Knowledge discovery and disaster chain network analysis of coal mine roof accidents
-
摘要:
为预防煤矿顶板事故的发生,发现事故致因的知识和致因灾害链的网络分析具有重要意义;以煤矿顶板事故报告中的事故原因为数据集,构建事故致因知识图谱,并基于复杂网络对致因灾害链进行网络分析。研究结果表明:顶板冒顶、空顶、支架变形等22个原因是导致煤矿顶板事故发生的主要致灾因子;灾害链网络分析中风险度最大的边为违章指挥→作业人员违章作业;对关键节点和边实施管控,能有效控制煤矿顶板事故的发生或衍生发展。
Abstract:In order to prevent the occurrence of coal mine roof accidents, it is of great significance to realize the knowledge discovery of accident causes and the network analysis of disaster chain. Taking the accident causes in the coal mine roof accident report as the data set, the knowledge map of accident causes was constructed. Further, based on the complex network theory, the causal disaster chain network is constructed and analyzed. The research results show that 22 factors, such as roof collapse, empty roof and support deformation, are the main disaster causing factors of coal mine roof accidents; the side with the greatest risk is from illegal command to illegal operation of operators in disaster chain network analysis. The control of key nodes and edges can effectively control the occurrence or derivative development of coal mine roof accidents.
-
Keywords:
- coal mine roof accident /
- complex network /
- knowledge discovery /
- disaster chain /
- safety management
-
矿井瓦斯是煤炭开采中常见的有害气体,对煤矿的安全生产构成严重威胁。根据国家矿山安监局的统计数据,从2001年到2020年的20年间共发生了484起煤与瓦斯突出事故,导致
3195 人死亡[1]。煤层气的开发不仅对矿井瓦斯防治具有重要意义,而且还具备一定的经济价值。因此,深入研究煤层致裂增透技术对于促进煤层气开发具有重要意义[2]。水力压裂技术提高了煤层气回收效率,但存在地层损害、大量水消耗等环境问题。鉴于水力压裂、水力冲孔等增透技术[3-4]的局限性,学者们开始广泛关注低温致裂增透技术[5-7]。最早是在1971年,由美国学者ALLEN and BAUER[8]提出向岩石裂缝中循环注入冷冻流体以达到增产效果;之后在1997年,学者MCDANIEL等[9]成功通过现场试验证明使用液氮重复致裂煤层能达到良好的增产效果,由此开始了低温致裂技术的研究。近些年来,李波等[10]研究发现不同冻结温度下,内部拉应力随温度下降而增加,大裂隙比小裂隙更早冻结;翟成等[11]研究发现随着冻融次数的增加,试样的中大孔占比、总孔隙度、有效孔隙度和渗透率显著增加,煤体透气性增强;张春会等[12]研究发现饱水程度对液氮溶浸致裂效果影响很大;张磊等[13]研究发现液氮溶浸处理可显著提升煤体的渗透率,但随着溶浸时间的增长,液氮的致裂增透效果是在逐渐衰减的。在低温致裂的研究中,冷冲击饱水煤体致裂的力学特性及渗流规律有待完善。
为此,开展了不同温度冷冲击饱水试样致裂的渗流试验研究,研究结果对煤层增透和煤层气抽采有一定的参考价值。
1. 试验方法
1.1 试样制备
试样取自新景矿[14],均来自同一块大尺寸原煤,其工业分析指标结果如下(质量分数):水分2.26%,灰分19.6%,挥发分10.64%,固定碳67.89%。根据实验室岩样制备标准,按GB/T 23561.7—2009《煤和岩石物理力学性质测定方法》[15],进行试样的制备。步骤如下:利用水力台式钻机钻取尺寸为ϕ50 mm$ \times $100 mm的标准试样;经水力钻孔所取试样表面较为粗糙,因此,手工打磨试样表面,目的是使试样表面受力均匀;用保鲜膜包裹,放置于实验室备用。
使用负压饱水装置对试样进行饱水处理,每隔12 h对试样称重,直至试样质量不再变化,得到最大含水率试样。
搭建冷冲击处理平台,平台包括温控系统、试验箱、声发射系统。其优点为恒温控制,精度高、易操作。温度控制范围为−45~190 ℃,温度控制精度±0.5 ℃,升降温速率0.7~1 ℃/min。具体试样参数及冲击温度见表1。
表 1 试样参数及冲击温度Table 1. Specimen parameters and impact temperatures试样编号 尺寸/(mm×mm) 冲击温度/ ℃ 1# ϕ50.22×102.01 — 2# ϕ50.12×101.53 0 3# ϕ50.34×101.23 −20 4# ϕ50.24×101.06 −40 1.2 试验内容和试验步骤
采用WYS-800微机控制电液伺服三轴试验装置进行测试,研究饱水试样冷冲击处理后渗流规律及力学特性。
试验采用氮气,瓦斯压力取0.5 MPa。在轴压为3.06、8.15、13.25、18.34、23.44、28.54 MPa时,测量试样的渗透率,直至试样被压碎。同时,使用声发射监测仪监测应力加载过程中试样的声发射信号。
试验步骤如下:
1)安装饱水处理试样1。将试样装入热缩管,用热风枪烘烤热缩管,以防止渗流测试气体从煤壁外侧通过,用704硅胶进行密封处理。安装环向引伸计,放入三轴室内。
2)将试样安装好后加油排空三轴压力室内的空气,对试样加载轴向力到预定值,加载围压到预定值(加载方式:围压轴压交替加载)。
3)吸附氮气。试验前需对试样进行抽真空处理,通入氮气到预定压力,待12 h后视为氮气吸附平衡。
4)渗透率测试。在试样吸附平衡后,打开出气管阀门,待流速稳定之后,测定流量数据,重复测量3次,然后取平均值,作为试样的初始渗透率。
5)在三轴室外壁安装声发射探头,监测应力加载过程的声发射信号。继续加载轴向应力,待吸附平衡后再次测定流量,直至试样被压碎,在试验过程中数据自动记录。
6)将低温饱水处理试样2~试样4,重复步骤1)~步骤4),完成其余试样的渗透率测试及三轴力学试验。试验流程图如图1所示。
2. 试验结果
2.1 冷冲击饱水试样三轴力学特性
对不同温度冲击的试样进行三轴压缩试验,得到的试样的峰值强度及应变见表2,常温及不同温度冷冲击处理试样应力−应变曲线如图1所示。在轴向应力的加载过程中试样的变形分为4个阶段:压密阶段、弹性阶段、屈服阶段、破坏阶段。
表 2 试样的峰值强度及应变Table 2. Peak strength and strain of specimens试样 峰值强度/MPa 峰值轴向应变 峰值径向应变 1# 35.128 66 0.017 34 −0.059 28 2# 29.859 87 0.018 67 −0.047 44 3# 24.866 24 0.011 76 −0.022 80 4# 24.570 70 0.013 43 −0.037 12 由图1可以看出,不同冷冲击温度处理下的试样及常温下未处理试样的轴向应力-轴向应变趋势基本一致。随着冷冲击处理温度的降低,试样的峰值强度降低。由煤体内外温差形成温度梯度导致矿物颗粒变形相互制约产生的热应力,计算如下[16-17]:
$$ {\sigma _{ij}} = {\alpha _{ij}}{E_{ij}}\Delta T{\delta _{ij}} $$ (1) 式中:${\sigma _{ij}}$为热应力,MPa;${\alpha _{ij}}$为煤体线膨胀系数;${E_{ij}}$为煤体弹性模量;$\Delta T$为温度变化, ℃;${\delta _{ij}}$为Kronecker符号,取值为$1$。
根据热应力公式,温差越大,产生的热应力越大,对煤体的造成的破坏越大,其抗压强度也会越低。因此随着冷冲击温度降低,试样的抗压强度逐渐降低。经冷冲击处理后,相同轴向应力下,试样的轴向应变大于未经处理的常温试样。
而与之相反,经冷冲击处理后,在相同轴向应力下,试样的径向应变小于未经处理的常温试样。
使用Origin软件对不同冷冲击处理试样的温度与峰值强度进行拟合,冷冲击温度与峰值强度拟合曲线如图2所示,由图2可以看出,随着冷冲击温度的降低,试样的峰值强度逐渐减小。
为进一步研究冷冲击处理饱水试样的结构劣化,在三轴压缩加载的同时,进行声发射定位,记录试样在不同应力阶段下声发射定位图,即峰值应力的20%、60%、100% 3个阶段。不同温度冷冲击煤体三轴压缩声发射定位图如图3所示,图中数字1~6分别为声发射探头的6个定位点。
将定位图分为初始期、活跃期、峰值前后3个阶段,在处理试样达到峰值应力的20%之前,应力−应变曲线处于弹性阶段前期,内部声发射事件较为稀疏,此过程产生的声发射事件主要来源于试样孔裂隙以及冷热冲击造成的结构弱面压缩破坏;继续加载过程中,当达到冷热冲击试样峰值应力的60%之前,应力−应变曲线处于弹性阶段后期,孔裂隙开始发育,声发射事件增加;峰值应力前后,冷热冲击试样裂隙已完全发育贯通,达到峰值应力后试样破坏。
在应力加载的各个阶段,冷冲击的温差越大,定位图中记录的声发射事件越密集。冷冲击温差大的试样,在应力加载到下一阶段时,定位图中的声发射事件增加的幅度大于温差小的试样。由此可见,冷冲击温差越大,对饱水试样造成的损伤越大。
2.2 冷冲击温度对煤体渗透率的影响
低渗透煤层主要依赖微孔裂隙通道,其渗透率较低。在气体渗透时,随着平均孔隙压力增加,渗透率减小,表现出滑脱效应。温度冲击会导致微孔裂隙重构,改变滑脱效应程度,反映煤体的损伤演化。在试验中,通过对煤样进行负压饱水处理,使水分充满煤体内部的孔裂隙。渗透率计算公式如下:
$$ K = \frac{{2Q{p_0}\mu L}}{{(p_1^2 - p_2^2)S}} $$ (2) 式中:$K$为渗透率,${{\mathrm{m}}^2}$;${p_0}$为标准大气压,${\mathrm{MPa}}$;$Q$为流量,cm3/s;$\mu $为动力黏度系数,${\mathrm{MPa}} \cdot {\mathrm{s}}$;$L$为试样长度,${\mathrm{cm}}$;${p_1}$为进气口压力,${\mathrm{MPa}}$;${p_2}$为出气口压力,${\mathrm{MPa}}$;$S$为试样渗流断面面积,${{\mathrm{cm}}^2}$。
常温及不同温度冷冲击处理试样轴压−渗透率曲线如图4所示;常温及不同温度冷冲击处理试样轴向应力−应变−渗透率曲线如图5所示;常温及不同温度冷冲击处理试样渗透率表见表3。
由图4可以看出:不同冷冲击温度处理的试样渗透率随轴向应力的变化趋势一致,随着轴向应力增大,试样的渗透率逐渐减小;试样的初始渗透率随着冷冲击温度的降低而增大,−40 ℃冷冲击下的试样初始渗透率达到最大,表明其增透效果最好;但是随着轴压增大,−40 ℃冷冲击处理的试样渗透率会有个明显的上升趋势。
煤体内部孔−裂隙结构影响着煤体渗透率,试样受载时,内部结构发生变化,导致渗透率改变。由图5可以看出,其变化趋势基本一致,渗透率随着应变的增加而减小。试样在受载阶段的渗透率变化规律为:在压密阶段,轴向载荷增大,孔−裂隙结构逐渐闭合,允许气体通过的通道变窄,渗透率减小。应力继续加载,试样进入弹性变形阶段,渗透率不断下降,下降的速率减小。当加载应力超过试样局部结构承受极限时,裂隙结构不断发育,各裂隙连接贯通,渗透率增大。在屈服阶段,裂纹快速发育,形成宏观破裂面,达到峰值强度,试样整体破碎。
由表3可知:经0 ℃冷冲击处理后,试样的初始渗透率较未经处理常温试样的渗透率增加了15.89%;−20 ℃冷冲试样较常温下试样的渗透率增加了39.19%,相比0 ℃冷冲击下的试样,初始渗透率的增长为20.1%,增长幅度不大;−40 ℃冷冲击处理的试样较常温试样渗透率增长为202.9%,较−20 ℃冷冲击,初始渗透率增长117.2%,增长幅度较大。饱水试样经冷冲击处理后,水分冻结成冰,导致体膨胀,扩大了瓦斯流动通道。当膨胀应力超过煤体的承受极限时,试样内部结构受损,裂隙扩展,渗透率增加。由于冷冲击处理引起煤体的不均匀变形,产生热应力,且温度差异引发温度场,煤体内外温差导致相互制约的热应力。当热应力超过煤体的抗拉强度时,内部裂隙扩展并连接,渗透率增大。冷冲击产生的热应力与冷冲击温度成正比,即温度越低,热应力越大,试样渗透率增大。因此,随着冷冲击温度的降低,试样的渗透率不断增加。
表 3 常温及不同温度冷冲击处理试样渗透率表Table 3. Permeability table of samples treated with cold shock at normal temperature and different temperatures试样 冷冲击
温度/ ℃初始渗透率/
10−15 m2较常温试样
渗透率增幅/%较上一档冷冲击
处理试样的
渗透率增幅/%1 常温 0.096 6 — — 2 0 0.111 9 15.89 15.89 3 −20 0.134 4 39.19 20.10 4 −40 0.292 5 202.90 117.62 3. 结 语
1)经冷冲击处理后,试样的峰值强度降低,随着冷冲击温度的降低,试样的峰值强度越低。不同温度冷冲击处理的试样的应力−应变曲线变化趋势基本一致。
2)试样的渗透率随冷冲击温度的降低而增大,与试样的峰值强度变化吻合。冷冲击的温度越低,温差越大,产生的热应力越大,对试样的损伤越强,煤体内部的裂隙发育扩张越剧烈,导致渗透率增大。
3)−40 ℃冷冲击下试样的初始渗透率最大,但随着应力加载,渗透率减小的速度最快。−20~−40 ℃之间,初始渗透率的增长幅度最大。
-
表 1 直接致因集群关键词表
Table 1 Keywords table of cluster of direct causes
集群 关键词 主要内容 1# 破碎、液压、支架、应力、地质、压死、岩石、综采 顶板破碎和支架 2# 导致、大面积、垮落、受伤、巷道 顶板垮落 3# 巷道、顶板、支护、强度、冒顶、压力、加固、维修 顶板支护和支护强度 4# 事故、致死、措施、规程、工人、违反、支柱 事故的致死和措施规程 5# 掘进、空顶、锚杆、锚索、失稳、离层、施工 空顶作业、违章掘进和锚杆锚索 6# 工作面、作业、现场、人员、违章、冒险、发生 工作面和违章作业 表 2 间接致因集群的关键词表
Table 2 Keywords table of indirect cause cluster
集群 关键词 主要内容 1# 工作面、顶板、事故、支护、支架、巷道、施工 工作面的支护 2# 针对性、措施、技术、规程、编制、不严 作业规程和技术措施 3# 作业、人员、职工、安全意识、风险、淡薄、教育、培训、 人员的教育培训和安全意识 4# 现场、管理、执行、管理人员、安全检查、领导、检查、职责 管理人员和现场管理 5# 制度、严格执行、敲帮、问顶、管理混乱、违反、安全措施 制度和措施的执行 6# 未能、消除、采取有效、事故隐患、防护 事故隐患和有效措施 7# 隐患、排查、落实、责任、到位、安全监管、整改、督促 隐患排查和责任落实 表 3 煤矿顶板事故的致灾因子
Table 3 Disaster causing factors of coal mine roof accidents
编号 致灾因子 编号 致灾因子 ${V_1}$ 顶板冒顶 ${V_{12}}$ 违章指挥 ${V_2}$ 顶板破碎 ${V_{13}}$ 锚杆锚索失效 ${V_3}$ 空顶 ${V_{14}}$ 安全措施不到位 ${V_4}$ 矸石离层掉落 ${V_{15}}$ 未遵守作业规程 ${V_5}$ 支护作业不可靠 ${V_{16}}$ 培训教育不到位 ${V_6}$ 顶板垮落 ${V_{17}}$ 管理人员作风不实 ${V_7}$ 液压支架工作阻力不足 ${V_{18}}$ 未制定针对性措施 ${V_8}$ 支架变形 ${V_{19}}$ 未能排查隐患 ${V_9}$ 工作面处理不当 ${V_{20}}$ 责任落实不到位 ${V_{10}}$ 作业人员违章作业 ${V_{21}}$ 管理混乱 ${V_{11}}$ 未敲帮问顶 ${V_{22}}$ 安全意识淡薄 表 4 节点的度中心性
Table 4 Degree centrality of nodes
节点 $D_{{\mathrm{C}}i}$ 节点 $D_{{\mathrm{C}}i}$ 节点 $D_{{\mathrm{C}}i}$ 节点 $D_{{\mathrm{C}}{{i}}}$ ${V_1}$ 0.286 ${V_7}$ 0.439 ${V_{13}}$ 0.190 ${V_{19}}$ 0.619 ${V_2}$ 0.286 ${V_8}$ 0.143 ${V_{14}}$ 0.524 ${V_{20}}$ 0.190 ${V_3}$ 0.238 ${V_9}$ 0.619 ${V_{15}}$ 0.238 ${V_{21}}$ 0.286 ${V_4}$ 0.381 ${V_{10}}$ 0.667 ${V_{16}}$ 0.476 ${V_{22}}$ 0.524 ${V_5}$ 0.619 ${V_{11}}$ 0.286 ${V_{17}}$ 0.381 ${V_6}$ 0.476 ${V_{12}}$ 0.190 ${V_{18}}$ 0.524 表 5 节点的综合重要度
Table 5 Comprehensive importance of nodes
节点 $I_{i}$ 节点 $I_{i}$ 节点 $I_{i}$ 节点 $I_{i}$ ${V_1}$ 0.137 ${V_7}$ 0.406 ${V_{13}}$ 0.048 ${V_{19}}$ 0.771 ${V_2}$ 0.137 ${V_8}$ 0.044 ${V_{14}}$ 0.587 ${V_{20}}$ 0.045 ${V_3}$ 0.091 ${V_9}$ 0.752 ${V_{15}}$ 0.091 ${V_{21}}$ 0.178 ${V_4}$ 0.227 ${V_{10}}$ 1.00 ${V_{16}}$ 0.403 ${V_{22}}$ 0.514 ${V_5}$ 0.766 ${V_{11}}$ 0.151 ${V_{17}}$ 0.344 ${V_6}$ 0.318 ${V_{12}}$ 0.098 ${V_{18}}$ 0.605 表 6 边介数及风险度
Table 6 Edge betweenness and risk degree
边:边介数&风险度 边:边介数&风险度 边:边介数&风险度 边:边介数&风险度 边:边介数&风险度 $\left( {12,10} \right)$:14.250&5.597 $\left( {19,5} \right)$:3.833&2.277 $\left( {10,14} \right)$:6.000&1.358 $\left( {22,19} \right)$:2.310&1.174 $\left( {20,18} \right)$:5.500&0.801 $\left( {21,9} \right)$:5.167&2.847 $\left( {21,18} \right)$:4.667&1.899 $\left( {22,9} \right)$:2.643&1.271 $\left( {18,22} \right)$:8.583&1.168 $\left( {15,5} \right)$:3.500&0.801 $\left( {15,10} \right)$:8.000&2.569 $\left( {21,14} \right)$:4.083&1.483 $\left( {16,10} \right)$:4.000&1.225 $\left( {8,5} \right)$:8.000&0.975 $\left( {9,5} \right)$:2.333&0.760 $\left( {20,19} \right)$:6.500&2.563 $\left( {21,19} \right)$:3.333&1.367 $\left( {14,5} \right)$:4.833&1.200 $\left( {16,22} \right)$:2.333&0.947 $\left( {10,6} \right)$:2.643&0.746 $\left( {22,10} \right)$:4.643&2.464 $\left( {18,14} \right)$:3.000&1.362 $\left( {17,19} \right)$:3.833&1.181 $\left( {10,5} \right)$:2.000&0.887 $\left( {19,7} \right)$:4.917&0.706 -
[1] 康红普,张镇,黄志增. 我国煤矿顶板灾害的特点及防控技术[J]. 煤矿安全,2020,51(10):24−33. KANG Hongpu, ZHANG Zhen, HUANG Zhizeng. Characteristics of roof disasters and controlling techniques of coal mine in China[J]. Safety in Coal Mines, 2020, 51(10): 24−33.
[2] 明崯崯,李亨英,李忠卫. 基于ISM法对煤矿顶板事故致因研究[J]. 煤矿安全,2013,44(9):238−240. MING Yinyin, LI Hengying, LI Zhongwei. The research on causing factors of mine roof accidents based on ISM method[J]. Safety in Coal Mines, 2013, 44(9): 238−240.
[3] 秦岩,盛武. 基于贝叶斯网络的煤矿顶板事故致因研究[J]. 矿业安全与环保,2022,49(3):136−142. QIN Yan, SHENG Wu. Study on causes of coal mine roof accidents based on Bayesian network[J]. Mining Safety & Environmental Protection, 2022, 49(3): 136−142.
[4] 邢艳冬,周建亮,孟凡旺,等. 基于时空的煤矿顶板事故中施工作业耦合网络研究[J]. 煤矿安全,2023,54(6):144−152. XING Yandong, ZHOU Jianliang, MENG Fanwang, et al. Research on construction operation coupling network in coal mine roof accidents based on spacetime[J]. Safety in Coal Mines, 2023, 54(6): 144−152.
[5] 陈晓勇,施式亮,李润求,等. 基于G1-DEMATEL法的煤矿顶板事故风险因素分析[J]. 湖南科技大学学报(自然科学版),2018,33(1):16−22. CHEN Xiaoyong, SHI Shiliang, LI Runqiu, et al. Study on risk factors of coal mine roof accident based on G1-DEMATEL method[J]. Journal of Hunan University of Science and Technology (Natural Science Edition), 2018, 33(1): 16−22.
[6] 田水承,周鹏辉. 基于文本挖掘的煤矿顶板事故致因网络分析[J]. 煤炭技术,2023,42(2):117−121. TIAN Shuicheng, ZHOU Penghui. Network analysis of coal mine roof accident causation based on text mining[J]. Coal Technology, 2023, 42(2): 117−121.
[7] 黄昕,靳健,林作忠,等. 基于复杂网络的城市深部空间火灾灾害链分析及智慧减灾方法研究[J]. 应用基础与工程科学学报,2021,29(5):1280−1291. HUANG Xin, JIN Jian, LIN Zuozhong, et al. Complex network-based analysis on fire disaster chain of urban deep underground space and intelligent disaster mitigation measures[J]. Journal of Basic Science and Engineering, 2021, 29(5): 1280−1291.
[8] 郭文亚,周志鹏. 地铁工程建设安全事故致因网络模型构建与解析[J]. 中国安全科学学报,2020,30(1):155−161. GUO Wenya, ZHOU Zhipeng. Subway construction accident causation network modeling and analysis[J]. China Safety Science Journal, 2020, 30(1): 155−161.
[9] 李丽,陈志平,刘新,等. 复杂通风网络矿井反风技术分析及火灾应急对策[J]. 煤矿安全,2019,50(10):177−180. LI Li, CHEN Zhiping, LIU Xin, et al. Analysis of anti-wind technology with complex ventilation network mine and fire emergency countermeasures[J]. Safety in Coal Mines, 2019, 50(10): 177−180.
[10] 杨应柳,晋良海,邵波,等. 基于复杂网络的煤矿火灾爆炸致因研究[J]. 中国安全科学学报,2023,33(1):145−151. YANG Yingliu, JIN Lianghai, SHAO Bo, et al. Research on causes of coal mine fire and explosion based on complex network[J]. China Safety Science Journal, 2023, 33(1): 145−151.