3D Anisotropic P- and S-Mode Wavefields Separation in 3D Elastic Reverse-Time Migration
Jiahui Zuo,Fenglin Niu,Lu Liu,Shuai Da,Houzhu Zhang,Jidong Yang,Lele Zhang,Yang Zhao
DOI: https://doi.org/10.1007/s10712-021-09688-8
2022-01-01
Surveys in Geophysics
Abstract:In 2D anisotropic media, non-stationary filters and low-rank approximation methods are classical strategies to compute the decomposition operators, but they suffer from expensive computational costs for 3D media. This study adopts the eigenform analysis into 3D vertical transverse isotropic (VTI) media and produces the separated vector P and S wavefields with the same amplitudes, phases, and physical units as the input elastic wavefields. We first built a 3D zero-order pseudo-Helmholtz decomposition operator by deriving the eigenvalues and eigenvectors of the 3D VTI wave equations in the wavenumber domain. The eigenvalues refer to the phase velocities of P-, SH-, and SV-wave, and the corresponding eigenvectors are pointing to their polarizations. Second, we use the pseudo-Helmholtz decomposition operator to construct a 3D anisotropic Poisson’s equation. Based on the laterally homogeneous assumption, Poisson’s equation is solved in the mixed domain z - k_x - k_y , where k_x , k_y , and z denote the horizontal wavenumbers and depth, respectively. Third, we obtain the vector P and S wavefields using the proposed 3D pseudo-Helmholtz decomposition operator in the space domain. Lastly, 3D PP and PS images are calculated with a dot-product imaging condition. The anisotropic amplitude versus offset (AVO) responses of the 3D elastic reverse-time migration (ERTM) images are also validated by analytical solutions (Ruger’s equations). Our proposed 3D pseudo-Helmholtz decomposition operator degrades to a gradient operator satisfying isotropic media conditions. In addition, the method is easy to extend into 3D due to its high-efficiency cost. Several numerical examples with large shear anisotropy are selected to demonstrate the feasibility of our proposed pseudo-Helmholtz decomposition method in 3D applications.