# Ray equation for velocity linear with depth

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 |

Suppose we have a 2D stratified medium with horizontal coordinate *x* and depth coordinate *y*. We work in the first quadrant, so the surface of the earth is the line . The *y* axis in the upward direction indicates depth. We assume that , *a* are two known positive constants. Let the velocity increase linearly with depth *y* but not vary with horizontal coordinate *x*. It follows that slowness decreases with depth. Because

**(**)

it follows that the gradient of the slowness is

**(**)

We see that in this case the gradient of slowness is a vector that always points in the negative direction on the depth axis. In other words, the gradient of slowness points vertically toward the surface of the earth.

In ray theory, it is conventional to let be the angle between the positive *y* direction and the raypath. Thus, the unit tangent vector to the ray is

**(**)

It seems odd to see the sine and cosine in this order, but it is a built-in pitfall of ray theory that one must get used to. The ray equation

**(**)

becomes

**(**)

Because , the first component of the ray equation **30** is

**(**)

which says that is a constant. This constant is designated by *p* and is called the Snell parameter. Thus, the first component 31 of the ray equation gives the so-called Snell equation,

**(**)

The sine refers to the horizontal component. Equation **32** says that the horizontal component of the raypath vector *n***u** is the same at all points along the raypath. At the shot point, the ray makes the angle with the vertical. Choose a specific point (*x*,*y*) on the raypath. Designate the angle at this point by . The Snell equation says that

**(**)

The second component of ray equation **30** is

**(**)

The cosine refers to the vertical component. If we carry out the differentiation on the left side of equation **34**, we obtain

**(**)

which gives

**(**)

or

**(**)

or

**(**)

It follows that the second component 34 of the ray equation gives

**(**)

The first component 31 of the ray equation (Snell’s law) says that *p* is constant. Thus, the second component 34 of the ray equation gives

**(**)

This equation is the differential equation of a circle with radius

**(**)

In conclusion, the solution of the ray equation tells us that the raypath is an arc of a circle of radius

## Continue reading

Previous section | Next section |
---|---|

Ray equation | Raypath 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
- 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