基于多线程的地震相干体属性提取算法
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
计 算 机 系 统 应 用
ｈｔ：ｗ ． Ｓ ．ｒ． ｔ ／ ｗｗｃ － ｏｇｃ ｐ／ —ａ ｎ
２１ 年 第 ２ 卷 第 １ ０２ ｌ １期
基 于 多线程 的地震相 干体属 性提取算法①
杨 尚琴， 自龙， 许 洪承煜
（ 中国石油化工股份有限公司 石油物探技术研究院， 京 ２ １０ ） 南 １１３
摘
要 ：为了充分发挥计算机的多核优势， 高地震数据相干体 的计算速度，通过研 究在多核上的多线程并行技 提
术，完成 了并行相干体算法 的设计 与实现，并分别对 串行和并行算法进行性 能 比较 测试． 测试 结果表 明： ｔｒａ Ｐｈｅｄ 多线程 技术可 以充 分利用 多核资源，取得 比较理想 的线性加速 比，且提 高了系统 的计算 效率，非常适合于大数据
ｒ ｓ ｌ ａｓ ｄ ｃ ｔｓｔ ａ ｉｇ Ｐｔｒ ａ ａ ｍｐ ｏ ｅｔ ｅ ｃ ｍｐ ｔｔｏ ａ ｆ ｃ ｅ ｃ ｆｔ ｅ ｓ ｔｍ ｎ ｅ ｖ ｒ ｕ ｔｂ ｅ ｆ ｒ ｅ ｕ ｔ ｌｏ ｉ ｉａ ｅ ｔ ｕｓ ｈ ｅ ｄ Ｃ ｉ ｒ ｖ ｏ ｕ ｉｎ ｌｅ ｉｎ ｙ ｏ ｈ ｙｓｅ ａ ｄ ｂ ｅｙ ｓ ｉ ｌ ｏ ｎ ｈ ｎ ｎ ｈ ａ ｉ ａ ｐ ｏ ｅ ｓ ｇ ｌ ｇ －ｃ ｅｓ ｉｍｉ ａ ． ｒ ｃ ｓ ｉ ｅ ｓ ａ ｅｓ ｃｄ ｔ ｎ ａ ｒ ｌ ａ Ｋｅ ｒ ｓ ｃ ｈ ｒ ｎ ｅ ｓ ｉｍｉ ａａｐ ｏ ｅ ｓｎ ； ｔ ｒ ａ ｍｕ ｔ－ｈ ｅ ｄ ｇ ｐ ｒ ｌ ｌ ｏ ｕａ ｉｎ ｙｗｏ ｄ ： ｏ ｅ ｅ ｃ ； ｅｓ ｃｄ ｔ ｒ ｃ ｓｉ ｇ Ｐ ｈ ｅ ｄ； ｌｉ ｒ ａ ｉ ； ａ ａ ｌ ｍｐ ｔｔｏ ｔ ｎ ｅｃ
ｉ ｐ ｅ ｅ ｔｄ． ｅｃ ｍ ｐ ｒｓ ｎｔ ｓ ｆｔ ｅｐ ｒｏ ｍ ａ ｃ ｅｗｅ ｎ ｔ ｅｓ ｒａ ｎ ｈ ａ ａｌｌ ｌ ｏ ｉ ｍ ａ ｅ ｏ ｐ ｅｅ ｍ ｌｍ ｎ ｅ Ｔｈ ｏ ａｉｏ ｔｏ ｅ ｒ ｎ ｅｂ ｔ ｅ ｈ ｅｉ ｌ ｄ ｔ ｅｐ ｒ ｌｅ ｇ ｒｔ ｅ ｈ ｆ ａ ａ ｈ ｈ ｓｂｅ ｎ ｃ ｍ ｌｔｄ．
Ｔ ｅｒ ｓ ｌ ｓ ｏ ａ ｓ ｇ Ｐ ｈ ｅ ｄＣ ｋ ｌ ｕ ｅ ｏ ｌ － ｏ ｅｒ ｓ ｕ ｃ ｓａ ｄ ａ ｈ ｅ ｅ ａ ｅ ｌ ｉｅ ｐ ｅ ｕ ． ｈ ｈ ｕ ｔ ｈ ｗｓ ｈ ｔ ｉ ｔ ｒ ａ ａ ｍａ ｅ ｆ ｌ ｓ ｆ ｅ ｔ ｕ ｎ ｎ ｕ ｍｕｔ ｃ ｒ ｅ ｏ ｒ ｅ ｃ ｉｖ ｉ ａ ｌ ａ ｓ ｅ ｄ ｐ Ｔ ｅ ｉ ｎ ｎ ｄ ｎ ｒ
意【 ．
的 Ｃ ＤＡ技术． Ｃ Ｄ Ｕ ｎ） Ｕ Ａ并行利用的 Ｇ Ｕ通常在显卡 Ｐ 上，需要 Ｎ Ｄ Ａ等专 门的显卡支持， ＶＩ Ｉ 代价 比较高，且 目前 Ｃ Ｄ 并行计 算的数据要经过 内存才 能加载到 Ｕ Ａ
ＧＰ 考 虑到本算法 的特 定性—— 没有大规模矩 阵数 Ｕ，
值运算，不能很好的解决 Ｉ ／ ０与计算并行问题，本文不 讨论此方案． ２Ｏ ｅＭＰ 通过对原有 的串行代码插入 （） ｐ ｎ
一
些指 导性 的注释，并进行 必要的修 改，可 以快速 的
随着物探 技术的发展，地震属 性技术在地震 资料 解释 中的作用越来越 大，例如在储层预测和表征方面，
利用地 震属性技术可 以预 测岩性和有利Fra Baidu bibliotek集 体， 描述
术，并根据现有 硬件 的特 性，充分利用 多核计算机 来
帮助我们完 成工作 ．目前 ， 以使用 的并行编程 技术 可
主要包括针对 Ｃ Ｕ的 Ｏ ｅＭＰ Ｐｈｅｄ Ｐ ｐｎ 、 ｔ ａ 等和针对 ＧＰ ｒ Ｕ
Ｙ ＡＮＧ Ｓ ａ ｇＱｉ， ｈ ｎ － ｎ ＸＵ ｉ ｏ ｇ ＨＯＮＧ Ｃｈ ｎ － ｕ Ｚ－ ｎ ， Ｌ ｅｇＹ （Ｉ Ｐ ＣＧ ｏ ｈ ｃｌ ｓａｃ ｓｔｔ， ｎｉｇ １１３ Ｃ ｉａ ＳＮＯ Ｅ ｅｐ ｙ ｉ ｅｒｈ ｎｔｕｅＮａｊ １０ ， ｈｎ ） ｓ ａＲｅ Ｉ ｉ ｎ２
Ａｂｓｒ ｃ ：Ｆ ｒ ｔｋｎｇ ｆ ｌ ｄ ａ ｔｇ ｆ ｍｕｔ－ ｏ ｅ ｃ ｍｐ ｔｒ ｄ ｉ ｒ ｖｎ ｏ ｕ ｉｇ ｓ ｅ ｆｃ ｈ ｒ ｎ ｕ ｆ ｔ ａ ｔ ｏ ａ ｉ ｕｌ ａ ｖ ｎａ ｅ ｏ ｌ ｃ ｒ ｏ ｕ ｅ ｓ ａ ｍｐ ｏ ｉｇ ｃ ｍｐ ｔ ｐｅ ｄ ｏ ｏ ｅ ｅ ｔｃ ｂｅ ｏ ｉ ｎ ｎ ｓ ｉｍｉ ｔ，ａｃ ｈ ｒ ｎ ｅ ｐ ｒ ｌ ｌａｇ ｒｔ ｍ ｔ ｈ ｅ ｅ ｒ ｈ ｏｆｍｕｌ － ｏ ｅｐ ａ ｌｌ ｅ ｈ ｌ ｇ ａ ｅ ｎ ｄ ｓｇ ｅ ｄ ｅｓ ｃｄａａ ｏ ｅｅ ｃ ａａ ｌ ｌ ｏ ｉｈ ｗｉ ｔ ｅ ｒ ｓ ａ ｃ ｔ ｃ ｒ ａ ｌ ｃ ｎｏ ｏ ｙ ｈ ｓｂ ｅ ｅ ｉ ｎ ｄ ａ ｅ ｈ ｉ ｒ ｅｔ ｎ
油藏特 征及 孔隙度变化 ， 找死油 区， 至于监测 流 寻 甚 体运 动和进行其他综 合研 究． 而地 震数据相干体 属性 是众 多地震属性 中相 当重要又常用 的一种，它是通 过 对三 维数据体 的不连 续性进 行分析，可识别构造和 断 层 的分布，使解释人员在解释之 前就能获得研 究区域 粗 略的构造几 何形态 和断层分布情况 ， 免解释 的随 避
量的地震数据 处理 的应用．
关键 词：相干体；地震数据处理； ｔ ｅｄ Ｐｈ ａ；多线程技术；并行计算 ｒ
Ａｔｒｂｕｔ ｔ ａ ｔｏ Ａｌ ｏ ｉｈ ｆＳ ｉｍ ｉ ｈｅ ｅ ｃ ｓ ｄ ｏ Ｍ ｕｌｉＴｈｒ ａ ｎｇ ｔ ｉ ｅＥｘ ｒ ｃ ｉ ｎ ｇ ｒ ｔ ｍ ｏ ｅ ｓ ｃ Ｃｏ ｒ ｎ ｅＢａ ｅ ｎ ｔ— ｅ ｄｉ
ｈｔ：ｗ ． Ｓ ．ｒ． ｔ ／ ｗｗｃ － ｏｇｃ ｐ／ —ａ ｎ
２１ 年 第 ２ 卷 第 １ ０２ ｌ １期
基 于 多线程 的地震相 干体属 性提取算法①
杨 尚琴， 自龙， 许 洪承煜
（ 中国石油化工股份有限公司 石油物探技术研究院， 京 ２ １０ ） 南 １１３
摘
要 ：为了充分发挥计算机的多核优势， 高地震数据相干体 的计算速度，通过研 究在多核上的多线程并行技 提
术，完成 了并行相干体算法 的设计 与实现，并分别对 串行和并行算法进行性 能 比较 测试． 测试 结果表 明： ｔｒａ Ｐｈｅｄ 多线程 技术可 以充 分利用 多核资源，取得 比较理想 的线性加速 比，且提 高了系统 的计算 效率，非常适合于大数据
ｒ ｓ ｌ ａｓ ｄ ｃ ｔｓｔ ａ ｉｇ Ｐｔｒ ａ ａ ｍｐ ｏ ｅｔ ｅ ｃ ｍｐ ｔｔｏ ａ ｆ ｃ ｅ ｃ ｆｔ ｅ ｓ ｔｍ ｎ ｅ ｖ ｒ ｕ ｔｂ ｅ ｆ ｒ ｅ ｕ ｔ ｌｏ ｉ ｉａ ｅ ｔ ｕｓ ｈ ｅ ｄ Ｃ ｉ ｒ ｖ ｏ ｕ ｉｎ ｌｅ ｉｎ ｙ ｏ ｈ ｙｓｅ ａ ｄ ｂ ｅｙ ｓ ｉ ｌ ｏ ｎ ｈ ｎ ｎ ｈ ａ ｉ ａ ｐ ｏ ｅ ｓ ｇ ｌ ｇ －ｃ ｅｓ ｉｍｉ ａ ． ｒ ｃ ｓ ｉ ｅ ｓ ａ ｅｓ ｃｄ ｔ ｎ ａ ｒ ｌ ａ Ｋｅ ｒ ｓ ｃ ｈ ｒ ｎ ｅ ｓ ｉｍｉ ａａｐ ｏ ｅ ｓｎ ； ｔ ｒ ａ ｍｕ ｔ－ｈ ｅ ｄ ｇ ｐ ｒ ｌ ｌ ｏ ｕａ ｉｎ ｙｗｏ ｄ ： ｏ ｅ ｅ ｃ ； ｅｓ ｃｄ ｔ ｒ ｃ ｓｉ ｇ Ｐ ｈ ｅ ｄ； ｌｉ ｒ ａ ｉ ； ａ ａ ｌ ｍｐ ｔｔｏ ｔ ｎ ｅｃ
ｉ ｐ ｅ ｅ ｔｄ． ｅｃ ｍ ｐ ｒｓ ｎｔ ｓ ｆｔ ｅｐ ｒｏ ｍ ａ ｃ ｅｗｅ ｎ ｔ ｅｓ ｒａ ｎ ｈ ａ ａｌｌ ｌ ｏ ｉ ｍ ａ ｅ ｏ ｐ ｅｅ ｍ ｌｍ ｎ ｅ Ｔｈ ｏ ａｉｏ ｔｏ ｅ ｒ ｎ ｅｂ ｔ ｅ ｈ ｅｉ ｌ ｄ ｔ ｅｐ ｒ ｌｅ ｇ ｒｔ ｅ ｈ ｆ ａ ａ ｈ ｈ ｓｂｅ ｎ ｃ ｍ ｌｔｄ．
Ｔ ｅｒ ｓ ｌ ｓ ｏ ａ ｓ ｇ Ｐ ｈ ｅ ｄＣ ｋ ｌ ｕ ｅ ｏ ｌ － ｏ ｅｒ ｓ ｕ ｃ ｓａ ｄ ａ ｈ ｅ ｅ ａ ｅ ｌ ｉｅ ｐ ｅ ｕ ． ｈ ｈ ｕ ｔ ｈ ｗｓ ｈ ｔ ｉ ｔ ｒ ａ ａ ｍａ ｅ ｆ ｌ ｓ ｆ ｅ ｔ ｕ ｎ ｎ ｕ ｍｕｔ ｃ ｒ ｅ ｏ ｒ ｅ ｃ ｉｖ ｉ ａ ｌ ａ ｓ ｅ ｄ ｐ Ｔ ｅ ｉ ｎ ｎ ｄ ｎ ｒ
意【 ．
的 Ｃ ＤＡ技术． Ｃ Ｄ Ｕ ｎ） Ｕ Ａ并行利用的 Ｇ Ｕ通常在显卡 Ｐ 上，需要 Ｎ Ｄ Ａ等专 门的显卡支持， ＶＩ Ｉ 代价 比较高，且 目前 Ｃ Ｄ 并行计 算的数据要经过 内存才 能加载到 Ｕ Ａ
ＧＰ 考 虑到本算法 的特 定性—— 没有大规模矩 阵数 Ｕ，
值运算，不能很好的解决 Ｉ ／ ０与计算并行问题，本文不 讨论此方案． ２Ｏ ｅＭＰ 通过对原有 的串行代码插入 （） ｐ ｎ
一
些指 导性 的注释，并进行 必要的修 改，可 以快速 的
随着物探 技术的发展，地震属 性技术在地震 资料 解释 中的作用越来越 大，例如在储层预测和表征方面，
利用地 震属性技术可 以预 测岩性和有利Fra Baidu bibliotek集 体， 描述
术，并根据现有 硬件 的特 性，充分利用 多核计算机 来
帮助我们完 成工作 ．目前 ， 以使用 的并行编程 技术 可
主要包括针对 Ｃ Ｕ的 Ｏ ｅＭＰ Ｐｈｅｄ Ｐ ｐｎ 、 ｔ ａ 等和针对 ＧＰ ｒ Ｕ
Ｙ ＡＮＧ Ｓ ａ ｇＱｉ， ｈ ｎ － ｎ ＸＵ ｉ ｏ ｇ ＨＯＮＧ Ｃｈ ｎ － ｕ Ｚ－ ｎ ， Ｌ ｅｇＹ （Ｉ Ｐ ＣＧ ｏ ｈ ｃｌ ｓａｃ ｓｔｔ， ｎｉｇ １１３ Ｃ ｉａ ＳＮＯ Ｅ ｅｐ ｙ ｉ ｅｒｈ ｎｔｕｅＮａｊ １０ ， ｈｎ ） ｓ ａＲｅ Ｉ ｉ ｎ２
Ａｂｓｒ ｃ ：Ｆ ｒ ｔｋｎｇ ｆ ｌ ｄ ａ ｔｇ ｆ ｍｕｔ－ ｏ ｅ ｃ ｍｐ ｔｒ ｄ ｉ ｒ ｖｎ ｏ ｕ ｉｇ ｓ ｅ ｆｃ ｈ ｒ ｎ ｕ ｆ ｔ ａ ｔ ｏ ａ ｉ ｕｌ ａ ｖ ｎａ ｅ ｏ ｌ ｃ ｒ ｏ ｕ ｅ ｓ ａ ｍｐ ｏ ｉｇ ｃ ｍｐ ｔ ｐｅ ｄ ｏ ｏ ｅ ｅ ｔｃ ｂｅ ｏ ｉ ｎ ｎ ｓ ｉｍｉ ｔ，ａｃ ｈ ｒ ｎ ｅ ｐ ｒ ｌ ｌａｇ ｒｔ ｍ ｔ ｈ ｅ ｅ ｒ ｈ ｏｆｍｕｌ － ｏ ｅｐ ａ ｌｌ ｅ ｈ ｌ ｇ ａ ｅ ｎ ｄ ｓｇ ｅ ｄ ｅｓ ｃｄａａ ｏ ｅｅ ｃ ａａ ｌ ｌ ｏ ｉｈ ｗｉ ｔ ｅ ｒ ｓ ａ ｃ ｔ ｃ ｒ ａ ｌ ｃ ｎｏ ｏ ｙ ｈ ｓｂ ｅ ｅ ｉ ｎ ｄ ａ ｅ ｈ ｉ ｒ ｅｔ ｎ
油藏特 征及 孔隙度变化 ， 找死油 区， 至于监测 流 寻 甚 体运 动和进行其他综 合研 究． 而地 震数据相干体 属性 是众 多地震属性 中相 当重要又常用 的一种，它是通 过 对三 维数据体 的不连 续性进 行分析，可识别构造和 断 层 的分布，使解释人员在解释之 前就能获得研 究区域 粗 略的构造几 何形态 和断层分布情况 ， 免解释 的随 避
量的地震数据 处理 的应用．
关键 词：相干体；地震数据处理； ｔ ｅｄ Ｐｈ ａ；多线程技术；并行计算 ｒ
Ａｔｒｂｕｔ ｔ ａ ｔｏ Ａｌ ｏ ｉｈ ｆＳ ｉｍ ｉ ｈｅ ｅ ｃ ｓ ｄ ｏ Ｍ ｕｌｉＴｈｒ ａ ｎｇ ｔ ｉ ｅＥｘ ｒ ｃ ｉ ｎ ｇ ｒ ｔ ｍ ｏ ｅ ｓ ｃ Ｃｏ ｒ ｎ ｅＢａ ｅ ｎ ｔ— ｅ ｄｉ