In elastic imaging, the extrapolated vector fields are decomposed into pure wave modes, such that the imaging condition produces interpretable images, which characterize reflectivity of different reflection types. Conventionally, wavefield decomposition in anisotropic media is costly as the operators involved is dependent on the velocity, and thus not stationary. In this abstract, we propose an efficient approach to directly extrapolate the decomposed elastic waves using lowrank approximate mixed space/wavenumber domain integral operators for heterogeneous transverse isotropic (TI) media. The low-rank approximation is, thus, applied to the pseudospectral extrapolation and decomposition at the same time. The pseudo-spectral implementation also allows for relatively large time steps in which the low-rank approximation is applied. Synthetic examples show that it can yield dispersionfree extrapolation of the decomposed quasi-P (qP) and quasi- SV (qSV) modes, which can be used for imaging, as well as the total elastic wavefields.
Bibliographical noteKAUST Repository Item: Exported on 2020-10-01
Acknowledgements: The first author would like to thank Sergey Fomel for sharing his experience in designing low-rank approximate algorithm, and Chenlong Wang and Tengfei Wang for helpful discussion. We thank KAUST for it’s support. We thank the authors of Madagascar for providing this software platform for reproducible computational experiments.