CN101709863B - 燃煤锅炉炉膛压力***混合控制方法 - Google Patents
燃煤锅炉炉膛压力***混合控制方法 Download PDFInfo
- Publication number
- CN101709863B CN101709863B CN2009101557929A CN200910155792A CN101709863B CN 101709863 B CN101709863 B CN 101709863B CN 2009101557929 A CN2009101557929 A CN 2009101557929A CN 200910155792 A CN200910155792 A CN 200910155792A CN 101709863 B CN101709863 B CN 101709863B
- Authority
- CN
- China
- Prior art keywords
- proportional
- model
- derivative controller
- furnace pressure
- integral
- 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
Links
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明涉及一种燃煤锅炉炉膛压力***混合控制方法。本发明方法首先基于燃煤锅炉炉膛压力实时过程数据建立过程模型,挖掘出基本的过程特性;然后基于该过程模型建立比例积分微分控制回路;最后通过计算预测PI控制器的参数,将比例积分微分控制与燃煤锅炉炉膛压力对象整体实施预测PI控制。本发明方法弥补了传统控制的不足,并有效地方便了控制器的设计,保证控制性能的提升,同时满足给定的生产性能指标。本发明提出的控制技术可以有效减少理想炉膛压力工艺参数与实际炉膛压力工艺参数之间的误差,进一步弥补了传统控制器的不足,同时保证控制装置操作在最佳状态,使生产过程的炉膛压力工艺参数达到严格控制。
Description
技术领域
本发明属于自动化技术领域,涉及一种燃煤锅炉炉膛压力***的预测比例积分控制(预测PI)与比例积分微分控制(PID)的混合控制方法。
背景技术
燃煤锅炉是电力生产部门的重要动力设备,其要求是供给合格的蒸汽,使燃煤锅炉发汽量适应负荷的需要。为此,生产过程的各个主要工艺参数必须严格控制。然而燃煤锅炉设备是一个复杂的被控对象,输入量与输出量之间相互关联。对于锅炉炉膛压力***来说:蒸汽负荷发生变化引起蒸汽压力和过热蒸汽温度变化的同时,也引起炉膛压力发生变化;燃料量的变化直接影响蒸汽压力,过热蒸汽温度、过剩空气和炉膛负压的变化;减温水的变化会导致过热蒸汽温度、蒸汽压力发生变化,进一步导致炉膛压力等的变化。这些不利因素导致传统的控制手段精度不高,又进一步导致后续生产控制参数不稳定,产品合格率低,锅炉效率低下。目前实际工业中燃煤锅炉的炉膛压力控制基本上采用传统的简单的控制手段,甚至手动操作,控制参数完全依赖技术人员经验,使生产成本增加,控制效果很不理想。我国燃煤锅炉控制与优化技术比较落后,能耗居高不下,控制性能差,自动化程度低,很难适应节能减排以及间接环境保护的需求,这其中直接的影响因素之一便是燃煤锅炉***的控制方案问题。
发明内容
本发明的目标是针对现有的燃煤锅炉炉膛压力***控制技术的不足之处,提供一种燃煤锅炉炉膛压力***混合控制方法,具体是基于预测比例积分与比例积分微分控制的混合控制方法。该方法弥补了传统控制方式的不足,保证控制具有较高的精度和稳定性的同时,也保证形式简单并满足实际工业过程的需要。
本发明方法首先基于燃煤锅炉炉膛压力实时过程数据建立过程模型,挖掘出基本的过程特性;然后基于该过程模型建立比例积分微分控制回路;最后通过计算预测PI控制器的参数,将比例积分微分控制与燃煤锅炉炉膛压力对象整体实施预测PI控制。
本发明的技术方案是通过数据采集、过程辨识、预测机理、数据驱动、优化等手段,确立了一种燃煤锅炉炉膛压力***的预测PI与比例积分微分控制的混合控制方法,利用该方法可有效提高控制的精度,,同时满足给定的生产性能指标。
本发明方法的步骤包括:
(1)利用燃煤锅炉炉膛压力实时过程数据建立过程模型,具体方法是:
首先建立燃煤锅炉炉膛压力实时运行数据库,通过数据采集装置采集N组实时过程运行数据,将采集的实时过程运行数据作为数据驱动的样本集合,表示为{Xi,y(i)}i=1 N,i=1,2,…,N,其中Xi表示第i组工艺参数的输入数据,y(i)表示第i组工艺参数的输出值。
然后以该炉膛压力实时过程运行数据集合为基础建立基于最小二乘法的离散差分方程形式的局部受控自回归滑动平均模型:
yL(k)=ΦTX,Φ=[a′1,a′2,…,a′n,b′0,b′1,…,b′m-1]T
X=[y(k-1),…,y(k-n),u(k-d-1),…,u(k-d-m)]T
其中,yL(k)表示当前时刻过程模型的工艺参数的输出值,X表示过程模型的工艺参数的过去时刻的输入和输出数据的集合,u(k)表示当前过程模型工艺参数对应的控制变量,k为当前的递推步数,Φ表示通过辨识得到的模型参数的集合,T表示矩阵的转置,n,m,d+1分别为对应实际过程的输出变量阶次、输入变量阶次、实际过程的时滞。
采用的辨识手段为:
(2)采用典型的响应曲线法设计炉膛压力过程模型的比例积分微分控制器,具体方法是:
a.将过程模型的比例积分微分控制器停留在手动操作状态,操作拨盘使其输出有阶跃变化,由记录仪表记录过程模型的输出值,将过程模型输出值yL(k)的响应曲线转换成无量纲形式yL *(k),具体是:
其中,yL(∞)是过程模型的比例积分微分控制器的输出有阶跃变化时的过程模型输出yL(k)的稳态值。
K=yL(∞)/q
T=2(k1-k2)
τ=2k1-k2
其中,q为过程模型的比例积分微分控制器输出的阶跃变化幅度。
c.计算过程模型的比例积分微分控制器的参数,具体是:
Kc=1.2T/Kτ
Ti=2τ
Td=0.5τ
其中Kc为比例积分微分控制器的比例参数,Ti为比例积分微分控制器的积分参数,Td分别为比例积分微分控制器的微分参数。
(3)设计预测比例积分比例积分微分控制器,具体步骤是:
d.将过程模型的比例积分微分控制器停留在自动操作状态,操作拨盘使其输入有阶跃变化,由记录仪表记录实时过程的输出,将过程输出值y(k)的响应曲线转换成无量纲形式y*(k),具体是:y*(k)=y(k)/y(∞)
其中,y(∞)是过程模型的比例积分微分控制器的输入有阶跃变化时的过程模型输出y(k)的稳态值。
e.选取满足y(k3)=0.39,y(k4)=0.63的另两个计算点k3和k4,依据下式计算预测比例积分比例积分微分控制器所需要的参数K1,T1和τ1:
K1=y(∞)/q1
T1=2(k3-k4)
τ1=2k3-k4
其中,q1为过程模型的比例积分微分控制器输入的阶跃变化幅度。
f.将步骤e得到的参数转化为拉普拉斯形式的局部受控传递函数模型:
其中,s为拉普拉斯变换算子,λ1为局部受控传递函数模型的时间常数,L1为局部受控传递函数模型的时滞,y(s)表示当前时刻过程模型的输出值的拉普拉斯变换,q1(s)表示过程模型的比例积分微分控制器输入的拉普拉斯变换。
λ1=T1
L1=τ1
g.依据步骤f计算出的模型参数整定预测比例积分比例积分微分控制器的参数,具体方法是:
①对该对象设计预测比例积分控制器。选期望的闭环传递函数模型为Gq2(s)
λ2为期望的闭环传递函数模型的时间常数,L2为期望的闭环传递函数模型的时滞,L2=L1;
②预测比例积分比例积分微分控制器的传递函数Gc1(s)可由下式来表示
③依据步骤②得到当前的预测比例积分比例积分微分控制器参数值u(s)
本发明提出的一种基于数据驱动的模型选取和预测PI-PID混合控制方法弥补了传统控制的不足,并有效地方便了控制器的设计,保证控制性能的提升,同时满足给定的生产性能指标。
本发明提出的控制技术可以有效减少理想炉膛压力工艺参数与实际炉膛压力工艺参数之间的误差,进一步弥补了传统控制器的不足,同时保证控制装置操作在最佳状态,使生产过程的炉膛压力工艺参数达到严格控制。
具体实施方式
以循环流化床锅炉***炉膛压力过程控制为例:
这里以该***炉膛压力回路的控制作为例子加以描述。炉膛压力不仅受到空气流量的影响,同时也受燃料流量,减温水流量和蒸汽流量的影响。调节手段采用进空气量,其余的影响作为不确定因素。
(1)建立该循环流化床锅炉***的炉膛压力过程模型。
通过数据采集装置采集实时过程炉膛压力运行数据,将采集的实时过程炉膛压力运行数据作为数据驱动的样本集合采用最小二乘法推理,建立基于最小二乘法的离散差分方程形式的炉膛压力过程模型。
其中,***调用推理机采用最小二乘法进行炉膛压力过程模型参数的辨识,这些参数包括元素Φ中变量的个数和具体数值。
其中y(k)是实际炉膛压力测量值,Φk TXk是炉膛压力过程模型的输出值。
这个过程是第一步推理过程。这个第一步推理是初步挖掘实际炉膛压力回路的基本特性。
(2)设计炉膛压力过程模型的比例积分微分控制器,具体方法是典型的响应曲线法。
第一步:将炉膛压力比例积分微分控制器停留在“手动操作”状态,操作进空气量的拨盘使进空气量控制器输出有个阶跃变化,由记录仪表记录炉膛压力过程模型的输出值,将炉膛压力过程模型输出值yL(k)的响应曲线转换成无量纲形式yL *(k):
其中,yL(∞)是炉膛压力过程模型输出yL(k)的稳态值。
K=yL(∞)/q
T=2(k1-k2)
τ=2k1-k2
其中,q为炉膛压力比例积分微分控制器输出的阶跃变化幅度。
第三步:依据第二步计算出的K,T和τ整定炉膛压力比例积分微分控制器的参数:
Kc=1.2T/Kτ
Ti=2τ
Td=0.5τ
其中Kc,Ti,Td分别为比例积分微分控制器的比例参数,积分参数,微分参数。
(3)设计炉膛压力过程的预测PI-PID控制器,具体方法是:
针对设计的炉膛压力比例积分微分控制器和过程模型组成的基本控制回路建立该锅炉炉膛压力实时运行过程数据库,通过数据采集装置采集炉膛压力实时过程运行数据,依据炉膛压力实时过程运行数据建立预测PI-PID控制所需的预测模型,基于该预测模型设计相应的炉膛压力实时过程预测PI-PID控制器,具体步骤是:
第一步:将炉膛压力比例积分微分控制器停留在“自动操作”状态,操作炉膛压力比例积分微分控制器的输入使炉膛压力比例积分微分控制器的输入有个阶跃变化,由记录仪表记录炉膛压力实时过程的输出,将炉膛压力实时过程输出值y(k)的响应曲线转换成无量纲形式y*(k):
y*(k)=y(k)/y(∞)
其中,y(∞)是炉膛压力实时过程输出y(k)的稳态值。
第二步:选取2个计算点,y(k3)=0.39,y(k4)=0.63,依据以下计算公式计算炉膛压力预测PI-PID控制器所需要的参数K1,T1和τ1:
K1=y(∞)/q1
T1=2(k3-k4)
τ1=2k3-k4
其中,q1为炉膛压力比例积分微分控制器输入的阶跃变化幅度。
第三步:将第二步得到的参数转化为拉普拉斯形式的局部受控传递函数模型:
其中,y(s)表示当前时刻炉膛压力过程模型输出值的拉普拉斯变换,q1(s)表示炉膛压力过程模型的比例积分微分控制器输入的拉普拉斯变换。
λ1=T1
L1=τ1
第四步:依据第三步计算出的模型参数整定炉膛压力预测PI-PID控制器的参数,具体方法是:
①对该对象设计预测比例积分控制器。选期望的闭环传递函数模型为Gq2(s)
λ2为期望的闭环传递函数模型的时间常数,L2为期望的闭环传递函数模型的时滞,L2=L1;
②预测比例积分比例积分微分控制器的传递函数Gc1(s)可由下式来表示
③依据步骤②得到当前的预测比例积分比例积分微分控制器参数值u(s)
Claims (1)
1.燃煤锅炉炉膛压力***混合控制方法,其特征在于该方法包括以下步骤:
(1)利用燃煤锅炉炉膛压力实时过程数据建立过程模型,具体方法是:
首先建立燃煤锅炉炉膛压力实时运行数据库,通过数据采集装置采集N组实时过程运行数据,将采集的实时过程运行数据作为数据驱动的样本集合,表示为{xi,y(i)}i=1 N,i=1,2,…,N,其中xi表示第i组工艺参数的输入数据,y(i)表示第i组工艺参数的输出值;
然后以该炉膛压力实时过程运行数据集合为基础建立基于最小二乘法的离散差分方程形式的局部受控自回归滑动平均模型:
X=[y(k-1),…,y(k-n),u(k-d-1),…,u(k-d-m)]T
其中yL(k)表示当前时刻过程模型的工艺参数的输出值,x表示过程模型的工艺参数的过去时刻的输入和输出数据的集合,u(k)表示当前过程模型工艺参数对应的控制变量,k为当前的递推步数,Φ表示通过辨识得到的模型参数的集合,T表示矩阵的转置,n,m,d+1分别为对应实际过程的输出变量阶次、输入变量阶次、实际过程的时滞;
采用的辨识手段为:
(2)采用典型的响应曲线法设计炉膛压力过程模型的比例积分微分控制器,具体方法是:
a.将过程模型的比例积分微分控制器停留在手动操作状态,操作拨盘使其输出有阶跃变化,由记录仪表记录过程模型的输出值,将过程模型输出值yL(k)的响应曲线转换成无量纲形式yL *(k),具体是:
其中,yL(∞)是过程模型的比例积分微分控制器的输出有阶跃变化时的过程模型输出yL(k)的稳态值;
K=yL(∞)/q
T=2(k1-k2)
τ=2k1-k2
其中q为过程模型的比例积分微分控制器输出的阶跃变化幅度;
c.计算过程模型的比例积分微分控制器的参数,具体是:
Kc=1.2T/Kτ
Ti=2τ
Td=0.5τ
其中Kc为比例积分微分控制器的比例参数,Ti为比例积分微分控制器的积分参数,Td分别为比例积分微分控制器的微分参数;
(3)设计预测比例积分比例积分微分控制器,具体步骤是:
d.将过程模型的比例积分微分控制器停留在自动操作状态,操作拨盘使其输入有阶跃变化,由记录仪表记录实时过程的输出,将过程输出值y(k)的响应曲线转换成无量纲形式y*(k),具体是:y*(k)=y(k)/y(∞)
其中,y(∞)是过程模型的比例积分微分控制器的输入有阶跃变化时的过程模型输出y(k)的稳态值;
e.选取满足y(k3)=0.39,y(k4)=0.63的另两个计算点k3和k4,依据下式计算预测比例积分比例积分微分控制器所需要的参数K1,T1和τ1:
K1=y(∞)/q1
T1=2(k3-k4)
τ1=2k3-k4
其中q1为过程模型的比例积分微分控制器输入的阶跃变化幅度;
f.将步骤e得到的参数转化为拉普拉斯形式的局部受控传递函数模型:
其中s为拉普拉斯变换算子,λ1为局部受控传递函数模型的时间常数,L1为局部受控传递函数模型的时滞,y(s)表示当前时刻过程模型的输出值的拉普拉斯变换,q1(s)表示过程模型的比例积分微分控制器输入的拉普拉斯变换;
λ1=T1
L1=τ1
g.依据步骤f计算出的模型参数整定预测比例积分比例积分微分控制器的参数,具体方法是:
①对该对象设计预测比例积分控制器;选期望的闭环传递函数模型为Gq2(s)
λ2为期望的闭环传递函数模型的时间常数,L2为期望的闭环传递函数模型的时滞,L2=L1;
②预测比例积分比例积分微分控制器的传递函数Gc1(s)可由下式来表示
③依据步骤②得到当前的预测比例积分比例积分微分控制器参数值u(s)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101557929A CN101709863B (zh) | 2009-12-18 | 2009-12-18 | 燃煤锅炉炉膛压力***混合控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101557929A CN101709863B (zh) | 2009-12-18 | 2009-12-18 | 燃煤锅炉炉膛压力***混合控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101709863A CN101709863A (zh) | 2010-05-19 |
CN101709863B true CN101709863B (zh) | 2012-01-11 |
Family
ID=42402662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101557929A Expired - Fee Related CN101709863B (zh) | 2009-12-18 | 2009-12-18 | 燃煤锅炉炉膛压力***混合控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101709863B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102840571B (zh) * | 2012-09-20 | 2015-08-12 | 贵州电力试验研究院 | 基于子空间辨识的火力电站锅炉输出过热蒸汽的预报方法 |
CN102915011B (zh) * | 2012-09-26 | 2015-06-10 | 中国神华能源股份有限公司 | 实现分布式控制***对五个参数取中操作的方法及装置 |
CN103090410B (zh) * | 2013-02-01 | 2015-07-08 | 莱芜钢铁集团电子有限公司 | 加热炉助燃风压力控制方法、装置及*** |
CN104317321A (zh) * | 2014-09-23 | 2015-01-28 | 杭州电子科技大学 | 状态空间预测函数控制优化的焦化炉炉膛压力控制方法 |
CN104964263B (zh) * | 2015-06-23 | 2018-04-10 | 大唐韩城第二发电有限责任公司 | 一种锅炉主汽压力的控制方法 |
CN105159097A (zh) * | 2015-10-10 | 2015-12-16 | 杭州电子科技大学 | 炼油加热炉炉膛压力的多变量预测控制pid控制方法 |
CN107991886B (zh) * | 2017-12-28 | 2020-08-28 | 杭州电子科技大学 | 一种废塑料气化炼油炉膛压力的预测优化控制方法 |
CN112180876B (zh) * | 2020-10-19 | 2022-07-19 | 广东省特种设备检测研究院 | 基于大数据的燃气锅炉的节能控制方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1844739A (zh) * | 2006-04-11 | 2006-10-11 | 中控科技集团有限公司 | 一种循环流化床锅炉负荷串级燃烧控制***及方法 |
-
2009
- 2009-12-18 CN CN2009101557929A patent/CN101709863B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1844739A (zh) * | 2006-04-11 | 2006-10-11 | 中控科技集团有限公司 | 一种循环流化床锅炉负荷串级燃烧控制***及方法 |
Non-Patent Citations (1)
Title |
---|
陈云云,宁芳青.锅炉蒸汽温度控制***中广义预测控制的引用研究.《冶金自动化》.2005,394-396. * |
Also Published As
Publication number | Publication date |
---|---|
CN101709863A (zh) | 2010-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101709869B (zh) | 燃煤锅炉过热蒸汽温度***混合控制方法 | |
CN101709863B (zh) | 燃煤锅炉炉膛压力***混合控制方法 | |
CN100545772C (zh) | 一种燃煤锅炉***混合控制方法 | |
CN101709867B (zh) | 燃煤锅炉汽包水位***混合控制方法 | |
CN101286044B (zh) | 一种燃煤锅炉***蒸汽温度混合建模方法 | |
CN102401371B (zh) | 一种基于多变量预测控制的再热汽温优化控制方法 | |
CN103225799B (zh) | 一种火电厂主蒸汽温度控制方法 | |
CN102053562B (zh) | 一种裂解炉出口温度混合控制方法 | |
CN102608911B (zh) | 一种基于多参数预测的火电机组协调控制方法 | |
CN102156496B (zh) | 一种反应釜温度混合控制方法 | |
CN100596325C (zh) | 一种循环流化床锅炉负荷串级燃烧控制***及方法 | |
CN102494336A (zh) | 一种循环流化床锅炉燃烧过程多变量控制方法 | |
CN104699064B (zh) | 一种火电机组协调控制***的煤质自校正预测控制方法 | |
CN104281057A (zh) | 一种应用于变压器冷却***的复合pid模糊控制方法 | |
CN113031565B (zh) | 一种火电机组一次调频功率响应性能预测方法及*** | |
WO2022105356A1 (zh) | 带有增量式调节功能的核电机组控制棒调节方法及*** | |
CN105180136A (zh) | 基于分数阶pi动态矩阵的火电厂锅炉主蒸汽温度控制方法 | |
CN102736640A (zh) | 用于除氧器水位控制的***和方法 | |
CN105202519A (zh) | 供热机组调频调峰全工况协调控制方法 | |
CN202677208U (zh) | 用于除氧器水位控制的*** | |
CN110579968A (zh) | 一种超超临界机组深度调峰协调***预测控制策略 | |
CN110673482A (zh) | 一种基于神经网络预测的电站燃煤锅炉智能控制方法和*** | |
CN210951273U (zh) | 一种超临界机组除氧器液位自动控制*** | |
CN104196640A (zh) | 一种基于重型燃气轮机模型的解耦控制方法和*** | |
CN103336438A (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: 20120111 Termination date: 20211218 |