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 a distance between predicted and recorded data. The predicted data is 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 interest 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 however to carefully implement this method to avoid prohibitive computational costs. First, this requires to work 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 payed to define the stopping criterion for the inner linear iterations associated with the computation of the Newton descent direction. We propos