摘要：有限差分方法是波场数值模拟的一个重要的方法,其在地震波正演模拟及偏移成像等方面得到了广泛应用.但常规的有限差分法本身存在着数值频散问题,会降低波场模拟的精度与分辨率,为了克服常规差分算子的数值频散,这里采用了25点有限差分算子,再根据最优化理论求取的优化系数,建立了频率空间域中弹性波波动方程的差分格式;为了消除边界反射,引入了最佳匹配层,构造带有最佳匹配层的差分方程.另外,针对有限差分算法存在计算效率低的问题,本文对有限差分进行了并行算法的分析与研究,给出了相关算法描述.最后由弹性波波动方程和边界条件,通过频率域有限差分法,分别利用不同震源对弹性波在均匀各向同性介质及层状介质中的传播过程进行了数值正演模拟.通过正演模拟,得到了单频波波场、时间切片和共炮点道集,为下一步的研究工作如成像、反演打下了研究基础.
