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

Рис. 2 Ступенчатая функция

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

Коррекция кинематических поправок

Скорректированные кинематические поправки определяют на основе разновременного анализа сейсмограммы ОГТ по вееру гипербол (парабол), пересекающихся на трассе с абсциссой . Схематически процесс такого анализа сводится к следующему. Задают набор из значений , в пределах которого заключено значение фиктивной скорости для искомой волны. Для каждого значения по формуле:

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

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

Рассматриваемый принцип принят за основу в различных рабочих алгоритмах, различающихся использованием более помехоустойчивых, нежели простое суммирование, операторов разновременного криволинейного анализа, который с целью минимизации времени счета выполняют не непрерывно повремени, как это делают при получении суммоленты, а в дискретных точках (вертикальные спектры) с шагом (обычно через 25—100 мс), либо на фиксированных для каждой сейсмограммы ОГТ временах , соответствующих линии горизонта на временном разрезе (горизонтальные спектры). Вместе с тем, все разновидности способов определения скорректированных кинематических поправок основаны на разновременном анализе по вееру гипербол (парабол).

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

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

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

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

Третий способ заключается в разновременном анализе по набору парабол второй степени. Если в сейсмограмму ОГТ были введены исходные кинематические поправки, годограф отраженной волны трансформируется в параболу второй степени. Следовательно, разновременной анализ сейсмограммы, скорректированной исходными (априорными) кинематическими поправками, может быть реализован по вееру парабол второй степени. Веер парабол задают таким образом, чтобы обеспечивался постоянный сдвиг на крайнем канале сейсмограммы между двумя соседними линиями анализа. Данный способ имеет серьезные ограничения, связанные с искажением сигналов в процессе ввода исходных кинематических поправок.

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

Определение кинематических поправок по профилю

Густота точек по профилю (площади), для которых находятся скорректированные кинематические поправки, определяется сейсмогеологическими условиями изучаемого района. В сложных условиях (криволинейные несогласно залегающие границы раздела, большие углы падения) скорость может резко менять свою величину на сравнительно небольших участках профиля. В этом случае шаг по профилю между точками OГT, для которых находится зависимость сокращается до 0,5—1,0 км. В более простых ситуациях шаг может равняться 1,5— 3,0 км. На основании найденных в дискретных точках профиля значений строят сеточную модель кинематических поправок, представляющую собой таблицу . В промежуточных точках по осям и значения определяют на основании линейной интерполяции.

Практическая часть

Пакет RadExPro+ предназначен для обработки многоканальных сейсмоакустических данных на компьютерах, работающих под управлением операционной системы MS Windows. По структуре и интерфейсу пакет близок к наиболее распространенным пакетам обработки, таким как PROMAX, GEOVECTUR, IXL, OMEGA и пр. Включенные в его состав процедуры позволяют осуществлять основные операции с данными, характерные для систем обработки данных:

• ввод данных, записанных в различных форматах, включая произвольный, задаваемый пользователем. При этом реализованы дополнительные возможности, позволяющие автоматически определять количество трасс в файле и шаг по профилю;

• интерполяцию данных на регулярную сеть наблюдений, что может быть использовано, например, для объединения различных файлов данных в один профиль;

• математическую обработку и анализ данных;

• визуализацию результатов с широким набором возможностей;

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

Параметры для процедуры NMO/NMI

Модуль позволяет, используя функцию скорости, рассчитать и ввести кинематические поправки в отсчеты трасс ОГТ путем линейной интерполяции. При активации модуля появляется окно, содержащее две вкладки: Velocity и NMO.

На вкладке NMO задаются параметры расчета кинематических поправок:

• NMO — выберите эту опцию, если необходимо ввести кинематические поправки в отсчеты трасс.

• NMI — выберите эту опцию, если необходимо произвести скоростную инверсию, т.е. к сейсмограммам с введенными кинематическими поправками применить обратный кинематический закон.

• Mute percent — параметр мьютинга в процентах. Растяжение трасс после применения NMO является нежелательным, но неизбежным результатом. Необходимо задать параметр мьютинга в процентах, для того, чтобы все данные, растянувшиеся более, чем на заданное количество процентов, были обнулены.

На вкладке Velocity сгруппированы параметры задания скоростного закона:

Скоростной закон можно задать тремя способами:

• активировав опцию Single velocity function (скоростной постоянный закон), задать его вручную. Порядок записи следующий:

время:скорость,время-время:скорость и т.д.

Скорости здесь задаются в м/с.

• активировав опцию Get from file, задать скоростной закон из файла. Для этого, нажав Browse…, выберите нужный файл в открывшемся стандартном диалоговом окне.

• активировав опцию Database — picks, задать скоростной закон, предварительно сохраненный в базе данных проекта. Для этого, нажав кнопку Browse…, выберите нужный объект базы данных в открывшемся стандартном диалоговом окне.

В поле Velocity type (тип скорости) необходимо указать тип скорости: