Полный текст

Введение При оценке достоверности построения геологической модели на основе трехмерного геологического моделирования наиболее распространены две технологии: перекрестной оценки и стохастического моделирования с оценкой неопределенности. Технология перекрестной оценки заключается в последовательном исключении скважин из набора, использовавшегося при построении модели, и измерении погрешности построения модели в точках скважин. Технология стохастического моделирования позволяет получить представительный ансамбль реализаций, который может учитывать неопределенность в структурных, литологических и петрофизических построениях. На основании этих данных определяются достоверность построения геологической модели, возможные диапазоны разброса параметров модели, гистограммы распределения запасов, зоны повышенного риска бурения и др. [1]. Остановимся на неопределенности геологических моделей, которая явилась причиной возникновения технологии стохастического моделирования с оценкой неопределенности. Высокая степень неопределенности существует в любой геологической модели месторождения. Пессимистические и оптимистические модели создаются для объективного понимания неопределенности и оценки возможных рисков. Имея множество возможных неопределенностей, зачастую невозможно узнать, как наилучшим образом создать эти модели и как сделать их наиболее информативными. С помощью стохастического моделирования петрофизических свойств появляется возможность создавать множество равновероятных реализаций, каждая из которых соответствует исходным данным и полевой информации [2]. Диапазон возможных результатов можно оценить в процессе моделирования. В итоге можно выбрать для прогнозирования различных свойств наиболее приемлемые модели. В данной работе представлено 4 различных варианта построения трехмерных геологических моделей. При оценке достоверности построения геологических моделей методом стохастического моделирования по каждому из вариантов был выполнен процесс многократного (32 реализации) перестроения модели на основе пересчета куба пористости, литологии и привлечения дополнительных исходных данных при задании настроек стохастического алгоритма моделирования. Все типы исходных данных являются дополнительными и используются для ограничения стохастического алгоритма, т.е. для получения более реалистичного результата построения трехмерных литолого-фациальных моделей. Рассмотрим наиболее подробно каждый из вариантов моделирования. Первый вариант моделирования В первом варианте путем стохастического моделирования выполнялось распределение коэффициента пористости Кп. Преимущества данного инструмента в том, что можно получать распределение на основе исходных данных и знаний об их трендах, а также оценивать неопределенность построений и связанный с ней риск. В основе стохастического моделирования лежит геостатистический метод – кригинг, который основывается на двух принципах: использование пространственных соотношений между неизвестными и известными данными, чтобы рассчитать неизвестные; и учет близко расположенных данных для корректного взвешивания каждой точки [3, 4]. Исходными данными для петрофизического моделирования являются: – трехмерная геологическая сетка; – скважинные кривые, которые осреднены до разрешения ячеек трёхмерной сетки; – результаты интерпретации ГИС (непрерывная кривая Кп), а также данные Рис. 1. Распределение Кп по 32 реализациям первого варианта моделирования о предельных значениях параметров пород продуктивных пластов [5]. В итоге по 32 реализациям первого варианта моделирования было получено распределение коэффициента пористости (рис. 1). В данном случае гистограмма коэффициента пористости, построенная по 32 реализациям первого варианта моделирования, имеет бимодальный вид распределения, причем выделенные моды по отношению друг к другу являются практически равнозначными. Следовательно, выдвигаем гипотезу, что данный вид распределения обусловлен присутствием в выборке двух типов фаций: протоки дельты и дельтового залива. Исходя из вышеизложенного необходимо отметить, что в пределах нижнетиманского циклита Кирилловского месторождения целесообразно выполнить работу по разделению рассматриваемого фациального комплекса на два типа фаций [6]. Таким образом, при расчете второго варианта моделирования для построения трехмерных кубов пористости в качестве дополнительных исходных данных были учтены результаты литолого-фациального анализа девонской терригенной толщи Кирилловского месторождения, где рассматриваемый в данной работе циклит представлен двумя основными фациальными комплексами типа протока дельты и дельтовый залив. Второй вариант моделирования Во втором варианте моделирования, помимо непрерывной кривой пористости, были учтены результаты литолого-фациального анализа (по А.П. Вилесову) и данные сейсмофациального анализа (по И.С. Путилову) продуктивной живетско-нижнефранской терригенной толщи Кирилловского месторождения нефти. По данным литолого-фациального анализа нижнетиманский циклит представлен двумя основными фациальными комплексами: дельтовые заливы и дельтовые протоки [7], которые были использованы при объектном моделировании песчаных тел. Из результатов сейсмо-фациального анализа были использованы границы зон распространения рассматриваемых фациальных комплексов. При обобщении непрерывных данных геофизических исследований скважин, а также результатов литолого- и сейсмофациального анализа была получена и преобразована в дискретный вид фациальная кривая Zones, в которой каждому фациальному типу соответствует определенный индекс, например, для дельтовых заливов был задан индекс 0, а для дельтовых проток – 1. Полученная кривая фаций необходима в качестве одного из исходных данных для выполнения фациального моделирования песчаных объектов, выделенных в пределах нижнетиманского циклита по результатам литолого-фациального анализа керна [5]. Моделирование песчаных тел отдельно для каждого типа фаций осуществлялось с помощью инструмента SedSeis. Этот метод наиболее подходит для моделирования рассматриваемых фациальных комплексов, когда форма объектов напоминает русло [8]. В качестве исходной информации, помимо полученной ранее дискретной фациальной кривой, были использованы границы зон распространения рассматриваемых типов фаций, которые были выделены по данным сейсморазведки методом 3D, а также осевые (цен- Рис. 2. Распределение Кп по 32 реализациям второго варианта моделирования отдельно для каждого типа фаций тральные) линии песчаных объектов. Сопоставление получившихся распределений фаций велось визуально по горизонтальным и вертикальным срезам, по всему объему трехмерных геологических моделей. Затем так же, как и в первом варианте, путем стохастического моделирования выполнялся расчет куба пористости с учетом непрерывной скважинной кривой Кп, данных литолого- и сейсмофациального анализа, а также результатов фациального моделирования песчаных объектов. В итоге по 32 реализациям второго варианта моделирования было получено распределение коэффициента пористости (рис. 2). Представленные на рис. 2 гистограммы коэффициента пористости имеют одномодальный вид распределения, который эмпирически установлен для большинства однородных типов пород, принадлежащих определенному типу фаций: дельтовая протока и залив дельты. Также на гистограммах представлено предельное значение коэффициента пористости (Кп* = 10 %). Анализируя представленные результаты, можно сделать вывод, что распределение частоты встречаемости коллекторов в песчаниках и алевролитах дельтовых проток имеет ярко выраженное смещение в сторону высокоемких пород, тогда как среди песчано-алевритовых пород дельтовых заливов в основном преобладают неколлекторы (значение Кп < 0,1 д. ед.), которые составляют 95 % от всей выборки [9]. Также необходимо отметить, что полученные результаты по второму варианту моделирования указывают на корректное разделение рассматриваемого в данной работе нижнетиманского циклита на отдельные фациальные комплексы типа протока дельты и дельтовый залив. Третий вариант моделирования В третьем варианте при построении трехмерных кубов пористости, помимо использованных в предыдущих вариантах дополнительных данных, были учтены отдельно для каждого типа фаций результаты статистической характеристики непрерывной скважинной кривой Кп по ряду месторождений нефти (Кирилловское, Кустовское, Мало-Усинское и Андреевское). Данные результаты представлены в виде комплексных гистограмм нормального распределения Кп и были использованы в качестве трендов при настройках стохастического алгоритма, необходимого для распределения коэффициента пористости, отдельно для каждого типа фаций. По результатам анализа комплексных гистограмм нормального распределения Кп можно сделать следующий вывод: трендовые гистограммы имеют одномодальное, близкое к нормальному виду распределение параметра пористости; наилучшими коллекторскими свойствами обладают песчаники и алевролиты, приуроченные к фациальному комплексу типа протока дельты; наблюдается практически полное отсутствие коллекторов в породах дельтовых заливов [5]. Затем, что так же, как и в предыдущих вариантах, путем задания определенных настроек стохастического алгоритма выполнялось распределение параметра пористости с учетом непрерывной скважинной кривой Кп, данных литолого- и сейсмофациального анализа, результатов фациального моделирования песчаных объектов, а также с учетом трендовых гистограмм непрерывной скважинной кривой Кп по ряду месторождений нефти (Кирилловское, Кустовское, Мало-Усинское и Андреевское) [10]. Анализ полученных результатов по третьему варианту моделирования показал, что значения коэффициента пористости так же, как и в предыдущем варианте, в фациальных комплексах типа протока дельты и залив имеют одномодальный вид распределения и подтверждают корректное разделение рассматриваемого в данной работе нижнетиманского циклита на основные фациальные комплексы. Четвертый вариант моделирования В четвертом варианте при построении трехмерных параметров пористости, помимо использованных в предыдущих вариантах исходных данных, дополнительно были учтены в качестве трендов прогнозные схемы коэффициентов пористости, которые получены при детальном изучении геологического строения Кирилловского месторождения нефти с помощью пространственной сейсморазведки методом 3D. В ячейках стратиграфической сетки были рассчитаны значения прогнозного коэффициента пористости. По данной сетке была построена промежуточная карта Кп для интервала, соответствующего нижнетиманским отложениям. Учитывая слабые стороны статистической модели, полученные данные были обработаны геостатистическими методами для повышения их достоверности. С использованием геостатистического алгоритма кокрайгинга была получена окончательная схема прогнозных Кп в продуктивных интервалах терригенного девона [11, 12]. В итоге для нижнетиманских терригенных отложений были построены карты прогнозного коэффициента пористости исходя из анализа количественных связей между интервальными атрибутами волнового поля и геолого-геофизическими данными. Эти прогнозные карты пористости были использованы в качестве трендов при дальнейшем моделировании [13]. Полученные в результате прогнозные схемы пористости в интервале нижнетиманских отложений в пределах скважин соответствуют скважинным параметрам (по РИГИС), а в межскважинном пространстве отражают неоднородность свойств. Результаты по четвертому варианту моделирования, как и в предыдущих вариантах, подтверждают корректное разделение рассматриваемого в данной работе нижнетиманского циклита на основные фациальные комплексы: протока дельты и дельтовый залив. Заключение В итоге по каждому из рассмотренных в данной работе вариантов моделирования было получено по 32 реализации распределения трехмерного параметра пористости по разрезу терригенных отложений нижнетиманского циклита (т.е. Кп распределен как в породах-коллекторах, так и в неколлекторах). Анализируя полученные результаты трехмерных параметров пористости, представленные по каждому варианту моделирования, с помощью математической статистики отдельно по каждому фациальному комплексу, рассчитали t-критерий и уровень значимости р. Данные значения получены путем последовательного сравнения друг с другом вариантов моделирования. Результаты сравнения вариантов отражены в таблице. По результатам, представленным в таблице, можно сделать вывод, что для фациальных комплексов типа протока дельты и дельтовый залив при каждом последующем варианте моделирования наблюдается увеличение значения t-критерия и снижение уровня значимости р, т.е. происходит уменьшение дисперсии и значения коэффициента пористости становятся ближе к среднему значению. Данная тенденция связана с привлечением к каждому после- Результаты последовательного сравнения вариантов моделирования Вариант Протока + залив Протока Залив t-критерий р t-критерий р t-критерий р 1_2 0,8697130 0,384457 – – – – 2_3 1,2826250 0,199624 0,266172 0,790107 0,113961 0,909268 3_4 1,0217000 0,306923 0,881625 0,377980 2,482120 0,01306 дующему варианту моделирования дополнительных исходных данных. Исходя из вышеизложенного, можно сделать вывод, что результаты, полученные по четвертому варианту моделирования, являются наиболее информативными и достоверными, а также подтверждают корректное разделение рассматриваемых в данной работе нижнетиманских терригенных отложений Кирилловского месторождения нефти на основные типы фаций (протока дельты и дельтовый залив). Затем по всем представленным в четвертом варианте моделирования реализациям проводился расчет трехмерного параметра литологии путем разделения полученного параметра пористости по его предельному значению (Кп*) [14, 15]. Данный трехмерный параметр наиболее точно характеризует пространственное положение коллекторов в пределах рассматриваемой площади.

Об авторах

Александр Борисович Дерюшев

Автор, ответственный за переписку.
Email: ABDeryushev@mail.ru
614066, г. Пермь, ул. Советской Армии, д. 29

инженер 1-й категории отдела геологического моделирования, филиал ООО «ЛУКОЙЛ-Инжиниринг» «ПермНИПИнефть» в г. Перми

Денис Владимирович Потехин

Email: Potekhin@permnipineft.ru
614066, г. Пермь, ул. Советской Армии, д. 29

начальник отдела геологического моделирования, филиал ООО «ЛУКОЙЛ-Инжиниринг» «ПермНИПИнефть» в г. Перми

Список литературы

Статистика

Просмотры

Аннотация - 273

PDF (Russian) - 49

Ссылки

• Ссылки не определены.