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

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

Florent Sourbier ¹ Stéphane Operto ¹ Jean Virieux ² Patrick Amestoy ³ Jean-Yves l'Excellent ⁴

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 : Friday, October 11, 2019 - 8:22:56 PM

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

Identifiers

Collections

Citation

Export

Metrics

683

Contact

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

Identifiers

Collections

Citation

Export

Share

Metrics

683

Contact