Comparative analysis of seismic P-wave and S-wave inversion capabilities
DAI Fucai1,2, ZHANG Feng1,2, QIAN Zhongping3,4, SUN Pengyuan3,4, ZOU Zhen3,4, LI Xiangyang1,2,3
1. State Key Laboratory of Oil and Gas Resources and Exploration, China University of Petroleum (Beijing), Beijing 102249, China; 2. CNPC Key Laboratory of Geophysical Exploration, China University of Petroleum (Beijing), Beijing 102249, China; 3. BGP Inc., CNPC, Zhuozhou, Hebei 072750, China; 4. National Engineering Research Center for Oil and Gas Exploration Computer Software, Beijing 100088, China
Abstract:At present, seismic exploration of oil and gas mainly uses the information of seismic P-wave, which has high uncertainty in the fields of multi-parameter inversion, complex and unconventional oil and gas reservoir prediction, and fluid identification. In recent years, the development of S-wave (SV-SV wave and SH-SH wave) vibroseis technology has greatly improved the quality of S-wave data, providing necessary conditions for the development and wide application of seismic S-wave exploration. However, compared with the traditional PP wave inversion and PP-PS wave joint inversion, the inversion ability of seismic S-wave is not clear. This paper first compares and analyzes the accuracy of the approximate formula of the reflection coefficient of P-wave, PSV wave, and S-wave represented by different elastic parameters and then uses the approximate formula of reflection coefficient to construct the linear inversion framework for different elastic parameters. Combining the inversion covariance matrix and the inversion condition number, this paper compares and analyzes the uncertainty and ill-posed nature of seismic P-wave and S-wave in the inversion of different elastic parameters. The results show that the uncertainty of shear modulus and density inversion using PP wave is the largest, and the uncertainty of PP wave inversion is still high even if using large-angle data. PSV wave inversion can effectively reduce the ill-posed nature of inversion, but it is difficult to be independently applied to actual data inversion. The joint inversion of PP wave and PSV wave can effectively reduce the uncertainty and ill-posed nature of inversion, but the joint inversion effect depends on the matching accuracy of P-wave and S-wave. SH-SH wave and SV-SV wave have obvious advantages in inversion of S-wave impedance, S-wave velocity, shear modulus, and density, which can effectively reduce the demand for input data angle. The results can provide a theoretical basis for the study of S-wave inversion methods.
郭强, 雒聪, 刘红达, 等.自适应优化参数模拟退火的叠前地震联合反演方法[J].石油地球物理勘探, 2023, 58(3):670-679.GUO Qiang, LUO Cong, LIU Hongda, et al.Prestack seismic hybrid inversion based on simulated annealing algorithm with adaptive optimization parameters[J].Oil Geophysical Prospecting, 2023, 58(3):670-679.
[2]
RUSSELL B H, HEDLIN K, HILTERMAN F J, et al.Fluid-property discrimination with AVO:A Biot-Gassmann perspective[J].Geophysics, 2003, 68(1):29-39.
[3]
RUSSELL B H, GRAY D, HAMPSON D P.Linearized AVO and poroelasticity[J].Geophysics, 2011, 76(3):C19-C29.
[4]
DE NICOLAO A, DRUFUCA G, ROCCA F.Eigenvalues and eigenvectors of linearized elastic inversion[J].Geophysics, 1993, 58(5):670-679.
[5]
WANG Y.Simultaneous inversion for model geometry and elastic parameters[J].Geophysics, 1999, 64(1):182-190.
[6]
BLOM N, BOEHM C, FICHTNER A.Synthetic inversions for density using seismic and gravity data[J].Geophysical Journal International, 2017, 209(2):1204-1220.
[7]
ROY B, ANNO P, GURCH M.Imaging oil-sand reservoir heterogeneities using wide-angle prestack seismic inversion[J].The Leading Edge, 2008, 27(9):1192-1201.
[8]
DOWNTON J E.Seismic Parameter Estimation from AVO Inversion[D].University of Calgary, Calgary, 2005.
[9]
曹俊兴, 薛雅娟, 田仁飞, 等.深层碳酸盐岩储层含气性检测方法技术研究[J].石油物探, 2019, 58(1):9-16.CAO Junxing, XUE Yajuan, TIAN Renfei, et al.Advances in hydrocarbon detection in deep carbonate reservoirs[J].Geophysical Prospecting for Petroleum, 2019, 58(1):9-16.
[10]
GRANLI J R, ARNTSEN B, SOLLID A, et al.Imaging through gas-filled sediments using marine shear-wave data[J].Geophysics, 1999, 64(3):668-677.
[11]
AVSETH P, MUKERJI T, MAVKO G.Quantitative Seismic Interpretation:Applying Rock Physics Tools to Reduce Interpretation Risk[M].Cambridge University Press, Cambridge, 2005.
[12]
CRAMPIN S.Evaluation of anisotropy by shear-wave splitting[J].Geophysics, 1985, 50(1):142-152.
[13]
DUFFAUT K, ALSOS T, ROGNØ H, et al.Shear-wave elastic impedance[J].The Leading Edge, 2000, 19(11):1222-1229.
[14]
GONG X, DU Q, ZHAO Q.SP-and SS-imaging for 3D elastic reverse time migration[J].Geophysics, 2018, 83(1):A1-A6.
[15]
KAZEI V, ALKHALIFAH T.Scattering radiation pattern atlas:What anisotropic elastic properties can body waves resolve?[J].Journal of Geophysical Research:Solid Earth, 2019, 124(3):2781-2811.
[16]
王赟, 孙丽霞, 李栋青, 等.勘探地震中的六分量观测[J].石油物探, 2021, 60(1):13-24.WANG Yun, SUN Lixia, LI Dongqing, et al.Six-component observation for exploration seismology[J].Geophysical Prospecting for Petroleum, 2021, 60(1):13-24, 33.
[17]
王波, 聂其海, 陈进娥, 等.四维多波地震在油藏动态监测中的应用[J].石油地球物理勘探, 2021, 56(2):340-345.WANG Bo, NIE Qihai, CHEN Jin'e, et al.Application of 4D multi-component seismic survey in dynamic reservoir monitoring[J].Oil Geophysical Prospecting, 2021, 56(2):340-345.
[18]
代福材, 张峰, 李向阳, 等.多波联合反演研究进展[J].石油物探, 2021, 60(1):57-69.DAI Fucai, ZHANG Feng, LI Xiangyang, et al.Pro-gress of multi-wave joint inversion[J].Geophysical Prospecting for Petroleum, 2021, 60(1):57-69.
[19]
李向阳, 张少华.勘探地震中横波分裂研究四十年回顾[J].石油物探, 2021, 60(2):190-209.LI Xiangyang, ZHANG Shaohua.Forty years of shear-wave splitting in seismic exploration:An overview[J].Geophysical Prospecting for Petroleum, 2021, 60(2):190-209.
[20]
BHATTACHARYA S N.Exact solutions of SH wave equation for inhomogeneous media[J].Bulletin of the Seismological Society of America, 1970, 60(6):1847-1859.
[21]
AKI K, RICHARDS P G.Quantitative Seismology:Theory and Methods[M].W.H.Freeman and Co, San Francisco, 1980.
[22]
VIRIEUX J.SH-wave propagation in heterogeneous media:Velocity-stress finite-difference method[J].Geophysics, 1984, 49(11):1933-1942.
[23]
姚振兴, 徐大方.关于一种新型横波震源(SH型)的探讨[J].地球物理学报, 1986, 29(4):382-389.YAO Zhenxing, XU Dafang.A pure shear wave source for reflection prospecting[J].Chinese Journal of Geophysics, 1986, 29(4):382-389.
[24]
FU L Y, BOUCHON M.Discrete wavenumber solutions to numerical wave propagation in piecewise heterogeneous media:I.Theory of two-dimensional SH case[J].Geophysical Journal International, 2004, 157(2):481-498.
[25]
王美霞, 杨顶辉, 宋国杰.二维SH波方程的半解析解及其数值模拟[J].地球物理学报, 2012, 55(3):914-924.WANG Meixia, YANG Dinghui, SONG Guojie.Semi-analytical solutions and numerical simulations of 2D SH wave equation[J].Chinese Journal of Geophysics, 2012, 55(3):914-924.
[26]
WON I J, CLOUGH J W.A new torsional shear-wave generator[J].Geophysics, 1981, 46(11):1607-1608.
[27]
詹正彬, 姚姚.多波及横波地震勘探[M].北京:地质出版社, 1994.ZHAN Zhengbin, YAO Yao.Multi Wave and Transverse Wave Seismic Prospecting[M].Geology Press, Beijing, 1994.
[28]
DASIOS A, MCCANN C, ASTIN T R, et al.Seismic imaging of the shallow subsurface:shear-wave case histories[J].Geophysical Prospecting, 1999, 47(4):565-591.
[29]
BUSH I, CRAMPIN S.Paris Basin VSPs:case history establishing combinations of fine-layer (or lithologic) anisotropy and crack anisotropy from modelling shear wavefields near point singularities[J].Geophysical Journal International, 1991, 107(3):433-447.
[30]
LI X, CRAMPIN S.Linear-transform techniques for processing shear-wave anisotropy in four-component seismic data[J].Geophysics, 1993, 58(2):240-256.
[31]
THOMSEN L.Converted-wave reflection seismo-logy over inhomogeneous, anisotropic media[J].Geophysics, 1999, 64(3):678-690.
[32]
LI X Y, YUAN J.Converted-wave moveout and conversion-point equations in layered VTI media:theory and applications[J].Journal of Applied Geophysics, 2003, 54(3/4):297-318.
[33]
李录明, 罗省贤.多波资料处理及解释方法的研究进展[J].石油地球物理勘探, 2006, 41(6):663-671.LI Luming, LUO Shengxian.Progress in study of multi-wave data processing and interpretation[J].Oil Geophysical Prospecting, 2006, 41(6):663-671.
[34]
ZHANG F, WANG Y.Wavelet-preserved PP-and PS-wave registration[J].Journal of Geophysics and Engineering, 2010, 7(4):395-403.
[35]
WANG Y, WANG C, LI J, et al.Frequency-dependent S-wave splitting parameters analysis:A case study from azimuthal PS data, Sanhu area of the Qaidam Basin, China[J].Geophysics, 2019, 84(6):B375-B386.
[36]
DAI H, LI X Y, CONWAY P.Imaging beneath gas clouds using 3D prestack Kirchhoff time migration of PS-converted waves:A case study from the North Sea[J].The Leading Edge, 2007, 26(4):522-529.
[37]
VEIRE H H, LANDRØ M.Simultaneous inversion of PP and PS seismic data[J].Geophysics, 2006, 71(3):R1-R10.
[38]
ZHANG F.Joint Inversion of Seismic PP-and PS-Waves in the Ray Parameter Domain[D].Imperial College London, London, 2010.
[39]
芦俊, 王赟, 石瑛.利用多波地震与测井数据联合反演预测煤岩的坚固性[J].地球物理学报, 2011, 54(11):2967-2972.LU Jun, WANG Yun, SHI Ying.Coal hardness prediction using joint inversion of multi-wave seismic data and logging[J].Chinese Journal of Geophysics, 2011, 54(11):2967-2972.
[40]
ALKAN E, HARDAGE B.Directs wave modes produced by vertical and horizontal vibrators[J].Journal of Seismic Exploration, 2013, 22(2):147-168.
[41]
GUPTA M, DEANGELO M V, HARDAGE B.PP and SS wave interpretation of a carbonate formation:A case study from the Arbuckle interval in Wellington field, Kansas[C].SEG Technical Program Expanded Abstracts, 2015, 34:2082-2087.
[42]
DENG Z, WU W, JING Y, et al.The joint PP and SH-SH data characterization of structures and reservoirs[C].Extended Abstracts of 81st EAGE Conference and Exhibition, 2019, 1-5.
[43]
ZHANG F, LI X.Inversion of the reflected SV-wave for density and S-wave velocity structures[J].Geophysical Journal International, 2019, 221(1):1635-1639.
[44]
王波, 张峰, 代福材, 等.VTI介质SH-SH波地震反演方法研究[J].地球物理学报, 2023, 66(5):2112-2122.WANG Bo, ZHANG Feng, DAI Fucai, et al.Study on seismic inversion method of SH-SH wave in VTI media[J].Chinese Journal of Geophysics, 2023, 66(5):2112-2122.
[45]
BORTFELD R.Approximations to the reflection and transmission coefficients of plane longitudinal and transverse waves[J].Geophysical Prospecting, 1961, 9(4):485-502.
[46]
RICHARDS P G, FRASIER C W.Scattering of elastic waves from depth-dependent in homogeneities[J].Geophysics, 1976, 41(3):441-458.
[47]
FATTI J L, SMITH G C, VAIL P J, et al.Detection of gas in sandstone reservoirs using AVO analysis:A 3-D seismic case history using the Geostack technique[J].Geophysics, 1994, 59(9):1362-1376.
[48]
ZOU K.S-Zero Stack:A converted wave processing to extract subsurface density information[J].Geophysics, 2014, 79(3):N1-N10.
[49]
周竹生.P-SV波和SH波的AVO分析[J].石油地球物理勘探, 1993, 28(4):430-438.ZHOU Zhusheng.AVO analysis of P-SV and SH waves[J].Oil Geophysical Prospecting, 1993, 28(4):430-438.
[50]
DAI F, ZHANG F, LI X.SH-SH wave inversion for S-wave velocity and density[J].Geophysics, 2022, 87(3):A25-A32.
[51]
ZHANG F.Simultaneous inversion for S-wave velocity and density from the SV-SV wave[J].Geophysics, 2021, 86(2):R187-R195.
[52]
ZHANG F, ZHANG T, LI X.Seismic amplitude inversion for the transversely isotropic media with vertical axis of symmetry[J].Geophysical Prospecting, 2019, 67(9):2368-2385.
[53]
WEI Z T, TANG X M.Numerical simulation of radiation, reflection, and reception of elastic waves from a borehole dipole source[J].Geophysics, 2012, 77(6):D253-D261.
[54]
JIN S, CAMBOIS G, VUILLERMOZ C.Shear-wave velocity and density estimation from PS-wave AVO analysis:Application to an OBS dataset from the North Sea[J].Geophysics, 2000, 65(5):1446-1454.