Abstract:As a very important forward computational tool, 3D traveltime computational scheme has been widely used in many 3D seismic data processing techniques such as traveltime inversion, tomography, migration, etc. In the literatures, finite-difference scheme and conventional bilinear interpolation method are frequently used for calculating the 3D traveltime. The finite-difference scheme has high-efficiency, but is not accurate enough. The conventional bilinear interpolation method has high-accuracy, but is not efficient enough. To build a method for computing the traveltime in 3D space with high-efficiency and high-accuracy, we firstly propose the analytic solutions of the binary minimum transcendental equation by making full use of the results of the plane wave bilinear assumption. Using the above analytic solutions and the deducing process of the conventional bilinear interpolation formulas, we derive some new local computational formulas of 3D traveltime. Secondly, we present a local computational strategy of 3D traveltime by introducing the upwind idea that was proposed in the fast marching method. The above new local computational formulas and strategy make up a new local computational method of 3D traveltime named upwind bilinear interpolation scheme. The new local computational scheme of 3D traveltime is more accuracy and efficient than the conventional bilinear interpolation scheme. We also prove that the new local scheme complies with the propagation law of seismic wave and the eikonal equation. Then, to calculate the traveltime distributed in the whole 3D model, we adapt the above local computational formulas and upwind bilinear interpolation strategy to the general narrow band technique of the fast marching method and build a global algorithm. This global algorithm is named fast marching upwind bilinear interpolation method. Finally, we analyze the accuracy and efficiency of the new method in the 3D homogeneous media and 3D layered media, and verify the stability and effectiveness of the new method with a modified 3D Marmousi model. Based on the above analysis and verification, we can make the following conclusions: (1)A constructive proof indicates that the new method proposed in this paper abides strictly by the Fermat principle and eikonal equation; (2) The new method presented here is accurate and efficient than the 1st-order finite-difference method; (3)For introducing the upwind idea into the conventional bilinear interpolation formulas, the new method proposed in this paper is concise and unconditionally stable; (4)The new method presented here is stable and flexible in 3D complex media.

3D Traveltime Computation Using Fast Marching Upwind Bilinear Interpolation Method

Traveltime computation using the upwind finite difference method with nonuniform grid spacing in a 3D undulating surface condition

Computation of multiple seismic traveltime in mountainous areas with complex 3D conditions using the multistage group marching upwind hybrid method

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

An Adaptive Finite-Difference Method for Seismic Traveltime Modeling Based on 3D Eikonal Equation

A multistage upwind interpolation method for computing traveltimes and raypaths of seismic waves including complex seabed topography

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

3-D Traveltime Computation Using Interpolation of Traveltime Perturbations

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

Traveltimes computation using linear interpolation and narrow band technique under complex topographical conditions

Linear Traveltime Perturbation Interpolation: a Novel Method to Compute 3-D Traveltimes

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

Travel-time computation based on linear interpolation and narrow-band technique

Finite-difference Calculation of Traveltimes Based on Rectangular Grid

Traveitimes computation using linear interpolation and narrow band technique under complex topographical conditions

The Comparison of Three Schemes for Computing Seismic Wave Traveltimes in Complex Topographical Conditions

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

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

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

An Optimized Three-Dimensional Time-Space Domain Staggered-Grid Finite-Difference Method

A fast algorithm for computing complex traveltime based on L-BFGS method