基于地震波和电磁波CT联合探测的采掘巷道冲击危险性评价方法

基于地震波和电磁波CT联合探测的采掘巷道冲击危险性评价
方法
LIU Shaohong;PAN Junfeng;WANG Hongtao;QI Qingxin;TANG Zhongyi;FENG Meihua;LOU Jinfu;ZHANG Yulong
【摘要】电磁波CT能够探测掘进工作面,但无法探测掘进工作面后方的大面积区域;地震波CT能够探测回采工作面内部,但是探测不了两侧的煤柱,即基于单一CT 探测手段的冲击危险性评价方法,在矿井采掘巷道存在评价盲区.针对上述实际问题,建立了一种基于地震波和电磁波CT联合探测的采掘巷道冲击危险性评价方法;将2种现场实测手段电磁波CT和地震波CT2相结合,研究了基于2种CT探测的冲击危险性评价方法的数学表达;以时间和地点为状态变量,以现场动力显现情况及CT 探测结果为驱动变量,确定了评价方法的内在逻辑关系及实施流程.通过分析基于地震波和电磁波CT探测的冲击危险性评价的力学基础,建立了冲击危险性评价模型;在2种CT探测的基础上,实现了采掘巷道冲击危险等级的确定及危险区域的划分;将该评价方法进行工程实践表明,通过将2种CT探测方式联合起来,实现了掘进工作面及其后方大面积区域,回采工作面及区段煤柱区的冲击危险性评价,验证了评价方法的有效性.该方法克服了基于单一CT探测手段的评价方法在探测尺度上的局限,几乎覆盖了矿井采掘巷道所有潜在冲击危险区域,为冲击地压治理提供了较为准确的依据,具有较好的应用效果.在时间上,该方法将CT探测与矿井的采掘接续相结合,初步确定探测顺序,并根据实际的动力显现情况和之前CT探测的结果加以调整.从空间来说,覆盖了掘进面、回采工作面和煤柱等矿井采掘巷道几乎所有潜在冲击危险区域,基本保证了无评价盲区.并且,该方法充分发挥了震动波与电磁波CT探测操作过程的便捷性,能够较好的适应井下现场动力显现的突发性.
【期刊名称】《煤炭学报》
【年(卷),期】2018(043)011
【总页数】12页(P2980-2991)
【关键词】地震波CT;电磁波CT;冲击危险性评价;采掘巷道;煤柱
【作者】LIU Shaohong;PAN Junfeng;WANG Hongtao;QI Qingxin;TANG Zhongyi;FENG Meihua;LOU Jinfu;ZHANG Yulong
【作者单位】;;;;;;;
【正文语种】中文
【中图分类】TD324
根据我国煤矿冲击地压案例，掘进工作面冲击地压有时发生在掘进工作面，有时发生在掘进工作面后方100 m以外区域，例如新疆硫磺沟煤矿“1·30”冲击地压事
故和义马千秋煤矿“11·3”冲击地压事故;当煤柱宽度留设不合理时，回采工作面
邻空巷道的煤柱帮常常出现冲击地压显现，例如鹤岗峻德矿“3·15”冲击地压事故。

因此，对采掘巷道的冲击危险性进行准确评价，进而指导冲击地压治理工作，具有重要的现实意义。

冲击危险性评价的传统方法包括综合指数法、动态权重法、冲击地压可能性指数诊断法等[1-3];近年来，基于现场监测及探测手段的冲击危险性评估方法被广泛应用[4-8];这些评价方法不同程度的指导了现场冲击地压防治工作。

但是，由于煤岩地
质条件的多变性、复杂性，基于CT探测的评价方法存在一定局限。

具体来说，电磁波CT能够探测掘进工作面，但无法探测掘进工作面后方的大面积区域;地震波
CT能够探测回采工作面内部，但是探测不了两侧的煤柱。

换而言之，采用单一的CT探测手段进行危险性评价时，无法覆盖矿井所有采掘巷道的围岩，存在评价盲区。

为了解决上述问题，将电磁波CT与地震波CT探测手段联合起来，建立矿井采掘巷道冲击危险性评价方法。

首先，研究基于地震波和电磁波CT探测的冲击危险性评价方法的数学表达，获得该方法的内在逻辑关系及流程;然后，通过分析基于地震波和电磁波CT探测的冲击危险性评价的力学基础，建立冲击危险性评价方法;将该方法在平庄矿区和彬长矿区进行工程实践，表明该方法几乎能够覆盖矿井采掘巷道的潜在冲击危险区域。

1 基于地震波和电磁波CT探测的冲击危险性评价方法的数学表达
为了明确基于地震波和电磁波CT探测的冲击危险性评价方法的内在逻辑关系及流程，对该方法的数学表达进行研究。

首先，对研究问题进行分析，问题中涉及的变量包括，探测时期:掘进卸压前T1、掘进卸压后T2、回采卸压前T3和回采卸压后T4;探测区域类型:掘进工作面区域X1，掘进工作面后方X2，煤柱区X3，回采工作面X4;探测方式:地震波CT探测Z1，电磁波CT探测Z2;现场实际动力显现情况:巷道围岩破坏性冲击显现为F1，强震动、瞬间围岩变形、煤体弹射等动力现象为F2，无冲击显现为F3。

问题分析中，需要进行一系列假设。

假设1:煤矿具备2种探测设备各1套，且均可正常使用;假设2:矿上现有的条件可以开展2种CT探测工作;假设3:根据2种CT探测手段的特点，探测方法与时期、空间是一一对应的关系，具体见式(1)，表1。

Zm =Z(Ti,Xj)(m=1,2;i=1,2,3,4;j=1,2,3,4)
(1)
表1 不同时期和地点的CT探测手段Table 1 CT detection methods at
different time and locations区域掘进卸压前T1掘进卸压后T2回采卸压前T3回采卸压后T4掘进工作面区域X1Z2Z2——掘进工作面后方X2Z1Z1——煤柱区X3Z2Z2Z2Z2回采工作面X4——Z1Z1
问题的目标是根据矿井的实际采掘布局接续及现场动力显现情况，得到不同时期、不同地点、由不同探测方式得到的危险性Wn。

探测流程为给定二维空间R2上的函数，其自变量为探测区域内的空间坐标xj，探测方式Zm，即
Wn{Fo→Fp,Wn-1}=W(xj,Zm)
(j=1,2,3,4;m=1,2;o,p=1,2,3;o≤p)
(2)
其中，n由CT探测中冲击危险等级的划分标准来确定。

在上述问题中，现场动力显现情况F以及上个区域的CT探测结果属于驱动变量，其中现场动力显现比较强的区域优先探测，并且上个区域的CT探测出的危险区域也为下个区域的CT探测提供依据，由此建立基于地震波和电磁波CT联合探测的冲击危险性评价的数学表达式(2)。

2种CT探测方式正是由于一种探测方式的结果能够对另一种产生影响而联合起来，具体表述详见第5节讨论部分。

在数学表达中，探测区域的空间坐标xj和探测方式Zm属于状态变量，探测方式Z由探测时期T和探测区域类型X确定，见式(1)，表1。

需要注意的是，每一个矿井的动力显现情况不同，动力显现情况F1，F2，F3不一定全部具备，同时满足动力显现情况F1，F2，F3的时间和地点也可能不惟一，由于F仅是驱动变量不影响整个函数(2)的一一对应性。

此外，2个驱动变量的优先级是F>Wn-1。

最终的评价结果是根据CT探测结果对探测区域内每个空间坐标都赋予一个冲击危险性。

不失一般性，举例说明了基于CT探测的采掘巷道冲击危险性评价方法的结果。

该方法的评价结果包括矿井采掘巷道的评价流程以及冲击危险性，如图1所示。

流程图中体现了以下逻辑关系:首先，矿井各采掘区域的动压显现程度作为“驱动变
量1”，确定了联合探测流程的先后顺序;其次，已完成的探测结果作为“驱动变量2”，能够对矿井未来的探测方式产生影响;最后，当上述两者确定后，再根据时间和地点确定探测方式，进而得到冲击危险等级及危险区域。

图1 基于地震波和电磁波CT探测的冲击危险性评价结果Fig.1 Evaluation of impact hazard based on seismic wave and electromagnetic wave CT detection
2 基于地震波和电磁波CT探测的冲击危险性评价的力学基础
根据已有研究成果[9]，由煤岩变形破坏过程应力-应变曲线，建立的煤岩体冲击危险性判别准则为
(3)
式中，Wε(t)为煤岩体破坏危险性指数，无量纲;t为时间，s;ε(t)为煤岩体当前时刻的应变值;ε0为煤岩体出现微破裂时的初始应变值;ε1为煤岩体完全破坏时的最终应变值。

2.1 电磁波CT的力学基础
电磁波CT法是利用无线电波在2个钻孔中分别发射和接收，根据不同位置上接收的场强大小，来确定地下不同介质分布的一种地下地球物理勘查方法，它的工作频率最高可达32 MHz。

下式为电磁波CT法中的场强观测值公式[10-11]:
(4)
(5)
其中，E为接收点的场强值，为初始辐射常数，无量纲;β为吸收系数，即介质中单位距离对电磁波的吸收值，dB/m;f(θ)为收发天线的方向因子函数;r为发射与接收
点之间的距离，m;ω为天线圆频率，为介电常数，F/m;μ为磁导率，H/m;Cσ为
电导率，1/(Ω·m)，σ=1/z，z为电阻率，Ω·m。

由式(5)，电磁波CT探测中煤的电磁波吸收系数主要随介电常数和电阻率的变化而变化。

由实验室实验结果[17]可知，煤的介电常数和电阻率与煤体承受应力相关，而煤体承受应力与冲击危险性关系密切，因此煤岩内电磁波CT探测与冲击危险性是相关的。

通过实验发现[12]，煤的电阻率z、介电常数ε与其承受应力σ间呈线性关系，如式(6)，由此可得电磁波吸收系数与应力近似呈线性关系，即与冲击危险性呈线性
关系。

(6)
式中，a0,b0,c0,d0均为常数，与煤体自身的性质有关。

2.2 地震波CT的力学基础
在应力传播路径上，根据动量守恒定律，沿传播方向，应力、振速和波速之间的相互关系为
(7)
式中，σm为P波作用产生的正应力，Pa;τ为S波作用产生的剪应力，Pa;ρ为传
播介质的密度，kg/m3;vP和vS为P，S波导致的介质振速，m/s;VP和VS为P，S波波速，m/s。

综上分析，电磁波吸收系数和地震波波速与冲击危险性间呈近似线性关系，因此电磁波与震动波探测的评价指标均可参照式(3)进行归一化处理[9]，并需根据参量的
物理意义加以调整。

3 基于地震波和电磁波CT探测的冲击危险性评价模型
根据地震波和电磁波CT的力学基础，分别建立基于2种CT探测手段的冲击危险性评价模型。

3.1 基于地震波CT探测的冲击危险性评价模型
3.1.1 波速异常系数
(8)
式中，为测区内围岩纵波波速平均值，为测区内围岩极限纵波波速值(岩层临界破坏时)，m/ms。

从震波波速角度表征围岩冲击危险性，其值可正，亦可负。

若AC为正，则表明该处围岩大概率处于应力集中状态，该值越大应力集中程度越高;若AC为负，应力集中程度则较低。

3.1.2 波速梯度系数
波速梯度系数GC从震波波速梯度角度表征围岩冲击危险性，其表达式为
(9)
式中，GP为测区内围岩某点的纵波波速梯度，为测区内围岩极限纵波波速梯度，即围岩处于临界破坏状态，s-1。

速度场梯度GP为波速变化率最大方向上的梯度值。

一般对离散数据周围8个节点求取一阶方向导数，取其最大值，可表示为
(10)
式中，d为网格边长，m;x，y分别为相邻网格的纵、横方向上的编号。

3.1.3 评价模型
综合考虑波速大小和波速梯度的影响，建立冲击地压危险性评价模型为
(11)
其中，C为地震波CT的冲击地压危险性指数;a1，b1为权重系数，均取0.5。

VP 可通过层析成像技术反演获取;GP则通过式(9)计算得出;可通过地震波CT反演结果估算:对于无明显动力现象煤岩层，取1.2 maxVP，1.2 maxGP，对于弱动力显现区域取1.1 maxVP，1.1 maxGP，对于强烈动力显现区域分别取maxVP，maxGP，式中maxVP为反演出的波速CT图像中最大波速值，maxGP为依据波速CT图像计算得出的最大波速梯度[8]。

3.2 基于电磁波CT探测的冲击危险性评价模型
3.2.1 吸收系数异常指数
考虑电磁波吸收系数对煤体冲击危险性的影响，构建吸收系数异常指数如下:
(12)
式中，β为吸收系数的实测值，dB/m;βmax为探测区域围岩对电磁波吸收系数的最大值，dB/m;β0为探测区域围岩对电磁波吸收系数的平均值，dB/m;α为动压特征参数。

根据以往实测案例及实验结果，α在无动压显现区域取1.2，强动压显现区域取1.0，弱动压显现区域取1.1[8]。

吸收系数异常变化与煤岩体裂隙发育程度之间的关系为:正值表征电磁波衰减异常与围岩破坏程度的关系，数值越大说明破坏程度越高，此时对于支护内围岩该值越大说明支护质量越差，而对于支护外围岩则为应力释放越充分;负值指数表征电磁波衰减异常与围岩完整程度的关系，数值越小则完整性越好,支护内围岩该值越小说明支护质量越好，而支护外围岩该值越小则应力越集中。

3.2.2 吸收系数梯度指数
考虑电磁波吸收系数变化梯度与冲击危险性的关系，构建吸收系数梯度指数的计算
式如下:
(13)
式中，Gβ为测区内围岩某节点电磁波吸收系数的梯度，dB/m2;Gβmax为现场条件下测区内电磁波吸收系数梯度的最大值，dB/m2;α为动压特征参数，取值同前。

吸收系数梯度Gβ取测区内某一节点相邻电磁波吸收系数变化率的最大值。

一般对离散数据周围8个节点求取一阶方向导数，最后取其最大值，中心网格的吸收系
数梯度可表示为
(14)
式中，d为网格边长，m;x，y分别为相邻每一个网格的纵、横方向编号。

3.2.3 评价模型
综上，近场围岩的冲击危险指数记作:
(15)
式中，Ds为近场围岩的冲击危险指数;h为探测深度，m;hs为支护结构的范围，m;本模型中k，l均取0.5[9,12-13]。

3.3 冲击危险等级的划分
根据统计学、模糊数学、实验室试验和现场实测等，模型中将冲击危险定量化分为4个等级，包括无冲击危险、弱冲击危险、中等冲击危险和强冲击危险，见表
2[8,14-15]。

危险性指数最大值为1，由于波速和吸收系数异常指数存在负值，相应的危险性指数同样存在负值[13]。

表2 冲击危险等级划分Table 2 Classification of danger grade围岩指标冲击危险等级C和Ds值① 无< 0.25② 弱(0.25,0.50)③ 中(0.50,0.75)④ 强(0.75,1.00)
4 工程实践
4.1 掘进工作面的工程实践
4.1.1 掘进工作面概况及探测方案
探测地点为内蒙古某矿西068掘进工作面，埋深平均541 m，煤层平均厚度17.2 m，倾角约30°，煤层顶板为抗压强度147 MPa的岩浆岩，底板岩性为砂岩;西068工作面长度99 m，阶段煤柱净宽8 m，巷道宽为4.5 m。

2015-08-
05T17:27时，西068运输巷掘进31 m时，本巷及其与石门交叉区域出现强烈冲击显现，微震监测震级为2.4级[20]。

针对上述情况，对掘进工作面开展电磁波CT探测，对掘进工作面后方大面积区域开展地震波CT探测。

由于掘进工作面附近区域动力显现较为明显，这里首先对其开展电磁波CT探测，然后是掘进工作面后方的地震波CT探测。

首先，电磁波CT探测孔均为仰角10°，孔径76 mm;1号测区为掘进工作面超前区域，2号和3号测区分别为掘进工作面后方煤壁0～20 m和20～36 m。

根据实验室实验及现场实测结果表明，煤岩变形破裂向外辐射的电磁波频谱范围一般在1 kHz～1 MHz，与被本文CT探测的电磁波信号不在一个频段，2种电磁波间不会相互干扰[16-19]。

电磁波探测区域及参数详见图3和表3。

图2 掘进面电磁波CT探测区域及测区布置Fig.2 Electromagnetic wave CT detection area and measuring area layout between two holes in heading face
表3 掘进工作面电磁波CT探测参数Table 3 Parameters of electromagnetic wave CT detection between two holes in heading face测区编号激发孔深度/m接收孔深度/m孔间距/m探测面积/m2激发次数电磁波频率/106 Hz1号13133.545.518212～222号131319280.021012～223号9916144.09012～22 电磁波CT探测技术在实际操作中，分别在探测区域两侧钻取激发孔和接收孔，将
激发设备和接收设备安装至孔中，激发点发射高频电磁波，位于另一侧的接收点进行接收。

采用完备测量，即激发孔和接收孔中从孔口开始向里间隔1 m设定观测点，所有观测点均有组合。

探测中均采用“右发左接”的方式。

采用地震波CT对掘进工作面后方50～338 m区域进行探测。

综合考虑西068掘进工作面动力显现区域分布情况及周围开采环境，根据设备探测能力，设计2轮
探测，两轮探测接收探头均布置运输巷，共实施爆破震源54个。

运输巷探测范围265 m，回风巷探测范围275 m，探测面积23 490 m2。

探测系统布置如图3所示。

图3 西068掘进工作面后方区域地震波CT探测系统布置Fig.3 Layout map of seismic wave CT detection system in west 068 heading face
4.1.2 掘进工作面的电磁波CT危险性评价
图4～6为西068运输巷掘进工作面冲击危险性指数分布云图及等值线图[13]。

依据表2的分类标准，划定冲击危险区域。

以下具体分析:
图4 掘进工作面超前区域冲击地压危险性指数Ds分布Fig.4 Distribution diagram of impact hazard index Ds in advance area of heading face
图5 掘进工作面后方0～20 m帮部冲击地压危险性指数Ds值分布Fig.5 Distribution diagram of impact hazard index Ds in the 0～20 m section of heading face back
图6 掘进工作面后方20～36 m帮部冲击地压危险性指数Ds分布Fig.6 Distribution diagram of impact hazard index Dsin the 20～36 m section of heading face back
(1)由图4可知，掘进工作面超前12 m范围内，局部区域具有弱冲击危险，其中10 m范围内危险区域面积相对较大，10 m范围外危险区域面积逐渐减小，说明10 m范围内煤体较为完整，在掘进面超前支承压力的影响下，形成了应力集中区。

(2)由图5和图6可知，在支护范围内，掘进面帮部在后方0～4 m，具有小面积的中等冲击危险区域，说明该区域在掘进扰动下煤体破碎性较高，裂隙发育，支护损伤程度较高。

而帮部在掘进工作面后方4～30 m，支护内围岩仅有小面积的弱冲击危险区，表明该范围支护内围岩具有较高的完整性，裂隙发育程度较低，现有支护方式对帮部围岩的控制较好，围岩具有较强的抗冲击能力。

支护外围岩，在掘进工作面后方0～30 m帮部，深度大于8 m具有弱冲击危险，说明该范围帮部支护外煤体完整性较高，卸压不充分，在巷道开挖引起的侧向支承压力影响下，产生了应力集中区，储存了较多弹性能。

4.1.3 掘进工作面后方区域的地震波CT危险性评价
通过掘进工作面附近区域的电磁波CT探测确定了危险区域及危险等级，在对冲击危险区域开展卸压措施后，对后方区域进行地震波CT探测。

图7为西068工作面已掘部分探测区域冲击危险性指数C分布图，图中以蓝色到红色从小到大来代表探测区域内冲击危险性指数，由图中可见测区内的煤岩体总体上处于弱冲击危险等级，这些区域多位于探测区域的中部;部分区域处于中等冲击危险等级，这些区域多与工作面运输巷和回风巷相邻。

图7 西068掘进工作面后方区域冲击危险性指数C分布Fig.7 Distribution diagram of impact hazard index C for rear area of west 068 heading face 图8 西068工作面巷道冲击危险区范围Fig.8 Scale of the impact danger zone in upper and lower roadway of west 068 working face
采掘空间冲击危险性的判定是冲击地压危险性评价的重点。

根据实测煤岩层冲击危险区域分布范围及至巷帮的距离，划定了西068工作面掘进巷道冲击危险区域，如图8所示，冲击危险等级:ⅲ-1，ⅲ-2为中等，ⅲ-1，ⅲ-2，ⅲ-3，ⅲ-4，ⅲ-5和ⅲ-6为弱。

对比震波CT探测结果发现，在采取措施后，危险性变弱，具有较好的卸压效果，保障了巷道的安全掘进。

4.2 回采工作面的工程实践
4.2.1 回采工作面概况及探测方案
试验地点为陕西长武亭南煤矿207工作面回风巷，探测区域埋深529 m，主采4号煤层平均厚度17.3 m，倾角7°，具有强冲击倾向，煤层顶板为抗压强度28.7 MPa的细粒砂岩，具有强冲击倾向，底板为砂质泥岩，具有弱冲击倾向;207工作面长度200 m，207回风巷为邻空巷道，巷宽为5.2 m，其一侧区段煤柱宽度30 m，采掘平面图如图9所示。

207回风巷掘出后在区段煤柱内布置巷道的区域，出现了动力显现，威胁了矿井的安全生产，防冲形势较为严峻。

针对上述情况，对煤柱区开展电磁波CT探测，对工作面内开展地震波CT探测，由于区段煤柱侧动力显现较为明显，这里首先对其开展电磁波CT探测，然后是工作面的地震波CT 探测。

图9 207工作面采掘工程平面图Fig.9 Map of mine working of 207 workingface
首先，对煤柱区开展电磁波CT探测，包括两类区域:① 区段煤柱内多巷道布置区域;② 区段煤柱内未布置巷道区域。

具体的探测参数如图10和表4所示。

然后，对工作面开展地震波CT探测。

根据矿井实际情况及设备探测能力将207工作面分6个区域。

共设计激发震源444个，获得震波数据8 957道，探测走向范围约为2 545 m，探测面积约509 295 m2。

现场探测时，放炮端布置在运输巷和开切眼区域，接收端布置在回风巷。

图10 煤柱区孔间电磁波CT探测区域及测区布置Fig.10 Electromagnetic wave CT detection area and layout in coal pillar area
表4 孔间电磁波CT探测参数Table 4 Parameters of electromagnetic wave CT detection between two holes in heading face测区编号激发孔深度/m接收孔深度/m孔间距/m探测面积/m2激发次数电磁波频率/(106 Hz)2号1313153001
6002～323号1313153001 6002～324号77121809002～32
4.2.2 回采面的地震波CT危险性评价
图15为207工作面1～3号探测区域、4～6号探测区域冲击危险性指数C分布图，由图中可见测区内的煤岩体总体上处于中等冲击危险等级。

根据实测煤岩层冲击危险区域(图15)分布范围，划定了207工作面探测区域回采
巷道冲击危险区域，如图16所示。

巷道冲击危险等级:ⅳ-1，ⅳ-2，ⅳ-3，ⅳ-4，ⅳ-5，ⅳ-6为强;ⅲ-1，ⅲ-2，ⅲ-3，ⅲ-4，ⅲ-5，ⅲ-6，ⅲ-7，ⅲ-8，ⅲ-9，
ⅲ-10为中等;ⅱ-1，ⅱ-2，ⅱ-3，ⅱ-4，ⅱ-5，ⅱ-6，ⅱ-7，ⅱ-8，ⅱ-9，ⅱ-10，ⅱ-11，ⅱ-12，ⅱ-13，ⅱ-14，ⅱ-15为弱。

4.2.3 区段煤柱的电磁波CT危险性评价
由地震波CT的探测结果，区段煤柱内布置巷道的区域以及巷道交叉口附近区域冲击危险性较高。

图11～16为207回风巷一侧区段煤柱的冲击危险性指数Ds分布云图及等值线图。

依据冲击危险等级分类标准，将各冲击危险区域划定。

以下对几个探测范围具体分析:
图11 1号区域(回风巷回采帮)冲击地压危险性指数Ds分布Fig.11 Distribution diagram of impact hazard index Ds in No.1 area
图12 2号区域(回风巷30 m煤柱帮)冲击地压危险性指数Ds分布Fig.12 Distribution diagram of impact hazard index Ds in No.2 area
图13 3号区域(回风巷回采帮)冲击地压危险性指数Ds分布Fig.13 Distribution diagram of impact hazard index Ds in No.3 area
图14 4号区域(回风巷18 m煤柱帮)冲击地压危险性指数Ds分布Fig.14 Distribution diagram of impact hazard index Ds in No.4 area
由图11可知，回风巷30 m煤柱区的回采帮支护区内仅具有较小的冲击危险区，而支护区外冲击危险区域至煤壁的距离大于9 m，说明回采帮支护结构完整性较
好，抗冲击能力较强，支护区外没有弹性能聚集，整体上冲击危险性较低;由图12可知，煤柱帮在支护区内和支护区外均分布冲击危险区域，说明煤柱帮支护结构的完整性较低，抗冲击能力较低，而支护区外煤体较为完整，聚集较多的弹性能。

由图13可知，回风巷18 m煤柱区回采帮的冲击危险区域多分布于支护区外，支护区内仅具有小面积的弱冲击危险区，说明回采帮冲击危险性较低，但与30 m煤柱区的回采帮相比冲击危险性有所升高;由图14可知，18 m煤柱帮在支护内和支护区外具有大面积的弱及中等冲击危险区域，说明相对于其他探测区域，该区域支护结构完整性更差，支护区外集聚的弹性能更高，且支护区外的冲击危险区域至煤壁的距离仅有4.5 m，由于距离较小，该危险区能量释放导致冲击地压的概率大幅升高。

综上，18 m区段煤柱不仅能够导致煤柱帮冲击危险性的大幅升高，还能导致回采帮冲击危险性的升高，煤柱帮冲击危险性的增幅更大。

4.3 应用效果
通过实践发现，地震波及电磁波CT联合探测的评价方法能够实现掘进工作面及其后方大面积区域，回采面及一侧区段煤柱区的危险程度确定及危险区域划分。

根据评价结果采取治理措施后，西068掘进工作面和207回采工作面危险区域附近发生了强矿震，虽然现场震感较为强烈，但巷道支护几乎没发生任何破损，验证了本文提出的评价方法的有效性。

图15 207工作面冲击危险性指数C分布Fig.15 Distribution diagram of impact hazard index C for 207 working face detection area
图16 207工作面回采巷道冲击危险区范围Fig.16 Scale of the impact danger zone in upper and lower roadway of 207 working face
5 讨论
(1)2种CT探测方式间的联合关系。

2种CT探测方式间联合关系是根据2种方法基于同一力学基础以及危险性划分标准，2种方法的探测结果可以相互比较，并且。