Elastic Wavefield Separation in Anisotropic Media Based on Eigenform Analysis and Its Application in Reverse-Time Migration




Journal Title

Journal ISSN

Volume Title


Oxford University Press



Separating compressional and shear wavefields is an important step in elastic reverse-time migration, which can remove wave-mode crosstalk artefacts and improve imaging quality. In vertical (VTI) and titled (TTI) transversely isotropic media, the state-of-the-art techniques for wavefield separation are based on either non-stationary filter or low-rank approximation. They both require intensive Fourier transforms for models with strong heterogeneity. Based on the eigenform analysis, we develop an efficient pseudo-Helmholtz decomposition method for the VTI and TTI media, which produces vector P and S wavefields with the same amplitudes, phases and physical units as the input elastic wavefields. Starting from the elastic VTI wave equations, we first derive the analytical eigenvalues and eigenvectors, then use the Taylor expansion to approximate the square-root term in the eigenvalues, and finally obtain a zero-order and a first-order pseudo-Helmholtz decomposition operator. Because the zero-order operator is the true solution for the case of ϵ = δ, it produces accurate wavefield separation results for elliptical anisotropic media. The first-order separation operator is more accurate for non-elliptical anisotropy. Since the proposed pseudo-Helmholtz decomposition requires solving an anisotropic Poisson's equation, we propose two fast numerical solvers. One is based on the sparse lower-upper (LU) factorization, which can be repeatedly applied to the input elastic wavefields once computing the lower and upper triangular matrices. The second solver assumes the model parameters are laterally homogeneous within a given migration aperture. This assumption allows us to efficiently solve the anisotropic Poisson's equation in the z k x domain, where k x and z denote the horizontal wavenumber and depth, respectively. Using the coordinate transform, we extend the pseudo-Helmholtz decomposition to the TTI media. The separated vector wavefields are used to produce PP and PS images by applying a dot-product imaging condition. Several numerical examples demonstrate the feasibility and applicability of the proposed methods. © The Author(s) 2019. Published by Oxford University Press on behalf of The Royal Astronomical Society.



Seismology--Computer programs, Wave equation--Numerical solutions, Anisotropy, Approximation theory, Eigenvalues, Fourier transformations, Poisson's equation, Seismology, Wave equation, Eigenfunctions, Eigenvectors



©2019 The Authors. Published by Oxford University Press on behalf of The Royal Astronomical Society