Gaussian beams in inhomogeneous media: A review

The paper outlines the most important results of the paraxial complex geometrical optics (CGO) in respect to Gaussian beams diffraction in the smooth inhomogeneous media and discusses interrelations between CGO and other asymptotic methods, which reduce the problem of Gaussian beam diffraction to the solution of ordinary differential equations, namely: (i) Babich’s method, which deals with the abridged parabolic equation and describes diffraction of the Gaussian beams; (ii) complex form of the dynamic ray tracing method, which generalizes paraxial ray approximation on Gaussian beams and (iii) paraxial WKB approximation by Pereverzev, which gives the results, quite close to those of Babich’s method. For Gaussian beams all the methods under consideration lead to the similar ordinary differential equations, which are complex-valued nonlinear Riccati equation and related system of complex-valued linear equations of paraxial ray approximation. It is pointed out that Babich’s method provides diffraction substantiation both for the paraxial CGO and for complex-valued dynamic ray tracing method. It is emphasized also that the latter two methods are conceptually equivalent to each other, operate with the equivalent equations and in fact are twins, though they differ by names. The paper illustrates abilities of the paraxial CGO method by two available analytical solutions: Gaussian beam diffraction in the homogeneous and in the lens-like media, and by the numerical example: Gaussian beam reflection from a plane-layered medium.     

Comparison of the anisotropic-ray-theory rays and anisotropic common S-wave rays with the SH and SV reference rays in a velocity model with a split intersection singularity

We describe the behaviour of the anisotropic–ray–theory S–wave rays in a velocity model with a split intersection singularity. The anisotropic–ray–theory S–wave rays crossing the split intersection singularity are smoothly but very sharply bent. While the initial–value rays can be safely traced by solving Hamilton’s equations of rays, it is often impossible to determine the coefficients of the equations of geodesic deviation (paraxial ray equations, dynamic ray tracing equations) and to solve them numerically. As a result, we often know neither the matrix of geometrical spreading, nor the phase shift due to caustics. We demonstrate the abrupt changes of the geometrical spreading and wavefront curvature of the fast anisotropic–ray–theory S wave. We also demonstrate the formation of caustics and wavefront triplication of the slow anisotropic–ray–theory S wave.Since the actual S waves propagate approximately along the SH and SV reference rays in this velocity model, we compare the anisotropic–ray–theory S–wave rays with the SH and SV reference rays. Since the coupling ray theory is usually calculated along the anisotropic common S–wave rays, we also compare the anisotropic common S–wave rays with the SH and SV reference rays.     

Spatial derivatives and perturbation derivatives of amplitude in isotropic and anisotropic media

Explicit equations for the spatial derivatives and perturbation derivatives of amplitude in both isotropic and anisotropic media are derived. The spatial and perturbation derivatives of the logarithm of amplitude can be calculated by numerical quadratures along the rays. The spatial derivatives of amplitude may be useful in calculating the higher-order terms in the ray series, in calculating the higher-order amplitude coefficients of Gaussian beams, in estimating the accuracy of zero-order approximations of both the ray method and Gaussian beams, in estimating the accuracy of the paraxial approximation of individual Gaussian beams, or in estimating the accuracy of the asymptotic summation of paraxial Gaussian beams. The perturbation derivatives of amplitude may be useful in perturbation expansions from elastic to viscoelastic media and in estimating the accuracy of the common-ray approximations of the amplitude in the coupling ray theory.     

Paraxial ray methods in inhomogeneous anisotropic media: Initial conditions

Paraxial ray methods have found broad applications in the seismic ray method and in numerical modelling and interpretation of high-frequency seismic wave fields propagating in inhomogeneous, isotropic or anisotropic structures. The basic procedure in paraxial ray methods consists in dynamic ray tracing. We derive the initial conditions for dynamic ray equations in Cartesian coordinates, for rays initiated at three types of initial manifolds given in a three-dimensional medium: 1) curved surfaces (surface source), 2) isolated points (point source), and 3) curved, planar and non-planar lines (line source). These initial conditions are very general, valid for homogeneous or inhomogeneous, isotropic or anisotropic media, and for both a constant and a variable initial travel time along the initial manifold. The results presented in the paper considerably extend the possible applications of the paraxial ray method.     

A note on dynamic ray tracing in ray-centered coordinates in anisotropic inhomogeneous media

Dynamic ray tracing plays an important role in paraxial ray methods. In this paper, dynamic ray tracing systems for inhomogeneous anisotropic media, consisting of four linear ordinary differential equations of the first order along the reference ray, are studied. The main attention is devoted to systems expressed in a particularly simple choice of ray-centered coordinates, here referred to as the standard ray-centered coordinates, and in wavefront orthonormal coordinates. These two systems, known from the literature, were derived independently and were given in different forms. In this paper it is proved that both systems are fully equivalent. Consequently, the dynamic ray tracing system, consisting of four equations in wavefront orthonormal coordinates, can also be used if we work in ray-centered coordinates, and vice versa. vcerveny@seis.karlov.mff.cuni.cz     

Paraxial ray methods for anisotropic inhomogeneous media

A new formalism of surface-to-surface paraxial matrices allows a very general and flexible formulation of the paraxial ray theory, equally valid in anisotropic and isotropic inhomogeneous layered media. The formalism is based on conventional dynamic ray tracing in Cartesian coordinates along a reference ray. At any user-selected pair of points of the reference ray, a pair of surfaces may be defined. These surfaces may be arbitrarily curved and oriented, and may represent structural interfaces, data recording surfaces, or merely formal surfaces. A newly obtained factorization of the interface propagator matrix allows to transform the conventional 6 × 6 propagator matrix in Cartesian coordinates into a 6 × 6 surface-to-surface paraxial matrix. This matrix defines the transformation of paraxial ray quantities from one surface to another. The redundant non-eikonal and ray-tangent solutions of the dynamic ray-tracing system in Cartesian coordinates can be easily eliminated from the 6 × 6 surface-to-surface paraxial matrix, and it can be reduced to 4 × 4 form. Both the 6 × 6 and 4 × 4 surface-to-surface paraxial matrices satisfy useful properties, particularly the symplecticity. In their 4 × 4 reduced form, they can be used to solve important boundary-value problems of a four-parametric system of paraxial rays, connecting the two surfaces, similarly as the well-known surface-to-surface matrices in isotropic media in ray-centred coordinates. Applications of such boundary-value problems include the two-point eikonal, relative geometrical spreading, Fresnel zones, the design of migration operators, and more.     

A note on two-point paraxial travel times

Recently, several expressions for the two-point paraxial travel time in laterally varying, isotropic or anisotropic layered media were derived. The two-point paraxial travel time gives the travel time from point S′ to point R′, both these points being situated close to a known reference ray Ω, along which the ray-propagator matrix was calculated by dynamic ray tracing. The reference ray and the position of points S′ and R′ are specified in Cartesian coordinates. Two such expressions for the two-point paraxial travel time play an important role. The first is based on the 4 × 4 ray propagator matrix, computed by dynamic ray tracing along the reference ray in ray-centred coordinates. The second requires the knowledge of the 6 × 6 ray propagator matrix computed by dynamic ray tracing along the reference ray in Cartesian coordinates. Both expressions were derived fully independently, using different methods, and are expressed in quite different forms. In this paper we prove that the two expressions are fully equivalent and can be transformed into each other.     

PARABOLIC AND HYPERBOLIC PARAXIAL TWO-POINT TRAVELTIMES IN 3D MEDIA1

The 4 × 4 T -propagator matrix of a 3D central ray determines, among other important seismic quantities, second-order (parabolic or hyperbolic) two-point traveltime approximations of certain paraxial rays in the vicinity of the known central ray through a 3D medium consisting of inhomogeneous isotropic velocity layers. These rays result from perturbing the start and endpoints of the central ray on smoothly curved anterior and posterior surfaces. The perturbation of each ray endpoint is described only by a two-component vector. Here, we provide parabolic and hyperbolic paraxial two-point traveltime approximations using the T -propagator to feature a number of useful 3D seismic models, putting particular emphasis on expressing the traveltimes for paraxial primary reflected rays in terms of hyperbolic approximations. These are of use in solving several forward and inverse seismic problems. Our results simplify those in which the perturbation of the ray endpoints upon a curved interface is described by a three-component vector. In order to emphasize the importance of the hyperbolic expression, we show that the hyperbolic paraxial-ray traveltime (in terms of four independent variables) is exact for the case of a primary ray reflected from a planar dipping interface below a homogeneous velocity medium.     

Common-ray tracing and dynamic ray tracing for S waves in a smooth elastic anisotropic medium

Anisotropic common S-wave rays are traced using the averaged Hamiltonian of both S-wave polarizations. They represent very practical reference rays for calculating S waves by means of the coupling ray theory. They eliminate problems with anisotropic-ray-theory ray tracing through some S-wave slowness-surface singularities and also considerably simplify the numerical algorithm of the coupling ray theory for S waves. The equations required for anisotropic-common-ray tracing for S waves in a smooth elastic anisotropic medium, and for corresponding dynamic ray tracing in Cartesian or ray-centred coordinates, are presented. The equations, for the most part generally known, are summarized in a form which represents a complete algorithm suitable for coding and numerical applications.     

3D multivalued travel time and amplitude maps

An algorithm for computing multivalued maps for travel time, amplitude and any other ray related variable in 3D smooth velocity models is presented. It is based on the construction of successive isochrons by tracing a uniformly dense discrete set of rays by fixed travel-time steps. Ray tracing is based on Hamiltonian formulation and includes computation of paraxial matrices. A ray density criterion ensures uniform ray density along isochrons over the entire ray field including caustics. Applications to complex models are shown.     

Comparison of Ray Methods with the Exact Solution in the 1-D Anisotropic “Simplified Twisted Crystal” Model

The exact analytical solution for the plane S-wave, propagating along the axis of spirality in the simple 1-D anisotropic simplified twisted crystal model, is compared with four different approximate ray-theory solutions. The four different ray methods are (a) the coupling ray theory, (b) the coupling ray theory with the quasi-isotropic perturbation of travel times, (c) the anisotropic ray theory, (d) the isotropic ray theory. The comparison is carried out numerically, by evaluating both the exact analytical solution and the analytical solutions of the equations of the four ray methods. The comparison simultaneously demonstrates the limits of applicability of the isotropic and anisotropic ray theories, and the superior accuracy of the coupling ray theory over a broad frequency range. The comparison also shows the possible inaccuracy due to the quasi-isotropic perturbation of travel times in the equations of the coupling ray theory. The coupling ray theory thus should definitely be preferred to the isotropic and anisotropic ray theories, but the quasi-isotropic perturbation of travel times should be avoided. Although the simplified twisted crystal model is designed for testing purposes and has no direct relation to geological structures, the wave-propagation phenomena important in the comparison are similar to those in the models of the geological structures.In additional numerical tests, the exact analytical solution is numerically compared with the finite-difference numerical results, and the analytical solutions of the equations of different ray methods are compared with the corresponding numerical results of 3-D ray-tracing programs developed by the authors of the paper.     

Errors Due to the Common Ray Approximations of the Coupling Ray Theory

The common ray approximation considerably simplifies the numerical algorithm of the coupling ray theory for S waves, but may introduce errors in travel times due to the perturbation from the common reference ray. These travel-time errors can deteriorate the coupling-ray-theory solution at high frequencies. It is thus of principal importance for numerical applications to estimate the errors due to the common ray approximation.We derive the equations for estimating the travel-time errors due to the isotropic and anisotropic common ray approximations of the coupling ray theory. These equations represent the main result of the paper. The derivation is based on the general equations for the second-order perturbations of travel time. The accuracy of the anisotropic common ray approximation can be studied along the isotropic common rays, without tracing the anisotropic common rays.The derived equations are numerically tested in three 1-D models of differing degree of anisotropy. The first-order and second-order perturbation expansions of travel time from the isotropic common rays to anisotropic-ray-theory rays are compared with the anisotropic-ray-theory travel times. The errors due to the isotropic common ray approximation and due to the anisotropic common ray approximation are estimated. In the numerical example, the errors of the anisotropic common ray approximation are considerably smaller than the errors of the isotropic common ray approximation.The effect of the isotropic common ray approximation on the coupling-ray-theory synthetic seismograms is demonstrated graphically. For comparison, the effects of the quasi-isotropic projection of the Green tensor, of the quasi-isotropic approximation of the Christoffel matrix, and of the quasi-isotropic perturbation of travel times on the coupling-ray-theory synthetic seismograms are also shown. The projection of the travel-time errors on the relative errors of the time-harmonic Green tensor is briefly presented.     

Errors due to the anisotropic-common-ray approximation of the coupling ray theory

The common-ray approximation eliminates problems with ray tracing through S-wave singularities and also considerably simplifies the numerical algorithm of the coupling ray theory for S waves, but may introduce errors in travel times due to the perturbation from the common reference ray. These travel-time errors can deteriorate the coupling-ray-theory solution at high frequencies. It is thus of principal importance for numerical applications to estimate the errors due to the common-ray approximation applied. The anisotropic-common-ray approximation of the coupling ray theory is more accurate than the isotropic-common-ray approximation. We derive the equations for estimating the travel-time errors due to the anisotropic-common-ray (and also isotropic-common-ray) approximation of the coupling ray theory. The errors of the common-ray approximations are calculated along the anisotropic common rays in smooth velocity models without interfaces. The derivation is based on the general equations for the second-order perturbations of travel time.     

Calculation of Synthetic Seismograms with Gaussian Beams

— In this paper, an overview of the calculation of synthetic seismograms using the Gaussian beam method is presented accompanied by some representative applications and new extensions of the method. Since caustics are a frequent occurrence in seismic wave propagation, modifications to ray theory are often necessary. In the Gaussian beam method, a summation of paraxial Gaussian beams is used to describe the propagation of high-frequency wave fields in smoothly varying inhomogeneous media. Since the beam components are always nonsingular, the method provides stable results over a range of beam parameters. The method has been shown, however, to perform better for some problems when different combinations of beam parameters are used. Nonetheless, with a better understanding of the method as well as new extensions, the summation of Gaussian beams will continue to be a useful tool for the modeling of high-frequency seismic waves in heterogeneous media.     

线性连续变化非弹性介质中地震波射线轨迹与走时研究

地球内部非弹性介质广泛存在,垂向线性连续变化非弹性介质模型具有重要理论与实际意义,本文利用复速度概念推导得到了复速度随深度线性连续变化的非弹性介质中地震波射线轨迹与等时线方程.理论结果表明:在线性连续变化非弹性介质中,射线轨迹呈圆形,等时面的形状为椭圆.最后对射线轨迹与等时线方程的适用性条件进行了讨论.     

Eigenrays in 3D heterogeneous anisotropic media,Part II: Dynamics

This paper is the second in a sequel of two papers and dedicated to the computation of paraxial rays and dynamic characteristics along the stationary rays obtained in the first paper. We start by formulating the linear, second‐order, Jacobi dynamic ray tracing equation. We then apply a similar finite‐element solver, as used for the kinematic ray tracing, to compute the dynamic characteristics between the source and any point along the ray. The dynamic characteristics in our study include the relative geometric spreading and the phase correction due to caustics (i.e. the amplitude and the phase of the asymptotic form of the Green's function for waves propagating in 3D heterogeneous general anisotropic elastic media). The basic solution of the Jacobi equation is a shift vector of a paraxial ray in the plane normal to the ray direction at each point along the central ray. A general paraxial ray is defined by a linear combination of up to four basic vector solutions, each corresponds to specific initial conditions related to the ray coordinates at the source. We define the four basic solutions with two pairs of initial condition sets: point–source and plane‐wave. For the proposed point–source ray coordinates and initial conditions, we derive the ray Jacobian and relate it to the relative geometric spreading for general anisotropy. Finally, we introduce a new dynamic parameter, similar to the endpoint complexity factor, presented in the first paper, used to define the measure of complexity of the propagated wave/ray phenomena. The new weighted propagation complexity accounts for the normalized relative geometric spreading not only at the receiver point, but along the whole stationary ray path. We propose a criterion based on this parameter as a qualifying factor associated with the given ray solution. To demonstrate the implementation of the proposed method, we use several isotropic and anisotropic benchmark models. For all the examples, we first compute the stationary ray paths, and then compute the geometric spreading and analyse these trajectories for possible caustics. Our primary aim is to emphasize the advantages, transparency and simplicity of the proposed approach.     

Nonlinear ray perturbation theory with its applications to ray tracing and inversion in anisotropic media

A comprehensive approach, based on the general nonlinear ray perturbation theory (Druzhinin, 1991), is proposed for both a fast and accurate uniform asymptotic solution of forward and inverse kinematic problems in anisotropic media. It has been developed to modify the standard ray linearization procedures when they become inconsistent, by providing a predictable truncation error of ray perturbation series. The theoretical background consists in a set of recurrent expressions for the perturbations of all orders for calculating approximately the body wave phase and group velocities, polarization, travel times, ray trajectories, paraxial rays and also the slowness vectors or reflected/transmitted waves in terms of elastic tensor perturbations. We assume that any elastic medium can be used as an unperturbed medium. A total 2-D numerical testing of these expressions has been established within the transverse isotropy to verify the accuracy and convergence of perturbation series when the elastic constants are perturbed. Seismological applications to determine crack-induced anisotropy parameters on VSP travel times for the different wave types in homogeneous and horizontally layered, transversally isotropic and orthorhombic structures are also presented. A number of numerical tests shows that this method is in general stable with respect to the choice of the reference model and the errors in the input data. A proof of uniqueness is provided by an interactive analysis of the sensitivity functions, which are also used for choosing optimum source/receiver locations. Finally, software has been developed for a desktop computer and applied to interpreting specific real VSP observations as well as explaining the results of physical modelling for a 3-D crack model with the estimation of crack parameters.     

Ray tracing and geodesic deviation of the SH and SV reference rays in a heterogeneous generally anisotropic medium which is approximately uniaxial

The coupling ray theory is usually applied to anisotropic common reference rays, but it is more accurate if it is applied to reference rays which are closer to the actual wave paths. If we know that a medium is close to uniaxial (transversely isotropic), it may be advantageous to trace reference rays which resemble the SH–wave and SV–wave rays. This paper is devoted to defining and tracing these SH and SV reference rays of elastic S waves in a heterogeneous generally anisotropic medium which is approximately uniaxial (approximately transversely isotropic), and to the corresponding equations of geodesic deviation (dynamic ray tracing). All presented equations are simultaneously applicable to ordinary and extraordinary reference rays of electromagnetic waves in a generally bianisotropic medium which is approximately uniaxially anisotropic. The improvement of the coupling–ray–theory seismograms calculated along the proposed SH and SV reference rays, compared to the coupling–ray–theory seismograms calculated along the anisotropic common reference rays, has already been numerically demonstrated by the authors in four approximately uniaxial velocity models.     

Transformation relations for second-order derivatives of travel time in anisotropic media

In the computation of paraxial travel times and Gaussian beams, the basic role is played by the second-order derivatives of the travel-time field at the reference ray. These derivatives can be determined by dynamic ray tracing (DRT) along the ray. Two basic DRT systems have been broadly used in applications: the DRT system in Cartesian coordinates and the DRT system in ray-centred coordinates. In this paper, the transformation relations between the second-order derivatives of the travel-time field in Cartesian and ray-centred coordinates are derived. These transformation relations can be used both in isotropic and anisotropic media, including computations of complex-valued travel times necessary for the evaluation of Gaussian beams.     

耗散大气中的声波射线追踪

基于分层大气中声波的局地色散关系方程,建立一种计入真实大气衰减效应的有耗大气声波射线追踪模型.在色散方程的虚部中导出声波在运动大气中的耗散系数和竖直方向上的增长因子,并利用真实大气中的衰减理论对所得到的耗散系数进行修正.利用Hamilton方程组解出大气声波在考虑耗散效应下的射线微分方程组.该有耗射线追踪模型的数值模拟...