Extended full waveform inversion in the time domain by the augmented Lagrangian method

Abstract

Extended full-waveform inversion (FWI) has shown promising results for accurate estimation of subsurface parameters when the initial models are not sufficiently accurate. Frequency-domain applications have shown that the augmented Lagrangian (AL) method solves the inverse problem accurately with a minimal effect of the penalty parameter choice. Applying this method in the time domain, however, is limited by two main factors: (1) The challenge of data-assimilated wavefield reconstruction due to the lack of an explicit time-stepping and (2) The need to store the Lagrange multipliers, which is not feasible for the field-scale problems. We show that these wavefields are efficiently determined from the associated data (projection of the wavefields onto the receivers space) by using explicit time stepping. Accordingly, based on the augmented Lagrangian, a new algorithm is proposed which performs in "data space" (a lower dimensional subspace of the full space) in which the wavefield reconstruction step is replaced by reconstruction of the associated data, thus requiring optimization in a lower dimensional space (convenient for handling the Lagrange multipliers). We show that this new algorithm can be implemented efficiently in the time domain with existing solvers for the FWI and at a cost comparable to that of the FWI while benefiting from the robustness of the extended FWI formulation. The results obtained by numerical examples show high-performance of the proposed method for large scale time-domain FWI.

References in corpus (6)

• A penalty method for PDE-constrained optimization in inverse problems
• Improving full waveform inversion by wavefield reconstruction with the alternating direction method of multipliers
• Compound Regularization of Full-waveform Inversion for Imaging Piecewise Media
• Implementing bound constraints and total-variation regularization in extended full waveform inversion with the alternating direction method of multiplier: application to large contrast media
• Full Waveform Inversion with Adaptive Regularization
• Wavefield Reconstruction Inversion: an example

Cited by in corpus (5)

• Full Waveform Inversion and Lagrange Multipliers
• Improving full-waveform inversion based on sparse regularisation for geophysical data
• Inexact Augmented Lagrangian Method-Based Full-waveform Inversion with Randomized Singular Value Decomposition
• Algorithmic analysis towards time-domain extended source waveform inversion
• Enhanced Geological Prediction for Tunnel Excavation Using Full Waveform Inversion Integrating Sobolev Space Regularization with a Quadratic Penalty Method