D0I:10.13374才.i8sn1001663.2013.02.013 第35卷第2期 北京科技大学学报 Vol.35 No.2 2013年2月 Journal of University of Science and Technology Beijing Feb.2013 等温固定床多球团还原模型 国宏伟2)，范德球，张建良12，曹英杰)，孙德鹏) 1)北京科技大学治金与生态工程学院，北京100083 2)北京科技大学钢铁冶金新技术国家重点实验室，北京100083 通信作者，E-mil:ghwustb@gmail.com 摘要针对单个球团未反应核模型在研究多球反应体系中存在的局限性，首先基于单球团未反应核动力学模型，构 建了等温固定床中多球体系下球团矿还原模型：然后通过引入无因次变量，得到方程的解析解，并且采用特征线法进行 求解得到了方程的数值解：最后在中温管式炉中设计了等温固定床球团矿还原试验，实验结果与模型计算值有较好的 吻合. 关键词球团矿：矿石还原：数学模型：等温：固定床 分类号TF325.1 Reduction model of multiple pellets in isothermal fixed beds GUO Hong-wei)FAN De-qiu),ZHANG Jian-liang2),CAO Ying-jie,SUN De-peng 1)School of Metallurgical and Ecological Engineering,University of Science and Technology Beijing,Beijing 100083,China 2)State Key Laboratory of Advanced Metallurgy,University of Science and Technology Beijing,Beijing 100083.China Corresponding author,E-mail:ghwustbagmail.com ABSTRACT For the limitation of the single-pellet unreacted shrinking core model in studying a multi-pellet system. a mathematical model of iron oxide pellet reduction in an isothermal fixed bed was constructed based on the reduction behavior of single pellets.The analytical solution of the model was obtained by introducing dimensionless variables, and the characteristic line method was used for numerical analysis.At last the reduction experiments of the fixed bed consisted of inultiple pellets were performed in a middle temperature tube furnace and the predicted results by the model agreed well with the measured ones by experiment. KEY WORDS iron oxide pellets;ore reduction;mathematical models;isothermal:fixed beds 球团矿是高炉炼铁的主要原料之一，研究球团 (3)阶段模型).模型既不能完全由均相模型来 矿在高炉中的反应热力学和动力学行为，对于高炉 描述，又不能由界面模型描述.在球团矿还原初期， 的高效、低成本治炼具有重要的意义四 反应在球团表面附近进行的速度很快（内扩散阻力 目前对于球团矿在高炉中的动力学行为2-3到研 小)，在球团靠近表面的一定区域内反应类型趋向于 究主要有以下模型. 均相反应，此过程在产物层形成之前的阶段称为第 (1)界面模型.该模型认为铁矿石的还原反应 一阶段，球团表面附近的固体反应物被消耗形成产 是在未反应核心(Fe2O3)与产物(Fe)之间的界面上 物层后，称为第二阶段. 进行. (4)微粒模型8).模型假设固体球团是由很多微 (2)均相模型5-).该模型认为气体反应物在反 小颗粒组成，微粒与微粒之间的空隙构成球团内部 应过程中通过固体反应物的扩散系数很大，且在整 的孔隙结构，反应气体正是通过微粒与微粒之间的 个反应过程中扩散系数将保持不变 空隙扩散到固体球团的内部参与球团的还原反应 收稿日期：2012-0907 基金项目：高等学校博士学科点专项科研基金资助项目(20090006120004)：国家自然科学基金资助项目(51204013)：中央高校青年 人才培养基金资助项目(FRF-TP-12-020A)
DOI ：10．13374／j ．issn1001－053x．2013．02．013
162 北京科技大学学报 第35卷 (5)多级区域反应模型9.该模型能够同时描述 H→M、M→W和W→Fe的多级反应情形（其中H 代表赤铁矿，M代表磁铁矿，W代表浮氏体)，并得 出了反应时间加和性特征 (6)孔隙模型1.该模型认为，在气固反应过 程中，气体反应物通过固体反应物内部孔隙扩散进 入固体反应物内部，同时在扩散过程中，气体反应 物又同时与孔壁处的固体反应物发生反应生成固体 产物.由于生成的固体产物与原本的固体反应物之 间具有不同的物理性质（如密度），这样在反应过程 中就会造成固体反应物孔隙结构的变化 建立模型对实际球团矿还原行为进行模拟，其 主要目的令其过程简单化.上述模型研究中，以未 CO 反应核1-12界面模型，应用最为广泛，但模型多 停滞于单球体系下的研究 图1固定床球团示意图 针对单个球团未反应核模型在研究多球反应 Fig.1 Schematic of the fixed bed with pellets 体系中存在的局限性，笔者首先基于单球团未反应 核动力学模型，构建了等温固定床中多球体系下球 (1-E)AdzPFe2O3 团矿还原模型 MFe203 +- 1 等温固定床球团矿还原模型 =-Reo Adtdz. (3) 1.1模型建立 化简式(3)得 1.1.1模型的假定 af ()床层内反应温度恒定，且反应过程中床层 R'co 0t=31-ea (4) 的空隙率保持不变： (2)床层内气流的速度恒定： 式中：∫为球团的转化率；t为时间，s:a=po/ (3)忽略径向上床内参数的变化，即转换率和 MFe203,其中ppe2O3为Fe2O3密度(kgm-3), 反应气浓度在径向上保持不变 MFezO3为Fe2O3的摩尔质量(kg-mol-). 1.1.2控制方程的建立 式(2)和(4)的偏微分方程为描述本模型基本 (1)CO物料衡算.设在反应时刻t,取如下图1 控制方程.R吃。是球团转化率∫和床内反应气CO 所示的床层内z~z+dz的控制体 的浓度Cco的函数，选取拟稳态下，界面反应为一 对控制体内CO进行物料衡算有 级可逆反应时，并将转化率的表达式∫=1-r3/ 代入（其中r表示未反应核半径，m:ro表示球团半 EuCCoA-EUA aCcodz Cc0+82 径，m),未反应核速率表达式如下： Roo Adz EAdz aCco (1) Ot Rco= 化简式(1)得 3(1-E)(Cco-Cco.)/To 8Cco=-Rco aCco+ (2) 言+发1-小n-+- (5) 式中：e为床层孔隙率：u为相对速度，ms-1；A 式中：D。表示扩散系数，m2.s-1:r表示气相 为固定床的面积，m;0=eu为反应气的空 传质系数，m:s1；k表示表观速率反应常数， 塔速度，ms1；Cc。为床内反应气中C0的浓 s-1:Cc0,e表示C0的化学反应平衡浓度，mol 度，molm3:R。为C0的消耗速度，molm-3s-1. m-3:K9表示化学反应平衡常数 (2)球团内Fe2O3的物料衡算.在反应时间间 此外，式(2)和(4)的初始条件为 隔t~t+dt内对以上所取的控制体内Fe2O3进行 物料衡算，有以下表达式： t=0:f=0,Cco=Cco.e. (6)
第2期 国宏伟等：等温固定床多球团还原模型 163· 边界条件为 对式(12)两边进行积分，并化简得 z=0:Cc0=Cc0,0 (7)） 1-x=品(+) (13) 在给定的初始条件（⑥）和边界条件(7)下，通 过求解偏微分方程(2)和(4)便可得到固定床内沿 着床层高度方向上球团转化率∫和CO浓度的分 x=品-(x+) (14) 布情况.以下从解析和数值两个方面对上述问题进 将式(14)代入到(9)中，且Y为常数，得 行求解. 1.2固定床内操作特性解析 -o-x+n (15) 在1.1中固定床内一维模型是由一组偏微分方 程(2)和(4)描述的，用解析的方法求解这样一组 在任意确定位置，对式(15)积分得 偏微分方程有很大的困难，通常情况下都是使用数 值的方法进行解析的.关于解析解，Moriyama13曾 (16) 通过一定合理的假设和简化得到近似的解析解.以 =9-，x- 下便是近似求解的过程. 对于一般的气固反应有，γ的数量级在103 首先，引入适当的无因次量 ~10-5,于是Y0Xdn<0fd:且随着反应不断 无因次浓度：X=(Cco-Cco.c)/(Cco.0- 地进行，有y0>y0xdn.因此式(16)可以在保 Cco.e). 证计算精度的情况下略去Y0Xdn项，简化为 无因次时间：0=ut/L. 无因次轴向距离：7=L. (17) 固定床内气固反应物的摩尔比：y=e(Cco.0- Cco,e)/i3a(1-e八： 以下可以通过从式(17)出发，得出固定床反应 无因次参数：a=3(1-e)Lk/roo,3=3(1- 器的操作特性 E)LDe/r3uo:6=3(1-e)Lk(1+1/Ke)/rouo. 将以上的无因次量代入到式(2)、(4)、(6)和(7) 1.2.1反应完结带高度ne和反应完全时间0： 中，得到以下的无因次方程： 在固定床中，反应完结带高度指的是在该高度 axax =-X(f)x, 上，颗粒填充床内的固体反应物已经反应完全，且 (8) 在低于此高度方向上床内颗粒物均己反应完全，如 图2所示. 器-A (9) 由固定床内反应特性知，反应完结高度是随着 式中， 反应时间的推移逐渐向固定床内部推进的，因此位 于固定床最底层的球团将最先被完全还原 Af)=1/1/a+(1/B){1-f)1/3-1} +(1/)(1-f)231. 初始条件(6)的无因次形式： 未反应带 0=0:f=0,X=0. (10) 反应带 边界条件(7)的无因次形式： 7=0:X=1. (11) 联立式(8)和(9)，消去(f)得 X+ a 反应完结带 品(x+)-0变形得 =品(x+) 图2不同区域的示意图 (12) Fig.2 Schematic diagram of different zones
.164 北京科技大学学报 第35卷 设c=0,即固定床内最底层球团的反应完结时 观动力学参数(r、D。、k、K曰等)、操作条件气流 间为9。，则在≤9。区间，反应完结的高度为零，由 速度uo以及床层参数e、L、o有关，但均与反应时 式(17)可得 间无关.这就说明随着反应时间的推移，床层内球 1f df 0e= 团的转化率以保持不变的分布向床层内部推移.球 7。万 (18) 团转化率的分布特性只与单个球团宏观动力学模型 当反应时间6>9。以后，反应完结高度前沿 的A(f)和式(24)有关. 由最底层球团向床层内部逐渐推进.由于在”< 下面讨论球团宏观动力学过程不同控速时，反 的反应完结带内，球团的转化率为=1，因此将式 应带内球团转化率的曲线分布 (17)在整个反应完结带内积分可得 (1)气膜内传质控速.气膜扩散控速时，式(24) 简化为 dn,ne=70- df (·(19) of =-f. (25) On 将式(18)代入式(19)，整理得 根据式(25)，可以得气膜传质控速时反应带内部球 7e=y9-y8e=y(0-0e). 团转化率的分布如图3(a)所示. (20) 由图3(a)可知当还原性气体C0的气流速度 根据式(20)可以计算得到全床反应完结(=1)的 o分别取5、10和15ms-1时，反应带的宽度变 时间为 宽.这主要是因为气流速度增大以后，：的绝对值 0:=8e+1/. (21) 减小，使得af/8m的变化趋势变缓，af/0m由初始 (8f/8m)f=1=-a变化到（⑦f/an)=n=0需要更长 由式(21)可以得出全床反应完结时间A只与 的距离，所以反应带的宽度变宽 0。和Y有关，因此只要给定化学反应速率中表达 (2)内扩散控速.内扩散控速时，式(24)简化为 式()和量纲一的固定床内气固反应物摩尔比便 可近似地求出9 af Bf 由式(20)有dne/d9=Y,当0>0。时，反应完 8m(1-f)5-1 (26) 结带高度将会以与反应时间成正比例的速度向床层 根据上式，可以得内扩散传质控速时反应带内部球 内部推进，推进的量纲一的速度为。=Y 团转化率的分布如图3(b)所示 1.2.2反应带内球团转化率分布 (3)化学反应控速.化学反应控速时，式(24)简 反应带指的是反应完结带前沿上方球团转化 化为 率在0<f<1范围之内区域，高度为.反应带内 of =8f(1-f)3 On (27) 的球团的转化率在轴向上的分布对于固定床内操作 有着很重要的意义， 根据上式，可以得化学反应控速时反应带内部球团 式(17)两边分别对”求偏导数，则方程式左边 转化率的分布如图3(c)所示. 变为 由图3(c)可看出，Of/Bn的绝对值的变化随着 ∫的增加先变大，当达到某一拐点后，又出现逐渐 df df 1 df 变小的趋势.同理，当提高还原气气流速度后，反 A(万 dn A(f)dn 应带的宽度也明显变宽.这由图3（©）中气流速度 (22) 提高1倍后，反应时间同为2时，红线代表的反 同理，方程式右边变为 应带宽度明显比黑线所代表的反应宽度要宽可以看 出.其主要原因和气膜传质控制时，气流速度增大 (23) 造成8f/⑦m变化趋势变缓类似. 则求出偏导之后，式(17)变为 (4)混合控速.当宏观动力学由混合控速作用 时，只能通过求解式(24)才能得到反应带内球团转 af=-jλf) O (24) 化率的分布，但式(17)中有非线性项很难得到解析 解，这时可采用四阶R-K法14得到数值解. 式(24)表示固定床内反应带中固体转化率的分布 由图3(d)可以看出，随着反应时间的推移，反 情况.由式子右端可知0f/8n的分布只受到球团转 应带内气团转化率以相同的分布曲线往床层内部推 化率∫和()分布的影响，而()与单个球团宏 必
第2期 国宏伟等：等温固定床多球团还原模型 ·165· 1.0 1.0 0.8 0.8 0.6 △nla=5 0.6 0.4 0.4 △=10」 0.2 △lu=15 0.2 0.0 0.0 (a) 6) 1.0 1.0 0.8 0 0.8 8<02<0 0.6 0.6 0.4 0.4 0.2 0.2 0.0 0.2 0.4 0.6 0.8 1.0 (c) (d) 图3不同控速环节时反应带内球团转化率分布.(a)气膜传质控速：(b)内扩散控速：(c)化学反应控速：(d)混合控速 Fig.3 Distribution of pellet conversion rate in the internal reaction zone under different control conditions:(a)gas film diffusion control condition;(b)internal diffusion control condition;(c)chemical reaction control condition;(d)mixed control condition 1.3模型数值求解 根据以上两条特征线，可以得到在边界条件 模型的控制方程与边界条件已经由式(8)~(11) 上，即=0时，边界n和0的分布. 描述.下面采用特征线的方法进行求解， n=8 刀=6-2△0 对式（⑧）和（⑨）中x和f求全微分得 3△ axae+ dx=604 axdm, (28) 20 721 /23 24 26 afan. 8fa0 2△ 22 25 df=5 (29) On 定义如下两条特征曲线： 10 11 12 14 15 /16 △7 d0 =1, (30) dn 00 01 02 03 04 /06 dn=0. 49 2△0 3△4△5△86△07△0 d (31) 图4特征曲线组成的网格 将以上两条特征曲线代入式(28)和(29)，则双曲型 Fig.4 Grid system composed of characteristic curves 偏微分方程组变成了以下的常微分方程组： 首先，7=0时，X=1,将其代入到式(33)中，得 dx =-X(f)x, (32) dn 影=m (34) af =Af)xT. (33) de 积分式(34)得0和f的关系： 积分式(30)和(31)可以得到其积分式如下： 0=f(1/a-1/3)+3[(1-(1-f)23)】/(23) =0+const, (30a) +3[(1-(1-f)/3)]/6. (35) n=const. (31a) 根据上式可以得到在给定的反应时间内边界 这样，在0和7组成的坐标体系内形成以(30a)和 上球团的转化率f.当球团的转化率为1时，反 (31a)表示的一系列特性线交叉所形成网格体系，如 应时间为0。=[1/a+1/(23)+1/(36)〗/y;当反应 图4所示. 时间0>6。时，转化率保持1不变；当反应时间