第四章循环流化床锅炉炉内传热计算.docx

1、第四章循环流化床锅炉炉内传热计算第四章 循环流化床锅炉炉内传热计算循环流化床锅炉炉膛中的传热是一个复杂的过程，传热系数的计算精度直接影响了受热面设计时的布置数量，从而影响锅炉的实际出力、蒸汽参数和燃烧温度。正确计算燃烧室受热面传热系数是循环流化床锅炉设计的关键之一，也是区别于煤粉炉的重要方面。 随着循环流化床燃烧技术的日益成熟，有关循环流化床锅炉的炉膛传热计算思想和方法的研究也在迅速发展。许多著名的循环流化床制造公司和研究部门在此方面也做了大量的工作，有的已经形成商业化产品使用的设计导则。但由于技术保密的原因，目前国内外还没有公开的可以用于工程使用的循环流化床锅炉炉膛传热计算方法，因此对它的研

2、究具有重要的学术价值和实践意义。清华大学对CFB锅炉炉膛传热作了深入的研究，长江动力公司、华中理工大学、浙江大学等单位也对CFB锅炉炉膛中的传热过程进行了有益的探索。根据已公开发表的文献报导，考虑工程上的方便和可行，本章根椐清华大学提出的方法，进一步分析整理，作为我们研究的基础。为了了解CFB锅炉传热计算发展过程，也参看了巴苏的传热理论和计算方法，浙江大学和华中理工大学的传热计算与巴苏的相近似。4.1 清华的传热理论及计算方法4.1.1 循环流化床传热分析CFB锅炉与煤粉锅炉的显著不同是CFB锅炉中的物料(包括煤灰、脱硫添加剂等)浓度Cp大大高于煤粉炉，而且炉内各处的浓度也不一样，它对炉内传热

3、起着重要作用。为此首先需要计算出炉膛出口处的物料浓度Cp，此处浓度可由外循环倍率求出。而炉膛不同高度的物料浓度则由内循环流率决定，它沿炉膛高度是逐渐变化的，底部高、上部低。近壁区贴壁下降流的温度比中心区温度低的趋势，使边壁下降流减少了辐射换热系数；水平截面方向上的横向搅混形成良好的近壁区物料与中心区物料的质交换，同时近壁区与中心区的对流和辐射的热交换使截面方向的温度趋于一致，综合作用的结果近壁区物料向壁面的辐射加强，总辐射换热系数明显提高。在计算水冷壁、双面水冷壁、屏式过热器和屏式再热器时需采用不同的计算式。物料浓度Cp对辐射传热和对流传热都有显著影响。燃烧室的平均温度是床对受热面换热系数的另

4、一个重要影响因素。床温的升高增加了烟气辐射换热并提高烟气的导热系数。虽然粒径的减小会提高颗粒对受热面的对流换热系数，在循环流化床锅炉条件下，燃烧室内部的物料颗粒粒径变化较小，在较小范围内的粒径变化时换热系数的变化不大，在进行满负荷传热计算时可以忽略，但在低负荷传热计算时，应该考虑小的颗粒有提高传热系数的能力。炉内受热面的结构尺寸，如鳍片的净宽度、厚度等，对平均换热系数的影响也是非常明显的。鳍片宽度对物料颗粒的团聚产生影响；另一方面，宽度与扩展受热面的利用系数有关。根据实验研究，可以归纳出循环流化床锅炉燃烧室受热面传热系数的计算方法。CFB锅炉炉膛受热面的吸热量按下式计算： (4-1)式中 Q传

5、热量，W；K基于烟气侧总面积的传热系数，W/m2K；T温差，K；H烟气侧总面积，m2。4.1.2 受热面结构尺寸对传热的影响传热系数K按式(4-2)计算，其中分母包括四部分热阻：烟气侧热阻；工质侧热阻和受热面本身热阻；以及附加热阻as。 (4-2)式中 烟气侧向壁面总表面的名义换热系数，W/m2K；f工质侧换热系数，W/m2K，可按苏1973年热力计算标准求取；Ht烟气侧总面积，m2；Hf工质侧总面积，m2；as附加热阻， m2K/W；1管子厚度，m；受热面金属导热系数，W/m2K； (4-3)式中 P鳍片面积系数，；Hfin鳍片面积，m2；Ht受热面外部面积，m2。 (4-4)图4-1 烧室

6、受热面结构简图s, d管子节距、外径，m，见图4-1。鳍片利用系数， (4-5)式中 与受热面受热情况、膜式壁鳍片结构尺寸和材料等有关，可表示为 (4-6)式中 N受热情况，单面受热N=1，双面受热N=2；h实际鳍片高度 图4-1 炉膛受热面结构简图 (4-7)鳍片厚度，m；s受热面污染系数，取为0.0005；h折算高度，m： (4-8)h”有效高度，m： (4-9)根据实验和运行数据，可得到鳍片宽度系数与结构尺寸的关系： (4-10)b烟气侧换热系数，见式(4-15)：as附加热阻，在计算耐火材料涂层受热面时考虑： (4-11)a受热面耐火层厚度，m；a受热面耐火层导热系数，W/mK，按式(

7、4-12)计算： (4-12)式中 a0、a1系数；耐火层平均温度，K，按式(4-13)计算： (4-13)式中 Tb烟气侧温度，K；Tw受热面壁面温度，K，见式(4-17)：受热面外内面积比为 (4-14)式中 1管壁厚度，m；s管节距，m；鳍片厚度，m。4.1.3 CFB锅炉烟气侧换热系数b炉膛烟气物料两相混合物向壁面的换热包括对流和辐射两部分，按两者的线性叠加，则有 (4-15)式中 r辐射换热系数，W/m2K，见式(4-16)：c对流换热系数，W/m2K，见式(4-26)： (4-16)式中 Boltzmann常数；Tw水冷壁管壁温度，按式(4-17)计算： (4-17)式中 Tf受热

8、面内工质温度，K。水冷壁管壁内外侧温差 (4-18)式中 Tb烟气侧温度，K；Tf受热面内工质温度，K；N受热情况，1或2；w导热影响系数，w=0.2+0.007 (4-19)式中 金属导热系数，W/m2K；壁面与烟气侧的系统黑度可写作式(4-20)的形式： (4-20)式中 b烟气侧黑度，按式(4-21)计算：壁面黑度，一般为0.50.8。在气固两相中，烟气侧黑度包括颗粒黑度和烟气黑度两部分： (4-21)式中 p固体物料黑度，由式(4-22)计算： (4-22)式中 B系数，各向同性反射时为0.5，漫反射颗粒为，本文中取为；物料表面平均黑度，与固体颗粒的浓度有关，可表示为 (4-23)式中

9、 C常数；C为0.10.2；Cp物料空间浓度，kg/m3。g烟气黑度，由式(4-24)计算： (4-24)烟气辐射减弱系数k可按下式简单计算： (4-25)式中，烟气中水蒸气份额；烟气中三原子气体份额；sg烟气辐射厚度，近似为下降流厚度，m。对流换热系数由烟气对流和颗粒对流两部分组成，即 (4-26)式中 烟气对流换热系数，W/m2K，计算见式(4-27)；颗粒对流换热系数，计算见式(4-28)。 (4-27)式中 烟气对流系数，46J/m3K；vf烟气速度，m/s。 (4-28)式中 vf烟气速度，m/s，该项为颗粒对流强度与颗粒粒径的直接修正；初始流态条件下颗粒对流理论换热系数，其值与颗粒

10、的粒度、温度、受热面布置有关；颗粒对流系数，按式(4-29)计算： (4-29)式中 Cpc颗粒系数，0.010.02；Cp炉膛局部物料浓度，kg/m3；n1常数，0.851.25。根据第二章中上部快速床的分析，则受热面所在位置的浓度与其高度位置密切相关，用于传热的平均浓度关联到受热面的平均高度，则双面水冷壁、屏过、屏再局部物料浓度Cp按式(4-30)计算： (4-30)图4-2 特征携带量式中 Cpp实际温度下炉膛出口处特征物料浓度，kg/m3；该数值可以根据图4-2选定，并根据经验予以修正。Hlt炉膛总高度，m；Hpz双面水冷壁屏再或屏过总高度，m。水冷壁物料浓度Cp按式(4-31)计算：

11、 (4-31)式中 h1炉膛下部冷灰斗锥体计算高度(从布风板算起)，m； dh梯形段上直段耐火层高度，m。4.1.4 按清华方法对一台440 t/h 贫煤CFB锅炉的计算用清华方法对按某国外引进程序设计的锅炉输入数据及计算结果进行了分析校核，以便了解影响传热的因素和影响关系。该炉为燃烧贫煤的440 t/h CFB锅炉，100%，、50%负荷的计算结果见表4-1表4-4。此外，按清华方法对一台440 t/h 无烟煤、440 t/h 烟煤、480 t/h褐煤CFB锅炉炉膛也进行了同样的传热计算，结果示于表4-10。4.1.5 100%负荷全炉膛传热量计算结果的校核在上节中已经求出水冷壁、双面水冷壁

12、、屏过、屏再四部分受热面所吸收的热量。其和应等于锅炉热平衡计算中在炉内的传热量。以新乡440 t/h锅炉主循环回路作为对象，热平衡炉内传热量Q1： (4-32)其中 (4-33)式中 Iff回料器及冷渣器反回风带入的热量，kJ/kg；主循环回路出口(分离器出口)烟气焓，烟温883 查温焓表，当过量空气系数=1.2时，。 Ifh离开主循环回路(分离器出口)的飞灰带走的热焓，kJ/kg kJ/kg (4-34)式中 fh 飞灰份额，%；Aar 燃料中灰份，%；Ash 加石灰石产生的灰份；Cfh 飞灰可燃物含量，%；(C) fh飞灰热烩，kJ/kg。将具体数据代入式(3-34)后得：kJ/kg炉膛传

13、热计算中炉内四种受热面总的吸热量为：129.29MW(水冷壁)+25.59MW(水冷屏)+46.94MW(屏过)+34.15MW(屏再)=235.97MW，该数值与炉内热平衡计算的传热量232.66 MW相差小于1.5%，故可以结束计算。4.1.6 低负荷传热计算一般的，煤粉炉当处于低负荷运行时，相对于正常负荷时，炉膛中的水冷壁受热面显得过大，导致炉内温度水平大大降低，炉膛出口温度也下降。为了维持低负荷时汽温仍保持在额定范围内，在设计锅炉时，除了额定工况的计算外，还必须进行70%、50%负荷的计算，这时一般要大大增加过热器及再热器受热面，以保证低负荷时温压大大降低的情况下仍能达到汽温的要求。但

14、对于循环流化床锅炉，低负荷时，烟气流速减小，烟气携带固体的能力下降，可使理论燃烧温度上升(参照下一节)，从而可以弥补由于在低负荷时相对于正常负荷时过大的水冷壁受热面而造成的烟气过度冷却。同时，也可以降低水冷壁的传热系数，使炉膛出口温度较少变化，从而维持过热汽温达到额定值。低负荷传热计算一般进行75%和50%额定负荷计算。下面讨论几个工况参数的变化情况。(1) 床层温度和炉膛出口温度100%负荷时由于内外物料循环流量较高，炉膛上下乃至于整个主循环回路的温度基本一致。但低负荷时炉内，物料循环流率显著降低，趋向于鼓泡床，故床层温度显著高于炉膛出口温度。这时为了求得床层温度，就得进行分段计算，进行密相

15、区传热计算。而为了求得炉膛出口温度仍可以进行全炉膛计算。(2) 密相区燃烧率为了进行分段计算，就需要知道密相区的燃烧率、上升和下降的物料量和物料温度。经分析，低负荷时燃烧工况向鼓泡床转化，故燃烧率a应大于正常运行时的m。以100 MWe级CFB锅炉为例，正常运行时取m0.47，低负荷时取a0.6。(3) 上升与下降循环物料的温差考虑循环物料量降低，故上升与下降物料的温差也应减小，取为3。(4) 烟气速度u0烟气速度受煤耗量Bj和烟气体积(由于增加，体积增加)和烟气温度Qpj的影响，一般低负荷时烟气速度下降。以100 MWe机组为例，100%负荷时u0=5.68 m/s；75%负荷时u0=3.8

16、1 m/s；50%负荷时u0=3.18 m/s。(5) 上升的循环物料量由于负荷降低，分离器效率降低，故循环物料量也相应比满负荷时要降低。降低多少可以通过校核计算求知。就是说，根据锅炉说明书给出低负荷时的床温cc或根据实际运行时测出的床温来反求循环物料量。至于下降和上升的循环物料量比m也只能通过校核计算求得。从50%负荷实际计算看出密相区燃烧率变化对物料浓度影响不大，而改变下降与上升的物料量比m值则对物料浓度影响很敏感。m减少，则物料浓度Cp减小很多。物料浓度除按上述校核计算求取外，可按式(4-35)计算。 (4-35)Gs可由资料根据烟气速度求取，例如图5-2。假定烟气速度为3.18 m/s

17、，则Gs7，则kg/m3。(6) 分离器分离效率低负荷时由于烟气量减少，则分离器进口烟气速度降低，因而使分离器效率降低，从而导致循环量GLC和物料浓度Cp减少。(7) 烟气辐射层厚度s烟气辐射层厚度sg随着负荷的下降而下降，可参照资料计算，但它对传热影响不是很大。以440 t/h锅炉为例所进行的50%负荷全炉膛计算结果见表4-1表4-4。其中4种受热面总计传热量为55.36+13.02+24.78+18.80=111.96 MW；而根据热平衡计算炉内传热量为119.55 MW，误差为6%。表4-5为相关的440 t/h锅炉50%负荷性能参数计算结果。由于床层温度是可控制量，因此计算中通常假定某

18、个低于满负荷的温度作为计算依据。以此为基础，进行炉膛传热计算，得到炉膛出口烟气温度。为便于计算，在积累了大量经验的基础上，低负荷计算可以根据经验确定床底温度，第五章表5-11给出了经验总结结果，是可以用于设计计算的。表4-1 某440t/h CFB锅炉100% 、50%负荷全炉膛水冷壁传热计算项 目符号单位100%负荷50%负荷烟气速度Vfm/s5.683.18床侧温度TbK11851012受热面内工质温度TfK613613管节距Sjm0.090.09管外径dm0.060.06鳍片厚度m0.0060.006管壁厚1m0.00650.0065物料浓度CPPkg/m31.720.73炉膛总高度Hl

19、tm39.4139.41炉膛下部计算高度Hpgm5.55.5梯形段上直段耐火层高度Hnhm0.450.45局部物料空间浓度Cpkg/m322.469.56颗粒对流理论换热系数W/m2K100100烟气中水蒸汽份额rH2O%0.0640.064烟气中三原子气体份额r%0.150.15烟气侧水冷壁总面积Htm212031203工质侧换热系数fW/m2K1500015000实际设计运行系数*Xiu10.211受热面受热情况N单面1、双面211烟气辐射厚度Sm0.20.1壁面黑度w0.80.8受热面金属导热系数W/m2K40.3940.39受热面壁面污染系数sm2K/W0.00050.0005受热面耐

20、火层厚度am100100涂层水冷壁面积m2360.7360.7常数B1/22/32/32/3Boltzmann常数W/m2K45.67E-085.67E-08烟气对流系数*Cgc55鳍片宽度系数*0.9423250.942325耐火材料系数Aa02.52.5耐火材料系数Ba10.000250.00025续表4-1 某440t/h CFB锅炉100%、50%负荷全炉膛水冷壁传热计算项 目符 号单 位100%负荷50%负荷颗粒对流系数0.236 0.108 颗粒对流理论换热系数W/m2K56.306 19.325 烟气对流换热系数W/m2K16.866 11.237 对流换热系数cW/m2K73.

21、172 30.563 烟气辐射减弱系数k0.08656 0.15146 物料表面平均黑度0.74161 0.53498 固体物料黑度p0.90490 0.80993 烟气黑度g0.01716 0.01503 床层黑度b0.90653 0.81279 系统黑度0.73904 0.67553 受热面管壁温差TwK13.767 9.604 管外壁温度TwK626.767 622.604 辐射换热系数rW/m2K136.432 88.390 换热系数bW/m2K209.604 118.953 鳍片高度hm0.01500 0.01500 折算高度*hm0.01592 0.01592 有效高度*h”m0.

22、01592 0.01592 鳍片厚度系数*v0.04333 0.04333 折算厚度*0.00026 0.00026 参数4.05466 3.11917 鳍片利用系数0.99861 0.99918 鳍片面积比(P)Hfin/Ht0.25370 0.25370 名义床侧换热系数W/m2K189.654 112.252 受热面内外面积比Ht/Hf1.602 1.602 壁面平均温度*Tw_K619.884 617.802 受热面内外温差*TK572.000 399.000 受热面耐火层平均温度*Ta_K905.884 817.302 受热面耐火层导热系数*aW/m2K2.657 2.635 附加热

23、阻as0.03814 0.03845 传热系数KW/m2K180.490 108.977 光管水冷壁受热面吸热量QggMW124.20 52.31 涂层水冷壁传热系数KW/m2K23.157 21.218 涂层水冷壁吸热量QtcMW4.778 3.054 水冷壁受热面总吸热量QMW128.976 55.362 表4-2 某440t/h CFB锅炉100%、 50% 负荷双面水冷壁全炉膛传热计算项 目单 位符 号100%负荷50%负荷烟气速度Vfm/s5.683.18床侧温度TbK11851058受热面内工质温度TfK613613管节距Sjm0.07270.0727管外径dm0.060.06鳍片

24、厚度m0.0060.006管壁厚1m0.00650.0065物料浓度CPPkg/m31.720.73炉膛总高度Hltm39.4139.41双面水冷壁总高度Hssm2727局部物料空间浓度Cpkg/m314.7376.272颗粒对流理论换热系数W/m2K100100烟气中水蒸汽份额rH2O%0.0640.064烟气中三原子气体份额r%0.150.15烟气侧水冷壁总面积Htm2260260工质侧换热系数fW/m2K1500015000实际设计运行系数*Xiu10.211受热面受热情况N单面1、双面222烟气辐射厚度Sm0.80.4壁面黑度w0.50.80.80.8受热面金属导热系数W/m2K40.

25、3940.39受热面壁面污染系数sm2K/W0.00050.0005受热面耐火层厚度am100100涂层水冷壁面积m214.814.8常数B1/22/30.50.5Boltzmann常数W/m2K45.67E-085.67E-08鳍片宽度系数*0.980.98耐火材料系数Aa02.52.5耐火材料系数Ba10.000250.00025烟气对流系数*W/m2K，4555续表4-2 某440t/h CFB锅炉100%、50%负荷双面水冷壁全炉膛传热计算项 目单 位符 号100%负荷50%负荷颗粒对流系数0.162 0.072 颗粒对流理论换热系数W/m2K38.630 12.929 烟气对流换热系

26、数W/m2K16.866 11.237 对流换热系数cW/m2K55.496 24.166 烟气辐射减弱系数k0.04022 0.06867 物料表面平均黑度0.64007 0.43907 固体物料黑度p0.86103 0.75637 烟气黑度g0.03167 0.02710 床层黑度b0.86543 0.76297 系统黑度0.71150 0.64075 受热面管壁温差TwK13.878 10.797 管外壁温度TwK626.878 623.797 辐射换热系数rW/m2K131.365 92.170 换热系数bW/m2K186.861 116.336 鳍片高度hm0.00635 0.00635 折算高度*hm0.00645 0.00645 有效高度*h”m0.00456 0.00456 鳍片厚度系数*v0.20472 0.20472 折算厚度*0.00123 0.00