Skip to Main content Skip to Navigation
Journal articles

FWT2D: A massively parallel program for frequency-domain full-waveform tomography of wide-aperture seismic data—Part 1: Algorithm

Florent Sourbier 1 Stéphane Operto 1 Jean Virieux 2 Patrick Amestoy 3 Jean-Yves l'Excellent 4
1 GEOAZUR 6526 - Géoazur
2 LGIT - Laboratoire de Géophysique Interne et Tectonophysique
3 ENSEEIHT - Ecole Nationale Supérieure d'Electrotechnique, d'Electronique, d'Informatique, d'Hydraulique et de Télécommunications
4 GRAAL - Algorithms and Scheduling for Distributed Heterogeneous Platforms
Inria Grenoble - Rhône-Alpes, LIP - Laboratoire de l'Informatique du Parallélisme
Abstract : This is the first paper in a two-part series that describes a massively parallel code that performs 2D frequency-domain full-waveform inversion of wide-aperture seismic data for imaging complex structures. Full-waveform inversion methods, namely quantitative seismic imaging methods based on the resolution of the full wave equation, are computationally expensive. Therefore, designing efficient algorithms which take advantage of parallel computing facilities is critical for the appraisal of these approaches when applied to representative case studies and for further improvements. Full-waveform modelling requires the resolution of a large sparse system of linear equations which is performed with the massively parallel direct solver MUMPS for efficient multiple-shot simulations. Efficiency of the multiple-shot solution phase (forward/backward substitutions) is improved by using the BLAS3 library. The inverse problem relies on a classic local optimization approach implemented with a gradient method. The direct solver returns the multiple-shot wavefield solutions distributed over the processors according to a domain decomposition driven by the distribution of the LU factors. The domain decomposition of the wavefield solutions is used to compute in parallel the gradient of the objective function and the diagonal Hessian, this latter providing a suitable scaling of the gradient. The algorithm allows one to test different strategies for multiscale frequency inversion ranging from successive mono-frequency inversion to simultaneous multifrequency inversion. These different inversion strategies will be illustrated in the following companion paper. The parallel efficiency and the scalability of the code will also be quantified.
Keywords : Seismic imaging Full-waveform inversion Parallel computation
Document type :
Journal articles
Complete list of metadatas
https://hal-insu.archives-ouvertes.fr/insu-00354703
Contributor : Pascale Talour <>
Submitted on : Tuesday, January 20, 2009 - 4:04:46 PM
Last modification on : Tuesday, December 8, 2020 - 10:51:13 AM

Identifiers

Collections

INSU | UNICE | CNRS | INRIA | UNIV-SAVOIE | OCA | GEOAZUR | LIP | UPMC | ENS-LYON | IFSTTAR | SIGAMM | UNIV-COTEDAZUR | UNIV-LYON1 | SORBONNE-UNIVERSITE | UNIV-TLSE3 | IRSTEA | UDL | INRAE | UGA

Citation

Florent Sourbier, Stéphane Operto, Jean Virieux, Patrick Amestoy, Jean-Yves l'Excellent. FWT2D: A massively parallel program for frequency-domain full-waveform tomography of wide-aperture seismic data—Part 1: Algorithm. Computers & Geosciences, Elsevier, 2009, 35 (3), pp.487-495. ⟨10.1016/j.cageo.2008.04.013⟩. ⟨insu-00354703⟩

Export

BibTeX TEI DC DCterms EndNote

Share

Metrics

Record views

765