论文详情
三维地震波动方程数值求解对计算机的内存大小和运算速度都有很高的要求。采用基于消息传递接口(Message Passing Interface,MPI)的并行算法对三维空间的瑞雷面波进行了交错网格有限差分正演模拟。该算法将待模拟区域划分为若干个子区域，各个进程互相协同，并行完成各个子区域的数值模拟过程，从而达到扩大模型规模、加快模拟速度的目的。数值模拟过程中，采用声学-弹性界面法处理自由地表边界。利用均匀各向同性介质模型模拟所得的单道地震记录与解析解的对比结果和波场快照验证了算法的可行性和正确性；通过3层速度递增模型数值模拟所得波场记录的频散曲线与解析解对比，进一步验证了算法的有效性。
Numerical simulation of three-dimensional wave equation requires large computer memory and high calculation speed.Parallel algorithm based on Message Passing Interface (MPI) was adopted to simulate Rayleigh wave by using staggered-grid finite difference method.The model region to be simulated was divided into several sub-regions.Then all processors collaborate each other and the simulation in sub-regions is completed by each processor in parallel.Accordingly,the purposes of accelerate the computing speed and expand the scale of the model are achieved.Acoustic-elastic boundary approach was adopted to implement the traction-free boundary in the simulation.Feasibility and correctness of the simulation method were certificated by comparing the single-channel seismic records from the heterogeneous isotropic medium model obtained by numerical simulation with the results of analytical solution and the snapshots.The method is further verified by the comparison of dispersion curves derived from wave records and analytical solution of a three-layer velocity incremented model.