A finite difference scheme for solving the eikonal equation including surface topography
A finite difference scheme for solving the eikonal equation including surface topography
复制标题
DOI：
10.1190/1.3580634
复制
发表时间：
2011-06
期刊：
影响因子：
3.3
通讯作者：
Jianguo Sun;Zhangqing Sun;Fuxing Han
中科院分区：
文献类型：
--
作者：
Jianguo Sun;Zhangqing Sun;Fuxing Han
为了求解曲面和曲面界面附近区域的程函方程，我们分别在曲面和界面附近区域采用非均匀网格间距的二阶迎风差分格式。具体而言，在地球表面和所考虑的接口的直接邻域中，我们用表面（接口）点和所考虑的网格点之间的垂直距离来替换垂直方向上的规则网格间距。然而，对于水平方向，仅使用规则网格点。将传统的迎风差分格式转化为非均匀网格格式。此外，为捕获和传播的局部波前附近的弯曲地球的表面（接口），我们适应快速推进方法，通过引入新的点类型，即表面点，点以上的表面，接口点，和点下的接口。如果我们……
For solving the eikonal equation in the regions near the curved earth’s surface and the curved interface, we find a second order upwind finite difference scheme that uses nonuniform grid spacing in the regions near the earth’s surface and the interface, respectively. Specifically, in the direct neighborhood of the earth’s surface and of the considered interface, we replace the regular grid spacing in the vertical direction by the vertical distance between the surface (interface) point and the grid point under consideration. For the horizontal direction, however, only the regular grid points are used. As a result, the conventional upwind finite difference formulas are changed into the ones with nonuniform grid spacing. Furthermore, for capturing and propagating the local wavefront near the curved earth’s surface (interface), we adapt the fast marching method by introducing new point types, namely the surface point, the point above the surface, the interface point, and the point under the interface. If we u...