Astronomy and Astrophysics – Astronomy
Scientific paper
May 1999
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=1999geoji.137..364t&link_type=abstract
Geophysical Journal International, Volume 137, Issue 2, pp. 364-380.
Astronomy and Astrophysics
Astronomy
12
Anisotropy, Body Waves, Finite Difference Methods, One-Way Wave Equations, Ray Theory
Scientific paper
A derivation is given for the factorization of the generally anisotropic, elastic wave equation into two operators which are first order with respect to a preferred direction of propagation. One factor gives the `one-way' equation for forward-going body waves. A 3 x 3 matrix notation is used which emphasizes the role of the Christoffel equation (or what might be called its `preferred coordinate projection') familiar from ray theory. To recognize this it helps to obtain first an exact factorization of the wave equation for homogeneous anisotropic media. When inhomogeneities exist it is necessary to use a Fourier representation for differential operators with respect to the transverse coordinates, but still the Christoffel equation is the key. The so-called `square root operator' required in the factorization becomes more correctly the `root of a matrix operator ratic equation' intimately connected to the algebraic Christoffel equation. The factorization does not assume narrow-angle propagation or involve an a priori reference phase. It includes the P and two S waves and the forward coupling between them. It remains valid at slowness-surface singularities such as conical points (also known as acoustic axes-a most stringent test) and for wave fronts which fold. In the frequency domain this wide-angle one-way equation involves Fourier synthesis with respect to transverse position/slowness and its time-domain analogue involves2-D slant stacking (Radon transformation), multiplication by a 3 x 3 `propagator' matrix and an inverse 2-D slant stack. The one-way equation can be solved `analytically' by using `path integrals'. These in turn may be reduced by stationary-phase arguments to standard ray theory, Maslov and Kirchhoff representations. In this sense, the `gap' between ray theory and the numerical solution of the full wave equation is bridged or filled in, at least for forward propagation. Alternatively, the one-way wave equation can be solved by numerical `forward-stepping' algorithms and here again a number of choices are available. Narrow-angle Taylor or rational approximations to the matrix propagator in the Fourier/slant-stack representation lead to the anisotropic analogues of the classic 15 deg and 45 deg partial differential equations. The integral wide-angle form is also reminiscent of the `split-step'and `elastic phase screen' approaches of Tappert, Wu and others, suggesting other implementations.
No associations
LandOfFree
The `gap' between seismic ray theory and `full' wavefield extrapolation does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with The `gap' between seismic ray theory and `full' wavefield extrapolation, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and The `gap' between seismic ray theory and `full' wavefield extrapolation will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-1125092