Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Теплопередача в скважинах

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
15.09 Mб
Скачать

что радиоактивные элементы распространены даже глубже земной коры.

Вкорне изменили представления о тепловой историй Земли новые

космогонические взгляды,

основоположником

которых

был акад.

О. Ю. Шмидт. Согласно

этим представлениям

Земля

и планеты

образовались путем постепенного объединения холодных частиц газопылевого протопланетного облака, вращавшегося вокруг Солнца.

Земля, образовавшаяся в результате объединения частиц, была вначале сравнительно однородным телом с одинаковой по объему концентрацией радиоактивных элементов, Стадия равномерного

распределения

источников

 

 

 

 

имела место примерно 4 млрд,

 

 

 

 

лет назад. В тот период

 

со­

 

 

 

 

держание радиоактивных эле­

 

 

 

 

ментов в несколько

раз

пре­

 

 

 

 

вышало

теперешнее,

и

по­

 

 

 

 

этому тогда

Земля

 

в

целом

 

 

 

 

разогревалась.

В определен­

 

 

 

 

ный

момент

создались усло­

 

 

 

 

вия, благоприятные для диф­

 

 

 

 

ференциации вещества в верх­

 

 

 

 

них слоях Земли.

Этому мо­

 

 

 

 

менту

соответствует

 

макси­

 

 

 

 

мальная температура

в верх­

 

Радиус Земли, км

 

ней

части оболочки

Земли.

 

 

 

 

Происходило

формирование

В

'

Е

 

земной коры,

и перераспре­

 

деление

радиоактивных

эле­

Рис. 14. Современное тепловое

состояние

ментов,

приведшее

к

повы­

 

 

Земли:

 

шенной

их

концентрации

в

1 — температура

 

плавления; 2 — распределение

земной коре. Генерация

теп­

температуры по радиусу Земли; 3 — зоны плава­

ния

вещества Землп

 

ла

в ней

возросла

 

в

 

не­

 

 

 

 

сколько

раз.

приведено

результирующее

распределение

темпера­

На

рис. 14

туры, характеризующее термическое состояние Земли в настоящее время. Как видно из графика, в температурном поле Земли можно выделить область аккумуляции тепла и область оттока тепла к по­ верхности. В глубоких недрах и сейчас происходит разогрев вещества, однако этот процесс может не оказывать влияния на тепловой поток, наблюдаемый у поверхности.

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

Заметим, что приведенная схема распределения температуры по радиусу Земли позволяет предположить наличие локальных магма­

тических очагов

в Земле на сравнительно небольших

глубинах

(до 700 км). Эти

очаги не являются сплошными поясами

в Земле

41

и могут создавать определенную тектоническую картину, локализо­ ванную на небольшой территории.

Гипотеза теплового режима Земли, основанная на представлениях

оее происхождении путем постепенного объединения холодных час­ тиц 1азопылевого протопланетного облака, удовлетворяет сведениям

отеплооттоке с поверхности Земли, содержании радиоактивных

элементов, возрасте Земли и т. д. Эта гипотеза наилучшим обра­ зом объясняет многие явления природы (землетрясения, тектониче­ ские перемещения, вулканическую деятельность).

§ 2. ТЕМПЕРАТУРНЫЙ РЕЖИМ ГЕЛИОТЕРМОЗОНЫ

Под гелиотермозоной понимают внешнюю оболочку Земли, из пределов которой за геологическое время происходит отток тепла к поверхности и в толще которой на термическом режиме, а также на физических свойствах породы и вод сказывается влияние Солнца в любой степени и форме.

По данным геофизики принимается, что мощность этой зоны ориен­ тировочно равна 200—1000 км. С точки зрения проходки и эксплуа­ тации скважин наибольший интерес представляют слои земной коры мощностью до 20—30 км.

При описании температурного поля гелиотермозопы принимают во внимание не только решающее влияние солнечной радиации на термический режим пород, но и любое проникновение температурных волн с поверхности Земли.

В отличие от искусственного температурного поля, формирующе­ гося под воздействием человека, температурный режим гелиотермо­ зоны является естественным.

Для верхней подзоны, где геотермическая зональность согласуется с климатической, при характеристике естественного температурного поля учитывают в основном закономерности распределения радиа­ ционного баланса у поверхности Земли и распределение температур по вертикали, определяемое геотермическим градиентом.

Температурные условия в пункте наблюдения характеризуются среднесуточными, среднемесячными и среднегодовыми температу­ рами. Значения среднесуточных и среднемесячных температур пре­ терпевают изменения в различные годы, в то время как среднегодо­

вая температура

для данного пункта изменяется не существенно

(не более десятых долей градуса).

Анализ естественного регионального температурного поля сво­

дится к

решению

дифференциального уравнения теплопроводности

и учету

влияния

основных для данного региона географических

игеологических факторов.

Вобщем случае температурное поле гелиотермозоны нестацио­ нарно, но расчеты с достаточной степенью приближения можно вести для стационарных условий, учитывая, что изменения темпера­ туры для рассматриваемого периода времени и сравнительно неболь­ шой глубины ничтожно малы.

42

Обычно исходят из наиболее простой физической модели: имеет место стационарный теплообмен между двумя полупространствами, одно из которых (например, поверхность Земли) находится при по­ стоянной температуре to, а другое — нижележащие горные породы — имеет температуру t (Z), изменяющуюся в зависимости от расстоя­ ния между источником тепла и рассматриваемым слоем. В этих усло­ виях дифференциальное уравнение теплопроводности имеет следу­ ющий вид:

 

d z2 ■= 0.

 

Поэтому dtldZ = Г и t = r Z -f- А, где А

постоянная интегри­

рования.

,

 

Используем граничные условия первого рода: при Z = 0 t (Z) =

=-- t о. Тогда

постоянная интегрирования А

= to и уравнение для

определения текущей температуры примет вид

 

t = tо -f* r Z .

(И1.1)

В случае слоистости рассматриваемой плоскопараллелыюй среды уравнение (III.1) можно привести к следующему виду:

 

Z

Z

 

2

“Ь ^ Г (Z) dZ = £0 +

Я^ X(Z)

(HI-2)

Если число слоев ограничено, то распределение температуры по

глубине можно определить из формулы

 

 

 

t

 

(Ш .З )

 

2=1

 

 

где п — количество слоев рассматриваемой среды;

AZt — коэффи­

циент теплопроводности и мощность i-того слоя.

 

При выводе формул (III.1), (III.2),

(Ш .З) предполагалось, что

в среде отсутствуют рассеянные источники тепла и поэтому плотность теплового потока постоянна по глубине. Влияние рассеянных источ­ ников тепла на естественное тепловое поле можно оценить, если считать, что радиоактивные элементы в рассматриваемой слоистой среде распределены равномерно [41]

' - « • + ! < iiu )

где qо — плотность теплового потока в верхнем слое гелиотермозоны; х — количество тепла, выделяемое единицей объема породы за счет радиоактивного распада.

Последняя формула показывает, что с учетом тепла радиоактив­ ного распада распределение температуры по глубине зависит не

43

только от теплопроводности горных пород, слагающих разрез, но

исвязано с уменьшением плотности теплового потока с глубиной. Выражения для определения стационарной температуры, предло­

женные^ другими исследователями [42*, 101, 104], получены при рас­ смотрении аналогичной физической модели и принципиально мало отличаются от формулы (III.4).

Использовать полученные аналитическим путем зависимости для определения расчетным путем стационарной температуры в сква­ жинах можно в том случае, если имеются надежные сведения о тепловых потоках вблизи поверхности Земли, о теплофизических свойствах горных пород, о температуре нейтрального слоя данного района и о количестве тепла, выделяемого горными породами за счет радиоактивного распада.

Особые трудности при прогнозировании естественной температуры возникают в тех случаях, когда температурное поле осложнено рядом геологических факторов, учесть которые расчетом не всегда воз­ можно. В таких случаях предпочтительны эмпирические зависимости и непосредственные измерения стационарной температуры в сква­ жинах.

Если предположить, что основной источник тепла расположен значительно глубже забоя рассматриваемой скважины, то, распола­ гая сведениями о градиенте температуры и теплофизических свой­ ствах вблизи поверхности Земли, можно приближенно определить плотность теплового потока, используя для этого закон Фурье

Я= Г<До» (III.5)

где Г о, Я'о -— соответственно геотермический градиент и коэффициент теплопроводности пород вблизи земной поверхности.

Учитывая тесную взаимосвязь между тепловыми потоками, гра­ диентами температуры, теплофизическими свойствами горных пород, и температурными условиями у поверхности Земли, рассмотрим более подробно каждый из этих параметров.

§ 3. РЕЗУЛЬТАТЫ ПРОМЫСЛОВЫХ ИЗМЕРЕНИЙ ТЕПЛОВЫХ ПОТОКОВ И ГРАДИЕНТОВ ТЕМПЕРАТУРЫ

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

Под удельным тепловым потоком из недр Земли понимают потери тепла в единицу времени с единицы поверхности путем теплопровод­ ности. Поэтому определение теплового потока должно быть основано на точных измерениях градиента температуры и теплопроводности именно для тех пород, для которых определен градиент температуры.

Поскольку величина градиента температуры изменяется по глу­ бине не только в результате изменения теплопроводности, но и под

44

влиянием других факторов, то нельзя считать удовлетворительным способ определения теплового потока путем обычного перемножения усредненного градиента но всей рассматриваемой глубине на сред­ нюю для данного района теплопроводность горных пород. Именно поэтому пока не представляется возможным использовать для опре­ деления теплового потока промысловый термокаротаж.

До последнего времени надежными методами определения тепло­ вого потока считают те, которые осуществляют на основе данных точечных измерений температуры по глубине с точностью до сотых долей градуса. При этом необходимо располагать большим числом определений теплопроводности кернов пород, отобранных из различ­ ных горизонтов.

Как в СССР, так и за рубежом проведены тысячи измерений темпе­ ратуры в скважинах и шахтах. Многие результаты обобщены и систе­ матизированы. Однако немногие из них могут быть использованы для определения тепловых потоков, так как при измерениях не всегда использовалась единая методика, применялись малочувствительные приборы, редко исследовались теплофизические свойства горных пород. Поэтому в настоящее время величину теплового потока оце­ нивают только по результатам специальных исследований, которых, к сожалению, проведено очень мало.

Такие специальные работы по определению теплового потока были проведены АН СССР [102] на Балтийском и Украинском щитах, в Средней Азии и на Северном Кавказе (табл. 3).

 

Т а б л и ц а 3

Результаты специальных измерений теплового потока

и геотермического градиента в СССР

 

 

Величина

Величина

 

геотермиче­

Пункт измерения

теплового

ского

потока,

градиента,

 

10-*

1Q-*

 

кал/см**с

градус/м

Украинский щ и т .........................................

0,7—0,8

0,6—1,52

Балтийский щит .........................................

0,72—0,86

0,8-1,42

Средняя Азия (на юге Казахстана) . .

1,2

1,36-2,0

Северный Кавказ (Краснодарский край)

1,8-2,1

3 ,0 -6 ,0

Как видно из табл. 3, в областях древних кристаллических щитов наблюдаются низкие градиенты температуры и малые величины тепловых потоков. Считают, что причиной этому служит обедненность кристаллических щитов радиоактивными элементами. Что касается больших величин тепловых потоков в Предкавказье, то существуют разноречивые мнения относительно причин, обусловивших эту ано­ малию. Возможно, что такое явление связано с наличием неглубо­ кого очага расплавленного вещества.

43

Аналогичные измерения теплового потока, проведенные на юге Северной Америки (рис. 15) и в Японии (рис. 16), показывают, что

Рис. 15. Распределение теплового потока на юге Северной Америки (по данным Roy, 1963).

Цифры показывают величину теплового потока в 10" б кал/см2*с

величина потока тепла из недр к поверхности Земли колеблется в значительных пределах. Аномалии пока не имеют однозначного объяснения, хотя установлено, что они приурочены к районам совре­ менных тектонических дви­

жений.

В настоящее время в боль­ шом объеме проводят иссле­ дования по определению удельных тепловых потоков в различных районах земного шара *.

* Для систематизации резуль­ татов исследований организован Международный комитет по тепло­ вому потоку (IHFC) под председа­ тельством Ф. Берча. Вице-предсе­ дателем комитета избрана Е. Лю­ бимова.

Рис. 16. Карта теплового потока в v Японии (по данным Уеда п Хоран).

Цифры показывают величину теп­ лового потока в 10_б кал/см2*с

46

Обобщение фактического материала позволило выявить средние значения геотермических градиентов и удельных тепловых потоков для различных районов мира (табл. 4).

 

 

 

 

 

Т а б л и ц а 4

 

 

Наибольшая

Геометрический гради­

Средний

 

 

ент, градус/м

Районы исследований

глубина

 

 

 

тепловой

исследова­

 

 

 

поток,

 

 

ния, м

от

до

средний

кал/см2-с

 

 

 

А ф р и к а .........................................

 

3048

0,007

0,022

0,012

1,1

Канада .........................................

 

2134

0,009

0,016

0,013

0,96

США .............................................

 

2743

0,0186 ; 0,039

0,026

1,47

И р а н .............................................

 

914

0,87

Англия .........................................

 

1200

0,0133

0,08

0,034

1,33

Тихий океан .............................

. . . .

1520

0,04

0,214

0,083

1,45

Атлантический океан

1440

0,024

0,055

0,039

0,98

Нетрудно заметить, что величина геотермического градиента по

земному шару колеблется в довольно широких пределах. Не

менее

существенные колебания

геотермического градиента наблюдаются

в пределах СССР. Так,

если его величина на Кавказе

0,04-f-

^ 0,08 градус/м, то в Кривом Роге и Белоруссии градиент температуры не превышает величины 0,011—0,012 градус/м (табл. 5). В пределах Предкавказья и Кавказа геотермический градиент колеблется в до­

вольно широких пределах

(от 0,032 до 0,083 градус/м).

 

 

Т а б л и ц а 5

 

Средняя

Средний геотер­

Район наблюдения

температура

мический гради­

на глубине

ент в интервале

 

1000 м,

100—1000 м,

 

°С

градус/м

Северный Кавказ .................................

90,7

0,083

Дагестан .................................................

55,6

0,047

Западный Кавказ .................................

50,4

0,04

Апшеронский полуостров

................ 47,2

0,036

Западное Предкавказье.........................

41,9

0,032

Донбасс ...................................................

39,6

0,031

Эмба ..........................................................

41,4

0,03

Западная Украина .............................

31,4

0,023

Нижнее П оволж ье.................................

28,6

0,02

Самарская Лука .................................

24,8

0,016

Башкирия .............................................

18,4

0,012

Белоруссия .............................................

23,4

0,012

Камское Приуралье .............................

17,4

0,011

Кривой Р о г .............................................

19,5

0,09

47

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

 

 

 

 

 

мости

от

изменения

теплопровод­

 

 

 

 

 

ности

пород.

Увеличение геотерми­

 

 

 

 

 

ческого градиента с глубиной обычно

 

 

 

 

 

наблюдается

в

местах перехода от

 

 

 

 

 

вышележащих пород

низкого тепло­

 

 

 

 

 

вого сопротивления к

породам с по­

 

 

 

 

 

вышенным тепловым сопротивлением.

 

 

 

 

 

В терригенных

отложениях

геотер­

 

 

 

 

 

мический

градиент

заметно

увели­

 

 

 

 

 

чивается при переходе от песчаных

 

 

 

 

 

пород

к

глинистым.

Наибольшее

 

 

 

 

 

увеличение

 

геотермического

гради­

 

 

 

 

 

ента наблюдается при переходе от

 

 

 

 

 

плотных карбонатных

и гидрохими­

500

1000

1500

2000

Н,м

ческих

осадков

к

рыхлым

терри-

 

 

 

 

 

генным отложениям.

 

 

 

Рис. 17. Изменение усредненных

Некоторую роль в геотермической

значений геотермического градиен­

характеристике

разреза играет плот­

та с глубиной в песчано-глини-

ность

пород

(рис. 17).

При прочих

стых

отложениях:

 

1 — Центральное

Предкавказье;

2

равных условиях с увеличением плот­

месторождение Махачкала; 3 — районы

ности пород геотермический градиент

Эмбы; 4 — Аишеронский

полуостров;

S — Юго-западная

Туркмения;

6 —

уменьшается.

 

 

 

 

 

Астраханское Поволжье

§ 4. О ТЕМПЕРАТУРЕ НЕЙТРАЛЬНОГО СЛОЯ

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

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

Нейтральный слой представляет собой условную границу между верхней зоной переменных температур и нижней зоной стационарных температур. В результате сезонных колебаний температуры термиче­ ское состояние грунта над нейтральным слоем периодически изме­ няется. Изменяется и тепловой поток как по величине, так и но направлению.

Переменное температурное поле в этой зоне не поддается точному описанию простыми аналитическими зависимостями, так как подвер­ жено чрезмерному влиянию процессов передвижения влаги, ее фазо­

'43

вых превращений и т. д. Однако для оценки температур на больших глубинах и выяснения роли влияния поверхностных факторов на геотермическое поле вполне достаточно приближенного его описания. В качестве такого приближения обычно используют экспоненциаль­ ную зависимость для затухания колебаний с глубиной, которая полу­ чена при условии, что теплопередача в породах верхней зоны проис­ ходит исключительно путем теплопроводности, которая принимается неизменной. Для соблюдения этого условия исключают из рассмотре­ ния самые поверхностные слои, теплофизические свойства которых сильно изменяются благодаря сезонным колебаниям влажности.

Таким образом, температуру грунта в первом приближении можно описать следующим уравнением:

 

 

 

 

 

 

 

 

 

(1н*б)

где

Th — переменная температура

на глубине

h в любой

момент

времени

т;

Тн0 — среднегодовая

температура

пород

на

глубине

hQ(ho <

h)\

A Q — амплитуда

колебаний

температуры

на

глубине

h0;

т0 — период циклических

колебаний

температуры.

 

 

 

Наиболее точно положение нейтрального слоя можно определить

путем систематических измерений температуры на различных глуби­ нах в зоне переменных температур. Определение глубины нейтраль­ ного слоя по данным таких измерений можно осуществить графиче­ ским методом. Более приближенно глубину нейтрального слоя опре­ деляют по формуле

(Ш.7)

где А н — допустимая погрешность изменений температуры (обычно принимают А н = 0,1° С).

Отсутствие прямых измерений естественной температуры пород нейтрального слоя в невскрытом массиве приводит к разноречивым мнениям относительно ее величины. Одни исследователи предлагают принимать эту величину равной среднегодовой температуре воздуха, другие — среднегодовой температуре почвы, третьи — величину, в 1,5—2 раза превышающую среднегодовую температуру воздуха.

Для оценки температуры нейтрального слоя следует исходить из среднегодовых температур, определенных для глубин, превышающих глубину сезонного промерзания или оттаивания.

В настоящее время большинство метеорологических станций

СССР проводит измерения температур грунтов до глубины 3,2 м. В большинстве районов эта глубина больше мощности слоя сезонного промерзания, и поэтому температура на этой глубине почти соответ­ ствует температуре нейтрального слоя.

В табл. 6 приведены обработанные данные метеослужбы г. Крас­ нодара о среднегодовой температуре и величине среднего градиента

4 Заказ 1249

49

 

 

 

 

Т а б л и ц а 6

Данные метеослужбы об изменении среднегодовой

температуры

 

в г. Краснодаре за 1967—1968 гг.

 

 

 

Температура, °С

 

 

Воздух

 

Грунт на глубине, м

 

 

0.4

0,8

1,6

3.2

11,7

12,75

12,9

13,08

13,4

 

Градиент температуры,

градус/м

 

2,6

0,375

 

0,225

0,2

температуры в малых интервалах глубин вблизи дневной поверх­ ности. Анализируя табл. 6, убеждаемся, что при данном темпе изме­ нения величины градиента температуры среднегодовая температура на глубине нейтрального слоя окажется на 0,3—0,5° С выше темпе­ ратуры грунта на глубине 3,2 м и составит для г. Краснодара вели­ чину 13,7—13,9° С. Величина градиента температуры на этой же глубине (10—30 м) уменьшится до 0,04—0,05 градус/м, что соответ­ ствует его среднему значению в интервале 100—1000 м.

Что касается соотношения между температурой нейтрального слоя и среднегодовой температурой воздуха, то разница между ними колеблется в сравнительно широких пределах (от 0,8 до 12,5° С). Превышение температуры грунта над температурой воздуха замечено давно, и многие исследователи пытались найти количественную связь между средней многолетней температурой воздуха и температурой слоя постоянных годовых температур. Однако эта задача оказалась неразрешимой по причине непостоянства разности между температу­ рой воздуха и пород для разных районов. Особенно резко отличаются указанные температуры в районе вечномерзлых пород. Если для евро­ пейской части СССР средняя разница между температурами пород и воздуха составляет около 2° С, то для Сибири она равна примерно 6° С. Поэтому использовать среднегодовую температуру воздуха в качестве температуры нейтрального слоя следует с особой осторож­ ностью.

Определенный научный и практический интерес представляет оценка стабильности температуры воздуха и верхних слоев гелиотермозоны с течением времени. Известно, что температура поверхности Земли изменяется во времени не только в результате изменения про­ цессов выделения глубинного тепла, но и под действием климатиче­ ских циклов. Возникает вопрос: нельзя ли по данным метеослужбы судить о кинетике глубинных тепловых процессов? Ведь климати­ ческие циклы приводят к согласованному изменению температуры

50