3D Traveltime Computation Using Fast Marching Upwind Bilinear Interpolation Method
Sun Zhang-Qing,Sun Jian-Guo,Yue Yu-Bo,Jiang Zhao-Nan
DOI: https://doi.org/10.6038/cjg20150616
2015-01-01
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.