Quick Search:    advanced search

GSW Home   GeoRef Home   My GSW Alerts   Contact GSW   About GSW   Journals List   Help

This Article Full Text (PDF) References Alert me when this article is cited Alert me if a correction is posted 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 CERVENy, V. Articles by HRON, F. Search for Related Content GeoRef GeoRef Citation

The ray series method and dynamic ray tracing system for three-dimensional inhomogeneous media

V. ERVEN^* and F. HRON

INSTITUTE OF EARTH AND PLANETARY PHYSICS DEPARTMENT OF PHYSICS UNIVERSITY OF ALBERTA, EDMONTON, ALBERTA, Canada, T6G 2J1

Abstract

The ray series method is used to study propagation of seismic waves in the three-dimensional media consisting of generally inhomogeneous layers separated by curved interfaces. The investigation is carried out with the help of the so-called ray centered corrdinate system which was proposed by Popov and Peník (1978a). It is shown that in this system the principal components of the amplitude coefficients in the ray series for S waves do not rotate about the ray with respect to the basis vectors when the wave progresses, even though they rotate with respect to the unit vectors n^ and b^ along the direction of the normal and binormal to the ray, respectively. This considerably simplifies the final expressions for the amplitude coefficients for S waves, whose two principal components are decoupled in the ray-centered coordinate system. The ray-centered coordinate system is also applied to the eikonal equation in order to produce a dynamic ray tracing system consisting of three nonlinear ordinary differential equations of the first order determining the second derivatives of the time field and, in this manner, even the basic geometrical properties of the wave fronts (e.g., principal curvatures and geometrical spreading) along the ray. Several different modifications of the dynamic ray tracing system are presented. It is demonstrated that in the case of generally inhomogeneous two-dimensional media the dynamic ray tracing system reduces, under certain not too restrictive conditions, to the single first order differential equation of the Riccati type. Finally, the phase matching method is used to determine discontinuities of individual quantities in the dynamic ray tracing system when the wave is impinging on a curved interface separating two generally inhomogeneous media. Since all basic equations are presented in a computationally convenient matrix formulation, they can be readily employed for any numberical evaluation of dynamic properties of seismic waves propagating through structurally complicated media. As the paper describes all basic features of asymptotic ray theory (the name under which the ray series method is known on this continent), it can serve as a starting point for anyone wishing to develop computer programs for the computation of synthetic seismograms.

Footnotes

^* On leave from Geophysical Institute, Charles University, Prague, Czechoslovakia.

This article has been cited by other articles:

				W. Kim and C.-E. Baag Rapid and Accurate Two-Point Ray Tracing Based on a Quadratic Equation of Takeoff Angle in Layered Media with Constant or Linearly Varying Velocity Functions Bulletin of the Seismological Society of America, August 1, 2002; 92(6): 2251 - 2263. [Abstract] [Full Text] [PDF]

				F. HRON and B. S. ZHENG On the longitudinal component of the particle motion carried by the shear PS wave reflected from the free surface at normal incidence Bulletin of the Seismological Society of America, October 1, 1993; 83(5): 1610 - 1616. [PDF]

				P. F. DALEY and F. HRON Reflection of an incident spherical P wave on a free surface (near-vertical incidence) Bulletin of the Seismological Society of America, June 1, 1987; 77(3): 1057 - 1070. [Abstract] [PDF]

				V. F. CORMIER and G. C. BEROZA Calculation of strong ground motion due to an extended earthquake source in a laterally varying structure Bulletin of the Seismological Society of America, February 1, 1987; 77(1): 1 - 13. [Abstract] [PDF]

				F. HRON and B. G. MIKHAILENKO Numerical modeling of nongeometrical effects by the Alekseev-Mikhailenko method Bulletin of the Seismological Society of America, August 1, 1981; 71(4): 1011 - 1029. [Abstract] [PDF]

				E. S. KREBES and F. HRON Synthetic seismograms for SH waves in layered anelastic medium by asymptotic ray theory Bulletin of the Seismological Society of America, December 1, 1980; 70(6): 2005 - 2020. [Abstract] [PDF]