Файл: Щербань, А. Н. Прогноз и регулирование теплового режима при бурении глубоких скважин.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.10.2024
Просмотров: 104
Скачиваний: 0
стве, определяемый как среднеарифметическое из расходов воздуха: на входе (С2) и выходе (£?х) из межтрубиого пространства [91].
Преобразование уравпеиий (4.92)—(4.94) с учетом выражений (4.95) и (4.97), известных из термодинамики влажного воздуха,, и (4.98), приводит к следующей системе, из которой нетрудно вы
числить искомые значения t2, t3, f4 |
при заданной температуре гх: |
if4 = Кх -\- ~\fК\ -}- D1{іх [Е9—E3(ES |
Е7Е5)] -\-Ег {Е7Ей -\-Es) -|- |
-}- ЕfjE4 — Е 0};
t-2= 7?х -f- E2t<i-j- E 3t^,
t3 = E3t2 T?4.
В выражениях (4.99)—(4.101)
Ел= -tJ ----- 0,5 |
АгЕ^ |
|
|
427c„ |
|
' cp |
|
E* |
0,5й 2 |
|
|
11 + I |
|
|
|
L— 0,5/1 a ( |
-ë=- Eb |
||
2 U |
|||
|
|
^cp |
) ' |
1 +0,5Й2СЕі2—1)
E3
1 -0,5 Л 2 ( я 2 + 1—
7 " [ ^ - ( і - 0 . 5 Л 2) + / '( а : з - ^ ) — ^ 7 - ] - ^ п . з а б
E, _ lp
1 +0.5^2
р1— 0,5А2
5~~ 1 +0.5Л2 ’
Е а ■ £ ( " тс р 4 - |
н |
|
|
|
427 |
|
|
||
Go |
|
0,5(Л1- 4 |
в)]; |
|
^ т= ^ і -[1 + |
||||
Тср |
|
|
|
|
Еа= Ег1(1 - |
0,5 Лх) + 0,542 (1 - Е и )\ |
|||
Es= Д12 (1 - |
0,5Aj) + О М 2 (1 - ^ 12); |
|||
Е10— 1 + 0,5 |
+ |
■ |
гаер,, |
|
Е.1 1 ' |
0,375<?2 (<?— |
1) . |
|
|
|
1c p |
|
|
|
|
|
|
|
|
En — 0,125Gg (g — |
1) |
|
||
|
|
Gcp |
|
|
Kx= —0,5.DX[Ea—E12 (Ee |
E 7Ee)]; |
D |
c p |
|
r l cp4 |
||
|
(4.99)
(4.100)
(4.101)
(4.102)
(4.103)
(4.104)
(4.105)
(4.106)
(4.107)
(4.108)
(4.109)
(4.110)
(4.111)
(4.112)
(4.113)
(4.114)
(4.115)
■Ml
Температура воздушной струи в бурящейся скважине, исходя из выражений (4.99)—(4.101), определяется в следующей последо вательности:
1)температура воздуха на выходе из скважины — по формуле
(4.99);
2)температура воздуха перед бурильным инструментом — по формуле (4.101);
3)температура воздуха иа выходе у призабойной зоны — по формуле (4.100).
Для обратного теплового расчета бурящейся скважины (рис. 22, 6) расчетная зависимость для определения температуры воздуха на выходе из охлаждающего устройства t0 при заданной температуре иа входе в забой t2, полученная исходя из конечно-разностных уравне ний теплового баланса, записанных последовательно для межтруб ного пространства, призабойной зоны и бурильной колонны, имеет впд
- В Д * 2 + |
( - f i - - Е 10) + Е 7Е4-I- а д } , |
(4.116) |
где
Остальные обозначения приведены выше.
При обратной схеме продувки количество воздуха по глубине скважины практически непзмепно, так как давление воздуха в бу рильной колонне ниже давления в межтрубном пространстве, что исключает возможность утечек через резьбовые соединения. В связи с этим расчеты для обратной схемы выполняются только по форму лам (4.85), (4.86) и относящимся к ним расчетным зависимостям.
Г Л А В А 5
ТЕОРЕТИЧЕСКИЕ ОСНОВЫ РАСЧЕТА ЕСТЕСТВЕННОЙ ТЕМПЕРАТУРЫ ГОРНЫХ ПОРОД ПО ТЕМПЕРАТУРНЫМ ИЗМЕРЕНИЯМ В БУРЯЩИХСЯ СКВАЖ ИНАХ
ОПРЕДЕЛЕНИЕ ТЕМПЕРАТУРЫ ГОРНЫХ ПОРОД ПРИ ДЛИТЕЛЬНОЙ ВЫ СТОЙКЕ БУРЯЩ ЕЙ С Я С КВ АЖ И Н Ы
При геологических и геофизических исследованиях объектом измерения и прогнозирования часто является естественная темпе ратура горных пород или геотермический градиент, выражающий зависимость между глубиной залегания горных пород и их есте ственной температурой.
Геотермический градиент является важной технической величи ной, используемой при выполнении тепловых расчетов рудничного воздуха в глубоких шахтах, разработке и осуществлении проектов сверхглубокого бурения, геотермальных систем и т. и. Эта величина определяется главным образом в результате измерения температуры в шпурах и скважинах различной глубины и конструкции, ио в ос новном в глубоких буровых скважинах.
Основной объем геотермических измерений приходится в настоя щее время на долю скважин, бурение которых окончено или прер вано на длительное время по различным причинам. При этом в ре зультате температурных исследований вводятся поправки на иска жающее влияние различных факторов, обусловленных технологи ческими и физическими процессами, имеющими место при бурении и выстойке скважины или являющимися их следствием (см. гл. 1).
Наиболее существенное влияние на качество измерения геотер мического градиента оказывают бурение и циркуляция промывочной жидкости, приводящие к нарушению естественного температурного поля в горном массиве, пересекаемом скважиной.
Однако при достаточно большом времени периода простоя сква жины температура раствора становится равной естественной
8 Заказ 060 |
ИЗ |
температуре окружающих горных пород, что позволяет для определе ния геотермического градиента ограничиться измерением температуры бурового раствора. Требуется вычислить продолжительность простоя скважины, при которой обеспечивается равепство температур внутри и вне скважины (прямая задача), а также рассчитать изменение температуры раствора вблизи стенки скважины в зависимости от заданного времепп простоя (обратная задача).
Аналитическое решение обеих задач можно получить, исходя из системы уравнений теплопроводности для среды, заполняющей скважину, и для окружающего ее бесконечного массива, как это сделано, например, в гл. 2 для температуры в любой точке массива или скважины. Однако решение такой задачи относительно времени простоя, по истечении которого происходит выравнивание темпера тур виутрп и впе скважины, чрезвычайно громоздко. В связп с этим авторами разработана діетодика численного решения с помощью ЭВМ для расчетной плоскости, перпендикулярной к оси скважины п достаточно удаленной от ее забоя. При такой постановке, как было показано выше, любой период простоя можно рассматривать, как первый. Расчеты выполнялись на электронно-вычислительной циф ровой машине. Выполнение програмдіы осуществлялось на несколь ких этапах в зависимости от величины времени выравнивания тем ператур твыр.
Расчет проводился по следующим формулам, являющимся пре
образованным |
решением системы |
(2.57)—(2.62): |
|
|||||||
|
|
|
JV |
|
|
|
|
|
|
|
|
Пш - |
V [ £ W Tnp)]"-1lE W T»biP)] |
(54) |
|||||||
|
|
|
JV |
|
|
|
||||
где |
и |
іпр ■ |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
||
ттш /■_ |
\ _ |
1 — фпр |
|
|
|
|
|
(5.2) |
||
и іп р ^TnpJ ' - |
1 фпр— Рац* [ ( С / ‘2У - 1 ) ( £ |
7 п р - |
1 ) - * п р ] ; |
|||||||
|
|
Фпр = |
ехР (z2Tnp) erfc (z 1/тПр); |
|
(5.3) |
|||||
|
■Unp = |
l — j / ' -^r- erfc |
r - R о |
\ |
(5.4) |
|||||
|
2 VflaTnp |
/ |
||||||||
|
|
|
|
|
|
|
||||
|
u & '= i — / 4 ^ |
erfc |
r - R o |
\ |
(5.5) |
|||||
|
2 Va2rц |
) ’ |
||||||||
|
|
|
|
|
|
|
||||
■^np |
e,p |
|
R о |
- z 2T,np |
erfc |
|
r — R о |
|
(5.6) |
|
Т Щ |
2 |
а2тПр |
||||||||
|
|
|
|
Величины тпр и твыр в выражении (5.1) представляют собой соответственно текущее время (продолжительность) простоя сква-
114
жииы и |
искомое значение времени выравнивания горных пород |
и среды, |
заполняющей скважину. |
Формулы для U[пр (тВЬІр ) аналогичны (5.2)—(5.6). При R —r0 = О
■формула |
(5.2) |
принимает вид (2.112), где |
|
|
|
|
|
||||||||||
|
|
|
|
Фпр = |
e x p |
(г 2т пр) erfc (z2 ] / т пр) ; |
|
|
(5.7) |
||||||||
|
|
|
znp = |
exp (z2xnp) erfc (z2 уЧ пр) = |
cpnp. |
|
|
(5.8) |
|||||||||
Аналогично для |
U%p (твыр) |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
Фоыр = |
exp (z2TBbIp) erfc (z2 ]/т выр); |
|
|
(5.9) |
|||||||||
|
|
|
^выр = |
exp (z2TBbIp) erfc (z2V твыр) = срвыр. |
|
(5.10) |
|||||||||||
Таким образом, |
вычисление |
£ДпР (тпр) и |
|
|
(тВЬІр) |
сводится |
|||||||||||
к вычислению |
функции |
вида |
1—exp (ж2)еііс (ж), |
или 1—exp (ж2) |
|||||||||||||
{1—erf (ж)]. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Исходные |
данные |
для |
расчета: |
|
|
|
|
|
|
|
|||||||
С^іпр (твыр) = |
0.99; |
тпр = 18 ч; |
z = |
0,128; |
R 0 = |
0,08; |
тц = |
6 ч; |
|||||||||
а 2 = 2 9 ,3 -ІО-4; |
|
N = |
1, |
2, |
3, |
4, |
5, |
6, |
. . ., |
300. |
|
|
|
|
|||
Для |
твы = |
1 |
+ |
1500 ч |
функция erfc (ж) |
определялась по |
таб |
||||||||||
лицам. |
Для |
твыр = |
1500 ч |
и |
более |
вычисления |
проводились |
по |
|||||||||
формуле |
|
|
|
|
|
|
|
|
1 |
|
1-3 |
1-3-5 |
|
|
|
|
|
|
erf (ж) = |
1 ----- ^ ( 1 |
|
|
|
|
(5.11) |
||||||||||
|
2хі |
|
(2 * 2)2 |
(2*2)3 |
|
|
|||||||||||
|
|
|
|
|
X у п, \ |
|
|
|
|
|
|
для больших X [24]. Тогда функция
1— ехр (ж2) [1 —erf (ж)]
приводится к виду
1 --------— |
f l -----— + ■1' 3- |
1-3-5 . |
(5.12) |
||
(2*2)3 “Г |
|||||
*]/Д |
V |
2*2 ^(2*2)2 |
|
Программа вычисления состоит из трех подпрограмм:
1) |
вычисление |
функции |
|
|
||
|
|
*7І п р ( т ВЫр ) |
= 1 — exp (z2rBblp) erfc (z ( ] / твыр) |
(5.13) |
||
для |
тВЬ1р < |
1500 ч, z = |
0,128; |
|
|
|
2) |
вычисление той же функции для тВЬ1р ^ 1500 ч, |
z = 0,128 |
(под |
|||
программа 1, табл. 9); |
|
|
|
|||
3) |
вычисление |
функции по формуле (5.1) для |
твыр = 1-^- |
18 ч, |
||
z = 0,128, |
тпр = |
18 ч (подпрограмма 2, табл. 9). |
|
|
S* |
115 |
|
|
|
|
|
|
|
Т а б л и ц а 9 |
|
Подпрограмма 1 |
для вычисления функции (5.1) при времени |
|||||
|
|
|
простоя свыше 1500 ч |
|
|
|
|
00 |
Чт |
02 |
25 |
Умн |
10 |
|
|
01 |
Y |
00 |
26 |
Зп |
10 |
|
|
02 |
Умн |
01 |
27 |
ФР |
00 |
|
|
03 |
З п |
08 |
28 |
Выч 2 |
03 |
|
|
04 |
Чт |
86 |
29 |
УП 1 |
17 |
|
|
05 |
Сл |
04 |
30 |
Чт |
81 |
Исходные данные |
|
06 |
Зп |
04 |
31 |
V |
00 |
|
|
07 |
Чт |
88 |
32 |
Зп |
И |
|
|
0S |
Дел |
08 |
33 |
Чт |
06 |
Адрес |
Величина |
09 |
Дел |
08 |
34 |
Дел |
08 |
|
|
10 |
Выч 2 |
00 |
35 |
Дел |
И |
|
|
11 |
Зп |
09 |
36 |
Выч 2 |
86 |
01 |
Z |
12 |
Зп |
10 |
37 |
Ост |
01 |
02 |
|
13 |
ФР |
00 |
38 |
Выч 1 |
07 |
03 |
е |
14 |
Выч 2 |
03 |
39 |
УП 1 |
41 |
04 |
0 |
15 |
УП 1 |
17 |
40 |
Ост |
00 |
05 |
|
16 |
БП |
30 |
41 |
Чт |
02 |
06 |
0,99 |
17 |
Чт |
06 |
42 |
Сл |
05 |
07 |
|
18 |
Сл |
10 |
43 |
Зп |
02 |
|
|
19 |
Зп |
06 |
44 |
Чт |
00 |
|
|
20 |
Чт |
04 |
45 |
Зп |
04 |
|
|
21 |
Сл |
87 |
46 |
Чт |
80 |
|
|
22 |
Зп |
04 |
47 |
Зп |
06 |
|
|
23 |
Чт |
09 |
48 |
БП |
00 |
|
|
24 |
Умн |
04 |
|
|
|
|
|
Подпрограмма 2 для вычисления функции £71пр (тВыр) прн временп простоя от 1 до 18 ч
00 |
Чт |
02 |
25 |
Зп |
03 |
50 |
Зп |
13 |
|
|
|
01 |
V |
00 |
26 |
Чт |
06 |
51 |
БП |
41 |
|
|
|
02 |
Умн |
01 |
27 |
Сл ф |
09 |
52 |
Чт 11 |
06 |
Нсходпыс данные |
||
03 |
Зп |
20 |
28 |
Зп |
06 |
53 |
Умн |
13 |
|||
04 |
Ост |
00 |
29 |
Чт |
03 |
54 |
Дел |
04 |
|
|
|
05 |
Чт |
20 |
30 |
Выч 1 |
10 |
55 |
Ост |
01 |
Адрес |
Величина |
|
06 |
Умн |
20 |
31 |
УП 1 |
И |
56 |
Чт |
04 |
|
|
|
07 |
ехр |
00 |
32 |
Чт |
07 |
57 |
Выч 1 |
08 |
|
|
|
08 |
Умн |
05 |
33 |
Зп |
06 |
58 |
УП 1 |
34 |
01 |
Z |
|
09 |
Выч 2 |
86 |
34 |
Чт |
86 |
59 |
Ост |
01 |
|||
02 |
ТПр |
||||||||||
10 |
Зп |
25 |
35 |
Зп |
12 |
60 |
Чт |
. 06 |
|||
03 |
ТвыР |
||||||||||
11 |
Чт |
03 |
36 |
Зп |
13 |
61 |
Сл Ф |
09 |
|||
12 |
V |
00 |
37 |
Чт |
04 |
62 |
Зп |
06 |
04 |
1 |
|
13 |
Умн |
01 |
38 |
Сл |
86 |
63 |
Чт |
86 |
05 |
erfc (х) |
|
06 |
0,05000 |
||||||||||
14 |
Зп |
20 |
39 |
Зп |
04 |
64 |
Зп |
04 |
|||
07 |
0,05000 |
||||||||||
15 |
Ост |
00 |
40 |
Зп |
15 |
65 |
Сл |
И |
|||
08 |
300 |
||||||||||
16 |
Чт |
20 |
41 |
Чт |
15 |
66 |
Зп |
И |
|||
09 |
0,00100 |
||||||||||
17 |
Умн |
20 |
42 |
Выч 1 |
86 |
67 |
Выч 1 |
10 |
|||
10 |
19 |
||||||||||
18 |
ехр |
00 |
43 |
Зп |
15 |
68 |
УП 1 |
34 |
|||
И |
1 |
||||||||||
19 |
Умн |
05 |
44 |
УП 2 |
52 |
69 |
Ост |
00 |
|||
12 |
1 |
||||||||||
20 |
Выч 2 |
86 |
45 |
Чт |
25 |
|
|
|
|||
|
|
|
13 |
1 |
|||||||
21 |
Зп И |
06 |
46 |
Умн |
12 |
|
|
|
|||
22 |
Ост |
01 |
47 |
Зп |
12 |
|
|
|
|
|
|
23 |
Чт |
03 |
48 |
Чт |
13 |
|
|
|
|
|
|
24 |
Сл |
86 |
49 |
Сл |
12 |
|
|
|
|
|
116