Fast 3D frequency-domain full waveform inversion with a parallel Block Low-Rank multifrontal direct solver: application to OBC data from the North Sea

Patrick R. Amestoy 1, 2 Romain Brossier 3 Alfredo Buttari 4, 2 Jean-Yves L'Excellent 5, 6 Théo Mary 7, 2 Ludovic Métivier 8, 3, 4 Alain Miniussi 9 Stéphane Operto 9
1 INP de Toulouse - Institut Nationale Polytechnique de Toulouse
2 IRIT - Institut de recherche en informatique de Toulouse
3 ISTerre - Institut des Sciences de la Terre
4 CNRS - Centre National de la Recherche Scientifique
5 LIP - Laboratoire de l'Informatique du Parallélisme
6 ROMA - Optimisation des ressources : modèles, algorithmes et ordonnancement
Inria Grenoble - Rhône-Alpes, LIP - Laboratoire de l'Informatique du Parallélisme
7 UPS - Université Toulouse III - Paul Sabatier
8 EDP - Equations aux Dérivées Partielles
LJK - Laboratoire Jean Kuntzmann
9 GEOAZUR - Observatoire de la Côte d'Azur
Abstract : Wide-azimuth long-offset OBC/OBN surveys provide a suitable framework to perform computationally-efficient frequency-domain full waveform inversion (FWI) with a few discrete frequencies. Frequency-domain seismic modeling is performed efficiently with moderate computational resources for a large number of sources with a sparse multifrontal direct solver (Gauss-elimination techniques for sparse matrices). Approximate solutions of the time-harmonic wave equation are computed using a Block Low-Rank (BLR) approximation, leading to a significant reduction in the operation count and in the volume of communication during the LU factorization as well as offering a great potential for reduction in the memory demand. Moreover, the sparsity of the seismic source vectors is exploited to speed up the forward elimination step during the computation of the monochromatic wavefields. The relevance and the computational efficiency of the frequency-domain FWI performed in the visco-acoustic VTI approximation is shown with a real 3D OBC case study from the North Sea. The FWI subsurface models show a dramatic resolution improvement relative to the initial model built by reflection traveltime tomography. The amplitude errors introduced in the modeled wavefields by the BLR approximation for different low-rank thresholds have a negligible footprint in the FWI results. With respect to a standard multifrontal sparse direct factorization, and without compromise on the accuracy of the imaging, the BLR approximation can bring a reduction of the LU factor size by a factor up to three. This reduction is not yet exploited to reduce the effective memory usage (ongoing work). The flop reduction can be larger than a factor of 10 and can bring a factor of time reduction of around three. Moreover, this reduction factor tends to increase with frequency, namely with the matrix size. Frequency-domain visco-acoustic VTI FWI can be viewed as an efficient tool to build an initial model for elastic FWI of 4-C OBC data.
Keywords : finite difference vertical transverse isotropy block low-rank frequency domain Wave equation multifrontal method
Type de document :
Article dans une revue
Geophysics, Society of Exploration Geophysicists, 2016, 81 (6), pp.R363-R383. 〈10.1190/geo2016-0052.1〉
Liste complète des métadonnées
https://hal.inria.fr/hal-01349119
Contributeur : Equipe Roma <>
Soumis le : mardi 26 juillet 2016 - 18:04:14
Dernière modification le : vendredi 5 octobre 2018 - 16:52:04

Identifiants

Collections

UNICE | OSUG | OCA | LIP | INRIA | ENS-LYON | ISTERRE | LJK | UNIV-SAVOIE | UNIV-TLSE3 | SMS | UNIV-TLSE2 | IFSTTAR | INSU | UGA | LJK_MAD | LJK_MAD_EDP | UNIV-COTEDAZUR | UNIV-LYON1

Citation

Patrick R. Amestoy, Romain Brossier, Alfredo Buttari, Jean-Yves L'Excellent, Théo Mary, et al.. Fast 3D frequency-domain full waveform inversion with a parallel Block Low-Rank multifrontal direct solver: application to OBC data from the North Sea. Geophysics, Society of Exploration Geophysicists, 2016, 81 (6), pp.R363-R383. 〈10.1190/geo2016-0052.1〉. 〈hal-01349119〉

Exporter

BibTeX TEI DC DCterms EndNote

Partager

Métriques

Consultations de la notice

729

Contact


archive-ouverte@inria.fr