Prediction of water inrush source of coal seam floor based on Fisher discriminant model
-
摘要: 针对传统矿井突水水源判别方法对矿井煤层底板突水水源判别准确率低的问题，以城郊煤矿二2煤层为例，建立了Fisher矿井底板突水水源判别分析模型。城郊煤矿二2煤层具有突水威胁的含水层分别为煤系砂岩含水层和底板太原组岩溶裂隙含水层，考虑到水化学离子的重要性及数据的有效性，采用煤层底板有突水威胁的砂岩水、太灰水和混合水3类水样分析资料作为样本，选取Ca2+，Mg2+，Na++K+，HCO3−，Cl−，SO42−这6种离子含量和矿化度作为矿井突水水源判别分析的变量。利用SPSS软件求得2个典型Fisher判别函数（第1和第2判别函数），计算出典型判别函数在3类水样类别的中心值，通过比较待判水样函数值与这3类水样类别的中心值距离即可判别样本归属哪一类别。利用回代估计法对Fisher矿井底板突水水源判别分析模型进行检验，结果表明：该模型的判别正确率达93.3%，判别结果可信度高。利用该模型对城郊煤矿二2煤层10个已知水样进行分类，得出10个水样的判别结果与实际情况吻合，判别正确率为100%。Abstract: In order to solve the problems of low accuracy mine water inrush source discriminant method for mine floor water inrush source discrimination, taking the second level coal seam of suburban coal mine as an example, the Fisher mine floor water inrush source discriminant model is established. The aquifers with the threat of water inrush in the second level coal seam of suburban coal mine are the coal-measure sandstone aquifer and the karst fractured aquifer of the Taiyuan Formation in the floor. Considering the importance of hydrochemical ions and the validity of the data, three kinds of water quality analysis data of sandstone water, limestone water and mixed water with water inrush threat in the coal seam floor are used as samples. The content and mineralization of six kinds of ions, Ca2+, Mg2+, Na++K+, HCO3−, Cl− and SO42−, are selected as the discriminant analysis variables for the identification of mine inrush water sources. Two typical Fisher discriminant functions (the first and the second discriminant functions) are obtained by SPSS software. The central values of the typical discriminant functions in the three water quality groups are calculated. By comparing the distance between the function values of the water samples to be discriminated and the central values of the three water quality groups, it is able to determine which group the samples belong to. The back substitution estimation method is used to test the Fisher mine floor water inrush source discriminant model. The results show that the discriminant accuracy rate of the model is 93.3%, and the discriminant results are highly reliable. The model is used to classify 10 known water samples in the second level of suburban coal mine. The results show that the discriminant effect of 10 water samples is consistent with the actual situation, and the discriminant accuracy rate is 100%.
-
0. 引言
快速有效地预测瓦斯浓度可及时发现井下灾害隐患、组织人员撤离，有效保障井下人员和设备安全。因此，提高对瓦斯浓度的预测预警效率，增强对瓦斯浓度变化的实时感知能力，始终是煤矿瓦斯监测需要解决的问题。
常用的瓦斯浓度预测方法包括灰色模型（Gray Model，GM）、混沌时间序列分析、支持向量机（Support Vector Machine，SVM）和神经网络4种。GM先对瓦斯浓度的变化规律进行动态理论分析，再利用灰色系统控制理论建立合适的模型对数据进行有效预测。混沌时间序列分析使用混沌时间序列理论，结合相空间重构找到影响因素和瓦斯浓度变化的关系，从而建立模型对瓦斯数据进行预测。SVM能够较好地解决应用中常见的过拟合、陷入局部极小值或模型泛化能力差等问题，采用小波算法、相空间重构等对瓦斯数据进行预处理。神经网络通过结合影响瓦斯浓度变化的因素指标，将这些因素指标或瓦斯自身数据或二者结合作为输入信息，通过建立合适的神经网络模型对瓦斯浓度进行预测。其中，常用于瓦斯浓度预测任务的神经网络包括反向传播神经网络（Back Propagation Neural Network，BPNN）、径向基神经网络、极限学习机等，多数预测模型通过结合不同的优化算法，优化模型的网络参数或模型结构，提高瓦斯浓度预测精度[1]。
为解决煤矿井下复杂环境下瓦斯浓度预测精度较低的问题，许多学者对传统的瓦斯浓度预测模型进行了优化。王雨虹等[2]将蝗虫优化算法（Improved Grasshopper Optimisation Algorithm，IGOA） 与长短时记忆网络（Long Short Term Memory，LSTM）相结合，提出一种基于蝗虫优化算法的瓦斯浓度预测模型，实验结果表明该模型具有更高的预测性能。戚昱[3]采用遗传算法（Genetic Algorithm，GA）对BP神经网络进行优化，提出了基于信息融合和GA−BP的煤矿瓦斯浓度预测模型，经试验验证发现该模型的预测误差较低，提高了瓦斯浓度预测的准确率。冉啟华等[4]在门控循环单元（Gated Recurrent Unit，GRU）的基础上，提出一种基于A−GRU的瓦斯浓度序列预测模型，通过调整网络隐藏层参数，使得误差损失最小化，提高了瓦斯浓度的预测精度。梁运培等[5]利用布谷鸟搜索（Cuckoo Search，CS）算法对LSTM进行优化，建立了基于CS−LSTM的瓦斯浓度预测模型，通过对隐藏层层数等4个超参数进行寻优，提升了全局寻优能力和模型预测精度。付华等[6]结合皮尔逊相关系数法提出一种基于多传感器的深度长短期记忆网络（Long Short-Term Memory，LSTM）瓦斯浓度预测模型，实验结果表明该模型在多步滚动预测等方面效果较好。上述文献从时间维度进行建模，虽然通过引入各种优化算法提高了瓦斯浓度预测精度，但都忽略了瓦斯浓度数据的空间特性。瓦斯运移和分布规律与时间、空间密切相关[7]，仅从时间角度处理易导致重要先验知识丢失，影响预测效果。
针对上述问题，本文将图卷积神经网络（Graph Convolutional Networks，GCN）[8]与门控循环单元(Gated Recurrent Unit，GRU)[9]有机结合，构建基于GCN−GRU的瓦斯浓度时空分布预测模型，并对瓦斯浓度进行预测。首先，分析监测区域内瓦斯浓度的时空相关性。其次，利用GCN对煤矿井下瓦斯监测节点图进行处理，得到相应的空间依赖关系，即瓦斯空间特征信息。再次，利用GRU对时间序列下各时刻组成的瓦斯空间特征信息进行处理，通过基于序列到序列模型和自动编码器，生成模型预测结果。最后，通过预测结果与井下瓦斯浓度实际值进行对比，验证模型的预测精度。
1. 瓦斯浓度数据的时空特性分析
随着时间变化而变化的数据称为时间序列数据，随着时间和空间2个维度同时发生变化的数据称为时空数据。时空数据具有动态变化、高维度和非线性等特征[10]。对瓦斯浓度序列数据特性进行研究是瓦斯时空序列分析与建模的基础，所以本文先对瓦斯数据的时间特性和空间特性分别进行研究。
1.1 时间特性
某煤矿采掘工作面瓦斯体积分数变化情况如图1所示。由图1（a）可看出，瓦斯体积分数整体随时间呈周期性波动，每个周期内的变化情况存在一定差异，如瓦斯波动范围不同等。由图1（b）可看出，在1个周期内，瓦斯体积分数随时间发生非线性变化，趋势性特征明显。结合不同的时间尺度来看，瓦斯体积分数变化均存在一定的规律性，受煤矿井下实际开采条件等因素的影响，某一时刻的瓦斯体积分数与前一段历史时间内的体积分数密切相关。因此，为了对瓦斯体积分数进行长期预测，需对瓦斯体积分数进行时间依赖建模，捕捉瓦斯体积分数在时间维度的依赖关系。
1.2 空间特性
开采工作面的瓦斯主要来源于工作面附近煤层、岩层、落煤区域等。在工作面巷道内，瓦斯在扩散运移过程中主要有2种运动形式：一种是来自通风系统的空气与瓦斯之间存在一定程度的浓度差和温度差，根据菲克定律，瓦斯会发生扩散运动；另一种是巷道内瓦斯气体随风流运动。
以U型巷道内瓦斯的扩散运动为例，对瓦斯运移规律进行分析。假设瓦斯气体由体积单元组成，设任意体积单元边长分别为$ \Delta x,\Delta y,\Delta {\textit{z}} $，如图2所示，x，y，z分别为水平方向、竖直方向、垂直方向。
为了便于分析，假设气体只向x轴正方向移动，令 dy =dz =1。单位时间内进入体积单元的瓦斯量为
$$ \Delta {Q_{{\rm{in}}}} = EJ\Delta y\Delta {\textit{z}}\Delta t - P\frac{{\partial J}}{{\partial x}}\Delta y\Delta {\textit{z}}\Delta t $$ (1) 式中：$ E $为瓦斯气体的平均流速；$ J $为体积单元中的平均瓦斯体积分数；$ \Delta t $为瓦斯扩散时间；P为瓦斯扩散系数。
单位时间内从体积单元流出的瓦斯量为
$$ \begin{split} \Delta {Q_{{\rm{out}}}} =& \left(EJ + \frac{{\partial EJ}}{{\partial x}}\Delta x\right)\Delta y\Delta {\textit{z}}\Delta t - \\& \left(P\frac{{\partial J}}{{\partial x}} + \frac{\partial }{{\partial x}}\left(P\frac{{\partial J}}{{\partial x}}\right)\Delta x\right)\Delta y\Delta {\textit{z}}\Delta t \end{split} $$ (2) 瓦斯量在体积单元中会发生一定量的衰减，单位时间内由于衰减而减少的瓦斯量为
$$ \Delta {Q_{\rm{J}}} = pJ\Delta x\Delta y\Delta {\textit{z}}\Delta t $$ (3) 式中p为衰减速率，是一个常数。
体积单元内最终的瓦斯量为
$$ \begin{split} \Delta Q =& \left(EJ - P\frac{{\partial J}}{{\partial x}}\right)\Delta y\Delta {\textit{z}}\Delta t - \Bigg(EJ + \frac{{\partial EJ}}{{\partial x}}\Delta x - \\& \left(P\frac{{\partial J}}{{\partial x}} + \frac{\partial }{{\partial x}}\left(P\frac{{\partial J}}{{\partial x}}\right)\Delta x\right)\Bigg)\Delta y\Delta {\textit{z}}\Delta t - pJ\Delta x\Delta y\Delta {\textit{z}}\Delta t \end{split} $$ (4) 对式（4）两边求导得
$$ \frac{{\partial J}}{{\partial t}} = P\frac{{{\partial ^2}J}}{{\partial {x^2}}} - E\frac{{\partial J}}{{\partial x}} - pJ $$ (5) 式（5）表达了工作面瓦斯浓度与时间$t$和空间位置$x$之间的关系，说明瓦斯浓度序列具有时空相关性。
1.3 瓦斯传感器的布置
根据文献[11]，瓦斯传感器应垂直悬挂于顶板、巷道侧壁附近等区域。以煤矿井下U型通风方式为例，将瓦斯传感器布设于同一采掘区域的上隅角A、工作面进风巷B、2处工作面回风巷C和D处，如图3所示。本文使用的瓦斯浓度数据来自这4个监测节点，采集区域空间结构连通。瓦斯浓度数据可以定义在节点图上，各监测节点构成的图结构如图4所示，整体由$m$个图结构数据帧组成，表示在时间维度上拓扑结构的变化。
2. 瓦斯浓度预测模型的建立
瓦斯浓度是关于时间和空间位置的函数，根据工作面瓦斯浓度序列的时空相关性，建立煤矿井下瓦斯浓度时空分布预测模型。
2.1 问题定义
定义$X \in {{\bf{R}}^{N \times M}}$为模型输入，即初始图信号，其中$N$为瓦斯监测节点数目，$M$为每个节点的特征数量，通过构建瓦斯监测节点图$ G = (v,e) $对节点间复杂的依赖关系进行建模，其中$ v $为瓦斯监测站点集，$ |v| = N $，$ e $为边集，表示传感器之间的距离。
瓦斯浓度预测的任务为建立瓦斯浓度预测模型$i( \cdot )$，基于节点图G，对未来时间段${\tau _{{\rm{out}}}}$所有监测节点的瓦斯浓度$O \in {{\bf{R}}^{N \times {\tau _{out}}}}$进行预测。
$$ \left\{ {X;G} \right\}\xrightarrow{{i( \cdot )}}O $$ (6) 2.2 框架分析
瓦斯浓度预测模型整体框架如图5所示。首先，对瓦斯浓度历史数据进行预处理，根据各采集节点间的空间距离，构建瓦斯浓度空间节点图，用于对节点间复杂的依赖关系进行建模。然后，在每个采样时间点，将瓦斯浓度和节点间的距离权重参数作为输入，在获取瓦斯的空间节点图结构之后，通过GCN进行空间特征自适应学习和图卷积运算，获得瓦斯的空间特征[12]，再将获取的空间特征信息转换为序列数据，输入到GRU。最后，通过基于GRU的编码器对提取到的节点信息进行编码，并将编码结果输入解码器，在解码器中使用与编码器相同的参数矩阵进行运算，在解码之后获得预测输出。
2.3 基于图的空间依赖建模
在图信号处理领域，卷积神经网络（Convolutional Neural Network，CNN）应用广泛，能够对图像、二维网格类型的欧氏空间数据进行局部特征提取，但是对于交通流、瓦斯这类非欧氏空间数据，CNN难以进行复杂的特征映射，因此无法对这类数据进行空间相关性捕捉。近年来，图神经网络（Graph Neural Network，GNN）以强大的图特征捕捉能力，成为处理复杂非结构化数据的重要方法之一[13]。由于图上平移不变性的缺失，难以将卷积直接应用在节点域，而GCN在GNN基础上引入卷积定义，通过学习到的图卷积核，在谱空间对输入信号进行图卷积[14]。目前，GCN模型已成功应用在图像分类[15]、无监督学习[16]等领域。综合考虑瓦斯浓度序列数据的时空特性，本文采用GCN对瓦斯浓度进行空间依赖建模。
在完成瓦斯空间节点图构建之后，通过GCN实现节点间的信息交互，以对节点间的空间依赖关系进行建模。
$$ {H^{\left( {l + 1} \right)}} = \sigma \left( {\tilde {\boldsymbol{D}} - \frac{1}{2}\tilde {\boldsymbol{A}} \tilde {\boldsymbol{D}} - \frac{1}{2}{H^l}{\theta ^l}} \right) $$ (7) 式中：$ {H^{(l + 1)}} $和$ {H^{(l)}} $分别为第$l + 1$层和第$l$层的输出；$ \sigma $为模型的Sigmoid激活函数；$\tilde {\boldsymbol{D}} $为归一化后的度矩阵；$\tilde {\boldsymbol{A}} $为特征矩阵；$ {\theta ^{(l)}} $为第$ l $层的参数集合。
通过反向传播方式进行整体训练，并选择2层GCN网络进行前向传播。经前向计算，将图表示的数据投影到谱域进行空间特征学习，如图6所示。在瓦斯浓度监测网络中，假设节点1是煤矿中心区域的监测点，GCN模型先得到中心节点与附近监测节点间的距离和权重，再构建节点分布拓扑图，对输入的图信号在谱域上进行投影，对瓦斯空间分布网络的拓扑结构和瓦斯节点属性进行计算，最终得到瓦斯浓度的空间特征，即捕捉到瓦斯浓度的空间依赖关系。
2.4 序列到序列模型
2.4.1 基于GRU的时间依赖建模
采用GRU对瓦斯浓度数据进行时间依赖建模。GRU模型主要由重置门和更新门组成，如图7所示，其中${{\textit{n}}_t}$为更新门的输出，${r_t}$为重置门的输出，${h_t}$为t时刻的隐藏状态，$ {h_{t - 1}} $为t−1时刻的隐藏状态。重置门通过对上一时刻的隐藏状态与当前时刻的输入数据进行拼接，再通过tanh激活函数将数据进行缩放，实现对上一时刻隐藏状态的选择性记忆；更新门通过对上一时刻和当前时刻的计算，决定信息的流入量。
具体计算过程如下：
$$ {{\textit{n}}_t} = \sigma \left( {{{\boldsymbol{W}}^{\left( {\textit{n}} \right)}}X_t' + {{\boldsymbol{U}}^{\left( {\textit{n}} \right)}}{h_{t - 1}}} \right) $$ (8) $$ {r_t} = \sigma \left( {{{\boldsymbol{W}}^{\left( r \right)}}X_t' + {{\boldsymbol{U}}^{\left( r \right)}}{h_{t - 1}}} \right) $$ (9) $$ {h}_{t}'=\mathrm{tanh}({\boldsymbol{W}}{X}_{t}'+{r}_{t}\odot {\boldsymbol{U}}{h}_{t-1}) $$ (10) $$ {h}_{t}={h}_{t}'\odot\left(1-{{\textit{n}}}_{{t}}\right)+{{\textit{n}}}_{t}\odot{h}_{t-1} $$ (11) 式中：$X_t'$为经多层图卷积计算后的图信号；$h_t'$为基于重置门计算的隐藏状态；${{\boldsymbol{W}}^{(n)}}$、$ {{\boldsymbol{W}}^{(r)}} $、${{\boldsymbol{U}}^{(n)}}$、$ {{\boldsymbol{U}}^{(r)}} $、$ {\boldsymbol{W}} $、${\boldsymbol{ U}} $均为可学习参数，以权值矩阵的形式参与计算；⊙为哈达玛乘积。
2.4.2 序列编解码
序列到序列模型在文本摘要[17]、回归预测、语义问答[18]和语音识别[19]等领域得到了广泛应用，是近年来序列数据处理领域产生重大影响的模型之一[20]。
本文采用序列到序列模型进行长期预测，并引入自动编码器生成预测结果。序列到序列模型的内部结构如图8所示。编码器内部主要由GCN和GRU 2种网络结构单元组成，通过自动编码器将瓦斯浓度数据映射为隐藏空间表征K，将K输入解码器进行数据重构，获得最终时空预测输出。通过数据压缩，对输入的数据进行降维去噪处理，提高训练效率。
假设从初始时刻t到最终时刻T输入到模型的瓦斯图信号为$\{ {X_t},{X_{t+1}},\cdots ,{X_T}\}$，将该图信号经过2层神经网络，得到经图卷积计算后的瓦斯图信号$\{ X_t', X_{t+1}',\cdots ,X_T'\}$，再经过GRU单元，得到瓦斯浓度序列的时间依赖关系。
$$ {Y_t} = {f_{{\rm{GRU}}}}({h_{t - 1}},X_t') $$ (12) 式中：$ {Y_t} $为编码器的输出，即瓦斯浓度序列的时间依赖关系；$ {f_{{\rm{GRU}}}}( \cdot ) $为GRU单元内部的计算过程。
通过编码器结构计算得到隐藏空间表征K，将空间表征序列作为站点编码向量，用来初始化解码GRU。在解码器中，GRU单元将每一时刻得到的隐藏状态输入到全连接层，最终得到每一步的预测结果。
3. 试验与结果分析
3.1 瓦斯浓度数据集
本试验所用的瓦斯浓度数据来自某煤矿2018年3−5月的工作面实时监测记录。数据的平均采集间隔为5 min，共采集26 481条数据。由于历史数据采集的时间窗口长度与瓦斯浓度数据存在时间关联性，不同的窗口大小影响着预测效果，所以选择大小分别为6，12，20，50的滑动窗口进行模型训练。当窗口大小为20时，损失函数收敛速度最快，因此设置数据采集窗口大小为20。试验设置采集步长为1，对未来6个时间步的瓦斯浓度进行预测，并按时间顺序以6∶2∶2的比例将数据集划分为训练集、验证集和测试集。
在采集过程中，大量不可预测因素（如灰尘、水汽、人为因素等）易造成数据不完整或不准确，导致瓦斯浓度数据出现异常值或缺失值。因此根据国家标准[11]，对采集到的瓦斯数据设定有效值范围，以排除异常值，所有不在有效值范围内的节点监测数据将被认为是异常值，在后续处理中视为缺失值。
本文选择3σ法进行异常值检测，当数值与平均值的差值超过3倍标准差时，将该数值判定为异常值。选择线性插值法进行缺失值填充，通过对包含缺失值的一段序列数值进行线性拟合，计算缺失位置的理论数值。
$$ k(a) = k({a_0}) + \frac{{k({a_1}) - k({a_0})}}{{{a_1} - {a_0}}}(a - {a_0}) $$ (13) 式中：$k(a)$为拟合出的线性变化函数；a为缺失值的数据样本序号；${a_0}$和${a_1}$分别为缺失值前、后的有效数据样本序号。
3.2 试验设置
模型采用Python语言在PyTorch深度学习框架下实现。训练的超参数设定：批数量大小为64，学习率为0.003。采用Adam优化器，对模型进行最多100个轮次的训练。为避免过拟合，选择早停策略，使模型直接使用训练中的最优参数。训练损失采用MSE函数计算。
在其他模型参数不变的情况下，分别使用不同的网络层数b和隐藏单元数c进行测试。测试结果如图9所示，其中b∈{1，2，3，4}，c∈{16，32，64，128}。可看出随着隐藏单元数增加，预测误差先减小后增大，当隐藏单元数为64时，误差最小；当网络层数为2时，预测误差最小。故本试验设置隐藏单元数为64，网络层数为2。
3.3 模型对比和结果分析
本文模型的损失函数曲线如图10所示。可看出损失曲线在不到20个训练轮次中开始收敛并趋于平稳，整体波动较小。
为验证本文模型的准确性，选用基于时间序列的基线模型与本文模型进行对比。历史平均（Historical Average，HA）模型通过历史时段内的平均瓦斯浓度信息进行预测；移动平均自回归模型（Autoregressive Integrated Moving Average Model，ARIMA）是一种经典的时间序列预测模型，通过对原始序列数据进行识别，将识别的数据与当前序列数据同时进行分析，预测未来值；支持向量回归（Support Vector Regression，SVR）模型通过在特征空间寻找最优超平面，使得特征空间中所有样本到该平面的距离最小。
选取HA模型、SVR模型与本文模型进行对比，不同模型在3月份的预测结果如图11所示。可看出本文模型预测结果与实际数据的拟合度优于HA模型和SVR模型。本文模型能够较为准确地预测瓦斯浓度的总体变化趋势，但是对于单位时间内瓦斯浓度波动较大的部分（如局部极大值处），模型难以进行精确预测，出现该现象的主要原因是GCN模型在傅里叶域中定义了一个平滑的滤波器，该滤波器通过移动计算分布点的空间特征，不断对图信号进行平滑操作，使得图信号缺少一定的可区分性。
为更加客观地了解各模型的预测精度，验证本文模型的性能，采用均方根误差、平均绝对误差、准确率作为预测精度评价指标，采用判定系数R2衡量模型的预测能力。各模型的性能指标对比见表1，其中*表示数值过小，故忽略不计。可看出HA模型的判定系数R2略高于本文模型，且均方根误差和平均绝对误差较本文模型的高。SVR模型的均方根误差和平均绝对误差较HA模型明显增大，相较其他模型，SVR的判定系数R2明显减小，说明SVR模型在实验过程中出现过拟合现象，导致模型鲁棒性降低，瓦斯浓度预测精度下降。ARIMA模型的均方根误差和平均绝对误差均高于本文模型，准确率低于本文模型。本文模型的均方根误差较HA模型、SVR模型、ARIMA模型分别降低了0.5%，71.4%，37.9%，平均绝对误差较HA模型、SVR模型、ARIMA模型分别降低了10.5%，82.4%，82.4%准确率较HA模型、SVR模型、ARIMA模型分别提高了0.06%，17.7%，13.8%，表明本文模型具有较强的鲁棒性，且泛化性能较好。
表 1 各模型性能指标Table 1. Performance indexes of each model模型 均方根误差 平均绝对误差 准确率/% R2 HA 0.039 3 0.023 2 79.89 0.794 7 SVR 0.067 0 0.061 7 65.76 0.405 1 ARIMA 0.053 9 0.038 3 68.91 * 本文模型 0.039 1 0.021 0 79.94 0.698 1 本文模型较HA模型、SVR模型、ARIMA模型更能关注到前序重要特征的影响。这主要是由于GRU的2个门关注数据的时间特征，而HA，ARIMA等方法无法对长期记忆进行选择，导致数据冗余，因此对复杂的、非平稳的时间序列数据，难以保证预测精度。作为循环神经网络的改进模型，GRU在保留门控功能的基础上，减少训练参数，在一定程度上提高了模型训练效率，降低了训练时长。
4. 结论
1） 瓦斯浓度预测是矿井瓦斯灾害防治的重要依据，研究瓦斯的时空分布预测问题，以GCN为基础建立瓦斯浓度的时空预测模型，具有更高的预测准确性。本文分析了瓦斯浓度的时间和空间特性，通过GCN对瓦斯浓度进行空间分布建模，获取瓦斯的空间分布特征，以GRU进行瓦斯浓度时间变化建模，提出一种预测煤矿井下瓦斯时空分布的深度神经网络模型。
2） 本文模型预测结果与实际数据的拟合度优于HA模型和SVR模型。本文模型较HA模型、SVR模型、ARIMA模型更能关注到前序重要特征的影响。
3） 本文模型的均方根误差较HA模型、SVR模型、ARIMA模型分别降低了0.5%，71.4%，37.9%，平均绝对误差分别降低了10.5%，82.4%，82.4%，准确率分别提高了0.06%，17.7%，13.8%，表明本文模型具有较强的鲁棒性，且泛化性能较好。
-
![]()
图 2 第1和第2判别函数分组结果
Figure 2. Grouping results of the first and second discriminant functions
表 1 训练样本
Table 1 Training sample
序号 判别指标含量/（mg·L−1） 类别 Ca2+ Mg2+ Na++ K+ HCO3− Cl− SO42− M 1 95.13 31.92 1 941.1 289.58 131.16 4 230.10 6 739.90 砂岩水 2 225.70 81.84 2 171.0 266.83 109.90 4 836.41 7 698.30 砂岩水 3 133.07 55.43 2 206.8 300.67 140.03 4 428.37 7 286.80 砂岩水 4 30.84 24.19 1 463.6 406.93 287.14 2 454.08 4 680.81 砂岩水 5 40.30 23.45 1 485.0 405.42 297.78 2 749.76 5 019.32 砂岩水 6 23.49 11.54 1 843.9 524.47 235.74 3 221.08 5 868.13 砂岩水 7 60.72 25.04 1 665.9 302.45 251.70 3 094.96 5 423.65 砂岩水 8 73.75 31.46 1 613.5 328.19 274.74 3 117.44 5 453.08 砂岩水 9 38.49 13.52 1 724.6 424.72 265.88 3 156.29 5 654.40 砂岩水 10 58.76 23.35 1 757.2 470.46 280.06 2 952.35 5 560.64 砂岩水 11 66.86 31.95 1 649.0 276.94 274.74 2 903.79 5 243.07 砂岩水 12 62.81 18.43 1 920.0 437.09 210.93 3 603.03 6 260.54 砂岩水 13 200.40 81.44 736.4 350.78 177.25 1 964.43 3 517.73 太灰水 14 173.15 28.69 688.0 0 202.06 1 469.72 2 605.20 太灰水 15 99.40 36.95 1 129.6 400.90 191.43 2 065.29 3 938.47 太灰水 16 136.27 34.03 686.0 0 198.52 1 421.69 2 515.40 太灰水 17 133.70 63.72 1 350.0 316.97 251.34 2 925.84 5 050.99 太灰水 18 229.26 0 1 398.0 255.46 246.38 2 065.29 4 347.13 太灰水 19 138.28 0 1 503.0 133.43 249.92 2 276.62 4 438.26 太灰水 20 192.38 85.09 919.9 336.47 203.84 2 046.08 3 794.84 太灰水 21 214.43 88.73 913.6 332.89 200.29 2 113.32 3 869.80 太灰水 22 320.64 106.96 489.6 304.25 219.79 1 719.47 3 168.75 太灰水 23 194.31 77.34 669.0 193.29 207.38 1 600.84 2 948.84 太灰水 24 184.37 80.71 906.0 322.15 209.16 1 988.44 3 694.47 太灰水 25 183.57 79.49 927.1 325.73 205.61 2 094.11 3 823.36 太灰水 26 186.21 79.80 901.2 318.57 207.38 1 969.52 3 667.63 太灰水 27 215.63 90.43 823.1 329.31 210.93 1 911.59 3 585.55 太灰水 28 199.60 86.79 929.1 322.15 209.16 2 074.90 3 830.01 太灰水 29 157.11 74.87 958.6 300.67 200.29 2 209.38 3 905.09 太灰水 30 198.36 85.44 896.2 332.89 203.84 1 998.62 3 554.19 太灰水 31 201.59 86.43 894.8 332.89 203.84 2 008.33 3 566.52 太灰水 32 202.40 85.94 889.7 340.05 200.29 2 037.43 3 590.25 太灰水 33 201.59 87.65 878.0 332.89 202.06 1 979.22 3 520.03 太灰水 34 203.21 87.41 889.2 340.05 198.52 2 027.73 3 580.54 太灰水 35 204.43 85.94 891.3 340.05 198.52 2 047.13 3 602.18 太灰水 36 202.40 85.44 888.8 332.89 202.06 1 988.92 3 539.53 太灰水 37 202.40 87.90 906.6 336.47 202.06 2 076.24 3 647.83 太灰水 38 206.45 85.94 886.6 336.47 202.06 2 095.64 3 649.20 太灰水 39 202.40 87.90 903.9 336.47 200.29 2 047.13 3 614.31 太灰水 40 331.94 116.26 428.05 345.42 215.36 1 474.71 2 742.29 混合水 41 246.93 168.43 408.50 338.88 218.90 1 455.31 2 671.18 混合水 42 332.75 116.38 429.60 309.69 220.68 1 541.66 2 799.16 混合水 43 344.09 117.36 381.0 311.18 221.56 1 523.22 2 746.15 混合水 44 331.94 114.42 434.30 311.18 222.27 1 538.75 2 800.68 混合水 45 238.03 102.39 1 068.40 287.36 176.36 2 580.75 4 315.95 混合水 
下载: 导出CSV
表 2 典型判别函数系数项和常数项
Table 2 Coefficient term and constant term of typical discriminant function
变量 第1判别函数 第2判别函数 Ca2+ −0.013 10 0.017 2 Mg2+ −0.018 80 0.006 9 Na++K+ −0.003 40 0.002 1 HCO3− 0.001 40 0.003 4 Cl− 0.020 80 0.041 7 SO42− 0.004 70 0.008 6 M −0.000 37 −0.005 2 常数项 −7.175 40 −14.165 5
下载: 导出CSV
表 3 判别函数方差
Table 3 Variance of discriminant function
判别函数 特征值 方差贡献率/% 累积方差贡献率/% 典型相关系数 第1判别 8.514 92.8 92.8 0.946 第2判别 0.660 7.2 100.0 0.631
下载: 导出CSV
表 4 典型判别函数在各类别的中心值
Table 4 Center value of typical discriminant function in each category
类别 判别函数 第1判别函数 第2判别函数 砂岩水 4.509 0.343 太灰水 −1.230 −0.542 混合水 −3.481 1.751
下载: 导出CSV
表 5 回代估计结果
Table 5 Back substitution estimation results
序号 第1判
别函数值距砂岩水
中心值
距离距太灰水
中心值
距离距混合水
中心值
距离回代
判别实际
类别序号 第1判
别函数值距砂岩水
中心值
距离距太灰水
中心值
距离距混合水
中心值
距离回代
判别实际
类别1 5.205 0.696 6.435 8.686 砂岩水 砂岩水 24 −1.269 5.778 0.039 2.212 太灰水 太灰水 2 3.831 0.678 5.061 7.312 砂岩水 砂岩水 25 −0.921 5.430 0.309 2.560 太灰水 太灰水 3 4.317 0.192 5.547 7.798 砂岩水 砂岩水 26 −1.382 5.891 0.152 2.099 太灰水 太灰水 4 3.536 0.973 4.766 7.017 砂岩水 砂岩水 27 −1.863 6.372 0.633 1.618 太灰水 太灰水 5 4.849 0.340 6.079 8.330 砂岩水 砂岩水 28 −1.301 5.810 0.071 2.180 太灰水 太灰水 6 4.902 0.393 6.132 8.383 砂岩水 砂岩水 29 −0.222 4.731 1.008 3.259 太灰水 太灰水 7 4.331 0.178 5.561 7.812 砂岩水 砂岩水 30 −1.504 6.013 0.274 1.977 太灰水 太灰水 8 4.824 0.315 6.054 8.305 砂岩水 砂岩水 31 −1.519 6.028 0.289 1.962 太灰水 太灰水 9 5.319 0.810 6.549 8.800 砂岩水 砂岩水 32 −1.438 5.947 0.208 2.043 太灰水 太灰水 10 4.188 0.321 5.418 7.669 砂岩水 砂岩水 33 −1.644 6.153 0.414 1.837 太灰水 太灰水 11 3.779 0.730 5.009 7.260 砂岩水 砂岩水 34 −1.553 6.062 0.323 1.928 太灰水 太灰水 12 5.027 0.518 6.257 8.508 砂岩水 砂岩水 35 −1.465 5.974 0.235 2.016 太灰水 太灰水 13 −1.598 6.107 0.368 1.883 太灰水 太灰水 36 −1.610 6.119 0.380 1.871 太灰水 太灰水 14 −2.083 6.592 0.853 1.398 太灰水 太灰水 37 −1.336 5.845 0.106 2.145 太灰水 太灰水 15 −0.055 4.564 1.175 3.426 太灰水 太灰水 38 −1.195 5.704 0.035 2.286 太灰水 太灰水 16 −1.961 6.470 0.731 1.520 太灰水 太灰水 39 −1.490 5.999 0.260 1.991 太灰水 太灰水 17 3.046 1.463 4.276 6.527 砂岩水 太灰水 40 −4.204 8.713 2.974 0.723 混合水 混合水 18 −1.177 5.686 0.053 2.304 太灰水 太灰水 41 −4.007 8.516 2.777 0.526 混合水 混合水 19 0.533 3.976 1.763 4.014 太灰水 太灰水 42 −3.867 8.376 2.637 0.386 混合水 混合水 20 −1.357 5.866 0.127 2.124 太灰水 太灰水 43 −3.919 8.428 2.689 0.438 混合水 混合水 21 −1.481 5.990 0.251 2.000 太灰水 太灰水 44 −3.814 8.323 2.584 0.333 混合水 混合水 22 −3.050 7.559 1.820 0.431 混合水 太灰水 45 −1.076 5.585 0.154 2.405 太灰水 混合水 23 −2.330 6.839 1.100 1.151 太灰水 太灰水
下载: 导出CSV
表 6 煤层底板突水水源判别分析模型指标数值
Table 6 Indicator value of mine floor water inrush source discrminant model
序号 判别指标含量/（mg·L−1） 实际类别 Ca2+ Mg2+ Na++ K+ HCO3− Cl− SO42− M 1 86.43 35.36 1706.6 373.70 214.47 3333.65 5768.76 砂岩水 2 76.76 31.04 1681.9 320.31 235.740 3085.40 5447.14 砂岩水 3 88.03 36.77 1536.0 347.00 245.67 3144.77 5412.92 砂岩水 4 206.45 85.94 900.6 336.47 200.29 2076.24 3647.26 太灰水 5 207.26 82.50 900.8 325.73 198.52 2037.43 3596.87 太灰水 6 209.69 83.97 888.1 336.47 200.29 2018.03 3572.50 太灰水 7 208.48 87.16 883.1 329.31 202.06 2008.33 3561.56 太灰水 8 210.5 85.94 868.8 336.47 200.29 2056.84 3594.60 太灰水 9 327.89 115.40 403.5 299.27 217.66 1534.87 2752.40 混合水 10 323.85 112.08 404.0 284.38 219.79 1533.90 2743.50 混合水
下载: 导出CSV
表 7 各水样的Fisher判别结果
Table 7 Fisher discrimination results of each water sample
序号 第1判别
函数值距砂岩水
中心值距离距太灰水
中心值距离距混合水
中心值距离判别结果 1 3.914 0.595 5.144 7.395 砂岩水 2 3.523 0.986 4.753 7.004 砂岩水 3 4.285 0.224 5.515 7.766 砂岩水 4 −1.422 5.931 0.192 2.059 太灰水 5 −1.584 6.093 0.354 1.897 太灰水 6 −1.632 6.141 0.402 1.849 太灰水 7 −1.674 6.183 0.444 1.807 太灰水 8 −1.441 5.950 0.211 2.040 太灰水 9 −3.830 8.339 2.600 0.349 混合水 10 −3.694 8.203 2.464 0.213 混合水
下载: 导出CSV
-
[1] 刘伟,刘晨君. 改革开放四十年来煤炭行业安全发展之路[J]. 煤炭经济研究,2018,38(11):34-42. LIU Wei,LIU Chenjun. The road to safe development of the coal industry in the past 40 years of reform and opening up[J]. Coal Economic Research,2018,38(11):34-42.
[2] MA Dan,DUAN Hongyu,CAI Xin,et al. A global optimization-based method for the prediction of water inrush hazard from mining floor[J]. Water,2018,10(11):1618. DOI: 10.3390/w10111618
[3] WANG Jing,LI Shucai,LI Liping,et al. Attribute recognition model for risk assessment of water inrush[J]. Bulletin of Engineering Geology & the Environment,2019,78(2):1057-1071.
[4] 彭苏萍. 我国煤矿安全高效开采地质保障系统研究现状及展望[J]. 煤炭学报,2020,45(7):2331-2345. PENG Suping. Current status and prospects of research on geological assurance system for coal mine safe and high efficient mining[J]. Journal of China Coal Society,2020,45(7):2331-2345.
[5] 汪洋,左文喆,王斌海,等. 矿井突水水源判别方法研究进展[J]. 现代矿业,2018,34(1):69-73. DOI: 10.3969/j.issn.1674-6082.2018.01.012 WANG Yang,ZUO Wenzhe,WANG Binhai,et al. Study progress of discriminanat method of the sources of mine water inrush[J]. Modern Mining,2018,34(1):69-73. DOI: 10.3969/j.issn.1674-6082.2018.01.012
[6] 李燕,徐志敏,刘勇. 矿井突水水源判别方法概述[J]. 煤炭技术,2010,29(11):87-89. LI Yan,XU Zhimin,LIU Yong. Summary on methods of distinguishing sources of mine water-invasion[J]. Coal Technology,2010,29(11):87-89.
[7] 穆文平. 北阳庄矿煤层底板断层突水机理与岩溶水疏降水量预测[D]. 北京: 中国矿业大学（北京）, 2018. MU Wenping. Mechanism of water inrush on faults of coal seam floor and prediction of dewatering rate from karst aquifers in Beiyangzhuang Mine[D]. Beijing: China University of Mining and Technology（Beijing）, 2018.
[8] 孙林华,桂和荣,陈松. 深层地下水稀土元素无机形态及其对稀土特征的影响−以皖北任楼矿煤系含水层为例[J]. 煤田地质与勘探,2011,39(3):38-43. DOI: 10.3969/j.issn.1001-1986.2011.03.008 SUN Linhua,GUI Herong,CHEN Song. Inorganic speciation of deep groundwater and its effect on the characteristics of rare earth elements:with aquifers in coal bearing masrues in Renlou Coal Mine in the north of Anhui province as example[J]. Coal Geology & Exploation,2011,39(3):38-43. DOI: 10.3969/j.issn.1001-1986.2011.03.008
[9] 乔元栋,孟召平,张村,等. 复杂构造井田含水层特征及其水力联系辨识[J]. 煤炭学报,2021,46(12):4010-4020. QIAO Yuandong,MENG Zhaoping,ZHANG Cun,et al. Characteristics of mine aquifer with complex structure and identification of its hydraulic connection[J]. Journal of China Coal Society,2021,46(12):4010-4020.
[10] 曾一凡, 梅傲霜, 武强, 等. 基于水化学场机器学习分析与水动力场反向示踪模拟耦合的矿井涌（突）水水源综合判识技术[J/OL]. 煤炭学报: 1-14[2022-03-22]. DOI: 10.13225/j.cnki.jccs.2021.1979. ZENG Yifan, MEI Aoshuang, WU Qiang, et al. Source discrimination of mine water inflow or inrush using a hydrochemical field machine learing analysis and hydrodynamic field reverse tracer simulation coupling technique[J/OL]. Journal of China Coal Society: 1-14[2022-03-22]. DOI: 10.13225/j.cnki.jccs.2021.1979.
[11] 翟明娟. 距离判别分析及其评价[J]. 长治学院学报,2012,29(2):34-36,55. DOI: 10.3969/j.issn.1673-2014.2012.02.011 ZHAI Mingjuan. Distance discretion analysis and evaluation[J]. Journal of Changzhi University,2012,29(2):34-36,55. DOI: 10.3969/j.issn.1673-2014.2012.02.011
[12] 徐星,李喆,曾珠. 煤层底板突水危险性AHP−Fisher判别分析模型[J]. 煤炭技术,2018,37(7):153-155. XU Xing,LI Zhe,ZENG Zhu. AHP-Fisher discrimination model of water inrush risk from coal seam floor[J]. Coal Technology,2018,37(7):153-155.
[13] 黄利文. 基于变量择优的Fisher逐步判别分析方法[J]. 系统科学与数学,2021,41(8):2338-2348. HUANG Liwen. Fisher stepwise discriminant analysis method based on selecting optimal variables[J]. Journal of Systems Science and Mathematical Sciences,2021,41(8):2338-2348.
[14] 温廷新,于凤娥. 基于KPCA−Fisher判别分析的煤炭自燃预测研究[J]. 矿业安全与环保,2018,45(2):49-53,58. DOI: 10.3969/j.issn.1008-4495.2018.02.011 WEN Tingxin,YU Feng'e. Research on prediction of coal spontaneous combustion based on KPCA-Fisher discriminant analysis[J]. Mining Safety & Environmental Protection,2018,45(2):49-53,58. DOI: 10.3969/j.issn.1008-4495.2018.02.011
-
期刊类型引用(5)
1. 李晓梅，黄建勇，张泽治. 基于改进权函数距离的机器人运动偏差补偿算法设计. 吉林大学学报(信息科学版). 2024(01): 86-92 . 
百度学术
2. 郑林超. 掘进机截齿截割的受力分析及仿真. 机械管理开发. 2024(02): 79-80+83 . 百度学术
3. 王豪，栗继祖. 掘进机二运防脱道断电装置的设计及应用分析. 机械管理开发. 2023(01): 134-135 . 百度学术
4. 赵书山. 镐型截齿不同角度对掘进机截割性能的影响研究. 机械管理开发. 2023(06): 28-29+32 . 百度学术
5. 韩永亮. 智能掘进技术在任家庄煤矿的应用及创新. 工矿自动化. 2022(S2): 30-32+41 . 本站查看
其他类型引用(0)
计量
- 文章访问数: 253
- HTML全文浏览量: 82
- PDF下载量: 32
- 被引次数: 5
