A multistage upwind interpolation method for computing traveltimes and raypaths of seismic waves including complex seabed topography

摘要

To compute traveltimes and raypaths of seismic waves of different types in ocean regions with complex seabed topography, a new method, called multistage upwind interpolation method, is presented. The key point of the method is to construct the linear interpolation formulas for the traveltimes and raypaths by using the basic ideas of the upwind finite difference method, the multistage fast marching method, and of the multistage backward ray computation scheme. A numerical evaluation and some numerical examples in some typical models with complex seabed topography show that the multistage upwind interpolation method is accurate, unconditionally stable, and flexible enough when it is used in ocean regions with the complex seabed topography.