Generalized impedance blocky inversion based on a-nalytic solution to wave equation
LI Yuanqiang1,2,3, HUO Zhizhou4, LI Jingye1,2,3, CHEN Xiaohong1,2,3, ZHANG Jian1,2,3, GENG Weiheng1,2,3
1. College of Geophysics, China University of Petroleum(Beijing), Beijing 102249, China; 2. State Key Laboratory of Petroleum Resources and Prospecting, Beijing 102249, China; 3. National Engineering Laboratory for Offshore Oil Exploration, Beijing 102249, China; 4. Sinopec Petroleum Exploration and Production Research Institute, Beijing 100083, China
Abstract:Since the pre-processing of pre-stack gathers is based on the assumption of acoustic medain in many cases,the the gathers tend to be with more acoustic AVO features.In addition,density inversion is unstable.This paper proposes a generalized impedance blocky inversion based on analytic solution to acoustic wave equation.The generalized acoustic impedance is inverted by a partially stacked profile,which varies with the angle of incidence; and on this basis,more accurate velocity and stable density are extracted.For the conventional impe-dance inversion method,transmission loss and inter-layer multiples are neglected.Based on the recursive formula of derivation,the one-dimensional acoustic wave equation is solved analytically to obtain the full-wavefield responses at different incident angles,and the Fréchet derivatives are analytically derived for gradient-descent inversion algorithm.Most of the inversion methods are based on smoothing constraints,which fundamentally lead to unfocused boundaries for inversion results.In order to improve the resolution of the inversion results,blocky constraints can be introduced based on the Bayesian inference framework to obtain stable and high resolution inversion results.According to the above theory,we first uses model data to analyze the influence of the incompleteness of the forward method on seismic responses,further verify the validity of the inversion method,and extract the accurate velocity and density.Then the ability to cha-racterize the boundary for blocky constraintis tested by adding noises.Both model and actual data prove that the inversion results from the new method have higher resolution,the boundary is clearer,and the extracted velocity and density profiles are stable and accurate.
马劲风.地震勘探中广义弹性阻抗的正反演[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.
[3]
Bosch M,Mukerji T,Gonzalez E F.Seismic inversion for reservoir properties combining statistical rock physics and geostatistics:A review[J].Geophysics,2010,75(5):A165-A176.
[4]
罗鑫,陈学华,张杰,等.基于依赖频率AVO反演的高含气饱和度储层预测方法[J].石油地球物理勘探,2019,54(2):356-364.LUO Xin,CHEN Xuehua,ZHANG Jie,et al.High gas-saturation reservoir prediction based on frequency-dependent AVO inversion[J].Oil Geophysical Prospecting,2019,54(2):356-364.
[5]
Lindseth R O.Synthetic sonic logs:A process for stratigraphic interpretation[J].Geophysics,1979,44(1):3-26.
[6]
Walker C,Ulrych T J.Autoregressive recovery of the acoustic impedance[J].Geophysics,1983,48(10):1338-1350.
[7]
Ma X Q.A constrained global inversion method using an over-parameterized scheme:Application to post-stack seismic data[J].Geophysics,2001,66(2):613-626.
[8]
吕铁良.波阻抗约束反演中的约束方法研究[D].山东青岛:中国石油大学(华东),2007.
[9]
Cooke D A,Schneider W A.Generalized linear inversion of reflection seismic data[J].Geophysics,1983,48(6):665-676.
[10]
Oldenburg D W,Scheuer T,Levy S.Recovery of the acoustic impedance from reflection seismograms[J].Geophysics,1983,48(10):1318-1337.
Oliveira S A M,Braga I L S,Lacerda M B,et al.Extending the useful angle range for elastic inversion through the amplitude-versus-angle full-waveform inversion method[J].Geophysics,2018,83(3):R213-R226.
[13]
杨文采.非线性地震反演方法的补充及比较[J].石油物探,1995,34(4):109-116.YANG Wencai.Supplement and comparison of nonlinear seismic inversion methods[J].Geophysical Prospecting for Petroleum,1995,34(4):109-116.
[14]
Gholami A.Nonlinear multichannel impedance inversion by total-variation regularization[J].Geophysics,2015,80(5):R217-R224.
[15]
王治强,曹思远,陈红灵,等.基于TV约束和Zoeppritz矩阵分解的波阻抗反演[J].石油地球物理勘探,2017,52(6):1193-1199.WANG Zhiqiang,CAO Siyuan,CHEN Hongling,et al.Wave impedance inversion based on TV regularization and Zoeppritz-sparse matrix factorization[J].Oil Geophysical Prospecting,2017,52(6):1193-1199.
[16]
林恬,孟小红,张致付.基于约束最小二乘与信赖域的储层参数反演方法[J].地球物理学报,2017,60(10):3969-3983.LIN Tian,MENG Xiaohong,ZHANG Zhifu.The petrophysical parameter inversion method based on constrained least squares and trust region domain[J].Chinese Journal of Geophysics,2017,60(10):3969-3983.
[17]
Tarantola A.Inversion of seismic reflection data in the acoustic approximation[J].Geophysics,1984,49(8):1259-1266.
[18]
Virieux J,Operto S.An overview of full-waveform inversion in exploration geophysics[J].Geophysics,2009,74(6):WCC1-WCC26.
[19]
杨午阳,王西文,雍学善,等.地震全波形反演方法研究综述[J].地球物理学进展,2013,28(2):766-776.YANG Wuyang,WANG Xiwen,YONG Xueshan,et al.The review of seismic full waveform inversion method[J].Progress in Geophysics,2013,28(2):766-776.
[20]
McAulay A D.Prestack inversion with plane-layer point source modeling[J].Geophysics,1985,50(1):77-89.
[21]
Oliveira S,Loures L,Moraes F,et al.Nonlinear impedance inversion for attenuating media[J].Geophy-sics,2009,74(6):R111-R117.
[22]
Ursin B,Stovas A.Reflection and transmission responses of a layered isotropic viscoelastic medium[J].Geophysics,2002,67(1):307-323.
[23]
Alemie W,Sacchi M D.High-resolution three-term AVO inversion by means of a Trivariate Cauchy probability distribution[J].Geophysics,2011,76(3):R43-R55.
[24]
Theune U,Jensås I Ø,Eidsvik J.Analysis of prior models for a blocky inversion of seismic AVA data[J].Geophysics,2010,75(3):C25-C35.
[25]
Charbonnier P,Blanc-Féraud L,Aubert G,et al.Deterministic edge-preserving regularization in computed imaging[J].IEEE Transactions on Image Processing,1997,6(2):298-311.