建立了铜闪速熔炼过程气粒两相间砷分配行为的数学模型,基于软件Fluent 6.3对其分配行为进行了数值仿真。仿真结果表明:Cu 3AsS 4的脱砷过程和As 2O 3的挥发过程均非常迅速,精矿脱砷过程在距离反应塔顶3 m的区域内基本完成,烟灰中As 2O 3的挥发过程在距离反应塔顶1 m的区域内基本完成;距离塔顶3 m以内的区域是调整砷在各相中分配的关键区域,对于所计算的工况,砷在熔体与气相中的分配比例分别为48.1%、51.9%。 The numerical simulation model for the arsenic distribution behavior between gas and particle phases in a copper flash smelting process was established, and the distribution behavior was nu-merically simulated based on the software Fluent 6.3. The simulating results show that the arsenic was separated rapidly from Cu 3AsS 4 and As 2O 3 was quickly volatilized. The arsenic separation from the concentrate was almost completed within 3 m from the top of the reaction tower. The volatilization of As 2O 3 in the soot was finished within 1 m from the top of the reaction tower. The zone within 3 m from the top of the tower is the key to adjust the distribution of arsenic in each phase. For the simulated case, the proportion of arsenic distribution is 48.1% in the melt and 51.9% in the gas phase.
龙 鹏1,2*,刘达兵1,陈 卓1
1中南大学能源科学与工程学院,湖南 长沙
2桂林航天工业学院能源与建筑环境学院,广西 桂林
Email: *longpeng@guat.edu.cn
收稿日期:2018年5月13日;录用日期:2018年6月5日;发布日期:2018年6月12日
摘 要
建立了铜闪速熔炼过程气粒两相间砷分配行为的数学模型,基于软件Fluent 6.3对其分配行为进行了数值仿真。仿真结果表明:Cu3AsS4的脱砷过程和As2O3的挥发过程均非常迅速,精矿脱砷过程在距离反应塔顶3 m的区域内基本完成,烟灰中As2O3的挥发过程在距离反应塔顶1 m的区域内基本完成;距离塔顶3 m以内的区域是调整砷在各相中分配的关键区域,对于所计算的工况,砷在熔体与气相中的分配比例分别为48.1%、51.9%。
关键词 :砷分配,闪速熔炼,数值仿真
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
砷为我国重点防控的五种重金属污染物之一,对环境和人体的危害已众所周知。铜冶炼企业作为有色金属冶炼行业的重要组成部分,一直以来是产砷排砷的大户,砷在冶炼工序中可转化为挥发性和毒性极大的三氧化二砷,腐蚀闪速炉炉壁耐火材料,可使SO2转化触媒五氧化二钒中毒,可转化为无色剧毒的砷化氢气体,在电解过程可导致阳极钝化,影响阴极铜的生产和质量 [
目前对有色金属冶炼过程中砷分配行为的研究大多以铜冶炼过程为代表,采用基于化学热力学和化学平衡理论的热力学建模法进行 [
本文以我国某炼铜企业的闪速熔炼炉为原型,建立数值仿真物理结构模型。由于研究中着重关注砷元素在悬浮熔炼条件下的反应与分配行为特点,因此计算区域仅包含闪速熔炼炉反应塔、沉淀池部分气相空间以及反应塔顶部的精矿喷嘴、天然气烧嘴等,如图1所示。
用于描述闪速熔炼过程的数学模型有 [
图1. 铜闪速炉物理模型示意图。1-精矿喷嘴;2-塔顶天然气烧嘴;3-反应塔;4-沉淀池气相空间
颗粒表面反应:
2CuFeS 2 + O 2 = Cu 2 S + 2FeS + SO 2 (1)
FeS 2 + O 2 = FeS + SO 2 (2)
3FeS + 5O 2 = Fe 3 O 4 + 3SO 2 (3)
2Cu 2 S + 3O 2 = 2Cu 2 O + 2SO 2 (4)
FeS + 3Fe 3 O 4 = 10 FeO + SO 2 (5)
3Cu 2 O + FeS = 6Cu + FeO + SO 2 (6)
4Cu 3 AsS 4 + 13O 2 = 6Cu 2 S + As 4 O 6 + 10 SO 2 (7)
2As 2 O 3 = As 4 O 6 (8)
4Fe 3 O 4 + 3As 4 O 6 + 7O 2 = 12FeAsO 4 (9)
体积反应:
CH 4 + 2O 2 = CO 2 + 2H 2 O (10)
As 4 O 6 = 4AsO + O 2 , T > 1073 K (11)
三种边界条件:1) 各气流入口均设置为质量流量入口,以某炼铜企业实际操作参数为依据,各股气流的参数如表1所示;2) 对于动量方程,计算区域内各壁面均设置为无滑移壁面,对于能量方程,壁面温度根据闪速熔炼炉炉内温度分布特点及与其接触的烟气温度进行设置;3) 烟气出口设置为压力出口,根据在线监测数据,出口处的压力设置为−42 Pa。
闪速熔炼炉投料量为60.56 kg/s,通过取样并进行相关实验分析,其配料、特征粒径、颗粒物相组成分别如表2~表4所示。
工况 | 工艺风 | 分散风流量(kg/s) | 中央氧流量(kg/s) | 中央天然气流量(kg/s) | 塔顶烧嘴(3个) | |||
---|---|---|---|---|---|---|---|---|
工艺风量(kg/s) | 氧浓(%) | 风速 (m/s) | 天然气量(kg/s) | 燃烧风量(kg/s) | ||||
基准工况 | 14.43 | 76.8 | 90 | 1.129 | 0.874 | 0.03 | 0.03 | 0.538 |
表1. 闪速炉配风参数
炉料 | 精矿 | 渣精矿 | 吹炼渣 | 石英 | 烟灰 | |
---|---|---|---|---|---|---|
质量百分数(%) | 71.42 | 3.17 | 4.76 | 12.39 | 8.26 | |
特征粒径(μm) | 24.642 | 37.0 | 926.873 | 429.548 | 24.642 | |
表2. 闪速炉配料参数及特征粒径
成分 | CuFeS2 | FeS2 | FeS | Fe3O4 | Cu2O | SiO2 | Fe2SiO4 | As2O3 | Cu3AsS4 |
---|---|---|---|---|---|---|---|---|---|
质量百分数(%) | 66.94 | 4.80 | 1.76 | 2.79 | 0.81 | 6.19 | 13.73 | 0.66 | 1.57 |
表3. 精矿、烟灰的物相组成
成分 | Cu2S | Fe3O4 | Cu2O | SiO2 | Fe2SiO4 | FeAsO4 |
---|---|---|---|---|---|---|
质量百分数(%) | 29.45 | 7.24 | 5.96 | 15.67 | 35.92 | 5.17 |
表4. 渣精矿的物相组成
闪速熔炼过程中,含砷物质Cu3AsS4、As2O3的反应过程及中间产物As4O6的转化过程决定着砷的分配行为,控制着砷的流向。下面将根据仿真结果对此进行逐一分析。
气相中As4O6和AsO的浓度分布可以清楚地反映出砷在气相空间内的转化过程,闪速炉中心对称面上As4O6和AsO的气相浓度分布如图2所示。
As4O6和AsO在炉内的总体反应过程如图3所示。根据反应速率相对大小的不同,As4O6的生成与消耗过程可以分为以下三个阶段。
1) 生成速率大于消耗速率阶段:这一过程发生在距离塔顶0.75 m的区域内,此时精矿和配风刚入炉,温度较低,主要以As4O6的生成为主,在离塔顶0.75 m处时As4O6的量到达最大值。此时,虽然精矿还未着火,但在中央天然气燃烧产生的高温区域的作用下,部分As4O6开始分解,从图3(b)可以看出,在离塔顶0.75 m处时,炉内已经生成了一定的AsO。
2) 消耗速率大于生成速率阶段:这一过程发生在距离塔顶0.75~2.5 m的区域内,伴随着精矿的着火和中央天然气的燃烧,气相的温度急剧升高,在达到分解温度后,As4O6的生成速率便将大大低于其消耗的速率。图3(a)表明在距离塔顶2.5 m时,炉内原来积累的As4O6已消耗殆尽。
3) 生成速率等于消耗速率阶段:这一过程发生在距离塔顶2.5~3.25 m的区域,此时由于精矿已经完全着火并开始燃烧,气相的温度远超过As4O6的分解温度,致使As4O6在生成的一瞬间即被分解和消耗。
图2. 闪速熔炼炉对称面上的气相浓度场分布云图
图3. 炉内As4O6和AsO的变化趋势:(a) 炉内As4O6的生成与消耗趋势;(b) 炉内AsO的生成趋势;(c) 炉内AsO的净流量
由图3(c)可以看出,在这一阶段由As4O6转化而来的AsO变化很小,可以认为As4O6的生成与转化过程在距离塔顶2.5 m的区域内已基本完成。
从图3(c)中可以看出,在距离塔顶2 m以内的区域,生成了超过95%的AsO,其中,距离塔顶1 m以内区域生成的AsO约占40%,距离塔顶1~2 m之间的区域生成的AsO约占60%,在距离塔顶3 m以下的区域,不再有AsO生成。据此可以推论:距离塔顶1~2 m之间的区域是生成AsO的主要区域,含砷物质的反应过程在距离塔顶3 m的位置处已基本完成。
混合精矿中Cu3AsS4在闪速熔炼过程中化学反应如式(9)所示,其反应过程如图4(a)、图4(b)所示。从图中可以看出,精矿中的Cu3AsS4入炉后立即发生化学反应,砷以As4O6的形式进入气相,铜以Cu2S
图4. Cu3AsS4、As2O3的反应过程。(a) 混合精矿、烟灰中Cu3AsS4质量分布;(b)Cu3AsS4的质量分数沿反应塔高度的变化;(c) 混合精矿、烟灰中As2O3质量分布;(d) As2O3的质量分数沿反应塔高度的变化;(e) 混合精矿、烟灰中FeAsO4质量分布;(f)FeAsO4的质量分数沿反应塔高度的变化
砷带入量(t/h) | 砷带出量(t/h) | |||
---|---|---|---|---|
干矿 | 烟灰 | 熔体e(炉渣+冰铜) | 气相(烟气、烟尘) | |
0.732 | 0.866 | 0.768 | 0.827 | |
质量百分数(%) | 45.8% | 54.2% | 48.1% | 51.9% |
表5. 仿真得到的砷在各相中的分配比
注e:熔体即为反应后落入沉淀池的颗粒相(精矿颗粒、烟灰、渣精矿、吹炼渣和石英)。
的形式继续存在于颗粒中;整个反应过程在精矿喷嘴附近的富氧环境下进行,在距离塔顶2 m的位置处已基本反应完全,转化过程迅速而稳定。
烟灰中的As2O3在塔内发生的化学反应如式(8)所示,转变过程如图4(c)、图4(d)所示。从图中可以看出,在炉内高温环境下,As2O3的挥发过程极为迅速,在距离塔顶1 m的位置处时已全部以As4O6的形式进入气相。
FeAsO4的生成反应如式(9)所示,生成过程如图4(e)、图4(f)所示。图4(e)表明,料柱中央的精矿颗粒FeAsO4的含量较高,而料柱外围的精矿颗粒FeAsO4的含量较低,这是因为位于塔中央的精矿颗粒更为集中,升温着火慢,粒径较小容易氧化造成的。从图4(f)可以看出:进入熔体的FeAsO4同进入气相的AsO一样,主要在距离塔顶3 m以上的区域内生成,而在3 m以下的区域,几乎不再有FeAsO4生成;精矿颗粒中FeAsO4含量的提高主要是由于精矿的进一步脱硫造成。
闪速熔炼过程的数值仿真计算得到的砷在熔体及气相中的分配比例列于表5中。这里,由于烟尘随同烟气一起进入后续的工序,并不进入位于沉淀池的熔体中,因此将烟尘的带走的砷算在气相中。从中可知,在熔炼过程,最终有51.9%的砷进入熔体,48.4%的砷进入气相,而这和实际生产过程中收集的数据相符。
建立了描述铜闪速熔炼过程中砷分配行为数学模型,并对其进行了数值仿真,揭示了铜闪速熔炼过程中含砷物质的转化过程。主要结论如下:
1) Cu3AsS4的脱砷过程和As2O3的挥发过程都非常迅速,精矿中含砷物质Cu3AsS4的脱砷过程在距离反应塔顶3 m的区域内已基本完成,而烟灰中的含砷物质As2O3的挥发过程在距离反应塔顶1 m的区域内已基本完成;
2) 距离塔顶3 m以内的区域是调整砷进入熔体或气相比例的关键区域,闪速熔炼炉出口处砷在熔体中的分配比例为48.1%,在气相中的分配比例为51.9%。
研究尚未就砷在铜闪速熔炼过程中分配行为影响因素进行深入分析,后续将就此进行进一步研究。
环保公益性行业科研专项经费项目(NO. 201509050)。
龙 鹏,刘达兵,陈 卓. 铜闪速熔炼砷分配行为的数值仿真 Numerical Simulation of Arsenic Distribution Behavior in Copper Flash Smelting[J]. 冶金工程, 2018, 05(02): 85-92. https://doi.org/10.12677/MEng.2018.52012