Finite difference calculation of traveltime on non-orthogonal grid |
| |
| Authors: | Xiangfang Zeng Sidao Ni |
| |
| Affiliation: | State Key Laboratory of Geodesy and Earth's Dynamics, Institute of Geodesy and Geophysics, Chinese Academy of Sciences, Wuhan 430077, China |
| |
| Abstract: | Finite difference methods have been widely employed in solving the eikonal equation so as to calculate traveltime of seismic phase. Most previous studies used regular orthogonal grid. However, much denser grid is required to sample the interfaces that are undulating in depth direction, such as the Moho and the 660 km discontinuity.Here we propose a new finite difference algorithm to solve the eikonal equation on non-orthogonal grid(irregular grid).To demonstrate its efficiency and accuracy, a test was conducted with a two-layer model. The test result suggests that the similar accuracy of a regular grid with ten times grids could achieve with our new algorithm, but the time cost is only about 0.1 times. A spherical earth model with an undulant660 km discontinuity was constructed to demonstrate the potential application of our new method. In that case, the traveltime curve fluctuation corresponds to topography. Our new algorithm is efficient in solving the first arrival times of waves associated with undulant interfaces. |
| |
| Keywords: | finite difference traveltime non-orthogonal grid topography |
| 本文献已被 CNKI 等数据库收录！ |
| 点击此处可从《地震科学（英文版）》浏览原始摘要信息 |
|
点击此处可从《地震科学（英文版）》下载免费的PDF全文 |
|
