Pseudo-bending method for three-dimensional seismic ray tracing in a spherical earth with discontinuities

We present a 'pseudo-bending' approach to 3-D ray tracing in a spherical earth with discontinuities. This method is based on a three-point perturbation associated with a first-order approximation, while Snell's law in curvilinear coordinates is applied at the discontinuities. We demonstrate the computational accuracy and efficiency of the pseudo-bending method in tracing rays for various velocity models by comparing results with analytical solutions and with results from the bending method. The improvement of efficiency is significant, but is reduced as the number of discontinuities increases. Since the bending approach may be computationally unstable in some situations, even though it is exact, the pseudo-bending approach is preferable for automatic calculation of rays.