CN108763651B - 一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法 - Google Patents
一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法 Download PDFInfo
- Publication number
- CN108763651B CN108763651B CN201810402101.XA CN201810402101A CN108763651B CN 108763651 B CN108763651 B CN 108763651B CN 201810402101 A CN201810402101 A CN 201810402101A CN 108763651 B CN108763651 B CN 108763651B
- Authority
- CN
- China
- Prior art keywords
- air
- secondary air
- baffle
- coal
- hearth
- 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
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000003245 coal Substances 0.000 claims abstract description 55
- 238000009423 ventilation Methods 0.000 claims abstract description 12
- 238000002485 combustion reaction Methods 0.000 claims abstract description 5
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 19
- 229910052799 carbon Inorganic materials 0.000 claims description 19
- 238000013519 translation Methods 0.000 claims description 19
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 11
- 239000003546 flue gas Substances 0.000 claims description 11
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 10
- 239000001301 oxygen Substances 0.000 claims description 10
- 229910052760 oxygen Inorganic materials 0.000 claims description 10
- 230000002093 peripheral effect Effects 0.000 claims description 8
- 238000000227 grinding Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000002956 ash Substances 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 239000010881 fly ash Substances 0.000 claims description 4
- 238000010298 pulverizing process Methods 0.000 claims description 4
- 239000002893 slag Substances 0.000 claims description 4
- 238000013461 design Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 abstract description 11
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 abstract description 9
- 238000012544 monitoring process Methods 0.000 abstract description 3
- 239000000463 material Substances 0.000 abstract 1
- 239000000446 fuel Substances 0.000 description 5
- 230000001133 acceleration Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- -1 radical hydrogen Chemical class 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 239000011593 sulfur Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Regulation And Control Of Combustion (AREA)
Abstract
本发明公开了一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法,无需进行专门的炉膛冷态通风试验,节约了锅炉点火启动时间,节省了人力和物力;本发明得到的配风挡板的阻力特性曲线可用于测量燃烧器喷口的空气流量,从而对整个炉膛的配风进行监控和调整,能够显著提高炉膛内煤粉的燃烧效率,降低氮氧化物的生成。
Description
技术领域
本发明属于锅炉燃烧器运行监测和控制领域,尤其涉及一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法。
背景技术
锅炉配风挡板的过流特性是指其开度与流过的空气流量的对应关系,通常采用流量系数或阻力系数随开度的变化曲线来描述,这种流量或阻力特性曲线是表征挡板调节性能的重要数据,不仅决定了送入炉膛的空气气流的分配和调节特性,还可以用来测量燃烧器喷口的空气流量,从而对整个炉膛的配风进行监控和调整,提高煤粉燃烧效率,降低氮氧化物的生成。
目前,锅炉配风挡板流量或阻力特性曲线通过炉膛冷态通风试验获得,这项试验是在维持配风挡板和炉膛总压降不变的情况下,测量不同挡板开度下的空气流量或流速。申请号为201510493413.2的中国专利公开了一种通过对多个运行状态下的风箱-炉膛压降和燃烧器喷口风速的测试,确定燃烧器二次风挡板阻力特性曲线的方法。由于大容量锅炉二次风配风挡板和燃烧器喷口的数量大,炉膛冷态通风试验工况多,风量或风速的测试工作费时费力,劳动强度较大;为便于测量燃烧器喷口风速,申请号为201520265921.0的中国专利公开了一种用于测试二次风挡板特性的装置,使喷口风速测量条件有一定改善,但炉膛冷态通风试验的测量项目、试验工况以及费时费力的特点并没有得到实质性改变。
发明内容
为解决现有技术的缺陷,本发明具体公开了一种从运行数据中提取锅炉配风挡板过流特性的方法,该方法无需进行炉膛冷态通风试验,而是从锅炉日常运行数据中通过学习算法获取锅炉配风挡板过流特性,特别适用于配备直吹式制粉系统的切向燃煤锅炉。
为实现上述目的,本发明的具体方案如下:
一种从运行数据中提取锅炉配风挡板过流特性的方法,包括以下步骤:
(1)采集n组锅炉运行工况的数据,这些运行工况应覆盖重要运行参数的典型调节范围;
(2)对锅炉二次风配风挡板进行分类,对于同类配风挡板,它们的阻力特性曲线相同;对于不同种类的配风挡板,阻力特性曲线存在差异;所有配风挡板的阻力特性曲线通过标准挡板阻力特性曲线的横坐标向左平移得到,每类挡板阻力曲线的横坐标平移量是唯一未知的参数,也是本发明要确定的参数,这些横坐标平移量参数通过如下试算方法求取;
(3)假定各类二次风配风挡板的横坐标向左平移量,根据标准挡板的阻力特性曲线,得到每类二次风配风挡板的阻力特性曲线;
(4)利用n组工况数据中各二次风配风挡板的开度,根据各挡板的阻力特性曲线,得到各工况下每个二次风挡板在所处开度下的阻力系数;
(5)利用n组工况数据中的二次风温度和压力,计算各组数据中的二次风密度;利用n组数据中的磨煤机出口一次风温度和压力,计算各工况下的一次风密度;
(6)利用n组工况数据中的运行磨煤机通风量和磨煤量,得到各工况下进入炉膛的总一次风流量和总燃煤煤流量;
(7)利用n组工况数据中二次风箱到炉膛出口压降数据,根据步骤(4)得到的二次风挡板的阻力系数,计算各工况下每个二次风喷口风速,再计算二次风喷口流量,进而得到总二次风流量;
(8)根据n组工况数据中的总一次风流量和总二次风流量,计算各工况下进入炉膛的总空气量;
(9)利用n组工况数据中入炉煤的元素分析数据、锅炉飞灰和炉渣含碳量,得到各工况下每千克入炉煤根据实际烧掉的碳计算的理论空气量;
(10)利用n组工况数据中进入炉膛的总空气量、总燃煤煤流量以及每千克入炉煤中按照实际烧掉的碳所确定的理论空气量,计算各工况下炉膛出口过剩空气系数,得到的结果称为过剩空气系数定义值;
(11)利用n组工况数据中炉膛出口烟气含氧量,计算各工况下炉膛出口过剩空气系数,得到的结果称为过剩空气系数测量值;
(12)对于每组工况数据,计算炉膛出口过剩空气系数的定义值和测量值之间偏差的平方值;计算n组工况数据的总偏差平方和;
(13)改变步骤(3)中假设的横坐标平移量,重复步骤(4)到步骤(12),寻找使过剩空气系数定义值和测量值偏差平方和最小值的横坐标平移量;
(14)使过剩空气系数定义值和测量值偏差平方值和最小的横坐标平移量,作为最终求得的横坐标平移量,根据标准挡板的阻力特性曲线,得到各个二次风配风挡板的阻力特性曲线。
进一步地,步骤(1)中,每组锅炉运行数据包括如下参数:机组负荷,所有二次风配风挡板的开度,二次风温度和压力,一次风温度和压力,运行磨煤机通风量和磨煤量,二次风箱到炉膛出口压降,入炉煤的元素分析,飞灰和炉渣含碳量,炉膛出口烟气含氧量;重要运行参数是指机组负荷、二次风配风挡板开度、二次风箱到炉膛出口压降和炉膛出口烟气含氧量,在所选择的n组运行工况数据中,这些参数的变化范围应覆盖日常运行中的调节范围。
进一步地,步骤(2)中,二次风配风挡板所分类的类数应少于步骤(1)中运行工况数据组数n;二次风配风挡板分类根据它们所对应的二次风喷嘴的功能来进行,将辅助风、周界风、紧凑燃尽风、分离燃尽风以及最底层火下风等二次风喷嘴所对应的挡板分作不同类挡板。
进一步地,步骤(2)中,标准挡板的阻力特性为:在开度为0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0时,阻力系数分别为200,65,26,16,8.3,4,1.8,0.85,0.48,0.3,该对应关系可拟合成函数曲线。
进一步地,步骤(5)中,一、二次风密度由对标准状态下空气的密度进行温度和压力修正得到。
进一步地,步骤(6)中,对于中速磨直吹式制粉系统,进入炉膛的总一次风流量等于运行磨煤机通风量之和;总燃煤煤流量等于运行磨煤机的磨煤量之和。
进一步地,步骤(7)中,二次风喷口风速根据下列关系计算:二次风箱到燃烧器喷口的阻力系数等于风道进、出口局部阻力系数和二次风挡板阻力系数之和,该阻力系数和燃烧器喷口动压的乘积为二次风箱到燃烧器喷口的压降,该压降与炉膛压降之和等于二次风箱到炉膛出口压降;
二次风喷口流量等于喷口风速、二次风密度和喷口面积的乘积;
总二次风流量等于所有二次风喷口流量之和。
步骤(8)中,进入炉膛的总空气量等于总一次风流量、总二次风流量和炉膛漏风量之和,炉膛漏风量根据炉膛设计漏风系数得到。
进一步地,步骤(9)中,实际烧掉的碳是指1kg煤含有的碳元素质量减去1kg煤所产生的灰渣中的碳的质量;对于n组工况数据,当入炉煤质稳定时,这些工况可采用相同的入炉煤元素分析数据;理论空气量是指完全燃烧所需要的空气量。
进一步地,步骤(10)中,进入炉膛的总空气流量除以总燃煤流量与每千克入炉煤按照实际烧掉的碳计算的理论空气量的乘积,所得的商即为炉膛出口过剩空气系数的定义值。
进一步地,步骤(12)中,n组工况数据的总偏差平方和是指每组工况数据中炉膛出口过剩空气系数的定义值和测量值之间偏差的平方值之和;
步骤(13)中,改变各挡板阻力特性曲线的横坐标平移量,寻找过剩空气系数偏差平方和最小值,可采用试算算法;为求解提高效率,也可采用目标函数求最小值的数学算法来进行,其中目标函数为过剩空气系数定义值和测量值偏差平方值和。
综上,本发明的上述技术方案的有益效果如下:
为得到锅炉配风挡板阻力特性曲线,采用传统的炉膛冷态通风试验法,测试工作费时费力,劳动强度大,本发明公开了一种从运行数据中提取锅炉配风挡板过流特性的方法,无需进行专门的炉膛冷态通风试验,节约了锅炉点火启动时间,节省了人力和物力;本发明得到的配风挡板的阻力特性曲线可用于测量燃烧器喷口的空气流量,从而对整个炉膛的配风进行监控和调整,能够显著提高炉膛内煤粉的燃烧效率,降低氮氧化物的生成。
附图说明
图1为标准挡板不同开度下的阻力系数值、标准挡板阻力系数的拟合曲线以及经过横坐标平移得到的实际挡板的阻力系数曲线;
图2为实例锅炉的燃烧器喷口布置、喷口面积及配置的二次风配风挡板;
图3为过剩空气系数偏差平方和函数f(x1,x2,......,x6)的计算流程;
图4为AA层二次风挡板阻力系数ζ1与开度γ1的坐标图;
图5为辅助二次风挡板阻力系数ζ2与开度γ2的坐标图;
图6为油枪二次风挡板阻力系数ζ3与开度γ3的坐标图;
图7为SOFA二次风挡板阻力系数ζ6与开度γ6的坐标图。
具体实施方式
以下结合附图对本发明的特征和原理进行详细说明,所举实施例仅用于解释本发明,并非以此限定本发明的保护范围。
本发明公开了一种从运行数据中提取锅炉配风挡板过流特性的方法,具体包括以下步骤:
(1)从电厂分散控制系统的数据库中选取n个工况的锅炉运行数据,这些工况的选取原则是,主要参数如机组负荷、二次风配风挡板开度、二次风箱到炉膛出口压降和炉膛出口烟气含氧量的变化范围应在日常运行调节区间内,工况的个数n应大于二次风喷口的种类数m。
(2)对于切向燃烧锅炉,二次风喷口按功能分为辅助风、周界风、紧凑燃尽风、分离燃尽风以及最底层火下风等,不同种类喷口对应的二次风配风挡板具有不同的阻力特性曲线,某类二次风配风挡板的阻力特性曲线可通过标准挡板阻力特性曲线的横坐标向左平移得到,图1是标准挡板不同开度下阻力系数值、标准挡板阻力系数的拟合曲线以及标准挡板阻力系数曲线向左平移x个单位横坐标后所得的曲线,平移x个单位横坐标后所得的曲线可作为实际二次风配风挡板的阻力特性曲线,其中,x为待定参数,与配风挡板的种类有关。
图1中标准挡板阻力系数拟合函数曲线如下:
ζ=ζ0+ae-γ/p+be-γ/q
式中γ为挡板开度,ζ为阻力系数,ζ0=-0.457,a=97.258,b=689.233,p=0.202,q=0.063。
图1中平移x个单位横坐标后得到的实际二次风配风挡板阻力特性曲线如下:
ζ=ζ0+ae-(γ+x)/p+be-(γ+x)/q
式中x个为向左平移的横坐标量。
(3)假定m类二次风配风挡板的横坐标向左平移量分别为x1,x2,……,xm,如图1所示,对标准挡板的阻力特性曲线进行横坐标平移,得到各类二次风配风挡板的阻力特性曲线。
(4)每个数据工况中,各个二次风配风挡板的开度分别为γ1,γ2,......γN,根据各个挡板的阻力特性曲线,计算每个二次风挡板在所处开度下的阻力系数值ζ1,ζ2,......ζN。
(5)根据每个数据工况中的一次风温度和压力、二次风温度和压力,分别按下式计算一次风密度、二次风密度。
ρ=0.003483(p0+p)/(273.15+t)
式中ρ为一、二次风密度,kg/m3;p为一、二次风压力,Pa;t为一、二次风温度,℃;p0为标准大气压力,Pa。
(6)根据每个数据工况中的运行磨煤机通风量Mfi和磨煤量Mci,计算各工况下进入炉膛的总一次风流量M1和总燃煤煤流量Mc。对于配有中速磨直吹式制粉系统的锅炉:
式中Mfi、Mci分别为第i台磨煤机的通风量和给煤量,kg/s。
(7)根据每个数据工况中二次风箱到炉膛出口压降DP,计算每个二次风喷口风速、风量以及总二次风流量。
第j个二次风喷口的风速uj按下式计算
式中ΔPb为风箱到喷口出口的压降,Pa;ζj为配风挡板阻力系数;ζ1为二次风箱到风道进口的局部阻力系数,ζ2为二次风道到炉膛入口的阻力系数,ρ2为二次风密度,kg/m3。
ΔPb等于二次风箱到炉膛出口压降ΔP与炉膛压降ΔPL之差,按下式计算:
ΔPb=ΔP-ΔPL
炉膛压降ΔPL按下式计算:
式中M2为二次风总流量,kg/s;其余符号同前。
单个二次风喷口的流量根据喷口风速uj、二次风密度ρ2和喷口面积A2j计算,二次风总流量M2等于所有二次风喷口流量之和:
式中β为喷口附加漏风修正系数。
二次风喷口风速uj的计算方程右边隐含了uj值,因此,可采用迭代法求解,迭代的变量选为风箱到喷口出口的压降ΔPb,迭代收敛后求得uj值,并最终求得二次风总流量M2。
(8)根据每个数据工况中得到的总一次风流量和总二次风流量,计算进入炉膛的总空气量,炉膛漏风系数为ΔαL,得到进入炉膛的总空气量Mt为:
Mt=(1+ΔαL)(M1+M2)
(9)每个数据工况中,每千克入炉煤根据实际烧掉的碳计算的理论空气量m0按下式计算:
m0=0.1151Cb+0.342Har+0.0431Sar-0.0432Oar
其中实际烧掉的碳Cb为
(10)每个数据工况中,过剩空气系数定义值αLd为
(11)每个数据工况中,过剩空气系数测量值αLm按下式计算:
αLm=21/(21-γO2)
式中γO2为炉膛出口烟气含氧量,%。
(12)每个数据工况中,计算炉膛出口过剩空气系数的定义值和测量值之间偏差的平方值,并得到n组数据工况的总偏差平方和为:
式中αLd(i)为第i个数据工况中过剩空气系数定义值,αLm(i)为第i个数据工况中过剩空气系数测量值。
(13)将步骤(3)中假设的横坐标平移量x1,x2,……,xm视为总偏差平方和的函数,即
可以用试算的方法确定使函数f达到最小时的x1,x2,......,xm值;为提高求解效率,也可用采用求函数f(x1,x2,......,xm)最小值的方法,搜索使函数f达到最小值时的x1,x2,......,xm的值。
(14)将求得的x1,x2,......,xm值代替图1中的x,得到每个二次风挡板的阻力特性曲线,该挡板的阻力特性曲线也可用下面的函数表示:
式中ζj为第j个二次风配风挡板的阻力系数,γj为第j个二次风配风挡板的开度,xj为第j个二次风配风挡板的横坐标平移量。
实施例一:
某300MW切向燃烧煤粉锅炉,采用直吹式制粉系统,配有5台中速磨煤机,同一台磨煤机出口的一次风分配到同一层的4只燃料喷口,总计有5层共20只燃料喷口,从下至上,燃料喷口依次记为A、B、C、D、E层;空气喷口共19层,每层4只,共76只;这19层空气喷口分为6类:①AA层二次风,位于最底层;②辅助二次风:共3层,分别是CD、DE和EE层二次风;③油枪二次风:共2层,分别是AB、BC层油枪二次风;④燃料喷口周界二次风:共5层,分别是A周界风、B周界风、C周界风、D周界风、E周界风,⑤贴壁二次风,共4层,分别是BC、DE、EE和OFA层贴壁二次风;⑥SOFA二次风,共4层,分别为SOFA1、SOFA2、SOFA3和SOFA4。
上述19层空气喷口由15层二次风配风挡板γAA、γA、γAB、γB、γBC、γC、γCD、γD、γDE、γE、γEE、γSOFA1、γSOFA2、γSOFA3、γSOFA4控制,见图2,这些挡板采用分层控方式,即同层4只挡板同步控制,因此,可作为1个挡板处理;同类空气喷口对应的风道尺寸和挡板结构相同,安装时预留相同的零位空隙,因此具有相同的过流特性,但不同类的配风挡板因结构和安装预留间隙不同,具有不同的过流特性。二次风喷口配风挡板共分6类。燃烧器喷口和二次风挡板的布置见图2。该锅炉基本数据如下:
同层四只喷口的面积之和见图2;
二次风挡板数量N=15;
二次风配风挡板分类数量m=6;
燃烧器中心线到炉膛出口的高度H=34.3m;
二次风箱到风道进口的局部阻力系数ζ1=0.5;
二次风道到炉膛入口的阻力系数ζ2=1.0;
炉膛阻力系数ζL=0.78;
炉膛漏风系数ΔαL=0.05;
炉内烟气密度ρL=0.24kg/m3;
喷口附加漏风修正系数β=1.277;
大气压力p0=100900Pa;
重力加速度g=9.8m/s2。
从电厂控制系统DCS中选择了12组数据,即n=12,其中机组负荷Me、二次风箱到炉膛出口压降△P、炉膛出口烟气含氧量γO2和二次风配风挡板的开度γ的变化范围较大,覆盖了日常运行调节区间,选取的数据工况个数n(=12)大于挡板分类数量m(=6),这n组数据见下表1:
表1
假定AA层二次风、辅助二次风、油枪二次风、燃料喷口周界二次风、贴壁二次风、SOFA二次风等六类配风挡板的横坐标平移量分别为x1,x2,...,x6,它们作为自变量,按照图3给出的框图计算采用迭代算法计算二次风喷口的风速uj、总二次风量M2以及12组数据的过剩空气系数偏差平方和函数f(x1,x2,......,x6),求函数f的最小值,并得到函数f取得最小值时的x1,x2,...,x6的值,结果为x1=0.173,x2=0.155,x3=0.25,x4=0.02,x5=0.169,x6=0.25,分别将这些值代入中,可得到第j类挡板阻力系数ζj与开度γj的函数关系式,图4到图7是主要二次风挡板阻力系数ζj随开度γj的变化,图中实线为本发明得到的挡板阻力系数与开度的函数曲线,图中实点是采用传统的炉膛冷态通风试验得到的挡板阻力系数与开度关系的数据,可以看出,本发明得到的结果与传统的炉膛冷态通风试验方法接近。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (10)
1.一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法,其特征在于:
(1)采集n组锅炉运行工况的数据,这些运行工况应覆盖重要运行参数的典型调节范围;重要运行参数是指机组负荷、二次风配风挡板开度、二次风箱到炉膛出口压降和炉膛出口烟气含氧量;
(2)对锅炉二次风配风挡板进行分类,对于同类配风挡板,它们的阻力特性曲线相同;对于不同种类的配风挡板,阻力特性曲线存在差异;所有配风挡板的阻力特性曲线通过标准挡板阻力特性曲线的横坐标向左平移得到,每类挡板阻力曲线的横坐标平移量是唯一未知的参数,也是本方法要确定的参数,这些横坐标平移量参数通过如下试算方法求取;
(3)假定各类二次风配风挡板的横坐标向左平移量,根据标准挡板的阻力特性曲线,得到每类二次风配风挡板的阻力特性曲线;
(4)利用n组工况数据中各二次风配风挡板的开度,根据各挡板的阻力特性曲线,得到各工况下每个二次风挡板在所处开度下的阻力系数;
(5)利用n组工况数据中的二次风温度和压力,计算各组数据中的二次风密度;利用n组数据中的磨煤机出口一次风温度和压力,计算各工况下的一次风密度;
(6)利用n组工况数据中的运行磨煤机通风量和磨煤量,得到各工况下进入炉膛的总一次风流量和总燃煤煤流量;
(7)利用n组工况数据中二次风箱到炉膛出口压降数据,根据步骤(4)得到的二次风挡板的阻力系数,计算各工况下每个二次风喷口风速,再计算二次风喷口流量,进而得到总二次风流量;
(8)根据n组工况数据中的总一次风流量和总二次风流量,计算各工况下进入炉膛的总空气量;
(9)利用n组工况数据中入炉煤的元素分析数据、锅炉飞灰和炉渣含碳量,得到各工况下每千克入炉煤根据实际烧掉的碳计算的理论空气量;
(10)利用n组工况数据中进入炉膛的总空气量、总燃煤煤流量以及每千克入炉煤中按照实际烧掉的碳所确定的理论空气量,计算各工况下炉膛出口过剩空气系数,得到的结果称为过剩空气系数定义值;
(11)利用n组工况数据中炉膛出口烟气含氧量,计算各工况下炉膛出口过剩空气系数,得到的结果称为过剩空气系数测量值;
(12)对于每组工况数据,计算炉膛出口过剩空气系数的定义值和测量值之间偏差的平方值;计算n组工况数据的总偏差平方和;
(13)改变步骤(3)中假设的横坐标平移量,重复步骤(4)到步骤(12),寻找使过剩空气系数定义值和测量值偏差平方和最小值的横坐标平移量;
(14)使过剩空气系数定义值和测量值偏差平方值和最小的横坐标平移量,作为最终求得的横坐标平移量,根据标准挡板的阻力特性曲线,得到各个二次风配风挡板的阻力特性曲线。
2.根据权利要求1所述的一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法,其特征在于,步骤(1)中,每组锅炉运行数据包括如下参数:机组负荷,所有二次风配风挡板的开度,二次风温度和压力,一次风温度和压力,运行磨煤机通风量和磨煤量,二次风箱到炉膛出口压降,入炉煤的元素分析,飞灰和炉渣含碳量,炉膛出口烟气含氧量;在所选择的n组运行工况数据中,这些参数的变化范围应覆盖日常运行中的调节范围。
3.根据权利要求1所述的一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法,其特征在于,步骤(2)中,二次风配风挡板所分类的类数应少于步骤(1)中运行工况数据组数n;二次风配风挡板分类根据它们所对应的二次风喷嘴的功能来进行,将辅助风、周界风、紧凑燃尽风、分离燃尽风以及最底层火下风这些二次风喷嘴所对应的挡板分作不同类挡板。
4.根据权利要求1所述的一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法,其特征在于,步骤(2)中,标准挡板的阻力特性为:在开度为0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0时,阻力系数分别为200,65,26,16,8.3,4,1.8,0.85,0.48,0.3,将开度与阻力系数的对应关系拟合成函数曲线。
5.根据权利要求1所述的一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法,其特征在于,步骤(5)中,一、二次风密度由对标准状态下空气的密度进行温度和压力修正得到。
6.根据权利要求1所述的一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法,其特征在于,步骤(6)中,对于中速磨直吹式制粉系统,进入炉膛的总一次风流量等于运行磨煤机通风量之和;总燃煤煤流量等于运行磨煤机的磨煤量之和。
7.根据权利要求1所述的一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法,其特征在于,步骤(7)中,二次风喷口风速根据下列关系计算:二次风箱到燃烧器喷口的阻力系数等于风道进、出口局部阻力系数和二次风挡板阻力系数之和,该阻力系数和燃烧器喷口动压的乘积为二次风箱到燃烧器喷口的压降,该压降与炉膛压降之和等于二次风箱到炉膛出口压降;
二次风喷口流量等于喷口风速、二次风密度和喷口面积的乘积;
总二次风流量等于所有二次风喷口流量之和;
步骤(8)中,进入炉膛的总空气量等于总一次风流量、总二次风流量和炉膛漏风量之和,炉膛漏风量根据炉膛设计漏风系数得到。
8.根据权利要求1所述的一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法,其特征在于,步骤(9)中,实际烧掉的碳是指1kg煤含有的碳元素质量减去1kg煤所产生的灰渣中的碳的质量;对于n组工况数据,当入炉煤质稳定时,这些工况可采用相同的入炉煤元素分析数据;理论空气量是指完全燃烧所需要的空气量。
9.根据权利要求1所述的一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法,其特征在于,步骤(10)中,进入炉膛的总空气流量除以总燃煤流量与每千克入炉煤按照实际烧掉的碳计算的理论空气量的乘积,所得的商即为炉膛出口过剩空气系数的定义值。
10.根据权利要求1所述的一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法,其特征在于,步骤(12)中,n组工况数据的总偏差平方和是指每组工况数据中炉膛出口过剩空气系数的定义值和测量值之间偏差的平方值之和;
步骤(13)中,改变各挡板阻力特性曲线的横坐标平移量,寻找过剩空气系数偏差平方和最小值,采用目标函数求最小值的数学算法来进行,其中目标函数为过剩空气系数定义值和测量值偏差平方值和。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810402101.XA CN108763651B (zh) | 2018-04-28 | 2018-04-28 | 一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810402101.XA CN108763651B (zh) | 2018-04-28 | 2018-04-28 | 一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108763651A CN108763651A (zh) | 2018-11-06 |
CN108763651B true CN108763651B (zh) | 2022-04-12 |
Family
ID=64008709
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810402101.XA Active CN108763651B (zh) | 2018-04-28 | 2018-04-28 | 一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108763651B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109681907B (zh) * | 2018-12-29 | 2024-02-27 | 浙江大学 | 一种燃煤电厂锅炉变负荷风煤匹配控制装置与方法 |
CN110081457B (zh) * | 2019-06-20 | 2023-11-28 | 吉林大学 | 一种预防锅炉水冷壁爆管的二次风风管装置及其控制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN204556242U (zh) * | 2015-04-29 | 2015-08-12 | 华电国际电力股份有限公司技术服务中心 | 锅炉二次风门挡板特性在线测量装置 |
CN104848248A (zh) * | 2014-02-19 | 2015-08-19 | 宝山钢铁股份有限公司 | 煤气混烧锅炉风门挡板开度的控制方法 |
CN105160158A (zh) * | 2015-08-12 | 2015-12-16 | 山东电力研究院 | 锅炉二次风挡板特性试验数据的处理方法 |
CN105783025A (zh) * | 2016-03-11 | 2016-07-20 | 国网山东省电力公司电力科学研究院 | 一种监测低NOx切向燃煤锅炉炉内风粉分布的方法 |
CN105841141A (zh) * | 2016-04-08 | 2016-08-10 | 国网山东省电力公司电力科学研究院 | 一种切向燃烧摆动喷嘴附加进风的运行评估方法 |
EP3239611A1 (en) * | 2014-12-25 | 2017-11-01 | Fuji Electric Co., Ltd. | Combustion control system, combustion control method, combustion control program, and computer-readable recording medium |
CN207264430U (zh) * | 2017-09-30 | 2018-04-20 | 中电华创电力技术研究有限公司 | 正压式中速磨煤机制粉装置煤粉细度实时监测系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130035794A1 (en) * | 2011-08-03 | 2013-02-07 | Behzad Imani | Method and system for controlling building energy use |
-
2018
- 2018-04-28 CN CN201810402101.XA patent/CN108763651B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104848248A (zh) * | 2014-02-19 | 2015-08-19 | 宝山钢铁股份有限公司 | 煤气混烧锅炉风门挡板开度的控制方法 |
EP3239611A1 (en) * | 2014-12-25 | 2017-11-01 | Fuji Electric Co., Ltd. | Combustion control system, combustion control method, combustion control program, and computer-readable recording medium |
CN204556242U (zh) * | 2015-04-29 | 2015-08-12 | 华电国际电力股份有限公司技术服务中心 | 锅炉二次风门挡板特性在线测量装置 |
CN105160158A (zh) * | 2015-08-12 | 2015-12-16 | 山东电力研究院 | 锅炉二次风挡板特性试验数据的处理方法 |
CN105783025A (zh) * | 2016-03-11 | 2016-07-20 | 国网山东省电力公司电力科学研究院 | 一种监测低NOx切向燃煤锅炉炉内风粉分布的方法 |
CN105841141A (zh) * | 2016-04-08 | 2016-08-10 | 国网山东省电力公司电力科学研究院 | 一种切向燃烧摆动喷嘴附加进风的运行评估方法 |
CN207264430U (zh) * | 2017-09-30 | 2018-04-20 | 中电华创电力技术研究有限公司 | 正压式中速磨煤机制粉装置煤粉细度实时监测系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108763651A (zh) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102425807B (zh) | 一种煤粉锅炉燃烧前馈反馈复合优化控制方法 | |
CN105276611B (zh) | 火电厂锅炉燃烧调整优化方法与系统 | |
CN108800191B (zh) | 一种切向燃煤锅炉二次风配风的动力优化方法 | |
CN110848733B (zh) | 一种基于煤质在线监测的燃烧优化方法 | |
CN112524637B (zh) | 一种基于风粉和co在线监测的锅炉燃烧优化方法和系统 | |
CN109084324B (zh) | 生物质锅炉的燃烧风量控制系统及控制方法 | |
WO2017133316A1 (zh) | 对冲燃烧进风量的确定方法、装置及自动控制系统 | |
CN109058969B (zh) | 一种超超临界二次再热塔式锅炉运行控制方法 | |
CN105783025B (zh) | 一种监测低NOx切向燃煤锅炉炉内风粉分布的方法 | |
CN108763651B (zh) | 一种从锅炉运行数据中提取燃烧器配风挡板过流特性的方法 | |
CN105841141B (zh) | 一种切向燃烧摆动喷嘴附加进风的运行评估方法 | |
CN105605609B (zh) | 一种火电厂锅炉燃烧氧量优化方法 | |
CN110135065B (zh) | 一种基于煤质软测量技术的燃烧配风优化方法 | |
CN105160158A (zh) | 锅炉二次风挡板特性试验数据的处理方法 | |
CN113609663B (zh) | 一种基于物理模型与历史数据分析的燃烧优化方法 | |
CN105240868A (zh) | 基于风煤比煤质校正的锅炉自动燃烧调整控制方法 | |
CN111853848B (zh) | 一种燃煤锅炉不同层燃烧器之间燃料量分配的优化方法 | |
CN2711539Y (zh) | 循环流化床锅炉燃烧自动调节系统 | |
CN113341713B (zh) | 一种燃煤机组煤风同步动态协控方法 | |
CN107143873A (zh) | 一种超超临界前后对冲锅炉氮氧化物燃烧优化调整方法 | |
CN113390099A (zh) | 一种燃煤热电联产机组自适应btu下的协调控制方法 | |
TW201837628A (zh) | 模擬結果之評價裝置及方法 | |
CN112945600B (zh) | 基于水冷壁高温腐蚀防治的煤粉锅炉一体化综合调整方法 | |
CN106895434A (zh) | 塔式炉空气分级燃烧方式下的低氮燃烧自动控制方法 | |
CN115095884A (zh) | 基于数字孪生的锅炉精细化配风方法和锅炉系统 |
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 |