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

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

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

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) псевдоакустический каротаж (РАК).

Раздел 3.2. Выбор способа решения прямой
динамической задачи

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

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

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

,

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

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

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