Geometrical theory of diffraction, evanescent waves, complex rays and Gaussian beams

Abstract

Summary. The Gaussian beam method has recently been introduced into synthetic seismology to overcome shortcomings of the ray method, especially in transition regions due to focusing or diffraction where ray theory fails. One proceeds by discretizing the initial data as a superposition of paraxial Gaussian beams, each of which is then traced through the seismic environment. Since Gaussian beam fields do not diverge in ray transition regions, they are ‘uniformly regular’ although the quality of this regularity depends on the beam parameters and on the ‘numerical distance’ which defines the extent of the transitional domain. However, when Gaussian beam patches are used to simulate non-Gaussian initial data, there arise ambiguities due to choice of patch size and location, beam width, etc., which are at the user's disposal. The effects of this arbitrariness have customarily been explored by trial and error numerical experiment but no quantitative recommendations have emerged as yet. As a step toward a priori predictive capability, it is proposed here to perform a systematic study on analytically tractable prototype models of how the parameters and location of a single beam affect the quality of the observed seismic field, especially in ray transition regions. The conversion of ordinary ray fields into beam fields in canonical configurations can be accomplished conveniently by displacing a real source point into a complex coordinate space. Thus, the desired beam solutions can be obtained directly from available ray, and even paraxial ray, fields. Complex ray theory and its implications are reviewed here, with an emphasis on improvements of beam tracking schemes employed at present.