К оглавлению журнала

 

УДК 550.834

© В. В. Корягин, 1990

Интерпретация аномалий сейсмической записи с помощью итеративного моделирования

В.В. КОРЯГИН (ВО ИГиРГИ)

При различном насыщении пластов-коллекторов водой, нефтью или газом происходит изменение их упругих свойств. Лабораторные исследования на керне и скважинные наблюдения показывают, что при пористости коллектора в пределах 15– 25 % и глубине его залегания до 2 км замещение воды нефтью вызывает уменьшение скорости продольных волн на 8–12, замещение газом – на 10–20 %. Акустическая жесткость уменьшается соответственно на 10–15 и 15–25 % [1].

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

Однако перечисленные задачи пока удается решать в условиях, когда пласт-коллектор имеет мощность не менее 10–15 м, является сравнительно однородным, а минимальная высота залежи около 10 м. В качестве примера успешного решения задачи по оконтуриванию нефтяной залежи, отвечающей указанным условиям, рассмотрим Суравинскую структуру в Куйбышевской области. Основная залежь нефти на ней вскрыта в бобриковском горизонте, в пласте Б2, общая мощность пласта-коллектора 21 м, нефтенасыщенной части в своде структуры 14. На временных разрезах по всем сейсмическим профилям, проходящим через свод структуры (рис. 1, а), с различной степенью отчетливости наблюдается аномалия записи в виде ослабления амплитуд по отражающему горизонту У, сопоставляемому с бобриковским. Пример такой аномалии представлен на фрагменте временного разреза по профилю 078363 на рис. 1, б. Моделирование убедительно показало, что ослабление амплитуд отражения У обусловлено нефтяной залежью в пласте Б2 (рис. 1, в, г).

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

В модели на рис. 1 потребовалось уменьшить скорость продольных волн в нефтенасыщенной части пласта Б2 до 3,8 км/с, т. е. на 17 % по отношению к скорости их в водонасыщенной части (4,6 км/с), чтобы амплитудные аномалии на синтетических разрезах соответствовали наблюденным (см. рис. 1, б, в). Высокий перепад скоростей при переходе в нефтенасыщенную часть пласта, не согласующийся с результатами некоторых лабораторных и скважинных сейсмометрических исследований [1, 4], можно объяснить двумя причинами. Во-первых, залежь нефти в пласте Б2 на Суравинской структуре характеризуется сравнительно высоким газовым фактором (около 100 м3/т), а добавление в жидкость хотя бы малого количества газа ведет к быстрому уменьшению скорости продольных волн [1]. Во-вторых, при формировании залежи Б2, по-видимому, особенно интенсивно происходило изменение петрофизических характеристик коллектора за счет растворения минералов на древних ВНК в условиях периодически дискретного поступления УВ в ловушку. Процесс такого изменения коллектора в объеме залежи подробно рассмотрен ранее [3].

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

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

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

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

Однако нередко характер задачи, решаемой с помощью моделирования, таков, что для интерпретации используется реальный временной разрез без миграции. В этом случае синтетический разрез должен вычисляться по нормальным лучам (в лучевом приближении) или по волновой теории [2]. Очевидно, процесс моделирования в данном случае становится более трудоемким.

Не менее важным вопросом является выбор импульса, моделирующего сейсмический сигнал. В практике моделирования широкое распространение получил импульс Пузырева, в котором начальную фазу j, преобладающую частоту w0 и затухание Р предлагается находить следующим образом: j принимается равной p/2 (симметричный импульс) на основании того, что в процессе обработки реальных сейсмических записей в результате применения всех видов фильтраций (деконволюция, полосовая фильтрация) стремятся получить элементарный сигнал симметричной формы (нульфазовый); w0 и Р определяются из теоретического амплитудного спектра импульса Пузырева, которым наилучшим образом (способом наименьших квадратов) аппроксимируется амплитудный спектр сигнала, наблюденного по осредненным ФАК реальных сейсмических записей. В дальнейшем наряду с коррекцией параметров модели найденные параметры импульса также уточняются.

Булатовская площадь. Данная площадь расположена в пределах Мелекесской впадины, в Куйбышевской области. По отложениям нижнего карбона она приурочена к осевой части Усть-Черемшанского прогиба ККСП, где мощность терригенных отложений нижнего карбона составляет около 400 м. Булатовское поднятие выявлено в 1983 г., детализировано в 1984 г. По отражающему горизонту У, сопоставляемому с кровлей тульского горизонта, поднятие характеризуется почти изометрической формой, очень контрастное, имеет амплитуду около 60 м, но очень малые размеры в плане – 1.5X1,5 км (рис. 2, а). На поднятии пробурено четыре скважины, открыто несколько залежей нефти, одна из них в пласте Б2 бобриковского горизонта. Залежь пластового типа, мощность нефтенасыщенной части в своде структуры равна 29 м, нефть характеризуется вязкостью и большой плотностью (0,942 г/см3).

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

Первым побуждением, естественно, было связать указанные аномалии с нефтяной залежью в пласте Б2. Однако сейсмическое моделирование позволило дать другую интерпретацию этим аномалиям. Как указывалось выше, нефть в пласте Б2 на Булатовской структуре характеризуется высокой вязкостью, большой плотностью и отсутствием растворенного газа, в связи с чем в модели на рис. 2, д скорость в нефтенасыщенной части пласта Б2 была незначительна, всего на 4 %, т. е. на 180 м/с, уменьшена по сравнению с водонасыщенной его частью. Такой слабый контраст в упругих свойствах при переходе из водонасыщенной части пласта (скорость – 4,52 км/с) к нефтенасыщенной (скорость 4,34 км/с) предопределил отсутствие заметных изменений в динамике записи на синтетических временных разрезах, полученных для модели с залежью (рис. 2, в) и без нее (рис. 2, г).

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

Модель 1. Пласты глин и алевролитов имеют линзовидный характер, распространены только в пределах поднятия, а в направлении к прогибам происходит их выклинивание, как это показано на рис. 2, д; начальное распределение сейсмических скоростей в данной модели было принято следующим: в песчанике пласта Б2 – 4,5 км/с, в чистых глинах – 3,5, в алевролитах – 3,8.

Модель 2. Пласты глин и алевролитов были заданы достаточно протяженными, т. е. не выклинивающимися на крыльях структуры, но в них происходит довольно резкое литологическое замещение в направлении к прогибам, где песчанистые фракции начинают преобладать и скорости продольных волн в последних (4,45 км/с) и во вмещающих песчаниках (4,52 км/с) мало отличаются.

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

Таким образом, представленные на рис. 2 результаты итеративного моделирования подтверждают предположение, что пласты глин и алевролитов в средней части пласта Б2 имеют линзовидный характер распространения и что именно с таким геологическим строением Булатовского поднятия связано появление аномалий в записи отражающего горизонта У. Залежь нефти не создает заметных аномалий эффектов. В пользу последнего вывода свидетельствует также тот факт, что уменьшение скорости в нефтенасыщенной части пласта Б2 до 4000 м/с, т. е. на 11 % ниже, чем в водонасыщенной части, вызывает резкое ухудшение сходства синтетического разреза с реальным (напомним, что в модели на рис. 2, д скорость в залежи уменьшена на 4 %).

Крюковская прощадь. Данная площадь расположена в юго-западной части Бузулукской впадины. Структура (Крюковская) выявлена в 1979 г., подготовлена по горизонтам карбона и девона в 1981 г. (рис. 3, а). Еще до бурения, т. е. по данным сейсморазведки, был ясен генезис Крюковского поднятия – типично эрозионного типа, связанного с останцами рельефа поверхности фундамента. В структуре осадочного чехла амплитуда поднятия последовательно уменьшается вверх по разрезу, о чем свидетельствуют такие данные: по поверхности фундамента амплитуда поднятия около 100, по кровле терригенного девона 60, бобриковского горизонта 35 м.

На поднятии пробурены две скважины, в своде и на северо-западном крыле, открыто две залежи нефти, в терригенном девоне (пласт Д1) и в бобриковском горизонте (пласт Б2).

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

На рис. 3, б приведен фрагмент временного разреза по профилю 048128, на котором в центре, на времени 0,92 с, благодаря достаточно корректной обработке с сохранением амплитуд аномалия проявилась особенно выразительно. Эта аномалия имеет классический вид, т. е. амплитуда первой фазы отражения У в области свода заметно ослаблена (темное пятно), а вторая фаза, наоборот, испытывает усиление более чем в два раза (яркое пятно). На времени 0,97 с прослеживается еще одно отражение, связанное с двумя глинистыми пластами в карбонатных отложениях турне. В записи этого отражения также наблюдается ослабление амплитуд в сводовой части поднятия.

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

Прежде всего, значительное увеличение амплитуд второй фазы отражения У, т. е. яркое пятно, объясняется тем, что в нижней части бобриковского горизонта на крыльях складки в направлении к своду происходит замещение пласта песчаника, имеющего мощность 10 м и скорость продольных волн 4,5 км/с, глинами со скоростью волн 3,5 км/с (рис. 3, д). В результате этого в сводовой части структуры с карбонатами турне контактируют глины с низкой скоростью продольных волн. Коэффициент отражения здесь почти в два раза выше, чем в песчаниках, залегающих на карбонатах турне в крыльевых частях и прогибах. В том, что залежь нефти не оказывает влияния на аномалию увеличенных амплитуд, можно убедиться при сравнении синтетических разрезов с залежью и без нее (рис. 3, в, г). На обоих разрезах данная аномалия имеет одинаковый вид, хотя скорость в залежи была уменьшена на 15 % по сравнению с водонасыщенной частью.

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

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

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

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

СПИСОК ЛИТЕРАТУРЫ

  1. Авербух А.Г. Изучение состава и свойств горных пород при сейсморазведке.– М.: Недра.– 1982.
  2. Корягин В.В., Сахаров Ю.П. Математическое моделирование в сейсморазведке.– М.: Наука.– 1988.
  3. Сахибгареев Р.С., Журавлева А.В. Литогенетические особенности коллекторов нефтяных залежей в связи с прямыми поисками // Советская геология.– 1987.– № 6.– С. 31–41.
  4. Трапезникова Н.А., Шушакова Н.С., Ворожцов Л.Н. Результаты математического моделирования сейсмических полей для нефтегазовых залежей Среднего Приобья // Прикладная геофизика.– 1979.– Вып. 94.– С. 45–56.

Abstract

It has been illustrated by the examples of the interpretation of anomalies of seismic record amplidudes that at comparatively simple geologic framework of a field these anomalies caused, generally by oil pool, may be used to solve a number of geological-field problems. Freguently encountered are the anomalies of amplidudes which are marked in the crestal parts of uplifts but which are not related to oil pools and appear to nesult from lithologic substitutions in terrigenous rocks, increased fracturing of carbonates above basement projections, etc. It has been demonstrated that seismic modeling is an efficient tool for revealing such false anomalies.

Рис. 1. Сейсмическое моделирование на Суравинской структуре с целью оконтуривания нефтяной залежи Б2.

а – структурная карта по отражающему горизонту У (контур залежи: 1 – по структурной основе, 2 – по динамике – амплитудам отражения У); б – реальный временной разрез; в – синтетический разрез; г – окончательная модель: 3 – известняки, 4 – аргиллиты, 5 – известняки глинистые, 6 – песчаники

Рис. 2. Результаты итеративного моделирования на Булатовском поднятии.

а – структурная карта по отражающему горизонту У: 1 – сейсмические профили, 2 – изогипсы отражающего горизонта У, 3 – приподнятый участок, 4 –глубокие скважины; 6 – реальный временной разрез по профилю 028434; синтетический разрез для модели: в – с залежью нефти, г – без залежи; д – окончательная модель (штриховкой показана залежь нефти в пласте Б2). Усл. обозн. см. на рис. 1

Рис. 3. Подбор модели на Крюковском поднятии.

а – структурная карта по отражающему горизонту У: 1 – сейсмические профили, 2 – изогипсы отражающего горизонта У, м; 3 – глубокие скважины; б – реальный временной разрез по профилю 048128, синтетический разрез для модели в – без залежи, г – с залежью; д – окончательная модель (штриховкой показана залежь нефти в пласте Б2). Усл. обозн. см на рис. 1