高海拔地区矿井风机状态动态评估
Dynamic Status Evaluation of Main Fans of Mine at High Altitude Region
通讯作者:
收稿日期: 2020-05-25 修回日期: 2020-07-18 网络出版日期: 2021-01-29
基金资助: |
|
Received: 2020-05-25 Revised: 2020-07-18 Online: 2021-01-29
作者简介 About authors
王利鹏(1994-),男,河南洛阳人,硕士研究生,从事风险分析与评估、安全评价技术研究工作
关键词:
Keywords:
本文引用格式
王利鹏, 闫放, 李孜军, 王方.
WANG Lipeng, YAN Fang, LI Zijun, WANG Fang.
矿产资源是国民经济建设的重要物质基础。中国的高海拔地区,特别是西部地区,由于其独特的地理位置而拥有大量的矿产资源。随着西部大开发的实施和“一带一路”倡议的提出,近年来我国西部地区有许多矿山已投入生产或正在规划建设[1]。在矿产资源的开采过程中,矿井风机通过提供新鲜空气和消除有毒有害气体,为矿产资源安全开采提供保障,是高原地区矿井安全生产的重要组成部分。由于高海拔地区具有低压缺氧、气候干燥和空气稀薄等特点,给高原地区矿产资源的安全有效开采带来了挑战[2-3]。特别是低压缺氧的环境,矿井风机的工作效率受到了严重的影响,高原地区矿井风机的降效问题威胁着矿产资源的安全开采[4-6]。因此,建立综合的评估指标体系,对高海拔地区的矿井风机运行状态进行科学的评估,不断提高矿井风机的工作效率,不仅是矿井风机安全运行的必然要求,也是矿产资源安全、有效开采的重要保证。
目前,国内外学者对矿井风机的研究取得了比较好的成果。马宁等[7]从人员因素、设备因素和环境因素构建评价指标体系,利用模糊层次分析法对高原矿井人机功效进行评估。宋品芳等[8]构建了矿井风机性能评价指标体系,然后利用熵权法确定权重,最后利用熵权模糊综合评价法对高海拔矿井风机性能进行了综合评估。李琦等[9]利用柴油机机械效率数学模型分析了高海拔地区的环境因素对施工机械效率的影响程度和规律。王瑜敏等[4]根据矿井通风学原理,对高原地区矿机风机通风降效特征进行了研究,认为低压缺氧是高原地区矿井作业人员和设备降效的根本原因。Villiers等[10]利用实地调研的方法对南非金矿辅扇进行了综合的性能评估。Hirano等[11]利用CFD对小轴流风机进行了性能评估。
上述评估方法在高原地区矿井风机的安全管理工作中有很大的应用空间,但是也存在一定的局限性。例如,所选用的方法只能对矿井风机进行静态评估,无法对高原地区矿井风机的运行状态进行动态预测评估,同时无法克服评估过程中的不确定性。因此,如何构建不仅能有效克服评估过程中产生的不确定性而且能够对矿井风机运行状态进行动态预测评估的数学评估模型,是高原地区矿井风机安全管理面临的重要挑战。基于此,本文在已有的研究基础之上,构建了高原地区矿井风机评估指标体系,引入集对分析理论对高原地区矿井风机运行状态进行了科学评估。该方法不仅可以有效地克服评估过程中产生的不确定性,而且可以从静态和动态2个方面对高原地区矿井风机进行综合评估。结合静态和动态的矿井风机状态评估结果,可以为高原地区矿井风机的安全管理提供有针对性的管理方案,不断提高高原地区矿井风机的工作效率。
1 集对理论
1.1 集对分析
式中:
1.2 集对势
集对势可以反映2个集合在特定背景下的联系程度,可表示为
式中:当
1.3 偏联系数
五元联系数的偏联系数如
一阶偏联系数:
式中:
二阶偏联系数:
式中:
三阶偏联系数:
式中:
四阶偏联系数:
式中:
计算一阶偏联系数时,令
2 基于SPA-IAHP的矿井风机状态评估模型构建
2.1 构建评估指标体系
图1
(1)环境因素:高海拔地区具有低压缺氧的特点,它是造成高海拔地区设备降效的根本原因[4]。因此本文将环境温度、环境湿度和大气压力作为评估高海拔地区矿井风机的评估指标。
(2)管理因素:良好的管理对高海拔地区矿井风机的正常运行至关重要,它决定着矿井风机运行的状况。完善的管理制度、精准的应急方案和合理的人员分配将优化风机的运行状况。
(3)风机因素:风机因素是决定风机运行状况的决定性因素。根据《金属非金属矿山在用主通风机系统安全检查规范》(AQ 2054-2016)[23],将外部结构、保护设施、监控系统、维修措施、能耗状态和通风效率作为考核风机运行状态的决定性因素。
2.2 确定指标权重
高原地区矿井风机所处环境具有复杂性和不确定性,导致矿井风机指标权重的分配具有很大的不确定性。因此,引入不确定性层次分析法进行权重的计算。由于不确定性层次分析法得到的是权重区间,具有确定性和不确定性的特征,因此可以利用集对分析中的确定度、差异度和对立度进行不确定性分析,进而可以将区间权重转换成精确权重。具体的步骤如下:
(1)建立区间比较矩阵
基于1~9标度法,各位专家对评价指标两两比较进行打分,形成区间比较矩阵A,如
基于区间数理论[25],AL和AU可以获得,如
(2)计算最大特征值和最大特征向量
分别计算矩阵AL和AU对应的最大特征值,即
(3)计算区间权重
区间权重由式(
(4)计算精确权重
区间权重可以由
2.3 构建五元联系数同异反评估模型
构建“优秀,良好,中等,差,较差”状态评语集,邀请专家根据实际情况对指标集
可简化为五元联系数:
式中:
3 案例研究
以云南省某金矿为实例,该矿山位于海拔高度为3 400 m的云南省迪庆藏族自治州。该金矿为地下开采,开采深度为标高4 000 m左右,年产量约30万t。
3.1 矿井风机状态评估结果
根据建立的矿井风机运行状态评价指标,邀请专家利用1~9比例标度法对环境温度A1、环境湿度A2和大气压力A3进行评估,建立区间判断矩阵:
通过
通过MATLAB计算得到
通过式(
因此,利用式(15)~
通过同样的方式,可以计算出矿井风机所有评估指标的权重,如表1所示。
表1 矿井风机指标权重
Table 1
二级指标 | 权重 | 三级指标 | 权重 |
---|---|---|---|
环境因素A | 0.4989 | 环境温度A1 | 0.1207 |
环境湿度A2 | 0.2282 | ||
大气压力A3 | 0.6511 | ||
管理因素B | 0.0148 | 管理制度B1 | 0.1447 |
应急方案B2 | 0.2911 | ||
人员分配B3 | 0.5642 | ||
风机因素C | 0.4677 | 外部结构C1 | 0.4401 |
保护设施C2 | 0.0452 | ||
监控系统C3 | 0.0579 | ||
维修措施C4 | 0.1073 | ||
能耗状态C5 | 0.1632 | ||
通风效率C6 | 0.1863 | ||
人的因素D | 0.0186 | 心理状况D1 | 0.4152 |
身体状况D2 | 0.2164 | ||
工作能力D3 | 0.2825 | ||
教育水平D4 | 0.0859 |
邀请专家根据所建立的指标体系进行评估。以环境温度A1、环境湿度A2和大气压力A3为例,可得到状态评估矩阵
然后,利用3.2所计算的指标权重,根据
表2 矿井风机评估结果
Table 2
二级指标 | 三级指标 | 静态联系数 | 集对势 |
---|---|---|---|
环境因素A | 环境温度A1 | 0.01+0.30i1+0.40i2+0.10[i3+0.1j | 均势 |
环境湿度A2 | 0.00+0.30i1+0.40i2+0.20i3+0.10j | 均势 | |
大气压力A3 | 0.00+0.30i1+0.40i2+0.20i3+0.10j | 均势 | |
0.01+0.30i1+0.40i2+0.19i3+0.10j | 均势 | ||
管理因素B | 管理制度B1 | 0.30+0.20i1+0.30i2+0.10i3+0.10j | 偏同势 |
应急方案B2 | 0.20+0.40i1+0.30i2+0.00i3+0.10j | 均势 | |
人员分配B3 | 0.10+0.30i1+0.30i2+0.30i3+0.00j | 均势 | |
0.16+0.31i1+0.30i2+0.18i3+0.04j | 偏同势 | ||
风机因素C | 外部结构C1 | 0.10+0.20i1+0.30i2+0.10i3+0.30j | 偏反势 |
保护设施C2 | 0.00+0.30i1+0.30i2+0.30i3+0.10j | 均势 | |
监控系统C3 | 0.30+0.30i1+0.10i2+0.10i3+0.20j | 均势 | |
维修措施C4 | 0.10+0.30i1+0.30i2+0.20i3+0.10j | 均势 | |
能耗状态C5 | 0.10+0.30i1+0.50i2+0.00i3+0.10j | 均势 | |
通风效率C6 | 0.10+0.20i1+0.40i2+0.20i3+0.10j | 均势 | |
0.11+0.24i1+0.34i2+0.12i3+0.19j | 偏同势 | ||
人的因素D | 心理状况D1 | 0.10+0.30i1+0.40i2+0.00i3+0.20j | 均势 |
身体状况D2 | 0.00+0.10i1+0.30i2+0.40i3+0.20j | 偏反势 | |
工作能力D3 | 0.20+0.10i1+0.40i2+0.10i3+0.20j | 均势 | |
教育水平D4 | 0.20+0.40i1+0.20i2+0.10i3+0.10j | 均势 | |
0.11+0.21i1+0.38i2+0.12i3+0.18j | 均势 | ||
0.06+0.27i1+0.37i2+0.16i3+0.14j | 均势 |
表3 偏联系数计算结果
Table 3
一级指标 | 二级指标 | 一阶偏联系数 | 趋势 | 二阶偏联系数 | 趋势 | 三阶偏联系数 | 趋势 | 四阶偏 联系数 | 趋势 |
---|---|---|---|---|---|---|---|---|---|
环境因素 A | 环境温度A1 | 0.25+0.43i1+0.80i2+0.50i3 | 下降 | 0.37+0.35i1+0.62i2 | 下降 | 0.51+0.36i1 | 上升 | 0.59 | 上升 |
环境湿度A2 | 0.00+0.43i1+0.67i2+0.67i3 | 下降 | 0.00+0.39i1+0.50i2 | 下降 | 0.00+0.44i1 | 下降 | 0.00 | 过渡 | |
大气压力A3 | 0.00+0.43i1+0.67i2+0.67i3 | 下降 | 0.00+0.39i1+0.50i2 | 下降 | 0.00+0.44i1 | 下降 | 0.00 | 过渡 | |
0.03+0.43i1+0.68i2+0.66i3 | 下降 | 0.07+0.39i1+0.51i2 | 下降 | 0.63+0.37i1 | 下降 | 0.26 | 上升 | ||
管理因素 B | 管理制度B1 | 0.60+0.40i1+0.75i2+0.50i3 | 上升 | 0.60+0.35i1+0.60i2 | 下降 | 0.63+0.37i1 | 上升 | 0.63 | 上升 |
应急方案B2 | 0.33+0.57i1+1.00i2+0.00i3 | 上升 | 0.37+0.36i1+1.00i2 | 下降 | 0.50+0.27i1 | 上升 | 0.65 | 上升 | |
人员分配B3 | 0.25+0.50i1+0.50i2+1.00i3 | 下降 | 0.33+0.50i1+0.33i2 | 下降 | 0.40+0.60i1 | 下降 | 0.40 | 上升 | |
0.34+0.51i1+0.63i2+0.82i3 | 下降 | 0.40+0.45i1+0.43i2 | 下降 | 0.47+0.51i1 | 下降 | 0.48 | 上升 | ||
风机因素 C | 外部结构C1 | 0.33+0.40i1+0.75i2+0.25i3 | 上升 | 0.45+0.35i1+0.75i2 | 下降 | 0.57+0.32i1 | 上升 | 0.64 | 上升 |
保护设施C2 | 0.00+0.50i1+0.50i2+0.75i3 | 下降 | 0.00+0.50i1+0.40i2 | 下降 | 0.00+0.56i1 | 下降 | 0.00 | 过渡 | |
监控系统C3 | 0.50+0.75i1+0.50i2+0.33i3 | 上升 | 0.40+0.60i1+0.60i2 | 下降 | 0.40+0.50i1 | 下降 | 0.44 | 上升 | |
维修措施C4 | 0.25+0.50i1+0.60i2+0.67i3 | 下降 | 0.33+0.45i1+0.47i2 | 下降 | 0.42+0.49i1 | 下降 | 0.46 | 上升 | |
能耗状态C5 | 0.25+0.38i1+1.00i2+0.00i3 | 上升 | 0.40+0.27i1+1.00i2 | 下降 | 0.59+0.21i1 | 上升 | 0.74 | 过渡 | |
通风效率C6 | 0.33+0.33i1+0.67i2+0.67i3 | 下降 | 0.50+0.33i1+0.50i2 | 下降 | 0.60+0.40i1 | 上升 | 0.60 | 上升 | |
0.31+0.41i1+0.74i2+0.39i3 | 下降 | 0.43+0.36i1+0.66i2 | 下降 | 0.55+0.35i1 | 上升 | 0.61 | 上升 | ||
人的因素 D | 心理状况D1 | 0.25+0.43i1+1.00i2+0.00i3 | 上升 | 0.37+0.30i1+1.00i2 | 下降 | 0.55+0.23i1 | 上升 | 0.70 | 上升 |
身体状况D2 | 0.00+0.25i1+0.43i2+0.67i3 | 下降 | 0.00+0.37i1+0.39i2 | 下降 | 0.00+0.48i1 | 下降 | 0.00 | 过渡 | |
工作能力D3 | 0.67+0.20i1+0.80i2+0.33i3 | 上升 | 0.77+0.20i1+0.71i2 | 下降 | 0.79+0.22i1 | 上升 | 0.78 | 上升 | |
教育水平D4 | 0.38+0.67i2+0.67i2+0.50i3 | 下降 | 0.33+0.50i1+0.57i2 | 下降 | 0.40+0.47i1 | 下降 | 0.46 | 上升 | |
0.34+0.36i1+0.76i2+0.40i3 | 下降 | 0.49+0.32i1+0.66i2 | 下降 | 0.61+0.33i1 | 上升 | 0.65 | 上升 | ||
0.18+0.42i1+0.7i2+0.53i3 | 下降 | 0.30+0.38i1+0.57i2 | 下降 | 0.44+0.40i1 | 上升 | 0.53 | 上升 |
从二级指标的情况来看,环境因素A、管理因素B、风机因素C和人的因素D的状态评估结果均处于中等水平,环境因素A和人的因素B的集对势为均势,管理因素C和风机因素D的集对势为偏同势。说明在矿井风机的安全管理中应该重点关注环境因素和人的因素。
从三级指标的评估结果来看,13个指标的集对势为均势,2个指标的集对势为偏反势,1个指标的集对势为偏同势。对于处于偏反势的2个指标,即外部结构C1和身体状况D2,应该采取相应的技术措施和管理措施以提高该指标目前的状态。
3.2 矿井风机状态发展趋势分析
偏联系数可以反映矿井风机状态动态发展趋势。根据表2的计算结果,风机状态总的评估结果为中等水平。尽管矿井风机的运行状态不是差或者较差的状态,但是评估结果的第一阶偏联系数和第二阶偏联系数呈现下降趋势,表明风机的运行状态存在向较差状态发展的趋势。因此,矿井风机的管理工作不能有任何的懈怠,需要投入更多的精力去改善风机的运行状态。
对于二级指标,环境因素A和管理因素B的第一阶偏联系数、第二阶偏联系数和第三阶偏联系数均处于下降状态,说明这些指标在未来的一段时间内向较差状态方向发展。因此,应该加强对环境的检测和管理,使矿井风机处于一个良好的运行环境,减少高原地区环境因素对矿井风机正常运行的影响。对于风机因素C和人的因素D,2个指标的第一阶偏联系数和第二阶偏联系数均为下降状态;此外,2个指标的第三阶偏联系数和第四阶偏联系数为上升状态。这些信息表明,人的因素C和风机因素D的状态在未来的一段时间里将呈现动态的变化趋势。因此,通过偏联系数可以从宏观上预测风机未来管理工作中的重点,以达到不断提高风机工作效率的目的。
为了能够更有针对性地对矿井风机的运行状态进行预测,笔者从微观方面对每个评估指标进行了详细分析。对于三级指标,如图2所示,这些因素(管理规定B1,应急管理程序B2,保护装置C3,维护措施C5,人员的心理状态D1和人员的知识水平D3)的一阶偏联系数均大于0,说明这些指标的状态在未来的一段时间内存在向更好状态发展的趋势;其中,人员的知识水平D3的值最大,说明该指标具有向更好状态发展的稳定趋势。从图2可以看出,这些指标(环境温度A1,环境湿度A2,大气压力A3,人员配备合理性B3,风扇的运行状态C1,外部结构状态C2,监视系统C4,能耗状态C6和人员身体状况D2)的第一阶偏联系数均小于0,说明这些指标存在向较差状态发展的趋势;其中,人员配备合理性B3和人员身体状况C2的状态呈现急剧向较差状态发展的趋势。因此,在矿井风机的管理工作中应该优先关注状态处于下降或者急剧下降的指标。
图2
基于图3,所有指标的第二阶偏联系数均小于0,因此,在未来的一段时间内将会持续呈现向较差状态发展的趋势。
图3
图4
图5
另外,从整体上分析,由于这些因素(环境湿度A2,外部结构状态C2和人员的身体状况D2)的第一偏联系数、第二偏联系数和第三偏联系数均小于0,表明这些指标的状态将持续呈现下降趋势,意味着这些指标持续向较差状态方向发展。此外,这4个因素的第四阶偏联系数全部处于过渡状态,表明这些指标的状况向较差状态方向发展的趋势较强。因此,这些因素应该列为高风险因素,作为重点管理对象。例如,采用无损检测技术检测矿井风机的外部结构;制定合理的休息计划以改善人员身体状况等。综上所述,这些动态的评估结果可以为高原矿山风机的管理提供重要思路,管理者可以从这些因素着手进行矿山风机的优化管理,以达到改善高原地区矿山风机状态的目的。
针对高原地区矿井风机状态评估方法中存在的局限性,引入了集对分析理论。从评估结果来看,该方法能够有效地克服评估过程中的不确定性,使评估结果更加准确。此外,利用偏联系数理论可以对矿井风机的运行状态进行准确预测,掌握矿井风机状态未来的发展趋势,为矿井风机管理者提供有针对性的管理方案。因此,本研究所提出的评估方法可以有效弥补现存评估方法的不足,评估结果可以为矿井风机的安全管理提供指导性建议。
4 结论
(1) 结合不确定性区间层次分析法和三元联系数,克服了无法得到矿井风机指标精确权重的问题,使获得的指标权重更加准确。
(2) 解决了高原矿山风机评估过程中的不确定问题。构建的五元联同异反评估模型能够更有效地处理不确定问题,使评估结果更加准确和科学,矿井风机的评估结果为中等水平。通过对比云南省迪庆藏族自治州高原地区矿井风机的实际情况,验证了该方法的有效性。
(3) 实现了矿井风机状态动态评估。引入集对分析理论中的偏联系数,综合评估了矿井风机每个指标未来的发展趋势,实现了矿井风机动态状态评估;结合静态评估,构建了静态和动态相结合的评估方法,使评估结果具有系统性和准确性。
http://www.goldsci.ac.cn/article/2020/1005-2518/1005-2518-2020-28-6-930.shtml
参考文献
高寒高海拔地区矿山工程设计要点思考
[J].
Key design consideration for the mine engineering in the alpine and high-altitude regions
[J].
高海拔矿山通风系统改造方案优选研究
[J]
Optimization of ventilation system transformation at high altitude mine
[J].
基于Ventsim的高原矿井通风系统优化
[J].
Optimization of plateau mine ventilation system based on Ventsim
[J].
高海拔矿井风机通风降效特征的研究
[J].
Research on characteristics of ventilation effect’s reduction of fan in high altitude mines
[J].
高原寒冷条件下南水北调西线工程施工人机效率的思考
[J].
On manpower and machine efficiency of water diversion from south to north for west line construction under the condition on cold plateau
[J].
高寒地区人工机械施工效率研究
[J].
Research on construction efficiency of artificial machinery in alpine region
[J].
基于模糊层次分析法的高原矿井人机功效评价
[J].
Evaluation of human-machine effectiveness of plateau mine based on fuzzy analytic hierarchy process
[J].
基于熵权模糊法的高海拔矿井风机性能影响因素分析
[J].
Analysis of factors influencing the performance of high altitude mine fan based on entropy weight fuzzy method
[J].
高海拔环境对施工设备机械效率的影响研究
[J].
Study on mechanical efficiency of construction equipment in high altitude environment
[J].
Evaluating the impact of auxiliary fan practices on localised subsurface ventilation
[J].
Study on performance evaluation of small axial fan
[J].
集对论——一种新的不确定性理论方法与应用
[J].
A new method and application of uncertainty theory
[J].
A set pair analysis based layer of protection analysis and its application in quantitative risk assessment
[J].
Application of a cloud model-set pair analysis in hazard assessment for biomass gasification stations
[J].
基于广义集对分析的系统危险性评价研究
[D].
Research on System Risk Assessment Based on Generalized Set Pair Analysis
[D].
Water resources carrying capacity evaluation and diagnosis based on set pair analysis and improved the entropy weight method
[J].
多元联系数法在高速公路PPP项目风险评价中的应用
[J].
Application of multivariate connection number method in risk assessment of expressway PPP project
[J].
基于集对分析疗效曲线与偏联系数的银屑病用药优选探讨
[J].
Optimization of psoriasis medication based on set pair analysis,efficacy curve and partial correlation coefficient
[J].
基于多元联系数的矿工不安全行为风险态势评估
[J].
Risk state assessment on unsafe behavior of miners based on multivariate connection number
[J].
基于联系数的地铁隧道施工邻近建筑物风险评价
[J].
Risk assessment of neighboring building in metro tunneling construction based on connection number
[J].
偏联系数在投标企业发展趋势分析中的应用
[J].
Apply of partial connection number in the analysis development trend of bidding enterprise
[J].
基于偏联系数的医院医疗质量发展趋势综合分析
[J].
The synthetic analysis of developing trend of medical quality of hospital based on partial connection number
[J].
高原矿井作业人员环境适应性分析
[J].
Environment adaptability analysis of high altitude miners
[J].
Comprehensive benefit evaluation of the power distribution network planning project based on improved IAHP and multi-level extension assessment method
[J].
/
〈 |
|
〉 |
