A temperature prediction model for coal spontaneous combustion based on PSO-SRU deep artificial neural networks
-
摘要: 针对传统煤自燃温度预测模型泛化能力不强、鲁棒性较差的问题，提出了一种基于改进粒子群(PSO)优化简单循环单元(SRU)的煤自燃温度预测模型(PSO−SRU模型)。首先，对煤自燃程序升温实验中采集的气体浓度数据进行预处理，选取与煤温相关性较强的O2，CO，CO2，CH4，C2H4作为煤温预测指标，并将预测指标划分为训练集和测试集；其次，构建SRU预测模型拟合训练集中煤自燃温度与气体指标间非线性规律，将平均绝对误差(MAE)作为适应度函数，利用改进的PSO算法优化SRU预测模型参数；最后，将测试集数据输入参数最优的SRU预测模型，利用SRU计算得到煤自燃温度预测值。实验结果表明：通过指标择优和参数寻优后，PSO−SRU模型在测试集上的MAE相较于基于支持向量回归(SVR)、随机森林（RF）和反向传播（BP）的煤自燃温度预测模型分别降低了12.58，7.65，5.91 ℃，表明PSO−SRU模型在一定程度上提高了预测精度；均方根误差（RMSE）分别降低了22.65，17.45，8.94 ℃，PSO−SRU模型在训练集和测试集上的决定系数(R2)仅相差0.03，表明PSO−SRU模型具有良好的泛化性和鲁棒性。Abstract: Traditional temperature prediction models for coal spontaneous combustion typically have low generality and robustness. This paper improves them by proposing a coal spontaneous combustion temperature prediction model based on particle swarm optimization and simple recurrent unit(PSO-SRU). It firstly pre-processes the gas concentration data collected from temperature programmed oxidation tests, selects the concentration data of O2, CO, CO2, CH4, C2H4 that highly relate to the coal temperature as the prediction indicators, and further separates the indicators into training and testing data sets. Then, a SRU based prediction model over the training data set is trained to learn the nonlinear relationship between the coal spontaneous combustion temperature and the indicators. Mean absolute error(MAE) forms the fitness function and PSO algorithms are involved to optimize the SRU prediction model's parameters. Finally, the PSO-SRU model with optimized parameters are applied over the testing data set to predict the coal spontaneous combustion temperature. Experiments show the PSO-SRU model can improve the prediction accuracy, as the model's MAE and root mean square error(RMSE), comparing with those generated by support vector regression(SVR), random forest(RF), and back propagation(BP), decreases by 12.58, 7.65, 5.91 ℃, and 22.65, 17.45, 8.94 ℃ respectively. The PSO-SRU model also demonstrates a good generality and robustness, as the difference of determination coefficient (R2) of the model over the training and testing data sets is only 0.03.
-
0. 引言
我国对煤炭的需求短期内不会发生根本性变化，煤炭在能源结构中仍将长期处于主体地位[1-2]。长期大量开采导致采空区遗煤增多，煤自燃灾害频发[3-4]。煤自燃是一个复杂的动态氧化过程，一旦发生会造成巨大资源浪费和环境污染[5-6]。因此，开展煤自燃的温度预测研究对防控煤自燃灾害发生具有重要意义[7]。
近年来，学者们围绕煤自燃温度预测相关问题提出了多种方法，主要有测温法、自燃实验预测法和气体分析法等[8-10]。其中气体分析法根据煤体氧化升温时不同温度状态会释放不同浓度的气体产物这一规律，通过监测煤自燃过程中气体指标浓度预测煤自燃温度[11-13]。气体分析法因规律性强、灵敏度高而被广泛应用[14-15]。邓军等[16]建立了支持向量回归(Support Vector Regression, SVR)的煤自燃预测模型，并采用粒子群(Particle Swarm Optimization, PSO)选取最佳的核函数和惩罚因子，提高了预测精度，但该模型对参数选取敏感，易陷入局部最优。刘宝等[17]提出了基于相关向量机(Relevance Vector Machine, RVM)的煤自燃预测方法，简化了模型参数选取，泛化性较好，但预测精度有待提高。昝军才等[18]采用反向传播(Back Propagation, BP)神经网络学习煤温与气体指标的非线性映射关系，实现煤自燃温度预测，但该方法易出现过拟合现象，泛化能力较差。郑学召等[19]建立了基于随机森林(Random Forest, RF)的煤自燃温度预测模型，优化了决策树深度和数量，参数优化简单，泛化性好，但预测精度和模型鲁棒性有待提高。
综上可知，现有的煤自燃温度预测方法存在预测模型泛化性不强、鲁棒性较差的问题。针对上述问题，本文考虑真实煤自燃监测数据前后具有强关联性，提出了基于PSO算法优化简单循环单元(Simple Recurrent Units, SRU)的煤自燃温度预测模型−PSO−SRU模型。将煤自燃多特征数据输入SRU预测模型，利用SRU神经单元挖掘特征数据间的非线性关系，根据挖掘到的规律和当前煤自燃特征的输入预测煤自燃温度；通过改进的PSO算法优化SRU预测模型参数，提高了模型的预测精度和鲁棒性。 该模型可为利用深度神经网络预测煤自燃温度提供解决思路。
1. 理论基础
1.1 SRU
SRU是循环神经网络(Recurrent Neural Network, RNN)的优秀变体[20]，与长短期记忆网络(Long Short Term Memory, LSTM)和门循环单元(Gate Recurrent Unit, GRU)等循环神经单元的结构一样，都是基于“门控”结构构建。SRU简化了单元结构，使RNN更加轻便，其内部结构和网络结构如图1所示。
![]()
图 1 SRU内部结构和网络结构Figure 1. Simple recurrent units(SRU) interior structure and network structureSRU结构计算过程如下：
$$ {{{{f}}}_{{{t}}}} = {\sigma_{\rm{F}}}\left({{{\boldsymbol{W}}_{\rm{F}}}{{\boldsymbol{X}}_{{{t}}}} + {{\boldsymbol{U}}_{\rm{F}}} \odot {{{{c}}}_{{{{t}} - 1}}} + {{{\boldsymbol{b}}}_{\rm{F}}}} \right)$$ (1) $$ {{{{c}}}_{{{t}}}} = {{{{f}}}_{{{t}}}} \odot {{{{c}}}_{{{{t}} - 1}}} + \left( {1 - {{{{f}}}_{{{t}}}}} \right) \odot \left( {{\boldsymbol{W}}{{\boldsymbol{X}}_{{{t}}}}} \right) $$ (2) $$ {{{{r}}}_{{{t}}}} = {\sigma _{\rm{R}}}\left( {{{\boldsymbol{W}}_{\rm{R}}}{{\boldsymbol{X}}_{{{t}}}} + {{\boldsymbol{U}}_{\rm{R}}} \odot {{{{c}}}_{{{{t}} - 1}}} + {{{\boldsymbol{b}}}_{\rm{R}}}} \right) $$ (3) $$ {{{{h}}}_{{{t}}}} = {{{{r}}}_{{{t}}}} \odot {{{{c}}}_{{{t}}}} + \left( {1 - {{{{r}}}_{{{t}}}}} \right) \odot {{\boldsymbol{X}}_{{{t}}}} $$ (4) 式中：ft为遗忘门，用来决定单元状态需要丢弃过去哪些信息，t为时刻；σF为遗忘门激活函数；WF为遗忘门细胞状态的参数矩阵；Xt为输入数据；UF , bF, UR, bR为需要在训练过程中学习的参数向量；⊙为点乘运算符号；ct−1 , ct为神经单元状态；W 为细胞状态参数矩阵；rt为重置门，控制前一单元状态有多少信息被写入到当前细胞ht中；σR为重置门激活函数；WR为重置门细胞状态的参数矩阵。
SRU网络结构由轻量循环（式（1）和式（2））和高速网络（式（3）和式（4））2个重要部分组成。轻量循环采用UF⊙ct−1 的计算方式，降低单元与单元间的依赖程度，提高运算速度。高速网络采用(1− rt)⊙Xt跳跃连接的方式，使梯度能够直接向前传播，缓解梯度消失和梯度爆炸问题。
1.2 PSO算法
PSO算法从拟态生物学的角度出发，将鸟群在空间内的觅食过程、状态及行为抽象为无质量粒子在多维空间内搜索满足最优适应度值的过程、状态和行为。其核心思想是信息共享，通过粒子群中记录的全局最优适应度和个体最优适应度更新粒子的速度和位置，计算自身的适应度，将信息分享给群体，各粒子调整自身的速度和位置，朝最优适应度粒子周围聚集，达到搜索终止条件的最优适应度粒子的每个维度值即为全局最优解[21]。
在D维搜索空间下，PSO算法的更新公式为
$$ v_i^{q + 1} = {\bar \omega _{}}v_i^q + {\alpha _1}\beta _1^q\left( {p_{{\rm{best}}i}^q - \lambda _i^q} \right) + {\alpha _2}\beta _2^q\left( {g_{{\rm{best}}i}^q - \lambda _i^q} \right) $$ (5) $$ \lambda _i^{q + 1} = \lambda _i^q + v_i^q $$ (6) 式中：νi，λi分别为第i个粒子的速度和位置，i=1，2，…，s，s为粒子数，通常取40；
$\bar \omega$ 为惯性权重，通常取0.9；q为当前迭代次数；α1，α2为学习因子，控制着第q−1次迭代的经验对第q次迭代的影响程度，通常α1=α2=0.5；β1，β2为[0,1]内的随机数；pbesti为个体粒子历史最优位置；gbesti为全局最优个体粒子位置。2. 基于PSO−SRU的煤自燃温度预测模型
2.1 PSO−SRU模型构建
PSO−SRU模型架构如图2所示。通过煤自燃程序升温实验采集到n个时刻煤自燃数据样本，令数据集合X={X1, X2, …, Xn}，Xt为第t个时刻样本，Xt=[xt1, xt2, …, xtm]，xtj为第t个时刻第j个预警指标气体浓度，t=1，2，…，n，j=1，2，…，m，m为预警气体指标个数。PSO−SRU模型构建过程如下：
![]()
图 2 PSO−SRU模型架构Figure 2. Temperature prediction model framework for coal spontaneous combustion based on particle swarm optimization and simple recurrent unit (PSO-SRU)(1) 通过实验采集煤自燃预警指标气体浓度，对数据进行清洗、归一化处理，计算气体与煤温相关性，筛选相关性气体指标并将数据集划分为训练集和测试集。
(2) 建立SRU煤自燃温度预测模型，将Xt输入SRU进行矩阵计算，通过式(1)—式(4)计算得到下一时刻煤自燃温度
$ {P'_{t + 1}} $ 。(3) 提升SRU预测模型性能，将隐藏层神经元数L1，L2，学习率ɛ及模型训练次数e作为PSO算法寻优对象，改进PSO算法，提升寻优算法性能。
(4) 用最优超参数建立SRU预测模型，将测试集输入参数最优的SRU预测模型，得到煤温预测结果。
2.1.1 数据采集及处理
实验选取某矿采集的气煤煤样进行程序升温实验。装煤总质量为1 kg，装煤高度为17.5 cm，混合煤样平均粒径为 4.18 mm，空隙率为0.48%，平均升温速率为0.3 ℃/min，供风量为120 mL/min。利用程序升温装置对煤样进行加热，并向程序升温装置送入预热空气，测定气体产物，当温度升高到400 ℃预定温度时停止加热。共得到625组样本数据，记录气煤煤样氧化升温过程中产生的多种气体(如CO，CO2，CH4，C2H4，O2)浓度和温度等。625组原始数据中包含缺失值和异常值的数据，会使建模过程陷入混乱，导致输出结果不可靠，通过剔除异常值并对缺失值进行插补，得到601组数据，采用MinMaxScaler方法将数值映射到[0,1]之间。
2.1.2 预警气体指标分析
不同气体指标与煤温随时间变化关系如图3所示。从图3(a)可看出，因为煤氧复合反应的发生，实验环境中的O2体积分数呈现随煤温升高而逐渐减小的规律。从图3(b)、图3(c)可看出，从煤温开始发生变化时，CO和CH4体积分数呈现随煤温升高而不断增大的趋势，并且与煤自燃温度变化有明显的强相关性关系，具有很好的规律性，可以用于煤自燃预测预警的整个周期。从图3(b)、图3(d)可看出，CO2和C2H4体积分数在煤温为100 ℃以前没有明显变化，而在煤温达到100 ℃以后，能够反映煤温的剧烈变化，因此，可以作为煤自燃发生后期的预警气体。
2.2 SRU预测模型参数寻优
在标准PSO算法中，惯性权重
$\bar \omega $ 代表粒子保持前一时刻运动状态的能力，通常为固定值，当$\bar \omega $ 较大时，粒子在全局范围内进行搜索过程中受前一次迭代速度的影响，粒子的飞行距离较大，即使到了最优粒子所在的区域，也会不停围绕最优粒子振荡，无法快速靠近进行精细化搜索；当$\bar \omega $ 较小时，粒子得到比前一次迭代更小的速度，只能在很小的范围内飞行，此时虽然拥有较好的局部搜索能力，但容易陷入局部极小值[22-23]。2.2.1 改进PSO算法
为了平衡全局搜索和局部搜索，提高PSO算法寻优的性能，根据[−4，4]之间双曲正切曲线变化规律对标准PSO算法权重进行改进，使惯性权重随粒子搜索迭代次数动态改变。改进后的PSO算法权重更新公式为
$$ \bar \omega =\frac{{ {{\bar \omega _{\max }} + {\bar \omega _{\min }}} }}{2} + \tanh \left( - 4 + \frac{{A ({q_{\max }} - q)}}{{{q_{\max }} \left( {1+{{\rm{exp}}{ (- q)}}} \right)}}\right) \frac{{{\bar \omega _{\max }} - {\bar \omega _{\min }}}}{2} $$ (7) 式中：
$\bar \omega $ max为惯性权重最大值，取值为0.9；$\bar \omega $ min为惯性权重最小值，取值为0.4；A为控制全局搜索的参数，经过测试，取值为14； qmax为最大迭代次数。改进后的权重随迭代次数变化趋势如图4所示。
为了进一步避免算法陷入局部最小值，提高种群多样性和跳出局部极值的能力，从拟态物理学角度出发，将粒子调整自身朝个体最优位置和全局最优位置飞行的过程看作是受到二者的引力作用。与此同时，所有粒子都避免自己成为全局最差粒子，纷纷远离当前全局最差粒子，将粒子远离全局最差粒子构造为一股排斥力，使粒子陷入局部极值时迅速跳出[24-25]。在前人研究的基础上，综合考虑粒子在飞行过程中受到的引力和排斥力，即粒子更新将同时受到2个引力和1个排斥力的影响，对粒子速度更新进行改进，改进后的粒子速度更新公式为
$$ \begin{split} v_i^{q + 1} =& {\bar \omega _{}}v_i^q + {\alpha _1}\beta _1^q\left( {p_{{\rm{best}}i}^q - \lambda _i^q} \right) + {\alpha _2}\beta _2^q\left( {g_{{\rm{best}}i}^q - \lambda _i^q} \right) \\ &- \beta _3^q\left( {g_{{\rm{worst}}i}^q - \lambda _i^q} \right) \end{split} $$ (8) 式中：gworsti为全局最差个体粒子位置；β3为[0，1]内的随机数。
第q+1次迭代时，粒子λi受到的个体最优粒子对其的引力为
$p_{{\mathrm{best}}i}-\lambda_{i} $ ，受到全局最优粒子对其的引力为$g_{{\text {best}}i}-\lambda_{i} $ ，受到全局最差粒子对其的排斥力为$g_{{\mathrm{worst}}i}-\lambda_{i} $ ，加入排斥力后，能使粒子更快地向全局范围内适应度小的位置聚集，提高算法的收敛速度。2.2.2 改进PSO优化SRU预测模型
（1） 将获取的煤自燃温度和气体指标原始数据进行数据预处理，将数据划分为训练集和测试集。
（2） 确定参数寻优范围，初始化粒子群参数。最大迭代次数设为1 000；寻优范围设置如下：2个隐藏层神经元数的取值范围均为[1，200]，学习率边界范围为[0.0001，0.01]，训练次数寻优范围为[1，500]。
（3） 利用当前粒子群中各粒子的参数取值建立SRU煤自燃温度预测模型，学习训练集中气体指标与煤自燃温度之间的非线性规律，再对测试集样本进行预测。
（4） 计算各粒子适应度，更新共享信息。记录当前取得最小适应度的粒子，更新粒子群的全局最优位置
$g_{{\text {best}}i} $ ；更新每个粒子从搜索开始后各自的历史最优位置$p_{{\text {best}}i} $ ；记录当前取得最大适应度的粒子，更新粒子群的全局最差位置$g_{{\text {worst}}i} $ 。适应度计算公式：$$ \psi = \frac{1}{n}\sum\limits_{t = 1}^n {\left| {{{P}'_t} - {P_t}} \right|} $$ (11) 式中：
$\psi $ 为计算出的训练样本真实值与预测值之间的误差适应度；$ {P'_t} $ 为第t个时刻样本的煤自燃温度预测值；Pt为第t个时刻样本的煤自燃温度真实值。（5） 判断是否达到终止条件。若未达到终止条件，根据迭代次数更新
$\bar \omega $ ，并依据式(8)、式(6)更新每个粒子的速度和位置，然后转步骤（2）继续执行；否则执行步骤（6）。（6） 输出全局最小适应度及全局最小适应度对应粒子的位置参数，该位置参数即为全局最优解。
3. 实验与结果分析
实验环境配置如下：CPU型号为i7−3612QM；GPU型号为RTX1080Ti；内存容量为12 GB；编程语言为Python−3.6；编程平台为PyCharm−2019.2.5(Professional Edition)；集成环境管理为Anaconda Navigator−1.3.1；深度学习框架为PyTorch−1.6；机器学习框架为Scikit−learn−0.23.2。
3.1 预警指标择优
煤氧复合作用下，煤体氧化升温过程中会释放不同浓度的气体产物，需要分析各气体指标与温度在数据上的紧密程度，筛选出与煤自燃温度相关性尽可能大的气体指标。对插补过的指标气体数据与煤样温度数据进行相关性分析，对数据的每一列均采用Kolmogorov−Smirnov检验，其P−value均小于0.001，采用Pearson相关系数法分析相关性。
$$\theta = \frac{{\displaystyle\sum\limits_{t = 1}^n {({{{x}}_t} - \bar x)({{{y}}_t} - \bar y)} }}{{\sqrt {\displaystyle\sum\limits_{t = 1}^n {{{({{{x}}_t} - \bar x)}^2}{{({{{y}}_t} - \bar y)}^2}} } }} $$ (12) 式中：θ 为线性相关的程度；xt, yt为任意2个特征指标的样本值；
$ \bar x $ ，$ \bar y $ 为对应特征指标的均值。选取与温度相关系数大于0.7的气体指标O2，CO，CO2，CH4，C2H4作为预测模型输入数据，煤温作为预测模型目标输出，将这些数据样本按照70%训练集、30%测试集进行划分。各指标间相关性系数见表1。
表 1 温度与气体指标间的相关性Table 1. Correlation between temperature and gas indexes名称 温度 O2 CO CO2 CH4 C2H4 温度 1.00 −0.74 0.80 0.82 0.81 0.72 O2 −0.74 1.00 −0.69 −0.85 −0.72 −0.65 CO 0.80 −0.69 1.00 0.83 0.68 0.61 CO2 0.82 −0.85 0.83 1.00 0.67 0.60 CH4 0.81 −0.72 0.68 0.67 1.00 0.92 C2H4 0.72 −0.65 0.61 0.60 0.92 1.00 3.2 评价指标计算
为了评估PSO−SRU模型性能，对模型的预测结果作出统一客观的评价，选取平均绝对误差(Mean Absolute Error, MAE)、均方根误差(Root Mean Square Error, RMSE)、决定系数R2作为评价指标。
(1) MAE计算公式为
$$ {{E_{\rm{MA}}}} = \frac{1}{n}\sum\limits_{t = 1}^n {\left| {{{P}'_t} - {P_t}} \right|} $$ (13) MAE是预测值与真实值之间绝对误差的平均值，能避免误差相互抵消。MAE越小，说明预测越精准。
(2) RMSE计算公式为
$$ {{E_{\rm{RMS}} }}= \sqrt {\frac{1}{n}\sum\limits_{t = 1}^n {{{\left( {{{P}'_t} - {P_t}} \right)}^2}} } $$ (14) RMSE能体现出预测值的离散程度，RMSE越大，说明有较多或较大的偏差值出现，RMSE越小，说明模型性能越好。
(3) R2计算公式为
$$ {{{R}}^2}{\text{ = }}{{\sum\limits_{t = 1}^n {{{\left( {{{P}'_t} - \bar P} \right)}^2}} } \mathord{\left/ {\vphantom {{\sum\limits_{t = 1}^n {{{\left( {{{P}'_t} - \bar P} \right)}^2}} } {\sum\limits_{t = 1}^n {{{\left( {{P_t} - \bar P} \right)}^2}} }}} \right. } {\sum\limits_{t = 1}^n {{{\left( {{P_t} - \bar P} \right)}^2}} }} $$ (15) 式中
$ \bar P $ 为真实值的平均值。R2可以判别模型拟合非线性关系的程度，取值范围为[0,1]，数值越接近1，表示模型预测性能越好。
3.3 实验结果分析
3.3.1 模型结构确定
为了选择合适的隐藏层数，本文利用改进的PSO算法分别优化隐藏层的SRU预测模型，改变不同隐藏层数模型的训练次数，进行MAE和运行时间对比，结果如图5所示。
![]()
图 5 不同隐藏层数SRU预测模型MAE与时间对比Figure 5. Comparison of mean absolute errors(MAE) and running time under various hidden layer of simple recurrent units(SRU) perdiction model从图5可看出，不同隐藏层SRU预测模型的MAE随训练次数增加逐渐减小，运行时间随隐藏层数增加几乎为线性增长。1层隐藏层模型的最小MAE始终大于2，3，4层模型；3层隐藏层模型的最小MAE略高于2，4层模型；4层隐藏层模型的MAE最小，虽然比2层模型的MAE略小，但运行时间是2层模型的2.1倍。因此，在综合考虑预测精度和时效性的前提下，本文采用含有2层隐藏层的SRU构建PSO−SRU模型。
3.3.2 模型参数优化
SRU预测模型由输入层、2层隐藏层和1层输出层构成。为了验证改进的PSO算法性能，分别采用改进的PSO算法和PSO算法对含有2层SRU隐藏层的模型进行优化对比。寻优过程中二者的适应度变化曲线如图6所示。从图6可看出，改进的PSO算法明显比PSO算法可更快求得全局最优解，并最终求解到全局最小适应度，说明改进的PSO算法全局寻优能力更强，算法收敛速度更快。基于改进的PSO算法最终寻找到的最优SRU预测模型参数如下：L1=18，L2=18，ɛ=0.001，e=21。
![]()
图 6 PSO算法改进前后适应度变化曲线Figure 6. Fitness value change curves before and after improving particle swarm optimization(PSO) algorithm3.3.3 性能对比
为了验证PSO−SRU模型的性能，在使用相同训练集和测试集数据前提下，利用改进的PSO算法对各对比模型参数进行优化。基于SVR算法的预测模型[16]惩罚因子为9.16，核参数为12.47；基于BP算法的预测模型[18]2层隐藏层神经元数分别为12和8；基于RF算法的预测模型[19]决策树数量为115，最大树深为5。采用MAE、RMSE和R2指标对不同预测模型的性能进行评价，结果见表2。
表 2 不同预测模型的性能预测结果对比Table 2. Comparison of predictions using various models模型 MAE/℃ RMSE/℃ R2 训练集 测试集 训练集 测试集 训练集 测试集 SVR 17.35 18.01 27.21 31.34 0.86 0.83 RF 10.81 13.08 14.86 26.14 0.91 0.86 BP 7.49 11.34 12.51 17.63 0.97 0.89 POS−SRU 3.15 5.43 6.27 8.69 0.99 0.96 由表2分析可知，PSO−SRU模型在测试集上的决定系数R2比基于SVR，RF和BP算法的预测模型分别提高了15.67%，11.63%和7.87%，且PSO−SRU模型的R2为0.96，说明其预测效果及性能最好，预测结果最接近真实值。在训练集和测试集上的预测结果表明：基于BP算法的预测模型容易出现过拟合；基于SVR和RF算法的预测模型在训练样本和测试样本上预测误差差距较小、泛化性较好，但预测精度不足；PSO−SRU模型在训练集和测试集上的R2差值为0.03，说明PSO−SRU模型具有良好的泛化性。相较于基于SVR，RF和BP算法的预测模型，PSO−SRU模型在测试集上的MAE分别降低了12.58，7.65，5.91 ℃，说明PSO−SRU模型在测试样本上的预测精度优于其他3个模型；同时，RMSE分别降低了22.65，17.45，8.94 ℃，说明PSO−SRU模型预测的煤自燃温度离散化程度低，具有良好的鲁棒性。
不同模型测试样本真实煤温与预测煤温对比如图7所示。
![]()
图 7 不同模型测试样本真实煤温与预测煤温对比Figure 7. Comparison of real and predicted temperatures of testing samples using different models从图7可看出，基于SVR算法的预测模型预测值偏离真实值程度最大，基于RF算法的预测模型预测值离散化程度较大，模型鲁棒性较差，基于BP算法的预测模型预测值曲线虽然比基于SVR，RF算法的预测模型更贴近真实值，但预测精度仍然低于PSO−SRU模型。PSO−SRU模型预测的煤温曲线平滑，验证了其用于提高煤自燃温度预测精度、鲁棒性和泛化能力的有效性。
4. 结论
（1） PSO−SRU模型利用SRU挖掘煤自燃温度与气体指标特征之间的非线性关系，根据历史监测数据预测出煤自燃温度，并利用改进的PSO算法对模型参数进行优化。相较于标准PSO算法，改进的PSO算法在优化SRU煤自燃温度预测模型超参数的过程中，降低了预测模型选取参数的难度和模型结构复杂程度，提高了参数寻优效率。
（2） 实验结果表明：与基于SVR，RF，BP算法的煤自燃温度预测模型相比，PSO−SRU模型在测试集上决定系数R2分别提高了15.67%，11.63%，7.87%，且R2为0.96，预测结果最接近真实值；MAE分别降低了12.58，7.65，5.91 ℃，预测精度得到了提升； RMSE分别降低了22.65，17.45，8.94 ℃，预测曲线在拐点处更平滑，说明PSO−SRU模型具有良好的鲁棒性。
-
![]()
图 1 SRU内部结构和网络结构
Figure 1. Simple recurrent units(SRU) interior structure and network structure
![]()
图 2 PSO−SRU模型架构
Figure 2. Temperature prediction model framework for coal spontaneous combustion based on particle swarm optimization and simple recurrent unit (PSO-SRU)
![]()
图 5 不同隐藏层数SRU预测模型MAE与时间对比
Figure 5. Comparison of mean absolute errors(MAE) and running time under various hidden layer of simple recurrent units(SRU) perdiction model
![]()
图 6 PSO算法改进前后适应度变化曲线
Figure 6. Fitness value change curves before and after improving particle swarm optimization(PSO) algorithm
![]()
图 7 不同模型测试样本真实煤温与预测煤温对比
Figure 7. Comparison of real and predicted temperatures of testing samples using different models
表 1 温度与气体指标间的相关性
Table 1 Correlation between temperature and gas indexes
名称 温度 O2 CO CO2 CH4 C2H4 温度 1.00 −0.74 0.80 0.82 0.81 0.72 O2 −0.74 1.00 −0.69 −0.85 −0.72 −0.65 CO 0.80 −0.69 1.00 0.83 0.68 0.61 CO2 0.82 −0.85 0.83 1.00 0.67 0.60 CH4 0.81 −0.72 0.68 0.67 1.00 0.92 C2H4 0.72 −0.65 0.61 0.60 0.92 1.00 
下载: 导出CSV
表 2 不同预测模型的性能预测结果对比
Table 2 Comparison of predictions using various models
模型 MAE/℃ RMSE/℃ R2 训练集 测试集 训练集 测试集 训练集 测试集 SVR 17.35 18.01 27.21 31.34 0.86 0.83 RF 10.81 13.08 14.86 26.14 0.91 0.86 BP 7.49 11.34 12.51 17.63 0.97 0.89 POS−SRU 3.15 5.43 6.27 8.69 0.99 0.96
下载: 导出CSV
-
[1] 谢和平,吴立新,郑德志. 2025年中国能源消费及煤炭需求预测[J]. 煤炭学报,2019,44(7):1949-1960. XIE Heping,WU Lixin,ZHENG Dezhi. Prediction on the energy consumption and coal demand of China in 2025[J]. Journal of China Coal Society,2019,44(7):1949-1960.
[2] 林柏泉,李庆钊,周延. 煤矿采空区瓦斯与煤自燃复合热动力灾害多场演化研究进展[J]. 煤炭学报,2021,46(6):1715-1726. LIN Baiquan,LI Qingzhao,ZHOU Yan. Research advances about multi-field evolution of coupled thermodynamic disaster in coal mine goaf[J]. Journal of China Coal Society,2021,46(6):1715-1726.
[3] 邓军,白祖锦,肖旸,等. 煤自燃灾害防治技术现状与挑战[J]. 煤矿安全,2020,51(10):118-125. DENG Jun,BAI Zujin,XIAO Yang,et al. Present situation and challenge of coal spontaneous combustion disasters prevention and control technology[J]. Safety in Coal Mines,2020,51(10):118-125.
[4] 程卫民,张孝强,王刚,等. 综放采空区瓦斯与遗煤自燃耦合灾害危险区域重建技术[J]. 煤炭学报,2016,41(3):662-671. CHENG Weimin,ZHANG Xiaoqiang,WANG Gang,et al. Reconstruction technology of gas and coal spontaneous combustion coupled hazard in fully mechanized caving goaf[J]. Journal of China Coal Society,2016,41(3):662-671.
[5] 张玉涛,李亚清,邓军,等. 煤炭自燃灾变过程突变特性研究[J]. 中国安全科学学报,2015,25(1):78-84. ZHANG Yutao,LI Yaqing,DENG Jun,et al. Study on catastrophe characteristics of coal spontaneous combustion[J]. China Safety Science Journal,2015,25(1):78-84.
[6] 秦波涛,仲晓星,王德明,等. 煤自燃过程特性及防治技术研究进展[J]. 煤炭科学技术,2021,49(1):66-99. QIN Botao,ZHONG Xiaoxing,WANG Deming,et al. Research progress of coal spontaneous combustion process characteristics and prevention technology[J]. Journal of Coal Science and Technology,2021,49(1):66-99.
[7] 邓军,徐精彩,陈晓坤. 煤自燃机理及预测理论研究进展[J]. 辽宁工程技术大学学报,2003,22(4):455-459. DENG Jun,XU Jingcai,CHEN Xiaokun. Perspectives on spontaneous combustion mechanism and prediction theory of coal[J]. Journal of Liaoning Technical University,2003,22(4):455-459.
[8] 朱红青,王海燕,王斐然,等. 煤堆测温技术研究进展[J]. 煤炭科学技术,2014,42(1):50-54. ZHU Hongqing,WANG Haiyan,WANG Feiran,et al. Research progress on coal stockpile temperature measuring technology[J]. Coal Science and Technology,2014,42(1):50-54.
[9] 陈欢,杨永亮. 煤自燃预测技术研究现状[J]. 煤矿安全,2013,44(9):194-197. CHEN Huan,YANG Yongliang. Research status of predicting coal spontaneous combustion[J]. Safety in Coal Mines,2013,44(9):194-197.
[10] 李林,陈军朝,姜德义,等. 煤自燃全过程高温区域及指标气体时空变化实验研究[J]. 煤炭学报,2016,41(2):444-450. LI Lin,CHEN Junchao,JIANG Deyi,et al. Experimental study on temporal variation of high temperature region and index gas of coal spontaneous combustion[J]. Journal of China Coal Society,2016,41(2):444-450.
[11] 郝宇,叶正亮. 不同甲烷气氛下煤自燃指标气体及活化能研究[J]. 工矿自动化,2019,45(11):65-69. HAO Yu,YE Zhengliang. Research on index gas and activation energy of coal spontaneous combustion under different methane atmosphere[J]. Industry and Mine Automation,2019,45(11):65-69.
[12] 周冬,刘贞堂,钱继发,等. 采空区内煤自燃气体特征及产生规律分析[J]. 工矿自动化,2019,45(3):18-22. ZHOU Dong,LIU Zhengtang,QIAN Jifa,et al. Analysis of gas characteristics and generation rules of coal spontaneous combustion in goaf[J]. Industry and Mine Automation,2019,45(3):18-22.
[13] 金永飞,郭军,文虎,等. 煤自燃高温贫氧氧化燃烧特性参数的实验研究[J]. 煤炭学报,2015,40(3):596-602. JIN Yongfei,GUO Jun,WEN Hu,et al. Experimental study on the high temperature lean oxygen oxidation combustion characteristic parameters of coal spontaneous combustion[J]. Journal of China Coal Society,2015,40(3):596-602.
[14] 高峰,王文才,李建伟,等. 浅埋煤层群开采复合采空区煤自燃预测[J]. 煤炭学报,2020,45(增刊1):336-345. GAO Feng,WANG Wencai,LI Jianwei,et al. Prediction of coal spontaneous combustion in compound gob of shallow seam group mining[J]. Journal of China Coal Society,2020,45(S1):336-345.
[15] LEI Changkui,DENG Jun,CAO Kai,et al. A random forest approach for predicting coal spontaneous combustion[J]. Fuel,2018,223(1):63-73.
[16] 邓军,雷昌奎,曹凯,等. 煤自燃预测的支持向量回归方法[J]. 西安科技大学学报,2018,38(2):175-180. DENG Jun,LEI Changkui,CAO Kai,et al. Support vector regression approach for predicting coal spontaneous combustion[J]. Journal of Xi'an University of Science and Technology,2018,38(2):175-180.
[17] 刘宝,穆坤,叶飞,等. 基于相关向量机的煤自燃预测方法[J]. 工矿自动化,2020,46(9):104-108. LIU Bao,MU Kun,YE Fei,et al. Prediction method of coal spontaneous combustion based on relevance vector machine[J]. Industry and Mine Automation,2020,46(9):104-108.
[18] 昝军才,魏成才,蒋可娟,等. 基于BP神经网络的煤自燃温度预测研究[J]. 煤炭工程,2019,51(10):113-117. ZAN Juncai,WEI Chengcai,JIANG Kejuan,et al. Prediction of coal spontaneous combustion temperature based on BP neural network[J]. Coal Engineering,2019,51(10):113-117.
[19] 郑学召,李梦涵,张嬿妮,等. 基于随机森林算法的煤自燃温度预测模型研究[J]. 工矿自动化,2021,47(5):58-64. ZHENG Xuezhao,LI Menghan,ZHANG Yanni,et al. Research on the prediction model of coal spontaneous combustion temperature based on random forest algorithm[J]. Industry and Mine Automation,2021,47(5):58-64.
[20] TAO Lei, YU Zhang, WANG Sidai, et al. Simple recurrent units for highly parallelizable recurrence[C]// Proceedings of the 2018 Conference on Empirical Methods in Natural Language Processing, Brussels, 2018: 4470-4481.
[21] KENNEDY J, EBERHART R C. Particle swarm optimization[C]//Proceedings of the IEEE International Conference on Neural Networks, Perth, 1995: 1942-1948.
[22] 闫群民,马瑞卿,马永翔,等. 一种自适应模拟退火粒子群优化算法[J]. 西安电子科技大学学报,2021,48(4):120-127. YAN Qunmin,MA Ruiqing,MA Yongxiang,et al. Adaptive simulated annealing particle swarm optimization algorithm[J]. Journal of Xidian University,2021,48(4):120-127.
[23] 敖永才,师奕兵,张伟,等. 自适应惯性权重的改进粒子群算法[J]. 电子科技大学学报,2014,43(6):874-880. AO Yongcai,SHI Yibing,ZHANG Wei,et al. Improve particle swarm optimization with adaptive inertia weight[J]. Journal of University of Electronic Science and Technology of China,2014,43(6):874-880.
[24] 汪雅文, 钱谦, 冯勇, 等. 融合吸引排斥和双向学习的改进粒子群算法[J/OL]. 计算机工程与应用: 1-9. [2021-07-20]. http://kns.cnki.net/kcms/detail/11.2127.TP.20210618.0839.002.html. WANG Yawen, QIAN Qian, FENG Yong, et al. An improved particle swarm optimization algorithm based on attraction-repulsion and bidirectional learning strategies[J/OL]. Computer Science and Application: 1-9. [2021-07-20]. http://kns.cnki.net/kcms/detail/11.2127.TP.20210618.0839.002.html.
[25] 姚成玉,赵哲谕,陈东宁,等. 有向动态拓扑混合作用力微粒群优化算法及可靠性应用[J]. 机械工程学报,2017,53(10):166-179. DOI: 10.3901/JME.2017.10.166 YAO Chengyu,ZHAO Zheyu,CHEN Dongning,et al. Unidirectional dynamic topology hybrid force PSO algorithm and its applications in reliability optimization[J]. Journal of Mechanical Engineering,2017,53(10):166-179. DOI: 10.3901/JME.2017.10.166
-
期刊类型引用(8)
1. 孙吉平. 基于SF_6质量浓度变化特征的煤矿火灾状态识别分析. 山西煤炭. 2025(01): 42-49 . 
百度学术
2. 李延河，万志军，于振子，苟红，赵万里，周嘉乐，师鹏，甄正，张源. 基于PSO-SVR的掘进工作面风温预测. 煤炭科学技术. 2025(01): 183-191 . 百度学术
3. 王斌，贾澎涛，郭风景，孙刘咏，林开义. 基于多特征融合的煤自燃温度深度预测模型. 中国矿业. 2024(02): 84-90 . 百度学术
4. 邓军，李鑫，王凯，王伟峰，闫军，汤宗情，康付如，任帅京. 矿井火灾智能监测预警技术近20年研究进展及展望. 煤炭科学技术. 2024(01): 154-177 . 百度学术
5. 曹富荣，吴学松，李军，付天予，刘佳伟，李志辉，杨小彬. 基于机器学习的多气体指标煤自燃温度预测. 煤矿安全. 2024(04): 106-113 . 百度学术
6. 薛凯隆，崔欣超，祁云，齐庆杰. 基于DBO-SVM的采空区煤自燃危险性预测. 沈阳理工大学学报. 2024(06): 85-90 . 百度学术
7. 贾澎涛，张杰，郭风景. 基于GAT-Informer的采空区煤自燃温度预测模型. 工矿自动化. 2024(11): 92-98+108 . 本站查看
8. 赵春光，周李兵，高鹏，刘军伟. 基于改进PSO算法的煤矿相机参数优化设计. 煤矿机电. 2023(05): 1-6 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 222
- HTML全文浏览量: 104
- PDF下载量: 39
- 被引次数: 12
