Математическое моделирование в сейсморазведке

Лекция 1

ВВЕДЕНИЕ

МЕТОД МАТЕМАТИчЕСКОГО МОДЕЛИРОВАНИя ПОЛУчИЛ РАСПРОСТРАНЕНИЕ В
СЕЙСМОРАЗВЕДКЕ ПРИМЕРНО С СЕРЕДИНЫ 60-Х ГОДОВ. ИСПОЛЬЗОВАНИЕ СИНТЕТИчЕСКИХ
СЕЙСМОГРАММ (СС), яВЛяВШИХСя РЕЗУЛЬТАТОМ РЕШЕНИя ОДНОМЕРНОЙ ДИНАМИчЕСКОЙ
ЗАДАчИ, ИМЕЛО СЛЕДУЮЩИЙ ЦЕЛИ:
. анализ процесса формирования поля отраженных волн в тонкослоистой среде(
. оценку роли многократных волн в этом поле(
. определение сейсмических эффектов, обусловленных изменением литологии или углеводородосодержания и др.

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

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

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

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

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

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

1 Системный анализ проблемы интерпретации данных сейсмических наблюдений

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

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

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

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

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

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

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

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

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

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

VI уровень в рассматриваемой схеме представляют процессы принятий по коррекции параметров в общем случае всех операций уровней II–V. В частности, при наименее "глубокой" обратной связи корректируются параметры сейсмомоделирования, т. е. сейсмогеологическая модель и модель импульса падающей волны. Исходными данными для принятия таких решении являются оценки сходства ("рассогласования"), получаемые на уровне V.

2 Теоретические вопросы автоматизированной интерпретации данных сейсморазведки

Лекция 2
ТАБЛИЦА 1. ВЛИяНИЕ ПАРАМЕТРОВ ДВУМЕРНОГО СЕЙСМОМОДЕЛИРОВАНИя

НА ХАРАКТЕРИСТИКИ ОТРАЖЕНИЙ

Способы построения сейсмических моделей геологических сред

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

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

1 Построение одномерных моделей

Исходная информация, т. е. значения детальных скоростей и плотностей, для построения одномерных тонкослоистых моделей может быть получена несколькими способами:
1. По данным акустического (АК), гамма-гамма (ГГК) или гравитационного каротажей после соответствующей их обработки; обработка АК обычно включает процедуры вычисления скоростей с учетом кавернометрии, коррекции полученных скоростей по сейсмическому каротажу (СК), осреднения и др.;

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

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

Из перечисленных способов предпочтение следует отдать использованию данных АК и ГГК.

Осреднение данных АК и ГГК

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

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

Способ осреднения с порогом. Применение его позволяет получить тонкослоистую модель в виде серии однородных слоев большей мощности по сравнению с исходным разрезом. Все границы в модели представляются границами первого рода. Сущность алгоритма осреднения в данном способе заключается в том, что по заданным (V – величине значимой скоростной дифференциации и ((min – минимальной временной мощности слоев из разреза исключаются тонкие слои, время пробега в которых (ti < ((min, объединяются слои с номерами i и i-1, если разница скоростей в них удовлетворяет условию

|Vi – Vi–1| ( (V

Значение скорости в объединенном слое вычисляется как среднее из Vi и
Vi-1. Пороговое значение скачка скорости (V может быть различным для разных частей разреза.

Изменяя (V, можно менять число слоев в модели N, так как оно тем меньше, чем больше (V. Это может быть использовано для автоматического поиска моделей с числом слоев, находящихся в заданных пределах Nmin – Nmax.

2 Построение двумерных моделей

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

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

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

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

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

1 Построение модели по данным бурения

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

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

2 ПОСТРОЕНИЕ МОДЕЛЕЙ ПО ДАННЫМ БУРЕНИя И СЕЙСМОРАЗВЕДКИ

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

(см. п. 1) как доминирующие при формировании отраженной волны. Если по какой-либо скважине получается невязка, то она "разбрасывается" по линейному закону в глубины сейсмической границы между скважинами.
3. На полученный в результате такой коррекции сейсмический разрез, который можно назвать базисной толстослоистой моделью, в точках расположения скважин наносятся тонкослоистые модели, соответствующие моделируемому интервалу. В пределах моделируемого интервала проводятся границы отдельных литологически однородных тонких слоев. При этом в зависимости от предполагаемой степени сложности двумерной модели подходы к ее построению могут быть различными. В зонах выдержанной корреляции сейсмических данных, которые, как правило, соответствуют согласному или близкому к нему залеганию пород, эти границы проводятся так, чтобы они соединяли отметки по скважинам и были параллельны сейсмическим границам между скважинами. Участки изменений сейсмических данных (схождение осей синфазности, изменения формы и интенсивностей колебаний, разрывы в корреляции) тщательно анализируются и с учетом данных по скважинам задаются возможные модели изменений мощности слоев, литолого-фациальных замещений, появления углеводородов и др. Нередки случаи, когда в пределах одного моделируемого интервала встречаются участки различной сложности.
4. Задаются упругие параметры (скорости и плотности) во всех слоях модели, при этом в точках между скважинами эти параметры находятся путем линейной интерполяции значений, полученных ранее в процессе формирования одномерных моделей в точках расположения скважин.

3 Построение моделей по данным сейсморазведки

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

Уральской провинции – чаще из интерполированных или экстраполированных сейсмокаротажных данных.
2. Интервал временного разреза, соответствующий моделируемому объекту, преобразуется во временной разрез волновых сопротивлений по методике псевдоакустического каротажа (ПАК).
3. В ряде точек профиля строятся одномерные модели волновых сопротивлений.

Затем от волновых сопротивлений с использованием формулы ( =аVb, где ( – плотность, V – скорость, переходят к оценкам скорости и плотности.

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

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

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

4 Влияние нефтегазонасыщенности на упругие свойства пород

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

Прямые измерения в скважинах с помощью сейсмического просвечивания и
СК выполнены в ограниченном объеме и полученные результаты не всегда достаточно точны. Обобщение данных показывает, что в нефтенасыщенных песчаных коллекторах при глубинах 1500–3000 м и средней пористости 20% скорость продольных волн уменьшается на 6–12%, в газонасыщенных коллекторах
– на 15–30% по сравнению с водонасыщенным коллектором.

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

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

, (2.1)

где Uп и (п – параметры, зависящие соответственно от упругости и плотности флюида; Uск и (ск – параметры, характеризующие упругость и плотность скелета (остова) породы.

Значения U и ( следующим образом выражаются через свойства твердого материала породы и насыщающего ее флюида:
1) (ск = (тв (1 – Kп), где (тв – плотность материала, слагающего твердую фазу породы, Kп – пористость;
2) (п = (ф Kп, где (ф – плотность флюида, т. е. плотность воды, нефти, газа или их смеси;
3) , где (ск – сжимаемость скелета (относительное изменение объема скелета при всестороннем упругом сжатии породы), Gск – модуль сдвига скелета;
4) где (тв – сжимаемость материала, слагающего скелет породы, (ф – сжимаемость флюида, величины (тв и (ск связаны соотношением (ск = (тв + Kп(п ((п – сжимаемость порового пространства).

При использовании формулы (2.1) основная трудность заключается в выборе величин (ск и Gск.

Для приближенных расчетов можно использовать уравнение среднего времени (уравнение Уилли)

, (2.2)

где Vп – скорость в коллекторе, заполненном флюидом; Vск – скорость в скелете; Vф – скорость во флюиде, Kп – коэффициент пористости. Формула
(2.2) справедлива для хорошо сцементированных пород. Величину плотности можно оценить по уравнению

(п = (ск (1 – Kп) + (фKп, (2.3)

где (п – плотность коллектора, заполненного флюидом, (ск – плотность скелета, (ф – плотность флюида.

Если поры заполнены несколькими компонентами, например газ–вода, нефть–вода и т. д., то имеет место уравнение

(п = (ск (1 – Kп) + (фKп + ((в – (ф)SвKп,

где (в – плотность воды, Sв – коэффициент водонасыщенности.

Методика интерпретации на основе итеративного моделирования

1 Особенности получения и обработки сейсмических данных, интерпретируемых на основе моделирования

Главное требование, предъявляемое к данным сейсмических наблюдений, которые интерпретируются с помощью итеративного моделирования, состоит в повышенном отношении сигнал/помеха. Опыт сейсмомоделирования показывает, что нижний предел отношения энергии сигнала к энергии помехи, равный 10 –
15, является достаточным для того, чтобы в процессе итеративного подбора параметров модели достичь достаточно высокую степень сходства СВР и реального временного разреза (РВР). Это предельное значение установлено на основе тестового моделирования и сопоставления СВР и РВР по нормированной функции взаимной корреляции (НФВК) и значений отношения сигнал/помеха на
РВР по одинаковым фрагментам временных разрезов. На рис. 4 показан пример такого сопоставления по профилю 39 Северо-Маркинской площади, из которого видно, что сходство СВР и РВР до 0,8 и выше удавалось получить только на участках, где отношение сигнал/помеха на РВР достигало 10 – 15 и выше.

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

ДОСТИЖЕНИЕ ПОДОБНОГО КАчЕСТВА РВР ЕСТЕСТВЕННО НАКЛАДЫВАЕТ БОЛЕЕ
ЖЕСТКИЕ ТРЕБОВАНИя НА МЕТОДИКУ ПОЛЕВЫХ НАБЛЮДЕНИЙ И ПОСЛЕДУЮЩУЮ ОБРАБОТКУ
СЕЙСМИчЕСКИХ ДАННЫХ.

1 Методика полевых наблюдений

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

На поисковом этапе исследований, в целях выявления рифогенных образований, грабенообразных прогибов, выступов кристаллического фундамента методика полевых работ может быть близка к производственной или отличаться от нее некоторым увеличением мощности интерференционных систем при возбуждении и приеме. Основные элементы такой методики следующие:
1) плотность сети профилей 1,5–2,0 пог. км на 1 км2;
2) схема наблюдения – в основном центральная;
3) кратность перекрытия 12 или 24;
4) максимальное расстояние взрыв – прибор Хmax = 1700–2500 м;
5) вынос 25–200 м;
6) расстояние между каналами 40–50 м;
7) группирование сейсмоприемников до 36 на канал, причем расположение приемников в одну или две линии на базе не более 50 м;
8) возбуждение – взрывы в одиночных скважинах с оптимальной глубины или из группы мелких (4–5 м) скважин на базе не более 40–50 м.

При детальных исследованиях требования к методике полевых наблюдений повышаются и сводятся к следующему.
1) плотность профилей должна быть не менее 3 пог. км на 1 км2, причем при детализации, например, грабенообразных прогибов большую часть профилей следует ориентировать вкрест прогиба с расстоянием между ними не более
500 м;
2) в целях повышения пространственной разрешенности расстояние между каналами не должно превышать 25–30 м;
3) группирование сейсмоприемников увеличивается до 48–60 элементов на канал, причем эти элементы располагаются по площади в виде 4–5 параллельных ниток; база группы должна быть не более 50 м.

2 Методика цифровой обработки

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

Выполнение указанных требований достигается при использовании усложненного графа обработки, содержащего следующие процедуры:
1) демультиплексация цифровых записей (DMXT);
2) редакция (REDX);
3) коррекция амплитуд за геометрическое расхождение и поглощение (RAMP);
4) вычитание среднескоростных волн-помех (RECON);
5) минимально-фазовая деконволюция исходных записей (DECVTX);
6) широкополосная фильтрация исходных записей (FILVTX);
7) коррекция амплитуд за неидентичность условий возбуждения и приема

(NORM);
8) коррекция статических поправок (SUMLAK);
9) коррекция кинематических поправок (сканирование или вертикальные спектры, KINVC);
10) автоматическая коррекция статических поправок (PAKS);
11) накапливание по ОГТ (SUMLC);
12) погоризонтный анализ скоростей (горизонтальные спектры скоростей,

HORSP);
13) независимая потрассовая коррекция остаточных фазовых сдвигов в нескольких временных окнах (WINCOR);
14) когерентная фильтрация (AMCOD);
15) нуль-фазовая деконволюция по разрезу (ZEDEC);
16) широкополосная фильтрация по разрезу (FILVTX);
17) когерентная фильтрация (AMCOD);
18) миграция (MIGFK);
19) псевдоакустический каротаж (РАК).

2 Выбор способа решения прямой

динамической задачи

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

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

Для оценки величины области формирования отраженного импульса обычно используется параметр первой зоны Френеля F, который рассчитывается по известной формуле(

где Н – глубина залегания отражающей границы; ( – длина волны. Если протяженность отражающего элемента, связанного с какой-либо неоднородностью в геологическом разрезе, составляет величину F зоны Френеля и более, то этот элемент отобразится на временном разрезе с максимальной амплитудой, соответствующей отражению от бесконечно длинной границы. При уменьшении горизонтальных размеров элемента (меньше F) он будет отображаться на временном разрезе с заметным уменьшением амплитуды, все меньше походить на отражение и все больше приобретать вид дифракции, соответствующей отражающей точке.

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

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

1 Пример 1. Моделирование микрограбенов

Данный пример (рис. 5) иллюстрирует отличие волновых полей от грабенообразных прогибов при различной их ширине. Последняя варьировалась, исходя из величины зоны Френеля, которая для модели на рис. 5, а при видимой длине волны ( = 160 м и глубине границы Н = 2400 м составляет F =
880 м. Поэтому ширина грабенов была задана следующей: l1 = 0,5F = 440 м, l2
= F = 880 м, l3 = 2F = 1760 м.

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

Вернуться