Abstract：
3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method (FMM) is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method (WFC), which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.
. Joint 3D traveltime calculation based on fast marching method and wavefront construction[J]. APPLIED GEOPHYSICS, 2017, 14(1): 56-63.
[1]
Alkhalifah, T., and Fomel, S., 2001, Implementing the fast marching eikonal solver: spherical versus Cartesian coordinates: Geophysical Prospecting, 49(2), 165-178.
[2]
Benaïchouche, A., Noble, M., and Gesret, A., 2015, First arrival traveltime tomography using the fast marching method and the adjoint state technique: 77th EAGE Conference and Exhibition 2015.
[3]
Coman, R., and Gajewski, D., 2005, Traveltime computation by wavefront-orientated ray tracing: Geophysical Prospecting, 53(1), 23-36.
[4]
Eskandari, M., and Safian, R., 2010, Inverse scattering method based on contour deformations using a fast marching method: Inverse Problems, 26(9), 095002.
[5]
Han, F., Sun, J., Sun, Z., et al., 2009, Positioning of grid points in wavefront construction: Applied Geophysics, 6(3), 248-258.
[6]
Huang, X., Sun, J., and Sun, Z., 2016, Local algorithm for computing complex travel time based on the complex eikonal equation: Physical Review E, 93(4), 043307.
[7]
Iversen, E., and Kaschwich, T., 2016, Twin wavefront construction approach to simulation of shear waves in transversely isotropic layered media: 78th EAGE Conference and Exhibition 2016.
[8]
Julian, B. R., and Gubbins, D., 1977, Three-dimensional seismic ray tracing: J. Geophys, 43(1), 95-114.
[9]
Lelièvre, P. G., Farquharson, C. G., and Hurich, C. A., 2011, Computing first-arrival seismic traveltimes on unstructured 3-D tetrahedral grids using the fast marching method: Geophysical Journal International, 184(2), 885-896.
[10]
Liu, S., Wang, H., Liu, T., et al., 2015, Kirchhoff PSDM angle-gather generation based on the traveltime gradient: Applied Geophysics, 12(1), 64-72.
[11]
Liu, S., Wang, H., Yang, Q., et al., 2014, Traveltime computation and imaging from rugged topography in 3D TTI media: Journal of Geophysics and Engineering, 11(1), 015003.
[12]
Motamed, M., and Runborg, O., 2015, A wavefront-based Gaussian beam method for computing high frequency wave propagation problems: Computers & Mathematics with Applications, 69(9), 949-963.
[13]
Rawlinson, N., and Sambridge, M., 2004, Wave front evolution in strongly heterogeneous layered media using the fast marching method: Geophysical Journal International, 156(3), 631-647.
[14]
Sethian, J. A., 1996, A fast marching level set method for monotonically advancing fronts: Proceedings of the National Academy of Sciences, 93(4), 1591-1595.
[15]
Sethian, J. A., and Popovici, A. M., 1999, 3-D traveltime computation using the fast marching method: Geophysics, 64(2), 516-523.
[16]
Sun, J., Sun, Z., and Han, F., 2011, A finite difference scheme for solving the eikonal equation including surface topography: Geophysics, 76(4), 53-63.
[17]
Sun, Y., and Fomel, S., 1998, Fast-marching eikonal solver in the tetragonal coordinates: 68th Annual International Meeting, SEG, Expanded Abstracts, 1949-1952.
[18]
Sun, Z., Sun, J., and Han, F., 2012, Traveltime computation using the upwind finite difference method with nonuniform grid spacing in a 3D undulating surface condition: Chinese J. Geophys. (in Chinese), 55(7), 2441-2449.
[19]
Szostek, K., and Piórkowski, A., 2016, Real-time simulation of ultrasound refraction phenomena using ray-trace based wavefront construction method: Computer Methods and Programs in Biomedicine, 135, 187-197.
[20]
Vidale, J. E., 1988, Finite-difference calculation of travel times: Bulletin of the Seismological Society of America, 78(6), 2062-2076.
[21]
Vidale, J. E., 1990, Finite-difference calculation of traveltimes in three dimensions: Geophysics, 55(5), 521-526.
[22]
Vinje, V., Iversen, E., Åstebøl, K., et al., 1996, Estimation of multivalued arrivals in 3D models using wavefront construction—Part I: Geophysical Prospecting, 44(5), 819-842.
[23]
Vinje, V., Iversen, E., and Gjøystdal, H., 1993, Traveltime and amplitude estimation using wavefront construction: Geophysics, 58(8), 1157-1166.
2017, Vol. 14 
