矿井瞬变电磁法参数优化试验及超前探测应用.pdf
第 48 卷第 6 期煤 炭 科 学 技 术Vol 48 No 6 2020 年6 月Coal Science and Technology Jun.2020 移动扫码阅读 牟 义ꎬ李江华ꎬ徐 慧ꎬ等.矿井瞬变电磁法参数优化试验及超前探测应用[J].煤炭科学技术ꎬ2020ꎬ486 184-190.doi10 13199/ j cnki cst 2020 06 024 MU Yiꎬ LI Jianghuaꎬ XU Huiꎬet al.Parameters optimization test of mine transient electromagnetic method and appli ̄ cation of advanced detection [ J] Coal Science and Technologyꎬ 2020ꎬ 48 6 184 - 190. doi 10 13199/ j cnki cst 2020 06 024 矿井瞬变电磁法参数优化试验及超前探测应用 牟 义1ꎬ2ꎬ3ꎬ李江华1ꎬ2ꎬ3ꎬ徐 慧1ꎬ2ꎬ3ꎬ马志超1ꎬ2ꎬ3ꎬ崔金亮1ꎬ2ꎬ3ꎬ张永超1ꎬ2ꎬ3 1.煤炭科学技术研究院有限公司 安全分院ꎬ北京 100013ꎻ2.煤炭科学研究总院 煤炭资源高效开采与洁净利用国家重点实验室ꎬ 北京 100013ꎻ3. 北京市煤矿安全工程技术研究中心ꎬ北京 100013 摘 要为了精准超前探测采空区及积水区ꎬ确保煤矿的安全掘进ꎬ通过对已知积水采空区进行矿井 瞬变电磁法探测对比试验ꎬ分析不同参数对积水采空区电场响应特征的影响ꎬ优化了一套适合积水采 空区探测的发射频率、叠加次数以及采集密度等采集参数ꎬ并将该优化参数应用到井下超前探测实践 中ꎬ成功发现前方采空区及积水区ꎮ 结果表明矿井瞬变电磁法采用 6.25 Hz 的发射频率、32 次叠加 以及中密度采集方式ꎬ对采空区及积水区超前探测较好ꎻ将该采集优化参数应用到晋城矿区井下超前 探测ꎬ以王台铺矿超前探测应用为例ꎬ发现掘进工作面前方存在明显的 1 处高阻异常体和 2 处低阻异 常体ꎬ且分界面明显ꎬ其与查明的小煤窑巷道分布区基本重合ꎬ后经过打钻验证ꎬ正前方 50 m 之外的 高阻异常区见空ꎬ通过钻孔透视发现ꎬ为小煤窑采空区且无积水ꎬ后退 30 mꎬ向右侧帮 1 处低阻异常 区打钻ꎬ打到深度 35 m 遇空出水ꎬ水量较大ꎬ为采空区积水ꎬ取得较好探测效果ꎻ根据在鄂尔多斯地方 煤矿、临汾地方煤矿、晋煤集团、朔州地方煤矿等井下探测总结的经验ꎬ一般视电阻率位于 30 Ωm 以下推断为含水异常ꎬ30 Ωm 以上推断为不含水正常岩层ꎬ更高则推断为高阻异常ꎮ 关键词超前探测ꎻ采空区ꎻ瞬变电磁法ꎻ发射频率ꎻ叠加次数ꎻ采样密度ꎻ矿井 中图分类号TD166ꎻ P631.3 文献标志码A 文章编号0253-2336202006-0184-07 Parameters optimization test of mine transient electromagnetic method and application of advanced detection MU Yi1ꎬ2ꎬ3ꎬ LI Jianghua1ꎬ2ꎬ3ꎬ XU Hui1ꎬ2ꎬ3ꎬ MA Zhichao1ꎬ2ꎬ3ꎬ CUI Jinliang1ꎬ2ꎬ3ꎬ ZHANG Yongchao1ꎬ2ꎬ3 1 Mine Safety Technology BranchꎬChina Coal Research Instituteꎬ Beijing 100013ꎬChinaꎻ2 State Key Laboratory of Coal Mining and Clean Utilizationꎬ China Coal Research Instituteꎬ Beijing 100013ꎬChinaꎻ3 Beijing Mine Safety Engineering Technology Research Centerꎬ Beijing 100013ꎬChina AbstractIn order to accurately detect the gob and the water-bearing area in advanceꎬ and ensure the safe excavation of the coal mineꎬ the mine transient electromagnetic method was used to detect and compare the test in the known water-bearing gobꎬ and the electric field response characteristics of different parameters to the mined-out gob were analyzed. A set of data acquisition parameters suitable for the detection of the water levelꎬ the number of superpositions and the acquisition density were optimized and the optimized parameters were ap ̄ plied to the practice of underground advance detectionꎬ the front gob and water-bearing gob area were successfully discovered. The results showed that the mine transient electromagnetic method used the emission frequency of 6.25 Hzꎬ 32 superposition and medium density ac ̄ quisition methodsꎬ and it was better to detect the gob and the water-bearing areaꎻ The acquisition optimization parameters were applied to the underground exploration in Jincheng Mining Area. Taking the advanced detection application of Wangtaipu mine as an exampleꎬ it was seen that there was an obvious high-resistance anomaly and two low-resistance anomalies in front of the faceꎬ and the interface was obvi ̄ ousꎬ which was basically the same as the identified distribution area of small coal mine roadway. After drilling and verificationꎬ the high- resistance anomaly area 50 m ahead of the front was emptyꎬ and it was found through the borehole fluoroscopy that it was a small coal mine gob 收稿日期2019-10-18ꎻ责任编辑赵 瑞 基金项目国家科技重大专项资助项目2016ZX05045001-004ꎻ国家自然科学基金青年基金资助项目51804162ꎬ51704162 作者简介牟 义1983ꎬ男ꎬ山东日照人ꎬ副研究员ꎮ Tel010-84264744ꎬE-mail61673275@ qq.com and there was no water. When retreating 30 m backwardꎬ a low-resistance anomaly area was drilled to the right sideꎬ hitting a depth of 35 m 481 牟 义等矿井瞬变电磁法参数优化试验及超前探测应用2020 年第 6 期 with effluentꎬ large amount of waterꎬ and the water was collected in the gobꎬ which has a good detection effect. According to the exploration and summary experience of local coal mine in Ordosꎬ Linfen local coal mineꎬ Jinmei Groupꎬ Shuozhou local coal mine and other underground minesꎬ the apparent resistivity value was generally below 30 Ωmꎬ which was inferred to be water-containing anomalyꎬ and 30 Ωm or more was inferred to be a normal rock formation without waterꎬ and higher than that would be inferred to be a high-resistance anomaly. Key words advance detectionꎻ gobꎻ transient electromagnetic methodꎻ transmitting frequencyꎻ number of superpositionꎻ sampling densityꎻmine 0 引 言 矿井 瞬 变 电 磁 法 Transient Electromagnetic MethodꎬTEM采用电磁波在煤矿井下全空间激发一 次场ꎬ在一次场脉冲的间隔期间断电后 ꎬ观测由 围岩体感应形成的二次场随时间的变化[1]ꎮ 其在 对采煤工作面内部以及巷道掘进前方隐伏小煤窑采 空区、断层、陷落柱、水文钻孔等地质异常体的导含 水性超前探测中取得了良好的应用效果[2]ꎮ 受井 下空间及防爆限制ꎬ矿井瞬变电磁探测回线边长尺 寸和发射功率较小ꎬ为了达到相应的探测深度及探 测精度ꎬ研究合适的数据采集及处理参数尤为重 要[3]ꎮ 国内不少专家学者进行了大量研究ꎬ姜志海 等[4]采用物理试验对矿井瞬变电磁观测系统的发 射功率、发射磁矩、关断时间等参数随发射线圈匝数 的变化关系ꎬ以及关断时间、接收信号随接收线圈匝 数的变化关系进行了定性研究ꎻ李飞等[5]通过数值 模拟和现场试验ꎬ分别研究了矿井半空间和全空间 条件下瞬变电磁法偶极装置观测信号随收发距的变 化规 律ꎻ 陈 明 生 等[6]应 用 时 域 有 限 差 分 方 法 FDTD分析地下和地面的瞬变响应ꎬ采用 35 ms 的延迟时间可探测较浅的异常体ꎻ陈健强[7]开展了 瞬变电磁法探测参数的现场试验研究ꎬ通过发射频 率、叠加次数、数据采集密度等采集参数及视电阻率 系数、深度系数等处理参数的对比试验ꎬ得出了适用 于试验矿区及相似地球物理场条件下的最优探测参 数ꎻ邱浩[8]在朝阳煤矿9101 工作面巷道掘进过程中 通过基于瞬变电磁法频率、叠加次数的参数优化试 验ꎬ确定了工作面内部老旧巷道探测的最优参数ꎻ文 献[9-12]分别从地面瞬变电磁法中心回线、小线 框、浅埋深采空区以及不同系统装置的角度对各类 参数进行了地面试验研究ꎬ优化了相关地面采集技 术ꎬ取得了一定成果ꎮ 在以上研究基础上ꎬ为了获取适合井下施工的 矿井瞬变电磁法参数ꎬ在晋煤集团朝阳矿已知积水 采空区进行了大量参数优化试验工作ꎬ并将优化结 果应用到晋煤集团王台铺矿超前探测采空区及积水 区ꎬ取得一定的成果ꎬ对榆林、鄂尔多斯、山西等地的 矿井瞬变电磁法探测工作可以起到一定的借鉴 作用ꎮ 1 试验区概况及试验设计 1.1 试验区概况 晋煤集团某矿 9101 工作面轨道巷位于 9 号煤 层中ꎬ探测位置东面为矿界ꎬ南部为未采区可能存 在采空区ꎮ 井田内可采煤层为山西组的 3 号煤层 及太原组的 9、15 号煤层ꎮ 2015 年 8 月 23 日在轨道巷 1 152 m 处施工超 前探测钻孔ꎬ在掘进工作面正前方 45 m 处与右帮方 位角 160方向 65 m 处见空ꎬ并有采空区积水涌出 图 1ꎮ 经测定ꎬ水压最大可达 0.2 MPaꎬ涌水量为 1020 m3/ hꎮ 采取水样进行水质化验后ꎬ发现非常 高的矿化度和 SO 2- 4 含量ꎬ具备采空区积水水文地质 特征ꎮ 超前探测钻孔揭露采空区以后ꎬ为了提高瞬变 电磁法探测精度ꎬ在 9101 工作面轨道巷 1 152 m 处 图 1ꎬ进行了矿井瞬变电磁法超前探测参数优化 试验的现场数据采集工作ꎮ 试验现场轨道巷宽度 4.1 mꎬ高度 2.2 mꎬ锚杆锚网支护ꎬ综掘机退后 27 mꎬ掘进工作面正中间 1 号钻孔与右帮 B5 号钻孔作 为放水孔对采空区积水进行疏放ꎬ外接排水钢管ꎬ长 度 1 m 左右ꎬ掘进工作面与右帮管路接口处有少量 涌水ꎮ 距掘进工作面 35 m 处有钢管等铁器ꎬ右帮 距掘进工作面 5 m 处钻孔有渗水ꎬ颜色为黄色ꎮ 图 图 1 物探位置钻孔示意 Fig.1 Schematic of geophysical location drilling 581 2020 年第 6 期煤 炭 科 学 技 术第 48 卷 1 中ꎬ1 号钻孔方位角 90ꎬ倾角+5ꎬ孔深 45 m 时 出水ꎬ出水量约为 10 m3/ hꎬ水压 0.2 MPaꎮ B4 号钻 孔方位角 180ꎬ倾角-1ꎬ孔深 127 mꎬ揭露全煤ꎮ B5 号钻孔方位角 160ꎬ倾角-0.9ꎬ孔深时 65 m 出 水ꎬ出水量约为 10 m3/ hꎬ水压 0.2 MPaꎮ 1.2 试验设计 矿井瞬变电磁法超前探测试验设计探测方向为 顺煤岩层向掘进工作面前方探测ꎬ共布置横向探测角 度 14 个图 2ꎬ探测角度分别是左侧帮180、165、 150、135、120、105、正前方90、 90、右侧帮 0、15、30、45、60、75ꎬ计 14 个测点数据[12]ꎮ 图 2 横向探测角度示意 Fig.2 Schematic of lateral detection angle 2 参数优化试验 2.1 发射频率 矿井瞬变电磁一般要求井下探测区域附近进行 断电处理ꎬ以防止受到生产用电的干扰ꎻ矿井瞬变电 磁仪的发射波形为方形波ꎬ发射频率为普通用电频 率的倍数或约数ꎬ通过滤波可完全对普通用电频率 的干扰进行剔除[13-14]ꎮ 频率越高ꎬ探测深度越小ꎬ 目标体分辨率越高ꎻ频率越低ꎬ探测深度大ꎬ目标体 分辨率越低ꎮ 在频率的选择上ꎬ一般以实际探测时 曲线的完整为标准ꎬ既要兼顾探测深度ꎬ又能兼顾探 测精度[15-16]ꎮ 本次试验选取发射频率为 2.5、6.25 和 25 Hz 时 进行试验叠加次数均为 32 次ꎬ采样密度均为中密 度ꎮ 1 号钻孔方向探测角度 90方向不同发射 频率的归一化电压-时间曲线和视电阻率-深度曲 线分别如图 3 和图 4 所示ꎬ由图 3 可以看出6.25 和 25 Hz 曲线衰减较为圆滑ꎬ2.5 Hz 曲线 40 ms 之前曲 线衰减较为圆滑ꎬ40 ms 之后曲线尾支出现较大的 跳跃性ꎬ曲线衰减不稳定ꎬ且随着频率从2.5、6.25 Hz 到 25 Hz 逐步增大ꎬ其感应二次场随衰减时间的增 加而减少ꎻ在图 4 中表现为随频率增大ꎬ其最大探测 深度减少ꎬ2.5、6.25、25 Hz 最大探测深度分别为 170、120、100 m 左右ꎮ 结合图 1 中在探测深度45 m 打钻出水位置ꎬ可以看出随着频率增大ꎬ视电阻率 也逐步增大ꎬ25 Hz 的曲线视电阻率显示较高ꎬ不能 很好地表现低阻异常ꎬ2.5 Hz 曲线视电阻率虽然较 低ꎬ但衰减不稳定ꎬ影响探测精度ꎬ而 6.25 Hz 曲线 衰减圆滑ꎬ且在出水位置表现出相吻合的低阻异常ꎬ 与现场实际最为符合ꎮ 图 3 1 号钻孔方向测点归一化电压-时间曲线 Fig.3 Normalized voltage-time curves of measuring points in No.1 borehole direction 图 4 1 号钻孔方向测点视电阻率-深度曲线 Fig.4 Resistivity-depth curves of measuring points in No.1 borehole direction 不同发射频率的瞬变电磁法顺层方向超前探测 视电阻率等值线拟断面如图 5 所示ꎬ图中横坐标 0 点位置为探测位置ꎬ横坐标负轴对应巷道左侧ꎬ正轴 对应巷道右侧ꎬ纵坐标为沿顺层方向探测的距离ꎻ图 中从冷色调到暖色调表示视电阻率不断升高ꎮ 从图 5 中可以明显看出2.5 Hz 和 6.25 Hz 的断面图中ꎬ1 号钻孔 45 m 处开始出现凹陷状冷色调低阻异常ꎬ与 现场打钻情况相符ꎬ而 25 Hz 的断面图则出现正常 的视电阻率值ꎬ冷色调低阻异常出现在 70 m 之后ꎬ 与打钻出水情况不符ꎻ考虑到现场掘进工作面左右 帮均有锚网及钻杆等铁器干扰ꎬ在图 5 中所有断面 图中左右帮均出现大片低阻干扰异常ꎬ与 B5 号钻 孔出水位置重合ꎬ不具备验证性ꎮ 2.2 叠加次数 矿井瞬变电磁法在数据采集时ꎬ可能遇到游散 681 牟 义等矿井瞬变电磁法参数优化试验及超前探测应用2020 年第 6 期 图 5 不同发射频率超前探测视电阻率等值线拟断面 Fig.5 Pseudo-section of apparent resistivity isoline for advanced detection with different emission frequencies 电流等随机干扰ꎬ通过多次数据采集并叠加到一起 可以减弱随机干扰ꎬ增强数据的真实性ꎻ但有时遇到 井下交流电缆等周期性干扰ꎬ通过多次数据采集叠 加ꎬ可能增强或突出周期性干扰ꎬ降低探测结果的可 靠性[17-18]ꎮ 另外ꎬ叠加次数越多ꎬ采集时间越长ꎬ施 工效率越低ꎬ所以不能简单增加叠加次数来提高数 据质量ꎬ不能简单从叠加次数的多少来说明探测结 果的质量ꎬ选择叠加次数也是以实际探测时曲线的 完整、圆滑为标准[19-20]ꎮ 本次选取叠加次数为 16、32 和 64 次进行试验 发射频率均为 6.25 Hzꎬ采样密度均为中密度ꎮ 1 号钻孔方向探测角度 90方向不同叠加次数的归 一化电压-时间曲线和视电阻率-深度曲线分别如 图 6 和图 7 所示ꎮ 由图 6 可知3 种叠加次数曲线 衰减均较为圆滑ꎬ只有 16 次叠加曲线尾支出现上 扬ꎬ不能有效压制晚期二次场干扰ꎬ不同之处在于在 图 6 1 号钻孔方向测点归一化电压-时间曲线 Fig.6 Normalized voltage-time curves of measuring points in No.1 borehole direction 图 7 1 号钻孔方向测点视电阻率-深度曲线 Fig.7 Resistivity-depth curves of measuring points in No.1 borehole direction 0.13.0 ms 时ꎬ32 次叠加次数曲线整体归一化电压 值最大ꎬ16 次和 64 次叠加次数曲线最小ꎬ16 次叠加 次数相比较于 64 次整体略大ꎻ图 7 中 32 次叠加次 数曲线整体视电阻率位于 30 Ωm 以下ꎬ而 16 次 和 64 次叠加次数曲线在 30 Ωm 以上ꎬ在 1 号钻 孔 45 m 深度也有类似的表现ꎮ 根据在井下探测总 结的经验ꎬ一般视电阻率值位于 30 Ωm 以下推断 为含水异常ꎬ30 Ωm 以上推断为不含水正常岩层 或高阻异常ꎮ 因此ꎬ32 次叠加次数在 1 号钻孔 45 m 深度表现为低阻含水异常ꎬ与出水位置相吻合ꎬ 与现场实际最为符合ꎬ为矿井瞬变电磁法最佳叠 加次数ꎮ 瞬变电磁法顺层方向不同叠加次数超前探测视 电阻率等值线拟断面如图 8 所示ꎬ可以明显看出32 次的断面图中 1 号钻孔 45 m 处开始出现大片蓝色 图 8 不同叠加次数超前探测视电阻率等值线拟断面 Fig.8 Pseudo-section of apparent resistivity isoline for advanced detection with different number of superpositions 781 2020 年第 6 期煤 炭 科 学 技 术第 48 卷 低阻异常ꎬ与现场打钻情况相符ꎬ而 16 次和 64 次叠 加次数的断面图均出现正常或偏高的视电阻率值ꎬ 与打钻出水情况不符ꎻ左右帮均存在铁器干扰ꎬ断面 图左右帮不具备试验验证性ꎮ 2.3 采样密度 一般采样密度越高ꎬ采集的数据点越多ꎬ数据质 量越可靠ꎬ但采样密度越高ꎬ采集时间越长ꎬ施工效 率就越低ꎬ且受限于计算机处理效率ꎬ不能无限增大 采样密度[21-22]ꎮ 另外ꎬ采样密度越高ꎬ采集的电性 信息越多ꎬ地质信息越丰富ꎬ分辨异常越精细ꎬ同时 采集到的干扰信息也可能增多ꎬ干扰增大ꎬ所以需要 试验选择合理的采样密度来提高数据质量ꎬ不能简 单从增加采样密度来提高异常分辨率ꎬ选择采样密 度也是以实际探测时曲线的衰减圆滑程度为标准ꎬ 更要兼顾实际ꎮ 本次试验选取低密度、中密度和高密度进行试 验发射频率均为6.25 Hzꎬ叠加次数均为32 次ꎮ 1 号钻孔方向探测角度 90方向不同采样密度的归 一化电压-时间曲线和视电阻率-深度曲线如图 9 和图 10 所示ꎬ由图 9 可知ꎬ3 种采样密度曲线衰减 均较为圆滑ꎬ且较为重合ꎬ高密度采样曲线尾支稍显 上翘ꎬ归一化感应电压未见明显差异ꎻ图 10 中表现 为较明显的视电阻率差异ꎬ中密度采样曲线整体视 电阻率最低ꎬ低密度采样曲线视电阻率整体居中ꎬ高 密度采样视电阻率值整体偏高ꎬ但均位于 30 Ωm 以下ꎬ属于低阻含水异常ꎬ低阻异常幅度有差异ꎮ 因 此ꎬ中密度采样曲线吻合程度最高ꎬ为矿井瞬变电磁 法最佳采样密度[23]ꎮ 图 9 1 号钻孔方向测点归一化电压-时间曲线 Fig.9 Normalized voltage-time curves of measuring points in No.1 borehole direction 图 10 1 号钻孔方向测点视电阻率-深度曲线 Fig.10 Resistivity-depth curves of measuring points in No.1 borehole direction 瞬变电磁法顺层方向低密度、中密度、高密度采 样的超前探测视电阻率等值线拟断面如图 11 所示ꎮ 可以看出低密度采样断面图中低阻异常位于正 图 11 不同采样密度的超前探测视电阻率等值线拟断面 Fig.11 Pseudo-section of apparent resistivity isoline for advanced detection with different sampling densities 前方及偏左方向ꎬ中密度低采样断面图中低阻异常 呈现在正前方一片ꎬ左右对称ꎬ而高密度采样断面图 中低阻异常位于正前方及偏右方向ꎬ中密度采样断 面图中的低阻异常可以较好地涵盖现场实际富水区 分布ꎮ 3 应用实例分析 按照以上优化的瞬变电磁法发射频率、叠加次 数以及采样密度等参数ꎬ对山西晋城矿区王台铺矿 进行矿井瞬变电磁法超前探测工作ꎬ以准确查明前 方采空区及积水情况ꎮ 超前探测位置位于 XV3308 工作面的 33083 巷 500 m 处距离开切眼 30 mꎬ探 测成果叠加如图 12 所示ꎬ图中从蓝色到红色冷色 调到暖色调表示视电阻率不断升高ꎮ 可以看出在 掘进工作面两帮随着距离增大ꎬ出现较大面积的低 阻异常区蓝色ꎬ推断为含水采空区ꎻ而在掘进工 作面前方较近区域存在大范围高阻异常区红色ꎬ 推断为不含水或弱含水采空区ꎮ 881 牟 义等矿井瞬变电磁法参数优化试验及超前探测应用2020 年第 6 期 图 12 顺层切片叠加及小窑采空区分布 Fig.12 Stratigraphic slice stacking and small kiln goaf distribution 由图 12 可以看出存在明显的 1 处高阻异常体 和 2 处低阻异常体ꎬ且分界面明显ꎬ其与查明的小窑 巷道分布区基本重合ꎮ 后经过矿方打钻验证ꎬ正前 方 50 m 之外的高阻异常区见空ꎬ通过钻孔透视发 现ꎬ为小煤窑采空区ꎬ无积水ꎬ后退 30 mꎬ向右侧帮 1 处低阻异常区打钻ꎬ打到深度 35 m 遇空出水ꎬ水量 较大ꎮ 可以看出ꎬ瞬变电磁法发射频率、叠加次数以 及采样密度优化结果在晋城矿区煤矿应用取得较好 验证效果ꎮ 4 结论与建议 1矿井瞬变电磁法确定了采用 6.25 Hz 的发射 频率、32 次叠加和中密度采集方式的采空区及积水 区超前探测优化参数ꎬ可取得较好探测效果ꎮ 2根据在鄂尔多斯地方煤矿、临汾地方煤矿、 晋煤集团、朔州地方煤矿等井下探测总结的经验ꎬ一 般视电阻率值位于 30 Ωm 以下推断为含水异常ꎬ 30 Ωm 以上推断为不含水正常岩层ꎬ更高则推断 为高阻异常ꎮ 3将采集优化参数应用到瞬变电磁法实例中ꎬ 成功发现前方采空区及积水区ꎮ 4本次试验由于防爆仪器本身局限性ꎬ部分常 用参数如发射电流等不能试验ꎬ并且尚未考虑资料 处理过程部分参数如视电阻率系数、深度系数等对 数据曲线的影响ꎬ未来试验中会加大这方面参数的 试验工作ꎮ 致谢感谢晋煤集团技术研究院物探工程分公 司以及煤炭科学技术研究院有限公司安全分院晋城 项目部的各位领导、同事对本次试验的重视和支持 参考文献References [1] 申宝宏ꎬ郑行周ꎬ弯效杰.煤矿隐蔽致灾因素普查技术指南 [M].北京煤炭工业出版社ꎬ2015. [2] 李宏杰ꎬ张 彬ꎬ李 文ꎬ等.煤矿采空区灾害综合防治技术与 实践[M].北京煤炭工业出版社ꎬ2016. [3] 闫 述ꎬ陈明生ꎬ傅君眉.瞬变电磁场的直接时域数值分析[J]. 地球物理学报ꎬ2002ꎬ452275-284. YAN Shuꎬ CHEN Mingshengꎬ FU Junmei. Direct time domain nu ̄ merical analysis of transient electromagnetic fields [J]. Chinese Journal of Geophysicsꎬ 2002ꎬ 45 2 275-284. [4] 姜志海ꎬ岳建华ꎬ刘树才.多匝重叠小回线装置形式的矿井瞬 变电磁观测系统[J].煤炭学报ꎬ2007ꎬ3211 1152-1156. JIANG Zhihaiꎬ YUE Jianhuaꎬ LIU Shucai. Mine transient electro ̄ magnetic observation system in the form of multi-turn and small- loop device [J] . Journal of China Coal Societyꎬ 2007ꎬ 32 11 1152-1156. [5] 杨撷民ꎬ霍俊杰ꎬ李彦荣.基于 MIDAS 软件对受采空区影响坡 体的稳定性研究[J].煤炭科学技术ꎬ2019ꎬ47889-95. YANG XieminꎬHUO JunjieꎬLI Yanrong.Study on stability of slope affected by goaf based on MIDAS software[J].Coal Science and Technologyꎬ2019ꎬ47889-95. [6] 李 文.煤矿采空区地面综合物探方法优化研究[J].煤炭科学 技术ꎬ2017ꎬ45 1194-199. LI Wen. Optimization study of surface comprehensive geophysical detection methods of coal mine goafs[J]. Coal Science and Tech ̄ nologyꎬ 2017ꎬ 451 194-199. [7] 姜国庆ꎬ贾春梅ꎬ牟 义ꎬ等.矿井瞬变电磁法全空间处理及三 维解释方法研究[J].地球物理学进展ꎬ2015ꎬ3052164- 2170. JIANG Guoqingꎬ JIA Chunmeiꎬ MU Yiꎬ et al. Study on full-space treatment and 3D interpretation method of mine transient Electro ̄ magnetic Method [J]. Progress in Geophysicsꎬ 2015ꎬ 30 5 2164-2170. [8] 牟 义.典型地质异常体电磁法响应特征研究[J].煤矿开采ꎬ 2017ꎬ2244-9. MU Yi. Study of electromagnetic response characteristics of typical geological anomalous body[J].Coal Mining Technologyꎬ 2017ꎬ 22 4 4-9. [9] 杨云见.不同观测参数条件下中心回线瞬变响应特征分析[J]. 物探装备ꎬ2006ꎬ163216-220. YANG Yunjian. Analyzing transient response of central loop under different conditions of observing parameters[J].Equipment for Ge ̄ ophysical Prospectingꎬ 2006ꎬ 163216-220. [10] 陈明生ꎬ闫 述ꎬ石显新ꎬ等.二维地质体的瞬变电磁场响应特 征[J].地震地质ꎬ 20012252-256. CHEN Mingshengꎬ YAN Shuꎬ SHI Xianxinꎬ et al. Transient electromagnetic response characteristics of two-dimensional geo ̄ logical bodies [J]. Seismology and Geologyꎬ 20012252-256. [11] 郭 纯ꎬ刘白宙ꎬ白登海.地下全空间瞬变电磁技术在煤矿巷 道掘进头的连续跟踪超前探测[J].地震地质ꎬ2006ꎬ 283 456-462. GUO Chunꎬ LIU Baizhouꎬ BAI Denghai. Study on continuous fol ̄ low-up further detection of underground all-time transient elec ̄ tromagnetic technique in coal mining tunnel driving head [J]. Seismology and Geologyꎬ 2006ꎬ 28 3 456-462. 981 2020 年第 6 期煤 炭 科 学 技 术第 48 卷 [12] 刘志新ꎬ岳建华ꎬ刘仰光.扇形探测技术在超前探测中的应用 研究[J].中国矿业大学学报ꎬ2007ꎬ366822-868. LIU Zhixinꎬ YUE Jianhuaꎬ LIU Yangguang. Application of fan- shaped detection technology in advanced detection [J]. Journal of China University of Mining and Technologyꎬ 2007ꎬ 36 6 822-868. [13] 程久龙ꎬ王玉和ꎬ于师建ꎬ等.巷道掘进中电阻率法超前探测原 理与应用[J].煤田地质与勘探ꎬ2000ꎬ28460-62. CHENG Jiulongꎬ WANG Yuheꎬ YU Shijianꎬ et al.Development of resistivity method for advance exploration in roadway excavation [J]. Coalfield Geology and Explorationꎬ 2000ꎬ 28 4 60-62. [14] 牟 义.切片技术在煤矿回采工作面瞬变电磁法探测中的研 究与应用[J].中国煤炭ꎬ2014ꎬ40335-39. MU Yi. Research and application of slicing technology in transient electromagnetic detection in coal mining face[J]. China Coalꎬ 2014ꎬ 403 35-39. [15] 陈健强ꎬ张俊英ꎬ李宏杰ꎬ等.掘进面瞬变电磁探测采空区及三 维成像技术研究[J].煤矿开采ꎬ2017ꎬ22118-21. CHEN Jianqiangꎬ ZHANG Junyingꎬ LI Hongjieꎬ et al. 3-D ima ̄ ging technology and goaf transient electromagnetic exploration of driving face[J]. Coal Mining Technologyꎬ 2017ꎬ 221 18- 21. [16] 窦文武ꎬ牟 义ꎬ邱 浩ꎬ等.基于 Voxler 的瞬变电磁勘探三维 解释技术研究[J].煤炭科技ꎬ2015437-39. DOU Wenwuꎬ MU Yiꎬ QIU Haoꎬet al. Study on 3D interpretation technology of transient electromagnetic exploration based on voxler [J]. Coal Science and Technologyꎬ 20154 37-39. [17] 夏春鹏ꎬ窦文武ꎬ樊林林ꎬ等.不同视电阻率导体对瞬变电磁法 的响应特征分析[J].中州煤炭ꎬ201511114-117. XIA Chunpengꎬ DOU Wenwuꎬ FAN Linlinꎬ et al.Characteristics of response of different apparent resistivity conductors to transient electromagnetic method[J]. Zhongzhou