Abstract:Density is difficult to reconstruct in multi-parameter full waveform inversion (FWI). This difficulty mainly comes from the coupling effects between velocity and density. In this study, we apply the truncated Gauss-Newton method to multi-parameter FWI for simultaneous estimation of velocity and density in acoustic media. The inverse Hessian whose off-diagonal blocks reflect the coupling effects between different parameters is incorporated to decouple velocity and density during the inversion procedure. In conventional truncated Newton algorithms, the gradients of the objective function with respect to model parameters are calculated with the first-order adjoint-state method, where the forward propagated source wavefields are zero-lag cross-correlated with the back-propagated residual wavefields. The model update is computed through a matrix-free conjugate gradient (CG) solution of the Newton linear system. This operation requires the computation of the Hessian-vector products, which can be acquired by the second-order adjoint-state formulas or the finite difference method. In such a way, two additional forward wave-propagation problems have to be solved in each internal CG loop, which increase the computational cost dramatically. In this study, we calculate the gradient based on the Born sensitivity kernels, where the kernel-vector products are transformed to the accumulation of vector-scalar products. Meanwhile, the Hessian vector products are implemented by means of two nested kernel-vector products, thus the forward simulations are no longer required. To demonstrate the effectiveness of the proposed method, we apply it to a canonical inclusion model and the SEG/EAGE Marmousi-2 model to simultaneously determine velocity and density. The results show that this method can obtain reasonable results for both velocity and density without using any inversion strategy. The vertical profiles extracted from the reconstructed models depict good consistency with the true ones. Accurate consideration of the inverse Hessian is of great importance in the context of multi parameter FWI, which helps to mitigate the coupling effects and rescale the magnitudes of different parameters. The Hessian-vector products can be efficiently computed using kernel-vector products without solving any additional forward simulation problems, thus the computational cost of the truncated Gauss-Newton FWI is reduced. Numerical experiments prove the capability of the proposed method on recovering multiple parameters simultaneously.

Finite bandwidth simultaneous inversion of three parameters in stratified elastic media

A Multi-parameter Full Waveform Inversion Strategy in Acoustic Media

Full Waveform Inversion Method for Horizontally Inhomogeneous Stratified Media

Simultaneous inversion of layered velocity and density profiles using Direct Waveform Inversion (DWI): 1D case

A Multi-Parameter Full Waveform Inversion Strategy for Acoustic Media with Variable Density

Application of the finite-difference contrast source inversion method to multiparameter reconstruction using seismic full-waveform data

Prestack waveform inversion based on analytical solution of the viscoelastic wave equation

Joint Laplace- and Fourier-domain Seismic Inversion for Orthorhombic Medium Parameter Estimates

Incorporating attenuation effects into frequency-domain full waveform inversion from zero-offset VSP data from the Stokes equation

Study on Waveform Inversion Method for Elastic Parameters of Stratified Medium and Its Applications

Simultaneous Estimation of Velocity and Density in Acoustic Multiparameter Full-Waveform Inversion Using an Improved Scattering-Integral Approach

Seismic nonlinear frequency-dependent inversion for viscoelastic media based on a modified scattering coefficient approximation

Elastic multi-parameter full-waveform inversion for anisotropic media

Full-Wavefield Surface Wave Method: Integrating Rigorous and Efficient Methods for Enhanced Subsurface Exploration

Acoustic- and elastic-waveform inversion with total generalized p-variation regularization

Multi-parameter full waveform inversion for velocity and density based on Born sensitivity kernels

P-P and Dynamic Time Warped P-SV Wave AVA Joint-Inversion with ℓ 1–2 Regularization

An Avo Inversion Method in the Frequency Domain Based on A Layered Model-A Bakken Shale Case Study

Elastic Model Low- to Intermediate-Wavenumber Inversion Using Reflection Traveltime and Waveform of Multicomponent Seismic Data

Inversion of Fracture Weakness Parameters Based on the 3CVSP Data?Part I: HTI Media Composed of A Single Fracture Set in the Isotropic Background Media

Regularized Seismic Velocity Inversion Based on Inverse Scattering Theory