First break traveltime tomographic imaging based on the eikonal equation in VTI medium

ZHANG Kai1,2 CHEN Junyi1,2 ZHU Baihang1 LI Zhenchun1,2 DING Pengcheng1

(1.School of Geosciences, China University of Petroleum (East China) , Qingdao, Shandong, China 266580)
(2.Laboratory of Marine Mineral Resources Evaluation and Detection Techniques, National Laboratory for Marine Science and Technology, Qingdao, Shandong, China 266061)
【Knowledge Link】Fast marching method

【Abstract】The first break traveltime tomography is an effective method for the inversion of anisotropic parameter. Its difficulties are the calculation of first break traveltime and joint inversion of anisotropic parameters in VTI medium. In this paper, we use, NMO velocity, and non-elliptic coefficient to characterize the anisotropy. The ray tracing based on Runge-Kutta algorithm is used to calculate Fréchet kernel. According to the group slowness approximate expressions, the Fréchet kernels of these three parameters are derived. In order to overcome the magnitude between the three parameters, a double-round strategy is employed. The least squares orthogonal decomposition (LSQR) algorithm is used to carry out the joint inversion of the two methods. The numerical simulations show that the proposed method has better inversion accuracy than the Thomsen Kutta algorithm is not only applicable to isotropic medium but also to VTI medium.

【Keywords】 VTI medium; eikonal equation; first break traveltime tomography; Fréchet kernel; Thomsen parameter;


【Funds】 National Science and Technology Major Project on Oil and Gas Development (2017ZX05032-003-002) Natural Science Foundation of Shandong Province (ZR2017MD014) Key Research and Development Program of Shandong Province (2018GHY115016) Innovative Scientific Research Program of China University of Petroleum (East China) (17CX02052)

Download this article


    [1] Qin Ning, Li Zhenchun, Zhang Kai et al. Traveltime velocity tomography based on ADCIGs of Gaussian beam. OGP, 2013, 48 (3): 366–371 (in Chinese).

    [2] Thomsen L. Weak elastic anisotropy. Geophysics, 1986, 51 (10): 1954–1966.

    [3] Tsvankin I. Inversion of reflection traveltime for transverse isotropic media. Geophysics, 1995, 60 (1): 268–284.

    [4] Alkhalifah T. Traveltime computation with the linearized eikonal equation for anisotropic Media. Geophysical Prospecting, 2002, 50 (4): 373–382.

    [5] Alkhalifah T. Acoustic approximations for processing in transversely isotropic media. Geophysics, 1998, 63 (2): 623–631.

    [6] Alkhalifah T. The offset-midpoint traveltime pyramid for transversely isotropic media. Geophysics, 2000, 65 (4): 1316–1325.

    [7] Alkhalifah T. An acoustic wave equation for anisotropic media. Geophysics, 2000, 65 (4): 1239–1250.

    [8] Alkhalifah T. Fomel S. Implementing the fast marching eikonal solver: spherical versus Cartesian coordinates. Geophysical Prospecting, 2001, 49 (2): 165–178.

    [9] Sethian J A. Fast marching methods. Siam Review, 1999, 41 (2): 199–235.

    [10] Sethian J A. Level Set Methods and Fast Marching Methods. Cambridge University Press, Cambridge, 1999.

    [11] Fomel S. On an elliptic approximations for qP velocities in transversally isotropic media. Geophysical Prospecting, 2004, 52 (1): 247–259.

    [12] Chapman C H, Pratt R G. Traveltime tomography in anisotropic media—I: Theory. Geophysical Journal International, 1992, 109 (1): 1–19.

    [13] Zhang Wensheng, He Qiaodeng. Statistical tomography in anisotropic media by using crosshole traveltimes of qP wave and qSV wave. Chinese Journal of Geophysics, 1999, 42 (S1): 195–203 (in Chinese).

    [14] Lu Minghui, Tang Jianhou, Yang Huizhu et al. P-wave traveltime analysis and Thomsen parameters inversion in orthorhombic media. Chinese Journal of Geophysics, 2005, 48 (5): 1167–1171 (in Chinese).

    [15] Zhou B, Greenhalgh S A, Green A. Nonlinear traveltime inversion scheme for crosshole seismic tomography in tilted transversely isotropic media. Geophysics, 2008, 73 (4): 17–33.

    [16] Huang Guangnan, Deng Juzhi, Li Hongxing et al. Nonlinear traveltime tomography of P-SV and P-SH waves. OGP, 2015, 50 (6): 1127–1133 (in Chinese).

    [17] Qin Ning, Li Zhenchun, Yang Xiaodong et al. Image domain travel-time tomography velocity inversion based on automatic picking. OGP, 2012, 47 (3): 392–398 (in Chinese).

    [18] Liu Yuzhu, Dong Liangguo. Analysis of influence factor of first-breaks tomography. OGP, 2007, 42 (5): 544–553 (in Chinese).

    [19] Liu Yuzhu, Ding Kongyun and Dong Liangguo. The dependence of first arrival travel time tomography on initial model. OGP, 2010, 45 (4): 502–511 (in Chinese).

    [20] Liu Yuzhu, Wang Guangyin, Dong Liangguo et al. Joint inversion of VTI parameters using nonlinear traveltime tomography. Chinese Journal of Geophysics, 2014, 57 (10): 3402–3410 (in Chinese).

    [21] Jiang F, Zhou H W. Traveltime inversion and error analysis for layered anisotropy. Journal of Applied Geophysics, 2011, 73 (2): 101–110.

    [22] Jiang F, Zhou H W. A strategy to estimate anisotropic parameters by error analysis of traveltime inversion. SEG Technical Program Expanded Abstracts, 2010, 29, 323–327.

    [23] Tsvankin I. P-wave signatures and notation for transversely isotropic media. Geophysics, 1996, 61 (2): 467–483.

    [24] Ding Pengcheng, Yang Guoquan, Li Zhenchun et al. Ray tracing based on 3D multi-stencils fast marching algorithm for complex near-surface model. GPP, 2016, 55 (4): 483–492 (in Chinese).

    [25] Yang Guoquan, Ding Pengcheng, Li Zhenchun et al. 3D first arrival traveltime tomography with linearized eikonal equation and shaping regularization. OGP, 2017, 52 (2): 264–272 (in Chinese).

    [26] Wang Qi, Zhu Pan, Ye Pei et al. Seismic wave travel time mixed grid linear interpolation ray tracing algorithm for relief surface. OGP, 2018, 53 (1): 35–46 (in Chinese).

This Article


CN: 13-1095/TE

Vol 53, No. 06, Pages 1218-1226+1112-1113

December 2018


Article Outline



  • 1 Introduction
  • 2 Principle
  • 3 Numerical simulation
  • 4 Conclusions
  • References