CN111521430B - 一种余热锅炉性能试验方法 - Google Patents

一种余热锅炉性能试验方法 Download PDF

Info

Publication number
CN111521430B
CN111521430B CN202010520009.0A CN202010520009A CN111521430B CN 111521430 B CN111521430 B CN 111521430B CN 202010520009 A CN202010520009 A CN 202010520009A CN 111521430 B CN111521430 B CN 111521430B
Authority
CN
China
Prior art keywords
flue gas
pressure
steps
heat boiler
following
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
CN202010520009.0A
Other languages
English (en)
Other versions
CN111521430A (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.)
Electric Power Research Institute of State Grid Anhui Electric Power Co Ltd
Anhui Xinli Electric Technology Consulting Co Ltd
Original Assignee
Electric Power Research Institute of State Grid Anhui Electric Power Co Ltd
Anhui Xinli Electric Technology Consulting 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 Electric Power Research Institute of State Grid Anhui Electric Power Co Ltd, Anhui Xinli Electric Technology Consulting Co Ltd filed Critical Electric Power Research Institute of State Grid Anhui Electric Power Co Ltd
Priority to CN202010520009.0A priority Critical patent/CN111521430B/zh
Publication of CN111521430A publication Critical patent/CN111521430A/zh
Application granted granted Critical
Publication of CN111521430B publication Critical patent/CN111521430B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M99/00Subject matter not provided for in other groups of this subclass
    • G01M99/005Testing of complete machines, e.g. washing-machines or mobile phones

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Steam Boilers And Waste-Gas Boilers (AREA)

Abstract

本发明公开了一种余热锅炉性能试验方法,其特征是首先基于目前测量水平和仪器精度测量相关可测准的物理参数,然后根据这些准确的测量结果进行计算,分别计算出余热锅炉热效率,出力以及汽水侧和烟气侧压降。本发明在相关试验标准的基础上考虑到现场实际情况、测量仪器以及测量水平等因素,给出了一种实际可操作的试验方法,简化了试验过程,提高了试验结果的准确性和可信度。

Description

一种余热锅炉性能试验方法
技术领域
本发明涉及燃气蒸汽联合循环机组余热锅炉性能试验技术领域,更具体地说是一种余热锅炉性能试验方法。
背景技术
新能源发电负荷不稳定,随着时间发生变化,比如风电随着风速变化而变化,光伏发电随着日照的强度而发生变化。因此,新能源发电具有间歇性和波动性的特征。大规模快速增长的新能源发电并网后必然对电网的电压、电流和频率产生较大冲击,影响电网的电能质量。为了消除这些不利影响,需要增加快速调峰机组的容量。调峰性能好、运行灵活的燃气蒸汽联合机组将会用于适应装机规模和发电量不断增长的新能源发电。燃气蒸汽联合机组包括压气机、燃气轮机、余热锅炉、汽轮机和发电机等主要设备。随着燃气蒸汽联合机组的发展,余热锅炉作为余热利用的关键设备,其性能指标是比较重要的技术参数。在机组建设之后或者检修前后,需要对余热锅炉的性能进行测量试验,对新建的余热锅炉进行性能考核或者对已投产的余热锅炉的检修工作提供参考与建议,并将检修前与检修后的余热锅炉性能指标进行比较,以确定检修的效果。
目前,针对联合循环余热锅炉性能试验,国家能源局2015年发布了DL/T 1427-2015《联合循环余热锅炉性能试验规程》,规程对联合循环余热锅炉性能试验过程做了详细规定,从试验前的准备,试验执行过程以及试验报告的格式均做了详细介绍,并对测量的仪器和测量方法也做了规定,给出了试验数据处理与计算的过程,通过算例进行演示。尽管如此,在目前的测量水平和测量装备下,规程中给出的试验方法更多的停留在理论阶段,不具有可操作性,并不适用于各种类型的余热锅炉性能试验的现场实际情况。根据目前的测量水平,规程中有些参数并不能直接准确获得,比如:余热锅炉的蒸汽流量无法准确测量,同轴燃气蒸汽联合循环机组的燃气轮机做功发电功率与蒸汽轮机做功发电功率不能完全分开计量等。正是由于相当一些关键参数在实际现场的条件下无法直接准确测量,因此现有规程中从能量守恒角度提出的试验方法在实际应用过程中不具有可操作性。
发明内容
本发明是为避免上述现有技术所存的不足,提供一种余热锅炉性能试验方法,根据现场实际可准确测量的数据,并依据质量守恒的原则给出针对现场测量情况的计算方法,使得对于余热锅炉性能试验具有可操作性,试验结果更加准确与可信。
本发明为解决技术问题采用如下技术方案:
本发明余热锅炉性能试验方法的特点是按如下步骤进行:
步骤1、确定余热锅炉烟道中进口烟气组分和摩尔分数,所述进口烟气组分包括N2、O2、CO2、H2O、Ar和SO2,一一对应的摩尔分数为
Figure GDA0003194156170000021
XAr
Figure GDA0003194156170000022
无补燃余热锅炉烟道中出口烟气组分和摩尔分数与进口相同;
步骤2、计算获得余热锅炉进口烟气焓值Hin
确定进口烟气温度Tin,并根据进口烟气温度Tin计算获得进口烟气中各组分的焓值,包括:
进口烟气中N2的焓值
Figure GDA0003194156170000023
为:
Figure GDA0003194156170000024
进口烟气中O2的焓值
Figure GDA0003194156170000025
为:
Figure GDA0003194156170000026
进口烟气中CO2的焓值
Figure GDA0003194156170000027
为:
Figure GDA0003194156170000028
进口烟气中H2O的焓值
Figure GDA0003194156170000029
为:
Figure GDA00031941561700000210
进口烟气中Ar的焓值HAr为:
Figure GDA00031941561700000211
进口烟气中SO2的焓值
Figure GDA00031941561700000212
为:
Figure GDA00031941561700000213
其中,
Figure GDA00031941561700000214
aj,Ar
Figure GDA00031941561700000215
均为经验系数,j=1,2,…,6;
计算获得余热锅炉的进口烟气焓值Hin为:
Figure GDA00031941561700000216
其中,
Figure GDA00031941561700000217
MAr
Figure GDA00031941561700000218
一一对应为N2、O2、CO2、H2O、Ar和SO2的摩尔质量;
步骤3、按照与步骤2相同的方式确定出口烟气温度Tout,计算获得出口烟气焓值Hout
步骤4、计算获得进口烟气质量流量Qin
按照网格法分别测量获得烟道中烟气进口网格点烟气动压ΔPi和网格点烟气静压Pi;则:通过计算平均值获得进口烟气静压Pin为:
Figure GDA00031941561700000219
其中,n为烟道中烟气进口网格点的数量;
进口烟气密度ρin为:
Figure GDA0003194156170000031
其中,L为理想气体的摩尔体积;
烟气进口网格点烟气流速Vi为:
Figure GDA0003194156170000032
通过计算平均值获得进口烟气流速Vin为:
Figure GDA0003194156170000033
计算获得进口烟气质量流量Qin为:Qin=ρinVinA
其中,A为余热锅炉烟道进口的截面积;
步骤5、计算获得余热锅炉热效率η
余热锅炉的输入热量Ein为:Ein=QinHin
针对密封状态良好,且无泄漏的无补燃余热锅炉,烟道中出口烟气质量流量Qout等于进口烟气质量流量Qin,则有:
余热锅炉的排烟热量Eout为:Eout=QinHout
余热锅炉的排烟热损失q2为:
Figure GDA0003194156170000034
若余热锅炉的输入热量Ein满足式(2):
1.055×108kJ/h≤Ein≤1.583×109kJ/h (2)
则余热锅炉的散热损失q5为:
Figure GDA0003194156170000035
若余热锅炉的输入热量Ein不满足式(2),按如下方式确定余热锅炉的散热损失q5
若0<De≤850t/h,则余热锅炉的散热损失q5由式(3a)计算获得:
Figure GDA0003194156170000036
若850<De<2000t/h,则余热锅炉的散热损失q5由式(3b)计算获得:
Figure GDA0003194156170000041
其中,D为余热锅炉给水流量的现场测量值,De为余热锅炉给水流量的额定值;
则计算获得余热锅炉热效率η为:
η=100-q2-q5-qoth
其中,qoth为不可计量热损失,qoth值包括连续排污热损失和定期排污热损失;
步骤6、计算获得余热锅炉出力,包括高压蒸发量Dhs、中压蒸发量Dms和低压蒸发量Dls
高压蒸发量Dhs为:Dhs=xDh
中压蒸发量Dms为:Dms=xDm
低压蒸发量Dls为:Dls=xDl
Dh为高压蒸发量设计值,Dm为中压蒸发量设计值,Dl为低压蒸发量设计值;
x是依据能量守恒原理由式(4)计算获得的蒸发量比例常数:
Figure GDA0003194156170000042
Qzr为再热减温水流量,Hzw为再热减温水焓值;
Hrz为热再热蒸汽焓值,Hlz为冷再热蒸汽焓值;
Hh为高压蒸汽焓值,Hm为中压蒸汽焓值,Hl为低压蒸汽焓值,Hw为总给水焓值;
若余热锅炉有再热减温水系统,则:再热减温水流量Qzr、再热减温水的压力Pzw和再热减温水的温度Tzw均为试验期间采集数据的算术平均值;再热减温水焓值Hzw根据再热减温水的压力Pzw和再热减温水的温度Tzw查表或者通过IAPWS-iF97公式计算获得;
若余热锅炉无再热减温水系统,则:再热减温水流量Qzr和再热减温水焓值Hzw均为零;
若余热锅炉有再热系统,则:热再热蒸汽的压力Prz和热再热蒸汽的温度Trz,以及冷再热蒸汽的压力Plz和冷再热蒸汽的温度Tlz均为试验期间采集数据的算术平均值;热再热蒸汽焓值Hrz和冷再热蒸汽焓值Hlz分别根据热再热蒸汽的压力Prz和热再热蒸汽的温度Trz,冷再热蒸汽的压力Plz和冷再热蒸汽的温度Tlz查表或者通过IAPWS-iF97公式计算获得;
若余热锅炉无再热系统,则:热再热蒸汽焓值Hrz和冷再热蒸汽焓值Hlz均为零;再热减温水流量Qzr和再热减温水焓值Hzw也均为零;
若余热锅炉为三压系统,则:高压蒸汽的压力Ph和高压蒸汽的温度Th、中压蒸汽的压力Pm和中压蒸汽的温度Tm、低压蒸汽的压力Pl和低压蒸汽的温度Tl,以及总给水的压力Pw和总给水的温度Tw均取为试验期间采集数据的算术平均值;高压蒸汽焓值Hh、中压蒸汽焓值Hm、低压蒸汽焓值Hl,以及总给水的焓值Hw分别根据Ph、Th、Pm、Tm、Pl、Tl、Pw以及Tw的值查表或者通过IAPWS-iF97公式计算获得;
若余热锅炉为双压系统,则:中压蒸发量Dms、中压蒸汽的压力Pm和中压蒸汽的温度Tm,以及中压蒸汽焓值Hm均为零;
若余热锅炉为单压系统,则:中压蒸发量Dms、中压蒸汽的压力Pm和中压蒸汽的温度Tm、中压蒸汽焓值Hm、低压蒸发量Dls、低压蒸汽的压力Pl和低压蒸汽的温度Tl,以及低压蒸汽焓值Hl均为零;
步骤7、计算获得高压汽水侧压降ΔPh、中压汽水侧压降ΔPm、低压汽水侧压降ΔPl,以及烟气侧压降ΔPg
高压汽水侧压降ΔPh为:ΔPh=Phw-Ph
Phw为高压给水泵出口压力,取为试验期间DCS采集数据的算术平均值;
中压汽水侧压降ΔPm为:ΔPm=Pmw-Pm
Pmw为中压给水泵出口压力,取为试验期间DCS采集数据的算术平均值;
低压汽水侧压降ΔPl为:ΔPl=Plw-Pl
Plw为低压给水泵出口压力,取为试验期间DCS采集数据的算术平均值;
烟气侧压降ΔPg为:ΔPg=Pin-Pout
Pout为余热锅炉出口烟气静压,按照网格法测量烟道中烟气出口各网格点的烟气静压,按照与式(1)相应的方式计算获得;若现场不具备直接测量的条件,则余热锅炉烟道中进口烟气静压Pin和出口烟气静压Pout均取为试验期间DCS采集数据的算术平均值。
本发明余热锅炉性能试验方法的特点也在于在所述步骤1中,采用红外烟气分析仪现场直接测量获得进口烟气组分和摩尔分数;或按如下方式计算获得进口烟气各组分的摩尔分数
Figure GDA0003194156170000051
XAr
Figure GDA0003194156170000052
的值:
步骤1.1、计算获得试验环境中空气各组分的摩尔分数,所述空气各组分包括H2O、N2、O2、Ar和CO2,一一对应的摩尔分数为
Figure GDA0003194156170000053
XAr,a
Figure GDA0003194156170000054
计算获得空气中水蒸汽分压
Figure GDA0003194156170000055
为:
Figure GDA0003194156170000056
HR为试验环境中的空气相对湿度,t0为试验环境中的空气温度;
计算获得空气中水蒸汽组分摩尔分数
Figure GDA0003194156170000061
为:
Figure GDA0003194156170000062
Pact为试验环境中的大气压力;空气中其它组分的摩尔分数
Figure GDA0003194156170000063
XAr,a
Figure GDA0003194156170000064
根据干空气各组分摩尔分数与空气中水蒸汽组分摩尔分数
Figure GDA0003194156170000065
一一对应计算获得;
步骤1.2、计算获得进入燃气轮机的空气体积流量Qa
单位体积的燃气完全燃烧所需空气中O2体积a为:
Figure GDA0003194156170000066
单位体积的燃气完全燃烧后体积的变化量b为:
Figure GDA0003194156170000067
其中,
Figure GDA0003194156170000068
Xr,CO
Figure GDA0003194156170000069
分别为进入燃气轮机的燃气组分H2、CO、CmHn、H2S和O2的摩尔分数,由在线燃气组分分析仪或取燃气样品分析获得;m和n为燃气组分CmHn的下角标系数;则,进入燃气轮机的空气体积流量Qa为:
Figure GDA00031941561700000610
其中,Qr为进入燃气轮机的燃料体积流量,由试验期间DCS采集数据求算术平均值获得;
Figure GDA00031941561700000611
为进口烟气中O2的摩尔分数,无补燃余热锅炉进口烟气中O2的摩尔分数与出口相同,由CEMS在线系统获得;
步骤1.3、计算获得进口烟气各组分的摩尔分数
Figure GDA00031941561700000612
XAr
Figure GDA00031941561700000613
的值;进口烟气中N2的摩尔分数
Figure GDA00031941561700000614
为:
Figure GDA00031941561700000615
进口烟气中CO2的摩尔分数
Figure GDA00031941561700000616
为:
Figure GDA00031941561700000617
进口烟气中H2O的摩尔分数
Figure GDA00031941561700000618
为:
Figure GDA0003194156170000071
进口烟气中Ar的摩尔分数XAr为:
Figure GDA0003194156170000072
进口烟气中SO2摩尔分数
Figure GDA0003194156170000073
为:
Figure GDA0003194156170000074
其中,
Figure GDA0003194156170000075
分别为进入燃气轮机的燃气组分N2、CO2和H2O的摩尔分数,由在线燃气组分分析仪或取燃气样品分析获得;无补燃余热锅炉出口烟气中各组分的摩尔分数与进口相同。
本发明余热锅炉性能试验方法的特点也在于:
当现场不具备条件直接测量进口烟气质量流量时,将权利要求1的步骤4中计算获得进口烟气质量流量Qin的方法替换为由式(10)计算获得进口烟气质量流量Qin
Figure GDA0003194156170000076
本发明余热锅炉性能试验方法的特点也在于:
步骤2中所述的进口烟气温度Tin按如下方式获得:采用热电偶按照网格法现场测量获得进口烟气网格点温度Ti,i表示进口烟气网格点,i=1,2,…,n,n为烟道中烟气进口网格点的数量;通过计算平均值获得进口烟气温度Tin为:
Figure GDA0003194156170000077
若现场不具备直接测量的条件,则进口烟气温度Tin取自DCS数据。
本发明余热锅炉性能试验方法的特点也在于:针对补燃余热锅炉,因补燃余热锅炉的出口烟气质量流量Qout大于进口烟气质量流量Qin,余热锅炉的出口烟气组分与进口烟气组分存在有不同,并有:
补燃燃气质量流量Qbr,m为:Qbr,m=Qout-Qin
进入余热锅炉的热量Ein为:Ein=QinHin+Qbr,mHbr
其中,Hbr为补燃燃气的低位热值,由在线燃气组分分析仪或取燃气样品分析获得;
出口烟气热量Eout为:Eout=QoutHout
若现场具备直接测量烟气组分、温度和流量的条件,则按步骤1到步骤4的相应方法通过测量和计算同样获得余热锅炉出口烟气组分的摩尔分数
Figure GDA0003194156170000081
Xc,Ar
Figure GDA0003194156170000082
出口烟气焓Hout,以及出口烟气质量流量Qout
在现场不具备直接测量烟气组分、温度和流量的条件时,在所述步骤1到步骤4中,按如下方式计算获得进口烟气各组分的摩尔分数
Figure GDA0003194156170000083
XAr
Figure GDA0003194156170000084
进口烟气质量流量Qin以及进口烟气焓值Hin;出口烟气组分的摩尔分数
Figure GDA0003194156170000085
Xc,Ar
Figure GDA0003194156170000086
出口烟气质量流量Qout以及出口烟气焓值Hout
步骤a、计算获得进入燃气轮机的空气体积流量Qa
单位体积的补燃燃气完全燃烧所需空气中O2的体积a2为:
Figure GDA0003194156170000087
单位体积的补燃燃气完全燃烧后体积的变化量b2为:
Figure GDA0003194156170000088
其中,
Figure GDA0003194156170000089
Xbr,CO
Figure GDA00031941561700000810
Figure GDA00031941561700000811
一一对应为进入余热锅炉的补燃燃气组分H2、CO、CmHn、H2S和O2的摩尔分数,是由在线燃气组分分析仪或取燃气样品分析获得;m和n为燃气组分CmHn的下角标系数;进入燃气轮机的空气体积流量Qa为:
Figure GDA00031941561700000812
其中,Qbr为补燃燃气的体积流量,由试验期间DCS采集数据求算术平均值获得;
Figure GDA00031941561700000813
为出口烟气中O2的摩尔分数,由CEMS在线系统获得;
步骤b、计算获得余热锅炉进口烟气各组分的摩尔分数
Figure GDA00031941561700000814
XAr
Figure GDA00031941561700000815
以及进口烟气质量流量Qin的值;余热锅炉进口烟气中O2的摩尔分数
Figure GDA00031941561700000816
为:
Figure GDA00031941561700000817
余热锅炉进口烟气中其它各组分的摩尔分数
Figure GDA00031941561700000818
XAr
Figure GDA00031941561700000819
以及进口烟气质量流量Qin按照式(5)-(10)分别计算获得;
步骤c、计算获得余热锅炉出口烟气各组分摩尔分数
Figure GDA0003194156170000091
Xc,Ar
Figure GDA0003194156170000092
以及出口烟气质量流量Qout的值,则有:
出口烟气中N2的摩尔分数
Figure GDA0003194156170000093
为:
Figure GDA0003194156170000094
出口烟气中CO2的摩尔分数
Figure GDA0003194156170000095
为:
Figure GDA0003194156170000096
出口烟气中H2O的摩尔分数
Figure GDA0003194156170000097
为:
Figure GDA0003194156170000098
出口烟气中Ar的摩尔分数Xc,Ar为:
Figure GDA0003194156170000099
出口烟气中SO2的摩尔分数
Figure GDA00031941561700000910
为:
Figure GDA00031941561700000911
余热锅炉出口烟气质量流量Qout为:
Figure GDA00031941561700000912
步骤d、余热锅炉的进口烟气焓值Hin和出口烟气焓值Hout按权利要求1中的步骤2的方式计算获得,余热锅炉的进口烟气温度Tin和出口烟气温度Tout取自DCS数据。
与已有技术相比,本发明有益效果体现在:
1、本发明根据现场实际情况给出的余热锅炉性能试验方法,在现场实际应用中更具有可操作性;
2、本发明通过实际测量条件与风险评估,选择相应的试验过程和计算方法,很大程度上保障了现场试验人员人身和设备的安全;
3、本发明依据目前测量水平和技术,从实际测量与理论计算相结合的方法根据实际可准确测量的数据计算获得余热锅炉的性能指标,从而保证了试验结果的准确性与可信性。
附图说明
图1本发明方法操作流程示意图;
具体实施方式
本实施例中余热锅炉性能试验方法按如下步骤进行:
步骤1、确定余热锅炉烟道中进口烟气组分和摩尔分数,进口烟气组分包括N2、O2、CO2、H2O、Ar和SO2,一一对应的摩尔分数为
Figure GDA0003194156170000101
XAr
Figure GDA0003194156170000102
无补燃余热锅炉烟道中出口烟气组分和摩尔分数与进口相同。
步骤2、计算获得余热锅炉进口烟气焓值Hin
确定进口烟气温度Tin,并根据进口烟气温度Tin计算获得进口烟气中各组分的焓值,包括:
进口烟气中N2的焓值
Figure GDA0003194156170000103
为:
Figure GDA0003194156170000104
进口烟气中O2的焓值
Figure GDA0003194156170000105
为:
Figure GDA0003194156170000106
进口烟气中CO2的焓值
Figure GDA0003194156170000107
为:
Figure GDA0003194156170000108
进口烟气中H2O的焓值
Figure GDA0003194156170000109
为:
Figure GDA00031941561700001010
进口烟气中Ar的焓值HAr为:
Figure GDA00031941561700001011
进口烟气中SO2的焓值
Figure GDA00031941561700001012
为:
Figure GDA00031941561700001013
其中,
Figure GDA00031941561700001014
aj,Ar
Figure GDA00031941561700001015
为查表获得的如表1-1和表1-2所示的经验系数,j=1,2,…,6;
计算获得余热锅炉的进口烟气焓值Hin为:
Figure GDA00031941561700001016
其中,
Figure GDA00031941561700001017
MAr
Figure GDA00031941561700001018
一一对应为查表获得的N2、O2、CO2、H2O、Ar和SO2的摩尔质量,如表2所示。
步骤3、按照与步骤2相同的方式确定出口烟气温度Tout,计算获得出口烟气焓值Hout
步骤4、计算获得进口烟气质量流量Qin
采用毕托管和微压计按照网格法分别测量获得烟道中烟气进口网格点烟气动压ΔPi和网格点烟气静压Pi;则:
通过计算平均值获得进口烟气静压Pin为:
Figure GDA0003194156170000111
其中,n为烟道中烟气进口网格点的数量;
进口烟气密度ρin为:
Figure GDA0003194156170000112
其中,L为查表获得的理想气体的摩尔体积;
烟气进口网格点烟气流速Vi为:
Figure GDA0003194156170000113
通过计算平均值获得进口烟气流速Vin为:
Figure GDA0003194156170000114
计算获得进口烟气质量流量Qin为:Qin=ρinVinA
其中,A为余热锅炉烟道进口的截面积。
步骤5、计算获得余热锅炉热效率η
余热锅炉的输入热量Ein为:Ein=QinHin
针对密封状态良好,且无泄漏的无补燃余热锅炉,烟道中出口烟气质量流量Qout等于进口烟气质量流量Qin,则有:
余热锅炉的排烟热量Eout为:Eout=QinHout
余热锅炉的排烟热损失q2为:
Figure GDA0003194156170000115
若余热锅炉的输入热量Ein满足式(2):
1.055×108kJ/h≤Ein≤1.583×109kJ/h (2)
则余热锅炉的散热损失q5为:
Figure GDA0003194156170000121
若余热锅炉的输入热量Ein不满足式(2),按如下方式确定余热锅炉的散热损失q5
若0<De≤850t/h,则余热锅炉的散热损失q5由式(3a)计算获得:
Figure GDA0003194156170000122
若850<De<2000t/h,则余热锅炉的散热损失q5由式(3b)计算获得:
Figure GDA0003194156170000123
其中,D为余热锅炉给水流量的现场测量值,De为锅炉设计说明书中记载的余热锅炉给水流量的额定值;则计算获得余热锅炉热效率η为:
η=100-q2-q5-qoth
其中,qoth为不可计量热损失,包括连续排污热损失、定期排污热损失,及其它不可测量热损失。
步骤6、计算获得余热锅炉出力,包括高压蒸发量Dhs、中压蒸发量Dms和低压蒸发量Dls高压蒸发量Dhs为:Dhs=xDh,中压蒸发量Dms为:Dms=xDm,低压蒸发量Dls为:Dls=xDl;Dh为高压蒸发量设计值,Dm为中压蒸发量设计值,Dl为低压蒸发量设计值,Dh、Dm和Dl均由锅炉设计说明书获得,x是依据能量守恒原理由式(4)计算获得的蒸发量比例常数:
Figure GDA0003194156170000124
Qzr为再热减温水流量,Hzw为再热减温水焓值;
Hrz为热再热蒸汽焓值,Hlz为冷再热蒸汽焓值;
Hh为高压蒸汽焓值,Hm为中压蒸汽焓值,Hl为低压蒸汽焓值,Hw为总给水焓值。
若余热锅炉有再热减温水系统,则:再热减温水流量Qzr、再热减温水的压力Pzw和再热减温水的温度Tzw均为试验期间采集数据的算术平均值;再热减温水焓值Hzw根据再热减温水的压力Pzw和再热减温水的温度Tzw查表或者通过IAPWS-iF97公式计算获得。
若余热锅炉无再热减温水系统,则:再热减温水流量Qzr和再热减温水焓值Hzw均为零。
若余热锅炉有再热系统,则:热再热蒸汽的压力Prz和热再热蒸汽的温度Trz,以及冷再热蒸汽的压力Plz和冷再热蒸汽的温度Tlz均为试验期间采集数据的算术平均值;热再热蒸汽焓值Hrz和冷再热蒸汽焓值Hlz分别根据热再热蒸汽的压力Prz和热再热蒸汽的温度Trz,冷再热蒸汽的压力Plz和冷再热蒸汽的温度Tlz查表或者通过IAPWS-iF97公式计算获得。
若余热锅炉无再热系统,则:热再热蒸汽焓值Hrz和冷再热蒸汽焓值Hlz均为零;再热减温水流量Qzr和再热减温水焓值Hzw也均为零。
若余热锅炉为三压系统,则:高压蒸汽的压力Ph和高压蒸汽的温度Th、中压蒸汽的压力Pm和中压蒸汽的温度Tm、低压蒸汽的压力Pl和低压蒸汽的温度Tl,以及总给水的压力Pw和总给水的温度Tw均取为试验期间采集数据的算术平均值;高压蒸汽焓值Hh、中压蒸汽焓值Hm、低压蒸汽焓值Hl,以及总给水的焓值Hw分别根据Ph、Th、Pm、Tm、Pl、Tl、Pw以及Tw的值查表或者通过IAPWS-iF97公式计算获得。
若余热锅炉为双压系统,则:中压蒸发量Dms、中压蒸汽的压力Pm和中压蒸汽的温度Tm,以及中压蒸汽焓值Hm均为零。
若余热锅炉为单压系统,则:中压蒸发量Dms、中压蒸汽的压力Pm和中压蒸汽的温度Tm、中压蒸汽焓值Hm、低压蒸发量Dls、低压蒸汽的压力Pl和低压蒸汽的温度Tl,以及低压蒸汽焓值Hl均为零。
步骤7、计算获得高压汽水侧压降ΔPh、中压汽水侧压降ΔPm、低压汽水侧压降ΔPl,以及烟气侧压降ΔPg
高压汽水侧压降ΔPh为:ΔPh=Phw-Ph
Phw为高压给水泵出口压力,取为试验期间DCS采集数据的算术平均值;
中压汽水侧压降ΔPm为:ΔPm=Pmw-Pm
Pmw为中压给水泵出口压力,取为试验期间DCS采集数据的算术平均值;
低压汽水侧压降ΔPl为:ΔPl=Plw-Pl
Plw为低压给水泵出口压力,取为试验期间DCS采集数据的算术平均值;
烟气侧压降ΔPg为:ΔPg=Pin-Pout
Pout为余热锅炉出口烟气静压,按照网格法测量烟道中烟气出口各网格点的烟气静压,按照与式(1)相应的方式、根据烟道中烟气出口网格点的烟气静压求算术平均值获得。若现场不具备直接测量的条件,则余热锅炉烟道中进口烟气静压Pin和出口烟气静压Pout均取为试验期间DCS采集数据的算术平均值。
具体实施中,在步骤1中是采用红外烟气分析仪现场直接测量获得进口烟气组分和摩尔分数;或按如下方式计算获得进口烟气各组分的摩尔分数
Figure GDA0003194156170000141
XAr
Figure GDA0003194156170000142
的值:
步骤1.1、计算获得试验环境中空气各组分的摩尔分数,空气各组分包括H2O、N2、O2、Ar和CO2,一一对应的摩尔分数为
Figure GDA0003194156170000143
XAr,a
Figure GDA0003194156170000144
计算获得空气中水蒸汽分压
Figure GDA0003194156170000145
为:
Figure GDA0003194156170000146
HR为试验环境中的空气相对湿度,单位%;
t0为试验环境中的空气温度,单位℃,HR和t0均为现场直接测量获得;
计算获得空气中水蒸汽组分摩尔分数
Figure GDA0003194156170000147
为:
Figure GDA0003194156170000148
Pact为试验环境中的大气压力,单位Pa,由现场直接测量获得;空气中其它组分的摩尔分数
Figure GDA0003194156170000149
XAr,a
Figure GDA00031941561700001410
根据干空气各组分摩尔分数与空气中水蒸汽组分摩尔分数
Figure GDA00031941561700001411
一一对应计算获得,干空气各组分摩尔分数为查表获得;
步骤1.2、计算获得进入燃气轮机的空气体积流量Qa
单位体积的燃气完全燃烧所需空气中O2体积a为:
Figure GDA00031941561700001412
单位体积的燃气完全燃烧后体积的变化量b为:
Figure GDA00031941561700001413
其中,
Figure GDA00031941561700001414
Xr,CO
Figure GDA00031941561700001415
分别为进入燃气轮机的燃气组分H2、CO、CmHn、H2S和O2的摩尔分数,由在线燃气组分分析仪或取燃气样品分析获得;m和n为燃气组分CmHn的下角标系数;则,进入燃气轮机的空气体积流量Qa为:
Figure GDA00031941561700001416
其中,Qr为进入燃气轮机的燃料体积流量,单位Nm3/h,由试验期间DCS采集数据求算术平均值获得;
Figure GDA00031941561700001417
为进口烟气中O2的摩尔分数,无补燃余热锅炉进口烟气中O2的摩尔分数与出口相同,由CEMS在线系统获得。
步骤1.3、计算获得进口烟气各组分的摩尔分数
Figure GDA0003194156170000151
XAr
Figure GDA0003194156170000152
的值;进口烟气中N2的摩尔分数
Figure GDA0003194156170000153
为:
Figure GDA0003194156170000154
进口烟气中CO2的摩尔分数
Figure GDA0003194156170000155
为:
Figure GDA0003194156170000156
进口烟气中H2O的摩尔分数
Figure GDA0003194156170000157
为:
Figure GDA0003194156170000158
进口烟气中Ar的摩尔分数XAr为:
Figure GDA0003194156170000159
进口烟气中SO2摩尔分数
Figure GDA00031941561700001510
为:
Figure GDA00031941561700001511
其中,
Figure GDA00031941561700001512
分别为进入燃气轮机的燃气组分N2、CO2和H2O的摩尔分数,由在线燃气组分分析仪或取燃气样品分析获得;无补燃余热锅炉出口烟气中各组分的摩尔分数与进口相同。
当现场不具备条件直接测量进口烟气质量流量时,将步骤4中计算获得进口烟气质量流量Qin的方法替换为由式(10)计算获得进口烟气质量流量Qin
Figure GDA00031941561700001513
具体实施中,步骤2中的进口烟气温度Tin按如下方式获得:采用热电偶按照网格法现场测量获得进口烟气网格点温度Ti,以i表示进口烟气网格点,i=1,2,…,n,n为烟道中烟气进口网格点的数量,通过计算平均值获得进口烟气温度Tin为:
Figure GDA0003194156170000161
若现场不具备直接测量的条件,则进口烟气温度Tin取自DCS数据。
具体实施中,针对补燃余热锅炉,因补燃余热锅炉的出口烟气质量流量Qout大于进口烟气质量流量Qin,余热锅炉的出口烟气组分与进口烟气组分存在有不同,并有:
补燃燃气质量流量Qbr,m为:Qbr,m=Qout-Qin
进入余热锅炉的热量Ein为:Ein=QinHin+Qbr,mHbr
其中,Hbr为补燃燃气的低位热值,由在线燃气组分分析仪或取燃气样品分析获得;
出口烟气热量Eout为:Eout=QoutHout
若现场具备直接测量烟气组分、温度和流量的条件,则按步骤1到步骤4的相应方法通过测量和计算同样获得余热锅炉出口烟气组分的摩尔分数
Figure GDA0003194156170000162
Xc,Ar
Figure GDA0003194156170000163
出口烟气焓Hout,以及出口烟气质量流量Qout
在现场不具备直接测量烟气组分、温度和流量的条件,在步骤1至步骤4中,按如下方式计算获得进口烟气各组分的摩尔分数
Figure GDA0003194156170000164
XAr
Figure GDA0003194156170000165
进口烟气质量流量Qin以及进口烟气焓值Hin;出口烟气组分的摩尔分数
Figure GDA0003194156170000166
Xc,Ar
Figure GDA0003194156170000167
出口烟气质量流量Qout以及出口烟气焓值Hout
步骤a:计算获得进入燃气轮机的空气体积流量Qa
单位体积的补燃燃气完全燃烧所需空气中O2的体积a2为:
Figure GDA0003194156170000168
单位体积的补燃燃气完全燃烧后体积的变化量b2为:
Figure GDA0003194156170000169
其中,
Figure GDA00031941561700001610
Xbr,CO
Figure GDA00031941561700001611
Figure GDA00031941561700001612
一一对应为进入余热锅炉的补燃燃气组分H2、CO、CmHn、H2S和O2的摩尔分数,由在线燃气组分分析仪或取燃气样品分析获得;m和n为燃气组分CmHn的下角标系数;进入燃气轮机的空气体积流量Qa为:
Figure GDA00031941561700001613
其中,Qbr为补燃燃气的体积流量,单位Nm3/h,由试验期间DCS采集数据求算术平均值获得;
Figure GDA0003194156170000171
为出口烟气中O2的摩尔分数,由CEMS在线系统获得。
步骤b、计算获得余热锅炉进口烟气各组分的摩尔分数
Figure GDA0003194156170000172
XAr
Figure GDA0003194156170000173
以及进口烟气质量流量Qin的值;余热锅炉进口烟气中O2的摩尔分数
Figure GDA0003194156170000174
为:
Figure GDA0003194156170000175
余热锅炉进口烟气中其它各组分的摩尔分数
Figure GDA0003194156170000176
XAr
Figure GDA0003194156170000177
以及进口烟气质量流量Qin按照式(5)-(10)分别计算获得。
步骤c、计算获得余热锅炉出口烟气各组分摩尔分数
Figure GDA0003194156170000178
Xc,Ar
Figure GDA0003194156170000179
以及出口烟气质量流量Qout的值,则有:
出口烟气中N2的摩尔分数
Figure GDA00031941561700001710
为:
Figure GDA00031941561700001711
出口烟气中CO2的摩尔分数
Figure GDA00031941561700001712
为:
Figure GDA00031941561700001713
出口烟气中H2O的摩尔分数
Figure GDA00031941561700001714
为:
Figure GDA00031941561700001715
出口烟气中Ar的摩尔分数Xc,Ar为:
Figure GDA00031941561700001716
出口烟气中SO2的摩尔分数
Figure GDA00031941561700001717
为:
Figure GDA00031941561700001718
余热锅炉出口烟气质量流量Qout为:
Figure GDA00031941561700001719
步骤d、余热锅炉的进口烟气焓值Hin和出口烟气焓值Hout按步骤2的相应方式计算获得,余热锅炉的进口烟气温度Tin和出口烟气温度Tout取自DCS数据。
如图1所示的操作流程,首先收集余热锅炉类型的资料和现场实际运行状态的信息,判断试验对象余热锅炉是否是补燃锅炉,实际现场是否进行补燃运行;评估测试风险确定实际现场是否具备直接测量余热锅炉进、出口烟气组分的摩尔分数、流量和温度的试验条件;最后根据锅炉类型、运行状态与试验测试条件,按照相应的步骤获得试验数据,并对试验数据进行处理,计算获得余热锅炉热效率、锅炉出力和锅炉汽水侧、烟气侧压降的结果。
表1-1:查表获得的根据温度求取N2、O2和CO2气体焓值的经验系数
气体 N<sub>2</sub> O<sub>2</sub> CO<sub>2</sub>
系数a<sub>1</sub> -6.5604875×10<sup>3</sup> 8.9009217×10<sup>3</sup> -9.3398089×10<sup>3</sup>
系数a<sub>2</sub> -1.1333376×10<sup>2</sup> 1.259436×10<sup>2</sup> -1.1834503×10<sup>2</sup>
系数a<sub>3</sub> 2.20796×10<sup>2</sup> -8.72569×10<sup>2</sup> 3.94483×10<sup>2</sup>
系数a<sub>4</sub> 1.8053856 0.29076177 1.0016303
系数a<sub>5</sub> -1.2660079×10<sup>-3</sup> 5.5785819×10<sup>-4</sup> 2.3651695×10<sup>-4</sup>
系数a<sub>6</sub> 1.3698985×10<sup>-6</sup> -5.9211059×10<sup>-8</sup> -1.3396751×10<sup>-8</sup>
系数a<sub>7</sub> -7.142453×10<sup>-10</sup> -1.3143737×10<sup>-10</sup> -3.6320845×10<sup>-11</sup>
系数a<sub>8</sub> 1.4957213×10<sup>-13</sup> 5.3996454×10<sup>-14</sup> 1.0767527×10<sup>-14</sup>
表1-2:查表获得的根据温度求取H2O、Ar和SO2气体焓值的经验系数
气体 H<sub>2</sub>O Ar SO<sub>2</sub>
系数a<sub>1</sub> 1.8220844×10<sup>4</sup> 0 6.8926679×10<sup>3</sup>
系数a<sub>2</sub> 2.6564163×10<sup>2</sup> 0 1.1797846×10<sup>2</sup>
系数a<sub>3</sub> -1.80771×10<sup>3</sup> -1.50225×10<sup>2</sup> -7.00143×10<sup>2</sup>
系数a<sub>4</sub> 0.4300414 0.52033331 -0.30588875
系数a<sub>5</sub> 1.6667329×10<sup>-3</sup> 0 1.4305209×10<sup>-3</sup>
系数a<sub>6</sub> -1.1295923×10<sup>-6</sup> 0 -1.0862046×10<sup>-6</sup>
系数a<sub>7</sub> 5.7171965×10<sup>-10</sup> 0 4.6926969×10<sup>-10</sup>
系数a<sub>8</sub> -1.2340574×10<sup>-13</sup> 0 -8.7450866×10<sup>-14</sup>
表2查表获得的各气体的摩尔质量
气体组分 N<sub>2</sub> O<sub>2</sub> CO<sub>2</sub> H<sub>2</sub>O Ar SO<sub>2</sub>
摩尔质量M 28.01348 31.9988 44.0098 18.01528 39.948 64.0648
应用实例:
某电厂某台燃气蒸汽联合循环余热锅炉为双压带自身除氧器、卧式、无补燃型、无再热系统的自然循环燃机余热锅炉。经现场勘查与交流,发现现场测量环境恶劣、测量风险较高,而且没有试验测点,因此按如下方式操作:
第一步:采用温湿度仪与大气压力表测量空气相对湿度HR为80%、环境温度t0为15.6℃,大气压Pact为101325Pa;
第二步:现场取燃气样品化验分析获得进入燃气轮机的燃气组分以及各组分的摩尔分数(体积分数),实际现场采用高炉煤气和焦炉煤气作为燃机的燃料,分别取样化验各燃气的组分,根据混合比例,计算出进入燃气轮机的燃气各组的平均分摩尔分数(体积分数)如表3:
表3
Figure GDA0003194156170000191
第三步:通过电厂DCS分别获得进入燃气轮机的燃料体积流量为267061Nm3/h(高炉煤气+焦炉煤气),余热锅炉进、出口烟气温度分别为575.1℃和136.3℃,高压蒸汽压力7.5MPa(表压)、温度505℃,低压蒸汽压力0.647MPa(表压)、温度278.6℃,除氧蒸发器进水压力1.173MPa,进水温度58℃,高压给水泵出口压力10MPa(表压),低压给水泵出口压力1.26MPa(表压),余热锅炉烟气进口压力3.14kPa(表压),出口压力0.1kPa(表压)。查设计参数,高压蒸汽量与低压蒸汽量之比为6.35:1。
第四步、根据CEMS在线系统获得出口烟气中O2(湿基)的摩尔分数为0.118(体积百分数为11.8%)。
第五步、按照步骤1.1至步骤1.3,步骤2至步骤7的相应公式对以上试验数据进行处理与计算,获得锅炉热效率、锅炉出力和锅炉汽水侧、烟气侧压降。的结果如表4
表4
Figure GDA0003194156170000201
Figure GDA0003194156170000211
本实施例中计算获得的锅炉热效率、锅炉出力和锅炉汽水侧、烟气侧压降如表4。

Claims (5)

1.一种余热锅炉性能试验方法,其特征是按如下步骤进行:
步骤1、确定余热锅炉烟道中进口烟气组分和摩尔分数,所述进口烟气组分包括N2、O2、CO2、H2O、Ar和SO2,一一对应的摩尔分数为
Figure FDA0003194156160000011
XAr
Figure FDA0003194156160000012
无补燃余热锅炉烟道中出口烟气组分和摩尔分数与进口相同;
步骤2、计算获得余热锅炉进口烟气焓值Hin
确定进口烟气温度Tin,并根据进口烟气温度Tin计算获得进口烟气中各组分的焓值,包括:
进口烟气中N2的焓值
Figure FDA0003194156160000013
为:
Figure FDA0003194156160000014
进口烟气中O2的焓值
Figure FDA0003194156160000015
为:
Figure FDA0003194156160000016
进口烟气中CO2的焓值
Figure FDA0003194156160000017
为:
Figure FDA0003194156160000018
进口烟气中H2O的焓值
Figure FDA0003194156160000019
为:
Figure FDA00031941561600000110
进口烟气中Ar的焓值HAr为:
Figure FDA00031941561600000111
进口烟气中SO2的焓值
Figure FDA00031941561600000112
为:
Figure FDA00031941561600000113
其中,
Figure FDA00031941561600000114
aj,Ar
Figure FDA00031941561600000115
均为经验系数,j=1,2,…,6;
计算获得余热锅炉的进口烟气焓值Hin为:
Figure FDA00031941561600000116
其中,
Figure FDA00031941561600000117
MAr
Figure FDA00031941561600000118
一一对应为N2、O2、CO2、H2O、Ar和SO2的摩尔质量;
步骤3、按照与步骤2相同的方式确定出口烟气温度Tout,计算获得出口烟气焓值Hout
步骤4、计算获得进口烟气质量流量Qin
按照网格法分别测量获得烟道中烟气进口网格点烟气动压ΔPi和网格点烟气静压Pi;则:通过计算平均值获得进口烟气静压Pin为:
Figure FDA00031941561600000119
其中,n为烟道中烟气进口网格点的数量;
进口烟气密度ρin为:
Figure FDA0003194156160000021
其中,L为理想气体的摩尔体积;
烟气进口网格点烟气流速Vi为:
Figure FDA0003194156160000022
通过计算平均值获得进口烟气流速Vin为:
Figure FDA0003194156160000023
计算获得进口烟气质量流量Qin为:Qin=ρinVinA
其中,A为余热锅炉烟道进口的截面积;
步骤5、计算获得余热锅炉热效率η
余热锅炉的输入热量Ein为:Ein=QinHin
针对密封状态良好,且无泄漏的无补燃余热锅炉,烟道中出口烟气质量流量Qout等于进口烟气质量流量Qin,则有:
余热锅炉的排烟热量Eout为:Eout=QinHout
余热锅炉的排烟热损失q2为:
Figure FDA0003194156160000024
若余热锅炉的输入热量Ein满足式(2):
1.055×108kJ/h≤Ein≤1.583×109kJ/h (2)
则余热锅炉的散热损失q5为:
Figure FDA0003194156160000025
若余热锅炉的输入热量Ein不满足式(2),按如下方式确定余热锅炉的散热损失q5
若0<De≤850t/h,则余热锅炉的散热损失q5由式(3a)计算获得:
Figure FDA0003194156160000026
若850<De<2000t/h,则余热锅炉的散热损失q5由式(3b)计算获得:
Figure FDA0003194156160000031
其中,D为余热锅炉给水流量的现场测量值,De为余热锅炉给水流量的额定值;
则计算获得余热锅炉热效率η为:
η=100-q2-q5-qoth
其中,qoth为不可计量热损失,qoth值包括连续排污热损失和定期排污热损失;
步骤6、计算获得余热锅炉出力,包括高压蒸发量Dhs、中压蒸发量Dms和低压蒸发量Dls高压蒸发量Dhs为:Dhs=xDh
中压蒸发量Dms为:Dms=xDm
低压蒸发量Dls为:Dls=xDl
Dh为高压蒸发量设计值,Dm为中压蒸发量设计值,Dl为低压蒸发量设计值;
x是依据能量守恒原理由式(4)计算获得的蒸发量比例常数:
Figure FDA0003194156160000032
Qzr为再热减温水流量,Hzw为再热减温水焓值;
Hrz为热再热蒸汽焓值,Hlz为冷再热蒸汽焓值;
Hh为高压蒸汽焓值,Hm为中压蒸汽焓值,Hl为低压蒸汽焓值,Hw为总给水焓值;
若余热锅炉有再热减温水系统,则:再热减温水流量Qzr、再热减温水的压力Pzw和再热减温水的温度Tzw均为试验期间采集数据的算术平均值;再热减温水焓值Hzw根据再热减温水的压力Pzw和再热减温水的温度Tzw查表或者通过IAPWS-iF97公式计算获得;
若余热锅炉无再热减温水系统,则:再热减温水流量Qzr和再热减温水焓值Hzw均为零;
若余热锅炉有再热系统,则:热再热蒸汽的压力Prz和热再热蒸汽的温度Trz,以及冷再热蒸汽的压力Plz和冷再热蒸汽的温度Tlz均为试验期间采集数据的算术平均值;热再热蒸汽焓值Hrz和冷再热蒸汽焓值Hlz分别根据热再热蒸汽的压力Prz和热再热蒸汽的温度Trz,冷再热蒸汽的压力Plz和冷再热蒸汽的温度Tlz查表或者通过IAPWS-iF97公式计算获得;
若余热锅炉无再热系统,则:热再热蒸汽焓值Hrz和冷再热蒸汽焓值Hlz均为零;再热减温水流量Qzr和再热减温水焓值Hzw也均为零;
若余热锅炉为三压系统,则:高压蒸汽的压力Ph和高压蒸汽的温度Th、中压蒸汽的压力Pm和中压蒸汽的温度Tm、低压蒸汽的压力Pl和低压蒸汽的温度Tl,以及总给水的压力Pw和总给水的温度Tw均取为试验期间采集数据的算术平均值;高压蒸汽焓值Hh、中压蒸汽焓值Hm、低压蒸汽焓值Hl,以及总给水的焓值Hw分别根据Ph、Th、Pm、Tm、Pl、Tl、Pw以及Tw的值查表或者通过IAPWS-iF97公式计算获得;
若余热锅炉为双压系统,则:中压蒸发量Dms、中压蒸汽的压力Pm和中压蒸汽的温度Tm,以及中压蒸汽焓值Hm均为零;
若余热锅炉为单压系统,则:中压蒸发量Dms、中压蒸汽的压力Pm和中压蒸汽的温度Tm、中压蒸汽焓值Hm、低压蒸发量Dls、低压蒸汽的压力Pl和低压蒸汽的温度Tl,以及低压蒸汽焓值Hl均为零;
步骤7、计算获得高压汽水侧压降ΔPh、中压汽水侧压降ΔPm、低压汽水侧压降ΔPl,以及烟气侧压降ΔPg
高压汽水侧压降ΔPh为:ΔPh=Phw-Ph
Phw为高压给水泵出口压力,取为试验期间DCS采集数据的算术平均值;
中压汽水侧压降ΔPm为:ΔPm=Pmw-Pm
Pmw为中压给水泵出口压力,取为试验期间DCS采集数据的算术平均值;
低压汽水侧压降ΔPl为:ΔPl=Plw-Pl
Plw为低压给水泵出口压力,取为试验期间DCS采集数据的算术平均值;
烟气侧压降ΔPg为:ΔPg=Pin-Pout
Pout为余热锅炉出口烟气静压,按照网格法测量烟道中烟气出口各网格点的烟气静压,按照与式(1)相应的方式计算获得;若现场不具备直接测量的条件,则余热锅炉烟道中进口烟气静压Pin和出口烟气静压Pout均取为试验期间DCS采集数据的算术平均值。
2.根据权利要求1所述的余热锅炉性能试验方法,其特征是在所述步骤1中,采用红外烟气分析仪现场直接测量获得进口烟气组分和摩尔分数;或按如下方式计算获得进口烟气各组分的摩尔分数
Figure FDA0003194156160000041
XAr
Figure FDA0003194156160000042
的值:
步骤1.1、计算获得试验环境中空气各组分的摩尔分数,所述空气各组分包括H2O、N2、O2、Ar和CO2,一一对应的摩尔分数为
Figure FDA0003194156160000043
XAr,a
Figure FDA0003194156160000044
计算获得空气中水蒸汽分压
Figure FDA0003194156160000045
为:
Figure FDA0003194156160000046
HR为试验环境中的空气相对湿度,t0为试验环境中的空气温度;
计算获得空气中水蒸汽组分摩尔分数
Figure FDA0003194156160000051
为:
Figure FDA0003194156160000052
Pact为试验环境中的大气压力;空气中其它组分的摩尔分数
Figure FDA0003194156160000053
XAr,a
Figure FDA0003194156160000054
根据干空气各组分摩尔分数与空气中水蒸汽组分摩尔分数
Figure FDA0003194156160000055
一一对应计算获得;
步骤1.2、计算获得进入燃气轮机的空气体积流量Qa
单位体积的燃气完全燃烧所需空气中O2体积a为:
Figure FDA0003194156160000056
单位体积的燃气完全燃烧后体积的变化量b为:
Figure FDA0003194156160000057
其中,
Figure FDA0003194156160000058
Xr,CO
Figure FDA0003194156160000059
分别为进入燃气轮机的燃气组分H2、CO、CmHn、H2S和O2的摩尔分数,由在线燃气组分分析仪或取燃气样品分析获得;m和n为燃气组分CmHn的下角标系数;则,进入燃气轮机的空气体积流量Qa为:
Figure FDA00031941561600000510
其中,Qr为进入燃气轮机的燃料体积流量,由试验期间DCS采集数据求算术平均值获得;
Figure FDA00031941561600000511
为进口烟气中O2的摩尔分数,无补燃余热锅炉进口烟气中O2的摩尔分数与出口相同,由CEMS在线系统获得;
步骤1.3、计算获得进口烟气各组分的摩尔分数
Figure FDA00031941561600000512
XAr
Figure FDA00031941561600000513
的值;
进口烟气中N2的摩尔分数
Figure FDA00031941561600000514
为:
Figure FDA00031941561600000515
进口烟气中CO2的摩尔分数
Figure FDA00031941561600000516
为:
Figure FDA00031941561600000517
进口烟气中H2O的摩尔分数
Figure FDA00031941561600000518
为:
Figure FDA0003194156160000061
进口烟气中Ar的摩尔分数XAr为:
Figure FDA0003194156160000062
进口烟气中SO2摩尔分数
Figure FDA0003194156160000063
为:
Figure FDA0003194156160000064
其中,
Figure FDA0003194156160000065
分别为进入燃气轮机的燃气组分N2、CO2和H2O的摩尔分数,由在线燃气组分分析仪或取燃气样品分析获得;无补燃余热锅炉出口烟气中各组分的摩尔分数与进口相同。
3.根据权利要求2所述的余热锅炉性能试验方法,其特征是:
当现场不具备条件直接测量进口烟气质量流量时,将权利要求1的步骤4中计算获得进口烟气质量流量Qin的方法替换为由式(10)计算获得进口烟气质量流量Qin
Figure FDA0003194156160000066
4.根据权利要求1所述的余热锅炉性能试验方法,其特征是:
步骤2中所述的进口烟气温度Tin按如下方式获得:采用热电偶按照网格法现场测量获得进口烟气网格点温度Ti,i表示进口烟气网格点,i=1,2,…,n,n为烟道中烟气进口网格点的数量;通过计算平均值获得进口烟气温度Tin为:
Figure FDA0003194156160000067
若现场不具备直接测量的条件,则进口烟气温度Tin取自DCS数据。
5.根据权利要求1所述的余热锅炉性能试验方法,其特征是:针对补燃余热锅炉,因补燃余热锅炉的出口烟气质量流量Qout大于进口烟气质量流量Qin,余热锅炉的出口烟气组分与进口烟气组分存在有不同,并有:
补燃燃气质量流量Qbr,m为:Qbr,m=Qout-Qin
进入余热锅炉的热量Ein为:Ein=QinHin+Qbr,mHbr
其中,Hbr为补燃燃气的低位热值,由在线燃气组分分析仪或取燃气样品分析获得;
出口烟气热量Eout为:Eout=QoutHout
若现场具备直接测量烟气组分、温度和流量的条件,则按步骤1到步骤4的相应方法通过测量和计算同样获得余热锅炉出口烟气组分的摩尔分数
Figure FDA0003194156160000071
Xc,Ar
Figure FDA0003194156160000072
出口烟气焓Hout,以及出口烟气质量流量Qout
在现场不具备直接测量烟气组分、温度和流量的条件时,在所述步骤1到步骤4中,按如下方式计算获得进口烟气各组分的摩尔分数
Figure FDA0003194156160000073
XAr
Figure FDA0003194156160000074
进口烟气质量流量Qin以及进口烟气焓值Hin;出口烟气组分的摩尔分数
Figure FDA0003194156160000075
Xc,Ar
Figure FDA0003194156160000076
出口烟气质量流量Qout以及出口烟气焓值Hout
步骤a、计算获得进入燃气轮机的空气体积流量Qa
单位体积的补燃燃气完全燃烧所需空气中O2的体积a2为:
Figure FDA0003194156160000077
单位体积的补燃燃气完全燃烧后体积的变化量b2为:
Figure FDA0003194156160000078
其中,
Figure FDA0003194156160000079
Xbr,CO
Figure FDA00031941561600000710
Figure FDA00031941561600000711
一一对应为进入余热锅炉的补燃燃气组分H2、CO、CmHn、H2S和O2的摩尔分数,是由在线燃气组分分析仪或取燃气样品分析获得;m和n为燃气组分CmHn的下角标系数;进入燃气轮机的空气体积流量Qa为:
Figure FDA00031941561600000712
其中,Qbr为补燃燃气的体积流量,由试验期间DCS采集数据求算术平均值获得;
Figure FDA00031941561600000713
为出口烟气中O2的摩尔分数,由CEMS在线系统获得;
步骤b、计算获得余热锅炉进口烟气各组分的摩尔分数
Figure FDA00031941561600000714
XAr
Figure FDA00031941561600000715
以及进口烟气质量流量Qin的值;余热锅炉进口烟气中O2的摩尔分数
Figure FDA00031941561600000716
为:
Figure FDA00031941561600000717
余热锅炉进口烟气中其它各组分的摩尔分数
Figure FDA00031941561600000718
XAr
Figure FDA00031941561600000719
以及进口烟气质量流量Qin按照式(5)-(10)分别计算获得;
步骤c、计算获得余热锅炉出口烟气各组分摩尔分数
Figure FDA0003194156160000081
Xc,Ar
Figure FDA0003194156160000082
以及出口烟气质量流量Qout的值,则有:
出口烟气中N2的摩尔分数
Figure FDA0003194156160000083
为:
Figure FDA0003194156160000084
出口烟气中CO2的摩尔分数
Figure FDA0003194156160000085
为:
Figure FDA0003194156160000086
出口烟气中H2O的摩尔分数
Figure FDA0003194156160000087
为:
Figure FDA0003194156160000088
出口烟气中Ar的摩尔分数Xc,Ar为:
Figure FDA0003194156160000089
出口烟气中SO2的摩尔分数
Figure FDA00031941561600000810
为:
Figure FDA00031941561600000811
余热锅炉出口烟气质量流量Qout为:
Figure FDA00031941561600000812
步骤d、余热锅炉的进口烟气焓值Hin和出口烟气焓值Hout按权利要求1中的步骤2的方式计算获得,余热锅炉的进口烟气温度Tin和出口烟气温度Tout取自DCS数据。
CN202010520009.0A 2020-06-09 2020-06-09 一种余热锅炉性能试验方法 Active CN111521430B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010520009.0A CN111521430B (zh) 2020-06-09 2020-06-09 一种余热锅炉性能试验方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010520009.0A CN111521430B (zh) 2020-06-09 2020-06-09 一种余热锅炉性能试验方法

Publications (2)

Publication Number Publication Date
CN111521430A CN111521430A (zh) 2020-08-11
CN111521430B true CN111521430B (zh) 2021-09-28

Family

ID=71913064

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010520009.0A Active CN111521430B (zh) 2020-06-09 2020-06-09 一种余热锅炉性能试验方法

Country Status (1)

Country Link
CN (1) CN111521430B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113487190A (zh) * 2021-07-09 2021-10-08 西安热工研究院有限公司 一种整体煤气化联合循环发电机组全厂性能评估的试验方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101576248A (zh) * 2009-06-04 2009-11-11 重庆新仪自控系统工程有限公司 注汽锅炉最优化运行控制系统
CN103984868A (zh) * 2014-05-21 2014-08-13 华中科技大学 一种富氧燃烧锅炉热力性能获取方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7039555B2 (en) * 1998-03-24 2006-05-02 Exergetic Systems Llc Method for detecting heat exchanger tube failures and their location when using input/loss performance monitoring of a recovery boiler
AU2007217133B2 (en) * 2006-02-27 2013-05-30 Highview Enterprises Limited A method of storing energy and a cryogenic energy storage system
CN108021782A (zh) * 2016-11-15 2018-05-11 赵志渊 一种确定单压余热锅炉热效率及排烟参数的计算方法
CN108019732A (zh) * 2016-11-23 2018-05-11 林振娴 一种确定双压余热锅炉高低压蒸汽流量的计算方法
CN107796851B (zh) * 2017-11-21 2020-08-11 江苏海事职业技术学院 高炉煤气锅炉入炉煤气热值与锅炉热效率在线监测方法
CN109899775A (zh) * 2019-01-04 2019-06-18 苏州西热节能环保技术有限公司 一种燃气-蒸汽联合循环机组余热锅炉效率的计算方法
CN110108336A (zh) * 2019-03-29 2019-08-09 苏州西热节能环保技术有限公司 燃气-蒸汽联合循环机组余热锅炉烟气流量的计算方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101576248A (zh) * 2009-06-04 2009-11-11 重庆新仪自控系统工程有限公司 注汽锅炉最优化运行控制系统
CN103984868A (zh) * 2014-05-21 2014-08-13 华中科技大学 一种富氧燃烧锅炉热力性能获取方法

Also Published As

Publication number Publication date
CN111521430A (zh) 2020-08-11

Similar Documents

Publication Publication Date Title
Si et al. Exergy analysis of a 1000 MW double reheat ultra-supercritical power plant
CN102494714B (zh) 电站锅炉效率与燃煤热值、灰分和水分的同步测算方法
CN103646176B (zh) 汽轮机汽封改造后节能效果的综合计算方法
CN108691585B (zh) 一种凝汽式汽轮机低压缸效率的计算方法
CN107201921B (zh) 一种汽轮机热耗率在线监测系统及测量方法
CN110608070B (zh) 汽轮机最小安全流量控制方法
CN109388844B (zh) 低压省煤器节能效果的修正计算方法
CN113095591B (zh) 一种用于火电机组运行参数自寻优的耗差分析方法
Henke et al. Validation of a T100 micro gas turbine steady-state simulation tool
CN106960113B (zh) 一种单轴联合循环发电机组性能验收中机岛与炉岛的责任划分方法
WO2023173624A1 (zh) 一种重型燃气轮机关键部件效能参数软测量方法
CN107340137A (zh) 一种重型燃气轮机透平效率在线监测系统装置及其方法
CN103776502B (zh) 火电机组中低压缸入口热再热蒸汽质量流量实时计量方法
CN108197723B (zh) 煤电机组供电煤耗与污染物排放的优化节能调度方法
CN111521430B (zh) 一种余热锅炉性能试验方法
CN106096310A (zh) 一种抽汽供热机组供热工况热力性能计算方法
CN108182553A (zh) 一种燃煤锅炉燃烧效率在线测量方法
CN104483152B (zh) 非再热回热复合循环机组的热耗率测定方法
CN103697958B (zh) 燃煤机组汽包出口饱和蒸汽质量流量的实时计量方法
CN109932649A (zh) 一种火电机组发电效率监测方法
CN207114187U (zh) 一种重型燃气轮机透平效率在线监测系统装置
CN103728055B (zh) 一种火电机组锅炉炉膛出口烟气能量的实时估计方法
CN114004159B (zh) 一种余热锅炉过热器变工况效能计算方法
Le et al. Comparison of model-driven soft measurement methods for compressor air flow in gas-steam combined cycle power units
Wang et al. Sensitive analysis of energy consumption of operating parameters for coal-fired unit

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