《基于MPI+CUDA转换波Kirchhoff叠前时间偏移并行算法研究及实现》.pdf

第 卷第 期 石 油 物 探 年 月 文章编号 基于 和 的转换波 叠前时间偏移并行计算 喻 勤 张少华 孔选林 中国石油化工股份有限公司西南油气分公司勘探开发研究院德阳分院四川德阳 中国石油化工股份有限公司多波地震技术重点实验室四川德阳 东方地球物理公司 科技信息处河北涿州 摘要转换波 叠前时间偏移可以实现全空间三维转换波资料的准确成像但转换波叠前偏移数据量 巨大而且需偏移迭代多次来寻找匹配的偏移速度模型导致偏移处理周期长效率低限制了转换波偏移技术在 生产上的应用规模 目前解决上述问题采用的方法主要是应用集群来提高计算效率但集群存在成本高 功耗大占用空间大维护成本高等缺点给出了一种基于 和 的转换波 叠前时间偏移并行算法将细粒度线程级的 并行计算融入粗粒度进程 级 并行编程模型利用实际转换波数据分别在 单核 单卡 和 个节点测试平 台上对算法性能进行了测试结果表明 和 个节点的计算速度是 单核的近 倍可以大 幅度提高转换波 叠前时间偏移的计算效率降低计算成本 关键词转换波地震数据处理 叠前时间偏移集群并行计算 中图分类号 文献标识码 转换波 叠前时间偏移是转换波地 等但因为 面元化不易确定转换点真实位 震数据处理中的一项关键技术可实现全空间三维 置 不能适应层间速度剧烈变化和大陡倾角 情况所以需要进行转换波叠前时间偏移该技术 转换波资料的准确成像 但转换波勘探数据量 不需要抽取 道集和 处理就能实现全空 巨大进行叠前偏移处理异常耗时为了提高生产 效率目前大多数偏移算法采用大规模 集群 间三维转换波资料的准确成像 进行计算但因集群设备存在着成本高功耗大占 根据各向异性双平方根方程叠前时间偏移方 用空间大以及维护成本高等缺点其发展空间受到 程 可以写成 限制近几年应用 高性能计算进行地震偏 移已成为一个新兴的研究热点国内外很多地球物 槡 理公司都开展了基于 的偏移算法研究 使叠前时间偏移叠前深度偏移单程波波动方程 槡 偏移逆时偏移等性能得到数十倍的提升我们在 环境下利用 式中 为散射波旅行时 是半源检偏移距 为 等工具实现了基于 和 炮检距 为各向异性参数 是炮点到 成像点的下行 波走时 是成像点到检波器的 的转换波 叠前时间偏移并行算法 上行 波走时 是下行 波速度 是下行 波 速度由于从实际转换波数据中很难获得 在 转换波 叠前时间偏移原 各向异性速度分析时也得不到 因此这两个参 理 量利用等效的转换波速度 转换得到 在转换波地震数据处理中由于下行 波和 上行 波的路径不对称需要做共转换点 收稿日期 改回日期 选排因此转换波资料处理流程一般为抽 作者简介喻勤 男硕士工程师主要从事转换波处理和 道集速度分析 叠加以及叠后偏移 并行计算方面的研究工作 石 油 物 探 第 卷 等多个版 本最常用和稳定的是 和 是目前超大规模并行计算应用最广效率最 图 为转换波成像示意图转换波叠前时间偏 高的并行编程模型几乎被所有并行环境共享和 移算法是计算出炮点到成像点的下行 波走时 分布式存储并行机集群系统等和多进程操作系 成像点到检波器的上行 波走时 依据散射 统 所支持基于 波旅行时公式将上下行波走时之和 上的 开发的应用程序具有可移植性适用性易用性可 能量累加到成像点处转换波 叠前时 扩展性强性能高等优点通过 可以方便地 间偏移能将任何道集形式的输入转换成共成像点 连接现有计算机进行高性能集群计算 道集任意一个检波点上接收的能量都可 与 的最大优势均在于并行计算 以依据炮检点间的空间关系和各向异性旅行时方 但是 面向的是主处理器 进程间的并行 程分配到空间所有可能的成像位置而所有炮检对 属于粗粒度级的并行而 的并行是协处理 的能量都可以依据方程给出的射线路径累加到图 器线程间的并行属于细粒度级的并行要 中散射点成像点的成像位置从而使转换波在 实现基于 和 的并行算法需要对计算 三维空间任何位置准确成像 任务进行划分偏向管理控制的部分采用进程间的 并行计算计算密度大的部分采用线程间的并行计 %% xP xS % 算通过主处理器 和协处理器 结合 可以最大限度地发挥计算机的处理能力 S P 基于 和 的 叠前时间偏 移并行算法设计主要包括以下几个方面 t 从硬件系统上考虑采取一个 绑定一 个 的方式这样可以形成最小计算粒度的节 点比较灵活地进行数值计算保证其单独享有 和 缓存避免多核 产生的数据竞争问题 % 减少资源需求从而降低内存负载保持 到 h h 显存 信道 到主存以及主存到硬盘总 图 转换波成像原理 线的高带宽 转换波叠前时间偏移从理论上规避了输入数 进程间的通信通过 来完成每个 据零炮检距的假设从而避免了校正叠加 进程绑定一个 主要用于完成均