Файл: Щербань, А. Н. Прогноз и регулирование теплового режима при бурении глубоких скважин.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 19.10.2024

Просмотров: 109

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Решение (6.17) ищем в

виде

 

 

 

СО

 

 

Ux =

S

\r I е{»ч>;

(6.21)

 

П=1

 

 

 

СО

 

 

г/2 =

S ^ 2

I'■ I е-(Пф.

(6.22)

 

п=1

 

 

Подставляя последовательно (6.21) и

(6.22) в уравнение

(6.17),

п о л у ч а е м дифференциальные уравнения

для

F 1 н F 2

 

FFl + r F l - F ,

г - b n

=

0;

(6.23)

 

FFl rF’F.y Iir-bn +

li11=

0.

(6.24)

Уравнения (6.23) и (6.24) представляют собой уравнения Бесселя, решения которых выражаются через модифицированные функции

Бесселя

FX= C1Jn(—i ]/і ynr) -f C2Kn

ynr);

 

 

 

 

(6.25)

где

 

Л = C3Jn (УІУпг) + CxK n (VIy,tr),

 

(6.26)

 

 

 

__

 

 

 

 

 

Vn= V K = ) f ~ .

 

 

(6-27)

Учитывая,

что Jn и K n

выражаются через функции

Кельвина

п-го порядка:

 

 

 

 

 

 

 

J n iV i Уп') =

Ьег„ (упг) і Ьеіл (ynr)

 

 

 

 

Jni—і У І упг) = Ьег„ (ynr) i bei„ (ynr)

 

^ 2g^

 

 

K n (Vi ynr) =

ker„ (ynr) + i kei„ (ynr)

 

 

 

 

K n (—i V i y nr) = ker„ (ynr) i kei„ (У,У

 

 

можем записать (6.25) и (6.26) в виде

 

 

 

F1=

Су [Ьегл (упг) — i bei,, (у„г)] -f С2[ker„ (уnr) — i kei„ (yr/)];

(6.29)

F2 =

C3[ber„ (ynr) + Ьеіл (у„г)] +

C4 [ker/t (ynr) —kein (y„?■)].

(6.30)

Тогда

частные решения

(6.21)

и (6.22)

получим в

виде

 

Uy = (cos ?гср

i sin /гср) {Сх[Ьег„ (упг) —i bein (у„г)] + С2[ker„ (улг )+

 

 

+

i kei„ (у„г)]};

 

 

(6.31)

и2=

(cos rap — i sin rap) {С3[Ьегл (ynr)

і Ьеі„ (ynr) ]

 

 

 

 

+ С4 [ker„ (ynr) + kein (Y»?')]}.

 

(6.32)

І28


В (6.31) и

(6.32) учтено, что

 

 

 

 

 

 

 

eimp _

C0S ?г(р

i sin /гср;

 

 

 

 

 

с-іп(р _

cos

Igin пф.

 

 

Решением,

удовлетворяющим

(6.17),

будет

 

 

и = их+ и3 =

(cos іщ + і sin ?гср) {Cx [ber„ (ynr) — i ber„ (у„г)] 4-

 

+ Co [ker„ (ynr) i kei„ (y,/)]} + (cos «rp — i sin ncp) {C3[ber„ (ynr) -|-

 

І bei„ (y„r)] 4- C4 [ker„ (ynr) +

i kei in(y„r)]}.

(6.33)

Введем

обозначения

 

 

 

 

 

 

 

 

bern(ynr) = xn;

 

 

 

 

 

 

bein(y„r)=y„;

 

 

 

 

 

 

ker„ (ynr) = zrl;

 

 

 

 

 

 

kein (Упг) = V

 

 

 

Тогда

(6.33) примет вид

 

 

 

 

 

и = (cos 7кр 4- і sin шр) {С4 (xniyn) — С, (Zn—іт„)} -f-

 

 

-г (cos„cp— i sin ncp) {C3 (x„+ iyn) -'г C4 (zn +

іт„)}.

(6.34)

Удовлетворяя

условию

(6.18),

получаем:

 

 

{С4 [хп (В0) ііу'п(Д0)] + С2 [z'n(Д0) -f ітп’ (Д0)]} (sin na i cos na +

i) -f

 

 

4- {C3 [x"n(Д0) + іУп(7?o)] 4- Ct [z'n(Д0) 4-

 

 

 

 

4- ix'n(Д0)]} (sin na 4- i cos na —i)

.

(6.35)

Из граничного условия

(6.19)

получим:

 

 

{Сг [х'п(Д0) —іуп’ (Д0)1 4- Со [z'n(Д0) —іСп (Д0)]} cos na —sin na i) 4- 4- {C3 [Xn(Д„) 4- iy'n(До)] 4- c t [z'n(Д0) + гт; (Д0)]} (i —sin na i cos na) =

=—h [{Ci [xn(Д0) iyn(Д0)] 4- C3 [z„ (Д0) ixn(Д0)]} X X (г cos na —sin na i) 4- {C3 [xn(Д0) 4- iyn(Д0)] 4-

4 -C4 [z„ (Д0) 4 -ітл (Д0)]} (i —sin na — z cos na)].

(6.36)

Из граничного условия (6.20) получим:

{Сі [хп’ (r0) — iy'n(/•„)] 4- Co [z'n (r0) —ix'n(r0)]} (cos ncp 4- i sin ncp) 4- + {C3 [x'n (r0) 4- iy'n (r0)] 4- C4 [z’n (r0) 4- ix'n (?■„)]} (cos 7zcp — i sin ncp) =

= — h [{Ci [xn(r0) 4- ІУп('‘о)] + CoJza (r0) — ixn(r0)]} (cos ncp 4- i sin cp)4~ 4- {Cg [xn(r0) 4- iynM + Ci [zn(r0) 4- *T„(r0)]} (cos ncp — г sin ncp)]. (6.37)

9 Заказ 660

129-



Введем

обозначения:

 

 

 

 

 

 

 

 

 

хп(Д0)

іУп (До) =~ — (Д0);

 

 

 

 

х п ( R 0) - f

iyfi

(R 0) =

А

4

(/?„);

 

 

 

 

^ ( д 0) - і т „ ( д 0) = д ' - ( Д о ) ;

 

 

 

 

(До) +

г'в> (Д 0) =

Д '

-ь (Д 0);

 

 

 

 

(Д 0) -

і у п(Д 0) =

^

— (Д 0);

 

 

 

 

х л (До) _г гі/« (До) = -4 н • (До);

 

 

 

 

4 (До) — ітя (Д0) = Д — (Д0);

 

 

 

 

-л (Д 0) +

':тл (До) = Д +

(Д 0);

 

 

 

 

і cos /га — sin /га — г = а;

 

 

 

 

г г cos /га — sin /га = а;

 

 

 

а; (/•„) cos /гср

^ (/■„)sin /гср =

Д ' — (/■„);

 

 

 

До) cos /гср -f х/л (/’о) sin /гср = D -f- (r0);

 

 

4

(r0) sin /гср — y'„ (/•„) cos /гср = Dh—(/•„);

 

 

4

До) sin /гср — //„ (/•„) cos /гср = D (/•„);

 

 

 

(r0) cos /гср0 -f Yn (/■„) sin /гср =

AT 4

Д0);

 

 

z’n(/•„) sin /гср — T„ (/•„) cos /гср =

M" (r0);

 

 

4

До) cos /гср -

t„ (/•„) sin /гср =

M

(/•„);

 

 

4

(/•„) sin /гср t„ (/■„) cos /гср = M (/■„).

 

После преобразований получим

 

 

 

 

 

Сз [ Д ' — (/'о) 4 4 P — Д о )] — ^2 (^2 — С 4) Д / ' — До) ~h сѴ/ — (г~о) 1

 

 

 

D ' — До) 4 4 - D — До)

 

 

 

 

=

С3- к 2(С2- С ^

Л1'-(г0)+М-(го)

(6.38)

И

 

 

 

D ' — Д о) -t- HD о)

 

 

 

 

 

 

 

 

 

 

2С3 [D'

До) +

D 4 До) ^2 ! — ^ 2

2 — Д4)

Л/ ' - ( , . 0) + л / - ( г - о )

 

■0' + До)+М >-До)

A

 

X [Д ' 4

(r0) -

h2D 4 (r0)l -

(C24 c 4) [M" Д0) -

/г,М - (r0)].

(6.39)

Для того, чтобы решение (6.32) было действительным, положим:

 

 

Cx — f — ix;

C2 = £ — i£;

 

(6.40)

 

 

C3 — f-\-ix\

Ci =

t,A- г|.

 

 

 

 

 

130


Тогда

 

СI -{-Со .

j-

 

 

 

t

 

 

 

 

У—

2

^ ~

2

 

® =

« (С, —С3);

1 =

г(С ,-С ,).

(6‘41)

С учетом последних соотношений решение будет иметь вид (6.21),

(6.22). Тогда общее

решение

будет

 

 

 

 

 

 

 

ОО

 

 

 

 

 

и =

Ъ и п.

 

 

(6.42)

 

 

 

 

п = 1

 

 

 

Вычислив U по формуле (6.42), можно с помощью (6.16) перейти к температуре в любой точке по толщине стенки бурильной трубы.

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕМПЕРАТУРНОГО ПОЛЯ КОНСТРУКЦИИ С КВ АЖ И Н Ы И ГОРНОГО М АССИВА ПРИ РОТОРНОМ БУРЕНИИ

Как указывалось выше, аналитически невозможно получить замкнутое решение сопряженной задачи о температуре бурильной и обсадной колони и массива в процессе бурения для общего случая (6.3)—(6.10). Полученное аналитическое решение (6.21)—(6.42) вы­ полнено при определенных упрощающих допущениях и требует работы на ЭВМ для выполнения расчетов. В связи с этим темпера­ турное поле конструкции бурящейся скважины и горного массива при роторном бурении применительно к типовой конструкции сверх­ глубокой скважины исследовалось методом электромоделпрования конечпо-разностиым выражением уравнения теплопроводности и за­ коном Кирхгофа для цепи, сходящейся в узел. Преимуществом этого метода является возможность полного воспроизведения в элек­ трических величинах-аналогах геометрических и физических пара­ метров данной модели, задание любых, в том числе смешанных, граничных условий и получение достаточно точных результатов дискретно во времени и пространстве.

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

Математическое моделирование поставленной задачи проводи­ лось на электроинтеграторе Киевского государственного универ­ ситета им. Т. Г. Шевченко, представляющем собой уникальную конструкцию универсальной электрической модели сеточного типа на 1200 узлов и 2400 магазинов сопротивлений, дающих возможность

9*

131