Бесплатный автореферат и диссертация по геологии на тему
Математическое моделирование термической эволюции осадочных бассейнов и условий реализации их углеводородного потенциала
ВАК РФ 04.00.12, Геофизические методы поисков и разведки месторождений полезных ископаемых

Автореферат диссертации по теме "Математическое моделирование термической эволюции осадочных бассейнов и условий реализации их углеводородного потенциала"



^^ На правах рукописи

Л V %v

ГАЛУШКИН ЮРИЙ ИВАНОВИЧ

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕРМИЧЕСКОЙ ЭВОЛЮЦИИ ОСАДОЧНЫХ БАССЕЙНОВ И УСЛОВИЙ РЕАЛИЗАЦИИ ИХ УГЛЕВОДОРОДНОГО ПОТЕНЦИАЛА

Специальность 04.00.12 - Геофизические методы поисков и разведки

месторождений полезных ископаемых 04.00.13 - Геохимические методы поисков месторождений полезных ископаемых

АВТОРЕФЕРАТ

диссертации на соискание ученой степени доктора технических наук

Москва 1998

Работа выполнена в Учебно-научном Музее землеведения Московского Государственного Университета им. М.В.Ломоносова и в лаборатории геологии и геохимии нефти и газа Всероссийского научно-исследовательского института геологических, геофизических и геохимических систем

Официальные оппоненты: доктор физико-математических наук.

заведующий лабораторией ВНИИГеосистем, профессор Полячснко Анатолий Львович

доктор геолого-минералогических наук, заведующий лабораторией ВНИГНИ профессор Шеин Виктор Степанович

доктор физико-математических наук, ■заведующий лабораторией института Океанологии РАН им. Ширшова профессор Сорохтин Олег Георгиевич

Ведущая организация: Институт проблем нефти и газа РАН и Министерства общего и профессионального образования РФ, Москва

Защита состоится " 1998г. в) 7 час

на заседании Диссертационного Совета ДО 71.10.01 при Всероссийском научно-исследовательском институте геологических, геофизических и геохимических систем по адресу: 113105, Москва, Варшавское шоссе, д. 8

С диссертацией можно ознакомиться в научной библиотеке ВНИИгеосистем

Автореферат разослан ____"__1998 г.

Ученый секретарь диссертационного совета - профессор

В.С. Лебедев

ВВЕДЕНИЕ

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

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

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

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

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

1. Расчет эволюции температурного режима и полей деформаций литосферы в процессе континентального рифтогенеза.

2. Численный анализ изменения температурного режима континентальной и океанической литосферы в бассейнах с разнообразной тектонической историей и нестационарным горизонтальным полем скоростей перемещения вещества литосферы.

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

4. Оценка изменения степени метаморфизма органического вещества, реализации нефтегазогенерационного потенциала бассейна и порога первичной эмиграции жидких УВ.

5. Восстановление спектра химико-кинетических реакций, управляющих процессами созревания ОВ осадков.

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

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

Работа выполнена в секторе геодинамики Музея землеведения Московского Государственного университета им. М.В. Ломоносова и в лаборатории геологии и геохимии нефти и газа во ВНИИГеосистем и суммирует основные научные результаты, полученные автором за последние 20 лет.

Личный научный вклад автора состоит в решении следующих вопросов:

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

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

- выявление на основе численных расчетов закономерностей и особенностей развития термического режима бассейнов на рифтозом и пострифтовом этапах их развития.

- численный анализ основных закономерностей и особенностей в истории созревания ОВ и реализации УВ потенциала осадочных бассейнов разного геодинамического типа

Автор является идеологом и основным исполнителем разработок. Теоретические и методические результаты получены непосредственно автором.

Основные защищаемые положения.

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

2) Численный анализ термического режима, рельефа поверхности и гравитационного поля осадочных бассейнов континентальных окраин и задугового спрединга предполагает учет таких процессов как латеральный теплообмен между разновозрастными или разнотиповыми (океаническим и континентальным) блоками литосферы, временные вариации в скоростях спрединга и перескоки осей спрединга. Первый из перечисленных процессов определяет амплитуду рельефа и размеры палеоспрединговых структур в

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

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

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

5) Сравнение измеренных значений отражательной способности витринита, 110%, со значениями, вычисленными по кинетической модели созревания витринита, сопоставление наблюдаемого и вычисленного глубинных профилей температур и анализ вариации амплитуды тектонического по!ружения фундамента являются действенными контролирующими факторами в системе численного моделировании бассейнов.

6) Численная модель упруго-пластичного растяжения литосферы с конечным пределом пластичности адекватно воспроизводит процесс утонения литосферы бассейнов при рифтогенезе с соблюдением временных и пространственных масштабов явления. Это моделирование вместе с численным анализом многих кривых тектонического погружения бассейнов подтверждают тот факт, что продолжительности периодов растяжения и тепловой активизации литосферы бассейнов могут составлять десятки млн. лет, существенно превосходя времена, предполагаемые в моделях мгновенного растяжения литосферы.

7) Разработанная модификация пакета программ по моделированию бассейнов позволяет анализировать изменения термического состояния осадочных пород (и пород фундамента), а также процесс формирования ореола зрелости вокруг интрузивных тел, внедрившихся в реальный осадочный разрез бассейна. Форма и размеры ореолов зрелости ОВ (профилей Яо%) в окрестности интрузивных тел согласуются с конечными временами их формирования (от минут до полугода) иногда и без прямого контакта расплавленных интрузивных и вмещающих пород. Резкие изменения глубинных градиентов степени катагенеза ОВ пород, наблюдаемые в глубоких горизонтах рифтогенных бассейнов, можно объяснить высоким тепловым потоком от гидротермальной и интрузивной деятельности на рифтовом этапе развития бассейнов.

8) Применение алгоритма с различными частотными факторами А1 в программном блоке восстановления кинетического спектра по данным экспериментального пиролиза образцов материнских пород позволяет существенно улучшить подбор кинетических параметров реакций по сравнению с традиционным методом, предполагающим одно значение А1 для всех кинетических реакций. Модификация алгоритма рассматриваемого программного блока путем включения в него члена, ответственного за геологический этап созревания исследуемого образца материнской породы, позволяет получать спектры химических реакций, управляющих созреванием керогена, заметно

измененые в области низких энерий и осуществлять восстановление спектра для образцов с исходной степенью зрелости ко < 0.7-0.8% .

Научная новизна диссертации состоит в разработке автором следующих проблем:

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

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

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

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

5. Используя модификацию пакета программ по моделированию бассейнов, автор первым численно показал, что форма и размеры ореолов зрелости ОВ (профилей 1*о%) в окрестности интрузивных тел согласуются с конечными временами формирования этих тел (от минут до полугода) иногда и без прямого контакта расплавленных интрузивных и вмещающих пород.

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

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

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

С 1975 по 1992 г. автор в качестве ведущего исполнителя вел совместные работы по темам с ПГО "Ссвморгеология", НПО "Союзморгео" и НПО "Южморгео" по изучению перспектив нефтегазоностости бассейнов переходных зон Арктики, Антарктики и Северной Атлантики, а также по термическим проблемам формирования различных структур океанического дна в связи с природой образования сульфидных месторождений. С 1981 по 1992 год автор в сотрудничестве с Я.Б. Смирновым изучал термическую эволюцию и условия созревания ОВ ряда бассейнов Восточной Европы. С 1988 года в сотрудничестве с Н.В.Лопатиным, и с 1993 по 1996 г. под его руководством как сотрудник лаборатории нефтяной геологии и геохимии ВНИИГеосистем, автор занимался проблемами геотермии инефтегазоносности Западно-Сибирского, ЮжноБаренцева, Печерского бассейнов, а также бразильского бассейна ГТарнаиба. Все рекомендации практического, методического и теорегического характера были своевременно переданы в ПГО"Ссвморгеология", в НПО "Союзморгео" и геологическим организациям "Архангельскгеология", "Уренгойгеология" и реализованы в этих организациях.

Аппробация диссертации. Основные разделы и отдельные положения работы докладывались автором на международных симпозиумах, совещаниях, конгрессах, посвященных геодинамике океанической и континентальной литосферы, термической эволюции Земли, проблемам моделирования бассейнов: Иркутск, 1975; Южно-Сахалинск,1976; Геленджик, 1990,1992; Звенигород, 1991; Аксаково,1993; Москва,1993; Фрайберг (Германия), 1987; Ставангер (Норвегия), 1991; Абердин (Англия), 1992; Париж,1992; Вена,1994; Хьюстон,1995; Абу-Даби (Кувейт), 1995), на всесоюзных и всероссийских совещаниях (Москва, 1976,1977,1982,1988,1991; Ялта, 1982; 1988;Гелснджик, 1984,1986,1988; Иркутск,1989; Владивосток,1976,1987,1989; Тбилиси,1980) и многократно - на семинарах и чтениях в МГУ.

Публикации. По теме диссертации автором опубликовано более 120 работ, среди которых 72 статьи в журналах и сборниках и 6 монографий (в соавторстве). Диссертация в целом обсуждалась на Ученом Совете Музея землеведения МГУ и Ученом Совете ВНИИГеосистем.

Благодарности. Автор приносит глубокую благодарность директору Музея землеведения МГУ, академику РАЕН, профессору С А Ушакову за постоянное внимание и поддержку работы и за большое участие в формировании научного мировоззрения соискателя в области геодинамическон природы формирования бассейнов.

Автор искренне благодарен д.г.-м.н. Е.П.Дубинину, совместно с которым решались и обсуждались многие из проблем, касавшихся океанической литосферы.

Чувство глубокой благодарности автор сохраняет к старшему научному сотруднику лаборатории геотермических исследований ГИНа АН РАН Я.Б. Смирнову, совместная работа с которым оказала большое влияние на формирование мировоззрения автора как геотермика.

Автор благодарен доктору наук, заведующему лаборатории нефтяной геологии и геохимии ВНИИГеосистем Н.В.Лопатину за предложение ряда научных тем, связанных с моделированием бассейнов ("Химико-кинетическое

моделирование процессов генерации углеводородов", "Моделирование процесса созревания OB в окрестности интрузивных тел" и др.), и плодотворное обсуждение научных результатов.

Автор благодарен за ценные замечания в области химико-кинетического моделирования процессов преобразования керогенов и моделирования бассейнов, высказанные в ходе научных дискуссий в 1994-1995 годах докторами Р. Шайфером, Р. Литтке, и X. Пульхау в институте Нефти и Органической геохимии в Юлихе и во время их пребывания в Москве, а также докторами Д. Бэнсом и Д. Кюри во время научного семинара рабочей группы по низкозрелым газам в в Исследовательском Центре компании "Эксон Продакшн" в Хьюстоне в декабре 1995 года.

Для автора оказались очень полезными совместная работа и научные дискуссии с докторами А. Бсртольдом, Мунзером Махусом, Г.Е.Яковлевым, A.B. Муравьевым, О.И. Симоненковой, Р.И., Кутасом, A.B. Ершовым, А.И. Шемендой, P.M. Деменицкой, А.Д. Дучковым , Л.С. Соколовой, В.М. Сугробовым, A.B. Ершовым, Е.А. Тищепко, Л.И. Лобковским, A.A. Ковалевым, A.M. Городницким, Б.Г. Поляком, М.Д. Хуторским и другими.

Автор благодарен A.B. Свешникову за помощь в оформлении диссертации. Завершающая стадия работы выполнялась при финансовой поддержке Международного (ISF grant MER 300, 1994-1995 гг.) и Российского (РФФИ 19961998) Фондов Фундаментальных Исследований.

СОДЕРЖАНИЕ РАБОТЫ

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

Глава 1. ГЕОДИНЛМИЧЕСКИЕ ПРЕДПОСЫЛКИ ФОРМИРОВАНИЯ РИФТОГЕННЫХ ОСАДОЧНЫХ БАССЕЙНОВ И СПЕЦИФИКА МОДЕЛИРОВАНИЯ ИХ ТЕРМО-МЕХАНИЧЕСКОЙ ИСТОРИИ 1.1 Эволюционный ряд рифтогенных бассейнов и разнообразие геодинамических обстановок их формирования

В разделе на основе имеющихся в литературе классификаций осадочных бассейнов (Высоцкий, Кучерук, 1978; Кучерук, Ушаков, 1985; Артюшков, 1993; и др.) показано, что рифтогенные бассейны составляют существенную часть нефтегазоносных бассейнов Земли и их эволюция - возникновение, развитие и преобразование или разрушение - составляет часть общей глобальной эволюции литосферы Земли. Большое разнообразие бассейнов, строения их осадочной толщи и фундамента объясняют на основе так называемого эволюционного ряда бассейнов, приуроченного к тектоническому циклу развития литосферы. Этот цикл начинается расколом континента, включает образование океанической впадины и заканчивается орогенезом в процессе поглощения океанической коры, закрытия океана, сближения и столкновения континентов (Ушаков, Галушкин, 1983; Кучерук, Ушаков, 1985; Sokolov, 1985; Khain, Sokolov and Shein, 1991; Хаин, Болуховский, 1993; Kleshov and Shein, 1991; Cloetingh et a!., 19%).

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

В процессе развития бассейна тектоническое погружение его фундамента можег обуславливаться как поверхностной нагрузкой осадков и воды, так и термическими и тектоническими изменениями в распределении плотности пород фундамента с глубиной, в том числе и вызыванным фазовыми переходами пород коры и мантии литосферы; MacKenzie 1978; Hegarty et al.1988; Rehault et al. 1990; Артюшков, 1993; Cloetingh et al.,1994,1996) . Однако, как следует из обсуждения вопроса в диссертации, вклад фазовых переходов пород нижних горизонтов континентальной коры в более плотные гранулитовые и эклогитовые фации в тектоническое погружение внутриконтинентальных бассейнов может быть оценен лишь качественно в силу слабой изученности кинетических параметров данного перехода (Соболев, 1978; Артюшков и др., 1979, 1982, 1998; Артюшков, Бэер, 1983; 1984; Hamdani et al.,. 1991;1994;Barid et al., 1995; Артюшков, Мернер, 1997).

1.2 Термо-механические модели образования рифтовых бассейнов

Утонение литосферы бассейна и астеносферный диапиризм под осевой зоной континентального рифта рассматриваются как основные процессы,

управляющие формированием рифтогенных бассейнов (Артюшков, 1993; Serahne and Mallavielle, 1994; Vejbeck, 1994). При этом анализ геологической и геофизической информации по формированию структур растяжения, приводит к выводу, что процесс рифтинга может занимать ог 16 до 70 млн. лет для континентальной и от 3 до 12 млн. лет для литосферы задуговых бассейнов (Takeshita and Yamaji, 1990; Anders and Sleep, 1992; Pedersen, 1994) и, следовательно, его нельзя считать мгновенным, как это предполагается, например, в модели (McKenzie, 1978;1980). Экспериментальное физическое моделирование процесса рифтогенеза континентальной литосферы, осуществленное в лаборатории моделирования геодинамических процессов МГУ А.И. Шемендой (Шеменда. 1985; Shemenda, 1993), продемонстрировало ведущую роль утонения литосферы снизу в процессе рифтогенеза при существенно меньших амплитудах деформаций приповерхностных слоев литосферы. Эти эксперименты показали также, что эволюция поля деформаций при рифтогенном растяжении литосферы хорошо описывается в рамках модели деформаций упруго-пластичной литосферы с конечным пределом ползучести, x¡ :

czjèl = 0 при а,л < х, и cejèt = [(ахх - ij /а)" при ахх > х, (1) где а имеет смысл эффективной вязкости среды: ахх = х, + а- (cejcl)"" В численной модели рифтогенеза при расчете деформаций мы использовали реологический закон (1) вместе с обычным упрощающим предположением о независимости скоростей деформаций éejet от глубины z (Alvares et al., 1984; Pederson, 1994). Тогда выполнение двух условий - постоянство силы, растягивающей литосферу, в каждом ее сечении (ахх-Н(x.t) =F(l)) и - выход

скоростей смещения пород литосферы у границы зоны пластичности x=XTS к заданному значению скорости растяжения Vo:

XTS XTS(t) u-TO/.v

Vo= J (Л,/*).«& = (%• J (2)

о a о "Vх' ч

дало возможность путем итераций найти границу области пластических деформаций XTS как функцию времени, если заданы реологические параметры х, а и п и если определена форма нижней границы литосферы H(x,t) на данный момент времени I. После этого поле скоростей смещения частиц литосферы определялось но формулам:

4,(x.t) = V0;V2(x,z,t) = 0 при х>XTS и

X

V„(x,t) = Jf(x,t)-dx;VI{x,z,t)= -g(x,t) z при 0 <х.< JO'S (3)

о

Найденное поле скоростей использовалось для определения нового положения нижней границы литосферы через шаг времени At и повторения цикла вычислений. Эта модель в соответствии с наблюдаемым в природе и в экспериментах обеспечивала плавное уменьшение скоростей деформаций в растягиваемой литосфере ¿ejcl с приближением к границе пластичности и объясняла ведущую роль утонения литосферы снизу. Полученное поле скоростей Vx(x,t) и Vz(x,z,l) использовалось в решении нестационарного, двумерного уравнения теплопроводности с конвективными членами, (см. раздел 4.1), которое определяло эволюцию температурного поля литосферы ь процессе рифтогенеза.

В диссертации приведен пример численного моделирования полей деформаций и температур литосферы в процесс«; растяжения последней с постоянной скоростью Уо=0.5 см/год и с разумными значениями реологических

параметров пород литосферы: a=4-I(f2 пуаз, i=l(f дин/см2, п-1 (Галушкин и др., 1988). Расчеты показали, что в результате такого растяжения подошва литосферы у оси рифта в течении 20 млн. лет поднялась с глубины около 145 км до 37 км. В процессе растяжения изотермы приближались к поверхности в области "шейки" и тепловой поток на оси "шейки" непрерывно рос от начального значения около 30 мВт/м2 до 37 мВт/м2 через 5 млн.лет рифтинга, затем до 57 мВт/м2 через 12 млн. лет и достигал 112 мВт/м2 после 20 млн. лет растяжения. Характерно, что рассчитанные значения теплового потока заметно отставали от своих стационарных величин, особенно в крыльях структуры (Галушкин и др., 1988). Таким образом, численное моделирование в рамках модели растяжения упруго-пластичной литосферы с конечным пределом пластичности неплохо передает временные и пространственные масштабы процесса утонения литосферы при рифтогенезе.

1.3 Кинематика литосферных плит и специфические типы осадочных бассейнов

Раздел 1.3 диссертации посвящен применениям разработанной автором компьютерной программы для расчета параметров модели мгновенной кинематики движения литосферных плит Земли (Галушкин, Ушаков, 1978, 1979; Ушаков, Галушкин, 1978), основанной на алгоритмах американских геофизиков (Minster and Jordan, 1974), Эта модель позволяет судить о природе формирования зон растяжения, сопутствующих образованию рифтогенных бассейнов, о современном напряженном состоянии их литосферы. Она дает возможность представить скорости формирования новой литосферы в центрах спрединга, поглощения плит на конвергентных границах, включая глобальные пояса сжатия. Сложный характер взаимодействия краев литосферных плит и блоков при их взаимных перемещениях, а также особенности в простирании их границ ведут к формированию разнообразных типов бассейнов, и обуславливают существование структур, в которых области растяжения будут соседствовать с областями сжатия, а районы высокого термического градиента с низкотемпературными областями. Изменения в относительных движениях литосферных плит и блоков могут служить причиной как зарождения новых, так и отмирания активных до этого времени рифтовых структур с последующим образованием на их месте авлакогенов или бассейнов палеоспрединговых хребтов. Эти изменения могут привести также к появлению осложнений на отдельных участках дивергентных, конвергентных и трансформных границ плит и формировать специфические для данных обстановок бассейны с локальным растяжением фундамента.

Глава 2. "МОДЕЛИРОВАНИЕ БАССЕЙНОВ" КАК СИСТЕМА ЧИСЛЕННОГО

АНАЛИЗА ТЕРМИЧЕСКОЙ ЭВОЛЮЦИИ ОСАДОЧНЫХ БАССЕЙНОВ И РЕАЛИЗАЦИИ ИХ УГЛЕВОДОРОДНОГО ПОТЕНЦИАЛА 2.1 Основные компоненты системы моделирования ГАЛО

Большая часть диссертации, включающая главы 2 и 3, посвящена принципам построения, исходному алгоритмическому базису и применениям разработанного автором программного комплекса по моделированию бассейнов, позволяющего численно реконструировать историю погружения и эволюцию термического режима осадочных пород и реализацию нефтегазогенерационного потенциала материнских свит в процессе эволюции бассейна. Общий принцип построения комплекса показан на рис.1 и обсуждается в диссертации. Одной из основных целей автора при создании пакета по моделированию бассейнов (ГАЛО), была по возможности полная интеграция процессов, определяющих

ВВОД ДАННЫХ ПО СТРОЕНИЮ И ЭВОЛЮЦИИ БАССЕЙНА

ИСХОДНЫЕ ПАРАМЕТРЫ МОДЕЛИРОВАНИЯ

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

ЧИСЛЕННОЕ МС ДЕЛИРОВАНИЕ

расчет пористости пород

реконструкция истории погружения бассейна

расчет теплопроводности, теплоемкости итеплогеие-рации пород

—► профиль

РЕШЕНИЕ УРАВНЕНИЯ ТЕПЛОПЕРЕНОСА: температуры

температура пород фундамента и осадков

вычисление

значений

Рис.1 Принцип построения и основные блоки системы моделирования бассейнов ГАЛО

термический режим в осадочной толще бассейна, с термической и тектонической историей подстилающей литосферы. Автор намеренно ограничился анализом "плоских" (одномерных) бассейнов и изменил подход к решению проблемы, предпринятый в вариантах известных пакетов по моделированию бассейнов MATOIL, GENEX, TEMISPACK, PDI (Welte and Yukler, 1981; Tissot et al.,1987; Nakayama and Lerche, 1987; Welte and Yalcin, 1988; Espitalie et al.,1988; Ungerer, 1990; Ungerer et al.,1990). в плане постановки граничных условий. Он включил в область численного поиска распределения температур, помимо осадочной толщи, подстилающую литосферу и часть астеносферы. Такой подход наряду с корректной трактовкой теплообмена между осадочными слоями и фундаментом давал возможность путем численного анализа теплового режима, распределения плотностей пород и вариаций амплитуд тектонического погружения фундамента проследить взаимосвязь процессов формирования осадочной толщи бассейнов и эволюции ее температурного режима с глубинными тепловыми и тектоническими процессами d астеносфере и подстилающей литосфере рифтогенных бассейнов и привлечь анализ вариаций амплитуд тектонического погружения фундамента для контроля исходных параметров моделирования бассейнов. Учет термической истории погружающегося образца нефтематеринских пород в процедуре восстановления параметров эффективных химико-кинетических реакций, управляющих процессами созревания OB в осадках, также составлял момент, отличающий пакет ГАЛО от разработанных в германском, французском и американском исследовательских центрах.

2.2 Реконструкция истории погружения осадочных слоев

При погружении осадки сжимаются с высвобождением части поровых вод. Процесс консолидации рассматривается в обычном предположении (Perrier and Qublier, 1974): изменение обьема твердой матрицы пород при их захоронении пренебрежимо по сравнению с уменьшением их порового пространства

[l-P(Zj] dZ, = fl-P(Z2)J dZ2 (4) где P(Z) - пористость на глубине Z и dZj и dZ2 - толщины слоя на глубинах Z/ and 7.2, соответственно. Известная процедура разуплотнения (backstripping) для элементарного слоя осадков dZ основана на соотношении (4) и законе уменьшения пористости с глубиной (Sclater and Christie, 1980; Deming and Chapman, 1989 и др.), из которых в пакете ГАЛО предусмотрено использование экспоненциального, линейного и задаваемого произвольной функцией глубины в зависимости от литологии и физического состояния изучаемых пород. Значения параметров Р(0)- средней пористости пород в верхних 100-200 м разреза, и В -масштаба изменения пористости с

глубиной, используемых в законах изменения пористости, определяются либо из анализа петрофизических свойств пород конкретного разреза, либо на основании среднемировых данных для основных литологических единиц и обсуждаются в соответствующем разделе диссертации. При этом пористость пород, представленных комбинацией литологических единиц, вычислялась, согласно (Doligezet al., 1986):

----= V ——— (5)

где С; - доля i-ой литологической единицы в породе, п - число разных фракций в породе и Pj(Z) - пористость пород i-ой литологической единицы на глубине Z.

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

отложенных толщ (2-3 раза) и помимо нагрузки лежащих выше осадочных пород существенно определяется литологией пород, их минеральным составом, а также историей погружения свиты (см., например, рис. 26). В наших моделях использовалось стандартное предположение о том, что разуплотнения пород при снятии нагрузки во время эрозии не происходит. Вариации пористости при растворении минералов, отложении цемента или процессов перекристаллизации учитывались в модели косвенно через эффективное изменение масштабного параметра В со временем погружения породы в бассейне.

Эволюция осадочного бассейна может сопровождаться появлением зон с аномально высоким норовым давлением (АВПД), и часть раздела 2.2 диссертации посвящена обсуждению механизмов его генерации. Анализ современной теоретической и геологической литературы но этому вопросу показал, что попреки традиционным представлениям. акватермальное давление и неравновесное уплотнение не являются распростракнеными причинами появления АВПД за исключением районов с лавинным осадконакоилением. Тектоническое сжатие в бассейнах современных тектонически активных областей, преобразование керогена и вторичный крекинг УВ в материнских свитах бассейнон, а также движение метеорных вод под действием гидравлического напора в неглубоких осадочных горизонтах с обильным обводнением могут рассматриваться как наиболее вероятные источники АВПД в современных бассейнах.

2.3.Восстановление термической истории бассейна

Распределение температуры в осадочной толще и подстилающей литосфере бассейна находится решением уравнения теплоиереноса, записанного в консерватвной форме (Паскоиов и др.,1988):

д(С\ Т)/а + 0{Cnv VT)/8Z = (?[K((7I7(7Z)]/(7i + Л (7., I) (6) где T- температура, К - теплопроводность, Cv=pCp - обьемная теплоемкость и Ср - 1еплоемкость на единицу массы, р - плотность, Л - генерация тепла в единице обьема, V - скорости г рунтовых вод и воды, выжатой при консолидации осадков , относительно "неподвижного" фундамента, Z - глубина и t - время. Уравнение решалось в системе координат, связанной с фундаментом. Воспроизводился природный процесс осадконакопления путем наращивания сверху каждый элементарный промежуток времени ill области счета на элемент dZ, который представлял количество осадков, отложенных на поверхности за интервал til при температуре, отвечающей климатическим условиям в момент отложения (или уменьшения области на этот элемент в периоды эрозии). Учитывалась консолидация осадков под действием веса вновь отложенного элемента и соответствующие деформационные смещения узловой сетки (Галушкин, 1990; Mukhous, Galushkin. and Lopatin, 1997). Конвективный член в уравнении (6) мог быть связан только с движением грунтовых вод в погружающемся бассейне, вызванным выжиманием норовых вод при консолидации осадков, градиентом гидростатического давления от перепада рельефа, вариациями солености и температуры грунтовых вод.

Одной из особеностей уравнения (6) являлось то, что теплофизические параметры осадочных пород, будучи функциями их литологическото состава, пористости и температуры, могли заметно изменяться с глубиной и временем. Плотность (ps), обьемная теплоемкость(Cvs), и теплопроводноеть(Кя) пород определялась через их значения для минерального скелета (pm. Cvm. Km) и воды (rw. О»'. Kw) Doligez et al.,1986; Doming and Chapman, 1989; Ungerer et al„ 1990) :

и 0-|

01 -

X 1 -

4J

Г 2-J-

U,

тектоническое погружение

фундамента

FACI

-----|__ТЛЮ

20-

40

ПО

500°С _

íoxo ^ „ - —

750°С

фазовый переход

1000°С

■^основание

л

'J--4-.

I 1 ч

I ' V

I 1 \ ■I I

еоо

ж)

о

100-

а)

Л >>

ей 50 О. 0> с Я о н

О 8 0.6 -0.4 О 2

-т-г-г-Т ' ' i I -500 -300 -100 Бремя (млн.лет)

1-00

О

-500 -300 -100

Бремя (млн.лет)

-400 -200

Время (млн.лет)

Рис. Реконструкция температурно-нременной истории бассейна Уэд эль Mita площади Такхухт (см. текст)

0

s.

Ps(Z) = Pm 1¡ - P(Z) 1 +Pw p(Z) Cvs(Z) = Cvm [ 1 - P(Z) ] + Cm P(Z) (7)

Ks(Z)=Km(i-P(Z)) KnP(Z) Нелинейная зависимость теплопроводности от пористости подтверждается большим числом экпериментальных исследований (Beck 1976; Oxburgh, AndrewsSpeed 1981 и др.). Для воды: Cpw = 4186.8 Дж/кг °С, pw = 1030 кг/м3 и Л'н определяется в мВт/м°С как (Deming and Chapman, 1989):

Kw= 0.565 + 0.00188 Т-0.00000723 Т2 ( 0<Т<137°С) Kw = 0.602 +0.00131 Т-0.00000514 Т2 ( 137 < Т< 300°С) (8) Теплопроводности минерального скелета и воды также менялись с температурой (Burrus,Andebert,1990; Cermak,1995; Kukkonen and Joeleht.1996):

Km = Km(T=0°C)/( 1+A1-T) (9) Наиболее сильные изменения теплопроводности с температурой, как следует из данных, приведенных в диссертации, характерны для пород типа галита, для которых теплопроводность уменьшается в полтора раза при возрастании температуры на 100°С.

Значения параметров Km. Cvm. pm, а также генерации тепла, А и температурного коэффициента матричной теплопроводности, AI, могли быть определены по данным изучения петрофизических свойств пород конкретного разреза, но чаще оценивались на основании среднемировых данных для основных литологических единиц осадочных пород (Sclater and Christie, 1980; Beaumont et al., 1982; GofT, 1983; Hutchinson, 1985; Stockmal et al., 1986; Burrus and Anderbert, 1990; Deming and Chapman, 1989; Nielsen and Balling, 1990; Vasseur et al., 1995; Ilkusk, 1995; Correira and Jones, 1996 и т.д.). Теплофизические параметры пород, представленных комбинациями литологических единиц, вычислись в нашей модели с помощью специальных процедур осреднения в соответствии с (Doligez et al.. 1986; Ungereret al., 1990):

Pm ~ Pm ¡ Q + Pm2 Q+... + pmn Cn

Pm Cpm = pmlCpl C¡ + p„,2Cp2 C2 +... + PmncPn Qi (Ю) с с с

Km=Km¡ ' K1112 :.... Kmn "

Коэффициент температурного расширения матрицы пород AI в (9) определялся для смеси пород эффективно, через значения Km и Кпц при 100°С:

Al = 0.011 EXP ( £1п С i Ln(l + 100-Alj))-1 ] (11) где С i есть доля i-ой литологической компоненты в рассматриваемой осадочной породе. В согласии с (7) изменение теплопроводности осадочных пород с глубиной может быть значительным. Например, для песчаников теплопроводность меняется от значения 1.256 вт/м°С на поверхности до 3.2 вт/м°С на больших глубинах. Для глин аналогичные пределы изменения составляют от 0.92 вт/м°С до 2.1 вт/м°С. Вариации в значениях теплопроводности обуславливают заметные изменения геотермического градиента в пределах осадочной толщи бассейнов. Расчеты автора показли, что изменение теплофизических параметров пород с глубиной в процессе их консолидации существенно влияет на результаты моделирования термической истории осадков, особенно в режимах интенсивного осадконакопления или эрозии (Галушкин и др.,1985; Galushkin et al., 1991; Галушкин, Кутас, 1995; 1996; Makhous, Galushkin, and Lopatin, 1997).

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

параметров пород (К, Су, р, А) определялось в соответствии с данными геолого-геофизического изучения района. В системе ГАЛО предусмотрен выбор из четырех типов литосферы, первые три из которых были "стандартными" и включали литосферы континентального щита, мирового океана и окраинного моря с характеристиками, взятыми из соответствующих литературных источников (Смирнов, 1980; Ваег, 1981; Deming and Chapman, 1989; Rybach, 1996), а четвертый тип был "нестандартным" и характеризовался произвольными значениями параметров литосферы (числом слоев, их мощностями, плотностями, теплопроводностями и объемными генерациями тепла в слагающих породах). В согласии с (Ваег 1981; Сегшак, 1996; Rybach, 1996), вклад радиогенной составляющей пород литосферы в поверхностный тепловой поток увеличивается от Qpad = 10.5 мвт/м2 для океанической до Qpad = 13.-15. мвт/м2 для литосферы окраинных морей и Qpad = 20-30 мвт/м2 для континентальной литосферы. Этот вклад уточнялся при рассмотрении каждого конкретного района.

В термических расчетах системы ГАЛО учитывался эффект скрытой теплоты плавления, который рассматривался в рамках известного энтальпийного приближения (Carislaw and Jaeger, 1959). Доля расплавленной фракции, / по предположению, линейно увеличивалась с ростом температуры, Т, в пределах интервала Ts < Т < TI

__ Т~ Ts Tl -Ts

и теплоемкость Ср в этом интервале температур заменялась на значение С.р':

су=Ср*~— (12)

/7- is

где L - скрытая теплота плавления (100 кал/г для перидотит овых пород, Forsyth and Press, 1971), Ts - температура солидуса, которая являлась функцией давления (кривая 3 на рис.2д) и определялась по кривой солидуса для перидотита с 0.2% Н20 (Wyllic 1979):

Ts(lOOirC) =1.03 + 0.00391 Р для 0<Р< 12.8 Кбар Ts(1000"C) =1.236 - 0.0122 Р для 12.8 < Р <20.16 Кбар (13) Ts(1000°C) =0.742 +0.0122 Р для Р< 20.16 близкой к кривой плавления сухого лерцолита с содержанием НгО менее 0.5% (Hirose, 1997). 77 - температура ликвидуса: 77 = Ts + 450 - 600°С для пород мантии и 77 = Ts + 75-100°С для базальтовых пород (Turcotte and Schubert, 1982). Кривая солидуса (13) использовалась в модели для определения глубины подошвы литосферы (см. рис.2д, ж). Выделение и поглощение скрытой теплоты плавления замедляет процессы остывания и нагревания литосферы в периоды ее тепловой релаксации и активизации. Они особенно важны при анализе термического воздействия интрузий, вмещающихся в осадочную толщу бассейна (Galushkin, 1997).

В процессе моделирования вместе с температурой на каждом шаге времени вычислялось распределение плотностей пород в осадках и фундаменте. При расчетах плотностей пород в фундаменте учитывались изменения плотностей, вызванные 1) различием в строении и составе литосфер различного типа, 2) скачками плотностей при фазовых переходах в пределах мантии литосферы перидотитового состава, а также 3) изменения за счет термического расширении и изотермического сжатия пород:

/>, (Z, t) = р0 (Z, 1) • [1 - а ■ 1(2, t) + fi ■ P(Z, t)). (14) где а - коэффициент термического расширения (Parsons,Sclater 1977), ß -изотермический коэффициент сжатия пород (Touloukian and Но 1981), р0 (Z,t) -распределение плотности пород с глубиной при стандартных условиях (Р= 1 атм. и

Т=20°С). Параметр р0 включал вариации плотности, обусловленные типом пород (кора, мантия, "гранитный" или "базальтовый" слои, упомянутыми фазовыми переходами в литосфере, а также изменениями в глубинном распределении плотностей, возникающих при растяжении литосферы. Из фазовых переходов в литосфере в системе ГАЛО учитывались два основных перехода в пределах мантии перидотитовой литосферы (Forsyth, Press 1971): плакиоглазового перидотита в пироксеновый, контролируемого границей фазового перехода: Т(°С) = 194.17 • Р(кбар) ■ 761.16

со скачком плотностей от 3.26 до 3.30 г/см3 и пироксенового перидотита в гранатовый (рис. 2ж) с границей фазового перехода: Т(°С) - 83.94 ■ Р(кбар)-646.38

и со скачком плотностей от 3.30 до 3.38 г/см3 при средней плотности пород мантии рга=3.3 г/см3. Изменение в глубинах границ фазовых переходов по мерс развития бассейна могло обуславливать дополнительные 200-400 м в погружении или воздымании поверхности фундамента. В пределах континентальной литосферы первый фазовый переход часто не имеет место, так как по Р-Т условиям приходится на кору, тогда как в океанической литосфере и литосфере краевых морей первый фазовый переход по вкладу в рельеф поверхности подчас оказывается более существенным, чем второй, несмотря на меньший контраст плотностей.

Температура на верхней границе области счета (Z-0) определялась палеоклиматическими условиями, характерными для эволюции бассейна (см..например, рис.2а). Ее влияние на численные результаты реконструкции теплового режима бассейна обсуждаются в диссертации. Проблема выбора условия на нижней границе области счета не является простой задачей и подробно обсуждается в диссертации, где, в частности, отмечается что при реконструкциях теплового режима бассейна имеются весомые физические основания для задания в основании области счета температуры (а не теплового потока). При этом температура Tlow, поддерживаемая в основании области счета, определялась вместе с глубиной нижней границы Zlow пересечением предполагаемой геотермы самого холодного за время эволюции бассейна состояния литосферы с кривой солидуса пород мантии литосферы (15). В случае, когда такого пересечения не достигалось (холодная континентальная литосфера), глубина нижней границы области счета принималась равной 200 км. Начальное распределение температуры строилось решением стационарного аналога уравнения (6) для фундамента бассейна (линия 2 на рис.2д

Уравнение (6) решалось с использованием неявной конечно-разностной схемы, близкой к рассмотренной в (Peacemont and Rtachford 1954), но приспособленной к меняющимся значениям термофизических параметров пород и переменным шагам по времени и глубине (Галушкин, 1990).

(р-Ср-Т)к" -(р-Ср-Т)" (fW• Cpw• V• Т)^' -(pw• Cpw• V■ T)k*j =

dt dZk+dZh(.

к k.i (l5)

" dZhll.(dZK+dZkJ dZk(dZk+dZkJ где коэффициенты теплопроводности Kk'k\, и Kk*k_, осреднены но соседним слоям

dZk,dZk>i и dZk,dZk , на (п+1)-ом шаге времени:

dZt+dZM

dZt +dZt.t

dZk

X"

.and..KtJ,.t= ' ,;-'- (16)

Тридиагональная система уравнений для искомых значений температур в узлах решетки Т^1,1, Т^1*1.. и.. Т"+*, Т^.*,1 , дополненная граничными условиями на поверхности и в основании области счёта, решалась методом прогонки (см. например, Самарский, Гулин, 1989). В результате получалось распределение температур на следующем, (п+1)-ом, шаге времени. Шаг Az менялся по глубине как непрерывная кусочно-линейная функция от значений 0.5 - 20 м у поверхности осадочной толщи до 1000 -3000 м в основании области счета (z= 100-200 км).

Корректность разностной схемы проверялась сравнением численных решений с аналитическими (Карслоу и Егер 1964; Тёркот и Шуберт 1985) и полуаналитическими (Гольмшток 1979,1980; Галушкин и Смирнов, 1987) аналогами, а также сопоставлением решений, полученных с использованием различных шагов At и Az н подробно обсуждается в разделе 2.3 диссертации. Во всех случаях подходящий выбор шагов At и Az позволял обеспечить совпадение температур и тепловых потоков с результатами аналитических н полуаналитических методов с точностью до десятых долей процента. Как отмечается в диссертации, анализ режима эрозии требуег особой осторожности в выборе шагов At и Az (Галушкнн, 1988, 1990; Makhous, Galushkin and Lopatin, 1997)).

2.4 Тектоническое погружение, растяжение и термическая активизация литосферы

Известно, что в модели изостатического отклика литосферы на нагрузку амплитуда тектонического погружения поверхности фундамента складывается из двух основных составляющих (Sclater and Christie, 1980; McKenzie, 1981):

ZT=ZTs + ZTb (17) Здесь ZTs определяет относительное изменение глубины поверхности фундамента при удалении нагрузки воды и осадков и может быть записано в виде:

ZTs(t) - ZTs(0) = Р' ■ S(t) + Р-°-~Р" .[Z„ (t) - Z„ (0)] (18) Р, Р,

Вторая часть тектонического погружения,Z77>, связана с процессами, изменяющими распределение плотности в колонке фундамента (например, при нагревании и охлаждении пород фундамента, его растяжении или при смещении границ фазовых переходов в литосфере):

ZMt)-.ZTb(0)3JMK°) т

р, g

(Makhous, Galushkin, and Lopatin, 1997). В уравнениях (18, 19): t - время, t = 0 -начало формирования бассейна, Zw(t) - палеоглубины водной колонки на время t, ZTs(i) - амплитуда тектонического погружения, S(t) - толщина осадочного покрова, ра и рур- плотности астеносферы и воды, ps(i)- средняя плотность пород в столбце осадков, g - ускорение силы тяжести и G - вес столбца фундамента некоторой фиксированной высоты /о:

G(t) = gfp,(Z,t)dZ (20)

Плотность пород фундамента р[(Х,1) зависела от температуры, давления и состава литосферы в согласии с (14).

В наших моделях поверхность изостазии совпадала с основанием области счета и приходилась на глубинные горизонты хорошо прогретых реологически ослабленных пород астеносферы или низов литосферы, способных течь под действием сравнительно слабых разностей напряжений. В диссертации рассматривается этот вопрос и отмечается, что отклонения от состояния локальной изостазии в осадочных бассейнах ожидаются незначительными как на этапах ослабления литосферы бассейна (при ее аномальном прогревании или растрескивании при растяжениях фундамента), так и на региональном этапе развития бассейна, когда характерные горизонтальные размеры осадочного покрова заметно превосходят эффективную "упругую" мощность подстилающей литосферы (Галушкин и др., 1991). Для таких изостатически уравновешенных областей воздымание поверхности фундамента в отсутствии эрозии может вызвываться уменьшением плотности пород в колонке фундамента в ответ на прогревание литосферы бассейна или/и на погружение границ фазовых переходов (обрамление Красного моря, Байкала и др.). И, напротив, при относительно слабом воздействии прогревания уменьшение мощности "нормальной" коры в ответ на растяжение фундамента или как результат действия фазовых переходов, переводящих породы нижних горизонтов континентальной коры в более плотные гранулитовые и эклогитовые фации должно вызывать погружение поверхности фундамента. Рельеф района Афарского треугольника представляет пример одновременного действия интенсивной тепловой активизации и тектонического растяжения. В бассейнах с локально-изостатическим откликом литосферы на нагрузку сопоставление вариаций амплитуды тектонического погружения фундамента, вычисленных альтернативными методами: I) удалением поверхностной нагрузки воды и осадков (процедура "backstripping") и 2) из анализа распределения плотностей в фундаменте, - даст дополнительный метод контроля за параметрами тектонической и тепловой истории бассейна. Если некоторый этап развития бассейна характеризуется состоянием горизонтального сжатия, когда отклонения от локальной изостазии становятся заметными, то в расчетах тектонического погружения необходимо вводить поправочный множитель у<1, учитывающий уменьшение отклика литосферы на нагрузку от изостатического за счет конечной жесткости литосферы бассейна, поддерживаемой в том числе и региональными напряжениями. Более сложная ситуация будет в районах проявления активной тектоники, типа предгорных прогибов, когда появляется заметная динамическая составляющая тектонического погружения.

Как показывают исследования, процессы тепловой активизации и растяжения фундамента (последние, как отмечается в диссертации, могут включать в себя и сокращение мощности коры при фазовых переходах в ее основании) могут неоднократно повторяться в истории бассейна. Метод численного анализа влияния этих процессов на тепловую историю бассейнов обсуждается в разделе 2.4 диссертации. Рис.2 показывает пример влияния двух периодов медленного растяжения с амплитудами около 1.2 и продолжительностью около 100 и 50 млн. лет (PACI и РАС2 на рис.2е), а также двух тепловых активизаций с эффективными теловыми потоками на поверхности около 100 и 60 мВт/мг (ТАК1 и ТАК2 на рис.2е, ж) на термическое состояние литосферы бассейна Уэд эль Миа в Восточном Алжире.

2.5 Методы контроля палеотемпературной истории бассейна

В данном разделе диссертации обсуждается вопрос о реперах палеотемпературной истории бассейна. Раздел начинается с рассмотрения отражательной способности

внтрннита Ro% как наиболее распространенного и обоснованного метода контроля степени катагенеза ОВ, обсуждаются физические и геохимические условия применения этого метода. К его преимуществам, помимо доступности и воспроизводимости, относится возможность теоретического расчега значений Ro по палеотемпературной истории образца. В системе ГАЛО для вычисления значений Ro использовалась кинетическая модель витринита, разработанная в Ламонтской лаборатории США (Sweeney and Burnham, 1990). В этой модели процесс созревания витринита описывается 20-ю реакциями первого порядка Аррениуса со значениями частотного фактора А и исходными потенциалами реакций Хю, приведенными в табл.7-2 диссертации. В процессе моделирования на каждом шаге времени и в каждом узле решетки в пределах осадочной толщи вычислялась степень преобразованности витринита, и затем значения отражательной способности витринита определялись из корреляционного соотношения (Sweeney and Burnham, 1990). Сравнение вычисленных значений Ro с величинами, измеренными в современном разрезе, использовалось в системе ГАЛО для корректировки исходных параметров модели (рис.2в). Значения температурно-временного индекса (ТВИ, Лопатин, 1971; Waples, 1980) также рассматривались в наших моделях, однако они использовались лишь как дополнительный, простой, но приближенный метод характеристики палеотемпературного режима осадков.

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

2.6 Моделирование истории реализации углеводородного потенциала нефтематеринских свит бассейна

Процессы созревания керогена и история генерации УВ рассматривались в системе ГАЛО с использованием алгоритмов, представленных в (Tissot and Espitalie, 1975; Tissot et al., 1987; Espitalie et al.,1988; Ungerer, 1990; Issler and Snowdon, 1990). Предполагалось, что созревание OB нефтегазоматеринских пород можно описать с помощью ряда независимых реакций Аррениуса первого порядка

Ki(t) = Ai exp ( -Ei/R T(t) ) (21) где Ki - скорость реакции, Ai - частотный фактор Аррениуса, Ei - энергия активации, R - газовая постоянная, t - время, и Т - температура в °К. Если napaMerps Ai, Ei и Xio (исходный потенциал генерации УВ для i-ой реакции) известны для всех N реакций, то количество генерированных УВ и скорость их генерации вычислялись по формулам (Tissot and Espitalie, 1975):

Q(t)= ^Xio-[l-exp(-JlKi(t')-dt'].

' N <0 (22)

S2(t) = ^ = S Xio • Ki(t) - exp(-|'o Ki(f) ■ df)

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

экспериментального пиролиза образцов нефтематеринских пород на установках типа Rock-Eval был разработан специальный программный модуль , включенный в систему ГАЛО. В этом модуле кинетические параметры эффективных реакций £/, Ai и Xio находятся с использованием нелинейного алгоритма оптимизации (метод Левенберга-Макуиртера, Press et al., 1986; Issler and Snowdon, 1990) в сочетании с линейным методом поиска коэффициентов. Решается обратная задача статистики, в которой минимизируется функция ошибок % ПРИ сравнении скоростей генерации углеводородов, S2, вычисленных по уравнению (22), и измеренных в экспериментальной установке открытого пиролиза для набора температур нагревания образца T¡ в пиролизаторе:

*2(a)=Ír-S*-S*(T"a)]2 (24) Tt

Здесь а, - ошибка /-ого измерения, N - число измерений, а - искомый вектор кинетических параметров, включающий М искомых исходных потенциалов реакций, Xi, и М искомых частотных факторов реакций, A¡, при заранее выбранных значениях энергий активации Ег . Размерность искомого вектора а могла достигать 80-ти. Применяя метод градиентного спуска к решению обратной статистической задачи (Press et al.,1986), мы получали отклонение 5ага ш-той компоненты вектора а от ее предыдущего приближения решением системы уравнений:

(25)

где

„ -о, I? V 1 ^г(Т„а)

" Km ~ "-J'. ¡'1 _ ' 1

¿hK hl a ¿fek

(26)

и вычисляется согласно (22). В отличие от зарубежных аналогов наш алгоритм восстановления параметров спектров реакций не ограничивался вариантом с одним значением частотного фактора А, и это позволяло осуществлять более качественный подбор кинетических параметров (сравни рис.За, б -один фактор А, с кривыми на рис Зв, г - различные А|).

Нередко встречается ситуация, когда образцы материнских пород, подвергаемых экспериментальному пиролизу, характеризуются заметным уровнем катагенеза ОВ (Яо»0.50-0.70%), приобретенным при погружении образца в геологическую историю бассейна. В нашем алгоритме учтен этот фактор и, в отличие от французского, германского и американского пакетов, он включает геологический этап созревания ОВ пиролизируемого образца. Для этого интеграл по времени в (22) представлялся в виде суммы двух интегралов:

(к«ГНГ= (27)

где ю < !'< // есть временной интервал истории погружения образца в бассейне и (/</'</- временной интервал программируемого нагревания образца в пиролизаторе. При известной из реконструкции температурно-временной истории материнских пород (рис.Зе или 2з) такой подход позволял восстанавливать спектры эффективных кинетических реакций по данным пиролиза образцов с исходными значениями Яо < 0.70-0.80% с заметным уточнением распределения энергий в низкоэнергстичной области спектра (сравни спектр на рис.Зв,

g 15-

& 10-

-0.3

35 45 55 85 75 -120 -t)0 -00 -ho'

" _ Энергия (Кхал/моль) Время (unH.mi)

Рис.О Кинетические спектры: a) j) —<5еэ учета геологической истории, д) — с учетом геологической истори'1 созревания образца, представленной на рис.Зе. На рис.36, Зг — скорости генерации У13 в эксперименте, (сплошные линии—вычисленные, символы — измеренные значения)

полученный без учета геологической истории образца, со спектром на рис.3 д, учитывающим эту историю).

Формулы (21,22), позволяющие восстанавливать историю реализации УВ потенциала пород по известным кинетическим параметрам реакций созревания ОВ, предполагают, что процесс крекинга керогена осуществляется через реакции первого порядка. При наличии реакций вторичного крекинга жидких УВ на газовые и коксовую составляющие, имеющих место при температурах Т > 140-150°С, выход различных УВ-составляющих (в отсутствии заметной эмиграции УВ из слоя) будет определяться уравнениями (Tissot et al.,1987; Espitalie et al., 1988):

dq/dt = С X + В q (28) с начальным условием: qi(0)=Copr(0)q,o, j=l,...,n. Матрица Q=Ai-EXP(-Ei/R-T)a,j описывает первичный крекинг, X - массив остаточных потенциалов первичных реакций (dxi/dt=-A, EXP(-iyR-T) x,(t)). Матрица В описывает процессы вторичного крекинга неустойчивых фракций УВ [B,j=A'j-EXP(-E'j/R-T) aji для i=I,...,m, j=l,...,p и Bjj=A'j EXP(-E'j/RT) для j=l,...,p и В^=0для i=l,...,m, j=p+l,...,m]. -а^- вклад i-ой фракции в продукте j-ой вторичной реакции, р - число неустойчивых фракций (число вторичных реакций), ш - полное число фракций, n-число первичных реакций, Ау,Е']-предэкспоненциальный фактор и энергия активации реакций вторичного крекинга. Для трехкомпонентной системы (нефть-Qou, ra3-Qgas и кокс-Qcoke) соответствующие выходы УВ как функции времени можно представить в замкнутом виде:

^EXPi-jKj,-сй'Л/к',,, C'k EXPijK^ df)-df] 11 11 1. t„ !„

QB„ ¿xg1.[l-EXP(-jKlgdt') (29)

t„ И и t.

Q«*. = ¿0 ~ xri)-JK^ -Q'ou -dt'

i

где C'h = xo, • EXP(-jK'd| dt')- концентрация i-ой составляющей керогена, >»

преобразующейся в жидкие УВ, K»ii' - скорость этой реакции, Kvr' - скорость крекинга i-ой составляющей жидких УВ, хп - доля выходи газовых УВ в этой реакции, xOj, хц - исходные потенциалы i-ой и j-ой реакций первичного крекинга керогена. Формулы (29) применялись для расчетов относительных вкладов жидкой, газовой и коксовой составляющих в полный выход УВ (Галушкин и др.,1995; Галушкин, Кутас, 1995; Makhous, Galushkin and Lopatin, 1997). В расчетах использовались кинетические параметры реакций крекинга керогенов стандартных типов I, II, III (Tissot et al.,1987; Espitalie et al.,1988) на 3 (нефть, газ и кокс) и 5 (С[, С2-С5, Cg-C[5, С15+ и кокс) фракций (см. раздел 2.6 диссертации).

Система моделирования бассейнов ГАЛО включала также оценку порога первичной эмиграции жидких УВ (рис.2 и) и количества эмигрировавших УВ (Лопатин, Галушкин и др.,1992; Lopatin, Galushkin and Makhous, 1996). В расчетах учитывались изменения пористости и плотности материнских пород с глубиной и временем, а также уменьшение содержания Сорг по мере созревания ОВ в них. Как отмечается в дисссергации, оценки порога эмиграции очень критичны к значению пороговой пористости, Plh, которая представляет долю порового пространства, заполненного жидкими УВ, прежде чем начнется их эмиграция из материнских пород. Величина Plh может меняться от 5 до 60% в зависимости от типа ОВ и литологии вмещающих пород (Espitalie et al.,1988; Ungerer, 1990).

2.7 Моделирование термического режима бассейнов в условиях похолодания климата в плиоцен-четвертичное время с учетом формирования криолитозон.

Вариации климата в плиоцен-четвртичное время вызвали существенные отклонения термического режима осадочных толщ от стационарного в бассейнах высоких и средних широт Северного полушария. Формирование и деградация горизонтов вечномерзлых пород неоднократно имели место в течение всего рассматриваемого периода. Методы, позволяющие оценить влияние этих процессов на термическое состояние осадочной толщи бассейнов высоких широт, обсуждаются в разделе 2.7 диссертации. Показано, что при таких обстоятельствах значительную роль в формировании температурного режима бассейнов играют изменения теплоемкости и теплопроводности пород, связанные с процессами замерзания воды и таяния льда в пористых осадочных породах, а также с вариациями литологии пород в реальном разрезе бассейна. Для мерзлых пород теплоемкость и теплопроводность в заметной степени зависят от содержания незамерзшей воды и льда в породе (Кудрявцев, 1981; Lachenbruch et al.,1982)

„ -, (1-ф(1)) „ M.W(T)Í ¿íz)-(l-W<TH

К = Km ■ Kw' ■ Ki

Cv = Cvm[l-<t,(z)]+Cvw<t>(z) W(T)+Cviifi(z) (l-W(T))+pLdW(T)/dT (30) -a для талых пород они определяются формулами (9). В (30): z-глубина, ф-пористость, Т-температура, Krn-теплопроводность минерального скелета осадочной породы, Kw-тештопроводность воды, изменяющаяся с температурй, • Кл=2.26 вт/м°К-теплопроводность льда, W-относительная доля воды, не замерзающей при температуре Т, TL-температура ликвидуса льда или начала замерзания воды, зависящая от давления и содержания солей в грунтовых водах (Nixon, 1986):

TL = 0°C - 0.073 Р -0.064 Cs (31) Р -давление поровои жидкости в МРа, предполагаемое гидростатическим, Cs -концентрация солей NaCl и КС1 в г/литр. Теплопроводность льда почти в 4 раза превосходит среднюю теплопроводность воды, поэтому для мерзлых пород теплопроводность сильно зависит от вида кривой льдистости, W(T). Последний член в формуле (30) для Cv описывает вклад скрытой теплоты плавления льда в энтальпийном приближении, когда температуры фазового перехода "вода-лед", а вместе с ними и тепловой эффект скрытой теплоты плавления, оказываются распределенными по интервалу температур, определяемому формой функции W(T) (Jame and Norum, 1980; Кудрявцев, 1981; Nixon, 1986; Ершов, 1990). Такой подход несомненно ближе к описанию реальных процессов в мерзлых породах, нежели более популярный классический метод Стефана с фиксированной температурой перехода воды в лед для всего обьема жидкости в породе. Эффективный температурный интервал перехода воды в лед (и обратно) определяется видом кривой W(T). В зависимости от формы кривой эффект скрытой теплоты плавления может увеличивать эффективную теплопроводность мерзлых пород на 1-2 порядка (Galushkin, 1997).

На основе рассмотренных алгоритмов автором был разработан пакет программ "ФРОСТ", позволяющий проводить численный анализ эволюции температурного режима бассейнов в условиях резких колебаний климата. Пакет создавался как один из прикладных вариантов общего пакета программ ГАЛО по моделированию палеотемпературной истории бассейнов и оценке их углеводородного потенциала. Существенно, что изменение теплового режима осадочной толщи в течении последних 3-4 млн. лет моделировалось как продолжение эволюции бассейна на более ранних этапах его развития, начиная с

момента зарождения бассейна. Соответствующие расчеты, проведенные на примерах отдельных площадей Западно-Сибирского и Печерского бассейнов (Галушкин, Лопатин, 1997, Galushkin, 1997), позволили уточнить оценки нестационарного термического эффекта вариаций климата, которые в литературе ограничивались в основном относительно короткими интервалами времени в позднечетвртичный период и осуществлялись в рамках относительно простых моделей, включавших два однородных слоя талых и мерзлых пород (Баулин, 1962; Кудрявцев, 1981; Lachenbruch et al., 1982; Сигунов, Фартышсв, 1991, 1995; Lebret ct al., 1994). Расчеты, представленные в диссертации, позволили оценить вариации мощности зон вечномерзлых пород, также как и зон устойчивого существования метановых газогидратов. Они показали, что в силу относительной кратковременности процессы формирования криолитозон не оказали заметного влияния на историю созревания ОВ осадков, но их влияние на современное распределение температур значительно. Последнее, как известно, служит важным контролирующим фактором при построении модели. По этой причине для бассейнов северных широт анализ температурного режима криолитозон является обязательной частью процедуры моделирования бассейнов.

2.8 Численный анализ термической эволюции и созревания ОВ осадков под влиянием тепла интрузий.

Этот раздел диссертации посвящен методам численной оценки размеров и интенсивности ореола термического влияния интрузии. Мировые данные согласуются с тем фактом, что термическое влияние интрузии простирается, как правило, на расстояния 30 - 60 % толщины интрузии и очень редко превосходит 100% толщины. Они также налагают ограничения на возможные максимальные температуры вмещающих пород в зоне контакта с интрузией значениями, которые могут быть на 200-500°С ниже исходной температуры жидкой магмы. Ограниченный ореол термического воздействия интрузий вместе с относительно низким уровнем метаморфизма пород контактной зоны могли быть обьяснены лишь в рамках модели, отличной от традиционной модели мгновенного вмещения расплавленного магматического тела во вмещающие породы.

Для решения проблемы была разработана специальная модификация пакета моделирования бассейнов, позволявшая проводить вычисления эволюции термического режима осадочной толщи и фундамента вместе с анализом температурного режима самих интрузий и вмещающих пород с одновременным расчетом степени метаморфизма ОВ вмещающих пород. Существенно, что вычисления проводятся для реального осадочного разреза, содержащего интрузию. Специальная процедура обеспечивала разбиение области счета в окрестности интрузии на достаточно мелкие шаги по глубине и времени, гарантируя устойчивость решения и выполнения условия, что смещение фронта плавления за шаг времени было в 2-3 раза меньше соответствующего шага по глубине. Модель позволяла учитывать тепловое и катагенетическое влияние гидротермального теплопереноса выше интрузии, а также тепловой эффект от реакций дегидратаци во вмещающих породах (Galushkin, 1997).

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

моделях мощность интрузивного тела наращивается в течении некоторого промежутка времени по механизму "дайка в дайке". Алгоритм вычисления был близок к использованному нами при анализе процесса формирования подосевой коровой магматической камеры срединно-океанических хребтов (Галушкин и Дубинин, 1993; 1994). Результаты расчетов по этой модели, представленные в диссертации, демонстрируют хорошее совпадение вычисленных и наблюдаемых термических ореолов интрузий (Са^Нкт, 1997).

Глава 3. ТЕРМИЧЕСКАЯ ЭВОЛЮЦИЯ ПЛАТФОРМЕННЫХ БАССЕЙНОВ И

УСЛОВИЯ РЕАЛИЗАЦИИ ИХ НЕФТЕГАЗОГЕНЕРАЦИОННОГО ПОТЕНЦИАЛА

3.1 Термический режим осадочной толщи и условия созревания ОВ на рифтовой стадии развития бассейнов

Высокий тепловой поток на рифтовой стадии развития бассейна способствует раннему созреванию ОВ в осадках. Степень преобразованности ОВ при этом существенно зависит от интенсивности осадконакопления. Этот факт иллюстрируется результатами численного моделирования процесса остывания литосферы с начальным тепловым потоком на поверхности <3=105 мВт/м2, типичным для современных рифтовых областей, и с одновременным отложением песчанисто-глинистых осадков со скоростями от 0.1 до 1.0 км/млн.лет в течении первых 10 млн лет остывания с последующей релаксацией термического режима при неизменной мощности осадочной толщи. Расчеты показывают, что при максимальной скорости осадконакопления У=1 км/млн.лет породы в основании осадочной толщи уже через 5 млн. лет после начала осадконакопления достигают температур свыше 150°С при степени катагенеза пород Ио > 1.0%. Такие же условия для скоростей осадконакопления V > 0.4 км/млн.лет достигаются в пределах первых 10 млн. лет. Моделирование, проведенное для трех стандартных типов керогена (кероген типа II с исходным потенциалом Н1=627 мг УВ/г Сорг, кероген типа II с более бедным потенциалом Н1=377 мг УВ/г Сорг и кероген типа III с исходным потенциалом Н1=!60 мг УВ/г С^,,), показывает, что достигнутые температуры и уровни зрелости соответствуют началу вторичного крекинга жидких УВ при условии, что эти УВ не имели возможность мигрировать в верхние горизонты с более умеренными температурами. В этих условиях жидкие УВ, генерированные в основании накопленной толщи осадков, разлагаются на газовые составляющие и кокс в течение первых 7-11 млн лет, а исходный потенциал генерации УВ этих пород будет реализован за это время почти полностью. Другая ситуация, как показано в диссертации, будет иметь место для более умеренных скоростей осадконакопления и соответственно меньших мощностей отложений.

Тесная связь между скоростью осадконакопления и степенью катагенеза ОВ в породах на рифтовом этапе развития бассейна подтверждается и результатами моделирования для осадочных разрезов Днепрово-Донецкого (Галушкин и др., 1985; Галушкин, 1988; ва^Икт е! а1.,1991), Западно-Сибирского (Галушкин и др., 1994, 1997), Северо-Германского (ВегНюШ ипс! Са^Ькт, 1986; 1988) и других бассейнов. Во всех случаях, когда мощность осадочного покрова, отложившегося в условиях повышенного теплового потока, превосходила 2 км, достигался значительный уровень уровень катагенеза ОВ и реализации исходного нефтегазогенерационного потенциала глубокопогруженных пород уже на рифтовом этапе развития бассейна.

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

миллионов лет) процессов тепловой активизации литосферы, сопровождающих рифтогенез. Так, для Нижнепурского прогиба Уренгойского вала ЗападноСибирского бассейна анализ предполагает заметную тепловую активизацию литосферы с эффективным тепловым потоком на поверхности 65-85 мвт/м2, которая охватывала не только триас, но и подавляющую часть нижней юры (Галушкин и др., 1994, 1997; Galushkin, 1997). Предполагаемая из анализа тектонического погружения тепловая активизация литосферы в течение 20 млн. лет рифтового этапа развития бассейна помогла обьяснить быстрое пострифтовое погружение Днепрово-Донецкого бассейна в каменоугольный период (Галушкин и Кутас, 1995), причины которого были предметом оживленных дискуссий (Артюшков, 1993; Nikishin et al., 1996, 1997). Продолжительная тепловая активизация на рифтовом этапе эволюции бассейна может сопровождаться внедрением силлов и гидротермальной циркуляцией в осадочной толще и верхних горизонтах фундамента. Учет этих процессов в моделировании позволяет обьяснить резкое увеличение отражательной способности внтринита в нижних горизонтах осадочных толщ, типичное для многих рифтогенных бассейнов мира (Галушкин и др.,1994, 1997; Galushkin, 1997).

3.2 Анализ эволюции температурного режима и условий созревания OB в пострифтовый этап развития бассейнов

Интенсивность осадконакопления остается важным фактором формирования термической истории осадочных пород и на пострифтовом этапе развития бассейна. Периоды медленного осадконакопления в отсутствии тепловой активизации фундамента характеризуются относительно умеренным ростом уровня катагенеза OB, как это демонстрируют численные реконструкции разрезов Ромашкинского месторождения Волго-Камской антеклизы (Смирнов и Галушкин, 1986), отдельных площадей Днепрово-Донецкого авлакогена (Галушкин и др.,1985; Galushkin et al.,1991), Пипятского бассейна (Berthold, Galushkin und Muller, 1986), ряда разрезов Северо-Германской впадины (Berthold und Galushkin, 1986; 1988) , Печерской синеклизы (Смирнов and Галушкин, 1986), бассейна Иллизи в Алжире (Munzer, Galushkin, and Lopatin, 1997) и других.

Напротив, интенсивное осадконакопление может перемещать глубинные осадки в область больших глубин и температур, и заметно увеличивать уровень катагенеза OB. Так, отложение 6 км терригенных песчано-алевролито-глинистых фракций в триасе в нашем примере моделирования бассейна Южно-Баренцевой впадины в районе Штокмановского месторождения имело итогом увеличение температуры пород нижнего триаса до 150°С и уровня катагенеза до Ro= 1.2-1.4% всего лишь за первые 40 млн.лет развития бассейна. Как следствие, УВ потенциал этих пород был реализован еще на юрском этапе развития бассейна со вторичным крекингом жидких УВ, не успевших эмигрировать из пород (Галушкин и др., 1994). Рост уровня катагенеза при высоких скоростях осадконакопления имел место и в Восточно-Черноморской впадине, где интенсивное осадконакопление в позднем кайнозое сопровождалось значительной тепловой активизацией фундамента (Смирнов и Галушкин, 1986; Galushkin et al.,1991; Spadini et al.,1996).

Эрозия относится к числу распространенных процессов в истории бассейнов. Как показало моделирование, интенсивная эрозия, перемещая породы к поверхности, в область более низких температур, может законсервировать уровень катагенеза OB в них на величине, близкой к доэрозионной (пермская эрозия силурийских пород на площади Такхухт бассейна Оуэд-эль-Миа в Алжире -рис.2б,з; Lopatin, Galushkin, and Makhous, 1996; Makhous, Galushkin, and Lopatin, 1997). Считают, что эрозия всегда имеет результатом повышение теплового

потока и температурного градиента в среде (Череменский, 1972). Однако, моделирование показало, что это справедливо лишь для эрозии однородной по теплофизическим свойствам среды, а для реального осадочного разреза с термическими характеристиками пород, меняющимися с глубиной, умеренная эрозия может приводить к уменьшению температурного градиента в осадочной толще (Смирнов и Галушкин, 1986).

3.3 Роль тепловой активизации н интрузивной деятельности в бассейнах

Тепловая реактивизация фундамента довольно распространенное явление в термической истории многих осадочных бассейнов. Так, в наших моделях эволюции бассейнов рассматривались: девонская тепловая активизация в Печорском бассейне (Смирнов и Галушкин, 1986; Galushkin et al.,1991), пермская активизация в Северо-Германском бассейне (Berthold und Galushkin, 1986; 1988), юрская и кайнозойская в Западно-Сибирском (Lopatin, Galushkin, and Makhous, 1996; Галушкин и др.,1994; Галушкин и др.,1997), пермско-нижнекарбоновая в Днепрово-Донецкой впадине (Галушкин и Кутас, 1995) и другие. Как правило, тепловая активизациия литосферы бассейна приводит к заметному увеличению температуры осадочных пород и уровня катагенеза ОВ в них. Исключения бывают, когда активизация сопровождается интенсивной эрозией, перемещающей глубинные породы к поверхности, в область низких температур, как это имело место во время пермской тепловой активизации в бассейне Уэд эль-Миа в Алжире. Здесь мощная эрозия привела к заметному уменьшению температуры силурийских материнских пород, что имело результатом лишь незначительное увеличение их катагенеза (рнс.2 б, ж; Lopatin, Galushkin, and Makhous, 1996; Makhous, Galushkin, and Lopatin, 1997). В то же время кайнозойская тепловая активизация литосферы бассейна не сопровождалась эрозией и привела к заметному росту температуры и уровня катагенеза пород (рис.2 б, ж).

На отдельных этапах своего развития литосфера бассейнов в результате действия региональных полей напряжений может подвергаться растяжению (рис.2е). Результаты моделирования, приведенные в главе 3 диссертации, показывают, что амплитуда растяжения, оцениваемая по вариациям кривой тектонического погружения, редко превышает значение 1.20, в то время как продолжительность процесса, как правило, превосходит 30-40 млн.лет (глава 3 диссертации). По этой причине скорости деформаций пород при растяжении литосферы бассейна остаются незначительными, а вместе с ними малы и термические следствия процесса растяжения (рис.2б, е, ж; Makhous, Galushkin, and Lopatin, 1997).

В некоторых бассейнах (Восточно-Сибирский в России, Парана и Парнаибо в Бразилии) тепловая реактивизация литосферы сопровождалась обильной интрузивной активностью с прогреванием пород и созреванием ОВ в окрестности интрузий. В разделе 3.3 диссертации методика моделирования термического воздействия интрузий, описанная в главе 2, применяется для анализа термической эволюции и процесса созревания ОВ во вмещающих породах ряда континентальных интрузий с мощностями тел от 0.9 до 118 м. На примере внедрения 5-ти юрских и меловых интрузий (мощностью до 200 м) в осадочный разрезе бассейна Парнаиба в Бразилии анализируются тепловые следствия внедрения роя интрузий в осадочный чехол бассейнов. Результаты моделирования подтверждают тот факт, что внедрение интрузий в осадочный чехол бассейнов может существенно повысить уровень катагенеза ОВ вмещающих пород на расстояниях от 0.5 до 1.0 ширины интрузии. Подтвердились также выводы раздела 2.8 диссертации о том, что вычисления в рамках традиционной модели мгновенного внедрения интрузии заметно переоценивают температуру и уровень

метаморфизма вмещающих пород. Рассмотренные для этих примеров модели, предполагающие наращивание интрузивных тел в течении определенного промежутка времени (от минут до полугода), и, возможно, без прямого контакта расплавленных и вмещающих пород, удовлетворительно объясняют наблюдаемые размеры ореолов зрелости интрузий. Процесс генерации УВ в окрестности интрузий оказывается весьма кратковременным (от минут до нескольких лет), поэтому остается актуальным вопрос о конкретном механизме взрывной эмиграции жидких УВ из зоны контакта с магматическим телом и формирования для них подходящих ловушек.

3.4 Реконструкция теплового режима и истории реализации потенциала генерации УВ типичных континентальных рифтогенных бассейнов

В разделе диссертации рассмотриваются реконструкции теплового режима и истории генерации УВ для осадочных разрезов Западно-Сибирского, Днепрово-Донецкого, Южно-Баренцева бассейнов России и бассейна Уэд эль Миа в Восточном Алжире, что позволило сделать выводы о влиянии некоторых специфичных процессов на формирование термического режима и эволюцию катагенеза ОВ осадочных бассейнов.

Так, моделирование на примере осадочного разреза Уренгойского месторожденя предполагает, что при содержании ОВ в рассмотренном примере Сорг=1-3% в меловых и юрских горизонтах эффект на тепловое состояние бассейна от снижения теплопроводности пород за счет органических примесей не превышает 3-5°С по температуре и 0.10-0.12% по Яо. Более значительным оказалось влияние уменьшения теплопроводности , вызванное появлением в осадочном разрезе значительных пропластов со свободными УВ. Для осадочного разреза скв.266 соответствующее повышение температуры пород достигало 10-13°С на глубинах 3000 - 5700 м, однако, этот эффект не превосходил 3°С для разреза скв. 411, где аккумуляция УВ была заметно меньшей.

Реконструкция термического режима Западно-Сибирского бассейна позволила нам восстановить химико-кинетический спектр генерации УВ пород покурской свиты, воспользовавшись данными экспериментального пиролиза образца покурской свиты на установке с хорошем контролем температур. Процедура восстановления спектра была аналогичной описанной в главе 2 и учитывала созревание ОВ образца на этапе его погружения в бассейне (рис.3д,е), Результирующий спектр предполагает, что термический крекинг ОВ пород покурской свиты может вносить заметный вклад в формирование больших месторождений метана Уренгойской площади.

Па примерах баженовской свиты Западно-Сибирского бассейна и силурийской свиты бассейна Уэд эль-Миа в Восточном Алжире моделирование подтвердило тот факт, что аномальное тепло радиоактивности, выделяющееся, в сравнительно узких слоях нефтематеринских пород (100-150 м) не оказывает заметного влияния на формирование температурного режима пород и процесс созревания в них ОВ.

В главе 3 на примере реконструкции теплового режима ЗападноСибирского бассейна было оценено также влияние резких изменений климата в плиоцен-четвертичное время на температурный режим осадков и фундамента бассейнов высоких широт. Результаты моделирования существенно зависели от конкретной литологии изучаемого разреза. Они предполагают, что вариации климата в течении последних 3.4 млн. лет в Западно-Сибирском осадочном бассейне привели к понижению температур в верхних 1500 м осадочного разреза на 15-20°С и на 8-10°С в его нижних горизонтах. В результате последнего

оледенения (23-15 тысяч лет назад) максимальная глубина погружения границы мерзлых пород составляла около 450 м и наблюдалась около 13 тысяч лет назад. В настоящий момент, согласно расчетам, зона нермафроста Уренгойской площади деградирует со скоростью около 6 м/1000 лет и ее мощность в местах Уренгойской площади, не подверженных действию аномальных факторов (озера, долины рек), должна достигать 350 м (Оа^яЬкт, 1997).

Вместе с зонами вечномерзлых пород в периоды резких колебаний климата происходило формирование и деградация зон, благоприятных для устойчивого существоания газогидратов. Моделирование предполагает, что за всю анализируемую историю климата на изучемой площади Уренгойского месторождения нижняя граница устойчивого существования газогидратов не опускалась ниже 800 м, всегда оставаясь выше кровли сеноманского газового резервуара. Тогда один из основных механизмов образования газогидратных залежей внутри зон скоплений свободного газа становится не актуальным для изучаемой площади, за исключением, быть может, сравнительно узких зон активной вертикальной миграции газа из сеноманских залежей (Галушкин, Лопатин, 1997).

Глава 4. ПРОСТРАНСТВЕННО-ВРЕМЕННЫЕ МОДЕЛИ ЭВОЛЮЦИИ ТЕРМИЧЕСКОГО РЕЖИМА ЛИТОСФЕРЫ РИФТОГЕННЫХ ОСАДОЧНЫХ

БАССЕЙНОВ

В главе 4 диссертации анализируются те аспекты формирования термического режима литосферы рифтогенных бассейнов, которые не могут быть изучены в рамках моделей плоского бассейна, применявшихся в главах 2 и 3, и требуют пространственно-временных методов численного анализа.

4.1 Методы анализа термической эволюции, теплового потока, рельефа поверхности и гравитационного поля подстилающей литосферы континентальных и океанических осадочных бассейнов

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

+ ^|(рСр-т.^) + 4(Р.ср.т.у2) = 4(к4т)+|:(к4т)+А (32>

¿Х м ¿2 ¿х м дг. дг.

где р- плотность, Ср - теплоемкость, К - теплопроводность, А -теплогенерация, Г-температура пород литосферы, V, - х и г-составляющие скоростей

перемещения вещества коры и мантии. Все параметры в уравнении являлись функциями координат (х, г) и времени (г)- Скрытая теплота плавления, игравшая заметную роль в формировании теплового режима литосферы, учитывалась при решении (32) энтальпнйным методом, обсуждавшимся в главе 2. Для задач, рассматриваемых в этой главе, граничные условия предполагали поддержание определенных температур на верхней (г=0.) и нижней (7= 100-200 км) границах области счета, условие сТ/сЬс- 0 на достаточном удалении от оси рифта (или спрединга) и поддержание заданного высокоградиентного распределения

температур на оси области х=0. во время континентального рифтогенеза или спрединга или условие симметрии здесь {дТ/дх-0) в периоды остывания литосферы. Уравнение (32) решалось по неявной конечно-разностной схеме близкой к рассмотренной в (Peacemont and Rtachford 1954), но приспособленной к изменяющимся значениям термофизических параметров пород и переменным шагам по времени и глубине и использующую схему с опережением для аппроксимации конвективного члена: ОСЬ X:

[ (Р-Ср-УхП 2-е, ] +

Ах, Дх,(Ах, + Дх,и) м,к

.(р-уС , (P-Cp-V.c | 2-Kgj, | 2-Kg?, п<| +

Atn„ Дх, Дх,(Дх,+Ах,(1) ДхИ1(Дх, + Дх,(1) *

I-iJlHil-]. тп;' = [--.]. Т", +

,(Р-Ср-); | (р-ср-Уг)Г, 2-К<"к_,___2-Ка„ jn + Ш)

At„„ Дzkt1 Azk(Azk+Azk<1) Azki,(Azk + Azkf,) lh

. , 2-km»I , Tn ,nn

Azktl +Azktl(Azk+AzkJ]T'-+A^ .............................................................for..2<i< JM -),..! <k< KM

ОСЬ Z:

2 • K0)

Azk(Azk+Azk<1)J tt"'

|f(PS-r (Р'УЧГ , 2-KS., | 2-Kg>t, ^

Atn+2 Azk), Azk(Azk + Azktl) Azkt,(Azk+A2k„)J w

|f(p-yv,C 2-k'',L, nt2 (pcp-v,n 2 -Kg', r„

(P-Cp-)";' (Р-ср-У,С 2-Kg!, 2-C, +

Atnt2 Ax, Ax,(Ax, + Ax,J Ax,f,(Ax,+AxM) "

9 K(k)

+[ ""—-] • T/:;; + A-2

for.2 < k < KM - l,..l < i < JM

к«> _ Лгк + А2к„ ,„ _ Агк + Агк_, (к) _ Лх, + Ах,„ (к) _ Ах, + Ах,,

= " кк"' = + Л^,.....= Дх, + Ах,,1 = Дх, + Дх,/

к,к4, К№_, К№ К,л К,_1к

Тридиагональные системы уравнений для решения (32) методом альтернативных направлений включали также граничные условия на оси и боковой границе

области Х=0 (1=1) и Х=ХМ у=.)М) по оси X и на поверхности Х=0 (К=1) и в основании области счёта 2=2М (К=КМ) по оси Ъ. После решения уравнений (39) методом прогонки получалось распределение температур на следующих (п+1)-ом и (п+2)-ом шагах времени. Точность используемой разностной схемы проверялась сравнением с аналитическими решениями для распределения температур в остывающем однородном полупространстве (Карслоу, Егер, 1964) и для теплообмена контактирующих блоков однородной литосферы с различными начальными температурами в областях X > 0 и X < 0 (Теркот, Шуберт, 1985), а также сравнением решений, полученных для разных шагов ¿Л, Ах, Ах.

Одновременно с распределением температур рассчитывался рельеф поверхности литосферы, Н^х.гДБсЫег е1 а!., 1981):

[МХГДгДЬ/Кх.г,!)]^

Н(хД) И-2-] (34)

/>с(ге!е10 + ре(тоИо) - рт(тоИо) - рщ

и аномалии гравитационного поля в свободном воздухе, Как и в

процедуре моделирования бассейнов, учитывалось изменение плотности пород литосферы с температурой и давлением в соответствие с формулой (16), где, как и выше, параметр р0 являлся функцией трех переменных: х, г, I, меняясь при смещении границ фазовых переходов в литосфере, и растяжении литосферы. При расчетах рельефа поверхности и гравитационных аномалий предполагалось состояние изостатического равновесия литосферы с уровнем изостазии на нижней границе области счета (гм= 100-200 км). Расчет гравитационных аномалий для глубинных разрезов литосферы бассейнов проводился с помощью стандартного алгоритма вычисления гравитационного поля от четырехугольной пирамиды произвольного сечения и позволял находить аномалии гравитационного поля от произвольного распределения плотностей в слоисто-блоковых средах. К преимуществам нашго программного блока относились возможность вычисления аномалий с использованием той же конечно-разностной сетки, что применялась при решении (32), (33), возможность проводить расчеты для произвольного (симметричного или несимметричного относительно центральной оси) распределения плотностей, а также возможность вычислять вклад в А%с, от произвольного рельефа поверхности дна или суши при произвольно выбранном уровне наблюдения гравики.

Рассмотренная схема решения применялась в главе 4 диссертации дня анализа термической эволюции, рельефа поверхности и гравитационного поля литосферы осадочных бассейнов переходных зон Атлантики (Галушкин и др.,1991; 1992), пассивных окраин Антарктиды (моря Беллинсгаузена, моря Уэдцелла, Австралийского сектора континентальной окраины Антарктиды; Галушкин, Дубинин, 1990; Дубинин, Галушкин, 1990; Галушкин и др.,1991, 1992), литосферы задугового спрединга и океанической литосферы в областях перескока осей спрединга (Паресе Вела в Филиппинском море, океанической литосферы хребта Математиков и Галапагосского поднятия в восточной части Тихого океана; Галушкин и Дубинин 1991, 1992; Галушкин и др. 1991; 1992; Муравьев и Галушкин, 1991; Смирнов, Галушкин и др.,1991; Бгшгпоу, БидгоЬоу, Оа^Ыст, Л а1.,1995),

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

формами соотношений рыхлых осадков и акустического фундамента необходимо учитывать топографическую поправку совместно с поправкой на нерегулярность теплопроводности и генерации радиоактивного тепла в разнородных блоках океанической и континентальной коры. В немногих случаях простейших геометрических форм эту задачу можно решить аналитическими, или полуаналитическими методами (Lachenbruch, 1968; Череменский, 1977; Любимова и др., 1973 и др.). Однако, в реальных условиях формы рельефа и залегающих тел оказываются довольно сложными. Специально для анализа подобных случаев диссертантом был разработан компьютерный пакет для решения стационарного уравнения теплопроводности:

1"<К|:Т> + |:(К4Т> + А = 0 (35)

ди ся. dz. dz.

в области с произвольной формой границ и с произвольным распределением теплопроводности К и генерации тепла радиоактивности А внутри нее, т.е. с произвольной зависимостью этих параметров от координат х, г. Как показано в диссертации, для верхнего слоя литосферы мощностью 3-10 км, которым и ограничивается влияние топографической поправки и структурного фактора, применение стационарного уравнения (35) оправдано, если отсутствуют динамические или какие-либо иные процессы, способные влиять на температурное распределение в породах с характерным временем т < 3 - 5 млн. лет. Уравнение (35) решалось при задании температуры и теплового потока на границе.

Решение уравнения (35) с произвольной резко меняющейся верхней границей области счета, отражающей реальный сильно расчлененный рельеф дна океана или поверхности суши, находилось методом конечных элементов (Зенкевич, 1975; Bathe and Wilson, 1976). Дискретизация области на конечные элементы достигалась с помощью вариационного метода с использованием функционала, математически эквивалентного дифференциальному уравнению. Таким функционалом для уравнения (35) будет

F = ii{*[K■(¿TJ'+K-^TJ'l-A.TJ-dx.dz + J^q.T-dl (36)

Диссертантом использовалась квадратичная аппроксимация функции Т(х, z) в пределах треугольного элемента, в соответствии с которой в каждом треугольнике рассматривались шесть узлов: три в вершинах и три на серединах сторон треугольника. Условия минимума функционала dF/dTk = 0 приводило к системе линейных алгебраических уравнений относительно искомых значений температур в узлах элементов, которая решалась прямым методом исключения (Гаусса). Последнее давало существенный выигрыш во времени счета по сравнению с традиционно употребляемыми итерационными методами. Использовалась одномерная компактная форма записи матрицы жесткости, заметно экономившая оперативную память компьютера (Bathe and Wilson, 1976). Тепловой поток в узлах конечно-элементной сетки определялся по значениям найденного поля температур методом Фишера (Fisher, 1976). Точность используемого метода решения проверялась сравнением с известными аналитическими решениями и сопоставлением решений с различными разбиениями на элементы.

Метод конечных элементов использовался в диссертации при анализе теплового режима верхних слоев литосферы Командорской котловины в южной части Берингова моря (Галушкин и др., 1986), нескольких профилей через литосферу Восточной Европы, профиля через Предкарпатскую депрессию с наличием соляных диапиров (Galushkin et al.,1991), а также в ряде суббассейнов в Филиппинском море (Муравьев и Галушкин, 1991; Смирнов, Галушкин и др., 1991; Smirnov, Sugrobov, Galushkin et al.,1995).

4.2 Термическая эволюция литосферы при континентальном рифтогенезе

Обсуждение термической эволюции литосферы при континентальном рифтогенезе в разделе 4.2 диссертации начинается с анализа распределения температур в окрестности трещины, пронизывающей однородную литосферу и заполненной расплавом, полученным решением нестационарного уравнения теплопроводности методом Фурье (Ушаков, Галушкин, 1977, 1978). Это решение, несмотря на приближение, связанное с неучетом скрытой теплоты плавления, позволяет получить грубую оценку пространственных размеров тепловых аномалий. Расчеты показывают, что при заполнении узкой трещины поднимающимся горячим мантийным веществом литосфера по периферии трещины за первые миллионы лет прогревается на расстояния порядка 10 км. Подьем поверхности литосферы, обусловленный таким прогревом, при полной изостазии может составлять 2.5-3 км, из них 2-2.5 км - результат термического расширения пород и 0.4-0.7 км - следствие подъема фазовых границ в окрестности трещины в пределах литосферной мантии (Ушаков и Галушкин 1977, 1978). Решение предполагает, что с увеличением времени прогревния по трещине на ее бортах формируется некоторая конечная (стационарная) форма термического рельефа приразломных гор, с характерной шириной около 100 км.

Однако, он не согласуется с характерными линейными размерами областей пририфтового континентального поднятия большинства районов прошлого и современного рифтогенеза, которые составляют от нескольких сотен до тысячи км (Эфиопия, Аравийский полуостров). Такой рельеф не объясняется в модели прогревания по трещине и для этого привлекается модель сводового пририфтового поднятия астеносферы или модель "горячих пятен" Уилсона-Моргана. Прогревание литосферы по широкой области (600-800 км), включающей рифтовую трещину, ведет к формированию здесь обширного сводового поднятия, характеризующегося повышенным тепловым потоком на поверхности (65-85 мВт/м2) и сокращенной мощностью литосферы (50-55 км; Ушаков и Галушкин, 1978; Артюшков, 1993; Wilson, 1993; Никишин и др.,1996)). Наличие сводового поднятия учитывалось в наших моделях эволюции термического режима континентальных рифтовых зон, описанных в главе 4 диссертации. Нами проводился численный анализ термической релаксации такого первоначального возмущения термического режима литосферы рифтовых зон для ЗападноСибирского континентального бассейна (Соколова, Галушкин и др. 1990), а также ряда бассейнов пассивных окраин (Пелотос, Сантос в Южно-Американском секторе атлантической окраины, трога Фильхнера, моря Уэддела и Австралийского сектора в переходной зоне Антарктиды: Галушкин и Дубинин 1990, Галушкин и др. 1991, 1992). Анализ показал, что релаксация тепловой аномалии, сопровождающей континентальный рифтогенез, может вносить заметный вклад в формирование теплового режима литосферы еще в течении 3050 млн. лет по завершении процесса рифтогенеза.

4.3 Моделирование эволюции теплового режима и условий созревания OB в бассейнах пассивных континентальных окраин

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

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

Термическая эволюция литосферы континентальных окраин анализировалась нами на примере формирования рифтогенной континентальной окраины атлантического типа и сопряженной океанической литосферы на участке от 115° до 138° в д. австралийского сектора Антарктиды. Предрифтовое прогревание континентальной литосферы по широкой области, латеральный обмен тепла между блоками с различным глубинным строением и тепловой историей, резкие изменения скорости спрединга океанической литосферы во времени (Veevers, 1986) - таковы основные факторы, формировавшие тепловой режим литосферы этой окраины на разных этапах ее развития. Характерные черты эволюции термического режима литосферы этой континентальной окраины, рассчитанные в нашей модели, представлены на рис. 4. Как следует из расчетов, даже по прошествии 43 млн. лет после раскола континента в рельефе поверхности и глубинах изотерм континентальной литосферы еще заметны платообразные участки - следы остывания аномально прогретой континентальной литосферы широкого околорифтового поднятия (рис.4а). Они исчезают при дальнейшем остывании литосферы и уже не заметны для времени 38 млн. лет назад (Галушкин и Дубинин,1990; Галушкин и др., 1991; 1992). В процессе спрединга зона перехода от термического режима океанической литосферы к режиму континентальной становилась шире и в настоящее время имеет характерные размеры 700 - 1000 км (рис.4б). Прогревающее влияние океанической литосферы на континентальный край сказывается в возростании теплового потока, термического рельефа поверхности и глубин изотерм с приближением к зоне контакта океан-континент со стороны континента. В свою очередь, охлаждающее влияние континентального блока литосферы на океанический выражается в резком погружении изотерм и подошвы океанической литосферы в сторону континента. Тепловой поток через поверхность океанической литосферы быстро падает от значений 500 мВт/мг в осевой зоне до 50-55 мВт/м2 для участков эоценового возраста (пунктирная линия в верхней части рис.4б).

Мы использовали данные глубоководного бурения для скважин DSDP 265 (возраст фундамента 12 млн лет, мощность осадочного чехла 440 м), 266 (возраст фундамента 22.5 млн. лет, мощность осадков 370 м) и 269 (возраст фундамента 65 млн. лет, мощность осадочного покрова 6900 м), а также данные сейсмического зондирования в районе скв. DSDP 269, чтобы сопоставить вычисленные нами амплитуды тектонического погружения поверхности литосферы с наблюдаемым положением поверхности фундамента в предположении локально-изостатического отклика литосферы на поверхностную нагрузку. Считалось, что наличие продольных разломов вблизи зоны контакта океанической и континентальной литосфер делает механически независимыми вертикальные подвижки континентального и океанического блоков литосферы и тогда "нетермическая" часть влияния континентального края литосферы на погружение океанического, сводилась лишь к поставке осадочного материала с континентального склона на поверхность океанической литосферы. Сравнение показало хорошее совпадение вычисленного в нашей модели и наблюдаемого (по сейсмическим данным) рельефов поверхности фундамента (Галушкин и Дубинин, 1990; Галушкин и др., 1991; 1992). Те же расчеты демонстрируют, что в рассматриваемом случае согласие аналитических и полуэмпирических (Сорохтин, 1973; Parson and Sclater, 1977) оценок погружения океанического фундамента наблюдается лишь до возраста литосферы 30-35 млн. лет. Далее контакт с относительно холодным

в) г)

Рис 4 Термическое состояние литосферы Австрало-Антарктического хребта а) - 53, б) - 38, в) - 0 млн. лет назад, т е. после 39, 53, и 92 млн. лет с начала спрединга. В верхних частях рисунков: пунктир-тепловой поток, сплошная линия-термический рельеф по верхности литосферы. В нижних частях рисунков: сплошные линии-изотермы, и Х¥2 - границы фазовых переходов, коротко-пунктирные линии-граниа МОХО в континентальной литосфере, 21. -оцениваемое положение подошвы литосферы. г) - Реконструкция температурно-временной истории развития осадочной толщи бассейна в районе скв 269 Ввверху-тепловой поток, сплошные линии-подошвы осадочных слоев, пунктирные-изотермы; заштрихованы области с ТВИ=7-160 ("окно генерации нефти" и ТВИ > 1600 ("окно генерации сухого газа")

континентальным блоком вызывает расхождение этих оценок, которое для возраста литосферы 50-65 млн лет составляет уже 1 - 1.5 км (см. точечно-пунктирную линию и линию на верхней части рис.4б).

Реконструкция палеотемпературных условий созревания ОВ, проведенная для полного поперечного сейсмического профиля бассейна Пелотес и одного из разрезов бассейна Сантос Южно-Американской континентальной окраины, показало, что максимально достижимые степени катагенеза ОВ осадочных пород в них сильно варнруют в зависимости от положения разрезов относительно зоны контакта континентальной и океанической литосфер (Галушкин и Дубинин, 1990; Галушкин и др., 1991; 1992).

4.4 Численный анализ термического режима литосферы и условий созревания ОВ в пассивных окраинах сложного генезиса

Соввместный анализ аномалий гравитационного поля, рельефа дна и эволюции термического состояния литосферы переходной эписубдукционной окраины Тихоокеанского сектора Антарктиды в районе моря Беллинсгаузена (раздел 4.4 диссертации; Дубинин, Галушкин, 1990; Галушкин и др.,19916 1992а) дает возможность оценить температурный режим и условия созревания ОВ в осадочной толще бассейнов со сложной историей развития, включавшей перескок оси спрединга, одновременное действие двух центров спрединга и отмирание старого центра по мере его приближения к континентальной окраине в процессе субдукции океанической литосферы. Термическое состояние литосферы, изменение термического рельефа ее поверхности и тепловой поток, получались решением уравнения теплопроводности (38) по методике, описанной в диссертации и в работах (Дубинин и Галушкин, 1990; Галушкин и др., 1991). Наши расчеты показали, что толщина литосферы в осевой зоне хребта Алук за 20 млн. лет остывания увеличилась с 6 до 42 км. Тепловой поток за то же время упал здесь со значений 500 до 100 мВт/м2. Вычисленная глубина дна океана увеличивалась с удалением от оси древнего хребта Алук и достигала максимума (около 2 км относительно уровня гребня хребта Алук) на расстоянии 1300 км от оси. С последующим удалением от оси глубина дна уменьшалась, что было связано с тепловым влиянием молодой литосферы ВТП. Такая картина хорошо согласуется с характером изменения современных глубин дна моря Беллинсгаузена.

Сопоставление рассчитанных гравитационных аномалий Д^ с наблюдаемыми говорит об изостатическом выравнивании структур палеожелоба за время, прошедшее с момента прекращения пододвигания, что подтверждается и анализом тектонического рель&фа поверхности фундамента, приведенным в диссертации и работах (Дубинин и Галушкин, 1990; Галушкин и др.,1991). Оценки степени катагенеза ОВ показали, что из-за малых мощностей осадочных толщ в рассмотренных разрезах бассейна ОВ осадков не смогло достигнуть стадии катагенеза, соответствовавшей началу генерации жидких УВ, за исключением, быть может, районов, прилегающих к палеожелобу, где мощности осадочной толщи могут превосходить 2 км.

4.5 Моделирование термического режима литосферы бассейнов краевых морей

Как уже отмечалось, измеренные значения теплового потока служат одним из основных источников информации о строении и истории формирования бассейнов задуговых морей наряду с сейсмическими, магнитными и гравитационными данными (Сугробов и др., 1983; Смирнов и др., 1991). Поэтому в первой части этого раздела диссертации рассматривается существенный для

геотермии краевых бассейнов вопрос об определении поправок к наблюдаемым значениям теплового потока с целью получения корректных фоновых значений потока, характеризующих тепловое и геодинамическое состояние коры и мантии литосферы бассейна. Среди этих поправок коррекция на скорость осадконакопления, рассмотренная в разделе 2.3.9 диссертации, и на контрасты рельефа и теплопроводности, обсуждавшаяся в разделе 4.1, являются самыми распространенными. Как отмечается в диссертации, осадконакопление может занижать среднюю величину теплового потока в задуговых бассейнах на 18% -30% (Смирнов, Сугробов и Галушкин., 1982; Галушкин и др., 1986; Галушкин и Смирнов, 1987), и даже более чем на 60%, как в современном центре задугового спрединга впадины Манус в Новогвинейском море с интенсивным современным осадконакоплением (Сугробов и др., 1983).

Но более существенные искажения теплового потока, часто превышающие поправки на осадконакопление, могут возникать в зонах резко расчлененного рельефа и сложных геометрических форм соотношений рыхлых осадков и акустического фундамента, имеющих разную теплопроводность. Такие поправки к потоку определялись нами для литосферы южной части Командорской котловины с помощью решения стационарного уравнения теплопроводности (41) методом конечных элементов (42). Этим методом оценивались искажения теплового режима от неоднородностей рельефа до глубин, в 3 - 7 раз превышающих амплитуду рельефа, тогда как эволюция режима в более глубоких слоях определялась из решения нестационарного двумерного уравнения (38). Расчеты показали, что поправки за рефракцию теплового потока вдоль границ выступа поднятия фундамента южной части Командорской плиты достигают 30 -40% и остаются заметными на расстояниях до 2 км от подножья уступов, а тешГопой поток в осевой части крупной грабенообразной структуры, разделяющей Командорскую плиту и Алеутскую гряду, за счет рефракции оказывается в 1.5 раза выше, чем на бортах грабена, и у подножий склонов этого грабена возможно локальное завышение фонового теплового потока в 4 раза (Галушкин и др., 1986). В диссертации отмечается также, что в областях со сложным рельефом дна океана (на бортах грабена, склонах уступа и краевых частях поднятий фундамента) могут возникать значительные горизонтальные составляющие теплового потока (до 20 - 40% от полной величины потока).

В осевых зонах молодых центров спрединга с ограниченным осадконакоплением большой разброс в значениях теплового потока объясняется гидротермальной деятельностью (Ymano and Uyeda, 1988; Смирнов и др., 1991). В дисертации отмечается, что в целом распределение фонового теплового потока в задуговых бассейнах в зависимости от возраста коры после введения указанных поправок согласуется с распределением, наблюдаемым в срединно-океанических хребтах (см. раздел 4.5.1 диссертации; Ymano and Uyeda, 1988; Смирнов и др., 1991). Анализ геофизической информации предполагает, что задуговой спрединг, часто сопровождающийся перескоком осей растяжения, является основным механизмом формирования большинства окраинных бассейнов (Karig, 1971; Uyeda, 1977; Смирнов и др., 1988).

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

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

Еще один пример численного анализа термической эволюции литосферы краевых бассейнов, рассмотренный в диссертации, касался эволюции осадочных бассейнов Филиппинского моря. Один из этих бассейнов занимает Западно-Филиппинскую котловину, которая по всей вероятности сформировалась от 60 до 35 млн лет назад в процессе двухстадийного спрединга с центром спрединга, расположенным в Центрально-Филиппинском разломе, второй - бассейн котловины Паресе Вела образован в результате спрединга, начавшегося около 32 млн. лет назад на литосфере возраста около 30 млн. лет и продолжавшегося до 17 млн. лет назад, и, наконец, третий бассейн Марианского трога является современным активным задуговым бассейном (Hilde and Lee, 1984; Mrozowski and Hayes, 1979; Муравьев и др.,1988; Смирнов, Галушкин и др.,1991). Численная реконструкция термической истории литосферы бассейнов Филиппинского моря, рассмотренная в разделе 4.5.3 диссертации, неплохо обьясняет современное распределение теплового потока, особенности батиметрии и термической структуры фундамента этих краевых бассейнов и в целом хорошо согласуются с концепцией Карига об эволюции Филиппинского моря в результате нескольких этапов задутового спрединга (Смирнов, Ямано, Уеда, Галушкин и др.,1991).

ЗАКЛЮЧЕНИЕ

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

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

1) Рифтогенные бассейны могут развиваться в самых разнообразных геодинамических обстановках, связанных с особенностями взаимодействия литосферных плит н конфигурацией их границ: в региональных областях

растяжения литосферы (Восточно-Африканская рифтовая система), зонах регионального сдвига (Байкальский рифт), внутри локальных центров растяжения возникающих в пределах трансформных границ плит (желоб Кайман), континентальных поясов сжатия (Рейнский грабен) и задуговых областей конвергентных границ океана (Командорская котловина).

2) Процесс утонения литосферы и эволюции сс термического режима на рифтовой стадии развития бассейна может быть численно воспроизведен с соблюдением пространственных и временных масштабов в рамках модели растяжения литосферной плиты с конечным пределом пластичности. Результаты этого моделирования, а также численный анализ вариаций амплитуд тектонического погружения фундамента, проведенный для ряда осадочных бассейнов, предполагает, что продолжительность тепловой активизации литосферы на рифтогенном этапе эволюции бассейна (и в периоды его реактивизации) заметно превосходит оцениваемую в модели мгновенного растяжения литосферы.

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

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

5) Численное моделирование показывает, что перераспределение жидких и газообразных продуктов реакций созревания ОВ в процессе вторичного крекинга жидких УВ может играть заметную роль в истории реализации нефтегазогенерационного потенциала осадочных свит бассейнов при попадании осадков в область влияния высоких температур (Т > 150°С), что может иметь место ка рифтовом этапе развития бассейнов или в периоды их тепловой реактивизации.

6) Учет геологического этапа погружения образца материнской породы в процедуре численного восстановления кинетического спектра реакций из данных открытого экспериментального пиролиза образцов материнских пород позволяет уточнить низкоэнергетичную часть эффективного спектра кинетических реакций, управляющих процессом созревания ОВ, и допускает восстановление спехтра по данным пиролиза образцов с ненулевой исходной зрелостью ОВ (Яо 5 0.60-0.70%). Использование алгоритма с различными частотными факторами А! при поиске эффективного спектра реакций существенно улучшает подбор кинетических параметров реакций по сравнению с традиционным методом, предполагающим одно значение частотного фактора для всех реакций.

7) При моделировании бассейнов высоких широт Северного и Южного полушарий анализ температурного режима осадочной толщи в условиях неоднократного формирования и деградации криолитозон является обязательной частью пакета программ по моделированию бассейнов, так как современный профиль температур используется как один из контролирующих факторов. Так, вариации климата в течении последних 3.4 млн. лет в Западно-Сибирском осадочном бассейне привели к понижению температур в верхних 1500 м осадочного разреза на 15-20°С и на 8-10°С в его нижних горизонтах.

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

9) Моделирование процесса эрозии предъявляет жесткие требования к выбору пространственно-временных шагов в конечно-разностной схеме решения уравнения теплопроводности. Использование грубых шагов V/. и VI может приводить к заметным искажениям в оценках теплового потока. Термическое воздействие осадочных толщ, удаленных эрозией, на степень катагенеза ОВ нижележащих осадков кроме амплитуды эрозии определяется в сильной степени до- и постэрозионной историей погружения бассейна. Продолжительная и медленная эрозия уплотняющихся осадков может приводить, вопреки устоявшемуся мнению, и к уменьшение градиента температур в осадочной толще бассейна.

10) Моделирование теплового воздействия интрузий показало, что их внедрение в осадочный чехол бассейнов может существенно повысить уровень катаг енеза ОВ вмещающих пород на расстояниях от 0.5 до 1.0 ширины интрузии. При этом заметные отклонения значений отражательной способности витринита, 11„%, измеренных в окрестности интрузивных тел, от вычисленных в рамках модели мгновенного внедрения интрузий, находят объяснение в модели с конечным временем формирования интрузивного тела (от минут до полугода), не предполагающей прямого контакта расплавленных и вмещающих пород.

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

12) Численное моделирование бассейнов краевых морей, океана и некоторых пассивных окраин сложного генезиса позволило установить существенную роль процессов перескока осей спрединга (в отдельных случаях и одновременного действия двух центров спрединга), а также тепловой релаксации

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

Полученные результаты использовались при составлении рекомендаций практического, методического и теоретического характера в отчетах МГУ и ВНИИГеосистем, переданы» для реализации в ПГО"Севморгеология", в НПО "Союзморгео" и геологическим организациям "Архангельскгеология", "Урснгойгеология" и др.

Основные положения диссертации опубликованы в следующих работах: Монографии:

1. Кинематика плит и океаническая литосфера. Физика Земли т.З. Литосфера Земли (по геофизическим данным). ч.1. М., ВИНИТИ. 1978. 272 с (соавтор С.А. Ушаков)

2. Гравитационное поле и рельеф дна Мирового океана. Л., Недра, 1979. 296 с (соавторы С.А. Ушаков, Е.П. Дубинин, Г.И. Гапоненко, О.П. Иванов, С.С. Иванов, K.M. Каверзнев, В.Н. Шимараев)

3. Конвергентные границы плит. Физика Земли т.4. Литосфера Земли (по геофизическим данным). ч.2. М., ВИНИТИ. 1979. 224 с (соавтор С.А. Ушаков)

4. Геофизический анализ палеотектоники литосферы Земли. Физика Земли т. 7. Литосфера Земли. ч.З. М., ВИНИТИ, 1983. 228 с (соавтор С.А. Ушаков)

5. Строение и развитие литосферы переходных зон Южного океана. Физика Земли" т. 11, М., ВИНИТИ, 1991,187 с (соавторы Е.П. Дубинин, Ю.И. Прозоров, С.А. Ушаков)

6. Геотраверс Северо-Китайская равнина-Филиппинский желоб М., Наука. 1991, 150 с (соавторы Я.Б. Смирнов, М. Ямано, С. Уеда, A.B. Муравьев)

Статьи:

7. Природа складчатости на дне Черного моря в зоне перехода к Крыму и Кавказу // ДАН СССР, 1977, т.233, N 5. (соавторы С.А. Ушаков, О.П. Иванов)

8. Геодинамический анализ эволюции литосферы на раздвигающихся границах в связи с развитием Байкальской рифтовой зоны. II В сб "Континентальный рифтогенез", Новосибирск, Наука. 1977, (соавтор С.А. Ушаков)

9. Глобальная картина мгновенной кинематики литосферных плит. II Вест.Москов ун- та. Сер.геолог. 1978. N4, с.20-33. (соавтор С.А. Ушаков)

Ю.Модель мгновенной кинематики литосферных плит Земли. IIВ сб.Теодинамика дна океана", N.2, Москва, Наука. 1979 (соавтор С.А. Ушаков)

11.Мгновенная кинематика относительных перемещений литосферных плит. // В сб. "Строение литосферных плит", М., Ин-т Океанологии. 1979 (соавтор С.А. Ушаков)

12.Геометрия перемещений литосферных плит на сфере. II В сб. "Геофизика океана", т.2. Геодинамика, М.. Наука, 1979 (соавтор С.А. Ушаков)

13.Нарушение изостазии и геодинамика Карибского региона. II В сб.'Теодинамика Карибского региона", М., Наука. 1980. с. 63-77 (соавтор С.А. Ушаков, Е.П. Дубинин, А.И. Авдеев )

14.Тепловой поток в зоне сочленения Алеутско-Камчатской островодужной системы. //Вулканология и сесмология, 1982. N 6, с.96-115. (соавтор Я.Б. Смирнов, В.М. Сугробов)

15.Теплопередача при обтекании скважины фильтрационным потоком. // Всстник Моск.ун-та, Сер.геол. 1983. N6, с.94-96. (соавтор A.A. Куваев)

16.Черное море. // В сб."Методические и экспериментальные основы геотермии" М., Наука. 1983 (соавтор Я.Б. Смирнов)

^.Характеристика термической истории Восточно-Эльбской части СевероЗападноевропейского бассейна и Днспрово-Донецкой впадины: механизм развития и оценка положения зон нефтегазоносности. IIВ сб.'Тектоника плит и палеозойские ископаемые", М., 1985, с.132-147. (соавторы Я.Б. Смирнов, А.Бертольд, Г. Ольсзак)

18.Исследование структуры геотермического поля литосферы южной части Командорской котловины II Вулканология и сейсмология, 1986, N 5, с. 3-16 (соавторы A.B. Муравьев, Я.Б. Смирнов, В.М. Сугробов)

19.Методы изучения термической структуры литосферы в связи с анализом глубинного строения и оценкой перспектив нефтегазоносности (на примере эталонного сейсмологического разреза: Полярный Урал - Черное море) II В сб.'Теотермическая карта Северной Евразии и методы анализа термической структуры литосферы ", 1986. М., с. 94-157 (соавтор Я. Б. Смирнов)

20.Термическая история осадочных бассейнов; экспресс-методы оценки теплового потока .11 Гелогия и Геофизика, 1987. N 11, с.105-112 (соавтор Я. Б. Смирнов)

21.Температура и положение зон генерации углеводородов в процессе развития осадочных бассейнов // "Геология, методы поисков и разведки месторождений нефти и газа" 1988. М., вып. 12, с.6-11.

22.Результаты физического и математического моделирования континентального рифтогеза II В сб. "Математические методы прогнозирования нефтегазоносности недр" 1988, М. ч.2, с.24. (соавторы А.И. Шемеида, С. А. Ушаков)

23.Температурныс условия и положение зон генерации углеводородов в процессе развития осадочных бассейнов: описание метода и программы расчета // Жизнь Земли. 1990. М. Изд. МГУ. с.102-108

24.Геотермическая модель литосферы вдоль профиля ГСЗ "Березово-Усть-Мая" в пределах Западно-Сибирской плиты // Геология и геофизика, 1990. N 9, с. 84-92 (соавторы JI.C. Соколова, А.Д. Дучков, JI.B. Смирнов)

25.Эволюция температурного поля осадочного чехла северной части ЗападноСибирской плиты // Геология и геофизика, 1990. N 10, с. 51-60 (соавторы А.Д. Дучков, Л.В.Смирнов, Л.С. Соколова

26.Исгория развития и термическая эволюция литосферы моря Беллинсгаузена // Геофизический журнал, 1990, т. 12, N 14, с.64-70 (соавтор Е.П. Дубинин)

27.Термическая эволюция литосферы и осадочной толщи бассейнов окраины Антарктиды (Австрийский сектор) II Океанология, 1990, т.30, N I, с. 86-92 (соавтор Е.П. Дубинин)

28.Эволюция рельефа дна и термического режима литосферы при перескоке от спрединга II Тихоокеанская геология, 1991. N 6, с. 123-138 (соавтор Е.П. Дубинин)

29.Модель термической эволюции литосферы котловины Паресе Вела (Филиппинское море) // В сб."Теоретическис и экспериментальные исследования геотермического режима акваторий", 1991. М., Наука, с.157-167 (соавтор A.B. Муравьев)

30.Термический режим литосферы при перескоке оси спрединга хребта Математиков // Известия РАН, "Физика Земли", 1992, N9,1992, с.59-69 (соавтор Е.П. Дубинин)

31.Эволюция литосферы переходных окраин Антарктиды // Бюл.Москв.о-ва испыт. природы, отдел геол., 1992, т. 67, вып.З, с.3-15 (соавторы Е.П. Дубинин, С.А.Ушаков)

32.Вариации глубин зон нефтегазогенерации в пределах бассейнов континентальных окраин: пример бассейна Пелотес // Жизнь Земли. 1992. М. Изд. МГУ, с.76-85 (соавторы Е.П. Дубинин, С.А. Ушаков)

33.Модель образования и развития магматической камеры рифтовых зон срединно- океанических хребтов // ДАН РАН, 1993, т.332, N 4, с.496-499 (соавтор Е.П. Дубинин)

34.Принципы численного анализа термической эволюции литосферы и осадочных бассейнов пассивных континентальных окраин различных геодинамических типов // Геология, 1993, ч. 1, изд. МГУ, с. 252-262 (соавторы Е.П. Дубинин, С.А. Ушаков)

ЗЗ.Термическая структура осевой зоны срединно-океанического хребта. Статья 1. Формирование и эволюция осевой магматической камеры II Изв. АН РАН. сер. Физика Земли, 1994, N 5, с.11-19 (соавторы Е.П. Дубинин, А.И. Шеменда)

36.Термическая структура осевой зоны срединно-океанического хребта. Статья 2. Влияние линзы расплава на форму и эволюцию магматической камеры // Изв. АН РАН. сер.Физика Земли, 1994, N 5, с.19-26 (соавторы Е.П. Дубинин, А.И. Шеменда)

37.Магматическая камера рифтовых зон срединно- океанических хребтов: термическая модель формирования и эволюции // Вулканология и сейсмология, 1994, N 4-5, с.90-98 (соавтор Е.П. Дубинин)

38.Численное моделирование термической эволюции осадочных бассейнов и история генерации углеводородов // В сб. "Нефтегазовые ресурсы", 1994, М. вып.!, с.45-51 (соавторы Е.П. Дубинин, С.А. Ушаков)

39.Дцепрово-Донецкий палеорифт: эволюция теплового режима и нефтегазоносность // Геофизический журнал, 1995, т.17, N 3, с.13-23 (соавтор Р. И. Кутас)

40.Термический режим нефтематсринских свит Штокмановского месторождения и генерация углеводородов керогеном континентального и морского типов // В сб. "Нефтегазовые ресурсы", 1995, М., вып.2, с. 41-46 (соавторы Н.В. Лопатин, С.А. Ушаков, Г.Е. Яковлев)

41. Численное моделирование эволюции катагенеза отложений юры и триаса. // В сб. "Тюменская сверхглубокая скважина" В.Б. Мазур (ред.), 1996, Пермь, с. 279286 (соавторы Н.В.Лопатин, Т.П. Емец).

42.Криолитозона в плейстоцен-голоценовой истории Уренгойского района Западной Сибири // ДАН РАН, 1997, т. 356, п. 3, с. 393-397 (соавтор Н.В. Лопатин)

43.Геохимические предпосылки поисков нефти и газа в глубокозалегающих юрских и триасовых отложениях Западной Сибири // Геология, геофизика и разработка нефтяных месторождении, 1997, N 4, 2-15. (соавторы Н.В. Лопатин, Т.П. Емец, О.И. Симоненкова)

44.06 источнике нефтен, обнаруженных в коре выветривания и кровле палеозойского фундамента на площадях среднего Прнобья // Геология, геофизика и разработка нефтяных месторождений, 1997, N 7, 6-22. (соавторы Н.В. Лопатин, Т.П. Емец, О.И. Симоненкова)

45,Оценка влияния процессов формиирования гигантских скоплений газа на термический режим осадочных бассейнов на примере Уренгойского месторождения Западной Сибири. Геохимия, 1998 (в печати; соавторы: Симоненкова О.И., Лопатин Н.В.)

46. Термический ореол интрузии и возможный механизм ее вмещения в осадочный разрез поднятия Кэйп-Верде. Геохимия, 1998 (в печати)

47. Mathematische Modelirungen der Senkenbildung am Beispiel der N-P Senke // Zeitschriftfur angewandte Geologie, 1986. Bd.32, IIf.10, s.262-267,(coauthor A. Berthold)

48.Geothermische Modellierungen am Bcispicl der Pripjat-Senke II Zcitschrift fur an gewandte Geologic, 1986, Bd.32, Hf.ll, s.283-286, (coauthors A. Berthold, A. Muller)

49.ModeIlierungen zum Temperaturfeld wahrend des Sedimentationsprozesscs // In "Methoden und Ergebnisse geothermischer Untersuchungen", Leipzig, 1988, s. 113122 (coauthor A. Berthold)

50.Heat flow an analysis of the thermal structure of the lithosphere in the Europenian part of the USSR // In : "Exploration of the Deep Continental Crust. Terrestrial Heat Flow and the Lithosphere Structure" 1991, Springer Verlag, p.206-238 (coauthors R.I. Kutas, Ya.B. Smirniv)

51.Terrestrial Heat Flow in the Transition Zone from Asia to the NW Pacific Ocean. In Gupta M.L. and Yamano M. (eds.), "Terrestrial heat flow and geothermal energy in Asia", Rotterdam, 1995, pp. 237-250 (coathors Smirnov Ya.B., V.M.Sugrobov, A.G.Rodnikov,A.V.Muraliev, N.I.Seliverstov, V.V.Soinov and F.A.Yanovsky).

52.Evolution of sedimentary basins and petroleum formation // In "Geodynamic evolution of sedimentary basins", 1996, ed. TECHNIP, 435-453 (coauthors N.V. Lopatin, M. Makhous)

53.Burial history and kinetic modelling for hydrocarbon generation. Part I: The GALO Model II AAPG Bull. 1997, v.81, n.10, p. 1660-1678 (coauthors M. Makhous, N.V. Lopatin)

54.Burial history and kinetic modelling for hydrocarbon generation. Part II: Application of the Model to Saharan Basins // AAPG Bull. 1997, v.81, n.10, p. 1679-1699 (coauthors M. Makhous, N.V. Lopatin)

55.Numerical simulation of permafrost evolution as a part of basin modeling: permafrost in Pliocene-Holocene climate history of Urengoy field in West Siberian basin. Canad. J. Earth Science, 1997 , v.34, n.7,935-948

56.The thermal effect of igneous intrusive bodies on maturity of organic matter - A possible mechanism of intrusion formation // Organic Geochemistry, 1997, v. 27, n. 11-12, p. 645-658.

Типография ордена "Знак Почета" издательства МГУ 119899, Москва, Ленинские горы Заказ № )206 • Тираж 100 экз.

Содержание диссертации, доктора технических наук, Галушкин, Юрий Иванович

Введение.

Глава 1. ГЕОДИНАМИЧЕСКИЕ ПРЕДПОСЫЛКИ ФОРМИРОВАНИЯ РИФТОГЕННЫХ ОСАДОЧНЫХ БАССЕЙНОВ И СПЕЦИФИКА МОДЕЛИРОВАНИЯ ИХ ТЕРМО-МЕХАНИЧЕСКОЙ ИСТОРИИ.

1.1 Эволюционный ряд рифтогенных бассейнов и разнообразие reo динамических обстановок их формирования.

1.2 Термо-механические модели образования рифтовых бассейнов.

1.3 Кинематика литосферных плит и специфические типы осадочных бассейнов

Введение Диссертация по геологии, на тему "Математическое моделирование термической эволюции осадочных бассейнов и условий реализации их углеводородного потенциала"

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

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

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

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

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

1. Расчет эволюции температурного режима и полей деформаций литосферы в процессе континентального рифтогенеза.

2. Численный анализ изменения температурного режима континентальной и океанической литосферы в бассейнах с разнообразной тектонической историей и нестационарным горизонтальным полем скоростей перемещения вещества литосферы.

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

4. Оценка изменения степени метаморфизма органического вещества, реализации нефтегазогенерационного потенциала бассейна и порога первичной эмиграции жидких УВ.

5. Восстановление спектра химико-кинетических реакций, управляющих процессами созревания OB осадков.

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

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

Заключение Диссертация по теме "Геофизические методы поисков и разведки месторождений полезных ископаемых", Галушкин, Юрий Иванович

ЗАКЛЮЧЕНИЕ

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

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

1) Анализ природы образования рифтогенных бассейнов и кинематики относительных движений литосферных плит показывает, что формирование рифтовых бассейнов часто связано с особенностями взаимодействия литосферных плит, конфигурации границ и относительных движений плит на их границах. Рифтогенные бассейны могут развиваться в самых разнообразных геодинамических обстановках: региональных областях растяжения литосферы (Восточно-Африканская рифтовая система), зонах регионального сдвига (Байкальский рифт), внутри локальных центров растяжения возникающих в пределах трансформных границ плит (желоб Кайман), континентальных поясов сжатия (Рейнский грабен) или задуговых областей конвергентных границ океана (Командорская котловина).

2) Процесс утонения литосферы и эволюции ее термического режима на рифтовой стадии развития бассейна может быть численно воспроизведен с соблюдением пространственных и временных масштабов в рамках модели растяжения литосферной плиты с конечным пределом пластичности. Результаты этого моделирования, а также численный анализ вариаций амплитуд тектонического погружения фундамента, проведенный для ряда осадочных бассейнов, предполагает, что продолжительность тепловой активизации литосферы на рифтогенном этапе эволюции бассейна (и в периоды его реактивизации) заметно превосходит оцениваемую в модели мгновенного растяжения литосферы.

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

4) Вычисление отражательной способности витринита, Яо%, по кинетической модели созревания витринита и сравнение вычисленных значений с измеренными, анализ вариаций амплитуды тектонического погружения фундамента в истории бассейна и сопоставление измеренных и вычисленных глубинных распределений температур служат основными инструментами для построения и корректировки исходных параметров моделирования бассейна, наряду с имеющейся геолого-геофизической информацией о его строении и эволюции. При этом анализ вариаций амплитуды тектонического погружения поверхности фундамента может использоваться для уточнения последовательности тектоно-термических событий, имевших место в истории рассматриваемого осадочного бассейна.

5) Численное моделирование показывает, что перераспределение жидких и газообразных продуктов реакций созревания ОВ в процессе вторичного крекинга жидких УВ может играть заметную роль в истории реализации нефтегазогенерационного потенциала осадочных свит бассейнов при попадании осадков в область влияния высоких температур (Т > 150°С), что может иметь место на рифтовом этапе развития бассейнов или в периоды их тепловой реактивизации.

6) Учет геологического этапа погружения образца материнской породы в процедуре численного восстановления кинетического спектра реакций из данных открытого экспериментального пиролиза образцов материнских пород позволяет уточнить низкоэнергетичную часть эффективного спектра кинетических реакций, управляющих процессом созревания ОВ, и допускает восстановление спектра по данным пиролиза образцов с ненулевой исходной зрелостью ОВ (Яо < 0.60-0.70%). Использование алгоритма с различными частотными факторами А1 при поиске эффективного спектра реакций существенно улучшает подбор кинетических параметров реакций по сравнению с традиционным методом, предполагающим одно значение частотного фактора для всех реакций.

7) При моделировании бассейнов высоких широт Северного и Южного полушарий анализ температурного режима осадочной толщи в условиях неоднократного формирования и деградации криолитозон является обязательной частью пакета программ по моделированию бассейнов, так как современный профиль температур используется как один из контролирующих факторов. Так, вариации климата в течении последних 3.4 млн. лет в Западно-Сибирском осадочном бассейне привели к понижению температур в верхних 1500 м осадочного разреза на 15-20°С и на 8-10°С в его нижних горизонтах.

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

9) Моделирование процесса эрозии предъявляет жесткие требования к выбору пространственно-временных шагов в конечно-разностной схеме решения уравнения теплопроводности. Использование грубых шагов Уг и VI может приводить к заметным искажениям в оценках теплового потока. Термическое воздействие осадочных толщ, удаленных эрозией, на степень катагенеза ОВ нижележащих осадков кроме амплитуды эрозии определяется в сильной степени до- и постэрозионной историей погружения бассейна. Продолжительная и медленная эрозия уплотняющихся осадков может приводить, вопреки устоявшемуся мнению, и к уменьшение градиента температур в осадочной толще бассейна.

10) Моделирование теплового воздействия интрузий показало, что их внедрение в осадочный чехол бассейнов может существенно повысить уровень катагенеза ОВ вмещающих пород на расстояниях от 0.5 до 1.0 ширины интрузии. При этом заметные отклонения значений отражательной способности витринита, 1^%, измеренных в

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

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

12) Численное моделирование бассейнов краевых морей, океана и некоторых пассивных окраин сложного генезиса позволило установить существенную роль процессов перескока осей спрединга (в отдельных случаях и одновременного действия двух центров спрединга), а также тепловой релаксации отмерших центров спрединга в формировании термического режима литосферы этих бассейнов. Расчеты показали, что при перескоках локальных осей спрединга в задуговых бассейнах процессы латерального теплообмена во многом определяют ширину области перехода между старым и молодым блоками литосферы, а также скорость релаксации рельефа палеоспрединговых хребтов.

Библиография Диссертация по геологии, доктора технических наук, Галушкин, Юрий Иванович, Москва

1. Аплонов C.B. Локальные палеооси спрединга: магнитные аномалии и геодинамика. Изв. АН СССР. Сер.Физика. Земли, 1989, п.1, с. 3-12.

2. Артюшков Е.В. Физическая тектоника. М., Наука, 1993, 457 с.

3. Артюшков Е.В., Бэер М.А., Соболев C.B., Яншин A.A. Механизм образования геосинклинальных поясов. Советская геология, 1982, п. 1, с. 14-27.

4. Артюшков Е.В., Бэер М.А. Механизм образования глубоких прогибов на континентальной коре Урала, Аппалачей и Скандинавских каледонид. Изв. АН СССР, сер. геолог., 1983, п.5, с. 25-52.

5. Артюшков Е.В., Бэер М.А. О роли горизонтальных и вертикальных движений в образовании прогибов на континентальной коре в складчатых поясах. Изв. АН СССР, сер. геолог., 1984, п.5, с. 3-28.

6. Артюшков Е.В., Мернер H.A. Быстрое образование крутых флексур под осадочными бассейнами указание на временную потерю прочности литосферного слоя. ДАН РАН, 1997, т. 356, п. 3, с. 382-386.

7. Артюшков Е.В., Шлезингер А.Е., Яншин A.A. Основные типы и механизмы образования структур на литосферных плитах. Статья 1. Континентальные платформы. Бюлл. МОИП, 1979, т.54, вып. 2, с. 8-30.

8. Архипов С.А. Палеогеография и хроностратиграфия плейстоцена севера Сибири (обзор новейших данных)//В сб."Четвертичный период",ред.Яншин А.Л., Кишинев, 1989, с.201-214.

9. Ю.Архипов С.А. Хроностратиграфия плейстоцена-основа для палеоклиматических реконструкций и периодизации биоты // В сб. "Эволюция климата, биоты и среды обитания человека в позднем кайнозое Сибири", ред.В.АЗахаров, Новосибирск, 1991, с. 17-30.

10. Балобаев В.Т. Геотермия мерзлой зоны литосферы Севера Азии. Новосибирск, Наука, 1991, 193 с.

11. Баулин В.В. Основные этапы истории развития многолетнемерзлых пород на территории Западно-Сибирской низменности // Труды ин-та мерзлотоведения им.В. А.Обручева, 1962, т.19, с.5-18.

12. Белонин М.Д.,Григоренко Ю.Н.,Киричкова А.И. и др. Геология нефти и газа, методы прогноза и поисков месторождений углеводородов (по материалам 29-ой сессии МГК), М., Геоиннформмарк, 1994, с.15-17.

13. Бяков В.М., Шиманов Г.Г.Б Степанова О.П. Эффект ионизирующей радиации на углефикацию органической материи в Земле. Химия высоких энергий, 1987, т.21, N 1, с. 45-49.

14. Буракович Т.К., Воронцов И.М., Гордиенко В.В. и др. О природе Тарханкутской аномалии в Крыму. Геофизич. журнал, 1987, т.9, п.5, с. 50-58.

15. Величко A.A. Структура вариаций климата в мезо-кайнозойское время по данным из Восточной Европы.// В сб"Климаты в геологической истории Земли", ред. Величко A.A., Чепалыга А.Л., М, Наука, 1987, с.5-53.

16. Величко А.А.(ред.). Западная Сибирь. Развитие ландшафтов и климата Северной Евразии. Поздний плейстоцен-голоцен, М.,Наука, 1993, вып.1, с.27-46.

17. Величко A.A., O.K. Борисова, С.А.Лаухин, Т.В. Светлицкая. Реконструкция палеоклиматов плиоцена Северной Евразии. // В сб'Ъсероссийское совещание по изучению четвертичного периода", ред. А.Л. Яншин, М., 1994, с.54.

18. Верба М.Л., Алексеева А.Б. Влияние интрузии на содержание битумов во вмещающих палеозойских карбонатных породах Норильского района. Сб. "Вопросы геологии и нефтегазоносности Тунгусской синеклизы", Ленинград, Труды ВСЕГЕИ, 1972, п. 308, с. 124-142.

19. Волкова B.C. Колебания климата в Западной Сибири в позднеплиоценовое и четвертичное время. //В сб. "Эволюция климата, биоты и среды обитания человека в позднем кайнозое Сибири", ред.В.АЗахаров, Новосибирск, 1991, с.30-40.

20. Волкова B.C. Динамика природной среды и климат казанцевского межледниковья в Западной Сибири (по данным палинологии). // В сб'Ъсероссийское совещание по изучению четвертичного периода", ред. А.Л. Яншин, Тез.докл., М, 1994, с.58.

21. Вотах М.Р. и В.А. Климатов. Растительность и климат Томского Приобья в голоцене. // Геология и геофизика, 1994, т.35, N 10, с.25-31.

22. В.Н.Гавриш, Г.Д.Забелло, Л.И.Рябчун и др. Геология нефтегазоносности Днепрово-Донецкой впадины. Глубинное строение и геотектоническое развитие.// Киев: Наук.Думка, 1987 208 с.

23. Галушкин Ю.И., Ушаков С.А. Глобальная картина мгновенной кинематики литосферных плит. Вест. МГУ, сер.геолог., №4, 1978.

24. Галушкин Ю.И., Ушаков С.А. Модель мгновенной кинематики литосферных плит Земли. Сб. "Геодинамика дна океана", N.2, Москва, Наука. 1979.

25. Галушкин Ю.И., Ушаков С.А. Геометрия перемещений литосферных плит на сфере. Кн. "Геофизика океана, т.2. Геодинамика", М. Наука, 1979.

26. Галушкин Ю.И., Куваев A.A. Теплопередача при обтекании скважины фильтрационным потоком. Вестник Моск. ун-та, Сер.геол., 1983, N6, с.94-96.

27. Галушкин Ю.И., A.B. Муравьев, Я.Б. Смирнов, В.М. Сугробов Исследование структуры геотермического поля литосферы южной части Командорской котловины. Вулканология и сейсмология, 1986, N 5, с.3-16.

28. Галушкин Ю.И., Смирнов Я.Б. Термическая история осадочных бассейнов; экспресс-методы оценки теплового потока. Гелогия и Геофизика, 1987. N11, с. 105-112.

29. Галушкин Ю.И., Шеменда А.И., Ушаков С.А. Результаты физического и математического моделирования континентального рифтогеза. Сб. "Математические методы прогнозирования нефтегазоносности недр" 1988, М. ч.2, с.24.

30. Галушкин Ю.И. Температура и положение зон генерации углеводородов в процессе развития осадочных бассейнов "Геология, методы поисков и разведки месторождений нефти и газа" 1988. М., вып. 12, с. 6-11.

31. Галушкин Ю.И., Прозоров Ю.И. Термические условия в осадочной толще ЮжноКаспийской впадины в связи с вопросами нефтегазообразования. Сб. "Геодинамические методы прогнозирования нефтегазоносности недр" 1988. М. ч.2, с.225.

32. Галушкин Ю.И. Температурные условия и положение зон генерации углеводородов в процессе развития осадочных бассейнов: описание метода и программы расчета "Жизнь Земли". 1990. М. Изд. МГУ. с. 102-108

33. Галушкин Ю.И., Дубинин Е.П. Термическая эволюция литосферы и осадочной толщи бассейнов окраины Антарктиды (Австрийский сектор). Океанология, 1990,t.30,N 1,с.86-92.

34. Галушкин Ю.И., Дубинин Е.П., Прозоров Ю.И., Ушаков С.А. Строение и развитие литосферы переходных зон Южного океана. Физика Земли", т. 11, М.,ВИНИТИ, 1991,187с.

35. Галушкин Ю.И., Муравьев A.B. Модель термической эволюции литосферы котловины Паресе Вела (Филиппинское море). Сб. "Теоретические и экспериментальные исследования геотермического режима акваторий", 1991. М., Наука, с. 157-167.

36. Галушкин Ю.И., Дубинин Е.П. Эволюция рельефа дна и термического режима литосферы при перескоке от спрединга. Тихоокеанская геология, 1991. N 6, с. 123-138.

37. Галушкин Ю.И., Дубинин Е.П. Термический режим литосферы при перескоке оси спрединга хребта Математиков. Известия РАН, "Физика Земли", 1992, N9,1992, с.59-69.

38. Галушкин Ю.И.„ Дубинин Е.П., Ушаков С.А. Эволюция литосферы переходных окраин Антарктиды. Бюлл.МОИП, отд.геол., т.67, вып. 3, 1992 (13 стр.)

39. Галушкин Ю.И. Дубинин Е.П., Ушаков С.А. Вариация глубин зон нефтегазогенерации в пределах бассейнов континентальных окраин. Сб. Жизнь Земли, изд.МГУ, 1992, с. 76-85

40. Галушкин Ю.И., Дубинин Е.П. Модель образования и развития магматической камеры рифтовых зон срединно- океанических хребтов // ДАН РАН, 1993, т.332, N 4, с.496-499.

41. Галушкин Ю.И., Е.П. Дубинин, С.А. Ушаков. Принципы численного анализа термической эволюции литосферы и осадочных бассейнов пассивных континентальных окраин различных геодинамических типов. Сб." Геология" , ч.1, изд. МГУ, 1993, с. 252-262.

42. Галушкин Ю.И., Дубинин Е.П., Шеменда А.И. Термическая структура осевой зоны срединно-океанического хребта. Статья 2. Влияние линзы расплава на форму и эволюцию магматической камеры // Изв. АН РАН. сер.Физика Земли, 1994, N 5, с. 19-26.

43. Галушкин Ю.И., Дубинин Е.П. Магматическая камера рифтовых зон срединно-океанических хребтов: термическая модель формирования и эволюции // Вулканология и сейсмология, 1994, N 4-5, с.90-98.

44. Галушкин Ю.И., Дубинин Е.П., Ушаков С.А. Численное моделирование термической эволюции осадочных бассейнов и история генерации углеводородов. Сб. "Нефтегазовые ресурсы", 1994, М. вып.1, с.45-51

45. Галушкин Ю.И., Кутас Р.И. Днепрово-Донецкий палеорифт: эволюция теплового режима и нефтегазоносность. Геофизический журнал, 1995, т. 17, N 3, с. 13-23.

46. Галушкин Ю.И., Н.В. Лопатин, С.А. Ушаков, Г.Е. Яковлев. Термический режим нефтематеринских свит Штокмановского месторождения и генерация углеводородов керогеном континентального и морского типов. В сб. "Нефтегазовые ресурсы", 1995, М., вып.2, с. 41-46.

47. Галушкин Ю.И., Н.В. Лопатин, Т.П. Емец. Численное моделирование эволюции катагенеза отложений юры и триаса. Сб. "Тюменская сверхглубокая скважина" В.Б. Мазур (ред.), 1996, Пермь, с. 279-286.

48. Галушкин Ю.И., Н.В. Лопатин. Криолитозона в плейстоцен-голоценовой истории Уренгойского района Западной Сибири // ДАН РАН, 1997, т. 356, п. 3, с. 393-397.

49. Галушкин Ю.И., Симоненкова О.И., Лопатин Н.В. Оценка влияния процессов формиирования гигантских скоплений газа на термический режим осадочных бассейнов на примере Уренгойского месторождения Западной Сибири. Геохимия, 1998 (в печати).

50. Галушкин Ю.И. Термический ореол интрузии и возможный механизм ее вмещения в осадочный разрез поднятия Кэйп-Верде. Геохимия, 1998 (в печати)

51. Глико А.О. Асимптотическое решение задачи о термическом утонении литосферы. Изв.АН СССР, сер. Физика Земли, 1986, п. 4, с.3-13.

52. Гольдберг A.B. Основы региональной палеоклиматологии. М., Недра, 1987,222 с.

53. Гольмшток А.Л. Эффект осадконакопления на глубинный тепловой поток. Океанология, 1979, т. 19, N 6, с. 1133-1138.

54. Гольмшток А.Л. Тепловой поток, вызванный генерацией тепла при отложении осадков. Океанология, 1981, т.21, N 6, с.1029-1033.

55. Горальчук М.И., Крицук Л.М. Особенности распределения температуры грунтов на левобережье р.Пур (междуречье рек Ево-Лха и Ягенетта)//Геокриологические исследования, М., Труды ВСЕГИНГЕО, 1980, N 138, с.22-34.

56. Гохман М.Р. Определение расчетных температур засоленных глинистых грунтов с учетом фазовых переходов грунтовой влаги в спектре отрицательных температур. Основания, фундаменты и механика грунтов, 1993, N 1, с. 19-21.

57. Грачев А.Ф. Рифтовые зоны Земли. М., Недра, 1987.

58. Грачев А.Ф., Николайчик В.В. Термические модели пострифтового погружения. Докл. АНСССР, 1985, т. 283, с. 1439-1442.

59. Дейли Дж, Хамман Д. Механика жидкости. М., ИЛ, 1971.

60. Дубинин Е.П. Трансформные разломы океанической литосферы. М., МГУ, 1987, 182 с.

61. Дубинин Е.П. Палеограницы плит океанической литосферы. Тихоокеанская геология, 1994, п.З, с. 3-21.

62. Дубинин Е.П., Галушкин Ю.И. История развития и термическая эволюция литосферы моря Беллинсгаузена. Геофизический журнал, 1990, т. 12, N. 14, с. 64 -70.

63. Дучков А.Д., Лысак C.B., Балобаев В.Т. и др. Тепловое поле недр Сибири. Новосибирск, Наука, 1987, 196 с.

64. Дучков А.Д., Соколова Л.С., Новиков Г.Н., и др. Тепловой поток юго-восточной части Западно-Сибирской плиты. Геология и Геофизика, 1988, п. 8, с. 77-85.

65. Дучков А.Д., Галушкин Ю.И., Л.В.Смирнов, Л.С. Соколова Эволюция температурного поля осадочного чехла северной части Западно-Сибирской плиты. Геология и геофизика, 1990. N 10, с. 51-60.

66. Ершов Э.Д (ред.) Геокриология СССР: Западная Сибирь. М., Недра, 1989, 456 с.

67. Зыкин B.C., B.C. Зажигин, А.Ю. Казанский. Поздний неоген юга Западно-Сибирской равнины: стратиграфия, палеомагнетизм, основные климатические события.// Геология и Геофизика, 1991, N 1, с.78-86.

68. Истомин В.А.,В.С. Якушев. Газовые гидраты в природных условиях.М,Недра, 1992,236 с.

69. Казаринов В.В., Хоменко A.B. Влияние траппов на палеозойские нефте и газогенерирующие породы района Тунгусски и Лены. Сб "Литология и геохимия нефтегазоносных толщ Сибирской платформы", Наука, М, 1981, с. 113-117.

70. Карслоу Г., Егер Д. Теплопроводность твердых тел. М., Наука, 1964, с.487.

71. Климанов В. А. Особенности изменения климата Северной Евразии в позднеледниковье и голоцене. //Бюлл.МОИП, 1994, т.69, N 1, с.58-63.

72. Климанов В.А. и В.В. Клименко. Колебания температуры в климатических оптимумах голоцена и плейстоцена// ДАН России, 1995, т.342, N 2, с.242-245.

73. Козлов В.П., Токарев ЛВ Масштабы образования газа в осадочных толщах (пример Донецкого бассейна). Советская Геология, 1961, п.7, с. 19-33.

74. Конторович А.Э., Нестеров И.И., Салманов Ф.К. и др. 1975. Геология нефти и газа Западной Сибири. М., Недра, с.250.

75. Конторович А.Э., Фотиади Э.Э., Демин В.И. и др. 1981. Прогноз месторождений нефти и газа. М., Недра, с.350.

76. Котляков В.Н. Глобальные изменения природы в "зеркале" ледяного керна. // Природа, 1992, N 7, с.59-68.

77. Кудрявцев В. А. Мерзлотоведение. М., изд. МГУ, 1981, 240 с.

78. Курчиков А. Гидро-геотермические критерии нефтегазоносности. М., Недра, 1992, 231 с.

79. Курчиков А.Р., Ставицкий Б.П. 1987. Геотермия нефтегазоносных областей Западной Сибири. М., Недра, с. 134.

80. Кучерук Е.В., Ушаков С.А. Тектоника плит и нефтегазоносность (геофизический анализ). Физика Земли т.8, М., ВИНИТИ, 1985, 200 с.

81. Кутас Р.И. Поле тепловых потоков и термическая модель земной коры. Киев, Наукова думка, 1978, 140 с.

82. Кутас Р.И., Цвещенко В.А., Корчигин И.Н. Моделирование теплового поля континентальной литосферы. Киев, Наукова думка, 1989, 191 с.

83. Лаврушин Ю.А. и М.Н. Алексеев. Позднекайнозойские палеоклиматические события в Арктическом океане. // Стратиграф. Геолог, корреляция, 1993, т.1, N 6, с.48-52.

84. Лаухин С.А. Изменения климата в плиоцен-плейстоцене Северо-Востока Азии. // Стратиграф. Геолог, корреляция, 1993, т.1, N 6, с.59-65.

85. Лопатин Н.В. Температура и геологическое время как фактор углефикации. Известия АН СССР, сер. геолог., 1971, N3, с. 95-106.

86. Лопатин Н.В., Емец Т.П. Пиролиз в нефтегазовой геохимии. Москва, Наука, 1987, 144 с.

87. Лопатин Н.В., Емец Т.П., Симоненкова О.И., Галушкин Ю.И. Геохимические предпосылки поисков нефти и газа в глубокозалегающих юрских и триасовых отложениях Западной Сибири.// Геология, геофизика и разработка нефтяных месторождений. 1997. N 4. стр. 2-16.

88. Лопатин Н.В., Емец Т.П., Симоненкова О.И., Галушкин Ю.И. Об источнике нефтей, обнаруженных в коре выветривания и кровле палеозойского фундамента на площадях среднего Приобья // Геология, геофизика и разработка нефтяных месторождений, 1997, N 7, 6-22.

89. Любимова Е.А. Термика Земли и Луны. М, Наука, 1968, 279 с.

90. Любимова Е.А., Александров А.Л., Дучков А.Д. Методика изучения тепловых потоков через дно океанов. М. Наука. ,1973, 176 с.

91. Ю4.Милановский Е.Е. Рифтогенез в истории Земли. М., Недра, 1983, 1987.

92. Милановский Е.Е. 1991. Рифтогенез и его роль в тектоническом строении Земли и ее мезо-кайнозойской геодинамики. Геотектоника, N 1, с. 3-21.

93. Юб.Муравьев A.B. Тепловой поток в южной части Командорской впадины. В сб.'Теотермические исследования на дне акватории" М., Наука, 1988, с. 438-443.

94. Муравьев A.B., Селивестров Н.И., Смирнов Я.Б., Сугробов В.М. Новые данные по тепловому потоку в районе подводного четвертичного вулканизма Командорской котловины. Докл. АНСССР, 1990, т.312, п.2, с. 438-443.

95. Никишин А.М. О зависимостти структуры континентальных рифтов от мощности и физического состояния дорифтовой литосферы. ДАН СССР, 1987, т. 293, п. 5, с. 12091213.

96. Ю9.Пасконов В.М., В.И. Полежаев, Л.А. Чудов. Численное моделирование тепло- и массо-переноса. М., Наука, 1984, 285 с.

97. Ровенская A.C., Горшков В.И., Немченко H.H. Ораздельном прогнозировании зон газо- и нефтенакопления на Севере Западной Сибири. Тр. Зап. СибНИГНИ, 1978, п. 138, с. 83-88.

98. ПЗ.Роднова E.H. Изменение коллекторских свойств осадков в контактной зоне траппов в центральной части Тунгусской синеклизы. Сб. "Вопросы геологии и нефтегазоносности Тунгусской синеклизы", Ленинград, Труды ВСЕГЕИ, 1976, п. 308, с. 118-133.

99. Сальников В.Е. Геотермический режим Южного Урала. М., Наука, 1984, 89 с.

100. Самарский A.A., Гулин A.B. Численные методы. М., Наука, 1989, 430 с.

101. Пб.Селивестров Н.И., Баранов Б.В., Егоров Ю.О., Шкира В.А. Новые данные о строенииюжной части Командорской котловины по результатам 25 и 26-го рейса НИС "Вулканолог". Вулканология и сейсмология, 1988, п.4, с. 3-20.

102. Сигунов Ю.А., Фартышев А.И. Исследования эволюции криолитозоны Арктического шельфа методами математического моделирования. Геология и геофизика, 1991,N 8,с24-31

103. Сигунов Ю.А., Фартышев А.И. Промерзание и протаивание Восточно-Арктического шельфа в позднем плейстоцене (вычислительный эксперимент). Геология и геофизика, 1995, т. 36, N9, с.36-41.

104. Смирнов Я.Б. Тепловое поле территории СССР : пояснительная записка к картам тепллового потока и глубинных температур в массштабе 1:10 ООО ООО. Москва, ГУГК, 1980, 150 с.

105. Смирнов Я.Б., Сугробов В.М., Галушкин Ю.И. Тепловой поток в зоне сочленения Алеутско-Камчатской островодужной системы. Вулканология и сейсмология, 1982. N 6, с.96-115.

106. Смирнов, Я.Б. Галушкин Ю.И., Черное море. Сб. "Методические и экспериментальные основы геотермии" М., Наука. 1983.

107. Смирнов Я.Б., Сугробов В.М. Земной тепловой поток в Курило-Камчатской и Алеутской провинциях, часть 2. Карта измеренного и фонового теплового потока. Вулканология и сейсмология, 1980, п.1, с. 16-31.

108. Смирнов Я.Б., М. Ямано, С. Уеда, Ю.И. Галушкин, A.B. Муравьев и др. Геотраверс Северо-Китайская равнина-Филиппинский желоб-Марианский желоб. М., Наука. 1991, 150 с.

109. Соболев C.B. Модели нижней части коры на континентах с учетом фазового перехода габбро-эклогит. Сб. "Проблемы петрологии земной коры и верхней мантии", Новосибирск, 1978, с. 347-355.

110. Соколова JI.C., Ю.И. Галушкин, А.Д. Дучков, Л.В. Смирнов Геотермическая модель литосферы вдоль профиля ГСЗ "Березово-Усть-Мая" в пределах Западно-Сибирской плиты // Геология и геофизика, 1990. N 9, с. 84-92.

111. Сорохтин О.Г., Ушаков С.А. Глобальная эволюция Земли. М., изд. МГУ, 1991, 446 с.

112. Сурков B.C., Смирнов Л.В., Жеро О.Г. Раннемезозойский рифтогенез и его влияние на структуру литосферы Западно-Сибирской плиты. Геология и геофизика, 1987, N 9, с.3-11.

113. Сурков B.C., Смирнов A.B. Тектонические события кайнозоя и фазовая дифференциация углеводородов в готерив-сеноманском комплексе Западно-Сибирского бассейна.//Геология нефти и газа. 1994. N11. с.3-43.

114. Сухорукова С.С. К вопросу о палеогляциоклиматических реконструкциях плейстоценовых оледенений Западной Сибири. //В сб. "Эволюция климата, биоты и среды обитания человека в позднем кайнозое Сибири", ред.В.А.Захаров, Новосибирск, 1991, с.78-84.

115. Теркот Д., Шуберт Дж. Геодинамика. М., Мир, 1985, т.1, 374 с

116. Ушаков С.А., Красс М.С. Сила тяжести и вопросы механики недр Земли. М., Недра, 1972, 157 с.

117. Ушаков С.А. Строение и развитие Земли. Физика Землит.1. М., ВИНИТИ, 1974, 269 с.

118. Ушаков С.А, Галушкин Ю.И. Геодинамический анализ эволюции литосферы на раздвигающихся границах в связи с развитием Байкальской рифтовой зоны. Сб "Проблемы рифтогенеза", Иркутск, 1975.

119. Ушаков С.А., Галушкин Ю.И. Геодинамический анализ эволюции литосферы на раздвигающихся границах в связи с развитием Байкальской рифтовой зоны. Сб "Континентальный рифтогенез", Новосибирск, Наука. 1977.

120. Ушаков С.А., Галушкин Ю.И., Городницкий А.М. Условия срыва подводных гор при их затягивании под островные дуги. Сб. "Тектоника литосферных плит", М., ИО АН СССР, 1977 (9 стр.)

121. Ушаков С.А., Галушкин Ю.И., Иванов О.П. Природа складчатости на дне Черного моря в зоне перехода к Крыму и Кавказу. ДАН СССР, 1977, т.233, N 5.

122. Ушаков С.А., Галушкин Ю.И., Дубинин Е.П., Иванов О.П. Планетарные разломы океанической литосферы (геодинамическая характеристика) // В сб " Разломы земной коры", М, Наука, 1977, (12 стр.)

123. Ушаков С.А., Галушкин Ю.И. Кинематика плит и океаническая литосфера. Физика Земли т.З. Литосфера Земли (по геофизическим данным). ч.1. М., ВИНИТИ, 1978. 272 с.

124. Ушаков CA., Галушкин Ю.И. Мгновенная кинематика относительных перемещений литосферных плит. Сб. Строение литосферных плит, М., ИО АН СССР, 1979

125. Ушаков С.А., Ю.И. Галушкин, Е.П. Дубинин, Г.И. Гапоненко, О.П. Иванов, С.С. Иванов, K.M. Каверзнев, В.Н. Шимараев. Гравитационное поле и рельеф дна Мирового океана. JL, Недра, 1979. 296 с.)

126. Ушаков С.А., Галушкин Ю.И. Конвергентные границы плит. Физика Земли т.4. Литосфера Земли (по геофизическим данным). ч.2. М., ВИНИТИ. 1979. 224 с

127. Ушаков С.А., Авдеев А.И., Галушкин Ю.И., Дубинин Е.П. Нарушение изостазии и геодинамика Карибского региона. Сб. "Геодинамика Карбского региона", М., Наука, 1980, 63-77.

128. Ушаков С.А., Галушкин Ю.И., Дубинин Е.П., Иванов О.П. Природа подводных вулканических гор // В сб."Комплексные исследования природы океана" М., 1980. N.7, с. 146-150.

129. Ушаков С.А., Галушкин Ю.И. Деформация осадков на конвергентных границах плит ( в связи с проблемой складчатости на дне Черного и Каспийского морей) . Сб. "Проблемы геодинамики Кавказа", М., Наука, 1982 (12 стр.)

130. Ушаков С.А., Галушкин Ю.И. Геофизический анализ палеотектоники литосферы Земли. Физика Земли т. 7. Литосфера Земли. ч.З. М., ВИНИТИ, 1983. 228 с.

131. Феоктистов Г.Д. Метаморфизм глинисто-песчаных пород вблизи от контакта с интрузией. Наука, М, 1972, 100 стр.

132. Федотов С.А. Поднятие основных магм в земной коре и механизм излияния базальта по трещинам. Изв. АН СССР, сер. геол., 1976, п. 10, стр. 5-23

133. Фролов Н.М., Лялько В.И., Митник М.М. Гидротермические методы исследований в нефтяной гидрогеологии. В сб. :Геотермические методы исследований в гидрогеологии. М.:Недра, 1979. с.210-212.

134. Хаин В.Е., Божко H.A., Сеславинский К.Б., Балуховский А.Н. Историческая геотектоника. М. Недра, 1988-1993.

135. Хаин В.Е., Ломизе М,Г. Геотектоника с основами геодинамики. М., МГУ, 1995, 479 с.

136. Череменский Г. А. Прикладная геотермия. Ленинград, Недра, 1977,225 с.

137. Шеменда А.И. Критерии подобия при механическом моделировании тектонических процессов. Геология и геофизика, 1983, п. 10, с. 10-19.

138. Шеменда А.И. Некоторые закономерности деформации литосферы при растяжении (по результатам физического моделирования). Докл.АНСССР, 1984, т.275, п.2, с. 346-348.

139. Шик С.М. Климатическая ритмичность в плейстоцене Восточно-Европейской платформы. // Стратиг.Геолог.коррел., 1993, т.1, N4, с. 12-24.

140. Яновский Ф.А., Сугробов В.Н., Селивестров Н.И. Тепловое поле и геотермическая модель Командорской котловины Берингова моря. Вулканология и сейсмология, 1997, п.2, с. 16-32.

141. Alvarez F., Virieux J., and Le Pichon X. Thermal consequences of lithosphere extension over continental margins the initial stretching phase. Geophys. J. Roy. astron. Soc. , 1984, v. 78, p. 389-411.

142. Anderson, D.L., 1979, The deep structure of continents: J.geophys.Res., v.84, N.B13, p.7555-7560.

143. Anderson, D.L., 1980, The temperature profile of the upper mantle: J.geophys.Res.,v.85, N.B12, p.7003-7010.

144. Barid D. J., Knapp J.H., Steer D.H., Brown L.D., and Nelson K.D. Upper mantle reflectivity beneath the Willinston basin, phase-change Moho, and the origin of intracratonic basins. Geology, 1995, v. 23, n. 5, p. 431-434.

145. Barker C., 1988. Generation of anomalous internal pressures in source rocks and its role in driving petroleum migration. Revue de I'nstitut Franqais de Petrole, v.43, p. 349-355.

146. Bathe K-J., and Wilson E.L. Numerical methods in finite element analysis, N-Y, 1976, 470 p.

147. Beck, E., 1976, An improved method of computing the thermal conductivity of fluid-filled sedimentary rocks: Geophysics, v.41, N.l, p. 133-144.

148. Berthold A. und Galushkin Yu.I. Mathematische Modelirungen der Senkenbildung am Beispiel der N-P Senke. Zeitschrift für angewandte Geologie, 1986. Bd.32, Hf.10, s.262-267.

149. Berthold A., Galushkin Yu.I., und Muller A. Geothermische Modellierungen am Beispiel der Pripjat-Senke. Zeitschrift fur an gewandte Geologie, 1986, Bd.32, Hf.l 1, s.283-286.

150. Berthold A. und Galushkin Yu.I. Modellierungen zum Temperaturfeld wahrend des Sedimentationsprozesses. In "Methoden und Ergebnisse geothermischer Untersuchungen", Leipzig, 1988, s. 113-122.

151. Bertotti G. and ter Voorde M. Thermal effects of normal faulting during rifted basin formation. 2. The Lungano-Val Grande normal fault and the role of pre-existing thermal anomalies. Tectonophysics, 1994, v. 240, n. 1-4, p. 145-157.

152. Bethke C.M. 1985. A numerical model of compaction-driven groundwater flow and heattransfer ant its application to the paleohydrology of intracratonic sedimentary basin J.Geophys.Res., v.90, p.6817-6828.

153. Bethke C.M. 1989. Modelling subsurface flow in sedimentary basins. Geologische Rundschau, v.78,Nl, p. 129-154.

154. Beaumont Ch., Ch.E.Keen, and R.Boutilier. 1982, On the evolution of rifted continental margins: comparison of models and observations for the Nova Scotian margin., Geophys.J.Roy.Astr.Soc.,v.70, p.667-715.

155. Bishop A.N. and Abbott C.D. The interrelationship of biological marker maturity parameters and molecular yields during contact metamorphism. Geochim. Cosmochim. Acta, 1993, 87, n. 15, p. 3661-3668.

156. Bishop A.N. and Abbott C.D. Vitrinite reflectance and molecular geochemistry of Jurassic sediments: the influence of heating by Tertiary dykes (northwest Scotland). Org. Geochem., 1995, 22, p. 165-177.

157. Bradley J.S., 1975. Abnormal formation pressure. AAPGBull. v.59, p. 957-973.

158. Bredehoeft J.D., R.D. Djevanshir, and K.R. Beilitz. 1988. Lateral fluid flow in a compacting sand-shale sequence, south Caspian Sea. AAPGBull. v.12, p. 416-424.

159. Bredehoeft J.D., and I.S. Papadopulos. 1965. Rates of vertical groundwater movementestimated from the Earth's thermal profile. Water Resour. Res. v.l, N 2, p.325-328.

160. Buntebrach G. and Stegena L. Methods in paleogeothermics. in: "Lectures Notes in Earth Sciences", Bhattacharjis et al. (eds.), Springer Verlag, 1986, v.5, p. 5-43.

161. Burnham and Braun 1990 Development of a detailed model of petroleum formation, destruction, and expulsion from lacustrine and marine source rocks. Org.Geochemistry, v. 16, N.l-3, p.27-39.

162. Burnham A.K. and Sweeney J.J. A chemical kinetic model of vitrinite maturation and reflectance. Geochim. Cosmochim. Acta, 1989, v. 53, n. 10, p. 2649-2657.

163. Burruss R.C. Paleotemperatures from fluid inclusions: advances in theory and technique. In: "Thermal history of sedimentary basins" Naeser N.D. and McCulok T.H. (eds.), Springer Verlag, 1989, p. 119-131.

164. Burrus J, Andebert F. 1990, Thermal and compaction processes in a young rifted basin containing evaporites: gulf of Lions, France. AAPGBull.v 74, N 9, p 1420-1440.

165. Cande S.C., Herron E.M., and Hall B.R. Cenozoic tectonic history of the southeast Pacific. Earth Planet. Sci. Letters, 1982, v.57, p.63-74.

166. Cermak V., Bodri L., Rybach L., and Buntebarth G. Relationship between seismic velocity and heat production: comparison of two sets of data and test of validity. Earth Planet. Sci. Letter, 1990, v. 99, p.48-57.

167. Chekunov A.V., Gavrish V.K., Kutas R.I., Ryabchun L.I. Dnieper-Donets paleorift// Tectonophysics, 1992, v.208, p.257-272.

168. Class C. and Goldstein S.L. Plume-lithosphere interactions in the ocean basins: constraints from the source mineralogy. Earth Planet. Sci. Letters, 1997, v. 150, p. 245-260.

169. Clauser C., and H. Villinger. 1990. Analysis of conductive and convective heat transfer in a sedimentary basin, demonstrated for the Rhein graben. Geophys. J. Int. v. 100, p. 393-414.

170. Clayton J.L. and Bostick N.H. Temperature effect on kerogen and on molecular and isotopic composition of organic matter in Pierre shale near an igneous dike. // Org. Geochem., 1986, 10, n. 1/3, p. 135-143.

171. Cloetingh S.,Ben-Abraham,Sass W., and Horvath F. Dynamics of basin formation and strike-slip tectonics. Tectonophysics, 1996, v.266, n.1-4. p. 1-10.

172. Colombo U., Gazzarrini F., Gonfiantini R., Kneuper G., Teichmuller M., Teichmuller R. Carbon isotope study on methane from German coal deposits. in "Advances in Organ. Geochem., 1966", 1970, 1-26.

173. Combarnous M. Natural convection in porous media and geothermal systems. // In Int. Heat Transfer Conf., 6th, 1978, p. 45-59.

174. Correira A. and Jones F.W. On the importance of measuring thermal conductivities for heat flow density estimates: an example from the Jeanne d'Arc Basin, offshore eastern Canada. Tectonophysics, 1996, v. 257, n. 1, p. 71-80.

175. Davey F.J. The Antarctic margin and its possible hydrocarbon potential. Tectonophysics, 1985, v. 114, p. 443-476.

176. Delaney P.T. and Pollard D.D. Solidification of basaltic magma during flow in a dike. Amer. J. Sci., 1982, 282, n. 6, p. 856-885.

177. Deming, D., and D.S.Chapman, 1989, Thermal histories and hydrocarbon generation: Example from Utah-Wyoming thrust belt: AAPG Bull., v.73, N. 12, p. 1455-1471.

178. Deming D., J.A. Nunn, and D.G. Evans, 1990. Thermal effect of compaction-driven groundwater flow from overthrust belts. J.Geophys. Res., v.95, N.B5, p.6669-6683.

179. Deming D., 1994. Factors neccessary to define a pressure seal. AAPG Bull, v.78, p. 10051009.

180. Dunbar J.A. and Sawyer D.S. Three-dimensional dynamical model of continental rift propagation and margine plateau formation. J.geophys.Res. v. 101, n.B 12, p. 27.845-27.863.

181. Durand B., Alpern B., Pittion J.L., and Pradier B. Reflectance of vitrinite as control of thermal history of sediments. In: "Thermal modelling in sedimentary basins", ed. Burrus J., Paris, 1986, p. 441-474.

182. Dykstra J. Compaction correction for burial history curves applications to Lopatin's method for source rock maturation determination. Geobyte, 1987, v.2, N 4, p. 16-23.

183. Dypvick H. Clay mineral transformations in Tertiary and Mesozoic sediments from North Sea. AAPG Bull., 1983, v. 67, n.l, p. 160-165.

184. England W., MacKenzie A., Mann D., and Quigley T., 1987. The movement and entrapment of petroleum fluids in the subsurface. Journal of the Geological Society, London, v. 144, p. 165-180.

185. Espitalie, J., P.Ungerer, I.Irvin, and E.Marquis, Primary cracking of kerogens. Experimenting and modelling CI, C2-C5, C6-C15 classes of hydrocarbons formed: Org.Geochemistry, 1988,. v.13, N.4-6, p.893-899.

186. Fischer K. On the calculation of higher derivatives in finite elements. Comput. Method Appl. Mech. Eng., 1976, v. 7, n.3, p. 320-330.

187. Fisher, A.T. and G. Zwart, 1996. Relation between permeability and effective stress along a plate-boundary fault, Barbadosaccretionary complex. Geology, v. 24, p.307-310.

188. Forbes P.L., Ungerer P., and Mudford B.S., 1992. A two-dimensional model of overpressure development and gas accumulation in Venture field, Eastern Canada. AAPG Bulletin, v.76, n.3, p. 318-338.

189. Forsyth, D.W. and F.Press, 1971, Geophysical tests of petrological models of the spreading lithosphere: J.Geophys.Res., v.76, p.7963-7972.

190. Fowler C.M.R. and Nisbet E.G. The subsidence of the Willinston basin. Canad. J. Earth Science, v. 22, p. 408-415.

191. Frakes L.A. Climates throughout geologic time.NY-Ams. Elsevier. 1979, 310 pp

192. Fredericks P.M., Warbrooke P., and Wilson M.A. A study of the effect of igneous intrusions on the structure of an Australian high volatile bituminous coal. // Org. Geochem., 1985, 8, n. 5, p. 329-340.

193. Gaffin S. Ridge volume dependence on seafloor generation rate and inversion using long term sealevel change. Amer. J. Science., 1987, v.287, p. 596-611.

194. Galushkin Yu.I. Numerical simulation of permafrost evolution as a part of basin modeling: permafrost in Pliocene-Holocene climate history of Urengoy field in West Siberian basin. Canad. J. Earth Science , 1997 , v.34, n.7, 935-948

195. Galushkin Yu.I. The thermal effect of igneous intrusive bodies on maturity of organic matter -A possible mechanism of intrusion formation. Organic Geochemistry, 1997, v. 27, n. 11-12, p. 645-658.

196. Gibbons,M.J., Williams A.K., Piggott N., and Williams G.M. Petroleum geochemistry of the Southern Santes Basin, offshore Brazil. J.geol.Soc., 1983, v. 140, p. 423-430.

197. Gliko A.O. and Maireshal J-C. Non-linear asymptotic solution to Stefan-like problems and the validity of the linear approximation. Geophys. J. Internal., 1989, v.99, p. 801-809.

198. GofF, J.C., 1983, Hydrocarbon generation and migration from Jurassic source rocks in the E-Shethland basin and Viking graben of the Nothern North Sea. J.geol.Soc.Lond., vl40, p.445-474.

199. Gretener, P.E., 1981, Geothermics: usiug temperature in hydrocarbon exploration Education Course, Note Series, 17 Short Course, San Francisco Annual Meeting, May 1981.

200. Gudmundsson A. Emplacement of dikes, sills, and crustal magma chambers at divergent plate boundaries. // Tectonophysics, 1990,176, n. 3/4, p. 257-275.

201. Guidish T.H., Lerche I., Kendall C.G., O'Brien J.J. Relationship between eustatic sea level changes and basement subsidence. AAPGBull., 1984, v. 68, n. 2, p. 164-177.

202. Guptasarma D., N. Purnachandra Rao., Ravikumar M., and Satyabala S.P. The need for geodetic observations to constrain the relative motion of the Indian tectonic plate. Phys. Earth and Planet. Inter., 1991, v. 68, p. 41-44.

203. Hamdani Y., Mareshal I.C., and Arcani-Hamed J. Phase change and thermal subsidence in intracontinental sedimentary basins. Geophys. J. Intemat., 1991, v.106, p. 657-665.

204. Hamdani Y., Mareshal I.C., and Arcani-Hamed J. Phase change and thermal subsidence of the Willinston basin. Geophys. J. Internat., 1994, v. 116, n. 3, p. 585-597.

205. Hanbaba P. and Juntgen H. Zur Ubertragbarkeit von Laboratoriums- Untersuchungen auf geochemische Prozesse der Gasbildung aus Steinkohle und über den Einfluss von Sauerstoff auf die Gasbildung. in"Advan. in Organ. Geochem.,1968", 1969, 459-471

206. Hanbaba P. Juntgen H., Peters W. Nicht-isotherme Reaktionskinetik der Kohlenpyrolyse. Teil II: Entweiterung der Teorie der Gasabspaltung und experimentelle Bestätigung an Steinkohlen. Brennstoff-Chemie, Bd.49, N 12, 368-376, (1968).

207. Hanson R.B. and Barton M.D. Thermal development of low-pressure metamorphic belts: results from two-dimensional numerical models. // J. Geophys. Res., 1989, 94, n. B8, p. 1036310377.

208. Hardee H.C. Permeable convection above magma bodies. // Tectonophysics, 1982, 84, p. 179195.

209. Haxby W.F., Turcotte D.L., and Bird I.M. Thermal and mechanical evolution of the Michigan basin. Tectonophysics, 1976, p. 57-75.

210. Heeremans M., Larsen B.T., and Stel H. Paleostress reconstruction from kinematic indicators in the Oslo Graben, southern Norway: new constraints on the mode of rifting. Tectonophysics, 1996, n. 1-4, p. 55-79:

211. Hegarty K.A., Weissei J.K., and Mutter J.C. Subsidence history of Australian's southern margin: constraints on basin models. AAPGBull., 1988, v. 72, p. 615-633.

212. Herron E.M. and Tucholke B.E. Sea-floor magnetic patterns and basement structure in the southern Pacific. Initial Reports of DSDP Leg.35., 1976.

213. Hilde T.W.C. and Lee C/S. Origin and evolution of the West Philippine basin: a new interpretation. Tectonophysics, 1984, v. 102, p. 85-104.

214. Hirose K. Melting experiments in lherzolite KLB-1 under hydrous conditions and generation of high-magnesian andesite melts. Geology, 1997, v.25, n.l, p. 42-44.

215. Hunt J.M., Lewan M.D., and Henuet R.J-C. Modelling oil generation with time-temperature index graphs based on the Arrhenius equation. AAPG Bull., 1991, v. 75, n. 4, p. 795-807.

216. Husson D.M. and Uyeda S. Tectonic processes and the history of the Mariana arc: a synthesis of the results of Deep See Drilling Project Leg 60. Init.Rep.Deep See Drill. Proj., 1981, v.60, p. 909-929.

217. Hutsberger T.L. and Lerche I. Determination of paleo heat-flux from fission scar tracks in apatite. J. Petrol. Geol., 1987, v. 10, n. 4, p. 365-394.

218. Hutchinson, I., 1985, The effects of sedimentation and compaction on oceanic heat flow Geophys.Roy.astron.Soc. v.82, p.439-459.

219. Ito K. and Kennedy G.C. An experimental study of the basalt-garnet-granulite-eclogite transformation. In "The structure and physical properties of the Earth's crust", Washington, AGU, Geophysic Monogr. 1971, n. 14, p. 303-314.

220. Jame Y-W and D.I. Norum Heat and mass transfer in a freezing unsaturated porous medium //Water resources research, 1980, v. 16, N4, p. 811-819.

221. Jowett E.C., L. M. Cathles III, and B.W. Davis, 1993. Predicting depthsof gypsum dehydratation in evaporitic sedimentary basins. AAPG Bull. v.77, p. 402-413.

222. Juntgen H. and Van Heek K.H. Reaktionsablaufe unter nicht-isothermen Bedingungen. Fortschr.Chem.Forsch., Bd.13. N 3/4, 601-699, (1970).

223. Juntgen H. and Klein J. Entstehung von Erdgas aus kohligen Sedimenten. Erdöl and Kohle, Bd. 28, N2, 65-73,(1975).

224. Juntgen H. and Klein J. Entstehung von Erdgas aus kohligen Sedimenten. Erdöl and Kohle. Ergänzungsband. Erdgas Petrochemie.,Bd. 1, 52-69, (1975).

225. Juntgen H. and Van Heek K.H. An update of German non-isithermal coal pyrolysis work. Fuel Processing Technology, v.2, N.4, 261-293 (1979).

226. Karner G.D. and Watts A.B. On isostasy at Atlantic type continental margins. J.geophys. Res., 1982, v. 87, n. B4, p. 2923-2948.

227. Keen C.E. The dynamics of rifting: deformation of the lithosphere by active and passive driving forces. Geophys. J. Roy. astron. Soc., 1985, v. 80, p. 95-120.

228. Konrad J.M., J.T.C. Seto. Freezing of a clayey silt contaminated within organic solvent//J. of Contaminant Hydrology, 1991, v.8, N 3/4, p.335-355.

229. Kroos B.M. and Leythauser D. Diffusion of methane and ethane through the reservoir cap rocks: implications for the timing and duration of catagenesis; Discussia. AAPG Bull., 1997, v. 81, n. 1, p. 155-161.

230. Kukkonen I.T. Thermal aspects of groundwater circulation in bedrock and its effect on crustal geothermal modelling in Finland, the central Fennoscandian shield. Tectonophysics, 1995, v. 244, n. 1-3, p. 119-136.

231. Kukkonen I.T., Joelecht A. Geothermal modelling of the lithosphere in the central Baltic Shield and its southern slope. Tectonophysics, 1996, v. 255, n. 1-2, p. 25-45.

232. Lachenbruch A.H. Rapid estimation of the topographic disturbance to superficial thermal gradients. Rev. Geophys., 1968, v. 6, n. 3, p. 365-400.

233. Lachenbruch A.H., Sass J.H., Marshall B.V., and Moses T.H. Permafrost, heat flow, and the geothermal regime at Prudhoe Bay, Alaska// J. geophys. Res.,1982, v.87, N B11, p.9301-9316.

234. Langford F.F. and Blaue-Valleron M.-M. Interpreting Rock-Eval pyrolysis data using graphs of pyrolizable hydrocarbons vs. total organic carbon. AAPG Bull., 1990, v. 74, n. 6, p. 799-804.

235. Lebret P., Dupas A., Clet M., et al. Modeling of permafrost thickness during the late glacial stage in France: preliminary results. CanJ.Earth Sci., 1994, v.31, N 6, p. 959-968.

236. Le Pichon X., Angelier J., and Sibuet J-C. Plate boundaries and extensional tectonics. Tectonophysics, 1982, v.81, p. 239-256.

237. Leyden R., Damuth J.A., Ongley L.K. et al., Salt diapirs on Sao Paulo Plateau, Southeastern Braziliancontinental margin. AAPG Bull., 1978, v.62, n.4, p. 657-669.

238. Leyden R., Ludwig W.J., Ewing M. Structure of continental margin of Punta del Este, Uruguay, and Rio de Janeiro, Brasil. AAPG Bull., 1971, v.55, n.12, p. 2161-2173.

239. Littke R., Leythaeuser D., Radke M., Schaefer R.G. Petroleum generation and migration in coal seams of the Carboniferous Ruhr Basin, Northwest Germany. Org.Geochem., v. 16, N.l-3, 247258, (1990).

240. Lister C.R.B. Estimators for heat flow and deep rock properties based on boundary layer theory. Tectonophysics, 1977, v.41, p. 157-171.

241. Lonsdale P. A transform continental margin rich in hydrocarbons, Galf of California. AAPG Bull., 1985, v.69, n.7, 1160-1180.

242. Lopatin N.V., Galushkin Yu.I., and Makhous M. Evolution of sedimentary basins and petroleum formation. In "Geodynamic evolution of sedimentary basins", 1996, ed. TECHNIP, p. 435-453.

243. Luo X. and G. Vasseur, 1992. Contributions of compaction and aquathermal pressuring to geopressure and the influence of environmental conditions. AAPG Bull, v.76, p. 1550-1559.

244. Lutz T.M. and Omar G. An inverse method of modelling thermal histories from apatite fission track data. Earth Planet. Sci. Letters, 1991, v. 104, p. 181-195.

245. MacKenzie A.S. and McKenzie D. Isomerization and aromatization of hydrocarbons in sedimentary basins formed by extension. Geol.Mag., 1983, v. 120, p. 417-470.

246. Magara K. Comparison of porosity-depth relationship of shale and sandstone. J. Petrol. Geology, 1980, v. 3, n. 2, p. 175-185.

247. Makhous M., Galushkin Yu.I., and Lopatin N.V. Modelling of tectonic subsidence and thermal histories in Saharan basins. Terra nova, 1995, v.7, Abstract supplement, N1, p.l 16.

248. Makhous M., Galushkin Yu.I. and N.V. Lopatin. Burial history and kinetic modelling for hydrocarbon generation. Part I: The GALO Model. AAPG Bull. 1997, v.81, n.10, p. 1660-1678.

249. Makhous M., Galushkin Yu.I. and N.V. Lopatin. Burial history and kinetic modelling for hydrocarbon generation. Part II: Application of the Model to Saharan Basins // AAPG Bull. 1997, v.81, n.10, p. 1679-1699.

250. Mann D.M. and A.S. MacKenzie, 1990. Prediction of pore fluid pressures in sedimentary basins. Marine and Petroleum Geology, v.7, p. 55-65.

251. Mareshal J-C. Mechanism of uplift preceding rifting. Tectonophysics, 1983, v. 94, p. 51-66.

252. Mareshal J-C. Uplift and heat flow following the injection of magmas into lithosphere. Geophys. J. Roy. astron. Soc., 1983, v.73, p. 109-127.

253. McKenzie D. P. Some remarks on heat-flow and gravity anomalies. J. geophys. Res., 1967, v.72, n. 24.

254. McKenzie D.P. Some remarks on the development of sedimentary basins. Earth Planet. Science Letters, 1978, v.40, p.25-32.

255. McKenzie, D.P.,The variation of temperature with time and hydrocarbon maturation in sedimentary basins formed by extension. Earth Planet.Science Letters, 1981, v.55, p.87-98.

256. Meissner F.F. Examples of abnormal fluid pressure producted by hydrocarbon generation. (Abs.) AAPGBull., 1980, v.64, p. 749.

257. Minster J.B., Jordan T.H., Molnar P., and Haines E. Numerical modelling of instantaneous plate tectonics. Geophys. J. Roy. astron. Soc., 1974, v.36, p. 541-576.

258. Minster J.B. and Jordan T.H. Present-day plate motions. J. geophys. Res., 1978, v. 83, p. 53315334.

259. Molnar P. and Tapponier P. Cenozoic tectonics of Asia: effect of a continental collision. Science, 1975, v. 189, p.419.

260. Morgan P. and Ramberg I.B. Physical changes in the lithosphere associated with thermal relaxation after rifting. Tectonophysics, 1987, v. 143, p. 1-11.

261. Morrow D.W. and Issler D.R. Calculation of vitrinite reflectance from thermal histories: a comparison of some methods. AAPG Bull., 1993, v. 77, n. 4, p. 610-624.

262. Mrozowski C. L. Lewis S.D., and Hayes D.E. Complexities in thre evolution of the West Phillippine Basin. Tectonophysics, 1982, v. 82, p. 1-24.

263. Mudford B.S., 1988. Modelling the occurrence of overpressures on the Scotian Shelf, Offshore Eastern Canada. J. geophys. Res., v. 93, n. B7, p. 7845-7855.

264. Mudford B.S. and M. Best, 1989. Venture gas field, Offshore Nova Scotia: case study of overpressuring in region of low sedimentation rate. AAPG Bulletin, v. 73, n. 11, p. 1383-1386.

265. Mudford B.S., 1990. A one-dimensional, two phase model of overpressure generation in the Venture gas field, Offshore Nova Scotia. Canad. Petrol. Geology Bulletin, v. 38, n. 2, p. 246258.

266. Nakayama, K. and I.Lerche, 1987, Basin analysis by model simulation: effects of geologic parameters on ID and 2D fluid flow systems with application to an oil field: Golf Coast Association of Geological Societies Transactions, v.37, p.175-184.

267. Neugebauer H.J. Mechanical aspects of continental rifting. Tectonophysics, 1983, v. 94, p. 91108.

268. Newman R. and White N. Rheology of the continental lithosphere inferred from sedimentary basin. Nature, 1997, v. 385, February, p. 621-624.

269. Nielsen, S.B., Balling N., 1990, Subsidence, heat flow, and hydrocarbon generation in extensional basins. First break, v 8, N 1, p.23-31.

270. Nikishin A.M., Cloetingh S., Ershov A.V. et al., Dnieper-Donets basin, Ukrainian shield and Voronezh high : Carbonifereous-Cenozoic movements as reaction to reorganisations of regional stress field. Tectonophysics, 1997 (in press).

271. Nixon J.F. Thermal simulation of subsea saline permafrost// Canad.J.Earth Sci., 1986, v.23, p.2039-2046.

272. Parson B. and Sclater I.C. An analysis of the variation of oceanic floor bathimetry and heat flow with age. J. geophys. Res., 1977, v. 82, p. 803-820.

273. Peck D.L., Hamilton M.S., and Shaw H.R. Numerical analysis of lava lake cooling models. Part H: Application to Alal lava lake, Hawaii. // Amer. J. Sci. 1977, 277, n. 4, p. 415-457.

274. Pedersen T. Some remarks on lithospheric forces and decompression magmatism. Tectonophysics, 1994, v. 240, n. 1-4, p. 11-19.

275. Peaceman, D.W. and H.H.Rachford, 1955,. The numerical solution of parabolic and elliptic differential equations. Joura. Society of Industrial and Applied Mathematics v. 3, N 1, p. 28-41.

276. Perregard J. and Schiener E.J. Thermal alteration of sedimentary organic matter by a basalt intrusive (Kimmeridgian shales, Milue Land, East Greenland). // Chemical Geology, 1979, 26, n. 3/4, p. 331-343.

277. Perrier, B., and J.Quiblier, 1974, Thickness changes in sedimentary layers during compaction history: methods for quantitative evaluation, AAPG Bull., v.58, N.3, p.507-520.

278. Person M. and G. Garven. 1992. Hydrologie constraints on petroleum generation withincontinental rift basins: theory and application to the Rhine Graben. AAPG Bull.v.76, N 4, p.468-488.

279. Peters K.E., Whelan J.K., Hunt J.M., and Tarafa H.F. Programmed pyrolysis of organic matter from thermally altered Cretaceous black shales. // AAPG Bull., 1983, 67, n. 11, p. 21372149.

280. Petit C., Burov E., and Deverckere J. On the structure and mechanical behaviour of the extending lithosphere in the Baikal Rift from gravity model. Earth Planet. Sci. Letters, 1997, v. 149, n. 1-4, p. 29-42.

281. Phipps Morgan J. and Chen Y.J. The genesis of oceanic crust: magma injection, hydrothermal circulation, and crustal flow. // J. Geophys. Res. 1993, 98, n. B4, p. 6283-6297.

282. Phipps Morgan J, Parmentier E.M., and Lin J. Mechanisms for the origin of mid-oceanic ridge axial topography: Implication for the thermal and mechanical structure at accreting plate boundaries. // J. Geophys. Res., 1987, 92, ser.B, p. 12823-12836.

283. Press, W.H., B.P.Flannery, S.A.Teukolsky, and W.T.Vetterling, 1986, Numerical recipes-The art of scientific computing: Cambridge, Cambridge University Press, 818 p.

284. Pytte M. and Reynolds R.C. The thermal transformation of smectite to illite. In: "Thermal history of sedimentary basins" Naeser N.D. and McCulok T.H. (eds.), Springer Verlag, 1989, p. 132-140.

285. Rabinovich P. and Cooper A. Structure and sediment distribution in the Western Bering Sea. Marine Geol., 1977, v. 24, p. 309-320.

286. Raymond A.C. and Murchison D.G. Effect of volcanic activity on level of organic maturation in Carbonifereous rocks of East Fife, Midland Valley of Scotland. Fuel, 1988a, 67, n. 8, p. 11641166.

287. Raymond A.C. and Murchison D.G. Development of organic maturation in the thermal aureoles of sills and its relation to sediment compaction. Fuel, 1988b, 67, n. 12, p. 1599-1608.

288. Raymond A.C. and Murchison D.G. Organic maturation and its timing in a Carbonifereous sequence in the central Midland Valley of Scotland: comparisons. Fuel, 1988,68,n. 3, p. 328-334.

289. Rehault J-P., Tisseau Ch., Brunet M-F, and Londen K.E. Subsidence analysis on the Sardinian margin and the Central Tyrrenian basin: thermal modelling and heat flow control; deep structure implications. J. of Geodynamics, 1990, v. 12, 269-310.

290. Rouden L. and Keen C.E. Rifting processes and thermal evolution of the continental margin of eastern Canada determined from subsidence curves. Earth Planet.Science Letters, 1980, v.51, p.343-361.

291. Roussel J. and Linger J.L. A review of deep structure and ocean-continent transition in the Senegal basin (West Africa). // Tectonophysics, 1983, 91, p. 183-211.

292. Rowley D.B. and Sahagian D. Depth-dependent stretching a different approach. Geology, 1986, v. 14, p. 32-35.

293. Rudke M., Horsfeld B., Littke R., and Bullkotter J. Maturatiion and petroleum generation. In: "Petroleum and basin evolution", Welte D.H., Horsfield B., and Baker D.R. (eds.),.Springer-Verlag, 1997, p. 171-229.

294. Ryan P.D. and Dewey J.F. Continental eclogites and Wilson cycle. J. geolog. Soc. London, 1997, v. 154, p. 437-442.

295. Rybach L. Heat sources, heat transfer and rock types in the lower continental crust-inference from Deep Drilling. Tectonophysics, 1996, v. 257, n. 1, p. 1-6.

296. Sajgo Cs. and Lefler J. A reaction kinetic aproach to the temperature-time history of sedimentary basins. In: "Paleogeothermics. Lecture Notes in Earth Sciences", Buntebarth F. and Stegena I. (eds.) Springer, 1986, v. 5, p. 119-151.

297. Saxby J.D., Bennett A.J.R., Corcoran J.F., Lambert D.E., RilleyK.W. Petroleum generation: simulation over six years of hydrocarbon formation from torbanite and brown cole in a subsiding basin. Org.Geochem., v.9, N.2, 69-81, (1986).

298. Schatz J.F. and Simmons G. Thermal conductivity of Earth materials at high temperatures. J. geophys. Res., 1972, v. 77, n. 35, p. 6966-6983.

299. Schaw H.R., Hamilton M.S., and Peck DL. Numerical analysis of lava lake cooling models. Part I: Description of the method. // Amer. J. Sci. 1977, 277, n. 4, p.384-414.

300. Sclater J.G. and Celerier B. Extensional models for the formation of sedimentary basins and continental margins. Norsk Geologisk. Tidsskrift, 1987, v. 67, p. 253-267.

301. Sclater, J.G., and P.AF.Christie, 1980, Continental stretching: an explanation of the Post-Midcretaceous subsidence of the central North sea basin: J.Geophys.Res.,v.85,N.B7, p.3711-3739.

302. Sclater J.G., Parsons B., and Jaupart C. Ocean and continent similarities and differences in the mechanisms of heat loss. // J. Geophys. Res., 1981, 86, n. B12, p. 11535-11552

303. Shemenda A.I. Subduction of lithosphere and back-arc dynamics: insight from physical modelling. J. Geophys. Res., 1993, v. 98 (Ser.B), 16167-16185.

304. Shi Ya. and Wang C-Yu., 1986. Pore generation in sedimentary basins: overloading versus aquathermal. J. geophys. Res. v. 91, n. B2, p. 2153-2162.

305. Simoneit B.R.T., Brenner S., Peters K.E., and Kaplan I.R. Thermal alteration of Cretaceous black shale by basaltic intrusions in the Eastern Atlantic. // Nature, 1978, 273, n. 5663, p. 501504.

306. Simoneit B.R.T., Brenner S, Peters K.E., and Kaplan I.R. Thermal alteration of Cretaceous black shale by diabase intrusions in the Eastern Atlantic. Part n. Effects on bitumen and kerogen. //Geochem. Cosmochim. Acta, 1981, 45, n. 9, p. 1581-1602.

307. Sleep N.H. and M.L. Blanpied, 1992. Creep, compaction and the weak rheology of major faults. Nature, v.359, 687-692.

308. Sloan, E.D. 1990. Clathrate hydrates of natural gases. Marcel Dekker, Inc., New York and Basel, 234 p.

309. Spadini G., Robinson A., and Cloetingh S. Western versus Eastern Black Sea tectonic evolution : pre-rift lithospheric controls of basin formation. Tectonophysics, 1996, v/266, n.1-4, p. 139154.

310. Spohn T. and Schubert G. Convective thinning of the lithosphere: a mechanism for rifting and midplate volcanism on Earth, Venus, and Mars. Tectonophysics, 1983, v. 94, p. 67-90. 1983

311. Stephenson RA. Continental rift development in Precembrian and Phanerozoic Europe: EUROPROBE and the Dnieper-Donets rift and Polish trough basin. Sedimentary Geology. 1993, v.86, p. 159-175.

312. Stillman C.J., Furnes H., LeBas M.J., Robertson A.H.F., and Zielonka. J. The geological history of Maio, Cape Verde Islands. // J. Geol. Soc. London, 1982,139, p.347-361.

313. Stockmal, G.S., C. Beaumont, R. Boutilier, 1986, Geodynamic models of convergent margin tectonics: transition from rifted margin to overthrust belt and consequences for foreland-basin development, AAPGBull, v 70, N 2, p 181-190.

314. Su D., White N., and McKenzie D. Extension and subsidence of the Pearl River Month basin, northern South China Sea. Basin Research, 1989, v. 2, p. 205-222.

315. Sweeney, J.J., and A.K.Burnham, 1990, Evolution of a simple model of vitrinite reflectance based on chemical kinetics: AAPG Bull, v.74, n.10, p. 1559-1570.

316. Sweeney J.J., Braun R.L., Burnham A.K., Talukdar S., and Vallejos C., 1995. Chemical kinetic model of hydrocarbon generation, expulsion, and destruction applied to the Maracaibo Basin, Venezuela. AAPG Bulletin, v. 79, n.10, p. 1515-1532.

317. Takeshita T. and A. Yamaji. Acceleration of continental rifting due to a thermomechanical instability. Tectonophysics, 1990, v. 181, p. 307-320.

318. Tapponier P., Pelzer G., and Armizo R. On the mechanics of the collision between India and Asia. In: "Collision Tectonics", eds. M.P. Coward and A.C. Ries, Geol. Soc. Spec. Publ. v. 19, p. 115-157.

319. Ter Voorde M and Bertotti G. Thermal effects of normal faulting during rifted basin formation. 1. A finite difference model. Tectonophysics, 1994, v. 240, n. 1-4, p. 133-144.

320. Tissot, B.P.,and J.Espitalie, 1975, L'évolution thermique de la matiere organique des sediments applications d'une simulation mathématique: Revue de l'InstitutFranccais du Petrole, v.30, p.743-777.

321. Tissot, B.P., RPelet, and P.Ungerer, 1987, Thermal history of sedimentary basins, maturation indices, and kinetics of oil and gas generation: AAPG Bull., v.71, N.12, p. 1445-1466.

322. Touloukian Y.S. and Ho C.Y. Physical properties of rocks and minerals. McGrew-Hill, 1981, 548 p.

323. Toupin D., Eadington P. J., Person M., et al. Petroleum hydrology of the Cooper and Eromanga basins, Australia: some insights from mathematical modelling and fluids inclusion data. AAPG Bull., 1997, v. 81, n. 4, p. 577-603.

324. Turcotte D.L. Magma migration. Ann. Rev. Earth Planet. Sci. 1982,10, n. 1557, p. 397-408,

325. Turcotte D.L. and Bernthal M.J. Synthetic coral-reef terraces and variations of Quaternary sea level. Earth.Planet.Sci.Letters, 1984, v. 70, p. 121-128.

326. Turcotte,D.L.,and G.Schubert, 1982, Geodynamics: Application of continuum physics to geological problems: New York, John Wiley & Sons,v.l, 376 p.

327. Ujie Y. Thermal alteration of kerogen as an indicator of contact metamorphism to sedimentary rocks: H-NMRTj and element composition.// Geochem. J., 1984,18, p. 163-166.

328. Ungerer, P., 1990, State of the art of research in kinetic modelling of oil formation and expulsion: Org.Geochemistry, v. 16, N.l-3, p. 1-27.

329. Ungerer, P., I.Burrus, B.Doligez, P.Chenet, and F.Bessis, 1990, Basin evolution by integrated two-dimensional modelling of heat transfer, fluid flow, hydrocarbon generation, and migration: AAPG Bull., v.74, N.3, p.309-335.

330. Ushupi E., Emery K.O., Bowin C.O., and Phillips J.D. (1976) Continental margin off Western Africa: Senegal to Portugal. American Association of Petroleum Geologist Bull, 60, n. 5, p. 809-879.

331. Vail P.R., Mitchum R.M., and Thompson S. Global cycles of relative changes of sea level. AAPG Mem. , 1977, n. 26, p. 83-97.

332. Van der Beek P. Flank uplift and topography at the central Baikal Rift (SE Siberia): a test of cinematic models for continental extension. Tectonics, 1997, v. 16, n. 1, p. 122-136.

333. Van der Linden, W.I.M. The crustal structure and evolution of the continental margin off Senegal and Gambia, from total-intensity magnetic anomalies. // Geologie en Mijnbouw, 1981, 60 , n. 2, p. 257-266.

334. Van Heek K.H., Juntgen H., Luft K-F., Teichmuller M. Aussagen zur Gasbildung in fruhen Ikohlungsstadian auf grand von Pyrolyseversuchen. Erdol, Kohle, Erdgas, Petrochemie. 1971, Bd.24, N.9, 566-572;

335. Vasseur G., Brigand F., and Demongodin L. Thermal conductivity estimation in sedimentary basins. Tectonophysics, 1995, v. 244, n. 1-3, 167-174.

336. Veevers J.J. Breakup of Australia and Antarctica estimated as mid-Cretaceous (95+. 5 Ma) from magnetic and seismic data at the continental margin. Earth Planet. Sci. Letters? 1986, v.77, p. 91-99.

337. Wang X., Lerche L., and Walters C. The effect of igneous intrusive bodies on sedimentary thermal maturity. // Org. Geochem., 1989, 74, n. 6, p. 571-584.

338. Waples D.W. 1980. Time and temperature in petroleum formation: application of Lopatin's method to petroleum exploration AAPG Bull., v.64, p.916-926.

339. Waples D.W. Geochemistry in petroleum exploration. N-Y, 1985, 456 p.

340. Waples D.W., Kamata H., and Suizu M. 1992. The art of maturity modelling. Part 1 : Finding of satisfactory geological model. AAPG Bull, v.76, N.l, p.30-46.

341. Ward S.N. Pacific-North America plate motions: new results from very long baseline interferometry. J. geophys. Res., 1990, v. 95, n. B13, p. 21,965-21,981.379.Watts, Stackler 1979378

342. Welte D.H. and Yukler M.A. Petroleum origin and accumulation in basin evolution a quantitative model. AAPGBull., 1981, v. 65, n. 8, p. 1387-1396.

343. Welte, D.H., and M.M.Yalcin, 1988, Basin modelling a new comprehensive method in petroleum geology: Advances in Organic Geochemistry, v. 13, p. 141-151.

344. Welte D.H., Schaefer R.G., Yalkin M.N. Gas generation from source rocks: aspects of a quantitative treatement. Chemical Geology, v.71, 105-116, (1988).

345. Welte D.H., Horsfield B., and Baker D.R. (eds.). Petroleum and basin evolution.Springer-Verlag, 1997, 536 p.

346. Wilkinson M., Secondary porosity generation during deep burial associated with overpressure leak-off: Fulmal formation, United Kingdom, Central Graben. AAPG Bull., 1997, v.81, n.5, p.803-813.

347. Wyllie, P.J., 1979, Magmas and volatile components: Am.Mineral, v.64, p.469-500.

348. Zorin Yu. and Lepina S.V. Geothermal aspects of development of asthenospheric upwelling beneath continental rift zones. Journal of Geodynamics, 1989, v.3, p. 1-22.

349. Yamano M., Uyeda S., Sibuet J.C., and Foucher J.P. Heat flow anomaly in the middle Okinawa trough. Tectonophysics, 1989, v. 159, p. 307-318.