荧光显微图像复原有着很多重要的应用,例如,天文成像,电子显微镜成像,单光子发射计算机断层成像术和正电子发射断针对采用机理分析方法得到的煤发热量存在响应速度迟缓、变负荷时动态误差大的问题,提出一种基于状态观测器的动态补偿方法。依据单元机组负荷–压力模型,由汽包压力和汽轮机一级压力构造锅炉总热量信号,采用一高阶惯性环节描述锅炉给煤量与总热量信号之间的动态关系。针对此高阶惯性环节设计状态观测器,得到具有较小惯性的中间状态信号,利用此信号除以给煤量得到煤发热量信号。实验结果表明:动态补偿后的煤发热量具有抗扰动性好、响应速度快的优点。 The coal calorific value calculated by using the method of mechanism analysis has slow response speed and big dynamic error under variable load conditions. In order to solve this problem, a dy-namic compensation method based on state observer is proposed. Based on the load-pressure model of unit, the total heat input into the boiler is constructed by the drum pressure and the tur-bine primary pressure. The dynamic relationship between the boiler fuel flow and the total heat is described by a high order inertia link. By designing a state observer for the high order inertia link, the intermediate state signal with smaller inertia is obtained, which is divided by the fuel flow to gain the signal of the coal calorific value. The signal of coal calorific value has small dynamic error under variable load conditions, Experiments show that: the signal of coal calorific value has better behaviors on disturbance restraint and dynamic response.
武现聪1,田亮2,杜文霞1
1河北师范大学职业技术学院,河北 石家庄
2华北电力大学计算机与控制工程学院,河北 保定
收稿日期:2016年7月26日;录用日期:2016年8月14日;发布日期:2016年8月17日
针对采用机理分析方法得到的煤发热量存在响应速度迟缓、变负荷时动态误差大的问题,提出一种基于状态观测器的动态补偿方法。依据单元机组负荷–压力模型,由汽包压力和汽轮机一级压力构造锅炉总热量信号,采用一高阶惯性环节描述锅炉给煤量与总热量信号之间的动态关系。针对此高阶惯性环节设计状态观测器,得到具有较小惯性的中间状态信号,利用此信号除以给煤量得到煤发热量信号。实验结果表明:动态补偿后的煤发热量具有抗扰动性好、响应速度快的优点。
关键词 :煤发热量,软测量,状态观测器,动态补偿
近年软测量在火电厂热工参数测量中得到广泛应用,如烟气流量、飞灰含碳量软测量等等。软测量定义为:通过一定的数学方法,寻找与待测量变量相关的辅助变量,建立辅助变量与待测量变量之间的数学模型,进而得到对待测变量的最佳估计。一些通过机理分析结合数据拟合得到的软测量结果在稳态工况下使用没问题,但在动态工况下往往存在巨大的动态误差。目前大型机组都参与电网调峰及一次调频,机组负荷变化频繁,有必要研究软测量动态补偿方法,以提高测量的快速性、减小动态误差。杜文莉等人针对静态软测量无法反映工业过程的动态信息这一问题,提出一种基于最小二乘支持向量机和自回归–滑动平均模型的软测量建模方法,对静态模型增加动态校正环节,实现了对静态模型的动态校正,以改善系统动态响应特性 [
热工对象往往具有大惯性、大迟延,利用状态观测器将对象中惯性小、响应速度快的中间状态进行重构,可以实现状态反馈控制,提高系统响应速度 [
正常工况下,煤发热量Qar (MJ/kg)计算方法 [
式中:Qb为锅炉总热量(MW);uB为机组总给煤量(kg/s)。
典型机组一般采用皮带称重式给煤机,uB由各个给煤机给煤量求和得到。根据能量平衡原理,采用以下公式得到Qb:
式中:Q1为锅炉汽水系统有效吸热量(MW);Q2为锅炉排烟总焓(MW);q3为可燃气体未完全燃烧热损失(%);q4为固体未完全燃烧热损失(%);q5为锅炉散热损失(%);q6为灰渣物理热损失(%)。
文献 [
公式(1)在静态工况下使用是没有问题的,但在机组变负荷时,直接用公式(1)会出现巨大的动态误差。例如升负荷时,煤发热量没变,给煤量首先增加,经过给煤机、磨煤机、一次风粉管道后进入炉膛燃烧释放热量,经汽水系统吸收表现为蒸汽流量、压力变化,锅炉总热量要经过一段时间才达到稳定,导致给煤量与燃料总热量在时间上的不一致,计算出来的煤发热量在这段时间降低。另外当煤质真正变化时,由于制粉惯性和蓄热惯性,计算的煤发热量响应时间较长,不能及时反映煤质变化。针对这些问题,研究一种基于状态观测器的动态补偿方法,以提高其变负荷时的抗扰动性及变煤质时的动态响应速度。
以二阶对象为例,设被控对象传递函数如下:
式中,U(s)为给煤量的拉氏变换(kg/s);Y(s)为总热量的拉氏变换(MJ);Tb为对象惯性时间(s);K为对象增益(MJ/kg),在实际过程中可能发生变化。
对象中存在两个状态变量,离输入最近的状态响应必然快于输出响应。给对象补偿一个惯性环节
根据式(4)设计状态观测器,将重构小惯性状态
以Tb = 100,Ts = 1,K = 16为例进行两种情况下的仿真:t = 250 s时输入变而增益不变如图2曲线a所示;在t = 1250 s时增益变而输入不变如曲线b所示。Y(s)如曲线c所示,
状态观测器误差校正环节的参数确定:首先采用极点配置方法确定Kp、g2、g3,得到最靠近输入且惯性小的重构状态
图1. 基于状态观测器的增益软测量动态补偿模型
图2. 增益软测量动态补偿仿真曲线
设计低通滤波器对动态补偿后的K进行滤波,要求滤波器能够消除输入量变化时造成的扰动。
研究选用的是中值限速滤波方法——先进行限速滤波再进行中值滤波。中值滤波和限速滤波均具有低通滤波特性,其中,中值滤波具有完全消除尖峰脉冲(脉冲长度小于滤波器宽度)和完全通过阶跃和斜坡的特点;限速滤波器由非线性速率限制环节构成,具有低通滤波性随输入幅值的变大而变强和快速消除尖峰脉冲的特点。
中值限速滤波器的滤波长度N = 3,一个信号为原始信号,另外两个为原始信号经过限速滤波后的信号,依靠经验选择速率限制分别为±0.25、±0.125。图2曲线g为滤波后的增益,Klv表示。可以看出,无论是输入变化还是增益变化,动态补偿后的增益经过滤波,抗输入扰动性和动态响应速度均有显著改善。
根据现有的机组负荷–压力简化非线性模型给出了Matlab中simulink模型图,如图3所示。图中,uB为给煤量(kg/s);τ为制粉过程的迟延时间(s);Kf为燃料动态时间(s);rB为锅炉燃烧率(kg/s);K1为燃料增益(MJ/kg);Qb为锅炉总热量(MW);μT为汽轮机调门开度(%);K2为过热器阻力系数;Cb为锅炉蓄热系数;Pd为汽包压力(MPa);Pt为机前压力(MPa);P1为汽轮机一级压力(MPa);K3为汽轮机调门指令增益;P0为蒸汽做功(MW);Kt为汽轮机动态时间(s);Ne为机组负荷(MW)。
文献 [
式中,Qar为煤发热量(MJ/kg),ηu为机组发电效率,可由机组设计的负荷–发电煤耗曲线确定,这样煤发热量仅与K1相关,求得K1即可知煤发热量。
由图3可知,机组负荷–压力简化的非线性模型主要分为三个动态部分 [
由图3得锅炉总热量信号构造公式:
式(6)得到的Qb静态精度高但动态性差,文献 [
式中:qa为锅炉总风量(Nm3/s);O2为锅炉排烟氧量(%);Kaq为煤空气热量比(Nm3/MJ),对于大部分煤种近似为0.263;T0为滤波器的惯性时间(s),决定了滤波截止频率和带宽,目前主要靠经验选取。
联立式(6)~(7)得到了更接近实际的锅炉总热量,煤的制粉及炉内燃烧过程就简化为单入–单出环节。
图3. 机组负荷–压力非线性简化模型
采用上节思路,对此环节设计状态观测器。
对于带有纯迟延环节的系统设计状态观测器时,往往是先用高阶惯性代替纯迟延。人们建立图3的模型时为了避免繁琐的高阶结构,进行了适当的简化,其中制粉动态特性简化为纯迟延和一阶惯性。实际上,用高阶惯性描述制粉系统的动态过程更能模拟实际动态过程。
结合图3,给出制粉及炉内燃烧释放热量的高阶惯性描述:
式中,T1表示煤在给煤机和一次风粉管道的动态时间(s);T2为煤粉在炉膛内停留及燃烧的动态时间(s);
据此设计的状态观测器模型及燃料增益软测量动态补偿模型如图4所示。图中
现以660 MW机组的负荷–压力非线性模型为例,进行基于状态观测器的煤发热量动态补偿方法验证。模型中各个参数为τ = 23,Kf = 189,K1 = 19.1,K2 = 0.000122,K3 = 0.537,Kt = 12,Cb = 7690。图3中制粉及炉内燃烧部分(虚线框)的锅炉总热量是可测的,那么虚线框部分就从模型中分离出来成为一个单入单出系统。利用响应曲线的方法,求得制粉系统高阶惯性描述中的参数为:T1 = 8.13,T2 = 1,
图5中
图4. 基于状态观测器的煤发热量软测量动态补偿模型
在某600 MW机组为例,利用其运行过程中的历史数据,计算出实际的煤发热量总热量。将其作为状态观测器模型中的输入,在DCS中以组态方式实现补偿逻辑,观察各种工况下经过补偿和未经补偿的煤发热量软测量结果。
图6为机组变负荷工况下软测量结果对比,可以看出,动态补偿后得到的煤发热量经过负荷修正,基本不受负荷变化的影响,具有较好的抗负荷扰动性能。而未经补偿的煤发热量计算值,升负荷时偏低,降负荷时偏高。
在现场进行煤质变化实验非常困难,这里采用等效模拟实验方法进行验证,其思路是:保持总给煤量不变,人为减少锅炉总热量,模拟煤发热量突然降低的情况。减少锅炉总热量的做法:以式(8)所示的
图5. 基于状态观测器动态补偿的煤发热量仿真曲线
图6. 变负荷工况下计算煤发热量
图7. 模拟煤质变化工况下计算煤发热量
图8. 启/停磨工况下计算煤发热量
制粉及炉内燃烧释放热量动态过程为模拟对象,根据发电负荷曲线选择一段稳定负荷运行的工况,在t = 500 s时改变对象中的燃料增益,将式(8)的输出加入噪声修正后作为模拟所需的锅炉总热量。由图7可以看出,与补偿前相比,补偿后的煤发热量能够快速、准确的反映煤质变化。
图8为机组负荷由380 MW下降至350 MW,期间停一台磨煤机,对应停一层燃烧器情况下计算得到的煤发热量。启/停磨会对锅炉燃烧状态造成强烈扰动,可发现动态补偿后的煤发热量基本不受磨启/停影响,具有很好的抗燃烧扰动能力。
本文进行了基于状态观测器的热工参数软测量动态补偿方法研究,并将其应用于煤发热量软测量中,通过实验验证了该方法的可行、有效性。并得到了以下结论:
1) 利用机组负荷–压力简化模型,由汽包压力和一级压力构造锅炉热量信号,并将制粉及炉内燃烧动态过程简化为一个单入单出高阶惯性环节,便于状态观测器设计的同时也更加符合实际情况。
2) 在设计状态观测器时引入了PID结构,实现了对象增益变化后,重构状态仍可以快速、无静差的逼近真实状态。利用观测到的小惯性状态对软测量值进行动态补偿,有效提高了软测量的速度,减小了动态误差。
3) 实验表明,机组负荷变化频繁时,不加动态补偿得到的软测量结果远不能满足机组要求。采用本文动态补偿方法得到的煤发热量结果具有很好的抗负荷扰动和抗燃烧扰动能力,在煤质变化时能够快速、准确的计算出新的煤发热量值,及时反映煤质变化。
河北省教育厅科研计划重点项目基金资助(ZD2016053),河北师范大学自然科学科研基金青年基金(L2015Q10)。
武现聪,田亮,杜文霞. 基于状态观测器的软测量动态补偿方法Dynamic Compensation Method for Soft Sensing Based on State Observer[J]. 建模与仿真, 2016, 05(03): 89-97. http://dx.doi.org/10.12677/MOS.2016.53012