锅炉过热汽温系统的控制

3. 广义预测控制（GPC ）[1,3,4,6]
3.1 原理简介 广义预测控制是随着自适应控制的研究而发展起来的一种预测控制算法，它在保持了最
小方差自校正的在线辨识，输出预测，最小方差控制的基础上，吸收了 DMC 和 MAC 中的 滚动优化策略，提出了兼具自适应控制和预测控制性能。GPC 基于参数模型，引入了不相 等反馈校正等特征，呈现出优良的控制性能和鲁棒性。 3.1.1 预测模型
高响应速度。
PID 控制器的相位角 ∠Gc = −90o + arctanτ1ω + arctanτ 2ω 。只要参数选择适当，PID
控制器可使相位裕度增加，不但有利于系统稳定，还可减弱系统振荡程度，改善动态性能。 PID 控制所具有的这些功能使得它广泛应用于现代生产中，它的鲁棒性好，即其控制品
质对被控对象特性的变化不太敏感，当被控对象简单或控制要求低时，采用 PID 控制效果 较好。
内扰动及副调节器的任务是快速消除内扰，一般可选用纯比例调节器，这里另
Wa1(s) = 25 。外回路及主调节器的任务是维持过热汽温稳定。在工业上一般采用 PI 和 PID
调节器。
对于锅炉过热汽温这类具有大惯性的对象，采用常规的 PID 调节器，在工作点附近的 小范围内，由于其动态特性近似于线性，有可能控制得较好；但当大范围改变给定值或受外
由于参数模型是最小化模型，需要已知模型结构，但模型参数远比非模型参数要小，减 少了预测控制算法的计算量。为了克服模型的参数失配对输入预测误差的影响，在基于参数 模型的预测控制算法中，引进了自适应控制的在线递推算法估计模型参数，并用估计的参数 取代原模型参数，从而可进行预测算法。由于将自适应控制与预测控制相结合，因而用于过 程参数慢时所引起的预测模型输出误差得以及时修正，从而改善了系统的动态性能。 3.1.2 滚动优化
预测控制是一种闭环控制算法，在通过优化确定了一系列未来的控制作用后，为了防止 模型失配或环境影响引起控制对理想状态的偏离，它通常不是把这些控制作用逐一全部实 施，而只是实现本时刻的控制作用。到下一采样时刻，则首先检测对象的实际输出，并利用 这一实时信息对基于模型的预测进行修正，然后在进行新的优化。
反馈校正的形式是多样的，可以在保持预测模型不变的基础上，并力图在优化时对系统 未来的动态行为做出较准确的预测。因此，预测控制中的优化不仅基于模型，而且应用了反 馈的信息，因而构成了闭环优化。
3.2 GPC 算法
广义预测控制是兼具自适应控制和预测控制性能的一种预测控制。则其原理具有预测控
制的一般特性，但具体实施过程有具有其特点
3.2.1 预测模型
该算法采用自回归积分滑动平均模型(CARIMA)。
GPC 采用具有随机阶跃扰动非平稳噪声的离散差分方程描述，即 CARIMA 模型：
A(z−1) y(k) = z−d B(z−1)u(k ) + C(z−1)ξ /+

锅炉过热汽温系统的控制
刘学智
武汉大学电气工程学院 武汉 430072
Email: newyouthliu@
摘要: 现代锅炉的过热器是在高温、高压条件下工作的，锅炉出口的过热蒸汽汽温对电厂的 安全经济运行有重大影响。本文以锅炉过热汽温控制系统为研究对象，通过仿真曲线比较了 两种 不同的控制 方案(PID 和 GPC)的控 制效果。采 用常规的 PID 控制时 ，本文使用 Matlab /Simulink 对其进行了仿真。 本文重点讨论和分析了 GPC 控制方案，从预测控制的 基本原理出发，以自回归积分滑动平均预测模型(CARIMA)为基础，详细分析了 GPC 控制算法 (预测模型、滚动优化、反馈校正)，并编写出相应的程序，绘制了采用 GPC 控制时的仿真曲 线和计算了动态性能（调整时间和超调量）。在此基础上，比较和分析了两种控制方案的控 制效果，得到了如下结论：GPC 的控制效果比 PID 好，在系统有扰动和被控对象的参数变化 较大(惯性或增益增大)时，其优势愈加凸显，呈现出优良的控制性能和鲁棒性。因此， GPC 所具有的自适应控制和预测控制的优点，使得其更适合于像锅炉过热蒸汽控制这样的大延 迟、大惯性和时变的系统。
界环境（包括工况）太大扰动时，就需要即时修正 PID 参数，否则将达不到控制要求。 预测控制对多容大惯性的对象具有较好的控制效果，它的滚动优化和反馈校正能力可以
动态的补偿模型与对象间的失配误差，具有一定的鲁棒性。
现对主调节器Wa2 分别采用 PID 和 GPC 控制方法，对控制过程进行仿真。
2. 比例积分微分控制（PID）[5]
r+
_
Wa2
u+ d+
_
Wa1
Wo1
Wo2
y
WH1
WH2 图 1 串级控制框图
1

已知 Wo1 ( s )
=
8 (1 + 15s)2
℃/mA
，
Wo 2
(s)
=
1.125 (1+ 25s)3
℃/mA，
WH1(s) = 0.1 ℃/mA ，
WH 2 (s) = 0.1℃/mA
(1.1)
式中，y(t), u(t), ξ (t)分别为系统的输出，输入和均值为零，方差为σ 2 的白噪声, + 为差分
nc
na
∑ ∑ 算子,且 + =1− z−1 , C(z−1) = 1+ ci z−i , A(z−1) = 1+ ai z−i ,
i =1
i =1
nb
nc
∑ ∑ B(z−1) = 1+ bi z−i , C(z−1) = 1+ ci z−i ，d 为系统的滞后。若 d>1，只需令 B(z−1) 多项
绝大多数电厂采用如图所示的具有导前信号的串级控制框图如图 1 所示，W 为过热汽
温设定值，y 为过热汽温，Wa2 (s) 和Wa1(s) 分别为过热汽温控制系统的主调节器和副调节
器；Wa1(s) ，Wo2 (s) 分别为调节对象惰性区及导前区的传递函数；WH1(s) ，WH 2 (s) 分别
为导前汽温和过热汽温的测量单元；d 为扰动。
G j (z−1)
=
g
jቤተ መጻሕፍቲ ባይዱ0
+
g
z −1
j ,1
+"+
g
z−
j, j+nb−1
j −nb+1
这里， deg(Rj) = j −1, deg(Sj) = na , deg(G j) = nb + j −1
现令
G j (z-1)=G j (z-1)+z-jFj (z-1)
Gj (z−1)
=
g
j,0
+
g
z −1
i =1
i =1
式中的前 d-1 项系数为零即可，即 b0 = b1 = b2 = ... = bd −2 = 0 ，现取 d=1，则(1.1)式可写为
3

A(z-1)y(k)=B(z-1)u(k-1)+C(z-1)ξ (t) /+
(1.2)
其可简化为
A(z-1)y(k)=B(z-1)+u(k-1)+C(z-1)ξ (t)
预测控制通过性能指标的最优来确定未来的控制作用。这一性能指标涉及到系统未来的 行为。例如，通常可取被控对象输出在未来的采样点上跟踪某一期望轨迹的方差为最小，但 也可取更广泛的形式，例如要求控制能量为最小而同时保持输出在某一给定范围内等等。性 能指标中涉及到的系统未来的行为，是根据预测模型由未来的控制策略决定的。
预测控制中的优化是一种有限时段的滚动优化。在每个采样时刻，优化性能指标只涉及 到从该时刻起未来有限的时间，而到下一采样时刻，这一优化时段同时向前推移。因此，预 测控制不是用一个对全局相同的优化性能指标，而是在每个时刻有一个相对该时刻的优化性 能指标。不同时间区域则是不同的。因此，在预测控制中，优化不是一次离线进行，而是反 复在线进行，这就是滚动优化的含义，也是预测控制区别于传统最优控制的根本点。 3.1.3 反馈校正
1. 问题背景
现代锅炉的过热器是在高温，高压条件下工作的，锅炉出口的过热蒸汽温度是整个汽水 行程中工质的最高温度，对电厂的安全经济运行有重大影响。目前大型机组过热汽温调节大 多采用喷水减温方案。采用喷水减温调节过热汽温时，一般把过热蒸汽分为两个区域：导前 区和惰性区。从喷水减温的工艺过程可知，以喷水量为输入，以过热蒸汽温度为输出，对象 具有大延迟，大惯性和时变性。过热器管道长度和蒸汽容积较大，当减温水流量发生变化时， 过热器出口蒸汽温度有较大的延迟；负荷变化时，主蒸汽的动态特性变化明显；且主蒸汽温 度还具有分布参数和扰动多的特点。
其中，满足如下的 Diophantine 方程：
C(z−1) = A(z-1)R j (z-1)+z-jSj (z-1)
其中
G j (z-1)=R j (z-1)B(z-1)
Rj (z−1)
=
rj,0
+
rj,1z −1
+"
+
rj,
z−
j −1
j +1
S j (z−1) = s j,0 + s j,1z−1 +" + sna z−na
j ,1
+"
+
g
j,
z−
j −1
j +1
Fj (z−1) =
f j,0
+
f j,1z−1 +" +
f z−nb+1 j ,nb −1
其中， deg(Gj ) = j −1 , deg(Fj) = nb 。 为了方便讨论，取 C(z−1) =1。