瑞雷面波凭借其衰减小、信噪比高、抗干扰能力强、分辨率高及在层状介质中具有频散特性,在地震勘察和工程检测领域有着广泛的应用,其数值模拟也日益受到学界重视。为了更加接近真实的地质模型,为实际的瑞雷面波勘探提供充足的理论支撑,通过交错网格有限差分数值模拟,对三维复杂介质中的瑞雷面波传播特性进行了研究。三维波动方程数值模拟对计算机内存和计算速度有很高的要求,为了克服单PC串行算法在三维数值计算时数据存储量小和计算速度慢的局限,采用基于消息传递接口(Message Passing Interface, MPI)的并行算法对瑞雷面波进行了交错网格有限差分正演模拟。将模拟区域划分为若干个子区域,各个进程互相协同,并行完成各个子区域的数值模拟,从而达到加速模拟速度,扩大模型规模的目的。数值模拟过程中,采用声学-弹性界面近似法(Acoustic-Elastic boundary Approach,AEA)处理水平自由地表,吸收边界通过指数衰减的方法处理。模拟所得的单道地震记录与解析解的对比结果和波场快照验证了模拟的可行性和正确性。本文对水平层状介质(三层速度递增、含软夹层)、地表裂缝、地下空洞、地表水平倾斜、地表隆起、地表凹陷等模型下的瑞雷面波进行了数值模拟。在模拟起伏地表模型中的瑞雷面波时,采用坐标映射的方法,先将地形函数及其导数引入波动方程后再通过交错网格有限差分求解。通过一系列数值模拟算例,分析了瑞雷面波高阶模式形成机理,直观地揭示了瑞雷面波在复杂介质中的传播过程,并研究了其传播特性。