Skip to Main content Skip to Navigation
Journal articles

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 Amestoy 1 Romain Brossier 2 Alfredo Buttari 1, 3 Jean-Yves l'Excellent 4, 5 Théo Mary 6 Ludovic Métivier 7, 2 Alain Miniussi 8 Stéphane Operto 8
1 IRIT-APO - Algorithmes Parallèles et Optimisation
IRIT - Institut de recherche en informatique de Toulouse
2 ISTerre - Institut des Sciences de la Terre
3 CNRS - Centre National de la Recherche Scientifique
4 LIP - Laboratoire de l'Informatique du Parallélisme
5 ROMA - Optimisation des ressources : modèles, algorithmes et ordonnancement
Inria Grenoble - Rhône-Alpes, LIP - Laboratoire de l'Informatique du Parallélisme
6 IRIT - Institut de recherche en informatique de Toulouse
7 EDP - Equations aux Dérivées Partielles
LJK - Laboratoire Jean Kuntzmann
8 GEOAZUR 7329 - Géoazur
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 : vertical transverse isotropy finite difference block low-rank multifrontal method Wave equation frequency domain
Document type :
Journal articles
Complete list of metadata
https://hal.inria.fr/hal-01349119
Contributor : Equipe Roma Connect in order to contact the contributor
Submitted on : Tuesday, July 26, 2016 - 6:04:14 PM
Last modification on : Saturday, January 15, 2022 - 3:52:21 AM

Identifiers

Collections

CNRS | OSUG | OCA | INRIA | ENS-LYON | ISTERRE | LJK | UNIV-SAVOIE | SMS | UNIV-TLSE2 | IFSTTAR | INSU | LJK_MAD | LJK_MAD_EDP | UNIV-COTEDAZUR | UNIV-LYON1 | UNIV-TLSE3 | UT1-CAPITOLE | UDL | SIGAMM | IRIT | IRIT-APO | GEOAZUR | UGA | IRIT-CISO | IRIT-CNRS

Citation

Patrick 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⟩

Export

BibTeX TEI DC DCterms EndNote

Share

Metrics

Les métriques sont temporairement indisponibles

 

Contact                    Legal Notice                    Personal Data