曙光4000A上的三维叠前深度偏移并行计算.pdfVIP

曙光4000A上的三维叠前深度偏移并行计算 孔祥宁 方伍宝 赵改善袁联生 中国石化石油勘探开发研究院南京石油物探研究所，南京，210014 摘要波动方程叠前深度偏移在地震勘探成像处理方面起着不可替代的作用，随着高性能大 规模并行计算机技术的发展，波动方程叠前深度偏移计算在地震勘探中的应用有了很大的进步。 波动方程叠前深度偏移处理中庞大的数据规模与海量计算对计算性能提出了很高的要求。曙光 4000A超级计算机系统是我国目前峰值速度最快的商用超级计算机系统，无论是硬件平台建设 还是应用软件的配置方面，都可获得良好的应用性能，基于该系统设计的三维波动方程叠前深 度偏移(炮域)PSDM软件，采用动态负载平衡并行计算模式，具有较高的计算效率、高度的 可扩展性和可靠性。 关键词曙光4000A超级计算机叠前深度偏移并行计算负载平衡 查明地下复杂构造和深层构造已成为目前石油勘探的重点之一，波动方程叠前深度偏移 因其对该类构造具有良好的成像效果而被认为是最为有效的方法之一，在地震勘探成像处理 方面起着不可替代的作用。庞大的计算量与数据量使得叠前深度偏移的理论研究及其应用推 广一直紧随着计算机技术的发展而进步。随着高性能大规模并行计算机和并行处理技术的发 展以及并行计算环境和并行编译器的出现，使得地球物理学家研究叠前深度偏移得到深化和 发展。近几年来，地震叠前深度偏移技术正逐步从克希霍夫算法向波动方程算法发展，从而 可以将叠前深度偏移技术的应用领域从复杂构造成像扩大到复杂构造条件下的岩性研究和油 藏成像。中国石化石油勘探开发研究院南京石油物探研究所和同济大学海洋与地球科学学院 等单位开展了技术合作，成功开发出一套具有自主知识产权的波动方程地震叠前偏移成像软 单炮波动方程叠前深度偏移计算的技术特点，采用一种粗粒度的、动态负载平衡的并行策 略，这样不仅数据流程明晰简洁，而且使叠前深度偏移程序具有很好的系统容错功能。 曙光4000A超级计算机系统是我国目前峰值速度最快的商用超级计算机系统，该系统 KVM维护网络共同组成。全机总共有2128个64位AMD2．4G的CPU，每个结点 Opteron 8．0 Linue 目前安装的操作系统为Turbo 64位服务器版本，向下兼容32位版本的应用软件。 3．0．1基本库。 系统软件包含有gcc3．2．2—5编译器、MPICH1．2．5并行编译环境和fftw 负载共享软件LSF(LoadSharing 软件，在使用范围上，LSF既可用于科学计算，也可用于企业的事务处理；在功能上，除了 一22；一 具有一般的作业管理特性外，在负载平衡、系统容错、检查点操作、进程迁移等方面都有一 定的独特之处和强大功能。网络系统具有高带宽、低延迟、可扩展及易实现等特点。 针对曙光4000A超级计算机系统64位集群计算平台的技术优势和波动方程叠前深度偏 移计算的技术特点，我们对相关程序进行了优化，充分发挥曙光4000A的技术优势和性能 优势，对实际三维地震数据进行计算，取得了较好的成果。通过实际数据的计算结果表明， 并行算法是高效的，同时具有良好的可扩展性，在大规模配置的集群计算机系统下运行，能 保持很高的并行加速比。 1 波动方程叠前深度偏移的并行算法 由于单炮波动方程叠前深度偏移巨大的计算量，因此必须采用并行计算的策略来实现其 实际的应用；而叠前深度偏移问题中所固有的可分解性和线性叠加性质，使我们拥有了天然 的并行能力，为我们实现波动方程叠前深度偏移的并行计算提供了保障。 进行基于单炮波动方程叠前深度偏移成像时，共炮域采用分步傅里叶法(SSF)，傅里 相应的成像点位置进行偏移成像。主要的循环体是一个炮集，在循环体中每一次偏移成像都 是一个相对独立的过程，计算时相互间不需要或很少需要进行数据交换，具有很强的并行 性。因此，按炮集划分进行并行计算是很容易理解的。将一个炮集的深度偏移作为一个并行 粒度，在一个CPU上