Accéder directement au contenu Accéder directement à la navigation
Article dans une revue

Two-dimensional elastic full waveform inversion using Born and Rytov formulations in the frequency domain

Abstract : We perform the full elastic waveform inversion in the frequency domain in a 2-D geometry. This method allows imaging of two physical seismic parameters, using vertical and horizontal field components. The forward problem is discretized using finite difference, allowing to simulate the full elastic wavefield propagation. Moreover, it is solved in the frequency domain, a fast approach for multisource and multireceiver acquisition. The non-linear inversion is based on a pre-conditioned gradient method, where Born and Rytov formulations are used to compute Fréchet derivatives. Parameter perturbations linearly depend on fields perturbations in the Born kernel, and on the generalized complex phases of fields in the Rytov kernel, giving different Fréchet derivatives. The gradient is pre-conditioned with the diagonal part of the inverse Hessian matrix, allowing to better estimate the stepping in the optimization direction. Non-linearity is taken into account by updating parameters at each iteration and proceeding from low to high frequencies. The latter allows as well to progressively introduce smaller wavelengths in parameter images. On a very simple synthetic example, we examine the way the inversion determines the Vp (P-wave velocity) and Vs (S-wave velocity) images. We highlight that, with a transmission acquisition, final parameter images weakly depend on the chosen formulation to compute Fréchet derivatives and on the inverted parameters choice. Of course, convergence strongly depends on the medium wavenumber illumination which is related somehow to the acquisition geometry. With a reflection acquisition, the Born formulation allows to better recover scatterers. Moreover, the medium anomalies are not well reconstructed when surface waves propagate in the medium. This may be due to the evanescent nature of surface waves. By selecting first body waves and then surface waves, we improve the convergence and properly reconstruct anomalies. This shows us that preparation of the seismic data before the inversion is as critical as the initial model selection.
Type de document :
Article dans une revue
Liste complète des métadonnées
Contributeur : Pascale Talour Connectez-vous pour contacter le contributeur
Soumis le : jeudi 25 mars 2021 - 17:55:34
Dernière modification le : jeudi 4 août 2022 - 16:59:17
Archivage à long terme le : : samedi 26 juin 2021 - 19:16:57


Fichiers éditeurs autorisés sur une archive ouverte




Céline Gélis, Jean Virieux, Gilles Grandjean. Two-dimensional elastic full waveform inversion using Born and Rytov formulations in the frequency domain. Geophysical Journal International, Oxford University Press (OUP), 2007, 168 (2), pp.605 à 633. ⟨10.1111/j.1365-246X.2006.03135.x⟩. ⟨insu-00355241⟩



Consultations de la notice


Téléchargements de fichiers