固定床煤气化过程的数学模型与模拟计算 1999.pdf

固定床煤气化过程的数学模型与模拟计算 项友谦 张利平 中国市政工程华北设计研究院 天津300074 武汉乾能燃气热力有限公司 武汉430056 摘要 根据煤与氧气 二氧化碳和水蒸汽反应的动力学数据 在传质和气 固相传热的基 础上 建立了固定床煤气化过程的数学模型 探讨模型的简化和有关参数的确定 并对常 用炉型进行了计算机模拟计算与分析 关键词 煤转化 固定床煤气化 数学模型 模拟计算 中图分类号 TU996 TQ541 TQ018 1 引 言 我国煤炭气化工艺中 使用最多的是固定床气 化工艺 如鲁奇加压气化炉 常压发生炉 发生炉型 两段炉 水煤气炉 水煤气型两段炉等 这些气化工 艺已很成熟 设计也比较完善 但操作技术还有待继 续提高 过去操作条件的选择 尤其是优化操作条 件的确定多凭经验判断 煤种及条件变化时 要用试 验的方法进行较长时间摸索 才能确定 试验法无 凝是最现实和可靠的方法 它是利用实际生产装置 改变操作条件 测定其热效率 以此来确定优化的操 作条件 此法比较直观 但试验次数多 人力物力消 耗大 测量仪表要求较高 有时还要影响生产 如果 能从气化过程机理出发 根据 三传一反 的原理建 立数学模型 在计算机上进行模拟计算 再根据计算 结果寻找优化的操作条件 这样 所需数据多数可 从反应器外得到 模拟计算方便 人力 物力消耗少 本文主要介绍计算机模拟所需数学模型的建立 数 值解法以及典型过程模拟计算的结果 2 数学模型及数值解法 2 1 非稳态模型 为了使模型的计算简便 失真又不严重 建模过 程特做如下基本假设 1 气化炉内气体及固体呈活塞运动 无径向 参数差别 2 气体的轴向扩散忽略不计 3 反应器内压力视为常数 4 气化段内只进行气化和燃烧过程 根据守恒定律 净流入速率 净发生速率 累积速率 今取气化炉内高度上x与dx之间的微元体对i 组分进行物料平衡得出 5Ci 5t 5 C i Ux 5x ij Rj C T 1 分别对固相和气相进行能量衡算则得到 5 Cp s Cs Ts 5t 5 Cp s Cs Ts Us 5x F Tg Ts K Rj C T HR 2 5 Cp g Cg Tg 5t 5 Cp g Cg Tg Ug 5x F Ts Tg 1 K Rj C T HR 3 式 1 3 构成的偏微分方程组为非稳态煤气 化过程数学模型 水煤气炉 水煤气两段炉 两段式 地下气化过程以及稳态煤气化过程的开停车过程就 属于此类模型 2 2 稳态煤气化过程模型 由于稳态煤气化过程参数不随时间而变化 上 述方程中的对时间的偏导数可取消 则上述方程变 为常微分方程组 8 煤气与热力 1999年9月 收稿日期 1999 05 14 5 C i Ux 5x ij Rj C T 4 5 T s Us 5x F Tg Ts K Rj C T HR Cps Cs 5 5 T g Ug 5x F Tg Ts 1 K Rj C T HR Cps Cs 6 2 3 数学模型参数的确定 煤气化炉主要由气化与燃烧反应段组成 主要 反应为碳的燃烧 碳与水蒸气和碳与二氧化碳的反 应 固体燃料为一多孔介质物料 其反应主要在孔 内表面进行 因此整个反应可设想由下述步骤组成 1 气体反应剂从主流经气膜扩散到固体外表 面 2 气体反应剂经孔扩散到固体内表面 3 气体反应剂与固体在孔内表面进行反应 4 气体产物由孔内扩散到固体外表面 5 气体产物由固体外表面扩散到气相主流 上述各步的速度往往差别很大 总速度取决于 速度最慢的阶段 即所谓 控制阶段 如果对 1 3 步进行考察 在一级反应时 总反应速度r可用 下式表示 r Keff Cin Keff 1 1 Kg 1 Kr 下面分别讨论本征反应速度常数Kr 孔效率 和气膜传质系数Kg的确定 1 本征反应速度常数Kr的确定 当物料颗粒直径小时 颗粒的比表面积大 温度 低时 反应速度低 这时扩散速度大于反应速度 内 表面的反应速度基本上代表整个反应的速度 这时 反应速度常数为本征反应速度常数 在动力学参数 测定时 物料的重量易于测定 因此常以重量为基准 的反应速度常数Km表示 Kr Km 2 孔效率 的计算 孔效率 表示孔扩散对反应速度的影响程度 它可通过x和x dx之间微元体内的反应速度与扩 散速度的平衡 建立微分方程求解得出 对于圆柱形孔 tanh 式中 梯勤模数 Rk kV De 对于锥形孔 1 1 tanh3 1 3 式中 梯勒模数 Dk 6 kV De 3 Kg的计算 Kg可用下式计算 Kg Sh De F DK 对于球形颗粒 Sh 2 0 0 6Re0 6 Sc0 33 对于填充床 Sh 2 1 1Re0 6 Sc0 33 对于燃烧反应 温度一旦超过燃点温度 反应速 度就会很快 以至扩散成为阻力 因此Keff Kg 燃点温度采用前述测得数据 1 6 气化炉内温度高 传热过程非常复杂 既有辐射 传热 又有对流传热 既有固体内部的导热 又有气 固之间的传热 由于气固之间传热系数的大小对整 个气化炉的热平衡影响不大 主要影响气体和固体 之间的温度差别 因此 气体与固体之间的传热系 数采用下式计算即可 Nu Dk 式中 Nu 2 0 1 8Re0 5 Pr0 33 3 气化过程的模拟计算 煤的气化工艺很多 下面就几种典型气化过程 进行模拟计算 3 1 水煤气型两段炉的模拟计算 水煤气型两段炉与水煤气炉差别不大 均属循 环制气过程 数学模型能否反映实际的生产装置 须与实际生产数据进行比较 对于稳态数学模型可 用正常操作数据来比较 对于非稳态数学模型就比 较困难 由于实际生产过程的动态数据很难测定 炉 出口的气体浓度只能测得某一阶段的平均值 对于 水煤气炉来讲 最重要的参数之一是生成煤气的组 成 因此 我们在操作条件相同的状况下以各阶段的 煤气组成作为主要参数进行比较 表1为水煤气型两段炉上吹阶段的实际生产数 据与模拟计算数据的比较示例 5 通过比较 我们 认为水煤气型两段炉的数学模型基本上能反映实际 气化炉的状况 9 第19卷 第5期 煤气与热力 3 2 发生炉型两段炉的模拟计算 发生炉型两段炉与发生炉的区别在于上面有干 馏段 干馏产品的产率取自煤干馏分析数据 表2 为发生炉型两段炉模拟计算与实测值的比较 7 表1 水煤气型两段炉模拟计算数据 O2COCO2H2CH4N2 鼓风 阶段 计算值0 636 1016 6776 60 实测值7 5414 7277 54 上吹 阶段 计 算 值 气化段34 8410 1055 06 干馏段12 9613 2629 1244 66 混 合31 8310 5351 506 13 实测值31 117 6646 945 808 29 上吹阶段循环时间 16112030 最高反应温度 K 14601446142714071387 水蒸气分解率 66 3665 6954 6663 1761 62 表2 发生炉型两段炉煤气组成的模拟值与实测值的比较 上段下段煤气总管 模拟值实测值模拟值实测值模拟值实测值 CO 25 725 028 825 427 425 4 CH4 8 07 71 82 04 53 5 CO2 5 64 06 26 05 94 7 H2 15 413 18 610 949 952 0 N2 44 549 254 355 149 952 0 O2 0 20 50 30 40 30 4 H2S 0 30 00 2 CnHm 0 20 50 20 20 20 3 煤气热值 M J 8 077 625 225 226 666 30 图1 典型原料在鲁奇炉内气化时的温度分布 3 3 鲁奇炉的模拟计算 鲁奇炉是加压固定床 也称移动床 最成熟的炉 型 其优点是 能用劣质煤 气化强度高 能耗低 生 成煤气压力高 便于后处理和远距离输送 图1为三种典型原料 活性炭 焦炭 石墨 在典 型操作条件 粒径0 02 m 压力2 5 MPa 气化强度 800 kg m2 h 汽氧比 6 0 下计算的炉内温度分 布 1 2 3 4 结 语 1 煤与氧气 二氧化碳和水蒸汽反应的动力 学数据 用于固定床气化过程数学模型的计算 基本 01 煤气与热力 1999年9月 可行 2 开发的数学模型已用于多种固定床炉 与 实测数据比较符合 3 固定床气化过程的数学模拟可以预测气化 过程不同条件下各种参数的变化规律 对于指导实 际生产 科研设计以及培训人员均有意义 4 固定床气化过程数学模拟计算 可用于过 程优化条件的确定比试验法安全 省时 省料 具有 广阔应用的前景 符 号 说 明 C 组分的浓度 kmol m3 Cp 比热 kJ kmol K De 扩散系数 m2 h dK 颗粒直径 m F 比表面积 m2 m3 G 吉布斯自由能 kJ kmol H 反应热 kJ kmol K 反应热分配系数 Keff 有效反应速度常数 1 h Kg 气膜传质速率常数 1 h Km 体积为基准的反应速度常数 m3 kg h Kr 有效反应速度常数 1 h n 反应级数 Nu Pr Re Sc Sh 特征数 R 气体常数 kJ kmol K T 温度 或K U 速度 m h 传热系数 kJ m2 h K 堆密度 kg m3 孔效率 导热系数 kJ m h K ij 化学计量系数 梯勤模数 参 考 文 献 1 Xiang Y ouqian Teoretische Berechnugen der G leichgewicht szusamensetzung Bestimung und Mathematisches Model fuer die K ohle vergasung Engler Bunte Institut der Universitaet Karlsruhe TH 1985 2 Xiang Y ouqian et Mathematical model for moving bed solid fuel gasifi cation under pressure The 4th Japan China symposium on coal c1 chemistry proceding Osaca Japan May 25 28 1993 167 172 3 项友谦 固体燃料加压气化过程的分析与模拟 煤气与热力 1987 1 3 11 4 步学朋 彭万旺 项友谦 煤加压气化过程数学模型的研究 煤化 工 1993 1 6 15 35 5 项友谦 循环制气过程的数学模型 城市煤气 1995 12 14 17 6 项友谦 煤气化过程反应活性及动力学研究 煤气与热力 1999 4 4 10 7 张乙明 陶瓷工业窑炉用煤气发生炉气化过程分析及计算的模 拟 硕士论文 天津 天津大学 1988 MATHEMATICAL MODEL AND SIMULATE CALCULATION FOR FIXED BED COAL G ASIFICATION Xiang Y ouqian North China Municipal Engineering Design Reseach Institute Tianjin 300074 China Zhang Liping Wuhan Qianneng Gas Heat Company LTD Wuhan430056 China ABSTRACT According to the reaction kinetics parameter of the gasification material with O2 CO2 H2O and based on the heat transfer mass transfer gas solid phases heat trans