Геотермическое поле земли
Геотермическое распределение температуры в массиве пород данного региона являестя стационарным и определяется региональным тепловым потоком из недр земли, условием теплообмена на поверхности и теплопроводностью пород.
В отсутствии латерального теплового потока и сильной анизотропии теплопроводности, это распределение может быть записано через вертикальный геотермический градиент, следующим образом (см. вывод): LaTeX Math Block |
---|
| T_g(x,y,z) = T_{ref}(x,y) + \int_{z_{ref}}^z G_T(x,y,z) dz |
где LaTeX Math Block |
---|
| G_T(z) =\frac{d T_g}{d z}= \frac{j_z}{\lambda_e} |
|
вертикальная проекция регионального геотермического градиента | | распределение вертикальной проекции регионального теплового потока из недр Земли на заданной глубине (обычно LaTeX Math Inline |
---|
body | j_e = 25 \div 55 \ mW/m^2 |
---|
|
) | | профиль теплопроводности в массиве горных пород | | опорная глубина, для которой известна стационарная температура пород | LaTeX Math Inline |
---|
body | T_{ref}(x,y) = T_g(x,y,z_{ref}) |
---|
|
| значение опорной температуры |
Anchor |
---|
| T_g_deduction |
---|
| T_g_deduction |
---|
|
Expand |
---|
title | Общий подход к моделированию геотермического поля |
---|
|
В общем случае, геотермическое распределение температуры в массиве пород данного региона является решением стационарного уравнения теплопроводности: LaTeX Math Block |
---|
anchor | nablaT |
---|
alignment | left |
---|
| \nabla {\bf j} = 0 \quad \rightarrow \quad \partial_{\alpha} \, j_{\alpha} = 0
|
где LaTeX Math Block |
---|
| j_{\alpha} = \lambda_{\, \alpha \beta} \: \partial_{\beta} T_g |
| плотность теплового потока Земли | LaTeX Math Inline |
---|
body | \hat \lambda = \lambda_{\, \alpha \, \beta} |
---|
|
| тензор теплопроводности |
Как правило, теплопроводность имеет небольшую анизотропию вдоль и поперек горизонта залегания и в главных осях имеет вид: LaTeX Math Block |
---|
| \hat \lambda= \begin{pmatrix}
\lambda_{\perp} & 0 & 0 \\
0 & \lambda_{\perp} & 0 \\
0 & 0 & \lambda_{||}
\end{pmatrix} |
где | теплопроводность пород вдоль горизонта залегания | | теплопроводность пород перпендикулярно горизонту залегания |
При анализе небольших по площади участков (десятки километров) можно пренебречь влиянием латерального теплового потока и анизотропией тензора теплопроводности и тогда, геотермическое распределение температуры в регионе можно записать как LaTeX Math Block |
---|
anchor | nablaT |
---|
alignment | left |
---|
| \partial j_z = 0
|
откуда LaTeX Math Block |
---|
| \lambda_z(x,y,z) \, \partial_z \, T_g = j_z(x,y)
|
где региональное распределение вертикально компоненты плотности теплового потока Земли.
Интегрируя уравнение LaTeX Math Block Reference |
---|
| получим уравнение LaTeX Math Block Reference |
---|
| .
|
Условия теплообмена на поверхности задаются опорной глубиной и температурой. Понятие опорной глубины и температуры связано с тем, что верхний слой грунта не является стационарным во времени и подвержен дневным и сезонным колебаниям температуры, вызванным переменным влиянием солнца и атмосферных явлений. Как правило, эти влияния проникают не глубже 20 – 30 м. Эта отметка называется глубиной залегания нейтрального слоя, а его температура LaTeX Math Inline |
---|
body | T_{n_0} = T_g(x,y,z_{n_0}) |
---|
| называется температурой нейтрального слоя.Именно ее и выбирают в качестве опорной для расчета температурного поля ниже опорной отметки.
Нейтральный слой
Математическая модель температуры грунта на произвольной глубине и в произвольный момент времени , от момента начала цикла температурных колебаний, может быть решена численно.Однако полезность этих расчетов на коротких временах (сутки и даже месяцы) ограничивается спонтанностью входящих в модель параметров: - волатильность степени солнечного потока в данном регионе в зависимости от облачности и степени поглощения/отражения поверхности земли
- волатильности атмосферных явлений (ветер и осадки)
- волатильности активности грунтовых вод
Для практических целей моделирования температуры пород выше нейтрального слоя вполне достаточно феноменологической корреляции [Kasuda, 1965]: LaTeX Math Block |
---|
| T(t, z) = T_{srf} + \frac{j_z}{\lambda_e} (z-z_{srf}) + T_A \, \exp \bigg[ \, {(z_{srf}-z}) \sqrt{\frac{\pi}{a_e \, \delta_T}} \, \bigg] \, \cos \bigg[ \, 2 \pi \frac{t - t_{min}}{\delta_T} + (z_{srf} -z) \sqrt {\frac{\pi}{a_e \, \delta_T}} \, \bigg] |
где | абсолютная отметка поверхности земли (обычно полагается ) | | средняя температура пород на поверхности (обычно по метеосводкам региона) | | среднее значение амплитуды циклических колебаний на поверхности земли (обычно по метеосводкам региона) | | период циклических колебаний температуры (обычно LaTeX Math Inline |
---|
body | \delta t = 1 \, {\rm год} |
---|
|
) | | смещение от начала цикла при котором отклонение температуры от среднего минимально | LaTeX Math Inline |
---|
body | a_e = \frac{\lambda_e}{\rho_e \, c_e} |
---|
|
| температуропроводность грунта | | плотность грунта | | удельная объемная теплоемкость грунта при постоянном давлении |
Волатильность теплообмена с атмосферой зашифрована в параметры и и определяется по региональным метеосводкам.
Формула LaTeX Math Block Reference |
---|
| применима как к сезонным, так и дневным колебаниям температуры, но с учетом вышесказанных оговорок о волатильности параметров теплообмена на коротких временах наблюдений.
По этой же формуле можно оценить глубину залегания нейтрального слоя (по пороговому значению погрешности температурных измерений): LaTeX Math Block |
---|
| z_{n_0} = z_{srf} + \sqrt{\frac{a_e \, \delta_T }{\pi}} \, \ln \frac{T_A }{\delta T_{err} } |
где | измерительная погрешность термометра (обычно LaTeX Math Inline |
---|
body | \delta T_{err} = 0.1 \, ^\circ \rm C |
---|
|
) |
Expand |
---|
title | Вывод формулы глубины залегания нейтрального слоя |
---|
|
Для вывода формулы LaTeX Math Block Reference |
---|
| запишем значение максимальной и минимальной температуры по формуле LaTeX Math Block Reference |
---|
| в зависимости от момента времени на хронологической шкале цикла:
LaTeX Math Block |
---|
| T_{max} = T_{srf} + \frac{j_z}{\lambda_e} (z-z_{ref}) + T_A \, \exp \bigg[ \, {(z_{ref}-z}) \sqrt{\frac{\pi}{a_e \, \delta_T}} \, \bigg] |
LaTeX Math Block |
---|
| T_{min} = T_{srf} + \frac{j_z}{\lambda_e} (z-z_{ref}) - T_A \, \exp \bigg[ \, {(z_{srf}-z}) \sqrt{\frac{\pi}{a_e \, \delta_T}} \, \bigg] |
откуда получается размах температурных колебаний на глубине : LaTeX Math Block |
---|
| \delta T_{max} = T_{max} - T_{min} = 2 T_A \, \exp \bigg[ \, {(z_{srf}-z}) \sqrt{\frac{\pi}{a_e \, \delta_T}} \, \bigg] |
Если метрологическая погрешность термометра составляет , то он в состоянии зарегистровать колебания свыше LaTeX Math Inline |
---|
body | \delta T_{max} = 2 \, \delta T_{err} |
---|
| то есть: LaTeX Math Block |
---|
| \delta T_{err} = T_A \, \exp \bigg[ \, {(z_{srf}-z_0}) \sqrt{\frac{\pi}{a_e \, \delta_T}} \, \bigg] |
откуда и вытекает формула LaTeX Math Block Reference |
---|
| . |
Геотерма скважины
Геотермическое распределение температуры вдоль траектории данной скважины называется геотермой данной скважины и являестя как функцией геотермического поля региона, так и траектории данной скважины и задается следующей моделью: LaTeX Math Block |
---|
| T_g(l) = T_{n_0}(x,y) + \int_{z_{n_0}}^{z(l)} G_T(x,y,z) dz = T_{n_0}(x,y) + \int_{l_{n_0}}^l G_T(x,y,z(l)) \cos \theta dl |
где | отметка нейтрального слоя вдоль траектории скважины (обычно так как начальные участки скважин не имеют сильного отклонения от вертикали) |
В регионах, где нейтральный слой остается постоянным по площади LaTeX Math Inline |
---|
body | T_{n_0}(x,y) = T_{n_0}= \rm const |
---|
| , а геотермический градиент остается постоянным по площади по глубине LaTeX Math Inline |
---|
body | G_T(x,y,z) = G_T = \rm const |
---|
| , геотермическое распределение температуры в породах принимает простой вид: LaTeX Math Block |
---|
anchor | T_g_const |
---|
alignment | left |
---|
| T_g(l) = T_{n_0} + (z(l) - z_{n_0}) \, G_T \, \sin \theta(l) |
Однако в большом количестве практических случаев это не так и применение среднего по всему разрезу значения геотермического градиента для оценки геотермического распределения температур по формуле LaTeX Math Block Reference |
---|
| может привести к значительным погрешностям.
Наиболее популярные причины отклонения геотермического профиля от LaTeX Math Block Reference |
---|
| являются: - большой контраст теплопроводности по горизонтам
- латеральные потоки флюида (как секулярные так и возникшие в процессе разработки месторождений)
- близкое расположение активного разлома
к котором следует еще добавить техногенные явления в процессе геотермических исследований:
- действующие перетоки в исследуемой скважине
- недостаточный срок стабилизации температуры после останова скважины для геотермических исследований
Наиболее популярными причинами дрейфа геотермы или е отдельных участков по площади являются: - переменный по площади профиль теплопроводности пластов (вызванный особенностями распределения фаций и вертикальными разломами)
- вертикальные перетоки в разломах и скважинах
- переменный тепловой поток Земли вызванный конвективными ячейками в мантии
- переменные условия теплообмена на поверхности земли (например суша-вода)
|