Fast full waveform inversion with source encoding and second-order optimization methods - EDP Accéder directement au contenu
Article Dans Une Revue Geophysical Journal International Année : 2015

Fast full waveform inversion with source encoding and second-order optimization methods

Clara Castellanos
  • Fonction : Auteur
Stéphane Operto
Romain Brossier
Jean Virieux
  • Fonction : Auteur
  • PersonId : 857422

Résumé

Full waveform inversion (FWI) of 3-D data sets has recently been possible thanks to the development of high performance computing. However, FWI remains a computationally intensive task when high frequencies are injected in the inversion or more complex wave physics (viscoelastic) is accounted for. The highest computational cost results from the numerical solution of the wave equation for each seismic source. To reduce the computational burden, one well-known technique is to employ a random linear combination of the sources, rather that using each source independently. This technique, known as source encoding, has shown to successfully reduce the computational cost when applied to real data. Up to now, the inversion is normally carried out using gradient descent algorithms. With the idea of achieving a fast and robust frequency-domain FWI, we assess the performance of the random source encoding method when it is interfaced with second-order optimization methods (quasi-Newton l-BFGS, truncated Newton). Because of the additional seismic modelings required to compute the Newton descent direction, it is not clear beforehand if truncated Newton methods can indeed further reduce the computational cost compared to gradient algorithms. We design precise stopping criteria of iterations to fairly assess the computational cost and the speed-up provided by the source encoding method for each optimization method. We perform experiment on synthetic and real data sets. In both cases, we confirm that combining source encoding with second-order optimization methods reduces the computational cost compared to the case where source encoding is interfaced with gradient descent algorithms. For the synthetic data set, inspired from the geology of Gulf of Mexico, we show that the quasi-Newton l-BFGS algorithm requires the lowest computational cost. For the real data set application on the Valhall data, we show that the truncated Newton methods provide the most robust direction of descent.
Fichier principal
Vignette du fichier
ggu427.pdf (10.47 Mo) Télécharger le fichier
Origine : Fichiers éditeurs autorisés sur une archive ouverte

Dates et versions

hal-01888507 , version 1 (08-11-2021)

Licence

Paternité

Identifiants

Citer

Clara Castellanos, Ludovic Métivier, Stéphane Operto, Romain Brossier, Jean Virieux. Fast full waveform inversion with source encoding and second-order optimization methods. Geophysical Journal International, 2015, 200 (2), pp.720 - 744. ⟨10.1093/gji/ggu427⟩. ⟨hal-01888507⟩
107 Consultations
18 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More