Full Waveform Inversion and the Truncated Newton Method

Abstract

Full waveform inversion (FWI) is a powerful method for reconstructing subsurface parameters from local measurements of the seismic wavefield. This method consists in minimizing the distance between predicted and recorded data. The predicted data are computed as the solution of a wave-propagation problem. Conventional numerical methods for the resolution of FWI problems are gradient-based methods, such as the preconditioned steepest descent, or more recently the $l$-BFGS quasi-Newton algorithm. In this study, we investigate the desirability of applying a truncated Newton method to FWI. The inverse Hessian operator plays a crucial role in the parameter reconstruction. The truncated Newton method allows one to better account for this operator. This method is based on the computation of the Newton descent direction by solving the corresponding linear system through an iterative procedure such as the conjugate gradient method. The large-scale nature of FWI problems requires us, however, to carefully implement this method to avoid prohibitive computational costs. First, this requires working in a matrix-free formalism and the capability of computing efficiently Hessian-vector products. To this purpose, we propose general second-order adjoint state formulas. Second, special attention must be paid to define the stopping criterion for the inner linear iterations associated with the computation of the Newton descent direction. We propose several possibilities and establish a theoretical link between the Steihaug--Toint method, based on trust regions, and the Eisenstat and Walker stopping criterion, designed for a method globalized by linesearch. We investigate the application of the truncated Newton method to two test cases: The first is a standard test case in seismic imaging based on the Marmousi II model. The second one is inspired by a near-surface imaging problem for the reconstruction of high-velocity structures. In the latter case, we demonstrate that the presence of large amplitude multiscattered waves prevents standard methods from converging while the truncated Newton method provides more reliable results.