Two-dimensional elastic full waveform inversion using Born and Rytov formulations in the frequency domain - BRGM - Bureau de recherches géologiques et minières Accéder directement au contenu
Article Dans Une Revue Geophysical Journal International Année : 2007

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

Résumé

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.
Fichier principal
Vignette du fichier
168-2-605.pdf (9.71 Mo) Télécharger le fichier
Origine : Fichiers éditeurs autorisés sur une archive ouverte

Dates et versions

insu-00355241 , version 1 (25-03-2021)

Identifiants

Citer

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, 2007, 168 (2), pp.605 à 633. ⟨10.1111/j.1365-246X.2006.03135.x⟩. ⟨insu-00355241⟩
365 Consultations
124 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More