三维电场直流电阻率边界单元法地形改正.pdfVIP

第 卷 第 期 ３６ ５ Ｖｏｌ．３６ Ｎｏ．５ ＿ ＿ 年 月 ＣＯMＰＵＴＩＮＧ ＴＥＣＨＮＩＱＵＥＳ ＦＯＲ ＧＥＯＰＨＹＳＩＣＡＬ ＡＮＤ ＧＥＯＣＨＥMＩＣＡＬ ＥＸＰＬＯＲＡＴＩＯＮ ２０１４ ９ Ｓｅ ｔ．２０１４ ｐ ＿ 文章编号： （ ） ＿ １００１-１７４９ ２０１４ ０５-０５２９-０６ 三维电场直流电阻率边界单元法地形改正 ， １ ３ １ ２ ２ ３ 段长生 ，王绪本 ，吕玉增 ，李长伟 ，黄修保 （ 成都理工大学，成都 ； 桂林理工大学，桂林 ； １． ６１００５９ ２． ５４１００４ ＿ ＿ ３．赣中南地质矿产勘查研究院，南昌 ３３００２９） ＿ 摘 要：电法观测易受到地形的影响，给资料解释造成困难。这里对数值计算中地面网格节点 ＿ 的立体角（）公式进行了严格地推导，利用边界单元法对三维直流电场的地形影响进行了计算 ω 并编写了程序。 角域山脊地形算例表明，解析解与数值解的电位及视电阻率相对误差小于 ９０° ，验证了算法及程序的正确性。实例计算了 山谷地形 观测剖面的地形效应，当旁侧剖 ５％ ３Ｄ ２Ｄ 面与主剖面距离大于２４ ｍ 时，地形效应最大值小于５％。 关键词：边界单元法；地形改正；立体角；球面三角形 中图分类号： ＋ 文献标志码： ： ／ Ｐ ６３１．３ ２２ Ａ DOI １０．３９６９ ．ｉｓｓｎ．１００１-１７４９．２０１４．０５．０３ ＿ ＿ ｊ 复杂性，这些方法的应用受到一定限制。随着计算 ０ 引言 机技术的发展和地球物理解释方法的进步，众多学 ＿ 者研究了适合于起伏地形的 、 数值模拟和解 ２Ｄ ３Ｄ 世纪 年代以来，电法勘探在矿产资源勘 ２０ ５０ 释技术，从而使地形影响问题的模拟、识别和校正成 查和工程勘察等方面的应用越来越广泛，但地形起 ［］６ 为可能 。主要的数值计算方法有：积分方程法、边 伏较大的地区，不仅采集难度加大，资料畸变也严 界元方法、有限差分法和有限元方法等。其中