Процедура DECONSC.
Подготовка данных для поверхностно-согласованной деконволюции сейсмических трасс.
Зарегистрированный сейсмический сигнал может рассматриваться как свертка сигнала в источнике с характеристикой аппаратуры, характеристикой сейсмоприемников и откликом Земли. Задачей деконволюции является оценка этих эффектов как линейных фильтров, а затем расчет и применение обратных фильтров с целью устранения неоднородностей в пунктах возбуждения и в пунктах приема.
Концепция поверхностно-согласованной деконволюции основана на идее, что сейсмический сигнал может быть разбит на компоненты, обусловленные влиянием источника, приемника, удаления, и ОГТ.
(1)
где a – зарегистрированный сигнал, s – компонента источника r – компонента приемника l – компонента удаления, c – компонента ОГТ (отклик среды). Символ означает свертку операторов.
Нашей задачей является определение и устранение факторов s , r и l.
В спектральной области уравление (1) преобразуется в произведение спектров.
(2)
Будем предполагать, что все компоненты уравнения (2) – минимально-фазовые. Это позволяет нам избавиться от фазовых спектров и оперировать только амлитудными спектрами (автокорреляционными функциями).
После отбрасывания фазовых спектров и логарифмирования уравнения (2) получаем систему линейных уравнений
(3)
Эта система очень похожа на постановку задачи в коррекции остаточных амплитудных поправок. Отличие состоит лишь в том, что там мы оперируем числами, а здесь – векторами. Это приведет к некоторому увеличению времени счета, что для современных компьютеров не принципиально.
Параметры:
PRED - интервал предсказания в дискретах,
Nois Percent - уровень "белого шума" в процентах,
Параметр можно указывать как целое число и как число с плавающей точкой.
NPTS - длина оператора фильтрации в дискретах.
Start Time, Stop Time - временное окно на нулевом удалении в мсек, в котором вычисляется автокорреляционная функция.
Значения должны быть умножены на делитель времен, указанный в Geom2D .
Velocity - кажущаяся скорость, которая определяет зависимость верхней границы окна расчета автокорреляционной функции от удаления ПП-ПВ. Задается в м/сек.
Start Time (L)= Start Time+L/V где L - расстояние ПП от ПВ.
Если V=0, то Start Time (L)= Start Time =const
Horizon Number - Номер горизонта в файле с описанием горизонтов.
Если этот параметр не равен нулю, окно расчета автокорреляционных функций будет отслеживаться с учетом времен горизонта.
Horizon Ext - расширение файла с описанием горизонтов.
Требуется только если Horizon Number не равен 0. Содержит прослеженные в программе FINES горизонты. Обработка с указанием горизонта возможна только при обработке профилей 2-Д. В 3-Д 8-й параметр пока воспринимается как ошибка.
Процедура записывает Автокорреляционные функции всех трасс в файл "LineName.SCD". Если такой файл был создан ранее, то можно дополнить его новыми данными (обработка фрагментов одного профиля). Тип сортировки сейсмотрасс не имеет значения.
Процедура применяет к входным трассам обычную предсказывающую деконволюцию с параметрами Pred, Noise, Npts. Это позволяет легко сопоставлять результат последующей поверхностно-согласованной деконволюции с результатом обычной предсказывающей деконволюции. Если выходной набор не нужен, то запись сейсмотрасс на выходе можно подавить установкой переключателя Dummy Write в окне Run SDPS.
Примечание: При обработке данных 3-Д пока не допускается задание горизонта.
Внимание: Горизонт не обязательно должен придерживаться фазы. Он просто задает переменное по профилю окно вычисления ФАК. При этом параметры Start Time и Stop Time воспринимаются как часть окна над горизонтов (Time Above) и под горизонтом (Time Below). Величина окна при отсутствии горизонта равна разности, а при наличии горизонта – сумме 4-го и 5-го параметров.
Файл амплитудных спектров обрабатывается процедурой DECONS2.
Процедура DECONS2 рассчитывает операторы деконволюции за ПП, ПВ и Удаление и применяет их с сейсмическим трассам.
Внимание: недопустимо задавать DECONSC и DECONS2 в одном задании.