Imaging of complex media with elastic wave equations - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Poster Année : 2014

Imaging of complex media with elastic wave equations

Résumé

Even if RTM has enjoyed the tremendous progresses of scientific computing, its performances can still be improved, in particular when applied to strong heterogeneous media. In this case, images have been mainly obtained by using direct arrivals of acoustic waves and the transition to elastic waves including multiples is not obvious essentially because elastic waves equations are still more computationnaly consuming. The accuracy of numerical wave fields is obviously of great importance. We have thus chosen to consider high-order Discontinuous Galerkin Methods which are known to be well-adapted to provide accurate solutions based upon parallel computing. Now one of the main drawback of RTM is the need of storing a huge quantity of information which is redhibitory when using elastic waves. However RTM procedure can be modified by giving up the idea of applying steps (i)-(iii) sequentially. For that purpose, we apply the Griewank algorithm following Symes' ideas for the acoustic RTM. The idea is to find a compromise between the number of wave equations to solve and the number of numerical waves that we have to store. This is the so-called Optimal Checkpointing. By reducing the occupancy of the memory, RTM should be efficient even when using elastic waves. By this way, the question of knowing if considering elastic waves including multiples allow to improve images of heterogeneous media could be considered. It must involved a careful numerical analysis including the evaluation of the impact of the imaging condition. It is thus necessary to derive accurate imaging conditions, which could takeadvantage of all the information contained in the wavefield. For acoustic media, Claerbout proposed an imaging condition which is widely used and turns out to be sufficient to accurately reproduce interfaces. But Claerbout conditions do not take wave conversions into account and it is not clear if conversions do or do not not contain interesting information to get accurate images of heterogeneous media. However, since P-wave and S-wave interact with each other, it might be relevant to use an imaging condition including these interactions. In fact, this has been done successfully by J.Tromp and C. Morency for seismology applications based upon the inversion of the global Earth. Their approach is based upon the state adjoint and it involves sensitivity kernels which are defined from the propagated and the backpropagated fields. Now it has been shown that full wave form inversions using these sensitivity kernels may be polluted by numerical artefacts. One solution is to use a linear combination of the sensitivity kernels to delete artefacts. In this work, we propose then a new imaging condition which construction is inspired from with some approximations required to keep admissible computational costs. We illustrate the properties of the new imaging condition on industrial benchmarks like the Marmousi model. In particular, we compare the new imaging condition with other imaging conditions by using as criteria the quality of the image and the computational costs required by the RTM.The authors acknowledge the support by the Inria-Total strategic action DIP (dip.inria.fr).
Fichier principal
Vignette du fichier
poster_doctoriales_en.pdf (2.57 Mo) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)

Dates et versions

hal-01096620 , version 1 (30-12-2014)

Identifiants

  • HAL Id : hal-01096620 , version 1

Citer

Jérôme Luquel, Hélène Barucq, Julien Diaz, Henri Calandra. Imaging of complex media with elastic wave equations. SIAM Conference on IMAGING SCIENCE, May 2014, Hong Kong, China. 2014. ⟨hal-01096620⟩
377 Consultations
105 Téléchargements

Partager

Gmail Facebook X LinkedIn More