Deriving a fast P/S decoupled operator from 3D anisotropic Christoffel equation with its application in elastic reverse time migration and angle-domain common-image gathers
Lele Zhang,Yang Zhao,Lu Liu,Houzhu Zhang,Ke Chen,Jiahui Zuo,Yixin Wang,Fuyu Zhu,Shichang Li
DOI: https://doi.org/10.1109/tgrs.2023.3308162
IF: 8.2
2023-01-01
IEEE Transactions on Geoscience and Remote Sensing
Abstract:A large-scale computation usually hinders the production of the existing 3-D anisotropic P/S wave-mode decomposition methods, such as low-rank approximation, lower-upper (LU) factorization, or local Fourier transformation (FT). To tackle this problem, we develop a fast decoupled operator in 3-D vertical transverse isotropic (VTI) media for P/S wave-mode decomposition and apply this operator to elastic reverse time migration (ERTM) and angle-domain common-image gathers (ADCIGs). To start with, we construct the 3-D VTI Christoffel equation in the wavenumber domain and obtain three elliptical solutions that represent the P-, SV-, and SH-polarization directions. We project the original elastic wavefields to these polarization vectors and derive the decoupled formulations of P-, SV-, and SH-waves. To improve the efficiency, these decoupled formulations are then converted to the space domain and depend on the model parameters and phase angles. Further removing the phase angle term from such space-domain decoupled formulations, the approximate P-, SV-, and SH-wavefield components are finally obtained. These approximate wave components are proven to compute an accurate phase angle as that of exact P-, SV-, and SH-waves and can be used for ERTM imaging and ADCIGs extraction. A fast decoupled operator is thusly developed and only requires gradient operations and once fast Fourier transformation (FFT, for SV-waves). Both simple and complex synthetic examples demonstrate the effectiveness and feasibility of our approaches.
imaging science & photographic technology,remote sensing,engineering, electrical & electronic,geochemistry & geophysics