2022/04/13 by Malovichko, M., Orazbayev, A., Khokhlov, N.
#FOS: Computer and information sciences #FOS: Mathematics #FOS: Physical sciences #Geophysics (physics.geo-ph) #Mathematical Software (cs.MS) #Numerical Analysis (math.NA)
paper · doi:10.48550/arxiv.2204.06489
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.