Abstract:
Seismic ray tracing algorithms play an important role in seismology and seismic exploration. The modified shortest path method provides more candidate ray directions by adding secondary nodes to the grid surface, which has the advantages of stability, efficiency, and controllable accuracy. However, the addition of secondary nodes significantly increases memory requirements. In this study, we first used a special node numbering method to dynamically store wavefront grids. Then, in conjunction with the modified shortest path method, we developed the shortest path ray tracing method based on dynamic storage of seismic wavefront. The dynamically stored shortest path method (DSPM) can be implemented by revising the modified shortest path method, and both have consistent computational accuracy. Numerical simulation results show that the DSPM algorithm can trace the ray trajectory and travel time in complex media while reducing memory requirements by approximately one order of magnitude and maintaining computational efficiency to a certain extent. This method is suitable for performing parallel ray tracing on a multicore computer and provides an efficient parallel forward algorithm for large-scale inversion tasks.