Файл: Щербань, А. Н. Прогноз и регулирование теплового режима при бурении глубоких скважин.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.10.2024
Просмотров: 91
Скачиваний: 0
А, |
4 заб^'заб (1 |
£2 ) . |
|
Gcp |
’ |
||
|
Ал _ k-J'ст (1 — Sa)
Gcn
(3.192)
(3.193)
Без учета эксцентриситета бурильной колонны решение системы (3.179)—(3.181) имеет следующий вид:
1(1 —0,5.13.) + 0,5,1., ( £ ,+ £ ,і)
1— 0,5,1! (Ei + Е3— 1) t 3= E j , 4- Eo',
h = Е зк + E i,
где
E _ |
1 -0 ,5 A, |
|
1 |
1+0,5 А 2 ’ |
|
,i , |
. 1 |
Озаб |
А 2^п. Заб ~Г |
/ 1 |
|
ТР _ |
|
(*СР . |
2 1+0,5Л2
Ея
[1 — 0,5 (Лг+ Л3)1 -+0,5Лг . 1 + 0,5 И 1 + Л3)
рі?
4 1 + 0,5 (А! + Л3) *
) — Е 2 [0,5 (А±-f- П3) — 1] + 0,5A-J.^,
Л _ |
^трТ'тр . |
|
A l |
Gc |
' |
л |
^ Ч заб ^ заб . |
|
+*о |
/ч |
, |
|
Gc |
’ |
А |
Ст^Ст |
■ |
Л з |
Gc |
(3.194)
(3.195)
(3.196)
(3.197)
(3.198)
(3.199)
(3.200)
(3.201)
(3.202)
(3.203)
(3.204)
При выполнении расчетов по формулам (3.182) — (3.184) и (3.194) —
(3.196) для периода бурения принимается |
А 2 = 0; А 3 = 0, для |
периода промывки Q3a6 = 0. |
|
Для оценки влияния принятого при |
составлении исходных |
уравнений (3.179)—(3.181) допущения о линейном характере изме нения температуры промывочной жидкости но глубине скважины на точность конечных результатов выполнено сопоставление резуль татов, полученных для периода промывки скважины по формулам
(3.182)—(3.184) и (3.194)—(3.196) и по аналитическим зависимо стям (3.146), (3.147) и (3.157), (3.158). Результаты расчетов при ведены в табл. 6.
На рис. 19 показано изменение температуры промывочной жидкости (глинистого раствора) для скважины глубиной 2000 м (рис. 19, а) и 3000 м (рис. 19, б), вычисленное по формулам (3-146), (3.147) (кривые 3 и 4) и (3.182)—(3.184) (кривые 1 и 2). Расчеты
6* |
S3 |
выполнялись при следующих исходных данных; |
d = |
0,155 м; |
D = |
||||
0,28 м; к = 58,2 Вт/м2 • °С; а = |
0,0364 °С/м; ех = |
0,361; е2 == |
0,2. |
||||
Коэффициент нестационарного теплообмена определялся по фор |
|||||||
муле (3.171) и |
составил |
для |
скважины глубиной |
2000 м к'т — |
|||
— 14 Вт/м2-°С; |
к'х = |
16,3 |
Вт/м2-°С; для скважины глубиной 3000 м |
||||
К\ = 11,6 Вт/м2 ■°С; |
Ay = |
14 Вт/м2 • °С. |
|
|
|
Рис. 19. Сравнение результатов тепловых расчетов.
Как видно из выполненного сопоставления, упрощенные способы тепловых расчетов для бурящихся скважин относительно небольшой глубины являются достаточно точными.
Таблица 6
Глубина |
|
Температура промывочной жидкости в бурильной колонне |
|
|||||
|
II межтрубном пространстве, в °С, вычисленная по формулам |
|||||||
забоя, м |
_______________________________________________________________ |
|||||||
|
43.182) -(3.184) |
(3.146) |
(3.147) (3.194)--(3.196) |
(3.157) |
(3.158) |
|||
0 |
18,0 |
32,4 |
18,0 |
34,6 |
18,0 |
32,2 |
18,0 |
34,4 |
1850 |
23,7 |
33,5 |
24,1 |
34,9 |
22,5 |
33,3 |
23,0 |
34,6 |
2000 |
18,0 |
18,0 |
||||||
0 |
18,0 |
44,0 |
35,9 |
18,0 |
43,6 |
35,2 |
||
2850 |
37,0 |
45,2 |
41,4 |
43,0 |
36,3 |
44,4 |
40,0 |
42,1 |
3000 |
|
|
МЕТОДИКА РАСЧЕТА ТЕМПЕРАТУРЫ ПРОМЫВОЧНОЙ Ж ИД КО СТИ С УЧЕТОМ КОНСТРУКЦИИ БУРЯЩ ЕЙСЯ С КВАЖ И Н Ы
Изложенные выше аналитические методы расчета температуры промывочной жидкости, основанные на решении системы двух дифференциальных уравнений теплового баланса бурильной ко лонны и межтрубиого кольцевого пространства, не позволяют учи тывать переменный характер коэффициента теплопередачи через
84
стенку бурильной колонны и коэффициента нестационарного тепло обмена на стенке скважины, обусловленный изменением диаметра ствола и бурильной колонны по глубине скважины и соответственно изменением скорости движения промывочной жидкости, от которой зависят указанные выше расчетные коэффи циенты.
В табл. 7 приведены значения коэффициента теплопередачи и коэффициента нестационарного теплообмена для различных участков скв. Вул- кановская-3. Как видно из табл. 7, разница в значениях коэффициентов иа отдельных участках равной протяженности достигает одного порядка. Очевидно, что усреднение расчетных коэффици ентов по глубине скважины при расчете тем пературы промывочной жидкости, если рассмат ривать скважину как выработку с равновеликим диаметром, равным некоторому среднему, а днаАгетр бурильной колонны также считать неизмен ным, влечет за собой внесение определенной
погрешности в результаты расчетов. |
|
|
Рис. |
20. |
Схема |
||||
Для более |
точного |
определения температуры |
|||||||
промывочной |
жидкости |
скважину |
необходимо |
разбивки |
скважи |
||||
ны на участки для |
|||||||||
разбивать |
на |
участки |
с |
неизАіенньш диаметром |
выполнения теп |
||||
ствола или бурильной колонны (рис. |
20), |
а тем |
ловых |
|
расчетов |
||||
пературу |
вычислять |
последовательно, |
начиная |
с учетом конструк |
|||||
с устья скважины, для каждого из этих участков. |
ции |
скважины. |
|||||||
|
|
|
Таким образом, фактически для каждой скважины необходимо решать систему дифференциальных уравнений теплового
баланса, содержащую количество уравнений, равное удвоенному числу расчетных участков, причем граничные условия должны быть записаны
Расчетный
участок
Устье
скважины
•1
2
3
4
5
Глубина конца участка *, м |
Коэффициент не стационарного теплообмена, &т, Вт/м2•°С |
Коэффициент теплопередачи, Вт/м’ •°С |
0 |
_ |
|
|
|
|
150 |
6,03 |
218 |
2200 |
5,6 |
754 |
2400 |
15,6 |
І490 |
3800 |
20,6 |
950 |
4800 |
72 |
2430 |
Таблица 7
Температура промывочной жидкости в конце участка*, в °С, при расчетах
с усреднени |
на ЭВМ |
по методу |
|||
ем А |
и В |
|
|
пени моделей |
|
и |
и |
Ц |
U |
и |
1, |
18 |
23,6 |
18 |
25,4 |
18 |
25,4 |
— |
— |
19,2 |
26,7 |
— |
— |
— |
— |
78,6 |
86,1 |
— |
— |
87,5- |
93,5 |
87,3 |
94,7 |
87,3 |
94,7 |
— |
— |
122,8 |
128,3 |
— |
— |
133,7 |
133,5 |
151,9 |
152 |
160,1 |
160,1 |
* «Конец участка» отсчитывается по направлению от устья к забою скважины.
8S
отдельно для каждого участкаПри общем числе участков N в сква жине глубиной II граничные условия такой системы имеют вид для первого участка
|
|
|
tx — tu |
при h ~ 0; |
(3.205) |
|||
|
|
|
— |
|
при h = h1 = l1; |
(3.206) |
||
для |
второго |
участка |
|
|
|
|
|
|
|
|
|
|
, |
при /г = |
|
(3.207) |
|
|
|
|
*2. = *12 + |
^ 2 |
ПРИ |
h = h.1 = l1-)-l2-, |
(3.208) |
|
.для любого |
t-того участка |
|
|
|
|
|||
|
П,-= |
при h = Ѵ і = Іл + |
1, -1- . . . + 1,.г\ |
(3.209) |
||||
|
t2i — tx.-{-Atl |
при |
h = |
/г,- = ^ -f- /о-І-. . . -{- lt; |
(3.210) |
|||
для |
последнего |
участка, |
который |
оканчивается забоем |
скважины, |
|||
|
|
|
h N = hN^ |
при |
h = H — lN] |
(3.211) |
||
|
|
|
hN = hN + At3a6 |
при |
h = H. |
(3.212) |
В выражениях (3.205)—(3.212) величина At представляет собой перепад температур между потоками промывочной жидкости в бу рильной трубе и затрубиом пространстве. Таким образом, точный расчет температуры в бурящейся скважине, который ввиду большого ■количества дифференциальных уравнений и относящихся к ним
граничных условий |
выполняется с помощью ЭВМ, заключается |
в нахождении таких |
значений функций tx (II), 12 (II) и At (II), |
при которых соблюдается равенство (3.212), являющееся контроль ным при заданных исходных параметрах. Такой расчет неудобен из-за больших затрат машинного времени.
Ниже излагается методика, позволяющая упростить выполнение тепловых расчетов бурящейся скважины без заметного ущерба для их точности. Это достигается заменой машинного расчета указанной выше системы с граничными условиями (3.205)—(3.212) расчетом цепи моделей расчетных участков по аналогии с методикой, пред ложенной в работах Б. И. Медведева для тепловых расчетов горных выработок в глубоких шахтах и основанной на следующих сообра жениях.
Приращение температуры потока ät на элементарном отрезке горной выработки dL может быть определено по уравнению
dt = ^ ( t u - t ) d L + ^ - . |
(3.213) |
где кх — коэффициент нестационарного теплообмена; |
G — расход |
жидкости или воздуха; U — периметр горной выработки; с — тепло |
емкость; гп и t — соответственно температура неохлажденного
■86
горного массива и температура потока на отрезке dL\ q' — тепло выделения от абсолютных истопников тепла на единицу длины выработки, отнесенные к единице веса протекающей по ней среды.
Приращение температуры иа элементарном отрезке модели горной выработки можно выразить уравнением
|
|
|
к |
к и к { к , |
k J J (i |
|
Л J r |
1 W b |
<ndL |
(3.214)' |
|
|
|
|
|
i |
|
||||||
|
|
|
|
к а к с |
Gc Ѵ'П |
1>а Ь 1 |
к е |
с |
|
||
где |
|
|
|
|
и ы • |
|
|
|
^П. м . |
|
|
1, |
|
м |
, |
“-I7 |
к |
|
|
|
|||
|
К |
|
и |
|
^ в |
|
^п |
|
|
||
|
|
|
|
|
|
|
|
||||
|
|
|
|
кв = |
Gm . |
|
см . 1. _9м |
|
|||
|
|
|
|
|
G ’ |
|
С ’ |
КЧ - |
д |
|
|
масштабы моделирования величин, |
входящих |
в уравнение |
(3-214),. |
а кХы, и ы, ta, tn м, Іы, GK, см, qu — значения этих величин для модели. Уравнение (3.214) определяет не только приращение темпера туры потока в модели dtM, но и приращение ее dt на элементарном
отрезке выработки.
Таким образом, цепь моделей подбирается в соответствии с рией подобия, чтобы всю ее можно было рассматривать как выработку.
Рассмотрим приращение температуры на элементарном отрезкескважины, которое для несжимаемой жидкости определяется систе мой уравнений, например типа (3.42), (3.43).
Для модели бурящейся скважины система уравнений, аналогич ная (3.42), (3.43), имеет вид
dt1 — |
kkjjkfki kjpUjp |
(3.215)' |
(ta — tj) dh; |
||
duM= kt dtz= |
h g r p . (t2 - tl) dh - |
|
kkjrkfki |
h tj |
(3.216), |
------F/T— |
(ln„-г oft — h)dh, |
где
|
ч |
, . |
к |
= |
’ |
|
к т |
|
к |
Uтр м ■• |
|
/1,тр и |
Р т р |
|
|
||
|
k l |
І и ; |
|
|
1 |
. |
. |
чзЗ К II |
|
к .
_
k G -“
м.
дтр
м_
*1
Gm .
G ’
k j j
^2М
t o
K =
= |
II |
|
^п0 м
гп0
СМ
С
Коэффициенты к, |
ku,t,i,G,o K p и |
kTpU в уравнениях (3.215),. |
(3.216) представляют |
собой масштабы |
моделирования физических |
величин кх, ктр, U, С/тр, t±, t2, tn<>, l, G, с, входящих в эти уравне ния, причем те же величины с индексом «м» относятся к модели.
8 Т