Оценка качественной и количественной сходимости результатов расчета с экспериментальными данными
Экспериментальную проверку методики математического моделирования процесса бурения лучше всего выполнять путем сравнения результатов расчета механической скорости и крутящего момента с данными, полученными непосредственно в процессе глубокого бурения, поскольку только в этом случае функционирует колонна бурильных труб в качестве динамического элемента системы. Однако в промысловых условиях весьма сложно получить надежные зависимости показатФдейыттроцесса бурения от параметров режима и других факторов. Поэтому с целью проверки этих зависимостей, полученных расчетным путем, приходится прибегать к сравнению их также и со стендовыми данными, многие из которых широко известны и опубликованы в литературе.
Для удобства проверки надежности расчетных характеристик бурения приведем сначала некоторые результаты расчетов тех показателей бурового процесса, которые получаются непосредственно путем математического моделирования без привлечения каких-либо экспериментальных данных, за исключением тех, которые характеризуют свойства разбуриваемых пород. Таковыми являются, как было изложено выше, углубление зубцов долота в породу гтах, скорость проходки неизношенным долотом v и крутящий момент Мср.
Для определения взаимосвязи между показателями эффективности и параметрами режима необходимо пользоваться конкретными экспериментальными кривыми P(z). Поэтому искомые зависимости не могут быть получены в общем виде. Был использован метод математического эксперимента, в котором изменяли нагрузку на долото, скорость его вращения и проводили соответствующие расчеты для выбранных вариантов. горных пород.
Расчеты углубления зубцов долота, крутящего момента и скорости бурения были проведены для следующих реальных пород: индианского известняка, мрамора «коелга» и черного среднезернистого гранита. Эти три породы охватывают довольно широкий диапазон механических свойств. Индианский известняк при высоком давлении является весьма пластичной горной породой, мрамор — хрупко-пластичной и гранит — хрупкой.
Расчетные условия были выбраны, по возможности, таким образом, чтобы наиболее полно оценить влияние параметров режима и свойств породы на показатели эффективности бурения.
Так, параметры режима для всех трех пород варьировались в следующих пределах: осевая нагрузка на долото G = 5, 10, 15, 20, 25, 30 тс; скорость вращения долота п = 40, 70, 100, 150, 200, 250, 400, 700, 1000 об/мин.
Был использован вариант модели, в котором долото представлено тремя цилиндрическими катками с разным числом зубцов на каждом. Диаметр долота 246 мм, вооружение — призматические зубцы для расчета бурения на известняке и мраморе и полусферические штыри для расчетов на граните. Колонна бурильных труб — полубесконечная постоянного сечения.
В настоящей серии расчетов для щения долота использовали одну и зависимость P(z) для данной породы. В действительности динамическая характеристика сопротивления породы может изменяться при изменении скорости взаимодействия зубцов долота с забоем. Характер этих изменений и вызываемые ими последствия будут представлены ниже. Здесь же расчеты показали, что величина углубления зубцов в породу Zmax И КруТЯЩИЙ MO — мент Мср сравнительно слабо зависят от скорости вращения долота при изменении п от 40 до 1000 об/мин. Определяющими факторами являются свойства горной породы и осевая нагрузка на долото. На рис. 21 представлены зависимости zmax(G) для всех трех пород. Зависимости эти нелинейны, как и должно быть, учитывая сложный характер кривых P(z). Углубление зубцов в породу резко увеличивается с ростом осевой нагрузки.
различных скоростей вра — ту же экспериментальную |
zmaxiMM Рис. 21. Расчетные зависимости максимальной глубины внедрения зубца от осевой нагрузки Порода: 1 — индианский известняк; 2—мрамор; 3— гранит |
При небольших нагрузках (5 —
10 гс) темп роста zmax сравнительно невелик, но в интервале от 10 до 20 тс zmax увеличивается более чем в 3 раза.
Полученные расчетные зависимости крутящего момента от осевой нагрузки имеют характер, аналогичный зависимостям 2max(G). Этого следовало ожидать, поскольку крутящий момент в значительной степени определяется величиной углубления зубцов. Для сравнения качественного характера полученных расчетных зависимостей Mcp(G) с экспериментальной были использованы данные Д. С. Роули [52]. Так как эксперименты Роули производились долотом другого типоразмера, то сравнение выполнено в относительных величинах. Сравниваемые зависимости (рис. 22) демонстрируют хорошее качественное совпадение. Данные стендовых и промысловых исследований В. В. Симонова и Г. Д. Бревдо [28], Л. Е. Симонянца и
А. М. Некрасова |30] и других также показывают линейный или близкий к нему характер зависимостей Mcp(G). Это свидетельствует о том, что физическая интерпретация формирования крутящего момента при моделировании выполнена достаточно обоснованно и непротиворечиво.
Полученные расчетные зависимости начальной механической скорости бурения от скорости вращения долота очень близки к линейным, но в некоторых случаях отклоняются от нее в сторону снижения темпа роста v (рис. 23, а). Это отклонение особенно отчетливо проявляется при высоких осевых нагрузках на долото. Такого рода зависимости v(n) достаточно хорошо известны в литературе и подтверждены многочисленными стендовыми и промысловыми экспериментальными исследованиями [13, 33, 50, 51].
ц% Рис. 24. Сравнение расчетной зависимости скорости проходки от осевой нагрузки с экспериментальной зависимостью по Бингхэму: 1 — расчетная зависимость; 2 — экспериментальная за висимость |
Часто экспериментально наблюдается более значительное снижение темпа роста скорости при увеличении числа оборотов. Причины этого явления будут подробно исследованы ниже.
а 6 Рис. 23. Расчетные зависимости скорости проходки от числа оборотов (а) и осевой нагрузки (б) |
Кривые зависимости механической скорости от осевой нагрузки на долото также стабильны для различных пород. Они характеризуются нарастающим темпом увеличения v при малых нагрузках и затем почти линейной зависимостью (рис. 23,6). Подобные зависимости были получены в стендовых и промысловых условиях Д. Спиром, П. Муром, М. Бингхэмом [13, 50, 51, 53] и рядом других исследователей. На рис. 24 дано сравнение в относительных величинах расчетных кривых с экспериментальными, полученными М. Бингэмом [13]. Расчетные данные качественно достаточно хорошо согла
суются с общеизвестными стендовыми и экспериментальными данными.
р>тс л /О ;1УУЛ/ |
Рис. 25. Экспериментальные зависимости силы сопротивления породы от глубины внедрения зубца: а — мелкозернистый слюдистый песчаник, глубина 2460—2600 м, зубец 10X1 мм; б — известняк верхнего карбона, глубина 770—970 м, зубец 5X1 мм; в — мелкокристаллический плотный известняк, глубина 1500—1800 м, зубец 5X1 мм; г — мелкокристаллический плотный известняк, глубина 1900—2200 м, зубец 5X1 мм |
1 Z г, мн |
Качественное совпадение расчетных и экспериментальных зависимостей дает основание полагать, что предложенная модель процесса бурения достаточно хорошо воспроизводит физические явления, происходящие при разрушении забоя скважины долотом.
Полезно, однако, попытаться сделать количественную проверку совпадения расчетной величины скорости бурения с реальными промысловыми данными, т. е. провести так . называемую пассивную проверку. Такая проверка была произведена по нескольким регионам, где имелась возможность получения кернового материала и более или менее достоверных данных по отработке буровых долот.
Керны были отобраны в нескольких различных регионах и с различных глубин: мелкозернистый слюдистый песчаник
(2460—2600 м), известняк верхнего карбона (770—970 м), мелкокристаллический плотный известняк (1500—1800 и 1900— 2200 м). Лабораторная обработка кернов производилась в условиях, соответствующих глубине залегания породы с учетом параметров глинистого раствора. Соответствующие зависимости P{z) и V(г) приведены на рис. 25 и 26. Для моделирования были отобраны конструкции долот, наиболее употребительные в каждом конкретном регионе. Конструктивные параметры этих долот были введены в модель.
Диапазон применявшихся на практике нагрузок колеблется в сравнительно незначительных пределах. При расчетах этот диапазон был существенно расширен. Таким образом, сравнение показателей могло быть выполнено только на отдельных участках расчетных зависимостей. Расчетные зависимости скорости проходки и от удельной осевой нагрузки на долото даны на рис. 27. В трех сравниваемых вариантах (рис. 27 а, в, г) применяли турбинный способ бурения и фактическое число оборотов определялось по моментным характеристикам турбобуров, в четвертом случае (рис. 27,6) —электробурение.
Если принять во внимание, что горизонты, в которых про-
Рис. 26. Экспериментальные зависимости объема разрушенной породы от глубины внедрения зубца: |
а — мелкозернистый слюдистый песчаник, зубец 10X1; б—-известняк верхнего карбона (!), мелкокристаллические плотные известняки (2, 3), зубец 5×1
Рис. 27. Сравнение расчетных данных по скорости проходки с промысловыми экспериментальными данными:
сплошная линия — расчетные зависимости, кружки — экспериментальные точки
изводилась отработка долот, не однородны по литологическому составу и часто содержат кроме исследованной еще и ряд других пород, то следовало ожидать, что только часть точек должна совпадать с расчетными кривыми. В этом смысле результаты сопоставления можно считать в высшей степени удовлетворительными. Несомненно, что порядок величин, получаемых расчетным путем, является вполне достоверным.
Полученные обнадеживающие результаты послужили основанием для активной и более детальной проверки детерминированной модели процесса бурения. На экспериментальной буровой— стенде была спланирована и осуществлена специальная программа для активной проверки функционирования детерминированной модели. Методика проведения экспериментов предусматривала бурение двух параллельных стволов. В одном из них проводился сплошной отбор керна, в другом — выполнялись буровые эксперименты при различных значениях осевой нагрузки и скорости вращения долота.
Бурили электробуром с проходкой 0,4 м при каждой комбинации нагрузки на долото и скорости его вращения. Одновременно непрерывно отбирали шлам. Исследования шлама и кернов позволили составить крупномасштабный классификационный разрез пробуренного интервала и выбрать керны для лабораторного исследования!.
Условия проведения экспериментов.
— Разбуриваемые породы:
известняк, категория твердости по штампу 4—5 и 5—6;
глина, категория твердости по штампу 2;
алевролито-песчано-глинистые разности, категория твердости по штампу 3.
— Забойное давление:
250 кгс/см2, соответствующее глубине 2000 м.
—Горное давление:
500 кгс/см2.
— Параметры глинистого раствора:
плотность 1,14—1,15 г/см3;
водоотдача 7 см3 за 30 мин.
Учитывая, что перечисленные выше породы слабопроницаемы и процесс динамической фильтрации при лабораторном исследовании кернов практически не развивается, поровое давление принято равным нулю.
Скорости внедрения зубца в породу при лабораторных экспериментах соответствовали скорости вращения долота п = = 225 и 670 об/мин. ни а
Количество кернов, использованных при проведении эке-
10 |
чериментов, определялось из условия необходимости 6—8 повторений.
_£ 1^ |
Выходными данными при лабораторном исследовании кернов являлись зависимости P(z) и V(.z).
°8 |
о О |
Входные параметры модели:
пользования различных методов проверки, сводится к тому, что качественная и количественная сходимость результатов расчета на базе математического моделирования процесса бурения является удовлетворительной.
Таблица 2
Интервал бурения, м |
п, об/мин |
G, тс |
V, эксперимен тальная |
м/ч расчетная |
Отклонение, и / 0 |
2000—2005 |
670 |
17,5 |
18,6; 13,7 |
16,7 |
3 |
11,8 |
7,3; 6,5 |
8,3 |
20 |
||
2074—2078,5 |
225 |
17,8 |
11,8; 10,5 |
12,1 |
8 |
10,2 |
5,1; 6,0 |
6,3 |
13 |