基于突变级数法的煤矿外因火灾危险等级评价.pdf
doi10. 3969/ j. issn. 1001-8352. 2018. 06. 011 基于突变级数法的煤矿外因火灾危险等级评价 ❋ 贾进章①② 董铭鑫①② 胡竞文③ ①辽宁工程技术大学安全科学与工程学院辽宁阜新, 123000 ②矿山热动力灾害与防治教育部重点实验室辽宁阜新,123000 ③中共呼伦贝尔市委统战部内蒙古呼伦贝尔,021000 [摘 要] 为了给煤矿外因火灾危险等级评价提供一种更加合理的方法,采用将突变理论与模糊数学相结合的突 变级数法,从明火致因、电火致因、其他原因等 3 个方面分析了煤矿外因火灾的影响因素,建立了煤矿外因火灾危 险评价指标体系,并运用突变级数法对煤矿外因火灾危险等级进行评价。 以沈阳焦煤集团林盛煤矿为例,基于突 变级数法,求得该煤矿外因火灾危险等级评价的总突变隶属矩阵,依据总突变隶属函数值越大越优的原则,得到该 煤矿外因火灾危险等级为二级,即较安全。 依据评价结果,提出了提高该煤矿外因火灾安全性的相关对策。 [关键词] 危险等级评价;煤矿外因火灾;突变级数法;评分 [分类号] TD75 +2. 3 Risk Grade uation of External Fire in Coal Mine Based on Catastrophe Progression JIA Jinzhang①②, DONG Mingxin①②, HU Jingwen③ ① College of Safety Science and Engineering, Liaoning Technical University Liaoning Fuxin, 123000 ② Key Laboratory of Mine Thermodynamic Disasters and Control of Ministry of Education Liaoning Fuxin, 123000 ③ The United Front Department of the CPC Hulun Buir Municipal Party Committee Inner Mongolia Hulun Buir, 021000 [ABSTRACT] In order to provide a more reasonable to uate the risk rank of external fire in coal mine, an uation named catastrophe progression , which combines catastrophe theory and fuzzy math , was applied. Potential causes include coal mine external fires, open fires, electric fires and other reasons, were identified. Hazard uation index system of coal mine external fire was populated, and the catastrophe progression was applied to uate the risk ranks. Linsheng Coal Mine of Shenyang Coking Coal Group was taken as an example, and catastrophe progression was used to determine the total mutation membership matrix of the fire risk grade assessment of the coal mine. According to the principle that bigger total mutations subordinate function values tend to be better rated, the fire danger grade is assessed as grade two, which means relatively safe. Accordingly countermeasures to improve the external causes of fire risk in coal mines were put forward. [KEYWORDS] risk grade uation; external fire of coal mine; catastrophe progression ; scoring 引言 煤矿火灾危险性高,事故后果严重,而且易造成 一系列次生灾害。 一般将矿井中由外来热源引发可 燃物质燃烧造成的火灾称为煤矿外因火灾。 煤矿外 因火灾具有突发性和随机性强、传播速度快的特点, 最快速度可达 50 100 m/ h[1]。 煤矿外因火灾会造 成大量的人员伤亡、严重的财产损失以及恶劣的社 会影响。 据统计,9 成以上的煤矿重大火灾事故是 由外因火灾所引起,而且外因火灾事故致死的人数 占火灾总死亡人数的 65% [2]。 目前,国外方面,美国对火灾的研究较为深入, 美国国家标准技术委员会NIST对性能化防火和 第 47 卷 第 6 期 爆 破 器 材 Vol. 47 No. 6 2018 年 12 月 Explosive Materials Dec. 2018 ❋ 收稿日期2017-12-04 基金项目国家自然科学基金矿井瓦斯爆炸冲击波与通风动力耦合传播机理研究51374121 作者简介贾进章1974 - ,男,博士,教授,研究方向为矿井通风与灾害防治及安全评价。 E-mailjiajinzhang@163. com 通信作者董铭鑫1991 - ,男,硕士研究生,研究方向为矿井瓦斯爆炸、建筑火灾安全及安全评价。 E-mail1148404056@ qq. com 超高建筑火灾评价进行过专题研究[3];此外,火灾 评价模型的开发方面,众所周知的有 BFSEM 模 型[4]和 FRECAM 模型[5]等。 国内方面,主要是利用 层次分析法、模糊综合评价法以及神经网络评价法 等对煤矿火灾危险性进行评价。 王旭等[6]在分析 了煤矿事故特点后,将模糊综合评价法运用到了煤 矿安全评价中;邓军等[7]深入研究了层次分析法, 并将其应用到煤矿安全生产能力指标体系以及火灾 危险性评价工作中;冯海明等[8]将神经网络应用到 煤炭自燃预测工作中,并且在 Matlab 上得到了实 现。 上述方法虽然能够对煤矿火灾危险性作出合理 的评价,但是其计算过程复杂,有些指标的权重的计 算工作量大,而且人为主观因素影响大。 通过查阅 文献,虽然突变级数法作为比较成熟的评价方法被 广泛应用于工程评价、经济学和社会学等领域,但未 见将突变理论应用于煤矿外因火灾安全评价方面的 研究。 煤矿外因火灾具有独特的突变性,火灾是否发 生以及发生的时间都很难测算,即使计算出了火灾 发生概率,也难以确定它是否会造成煤矿事故风险。 将突变级数法[9]引入到煤矿外因火灾危险等级评 价中,更好地凸显了突变级数法处理突变事件的优 势。 该方法省去对各指标采取权重的工作,只需考 虑各指标对煤矿外因火灾影响的相对重要性[10],不 仅保证了评价结果的真实性、客观性,而且减轻了工 作量。 1 基于突变级数法的煤矿外因火灾危 险等级评价 1. 1 突变级数法的基本原理 法国数学家 Rene Thom 于 20 世纪 70 年代提出 了较为典型的突变理论。 后来,Zeeman 等在此基础 上进行了改进。 控制变量的个数 r 主要决定突变类 型当 r≤4 时,只有折叠突变、尖点突变、燕尾突变、 蝴蝶突变、椭圆突变、双曲突变以及抛物型突变 7 种 基本的突变模型和与之相对应的势函数形式[11-12]。 表 1 中是 4 种较为常见的突变系统。 其中,x 是 f x的状态变量;a、b、c 和 d 是函数的控制变量。 表 1 中具有高度概括性和普遍性的突变模型是 由 Rene Thom 经过严谨的计算推导得出的,尽管其 证明过程非常困难,但对其证明成果的运用则相对 比较容易。 现以尖点突变模型r =2为例进行相关 的计算推导。 求取尖点突变模型的势函数fx = x4+ ax2+ bx中状态变量 x 的一阶和二阶导数;首先, 通过 f ′x = 0 计算得到平衡曲面 M4x3+ 2ax + b =0;其次,通过 f ″x = 0 计算得到奇点集 S 12x2+2a = 0;最后,通过消去 x 计算得到分歧方程 B8a3+27b2=0。 奇点集 S 在控制变量 a、b 上的投 影构成了分歧集 B。 系统突变发生的条件是控制 变量 a、b 的取值,使方程 8a3+27b2=0 成立。 表 1 突变级数法常用的 4 种突变类型[13-16] Tab. 1 Four kinds of mutation types commonly used in the catastrophe progression 突变 模型 状态 变量 控制 变量 势函数分歧方程 折叠 突变 11fx = x3+ axa =0 尖点 突变 12fx = x4+ ax2+ bx a = -6x2, b =8x3 燕尾 突变 13 fx = x5+ ax3+ bx2+ cx a = -6x2, b =8x3,c = -3x4 蝴蝶 突变 14 fx = x6+ ax4+ bx3+ cx2+ d a = -10 x2,b =20 x3, c = -15x4,d =4x5 图 1 是突变级数法常见的 4 种突变模型。 利用奇点方程推导得到归一化公式,从而将控 制变量和状态变量的取值控制在 0 1 之间,满足传 统的模糊隶属函数取值规律,进而求出系统总突变 隶属函数值。 与表 1 中各突变模型相对应的归一化公式又 称突变级数如式1 式4所示[10] 折叠突变 a折叠突变 b尖点突变 c燕尾突变 d蝴蝶突变 图 1 突变级数法常见的 4 种突变模型示意图 Fig. 1 Schematic diagram of the four mutation models commonly used in catastrophe progression 06 爆 破 器 材 第 47 卷第 6 期 Xa= a 1 2; 1 尖点突变 Xa= a 1 2,Xb= b 1 3; 2 燕尾突变 Xa= a 1 2,Xb= b 1 3,Xc= c 1 4; 3 蝴蝶突变 Xa= a 1 2,Xb= b 1 3,Xc= c 1 4,Xd= d 1 5。 4 1. 2 突变级数法的一般评价步骤 突变级数法主要步骤包括[17] 1确定评价的总指标,对评价目标进行多层次 矛盾分解; 2把突变理论与模糊数学相结合,通过计算处 理得到突变模糊隶属函数; 3运用归一化公式进行综合量化运算,最后归 一为一个参数即求出总的隶属函数,根据所得结 果对评价目标进行排序分析。 1. 3 煤矿外因火灾危险等级评价 1. 3. 1 煤矿外因火灾的危险源分析 燃烧的发生与发展需要具备 3 个基本要素,即 可燃物、引火源和氧化物氧气。 分析矿井火灾发 生的实际情况和燃烧的三要素,消除可燃物和氧气 这两个条件的难度高、工程量大;因此,防止矿井外 因火灾发生的最直接、有效方法是杜绝引火源。 外因火灾是由外来热源引起的。 通过分析、查 阅相关资料,参考煤矿安全规程2016 修订,可 知造成煤矿外因火灾事故的火源主要包括 1井下明火,如吸烟、安全灯使用不规范、电焊 和气焊焊接起火、电弧火花等。 2电器设备起火,如井下各种机电设备由于长 期连续运转,在检修工作不到位的情况下非常容易 出现接触不良、漏电、短路、接地故障、电器故障火 花、电器照明故障、电线老化、脱皮等故障,造成设备 过流或摩擦过热,引发设备燃烧,引燃设备周围的可 燃物,设备一旦发生燃烧,极易造成火势的蔓延。 3其他原因起火,如常见的有违规放炮、运输 系统起火、瓦斯和煤尘爆炸、物体摩擦起火等[18]。 1. 3. 2 建立评价指标体系 煤矿火灾作为一个独特的系统工程,自身的危 险性特点以及动态特性决定了煤矿火灾发生的原因 多种多样[10]。 根据科学性原则、系统性原则、可操作性原则、 相对独立性原则、规范性原则、代表性原则,对引发 煤矿外因火灾的一般因素进行了分析,如明火致因 中的吸烟、电焊、火焊等;电火致因中的短路、电线老 化等;其他原因中的机械摩擦等。 并结合煤矿外因 火灾发生的特殊条件,如矿井外因火灾主要发生在 采空区、遗留的煤柱、破裂的煤壁、煤巷的高冒处以 及浮煤堆积的地点,它的发生和发展相对较慢、火源 潜伏期长、火源相对比较隐蔽等特点,而且井下火灾 一般是在空气极其有限的情况下发生的,特别是采 空区火灾和煤柱内火灾;即使发生在风流畅通的地 点,其空间和供氧条件也是有限的。 在文献[19]的 基础上,主要针对引火源这一火灾发生的必要条件, 建立煤矿外因火灾的安全评价指标体系。 指标体系 共有 15 个指标,构成 4 个突变系统。 其中 3 个二级 指标构成燕尾型突变系统,12 个三级指标构成 3 个 蝴蝶突变系统[10,18],见图 2。 图 2 煤矿外因火灾危险等级评价指标体系 Fig. 2 Index system for external fire risk grade assessment of coal mine 162018 年 12 月 基于突变级数法的煤矿外因火灾危险等级评价 贾进章,等 万方数据 1. 3. 3 评价集的建立[10] 将评价集划分为 V = {安全 V1较安全 V2一般 安全 V3不安全 V4}。 将其赋予相应的分值后,V = {90 80 70 60}。 显然,分数越高表示越安全。 1. 3. 4 隶属度及隶属函数的确定 确定因素集中各因素隶属于评价集 V 的程度, 借鉴文献[15] 的研究成果。 确定隶属函数如式 5 式8 f1u = 1,90≤u 100; u -80 90 -80, 80≤u 90; 0,u 80。 5 f2u = 100 - u 100 -90, 90≤u 100; 1,80≤u 90; u -70 80 -70, 70≤u 80; 0,u 70。 6 f3u = 0,90≤u 100; 90 - u 90 -80, 80≤u 90; 1,70≤u 80; u -60 70 -60, u 70。 7 f4u = 0,80≤u 100; 80 - u 80 -70, 70≤u 80; 1,u 70。 8 1. 3. 5 高级评价[10] 邀请相关领域的专家,对煤矿外因火灾的评价 指标进行评分,通过隶属函数,得到评价体系指标层 初始模糊隶属函数值矩阵 RU1、RU2、RU3。 以明火致 因 U1为例,将 U1下的三级指标 U11、U12、U13、U14的 评测值代入隶属函数 f1u、f2u、f3u、f4u,可 通过计算得到 RU1。 同理,将 U2、U3依次代入隶属 函数,计算可得 RU2、RU3。 以 U1为例进行计算。 初始模糊隶属函数值矩阵 RU1= f1U11f2U11f3U11f4U11 f1U12f2U12f3U12f4U12 f1U13f2U13f3U13f4U13 f1U14f2U14f3U14f4U14 。 9 用蝴蝶型的归一化公式对所有三级指标进行计 算,得到指标层的突变级数矩阵 R′U1、R′U2、R′U3。 指标层突变级数值矩阵见式10。 R′U1= f1U11f2U11f3U11f4U11 3 f1U12 3 f2U12 3 f3U12 3 f4U12 4 f1U13 4 f2U13 4 f3U13 4 f4U13 5 f1U14 5 f2U14 5 f3U14 5 f4U14 。 10 将突变级数的平均值作为状态变量 x 的值进行 递归计算,可得到二级指标 U1、U2、U3对应评语中 “安全” 的隶属函数值分别为 rf1 U1 、rf1 U2 、rf1 U3 ,如式 11 式13所示。 rf1 U1 = f1U11 + 3 f1U12 + 4 f1U13 + 5 f1U14 4 ; 11 rf1 U2 = f1U21 + 3 f1U22 + 4 f1U23 + 5 f1U24 4 ; 12 rf1 U3 = f1U31 + 3 f1U32 + 4 f1U33 3 。 13 将 U1、U2、U3隶属于其他评语的函数值进行排 列,可得二级指标准则层的模糊隶属函数值矩阵 RU= rf1 U1 rf2 U1 rf3 U1 rf4 U1 rf1 U2 rf2 U2 rf3 U2 rf4 U2 rf1 U3 rf2 U3 rf3 U3 rf4 U3 。14 根据燕尾型突变模型归一化公式对 U1、U2、U3 进行处理,可得准则层突变级数值矩阵 R′U= rf1 U1 rf2 U1 rf3 U1 rf4 U1 3 rf1 U2 3 rf2 U2 3 rf3 U2 3 rf4 U2 4 rf1 U3 4 rf2 U3 4 rf3 U3 4 rf4 U3 。15 将 R′U的第一列的平均值作为隶属于评价集中 “安全”的总突变隶属函数值,可得 r1= rf1 U1 + 3 rf1 U2 + 4 rf1 U3 3 。16 同理,可得隶属于评价集中其他评语的总突变 隶属函数值 r2、r3、r4。 总突变隶属函数值矩阵 R = [r1r2r3r4]。 依据总突变隶属函数值越大越优原 则,得出评价结论。 2 应用实例 以沈阳焦煤集团林盛煤矿为例,用基于突变级 数法的煤矿外因火灾危险性模型对其进行分析。 根 26 爆 破 器 材 第 47 卷第 6 期 万方数据 据该煤矿的实际情况,邀请辽宁工程技术大学矿井 火灾安全领域、东北大学矿井火灾安全领域以及中 国矿业大学煤矿安全领域的 10 位专家对可能引起 外因火灾的各项指标打分。 为使专家打出的评价分数尽可能准确反映评价 对象的重要程度,特意邀请的专家中,5 位长期从事 现场工作,5 位长期从事科研工作,取评分结果的算 术平均值作为评测数据[18],结果见表 2。 表 2 煤矿外因火灾评测参数 Tab. 2 uation parameters for external five of coal mine 目标层 二级指 标层 三级指标层 评测 数据 煤 矿 外 因 火 灾 危 险 性 评 价 指 标 体 系 U 明火致 因 U1 吸烟 U11 电焊、火焊、喷灯焊 U12 大灯泡取暖 U13 电弧火花 U14 75 77 80 68 电火致 因 U2 短路、过负荷、接地故障 U21 接触不良、漏电、照明故障 U22 电器故障火花 U23 电线老化、脱皮 U24 90 88 95 78 其他原 因 U3 违规放炮 U31 运输系统起火 U32 瓦斯、煤尘爆炸 U33 机械摩擦及物体碰撞 U34 71 61 85 87 将表 2 中各个三级指标的评测分值代入隶属函 数,利用相应模型的归一公式进行计算,得到 1指标层隶属函数值矩阵。 RU1= 00. 510. 5 00. 710. 3 010. 10 000. 81 ; RU2= 10. 100 0. 810. 20 10. 500 00. 810. 2 ; RU3= 00. 110. 9 0000 0. 510. 50 0. 710. 30 。 2指标层突变级数值矩阵 。 R′U1= 00. 70710. 707 00. 88810. 669 010. 5620 000. 9561 ; R′U2= 10. 31600 0. 92810. 5850 10. 84100 00. 95610. 725 ; R′U3= 00. 31610. 949 0000 0. 84110. 8410 0. 4430. 5790. 6570. 237 。 3准则层隶属函数值矩阵。 RU= 00. 6490. 8800. 594 0. 7320. 7780. 3960. 181 0. 4430. 5790. 6570. 237 。 4准则层突变级数值矩阵。 R′U= 00. 8060. 9380. 771 0. 9010. 9200. 7350. 566 0. 8160. 8720. 9000. 698 。 由此,得到该煤矿外因火灾危险性评价的总突 变隶属矩阵 R = [0. 572 0. 866 0. 858 0. 678] 。 依据总突变隶属函数值越大越优的原则,总突变隶 属函数值隶属于评语集中第二级的值最大,可知该 煤矿外因火灾危险等级为二级,即较安全。 表明该 矿火灾危险等级相对较好,可接受。 但隶属于第三 级的总突变隶属函数值为 0. 858,火灾危险等级有 升高的风险。 为了进一步降低煤矿外因火灾危险等级,同时 保障煤矿的安全生产和正常运营,根据评价的结果 及专家打分的数据提出相关对策定期进行安全检 查,针对电弧火花、违规放炮、运输系统起火及物体 碰撞的安全性相对较差的情况,需要加强安全管理 工作,加强对职工的培训,严格按照操作规程施工作 业;针对电焊、火焊、喷灯焊违规使用情况以及电器 设备的老化和脱皮情况,需要定期更换受损、老化的 设备;此外,还要尽量避免瓦斯、煤尘爆炸的发生。 3 结论 1从明火致因、电火致因、其他原因 3 个方面 出发,并根据实地调查分析以及相关国家法规和行 业标准,建立了煤矿外因火灾危险性评价指标体系, 指标体系共有 15 个指标,构成 4 个突变系统,其中 3 个二级指标构成燕尾型突变系统,12 个三级指标 构成 3 个蝴蝶突变系统。 2运用突变级数法对各级指标进行评测,能较 好地反映煤矿外因火灾事故过程中的突变特性。 利 用基于突变级数法建立的煤矿外因火灾危险性评价 362018 年 12 月 基于突变级数法的煤矿外因火灾危险等级评价 贾进章,等 万方数据 体系,以沈阳焦煤集团林盛煤矿为例进行了计算分 析,得到该煤矿外因火灾危险等级为较安全,并根据 评价结果,提出了提高煤矿外因火灾安全水平的技 术对策。 参 考 文 献 [1] 环境保护部环境应急与事故调查中心. 环境风险及化 学品检查技术报告[R]. 北京环境保护部,2012. [2] 刘铁民,张兴凯,刘功智. 安全评价方法应用指南 [M]. 北京化学工业出版社,2005. [3] DUNJ J,FTHENAKIS V, VLCHEZ J A, et al. Hazard and operability HAZOP analysis. A literature review [J]. Journal of Hazardous Materials,2010,1731/2/3 19-32. [4] CHUVIECO E, AGUADO I, YEBRA M, et al. Develop- ment of a framework for fire risk assessment using remote sensing and geographic ination system technologies [J]. Ecological Modelling,2010,221146-58. [5] WESSLING S, KUENZER C, KESSELS W, et al. Nu- merical modelling for analyzing thermal surface anomalies induced by underground coal fires[J]. International Jour- nal of Coal Geology,2008,743/4175-184. [6] 王旭,霍德利. 模糊综合评价法在煤矿安全评价中的 应用[J]. 中国矿业,2008,17575-78. WANG X, HUO D L. Application of fuzzy comprehensive uation in coal safety assessment[J]. China Mining Magazine,2008,17575-78. [7] 邓军,杨一帆,翟小伟. 基于层次分析法的矿井煤自燃 危险性模糊评价[J]. 中国安全生产科学技术, 2014, 104120-125. DENG J,YANG Y F, ZHAI X W. Fuzzy uation on risk of coal spontaneous combustion based on AHP[J]. Journal of Safety Science and Technology,2014,104 120-125. [8] 冯海明,金龙哲,张春芝. 基于神经网络的煤炭自燃预 测及在 MATLAB 上的实现[J]. 中国煤炭, 2008,34 582-84. [9] 蒋军成. 突变理论及其在安全工程中的应用[J]. 南京 化工大学学报,1999,21124-28. JIANG J C. Catastrophe theory and its application in safe- ty engineering [ J]. Journal of Nanjing University of Chemical Technology,1999,21124-28. [10] 张旺,冯涛,武剑,等. 基于突变级数法的煤矿瓦斯爆 炸危险性评价[J]. 矿业工程研究,2016,31141- 45. ZHANG W, FENG T, WU J, et al. Gas explosion risk assessment in coal mines based on catastrophe progres- sion [J]. Mineral Engineering Research, 2016, 31141-45. [11] 田冷,申智强,刘柳,等. 多层叠合致密砂岩气藏产量 劈分方法研究[J]. 天然气勘探与开发,2016,393 41-44. TIAN L, SHEN Z Q, LIU L, et al. Production split for multilayer-superimposed tight sandstone gas reservoirs[ J]. Natural Gas Exploration and Develop- ment, 2016,39341-44. [12] LIU M H. The research on catastrophe theory’s applica- tion in the re of sports curriculum[C] / /2010 Se- cond International Workshop on Education Technology and Computer Science. Wuhan, 2010790-793. [13] 邹义怀,江成玉,李春辉. 基于突变理论的煤与瓦斯 突出危险性预测的研究 [J]. 中国煤炭,2011,37 10105-108. ZOU Y H, JIANG C Y, LI C H. Research on risk prediction of coal and gas outburst based on catastrophe theory[J]. China Coal,2011,3710105-108. [14] 丁元春,翁发禄,唐敏康. 基于突变理论的高层建筑 火灾安全综合评价[J]. 建筑科学,2012,28122- 25. DING Y C, WENG F L, TANG M K. The research on catastrophe theory’s application in the fire risk assess- ment of high-rise building[J]. Building Science,2012, 28122-25. [15] 张俊朋,杨晓雯. 煤矿井外因火灾风险的灰类-FAHP 评价[J]. 安全与环境工程,2017,24489-93. ZHANG J P,YANG X W. Risk uation of exogenous mine fire based on grey clustering-FAHP [J]. Safety and Environmental Engineering,2017,24489-93. [16] 徐阿猛,齐黎明. 基于突变理论的马家沟矿煤与瓦斯 突出过程研究[J]. 煤炭工程,20111271-73. XU A M, QI L M. Study on coal and gas outburst process of Majiagou Mine base on catastrophe theory [J]. Coal Engineering,20111271-73. [17] 袁文瑞. 基于产业集群的工业园区竞争力研究[D]. 哈尔滨哈尔滨工业大学,2009. YUAN W R. Research on the competitiveness of indus- trial park based on industrial clusters[D]. HarbinHar- bin Institute of Technology,2009. [18] 贾进章,董晓雷. 基于模糊综合评价-集值统计法的煤 矿外因火灾危险性分析[J]. 安全与环境学报,2015, 15211-14. JIA J Z, DONG X L. Fuzzy uation and statistically oriented in analyzing the external causes of the coal mine fire[J]. Journal of Safety and Environment, 2015,15211-14. [19] 闫明仁,张人伟,段汉文,等. 煤矿外因火灾危险性的 模糊综合评价[J]. 煤矿安全,2009690-92. 46 爆 破 器 材 第 47 卷第 6 期 万方数据