边界元法在中梯装置视电阻率二维地形改正中的应用 - 矿产勘查.pdfVIP

第 １卷　第４期 矿 产 勘 查 Ｖｏｌ１　Ｎｏ４ ２０１０年７月 　　　　　　　　　　　　　　　ＭＩＮＥＲＡＬＥＸＰＬＯＲＡＴＩＯＮ　　　　　　　　　　　　　　　 Ｊｕｌｙ，２０１０ 边界元法在中梯装置视电阻率二维 地形改正中的应用 梁芳敏，张小路，匡应平，魏继祖 （桂林理工大学地球科学学院，桂林　５４１００４） 摘　要　文章应用边界单元法进行中间梯度装置视电阻率法二维任意地形改正。在此过程中，应用边界单元法数 学理论基础，编制求解方程组的计算机程序，并计算了纯地形的视电阻率。通过与模拟实验结果对比，证 明其正确性，然后将实际项目的实测电阻率与其纯地形视电阻率进行比值改正，有效地解决实际地球物 理问题，验证了方法的实用性和有效性。 关键词　电阻率法勘探　视电阻率　地形改正　奇异值分解　边界单元法 中图分类号：　Ｐ６３１．３　文献标识码：Ａ　　　文章编号：１６７４－７８０１（２０１０）０４－０３８９－０４ 方法。 ０　引言 ［１］ 根据前人 的模拟设计及理论推导得到的线 在电法勘探的实际工作中，常常受到地形因素 性代数方程组： 的干扰影响，如在山区及丘陵地带的起伏地形等，使 Ｗｉ １ －ＦＵ＝ Ｕ （１１） 视电阻率曲线发生畸变。这些干扰使视电阻率曲线 ( ２π ) ２ ０ 　　 其中Ｗ ＝ｄｉａｇ（ ），Ｆ＝（Ｆ），Ｕ＝（ｕ）， 变化复杂，或者使异常曲线特征点变得不明显。由 ｉ ωｉ ｉｊ ｉ 于地形影响因素的存在，直接影响和降低了电法勘 Ｕ ＝（ｕ）；ｉ，ｊ＝１，···，ｎ ０ ０ｉ 探的地质效果，忽略这些影响将使推断解释工作难 方程（１１）含有 ｎ个方程和 ｎ个未知数 ｕ， １ 以进行，甚至导致错误的地质结论。因此，弄清楚地 ···，ｕ．用奇异值分解法解此线性代数方程组。通 ｎ 形对电剖面 异常曲线的影响，并采取有效的措施 过编制 ＶＣ＋＋语言程序求解，得到边界上各节的函 ρ ｓ 即地形改正，对电法实测数据进行预处理是电法资 数值；通过求得各点的电位值，再根据电阻率装置类 料解释前非常重要的工作。边界单元法在地形改正 Δｕ 型、装置大小，按视电祖率公式 ＝Ｋ ，逐点计算 ρ ｓ Ｉ 问题上可发挥其优越性（徐世浙、楼云菊等，１９８５；