Relative non-high-frequency leakage equivalent N-drop smooth spectrum analog deconvolution method

Abstract

The invention belongs to the seismic data processing and relates to a method by utilizing an equivalent N point smooth spectrum for the simulation of deconvolution. The implementation steps are that the seismic data is collected; the Fourier transform is made for the seismic data at a certain time window; a low pass filter operator of the time domain with an interception frequency of FE and a sampling rate of DT is calculated, which is taken as the equivalent N point smooth operator; FE and DT satisfy the equation that FE*DT is equal to X*DF*1.024; the estimate value of the amplitude spectrum of the seismic wavelet and the reverse wavelet are obtained, which are convoluted with the original seismic data; and the deconvolution seismic data of the time window is obtained. The single peak hypothesis is not applied for the amplitude spectrum of the seismic wavelet in the invention, which suits for complex amplitude spectrum of the seismic wavelet; compared with the least square polynomial fitting, the shape of the simulated amplitude spectrum curve of the seismic record is not affected by the global influence; the invention has the advantages of high calculation efficiency, better stability and improved resolution of seismic profiles.

Description

technical field [0001] The invention relates to geophysical exploration technology, belongs to the deconvolution of seismic data processing, in particular to an equivalent N-point smoothing spectrum simulation deconvolution method capable of obtaining high-resolution seismic sections with relatively no high-frequency leakage. Background technique [0002] Geophysical exploration needs to further process seismic data after obtaining observation results through conventional methods. Generally, there are three main processing procedures: deconvolution, stacking, and migration. The seismic profile is actually the result of the convolution of the subsurface reflection coefficient sequence and the seismic wavelet. Since the seismic wavelet has a certain length of time and is not a sharp pulse, its convolution effect makes the reflection layer (or event axis) displayed on the seismic section thicker, and the reflection layers adjacent to each other in time are connected or overlapp...

AI Technical Summary

Problems solved by technology

[0006] The conventional spectral simulation deconvolution method limits the seismic wavelet amplitude spectrum to a single-peak curve, which is not enough to describe the complexity of the seismic wavelet amplitude spectrum; the polynomial least squares fitting algorithm has statistical characteristics, so that the shape of the final fitted curve depends on For each part of the original data, the curve shapes of different frequency bands are related to each other; solving the polynomial needs to solve the matrix equation, the algorithm is relatively complicated, the data sensitivity is strong, and it is difficult to obtain high-resolution seismic sections

Images