RU2014118825A - METHOD FOR PROCESSING AND INTERPRETING SEISMIC DATA - Google Patents

Abstract

1. Способ обработки и интерпретации сейсмических данных, включающий получение временного сейсмического разреза в исследуемой толще земли, отображение сейсмического разреза в виде набора сейсмических трасс, определение спектральных и энергетических характеристик отраженных и рассеянных волн по выборкам данных на сейсмическом разрезе в перемещающемся скользящем окне и вынесение суждения об отражающих свойствах и локальных неоднородностях объектов геологической среды, отличающийся тем, что обработка данных повыделению энергии рассеянных волн сейсмического поля включает следующие последовательные операции: а) определение дисперсии по отдельным трассам временного разреза в скользящем окне; б) разложение (декомпозиция) полученных значений дисперсии на трендовую и локальную составляющие, обусловленные соответственно энергией отраженных и рассеянных волн, путем одномерной адаптивной энергетической фильтрации; в) разложение трендовой и локальной составляющих дисперсии посредством вейвлет-преобразования; г) привязка энергии вейвлет-преобразования трендовой составляющей дисперсии и энергии локальной составляющей дисперсии по времени, а интерпретация данных по энергии рассеянных волн включает суждение о трещинно-кавернозных зонах неоднородностей геологической среды, связанных с зонами коллекторов и месторождений углеводородов.2. Способ по п.1, отличающийся тем, что определение дисперсии D по отдельным трассам временного разреза - операция а) - производят в скользящем окне, по крайней мере, из трех точек путем расчета по амплитудам A(t) отраженных волн для, по крайней мере, трех дискретов по времени I) в �1. A method of processing and interpreting seismic data, including obtaining a temporary seismic section in the studied body of the earth, displaying a seismic section in the form of a set of seismic traces, determining the spectral and energy characteristics of reflected and scattered waves from data samples on a seismic section in a moving sliding window and making a judgment about reflective properties and local heterogeneities of objects of the geological environment, characterized in that the data processing by the release of energy nnyh seismic wave field comprises the following successive steps: a) determining the variance on selected tracks time section in the moving window; b) decomposition (decomposition) of the obtained dispersion values into trend and local components, due to the energy of reflected and scattered waves, respectively, by one-dimensional adaptive energy filtering; c) the expansion of the trend and local components of the variance through the wavelet transform; d) the binding of the energy of the wavelet transform of the trend component of the dispersion and the energy of the local component of the dispersion in time, and the interpretation of the data on the energy of the scattered waves includes a judgment on the fracture-cavernous zones of heterogeneities of the geological environment associated with the zones of reservoirs and hydrocarbon deposits. 2. The method according to claim 1, characterized in that the determination of the dispersion D by individual tracks of the time section — operation a) —is performed in a sliding window from at least three points by calculating the amplitudes A (t) of the reflected waves for at least , three discrete time I) in �

Claims (5)

1. Способ обработки и интерпретации сейсмических данных, включающий получение временного сейсмического разреза в исследуемой толще земли, отображение сейсмического разреза в виде набора сейсмических трасс, определение спектральных и энергетических характеристик отраженных и рассеянных волн по выборкам данных на сейсмическом разрезе в перемещающемся скользящем окне и вынесение суждения об отражающих свойствах и локальных неоднородностях объектов геологической среды, отличающийся тем, что обработка данных повыделению энергии рассеянных волн сейсмического поля включает следующие последовательные операции: а) определение дисперсии по отдельным трассам временного разреза в скользящем окне; б) разложение (декомпозиция) полученных значений дисперсии на трендовую и локальную составляющие, обусловленные соответственно энергией отраженных и рассеянных волн, путем одномерной адаптивной энергетической фильтрации; в) разложение трендовой и локальной составляющих дисперсии посредством вейвлет-преобразования; г) привязка энергии вейвлет-преобразования трендовой составляющей дисперсии и энергии локальной составляющей дисперсии по времени, а интерпретация данных по энергии рассеянных волн включает суждение о трещинно-кавернозных зонах неоднородностей геологической среды, связанных с зонами коллекторов и месторождений углеводородов.1. A method of processing and interpreting seismic data, including obtaining a temporary seismic section in the studied body of the earth, displaying a seismic section in the form of a set of seismic traces, determining the spectral and energy characteristics of reflected and scattered waves from data samples on a seismic section in a moving sliding window and making a judgment about reflective properties and local heterogeneities of objects of the geological environment, characterized in that the data processing by the release of energy nnyh seismic wave field comprises the following successive steps: a) determining the variance on selected tracks time section in the moving window; b) decomposition (decomposition) of the obtained dispersion values into trend and local components, due to the energy of reflected and scattered waves, respectively, by one-dimensional adaptive energy filtering; c) the expansion of the trend and local components of the variance through the wavelet transform; d) the binding of the energy of the wavelet transform of the trend component of the dispersion and the energy of the local component of the dispersion over time, and the interpretation of the data on the energy of the scattered waves includes a judgment on the fracture-cavernous zones of heterogeneities of the geological environment associated with the zones of reservoirs and hydrocarbon deposits. 2. Способ по п.1, отличающийся тем, что определение дисперсии D по отдельным трассам временного разреза - операция а) - производят в скользящем окне, по крайней мере, из трех точек путем расчета по амплитудам A_i(t_i) отраженных волн для, по крайней мере, трех дискретов по времени I) в соответствии с выражением2. The method according to claim 1, characterized in that the determination of the dispersion D for individual tracks of the time section — operation a) —is performed in a sliding window from at least three points by calculating the amplitudes A _i (t _i ) of the reflected waves for at least three time samples I) in accordance with the expression

где A_i - значения измеренных амплитуд в скользящем окне из n значений (i=1,2,…n);where A _i are the values of the measured amplitudes in a sliding window of n values (i = 1,2, ... n);

- среднее арифметическое значение измеренных амплитуд для n дискретных значений по времени t_i.

- the arithmetic average of the measured amplitudes for n discrete values in time t _i . 3. Способ по п.1, отличающийся тем, что разложение значений дисперсии на трендовую и локальную составляющие - операция б) - путем одномерной адаптивной энергетической фильтрации реализуют посредством расчета автокорреляционной функции, размер n_δ=1,2 r₀ базового окна адаптивного фильтра выбирают по величине радиуса r₀ корреляции автокорреляционной функции, а размер n_i текущего (i-того) окна фильтрации и весовые коэффициенты h_i фильтра определяют из уравнения3. The method according to claim 1, characterized in that the decomposition of the dispersion values into a trend and local components - operation b) - is implemented by means of a one-dimensional adaptive energy filtering by calculating the autocorrelation function, the size n _δ = 1.2 r ₀ of the adaptive filter base window is selected by the value of the radius r _{0 of the} correlation of the autocorrelation function, and the size n _{i of the} current (i-th) filter window and the weight coefficients h _{i of the} filter are determined from the equation [R(m-i)-λ_max·I(m-i)]h_i=0,[R (mi) -λ _max · I (mi)] h _i = 0, где R(m-i) - корреляционная оценка корреляционной матрицы размерностью m;where R (m-i) is the correlation estimate of the correlation matrix of dimension m; λ_max - максимальное собственное значение матрицы R(m-i);λ _max is the maximum eigenvalue of the matrix R (mi); I(m-i) - единичная матрица.I (m-i) is the identity matrix. 4. Способ по п.1, отличающийся тем, что разложение трендовой и локальной составляющих дисперсии посредством вейвлет-преобразования -операция в) - проводят в соответствии с формулой4. The method according to claim 1, characterized in that the decomposition of the trend and local components of the dispersion by means of the wavelet transform-operation c) is carried out in accordance with the formula

где W(a,b) - свертка амплитуд дисперсии D(t),характеризующая энергию волнового поля;where W (a, b) is the convolution of the dispersion amplitudes D (t), which characterizes the energy of the wave field; a, b - вещественные значения, характеризующие соответственно частоту импульса Морле-вейвлета и его перемещение вдоль трассы;a, b are real values characterizing, respectively, the frequency of the Morlet wavelet impulse and its movement along the path; ψ(t) - заданная форма вейвлета, в качестве которой используется импульс в виде Морле-вейвлета.ψ (t) is the given wavelet form, which is used as a momentum in the form of a Morlet wavelet. 5. Способ по п.1, отличающийся тем, что этапы а) - г) обработки данных реализуют с использованием цифрового программируемого процессора или персонального компьютера. 5. The method according to claim 1, characterized in that steps a) to g) of data processing are implemented using a digital programmable processor or a personal computer.