Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Основы механики глубокого бурения

..pdf
Скачиваний:
5
Добавлен:
12.11.2023
Размер:
7.62 Mб
Скачать

Рис. 7.4. К выявлению зависимости скорости

л

бурения от проходки

и

О

АН

Ун

Начальное условие: о(0) = VQ. Решение этого уравнения (7.7) найдено выше.

Рассмотрим уравнение (7.6) несколько с иных позиций. Пред­ ставим его левую часть в виде

dv _ dH d v _ dv dt dt dH dH'

что после подстановки в (7.6) дает:

(7.13)

Здесь Я г- текущая глубина скважины, а через Но обозначим на­ чальную глубину бурения. Тогда начальные условия для уравне­ ния (7.13) при t - 0 запишутся: Я = Яо, v = Vo-

Поместим новое начало координатной оси в точку Н - Но (см. рис. 7.4). Тогда начальные условия примут следующий вид: Я = 0, v = VQ. Уравнение (7.13) - уравнение с разделяющимися переменными и его решение после несложных преобразований можно представить так:

(7.14)

71

Значение же механической скорости бурения v(AН) в конце интервала АНсогласно (7.14) запишется:

1

»(ДЯ) = Ц 1 - ! Ь ^ ! Й Д Я | !~\

(7.15)

Рассмотрим теперь несколько частных случаев зависимости

(7.15).

Очевидно, что при р. = О имеем v(AH) = vo, то есть, как и сле­ довало ожидать в данном случае, износ долота отсутствует и зна­ чение v = vo сохраняется на всем интервале АН.

Пусть теперь р. Ф0. Рассмотрим решение (7.15) для ряда зна­ чений параметра ц:

Л= —1, v(AH) =v0

Л= 0, o(AH) = n0| l - | - A t f j 2;

11 = 1, »(Aff) = q,-(l-b -A ffJ;

11 = 2, о(АЯ) = о0 -е-'‘дН;

ц = 3, v(АН) = v01

1 + |к»0ДЯ ’

А теперь проанализируем полученные зависимости.

При значениях Т| = -1, 0, 1 (три верхние зависимости), видно, что в случае больших значений скоростей vQ(и малых интерва­ лах АН), как и следовало ожидать (после пренебрежения вторы­ ми членами в скобках при о0 -» °°), скорость на нижней границе интервала v(AH) = v0. В случае значения ц = 2 при больших зна­ чениях vo скорость t>(AН) на нижней границе интервала будет также возрастать, однако всегда будет меньше скорости г>0 в ецЛ// (сказывается износ долота при прохождении интервала АН). На­ конец, для случая т] = 3 имеем:

limn (ЛЯ) = lim -— Зй-

1

ци0ДЯ

идя ‘

72

Мы видим, что в двух последних случаях, несмотря на весьма большую начальную скорость бурения, на нижней границе ин­ тервала ДН скорость бурения отличается от начальной, причем, по логике вещей, и в последнем случае = 3) она должна быть меньше значения механической скорости бурения Vo в силу изно­ са долота на данном интервале: v(AH) < г>0.

Рассмотрим теперь изменение глубины интервала Н в зави­ симости от времени. С этой целью проинтегрируем выражение (7.7):

2-д'

 

 

11

Я (0 =

И (2 -Ч)

(7.16)

 

Очевидно, что при р = 0 износ долота отсутствует, скорость механического бурения v = оо, потому углубление интервала за­ пишется как

H(t) = v0-t.

Следовательно, в данном случае теоретически одним долотом можно пробурить интервал сколь угодно большой величины.

Рассмотрим решение (7.16) при р * 0 для ряда значений па­ раметра ц:

Т1 = Ч Я ( 0 =

л = 1, ^ ( 0 = ^- *(! ~

);

т| = 2, Я (0 = - ■In (1+рг>0 t);

i i = 3 . H ( t ) = ^ - ( , / T T W f - 4

Из данных выражений следует, что при ц = -1 максимально возможное время бурения fmax найдется при обращении в нуль

73

разности в круглых скобках (при ее отрицательном значении под знаком квадратного корня глубина интервала будет комплексным числом, что в данном случае не имеет физического смысла), и, естественно, при этом значении времени найдется и максимально возможная глубина интервала Ятах:

= 1 1

зц'

Рассмотрим теперь случай ц - 0. Здесь так же очевидно, что Ятах будет достигаться при обращении в нуль разности в круг­ лых скобках (при дальнейшем возрастании времени бурения глубина интервала вместо увеличения будет уменьшаться, что лишено физического смысла). Отсюда

tn

=

я

= -3L

 

ц

,пах

2ц

Случай Т| = 1 дает нам теоретическое время бурения, равным бесконечности, но при этом величина интервала имеет конечное значение:

И, наконец, при значениях параметра ц = 2 и "П= 3 согласно математической модели углубления забоя (7 .6 ) время бурения и глубина пройденного интервала могут быть сколь угодно большими. Но вернемся к скоростям механического бурения. В случаях значений параметра т| = - 1 , 0 , 1 , как было уста­ новлено выше, падение механической скорости из-за износа долота меньше, чем в случаях значений т| = 2, 3. Однако, проходки на долото в первых трех случаях ограничены, в то время как в последних двух бурение, теоретически, может вес­ тись одним долотом неограниченное время на неограниченную глубину. Отсюда следует парадоксальный вывод: при менее из­ нашиваемом долоте проходка меньше, чем при более изнаши­ ваемом.

Итак, проведенный анализ показал, что в области математиче­ ского моделирования процессов углубления забоя царит, вообще говоря, изрядная неразбериха. Следствием этого является тот факт, что математических моделей углубления забоя к настоя­ щему времени опубликовано уже несколько десятков. И что ещё необходимо заметить: в подавляющем большинстве предлагае­ мые зависимости игнорируют механические свойства буриль­

74

ной колонны (это можно увидеть из продемонстрированных в настоящей лекции примеров). К чему приводит игнорирование механических свойств БК мы и попытаемся разобраться в сле­ дующей лекции.

ЛИТЕРАТУРА

1.Бреодо ГД. Проектирование режима бурения. - М.: Недра, 1988.

2.By Тхьсн Лыонг. Исследование динамической устойчивости процесса углуб­ ления скважины в породах кристаллического фундамента месторождений Боль­ шой Тигр и Дракон Вьетнама (синергетический подход): Авторсф. дне. канд-та техн. наук. -- М„ 1999.

3.Григулецкиы В.Г. Оптимальное управление при бурении скважин. - М.: Не­ дра, 1988.'

4.Меджидов Г.Н. Повышение эффективности процесса бурения скважин с

учётом закономерностей разрушения горных пород в осложнённых условиях:

Автореф. дне. канд-та техн. наук. - Баку, 1999.

о. Середа Н.Г., Соловьёв Е.М. Бурение нефтяных и газовых скважин. - М: Недра, 1988.

Лекция 8

ТРАНСФОРМАЦИЯ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ УГЛУБЛЕНИЯ ЗАБОЯ ПРИ НАЛИЧИИ КРУТИЛЬНЫХ АВТОКОЛЕБАНИЙ БУРИЛЬНОЙ КОЛОННЫ

При создании методов прогнозирования эффективности рабо­ ты породоразрушающего инструмента на забое скважины и оп­ тимизации режимов бурения, как правило, используются эмпи­ рические зависимости, связывающие показатели отработки долот с режимными параметрами, без учета механических свойств БК. Подобный подход, как мы увидели в предыдущей лекции, поро­ дил большое количество указанных зависимостей, довольно час­ то не согласующихся друг с другом. Отмеченные факты, при прочих равных условиях, имеют место из-за пренебрежения вза­ имным влиянием долота и колонны друг на друга, что приводит к ложной трактовке результатов исследований закономерностей процесса бурения нефтяных и газовых скважин со всеми выте­ кающими отсюда последствиями для теории и практики этой области техники. Рассмотрению указанных проблем и посвящен излагаемый ниже материал.

Многие из зависимостей, подобных приведенным выше, по­ лучаются в лабораторных условиях по схеме, изображенной на рис. 8.1 (учет осевой нагрузки на долото Р, скорости его враще­ ния п, и иногда расхода промывочной жидкости Qp при полном игнорировании остальных факторов, которые «загоняются» в эмпирические коэффициенты). Тем не менее сильное влияние механических характеристик БК на процесс разрушения горной породы отмечается во многих исследованиях. Однако, ни в одной из приведенных в предыдущей лекции зависимостей, используе­ мых при проектировании режимов бурения скважин, не нашли отражения геометрические и механические характеристики ком­ поновок бурильных колонн, которыми в каждом конкретном случае ведется бурение скважины. Короче говоря, представлен­ ные выше математические зависимости процесса углубления за­ боя не отражают влияния на него такого важнейшего элемента, как бурильная колонна, осуществляющая связь между забоем скважины и дневной поверхностью.

76

Случай 2. Точка А находится в зоне РВ, а точки В* и С - в зоне НВ; сочетания параметров (Рь «,), (Р3, щ), (Р3, п2).

Случай 3. Все три опытные точки А*, В', С лежат в зоне НВ, что дает сочетания (Р2, щ), (Р3, щ), (Р3, п2).

Случай 4. Две точки - точка А и точка С - находятся в зоне РВ, а точка А* - в зоне НВ; соответствующие сочетания здесь (Pi, «О, (Ръ Щ), (Р2, п2).

Представленные случаи описывают основные сочетания па­ раметров (Ph tij) при проведении эксперимента в процессе буре­ ния вертикальной скважины на заданной глубине Н. Пусть про­ водится отработка долот. Тогда для определения в формуле (8.1) неизвестных постоянных А, ос, Р необходимы три уравнения. По­ лучить их можно опытным путем, замерив три значения v при

различных комбинациях параметров Р и щ:

 

Ц = АР?п\, v2 = АР?п\, v3= AP?nl

(8.3)

Прологарифмировав эти три соотношения, получаем:

 

In Л + alnP, -I-РIn я, = 1пг>,,

 

In А + aln Р2 +р Inщ = Inv2,

(8.4)

1пЛ + а1пР3 +р1пяз = 1па3.

 

Система уравнений (8.4) является линейной системой с неиз­ вестными величинами In А, а и р. Согласно схеме значений па­ раметров Р и п, представленной на рис. 8.2, мы будем использо­ вать следующие сочетания (Р„ nj)\

(Рь rii), (Рь wj), (Pj, и,) V г, 7 = 1 , 2 при условии i * j.

В этом случае решение системы (8.4) принимает наиболее простую форму:

In Ч

in^L

 

a =—

р-----Vi = 1,2. j = 1,2, i * j - (8.5)

lni l

ln-i-

pi nj

Pj

Tij

 

Записанные соотношения нужно понимать следующим обра­

зом.

Если в первом опыте (из серии трех экспериментов) берется пара значений (Pi, ni), то во втором - (Pi, п2) (то есть осевая

79

нагрузка остается прежней, а изменяется скорость вращения до­ лота со значения щ на значение п2), а в третьем опыте - (Р2, п2) (здесь остается без изменения скорость вращения породоразру­ шающего инструмента, но изменяется осевая нагрузка с Pi на Р2). Подобный способ выбора изменения параметров Р и и0 при про­ ведении экспериментов гарантирует выполнение условия нера­ венства нулю главного определителя системы (8.4) и существен­ но упрощает процедуру вычисления опытных констант А, а и р. А теперь будем рассуждать следующим образом.

Пусть результаты стендовых экспериментов для некоторой пары «долото - горная порода» представлены формулой (8.1). В случае 1 разница между механической скоростью в стендовых (равномерное вращение) и механической скоростью в промысло­ вых условиях для сходных типов долота и разбуриваемой поро­ ды будет минимальной, а потому примем в виде стендовой зави­ симости (8 .1 ) механическую скорость бурения при равномерном вращении бурильной колонны. Соотношения (8.3) в случае 1 запишутся как

у,, = APfnf, у12 = APfrtf, у22 = AP?rtf,

и экспериментальные параметры а, р и А согласно соотношений (8.5) суть:

\пМ

lnfk

 

 

 

 

a = _% L, p = - 3 i , A = -5 V .

J f

(8 .6 )

b i

p

ina.

»5

 

P2

n2

 

 

 

 

Поставим теперь следующий вопрос: как изменятся парамет­ ры а, Р, А в случаях 2, 3 и 4 сочетаний(Р„ щ)? Вспомним, что при автоколебаниях БК механическая скорость бурения у уменьшается, а ее значение в зоне НВ оценивается согласно

(8.2).

Рассмотрим случай 2. Если бы здесь имело место равномер­ ное вращение БК при всех сочетаниях (Р,, п}), то соотношения (8.3) записались бы как

уи = APfrtf, у31 = АР*rtf, у32 = АР*п\,

причем числовые значения ос, Р, А, определенные из данных уравнений, должны были бы совпадать с параметрами (8 .6 ). Од­ нако в этом случае при сочетаниях (Р3, щ) и (Р3, га2) возникают