基于COMSOL掘进巷道直流电法超前探测研究_郝耀军.pdf

煤矿现代化2020 年第 6 期总第 159 期 0引言 煤矿巷道掘进时， 如果存在水害威胁， 一定要采 取先探测后采掘的措施。 当异常体存在于巷道迎头前 方时， 要进行超前探测来发现异常体。矿井直流电法 超前探测技术与一般的电测深的原理有所不同， 因此 所选择的装置也有所不同， 直流电法超前探测所选的 测量装置一般是单点供电的三极 AMN 装置， 该装置 在超前探测中对含水异常体有较好的探测效果和较 高的异常分辨率。 本文通过有限元分析软件进行数值 模拟与生产实例相结合， 总结出巷道掘进迎头前方地 质异常体对测量结果的影响规律， 为矿井直流电法超 前探测资料的解释提供理论依据， 也能够对实际测量 中采集数据的反演解释工作提供指导[1]。 1勘探方法原理 矿井直流电法与地面电阻率法原理相同， 虽然选 用的测量装置不尽相同， 都会有供电电极对大地进行 供电， 形成稳定电流场， 在供电电极所在的测线 （其余 位置 ） 设置测量电极， 测得电位数据将用于计算地层 的视电阻率。对于矿井直流电法来说， 按照事先设定 好的测线将供电以及测量电极布设在被测巷道底板 或者侧帮上， 通过改变测量电极位置从而达到观测供 电电极形成稳定电流场的分布以及变化特征， 利用上 述的测量数据来超前预测巷道前方的构造异常情况， 从而达到超前探测的目的。 三极 AMN 装置是最早应用于矿井直流电法超 前探测的装置， 如图 1 所示。目前广泛应用于超前探 测巷道迎头前方的地质异常体。 在三维空间介质中利 用单点电源 A 进行供电（另一供电电极 B 设置在相 对无穷远处 ） ， 并且将 M、 N用作测量电势差。超前探 测与电测深原理不同，因此相应的操作方法有所不 同， 超前探测电极布设方法 将供电电极 A 固定到巷 道迎头上， M、 N 测量电极处于测线上，并且位于 A 点后方， 所测电位差为△UMN， 记录点处的电场强度 与供电点 A 和测量电极 M、 N 的距离有很大关系， 记 录点为两测量电极的中点。 在巷道迎头处进行超前探 测时，当巷道的垂直截面尺寸都小于记录点到 A 点 距离时， 巷道对测量结果影响较小； 当巷道的垂直截 面尺寸都大于记录点到测量点距离时， 巷道对测量结 基于 COMSOL 掘进巷道直流电法超前探测研究 郝耀军 ， 李浩 ， 武 艺 （西山煤电集团西曲矿 ， 山西 太原 030200 ） 摘要 煤矿巷道掘进过程中存在突水隐患时， 利用矿井直流电法超前探测技术， 对巷道前方煤岩体 富水性进行探查， 对矿井防治水工作具有指导意义。本文通过 COMSOL Multiphysics 数值模拟与工 程实例相结合的方法研究了矿井直流电法超前探测技术影响因素与视电阻率曲线变化之间的规律， 并利用这些规律指导矿井直流电法超前探测工程实践， 以及提供理论支持。 关键词 直流电法 ； 超前探测 ； 数值模拟 ； 视电阻率 中图分类号 TD166 ; TD745.2文献标识码 A 文章编号 1009-0797 （2020 ） 06-0104-03 Research on Advance Detection of Direct Current in Driving Roadway Based on COMSOL HAO Yao jun , LI Hao , WU Yi （Xi Shan Coal Electricity Group Xi Qu Coal Mine ，Taiyuan 030200 ， China ） Abstract When there are hidden dangers of water inrushed in the process of mine roadway excavation, the water-rich of coal and rock mass in front of the roadway is explored by using the advanced detection technology of mine direct current , which has guiding significance for mine water prevention and control work.In this paper, the law between the influence factors of mine direct current advanced detection technology and the change of apparent resistivity curve is studied by COMSOL Multiphysics numerical simulation and engineering example, and these laws are used to guide the practice of mine direct current advanced detection engineering and provide theoretical support. Keywords direct current （DCM ）； advanced detection ； numerical simulation ； apparent resistivity 104 ChaoXing 煤矿现代化2020 年第 6 期总第 159 期 果影响很大， 当这种影响不可被忽略时， 按照视电阻 率公式计算视电阻率值，然后进行消除巷道影响处 理， 就可以得出沿巷道的视电阻率剖面曲线[3]。 图 1AMN 装置超前探测原理图 2矿井直流电法超前探测异常体数值模拟 2.1有限元数值模拟原理及流程 直流电法正演模拟是将事先设定好的的地电模 型利用数值计算来求解稳定场的分布情况。利用 COMSOL Multiphysics 进行地电模型正演数值模拟过 程， 结合下述的流程图分析， 详细步骤如图 2 所示 1 ）首先根据分析问题需要选择三维空间模型， 选择软件的电流ec模块， 研究状态为稳态。 2 ）进行全局变量设置和几何模型构建，然后全 部构建几何模型形成联合体，根据地电模型需要， 对 几何联合体进行点、 边界和域的定义， 便于对地电模 型各部分赋值等工作。 3 ）根据实际情况，进行地电模型的材料设置， 然后对静电场的点电源进行赋值， 以及接地边界的 定义。 4 ）地电模型的网格剖分，地电模型的主体划分 为自由剖分四面体网格， 无穷远层状介质， 网格剖分 成薄层扫掠形式。 由于建立的模型较大需要对重点区 域进行细致研究， 就要 将重点部位网格细化剖分， 细化的重点为巷道周 围以及异常体周围。供电电极周围网格超细化剖分。 5 ）地点模型数值运算，运算速度与模型划分的 总单元数有关， 模型的单元数越多， 数值计算的速度 越缓慢。 6 ）地电模型的后处理过程，可以根据需要对正 演模拟的数据集进行三维整体输出， 也可以进行截面 数据输出或者截点数据输出。也可以利用 COMSOL Multiphysics 软件进行地电模型的三维切片、 一维、 二 维、 三维图组处理， 然后输出成果图。 图 2COMSOL Multiphysics 数值模拟流程图 2.2多异常体并存对视电阻的影响 在实际生产中， 往往遇到复杂的地层存在， 有可 能是高阻和低阻同时出现， 通过矿井直流电法数值模 拟总结复杂地层下视电阻率曲线变化规律， 对于超前 探测具有较好的指导作用， 具有实际意义。 模型设定围岩长度 1100m，宽度 500m，高度 500m， 供电点 A 在迎头底板处， 巷道长度 500m， 截面 为 6m 的正方形， 巷道电阻率为 108 Ω m， 供电电流 设置为 1A， 电极间距为 10m， 测量装置有三极 AMN， 下面 3 个模型设定了不同异常体存在情况。 图 3两个异常体分布示意图 图 3 为两个异常体共同影响下的电位分布图， 图 中上边的是高阻异常体， 下边是低阻异常体。 模型 1 单独低阻体阻值 10 Ω m， 异常体截面平 行于巷道， 异常体截面为 6060m， 厚度为 6m， 位于 巷道迎头侧前方 1030m处。 模型 2 单独高阻体阻值 3000 Ω m， 异常体截面 平行于巷道， 异常体截面为 6060m， 厚度为 6m， 位 于巷道迎头侧前方 1030m处。 模型 3高阻低阻体共同存在，低阻体阻值 10 Ω m， 高阻体阻值 3000 Ω m， 异常体截面都平行于 巷道， 异常体截面为 6060m， 厚度为 6m， 位于巷道 迎头侧前方 1030m处。 105 ChaoXing 煤矿现代化2020 年第 6 期总第 159 期 图 4多异常体存在时视电阻率曲线对比图 图 4 为不同异常体影响下的视电阻率曲线， 从图 中可以看出 1 ）当巷道侧前方只有低阻体时，视电阻率曲线 所显示的低阻异常较明显； 当巷道侧前方只有高阻体 时， 视电阻率曲线显示的高阻异常幅度不大。 2 ）当巷道侧前方高阻体和低阻体共同存在时， 通过分析三条曲线异常变化趋势可知， 高阻体和低阻 体共同作用时的视电阻率曲线， 基本符合单独高阻和 单独低阻视电阻率曲线异常叠加效果。 3煤矿探测实例 3.1南六盘区地质概况 西曲矿南六采区面积共 1.6km2，区内沟谷枞横、 地形复杂， 小煤窑众多。 区内上组煤有 2.3、 4 煤， 4 煤属不稳定局部可采煤层， 且 4 煤层厚度变化较大。 区内原勘探程度较低， 矿井后续生产衔接任务无法开 展。西曲矿现计划对南六采区下组煤层进行回采， 而 该采区上覆煤层已基本全部被破坏。 南六采区小煤矿 分布较多， 区内上组煤大部分被小窑破坏， 小窑破坏 程度、 巷道位置及巷道内积水情况均不确定， 为保证 煤矿安全生产， 提前做好预测预报工作， 现利用直流 电法对该区域皮带下山巷道进行超前探测。 3.2数据采集及成果解释 本次井下超前探测数据采集使用的是 WJDJ- 4 型高密度电阻率系统， 沿巷道中心布置测线， 电极间 距为 10m， 测量电压为 144V， 供电时间为 500ms， 采 用三极 AMN 装置进行数据采集， 皮带下山巷道内共 布设了 29 个电极，本次井下超前探测工作巷道内停 止了施工并关闭了用电设备，采集的数据较为稳定， 采集数据质量好。 根据直流电法超前探测数值模拟二维视电阻率 平面图可知， 当巷道前方有低阻异常时， 所形成的二 维视电阻率平面图会显示低阻异常， 并且巷道前方的 异常体规模越大，低阻异常范围在图中显示的也越 大。 通过此规律分析解释西曲矿南六采区皮带下山的 二维视电阻率分布情况[4]。 数据处理后得到电阻率色谱平面图见图 5， 平面 图宽度约 50m， 巷道位置如图所示， 实线为已掘巷道， 用正数表示， 虚线为未掘巷道， 用负数表示， 平面图等 值线为地层电阻率的对数值。 图 5南六采区皮带下山巷道超前探测电阻率色谱平面图 迎头前方 85m 范围内地层电阻率基本为高阻显 示， 电阻率对数值大于 1， 说明 85m 范围内地层富水 性不强，巷道能够向前掘进；在迎头前方 85～105m 范围地层电阻率有所下降，电阻率对数值 0.4～1， 说 明该范围内地层裂隙相对发育，富水性有所增强； 在 迎头前方 105～130m 范围地层电阻率最低，电阻率 对数值小于 0.4, 说明该位置范围断裂、裂隙比较发 育， 地层富水性强； 迎头前方 130m 以后范围地层电 阻率为高阻显示， 电阻率对数值大于 1， 说明该位置 岩性变化， 为奥灰反映[3]。后经钻探验证， 与物探结果 基本吻合。 4结论 1 ）当有多个异常体存在巷道迎头前方时，所测 得视电阻率曲线异常是由所有异常体影响的叠加结 果。通过实例验证数值模拟， 利用数值正演模拟形成 二维视电阻率平面图， 总结规律指导反演解释工作。 2 ）直流电法超前探测技术在实际探查工作中， 测量条件虽复杂， 数据规律起伏较大， 但该方法在掘 进前方富水性探测的应用中， 仍有较为明显的效果。 参考文献 [1] 郝耀军, 朱鲁, 高卫富. 基于矿井直流电法的底板岩层富 水性研究[J]. 山东煤炭科技,2016 （05） 136- 138. [2] 牟义. 典型地质异常体电磁法响应特征研究 [J]. 煤矿开 采,2017,22 （04） 4- 9. [3] 张三敏,王娟,李增涛,魏晓航. 瞬变电磁法探测煤矿采空 区的应用研究 [J]. 工程地球物理学报,2014,11（05） 683- 687. 作者简介 郝耀军 （1989-） ， 山西太原人， 硕士研究生， 助理工程师， 主要从事矿井电法勘探及其应用。 （收稿日期 2020- 2- 26） 106 ChaoXing