1. 引言
顾北煤矿地处淮南煤田中部,该井田为全隐伏煤田,含煤系为石炭二叠系,其主要为二叠系山西组煤层段 [1] [2] [3],煤层厚度平均69 m,含煤1层区内1煤厚2.6~10.2 m,平均7.6 m;上距4-1煤层66.5~75.3 m,平均70.6 m;下距1灰14.6~25.4 m,平均约19.1 m;采区中部较薄,两侧逐渐增厚;一般含1~3层泥岩夹矸;煤层发育稳定,属稳定型煤层 [4] [5]。为保证资源及时供应,现需要继续开采1煤层。但1煤层开采受底板灰岩水害威胁,水文地质条件复杂。为此,建立了地面–井下注浆系统。利用导管进行注浆加固和含水层改造从而达到可安全开采的目的。
13521采区工作面主要受底部太原组灰岩岩溶裂隙水影响,其中C3I组灰岩含水层为1煤层底板直接充水含水层,厚度为29.6~35 m,平均厚度为32.3 m。通过底板的构造裂隙与采动裂隙影响工作面,因此采用注浆加固来进行水害防治。本文以顾北煤矿为例,结合AHP对煤层底板区注浆区域综合治理进行模拟评价 [6]。
2. 注浆孔区域概况
13521工作面主体为向东倾斜的单斜构造 [7]。根据工作面的综合历次三维地震数据、,工作面发育断层14条,其中对采掘影响较大的断层2条,分别为FS854和FS960 (如下表1)。13521工作面煤层发育稳定,属稳定型煤层。
![](Images/Table_Tmp.jpg)
Table 1. Three-dimensional seismic exploration and interpretation of faults in 13521 working face
表1. 13521工作面三维地震勘查解释断层
13521工作面地面区域探查治理范围内施工3个钻探主孔、26个钻探分支孔,钻孔钻进及注浆工作于2021年2月27日全部完成。本次探查治理共施工3个钻探主孔,26个钻探分支孔,共完成钻探工程量17,576 m,其结构图如图1所示。
![](//html.hanspub.org/file/5-2890536x7_hanspub.png?20211019085553434)
Figure 1. Schematic diagram of drilling well structure for exploration and treatment of surface area in Nanyi 1 coal mining area
图1. 南一1煤采区地面区域探查治理钻孔井身结构示意图
3. 评价指标的设计
层次分析法,简称AHP。层次分析法通常分为四个步骤:1) 建立层次结构模型;2) 构造比较判断矩阵;3) 计算矩阵的最大特征根及其对应的特征向量;4) 层次排序和其一致性检验 [8] [9]。此次采用分析法综合分析16个因素影响下的注浆效果评价。
根据前述,将13521工作面范围区域探查治理效果评价指标分为钻孔质量、注浆质量和效果检验三大方面。分别对应总目标层、分目标层和准则层 [10] 具体如图2所示。
![](//html.hanspub.org/file/5-2890536x8_hanspub.png?20211019085553434)
Figure 2. Hierarchical analysis structure diagram
图2. 层次分析结构图
3.1. 判断矩阵构造与计算
权重计算以指标层开始逐步进行探究 [11]。为掌握指标层中有16个因素在区域治理综合评价中的所占权重,和各因素之间的关系,引入1~9的区分标度(见表2)。对各个因素进行比较,采用层次分析法,确定各因素的权重,构造出判断矩阵(见表3~7)。后采用AHP软件进行矩阵的判断与求解 [12]。
![](Images/Table_Tmp.jpg)
Table 2. Discrimination scale of each element in the judgment matrix
表2. 判断矩阵中各元素的区分标度
![](Images/Table_Tmp.jpg)
Table 3. Judgment matrix and weights for comprehensive evaluation of regional governance
表3. 区域治理综合评价的判断矩阵与权重
![](Images/Table_Tmp.jpg)
Table 4. Judgment matrix and weights for borehole quality evaluation
表4. 钻孔质量评价的判断矩阵与权重
![](Images/Table_Tmp.jpg)
Table 5. Judgment matrix and weights for grouting quality evaluation
表5. 注浆质量评价的判断矩阵与权重
![](Images/Table_Tmp.jpg)
Table 6. Judgment matrix and weights of effect test evaluation
表6. 效果检验评价的判断矩阵和权重
![](Images/Table_Tmp.jpg)
Table 7. The value table of random consistency index R.I. with different matrix betweennesses
表7. 矩阵介数不同的随机一致性指标R.I.取值表
假定分目标层(B)指标中Bi与准则层中
有联系,建立判断矩阵,如下:
通过各指标的比较确定各指标的优先权重等级,表示指标关系的相对重要性。其评价结果采用如下矩阵表示。
指标权重的确定:
1) 计算矩阵每行所有元素的几何平均数
,可得,
(3.1)
2) 进行归一化
,可得,
(3.2)
3) 计算判断矩阵最大特征值
,其中
为向量Aw的第i个元素。 (3.3)
4) 一致性检验
(3.4)
(3.5)
其中,C.I与R.I是平均随机一致性指标,R.I.取值见上表7所示。
如表8,判断矩阵的C.R.分别为0、0、0.097和0.091均小于0.1,符合上述层次分析法所要求的一致性检验结果标准。
根据以上计算所得各目标层次指标的所占权重如表2~8所示。为了实现更加总体客观的评价,还需要进行总目标层次排序。其计算方式如下所示:
由此计算所得组合权重判断矩阵见表9。
![](Images/Table_Tmp.jpg)
Table 9. Priority weight of grouting effect evaluation index
表9. 注浆效果评价指标优先级权重
基于上述各指标的权重性分析,对各指标的分析如下:
3.2. 计算结果及其评价
1) 钻孔质量评价(B1)
跟层率(C1)其中各分支水平孔实钻跟层率>89%,符合设计中要求。
钻孔孔深(C2)其中各分支水平孔钻孔孔深均>设计孔深。
轨迹偏差量(C3)各水平分支孔实际轨迹与设计轨迹的偏差均<2 m的要求。
2) 注浆质量评价(B2)
注浆材料(C4)其中各指标质量均符合国家GB175-2007标准。
注浆压力(C5)均大于钻孔水平段所在位置奥灰含水层静水压力的1.5倍 [13] [14]。
注浆用量(C6)均小于设计注浆结束吸水率(0.010 L/min·m·m)。
注浆异常情况次数(C7),注浆期间发生严重漏失异常2次。
异常情况及处理(C8)各异常情况均得到及时有效处理。
浆液扩散有效范围(C9),各分支孔浆液设计间距为60 m,结合注浆模拟可知,实际浆液扩散有效半径>30 m。
3) 效果检验评价(B3)
验证孔水压(C10)各验证孔终孔出水水压为0~1.0 MPa,水压值<深部灰岩水压。
验证孔水量(C11)10个验证钻孔,共有6个出水,水量0.2~0.85 m3/h,钻孔出水量较小。
验证孔水温(C12)验证钻孔实测水温值在32℃~37℃之间,比奥灰水温低。
验证孔水质(C13)各验证孔的常规离子含量、TDS及舒卡列夫水化学分类均与奥灰水存在差距。
物探异常区数(C14)物探异常区面积(C15通过槽波地震和无线电波坑透探查,综合确定13521底板有5处低阻区。
低阻区内钻孔出水量(C16)低阻异常区验证钻孔水量0.4~0.85 m3/h。
根据实际情况将区域探查治理效果等级水平划分为四个等级 [15],并作了定性和定量的描述,具体见表10。
![](Images/Table_Tmp.jpg)
Table 10. Classification of grouting effect grade level
表10. 注浆效果等级水平划分
利用以上16个分目标定量化,对注浆评价数值V进行计算,其计算方式如下:
(3.6)
式中:Zi为各指标无量纲化评价值;Ci为各指标的权重值。
结合各指标与实际情况,对区域综合探查治理效果影响指标评价数值计算,V < 0.25,对照表9可知,确定其属于等级I,注浆效果良好,区域探查治理取得了很好效果。
4. 结论
1) 根据实际现场采用物探方法检测注浆效果与模拟效果评价大致一致可以确定AHP分析法在注浆效果评价中是可以实现一定的预测效果。
2) 对于注浆加固改善采区工作面灰岩水害防治取得了很好的效果。注浆效果明显,探查治理效果良好。
3) 利用各因素关系、专家打分、MATLAB计算,通过建立层次分析结构、构造与计算判断矩阵、确定各因素的权重以及建立评价等级等步骤,分步完成对于多因素影响评价模型。根据工作面实际评价效果从而推动模型的优化与分析。