基于神经网络模糊多模型软测量法的球磨机存煤量测量方法

基于神经网络模糊多模型软测量法的球磨机存煤量测量方法

近年来，正确的煤炭储存能力（以下简称差动）的存煤量测量一直是一个尚未解决的问题。进行准确的存煤量测量不仅有助于提高球磨机的控制质量，而且有助于降低能源消耗，提高效率。这对球磨机的安全运行和经济效益具有重要的现实意义。国外大多采用直接测量存煤量的方法,如在球磨机轴颈加秤、使用称重给煤机或在原煤斗和给煤机间加装计量设备等.这些方法的投资、改装和维护的工作量大,而且一旦钢球装载量或煤的湿度发生变化时会直接影响存煤量的测量.国内多采用间接测量的方法,如差压法、功率电流法、音频法、振动法、物位法、油压法、气压差动法和超声波法等,这些通过单一信号来表征球磨机内存煤量的方法均由于一定的局限性而无法实现准确有效的测量.目前,球磨机存煤量测量的发展趋势主要有两种:一是充分利用信号分析技术,对特征因素进行深入挖掘并找出所包含的潜在信息;二是多信号融合软测量,通过多信号联合建立更准确的软测量模型.近年来第二种趋势正为越来越多的人所青睐,成为发展的主流.1球磨机的立过程为了深入研究球磨机存煤量测量的问题,需先建立球磨机的机理模型.在模型建立过程中,作如下假设:不深入探究磨筒内复杂的磨煤过程,仅从总体上反映球磨机的特性;不考虑制粉系统内散热损失,漏风集中在球磨机入口;磨出口温度是球磨机内煤、风混合物温度;将磨筒、下降干燥管和回粉管合并在一起考虑.1.1钢球磨损的最佳工艺参数钢球磨损和钢球装载量按下式求得mr=Ρ(t)m+S(t)x(1)mr=P(t)m+S(t)x(1)式中:mr为磨球经过磨损后的剩球质量配比向量,mr=[mr1,mr2,…,mrq]T;m为磨球的初始质量配比向量,m=[m1,m2,…,mq]T;x为球磨机实际运行工况的参数向量,x=[ξ,η,δ,Bg,w(M)]T,其中ξ为煤的可磨度,η为磨球的耐磨性系数,δ为煤的细度,Bg为给煤量,w(M)为煤的含水量;P(t)为磨球磨损分布概率矩阵;S(t)为不同工况下运行参数对磨球磨损的影响系数矩阵.当不考虑工况运行参数对磨损的影响时,钢球磨损的经验公式为:mi-mr,i=990.7e-0.0763√(mi/π)(2)mi−mr,i=990.7e−0.076(mi/π)√3(2)式中:mi为i尺寸级别磨球的初始装球质量,kg;mr,i为i尺寸级别磨球经过磨损后的剩球质量,kg;i=1,2,…,q.由于存在磨损,钢球装载量Mg,q随着时间和工况的改变而改变,明显地反映了球磨机的时变性.Μg‚q=q∑i=1mr,i(3)1.2煤调机产粉量dwm/d球磨机进出口质量平衡方程为(Bg-Bm)/3.6=dWm/dτ(4)式中:Bg为给煤机的给煤量,t/h;Bm为球磨机出粉量,t/h;Wm为球磨机内的存煤量,kg.1.3机械热-磨煤过程热容和tg、tf、ctf的计算球磨机进出口能量平衡方程为ddτ[(CmWm+CgqΜg‚q)Τ]=∑Qin-∑Qout(5)式中:∑Qin、∑Qout分别为进入和流出球磨机的热量,kJ/s;Cm、Cgq为煤粉比热容和钢球比热容,J/(g·K);T为球磨机出口温度,℃.∑Qin=(CgBgtg/3.6)+CzGztz+ClGltl+Q0(6)式中:Cg、Cz、Cl分别为原煤比热容、干燥剂比热容和冷风比热容,J/(g·K);Gz、Gl分别为干燥剂流量和冷风流量,kg/s;tg、tz、tl分别为原煤温度、干燥剂温度和冷风温度,℃;Q0为磨煤过程中产生的机械热-磨煤散热,kJ/s.∑Qout=CmBmΤ/3.6+Ctf(Gz+Glf)Τ+BgΔw(2491+1.884Τ-4.19tg)/3.6(7)式中:Ctf为通风比热容,J/(g·K);Glf为漏风量,kg/s;Δw为磨煤过程中煤的水分蒸发量.1.4旋转磨损机的冷却公式球磨机的出口温度方程为:Τ=1CmWm+CgqΜg‚q∫(∑Qin-∑Qout)dτ(8)1.5球磨机出粉量计算建立球磨机的出粉量模型如下Bm=B0kΤkWmktfkcf(9)式中:B0是根据球磨机型号计算的出粉量,t/h;kT、kWm、ktf、kcf分别为球磨机出口温度、球磨机内存煤量、球磨机内通风量及粗粉分离器折向门开度变化对球磨机出粉量的影响系数.1.6球磨机阻力系数fm建立球磨机的进出口压差-流量模型如下:Δpm=fm(1+0.8um)G2tf/ks(10)式中:fm为球磨机的阻力系数;um为通过球磨机的煤粉质量浓度;ks为筒内空气流通面积对磨煤机进出口压差Δpm的影响系数.1.7球磨机出口压力为便于建模将粗粉分离器阻力、细粉分离器阻力、煤粉提升阻力视为常数,且将热风门入口压力、排粉机出口压力也视为常数,则球磨机的出口压力pm0为:pm0=ζG2tf+pzx-p0(11)式中:ζ为粗粉、细粉分离器阻力和煤粉提升阻力三者的综合阻力系数;pzx为排粉机出口压力,Pa;p0为排粉风机的零位压力,Pa.球磨机入口负压方程:pm=pm0+Δpm=ζG2tf+pzx-p+Δpm=ζG2tf+pzx-p+fm(1+0.8um)G2tf/ks(12)1.8m的2.5m将某DTM350/600型球磨机的相关参数代入式(3)、式(4)、式(8)、式(9)、式(10)和式(12),建立的机理模型如下:{Μgq=q∑i=1(mi-990.7e-0.0763√(mi/π))Wm=∫{0.278Bg-[39.4071+e-(Τ-530)2-19.703]×(3.333×10-4Wm-2.778×10-8W2m)×1.5Gz+0.024Bg1.05Gz+0.017Bg+10.3}dτΤ=∫[(376.253Gz-1.293GzΤ-39.18Bg-0.047BgΤ-0.467BmΤ)/(0.502Μgq+1.533Wm)]dτΔpm=[0.196×(1.5Gz+0.024Bg)2+0.086Bm(1.5Gz+0.024Bg)]/(0.55-3.54×10-6Μgq-2.606×10-5Wm)pm=5.823×(1.5Gz+0.024Bg)2-8740+[0.196×(1.5Gz+0.024Bg)2+0.086Bm(1.5Gz+0.024Bg)]/(0.55-3.54×10-6Μgq-2.606×10-5Wm)(13)2模糊聚类与径向基函数rbf神经网络为了更好地解决球磨机存煤量测量的问题,首先通过灰色关联分析选取合适的辅助变量,然后采用神经网络模糊多模型软测量的方法,即将模糊C-均值(FuzzyC-Means,FCM)聚类与径向基函数RBF神经网络多模型相结合,以实现存煤量的准确测量.2.1邓氏关联度算法灰色关联分析的基本思想是根据参考序列和比较序列变化曲线间的相似程度来判断各因素之间的关联程度,即灰关联度.它具有简单直观,对样本量的大小要求不高,定量分析与定性分析的结果一般能够吻合等优点.尤其在分析包含未知信息的灰色系统时,其优势更加明显.与其他的灰关联度相比,邓氏关联度应用最早,算法相对成熟,使用范围也较为广泛.邓氏关联度算法简要描述:首先确定分析的序列矩阵,包括参考序列和比较序列;其次对变量序列进行无量纲化;再次求绝对差值序列(Δ1,Δ2,…,Δn)、最大差Δmax和最小差Δmin;最后求关联系数ξi及灰关联度γi,i=1,2,…,n.ξi(k)=Δmin+ρΔmaxΔi(k)+ρΔmax(14)γi=1mm∑k=1ξi(k)(15)2.2多模型多模型软测量方法神经网络多模型理论的基本思想是分别建立多个独立的神经网络模型,然后将各个模型的输出以某种方式进行联接作为系统最终输出,从而改进模型的预测能力.FCM聚类算法是根据指定的聚类数,由样本与聚类中心的距离关系进行分类.将二者相结合,就形成了神经网络模糊多模型软测量的方法.考虑一个多输入单输出(MISO)的系统,设样本的输入数据矩阵为X=[x1,x2,…,xn]T,样本的输出数据矩阵为Y=[y1,y2,…,yn]T.将X分为训练样本集Xtrain和检验样本集Xtest.神经网络模糊多模型软测量方法见图1,其思路如下:(1)采用FCM聚类算法对训练样本集Xtrain中的数据进行聚类,生成隶属度矩阵U.(2)将训练样本集Xtrain作为神经网络的输入,由(1)得到的隶属度矩阵U作为输出,通过这些输入输出数据训练一个神经网络NN0,从而可以应用此网络来构造样本的隶属度,即当输入检验样本集Xtest时,就可以生成对应的隶属度矩阵U′=[U′1,U′2,…,U′c].(3)根据最大隶属度原则将训练样本分成C类样本数据组.训练样本集Xtrain中的任一样本xk的隶属度矩阵为uk=[u1k,u2k,…,uck],其中c为聚类数.令ujk=max{uik},i=1,2,…,c,ujk为训练样本xk对第j个聚类的隶属度.根据最大隶属度原则,将xk归于第j个聚类.(4)对于由(3)得到的C类样本数据组,分别针对每类数据组及其对应的输出建立RBF神经网络模型,则生成了C个神经网络NN1,NN2,…,NNc.当输入检验样本集Xtest时,所有神经网络输出的集合为Y′={Y′1,Y′2,…,Y′c}.(5)对于任一输入样本xk∈X,其多模型软测量输出为:yk=c∑i=1uikyik(16)式中:uik为训练样本xk对第i个聚类的隶属度;yik为xk输入到第i个模型后得到的输出;k=1,2,…,n,n为所含的样本数.对于检验样本集Xtest,其多模型软测量输出为:Y=c∑i=1(U′i)ΤY′i(17)3模拟研究3.1钢球磨煤机主要运行状态的关联度分析辅助变量过多会导致过学习现象,辅助变量过少会导致欠学习现象,因此选择合适的辅助变量可以大大提高软测量模型的测量精度.原始可测变量数目很多,所以应当根据灵敏、准确以及维数尽量少的原则,对原始变量进行选择以达到降维目的.首先进行粗选,利用机理分析的方法选择一些对被测变量影响较大的相关变量作为预选变量.它们与存煤量之间的关系见表1.然后进行细选,应用适当算法从预选变量集中选取辅助变量,剔除冗余的变量.钢球磨煤机属于部分信息可知的灰色系统,故采用灰色关联分析中的邓氏关联度方法来计算各预选变量与存煤量之间的关联度.除磨电流之外的7个预选变量的关联度分析结果见表2.在指定满意阈值为0.355的前提下,预选变量集中的进出口压差、入口负压、入口风温、出口温度及通风量这5个变量,与存煤量之间的关联度满足要求,因此,选取它们作为软测量的辅助变量.为分析磨电流与存煤量之间的关联度,笔者根据球磨机的运行状态将其划分为四个工作区域:低负荷区、习惯工作区、优化工作区和堵磨区(下文3.2节中将会详细阐述).由于磨电流在低负荷区和习惯工作区内变化不明显,而在优化工作区和堵磨区内变化明显,所以在分析磨电流与存煤量之间的关联度时分为两个阶段讨论,分析结果见表3.在指定满意阈值为0.355的前提下,选取第二阶段(优化工作区、堵磨区)中的磨电流作为软测量的辅助变量,这样可以显著提高存煤量的检测灵敏度.综上所述,当球磨机运行在低负荷区和习惯工作区时,选取了5个辅助变量;当其运行在优化工作区和堵磨区时,选取了6个辅助变量.而且这些辅助变量均为工业现场易采集的信号,不需要加装检测振动或噪声等信号的特殊变送器,因此具有更广泛的实用性.3.2训练样本的选择根据式(13)所建立的机理模型选取样本,并进行信号预处理(包括滤波和归一化处理),最终得到900个样本.每个样本中含有6个特征参数,分别为进出口压差、入口负压、入口风温、出口温度、通风量及磨电流.从总体样本中选取700个作为训练样本,记作Xtrain;其余的200个作为检验样本,记作Xtest.采用FCM算法对训练样本集Xtrain中的数据进行聚类,生成四个聚类中心,每种聚类所包含的样本个数分别为:286,170,135,109.聚类结果见图2.值得注意的是,FCM聚类时,样本所包含的6个特征参数均作为输入变量参与了仿真试验,但绘图时最多只能显示3个.通过模糊聚类将训练样本划分为四类:低负荷区、习惯工作区、优化工作区和堵磨区.对于不同的聚类,磨煤机存煤量与出力的关系见图3.3.3入口通风压、出口温度及通风量对模糊聚类分成的各类训练样本建立RBF神经网络模型.当构建球磨机运行在低负荷区和习惯工作区内的RBF神经网络模型时,输入节点(辅助变量)分别为:进出口压差、入口负压、入口风温、出口温度及通风量;构建优化工作区和堵磨区内的网络模型时,输入节点分别为:进出口压差、入口负压、入口风温、出口温度、通风量及磨电流.由此得到代表球磨机不同运行状态的4个RBF神经网络,见图4.对于任一样本x,由FCM聚类生成的隶属度为u=[u1,u2,u3,u4],4个不同RBF神经网络输出的集合为{y1,y2,y3,y4},综合此样本的各网络输出及其对应的隶属度,由式(16)得到球磨机存煤量的软测量输出为y=u1y1+u2y2+u3y3+u4y4.神经网络多模型结构见图5.3.4个rbf神经网络模型用训练样本集Xtrain中的700个样本训练一个RBF神经网络模型,见图6.这是一个MI