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

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

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

Как видно из выражения (89), максимальное значение т получается при а = 45°, тогда имеем

ог = а, =

рт

 

 

 

2 ( т ■ ■1)

 

 

(90)'

 

 

 

 

 

,Р(т — 2)

 

 

 

ал

2 (т — 1)

 

 

 

Определив значение вертикального давления р из формулы

(87) и подставив его в выражение (90), для

тт ах найдем

_

— 2) куа

-ctgcp.

(91)

4 ( т — 1) (0,022 + /с2)

 

 

 

Из совместного решения уравнений (88) и (91) определяем

радиус податливой плоскости

 

 

 

 

4 — 1) (0,022 -f- к2) тс

 

(92)

— 2) ку ctg ср

 

 

 

Определим значение а

при следующих

данных:

ф = 30°;

у= 2 тс/м3; к = 0,42; тс=1 кгс/см2.

Сначала подсчитаем величину т

 

 

1 — sin ф

1 — sin 30°

 

 

При

принятых

значениях

<р,- у,

/с, т,

тс из формулы

(92)

получим

а = 8,2 м.

Следовательно,

до значения радиуса подат­

ливой плоскости, равного 8,2

м, не

будет

образовываться

свод

и не будет нарушений обсадных колонн, связанных с воздей­ ствием осевой нагрузки.

Величина радиуса податливой плоскости (как видно из рис. 15) связана с радиусом выработки и длиной (мощность) фильтра следующим образом

a = e + h tg (4 5

°---- 2. У

(93)

Как видно из уравнения (93),

максимальное

значение а

получится за счет или длины фильтра, или радиуса выработки. Отсюда заключаем, что при коротких фильтрах (h < а) обра­ зование свода будет за счет радиуса выработки, а при длинных фильтрах ( h ^ а) — длины фильтра, т. е. е= 0.

Подставляя значение а из выражения (92) в формулы (84) и (87), получаем величины осевой нагрузки и бокового давле­ ния от «свободного тела» при образовании свода

41,6 0,1) (т — I)2 (0,022 + к2) л п 2с

Q =

— 2)2 ку ctg ф

Ph(mах)

( 94)

2 — 1) ктс

т — 2

60

Подсчитаем величины Q в тс и рцmax) в кгс/см2 для обсад­ ной колонны диаметром 168 мм при принятых выше значениях Ф, у, к, тс

Q =

41,6(0,42 — 0,1) (4 — I)2 (0,022 + 0,422) 3,14 • 0,084 • 102

104;

(4 — 2)2 • 0,42 - 2 - 1 , 7 3

 

 

 

 

 

 

 

Ph(max)

2 (4 — 1) 0,42 •

1

1,26.

 

 

4 — 2

 

 

 

 

 

 

 

Как видно из полученных результатов, величина бокового давления от «свободного тела» незначительна. Поэтому при дальнейшем исследовании его действием можно пре­ небречь. Что касается ве­ личины осевой нагрузки, то она способна привести к продольному изгибу фильтровой части колонн или разрыву муфтовых соединений выше фильт­ ра.

Осевая сила была оп­

 

ределена

без учета

угла

 

падения пласта,

т.

е. в

 

случае

горизонтального

 

расположения пласта.

 

В случае горизонталь­

 

ного

расположения

пла­

 

ста

масса

горных

пород,

 

заключенная вутри свода,

Рис. 16. Схема образования свода над

представляющего

полови­

фильтром при большом угле падения

ну эллипсоида, имеет сим­

пласта

метричное

расположение

 

по отношению к оси

колонны, и благодаря этому в зоне фильт­

ра возникает только

осевая

нагрузка (см. рис. 15). Когда угол

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

Определим величины изгибающего момента и осевого уси­ лия с учетом угла падения пласта.

На рис. 16 представлена схема свода, образованного над фильтром, без окружающей массы горных пород. Как видно из схемы, масса горных пород, заключенная внутри свода, в за­ висимости от угла падения пласта распределяется неравномер­ но по отношению к оси колонны.

61

Величина изгибающего момента при неравномерном рас­ пределении массы определится по формуле

 

M = (F1 — Ft)r,

 

(95)

где Fь

F2— силы трения,

возникающие между трубами

эксплуатационной колонны

и

горной

породой,

заключенной

внутри

свода, соответственно

слева и

справа;

г — наружный

радиус трубы.

 

 

 

 

Для определения величины F\ и F2 можно пользоваться вы­

ражением для определения силы трения, т. е.

 

 

dF =

2nrKpvdz tg ф,

(96)

где pv определяется по формуле (81).

Формулы (81) и (96) получены при условии отсутствия угла падения пласта, т. е. в случае, когда сс = 0.

Для возможности применения формул (81) и (96) к рас­ сматриваемому случаю, когда афб, из каждой половины сво­ да мысленно отбросим одинаковые объемы, которые соответ­ ствуют почти равным треугольникам AEF и BCD (см. рис. 16). Тогда рассматриваемые своды с наклонными основаниями превратятся в своды с горизонтальными основаниями, и для определения сил трения F\ и F2 в этом случае можно пользо­ ваться следующими выражениями:

 

 

 

ь

 

F ,= ±

ь

 

 

 

 

F ^ ^ - j d F ;

f dF,

 

(97)

 

 

 

h1

 

 

h,-\-h2

 

 

где

 

 

 

 

 

 

 

 

 

 

= t a r) tg a;

h2 = (a + r) tg a.

 

 

Из

совместного решения

уравнений

(81), (96) и (97) най­

дем

 

 

__

mi.

 

 

 

 

 

 

 

 

 

 

 

=

1 8

-------

-e

4 ---- ^---- 2г|>2----- Д [arcsinjAn —

 

4 + ri2

 

 

 

 

 

— (1 — 2^) |/>!(1 —

-

4 -f

- [2\p! — ц

— T|?i)

— 1] X

 

 

 

 

if

 

 

 

 

 

у

e -T) arc sin

_[_ 2 t|)1| ;

 

 

 

F« = A{

2T1

 

4

---- 2 (фх +

ф2)2 —

 

 

4 + t}2

 

 

 

 

 

 

 

 

----- -

arc sin \ f

+.ф2— [1—2^ i+ T 2)lK(Ti+T2) [1—Л +

Ф2)]]—

 

■[2 vTi + ^

—Л V СФ1 + Ф2' [ i —(Ti +

 

1] X

62

 

 

у е - л arc sin

V Ч),+ф! - f

2 (фх +

ф2)

,

где

 

 

 

 

 

 

 

 

 

\ _

4я"су6г *8 Ф

^

'

(о— r-)tgg .

 

(a + r) tg a

 

4 +

T)2 ’

Vl

 

 

2b

 

2b

Подставляя

значения Fi

и

F2 в

формулу

(95), получаем

м

= А г {-%- [arcsin \/ ф7ф-фГ — 11—2(фх + ф2) X

X КОИ + ФгШ — (Фх + W] ] ---- [arc sin ГфГ— (1 — 2фх) X

X / ф ] (1 — Фх)------ л- - 4

[2фх — т) )/ % (1 — Ф1) — 1 ] X

X е-ч « s>" ^

,

[2 (фх +

ф2) -т^Х фх+Ф г'П -Х Ф х+Ф г']-

 

 

4 + г)2

 

 

 

 

 

 

 

l] t -r\a rc sin УфГРфГ —

2 % -f 2 [(Яр! +

Ч5^ 2 — ^?] } .

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

|[arcsin р ф , +

ф2 — [1 — 2 (фх +

q>2) X

X V (Ф х + ф у [1 —

(Ф х + Ф[arcа '1 —sin V b —

(1 — 2ф,) v 4 i(l —Ф11 ! | .

(98)

Величина осевой силы в этом случае определится как сум­ ма сил трения

Q

Чтобы можно было пользоваться выражениями (81) и (96), из первого члена правой стороны формулы (99), который пред­ ставляет собой суммарную силу трения массы горных пород левой стороны свода, отброшена часть объема, соответствую­ щая треугольнику AFF, а ко второму члену, представляющему собой суммарную силу трения массы горных пород правой стороны свода, прибавлена часть такого же объема. В резуль­ тате получим

 

А ( лт]

_ iHL

Q

i L e 4 1 - 2 ( ф 2 + ф , 2 ) _

2~ (~Т~

 

4 + rf

63

------[arc sin Yг|>1 +

arc sin l/if>2

(1 —

 

'M l — Ф1 ) —

(1 1|)2) l i|>2 (1 —

U V l

-----7 -7— -

[2 -ф! —

г) V 4>i (1 —

4'i) — 1] X

X e“ 11 arcsin V ¥ --------

[2Ц», - 1 1

l A

h l l - i W

-

 

 

4 + T]2

 

 

 

 

 

__ J J 0 — ri arc sin ^

 

+ 2;ii31

+

ii)2'J .

 

В этом случае влиянием второго, четвертого, пятого, шесто­ го и седьмого членов можно пренебречь; тогда

0 = — {( — --Г ]------—[arc sin]Api -f arc sin V ф2

2 4 J 2

— (1 — 'M K 'M l — Ф1) — (1 — 'ЫУЧгО —Ф2Ч I • (1°°)

Отметим, что для практических целей можно упростить фор­ мулы (98) и (100). Для этого нужно отбросить величину г, входящую в эти формулы, так как она по сравнению с величи­ ной а очень незначительна. Тогда для расчета будем иметь

M =

4 -j- т]2

)t

 

( 101)

 

 

 

Anrkyb2

Г

— 1,28v)

ЛГ]

4 + +

(1

- 1 ] tg (P,

!

 

4

где

v = arcsin

vx = arc sin и

a tg а

( 1

2 Ъ

\

a tg a

t^

a t g a \ ■ /

a tg a / j

a tg a \ ,

b J

V

2b V

2b J

 

 

( 102)

2 a tg a \

 

/ a tg a / j

a tg a \

Подставляя значения b и rj из выражений (75), (80) в урав­ нение (101) и производя некоторые преобразования, получаем

лл

1,65л (£га)2 у

N ,

 

М

= ---------- — (Vi — v) ctg cp;

 

 

0 ,0 2 2 + /г2

v

 

 

[kya^^ __ 1 ,28v)/e — 0,1] ctg tp.

(103)

0,022 + k"-

 

 

Формулы (ЮЗ) позволяют приблизительно определить величины изгибающего момента и осевой силы, действующих

64

на обсадную колонну в зоне фильтра при наличии угла падения

пласта.

В случае, когда а = 0, будем иметь М = 0, а выражение для Q в этом случае превратится в формулу (94) для определения величины осевой нагрузки в случае отсутствия угла падения пласта.

П р и м ер .

Определить М и Q для

обсадной колонны диаметром 168 мм

при

а = 30°|; тс = 1,1

кгс/см2; ср = 30°; /с =

1,42; у= 2 тс/м3.

 

Для

облегчения

расчета

на

 

рис. 17, 18 представлены графики за­

 

висимостей между радиусом основа­

 

ния свода а и углом внутреннего

 

трения

грунта ср (при тс= 1 кгс/см2),

 

подсчитанная по формуле (92), а так­

 

же

между

коэффициентом v и углом

 

падения пласта а для различных зна­

 

чений

ф по формуле (102).

 

 

 

Из рис. 17 находим а. При зна­

 

чении

ф = 30°,

его

величина

равна

 

8,2 м. Ввиду того, что нам нужно

 

знать

величину

а

для

случая,

когда

 

тс =

1,1, найденное значение умножа­

 

ем

на

1,1.

Тогда получим а = 9,02 м.

 

На

рис. 18

значение коэффициента v

 

Рис. 17. Зависимость радиуса основа­

Рис. 18. Зависимость коэффициен­

ния свода от угла внутреннего тре­

та v от угла падения пласта

ния грунта

 

при а = 30° с углом внутреннего трения грунта 30° равно 26-10—3 (третья кри­

вая). Величину vi

при указанных значениях а и ф находим по формуле (102),

равной 67,3-10_3.

Таким образом

 

 

1,65-3,14 (0,42-0,084-9,02)2-2

- Ю-з) 1,733 = 0,422;

 

М =

 

- (67 -10 - з _ 26

 

0 ,0 2 2 +

0,422

 

 

2,6-3,14-0,084-0,42-2-9,022

 

^ =

0,022 + 0,422

1(1 ' 1,28-26-10-3)0,42 — 0,1] 1,732=122,6.

Осевая сила без учета угла падения пласта, т. е. по формуле (94), при вышеуказанных данных равна 130,2 тс.

3 Зак. 1002

65

Как видно из полученных результатов, величина изгибающегомомента с учетом угла падения пласта незначительна и при расчетах ее влиянием можно пренебречь. Что касается влияния угла падения пласта на величину осевой силы, то она заметно уменьшается с увеличением угла падения. В рассматриваемом примере это уменьшение составляет 6,2%.

Как видно из рис. 15, при наличии муфты в сечении А-Ау резьбовое соединение может разрушиться от осевой силы, вы­ зывая последующий отвод колонны. Поэтому для сечения А-А прочностной характеристикой является сопротивляемость резьбы на растяжение. Для определения страгивающей нагрузки в резь­ бовом соединении в сечении А-А можно использовать формулу Ф. И. Яковлева.

Длина колонны между сечениями В-В и С-С будет претер­ певать продольный изгиб от сжимающей силы. Так как в се­ чениях В-В и С-С колонна имеет крепкую связь с окружающей породой, то эту длину колонны можем приблизительно рассмат­ ривать как стержень, работающий на продольный изгиб с за­ щемленными концами.

При длинных фильтрах опасной зоной будет участок колон­ ны, заключенный между сечениями В-В и С-С, так как в этой зоне после выноса определенной массы песка из призабойной зоны создаются более благоприятные условия для продольногоизгиба.

Продольная устойчивость фильтровой части обсадной колонны

Продольная устойчивость стержня с различными концевыми условиями достаточно освещена в технической литературе. Существующие методы решения продольной устойчивости стерж­

ня можно с успехом применять при

решении задач, связанных

с продольной устойчивостью обсадной колонны.

Определим критическую силу для

фильтровой части колонн

с учетом перфорированных отверстий. Для определения крити­ ческой силы воспользуемся энергетическим методом.

Уравнение

изогнутой

оси

стержня возьмем в следующем

виде:

 

 

 

 

 

 

 

y =

f sin2^

,

(104)

где / — стрела

прогиба;

/ — длина перфорированного участка

обсадной колонны.

 

 

 

 

Уравнение (104) соответствует условию, когда оба конца

стержня защемлены, при

этом граничные условия следующие:

 

при х = 0

у = 0,

у' =

0;

 

 

 

 

 

(105)

 

при х = I

у = 0,

у' =

0.

66

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

U — A=± 0,

(106)

где U потенциальная энергия

деформации изгиба

стержня,

равная

 

 

о

 

 

А — работа внешней силы Р в

процессе деформации

стержня

в вертикальном перемещении

 

 

о

При определении критической силы возьмем неблагоприятное положение отверстий для устойчивости колонн

(рис. 19).

Предположим, что жесткость ко­ лонны, где имеются перфорированные отверстия, уменьшается на величину АВ(х), тогда жесткость в любом сечении, лежащем на расстоянии х от нижнего конца, будет равна

где

В — жесткость

сечения колонны

без

отверстия.

энергия изгиба при

 

Потенциальная

этом

о

или, учитывая выражение (104),

о

Определим А

о

х

г~1—-

nojnn

i

Рис. 19. Схема определения критической осевой силы для фильтровой части об­ садной колонны с учетом перфорированных отверстий

U

3* 67

Подставляя значения U и А в уравнение (106), для Р полу­ чаем

 

 

 

Р =

4п*В

 

 

 

2 (* ДБ (х)

о 2ях dx

(Ю7)

 

 

 

 

 

Р

 

 

 

/

J

В

 

I

 

 

Как

видно

из

уравнения

(107),

при

наличии

отверстий на

трубе критическая сила уменьшается на величину

 

 

 

 

 

 

 

 

 

i

 

 

 

» 2 ях ,

 

 

 

 

 

 

 

а =

 

2 Г &В (х)

 

 

 

 

 

 

 

— \

---- cos2*----------dx.

 

 

 

 

 

 

 

 

 

I

о.)

В

 

l

 

 

 

Или,

интегрируя по участкам фильтра, получим

 

 

 

 

 

ДБ (х)

 

 

 

 

 

М-р

 

 

 

 

 

 

 

 

 

l ^ L d x + Г

в

 

+

 

 

 

 

 

 

 

 

 

i

 

 

J

 

I

,

2А,+р

 

 

* 2 лх

 

,

 

2 (М~Р)

 

 

2ях

I

С

ДВ(х)

 

 

 

С

ДБ (x)

2

+

j

— ^ - c o s 2 —

 

d x +

|

--------—

COS2 ------------ b •

 

M-P

 

 

 

 

 

 

 

 

2X+p

Б

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегралы, стоящие на нечетных местах, все уничтожаются*

так

как

на

соответствующих участках £ (х )= 0; тогда получим

 

 

х+р

 

2 2лх 1

2 (Х+р)

 

2 2лх I

 

 

С ДБ (х)

Г*

ДБ (х)

 

 

\

-----— cos2----- dx

-f

J

-----— cos2------ b

 

 

J

Б

 

 

 

/

 

 

 

Б

 

/

 

 

 

 

 

 

 

 

 

 

 

 

2A,-f-p

 

 

 

 

 

 

 

 

 

AT

 

 

я (М-Р)

B

i c o i d

, ,

 

 

 

 

-

f

S

 

 

I

 

 

 

 

 

 

 

 

в

 

i

 

 

 

 

 

 

 

n= 1nX+ (n—1) p

 

 

 

 

 

где N — число отверстий.

 

 

 

 

 

 

 

 

Момент инерции ДВ(х) определяется из рис. 19

 

 

 

 

 

 

АВ (х) = к (Б4 — г4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

где

к — число отверстий

по окружности

трубы (на рис. 19 по­

казано одно отверстие); R — наружный радиус трубы; г — внут­

ренний радиус трубы.

 

 

 

 

 

 

находим ср; тогда

Из соотношения а(х)/2я/? = ф/2я

 

 

 

 

 

АВ (х) =

к (R4 — г4)

а{х)

 

 

 

 

 

 

 

 

8

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

68

Как видно из рис. 18, а(х) меняется в интервале (0, р). Для расчета возьмем устойчивое состояние колонн с минимальной жесткостью, т. е. максимальное значение а(х):

AB(x) = ~KiR1~ r,)

8 2

Учитывая, что £ = л(Р 4—г4)/4, для а получим

 

 

р/с

N

 

п (A,-f-p)

 

 

 

 

 

 

 

2пх

 

 

 

 

пШ S

 

I

~т~ dx.

 

 

 

 

 

 

п<=1пК-|- (п—1) р

 

 

Выполняя вычисления для интегралов, находим

 

 

 

 

 

 

 

 

N

 

 

 

|Л/р + +

 

sin

l

V

cos + [2n (k + p) — p] .

 

 

 

/

 

 

 

 

 

 

 

 

n =

l

 

 

Так как диаметр отверстий по сравнению с длиной фильтра

незначителен,

 

можно написать

 

 

 

 

а ■

I

 

 

N

 

\

 

 

2лШ

N + У ’ соэ— {2п:Х + р> — р]

.

(108)

 

 

 

 

JamJ

 

I

 

 

 

 

 

 

 

п~ 1

)

 

 

Подсчитаем величину а для фильтра длиной 10 м из 168-мм

обсадной трубы с отверстиями диаметром 1,27 см.

1 м примем

Среднюю частоту отверстий по длине трубы на

равной 20. Тогда при длине 10 м будем иметь 200 отверстий. Расстояние К (в см) между отверстиями определится из соот­ ношения:

 

l — Np

_ 1000-200-1,27

: 3,72.

 

N + l

 

200+ 1

 

Абсолютное

значение

2

 

cos — [2/г(^+ р)—р]= 128,55. Под-

ставляя все данные в выражение (108), получим

а =

-----------------(200 + 128,55) « 0,01 к.

 

2-3,14.8,4.1000v

 

Зная а, для критической нагрузки получим

 

 

4л2В

 

 

Лер —

/2

- (1 — 0,01/с).

 

При к= 1, 2, 3, ..., Ркр соответственно

равно 0,99Р; 0,98Р;

0,97Р и т. д.

 

 

 

 

 

 

 

4л2В

(108а)

69

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