The truncated Newton method for Full Waveform Inversion
- Others:
- Institut des Sciences de la Terre (ISTerre) ; Université Joseph Fourier - Grenoble 1 (UJF)-Institut Français des Sciences et Technologies des Transports, de l'Aménagement et des Réseaux (IFSTTAR)-Institut national des sciences de l'Univers (INSU - CNRS)-Institut de recherche pour le développement [IRD] : UR219-PRES Université de Grenoble-Université Savoie Mont Blanc (USMB [Université de Savoie] [Université de Chambéry])-Centre National de la Recherche Scientifique (CNRS)
- Géoazur (GEOAZUR 6526) ; Institut de Recherche pour le Développement (IRD)-Université Pierre et Marie Curie - Paris 6 (UPMC)-Université Nice Sophia Antipolis (1965 - 2019) (UNS) ; COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire de la Côte d'Azur ; COMUE Université Côte d'Azur (2015-2019) (COMUE UCA)-Université Côte d'Azur (UCA)-Université Côte d'Azur (UCA)-Centre National de la Recherche Scientifique (CNRS)
Description
Full Waveform Inversion (FWI) is a promising seismic imaging method. It aims at computing quantitative estimates of the subsurface parameters (bulk wave velocity, shear wave velocity, rock density) from local measurements of the seismic wavefield. Based on a particular wave propagation engine for wavefield estimation, it consists in minimizing iteratively the distance between the predicted wavefield at the receivers and the recorded data. This amounts to solving a strongly nonlinear large scale inverse problem. This minimization is generally performed using gradient-based methods. We investigate the possibility of applying the truncated Newton (TrN) method to this problem. This is done through the development of general second-order adjoint state formulas that yield an efficient algorithm to compute Hessian-vector products, and the design of an adaptive stopping criterion for the inner conjugate gradient (CG) iterations. Numerical results demonstrate the interest of using the TrN method when multi-scattered waves dominate the recorded data.
Abstract
Extended abstract for the NCMIP 2012 Conference in Cachan
Abstract
International audience
Additional details
- URL
- https://hal.science/hal-00763730
- URN
- urn:oai:HAL:hal-00763730v1
- Origin repository
- UNICA