Развитие метода AVO (метод AVO – Amplitude Variation with Offset) началось с конца 1980-х гг. и связано с началом регистрации на временных сейсмических разрезах «ярких пятен» («bright spot»). Было замечено, что резкое изменение амплитуд отраженных волн на сейсмических разрезах может быть связано с газовыми залежами. Дальнейшие наблюдения показали, что яркие амплитудные аномалии могут быть вызваны также литологическими замещениями, и стал развиваться новый подход к изучению «ярких пятен» – анализ изменения амплитуд с удалением от источника возбуждения упругих волн (analysis of Amplitude Variation with Offset) [1, 2]. В основных алгоритмах AVO-анализа используется аппроксимация уравнения Шуэ [1, 2] вида: R(θ) = R0 + Grsin2θ, в которой R0 определяется только контрастностью акустической жесткости (AI = Vp*ρ) для продольных волн на границе двух сред, градиент GR – величина изменения амплитуд с удалением, определяется соотношением продольных (Vp) и поперечных (Vs) волн и является параметром, наиболее чувствительным к насыщению флюидом [1, 2].
В настоящее время существует много разнообразных направлений развития и применения данного метода, в которых величины R0 и GR, являясь основными атрибутами анализа, образуют вторичные атрибуты.
В рамках данной статьи рассматривается технология построения трехмерных кроссплотов AVO – кроссплотов R0 и GR с использованием третьего параметра в цветокодированном виде. Данная технология используется для решения задач прогнозирования газоносных объектов на месторождениях Вилюйской синеклизы.
В пределах рассматриваемых участков эффект аномального усиления амплитуд (эффект «ярких пятен») на реальных сейсмических данных дают как газовые залежи, так и угольные пропластки. Анализ каротажных кривых (АК и ГГКп) показал, что для изучаемой территории упругие характеристики углей и газонасыщенных песков имеют резкий контраст и их не сложно разделить по параметру AI (Vp*ρ), но вмещающие породы, представленные высокоскоростными глинами и аргиллитами, имеют более широкий диапазон изменений значений упругих характеристик, внося тем самым неоднозначность в интерпретацию аномалий.
Проблему неопределенности удалось решить путем построения теоретических петрофизических трендов, основанных на данных каротажа. В качестве анализируемых параметров использовались AI и Vp/Vs. Деление типов литологии по тому или другому параметру позволило оценить информативность последнего для интерпретации эффектов «ярких пятен» в сейсмическом волновом поле.
Информативный атрибут включался в цветокодированном виде в качестве третьего параметра для построения кроссплотов AVO, что позволило получить многомерную информацию и интерпретировать эффекты «яркого пятна» в пределах изучаемой территории.
Краткая характеристика участка работ
В административном отношении площадь работ находится на территории Вилюйского и Верхневилюйского районов Республики Саха (рис. 1).
Рис. 1. Республика Саха (Якутия) (Научно-технический сборник № 3(19) /2014) [4]
В нефтегазоносном отношении участок исследований расположен в северо-западной части Вилюйской газоносной области, входящей в состав Хатанго-Вилюйской нефтегазоносной провинции. В разрезе данного участка ранее проведенными работами были выявлены два продуктивных комплекса: верхнепермский и нижнетриасовый. На соседних площадях были получены притоки газа из верхнеюрских и нижнемеловых отложений [4].
В тектоническом отношении площадь расположена на северо-западном склоне Вилюйской гемисинеклизы. Бурением вскрыты породы палеозойского (среднего – верхнего) и мезозойского возраста.
Фундамент в районе исследований и к востоку от него образует Вилюйскую депрессионную зону, разделяющую Алданское, Анабарское и Оленекское поднятия. В пределах исследуемого участка рельеф фундамента представляет собой выступ, осложненный серией складок северо-восточного простирания. Доминирующая северо-восточная ориентировка тектонических элементов сохраняется на всех вышезалегающих этажах, прерываясь лишь в центральной части площади сдвиговыми деформациями, формирующими в центральной части площади обособленные блоки. Вся сложившаяся система трещин и деформаций является благоприятным фактором для миграции УВ (рис. 2, а, б, в, г).
Рис. 2. Тектоническое строение территории: а, б – схема геодинамической эволюции сдвига; в, г – внутреннее строение разломной зоны в интервале триаса и нижнего мела; д, е – сейсмоимидж аномалий нижнего мела с интерпретацией и без интерпретации
Анализ сейсмических данных
На площади исследований проведены сейсморазведочные работы 3D. На сейсмических разрезах в пределах исследуемого участка в интервалах продуктивных горизонтов и в тех интервалах, где по материалам ранее проведенных сейсморазведочных работ 2D перспективные объекты не выделялись, фиксируются амплитудные аномалии типа «яркое пятно».
Сейсмоимидж аномалий верхней части разреза, в интервале отложений бергеинской свиты, соответствует образу тектонически экранированных ловушек (рис. 2, д, е). Подобный вывод сформировался по совокупности внешних признаков, таких как четкие очертания динамических аномалий, вписывающихся в блоки системы в купольной части поднятия, при наличии проводящей системы трещин и сдвиговых деформаций.
Сейсмоимидж объекта в низах разреза, в интервале отложений хомустанской свиты, по протяженной форме, окаймляющей поднятие с юго-запада и четким границам, повторяющим контур склона, соответствуют образу стратиграфической ловушки УВ.
Выделение газонасыщенных аномалий по сейсмическим данным традиционно проводится по технологии AVO, в частности, путем построения кроссплотов GR как функция R0 («интерцепт-градиент»). В данном случае в анализ привлекалась третья компонента, которая в цветокодированном виде содержала информацию о литологии: код раскраски точек кроссплота соответствовал тому или иному типу литологии.
Анализ каротажных диаграмм
Интерпретацию литологического состава и флюидонасыщения эффективно проводить по кроссплотам каротажных диаграмм [3, 5]. В качестве анализируемых параметров были выбраны упругие характеристики AI и Vp/Vs (Vp/Vs как функция AI). Код раскраски петрофизических трендов основан на выделении четырех литотипов в исследуемом интервале разреза: «уголь», «газоносный песчаник», «глина», «вмещающие породы». Данная классификация основана на замерах в разрезе скважины потенциалов собственной поляризации, естественной радиоактивности пород и водородосодержания. Кроссплот включает в себя тренды распределения акустических свойств для углей и осадочных пород в зависимости от содержания в них глинистой компоненты и насыщения (рис. 3) [6].
Рис. 3. Тренды распределения акустических свойств для пород: а) бергеинской свиты; б) хомустанской свиты; в) идеализированные теоретические тренды (Одегаард и Авсет, 2004)
На построенных по сейсмическим данным кроссплотах AVO локализуются отдельные группы аномалий, что говорит о чувствительности исследуемого волнового поля к контрастности скоростей поперечных волн и соотношений плотности на границах двух контактирующих сред. Классификацию AVO аномалий по литотипам можно осуществить путем введения в AVO кроссплоты в цветокодированном виде третьей компоненты.
Технологическая схема использования анализа AVO сводится к следующей последовательности:
1. Типизация разреза по диаграммам каротажа с назначением цветовых кодов.
2. Построение теоретических петрофизических трендов, выбор информативного параметра.
3. Построение трехмерных кроссплотов AVO, классификация аномалий AVO по типу литологии.
4. Выделение классов аномалий на временных сейсмических разрезах и на картах изохрон.
Построение петрофизических трендов и интерпретация результатов AVO-анализа в интервале нижнемеловых отложений
Исследованиями кернового материала в пределах этой и соседних площадей в отложениях нижнего мела, в бергеинской свите, установлены мощные песчаные пласты с хорошими коллекторскими свойствами. Покрышками песчаных тел являются слои высокоскоростных глин, залегающих под угленосной толщей. Поскольку на данном месторождении газонасыщенность нижних этажей разреза установлена по данным бурения, то наличие проводящей системы трещин и разломов, с одной стороны, и наличие покрышек, сложенных высокоскоростными глинистыми пропластками, с другой, могут служить благоприятными факторами для скоплений УВ в ловушках тектонического типа в бергеинской свите.
Анализ каротажных кривых и типизация разреза проводились по скважине, хорошо охарактеризованной керновым материалом. Код раскраски, присвоенный литотипам на диаграммах каротажа, перенесен в область кроссплота для интерпретации петрофизических трендов.
Из характера распределения данных на графике следует, что литотип «насыщенный песчаник» и литотип «уголь» контрастны по своим акустическим параметрам. Это означает, что выбор параметра AI для включения в кросс-плот AVO в качестве третьей переменной позволит с достаточной долей уверенности определить на кроссплотах AVO и в сейсмическом волновом поле аномалии, вызванные угольными пропластками.
Распределение данных каротажа по такому параметру, как «контрастность скоростей», показывает, что тренды анализируемых литотипов локализованы в разных диапазонах значений Vp/Vs. Из этого следует, что на кроссплотах AVO аномалии, вызванные насыщенными песчаниками, и аномалии, вызванные угольными пропластками, также локализуются в отдельные группы аномалий относительно фоновой линии, что позволит их идентифицировать в сейсмическом волновом поле (рис. 4).
Рис. 4. Интерпретация результатов AVO-анализа для интервала отложений бергеинской свиты: слева – по коду раскраски кроссплота в сейсмическом волновом поле выделен литотип «газоносный песчаник»; справа – то же для литотипа «уголь»
Построение петрофизических трендов и интерпретация результатов AVO-анализа в интервале отложений верхней перми
На площади исследования бурением вскрыта толща терригенных отложений верхней перми, где среди мощных плотных пластов аргиллитов, глин и углей залегают пласты песчаников различной глинистости. В скважинах 230 и 231 данный интервал представлен газоносными песчаниками с низкой пористостью и слабой проницаемостью, ввиду чего к испытанию не рекомендовался. В скв. 238, расположенной в купольной части структуры, из этого же пласта при испытании был получен промышленный приток газа.
Эти газоносные залежи производят во временном сейсмическом поле эффект «ярких пятен». Ниже пласта на 30 м залегает мощный угольный пропласток, который дает аналогичный эффект «яркого пятна» (рис. 5). Контрастность литотипов «уголь» – «газоносный песчаник» по параметру AI можно установить по характеру распределения данных на диаграммах петрофизических трендов для данного интервала. Значения по этому параметру не перекрываются, что является благоприятным фактором для включения AI в кроссплот AVO в качестве третьей переменной.
Рис. 5. Интерпретация результатов AVO-анализа для интервала отложений хомустанской свиты (верхняя пермь): а) сейсмоимидж ловушки; б) по коду раскраски кроссплота в сейсмическом волновом поле выделен литотип «газоносный песчаник»; в) то же для литотипа «уголь»
Перекрытие значений по параметру Vp/Vs говорит о том, что можно ожидать, что группы аномалий на кроссплотах AVO также будут перекрываться. Тем не менее в данной работе удалось идентифицировать в сейсмическом волновом поле амплитудные аномалии, связанные с газоносностью ловушки неструктурного типа в верхнепермских отложениях, и отличить их от аномалий, вызванных мощным пластом угля, залегающим ниже ловушки на 30 м (рис. 5).
Выводы
По результатам проведенных сейсморазведочных работ 3D выявлены тектонические ловушки в отложениях нижнего мела и уточнены контуры литологической ловушки верхней перми, детализировано тектоническое строение площади и установлена проводящая роль тектонических элементов. Построены петрофизические тренды и проведена интерпретация литологического состава для отложений бергеинской свиты нижнего мела и хомустанской свиты верхней перми. Проведен AVO-анализ с включением в технологическую цепочку результатов интерпретации петрофизических трендов, построенных по данным каротажных диаграмм, что позволило сократить риски при оценке перспектив газоносности исследуемых интервалов.
Данная работа проводилась на базе материалов, систематизированных и обработанных в ООО «ИНГЕОСЕРВИС» (г. Тюмень). Автор выражает свою благодарность сотрудникам и руководству компании за предоставленные материалы.