2. 辽宁石油化工大学 石油天然气工程学院, 辽宁 抚顺 113001
2. College of Petroleum Engineering, Liaoning Shihua University, Fushun 113001, China
目前, 钢铁企业烧结工序是最主要的二氧化硫 (SO2) 排放工序, 同时也是仅次于炼铁工序的第二大耗能工序.烧结工序本质上是将焦炭、氧化钙等粉状物料和铁矿石细粒料在不完全溶化的条件下烧结成块的过程, 在烧结过程中伴随产生大量烟气.现阶段, 国内大中型钢铁企业将大量烧结烟气统一回收进行脱硫处理, 给烟气脱硫装置带来了巨大负荷[1].烧结烟气中SO2主要来源于烧结混合料中的硫化物 (FeS2, CuFeS2, Cu2S等) 和硫酸盐 (BaSO4, CaSO4, MgSO4等), 在烧结过程中90%以上的硫化物和有机硫被氧化或燃烧生成硫的气态化合物 (主要为SO2) 而释放[2-3].国内外烧结烟气中SO2排放的研究主要关注于SO2迁移和富集机理[4-6]、烟气脱硫方法与工艺[7-9]、烟气中SO2循环模拟[10-12]、烟气脱硫吸收机理[13]等问题的研究, 而鲜有文献针对烧结机水平带长方向烟气中SO2产生机理、浓度分布特点及其影响因素进行研究与分析.
基于此, 本文拟以烧结过程中SO2产生和吸附机理为基础, 针对烧结过程中不同位置化学反应特点, 分别建立了不同区域SO2浓度数学计算模型; 结合国内某大型烧结机运行参数, 对该模型计算结果进行了误差分析与数据拟合; 分析与讨论了料层高度、台车运行速率、混合料配碳量和含水量对烧结烟气中SO2浓度的影响及其变化规律.以期为国内钢铁行业实现“节能减排”目标提供一定的理论依据.
1 烧结烟气SO2浓度分布计算模型 1.1 基本假设烧结过程SO2的产生释放过程是个复杂的化学反应和传质过程, 为了建立烧结烟气SO2浓度计算的数学模型, 需要对产生和传递过程做出相应的简化, 因此本模型将烧结机作为固定床的传质进行解析, 并作如下假设来简化模型:
1) 烧结矿在台车内铺料均匀, 在各个方向属性相同;
2) 忽略纵向、横向之间的传质;
3) 气固之间的传质瞬间达到平衡;
4) 由于烧结矿带的吸附能力远远小于混合料带的吸附能力, 假设其不具有吸附能力;
5) 由于烧结混合料颗粒较细, 其内部扩散远远大于外部扩散, 忽略其内部扩散过程, 所以把在燃烧带和预热带SO2的产生释放模型简化成内释放源模型, 在烧结矿带不存在SO2扩散析出过程, 其产生量以内释放源的形式在分层中平均化处理.
1.2 数学模型的建立 1.2.1 计算模型基本方程式在本计算模型中, 需要对烧结机竖向方向构成的平面内各带的SO2气体的释放—吸收—解吸—扩散排放行为分别逐层进行物质平衡的分析计算[14], 在根据烧结过程SO2的产生机理的基础上做出合理模型假设,已将各带的SO2气体简化为内释放源q, 而吸收过程在本模型的计算过程中需要根据传质理论进行分析, 推导出其传质的控制方程.
基于上文的假设条件, 在层内取单元体积dV=dx·dy·1= dx·dy(宽度为1), 在dt的时间间隔里, 如果该单元体积中存在烧结混合料对SO2的吸收作用, 可以把单元体积看作处于静止状态的不稳定传质来分析, 通过气体和固体颗粒进行SO2的质量平衡计算, 根据模型假设简化可以得到本计算模型中烧结混合料对SO2的吸收作用的3个方程式.
根据单位体积内的质量平衡, 即可得到描述床层吸附剂颗粒内外吸附质浓度变化的连续性方程:
(1) |
传质速率方程表示为
(2) |
式中:ε为孔隙率; c为气体中SO2的质量浓度, mg/m3; ρs为烧结矿密度,kg/m3; t为时间, s; c*为达到平衡时气体中SO2的质量浓度, mg/m3; q为每吨烧结矿的点火热耗, kJ/t; K0为固体对气体中SO2的吸附速度, m3/(mg·s); u为气体的流速, m/s; z为烧结料层高度, mm.
当气体中组分的体积分数较低时 (低于3%), 可以假设吸附等温线为线性关系, 进一步由吸附平衡浓度和平衡压力之间的关系可知吸附平衡质量浓度与吸附量之间也可以表达为线性关系.即
(3) |
式中:kf为吸附量与吸附平衡浓度之间的比例系数, mg/m3; s为固体吸附量, mg/g.
1.2.2 计算基本方程式的建立基于上述假设得到的内释放源和推导出的烧结混合料对SO2的吸收作用的方程式, 结合整个烧结料层与气相之间发生的物理化学变化, 在结合烧结过程SO2的产生机理和烧结料层对SO2的吸收—解吸机理的基础上, 将各带的SO2气体的释放—解吸—扩散排放行为的具体特性以具体求解方程式计算烧结料层的吸收量与气体浓度的变化表示出来.由于目前烧结混合料吸收SO2的研究很少, 吸附速率常数未有定量研究结果[14].通常以烧结料层温度的变化以及烧结过程中所发生的物理化学反应变化为依据, 将正在烧结的料层自上而下划分为烧结矿带、燃烧带、干燥预热带和过湿带.
1) 烧结矿带.焦粉与铁矿石充分燃烧后, 高温烧结矿表面与烧结机上方冷空气接触换热, 形成烧结矿带.在该带中, 燃料燃烧已经结束, 高温熔融物凝固形成多孔的烧结矿饼, 料层中熔融物 (液相) 冷凝, 矿物析晶.因此, 本模型认为在该带中烧结饼不释放SO2, 对SO2也无吸附作用, 即在该带中不存在SO2气体的释放—解吸—扩散排放行为.
2) 燃烧带:铁矿石与焦粉充分接触并发生燃烧反应, 形成燃烧带.在烧结矿带被预热的空气进入燃烧带, 与焦粉接触发生燃烧反应, 放出大量的热 (温度可达1 300~1 500 ℃), 产生高温, 在此条件下, 伴随各种物理化学反应, 低熔点物质继续熔化, 形成一定数量的液相; 部分氧化物分解、还原、氧化, 硫化物、硫酸盐和碳酸盐分解等.在化学反应过程中, 有机硫和碳酸盐中S元素经反应生成SO2气体释放.本模型认为燃烧带产生并释放SO2, 不存在吸收作用, 对于SO2的释放量根据各种硫份的含量及一般转化率采用平均内释放源的方式进行计算.根据质量平衡,烧结料层的吸收量与气体浓度变化的具体求解方程式分别由式 (4) 和式 (5) 表示.
烧结料层的吸收量求解方程式:
(4) |
气体SO2浓度变化的求解方程式:
(5) |
式中:Qx为标态烟气流量, Nm3/h; Qf为单元体积烧结料层硫份转变为SO2的质量, Nm3/h.
在燃烧带中, 在加热到700 ℃左右时焦粉燃烧, 有机硫燃烧产生SO2, 而硫酸盐中的硫元素则在900 ℃以上才分解出来, 燃烧带由于燃烧放出大量的热, 升温极快, 可以认为有机硫和硫酸盐里面的硫元素同时转变成SO2释放出来.
3) 干燥预热带.高温烧结烟气向烧结机底部流动的过程中作为热源与未燃烧的混合料接触, 导致混合料表面温度随之升高, 形成干燥预热带.经过燃烧带产生的高温烟气与混合料进行换热, 在烟气的加热下, 料层中的游离水迅速蒸发, 温度迅速升高, 各种硫化物在被加热的过程中发生氧化反应产生SO2, 同时在此过程中, 由于料层中存在着水分、碱性物质, 对SO2存在着吸附作用.因此在该带中既有SO2的释放又有吸收.本模型则同时考虑SO2的释放和吸收, SO2的释放量是根据硫化物的含量及一般转化率采用平均内释放源的方式进行计算, SO2的吸收量是根据固定床的吸附机理进行计算.根据质量守恒原理, 烧结料层的吸收量与气体浓度的变化的具体求解方程式分别由式 (6) 和式 (7) 表示.
烧结料层的吸收量求解方程式:
(6) |
气体SO2浓度变化的求解方程式:
(7) |
各种硫化物在不同温度发生化学反应产生SO2, 但其前后温度相差不是很大, 而且该带中发生反应的区域比较密集, 本模型假设硫化物产生SO2在某一温度同时释放来计算得到单元体积产生的SO2质量.
4) 过湿带.上层高温烟气中水蒸气进入下层冷料发生冷凝并转移至烧结料中, 导致烧结料的含水量超过混合料原始水分, 形成过湿带.因为过湿带含水量较高, 所以对SO2具有极强的吸收作用, 且没有SO2的产生和释放, SO2的吸收量则根据固定床的吸附机理进行计算.根据质量平衡,烧结料层的吸收量与气体浓度变化的具体求解方程式分别由式 (6) 和式 (8) 表示.
气体SO2浓度变化的求解方程式:
(8) |
在此带中, 由于水分的冷凝导致烧结料层对SO2的吸收作用增强, 本模型对吸收速率常数K0做出了相应的修正, 使其与该带的吸收强度相符.
1.3 计算结果及分布规律以国内某360 m2烧结机运行参数为例 (见表 1), 运用烧结烟气SO2质量浓度分布计算模型, 计算得到该烧结机烧结烟气SO2质量浓度分布规律, 并与现场实测数据进行对比分析, 如图 1所示.
由图 1分析可知, 烧结烟气SO2质量浓度分布模型计算结果与实验测试的结果之间相对误差小于5%.因此, 认为该模型计算结果真实可靠, 可应用于烧结烟气SO2质量浓度分布规律及其影响因素等问题研究.
为研究烧结烟气SO2质量浓度分布规律, 现对模型计算结果进行非线性拟合 (图 1中实线).拟合结果表明:烧结烟气SO2浓度分布随烧结机带长的增加呈Cubic函数分布 (即呈先逐渐增加达到峰值点, 而后逐渐下降的趋势), 其拟合决定系数R2约为0.948, 拟合精度较高.
根据《钢铁烧结、球团工业大气污染物排放标准 (GB 28662—2012)》中要求, 当烧结烟气中SO2质量浓度低于200 mg/m3时, 烟气可直接排放.因此, 该烧结机最前端约至带长方向1/3处烧结烟气满足直接排放要求, 可减轻烧结烟气脱硫设备工作负荷.
2 烧结烟气SO2浓度影响因素分析以某360 m2烧结机运行参数为例, 计算了不同操作参数 (料层高度和台车运行速率) 以及原料参数 (混合料配碳量和含水量) 情况下,烧结烟气SO2浓度沿台车运行方向上的变化曲线.
2.1 料层高度的影响图 2为600, 650, 700和750 mm料层高度下烧结烟气SO2浓度沿台车运行方向上的变化曲线.由图 2可以看出, 烟气SO2浓度随着烧结机带长的增加而呈Cubic函数规律, 且料层高度越低, 烟气SO2浓度峰值越小, 位置越靠前.当料层高度为600 mm时, 烟气SO2质量浓度峰值 (约为661×10-4 mg/m3) 出现在17#风箱位置; 当料层高度为750 mm时, 烟气SO2浓度峰值 (约为842×10-4 mg/m3) 出现在22#风箱位置.这是因为在烧结过程中, 燃烧带所释放的烟气中的部分SO2被下方的烧结混合料所吸附, 但随着燃烧的继续, 下方的烧结混合料最终燃烧殆尽, 造成了烟气SO2浓度沿带长方向先增加而后降低的趋势.此外, 在料层高度较高的情况下, 其底部吸收的SO2相对较多, 当烧结燃烧带到达烧结底层时, 底层所吸收的硫分再次释放, 所以产生的SO2也越多.在烧结机尾部, 烟气SO2浓度随着料层高度的降低而迅速下降, 当料层高度小于650 mm时, 烧结机尾部SO2质量浓度低至0, 这是由于烧结机尾部烧结混合料中的硫分已经反应完全所致.
图 3为台车运行速率分别为2.1, 2.2, 2.3和2.4 m/min时烧结烟气SO2质量浓度沿台车运行方向上的变化曲线.
由图 3可以看出, 台车运行速率越小, 烟气SO2浓度峰值点距离烧结机尾端越远, 且台车尾端处的烟气SO2浓度越低.其中, 当台车运行速率分别为2.1和2.4 m/min时, 烟气SO2浓度峰值分别出现在18#和20#风箱位置.这是因为在台车运行长度、料层高度、混合料配碳量和含水量一定的情况下, 台车运行速率越小, 到达同一位置风箱所需的时间越长, 烧结时间也就越长, 其烧结混合料厚度越小, 吸附烟气SO2作用就越小, 从而烟气SO2质量浓度峰值出现越早, 尤其是当台车速率小于2.3 m/min时, 会造成过烧的情况, 影响烧结矿质量.
2.3 混合料配碳量的影响图 4为混合料配碳量分别为2.6%, 2.8%, 3.0%和3.2%时烧结烟气SO2质量浓度沿台车运行方向上的变化曲线.
由图 4可以看出, 烟气SO2浓度受混合料配碳量的影响明显, 混合料配碳量越大, 烟气SO2浓度峰值越大, 越提前出现.其中, 当混合料配碳量为3.2%时, 烟气SO2浓度峰值 (约为862×10-6) 出现在19#风箱位置.而当烟气SO2浓度峰值点过后, 烟气SO2浓度下降速率随着混合料配碳量的增大而增大.这是因为混合料配碳量越高, 烧结速率越快, 在同一风箱位置, 烧结混合料吸附SO2作用越小, 烟气SO2浓度就越高, 同时由于烧结速率越快, 到达烧结机终点的时间就越短, 烟气中SO2浓度峰值点也就越先出现.
2.4 混合料含水量的影响图 5为混合料含水量分别为5%, 6%, 7%和8%时烧结烟气SO2浓度沿台车运行方向上的变化曲线.
由图 5可以看出, 烧结烟气SO2浓度受混合料含水量的影响不明显, 其峰值均出现在20#风箱位置.在同一风箱所对应的位置, 烟气SO2浓度随着混合料含水量的降低而增高, 其中, 当混合料含水量分别为5%和8%时, 烧结烟气SO2浓度峰值分别约为785×10-6和723×10-6.这是因为混合料含水量越高, 单位烧结混合料的吸附作用就越大, 其烟气中SO2浓度就越高.由于吸附的SO2主要富集在料层底部, 所以吸附作用越大, 则料层底部燃烧后释放的SO2就比较大, 且峰值就会相对更高[15].相对较小的含水量, 其消耗的热量也相对较少, 烧结速度就相对较快, 烟气SO2就释放较快, 尾部就先出现低SO2浓度, 末端风箱位置的SO2浓度就越低.
综上所述, 烟气SO2浓度沿着烧结机带长方向上的分布呈先增大后减少的趋势, 其中, 料层高度、台车运行速率和混合料含水量越低, 及混合料配碳量越高, 烟气SO2浓度峰值出现的位置也就越靠前.
3 结论1) 以某360 m2烧结机为例, 验证了烧结烟气SO2浓度分布计算模型, 相对误差小于5%, 拟合得到烧结机带长方向上烟气SO2浓度分布呈Cubic函数分布, 拟合决定系数R2约为0.948.
2) 料层高度、台车运行速率以及混合料配碳量的变化对烧结烟气SO2浓度分布影响较明显, 而混合料含水量的影响较小, 烟气SO2峰值出现的位置随着料层高度、台车运行速率和混合料含水量的降低以及混合料配碳量的升高而提前.
3) 当混合料含水率不变 (7%) 时, 保持料层高度在650~700 mm之间, 台车运行速率约2.3 m/min, 混合料配碳量约为3%, 有助于保证混合料的烧结质量, 同时, 烧结机最前端至带长方向1/3处烧结烟气SO2浓度均小于200 mg/m3.
[1] | Long H, Chen X, Chun T, et al. Sulfur balance calculation of new desulfurization technology in the iron ore sintering process[J]. Metallurgical Research & Technology, 2016, 113(1): 107. |
[2] | Liu Z J, Zhang J L, Yang T J. Research and development of sintering flue gas desulphurization technology[J]. China Metallurgy, 2009, 2: 1–5, 9. |
[3] | Liu W Q. Development of sintering flue gas desulfurization technology in iron and steel industry[J]. China Environmental Protection Industry, 2009, 5: 24–27. |
[4] | Zhang H, Rao M, Fan Z, et al. Effects of circulated flue gas components on iron ore sintering[J]. ISIJ International, 2012, 52(12): 2139–2144. DOI:10.2355/isijinternational.52.2139 |
[5] | Li G, Liu C, Rao M, et al. Behavior of SO2 in the process of flue gas circulation sintering (FGCS) for iron ores[J]. ISIJ International, 2014, 54(1): 37–42. DOI:10.2355/isijinternational.54.37 |
[6] |
董辉, 李磊, 蔡九菊, 等.
烧结余热回收竖罐内料层传热过程数值计算[J]. 东北大学学报 (自然科学版), 2012, 33(9): 1299–1302.
( Dong Hui, Li Lei, Cai Jiu-ju, et al. Numerical simulation of heat exchange in vertical tank of waste heat recovery[J]. Journal of Northeastern University (Natural Science), 2012, 33(9): 1299–1302. ) |
[7] | Wang S J, Zhu P, Zhang G, et al. Numerical simulation research of flow field in ammonia-based wet flue gas desulfurization tower[J]. Journal Energy Institute, 2014, 88(3): 284–291. |
[8] | Wang X, Li Y, Zhu T, et al. Simulation of the heterogeneous semi-dry flue gas desulfurization in a pilot CFB riser using the two-fluid model[J]. Chemical Engineering Journal, 2015, 264: 479–486. DOI:10.1016/j.cej.2014.11.038 |
[9] | Ooi T C, Lu L. Formation and mitigation of PCDD/Fs in iron ore sintering[J]. Febs Letters, 2011, 85(3): 291–299. |
[10] | Wang G, Wen Z, Lou G, et al. Mathematical modeling of and parametric studies on flue gas recirculation iron ore sintering[J]. Applied Thermal Engineering, 2016, 102: 648–660. DOI:10.1016/j.applthermaleng.2016.04.018 |
[11] | Wang G, Wen Z, Lou G, et al. Mathematical modeling and combustion characteristic evaluation of a flue gas recirculation iron ore sintering process[J]. International Journal of Heat and Mass Transfer, 2016, 97: 964–974. DOI:10.1016/j.ijheatmasstransfer.2016.02.087 |
[12] | Yamaoka H, Kawaguchi T. Development of a 3-D sinter process mathematical simulation model[J]. ISIJ International, 2005, 45: 522–531. DOI:10.2355/isijinternational.45.522 |
[13] | Wang J, Zhong W. Simultaneous desulfurization and denitrification of sintering flue gas via composite absorbent[J]. Chinese Journal of Chemical Engineering, 2016, 24(8): 1104–1111. DOI:10.1016/j.cjche.2016.04.005 |
[14] |
潘建. 铁矿烧结烟气减量排放基础理论与工艺研究[D]. 长沙: 中南大学, 2007.
( Pan Jian.Theoretical and process studies of the abatement of flue gas emissions during iron ore sintering[D].Changsha:Central South University, 2007. ) |
[15] | Chen X L, Huang Y S, Gan M, et al. Effect of H2O (g) content in circulating flue gas on iron ore sintering with flue gas recirculation[J]. Journal of Iron & Steel Research International, 2015, 22(12): 1107–1112. |