基于拓扑网络算法的煤矿瓦斯爆炸风险度量
2022-05-06成连华郭阿娟刘黎郭慧敏
成连华 郭阿娟 刘黎 郭慧敏
摘要：為有效预防煤矿瓦斯爆炸事故，对瓦斯爆炸风险进行量化分析。首先，通过分析2010—2020年152起瓦斯爆炸事故调查报告，基于层次全息建模方法识别出67个风险因素，提取出214条事故链，构建瓦斯爆炸风险演化路径;然后将迭代加权思想与设计结构矩阵相结合，提出改进集成风险传播法，计算风险演化过程中风险的直接和间接传播概率;在此基础上，将风险因素作为网络节点，因素之间的影响关系作为网络边，构建瓦斯爆炸风险拓扑网络度量模型，并运用Gephi 0.9.2软件分析计算网络拓扑参数，实现风险演化过程中的累积风险可视化。结果表明：煤矿瓦斯爆炸风险拓扑网络呈现明显小世界网络特征，聚集程度偏低而演化性较强;识别出最短风险演化路径7条，关键性风险节点包括施工操作不当、风量不足、瓦斯异常涌出;初始风险因素不同，在风险演化过程中路径不同，风险权值也存在差异，且在风险演化过程中存在风险损失现象。研究结果为煤矿瓦斯爆炸风险评估提供理论依据，对预防瓦斯爆炸事故发生具有现实意义。
关键词：瓦斯爆炸;风险演化;拓扑网络;风险度量;网络参数
中图分类号：X 936文献标志码：A
文章编号：1672-9315（2022）02-0268-08
DOI：10.13800/j.cnki.xakjdxxb.2022.0210开放科学（资源服务）标识码（OSID）：
Risk measurement of coal mine gas explosion based
on topological network algorithmCHENG Lianhua GUO Ajuan LIU Li GUO Huimin
（1.College of Safety Science and Engineering，Xian University of Science and Technology，Xian 710054，China;
2.West Mine Gas Intelligent Drainage Engineering Research Center，
Xian University of Science and Technology，Xian 710054，China）Abstract：In order to effectively prevent gas explosion accidents，the risk of gas explosion is quantitatively analyzed.Firstly，by analyzing the investigation reports of 152 gas explosion accidents from 2010 to 2020，67 risk factors were identified based on the hierarchical holographic model method，214 accident chains were extracted，and the gas explosion risk evolution path was constructed.Then the iterative weighting idea was combined with the design structure matrix with an improved integrated risk propagation method proposed to calculate the direct and indirect propagation probability of risk in the process of risk evolution.The risk factors，accordingly，as network nodes and the influence relationships between factors as network edges to build a topological network measurement model of gas explosion risk.The Gephi 0.9.2 software was used to analyze and calculate the network topological parameters，so as to realize the visualization of cumulative risks in the risk evolution process.The results show that the topological network of coal mine gas explosion risk presents small-world network characteristics，with low aggregation degree and strong evolution.Seven shortest risk evolution paths are identified，and the key risk nodes involve improper construction operation，insufficient air volume and abnormal gas emission.With the initial risk factors and the paths of risk evolution being different，the risk weights are also different.There arises risk loss in the risk evolution process.The research results provide a theoretical basis for the risk assessment of coal mine gas explosion，and have practical significance for the prevention of gas explosion accidents.
Key words：gas explosion;risk evolution;topological network;risk measurement;network parameter
0引言
近年来，中国煤炭行业安全生产状况持续好转，但随着煤炭开采向深部转移，煤层瓦斯含量急剧升高，瓦斯爆炸风险依然存在。国家矿山安全监察局统计数据显示，2010—2020年煤矿事故死亡10 203人，因瓦斯爆炸事故死亡2 762人，约占煤矿事故死亡人数的27.10%，瓦斯爆炸风险严重威胁工人的生命安全。因此，对瓦斯爆炸风险演化过程定量分析，对制定预防措施具有重要现实意义。
目前国内外学者围绕煤矿瓦斯爆炸风险评价展开了大量研究。LI等基于改进层次分析法、贝叶斯网络等方法，评估了煤矿瓦斯爆炸风险[1-2];TONG等运用概率风险评估模型、蒙特卡罗法分析煤矿瓦斯爆炸事故中矿工不安全行为所造成的事故风险[3-4];田水承等运用文献分析法和熵权法确定矿工不安全状态评价指标综合权重，并进行实例分析[5];王婉青和李红霞分别基于IAHP方法、模糊综合评价法分析了瓦斯爆炸事故风险中的关键风险，同时对事故后果提出预防措施[6-7]。但大多数风险评价方法集中于对事故结果进行分析，风险演化过程的定量分析有待进一步研究。拓扑网络是一种用来分析复杂系统相互作用及演化的工具，由于其契合复杂系统的特点，在城市交通[8-9]、电力系统[10]、房建事故[11]等诸多领域得到应用。如ZHOU等构建地铁施工事故网络（SCAN）探究26种事故之间的联系[12];丁明等运用动态故障树理论，提出停电事故链网络模型，并评估支路的关键风险程度[13];周新宇等构建房建事故行为致险网络模型，解析了人的不安全行为与事故间的作用机理[14]。然而当前拓扑网络应用于煤矿瓦斯爆炸风险分析研究主要集中于无向无权网络，对于有向加权网络研究还处于探索阶段。
鉴于此，运用拓扑网络，将影响瓦斯爆炸事故的风险因素作为节点，因素之间的演化关系作为连接边，通过计算风险演化过程中直接和间接传播概率，构建煤矿瓦斯爆炸风险拓扑网络度量模型。同时对节点出入度、聚类系数及最短路径等拓扑参数进行分析，确定事故演化最短路径和关键风险因素，以期为煤矿瓦斯爆炸风险量化分析提供参考。
1瓦斯爆炸风险演化分析
1.1瓦斯爆炸事故案例收集
以国家矿山安全监察局官网中统计的煤矿瓦斯爆炸事故调查报告作为数据来源，通过搜集2010—2020年的瓦斯爆炸事故调查报告，梳理了152起瓦斯爆炸事故案例，主要包含事故基本情况、事故经过、事故原因等基本信息。
1.2瓦斯爆炸风险因素分析
基于层次全息建模（hierarchical holographic modeling，HHM）方法[15]从多维度、多层次进行事故过程识别。结合安全系统工程理论从人、物、管理、环境、技术方法5个维度进行分类，归纳分析出67个风险因素、2个风险事件。
统计瓦斯爆炸风险出现频次，按公式（1）计算发生概率Pi。瓦斯爆炸事故Ck中如果出现风险因素类别i，则Ck=1，否则Ck=0。其中k为某个瓦斯爆炸事故案例;k={1，2，…，n}，n=152。
由于在一次瓦斯爆炸事故中可能存在多个风险，因此概率和大于1。统计分析瓦斯爆炸风险结果见表1。
1.3瓦斯爆炸风险演化路径构建
任何一起事故过程都可用一个或几个事故链描述。在瓦斯爆炸风险演化过程中，当风险节点状态偏离时，风险累积达到风险阀值，触发整个系统发生连锁反应，进而导致事故发生。
依照事故发展流程，确定风险节点的时序关系和逻辑顺序，将风险节点依序相连形成事故链[16-17]，实现对瓦斯爆炸风险演化过程的描述。以山西平遥二亩沟煤业瓦斯爆炸事故为例进行分析，对事故直接原因①和②进行分析，可先后提取出D17与B22的2个风险因素。事故经过表明：炮采组工人张某（无爆破工特种作业人员证件）未执行“一炮三检”和“三人连锁爆破”制度违章爆破，爆破明火引爆9103工作面采空区涌入瓦斯，发生瓦斯爆炸，可先后提取出A16，A17和D21的3个风险因素和F11这1个风险事件。最终从该案例中提取出8个风险因素、1个风险事件和4条事故链，见表2。
分析152起煤矿瓦斯爆炸事故案例，將含有共同风险因素的事故链融合到一个全局网络中[18]，运用系统动力学软件构建煤矿瓦斯爆炸风险演化路径，如图1所示。
2.1风险网络连边计算
引入交互作用矩阵（factor interaction matrix，FIM）理论计算风险连边权重，以分析系统影响因素的交互作用关系，如图2所示。
因子交互矩阵共有n个风险因素，主对角上各风险因素表示为Pi（i=1，2，…，n），非对角上方区域表示对角线上其他风险因素对相应Pi的影响xij，非对角下方区域表示对角线上Pi对其他各风险因素的影响xji。每一行对应的数值相加表示为Ci，每一列的对应数值相加表示为Ei，则风险网络连边权重wi的计算表达式为
由于文章篇幅所限，以安全检查不到位→线头外露→电路老化→电气火花→引爆瓦斯（D19→B14→B11→B22→F11）风险演化路径为例进行计算，由FIM理论计算可得风险网络连边权重wi=（0.55，0.48，0.35，1.08）。
2.2风险网络迭代加权
设计结构矩阵（design structure matrix，DSM）是一种用来表示系统中元素及相互作用的结构化建模方法，反映风险传播网络中一个元素通过所有可能路径影响另一个元素的集成风险传播概率，可直观对复杂系统进行可视化分析[19-20]。笔者将迭代加权思想与DSM相结合，构建集成风险传播DSM矩阵。
2.3風险拓扑网络形成
将风险迭代加和权重转化成CSV文件，并运用Gephi 0.9.2软件构建煤矿瓦斯爆炸风险加权拓扑网络，如图4所示。
从图4可以看出，煤矿瓦斯爆炸风险拓扑网络模型是一个完整的连通网络。圆圈代表67个风险节点，连线代表214条网络连边，连线的粗细程度代表节点间的集成传播概率。如B25线条最粗，可直观看出与该风险链接路径最多，风险权值变化最大。通过定量描述瓦斯爆炸风险加权网络中风险因素之间的复杂关系，为煤矿瓦斯爆炸风险度量提供了一个新的视角。
3拓扑网络风险度量分析
3.1风险节点度分析
节点度分析是网络拓扑结构最重要的指标之一，反映了节点在网络中的直接影响力和重要性[22]。节点的度值越大，表明节点越重要。节点的度可表示为
通过Gephi 0.9.2软件分析得到节点度分布如图5所示，横坐标为度贡献值，纵坐标为度分布数量。平均度为2.315，表明每个风险因素至少与2个不安全因素存在链接关系。
选取风险节点中出、入度数排前的节点共20个，度分布如图6所示。从图6可以看出，B20，B34的入度和出度最高，表明诱发该风险的途径较多，因此对该风险进行有效控制，抑制瓦斯爆炸发生的效果最佳。同理，A12，B17，B24等度分布较高的风险节点也是风险防控中需重点关注对象。
3.2风险聚类分析
聚类系数是网络节点中邻接节点之间相互连接的比例，反映了网络节点集团化程度。对于与k个节点相连的节点i而言，它最多可与这些邻居节点构成k（k-1）/2条连边，则实际存在连边Ei与最多可能连边个数之比是节点i的聚类系数Ci。
通过Gephi 0.9.2软件分析得到风险网络模型的聚类系数为0.212，表明在整个网络演化中，导致瓦斯爆炸的成因复杂，但大部分风险因素之间的相互影响并不明显。
筛选节点加权聚类系数排前的节点共20个，如图7所示。从图7可以看出，B20，B33聚类系数最高，表明该节点集团化严重，若对该节点及邻接节点构成的子系统采取断链控制措施，可有效减缓事故发生。
3.3最短路径分析
节点度分布及聚类系数分别分析了风险节点的重要程度及聚类程度，但未考虑上一节点引发后续节点的风险传递。通过Gephi 0.9.2软件分析得到平均最短路径为2.506，表明瓦斯爆炸风险拓扑网络中一个节点状态发生变化平均仅通过2506步就能够引起邻接其它节点发生变化。
Dijkstra算法是荷兰著名计算机科学家E.W.Dijkstra于1959年提出关于图论中求最短路径常用算法[23]。借鉴Dijkstra算法，编制计算程序，求解瓦斯爆炸风险最短路径，得到风险节点的最短路径及其最短路径值（风险权值），见表3。
通过分析表3可知，风险演化路径越短，权重值越大，但在风险路径传递过程中风险权值逐渐减少，说明该过程存在风险损失;B31→F12风险路径最短，风险权重变化显著，表明B31（冲击波蔓延）对瓦斯爆炸风险演化影响较大;虽然B11→B16→B25、A20→B20→F11和B17→B16→B25这3条最短演化路径都是由3个风险因素构成，但初始风险因素不同，风险演化路径不同，所以风险权值也不同。此外，风险演化最短路径也从侧面反映了瓦斯爆炸事故发生的可能性。在实际瓦斯爆炸风险管控工作中，优先对最短风险路径上的风险因素加以控制，可以提高解决事故的效率。
4结论
1）通过统计分析2010—2020年共计152起瓦斯爆炸事故案例，采用HHM方法从人、物、管理、环境、技术方法5个维度进行风险分析，得到67个风险因素，形成214条瓦斯爆炸风险演化路径链。
2）以拓扑网络算法为基础，引入DSM对瓦斯爆炸风险进行迭代加权计算，并运用Gephi 0.9.2软件构建瓦斯爆炸风险拓扑网络度量模型，实现瓦斯爆炸风险演化过程的量化分析，以揭示事故发生的潜在特征。
3）通过分析瓦斯爆炸风险有向带权网络参数，发现风险节点聚集程度较低，呈现明显小世界网络特征，表明在整个风险演化过程中大部分风险因素之间相互影响不明显，存在显著传递关系。
4）运用Dijkstra算法搜寻风险演化最短路径，得到7条最短风险演化路径和3个关键风险因素，优先对瓦斯爆炸最短风险路径上关键风险加以管控，对预防瓦斯爆炸事故发生、保障企业安全运行具有现实意义。
参考文献（References）：
[1]LI M，WANG H T，WANG D M，et al.Risk assessment of gas explosion in coal mines based on fuzzy AHP and Bayesian Network[J].Process Safety and Environmental Protection，2020，135：207-218.
[2]SHI L，WANG J H，ZHANG G M，et al.A risk assessment method to quantitatively investigate the methane explosion in underground coal mine[J].Process Safety & Environmental Protection，2017，107（8）：317-333.
[3]TONG R，YANG Y G，MA X F，et al.Risk assessment of miners unsafe behaviors：A case study of gas explosion accidents in Coal Mine China[J].International Journal of Environmental Research and Public Health，2019，16（10）：1765.
[4]许素睿.基于蒙特卡罗法的煤矿工人不安全行为风险评估[J].中国安全科学学报，2020，30（4）：172-178.XU Surui.Risk assessment of unsafe behaviors of coal miners based on Monte Carlo method[J].China Safety Science Journal，2020，30（4）：172-178.
[5]田水承，孙璐瑶，唐艺璇，等.矿工不安全状态评价指标体系的构建与分析[J].西安科技大学学报，2021，41（3）：402-409.TIAN Shuicheng，SUN Luyao，TANG Yixuan，et al.Construction and analysis of evaluation index system for miners unsafe state[J].Journal of Xian University of Science and Technology，2021，41（3）：402-409.
[6]王婉青，吕淑然，张宇栋.基于IAHP方法的煤矿瓦斯爆炸事故分析[J].煤炭技术，2019，38（5）：124-126.WANG Wanqing，LYU Shuran，ZHANG Yudong.Analysis of coal mine gas explosion accident based on IAHP method[J].Coal Technology，2019，38（5）：124-126.
[7]李红霞，郑佳敏.煤矿瓦斯爆炸模糊Bow-tie模型分析[J].矿业安全与环保，2020，47（5）：119-126.LI Hongxia，ZHENG Jiamin.Analysis of fuzzy Bow-tie model for gas explosion in coal mine[J].Mining Safety & Environmental Protection，2020，47（5）：119-126.
[8]王红勇，温瑞英.基于复杂网络的空中交通态势风险评估方法[J].中国安全科学学报，2018，28（5）：172-178.WANG Hongyong，WEN Ruiying.Risk assessment method for air traffic situation based on complex network[J].China Safety Science Journal，2018，28（5）：172-178.
[9]胡立伟，范仔健，张苏航，等.基于复杂网络的城市交通拥塞因子风险传播机理及其应用研究[J].交通运输系统工程与信息，2021，21（2）：224-230.HU Liwei，FAN Zaijian，ZHANG Suhang，et al.Research on risk transmission mechanism and application of urban traffic congestion factor based on complex network[J].Transportation Systems Engineering and Information，2021，21（2）：224-230.
[10]梅丹，王公寶，叶志浩，等.舰船电力系统复杂网络特性研究[J].海军工程大学学报，2016，28（6）：10-14.MEI Dan，WANG Gongbao，YE Zhihao，et al.Research on complex network characteristics of naval power system[J].Journal of Naval University of Engineering，2016，28（6）：10-14.
[11]张子琛，郭聖煜，周新宇，等.工程不安全行为风险后果量化评价[J].土木工程与管理学报，2021，38（5）：125-129，139.ZHANG Zichen，GUO Shengyu，ZHOU Xinyu，et al.Quantitative evaluation of consequences of behavioral risk in construction[J].Journal of Civil Engineering and Management，2021，38（5）：125-129，139.
[12]ZHOU Z，IRIZARRY J，LI Q，et al.Using network theory to explore the complexity of subway construction accident network（scan）for promoting safety management[J].Safety Science，2014，64（64）：127-136.
[13]丁明，肖遥，张晶晶，等.基于事故链及动态故障树的电网连锁故障风险评估模型[J].中国电机工程学报，2015，35（1）：821-829.DING Ming，XIAO Yao，ZHANG Jingjing，et al.Risk assessment model of power grid cascading faults based on fault chain and dynamic fault tree[J].Proceedings of the CSEE，2015，35（1）：821-829.
[14]周新宇，郭聖煜.房建工程事故行为致险链作用机理研究[J].工程管理学报，2019，33（1）：112-116.ZHOU Xinyu，GUO Shengyu.Mechanism of behavioral risk chain for building construction accident[J].Journal of Engineering Management，2019，33（1）：112-116.
[15]刘家国，崔进，周欢，等.基于HHM-RFRM的船舶航行风险评估方法研究[J].中国管理科学，2019，27（5）：174-183.LIU Jiaguo，CUI Jin，ZHOU Huan，et al.Research on navigation risk assessment method based on HHM-RFRM[J].Chinese Journal Management Science，2019，27（5）：174-183.
[16]LAM C Y，FUSE M.Risk analysis for consumer-level utility gas and liquefied petroleum gas incidents using probabilistic network modeling：A case study of gas incidents in Japan[J].Reliability Engineering & System Safety，2019，185：198-212.
[17]LAM C Y，CRUZ A M.Assessment of risk factors and effects in hydrogen logistics incidents from a network modeling perspective[J].International Journal of Hydrogen Energy，2019，44：20572-20586.
[18]阳雨平，黄丕森，陈国国.基于改进FIM-未确知测度的尾矿库风险评价模型及应用[J].安全与环境学报，2021，21（3）：996-1004.YANG Yuping，HUANG Pisen，CHEN Guoguo.Safety evaluation of tailings pond based on FIM-optimization unascertained measure[J].Journal of Safety and Environment，2021，21（3）：996-1004.
[19]YANG Q，YAO T，ZHANG B.An overlapping-based design structure matrix for measuring interaction strength and clustering analysis in product development project[J].IEEE Transactions on Engineering Management，2014，61（1）：159-170.
[20]楊青，吕佳芮，索尼亚.基于设计结构矩阵（DSM）的复杂研发项目建模与优化研究进展[J].系统工程理论与实践.2016，36（4）：989-1002.YANG Qing，LYU Jiarui，Sonia.Progress in modeling and optimization for complex R&D based on design structure matrix[J].Systems Engineering Theory & Practice，2016，36（4）：989-1002.
[21]HAMRAZ B，CALAWEL N H，CLARKSON P J.A matrix-calculation-based algorithm for numerical change propagation analysis[J].IEEE Transactions on Engineering Management，2013，60（1）：186-198.
[22]薛锋，何传磊，黄倩.成都地铁网络的关键节点识别方法及性能分析[J].中国安全科学学报，2019，29（1）：93-99.XUE Feng，HE Chuanlei，HUANG Qian.Identification of key nodes in Chengdu metro network and analysis of network performance[J].China Safety Science Journal，2019，29（1）：93-99.
[23]王特起，谢亚琴.基于Dijkstra算法的校园导航系统的设计与实现[J].通信技术，2019，52（8）：1937-1943.WANG Teqi，XIE Yaqin.Design and implementation of campus navigation system based on Dijkstra algorithm[J].Communications Technology，2019，52（8）：1937-1943.
