基于MPI的三维瑞雷面波有限差分并行模拟.pdfVIP

第 52卷第 4期 石 油 物 探 Vo1．52，No．4 2013年 7月 GEOPHYSICALPR0SPE rING F0R PETR0LEUM Ju1．，2013 文章编号：1000—1441(2013)04—0354—09 基于MPI的三维瑞雷面波有限差分并行模拟 张明财 ，熊章强 ，张大洲 (1．中国水电工程顾问集团公司西北勘测设计研究院，陕西西安710043；2．中南大学地球科学与 信息物理学院，湖南长沙41O083) 摘要：三维地震波动方程数值求解对计算机的内存大小和运算速度都有很高的要求。采用基于消息传递接VI (MessagePassingInterface，MPI)的并行算法对三维空间的瑞雷面波进行了交错网格有限差分正演模拟。该算 法将待模拟区域划分为若干个子区域，各个进程互相协同，并行完成各个子区域的数值模拟过程，从而达到扩大 模型规模、加快模拟速度的目的。数值模拟过程中，采用声学一弹性界面法处理自由地表边界。利用均匀各向同 性介质模型模拟所得的单道地震记录与解析解的对比结果和波场快照验证了算法的可行性和正确性 ；通过 3层 速度递增模型数值模拟所得波场记录的频散曲线与解析解对比，进一步验证了算法的有效性。 关键词：MPI；瑞雷面波；有限差分 ；并行模拟 DOI：10．3969／j．issn．1000—1441．2013．．04．004 、 中图分类号：P631 文献标识码：A 瑞雷面波凭借其衰减小、信噪比高、抗干扰能 种：以真空形式处理应力的自由边界法[、应力镜 力强、分辨率高以及在层状介质中所具有的频散特 像法[1o]、横向各向同性介质代替自由边界法[11和 性，在地震勘探中有着广泛的应用，其正演模拟也 声学一弹性介质边界代替 自由边界法[12]。由于声 成为国内外学者的研究热点。HestholmE1]通过模 学一弹性介质边界代替 自由边界法符合物理规律， 拟对比二维、三维空间中的瑞雷面波发现，二维介 易于编程实现，应用效果 良好，我们利用该方法处 质难以模拟实际的面波传播特性。董 良国等[2对 理 自由地表 。 一 阶弹性波方程交错网格高阶差分解法进行了研 基于单PC机的串行算法势必会受到当前计 究，给出了差分方程的稳定性条件。Igel等[4]给出 算机数据存储量小、运算速度慢的限制，采用并行 了三维球坐标系下的交错网格有限差分公式，并对 计算技术增加数据存储量、减少正演模拟时间是一 地震波场进行了数值模拟。Peter等[5]利用四阶有 条非常有效的途径。王月英r1。]采用基于 MPI的 限差分交错网格对三维介质中的瑞雷面波进行了 并行算法对三维波动方程进行了有限单元正演模 数值模拟，研究了镜像法处理 自由地表时瑞雷面波 拟，并验证了该算法的可行性和有效性；Minkoff：¨] 的精度与最小波长 内网格点数的关系。熊章强[6] 研究了三维速度一应力有限差分地震波数值模拟方 通过正演模拟，对复杂介质下的瑞雷面波传播特征 法，通过数值模拟提高了运算规模和效率；何兵寿 进行研究，利用高阶交错网格有限差分对垂直界 等[1]研究了利用有限差分法并行求解该方程的基 面、斜界面、垂直低速带等复杂模型进行了正演模 本思路与方法，极大地提高了数值求解弹性波方程 拟计算，并提取相应的频散 曲线进行了分析。周竹 的计算效率。 生等[7]利用交错网格有限差分对弹性介质中的瑞 MPI是一种比较著名的应用于并行环境的消 雷面波进行正演模拟研究，在人工边界上采用了变 息传递标准，其