Файл: Щербань, А. Н. Прогноз и регулирование теплового режима при бурении глубоких скважин.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 (xn—iyn) — С, (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 |