Abstract:A fractured area, such as a fault area, usually induces orthorhombic anisotropy. Ignoring orthorhombic anisotropy may degrade the subsurface image by creating a well mistie and blurring the image. Traveltime computation is essential for many processing techniques, such as depth imaging and tomography. Solving the ray-tracing system and eikonal equation are two popular methods for traveltime computation in isotropic media. However, because the ray-tracing system becomes complex and the eikonal equation becomes highly nonlinear, their applications in orthorhombic media become complex. We have developed an alternative 3D traveltime computation method in orthorhombic media based on dynamic programming. To avoid solving the complex ray-tracing system and nonlinear eikonal equation, it adopts an explicitly expressed group velocity from the moveout approximation to describe the propagation of the wavepath and computes the traveltime by Fermat’s principle. Similar to depth extrapolation, it computes the traveltime from one depth to the next depth and does not suffer from a shadow zone. Besides, three strategies of traveltime computation are proposed to deal with different geologic scenarios. Because classic dynamic programming (i.e., the first strategy) computes all possible wavepaths (i.e., 24 wavepaths) across one spatial location, it may be computationally intensive. Based on the idea of wavefield decomposition (e.g., upgoing and downgoing), the second and third strategies simplify the traveltime computation and reduce the computational cost. Numerical examples on the vertical and tilted orthorhombic models indicate that the traveltime contour obtained by our method matches well with the wavefront extrapolated from the wave equation. Our method can be applied in depth imaging and tomography.

High-efficient raytracing method in 1D horizontal layered VTI medium

A Fast and Robust Two-Point Ray Tracing Method in Layered Vertical Transversely Isotropic Media with Strong Anisotropy

Wave Equation Migration for 3-D VSP Data

3D Traveltime Computation in Orthorhombic Media Using Layer-by-layer Wavefront Marching

Fast and Accurate 3-D Ray Tracing Using Bilinear Traveltime Interpolation and the Wave Front Group Marching

3D Traveltime Computation for Quasi-P-wave in Orthorhombic Media Using Dynamic Programming

Finite-difference Calculation of Traveltimes Based on Rectangular Grid

Traveltime calculations for qP, qSV, and qSH waves in 3D tilted transversely isotropic media using fast sweeping method with a Newton method local solver

Seismic complex ray tracing in 2D/3D viscoelastic anisotropic media by a modified shortest-path method

3-D Traveltime Calculation In Arbitrary Velocity Distribution With Dynamic Programming Approach

A Fast Sweeping Method for Calculating Qp-Wave Traveltimes in 3-D Vertical Transversely Isotropic Media Using a Quadratic Equation

A 3D Reflection Ray-Tracing Method Based on Linear Traveltime Perturbation Interpolation

A New Ray Tracing Technique for Crosshole Radar Traveltime Tomography Based on Multistencils Fast Marching Method and the Steepest Descend Method

First-arrival Traveltime Computation for Quasi-P Waves in 2D Transversely Isotropic Media Using Fermat's Principle-Based Fast Marching

A Novel 3d Ray-Tracing Acceleration Technique Based On Kd-Tree Algorithm For Radio Propagation Predictionin Complex Indoor Environment

A search algorithm for first break traveltime in 2D TTI medium

Local angle-domain ray tracing and prestack depth migration in TI medium

An Efficient and Accurate OpenMP-MPI Based 3D Beam Tracing Method for Radio Propagation Simulation in Large Scale Urban Environments

Application of grid successive subdivision algorithm in 3-D seismic ray tracing

Study of Shortest Path Method of Ray Tracing Algorithm

THE APPLICATION OF LINEAR TRAVELTIME INTERPOLATION RAY TRACING TO THE CALCULATION OF FIRST-BREAK TRAVELTIME FOR THE NEAR-SURFACE MODEL