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

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

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

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

Добавлен: 19.10.2024

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

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

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

стве, определяемый как среднеарифметическое из расходов воздуха: на входе (С2) и выходе (£?х) из межтрубиого пространства [91].

Преобразование уравпеиий (4.92)—(4.94) с учетом выражений (4.95) и (4.97), известных из термодинамики влажного воздуха,, и (4.98), приводит к следующей системе, из которой нетрудно вы­

числить искомые значения t2, t3, f4

при заданной температуре гх:

if4 = Кх -\- ~\fК\ -}- D1{іх [Е9E3(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[EaE12 (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