论文部分内容阅读
地震波场数值模拟一般基于水平地表条件,而当前近地表地震方法研究中急需起伏地表条件下的正演模拟数据.起伏地表条件下的二维地震波场波动方程高阶有限差分数值模拟方法研究中面临的主要问题是边界条件的处理.从以速度-应力表示的一阶波动方程出发,导出了起伏边界情况下的边界条件.该条件与波动方程具有相同的形式,可以用同样的方法来处理正演过程中的网格内部点和边界点.用交锗网格的高阶有限差分方法解波动方程,在满足稳定性要求时,可获得时间和空间都是高阶精度的结果.本算法在高性能计算机集群MPI环境中并行实现时,采用主从模式设