УДК 553.98:550.8:551.76:552.5(571.62)
ГЕОТЕРМИЯ И ОЦЕНКА НЕФТЕГАЗОВОГО ПОТЕНЦИАЛА БУРЕИНСКОГО БАССЕЙНА
(ДАЛЬНИЙ ВОСТОК РОССИИ)
Прохорова Полина Николаевна1,
prokhorova.polina1988@gmail.com
Развозжаева Елена Петровна2,
rep@itig.as.khb.ru
Исаев Валерий Иванович1,
isaevvi@tpu.ru
1 Национальный исследовательский Томский политехнический университет, Россия, 634050, г. Томск, пр. Ленина, 30.
2 Институт тектоники и геофизики им. Ю.А. Косыгина ДВО РАН, Россия, 680000, г. Хабаровск, ул. Ким Ю Чена, 65.
Актуальность. Буреинский осадочный бассейн - наиболее интересная в нефтегазоносном отношении межгорная впадина слабо-изученной материковой части Дальнего Востока России. Положительная оценка в отношении перспектив обнаружения в отложениях бассейна месторождений нефти и газа обоснована результатами ранее проведённых производственных и научных работ. Цель:дальнейшая аргументация перспектив нефтегазоносности верхнеюрско-нижнемеловых отложений Буреинского бассейна на основе применения отечественного программного обеспечения бассейнового моделирования TeploDialog, обладающего оригинальными особенностями.
Объект: осадочные комплексы средне-верхнеюрский, верхнеюрско-нижнемеловой и меловой Кындалского грабена, наиболее хорошо изученной структуры в составе Буреинского краевого прогиба. Ранее проведенные исследования показали, что в процесс нефтегазообразования вовлечён весь юрско-меловой разрез Буреинского бассейна.
Методы. В комплексе ТeploDialog реализован метод палеотемпературного моделирования, основанный на численном решении уравнения теплопроводности горизонтально-слоистого твердого тела с подвижной верхней границей. В математическую модель непосредственно включены климатический вековой ход температур на земной поверхности, как краевое условие, и пале-отемпературы из определений отражательной способности витринита, как «наблюденные». Метод не требует априорных сведений о природе и величинах глубинного теплового потока, поток определяется решением обратной задачи геотермии в рамках параметрического описания седиментационной истории и истории теплофизических свойств осадочной толщи. Результаты. Установлено, что геотермия грабена, начиная с юрского времени, благоприятствовала процессам генерации жидких углеводородов в породах осадочного чехла, которые могли начаться около 107 млн лет назад, в середине альба. Наиболее интенсивно нефтегенерация происходила в период 107-40 млн лет назад, до середины эоцена. Температурные условия генерации жидких углеводородов сохраняются для талынджанского очага. Наибольшая продолжительность температурных периодов, благоприятных для газообразования, наблюдается в чемчукинском очаге. Условия генерации газа сохраняются для ургальского и чагдомынского очагов. Результаты, полученные с использованием отечественной системы TeploDialog и ранее полученные в системе PetroMod, во многом совпадают. Названные системы характеризуются разными подходами к рассмотрению тектоно-се-диментационных и термических условий генерации УВ. Поэтому результаты настоящих исследований являются дополнительной аргументацией положительных перспектив нефтегазоносности верхнеюрско-нижнемеловых отложений Буреинского бассейна.
Ключевые слова:
Верхнеюрско-нижнемеловые отложения, геотермический режим, очаги генерации углеводородов, Буреинский осадочный бассейн, Дальний Восток России.
Введение
Верхнебуреинская впадина (Буреинский осадочный бассейн) - одна из наиболее интересных в нефтегазоносном отношении межгорных впадин южной части материка в Дальневосточном регионе [1, 2].
Положительная оценка в отношении перспектив обнаружения в отложениях Буреинского бассейна месторождений нефти и газа обоснована результатами ранее проведённых работ по изучению его геологического строения, нефтегазоносности, газогеохимии и подземных вод [3].
В настоящее время при решении задач, связанных с прогнозом нефтегазоносности, широко применяется бассейновое моделирование [4-6]. Этот процесс можно разделить на две относительно независимые части: моделирование истории погру-
жения, тектоники, термического развития (первая часть) и моделирования процессов генерации, миграции, аккумуляции и сохранности углеводородов (УВ) в бассейне (вторая часть). Первая часть в зарубежной литературе определена термином «бассейновое моделирование» («basin modeling»), вторая - «моделирование углеводородных систем» («petroleum system modeling»). Можно придерживаться этой терминологии.
Количество генерированных УВ рассчитывается на основе реконструкции геотемпературного режима нефтематеринских отложений. На региональном этапе геологоразведочных работ технология бассейнового моделирования позволяет выделить очаги генерации УВ и проследить их динамику на протяжении эволюции бассейна. Существую-
DOI 10.18799/24131830/2019/1/51
65
щие на сегодняшний день программно-математические комплексы бассейнового моделирования (ГАЛА, Теш1в и др.) характеризуются разными подходами к рассмотрению тектоно-седимента-ционных и термических условий генерации УВ в бассейне [7, 8].
Первая и единственная оценка реализации генерационного потенциала нефтематеринских отложений Буреинского бассейна в пределах Кын-далского грабена с помощью программного комплекса бассейнового моделирования представлена в работе [9]. Моделирование выполнено с применением программного обеспечения PetroMod Ш компании Шлюмберже.
Цель настоящей работы - дальнейшая аргументация перспектив нефтегазоносности верх-неюрско-нижнемеловых отложений Буреинского бассейна на основе применения отечественного программного обеспечения бассейнового моделирования TeploDialog, обладающего оригинальными превосходными особенностями.
Достижение поставленной цели осуществлялось путем решения следующих задач: 1) актуализация геолого-геофизических, геохимических и нефтегеологических данных в пределах Кындал-ского грабена; 2) одномерное моделирование глубинного теплового потока и термической истории потенциально нефтегазоматеринских свит; 3) выделение по геотемпературному критерию очагов генерации нефти и газа; 4) оценка времени и интенсивности реализации фаз нефте- и газообразования органическим веществом нефтегазомате-ринских отложений.
О методике исследований
Оценка реализации генерационного потенциала нефтегазоматеринских пород в пределах Кын-далского грабена выполнена с помощью комплекса ТeploDialog. В комплексе реализован метод палео-температурного моделирования, основанный на численном решении уравнения теплопроводности горизонтально-слоистого твердого тела с подвижной верхней границей [10]. Метод опирается на решение прямой задачи геотермии, как частный случай найденных решений для набора горизонтальных бесконечных плоскопараллельных слоев, расположенных в нижнем полупространстве [11].
Расчет палеотемператур состоит из двух этапов. На первом, по распределению температур, «наблюденных» в скважине, рассчитывается тепловой поток через поверхность основания осадочного чехла, т. е. решается обратная задача геотермии. На втором этапе, с известным значением теплового потока, решаются прямые задачи геотермии - рассчитываются температуры в заданных точках осадочной толщи (в том числе в материнских отложениях) на заданные моменты геологического времени.
В математическую модель непосредственно включены палеотемпературы из определений
отражательной способности витринита (ОСВ), как «наблюденные». Для перехода от ОСВ к соответствующей геотемпературе используется диаграмма «Линии значений отражательной способности витринита, нанесенные на измененную схему Кон-нона» [12].
Отличия TeploDialog и PetroMod заключаются в разных подходах к определению граничных условий при моделировании тектонической и тепловой истории отложений осадочного бассейна.
Нижним граничным условием является значение теплового потока из основания бассейна осад-конакопления. При моделировании в PetroMod предусматриваются переменные значения плотности теплового потока из основания осадочного чехла, определяемые на основе модели растяжения двухслойной литосферы Маккензи [13]. При пале-отемпературном моделировании в TeploDialog предполагается квазипостоянство значения плотности теплового потока из основания. В TeploDia-log глубинный тепловой поток определяется решением обратной задачи геотермии, в рамках параметрического описания седиментационной истории и истории теплофизических свойств только осадочной толщи, без привлечения сведений о геодинамике ниже основания осадочного разреза.
Верхнее граничное условие - вековой ход поверхностных температур на протяжении моделируемого отрезка истории развития бассейна. В отличие от TeploDialog, в котором используются экспериментальные данные о вековом ходе температур на земной поверхности [14, 15], в PetroMod реализован теоретический расчёт поверхностных па-леотемператур на основе значений палеошироты территории исследований, определяющей интенсивность солярного источника тепла [4].
В методическом подходе, использующем Teplo-Dialog, применяется простой способ (экспересс-расчет) интегрального показателя генерации УВ. Расчет реализуется следующим образом. Решение прямых задач геотермии выполняется на моменты геологического времени, включающие времена начала/завершения формирования каждой свиты, перекрывающей материнскую. Балансовая модель процессов нефтегазообразования [16] позволяет по геотемпературному критерию выполнить выделение очагов интенсивного образования УВ из рассеянного органического вещества (РОВ) материнских отложений. Далее для материнской свиты рассчитывается интегральный показатель плотности генерации углеводородов Я [17]. Расчетное значение Я напрямую зависит от времени нахождения материнской свиты в «верхней зоне газообразования» (ВЗГ), «главной зоне нефтеобразования» (ГЗН), «нижней зоне газообразования» (НЗГ) и от геотемператур этих зон. Применяемый подход позволяет кумулятивно учитывать динамику геотемператур материнских отложений и достаточно просто выполнять пространственно-временную локализацию очагов генерации УВ.
Краткая характеристика территории исследований
В строении осадочного чехла Буреинского бассейна принимают участие четыре комплекса: ни-жне-среднеюрский (дешская ds, синкальтинская sn, эпиканская ep, эльгинская el свиты), средне-верхнеюрский (чаганыйская cg, талынджанская tl свиты), верхеюрско-нижнемеловой (ургальская ur, чагдамынская cg, чемчукинская cm, йорекская ir свиты) и меловой (кындалская kn свита и песчаниковая ps толща).
Кындалская грабен-синклиналь является наиболее хорошо изученной структурой в составе Буреинского краевого прогиба (рис. 1).
Кындалский грабен расположен в центральной части Кындалской грабен-синклинали, имеет север-северовосточную протяженность, около 56 км в длину при максимальной ширине до 17 км. В целом осадочный разрез грабена представлен ритмичным переслаиванием песчаников, алевролитов, аргиллитов, реже конгломератов, с редкими пластами углей и прослоями туфов.
132°30'0"Е
132°45'0"Е
5ri5'0"N-
51OTN-
50°45'0"N-
-51°15'0"N
-5Г0'0"Ы
-50°45'0"N
132С30'0"Е
132М50Е
-1200'
®1А
Рис. 1. Структурная схема Кындалского грабена (по подошве йорекской свиты): 1 - изоглубины подошвы йорекской свиты, м; 2-разломы: а) главные и их буквенное обозначение (К - Кындалский, У - Ургальский, А - Адниканский), б) второстепенные; 3 - глубокие скважины, по которым проводилось моделирование; 4 - территория Кындалской грабен-синклинали; 5 - гидросеть
Fig. 1. Structure map of Kyndal graben (from the bottom of iorekskaya suite): 1 - isodepth of iorekskaya suite bottom, m; 2 - faults: a) main and their letter identification (K - Kindalskiy, U - Urgalskiy, A - Adnikanskiy), b) affluent; 3 - deep wells that were modeled; 4-Kyndal Graben-synclinal area; 5 - hydrographie network
В период с 1965 по 1972 гг. в пределах Кын-далского грабена проведено структурно-картиро-вочное, параметрическое и поисковое бурение шести скважин, в том числе одной параметрической 1ПР на Восточном поднятии. Почти по всему осадочному разрезу были зафиксированы многочисленные проявления УВ, но промышленных скоплений не обнаружено. В 1980-1990 гг. пробурены две глубокие скважины: 1С в центральной части и 1А в его южной части грабена. В скважине 1А получены промышленные притоки газа, что послужило открытию Адниканского газового месторождения.
В 1990-х гг. в разрезе скважины 1А, вскрывшей чагдамынскую, чемчукинскую, йорекскую, кын-далскую свиты и песчаниковую толщу, проведено детальное изучение литологии, фильтрационно-ем-костных свойств пород и геохимии органического вещества. По результатам исследований выявлено, что в процесс нефтегазообразования вовлечён весь юрско-меловой разрез Буреинского осадочного бассейна, при этом исходное органическое вещество имеет преимущественно гумусовый состав [3].
Глубокие скважины 1А, 1С и 1ПР расположены соответственно в южной (Адниканская антиклинальная структура), центральной и северо-восточной (Восточное поднятие) частях грабена (рис. 1).
Далее приводятся краткие сведения о мощности вскрытого разреза, возрастной привязке, наличии размыва, измеренных значениях показателя ОСВ, а также результатах испытаний верхнеюр-ско-меловых пластов. Данные испытаний скважин 1А и 1С изучены и сведены из первичных «дел скважин» (материалы Хабаровского филиала ФБУ «Территориальный фонд геологической информации по Дальневосточному федеральному округу»), данные испытаний скважины 1ПР приняты из монографии [3], значения ОСВ (Щ) приняты из названной монографии.
Скважина 1А. На забое (2998 м) ургальская свита ^-К^, вскрытая на глубине 2600 м. Фиксируется перерыв осадконакопления в раннем мелу, размыв позднемеловых и палеоген-неогеновых отложений. Результаты испытаний кындалской свиты (К1-2): 438-455 м, фонтанный газ, 53 тыс. м3/сут; 993-1010 м, газ, 82 тыс. м3/сут; 1118-1140 м, вода+растворенный газ, 11 м3/сут. Результаты испытаний йорекской свиты (К^: 1202-1215 м, вода+растворенный газ, 7 м3/сут; 1218-1225 м, вода+растворенный газ, 7 м3/сут. Температуры, определенные по ОСВ: 1888 м, 99 °С (£0=0,64); 2196 м, 114 °С (Д°=0,75).
Скважина 1С. На забое (3459 м) талынджан-ская свита (^), вскрытая на глубине 2863 м. Фиксируются перерывы осадконакопления в поздней юре и раннем мелу, размыв позднемеловых и палеоген-неогеновых отложений. Результаты испытаний чемчукинской свиты (К!): 1638-1649 м, слабогазированный фильтрат бурового раствора с пленкой нефти. Результаты испытаний ургаль-ской свиты: 2680-2792 м, слабогазированный
фильтрат бурового раствора. Температуры, определенные по ОСВ: 1211м, 83 °С (Д°=0,52); 1263 м, 83 °С (Д»=0,52); 1355 м, 84 °С (Д£=0,53).
Скважина 1ПР. На забое (3012 м) ургальская свита (^-К^, вскрытая на глубине 2570 м. Фиксируется перерыв осадконакопления в раннем мелу, размыв позднемеловых и палеоген-неогеновых отложений. Имели место повышенные газопоказания при проходке кындальской, чемчукинской и чагдамынской свит. При проходке ургальской свиты в интервале 2720-3002 м наблюдались выбросы нафтено-метановой парафинистой нефти. Результаты испытаний ургальской свиты: 2941-2966 м, газ, 500 м3/сут.
Параметризация и моделирование
Палеореконструкции истории погружения и геотемператур выполнены для трех глубоких скважин: 1А, 1С и 1ПР Кындалского грабена (рис. 1). Пример параметризации тектоно-седиментацион-ной истории и теплофизических свойств осадочного разреза показан в табл. 1.
Термическая история осадочных комплексов восстанавливалась на основе использования фактических данных об ОСВ. Отражательная способность витринита - природного максимального термометра, наиболее предпочтительна для восстановления тепловой истории отложений, так как отражает кумулятивный эффект воздействия температур на РОВ [19, 20].
Скважина 1А пересекает аллохтонный и автохтонный блоки Адниканского надвига на глубине около 1850 м. Предварительный анализ показал, что значения показателя ОСВ для отложений висячего блока являются некондиционными. Для моделирования были отобраны значения показателя ОСВ, измеренные в образцах керна скважины с глубины ниже 1850 м.
Для скважины 1С для моделирования были выбраны значения показателя ОСВ по двум структур-но-картировочным скважинам 14-ск и 16-ск, расположенным в аналогичных структурно-тектонических условиях.
Для скважины 1ПР измеренные значения показателя ОСВ отсутствуют. Построение модели геотемператур выполнялось на основе экстраполяции рассчитанных величин теплового потока, полученных в скважинах 1А и 1С.
Применение вариативности [21] сценариев тек-тоно-седиментационной истории и реконструкций геотермического режима осадочного разреза позволили количественно определить мощность денудации в меловое время (табл. 1).
Контроль результатов моделирования осуществлялся сопоставлением «измеренных» и расчетных геотемператур в скважинах (табл. 2). Сред-неквадратические отклонения фактических значений от расчетных составили ±2-3 °С. Такие «невязки» оптимальны [22], т. к. «наблюденные» геотемпературы (ОСВ) имеют погрешность определения порядка ±2 °С [10].
Таблица 1. Пример параметризации осадочной толщи для моделирования (скважина 1А) Table 1. Example of parameterization of sedimentary strata for modeling (well 1A)
Свита, толща (стратиграфия)* Suite, strata (stratigraphy)* Мощность, м* Thickness, m* Возраст, млн лет назад** Age, million years ago** Время накопления, млн лет Time, million years Плотность, г/см3*** Density, g/cm3*** Теплопроводность, Вт/м-град Thermal conductivity, W/m-deg Температуропроводность, м2/с Heat diffusivity, m2/s Тепловыделение, Вт/м3 Heat production, W/m3
Квартер Quaternary system Q 30 0-1,81 1,81 1,6 0,91 6,5е-007 1,1е-006
P1-N2 -200 1,8-61,7 59,89 1,9 1,17 6,5е -007 1,1е-006
Песчаная Peschanaya p K2-P1 480 61,7-70,6 8,9 1,9 1,17 6,5е -007 1,1е-006
K2 -1100 70,6-98,2 27,6 2,3 1,53 8е-007 1,25е-006
Кындалская Kyndalskaya kn K1-2 1990 98,2-108,4 10,2 2,3 1,53 8е-007 1,25е-006
Иорекская Iorekskaya ir K1 437 108,4-116 7,6 2,5 1,69 8е-007 1,3е-006
Ki - 116-118,2 2,2 - - - -
Чемчукинская Chemchukinskaya cm K1 581 118,2-127,5 9,3 2,5 1,69 8е-007 1,3е-006
Чагдамынская Chagdaminskaya cg K1 382 127,5-136,4 8,9 2,58 1,76 8е-007 1,3е-006
Ургальская Urgalskaya ur J 3-K1 398 136,4-150,8 14,4 2,52 1,71 8е-007 1,3е-006
* - данные литолого-стратиграфическихразбивок взяты из монографии [3]; ** - в качестве основы использована Шкала геологического времени [18]; *** - плотность пород выделенных свит и толщ принята по материалам монографии [3]. Светло-серой заливкой показаны времена накопления потенциально нефтегазоматеринских свит и их параметрическое описание, серой - перерывы осадконако-пления и размывы меловых и палеоген-неогеновых отложений.
* - the data of lithological and stratigraphie differentiation are taken from the monograph [3]; ** - the geological time Scale is used as a basis [18]; *** - the rock density of the selected suites and strata are taken from the monograph [3]. Light-gray fill shows the accumulation times of potentially oil and gas source suites and their parametric description, gray fill shows sedimentation breaks and erosion of Cretaceous and Paleo-gene-Neogene deposits
ных значений плотности теплового потока на скважинах В1, Е1 и Е2 Переяславского грабена Сред-неамурской впадины (48-56 мВт/м2) [23]. Результаты моделирования не противоречат трем значениям теплового потока, измеренным в пределах Буреинского бассейна (40, 56 и 90) мВт/м2 [24], и диапазону значений теплового потока в бассейне Bide-Santang (46-59 мВт/м2) [25].
Результирующие графики истории погружения и геотемператур осадочных комплексов Кын-далского грабена представлены на рис. 2.
Интересно отметить кратковременное, но характерное охлаждение осадочного разреза, начиная с 108 млн лет назад (пиковое «провисание» изотерм, рис. 2), обусловленное как высокой скоростью накопления кындалской свиты, так и ее последующей мощной денудацией.
Рассеянное органическое вещество (РОВ) изучаемых осадочных комплексов преимущественно гумусового типа, поэтому нижняя температурная граница вхождения в верхнюю зону газообразования (ВЗГ) принята 60 °С, в главную зону нефтеоб-разования (ГЗН) - 95 °С, в нижнюю зону газообразования (НЗГ) - 130 С.
В соответствии с температурной градацией зон интенсивной генерации УВ, используя полученную табличную базу термической истории свит и ее графическое представление (рис. 2), проведена оценка степени реализации генерационного потенциала талынджанской, ургальской, чагдамын-ской, чемчукинской и йорекской свит. Выполнена
Таблица 2. Сопоставление измеренных и расчетных геотемператур в скважинах Table 2. Comparison of measured and calculated temperatures in wells
Глубина, м Depth, m Измеренные температуры, °С Measured temperatures, °С Способ измерения Measurement method Расчетные температуры, °С Calculated temperatu-res,°С Разница расчетных и измеренных температур, °С Difference between calculated and measured temperatures, °С
Скважина 1A/Well 1A
1888 99 по ОСВ 102 +3
2196 114 111 -3
Среднеквадратическое отклонение(«невязка»), °С Mean-square deviation («residual error»), °С ±3
Расчетная плотность теплового потока, мВт/м2 Calculated heat flow, mW/m2 49
Скважина 1C/Well 1С
1211 83 по ОСВ 81 -2
1263 83 83 0
1355 84 86 +2
Среднеквадратическое отклонение(«невязка»), °С Mean-square deviation («residual error»), °С ±2
Расчетная плотность теплового потока, мВт/м2 Calculated heat flow, mW/m2 45
Достоверность результатов моделирования, выполненного на скважинах 1А и 1С Кындалского грабена, подтверждается согласованностью полученных расчетных значений плотности теплового потока 45 и 49 мВт/м2 (табл. 2) с порядком расчет-
t, град. С
.,Млн.л?т / ^ ^ £ £ & £ ^ J 1. «пи. лаг/
«з
У' f>
4? / # #
,„30— A 1 J 2-3 tl 2 | ,'-95--.
ш ^
UJ "
ОТ К
О л
UJ А
* °
8 °
ГО =1
ш о
(D
X
?? г
О) -С
а> <D от ^ са о
а>
= S о ■
ы ф О О ОТ "О
= 2 Ш 1£>
Рис. 2. Палеореконструкции тектонической и термической истории в разрезах скважин: А - 1А,Б - 1С, В - 1ПР. 1 - изотермы, "С; 2 - стратиграфическая приуроченность отложений; 3 - изотерма пороговой температуры ГЗН. В верхней части рисунка (А) приведен график палеоклиматического векового хода температур на земной поверхности
Fig 2. Paleoreconstruction of tectonic and thermal history in well sections: A - 1А.Б- 1C,B- 1PR. 1 - isotherms, °C; 2 - stratigraphic confinement of sediments; 3 - isotherm of threshold temperature of main oil window. The upper part of the figure (A) shows a graph of the paleoclimatic secular trend of temperatures on the earth surface
экспресс-оценка плотности генерации УВ - рассчитан интегральный показатель R.
Результаты исследований
Сводные данные экспресс-оценки по интегральному температурно-временному показателю генерации УВ (R) талынджанской, ургальской, чагда-мынской, чемчукинской и йорекской свит приведены в табл. 3.
Анализ истории погружения и геотемператур (рис. 2), интегральных оценок генерации нефти и газа (табл. 3) показывает, что на протяжении истории погружения Буреинского бассейна в пределах Кындалского грабена генерация УВ могла происходить в талынджанской, ургальской, чагдамын-ской, чемчукинской и йорекской свитах. Кындал-ская свита за всю свою историю практически не
входила в ГЗН, лишь эпизодически, незначительно, подошвой достигала верхов ВЗГ.
Около 126,0 млн лет назад породы талын-джанской свиты вошли в ВЗГ. Благоприятные геотемпературные условия для генерации нефти отложениями талынджанской свиты (ГЗН) наступили 107 млн лет назад. Затем, начиная с 98,2 млн лет назад, свита вошла в НЗГ. Результаты расчётов свидетельствуют о том, что с 87,5 млн лет назад по настоящее время в отложениях талынджанской свиты генерируется нефть (ГЗН).
Около 114,0-112,0 млн лет назад породы ур-гальской свиты вошли в ВЗГ. Затем, приблизительно 103,0-101,0 млн лет назад, ургальская свита достигла ГЗН. Наибольшая длительность фазы генерации нефти ургальской свитой наблюдается в районе скважины 1А, где отложения продолжали
Потенциально материнская свита Potentially source suite Экспресс-расчет плотности генерации газа (R), усл. ед Rapid calculation of gas generation density, c.u. Период работы палеоочага генерации газа, млн лет назад; максимальные геотемпературы, "С Working period of gas paleokitchen, million years ago; the maximum geotemperatures, °С Время работы палеоочагов генерации газа, млн лет Working time of gas paleo-kitchen, million years Экспресс-расчет плотности генерации нефти (R), усл. ед Rapid calculation of oil generation density, c.u. Период работы палеоочага генерации нефти, млн лет назад; максимальные геотемпературы, °С Working period of oil paleo-kitchen, million years ago; maximum geotemperatures, °С Время работы палеоочагов генерации нефти, млн лет Working time of oil paleokitchen, million years
Скважина 1A/Well 1A
Ургальская Urgalskaya 17 112,0-103,0; 94 96,0-94,0; 132(ГОГ) 6,0-0,0; 94 17 128 103,0-96,0; 130 94,0-6,0; 120 96
Чагдамынская Chagdaminskaya 47 107,0-102,0; 83 77,5-67,0; 92 35,0-0,0; 94 50 63 102,0-77,5; 121 67,0-35,0; 107 56
Чемчукинская Chemchukinskaya 87 105,0-99,0; 93 87,5-0,0; 93 94 12 99,0-87,5; 107 12
Иорекская Iorekskaya 60 103,0-74,0; 91 71,0-28,0; 78 72 - - -
Скважина tt/Well 1С
Талынджанская Talindzhanskaya 31 126,0-107,0; 87 98,2-87,5; 139(ГОГ) 30 125 107,0-98,2; 130 87,5-0,0; 130 96
Ургальская Urgalskaya 42 114,0-103,0; 92 71,0-69,0; 94 32,0-0,0; 94 45 80 103.0-71.0; 123 69.0-32.0; 110 69
Чагдамынская Chagdaminskaya 72 107,0-101,0; 87 84,5-61,7; 94 48,0-0,0; 94 77 31 101.0-84.5; 109 61.7-48.0; 95 30
Чемчукинская Chemchukinskaya 93 105,0-7,0; 95 98 - - -
Иорекская Iorekskaya 40 102,0-77,5; 86 65,0-36,0; 67 54 - - -
Скважина ШР/Well 1PR
Ургальская Urgalskaya 50 112,0-101,0; 92 76,0-68,0; 94 35,0-0,0; 94 54 62 101.0-76.0; 109 68.0-35.0; 106 58
Чагдамынская Chagdaminskaya 80 107,0-99,0; 89 91,0-62,0; 94 48,0-0,0; 94 85 21 99.0-91.0; 98 62.0-48.0; 96 22
Чемчукинская Chemchukinskaya 88 105,0-3,0; 86 102 - - -
Иорекская Iorekskaya 37 101,0-80,0; 71 67,0-35,0; 69 53 - - -
Таблица 3. Экспресс-оценка генерации углеводородов (температурно-временные интервалы ВЗГ, НЗГ и ГЗН) Table 3. Rapid assessment of hydrocarbon generation (temperature and time intervals of upper gas formation zone, lower gas formation zone and main oil generation zone)
генерировать нефть в период с 103 до 6 млн лет назад. Раньше ургальская свита вышла из ГЗН в районе скважины 1ПР, около 35 млн лет назад. В районе скважины 1С главная фаза нефтегенерации ургальской свитой завершилась 32 млн лет назад. В настоящее время отложения ургальской свиты в пределах территории исследований являются газопроизводящими (ВЗГ).
Около 107,0 млн лет назад чагдамынская свита вошла в ВЗГ. Благоприятные условия для генерации нефти отложениями чагдамынской свиты сформировались приблизительно 102-99 млн лет назад. Затем, около 91,0-77,5 млн лет назад, под влиянием тектонического (общее воздымание территории и последующий размыв) и палеоклимати-ческого факторов, чагдамынская свита вышла из ГЗН. Вторая волна входа в ГЗН, обусловленная нарастанием мощности перекрывающей осадочной толщи, наступила приблизительно 69-62 млн лет назад. Раньше чагдамынская свита вышла из «нефтяного окна» в районе скважин 1ПР и 1С, около 48 млн лет назад. В районе скважины 1А главная фаза нефтегенерации отложениями чагда-мынской свиты завершилась 35 млн лет назад. В настоящее время породы чагдамынской свиты в пределах территории исследований являются газопроизводящими.
Около 105,0 млн лет назад чемчукинская свита вошла в ВЗГ. В период 99,0-87,5 млн лет назад в районе скважины 1А отложения чемчукинской свиты генерировали нефть (ГЗН). Фаза нефтеобра-зования в отложениях чемчукинской свиты в районе скважин 1С и 1ПР места не имела. Приблизительно 105,0-87,5 млн лет назад чемчукинская свита начала генерировать газ. Раньше чемчукин-ская свита вышла из ВЗГ в районе скважины 1С, 7 млн лет назад. В районе скважины 1ПР фаза газообразования завершилась 3 млн лет назад. В настоящее время в районе скважины 1А отложения чемчукинской свиты продолжают генерировать газ.
Около 103,0-101,0 млн лет назад отложения йорекской свиты достигли температуры 60 °С (ВЗГ). Затем, около 80,0-74,0 млн лет назад, йо-рекская свита вышла из ВЗГ. Вторая волна входа в ВЗГ наступила приблизительно 71-65 млн лет назад. Раньше йорекская свита вышла из ВЗГ в районе скважин 1С и 1ПР, около 36-35 млн лет назад. Наибольшая длительность фазы генерации газа йорекской свитой наблюдается в районе скважины 1А, где отложения продолжали генерировать газ в период с 71 до 28 млн лет назад.
Обсуждение результатов
Наибольшая продолжительность фазы нефте-образования (ГЗН) в Кындалском грабене наблюдается в талынджанском очаге генерации нефти -96 млн лет, очаг характеризуется максимальным значением плотности генерации нефти -Е=125 усл. ед (табл. 3). Соизмеримо долго, по геотемпературному критерию, «работает» в пределах
грабена ургальский очаг генерации нефти, продолжительностью от 58 до 96 млн лет, R=62-128 усл. ед. Продолжительность фазы генерации жидких УВ чагдомынским очагом составила от 22 до 56 млн лет, R=21-63 усл. ед. Вклад в генерацию нефти чемчукинским очагом мал, продолжительность генерации жидких УВ - 12 млн лет, R=12 усл. ед. На сегодняшний день геотемпературные условия генерации жидких УВ сохраняются для талынджанского очага.
Наибольшая фаза газообразования (ВЗГ) наблюдается в чемчукинском очаге генерации - от 94 до 102 млн лет, очаг характеризуется максимальным значением плотности генерации газа -R=87-93 усл. ед. Продолжительность периодов газогенерации, выделенных по геотемпературному критерию, чагдомынским, йорекским и ургаль-ским очагами составила 50-85 млн лет ^=47-80), 53-72 млн лет ^=37-60) и 17-54 млн лет ^=17-50), соответственно. Вклад в генерацию газа талынджанским очагом невелик, продолжительность периода газообразования 30 млн лет, R=31 усл. ед.
Вероятность очагов генерации УВ, выделенных по геотемпературному критерию, не противоречит результатам испытаний глубоких скважин на исследуемой территории. Непромышленные притоки нефти получены при испытании ургальской свиты, плёнки нефти - при испытании чемчукинской свиты.
Повышенные газопоказания зафиксированы при испытании чемчукинской, чагдомынской и ур-гальской свит. По данным бурения скважины 1А получен промышленный приток газа из кындал-ской свиты. Вполне вероятно, что углеводороды, генерируемые нижнемеловыми комплексами, могли мигрировать и аккумулировать в тектонически экранированных ловушках, образовавшихся во фронте Адниканского надвига (рис. 1).
Сопоставительным анализом результатов настоящей работы, выполненной с использованием российского комплекса палеотемпературного моделирования TeploDialog, с предшествующей оценкой территории исследований, выполненной ранее с применением импортной системы Ре^о-Mod, установлено, что термический режим прогнозируемых очагов генерации УВ (продолжительность периодов нефте- и газообразования, максимальные температуры периодов интенсивной генерации нефти и газа) отличается несущественно. Для примера сопоставления ниже представлена оценка термического режима очагов генерации УВ в районе скважины 1А, полученная ранее с использованием системы Petromod [9].
Для отложений ургальской свиты: период ВЗГ - 128-110 млн лет назад; периоды ГЗН -110-100 млн лет назад, с 91 млн лет назад по настоящее время; период НЗГ - 100-91 млн лет назад.
Для отложений чагдомынской свиты: периоды ВЗГ - 119-106 млн лет назад, 75-67 млн лет на-
зад, с 35 млн лет назад по настоящее время; периоды ГЗН - 106-75 млн лет назад, 67-35 млн лет назад.
Для отложений чемчукинской свиты: периоды ВЗГ - 110-100 млн лет назад; с 85 млн лет назад по настоящее время; период ГЗН - 100-85 млн лет назад.
Для отложений йорекской свиты: период ВЗГ -105-20 млн лет назад. Фаза нефтеобразования места не имела.
Если провести сопоставление результатов Те-ploDialog и PetroMod по всем трем скважинам (1А, 1С и 1ПР), то перечни выделенных по геотемпературному критерию очагов генерации нефти (та-лынджанский, ургальский, чагдомынский) и очагов газогенерации (ургальский, чагдомынский, чемчукинский, йорекский) полностью совпадают.
Заметное отличие результатов моделирования TeploDialog и PetroMod состоит в том, что в системе PetroMod получены периоды ВЗГи НЗГ более широкие. По-видимому, это обусловлено разными подходами к определению граничных условий при моделировании тектонической и тепловой истории отложений осадочного бассейна, о чем говорилось в разделах введения и методики исследований. Вместе с тем, можно отметить, что рассчитанные в Te-ploDialog значения плотности глубинного теплового потока (45-49 мВт/м2)согласуются со значениями современной плотности теплового потока, полученными с помощью PetroMod (48-49 мВт/м2).
СПИСОК ЛИТЕРАТУРЫ
1. Маргулис Л.С. Нефтегазоперспективные толщи Дальнего Востока России // Нефтегазовая геология. Теория и практика. -2012. - Т. 7. - № 4. URL: http://www.ngtp.ru/rub/2/65_ 2012.pdf (дата обращения 01.08.2018)
2. Razvozjaeva 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.
3. Осадочные бассейны Востока России. Т. 4. Буреинский осадочный бассейн: геолого-геофизическая характеристика, геодинамика, топливно-энергетические ресурсы / гл. ред. А.И. Хан-чук, отв. ред. Г.Л. Кириллова. - Владивосток: Дальнаука, 2012. - 360 с.
4. Hantschel T., Kauerauf A.I. Fundamentals of basin and petroleum systems modeling. - Heidelberg: Springer, 2009. - 476 p.
5. Galushkin Yu. Non-standard Problems in Basin Modelling. -Switzerland: Springer, 2016. - 274 p.
6. Reconstruction of the Cenozoic History of Hydrocarbon Fluids from Rifting Stage to Passive Continental Margin Stage in the Huizhou Sag, the Pearl River Mouth Basin / Yajun Li, Shu Jiang, Zhenglong Jiang, Hao Liu and Bingxi Li // Geofluids. - 2017. -UNSP 4358985. - 32 p.
7. Galushkin Yu.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.
8. 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,
Выводы
1. Геотермия Кындалского грабена, начиная с юрского времени, благоприятствовала процессам генерации жидких углеводородов в породах осадочного чехла, которые могли начаться около 107 млн лет назад, в середине альба. Наиболее интенсивно нефтегенерация происходила в период 107-40 млн лет назад, до середины эоцена. На сегодняшний день температурные условия генерации жидких углеводородов сохраняются для талынджанского очага.
2. Наибольшая продолжительность температурных периодов, благоприятных для газообразования, наблюдается в чемчукинском очаге. В настоящее время условия генерации газа сохраняются для ургальского и чагдомынского очагов.
3. Результаты выделения по геотемпературному критерию очагов генерации нефти и очагов газогенерации, полученные с использованием систем TeploDialog и PetroMod, во многом и в основном совпадают.
4. Названные системы характеризуются разными подходами к рассмотрению тектоно-седимента-ционных и термических условий генерации УВ. Поэтому результаты настоящих исследований являются дополнительной аргументацией положительных перспектив нефтегазоносности верхнеюрско-нижнемеловых отложений Буре-инского бассейна.
P.I. Safronov, S.A. Gus'kov, 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.
9. Razvozzhaeva E.P., Prokhorova P.N., Lapkovskii V.V. Numerical Modeling of the Tectonic and Thermal History of the Kyndal Graben of the Bureya Basin (Far East of Russia) // Russian Journal of Pacific Geology. - 2017. - V. 11. - № 3. - P. 205-222.
10. Mesozoic-Cenozoic Climate and Neotectonic Events as Factors in Reconstructing the Thermal History of the Source-Rock Bazhe-nov Formation, Arctic Region, West Siberia, by the Example of the Yamal Peninsula / V.I. Isaev, A.A. Iskorkina, G.A. Lobova, V.I. Starostenko, S.A. Tikhotskii, A.N. Fomin // Izvestiya, Physics of the Solid Earth. - 2018. - V. 54. - № 2. - P. 310-329.
11. Generalization of the Rayleigh-Tikhonov stationary geothermal problem for a horizontal layer / V.I. Starostenko, R.I. Kutas, V.N. Shuman, O.V. Legostaeva // Izvestiya, Physics of the Solid Earth. - 2006. - V. 42. - № 12. - P. 1044-1050.
12. Хант Дж. Геохимия и геология нефти и газа. - М.: Мир, 1982. - 704 с.
13. Royden L., Keen C.E. Rifting process and thermal evolution of the continental margin of eastern Canada determined from subsidence curves // Earth and Planet. Sci. Letters. - 1980. - V. 51. -P. 342-361.
14. Геотермия арктических морей / М.Д. Хуторской, В.Р. Ахмед-зянов, А.В. Ермаков, Ю.Г. Леонов, Л.В. Подгорных, Б.Г. Поляк, Е.А. Сухих, Л.А. Цыбуля. - М.: ГЕОС, 2013. - 232 с.
15. Искоркина А.А. Палеоклиматические факторы реконструкции термической истории нефтематеринской баженовской свиты арктического региона Западной Сибири // Известия Томского политехнического университета. Инжиниринг георесурсов. - 2016. - Т. 327. - № 8. - C. 59-73.
16. Модель катагенеза органического вещества (на примере баже-новской свиты) / Л.М. Бурштейн, Л.В. Жидкова, А.Э. Конто-рович, В.Н. Меленевский // Геология и геофизика. - 1997. -Т. 38. - №6. - С. 1070-1078.
17. Лобова Г.А., Попов С.А., Фомин А.Н. Локализация прогнозных ресурсов нефти юрско-меловых НГК Усть-Тымской мегав-падины // Нефтяное хозяйство. - 2013. - № 2. - С. 36-40.
18. Шкала геологического времени / У.Б. Харленд, А.В. Кокс, П.Г. Ллевеллин, К.А.Г. Пиктон, А.Г. Смит, Р. Уолтерс. - М.: Мир, 1985. - 140 с.
19. Yalcin M.N., Littke R., Sachsenhofer R.F. Thermal history of sedimentary basins. - Berlin: Springer Verlag, 1997. - 167 p.
20. Фомин А.Н. Катагенез органического вещества и нефтегазо-носность мезозойских и палеозойских отложений Западно-Сибирского мегабассейна. - Новосибирск: ИНГГ СО РАН, 2011. -331 с.
21. Тектоно-седиментационная интерпретация данных геотермии при выявлении и оценке позднеэоценовой эрозии на арктических месторождениях углеводородов (п-ва Ямал) / В.И. Исаев, В.И. Старостенко, Г.А. Лобова, А.Н. Фомин, А.К. Исагалиева // Известия Томского политехнического университета. Инжиниринг георесурсов. - 2017. - Т. 328. - № 7. - С. 19-31.
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. 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.
24. Thermal structure of lithosphere in Central Asian and Pacific belts and their adjacent cratons, from data of geoscience transects / P.Yu. Gornov, M.V. Goroshko, Yu.F. Malyshev, V.Ya. Podgornyi // Russian Geology and Geophysics. - 2009. -V. 50. - № 5. - P. 485-499.
25. Guo C., Qin Y., Lu L. Terrestrial heat flow and geothermal field characteristics in the Bide-Santang basin, western Guizhou, South China // Energy Exploration & Exploitation. - 2018. -V. 36. - №5. - P. 1114-1135.
Поступила 31.10.2018 г.
Информация об авторах
Прохорова П.Н., аспирант отделения геологии Инженерной школы природных ресурсов Национального исследовательского Томского политехнического университета.
Развозжаева Е.П., кандидат геолого-минералогических наук, заведующая лабораторией тектоники осадочных бассейнов Института тектоники и геофизики им. Ю.А. Косыгина ДВО РАН.
Исаев В.И., доктор геолого-минералогических наук, профессор отделения геологии Инженерной школы природных ресурсов Национального исследовательского Томского политехнического университета.
UDC 553.98:550.8:551.76:552.5(571.62)
GEOTHERMY AND ESTIMATION OF HYDROCARBON POTENTIAL OF THE BUREYA BASIN (RUSSIAN FAR EAST)
Polina N. Prokhorova1,
prokhorova.polina1988@gmail.com
Elena P. Razvozzhaeva2,
rep@itig.as.khb.ru
Valeriy I. Isaev1,
isaevvi@tpu.ru
1 National Research Tomsk Polytechnic University, 30, Lenin Avenue, Tomsk, 634050, Russia.
2 Institute of Tectonics and Geophysics named after Yu.A. Kosygin, 65, Kim Y. Chen street, Khabarovsk, 680000, Russia.
The relevance. The Bureya sedimentary basin is the most interesting in respect of oil and gas intermountain basin of frontier mainland of Russian Far East. An upbeat assessment of prospects of oil and gas field discovery in the basin is proved by the results of previous prospecting work and scientific research.
The aim of this paper is further argumentation of hydrocarbon prospects of upper Jurassic-lower Cretaceous deposits in the Bureya basin based on the use of domestic basin modeling software TeploDialog, having original features.
Study object is middle-upper Jurassic, upper Jurassic-lower Cretaceous and Cretaceous sedimentary sequences of the Kyndal graben, the most studied structure in the Bureya fore deep. The previous investigations have revealed that the whole Jurassic-Cretaceous section of the Bureya basin is involved in oil-and-gas formation.
Methods. In the software TeploDialog the paleotemperature modeling method is implemented based on the numerical solution of heat-transfer equation of a horizontally layered solid with a movable upper boundary. Mathematical model includes climatological secular trend of temperatures on the earth surface (boundary condition) andpaleotemperatures from vitrinite reflectance evaluation (observational). The method does not require a priori information about the nature and magnitude of the deep heat flow, the heat flow is determined by the solution of the geothermal inversion within the framework of parametric description of the sedimentation history and the history of the thermophysical properties of sedimentary strata.
Results. The study found that geothermy of graben since Jura permitted formation of liquid hydrocarbons in the sedimentary cover rocks, which could begin about 107 million years ago, in the middle of the Alb. The most intensive oil generation took place in the period of 107-40 million years ago, until the middle of the Eocene. Temperature conditions for the generation of liquid hydrocarbons are still the same for talyndzhansky kitchen. Maximal duration of temperature periods, favourable for gas formation is observed in chemcukin-sky kitchen. The conditions for gas generation are still the same for urgalsky and chegdomynsky kitchens. The results obtained using the domestic software TeploDialog and previously obtained in software PetroMod, strongly aligned. These computer programs are characterized by different approaches to consideration of tectonic-sedimentation and thermal conditions of hydrocarbon generation. Therefore, the results of this study are an additional argument for the positive prospects of oil and gas potential of the upper Jurassic-lower Cretaceous deposits of the Bureya basin.
Key words:
Upper Jurassic-lower Cretaceous deposits, geothermal mode, hydrocarbon kitchens, Bureya sedimentary basin, Russian Far East.
REFERENCES
1. Margulis L.S. Oil and gas prospective strata of Russian far East. Oil and gas geology. Theory and practice, 2012, vol. 7, no. 4. InRus. Available at: http://www.ngtp.ru/rub/2/65_2012.pdf (accessed 1 August 2018).
2. Razvozjaeva 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.
3. Osadochnye basseyny Vostoka Rossii. Vol. 4. Bureinskiy osadochny basseyn: geologo-geofizicheskaya kharakteristika, geodinamika, toplivno-energeticheskie resursy [Series «Sedimentary basins of East Russia». Vol. 4. The Bureya sedimentary basin: geological-geophysical characteristics, geodinamics, and fuel and energy resources]. Editor-in-chef academician A.I. Khanchuk, ed. by G.L. Kirillova. Vladivostok, Dalnauka Publ., 2012. 360 p.
4. Hantschel T., Kauerauf A.I. Fundamentals of basin and petroleum systems modeling. Heidelberg, Springer, 2009. 476 p.
5. Galushkin Yu. Non-standard Problems in Basin Modelling. Switzerland, Springer, 2016. 274 p.
6. Yajun Li, Shu Jiang, Zhenglong Jiang, Hao Liu and Bingxi Li. Reconstruction of the Cenozoic History of Hydrocarbon Fluids from Rifting Stage to Passive Continental Margin Stage in the Huizhou Sag, the Pearl River Mouth Basin. Geofluids, 2017. UNSP 4358985. 32 p.
7. Galushkin Yu.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.
8. Kontorovich A.E., Burshtein L.M., Malyshev N.A, Safronov P.I., Gus'kov S.A., Ershov S.V., Kazanenkov V.A., Kim N.S., Kontorovich V.A., Kostyreva E.A., Melenevsky V.N., Livshits V.R., Polyakov A.A., Skvortsov M.B. Historical-geological 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.
9. Razvozzhaeva E.P., Prokhorova P.N., Lapkovskii V.V. Numerical Modeling of the Tectonic and Thermal History of the Kyndal
Graben of the Bureya Basin (Far East of Russia). Russian Journal of Pacific Geology, 2017, vol. 11, no. 3, pp. 205-222.
10. Isaev V.I., Iskorkina A.A, Lobova G.A., Starostenko V.I., Tikhot-skii S.A., Fomin A.N. Mesozoic-Cenozoic Climate and Neotecto-nic Events as Factors in Reconstructing the Thermal History of the Source-Rock Bazhenov Formation, Arctic Region, West Siberia, by the Example of the Yamal Peninsula. Izvestiya, Physics of the Solid Earth, 2018, vol. 54, no. 2, pp. 310-329.
11. Starostenko V.I., Kutas R.I., Shuman V.N., Legostaeva O.V. Generalization of the Rayleigh-Tikhonov stationary geothermal problem for a horizontal layer. Izvestiya, Physics of the Solid Earth, 2006, vol. 42, no. 12, pp. 1044-1050.
12. Khant Dzh. Geokhimiya i geologiya nefti i gaza [Geochemistry and geology of oil and gas]. Moscow, Mir Publ., 1982. 704 p.
13. Royden L., Keen C.E. Rifting process and thermal evolution of the continental margin of eastern Canada determined from subsidence curves. Earth and Planet. Sci. Letters, 1980, vol. 51, pp. 342-361.
14. Khutorskoy M.D., Akhmedzyanov V.R., Ermakov A.V., Leo-nov Yu.G., Podgornykh L.V., Polyak B.G., Sukhikh E.A., Tsybu-lya L.A. Geotermiya arkticheskikh morey [Geochemistry of the Arctic seas]. Moscow, GEOS Publ., 2013. 232 p.
15. Iskorkina A.A. Paleoclimate factors of reconstructing thermal history of the petromaternal bazhenov suite of the Arctic region in Western Siberia. Bulletin of the Tomsk Polytechnic University. Geo Аssets Engineering, 2016, vol. 327, no. 8, pp. 59-73. In Rus.
16. Burshteyn L.M., Zhidkova L.V., Kontorovich A.E., Melenev-skiy V.N. Model of organic matter catagenesis (on the example of Bazhenov formation). Geology and Geophysics, 1997, vol. 38, no. 6, pp. 1070-1078. In Rus.
17. Lobova G.A., Popov S.A., Fomin A.N. Localization of predicted oil resources in the Jurassic and Cretaceous hydrocarbon play of Ust'-Tym megadepression. Oil Industry, 2013, no. 2, pp. 36-40. In Rus.
18. Kharlend U.B., Koks A.V., Llevellin P.G., Pikton K.A.G., Smit A.G., Uolters R. Shkala geologicheskogo vremeni [Geologic time scale]. Moscow, Mir Publ., 1985. 140 p.
19. Yalcin M.N., Littke R., Sachsenhofer R.F. Thermal history of sedimentary basins. Berlin, Springer Verlag, 1997. 167 p.
20. Fomin A.N. Katagenez organicheskogo veshchestva i neftegazo-nosnost mezozoyskikh i paleozoyskikh otlozheniy Zapadno-Sibir-skogo megabasseyna [Organic matter catagenesis and petroleum potential of Mesozoic and Paleozoic deposits of the West Siberian megabasin]. Novosibirsk, IPGG SB RAS Publ., 2011. 331 p.
21. Isaev V.I., Starostenko V.I., Lobova G.A., Fomin A.N., Isagalie-va A.K. Tectonic-sedimentation interpretation of the geother-mics data when identifying and assessing the late Eocene erosion on the Arctic hydrocarbon fields (Yamal peninsula). Bulletin of the Tomsk Polytechnic University. Geo Аssets Engineering, 2017, vol. 328, no. 7, pp. 19-31. In Rus.
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. Gulenok R.Yu., Isaev V.I., Kosygin V.Yu., Lobova G.A., Starostenko 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.
24. Gornov P.Yu., Goroshko M.V., Malyshev Yu.F., Podgornyi V.Ya. Thermal structure of lithosphere in Central Asian and Pacific belts and their adjacent cratons, from data of geoscience transects. Russian Geology and Geophysics, 2009, vol. 50, no. 5, pp. 485-499.
25. Guo C., Qin Y., Lu L. Terrestrial heat flow and geothermal field characteristics in the Bide-Santang basin, western Guizhou, South China. Energy Exploration & Exploitation, 2018, vol. 36, no. 5, pp. 1114-1135.
Received: 31 October 2018.
Information about the authors
Polina N. Prokhorova, postgraduate, National Research Tomsk Polytechnic University.
Elena P. Razvozzhaeva, Cand. Sc., head of the Laboratory of Sedimentary Basin Tectonics, Institute of Tectonics and Geophysics named after Yu.A. Kosygin.
Valeriy I. Isaev, Dr. Sc., professor, National Research Tomsk Polytechnic University.