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 . As before, let u be the unit vector along the ray. The ray in general will follow a curved path, and the raypath vector will be tangent to this curved raypath.
Let us derive the ray equation. This equation tells us how the raypath vector changes along the curved raypath. The equation for the unit tangent vector is
We also know that
Thus, the unit tangent vector u is
The slowness 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 is called the partial derivative of slowness with respect to x. The partial derivative is defined similarly. The gradient of the slowness surface is the vector with these partial derivatives as components; that is,
On the right side, we interchange the two operations to obtain
If we put together the above equations, we obtain the ray equation
Ray equation 25 says that the rate of change of the raypath vector 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.
|Ray equation for velocity linear with depth
Also in this chapter
- Reflection seismology
- Digital processing
- Signal enhancement
- The unit tangent vector
- 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