CN103576553B - 一种燃煤锅炉蒸汽温度的分数阶自整定控制方法 - Google Patents

一种燃煤锅炉蒸汽温度的分数阶自整定控制方法 Download PDF

Info

Publication number
CN103576553B
CN103576553B CN201310541805.2A CN201310541805A CN103576553B CN 103576553 B CN103576553 B CN 103576553B CN 201310541805 A CN201310541805 A CN 201310541805A CN 103576553 B CN103576553 B CN 103576553B
Authority
CN
China
Prior art keywords
fractional order
temperature
steam
vapor
coal
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.)
Expired - Fee Related
Application number
CN201310541805.2A
Other languages
English (en)
Other versions
CN103576553A (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.)
North China Electric Power University
Original Assignee
North China Electric Power 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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN201310541805.2A priority Critical patent/CN103576553B/zh
Publication of CN103576553A publication Critical patent/CN103576553A/zh
Application granted granted Critical
Publication of CN103576553B publication Critical patent/CN103576553B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种燃煤锅炉蒸汽温度的分数阶自整定控制方法。它属于自动控制领域,特别涉及一种基于分数阶模型的自整定控制方法。包括以下步骤:以采样间隔T对当前时刻t的蒸汽温度信号和调节量信号采样,并记录到历史数据库;计算控制系统的性能指标;确定是否更新控制器参数;从历史数据库中获取建模数据序列;采用成熟的群智能算法优化得到被控系统的分数阶数学模型;采用成熟的群智能算法优化控制器参数;最后计算控制器输出,施加到被控对象。本发明适合燃煤锅炉蒸汽温度系统的数字控制,具有步骤明确、概念清晰、鲁棒性强的特点,无需对控制对象做专门的扰动试验,也不需要在控制信号中加入伪随机序列,能够明显提高蒸汽温度的调节品质。

Description

一种燃煤锅炉蒸汽温度的分数阶自整定控制方法
技术领域
本发明涉及一种燃煤锅炉蒸汽温度的控制方法,具体是一种用于燃煤锅炉蒸汽温度调节的分数阶自整定控制方法,属于自动控制技术领域。
背景技术
虽然风力发电、太阳能发电、生物质发电,以及核能发电和水力发电等新能源发电技术近年来取得了显著的技术进步,并得到广泛的应用,但是燃煤锅炉发电机组仍然是我国相当长一段时间内的主力发电方式。燃煤锅炉发电厂的蒸汽温度(包括过热蒸汽温度和再热蒸汽温度)需要维持在相当精度的范围,温度过高会大大影响机组的使用寿命,温度过低又会大大降低锅炉效率。同时即使温度被控制在设计值附近但频繁波动,也同样会影响机组的使用寿命以及无谓的消耗能源。因此严格控制蒸汽温度在设计值附近是延长机组使用寿命、提高锅炉效率的重要手段之一。
经对现有技术的公开文献检索发现,A.Sanchez-Lopez,G.Arroyo-Figueroa,A.Villavicencio-Ramirez.Advanced control algorithms for steam temperatureregulation of thermal power plants.Electrical Power and Energy Systems,2004,26:779-785(热力发电厂蒸汽温度调节的先进控制算法,国际期刊:电力与能源系统,2004,26:779-785),虽然使用了先进的动态矩阵控制和模糊控制算法可以在一定程度上改善系统的调节品质,但是对于蒸汽温度系统,其特性随机组发电负荷的变化而变化,固定参数的控制器很难适应所有运行工况。文献(王万召,赵兴涛,宋艳萍.模糊RBF自整定PID控制器在过热汽温控制中应用.电力自动化设备,2007,27(11):48-50)以及文献(李致远.基于单神经元自适应PID的过热汽温控制系统研究.工业控制计算机,2010,23(8):65-66),采用了控制器参数的自整定(自适应)技术,在一定程度上可以使控制器适应蒸汽温度系统特性随机组发电负荷以及各种扰动的变化,但是由于蒸汽温度系统本质上具有严重的分布参数特性,在控制系统的属性描述上表现为分数阶微分方程和分数阶传递函数,采用常规的整数阶控制则不能得到很好的调节品质。
因此,本发明基于分数阶微积分理论,采用分数阶模型和分数阶控制器以及现代智能优化技术,提供一种新的燃煤锅炉蒸汽温度的自整定控制方法。
发明内容
本发明的目的在于提供一种基于分数阶模型、用于燃煤锅炉蒸汽温度调节的自整定控制方法,所发明的自整定控制方法能适用于各种燃煤锅炉蒸汽温度的自动调节。
为了实现上述目的,本发明的技术方案是根据以下步骤实现的:
步骤(1).以采样间隔T对当前时刻t的蒸汽温度信号y(t)和调节量(一般为减温水)信号u(t)采样,并记录到历史数据库;
步骤(2).计算反映控制系统调节效果的均方误差性能指标J(t);
步骤(3).确定是否更新控制器参数,如果是则进入步骤(4),否则转步骤(7);
步骤(4).从历史数据库中获取建模数据列;
步骤(5).采用成熟的群智能优化算法优化得到被控系统的分数阶数学模型;步骤(6).采用成熟的群智能优化算法优化控制器参数;
步骤(7).计算控制器输出,施加到被控对象,返回步骤(1)继续下一次采样;
其中所述蒸汽温度信号为过热蒸汽温度或再热蒸汽温度,所述调节量为减温水。
作为本发明的进一步改进,所述步骤1中采样间隔T取值范围5~20秒。
作为本发明的进一步改进,所述步骤2中反映控制系统调节效果的均方误差性能指标J(t)是在过去一段时间内被控变量蒸汽温度信号y(t)与期望值R(t)的均方差,即
J ( t ) = Σ i = 0 N - 1 [ R ( t - i ) - y ( t - i ) ] 2 / N
式中,N是计算性能指标所用的数据长度,取50~200;
y(t-i)和R(t-i)分别表示蒸汽温度信号y(t)与期望值R(t)在t-i时刻的值i=0,1,2,...,N-1。
作为本发明进一步的改进,所述步骤3中是否更新控制器参数的依据是均方误差性能指标J(t)是否达到给定的某个阈值J0,如果J(t)≥J0,则需要更新控制器参数。
作为本发明进一步的改进,所述步骤4中建模数据列取最近时间段内的200~600组数据。
作为本发明进一步的改进,所述步骤(5)和步骤(6)中群智能算法是混合微分进化与量子粒子群算法。
作为本发明进一步的改进,所述步骤(5)中包含内回路的广义被控系统的分数阶数学模型为如下形式
K ( 1 + T 2 S ) 2 ( 1 + T 1 S y ) ,
其中,s是拉氏变换的算子;K是被控系统模型的增益;T2是被控系统模型的整数阶部分的时间常数;T1是被控系统模型的分数阶部分的时间常数;n是被控系统模型的整数阶阶次;γ是被控系统模型的分数阶阶次。
作为本发明进一步的改进,所述步骤(6)中控制器为如下形式的分数阶PIIDP控制器
G c ( s ) = K p + K 1 s i + K D s P , λ , μ > 0 ,
其中,KP、KI和KD分别是分数阶PIIDP控制器的比例系数、积分系数和微分系数;λ、μ分别是分数阶积分阶次和微分阶次;s是拉氏变换的算子。参数λ、μ的不同取值可以获得比常规PID控制更宽的调节范围,常规PID控制只是平面内的一个点,而分数阶PIIDP控制器则占据一个矩形区域的范围。
本发明的有益效果是:控制器参数的整定无需对控制对象做专门的扰动试验,也不需要在控制信号中加入伪随机序列,而是直接根据性能指标判断获得模型辨识所需的建模数据列。同时由于使用了分数阶传递函数和分数阶PIIDP控制器,控制系统设计的自由度得到提高,具有比常规PID控制更强的鲁棒性,从而可以提高蒸汽温度的调节品质。
附图说明
图1为蒸汽温度的分数阶传递函数辨识方框图;
图2为蒸汽温度的分数阶PIIDP控制方框图;
图3为分数阶PIIDP控制器的参数λ、μ取值范围图;
图4为某过热蒸汽温度系统的分数阶与整数阶控制效果图;
图5为某过热蒸汽温度系统在对象时间常数增加20%后的分数阶与整数阶控制效果图;
图6为某过热蒸汽温度系统在对象增益减小20%后的分数阶与整数阶控制效果图。
具体实施方式
以下结合附图及并用某电厂过热蒸汽温度控制的实施例来说明本发明的实施方法和步骤,但不用来限制本发明的范围。
本发明的模型描述和分数阶控制器分别如附图1和附图2所示,该模型的明显优点是:对于具有本质分数阶特性的锅炉蒸汽温度系统,模型的精度明显高于使用传统的基于整数阶微积分的整数阶模型。本发明PIIDP控制器的明显优点是:对于具有本质分数阶特性的锅炉蒸汽温度系统,采用分数阶控制器的控制效果明显优于使用传统的基于整数阶微积分的整数阶控制器。如附图3所示,分数阶PIIDP控制器的参数λ、μ取值范围明显大于使用传统的基于整数阶微积分的整数阶PID控制器,因而对控制器的参数调整自由度更大,从而能够取得更好的控制效果。常规PID控制只是平面内的一个点,而分数阶PIIDP控制器则占据一个矩形区域的范围。
下面结合本发明的具体实施方式对一种燃煤锅炉蒸汽温度的分数阶自整定控制方法进行说明。
具体步骤如下:
1.以采样间隔10秒对当前时刻t的蒸汽温度(过热蒸汽温度或再热蒸汽温度)信号y(t)和调节量(为减温水)信号u(t)采样,并记录到历史数据库;
2.计算反映控制系统调节效果的均方误差性能指标J(t),
J ( t ) = Σ i = 1 N - 1 [ R ( t - i ) - y ( t - i ) ] 2 / N ,
式中,N是计算性能指标所用的数据个数,N=100,R是过热蒸汽温度控制的目标值,R=540℃;
3.确定是否更新控制器参数。给定阈值J0=0.5,如果J(t)<J0,则表明当前系统控制效果良好,不需要更新控制器,从而转步骤7按照已有的控制器继续执行;如果J(t)≥J0,则表明当前系统控制效果不够好,需要更新控制器参数。本实施例计算结果为J(t)=1.126,从而需要继续步骤4,以便更新控制器参数;
4.从历史数据库中获取建模数据列400组,即获取蒸汽温度数据[y(t),y(t-1),……,y(t-399)]以及调节量数据[u(t),u(t-1),……,u(t-399)];
5.采用群智能算法优化得到被控系统的分数阶数学模型。采用群智能算法优化得到减温水量变化对主蒸汽温度的分数阶数学模型为如下形式(式中分子的负号表示减温水量增加时主蒸汽温度降低)
- 1.237 ( 1 + 222.5 s ) 2 ( 1 + 865.6 s 1.26 ) ;
6.采用群智能算法优化控制器参数。采用群智能算法优化得到分数阶PIIDP控制器为
G c ( s ) = 0.263995 + 0.001305 s 1.11 + 0.885212 s 0.18 ;
7.计算控制器输出,施加到被控对象,返回1继续下一次采样。
采用本发明方法将蒸汽温度从540℃调节至543℃的控制效果如图4,图中还对比性的给出了常规整数阶控制效果图。为说明本发明可以获得更好的蒸汽温度控制品质,图5和图6分别给出了对象时间常数增加20%以及对象增益减小20%后的鲁棒性验证控制结果。值得说明的是:(1)其中用于对比的常规整数阶PID控制器参数也是优化得到的最好结果,其参数为KP=0.5035,KI=0.0013,KD=0.7916;(2)所有的控制均是在考虑执行机构的物理限制条件下进行的,这一点对于工业应用具有重要的实际意义。

Claims (7)

1.一种燃煤锅炉蒸汽温度的分数阶自整定控制方法,其特征在于,包括以下步骤:
步骤1:以采样间隔T对当前时刻t的蒸汽温度信号y(t)和调节量信号u(t)采样,并记录到历史数据库;
步骤2:计算反映控制系统调节效果的均方误差性能指标J(t);
步骤3:确定是否更新控制器参数,如果是则进入步骤4,否则转步骤7;
步骤4:从历史数据库中获取建模数据列;
步骤5:采用成熟的群智能算法优化得到被控系统的分数阶数学模型;
步骤6:采用成熟的群智能算法优化控制器参数;
步骤7:计算控制器输出,施加到被控对象,返回步骤1继续下一次采样。
2.根据权利要求1所述的一种燃煤锅炉蒸汽温度的分数阶自整定控制方法,其特征在于,所述步骤1中采样间隔T取值范围5~20秒。
3.根据权利要求1所述的一种燃煤锅炉蒸汽温度的分数阶自整定控制方法,其特征在于,所述步骤2中反映控制系统调节效果的均方误差性能指标J(t)是在过去一段时间内被控变量蒸汽温度信号y(t)与期望值R(t)的均方差,即
J ( t ) = Σ i = 0 N - 1 [ R ( t - i ) - y ( t - i ) ] 2 / N
式中,N是计算性能指标所用的数据长度,取50~200;
y(t-i)和R(t-i)分别表示蒸汽温度信号y(t)与期望值R(t)在t-i时刻的值i=0,1,2,...,N-1;所述步骤3中是否更新控制器参数的依据是均方误差性能指标J(t)是否达到给定的某个阈值J0,如果J(t)≥J0,则需要更新控制器参数。
4.根据权利要求1所述的一种燃煤锅炉蒸汽温度的分数阶自整定控制方法,其特征在于,所述步骤4中建模数据列取最近时间段内的200~600组数据。
5.根据权利要求1所述的一种燃煤锅炉蒸汽温度的分数阶自整定控制方法,其特征在于,所述步骤5和步骤6中群智能算法是混合微分进化与量子粒子群算法。
6.根据权利要求1所述的一种燃煤锅炉蒸汽温度的分数阶自整定控制方法,其特征在于,所述步骤5中被控系统的分数阶数学模型是基于分数阶微积分理论得到的如下形式的分数阶传递函数
K ( 1 + T 2 S ) n ( 1 + T 1 S y )
其中,s是拉氏变换的算子;K是被控系统模型的增益;T2是被控系统模型的整数阶部分的时间常数;T1是被控系统模型的分数阶部分的时间常数;n是被控系统模型的整数阶阶次;y是被控系统模型的分数阶阶次。
7.根据权利要求1所述的一种燃煤锅炉蒸汽温度的分数阶自整定控制方法,其特征在于,所述步骤6中控制器为如下形式的分数阶PIλDμ控制器
G c ( s ) = K p + K I s λ + K D s μ , λ , μ > 0 ,
其中,KP、KI和KD分别是分数阶PIλDμ控制器的比例系数、积分系数和微分系数;λ、μ分别是分数阶积分阶次和微分阶次;s是拉氏变换的算子。
CN201310541805.2A 2013-11-06 2013-11-06 一种燃煤锅炉蒸汽温度的分数阶自整定控制方法 Expired - Fee Related CN103576553B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310541805.2A CN103576553B (zh) 2013-11-06 2013-11-06 一种燃煤锅炉蒸汽温度的分数阶自整定控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310541805.2A CN103576553B (zh) 2013-11-06 2013-11-06 一种燃煤锅炉蒸汽温度的分数阶自整定控制方法

Publications (2)

Publication Number Publication Date
CN103576553A CN103576553A (zh) 2014-02-12
CN103576553B true CN103576553B (zh) 2016-08-24

Family

ID=50048570

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310541805.2A Expired - Fee Related CN103576553B (zh) 2013-11-06 2013-11-06 一种燃煤锅炉蒸汽温度的分数阶自整定控制方法

Country Status (1)

Country Link
CN (1) CN103576553B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105180136A (zh) * 2015-10-08 2015-12-23 南京信息工程大学 基于分数阶pi动态矩阵的火电厂锅炉主蒸汽温度控制方法
CN105318312B (zh) * 2015-11-26 2017-09-26 广东省粤电集团有限公司 一种电站锅炉蒸汽温度的全工况多步超前预测控制方法
CN107783423B (zh) * 2017-10-25 2020-03-27 珠海格力电器股份有限公司 基于机器学习的pid参数自整定方法及其装置
CN109001975B (zh) * 2018-08-02 2021-06-15 杭州电子科技大学 一种工业加热炉多模型分数阶控制方法
JP6943826B2 (ja) * 2018-09-28 2021-10-06 三菱パワー株式会社 熱機器の内部状態推定方法及び装置
CN110879576B (zh) * 2019-10-15 2020-09-15 杭州电子科技大学 水泥熟料冷却过程的分数阶模型广义二自由度控制方法
CN112000003B (zh) * 2020-08-31 2022-09-16 新疆大学 一种基于分数阶控制器的氧化槽的温度控制方法
CN112925198B (zh) * 2021-01-23 2022-06-14 西安热工研究院有限公司 一种一维多象限燃煤机组主汽温模糊控制方法
CN113244647A (zh) * 2021-05-25 2021-08-13 扬州工业职业技术学院 基于矩阵解耦的精馏塔vofflc温控方法及系统
CN116610169B (zh) * 2023-07-19 2023-10-20 四川公路桥梁建设集团有限公司 一种目标对象智能化养护系统及方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101286044A (zh) * 2008-05-12 2008-10-15 杭州电子科技大学 一种燃煤锅炉系统混合建模方法
CN101446997A (zh) * 2008-12-31 2009-06-03 中山大学 基于粒子群算法的功率电子电路优化方法
CN101498458A (zh) * 2009-03-02 2009-08-05 杭州电子科技大学 一种锅炉在线燃烧优化的模型更新方法
CN101498457A (zh) * 2009-03-02 2009-08-05 杭州电子科技大学 一种锅炉燃烧优化的方法
CN102123115A (zh) * 2011-04-15 2011-07-13 南京信息工程大学 基于粒子群优化的正交小波盲均衡方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7006900B2 (en) * 2002-11-14 2006-02-28 Asm International N.V. Hybrid cascade model-based predictive control system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101286044A (zh) * 2008-05-12 2008-10-15 杭州电子科技大学 一种燃煤锅炉系统混合建模方法
CN101446997A (zh) * 2008-12-31 2009-06-03 中山大学 基于粒子群算法的功率电子电路优化方法
CN101498458A (zh) * 2009-03-02 2009-08-05 杭州电子科技大学 一种锅炉在线燃烧优化的模型更新方法
CN101498457A (zh) * 2009-03-02 2009-08-05 杭州电子科技大学 一种锅炉燃烧优化的方法
CN102123115A (zh) * 2011-04-15 2011-07-13 南京信息工程大学 基于粒子群优化的正交小波盲均衡方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
锅炉-汽轮机系统的分数阶控制器设计;王东风等;《中国电机工程学报》;20100215;第30卷(第5期);论文第114-117页 *

Also Published As

Publication number Publication date
CN103576553A (zh) 2014-02-12

Similar Documents

Publication Publication Date Title
CN103576553B (zh) 一种燃煤锅炉蒸汽温度的分数阶自整定控制方法
Jiang et al. MPC-based frequency control with demand-side participation: A case study in an isolated wind-aluminum power system
CN108490790A (zh) 一种基于多目标优化的过热汽温自抗扰串级控制方法
CN101709869B (zh) 燃煤锅炉过热蒸汽温度系统混合控制方法
CN106919053A (zh) 一种基于变结构预测控制算法的火电机组协调控制系统
CN100545772C (zh) 一种燃煤锅炉系统混合控制方法
CN101709863B (zh) 燃煤锅炉炉膛压力系统混合控制方法
CN112803494B (zh) 一种含风光水火的多目标agc协调优化方法及系统
CN105573123A (zh) 一种基于改进的t-s模糊预测建模的火电机组机炉协调控制方法
CN112632774A (zh) 一种基于动态模态分解的数据驱动风电场频率控制方法
Bhati et al. Genetic optimization tuning of an automatic voltage regulator system
CN102854797B (zh) 火电机组先进控制多模型切换方法
Pain et al. Multiobjective optimization of load frequency control using PSO
CN110912185B (zh) 一种含风力发电电网自动发电控制系统pid控制器设计方法
CN111245032B (zh) 一种计及风电场集电线路降损优化的电压预测控制方法
CN105490266B (zh) 基于多变量拟合的发电机调速系统参数优化建模方法
CN115313416B (zh) 一种适用于储能电站辅助调频系统的多目标优化控制方法
CN107453366B (zh) 计及风电决策风险的含upfc多目标最优潮流计算方法
Wu et al. Parameter optimization for AP1000 steam generator feedwater control system using particle swarm optimization algorithm
CN108984979B (zh) 基于多变量频域法和启发式搜索算法结合的超超临界机组深度调峰控制器设计方法
Novikov et al. SELECTING A METHOD FOR THE PARAMETRIC ADAPTATION OF PI-CONTROLLER IN THE CONTROL SYSTEMS OF BOILER ASSEMBLIES AT THERMAL POWER STATIONS WITH SUPERCRITICAL PARAMETERS.
Hou et al. Parameters optimization of ADRC based on adaptive CPSO algorithm and its application in main-steam temperature control system
CN111082442B (zh) 一种基于改进fpa的储能容量优化配置方法
Zheng et al. Double fuzzy pitch controller of wind turbine designed by genetic algorithm
CN114263899B (zh) 一种热电联产母管制供热机组的锅炉负荷分配方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for 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

Granted publication date: 20160824

Termination date: 20191106

CF01 Termination of patent right due to non-payment of annual fee