摘要

The shortest path raytracing algorithm in networks is widely applied to calculate seismic ray paths and traveltimes, in which the beelines between network nodes are used to approximate the ray paths. When using a few nodes the ray usually follows a zig-zag path that is longer than the real raypath, and the calculated traveltime at a receiver is also larger than the real one. To solve this problem, we develop a shortest path raytracing method with dynamic networks based on the shortest path algorithm. In the process of expanding wavefront away from a source throughout the entire model, the wavefront times are sampled at all nodes on the cell edges. In each cell the traveltime at any point on an edge is expressed with linear interpolation between calculated traveltimes at adjacent nodes on the same cell edge, and then Fermat's Principle is used to determine the first arrival times at the nodes on the other edges of this cell. After the first arrival times at all nodes are calculated, trace raypath from a receiver to a source in the same way as the forward sampling for all pairs of receives and sources is performed. So the drawbacks of conventional shortest path raytracing methods are overcome, and the precision of calculated minimum traveltimes and the correspondent raypaths are improved greatly.