Register      Login
Exploration Geophysics Exploration Geophysics Society
Journal of the Australian Society of Exploration Geophysicists
RESEARCH ARTICLE

2D efficient ray tracing with a modified shortest path method

Meigen Zhang 1 2 4 Liyun Fu 2 Xinfu Li 3 Xiaofan Li 2
+ Author Affiliations
- Author Affiliations

1 Key Laboratory of Engineering Geomechanics, Chinese Academy of Sciences, 19 Beitucheng Xilu, Beijing 100029, China.

2 Key Laboratory of the Earth’s Deep Interior, Institute of Geology and Geophysics, Chinese Academy of Sciences, 19 Beitucheng Xilu, Beijing 100029, China.

3 China University of Geosciences, Beijing 100081, China.

4 Corresponding author. Email: mgzhang@mail.iggcas.ac.cn

Exploration Geophysics 40(4) 301-307 https://doi.org/10.1071/EG09016
Submitted: 6 March 2009  Accepted: 13 November 2009   Published: 7 December 2009

Abstract

The computation effort of ray tracing with the shortest path method (SPM) is strongly dependent on the number of the discretized nodes in a model and the number of ray directions emanating from a secondary source node. In the traditional SPM, a secondary source emanates rays to all the surrounding nodes. Obviously, most of them are not minimal traveltime raypaths. As a result, the efficiency of SPM can be greatly improved if some measures are taken to avoid those unnecessary computations. In the current study, we apply the traveltime information of neighbouring nodes and the incident rays to determine the effective target propagation directions of secondary source nodes in 2D case. Generally, the effective propagation directions are narrow bands with few surrounding nodes. Thus, most unnecessary ray directions of secondary source nodes are avoided. 2D model tests show that the computational speed of the improved method is about several to tens of times of that of the traditional SPM with the increase of network nodes and cells.

Key words: shortest path method, efficient ray tracing, effective propagation direction.


Acknowledgements

This work is sponsored by Key Laboratory of Engineering Geomechanics of Chinese Academy of Sciences, CNPC Innovation Fund and National Natural Science Foundation of China (grants 40437018 and 40874024). We thank the editors and two referees for their useful comments, which were of great help in improving the quality of this manuscript.


References

Asakawa, E., and Kawanaka, T., 1993, Seismic ray tracing using linear traveltime interpolation: Geophysical Prospecting 41, 99–111.
Crossref | GoogleScholarGoogle Scholar | Saito H. , 1989, Traveltimes and raypaths of first arrival seismic waves: computation method based on Huygens’ principle: 59th Annual International Meeting, SEG Expanded Abstracts, 244–247.

Schneider, W. A., Fanzingr, K. A., Balch, A. H., and Kruse, C., 1992, A dynamic programming approach to first arrival traveltime computation in media with arbitrarily distributed velocities: Geophysics 57, 39–50.
Crossref | GoogleScholarGoogle Scholar |

Thurber, C. H., and Ellsworth, W. L., 1980, Rapid solution of ray tracing problems in heterogeneous media: Bulletin of the Seismological Society of America 70, 1137–1148.


Van Avendonk, H. J. A., Harding, A. J., Orcutt, J. A., and Holbrook, W. S., 2001, Hybrid shortest path and ray bending method for traveltime and raypath calculations: Geophysics 66, 648–653.
Crossref | GoogleScholarGoogle Scholar |

van Trier, J., and Symes, W. W., 1991, Upwind finite-difference calculation of traveltime: Geophysics 56, 812–821.
Crossref | GoogleScholarGoogle Scholar |

Vesnaver, A. L., 1996, Ray tracing based on Fermat’s principle in irregular grids: Geophysical Prospecting 44, 741–760.
Crossref | GoogleScholarGoogle Scholar |

Vidale, J., 1988, Finite-difference calculation of traveltimes: Bulletin of the Seismological Society of America 78, 2062–2076.


Wang, Y., 1999, Simultaneous inversion for model geometry and elastic parameters: Geophysics 64, 182–190.
Crossref | GoogleScholarGoogle Scholar |

Wapenaar, C. P. A., Kinneging, N. A., and Berkhout, A. J., 1987, Principle of prestack migration based on the full elastic two-way wave equation: Geophysics 52, 151–173.
Crossref | GoogleScholarGoogle Scholar |

Zhang, J., and Toksoz, M. N., 1998, Nonlinear refraction traveltime tomography: Geophysics 63, 1726–1737.
Crossref | GoogleScholarGoogle Scholar |

Zhang, M., Chen, B., Li, X., and Wang, M., 2006, A fast algorithm of shortest path ray tracing: Chinese Journal of Geophysics 49, 1467–1474.


Zhao, A., Zhang, Z., and Teng, J., 2004, Minimum travel time tree algorithm for seismic ray tracing: improvement in efficiency: Journal of Geophysics and Engineering 1, 245–251.
Crossref | GoogleScholarGoogle Scholar |

Zhou, B., and Greenhalgh, S. A., 2005, ‘Shortest path’ ray tracing for most general 2D/3D anisotropic media: Journal of Geophysics and Engineering 2, 54–63.
Crossref | GoogleScholarGoogle Scholar | CAS |