Iterative PDE-constrained optimization for seismic full-waveform inversion

DOI10.1134/S0965542524700192MaRDI QIDQ6396447zbMATH OpenOpenAlexFDO

Authors M. S. Malovichko, A. Orazbayev, Nikolai I. Khokhlov, I. B. Petrov

Publication date 13 April 2022

Full work available at URL https://doi.org/10.1134/s0965542524700192

Mathematics Subject Classification ID

Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Inverse problems for PDEs (35R30) Newton-type methods (49M15) Numerical optimization and variational techniques (65K10) Finite difference methods for boundary value problems involving PDEs (65N06) Seismology (including tsunami modeling), earthquakes (86A15)

Abstract: This paper presents a novel numerical method for the Newton seismic full-waveform inversion (FWI). The method is based on the full-space approach, where the state, adjoint state, and control variables are optimized simultaneously. Each Newton step is formulated as a PDE-constrained optimization problem, which is cast in the form of the Karush-Kuhn-Tucker (KKT) system of linear algebraic equitations. The KKT system is solved inexactly with a preconditioned Krylov solver. We introduced two preconditioners: the one based on the block-triangular factorization and its variant with an inexact block solver. The method was benchmarked against the standard truncated Newton FWI scheme on a part of the Marmousi velocity model. The algorithm demonstrated a considerable runtime reduction compared to the standard FWI. Moreover, the presented approach has a great potential for further acceleration. The central result of this paper is that it establishes the feasibility of Newton-type optimization of the KKT system in application to the seismic FWI.

Recommendations

• An adjoint-based Jacobi-type iterative method for elastic full waveform inversion problem
• A class of preconditioned conjugate gradient methods for solving full-waveform inversion problem
• Full-wavefield inversion: an extreme-scale PDE-constrained optimization problem
• A Newton-CG method for large-scale three-dimensional elastic full-waveform seismic inversion
• Full waveform inversion and the truncated Newton method

Cited in

(1)

• A general approach to computing derivatives for Hessian-based seismic inversion

This page was built for publication: Iterative PDE-constrained optimization for seismic full-waveform inversion

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6396447)