CN113283091B - 一种助燃干空气量和湿烟气量的实时计算方法 - Google Patents
一种助燃干空气量和湿烟气量的实时计算方法 Download PDFInfo
- Publication number
- CN113283091B CN113283091B CN202110598225.1A CN202110598225A CN113283091B CN 113283091 B CN113283091 B CN 113283091B CN 202110598225 A CN202110598225 A CN 202110598225A CN 113283091 B CN113283091 B CN 113283091B
- Authority
- CN
- China
- Prior art keywords
- coal
- unit
- wet
- combustion
- dry air
- 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
- 239000000779 smoke Substances 0.000 title claims abstract description 80
- 238000004364 calculation method Methods 0.000 title claims description 9
- 239000003245 coal Substances 0.000 claims abstract description 211
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims abstract description 36
- 239000003546 flue gas Substances 0.000 claims abstract description 36
- 238000002485 combustion reaction Methods 0.000 claims abstract description 35
- 238000000034 method Methods 0.000 claims abstract description 34
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 33
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 33
- 239000001301 oxygen Substances 0.000 claims abstract description 33
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 16
- 229920006395 saturated elastomer Polymers 0.000 claims description 6
- 238000000611 regression analysis Methods 0.000 claims description 4
- 238000010438 heat treatment Methods 0.000 claims description 2
- 108700041286 delta Proteins 0.000 claims 1
- 238000004458 analytical method Methods 0.000 description 10
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 6
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 6
- 239000007789 gas Substances 0.000 description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 3
- 229910021529 ammonia Inorganic materials 0.000 description 3
- 229910052799 carbon Inorganic materials 0.000 description 3
- 229910052739 hydrogen Inorganic materials 0.000 description 3
- 239000001257 hydrogen Substances 0.000 description 3
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 229910052757 nitrogen Inorganic materials 0.000 description 3
- 229910052717 sulfur Inorganic materials 0.000 description 3
- 239000011593 sulfur Substances 0.000 description 3
- 238000000921 elemental analysis Methods 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 229910052729 chemical element Inorganic materials 0.000 description 1
- 239000002817 coal dust Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004043 responsiveness Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
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)
- Control Of Steam Boilers And Waste-Gas Boilers (AREA)
Abstract
本发明涉及一种助燃干空气量和湿烟气量的实时计算方法,包括如下步骤:分别确定机组负荷与设计给煤量和设计烟气运行湿基氧量的函数关系;拟合出单位质量燃煤理论湿烟气量和助燃干空气量的回归函数,并判断其准确性,计算单位质量燃煤实际助燃干空气量;分别计算每小时燃煤总量的所需助燃干空气量和产生湿烟气量的预测值与设计值,验证相对偏差是否在可接受范围内;计算实际给煤低位热值;计算进入锅炉的助燃干空气量和省煤器出口湿烟气量。本发明提出了一种以燃煤低位热值为自变量的助燃干空气量和湿烟气量的计算方法,可在不需燃煤元素组份的情况下准确计算出助燃干空气量和湿烟气量,解决了锅炉湿烟气量和助燃干空气量无法实时准确获得的问题。
Description
技术领域
本发明属于燃煤机组技术领域,具体涉及一种助燃干空气量和湿烟气量的实时计算方法。
背景技术
长期以来,对于燃煤机组的助燃干空气量和出口湿烟气量数据的准确实时获取是目前燃煤电厂都存在的问题,而对于助燃干空气量和出口湿烟气量数据的准确获取也直接影响到燃烧工况和SCR喷氨的调节响应性,进而影响到整个燃煤锅炉和SCR脱硝的性能状态。
参见中国专利CN104021290A公开了一种煤粉与高炉煤气混烧条件下的燃料燃烧计算方法,该方法依据燃料具体化学元素组份含量(如水分、灰分、元素碳、氢、氧、氮、硫等),提出了燃烧烟气量和助燃干空气量的计算方法。但由于在燃煤机组运行过程无法实时获得燃煤的化学组成,因此该计算方法无法实现现场实际应用。
发明内容
本发明的目的是提供一种助燃干空气量和湿烟气量的实时计算方法,用于解决锅炉湿烟气量和助燃干空气量无法实时计算的问题,为通过预置锅炉运行湿基氧量和SCR喷氨量提高其调节响应性。
为达到上述目的,本发明采用的技术方案是:
一种助燃干空气量和湿烟气量的实时计算方法,包括如下步骤:
S1:在机组稳态运行模式下,获取不同负荷下的设计给煤量和烟气运行湿基氧量,进而确定机组负荷L与设计给煤量Wcoal,b的对应函数关系、机组负荷L与设计烟气运行湿基氧量O2,b的对应函数关系;
S2:采用回归分析方法,拟合得到基于低位热值的单位质量燃煤理论湿烟气量Wflue,th,wet和单位质量燃煤理论助燃干空气量Wair,th,dry的回归函数,并判断回归函数的准确性,计算得到单位质量燃煤实际助燃干空气量Wair,act,dry,根据空气温度、相对湿度和大气压力,计算空气饱和蒸汽压和单位体积干空气中的水蒸气占比/>进而计算单位质量燃煤实际助燃湿空气量Wair,act,wet和单位质量燃煤实际产生的湿烟气量Wflue,act,wet;
S3:计算不同负荷下每小时燃煤总量的所需助燃干空气量预测值Wair,h和产生湿烟气量预测值Wflue,h,通过比较不同负荷下每小时燃煤总量的所需助燃干空气量预测值Wair,h与其设计值Wair,b、产生湿烟气量预测值Wflue,h与其设计值Wflue,b,验证相对偏差δ是否在可接受范围内;
S4:获取机组的实际运行负荷Li和实际给煤量Wcoal,act,通过机组负荷L与设计给煤量Wcoal,b的对应函数关系计算实际运行负荷Li下的设计给煤量Wcoal,b,i,并计算实际给煤低位热值Qnet,act;
S5:根据机组的实际给煤量Wcoal,act、实际运行湿基氧量O2,act以及实际给煤低位热值Qnet,act,计算出实时进入锅炉的助燃干空气量BWair,act,dry和助燃湿空气量BWair,act,wet,锅炉省煤器出口烟气量BWflue,act,wet。
优选地,所述的机组负荷L与设计给煤量Wcoal,b的对应函数关系、机组负荷L与设计烟气运行湿基氧量O2,b的对应函数关系分别为:
Wcoal,b=f(L,Qnet,b),
O2,b=g(L),
式中:
L为机组负荷,单位为MW,
Wcoal,b为设计给煤量,单位为t/h,
Qnet,b为设计煤低位热值,单位为MJ/kg,
O2,b为设计烟气运行湿基氧量,单位为%。
优选地,所述的单位质量燃煤理论湿烟气量Wflue,th,wet和所述的单位质量燃煤理论助燃干空气量Wair,th,dry通过拟合多组已有的电站锅炉煤种样品的元素分析和工业分析数据获得。
进一步优选地,所述的工业分析数据包括收到基灰、收到基水、低位热值,所述的元素分析数据包括碳、氢、氧、氮、硫成分。
优选地,所述的单位质量燃煤理论湿烟气量Wflue,th,wet和所述的单位质量燃煤理论助燃干空气量Wair,th,dry分别为
Wflue,th,wet=α1×Qnet+β1,
Wair,th,dry=α2×Qnet+β2,
式中:
Wflue,th,wet为单位质量燃煤理论湿烟气量,单位为m3/kg,
Wair,th,dry为单位质量燃煤理论助燃干空气量,单位为m3/kg,
Qnet为煤低位热值,单位为MJ/kg,
α1、α2、β1、β2均为常数。
优选地,所述的回归函数的准确性通过方差进行判断,所述的方差应大于0.9。
优选地,所述的空气饱和蒸气压单位体积干空气中的水蒸汽体积占比为/>单位质量燃煤实际助燃干空气量Wair,act,dry、单位质量燃煤实际助燃湿空气量Wair,act,wet、单位质量燃煤实际产生的湿烟气量Wflue,act,wet分别为:
式中:
为空气饱和蒸汽压,单位为Pa,
Tair为空气温度,单位为℃,
ф为空气相对湿度,单位为%,
Patm为大气压力,单位为Pa,
为单位体积干空气中的水蒸汽体积占比,
Wair,act,dry为单位质量燃煤实际助燃干空气量,单位为m3/kg,
Wair,th,dry为单位质量燃煤理论助燃干空气量,单位为m3/kg,
Wflue,th,wet为单位质量燃煤理论湿烟气量,单位为m3/kg,
Wflue,act,wet为单位质量燃煤湿基产生的湿烟气量,单位为m3/kg,
O2为烟气运行湿基氧量,单位为%。
优选地,所述的每小时入炉燃煤总量的所需助燃干空气量预测值Wair,h和产生的湿烟气量预测值Wflue,h分别为:
Wair,h=Wair,act,dry×Wcoal,b×(1-γ)×1000,
Wflue,h=Wflue,act,wet×Wcoal,b×(1-γ)×1000,
式中:
Wair,h为每小时燃煤总量的所需助燃干空气量预测值,单位为m3/h,
Wflue,h为每小时燃煤总量的产生湿烟气量预测值,单位为m3/h,
Wair,act,dry为单位质量燃煤实际助燃干空气量,单位为m3/kg,
Wflue,act,wet为单位质量燃煤实际产生的湿烟气量,单位为m3/kg,
Wcoal,b为设计给煤量,单位为t/h,
γ为入炉煤的未燃烧损失比例。
优选地,所述的相对偏差δ的可接受范围为-5%~+5%,所述的相对偏差δ包括每小时燃煤总量的所需助燃干空气量预测值Wair,h与其设计值Wair,b的相对偏差δ1、产生湿烟气量预测值Wflue,h与其设计值Wflue,b的相对偏差δ2,所述的相对偏差δ1和所述的相对偏差δ2分别为:
式中:
Wair,h为每小时燃煤总量的所需助燃干空气量预测值,单位为m3/h,
Wflue,h为每小时燃煤总量的产生湿烟气量预测值,单位为m3/h,
Wair,b为每小时燃煤总量的所需助燃干空气量设计值,单位为m3/h,
Wflue,b为每小时燃煤总量的产生湿烟气量设计值,单位为m3/h。
优选地,所述的实际给煤低位热值Qnet,act为:
式中:
Qnet,act为实际给煤低位热值,单位为MJ/kg,
Qnet,b为设计给煤低位热值,单位为MJ/kg,
Wcoal,act为实际运行负荷Li下的实际给煤量,单位为t/h,
Wcoal,b,i为实际运行负荷Li下的设计给煤量,单位为t/h。
优选地,所述的助燃干空气量BWair,act,dry和助燃湿空气量BWair,act,wet分别为:
式中:
BWair,act,dry为助燃干空气量,单位为m3/h,
BWair,act,wet为助燃湿空气量,单位为m3/h,
Wflue,th,wet为单位质量燃煤理论湿烟气量,单位为m3/kg,
Wair,th,dry为单位质量燃煤理论助燃干空气量,单位为m3/kg,
Wcoal,act为实际运行负荷Li下的实际给煤量,单位为t/h,
O2,act为实际烟气运行湿基氧量,单位为%,
为单位体积干空气中的水蒸汽体积占比。
优选地,所述的出口湿烟气量BWflue,act,wet为:
式中:
BWflue,act,wet为出口湿烟气量,单位为m3/h,
Wflue,th,wet为单位质量燃煤理论湿烟气量,单位为m3/kg,
Wair,th,dry为单位质量燃煤理论助燃干空气量,单位为m3/kg,
Wcoal,act为实际运行负荷Li下的实际给煤量,单位为t/h,
O2,act为实际烟气运行湿基氧量,单位为%,
为单位体积干空气中的水蒸汽体积占比。
优选地,所述的实时的实际运行负荷Li、实际给煤量Wcoal,act、实际烟气运行湿基氧量O2,act、实际给煤低位热值Qnet,act均通过机组的DCS系统获取。
由于上述技术方案运用,本发明与现有技术相比具有下列优点:
本发明基于回归分析方法,并结合机组负荷与燃煤量函数关系计算出的实际给煤低位热值,提出了一种以燃煤低位热值为自变量的助燃干空气量和湿烟气量的计算方法,可在不需给煤元素组份的情况下准确计算出机组当前负荷下的助燃干空气量和湿烟气量,解决了锅炉烟气量和助燃干空气量无法实时准确获得的问题,可通过预置锅炉燃烧运行氧量和SCR喷氨量提高调节响应性。
具体实施方式
下面对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种助燃干空气量和湿烟气量的实时计算方法,包括如下步骤:
S1:在机组稳态运行模式下,获得不同负荷下的设计给煤量和烟气运行湿基氧量,进而确定以机组负荷L为自变量的设计给煤量Wcoal,b折线函数、以机组负荷L为自变量的设计烟气运行湿基氧量O2,b折线函数,具体而言:
Wcoal,b=f(L,Qnet,b),
O2,b=g(L),
式中:
L为机组负荷,单位为MW,
Wcoal,b为设计给煤量,单位为t/h,
Qnet,b为设计煤低位热值,单位为MJ/kg,
O2,b为设计烟气运行湿基氧量,单位为%。
S2:采用回归分析方法,通过拟合已有的500多组电站锅炉煤种样品的元素分析和工业分析数据,得到基于低位热值的单位质量燃煤理论湿烟气量Wflue,th,wet和单位质量燃煤理论助燃干空气量Wair,th,dry,并通过方差判断回归函数的准确性,方差大于0.9则表明拟合的准确性高,进而计算单位质量燃煤实际助燃干空气量Wair,act,dry,并根据实时测得的空气温度、相对湿度和大气压力,计算空气饱和蒸汽压和单位体积干空气中的水蒸气体积占比/>进而计算单位质量燃煤实际助燃湿空气量Wair,act,wet和单位质量燃煤实际产生的湿烟气量Wflue,act,wet,具体而言:
Wflue,th,wet=α1×Qnet+β1,
Wair,th,dry=α2×Qnet+β2,
式中:
Wflue,th,wet为单位质量燃煤理论湿烟气量,单位为m3/kg,
Wair,th,dry为单位质量燃煤理论助燃干空气量,单位为m3/kg,
Wair,act,dry为单位质量燃煤实际助燃干空气量,单位为m3/kg,
Wflue,act,wet为单位质量燃煤湿基产生的湿烟气量,单位为m3/kg,
Qnet为煤低位热值,单位为MJ/kg,
O2为烟气运行湿基氧量,单位为%,
α1、α2、β1、β2均为常数,分别为0.2467、0.2496、0.718、0.3125,
为空气饱和蒸汽压,单位为Pa,
Tair为空气温度,单位为℃,
ф为空气相对湿度,单位为%,
Patm为大气压力,单位为Pa,
为单位体积干空气中的水蒸汽体积占比。
S3:根据工业分析(收到基灰、收到基水、低位热值)和元素分析(碳、氢、氧、氮、硫)结果,得到不同负荷下每小时燃煤总量的所需助燃干空气量设计值Wair,b和产生湿烟气量设计值Wflue,b;
计算不同负荷下每小时燃煤总量的所需助燃干空气量预测值Wair,h和产生湿烟气量预测值Wflue,h,具体而言:
Wair,h=Wair,act,dry×Wcoal,b×(1-γ)×1000,
Wflue,h=Wflue,act,wet×Wcoal,b×(1-γ)×1000,
式中:
Wair,h为每小时燃煤总量的所需助燃干空气量预测值,单位为m3/h,
Wflue,h为每小时燃煤总量的产生湿烟气量预测值,单位为m3/h,
Wair,act,dry为单位质量燃煤实际助燃干空气量,单位为m3/kg,
Wflue,act,wet为单位质量燃煤湿基产生湿烟气量预测值,单位为m3/kg,
Wcoal,b为设计给煤量,单位为t/h,
γ为入炉煤的未燃烧热损失比例,其值为0.42%;
计算不同负荷下每小时燃煤总量的所需助燃干空气量预测值Wair,h与其设计值Wair,b之间的相对偏差δ1、产生湿烟气量预测值Wflue,h与其设计值Wflue,b之间的相对偏差δ2,若相对偏差δ1和δ2均在-5%~+5%内,则拟合后的偏差在可接受范围内,具体而言:
式中:
Wair,h为每小时燃煤总量的所需助燃干空气量预测值,单位为m3/h,
Wflue,h为每小时燃煤总量的产生湿烟气量预测值,单位为m3/h,
Wair,b为每小时燃煤总量的所需助燃干空气量设计值,单位为m3/h,
Wflue,b为每小时燃煤总量的产生湿烟气量设计值,单位为m3/h。
S4:获取机组DCS系统显示的实际运行负荷Li和实际给煤量Wcoal,act,通过机组负荷L与设计给煤量Wcoal,b的对应函数关系计算实际运行负荷Li下的设计给煤量Wcoal,b,i,并计算实际给煤低位热值Qnet,act,具体而言:
式中:
Qnet,act为实际给煤低位热值,单位为MJ/kg,
Qnet,b为设计给煤低位热值,单位为MJ/kg,
Wcoal,act为实际运行负荷Li下的实际给煤量,单位为t/h,
Wcoal,b,i为实际运行负荷Li下的设计给煤量,单位为t/h。
S5:根据机组DCS系统显示的实际给煤量Wcoal,act、实际烟气运行湿基氧量O2,act以及实际给煤低位热值Qnet,act,计算出实际运行负荷Li下的进入锅炉的助燃干空气量BWair,act,dry和助燃湿空气量BWair,act,wet,锅炉省煤器出口湿烟气量BWflue,act,wet,具体而言:
式中:
BWair,act,dry为助燃干空气量,单位为m3/h,
BWair,act,wet为助燃湿空气量,单位为m3/h,
BWflue,act,wet为锅炉出口湿烟气量,单位为m3/h,
Wair,th,dry为单位质量燃煤理论助燃干空气量,单位为m3/kg,
Wflue,th,wet为单位质量燃煤理论湿烟气量,单位为m3/kg,
O2,act为实际烟气运行湿基氧量,单位为%,
为单位体积干空气中的水蒸汽体积占比。
Wcoal,act为实际运行负荷Li下的实际给煤量,单位为t/h。
以下具体给出实施例进行详细阐述:
选择某350MW的机组,计算负荷为300MW下的实时助燃干空气量和湿烟气量,具体包括如下步骤:
S1:在机组稳态运行模式下,分别获取负荷367.5MW、350MW、367.5MW、175MW、87.5MW下的设计给煤量(设计给煤低位热值Qnet,b为21.652MJ/kg)和设计烟气运行湿基氧量,如表1所示,则机组负荷L与设计给煤量Wcoal,b的具体函数关系为:
Wcoal,b=k×L+b
式中:
L为机组负荷,单位为MW;
Wcoal,b为设计给煤量,单位为t/h;
k和b均为常数,分别为0.3545和9.7697。
表1:不同负荷下的给煤量和烟气运行湿基氧量
S2:已知设计烟气运行湿基氧量O2,b为3.6%、设计给煤低位热值Qnet,b为21.652MJ/kg,则可计算得出Wflue,th,wet=6.059548m3/kg,Wair,th,dry=5.716839m3/kg,且通过拟合已有的500多组电站锅炉煤种样品的元素分析和工业分析数据得到的单位质量燃煤理论湿烟气量Wflue,th,wet和单位质量燃煤理论助燃干空气量Wair,th,dry方差分别为0.9763和0.9858,均大于0.9,因此拟合的准确性高,进而计算出单位质量燃煤实际助燃干空气量为Wair,act,wet=6.97741m3/kg,已知空气温度为20℃,相对湿度为55%,大气压力为101000Pa,计算得到单位体积干空气中的水蒸汽占比为
S3:当机组负荷为367.5MW时,根据工业分析和元素分析结果,得到机组负荷367.5MW下每小时燃煤总量的所需助燃干空气量设计值Wair,b为952746m3/h、产生湿烟气量Wflue,b设计值为1030231m3/h;已知设计给煤量Wcoal,b为139.23t/h,则每小时燃煤总量的所需助燃干空气量预测值Wair,h=6.97741×139.23×(1-0.42%)×1000=967385m3/h,每小时燃煤总量的产生湿烟气量预测值Wflue,h=7.3552×139.23×(1-0.42%)×1000=1019768m3/h;则助燃干空气量相对偏差产生湿烟气量相对偏差/>
其他负荷下的计算过程与367.5MW负荷下的计算过程相同,在此不再赘述,具体计算结果如表2中所示。
根据表2中的计算结果可知,不同负荷下每小时燃煤总量的所需助燃干空气量预测值Wair,h与其设计值Wair,b之间的相对偏差为1.54%~2.61%,产生湿烟气量预测值Wflue,h与其设计值Wflue,b之间的相对偏差为-0.1.02%~+0.63%,相对偏差均在-5%~+5%之间,则由此可判断拟合后的偏差在可接受范围内。
表2:每小时燃煤总量的所需助燃干空气量和产生烟气量的设计值与预测值
S4:实际运行负荷Li为300MW,机组DCS系统显示的实际给煤量Wcoal,act为121.5t/h,通过机组负荷L与设计给煤量Wcoal,b的对应函数关系,得到实际运行负荷300MW下的设计给煤量Wcoal,b,i=0.3545×300+9.769=116.12t/h,则可计算出实际给煤低位热值
S5:实际运行负荷Li为300MW时,获取实际运行湿基氧量O2,act为3.2%,
则可计算出进入锅炉的助燃干空气量为:
计算出锅炉省煤器出口湿烟气量为:
上述实施例只为说明本发明的技术构思及特点,其目的在于让熟悉此项技术的人士能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡根据本发明精神实质所作的等效变化或修饰,都应涵盖在本发明的保护范围之内。
Claims (8)
1.一种助燃干空气量和湿烟气量的实时计算方法,其特征在于:包括如下步骤:
S1:在机组稳态运行模式下,获取不同负荷下的设计给煤量和烟气运行湿基氧量,进而确定机组负荷L与设计给煤量Wcoal,b的对应函数关系、机组负荷L与设计烟气运行湿基氧量O2,b的对应函数关系,
S2:采用回归分析方法,拟合得到基于低位热值的单位质量燃煤理论湿烟气量Wflue,th,wet和单位质量燃煤理论助燃干空气量Wair,th,dry的回归函数,判断回归函数的准确性,计算得到单位质量燃煤实际助燃干空气量Wair,act,dry,根据空气温度、相对湿度和大气压力,计算空气饱和蒸汽压和单位体积干空气中的水蒸气体积占比/>进而计算单位质量燃煤实际助燃湿空气量Wair,act,wet和单位质量燃煤实际产生的湿烟气量Wflue,act,wet,
所述的空气饱和蒸气压单位体积干空气中的水蒸汽体积占比为/>单位质量燃煤实际助燃干空气量Wair,act,dry、单位质量燃煤实际助燃湿空气量Wair,act,wet、单位质量燃煤实际产生的湿烟气量Wflue,act,wet分别为:
式中:
为空气饱和蒸汽压,单位为Pa,
Tair为空气温度,单位为℃,
φ为空气相对湿度,单位为%,
Patm为大气压力,单位为Pa,Wair,act,dry为单位质量燃煤实际助燃干空气量,单位为m3/kg,
Wair,th,dry为单位质量燃煤理论助燃干空气量,单位为m3/kg,
Wflue,th,wet为单位质量燃煤理论湿烟气量,单位为m3/kg,
Wflue,act,Wet为单位质量燃煤湿基产生的湿烟气量,单位为m3/kg,
O2为烟气运行湿基氧量,单位为%,
为每体积干空气中的水蒸汽体积占比,
S3:计算不同负荷下每小时燃煤总量的所需助燃干空气量预测值Wair,h和产生湿烟气量预测值Wflue,h,通过比较不同负荷下每小时燃煤总量的所需助燃干空气量预测值Wair,h与其设计值Wair,b、产生湿烟气量预测值Wflue,h与其设计值Wflue,b,验证相对偏差δ是否在可接受范围内,所述的每小时燃煤总量的所需助燃干空气量预测值Wair,h和产生湿烟气量预测值Wflue,h分别为:
Wair,h=Wair,act,dry×Wcoal,b×(1-γ)×1000,
Wflue,h=Wflue,act,wet×Wcoal,b×(1-γ)×1000,
式中:
Wair,h为每小时燃煤总量的所需助燃干空气量预测值,单位为m3/h,
Wflue,h为每小时燃煤总量的产生湿烟气量预测值,单位为m3/h,
Wair,act,dry为单位质量燃煤实际助燃干空气量,单位为m3/kg,
Wflue,act,wet为单位质量燃煤实际产生的湿烟气量,单位为m3/kg,
Wcoal,b为设计给煤量,单位为t/h,
γ为入炉煤的未燃烧损失比例,
S4:获取机组的实际运行负荷Li和实际给煤量Wcoal,act,通过机组负荷L与设计给煤量Wcoal,b的对应函数关系计算实际运行负荷Li下的设计给煤量Wcoal,b,i,并计算实际给煤低位热值Qnet,act,
S5:根据机组的实际给煤量Wcoal,act、实际运行湿基氧量02,act以及实际给煤低位热值Qnet,act,计算出实时进入锅炉的助燃干空气量BWair,act,dry和助燃湿空气量BWair,act,wet,锅炉省煤器出口湿烟气量BWflue,act,wet。
2.根据权利要求1所述的助燃干空气量和湿烟气量的实时计算方法,其特征在于:所述的机组负荷L与设计给煤量Wcoal,b的对应函数关系、机组负荷L与设计烟气运行湿基氧量O2,b的对应函数关系分别为:
Wcoal,b=f(L,Qnet,b),
O2,b=g(L),
式中:
L为机组负荷,单位为MW,
Wcoal,b为设计给煤量,单位为t/h,
Qnet,b为设计煤低位热值,单位为MJ/kg,
O2,b为设计烟气运行湿基氧量,单位为%。
3.根据权利要求1所述的助燃干空气量和湿烟气量的实时计算方法,其特征在于:所述的单位质量燃煤理论湿烟气量Wflue,th,wet和所述的单位质量燃煤理论助燃干空气量Wair,th,dry分别为:
Wflue,th,wet=α1×Qnet+β1,
Wair,th,dry=α2×Qnet+β2,
式中:
Wflue,th,wet为单位质量燃煤理论湿烟气量,单位为m3/kg,
Wair,th,dry为单位质量燃煤理论助燃干空气量,单位为m3/kg,
Qnet为煤低位热值,单位为MJ/kg,
α1、α2、β1、β2均为常数。
4.根据权利要求1所述的助燃干空气量和湿烟气量的实时计算方法,其特征在于:所述的回归函数的准确性通过方差进行判断,所述的方差应大于0.9。
5.根据权利要求1所述的助燃干空气量和湿烟气量的实时计算方法,其特征在于:所述的相对偏差δ的可接受范围为-5%~+5%,所述的相对偏差δ包括每小时燃煤总量的所需助燃干空气量预测值Wair,h与其设计值Wair,b的相对偏差δ1、产生湿烟气量预测值Wflue,h与其设计值
Wflue,b的相对偏差δ2,所述的相对偏差δ1和所述的相对偏差δ2分别为:
式中:
Wair,h为每小时燃煤总量的所需助燃干空气量预测值,单位为m3/h,
Wflue,h为每小时燃煤总量的产生湿烟气量预测值,单位为m3/h,
Wair,b为每小时燃煤总量的所需助燃干空气量设计值,单位为m3/h,
Wflue,b为每小时燃煤总量的产生湿烟气量设计值,单位为m3/h。
6.根据权利要求1所述的助燃干空气量和湿烟气量的实时计算方法,其特征在于:所述的实际给煤低位热值Qnet,act为:
式中:
Qnet,act为实际给煤低位热值,单位为MJ/kg,
Qnet,b为设计给煤低位热值,单位为MJ/kg,
Wcoal,act为实际运行负荷Li下的实际给煤量,单位为t/h,
Wcoal,b,i为实际运行负荷Li下的设计给煤量,单位为t/h。
7.根据权利要求1所述的助燃干空气量和湿烟气量的实时计算方法,其特征在于:所述的助燃干空气量BWair,act,dry和助燃湿空气量BWair,act,wet分别为:
式中:
BWair,act,dry为助燃干空气量,单位为m3/h,
BWair,act,wet为助燃湿空气量,单位为m3/h,
Wflue,th,wet为单位质量燃煤理论湿烟气量,单位为m3/kg,
Wair,th,dry为单位质量燃煤理论助燃干空气量,单位为m3/kg,
Wcoal,act为实际运行负荷Li下的实际给煤量,单位为t/h,
O2,act为实际烟气运行湿基氧量,单位为%,
为单位体积干空气中的水蒸汽体积占比。
8.根据权利要求1所述的助燃干空气量和湿烟气量的实时计算方法,其特征在于:所述的出口湿烟气量BWflue,act,wet为:
式中:
BWflue,act,wet为出口湿烟气量,单位为m3/h,
Wflue,th,wet为单位质量燃煤理论湿烟气量,单位为m3/kg,Wair,th,dry为单位质量燃煤理论助燃干空气量,单位为m3/kg,Weoal,act为实际运行负荷Li下的实际给煤量,单位为t/h,
O2,act为实际烟气运行湿基氧量,单位为%,为单位体积干空气中的水蒸汽体积占比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110598225.1A CN113283091B (zh) | 2021-05-31 | 2021-05-31 | 一种助燃干空气量和湿烟气量的实时计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110598225.1A CN113283091B (zh) | 2021-05-31 | 2021-05-31 | 一种助燃干空气量和湿烟气量的实时计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113283091A CN113283091A (zh) | 2021-08-20 |
CN113283091B true CN113283091B (zh) | 2023-11-17 |
Family
ID=77282547
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110598225.1A Active CN113283091B (zh) | 2021-05-31 | 2021-05-31 | 一种助燃干空气量和湿烟气量的实时计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113283091B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109655489A (zh) * | 2019-01-23 | 2019-04-19 | 国电科学技术研究院有限公司 | 一种燃煤机组标准煤耗的计量方法及装置 |
CN109712676A (zh) * | 2018-12-17 | 2019-05-03 | 江苏海事职业技术学院 | 基于燃烧系统运行参数的煤气热值软测量方法 |
CN110162918A (zh) * | 2019-05-31 | 2019-08-23 | 上海电力学院 | 一种高炉煤气燃气机组直流锅炉效率的获取方法和系统 |
CN110345507A (zh) * | 2019-06-11 | 2019-10-18 | 张家港优赛电子商务有限公司 | 一种燃煤锅炉燃烧自动优化控制装置 |
CN111396919A (zh) * | 2020-04-20 | 2020-07-10 | 西安热工研究院有限公司 | 一种基于煤质实时在线分析的单元机组控制系统及方法 |
EP3730938A1 (en) * | 2019-04-23 | 2020-10-28 | Nestlé Skin Health SA | New method to visualize distribution profile of penetration of compound(s) through the nail |
CN112131517A (zh) * | 2020-09-02 | 2020-12-25 | 苏州西热节能环保技术有限公司 | 一种垃圾焚烧电厂入炉垃圾低位热值的测算方法 |
-
2021
- 2021-05-31 CN CN202110598225.1A patent/CN113283091B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109712676A (zh) * | 2018-12-17 | 2019-05-03 | 江苏海事职业技术学院 | 基于燃烧系统运行参数的煤气热值软测量方法 |
CN109655489A (zh) * | 2019-01-23 | 2019-04-19 | 国电科学技术研究院有限公司 | 一种燃煤机组标准煤耗的计量方法及装置 |
EP3730938A1 (en) * | 2019-04-23 | 2020-10-28 | Nestlé Skin Health SA | New method to visualize distribution profile of penetration of compound(s) through the nail |
CN110162918A (zh) * | 2019-05-31 | 2019-08-23 | 上海电力学院 | 一种高炉煤气燃气机组直流锅炉效率的获取方法和系统 |
CN110345507A (zh) * | 2019-06-11 | 2019-10-18 | 张家港优赛电子商务有限公司 | 一种燃煤锅炉燃烧自动优化控制装置 |
CN111396919A (zh) * | 2020-04-20 | 2020-07-10 | 西安热工研究院有限公司 | 一种基于煤质实时在线分析的单元机组控制系统及方法 |
CN112131517A (zh) * | 2020-09-02 | 2020-12-25 | 苏州西热节能环保技术有限公司 | 一种垃圾焚烧电厂入炉垃圾低位热值的测算方法 |
Non-Patent Citations (1)
Title |
---|
水泥窑富氧燃烧技术;王广东等;《节能减排》;第79-81页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113283091A (zh) | 2021-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Energy-exergy analysis and energy efficiency improvement of coal-fired industrial boilers based on thermal test data | |
CN102012968B (zh) | 煤粉锅炉热效率和煤质数据实时监测方法 | |
CN104008297B (zh) | 煤粉与高炉煤气混烧锅炉的热效率计算方法 | |
CN105181926A (zh) | 基于热平衡的煤粉掺烧煤气锅炉燃煤热值软测量方法 | |
CN109655488B (zh) | 基于混合煤气预热燃烧的煤气热值软测量方法 | |
CN113341713B (zh) | 一种燃煤机组煤风同步动态协控方法 | |
CN112664975B (zh) | 一种适用于煤粉锅炉的风量控制方法 | |
CN112710487B (zh) | 一种脱硫废水干燥系统对锅炉热效率影响的测评方法 | |
CN109519960B (zh) | 基于含氧量和飞灰含碳量在线监测的煤粉炉燃烧调控方法 | |
CN104598761A (zh) | 混烧锅炉运行参数变化对机组发电煤耗影响的分析方法 | |
CN111931346A (zh) | 一种燃煤锅炉入炉煤低位热值的实时计算方法及应用 | |
CN110864855A (zh) | 一种在中低负荷下测量电站锅炉空预器漏风率的装置及方法 | |
CN115187012A (zh) | 一种基于日数据验证燃煤电厂燃烧产生的碳排放方法 | |
CN113283091B (zh) | 一种助燃干空气量和湿烟气量的实时计算方法 | |
CN112580890A (zh) | 用于发电的锅炉混煤掺烧的发电可变成本预测方法及系统 | |
CN112417684B (zh) | 一种燃煤机组煤质成分动态计算方法 | |
CN109635464A (zh) | 一种钢厂煤气含硫量软测量方法 | |
CN110222351A (zh) | 煤气混烧锅炉入炉煤气热值与锅炉热效率在线监测方法 | |
CN109580711B (zh) | 高炉煤气与转炉煤气混烧条件下的煤气热值软测量方法 | |
CN112348703A (zh) | 一种基于供电煤耗最低的最优运行氧量简化分析方法 | |
Lala et al. | Theoretical Air Requirement and Combustion Flue Gases Analysis for Indigenous Biomass Combustion | |
CN104614197A (zh) | 煤粉与高炉煤气混烧锅炉热效率的在线监测方法 | |
CN114091275B (zh) | 缺空预器工质侧进/出口参数的锅炉热效率在线计算方法 | |
CN109489766A (zh) | 一种火力发电机组燃料燃烧碳氧化因子在线计量方法 | |
Gordin et al. | To calculate heat losses with flue gases of coal-fired boilers |
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 |