Seismic Blind Deconvolution Method Based on Generalized Gaussian Distribution
LI Hai-shan,WU Guo-chen,YIN Xing-yao
DOI: https://doi.org/10.6038/j.issn.1004-2903.2012.03.014
2012-01-01
Abstract:Deconvolution is an effective and essential method to improve seismic resolution.Traditional deconvolution techniques assume that the wavelet is minimum phase and the reflection coefficients have a white Gaussian noise distribution.The assumption of a Gaussian distribution means that recovery of the true phase of the wavelet is impossible;however,a non-Gaussian distribution in theory allows recovery of the phase.It is generally recognized that primary reflection coefficients typically have a non-Gaussian amplitude distribution,and the purpose of deconvolution is to recover the non-Gaussian distribution feature of the reflectivity.Blind deconvolution differs from conventional deconvolution in that the reflection coefficients have a non-Gaussian amplitude distribution and we do not make a priori assumption about the wavelet phase.Bussgang algorithm was proposed for convolutive blind source separation problem,in this paper,it was applied to seismic blind deconvolution process.Because the generalized Gaussian probability density function has the capacity for approximating any probability density function,in accordance with the statistical characteristics of the reflectivity,it was introduced in to reflect the super-Gaussian distribution characteristics of the reflectivity.On the basis of the statistical characteristics of the reflectivity and the principle of Bussgang algorithm,a corresponding objective function which utilizing Kullback-Leibler distance as non-Gaussian measure was established,and then,the memoryless nonlinearity involved in Bussgang algorithm was derived.finally,the seismic blind deconvolution technique based on generalized Gaussian distribution was realized.The results of model test showed that the conventional deconvolution methods,such as spiking deconvolution,estimate the wavelet characteristics from the autocorrelation,unfortunately,there is no information about the phase of the wavelet in the autocorrelation,so the deconvolution result is incorrect.However,testing results showed that the deconvolution method in this paper can be better adapt to non-minimum phase systems,the seismic wavelet and reflectivity can be estimated at the same time,the phase of wavelet can be correctly recovered,and the method is robust in the presence of noise.Then,the method was used to process the field data,in contrast with the spiking deconvolution result,the blind deconvolution technique not only improved the resolution of the profile,the phase axis superimposed together were separated,but also highlights the weak reflection information between the strong reflection axis.As an effective means to improve seismic resolution,the blind deconvolution technique can broaden the effective band of seismic records,and will play an important role in seismic exploration and exploitation.