科研之友
微信
新浪微博
Facebook
分享链接
摘要
基于弹性波传播方程,发展了一种高精度低数值频散的八阶ONAD(optimal nearlyanalytic discrete)方法,该方法利用八阶精度的近似解析离散算子对空间高阶偏导数进行离散,采用四阶精度的截断泰勒展开式离散时间高阶导数。八阶ONAD方法被用于模拟地震波在VTI介质模型和2个复杂层状介质模型中的传播。计算效率结果表明,该方法在运算速度和存储量上明显优越于八阶LWC方法。波场模拟结果显示,八阶ONAD方法在粗网格条件下可有效消除由速度强间断所造成的数值频散,有利于在强间断介质中使用粗网格进行波场模拟,是一种在地震勘探领域有着巨大应用潜力的数值方法。
- 单位
