Spectral decomposition of seismic data with reassigned smoothed pseudo Wigner–Ville distribution

Abstract Seismic signals are nonstationary mainly due to absorption and attenuation of seismic energy in strata. Referring to spectral decomposition of seismic data, the conventional method using short-time Fourier transform (STFT) limits temporal and spectral resolution by a predefined window length. Continuous-wavelet transform (CWT) uses dilation and translation of a wavelet to produce a time-scale map. However, the wavelets utilized should be orthogonal in order to obtain a satisfactory resolution. The less applied, Wigner–Ville distribution (WVD) being superior in energy distribution concentration, is confronted with cross-terms interference (CTI) when signals are multi-component. In order to reduce the impact of CTI, Cohen class uses kernel function as low-pass filter. Nevertheless it also weakens energy concentration of auto-terms. In this paper, we employ smoothed pseudo Wigner–Ville distribution (SPWVD) with Gauss kernel function to reduce CTI in time and frequency domain, then reassign values of SPWVD (called reassigned SPWVD) according to the center of gravity of the considering energy region so that distribution concentration is maintained simultaneously. We conduct the method above on a multi-component synthetic seismic record and compare with STFT and CWT spectra. Two field examples reveal that RSPWVD potentially can be applied to detect low-frequency shadows caused by hydrocarbons and to delineate the space distribution of abnormal geological body more precisely.

[1]  Naim Baydar,et al.  A comparative study of acoustic and vibration signals in detection of gear failures using Wigner-Ville distribution. , 2001 .

[2]  K. Kodera,et al.  A new method for the numerical analysis of nonstationary signals , 1976 .

[3]  J. E. Moyal Quantum mechanics as a statistical theory , 1949, Mathematical Proceedings of the Cambridge Philosophical Society.

[4]  Patrice Abry,et al.  Wavelet-based spectral analysis of 1/f processes , 1993, 1993 IEEE International Conference on Acoustics, Speech, and Signal Processing.

[5]  P. Xu,et al.  Crustal Velocity Structure of the Deep Continental Subduction Zone — A Wide Angle Reflection/Refraction Seismic Study on the Eastern Dabie Orogen , 2003 .

[6]  Dan Ebrom,et al.  The low-frequency gas shadow on seismic sections , 2004 .

[7]  F. Hlawatsch,et al.  Linear and quadratic time-frequency signal representations , 1992, IEEE Signal Processing Magazine.

[8]  M. Rauch-Davies,et al.  Spectral Decomposition – Transform Methods and Fluid and Reservoir Prediction Case Study , 2005 .

[9]  T. Gang,et al.  Studying surface structure in desert areas using multiple kinds of surface wave data , 2007 .

[10]  Xiaodong Zheng,et al.  Spectral decomposition using Wigner-Ville distribution with applications to carbonate reservoir characterization , 2008 .

[11]  M. Taner,et al.  Complex seismic trace analysis , 1979 .

[12]  William J. Williams,et al.  Improved time-frequency representation of multicomponent signals using exponential kernels , 1989, IEEE Trans. Acoust. Speech Signal Process..

[13]  Patrick Flandrin,et al.  Improving the readability of time-frequency and time-scale representations by the reassignment method , 1995, IEEE Trans. Signal Process..

[14]  Kurt J. Marfurt,et al.  Matching Pursuit Decomposition Using Morlet Wavelets , 2005 .

[15]  C. Robert Pinnegar,et al.  The S-transform with windows of arbitrary and varying shape , 2003 .

[16]  Xingong Li,et al.  Time‐Frequency decomposition based on Ricker wavelet , 2004 .

[17]  MansinhaL. Localization of the complex spectrum , 1996 .

[18]  G. Partyka,et al.  Interpretational applications of spectral decomposition in reservoir characterization , 1999 .

[19]  Yanghua Wang,et al.  Seismic time-frequency spectral decomposition by matching pursuit , 2007 .

[20]  O. Rioul,et al.  Wavelets and signal processing , 1991, IEEE Signal Processing Magazine.

[21]  Stéphane Mallat,et al.  Matching pursuits with time-frequency dictionaries , 1993, IEEE Trans. Signal Process..

[22]  Jerome Mars,et al.  Improving surface-wave group velocity measurements by energy reassignment , 2003 .

[23]  Xiao-Dong Zheng,et al.  Wigner-Ville distribution and its application in seismic attenuation estimation , 2007 .

[24]  Khaled H. Hamed,et al.  Time-frequency analysis , 2003 .

[25]  Keith Worden,et al.  TIME–FREQUENCY ANALYSIS IN GEARBOX FAULT DETECTION USING THE WIGNER–VILLE DISTRIBUTION AND PATTERN RECOGNITION , 1997 .

[26]  Bin Wang,et al.  Kirchhoff Migration of Seismic Data Compressed By Matching Pursuit Decomposition , 1996 .

[27]  Jechang Jeong,et al.  Kernel design for reduced interference distributions , 1992, IEEE Trans. Signal Process..

[28]  P. Anno,et al.  Spectral decomposition of seismic data with continuous-wavelet transform , 2005 .

[29]  L. Cohen,et al.  Time-frequency distributions-a review , 1989, Proc. IEEE.

[30]  Jing-Hua Gao,et al.  Generalized S Transform and Seismic Response Analysis of Thin Interbedss Surrounding Regions by Gps , 2003 .

[31]  Lalu Mansinha,et al.  Localization of the complex spectrum: the S transform , 1996, IEEE Trans. Signal Process..

[32]  E. Wigner On the quantum correction for thermodynamic equilibrium , 1932 .