Loading…

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 prob...

Full description

Saved in:
Bibliographic Details
Published in:Computational mathematics and mathematical physics 2024, Vol.64 (5), p.954-966
Main Authors: Malovichko, M. S., Orazbayev, A., Khokhlov, N. I., Petrov, I. B.
Format: Article
Language:English
Subjects:
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary: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.
ISSN:0965-5425
1555-6662
DOI:10.1134/S0965542524700192