HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS		QUICK SEARCH:   [advanced] Author: Keyword(s): Year:  Vol:  Page:

This Article Full Text (PDF) References Citation Map Services Email this article to a friend Similar articles in this journal Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by PEREYRA, V. Articles by KELLER, H. B. Search for Related Content GeoRef GeoRef Citation

Solving two-point seismic-ray tracing problems in a heterogeneous medium

Part 1. A general adaptive finite difference method

DEPARTAMENTO DE COMPUTACION UNIVERSIDAD CENTRAL DE VENEZUELA, CARACAS 104, Venezuela
OFFICE OF EARTHQUAKE STUDIES U.S. GEOLOGICAL SURVEY, MENLO PARK, CALIFORNIA 94025
APPLIED MATHEMATICS DEPARTMENT CALIFORNIA INSTITUTE OF TECHNOLOGY, PASADENA, CALIFORNIA 91125

Abstract

A study of two-point seismic-ray tracing problems in a heterogeneous isotropic medium and how to solve them numerically will be presented in a series of papers. In this Part 1, it is shown how a variety of two-point seismic-ray tracing problems can be formulated mathematically as systems of first-order nonlinear ordinary differential equations subject to nonlinear boundary conditions. A general numerical method to solve such systems in general is presented and a computer program based upon it is described. High accuracy and efficiency are achieved by using variable order finite difference methods on nonuniform meshes which are selected automatically by the program as the computation proceeds. The variable mesh technique adapts itself to the particular problem at hand, producing more detailed computations where they are needed, as in tracing highly curved seismic rays.

A complete package of programs has been produced which use this method to solve two- and three-dimensional ray-tracing problems for continuous or piecewise continuous media, with the velocity of propagation given either analytically or only at a finite number of points. These programs are all based on the same core program, PASVA3, and therefore provide a compact and flexible tool for attacking ray-tracing problems in seismology.

In Part 2 of this work, the numerical method is applied to two- and three-dimensional velocity models, including models with jump discontinuities across interfaces.

This article has been cited by other articles:

				T. Xu, Z. Zhang, E. Gao, G. Xu, and L. Sun Segmentally Iterative Ray Tracing in Complex 2D and 3D Heterogeneous Block Models Bulletin of the Seismological Society of America, April 1, 2010; 100(2): 841 - 850. [Abstract] [Full Text] [PDF]

				S. Ballard, J. R. Hipp, and C. J. Young Efficient and Accurate Calculation of Ray Theory Seismic Travel Time through Variable Resolution 3D Earth Models Seismological Research Letters, November 1, 2009; 80(6): 989 - 999. [Full Text] [PDF]

				T. J. MOSER, G. NOLET, and R. SNIEDER Ray bending revisited Bulletin of the Seismological Society of America, February 1, 1992; 82(1): 259 - 288. [Abstract] [PDF]

				W. A. PROTHERO, W. J. TAYLOR, and J. A. EICKEMEYER A fast, two-point, three-dimensional raytracing algorithm using a simple step search method Bulletin of the Seismological Society of America, June 1, 1988; 78(3): 1190 - 1198. [Abstract] [PDF]

				J. UM and C. THURBER A fast algorithm for two-point seismic ray tracing Bulletin of the Seismological Society of America, June 1, 1987; 77(3): 972 - 986. [Abstract] [PDF]

				R. P. COMER Rapid seismic ray tracing in a spherically symmetric earth via interpolation of rays Bulletin of the Seismological Society of America, April 1, 1984; 74(2): 479 - 492. [Abstract] [PDF]

				C. H. THURBER and W. L. ELLSWORTH Rapid solution of ray tracing problems in heterogeneous media Bulletin of the Seismological Society of America, August 1, 1980; 70(4): 1137 - 1148. [Abstract] [PDF]