Approximations for traveltime, slope, curvature, and geometrical spreading of elastic waves in layered transversely isotropic media

Each seismic body wave, including quasi compressional, shear, and converted wave modes, carries useful subsurface information. For processing, imaging, amplitude analysis, and forward modeling of each wave mode, we need approximate equations of traveltime, slope (ray-parameter), and curvature as a function of offset. Considering the large offset coverage of modern seismic acquisitions, we propose new approximations designed to be accurate at zero and infinitely large offsets over layered transversely isotropic media with vertical symmetry axis (VTI). The proposed approximation for traveltime is a modified version of the extended generalized moveout approximation that comprises six parameters. The proposed direct approximations for ray-parameter and curvature use new, algebraically simple, equations with three parameters. We define these parameters for each wave mode without ray tracing so that we have similar approximate equations for all wave modes that only change based on the parameter definitions. However, our approximations are unable to reproduce S-wave triplications that may occur in some strongly anisotropic models. Using our direct approximation of traveltime derivatives, we also obtain a new expression for the relative geometrical spreading. We demonstrate the high accuracy of our approximations using numerical tests on a set of randomly generated multilayer models. Using synthetic data, we present simple applications of our approximations for normal moveout correction and relative geometrical spreading compensation of different wave modes.

Homotopy solutions of the acoustic eikonal equation for strongly attenuating transversely isotropic media with a vertical symmetry axis

Many applications in seismology involve the modeling of seismic wave traveltimes in anisotropic media. We present homotopy solutions of the acoustic eikonal equation for P-waves traveltimes in attenuating transversely isotropic media with a vertical symmetry axis. Instead of the commonly used perturbation theory, we use the homotopy analysis method to express the traveltimes by a Taylor series expansion over an embedding parameter. For the derivation, we first perform homotopy analysis of the eikonal equation and derive the linearized ordinary differential equations for the coefficients of the Taylor series expansion. Then, we obtain the homotopy solutions for the traveltimes by solving the linearized ordinary differential equations. Results on approximate formulae investigations demonstrate that the analytical expressions are efficient methods for the computation of traveltimes from the eikonal equation. In addition, these formulas are also effective methods for benchmarking approximated solutions in strongly attenuating anisotropic media.

Large-offset P-wave traveltime in layered transversely isotropic media

Large-offset seismic data processing, imaging, and velocity estimation require an accurate traveltime approximation over a wide range of offsets. In layered transversely isotropic media with vertical symmetry axis (VTI), the accuracy of traditional traveltime approximations is limited to near offsets. Herein, we propose a new traveltime approximation that maintains the accuracy of the classical equations around zero offset, and exhibits the correct curvilinear asymptote at infinitely large offsets. Our approximation is based on the conventional acoustic assumption. Its equation incorporates six parameters. To define them, we use the Taylor series expansion of the exact traveltime around zero offset, and a new asymptotic series for infinite offset. Our asymptotic equation shows that the traveltime behavior at infinitely large offsets is dominated by the properties of the layer with the maximum horizontal velocity in the sequence. The parameters of our approximation depend on: the effective zero offset traveltime, the normal moveout velocity, the anellipticity, a new large-offset heterogeneity parameter, and the properties of the layer with the maximum horizontal velocity in the sequence. We apply our traveltime approximation: (1) to directly calculate traveltime and ray parameter at given offsets, as analytical forward modeling; and (2) to estimate the first four of the aforementioned parameters for the layers beneath a known high-velocity layer. Our large-offset heterogeneity parameter includes the layering effect on the reflections traveltime.