CN105700357B - 基于多变量pid-pfc的锅炉燃烧***控制方法 - Google Patents
基于多变量pid-pfc的锅炉燃烧***控制方法 Download PDFInfo
- Publication number
- CN105700357B CN105700357B CN201610109217.5A CN201610109217A CN105700357B CN 105700357 B CN105700357 B CN 105700357B CN 201610109217 A CN201610109217 A CN 201610109217A CN 105700357 B CN105700357 B CN 105700357B
- Authority
- CN
- China
- Prior art keywords
- boiler combustion
- control
- pfc
- combustion system
- pid
- 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.)
- Active
Links
- 238000002485 combustion reaction Methods 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000004044 response Effects 0.000 claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims description 30
- 239000013598 vector Substances 0.000 claims description 29
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 15
- 239000001301 oxygen Substances 0.000 claims description 15
- 229910052760 oxygen Inorganic materials 0.000 claims description 15
- 239000007983 Tris buffer Substances 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 238000013178 mathematical model Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 abstract description 8
- 238000011217 control strategy Methods 0.000 abstract description 4
- 230000007812 deficiency Effects 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 27
- 230000008569 process Effects 0.000 description 11
- 230000007547 defect Effects 0.000 description 7
- 239000007789 gas Substances 0.000 description 7
- 230000001276 controlling effect Effects 0.000 description 6
- 239000000446 fuel Substances 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 230000001105 regulatory effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/048—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators using a predictor
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Steam Boilers And Waste-Gas Boilers (AREA)
- Regulation And Control Of Combustion (AREA)
Abstract
本发明公开了基于多变量PID‑PFC的锅炉燃烧***控制方法,本发明将多变量PID和预测函数控制相结合,得到一种可应用于多输入多输出***的新型控制方法,将该方法引入到锅炉燃烧控制***中来替代传统的PID控制器,提供一种新型的控制策略。本方法克服了传统PID控制器存在的跟踪性能差、超调量大、受干扰后速降大的缺点,同时本方法相对传统PFC控制,克服了其稳态误差大和多个基函数造成的计算量大的不足。本方法不仅能够在锅炉燃烧***稳态运行中提高***运行效率,也可以在进行效率优化的同时提高***的响应速度,使得锅炉燃烧***在整个运行过程中都能兼顾效率与响应性能。在实际应用上,本发明针对锅炉燃烧***提供了一种新型有效的控制策略。
Description
技术领域
本发明涉及锅炉燃烧***控制的技术领域,特别是基于多变量PID-PFC的锅炉燃烧***控制方法。
背景技术
当今社会科学技术水平越来越高,燃气锅炉本身具有众多的优越特性,而且其控制简单,方式相当灵活,运送与装配都方便易行,因此其具有相当明朗的应用前景,是必然的未来趋势。锅炉***可以看作是多输入多输出(MIMO)的复杂模型,文献[张振兴,丁宝.燃气热水锅炉智能控制***的研究[D].哈尔滨:哈尔滨工业大学,2006:14-21.Karnopp D C,Margolis D L,Rosenberg R C.System d-ynamics;modeling and simulation ofmechatronic systems[M].New York:John Wiley and Sons Inc,2000:46-52.]分析了其动态特性,结果表明锅炉***可以相对的分为三个分别独立的***模型:1)锅炉燃烧***,2)蒸汽发生***,3)蒸汽过热***。
在如今实际应用领域,后两个***模型的控制方法已经相对成熟,趋于完善,但是锅炉燃烧***的控制却不尽如人意,因为它是一非线性时变多变量***,且其强耦合、扰动剧烈并幅值变化相当大。文献[谷洋洋,李来春,张绍娟,等.基于智能PID控制的燃气锅炉燃烧控制***研究[J].热能动力工程,2015,30(3):413-416.]的PID控制方法是比较早期的理论方法,相对难实现有效控制需求。文献[Satoru Gotoa,Masatoshi Nakamura,ShiroMatsumura.Aut-omatic realization of human experience for controlling var-iable pressure boilers[J].Control Engineering Practice,2002,10:15-22.]进一步有针对性设计了混合式智能控制方法,将PID、专家和前馈三种控制结合起来,融合三者优势于一身,发挥各自优势,但是该方法仍存在计算过程复杂、耗时长等的缺点。
预测控制是近年来发展起来的一类新型的计算机控制算法。其适用于不易建立精确数字模型且动态过程复杂的工业生产,所以它一出现就受到国内外工程界的重视,并已在石油、化工、冶金、机械等工业部门的控制***中得到了成功的应用。电机***是个典型的多变量、强耦合、动态过程复杂的控制***,同样难于建立精确的数学模型,但同时又是一个快速***,传统的预测控制在线计算量大,实时性差,可能伴随不明的控制规律,不适用于感应电机的转速控制。在这种背景下,预测函数(PFC)控制方法应运而生,其基于预测控制的基本原理发展而来,其详细内容可参见文献[王树青,金晓明.先进控制技术应用实例[M].北京,化学工业出版社,2005.]。预测函数与预测控制方法的基本原理基本相同:模型预测、滚动优化、反馈校正。其与预测控制的最大区别是注重控制量的结构形式,认为控制量是一组预先选定的基函数的线性组合。在国外,PFC已经在工业机器人的快速高精度跟踪、军事领域的目标跟踪等快速***中得到了成功的运用。但目前尚未发现将PID和预测函数控制结合并推广到多输入、多输出***,并将方法应用于锅炉燃烧***控制的文献、报导。针对PID控制的不足,文献[Sheng Wu,Ridong Zhang,Renquan Lu,Furong Gao.Designof dynamic matrix control based PID for residual oil outlet temperature in acoke furnace[J].Chemometrics and Intell-igent Laboratory Systems,2014(134):110-117.]提出一种基于PID控制的DMC控制算法,其结构相对简单,易于控制,其较好弥补PID算法不足,但是仍有参数量较大,造成计算量大且计算耗时也较长的不足。文献[Havlena V,Findejs J.Application of model predictive control to advancedcombustion control[J].Control Engi-neering Practice,2005,13(6):671-680.]进一步提出预测控制算法解决锅炉燃烧***的蒸汽压力控制问题,而且较好应对了空燃比的问题,该文献中控制算法将自回归模型(ARX)应用其中。现实应用中找到较为理想空燃比是很难的,因此这种传统控制方法在燃料热设定值发生剧烈跳变时无法实现的。同时,文献[朱玉璧,程相利,陶新建,等.智能控制在锅炉燃烧优化中的应用[J].中国电机工程学报,2008,28(11):82-86.]针对锅炉燃烧***有效控制燃烧的过程,设计了一种非线性模型预测控制算法,优化了算法。
发明内容
本发明所要解决的技术问题是克服现有技术的不足而提供基于多变量PID-PFC的锅炉燃烧***控制方法,本发明将多变量PID和预测函数控制相结合,得到一种可应用于多输入多输出***的新型控制方法,将该方法引入到锅炉燃烧控制***中来替代传统的PID控制器;本方法不仅能够在锅炉燃烧***稳态运行中提高***运行效率,也可以在进行效率优化的同时提高***的响应速度,使得锅炉燃烧***在整个运行过程中都能兼顾效率与响应性能。
本发明为解决上述技术问题采用以下技术方案:
根据本发明提出的基于多变量PID-PFC的锅炉燃烧***控制方法,包括以下步骤:
步骤1、将基于锅炉燃烧***的数学模型转化为状态空间方程,得出系数矩阵Am、Bm、Cm;
步骤2、根据下式计算控制量向量u(n):
u(n)=(Rp+Ri+Rd)Tfn(0)
其中:
Rp=(KiGTQG+fTRf)-1(Kp+Ki+Kd)DTQG
Ri=(KiGTQG+fTRf)-1(-Kp-2Kd)(q-1DTQG)
Rd=(KiGTQG+fTRf)-1Kd(q-2DTQG)
f=[fn1(0),fn2(0),···,fnJ(0)]T
Gnj(i)=CmAm i-1Bmfnj(0)+CmAm i-2Bmfnj(1)+···+CmBmfnj(i-1)
其中,上标T是矩阵转置,fn(0)表示基函数在时间为0时的值,Kp、Ki、Kd分别为广义比例项系数、积分项系数和微分项系数;fnj(i)为选定的基函数,i是整数且1≤i≤ns,f为基函数的值构成的矩阵,j为1到J之间整数,J是基函数的阶数;Q和R分别表示误差加权矩阵和控制加权矩阵;q-1和q-2为延时算子;ns为优化时域拟合点的个数,hi为第i个拟合点上的数值;yP(n)为当前时刻该锅炉燃烧***的三个输出量构成的向量,三个输出量分别为:主蒸汽压力、氧含量和炉膛负压;C(n)为工程中设定的该锅炉燃烧***的三个输出量构成的向量,三个输出量分别为主蒸汽压力、氧含量和炉膛负压;To是采样时间,Tr是参考轨迹的期望响应时间,I是单位矩阵,xm(n)是第n时刻模型状态值;
步骤3、将控制量u(n)以可执行文件的形式加载到DSP的RAM中,DSP的CAP口捕获单元读取位置信号,控制量经过多变量PID-PFC控制器后得到锅炉燃烧***的三个输出量的实际值,这三个输出量分别为:主蒸汽压力、氧含量和炉膛负压;将这三个输出量的实际值与预设的这三个输出量的参考值进行比较得到偏差,通过其偏差反馈调整控制量,从而控制锅炉燃烧***的运行。
作为本发明所述的基于多变量PID-PFC的锅炉燃烧***控制方法进一步优化方案,所述状态空间方程为:
其中,Ym(k)为k时刻模型预测输出向量,Xm(k)为k时刻模型状态值向量,U(k-1)为(k-1)时刻控制输入向量,Am、Bm、Cm为矩阵方程系数矩阵。
作为本发明所述的基于多变量PID-PFC的锅炉燃烧***控制方法进一步优化方案,步骤2中基函数fnj(i)为单位阶跃函数。
作为本发明所述的基于多变量PID-PFC的锅炉燃烧***控制方法进一步优化方案,单位阶跃函数的阶数J为1。
作为本发明所述的基于多变量PID-PFC的锅炉燃烧***控制方法进一步优化方案,所述ns为5。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
(1)在燃用混合煤气锅炉正常运行下,多变量PID预测函数控制相比传统的预测函数控制精度高,跟踪速度较快,稳态误差小,抗干扰能力强,既能保证锅炉燃烧***具有较好的稳定性与动态性能,又能提高锅炉动态时的运行效率;
(2)与包含多个基函数的预测函数控制相比,多变量PID预测函数算法计算量小,可根据实际情况调节Kp,Ki,Kd三个参数矩阵,控制灵活方便,无需重新修改控制程序,同时解决了预测函数控制模型不匹配导致的控制效果变差的问题;
(3)本发明有着跟踪速度较快,静差较小,抗干扰能力强,兼顾了动态响应速度和效率优化的优点。
附图说明
图1是锅炉燃烧***原理示意图。
图2是锅炉燃烧***输入输出关系图。
图3是锅炉燃烧***整体示意图。
图4是多变量PID-PFC输出量曲线。
图5是多变量PID-PFC控制量曲线。
图6是多变量PID-PFC、多变量PID与多变量PFC输出波形对比。
图7是多变量PID-PFC、多变量PIDMAC输出波形对比。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明的目的在于将PID和预测函数控制相结合,得到一种可应用于多输入多输出***的新型控制方法,将该方法引入到锅炉燃烧控制中来替代传统的预测函数控制,提供一种新型的控制策略。
1、选取基函数和参考轨迹
预测函数控制把控制输入结构看作影响***性能的关键。而在预测函数控制中在输入信号频谱有限的情况,控制输入仅属于一组与参考轨迹和对象性质有关的特定基函数族,基函数的选取的重要性可想而知。特别的,对于线性***的输出将是上述基函数作用于对象模型响应的加权组合。控制输入被表示为一系列已知基函数{fj}的线性组合,即
在上式中:U(k+i)为在k+i时刻的控制量向量,μj(k)为基函数加权系数向量,fj(i)为基函数在(k+i)Ts时的取值,J为基函数的阶数,P为预测步长。
在PFC(预测函数)中,为了使***的输出能够平缓地逐渐达到设定值,避免出现超调,根据预测输出值和过程输出值,我们可以规定一条渐进趋向于未来设定值的曲线,称为参考轨迹。其选定完全取决于设计者对***的要求。常见的参考轨迹如下:
Yr(k+i)=c(k+i)-αi[c(k)-Yp(k)] (2)
上式中:Yr(k+i)为(k+i)时刻的参考轨迹向量,YP(k)为k时刻的过程实际输出值向量,
c(k)为k时刻的设定值组成的向量,cn(k)为k时刻第n个变量的设定值,c(k)=[c1(k) c2(k) … cN(k)]T,n=1,2,…,N,αi为第i时刻的参考轨迹衰减因子,表征了参考轨迹趋于设定值的快慢程度,一般取其中Ts是采样时间,Tr是参考轨迹的期望响应时间,n=1,2,…,N。
2、建立基于混合煤气锅炉控制***空间模型
本发明控制对象为典型的混合煤气锅炉,根据现场采集的大量动态数据,按2s为一周期,经多变量预测误差模型辨识方法并可观规范型实施,可得被控***的空间模型:
式中,Ym(k)---k时刻模型预测输出向量;Xm(k)---k时刻模型状态值向量;
U(k-1)---(k-1)时刻控制输入向量;Am、Bm、Cm---矩阵方程系数矩阵。
对应输出量的三个控制目标:1.主蒸汽压力:控制在额定值1~3MPa左右;2.氧含量:控制在3%~4%左右的经济范围之内;3.炉膛负压:控制在-40~-20Pa左右的安全范围之间。
3、计算预测模型的模型输出
对于(k+i)时刻的模型状态值Xm(k+i),由上式(3)递推得到
由此可知,(k+i)时刻的模型预测输出为
其中,
4、计算补偿后的模型预测输出
在实际工业过程中,由于模型失配、噪声等原因,模型输出与过程输出之间存在一定的误差,即:
对于未来(k+i)时刻误差的预测,在控制***中可认为:
其中:为k时刻的误差向量,en(k)为第n个模型输出与过程输出之间的误差,n=1,2,…,N;
Yp(k)为k时刻的过程实际输出向量;
Ym(k)为k时刻的模型预测输出向量。
则未来P时刻预测模型被修正为
实际过程预测输出表达式为:
5、基于二次型PID目标函数求解出控制量
为了使控制***具有更好的控制品质,在多变量控制***中把PID控制和PFC控制结合起来,采用加入比例、积分,微分的新的目标函数,使推导的控制器具有广义上的比例、积分的结构特性。利用PID算法对PFC算法的目标函数进行改进,推导出的新型多变量PID预测函数算法不仅具有PID与PFC算法的优点,还能克服它们的缺点。
上式里,Kp,Ki,Kd分别为比例系数矩阵、积分系数矩阵和微分系数矩阵,Q和R分别为误差加权与控制量加权因子,且为正定矩阵,为预测误差,为预测误差的增量,为预测误差增量平方。且有,
U(k+i)=[u1(k+i),…,um(k+i)]T,i=0,1,2…,P-1
由于模型存在失配及噪音的原因,对象实际输出与模型预测输出不可避免的会有一定的误差:
上式里的未来p时刻预测模型被修改成:
令k+i时刻误差可表示为:
式中,为设定值矩阵,
定义
式中,D(k)=[D1 T(k),…,DP T(k)]T
由递推原理可得:
由式(8),多变量PID-PFC目标函数得
在此,引入后移算子q-1,则
ΔD(k)=(1-q-1)D(k),
Δ2D(k)=(1-q-1)2ΔD(k),
将式(10)代入后移算子,得
令W=Kp(1-q-1)2+Ki+Kd(1-q-1)4
则最优控制量可简化为:
令可得:
u(n)=(Rp+Ri+Rd)Tfn(0)
其中:
Rp=(KiGTQG+fTRf)-1(Kp+Ki+Kd)DTQG
Ri=(KiGTQG+fTRf)-1(-Kp-2Kd)(q-1DTQG)
Rd=(KiGTQG+fTRf)-1Kd(q-2DTQG)
f=[fn1(0),fn2(0),···,fnJ(0)]T
Gnj(i)=CmAm i-1Bmfnj(0)+CmAm i-2Bmfnj(1)+···+CmBmfnj(i-1)
其中,u(n)是***第n个时刻的控制量输出;Kp、Ki、Kd分别为广义比例项系数、积分项系数和微分项系数;fnj(i)为选定的基函数,f为基函数的值构成的矩阵,下标J表示基函数的阶数,j为1到J之间整数;Q和R分别表示误差加权矩阵和控制加权矩阵;q-1和q-2为延时算子;ns为优化时域拟合点的个数,hi为第i个拟合点上的数值;yP(n)为当前时刻该锅炉燃烧***的三个输出量(主蒸汽压力、氧含量和炉膛负压)构成的向量;C(n)为工程中设定的该锅炉燃烧***的三个输出量(主蒸汽压力、氧含量和炉膛负压)构成的向量;To是采样时间,Tr是参考轨迹的期望响应时间;
可得到最终的控制量:u(n)=(Rp+Ri+Rd)Tfn(0)
具体按照以下方法:
步骤1、将基于锅炉燃烧***的数学模型转化为状态空间方程,得出系数矩阵Am、Bm、Cm。
步骤2、根据下式计算控制量向量u(n):
u(n)=(Rp+Ri+Rd)Tfn(0)
其中:
Rp=(KiGTQG+fTRf)-1(Kp+Ki+Kd)DTQG
Ri=(KiGTQG+fTRf)-1(-Kp-2Kd)(q-1DTQG)
Rd=(KiGTQG+fTRf)-1Kd(q-2DTQG)
f=[fn1(0),fn2(0),···,fnJ(0)]T
Gnj(i)=CmAm i-1Bmfnj(0)+CmAm i-2Bmfnj(1)+···+CmBmfnj(i-1)
其中,u(n)是***第n个时刻的控制量输出;Kp、Ki、Kd分别为广义比例项系数、积分项系数和微分项系数;fnj(i)为选定的基函数,f为基函数的值构成的矩阵,下标J表示基函数的阶数,j为1到J之间整数;Q和R分别表示误差加权矩阵和控制加权矩阵;q-1和q-2为延时算子;ns为优化时域拟合点的个数,hi为第i个拟合点上的数值;yP(n)为当前时刻该锅炉燃烧***的三个输出量(主蒸汽压力、氧含量和炉膛负压)构成的向量;C(n)为工程中设定的该锅炉燃烧***的三个输出量(主蒸汽压力、氧含量和炉膛负压)构成的向量;To是采样时间,Tr是参考轨迹的期望响应时间;
步骤3、所诉锅炉燃烧控制为多变量PID预测函数控制,其控制参数为步骤2计算所得的控制量向量u(n);
将控制量u(n)以可执行文件的形式加载到DSP的RAM中,DSP的CAP口捕获单元读取位置信号,控制量经过多变量PID-PFC控制器后得到锅炉燃烧***的三个输出量(主蒸汽压力、氧含量和炉膛负压)。将三个输出量的参考值与实际输出量反馈值进行比较得到偏差,通过其偏差反馈,从而控制锅炉燃烧***的运行。
优选地,步骤2中所述基函数fnj(i)为单位阶跃函数,其阶数J的取值为1。
优选地,步骤2中所述基于PID预测函数的制方法,其特征在于,步骤2中所述优化时域拟合点的个数ns的取值范围为5。
为了验证本发明方法的效果,进行了下述实验:
经反复比较,多变量PID-PFC的调节参数以Kp=14,Ki=1000,Kd=6,P=5,M=2为宜。
图1是本发明基于PID-PFC的控制方法的原理示意图,三个控制量(燃料量、引风量、送风量)经过本设计出的控制策略,得到预期的输出量。
从控制的角度来看,锅炉是燃烧过程一个复杂的调节对象,其中有着许多调节参数和被调节的参数,同时存在着难以精确预料的扰动参数。这些参数间相互作用如图2所示。
理想锅炉燃烧控制***是多回路的控制***,但其控制相对复杂。因此,目前实际解决锅炉燃烧控制的方法是将锅炉看成几个相对独立的调节对象,将锅炉燃烧***视为三输入三输出***,输入量为送风量、燃料量和引风量;输出量为主蒸汽量、含氧量和炉膛压力。锅炉燃烧***结构如图3所示。
从图4分析可知:氧含量、主蒸汽压力和炉膛负压的上升时间分别为39s、34s和20s;三者可以看出几乎不存在超调量,其稳态误差也近于零;当主蒸汽压力在400s由1跳变到2时,炉膛负压和氧含量的相应变化很小,其控制效果较为理想。
针对图5的控制量波形图,多变量PID-PFC算法对锅炉燃烧***的三个控制量在开始时能做出快速反应,且时间都在100s之内。同时在第400s出现跳变时,三个控制量也可以做出迅速响应,时间不超过55s。从图5可以看出除去开始和跳变时间,其他情况下,三个控制量能较好的平缓运行,控制品质良好,***运行状态理想。
本发明在对多变量PID-PFC进行仿真,同时和多变量PID-PFC与多变量PID与多变量PFC控制效果进行整体的对比,分析其控制性能。
经多次对比,多变量PID-PFC参数选取如下Kp=14,Ki=1000,Kd=6,P=5,M=2。
经过多次对比,多变量PFC与PID参数选取如下:P=6,M=2。
图6中,实线、虚线和点实线分别代表PIDPFC的三个被控对象的输出响应、PFC的三个被控制对象的输出响应和PID的三个被控制对象的输出响应。由上述数据分析表明:针对三个被控量的输出响应来看,实线所代表的PIDPFC控制在整体上优越于PID和PFC控制,其近乎没有超调量,稳态误差也趋于零,控制性能良好。在个别方面,点实线代表的PID控制确有其优势,如PID控制的三个被控量曲线的上升时间明显快于其他两种算法,但其超调量过于大,对现实应用有明显的缺陷,因此整体看来,本算法的有效性高于其他两种算法。
本多变量PID-PFC算法将和文献[郭伟,王汉杰,夏友亮,周丽.基于状态空间方程的多变量PID-MAC在锅炉燃烧控制***中的应用[J].热力发电,2014,43(9):48-53.]的算法进行比较仿真,两者的参数见下表1和表2。
表1本算法参数
表2文献算法参数
图7可见,实线和虚线分别代表本发明PIDPFC的三个被控制对象输出响应和文献[14]PIDMAC的三个被控对象输出响应。在第400s时,针对主蒸汽压力的跟踪值突然跳变,由1跳变到2;第500s时,炉膛负压由-3跳变到-2的跟踪值;通过比较分析,明显看出本发明算法与文献算法在响应速度、超调量和调节时间这三个方面相比,本算法都具有优越性,其响应快速,并且没有超调量,动态性能和稳态性能更是好于文献中的仿真结果,体现其优越性。
Claims (5)
1.基于多变量PID-PFC的锅炉燃烧***控制方法,其特征在于,包括以下步骤:
步骤1、将基于锅炉燃烧***的数学模型转化为状态空间方程,得出系数矩阵Am、Bm、Cm;
步骤2、根据下式计算控制量向量u(n):
u(n)=(Rp+Ri+Rd)Tfn(0)
其中:
Rp=(KiGTQG+fTRf)-1(Kp+Ki+Kd)DTQG
Ri=(KiGTQG+fTRf)-1(-Kp-2Kd)(q-1DTQG)
Rd=(KiGTQG+fTRf)-1Kd(q-2DTQG)
f=[fn1(0),fn2(0),…,fnJ(0)]T
Gnj(i)=CmAm i-1Bmfnj(0)+CmAm i-2Bmfnj(1)+…+CmBmfnj(i-1)
其中,上标T是矩阵转置,fn(0)表示基函数在时间为0时的值,Kp、Ki、Kd分别为广义比例项系数、积分项系数和微分项系数;fnj(i)为选定的基函数,i是整数且1≤i≤ns,f为基函数的值构成的矩阵,j为1到J之间整数,J是基函数的阶数;Q和R分别表示误差加权矩阵和控制加权矩阵;q-1和q-2为延时算子;ns为优化时域拟合点的个数,hi为第i个拟合点上的数值;yP(n)为当前时刻该锅炉燃烧***的三个输出量构成的向量,三个输出量分别为:主蒸汽压力、氧含量和炉膛负压;C(n)为工程中设定的该锅炉燃烧***的三个输出量构成的向量,三个输出量分别为主蒸汽压力、氧含量和炉膛负压;To是采样时间,Tr是参考轨迹的期望响应时间,I是单位矩阵,xm(n)是第n时刻模型状态值;
步骤3、将控制量u(n)以可执行文件的形式加载到DSP的RAM中,DSP的CAP口捕获单元读取位置信号,控制量经过多变量PID-PFC控制器后得到锅炉燃烧***的三个输出量的实际值,这三个输出量分别为:主蒸汽压力、氧含量和炉膛负压;将这三个输出量的实际值与预设的这三个输出量的参考值进行比较得到偏差,通过其偏差反馈调整控制量,从而控制锅炉燃烧***的运行。
2.根据权利要求1所述的基于多变量PID-PFC的锅炉燃烧***控制方法,其特征在于,所述状态空间方程为:
其中,Ym(k)为k时刻模型预测输出向量,Xm(k)为k时刻模型状态值向量,U(k-1)为(k-1)时刻控制输入向量,Am、Bm、Cm为矩阵方程系数矩阵。
3.根据权利要求1所述的基于多变量PID-PFC的锅炉燃烧***控制方法,其特征在于,步骤2中基函数fnj(i)为单位阶跃函数。
4.根据权利要求3所述的基于多变量PID-PFC的锅炉燃烧***控制方法,其特征在于,单位阶跃函数的阶数J为1。
5.根据权利要求1所述的基于多变量PID-PFC的锅炉燃烧***控制方法,其特征在于,所述ns为5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610109217.5A CN105700357B (zh) | 2016-02-25 | 2016-02-25 | 基于多变量pid-pfc的锅炉燃烧***控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610109217.5A CN105700357B (zh) | 2016-02-25 | 2016-02-25 | 基于多变量pid-pfc的锅炉燃烧***控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105700357A CN105700357A (zh) | 2016-06-22 |
CN105700357B true CN105700357B (zh) | 2018-06-26 |
Family
ID=56222480
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610109217.5A Active CN105700357B (zh) | 2016-02-25 | 2016-02-25 | 基于多变量pid-pfc的锅炉燃烧***控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105700357B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107065541A (zh) * | 2017-03-22 | 2017-08-18 | 杭州电子科技大学 | 一种焦化炉炉膛压力***模糊网络优化pid‑pfc控制方法 |
CN107390515A (zh) * | 2017-08-01 | 2017-11-24 | 东华大学 | 预测pi算法和多变量解耦相结合的锅炉燃烧控制*** |
CN109062030A (zh) * | 2018-07-27 | 2018-12-21 | 南京信息工程大学 | 基于拉盖尔函数模型的火电单元机组负荷预测pid控制方法 |
CN111429010B (zh) * | 2020-03-27 | 2024-01-09 | 上海圆曦电力科技有限公司 | 一种基于内模控制结构的火力发电机组协调控制*** |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101211162B (zh) * | 2007-12-24 | 2012-05-16 | 庞国仲 | 用于燃煤机组和工业锅炉燃烧过程多变量时滞控制方法 |
CN101498457B (zh) * | 2009-03-02 | 2010-08-25 | 杭州电子科技大学 | 一种锅炉燃烧优化的方法 |
CN202512382U (zh) * | 2012-04-20 | 2012-10-31 | 青岛科技大学 | 链条锅炉燃烧滚动自寻优-pid复合控制*** |
JP6102475B2 (ja) * | 2013-05-01 | 2017-03-29 | 三浦工業株式会社 | ボイラシステム |
JP6164064B2 (ja) * | 2013-11-27 | 2017-07-19 | 三浦工業株式会社 | ボイラシステム |
CN103676651B (zh) * | 2013-12-02 | 2018-03-06 | 国家电网公司 | 基于状态观测模型的锅炉汽温预测控制方法 |
-
2016
- 2016-02-25 CN CN201610109217.5A patent/CN105700357B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105700357A (zh) | 2016-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105700357B (zh) | 基于多变量pid-pfc的锅炉燃烧***控制方法 | |
CN107479389B (zh) | 一种火电机组过热汽温预测模糊自适应pid控制方法 | |
CN108227488B (zh) | 基于滑模预测控制的超超临界火电机组协调控制方法 | |
CN104699064B (zh) | 一种火电机组协调控制***的煤质自校正预测控制方法 | |
CN105425612B (zh) | 一种水轮机调节***控制参数的优选方法 | |
CN102494336A (zh) | 一种循环流化床锅炉燃烧过程多变量控制方法 | |
CN104122797A (zh) | 一种新型火电机组负荷多变量预测控制方法 | |
CN107270283B (zh) | 一种基于循环流化床机组的多变量约束预测控制方法 | |
CN106707756B (zh) | 融合扩张观测器的超临界火电机组机炉协调控制方法 | |
CN110361974B (zh) | 基于bp-foa混合算法的水轮机调速***优化方法 | |
Lei et al. | Online optimization of fuzzy controller for coke-oven combustion process based on dynamic just-in-time learning | |
CN113835342B (zh) | 一种过热汽温***的抗扰预测控制方法 | |
CN100547504C (zh) | 基于模糊目标与模糊约束的工业过程预测控制方法 | |
CN106873377B (zh) | 不确定扰动下的焦炉加热燃烧过程火道温度自动调节方法 | |
CN107102550B (zh) | 一种超超临界火电机组控制分离器温度的预测控制方法 | |
CN106127351A (zh) | 一种高炉热风炉初始烧炉煤气流量优化*** | |
Tian et al. | The tuning principle of adaptive fuzzy fractional-order PID controller parameters | |
Duan | The design of predictive fuzzy-PID controller in temperature control system of electrical heating furnace | |
CN104850008B (zh) | 一种层燃锅炉的多变量控制与实时优化的方法 | |
CN114740713A (zh) | 一种湿法烟气脱硫过程的多目标优化控制方法 | |
CN113885326A (zh) | 考虑前馈的基于设定值有效激励的闭环***辨识方法 | |
CN113282043A (zh) | 基于多变量状态空间模型的超超临界机组协调控制方法 | |
CN107065538B (zh) | 汽包锅炉-汽轮机组的模糊跟踪控制方法 | |
Huang et al. | Simulation and analysis of main steam control system based on heat transfer calculation | |
CN110970936A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20190415 Address after: 257000 Dongying Port Economic Development Zone, Dongying City, Shandong Province Patentee after: SHANDONG SHENCHI PETROCHEMICAL Co.,Ltd. Address before: 210000 69 Olympic Sports street, Jianye District, Nanjing, Jiangsu. Patentee before: Nanjing University of Information Science and Technology |
|
PP01 | Preservation of patent right | ||
PP01 | Preservation of patent right |
Effective date of registration: 20240125 Granted publication date: 20180626 |