High-order Accurate Solution of the TTI Eikonal Equation

Summary High-frequency asymptotic methods, based on solving the eikonal equation, are widely used in many seismic applications including Kirchhoff migration and traveltime tomography. Finite-difference methods to solve the eikonal equation are computationally more efficient and attractive than ray tracing. But, finite-difference solution of the eikonal equation for a point-source contains inaccuracies due to source-singularity. Compared to the several proposed approaches to tackle source-singularity, factorization of the unknown traveltime is computationally efficient and simpler to implement. Recently, a factorization algorithm has been proposed to obtain clean first-order accuracy for tilted transversely isotropic (TTI) media. However, high-order accuracy of traveltimes is needed for computation of quantities that require traveltime derivatives, such as take-off angle and amplitudes. I propose an iterative fast sweeping algorithm to obtain high-order accuracy using factorization and a high-order finite-difference stencil. Numerical test shows improvements in accuracy of the TTI eikonal solution. This shows that once the source-singularity problem is tackled, high-order accurate solutions can be constructed easily. The method can be easily extended to media with lower anisotropic symmetries.