Seismic ray-tracing calculation based on parabolic travel-time interpolation

Abstract

A new seismic ray-tracing method is put forward based on parabolic travel-time interpolation (PTI) method, which is more accurate than the linear travel-time interpolation (LTI) method. Both PTI method and LTI method are used to compute seismic travel-time and ray-path in a 2-D grid cell model. Firstly, some basic concepts are introduced. The calculations of travel-time and ray-path are carried out only at cell boundaries. So, the ray-path is always straight in the same cells with uniform velocity. Two steps are applied in PTI and LTI method, step 1 computes travel-time and step 2 traces ray-path. Then, the derivation of LTI formulas is described. Because of the presence of refraction wave in shot cell, the formula aiming at shot cell is also derived. Finally, PTI method is presented. The calculation of PTI method is more complex than that of LTI method, but the error is limited. The results of numerical model show that PTI method can trace ray-path more accurately and efficiently than LTI method does.

Keywords

ray-tracing / parabolic travel-time interpolation method / linear travel-time interpolation method / travel-time / ray-path / grid cell

References

[1] VidaleJ. Finite-difference calculation of travel-times [J]. Bulletin of the Seismological Society of America, 1988, 78(6): 2 062-2 076

[2] AsawakaE. Seismic ray tracing using linear travel-time interpolation[J]. Geophysical Prospecting, 1993, 41(1): 99-112

[3] VidaleJ. Finite-difference calculation of travel-times in three dimensions [J]. Geophysics, 1990, 55(5): 521-526

[4] MoserT J. Shortest path calculation of seismic rays [J]. Geophysics, 1991, 56(1): 59-67

[5] ZhangSai-minCross Hole Seismic Tomography [D], 2003, Changsha, School of Info-Physics and Geomatics Engineering, Central South University(in Chinese)

[6] QinF, OlsenK. Finite-difference solution of the eikonal equation along expanding wavefronts [J]. Geophysics, 1992, 57(3): 478-487

[7] CaoS, GreenhalghS. Finite-difference solution of the eikonal equation using an efficient, first-arrival wavefront tracking scheme [J]. Geophysics, 1994, 59(4): 632-643

[8] VirieuxJ. Fast and accurate ray tracing by Hamiltonian perturbation [J]. Geophys Resource, 1991, 96(5): 579-594

[9] ChenJing-bo, QinMeng-zhao. Raytracing, symplectic algorithms and wave field simulation [J]. Chinese Journal of Computational Geophysics, 2001, 18(6): 481-486(in Chinese)

[10] LiuSheng-dong, LiCheng-hua, YangWei. Ray tracing method in seismic chromatography CT [J]. Journal of Huai-nan Institute of Technology, 1999, 19(1): 1-4(in Chinese)

[11] ZhangPo, LiuHong, LiYou-ming. The situation and progress of ray tracing method research [J]. Progress of Geophysics, 2000, 15(1): 37-45(in Chinese)

[12] XuKun, WuLu, WangMiao-yue. Improved moser ray-tracing [J]. Progress of Geophysics, 1998, 13(4): 60-66(in Chinese)

[13] LiJia-kang. Gain of near surface velocity in mountain regions [J]. Petroleum Exploration and Development, 2000, 27(2): 60-61(in Chinese)

[14] ZhaoYi, GaoEr-gen, LiuTong-qing, et al.. Total iterative ray-tracing in 3-D layered structure [J]. Computing Techniques for Geophysical and Geochemical Exploration, 2001, 23(2): 105-108(in Chinese)

[15] LiJia-kang, YuQin-fan. Constrained tomographic inversion for near-surface velocity [J]. Oil Geophysics Prospecting, 2001, 36(2): 135-140(in Chinese)

[16] GaoEr-gen, XuGuo-ming. Iterative ray-tracing method segment by segment under 3-D construction [J]. Oil Geophysics Prospecting, 2002, 37(1): 11-16(in Chinese)