Geological Structure-Guided Initial Model Building for Prestack AVO/AVA Inversion

Reconstructing an accurate and high-resolution subsurface model is attractive in the fields of both geology and seismology. However, due to the band-limited characteristics of seismic data, the inversion greatly depends on the reliability of the initial model. A fairly acceptable initial model could lay a good foundation for seismic inversion. In this article, we first introduce a well-log interpolation method with the local slope as a constraint for building a high-fidelity starting model in prestack amplitude versus offset/angle (AVO/AVA) inversion. First, we briefly review the basic theory of general seismic inversion. Then, instead of using the conventional preconditioned least-squares method, we introduce shaping regularization theory into the geological structure-guided well-log interpolation to accelerate the convergence. We use the plane-wave destruction (PWD) algorithm to extract the slope attribute from seismic data, images, or velocity models. The slope is used as the constraint to solve the inverse problem based on the shaping regularization method. Numerical examples demonstrate that the proposed initial model building method performs better than the conventional ones. It greatly improves the accuracy of inversion results. Furthermore, we apply the proposed model building method to the inverse problems of AVO/AVA inversion and reservoir parameter estimation of several field data sets for the first time, which demonstrate encouraging performance.

[1]  Ruben D. Martinez,et al.  Practical approaches for subsalt velocity model building , 2008 .

[2]  Xianhuai Zhu,et al.  Tomostatics: Turning‐ray tomography + static corrections , 1992 .

[3]  Konstantin Osypov,et al.  Robust refraction tomography , 2000 .

[4]  Colin A. Zelt,et al.  A Case History: Application of Frequency-Dependent Traveltime Tomography and Full Waveform Inversion to a Known Near-Surface Target , 2013 .

[5]  M. Warner,et al.  Anisotropic 3D full-waveform inversion , 2013 .

[6]  Dave Hale,et al.  Velocity analysis using weighted semblance , 2012 .

[7]  Qiang Guo,et al.  Prestack Seismic Inversion Based on Anisotropic Markov Random Field , 2018, IEEE Transactions on Geoscience and Remote Sensing.

[8]  Jean Virieux,et al.  Velocity model building by 3D frequency-domain, full-waveform inversion of wide-aperture seismic data , 2008 .

[9]  Siwei Li Wave-Equation Migration Velocity Analysis by Non-Stationary Focusing , 2013 .

[10]  Yangkang Chen,et al.  Seismic imaging of incomplete data and simultaneous-source data using least-squares reverse time migration with shaping regularization , 2016 .

[11]  Bin Wang,et al.  Advances in velocity model-building technology for subsalt imaging , 2008 .

[12]  Paul Sava,et al.  Wave‐equation migration velocity analysis. II. Subsalt imaging examples , 2004 .

[13]  Hicks,et al.  Gauss–Newton and full Newton methods in frequency–space seismic waveform inversion , 1998 .

[14]  Shuangquan Chen,et al.  Amplitude variation with angle inversion using the exact Zoeppritz equations — Theory and methodology , 2016 .

[15]  Xiaohong Chen,et al.  Prestack Waveform Inversion by Using an Optimized Linear Inversion Scheme , 2019, IEEE Transactions on Geoscience and Remote Sensing.

[16]  Min Bai,et al.  Obtaining free USArray data by multi-dimensional seismic reconstruction , 2019, Nature Communications.

[17]  M. Turhan Taner,et al.  Velocity spectra-digital computer derivation and applications of velocity functions , 1969 .

[18]  William W. Symes,et al.  Migration velocity analysis and waveform inversion , 2008 .

[19]  Florian A. Belina,et al.  Waveform Inversion of Crosshole Georadar Data: Influence of Source Wavelet Variability and the Suitability of a Single Wavelet Assumption , 2012, IEEE Transactions on Geoscience and Remote Sensing.

[20]  H. C. Larsen,et al.  Crustal structure of the southeast Greenland margin from joint refraction and reflection seismic tomography , 2000 .

[21]  Mike Warner,et al.  Next-generation seismic experiments: wide-angle, multi-azimuth, three-dimensional, full-waveform inversion , 2013 .

[22]  Jean Virieux,et al.  An overview of full-waveform inversion in exploration geophysics , 2009 .

[23]  Siyuan Cao,et al.  One-Step Slope Estimation for Dealiased Seismic Data Reconstruction via Iterative Seislet Thresholding , 2016, IEEE Geoscience and Remote Sensing Letters.

[24]  Marnik Vanclooster,et al.  Soil Surface Water Content Estimation by Full-Waveform GPR Signal Inversion in the Presence of Thin Layers , 2010, IEEE Transactions on Geoscience and Remote Sensing.

[25]  Raymond D. Mindlin,et al.  Compliance of elastic bodies in contact , 1949 .

[26]  Henri Calandra,et al.  High-performance 3D first-arrival traveltime tomography , 2010 .

[27]  Paul Sava,et al.  Wave-equation migration velocity analysis. I. Theory , 2004 .

[28]  Sergey Fomel,et al.  Velocity analysis using AB semblance , 2009 .

[29]  S. Qu,et al.  Velocity analysis of simultaneous-source data using high-resolution semblance—coping with the strong noise , 2016 .

[30]  Yangkang Chen,et al.  Geological structure guided well log interpolation for high-fidelity full waveform inversion , 2016 .

[31]  Børge Arntsen,et al.  Initial Velocity Models for Full Waveform Inversion , 2012 .

[32]  John R. Miller,et al.  Scaling-up and model inversion methods with narrowband optical indices for chlorophyll content estimation in closed forest canopies with hyperspectral data , 2001, IEEE Trans. Geosci. Remote. Sens..

[33]  Qing Huo Liu,et al.  Joint Inversion of Electromagnetic and Seismic Data Based on Structural Constraints Using Variational Born Iteration Method , 2018, IEEE Transactions on Geoscience and Remote Sensing.

[34]  Irena Hajnsek,et al.  Inversion of surface parameters from polarimetric SAR data , 2000, IGARSS 2000. IEEE 2000 International Geoscience and Remote Sensing Symposium. Taking the Pulse of the Planet: The Role of Remote Sensing in Managing the Environment. Proceedings (Cat. No.00CH37120).

[35]  A. Tarantola Inversion of seismic reflection data in the acoustic approximation , 1984 .

[36]  Stefano Tinti,et al.  Tsunami waveform inversion by numerical finite-elements Green’s functions , 2001 .

[37]  Chao Wang,et al.  High-Dimensional Waveform Inversion With Cooperative Coevolutionary Differential Evolution Algorithm , 2012, IEEE Geoscience and Remote Sensing Letters.

[38]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[39]  Richard Bamler,et al.  Tomographic SAR Inversion by $L_{1}$ -Norm Regularization—The Compressive Sensing Approach , 2010, IEEE Transactions on Geoscience and Remote Sensing.

[40]  Antoine Guitton,et al.  Constrained full-waveform inversion by model reparameterization1 , 2012 .

[41]  Yangkang Chen,et al.  Velocity analysis using similarity-weighted semblance , 2014 .

[42]  Hanming Chen,et al.  Modelling Elastic Wave Propagation Using A New Wave Equation and Temporal Fourth-order Finite-difference Method , 2016 .

[43]  Jinghuai Gao,et al.  A New Highly Efficient Differential Evolution Scheme and Its Application to Waveform Inversion , 2014, IEEE Geoscience and Remote Sensing Letters.

[44]  Jincai Zhang,et al.  Pore pressure prediction from well logs: Methods, modifications, and new approaches , 2011 .

[45]  Colin A. Zelt,et al.  Three‐dimensional seismic refraction tomography: A comparison of two methods applied to data from the Faeroe Basin , 1998 .

[46]  Wei Xiao,et al.  Hybrid Seismic Inversion Based on Multi-Order Anisotropic Markov Random Field , 2020, IEEE Transactions on Geoscience and Remote Sensing.

[47]  Thorkild B. Hansen,et al.  Inversion scheme for ground penetrating radar that takes into account the planar air-soil interface , 2000, IEEE Trans. Geosci. Remote. Sens..

[48]  R. Oezsen Velocity modelling and prestack depth imaging below complex salt structures: a case history from on-shore Germany , 2004 .

[49]  Sergey Fomel,et al.  Applications of plane-wave destruction filters , 2002 .

[50]  Jacques R. Ernst,et al.  A New Vector Waveform Inversion Algorithm for Simultaneous Updating of Conductivity and Permittivity Parameters From Combination Crosshole/Borehole-to-Surface GPR Data , 2010, IEEE Transactions on Geoscience and Remote Sensing.

[51]  Xiaohong Chen,et al.  Application of an adaptive acquisition regularization parameter based on an improved GCV criterion in pre-stack AVO inversion , 2017 .

[52]  Shaohuan Zu,et al.  Interpolating Big Gaps Using Inversion With Slope Constraint , 2016, IEEE Geoscience and Remote Sensing Letters.

[53]  Wei Chen,et al.  An open-source Matlab code package for improved rank-reduction 3D seismic data denoising and reconstruction , 2016, Comput. Geosci..

[54]  Yangkang Chen,et al.  A High-resolution Weighted Semblance for Dealing with AVO Phenomenon , 2016 .

[55]  R. Pratt Seismic waveform inversion in the frequency domain; Part 1, Theory and verification in a physical scale model , 1999 .

[56]  Jean Virieux,et al.  Building starting models for full waveform inversion from wide‐aperture data by stereotomography , 2013 .

[57]  Cong Luo,et al.  Pre-stack AVA Inversion by Using Propagator Matrix Forward Modeling , 2019, Pure and Applied Geophysics.

[58]  Harald Johnsen,et al.  SAR-ocean wave inversion using image cross spectra , 1995, IEEE Trans. Geosci. Remote. Sens..

[59]  William W. Symes,et al.  Automatic velocity analysis via shot profile migration , 2008 .

[60]  Xiaohong Chen,et al.  Application of Optimal Transport to Exact Zoeppritz Equation AVA Inversion , 2018, IEEE Geoscience and Remote Sensing Letters.

[61]  Yangkang Chen,et al.  Dealiased Seismic Data Interpolation Using Seislet Transform With Low-Frequency Constraint , 2015, IEEE Geoscience and Remote Sensing Letters.

[62]  Junjie Wu,et al.  Adaptive Differential Evolution by Adjusting Subcomponent Crossover Rate for High-Dimensional Waveform Inversion , 2015, IEEE Geoscience and Remote Sensing Letters.

[63]  Shaohuan Zu,et al.  Seismic imaging of simultaneous-source data using constrained least-squares reverse time migration , 2015 .