矩形棱柱体重力梯度张量异常正演计算公式-物探与化探.PDFVIP

　 第３９卷第６期 物　 探　 与　 化　 探 Ｖｏｌ．３９，Ｎｏ．６　 　 ２０１５年１２月 ＧＥＯＰＨＹＳＩＣＡＬ ＆ ＧＥＯＣＨＥＭＩＣＡＬ ＥＸＰＬＯＲＡＴＩＯＮ Ｄｅｃ．，２０１５　 ｄｏｉ：１０．１１７２０／ ｗｔｙｈｔ．２０１５．６．２０ － 舒晴，朱晓颖，周坚鑫，等．矩形棱柱体重力梯度张量异常正演计算公式［Ｊ］．物探与化探，２０１５，３９（６）：１２１７ １２２２．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０．１１７２０／ ｗｔｙｈｔ． ２０１５．６．２０ Ｓｈｕ Ｑ，ＺｈｕＸＹ，ＺｈｏｕＪＸ，ｅｔａｌ．Ｆｏｒｗａｒｄｍｏｄｅｌｉｎｇｅｘｐｒｅｓｓｉｏｎｓｆｏｒｒｉｇｈｔｒｅｃｔａｎｇｕｌａｒｐｒｉｓｍｗｉｔｈｃｏｎｓｔａｎｔｄｅｎｓｉｔｙ［Ｊ］．ＧｅｏｐｈｙｓｉｃａｌａｎｄＧｅｏｃｈｅｍｉｃａｌＥｘｐｌｏｒａ⁃ － ｔｉｏｎ，２０１５，３９（６）：１２１７ １２２２．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０．１１７２０／ ｗｔｙｈｔ．２０１５．６．２０ 矩形棱柱体重力梯度张量异常正演计算公式 １，２ ２ ２ ２ ２ 舒晴 ，朱晓颖 ，周坚鑫 ，高维 ，尹航 （１．吉林大学地球探测科学与技术学院，吉林 长春　 １３００２１；２．中国国土资源航空物探遥感中心， 北京　 １０００８３） 摘 要：矩形棱柱体重力梯度张量异常常用的正演计算公式含有解析奇点。 结合理论推导过程详细分析了“奇点” 存在的原因，并基于前人提出的长方体ΔＴ场及其梯度场无解析奇点理论表达式，建立了矩形棱柱体重力梯度张量 无解析“奇点”公式，通过模型计算结果对比和理论分析证明了其正确性。 关键词：重力梯度张量异常；矩形棱柱体；解析奇点 － － － 中图分类号：Ｐ６３２　 　 　 文献标识码：Ａ　 　 　 文章编号：１０００ ８９１８（２０１５）０６ １２１７ ０６ 　 　 矩形棱柱体作为最常用的模型构建单元被广泛 应用于重力梯度正演数值模拟之中。 笔者研究发 现，目前被广泛引用的矩形棱柱体重力梯度张量异 － ［１ ５］ ［６］ 常计算公式 均源自Ｄｏｎａｌｄ Ｐｌｏｕｆｆ ，该公式虽然 简洁、对称，但文献中仅给出了结果，没有详细地推 ［７］ 导过程，且存在解析奇点 ，在实际应用中需要对 观测平面进行特殊剖分。 因此，通过详细推导，剖析 解析奇点存在原因，并寻找新途径构建新的无解析 奇点公式显得十分必要。 图１　 矩形棱柱体与外部场点相互关系示意 １　 问题的提出 ２ ２ ２ ＝ － ＋ － ＋ － ＝ × 其中，ｒ （ξ ｘ） （η ｙ） （ζ ｚ） ，Ｇ ６．６７２ 在众多重力梯度张量分量正演数值模拟和反演