CC BY
27
УДК 553.98(571)
КОМПЛЕКСНАЯ ОЦЕНКА ПАЛЕОКЛИМАТИЧЕСКИХ ФАКТОРОВ РЕКОНСТРУКЦИИ ТЕРМИЧЕСКОЙ ИСТОРИИ НЕФТЕМАТЕРИНСКОЙ БАЖЕНОВСКОЙ СВИТЫ АРКТИЧЕСКИХ РАЙОНОВ ЗАПАДНОЙ СИБИРИ
Исаев Валерий Иванович1,
isaevvi@tpu.ru
Искоркина Альбина Альбертовна1,
iskorkina.a@mail.ru
Косыгин Владимир Юрьевич2,
kosyginv@inbox.ru
Лобова Галина Анатольевна1,
lobovaga@tpu.ru
Осипова Елизавета Николаевна1,
osipovaen@list.ru
Фомин Александр Николаевич3,
fominan@ipgg.sbras.ru
1 Национальный исследовательский Томский политехнический университет, Россия, 634050, г. Томск, пр. Ленина, 30.
2 Вычислительный центр ДВО РАН,
Россия, 680000, г. Хабаровск, ул. Ким Ю Чена, 65.
3 Институт нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН, Россия, 630090, г. Новосибирск, пр. Ак. Коптюга, 3.
Актуальность. Арктические районы Западной Сибири имеют уникальные палеоклиматические особенности, инверсионную тектоническую историю в палеоген-неогене, значительные вариации геохимических характеристик нефтематеринских отложений. Это обусловливает необходимость совершенствовать схемы и параметры количественной оценки ресурсов углеводородов объемно-генетическим методом, основанным на палеореконструкциях геотермического режима материнских отложений. Цель исследований: оценить влияние факторов палеоклимата - векового хода температур на земной поверхности и неоплейстоценовых толщ мерзлоты, ледниковых покровов, в комплексе с фактором эрозионных процессов, - на расчетный геотермический режим нефтематеринской баженовской свиты.
Объект исследований: баженовские отложения мезозойско-кайнозойского разреза, вскрытые глубокими скважинами на Арктической, Среде-Ямальской и Ростовцевской площадях (п-ов Ямал).
Методика исследований базируется на оригинальном компьютерном палеотемпературном моделировании, учитывающем параметры седиментационной истории и истории теплофизических свойств осадочной толщи, включающей вечномерзлые породы и ледники, и не требующем априорных сведений о величинах и природе глубинного теплового потока. Результаты исследований по выявлению и учету параметров денудаций, в комплексе с оценкой роли и учетом палеоклиматических характеристик, позволили сделать выводы и дать рекомендации. При определении ресурсов углеводородов объёмно-генетическим методом на землях арктического региона рекомендуется применять «арктический» вековой ход температур и учитывать динамику толщи неоплейстоценовой мерзлоты мощностью 300-600 м. В случае неучета вечной мерзлоты и палеоклима-тического хода температур расчетные ресурсы УВ могут быть занижены до 40-50 %. Сформулированы задачи предстоящих исследований, среди которых оценка роли палеоклиматических факторов в термической истории нефтематеринской китербют-ской свиты и разработка формализованных сценариев выявления и количественной оценки перерывов осадконакопления и размывов по данным геотермии.
Ключевые слова:
Палеоклимат, инверсионная тектоника, геотермический режим, нефтематеринские баженовские отложения, п-ов Ямал.
Введение
В рамках новой парадигмы развития сырьевой базы углеводородов РФ главными объектами прогнозирования и поисков становятся арктические районы Западной Сибири и сланцевые ресурсы баженовской свиты. Для реализации новой парадиг-
мы требуется создание новой фундаментальной теоретической базы и технологий поисков [1].
Оценка перспектив нефтегазоносности - оценка плотности ресурсов углеводородов и районирование территорий, зон нефтегазонакопления - выполняется объемно-генетическим методом (бассей-
новое моделирование). Количество генерированных углеводородов (УВ) рассчитывается на основе реконструкции геотемпературного режима нефте-материнских отложений [2, 3].
Новейшими исследованиями методом палео-температурного моделирования, на примере нефте-материнских баженовских отложений, вскрытых скважиной Арктическая 11 на п-ове Ямал, установлено [4], что неучет векового хода температур на поверхности Земли и толщи неоплейстоценовой мерзлоты не позволяет адекватно восстановить термическую историю материнских отложений. При определении ресурсов УВ объёмно-генетическим методом на землях арктического региона Западной Сибири рекомендовано применять «арктический» вековой ход температур и учитывать динамику толщи неоплейстоценовой мерзлоты.
В публикации Арктической экспедиции IODP 302 [5] приводятся результаты построения и анализа объединенной литосферно-бассейновой термальной модели (система PetroProb) в пределах Хребта Ломоносова. Для учета палеоклиматического фактора авторами построен «местный» (для района Хребта Ломоносова) вековой ход температур на земной поверхности, начиная со 100 млн лет назад. Для этого использовались результаты палинологического и изотопного анализа углерода органического вещества. Здесь следует отметить, что «местный» вековой ход температур на земной поверхности [5], в интервале геологического времени 100-4 млн лет назад, идеально совпадает с «арктическим» вековым ходом температур на земной поверхности [4], независимо построенным в результате сводки и обобщения экспериментальных определений и палеоклиматических реконструкций для севера Западно-Сибирской низменности.
В заключении S. Nelskamp с соавторами отметили, что эволюционирование температур на земной поверхности оказывает большое влияние на зрелость нефтематеринской породы: в зависимости от временных вариаций поверхностных температур могут быть большие или меньшие объемы получаемых УВ.
Основная цель настоящих исследований -дальнейшая аргументация существенного влияния 2-х факторов палеоклимата - векового хода температур на земной поверхности и неоплейстоценовых толщ мерзлоты - на расчетный геотермический режим нефтематеринских свит, определяющий подсчет плотности ресурсов. Для этого выполнены палеотемпературные исследования ба-женовских отложений, вскрытых глубокими скважинами на Ростовцевском и Средне-Ямальском месторождениях п-ова Ямал.
Вместе с тем представляет интерес количественная оценка влияния 3-го фактора палеоклимата - позднечетвертичных ледниковых покровов - на расчетный геотермический режим и степень реализации генерационного потенциала материнских отложений Западной Сибири. Опреде-
ление роли ледниковых покровов приобретает особую актуальность при оценке перспектив нефтега-зоносности арктических территорий Западной Сибири [3]. Для этого выполнено построение и приложение векового хода мощностей ледниковых покровов в четвертичное время.
Построение векового хода мощностей мерзлых пород в неоплейстоцене (северная палеоклиматическая зона Западной Сибири)
Вековой ход мощностей ледниковых покровов (северные широты 66-74°) построен на основе обобщения результатов детальных палеогеографических, стратиграфических исследований и палео-криологических реконструкций для севера Западно-Сибирской низменности, начиная с развития процессов оледенения в верхнем плейстоцене.
Отдел Подотдел Мощности ледников, м 50010001500 20002500 млн. лет
Плейстоцен~~~~—[олоткн Неоплейстоцен >а 0,02 0,04 0,06 0,08 0,10 0,12 0,14 0,16 0,18 0,20
1
< ¡г-. 18 —.
и 8
•
00 i
ч ( | )
0
0 1 020 30405 0607080 9
Рис. 1. Мощности ледников северной палеоклиматической зоны Западной Сибири в неоплейстоцене: 1-8 значение мощности ледников: 1 - по Глобальные изменения ... [7], 2 - по С.А. Архипову и др. [8], 3 - по С.В. То-мирдиаро [9], 4 - по С.А. Лаухину [10], 5 - по Четвертичные оледенения... [11]; 6 - по Е.В. Артюшкову [12];
7 - по В.В. Баулину [6], 8 - по Ю.Б. Баду [13]; 9 - кусочно-линейная аппроксимация векового хода мощностей ледников
Fig. 1. Capacities of glaciers of a northern paleoclimatic zone of Western Siberia in Pleistocene: 1-8 the value of power of glaciers: 1 - on Global changes... [7], 2 - after S.A. Arkhi-pov et al. [8], 3 - after S.V. Tomirdiaro [9], 4 - after S.A. Laukhin [10], 5 - on the Quaternary freezing... [11]; 6 - after E. V. Artyushkov [12]; 7 - after V. V. Baulin [6],
8 - after Yu. B. Badu [13]; 9 is the piecewise and linear approximation of the secular variation of glaciers capacity
На рис. 1 изображен вековой ход мощностей ледниковых покровов в четвертичное время для северной палеоклиматической зоны Западной Сиби-
ри, начиная с развития процессов оледенения в са-маровское время [6] и заканчивая деградацией ледников в конце сартанского времени [7].
66°00' 72°00'в.д.
Сеяха ® 1 • 2
Арктическое 4 5 Арк-11 •
Рис. 2. Обзорная схема территории исследований: 1 - населенный пункт и его название; 2 - поисково-разведочная скважина; 3 - сейсмический профиль работ МОГТ 2D; 4 - контур месторождения и его название (объекты экспериментальных исследований); 5 - гидрография и береговая линия; 6 - моделируемая скважина и ее индекс: Арк-11 - Арктическая 11, С-Я-14 - Средне-Ямальская 14, Рос-64 - Ростовцевская 64
Fig. 2. Review scheme of the studied area: 1 is the locality and its name; 2 is the exploration well; 3 is the seismic profile of works MOGT 2D; 4 is the contour of the field and its name; 5 is the hydrography and coastline; 6 is the well simulated and its index: Арк-11 - Arkticheskaya 11, С-Я-14 - Sredne-Yamalskaya 14, Рос-64 - Rostovtsev-skaya 64
Объект исследований
Моделирование палеогеотемпературных условий нефтематеринских баженовских отложений выполнено для палеоклиматических и геолого-геофизических условий арктических районов Западной Сибири в ходе исследований на площадях глубокого бурения месторождений углеводородов п-ова Ямал, расположенных на разных широтах (рис. 2, табл. 1).
Осадочный мезозойско-кайнозойский чехол территории исследования начинает формироваться в ранней юре. К концу волжского века трансгрессия моря расширилась, идет накопление баже-новской свиты обогащенной органиче-
ским веществом.
Нефтематеринская баженовская свита является основным источником формирования залежей УВ в ловушках верхнеюрского и мелового нефтегазоносных комплексов (НГК), а также приоритетной сланцевой формацией. Однако, в отличие от юго-западных районов Западно-Сибирской нефтегазоносной провинции [14], баженовские отложения в арктической зоне существенно варьируют как по концентрациям рассеянного органического вещества (РОВ), нередко уменьшаясь до 1-2 %, так и по типу РОВ, переходя к гумусово-сапропе-левому типу.
Начиная с апт-сеномана, морской режим господствует до начала эоцена. Раскрытие котловины Арктического бассейна приводит к смене знака вертикальных тектонических движений и наступает позднеэоценовая регрессия. Анализ мощностей палеоген-неогена [15] показывает, что кровля люлинвора (ирбита) могла быть подвергнута денудации. При этом размытый слой мог составить порядка 500-700 м.
Дискуссия о перерывах в осадконакоплении, эрозиях на рубеже позднего палеогена и эоцена в арктической зоне Западной Сибири продолжается по сей день [16-18]. Наличие перерывов в осадко-накоплении и размывов установлено однозначно, однако их объем и временные привязки однозначно не определены.
Анализ предварительно выполненных нами па-леотектонических и палеотемпературных реконструкций свидетельствует о том, что накопление шло до середины миоцена на протяжении 31,9 млн лет (нюрольская, тавдинская, атлымская, новоми-хайловская, туртасская, абросимовская свиты) в объеме 335-535 м и за 4 млн лет, в раннеби-щеульское время, эти отложения были размыты.
Ингрессиям бореального моря в среднем миоцене - раннем плиоцене, с конца бищеульского времени и до конца новопортовского, обязаны накопления осадков толщиной 113 м, которые в последующий этап положительных тектонических движений, за 1,3 млн лет, денудируются. С началом позднего миоцена идет накопление плиоцен-четвертичных озерно-аллювиальных осадков.
Таблица 1. Характеристика разрезов скважин Арктического, Средне-Ямальского и Ростовцевского месторождений Table 1. Characteristic of the sections of Arcticheskaya 11, Sredne-Yamalskaya 14 and Rostovtsevskaya 64 wells
Значение/Value
Характеристики/Characteristics Скважина/Well
Арктическая 11 Arcticheskaya 11 Средне-Ямальская 14 Sredne-Yamalskaya 14 Ростовцевская 64 Rostovtsevskaya 64
Забой, м/Bottom, m 3624 3383 3485
Отложения на забое (свита) Sediments at the bottom (suite) Левинская (J,lv) Levinskaya Надояхская (J1+J2nd) Nadoyahskaya Леонтьевская (J2ln) Leontjevskaya
Кровля баженовской свиты (J3+K,bg), м Roof of the Bazhenov suite (J3+K,bg), m 2792 2855 3162
Мощность баженовской свиты, м Capacity of the Bazhenov suite, m 16 35 33
Мощность палеогеновых отложений в современном разрезе, м Capacity of the Paleogene deposits in the current section, m 223 350 400
Мощность неоген-четвертичных отложений в современном разрезе, м Capacity of the Neogene-Quarternary deposits in the current section, m 280 200 266
Размыв палеоген-неогеновых отложений (14,5-18,5 млн лет назад), м Scour of the Paleogene-Neogene deposits (14,5-18,5 million years ago), m 535 435 335
Размыв неогеновых отложений (4,1-5,4 млн лет назад), м Scour of the Neogene deposits (4,1-5,4 million years ago), m 113 113 113
Мощность вечномерзлых пород в плиоцен-квартере (0,52-0,18 млн лет назад), м Capacity of the permafrost rocks in the Pliocene-Quarternary (0,52-0,18 million years ago), m 600 600 600
Мощность вечномерзлых пород в плиоцен-квартере (0,18-0,0 млн лет назад), м Capacity of permafrost rocks in the Pliocene-Quarternary (0,18-0,0 million years ago), m 300 300 300
Мощность ледникового покрова в неоплейстоцене (0,18-0,13 млн лет назад), м Capacity of glacial cover in the Neo-Pleistocene (0,18-0,13 million years ago), m 500 500 500
Мощность ледникового покрова в неоплейстоцене (0,13-0,05 млн лет назад), м Capacity of glacial cover in the Neo-Pleistocene (0,13-0,05 million years ago), m 1500 1500 1500
Мощность ледникового покрова в неоплейстоцене (0,05-0,015 млн лет назад), м Capacity of glacial cover in the Neo-Pleistocene (0,05-0,015 million years ago), m 500 500 500
Ахская; БЯ9-10; газ, 210 тыс.;
конденсат, 35,2; вода, 1,5 Akhskaya; BYa9-10; gas, 210 thousand; condensate, 35,2; water, 1,5
Ахская; БЯ9-10; вода, 3,9 Akhskaya; BYa9-10; water, 3,9 Ахская; НП2-3; газ, 6,4 тыс.
Результаты испытаний верхнеюрских-нижнемеловых пластов (свита; пласт; тип флюида; дебит, м3/сут) Results of the test of the Upper Jurassic-Lower Cretaceous layers (suite; layer; fluid type; output, m3/d) Не испытывались Not tested Akhskaya; NP2-3; gas, 6,4 thousand Ахская; НП5-6; газ, 11,5 тыс.; нефть, 6,5 Akhskaya; NP5-6; gas, 11,5 thousand; oil, 6,5 Ахская; НП5-7; нефть, 4,2 Akhskaya; NP5-7; oil, 4,2 Ахская; Ач-1; сухо Akhskaya; Ach-1; dry Баженовская; сухо Bazhenov; dry
Измеренные пластовые температуры (свита; глубина замера, м; температура, °С) Measured reservoir temperatures (suite; measurement depth; temperature, °C) Левинская/Levinskaya; 3533; 125 Левинская/Levinskaya; 3560; 126 Марресалинская Marresalinskaya; 846; 17 Ахская/Akhskaya; 2470;75 Ахская/Akhskaya; 2650; 81 Ахская/Akhskaya; 2660; 81
«Измеренные» температуры по ОСВ (глубина отбора, м; (R0t); температура, °С) Measured temperatures by reflective ability of a vitrinit (selection depth; (fl°t); temperature, °C) 2000; (0,65); 100 2500; (0,80); 120 1700;(0,52); 83 2200;(0,65); 100 3000; (0,80); 120 2096 м;(0,53); 84 2600; (0,63); 98 2827; (0,73); 111
В случае учета эрозионных процессов при моделировании по критерию «невязки» измеренных (наблюденных) и расчетных геотемператур результаты являются оптимальными: «невязки» составляют порядка ±2 °С, размывы палеоген-неогено-
вых отложений (335-535 м) и неогеновых отложений (113 м) подтверждаются. А результаты моделирования без учета эрозионных процессов однозначно неприемлемы, «невязки» составляют ±(5-12-16)°С.
О методике исследований
Для исследований применен метод палеотемпе-ратурного моделирования, основанный на численном решении уравнения теплопроводности горизонтально-слоистого твердого тела с подвижной верхней границей, реализованный в программном обеспечении TeploDialog [14, 19, 20].
Расчет палеотемператур состоит из двух этапов. На первом, по распределению температур Т, «наблюденных» в точках разреза скважины, рассчитывается тепловой поток q через поверхность основания осадочного чехла, т. е. решается обратная задача геотермии - классическая обратная задача геофизики. На втором этапе, с известным значением q решаются прямые задачи геотермии - непосредственно рассчитываются температуры и в заданных точках осадочной толщи 1 (в том числе в материнской свите) на любые заданные моменты геологического времени Ь.
Для решения обратной задачи геотермии используем в качестве «наблюденных» как измерения пластовых температур, полученные при испытаниях скважин, так и палеотемпературы, определенные по отражательной способности витринита (ОСВ) (табл. 1).
Первое краевое условие модели определяется температурой поверхности осадконакопления, т. е. фактором палеоклимата, и задается в виде кусочно-линейной функции «арктического» векового хода температур на поверхности Земли [4].
Параметризация осадочного разреза, вскрытого скважиной, определяющая параметры седимен-тационной и теплофизической модели, принимается в соответствии с литолого-стратиграфической разбивкой скважины (табл. 2).
Наличие перерывов в осадконакоплении и величины денудации подтверждены и оценены путем многовариантных расчетов при различных сценариях тектоно-седиментационной истории и последующего выбора сценария, отвечающего критерию «невязки».
Формирование, существование, деградация толщи вечномерзлых пород [4] и ледникового покрова (рис. 1) учитываются как своеобразные динамичные литолого-стратиграфические комплексы, обладающие аномально высокими значениями теплопроводности X и температуропроводности а. А накапливающийся материал ледника представлен аномально низкой плотностью а.
Формализованный учет толщи мерзлоты осуществляется начиная с 0,52 млн лет назад «мгновенной» (по меркам геологического времени, за 0,3+3,0 тыс. лет) заменой «нормальных» осадочных отложений 600-метровой толщей мерзлых пород со своими теплофизическими характеристиками. Затем эта толща мерзлых пород перекрывает осадочный чехол в течение 334 тыс. лет. Далее «мгновенно» (0,3+0,15+0,15 тыс. лет) толща мерзлых пород деградирует в объеме 300 м. И далее
вечномерзлые породы существуют в объеме 300 м до настоящего времени, в последние 182 тыс. лет.
Учет толщи ледникового покрова осуществляется начиная с 0,182 млн лет назад. За 5 тыс. лет сформировался 500-метровой ледник со своими те-плофизическими характеристиками и просуществовал в таком объеме 47 тыс. лет. Затем, за 10 тыс. лет, мощность ледника выросла до 1500 м, и в таком объеме он существует 70 тыс. лет. Далее, за 10 тыс. лет, ледник сокращается до 500 м и существует в течение 20 тыс. лет. К моменту 15 тыс. лет назад (конец сартанского времени) ледник полностью деградирует.
Решение прямых задач геотермии выполнялось на 52 ключевых момента геологического времени, соответствующих временам начала/завершения формирования каждой свиты, перекрывающей материнскую баженовскую, а также точкам «излома» векового хода температур на земной поверхности и «переломным» моментам формирования и деградации неоплейстоценовой мерзлоты, ледникового покрова.
Модель процессов нефтегазообразования [21] позволяет по геотемпературному критерию выполнить выделение очагов интенсивного образования нефтей из РОВ материнских отложений: пороговая геотемпература вхождения в главную зону нефтеобразования (ГЗН) и>95 °С для РОВ гумус-ово-сапропелевого типа.
Количественное определение влияния палео-климата на расчетный геотермический режим и на оценку степени реализации генерационного потенциала материнских отложений выполняется на основе анализа вариабельности результатов пяти вариантов палеотемпературных реконструкций. Вариант 1 - без учета факторов палеоклимата. Вариант 2 - учет «стандартного» векового хода температур, без учета неоплейстоценовой мерзлоты и ледников. Вариант 3 - учет «арктического» векового хода температур, без учета неоплейстоценовой мерзлоты и ледников. Вариант 4 - учет «арктического» векового хода температур, учет динамики неоплейстоценовой мерзлоты, без учета ледников. Вариант 5 - учет «арктического» векового хода температур, учет динамики неоплейстоценовой мерзлоты и ледников.
Основным критерием адекватности и предпочтительности результатов палеотемпера-турного моделирования выступает оптимальная согласованность («невязка») максимума расчетных геотемператур с «наблюденными» температурами «максимального палеотермометра» - с температурами, определенными по ОСВ. В той же степени важна оптимальная «невязка» расчетных геотемператур и с «наблюденными» пластовыми температурами. Оптимальная «невязка», принятая в классической разведочной геофизике, - это средняя квадратичная разность расчетных и наблюденных значений, равная погрешности наблюдений [22]. Эта погрешность порядка ±2 °С.
Таблица 2. Пример параметрического описания седиментационной истории и теплофизических свойств осадочной толщи, вскрытой скважиной Ростовцевская 64
Table 2. Parametrical description of history sedimentation and thermal properties of sedimentary thickness opened with the Rostov-skaya well 64
Свита, толща (стратиграфия) Suite, strata (stratigraphy) Мощность, м Thickness, m Возраст, млн лет назад Age, million years ago Время накопления, млн лет Accumulation time, million years Плотность, г/см3 Density, g/cm3 Теплопроводность, Вт/м-град Thermal conductivity, W/m-deg Температуропроводность, м2/с Thermal diffusivity, m2/s Тепловыделение, Вт/м3 Heat generation W/m3
Неоплейстоцен/Neo-Pleistocene Q-N2 - 0,00-0,015 0,015 - - - -
-500 0,015-0,02 0,005 - - - -
- 0,02-0,04 0,02 - - - -
-1000 0,04-0,05 0,01 - - - -
- 0,050-0,120 0,070 - - - -
+1000 0,120-0,130 0,010 0,92 2,25 1,2е-006 1,22е-007
- 0,130-0,177 0,047 - - - -
+500 0,177-0,182 0,005 0,92 2,25 1,2е-006 1,22е-007
300 0,1820-0,18215 0,00015 2,10 2,09 1,05е-006 1,22е-006
300 0,18215-0,1823 0,00015 2,10 1,3 7е-007 1,22е-006
-600 0,1823-0,1826 0,0003 - - - -
- 0,1826-0,5167 0,3341 - - - -
600 0,5167-0,5197 0,003 2,10 2,09 1,05е-006 1,22е-006
-600 0,5197-0,520 0,0003 - - - -
Квартер+плиоцен Q-N2/Quarter+Pliocene 266 0,520-4,1 3,58 2,04 1,29 6,5е-007 1,1е-006
N-2 -113 4,1-5,4 1,3 - - - -
Новопортовская N,-2/Novoportovskaya 50 5,4-8,4 3 2,08 1,33 7е-007 1,2е-006
Таволжанская N/Tavolzhanskaya 25 8,4-12,5 4,1 2,08 1,33 7е-007 1,2е-006
Бищеульская bsch N/Bishcheulskaya 38 12,5-14,5 2 2,08 1,33 7е-007 1,2е-006
N, -335 14,5-18,5 4 - - - -
Абросимовская N,/Abrosimovskaya 20 18,5-23,0 4,5 2,08 1,33 7е-007 1,2е-006
Туртасская tur P/Turtasskaya 50 23,0-28,0 5 2,08 1,33 7е-007 1,2е-006
Новомихайловская nvm P/Novomikhaylovskaya 35 28,0-30,0 2 2,08 1,33 7е-007 1,2е-006
Атлымская atl P3/Atlymskaya 65 30,0-34,0 4 2,08 1,33 7е-007 1,2е-006
Тавдинская tv P2/Tavdinskaya 110 34,0-42,6 8,6 2,08 1,33 7е-007 1,2е-006
Нюрольская nl P2/Nyurolskaya 55 42,6-50,4 7,8 2,08 1,33 7е-007 1,2е-006
Ирбитская P2/r/Irbitskaya 180 50,4-55,0 4,6 2,09 1,35 7е-007 1,2е-006
Серовская P, sr/Serovskaya 100 55,0-58,0 3 2,09 1,35 7е-007 1,2е-006
Тибейсалинская P, tb/Tibeysalinskaya 40 58,0-63,7 5,7 2,09 1,35 7е-007 1,2е-006
Ганькинская K2+P,gn/Gankinskaya 80 63,7-73,0 9,3 2,11 1,37 7е-007 1,25е-006
Березовская K2b/Berezovskaya 338 73,0-89,0 16 2,15 1,41 7,5е-007 1,25е-006
Кузнецовская K2kz/Kuznetsovskaya 18 89,0-92,0 3 2,18 1,43 8е-007 1,25е-006
Марресалинская K—K, mr/Marresalinskaya 466 92,0-102,0 10 2,26 1,49 8е-007 1,25е-006
Яронгская Kjar/Yarongskaya 206 102-108,5 6,5 2,39 1,6 8е-007 1,25е-006
Танопчинская K, tn/Tanopchinskaya 718 108,5-133,2 24,7 2,44 1,62 8е-007 1,25е-006
Ахская K, ah/Akhskaya 750 133,2-142,7 9,5 2,44 1,64 8е-007 1,25е-006
Баженовская J3+K, bg/Bazhenovskaya 33 142,7-149,3 6,6 2,42 1,62 8е-007 1,3е-006
Нурминская J2nr/Nurminskaya 72 149,3-161,7 12,4 2,42 1,62 8е-007 1,3е-006
Малышевская J2ml/Malyshevskaya 105 161,7-171,0 9,3 2,45 1,63 8е-007 1,3е-006
Леонтьевская J2ln/Leontyevskaya 113 171,0-173,0 2 2,47 1,65 8е-007 1,3е-006
Мощность разреза, м/Section capacity, m 3485
Примечание. Коричневой заливкой показаны времена накопления нефтематеринской баженовской свиты и ее параметрическое описание; серой - размывы палеоген-неогеновых отложений; синей - времена формирования, существования и деградации толщи мерзлоты; светло-синей - времена формирования, существования и деградации толщи ледников.
Note. Accumulation time of Bazhenovpetromaternal suites and its parametric description are filled with brown. Erosion of the Paleo-gene-Neogene deposits is filled with grey. Times of «instant» formation and «instant» degradation of the Neo-Pleistocene permafrost thickness are filled with blue. Time of existence and degradation of glaciers thickness is filled with dark blue.
В качестве второго основного критерия адекватности и предпочтительности результатов принята степень согласованности очагов интенсивной генерации УВ, выделяемых по геотемпературному критерию в материнских свитах, с установленной геологоразведкой нефтегазоносностью недр.
Важным критерием достоверности результатов палеотемпературного моделирования является согласованность расчетных значений плотности теплового потока q с данными экспериментального определения плотности теплового потока на территории исследований.
Для сопоставительной оценки вариантов моделирования для материнской свиты выполняется экспресс-расчет интегрального показателя плотности ресурсов генерированных нефтей (R, усл. ед.) по формуле [23]:
R =£ (Ut • 10Л
i=i
где Ul - расчетная геотемпература очага генерации нефти, С; tl - интервальное время действия очага -нахождения материнских отложений в ГЗН, млн лет; количество временных интервалов n определено числом интервалов геологического времени нахождения материнских отложений в ГЗН. Как следует из принятой формулы, расчетное значение плотности генерированных ресурсов (на участке скважины) напрямую зависит от времени нахождения материнской свиты в ГЗН и от геотемператур ГЗН.
Эффективность метода исследований заключается в следующем:
1. В отличие от известных систем бассейнового моделирования (ГАЛА - российская, Temis, PetroProb, PetroMod - зарубежные, реализующие палеотемпературное моделирование), используемая модель в программном обеспечении палеотемпературного моделирования TeploDia-log не требует априорных сведений о природе и величине глубинного теплового потока q из основания осадочного разреза. Глубинный тепловой поток определяется решением обратной задачи геотермии - классической обратной задачи разведочной геофизики. Решение обратной задачи выполняется в рамках параметрического описания седиментационной истории и истории теплофизических свойств только осадочной толщи, без привлечения сведений о геодинамике ниже основания осадочного разреза. В то время как известна сложность и неоднозначность определения теплового потока из основания осадочной толщи, базирующегося в системах ГАЛА, Temis, PetroProb и PetroMod на моделях рифтинга литосферы («defined rift phases») [3, 5, 24, 25].
2. Для геодинамических условий Западной Сибири, характеризующихся начиная с юрского времени, квазистационарностью глубинного
теплового q [26, 27], решение обратной задачи геотермии - определения q - в системе Tep1oDi-а1о§ выполняется однозначно. 3. В математическую модель, реализованную в Tep1oDia1og, непосредственно, в строгой математической форме, включены палеотемперату-ры из определений отражательной способности витринита, как «наблюденные». Геотемпературы из ОСВ используются аналогичным образом, как и измерения пластовых температур, полученные при испытаниях скважин. Для геотемператур по ОСВ учитывается время срабатывания «максимального палеотермометра». Никаких отдельных «калибровок» по температурам ОСВ выполнять не требуется.
Анализ влияния факторов палеоклимата на расчетный геотермический режим и оценку степени реализации генерационного потенциала баженовских отложений
В работе [4] выполнен анализ результатов пале-отектонических и палеотемпературных реконструкций осадочного разреза, включающего баже-новскую свиту, для скважины Арктическая 11. Отметим, что в указанной работе анализ выполнен без учета влияния ледниковых покровов. В настоящей статье выполнен детальный анализ влияния более полного комплекса палеоклиматических факторов (мезозойско-кайнозойский вековой ход температур на поверхности Земли, вековой ход мощностей неоплестоценовой мерзлоты и вековой ход мощностей позднечетвертичных ледниковых покровов), на примере скважины Ростовцевская 64, с последующим сводным анализом по всем трем скважинам (рис. 2).
Анализ расчетных значений плотности теплового потока q из основания осадочного разреза (табл. 3) показывает следующее. В вариантах 1, 3, 4 и 5 тепловой поток увеличивается на 2,4-2,6-3,6-2,4 мВт/м2по отношению к расчетному значению теплового потока варианта 2 -46,7 мВт/м2.
Анализ термической истории баженовской свиты (табл. 3) в разрезе скважины Ростовцевская 64 свидетельствует о том, что в варианте 1 (без учета всех факторов палеоклимата) материнская свита «пережила» самую короткую и самую «холодную» главную фазу нефтеобразования (ГФН).
В вариантах 2, 3, 4 и 5 (с учетом факторов па-леоклимата) баженовская свита имеет «богатые» термические истории ГФН. Главные фазы нефте-образования этих вариантов имеют разные значения абсолютных максимумов палеотемператур.
В варианте 4 присутствие толщи вечномер-злых пород, обладающих высокими значениями теплопроводности Я и температуропроводности a, приводит к увеличению расчетных значений плотности теплового потока q, что, в свою очередь, увеличивает расчетные геотемпературы материнских отложений.
Таблица 3. Расчетные геотемпературы баженовской свиты в разрезе скважины Ростовцевская 64 Table 3. Calculated geotemperatures of the Bazhenov suite in the section of Rostovskaya 64 well
Время, млн лет назад Time, million years ago «Стандартный» вековой ход, °С «Standard» century course, °С «Арктический» вековой ход, °С «Arctic» secular variation, °С Мощность мерзлых пород,м Capacity of permafrost rocks, m Мощность ледников, м Capacity of glaciers, m Глубина положения баженовской свиты, м Basement depth of Bazhenov suite, m Геотемпературы баженовской свиты, "С Geotemperatures of Bazhenov suite, "С
Вариант/Variant
1 2 3 4 Б
Геотемпературы Geotem-peratures Глубина положения баженовской свиты, м (с учетом ледникового покрова) Basement depth of Bazhenov suite, m (taking into account a glacial cover)
G G -4 3GG - З1У9 106 101 98 96 98 3179
0,015 -2 -1G 3GG - З1У8 106 100 98 96 9У 3179
G,G2 -3 -8 3GG 500 З1УУ 106 100 98 96 96 3679
G,G3 -4 -Б 3GG 500 З1УУ 106 100 98 96 95 3679
G,G4 -2 -6 3GG 500 З1УУ 106 100 98 96 95 3679
0,05 -1 -У 3GG 1Б00 З1У6 106 101 98 95 94 4679
0,0У -4 -4 3GG 1500 З1У4 106 101 9У 95 93 4679
G,11 -4 -Б 3GG 1500 З1У2 106 101 9У 94 92 4679
G,12 -2 -6 3GG 1500 З1У1 106 101 9У 93 91 4679
G,13 -1 -У 3GG 500 З1У1 106 101 9У 93 91 3679
С,1Б -4 -6 3GG 500 3169 106 101 96 93 91 3679
G,1УУ -6 -У 3GG 500 З16У 106 101 96 93 91 3679
G,182 -6 -У 3GG - З16У 106 101 96 93 91 3179
G,1826 -У -У 6GG - З16У 105 101 96 93 91 3167
G,2G -У -8 6GG - 3166 105 102 96 93 91 3166
G,24 -1G -9 6GG - 3163 105 102 96 93 91 3163
G,Б16У -6 -1G 6GG - 3145 104 102 94 98 95 3145
0,Б19У -Б -11 - - 3145 104 102 94 98 95 3145
G,Б2 -Б -11 - - 3145 104 102 94 98 95 3145
1,8 -3 -13 - - 3062 101 99 93 95 93 3062
3,2 -2 +Б - - 29У1 98 9У 103 105 103 2971
4,1 +3 +4 - - 2913 9У 96 102 103 101 2913
4,9 +Б +4 - - 2983 99 9У 104 106 103 2983
Б,4 +Б +4 - - 3026 100 100 105 10У 105 3026
8,4 +9 +Б - - 29У6 98 100 104 106 104 2976
1G +9 +6 - - 2966 98 101 105 106 104 2966
12,Б +1G +6 - - 2951 9У 101 104 106 104 2951
14,Б +1G +6 - - 2913 9У 102 104 106 104 2913
18,Б +11 +У - - 3248 109 113 11У 119 116 3248
23 +4 +8 - - 3228 108 108 11У 119 116 3228
28 +8 +8 - - З1У8 106 109 114 11У 115 3178
3G +1G +9 - - 3143 105 108 114 116 114 3143
34 +11 +9 - - З0У8 102 109 112 114 111 3078
ЗБ +14 +9 - - 3065 102 111 112 113 111 3065
42,6 +2G +12 - - 2968 98 114 111 112 111 2968
50 +21 +1Б - - 2916 96 112 111 113 111 2916
БG,4 +21 +1Б - - 2913 96 112 111 113 111 2913
ББ +21 +1Б - - 2УЗЗ 88 105 104 106 104 2733
Б8 +2G +16 - - 2633 85 101 101 102 101 2633
6З,У +19 +16 - - 2593 85 101 101 102 101 2593
УG +19 +16 - - 2539 83 98 99 101 99 2539
УЗ +19 +1Б - - 2513 82 9У 9У 99 9У 2513
8Б +19 +13 - - 2260 У2 89 86 8У 86 2260
89 +2G +13 - - 21УБ 69 86 83 85 83 2175
92 +2G +13 - - 21БУ 69 86 83 84 82 2157
1GG +22 +1Б - - 1У84 55 У4 У1 У2 У1 1784
1G2 +22 +1Б - - 1691 52 У1 68 У0 68 1691
108,5 +22 +1Б - - 1485 46 65 62 62 61 1485
12G +22 +16 - - 1151 35 55 51 52 51 1151
134 +21 +1Б - - У04 21 42 36 ЗУ ЗУ 704
1ЗБ +22 +1Б - - 625 19 40 34 35 34 625
142,Б +22 +1Б - - 33 1 23 16 16 16 33
Расчетный тепловой поток из основания, мВт/м2 Calculated heat flow from the basement, mW/m2 49,1 46,У 49,3 50,3 49,1
Отметим, что учет ледникового покрова (вариант 5) мало повлиял и на величину расчетного значения плотности теплового потока из основания q, и на интенсивность и продолжительность ГФН.
Сопоставление расчетных и «наблюденных» геотемператур для 3-х скважин приведено в табл. 4. Так как «наблюденные» (измеренные) температуры (включая определенные по ОСВ) имеют погрешность порядка ±2 °С, то варианты 1 и 2 решений нельзя признать приемлемыми. В этих вариантах «невязки» превышают оптимальную в 2 раза и много больше, а в варианте 1 разница с ОСВ достигает 11-12 °С.
В случае учета палеоклимата (варианты 3, 4и 5) как «невязки» для пластовых температур, так и сходимость с «максимальным палеотермометром» оптимальны и примерно равноценны. Несколько иные результаты получены для скважины Средне-Ямальская 14. Здесь оптимальным (приемлемым) можно признать, пожалуй, только вариант 4.
Расчет интегрального показателя Я (табл. 5) -экспресс-расчет плотности ресурсов генерированных ресурсов баженовских нефтей - дает максимальное или несколько большее значение для варианта 4, приемлемого по критерию «невязки». В этом варианте, помимо учета «арктического» векового хода температур на дневной поверхности, учтена динамика неоплейстоценовой мерзлоты.
Максимальные значения для скважин Средне-Ямальская 14 и Ростовцевская 64 (вариант 4: 73, 76 усл. ед.) примерно те же, что и в варианте 3 (72, 77 усл. ед.), в котором учтен только один фактор палеоклимата - вековой ход температур, но на 40-50 % больше, чем в варианте 1 (53, 51 усл. ед.), в котором не учтен ни один фактор палеоклимата.
Таблица 4. Сопоставление измеренных и расчетных геотемператур в скважинах Арктического, Средне-Ямальского и Ростовцевского месторождений Table 4. Comparison of the measured and calculated geotem-peratures in the wells of Arctic, Sredne-Yamal and Rostovtsev fields
и Вариант/Variant
E о U 3 £ CP CD d 1 2 3 4 5
Cp CD Q га Е CP CD ш ш
d - CD м nt e и U e и U e U U e U U e U U
5 m cd a p a p a p a p a p
та X CD £ Ю CD 8 з /V е/ CD и > е/ CD U > е/ CD U > е/ CD U > е/ s CD U
Ю ¥ ш ° ° га О О О О О
СРО^ и ^ ен а/ ен а/ ен а/ ен а/ ен а/
5 "О ач s ач s ач s ач s ач s
m з m з m з m з m з
га Q_ Q_ Q_ Q_ Q_
Me
Скважина Арктическая 11/Arcticheskaya 11 well
2000 100 по ОСВ 91 -9 97 -3 99 -1 102 +2 100 0
2500 120 by OSV 109 -11 115 -5 117 -3 121 +1 119 -1
3533 125 пласто- 133 +8 128 +3 126 +1 124 -1 126 +1
3560 126 in-place 133 +7 129 +3 127 +1 125 -1 126 0
Среднеквадрати-
ческое отклонение
(«невязка»), °С ±9 ±4 ±2 ±1 ±1
Mean squared error («true error»), °С
Скважина Средне-Ямальская 14/Sredne-Yamalskaya 14 well
1700 83 по ОСВ by OSV 78 -5 80 -3 80 -3 81 -2 80 -3
2200 100 95 -5 96 -4 96 -4 97 -3 96 -4
3000 120 122 +2 121 +1 122 +2 122 +2 121 +1
пласто-
846 17 вые in-place 36 +19 32 +15 28 +11 23 +6 29 +12
Среднеквадрати-
ческое отклонение
(«невязка»), °С ±10 ±8 ±6 ±4 ±6
Mean squared error («true error»), °С
Скважина Ростовцевская 64/Rostovtsevskaya 64 well
2470 75 пласто- 84 +9 80 +5 77 +2 75 0 77 +2
2650 81 вые 90 +9 85 +4 82 +1 81 0 82 +1
2660 81 in-place 90 +9 86 +5 83 +2 81 0 83 +2
2096 84 по ОСВ by OSV 76 -8 81 -3 84 0 85 +1 83 -1
2600 98 92 -6 96 -2 99 +1 101 +3 99 +1
2827 111 99 -12 103 -8 106 -5 108 -3 106 -5
Среднеквадрати-
ческое отклонение
(«невязка»), °С ±9 ±5 ±2 ±2 ±2
Mean squared error («true error»), °С
Примечание. Коричневой заливкой показаны варианты оптимальные (приемлемые) по критерию «невязки».
Note. The optimal variants (accepted) by «residual» criterion are shown in brown filling.
Достоверность результатов палеотемпера-турного моделирования, выполненного на Арктической, Средне-Ямальской и Ростовцевской площадях, подтверждается хорошей согласованностью по-
Примечание. Вариант 1 - без учета факторов палеоклимата. Вариант 2 - учет «стандартного» векового хода температур [28, 29], без учета неоплейстоценовой мерзлоты и ледников. Вариант 3 - учет «арктического» векового хода температур [4], без учета неоплейстоценовой мерзлоты и ледников. Вариант 4 - учет «арктического» векового хода температур и динамики неоплейстоценовой мерзлоты. Вариант 5 - учет «арктического» векового хода температур, динамики неоплейстоценовой мерзлоты и льда. Коричневой заливкой показаны температуры главной фазы нефтеобразования (ГФН), темно-коричневой заливкой - палеотемпературный максимум ГФН. Серой заливкой обозначены времена размыва палеоген-неогеновых отложений.
Note. Variant 1 - excludingpaleoclimate factors. Variant 2 - considering «standard» secular variation of temperatures [28, 29], without Neo-Pleistocene permafrost and glaciers. Variant 3 -considering «arctic» secular variation of temperatures [4], without Neo-Pleistocene permafrost and glaciers. Variant 4 - considering «arctic» secular variation of temperatures and dynamics of Neo-Pleistocene permafrost. Variant 5 - considering «arctic» secular variation of temperatures and dynamics of Neo-Pleistocene permafrost and the ice. Brown shaded areas indicate temperatures of the major oil generation zone (GFN), dark-brown colour shading indicates paleotemperature maxima GFN, grey shading shows the erosion times of Paleogene-Neogene sediments.
лученных расчетных значений плотности теплового потока (50-51-58 мВт/м2) с экспериментальными определениями плотности теплового потока для п-ва Ямал: 50-55 мВт/м2 [27], 49-56 мВт/м2 [30].
Как сообщалось выше, баженовская свита является основным источником формирования залежей УВ в ловушках верхнеюрского и мелового НГК, а также возможной сланцевой формацией. В этой связи важно оценить согласованность очагов интенсивной генерации баженовских нефтей, выделенных по геотемпературному критерию в разрезах скважин, с результатами испытаний верхнеюрских и нижнемеловых пластов (табл. 1). Но, к сожалению, в этих скважинах верхнеюрский НГК не испытывался, а нижнемеловые пласты испытаны только в скважине Ростовцевская 64. Результаты испытаний пластов ахской свиты в целом положительные и подтверждают очаг генерации нефти в разрезе скважины Ростовцевская 64. А для интенсификации притока при испытании непосредственно баженовской свиты вероятно требуются специальные работы - возможно, горизонтальный ствол и гидроразрыв пласта.
Сопоставление и обсуждение
результатов исследований
Сопоставление измеренных и расчетных геотемператур (табл. 4) позволяет заключить, что по классическому геофизическому критерию - критерию «невязки» - результаты только вариантов 3, 4 и 5 оптимальны (приемлемы) и равноценны. Учет «арктического» векового хода температур, неоплейстоценовой мерзлоты и ледникового покрова позволило достаточно корректно восстановить термическую историю баженовских отложений. Учет палеоклимата обуславливает увеличение расчетного палеотемпературного максимума в истории материнских отложений на площадях п-ва Ямал, на 10-13 С.
Полученные результаты по арктическому району Западной Сибири вполне согласуются с характером ранее полученных оценок существенного влияния мезозойско-кайнозойского климата, включая резкое похолодание в конце плиоцена, на геотермический режим осадочного чехла юго-востока Западной Сибири [31].
Вариант 4, учитывающий «арктический» вековой ход температур на земной поверхности и неоплейстоценовую мерзлоту - факторы палеоклима-та, оптимальный по сходимости («невязке») измеренных и расчетных геотемператур, представляет наиболее «богатую» термическую историю материнских баженовских отложений, а следовательно, обеспечивает наибольшую расчетную плотность ресурсов генерированных нефтей (табл. 5).
Учет ледникового покрова (вариант 5) мало повлиял (по отношению к варианту 4) и на величину расчетного значения плотности теплового потока из основания q, и на интенсивность и продолжительность ГФН, и на экспресс-оценку генерирован-
Таблица 5. Расчет интегрального показателя R, дающего экспресс-оценку плотности ресурсов генерированных баженовских нефтей Table 5. Calculation of the integrated indicator R giving the express assessment of density of resources of the generated Bazhenov oil
Вариант палеотемпературного моделирования Variant of paleo temperature modeling Экспресс-расчет ресурсов (R), усл. ед Express calculation of resources (R), cu. Период работы палеоочага млн лет назад Period of the work of the paleohearth, million years ago Время работы палеоочага, млн лет Lifetime of the paleohearth, million years Расчетная плотность теплового потока из основания осадочного разреза, мВт/м2 Settlement density of a thermal stream from the basis of a sedimentary section, mW/m2 Максимальные геотемпературы палеоочага генерации нефти, °С Maximum geotemperatures of the paleohearth oil generation, °С
Скважина Арктическая 11/Arcticheskaya 11 well
Вариант 1 Variant 1 52 50,4-0,0 50,4 56 119
Вариант 2 Variant 2 97 92,0-0,0 92,0 54 125
Вариант 3 Variant 3 69 85,0-3,2; 0,52-0,0 82,3 56 128
Вариант 4 Variant 4 102 92,0-0,2; 0,12-0,0 91,9 58 132
Вариант 5 Variant 5 98 89,0-0,5167; 0,070-0,0 88,6 58 130
Скважина Средне-Ямальская 14/Sredne-Yamalskaya 14 well
Вариант 1 Variant 1 53 50,4-0,0 50,4 54 117
Вариант 2 Variant 2 77 73,0-0,0 73,0 50 117
Вариант 3 Variant 3 72 70,0-3,2 66,8 51 118
Вариант 4 Variant 4 73 70,0-3,2 66,8 51 118
Вариант 5 Variant 5 71 70,0-3,2 66,8 51 117
Скважина Ростовцевская 64/Rostovtsevskaya 64 well
Вариант 1 Variant 1 51 50,4-0,0 50,4 49 109
Вариант 2 Variant 2 77 73,0-0,0 73,0 47 113
Вариант 3 Variant 3 77 73,0-3,2; 0,24-0,0 70,0 49 117
Вариант 4 Variant 4 76 73,0-0,5167; 0,07-0,0 72,6 50 119
Вариант 5 Variant 5 77 73,0-3,2; 0,52-0,5167; 0,04-0,0 69,8 49 116
Примечание. Коричневой заливкой показаныi варианты, оптимальные (приемлемые) по критерию «невязки».
Note. The optimal variants (accepted) by «residual» criterion are showed in brown filling.
ных ресурсов УВ (табл. 5). Такой результат объясняется тем, что «возмущения» тепловой модели (температур, теплофизических свойств) в верхней части разреза «доходят» в полной мере на глубины порядка 3000 м не ранее чем через 0,2-0,3 млн лет. Такое «запаздывание» было установлено нами и ранее [19, 32]. Очевидно, что это «запаздывание» в большей степени касается влияния ледниковых покровов, возраст которых принят 0,18 млн лет.
Вместе с тем можно отметить (табл. 3), что включение в модель динамики вечномерзлых пород (возраст 0,52 млн лет) закрывает «холодное» окно геотемператур баженовской свиты в период 1,80-0,52 млн лет назад, присутствующее в варианте 3 и обусловленное резким похолоданием на земной поверхности в этот период. Но, в свою очередь, с началом формирования вечномерзлых пород открывается «холодное» окно начиная с 0,24 млн лет назад (вариант 4), увеличиваясь в продолжительности при учете ледникового покрова (вариант 5).
Если принять во внимание полученные нами результаты, то можно констатировать несущественное влияние ледникового покрова (на землях п-ва Ямал) на термический режим материнских баженовских отложений, имеющий ключевое значение в процессах генерации УВ. Однако остается невыясненной роль ледников как фактора дополнительного давления на осадочный чехол, влияющего на генерацию и формирование залежей УВ.
В рамках принятой нами модели среды - твердое тело, горизонтальная слоистость - вполне адекватно описываются процессы теплопереноса для «плитной» осадочной толщи, при условии существенного превалирования кондуктивного механизма переноса тепла [33]. Названными условиями характеризуется Западно-Сибирская плита. Здесь магматические и гидротермальные процессы незначительны, отсутствуют или затухли 160-190 млн лет назад. Но в части механизма передачи давления для твердого тела дело обстоит не так. Дополнительно возникающее давление передается почти мгновенно, со скоростью сейсмической волны [34]. Для осадочного разреза п-ва Ямал дополнительное давление на баженовскую свиту от слоя льда мощностью 500-1500 м составит добавку от 7 до 20 %. Вместе с тем, поскольку осадочный разрез преимущественно горизонтально слоистый, каких-то аномальных (отличных от ли-тостатического) давлений, напряжений и смещений среды в гравитационном поле не будет.
Выводы
1. Рассмотрение проблемы совершенствования технологий прогнозирования нефтегазоносно-сти в арктических районах Западной Сибири, выполняемого объемно-генетическим методом, показало актуальность количественной оценки и учета роли мезозойско-кайнозойского климата, в частности вечномерзлых пород и ледников в позднечетвертичное время.
2. Выделены три основных палеоклиматических фактора Севера Западной Сибири, вероятно влияющих на реконструируемый геотермический режим осадочного разреза, включая неф-тематеринские отложения: 1-й фактор - вековой ход температур на поверхности Земли;
2-й фактор - формирование и деградация неоплейстоценовой толщи вечномерзлых пород;
3-й фактор - позднечетвертичные ледниковые покровы.
3. Сформулирована задача дальнейших исследований роли 1-го и 2-го факторов палеоклимата, впервые с комплексным учетом 3-го фактора, в реконструкциях геотермического режима неф-тематеринских баженовских отложений и в расчете степени реализации их генерационного потенциала.
4. Реализована методика исследований, которая основывается: 1) на палеотемпературном моделировании, наиболее полно учитывающем параметры тектоно-седиментационной истории и истории теплофизических свойств осадочной толщи, измеренные пластовые температуры и пале-отемпературы, определенные по данным ОСВ; 2) многовариантности палеотектонических и пале-отемпературных реконструкций и на анализе вариабельности результатов; 3) оценке соответствия результатов общепринятым критериям оптимальности решения обратной задачи геофизики; 4) согласованности расчетных значений плотности теплового потока с данными экспериментального определения плотности теплового потока на территории исследований.
5. На представительных мезозойского-кайнозой-ских разрезах месторождений п-ва Ямал (скважины Арктическая 11, Средне-Ямальская 14, Ростовцевская 64) установлено, что неучет «арктического» векового хода температур на поверхности Земли и толщи неоплейстоценовой мерзлоты не позволяет адекватно восстановить термическую историю баженовских отложений.
6. Достоверность результатов палеотемператур-ного моделирования уверенно контролируется сопоставлением с экспериментальными данными о тепловом потоке на территории исследований.
7. При определении ресурсов УВ объёмно-генетическим методом на землях арктического региона Западной Сибири рекомендуется применять «арктический» вековой ход температур и учитывать динамику толщи неоплейстоценовой мерзлоты мощностью 300-600 м. В случае неучета толщ вечной мерзлоты и палеоклимати-ческого хода температур расчетные ресурсы УВ могут быть занижены до 40-50 %.
8. Полученные нами первые результаты по оценке роли позднечетвертичных ледниковых покровов (в районе п-ва Ямал) позволили отметить несущественное влияние ледникового по-
крова на термический режим материнских ба-женовских отложений, имеющий ключевое значение в процессах генерации УВ.
Заключение
Объектом проведенных исследований был геотермический режим нефтематеринской баженов-ской свиты, которая является основным источником формирования залежей УВ в ловушках верхнеюрского и мелового НГК, а также возможной сланцевой формацией арктического региона Западной Сибири. Вместе с тем представляет несомненный интерес для исследований геотермический режим китербютской свиты (Jjkt), являющейся источником формирования залежей УВ в ловушках нижнеюрского и, возможно, доюрского НГК. Глинистая толща китербютская, обладающая нефте-материнским потенциалом, формируется во времена бореальных трансгрессий в ранней юре [35].
Как уже отмечалось, добавка к литостатиче-скому давлению за счет ледниковых покровов может составить 7-20 %, чем, возможно, нельзя пренебрегать. Кроме того, материал горных пород обладает диффузионной и дислокационной ползучестью [36]. При небольших глубинах преобладает диффузионная ползучесть. При такой ползучести среда ведет себя как ньютоновская несжимаемая сильно вязкая жидкость: вязкость порядка 1019—1021 Пуаз. В этом случае смещения точек среды могут составить несколько см/год. Учитывая скорость диффузионного смещения, «волна» аномальных диффузионных деформаций, обусловленная ледниковыми покровами, достигает глубин залегания баженовской свиты не позже, чем через
0.3.млн лет. Таким образом, дополнительное давление ледников по механизму диффузионной ползучести вероятно и может влиять на генерацию УВ.
Возможность выявления перерывов осадкона-копления и денудации методом геотермии, как
СПИСОК ЛИТЕРАТУРЫ
1. Конторович А.Э. Проблемы реиндустриализации нефтегазового комплекса России // Нефтяное хозяйство. - 2016. - № 3. - С. 14-15.
2. Tissot В.Р. Preliminary Data on the Mechanisms and Kinetics of the Formation of Petroleum in Sediments. Computer Simulation of a Reaction Flowsheet // Oil & Gas Science and Technology -Rev. IFP. - 2003. - V. 58. - № 2. - P. 183-202.
3. Historical-geological modeling of hydrocarbon generation in the Mesozoic-Cenozoic sedimentary basin of the Kara sea (basin modeling) / A.E. Kontorovich, L.M. Burshtein, N.A. Malyshev, P.I. Safronov, S.A. Gusfkov, S.V. Ershov, V.A. Kazanenkov, N.S. Kim, V.A. Kontorovich, E.A. Kostyreva, V.N. Melenevsky, V.R., Livshits A.A. Polyakov, M.B. Skvortsov // Russian Geology and Geophysics. - 2013. - V. 54. - № 8. - P. 1179-1226.
4. Искоркина А.А. Палеоклиматические факторы реконструкции термической истории нефтематеринской баженовской свиты арктического региона Западной Сибири // Известия Томского политехнического университета. Инжиниринг георесурсов. - 2016. - Т. 327. - № 8. - C. 59-73.
5. Influence of Surface Temperatures on Source Rock Maturity: an Example from the Russian Artic / S. Nelskamp, T. Donders, J.-D. van Wess, O. Abbink // ROGTEC. - 2014. - № 18. - P. 26-35.
методом, иным по физическим основам, наряду с геологическими методами и сейсморазведкой, существенно повышает достоверность знаний об истории геологического развития регионов Западной Сибири. Такая возможность особенно актуальна для северных районов Западно-Сибирской НГП, где по первым нашим оценкам, приведенным выше, размывы толщ палеоген-неогеновых пород достигают 400-600 м. В то время как для южного региона Западной Сибири размывы либо вообще не фиксируются, либо оцениваются в пределах 50-100 м.
Полученные первые результаты исследований по выявлению и учету параметров перерывов осад-конакопления и денудаций в комплексе с уникальными палеоклиматическими характеристиками Севера Сибири (вековой ход температур на земной поверхности, динамика плиоцен-четвертичных толщ вечномерзлых пород и динамика четвертичных ледниковых покровов) позволили сделать выводы и дать первые рекомендации. А именно, как адекватно восстановить термическую историю нефтематеринских баженовских отложений, а следовательно, повысить достоверность подсчета ресурсов УВ объемно-генетическим методом.
Таким образом, на очереди стоят следующие исследования для условий арктических районов Западной Сибири: 1) оценка роли палеоклимати-ческих факторов реконструкции термической истории китербютской свиты; 2) разработка формализованных сценариев и критериев выявления и количественной оценки параметров перерывов ос-адконакопления и размывов по данным геотермии; 3) аргументированная оценка влияние роли ледников как фактора дополнительного давления на осадочный чехол, влияющего на генерацию и формирование залежей УВ.
Статья подготовлена при частичной финансовой поддержке РФФИ в рамках научного проекта № 16-35-00080 мол_а.
6. Баулин В.В. Многолетнемерзлые породы нефтегазоносных районов СССР. - М.: Недра, 1985. - 176 с.
7. Глобальные изменения климата и природной среды позднего кайнозоя в Сибири / отв. ред. А.П. Деревянко. - Новосибирск: Изд-во СО РАН, 2008. - 511 с.
8. Палеогеография Западно-Сибирской равнины в максимум поздне-зырянского оледенения / С.А. Архипов, В.И. Астахов, И.А. Волков, В.С. Волкова, В.А. Панычев. - Новосибирск: Наука, 1980. - 103 с.
9. Томирдиаро С.В. Оледенения Арктического бассейна в плейстоцене и его связь с наземным и подземным оледенением суши // Гляциологические исследования в полярных странах. -Л.: Гидрометеоиздат, 1970. - С. 172-182.
10. Лаухин С.А., Пушкарь B.C., Черепанова М.В. Современное состояние реконструкций природной среды на севере Сибири в каргинское время (поздний плейстоцен) // Бюллетень московского общества испытателей природы. Отдел геологический. -2012. - Т. 87. - № 6. - С. 37-48.
11. Четвертичные оледенения Западной Сибири и других областей Северного полушария / под ред. С.А. Архипова, А.А. Величко. -Новосибирск: Наука, 1981. - 201 с.
12. Артюшков Е.В. Четвертичные оледенения и трансгрессии в Западной Сибири // Изв. АН СССР, сер. геол. - 1969. - № 7. - С. 98-114.
13. Баду Ю.Б. Влияние газоносных структур на мощность криогенной толщи Ямала // Криосфера Земли. - 2014. - Т. XVIII. -№3. - С. 11-22.
14. Районирование баженовской свиты и клиноформ неокома по плотности ресурсов сланцевой и первично-аккумулированной нефти (на примере Нюрольской мегавпадины) / В.И. Исаев, Г.А. Лобова, А.К. Мазуров, А.Н. Фомин, В.И. Старостенко // Геофизический журнал. - 2016. - Т. 38. - № 3. - С. 29-51.
15. Стратиграфия нефтегазоносных бассейнов Сибири. Кн. 9: Кайнозой Западной Сибири. - Новосибирск: Изд-во СО РАН, 2002.- 246 с.
16. Tectonic structure and history of evolution of the West Siberian geosyneclise in the Mesozoic and Cenozoic / V.A. Kontorovich, S.Yu. Belyaev, A.E. Kontorovich, V.O. Krasavchikov, A.A. Kontorovich, O.I. Suprunenko // Russian Geology and Geophysics. -2001. - V. 42. - P. 1832-1845.
17. История тектонического развития арктических районов Западно-Сибирской геосинеклизы в кайнозойское время / С.Ю. Беляев, С.А. Гуськов, В.С. Волкова, А.В. Истомин // Интерэкспо Гео-Сибирь. - 2013. - Т. 2. - № 1. URL: http://cyber-leninka.ru/article/n/istoriya-tektonicheskogo-razvitiya-ark-ticheskih-rayonov-zapadno-sibirskoy-geosineklizy-v-kaynozoy-skoe-vremya (дата обращения 26.11.2016).
18. Volkova V.S. Geologic stages of the paleogene and neogene evolution of the Arctic shelf in the Obi region (West Siberia) // Russian Geology and Geophysics. - 2014. - V. 55. - № 4. - P. 619-633.
19. Isaev V.I., Fomin A.N. Loki of generation of bazhenov- and togur-type oils in the southern NyuroFka megadepression // Russian Geology and Geophysics. - 2006. - V. 47. - № 6. - P. 734-745.
20. Estimation of the Oil-and-Gas Potential of Sedimentary Depression in the Far East and West Siberia Based on Gravimetry and Geothermy Data / R.Yu. Gulenok, V.I. Isaev, V.Yu. Kosygin, G.A. Lobova, V.I. Starostenko // Russian Journal of Pacific Geology. - 2011. - V. 5. - № 4. - P. 273-287.
21. Модель катагенеза органического вещества (на примере баженовской свиты) / Л.М. Бурштейн, Л.В. Жидкова, А.Э. Конто-рович, В.Н. Меленевский // Геология и геофизика. - 1997. -Т. 38. - № 6. - С. 1070-1078.
22. Strakhov V.N., Golizdra G.Ya., Starostenko V.I. Theory and practice of interpreting potential fields: Evolution in the 20th century // Izvestiya - Physics of the Solid Earth. - 2000. - V. 36. -№ 9. - P. 742-762.
23. Isaev V.I., Lobova G.A., Osipova E.N. The oil and gas contents of the Lower Jurassic and Achimovka reservoirs of the Nyurol'ka megadepression // Russian Geology and Geophysics. - 2014. -V. 55.- P. 1418-1428.
24. Galushkin Y.I., Sitar K.A., Kunitsyna A.V. Numerical modeling of the organic matter transformation in the sedimentary rocks of the northeastern Sakhalin shelf // Oceanology. - 2011. - V. 51. -№ 3. - P. 491-501.
25. Razvozzhaeva E.P., Kirillova G.L., Prokhorova P.N. Comparative analysis of fragments of the Mesozoic East Asian continental margin: the Kyndal (Bureya Basin, Russia) and Suibin (Sanjiang Basin, China) troughs // Russian Journal of Pacific Geology. -2014. - V. 8. - №6. - P. 404-422.
26. Эволюция температурного поля осадочного чехла Западно-Сибирской плиты / А.Д. Дучков, Ю.И. Галушкин, Л.В. Смирнов, Л.С. Соколова // Геология и геофизика. - 1990. - № 10. -С. 51-60.
27. Kurchikov A.R. The geothermal regime of hydrocarbon pools in West Siberia // Russian Geology and Geophysics. - 2001. -V. 42. - № 11-12. - P. 678-689.
28. Лопатин Н.В. Концепция нефтегазовых генерационно-акку-муляционных систем как интегрирующее начало в обосновании поисково-разведочных работ // Геоинформатика. -2006.- №3. - С. 101-120.
29. Галушкин Ю.И. Моделирование осадочных бассейнов и оценка их нефтегазоносности. - М.: Научный Мир, 2007. - 456 с.
30. Геотермия арктических морей / М.Д. Хуторской, В.Р. Ахмед-зянов, А.В. Ермаков, Ю.Г. Леонов, Л.В. Подгорных, Б.Г. Поляк, Е.А. Сухих, Л.А. Цыбуля. - М.: ГЕОС, 2013. - 232 с.
31. Палеоклиматические факторы реконструкции термической истории баженовской и тогурской свит юго-востока Западной Сибири / В.И. Исаев, А.А. Искоркина, Г.А. Лобова, А.Н. Фомин // Геофизический журнал. - 2016. - Т. 38. - № 4. - С. 3-25.
32. Нефтегазоносность нижнемеловых резервуаров Нюрольской мегавпадины / Е.Н. Осипова, Г.А. Лобова, В.И. Исаев, В.И. Старостенко // Известия Томского политехнического университета. - 2015. - Т. 326. - № 1. - С. 14-33.
33. Хуторской М.Д. Введение в геотермию. - М.: Изд-во РУДН, 1996. - 156 с.
34. Nadai A. Theory of Flow and Fracture of Solids. V. 2. - New York: McGraw-Hill Book Company, 1963. - 705 p.
35. Богоявленская В.И., Полякова И.Д. Перспективы нефтегазо-носности больших глубин Южно-Карского региона // Арктика: экология и экономика. - 2012. - № 3 (7). - С. 92-103.
36. Cserepes L., Yuen D.A., Schroeder B.A. Effect of the mid-mantle viscosity and phase-transition structure on 3D mantle convection // Physics of the Earth and Planetary Interiors. - 2000. - V. 118. -№1-2. - P. 125-148.
Поступила 28.11.2016 г.
Информация об авторах
Исаев В.И., доктор геолого-минералогических наук, профессор кафедры геофизики Института природных ресурсов Национального исследовательского Томского политехнического университета.
Искоркина А.А., аспирант кафедры геофизики Института природных ресурсов Национального исследовательского Томского политехнического университета.
Косыгин В.Ю., доктор геолого-минералогических наук, ведущий научный сотрудник Вычислительного центра ДВО РАН.
Лобова Г.А., доктор геолого-минералогических наук, доцент кафедры геофизики Института природных ресурсов Национального исследовательского Томского политехнического университета.
Осипова Е.Н., кандидат геолого-минералогических наук, ассистент кафедры геофизики Института природных ресурсов Национального исследовательского Томского политехнического университета.
Фомин А.Н., доктор геолого-минералогических наук, заведующий лабораторией геохимии нефти и газа Института нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН.
UDC 553.98(571)
INTEGRATED ASSESSMENT OF PALEOCLIMATE FACTORS OF RECONSTRUCTING THERMAL HISTORY OF PETROMATERNAL BAZHENOV SUITE IN ARCTIC REGIONS OF WESTERN SIBERIA
Valeriy I. Isaev1,
isaevvi@tpu.ru
Albina A. Iskorkina1,
iskorkina.a@mail.ru
Vladimir Yu. Kosygin2,
kosyginv@inbox.ru
Galina A. Lobova1,
lobovaga@tpu.ru
Elizaveta N. Osipova1,
osipovaen@list.ru
Aleksandr N. Fomin3,
fominan@ipgg.sbras.ru
1 National Research Tomsk Polytechnic University, 30, Lenin Avenue, Tomsk, 634050, Russia.
2 Computer center of the FEB RAS,
65, Kym St. Yu Chen street, Khabarovsk, 680000, Russia.
3 Institute of Petroleum Geology and Geophysics named after A.A. Trofimuk SB RAS, 3, Ac. Koptueg Avenue, Novosibirsk, 630090, Russia.
The relevance. The Arctic regions of Western Siberia have the unique paleoclimate features, the inversion tectonic history in Paleogene-Neogene and the considerable variations of geochemical characteristics of petromaternal deposits. This causes the necessity to improve schemes and parameters of quantitative assessment of hydrocarbon resources by volume and genetic method based on paleorecon-struction of geotemperature mode of petromaternal deposits.
The main aim of research is to estimate the impact of paleoclimate factors - temperature secular variation of the Earth surface and Neo-Pleistocene permafrost thickness, glacial sheets, in a complex with a factor of erosive processes, on the settlement geothermal mode of petromaternal of Bazhenov suite.
Object of researches: Bazhenov deposits of Mesozoic and Cenozoic sections, opened deep wells on Arctic, Sredne-Yamal and Rostov areas (Yamal Peninsula).
The method of research is based on original computer paleotemperature modeling considering the parameters of sedimentation history and the history of thermal properties of sedimentary thickness, including permafrost rocks and glaciers, without aprioristic data on sizes and the nature of a deep heat flow.
Research results on identification and registration of denudation parameters in conjunction with the assessment of the role and taking into account the characteristics of paleoclimatic features allowed making the conclusions and recommendations. It is recommended to apply «Arctic» secular variation of temperature and take into account the dynamics of Neo-Pleistocene permafrost strata of about 300-600 meters when determining the hydrocarbon resources by the volumetric-genetic method in the territory of Arctic region. Calculation hydrocarbon resources may be understated by 40-50 % in case if not taking into account permafrost and paleoclimatic secular variation of temperatures. The authors have stated the tasks of forthcoming researches, including the estimate of the role of paleoclimatic factors in thermal history of Kiterbyut source formation and elaboration of formalized scenarios for identifying and estimating the interruptions in sedimentation and erosion according to the geothermic data.
Key words:
Paleoclimate, inversion tectonic, geotemperature mode, petromaternal Bazhenov deposits, Yamal Peninsula.
The paper was partially financially supported by the RFBR within the scientific project no. 16-35-00080 Mou_a.
REFERENCES
1. Kontorovich A.E. Problems of reindustrialization of oil and gas complex of Russia. Neftyanoe xozyaystvo - Oil Industry, 2016, no. 3, pp. 14-15. In Rus.
2. Tissot B. Preliminary Data on the Mechanisms and Kinetics of the Formation of Petroleum in Sediments. Computer Simulation of a
Reaction Flowsheet. Oil & Gas Science and Technology -Rev.IFP, 2003, vol. 58, no. 2, pp. 183-202. 3. Kontorovich A.E., Burshtein L.M., Malyshev N.A., Safro-nov P.I., Gusfkov S.A., Ershov S.V., Kazanenkov V.A., Kim N.S., Kontorovich V.A., Kostyreva E.A., Melenevskiy V.N., Livshits V.R., Polyakov A.A., Skvortsov M.B. Historical-geologi-
cal modeling of hydrocarbon generation in the Mesozoic-Cenozoic sedimentary basin of the Kara sea (basin modeling). Russian Geology and Geophysics, 2013, vol. 54, no. 8, pp. 1179-1226.
4. Iskorkina A.A. Paleoclimate factors of reconstruction of thermal history of the petromaternal Bazhenov shale of the Arctic region of Western Siberia. Bulletin of the Tomsk Polytechnic University. Geo Assets Engineering, 2016, vol. 327, no. 8, pp. 59-73. In Rus.
5. Nelskamp S., Donders T., van Wess J.-D., Abbink O. Influence of Surface Temperatures on Source Rock Maturity: an Example from the Russian Artic. ROGTEC, 2014, no. 18, pp. 26-35.
6. Baulin V.V. Mnogoletnemerzlye porody neftegazonosnykh ray-onov SSSR [Permofrost rock of oil-and-gas regions of the USSR]. Moscow, Nedra Publ., 1985. 176 p.
7. Derevyanko A.P. Globalnye izmeneniya klimata iprirodnoy sredy pozdnego kaynozoya v Sibiri [Global changes of climate and environment of the late Cenozoic in Siberia]. Novosibirsk, SO RAN Publ., 2008. 511 p.
8. Arkhipov S.A., Astakhov V.I., Volkov I.A., Volkova V.S., Pa-nychev V.A. Paleogeografiya Zapadno-Sibirskoy ravniny v maksi-mum pozdnezyryanskogo oledeneniya [Paleogeography of the West Siberian Plain in a maximum glaciation of the late Zyrian]. Novosibirsk, Nauka Publ., 1980. 103 p.
9. Tomirdiaro S.V. Oledeneniya Arkticheskogo basseyna v pleysto-cene i ego svyaz s nazemnym i podzemnym oledeneniem sushi [Freezing of the Arctic basin in a Pleistocene and its relation to the land and underground freezing of land]. Glyatsiologicheskie issledovaniya v polyarnykh stranakh [Glaciological research in the Polar Regions]. Leningrad, Gidrometeoizdat Publ., 1970. pp. 172-182.
10. Laukhin S.A., Pushkar B.S., Cherepanova M.V. Current state of reconstruction of the environment in the north of Siberia in Kar-ginsk time (late Pleistocene). Bulletin of the Moscow society of testers of the nature. Geological department, 2012, vol. 87, no. 6, pp. 37-48. In Rus.
11. Arkhipova S.A., Velichko A.A. Chetvertichnye oledeneniya Za-padnoy Sibiri i drugikh oblastey Severnogo polushariya [Quaternary glaciation of Western Siberia and other regions of the Northern Hemisphere]. Novosibirsk, Nauka Publ., 1981. 201 p.
12. Artyushkov E.V. Chetvertichnye oledeneniya i transgressii v Za-padnoy Sibiri [Quaternary glaciations and transgressions in Western Siberia]. Izvestiya AN SSSR, ser. Geologiya, 1969, no. 7, pp. 98-114.
13. Badu Yu.B. Influence of gas-bearing structures on the capacity of Yamal cryogenic strata. Cryosphere of Earth, 2014, vol. XVIII, no. 3, pp. 11-22. In Rus.
14. Isaev V.I., Lobova G.A., Mazurov A.K., Fomin A.N., Starosten-ko V.I. Zoning of the Bazhenov formation and clinoforms of the Neocom by density of resources of shale and primary accumulated oil (for example, Nurolfka megahollow). Geofizicheskiy zhurnal -Geophysical journal, 2016, vol. 38, no. 3, pp. 29-51. In Rus.
15. Stratigrafiya neftegazonosnykh basseynov Sibiri. Kn. 9: Kayno-zoy Zapadnoy Sibiri [Stratigraphy of oil and gas bearing basins of Siberia. B. 9: The Cenozoic of Western Siberia]. Novosibirsk, SO RAN Publ., 2002. 246 p.
16. Kontorovich V.A., Belyaev S.Yu., Kontorovich A.E., Krasavchi-kov V.O., Kontorovich A.A., Suprunenko O.I. Tectonic structure and history of evolution of the West Siberian geosyneclise in the Mesozoic and Cenozoic. Russian Geology and Geophysics, 2001, vol. 42, pp. 1832-1845.
17. Belyaev S.Yu., Guskov S.A., Volkova V.S., Istomin A.V. History of tectonic development of Arctic districts of the West Siberian geosyneclise in Cenozoic. Interekspo Geo-Siberia, 2013, vol. 2, no. 1. In Rus. Available at: http://cyberleninka.ru/article/n/is-toriya-tektonicheskogo-razvitiya-arkticheskih-rayonov-zapadno-sibirskoy-geosineklizy-v-kaynozoyskoe-vremya (accessed: 26 November 2016).
18. Volkova V.S. Geologic stages of the Paleogene and Neogene evolution of the Arctic shelf in the Obf region (West Siberia). Russian Geology and Geophysics, 2014, vol. 55, no. 4, pp. 619-633.
19. Isaev V.I., Fomin A.N. Loki of generation of bazhenov- and togur-type oils in the southern Nyurol'ka megadepression. Russian Geology and Geophysics, 2006, vol. 47, no. 6, pp. 734-745.
20. Gulenok R.Yu., Isaev V.I., Kosygin V.Yu., Lobova G.A., Staro-stenko V.I. Estimation of the Oil-and-Gas Potential of Sedimentary Depression in the Far East and West Siberia Based on Gravime-try and Geothermy Data. Russian Journal of Pacific Geology, 2011, vol. 5, no. 4, pp. 273-287.
21. Burshteyn L.M., Zhidkova L.V., Kontorovich A.E., Melenev-skiy V.N. Model katageneza organicheskogo veshchestva (na pri-mere bazhenovskoy svity) [Model of katagenesis of organic matter (by the example of the Bazhenov formation)]. Russian Geology and Geophysics, 1997, vol. 38, no. 6, pp. 1070-1078.
22. Strakhov V.N., Golizdra G.Ya., Starostenko V.I. Theory and practice of interpreting potential fields: Evolution in the 20th century. Izvestiya - Physics of the Solid Earth, 2000, vol. 36, no. 9, pp. 742-762.
23. Isaev V.I., Lobova G.A., Osipova E.N. The oil and gas contents of the Lower Jurassic and Achimovka reservoirs of the Nyurolfka megadepression. Russian Geology and Geophysics, 2014, vol. 55, pp. 1418-1428.
24. Galushkin Y.I., Sitar K.A., Kunitsyna A.V. Numerical modeling of the organic matter transformation in the sedimentary rocks of the northeastern Sakhalin shelf. Oceanology, 2011, vol. 51, no. 3, pp. 491-501.
25. Razvozzhaeva E.P., Kirillova G.L., Prokhorova P.N. Comparative analysis of fragments of the Mesozoic East Asian continental margin: the Kyndal (Bureya Basin, Russia) and Suibin (Sanjiang Basin, China) troughs. Russian Journal of Pacific Geology, 2014, vol. 8, no. 6, pp. 404-422.
26. Duchkov A.D., Galushkin Yu.I., Smirnov L.V., Sokolova L.S. Evoliutsiya temperaturnogo polya osadochnogo chekhla Zapadno-Sibirskoy plity [Evolution of the temperature field of the sedimentary cover of the West Siberian plate]. Russian Geology and Geophysics, 1990, vol. 10, pp. 51-60.
27. Kurchikov A. R. The geothermal regime of hydrocarbon pools in West Siberia. Russian Geology and Geophysics, 2001, vol. 42, no. 11-12, pp. 678-689.
28. Lopatin N.V. Concept of oil and gas generative and accumulative systems as the integrating beginning in justification of exploration. Geoinformatika, 2006, no. 3, pp. 101-120. In Rus.
29. Galushkin Yu.I. Modelirovanie osadochnykh basseynov i otsenka ikh neftegazonosnosti [Modeling decantation basins and assessment of their petroleum potential]. Moscow, Nauchny mir Publ., 2007. 456 p.
30. Khutorskoy M.D., Akhmedzyanov V.R., Ermakov A.V., Leo-nov Yu.G., Podgornykh L.V., Poliak B.G., Sukhikh E.A., Cybu-lia L.A. Geotermiya arkticheskikh morey [Geothermic of the Arctic seas]. Ed. by Yu.G. Leonov. Moscow, GEOS Publ., 2013. 232 p.
31. Isaev V.I., Iskorkina A.A., Lobova G.A., Fomin A.N. Paleoclima-tic factors of reconstruction of thermal history of the Bazhenov and Togur suites in the southeast of Western Siberia. Geofizicheskiy zhurnal - Geophysical journal, 2016, vol. 38, no. 4, pp. 3-25. In Rus.
32. Osipova E.N., Lobova G.A., Isaev V.I., Starostenko V.I. Oil and gas reservoirs of Lower Nurolfka megahollow. Bulletin of the Tomsk Polytechnic University, 2015, vol. 326, no. 1, pp. 14-33. In Rus.
33. Khutorskoy M.D. Vvedenie v geotermiyu [Introduction into geothermic]. Moscow RUDN Publ., 1996. 156 p.
34. Nadai A. Theory of Flow and Fracture of Solids. Vol. 2. New York, McGraw-Hill Book Company, 1963. 705 p.
35. Bogoyavlenskaya V.I., Polyakova I.D. Petroleum potential of the great depths of the South Kara region. Arctic: ecology and economy, 2012, no. 3 (7), pp. 92-103. In Rus.
36. Cserepes L., Yuen D.A., Schroeder B.A. Effect of the mid-mantle viscosity and phase-transition structure on 3D mantle convec-
tion. Physics of the Earth and Planetary Interiors. 2000, vol. 118, no. 1-2, pp. 125-148.
Received: 28 November 2016.
Information about the authors
Valeriy I. Isaev, Dr. Sc., professor, National Research Tomsk Polytechnic University.
Albina A. Iskorkina, postgraduate, National Research Tomsk Polytechnic University.
Vladimir Yu. Kosygin, Dr. Sc., leading researcher, Computer center of the FEB RAS.
Galina A. Lobova, Dr. Sc., associate professor, National Research Tomsk Polytechnic University.
Elizaveta N. Osipova, Cand. Sc., assistant, National Research Tomsk Polytechnic University.
Aleksandr N. Fomin, Dr. Sc., head of the laboratory, Institute of Petroleum Geology and Geophysics named after A.A. Trofimuk SB RAS.
