本帖最后由 jayzone987 于 2025-9-28 13:01 编辑
1.前言 伴随着我国工业的迅猛发展，环境问题日益凸显。如何有效利用能源，尽可能的节能减排是我国经济发展需要解决的首要问题。据统计，工业生产中消耗的大量燃料，其热效率却很低，一般在30%左右，而被高温烟气等带走的热量缺高达40%~60%[1]。余热锅炉可利用高温烟气携带的热量产生蒸汽，被广泛应用于合成氨、乙烯、硫酸等行业。而余热锅炉的工艺参数设计的好坏，直接决定着是否能设计出经济性好、能耗低、传热效率高的余热锅炉。本文以某项目余热锅炉设计为例，从物性模拟、热力计算及水动力计算三个方面全面地阐述余热锅炉的工艺参数设计。 2.物性模拟 在余热锅炉的热力计算中，烟气的物性参数（如密度、定压比热、动力粘度、导热系数等）是必不可少的，因此准确的物性模拟就至关重要了。 一般情况下，烟气是各种气体的混合物，其成分主要是氮气、一氧化碳、二氧化碳、氧气、水蒸气等。可通过基于状态方程和活度系数方法的热力学模型，求解烟气的物性参数。工程上多用Aspen Plus进行物性模拟，其数据库收集的气液平衡、液液平衡数据十分完备。在硫磺回收装置中，工艺气还含有S的同素异形体，在余热锅炉内会发生相互转化，产生一定的热效应，这时可通过VMGSim软件进行物性模拟，计算出准确的物性参数。图1、图2分别为在Aspen Plus和VMGSim建立的热力学模型。 | | 图1 Aspen Plus热力学模型 | 图2 VMGSim热力学模型 |
热力学模型建立后，需要选择合适的物性方法进行模拟计算。常用的主要有基于状态方程的物性方法（如PENG-ROB、RK-SOAVE、PR-BM等）和基于活度系数的物性方法（如NRTL、UNIFAC、WILSON等）。极性的物系采用活度系数法进行模拟计算，其中又可按照其是否为电解质选择不同方法进行模拟；非极性或极性非常弱的物系则采用状态方程法进行模拟计算。 根据上述物性计算方法，对大气压力（101.325kPa）下的标准烟气（其组成如表1所示）物性进行计算，将计算结果与文献[2]中值进行比较，绘制曲线如图3、图4、图5、图6所示。 表1 大气压力（101.325kPa）下标准烟气组成
通过以上各图中计算值与文献值比较发现，在大气压力下标准烟气的密度、导热系数、动力粘度及定压比热曲线反映出的趋势均一样。计算值与文献值的密度曲线、动力粘度及定压比热曲线吻合较好，导热系数及比热曲线虽略有误差，但误差均在可接受范围内。同时由于标准烟气的物性组成固定，若考虑烟气成分的变化采用文献中值进行计算，则需要进行修正[3]，过程十分繁琐复杂。因此，本文以Aspen Plus进行烟气的物性模拟，可以较快的得出准确可信的结果。 下面以某项目的变换气余热锅炉为例，进行物性模拟，变换气组成及工艺条件如表2、表3所示，要求的副产压力为4.4MPa(A)的饱和蒸汽。 表2 变换气组成 表3 变换气工艺条件 采用Aspen Plus建立热力学模型，物性方法选用PENG-ROB，计算得到变换气的物性参数如表4所示。 表4 变换气物性参数 | | | | | | | 1 | 450 | -7704.274 | 1.753039 | 3.23E-02 | 1.02E-01 | | 2 | 442 | -7718.285 | 1.750272 | 3.20E-02 | 1.01E-01 | | 3 | 434 | -7732.274 | 1.747537 | 3.17E-02 | 1.00E-01 | | 4 | 426 | -7746.241 | 1.744836 | 3.14E-02 | 9.89E-02 | | 5 | 418 | -7760.186 | 1.742172 | 3.11E-02 | 9.78E-02 | | 6 | 410 | -7774.11 | 1.739547 | 3.09E-02 | 9.67E-02 | | 7 | 402 | -7788.013 | 1.736963 | 3.06E-02 | 9.56E-02 | | 8 | 394 | -7801.896 | 1.734422 | 3.03E-02 | 9.44E-02 | | 9 | 386 | -7815.758 | 1.731928 | 3.00E-02 | 9.33E-02 | | 10 | 378 | -7829.601 | 1.729483 | 2.97E-02 | 9.22E-02 | | 11 | 370 | -7843.424 | 1.72709 | 2.94E-02 | 9.11E-02 | | 12 | 362 | -7857.227 | 1.724752 | 2.91E-02 | 9.00E-02 | | 13 | 354 | -7871.013 | 1.722472 | 2.88E-02 | 8.88E-02 | | 14 | 346 | -7884.78 | 1.720255 | 2.85E-02 | 8.77E-02 | | 15 | 338 | -7898.53 | 1.718103 | 2.82E-02 | 8.66E-02 | | 16 | 330 | -7912.262 | 1.716022 | 2.79E-02 | 8.55E-02 | | 17 | 322 | -7925.978 | 1.714015 | 2.76E-02 | 8.43E-02 | | 18 | 314 | -7939.679 | 1.712086 | 2.73E-02 | 8.32E-02 | | 19 | 306 | -7953.364 | 1.710242 | 2.70E-02 | 8.21E-02 | | 20 | 298 | -7967.034 | 1.708487 | 2.67E-02 | 8.09E-02 | | 21 | 290 | -7980.691 | 1.706828 | 2.64E-02 | 7.98E-02 | | 22 | 282 | -7994.334 | 1.70527 | 2.61E-02 | 7.87E-02 | | 23 | 274 | -8007.966 | 1.703821 | 2.58E-02 | 7.75E-02 | | 24 | 266 | -8021.586 | 1.702489 | 2.55E-02 | 7.64E-02 | | 25 | 258 | -8035.196 | 1.701281 | 2.52E-02 | 7.52E-02 | | | | | | |
3.热力计算 一般的管壳式余热锅炉与管壳式换热器的主要区别为烟气温度较高，其入口处及出口处可能设有热防护结构。但余热锅炉的热力计算（传热计算和压降计算）与管壳式换热器类似，其中传热计算尤为复杂。余热锅炉的流体膜传热系数可按式1（传热关联式）进行计算。 式中， Nu——努塞尔数； B——经验系数； Re——雷诺数； Pr——普兰特常数； a、b——幂指数； file:///C:/Users/chenj/AppData/Local/Temp/msohtmlclip1/01/clip_image016.gif——壁温校正系数。 使用上式计算管程和壳程的流体膜传热系数时，可分别按管程和壳程中定性温度下的流体物性参数进行计算。 通过计算得到的管程、壳程的膜传热系数，则可通过下式计算得到总传热系数。 式中， K——总传热系数，W/(m2·K)； Ao——管外表面积，m2； Ai——管内表面积，m2； ho——管外流体管内膜传热系数，W/(m2·K)； hi——管内流体管内膜传热系数，W/(m2·K)； ro——管外污垢热阻，m2·K / W； ri——管内污垢热阻，m2·K / W； rp——金属热阻，m2·K / W。 最后可根据式3计算得到需要的换热面积。 式中， Q——热负荷，W； A——换热面积，m2； EMTD——有效平均温差，℃。 根据以上各式进行传热计算时，由于流体温度是沿换热管长度方向不断变化的，而流体的物性参数又是温度的函数。因此，如何确定定性温度下的流体物性参数则决定着传热计算得准确度。目前，普遍采用分段计算得方法进行计算，以解决这一问题。用合理的分段方式确定的定性温度下流体物性参数，可使计算结果更准确更接近实际。 工程上常用HTRI、HTFS软件进行传热计算，以代替繁琐的手工计算。这些软件经过几十年的不断完善，其计算结果较为准确，更接近实际工况。本文则利用HTRI软件，进行某项目余热锅炉的热力计算。 下面对某项目余热锅炉进行热力计算，根据本设备特点，余热锅炉的型式按TEMA[4]规则选择BXU，按正常工况的设计裕量30%左右进行设计，并校核最大工况，筛选具有代表性的设计方案如下表所示。 表5 设计方案 | | | | | | | | 换热管外径/mm | 25 | 25 | 32 | 32 | 38 | 38 | | 换热管壁厚/mm | 2.5 | 3 | 2.5 | 3 | 3 | 3.5 | | 换热管数量/个 | 1250 | 1400 | 740 | 800 | 550 | 550 | | 换热管有效长度/mm | 3800 | 3500 | 5400 | 5150 | 6600 | 6650 | | 换热面积/m2 | 866 | 910 | 894 | 931 | 947 | 953 | | 管束质量/kg | 22601 | 26890 | 23264 | 27697 | 27531 | 30908 | | 长径比 | 2.0 | 1.8 | 2.8 | 2.6 | 3.5 | 3.5 | | 正常工况 | Q/MW | 20.9 | | EMTD/℃ | 68.8 | | K/W/(m2·K) | 463.0 | 438.8 | 447.5 | 430.0 | 421.8 | 420.1 | | 设计裕量/% | 31.9 | 31.3 | 31.7 | 31.7 | 31.3 | 31.66 | | 压降/kPa | 17.4 | 17.2 | 14.7 | 14.8 | 13.5 | 15.21 | | Q/MW | 30.6 | | EMTD/℃ | 75.6 | | 设计裕量/% | 10.5 | 9.9 | 10.8 | 10.7 | 10.8 | 10.8 | | K/W/(m2·K) | 515.8 | 488.4 | 500.8 | 480.8 | 473.0 | 469.7 | | | | | | | |
比较上述方案，得出如下结论： (1) 方案1采用的Φ25×2.5mm的换热管，其换热面积和管束质量均最小，满足正常工况下的各个工艺要求，但在最大工况下，其压降不满足允许压降的要求； (2) 方案2相比方案1采用了相同外径的换热管，增加了壁厚。理论上在相同管数和有效管长的情况下，增加壁厚可增大流速，强化传热，可降低需要的换热面积，但设计过程中仍要受到允许压降等其他工艺条件的限制。考虑这些因素设计出的方案2，虽然在正常工况下的设计裕量、压降均与方案1接近，但其经济性较方案1差，同时仍然不满足最大工况的允许压降要求； (3) 方案3采用的Φ32×2.5mm的换热管，其换热面积和管束质量也相对较小，长径比较优。在正常工况下的压降为允许压降的73.5%，最大工况下的压降也满足允许压降的要求； (4) 虽然方案4可以满足正常工况和最大工况下的工艺条件要求，但较之方案3，其经济性较差； (5) 方案5和方案6均采用外径38mm的换热管，理论上在相同管数和有效管长的情况下，增加管径可增加换热面积，同时在相同壁厚下，增大管径可降低压降。受到经济性的制约，设计出的方案5方案6虽然可满足各个工况工艺条件的要求，但是其并不是最优方案。 综上所述，方案3为最优设计方案。因此，确定本项目余热锅炉的换热管布管方案为Φ32×2.5mm，有效管长为5400mm，其换热管间距、布管限定圆与筒体内壁距离等具体的布管方式，需考虑呼吸空位及热负荷等因素，因不涉及热力计算，故在此不再赘述。 4.水动力计算 水循环分为强制循环和自然循环，余热锅炉一般采用自然循环，而水动力计算则是为了保证余热锅炉水循环可靠性，典型的余热锅炉水循环模型如图7所示。
图7 典型的余热锅炉水循环模型 水动力计算的本质是求解一个结构特性和受热情况给定的流体管网的流量分配问题。各支路中的流量与其中的流动阻力、吸热量和重位压差密切相关。阻力、热量、流量相互影响、相互联系。 余热锅炉作为锅炉的一种，其水动力计算亦可用锅炉的水动力计算方法。锅炉的水动力计算方法主要由图解法和平差流量调节算法，其中图解法主要依据“串联管阻流量相同、压差相加”及“并联管组流量相加、压差相同”的原理合成各循环回路的水动力特性曲线，确定总体水动力特性，求解总压差及总循环流量，并借此确定各回路的工作点及各回路的压差、流量。 平差流量调节算法作为一种新算法，主要依据流体力学原理和锅炉水循环计算基本原则，不受系统结构限制，可实现对任意结构锅炉水循环特性的计算。 本文以基于平差流量调节算法的软件BESS进行余热锅炉的水动力计算。由热力计算结果，得到余热锅炉副产蒸汽在正常工况下蒸汽产量为31.6t/h，最大工况下的蒸汽产量为41.2t/h，从文献[5]初选循环倍率为20倍进行试算。设计汽包直径为2000mm，直段长度为9000mm，上升管为DN300mm，下降管为DN250mm，以此建立水动力计算模型如图8所示。
图8 水动力计算模型 水动力计算时需初步确定汽包布置高度，经过反复试算后，确定保证循环倍率最终高度，计算结果如下表所示。 表6 水动力计算结果
| | | | | | | | 正常工况 | 上升管组 | 1.054 | 258.9 | 2.9 | 29.3 | -54.7 | | 余热锅炉 | 1.47 | 258.9 | 0.1 | 31.2 | -10.3 | | 下降管组 | -9.086 | 258.9 | 2.2 |
| 65.1 | | 上升管组 | 1.054 | 259.1 | 3.376 | 20.1 | -55.2 | | 余热锅炉 | 1.47 | 259.1 | 0.12 | 20.9 | -10.1 | | | | |
| |
注：总阻力是正值时，为循环动力值；总阻力是负值时，为循环阻力值。 由以上计算结果，确定了整个水循环的布置形式如图9所示，该布置形式下各个工况的循环倍率均在20倍以上，满足水循环可靠性要求。
图9 余热锅炉水循环布置图 5.结论 (1) 本文采用AspenPlus软件对标准烟气的物性进行模并与文献值比较，验证了此种方法的可靠性及准确性，可采用此种方法进行多组分烟气或工艺气的物性模拟； (2) 本文以某项目的余热锅炉热力计算为例，设计了多种方案，并分析得到了最优方案； (3) 本文以基于平差流量调节算法的BESS软件进行余热锅炉的水动力计算，得到了满足循环倍率要求的水循环的锅炉布置形式。 综上所述，本文以某项目的余热锅炉工艺参数设计为例，按照物性模拟、热力计算及水动力计算分步计算的方法，将余热锅炉的工艺参数设计进行了全面阐述，介绍了一种余热锅炉的工艺参数设计方法可供借鉴。
| 
发表于 2025-9-28 11:27:39
显示全部楼层
IP卡
狗仔卡
楼主
提升卡
置顶卡
喧嚣卡
变色卡
千斤顶
显身卡