ADMM-based multiparameter wavefield reconstruction inversion in VTI acoustic media with TV regularization

Abstract

Full waveform inversion (FWI) is a nonlinear waveform matching procedure, which suffers from cycle skipping when the initial model is not kinematically-accurate enough. To mitigate cycle skipping, wavefield reconstruction inversion (WRI) extends the inversion search space by computing wavefields with a relaxation of the wave equation in order to fit the data from the first iteration. Then, the subsurface parameters are updated by minimizing the source residuals the relaxation generated. Capitalizing on the wave-equation bilinearity, performing wavefield reconstruction and parameter estimation in alternating mode decomposes WRI into two linear subproblems, which can solved efficiently with the alternating-direction method of multiplier (ADMM), leading to the so-called IR-WRI. Moreover, ADMM provides a suitable framework to implement bound constraints and different types of regularizations and their mixture in IR-WRI. Here, IR-WRI is extended to multiparameter reconstruction for VTI acoustic media. To achieve this goal, we first propose different forms of bilinear VTI acoustic wave equation. We develop more specifically IR-WRI for the one that relies on a parametrisation involving vertical wavespeed and Thomsen's parameters δ and . With a toy numerical example, we first show that the radiation patterns of the virtual sources generate similar wavenumber filtering and parameter cross-talks in classical FWI and IR-WRI. Bound constraints and TV regularization in IR-WRI fully remove these undesired effects for an idealized piecewise constant target. We show with a more realistic long-offset case study representative of the North Sea that anisotropic IR-WRI successfully reconstruct the vertical wavespeed starting from a laterally homogeneous model and update the long-wavelengths of the starting model, while a smooth δ model is used as a passive background model. VTI acoustic IR-WRI can be alternatively performed with subsurface parametrisations involving stiffness or compliance coefficients or normal moveout velocities and η parameter (or horizontal velocity).

Graph‐Space Optimal Transport Concept for Time‐Domain Full‐Waveform Inversion of Ocean‐Bottom Seismometer Data: Nankai Trough Velocity Structure Reconstructed From a 1D Model

“…MCS data can be utilized to mitigate this issue in the manner mentioned above. One can also consider more robust gradient regularizations, which are able to produce a broadband wavenumber update of the model, or which can utilize a priori information about the underlying structure (H. S. Aghamiry et al., 2019; Peters & Herrmann, 2017; Trinh et al., 2017).…”

Receiver-extension strategy for time-domain full-waveform inversion using a relocalization approach

“…This identity shows that the wave propagation problem can be rewritten equally as a linear operator A(m) acting on u or a linear operator B(u) acting on m. This identity is useful in the following developments. This property is true for general elastic and visco-elastic wave propagation, up to the choice of the parameterization for m, as is discussed in Aghamiry et al (2019a).…”

Ensemble-based uncertainty estimation in Full Waveform Inversion

“…Besides, this method does seem to display challenging memories limitation as it requires to store the optimal wavefields in memory for each frequency bands, which may become challenging for large scale 3D application. The extension to uncertainty estimation of multiparameter inversion also seems to be non-trivial in this extended domain FWI application, as only recent publications are addressing the multiparameter aspect of wavefield reconstruction inversion (Aghamiry et al 2019).…”