一种基于共享存储的叠前深度偏移并行算法A-软件学报.PDFVIP

1000-9825/2002/13(12)2231-07 ©2002 Journal of Software 软 件 学 报 Vol.13, No.12 一种基于共享存储的叠前深度偏移并行算法 1 1 2 1 1 李建江 , 舒继武 , 王有新 , 王鼎兴 , 郑纬民 1(清华大学 计算机科学与技术系,北京 100084); 2( 中国石油天然气集团公司 石油地球物理勘探局,河北 涿州 072751) E-mail: jj-li00@ 摘要: 为了解决叠前深度偏移计算量巨大这个问题,人们投入了很大精力来开发高效的并行算法.就此,在对一些相 关并行算法进行分析之后,根据三维叠前 Kirchhoff 深度偏移的特点,提出了一种基于共享存储的简化并行算法.其 中 Slave 进程直接存储和读取射线旅行时,消息传递总量显著减少; 同时合理地结合“任务池”技术与粗粒度级并行, 前者确保动态负载平衡,后者使得Slave 与Master 之间的通信开销尽可能小.实际数据的测试结果表明,此并行算法 高效且可扩展性较好. 关 键 词: 三维叠前Kirchhoff 深度偏移;共享存储;并行算法;消息传递;动态负载平衡;效率;可扩展性 中图法分类号: TP391 文献标识码: A 当前,查明地下复杂构造和深层构造已成为石油勘探的重点之一.深度偏移因其对该类构造具有良好的成 像效果[1~6] 而被公认为目前最为有效的方法之一,但是固有的巨大计算量阻碍了其实用化进程.随着并行计算 机和并行处理技术的发展以及并行计算环境和并行编译器的出现,其实用化成为可能. 目前已有研究针对其并 行实现[1~6],其中文献[1]提出了广义相位屏法单炮叠前深度偏移及其并行实现,文献[2]提出了三维叠后有限差 分深度偏移的并行算法,文献[3]提出了二维叠前有限差分法、裂步傅立叶方法和傅立叶有限差分法深度偏移的 并行实现.然而, 由于目前在石油勘探界使用最为普遍的叠前深度偏移技术是克希霍夫(Kirchhoff)积分法,所以 有必要对其并行算法进行研究.文献[4]提出了一种二维叠前 Kirchhoff 深度偏移的并行算法,在 SGI 和 IBM SP2 上获得了较高的并行效率.文献[5]提出了一种三维叠前Kirchhoff 深度偏移的并行实现方法,在IBM SP2 上 也获得了不错的并行效率.但其并行方案较为复杂、控制较多, 引入了一定的额外开销. 由于该方法在对地震数 据进行成像的过程中,是根据实际情况采用读取本地磁盘、通过消息传递、实时求取这3 种方式之一来获得所 需的射线旅行时, 因此不同 Slave 之间会存在大量的消息传递.文献[6]提出通过动态计算射线旅行时来实现三 维叠前Kirchhoff 深度偏移的并行化,但其对原始地震数据进行了静态分割,负载平衡不是很好, 同时通过设置同 步点来实现不同进程间的同步,这些都影响到算法的并行效率.本文针对三维叠前Kirchhoff 积分法深度偏移的 具体特点,在共享存储的情况下简化其并行实现方案.简单地说,采用 Slave 进程直接存储/读取射线旅行时的策 略,使得消息传递总量显著减少. 同时,合理结合“任务池”技术与粗粒度级并行,既保证了动态负载平衡,又使得 Slave 与 Master 之间的通信开销尽可能地小.实际数据的测试结果表明,本文提出的并行算法是高效的, 同时具 有较好的可扩展性. 收稿日期: 2001-04-20; 修改日期: 2001-06-26 基金项目: 国家自然科学基金资助项目 作者简介: 李建江(1971 －), 男, 四川广安人,博士生,主要研究领域为并行处理技术,并行编译;舒继武(1969 －), 男,湖北孝感人, 博士,副教授,主要研究领域为大规模科学与工程计算中并行算法,并行处理技术,并行应用软件;王有新( 1958－),男, 山东烟台人,博士, 教授级高级工程师,主要研究领域为地震资料数值处理方法研究及其软件开发;王鼎兴(1937 －), 男,江苏吴江人,教授,博士生导师,主 要研究领域为并行计算机体系结构,并行计算;郑纬民(1946 －), 男,浙江鄞县人,教授,博士生导师,主要研究领域为并行计算机体