发明名称 一种基于数据的冶金企业高炉煤气动态预测方法
摘要 一种基于数据的冶金企业高炉煤气动态预测方法。首先利用经验模态分解的消噪方法对已知的历史数据进行除噪处理,然后利用基于回声状态网络的时间序列预测方法对高炉煤气系统的产消量进行动态预测;其次在回声状态网络中引入贝叶斯推理方法建立各主要煤气用户煤气流量与煤气柜柜位之间的关系模型用来预测高炉煤气柜柜位在当前时刻后指定时间长度内的变化;本发明能够准确地预测高炉煤气系统的产消量和煤气柜柜位变化,为现场煤气平衡调度人员提供在线决策支持。
申请公布号 CN102147273B 申请公布日期 2012.11.28
申请号 CN201010103317.X 申请日期 2010.01.29
申请人 大连理工大学;上海宝信软件股份有限公司 发明人 赵珺;王伟;盛春阳;刘颖;张晓平;从力群;冯为民;吴毅平;沈兵
分类号 G01F1/66(2006.01)I;G01F23/296(2006.01)I 主分类号 G01F1/66(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 梅洪玉
主权项 1.一种基于数据的冶金企业高炉煤气动态预测方法,其特征在于具体步骤如下:(1)从冶金企业现场的实时数据库读取预测所需的高炉煤气系统各用户煤气流量数据和煤气柜柜位数据;(2)采用经验模态分解的消噪方法对高炉煤气系统各用户煤气流量数据进行消噪处理①记待处理的信号为u,找出信号u的所有局部极大、极小值点,用三次样条插值函数拟合出信号的上包络线和下包络线,求出两包络线的平均值m<sub>1</sub>,将u与m<sub>1</sub>之差记为h<sub>1</sub>,h<sub>1</sub>=u-m<sub>1</sub>,判断h<sub>1</sub>是否为IMF,若不是,则重复上述分解过程对其进一步分解;若是,第k步的结果h<sub>1k</sub>写成h<sub>1(k-1)</sub>-m<sub>1k</sub>=h<sub>1k</sub>,定义h<sub>1k</sub>=c<sub>1</sub>,则c<sub>1</sub>即为从原始数据中得到的第一个IMF;②计算r<sub>1</sub>=u-c<sub>1</sub>,将r<sub>1</sub>作为分解数据,重复上述分解过程,当剩余分量小于预定值或为单调<img file="FDA00001774850200011.GIF" wi="1896" he="205" />尺度信号的标准差;④计算第i个尺度信号相邻信号的标准差<img file="FDA00001774850200012.GIF" wi="61" he="49" />若<img file="FDA00001774850200013.GIF" wi="135" he="49" />则采用低通滤波器去噪;⑤将消噪后的IMF分量和其余的n-k个分量重构,得到消噪后的时间序列信号;(3)高炉煤气系统产消量预测①利用消噪后的数据构建训练样本集:将先前m个消噪后的历史数据[u(t),u(t+1),...,u(t+m-1),t=1,2,...,n-m]作为输入样本,用来预测第m+1个数据样本,训练样本集表示为S={(u(k),y(k))|k=1,2,…,n-m}<img file="FDA00001774850200014.GIF" wi="916" he="298" /><img file="FDA00001774850200015.GIF" wi="421" he="288" />其中u(k)∈R<sup>m</sup>,y(k)∈R分别为输入和输出数据;②建立基于回声状态网络的时间序列预测模型:<img file="FDA00001774850200016.GIF" wi="898" he="145" />其中,f<sup>out</sup>为线性函数;x(k)为k时刻DR的状态变量,初始值随机给定;W<sup>in</sup>为输入权值矩阵;W为DR内部神经元的连接矩阵,为使DR具有动态记忆能力,W保持1%~5%的稀疏连接,且谱半径设为0.75;W<sup>out</sup>为输出权值矩阵;③求解回声状态网络模型的权值:W<sup>in</sup>、W两者在网络学习前构建,而W<sup>out</sup>是在网络学习后,通过搜集到的状态矩阵M,M由u(k+1),x(k+1),y(k)组成,和搜集到的给定输出信号T,T由y(k+1)组成,利用公式(W<sup>out</sup>)<sup>T</sup>=M<sup>-1</sup>T求得;④预测高炉煤气系统各用户煤气流量:利用当前时间点前m个时间点的指定用户煤气流量值作为输入,预测当前时间点后n-m个时间点的指定用户煤气流量值,m和n的值由调度人员来指定;⑤预测高炉煤气系统在n-m个时间点内的产消量:高炉煤气在n-m个时间点内的产生量等于所有发生源用户在此n-m个时间点内的煤气流量之和,高炉煤气在n-m个时间点内的消耗量等于所有消耗源用户在此n-m个时间点内的煤气流量之和;(4)高炉煤气柜柜位预测建模①对从数据库中获得的高炉煤气系统各用户煤气流量数据和高炉煤气柜柜位数据进行量纲统一和归一化处理;②构建柜位预测训练样本集:利用影响柜位变化的P个高炉煤气系统用户的煤气流量历史数据和高炉煤气柜柜位历史数据构造模型的训练样本集<img file="FDA00001774850200021.GIF" wi="1271" he="286" /><img file="FDA00001774850200022.GIF" wi="399" he="282" />其中u(k)∈R<sup>P+1</sup>,y(k)∈R分别为输入和输出数据,u<sub>j</sub>(t)表示影响柜位变化的P个用户中的第j个在t时刻的煤气流量数据,H(t)表示t时刻高炉煤气柜的柜位数据,H(t+1)表示t+1时刻高炉煤气柜的柜位数据;③建立基于回声状态网络的柜位预测模型:<img file="FDA00001774850200023.GIF" wi="911" he="233" />其中,f<sup>out</sup>为线性函数;x(k)为k时刻DR的状态变量,初始值随机给定;y(x;W<sup>out</sup>)为网络输出值;D为实际输出值;ε为均值为0,精度为β的高斯白噪声;W<sup>in</sup>为输入权值矩阵;W为DR内部神经元的连接矩阵,为使DR具有动态记忆能力,W保持1%~5%的稀疏连接,谱半径选为0.75;W<sup>out</sup>为输出权值矩阵;④利用贝叶斯推理方法求解柜位预测模型的权值:W<sup>in</sup>、W在网络学习前随机构建,而W<sup>out</sup>是通过在回声状态网络中引入贝叶斯方法求得;因为ε是均值为0,精度为β的高斯白噪声,ε=D-y(x;W<sup>out</sup>),所以D服从均值为y(x;W<sup>out</sup>),精度为β的高斯分布,有<img file="FDA00001774850200031.GIF" wi="681" he="129" />引入贝叶斯原理<img file="FDA00001774850200032.GIF" wi="727" he="134" />用一个很宽的分布近似先验分布,于是有W<sup>out</sup>的先验分布<img file="FDA00001774850200033.GIF" wi="672" he="137" />由贝叶斯规则p(D)=∫p(D|W<sup>out</sup>)p(W<sup>out</sup>)dW<sup>out</sup>,将p(D|W<sup>out</sup>),p(D)和p(W<sup>out</sup>)的表达代入贝叶斯公式中得到<img file="FDA00001774850200034.GIF" wi="1232" he="129" />由于Z<sub>M</sub>与W<sup>out</sup>无关,所以将求解后验分布p(W<sup>out</sup>|D)取最大值时的<img file="FDA00001774850200035.GIF" wi="105" he="66" />转化为求解M(W<sup>out</sup>)取最小值时的<img file="FDA00001774850200036.GIF" wi="127" he="66" /><img file="FDA00001774850200037.GIF" wi="105" he="66" />为所求权值;将步骤(3)预测得到的高炉煤气系统各用户的n-m个煤气流量数据和从数据库中获得的高炉煤气柜柜位历史数据作为步骤(4)所建立的柜位预测模型的输入,预测当前时间点后n-m个时间点对应的柜位值。 
地址 116024 辽宁省大连市甘井子区凌工路2号