矿井直流电法三维正演计算的若干问题.pdfVIP

! 第# 卷第 期 物! 探! 与! 化! 探 $%’ # ，(%’ ! ! ))* 年* 月 +,-./012345 6 +,-3/,72345 ,8.5-94:2-( 4;’ ，))*! 矿井直流电法三维正演计算的若干问题 刘树才，刘志新，姜志海，岳建华 （中国矿业大学，江苏徐州! =))# ） 摘要：介绍了用于矿井直流电法三维正演模拟不同数值计算方法的特点及其适用范围，提出了用对数等间距和算 术等间距相结合的网格剖分方式，说明了电极互换原理所满足的地质条件。针对合成系数矩阵为非对角占优的特 点，采用了大数处理方法，使其线性方程组适用于现有的求解方法。最后通过对具体问题的探讨，阐述了地质模型 的计算精度和计算时间之间的相互关系，为模型参数的选择提供了一定的理论基础。 关键词：矿井直流电法；正演模拟；网格剖分；矩阵压缩 中图分类号： .?=’ ?! ! ! 文献标识码：4! ! ! 文章编号：=))) @ #A=#（))* ）) @ )=B) @ )? ! ! 矿井直流电法勘探中经常遇到三维地电构造， 数值计算方法。其基本原理是：首先把矿井直流电 大多为层状介质中含有三维异常体情况，目前多用 法所满足偏微分方程的边值问题，转化为相应泛函 一维、二维数值计算方法来进行分析，但解释结果与 的极值问题（变分问题）；再经过离散化，得到由空 实际地质模型有较大的差异。而三维空间能够真实 间各点未知电位值组成的高阶线性方程组；最后求 地反映实际异常体形态的变化情况，能更精确地反 解该方程组，可确定各点的电位分布。该方法适合 映地下电场的分布规律。三维数值模拟在矿井地质 于复杂矿井地质模型的计算（如顶底板内的局部岩 条件下的技术要点为边界条件设置、研究区域网格 溶、小断层或煤层内陷落柱等），只要适当地选择网 剖分及节点编码技术、系数矩阵压缩存储技术、大型 格单元体的类型及剖分步长的大小，即可使数值解 线性方程组的解法等。下面对矿井直流电法勘探正 达到较高的计算精度。 演计算中所面临的一些关键问题进行探讨。 积分方程法是从电场基本定律及场所满足的微 分方程边值问题出发，求得有关参数（电位! 或积 =! 计算方法的选择 累电荷密度 ），来满足积分方程，然后用近似计算 目前矿井直流电法勘探三维数值模拟技术主要 方法求此积分方程的数值解。其基本原理是：先将 包括三维有限差分法、三维有限单元法、积分方程法 三维地质模型表面用小面单元（三角形或矩形单 和边界单元法等。 元）进行离散，并设小单元内电荷密度 为常数；由 有限差分法又称网格法，其计算过程为：首先将 此将电荷密度满足的方程离散化，形成线性方程组， 求解区域离散成许多小正方形或长方形的网格，并 解此方程组即可得模型表面小单元的积累电荷密 以网格节点上的参数值来表征电场的空间分布，然 度；然后，对所有小单元在某一点# 处产生的电位 后用网格节点上电位函数的差商来近似代替该点的 求和，便得到表面积累电荷产生的异常场。该方法 偏导数（或微商），由此得到一个关于网格节点电位 易于处理三维局部异常的矿井地质问题。 值的高阶线型方程组，解出此方程组就可算出网格 边界单元法的基本思想与积分方程法相似，把 节点上的场参数值。有限差分法计算原理和程序设 研究区域的边界划分成一系列小单元，通过边界积 计比较简单，易于解决由二维过渡到三维地球物理