Abstract The Gabor and S transforms are frequently used in time-frequency decomposition methods. Constrained by the uncertainty principle, both transforms produce low-resolution time-frequency decomposition results in the time and frequency domains. To improve the resolution of the time-frequency decomposition results, we use the instantaneous frequency distribution function (IFDF) to express the seismic signal. When the instantaneous frequencies of the nonstationary signal satisfy the requirements of the uncertainty principle, the support of IFDF is just the support of the amplitude ridges in the signal obtained using the short-time Fourier transform. Based on this feature, we propose a new iteration algorithm to achieve the sparse time-frequency decomposition of the signal. The iteration algorithm uses the support of the amplitude ridges of the residual signal obtained with the short-time Fourier transform to update the time-frequency components of the signal. The summation of the updated time-frequency components in each iteration is the result of the sparse time-frequency decomposition. Numerical examples show that the proposed method improves the resolution of the time-frequency decomposition results and the accuracy of the analysis of the nonstationary signal. We also use the proposed method to attenuate the ground roll of field seismic data with good results.
This research is funded by the National Basic Research Program of China (973 Program) (No. 2011 CB201002), the National Natural Science Foundation of China (No. 41374117) and the great and special projects (2011ZX05005–005 - 008HZ and 2011ZX05006 - 002).
Cite this article:
WANG Xiong-Wen,WANG Hua-Zhong. Application of sparse time-frequency decomposition to seismic data[J]. APPLIED GEOPHYSICS, 2014, 11(4): 447-458.
[1]
Askari, R., and Siahkoohi, H. R., 2008, Ground roll attenuation using the S and x-f-k transforms: Geophysical Prospecting, 56, 105-114.
[2]
Bonar, D. C., and Sacchi, M. D., 2010, Complex spectral decomposition via inversion strategies: 80th SEG Annual International Meeting, Expanded Abstracts, 1408-1412.
[3]
Castagna, J., Sun, S., and Siegfried, R., 2003, Instantaneous spectral analysis: Detection of low-frequency shadows associated with hydrocarbons: The Leading Edge, 22, 120-127C.
[4]
Castagna, J., and Sun, S., 2006, Comparison of spectral decomposition methods: First Break, 24, 75-79.
[5]
Chuang, H., and Lawton, D. C., 1991, Some properties of thin beds: 61st Annual International Meeting, SEG, Expanded Abstracts, 224-227.
[6]
Ebrom, D., 2004, The low-frequency gas shadow on seismic sections: The Leading Edge, 23, 772,
[7]
Gao, J. H., 2003, Generalized S transform and seismic response analysis of thin interbeds: Chinese J. of Geophysics, 46(4), 526-532.
[8]
Gabor, D., 1946, Theory of communication: Journal of the Institute of Electrical Engineers, 93, 429-457.
[9]
Gholami, 2012, Sparse time-frequency decomposition and some applications: IEEE Transactions on Geoscience and Remote Sensing, 51, 3598-3604,
[10]
Huang, N. E., Shen, Z., Long, S. R., Wu, M. C., Shih, H. H., Zheng, Q., Yen, N. C., Tung, C. C., and Liu, H. H., 1998. The empirical mode decomposition and the Hilbert Spectrum for Nonlinear and Nonstationary Time Series Analysis: Proceedings of the Royal Society of London A, 454, 903-995.
[11]
Huang, N. E., and Shen, S. S. P., 2005,Hilbert-Huang transform and its applications: World Scientific,.ondon.
[12]
Li, X. Y., Chen, S. M., Wang, J. M., Pei, J. Y., and Wang, Y. B., 2012, Forward modeling studies on the time-frequency characteristics of thin layers: Chinese J. of Geophysics, 55(10), 3410-3419.
[13]
Liu, G., Fomel, S., and Chen, X., 2009, Time-frequency characterization of seismic data using local attributes: 79th SEG Annual International Meeting, Expanded Abstracts, 1825-1829.
[14]
Liu, J., and Marfurt, K. J., 2006, Thin bed thickness prediction using peak instantaneous frequency: 76th SEG Annual International Meeting, Expanded Abstracts, 968-972.
[15]
Liu, Y., and Fomel, S., 2010, Local time-frequency transform and its application to ground-roll noise attenuation: 80th SEG Annual International Meeting, Expanded Abstracts, 3711-3716.
[16]
Mallat, S., 2009, A wavelet tour of signal processing-The sparse way (Third Edition).
[17]
Marfurt, K. J., and Kirlin, R. L., 2001, Narrow-band spectral analysis and thin-bed tuning: Geophysics, 66, 1274-1283
[18]
Partyka, G., Gridley, A. J., and Lopez, J., 1999, Interpretational applications of spectral decomposition in reservoir characterization: The Leading Edge, 18, 353-360,
[19]
Ren, H., Goloshubin, G., and Hilterman, F., 2007, Spectra cross plot: 77th Annual International Meeting, SEG, Expanded Abstracts, 199-203.
[20]
Rioul, O., and Vetterli, M., 1991, Wavelets and signal processing: IEEE Signal Processing Magazine, 8(4), 14-38,
[21]
Sinha, S., Routh, P. S., Anno, P. D., and Castagna, J. P., 2005, Spectral decomposition of seismic data with continuous-wavelet transform: Geophysics, 70(6), 19-25,
[22]
Steeghs, P., and Drijkoningen, G., 2001, Seismic sequence analysis and attribute extraction using quadratic time-frequency representations: Geophysics, 66, 1947-1959,
[23]
Stockwell, R.G., Mansinha, L., and Lowe, R.P., 1996, Localization of the complex spectrum: the S transform: IEEE Transactions on Signal Processing, 44(4), 998-1001.
2014, Vol. 11 
