褐铁矿粉磁化焙烧的实验与模拟研究.pdf
褐铁矿粉磁化焙烧的实验与模拟研究 ① 朱 超, 杨 彪, 龚志军, 李保卫, 武文斐 (内蒙古科技大学 内蒙古自治区白云鄂博矿多金属资源综合利用重点实验室, 内蒙古 包头 014010) 摘 要 针对脱水后的多孔褐铁矿,利用同步热分析仪研究了 CO 在 823~873 K 之间还原褐铁矿的过程,利用传热、传质与化学反 应耦合的动力学模型模拟了微粉磁化焙烧的过程,同时还模拟了孔隙率对焙烧反应的影响。 实验结果表明,细颗粒所需的反应活化 能更低且反应时间明显缩短。 数值计算结果表明,当环境 CO 浓度由 1 mol/ m3增大到 3 mol/ m3时,褐铁矿磁化还原时间减少了 73%; 褐铁矿颗粒粒径越大,气体越难进入到颗粒内部;模拟和实验结果表明,425 μm 粗颗粒磁化焙烧所需时间约为 75 μm 颗粒的 2.5 倍。 关键词 褐铁矿; 磁化焙烧; 数值模拟 中图分类号 O551文献标识码 Adoi10.3969/ j.issn.0253-6099.2015.05.027 文章编号 0253-6099(2015)05-0103-04 Experiment and Numerical Simulation of Magnetic Roasting of Limonite Powders ZHU Chao, YANG Biao, GONG Zhi⁃jun, LI Bao⁃wei, WU Wen⁃fei (State Key Laboratory of Integrated Exploitation of Bayan Obo Multi⁃metal Resources, Inner Mongolia University of Science and Technology, Baotou 014010, Inner Mongolia, China) Abstract Aiming at the porosity of dehydrated limonite, reduction of limonite at 823~873 K with CO was studied by STA. Then, the magnetic roasting process of powder was simulated with a kinetics model that coupled heat and mass transfer with chemical reaction, and the effect of porosity on the roasting reaction was also simulated. The experimental results showed that the activation energy of fine particle reactions was lower and the reaction time was significantly shortened. And the numerical calculation indicated that the reduction time of limonite reduced by 73% when the ambient CO concentration increased from 1 mol/ m3to 3 mol/ m3. Furthermore, limonite of larger particle size brought in more difficulty for the gas into the interior of particles. It can be seen from the testing and simulation results that the magnetic roasting time required for the limonite of particle size of 425 μm is 2.5 times as much as that for limonite with size of 75 μm. Key words limonite; magnetic roasting; numerical simulation 磁化焙烧⁃磁选工艺是有效处理难选弱磁性氧化 铁矿最有效的方法之一,对于难选赤铁矿、褐铁矿及菱 铁矿等矿物,相比其他选矿方法,磁化焙烧可以得到更 好的选矿指标[1]。 目前为止对褐铁矿磁化焙烧的研 究基本只集中于实验研究[2-3]。 文献[4]研究了铁矿 磁化焙烧的热力学和动力学,发现鮞状赤铁矿焙烧效 果远不如赤铁矿和褐铁矿,得出减少颗粒尺寸、增加颗 粒与还原气体 CO 的接触面积可有效提高还原焙烧效 果。 文献[5]对褐铁矿直接还原焙烧的机理进行了研 究,结果表明还原时间和还原温度对还原效果都有着 重要影响,还原时间过长反而导致铁浓度降低,而温度 则正面影响着气体扩散和化学反应。 这些研究都只是 基于实验研究,而磁化焙烧的许多内在信息只能通过 数值模拟得到。 文献[6]通过实验研究了脱水后褐铁 矿的多孔性,并建立了多孔铁矿气固反应的动力学模 型,模拟结果与实验结果吻合得较好,但只讨论了温度 和浓度对磁化焙烧的影响。 本文利用同步热分析仪对不同尺寸褐铁矿颗粒磁 化焙烧动力学进行了分析,利用传热、传质与化学反应 耦合的动力学模型,来描述单个褐铁矿微粒的磁化焙 烧行为。 通过数值模拟,来说明颗粒尺寸、反应气体浓 度及颗粒孔隙率对磁化焙烧的影响。 ①收稿日期 2015-03-19 基金项目 内蒙古应用技术研究与开发资金计划(20130310);内蒙古高校创新团队研究计划(NMGIRT1406);白云鄂博矿多金属资源综合 利用国家重点实验室资助项目(BO-13-006) 作者简介 朱 超(1989-),男,江苏淮安人,硕士研究生,主要研究方向为流化床磁化焙烧。 通讯作者 武文斐(1964-),男,内蒙古包头人,教授,博士,主要研究方向为流化床微波磁化焙烧。 第 35 卷第 5 期 2015 年 10 月 矿矿 冶冶 工工 程程 MINING AND METALLURGICAL ENGINEERING Vol.35 №5 October 2015 1 实 验 1.1 实验材料 实验所用材料为褐铁矿粉,经测试,矿粉主要成分 如表 1 所示。 通过筛分得平均粒径分别为 75 μm 和 425 μm 左右的褐铁矿粉,备用。 实验所用气体分别为 高纯 CO(99.999%)和高纯氮气(99.999%)。 表 1 褐铁矿主要组成(质量分数) / % TFeFeOSiO2Na2OMgOAl2O3 32.801.8031.721.660.404.64 PSF K2O CaO 0.102.060.031.840.72 1.2 实验装置 还原实验装置采用耐驰(NETZSCH) 出品的型 号为 STA 449 F3 Jupiter 的同步热分析仪,分辨率为 0.1 μg。 1.3 实验方法 实验过程中,将实验样品尽可能平铺于坩埚中,样 品经前期加热脱水达到预定温度时开始通入 CO,反应 过程中通过电子天平来记录样品的失重量,反应全程 用氮气保护(流速 50 mL/ min)。 褐铁矿粉磁化焙烧程 度用还原率 β 来表示 β = m0 - m m0 - m 1 (1) 其中 m,m0,m1分别为反应过程中、反应前和反应结束 时的样品质量。 图 1 给出了较小颗粒(75 μm)褐铁矿在不同升温 速率下还原度与反应时间的关系曲线。 图 1 中 80 min 之后曲线突然发生明显下降,当总失重率达到 16.38% 左右时出现一个转折点,它对应 Fe2O3还原成 Fe3O4 的过程[7]。 因此将取转折点前的数据来分析褐铁矿 磁化焙烧过程中样品的失重情况。 时间/min 100 95 90 85 800 20406080100120 TG / 10 K / min 20 K / min 30 K / min 回 ① A 图 1 褐铁矿热重曲线 2 数学模型 褐铁矿的磁化焙烧是为了使 Fe2O3转变成具有磁 性的 Fe3O4的过程。 在还原过程中,褐铁矿颗粒通过 周围流动的气体升温,CO 经过孔隙扩散后吸附于褐铁 矿表面,当颗粒温度达到所需反应温度时,温度和活化 能同时影响着化学反应速率。 所以用于描述褐铁矿磁 化焙烧过程的模型中必须包含颗粒的温度方程、反应 气体的传质方程、与温度和气体浓度相耦合的反应速 率方程和描述孔隙结构变化的方程。 模型假设 1) 反应物 Fe2O3与其他物质(杂质等)彼此分离。 2) 化学反应是一级的、不可逆的。 3) 褐铁矿颗粒为球对称结构,且各向同性。 4) 忽略某些杂质对还原反应的催化影响。 5) 忽略少部分过还原现象。 2.1 颗粒的传热方程 由于颗粒的固相占据着很大一部分,且气固相接 触充分,故只考虑固相的传热,气固反应放出的热量作 为方程的内热源。 ρcp ∂T ∂t = 1 r2 ∂ ∂r λeffr2 ∂T ∂r - R COΔH (2) 式中 ρ 为固相密度,kg/ m3; cp为固相热容,kJ/ (kgK); r 为褐铁矿颗粒半径,m;T 为颗粒温度,K;t 为时间,s; λeff为多孔铁矿颗粒的有效导热系数,W/ (mK);RCO 为化学反应速率,mol/ (m3s);ΔH 为反应热,J/ mol。 2.2 气体的扩散方程 在气固反应过程中,反应速率受到气体扩散和浓 度的影响,忽略 CO2在反应过程中解吸附后向颗粒外 的扩散,得出颗粒的非稳态气相传质方程 εC ∂YCO ∂t = 1 r2 ∂ ∂r r2DinC ∂YCO ∂r + R CO (3) 式中ε 为褐铁矿颗粒的孔隙率;YCO为铁矿颗粒内CO 反 应气体的质量分数;C 为环境 CO 浓度,mol/ m3;Din为褐 铁矿内有效扩散系数,m2/ s。 考虑 Knudsen 扩散的影响,则其内部的有效扩散 系数表达式为 Din= ε t(Dk -1 + D m -1 ) (4) 式中 Dk为 Knudsen 扩散系数,m2/ s;Dm为分子扩散系 数,m2/ s。 2.3 反应速率方程 多孔褐铁矿的气固反应过程采用随机孔模型[7],反 应过程中颗粒的反应速率与其反应比表面积成正比 g = KρgYCOsp(5) 401矿 冶 工 程第 35 卷 K = max k0(1 - xFe2O3) 1+ KFe 2O3 K (1 - x) 2KFe 2O3 3K 1 - ψln(1 - ψ)e - E RT,KFe 2O3 (6) Fe2O3转化率 xFe2O3为 xFe2O3= 1 - exp[- τYCO(1 + ψτYCO/4)](7) 式中 g 为褐铁矿还原速率,kg/ s; ρg为气相密度,kg/ m3; K 和 KFe2O3分别为褐铁矿和 Fe2O3反应速率常数,m/ s; sp为褐铁矿实际反应表面积,m2;k0为褐铁矿还原频率 因子,m/ s;x 为褐铁矿的转化率;ψ为孔的结构参数; τ 为无因次时间;E 为反应活化能,J/ mol;R 为通用气 体常数,J/ (molK)。 2.4 孔隙结构的变化 颗粒孔隙率 ε 由褐铁矿脱水后形成的大孔孔隙率 εmac和反应物自身所具有的小孔孔隙率 εmic组成 ε = εmac + ε mic (8) εmic = ε Fe2O3(1 - εmac - ε ash) (9) 其中 εash为灰分的孔隙率。 2.5 初始和边界条件 初始条件 t=0 时,x=0, T=T0,YCO =Y 0。 边界条件 r=0 时,∂T/ ∂r=0, ∂YCO/ ∂r=0。 r=r0时,-λeff∂T/ ∂r=(hc +h r)(T-Tsur), -D in∂YCO/ ∂r=Kex(YCO -Y CO, sur)。 考虑外界对颗粒的辐射换热,外部表面的对流传 质系数为 Kex = D inSh(1 - ε)(1 - βash) / dp (10) 式中 hc为自然对流换热系数,W/ (m2K);hr为辐射 换热系数,W/ (m2K);Kex为颗粒外气体传质系数, m/ s;Sh 为 Sherwood 数;Tsur为环境温度,K;βash为褐铁 矿中灰的质量分数;YCO,sur为外部环境 CO 质量分数; dp为褐铁矿的表观直径,m。 3 模拟与实验结果及讨论 表观活化能表征着反应堆规模和能量消耗,它是 褐铁矿还原过程中最重要的参数之一。 影响活化能的 因素有许多,它的值取决于反应初始物质组成、还原气 体特性、温度范围、反应级数、含有的杂质、物理形状 等[8]。 本文使用 NETZSCH 热力学软件处理反应测试 结果,以自由模型的概念开始动力学测评,先采用“无 模型”拟合模式来得到初步的反应活化能。 采用 3 种 不同升温速率的数据计算得到动力学模型,模型与实验 测试结果吻合得较好,其相关系数基本达到了 99%。 具 体动力学参数见表 2。 其中 lgA 表示一级反应指前因 子的对数值。 数值模拟中,活化能对还原过程影响非 常大,由表 2 可以看出,细颗粒的活化能要明显小于较 粗颗粒的活化能,这一点也可以说明文献[4]得出的 减小颗粒尺寸可以提高焙烧效果的原因。 表 2 动力学参数 粒度 / μm E / (kJmol -1 ) lgA 一级反应 级数 相关 系数 75131.298025.82270.509900.989878 425151.289715.490.304640.992095 3.1 颗粒粒径对还原过程的影响 图 2 给出了在 840 K 下,不同颗粒尺寸下还原褐 铁矿的模拟和实验对比图。 由图可见,实验与模拟数 据吻合得较好。 随着粒径增大,还原时间也相应增大。 425 μm 粗颗粒磁化焙烧所需时间(249 s)约为 75 μm 颗粒(96 s)的 2.5 倍,这是因为当颗粒增大时,颗粒反 应所需活化能增大。 在同样的孔隙结构和温度条件 下,褐铁矿在初期反应时,在初步形成的孔隙中,孔隙 表面一些活化能较低的活性位暴露于 CO 中,从而初 期无论颗粒大小都可以在较低的活化能下进行反应。 反应过程中总是活化能低的活性位先反应,然后活化 能高的依次反应[9]。 而大粒径褐铁矿在反应后期需 要更高的活化能,从而导致实验所测得的表观活化能 增加。 这应该是由于在反应后期,大颗粒中大量灰分 和反应产物包覆着反应物,导致更少的还原气体渗透 进入褐铁矿反应物表面,同时由于产物气体向外扩散 的行程变长,进一步导致还原气体进入颗粒内的阻力变 大。 考虑到流化床中气体流速更高,若在实际流化床中 进行颗粒还原,小颗粒对缩短还原反应时间的效果将更 加明显。 可以看出小颗粒所需的反应活化能更小,同时 与还原气氛直接接触的颗粒初始反应比表面也相应增 加,因此也更有利于还原反应的进行。 但是过细的颗粒 同样也会带来经济方面的问题。 1.2 1.0 0.8 0.6 0.4 0.2 0.00 50100150200250 425 μm 实验值 425 μm 计算值 75 μm 实验值 75 μm 计算值 时间/s 还原率/ 图 2 840 K 下还原褐铁矿模拟和实验对比图 501第 5 期朱 超等 褐铁矿粉磁化焙烧的实验与模拟研究 当环境温度、初始孔隙率以及气体环境全部相 同时,褐铁矿颗粒粒径大小是影响还原过程的重要因 素。 颗粒大小对还原过程的影响首先将体现在颗粒内 部一氧化碳的浓度分布上。 如图 3 所示,环境温度为 843 K,初始孔隙率均为 0.4,时间步长和迭代次数相同 时,不同粒径褐铁矿颗粒内部存在浓度梯度,当褐铁矿 颗粒粒径为 50 μm 时,气体可以很好地扩散到褐铁矿 颗粒内部;当粒径为 150 μm 以及 450 μm 时,气体到 达中心的浓度明显低于粒径 50 μm 的颗粒,这说明褐 铁矿颗粒粒径越大,气体越难到达颗粒中心。 这是因 为,其他条件相同时,褐铁矿颗粒粒径越大,气体在颗 粒内部传输所需的时间越长,导致气体在褐铁矿颗粒 内部消耗的量变多,而随着反应进行,内部孔隙结构不 断发生变化,使得大颗粒的内部阻力变大,所以褐铁矿 颗粒粒径越大,气体越难进入到颗粒内部,从而进一步 减缓了颗粒的还原过程。 浓度/ mol m-3 r / R 回回 回 回 回 回 回 @@@@ @ @ @ AAAAA 回 @ A A A 2.0 1.6 1.2 0.8 0.4 0.0 0.00.20.40.60.81.0 50 μm 150 μm 450 μm 图 3 不同粒径下沿褐铁矿颗粒径向方向的 CO 浓度分布 3.2 还原气体浓度对还原过程的影响 粒径 75 μm 的褐铁矿颗粒,温度 840 K,环境 CO 浓度 1~3 mol/ m3,模拟不同环境 CO 浓度下与磁化反 应完全所需要时间的关系如图 4 所示。 CO浓度/mol m-3 回 回 回 回 回 250 200 150 100 50 0 1.01.52.02.53.0 反应时间/s 图 4 反应完全时间随 CO 浓度的变化 由图 4 可知,CO 与褐铁矿颗粒外表面和内部孔隙 壁面进行还原反应,褐铁矿的还原率是初始环境 CO 浓 度的函数,当 CO 浓度为 1.5~3.0 mol/ m3时,随着初始 环境 CO 浓度升高,褐铁矿的还原时间近似呈单调线性 下降。 当 CO 浓度由 1.5 mol/ m3减少到 1 mol/ m3时, 还原时间大幅度增加,说明在 840 K,CO 浓度低于 1.5 mol/ m3的情况下,还原气体浓度对还原反应已起 到了关键作用。 随着环境 CO 浓度增加,CO 分压力增 大,更有利于颗粒内部的还原反应。 当环境 CO 浓度 由 1 mol/ m3增大到 3 mol/ m3时,还原时间减少了 73%。 但需要指出的是,由于 Fe2O3向具有磁性的 Fe3O4转变非常容易发生,过高的 CO 浓度反而会导致 Fe3O4的过还原[4]。 3.3 孔隙率的影响 图 5 为褐铁矿颗粒在不同初始孔隙率下转化率随 时间的变化特性。 随着孔隙率增大,还原完成所需的 时间也相应缩短。 这是因为孔隙率增大,更有利于还 原气体向内部传递,即孔隙率增大,导致内扩散阻力减 小,从而使反应速率增大[10-11]。 1.2 1.0 0.8 0.6 0.4 0.2 0.00 50100150200250 0.3 0.4 0.5 时间/s 还原率/ 图 5 不同初始孔隙率下转化率随时间的变化 图 6 显示了颗粒内部一氧化碳浓度分布随孔隙率 的变化。 由图 6 可见,当颗粒孔隙率增大时,颗粒内部 还原气体的浓度梯度减小,即颗粒孔隙率增大,使得颗 粒内部的扩散阻力减小,从而使还原气体更易进入颗 粒内部,也使反应速率相应增大。 这与前述还原时间 分析结果相对应。 CO浓度/ mol m-3 r / R 1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.00.20.40.60.81.0 0.3 0.35 0.4 0.45 0.5 图 6 一氧化碳浓度分布随孔隙率的变化 (下转第 109 页) 601矿 冶 工 程第 35 卷 含铁、镁的稳定性较强。 由于浸出时间直接决定工艺 成本及效率,综合考虑,浸出时间控制为 180 min。 浸出率/ 回 回 回 回 回 回 @ @ A A A A @A A @A A @A A @A A AA A A AA AA A A A A 旁 旁 旁 旁 旁 旁 100 80 60 40 20 4 2 0 60120180240300 浸出时间/min Co Fe Cu Mg Ni Al Mn 1 图 2 浸出时间对浸出过程的影响 3 结 论 1) 对高品位水钴矿焙烧预处理后,其适宜的高压 酸浸条件为浸出剂硫酸浓度 0.2 mol/ L,浸出液固比 10∶1,浸出时间 180 min。 在此条件下,钴、铜和镍的浸 出率分别为 92.8%、96.7%和 95.3%。 2) 高压酸浸条件下,铁、镁浸出率低,可实现水钴 矿中有价元素钴、铜、镍与杂质铁、镁的初步分离。 参考文献 [1] 郭学益,姚 标,李晓静,等. 水钴矿中选择性提取铜和钴的新工 艺[J]. 中国有色金属学报,2012,22(6)1778-1784. [2] 刘忠胜,邢 飞,段英楠. 某钴铜精矿硫酸化焙烧试验研究[J]. 矿冶工程,2014(5)108-112. [3] Michel Shengo Lutandula, Banza Maloba. Recovery of cobalt and cop⁃ per through reprocessing of tailings from flotation of oxidised ores[J]. Journal of Environmental Chemical Engineering,2013(1)1085-1090. [4] Clotilde Apua M, Mulaba⁃Bafubiandi A F. Dissolution of oxidised Co⁃ Cu ores using hydrochloric acid in the presence of ferrous chloride[J]. Hydrometallurgy,2011,108233-236. [5] Mohammad Sadegh Safarzadeh, Nikhil Dhawan, Mustafa Birinci, et al. Reductive leaching of cobalt from zinc plant purification residues[J]. Hydrometallurgy,2011,10651-57. [6] 刘 俊,李林艳,徐盛明,等. 还原酸浸法从低品位水钴矿中提取 铜和钴[J]. 中国有色金属学报,2012,22(1)304-309 [7] 李淑梅. 刚果(金)水钴矿浸出新方法及热力学分析[J]. 中国有 色冶金,2014(2)29-32. [8] 马文强,徐盛明,李林艳,等. 用 Lix984 萃取分离水钴矿浸出液中 的铜钴[J]. 矿物学报,2011,31(3)560-565. [9] 刘建华,张焕然,王瑞祥,等. 氨法加压浸出钴铜氧化矿工艺[J]. 稀有金属,2012,36(1)149-153. [10] Xueyi Guo, Dong Li, Kyung⁃Ho Park, et al. Leaching behavior of metals from a limonitic nickel laterite using a sulfation⁃roasting⁃leac⁃ hing process[J]. Hydrometallurgy,2009,99144-150. �������������������������������������������������������������������������������������������������� (上接第 106 页) 4 结 论 利用同步热分析仪研究了 CO 在 823~873 K 之间 还原不同颗粒尺寸褐铁矿的过程,利用传热、传质与化 学反应耦合的动力学模型模拟分析了反应颗粒尺寸、 反应气体浓度及颗粒孔隙率对磁化焙烧的影响,实验 与模拟结果吻合得较好。 得到如下结论 1) 相对较粗颗粒,细颗粒的表观活化能减少了近 20 kJ/ mol,从反应活性方面证明了细颗粒还原焙烧效 果更好。 2) 颗粒粒径越大,气体越难进入到颗粒内部,结 果表明,425 μm 粗颗粒磁化焙烧所需时间约为 75 μm 细颗粒的 2.5 倍。 3) CO 浓度在 1.5~3.0 mol/ m3之间时,随 CO 浓 度升高,褐铁矿还原时间近似呈单调线性下降;CO 浓 度由 1.5 mol/ m3减少到 1 mol/ m3时,还原时间大幅度 增加。 当环境 CO 浓度由 1 mol/ m3增大到 3 mol/ m3 时,还原时间减少了 73%。 4) 孔隙率增大,颗粒内部还原气体的浓度梯度减 小,导致内扩散阻力减小,从而使反应速率增大。 参考文献 [1] 朱俊士. 选矿试验研究与产业化[M]. 北京冶金工业出版社, 2004. [2] 张 茂,王 东,陈启平,等. 云南某褐铁矿磁化焙烧⁃磁选工艺试 验研究[J]. 矿冶工程,2011(6)51-53. [3] 刘兴华,罗良飞,刘 卫,等. 某低品位褐铁矿的选矿工艺研究 [J]. 矿冶工程,2013(4)70-73. [4] Yu Y F,Qi C Y. Magnetizing roasting mechanism and effective ore dressing process for oolitic hematite ore[J]. Journal of Wuhan Uni⁃ versity of Technology⁃Materials Science Edition, 2011,26(2)176- 181. [5] Chen J G, Sun T C, Lu Q. Study of the mechanism of the direct re⁃ duction roasting of the limonite in Jiangxi[J]. Journal of Chemical and Pharmaceutical Research,2014,6(3)671-678. [6] 朱 超,龚志军,李保卫,等. 多孔褐铁矿磁化还原特性的模型构 建及数值模拟[J]. 过程工程学报,2014,14637-642. [7] Pineau A, Kanari N, Gaballah I. Kinetics of reduction of Iron oxides by H2 Part I Low temperature reduction of hematite[J]. Thermochim Acta,2006,447(1)89-100. [8] Wang X H, Zeng X J, Yang H L, et al. General modeling and nu⁃ merical simulation of the burning characteristics of porous chars[J]. Combust Flame,2012,159(7)2457-2465. [9] 廖 坚. 煤焦颗粒燃烧特性实验和数值模拟研究[D]. 哈尔滨哈 尔滨工业大学能源科学与工程学院,2010. [10] 郭 明,商志才,俞庆森,等. 化合物分子体积计算方法的研究 [J]. 浙江大学学报(理学版),2003(5)554-560. [11] 严继民. 吸附与凝聚⁃固体的表面与孔(第二版)[M]. 北京科 学出版社,1986. 901第 5 期郭秋松等 刚果(金)水钴矿高压酸浸实验研究