波动方程三维叠前深度偏移并行计算流程探索.pdfVIP

尽量利用Linux系统提供的操作系统命令予以实现。节点之间的通讯主要含数据传递、请求 传递、优先级排队等三种基本操作。通过这三种基本操作的有效组合，则可以完成节点任务 静态、动态分配、节点流程监控等任务。图4是节点A偏移进度落后，从而向其它节点求 助时节点任务动态分配示意图。 3．2．3 恢复点设置三维叠前深度偏移处理中耗费大量机时和计算机资源，计算时，将 偏移过程中单炮偏移结果或单频率偏移结果写人本节点磁盘，作为恢复点。设置恢复点后， 程序计算过程中出现异常中断后，可以随时恢复作业的运行。此外，如果运算过程中某一 节点出现故障，只需对这一节点进行调换或修理：同时，启动节点动态分配进程可在恢复 点后继续完成偏移作业，故障排除后，该节点可视需要随时加入系统工作，这正是集群式 并行机的灵活性所在。 3．3并行性能分析 衡量并行算法性能通常采用以下三个参数pJ： 1) 序中串行部分所耗的时间，P为处理器个数。由此可见，串行部分所占比例随着问题规模 增大而缩小，因而，求解问题的规模随着处理器增多而扩大时，并行度随之增加，仍可获得 理想的加速比。 I p， 2)并行效率Ep=Sp 3) 可扩展性：可扩展性是并行算法在多大程度上能够有效利用节点数增加的能力的一 种度量。随着节点数目的增加，如果并行效率基本保持不变，或者只是略有下降，则称该算 法在所用并行机上可扩展性良好。 针对上面所设计的炮集域并行方案，在集群式并行机Cluster上进行了测试。所用SW地 600 采样，每道采样点2500个．每炮512道，共约7200炮。成像数据体样点数约为1000X ×200。由于波动方程三维叠前深度偏移的串行算法需要大量频繁的I／O操作，同时对计算 机内存要求较高，而炮集域并行方案实质上是一个数据区域分解、各节点独立承担的过程， 当炮集数据分解到多个节点后，同时降低了对磁盘I／O与计算机内存的要求，正因此，炮集 域并行方案在集群式并行机Cluster上获得了线形加速比。从而表明，该方案具有甚优的可 扩展性。图5给出SW地区中一束测线计11．5kin2的三维地震资料计算用时比较。图中左侧 为三维叠前深度偏移单节点与32节点计算用时比较，右侧为用于快速提供精细速度一深度 模型的2．5维叠前深度偏移计算用时比较。可以看出，本文提出的并行算法效率较高。图6 为SW地区Line211经深时转换后的叠前深度偏移剖面。该剖面在信噪比、分辨率，特别是 深层复杂构造成像等诸多方面都获得了令人满意的效果。 4小结 本文针对波动方程三维叠前深度偏移大计算量问题，借助其算法内在并行性，运用集 群式并行机结构的软硬件特点，提出了两种算法处理流程对策，获得以甚低的资金投入， 而成功地实现了工业生产规模的计算。实际资料处理表明，所提出的算法分解流程在集群 式并行机Cluster上获得了线形加速比和甚优的可扩展性。 参考文献 【ll马在田．反射地震掌的现在与未来．篼；同济大学海洋地质与地球物理系主编反射地震学论文集．上海：尉F大学出版社’2I∞0．I．s． theEarth’SInt口ior,NewYork：BlackwdiScientificPublieatioas，1985． [21aa盯b伽峨I．F．血嗡ging 【3】 ：礼农．刘洪．李幼铭等．波动方程辛几何算法三维叠前深度偏移漉程及应用．中国科学院研究生院学报’2∞2．19(2)：134-140． ParallelWorkstationforSoemific ofthc ：4】tSterling,D．Bgc．kcl，D．Savatese．,end．BEOWULF：ACtmltmmtioll,：Proceedings 1995latematioalCoIlfCfgllccOnParallel 1995，Voi．1．1l-14． Pro∞醛ingOCPP)，August 嗍孙家昶，张#渡．迟学斌等．同络并行计算与分布式绾程环境_|匕京：科学出版社．1997． 152 B