Laplace-domain waveform modeling and inversion for the 3D acoustic–elastic coupled media

Abstract Laplace-domain waveform inversion reconstructs long-wavelength subsurface models by using the zero-frequency component of damped seismic signals. Despite the computational advantages of Laplace-domain waveform inversion over conventional frequency-domain waveform inversion, an acoustic assumption and an iterative matrix solver have been used to invert 3D marine datasets to mitigate the intensive computing cost. In this study, we develop a Laplace-domain waveform modeling and inversion algorithm for 3D acoustic–elastic coupled media by using a parallel sparse direct solver library (MUltifrontal Massively Parallel Solver, MUMPS). We precisely simulate a real marine environment by coupling the 3D acoustic and elastic wave equations with the proper boundary condition at the fluid-solid interface. In addition, we can extract the elastic properties of the Earth below the sea bottom from the recorded acoustic pressure datasets. As a matrix solver, the parallel sparse direct solver is used to factorize the non-symmetric impedance matrix in a distributed memory architecture and rapidly solve the wave field for a number of shots by using the lower and upper matrix factors. Using both synthetic datasets and real datasets obtained by a 3D wide azimuth survey, the long-wavelength component of the P-wave and S-wave velocity models is reconstructed and the proposed modeling and inversion algorithm are verified. A cluster of 80 CPU cores is used for this study.

[1]  C. Shin,et al.  Two‐dimensional waveform inversion of multi‐component data in acoustic‐elastic coupled media , 2008 .

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

[3]  P. Mora Nonlinear two-dimensional elastic inversion of multioffset seismic data , 1987 .

[4]  Qing Huo Liu,et al.  Higher-Order Numerical Methods for Transient Wave Equations , 2003 .

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

[6]  J. Virieux P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[7]  D. Komatitsch,et al.  Wave propagation near a fluid-solid interface : A spectral-element approach , 2000 .

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

[9]  Changsoo Shin,et al.  Waveform inversion using a logarithmic wavefield , 2006 .

[10]  Changsoo Shin,et al.  Why do Laplace-domain waveform inversions yield long-wavelength results? , 2013 .

[11]  Gilles Lambaré,et al.  Velocity macro‐model estimation from seismic reflection data by stereotomography , 1998 .

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

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

[14]  Sean C. Solomon,et al.  Tomographic imaging of the shallow crustal structure of the East Pacific Rise at 9°30′N , 1994 .

[15]  K. Marfurt Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations , 1984 .

[16]  Changsoo Shin,et al.  3D acoustic waveform inversion in the Laplace domain using an iterative solver Sukjoon Pyun*, Inha University, Woohyun Son, and Changsoo Shin, Seoul National University , 2010 .

[17]  Paul Sava,et al.  Wave‐equation migration velocity analysis with time‐shift imaging , 2011 .

[18]  Patrick R. Amestoy,et al.  Multifrontal parallel distributed symmetric and unsymmetric solvers , 2000 .

[19]  R. Madariaga Dynamics of an expanding circular fault , 1976, Bulletin of the Seismological Society of America.

[20]  Changsoo Shin,et al.  Proof of the existence of both zero- and low-frequency information in a damped wavefield , 2012 .

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

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

[23]  Changsoo Shin,et al.  2D Elastic Waveform Inversion in the Laplace Domain , 2010 .

[24]  Gilles Lambaré,et al.  Practical aspects and applications of 2D stereotomography , 2003 .

[25]  C. Shin,et al.  Laplace-domain wave-equation modeling and full waveform inversion in 3D isotropic elastic media , 2014 .

[26]  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.

[27]  R. Shipp,et al.  Two-dimensional full wavefield inversion of wide-aperture marine seismic streamer data , 2002 .

[28]  A. Levander Fourth-order finite-difference P-SV seismograms , 1988 .

[29]  Changsoo Shin,et al.  Comparison of waveform inversion, part 1: conventional wavefield vs logarithmic wavefield , 2007 .

[30]  Changsoo Shin,et al.  2D acoustic‐elastic coupled waveform inversion in the Laplace domain , 2010 .

[31]  C. Shin,et al.  Waveform inversion in the Laplace domain , 2008 .

[32]  Source estimation and direct wave reconstruction in Laplace-domain waveform inversion for deep-sea seismic data , 2011 .

[33]  Oz Yilmaz,et al.  Migration velocity analysis by wave-field extrapolation , 1984 .

[34]  Full waveform inversion comparison of conventional and broadband marine seismic streamer data, NW Shelf Australia , 2015 .

[35]  Jean Virieux,et al.  Quantitative imaging of complex structures from dense wide‐aperture seismic data by multiscale traveltime and waveform inversions: a case study , 2004 .

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

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

[38]  O. C. Zienkiewicz,et al.  The Finite Element Method: Its Basis and Fundamentals , 2005 .

[39]  A. Tarantola A strategy for nonlinear elastic inversion of seismic reflection data , 1986 .