期刊
GEOPHYSICAL JOURNAL INTERNATIONAL
卷 229, 期 3, 页码 1942-1963出版社
OXFORD UNIV PRESS
DOI: 10.1093/gji/ggac037
关键词
Numerical modelling; Numerical solutions
资金
- MURI Center for Material Failure Prediction through Peridynamics at the University of Arizona (AFOSR) [FA9550-14-1-0073]
The computation of travel times for P and S waves is crucial for earthquake and exploration seismology applications. However, solving the eikonal equation in anisotropic media with a complex phase velocity field is challenging. This study presents a new method using the peridynamic differential operator (PDDO) to solve the eikonal equation, which provides a non-local form and is immune to discontinuities.
Traveltimes of compressional (P) and shear (S) waves have proven essential in many earthquake and exploration seismology applications. An accurate and efficient traveltime computation for P and S waves is crucial for the success of these applications. However, solving the eikonal equation with a complex phase velocity field in anisotropic media is challenging. The eikonal equation is a first-order nonlinear hyperbolic partial differential equation. It represents the high-frequency asymptotic approximation of the wave equation. The fast marching and sweeping methods are commonly used due to their efficiency in numerically solving the eikonal equation. However, these methods suffer from numerical inaccuracy in anisotropic media with sharp heterogeneity, irregular surface topography and complex phase velocity fields. This study presents a new method for the solution of the eikonal equation by employing the peridynamic differential operator (PDDO). The PDDO provides the non-local form of the eikonal equation by introducing an internal length parameter (horizon) and a weight function with directional non-locality. The operator is immune to discontinuities in the form of sharp changes in field or model variables and invokes the direction of traveltime in a consistent manner. The weight function controls the strength of association among points within the horizon. Solutions are constructed in a consistent manner without upwind assumptions through simple discretization. The robustness of this approach is established by considering different types of eikonal equations on complex velocity models in anisotropic media. The examples demonstrate its unconditional numerical stability and results compare well with the reference solutions.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据