技术领域
[0001] 本
发明涉及一种利用神经网络预测粉煤覆盖煤堆安全堆储时间的方法，属于煤自燃预测技术领域。
背景技术
[0002] 我国的煤炭资源十分丰富，煤炭产量和消费量均居于世界前列，约占国内一次
能源生产和消费总量的85％以上，但是我国煤炭自燃火灾十分严重。从新疆、宁夏、内蒙古大面积煤田露头火灾和灵武、神府、大同等矿区地面煤场和秦皇岛等码头的煤堆自燃火灾到兖州、义
马、靖远、抚顺等矿区井下采煤区域的
煤层自燃火灾随处可见。煤的自燃是一种非常严重的
自然灾害，严重影响了煤炭的生产，造成了重大的经济和资源的损失。
[0003] 当前对于煤堆自燃倾向性和自然发火期的预测方法主要有现场检测法以及数学模型预测法。现场检验法主要是测量煤堆内部的
温度以及对于煤堆中气体进行浓度分析来对于自燃倾向性进行判断，这种方法具有一定的滞后性不能对煤堆的自燃进行很好的
预防；建立数学模型能够通过总结煤自燃的规律对自燃倾向性进行提前预测，建立数学模型的方式对自燃倾向性进行预测，主要有两种方法，传统多项式回归以及神经网络预测，此类方法多用于判断煤自燃危险等级，无法预测确定的自燃发火期。
发明内容
[0004] 为了解决当前预防煤自燃方法的预测滞后性以及不确定性的问题，本发明提供一种利用神经网络预测粉煤覆盖煤堆安全堆储时间的方法，预防煤堆的自燃。
[0005] 为实现上述目的，本发明采用的技术方案如下：一种利用神经网络预测粉煤覆盖煤堆安全堆储时间的方法，包括以下步骤：
[0006] 将粉煤覆盖厚度、空气流量、堆储时间、煤活化能、煤堆
空隙率设定为煤堆安全堆储时间的影响因素，将放热强度和温度作为评判煤堆自燃倾向性指标；
[0007] 选取不同煤种、不同空隙率的煤堆，按照统一直径堆放在不同空气流量的环境中，并在煤堆的中部以下覆盖不同厚度的粉煤，在不同的堆储时间测得温度，并在煤堆测点
位置取气体样品分析计算得到放热强度，从而获得N组包括粉煤覆盖厚度、空气流量、堆储时间、煤活化能、煤堆空隙率以及放热强度和温度的输入输出一一对应的数据组；
[0008] 将包括粉煤覆盖厚度、空气流量、堆储时间、煤活化能、煤堆空隙率的N组影响因素数据组进行归一化预处理；将其中的n组数据组作为训练样本集，其余N-n组数据组作为测试样本集；
[0009] 构建神经网络模型，将归一化预处理后的影响因素数据组作为神经网络模型的输入，将评判煤堆自燃倾向性指标作为神经网络模型的输出，建立煤堆安全堆储时间影响因素和煤堆自燃倾向性指标的非线性映射；
[0010] 设置神经网络模型的
最小均方误差，将训练样本集作为神经网络模型的输入和预期输出，对所述神经网络模型进行训练，并利用所述神经网络模型进一步设计样本，使之达到学习记忆功能；当所述神经网络模型学习
精度满足最小均方误差后即得到训练好的神经网络模型；
[0011] 调整神经网络中的隐含层
节点数目，将测试样本集输入训练好的神经网络模型，基于最小的训练误差确定神经网络的隐含层节点数目，得到用于预测粉煤覆盖煤堆安全堆储时间的BP神经网络模型；
[0012] 对待测煤堆进行测试，获得包括粉煤覆盖厚度、空气流量、堆储时间、煤活化能、煤堆空隙率的待测煤堆影响因素实测值，将所述待测煤堆影响因素实测值数据组输入BP神经网络模型，获得放热强度和温度的预测值；
[0013] 将放热强度和温度预测值用来计算待测煤堆的最短自燃发火期，得到粉煤覆盖煤堆安全堆储时间的预测结果。
[0014] 上述方法还包括：
[0015] 构建神经网络模型，将归一化预处理后的影响因素数据组作为神经网络模型的输入，将评判煤堆自燃倾向性指标温度与放热强度的归一化数值作为神经网络模型的输出结果，所述输出结果经过反归一化处理，得到放热强度与温度的转化值。
[0016] 所述归一化预处理方法为：
[0017]
[0018] 其中，xk为归一化后数据，x为原始数据，xmin、xmax分别为同一影响因素数据组中的最小值和最大值；所述数据包括粉煤覆盖厚度、煤堆空隙率、煤活化能、空气流量、堆储时间中的任一影响因素数据；经过归一化预处理后的数据xk分布在[0,1]区间内。
[0019] 所述反归一化预处理方法为：
[0020] x＝x输出×(xmax-xmin)
[0021] 其中，x为转化值，x输出为神经网络模型输出的归一化值，xmin、xmax分别为同一影响因素数据组中的最小值和最大值。
[0022] 上述方法还包括：
[0023] 气体取样分析，获得包括O2、CO和CO2的浓度，对应堆储时间和温度，得到O2的消耗速率以及CO和CO2的生成速率；
[0024] 通过以下公式计算得到最大放热强度qmax和最小放热强度qmin，并通过qmax和qmin的算术平均值计算得到放热强度q：
[0025]
[0026]
[0027] 其中，qmax、qmin为放热强度的最大值与最小值，单位J/s； 为
氧气的消耗速率，vCO、 分别为CO、CO2的生成速率，单位mol/s；ΔHCO、 分别为生成CO和CO2的平均
焓变，单位J/mol；qa为气体化学
吸附热，单位J/mol。
[0028] 上述技术方案中，所述最小均方误差选用1×10-5。
[0029] 上述技术方案中，所述粉煤为通过选取煤堆的适量煤样经过
粉碎后获得的粒径小于6mm的
煤粉。
[0030] 上述技术方案中，所述最短自燃发火期通过公式计算：
[0031]
[0032] 其中，ρc为煤体
密度，单位g/cm3；cc为煤体
比热容，单位kJ/(kg·K)；τ为煤自燃时间，单位d；T为温度，单位K；q(T)为放热强度，单位J/s。
[0033] 与
现有技术相比，本发明具有如下有益效果：通过利用BP神经网络对煤堆自燃发火期的预测，能够提前应对煤堆自燃；同时，能够对不同煤种、不同堆放条件的煤堆进行自燃发火期的预测，方法简单、可靠、适用范围广。
附图说明
[0034] 图1为本发明的神经网络模型结构示意图。
[0035] 图2为本发明
实施例一预测温度值和实测温度值对比图。
[0036] 图3为本发明实施例二预测温度值和实测温度值对比图。
[0037] 图4为本发明实施例三预测温度值和实测温度值对比图。
具体实施方式
[0038] 下面结合附图和具体实施例对本发明作进一步详细说明。
[0039] 为了实现对粉煤覆盖的煤堆安全堆储时间的预测，首先通过调研分析识别得到与安全堆储时间相关的因素，将粉煤覆盖厚度、空气流量、堆储时间、煤活化能、煤堆空隙率设定为煤堆安全堆储时间的影响因素，将放热强度与温度作为评判煤堆自燃倾向性指标，并通过放热强度和
温度计算出煤堆最短自燃发火期。
[0040] 采用中国榆家梁(YJL)、瓷窑湾(CYW)和上湾(SW)三种易自燃煤建立煤堆，进行了粉煤覆盖下煤堆低温氧化的模拟试验。首先测试出与煤种相关的数据，比如：
[0041] 煤的空隙率：
[0042]
[0043] 其中，ε为空隙率，单位cm3/g；ρ、ρα分别为煤的堆积密度以及煤的表观密度，单位kg/m3、g/cm3。
[0044] 煤的活化能：
[0045]
[0046] 其中，Ea为活化能，单位J/mol；K为速率常数；T为温度，单位K；R为摩尔气体常数，单位J/(mol*K)。
[0047] 将煤堆按照统一直径堆放，改变粉煤覆盖厚度和空气流量，随着堆储时间的变化，在统一高度和深度位置测定温度，并使用气体测量仪以及
传感器测量此处气体中主要组分的浓度，如O2、CO和CO2等。本实施例选取的是直径为35m且中下部分别覆盖30、70、100cm厚度粉煤的不同煤种不同空隙率的煤堆，堆放于
风速不同(0.08、0.10、0.12cm/s)环境中，选用测点位置为煤堆2m高、1.5m深处的煤氧化层。根据时间变化和气体组分浓度变化，得到O2的消耗速率以及CO和CO2的生成速率，依据如下公式计算煤堆不同堆储时间下的放热强度最大值和最小值：
[0048]
[0049]
[0050] 其中，qmax、qmin为放热强度的最大值与最小值，单位J/s； 为氧气的消耗速率，vCO、 分别为CO、CO2的生成速率，单位mol/s；ΔHCO、 分别为生成CO和CO2的平均焓变，单位J/mol；qa为气体
化学吸附热，单位J/mol。
[0051] 为了能够对不同煤种的煤堆都进行预测，直接通过qmax和qmin求取算术平均值近似认为是煤体的特征放热强度，作为上述数据组中的放热强度。
[0052] 记录下包括粉煤覆盖厚度、煤堆空隙率、煤活化能、空气流量、堆储时间以及放热强度和温度的输入输出一一对应的100组数据。
[0053] 构建如图1所示的包含
输入层、隐含层以及
输出层的神经网络模型，其中选用tansig、purelin作为传递函数，trainscg作为训练函数。将包括粉煤覆盖厚度、空气流量、堆储时间、煤活化能、煤堆空隙率的影响煤低温氧化的因素作为神经网络模型的输入，再将评判煤自燃倾向性的指标(放热强度与温度)作为神经网络模型的输出，建立影响因素与煤自燃倾向性的非线性映射。形成包含5个输入层节点、若干个隐含层节点和2个输出层节点的神经网络模型。
[0054] 然后使用100组数据的前70组数据作为神经网络模型的训练样本集，后30组数据-5作为神经网络模型的测试样本集，建立并获得符合精度(最小均方误差为1×10 )的BP神经网络模型。
[0055] 训练神经网络模型前为了避免数据绝对大小对结果的不良影响，采用如下办法对输入和输出数据进行归一化处理：(实际值-最小值)/(最大值-最小值)＝归一值，即[0056]
[0057] 其中，xk为归一化后数据，x为原始数据(某一影响因素的实际值)，xmin、xmax分别为同一因素原始数据中的最小值和最大值；所述数据包括粉煤覆盖厚度、煤堆空隙率、煤活化能、空气流量、堆储时间中的任一影响因素数据。
[0058] 通过神经网络模型计算得到输出值结果后，再采用反归一化推导出实际值：
[0059] x＝x输出×(xmax-xmin)
[0060] 其中，x为转化值，x输出为神经网络模型输出的归一化值，xmin、xmax分别为同一影响因素数据组中的最小值和最大值。
[0061] 在matlab中，归一化和反归一化分别用premnmx和postmnmx函数实现。
[0062] 最终得到放热强度与温度的实际值。
[0063] 设置神经网络模型的最小均方误差为1×10-5，将训练样本集作为神经网络模型的输入和预期输出，对神经网络模型进行训练，并在训练过程通过重复滚动输出利用神经网络模型进一步设计样本，使之达到学习记忆功能。当神经网络模型学习精度满足最小均方误差后即认为完成了训练，得到训练好的神经网络模型。
[0064] 设定神经网络模型隐含层的初始节点数，将测试样本集输入训练好的神经网络模型对其进行测试，并不断调整隐含层节点数目，重复上述过程，最后基于最小的训练误差确定确定神经网络的隐含层节点数目，得到用于预测粉煤覆盖煤堆安全堆储时间的BP神经网络模型。
[0065] 对待测定的煤堆进行测试，获得包括粉煤覆盖厚度、煤堆空隙率、煤活化能、空气流量、堆储时间的待测煤堆影响因素实测值数据组；将所述待测煤堆影响因素实测值数据组输入BP神经网络模型，预测煤堆的放热强度和温度
[0066] 利用神经网络预测的温度与放热强度的数值，代入下述公式计算煤堆的最短自燃发火期：
[0067]
[0068] 其中，ρc为煤体密度，单位g/cm3；cc为煤体比
热容，单位kJ/(kg·K)；τ为煤自燃时间，单位d；T为温度，单位K；q(T)为放热强度，单位J/s。
[0069] 通过上述分析过程可以看出，煤的放热强度测量和计算是一个复杂的过程，而利用训练好的BP神经网络模型来预测煤堆堆砌过程的温度和放热强度，从而预测出其最短自燃发火期，大大简化了其过程。
[0070] 实施例一
[0071] 选取直径为35m、中下部覆盖70cm厚度的粉煤的上湾煤，煤堆堆储于新鲜风流(0.08cm/s)环境中，煤堆的平均空隙率为0.40，测点选在距离地面2m高、深度1.5m的煤氧化层作为温度与放热强度测点。对现有煤堆输入数据得到输出的温度以及放热强度的数值。表1为放热强度的BP神经网络预测值和文献值对比表。图2为温度的BP神经网络模型预测值和实际煤堆温度实测值对比图。由表1可以看出BP神经网络预测出的放热强度的值在查阅资料得到的放热强度的平均值上下
波动，误差在允许范围内。由图2可以看出煤堆内部测点处的温度的值与BP神经网络的预测值相接近。因此本发明所构建的BP神经网络的值结果较为可靠。
[0072] 表1放热强度的BP神经网络预测值和文献值对比
[0073]
[0074] 注：放热强度文献值查取条件为：新鲜风流(0.08cm/s)中的松散煤体(
烟煤)在低温氧化时的放热强度平均值。
[0075] 根据BP神经网络模型预测的温度值T以及放热强度q求煤堆的最短自燃发火期。代入数据，得到在粉煤覆盖厚度为70cm时，煤堆的最短自燃发火期为109～139天。
[0076] 对于当前煤堆，煤堆覆盖70cm厚的粉煤煤堆堆放120天时温度快速增加，有自燃趋势。证明BP神经网络模型的预测结果可靠。
[0077] 实施例二
[0078] 选取直径为35m、中下部覆盖80cm厚度的粉煤的榆家梁煤，煤堆堆储于新鲜风流(0.08cm/s)环境中，煤堆的平均空隙率为0.45，测点选在距离地面2m高、深度1.5m的煤氧化层作为温度与放热强度测点。对现有煤堆输入数据得到输出的温度以及放热强度的数值。表2为放热强度的BP神经网络预测值和文献值对比表。图3为温度的BP神经网络模型预测值和实际煤堆温度实测值对比图。由表2可以看出BP神经网络预测出的放热强度的值在查阅资料得到的放热强度的平均值上下波动，误差在允许范围内。由图3可以看出煤堆内部测点处的温度的值与BP神经网络的预测值相接近。因此本发明所构建的BP神经网络的值结果较为可靠。
[0079] 表2放热强度的BP神经网络预测值和文献值对比
[0080]
[0081] 注：放热强度文献值查取条件为：新鲜风流(0.08cm/s)中的松散煤体(烟煤)在低温氧化时的放热强度平均值。
[0082] 根据BP神经网络模型预测的温度值T以及放热强度q求煤堆的最短自燃发火期。代入数据，得到在粉煤覆盖厚度为80cm时，煤堆的最短自燃发火期为119～145天。
[0083] 对于当前煤堆，煤堆覆盖80cm厚的粉煤煤堆堆放130天时温度快速增加，有自燃趋势。证明BP神经网络模型的预测结果可靠。
[0084] 实施例三
[0085] 选取直径为35m、中下部覆盖90cm厚度的粉煤的瓷窑湾煤，煤堆堆储于新鲜风流(0.08cm/s)环境中，煤堆的平均空隙率为0.48，测点选在距离地面2m高、深度1.5m的煤氧化层作为温度与放热强度测点。对现有煤堆输入数据得到输出的温度以及放热强度的数值。表3为放热强度的BP神经网络预测值和文献值对比表。图4为温度的BP神经网络模型预测值和实际煤堆温度实测值对比图。由表3可以看出BP神经网络预测出的放热强度的值在查阅资料得到的放热强度的平均值上下波动，误差在允许范围内。由图4可以看出煤堆内部测点处的温度的值与BP神经网络的预测值相接近。因此本发明所构建的BP神经网络的值结果较为可靠。
[0086] 表3放热强度的BP神经网络预测值和文献值对比
[0087]
[0088] 注：放热强度文献值查取条件为：新鲜风流(0.08cm/s)中的松散煤体(烟煤)在低温氧化时的放热强度平均值。
[0089] 根据BP神经网络模型预测的温度值T以及放热强度q求煤堆的最短自燃发火期。代入数据，得到在粉煤覆盖厚度为90cm时，煤堆的最短自燃发火期为129～150天。
[0090] 对于当前煤堆，煤堆覆盖90cm厚的粉煤煤堆堆放140天时温度快速增加，有自燃趋势。证明BP神经网络模型的预测结果可靠。
[0091] 以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何
修改、等同替换、改进等，均应包含在本发明的保护范围之内。
