适于大规模数据的逆时偏移实现方案-物探与化探.pdfVIP

　 第３９卷第５期 物　 探　 与　 化　 探 Ｖｏｌ．３９，Ｎｏ．５　 　 ２０１５年１０月 ＧＥＯＰＨＹＳＩＣＡＬ ＆ ＧＥＯＣＨＥＭＩＣＡＬ ＥＸＰＬＯＲＡＴＩＯＮ Ｏｃｔ．，２０１５　 ｄｏｉ：１０．１１７２０／ ｗｔｙｈｔ．２０１５．５．１７ － 张慧宇，王立歆，孔祥宁，等．适于大规模数据的逆时偏移实现方案［Ｊ］．物探与化探，２０１５，３９（５）：９７８ ９８４．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０．１１７２０／ ｗｔｙｈｔ．２０１５．５． １７ Ｚｈａｎｇ ＨＹ，Ｗａｎｇ ＬＸ，ＫｏｎｇＸ Ｎ，ｅｔ ａｌ．Ｔｈｅｉｍｐｌｅｍｅｎｔａｔｉｏｎ ｓｃｈｅｍｅｏｆ ｒｅｖｅｒｓｅ⁃ｔｉｍｅｍｉｇｒａｔｉｏｎｆｏｒｍａｓｓｓｅｉｓｍｉｃｄａｔａ［Ｊ］．Ｇｅｏｐｈｙｓｉｃａｌ ａｎｄ ＧｅｏｃｈｅｍｉｃａｌＥｘ⁃ － ｐｌｏｒａｔｉｏｎ，２０１５，３９（５）：９７８ ９８４．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０．１１７２０／ ｗｔｙｈｔ．２０１５．５．１７ 适于大规模数据的逆时偏移实现方案 张慧宇，王立歆，孔祥宁，李博 （中石化石油物探技术研究院，江苏 南京　 ２１１１０３） 摘 要：针对海量地震数据逆时偏移所面临的问题和挑战，解决逆时偏移海量计算和海量存储需求在ＣＰＵ／ ＧＰＵ异 构集群平台上实现时的技术瓶颈，重点介绍了基于ＧＰＵ⁃ＲＤＭＡ直接通信技术和ＧＰＵ⁃Ｓｔｒｅａｍ流技术的并行算法，形 成适合于ＧＰＵ集群的大规模分布式并行计算，即高密度采集地震数据的三维叠前逆时偏移计算技术；应用高密度 实际地震数据测试，计算效率可提升１５～４０倍。 关键词：逆时偏移；直接通信；流技术；计算效率；并行算法 － － － 中图分类号：Ｐ６３１．４　 　 　 文献标识码：Ａ　 　 　 文章编号：１０００ ８９１８（２０１５）０５ ０９７８ ０７ 　 　 随着勘探开发程度的不断深入，地震资料处理 存在剧烈变化的地下介质。 该技术具有相位准确、 对偏移成像的要求越来越高。 目前，叠前深度偏移 成像精度高、对介质横向速度变化和高陡倾角适应 是成像精度最高的一类方法，在工业界应用最多的 性强，甚至可以利用回转波、多次波正确成像等优 主要有基于射线理论的Ｋｉｒｃｈｈｏｆｆ 积分法和基于单 点，是目前理论最先进、成像精度最高的地震偏移成 程波动理论的波动方程两类叠前深度偏移方法。 像方法。 由于该技术计算量、存储需求量巨大这些 Ｋｉｒｃｈｈｏｆｆ 积分偏移算法的理论是基于波场高频近 瓶颈问题不易解决，导致 ＲＴＭ一直没有应用于生 似，而它的实现通常应用了单走时路径的假设，从而 产。 直到近期，随着高性能计算机及其应用技术的 使得该技术不适用于速度横向变化剧烈的介质的高 发展，许多学者对该方法进行了大量卓有成效的研 精度成像；而单程波动方程偏移算法的理论来源于 ［２］ 究工作，使其逐步走 向工业化。 Ｓｙｍｅｓ 提出 对全声波方程的单程波逼近，其有效的数值实现需 ＣｈｅｃｋＰｏｉｎｔｉｎｇ方法，在增加一定计算量的基础上减 要对描述单程波的拟微分方程实行进一步地简化， ［３］ 少波场的存储需求；Ｒｏｂｅｒｔ 提出用随机边界条件 进而导致了该偏移算法存在偏移倾角限制，不适应 代替吸收边界条件，通过增加一次波场反传重建波