1. 引言
煤矿瓦斯爆炸是在多种因素的共同作用下形成的极为复杂的灾害事故,同时也是威胁我国煤矿安全生产的重大灾害事故之一[1]。据不完全统计,近十年来我国煤矿每年都会发生重大及以上级别的瓦斯爆炸事故,且在现阶段中,90%以上的煤矿都存在瓦斯爆炸安全隐患[2]。随着高新技术的发展,煤矿的开采能力不断得以提升,煤层的开采深度随之增加,开采范围也在持续扩大[3]-[6],这对瓦斯爆炸风险的有效管控带来了巨大挑战,同时严重威胁井下作业人员的生命安全。由此可见,积极地开展瓦斯爆炸风险指标识别及风险评估,及时采取有效控制措施,从源头遏制事故发生,对我国煤矿安全生产有着重要意义。
诸多学者针对煤矿瓦斯爆炸风险评估进行了大量研究,例如:李润求等[7]运用粒子群优化算法对SVM模型参数进行优化,构建出瓦斯爆炸灾害风险识别的PSO-SVM模型。Zhang Qing [8]基于灾害系统理论,将煤矿瓦斯爆炸灾害分为孕灾环境子系统、致灾因子子系统及承灾体子系统,同时对其建立了随机森林模型。皮子坤等[9]结合前景理论与模糊综合评价理论以及区间数排序理论构建煤矿瓦斯爆炸风险评价模型,并通过实例对该模型进行验证。李敏等[10]基于模糊贝叶斯网络提出一种煤矿瓦斯爆炸风险评估方法,根据三角模糊数评估风险指标的条件概率,利用贝叶斯计算瓦斯爆炸发生概率,推理分析出瓦斯爆炸的事故成因。鲁锦涛等[11]从经典瓦斯爆炸事故案例中提取影响因素,基于集成灰色系统理论,构建灰色–物元风险评估模型,明确风险因素之间的作用路径。景国勋等[12]通过SPA-VFS方法,从人、机、环、管四个方面对煤矿瓦斯爆炸风险进行了全面评估,耦合模型评价结果与实际相符。Tauziede等[13]结合矿井瓦斯和作业环境的各方面调查数据,对矿井瓦斯爆炸事故风险进行全面评估。Pejic等[14]采用精细化方法构建矿井瓦斯爆炸风险指数,利用概率计算方法得出煤矿瓦斯爆炸的风险概率及危险程度。上述方法在实际工程中都取得了一定的效果,但由于导致瓦斯爆炸的因素众多,且煤矿安全生产系统中的大量信息具有模糊性和不确定性[15],加之无法对一些定性指标进行定量评价,这就使得所构建的模型对此类信息进行评价具有局限性。因此,需要探寻出一种既适用于煤矿前期勘察设计又能服务于巷道掘进施工过程,同时具有较好准确性、可靠性和实用性的多因素瓦斯爆炸风险评估模型。
王光远教授在1990年首次提出未确知测度理论,在分析不确定因素信息的同时,还能进行量化处理,在解决有序切割的问题上有着显著优势[16]。目前未确知测度理论广泛应用于巷道围岩稳定、冲击地压、隧道坍塌[17]-[22]等领域,并取得了较好的效果。鉴于此,本文将能够充分考虑瓦斯爆炸元素相关性的区间层次分析法(Interval Analytic Hierarchy Process, IAHP)和未确知测度理论相结合,从瓦斯爆炸事故的致灾机理和形成条件出发,基于危险源理论和事故致因理论,建立瓦斯爆炸风险评估指标体系。通过IAHP计算出指标权重,运用未确知测度的基本理论和计算规则,构建一种煤矿瓦斯爆炸风险评估模型,并将该模型应用于实际煤矿的风险评估,以校验其合理性与有效性,以期为煤矿瓦斯爆炸风险评估提供新的方法和思路。
2. 风险评估理论
2.1. IAHP
IAHP是AHP的一种改进和扩展,该方法通过区间数对比两两因素之间的重要程度,在建立判断矩阵中采用数值区间替代单数值,有效解决了决策者因信息掌握不全而造成判断矩阵不准确的问题。IAHP的计算过程与AHP一致,即首先确定与分析问题有关的影响因素,然后根据是否具有某些特征把这些因素聚集成组,将它们的共同特征作为系统中新层次的因素,同时这些因素也根据另外一组的特性构成更高层次的因素,最终形成单一的最高因素即分析目标,由此构建层次分析结构模型[23]。本文采用主观赋权法——IAHP确定煤矿瓦斯爆炸各影响因素权重,步骤如下[24]。
1) 根据煤矿瓦斯爆炸建立层次分析结构。层次分析模型应包含目标层、准则层和指标层。
2) 确定判断矩阵。通过多位专家评判,采用二元对比法,明确每个层次中各元素的相对重要性。
(1)
其中:A为区间矩阵,
和
为A所确定的矩阵;
和
分别为矩阵
和矩阵
的元素,且
,
;
表示利用二元对比法得到的元素
与元素
之间的相对重要性;n表示指标个数。
3) 区间判断矩阵的一致性检验。
(2)
(3)
(4)
式中:
为最大特征值;
为特征向量M的元素;k和
均为修正系数;
和
分别为矩阵
与矩阵
的最大特征根所对应的特征向量;Q为所计算出的权重范围。
4) 确定各影响因素权重。
(5)
式中
为层次结构模型中各指标的权重,通过式(2)~式(5)检查权重的分布是否在合理范围。
2.2. 未确知测度理论
假设有n个待评价的对象
,则评价对象空间
,每个评价对象包含有m个评价指标
,则指标空间
,第i个评价对象关于第j个评价指标的测度值为
。对每个测度值设定c个评价等级
,则等级空间
,若满足
或
,则称
是评价等级空间U的一个有序分割类[25]。
2.2.1. 单指标测度
令
表示指标测度值
属于第k个评价等级
的程度,当
满足:
(6)
(7)
(8)
则称
为未确知测度。其中式(8)表示“归一性”,式(9)表示“可加性”。所有测度值所构成的矩阵称为单指标未确知测度评价矩阵。
(9)
2.2.2. 多指标综合测度
令
表示
属于第k个评价类
的程度,则
(10)
由于
,且
,所以
是未确知测度,矩阵
(11)
称为多指标综合测度评价矩阵,
称为评价对象
的综合测度评价向量。
2.2.3. 置信度识别准则
在求出多指标综合测度评价矩阵的条件下,如果
是评价等级空间U的一个有序分割,即可引入置信度识别准则。假设置信度为
(
,且常取0.6或0.7),
,令
(12)
则认为评价对象
属于第
个评价类
。
3. IAHP-未确知测度综合评价
3.1. 评价指标体系的建立
建立瓦斯爆炸风险指标体系是对瓦斯爆炸进行风险评估的首要条件,也是评估结果可靠程度的主要影响因素。基于煤矿瓦斯爆炸事故的致灾机理和形成条件,将危险源理论和事故致因理论相结合,考虑“人–机–环境–管理”四个因素建立瓦斯爆炸风险评估指标体系[26] [27],见图1。
Figure 1. Hierarchical structure model of risk assessment index system
图1. 风险评估指标体系层次结构模型
3.2. 评价等级的确定
根据常用的分级标准与专家经验及评价指标的意义,同时考虑具体的实际情况,将煤矿瓦斯爆炸风险等级划分为安全、较安全、一般安全、较不安全和不安全5个评价等级,分别用P1,P2,P3,P4,P5表示。由于煤矿安全生产的特殊性和井下环境的复杂性,加之瓦斯爆炸受多种因素的综合影响,依据各指标的不同属性,其分级标准见表1。
Table 1. Evaluation index grading standard
表1. 评价指标分级标准
指标 |
等级区间 |
P1 (安全) |
P2 (较安全) |
P3 (一般安全) |
P4 (较不安全) |
P5 (不安全) |
X11; X21; X22; X23; X24; X25; X32; X41; X42; X43; X44 |
[90, 95) |
[85, 90) |
[75, 85) |
[65, 75) |
[50, 65) |
X12 |
[15, 16) |
[12, 15) |
[9, 12) |
[6, 9) |
[3, 6) |
X13 |
[120, 180) |
[100, 120) |
[60, 100) |
[30, 60) |
[20, 30) |
X14 |
[16, 20) |
[12, 16) |
[8, 12) |
[4, 8) |
[2, 4) |
X31 |
[2, 4) |
[4, 5) |
[5, 7) |
[7, 10) |
[10, 12) |
X33 |
[4, 5) |
[5, 6) |
[6, 8) |
[8, 10) |
[10, 12) |
X34 |
[95, 100) |
[90, 95) |
[85, 90) |
[80, 85) |
[75, 80) |
X35 |
[9, 12) |
[6, 9) |
[5, 6) |
[4, 5) |
[2, 4) |
X36 |
[8, 10) |
[10, 14) |
[14, 15) |
[15, 18) |
[18, 20) |
X37 |
[0.4, 0.5) |
[0.5, 0.6) |
[0.6, 0.7) |
[0.7, 0.8) |
[0.8, 0.9) |
3.3. 评价指标权重的确定
基于评价指标分级标准,同时咨询多位专家及经验丰富的煤矿作业人员,综合考虑其所给意见,采用二元对比法构造出区间判断矩阵,一级指标区间判断矩阵如表2所示,同理,可构造二级指标区间判断矩阵并确定权重,在此不再叙述。根据式(2)~式(5)分别计算出各层次的指标权重,如表3所示。根据各指标的单权重以及所属层次的权重,得出20个评价指标的综合权重,如表4所示。
Table 2. First-level index judgment matrix
表2. 一级指标判断矩阵
一级指标 |
X1 |
X2 |
X3 |
X4 |
X1 |
[1, 1] |
[1, 2] |
[1, 2] |
[1/2, 1] |
X2 |
[1/2, 1] |
[1, 1] |
[1, 2] |
[1/2, 1] |
X3 |
[1/2, 1] |
[1/2, 1] |
[1, 1] |
[1/2, 1] |
X4 |
[1, 2] |
[1, 2] |
[1, 2] |
[1, 1] |
Table 3. First-level index weight
表3. 一级指标权重
一级指标 |
|
X1 |
0.2673 |
X2 |
0.2248 |
X3 |
0.1890 |
X4 |
0.3179 |
一致性检验 |
k = 0.8715 < 1,β = 1.1265 > 1,满足一致性检验 |
Table 4. Comprehensive index weight
表4. 综合指标权重
一级指标 |
一级指标权重 |
二级指标 |
二级指标权重 |
综合权重 |
X1 |
0.2673 |
X11 |
0.1570 |
0.0420 |
X12 |
0.3234 |
0.0864 |
X13 |
0.1320 |
0.0353 |
X14 |
0.3845 |
0.1028 |
X2 |
0.2248 |
X21 |
0.4211 |
0.0947 |
X22 |
0.2756 |
0.0620 |
X23 |
0.1644 |
0.0370 |
X24 |
0.0849 |
0.0191 |
X25 |
0.0514 |
0.0116 |
X3 |
0.1890 |
X31 |
0.2124 |
0.0402 |
X32 |
0.0944 |
0.0178 |
X33 |
0.1604 |
0.0303 |
X34 |
0.2124 |
0.0402 |
X35 |
0.0238 |
0.0045 |
X36 |
0.0295 |
0.0056 |
X37 |
0.2647 |
0.0500 |
X4 |
0.3179 |
X41 |
0.1204 |
0.0383 |
X42 |
0.1801 |
0.0572 |
X43 |
0.4172 |
0.1326 |
X44 |
0.2818 |
0.0896 |
3.4. 建立未确知测度函数
根据表1中各指标的分级标准,结合单指标函数的定义选用直线法建立指标测度函数,本文给出通风系统、煤层瓦斯含量和风量供需比测度函数见图2。
(a) 通风系统测度函数
(b) 煤层瓦斯含量测度函数
(c) 风量供需比测度函数
Figure 2. Evaluation index measure function diagram
图2. 评价指标测度函数图
4. 实例应用分析
运用IAHP-未确知测度评价模型对贾宝山[28]等所作研究中的数据样本(表5)进行评价,以验证该模型的准确性。
Table 5. Evaluation of sample data
表5. 评价样本数据
工作面编号 |
X11 |
X12 |
X13 |
X14 |
X21 |
X22 |
X23 |
X24 |
X25 |
X31 |
X32 |
X33 |
X34 |
X35 |
X36 |
X37 |
X41 |
X42 |
X43 |
X44 |
1# |
95 |
12.18 |
183 |
26.05 |
95 |
95 |
95 |
95 |
95 |
11.21 |
95 |
6.36 |
105 |
6.5 |
22.34 |
0.51 |
95 |
95 |
95 |
90 |
2# |
95 |
11.56 |
163 |
21.32 |
90 |
90 |
90 |
90 |
90 |
10.16 |
90 |
5.10 |
100 |
7.2 |
21.12 |
0.56 |
90 |
95 |
95 |
90 |
3# |
90 |
10.84 |
151 |
19.58 |
95 |
95 |
95 |
90 |
90 |
12.45 |
85 |
9.53 |
100 |
6.6 |
20.51 |
0.76 |
90 |
95 |
90 |
90 |
4# |
90 |
13.31 |
193 |
24.50 |
95 |
95 |
95 |
95 |
90 |
16.01 |
95 |
7.93 |
105 |
8.8 |
22.34 |
0.41 |
95 |
95 |
95 |
90 |
5# |
95 |
11.85 |
171 |
23.34 |
90 |
95 |
95 |
95 |
95 |
6.12 |
85 |
4.01 |
98 |
9.5 |
19.89 |
0.68 |
90 |
95 |
90 |
90 |
以1#工作面为例,将各因素取值代入测度函数图中,根据式(9)可得到1#工作面人为因素的单指标测度矩阵,同理可得到设备因素、环境因素和管理因素的测度矩阵。
结合区间层次分析法所确定的各指标权重,由式(10)求出1#工作面的多指标综合评价向量,代入式(12),其中取
,求出k值,得到煤矿瓦斯爆炸风险等级评估结果见表6。
Table 6. Gas explosion risk level assessment results table
表6. 瓦斯爆炸风险等级评估结果表
工作面编号 |
综合评价指标测度 |
k值 |
等级 |
1# |
0.7355 |
0.1246 |
0.0905 |
0.0160 |
0.0305 |
1 |
P1 |
2# |
0.4843 |
0.3778 |
0.0764 |
0.0497 |
0.0088 |
2 |
P2 |
3# |
0.3643 |
0.3988 |
0.1016 |
0.0867 |
0.0457 |
2 |
P2 |
4# |
0.6869 |
0.1861 |
0.0782 |
0.0000 |
0.0457 |
2 |
P2 |
5# |
0.4167 |
0.4079 |
0.1625 |
0.0074 |
0.0025 |
2 |
P2 |
由表6可知,1#工作面的瓦斯爆炸风险等级为安全,其余4个工作面评估等级均为较安全,在样本数据中,1#工作面的指标值明显高于其他工作面,因此安全水平相对较高。在煤矿开采过程中,5个采煤工作面均未发生瓦斯爆炸事故,在人、设备、环境和管理方面也未发现
导致事故的危险因素,实际情况与评估结果基本吻合,说明运用IAHP-未确知测度模型对煤矿瓦斯爆炸风险进行等级评估是可行的。
5. 结论
1) 本文提出了一种基于IAHP-未确知测度的煤矿瓦斯爆炸风险评估方法,该方法利用IAHP确定煤矿安全生产系统中影响瓦斯爆炸的各评价指标权重,模型简单且操作容易,通过区间数的模糊性及不确定性,使所计算的指标权重更加客观、准确。同时,基于未确知测度函数对定性与定量指标进行评价,均满足“归一性”和“可加性”。两种方法结合,使评价过程更为科学严谨,为煤矿瓦斯爆炸风险评估提供了研究思路和科学依据。
2) 针对瓦斯爆炸事故本身的不确定性及矿井信息的未确知性,从人、机、环、管4个方面构建评价体系,将煤矿瓦斯爆炸风险划分为5个等级,为瓦斯爆炸风险水平的测度提供一定的基础。通过IAHP所计算权重与多指标综合测度评价矩阵的乘积可确定5个工作面的评价结果分别为安全、较安全、较安全、较安全和较安全,该评价结果符合现场工程实际,说明评价方法较为合理,对煤矿瓦斯爆炸风险评估具有一定的应用价值。
3) 在后续研究中,可将该模型结合到更多矿区及不同工作面的现场实际数据中,进一步验证其可靠性,以提高该模型的精确度和普适性,提高所述方法的评价和预测水平。
基金项目
国家自然科学基金项目“矿井多爆源瓦斯爆炸传播特性及热冲击动力学机制研究”(编号:52174178)、“矿井瓦斯爆炸能量释放转化特性及热冲击能量耗损研究”(编号:52174177)。