电力工程技术  2024, Vol. 43 Issue (3): 254-261  
0

引用本文  

范永胜, 陈卓, 郝勇生, 等. 基于相位补偿的过热汽温自抗扰控制[J]. 电力工程技术, 2024, 43(3): 254-261.
FAN Yongsheng, CHEN Zhuo, HAO Yongsheng, et al. Phase compensation based active disturbance rejection control for superheated steam temperature[J]. Electric Power Engineering Technology, 2024, 43(3): 254-261.

基金项目

国家自然科学基金资助项目(52076037)

作者简介

范永胜(1969), 男, 博士,高级工程师,从事火力发电厂技术研究、基建/生产技术管理工作(E-mail:16010163@ceic.com); 陈卓(1996),男,博士在读,研究方向为热工过程建模与控制; 郝勇生(1978),男,博士,教授,研究方向为企业、园区和大型公建能源管控理论与技术开发、能源动力系统先进测控理论与技术开发。

文章历史

收稿日期:2023-11-30
修回日期:2024-02-27
DOI: 10.12158/j.2096-3203.2024.03.027
文章编号: 2096-3203(2024)03-0254-08   中图分类号: TM62   
基于相位补偿的过热汽温自抗扰控制
范永胜1,2, 陈卓2, 郝勇生2, 苏志刚2, 魏洪轩2, 赵刚2    
1. 国家能源集团江苏电力有限公司,江苏 南京 210036;
2. 东南大学能源与环境学院,江苏 南京 210096
摘要:随着可再生能源接入电网比例的逐步增大,热力发电厂需要应对更加频繁、更大范围的负荷变化,给电厂的高阶大惯性过热汽温过程的控制带来严峻的挑战。为此,文中针对一类高阶大惯性过热汽温过程,提出一种基于相位补偿的自抗扰控制(phase compensation based active disturbance rejection control, PC-ADRC)方法。首先,阐述过热汽温系统的工作原理和控制难点。然后,采用低频近似法详细推导相位补偿(phase compensation, PC)网络模型,提出采用PC网络对模型动态特性进行补偿,得到等效降阶模型的简化思路。为便于工程应用,给出PC-ADRC系统的简单实现方法和等效模型分析。最后,对PC-ADRC系统的稳定性和鲁棒性进行研究。理论分析和仿真结果表明,所提出的PC-ADRC系统能有效提升高阶过程控制系统的鲁棒性和快速响应能力。
关键词过热汽温    火力发电    基于相位补偿的自抗扰控制(PC-ADRC)    高阶过程    低频近似    鲁棒性    
0 引言

近年来,由于可再生能源的清洁和可持续利用优势,其在电网中应用程度显著提高。然而,可再生能源,如风能和太阳能固有的不确定性和间歇性,对电网的电力平衡构成严峻挑战。因此,通过提高传统火电机组的调频能力,实时平衡电网功率偏差已成为一种新的能源系统发展趋势[1-2]。过热汽温系统参数作为火电机组工质的峰值温度和重要过程参数,对机组运行的安全性和经济性有着重要影响。通常过热汽温的控制偏差要求在±5 ℃[3]。过高的过热汽温会导致金属受损甚至机组异常停机[4],而较低的汽温会降低机组循环效率,增加机组耗煤[5]。一般来说,大范围的负荷变化会导致过大的过热汽温偏差[6]。因此,在机组宽负荷灵活运行的背景下,过热汽温系统的控制性能需要进一步提高。

然而,过热汽温过程是一个典型的高阶大惯性过程,这也是许多工业过程共同具有的特点,如再热汽温[7]和脱硫系统[8]。这个特点使得过热汽温的控制性能在工程实践中难以提高。为此,研究者提出许多先进的控制方法,如模糊控制[9]、分数阶比例积分微分(proportional integral derivative,PID)控制[10]和模糊模型预测控制(model predictive control,MPC)[11]。然而,这些先进的控制方法过多依赖精确的数学模型或计算复杂度较高,使得其在实际应用中还存在一定的限制。

自抗扰控制(active disturbance rejection control, ADRC)方法具有处理模型不确定性、非线性和外部扰动的强大能力[12],因此在控制界引起广泛关注[13]。ADRC继承PID的许多优点,如建模要求低、结构简单等,因此在工程实践中得到许多应用,如工业过程控制[5]、飞行器控制等[14]。一些现场应用中存在的实际问题也得到有效解决,如无扰切换、抗积分饱和[15]以及基于传递函数的离散实现[16]等。

然而,一般ADRC器中阶次和模型的阶次严格相关,一般ADRC器中阶次高于模型阶次一阶[17]。因此,对于高阶过程,会得到高阶的复杂控制器,不利于实际应用。在这种情况下,降阶ADRC器和低阶ADRC器成为有前景的研究方向,并得到广泛的研究。文献[18]证明低阶ADRC器可以有效镇定高阶开环稳定对象。文献[19]从考虑控制器的相位裕度和增益裕度的角度出发,为降阶自抗扰器开发基于编程的软件调参方法。文献[20]细致研究了低阶ADRC器控制高阶过程时的鲁棒调参问题。针对高阶系统,研究人员也提出一些改进的ADRC器。文献[21]通过增加输入通道补偿,降低主控制器阶数。文献[22]提出在一阶ADRC器中加入一个微分环节以提高高阶系统对输入信号的响应速度。然而,这些改进的控制器几乎不可避免地会导致新的问题,例如参数重调。事实上,从相位滞后系统中获得更多实时准确信息是提高控制质量的关键[23-24]

为此,文中提出一种基于相位补偿的ADRC(phase compensation based ADRC, PC-ADRC)器,该控制器适用于高阶大惯性过热汽温系统。文中从过程模型动态补偿的角度推导相位补偿(phase compensation, PC)网络模型,给出控制器参数定量整定方法。在PC-ADRC的控制系统稳定性和鲁棒性分析基础上,通过仿真验证所提出方法的有效性。

1 系统模型和ADRC器 1.1 过热汽温控制系统

图 1为典型燃煤机组锅炉系统结构,该系统可分为风烟系统和蒸汽系统。在烟气侧,原煤经给煤机输送至磨煤机磨成煤粉,煤粉由一次风输送至锅炉。在锅炉内,煤粉进一步与二次风混合燃烧,产生高温烟气。烟气与金属壁面进行热交换,加热其中的水或蒸汽。在蒸汽侧,给水经省煤器加热后进入汽包,在锅炉的水冷壁进一步吸热蒸发为水蒸汽。蒸汽在过热器进一步加热为主蒸汽,最后进入汽轮机膨胀做功并驱动发电机发电。

图 1 典型燃煤机组锅炉系统结构 Fig. 1 Structure of a typical coal-fired unit boiler system

过热汽温控制系统结构如图 2所示,燃煤机组通常采用减温器喷水方式来调节过热汽温。给水泵的减温喷水与高温蒸汽直接混合,可以有效降低蒸汽温度。通过过热汽温的测量和控制机构形成一个闭环控制系统,在该系统中由阀门开度来控制喷淋水的流量从而控制过热汽温。

图 2 过热汽温控制系统结构 Fig. 2 Structure of superheated steam temperature control system
1.2 控制难点

过热汽温控制系统是火电机组中最重要的控制系统之一,也是典型的高阶大惯性过程[20],常采用串级控制结构。由文献[5]中研究可知串级系统的副回路可以整定为快速随动系统,因此控制系统可以等效为主控制器和过热汽温惰性区的单回路系统。文中过热汽温惰性区传递函数模型[25]如式(1)所示。

$ G_{\mathrm{p}}(s)=\frac{1.276}{(1+18.4 s)^6} $ (1)

式中:s为传递函数模型中的拉普拉斯算子。

首先,由于过热汽温系统的复杂性,其精确的过程模型难以获取。在实际工程中,一般通过阶跃响应实验来确定其近似数学模型。其次,过热汽温系统受到负荷变化、煤质变化和燃烧不稳定性等多种扰动的影响。此外,过热汽温过程是一个典型的高阶大惯性过程。虽然减温喷水是最常用的过热汽温调节方式,但在喷水阀门开度变化后,出口温度仍需5 min左右才能达到稳定。因此,文中旨在提出一种实用的控制方案来处理过热汽温控制系统内部和外部的不确定性,以提高其系统的快速响应能力。

1.3 典型高阶过程ADRC器

一个典型ADRC器六阶过程模型的表示如式(2)所示。

$ G_{\mathrm{p}}(s)=\frac{K}{(1+T s)^6} $ (2)

式中:K为增益系数;T为时间常数。过程模型的状态空间表示如式(3)所示。

$\left\{\begin{array}{l} \dot{x}_1=x_2 \\ \dot{x}_2=x_3 \\ \dot{x}_3=x_4 \\ \dot{x}_4=x_5 \\ \dot{x}_5=x_6 \\ \dot{x}_6=f+b u \\ y=x_1 \end{array}\right. $ (3)

式中:uxi(i=1, 2, 3, 4, 5, 6)、y分别为ADRC器控制输出量、状态量和过程输出量;$ \dot{x}_i$xi的导数;b为输入通道增益;f为总扰动,不仅包括外部扰动,还包括未知的内部动态。

因此,基于式(3),扩张状态观测器(extended state observer,ESO)被设计为:

$\left\{\begin{array}{l} \dot{z}_1=z_2+\beta_1\left(y-z_1\right) \\ \dot{z}_2=z_3+\beta_2\left(y-z_1\right) \\ \dot{z}_3=z_4+\beta_3\left(y-z_1\right) \\ \dot{z}_4=z_5+\beta_4\left(y-z_1\right) \\ \dot{z}_5=z_6+\beta_5\left(y-z_1\right) \\ \dot{z}_6=z_7+b_0 u+\beta_6\left(y-z_1\right) \\ \dot{z}_7=\beta_7\left(y-z_1\right) \end{array}\right. $ (4)

式中:zi$ \dot{z}_i$分别为观测器各个状态变量的估计和其导数;z7$ \dot{z}_7$分别为观测器总扰动的估计值和其导数;β1β7为观测器增益;b0为输入通道增益b的估计值。

已有文献证明当观测器增益取值合理时,各个估计值可以渐进跟踪系统状态量和总扰动[26-27]。通过对总扰动进行前馈估计补偿,过程模型可以被近似为串联积分器。将状态反馈控制器设计为:

$ u=\frac{\overline{\boldsymbol{K}}(\overline{\boldsymbol{r}}-\boldsymbol{z})}{b_0} $ (5)

式中:K=[k1 k2 k3 k4 k5 k6 1],ki为控制器增益;r=[r 0 0 0 0 0 0]Tr为设定值;z=[z1 z2 z3 z4 z5 z6 z7]T

可以看出,常规设计所得的ADRC器为复杂的高阶控制器,而低阶控制器因其结构简单、易于实现而在实际应用中更为广泛。因此,文中将设计一种PC-ADRC器。

2 PC-ADRC 2.1 PC

在高阶大惯性过热汽温系统的控制实践中,应先对原系统动态进行估计或补偿,进而对补偿后的简化系统设计低阶控制器,从而提升原复杂系统的控制性能,例如著名的Smith预估器。基于此,文中提出一种PC-ADRC,其控制原理如图 3所示。其中d为外部扰动;y′为PC输出;$ \hat{f}$为总扰动的估计值;u0为状态反馈输出;kpkd分别为控制器比例和微分增益。

图 3 PC-ADRC原理 Fig. 3 PC-ADRC principle

该控制结构与传统的二阶ADRC器的主要区别在于PC网络,PC网络模型可以表示为:

$ P(s)=\frac{a^2 \tau^2 s^2+2 \xi a \tau s+1}{\tau^2 s^2+2 \xi \tau s+1} $ (6)

式中:a为PC网络的初始放大系数;ξ为阻尼比;τ为补偿器时间常数,是最重要的参数,与系统的阶次以及时间常数T相关。

PC的主要目的是对过程动态进行补偿,即在低频段进行等效来降低系统阶次:

$ \begin{aligned} & G_{\mathrm{p}}(s) P(s)=\frac{K}{(1+T s)^n} \times \frac{a^2 \tau^2 s^2+2 \xi a \tau s+1}{\tau^2 s^2+2 \xi \tau s+1} \underset{s \rightarrow 0}{\approx} \\ & \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ G_{\mathrm{p}}^{\prime}(s)=\frac{K^{\prime}}{\left(1+T^{\prime} s\right)^{n-2}} \end{aligned} $ (7)

式中:Gp(s)为补偿后的过程模型(等效模型);K′、T′分别为补偿后模型的增益系数和时间常数;n为原过程模型阶次。可以在简化模型Gp(s)的基础上设计低阶ADRC器,与1.3节类似,二阶ADRC器可以设计如式(8)所示。

$ u=\frac{u_0-z_3}{b_0}=\frac{k_{\mathrm{p}}\left(r-z_1\right)-k_{\mathrm{d}} z_2-z_3}{b_0} $ (8)

值得注意的是,由于精确的低频近似,二阶ADRC器在低频范围内可以有效控制相应阶次的模型。因此,可以有效提高系统的控制性和鲁棒性。

至此,需要整定的控制器参数包括PC网络参数aξτ,二阶ADRC器的观测器增益β1、β2β3和控制器增益kpkd。在进行ADRC器调参前,需要先设计PC网络。

2.2 PC网络设计

如2.1节所述,PC网络有3个待设计参数aξτ,并且其都具有特定的物理含义。经过大量的测试,a建议取值区间为[2.6,3.5],以获得合理的放大性能与噪声,文中为推导简明,a取3;ξ值越小则系统振荡越明显,因此ξ取0.5,其相应的衰减率为0.973(其他aξ取值并不显著影响后续推导过程)。

由式(6)、式(7)以及a=3、ξ=0.5可得:

$ \frac{K}{(1+T s)^n} \times \frac{a^2 \tau^2 s^2+2 \xi a \tau s+1}{\tau^2 s^2+2 \xi \tau s+1} \approx \frac{K^{\prime}}{\left(1+T^{\prime} s\right)^{n-2}} $ (9)

其静态增益为:

$ \begin{gathered} \lim _{s \rightarrow 0} \frac{K}{(1+T s)^n} \times \frac{a^2 \tau^2 s^2+2 \xi a \tau s+1}{\tau^2 s^2+2 \xi \tau s+1}= \\ \lim _{s \rightarrow 0} \frac{K^{\prime}}{\left(1+T^{\prime} s\right)^{n-2}} \end{gathered} $ (10)

由式(10)可得K=K′。进一步有:

$ \begin{gathered} \left(1+T^{\prime} s\right) n-2\left(9 \tau^2 s^2+3 \tau s+1\right) \approx \\ (1+T s) n\left(\tau^2 s^2+\tau s+1\right) \end{gathered} $ (11)

考虑二阶泰勒近似可得:

$ \begin{gathered} {\left[1+(n-2) T^{\prime} s+\frac{(n-2)(n-3)}{2} T^{\prime 2} s^2\right] \times} \\ \left(9 \tau^2 s^2+3 \tau s+1\right) \approx \\ {\left[1+n T s+\frac{n(n-1)}{2} T^2 s^2\right] \times\left(\tau^2 s^2+\tau s+1\right)} \end{gathered} $ (12)

由前2项系数相等推导可得:

$ \left\{\begin{array}{l} 3 \tau+(n-2) T^{\prime}=\tau+n T \\ \left.9 \tau^2+\frac{(n-2)(n-3)}{2} T^{\prime 2}+3 \tau(n-2) \right\rvert\, T^{\prime}= \\ \ \ \ \ \ \ \ \ \ \ \ \tau^2+\frac{n(n-1)}{2} T^2+n \tau T \end{array}\right. $ (13)

求解τT′可得:

$\left\{\begin{array}{l} \tau=\frac{T}{2(2 n-5)}[\sqrt{5 n(n-2)}-n] \\ T^{\prime}=\frac{n T-2 \tau}{n-2} \end{array}\right. $ (14)
2.3 PC网络实现

文中固定式(6)所示PC模型的一个重要原因是其可以通过组态软件轻松实现,具有可移植性和广泛适用性。2.2节中得到的PC网络可以通过简单的增益块和积分块来实现和封装。

PC网络如图 4所示,Ki为PC网络增益系数;T1T2为积分器的时间常数;Ein(s)为输入;Eout(s)为输出;E1(s)、E2(s)均为中间变量。

图 4 PC网络 Fig. 4 PC network

图 4可得:

$\left\{\begin{array}{l} E_1(s)=K_1 K_3 E_{\text {out }}(s)-K_2 K_3 E_{\text {in }}(s) \\ E_2(s)=E_{\text {out }}(s)\left(\frac{1}{T_1 T_2 s^2}+\frac{K_1 K_3}{T_2 s}\right)- \\ \ \ \ \ \ \ \ E_{\text {in }}(s)\left(\frac{1}{T_1 T_2 s^2}+\frac{K_2 K_3}{T_2 s}\right) \\ E_{\text {out }}(s)=K_4 K_6 E_{\text {in }}(s)-K_5 K_6 E_2(s) \end{array}\right. $ (15)

进一步可得:

$ P(s)=\frac{E_{\text {out }}(s)}{E_{\text {in }}(s)}=\frac{\frac{K_4}{K_5} T_1 T_2 s^2+K_2 K_3 T_1 s+1}{\frac{1}{K_5 K_6} T_1 T_2 s^2+K_1 K_3 T_1 s+1} $ (16)

$ \frac{K_4}{K_5} $=9、K5K6=1、K2K3=3、K1K3=1、T1=T2=τ,可得a=3、ξ=0.5时的PC模型。

3 鲁棒性分析与参数整定 3.1 鲁棒性分析

PC-ADRC的等效二自由度结构如图 5所示,图中,Gc(s)为等效反馈控制器模型;r(s)、e(s)、d(s)、u(s)、y(s)分别为参考信号r、误差信号e、扰动信号d、输出信号u、被控变量y的拉普拉斯变换;GF(s)为等效前馈控制器传递函数。控制系统的开环传递函数Gop(s)应为:

$\begin{gathered} G_{\mathrm{op}}(s)=G_{\mathrm{c}}(s) G_{\mathrm{p}}(s) P(s)= \\ \frac{\left(k_{\mathrm{p}} \beta_1+k_{\mathrm{d}} \beta_2+\beta_3\right) s^2+\left(k_{\mathrm{p}} \beta_2+k_{\mathrm{d}} \beta_3\right) s+k_{\mathrm{p}} \beta_3}{b_0\left[s^3+\left(k_{\mathrm{d}}+\beta_1\right) s^2+\left(k_{\mathrm{p}}+\beta_2+k_{\mathrm{d}} \beta_1\right) s\right]} \times \\ \frac{K}{(1+T s)^6} \times \frac{9 \tau^2 s^2+3 \tau s+1}{\tau^2 s^2+\tau s+1} \end{gathered} $ (17)
图 5 PC-ADRC的等效二自由度结构 Fig. 5 The equivalent two-degree freedom structure of PC-ADRC

可得闭环控制系统的灵敏度函数S(s)和互补灵敏度函数T(s):

$ \left\{\begin{array}{l} S(s)=1 /\left(1+G_{\text {op }}(s)\right) \\ T(s)=G_{\text {op }}(s) /\left(1+G_{\text {op }}(s)\right) \end{array}\right. $ (18)

为便于使用,选取灵敏度函数和互补灵敏度函数的峰值函数MSMT作为鲁棒性指标[28]

$ \left\{\begin{array}{l} M_{\mathrm{S}}=\max |S(s)| \\ M_{\mathrm{T}}=\max |T(s)| \end{array}\right. $ (19)
3.2 控制器参数整定

文中采用带宽参数法[29]以进一步减少控制器调参数量:

$\left\{\begin{array}{l} \beta_1=3 \omega_0 \\ \beta_2=3 \omega_0^2 \\ \beta_3=\omega_0^3 \\ k_{\mathrm{p}}=\omega_{\mathrm{c}}^2 \\ k_{\mathrm{d}}=2 \xi_{\mathrm{c}} \omega_{\mathrm{c}} \end{array}\right. $ (20)

式中:ω0ωc分别为观测器和闭环系统带宽;ξc为系统阻尼比,常取1。至此,PC-ADRC的待调参数减少为3个,即ω0ωcb0

ωc的选择应结合过程信息,如系统阶次和时间常数。参考文献[30-31],文中ωc取值为:

$\omega_{\mathrm{c}}=\frac{10}{\sigma n^{\prime} T^{\prime}} $ (21)

式中:σ为期望的设定值跟踪时间因子,对于高阶过程,常取1~4[20]n′为补偿后的系统阶次。

观测器带宽ω0决定了被观测和补偿的总扰动速度。因此,从控制性能提升角度而言,较大的ω0值更好。但在实际过程中ω0的取值是受限的,如受采样率的限制。采用经典的比例关系,文中ω0的取值为:

$ \omega_0=10 \omega_{\mathrm{c}} $ (22)

如前所述,对于ADRC器,b0是一个重要参数,尤其在采用低阶ADRC高阶过程时,b0过小会导致控制系统不稳定,而b0太大又会使得控制性能下降[31]。文献[20]针对K/(1+Ts)n型高阶过程ADRC器提出关于b0的参数调整方法。但是,由于新引入的PC网络,这种定量关系不能直接使用。文中将使用鲁棒回路成形法确定b0

图 3可知,PC-ADRC系统的闭环传递函数可以近似为:

$ \begin{gathered} G_{\mathrm{cl}}(s)=\frac{Y(s)}{R(s)} \approx \frac{k_{\mathrm{p}}}{s^2+k_{\mathrm{d}} s+k_{\mathrm{p}}} \times \frac{\tau^2 s^2+\tau s+1}{9 \tau^2 s^2+3 \tau s+1}= \\ \frac{\omega_{\mathrm{c}}^2}{s^2+2 \xi_{\mathrm{c}} \omega_{\mathrm{c}} s+\omega_{\mathrm{c}}^2} \times\left(\frac{1}{9} s^2+\frac{1}{9 \tau} s+\frac{1}{9 \tau^2}\right) \div \\ \left(s^2+\frac{1}{3 \tau} s+\frac{1}{9 \tau^2}\right) \end{gathered} $ (23)

式中:R(s)、Y(s)分别为设定值和闭环输出。可见,最终等式中的第二项完全由PC网络决定。因此,可以通过选择ωcξc来调节闭环系统响应。例如,当ωc和1/(3τ)比较接近时,可以选择较大的阻尼比ξc,而当ωc为主导极点时,第二项的动态影响很小,可以忽略不计。

4 仿真试验及结果分析

文中将通过仿真验证过热汽温系统低频近似的准确性。等效模型Gp(s)P(s)和理想模型Gp(s)频率响应见图 6,可以看出,等效模型和理想模型在低频段的响应高度一致,即很好地反映其低频动态。

图 6 等效模型和理想模型频率响应比较 Fig. 6 Frequency response comparison between the equivalent model and the ideal model

通过过热汽温系统控制仿真验证文中所提出方法的有效性,仿真中与以下3种不同的控制器进行设定值跟踪和负荷扰动仿真对比,分别是:

(1) 基于鲁棒调参方法的二阶ADRC器[20]

(2) 简单内模控制-比例积分(simple internal model control-proportion integral, SIMC-PI)[32]控制器;

(3) PC-比例积分(proportion integral,PI),即PC-PI控制器,即在SIMC-PI控制器中加入PC,以验证PC网络的适用性。

仿真中,在150 s和1 500 s分别加入设定值单位阶跃扰动和阶跃负载扰动。

为对比公平,PC-ADRC、ADRC和SIMC-PI采用相同的鲁棒性约束,具体控制器参数和MSMT指标见表 1,其中,KPKI分别为PI控制器的比例和积分增益。而PC-PI则是在SIMC-PI控制系统输出通道直接加入PC以验证其适用性。PC-ADRC中ξc取值为1.1,b0图 7所示的鲁棒回路成形方法确定。

表 1 控制器参数与鲁棒性指标 Table 1 Controller parameters and robustness metrics

图 7 鲁棒回路成形法整定b0示意 Fig. 7 obustness loop shaping method for tuning b0

仿真结果如图 8所示,其中,被控输出为过热汽温;控制信号为阀门开度;Set-point为参考信号。结果表明PC-ADRC与SIMC-PI有相近的跟踪性能。单位阶跃扰动跟踪过程的误差绝对值积分(integral absolute error,IAE)如表 2所示(IAE1),与其他3种控制器相比,PC-ADRC器显著提高控制系统的抗扰性能。阶跃负载扰动过程的IAE2指标也如表 2所示,这也表明PC可以有效地补偿过热汽温系统的大惯性动态。为进一步验证所提出的控制方法对一般扰动的有效性,考虑加入斜率为0.01斜坡形式和周期为50 s正弦形式的扰动,如图 9所示,仿真结果如图 10所示。可以看出,PC-ADRC器对这2类扰动也表现出很好的抗扰性能,如更快的响应和更小的扰动偏差峰值。斜坡和正弦扰动过程的IAE指标如表 2所示(IAE3和IAE4)。

图 8 不同控制器控制效果对比 Fig. 8 Control effect comparison of different controllers

表 2 扰动过程IAE指标对比 Table 2 IAE indicators comparison of disturbance process

图 9 斜坡和正弦扰动示意 Fig. 9 Slope and sinusoidal disturbance schematic

图 10 斜坡和正弦扰动控制效果对比 Fig. 10 Control effect comparison of slope and sinusoidal disturbance

相同参数下PC-ADRC与ADRC鲁棒性指标对比如图 11所示,其中Gm, PC-ADRCGm, ADRC分别为PC-ADRC器和ADRC器的幅值裕度绝对值;Pm, PC-ADRCPm, ADRC分别为PC-ADRC器和ADRC器的相位裕度绝对值。当采用同阶ADRC器和相同的控制器参数时(b0=-0.045、ω0=0.423、ωc=0.042 3),PC-ADRC器有着更好的鲁棒性。这说明PC可以直接加入到原ADRC系统中,作为控制器补丁,而不会使控制器不稳定。

图 11 相同参数下PC-ADRC与ADRC鲁棒性指标对比 Fig. 11 Robustness indicators comparison for the PC-ADRC and ADRC under the same parameters
5 结论

文中针对高阶大惯性过热汽温过程,提出一种PC-ADRC策略,可显著提高控制系统的抗扰性能。在介绍过热汽温系统的基础上,阐述模型不确定性、外部扰动和高阶大惯性是过热汽温系统的控制难点所在。

为提高控制性能,从理论上推导PC模型,分析PC后控制系统的动态特性,并给出其参数的整定方法。通过仿真实验验证所提出控制方法的有效性。结果表明,PC-ADRC器相较传统PI控制器和ADRC器有更好的控制性能。具体而言,在相近的鲁棒性约束下,相较于SIMC-PI器,所提出的控制方法在单位阶跃扰动值、阶跃负载扰动、斜坡负载扰动和正弦扰动下的IAE1—IAE4性能指标分别提升12%、51%、51%和8.5%,表明其具有更好的设定值跟踪和扰动抑制性能。

进一步地,仿真中通过直接在原PI控制器回路中加入PC模块,可以有效提升控制系统的抗扰性能,相较于SIMC-PC控制器,其在阶跃负载扰动、斜坡负载扰动和正弦扰动下IAE2—IAE4性能指标分别提升5.5%、5.1%和3.4%。表明PC模型不仅适用于文中的ADRC器,而且适用于其他基于误差的控制器,如PID控制器。

参考文献
[1]
李卫华. 提高燃煤机组调频调峰性能的综合控制技术[J]. 热力发电, 2014, 43(11): 6-11, 23.
LI Weihua. Integrated control technology for improving PFC and AGC performance of coal-fired units: research and application[J]. Thermal Power Generation, 2014, 43(11): 6-11, 23. DOI:10.3969/j.issn.1002-3364.2014.11.02 (0)
[2]
陈作舟, 余浩, 王盼盼, 等. 海上风电集群与火电打捆外送系统短路比定义及影响因素分析[J]. 发电技术, 2022, 43(2): 207-217.
CHEN Zuozhou, YU Hao, WANG Panpan, et al. Definition and influencing factors of short-circuit ratio between offshore wind power cluster and thermal power bundling system[J]. Power Generation Technology, 2022, 43(2): 207-217. (0)
[3]
MAZALAN N A, MALEK A A, WAHID M A, et al. Review of control strategies employing neural network for main steam temperature control in thermal power plant[J]. Jurnal Teknologi, 2014, 66(2): 73-76. (0)
[4]
LIU J Z, YAN S, ZENG D L, et al. A dynamic model used for controller design of a coal fired once-through boiler-turbine unit[J]. Energy, 2015, 93: 2069-2078. DOI:10.1016/j.energy.2015.10.077 (0)
[5]
SUN L, HUA Q S, SHEN J, et al. Multi-objective optimization for advanced superheater steam temperature control in a 300 MW power plant[J]. Applied Energy, 2017, 208: 592-606. DOI:10.1016/j.apenergy.2017.09.095 (0)
[6]
CHAN K H, DOZAL-MEJORADA E J, CHENG X, et al. Predictive control with adaptive model maintenance: application to power plants[J]. Computers & Chemical Engineering, 2014, 70: 91-103. (0)
[7]
孙明, 董泽. 二次再热机组汽温系统的串级自抗扰预测控制[J]. 中国电机工程学报, 2019, 39(18): 5415-5423, 5589.
SUN Ming, DONG Ze. Cascade active disturbance rejection predictive control for superheated steam temperature of ultra-supercritical double-reheat unit[J]. Proceedings of the CSEE, 2019, 39(18): 5415-5423, 5589. (0)
[8]
YI H, SU Z G, WANG P H. Closed-loop shaping-based active disturbance rejection control with delay robustness and application to wet desulfurization process[J]. IEEE Transactions on Industrial Electronics, 2024, 71(6): 6117-6127. DOI:10.1109/TIE.2023.3290252 (0)
[9]
杨新民. 智能控制技术在火电厂应用研究现状与展望[J]. 热力发电, 2018, 47(7): 1-9.
YANG Xinmin. Application status and prospect of intelligent control technology in thermal power plants[J]. Thermal Power Generation, 2018, 47(7): 1-9. (0)
[10]
郭伟, 程晓冲. FOPID-GPC算法在电厂过热蒸汽温度控制中的应用[J]. 热力发电, 2012, 41(3): 37-41.
GUO Wei, CHENG Xiaochong. Application of fractional-order PID type generalized predictive control algorithm in superheated steam control system of power plant[J]. Thermal Power Generation, 2012, 41(3): 37-41. DOI:10.3969/j.issn.1002-3364.2012.03.037 (0)
[11]
WU X, SHEN J, LI Y G, et al. Fuzzy modeling and predictive control of superheater steam temperature for power plant[J]. ISA Transactions, 2015, 56: 241-251. DOI:10.1016/j.isatra.2014.11.018 (0)
[12]
XUE W C, HUANG Y. Performance analysis of active disturbance rejection tracking control for a class of uncertain LTI systems[J]. ISA Transactions, 2015, 58: 133-154. DOI:10.1016/j.isatra.2015.05.001 (0)
[13]
韩京清. 自抗扰控制技术: 估计补偿不确定因素的控制技术[M]. 北京: 国防工业出版社, 2008.
HAN Jingqing. Active disturbance rejection control technique[M]. Beijing: National Defense Industry Press, 2008. (0)
[14]
付京博, 邵会兵, 詹韬. 基于深度强化学习的飞行器自抗扰控制技术[J]. 计算机仿真, 2022, 39(10): 54-59.
FU Jingbo, SHAO Huibing, ZHAN Tao. Aircraft active disturbance rejection control technology based on deep reinforcement learning[J]. Computer Simulation, 2022, 39(10): 54-59. (0)
[15]
SUN L, LI D H, HU K T, et al. On tuning and practical implementation of active disturbance rejection controller: a case study from a regenerative heater in a 1 000 MW power plant[J]. Industrial & Engineering Chemistry Research, 2016, 55(23): 6686-6695. (0)
[16]
HERBST G. Transfer function analysis and implementation of active disturbance rejection control[J]. Control Theory and Technology, 2021, 19(1): 19-34. DOI:10.1007/s11768-021-00031-5 (0)
[17]
XUE W C, HUANG Y. On performance analysis of ADRC for nonlinear uncertain systems with unknown dynamics and discontinuous disturbances[C]//Proceedings of the 32nd Chinese Control Conference. Xi'an, China. IEEE, 2013: 1102-1107. (0)
[18]
ZHAO C Z, HUANG Y. Design ADRC for two special kinds of systems[C]//Proceedings of the 30th Chinese Control Conference. Yantai, China. IEEE, 2011: 229-234. (0)
[19]
NOWAK P, CZECZOT J, KLOPOT T. Robust tuning of a first order reduced active disturbance rejection controller[J]. Control Engineering Practice, 2018, 74: 44-57. DOI:10.1016/j.conengprac.2018.02.001 (0)
[20]
HE T, WU Z L, LI D H, et al. A tuning method of active disturbance rejection control for a class of high-order processes[J]. IEEE Transactions on Industrial Electronics, 2020, 67(4): 3191-3201. DOI:10.1109/TIE.2019.2908592 (0)
[21]
WU Z L, HE T, LI D H, et al. Superheated steam temperature control based on modified active disturbance rejection control[J]. Control Engineering Practice, 2019, 83: 83-97. DOI:10.1016/j.conengprac.2018.09.027 (0)
[22]
WU Z L, LI D H, HE T, et al. A comparison study of a high order system with different ADRC control strategies[C]//2018 37th Chinese Control Conference (CCC). Wuhan, China. IEEE, 2018: 1-6. (0)
[23]
高志强. 浅谈工程控制的信息问题[J]. 系统科学与数学, 2016, 36(7): 908-923.
GAO Zhiqiang. On the problem of information in engineering cybernetics[J]. Journal of Systems Science and Mathematical Sciences, 2016, 36(7): 908-923. (0)
[24]
WEI W, ZHANG Z Y, ZUO M. Phase leading active disturbance rejection control for a nanopositioning stage[J]. ISA Transactions, 2021, 116: 218-231. DOI:10.1016/j.isatra.2021.01.004 (0)
[25]
常太华, 江清潘, 朱红路. 火电厂过热汽温系统内模控制研究[J]. 华东电力, 2009, 37(12): 2077-2080.
CHANG Taihua, JIANG Qingpan, ZHU Honglu. Study of internal model control for power plant superheated steam temperature[J]. East China Electric Power, 2009, 37(12): 2077-2080. (0)
[26]
陈增强, 孙明玮, 杨瑞光. 线性自抗扰控制器的稳定性研究[J]. 自动化学报, 2013, 39(5): 574-580.
CHEN Zengqiang, SUN Mingwei, YANG Ruiguang. On the stability of linear active disturbance rejection control[J]. Acta Automatica Sinica, 2013, 39(5): 574-580. (0)
[27]
WU Z H, ZHOU H C, GUO B Z, et al. Review and new theoretical perspectives on active disturbance rejection control for uncertain finite-dimensional and infinite-dimensional systems[J]. Nonlinear Dynamics, 2020, 101(2): 935-959. DOI:10.1007/s11071-020-05845-7 (0)
[28]
SU Z G, SUN L, XUE W C, et al. A review on active disturbance rejection control of power generation systems: Fundamentals, tunings and practices[J]. Control Engineering Practice, 2023, 141: 105716. DOI:10.1016/j.conengprac.2023.105716 (0)
[29]
GAO Z Q. Scaling and bandwidth-parameterization based controller tuning[C]//Proceedings of the 2003 American Control Conference. Denver, CO, USA. IEEE, 2003: 4989-4996. (0)
[30]
SUN Y G, SU Z G, SUN L, et al. Time-delay active disturbance rejection control of wet electrostatic precipitator in power plants[J]. IEEE Transactions on Automation Science and Engineering, 2023, 20(4): 2748-2760. DOI:10.1109/TASE.2022.3201539 (0)
[31]
ZHAO C Z, LI D H. Control design for the SISO system with the unknown order and the unknown relative degree[J]. ISA Transactions, 2014, 53(4): 858-872. DOI:10.1016/j.isatra.2013.10.001 (0)
[32]
SKOGESTAD S. Simple analytic rules for model reduction and PID controller tuning[J]. Journal of Process Control, 2003, 13(4): 291-309. DOI:10.1016/S0959-1524(02)00062-8 (0)
Phase compensation based active disturbance rejection control for superheated steam temperature
FAN Yongsheng1,2, CHEN Zhuo2, HAO Yongsheng2, SU Zhigang2, WEI Hongxuan2, ZHAO Gang2    
1. China Energy Jiangsu Power Co., Ltd., Nanjing 210036, China;
2. School of Energy and Environment, Southeast University, Nanjing 210096, China
Abstract: With the continuous penetration of renewable energies, thermal power plants are facing more frequent and wider variable load operation. The control of high order, large inertia superheated steam temperature processes in power plants is facing a great challenge. Therefore, a phase compensation based active disturbance rejection control (PC-ADRC) strategy for a class of high-order, large-inertia processes is proposed. Firstly, the working principles and control problems of SST system are illustrated. Then, the phase compensation (PC) network model is derived carefully by low-frequency approximation. A model simplification method is presented, in which the PC is adopted to compensate process dynamics and obtain reduced-order plant equivalently. To facilitate engineering applications, its simple implementation and equivalent model analysis for the PC-ADRC system are also given. At last, the stability and robustness of the PC-ADRC control system are theoretically studied. Theoretical analysis and simulation results show that the proposed controller can effectively improve the robustness and fast response ability of the high-order process control system.
Keywords: superheated steam temperature    thermal power generation    phase compensation based active disturbance rejection control (PC-ADRC)    high order process    low-frequency approximation    robustness