Nonlinear multi-parameter hybrid inversion of pre-stack and post-stack seismic data based on Zoeppritz equation
ZHANG Lingyuan1, ZHANG Hongbing1, SHANG Zuoping2, YAN Lizhi1, REN Quan1
1. College of Earth Science and Engineering, Hohai University, Nanjing, Jiangsu 211100, China; 2. College of Mechanics and Materials, Hohai University, Nanjing, Jiangsu 211100, China
Abstract:At present, Shuey or Gray approximate formula and basis pursuit theory are widely applied in inversion of Poisson ratio. The same is true for the velocity ratio which is calculated from P- and S- waves velocity obtained from pre-stack inversion. But it is not avoidable to reduce the accuracy of inversion parameters. We try to apply Zoeppritz equation in inversion for velocity ratio and Poisson ratio. First we developed a hybrid pre-stack and post-stack seismic inversion workflow by the convolution model based on the reflection coefficient of vertical incidence and Zoeppritz equation. Then we constructed a new objective function which includes edge-preserving regularization to reduce the adverse effect of inversion unsuitability. We applied a fast simulated annealing algorithm to achieve global non-linear optimization. To reduce the influence from the magnitude difference of multiple parameters and the poor stability of multi-parameter simultaneous random search, we used an improved Fatti equation in inversion, which introduces two fitting formulas, a density and P-wave velocity formula, and a S-velocity and P-wave velocity formula, to improve the accuracy and stability of inversion. The inverted results of field data indicate that the velocity ratio obtained by direct Zoeppritz equation inversion is better than that obtained by approximate formulas. The inversion results from the new workflow have good consistency with the logging data. In a word, post-stack inversion obtained more accurate P-wave velocity and density, while the accuracy of velocity ration obtained by pre-stack inversion is higher than that obtained by post-stack inversion, and the inversion of 18°~24° angular trace gathers is better than that of 3°~9° and 33°~39° gathers.
王秀娟,吴时国,余鹏,等.混合地震反演技术及其在水合物勘探中的应用[J].石油物探,2007,46(3):278-282.WANG Xiujuan,WU Shiguo,YU Peng,et al. Application of hybrid seismic inversion in hydrate exploration[J].Geophysical Prospecting for Petroleum, 2007,46(3):278-282.
[2]
李磊,王英民,杨绍国,等.混合地震反演技术及其在东海南部陆架盆地中的应用[J].中国石油大学学报,2007,31(5):28-34.LI Lei,WANG Yingmin,YANG Shaoguo,et al.Application of hybrid seismic inversion:a case study from the southern shelf basin of east China Sea[J].Journal of China University of Petroleum,2007,31(5):28-34.
[3]
杨文采.地球物理反演的理论与方法[M].北京:地质出版社,1997.
[4]
Pratt R G,Shin C,Hicks G J.Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion[J].Geophysical Journal International,1998, 133(2):341-362.
[5]
Virieux J,Operto S.An overview of full-waveform inversion in exploration geophysics[J]. Geophysics,2009,74(6):WCC1-WCC26.
[6]
张衡,刘洪,刘璐,等.基于平均导数方法的声波方程频率域高阶正演[J].地球物理学报,2014,57(5):1599-1611.ZHANG Heng,LIU Hong,LIU Lu,et al.Frequency domain acoustic equation high-order modeling based on an average-derivative method[J].Chinese Journal of Geophysics,2014,57(5):1599-1611.
[7]
陈生昌,陈国新.多主频波场的时间阻尼全波形反演[J].地球物理学报,2017, 60(8):3229-3237.CHEN Shengchang,CHEN Guoxin.Time-damping full waveform inversion of multi-dominant-frequency wavefields[J].Chinese Journal of Geophysics,2017,60(8):3229-3237.
[8]
Aki K,Richards P G.Quantitative Seismology:Theory and Methods[M].W.H. Freeman and Co.,1980.
[9]
Shuey R T.A simplification of the Zoeppritz equations[J]. Geophysics,1985,50(4):609-614.
[10]
Fatti J L,Smith G C,Vail P J.Detection of gas in sandstone reservoirs using AVO analysis:A 3-D seismic case history using the Geostack technique[J].Geo-physics,1994,59(9):1362-1376.
[11]
刘畅,李超,朱振宇,等.杨氏模量和泊松比基追踪反演[J]. 石油地球物理勘探,2019,54(6):1310-1315.LIU Chang,LI Chao,ZHU Zhenyu,et al.Basis pursuit inversion for Young's modulus and Poisson's ratio[J]. Oil Geophysical Prospecting, 2019,54(6):1310-1315.
[12]
张宏兵,尚作萍,杨长春,等.波阻抗反演正则参数估计[J]. 地球物理学报,2005, 48(1):181-188.ZHANG Hongbing,SHANG Zuoping,YANG Changchun,et al.Estimation of regular parameters for the impedance inversion[J].Chinese Journal of Geophy-sics, 2005,48(1):181-188.
[13]
吉洪诺夫,阿尔先宁(著);王秉忱(译).不适定问题的解法[M].北京:地质出版社,1979.
[14]
Sen M K,Roy I G.Computation of differential seismograms and iteration adaptive regularization in prestack waveform inversion[J].Geophysics,2003,68(6):2026-2039.
Zhang H,Guo Q,Liang L,et al.A nonlinear method for multiparameter inversion of pre-stack seismic data based on anisotropic Markov random field[J].Geophysical Prospecting,2018,66(3):461-477.
[19]
Liang L F,Zhang H B,Dan Z W,et al.Prestack density inversion using the Fatti equation constrained by the P- and S-wave impedance and density[J].Applied Geophysics,2017,14(1):133-141.
[20]
王保丽,印兴耀,张繁昌,等.基于Fatti近似的弹性阻抗方程及反演[J].地球物理学进展,2008,23(1):192-197.WANG Baoli,YIN Xingyao,ZHANG Fanchang,et al. Elastic impedance equation based on Fatti approximation and inversion[J].Progress in Geophysics,2008,23(1):192-197.
[21]
Tarantola A,Valette B.Generalized nonlinear inverse problems solved using the least squares criterion[J].Review of Geophysics and Space Physics,1982,20(2):219-232.
[22]
张宏彬,何礁登.宽带约束反演[J].石油物探,1995,34(1):1-10.ZHANG Hongbin,HE Qiaodeng.Broad-band constrained inversion[J].Geophysical Prospecting for Petroleum,1995, 34(1):1-10.
[23]
周竹生,何继善,赵荷晴.利用广义共轭梯度算法求解地震道反演问题[J].石油地球物理勘探,1998,33(4):439-447.ZHOU Zhusheng,HE Jishan,ZHAO Heqing.Solving a seismic trace inversion problem by using generalized conjugate gradient algorithm[J].Oil Geophysical Prospecting,1998,33(4):439-447.
[24]
刘璐,刘洪,张衡,等.基于修正拟牛顿公式的全波形反演[J].地球物理学报, 2013, 56(7):2447-2451.LIU Lu,LIU Hong,ZHANG Heng,et al.Full waveform inversion based on modified quasi-Newton equation[J]. Chinese Journal of Geophysics, 2013,56(7):2447-2451.
[25]
张霖斌,姚振兴,纪晨,等.快速模拟退火算法及应用[J]. 石油地球物理勘探,1997,32(5):645-660.ZHANG Linbin,YAO Zhenxing,JI Chen,et al.Fast simulated annealing algorithm and its application[J].Oil Geophysical Prospecting,1997,32(5):654-660.
[26]
Guo Q,Zhang H B,Tian J B,et al.A nonlinear multiparameter prestack seismic inversion method based on hybrid optimization approach[J].Arabian Journal of Geosciences,2018,11(3):48-60.
[27]
方中于,王丽萍,杜家元,等.基于混合智能优化算法的非线性AVO反演[J].石油地球物理勘探,2017,52(4):797-804.FANG Zhongyu,WANG Liping,DU Jiayuan,et al. Nonlinear AVO inversion based on hybrid intelligent optimization algorithm[J].Oil Geophysical Prospecting,2017,52(4):797-804.
[28]
Charbonnier P,Blanc-Féraud L,Aubert G.Determi-nistic edge-preserving regularization in computed imaging[J]. IEEE Transactions on Image Processing,1997,6(2):298-311.
[29]
郭强,张宏兵,曹呈浩,等.Zoeppritz方程叠前多参数反演及密度敏感性分析[J].石油地球物理勘探,2017,52(6):1218-1225.GUO Qiang,ZHANG Hongbing,CAO Chenghao,et al. Density-sensitivity analysis about prestack multi-parameter inversion based on the exact Zoeppritz equation[J].Oil Geophysical Prospecting, 2017,52(6):1218-1225.