CN104102134A - 一种通过性能指标实现再热汽温多变量协调预测控制的方法 - Google Patents

一种通过性能指标实现再热汽温多变量协调预测控制的方法 Download PDF

Info

Publication number
CN104102134A
CN104102134A CN201410351300.4A CN201410351300A CN104102134A CN 104102134 A CN104102134 A CN 104102134A CN 201410351300 A CN201410351300 A CN 201410351300A CN 104102134 A CN104102134 A CN 104102134A
Authority
CN
China
Prior art keywords
reheat steam
steam temperature
control
moment
spray valve
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410351300.4A
Other languages
English (en)
Other versions
CN104102134B (zh
Inventor
雎刚
钱晓颖
陈欢乐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201410351300.4A priority Critical patent/CN104102134B/zh
Publication of CN104102134A publication Critical patent/CN104102134A/zh
Application granted granted Critical
Publication of CN104102134B publication Critical patent/CN104102134B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Feedback Control In General (AREA)

Abstract

本发明公开了一种通过性能指标实现再热汽温多变量协调预测控制的方法,该方法将火电机组再热汽温系统看出一个2×1的多变量系统,基于新型性能指标,设计多变量预测控制器,同时控制再热汽温减温喷水阀门开度和烟气挡板开度。在该性能指标中,一方面基于对单步预测偏差及预测偏差变化率的要求,简化算法,实现快速控制,还能通过调整预测长度和预测偏差变化率权值有效调整控制系统的控制性能,另外一方面在指标中加入了减温喷水阀门开度的动态限制,该部分性能参数在控制过程中自适应变化,在提升再热汽温控制质量的同时,可减少喷入再热蒸汽的减温水量,有效提高机组的循环效率。

Description

一种通过性能指标实现再热汽温多变量协调预测控制的方法
技术领域
本发明属于热工自动控制技术领域,具体涉及一种通过性能指标实现再热汽温多变量协调预测控制的方法。
背景技术
再热汽温是火电厂热力系统中的重要参数,再热汽温控制品质的优劣直接影响到整个机组的安全和经济运行,再热汽温控制系统是机组的重要控制系统之一。然而无论采用蒸汽侧控制还是烟气侧控制,再热汽温对象都具有较大的延迟和惯性,再热汽温较难控制。目前,考虑到烟气挡板调节的可控性和精度都较差,一般采用两个独立的PID控制器分别实现蒸汽侧和烟气侧的调节,其中蒸汽侧的调节作为辅助调节。这种控制方法,一方面针对大惯性的再热汽温对象控制效果较差,同时蒸汽侧和烟气侧无法实现协调控制,蒸汽侧的减温喷水会影响机组的经济性。
也有一些先进控制策略的再热汽温控制方案,比如单变量预测控制、普通多变量预测控制等。但是这些方法一方面算法复杂,计算量大,无法实现工程在线应用,另外一方面控制算法中缺乏有效的控制参数对系统控制性能进行调节,不利于工程应用。同时,这些方法没有考虑机组的经济性,蒸汽侧减温喷水量无法优化,不能够实现蒸汽侧和烟气侧的协调控制。
发明内容
发明目的:针对上述问题,提出一种通过性能指标实现再热汽温多变量协调预测控制的方法,可以解决再热汽温调节动静态偏差大,蒸汽侧喷水量无法优化,影响机组经济性,以及普通预测控制计算量大,无法实现在线控制,缺乏有效控制参数对控制系统性能进行调节等问题。
技术方案:本发明提出一种新型性能指标,据此设计多变量预测控制器,并将其应用于再热汽温的控制系统中。在该指标中以烟气侧调节为主,协调控制烟气侧挡板开度和蒸汽侧减温喷水阀门开度量;指标中主要考虑再热汽温的单步预测偏差和预测偏差变化率的关系,同时在指标中添加蒸汽侧减温喷水阀门开度量的自适应限制,在控制过程中动态优化蒸汽侧减温喷水量,在保证再热汽温控制品质的同时,实现机组的经济运行。
本发明的实现步骤如下:
步骤1:获取再热汽温对象的阶跃响应模型,在稳态的工况下,分别作减温喷水阀门开度和烟气挡板开度的开环阶跃响应试验,获取系统响应数据,使用阶跃响应模型辨识方法,辨识得到减温喷水阀门开度—再热汽温和烟气挡板开度—再热汽温的传递函数模型;
步骤2:选择采样周期T,T首先满足香农定理,然后使T95/T取值在10到25之间,T95为再热汽温过渡过程上升到95%的调节时间,通过对步骤1中的传递函数模型进行离散化,得到再热汽温对象的CARIMA模型为: A ( z - 1 ) y ( k ) = B 1 ( z - 1 ) B 2 ( z - 1 ) · u 1 ( k - 1 ) u 2 ( k - 1 ) + 1 Δ ϵ ( k ) , 式中 A ( z - 1 ) = 1 + Σ i = 1 n a a i z - i , B 1 ( z - 1 ) = Σ i = 0 n b 1 b 1 , i z - i , B 2 ( z - 1 ) = Σ i = 0 n b 2 b 2 , i z - i , z-1为后移算子,Δ=1-z-1,y(k)为k时刻的再热汽温测量值,u1(k-1)为k-1时刻减温喷水阀门开度控制量,u2(k-1)为k-1时刻烟气挡板开度控制量,ε(k)为白噪声,ai、b1,i、b2,i分别为多项式A(z-1),B1(z-1)和B2(z-1)中z-i的系数,na、nb1、nb2分别为多项式A(z-1),B1(z-1)和B2(z-1)的阶次;
步骤3:根据丟番图方程1=Ej(z-1)A(z-1)Δ+z-jFj(z-1)解得控制变量多项式: E j ( z - 1 ) = Σ i = 0 j - 1 e j , i z - i , F j ( z - 1 ) = Σ i = 0 n a f j , i z - i , 其中ej,i、fj,i分别为多项式Ej(z-1),Fj(z-1)中z-i的系数,j=N-1,N,N为预测时域长度,其选择需大于再热汽温系统的纯延迟时间;
步骤4:定义多项式 G j 1 ( z - 1 ) = E j ( z - 1 ) B 1 ( z - 1 ) , G j 2 ( z - 1 ) = E j ( z - 1 ) B 2 ( z - 1 ) , 求得: G j 1 ( z - 1 ) = Σ i = 0 j + n b 1 - 1 g j , i 1 z - i , G j 2 ( z - 1 ) = Σ i = 0 j + n b 2 - 1 g j , i 2 z - i , 式中,分别为多项式中z-i项的系数;定义 f i = Σ i = 0 n a f j , i y ( k - i ) + Σ i = 1 n b 1 g j , j - 1 + i 1 Δ u 1 ( k - i ) + Σ i = 1 n b 2 g j , j - 1 + i 2 Δ u 2 ( k - i ) , 式中,Δu1(k-i),i=1…nb1,Δu2(k-i),i=1…nb2分别表示在k-i时刻减温喷水阀门开度和烟气挡板开度的输出控制量的增量,j=N-1,N;
步骤5:按下式对再热汽温系统未来第k+N-1和k+N时刻的输出进行预测:y(k+j|k)=Gj·ΔU+fj(j=N-1.N),式中,y(k+j|k)为未来k+j时刻再热汽温预测值,Gj为一维行向量,其中 G N - 1 = [ g N - 1 , N - 2 1 , . . . , g N - 1 , N - 1 - Nu 1 1 , g N - 1 , N - 2 2 , . . . , g N - 1 , N - 1 - Nu 2 2 ] , G N = [ g N , N - 1 1 , . . . g N , N - Nu 1 1 , g N , N - 1 2 , . . . , g N , N - Nu 2 2 ] , Nu1、Nu2分别为减温喷水阀门开度和烟气挡板开度的控制时域长度;ΔU=[Δu1(k),…,Δu1(k+Nu1-1),Δu2(k),…,Δu2(k+Nu2-1)],其中Δu1(k+i),i=0…Nu1-1,Δu2(k+i),i=0…Nu2-1分别表示在k时刻对未来k+i时刻减温喷水阀门开度和烟气挡板开度输出控制量增量的估计值;
步骤6:根据新型性能指标: J = [ e ( k + N ) + λ 1 ce ( k + N ) ] 2 + λ 2 Σ i = 0 Nu 1 - 1 u 1 ( k + i ) , 计算减温喷水阀门开度控制量增量Δu1(k)和烟气挡板开度控制量增量Δu2(k);
式中,e(k+N)=y(k+N|k)-yr(k)为第k+N时刻的预测偏差,式中yr(k)为再热汽温在k时刻的设定值;对应e(k+N-1)为k+N-1时刻的预测偏差;ce(k+N)=e(k+N)-e(k+N-1)为k+N时刻的预测偏差变化率;u1(k+i)为k+i时刻减温喷水阀门开度控制量的估计值;
λ1和λ2为控制参数,λ1一般取值范围为[0,5],λ2按下式在过程控制中自适应变化:k1为系数,取值范围为[0,1],将步骤5中的未来预测值带入上述性能指标中,并通过极值必要条件求得:ΔU=-(PT·P+QT·RT·R·Q)-1·[h·P+u1(k-1)·LT·RT·R·Q]T,式中,P=GN1·GN1·GN-1,h=fN-yr(k)+λ1·(fN-fN-1), Q2为Nu1×Nu2零矩阵,Q=[Q1,Q2], R = diag [ λ 2 , . . . , λ 2 ] Nu 1 ;
步骤7:计算减温喷水阀门开度控制量u1(k)=u1(k-1)+Δu1(k),和烟气挡板开度控制量u2(k)=u2(k-1)+Δu2(k),式中,u1(k-1)和u2(k-1)分别为减温喷水阀门开度和烟气挡板开度在k-1时刻的对应控制量;
步骤8:根据执行机构限制及实际运行情况,设定控制量约束,减温喷水阀门开度约束上下限为[u1,min,u1,max],烟气挡板开度约束上下限为[u2,min,u2,max],如果ui(k)>ui,max,则令ui(k)=ui,max,Δui(k)=ui,max-ui(k-1);如果ui(k)<ui,min,则令ui(k)=ui,min,Δui(k)=ui,min-ui(k-1),i=1,2;
步骤9:控制系统输出u1(k)和u2(k)到相应执行机构,在之后的每个采样周期内重复执行步骤5到步骤9。
有益效果:基于新型性能指标,设计再热汽温多变量预测控制算法,该算法能够实现再热汽温的快速稳定控制,在控制过程中,能保证减温喷水量和烟气挡板开度的动态协调,系统稳态时减温水量为零,保证机组运行经济性。同时,该算法中的预测时域长度N选取上只需大于纯延迟时间,控制时域长度一般取1即可,算法单步预测计算,计算量很小,有利于在线控制;可以通过调节预测时域长度N来有效调节控制系统的快速性,通过调节控制参数λ1有效调节控制系统的稳定性。
附图说明
图1是通过性能指标实现再热汽温多变量协调预测控制方法的控制结构图;
图2(a)-2(c)是通过性能指标实现再热汽温多变量协调预测控制系统在设定值阶跃扰动下的响应曲线;
图3(a)-(3c)是通过性能指标实现再热汽温多变量协调预测控制系统在外部阶跃扰动下的响应曲线。
具体实施方式
下面结合附图和具体实施例,对本发明做更进一步的解释。
步骤1:获取再热汽温对象的阶跃响应模型,在稳态的工况下,分别作减温喷水阀门开度和烟气挡板开度的开环阶跃响应试验,获取系统响应数据,使用阶跃响应模型辨识方法,辨识得到减温喷水阀门开度—再热汽温和烟气挡板开度—再热汽温的传递函数模型;
本实施例中再热器出口汽温y(℃)对减温喷水阀门开度u1(%)响应的传递函数为:
G 1 ( s ) = - 0.5889 ( 2 s + 1 ) ( 30 s + 1 ) 2 ( 145 s + 1 ) 2 e - 180 s
再热器出口汽温y(℃)对烟气挡板开度u2(%)响应的传递函数为:
G 2 ( s ) = 0 . 8155 ( 50 s + 1 ) 2 ( 110 s + 1 ) 2 e - 140 s
步骤2:选择采样周期T=10s。通过对步骤1中的传递函数模型进行离散化,得到再热汽温对象的CARIMA模型为: A ( z - 1 ) y ( k ) = B 1 ( z - 1 ) B 2 ( z - 1 ) · u 1 ( k - 1 ) u 2 ( k - 1 ) + 1 Δ ϵ ( k ) , 式中:
A(z-1)=1-6.77z-1+20.03z-2-33.82z-3+35.67z-4-24.06z-5+10.15z-6-2.46z-7+0.27z-8-0.0017z-9
B1(z-1)=z-14(1×10-5+6.5×10-5z-1+1.96×10-4z-2+9.63×10-5z-3+1.18×10-4z-4-1.21×10-4z-5+2.43×10-5z-6+3×10-6z-7-2.1310-8z-8)
B2(z-1)=z-18(-5.91×10-6-6.21×10-5z-1+1.55×10-4z-2-1.15×10-5z-3-2.02×10-4z-4+1.51×10-4z-5-1.71×10-5z-6-8.63×10-6z-7-7.79×10-8z-8)
式中,z-1为后移算子,Δ=1-z-1;ε(k)为白噪声;
步骤3:预测时域长度N选取为40,根据丟番图方程1=Ej(z-1)A(z-1)Δ+z-jFj(z-1)解得控制变量多项式: 式中,ej,i、fj,i分别为多项式Ej(z-1),Fj(z-1)中z-i的系数,j=N-1,N;
本实施例中,FN-1(z-1)和FN(z-1)用向量形式表示为:
FN-1(z-1)=[2160267.2,-14462687.3,42333166.4,-70765975.2,73907169.9,-49404550.7,20666166.9,-4963137.2,532952.8,-3371.9];
FN(z-1)=[2322980.6,-15560238.2,45568406.9,-76209835.5,79627836,-53250803.3,22283817.2,-5353602.2,575077.9,-3638.5];
步骤4:定义多项式 G j 1 ( z - 1 ) = E j ( z - 1 ) B 1 ( z - 1 ) , G j 2 ( z - 1 ) = E j ( z - 1 ) B 2 ( z - 1 ) , 求得: G j 1 ( z - 1 ) = Σ i = 0 j + n b 1 - 1 g j , i 1 z - i , G j 2 ( z - 1 ) = Σ i = 0 j + n b 2 - 1 g j , i 2 z - i , 式中分别为多项式中z-i项的系数;定义 f i = Σ i = 0 n a f j , i y ( k - i ) + Σ i = 1 n b 1 g j , j - 1 + i 1 Δ u 1 ( k - i ) + Σ i = 1 n b 2 g j , j - 1 + i 2 Δ u 2 ( k - i ) , 式中Δu1(k-i),i=1…nb1,Δu2(k-i),i=1…nb2分别表示在k-i时刻减温喷水阀门开度和烟气挡板开度的输出控制量的增量,j=N-1,N;
步骤5:减温喷水阀门和烟气挡板开度的控制时域长度Nu1、Nu2都取为1,按下式对再热汽温系统未来第k+N-1和k+N时刻的输出进行预测:y(k+j|k)=Gj·ΔU+fj(j=N-1.N),式中y(k+j|k)为未来k+j时刻再热汽温预测值,Gj为一维行向量,
本实施例中,GN-1=[0.3251,-0.1619],GN=[0.3471,-0.176];
步骤6:根据新型性能指标: J = [ e ( k + N ) + λ 1 ce ( k + N ) ] 2 + λ 2 Σ i = 0 Nu 1 - 1 u 1 ( k + i ) , 计算减温喷水阀门开度控制量增量Δu1(k)和烟气挡板开度控制量增量Δu2(k);
式中,e(k+N)=y(k+N|k)-yr(k),为第k+N时刻的预测偏差,其中yr(k)为再热汽温在k时刻的设定值;对应e(k+N-1)为k+N-1时刻的预测偏差;ce(k+N)=e(k+N)-e(k+N-1)为k+N时刻的预测偏差变化率;u1(k+i)为k+i时刻减温喷水阀门开度控制量的估计值;
λ1和λ2为控制参数,本实施例中λ1取1,λ2按下式在过程控制中自适应变化:k1为系数,本实施例中取0.1;
将步骤5中的未来预测值带入上述性能指标中,并通过极值必要条件 dJ dΔU = 0 求得:
ΔU=-(PT·P+QT·RT·R·Q)-1·[h·P+up(k-1)·LT·RT·R·Q]T
式中,P=GN1·GN1·GN-1,h=fN-yr(k)+λ1·(fN-fN-1), Q2为Nu1×Nu2零矩阵,Q=[Q1,Q2], R = diag [ λ 2 , . . . , λ 2 ] Nu 1 ;
步骤7:计算减温喷水阀门开度控制量u1(k)=u1(k-1)+Δu1(k),和烟气挡板开度控制量u2(k)=u2(k-1)+Δu2(k),式中,u1(k-1)和u2(k-1)分别为减温喷水阀门开度和烟气挡板开度在k-1时刻的对应控制量;
步骤8:根据执行机构限制及实际运行情况,设定控制量约束,减温喷水阀门开度约束上下限为[u1,min,u1,max],烟气挡板开度约束上下限为[u2,min,u2,max]。如果ui(k)>ui,max,则令ui(k)=ui,max,Δui(k)=ui,max-ui(k-1);如果ui(k)<ui,min,则令ui(k)=ui,min,Δui(k)=ui,min-ui(k-1),i=1,2;
本实施例中,减温喷水阀门开度约束上下限为[0,100];烟气挡板开度约束上下限为[0,100];
步骤9:控制系统输出u1(k)和u2(k)到相应执行机构,在之后的每个采样周期内重复执行步骤5到步骤9。
上述实施例的仿真结果如图2(a)至图3(c)所示,图2(a)至图2(c)为再热汽温控制系统在设定值阶跃扰动下的控制情况;图3(a)至图3(c)为再热汽温控制系统在外部阶跃扰动下的控制情况。图2(a)和图3(a)表明当再热汽温系统发生扰动时,本发明算法的控制系统能够快速调节再热汽温跟踪设定值,具有良好的控制性能。同时控制过程中,保证减温喷水量和烟气挡板开度的动态协调,系统稳态时减温喷水量为零,保证机组运行经济性。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (3)

1.一种通过性能指标实现再热汽温多变量协调预测控制的方法,其包括以下步骤:
步骤1:获取再热汽温对象的阶跃响应模型,在稳态的工况下,分别作减温喷水阀门开度和烟气挡板开度的开环阶跃响应试验,获取系统响应数据,使用阶跃响应模型辨识方法,辨识得到减温喷水阀门开度—再热汽温和烟气挡板开度—再热汽温的传递函数模型;
步骤2:选择采样周期T,T首先满足香农定理,然后使T95/T取值在10到25之间,T95为再热汽温过渡过程上升到95%的调节时间,通过对步骤1中的传递函数模型进行离散化,得到再热汽温对象的CARIMA模型为: A ( z - 1 ) y ( k ) = B 1 ( z - 1 ) B 2 ( z - 1 ) · u 1 ( k - 1 ) u 2 ( k - 1 ) + 1 Δ ϵ ( k ) , 式中 A ( z - 1 ) = 1 + Σ i = 1 n a a i z - i , B 1 ( z - 1 ) = Σ i = 0 n b 1 b 1 , i z - i , B 2 ( z - 1 ) = Σ i = 0 n b 2 b 2 , i z - i , z-1为后移算子,Δ=1-z-1,y(k)为k时刻的再热汽温测量值,u1(k-1)为k-1时刻减温喷水阀门开度控制量,u2(k-1)为k-1时刻烟气挡板开度控制量,ε(k)为白噪声,ai、b1,i、b2,i分别为多项式A(z-1),B1(z-1)和B2(z-1)中z-i的系数,na、nb1、nb2分别为多项式A(z-1),B1(z-1)和B2(z-1)的阶次;
步骤3:根据丟番图方程1=Ej(z-1)A(z-1)Δ+z-jFj(z-1)解得控制变量多项式: E j ( z - 1 ) = Σ i = 0 j - 1 e j , i z - i , F j ( z - 1 ) = Σ i = 0 n a f j , i z - i , 其中ej,i、fj,i分别为多项式Ej(z-1),Fj(z-1)中z-i的系数,j=N-1,N,N为预测时域长度,其选择需大于再热汽温系统的纯延迟时间;
步骤4:定义多项式 G j 1 ( z - 1 ) = E j ( z - 1 ) B 1 ( z - 1 ) , G j 2 ( z - 1 ) = E j ( z - 1 ) B 2 ( z - 1 ) , 求得: G j 1 ( z - 1 ) = Σ i = 0 j + n b 1 - 1 g j , i 1 z - i , G j 2 ( z - 1 ) = Σ i = 0 j + n b 2 - 1 g j , i 2 z - i , 式中,分别为多项式中z-i项的系数;定义 f i = Σ i = 0 n a f j , i y ( k - i ) + Σ i = 1 n b 1 g j , j - 1 + i 1 Δ u 1 ( k - i ) + Σ i = 1 n b 2 g j , j - 1 + i 2 Δ u 2 ( k - i ) , 式中,Δu1(k-i),i=1…nb1,Δu2(k-i),i=1…nb2分别表示在k-i时刻减温喷水阀门开度和烟气挡板开度的输出控制量的增量,j=N-1,N;
步骤5:按下式对再热汽温系统未来第k+N-1和k+N时刻的输出进行预测:y(k+j|k)=Gj·ΔU+fj(j=N-1.N),式中,y(k+j|k)为未来k+j时刻再热汽温预测值,Gj为一维行向量,其中 G N - 1 = [ g N - 1 , N - 2 1 , . . . , g N - 1 , N - 1 - Nu 1 1 , g N - 1 , N - 2 2 , . . . , g N - 1 , N - 1 - Nu 2 2 ] , G N = [ g N , N - 1 1 , . . . g N , N - Nu 1 1 , g N , N - 1 2 , . . . , g N , N - Nu 2 2 ] , Nu1、Nu2分别为减温喷水阀门开度和烟气挡板开度的控制时域长度;ΔU=[Δu1(k),…,Δu1(k+Nu1-1),Δu2(k),…,Δu2(k+Nu2-1)],其中Δu1(k+i),i=0…Nu1-1,Δu2(k+i),i=0…Nu2-1分别表示在k时刻对未来k+i时刻减温喷水阀门开度和烟气挡板开度输出控制量增量的估计值;
步骤6:根据新型性能指标: J = [ e ( k + N ) + λ 1 ce ( k + N ) ] 2 + λ 2 Σ i = 0 Nu 1 - 1 u 1 ( k + i ) , 计算减温喷水阀门开度控制量增量Δu1(k)和烟气挡板开度控制量增量Δu2(k);
式中,e(k+N)=y(k+N|k)-yr(k)为第k+N时刻的预测偏差,式中yr(k)为再热汽温在k时刻的设定值;对应e(k+N-1)为k+N-1时刻的预测偏差;ce(k+N)=e(k+N)-e(k+N-1)为k+N时刻的预测偏差变化率;u1(k+i)为k+i时刻减温喷水阀门开度控制量的估计值;
λ1和λ2为控制参数,λ1一般取值范围为[0,5],λ2按下式在过程控制中自适应变化:k1为系数,取值范围为[0,1],将步骤5中的未来预测值带入上述性能指标中,并通过极值必要条件求得:
ΔU=-(PT·P+QT·RT·R·Q)-1·[h·P+u1(k-1)·LT·RT·R·Q]T
式中,P=GN1·GN1·GN-1,h=fN-yr(k)+λ1·(fN-fN-1), Q2为Nu1×Nu2零矩阵,Q=[Q1,Q2], R = diag [ λ 2 , . . . , λ 2 ] Nu 1 ;
步骤7:计算减温喷水阀门开度控制量u1(k)=u1(k-1)+Δu1(k),和烟气挡板开度控制量u2(k)=u2(k-1)+Δu2(k),式中,u1(k-1)和u2(k-1)分别为减温喷水阀门开度和烟气挡板开度在k-1时刻的对应控制量;
步骤8:根据执行机构限制及实际运行情况,设定控制量约束,减温喷水阀门开度约束上下限为[u1,min,u1,max],烟气挡板开度约束上下限为[u2,min,u2,max],如果ui(k)>ui,max,则令ui(k)=ui,max,Δui(k)=ui,max-ui(k-1);如果ui(k)<ui,min,则令ui(k)=ui,min,Δui(k)=ui,min-ui(k-1),i=1,2;
步骤9:控制系统输出u1(k)和u2(k)到相应执行机构,在之后的每个采样周期内重复执行步骤5到步骤9。
2.根据权利要求1所述通过性能指标实现再热汽温多变量协调预测控制的方法,其特征在于:该预测控制系统预测时域长度N选择上只需大于纯延迟时间,同时控制时域长度Nu取[1,5]之间内的数值。
3.根据权利要求1所述的的一种通过性能指标实现再热汽温多变量协调预测控制的方法,其特征在于:通过调节预测时域长度N来有效调节控制系统的快速性,通过调节控制参数λ1有效调节控制系统的稳定性。
CN201410351300.4A 2014-07-22 2014-07-22 一种通过性能指标实现再热汽温多变量协调预测控制的方法 Expired - Fee Related CN104102134B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410351300.4A CN104102134B (zh) 2014-07-22 2014-07-22 一种通过性能指标实现再热汽温多变量协调预测控制的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410351300.4A CN104102134B (zh) 2014-07-22 2014-07-22 一种通过性能指标实现再热汽温多变量协调预测控制的方法

Publications (2)

Publication Number Publication Date
CN104102134A true CN104102134A (zh) 2014-10-15
CN104102134B CN104102134B (zh) 2016-08-31

Family

ID=51670364

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410351300.4A Expired - Fee Related CN104102134B (zh) 2014-07-22 2014-07-22 一种通过性能指标实现再热汽温多变量协调预测控制的方法

Country Status (1)

Country Link
CN (1) CN104102134B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104571022A (zh) * 2014-12-16 2015-04-29 华润电力湖北有限公司 基于煤耗与可控因子关系的耗差分析模型实验系统及方法
CN104699064A (zh) * 2015-03-20 2015-06-10 东南大学 一种火电机组协调控制系统的煤质自校正预测控制方法
CN104866140A (zh) * 2015-05-27 2015-08-26 小米科技有限责任公司 屏幕校准方法及装置
CN105889910A (zh) * 2016-05-04 2016-08-24 东南大学 一种循环流化床锅炉的新型agc控制方法
CN106287659A (zh) * 2016-08-02 2017-01-04 中国神华能源股份有限公司 再热汽温度控制方法及装置
CN106439786A (zh) * 2016-11-21 2017-02-22 华北电力大学(保定) 电站锅炉再热蒸汽温度的烟气侧和蒸汽侧协调预测控制方法
CN108803342A (zh) * 2018-07-05 2018-11-13 东南大学 一种单元机组负荷快速响应预测控制方法
CN111123704A (zh) * 2019-12-09 2020-05-08 国家能源集团泰州发电有限公司 一种二次再热机组再热汽温的节能型预测控制方法
CN112212712A (zh) * 2020-09-22 2021-01-12 东南大学 一种湿式冷却塔预热式防冻系统的温度控制方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000171003A (ja) * 1998-12-01 2000-06-23 Ishikawajima Harima Heavy Ind Co Ltd 石炭焚火力発電プラントの再熱蒸気温度制御方法
JP2009156553A (ja) * 2007-12-27 2009-07-16 Chugoku Electric Power Co Inc:The 再熱蒸気温度制御方法
JP2011027036A (ja) * 2009-07-27 2011-02-10 Hitachi Ltd 複合発電プラント及びその制御方法
CN102401371A (zh) * 2011-12-06 2012-04-04 东南大学 一种基于多变量预测控制的再热气温优化控制方法
JP2013133988A (ja) * 2011-12-26 2013-07-08 Mitsubishi Heavy Ind Ltd 石炭焚き火力発電プラントの制御装置及び制御方法
US20130333386A1 (en) * 2009-09-08 2013-12-19 General Electric Company Method and apparatus for controlling moisture separator reheater

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000171003A (ja) * 1998-12-01 2000-06-23 Ishikawajima Harima Heavy Ind Co Ltd 石炭焚火力発電プラントの再熱蒸気温度制御方法
JP2009156553A (ja) * 2007-12-27 2009-07-16 Chugoku Electric Power Co Inc:The 再熱蒸気温度制御方法
JP2011027036A (ja) * 2009-07-27 2011-02-10 Hitachi Ltd 複合発電プラント及びその制御方法
US20130333386A1 (en) * 2009-09-08 2013-12-19 General Electric Company Method and apparatus for controlling moisture separator reheater
CN102401371A (zh) * 2011-12-06 2012-04-04 东南大学 一种基于多变量预测控制的再热气温优化控制方法
JP2013133988A (ja) * 2011-12-26 2013-07-08 Mitsubishi Heavy Ind Ltd 石炭焚き火力発電プラントの制御装置及び制御方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
孙德立: "具有Smith预估器再热汽温联合控制系统的研究", 《华北电力学院学报》, no. 2, 31 December 1988 (1988-12-31), pages 87 - 91 *
雎刚,等: "单元机组负荷多变量模型预测控制", 《中国电机工程学报》, vol. 22, no. 4, 30 April 2002 (2002-04-30), pages 144 - 148 *
雎刚,等: "火电机组锅炉过热汽温的约束模型预测控制研究", 《热能动力工程》, vol. 16, no. 6, 30 November 2001 (2001-11-30), pages 641 - 643 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104571022A (zh) * 2014-12-16 2015-04-29 华润电力湖北有限公司 基于煤耗与可控因子关系的耗差分析模型实验系统及方法
CN104699064B (zh) * 2015-03-20 2017-03-08 东南大学 一种火电机组协调控制系统的煤质自校正预测控制方法
CN104699064A (zh) * 2015-03-20 2015-06-10 东南大学 一种火电机组协调控制系统的煤质自校正预测控制方法
CN104866140A (zh) * 2015-05-27 2015-08-26 小米科技有限责任公司 屏幕校准方法及装置
CN105889910A (zh) * 2016-05-04 2016-08-24 东南大学 一种循环流化床锅炉的新型agc控制方法
CN105889910B (zh) * 2016-05-04 2017-11-03 东南大学 一种循环流化床锅炉的新型agc控制方法
CN106287659A (zh) * 2016-08-02 2017-01-04 中国神华能源股份有限公司 再热汽温度控制方法及装置
CN106439786A (zh) * 2016-11-21 2017-02-22 华北电力大学(保定) 电站锅炉再热蒸汽温度的烟气侧和蒸汽侧协调预测控制方法
CN106439786B (zh) * 2016-11-21 2018-05-18 华北电力大学(保定) 电站锅炉再热蒸汽温度的烟气侧和蒸汽侧协调预测控制方法
CN108803342A (zh) * 2018-07-05 2018-11-13 东南大学 一种单元机组负荷快速响应预测控制方法
CN111123704A (zh) * 2019-12-09 2020-05-08 国家能源集团泰州发电有限公司 一种二次再热机组再热汽温的节能型预测控制方法
CN112212712A (zh) * 2020-09-22 2021-01-12 东南大学 一种湿式冷却塔预热式防冻系统的温度控制方法
CN112212712B (zh) * 2020-09-22 2022-03-04 东南大学 一种湿式冷却塔预热式防冻系统的温度控制方法

Also Published As

Publication number Publication date
CN104102134B (zh) 2016-08-31

Similar Documents

Publication Publication Date Title
CN104102134A (zh) 一种通过性能指标实现再热汽温多变量协调预测控制的方法
CN103225799B (zh) 一种火电厂主蒸汽温度控制方法
CN102401371B (zh) 一种基于多变量预测控制的再热汽温优化控制方法
CN104122797B (zh) 一种新型火电机组负荷多变量预测控制方法
CN111443681B (zh) 超临界火电机组协调控制系统多模型预测控制设计方法
CN102494336B (zh) 一种循环流化床锅炉燃烧过程多变量控制方法
Prasad et al. A local model networks based multivariable long-range predictive control strategy for thermal power plants
CN102841539B (zh) 基于多模型预测控制的亚临界协调控制方法
CN106019939B (zh) 基于非线性智能预测控制技术的超临界机组再热汽温优化控制系统
CN102841540A (zh) 基于多模型预测控制的超临界机组协调控制方法
CN104199299B (zh) 燃气轮机负荷调节性能的多变量受限广义预测控制方法
CN102444784A (zh) 基于动态矩阵控制的钢铁企业蒸汽管网压力控制系统
CN105676647A (zh) 一种基于状态观测器的冷热电三联供系统预测控制方法
Zhang et al. Zone economic model predictive control of a coal-fired boiler-turbine generating system
CN105786035A (zh) 基于试探式预测控制技术的火电机组scr脱硝优化控制系统
CN103134046A (zh) 一种火电机组过热汽温两级协调预测控制方法
CN103345153B (zh) 通用网络化预测模糊控制方法
CN101761917A (zh) 一种锅炉过热汽温模糊控制方法
CN107270283B (zh) 一种基于循环流化床机组的多变量约束预测控制方法
CA2943023A1 (en) Feedforward control with intermittent re-initialization based on estimated state information
CN110515304A (zh) 基于ARX-Laguerre函数模型的过热汽温PID预测控制方法
CN104181958A (zh) 一种过热汽温快速节能预测控制方法
CN113835342B (zh) 一种过热汽温系统的抗扰预测控制方法
CN105700357A (zh) 基于多变量pid-pfc的锅炉燃烧系统控制方法
CN106287659A (zh) 再热汽温度控制方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160831

Termination date: 20190722