1.一种基于增量式极限学习机的煤矿安全风险预警方法，其特征在于，包括：构建基于增量式极限学习机的煤矿安全风险预测模型；实时获取煤矿数据，对该数据进行预处理，将预处理后的数据输入到训练好的基于增量式极限学习机的煤矿安全风险预测模型中，得到该煤矿安全风险预测结果；根据煤矿安全风险预测结果对该煤矿进行预警；

对基于增量式极限学习机的煤矿安全风险预测模型进行训练的过程包括：S1：获取历史煤矿数据，对获取的数据进行预处理，得到规整的隐患时间序列数据；煤矿数据包括煤矿每日隐患数量、隐患类型、隐患级别、矿井历史事故记录以及矿井基础数据；

S2：设置时间粒度t，根据设置的时间粒度对历史煤矿数据进行划分，得到训练样本；

S3：基于极限学习机构建非线性的煤矿生产安全风险模型，设置模型训练时的MSE阈值；

S4：初始化模型的超参数，并采用遗传算法筛选出最优的模型超参数；

S41：初始化种群数、种群大小以及迭代次数；并对权重W与偏置bias进行随机取值；

S42：计算种群的适应度，计算适应度函数的表达为：

F＝min||MSETrainingset,MSETestingset||其中MSETrainingset为训练集的均方误差，MSETestingset为测试集的均方误差；

S43：对种群进行选择、交叉以及变异操作，产生新的子种群；

S44：计算新的子种群的适应度，若适应度函数值小于上一次适应度值，则更新模型权重W和隐藏层偏置bias；

S45：当达到最大迭代次数时，得到最优的W和bias；

S5：将数据集放入煤矿生产安全风险模型，根据筛选出的超参数对模型进行初训练，并通过广义逆矩阵获取模型输出层权重，并根据初训练结果选取最低MSE下的时间粒度t；

2

S51：定义训练损失函数Function(H,T)＝min‖Hβ‑T‖,其中H为输出矩阵，β为输出层权重，T为实际值，g表示激活函数， 表示第 个输入权重，xN表示第N个变量， 表示第 个偏置， 表示第 个输出权重， 表示第N个标签值，N表示训练集的个数， 表示隐藏层的数量；

S52：模型训练的过程即求解使得损失函数达到最小值时的输出层权重β，通过矩阵形式的最小二乘法计算模型的输出权重，其表达式为：T ‑1 T

β＝(HH) HT

S53：通过T＝Hβ对测试集进行预测，并根据均方误差MSE选取最优时间粒度t；

S6：根据最低MSE下的时间粒度t对模型进行增量训练，当预测MSE高于步骤S3中设置的MSE阈值时，对模型参数进行增量更新；

S7：对模型隐藏层数量进行更新，完成模型的训练过程。

2.根据权利要求1所述的一种基于增量式极限学习机的煤矿安全风险预警方法，其特征在于，对数据进行预处理的过程包括：对于一个矿区每天的隐患检查记录，按照隐患类型，隐患级别，检查单位提取同一时间段非缺失隐患数据，计算得到均值与方差，并利用计算得到的均值与方差对所有缺失值进行填充，得到规整的多粒度隐患度监测时间序列。

3.根据权利要求1所述的一种基于增量式极限学习机的煤矿安全风险预警方法，其特征在于，煤矿安全风险预测模型超参数包括：模型隐藏层数L，最优时间粒度t、模型权重W、隐藏层偏置bias以及模型输出层权重βi。

4.根据权利要求1所述的一种基于增量式极限学习机的煤矿安全风险预警方法，其特征在于，对模型进行增量训练的过程包括：S61：当预测MSE高于设定的MSE阈值时，定义增加训练数据后的损失函数H0为上一次的隐藏层输出矩阵，H1为新增数据的输出矩阵；T0表示上一次的标签值，T1表示新增的标签值，β1表示需要求解的新的输出权重；

S6 2：根据 广 义逆 矩阵 的计 算方 法可 得 其中 ：K1表示输出矩阵与其转置矩阵的乘积；

S63：将增加样本后的公式转换成由已知常量β0,K1,H1和T1表示的函数，得到增量式极限学习机训练过程的递推公式；其表达式为：其中，Kk+1表示第k+1次增量训练时隐层输出矩阵H转置相乘的结果；βk+1表示第k+1次增量训练时新的输出层权重的结果；Tk+1表示第k+1次增量训练时的训练标签；

S64：基于更新的输出层权重的结果β对测试集进行预测，预测结果为T＝Hβ。

5.根据权利要求1所述的一种基于增量式极限学习机的煤矿安全风险预警方法，其特征在于，对模型隐藏层数量进行更新的过程包括：S71：更新模型的权重W和偏置的值，权重更新公式为：Wl+1＝Wlw，偏置更新公式为bl+1＝blb，其中Wl和bl是更新之前含有l个隐含层节点数的权重和偏置，w为新增权重，b为新增偏置，且w和b为任意值；

S72：根据更新后的模型权重和偏置计算输出层矩阵；计算的表达式为：l+1 l l+1 l l+1 l

H ＝g([WW ]X+[bb ])＝[Hh]

l l+1 l

其中，g表示激活函数，X表示特征值，W表示更新之前的权重，W 表示新增的权重，b表l+1 l示更新之前的偏置，b 表示新增偏置，H表示隐层输出矩阵，h表示新增的输出矩阵；

S73：根据输出层矩阵更新隐藏层后的输出层权重，其表达式为：l+1 l+1 + l T l ‑1 l T

β ＝(H ) Y＝([Hh][Hh]) [Hh]Yl+1 + l+1 l

其中，(H ) 表示输出矩阵H 的可逆矩阵，Y表示标签值，H表示原有的输出矩阵，h表示新增的输出矩阵。

6.根据权利要求1所述的一种基于增量式极限学习机的煤矿安全风险预警方法，其特征在于，基于增量式极限学习机的煤矿安全风险预测模型处理数据的过程包括：其中，L表示隐藏层个数，g为激活函数，Wi表示模型中的权重，bias为隐藏层偏置，t表示时间窗口，xt表示t时间窗口下相关特征值，n表示为煤矿中矿井数量，βi表示输出层权重。

7.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行，以实现权利要求1至6中任一项基于增量式极限学习机的煤矿安全风险预警方法。

8.一种基于增量式极限学习机的煤矿安全风险预警装置，其特征在于，包括处理器和存储器；所述存储器用于存储计算机程序；所述处理器与所述存储器相连，用于执行所述存储器存储的计算机程序，以使所述一种基于增量式极限学习机的煤矿安全风险预警装置执行权利要求1至6中任一项基于增量式极限学习机的煤矿安全风险预警方法。