SeisElastic2D: An open-source package for multiparameter full-waveform inversion in isotropic-, anisotropic- and visco-elastic media

Abstract Full-waveform inversion (FWI) has emerged as a powerful technique to obtain high-resolution subsurface elastic properties. However, several complicating features, including interparameter trade-offs, cycle-skipping, and high computational cost, motivate a careful assessment and validation of candidate versions of FWI prior to use. Field data application of elastic FWI remains a challenging task, and widely available tools for inversion and appraisal do not exist. SeisElastic2D is an open-source package for multiparameter FWI in isotropic-, anisotropic- and visco-elastic media consisting of a set of FORTRAN 90 routines and SHELL scripts. This package has a modular structure built upon several existing open-source packages including upgraded SPECFEM2D, seisDD, and Seismic Unix. Various model parameterizations for isotropic- and anisotropic-elastic FWI, designed to reduce interparameter trade-offs, are available, as are different misfit functions designed to mitigate cycle-skipping and source-structure trade-offs. New misfit functions, measuring amplitude variation and central-frequency shift, are provided to estimate attenuation models in visco-elastic FWI. This package can be installed and used on high-performance computing cluster and local workstation supporting massively parallel interface, in order to satisfy the high computational requirements of elastic FWI. These features of SeisElastic2D make it a powerful tool for elastic FWI to bridge the gap between academic studies and industrial applications.

[1]  Jean-Yves L'Excellent,et al.  FWT2D: A massively parallel program for frequency-domain full-waveform tomography of wide-aperture seismic data - Part 2: Numerical examples and scalability analysis , 2009, Computational Geosciences.

[2]  I. Tsvankin.  P-wave signatures and notation for transversely isotropic media: An overview , 1996 .

[3]  Jeroen Tromp,et al.  Seismic waveform inversion best practices: regional, global, and exploration test cases , 2016 .

[4]  M. Warner,et al.  Semiglobal viscoacoustic full-waveform inversion , 2019, GEOPHYSICS.

[5]  Jeroen Tromp,et al.  Anelastic sensitivity kernels with parsimonious storage for adjoint tomography and full waveform inversion , 2016, 1604.05768.

[6]  Dmitry Borisov,et al.  Three-dimensional elastic full waveform inversion in a marine environment using multicomponent ocean-bottom cables: a synthetic study , 2015 .

[7]  B. Biondi,et al.  Simultaneous Inversion of Full Data Bandwidth by Tomographic Full Waveform Inversion , 2014 .

[8]  R. Plessix.  A review of the adjoint-state method for computing the gradient of a functional with geophysical applications , 2006 .

[9]  Gabriel Fabien-Ouellet,et al.  Time-domain seismic modeling in viscoelastic media for full waveform inversion on heterogeneous computing platforms with OpenCL , 2017, Comput. Geosci..

[10]  W. Pan,et al.  Amplitude-based misfit functions in viscoelastic full-waveform inversion applied to walk-away vertical seismic profile data , 2019, GEOPHYSICS.

[11]  J. Nocedal.  Updating Quasi-Newton Matrices With Limited Storage , 1980 .

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

[13]  D. Komatitsch,et al.  Introduction to the spectral element method for three-dimensional seismic wave propagation , 1999 .

[14]  Tariq Alkhalifah,et al.  A new full waveform inversion method based on shifted correlation of the envelope and its implementation based on OPENCL , 2019, Comput. Geosci..

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

[16]  Xiao‐Bi Xie,et al.  Full-intensity waveform inversion , 2018, GEOPHYSICS.

[17]  Ludovic Métivier,et al.  A guided tour of multiparameter full-waveform inversion with multicomponent data: From theory to practice , 2013 .

[18]  Wim A. Mulder,et al.  A correlation-based misfit criterion for wave-equation traveltime tomography , 2010 .

[19]  Qinya Liu,et al.  Tomography, Adjoint Methods, Time-Reversal, and Banana-Doughnut Kernels , 2004 .

[20]  Joakim O. Blanch,et al.  Modeling of a constant Q; methodology and algorithm for an efficient and optimally inexpensive viscoelastic technique , 1995 .

[21]  Brian Russell,et al.  Fluid-property discrimination with AVO: A Biot-Gassmann perspective , 2003 .

[22]  Ludovic Métivier,et al.  The SEISCOPE optimization toolbox: A large-scale nonlinear optimization library based on reverse communication , 2016 .

[23]  Robin P. Fletcher,et al.  Reverse time migration in tilted transversely isotropic "TTI… media , 2009 .

[24]  Marwan Charara,et al.  Feasibility study for an anisotropic full waveform inversion of cross‐well seismic data , 2008 .

[25]  Jianlei Zhang,et al.  Elastic full-waveform inversion based on GPU accelerated temporal fourth-order finite-difference approximation , 2020, Comput. Geosci..

[26]  J. Castagna,et al.  Relationships between compressional‐wave and shear‐wave velocities in clastic silicate rocks , 1985 .

[27]  J. Tromp,et al.  Finite-Frequency Kernels Based on Adjoint Methods , 2006 .

[28]  Liangguo Dong,et al.  Correlation-based reflection full-waveform inversion , 2015 .

[29]  Qingchen Zhang,et al.  Crosstalk-free simultaneous-source full waveform inversion with normalized seismic data , 2020, Comput. Geosci..

[30]  Felix J. Herrmann,et al.  Mitigating local minima in full-waveform inversion by expanding the search space , 2013 .

[31]  Separation of Migration and Tomography Modes of Full‐Waveform Inversion in the Plane Wave Domain , 2018 .

[32]  Gerard T. Schuster,et al.  Q-Least Squares Reverse Time Migration with Viscoacoustic Deblurring Filters , 2017 .

[33]  Use of prismatic waves in full-waveform inversion with the exact Hessian , 2020 .

[34]  Kristopher A. Innanen,et al.  Seismic AVO and the inverse Hessian in precritical reflection full waveform inversion , 2014 .

[35]  J. Korringa,et al.  On the Dependence of the Elastic Properties of a Porous Rock on the Compressibility of the Pore Fluid , 1975 .

[36]  K. Innanen,et al.  Estimating P- and S-wave inverse quality factors from observed seismic data using an attenuative elastic impedance , 2018 .

[37]  C. Shin,et al.  Improved amplitude preservation for prestack depth migration by inverse scattering theory , 2001 .

[38]  Dong-Joo Min,et al.  Frequency-domain elastic full waveform inversion for VTI media , 2010 .

[39]  J. Tromp,et al.  Misfit functions for full waveform inversion based on instantaneous phase and envelope measurements , 2011 .

[40]  D. Köhn,et al.  Time Domain 2D Elastic Full Waveform Tomography , 2011 .

[41]  José M. Carcione,et al.  Viscoacoustic wave propagation simulation in the earth , 1988 .

[42]  Romain Brossier,et al.  On the footprint of anisotropy on isotropic full waveform inversion: the Valhall case study , 2011 .

[43]  R. G. Pratt,et al.  Inversion strategies for visco-acoustic waveform inversion , 2013 .

[44]  J. Maddison,et al.  Full wavefield inversion of ambient seismic noise , 2018 .

[45]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[46]  Kristopher A. Innanen.  Inversion of the seismic AVF/AVA signatures of highly attenuative targets , 2011 .

[47]  Wei Zhou,et al.  A Time-Domain Preconditioned Truncated Newton Approach to Visco-acoustic Multiparameter Full Waveform Inversion , 2018, SIAM J. Sci. Comput..

[48]  S. Operto,et al.  Fast full waveform inversion with source encoding and second-order optimization methods , 2015 .

[49]  Jeroen Tromp,et al.  Double-difference adjoint seismic tomography , 2016, 1607.01811.

[50]  D. Nichols,et al.  Integrated VTI model building with seismic data, geologic information, and rock-physics modeling — Part 1: Theory and synthetic test , 2016 .

[51]  C. Bunks,et al.  Multiscale seismic waveform inversion , 1995 .

[52]  Laurent Demanet,et al.  Full-waveform inversion with extrapolated low-frequency data , 2016 .

[53]  Denes Vigh,et al.  3D prestack plane-wave, full-waveform inversion , 2008 .

[54]  W. Pan,et al.  Interparameter trade-off quantification for isotropic-elastic full-waveform inversion with various model parameterizations , 2019, GEOPHYSICS.

[55]  Wenyuan Liao,et al.  Total variation regularization for seismic waveform inversion using an adaptive primal dual hybrid gradient method , 2018 .

[56]  Research Note: Insights into the data dependency on anisotropy: an inversion prospective , 2016 .

[57]  Joakim O. Blanch,et al.  Viscoelastic finite-difference modeling , 1994 .

[58]  L. Thomsen.  Weak elastic anisotropy , 1986 .

[59]  Mauricio D. Sacchi,et al.  Elastic least-squares reverse time migration via linearized elastic full-waveform inversion with pseudo-Hessian preconditioning , 2017 .

[60]  Youli Quan,et al.  Seismic attenuation tomography using the frequency shift method , 1997 .

[61]  B. Biondi,et al.  Target-oriented elastic full-waveform inversion through extended-migration redatuming , 2018, SEG Technical Program Expanded Abstracts 2018.

[62]  R. Tonn,et al.  THE DETERMINATION OF THE SEISMIC QUALITY FACTOR Q FROM VSP DATA: A COMPARISON OF DIFFERENT COMPUTATIONAL METHODS1 , 1991 .

[63]  Yanhua O. Yuan,et al.  Multiscale adjoint waveform tomography for surface and body waves , 2015 .

[64]  Gerard T. Schuster,et al.  Wave-equation traveltime inversion , 1991 .

[65]  Frederik J. Simons,et al.  Multiscale adjoint waveform-difference tomography using wavelets , 2014 .

[66]  George A. McMechan,et al.  Comparison of two viscoacoustic propagators for Q-compensated reverse time migration , 2016 .

[67]  Albert Tarantola,et al.  Theoretical background for the inversion of seismic waveforms including elasticity and attenuation , 1988 .

[68]  Jeroen Tromp,et al.  SeisFlows - Flexible waveform inversion software , 2018, Comput. Geosci..

[69]  W. Pan,et al.  Interparameter trade-off quantification and reduction in isotropic-elastic full-waveform inversion: synthetic experiments and Hussar land data set application , 2018 .

[70]  René-Édouard Plessix,et al.  A parametrization study for surface seismic full waveform inversion in an acoustic vertical transversely isotropic medium , 2011 .

[71]  F. Dahlen,et al.  Surface-wave group-delay and attenuation kernels , 2006 .

[72]  Hao Hu,et al.  Prestack correlative least-squares reverse time migration , 2017 .

[73]  Yanfei Wang,et al.  Parameterization analysis and field validation of VTI-elastic full-waveform inversion in a walk-away vertical seismic profile configuration , 2020 .

[74]  Romain Brossier,et al.  Two-dimensional frequency-domain visco-elastic full waveform inversion: Parallel algorithms, optimization and performance , 2011, Comput. Geosci..

[75]  David C. Henley,et al.  Gabor deconvolution: Estimating reflectivity by nonstationary deconvolution of seismic data , 2011 .

[76]  Lianjie Huang,et al.  Acoustic- and elastic-waveform inversion with total generalized p-variation regularization , 2018, Geophysical Journal International.

[77]  Ru-Shan Wu,et al.  Seismic envelope inversion and modulation signal model , 2014 .

[78]  Yanfei Wang,et al.  On the influence of different misfit functions for attenuation estimation in viscoelastic full-waveform inversion: synthetic study , 2020 .

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

[80]  Don L. Anderson,et al.  Velocity dispersion due to anelasticity; implications for seismology and mantle composition , 1976 .

[81]  W. Pan,et al.  Elastic full-waveform inversion and parametrization analysis applied to walk-away vertical seismic profile data for unconventional (heavy oil) reservoir characterization , 2018 .

[82]  J. Tromp,et al.  3D elastic full-waveform inversion of surface waves in the presence of irregular topography using an envelope-based misfit function , 2018 .

[83]  B. Biondi,et al.  Preconditioned elastic full-waveform inversion with approximated Hessian , 2017 .