频率域电磁法在煤矿采空区调查中的应用
姜国庆1,2, 贾春梅3, 谭强4, 尚通晓1,2, 朱首峰1,2
1. 国土资源部地裂缝地质灾害重点实验室,江苏 南京 210018
2. 江苏省地质调查研究院 基础地质研究所,江苏 南京 210018
3. 中国石油化工股份有限公司石油物探技术研究院,江苏 南京 211103
4. 西北核技术研究所,陕西 西安 710024

作者简介: 姜国庆(1986-),男,物探工程师,硕士研究生,2011年毕业于山东科技大学,主要从事电磁法勘探研究及应用工作。E-mail:jiangguoqing220@126.com

摘要

采取有效技术手段进行煤矿采空区调查已成为矿区生态环境治理、矿井安全生产和重大工程建设项目选址中亟需解决的问题。本文以频率域电磁法基本原理为基础,对比分析EH-4电导率成像系统和可控源音频大地电磁法(CSAMT)的方法特点,通过数值模拟探讨不同介质充填采空区的电性响应特征,并结合典型工程实例对EH-4电导率成像系统和可控源音频大地电磁法探测煤矿采空区的实际应用效果进行分析。数值模拟结果和工程实践表明,煤矿采空区相对围岩具有显著的电性差异,以此为物性基础采用频率域电磁法进行采空区探测可以取得较为理想的效果。频率域电磁法在探测煤层采空区方面是一种行之有效的方法,可以成为煤矿防治水和采空区地质灾害调查的有效技术方法。

关键词: 煤矿采空区; 地球物理特征; 频率域电磁法; 探测效果
中图分类号:P631 文献标志码:A 文章编号:1000-8918(2015)03-0646-05
Application of fdem in detecting coal mined-out areas
JIANG Guo-Qing1,2, JIA Chun-Mei3, TAN Qiang4, SHANG Tong-Xiao1,2, ZHU Shou-Feng1,2
1. Key Laboratory of Earth Fissures Geological Disaster, Ministry of Land and Resources, Nanjing 210018, China
2. Basic Geological Research Institute, Geological Survey of Jiangsu Province, Nanjing 210018, China
3. Research Institute of Geophysical Prospecting, SinoPEC, Nanjing 211103, China
4. Northwest Institute of Nuclear Technology, Xi’an 710024, China;
Abstract

Detecting coal mined-out areas with effective measures has been major issue to be addressed in mining area environment management, mine safe production and major construction projects site selection. Based on the fundamental principle of frequency domain electromagnetic method (FDEM), method characteristics of EH-4 and CSAMT was compared and analyzed, electrical characteristic of mined-out area filled with different materials was analyzed through numerical modeling, and practical detecting effect was discussed by typical project cases of EH-4 and CSAMT. Numerical modeling result and engineering practice prove that compared with surrounding rock mined-out area has prominent electrical difference, and base on this physical property mined-out area detecting with FDEM can obtain ideal effect. FDEM is an effective method in detecting coal mined-out area and it has already been an effective technological means in coal mine water prevention and control and geological disaster survey of mined-out area.

Keyword: coal mined-out area; geophysical characteristic; FDEM; detecting effect

煤矿开采活动破坏了地下岩体的原始应力平衡, 引起上覆岩层冒落、断裂、弯曲变形等, 从而形成形态、规模各异的采空区。采空区的存在破坏了矿区的生态环境并对矿井安全生产和人民生命财产安全造成极大威胁, 同时对地面建筑物、构筑物的结构稳定性造成破坏。因此, 采取有效技术手段查明采空区的埋藏深度、大小、形态等特征, 并对采空区附近地质稳定性进行评价已成为矿区生态环境治理、矿井安全生产和重大工程建设项目选址中亟需解决的问题。

频率域电磁法中应用最早的是大地电磁法(MT), 其探测深度可达数十千米到数百千米[1], 主要用于研究深部地质构造问题[2, 3, 4]。后来, 随着能检测微弱音频大地电磁场信号的高灵敏度磁传感器和接收设备的出现, 音频大地电磁法(AMT)得以实现[5]。MT法和AMT法共同的特点是采用天然场源, 由于天然场源具有信号弱的特点, 其缺点为信噪比低、抗干扰能力差。为了弥补天然信号的不足, 人们将频率和功率都可控的人工场源引入到频率域电磁法勘探中, 出现了人工场源的EH-4电导率成像系统和可控源音频大地电磁法(CSAMT)。由于煤矿采空区调查时, 施工区域往往接近工矿区和居民区, 电磁干扰较为严重, 因此, EH-4电导率成像系统和CSAMT法在煤矿采空区调查中具有一定优势。

1 两种频率域电磁法对比

EH-4电导率成像系统是美国EMI和Geometrics公司联合生产的部分可控源的大地电磁测深仪[6], 该系统与CSAMT法[7]均属于引入人工场源的频率域电磁方法, 因此, 都具有抗干扰能力强、信噪比高的特点。但是由于EH-4系统仅在高频段采用人工场源, 而CSAMT法由低频到高频全部采用人工场源(例如美国GDP-32系统标准频率范围为1/64 Hz~8 192 Hz), 因此, CSAMT的探测深度大于EH-4系统, 但是EH-4系统的浅部分辨率要优于CSAMT法。同时, CSAMT法一般采用赤道偶极装置进行标量测量, 一次供电可以完成7个测深点的测量工作, 观测效率较高, 而EH-4电导率成像系统一般进行张量测量, 工作效率降低, 但可获得地下地质体更为丰富的地电信息。

EH-4系统和CSAMT法均受到场源效应的影响, 为了有效减小场源效应, 一方面通过选择合理的收发距, 另一方面需要在数据处理中采用非远区校正算法进行校正。由于CSAMT法人工供电的低频信号频率较低, 因此要求比较大的收发距(通常5~10 km) , 相应的发射装置的功率要求也比较高, 发射装置比较笨重, 而EH-4系统则减少了低频发射的笨重设备, 收发距缩短(400~500 m), 便于野外工作。但是当工作区面积较大或者测线较长时, EH-4系统需要频繁的移动发射站点。

在测线布置上, CSAMT方法测线布置受到发射源的限制, 测点分布平行于供电电极AB的方向, 而EH-4系统由于是连续单点观测, 测线既可以沿直线布置, 也可以根据地形或地质任务要求布置成弯曲测线, 因此, EH-4系统可以非常灵活地适用于各种不利的地形条件, 野外适应性更强。

因此, 针对具体的探测任务, 应根据施工条件、地质任务要求, 在兼顾探测效果和工作效率的前提下合理选择探测方法。

2 采空区电性特征数值模拟
2.1 煤矿采空区顶板“ 三带”

地下煤层被采出后, 在地下形成了一个空洞, 使原本平衡的地下应力系统遭到破坏, 岩层在上覆岩体的应力作用下发生变形、断裂、位移和冒落。顶板覆岩在垂直方向上的破坏可分为冒落带、裂隙带、弯曲变形带, 简称顶板“ 三带” , 如图1所示。

图1 煤矿采空区顶板“ 三带” 示意

采空区顶板“ 三带” 中冒落带和裂隙带高度的研究方法主要有经验公式法、物理模拟、数值模拟和现场实测等。研究成果表明, 顶板覆岩冒落带和裂隙带的发育高度是煤层厚度的几倍甚至几十倍, 冒落带和裂隙带的存在, 改变了上覆围岩的物理性质, 相对增大了采深比, 扩大了采空区的异常范围, 因此, 利用物探方法探测采空区, 并不仅仅是探测与采出煤层等体积的空洞, 而是探测包括空洞、冒落带和裂隙带在内的整个区域。

煤层开采以后, 采空区附近短期内形成一定规模的充气空间, 其地层电阻率与围岩具有显著差异[8], 但是由于此时充气空间规模有限, 当采深比较大时利用物探方法很难有效探测。开采结束一段时间后, 当采空区上方岩层在重力作用下发生坍塌变形, 冒落带、裂隙带发育成熟时, 若地下水通过裂隙向采空区汇集使采空区充水, 则由于水中溶解有大量的电解质, 采空区地层表现为低阻高极化的电性特征, 相反, 如果采空区未充水, 则由于地层中空洞、离层和裂隙的存在使岩矿石导电性变差, 此时, 采空区地层相对于围岩表现为高阻特征。

2.2 数值模拟分析

为了直观地讨论采空区的电性特征, 分别建立充水采空区(低阻)和未充水采空区(高阻)模型进行二维频率域电磁法有限元数值模拟。

如图2a为充水的低阻采空区模型, 第四系厚度20 m, 电阻率20 Ω · m; 煤系地层厚400 m, 电阻率100 Ω · m; 煤层厚5 m, 电阻率500 Ω · m; 采空区顶板覆岩“ 三带” ( 冒落带、裂隙带、弯曲变形带)厚度分别为20 m、60 m、20 m, 其电阻率分别为20 Ω · m、50 Ω · m、80 Ω · m, 导水裂隙带形状为两头高、中间低的马鞍状。图2b为采煤工作面斜长L变化时的电阻率断面, 由断面图可以看出, 受采空区充水、地层中岩矿石导电性增强的影响电阻率等值线呈现明显的下凹特征, 并且随着开采斜长L的增加下凹特征更加明显。

图2 低阻采空区模型数值模拟结果

图3 高阻采空区模型数值模拟结果

图3a为采空区顶板“ 三带” 充分发育但未充水的高阻采空区模型, 模型中顶板三带发育厚度与上述低阻采空区模型相同, 只是顶板“ 三带” 电阻率相对于围岩为高阻, 冒落带、裂隙带和弯曲变形带电阻率分别为1 000 Ω · m、800 Ω · m和150 Ω · m。图3b为采煤工作面斜长L变化时的电阻率断面, 由于

高阻采空区的存在, 电阻率等值线出现明显上凸形状, 并且随着采煤工作面斜长L的增加上凸特征更加明显。

数值模拟结果说明, 采空区的存在引起覆岩电阻率的显著变化, 当煤层采空区和顶板“ 三带” 的电阻率值与围岩电阻率具有显著差异, 并且采空区具有一定的规模时, 应用频率域电磁法进行煤层采空区探测在理论上是完全可行的。

3 工程实例
3.1 EH-4在山东某地块探测效果

山东某棚户区改造项目拟建场区毗邻一废弃的小煤窑, 为了对区域上方建筑的安全性进行评价, 必须首先查明小煤窑采空区的分布范围。本区含煤地层含上、中、下三层主要可采煤层, 三层煤总平均厚度为5.50 m, 煤层顶板埋深-100 m左右。根据已知地质资料和地质任务要求, 并考虑到地表实际情况, 测区范围内设计测线7条, 由北向南测线编号L1~L7。测线点距10 m, 控制目的层深度为± 0 m~-250 m。仪器设备采用美国EH-4电导率成像系统, 数据处理采用IMAGEM和EMAGE-2D数据处理系统, 正反演相结合, 反演采用BOSTIK反演。

通过EH-4电导率法探测结果可知, 在测区北部的L1~L5五条测线的反演电阻率断面上均有不同程度的采空区电性特征显示, 而在测区南部的L6和L7两条测线的反演电阻率断面上电性变化较为均匀, 无采空区特征显示。将L1~L5五条剖面上解译的采空区位置标定在平面图上可以看出测区范围内地下存在的采空区为一连续的采空区域, 本文选择两条具有代表性的剖面进行对比分析。

图4 山东某地EH-4探测反演电阻率断面

图4a为北部L1线反演电阻率断面。断面图清晰地反映出地下煤层采空区及上覆岩层移动、破坏的电阻率异常特征。测线方向30~110 m桩号段, 埋深70~100 m范围内, 地层电阻率等值线呈现下凹形状, 表明该段煤层开采后冒落带、裂隙带充水, 致使地层电阻率明显降低, 这是典型的采空区异常电性特征, 推断该区段煤层已被开采, 地层稳定性较差。

图4b为南部L6线反演电阻率断面, 可以看出测线下方地层电阻率具有很好的连续性, 未发现明显的煤层采空区低阻或高阻电阻率异常, 推断该测线地下无采空区存在, 该段地层稳定性好。

3.2 CSAMT在江苏某矿区探测效果

矿区位于灵山— 龙王山次级倾伏背斜, 开采煤层为二叠系龙潭组上1煤层。矿区范围内基岩裸露, 地表出露以三叠系青龙组灰岩为主。如图5为矿井采掘平面图, 从图中可以看出该废弃矿井为一对双斜井开拓, 在井筒的东西两侧分别为东、西两个采区, 图中阴影部分为已知采空区的范围。现场实地踏勘依然可以看到遗留的双斜井井口, 井筒内充水。为了查明煤矿后期开采时是否已向北部即深部煤层延伸, 本次探测在已知采空区的北部布置4条CSAMT法测线, 测线自南东向北西编号2线、4线、6线和8线, 测线点距20 m。仪器设备采用美国GDP-32多功能电法仪, CSAMT赤道偶极标量法测量, 频率范围8~8 192 Hz。野外采集的原始数据首先进行干扰剔除、去噪、近场校正及静态效应校正等预处理, 然后进行反演成像得到各测线反演电阻率断面。

图5 江苏某矿区CSAMT探测测线布置

图6 江苏某矿区CSAMT探测2线反演电阻率断面

如图6为2线反演电阻率断面, 由图中可以看出由浅到深地层电阻率表现为“ 高— 低— 高” 的变化特征。由已有地质资料可知, 2线断面中龙潭组上1煤层顶板埋深为-550 m左右。由断面图可以看出1 040 m ~1 260 m桩号段, 上1煤层位置附近地层电阻率等值线呈现明显的“ V” 字型低阻异常, 表明该段煤层开采后上覆岩层发生冒落、弯曲变形破坏, 改变了原来地层的连续性, 产生了大量的离层、空洞、裂缝且具有一定的含水性, 致使地层电阻率明显变化, 推断该区段煤层已被开采。而测区北部的4线、6线和8线电阻率断面图上无明显采空区电性特征显示。因此, 推断该废弃矿井西采区在开采过程中已经向深部延伸, 已知采空区北部的煤层也已经被开采。该探测结果已经被后期收集到的接近闭坑时的矿井资料所验证, 探测结果与实际开采情况完全吻合, 这充分说明了可控源音频大地电磁法应用于煤矿采空区探测的可行性和有效性。

4 结论

煤层开采后, 采空区附近应力场的变化将引起覆岩电阻率的显著变化, 该变化可以作为频率域电磁法探测采空区的地球物理前提。

数值模拟和工程实践表明, 频率域电磁法在探测煤层采空区方面是一种行之有效的方法, 该方法具有勘探精度高、效率高、成本低等优势, 已成为煤矿防治水和采空区地质灾害调查的有效技术手段。

The authors have declared that no competing interests exist.

参考文献
[1] 李金铭. 地电场与电法勘探[M]. 北京: 地质出版社, 2005. [本文引用:1]
[2] 张振宇, 王绪本, 方慧. 龙门山构造带中段大地电磁测深研究[J]. 物探与化探, 2012, 36(3): 377-381. [本文引用:1]
[3] 李爱勇, 柳建新, 朱春生, . 大地电磁测深在桂中坳陷油气勘探中的应用[J]. 物探与化探, 2012, 36(1): 8-12. [本文引用:1]
[4] 高景宏, 佟铁钢, 强建科, . 开封凹陷区地热资源大地电磁测深研究[J]. 物探与化探, 2010, 34(4): 440-443. [本文引用:1]
[5] 陈文华. 音频大地电磁法及其今后的研究方向[J]. 石油地球物理勘探, 1981, (3): 62-74. [本文引用:1]
[6] 郭建强, 武毅, 邵汝君, . StratagemTMEH-4电导率成像系统简介及应用[J]. 物探与化探, 1998, 22(6): 458-464. [本文引用:1]
[7] 何继善. 可控源音频大地电磁法[M]. 长沙: 中南工业大学出版社, 1990. [本文引用:1]
[8] 程久龙, 于师建. 覆岩变形破坏电阻率响应特征的模拟实验研究[J]. 地球物理学报, 2000, 43(5): 699-706. [本文引用:1]
[9] 刘菁华, 王祝文, 朱士, . 煤矿采空区及塌陷区的地球物理勘查[J]. 煤炭学报, 2005, 30(6): 715-719. [本文引用:1]
[10] 程建远, 孙洪星, 赵庆彪, . 老窑采空区的探测技术与实例研究[J]. 煤炭学报, 2008, 33(3): 251-255. [本文引用:1]
[11] 杨建军, 吴汉宁, 冯兵, . 煤矿采空区探测效果研究[J]. 煤田地质与勘探, 2006, 34(1): 67-70. [本文引用:1]
[12] 李娟娟, 潘冬明, 胡明顺, . 煤矿采空区探测的几种工程物探方法的应用[J]. 工程地球物理学报, 2009, 6(6): 728-732. [本文引用:1]
[13] 张莹, 张胜业. EH-4资料处理解释系统的研究[J]. 工程地球物理学报, 2005, 2(4): 311-315. [本文引用:1]
[14] 王玉和, 王厚臣, 程久龙, . EH-4电磁方法在断层定位及导水性探测中的应用[J]. 山东科技大学学报: 自然科学版, 2010, 29(4): 35-39. [本文引用:1]
[15] 吴璐苹, 石昆法, 李荫槐, . 可控源音频大地电磁法在地下水勘查中的应用研究[J]. 地球物理学报, 1996, 39(5): 712-717. [本文引用:1]
[16] 杨长福, 徐世浙. 国外大地电磁研究现状[J]. 物探与化探, 2005, 29(3): 243-247. [本文引用:1]
[17] 穆海杰, 王红兵. CSAMT法在南水北调中线采空区探测中的应用[J]. 工程地球物理学报, 2008, 5(3): 321-325. [本文引用:1]
[18] 刘红涛, 杨秀瑛, 于昌明, . 用VLF、EH-4和CSAMT方法寻找隐伏矿—以赤峰柴胡栏子金矿床为例[J]. 地球物理学进展, 2004, 19(2): 276-285. [本文引用:1]