Строение и история развития литосферы

На этом шаге было произведено изменение параметров модели. Если на первых двух временных шагах разрез включал всю литосферу, то для возраста 200 млн лет нижняя граница разреза была ограничена глубиной 70 км, установленной по данным ГСЗ. На этой глубине задавался тепловой поток, определенный на предыдущем этапе расчета. Он изменялся вдоль профиля от 34 до 40 мВт/м².

На разрезе коры после изменения геометрии модели геологическое строение отражено более детально. В модели учитывалась блоковая структура разреза, выраженная в наличии субвертикальных изгибов пластов на контакте пород гранитно-метаморфического комплекса с метаморфизованными осадочными отложениями палеозоя и базальтов с гипербазитами в юго-западной части профиля (интервал 100–200 км). Эти блоки контактируют по разломным или флексурным зонам, возникшим, судя по возрасту осадочных пород, не позже среднего палеозоя. Такой возраст дислокаций позволяет предположить отсутствие в них современных проницаемых зон, по которым осуществляется тепломассоперенос, поэтому разломы и флексуры моделировались только как границы контрастных теплофизических сред.

В результате расчета данной теплофизической модели с шагом дискретизации 80 млн лет было получено распределение температур в земной коре Баренцевоморского региона для времени 120 млн лет (граница раннего и позднего мела) (рис. 3). Пониженная теплопроводность гранитов по сравнению с базальтами и ультрабазитами и сложная геометрия слоев разреза обусловливают рефракцию глубинного теплового потока. Его плотность уменьшается на краях тел относительно пониженной теплопроводности из-за «обтекания теплового потока», но увеличивается в пространстве между ними, т. к. энергия идет «по пути наименьшего сопротивления». Это явление рефракции хорошо известно в теории потенциала и не требует дальнейших пояснений. Здесь уместно заметить, что практически все искажения глубинного теплового потока, отмеченные по результатам моделирования, связаны со структурно-теплофизическими неоднородностями разреза. Влияние рельефа морского дна как фактора, вызывающего искривление изотерм и искажение глубинного теплового потока, для данного профиля пренебрежимо мало по сравнению с вышеописанным эффектом.

Рис. 3. Температурный разрез (°С, сплошные линии) и распределение теплопроводности (Вт/(м·К), пунктирные линии) вдоль профиля ГСЗ-76 для времени 120 млн лет назад.

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

Следующий временной шаг рассчитан для возраста 60 млн лет (эоцен). В это время не происходило каких-либо структурных перестроек, кроме продолжающегося осадконакопления терригенных толщ. Но скорость осадконакопления для этого времени настолько низка, что сам процесс, как показали расчеты, не искажает глубинный тепловой поток. Следовательно, тепловое поле для этого момента отражало только релаксационные явления после воздействия позднепалеозойского термического источника. Распределение температур и тепловых потоков для этого времени близко к тому, которое было получено для 120 млн лет. Это говорит о том, что геотермическое поле становится квазистационарным, т. е. динамика его изменения во времени прекращается.

Положение изотерм катагенетического температурного интервала (110°С и 160°С) осталось на той же глубине, что и для предыдущего временного шага, т. е. от 4 до 5,5 км и от 4,5 до 6,5 км соответственно. Таким образом, приблизительно в конце позднего мезозоя в данном регионе установилось стационарное тепловое поле. На это же время указывают и геологические данные (Верба, Шаров, 1998), свидетельствующие о том, что в юре и позднее, с уходом рифтогенной активности за пределы Баренцево-Карской плиты, на всей ее площади устанавливается пострифтовый этап эволюции.

Современная термическая структура получена для следующего временного шага (0 млн лет). Профиль рассчитанного теплового потока на поверхности повторяет изгиб изотерм (см. рис. 2). Он достигает максимальных значений 65–67 мВт/м²

в интервалах L=500–700 км и приурочен к выступам гранитно-метаморфических пород, кровля которых по сейсмическим данным расположена на глубине 4–6 км под поверхностью дна. Полученный по результатам моделирования фоновый тепловой поток составляет 54 мВт/м². Это значение оценивается путем осреднения рассчитанных по профилю тепловых потоков (Хуторской, 1996).