Ray equation
 | |
| Series | Geophysical References Series |
|---|---|
| Title | Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing |
| Author | Enders A. Robinson and Sven Treitel |
| Chapter | 2 |
| DOI | http://dx.doi.org/10.1190/1.9781560801610 |
| ISBN | 9781560801481 |
| Store | SEG Online Store |
We need to understand how a single ray moving along a particular path - say, a seismic ray - can know what is an extremal path in the variational sense. To illustrate the problem, consider a seismic wave traveling through a medium whose slowness n increases in the direction of travel. Because the raypath is parallel to the gradient of the slowness, it undergoes no bending. However, if the contour lines of slowness are at an angle to the ray, then the ray will bend, even though the slowness at each point along the path is identical to what it was before when there was no bending. This shows that a ray’s path cannot be explained solely in terms of the value of the slowness on that path. We also must consider the slowness along neighboring paths - hat is, along paths not taken.
The classical wave explanation, proposed by Huygens, resolves this problem by saying that light does not propagate in the form of a single ray. According to the wave interpretation, light propagates as a wavefront possessing transverse width. If a small section of a propagating wavefront encounters the same value of slowness along its width, then the ray will not bend. If a small section of a propagating wavefront encounters different values of slowness along its width, then the ray will bend. The amount of bending depends on the gradient of the slowness. The wavefront propagates more rapidly on the side where the slowness is low (i.e., where the velocity is high) than it does on the side where the slowness is high. As a result, the wavefront naturally turns in the direction of increasing slowness.
In this section, the position vector r always represents a point on a specific raypath and not any arbitrary point in space. As time increases, r traces out the particular raypath in question. The seismic ray at any given point follows the direction of the gradient of the traveltime field Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): t\mathbf(r) . As before, let u be the unit vector along the ray. The ray in general will follow a curved path, and the raypath vector Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): n\mathbf{\ u} will be tangent to this curved raypath.
Let us derive the ray equation. This equation tells us how the raypath vector Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): n\mathbf{u} changes along the curved raypath. The equation for the unit tangent vector is
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \begin{align} \mathbf{u=}&\frac{d\mathbf{r}} {ds}. \end{align} ()
We also know that
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \begin{align} \frac{d\mathbf{r}} {ds}&{\rm =}\left(\frac{dx}{ds},\frac{dy}{ds}\right) . \end{align} ()
Thus, the unit tangent vector u is
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \begin{align} \mathbf{u=}&\left(\frac{dx}{ds},\frac{dy}{ds}\right) . \end{align} ()
The slowness Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): n\left(x,y\right) is a scalar function that depends on coordinates x and y. However, we can hold y constant and consider the curve that gives the variation of slowness with x alone. The slope of this curve Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \partial n/\partial x is called the partial derivative of slowness with respect to x. The partial derivative Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \partial n/\partial y is defined similarly. The gradient of the slowness surface $ n\left(x,y\right) $ is the vector with these partial derivatives as components; that is,
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \begin{align} \mathrm{grad}\ \textit{n}&{\ \rm =}\left(\frac{\partial n}{\partial x}{\rm ,\ }\frac{\partial n}{\partial y}\right). \end{align} ()
We know that the gradient of the traveltime is
$ {\begin{aligned}\mathrm {grad} \ {\textit {t}}&{\ {\rm {=}}}\left({\frac {\partial t}{\partial x}},{\frac {\partial t}{\partial y}}\right).\end{aligned}} $ ()
We recall eikonal equation 14. Now we will do the mathematics. We take the derivative of the eikonal equation with respect to path length s. We obtain
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \begin{align} \frac{d}{ds}\ \left(n\mathbf{u}\right)&{\rm =}\frac{d}{ds}{\rm (}\mathrm{grad}\ \textit{t}{\rm )}. \end{align} ()
On the right side, we interchange the two operations to obtain
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \begin{align} \frac{d}{ds}&{\rm (}\mathrm{grad}\ \textit{t}{\rm )=}\mathrm{grad}\ (\frac{dt}{ds}{\rm )}. \end{align} ()
We recognize Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \partial t{\rm /}\partial s as the slowness n. Thus, the right side of equation 23 is the gradient of the slowness, so we can write
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \begin{align} \frac{d}{ds}&{\rm (}\mathrm{grad}\ \textit{t}{\rm )=}\mathrm{grad}\ (n{\rm )}. \end{align} ()
If we put together the above equations, we obtain the ray equation
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): \begin{align} \frac{d}{ds}\ \left(n\mathbf{u}\right)&{\rm =}\mathrm{grad}\ \textit{n}. \end{align} ()
Ray equation 25 says that the rate of change of the raypath vector Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): n\mathbf{u} is equal to the gradient of the slowness. We know that eikonal equation 14 says that the raypath vector is a flow line on the traveltime surface. The ray equation says that the rate of change of the ray-path vector is a flow line on the slowness surface.
Continue reading
| Previous section | Next section |
|---|---|
| Snell’s law | Ray equation for velocity linear with depth |
| Previous chapter | Next chapter |
| Wave Motion | Visualization |
Also in this chapter
- Reflection seismology
- Digital processing
- Signal enhancement
- Migration
- Interpretation
- Rays
- The unit tangent vector
- Traveltime
- The gradient
- The directional derivative
- The principle of least time
- The eikonal equation
- Snell’s law
- Ray equation for velocity linear with depth
- Raypath for velocity linear with depth
- Traveltime for velocity linear with depth
- Point of maximum depth
- Wavefront for velocity linear with depth
- Two orthogonal sets of circles
- Migration in the case of constant velocity
- Implementation of migration
- Appendix B: Exercises
External links
| find literature about Ray equation |
|
|
|
|
|
|
|