起伏地表叠前时间偏移的多级并行优化技术.pdfVIP

第 52卷第3期 石 油 物 探 Vo1．52，No．3 2013年 5月 GEOPHYSICALPROSPECTING F0RPETR0LEUM May，2013 文章编号：1000—1441(2013)03—0280一O8 起伏地表叠前时间偏移的多级并行优化技术 马召贵 ，赵改善 ，武港山。，孙成龙 ，亢永敢 ，杨祥森 ，曹永生 (1．中国石油化工股份有限公司石油物探技术研究院，江苏南京211103；2．南京大学软件新技术 国家重点实验室，江苏南京210093) 摘要：Kirchhoff积分法叠前时间偏移并行实现方案除了要考虑数据体规模外 ，还必须考虑大规模节点部署条件 下的并行执行效率。而在 目前的数据处理规模及集群配置条件下，利用多进程进行成像域并行的实现策略存在 内存占用率高、不能充分发挥并行文件系统优势等问题，特别是在多节点大规模计算时，进程间的网络通讯、同 步开销及I／O问题严重影响并行效率。为了提高起伏地表叠前时间偏移算法在大规模部署条件下的偏移效率 ， 分析了该算法的技术特色和计算特点，提出了一种适用于积分法偏移的多级并行计算架构，采用偏移距域的多 进程数据域并行、地震数据 ]／O与偏移计算的异步并行以及单进程内的多线程成像域并行技术对原始算法进行 了架构升级及算法优化。利用34个节点的SMP集群系统和Panasas并行集群存储系统对优化前、后的算法进 行测试，某工区实际叠前地震数据的测试结果表明，该并行架构能够有效降低网络通讯和数据 I／O对计算性能 的影响，使算法的并行效率和大规模集群环境下部署的可扩展能力同时得到提升。 关键词：叠前时间偏移；起伏地表；地震数据处理；高性能计算；并行效率 DOI：10．3969／j．issn．1000—1441．2013．03．009 中图分类号：P631．4 文献标识码 ：A 起伏地表条件下的叠前时间偏移一般先进行 等 ]于2010年提出了一种适用于大规模数据体的 静校正处理，将叠前数据校正到一个固定面上再进 三维Kirchhoff积分法偏移实现方案，该方案将叠 行固定面偏移_1]。然而，现有的静校正方法都是基 前地震数据按共偏移距道集形式组织，逐个偏移距 于某种理论上的假设来进行的，当这种假设得不到 进行偏移，很大程度上降低了偏移计算对内存的依 满足时，应该抛弃传统的静校正处理，直接进行叠 赖，同时，将成像结果按时问维度进行分块处理，并 前偏移处理，并在偏移过程中解决静校正问题_2I3]。 充分利用本地盘存储数据，很好地解决了大工区大 叠前时间偏移在生产中发挥着极其重要的作用，但 规模数据偏移的内存不足问题，具有很好的工业价 随着勘探工区规模越来越大，以及石油勘探开发对 值。2012年，王华忠等Elo]对以上实现方案进行了 于地震勘探精度的要求越来越高，地震勘探数据量 优化，提出了一种适用于Kirchhoff积分法叠前深 与计算量都有了指数级的增长[4]。实践表明，在 度偏移的并行实现方案，该方案利用 OpenMP多 常规集群系统上进行三维海量地震数据的叠前时 线程技术实现了成像域沿时间切片的任务级并行， 间偏移甚费机时，时常使计算成果的产出时间与生 结合多线程共享内存的优势提高了多核处理器的 产进度需求之间产生严重脱节的局面_6]。 使用效率。偏移计算时的数据 I／O会严重影响偏 起伏地表 Kirchhoff叠前时间偏移具有计算 移效率，为了降低 I／O对计算效率的影响，王霖 等_1提出了一种叠前时间偏移并行模式的流水线 密集和I／0密集的特点，并行计算方案需要针对 改进方法，该方法通过建立数据缓