基于WNAD方法的非一致网格算法及其弹性波场模拟

详细信息   

• 责任者：宋国杰^c清华大学^d1978^e博士研究生 ; 杨顶辉 ; 陈亚丽^c西南石油大学 ; 马啸
• 刊名：地球物理学报
• 出版年：2010
• 卷期：53^b8
• 关键词：地震勘探20^a有限差分法 ; 算法 ; 数值模拟30^a非一致网格 ; 加权的近似解析离散化方法 ; 数值频散
• 页码：1985～1992
• 索书号：P260.66/141
• 图表：^c表1参15插图
• 分类号：2100
• 文摘期号：2011/03

摘要

加权近似解析离散化(WNAD)方法是近年发展的一种在粗网格步长条件下能有效压制数值频散的数值模拟技术。在地震勘探的实际应用中，不是所有情况都适合使用空间大网格步长。为适应波场模拟的实际需要，该文给出了求解波动方程的非一致网格上的WNAD算法。这种方法在低速区、介质复杂Ⅸ域使用细网格，在其他区域采用粗网格计算。在网格过渡区域，根据近似解析离散化方法的特点，采用了新的插值公式，使用较少的网格点得到较高的插值精度。数值算例表明，非一致网格上的WNAD方法能够有效压制数值频散，显著减少计算内存需求景和计算时间，进一步提高了地震波场的数值模拟效率。