1. Engineering Research Center for Seismic Disaster Prevention and Engineering Geological Disaster Detection of Jiangxi Province, East China University of Technology, Nanchang, Jiangxi 330013, China; 2. School of Geophysics and Measurement-Control Technology, East China University of Technology, Nanchang, Jiangxi 330013, China; 3. College of Mining Engineering, Taiyuan University of Technology, Taiyuan, Shanxi 030024, China; 4. Institute of Exploration and Development, SINOPEC Shanghai Offshore Oil & Gas Company, Shanghai 200120, China
Abstract:Time-varying property of the frequency components, i.e., the time-frequency characteristic, in seismic signals plays a significant role in hydrocarbon detection. The simplest time-frequency analysis technique,Gabor transform,has been widely used in seismic data interpretation. However,the time resolution of time-frequency spectra obtained via traditio-nal Gabor transform is poor because the spectra of adjacent reflection wavelets are seriously aliased with each other. This defect is unfavorable for high-resolution seismic data interpretation. To improve the time resolution of Gabor time-frequency spectra,we first express the calculation of a Gabor time-frequency spectrum as inversion problem solving according to the definition of inverse Gabor transform. Then,a group-sparse regularization constraint is introduced with the strategy of grouping the spectra of the same moment. Finally,the objective function is solved using the projected fast iterative soft-thresholding algorithm,and the Gabor time-frequency spectrum can be obtained. In addition,an adaptive method for constructing Gaussian window functions required by Gabor time-frequency analysis is developed with the instantaneous centroid frequency of the seismic signal. Theoretical signal tests and field data applications show that the proposed method can significantly compress the energy groups of Gabor time-frequency spectra in the time direction,simultaneously improve the time resolution of both low- and high-frequency components. Therefore,the top and bottom interfaces of thin reservoirs can be clearly revealed on both low- and high-frequency sections,which facilitates fine comparison and interpretation.
DENHAM L R.Seismic interpretation[J].Procee-dings of the IEEE,1984,72(10):1255-1265.
[2]
SINHA S,ROUTH P S,ANNO P D,et al.Spectral decomposition of seismic data with continuous-wavelet transform[J].Geophysics,2005,70(6):P19-P25.
[3]
LIU J,MARFURT K J.Instantaneous spectral attri-butes to detect channels[J].Geophysics,2007,72(2):P23-P31.
[4]
李思源,徐天吉.基于Wigner-Ville分布与Chrip-Z变换的高分辨时频分析方法[J].石油地球物理勘探,2022,57(1):168-175,211.LI Siyuan,XU Tianji.A new high-resolution time-frequency analysis method based on Wigner-Ville distribution and Chrip-Z transform[J].Oil Geophysical Prospecting,2022,57(1):168-175,211.
[5]
张付瑷,陈学华,罗鑫,等.改进的窗参数优化S变换及其在河道检测中的应用[J].石油地球物理勘探,2021, 56(4):809-814,881.ZHANG Fuai,CHEN Xuehua,LUO Xin,et al.Improved window parameter optimized S-transform and its application in channel detection[J].Oil Geophysical Prospecting,2021,56(4):809-814,881.
MARFURT K J,KIRLIN R L.Narrow-band spectral analysis and thin-bed tuning[J].Geophysics,2001,66(4):1274-1283.
[8]
高静怀,陈文超,李幼铭,等.广义S变换与薄互层地震响应分析[J].地球物理学报,2003,46(4):526-532.GAO Jinghuai,CHEN Wenchao,LI Youming,et al.Generalized S transform and seismic response analysis of thin interbeds[J].Chinese Journal of Geophysics,2003,46(4):526-532.
[9]
CASTAGNA J P,SUN S,SIEGFRIED R W.Instantaneous spectral analysis:Detection of low-frequency shadows associated with hydrocarbons[J].The Lea-ding Edge,2003,22(2):120-127.
[10]
WANG Y.Seismic time-frequency spectral decomposition by matching pursuit[J].Geophysics,2007,72(1):V13-V20.
[11]
GHOLAMI A.Sparse time-frequency decomposition and some applications[J].IEEE Transactions on Geoscience and Remote Sensing,2013,51(6):3598-3604.
[12]
王小杰,印兴耀,吴国忱.基于叠前地震数据的地层Q值估计[J].石油地球物理勘探,2011,46(3):423-428.WANG Xiaojie,YIN Xingyao,WU Guochen.Estimation of stratigraphic quality factors on pre-stack seismic data[J].Oil Geophysical Prospecting,2011,46(3):423-428.
[13]
HAO Y,WEN X,ZHANG B,et al.Q estimation of seismic data using the generalized S-transform[J].Journal of Applied Geophysics,2016,135:122-134.
[14]
郝亚炬,黄捍东,文晓涛,等.广义S域Q值估计方法及其在油气检测中的应用[J].石油地球物理勘探,2017, 52(5):1059-1066.HAO Yaju,HUANG Handong,WEN Xiaotao,et al.Q estimation in the generalized S domain and its application in the hydrocarbon detection[J].Oil Geophysical Prospecting,2017,52(5):1059-1066.
[15]
ASKARI R,SIAHKOOHI H R.Ground roll attenuation using the S and x-f-k transforms[J].Geophysical Prospecting,2008,56(1):105-114.
[16]
GABOR D.Theory of communication[J].Journal of the Institution of Electrical Engineers,1946,93(26):429-441.
[17]
STOCKWELL R G,MANSINHA L,LOWE R P.Localization of the complex spectrum:The S transform[J].IEEE Transactions on Signal Processing,1996,44(4):998-1001.
[18]
MANSINHA L,STOCKWELL R G,LOWE R P,et al.Local S-spectrum analysis of 1-D and 2-D data[J].Physics of the Earth and Planetary Interiors,1997,103(3-4):329-336.
[19]
陈学华.时频域高精度储层结构分析:理论、算法及应用[D].四川成都:成都理工大学,2009.CHEN Xuehua.High Resolution Reservoir Architecture Analysis in Time Frequency Domain:Theory,Algorithm and Application[D].Chengdu University of Technology,Chengdu,Sichuan,2009.
[20]
PURYEAR C I,PORTNIAGUINE O N,COBOS C M,et al.Constrained least-squares spectral analysis:Application to seismic data[J].Geophysics,2012,77(5):V143-V167.
[21]
SATTARI H,GHOLAMI A,SIAHKOOHI H R.Seismic data analysis by adaptive sparse time-frequency decomposition[J].Geophysics,2013,78(5):V207-V217.
[22]
SATTARI H.High-resolution seismic complex trace analysis by adaptive fast sparse S-transform[J].Geophysics,2017,82(1):V51-V67.
[23]
LIU Y,ZHAN Z,CAI J,et al.Projected iterative soft-thresholding algorithm for tight frames in compressed sensing magnetic resonance imaging[J].IEEE Tran-sactions on Medical Imaging,2016,35(9):2130-2140.
[24]
HAO Y,HUANG H,LUO Y,et al.Nonstationary acoustic-impedance inversion algorithm via a novel equivalent Q-value estimation scheme and sparse regularizations[J].Geophysics,2018,83(6):R681-R698.
[25]
KOWALSKI M,TORRÉSANI B.Sparsity and persistence:mixed norms provide simple signal models with dependent coefficients[J].Signal,Image and Video Processing,2009,3(3):251-264.