Для выполнения всех видов геологического моделирования необходимым условием является умение фиксировать результаты стратиграфических построений. Традиционно сохранялась такая информация в каталогах (таблицах) положения границ в разрезах скважин или обнажений [1]. Каталоги сопровождались пояснительными записками и схемами. Поиск, идентификация, согласование информации в каталогах и способ ее использования в геологических построениях были заботой геолога-интерпретатора.
При переходе к массовой обработке геолого-геофизических данных на ЭВМ, особенно с появлением первых СУБД (систем управления базами данных) появилась возможность иметь автоматизированные каталоги. А это потребовало формализации необходимого минимума базовых стратиграфических понятий (разработки модели данных), реализации специализированного программного обеспечения для заполнения и сопровождения каталога, непосредственного заполнения баз данных каталога и, наконец, подготовки технологической среды для его использования.
Принципы построения и модель данных автоматизированного каталога
База литолого-стратиграфических данных является частью интегрированной базы данных Научно-аналитического центра рационального недропользования (ИБД ЦРН). ИБД ЦРН выполнена на основе СУБД «Oracle-8» и эксплуатируется на протяжении семи лет. Одним из центральных объектов ИБД ЦРН является фонд поисково-разведочных скважин. Для каждой из скважин в ИБД хранятся её паспортные данные (наименование, координаты и т.п.). Всего ИБД содержит информацию по 14 тысячам поисково-разведочных скважин, пробуренных на территории Тюменской области. Для внутренней идентификации в ИБД каждой из скважин присваивается внутрисистемный уникальный идентификатор (номер UN). Вся проблемная информация по стволу скважин идентифицируется со скважиной по её уникальному номеру (В.М. Яковлев, 1999).
На рис. 1 представлена упрощенная структура стратиграфического раздела БД.
Рис. 1. Упрощенная структура стратиграфического раздела БД.
Геологическое представление разрезов, вскрытых в стволах скважин, в ИБД хранится в виде литологических колонок (блоки «Литология» и «Породы») и каталога литолого-стратиграфических разбивок (блоки «Разбивки», «Районы», «Номенклатура», «Облик»).
Литологические колонки скважин кодируются при помощи словаря литологических кодов, содержащего наименования двадцати типов пород, наиболее часто встречающихся в описаниях продуктивных отложений ХМАО (табл. 1). При необходимости таблица кодов расширяется.
Таблица 1.
Для каждого проинтерпретированного пропластка в таблице «Литология» заводится отдельная запись. Запись содержит ссылку на скважину (ее UN), глубину залегания кровли пропластка и его толщину (в метрах) и двух-трехбуквенный код приписанной разновидности пород. Кроме того, для слоя (пропластка) в ИБД могут храниться его основные петрофизические параметры (альфа ПС, сопротивление пород, пористость, проницаемость, продуктивность, характер насыщения и т.п.). Результаты послойной интерпретации могут объединяться в логические группы – варианты интерпретации, выполненные в рамках отдельных проектов.
При разработке структуры каталога мы пошли по пути естественного разделения общестратиграфических знаний (семантической модели предметной области) и сведений (модели седиментационного бассейна) от «фактуры» – конкретных глубин залегания границ [2]. Представленное в базе данных формализованное описание стратиграфических таксонов основывается на текущих версиях «Стратиграфического кодекса» [3] и «Региональных стратиграфических схем…» [4].
В настоящее время основой всех стратиграфических построений является Международная стратиграфическая шкала (МСШ). Внутри провинции Стратиграфический кодекс предписывает считать свиты и подсвиты основными таксономическими единицами местных стратиграфических шкал. И, наконец, имеются литолого-стратиграфические подразделения: маркирующие горизонты, пачки, пласты, пропластки и слои, используемые в качестве вспомогательных по отношению к стратонам местных стратиграфических шкал. В практике геологоразведочных работ на нефть и газ основным объектом анализа являются проницаемые пласты или их группы, именуемые резервуарами и комплексами. Согласование всех этих шкал мы производим через шкалу абсолютного возраста Г.П. Мясниковой и В.И. Шпильмана [5] с использованием элементов анализа седиментологической цикличности отложений [6].
Наиболее простую иерархическую структуру имеет Международная стратиграфическая шкала. В ряду: группа, система, отдел, ярус, подъярус – возрастной диапазон таксона верхнего уровня полностью слагается из суммы возрастных диапазонов таксонов старшего из подчиненных уровней.
Так как в «Региональных стратиграфических схемах…» местные таксоны датируются через возраст таксонов МСШ, то все местные шкалы региона, основные и вспомогательные, являются частными случаями детализации и уточнения МСШ.
В связи с тем, что границы местных шкал обычно приурочиваются к литологическим маркерам, необходимо оговорить свое понимание возрастного скольжения литологических границ.
При описании циклостратиграфической модели осадочного бассейна обычно указывается набор стратиграфических уровней (возрастов), ограничивающих во времени интервалы однонаправленного изменения условий осадконакопления [6]. Далее, согласно принципу Головкинского-Вальтера, предполагается, что внутри каждого из таких интервалов в результате осадконакопления образуется ансамбль взаимосогласованных литологических границ, названия которых меняются в зависимости от авторских предпочтений: сиквенс, циклит и т.п. – суть остается.
Для Западно-Сибирского мезозойского осадочного бассейна существуют несколько циклостратиграфических схем [6,7 и др.]. При составлении циклической шкалы мы, учитывая прикладной характер своей работы, ориентировались на существующую прак тику выделения нефтегазоносных комплексов и резервуаров [8].
Для нефтегазоносных отложений была выработана шкала требуемой детальности. В качестве примера опишем процедуру выделения временных интервалов субрегиональных клиноциклитов валанжина. Считается, что граница валанжина и готерива проходит внутри сармановской пачки, залегающей над пластом БС8 Сургутского района. Граница же берриаса и валанжина фаунистически не установлена и условно принята нами проходящей внутри пласта БВ14. В возрастном интервале валанжина таким образом умещаются девять циклитов – нефтегазоносных резервуаров от БС8 до БВ13 и частично резервуар БВ14. То есть, нам необходимо разделить временной интервал валанжина 139.5–130=9.5 млн. лет на десять частей. Соответственно, мы принимает временной шаг циклитов валанжинских резервуаров равным 0.9 млн. лет. Аналогично мы поступили и с другими нефтеносными горизонтами, но с учетом общей динамики изменения скоростей осадконакопления в интересующих нас районах. Так если для неокома мы приняли временной шаг циклитов, равным 0.9 млн. лет, то в нижней и средней юре аналогичные циклы (с учетом границ таксонов МСШ) были приняты равными 2 миллионам лет.
Мы отчетливо понимаем, что построение подобных шкал – большая и самостоятельная научная задача. Несколько легковесная самоуверенность нашей реализации объясняется критерием практической потребности. Для наших целей кодирования, упорядочения и синхронизации таксонов местных шкал точность определения возрастов не имеет принципиального значения – главное, чтобы они были равными у синхронных объектов и возрастали при переходе от молодых границ к более древним. Единственная цель, ради которой мы постарались достичь разумной аккуратности и точности, это возможность построения для отдельных скважин или участков территории графиков «глубина залегания – абсолютный возраст пород» [9], используемых в методике оценки потенциальных ресурсов В.И. Шпильмана [10]. На рис. 2 приведен такой график для Приразломной скв. 184-Р. Редкие квадраты (1) – отражают положение границ ярусов МСШ, густо посаженные ромбы (0) – уточнение графика за счет границ местных таксонов.
Для субрегиональных границ кровли проницаемых песчано-алевритовых отложений резервуаров (кровли циклитов) интервал их возрастного скольжения ограничивался «сверху», принимался равным точ, но возрасту одной из границ циклостратиграфической шкалы региона. При описании конкретной местной (районной) шкалы ее циклиты описывались ограниченными сверху и снизу граничными возрастами обобщенной шкалы. Таким образом, обеспечивалась синхронизация по абсолютному возрасту местных шкал через общее описание всех региональных и субрегиональных циклических объектов.
Рис. 2. График «глубина — абсолютный возраст» для Приразломной скважины 184-Р
Согласные границы свит и подсвит, там где это было возможно, привязывались к границам циклических объектов: кровлям или подошвам пластов. Там, где циклические объекты не выделялись, мы условно приняли границы свит синхронными, приписав им возрасты соответствующих таксонов МСШ прямо, в случае совпадения их возрастов по «Региональным стратиграфическим схемам…», либо с минимальным пересчетом и корректировкой в случае их несовпадения.
При наличии практической потребности каждая из местных шкал была детализирована до уровня возрастных интервалов, вмещающих основные продуктивные пласты и их пропластки. Критерий дробления возрастных интервалов – равномерность, удобство и возможность локального (по разрезу) уточнения, углубления и расширения местных шкал. Требование синхронизации местных шкал не только по региональным и субрегиональным циклическим объектам, но и по пластам, и пропласткам нами не выдвигалось. Совпадение возрастов пластов и пропластков возникало из стремления к минимизации общего числа ключевых возрастов, используемых при описании литостратиграфических объектов.
Описание возрастных соотношений и соподчиненности подразделений МСШ и всех литолого-стратиграфических объектов местных шкал от свит до пластов и пропластков содержится в таблице номенклатуры.
Для каждого из 46 нефтегазоносных районов Западно-Сибирской плиты описано множество границ местных таксонов, всего более трех тысяч записей. В таблице приводится пример описания границ, а именно — две записи номенклатуры Вартовского района. Две горизонтальные подстроки содержат имена таксонов одного ранга, соприкасающихся на данной границе.
Единичная запись таблицы «Номенклатура» содержит код района, условные абсолютные возрасты по уточненной и детализированной шкале Г.П. Мясниковой и В.И. Шпильмана [5], укороченное имя границы (NK) и имена свит, подсвит, пачек, пластов или пропластков, кровлей (нижняя подстрока) или подошвой (верхняя подстрока) которых является эта граница. Так, граница алымской и ванденской свит проходит по подошве пласта АВ13 (кровле покрышки пласта АВ2) и имеет возраст 118 миллионов лет.
Таблица 2.
Укороченное имя границы (NK) состоит из восьми первых символов имени старшего таксона, ограничиваемого этой границей сверху. Так, укороченное имя «ВАНДЕНСК» отражает тот факт, что возраст 118 млн. лет в Вартовском районе является кровлей ванденской сви ты, а следующая граница — всего лишь кровлей пласта АВ2.
Всего в базе данных описано более трех тысяч границ литолого-стратиграфических объектов осадочного чехла Западно-Сибирского осадочного бассейна на территории Тюменской области.
Для основных таксонов местных шкал (свит и подсвит) в таблице «Облик» может быть создано их стандартное описание по схеме, принятой в геологических словарях, включая автора, первоисточник, год введения, облик пород, возраст, толщины, характер границ, топономический смысл названия и т.п. Эти описания могут использоваться в справочном режиме или выдаваться в форме стандартных отчетов. Ниже приводится пример статьи такого отчета с кратким описанием новомихайловской свиты.
Таблица 3.
Кроме возможности использования в справочном режиме и в отчетах, эти хранимые в БД сведения программно помещаются в описания соответствующих таксонов при генерации сводных и гипостратотипических разрезов. Реально в ИБД хранится только шестьсот таких описаний.
Каталог разбивок построен в базе данных как таблица глубин залегания кровли стратиграфических под разделений. Каждая запись каталога содержит уникальный номер скважины, глубину залегания границы, код ее типа, код авторского подкаталога и укороченное имя таксона (NK), к которому относится эта граница. Список используемых типов границ приводится в табл. 2.
Таблица 2.
Технически собственно каталог разбивок построен из двух таблиц. В одной их них (RAZA) хранятся рабочие варианты авторских каталогов, в другой (RAZO) – одновариантный сводный каталог. Сводный каталог охватывает в идеале полный фонд скважин. Он составлен в результате сведения и переинтерпретации всех доступных документов – каталогов ручного обращения. (Каталог содержит разбивки из 56 первоисточников с 1964 по 2000 год, всего более 300 тысяч записей. Анализ результатов заполнения был дан в [11].) Сводный каталог автоматически пополняется и корректируется при работе кураторов со своими авторскими каталогами пользователей с особыми полномочиями, курирующих определенную территорию.
Авторский каталог может содержать самостоятельный вариант корреляции разрезов скважин, либо локальное (по разрезу или площади) дополнение к сводному варианту каталога. В последнем случае, при обращении к авторскому каталогу авторские разбивки могут дополнять или замещать собой соответствующие разбивки сводного каталога (см.табл. 3).
Таблица 3.
Пользователь (оператор или программа) видит сводный ката лог, на который «сверху» наложены разбивки авторского каталога. Проверка согласования авторского каталога со сводным каталогом при изменении авторского или сводного каталогов осуществляется программными средствами синхронно с внесением изменений.
Краткое описание реализованных информационных технологий
Администрирование и работа со стратиграфическим разделом ИБД ведется в основном в диалоговом режиме. Программное обеспечение написано в среде «Oracle-Forms» и гарантирует сохранение логической целостности базы. Массовый ввод информации может осуществляться и в пакетном режиме, также при жестком контроле вводимых данных.
Существуют программы для разнообразной выдачи стандартных форм каталогов и результатов интерпретации ГИС в виде файлов и структурированных документов ручного обращения. Пример выдачи каталога разбивок скважин по лицензионному участку приводится выше.
В заголовке таблицы приводится номер скважины по паспорту и ее уникальный внутрисистемный номер (UN). Глубины залегания границ извлечены из авторского и сводного каталогов. Границы в каталоге упорядочены по условному абсолютному возрасту, полные имена границ, извлеченных из каталогов, и их ранги (отражены отступом) расшифрованы в таблице «Номенклатура» по коротким именам границ, присутствующих в записях каталога.
Всего за время эксплуатации ИБД конечными пользователями было сгенерировано более тысячи каталогов ручного обращения, различной структуры и назначения.
Примером использования технологий подготовки документов ручного обращения является опубликование каталога посвитных разбивок [11]. Каталог содержит литолого-стратиграфические разбивки по восьми тысячам поисково-разведочных скважин, пробуренных на территории Ханты-Мансийского автономного округа Тюменской области. Представленные разбивки проинтерпретированы в рамках действующей (1991 года) региональной стратиграфической схемы.
В ГУП НАЦ РН ХМАО выполняется весь комплекс региональных и локальных геологических построений на основании информации поисково-разведочных скважин. В большинстве из них участвуют и стра тиграфические данные.
Одной из наиболее массовых задач является картирование поверхностей и свойств геологических тел. Исходными данными процедур картирования являются свойства картируемого объекта, определенные в некоторой точке пространства с прямоугольными декартовыми координатами. Обычно пространство двумерное (плоскость), а свойство замеряется по арифметической шкале (числом).
При создании технологии подготовки данных для картирования мы исходили из того, что, во-первых, непосредственное картирование может осуществляться различными программными средствами как в диалоговом, так и в пакетном режиме и, во-вторых, результаты картирования могут графически оформляться или использоваться в независимых программных средствах. Такова технологическая реальность коллективной работы в многопользовательском режиме. Это выдвигает задачу удобства подготовки, сохранения и передачи логически целостных наборов данных между независимыми программными средствами. Для обеспечения минимального выполнения этих требований мы пошли по пути создания в базе данных проектов, в рамках которых описываются наборы картируемых объектов, логическая структура характеризующих их наборов данных и непосредственно сохраняются сами подготовленные наборы.
Описание проекта содержит координаты прямоугольной рамки картирования, ссылки на используемые авторский и сводный вариант каталога и указание в рамках номенклатуры какого из районов описываются объекты данного проекта.
В проекте может быть описано множество поверхностей и объектов для картирования. Допустимыми являются четыре типа поверхностей:
1. Простые (указывается укороченное имя кров ли картируемого таксона, например – «БВ8») – в картировании участвуют разбивки, в которых есть имя данного таксона.
2. Дополненные аналогами — в картировании участвуют разбивки поименованного таксона и его возрастные аналоги, определяемые по условным абсолютным возрастам номенклатуры (например, кровля баженовской свиты –143 млн. лет и ее возрастных аналогов с тем же значением абсолютного возраста).
3. Составные границы описанных упорядоченным перечнем имен границ – фрагментов составной поверхности (например, поверхность доюрского основания описывается последовательностью укороченных имен «Р-Т» «КОРА ВЫВ» «ФУНДАМЕН», означающей, что по скважине на картирование приоритетно поступают разбивки по кровле промежуточного комплекса («Р-Т»), в случае их отсутствия берутся разбивки по кровле коры выветривания («КОРА ВЫВ») и только в случае отсутствия в скважине разбивок и по промежуточному комплексу и по коре выветривания в набор включается положение поверхности невыветрелого фундамента («ФУНДАМЕН»).
Картируемый объект проекта может содержать в своем описании одну или две поверхности: его кровлю и подошву. Каждая из них может быть простой, дополненной или составной границей.
С каждым из картируемых объектов проекта могут быть связаны один или несколько рабочих наборов данных для картирования. Рабочий набор представляет собой описатель его структуры и множество (последовательность) точек пересечения скважин и картируемого объекта с приписанными к ним атрибу тами (признаками объекта).
Реализовано несколько типов рабочих наборов данных для картирования. Их структура приводится в табл. 4.
Таблица 4.
Для каждого из них написана автономная каталогизированная на ядре СУБД утилита, осуществляющая заполнение рабочего набора данными, отражающими текущее состояние базы. Таким образом, рабочий набор представляет собой отражение, снимок базы данных. Он непротиворечив настолько, насколько согласованным было текущее содержимое базы и насколько точна логика проверки данных, прописанная в утилитах. Все это подвергается постоянной проверке и уточнению в ходе практической эксплуатации. Остается добавить, что с использованием этой технологии построено несколько тысяч карт от локальных до региональных на всю территорию ХМАО или Западно-Сибирской провинции.
Данные рабочего набора могут быть заполнены посредством активизации утилиты и запрошены в диалоговом режиме непосредственно из картирующих программ, например, программа «Medium» А.Г. Плавника или подготовлены также посредством активизации утилиты и выгружены в виде файлов из программ, написанных в «OracleForms». Пример фрагмента рабочего набора данных для картирования структурной поверхности приводится в табл. 5.
Таблица 5.
Второй тип массовых задач – построение моделей геологического пространства: сводных разрезов, профилей, трехмерных моделей. Действующие стандарты оформления пакетов геологической документации требуют обязательного представления сводных, типовых и (или) гипостратотипических геолого-геофизических разрезов и геологических профилей, фиксирующих принятую модель строения территории. Поэтому мы разработали и с 1994 года эксплуатируем специальную компьютерную технологию их построения, опирающуюся на существующую в ГУП ХМАО НАЦ РН интегрированную базу технико-экономических и геолого-геофизических данных и специализированную программу PLOTLOG (автор С.А. Предеин).
PLOTLOG — это программа построения геологических разрезов, схем корреляции, сводных геологических разрезов, схем опробования, блок диаграмм, планшетов интерпретации геолого-геофизической информации, профилей выравнивания и других геологических приложений. Программа позволяет визуализировать и корректировать разрез на экране компьютера и воспроизводить его твердую копию в любом вертикальном и горизонтальном масштабе на выводное устройство, поддерживаемое операционной средой MS WINDOWS 95–98 и WINDOWS NT.
Программа построения геологических разрезов оперирует с различными типами геологической информации:
* данными, задающими структурные взаимоотношения стратиграфических, литологических и тектонических границ, на основании которых производится построение структурного каркаса геологического разреза;
* производит учет структурных карт или сейсмических временных карт, представленных в числовом виде;
* данными, определяющими литологическую изменчивость изучаемого объекта;
* данными о характере насыщения изучаемого объекта флюидами;
* текстовой информацией, отображаемой на разрезе.
Все вышеперечисленные типы данных рассматриваются совместно для получения адекватной и непротиворечивой геологической модели.
Построение геологических моделей программой PLOTLOG осуществляется с использованием информации, выгружаемой в фиксированных форматах из литолого-стратиграфического раздела ИБД с помощью специализированного и достаточно сложного программного обеспечения.
Описания геологических профилей хранятся в ИБД в виде списка скважин, комплекса отображаемых геологических границ и ссылок на варианты корреляции и интерпретации ГИС.
Для построения профиля программой PLOTLOG из ИБД выгружается хранимая информация: набор используемых границ и описания скважин, а для каждой из скважин – глубины залегания границ, результаты послойной интерпретации ГИС, послойное описание литологической колонки и (при необходимости), привязанное к глубинам иерархическое описание взаимоотношения таксонов МСШ и местных литолого-стратиграфических шкал.
Рис. 3. Типичный локальный геологический разрез 1 – тонкие аргиллиты, 2 – аргиллиты битуминозные, 3 — глины, 4 – песчаники, 5 – нефтяные залежи.
Рис. 4. Типичная схема корреляции с литологической изменчивостью
Рис. 5. Фрагмент сводного геолого-геофизического разреза 1 – нижний мел К1, 2 – верхняя юра J3, 3 – средняя юра J2, 4 – нижняя юра J1, 5 – триас Т, 6 – палеозой PZ, 7 – песчаники, 8 – алевролиты, 9 – глины, 10 – аргиллиты, 11 – глины опоковидные, 12 – глины диатомовые, 13 – аргиллиты глауконитовые, 14 – аргиллиты битуминозные, 15 – известняк, 16 – базальты, 17 – глинистые сланцы
Всего в ИБД хранится несколько сотен таких описаний по построенным профилям. Примеры построения профилей на основе информации ИБД с использованием программы PLOTLOG часто приводились в публикациях «Вестника недропользователя». Построение сводных (типовых и гипостратотипических) геолого-геофизических разрезов в программе PLOTLOG аналогично построению схемы корреляции для одной искусственной скважины, которая составлена из фрагментов (интервалов) реальных поисково-разведочных скважин (рис. 5).
Интервалы эти берутся, как правило, между двумя регионально выдержанными границами и характеризуют весь вскрытый разрез от четвертичных отложений до фундамента. Такое пофрагментное описание сводного разреза хранится в ИБД и при необходимости выгружается в форматах программы PLOTLOG. Дополнительно выводятся привязанные к глубинам текстовые описания стратонов и отражающих сейсмических горизонтов. Для локальных участков территории построение сводных разрезов может осуществляться непосредственным описанием его интервалов в ИБД с последующей выгрузкой из нее данных для программы. Другая возможность построения сводного разреза — корректировка в программе PLOTLOG ранее построенных зональных гипостратотипов.
Для каждой зоны литолого-стратиграфического районирования на территории ХМАО были описаны в ИБД, подготовлены и построены в программе PLOTLOG типовые сводные геолого-геофизические разрезы (зональные гипостратотипы). Для переоценки потенциальных ресурсов УВ в верхнеюрском комплексе для всех эталонных участков на территории ХМАО по данной технологии было построено несколько десятков типовых разрезов и геологических профилей.
Хранение всей геолого-геофизической информации и описания профилей в ИБД избавляет пользователя от трудной и бессмысленной работы повторного сбора и подготовки данных. Весьма важным результатом такой технологии является и корректировка информации в ИБД, постоянное улучшение логической согласованности и целостности данных.
А это дает возможность многократно использовать как в программе PLOTLOG, так и в других приложениях, хранимую в ИБД информацию и верифицировать множество моделей геологического строения территории. Широкий спектр способов графического представления моделей от схем корреляции до блок-диаграмм позволяет интерпретатору-геологу надежно контролировать непротиворечивость рассматриваемой модели по всему комплексу доступной геолого-геофизической информации: структурным планам опорных сейсмических горизонтов, кривым и результатам интерпретации ГИС, литологическим колонкам и результатам испытания предположительно продуктивных пластов.
Таким образом, данные стратиграфического раздела ИБД ЦРН используются в производстве геологической документации, в частности, при подготовке лицензионных соглашений на недропользование, оценке ресурсов и проектировании геологоразведочных работ.
И, наконец, номенклатура литолого-стратиграфических объектов и каталог разбивок применяются для идентификации, поиска и упорядочения всей геолого-геофизической информации по стволу скважин.
Литература
1. Каталог стратиграфических разбивок / Под редакцией Н.Н. Ростовцева и А.А. Булынниковой
— Тр. ЗапСибНИГНИ.
— Вып. 67. -Тюмень.
— 1971.- 314 с.
2. Гришкевич В.Ф., Лагутина С.В. Информационно-функциональная модель автоматизированного каталога литолого-стратиграфических разбивок.//Математические методы прогнозирования нефтегазоносности в Западной Сибири/ Тр. ЗапСибНИГНИ.- 1988.
— Тюмень. -С. 131–140.
3. Стратиграфический кодекс.
— С.-Пб. -ВСЕГЕИ. -1992.-120 с.
4. Региональные стратиграфические схемы мезозойских отложений Западно-Сибирской равнины.
— Тюмень. – Изд. ЗапСибНИГНИ
— 1991.
5. Территория деятельности предприятия «Хантымансийскнефтегазгеология»/ Под редакцией А.М. Брехунцова, В.И. Шпильмана.- Тюмень
— Ханты-Мансийск.
— 1992.
— 96 с.
6. Карогодин Ю.Н.. Ершов С.В., Сафонов В.С и др. Приобская нефтеносная зона Западной Сибири: системно-литмологический аспект. Новосибирск. Изд-во СО РАН, НИЦ ОИГГМ.- 1996. -252 с.
7. Нежданов А.А., Огибенин В.В., Куренко М.И. и др. Региональная литмостратиграфическая схема мезозоя и кайнозоя Западной Сибири и основные закономерности размещения неантиклинальных ловушек углеводородов // Литмологические закономерности размещения резервуаров и залежей углеводородов. – Новосибирск.- Наука.
— 1990.-С. 80–108.
8. Шпильман В.И., Мясникова Г.П., Плавник Г.И. Выделение объектов анализа для прогнозирования и их предварительная обработка // Методика оценки прогнозных и перспективных запасов и обоснование их подсчетных параметров./Тр. ЗапСибНИГНИ.
— 1970. – Тюмень.
— Вып. 53.- С. 28–49.
9. Шпильман В.И., Гришкевич В.Ф. Характеристика истории тектонического развития резервуара //Применение математических методов при поисково- разведочных работах на нефть и газ в Западной Сибири/ Тр.ЗапСибНИГНИ.
— 1975.-Тюмень.
— Вып. 99.
— С.152–154.
10. Шпильман В.И., Мясникова Г.П., Плавник Г.И. и др. Метод количественной оценки потенциальных, прогнозных и перспективных ресурсов.
— ГУП НАЦ РН ХМАО.- Тюмень.
— 1999.- 99 с.
11. Каталог литолого-стратиграфических разбивок разрезов поисково-разведочных скважин. Ханты-Мансийский автономный округ //Под редакцией Гришкевича В.Ф, Теплякова Е.А. – Ханты-Мансийск.- 2000.-300с.
Автор: Болотова И.Л., Гришкевич В.Ф., Колосова Л.А., Лагутина С.В., Предеин С.А., Рыбьякова Е.В. (ГУП ХМАО НАЦ РН им. В.И.Шпильмана)