CN113378394B - 一种基于古尔维奇热平衡的智能吹灰算法 - Google Patents

一种基于古尔维奇热平衡的智能吹灰算法 Download PDF

Info

Publication number
CN113378394B
CN113378394B CN202110681637.1A CN202110681637A CN113378394B CN 113378394 B CN113378394 B CN 113378394B CN 202110681637 A CN202110681637 A CN 202110681637A CN 113378394 B CN113378394 B CN 113378394B
Authority
CN
China
Prior art keywords
temperature
flue gas
inlet
working medium
outlet
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
CN202110681637.1A
Other languages
English (en)
Other versions
CN113378394A (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.)
Zhongnan Electric Power Test and Research Institute of China Datang Group Science and Technology Research Institute Co Ltd
Original Assignee
Zhongnan Electric Power Test and Research Institute of China Datang Group Science and Technology Research Institute Co Ltd
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 Zhongnan Electric Power Test and Research Institute of China Datang Group Science and Technology Research Institute Co Ltd filed Critical Zhongnan Electric Power Test and Research Institute of China Datang Group Science and Technology Research Institute Co Ltd
Priority to CN202110681637.1A priority Critical patent/CN113378394B/zh
Publication of CN113378394A publication Critical patent/CN113378394A/zh
Application granted granted Critical
Publication of CN113378394B publication Critical patent/CN113378394B/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
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

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)
  • Control Of Steam Boilers And Waste-Gas Boilers (AREA)

Abstract

本发明涉及一种基于古尔维奇热平衡的智能吹灰算法,基于局部热量平衡原理的软测量技术,通过迭代算法,在已知出口烟气温度、进口工质温度、出口工质温度、及工质质量流量的情况下,当吸热量与放热量的校核误差<2%时,推算得到烟气进口温度,进而计算得到实际传热系数、理论传热系数、清洁因子。该算法解决了测点不全等问题,节约了成本。该算法能有效监测锅炉受热面的污染状态,确定最佳吹灰时机,进而维持受热面清洁度,提高锅炉运行效率。同时也节约了吹灰介质,减少了受热面管壁的腐蚀和磨损,间接延长了锅炉各受热面的使用效率。

Description

一种基于古尔维奇热平衡的智能吹灰算法
技术领域
本发明涉及一种基于古尔维奇热平衡的智能吹灰算法,属于智能算法领域。
背景技术
电能是现在人类文明中最重要的能源之一。现阶段我国电能的生产很大部分都来自于燃煤发电,然而燃煤在燃烧过程中,燃煤发电机组在运行过程中,燃煤所含的灰分会造成锅炉受热面积灰和结渣,使受热面的传热系数下降,影响了受热面的传热效率,使得主蒸汽、再热蒸汽温度下降,而排烟温度升高,从而使锅炉效率下降,造成能源的浪费,严重时甚至会危害机组的安全运行。因此,为了达到煤燃烧能量的利用的最大化,必须对锅炉管道的受热面进行及时的清扫吹灰。制定最合理、最科学的吹灰优化策略,才能保证一段时间内受热面传热效率的最大化,实现节能减排和机组的安全运行。
目前我国大型电站锅炉吹灰操作一般采用定时定量的程序控制方式。由于这种运行方式是在不了解受热面实际灰污状态的情况下进行的,不可避免地会产生吹灰不足或吹灰过度,且吹灰时间、部位大多凭经验决定,吹灰过程存在人为的主观性,这就造成了更多煤耗。现阶段节能减排给吹灰优化提出了更高的要求,只有对未来状态进行预测,并且提前准备,才能让燃煤电站有更好经济收益。
发明内容
针对现有技术的不足,本发明的目的是提供一种基于古尔维奇热平衡的智能吹灰算法。
为了实现上述目的,本发明所采用的技术方案是:
一种基于古尔维奇热平衡的智能吹灰算法,包括以下步骤:
(1)计算理论传热系数K0
Figure GDA0004085273260000011
其中:
a1=ad+af
Figure GDA0004085273260000012
Figure GDA0004085273260000021
Figure GDA0004085273260000022
式中:a1为烟气对管壁的放热系数;
a2为工质对管壁的放热系数;
ad为表面放热系数;
af为辐射传热系数;
Cz、Cs分别为气流流动方向上纵向和横向管排修正系数;
λ为烟气平均温度下的导热系数,W/(m·℃);
d为管子外径,m;
Re、Pr分别为烟气的雷洛数和普朗特数;
Cl为冲刷受热面相对长度的修正系数;
Ct为流体温度和管壁温度Tgb的修正系数。
σ0为玻尔兹曼常数;
εs为***黑度;
Tgb、T分别为管壁温度和***温度,K;
(2)通过古尔维奇热平衡智能算法确定换热器进口烟气温度;
通过出口烟气温度、进口工质温度、出口工质温度,工质质量流量,假设的进口烟气温度值,其中各温度可换算为对应的焓值,计算得到工质吸热量、烟气放热量,通过迭代算法,基于局部热量平衡原理,当吸热量与放热量的校核误差<2%,此时假定的烟气进口温度即为最终确定的进口烟气温度;
(3)计算实际传热系数Ksj
Figure GDA0004085273260000023
其中:
Qz=Dz(h”-h')
Figure GDA0004085273260000024
式中:Qz为工质吸热量,kJ/kg;
Dz为进口工质质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
A为传热面积,m2
Δtmax和Δtmin分别代表大温压和小温压;
其中,逆流时,大温压Δtmax=J′-t″;顺流时,大温压Δtmax=J′-t′;
逆流时,小温压Δtmin=J″-t′;顺流时,小温压Δtmin=J″-t″;
J′:烟气进口温度,℃;t″:工质出口温度,℃;J″:烟气出口温度,℃;t′:工质进口温度,℃;
(4)计算清洁因子CF
Figure GDA0004085273260000031
当CF=1时,受热面处于理想的洁净状态,即光滑管子的洁净因子;CF小于1则表示受热面受到了灰污的污染,越小则污染越严重;
步骤(2)的具体方法为:
1)省煤器进口烟气温度的算法
省煤器通过烟气挡板分为过热器侧和再热器侧;
①再热器侧省煤器
烟气温度和烟气焓值是一一对应关系,已知温度通过查表可以得出焓值,已知焓值可以求出温度;假设再热器侧省煤器烟气进口温度为J′;
再热器侧省煤器工质吸热量:Qz=Dz(h”-h')
再热器侧省煤器烟气放热量:Qy=mz(H'-H”)
其中:Dz为省煤器进口工质的质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
mz为再热器侧燃煤量,kg/s;
H'、H”分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qz-Qy
如果Qz-Qy<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
②过热器侧省煤器
假设过热器侧省煤器烟气进口温度为J′;
过热器侧省煤器工质吸热量:Qz=Dz(h”-h')
过热器侧省煤器烟气放热量:Qy=mg(H'-H”)
其中:Dz为省煤器进口工质的质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
mg为过热器侧燃煤量,kg/s;
H'、H”分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qz-Qy
如果Qz-Qy<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
2)低温过热器水平段进口烟气温度的算法
低温过热器水平段烟气出口温度=过热器侧省煤器进口烟气温度,假设低温过热器水平段烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=mg(H'-H”)
其中:Dz为低温过热器进口工质的质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
mg:过热器侧燃煤量,kg/s;
H'、H”分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.02Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
3)低温过热器垂直段进口烟气温度的算法
低温过热器垂直段烟气出口温度=低温过热器水平段烟气进口温度,假设低温过热器垂直段烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中,Dz:低温过热器进口工质的质量流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.05Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
4)低温再热器水平段进口烟气温度的算法
低温再热器水平段烟气出口温度=再热器侧省煤器烟气进口温度,假设低温再热器水平段烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气的放热量:Qy=mz(H'-H”)
其中,Dz:低温再热器进口工质的质量流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
mz:再热器侧燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.015Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
5)低温再热器垂直段进口烟气温度的算法
低温再热器垂直段烟气出口温度=低温再热器水平段烟气进口温度,假设低温再热器垂直段烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中,Dz:低温再热器进口工质的质量流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.04Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
6)高温再热器进口烟气温度的算法
高温再热器烟气出口温度=低温再热器垂直段烟气进口温度,假设高温再热器烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中:Dz:再热器蒸汽流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.06Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
7)高温过热器进口烟气温度的算法
高温过热器烟气出口温度=高温再热器烟气进口温度,假设高温过热器烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中,Dz:主蒸汽流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.06Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
8)屏式过热器进口烟气温度的算法
屏式过热器烟气出口温度=高温过热器烟气进口温度,假设屏式过热器烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
屏吸收的总热量:Qp=Qd p+Qf p
受热面的对流吸热量:Qd p=Qy-Qd fj
附加受热面对流吸热量:Qd fj=0.15Qz
屏吸收的辐射热量:Qf p=Qf-Qf fj
屏吸收的炉内直接辐射:Qf=Qf′-Qf
附加受热面获得的炉内直接辐射:
Figure GDA0004085273260000071
落入屏区的炉内直接辐射:
Figure GDA0004085273260000072
穿透屏区的炉内直接辐射:
Figure GDA0004085273260000081
落入屏区的炉膛出口热流:
Figure GDA0004085273260000082
炉膛的吸热量:
Figure GDA0004085273260000083
式中,Dz:进入屏式过热器的蒸汽流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:烟气进口和出口的焓值,kJ/kg;
F”:炉膛出口烟窗面积,㎡;
a:烟气黑度;
x'p:出口截面的角系数;
Ach:屏出口烟窗的面积,㎡;
Tpj:烟气的平均温度,K;
Bj:总燃煤量,kg/s;
β:炉膛与屏相互换热系数;
η:沿炉膛高度不均匀系数;
H1:炉膛有效辐射受热面积,㎡;
Figure GDA0004085273260000084
保温系数;
Q1:低位发热量,kJ/kg;
I1”:炉膛出口烟焓,kJ/kg;
Afj:附加受热面积,㎡;
S:总受热面积,㎡;
计算校核误差:Qp-Qz
如果Qp-Qz<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差。
本发明有益效果:
清洁因子CF(Cleanliness Factor)可以表征对流受热面的污染状态,是受热面洁净程度的标志,定义为实际传热系数和理论传热系数的比值。通常通过计算清洁因子来监测各受热面的积灰状况,然而计算清洁因子需要知道各受热面的烟气进出口温度,而很多电厂都没有测点,使清洁因子的计算受到了限制。
本发明利用古尔维奇热平衡的智能吹灰算法,基于局部热量平衡原理的软测量技术,通过迭代算法,在已知出口烟气温度、进口工质温度、出口工质温度、及工质质量流量的情况下,当吸热量与放热量的校核误差<2%时,推算得到烟气进口温度,进而计算得到实际传热系数、理论传热系数、清洁因子。该算法解决了测点不全等问题,节约了成本。
本发明基于古尔维奇热平衡的智能吹灰算法能有效监测锅炉受热面的污染状态,确定最佳吹灰时机,进而维持受热面清洁度,提高锅炉运行效率。同时也节约了吹灰介质,减少了受热面管壁的腐蚀和磨损,间接延长了锅炉各受热面的使用效率。
附图说明
图1为省煤器(再热器侧)进口烟气温度算法框图;
图2为省煤器(过热器侧)进口烟气温度算法框图;
图3为低温过热器水平段进口烟气温度算法框图;
图4为低温过热器垂段进口烟气温度算法框图;
图5为低温再热器水平段进口烟气温度算法框图;
图6为低温再热器垂直段进口烟气温度算法框图;
图7为高温再热器进口烟气温度算法框图;
图8为高温过热器进口烟气温度算法框图;
图9为屏式过热器进口烟气温度算法框图;
图10为本发明古尔维奇热平衡智能算法流程图。
具体实施方式
以下结合实施例对本发明的具体实施方式作进一步详细说明。
本发明公式中涉及的热量及焓值的单位均为kJ/kg。
实施例
一种基于古尔维奇热平衡的智能吹灰算法,包括以下步骤:
(1)计算理论传热系数K0
Figure GDA0004085273260000101
其中:
a1=ad+af
Figure GDA0004085273260000102
Figure GDA0004085273260000103
Figure GDA0004085273260000104
式中:a1为烟气对管壁的放热系数;
a2为工质对管壁的放热系数;
ad为表面放热系数;
af为辐射传热系数;
Cz、Cs分别为气流流动方向上纵向和横向管排修正系数;
λ为烟气平均温度下的导热系数,W/(m·℃);
d为管子外径,m;
Re、Pr分别为烟气的雷洛数和普朗特数;
Cl为冲刷受热面相对长度的修正系数;
Ct为流体温度和管壁温度Tgb修正系数。
σ0为玻尔兹曼常数;
εs为***黑度;
Tgb、T分别为管壁温度和***温度,K;
(2)通过古尔维奇热平衡智能算法确定换热器进口烟气温度;
通过出口烟气温度、进口工质温度、出口工质温度,工质质量流量,假设的进口烟气温度值,其中各温度可换算为对应的焓值,计算得到工质吸热量、烟气放热量,通过迭代算法,基于局部热量平衡原理,当吸热量与放热量的校核误差<2%,此时假定的烟气进口温度即为最终确定的进口烟气温度;
(3)计算实际传热系数Ksj
Figure GDA0004085273260000111
其中:
Qz=Dz(h”-h')
Figure GDA0004085273260000112
式中:Qz为工质吸热量,kJ/kg;
Dz为进口工质质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
A为传热面积,m2
Δtmax和Δtmin分别代表大温压和小温压;
其中,逆流时,大温压Δtmax=J′-t″;顺流时,大温压Δtmax=J′-t′;
逆流时,小温压Δtmin=J″-t′;顺流时,小温压Δtmin=J″-t″;
J′:烟气进口温度,℃;
t″:工质出口温度,℃;
J″:烟气出口温度,℃;
t′:工质进口温度,℃;
(4)计算清洁因子CF
Figure GDA0004085273260000113
当CF=1时,受热面处于理想的洁净状态,即光滑管子的洁净因子;CF小于1则表示受热面受到了灰污的污染,越小则污染越严重。
步骤(2)通过古尔维奇热平衡智能算法确定换热器进口烟气温度的具体方法如下:
1)省煤器进口烟气温度的算法
省煤器通过烟气挡板分为过热器侧和再热器侧。
①省煤器(再热器侧)
烟气温度和烟气焓值是一一对应关系,已知温度通过查表可以得出焓值,已知焓值可以求出温度。假设省煤器(再热器侧)烟气进口温度为J′。
省煤器(再热器侧)工质吸热量:Qz=Dz(h”-h')
省煤器(再热器侧)烟气放热量:Qy=mz(H'-H”)
其中:Dz为省煤器进口工质的质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
mz为再热器侧燃煤量,kg/s;
H'、H”分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qz-Qy
如果Qz-Qy<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差(如图1所示)。
②省煤器(过热器侧)
假设省煤器(过热器侧)烟气进口温度为J′。
省煤器(过热器侧)工质吸热量:Qz=Dz(h”-h')
省煤器(过热器侧)烟气放热量:Qy=mg(H'-H”)
其中:Dz为省煤器进口工质的质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
mg为过热器侧燃煤量,kg/s;
H'、H”分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qz-Qy
如果Qz-Qy<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差(如图2所示)。2)低温过热器水平段进口烟气温度的算法
低温过热器水平段烟气出口温度=省煤器(过热器侧)进口烟气温度,假设低温过热器水平段烟气进口温度为J′。
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=mg(H'-H”)
其中:Dz为低温过热器进口工质的质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
mg:过热器侧燃煤量,kg/s;
H'、H”分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.02Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差(如图3所示)。3)低温过热器垂直段进口烟气温度的算法
低温过热器垂直段烟气出口温度=低温过热器水平段烟气进口温度,假设低温过热器垂直段烟气进口温度为J′。
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中,Dz:低温过热器进口工质的质量流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.05Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差(如图4所示)。4)低温再热器水平段进口烟气温度的算法
低温再热器水平段烟气出口温度=省煤器(再热器侧)烟气进口温度,假设低温再热器水平段烟气进口温度为J′。
工质吸热量:Qz=Dz(h”-h')
烟气的放热量:Qy=mz(H'-H”)
其中,Dz:低温再热器进口工质的质量流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
mz:再热器侧燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.015Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差(如图5所示)。5)低温再热器垂直段进口烟气温度的算法
低温再热器垂直段烟气出口温度=低温再热器水平段烟气进口温度,假设低温再热器垂直段烟气进口温度为J′。
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中,Dz:低温再热器进口工质的质量流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.04Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差(如图6所示)。6)高温再热器进口烟气温度的算法
高温再热器烟气出口温度=低温再热器垂直段烟气进口温度,假设高温再热器烟气进口温度为J′。
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中:Dz:再热器蒸汽流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.06Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差(如图7所示)。7)高温过热器进口烟气温度的算法
高温过热器烟气出口温度=高温再热器烟气进口温度,假设高温过热器烟气进口温度为J′。
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中,Dz:主蒸汽流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.06Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差(如图8所示)。8)屏式过热器进口烟气温度的算法
屏式过热器烟气出口温度=高温过热器烟气进口温度,假设屏式过热器烟气进口温度为J′。
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
屏吸收的总热量:Qp=Qd p+Qf p
受热面的对流吸热量:Qd p=Qy-Qd fj
附加受热面对流吸热量:Qd fj=0.15Qz
屏吸收的辐射热量:Qf p=Qf-Qf fj
屏吸收的炉内直接辐射:Qf=Qf′-Qf
附加受热面获得的炉内直接辐射:
Figure GDA0004085273260000161
落入屏区的炉内直接辐射:
Figure GDA0004085273260000162
穿透屏区的炉内直接辐射:
Figure GDA0004085273260000163
落入屏区的炉膛出口热流:
Figure GDA0004085273260000164
炉膛的吸热量:
Figure GDA0004085273260000165
式中,Dz:进入屏式过热器的蒸汽流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:烟气进口和出口的焓值,kJ/kg;
F”:炉膛出口烟窗面积,㎡;
a:烟气黑度;
x'p:出口截面的角系数;
Ach:屏出口烟窗的面积,㎡;
Tpj:烟气的平均温度,K;
Bj:总燃煤量,kg/s;
β:炉膛与屏相互换热系数;
η:沿炉膛高度不均匀系数;
H1:炉膛有效辐射受热面积,㎡;
Figure GDA0004085273260000171
保温系数;
Q1:低位发热量,kJ/kg;
I1”:炉膛出口烟焓,kJ/kg;
Afj:附加受热面积,㎡;
S:总受热面积,㎡;
计算校核误差:Qp-Qz
如果Qp-Qz<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差(如图9所示)。

Claims (1)

1.一种基于古尔维奇热平衡的智能吹灰算法,其特征在于,包括以下步骤:(1)计算理论传热系数K0
Figure FDA0004085273220000011
其中:
a1=ad+af
Figure FDA0004085273220000012
Figure FDA0004085273220000013
Figure FDA0004085273220000014
式中:a1为烟气对管壁的放热系数;
a2为工质对管壁的放热系数;
ad为表面放热系数;
af为辐射传热系数;
Cz、Cs分别为气流流动方向上纵向和横向管排修正系数;
λ为烟气平均温度下的导热系数,W/(m·℃);
d为管子外径,m;
Re、Pr分别为烟气的雷洛数和普朗特数;
Cl为冲刷受热面相对长度的修正系数;
Ct为流体温度和管壁温度Tgb的修正系数。
σ0为玻尔兹曼常数;
εs为***黑度;
Tgb、T分别为管壁温度和***温度,K;
(2)通过古尔维奇热平衡智能算法确定换热器进口烟气温度;
通过出口烟气温度、进口工质温度、出口工质温度,工质质量流量,假设的进口烟气温度值,其中各温度可换算为对应的焓值,计算得到工质吸热量、烟气放热量,通过迭代算法,基于局部热量平衡原理,当吸热量与放热量的校核误差<2%,此时假定的烟气进口温度即为最终确定的进口烟气温度;
(3)计算实际传热系数Ksj
Figure FDA0004085273220000021
其中:
Qz=Dz(h”-h')
Figure FDA0004085273220000022
式中:Qz为工质吸热量,kJ/kg;
Dz为进口工质质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
A为传热面积,m2
Δtmax和Δtmin分别代表大温压和小温压;
其中,逆流时,大温压Δtmax=J′-t″;顺流时,大温压Δtmax=J′-t′;
逆流时,小温压Δtmin=J″-t′;顺流时,小温压Δtmin=J″-t″;
J′:烟气进口温度,℃;t″:工质出口温度,℃;J″:烟气出口温度,℃;t′:工质进口温度,℃;
(4)计算清洁因子CF
Figure FDA0004085273220000023
当CF=1时,受热面处于理想的洁净状态,即光滑管子的洁净因子;CF小于1则表示受热面受到了灰污的污染,越小则污染越严重;
步骤(2)的具体方法为:
1)省煤器进口烟气温度的算法
省煤器通过烟气挡板分为过热器侧和再热器侧;
①再热器侧省煤器
烟气温度和烟气焓值是一一对应关系,已知温度通过查表可以得出焓值,已知焓值可以求出温度;假设再热器侧省煤器烟气进口温度为J′;
再热器侧省煤器工质吸热量:Qz=Dz(h”-h')
再热器侧省煤器烟气放热量:Qy=mz(H'-H”)
其中:Dz为省煤器进口工质的质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
mz为再热器侧燃煤量,kg/s;
H'、H”分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qz-Qy
如果Qz-Qy<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
②过热器侧省煤器
假设过热器侧省煤器烟气进口温度为J′;
过热器侧省煤器工质吸热量:Qz=Dz(h”-h')
过热器侧省煤器烟气放热量:Qy=mg(H'-H”)
其中:Dz为省煤器进口工质的质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
mg为过热器侧燃煤量,kg/s;
H'、H”分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qz-Qy
如果Qz-Qy<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
2)低温过热器水平段进口烟气温度的算法
低温过热器水平段烟气出口温度=过热器侧省煤器进口烟气温度,假设低温过热器水平段烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=mg(H'-H”)
其中:Dz为低温过热器进口工质的质量流量,kg/s;
h”、h'分别为工质出口和进口的焓值,kJ/kg;
mg:过热器侧燃煤量,kg/s;
H'、H”分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.02Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
3)低温过热器垂直段进口烟气温度的算法
低温过热器垂直段烟气出口温度=低温过热器水平段烟气进口温度,假设低温过热器垂直段烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中,Dz:低温过热器进口工质的质量流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.05Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
4)低温再热器水平段进口烟气温度的算法
低温再热器水平段烟气出口温度=再热器侧省煤器烟气进口温度,假设低温再热器水平段烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气的放热量:Qy=mz(H'-H”)
其中,Dz:低温再热器进口工质的质量流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
mz:再热器侧燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.015Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
5)低温再热器垂直段进口烟气温度的算法
低温再热器垂直段烟气出口温度=低温再热器水平段烟气进口温度,假设低温再热器垂直段烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中,Dz:低温再热器进口工质的质量流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.04Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
6)高温再热器进口烟气温度的算法
高温再热器烟气出口温度=低温再热器垂直段烟气进口温度,假设高温再热器烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中:Dz:再热器蒸汽流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.06Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
7)高温过热器进口烟气温度的算法
高温过热器烟气出口温度=高温再热器烟气进口温度,假设高温过热器烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
其中,Dz:主蒸汽流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:分别为烟气进口和出口的焓值,kJ/kg;
计算校核误差:Qy-Qz-Qfj
其中,Qfj为附加受热面对流吸热量,Qfj=0.06Qz
如果Qy-Qz-Qfj<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差;
8)屏式过热器进口烟气温度的算法
屏式过热器烟气出口温度=高温过热器烟气进口温度,假设屏式过热器烟气进口温度为J′;
工质吸热量:Qz=Dz(h”-h')
烟气放热量:Qy=m(H'-H”)
屏吸收的总热量:Qp=Qd p+Qf p
受热面的对流吸热量:Qd p=Qy-Qd fj
附加受热面对流吸热量:Qd fj=0.15Qz
屏吸收的辐射热量:Qf p=Qf-Qf fj
屏吸收的炉内直接辐射:Qf=Qf′-Qf
附加受热面获得的炉内直接辐射:
Figure FDA0004085273220000071
落入屏区的炉内直接辐射:
Figure FDA0004085273220000072
穿透屏区的炉内直接辐射:
Figure FDA0004085273220000073
落入屏区的炉膛出口热流:
Figure FDA0004085273220000074
炉膛的吸热量:
Figure FDA0004085273220000075
式中,Dz:进入屏式过热器的蒸汽流量,kg/s;
h”、h':分别为工质出口和进口的焓值,kJ/kg;
m:总燃煤量,kg/s;
H'、H”:烟气进口和出口的焓值,kJ/kg;
F”:炉膛出口烟窗面积,㎡;
a:烟气黑度;
x'p:出口截面的角系数;
Ach:屏出口烟窗的面积,㎡;
Tpj:烟气的平均温度,K;
Bj:总燃煤量,kg/s;
β:炉膛与屏相互换热系数;
η:沿炉膛高度不均匀系数;
H1:炉膛有效辐射受热面积,㎡;
Figure FDA0004085273220000076
保温系数;
Q1:低位发热量,kJ/kg;
I1”:炉膛出口烟焓,kJ/kg;
Afj:附加受热面积,㎡;
S:总受热面积,㎡;
计算校核误差:Qp-Qz
如果Qp-Qz<2%,则此时假定的烟气进口温度即为最终确定的进口烟气温度,否则重新假设烟气进口温度进行计算,直到符合计算误差。
CN202110681637.1A 2021-06-19 2021-06-19 一种基于古尔维奇热平衡的智能吹灰算法 Active CN113378394B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110681637.1A CN113378394B (zh) 2021-06-19 2021-06-19 一种基于古尔维奇热平衡的智能吹灰算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110681637.1A CN113378394B (zh) 2021-06-19 2021-06-19 一种基于古尔维奇热平衡的智能吹灰算法

Publications (2)

Publication Number Publication Date
CN113378394A CN113378394A (zh) 2021-09-10
CN113378394B true CN113378394B (zh) 2023-04-18

Family

ID=77577893

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110681637.1A Active CN113378394B (zh) 2021-06-19 2021-06-19 一种基于古尔维奇热平衡的智能吹灰算法

Country Status (1)

Country Link
CN (1) CN113378394B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05264005A (ja) * 1992-01-21 1993-10-12 Babcock Hitachi Kk ボイラ装置の演算装置
CN105276563A (zh) * 2015-11-11 2016-01-27 东南大学 一种基于炉膛实时结渣情况的炉膛出口烟温软测量方法
CN106322412A (zh) * 2016-08-30 2017-01-11 上海交通大学 基于二维寻优的燃煤机组对流受热面智能吹灰方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101034009A (zh) * 2007-04-09 2007-09-12 上海发电设备成套设计研究院 一种大型燃煤锅炉在线检测、吹灰优化节能方法
US20130192541A1 (en) * 2010-04-29 2013-08-01 Siemens Aktiengesellschaft Method and device for controlling the temperature of steam in a boiler
CN102734783B (zh) * 2012-06-12 2014-09-03 广东电网公司电力科学研究院 一种电站超临界锅炉各级受热面监测数据参数的校验方法
CN104896457A (zh) * 2015-06-03 2015-09-09 华北电力科学研究院有限责任公司 锅炉各级受热面烟气温度确定方法及装置
CN107525064A (zh) * 2017-08-29 2017-12-29 国网浙江省电力公司电力科学研究院 基于温压计算的煤粉锅炉省煤器吹灰***优化操作方法
CN109300512B (zh) * 2018-10-11 2022-02-18 中国石油化工股份有限公司 一种适用于克劳斯炉的整套热力计算方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05264005A (ja) * 1992-01-21 1993-10-12 Babcock Hitachi Kk ボイラ装置の演算装置
CN105276563A (zh) * 2015-11-11 2016-01-27 东南大学 一种基于炉膛实时结渣情况的炉膛出口烟温软测量方法
CN106322412A (zh) * 2016-08-30 2017-01-11 上海交通大学 基于二维寻优的燃煤机组对流受热面智能吹灰方法

Also Published As

Publication number Publication date
CN113378394A (zh) 2021-09-10

Similar Documents

Publication Publication Date Title
CN103968401B (zh) 低pm2.5排放的除尘-脱硫-余热利用一体化***
CN112283689B (zh) 燃煤电站锅炉受热面积灰在线监测***及其检测方法
CN115374636B (zh) 一种锅炉数值模拟与性能计算耦合的壁温预测方法
CN105091932A (zh) 一种三分仓回转式预热器分段积灰监测方法
CN103939941A (zh) 一种融入了不可逆热力学的锅炉燃烧优化方法
CN113378394B (zh) 一种基于古尔维奇热平衡的智能吹灰算法
CN112833409A (zh) 一种基于动态损失预测的炉膛吹灰优化方法
CN111637478A (zh) 一种炉膛吹灰方法
CN111578723A (zh) 一种冶金高温含尘废气处理***
CN211011433U (zh) 一种锅炉吹灰优化控制***
CN105299894B (zh) 一种燃气热水锅炉空气循环节能***
CN205482504U (zh) 一种可水平放置的氟塑料管式换热器结构
CN105352174B (zh) 一种大型燃气热水锅炉
CN112097287B (zh) 一种锅炉节能与烟气脱白***、工艺、应用
CN205208690U (zh) 一种燃气热水锅炉水路循环节能***
CN205156329U (zh) 一种大型燃气热水锅炉
CN205156347U (zh) 一种燃气热水锅炉空气循环节能***
CN111271718A (zh) 一种用于低低温省煤器利用烟气吹灰的装置和方法
CN205940178U (zh) 一种转底炉烟气余热利用***
CN220355387U (zh) 带有倾斜式蛇形翅片管的活性炭余热锅炉
CN205505051U (zh) 一种锅炉烟热回收***
CN110362873A (zh) 一种吹灰控制方法、***、介质及设备
Yuanhao et al. On-line monitoring of ash fouling for air pre-heaters in coal-fired power plant boiler
CN105299726A (zh) 一种燃气热水锅炉水路循环节能***
CN202844882U (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