A novel nonstationary deconvolution method based on spectral modeling and variable-step sampling hyperbolic smoothing

Abstract Deconvolution is an important part of seismic processing tool for improving the resolution. One of the key assumptions made in most deconvolutional methods is that the seismic data is stationary. However, due to the anelastic absorption, the seismic data is usually nonstationary. In this paper, a novel nonstationary deconvolution approach is proposed based on spectral modeling and variable-step sampling (VSS) hyperbolic smoothing. To facilitate our method, firstly, we apply the Gabor transform to perform a time-frequency decomposition of the nonstationary seismic trace. Secondly, we estimate the source wavelet amplitude spectrum by spectral modeling. Thirdly, smoothing the Gabor magnitude spectrum of seismic data along hyperbolic paths with VSS can obtain the magnitude of the attenuation function, and can also eliminate the effect of source wavelet. Fourthly, by assuming that the source wavelet and attenuation function are minimum phase, their phases can be determined by Hilbert transform. Finally, the estimated two factors are removed by dividing them into the Gabor spectrum of the trace to estimate the Gabor spectrum of the reflectivity. An inverse Gabor transform gives the time-domain reflectivity estimate. Tests on synthetic and field data show that the presented method is an effective tool that not only has the advantages of stationary deconvolution, but also can compensate for the energy absorption, without knowing or estimating the quality factor Q.

[1]  Einar Kjartansson,et al.  Constant Q-wave propagation and attenuation , 1979 .

[2]  M. Turhan Taner,et al.  The use of the conjugate-gradient algorithm in the computation of predictive deconvolution operators , 1985 .

[3]  N. Hargreaves,et al.  Inverse Q filtering by Fourier transform , 1991 .

[4]  David C. Henley,et al.  Z-99 Gabor deconvolution: theory and practice , 2004 .

[5]  Tadeusz J. Ulrych,et al.  Seismic absorption compensation: A least squares inverse scheme , 2007 .

[6]  David C. Henley,et al.  Gabor deconvolution: Estimating reflectivity by nonstationary deconvolution of seismic data , 2011 .

[7]  Gary F. Margrave,et al.  Theory of nonstationary linear filtering in the Fourier domain with application to time‐variant filtering , 1998 .

[8]  Tadeusz J. Ulrych,et al.  Processing via spectral modeling , 1991 .

[9]  L. J. Griffiths,et al.  Adaptive deconvolution; a new technique for processing time-varying seismic data , 1977 .

[10]  G. Margrave,et al.  Phase correction in Gabor deconvolution , 2005 .

[11]  Yanghua Wang,et al.  Quantifying the effectiveness of stabilized inverse Q filtering , 2003 .

[12]  John W. C. Robinson,et al.  A technique for the continuous representation of dispersion in seismic data , 1979 .

[13]  Walter I. Futterman,et al.  Dispersive body waves , 1962 .

[14]  Dennis Gabor,et al.  Theory of communication , 1946 .

[15]  G. Clarke Time-Varying Deconvolution Filters , 1968 .

[16]  Yin Ji-feng Two key technologies about exploration of stratigraphic/lithological reservoirs , 2004 .

[17]  Dave Hale,et al.  Q‐adaptive deconvolution , 1982 .

[18]  Yanghua Wang,et al.  Inverse Q-filter for seismic resolution enhancement , 2006 .

[19]  Shoudong Wang Attenuation compensation method based on inversion , 2011 .

[20]  Yanghua Wang,et al.  A stable and efficient approach of inverse Q filtering , 2002 .

[21]  R. Stolt MIGRATION BY FOURIER TRANSFORM , 1978 .