Application of a preconditoned truncated Newton method to full waveform inversion

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

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

Full Waveform Inversion (FWI) is a powerful seismic imaging method, based on the iterative minimization of the distance between simulated and recorded wavefields. The inverse Hessian operator related to this misfit function plays an important role in the reconstruction scheme. As conventional methods use direct approximations of this operator, we investigate an alternative optimization scheme: the truncated Newton method. This two-nested-loops algorithm is based on the resolution of the Newton linear system through a matrix-free iterative solver at each outer iteration. On the 2D BP 2004 model widely used as a benchmark for FWI, the contrasts in wave velocities between salt structures and the upper water layer generate high amplitude multiple reflections. These multiple reflections strengthen the need for quite accurate approximation of the inverse Hessian operator and the truncated Newton method is shown to outperform than more conventional algorithms (l-BFGS, nonlinear conjugate gradient).

Mots clés: numerical optimization; Newton method; seismic imaging; Full Waveform Inversion

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

Romain Brossier, Ludovic Métivier, Stéphane Operto, Jean Virieux. Application of a preconditoned truncated Newton method to full waveform inversion. WAVES 2013 - 11th International Conference on Mathematical and Numerical Aspects of Waves, Jun 2013, Tunis, Tunisia. ⟨hal-00826610⟩