大地电磁三维正反演多核并行计算的设计与实现-物探与化探.PDFVIP

　 第３８卷第３期 物　 探　 与　 化　 探 Ｖｏｌ．３８，Ｎｏ．３　 　 ２０１４年６月 ＧＥＯＰＨＹＳＩＣＡＬ ＆ ＧＥＯＣＨＥＭＩＣＡＬ ＥＸＰＬＯＲＡＴＩＯＮ Ｊｕｎ．，２０１４　 ＤＯＩ：１０．１１７２０／ ｗｔｙｈｔ．２０１４．３．３３ － 顾观文，吴文鹂，梁萌．大地电磁三维正、反演多核并行计算的设计与实现［Ｊ］．物探与化探，２０１４，３８（３）：６０１ ６０６．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０．１１７２０／ ｗｔｙｈｔ． ２０１４．３．３３ Ｇｕ Ｇｕａｎｗｅｎ，ＷｕＷｅｎｌｉ，ＬｉａｎｇＭｅｎｇ．Ｔｈｅｍｕｌｔｉ⁃ｃｏｒｅｐａｒａｌｌｅｌｅｄｓｔｒａｔｅｇｉｅｓａｎｄｅｆｆｉｃｉｅｎｃｙａｎａｌｙｓｉｓｏｆ３Ｄｍａｇｎｅｔｏｔｅｌｌｕｒｉｃｆｏｒｗａｒｄｍｏｄｅｌｉｎｇａｎｄｉｎｖｅｒｓｉｏｎ［Ｊ］． － Ｇｅｏｐｈｙｓｉｃａｌ ａｎｄ Ｇｅｏｃｈｅｍｉｃａｌ Ｅｘｐｌｏｒａｔｉｏｎ，２０１４，３８（３）：６０１ ６０６．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０．１１７２０／ ｗｔｙｈｔ．２０１４．３．３３ 大地电磁三维正、反演多核并行计算的设计与实现 顾观文，吴文鹂，梁萌 （中国地质科学院地球物理地球化学勘查研究所，河北 廊坊　 ０６５０００） 摘 要：首先以频点为并行粒度，完成了ＭＴ三维正演并行计算，然后将该算法引入共轭梯度法ＭＴ三维反演过程中 的正演求数据的残差、反演方程右端项的求取、共轭梯度求解模型增量这三个计算环节中，使得ＭＴ三维反演中的 主要计算量实现并行化。 对编制的ＭＴ三维正、反演并行计算程序，在多核工作站上通过理论地电模型进行了正演 试算和反演验证，并在串行和并行两种计算方式下进行了计算效率对比。 关键词：大地电磁；多核并行；三维正演；三维反演 － － － 中图分类号：Ｐ６３１　 　 　 文献标识码：Ａ　 　 　 文章编号：１０００ ８９１８（２０１４）０３ ０６０１ ０６ 　 　 实现大地电磁法三维反演主要采用两种方式： 旋度，再将式（２）代入到式（１），则将电磁场满足的 一种方式是采用节省内存和节约时间的串行计算方 一阶微分方程变为二阶微分方程，可以得到电场Ｅ 法；另一种方式是采用性能高、内存大的小型计算 满足的方程 － 机、大型计算机或集群计算机的并行算法［１ １１］。 基 æ １ ö × ç × ÷ ＝ ∇ ∇ Ｅ σＥ ， （５） 于串行计算的大地电磁三维正、反演，虽然可能节省 èｉωμ ø 内存和节约时间，但三维大地电磁数据量巨大和地 求解这个偏微分方程组，得到电场分量Ｅ 、Ｅ 、Ｅ 的ｘ ｙ ｚ 质体本身的三维复杂结构，对计算资源的需求很高 分布后，可以根据式（１）求得磁场分量Ｈ 、Ｈ 、Ｈ 的ｘ ｙ ｚ 而基于并行运算的大地电磁三维正反演，应用了微 分布。 机和工作站多核计算技术，使得存储效率和计算速