We present a new ray tracing software with adaptive step size control as an efficient tool for modeling wave propagation in inhomogeneous anisotropic media. After a short review of the theoretical and numerical basics of ray tracing algorithms, we apply initial-value and two-point ray tracing to two strongly anisotropic velocity models. We compare the efficiency and accuracy of integrating ray tracing systems numerically with three different Runge-Kutta methods. These are the classical Runge-Kutta method with constant step size and two embedded Runge-Kutta methods, the Cash-Karp method and the Dormand-Prince method, with adaptive step size control. The results show that ray tracing with adaptive step size control is significantly more efficient than standard ray tracing algorithms with constant step size. Additionally, our software is more user-friendly because adaptive step size control provides relief from the necessity of having a priori knowledge about the spatial and temporal behavior of ray propagation and the maximum allowed step size in a given medium.


Article metrics loading...

Loading full text...

Full text loading...


  1. Butcher, J.C.
    [2016] Numerical Methods for Ordinary Differential Equations. John Wiley and Sons Ltd, 3rd edn.
    [Google Scholar]
  2. Cash, J.R. and Karp, A.H.
    [1990] A variable order Runge-Kutta method for initial value problems with rapidly varying right-hand sides. ACM Transactions on Mathematical Software, 16(3), 201–222.
    [Google Scholar]
  3. Červený, V.
    [2001] Seismic Ray Theory. Cambridge University Press.
    [Google Scholar]
  4. Dormand, J.R. and Prince, P.J.
    [1980] A family of embedded Runge-Kutta formulae. Journal of Computational and Applied Mathematics, 6(1), 19–26.
    [Google Scholar]
  5. Gajewski, D. and Pšenčǭk, I.
    [1987] Computation of high-frequency seismic wavefields in 3-D laterally inhomogeneous anisotropic media. Geophysical Journal International, 91(2), 383–411.
    [Google Scholar]
  6. Oelke, A., Krüger, O.S., and Shapiro, S.A.
    [2014] Modeling Wave Propagation of Microseismic Events in Anisotropic Media. Extendend abstract, 76th EAGE Conference & Exhibition.
    [Google Scholar]
  7. Reshetnikov, A. and Shapiro, S.A.
    [2015] Using multipath shear wave arrivals for enhanced microseis-mic location in strongly anisotropic shale. Extendend abstract, 77th EAGE Conference & Exhibition.
    [Google Scholar]
  8. Roeser, A. and Shapiro, S.A.
    [2016] 3D Anisotropic Full Waveform Modeling with an Enhanced OASES Workflow for Complex Source-receiver Geometries. Extendend abstract, 78th EAGE Conference & Exhibition.
    [Google Scholar]
  9. Vavryčuk, V.
    [2001] Ray tracing in anisotropic media with singularities. Geophysical Journal International, 145(1), 265–276.
    [Google Scholar]
  10. Virieux, J. and Farra, V.
    [1991] Ray tracing in 3D complex isotropic media: An analysis of the problem. Geophysics, 56(12), 2057–2069.
    [Google Scholar]

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error