Double deconvolution method in time-frequency domain for the non-stationary seismogram

The conventional seismic deconvolution methods almost assume that the seismogram is stationary, and the seismic wavelet is of zero-phase or minimum-phase. However, the actual seismic data can't meet the above assumptions. Thus, this paper proposes a double deconvolution method in time-frequency domain to improve the resolution of the non-stationary seismogram. Firstly, quadratic spectrum modeling method combined with higher order bispectrum method is used for the non-stationary seismogram to extract a single mixed-phase wavelet, and the first deconvolution is implemented by applying the spectrum division method for the entire of seismogram; then quadratic spectrum modeling method combined with higher order bispectrum in time-frequency domain is used to extract time-variant mixed-phase wavelet, and the second deconvolution is implemented by applying spectrum division method on every point spectrum in the time-frequency domain. Simulation experiments prove that the proposed method effectively overcomes the interference of adjacent strata, and greatly improve the resolution of the non-stationary seismogram.