CN105157055A - 基于最小二乘的循环流化床燃烧模型辨识方法 - Google Patents
基于最小二乘的循环流化床燃烧模型辨识方法 Download PDFInfo
- Publication number
- CN105157055A CN105157055A CN201510357392.1A CN201510357392A CN105157055A CN 105157055 A CN105157055 A CN 105157055A CN 201510357392 A CN201510357392 A CN 201510357392A CN 105157055 A CN105157055 A CN 105157055A
- Authority
- CN
- China
- Prior art keywords
- square
- model
- cfbc
- squares
- estimation
- 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.)
- Pending
Links
Landscapes
- Fluidized-Bed Combustion And Resonant Combustion (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了基于最小二乘的循环流化床燃烧模型得辨识方法,其特征在于:所述控制方法利用基本最小二乘参数估计的方法,具体步骤包括:(a)模型建立;(b)置初值:对循环流化床燃烧过程建立参数估计,设置算法的初始值。(c)采样数据:导入循环流化床燃烧过程的供给锅炉的燃料量(输入)和汽压(输出)数据;(d)参数递推估计;(e)迭代收敛。本发明通过基本最小二乘的估计方法,充分考虑锅炉燃烧的大时滞、非线性、强耦合的特性,建立递推的模型,有利于计算机的实时控制与在线校正,使得所建立的模型精度大大提高,有效地提高了过程系统的反应速度和控位精度,同时易于建立其他工业对象的输入输出模型。
Description
【技术领域】
本发明涉及自动控制技术领域,特别是基于最小二乘的循环流化床燃烧模型辨识的技术领域。
【背景技术】
目前,煤炭是我国的主要一次能源,煤电占总发电量的70%以上,还有大量的工业锅炉每年要消耗几亿吨煤炭。从我国的自然资源看,这种状况在相当长的时间里都难以有根本的改变。同时,煤炭作为燃料,长期以来存在着污染严重和利用率不高的问题。正因为如此,早在20世纪60年代初我国就出现了沸腾床(鼓泡流化床)锅炉,80年代末90年代初循环流化床锅炉又得到了飞速发展。循环流化床的的主要优点有:①燃料适应性广:对难燃煤种甚至废弃物也可高效燃用;②截面热强度高,热交换效率高:对改造现有锅炉具有现实意义;③炉膛温度对应于脱硫反应的温度区间:脱硫效果好,污染物排放少;④床料是高温循环灰:负荷适应性好;⑤燃料制备相对简单:只需干燥、破碎等等。循环流化床燃烧的主要任务是使燃料所提供的热量适应负荷的需要,并同时保证经济燃烧和安全运行,即维持汽压恒定,保证燃烧经济性(过剩空气系数)和维持炉膛负压一定。其中,维持汽压恒定是衡量负荷的标准,通常通过调节供给锅炉的燃料量来调节汽压。
最小二乘法,又称最小平方法,是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据,并使得这些求得的数据与实际数据之间误差的平方和为最小。
利用最小二乘的方法对循环流化床燃烧过程进行模型的辨识,并考虑各类噪声的影响,从而使锅炉的控制更加精确和有效。
【发明内容】
本发明的目的就是解决现有技术中的问题,提出基于最小二乘的循环流化床燃烧模型的辨识方法,可以明显减小锅炉燃烧模型建立的误差,从而有效地提高过程系统的控制精度,并且易于实现计算机控制。
为实现上述目的,本发明提出了基于最小二乘的循环流化床燃烧模型辨识方法,所述控制方法充分考虑循环流化床燃烧过程的大时滞、非线性、强耦合的特性,利用基本最小二乘参数估计的方法,具体步骤包括:
(a)模型建立:考虑随机过程的CARMA模型,并将模型转换为最小二乘结构;
(b)置初值:根据基本最小二乘的一步与递推算法,对循环流化床燃烧过程建立参数估计,并设置算法的初始值。
(c)采样数据:确定采样周期T,导入循环流化床燃烧过程的供给锅炉的燃料量(输入)和汽压(输出)数据;
(d)参数递推估计:根据基本最小二乘估计的递推算法,计算和P(k);
(e)迭代收敛:返回步骤(c),再次导入输入和输出的数据,直至算法收敛或者满足要求。
作为优选,所述CARMA模型在循环流化床燃烧过程的应用中,最小二乘的结构中考虑的噪声影响为白噪声,有色噪声的影响可以看成白噪声序列驱动的线性环节的输出。
作为优选,所述最小二乘的结构模型只考虑循环流化床的主要任务,即燃料与负荷的相关关系,从而简化了辨识的模型,提高计算机辨识的效率。
作为优选,所述最小二乘类方法采用的是基本最小二乘估计方法,因为引风量、给水量等因素而对模型辨识造成的干扰可通过最小二乘结构模型中的噪声项来表述。
作为优选,所述迭代收敛过程的结束可以通过满足目标准则的要求来确定,对循环流化床燃烧模型的辨识精度要求范围更加大,自由度也更大。
本发明的有益效果:本发明通过基本最小二乘的估计方法,充分考虑锅炉燃烧的大时滞、非线性、强耦合的特性,建立递推的模型,有利于计算机的实时控制与在线校正,使得所建立的模型精度大大提高,对于之后的锅炉液位控制过程、主蒸汽压力控制过程等更加有效,明显减小被控量的超调量与方差,有效地提高过程系统的反应速度和控位精度,同时具有很强的推广性,易于建立其他工业对象的输入输出模型。
本发明的特征及优点将通过实施例进行详细说明。
【具体实施方式】
本发明,具体步骤包括:
步骤一、考虑随机过程的CARMA模型,并将模型转换为最小二乘结构;
步骤二、根据基本最小二乘的一步与递推算法,对循环流化床燃烧过程建立参数估计,并设置算法的初始值;
步骤三、确定采样周期T,导入循环流化床燃烧过程的供给锅炉的燃料量(输入)和汽压(输出)数据;
步骤四、根据基本最小二乘估计的递推算法,进行参数递推估计,计算和P(k);
步骤五、返回步骤(c),再次导入输入和输出的数据,直至算法收敛或者满足要求。
所述CARMA模型在循环流化床燃烧过程的应用中,最小二乘的结构中考虑的噪声影响为白噪声,有色噪声的影响可以看成白噪声序列驱动的线性环节的输出。所述最小二乘的结构模型只考虑循环流化床的主要任务,即燃料与负荷的相关关系,从而简化了辨识的模型,提高计算机辨识的效率。所述最小二乘类方法采用的是基本最小二乘估计方法,因为引风量、给水量等因素而对模型辨识造成的干扰可通过最小二乘结构模型中的噪声项来表述。所述迭代收敛过程的结束可以通过满足目标准则的要求来确定,对循环流化床燃烧模型的辨识精度要求范围更加大,自由度也更大。
步骤一中,建立CARMA模型为随机过程的CARMA模型:
其中,y(k)为主汽压力的实时测定离散值,u(k)为电机输出的控制作用,d为纯时延,为包括测量噪声、引风量、给水量等各类干扰的过程噪声影响,且
na,nb,nc分别为输出、输入和噪声的阶次。
假定可观测,将上述模型转换为最小二乘结构:
步骤二中,最小二乘的一步与递推算法为:
算法初始值取
步骤四中,基本最小二乘估计的递推算法进行参数递推估计:
其中,计算和P(k)。
本发明工作过程:
本发明通过基本最小二乘的估计方法,充分考虑锅炉燃烧的大时滞、非线性、强耦合的特性,建立递推的模型,有利于计算机的实时控制与在线校正,使得所建立的模型精度大大提高,对于之后的锅炉液位控制过程、主蒸汽压力控制过程等更加有效,明显减小被控量的超调量与方差,有效地提高过程系统的反应速度和控位精度,同时具有很强的推广性,易于建立其他工业对象的输入输出模型。
上述实施例是对本发明的说明,不是对本发明的限定,任何对本发明简单变换后的方案均属于本发明的保护范围。
Claims (5)
1.基于最小二乘的循环流化床燃烧模型辨识方法,其特征在于:所述控制方法利用充分考虑循环流化床燃烧过程的大时滞、非线性、强耦合的特性,利用基本最小二乘参数估计的方法,具体步骤包括:
(a)模型建立:考虑随机过程的CARMA模型,并将模型转换为最小二乘结构;
(b)置初值:根据基本最小二乘的一步与递推算法,对循环流化床燃烧过程建立参数估计,并设置算法的初始值。
(c)采样数据:确定采样周期T,导入循环流化床燃烧过程的供给锅炉的燃料量(输入)和汽压(输出)数据;
(d)参数递推估计:根据基本最小二乘估计的递推算法,计算和P(k);
(e)迭代收敛:返回步骤(c),再次导入输入和输出的数据,直至算法收敛或者满足要求。
2.如权利要求1所述的基于最小二乘的循环流化床燃烧模型辨识方法,其特征在于:(a)步骤中,所述CARMA模型在循环流化床燃烧过程的应用中,最小二乘的结构中考虑的噪声影响为白噪声,有色噪声的影响可以看成白噪声序列驱动的线性环节的输出。
3.如权利要求1所述的基于最小二乘的循环流化床燃烧模型辨识方法,其特征在于:(a)步骤中,所述最小二乘的结构模型只考虑循环流化床的主要任务,即燃料与负荷的相关关系,从而简化了辨识的模型,提高计算机辨识的效率。
4.如权利要求1所述的基于最小二乘的循环流化床燃烧模型辨识方法,其特征在于:(b)步骤中,所述最小二乘类方法采用的是基本最小二乘估计方法,因为引风量、给水量等因素而对模型辨识造成的干扰可通过最小二乘结构模型中的噪声项来表述。
5.如权利要求1所述的基于最小二乘的循环流化床燃烧模型辨识方法,其特征在于:(e)步骤中,所述迭代收敛过程的结束可以通过满足目标准则的要求来确定,对循环流化床燃烧模型辨识的精度范围和自由度更加大。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510357392.1A CN105157055A (zh) | 2015-06-24 | 2015-06-24 | 基于最小二乘的循环流化床燃烧模型辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510357392.1A CN105157055A (zh) | 2015-06-24 | 2015-06-24 | 基于最小二乘的循环流化床燃烧模型辨识方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105157055A true CN105157055A (zh) | 2015-12-16 |
Family
ID=54798014
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510357392.1A Pending CN105157055A (zh) | 2015-06-24 | 2015-06-24 | 基于最小二乘的循环流化床燃烧模型辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105157055A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106773669A (zh) * | 2016-11-21 | 2017-05-31 | 华北电力大学(保定) | 一种燃料热值实时自适应校正的火电机组协调控制方法 |
CN106873370A (zh) * | 2017-03-07 | 2017-06-20 | 云南电网有限责任公司电力科学研究院 | 一种循环流化床机组协调控制方法 |
CN106908248A (zh) * | 2017-03-22 | 2017-06-30 | 哈尔滨工程大学 | 自识别单双韦伯燃烧规则经验参数自动校准方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011122800A (ja) * | 2009-12-14 | 2011-06-23 | Sumitomo Heavy Ind Ltd | 流動床燃焼炉及び流動床燃焼炉の運転方法 |
CN102494336A (zh) * | 2011-12-16 | 2012-06-13 | 浙江大学 | 一种循环流化床锅炉燃烧过程多变量控制方法 |
CN102661759A (zh) * | 2012-05-11 | 2012-09-12 | 西安科技大学 | 循环流化床锅炉非线性多变量关键参数的辨识及预测方法 |
CN103115356A (zh) * | 2013-02-07 | 2013-05-22 | 华北电力大学 | 超临界cfb锅炉燃烧信号的监测方法及优化控制方法 |
CN104696952A (zh) * | 2015-03-25 | 2015-06-10 | 济南大学 | 一种基于模型在线辨识的流化床锅炉引风量控制方法 |
-
2015
- 2015-06-24 CN CN201510357392.1A patent/CN105157055A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011122800A (ja) * | 2009-12-14 | 2011-06-23 | Sumitomo Heavy Ind Ltd | 流動床燃焼炉及び流動床燃焼炉の運転方法 |
CN102494336A (zh) * | 2011-12-16 | 2012-06-13 | 浙江大学 | 一种循环流化床锅炉燃烧过程多变量控制方法 |
CN102661759A (zh) * | 2012-05-11 | 2012-09-12 | 西安科技大学 | 循环流化床锅炉非线性多变量关键参数的辨识及预测方法 |
CN103115356A (zh) * | 2013-02-07 | 2013-05-22 | 华北电力大学 | 超临界cfb锅炉燃烧信号的监测方法及优化控制方法 |
CN104696952A (zh) * | 2015-03-25 | 2015-06-10 | 济南大学 | 一种基于模型在线辨识的流化床锅炉引风量控制方法 |
Non-Patent Citations (1)
Title |
---|
马强: "《循环流化床锅炉自动控制与自抗扰控制器》", 31 May 2015 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106773669A (zh) * | 2016-11-21 | 2017-05-31 | 华北电力大学(保定) | 一种燃料热值实时自适应校正的火电机组协调控制方法 |
CN106873370A (zh) * | 2017-03-07 | 2017-06-20 | 云南电网有限责任公司电力科学研究院 | 一种循环流化床机组协调控制方法 |
CN106908248A (zh) * | 2017-03-22 | 2017-06-30 | 哈尔滨工程大学 | 自识别单双韦伯燃烧规则经验参数自动校准方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105157055A (zh) | 基于最小二乘的循环流化床燃烧模型辨识方法 | |
Mousafarash et al. | Exergy and exergo-economic based analysis of a gas turbine power generation system. | |
CN103853881A (zh) | 基于自适应混沌和差分进化粒子群算法的水轮机参数辨识方法 | |
CN104879750B (zh) | 一种用于循环流化床锅炉的燃烧优化设备、系统以及方法 | |
CN104504292A (zh) | 基于bp神经网络预测循环流化床锅炉最佳工作温度的方法 | |
CN108278625B (zh) | 一种燃煤锅炉高温再热器吹灰方法 | |
CN105512388A (zh) | 一种基于labview的亚临界火电机组增强激励仿真建模方法 | |
CN103440528A (zh) | 基于耗差分析的火电机组运行优化方法及装置 | |
CN105259758A (zh) | 基于海量历史数据火电机组运行参数智能在线寻优方法 | |
CN105066481A (zh) | 基于trnsys的大型太阳能集热系统模型及建模方法 | |
CN108182553A (zh) | 一种燃煤锅炉燃烧效率在线测量方法 | |
CN110298502B (zh) | 基于能效最优的锅炉最佳氧量计算方法 | |
CN110929446A (zh) | 锅炉屏式受热面的热流密度和蒸汽流量分配的预测方法 | |
CN110888401B (zh) | 火力发电机组燃烧控制优化方法、装置及可读存储介质 | |
Djayanti | Energy Efficiency Improvement Strategies for Boilers: A Case Study in Pharmacy Industry | |
CN102298321A (zh) | 基于序列模式挖掘的火电厂球磨机制粉系统节能优化方法 | |
CN104806995B (zh) | 一种用于煤质多变情况下优化锅炉运行的方法 | |
Chu et al. | Evaluating CCS Investment of China by a Novel Real Option‐Based Model | |
CN105116721A (zh) | 基于最小二乘的水箱系统辨识方法 | |
Ziemele et al. | Decomposition analysis of district heating system based on complemented Kaya identity | |
CN112036091A (zh) | 一种利用回归算法预测四角切圆燃烧器喷口风量的方法 | |
CN110848662B (zh) | 一种基于参数自适应提高暂态性能的燃煤机组汽包锅炉水位控制方法 | |
Zhou et al. | Technical performance optimization of a novel geothermal hybrid power generation system | |
Zhang et al. | Analysis and Prediction of Factors Influencing Carbon Emissions of Energy Consumption Under Climate Change | |
Xi et al. | Robust Nonlinear Adaptive Backstepping Coordinated Control for Boiler-Turbine Units |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20151216 |
|
WD01 | Invention patent application deemed withdrawn after publication |