靳红梅, 党琪, 李洪安, 李占利
(西安科技大学 计算机科学与技术学院, 陕西 西安 710600)
国家能源局出台的相关政策指出,通风系统的安全稳定是矿井生产的必备前提,必须做到通风质量合格,工作环境合理,从而保证煤矿安全生产[1]。对于矿井通风系统而言,风速过大虽利于防止瓦斯积聚、稀释有毒有害气体,但会造成电能浪费,还可能造成漏风而引起火灾,或扬起地面粉尘等,危及工作人员健康或发生风源性事故。而风量不足会造成有毒有害气体浓度超标。因此,需对通风系统的通风质量进行动态评价,从而保证通风的安全持久化[2]。
针对煤矿通风的安全评价,国内外学者从全局角度出发,在评价指标和评价方法2个方面进行了深入研究。在评价指标方面,苏盈盈[3]用SVM(支持向量机)对不同指标组合的评价效果进行筛选,将最佳评价效果的指标组合作为精简后的指标体系,完成了小样本数据下矿井通风系统的16项评价指标体系的建立。芦志刚[4]从安全稳定、经济合理性角度出发建立了一套完整的通风系统评价指标体系,共建立了40个指标,完成了对通风系统的综合评价。张聪聪等[5]在全面分析煤矿通风系统可靠性影响因素的基础上,采用主成分分析法将13个评价指标简化为4个综合指标,完成了煤矿通风系统的可靠性评价。李强[6]为进一步提高煤矿企业通风安全监管,通过层次分析法建立了包含18项指标的通风系统安全稳定性评价指标体系。在评价方法方面,王海青[7]应用FMEA(Failure Mode and Effects Analysis,故障模式和影响分析)分析法对矿井通风系统进行了整体评价。孙佳保[8]采用主成分分析法将各种评价指标转变为相应的线性组合,进而利用构建的综合评价函数的方法评价矿井通风系统。刘洋[9]建立了基于层次分析法的煤矿通风系统评价体系,从环境、设施、管理3个方面提出了通风网络结构合理性、通风地点的温度、通风设备的可靠性、人员素质等10个具体指标。Jia Tinggui等[10]应用多元回归分析法分析了矿井通风系统的稳定性。毕娟等[11]利用集成模糊分析法和灰靶决策方法的评价方法对矿井通风系统进行可靠性评价。
综上可知,对矿井通风系统进行评价,有助于促进煤矿安全生产的持续推进。但现有的矿井通风系统评价方法是从矿井通风系统全局视角出发,对其进行整体性的评价,所得出的评价结论往往是静态的[12]、综合性的,得出的结果太过笼统,无法精准定位出矿井中用风范围内的通风情况并发现其存在的问题,也不能够实时反映当前时刻的通风质量状态。针对以上问题,本文建立了基于FCE-AHP(Fuzzy Comprehensive Evaluation Method-Analytic Hierarchy Process,模糊综合评价法-层次分析法)的矿井通风质量动态评价模型,用以对煤矿的通风质量进行实时评价。首先根据通风需求和作用将通风区域划分为进风区段、回风区段、采煤工作面、掘进工作面、硐室5个区域,根据每一种用风区域类型确定其对应的评价指标体系;然后根据评价指标体系获取对应测风点传感器实时监测到的数据;最后通过所建模型,依据实时监测到的数据对通风质量进行实时评价,确定煤矿通风质量的安全级别。该模型的评价结果是基于实时监测数据计算得出,能够及时反映通风质量状态,实时监测通风质量,保证矿井通风时刻安全。
为了能够从通风系统的评价结果中得到更多有效信息,更加精准地确定危险所在地,本文从通风系统各测风点出发展开研究,即对各个传感器所在位置进行划分,每个测风点所在区域为风流较为平稳的一个巷道或一片区域,其所测数据能反映该区域的整体通风情况,可作为一个原子用风区域,通过分析每一个原子用风区域的通风用途和需求,确定用风区域类型,将通风用途和需求相同的归为一类。在确定用风区域类型之后,根据每一种用风区域类型确定其对应的评价指标体系,可以更全面实现煤矿通风质量的评价。
根据通风需求和作用将通风区域划分为进风区段、回风区段、采煤工作面、掘进工作面、硐室5个区域。其中由于采煤工作面中上隅角区域的通风不太顺畅,其通风指标含量与采煤工作面中其他工作地点的要求准则有所不同,故该区域的危险评判界限不同,需将其进行单独分析,这里将采煤工作面分为2个子区域:一个为上隅角区域,一个为回风巷区域。基于不同用风区域的确定,分别考虑各区域对通风质量的不同影响因素,进而建立各区域相应的评价指标体系。具体的通风系统通风质量划区域评价层次指标体系如图1所示。
图1 通风系统通风质量划区域评价层次指标体系Fig.1 Ventilation quality region evaluation hierarchy index system of ventilation system
矿井通风系统是一个复杂随机的动态系统,对其通风质量评价所涉及的因素繁多,且通风质量是随煤矿生产推进而变化的。只依据考核评级办法完成对通风质量的打分不能满足对通风质量评价的实时性,需通过对矿井通风进行动态监控,根据实时监测数据进行动态评价,才能真实反映生产过程中煤矿的实时通风质量是否合格,进而有效促进矿井通风安全的方案调整,做到防患于未然。AHP、FCE等方法能够处理含有定性因素、定量因素的评价问题,使评价过程规范化,以提高评价的科学性。AHP能够将复杂对象的影响因素进行系统化、层次化的分析,计算简单、适用性广泛。FCE可以对多因素、多变量、多层次的复杂系统进行评价,能够对模糊安全状态的评价对象进行合理评价。故本文选择AHP与FCE相结合的方法完成通风质量的动态评价。
权重是指各评价因子对评价目标所起作用的重要程度,运用AHP可有效将原本复杂的定性评价过程进行合理量化,将无序的要素按照目标层、准则层和方案层归纳出有序的递阶层次结构[13]。本文采用AHP进行各指标权重的确定,根据9级标度确定相对重要程度[14],构建成对比较判断矩阵。9标度指标重要程度判断见表1。
表1 9标度指标重要程度判断Table 1 9 scale index importance judgment
为了确保利用AHP所得评价指标权重的基本合理性,需对成对比较矩阵进行一致性检验,即对于成对比较矩阵A,有
AW=λmaxW
(1)
式中:W为所对应的特征向量;λmax为最大特征值。
一致性指标CI为
(2)
式中N为矩阵阶数。
一致性比例CR为
(3)
式中RI为平均一致性指标,具体取值见表2。
表2 平均一致性指标值Table 2 Average consistency index values
当CR<0.1时,认为成对比较矩阵的一致性是可接受的。否则,需对成对比较矩阵作适当修改。
在确定各区域指标权重之后,本文按照优秀、良好、合格、危险4个等级对指标进行划分。基于矿井通风安全的研究,由于大部分通风质量指标范围具有上下限,即认为监测数据超出上下限范围的,视为危险级别。通风质量中各指标的量纲不同,含量标准不同,故其4个等级的确定标准范围也各不相同,又因为煤矿安全的相关文件中未明确规定出煤矿通风质量的优良级别界限,故本文通过函数分布划分通风质量的不同指标的各级别界限。基于对通风质量各指标监测数据的研究,可近似认为指标界限分布满足正态型分布,可设函数为
(4)
式中:x为指标变量值;a,b为常数,a≥0,b>0。
在等级划分函数确立的基础上,完成对各个指标的级别划分。根据通风质量各动态指标的危险界限,确定每个指标不同等级的界限。具体级别结合《煤矿安全规程》中的危险界限进行计算,完成各评价指标的标准分级。
隶属度函数的选择直接影响隶属度的大小,进而影响评价结果的合理性。由于通风质量的评价指标均为定量指标,柯西分布隶属度函数的计算过程简单,函数图像变化合理[15],故隶属度函数选取柯西分布隶属度函数。隶属度函数主要是用来量化不同等级下的指标,其函数类型的选择与指标的特点有关,需要依据具体因素进行分析。一般情况下,评价为极差和极好的2种情况会选用偏小或偏大型隶属度函数,其余情况可选择居中型隶属度函数。
偏小型隶属度函数为
(5)
居中型隶属度函数为
(6)
偏大型隶属度函数为
(7)
式中:f(x)为输出隶属度;α,β为调整函数的参数,其中,β取值为2,α一般由实践经验确定。
将类似柯西分布的隶属度函数作为确定元素隶属度的基础,对有区间范围的指标选用式(6),两极分布的选用式(5)或式(7)。
对于居中型隶属度函数有
(8)
式中xmax,xmin为对应级别相应参数的上下边界值。
对于偏大型或偏小型分布隶属度函数,结合煤矿通风的危险等级,当指标达到危险级别时,即具有“一票否决权”,也就是对于超出上下界限范围的监测数据,直接认定其危险级别的隶属度值为1,即f1(x)=1,为了不削弱危险程度的监测值的评价结果,故将煤矿通风质量的两极分布隶属度函数设为
f1(x)=1(x>xmax,x (9) 本文通过FCE对通风质量进行综合评价,其计算步骤如下: Step1 确定评级对象因素集:U={u1,u2,…un}(n=1,2,…),un为指标因素。 Step2 确定评语集:V={v1,v2,…vm}(m=1,2,…),vm为评语。 Step3 单因素评价:讨论单个因素对于评语集V的隶属度,得到第i个区域的单因素评价矩阵Ri为 (10) 式中rinm为计算得到的权重。 Step4 综合评价:综合考虑各指标隶属度完成评价对象的等级量化,再通过置信度识别原则得到最终的综合评价结果。 Q=wi·Ri= (11) 式中:wi为第i个区域的权重向量;·为模糊算子;q为矩阵乘积计算出来的综合权重。 东滩煤矿位于山东省兖州市,井田位于兖州煤田东部,井田范围大,通风方式为对角式通风,主要通风方法为抽出式通风。本文结合该煤矿数据,选择63上06工作面的采煤工作面回风巷进行通风质量评价,该通风线路各测点所在区域风流稳定,巷道断面规整,均符合测点所应遵循的原则。 63上06工作面通风系统的具体路线:进风区段(副井→井底车场→南翼轨道大巷→南翼辅运下山→南翼辅运巷)→其他用风区域(南翼第二胶带运输巷南段→辅助回风巷联络巷)→工作面回风巷(南翼回风巷南段→南翼回风巷)→掘进工作面(63上06采区轨道巷)→回风区段(支架联络巷→南翼总回风巷→西翼总回风巷→西风井底→西风井)。本文以南翼回风巷数据为例进行评价,具体数据见表3。 表3 南翼回风巷某一时刻通风质量指标数据Table 3 Ventilation quality index data at a certain moment in the south wing return airway 对于通风质量的评价,首先需完成指标权重的计算,即由多个煤矿通风安全专家对各区域指标进行重要度评判,取其均值,得到各指标间两两对比的重要度判断,再通过AHP计算该区域各指标权重。基于图1所建立的通风系统通风质量划区域评价层次指标,完成各区域的指标权重计算,判断矩阵记为Bi(i=1,2,…,6,分别表示采煤工作面的上隅角、回风巷、掘进工作面、硐室、进风区段、回风区段)。同时验证各判断矩阵Bi是否通过一致性检验,通过则计算权重,否则对判断矩阵进行修改,直到通过一致性检验,计算得出最终权重结果wBi,见表4。 在得到表4中的各指标的指标权重后,根据式(1)—式(3)对判断矩阵进行一致性检验,矩阵均通过了一致性检验,所得权重合理。从表4可得到采煤工作面回风巷各指标权重为wB2=[0.155 5,0.024 7,0.018 3,0.222 3,0.312 1,0.035 0,0.050 7,0.107 5,0.073 9]。 结合表3中东滩煤矿南翼回风巷中当前时刻所检测到的通风质量指标数据,计算当前时刻通风质量各指标在评价等级中的隶属度矩阵,完成该区域当前时刻的通风质量评价。 表5 采煤工作面回风巷危险级别划分Table 5 Classification of danger levels of return air roadway in coal working face 在确定该区域指标分级的基础上,确定隶属度权重,根据划分指标的级别界限,结合隶属度函数确定理论求得矿井通风采煤工作面回风巷区域的各个通风质量指标的隶属度函数。 风速的隶属度函数为 (12) 温度的隶属度函数为 (13) 湿度的隶属度函数为 (14) 风量供需比的隶属度函数为 (15) 瓦斯浓度的隶属度函数为 (16) CO2浓度的隶属度函数为 (17) CO浓度的隶属度函数为 (18) O2浓度的隶属度函数为 (19) 粉尘浓度的隶属度函数为 (20) 结合表3中传感器检测到的当前时刻数据,根据所建立的隶属度函数计算该区域指标在4个级别中的隶属度,归一化后所得各个指标的隶属度级别矩阵为 (21) 针对采煤工作面回风巷区域,此时w1即为wB2,计算w1·w2得 w=[0.589 9,0.204 4,0.152 0,0] (22) 根据置信度识别原则,选择置信度为μ=0.6,即有0.589 9+0.204 4=0.794 3>0.6,可得评价集V为良好,故南翼回风巷此刻的通风质量评级为良好,与实际情况相符。 (1) 分析了煤矿通风质量的影响因素,通过对《煤矿安全规程》及煤矿通风资料的研究分析,从通风需求和作用2个方面出发,将煤矿通风区域划分为进风区段、回风区段、采煤工作面、掘进工作面、硐室5个区域;全面考虑不同区域通风质量的影响因素,建立了不同区域的矿井通风质量动态评价指标体系;通过分析煤矿测风点实时监测数据的规律,结合《煤矿安全规程》中的指标界限,用符合监测数据分布的正态型函数对评价等级范围进行划分,提出了基于FCE-AHP的动态评价模型,对区域通风质量进行动态评价。该模型选取柯西分布隶属度函数,选择置信度判别原则确定评价级别,提高了通风质量评价的客观性。区域通风质量的动态评价有利于精准判断通风系统中的通风质量的危险源所在地。 (2) 利用基于FCE-AHP的动态评价模型对东滩煤矿的采煤工作面回风巷区域通风质量进行评价,评价结果为良好,煤矿经专家评定得到的通风质量为良好,评价结果与实际情况相符,说明评价模型的建立合理,能够对煤矿不同通风区域的通风质量进行动态评价,进而实现对矿井通风质量的综合评价,为改善煤矿通风系统质量提供依据,对煤矿通风安全生产的推进具有实用价值。2.4 综合评价
3 应用与分析
3.1 东滩煤矿概况
3.2 东滩煤矿南翼回风巷通风质量动态评价
4 结论