Astronomy and Astrophysics – Astronomy
Scientific paper
Feb 2007
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2007geoji.168..605g&link_type=abstract
Geophysical Journal International, Volume 168, Issue 1, pp. 605-633.
Astronomy and Astrophysics
Astronomy
11
Born And Rytov Formulations, Diffraction Tomography, Finite Difference Methods, Medium Wavenumber Illumination, Seismic Imaging, Waveform Inversion
Scientific paper
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.
Gélis Céline
Grandjean G.
Virieux Jean
No associations
LandOfFree
Two-dimensional elastic full waveform inversion using Born and Rytov formulations in the frequency domain does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Two-dimensional elastic full waveform inversion using Born and Rytov formulations in the frequency domain, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Two-dimensional elastic full waveform inversion using Born and Rytov formulations in the frequency domain will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-738552