Quality factor Q estimation based on S transform andvariational method
XU Li'nan1, GAO Jinghuai1, YANG Yang1, GAO Zhaoqi1, WANG Qian2
1. School of Information and Communications Engineering, Xi'an Jiaotong University, Xi'an, Shaanxi 710049, China; 2. School of Mathematics and Statistics, Hubei University of Arts and Science, Xiangyang, Hubei 441053, China
Abstract:Quality factor Q is an important parameter to quantitatively describe the viscoelastic attenuation. Accurate Q estimation is beneficial to reservoir identification and hydrocarbon detection. It can also be used for inverse Q filtering to improve the resolution of seismic data. The traditional Q estimation methods include logarithmic spectrum ratio (LSR), center frequency shift (CFS), and peak frequency shift (PFS), etc. LSR has poor noise immunity. Both CFS and PFS depend on the type of seismic wavelet. In response to these problems, this paper proposes a robust Q estimation method based on the S transform and variational method. Firstly, by studying the nonstationary convolution model, we derived the approximate representation of nonstationary seismic data in the S domain. Seco-ndly, on the basis of approximate representation, the objective function of quality factor Q and seismic wavelet is established and minimized based on the variational method, thereby obtaining the expression of Q estimation. Finally, we designed an adaptive selection scheme of an integral interval to improve the accuracy and noise resistance of the method. This scheme can automatically calculate the frequency parameters of the integration area based on the time-frequency spectrum of seismic data. Model examples demonstrate that the proposed method does not rely on the wavelet type and the length of the window function and shows good robustness to noise. The real data further verifies the effectiveness of the method in Q estimation.
RICKER N.The form and laws of propagation of seismic wavelets[J].Geophysics, 1953, 18(1):10-40.
[2]
马昭君, 刘洋. 地震波衰减反演研究综述[J]. 地球物理学进展, 2005, 20(4):1074-1081.MA Zhaojun, LIU Yang.A summary of research on seismic attenuation[J].Progress in Geophysics, 2005, 20(4):1074-1081.
[3]
苑书金, 董宁, 于常青.地震波衰减技术在鄂尔多斯盆地储层预测中的应用[J].石油物探, 2006, 45(2):182-185.YUAN Shujin, DONG Ning, YU Changqing.Application of seismic wave attenuation technology in reservoir prediction in Ordos Basin[J].Geophysical Prospecting for Petroleum, 2006, 45(2):182-185.
[4]
苏勤, 曾华会, 田彦灿, 等. 表层Q值确定性求取与空变补偿方法[J]. 石油地球物理勘探, 2019, 54(5):988-996.SU Qin, ZENG Huahui, TIAN Yancan, et al. Near-surface Q value estimation and quantitative amplitude compensation[J]. Oil Geophysical Prospecting, 2019, 54(5):988-996.
[5]
WANG Y H. A stable and efficient approach to inverse Q filtering[J].Geophysics, 2002, 67(2):657-663.
[6]
YANG Y, GAO J H, ZHANG G W, et al. An efficient attenuation correction method using the synchrosqueezing transform[C].SEG Technical Program Expanded Abstracts, 2017, 3376-3380.
[7]
赵秋芳, 云美厚, 朱丽波, 等. 近地表Q值测试方法研究进展与展望[J]. 石油地球物理勘探, 2019, 54(6):1397-1418.ZHAO Qiufang, YUN Meihou, ZHU Libo, et al. Progress and outlook of near-surface quality factor Q measurement and inversion[J]. Oil Geophysical Prospecting, 2019, 54(6):1397-1418.
[8]
LUO C, HUANG G T, LI X Y, et al. Q estimation by combining ISD with LSR method based on shaping-regularized inversion[J].IEEE Geoscience and Remote Sensing Letters, 2019, 16(9):1457-1461.
[9]
QUAN Y, HARRIS J M. Seismic attenuation tomography using the frequency shift method[J].Geophysics, 1997, 62(3):895-905.
[10]
高静怀, 杨森林, 王大兴. 利用VSP资料直达波的包络峰值处瞬时频率提取介质品质因子[J]. 地球物理学报, 2008, 51(3):853-861.GAO Jinghuai, YANG Senlin, WANG Daxing. Quality factor extraction using instantaneous frequency at envelope peak of direct waves of VSP data[J].Chinese Journal of Geophysics, 2008, 51(3):853-861.
[11]
GROSSMAN J P, MARGRAVE G F, LAMOUREUX M P, et al. A robust algorithm for constant-Q wavelet estimation using Gabor analysis[C].SEG Technical Program Expanded Abstracts, 2002, 21:2210-2213.
[12]
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.
[13]
高静怀, 张兵.基于时频滤波的吸收衰减参数估算[J].石油地球物理勘探, 2012, 47(6):931-936.GAO Jinghuai, ZHANG Bing. Seismic attenuation parameter estimation based on the time-frequency filtering[J]. Oil Geophysical Prospecting, 2012, 47(6):931-936.
[14]
MARGRAVE G F, LAMOUREUX M P, HENLEY D C. Gabor deconvolution:Estimating reflectivity by nonstationary deconvolution of seismic data[J].Geophysics, 2011, 76(3):W15-W30.
[15]
KOLSKY H. The propagation of stress pulses in visco-elastic solids[J].Philosophical Magazine, 1956, 1(8):693-710.
[16]
FUTTERMAN W. Dispersive body waves[J].Journal of Geophysical Research, 1962, 67(4):5279-5291.
[17]
DONOHO D L. Ideal spatial adaptation by wavelet shrinkage[J].Biometrika, 1994, 81(3):425-455.
[18]
夏洪瑞, 葛川庆, 彭涛. 小波时空变阈值去噪方法在可控震源资料处理中的应用[J]. 石油地球物理勘探, 2010, 45(1):23-27.XIA Hongrui, GE Chuanqing, PENG Tao. Application of wavelet time-space-varying threshold denoising method in vibroseis seismic data processing[J]. Oil Geophysical Prospecting, 2010, 45(1):23-27.
[19]
WANG Y H.Generalized seismic wavelets[J].Geo-physical Journal International, 2015, 203(2):1172-1178.
[20]
WANG L L, GAO J H, ZHANG W, et al. Enhancing resolution of nonstationary seismic data by molecular-Gabor transform[J]. Geophysics, 2013, 78(1):V31-V41.