СПОСОБ ОБРАБОТКИ И ИНТЕРПРЕТАЦИЙ СЕЙСМИЧЕСКИХ ДАННЫХ

1. Способ обработки и интерпретации сейсмических данных, включающий получение временного сейсмического разреза в исследуемой толще земли, отображение сейсмического разреза в виде набора сейсмических трасс, определение спектральных и энергетических характеристик отраженных и рассеянных волн по выборкам данных на сейсмическом разрезе в перемещающемся скользящем окне и вынесение суждения об отражающих свойствах и локальных неоднородностях объектов геологической среды, отличающийся тем, что обработка данных повыделению энергии рассеянных волн сейсмического поля включает следующие последовательные операции: а) определение дисперсии по отдельным трассам временного разреза в скользящем окне; б) разложение (декомпозиция) полученных значений дисперсии на трендовую и локальную составляющие, обусловленные соответственно энергией отраженных и рассеянных волн, путем одномерной адаптивной энергетической фильтрации; в) разложение трендовой и локальной составляющих дисперсии посредством вейвлет-преобразования; г) привязка энергии вейвлет-преобразования трендовой составляющей дисперсии и энергии локальной составляющей дисперсии по времени, а интерпретация данных по энергии рассеянных волн включает суждение о трещинно-кавернозных зонах неоднородностей геологической среды, связанных с зонами коллекторов и месторождений углеводородов.
i
i
2. Способ по п.1, отличающийся тем, что определение дисперсии D по отдельным трассам временного разреза - операция а) - производят в скользящем окне, по крайней мере, из трех точек путем расчета по амплитудам A(t) отраженных волн для, по крайней мере, трех дискретов по времени I) в соответствии с выражением
i
где A - значения измеренных амплитуд в скользящем окне из n значений (i=1,2,…n);
i
- среднее арифметическое значение измеренных амплитуд для n дискретных значений по времени t.
δ
0
0
i
i
3. Способ по п.1, отличающийся тем, что разложение значений дисперсии на трендовую и локальную составляющие - операция б) - путем одномерной адаптивной энергетической фильтрации реализуют посредством расчета автокорреляционной функции, размер n=1,2 r базового окна адаптивного фильтра выбирают по величине радиуса r корреляции автокорреляционной функции, а размер n текущего (i-того) окна фильтрации и весовые коэффициенты h фильтра определяют из уравнения
max
i
[R(m-i)-λ·I(m-i)]h=0,
где R(m-i) - корреляционная оценка корреляционной матрицы размерностью m;
max
λ - максимальное собственное значение матрицы R(m-i);
I(m-i) - единичная матрица.
4. Способ по п.1, отличающийся тем, что разложение трендовой и локальной составляющих дисперсии посредством вейвлет-преобразования -операция в) - проводят в соответствии с формулой
где W(a,b) - свертка амплитуд дисперсии D(t),характеризующая энергию волнового поля;
a, b - вещественные значения, характеризующие соответственно частоту импульса Морле-вейвлета и его перемещение вдоль трассы;
ψ(t) - заданная форма вейвлета, в качестве которой используется импульс в виде Морле-вейвлета.
5. Способ по п.1, отличающийся тем, что этапы а) - г) обработки данных реализуют с использованием цифрового программируемого процессора или персонального компьютера.