一种煤矿井下涌水量预测方法

技术领域

本发明涉及煤炭开采水害防治技术领域，更具体的涉及一种基于水泵性能的煤矿井下每日涌水量预测方法。

背景技术

随着东部区煤炭资源的逐渐枯竭，煤炭资源的开发利用逐步向西部转移，西部深部煤矿工作面高强度回采过程中，面临工作面顶板大量涌水的威胁。为了消除井下大量涌水对生产工作的干扰及可能产生的水害威胁，需要提前预测工作面回采过程中的涌水量，然后依据预测的涌水量设计、建设工作面排水系统。

现有针对煤矿井下涌水量预测的常规方法，需要大量人工测量采集数据，不仅消耗大量的财力、人力及物力，还会受到短期天气等原因的影响。例如针对于煤矿井下涌水量预测，使用的方法普遍为依据每月/每季/每年人工测量煤矿井下含水层厚、引用半径、影响半径、渗透系数、采掘深度等参数统计计算出的涌水量，使用灰色模型预测、BP网络模型等预测模型进行月/季/年维度下的涌水量预测，无法依据预测出来的涌水量与水泵性能进行联系，导致涌水量预测的精度差，进而不能准确对每日水泵的开泵计划进行评估。

发明内容

针对上述领域中存在的问题，本发明一种煤矿井下涌水量预测方法，能够解无法依据预测出来的涌水量与水泵性能进行联系，导致涌水量预测的精度差，进而不能准确对每日水泵的开泵计划进行评估的技术问题。

为解决上述技术问题，本发明公开了一种煤矿井下涌水量预测方法，包括以下步骤：

根据待预测的井下各个泵房的每日排水量、水仓液位、以及涌水速率，计算每日的井下各个泵房的涌水量、日水仓液位变化、涌水速率变化及各个泵房的涌水量与总涌水量的占比；

依据每日的井下各个泵房的涌水量，使用一次指数平滑算法进行第一次明日涌水量预测，获得第一次明日涌水量预测值；根据每日的涌水速率变化、日水仓液位变化及各个泵房的涌水量与总涌水量的占比，对获得的第一次明日涌水量预测值进行矫正，作为第二次明日涌水量预测值；

根据井下各个泵房的每日排水量、日水仓液位变化、涌水速率、前一日涌水量、第一次明日涌水量预测值、第二次明日涌水量预测值，对选取的泵房的涌水量进行单独分析，确定自变量因子和因变量因子，使用梯度提升树回归算法进行涌水量预测，对应获得选取的各个泵房涌水量预测值；

根据井下各个泵房的每日排水量、涌水速率、日水仓液位变化、前一日涌水量、第一次明日涌水量预测值、第二次明日涌水量预测值以及获得的选取的各个泵房涌水量预测值，使用梯度提升树回归算法，对井下各个泵房的涌水量进行二次涌水量预测，获得井下各个泵房的最终涌水量预测值；

根据获得的井下各个泵房的最终涌水量预测值与各个泵房的各水泵性能参数，分析计算获得每个泵房的水泵开泵时长预测值及开泵耗能预测值。

优选地，所述计算每日的井下各个泵房的涌水量、日水仓液位变化、涌水速率变化及各个泵房的涌水量与总涌水量的占比包括通过自动化系统采集井下各个泵房的每日排水量、水仓液位、以及涌水速率数据并存入sqlserver数据库中，利用python对采集的数据按照时间序列的顺序进行排序，根据采集来的数据计算出每日的井下各个泵房的涌水量及总涌水量。

优选地，所述获得第一次明日涌水量预测值的计算公式为：

F

其中，F

优选地，所述平滑常数α的取值范围为0.01～0.25。

优选地，所述对应获得选取的各个泵房涌水量预测值包括以下步骤：

分别对井下中央泵房、井下一号采区清水泵房、井下一号采区污水泵房、井下二号采区污水泵房的涌水量预测进行单独分析；

以井下中央泵房、井下一号采区清水泵房、井下一号采区污水泵房、井下二号采区污水泵房每个泵房的每日排水量、涌水速率、日水仓液位变化、前一日涌水量、第一次明日涌水量预测值、第二次明日涌水量预测值为自变量因子；

各个泵房的涌水量预测为因变量因子；

使用梯度提升树回归算法分别预测每个泵房的涌水量预测值，分别作为井下中央泵房涌水量预测值、井下一号采区清水泵房涌水量预测值、井下一号采区污水泵房涌水量预测值、井下二号采区泵房污水泵涌水量预测值。

优选地，所述获得井下各个泵房的最终涌水量预测值包括以下步骤：

以井下各个泵房的每日排水量、涌水速率、日水仓液位变化、前一日涌水量、第一次明日涌水量预测值、第二次明日涌水量预测值以及中央泵房涌水量预测值、井下一号采区清水泵涌水量预测值、井下一号采区污水泵涌水量预测值、井下二号采区泵房污水泵涌水量预测值作为自变量因子；

涌水量预测值为因变量因子；

使用梯度提升树回归算法进行二次涌水量预测，获得井下各个泵房的最终涌水量预测值。

优选地，所述梯度提升树回归算法包括以下步骤：

初始化

迭代建立m棵决策树，以训练数据集

对r

利用线性搜索估计叶节点区域的值

G

以各个泵房的每日的涌水量数据作为训练集训练模型。

优选地，所述梯度提升树回归算法中的决策树个数m为20。

优选地，所述使用梯度提升树回归算法进行二次涌水量预测的训练集为井下各个泵房的每日排水量、涌水速率、日水仓液位变化、第一次明日涌水量预测值、第二次明日涌水量预测值、前一日涌水量以及获得的选取的各个泵房涌水量预测值。

与现有技术相比，本发明具有如下有益效果：

本发明根据计算的每日的井下各个泵房的涌水量，使用一次指数平滑法预测涌水量，得到明日涌水量预测值，并对其进行矫正，可以对明日涌水量精准预测；根据明日涌水量预测值、矫正后的明日涌水量预测值，以及增加多个自变量参数，两次使用梯度提升树回归算法对井下各个泵房的涌水量进行预测，该方法将一次指数平滑法和两次梯度提升树回归算法相结合，提高了涌水量的预测精准度。该方法通过对明日涌水量的精准预测，并结合水泵性能参数，对各个泵房水泵的开泵时长与能耗进行预测，进而能够对每日水泵的开泵计划进行准确评估。

附图说明

图1为本发明的整体方法流程图；

图2为本发明的明日涌水量预测图；

图3为本发明的中央泵房涌水量预测图；

图4为本发明的最终涌水量预测图。

具体实施方式

下面将结合本发明实施例中的附图1-4，对本发明实施例中的技术方案进行清楚、完整地描述。应理解本发明中所述的术语仅仅是为描述特别的实施方式，并非用于限制本发明。

实施例

如图1所示，本发明实施例提供一种煤矿井下涌水量预测方法，包括以下步骤：

S1：更改涌水量预测的数据源，通过自动化系统采集井下各个泵房的每日排水量、水仓液位、以及涌水速率数据并存入sqlserver数据库中，如表1-2所示。

表1自动化系统采集存储于sqlserver中数据表中的数据

表2采集每个泵房每台泵的日运行数据

利用python对采集的数据按照时间序列的顺序进行排序，并依据各泵房的每日排水量、水仓液位以及涌水速率，计算出每日的井下各个泵房的涌水量及总涌水量，计算结果如表3所示，根据每日的井下各个泵房的涌水量及总涌水量，进一步计算日水仓液位变化、涌水速率变化及各个泵房的涌水量与总涌水量的占比。

表3各个泵房的水仓涌水速率及总涌水量数据

指数平滑法是一种性能优良、适应性强的方法，在各个方面都有广泛的应用。其作用主要体现在两个方面：一是用于预测，二是用于修匀历史数据，以测定时间数列的长期趋势。指数平滑法具有显著的特征，它有效地利用了全部历史数据且操作简单易行。

一次指数平滑法是一种加权预测，权数为α。它既不需要存储全部历史数据，也不需要存储一组数据，从而可以大大减小数据存储问题。甚至有时只需要一个最新观察值、最新预测值和α值，就可以进行预测。它提供的预测值是前一期预测值加上前期预测值中产生的误差的修正值。所以，将此算法应用于井下各个泵房的涌水量预测中，不仅可以解决数据非线性化的问题，同时，可以解决数据的延迟性问题。

S2：依据每日的井下各个泵房的涌水量，使用一次指数平滑算法进行第一次明日涌水量预测，计算公式为：

F

其中，F

获得第一次明日涌水量预测值，如表4所示。

表4第一次明日涌水量预测数据

如表4所示，为第一次明日涌水量预测值的部分数据，其中，真实值为当天的涌水实际数据，预测值为依据指数平滑预测出来的涌水值，预测上限与预测下限为预测出来的上限值与下限值，通过计算均方根误差、平均绝对误差和平均相对误差，得出均方根误差为129.7512，平均绝对误差为66.436，平均相对误差为0.0083。

通过分析可知，平均相对误差较小，平均绝对误差偏大，所以需进行调整减小平均绝对误差。

如图2所示，平均相对误差为0.0083，平均绝对误差为66.436。

S3：计算出井下各个泵房涌水量与总涌水量的占比情况，各个水仓的日液位变化，因为每日涌水速率的变化率也会影响涌水预测值，所以通过给定一个涌水速率变化率一个系数，对预测的第一次明日涌水量进行矫正，作为第二次明日涌水量预测值，在本发明中系数取值为0.02，获得第二次明日涌水量预测值如表5所示。

表5第二次明日涌水量预测数据

通过考虑涌水速率变化率的影响，对预测出来的第一次明日涌水量进行矫正，根据表5的第二次明日涌水量预测数据显示，由于近三天的涌水速率变化率过小，所以通过矫正得到的第二次明日预测值没有变化，在各阶段的时间中，涌水速率发生变化时，可对涌水量预测进行矫正。

梯度提升树回归算法GBDT(Gradient Boosting Decision Tree)是一种常用的机器学习算法，主要用于回归和分类问题。它通过组合多个决策树模型来建立一个强大的预测模型。GBDT具有很高的预测准确性，在各种机器学习竞赛和实际应用中都表现出色。对于缺失值和异常值具有一定的鲁棒性，不容易被极端值影响，能够处理高维稀疏数据。可以通过特征重要性评估，选择对预测结果影响较大的特征，提高模型的解释能力和泛化能力。GBDT中的每个决策树都是基于可解释的特征进行构建的，因此GBDT具有很好的可解释性，能够帮助我们理解数据和预测结果之间的关系。

本发明使用两次梯度提升树回归算法进行涌水量预测。

使用梯度提升树回归算法进行第一次涌水量预测，包括以下步骤：

分别选取井下中央泵房、井下一号采区清水泵房、井下一号采区污水泵房、井下二号采区污水泵房的涌水量进行单独分析；

以选取的每个泵房的每日排水量、涌水速率、日水仓液位变化、前一日涌水量、第一次明日涌水量预测值、第二次明日涌水量预测值为自变量因子；

以各个泵房的涌水量预测为因变量因子；

使用梯度提升树回归算法分别预测选取的每个泵房的涌水量预测值，分别作为井下中央泵房涌水量预测值、采区一号采区清水泵房涌水量预测值、采区一号采区污水泵房涌水量预测值、采区二号采区泵房污水泵涌水量预测值；

其中，梯度提升树回归算法包括以下计算步骤：

初始化

迭代建立m棵决策树，以训练数据集

对r

利用线性搜索估计叶节点区域的值

G

以各个泵房的每日的涌水量数据作为训练集训练模型，其中，迭代步长为0.1，最大深度为5，最大份数32，决策树个数为20，预测中央泵房示例结果，如表6所示。

表6中央泵房涌水量预测数据

如表6所示，为中央泵房涌水量预测数据的部分数据，通过使用梯度提升树回归算法对中央泵房的涌水量预测，其中，真实值为中央泵房实际涌水量值，预测值为中央泵房预测出来的值，均方根误差为41.6439，平均绝对误差为32.1247，平均相对误差0.0695，通过分析可知，平均相对误差偏大，平均绝对误差较相比单一指数平滑法预测有所减小。

如图3所示，平均相对误差为0.0695，平均绝对误差为32.1247。

S4：使用梯度提升树回归算法进行第二次涌水量预测，包括以下步骤：

以井下各个泵房的每日排水量、涌水速率、日水仓液位变化、前一日涌水量、第一次明日涌水量预测值、第二次明日涌水量预测值以及中央泵房涌水量预测值、井下一号采区清水泵涌水量预测值、井下一号采区污水泵涌水量预测值、井下二号采区泵房污水泵涌水量预测值作为自变量因子；

以涌水量预测值为因变量因子；

使用梯度提升树回归算法对井下各个泵房的涌水量进行二次涌水量预测，获得井下各个泵房的最终涌水量预测值。

以井下各个泵房的每日排水量、涌水速率、日水仓液位变化、第一次明日涌水量预测值、第二次明日涌水量预测值、前一日涌水量以及获得选取的各个泵房涌水量预测值数据作为训练集，使用梯度提升数回归算法预测出最终涌水量，其中迭代步长为0.1，最大深度为5，最大份数32，决策树个数为20，预测示例结果如表7所示。

表7最终的井下涌水量预测数据

如表7所示，为井下各个泵房的最终涌水量预测数据的部分数据，通过叠加使用梯度提升树回归算法，预测出最终的全矿涌水量，其中，最终预测出来的均方根误差为39.2095，平均相对误差为0.0039，平均绝对误差为30.6189。

相比指数平滑预测出来的误差与单一使用梯度提升树回归算法，叠加使用后的算法预测的结果均方根误差、平均相对误差、平均绝对误差都有明显的减小。

如图4所示，平均相对误差为0.0039，平均绝对误差为30.6189。

S5：根据井下各个泵房的最终涌水量预测值与各个泵房的各水泵性能参数，分析计算获得每个泵房的水泵开泵时长预测值及开泵耗能预测值，预测结果如表8-10所示。

表8中央泵房水泵开泵时长预测值及开泵耗能预测值

如表8所示，为中央泵房水泵开泵时长预测值及开泵耗能预测值的部分数据，通过最终预测出来的中央泵房的涌水量预测值、中央泵房的涌水量与总涌水量的占比以及中央泵房各水泵性能，对中央泵房各个水泵的开泵时长预测与耗电情况进行评估。

表9采区一号泵房污水泵开泵时长预测值及开泵耗能预测值

如表9所示，为采区一号泵房污水泵开泵时长预测值及开泵耗能预测值的部分数据，通过最终预测出来的涌水量、采区一号泵房污水泵涌水占比以及采区一号泵房污水泵各水泵性能，对采区一号泵房污水泵各个水泵的开泵时长预测与耗电情况进行评估。

如表10所示，为采区一号泵房清水泵开泵时长预测值及开泵耗能预测值的部分数据，通过最终预测出来的涌水量、采区一号泵房清水泵涌水占比以及采区一号泵房清水泵各水泵性能，对采区一号泵房清水泵各个水泵的开泵时长预测与耗电情况进行评估。

表10采区一号泵房清水泵开泵时长预测值及开泵耗能预测值

根据上述预测数据分析可知，煤矿井下的大部分涌水都汇入到各个井下泵房中，所以通过自动化系统采集的井下泵房的每日排水量、涌水速率以及水仓液位参数能够代替需要专业人员定期手动的测量采集含水层厚、引用半径、影响半径、渗透系数、采掘深度等参数来计算涌水量。

采用一次指数平滑法和两次梯度提升树回归算法，增加多个自变量参数预测涌水量，提高了涌水量的预测精准度，使得R方值为0.9995。通过对明日涌水量的精准预测，结合水泵性能参数，对各个泵房水泵的开泵时长与能耗进行预测，进而能够准确实现对每日水泵的开泵计划进行评估。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。

另外，除非另有说明，否则本发明使用的所有技术和科学术语具有本发明所属领域的常规技术人员通常理解的相同含义。本说明书中提到的所有文献通过引用并入，用以公开和描述与所述文献相关的方法。在与任何并入的文献冲突时，以本说明书的内容为准。