CN113091043B - 一种余热锅炉汽包水位全程自动控制的方法 - Google Patents
一种余热锅炉汽包水位全程自动控制的方法 Download PDFInfo
- Publication number
- CN113091043B CN113091043B CN202110231345.8A CN202110231345A CN113091043B CN 113091043 B CN113091043 B CN 113091043B CN 202110231345 A CN202110231345 A CN 202110231345A CN 113091043 B CN113091043 B CN 113091043B
- Authority
- CN
- China
- Prior art keywords
- lvl
- water level
- liquid level
- controller
- steam drum
- 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
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F22—STEAM GENERATION
- F22D—PREHEATING, OR ACCUMULATING PREHEATED, FEED-WATER FOR STEAM GENERATION; FEED-WATER SUPPLY FOR STEAM GENERATION; CONTROLLING WATER LEVEL FOR STEAM GENERATION; AUXILIARY DEVICES FOR PROMOTING WATER CIRCULATION WITHIN STEAM BOILERS
- F22D5/00—Controlling water feed or water level; Automatic water feeding or water-level regulators
- F22D5/26—Automatic feed-control systems
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E20/00—Combustion technologies with mitigation potential
- Y02E20/16—Combined cycle power plant [CCPP], or combined cycle gas turbine [CCGT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Thermal Sciences (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Control Of Steam Boilers And Waste-Gas Boilers (AREA)
Abstract
本发明公开了一种余热锅炉汽包水位全程自动控制的方法,涉及余热锅炉汽包水位控制技术,该方法步骤包括:步骤1:控制参数自适应模块将实际给水量输入至主控制器参数函数f1(x)和副控制器参数函数f2(x),监测运行工况,获取主、副控制器的控制参数;步骤2:前馈回路自适应模块根据PV‑SP值,生成7种液位模式,再根据不同的液位模式选择对应的前馈回路,前馈回路包括:FF1、FF2、FF3、FF4、FF5、FF6和FF7;步骤3:执行机构流量特性自适应模块实时计算现有执行机构流量特性的偏差,当偏差出现次数时间和次数满足要求时,触发流量特性曲线生成机制,并将生成的曲线更新至现有的特性曲线中。本发明的有益效果:汽包水位能全程投入自动。
Description
技术领域
本发明涉及余热锅炉汽包水位控制技术,具体的,涉及一种余热锅炉汽包水位全程自动控制的方法。
背景技术
现在燃气电厂几乎都采用燃气-蒸汽联合循环发电方式,以追求高发电效率。余热锅炉回收燃气轮机的余热,产生高压蒸汽,冲转汽轮机带动汽轮发电机组做做功。余热锅炉中一般有高压汽包、中压汽包、低压汽包水位,当其中一个水位过高和过低均会影响机组的安全,因此汽包水位是余热锅炉重要的控制参数之一。
目前余热锅炉的水位控制***在启机、停机、低负荷、变负荷阶段一般以人工操作为主,只有在机组负荷稳定的情况下采用自动控制。采用的控制技术大多基于三冲量控制、PID串级控制等,实际运行过程中,只有在***相对稳定的情况下实现自动控制。尤其是在出现虚假水位时,水位测量都会存在一定偏差。现有的控制方法难以实现余热锅炉汽包(汽包较小)水位的全程自动。
申请号为“CN201710457839.1”的专利公开了汽包锅炉启动阶段汽包水位控制方法及控制***,虽然该***可以控制汽包水位,但是该***只能控制汽包锅炉启动时的汽包水位,不能实现锅炉汽包水位全程自动调节。
发明内容
本发明的目的在于,针对上述问题,提出一种余热锅炉汽包水位全程自动控制的方法。
一种余热锅炉汽包水位全程自动控制的方法,包括如下步骤:
步骤1:控制参数自适应模块将实际给水量输入至主控制器参数函数f1(x)和副控制器参数函数f2(x),监测运行工况,获取主、副控制器的控制参数;
步骤2:前馈回路自适应模块根据PV-SP值,生成7种液位模式,再根据不同的液位模式选择对应的前馈回路,前馈回路包括:FF1、FF2、FF3、FF4、FF5、FF6和FF7;
步骤3:执行机构流量特性自适应模块实时计算现有执行机构流量特性的偏差,当偏差出现次数时间和次数满足要求时,触发流量特性曲线生成机制,并将生成的曲线更新至现有的特性曲线中。
所述步骤1包括如下子步骤:
步骤11:根据对象的动态特性,通过经验整定,输出同一控制器不同阶段的Kp,Ti控制参数组,
汽包水位主控制器控制参数组:
C1={(Kp1,Ti1),(Kp2,Ti2),…,(Kpn,Tin)},n∈R,
n代表控制器中参数的个数,(Kpn,Tin)代表不同阶段主控器比例和积分参数组;
汽包水位副控制器控制参数组:
C2={(Kp1,Ti1),(Kp2,Ti2),…,(Kpn,Tin)},n∈R,
n代表控制器中参数的个数,(Kpn,Tin)代表不同阶段主控器比例和积分参数组;
步骤12:利用二次函数进行平滑控制器参数组C1和C2中的Kp,Ti与给水流量之间的对应函数:
式中a2,a1,a0分别为待求系数,FL为给水流量,f(x)为在相应给水流量下的参数(Kp,Ti);
步骤13:利用最小二乘拟合公式,计算a2,a1,a0的数值,计算公式为:
Y=X*θ+E,
式中:X为n×k的矩阵,Y为n×1的列向量,θ为n×1的待求方程系数列向量,E为n×1的残差列向量,n为控制器参数的个数,k为待求方程系数的个数;
步骤14:令目标函数为:Q=ET*E,求解θ的最小二乘估计:
所述步骤2中的液位模式包括:
液位高危急:LVL_HHH;
液位高高:LVL_HH;
液位高:LVL_H;
液位正常:LVL_NORMAL;
液位低:LVL_L;
液位低低:LVL_LL;
液位低危急:LVL_LLL。
所述液位模式与前馈回路的关系为:
当触发液位高危急(LVL_HHH)时,前馈回路选择:
当触发液位高高(LVL_HH)时,前馈回路选择:
当触发液位高(LVL_H)时,前馈回路选择:
当触发液位正常(LVL_NORMAL)时,前馈回路选择:
当触发液位低(LVL_L)时,前馈回路选择:
当触发液位低低(LVL_LL)时,前馈回路选择:
当触发液位低危急(LVL_LLL)时,前馈回路选择:
上述公式中,Ki1,KRi,Ki2,K0i(i=1,2,…7)为相关系数,Ki1∈(-3,3),KRi∈(-300,300),Ki2∈(-10,10),K0i∈(-15,15);PV为汽包水位的测量值;SP为运行人员设定汽包水位的控制值。
所述步骤3包括如下子步骤:
步骤31:监测执行机构的指令YD和实际给水流量FL的正常运行数据,实时计算现有执行机构流量特性的偏差,当|YD-FL|的绝对值连续2分钟内大于2时,这样的情况出现3次,对现有的执行机构流量特性曲线进行修正;
步骤32:根据公式计算流量特性曲线y(x,w):
式中:x为实际给水流量,y为计算得到执行机构的指令,n为数据样本个数,wn代求多项式的系数;
误差函数为:
式中:N为x样本个数;
目标优化函数:
式中:N为x样本个数;
步骤33:根据步骤32中得到拟合的执行机构流量特性曲线y(x,w),选择x={0,10,20,30,40,50,55,60,75,80,85,87,90,91,93,95,97,100},分别将x输入至y(x,w)函数获取对应的y值;再将获得的(x,y)值导入值执行机构流量特性曲线模块F(x)中。
本发明的有益效果:控制参数自适应,余热锅炉启、停全过程只需要一套三冲量PID调节即可,减少逻辑冗余程度,提高逻辑运行效率;前馈回路自适应,提高在特殊工况下(如并泵,启、停机组,虚假水位等)的调节能力,确保汽包水位参数在“安全”范围内;给水执行机构流量特性自适应,提高执行机构的线性化,并减少运维工作量;汽包水位能全程投入自动,并在各种复杂工况下汽包水位均能控制在安全运行范围之内。
附图说明
图1为本发明的控制流程图。
具体实施方式
下面结合附图对本发明作进一步的说明。
如图1所示,一种余热锅炉汽包水位全程自动控制的方法,通过控制参数自适应模块、前馈回路自适应模块和执行机构流量特性自适应模块进行汽包水位控制。
控制参数自适应模块,根据生产工艺对控制品质的约束,如如衰减率φ,超调Mp等参数。依据对象的动态特性,通过经验整定,给出同一控制器不同阶段的Kp,Ti控制参数组,其中:
汽包水位主控制器控制参数组:
C1={(Kp1,Ti1),(Kp2,Ti2),…,(Kpn,Tin)},
n∈R,n代表控制器中参数的个数,(Kpn,Tin)代表不同阶段主控器比例和积分参数组。
汽包水位副控制器控制参数组:
C2={(Kp1,Ti1),(Kp2,Ti2),…,(Kpn,Tin)},
n∈R,n代表控制器中参数的个数,(Kpn,Tin)代表不同阶段主控器比例和积分参数组。
再利用二次函数进行平滑控制器参数组C1和C2中的Kp,Ti与给水流量之间的对应函数,即:
式中a2,a1,a0分别为待求系数,FL为给水流量,f(x)为在相应给水流量下的参数(Kp,Ti)。
再利用最小二乘拟合公式,计算a2,a1,a0数值,计算公式为:
Y=X*θ+E;
式中:X为n×k的矩阵,Y为n×1的列向量;Θ为n×1的待求方程系数列向量;E为n×1的残差列向量;而n为控制器参数的个数,k为待求方程系数的个数。
再令目标函数Q=ET*E,则可以得到求解θ的最小二乘估计:
通过上述方法将主控器参数组C1和副控制器C2转换为主控制器参数函数f1(x)和副控制器参数函数f2(x),当随着给水流量的不同,将得到不同的主、副控制器参数组,以适应机组不同的工况,取得满足生产工艺要求的控制品质。
前馈回路自适应模块,通过通过DCS实时采集汽包水位的信号,并三选、压力补偿等处理,作为汽包水位的测量值PV,运行人员设定汽包水位的控制值SP。根据PV-SP值的不同,生成7种液位模式。分别是液位高危急:LVL_HHH,液位高高LVL_HH,液位高LVL_H,液位正常:LVL_NORMAL,液位低:LVL_L,液位低低:LVL_LL,液位低危急:LVL_LLL。
当触发液位高危急LVL_HHH时,前馈回路选择:
当触发液位高高(LVL_HH)时,前馈回路选择:
当触发液位高(LVL_H)时,前馈回路选择:
当触发液位正常(LVL_NORMAL)时,前馈回路选择:
当触发液位低(LVL_L)时,前馈回路选择:
当触发液位低低(LVL_LL)时,前馈回路选择:
当触发液位低危急(LVL_LLL)时,前馈回路选择:
上述公式中,Ki1,KRi,Ki2,K0i(i=1,2,…7)为相关系数,Ki1∈(-3,3),KRi∈(-300,300),Ki2∈(-10,10),K0i∈(-15,15),需要根据不同的控制对象进行调整。
执行机构流量特性自适应模块,依据执行机构的指令YD和实际给水流量FL的正常运行数据,实时计算现有执行机构流量特性的偏差,当|YD-FL|的绝对值连续2分钟内大于2时,这样的情况出现3次,需要对现有的执行机构流量特性曲线进行修正。
流量特性曲线计算如下:
依据执行机构的指令(YD)和实际给水流量(FL)正常运行的数据,产生特性曲线,通过以下多项式拟合其曲线:
式中x为实际给水流量,y为计算得到执行机构的指令,n为数据样本个数,wn代求多项式的系数。
误差函数为:
式中:N为x样本个数;
目标优化函数:
式中:N为x样本个数。
通过上述算法得到拟合的执行机构流量特性曲线y(x,w)。选择x={0,10,20,30,40,50,55,60,75,80,85,87,90,91,93,95,97,100},分别将x输入至y(x,w)函数获取对应的y值。
最后,将获得的(x,y)值导入值执行机构流量特性曲线模块中。
Claims (5)
1.一种余热锅炉汽包水位全程自动控制的方法,其特征在于,包括如下步骤:
步骤1:控制参数自适应模块将实际给水量输入至主控制器参数函数f1(x)和副控制器参数函数f2(x),监测运行工况,获取主、副控制器的控制参数;
步骤2:前馈回路自适应模块根据PV-SP值,生成7种液位模式,再根据不同的液位模式选择对应的前馈回路,PV是汽包水位测量值,SP是设定的汽包水位控制值;
步骤3:执行机构流量特性自适应模块实时计算现有执行机构流量特性的偏差,当偏差出现次数时间和次数满足要求时,触发流量特性曲线生成机制,并将生成的曲线更新至现有的特性曲线中。
2.如权利要求1所述的一种余热锅炉汽包水位全程自动控制的方法,其特征在于,所述步骤1包括如下子步骤:
步骤11:根据对象的动态特性,通过经验整定,输出同一控制器不同阶段的Kp,Ti控制参数组,
汽包水位主控制器控制参数组:
C1={(Kp1,Ti1),(Kp2,Ti2),…,(Kpn,Tin)},n∈R,
n代表控制器中参数的个数;(Kpn,Tin)代表不同阶段主控器比例和积分参数组;
汽包水位副控制器控制参数组:
C2={(Kp1,Ti1),(Kp2,Ti2),…,(Kpn,Tin)},n∈R,
n代表控制器中参数的个数;(Kpn,Tin)代表不同阶段主控器比例和积分参数组;
步骤12:利用二次函数进行平滑控制器参数组C1和C2中的Kp,Ti与给水流量之间的对应函数:
式中a2,a1,a0分别为待求系数,FL为给水流量,f(x)为在相应给水流量下的参数(Kp,Ti);
步骤13:利用最小二乘拟合公式,计算a2,a1,a0的数值,计算公式为:
Y=X*θ+E,
式中:X为n×k的矩阵,Y为n×1的列向量,θ为n×1的待求方程系数列向量,E为n×1的残差列向量,n为控制器参数的个数,k为待求方程系数的个数;
步骤14:令目标函数为:Q=ET*E,求解θ的最小二乘估计:
3.如权利要求1所述的一种余热锅炉汽包水位全程自动控制的方法,其特征在于,所述步骤2中的液位模式包括:
液位高危急:LVL_HHH;
液位高高:LVL_HH;
液位高:LVL_H;
液位正常:LVL_NORMAL;
液位低:LVL_L;
液位低低:LVL_LL;
液位低危急:LVL_LLL。
4.如权利要求3所述的一种余热锅炉汽包水位全程自动控制的方法,其特征在于,所述液位模式与前馈回路的关系为:
当触发液位高危急(LVL_HHH)时,前馈回路选择:
当触发液位高高(LVL_HH)时,前馈回路选择:
当触发液位高(LVL_H)时,前馈回路选择:
当触发液位正常(LVL_NORMAL)时,前馈回路选择:
当触发液位低(LVL_L)时,前馈回路选择:
当触发液位低低(LVL_LL)时,前馈回路选择:
当触发液位低危急(LVL_LLL)时,前馈回路选择:
上述公式中,Ki1,KRi,Ki2,K0i(i=1,2,…7)为相关系数,Ki1∈(-3,3),KRi∈(-300,300),Ki2∈(-10,10),K0i∈(-15,15);PV为汽包水位的测量值;SP为运行人员设定汽包水位的控制值。
5.如权利要求1所述的一种余热锅炉汽包水位全程自动控制的方法,其特征在于,所述步骤3包括如下子步骤:
步骤31:监测执行机构的指令YD和实际给水流量FL的正常运行数据,实时计算现有执行机构流量特性的偏差,当|YD-FL|的绝对值连续2分钟内大于2时,这样的情况出现3次,对现有的执行机构流量特性曲线进行修正;
步骤32:根据公式计算流量特性曲线y(x,w):
式中:x为实际给水流量,y为计算得到执行机构的指令,n为数据样本个数,wn为代求多项式的系数;
误差函数为:
式中:N为x样本个数;
目标优化函数:
式中:N为x样本个数;
步骤33:根据步骤32中得到拟合的执行机构流量特性曲线y(x,w),选择x={0,10,20,30,40,50,55,60,75,80,85,87,90,91,93,95,97,100},分别将x输入至y(x,w)函数获取对应的y值;再将获得的(x,y)值导入值执行机构流量特性曲线模块F(x)中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110231345.8A CN113091043B (zh) | 2021-03-02 | 2021-03-02 | 一种余热锅炉汽包水位全程自动控制的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110231345.8A CN113091043B (zh) | 2021-03-02 | 2021-03-02 | 一种余热锅炉汽包水位全程自动控制的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113091043A CN113091043A (zh) | 2021-07-09 |
CN113091043B true CN113091043B (zh) | 2023-03-21 |
Family
ID=76666481
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110231345.8A Active CN113091043B (zh) | 2021-03-02 | 2021-03-02 | 一种余热锅炉汽包水位全程自动控制的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113091043B (zh) |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030075158A1 (en) * | 2000-02-09 | 2003-04-24 | Leonhard Milos | Method and device for a mass flow determination via a control valve and for determining a modeled induction pipe pressure |
JP4414269B2 (ja) * | 2004-04-23 | 2010-02-10 | シーケーディ株式会社 | 電磁比例弁の流量特性調整機構及びそれを用いた流量特性調整方法 |
RU2426943C2 (ru) * | 2009-09-08 | 2011-08-20 | Государственное образовательное учреждение высшего профессионального образования "Ивановский государственный энергетический университет имени В.И. Ленина" (ИГЭУ) | Система регулирования питания водой барабанного парогенератора |
CN102997222B (zh) * | 2012-11-27 | 2015-04-22 | 华北电力大学 | 一种锅炉汽包水位自适应pid控制方法 |
CN103558757B (zh) * | 2013-11-07 | 2016-09-21 | 中冶南方工程技术有限公司 | 热电锅炉汽包液位控制方法 |
US9476584B2 (en) * | 2013-12-12 | 2016-10-25 | General Electric Company | Controlling boiler drum level |
CN105159337B (zh) * | 2015-06-17 | 2017-11-17 | 国网天津市电力公司 | 电站锅炉汽包水位调节***快速投入自动方法 |
UA106702U (uk) * | 2015-08-12 | 2016-05-10 | Національний Аерокосмічний Університет Ім. М.Є. Жуковського "Харківський Авіаційний Інститут" | Адаптивна система автоматичного керування рівнем пароводяної суміші в котлі барабанного типу |
CN107016239B (zh) * | 2017-03-30 | 2020-01-31 | 贵州电网有限责任公司电力科学研究院 | 一种汽轮机阀门流量特性分析方法 |
CN107219758A (zh) * | 2017-07-03 | 2017-09-29 | 安徽马钢自动化信息技术有限公司 | 一种汽包液位调节方法 |
-
2021
- 2021-03-02 CN CN202110231345.8A patent/CN113091043B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113091043A (zh) | 2021-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5916043B2 (ja) | 湿分分離再加熱器を制御する方法及び装置 | |
WO2019165839A1 (zh) | 一种锅炉过热蒸汽的温度控制方法、装置和*** | |
CN102374518B (zh) | 使用动态矩阵控制的蒸汽温度控制 | |
CN104714526B (zh) | 基于凝结水节流调节预估的负荷控制***及方法 | |
JP5108644B2 (ja) | ボイラ制御装置ならびにボイラ制御方法 | |
CN110824926B (zh) | 一种火电机组深度调峰一次调频控制方法 | |
CN107780982B (zh) | 一种在线的间接空冷高背压供热机组背压控制***及方法 | |
US10126810B2 (en) | Method for controlling power generation unit to desired output as specified by load demand signal by using modified control signal | |
CN111637444B (zh) | 一种基于q学习的核电蒸汽发生器水位控制方法 | |
CN114089795B (zh) | 一种基于事件触发的模糊神经网络温度控制***及方法 | |
JP6608324B2 (ja) | 原子力プラントの出力制御装置及び方法 | |
US10240775B2 (en) | Multi-objective steam temperature control | |
EP1533482B1 (en) | Rapid power producing system and method for steam turbine | |
CN110376895A (zh) | 一种基于分层受限预测控制的火电机组协调控制方法 | |
CN112051810A (zh) | 一种供热机组的多能源分配和调度*** | |
KR100968796B1 (ko) | 보일러 및 터빈 응답지연 보상에 의한 발전소 출력 제어 안정도 향상 제어 방법 | |
CN110985139A (zh) | 一种基于应用备用汽源的小机进汽控制装置及方法 | |
CN113091043B (zh) | 一种余热锅炉汽包水位全程自动控制的方法 | |
CN104685426B (zh) | 用于热电厂中的控制和故障分析的优化的方法 | |
JP3666036B2 (ja) | 火力発電プラント起動制御システム及び起動制御方法 | |
CN111486441B (zh) | 一种锅炉过热蒸汽温度的超驰回路控制方法 | |
CN114562713A (zh) | 发电锅炉主蒸汽温度控制方法及*** | |
Emoto et al. | Integrated large-scale multivariable control and real-time optimization of a power plant | |
Shu et al. | Load regulation of the heat supply unit assisted by the energy storage in thermal networks via a reference model-based controller | |
Li et al. | Research and Treatment of Periodic Fluctuation of Steam Feed Pump Speed Regarding 900MW Nuclear Power Plant |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |