Preliminary investigations on computational methods for solving the two-point seismic ray-tracing problem in a heterogeneous and isotropic medium

Open-File Report 76-707
By:  and 

Links

Abstract

Two-point seismic ray tracing is an important problem in seismology. In this report, we study this problem by first deriving the differential equations governing seismic wave propagation between two end points in a 3-dimensional heterogeneous and isotropic medium. We then formulate three numerical schemes to solve the two-point ray-tracing problem. In the first method, introduced originally by R. L. Wesson, the second order ray equations are reduced by the central difference approximation to a set of simultaneous non-linear equations, which are then solved by Newton's method. In the second method, the second order ray equations are reduced to a set of first order equations which are solved by using an adaptive finite difference method introduced by M. Lentini and V. Pereyra. The third method, due originally to L. Euler uses a sum to approximate the integral for the travel time which are solved directly for the minimum time path. This formulation may be shown to be equivalent to Wesson's method.

We programed Wesson's method for a 3-dimensional inhomogeneous velocity medium, and also applied Lentini and Pereyra's method for a 2-dimensional linear velocity model. Our programs are tested for constant velocity, and linear velocity in one, and two, space variable models. Both programs give satisfactory answers in comparison with known analytic solutions, and with respect to each other.

Publication type Report
Publication Subtype USGS Numbered Series
Title Preliminary investigations on computational methods for solving the two-point seismic ray-tracing problem in a heterogeneous and isotropic medium
Series title Open-File Report
Series number 76-707
DOI 10.3133/ofr76707
Year Published 1976
Language English
Publisher U.S. Geological Survey
Description 66 p.
Google Analytic Metrics Metrics page
Additional publication details