一种具有系统非线性特性的汽包锅炉动态模型.docxVIP

一种具有系统非线性特性的汽包锅炉动态模型 1 单元机组模型的建立 单元装置的协调管理系统将炉设备设计作为一个整体进行控制，控制过程非常复杂。影响机组动态特性的因素除了其内在的物理结构属性以外,还与机组的运行方式、运行工况、外部环境等因素有关。一般说来,单元机组的动态特性是本质非线性的,并具有分布参数和时变特性。对单元机组动态特性的精确描述目前还难以得到。只有通过合理的简化与近似,采用机理分析或系统辨识方法才能够建立起来满足一定精度要求的动态特性数学模型。 常用的协调控制系统动态模型一般采用对非线性模型的某一特定工作点进行线性化而得到。但是这种线性化模型无法体现系统的非线性特性,因此利用这种工作点线性化的模型来验证协调控制系统的有效性是片面的。近年来较多的研究集中在构造一种通用的非线性协调控制系统动态仿真模型上。 2 锅炉动态特性 单元机组协调控制系统的被控对象是一种存在强烈耦合特征的、复杂的多变量系统。图1表示了协调控制系统动态模型的结构。它实际上是一种低阶的非线性模型。系统的结构一方面反映了能量平衡关系,另一方面反映了系统中存在的本质非线性特征。 系统的动态模型结构体现了2种相关的能量平衡关系:①汽包压力Pb反映了锅炉吸热量Dq和汽包出口蒸汽发热量Dk的平衡。汽包蓄热系数Cb反映了汽包蓄能的大小。② 主蒸汽压力Pt反映了汽包出口蒸汽发热量Dk与主蒸汽发热量Dt的平衡。主蒸汽管道蓄热系数Ct反映了主蒸汽管道蓄能的大小。 假设炉内燃烧与传热过程可简化为具有纯时延的一阶惯性特性,其传递函数可表示为 ΔDq=1Τ2s+1?k1e-τ1sΤ1s+1ΔB(1) 式中B为进入炉膛的燃料量;Dq为锅炉总有效吸热量。 考虑到锅炉汽包的蓄能,可得到第1个能量平衡方程: ΔDq-ΔDk=CddΔΡbdt(2) 式中Dk为汽包出口蒸汽对应的发热量。 考虑到主蒸汽管道的蓄能,可得到第2个能量平衡方程: ΔDk-ΔDt=CtdΔΡtdt(3) 式中Dt为主蒸汽流量对应的发热量。 锅炉过程动态中存在的非线性特征主要反映在2个方面: (1) 汽包压力Pb和主蒸汽压力Pt的压力降同汽包出口蒸汽流量Dk之间存在着平方根关系; (2) 主蒸汽流量Dt同汽轮机调节阀通流面积μ和主蒸汽压力Pt的乘积成比例关系。 可用下面2个方程描述: Pb-Pt=kD2k(4) Dt=μPt(5) 3 阀门开度扰动试验 图2显示了协调控制系统的动态仿真模型。限于篇幅,图中省略了相对简单的磨煤机动态过程。为了建立上述动态仿真模型,需要确定几个未定参数K1、K2等。K1为汽包的蓄热系数(量纲为s),K3为主蒸汽管道蓄热系数(量纲为s)。K4、K5、K6为内部系数,用于仿真模型的内部平衡。 在确定这些系数之前,需要获得机组的额定参数。假设机组汽包压力额定参数为Pbe,主蒸汽压力额定参数为Pte,主蒸汽流量额定参数为Dte。汽包压力的额定参数设置为积分器1的初始值,主蒸汽压力的额定参数设置为积分器2的初始值。 假设锅炉总有效吸热量Dq和汽轮机阀门开度指令均已标么化,则 Κ6=DteΡte(6) 那么由式(4)可得 Κ2√Ρbe-Ρte=Dke(7) 对上式进行标么化处理 Κ2=1.0√(Ρbe-Ρte)/Ρte(8) 为了使得加法器1的入口平衡,则必须满足 Κ4?B=√Ρbe-ΡteΚ2(9) 为了使得加法器3的入口平衡,必须满足 Ρte?μΚ5=√Ρbe-ΡteΚ2(10) 为得到K1和K3,需要作汽轮机阀门开度扰动试验。保持燃烧率不变,作阀门开度扰动试验,虽然主蒸汽流量和汽包压力都在变化,但是热量信号不会改变。汽包蓄热系数可由下列方程得到: ?Cb=t1∫0[Dt(t)-Dt(0)]dtΡb(0)-Ρb(t1)t/h?sΜΡa(11) 上式求出的蓄热系数是具有工程量纲意义上的数。要将上式求出的蓄热系数转换为数学模型中的K1和K3,一般地 Cb=?Cb?Κ5/Κ6(12) 式中Cb为蓄热时间常数,s。 上述方法求出的是整个锅炉的蓄热系数。而协调系统模型是需要将汽包的蓄热同主蒸汽管道的蓄热分开。理论分析和工程实践已经证明,汽包的蓄热系数通常占90%,而主蒸汽管道的蓄热系数通常占10%。那么 Cb=Κ1+Κ2Κ1≈90%?Cb?Κ3≈10%?Cb(13) 式(6)至式(13)给出了图2中未定参数的求取方法。K1代表了锅炉汽包的蓄热系数,K2反映了主蒸汽管道的沿程阻力,K3代表了主蒸汽管道的蓄热系数,K4为锅炉总有效吸热量的标定系数,K5为主蒸汽流量与对应发热量之间的标定系数,K6为主蒸汽流量标定系数。 4 负荷压力的渐进式预测 4.1 主蒸汽流量模型 从多变量控制系统的角度考虑,在机组的动态过程中,通常可以将主蒸汽压力和主蒸汽流量的响应认为是锅炉和汽轮机共同作用的结果。 假设负荷和压力的增