Full Waveform Inversion and the truncated Newton method

Ludovic Métivier (1, 2) , Romain Brossier (2) , Jean Virieux (2) , Stéphane Operto (3)

1 EDP - Equations aux Dérivées Partielles
2 ISTerre - Institut des Sciences de la Terre
3 GEOAZUR 6526 - Géoazur

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 propose several possibilities and establish a theoretical link between the Steihaug-Toint method, based on trust-regions, and the Eisenstat stopping criterion, designed for 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 multi-scattered waves prevents standard methods from converging while the truncated Newton method provides more reliable results.

Mots clés: large-scale inverse problems; Seismic imaging; Full Waveform Inversion; numerical optimization

Domaines: Optimisation et contrôle [math.OC]; Géophysique [physics.geo-ph]; Géophysique [physics.geo-ph]; Milieux et Changements globaux

Ludovic Métivier, Romain Brossier, Jean Virieux, Stéphane Operto. Full Waveform Inversion and the truncated Newton method. SIAM Journal on Scientific Computing, 2013, 35 (2), pp.B401-B437. ⟨10.1137/120877854⟩. ⟨hal-00763702⟩