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

книги из ГПНТБ / Мамедов, А. А. Нарушения обсадных колонн при освоении и эксплуатации скважин и способы их предотвращения

.pdf
Скачиваний:
13
Добавлен:
22.10.2023
Размер:
17.87 Mб
Скачать

ны для выявления в теле трубы того или иного дефекта завод­ ского происхождения.

Исследованиями установлено, что в зависимости от положе­ ния оси колонны величины осевой силы будут различны [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

b2a2

 

 

 

2b2

0,05lyH;

a' = —p; a''

= — 0,05/yB.

(120)

b2a2

 

 

 

 

 

 

Решая совместно уравнения

(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

Рис. 20. Схема изме­ нения температуры в скважине в осевом направлении
83

рована до глубины /. Функции кривой распределения начальной температуры спущенной колонны в скважину обозначим 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), получаем

 

 

,

Uh \ ,

г , Л . *о

 

к (г, х) = (- + Г ~ х)

Т +

(/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.

Соседние файлы в папке книги из ГПНТБ