求解3D弹性波方程的并行WNAD方法及其TI介质中的波场模拟

准确模拟TTI介质中弹性波的传播是研究地震各向异性、AVO反演的基础.在二维加权近似解析离散化(WNAD)算法的基础上,本文发展的并行WNAD算法是一种研究三维横向各向同性(TI)介质中弹性波传播的、快速高效的数值模拟方法.我们首先介绍三维WNAD方法的构造过程,然后与经典的差分格式——交错网格(SG)算法进行了比较.理论分析和数值算例表明,WNAD算法比交错网格算法更适合在高性能计算机上进行大规模弹性波场模拟.同时,本文利用并行的WNAD方法研究了弹性波在TTI介质中的传播规律,观测了TI介质中弹性波传播的重要特征:横波分离、体波耦合和速度各向异性等.在TTI介质分界面处,弹性波产生更加复杂的折射、反射和波型转化,使得波场非常复杂,研究和辨别不同类型的波能够加深我们对由裂隙诱导的各向异性介质的认识.