Multi-phase arrival tracking using tetrahedral cells within a 3D layered titled transversely isotropic anisotropic model involving undulating topography and irregular interfaces

To overcome a major problem in current ray tracing methods, which are only capable of tracing first arrivals, and occasionally primary reflections (or mode conversions) in regular cell models, we extend in this paper the multistage triangular shortest-path method (SPM) into 3D titled transversely isotropic (TTI) anisotropic media. The proposed method is capable of tracking multi-phase arrivals composed of any kind of combinations of transmissions, mode conversions and reflections. In model parameterization, five elastic parameters, plus two angles defining the titled axis of symmetry of TTI media are sampled at the primary nodes of the tetrahedral cell, and velocity value at secondary node positions are linked by a tri-linear velocity interpolation function to the primary node velocity value of that of a tetrahedral cell, from which the group velocities of the three wave modes (qP, qSV and qSH) are computed. The multistage triangular SPM is used to track multi-phase arrivals. The uniform anisotropic test indicates that the numerical solution agrees well with the analytic solution, thus verifying the accuracy of the methodology. Several simulations and comparison results for heterogeneous models show that the proposed algorithm is able to efficiently and accurately approximate undulating surface topography and irregular subsurface velocity discontinuities. It is suitable for any combination of multi-phase arrival tracking in arbitrary tilt angle TTI media and can accommodate any magnitude of anisotropy.