Traveltime approximation for strongly anisotropic media using the homotopy analysis method

Traveltime approximation plays an important role in seismic data processing, for example, anisotropic parameter estimation and seismic imaging. By exploiting seismic traveltimes, it is possible to improve the accuracy of anisotropic parameter estimation and the resolution of seismic imaging. Conventionally, the traveltime approximations in anisotropic media are obtained by expanding the anisotropic eikonal equation in terms of the anisotropic parameters and the elliptically anisotropic eikonal equation based on perturbation theory. Such an expansion assumes a small perturbation and weak anisotropy. In a realistic medium, however, the assumption of small perturbation likely breaks down. We present a retrieved zero-order deformation equation that creates a map from the anisotropic eikonal equation to a linearized partial differential equation system based on the homotopy analysis method. By choosing the linear and nonlinear operators in the retrieved zero-order deformation equation, we develop new traveltime approximations that allow us to compute the traveltimes for a medium of arbitrarily strength anisotropy. A comparison of the traveltimes and their errors from the homotopy analysis method and from the perturbation method suggests that the traveltime approximations provide a more reliable result in strongly anisotropic media.