CN104296131A - 一种双炉膛循环流化床机组的多变量协同控制方法 - Google Patents
一种双炉膛循环流化床机组的多变量协同控制方法 Download PDFInfo
- Publication number
- CN104296131A CN104296131A CN201410573371.9A CN201410573371A CN104296131A CN 104296131 A CN104296131 A CN 104296131A CN 201410573371 A CN201410573371 A CN 201410573371A CN 104296131 A CN104296131 A CN 104296131A
- Authority
- CN
- China
- Prior art keywords
- delta
- formula
- overbar
- input
- represent
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000004044 response Effects 0.000 claims description 31
- 230000008859 change Effects 0.000 claims description 18
- 239000000463 material Substances 0.000 claims description 15
- 238000012937 correction Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000005457 optimization Methods 0.000 claims description 5
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 4
- 229910052760 oxygen Inorganic materials 0.000 claims description 4
- 239000001301 oxygen Substances 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 230000000903 blocking effect Effects 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 claims description 2
- 238000001914 filtration Methods 0.000 claims description 2
- 238000002485 combustion reaction Methods 0.000 abstract description 15
- 230000006872 improvement Effects 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 3
- 230000033228 biological regulation Effects 0.000 abstract description 2
- 230000001276 controlling effect Effects 0.000 description 14
- 230000003134 recirculating effect Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 239000003245 coal Substances 0.000 description 4
- 239000000446 fuel Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 3
- 238000006477 desulfuration reaction Methods 0.000 description 3
- 239000003546 flue gas Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000003009 desulfurizing effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000002893 slag Substances 0.000 description 2
- 235000019738 Limestone Nutrition 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000023556 desulfurization Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000010881 fly ash Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000002440 industrial waste Substances 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 230000009063 long-term regulation Effects 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000010802 sludge Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F23—COMBUSTION APPARATUS; COMBUSTION PROCESSES
- F23C—METHODS OR APPARATUS FOR COMBUSTION USING FLUID FUEL OR SOLID FUEL SUSPENDED IN A CARRIER GAS OR AIR
- F23C10/00—Fluidised bed combustion apparatus
- F23C10/18—Details; Accessories
- F23C10/28—Control devices specially adapted for fluidised bed, combustion apparatus
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F23—COMBUSTION APPARATUS; COMBUSTION PROCESSES
- F23C—METHODS OR APPARATUS FOR COMBUSTION USING FLUID FUEL OR SOLID FUEL SUSPENDED IN A CARRIER GAS OR AIR
- F23C2206/00—Fluidised bed combustion
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种双炉膛循环流化床机组的多变量协同控制方法,它是构建循环流化床锅炉和汽轮机机组的一个多变量系统,采用基于跟踪误差变化率的多变量约束DMC算法。所采用的算法针对传统的DMC算法做出了改进,在目标函数中引入了跟踪误差变化率,提高了对设定值的跟踪效果,同时对输入输出附加了约束,使得求解的控制量在其物理范围内。由于控制系统同步考虑了左右两侧炉膛的床温和床压,在变工况调节时既能保证负荷跟踪能力,又能实现两个炉膛的床温床压平衡稳定和均衡燃烧。
Description
技术领域
本发明属于循环流化床发电厂控制领域,涉及循环流化床锅炉负荷、床温和床压的控制方法。
背景技术
循环流化床(Circulating Fluidized Bed,CFB)的锅炉燃烧技术作为高效能、低污染、燃烧煤泥矸石等劣质煤、适应性广的洁净燃煤技术,在全世界越来越受到广泛重视,正在成为燃煤技术的主力军。国家“十二五”规划明确提出“加大节能力度;推进粉煤灰、煤矸石、冶金和化工废渣及尾矿等工业废物利用”,给大型循环流化床锅炉的技术研究和产业化发展带来了的机遇,目前300MW循环流化床锅炉已逐步成为国内重要的发电锅炉形式,600MW超临界循环流化床锅炉也已有试点运行。
CFB锅炉炉膛内含有大量的物料,在燃烧过程中大量的物料被烟气携带到炉膛上部,经过布置在炉膛出口的分离器,将物料与烟气分开,并经过非机械式回料阀将物料直接回送至炉内或经外置床受热面后送至炉内,多次循环燃烧。由于物料浓度高,具有很大的热容量和良好的物料混合,一般每公斤烟气可携带若干公斤的物料,这些循环物料带来了高传热系数,使锅炉热负荷调节范围广,对燃料的适应性强。循环流化床锅炉的循环流化方式增加了燃料在炉膛内的停留时间,因此具有较高的燃烧效率,在低负荷下能稳定运行,而无需增加辅助燃料。
CFB锅炉运行温度通常在800~900℃之间,这是一个理想的脱硫温度区间,采用炉内脱硫技术,向床内加入石灰石和脱硫剂,燃料及脱硫剂经多次循环,反复进行低温燃烧和脱硫反应,加之炉内湍流运动剧烈,Ca/S摩尔比约为2时,可以使脱硫效率达到90%左右,SO2的排放量大大降低。同时循环流化床采用低温分级送风燃烧,使燃烧始终在低过量空气下进行,从而大大降低了NOx的生成和排放。
循环流化床锅炉的控制点多,燃烧系统耦合性强,燃烧响应滞后性大是循环流化床控制的难点,也是循环流化床锅炉本体结构和燃烧煤种所决定的。循环流化床控制难的问题是影响循环流化床锅炉又好又快发展的一个重要因素,如何解决好循环流化床特别是大容量循环床锅炉燃烧控制的稳定性、鲁棒性问题是一个迫切需要解决的问题。循环流化床锅炉的主要控制问题需要解决变工况时的床温和床压控制的稳定性,负荷控制的升降速率保证,以及国内主流300MW循环流化床锅炉双炉膛的床压平衡和均衡燃烧问题。
目前针对循环流化床锅炉建模和控制领域的研究工作主要体现在两个层面:一是通过分析循环流化床锅炉运行机理,建立锅炉机理模型,或基于现场试验所测数据,对不同工况下各变量参数之间的关系进行分析,建立控制系统设计所需要的控制模型。二是针对循环流化床锅炉燃烧系统的控制特性,研究和运用常规PID改进方法或先进控制方法进行控制。
对于循环流化床机组,床温和床压是两个重要的运行参数。当机组的负荷指令发生变化时,一、二次风量和给煤量均发生变化,从而对床温床压产生较大影响。传统的单回路PID控制方法是通过外置床回料阀调节单侧床温,负荷变化时时床温控制是难题;另外通过调节冷渣器转速调节单侧床压,但这种调节方式十分缓慢,运行人员通常只把冷渣器转速作为稳态长期调节床压的手段,不能较好满足动态调节需求。
同时300MW主流循环流化床机组的双炉膛结构设计,两侧炉膛的均衡燃烧和床温床压平衡是实际运行中的难题,特别是两侧床压受一次风的影响很大,变工况时容易产生失稳翻床。因此传统的单侧炉膛单回路控制附加一定的双向偏差补偿的控制方法只能满足稳态和小范围变化的调节需求。
发明内容
要解决的技术问题:针对现有技术的不足,本发明提出一种双炉膛循环流化床机组的多变量协同控制方法,用于解决变工况时循环流化床的床温、床压以及负荷协调控制的稳定性差以及双炉膛均衡燃烧控制困难的技术问题。
技术方案:为解决上述技术问题,本发明采用以下技术方案:
一种双炉膛循环流化床机组的多变量协同控制方法,以左炉一次风量、左炉二次风量、左炉给煤量、左炉回料阀开度、主汽调门开度、右炉一次风量、右炉二次风量、右炉给煤量和右炉回料阀开度共9个参数作为输入量,同时以实发功率、机前压力、氧量、左炉床温、右炉床温、左炉床压和右炉床压共7个参数作为输出量,基于跟踪误差变化率的多变量约束DMC算法进行协同控制。
DMC(动态矩阵控制)算法是一种基于对象阶跃响应的预测控制算法,它适用于渐进稳定的线性对象。DMC算法不要求对模型的结构有先验知识,适用于竖向模型难以精确建立的复杂对象,因而该算法在工过程中比较实用。一般的DMC控制都包括以下三部分:DMC算法的模型预测、DMC算法的滚动优化以及DMC算法的反馈校正。
进一步,在本发明中,针对双炉膛循环流化床机组的多变量协同控制,对传统的DMC算法进行了特别的优化和改进,以u1,u2,...,u9对应表示9个输入量并作为控制量,以y1,y2,...,y7对应表示7个输出量,在目标函数中考虑跟踪误差变化率,并考虑输入和输出约束,整个控制方法具体包括如下步骤:
步骤一、获取双炉膛循环流化床机组的阶跃响应模型:
在稳态工况下,以9个输入量为输入对7个输出量进行阶跃响应实验,经滤波平滑后,得到7个输出量相对于9个输入量的阶跃响应系数分别为 并满足 下标Nij表示第j个输入对第i个输出的阶跃响应模型的截断长度,截断长度的选择满足保证输出量的阶跃响应值已经达到稳态;
以上述阶跃响应系数的组合构成第j个输入与第i个输出之间的阶跃响应模型;
步骤二、设置控制器的相关参数,包括采样时间Ts(采样时间一般取10~15s)、各输出量的预测时域Pi(1≤i≤7)、各输入量的控制时域Mj(1≤j≤9)、误差校正系数跟踪误差权矩阵控制增量权矩阵以及跟踪误差变化率权矩阵其中,预测时域Pi、控制时域Mj、误差校正系数和三个权矩阵需要在仿真中进行调试,并根据现场的运行数据再做调整;
步骤三、将步骤一得到的阶跃响应模型带入到下列(1)式中,求解出7个输出数量的预测值
(1)式中,
预测值表示为
其中每个输入量的预测值表示为
Y0(k)是在控制作用不变下的自由响应,并且有
其中每个输入量的自由响应表示为
是动态矩阵,并且
动态矩阵中每个子矩阵表示如下
控制增量ΔU(k)表示为
其中每个控制增量表示为
求解时,在0时刻,所有的预测值初值取0;
步骤四、求解自由响应Y0(k)
首先计算出预测误差并进行反馈校正:
由于第i个输出量在k+1时刻的预测误差ei(k+1)由下列(2)式表示
(2)式中,yi(k+1)表示第i个输出量在k+1时刻的实际测量值;
所以第i个输出量在k+1时刻校正后的预测值为
(3)式中,
误差校正系数表示为
因此求解出所有输出量在k+1时刻校正后的预测值为
(4)式中,
E(k)=[e1(k+1),e2(k+1),...,e7(k+1)]T
然后求解自由响应:
用校正后的预测值表示第i个输出量的自由响应的第a个分量,即yi0(k+a)=yicor(k+a+1)a∈[1,Pi-1] (2)
特别的,用yicor(k+Pi)来近似的最后一项yi0(k+Pi),结合(5)式,表示为(6)式
(6)式中
因此根据(6)式将所有输出量的自由响应表示为(7)式
(7)式中
步骤五、设定目标函数和约束,求解控制增量ΔU(k),具体步骤如下:
首先将目标函数设定为
(8)式中
Epf(k)表示跟踪误差;ΔEpf(k)表示跟踪误差变化率;
表示跟踪误差的惩罚权重;表示跟踪误差变化率的惩罚权重;表示控制增量的惩罚权重;它们的具体形式如下
由于
并且对有Δuj(k-1)=0,所以
对所有的输出量,根据上式得到(10)式
ΔU(k-1)=AΔU(k) (10)
(10)式中
A=diag(α1,α2,α3,...α9)
结合(9)式和(10)式,并令则将转化为用ΔU(k)表示的形式
因此将Epf(k)和ΔEpf(k)转化为用ΔU(k)表示的形式
其中,W(k)表示各个输出量的设定值,前缀Δ表示各个量的变化值;
然后对输入量和输出量分别设置约束如下
输入量约束包括式(11)和式(12)
控制增量ΔU(k)满足
控制量U(k)满足
式(11)中,ΔU表示控制增量的下限,表示控制增量的上限;
式(12)中,U表示控制量的下限,表示控制量的上限;
具体的控制增量和控制量的上下限由实际需求和现场设备的物理极限决定,没有限制时可取为无穷;
下面将式(12)中控制量U(k)的约束转化为(15)式中用控制增量ΔU(k)来表示的形式:
由于对于第j个输入量的控制量uj有
uj(k)=uj(k-1)+Δuj(k)
uj(k+1)=uj(k-1)+Δuj(k)+Δuj(k+1)
.
.
.
uj(k+Mj-1)=uj(k-1)+Δuj(k)+Δuj(k+1)+…+Δuj(k+Mj-1)
将上述控制量uj的各式综合表示成(13)式的形式
(13)式中
所以对所有输入量,基于(13)式可得
(14)式中
结合(12)式和(14)式得到(15)式
利用模型的预测值代替实际的输出量,(16)式表示输出量的约束
(16)式中,Y表示实际输出量的下限,表示实际输出量的上限
结合(1)式和(16)式,将输出量的约束转化成为(17)式中用控制增量ΔU(k)来表示的形式:
结合式(11)、(15)和(17),将输入量的约束和输出量的约束综合表示成(18)式的形式
(18)式中
将(8)式变形为如下二次型的形式
二次型中
f0=(Y0-W)T(Y0-W)+(ΔY0-ΔW)T(ΔY0-ΔW)
由于f0与控制增量ΔU(k)无关,在优化时忽略;将ΔU(k)的求解转化为标准的二次优化问题,形式如下
利用二次规划问题的有效集算法对(19)式进行求解,得到控制增量ΔU(k);
步骤六、利用U(k)=U(k-1)+ΔU(k)分别计算9个控制量的值;
步骤七、输出新的控制量,并更新(1)式中的预测值,并重复执行步骤三到步骤七。
有益效果:
本发明将一种多变量预测控制算法应用在双炉膛循环流化床机组上,创新过程包含两个方面:(1)提取构建多变量对象模型,并在DMC算法的目标函数引入跟踪误差变化率,同时引入输入输出约束,提出了基于跟踪误差变化率的多变量约束DMC算法;(2)首次将这样一种控制方法应用到300MW双炉膛循环流化床机组的研究。具体的,主要包括以下4个方面的优点:
1、提取构建一个循环流化床锅炉和汽轮机整个机组的9个输入和7个输出的对象模型,采用基于跟踪误差变化率的多变量约束DMC算法,同时控制9个输入量,实现了床温、床压和负荷的有机协调,克服了传统PID单回路控制分别调节解藕难的缺陷;
2、所采用的算法针对传统的DMC算法做出了改进,在目标函数中引入了跟踪误差变化率,提高了跟踪效果,因为传统的DMC算法目标函数中只包含跟踪误差,这样在求解的过程中只能保证最终稳定时,输出量和设定值没有偏差,而不能保证在调节过程中输出量紧紧跟随设定值,但引入跟踪误差变化率可解决这一问题;
3、对输入量、输出量附加了约束,并用二次优化算法进行求解,使得控制系统能在保证跟踪效果,使控制量波动最小的基础上,让求解的控制量在其物理范围内,避免了直接附加约束造成的经济性损失;
4、同步考虑了左右两个炉膛的床温和床压,能有效保证两个炉膛床温和床压平衡,实现双炉膛的均衡燃烧,并将功率、机前压力、氧量控制和床温、床压控制协同考虑,构建多变量预测控制模型,能有效保证机组变工况时,床温床压协同负荷跟踪变化。
附图说明
图1为本发明的控制系统的示意图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
在图1中,系统的输出Y中一共有7个输出量,分别是实发功率NE、机前压力PT、氧量O、左炉床温TLF、右炉床温TRF、左炉床压PLF和右炉床压PRF,前缀r表示他们的设定值;系统的输入U一共有输入量有9个,分别是左炉一次风量F1,LF、左炉二次风量F2,LF、左炉给煤量BLF、左炉回料阀开度μv,LF、主汽调门开度μT、右炉一次风量F1,RF、右炉二次风量F2,RF、右炉给煤量BRF和右炉回料阀开度μv,RF。这样,系统的输入可以表示为U=[F1,LF,F2,LF,BLF,μv,LF,μT,F1,RF,F2,RF,BRF,μv,RF]T,系统的输出可表示为Y=[NE,PT,O,TLF,TRF,PLF,PRF]T,系统的输出设定值可表示为W=[r_NE,r_PT,r_O,r_TLF,r_TRF,r_PTF,r_PRF]T,之后,按照前面所述从步骤一到步骤七,便是本发明的具体工作流程。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (2)
1.一种双炉膛循环流化床机组的多变量协同控制方法,其特征在于:以左炉一次风量、左炉二次风量、左炉给煤量、左炉回料阀开度、主汽调门开度、右炉一次风量、右炉二次风量、右炉给煤量和右炉回料阀开度共9个参数作为输入量,同时以实发功率、机前压力、氧量、左炉床温、右炉床温、左炉床压和右炉床压共7个参数作为输出量,基于跟踪误差变化率的约束DMC算法进行协同控制。
2.根据权利要求1所述的双炉膛循环流化床机组的多变量协同控制方法,其特征在于:以u1,u2,…,u9对应表示9个输入量并作为控制量,y1,y2,…,y7对应表示7个输出量,在目标函数中考虑跟踪误差变化率,并考虑输入和输出约束,具体包括如下步骤:
步骤一、获取双炉膛循环流化床机组的阶跃响应模型
在稳态工况下,以9个输入量为输入对7个输出量进行阶跃响应实验,经滤波平滑后,得到7个输出量相对于9个输入量的阶跃响应系数分别为 并满足 下标Nij表示第j个输入对第i个输出的阶跃响应模型的截断长度,截断长度的选择满足保证输出量的阶跃响应值已经达到稳态;
以上述阶跃响应系数的组合构成第j个输入与第i个输出之间的阶跃响应模型;
步骤二、设置控制器的相关参数
包括采样时间Ts、各输出量的预测时域Pi(1≤i≤7)、各输入量的控制时域Mj(1≤j≤9)、误差校正系数跟踪误差权矩阵控制增量权矩阵以及跟踪误差变化率权矩阵
步骤三、将步骤一得到的阶跃响应模型带入到下列(1)式中,求解出7个输出数量的预测值
(1)式中,
预测值表示为
其中每个输入量的预测值表示为
Y0(k)是在控制作用不变下的自由响应,并且有其中每个输入量的自由响应表示为
其中所有自由响应的初值取0;
是动态矩阵,并且
动态矩阵中每个子矩阵表示如下
控制增量ΔU(k)表示为
其中每个控制增量表示为
求解时,在0时刻,所有的预测值初值取0;
步骤四、求解自由响应Y0(k)
首先计算出预测误差并进行反馈校正:
由于第i个输出量在k+1时刻的预测误差ei(k+1)由下列(2)式表示
(2)式中,yi(k+1)表示第i个输出量在k+1时刻的实际测量值;
所以第i个输出量在k+1时刻校正后的预测值为
(3)式中,
误差校正系数表示为
因此求解出所有输出量在k+1时刻校正后的预测值为
(4)式中,
E(k)=[e1(k+1),e2(k+1),...,e7(k+1)]T
然后求解自由响应:
用校正后的预测值表示第i个输出量的自由响应的第a个分量,即yi0(k+a)=yicor(k+a+1) a∈[1,Pi-1] (1)
特别的,用yicor(k+Pi)来近似的最后一项yi0(k+Pi),结合(5)式,表示为(6)式
(6)式中
因此根据(6)式将所有输出量的自由响应表示为(7)式
(7)式中
步骤五、设定目标函数和约束,求解控制增量ΔU(k),具体步骤如下
首先将目标函数设定为
(8)式中
Epf(k)表示跟踪误差;ΔEpf(k)表示跟踪误差变化率;
表示跟踪误差的惩罚权重;表示跟踪误差变化率的惩罚权重;表示控制增量的惩罚权重;它们的具体形式如下
由于
并且对有Δuj(k-1)=0,所以
对所有的输出量,根据上式得到(10)式
ΔU(k-1)=AΔU(k) (10)
(10)式中
A=diag(α1,α2,α3,...α9)
结合(9)式和(10)式,并令则将转化为用ΔU(k)表示的形式
因此将Epf(k)和ΔEpf(k)转化为用ΔU(k)表示的形式
其中,W(k)表示各个输出量的设定值,前缀Δ表示各个量的变化值;
然后对输入量和输出量分别设置约束如下
输入量约束包括式(11)和式(12)
控制增量ΔU(k)满足
控制量U(k)满足
式(11)中,ΔU表示控制增量的下限,表示控制增量的上限;
式(12)中,U表示控制量的下限,表示控制量的上限;
下面将式(12)中控制量U(k)的约束转化为(15)式中用控制增量ΔU(k)来表示的形式:
由于对于第j个输入量的控制量uj有
uj(k)=uj(k-1)+Δuj(k)
uj(k+1)=uj(k-1)+Δuj(k)+Δuj(k+1)
·
·
·
uj(k+Mj-1)=uj(k-1)+Δuj(k)+Δuj(k+1)+…+Δuj(k+Mj-1)
将上述控制量uj的各式综合表示成(13)式的形式
(13)式中
所以对所有输入量,基于(13)式可得
(14)式中
结合(12)式和(14)式得到(15)式
利用模型的预测值代替实际的输出量,(16)式表示输出量的约束
(16)式中,Y表示实际输出量的下限,表示实际输出量的上限
结合(1)式和(16)式,将输出量的约束转化成为(17)式中用控制增量ΔU(k)来表示的形式:
结合式(11)、(15)和(17),将输入量的约束和输出量的约束综合表示成(18)式的形式
(18)式中
将(8)式变形为如下二次型的形式
二次型中
f0=(Y0-W)T(Y0-W)+(ΔY0-ΔW)T(ΔY0-ΔW)
由于f0与控制增量ΔU(k)无关,在优化时忽略;将ΔU(k)的求解转化为标准的二次优化问题,形式如下
利用二次规划问题的有效集算法对(19)式进行求解,得到控制增量ΔU(k);
步骤六、利用U(k)=U(k-1)+ΔU(k)分别计算9个控制量的值;
步骤七、输出新的控制量,并更新(1)式中的预测值,并重复执行步骤三到步骤七。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410573371.9A CN104296131B (zh) | 2014-10-23 | 2014-10-23 | 一种双炉膛循环流化床机组的多变量协同控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410573371.9A CN104296131B (zh) | 2014-10-23 | 2014-10-23 | 一种双炉膛循环流化床机组的多变量协同控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104296131A true CN104296131A (zh) | 2015-01-21 |
CN104296131B CN104296131B (zh) | 2015-09-30 |
Family
ID=52315973
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410573371.9A Expired - Fee Related CN104296131B (zh) | 2014-10-23 | 2014-10-23 | 一种双炉膛循环流化床机组的多变量协同控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104296131B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104932451A (zh) * | 2015-04-17 | 2015-09-23 | 王玉娇 | 一种多设备协同控制方法及装置 |
CN105135425A (zh) * | 2015-09-17 | 2015-12-09 | 山西大学 | 带液阻调速器的一次风机控制系统及其控制方法 |
CN107270283A (zh) * | 2017-07-17 | 2017-10-20 | 安徽电气工程职业技术学院 | 一种基于循环流化床机组的新型多变量约束预测控制方法 |
CN107831651A (zh) * | 2017-10-12 | 2018-03-23 | 国网河北能源技术服务有限公司 | 煤粉电站锅炉scr与sncr脱硝系统协同控制方法 |
CN108897215A (zh) * | 2018-06-25 | 2018-11-27 | 西北工业大学 | 基于分布式模型预测控制的多海洋机器人协同环形扫描方法 |
CN110887038A (zh) * | 2019-12-26 | 2020-03-17 | 哈尔滨工业大学中远工控有限公司 | 一种循环流化床锅炉燃烧自适应控制系统及方法 |
CN112684699A (zh) * | 2020-11-10 | 2021-04-20 | 北京京能科技有限公司 | 变量幅值约束的超临界机组运行的控制方法及超临界机组 |
CN112742187A (zh) * | 2020-12-10 | 2021-05-04 | 山西漳山发电有限责任公司 | 一种脱硫系统中pH值的控制方法及装置 |
CN112797402A (zh) * | 2020-12-31 | 2021-05-14 | 华电国际电力股份有限公司天津开发区分公司 | 一种裤衩腿型循环流化床锅炉机组及床压预测系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103605284A (zh) * | 2013-11-14 | 2014-02-26 | 杭州电子科技大学 | 动态矩阵控制优化的废塑料裂解炉炉膛压力控制方法 |
-
2014
- 2014-10-23 CN CN201410573371.9A patent/CN104296131B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103605284A (zh) * | 2013-11-14 | 2014-02-26 | 杭州电子科技大学 | 动态矩阵控制优化的废塑料裂解炉炉膛压力控制方法 |
Non-Patent Citations (2)
Title |
---|
周武: "循环流化床锅炉燃烧建模与优化控制", 《上海交通大学硕士学位论文》, 1 December 2011 (2011-12-01), pages 46 - 49 * |
张强 等: "基于DMC-PID的循环流化床锅炉床层温度控制设计与实现", 《化工自动化及仪表》, vol. 36, no. 6, 10 December 2009 (2009-12-10) * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104932451A (zh) * | 2015-04-17 | 2015-09-23 | 王玉娇 | 一种多设备协同控制方法及装置 |
CN105135425A (zh) * | 2015-09-17 | 2015-12-09 | 山西大学 | 带液阻调速器的一次风机控制系统及其控制方法 |
CN107270283A (zh) * | 2017-07-17 | 2017-10-20 | 安徽电气工程职业技术学院 | 一种基于循环流化床机组的新型多变量约束预测控制方法 |
CN107831651A (zh) * | 2017-10-12 | 2018-03-23 | 国网河北能源技术服务有限公司 | 煤粉电站锅炉scr与sncr脱硝系统协同控制方法 |
CN108897215A (zh) * | 2018-06-25 | 2018-11-27 | 西北工业大学 | 基于分布式模型预测控制的多海洋机器人协同环形扫描方法 |
CN110887038A (zh) * | 2019-12-26 | 2020-03-17 | 哈尔滨工业大学中远工控有限公司 | 一种循环流化床锅炉燃烧自适应控制系统及方法 |
CN112684699A (zh) * | 2020-11-10 | 2021-04-20 | 北京京能科技有限公司 | 变量幅值约束的超临界机组运行的控制方法及超临界机组 |
CN112742187A (zh) * | 2020-12-10 | 2021-05-04 | 山西漳山发电有限责任公司 | 一种脱硫系统中pH值的控制方法及装置 |
CN112797402A (zh) * | 2020-12-31 | 2021-05-14 | 华电国际电力股份有限公司天津开发区分公司 | 一种裤衩腿型循环流化床锅炉机组及床压预测系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104296131B (zh) | 2015-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104296131B (zh) | 一种双炉膛循环流化床机组的多变量协同控制方法 | |
Sun et al. | Enhanced decentralized PI control for fluidized bed combustor via advanced disturbance observer | |
CN101750964B (zh) | 烧成系统基于脉冲响应模型的多变量预测控制方法 | |
CN103339441B (zh) | 用于氧燃料燃烧动力设备的优化综合控制 | |
CN108227488B (zh) | 基于滑模预测控制的超超临界火电机组协调控制方法 | |
CN100555125C (zh) | 气化炉的预测型比例积分微分控制方法 | |
CN102053562B (zh) | 一种裂解炉出口温度混合控制方法 | |
CN105955210B (zh) | 余热锅炉与工业锅炉联合发电系统的动态优化方法 | |
CN102425790B (zh) | 一种循环流化床锅炉在线优化自学习控制方法 | |
CN104765350A (zh) | 基于组合模型预测控制技术的水泥分解炉控制方法及系统 | |
CN109224815A (zh) | 一种基于多变量约束区间预测控制的氨法脱硫优化控制方法 | |
Jin et al. | Dynamic modeling and control for pulverized-coal-fired oxy-combustion boiler island | |
He et al. | Development and implementation of advanced control strategies for power plant cycling with carbon capture | |
CN107016176A (zh) | 一种混合智能锅炉综合燃烧优化方法 | |
CN101556038A (zh) | 循环流化床锅炉稳定运行与经济燃烧优化控制系统 | |
CN101498934B (zh) | 整体煤气化联合循环电站的协调控制方法 | |
CN106292277A (zh) | 基于全局滑模控制的亚临界火电机组协调控制方法 | |
Kocaarslan et al. | A fuzzy logic controller application for thermal power plants | |
CN105180136A (zh) | 基于分数阶pi动态矩阵的火电厂锅炉主蒸汽温度控制方法 | |
Niva et al. | Self-optimizing control structure design in oxy-fuel circulating fluidized bed combustion | |
CN102269973B (zh) | 面向烧结烟气脱硫过程的塔顶温度模糊预测补偿控制方法 | |
Liu et al. | Recent developments in applied thermal engineering: Process integration, heat exchangers, enhanced heat transfer, solar thermal energy, combustion and high temperature processes and thermal process modelling | |
Tang et al. | Distributed predictive control guided by intelligent reboiler steam feedforward for the coordinated operation of power plant-carbon capture system | |
CN102021254B (zh) | 一种高炉系统与trt装置分布式协调控制方法 | |
CN207046891U (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: 20150930 |