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

Автореферат диссертации по теме "Учет литологической макронеоднородности при цифровом геологическом моделировании месторождений нефти и газа"

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

ИСТОМИН СЕРГЕЙ БОРИСОВИЧ

УЧЕТ ЛИТОЛОГИЧЕСКОЙ МАКРОНЕОДНОРОДНОСТИ ПРИ ЦИФРОВОМ ГЕОЛОГИЧЕСКОМ МОДЕЛИРОВАНИИ МЕСТОРОЖДЕНИЙ НЕФТИ И ГАЗА

Специальность 25.00.10. — геофизика, геофизические методы поисков полезных ископаемых

АВТОРЕФЕРАТ

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

Москва 2009

СЮ34Э2 1

003492119

Работа выполнена на кафедре сейсмических и скважинных методов Российского государственного геологоразведочного университета и в ОАО «Центральная геофизическая экспедиция», г.Москва.

Научный руководитель:

Кандидат технических наук Мараев Игорь Алексеевич

Официальные оппоненты:

Доктор технических наук Денисов Сергей Борисович Кандидат геолого-минералогических наук Закревский Константин Евгеньевич

Ведущее предприятие:

РГУНГ им И.М. Губкина

Защита состоотся «И» марта 2010 г. в 15 часов на заседании Диссертационного совета Д.212.121.07 в Российском Государственном Геологоразведочном Университете им. Серго Орджоникидзе по адресу 117997, Москва, ул. Миклухо-Маклая, 23, РГГРУ. ауд. 6-38.

С диссертационной работой можно ознакомиться в библиотеке Российского Государственного Геологоразведочного Университета.

Автореферат разослан "И" февраля 2010 г.

Доктор технических наук, ученый секретарь диссертационного совета, профессор Л —т^А.Д. Каринский

Общая характеристика работы

Актуальность темы

Для прогноза степени вытеснения нефти и газа из продуктивных пластов, нефтеотдачи или коэффициента нефтеизвлечения (КИН) широко используются геолого-гидродинамические трехмерные цифровые модели. На основе этих моделей осуществляется оценка геологических и извлекаемых запасов углеводородов, обосновывается оптимальная стратегия доразведки и разработки продуктивных пластов. Создание и применение таких постоянно действующих геолого-технологических моделей (ПДГТМ) является одним из главных направлений повышения качества проектирования, управления и контроля разработки нефтяных и газонефтяных месторождений.

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

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

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

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

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

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

Основные задачи исследований

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

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

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

4. Алгоритмизация учета неоднородностей геологических объектов при моделировании на основе установленных взаимосвязей между критериями неоднородности и параметрами стохастического моделирования.

5. Разработка методики трехмерного геологического моделирования с учетом заданных критериев латеральной и вертикальной неоднородности с одновременным учетом результатов послойного моделирования и атрибутного анализа сейсмических данных.

Защищаемые положения

1. Необходимый набор критериев, достаточный для описания макронеоднородносги пластов при трехмерном геологическом моделировании, включает:

о Для вертикальной неоднородности: коэффициент песчанисгости и параметр расчлененности; о Для латеральной неоднородности: коэффициент распространения и параметр плотности ли-тологических границ.

Использование данного набора критериев обеспечивает отражение в модели геологической макронеоднородности пластов.

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

ронеоднородности пластов, что позволяет выполнять моделирование с учетом установленных взаимосвязей и получать предсказуемые результаты.

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

Научная новизна

1. Установлен информативный набор критериев геологической макронеоднородности для целей геолого-технологического моделирования.

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

3. Установлены количественные взаимосвязи между критериями макронеоднородности и параметрами стохастических индикаторных методов литофациального моделирования.

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

Практическая ценность

• Методика количественного анализа литологической макронеоднородности продуктивных пластов нефтяных и нефтегазовых месторождений и учета неоднородности при построении трехмерных геологических моделей является важным этапом в создании геолого-технологических моделей и используется в практике геологического моделирования в ОАО «ЦГЭ».

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

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

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

Фактический материал и личный вклад. Использованы материалы работ по подсчетам шасов и ТЭО КИН, построению сейсмогеологических моделей месторождений, выполненных в IAO «ЦГЭ» при непосредственном участии автора. С учетом проведенных исследований с ис-ользованием программного комплекса «DV-GEO» при непосредственном участии автора и под

его руководством определенных этапов в ОАО «ЦГЭ» были выполнены работы по построению и мониторингу геологических моделей более чем по двадцати месторождениям.

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

Результаты методических разработок опробованы и используются при выполнении построений трехмерных геолого-технологических моделей в отделении геоинформационных технологий ОАО «Центральная геофизическая экспедиция» г. Москва.

Объем и структура работы. Диссертация состоит из введения, четырех глав и заключения. Объем работы - 94 страницы текста, 16 рисунков, 5 таблиц. Список литературы содержит 57 наименований.

Благодарности. Диссертационная работа выполнена в период работы в ОАО «ЦГЭ», а также учебы в аспирантуре Российского государственного геологоразведочного университета. Основной объем работы выполнен в ОАО «ЦГЭ», где соискатель работает в Отделении Геоинформационных технологий - главным специалистом по геологии и геологическому моделированию.

Автор выражает благодарность начальнику ОАО «ЦГЭ» и научному руководителю работ по разработке программного комплекса «DV-GEO» д.т.н. Кашику A.C., глубокую признательность Дьяконовой Т.Ф. и Билибину С.И., без которых данная работа не могла состояться; своим научным руководителям профессору |Генри Николаевичу Боганик>1 и доценту Игорю Алексеевичу Ма-раеву.

Большое влияние на направление исследований в разное время оказали совместная работа и творческие контакты с Денисовым С.Б., Гавриловой Е.В, Кожевниковым Д.А., Коваленко К.В. Автор приносит глубокую благодарность им и всем, кто содействовал выполнению этой работы.

Содержание работы

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

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

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

Г.МЗолоевой. Классификацию геологических тел с выделением различных уровней структурной организации проводили Косыгин Ю.А., Дементьев Л.Ф., Токарев М.А. В соответствии со структурной организацией объектов могут быть выделены типы неоднородности, представленные в таблице 1.1.

Таблица 1. Иерархия типов геологической неоднородности

Номер уровня Тип неоднородности Иерархический уровень неоднородности

1 Ультрамикронеоднородность Минеральное зерно

2 Микронеоднородность Конгломерат минеральных зерен, образующих элементарный объем

3 Мезонеоднородность Условное геологическое тело единого вещественного состава, разделенное граничными свойствами пород (коллектор/неколлектор, высоко-, средне- и низкопроницаемые песчаники и т. п.)

4 Макронеоднородность Отдельные линзы или пропластки, сложенные единым типом пород, взаимное расположение этих геологических тел внутри пласта

5 Метанеоднородность Система из пластов или более крупных частей (зон)

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

Выделяются два основных проявления макронеоднородности: расчлененность продуктивного горизонта на разобщенные пласты и прерывистость отдельных пластов и прослоев по площади. Изучение макронеоднородности проводится в соответствии с данными видами ее проявления: по вертикали (по толщине горизонта) и по простиранию пластов (по площади).

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

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

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

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

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

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

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

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

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

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

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

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

Во второй главе «Критерии оценки литологической макронеоднородности» проводится обзор критериев, применяемых для оценки литологической макронеоднородности. Показывается взаимосвязь между показателями неоднородности. Приводятся способы оценки литологической макронеодиородности по исходным данным. Выводятся аналитические зависимости между показателями макронеоднородности и параметрами индикаторного стохастического моделирования.

Обзор критериев оценки литологической макронеоднородности. Всю совокупность показателей литологической макронеоднородности можно разделить на показатели вертикальной и показатели латеральной макронеоднородности.

Для оценки вертикальной неоднородности по данным скважин используются:

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

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

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

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

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

Коэффициент выклинивания Кл, показывающий долю выклинивающихся прослоев коллекторов. Он определяется отношением толщины выклинивающихся слоев коллектора к эффективной толщине. В отсутствии выклинивания Кл = О, при выклинивании всех прослоев Кл = 1.

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

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

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

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

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

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

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

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

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

К

Релс = -Ш-, где Нее - средняя толщина прослоев коллектора. То есть коэффициент песчанистогости, параметр расчлененности и средняя толщина прослоев коллектора - взаимосвязанные показатели вертикальной литологической макронеодиородносги пласта.

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

4^

ществует аналогичная связь: Рпл =--.

Ксг

Две пары параметров — два параметра вертикальной литологической макронеоднородносги и два параметра латеральной литологической макронеоднородности - наиболее полно описывают количественно и геометрически содержание в объеме пласта тел коллектора. Данные параметры являются базовыми для расчета большинства других критериев оценки литологической макронеоднородносги.

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

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

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

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

ном плотной сеткой скважин, можно приближенно оценить этот параметр аналитически: Рпл = 8,2347-^ К,.

Таблица 2. Параметры макронеоднородности для террнгенных пластов месторожде-_нпй Западной Сибири_

нго Условия седиментации Вертикальная макронеоднородность Латеральная макронеоднородность Средняя общая толщина геол. объекта, м Кол-во скв-н в анализе Площадь иссл-я, км' Средняя плотность сетки скважин, га/скв

Кпес, ср (м/м) Ррас, ср (1/м) к» верт К», ср (м2/м2) Рпл, ср (м(м2) Ки, лат

Среднеобсхая открытый шельф, каналы 0.43 0.20 0.47 0.41 0.001573 3.79 26.0 300 48.4 16.1

Среднеобсхая Мелководный шельф, каналы 0.76 0,16 0.21 0.77 0.001272 1.65 40.6 137 23.2 16.9

Среднеобсхая открытый шельф, бары, барьерные остров* 0 28 0.16 0.57 0.33 0.001486 456 10.8 138 23.2 16.8

Среднеобсхая открытый шельф, каналы 0.39 0.13 0.33 0.45 0.001443 3.23 25.4 399 56.7 14.2

Приуральская мелководные, намывные песчаные тела 0.50 0.13 0.26 0.63 0.000761 1.21 10.1 117 38.8 33.2

Приуральсхая Континентальные 0.47 0.22 0.47 0.48 0.001322 2.74 16.9 117 38.8 33.2

Приуральская континентальные 0.49 0.17 0.35 0.47 0.001105 2.34 24.3 117 38.8 33.2

Среднеобсхая мелководно-морские 0.48 0.08 0.17 0.66 0.000762 1.16 24.0 326 61.3 18.8

Среднеобсхая мелководно-морсхие 0.20 0.06 0.30 0.20 0.000689 3.51 14.0 325 61.3 18.9

Среднеобсхая мелководно-иорсхие 0.31 0.14 0.45 0.32 0.001393 4.31 23.0 325 61.3 18.9

Среднеобсхая континентальные 0.16 0.08 0.50 0.25 0.001034 4.11 28.7 140 44.1 31.5

Среднеобсхая мелководно-морские 0.29 0.07 024 0.45 0.000907 2.03 21.8 37 11.8 31.9

Среднеобсхая мелководно-морские 0.28 о.ов 0.21 0.48 0.000714 1.47 29.9 37 11.8 31.9

Среднеобсхая мелководно-морские 0,34 0.10 0.29 0.57 0.001375 2.40 25.6 31 6.2 20.0

Среднеобсхая мелководно-морские 0.28 0.09 0.32 0.41 0.001406 3.45 17.2 28 6.2 22.1

Фролоесхая мелководный шельф 0.49 0.10 0.20 0.52 0 001301 2.52 55.3 174 24.2 13.9

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

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

теральные (горизонтальные) вариограммы, походят наилучшим образом. Рассмотрим связь варио-грамм с параметрами макронеоднородности.

Исследователями из Университета Калифорнии (США) Б.Р.Саг1е и G.E.Fogg, занимавшимися проблемой взаимосвязи геологических неоднородностей с параметрами индикаторной вариограммы, были выведены следующие зависимости для разрезов скважин. Порог вариограммы напрямую связан с долей пространства, занимаемой одной из фаций (Кпес) и равен: с = к пес

Средний размер тел фации (Нср) равен произведению песчанистости на тангенс начального угла наклона вариограммы: НСР = КПЕС ■tg{a).

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

т = С-\~-~) при 02А5а у2-а 2-а )

у(Ы) = С при А > а

где у(К) - дисперсия, откладываемая на вариограмме по оси ординат; А - расстояние, откладываемое на вариограмме по оси абсцисс; а - радиус вариограммы; С - порог вариограммы.

Тангенс начального угла вариограммы равен отношению расстояния А к дисперсии у{И)

при стремящемся к нулю расстоянии А: tg(a) =- при А -> 0. Подставляя это тождество в

Л А)

уравнение вариограммы получим: tg{a) - . Из полученного уравнения мы можем получить

зависимость радиуса вариограммы от порога вариограммы и ее начального угла: а = 1.5 • С ■ tg{a).

Теперь, имея в наличии взаимосвязь всех параметров вариограммы и связи порога и начального угла вариограммы с показателями неоднородности, мы можем показать основные связи между вариограммой и вертикальной макронеоднородностью:

С = КПЕС ■ (1 -КПЕС), а = 1,5-С. = .

пес рас

где КПЕС - коэффициент песчанистости; Нср - средняя толщина прослоев коллектора, ¡'гас - параметр расчлененности; С и а - порог и радиус вариограммы.

Для трансформации критериев оценки латеральной неоднородности в параметры индикаторной вариограммы следует принять во внимание, что долю коллектора (аналог песчанистости)

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

пл

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

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

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

Два основных класса методов для построения трехмерных стохастических моделей дискретных параметров, существующих в настоящее время, - это методы на основе объектов и пиксельные (индикаторные) методы.

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

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

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

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

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

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

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

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

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

Таблица 3. Результаты оценки неоднородности пластов при различной плотности сетки скважин _

Возраст отложений пласта * * Верт. неоднородность Лат. неоднородность Верт. варио-грамма Лат. варио-грамма

2 2 Секта Условия с диментац* Плот-ть сетки, га/скв Плот-тъ сетки, скв/км3 Клее, ср (х/м) Ррас, сред (1/м) Км, верт К», ср (м/м ) Рпл,ср (км/ юг) Км, лат Сверт ср 1?верт ср(м) Сгор, ср Ягор, ср (м)

13.9 7.2 0.49 0.10 0.20 0.52 1.301 252 0.23 4.08 0.21 1079

и * 23.5 43 0.47 0.10 0.21 0.50 1.071 2.12 0.23 4.08 021 1278

3 1 310 32 0.47 0.10 0.21 0.50 0.999 2.00 0.23 3.95 0.21 1379

в о 80.7 1.2 0.47 0.10 0.21 0.54 0.627 1.16 0.24 4.14 0.21 2135

1 6 X к 8| 33.2 3.0 0.50 0.13 0.26 0.63 0.761 1.21 0.20 2.64 0.17 1405

2 ¿с о в а 74.6 1.3 0.55 0.13 0.24 0.61 0.359 0.59 0.21 2.81 0.15 2491

а ч !в * а 138.6 07 0.52 0.14 0.27 0.62 0.281 0.45 0.21 2.59 0.15 4011

332 3.0 0.47 0.22 0.47 0.48 1.322 2.74 0.22 1.63 0.23 1076

3 <N1 _ ! 74.6 1.3 0.46 0.24 0.52 0.44 0.695 1.57 022 1.45 0.23 1997

¥ ¡3 £ 138.6 0.7 0.46 0.22 0.48 0.46 0.514 1.12 0.22 1.53 0.22 2689

К Ф X 332 3.0 0.49 0.17 035 0.47 1.105 2.34 0.20 2.11 0.23 1315

4 сч к 74.6 1.3 046 0.15 0.33 0.44 0.554 126 0.21 2.35 0.20 2279

138.6 0.7 0.51 0.17 0.33 0.47 0.430 0.91 0.21 2.10 0.20 3195

к я X 16.1 6.2 0.43 0.20 0.47 0.41 1.573 3.79 021 1.71 0.21 893

О) в 235 4.3 0.43 0.20 0.47 0.43 1.041 244 0.21 1.71 0.20 1258

5 ее а 31.6 3.2 0.45 0.20 0.44 0.43 0.895 2.07 021 1.69 0.20 1478

* п 69.1 1.4 0.47 0.19 0.40 0.44 0.644 1.46 0.21 1.72 0.20 2028

а> ч 201.7 0.5 0.38 0.19 0.50 0.48 0.408 0.84 0.21 1.84 022 3438

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

1. Количественная оценка показателей макронеоднородности на основе исходных геолого-геофизических материалов и результатов их интерпретации.

1.1. Оценка показателей вертикальной макронеоднородности проводится с использованием результатов интерпретации материалов ГИС.

1.2. Оценка показателей латеральной неоднородности. Результаты такой оценки могут быть получены (ранжировано в порядке приоритета):

1.2.1. На участке опытно-эксплуатационного разбуривания залежи.

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

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

1.2.4. Аналитическим путем на основе зависимости между параметрами вертикальной и латеральной макронеоднородности.

2. С учетом полученных результатов оценки макронеоднородности рассчитываются параметры для вариограммного моделирования и строятся вертикальная и латеральная вариограммы.

3. Выполняется распределение коллекторов в объеме геологической модели с использованием индикаторных методов. Исходными данными для моделирования являются:

3.1. результаты интерпретации материалов ГИС в скважинах;

3.2. вертикальная и латеральная вариограммы;

3.3. ГСР по коэффициенту распространения;

3.4. тревдовые карты эффективных толщин (песчанистости), полученные в результате интерпретации сейсмических материалов или в результате построения концептуальной седимента-ционной модели.

4. Производится верификация построенной литологической модели резервуара. Первоочередным критерием верности построенной модели является корректность учета скважинных данных. Далее выполняется построение из куба «коллектор - неколлектор» карт песчанистости и гео-лого-статисгических разрезов по коэффициенту распространения и параметру плотности лито-логических границ в каждом слое модели. Полученные карты и ГСР сопоставляются с трендо-выми. Расхождение по ГСР для каждого слоя модели не должно превышать 5%.

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

Верификация геологических моделей. Широкое использование геолого-технологических моделей вызывает необходимость оценивать достоверность получаемых результатов и их соответствие исходной геолого-геофизической информации. В методических рекомендациях по проектированию за 2007г в разделе «Построение литологических моделей залежей и распределения фильтрационно-емкосгных свойств пластов» говорится о том, что для оценки достоверности кубов литологии используются построенные по этим кубам карты эффективных толщин, пористости и проницаемости, которые сравниваются с двухмерными картами, построенными по данным скважин. Отклонения в определении параметров не должны превышать ±5%.

Построение и сопоставление геолого-сгатистических разрезов (ГСР) по скважинам в интервале моделируемого пласта и по кубу литологии является одним из способов проверки правильности литологического моделирования. Этим способом осуществляется проверка учета латеральной неоднородности, поскольку ГСР отражает долю коллектора в слоях пласта (геологической моде-

ли). Фактически, значения каждого из отсчетов ГСР являются показателем коэффициента распространенности коллектора в отдельных слоях пласта. Однако даже полное совпадение коэффициента распространенности в модели со скважинными данными не гарантирует воспроизведения латеральной литологической неоднородности, которая, как показано, определяется двумя показателями: коэффициентом распространения и параметром плотности литологических границ.

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

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

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

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

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

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

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

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

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

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

Условия седиментации осадочных образований Западной Сибири наиболее системно рассматривались в работах А.Э.Конторовича, Ф.Г.Гурари, В.С.Суркова, М.В.Коржа и А.Г.Мухер и других исследователей. Особенности строения пластов Самотлорского месторождения изучались Денисовым С.Б. и др.

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

повышения градиента карты общих толщин на северо-западе территории совпадает с палеолинией бровки шельфа.

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

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

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

При оценке параметров литологической неоднородности по выделенным зонам наблюдается снижение неоднородности при приближении к источнику сноса материала (табл. 4).

Как сказано выше, основная Самотлорская залежь пласта БСв1"3 расположена в пределах третьей и четвертой зон. Оценка распределения удельных дебитов скважин основной Самотлор-ской залежи по фациальным зонам показывает, что количество относительно низкодебитных скважин в третьей зоне значительно превышает количество таких скважин в четвертой зоне. Это хорошо согласуется с распределением показателей макронеоднородности строения пласта.

Таблица 4. Показатели вертикальной макронеоднородности пласта БВ»'"3 Самотлорского региона по фациальным зонам

Показатели макронеоднородносги Зона 1 Зона 2 ЗонаЗ Зона 4

Коэффициент песчанистости - Клее, м/м 0,27 0,41 0,58 0,59

Параметр расчлененности - Ррас, 1/м 0,068 0,079 0,082 0,010

Коэффициент вертикальной макронеоднородности -Км,верт 0,252 0,193 0,141 0,017

Коэффициент распространенности - Кб, м7м^ 0Д7 0,41 0,58 0,59

Параметр плотности -Рпл, км/км'1 0,560 0,651 0,675 0,083

Коэффициент латеральной макронеоднородности -Км,лат 2,07 1,59 1,16 0,14

Пласт АВ1''3. Образование отложений пласта АВ1'"2 происходило в условиях постепенного углубления морского бассейна, сопровождавшегося снижением динамики процессов седиментации. Такие застойные условия и привели к образованию пород с характерной текстурой типа «рябчик», давшей название всему пласту АВ1'"2. Для пласта характерно частое тонкое переслаивание песчаников, алевролитов и глин, пятнистость и сероцветностъ. Пласт отличается от нижележащих пластов АВ ухудшенными ФЕС и площадным развитием глинистого типа разреза.

При переносе кластического материала с востока на запад через такие формы рельефа происходило в первую очередь осаждение и накопление наиболее крупнозернистых обломков. Этим фактом обуславливается измельчение обломочных фракций в интервале пласта в западном направлении. Конечным результатом процесса осадкообразования явилось постепенное замещение песчаников глинистыми разностями при продвижении на запад рассматриваемой территории. Описанные закономерности нашли отражение в выделенных фациальных областях АВ11"2.

При оценке параметров макронеоднородности по выделенным зонам наблюдается уменьшение неоднородности при приближении к источнику сноса терригенного материала (табл. 5).

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

Таблица 5. Показатели вертикальной макронеоднородностн пласта АВ!1"1 Самотлорского региона по фациальным зонам

Показатели макронеоднородности Зона 1 Зона 2 ЗонаЗ Зона 4

Коэффициент песчанистости - Кпес, м/м 0,14 0,16 0,52 0,85

Параметр расчлененности-Ррас, 1/м 0,096 0,109 0,197 0,111

Коэффициент вертикальной макронеоднородности - Км,верт 0,686 0,681 0,379 0,131

Коэффициент распространенности - Кв, м7м* 0,14 0,16 0,52 0,85

Параметр плотности - Рпл, км/км^ 0,791 0,898 1,622 0,914

Коэффициент латеральной макронеоднородности -Км,лат 5,65 5,61 3,12 1,08

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

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

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

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

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

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

Основные работы автора по теме диссертации:

1. Учет литологической неоднородности терригенных пластов при построении геологических моделей. Журнал «Геофизика», №4,2007г, с.186-190.

2. Особенности современного подхода к подсчету геологических запасов месторождений нефти и газа на основе трехмерных моделей. Журнал «Нефтяное хозяйство», № 10, 2006г. (Совместно с С.И. Билибиным, Т.Ф. Дьяконовой, Е.В. Гавриловой, Т.Г. Исаковой, Е.А. Юкановой).

3. Прогнозирование распространения коллекторов отложений вогулкинской толщи с использованием палеотектонического анализа на примере пласта П2 Тальникового месторождения. Журнал «Каротажник», №3-4 (116-117), 2004г, с.246-249.

4. Некоторые практические аспекты анализа межфлюидных контактов при геологическом моделировании залежей нефти и газа. Журнал «Геофизика», №4, 2007г, с.201-206. (Совместно с Кавуном М.М., Степановым A.B.).

5. Трехмерная геологическая модель - необходимый и обязательный этап изучения нефтегазового месторождения. «Недропользование XXI век» №4 Москва 2007г. Авторы: С.И.Билибин -к.т.н., Т.Ф.Дьяконова-д.г.м.н., Т.Г.Исакова, С.Б.Истомин., Е.А.Юканова. с. 38- 42.

6. Необходимость использования критериев неоднородности при геологическом моделировании неравномерно разбуренных площадей. Тезисы докладов VII Международной НПК «Геомодель-2005», г. Геленджик, 11-17 сентября, 2005г.

7. Особенности современного подхода к подсчету геологических запасов многопластовых месторождений нефти и газа на основе трехмерных моделей. V Международный технологический симпозиум «Новые ресурсосберегающие технологии недропользования и повышения нефтеотдачи» - тезисы докладов. 21-23 марта 2006 г. (в соавторстве с С.И. Билибиным, Т.Ф. Дьяконовой, Е.В. Гавриловой, Т.Г. Исаковой, Е.А. Юкановой).

8. Прогнозирование распространения коллекторов отложений вогулкинской толщи с использованием палеотектонического анализа на примере пласта П2 Тальникового месторождения. НПК «Геофизические и нефтепромысловые методы исследования скважин в комплексе с сейсморазведкой для построения и сопровождения геологических моделей залежей нефти и газа» - тезисы докладов, Москва, ОАО «ЦГЭ», 2004г.

9. Некоторые практические аспекты анализа межфлюидных контактов при геологическом моделировании залежей нефти и газа. Тезисы докладов IX Международной НПК «Геомодель-2007», г. Геленджик, 15-20 сентября, 2007г. (в соавторстве с Кавуном М.М., Степановым A.B.).

10. Выделение тектонических нарушений в волновом поле и их роль в формировании залежей углеводородов на Талинской площади (Красноленинский свод). Тезисы докладов IX Международной НПК «Геомодель-2007», г. Геленджик, 15-20 сентября, 2007г. (Совместно с Бричиковой М.П.).

11. Методика выявления тектонических нарушений по данным сейсморазведки 3D и их роль в формировании залежей углеводородов на Талинской площади. VI Международная научно-практическая конференция «ГЕОФИЗИКА 2007» - тезисы докладов, г. Санкт-Петербург, 15 октября, 2007г. (Совместно с Бричиковой М.П.).

12. Экспресс оценка разреза в скважинах статистическими методами с целью фациального зонирования терригенных пластов. «XVin Губкинские чтения» - тезисы докладов, г.Москва, 2009г.

Подписано к печати 10.02.10 Формат 60x90 1/16. объем 1 пл. Тираж 100 экз. заказ № 38

Отпечатано в ОАО «Центральная Геофизическая Экспедиция» Россия, 123289,Москва,ул.Народного Ополчения,38/3

Содержание диссертации, кандидата технических наук, Истомин, Сергей Борисович

Список рисунков.

Список таблиц.

Введение.

Глава 1. Современное состояние проблемы учета литологической макронеоднородности при построении геологических моделей.

1.1. Виды геологической неоднородности.

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

1.3. Способы изучения и учета литологической макронеоднородности.

1.4. Учет неоднородности при создании цифровых трехмерных геологических моделей.

Глава 2. Критерии оценки литологической макронеоднородности

2.1. Обзор критериев оценки литологической неоднородности

2.2. Взаимосвязь ме>еду показателями неоднородности и выбор показателей для геологического моделирования.

2.3. Оценка неоднородности по исходным данным.

2.4. Связь показателей неоднородности с параметрами стохастического моделирования.

Глава 3. Построение трехмерных компьютерных геолого-технологических моделей с учетом заданных показателей неоднородности.

3.1. Обзор методик построения трехмерных геологических моделей при литологическом моделировании.

3.2. Алгоритм выполнения построений трехмерной литологической модели с учетом заданных критериев макронеоднородности.

3.3. Влияние недоучета неоднородности при построении геологических моделей.

3.4. Верификация геологических моделей.

Глава 4. Влияние фациальных условий на литологическую макронеоднородность.

4.1. Фациальные предпосылки проявления неоднородности

4.2. Анализ фациального строения продуктивных пластов Самотлорского месторождения и оценка их макронеоднородности.

Введение Диссертация по наукам о земле, на тему "Учет литологической макронеоднородности при цифровом геологическом моделировании месторождений нефти и газа"

Актуальность темы

Для прогноза степени вытеснения нефти и газа из продуктивных пластов, нефтеотдачи или коэффициента нефтеизвлечения (КИН) в настоящее время широко используются геолого-гидродинамические трехмерные цифровые модели (геолого-технологические модели). На основе этих моделей осуществляется оценка геологических и извлекаемых запасов углеводородов, обосновывается оптимальная стратегия доразведки и разработки продуктивных пластов. Создание и применение таких постоянно действующих геолого-технологических моделей (ПДГТМ) является одним из главных направлений повышения качества проектирования, управления и контроля разработки нефтяных и газонефтяных месторождений.

Актуальность данного направления отражена в регламенте по созданию ПДГТМ нефтяных и газонефтяных месторождений, разработанном под эгидой Минтопэнерго РФ широким кругом специалистов (2000г.). Этот регламент является основным документом, описывающим требования, предъявляемые к компьютерным моделям нефтяных и газонефтяных пластов. К моменту написания данной работы в этом документе не были четко отражены требования к учету геологической неоднородности в процессе моделирования. Поэтому разработка методических и технологических решений для учета неоднородности литологического строения пластов при построении трехмерных моделей является в настоящее время одной из наиболее актуальных задач моделирования нефтяных и газовых месторождений.

Залежи многих месторождений нефти и газа приурочены к литологически неоднородным пластам. Неоднородность, выражающаяся в изменчивости формы залегания, минерального состава и фильтрационно-емкостных свойств коллекторов, 5 свойственна в той или иной мере всем продуктивным пластам. Цитологическая невыдержанность пластов даже на незначительных расстояниях, замещение коллекторов глинистыми породами, наличие слияний ме>еду прослоями коллекторов создают большие трудности при их разработке. Недоучет неоднородного строения продуктивных пластов приводит к существенным ошибкам при оценке запасов УВ и добычных возможностей залежей.

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

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

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

Цель работы

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

Основные задачи работы

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

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

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

4. Алгоритмизация учета неоднородностей геологических объектов при моделировании на основе установленных взаимосвязей между критериями неоднородности и параметрами стохастического моделирования.

5. Разработка методики трехмерного геологического моделирования с учетом заданных критериев латеральной и вертикальной неоднородности с одновременным учетом результатов послойного моделирования и атрибутного анализа сейсмических данных.

Защищаемые положения

1. Необходимый набор критериев, достаточный для описания макронеоднородности пластов при трехмерном геологическом моделировании, включает:

• Для вертикальной неоднородности: коэффициент песчанистости и параметр расчлененности;

• Для латеральной неоднородности: коэффициент распространения и параметр плотности литологических границ.

Использование данного набора критериев обеспечивает отражение в модели геологической макронеоднородности пластов.

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

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

Научная новизна

1. Установлен информативный набор критериев геологической макронеоднородности для целей геолого-технологического моделирования.

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

3. Установлены количественные взаимосвязи между критериями макронеоднородности и параметрами стохастических индикаторных методов литофациального моделирования.

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

Основные результаты

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

• Предложены алгоритмы оценки и учета при построении трехмерных геологических моделей показателей литологической макронеоднородности.

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

Практическая ценность работы

Методика количественного анализа литологической макронеоднородности продуктивных пла-стов нефтяных и нефтегазовых месторождений и учета неоднородности при построении трех-мерных геологических моделей является важным этапом в создании геолого-технологических моделей и используется в практике геологического моделирования в ОАО «ЦГЭ».

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

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

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

Апробация результатов

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

Результаты методических разработок опробованы и используются при выполнении построений трехмерных геолого-технологических моделей в отделении геоинформационных технологий ОАО «Центральная геофизическая экспедиция» г. Москва. На основе трехмерных геологических моделей, построенных с использованием алгоритмов учета литологической неоднородности, выполняются подсчеты геологических и извлекаемых запасов нефти и газа. Обосновываются коэффициенты извлечения нефти.

Благодарности

Диссертационная работа выполнена в период работы в ОАО «ЦГЭ», а также учебы в аспирантуре Российского государственного геологоразведочного университета. Основной объем работы выполнен в ОАО «ЦГЭ», где соискатель работает в Отделении Геоинформационных технологий - главным специалистом по геологии и геологическому моделированию.

Автор выражает благодарность начальнику ОАО «ЦГЭ» и научному руководителю работ по разработке программного комплекса «DV-GEO» д.т.н. Кашику A.C., глубокую признательность Дьяконовой Т.Ф. и Билибину С.И., без которых данная работа не могла состояться; своим научным руководителям профессору ¡Генри Николаевичу

Боганику] и доценту Игорю Алексеевичу Мараеву.

Большое влияние на направление исследований в разное время оказали совместная работа и творческие контакты с Денисовым С.Б., Гавриловой Е.В, Кожевниковым Д.А., Коваленко К.В., Шаховым П.А. Автор приносит глубокую благодарность им и всем, кто содействовал выполнению этой работы.

1. СОВРЕМЕННОЕ СОСТОЯНИЕ ПРОБЛЕМЫ УЧЕТА ЛИТОЛОГИЧЕСКОЙ МАКРОНЕОДНОРОДНОСТИ ПРИ ПОСТРОЕНИИ ГЕОЛОГИЧЕСКИХ МОДЕЛЕЙ

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

B.В.Воинова, Ш.К.Гиматудинова, Л.Ф.Дементьева, А.И. Деминой,

C.Б.Денисова, Г.М.Золоевой, М.М.Ивановой, И.М.Климушина, М.Л.Лейбина, В.Д.Лысенко, В.П.Меркулова, З.К.Рябининой, Е.И.Семина, В.В.Стасенкова, М.А.Токарева, О.С.Черновой, И.П.Чоловского, и др.

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

ЗАКЛЮЧЕНИЕ

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

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

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

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

86

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

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

1. Борисов Ю.П., Воинов В.В., Рябинина З.К. Влияние неоднородности пластов на разработку нефтяных месторождений. М., Недра, 1970.

2. Боровко H.H. Статистический анализ пространственных геологических закономерностей. Недра, Л.: 1971

3. Бриллиант Л.С., Аржиловский A.B., Цой В.Е. Экспертиза геологических и гидродинамических моделей основа качества исполнения проектной документации. Вестник российской академии естественных наук. - 2006. - Томб, №6.

4. Булыгин Д.В., Булыгин Д.Я. Геология и имитация разработки залежей нефти. 1996. - М.: Недра. - 382с.

5. Гавура В.Е. Геология и разработка нефтяных и газонефтяных месторождений. 1995. - М.: ВНИИОЭНГ. -496с.

6. Дементьев Л.Ф., Шурубор Ю.В., Азаматов В.И. и др. Оценка промышленных запасов нефти, газа и конденсата. М.: Недра, 1981.

7. Дементьев Л.Ф. Математические методы и ЭВМ в нефтегазовой геологии: Учеб. пособие для вузов. М.: Недра, 1983.

8. Демина А.И. Методика анализа прерывистости продуктивных пластов на основе трехмерных геологических моделей. // Газовая промышленность. 2006. - №12. - С.64-66.

9. Дьяконова Т.Ф., Билибин С.И., Денисов С.Б. Прогноз параметров коллекторов по данным комплексной интерпретации ЗД сейсморазведки и ГИС при построении цифровых геологических моделей. / Нефтяное хозяйство. 2000. - №10. - с.49-56.

10. Дюбрул Оливье. Использование геостатистики для включения в геологическую модель сейсмических данных. М.: SEG, EAGE, 2005.

11. Закревский К.Е., Геологическое 3D моделирование, М.: ООО «ИПЦ Маска», 2009.

12. Закревский К.Е., Майсюк Д.М., Сыртланов В.Р., Оценка качества 3D моделей, М.: ООО «ИПЦ Маска», 2008.

13. Золоева Г.М., Денисов С.Б, Билибин СИ. Геолого-геофизическое моделирование залежей нефти и газа: Учебное пособие. М.: ФГУП Изд-во «Нефть и газ» РГУ нефти и газа им. И.М. Губкина, 2005. - 172 с.

14. Золоева Г.М. Оценка неоднородности и прогноз нефтеизвлечения по ГИС. М.: Недра, 1995. - 212 с.

15. Иванова M.M., Дементьев Л.Ф., Чоловский И.П. Нефтегазопромысловая геология и геологические основы разработки месторождений нефти и газа. М., Недра, 1992.

16. Изучение геологической неоднородности продуктивных пластов / В.В. Воинов, М.Л. Лейбин, Е.И. Семин, З.А. Хрусталева. Науч.-техн. сб. по добыче нефти. М.: Недра, 1966. - № 14. - С. 3-7.

17. Истомин С.Б., Необходимость использования критериев неоднородности при геологическом моделировании неравномерно разбуренных площадей. VII Международная НПК «Геомодель-2005», г. Геленджик, 11-17 сентября, 2005г.

18. Истомин С.Б., Учет литологической неоднородности терригенных пластов при построении геологических моделей. Журнал «Геофизика», №4, 2007г, с. 186-190.

19. Истомин С.Б., Прогнозирование распространения коллекторов отложений вогулкинской толщи с использованием палеотектонического анализа на примере пласта П2 Тальникового месторождения. Журнал «Каротажник», №3-4 (116-117), 2004г, с.246-249.

20. Кавун М.М., Степанов A.B., Истомин С.Б., Некоторые практические аспекты анализа межфлюидных контактов при геологическом моделировании залежей нефти и газа. Журнал «Геофизика», №4, 2007г, с.201-206.

21. Каменский Г.Н., Корчебоков H.A., Разин К. И. Движение подземных вод в неоднородных пластах. М., Л.: Об'единенное научно-техническое издательство, 1935. - 168 с.

22. Каналин В.Г., Вагин С.Б., Токарев М.А., Ланчаков Г.А., Тимофеев В.А. Нефтегазопромысловая геология и гидрогеология. 1997. - М.: ОАО «Издательство «Недра». - 366с.

23. Кашик A.C., Билибин С.И., Гогоненков Г.Н., Кириллов С.А., Новые технологии при построении цифровых геологических моделейместорождений углеводородов / Технологии ТЭК. 2003. - № 3. - стр. 12-17

24. Корж М.В. Палеогеографические критерии нефтегазоносности юры Западной Сибири. Наука. М., 1978. 135 с.

25. Косыгин Ю.А. Основы тектоники. М.: Недра, 1974.

26. Лысенко В.Д., Мухарский Э.Д., Хамзин Р.Г. О неоднородности продуктивных пластов // Тр. ТатНИИ. Л.: Недра. - Вып. 6 -1964.

27. Медведев Е.А., Пьянков В.Н. Решение задач построения геолого-технологических моделей с использованием программного комплекса «БАСПРО Оптима». II Нефтяное хозяйство. 2005. - №10. - с.51-53.

28. Меркулов В.П., Чернова О.С. Закономерности пространственной изменчивости фильтрационно-емкостных параметров верхнеюрских коллекторов в связи с условиями седиментации / Сб. Научн. тр. школы-семинара НОЦ «ЮКОС-НОВОСИБИРСК», 2002. С. 182-188.

29. Методические рекомендации по проектированию разработки нефтяных и газонефтяных месторождений. 2007. - М.: Министерство природных ресурсов Российской Федерации.

30. Миллер Р. и Дж.Кан. Статистический анализ в геологических науках. Изд-во «МИР», М.: 1965.

31. Обстановки осадконакопления и фации: В 2-х т.: Пер. с англ./Под ред. Х.Рединга. М.: Мир, 1990. - 352 с.

32. Особенности современного подхода к подсчету геологических запасов месторождений нефти и газа на основе трехмерных моделей. С.И.Билибин, Т.Ф.Дьяконова, Е.В.Гаврилова, Т.Г.Исакова, С.Б.Истомин, ЕАЮканова. Журнал «Нефтяное хозяйство», № 10, 2006г.

33. Нестеров В.Н., Харахинов В.В., Семянов А.А., Шленкин С.И., Глебов А.Ф., Геологическая доразведка нефтяных месторождений Нижневартовского Приобья, Москва, Научный мир, 2006г

34. Повышение эффективности разработки нефтяных месторождений Башкирии / Э.М. Халимов, М.М. Саттаров, Ю.З. Зайнетдинов, М.Н. Галлямов. Уфа: Башк. книжн. изд-во, 1972.

35. Подсчет запасов нефти, газа, конденсата и содержащихся в них компонентов. Справочник. Под редакцией В.В.Стасенкова, И.С.Гутмана. М.: Недра, 1989. -270с.

36. Применение математических методов при поисках и разведке месторождений твердых полезных ископаемых. Секретариат СЭВ, М.: 1987г.

37. Рыбников A.B., Саркисов Г.Г. Стохастические геологические модели методы, технологии, возможности. / Нефтяное хозяйство. -2001. - №6. - с.22-25.

38. Регламент по созданию постоянно действующих геолого-технологических моделей нефтяных и газонефтяных месторождений. РД 153-39.0-047-00. М.: Минтопэнерго, 2000. - 150 с.

39. Семин Е.И. Геологическая неоднородность нефтеносных пластов и некоторые способы ее изучения // Тр. ВНИИ. М.: Гостоптехиздат, 1962. - Вып. 34. - С. 3^43.

40. Справочное руководство по проектированию разработки и эксплуатации нефтяных месторождений. Добыча нефти. Под общ. ред. Ш. К. Гиматудинова/ Р. С. Андриасов, И. Т. Мищенко,А- И. Петров и др. М., Недра, 1983, 455 с.

41. Справочное руководство по проектированию разработки и эксплуатации нефтяных месторождений. Проектирование разработки. Ш.К. Гиматудинов, Ю:Г. Борисов, М.Д. Розенберг и др. М., Недра, 1983, 463 с.

42. Спутник нефтегазопромыслового геолога. Справочник. Под редакцией Чоловского И.П. 1989. - М.: Недра. - 376с.

43. Стасенков В.В., Климушин ИМ., Бреев В.А. Методы изучения геологической неоднородности нефтяных пластов. М.: Недра, 1972.

44. Сыртланов В.Р., Майсюк Д.М., Трофимова И.А. Методика проведения экспертизы трехмерных моделей залежей нефти (на примере пласта БВ8° Самотлорского месторождения). Нефтепромысловое дело. 2003. - №10. - с. 14-17.

45. Токарев М.А. Комплексный геолого-промысловый контроль за текущей нефтеотдачей при вытеснении нефти водой. М.: Недра, 1990. -200с

46. Токарев М.А. Оценка и использование характеристик геологической неоднородности продуктивного пласта. 1983. - Уфа: изд-во УНИ.

47. Ф.Г.Гурари, Строение и условия образования клиноформ Западно-Сибирской плиты (история становления представлений): Монография. Новосибирск, СНИИГГиМС, 2003.

48. Чернова О.С. Седиментология резервуара. Учебное пособие по короткому курсу. 2004. - Томск. - 453с.

49. Чоловский И.П. Методы геолого-промыслового анализа при разработке крупных нефтяных месторождений. М.: Недра, 1966.

50. Carle S.F. and G.E. Fogg. 1996. Transition Probability-Based Geostatistics. Mathematical Geology, Vol.28, No.4, May 1996, p.453-476.

51. Journel A.G. and Gomez Hernandez, J.J., 1989, Stochastic Imaging of the Wilmington Clastic Sequence, SPE 19857.