- •Численная реконструкция термической истории осадочных бассейнов
- •Алгоритм решения воспроизводит природный процесс осадконакопления путем наращивания сверху области счета на элемент
- •Уравнение
- •Используется консервативная форма записи уравнения теплопроводности, которая является следствием закона сохранения энергии для
- •Конвективный фактор теплопроводности
- •Соответствующее выражение для кондуктивной составляющей теплового потока имеет вид:
- •При значениях параметра Пекле Pe 0.1
- •Другой пример: бассейн Уэд эль-Миа, где скорости осадконакопления не превосходили 145 м/млн.лет, а
- •Фильтрационное течение грунтовых вод вдоль почти горизонтальных осадочных слоев протяженностью сотни км практически
- •Термофизические параметры
- •Термофизические параметры осадочных пород
- •Табл. 2-1 Среднемировые значения петрофизических параметров основных осадочных фаций
- •Тепловые характеристики для смеси литологических единиц
- •Свойства поровых вод
- •Теплопроводность и теплогенерация осадочных пород заметно меняются с глубиной в зависимости от пористости,
- •Оценки теплопроводности пород через пористость, плотность и скорости сейсмических волн
- •Теплогенерация осадочных пород
- •Теплогенерация в осадочных породах
- •В моделировании полагают, что по мере погружения осадков генерация тепла в них увеличивается
- •Фундамент бассейна
- •Бассейны пассивных окраин характеризуются переходным типом фундамента: от континентального к океаническому
- •Фундамент в системе моделирования бассейнов
- •Теплопроводность мантии
- •Вклад радиационной составляющей в теплопроводность мантийных пород оказался заметно меньше предполагаемого в (Schatz
- •Теплоёмкость пород мантии
- •Скрытая теплота плавления
- •Выше - L - скрытая теплота плавления (90-100 кал/г для перидотитовых пород), Ts
- •1 и 3 – TL и Ts сухого гранатового перидотита (McKenzie, Bickle, 1988).
- •Основание литосферы и её термическая толщина определяется пересечением кривой солидуса пород мантии Ts(z)
- •Граничные условия в уравнении теплопроводности
- •Лапландский ледниковый горизонт (670 – 630 млн.л. назад). Обнаружен в Европе, Азии, Западной
- •Возмущения в температуре пород, вызванные вариациями поверх- ностной температуры, имеют тен- денцию затухать
- •Для морских осадков это температура дна моря. Она рассчитывается с учётом уменьшения температуры
- •Условия на нижней границе области счета
- •В основании области счёта поддерживается
- •Начальное распределение температур в литосфере
- •Для высоких значений Qo в интервале между температурами солидуса и ТМ (в области
- •Возраст окенической коры
- •Разностная схема решения уравнения теплопроводности
- •Шаги z, t
- •Корректность разностной схемы
- •Аналитические и полуаналитические решения для распределения температуры и теплового потока, используемые для проверки
- •Численное воспроизведение
- •2) Эрозия однородного полупространства
- •Сравнения с другими решениями
В моделировании полагают, что по мере погружения осадков генерация тепла в них увеличивается по закону:
A(z)=Am [1- (z)]
(см. слайд 10)
Вклад радиогенного тепла в тепловой поток
растет с толщиной осадочного покрова и может достигать 30% от теплового потока на поверхности в бассейнах с мощным осадочным заполнением (Баренцево море (Вержбицкий, 2002) и Мексиканский залив (McKenna et al., 1998)).
Материнские породы имеют высокое содержание U, Th, K
Cилурийские глины (Алжир; 70 м): A 6.7 мкВт/м3 ( Т 0.45°С) Баженовские глины (Уренгой; 15-50 м): A 7 мкВт/м3
( Т 0.5°С). Т.е. тепловой эффект выделения радиогенного тепла в узких слоях материнских пород ( 100-150 м) не имеет
заметного влияния на температуру и зрелость этих пород.
Фундамент бассейна
Осадочные бассейны могут развиваться на литосфере самого разного типа в соот- ветствии с глобальным цик- лом эволюции литосферы Земли.
Более того:
разные площади одного бас- сейна могут отвечать раз- ным типам литосферы (пас- сивные окраины, слайд 17). Тип фундамента бассейна может меняться также и в процессе его развития.
Бассейны пассивных окраин характеризуются переходным типом фундамента: от континентального к океаническому
Бассейн Сантос, Бразилия
Фундамент в системе моделирования бассейнов
При описании фундамента система моделирования использует 4 типа литосферы : 1) Континентальная; 2) Океаническая; 3) Задугового (краевого) моря и 4) с произвольными толщинами и характеристиками слоёв.
Вклад радиогенной составляющей пород литосферы в поверхностный тепловой поток увеличивается от Qрад = 10.5 мвт/м2 для океанической и Qрад = 13.-15. мвт/м2 для литосферы окраинных морей до Qрад = 23-25.мвт/м2 для континентальной литосферы.
Но различие в величине Qрад в различныx тектоническиx провинцияx может быть значительнымв
|
|
|
|
|
тип коры |
гранитные |
метаморфизи- |
базальтовая |
докембрийская |
|
плутоны |
рованная кора |
кора |
кора |
Q (мВт/м2) |
50-150 |
20-50 |
10-30 |
10-30 |
rad |
|
|
|
|
|
КОНТИНЕНТАЛЬНАЯ ЛИТОСФЕРА |
|
|||
интервал |
теплопроводность |
плотность |
генерация тепла |
||
глубин (км) |
(Вт/м°К) |
(г/см3) |
(мкВт/м3) |
||
|
кора |
“гранитный слой” |
|
||
0 Z 15 |
2.72 |
|
2.75 |
1.465 (0 Z 5) |
|
|
|
|
|
|
|
|
кора |
“базальтовый слой” |
0.840 (5 Z 15) |
||
|
|
||||
15 Z 35 |
1.88 |
мантия |
2.90 |
0.21 |
|
|
|
|
|
||
Z 35 |
3.56 или F(T)* |
3.30 |
0.0042 |
||
|
ОКЕАНИЧЕСКАЯ ЛИТОСФЕРА |
|
|||
0 Z 6.5 |
1.88 |
Кора |
2.85 |
0.630 |
|
Мантия |
|||||
|
|
|
|
||
Z 6.5 |
3.56 или F(T)* |
3.30 |
0.0042 |
||
|
ЛИТОСФЕРА КРАЕВЫХ МОРЕЙ |
|
|||
0 Z 16 |
2.51 |
Кора |
2.85 |
0.630 |
|
|
|
|
|
|
|
|
|
Мантия |
|
|
|
Z 16 |
3.56 или F(T)* |
3.30 |
0.0042 |
Вклад радиогенной составляющей пород литосферы в поверхностный тепловой поток увеличивается от Qрад = 10.5 мвт/м2 для океанической и Qрад = 13.-15. мвт/м2 для литосферы окраинных морей до Qрад = 23- 25.мвт/м2 для континентальной литосферы.
Такие значения генерации тепла Qрад означают, что, например, для потока тепла на поверхности 40 - 50 мВт/м2 температура пород океанической мантии будет расти с глубиной в полтора - два раза быстрее, чем в континентальной мантии.
Теплопроводность мантии
Теплопроводность мантии изменяется с глубиной как функция температуры ((Hofmeister, 1999; McKenzie et al.,2005) :
С ростом температуры увели- чивается плотность фононов, но падает длина их свободного пробега, что ведёт к падению теплопроводности. Но одновре- менно с ростом Т растёт интен- сивность радиационного пере- носа тепла фотонами.
kH (T ) |
|
b |
03 dm (T 273)m |
||
1 c T |
|||||
|
|
|
|||
Здесь Т в °С, b = 5.3, c = 0.0015, d |
= 1.753 10-2, d = - 1.0365 10- |
||||
4, |
|
|
0 |
1 |
|
|
|
|
|
d2 = 2.245 10-7, d3 = - 3.4071 10-11. Вычисленные k(T) показаны
Вклад радиационной составляющей в теплопроводность мантийных пород оказался заметно меньше предполагаемого в (Schatz and Simmons, 1972) из теоретических соображений.
Соответственно, теплопроводность мантийных пород, вычисленная по алгоритму (Schatz and Simmons, 1972), заметно превышает измеренные значения при температурах 900 – 1000ºС и выше.
Теплоёмкость пород мантии
Вариации теплоёмкости пород мантии (диопсид – смесь 0.89 форстерита и 0.11 фаялита) с температурой, Ср (Т), определяются согласно (Berman, Aranovich, 1996) и
(McKenzie et al.,2005) по формуле:
Ср = 1618.27 – 12538.9 / SQRT(T°K) – 173.7 107 /
T(°K)3 (Дж/кг°K)
Они согласуются с литературными данными для диопсида (Cawthorn, 1975; Robie et al,1978; Bottinga, Steinmetz, 1979). Изменение Ср от Т согласно этой формуле:
Т (°К) |
= 298.15 |
600 |
1000 |
1600 |
Т (°С) |
= 25 |
226.85 |
726.85 |
1326.85 |
Ср (Дж/кг°K) |
= 826.55 |
1098.3 |
1220.0 |
1304.4 |
Скрытая теплота плавления
Термическая эволюция бассейна на рифтовом этапе его развития или в периоды его тепловой реактивизации характеризуются заметными изменениями толщины литосферы в ответ на плавление пород в её основании или на затвердевание пород прилегающей астеносферы.
Выделение или поглощение скрытой теплоты плавления, заметно увеличивают «тепловую инерционность литосферы бассейна», замедляя процессы её остывания или разогревания.
Энтальпийное приближение: (Карслоу, Егер, 1964).
Доля расплавленной фракции, f, в породе линейно увеличивается с ростом температуры, T, в пределах интервала температур между солидусом и ликвидусом породы:
f |
|
T Ts |
при Ts < T < Tl |
|
Tl Ts |
||||
|
В этом интервале температур в уравнении теплопроводности в члене c |
||
( / t) теплоёмкость Ср заменялась на Cp‘ : |
L |
|
Cp' Cp |
||
Tl Ts |