一阶速度-应力弹性波方程的有限差分格式.docx

PAGE1

摘要

正演模拟是在利用地球物理方法进行勘探的研究中的一个十分重要的问题，地球物理正演模拟作为地球物理反演的基础，在一定程度上起着非常重要的指导性作用。地震波正演模拟方法可以根据波动方程的种类、模拟技术、网格种类等进行分类。本文选择的一阶速度-应力弹性波方程交错网格有限差分法。在简单均匀介质中对比理论和正演结果验证了方法和自己编写的代码的正确性；然后将精度从时间二阶、空间二阶提高到了时间二阶、空间八阶精度，对比了不同精度条件下的频散结果，发现精度越高，压制频散效果越好；在介质边界的处理上，引入了PML吸收边界，对比了无边界情况下和引入PML边界条件的情况下弹性波在边界上的效