背景技术
国家标准《GB10184-88电站锅炉性能试验规程》公开了煤粉锅炉热效率的计算方法如下:
煤粉锅炉热效率与热损失的关系
锅炉热效率ηb(%)由式(1)计算
ηb=100-(q2+q3+q4+q5+q6) (1)
其中,q2、q3、q4、q5、q6分别为煤粉锅炉的排烟热损失(%)、化学未燃尽热损失(%)、机械未燃尽热损失(%)、锅炉散热损失(%)、灰渣物理热损失(%)。
排烟热损失q2(%)由式(2)计算
其中,V
gy为1kg煤炭燃烧生成的干烟气体积(Nm
3/kg);
为1kg煤炭的烟气中所含水蒸气容积(Nm
3/kg);θ
py、t
0为锅炉排烟温度和空气预热器进口空气温度(℃);c
gy、
为烟气和水蒸气比热(kJ/(℃Nm
3)),根据温度θ
py、t
0从物理性质参数表达式得到。
煤炭燃烧生成的干烟气体积Vgy(Nm3/kg)由式(3)计算
其中,
为煤炭中实际烧掉碳的收到基质量含量百分数(%);S
ar为煤炭中硫的收到基质量含量百分数(%);N
ar为煤炭中氮的收到基质量含量百分数(%);V
0为实际燃烧理论干空气量(Nm
3/kg);α为过剩空气系数。
实际燃烧理论干空气量V0由式(4)计算
其中,Har为收到基氢质量百分数(%);Oar为收到基氧质量百分数(%)。
其中,War为收到基氢质量百分数(%);dk为空气绝对湿度,即每公斤干空气中含水蒸气的千克数(kg/kg干空气),由式(6)计算
其中,φ为按干、湿球温度查得的空气相对湿度,%;pact为就地大气压,Pa;(pb)0为在t0温度下的水蒸气饱和压力(Pa),在0~50℃范围内,按式(7)计算
化学未燃尽热损失q3(%)由式(8)计算
机械未燃尽热损失q4(%)由式(9)计算
其中,33727为碳燃烧生成二氧化碳的发热值(kJ/kg),即1克摩尔碳完全燃烧生成二氧化碳放出393.5kJ的热量;为收到基机械未燃尽碳百分数(%);Qnet为燃料低位发热值(kJ/kg);Aar为收到基灰质量百分数(%);Chz为灰渣含碳率。
锅炉散热损失q5(%)由式(10)计算
其中,
为额定蒸发量下的散热损失(%);D
e为锅炉的额定蒸发量(t/h);D为锅炉热效率测定时的实际蒸发量(t/h)。
灰渣物理热损失q6(%)由式(11)计算
其中,t
lz为由炉膛排出的炉渣温度(℃);t
cjh为由烟道排出的沉降灰温度,取为沉降灰斗上部空间的烟气温度(℃);c
lz、c
fh、c
cjh分别为炉渣、飞灰及沉降灰的比热kJ/(kg℃),由对应的温度利用灰渣比热特性表达式计算得到;α
lz、α
fh、α
cjh分别为炉渣、飞灰及沉降灰占总灰量的份额;
分别为炉渣中碳质量百分数(%)、飞灰中碳质量百分数(%)、沉降灰中碳质量百分数(%)。
锅炉运行调整将直接引起锅炉热效率的变化。例如,适量地增加送风量,碳颗粒周围氧气浓度增加,燃烧更加充分,未燃尽碳
减少,即q
4降低,锅炉热效率η
b上升;同时,由于风量增加,锅炉排烟体积V
gy增大使得q
2上升,锅炉热效率η
b因此又下降。故此,存在着一个最佳的送风量,使得锅炉热效率达到极大值。实现锅炉最佳运行调整的必要条件之一是:能够随时看到每个运行调整动作后的锅炉热效率值。但是,在锅炉热效率计算用到的参数中,有的需要长达4~5个小时的取样化验过程才能得到,例如
H
ar、O
ar、N
ar、S
ar、A
ar、W
ar、C
hz、Q
net,所以,现有技术存在着尚没有解决的困难,致使对锅炉最佳运行调整的实时监控至今尚未实现。
发明内容
本发明基于当代科学长期实践总结出的一个普遍规律:只要某系统是现实存在的实际系统,就必然存在对应的数学模型;本发明的目的是,对现有技术进行改进,提供一种科学合理,节省时间,监测准确的煤粉锅炉热效率和煤质数据实时监测方法。
为实现本发明目的所采取的技术方案是:一种煤粉锅炉热效率和煤质数据实时监测方法,其特征是,它包括以下内容:
〔1〕建立煤粉锅炉系统热效率实时监测数学模型:
煤粉锅炉热效率与热损失的关系
锅炉热效率ηb(%)
ηb=100-(q2+q3+q4+q5+q6) (1)
其中,q2、q3、q4、q5、q6分别为煤粉锅炉的排烟热损失(%)、化学未燃尽热损失(%)、机械未燃尽热损失(%)、锅炉散热损失(%)、灰渣物理热损失(%);
排烟热损失q2(%)
其中,V
gy为1kg煤炭燃烧生成的干烟气体积(Nm
3/kg),
为1kg煤炭的烟气中所含水蒸气容积(Nm
3/kg),θ
py、t
0为锅炉排烟温度和空气预热器进口空气温度(℃),c
gy、
为烟气和水蒸气比热(kJ/(℃Nm
3)),根据温度θ
py、t
0从物理性质参数表达式得到;
化学未燃尽热损失q3(%)
机械未燃尽热损失q4(%)
其中,33727为碳燃烧生成二氧化碳的发热值(kJ/kg),即1克摩尔碳完全燃烧生成二氧化碳放出393.5kJ的热量;
为收到基机械未燃尽碳百分数(%);Q
net为燃料低位发热值(kJ/kg);A
ar为收到基灰质量百分数(%);C
hz为灰渣含碳率;
锅炉散热损失q5(%)
其中,
为额定蒸发量下的散热损失(%),D
e为锅炉的额定蒸发量(t/h),D为锅炉热效率测定时的实际蒸发量(t/h);
灰渣物理热损失q6(%)
其中,t
lz为由炉膛排出的炉渣温度(℃);t
cjh为由烟道排出之沉降灰温度,取为沉降灰斗上部空间的烟气温度(℃);c
lz、c
fh、c
cjh分别为炉渣、飞灰及沉降灰的比热kJ/(kg℃),由对应的温度利用灰渣比热特性表达式计算得到;α
lz、α
fh、α
cjh分别为炉渣、飞灰及沉降灰占总灰量的份额;
分别为炉渣中碳质量百分数(%)、飞灰中碳质量百分数(%)、沉降灰中碳质量百分数(%);
理论空气量
每公斤煤炭的理论空气容积V0(Nm3/kg)
V0=0.0889927(Car+0.375344Sar)+0.265128Har-0.0334029Oar (12)
其中,Car、Har、Oar、Sar分别为煤炭收到基碳、氢、氧、硫的质量百分数(%),0.0889927为式(17)的推导结果,0.0889927×0.375344为式(23)的推导结果,0.265128为式(29)的推导结果,0.0334029为每kg氧的空气标准状态容积的百分之一;
烟气量
其中,Nar为煤炭收到基氮的质量百分数(%);
每公斤煤炭中燃烧碳生成二氧化碳的标准状态0℃760mmHg下的容积
(Nm
3/kg)
其中,
为煤炭收到基生成二氧化碳的碳质量百分数(%),二氧化碳标准状态密度为
每公斤煤炭中的碳燃烧生成一氧化碳的标准状态(0℃760mmHg)容积
(Nm
3/kg)
其中,
为煤炭收到基生成一氧化碳的碳质量百分数(%),一氧化碳标准状态(0℃760mmHg)密度为
每公斤煤炭中的硫燃烧生成二氧化硫的标准状态容积
(Nm
3/kg)
其中,Sar为煤炭收到基硫的质量百分数(%),二氧化硫标准状态密度为
燃烧生成二氧化碳的收到基碳质量百分数
每公斤煤炭燃烧前后的二氧化碳物质守恒关系
或者带入式(15) (25)
其中,CO2k、CO2分别为空气中二氧化碳气容积百分数(%)、烟气中二氧化碳气容积百分数(%);
燃烧生成一氧化碳的收到基碳质量百分数
这里忽略了空气中极少量的一氧化碳,每公斤煤炭燃烧前后的一氧化碳物质守恒关系
或者带入式(19) (26)
其中,CO为烟气中一氧化碳气容积百分数(%);
每公斤煤炭完全燃烧理论需氧量与实际燃烧耗氧量的容积差
其中,Oar为煤炭收到基氧质量百分数(%),CO为烟气中一氧化碳气容积百分数(%),CO2为烟气中二氧化碳气容积百分数(%),SO2为烟气中二氧化硫气容积百分数(%);
每公斤煤炭完全燃烧理论需氧量与实际燃烧耗氧量的容积差
还等于一氧化碳和机械未燃尽碳
的燃烧需氧量,从而得到机械未燃尽碳
的显函数
其中,1.8643978为式(16)的推导结果;
煤炭收到基碳质量百分数Car,碳燃烧生成二氧化碳、一氧化碳和没有燃烧部分的数量关系
其中,
分别为煤炭收到基生成二氧化碳的碳质量百分数(%)、生成一氧化碳的碳质量百分数(%)、机械未燃尽碳质量百分数(%)、煤炭收到基碳质量百分数(%);
收到基氢和水质量百分数Har、War,烟气中的水份,即氢的燃烧产物、煤炭中含有的水份和空气中含有的水份
其中,H2O为烟气中水蒸气容积百分数(%),H2Ok为空气中水蒸气容积百分数(%),War为煤炭收到基水质量百分数(%),Har为煤炭收到基氢质量百分数(%),0.11116039为式(30)的推导结果;
煤炭收到基氧质量百分数Oar根据氧的平衡方程式得到
或者
或者
煤炭收到基氮质量百分数Nar根据氮的平衡方程式得到
其中,
为烟气中氮氧化物容积(Nm
3/kg),
为烟气中氮氧化物含量(mg/Nm
3),1.59kg/Nm
3为二氧化氮(NO
2)标准状态密度,仅考虑了氮氧化物中的二氧化氮;
煤炭收到基硫质量百分数Sar根据硫的平衡方程式得到
其中,SO2k为空气二氧化硫体积百分数(%);
煤炭收到基灰质量百分数Aar由收到基定义计算
Aar=100-(Car+Har+Oar+Nar+Sar+War) (39)
空气预热器阻力,通过实测烟气在空气预热器中的阻力可以建立锅炉燃煤量、机械未燃尽热损失与煤质数据之间的关系;
烟气密度ρy(kg/Nm3)
其中:μy为烟气含灰浓度(kg/Nm3),为烟气的水蒸气容积(Nm3/kg),O2为烟气中氧气容积百分数(%),CO2为烟气中二氧化碳气容积百分数(%),H2O为烟气中水蒸气容积百分数(%),SO2为烟气中二氧化硫气容积百分数(%),N2为烟气中氮气容积百分数(%),这些参数均由实测得到;
烟气在空气预热器中的流速W(m/s)
其中,Bb为锅炉燃煤量(t/h),q4为机械未完全燃烧热损失(%),tkr为空气预热器烟气平均温度(℃),Ak为空气预热器烟气流通断面积(m2);
烟气在空气预热器中的阻力Δh(Pa)
其中,Δh为烟气在空气预热器中的阻力,现场实测;l为空气预热器烟气流道长度(m),由制造厂提供;dd为烟气在空气预热器中的流道当量直径(m),由制造厂提供;v为烟气粘度(m2/s),根据烟气温度由烟气特性程序计算;(a+b)为计入波纹高度的平均公差(mm),由制造厂提供;s为空气预热器板型结构特性参数(mm),由制造厂提供;
煤炭低位发热值Qnet
Qnet=339Car+1105.1Har+108.8Sar-108.8Oar-25.1War (43)
其中,Car、Har、Sar、Oar、War分别为煤炭收到基成分碳、氢、硫、氧、水(%);
锅炉燃煤量Bb
BbQnetηb=Qre (44)
其中,Bb为锅炉燃煤量(kg/h);Qre(kJ/h)为锅炉有效热量,通过实测蒸汽流量和焓、给水流量和焓、排污流量和焓计算得到;
过量空气系数α
其中,O2、CO为烟气的氧气容积百分数、烟气的一氧化碳容积百分数;
或者
其中,对于煤粉锅炉α
lz、α
fh、α
cjh分别为0.1、0.85、0.05,基本稳定;炉渣中碳质量百分数
沉降灰中碳质量百分数
取为定期取样化验值,由于这两部分灰渣份额较少,所以误差很少;
〔2〕使用计算机C语言程序进行原始数据的读取和数学模型的求解,其解算顺序是:根据烟气分析容积百分数CO、O
2用式(45)计算过量空气系数α,根据实际蒸发量(t/h)用式(10)计算散热热损失q
5(%);根据能够实时检测的锅炉运行数据求解式(1)、(2)、(8)、(9)、(11)、(12)、(13)、(15)、(19)、(22)、(24)、(25)、(26)、(31)、(33)、(34)、(35)、(36)、(37)、(38)、(39)、(40)、(41)、(42)、(43)、(44)和(46)数学模型组成的方程组,得到每公斤燃料的理论空气容积V
0(Nm
3/kg)、燃烧产物理论氮气容积
(Nm
3/kg)、烟气的二氧化碳容积
(Nm
3/kg)、一氧化碳容积
(Nm
3/kg)、二氧化硫容积
(Nm
3/kg)、每公斤燃料燃烧生成的干烟气容积
(Nm
3/kg)、每公斤煤炭完全燃烧理论需氧量与实际燃烧耗氧量的容积差
(Nm
3/kg)、煤炭收到基机械未燃尽碳质量百分数
(%)、煤炭生成二氧化碳的收到基碳质量百分数
(%)、煤炭生成一氧化碳的收到基碳质量百分数
(%)、煤炭收到基质量百分数碳C
ar(%)、氢H
ar(%)、氧O
ar(%)、氮N
ar(%)、硫S
ar(%)、灰A
ar(%)、水W
ar(%)、烟气密度ρ
y(kg/Nm
3)、烟气在空气预热器中的流速W(m/s)、煤炭低位发热值Q
net(kJ/kg)、排烟热损失q
2(%)、化学未完全燃烧热损失q
3(%)、机械未完全燃烧热损失q
4(%)、灰渣物理热损失q
6(%)、锅炉热效率η
b(%)、飞灰中碳质量百分数
(%)和燃煤量B
b(t/h)不能够实时检测的锅炉运行参数,其中,锅炉热效率η
b是锅炉最佳运行调整实时监控需要的关键参数,燃煤量B
b是锅炉运行重要的经济指标。其他方程为中间计算过程式或化学反应方程式。能够实时检测的锅炉运行数据包括:烟气分析一氧化碳容积百分数CO(%)、烟气分析二氧化碳容积百分数CO
2(%)、烟气分析二氧化硫容积百分数SO
2(%)、烟气分析氮气容积百分数N
2(%)、烟气分析水蒸气容积百分数H
2O(%)、烟气分析氧气容积百分数O
2(%)、锅炉排烟温度θ
py(℃)、空气预热器进口空气温度t
0(℃)、由炉膛排出的炉渣温度t
lz(℃)、由烟道排出之沉降灰温度t
cjh(℃)、空气中水蒸气容积百分数H
2O
k(%)、烟气中氮氧化物含量
(mg/Nm
3)、空气二氧化硫体积百分数SO
2k(%)、烟气含灰浓度μ
y(kg/Nm
3)、空气预热器烟气平均温度t
kr(℃)、实际蒸发量(t/h)、烟气在空气预热器中的阻力Δh(Pa)和锅炉有效热量Q
re(kJ/h)。烟气分析容积百分数以干烟气为基数,即CO+CO
2+SO
2+O
2+N
2=100。
本发明的煤粉锅炉热效率和煤质数据实时监测方法科学合理,节省时间,监测准确;实现了煤粉锅炉热效率和煤质数据实时监测,即,通过燃烧化学反应物质守恒定律和流体力学原理建立的映射关系,将能够实时检测的锅炉运行数据映射为不能够实时检测锅炉运行参数,实现了煤粉锅炉热效率和煤质数据实时测量;可操作性强,经济效益和社会效益显著。
具体实施方式
下面利用具体实施方式对本发明作进一步说明。
一种煤粉锅炉热效率和煤质数据实时监测方法,它包括以下内容:
〔1〕建立煤粉锅炉系统热效率实时监测数学模型:
煤粉锅炉热效率与热损失的关系
锅炉热效率ηb(%)由式(1)计算
ηb=100-(q2+q3+q4+q5+q6) (1)
其中,q2、q3、q4、q5、q6分别为煤粉锅炉的排烟热损失(%)、化学未燃尽热损失(%)、机械未燃尽热损失(%)、锅炉散热损失(%)、灰渣物理热损失(%)。
排烟热损失q2(%)由式(2)计算
其中,V
gy为1kg煤炭燃烧生成的干烟气体积(Nm
3/kg);
为1kg煤炭的烟气中所含水蒸气容积(Nm
3/kg);θ
py、t
0为锅炉排烟温度和空气预热器进口空气温度(℃);c
gy、
为烟气和水蒸气比热(kJ/(℃Nm
3)),根据温度θ
py、t
0从物理性质参数表达式得到。
化学未燃尽热损失q3(%)由式(8)计算
机械未燃尽热损失q4(%)由式(9)计算
其中,33727为碳燃烧生成二氧化碳的发热值(kJ/kg),即1克摩尔碳完全燃烧生成二氧化碳放出393.5kJ的热量;为收到基机械未燃尽碳百分数(%);Qnet为燃料低位发热值(kJ/kg);Aar为收到基灰质量百分数(%);Chz为灰渣含碳率。
锅炉散热损失q5(%)由式(10)计算
其中,
为额定蒸发量下的散热损失(%);D
e为锅炉的额定蒸发量(t/h);D为锅炉效率测定时的实际蒸发量(t/h)。
灰渣物理热损失q6(%)由式(11)计算
其中,t
lz为由炉膛排出的炉渣温度(℃);t
cjh为由烟道排出之沉降灰温度,取为沉降灰斗上部空间的烟气温度(℃);c
lz、c
fh、c
cjh分别为炉渣、飞灰及沉降灰的比热kJ/(kg℃),由对应的温度计算得到;α
lz、α
fh、α
cjh分别为炉渣、飞灰及沉降灰占总灰量的份额;
分别为炉渣中碳质量百分数(%)、飞灰中碳质量百分数(%)、沉降灰中碳质量百分数(%)。
理论空气量
每公斤煤炭的理论空气容积V0(Nm3/kg)由式(12)计算。
V0=0.0889927(Car+0.375344Sar)+0.265128Har-0.033402Oar (12)
其中,Car、Har、Oar、Sar分别为煤炭收到基碳、氢、氧、硫的质量百分数(%),0.0889927为式(17)的推导结果,0.0889927×0.375344为式(23)的推导结果,0.265128为式(29)的推导结果,0.0334029为每kg氧的空气标准状态容积的百分之一。
烟气量
燃烧产物理论氮气容积
(Nm
3/kg)由式(13)计算。
其中,Nar为煤炭收到基氮的质量百分数(%)。
碳燃烧生成二氧化碳,化学反应方程式如式(14)所示
C+O2=CO2 (14)
碳、氧、二氧化碳的分子量分别为12.011、32.000、44.011,从而,反应过程各物质的质量比例为12.011∶32.0∶44.011,也就是每公斤碳燃烧生成二氧化碳44.011/12.011=3.66422446公斤。每公斤煤炭中燃烧碳生成二氧化碳的标准状态(0℃760mmHg)下的容积
(Nm
3/kg)由式(15)表示
其中,
为煤炭收到基生成二氧化碳的碳质量百分数(%);二氧化碳标准状态密度为
从反应过程各物质的质量比例关系也得到,每公斤碳燃烧生成二氧化碳需要氧32.000/12.011=2.66422446公斤。每公斤煤炭中的机械未燃尽碳,如果燃烧生成二氧化碳的话,需要的氧气标准状态(0℃760mmHg)容积
(Nm
3/kg)由式(16)表示
其中,
为煤炭收到基机械未燃尽碳的质量百分数(%);氧气标准状态密度为
碳燃烧生成二氧化碳需要的空气量
(Nm
3/kg)由式(17)表示
根据化学反应方程式和摩尔定律可知,在二氧化碳反应过程中,二氧化碳容积在0.58%的相对误差下等于氧气容积,参见式(15)和式(16)。
碳燃烧生成一氧化碳,化学反应方程式如式(18)所示
2C+O2=2CO (18)
碳的分子量为12.011,氧的分子量为32.000,一氧化碳的分子量为28.011,从而反应过程各物质的质量比例为2×12.011∶32.0∶2×28.011,也就是每公斤碳燃烧生成一氧化碳28.011/12.011=2.33211223公斤,每公斤煤炭中的碳燃烧生成一氧化碳的标准状态(0℃760mmHg)容积
(Nm
3/kg)由式(19)表示
其中,
为煤炭收到基生成一氧化碳的碳质量百分数(%);一氧化碳标准状态密度为
从分子量关系也得到,每公斤碳燃烧生成一氧化碳需要氧32.000/(2×12.011)=1.33211公斤。在每公斤煤炭中,生成一氧化碳需要的氧气标准状态(0℃760mmHg)容积
(Nm
3/kg)由式(20)表示
根据化学反应方程式和摩尔定律可知,在一氧化碳反应过程中,一氧化碳容积在0.069%的相对误差下是氧气容积的两倍,参见式(19)和式(20)。
硫燃烧生成二氧化硫,化学反应方程式如式(21)所示
S+O2=SO2 (21)
硫的分子量为32.06,氧的分子量为32.000,二氧化硫的分子量为64.06,反应过程各物质的质量比例为32.06∶32.00∶64.06,也就是每公斤硫燃烧生成二氧化硫64.06/32.06=1.9981公斤,每公斤煤炭中的硫燃烧生成二氧化硫的标准状态容积(Nm3/kg)由式(22)表示
其中,S
ar为煤炭收到基硫的质量百分数(%),二氧化硫标准状态密度为
每公斤硫燃烧生成二氧化硫需要氧32.000/32.06=0.9981285公斤,即需要空气容积
(Nm
3/kg)由式(23)表示
根据化学反应方程式和摩尔定律可知,在二氧化硫反应过程中,二氧化硫气体容积等于氧气容积。
每公斤煤炭燃烧生成的干烟气量(Nm3/kg)由式(24)表示
二氧化碳的物质守恒为:空气中的二氧化碳+燃烧产生的二氧化碳=烟气中的二氧化碳,如式(25)表示了每公斤煤炭燃烧前后的二氧化碳物质守恒关系
或者带入式(15) (25)
其中,CO2k、CO2分别为空气中二氧化碳气容积百分数(%)、烟气中二氧化碳气容积百分数(%)。
考虑到空气中一氧化碳含量极少可以忽略,一氧化碳的物质守恒为:燃烧产生的一氧化碳量=烟气中的一氧化碳量,如式(26)表示了每公斤煤炭燃烧前后的一氧化碳物质守恒关系
或者带入式(19) (26)
其中,CO为烟气中一氧化碳气容积百分数(%)。
氢燃烧生成水,化学反应方程式如式(27)所示
2H2+O2=2H2O (27)
氢的分子量为2.0158,氧的分子量为32.000,水的分子量为18.0158,反应过程各物质的质量比例为2×2.0158∶32.00∶2×18.0158,也就是每公斤氢燃烧生成水需要消耗氧32.000/(2×2.0158)=7.9372954公斤,每公斤氢燃烧生成水2×18.0158/(2×2.0158)=8.9372954公斤。
每公斤煤炭氢燃烧消耗氧容积如式(28)所示
每公斤煤炭氢燃烧消耗空气容积如式(29)所示
其中,H
ar为煤炭收到基氢质量百分数(%),
为水蒸气标准态密度,注:一摩尔水为18克,变为蒸气在标准状况下是22.4升,所以密度为18/22.4=0.80357克/升。
每公斤煤炭完全燃烧理论需氧量与实际燃烧耗氧量的容积差
(Nm
3/kg)用式(31)表示
其中,Oar为煤炭收到基氧质量百分数(%),CO为烟气中一氧化碳气容积百分数(%),CO2为烟气中二氧化碳气容积百分数(%),SO2为烟气中二氧化硫气容积百分数(%)。
式(31)等号右边,第一项是煤炭自带的氧气,第二项是风机送入锅炉的氧气。第三、四项分别是燃烧生成二氧化碳、二氧化硫消耗的氧气。第五项是燃烧生成一氧化碳消耗的氧气。根据摩尔定律,相同摩尔数的不同气体在标准状态下体积相同,故而采用生成物容积表示消耗氧气容积。第六项是氢燃烧生成水消耗的氧气。
完全燃烧理论需氧量与实际燃烧耗氧量的容积差
还可以用机械未燃尽碳
和一氧化碳气体的燃尽需氧量表示,带入式(16),如式(32)表示
煤炭收到基碳质量百分数Car
进入锅炉的碳经过炉膛的燃烧过程分成三部分:其一,燃烧生成二氧化碳,占绝大部分;其二,燃烧生成一氧化碳,占很少部分;其三,还有少量没有燃烧。这三部分的数量关系用式(34)表示
其中,
分别为煤炭收到基生成二氧化碳的碳质量百分数(%)、生成一氧化碳的碳质量百分数(%)、机械未燃尽碳质量百分数(%)、煤炭收到基碳质量百分数(%)。
收到基氢和水质量百分数Har、War
烟气中的水份有三个来源,即氢的燃烧产物如式(30)、煤炭中含有的水份和空气中含有的水份,如式(35)所示
其中,H2O为烟气中水蒸气容积百分数(%),H2Ok为空气中水蒸气容积百分数(%),War为煤炭收到基水质量百分数(%),Har为煤炭收到基氢质量百分数(%),0.11116039为式(30)的推导结果。
收到基氧质量百分数Oar
煤炭收到基氧质量百分数根据氧的平衡方程式得到,如式(36)所示。氧平衡方程式的建立原则是:向锅炉提供的氧量等于燃烧消耗的氧量与烟气中的氧量之和。
或者
或者
在式(36)中,等号左边为向锅炉提供的氧量,由两部分组成,即煤炭自带的氧量和送风机送入锅炉的氧量;等号右边,第一、二、三、四项为燃烧消耗的氧量,第五项为烟气中的氧量。氮氧化物的生成也消耗氧,但是相比之下很少,对方程组的求解结果几乎没有影响,这里忽略了。2.66422446为式(16)的推导结果;1.33211223为式(20)的推导结果;0.9981285为式(21)的推导结果;7.9372953666为式(28)的推导结果。
收到基氮质量百分数Nar
煤炭收到基氮质量百分数Nar根据氮的平衡方程式得到,如式(37)所示
(37)
其中,
为烟气中氮氧化物容积(Nm
3/kg),
为烟气中氮氧化物含量(mg/Nm
3),1.59kg/Nm
3为二氧化氮(NO
2)标准状态密度,仅考虑了氮氧化物中的二氧化氮。
收到基硫质量百分数Sar
煤炭收到基硫质量百分数Sar根据硫的平衡方程式得到,利用式(22),如式(38)所示
其中,SO2k为空气二氧化硫体积百分数(%)。
收到基灰质量百分数Aar
煤炭收到基灰质量百分数由收到基定义计算,如式(39)所示
Aar=100-(Car+Har+Oar+Nar+Sar+War) (39)
空气预热器阻力
通过实测烟气在空气预热器中的阻力可以建立锅炉燃煤量、机械未燃尽热损失与煤质数据之间的关系。
烟气密度ρy(kg/Nm3)由式(40)计算
其中,μ
y为烟气含灰浓度(kg/Nm
3),
为烟气的水蒸气容积(Nm
3/kg),O
2为烟气中氧气容积百分数(%),CO
2为烟气中二氧化碳气容积百分数(%),H
2O为烟气中水蒸气容积百分数(%),SO
2为烟气中二氧化硫气容积百分数(%),N
2为烟气中氮气容积百分数(%)。这些参数可以实测得到。
烟气在空气预热器中的流速W(m/s)由式(41)计算
其中,Bb为锅炉燃煤量(t/h),q4为机械未完全燃烧热损失(%),tkr为空气预热器烟气平均温度(℃),Ak为空气预热器烟气流通断面积(m2)。
烟气在空气预热器中的阻力Δh(Pa)由式(42)
其中,Δh为烟气在空气预热器中的阻力,现场实测;l为空气预热器烟气流道长度(m),由制造厂提供;dd为烟气在空气预热器中的流道当量直径(m),由制造厂提供;v为烟气粘度(m2/s),根据烟气温度由烟气特性程序计算;(a+b)为计入波纹高度的平均公差(mm),由制造厂提供;s为空气预热器板型结构特性参数(mm),由制造厂提供。
煤炭低位发热值Qnet
煤炭低位发热值Qnet由式(43)计算
Qnet=339Car+1105.1Har+108.8Sar-108.8Oar-25.1War (43)
其中,Car、Har、Sar、Oar、War分别为煤炭收到基成分碳、氢、硫、氧、水(%)。
锅炉燃煤量Bb
锅炉燃煤量Bb由式(44)计算
BbQnetηb=Qre (44)
其中,Bb为锅炉燃煤量(kg/h);Qre(kJ/h)为锅炉有效热量,通过实测蒸汽流量和焓、给水流量和焓、排污流量和焓计算得到。
过量空气系数α
过量空气系数α由式(45)计算
其中,O2、CO为烟气的氧气容积百分数、烟气的一氧化碳容积百分数。
飞灰中碳质量百分数(%)
或者
其中,对于煤粉锅炉α
lz、α
fh、α
cjh分别为0.1、0.85、0.05,基本稳定;炉渣中碳质量百分数
沉降灰中碳质量百分数
取为定期取样化验值,由于灰渣份额较少,所以误差很少。
上述数学模型即包含了能够实时检测的锅炉运行数据,又包含了不能够实时检测的锅炉运行参数:
〔2〕使用计算机C语言程序进行原始数据的读取和数学模型的求解,其解算顺序是:根据烟气分析容积百分数CO、O
2用式(45)计算过量空气系数α,根据实际蒸发量(t/h)用式(10)计算散热热损失q
5(%);根据能够实时检测的锅炉运行数据求解由煤粉锅炉系统热效率实时监测数学模型组成的方程组,方程组的方程个数等于未知数个数,属于适定方程组,有唯一解,得到每公斤燃料的理论空气容积V
0(Nm
3/kg)、燃烧产物理论氮气容积
(Nm
3/kg)、烟气的二氧化碳容积
(Nm
3/kg)、一氧化碳容积
(Nm
3/kg)、二氧化硫容积
(Nm
3/kg)、每公斤燃料燃烧生成的干烟气容积
(Nm
3/kg)、每公斤煤炭完全燃烧理论需氧量与实际燃烧耗氧量的容积差
(Nm
3/kg)、煤炭收到基机械未燃尽碳质量百分数
(%)、煤炭生成二氧化碳的收到基碳质量百分数
(%)、煤炭生成一氧化碳的收到基碳质量百分数
(%)、煤炭收到基质量百分数碳C
ar(%)、氢H
ar(%)、氧O
ar(%)、氮N
ar(%)、硫S
ar(%)、灰A
ar(%)、水W
ar(%)、烟气密度ρ
y(kg/Nm
3)、烟气在空气预热器中的流速W(m/s)、煤炭低位发热值Q
net(kJ/kg)、排烟热损失q
2(%)、化学未完全燃烧热损失q
3(%)、机械未完全燃烧热损失q
4(%)、灰渣物理热损失q
6(%)、锅炉热效率η
b(%)、飞灰中碳质量百分数
(%)和燃煤量B
b(t/h)不能够实时检测的锅炉运行参数,其中,锅炉热效率η
b是锅炉最佳运行调整实时监控需要的关键参数,燃煤量B
b是锅炉运行重要的经济指标。能够实时检测的锅炉运行数据包括:烟气分析一氧化碳容积百分数CO(%)、烟气分析二氧化碳容积百分数CO
2(%)、烟气分析二氧化硫容积百分数SO
2(%)、烟气分析氮气容积百分数N
2(%)、烟气分析水蒸气容积百分数H
2O(%)、烟气分析氧气容积百分数O
2(%)、锅炉排烟温度θ
py(℃)、空气预热器进口空气温度t
0(℃)、由炉膛排出的炉渣温度t
lz(℃)、由烟道排出之沉降灰温度t
cjh(℃)、空气中水蒸气容积百分数H
2O
k(%)、烟气中氮氧化物含量
(mg/Nm
3)、空气二氧化硫体积百分数SO
2k(%)、烟气含灰浓度μ
y(kg/Nm
3)、空气预热器烟气平均温度t
kr(℃)、实际蒸发量(t/h)、烟气在空气预热器中的阻力Δh(Pa)和锅炉有效热量Q
re(kJ/h)。烟气分析容积百分数以干烟气为基数,即CO+CO
2+SO
2+O
2+N
2=100。
将能够实时检测的锅炉运行数据通过变送器和数字转换器送入锅炉监控计算机指定地址;采用计算机C语言将本发明提供的数学模型解算过程编制成应用程序为本领域技术人员所熟悉的技术。
一个具体实施例的数学模型计算与对比结果的验证:
这里针对配200MW汽轮发电机组的670t/h煤粉锅炉进行了电站锅炉热效率实时测量数学模型验证计算,其结果列于表1中。从表1中的锅炉热效率对比情况可以看到,本发明的方法可以投入到锅炉运行中实际应用,实时获取锅炉热效率,解决了实现锅炉最佳运行调整实时监控的技术难题之一。
在表1中,由数学模型计算的参数与热效率试验得到的参数有一定的差别,经分析,差别来自以下原因:
热效率试验参数是一定时间内的平均参数,数学模型计算参数是瞬时参数,是不具备对比充分条件的理由之一。
采用平均实测参数求取瞬时参数,是不具备对比充分条件的理由之二。
由于热效率试验实测参数不全,若干原始数据需要估计,如空气二氧化碳体积百分数、空气二氧化硫体积百分数、空气水蒸气体积百分数、氮氧化物。虽然估计,基本不影响数学模型的解算结果,因为空气成分在大部分地区基本上是常数。如果本方法使用在青藏高原,自然不可估计空气成分。
表1锅炉热效率在线测量程序验证计算对比表
利用上述数学模型的解算方法,采用C语言程序能够在0.4秒内完成原始数据读取和数学模型的求解,计算结果与根据《电站锅炉性能试验规程GB10184-88》进行的锅炉热效率试验结果一致。