一种基于极限学习机的水泥磨系统电耗指标预测方法技术方案

本发明专利技术公开一种基于极限学习机的水泥磨系统电耗指标预测方法，首先对相关变量进行收集，通过对水泥磨系统的运行分析，选取与水泥磨电耗相关的八个变量作为输入变量，吨电耗作为输出变量，构建ELM模型的输入输出层；使用IPSO算法对ELM的神经元个数L、输入层权重w及隐含层偏置b进行优化；IPSO算法在参数寻优过程中将模型的均方误差作为适应度函数；使用IPSO优化得到的ELM相关参数结合样本数据完成IPSO‑ELM预测模型的训练，将工业现场实际数据带入到训练好的模型完成水泥磨系统电耗的在线预测；本发明专利技术通过使用样本数据对IPSO‑ELM模型进行训练得到电耗预测模型，将实际水泥生产现场的变量数据输入到训练好的模型中，实现水泥磨电耗指标在线预测。

【技术实现步骤摘要】
一种基于极限学习机的水泥磨系统电耗指标预测方法
本专利技术属于水泥磨系统电耗指标预测领域，特别涉及一种基于极限学习机的水泥磨系统电耗指标预测方法。
技术介绍
水泥工业是我国传统工业，随着国家的发展，大型建筑不断出现，水泥使用量大大增加，水泥行业是我国主要的高耗能、高排放产业。实现水泥磨电耗指标的在线预测，有利于指导水泥磨研磨过程各个生产指标的调度优化，有利于降低水泥磨生产过程电耗，从而能达到节能减排的目的。但是由于水泥研磨系统各个生产指标之间存在滞后性，并且耦合性较高，因此很难通过建立机理模型模拟水泥磨系统的能耗模型；由于不同操作人员现场经验的差异，很难保证水泥磨处于高效，稳定的运行状态，从而导致能耗浪费。另外，在采用OPC技术对水泥工业现场DCS设备的数据提取过程中，由于人员操作不当，易会产生数据异常，如激增等，为后续建立精确地电耗预测模型增加难度。针对上述问题，赵辉等人采用主成分分析法得到影响水泥生产电耗的关键因素，降低了回归预测模型的复杂度，并提出基于改进多元非线性算法的水泥电耗预测模型，但该方法未解决水泥生产过程的时滞问题。谢武明等人将33个指标作为输入建立基于BP神经网络的污水处理厂电耗预测模型，对于电耗预测这类问题求解复杂非线性函数的全局极值，BP神经网络容易陷入局部最优。王昭旭等人将一种隐含时间序列深度信念网络用到水泥生产电耗预测，消除了时变时延对水泥能耗预测的影响，但深度信念网络预训练采用无监督学习，学习过程缓慢，难以进行水泥磨系统电耗在线预测。
技术实现思路
针对上述问题本专利技术提供一种基于极限学习机的水泥磨系统电耗指标预测方法，有效的解决了水泥磨系统各个电耗指标之间的滞后性问题以及多个电耗指标之间强耦合的问题，且相对于其他神经网络模型，极限学习机(ExtremeLearningMachine，ELM)训练速度及计算速度快，能达到实时训练-预测；极限学习机在训练之前随机产生权重w和偏置b，只需确定隐含层神经元的个数及隐含层神经元的能无限微分的激活函数，所以极限学习机的训练过程，即为求解的过程；然而根据经验设置极限学习机隐含层神经元的个数及随机生成权重w和偏置b会使预测准确度不稳定，在预测模型中加入改进粒子群优化算法(ImproveParticleSwarmOptimization，IPSO)对上述参数进行寻优，提高预测的准确度及稳定性，增强泛化能力。为实现上述目的，本专利技术所采用的技术方案是提供一种基于极限学习机的水泥磨系统电耗指标预测方法，其步骤包括：步骤1：确定输入、输出变量；通过对水泥磨系统的运行分析，确定与水泥磨电耗相关的输入变量为喂料量反馈、出磨斗提电流反馈、选粉机转速反馈、喂料斗提电流反馈、入库提升机电流反馈、比表面积、主机电流及选粉机进口压力共八个工业参量，确定输出变量为水泥磨吨电耗；步骤2：提取输入、输出变量；将水泥磨系统设备中的历史数据存入数据库，从数据库中提取所述输入变量和输出变量；步骤3：对PSO算法改进得到IPSO算法；改进的IPSO算法中的位置更新模型为：其中，wi为惯性权重，fiti为当前粒子第i次迭代的适应度值，Pgd为算法的全局最优值，k，k+1分别表示第k次、第k+1次迭代；c1和c2为加速因子，令c1＝c2＝2；r1和r2为分布于[0,1]之间的随机数；Mi为第i个粒子为S维的向量，Mi＝(Mi1,Mi2,…,Mid,…,Mis)T,Mid为第i个粒子在d维空间中的位置；Vi＝(Vi1,Vi2,…,Vid,…,Vis)T为第i个粒子的速度,Vid为第i个粒子在d维空间中的速度；Pi＝(Pi1,Pi2,…,Pid,…,Pis)T为第i个粒子的个体极值，Pid为第i个粒子在d维空间的个体极值；Pg＝(Pg1,Pg2,…,Pgd,…,Pgs)T为种群全局极值，Pgd为种群g在d维空间的全局极值；步骤4：利用IPSO算法及所述提取的输入、输出变量值对ELM模型参数进行寻优，得到ELM模型的寻优变量：神经元个数L、输入层权重w和隐含层偏置b；步骤5：利用所述寻优变量及所述提取的输入、输出变量值对ELM预测模型进行训练，得到输出权重步骤6：构建基于IPSO-ELM的水泥磨系统电耗指标预测模型为：其中，j＝1，2，…，L，L为隐含层神经元的个数，为输出权重，为激活函数，Wj为隐含层第j个神经元的输入权重矩阵，X＝[x1,x2,…xi,…,x7,x8]为输入变量时间序列数据集，xi＝[xi1,xi2,…,xit]i＝1,2…8，t为时间序列的长度，bj为隐含层第j个神经元的偏置；步骤7：进行在线预测：将实时输入变量数据输入基于IPSO-ELM的水泥磨系统电耗指标预测模型，得到实时的水泥磨吨电耗数值，实现水泥磨系统电耗指标的在线预测。进一步地，所述步骤4中寻优过程如下：步骤41：对种群进行初始化，设置粒子群的初始位置和初始速度；步骤42：根据适应度函数，计算每个粒子的适应度值；步骤43：对每个粒子进行位置比较，即对经过的位置和上一次的位置进行比较，如果适应度低于上一次的位置，就取代上一次的位置；步骤44：对个体最优位置和群体最优位置进行比较，如优于群体最优位置，就取代群体最优位置；步骤45：根据公式更新粒子速度和位置；步骤46：判断是否满足终止条件，若满足则停止运算，否则返回第二步重新开始。进一步地，所述步骤5中模型训练的具体步骤如下：首先利用构建好的输入输出层数据结合IPSO-ELM算法，进行IPSO-ELM水泥磨系统电耗预测模型的训练：对N个任意样本(Xi,ti)，其中，Xi＝[xi1,xi2,xi3,…,xin]T，ti＝[ti1,ti2,ti3,…,tim]T∈Rm，则存在L个隐含层节点的ELM预测模型表示为：其中，oi是ELM预测模型的输出，n为输入变量的个数，m为输出变量的个数，βj为隐含层第j个神经元的输出权重，为激活函数，L为隐含层神经元的个数，Wj为隐含层第j个神经元的输入权重矩阵，Xi＝[xi1,xi2,xi3,…,xin]T为输入变量值，bj为隐含层第j个神经元的偏置；N为样本总数；ELM网络训练的目标就是使预测值与真实值之间的误差最小，表达式为：即存在Wj，bj，βj使得：上述表达式的矩阵表达形式为：Hβ＝T(5)其中：其中，H的第i列是对输入Xi第i个隐藏节点的输出，其中T表示样本的实际值，神经元个数L、输入权重w和隐含层偏置b由IPSO算法确定；训练基于极限学习机的电耗预测模型即寻找系统的最小二乘解求得输出权值矩阵：称为H矩阵的摩尔-彭罗斯广义逆，并且该等式对应的解具有唯一性，使其训练误差达到最小。进一步地，所述步骤6中，t＝50，即输入数据为当前时点之前50分钟的输入变量数据构成的时间序本文档来自技高网...

【技术保护点】
1.一种基于极限学习机的水泥磨系统电耗指标预测方法，其特征在于，其步骤包括：/n步骤1：确定输入、输出变量；/n通过对水泥磨系统的运行分析，确定与水泥磨电耗相关的输入变量为喂料量反馈、出磨斗提电流反馈、选粉机转速反馈、喂料斗提电流反馈、入库提升机电流反馈、比表面积、主机电流及选粉机进口压力共八个工业参量，确定输出变量为水泥磨吨电耗；/n步骤2：提取输入、输出变量；/n将水泥磨系统设备中的历史数据存入数据库，从数据库中提取所述输入变量和输出变量；/n步骤3：对PSO算法改进得到IPSO算法；/n改进的IPSO算法中的位置更新模型为：/n

【技术特征摘要】
1.一种基于极限学习机的水泥磨系统电耗指标预测方法，其特征在于，其步骤包括：
步骤1：确定输入、输出变量；
通过对水泥磨系统的运行分析，确定与水泥磨电耗相关的输入变量为喂料量反馈、出磨斗提电流反馈、选粉机转速反馈、喂料斗提电流反馈、入库提升机电流反馈、比表面积、主机电流及选粉机进口压力共八个工业参量，确定输出变量为水泥磨吨电耗；
步骤2：提取输入、输出变量；
将水泥磨系统设备中的历史数据存入数据库，从数据库中提取所述输入变量和输出变量；
步骤3：对PSO算法改进得到IPSO算法；
改进的IPSO算法中的位置更新模型为：









其中，wi为惯性权重，fiti为当前粒子第i次迭代的适应度值，Pgd为算法的全局最优值，k，k+1分别表示第k次、第k+1次迭代；c1和c2为加速因子，令c1＝c2＝2；r1和r2为分布于[0,1]之间的随机数；Mi为第i个粒子为S维的向量，Mi＝(Mi1,Mi2,…,Mid,…,Mis)T,Mid为第i个粒子在d维空间中的位置；Vi＝(Vi1,Vi2,…,Vid,…,Vis)T为第i个粒子的速度,Vid为第i个粒子在d维空间中的速度；Pi＝(Pi1,Pi2,…,Pid,…,Pis)T为第i个粒子的个体极值，Pid为第i个粒子在d维空间的个体极值；Pg＝(Pg1,Pg2,…,Pgd,…,Pgs)T为种群全局极值，Pgd为种群g在d维空间的全局极值；
步骤4：利用IPSO算法及所述提取的输入、输出变量值对ELM模型参数进行寻优，得到ELM模型的寻优变量：神经元个数L、输入层权重w和隐含层偏置b；
步骤5：利用所述寻优变量及所述提取的输入、输出变量值对ELM预测模型进行训练，得到输出权重
步骤6：构建基于IPSO-ELM的水泥磨系统电耗指标预测模型为：



其中，j＝1，2，…，L，L为隐含层神经元的个数，为输出权重，为激活函数，Wj为隐含层第j个神经元的输入权重矩阵，X＝[x1,x2,…xi,…,x7,x8]为输入变量时间序列数据集，xi＝[xi1,xi2,…,xit]i＝1,2…8，t为时间序列的长度，bj为隐含层第j个神经元的偏置；
步骤7：进行在线预测；
将实时输入变量数据输入基于IPSO-ELM的水泥磨系统电耗指标预测模型，得到实时的水泥磨吨电耗数值，实现水泥磨系统电耗指标的在线预测。


2.根据权利要求1所述的基于极限学习机的水泥磨系统电耗指标预测方法，其特征在于，所述步骤4中寻优过程如下：

【专利技术属性】
技术研发人员：郝晓辰，李东栩，张志鹏，赵彦涛，冀亚坤，徐清泉，
申请(专利权)人：燕山大学，
类型：发明
国别省市：河北;13