
книги из ГПНТБ / Мамедов, А. А. Нарушения обсадных колонн при освоении и эксплуатации скважин и способы их предотвращения
.pdfны для выявления в теле трубы того или иного дефекта завод ского происхождения.
Исследованиями установлено, что в зависимости от положе ния оси колонны величины осевой силы будут различны [65, 73, 74]. В отличие от вертикального участка колонны, в искривлен ных участках ее от действия внутреннего и наружного избыточ ного давлений кроме осевых сил еще возникает изгибающий момент, в результате чего на этих участках величина и направ ление осевой силы меняются.
Поэтому при определении осевой силы нужно рассмотреть два случая.
1. В колонне отсутствуют искривленные участки — такое состояние колонны соответствует положению ее после натяжки.
Напишем относительное осевое удлинение колонны, исходя из обобщенного закона Гука:
е * = у К — Ц (сг< + 0Г,)Ь |
(117) |
Ввиду того, что колонна закреплена на устье и у границы цементного камня и возможность перемещения концов ее ис ключена, имеем условие ez = 0. Тогда получим
= |
Ц ( в / + |
0 , ) - |
Так как избыточное наружное гидростатическое давление жидкости возникает от разности наружного и внутреннего гид ростатического давлений, величина ot будет состоять из трех, а величина ог из двух слагаемых, т. е.
°t — ° t + |
|
°t+ |
у |
°r — Gr + |
°r У |
(119) |
|
где |
b2 + |
a2 |
|
|
Ь 2 —|- CL2 p. |
Aj-, |
|
o. = |
a, = |
|
|||||
-----!----- p |
---------0,05/yB; |
|
|||||
1 |
b2 — a2 |
|
1 |
b2— a2 |
|
|
|
|
2b2 |
0,05lyH; |
a' = —p; a'' |
= — 0,05/yB. |
(120) |
||
b2— a2 |
|
|
|
|
|
|
|
Решая совместно уравнения |
(118), (119) и (120), находим |
||||||
a, = |
- Д — [a2p - |
0,05/ (62у„ - a2yB)J. |
|
||||
|
b2 — a2 |
|
|
|
|
|
|
Или, заменяя а, b внутренним |
и наружным диаметрами, |
||||||
получаем |
|
|
|
|
|
|
|
= |
|
l d 2 p |
~ |
0 , 0 5 1 ( ° 2 у и ~ |
й*у Л |
(121) |
Умножая oz на величину площади поперечного сечения тру бы, находим осевую силу, возникающую в трубах колонны от
80
воздействия внутреннего рабочего и наружного избыточного гид ростатического давлений:
|
Р = i^ [d V - 0 |
,0 5 /(D 2v „ - d 27e)]. |
(122) |
2. |
В колонне имеются |
искривленные участки — это |
состоя |
ние соответствует моменту, когда колонна еще не натянута.
В этом случае на искривленном участке колонны от воздей ствия внутреннего рабочего и наружного избыточного гидро статического давления возникает изгибающий момент. Способ определения величины изгибающего момента описан в работах [66, 75]. В середине первой полуволны (с устья), где прогиб имеет максимальное значение, величина изгибающего момента определится формулой
М = f [d*/> - |
0,1 ‘{Н + - ^ ) (D*y. - d * Y . ) J |
U |
(123) |
где Н — прямолинейная |
часть колонны; L — длина |
полуволны. |
|
Как видно из формулы (123), с увеличением Н |
величина |
||
изгибающего момента М уменьшается. |
|
|
Это говорит о том, что наименьшая величина изгибающего момента будет у границы цементного камня и при этом состоя ние эксплуатационной колонны будет более устойчивым.
Можно определить, что осевая сила, которая приводит к возникновению изгибающего момента в искривленных участках колонны, имеет противоположное направление с осевой силой, возникающей в прямолинейном участке колонны.
Таким образом, представив момент как произведение осевой
силы на прогиб, для искривленного участка |
колонны |
получим |
|
Рг = f [*Р - 0 , 1 |
( # + у ) (D2Yh ~ d 2yB)] . |
(124) |
|
где Р1 — осевая сила в искривленной части колонны. |
колонны |
||
Величина осевой силы Рг в прямолинейной части |
|||
Н определится по формуле (122) |
|
|
|
Р2 = 3 L [<Рр- |
0,05Н (D*y„ - |
d V ]• |
(125) |
Осевая сила, действующая на всю длину колонны,
Р = Р1 + Р2,
или, пренебрегая величиной L/2 как малой и заменив Н на I—^кр, будем иметь
р = — f (! - ад ■d2P + f ( l - liV0Л (I - /кр) (D*yH- d*yB\ (126)
где /Кр — критическая длина искривленной части колонны.
81
Как видно из формул (122) и (126), осевые силы, возникаю щие в трубах эксплуатационной колонны в первом и во втором случаях, отличаются друг от друга как по величине, так и по направлению.
Так, если в первом случае внутреннее рабочее давление при водит к возникновению в трубах осевой растягивающей силы, то во втором случае — к осевой сжимающей силе. Обратное явление получается при воздействии на колонну наружного из быточного гидростатического давления жидкости.
РАБОТА ЭКСПЛУАТАЦИОННОЙ КОЛОННЫ В УСЛОВИЯХ ИЗМЕНЕНИЯ ТЕПЛОВОГО РЕЖИМА В СКВАЖИНЕ
При освоении и эксплуатации скважины в трубах эксплуа тационной колонны возникает осевая сила не только от воз действия внутреннего и наружного давлений, но и в результате изменения теплового режима в скважине.
В технической литературе величину осевой силы от воздей ствия теплового поля определяют, исходя из среднего значения температуры в трубах эксплуатационной колонны
Pt =aEFM, |
(127) |
где Pt — осевая сила, возникающая в колонне в |
результате |
нагрева или охлаждения; а — коэффициент линейного расши
рения; Е — модуль упругости |
материалов |
трубы; |
F — средняя |
площадь сечения трубы; At — средняя температура |
нагрева или |
||
охлаждения колонны |
|
|
|
_ {h — ^i) ~h (h — ^з) |
’ |
|
|
|
2 |
|
|
здесь tu t3— температура у |
устья и у |
границы |
цементного |
кольца до эксплуатации, принимаемая обычно по геотермиче скому градиенту; t2, /4— температура жидкости, движущейся по колонне в этих же местах.
Как видно из формулы (127), величина Pt определяется в предположении, что колонна по длине нагревается или охлаж дается равномерно. Такое допущение не соответствует реаль ному положению температурного состояния колонны в скважине
иприводит к получению завышенного значения осевой силы. Ввиду того, что с увеличением глубины скважины темпе
ратурный градиент возрастает и его воздействие на колонну может быть одной из основных причин нарушения герметичности эксплуатационной колонны, определение величины с учетом не равномерного изменения температуры в осевом направлении имеет большое значение.
Предположим, что обсадная колонна, состоящая из секций длиной /1, /2, ..., 1п, с площадью сечения трубы в каждой секции соответственно Fu F2, .... Fn спущена на глубину L и зацементи
82
рована до глубины /. Функции кривой распределения начальной температуры спущенной колонны в скважину обозначим ti(x)t а при работе скважины — t2(x).
На расстоянии х от устья скважины выделим в колонне бесконечно малый элемент dx. Полагая на участке dx разность изменения температуры постоянной, определим удлинение эле мента
dAlt = adx \t2 (х) —tx(*}].
Общее удлинение незацементированной части колонны со ставит
i
Alt = а [t2 (х) — tx(я)] dx.
Pis условия совместности деформации для колонны перемен ного сечения в работе [37] получена формула
i
a . [ [ t 2(x) — t1( x ] d x = ^ - + |
- ^ - + . . |
, + Д ^ - , |
||||
J l2 W |
^ |
|
EF, |
|
|
EFn |
откуда |
|
|
|
|
|
|
|
|
аЕ f [t2(х) — h (х)] dx |
|
|||
|
Pt = |
|
|
|
|
( 128> |
|
|
|
|
|
|
|
|
— |
+ — |
+ |
• • . + — |
|
|
|
|
Л |
Л |
|
Fn |
|
Как видно из формулы (128), для нахождения значения Pt |
||||||
функции t\(x) |
и t2(x) |
подлежат определению. |
При этом исхо |
|||
дим из следующих соображений [51]. |
tf ' t2 |
|||||
В период |
между |
окончанием |
цементи |
|||
ровки скважины и началом вскрытия про |
|
|||||
дуктивного пласта, т. |
е. |
до |
перфорации |
|
||
колонны, проходит достаточное время для |
|
|||||
уравновешивания температуры |
в |
системе |
|
|||
окружающая среда (горные породы — буро |
|
|||||
вой раствор)— обсадная |
колонна. Поэтому |
|
||||
изменение температуры в теле трубы обсад |
|
|||||
ной колонны в осевом направлении в этот |
|
|||||
период можно принять одинаковым с зако |
|
|||||
ном изменения температурного |
|
градиента |
|
|||
в горных породах в зависимости от глуби |
|
|||||
ны, как это показано на рис. 20, т. е. |
|
|||||
к(х) = ^ + |
1±=± -х , |
|
(129) |
|
где 11 — температура горных пород у устья скважины; to— температура на забое.
После вскрытия продуктивного пласта начинается период освоения и эксплуатации скважины. В этот период с выходом нефти и газа из пласта в колонне (при эксплуатации) или в связи с закачкой воды в пласт (при освоении) нарушается по стоянность теплового режима в радиальном направлении по всей ее длине.
Так как температура жидкости или газа в колонне U(x) отличается от температуры самой колонны ti(x), между ними происходит теплообмен.
В результате этого в случае 4(*)>4(*) колонна обсадных труб нагревается, а в случае 4 (* )< 4 (* )— охлаждается.
Для определения закона изменения температуры в теле тру бы, соответствующего этому периоду [функция t2(x)], рассмот рим кольцевой слой единичной длины, внутренний радиус кото
рого г, наружный r + dr. |
Допустим, что |
температура на |
его |
внутренней поверхности |
равна 4 (г), а |
на наружной — 4 |
(0 + |
+ dt2(r). |
|
|
|
Если считать, что количество тепла, протекающее в теле через какой-либо его слой, пропорционально поверхности этого слоя, тепловому напору, коэффициенту внутренней теплопровод ности и обратно пропорционально толщине слоя, тогда количе ство тепла, протекающее через рассмотренный нами кольцевой
слой в единицу времени, определится зависимостью |
|
q = K 2nrdt,(r) _ |
(130) |
dr |
|
Вследствие предполагаемой стационарности теплового со стояния q должно быть во всех концентрических слоях трубы одинаково.
Обозначив постоянную через m = q/2nK, можно дифференци альному уравнению (130) придать вид
dt.2(r) = т — .
Г
Интегрируя в интервалах от t2(r) до t2(x) и от г до а, получаем
/2 (г) — 4 (х) = т In — . |
(131) |
а |
|
Величину т следует определить из условия, что при г = Ь
4 {г) = 4 (х).
Тогда будем иметь
т _ 4 (*) — 4 (х)
Подставляя значение т в уравнение (131), находим
tl (*) In |
г |
Ь |
|
---- + |
t, (*) In---- |
|
|
h (Г, x) = -------- |
-— ----------- |
— • |
(132) |
in-5- a
Если закон изменения температуры жидкости или газа в ко лонне в период освоения и эксплуатации скважины считать прямолинейным, как это принято в работах [65, 73], то тогда для функции t3(x) будем иметь (см. рис. 20)
|
к ( х) = к~\— ° i |
~ х' |
(133) |
|
Подставляя значения функций ti(x) и t3(x) из формул |
(129) |
|||
и (133) в формулу (132), получаем |
|
|
||
, |
U— h \ , |
г , Л . *о |
|
|
к (г, х) = (- + —Г ~ х) |
Т + |
(/2+ |
(134) |
|
|
|
1п- |
|
|
Врезультате получили закон изменения температуры в тру бах обсадной колонны при нарушении постоянности теплового режима в скважине в радиальном и осевом направлении.
Вусловиях работы обсадной колонны изменение темпера туры в радиальном направлении (по толщине стенки трубы) очень мало влияет на ее напряженное состояние по сравнению
сизменением последней в осевом направлении. Кроме того, на
рушение незацементированной части обсадной колонны практи чески может произойти вследствие продольного изгиба нижней части этого участка, под действием осевой силы, вызванной из менением температуры в осевом направлении. Поэтому в даль нейшем величину температуры по толщине стенки трубы можно брать постоянной. Тогда для средней окружности сечения трубы, т. е. для r=(a + b)l2, зависимость изменения температуры по длине колонны примет вид
{ J. , ^0 — \ , а + Ь , Л , |
— h \ , |
2Ь |
|
|
||||
V 1 + ~ l ~ x ) ]n^ ~ |
+ {h + ~ l ~ xJ ] |
а + Ь |
(135) |
|||||
|
|
|
|
|
|
|
||
|
|
In — |
|
|
|
|
|
|
Подставляя значения |
ti(x) |
и |
^(я) |
из |
выражений (129) |
и |
||
(135) в формулу (128), |
находим |
величину |
осевой |
силы от |
из |
менения теплового режима в скважине при освоении и экс плуатации
85
a E |
, , + ^ |
, ) 1п^ |
+ ( „ + ^ |
, |
2b |
|
К + b |
||||||
|
|
|
|
|
||
p , = |
- > ti+l±i ~ * ) in“ ! djc |
|
|
|||
— + — + . . . + — V n — |
|
|||||
|
|
|||||
|
Л |
Л |
^ Fn J |
a |
|
или после решения интеграла и некоторых преобразований по лучаем
f Г L — h I i L — h i )
Pt = |
|
h |
|
|
1-2 |
|
In |
|
(136) |
|
|
“Ь |
• *H~” |
|
|
||||
|
|
Г |
|
г ~4” • |
p |
|
|
||
|
|
Fl |
|
F2 |
|
|
|
|
|
|
|
In ----- |
|
CL-f- Ь |
|
|
|||
|
|
|
2b |
|
|
||||
Здесь |
A |
— |
|
fl |
+ ft . |
В - |
|
|
|
= |
|
|
|
In - 2b |
|
|
|||
|
|
|
|
In — |
|
|
|
||
|
|
|
|
|
|
|
a + b |
|
|
a\ "f fl->+ |
• • |
|
• 4~ Gfl . |
b = |
bi -+- b2 -f- ■ • |
• 4~ bn |
(137) |
||
где au bx\ a2, b2\ |
|
an, bn — внутренние и наружные радиусы |
|||||||
трубы секции |
первой, |
|
второй и т. д.; |
п — число |
секций; |
h — |
высота подъема цемента за колонной.
Отметим, что для практического расчета можно упростить формулу (136). Для этого, во-первых, нужно рассматривать колонну односекционной, ввиду незначительной разницы в ве личинах площадей поперечных сечений отдельных секций (в расчете нужно брать площадь сечения трубы, имеющую сред нюю величину), во-вторых, для всех толщин стенок брать сред нее значение коэффициентов А и В.
Расчеты показывают, что коэффициенты А и В для труб диа метром 146 (толщина стенки от 6 до 12 мм) и 168 мм (толщина стенки от 6 до 14 мм) меняются соответственно от 0,445 до 0,518 и от —0,936 до —1,097, при этом средние значения А и В для всех размеров составляют 0,478 и —1,024. Учитывая это и
принимая для материала труб сх = 12-10~6 ---------- и £ = 2- 1х |
|
градус |
|
X Ю6 кгс/см2, получаем |
|
P t = Y ~ 14 — 1,0244) (L + h) — 0,0244 (L — К]. |
(138) |
Для сравнения результатов формул (127) и (138) определим величину Pt при следующих условиях: глубина скважины L= = 3500 м; высота подъема цемента h— 700 м; температура на
86
забое f0=90°C; температура жидкости на выкиде/2 = 60° С; тем пература горных пород у устья скважины ^ = 15° С. Произведем расчет эксплуатационной колонны диаметром 168 мм с толщи ной стенки 10 мм.
Следовательно,
p t = ± ^ о61 [(60 — 1,024 • 15) (3500 + 700) — 0,024 • 90 (3500 —
— 700,1 = 15422,4 кгс.
Теперь подсчитаем величину Pt по формуле (136). Величина At определится из выражения
(*«-*!)+ (*4-*»)
2
где
<. = |
'! + |
Л - |
У ^ |
= |
15 + |
(90 - |
1 5 ) - |g j - = |
75 °С; |
U = |
h + |
(to - |
to) ^ |
= |
60 + |
(93 - |
60) -g jjj- = |
84 °С; |
|
|
ы |
= ( 6 0 - 1 5 ) + ( 8 4 - 75) |
= 2 7 оС . |
|
|||
|
|
|
|
|
2 |
|
|
|
Pt = aEFAt = 12-10-6-2,1 • 106 *49,61 -27 = 33754,6.
Как видно из полученных результатов, при одних и тех же условиях величина Pt, подсчитанная по формуле (138), пример но в 2 раза меньше, чем по формуле (127).
Такая большая разница в значениях Pt, полученных по фор мулам (127) и (138), говорит о том, как велика роль неравно мерности распределения тепла в осевом направлении при опре делении осевой силы и что влиянием его пренебречь нельзя. Неравномерный нагрев или охлаждение эксплуатационной ко лонны по длине вызывает в сечениях трубы кроме осевой силы еще изгибающий момент.
Для определения величины изгибающего момента, возника ющего от неравномерного распределения тепла в осевом на правлении, допустим, что Т(х) представляет собой закон изме нения температуры, отсчитываемой от устья скважины. Примем, что все трубы эксплуатационной колонны разбиты плоскостями, перпендикулярными к оси скважины, на бесконечно тонкие кольца.
Вследствие температурного изменения, радиус любого коль ца колонны будет иметь увеличение (при нагревании) или уменьшение (при охлаждении). Причем, изменение радиуса кольца в разных сечениях колонны будет различным, и такая
87
неравномерная деформация стенки трубы по длине колонны бу дет вызывать изгибающий момент.
Такая же картина будет иметь место при взаимодействии боковой нагрузки от столба жидкости внутри и снаружи ко лонны в результате неравномерного распределения последней по стволу скважины.
Величина изгибающего момента может быть определена из дифференциального уравнения изгиба цилиндрического сосуда
7 7 + 4 р Ъ = -£ -, |
(139) |
где q — поперечная сила; D — цилиндрическая |
жесткость, рав- |
ная ——EZ3 — ; Е — модуль упругости материалов в кгс/см2;
р — коэффициент Пуассона материала трубы; г — средний ра диус трубы в см; б — толщина стенки трубы в см.
В нашей задаче q будет состоять из двух слагаемых, т. е.
<7= |
91 + 92. |
(140) |
где q1 — поперечная сила от |
столба |
жидкости внутри и снару |
жи колонны; q2— поперечная сила от изменения температуры. Величина qi для сечений, находящихся на расстоянии х от
начала координат (устья |
скважины), будет равна |
|
|
|
9i = |
0,l(Y„ — |
(141) |
где ун — удельный |
вес жидкости за колонной в гс/см3; |
||
ув — удельный |
вес жидкости внутри колонны |
в гс/см3. |
Для определения величины q2 предположим, что вследствие температурного изменения средний радиус колец колонны, на ходящихся на расстоянии * от начала координат, увеличивается на величину агТ(х). Это увеличение можно компенсировать, вернув кольцо к его первоначальному радиусу путем наложе ния внешнего давления интенсивности q2.
Радиальное перемещение от давления q2 равно r2q2/E6. При равнивая величины радиальных перемещений от температурного
расширения |
и внешнего давления, |
для q2 получаем |
|
|
|
|||
|
'9. = ~ - Т ( х ) , |
|
|
|
|
|
(142) |
|
где Т(х) определяется из зависимостей |
(129) |
и (135) |
|
|||||
|
t0- U |
а+ b |
|
. |
u - t 1 \ |
2Ъ |
||
|
(к + 1 |
" “ia |
+ ( h |
+ |
—— |
х |
I |
In------- |
Т (х) = t2 (х) |
— к (х)= ------- |
|
L |
|
д Ь |
|||
|
|
|
|
|
|
|
In —
а
88
или после соответствующих преобразований |
|
|
|
|||||
|
т(х) = А\в (Ч + |
|
^ |
+ |
Jl=Jl |
. |
(143) |
|
Из |
совместного решения |
уравнений |
(139), |
(140), |
(141), |
|||
(142) |
и (143) |
находим |
|
|
|
|
|
|
|
|
-^ - + 4Р4ц = |
Эх + |
ц, |
|
|
(144) |
|
где |
|
dx4 |
|
|
|
|
|
|
|
1 2 (1 -^ )0 4 4 ^ + ^) |
|
|
|
||||
|
|
|
|
|
||||
|
|
1 |
гб2 |
|
|
|
|
|
6 = "12£6з И |
{ ^ 7 Г [В 00 - |
V + 00 - VI + |
О'1(V. - |
Ц • |
(145) |
|||
Решение уравнения (144) |
представится в виде |
|
|
|||||
и == ер* (Сх cos $х + С2sin рх) + |
(С3cos рх + С4 sin рх) + |
11. |
Так как при увеличении х величина ер* возрастает, то по стоянные Сj и С2 должны быть равны нулю. Тогда имеем
и = |
(С3 cos Рх + £4 sin рх) + |
. |
(146) |
Постоянные С3 и С4 определяются из граничных условий, которые в нашем случае представятся в следующем виде:
при х = 0 |
и = 0; |
|
при х = I |
и = 0, |
(147) |
где / — длина незацементированной части колонны. Учитывая условия (147), находим
р _ |
Л |
^ |
_ |
Т| COS Р/ — (0/ И- Т]) fil |
(148) |
|
Сз |
4J34 |
* |
С 4 _ " |
4p4sinp/ |
||
|
Подставляя значения С3 и С4 в формулу (146), окончатель но определяем величину перемещения
4___9*+ ,> |
(149) |
4р4 sin р/ |
|
Изгибающий момент в сечении х определится по следующей формуле:
М = — Е1 d2U |
( 150) |
dx2 |
|
где / — момент инерции, равный б3/ 12.