3D acoustic modelling and waveform inversion in the Laplace domain for an irregular sea floor using the Gaussian quadrature integration method

Abstract This paper explores the concept of 3D acoustic modelling and inversion in the Laplace domain. We used the Gaussian quadrature integration method to assemble the element mass matrix of a coarse regular rectangular element containing both water and subsurface parts, approximating the properties of these two parts without refining the element into smaller components. The purpose of this investigation was to model an irregular sea floor using 3D Laplace-domain modelling and inversion without increasing the computational cost incurred when using a small grid interval to precisely delineate the sea floor. The algorithm was verified by comparing the analytical solution to the modelled wavefield obtained using the conventional Laplace-domain modelling method and the modelled wavefield obtained using the Gaussian quadrature integration in a two-layer model (with 50-m, 100-m, and 200-m grid spacing and 5, 10, 20, and 40 Gaussian quadrature integration points) where an oblique water–subsurface interface is given. The results show that, under certain conditions, the Gaussian quadrature integration method can generate relatively accurate and consistent wavefields compared to the conventional modelling method. Moreover, we performed 3D Laplace-domain inversion with this Gaussian quadrature modelling technique to obtain synthetic data with a flat and irregular interface using a 100-m and 200-m grid. The results obtained under the conditions herein show that 3D Laplace-domain waveform inversion with Gaussian quadrature integration can generate a more accurate and cost-efficient velocity model compared to that obtained using conventional Laplace-domain waveform inversion.

[1]  William W. Symes,et al.  A differential semblance algorithm for the inverse problem of reflection seismology , 1991 .

[2]  R. Higdon.  Absorbing boundary conditions for elastic waves , 1991 .

[3]  Equivalent source distribution for efficient 3-D acoustic wave equation modelling in the Laplace domain , 2011 .

[4]  William W. Symes,et al.  Velocity inversion by differential semblance optimization , 1991 .

[5]  Gary Cohen.  Higher-Order Numerical Methods for Transient Wave Equations , 2001 .

[6]  C. Officer,et al.  Introduction to the theory of sound transmission : with application to the ocean , 1958 .

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

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

[9]  R. Gerhard Pratt,et al.  Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies , 2004 .

[10]  J. Sochacki.  Absorbing boundary conditions for the elastic wave equations , 1988 .

[11]  W. Symes.  A Differential Semblance Criterion for Inversion of Multioffset Seismic Reflection Data , 1993 .

[12]  Mrinal K. Sen,et al.  Nonlinear multiparameter optimization using genetic algorithms; inversion of plane-wave seismograms , 1991 .

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

[14]  R. Schwarzenberger,et al.  The finite element method : a first approach , 1981 .

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

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

[17]  O. Zienkiewicz,et al.  Finite elements and approximation , 1983 .

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

[19]  Mrinal K. Sen,et al.  Nonlinear one-dimensional seismic waveform inversion using simulated annealing , 1991 .

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