CN111881554B - 一种锅炉随气温变化的优化控制方法 - Google Patents

一种锅炉随气温变化的优化控制方法 Download PDF

Info

Publication number
CN111881554B
CN111881554B CN202010603065.0A CN202010603065A CN111881554B CN 111881554 B CN111881554 B CN 111881554B CN 202010603065 A CN202010603065 A CN 202010603065A CN 111881554 B CN111881554 B CN 111881554B
Authority
CN
China
Prior art keywords
load
interval
air temperature
dimensional
data
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
Application number
CN202010603065.0A
Other languages
English (en)
Other versions
CN111881554A (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.)
Jiangsu Hanguang Intelligent Technology Co ltd
China University of Mining and Technology CUMT
Northeast Electric Power University
Original Assignee
Jiangsu Hanguang Intelligent Technology Co ltd
China University of Mining and Technology CUMT
Northeast Dianli 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 Jiangsu Hanguang Intelligent Technology Co ltd, China University of Mining and Technology CUMT, Northeast Dianli University filed Critical Jiangsu Hanguang Intelligent Technology Co ltd
Priority to CN202010603065.0A priority Critical patent/CN111881554B/zh
Publication of CN111881554A publication Critical patent/CN111881554A/zh
Application granted granted Critical
Publication of CN111881554B publication Critical patent/CN111881554B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24Querying
    • G06F16/245Query processing
    • G06F16/2458Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2462Approximate or statistical queries
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00Energy efficient computing, e.g. low power processors, power management or thermal management

Abstract

本发明提出一种锅炉随气温变化的优化控制方法,其特点是:包括如下步骤:1)建立按负荷、气温分布的优化控制规律:首先确定提取参数类型、设置数据包时间间隔,其次划分二维区间,计算所有二维区间各参数均值,然后在所有二维区间内筛选同时优于该区间的锅炉效率均值和NOx排放量均值的数据包,得到按负荷、气温分布的优化控制规律;2)生成优化控制指令:在锅炉变负荷调控中,按目标负荷、送风机入口风温匹配对应负荷、气温下的优化控制规律,最后根据优化控制规律与DCS控制信号生成修正后的优化指令进行精准调控,以实现锅炉的优化控制。本发明考虑气温对锅炉优化控制的影响,控制方法进一步贴合实际运行工况,具有快速响应、节能减排的优点。

Description

一种锅炉随气温变化的优化控制方法
技术领域
本发明属于锅炉运行优化控制技术领域,具体涉及一种锅炉随气温变化的优化控制方法,适用于燃煤锅炉变负荷时的闭环控制。
背景技术
实际运行的锅炉,其工况总是处于不断变化之中。随着电力系统的发展、电网容量的扩大以及新能源占比的持续增加,以及新能源发电固有的波动性,造成电网的峰谷负荷差进一步拉大,这要求各容量机组都需要参与电网调峰。而就目前国内电网组成结构而言,主要靠燃煤火电机组来承担调峰任务,许多无法满足调峰要求的燃煤电厂,濒临停机困局,因此亟待解决燃煤火电机组锅炉变负荷运行的灵活性问题,但是,当前对负荷的优化控制依然停留在非常落后的总量参数控制方法上,导致锅炉变负荷调控存在以下问题:
1.负荷响应速度不达标
燃煤机组增减负荷日趋频繁,而锅炉系统本身是一个具有大滞后特性的对象,这使得传统“水、煤、风”反馈调控具有延迟叠加现象,导致频繁增减负荷的响应时间无法满足当前电网要求。
2.响应负荷要求时,机组节能减排的潜力未被挖掘
燃煤发电机组调节的首要任务是快速响应负荷要求,而锅炉在变负荷调节中无法保持高发电效率和低NOx排放,即发电效率和NOx排放的调节具有滞后性。寻找一种可同时满足负荷要求又能提高锅炉效率,减少NOx排放的优化控制方法不失一种好的选择。
3.预测模型灵活性和适应性差
锅炉实际运行工况复杂多样,建立模型应贴合实际运行工况,未按实际运行工况划分的预测模型,其模型适用性差,导致无法灵活应对越来越大的电网峰谷负荷差。
现有技术对以上问题进行的研究,并给出多种解决方案。如亚临界循环流化床机组基于机理的燃机系统模型、基于深度信任网络的燃烧效率和NOx排放模型、描述锅炉热效率和NOx分布的物理模型、应用解析建模和人工智能相结合的锅炉混合模型等。
也有针对火电机组运行优化过程中存在工况变化大和目标参数最优值冲突的问题,提出基于工况划分的火电机组运行多目标优化方法,并验证了其有效性。所以分工况建立计算模型,使预测工况与实际工况一一对应,可达到锅炉机组控制精度的实质性提高。
目前,尚无有关针对锅炉根据气温(季节)变化的优化控制方法的相关报导。一年四季气温变化很大,我国北方城市冬季最低气温可达-35℃,而夏季最高气温可达35℃,一年四季最大温差可达70℃,南方城市温差相对较小但也有40℃,而一天之内气温变化可达一年气温变化幅度的1/4。因此,锅炉的优化控制必须考虑送风温度的影响,本发明在负荷工况划分的基础上,开展多维度工况划分,提供一种锅炉随气温(季节)变化的优化控制方法,对锅炉的节能减排具有实际的科学意义。
发明内容
本发明要解决的技术问题是:克服现有技术的缺点和不足,提供一种节能减排的锅炉随负荷、气温变化的优化控制方法,节能体现在经济性高,减排体现在NOx排放低。将机组DCS历史数据按负荷与当地气温(送风机入口风温)划分二维工况区间;对所有二维区间进行数据挖掘,得到按负荷、气温分布的优化控制规律;依据对应工况区间的控制规律开展精准调控,以此达到锅炉的优化控制和节能减排。
特别说明一下,本发明仅仅以DCS记录的实际负荷/总燃料量作为经济性评价的依据,没有考虑煤质变化的影响,得到的结果将是在实际煤质波动条件下的平均煤质下的经济性和排放水平的相对变化,并不是追求绝对准确的经济性评价指标。
本发明解决技术问题的技术方案是:一种锅炉随气温变化的优化控制方法,其特征在于:包括建立按负荷、气温分布的优化控制规律和生成优化控制指令,具体步骤如下:
1)建立按负荷、气温分布的优于优化目标的优化控制规律
①获取机组历史运行数据
a确定从DCS历史数据库中提取的参数类型xk,参数类型包括:实际负荷、总燃料量、各层给煤机给煤量、各层磨煤机一次风量、二次风总量、各层二次风门开度、各层给煤机转速、总风量、给水流量、主蒸汽压力、调节级压力、烟气含氧量、送风机入口风温、NOx浓度和过热器壁温,其中参数类别以k划分,k=1,2,…,K个;
b设定数据的时间间隔值;
c根据步骤a中列出的参数和步骤b设定的数据时间间隔值,以同一时刻所述参数的数据单元作为数据包,提取燃煤机组过去一年DCS历史运行数据;
②按负荷、气温划分二维区间,将历史数据按二维区间进行分类
a根据全年负荷主要运行区段,确定负荷优化区间:设机组在主要运行区段内的最高运行负荷为Smax,最低运行负荷为Smin,则负荷优化区间为[Smin,Smax];
b确定负荷区间数:设负荷划分间隔为LS,i代表划分的负荷区间数,i=1,2,…,m个,则负荷区间数m通过公式(1)获得;
m=(Smax-Smin)/LS (1)
c根据全年送风机入口风温变化,确定气温优化区间:设全年送风机入口最高风温为Tmax,最低风温为Tmin,则气温优化区间为[Tmin,Tmax];
d确定气温区间数:设气温划分间隔为LT,j代表划分的气温区间数,j=1,2,…,n个,则气温区间数n通过公式(2)获得:
n=(Tmax-Tmin)/LT (2)
e由步骤a~d获得第i负荷段、第j气温段的二维区间为公式(3):
{[Smin+(i-1)×LS,Smin+i×LS],[Tmin+(j-1)×LT,Tmin+j×LT]} (3)
③计算各二维区间的各类参数均值、锅炉效率均值和NOx排放均值
a将步骤①获取的机组历史运行数据的数据包按照步骤②划分的二维区间进行分类,超出负荷、气温二维区间的数据包予以舍弃;
b统计步骤a中各二维区间内数据包数量l,l=1,2,…,Lij个;则式(3)所表述的二维区间内的运行数据表示为:xi,j,k,l,所述二维区间的各类参数均值为公式(4)所示:
Figure GDA0003849087200000031
式中:
Figure GDA0003849087200000032
表示第i负荷段、第j气温段内第k个参数的均值;
c定义锅炉效率均值=实际负荷均值/总燃料量均值,锅炉效率均值为公式(5)所示:
Figure GDA0003849087200000033
式中:
Figure GDA0003849087200000034
Figure GDA0003849087200000035
为该二维区间的总燃料量均值,
Figure GDA0003849087200000036
Figure GDA0003849087200000037
为该二维区间的实际负荷均值,其中
Figure GDA0003849087200000038
d计算NOx排放均值
Figure GDA0003849087200000039
Figure GDA00038490872000000310
e二维区间的各类参数出现无效数据的处理;
④根据优化目标筛选二维区间内的数据包
在第i负荷段、第j气温段组成的二维区间内,以优于锅炉效率均值
Figure GDA00038490872000000311
和NOx排放均值
Figure GDA00038490872000000312
为优化目标筛选数据包,筛选条件为公式(6)所示:
Figure GDA00038490872000000313
式中:xi,j,E,l和xi,j,NOx,l表示第i负荷段、第j气温段、第l个数据包的锅炉效率和NOx排放,其中
Figure GDA00038490872000000314
⑤各二维区间优于锅炉效率均值和NOx排放均值的数据包的处理
a符合优化目标的数据包经筛选后保存在原二维区间内,并统计该区间内符合优化目标的数据包数量l′,其中l′=1,2,…,Li,j′,不符合优化目标的数据包予以淘汰;
b优选后的二维区间的数据包个数少于限定值的处理;
⑥从各二维区间优于优化目标的数据集合中计算该区间各类参数均值
对于符合优化目标的数据包再次计算该二维区间各类参数均值为公式(7)所示:
Figure GDA0003849087200000041
式中:
Figure GDA0003849087200000042
表示筛选后的第i负荷段、第j气温段、第k个参数的均值;x′i,j,k,l′表示筛选后第i负荷段、第j气温段、第k个参数的第l′个数据,由此得到按负荷、气温分布的优于优化目标的优化控制规律;
⑦二维区间内异常数据的处理
对于二维区间内的异常数据,依据相邻工况运行数据的渐变性,采用正则化矩阵对突变峰值进行均值滤波;
2)生成优化控制指令
①由AGC负荷指令生成目标负荷;
②按目标负荷、送风机入口风温匹配对应负荷、气温下的优化控制规律;
③若目标负荷、送风机入口风温与对应负荷、气温下的优化控制规律存在偏差,进行负荷、气温的二维差值校正,假设校正差值为Δi,j,k,则经差值校正后的优化控制信号为
Figure GDA0003849087200000043
④计算优化控制信号
Figure GDA0003849087200000044
与对应工况下DCS实际控制信号Sin,k的绝对差值
Figure GDA0003849087200000045
为公式(8)所示:
Figure GDA0003849087200000046
⑤依据优化控制指令的调节幅度,生成修正后的优化指令
若差值
Figure GDA0003849087200000047
在DCS实际控制信号Sin,k的设定幅度值以内,即
Figure GDA0003849087200000048
Figure GDA0003849087200000049
则生成修正后的优化指令
Figure GDA00038490872000000410
若差值
Figure GDA00038490872000000411
超出或等于DCS实际控制信号Sin,k的设定幅度值,即
Figure GDA00038490872000000412
则生成修正后的优化指令Sin,k·(1±设定幅度值)。
对于步骤1)的①所述的设定数据的时间间隔值为1~3分钟。
对于步骤1)的③二维区间的各类参数出现无效数据的处理:在所划分的不同负荷、不同气温的工况下,并不能保证二维区间的各类参数全部为有效运行数据,此时,将采用同一负荷段内其他二维区间的同类参数的平均值予以替代。
对于步骤1)的⑤优选后的二维区间的数据包个数少于限定值的处理:在数据包个数少于限定值的情况,认定该二维区间所含有效数据量不全,淘汰该区间已优选的数据包,并采用同一负荷段内其他二维区间的同类参数的平均值予以替代。
所述步骤1)的⑤的限定值为500。
对于所述步骤1)的⑦二维区间内异常数据的处理,依据相邻工况运行数据的渐变性,采用正则化矩阵对突变峰值进行均值替代,具体如下:
a)利用步骤②中按负荷、气温分布的二维区间建立二维数组A,数组元素表示为Ai,j,其中负荷段i=1,2,…,m,气温段j=1,2,…,n,按数组元素位于数组边界、四个顶角和内部共三种位置分类处理;
b)若区间位于数组边界,利用相邻三个区间的数据求均值,如公式(9)所示:
A1,2=(A1,1+A2,2+A1,3)/3 (9)
c)若区间位于数组顶角,利用相邻两个区间的数据求均值,如公式(10)所示:
Am,n=(Am-1,n+Am,n-1)/2 (10)
d)若区间位于数组内部,利用相邻四个区间的数据求均值,如公式(11)所示:
Ai,j=(Ai-1,j+Ai,j-1+Ai+1,j+Ai,j+1)/4 (11)
对于步骤2)的⑤所述的设定幅度值为5%。
本发明的一种节能减排的锅炉随气温变化的优化控制方法的优点体现在:
一是考虑气温(送风机入口风温)对锅炉优化控制的影响,将历史运行数据按负荷、气温进行精细划分,优化方案最大限度贴合实际运行工况,保证调控质量;
二是在精准调控的基础上,筛选二维工况区间内优于锅炉效率均值(实际负荷均值/总燃料量均值)和NOx排放均值的控制规律,以达到节能减排的目的;
三是优化方案安全可靠,调节幅度控制在原有DCS控制信号的±5%以内,保证机组稳定运行。
附图说明
图1为本发明的流程图;
图2为本发明的负荷、气温二维区间分布图;
图3为本发明的优化控制指令逻辑图
图4为实施例2的燃煤机组送风机入口风温全年变化曲线;
图5为实施例2的优化燃料量随负荷、气温的变化;
图6为实施例2的优化进风量随负荷、气温的变化。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
参见图1-3,实施例1,本实施例一种锅炉随气温变化的优化控制方法,包括建立按负荷、气温分布的优化控制规律和生成优化控制指令,具体步骤如下:
1)建立按负荷、气温分布的优于优化目标的优化控制规律
①获取机组历史运行数据
a确定从DCS历史数据库中提取的参数类型xk,参数类型包括:实际负荷、总燃料量、各层给煤机给煤量、各层磨煤机一次风量、二次风总量、各层二次风门开度、各层给煤机转速、总风量、给水流量、主蒸汽压力、调节级压力、烟气含氧量、送风机入口风温、NOx浓度和过热器壁温,其中参数类别以k划分,k=1,2,…,K个;
b设定数据的时间间隔值为1分钟;
c根据步骤a中列出的参数,按照步骤b设定的数据时间间隔,以同一时刻所述参数的数据单元作为数据包,提取燃煤机组过去一年DCS历史运行数据;
②按负荷、气温划分二维区间,将历史数据按二维区间进行分类
a根据全年负荷主要运行区段,确定负荷优化区间:设机组在主要运行区段内的最高运行负荷为Smax,最低运行负荷为Smin,则负荷优化区间为[Smin,Smax];
b确定负荷区间数:设负荷划分间隔为LS,i代表划分的负荷区间数,i=1,2,…,m个,则负荷区间数m通过公式(1)获得;
m=(Smax-Smin)/LS (1)
c根据全年送风机入口风温变化,确定气温优化区间:设全年送风机入口最高风温为Tmax,最低风温为Tmin,则气温优化区间为[Tmin,Tmax];
d确定气温区间数:设气温划分间隔为LT,j代表划分的气温区间数,j=1,2,…,n个,则气温区间数n通过公式(2)获得:
n=(Tmax-Tmin)/LT (2)
e由步骤a~d获得第i负荷段、第j气温段的二维区间为公式(3)
{[Smin+(i-1)×LS,Smin+i×LS],[Tmin+(j-1)×LT,Tmin+j×LT]} (3)
③计算各二维区间的各类参数均值、锅炉效率均值和NOx排放均值
a将步骤①获取的机组历史运行数据的数据包按照步骤②划分的二维区间进行分类,超出负荷、气温二维区间的数据包予以舍弃;
b统计步骤a中各二维区间内数据包数量l,l=1,2,…,Li,j个;则式(3)所表述的二维区间内的运行数据表示为:xi,j,k,l,所述二维区间的各类参数均值为公式(4)所示:
Figure GDA0003849087200000071
式中:
Figure GDA0003849087200000072
表示第i负荷段、第j气温段内第k个参数的均值;
c定义锅炉效率均值=实际负荷均值/总燃料量均值,锅炉效率均值为公式(5)所示:
Figure GDA0003849087200000073
式中:
Figure GDA0003849087200000074
Figure GDA0003849087200000075
为该二维区间的总燃料量均值,
Figure GDA0003849087200000076
Figure GDA0003849087200000077
为该二维区间的实际负荷均值,其中
Figure GDA0003849087200000078
d计算NOx排放均值
Figure GDA0003849087200000079
Figure GDA00038490872000000710
e二维区间的各类参数出现无效数据的处理
在所划分的不同负荷、不同气温的工况下,并不能保证二维区间的各类参数全部都出现有效运行数据,此时,将采用同一负荷段内其他二维区间的同类参数的平均值予以替代。
④根据优化目标筛选二维区间内的数据包
在第i负荷段、第j气温段组成的二维区间内,通过优于锅炉效率均值
Figure GDA00038490872000000711
和NOx排放均值
Figure GDA00038490872000000712
来优选数据包,筛选条件为公式(6)所示:
Figure GDA00038490872000000713
式中:xi,j,E,l和xi,j,NOx,l表示第i负荷段、第j气温段、第l个数据包的锅炉效率和NOx排放,其中
Figure GDA00038490872000000714
⑤各二维区间优于锅炉效率均值和NOx排放均值的数据包的处理
a符合优化目标的数据包经筛选后保存在原二维区间内,并统计符合优化目标的数据包数量l′,其中l′=1,2,…,Li,j′,不符合优化目标的数据包予以淘汰;
b优选后的二维区间的数据包个数少于限定值的处理
在数据包个数少于限定值为500的情况,认定该二维区间所含有效数据量不全,淘汰该区间已优选的数据包,并采用同一负荷段内其他二维区间的同类参数的平均值予以替代。
⑥从各二维区间优于优化目标的数据集合中计算该区间各类参数均值
对于符合优化目标的数据包再次计算该二维区间各类参数均值为公式(7)所示:
Figure GDA0003849087200000081
式中:
Figure GDA0003849087200000082
表示筛选后的第i负荷段、第j气温段、第k个参数的均值;x′i,j,k,l′表示筛选后第i负荷段、第j气温段、第k个参数的第l′个数据,由此得到按负荷、气温分布的优于优化目标的优化控制规律;
⑦二维区间内异常数据的处理
对于二维区间内的异常数据,依据相邻工况运行数据的渐变性,采用正则化矩阵对突变峰值进行均值滤波;
2)生成优化控制指令
①由AGC负荷指令生成目标负荷;
②按目标负荷、送风机入口风温匹配对应负荷、气温下的优化控制规律;
③若目标负荷、送风机入口风温与对应负荷、气温下的优化控制规律存在偏差,进行负荷、气温的二维差值校正,假设校正差值为Δi,j,k,则经差值校正后的优化控制信号为
Figure GDA0003849087200000083
④计算优化控制信号
Figure GDA0003849087200000084
与对应工况下DCS实际控制信号Sin,k的绝对差值
Figure GDA0003849087200000085
为公式(8)所示:
Figure GDA0003849087200000086
⑤依据优化控制指令的调节幅度,生成修正后的优化指令
若差值
Figure GDA0003849087200000087
在DCS实际控制信号Sin,k的5%以内,即
Figure GDA0003849087200000088
由图3所示此时S3=0,信号指令S1有效,修正后的优化指令即为
Figure GDA0003849087200000089
若差值
Figure GDA00038490872000000810
超出或等于DCS实际控制信号Sin,k的5%,即
Figure GDA00038490872000000811
此时S3=1,信号指令S2有效,修正后的优化指令即为Sin,k·(1±5%),此举以防止调控幅度过大,影响机组平稳运行。通过以上信号修正方法,最终达到节能减排的优化调控目标;
对于二维区间内的异常数据,依据相邻工况运行数据的渐变性,采用正则化矩阵对突变峰值进行均值滤波,具体如下:
a)利用步骤②中按负荷、气温分布的二维区间建立二维数组A,数组元素表示为Ai,j,其中负荷段i=1,2,…,m,气温段j=1,2,…,n,按数组元素位于数组边界、四个顶角和内部共三种位置分类处理;
b)若区间位于数组边界,利用相邻三个区间的数据求均值,如公式(9)所示:
A1,2=(A1,1+A2,2+A1,3)/3 (9)
c)若区间位于数组顶角,利用相邻两个区间的数据求均值,如公式(10)所示:
Am,n=(Am-1,n+Am,n-1)/2 (10)
d)若区间位于数组内部,利用相邻四个区间的数据求均值,如公式(11)所示:
Ai,j=(Ai-1,j+Ai,j-1+Ai+1,j+Ai,j+1)/4 (11)
参见图1-图6,实施例2,本实施例以我国南方某电厂6号600MW燃煤发电机组为研究对象,首先确定提取的参数类型如表1所示,其中包括:实际负荷、总燃料量、各层给煤机给煤量、各层磨煤机一次风量、二次风总量、各层二次风门开度、各层给煤机转速、总风量、给水流量、主蒸汽压力、调节级压力、烟气含氧量、送风机入口风温、NOx浓度、过热器壁温等58类参数,即K=58,数据包的提取时间间隔设置为1分钟;
表1提取DCS参数列表
Figure GDA0003849087200000091
Figure GDA0003849087200000101
根据以上设置规则,提取DCS历史数据,分析其送风机入口风温全年温度变化曲线,如图4所示。其温度变化范围为-6~+38℃,为保证区间内历史数据的数据包个数满足需求,温度区间取[-4,36],气温间隔取LT=2℃,则气温段个数n=20;分析机组负荷主要运行区间为[300~600],负荷间隔取LS=10MW,则负荷段个数m=30。
本实施例按照实施例1的的流程图经建立按负荷、气温分布的优化控制规律及输出优化控制指令后,分析优化控制效果,提取燃料量和送风量的优化控制指令。图5为优化燃料量随负荷、气温的变化,从图中可以看出相同负荷下,优化燃料量随气温的升高略有下降,显示机组夏季的运行经济性略高于冬季。图6为优化进风量随负荷、气温的变化,从图中可以看出进风量随气温变化非常显著,中高气温段,随着气温的上升,锅炉优化送风量明显增加;中低气温段,存在一个居中的温度点,优化送风量最低;高于或者低于这个居中温度,优化送风量不同程度增加。
本发明的实施例并非穷举,本领域技术人员不经过创造性劳动的简单复制和改进,仍属于本发明权利保护的范围。

Claims (7)

1.一种锅炉随气温变化的优化控制方法,其特征在于:包括建立按负荷、气温分布的优化控制规律和生成优化控制指令,具体步骤如下:
1)建立按负荷、气温分布的优于优化目标的优化控制规律
①获取机组历史运行数据
a确定从DCS历史数据库中提取的参数类型xk,参数类型包括:实际负荷、总燃料量、各层给煤机给煤量、各层磨煤机一次风量、二次风总量、各层二次风门开度、各层给煤机转速、总风量、给水流量、主蒸汽压力、调节级压力、烟气含氧量、送风机入口风温、NOx浓度和过热器壁温,其中参数类别以k划分,k=1,2,...,K个;
b设定数据的时间间隔值;
c根据步骤a中列出的参数和步骤b设定的数据时间间隔值,以同一时刻所述参数的数据单元作为数据包,提取燃煤机组过去一年DCS历史运行数据;
②按负荷、气温划分二维区间,将历史数据按二维区间进行分类
a根据全年负荷主要运行区段,确定负荷优化区间:设机组在主要运行区段内的最高运行负荷为Smax,最低运行负荷为Smin,则负荷优化区间为[Smin,Smax];
b确定负荷区间数:设负荷划分间隔为LS,i代表划分的负荷区间数,i=1,2,…,m个,则负荷区间数m通过公式(1)获得;
m=(Smax-Smin)/LS (1)
c根据全年送风机入口风温变化,确定气温优化区间:设全年送风机入口最高风温为Tmax,最低风温为Tmin,则气温优化区间为[Tmin,Tmax];
d确定气温区间数:设气温划分间隔为LT,j代表划分的气温区间数,j=1,2,…,n个,则气温区间数n通过公式(2)获得:
n=(Tmax-Tmin)/LT (2)
e由步骤a~d获得第i负荷段、第j气温段的二维区间为公式(3)
{[Smin+(i-1)×LS,Smin+i×LS],[Tmin+(j-1)×LT,Tmin+j×LT]} (3)
③计算各二维区间的各类参数均值、锅炉效率均值和NOx排放均值
a将步骤①获取的机组历史运行数据的数据包按照步骤②划分的二维区间进行分类,超出负荷、气温二维区间的数据包予以舍弃;
b统计步骤a中各二维区间内数据包数量l,l=1,2,…,Li,j个;则式(3)所表述的二维区间内的运行数据表示为:xi,j,k,l,所述二维区间的各类参数均值为公式(4)所示:
Figure FDA0003849087190000011
式中:
Figure FDA0003849087190000021
表示第i负荷段、第j气温段内第k个参数的均值;
c定义锅炉效率均值=实际负荷均值/总燃料量均值,锅炉效率均值为公式(5)所示:
Figure FDA0003849087190000022
式中:
Figure FDA0003849087190000023
Figure FDA0003849087190000024
为该二维区间的总燃料量均值,
Figure FDA0003849087190000025
Figure FDA0003849087190000026
为该二维区间的实际负荷均值,其中
Figure FDA0003849087190000027
d计算NOx排放均值
Figure FDA0003849087190000028
Figure FDA0003849087190000029
e二维区间的各类参数出现无效数据的处理;
④根据优化目标筛选二维区间内的数据包
在第i负荷段、第j气温段组成的二维区间内,通过优于锅炉效率均值
Figure FDA00038490871900000210
和NOx排放均值
Figure FDA00038490871900000211
来优选数据包,筛选条件为公式(6)所示:
Figure FDA00038490871900000212
式中:xi,j,E,l和xi,j,NOx,l表示第i负荷段、第j气温段、第l个数据包的锅炉效率和NOx排放,其中
Figure FDA00038490871900000213
⑤各二维区间优于锅炉效率均值和NOx排放均值的数据包的处理
a符合优化目标的数据包经筛选后保存在原二维区间内,并统计符合优化目标的数据包数量l′,其中l′=1,2,…,Li,j′,不符合优化目标的数据包予以淘汰;
b优选后的二维区间的数据包个数少于限定值的处理;
⑥从各二维区间优于优化目标的数据集合中计算该区间各类参数均值
对于符合优化目标的数据包再次计算该二维区间各类参数均值为公式(7)所示:
Figure FDA00038490871900000214
式中:
Figure FDA00038490871900000215
表示筛选后的第i负荷段、第j气温段、第k个参数的均值;x′i,j,k,l′表示筛选后第i负荷段、第j气温段、第k个参数的第l′个数据,由此得到按负荷、气温分布的优于优化目标的优化控制规律;
⑦二维区间内异常数据的处理
对于二维区间内的异常数据,依据相邻工况运行数据的渐变性,采用正则化矩阵对突变峰值进行均值滤波;
2)生成优化控制指令
①由AGC负荷指令生成目标负荷;
②按目标负荷、送风机入口风温匹配对应负荷、气温下的优化控制规律;
③若目标负荷、送风机入口风温与对应负荷、气温下的优化控制规律存在偏差,进行负荷、气温的二维差值校正,设校正差值为Δi,j,k,则经差值校正后的优化控制信号为
Figure FDA0003849087190000031
Figure FDA0003849087190000032
④计算优化控制信号
Figure FDA0003849087190000033
与对应工况下DCS实际控制信号Sin,k的绝对差值
Figure FDA0003849087190000034
为公式(8)所示:
Figure FDA0003849087190000035
从而生成优化控制指令;
⑤依据优化控制指令的调节幅度,生成修正后的优化指令
若差值
Figure FDA0003849087190000036
在DCS实际控制信号Sink的设定幅度值以内,即
Figure FDA0003849087190000037
设定幅度值,则生成修正后的优化指令
Figure FDA0003849087190000038
若差值
Figure FDA0003849087190000039
超出或等于DCS实际控制信号Sin,k的设定幅度值,即
Figure FDA00038490871900000310
设定幅度值,则生成修正后的优化指令Sin,k·(1±设定幅度值)。
2.如权利要求1所述锅炉随气温变化的优化控制方法,其特征在于:对于步骤1)的①所述的设定数据的时间间隔值为1~3分钟。
3.如权利要求1所述锅炉随气温变化的优化控制方法,其特征在于:对于步骤1)的③二维区间的各类参数出现无效数据的处理:在所划分的不同负荷、不同气温的工况下,并不能保证二维区间的各类参数全部出现有效运行数据,此时,将采用同一负荷段内其他二维区间的同类参数的平均值予以替代。
4.如权利要求1所述锅炉随气温变化的优化控制方法,其特征在于:对于步骤1)的⑤优选后的二维区间的数据包个数少于限定值的处理:在数据包个数少于限定值的情况,认定该二维区间所含有效数据量不全,淘汰该区间已优选的数据包,并采用同一负荷段内其他二维区间的同类参数的平均值予以替代。
5.如权利要求1或4所述锅炉随气温变化的优化控制方法,其特征在于:所述步骤1)的⑤的限定值为500。
6.如权利要求1所述锅炉随气温变化的优化控制方法,其特征在于:对于所述步骤1)的⑦二维区间内异常数据的处理,依据相邻工况运行数据的渐变性,采用正则化矩阵对突变峰值进行均值替代,具体如下:
a)利用步骤②中按负荷、气温分布的二维区间建立二维数组A,数组元素表示为Ai,j,其中负荷段i=1,2,…,m,气温段j=1,2,…,n,按数组元素位于数组边界、四个顶角和内部共三种位置分类处理;
b)若区间位于数组边界,利用相邻三个区间的数据求均值为公式(9)所示:
A1,2=(A1,1+A2,2+A1,3)/3 (9)
c)若区间位于数组顶角,利用相邻两个区间的数据求均值为公式(10)所示:
Am,n=(Am-1,n+Am,n-1)/2 (10)
d)若区间位于数组内部,利用相邻四个区间的数据求均值为公式(11)所示:
Ai,j=(Ai-1,j+Ai,j-1+Ai+1,j+Ai,j+1)/4 (11)。
7.如权利要求1所述锅炉随气温变化的优化控制方法,其特征在于:对于步骤2)的⑤所述的设定幅度值为5%。
CN202010603065.0A 2020-06-29 2020-06-29 一种锅炉随气温变化的优化控制方法 Active CN111881554B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010603065.0A CN111881554B (zh) 2020-06-29 2020-06-29 一种锅炉随气温变化的优化控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010603065.0A CN111881554B (zh) 2020-06-29 2020-06-29 一种锅炉随气温变化的优化控制方法

Publications (2)

Publication Number Publication Date
CN111881554A CN111881554A (zh) 2020-11-03
CN111881554B true CN111881554B (zh) 2022-11-25

Family

ID=73157195

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010603065.0A Active CN111881554B (zh) 2020-06-29 2020-06-29 一种锅炉随气温变化的优化控制方法

Country Status (1)

Country Link
CN (1) CN111881554B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112835950B (zh) * 2020-12-09 2023-03-28 华能陕西发电有限公司 一种dcs数据挖掘的湿法脱硫系统达标排放运行曲线的获取系统及方法
CN113834091B (zh) * 2021-10-12 2022-11-11 中国矿业大学 一种燃气锅炉燃烧优化供风系统控制方法
CN113834092B (zh) * 2021-10-12 2023-02-14 中国矿业大学 一种燃气锅炉燃烧优化控制方法
CN116979106B (zh) * 2023-09-22 2023-12-15 上海氢晨新能源科技有限公司 基于工作参数组的参数控制方法、系统和用电设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103244964A (zh) * 2013-04-27 2013-08-14 国家电网公司 基于声波测量炉膛温度场的燃烧优化控制系统及控制方法
JP2014126305A (ja) * 2012-12-27 2014-07-07 Hitachi Ltd ボイラ制御装置
CN106500128A (zh) * 2017-01-12 2017-03-15 东北电力大学 一种变负荷工况下火电厂锅炉NOx排放控制方法
CN108594663A (zh) * 2018-05-22 2018-09-28 马鞍山当涂发电有限公司 一种深度调峰下660mw超临界机组agc控制方法
CN108954375A (zh) * 2018-07-18 2018-12-07 厦门邑通软件科技有限公司 锅炉节煤控制方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7756591B2 (en) * 2006-04-25 2010-07-13 Pegasus Technologies, Inc. System for optimizing oxygen in a boiler

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014126305A (ja) * 2012-12-27 2014-07-07 Hitachi Ltd ボイラ制御装置
CN103244964A (zh) * 2013-04-27 2013-08-14 国家电网公司 基于声波测量炉膛温度场的燃烧优化控制系统及控制方法
CN106500128A (zh) * 2017-01-12 2017-03-15 东北电力大学 一种变负荷工况下火电厂锅炉NOx排放控制方法
CN108594663A (zh) * 2018-05-22 2018-09-28 马鞍山当涂发电有限公司 一种深度调峰下660mw超临界机组agc控制方法
CN108954375A (zh) * 2018-07-18 2018-12-07 厦门邑通软件科技有限公司 锅炉节煤控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Mathematical modeling and experimental validation of ash deposition in a pulverized-coal boiler;Shu Zheng等;《Applied Thermal Engineering》;20170105;720-729 *
根据气温对锅炉运行的影响构造闭环燃烧优化控制的试验研究;刘瑞东 等;《节能技术》;20220331;第40卷(第232期);99-104 *

Also Published As

Publication number Publication date
CN111881554A (zh) 2020-11-03

Similar Documents

Publication Publication Date Title
CN111881554B (zh) 一种锅炉随气温变化的优化控制方法
CN103576655B (zh) 一种电站锅炉燃烧子空间建模及多目标优化方法和系统
CN110486749B (zh) 一种火电机组锅炉燃烧优化控制方法及系统
CN111853848B (zh) 一种燃煤锅炉不同层燃烧器之间燃料量分配的优化方法
CN104534507A (zh) 一种锅炉燃烧优化控制方法
CN108490794B (zh) 一种深度调峰下660mw超临界机组agc控制系统
CN107016176A (zh) 一种混合智能锅炉综合燃烧优化方法
CN110260356B (zh) 一种流化床锅炉的节能控制方法
CN111623369B (zh) 一种利用烟气含氧量信号调节锅炉入炉燃料量的控制方法
CN102750424B (zh) 一种生物质炉燃烧优化的方法
CN110263452B (zh) 一种烟道内烟气时间分布特性分析方法、系统及脱硝系统
CN112628712A (zh) 一种基于风门阻力系数的二次风闭环优化控制系统
CN106801890B (zh) 减小锅炉在升降负荷过程中主再热汽温波动的方法
CN105605609A (zh) 一种火电厂锅炉燃烧氧量优化方法
CN110986085B (zh) 基于配风方式优化的四角切圆锅炉燃烧优化控制方法
CN113390099B (zh) 一种燃煤热电联产机组自适应btu下的协调控制方法
CN1249211A (zh) 中储式球磨机制粉系统自动控制方法
CN115145152A (zh) 一种锅炉燃烧与脱硝过程协同优化控制的方法
CN107831656A (zh) 一种火电机组协调控制系统节能优化技术
CN110594780B (zh) 一种燃煤电厂锅炉在线实时燃烧优化技术方法
CN107679351B (zh) 一种汽轮机agc系统的优化方法
Zhao et al. Modeling and optimization of efficiency and NOx emission at a coal-fired utility boiler
CN112611234A (zh) 一种掺烧高炉煤气的煤粉炉智能燃烧优化控制方法
Wang et al. Expert PID control of grinding condition for cement vertical roller mill
CN116822372A (zh) 一种基于igwo锅炉运行参数的优化方法

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
TA01 Transfer of patent application right

Effective date of registration: 20210804

Address after: 132012, Changchun Road, Jilin, Jilin, 169

Applicant after: NORTHEAST DIANLI University

Applicant after: China University of Mining and Technology

Applicant after: Jiangsu Hanguang Intelligent Technology Co.,Ltd.

Address before: 132012, Changchun Road, Jilin, Jilin, 169

Applicant before: NORTHEAST DIANLI University

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant