
Frequency‐domain elastic full waveform inversion for VTI media
Author(s) -
Lee HoYong,
Koo June Mo,
Min DongJoo,
Kwon ByungDoo,
Yoo Hai Soo
Publication year - 2010
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1111/j.1365-246x.2010.04767.x
Subject(s) - transverse isotropy , inversion (geology) , geology , isotropy , waveform , frequency domain , anisotropy , time domain , wave equation , inverse theory , mathematical analysis , geophysics , geometry , structural basin , physics , mathematics , optics , computer science , geomorphology , oceanography , deformation (meteorology) , voltage , quantum mechanics , computer vision
SUMMARY To describe subsurface structures in anisotropic media properly, particularly in transversely isotropic media with a vertical symmetry axis (VTI), which frequently appear in sedimentary basin environments, we develop a frequency‐domain elastic full waveform inversion algorithm for 2‐D VTI media. The inversion algorithm is based on the cell‐based finite‐difference modelling method and the adjoint state of the wave equation. Because the anisotropic inversion for VTI media deals with more elastic constants than the isotropic inversion, it is more prone to obtain local minimum solutions. For this reason, we may not succeed in properly describing the elastic constants of subsurface media if we only apply the standard inversion techniques to anisotropic waveform inversion. To compensate for the ill‐posedness of the anisotropic waveform inversion, we couple elastic constants C 11 and C 33 based on Thomsen's relationship, which is also supported by the sensitivity analysis with respect to the parameters. To enhance the inversion results, we apply the frequency‐selection strategy, moving from lower to higher frequencies and we carry out the inversion process over two stages. In both stages, all of the elastic constants are simultaneously optimized, as is done in the conventional waveform inversion. However, we only accept the inversion results for C 11 , C 33 and C 44 at the first stage, which will be used as the starting models for the second stage and C 13 is reinitialized as a linearly increasing model. We apply our waveform inversion algorithm to a simple 3‐layered model and a part of the overthrust model. For the 3‐layered model, the first inversion stage is enough to yield reasonable inversion results for all of the elastic constants. For the overthrust model, the second stage is needed to enhance the inversion results for C 13 .