The truncated Newton method for Full Waveform Inversion

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

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) methods use generally gradient based method, such as the nonlinear conjugate gradient method or more recently the l-BFGS quasi-Newton method. Several authors have already investigated the possibility of accounting more accurately for the inverse Hessian operator in the minimization scheme through Gauss-Newton or exact Newton algorithms. We propose a general framework for the implementation of these methods inside a truncated Newton procedure. We demonstrate that the exact Newton method can outperform the standard gradient-based methods in a near-surface application case for recovering high-velocity concrete structures. In this particular configuration, large amplitude multi-scattered waves are generated, which are better taken into account using the exact-Newton method.

Mots clés: nonlinear 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

Ludovic Métivier, Romain Brossier, Jean Virieux, Stéphane Operto. The truncated Newton method for Full Waveform Inversion. SEG 2012, Nov 2012, Las Vegas, United States. ⟨hal-00826613⟩