Процедура INTPRL
Процедура предназначена для выбора способа интерполяции трасс с целью качественной миграции исходных сейсмограмм процедурой MIGLPW.
Также возможно с помощью этой процедуры "восстанавливать" пропущенные Пункты возбуджения. Но следует иметь в виду, что такое применение процедуры является подлогом и может представлять исключительно любознательнй интерес.
На вход подаются сейсмограммы равных удалений. В 2-Д это сортировка LP. В 3-Д - сортировка LP (вдоль линий InLine) и CP (вдоль линий CrossLine)
Весьма желателен ввод всех статических поправок.
Параметры программы:
Mode – Способ интерполяции.
0 – Без интерполяции. Пропущенные трассы заполняются нулями.
1 - Линейная интерполяция пропущенных ПВ.
2 – Билинейная интерполяция
3 – Интерполяция в FX области.
Система уравнений решается методом исключений Гаусса.
4 - Интерполяция в FX области.
Система уравнений решается рекурсией Левинсона в комплексных числах.
Метод срабатывает намного быстрее, чем при Mode Intprl =3
5 - Интерполяция в FK- области.
При такой интерполяции программа запрашивает дополнительные параметры.
Внимание: В 3-Д эти три параметра (StartSP, StopSP и StepSP) рассчитываются автоматически. Параметры пользователя игнорируются.
При интерполяции в FX области требуются дополнительные парамеры.
Taper (ms)– длина конусного среза на концах трассы перед преобразованием Фурье в мсек.
Рекомендуется 1/3 длины окна Twlen.
Fmin (hz)– Минимальная частота (Гц) для обработки. Все частоты, меньшие Fmin, будут подавлены.
Рекомендуется 6-10 гц.
Fmax (hz)– Максимальная частота (Гц) для обработки. Все частоты, большие Fmax, будут подавлены.
Рекомендуется Fmax=0.6 частоты Найквиста =0.6/(2*dt)
Twlen (ms)– Окно для расчета фильтров. Рекомендуется не менее 300 мсек.
Должно быть не короче вертикально приращения времени для максимального наклона в пределах оператора фильтра Ntrf.
Если все горизонты имеют одинаковые наклоны, Это окно может быть равно полной длине сейсмической трассы.
Ntrw – Размер окна расчета фильтров по пространственной координате в трассах.
Это значение должно быть большим, чем число точек фильтра Ntrf.
Также, это значение не должно превышать максимальное количество трасс в сейсмограмме равных удалений.
Ntrf – Число пространственных точек (трасс) в комплексном фильтре предсказания.
Должно превышать размеры максимальной «дыры» пикетов ПВ на профиле.
Adding Noise Level - параметр регуляризации.
Определяет уровень добавляемого белого шума. При нулевом значении этого параметра возможно появление вырожденных автокорреляционных матриц.
Дополнительные параметры при интерполяции в области FK.
Start SP – начальный ПВ в метрах. Обычно это наименьший пикет ПВ.
Stop SP – конечный ПВ в метрах. Обычно это наибольший пикет ПВ.
Step SP – Желаемый шаг ПВ. Обычно равет шагу ПП.
Заголовки "воccтановленных" трасс формируются следующим образом:
Полевой номер сейсмограммы - наращивается от максимального значения во входном сейсмическом наборе данных.
Порядковый номер сейсмограммы - наращивается от максимального значения во входном сейсмическом наборе данных.
Порядковый номер ПВ - наращивается от максимального значения во входном сейсмическом наборе данных.
Удаление ПП-ПВ - Offset берется из паспорта профиля как удаление с номером текущего удаления.
Пикет ПВ наращивается от левого актуального значения до правого с шагом Step SP.
Пикет ПП PKпп=PKпв + Offset.
Пикет ОГТ PKогт=(PKпп=PKпв)/2
Такой подход позволяет получать по данным с интерполяцией суммы ОГТ, ОТВ и ОТП.
Пример интерполяции в F-X области.
