Iterative PDE-constrained optimization for seismic full-waveform inversion

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.