Abstract
The seismic rays and wavefront curvatures are determined by solving a system of nonlinear ordinary differential equations. For media with constant velocity and for media with constant velocity gradient, simplified solutions exist. In a general inhomogeneous medium these equations must be solved by numerical approximations. The integration of the ray‐tracing and wavefront curvature equations is then performed by a modified divided difference form of the Adams PECE (Predict‐Evaluate‐Correct‐Evaluate) formulas and local extrapolation. The interfaces between the layers are represented by bicubic splines. The changes in ray direction and wavefront curvature at the interfaces are computed using standard formulas. For three‐dimensional media, two quadratic traveltime approximations have been proposed. Both are based on a Taylor series expansion with reference to a ray from a reference source point to a reference receiver point. The first approximation corresponds to expanding the square of the traveltime in a Taylor series and taking the square root of the result. The second approximation corresponds to expanding the traveltime in a Taylor series. The two traveltime approximations may be expressed in source‐receiver coordinates or in midpoint‐half‐offset coordinates. Simplified expressions are obtained when the reference source and receiver coincide, giving zero‐offset approximations, for which the reference ray is a normal‐incidence ray. A new method is proposed for computing the second derivatives of the normal‐incidence traveltime with respect to the source‐receiver midpoint coordinates. By considering a beam of normal‐incidence rays it is shown that the second‐derivative matrix may be found by computing the wavefront curvature along a reference normal‐incidence ray starting at the reflection point with the wavefront curvature equal to the curvature of the reflecting interface. From this second‐derivative matrix the normal moveout velocity can be computed for any seismic line through the reference source‐receiver midpoint. It is also shown how a reverse wavefront curvature calculation may be used, in a time‐to‐depth migration scheme, to compute the curvature of the reflecting interface from the estimated second derivatives of the normal‐incidence traveltime. Numerical results for different three‐dimensional models indicate that the first traveltime approximation, based on an expansion of the square of the traveltime, is the most accurate for shallow reflectors and for simple models. For deeper reflectors the two approximations give comparable results, and for models with complicated velocity variations the second approximation may be slightly better than the first one, depending on the particular model chosen. A simplified traveltime approximation may be used in a three‐dimensional seismic velocity analysis. Instead of estimating the stacking velocity one must estimate three elements in a [Formula: see text] symmetric matrix. The accuracy and range of validity of the simplified traveltime approximation are investigated for different three‐dimensional models.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.