Direct inversion of P-wave to S-wave velocity ratio by Lp quasi-norm sparse constraints
ZHANG Tianyue1,2, LIN Kai1,2, WEN Xiaotao1,2, ZHAO Lian1,2, ZHANG Yuqiang1,2, LEI Yang3
1. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Chengdu University of Technology, Chengdu, Sichuan 610059, China; 2. College of Geophysics, Chengdu University of Technology, Chengdu, Sichuan 610059, China; 3. South Branch of SINOPEC Petroleum Engineering Geophysical Co., Ltd, Chengdu, Sichuan 610041, China
Abstract:The P-wave to S-wave velocity ratio(vP/vS) is a vital tool for gas reservoir identification,reservoir characterization,and lithology recognition. At present,the P-wave velocity and S-wave velocity are mainly obtained through the inversion of the reflection coefficient approximate equation,and then the vP/vS is calculated. However, this indirect calculation method creates a cumulative error. To obtain the vP/vS directly from pre-stack seismic data,this paper proposes a new generalized elastic impedance equation and further derives an approximate equation of the P-wave reflection coefficient,which is linked to the vP/vS,P-wave velocity,and density. To achieve high-precision inversion results,this paper proposes a prestack seismic inversion method based on the sparse constraint of the Lp quasi-norm utilizing the derived approximation equation of reflection coefficient,which is solved by the alternating direction multiplier algorithm. The proposed direct inversion method is applied to theoretical models and practical data and compared with the indirect inversion method. The results demonstrate that the direct inversion method exhibits higher inversion accuracy and clearer boundary characterization of gas-bearing reservoirs.
CHEN F B,ZONG Z Y,YANG Y M. Amplitude-variation-with-offset inversion using P- to S-wave velocity ratio and P-wave velocity[J]. Geophysics,2022,87(4):N63-N74.
[7]
WANG X Q,SCHUBNEL A,FORTIN J,et al. High vP/vS ratio:Saturated cracks or anisotropy effects?[J]. Geophysical Research Letters,2012,39(11):L11307.
[8]
YIN X Y,ZHANG S X. Bayesian inversion for effective pore-fluid bulk modulus based on fluid-matrix decoupled amplitude variation with offset approximation[J]. Geophysics,2014,79(5):R221-R232.
[9]
张金强,刘振峰,刘喜武,等. 致密砂岩储层自适应岩石物理建模方法[J]. 石油地球物理勘探,2021,56(4):792-800,808.ZHANG Jinqiang,LIU Zhenfeng,LIU Xiwu,et al. Self-adaptive rock physics modeling method for tight sandstone reservoirs[J]. Oil Geophysical Prospecting,2021,56(4):792-800,808.
[10]
程冰洁,徐天吉,梁群,等. 多波速度比参数含气性识别研究与应用[J]. 石油地球物理勘探,2014,49(2):307-315.CHENG Bingjie,XU Tianji,LIANG Qun,et al. Gas identification using multi-wave velocity ratio parameters[J].Oil Geophysical Prospecting,2014,49(2):307-315.
[11]
张凌远,张宏兵,尚作萍,等. 基于Zoeppritz方程的叠前和叠后混合多参数非线性地震反演[J]. 石油地球物理勘探,2021,56(1):164-171.ZHANG Lingyuan,ZHANG Hongbing,SHANG Zuoping,et al. Nonlinear multi-parameter hybrid inversion of pre-stack and post-stack seismic data based on Zoeppritz equation[J]. Oil Geophysical Prospecting,2021,56(1):164-171.
[12]
CONNOLLY P. Elastic impedance[J]. The Leading Edge,1999,18(4):438-452.
[13]
WHITCOMBE D N. Elastic impedance normalization[J]. Geophysics,2002,67(1):60-62.
[14]
WHITCOMBE D N,CONNOLLY P A,REAGAN R L,et al. Extended elastic impedance for fluid and lithology prediction[J]. Geophysics,2002,67(1):63-67.
[15]
马劲风. 地震勘探中广义弹性阻抗的正反演[J]. 地球物理学报,2003,46(1):118-124.MA Jinfeng. Forward modeling and inversion method of generalized elastic impedance in seismic exploration[J]. Chinese Journal of Geophysics, 2003, 46(1):118-124.
[16]
印兴耀,刘晓晶,吴国忱,等. 模型约束基追踪反演方法[J]. 石油物探,2016,55(1):115-122.YIN Xingyao,LIU Xiaojing, WU Guochen, et al. Basis pursuit inversion method under model constraint[J]. Geophysical Prospecting for Petroleum, 2016,55(1):115-122.
[17]
刘晓晶,印兴耀,吴国忱,等. 基于基追踪弹性阻抗反演的深部储层流体识别方法[J]. 地球物理学报,2016,59(1):277-286.LIU Xiaojing,YIN Xingyao,WU Guochen,et al. Identification of deep reservoir fluids based on basis pursuit inversion for elastic impedance[J]. Chinese Journal of Geophysics,2016,59(1):277-286.
[18]
SHE B,WANG Y J,ZHANG J S,et al. AVO inversion with high-order total variation regularization[J]. Journal of Applied Geophysics,2019,161:167-181.
WOODWORTH J,CHARTRAND R. Compressed sensing recovery via nonconvex shrinkage penalties[J]. Inverse Problems,2016,32(7):075004.
[21]
CHEN Y P,PENG Z M, GHOLAMI A, et al. Seismic signal sparse time-frequency representation by Lp-quasi norm constraint[J]. Digital Signal Processing,2019,87(4):43-59.
[22]
张雨强,文晓涛,吴昊等. 基于Lp拟范数稀疏约束和交替方向乘子算法的波阻抗反演[J]. 石油物探,2022,61(5):856-864.ZHANG Yuqiang,WEN Xiaotao,WU Hao,et al. Seismic acoustic impedance inversion using Lp quasi-norm sparse constraint and alternating direction multiplier algorithm[J]. Geophysical Prospecting for Petroleum,2022,61(5):856-864.
[23]
GOODWAY B,CHEN T W,DOWNTON J. Improved AVA fluid detection and lithology discrimination using Lamé petrophysical parameters "λρ", "μρ", & "λ/μ fluid stack" from P and S inversions[C]. SEG Technical Program Expanded Abstracts,1997,16:183-186.
[24]
AKI K,RICHARDS P G. Quantitative Seismology:Theory and Methods[M]. W H Freeman and Company,San Francisco,1980.
[25]
BOYD S,PARIKH N,CHU E,et al. Distributed optimization and statistical learning via the alternating direction method of multipliers[J]. Foundations and Trends in Machine Learning,2011,3(1):1-122.