有色金属科学与工程  2019, Vol. 10 Issue (2): 6-12
文章快速检索     高级检索
树脂吸附铜锌离子的数值模拟[PDF全文]
陈尚1 , 刘青1 , 赵世强1 , 郑其2     
1. 北京科技大学钢铁冶金新技术国家重点实验室,北京 100083;
2. 北京有色金属研究总院稀有金属冶金材料研究所,北京 100088
摘要:镀金原料氰化亚金钾的生产过程含氰废液中锌(Ⅱ)离子和铜离子(Ⅱ)浓度较高,为探究树脂吸附铜、锌离子的规律,运用Aspen Adsorption数值模拟软件对Φ10 mm×200 mm的树脂固定床吸附过程进行数值模拟,采用一阶迎风差分法作为偏微分方程的离散方法,研究不同条件下的吸附穿透曲线和吸附负荷分布曲线.结果表明:增加进料流量和进料浓度均可提高吸附效率,传质系数(KS)大于临界值时,吸附穿透曲线不再受其影响,KS(Zn2+)=4.3×10-3s-1, KS(Cu2+)=1.00×10-2s-1.
关键词树脂    吸附柱    锌离子    铜离子    穿透曲线    
Mathematical simulation on adsorption of zinc ion and copper ion with resin
CHEN Shang1 , LIU Qing1 , ZHAO Shiqiang1 , ZHENG Qi2     
1. State Key Laboratory of Advanced Metallurgy, University of Science and Technology Beijing, Beijing 100083, China;
2. Rare Metals and Metallurgy Materials Research Institute, General Research Institute for Nonferrous Metals, Beijing 100088, China
Abstract: Zinc(Ⅱ) ion and Copper(Ⅱ) ion come from the process of producing potassium aurocyanide, gold-plated material. Adsorption process of zinc ion and Copper ion resin in a φ10 mm×20 mm adsorption column was simulated and studied with the aid of Aspen Adsorption. Ion concentration at exit of adsorption column vs time curve and bed axial loading distribution of ion were obtained from simulation results with discrete partial differential equations of first order upwind difference method. The time of breakthrough was shortened, as feed flow rate and feed concentration increase respectively. The adsorption breakthrough curves did not change with mass transferring coefficient any more When the mass transfer coefficient is greater than the critical value, KS(Zn2+)=4.3×10-3 s-1, KS(Cu2+)=1.00×10-2 s-1. The modeling study provides reference data for scale-up experiment and industrial application.
Keywords: resin    adsorption column    zinc ion    copper ion    breakthrough curve    

氰化亚金钾是集成线路板或工艺品的主要镀金原料,制备过程产生的废液经活性炭吸附回收氰化亚金离子后,废液中的有价金属锌的含量高达0.310 g/L、铜的含量高达0.006 g/L.树脂对提金氰化废液进行综合处理是主流方法之一,具有重复使用率高、处理容量大、能够处理多种重金属离子和酸根等优势[1-7].实验中的静态吸附与动态实验结合研究的方法不便考察传质系数对吸附回收效率的作用,而数值模拟的方法可以不受此限制,模拟吸附过程,从而节约成本,提高效率[8-14].

彭松水[15]建立CO2待净化体系的双塔变温吸附模型,借助Aspen Adsorption软件动态模拟各组分的吸附、热再生、冷再生过程,对不同类型的吸附剂做出有效评价,同时对吸附剂制定较优的时序组合.李建伟等[16]运用Aspen Adsorption模拟煤层气的吸附非等温过程,研究床层温度、甲烷和氮气的吸附量以及气速等参数的分布规律.杨中山等[17]建立Φ0.2 m×0.5 m的吸附柱模型,运用Aspen Adsorption考察不同条件下氰化亚金离子的吸附柱穿透曲线.刘金昌等[18]依靠Aspen Adsorption软件对单塔活性炭吸附锌离子过程,研究浓度、流量和传质系数对吸附过程的影响.杨金垚等[19]运用Aspen Adsorption模拟烟气中CO2在固定床中的动态变温吸附过程,探究温度、时间对CO2捕捉率、产品纯度和分离耗散能的影响.薛才红等[20]以乙醇蒸汽为原料,运用Aspen Adsorption,探究了吸附压力、吸附剂体积对乙醇变压吸附过程的影响.竹涛等[21]运用Aspen Adsorption探究煤层气的变压吸附分离过程中原料气浓度、吸附压力对解吸气体甲烷浓度的影响.

文中运用Aspen Adsorption软件,以国内某金盐生产企业的提金氰化废液处理装置为研究对象,建立实验室规模的吸附柱,并对吸附过程进行数学建模,比较树脂吸附柱出口处锌离子的动态浓度数据与Aspen Adsorption模拟数据,验证了模型的可靠性.在此基础上,对树脂吸附铜、锌离子过程,分别从进料流量、进料浓度、传质系数及吸附柱的轴向负载等方面对吸附规律进行研究.

1 吸附数学模型的建立

为建立吸附过程的数学模型,首先应对实际吸附过程进行假设以简化问题,然后再对传质过程建立数学方程,选择合理的数学模型,最后借助Aspen Adsorption软件并给出初始化条件,最终对树脂吸附铜、锌离子的过程进行数值解析.

1.1 模型假设

针对树脂吸附铜、锌离子的物理过程,对固定床吸附铜、锌离子进行如下假设[17]:①液相均没有黏性且为不可压缩流体;②无温度变化,吸附过程属于绝热状态;③各吸附床层内的液相流动均为扩散活塞流;④不考虑轴向的物质扩散和压力降.

1.2 数学模型 1.2.1 物料平衡方程

采用离散方法求解吸附过程,将吸附柱离散化,分割成无数个单元微元体,对于任意微元体,吸附质的输入速率与输出速率的差值应与吸附质组分在利用微元体内的积累速率相等,利用吸附质浓度的变化关系可得到整个吸附柱的物料平衡方程:

(1)

整理得:

(2)

由模型假设可知,此流动为理想活塞流模型,即可忽略轴向扩散的影响,式(2)中的等号右侧等于零,物料平衡方程简化[17]为:

(3)

式(3)中,v是流体的线速度(m/s),C是组分浓度(kmol/m3),x是轴向距离(m),ε是孔隙率,t是吸附时间(s),ρ是堆积密度(kg/m3),q为吸附量(kmol/kg).初始条件:

(4)
(5)
(6)
1.2.2 吸附速率方程

固体床吸附采用线性推动力模型(LDF)描述,其吸附速率方程[22]为:

(7)

式(7)中,ke是集成扩散系数(s-1),q0是饱和吸附量(kmol/kg).

当表面扩散为内扩散的限制环节时,

(8)

式(8)中,Ds是表面扩散系数(m2/s);Rp是吸附剂颗粒半径(m).

1.2.3 数值方法

固体床吸附过程的离散方法采用快速、稳定、易收敛的一阶上风差分法(USD1,Upwind Differencing Scheme 1),其基于泰勒的一节展开式的推导方程[15]为:

(9)
2 参数设置 2.1 参数确定

Langrmuir方程和Freundlich方程均可用来描述吸附剂和吸附质间的吸附行为,其中

(10)
(11)

式(10)、式(11)中:Ce为平衡浓度(mg/L); qe为平衡吸附量(mg/g); qm为饱和吸附容量(mg/g); KL为Langrmuir吸附平衡常数(L/mg); KF为Freundlich吸附平衡常数(L1/n ·mgn-1/n· g-1).

使用树脂对不同浓度的锌离子溶液进行静态吸附实验,获得平衡浓度和平衡吸附量,然后分别用Langmuir方程和Freundlich方程对静态吸附数据进行线性拟合,RF2(0.998)>RL2(0.912),如图 1图 2所示.因此,本研究吸附等温线测定实验的有限实验数据说明了与树脂吸附锌离子过程较为吻合的是Freundlich吸附等温模型,而Freundlich吸附等温模型主要拟合物理吸附过程,可知树脂吸附锌离子过程可能主要依靠物理吸附进行,其吸附等温方程为:

图 1 Freundlich方程拟合吸附等温线 Fig. 1 Adsorption isotherm of Freundlich equation

图 2 Langmuir方程拟合吸附吸附等温线 Fig. 2 Adsorption isotherm of Langmuir equation

(12)

同理,将树脂对铜离子的静态吸附数据用2种吸附模型进行拟合. Langmuir方程的拟合系数RL2=0.9871,Freundlich方程的拟合系数RF2=0.9705.因此选择Langmuir模型,树脂对铜离子的吸附为单分子层吸附[23].而Langmuir模型主要拟合化学吸附或氢键物理吸附,树脂吸附铜离子时可能发生了离子交换的化学吸附.

2.2 吸附柱相关参数的确定

进料口中参数主要包括进料流量、进料浓度、温度和压强,等.表 1列出了树脂吸附铜、锌离子时的进料口参数值.

表 1 铜、锌离子进料参数 Table 1 Feeding parameters of zinc ion and copper ion
点击放大

吸附床层的参数主要包括吸附床层的高度、直径,吸附剂的堆积密度、空隙率,吸附等温参数,等.如表 2所列.

表 2 树脂吸附铜、锌离子时吸附床层主要参数 Table 2 Parameters of adsorbent bed filled with resin for adsorption of zinc ion and copper ion
点击放大

3 模拟结果与讨论

本研究通过数值模拟数据与已有的实验数据相对比,来验证模型的准确性.然后,利用Aspen Adsorption软件分别从进料流量、进料浓度、传质系数方面研究树脂吸附铜、锌离子的规律.

3.1 模型验证

运用Aspen Adsorption软件进行计算,并与实验数据对比,如图 3所示.当模型对进料锌离子浓度为4.73×10-3 kmol/m3,进料流量为4.6×10-5 m3/h,模拟结果能够正确预报吸附柱出料处锌离子浓度变化趋势, 与实验结果比较吻合,验证了树脂吸附模型的可靠性.因此,可以用来模拟预测树脂吸附铜、锌离子的情况.

图 3 数值模拟结果与实验结果的吸附穿透曲线的对比 Fig. 3 Comparison of adsorption breakthrough curves between numerical simulation results and experimental results

3.2 进料浓度对吸附过程的影响

固定床进料流量为4.6×10-5 m3/h时,锌离子含量分别为1.73×10-3 kmol/m3、4.73×10-3 kmol/m3、7.73×10-3 kmol/m3时的吸附穿透曲线如图 4所示.

图 4 不同进料浓度下的锌离子吸附穿透曲线 Fig. 4 Breakthrough curves of zinc ion at different feed concentrations

进料中锌离子的浓度从1.73×10-3 kmol/m3增至4.73×10-3 kmol/m3时,吸附床层开始穿透的时间从1.06 h缩短至1.00 h;当进料浓度从4.73×10-3 kmol/m3增加到7.73×10-3 kmol/m3,吸附床层开始穿透的时间从1.00 h缩短至0.98 h.因此,固定床层的开始穿透时间缩短幅度随着浓度的升高而降低.锌离子进料浓度从1.73×10-3 kmol/m3增至4.73×10-3 kmol/m3时,完全穿透的时间由13.11 h缩短至12.50 h;进料浓度从4.73×10-3 kmol/m3增至7.73×10-3 kmol/m3,吸附床层完全穿透的时间从12.5 h缩短至12.28 h.

当铜离子进料流量为4.6×10-4 m3/h时,取进料中离子的浓度分别取2.82×10-6 kmol/m3、4.82×10-6 kmol/m3、6.82×10-6 kmol/m3、8.82×10-6 kmol/m3、1.08×10-5 kmol/m3时,吸附穿透曲线如图 5所示.

图 5 不同进料浓度下铜离子吸附穿透曲线 Fig. 5 Breakthrough curves of copper ion at different feed concentrations

随着进料中铜离子浓度的增加,吸附开始穿透的时间缩短.以2.82×10-6 kmol/m3、6.82×10-6 kmol/m3、1.08×10-5 kmol/m3为例进行分析,当进料浓度由2.82×10-6 kmol/m3增至6.82×10-6 kmol/m3时,开始穿透的时间从2.00 h缩至1.72 h;当进料浓度由6.82×10-6 kmol/m3增至1.08×10-5 kmol/m3时,开始穿透的时间从1.72 h缩至1.58 h.

因此,进料浓度的增加会缩短完全穿透的时间.随着浓度的增加,吸附穿透曲线的斜率和吸附床层被完全穿透的速率均会增大,原因是:随着进料中锌离子的浓度的提高,单位时间进入床层的离子增大,而吸附床层饱和吸附量增大的幅度不够吸附增加的吸附质.部分锌、铜离子未被吸附而进入下一阶段床层,传质区迅速向床层底部移动,床层被穿透的速度增大.

3.3 进料流量对吸附过程的影响

进料中锌离子的含量为4.73×10-3 kmol/m3,进料流量分别为2.6×10-5 m3/h、4.6×10-5 m3/h和6.6×10-5 m3/h时,绘制吸附穿透曲线,如图 6所示.

图 6 不同进料流量下锌离子的吸附穿透曲线 Fig. 6 Breakthrough curves of zinc ion at different feeding flow rates

当进料流量由2.6×10-5 m3/h增至4.6×10-5 m3/h时,吸附床层开始穿透时间从2.83 h缩短至1.00 h.完全穿透时间从19.5 h缩短至12.5 h;当进料流量由4.6×10-5 m3/h增至6.6×10-5m3/h时,吸附床层开始穿透时间从1.00 h缩短至0.36 h.完全穿透时间从12.50 h缩短至9.50 h.

铜离子进料浓度是9.44×10-5 kmol/m3,进料流量分别为2.6×10-5 m3/h、4.6×10-5 m3/h、6.6×10-5 m3/h、8.6×10-5 m3/h时的吸附穿透曲线,如图 7所示.当进料流量由2.6×10-5 m3/h增加到4.6×10-5 m3/h时,吸附床层开始穿透的时间从3.08 h缩短到1.39 h;完全穿透时间从9.50 h缩短至6.00 h;当进料流量从4.6×10-5 m3/h增加至6.6×10-5 m3/h时,开始穿透时间由1.39 h缩短到0.75 h,完全穿透时间从6.00 h缩短至4.50 h.

图 7 不同流量下铜离子吸附穿透曲线 Fig. 7 Breakthrough curves of copper ion at different feeding flow rates

数据表明:在其他条件不变的情况下,随着进料流量的增加,开始、完全穿透时间均提前.原因是,随着进料流量的增加,部分锌离子未被树脂吸附就进入下一床层,吸附质与吸附剂接触的时间在减少,传质时间缩短,传质区移动到床层底部的速度增加.

3.4 传质系数对吸附过程的影响

当锌离子的进料流量为4.6×10-5 m3/h,进料浓度为4.73×10-3 kmol/m3时,研究传质系数对吸附穿透曲线的影响,如图 8所示.

图 8 不同传质系数下锌离子吸附穿透曲线 Fig. 8 Breakthrough curves of zinc ion at different mass transfer coefficient values

当传质系数KS从4.3×10-4 s-1增加至8.0×10-4 s-1时,开始穿透的时间从0.98 h增加到1.00 h,但传质系数KS为4.3×10-3 s-1、4.3×10-2 s-1穿透曲线几乎重合.数据表明,随着传质系数的增加,吸附床层穿透所需时间增加,但当传质系数KS≥4.3×10-3 s-1时,穿透曲线基本不改变.

当铜离子的进料流量为4.6×10-5 m3/h,进料浓度为9.44×10-5 kmol/m3时,研究传质系数对铜离子吸附的穿透曲线的影响,如图 9所示.

图 9 不同传质系数下铜离子吸附穿透曲线 Fig. 9 Breakthrough curves of copper ion at different mass transfer coefficient values

当传质系数KS从5.00×10-4 s-1增大到1.00×10-3 s-1时,两者曲线相近.当传质系数KS从1.00×10-3 s-1增大到5.00×10-3 s-1时,开始穿透时间从1.36 h延迟到1.61 h.但传质系数KS分别为0.01 s-1、0.10 s-1、1.00 s-1时,穿透曲线基本重合.这些表明,传质系数KS≤0.01 s-1时,穿透时间随着传质系数的减小而增加.但当传质系数KS≥0.01 s-1时,穿透曲线基本保持改变.

3.5 吸附柱的轴向荷载分布

当固体床进料流量为4.6×10-5 m3/h,进料锌离子浓度为4.73×10-3 kmol/m3时,进料锌离子浓度为9.44×10-3 kmol/m3时,以吸附床层高度为横坐标,以与床层高度对应的溶液浓度为纵坐标,绘制出不同时间点的溶液浓度与床层高度的变化关系曲线图,得到吸附柱轴向荷载分布曲线,如图 10图 11所示.

图 10 锌离子吸附荷载轴向分布 Fig. 10 Axial distributions of zinc ion adsorption capacity

图 11 铜离子吸附荷载轴向分布 Fig. 11 Axial distributions of copper ion adsorption capacity

吸附初期,吸附床层中离子被吸附,随着吸附时间的增加,吸附床层的进料口端开始形成传质区,并向床层出料口端移动.当吸附持续时间达t(Zn2+)为1.00 h、t(Cu2+)为1.36 h时,传质区的下表面到达床层底部,此时吸附床层开始被穿透;当吸附持续时间达t(Zn2+)为12.5 h、t(Cu2+)为6.00 h时,传质区的上表面到达床层底部,此时床层已被完全穿透,吸附柱丧失了吸附能力,床层出料口离子浓度与进料口浓度相等.此外,传质区在吸附过程沿轴向床层底部移动的过程中,传质区的长度并没有发生变化.

4 结论

本研究在液相吸附基本理论的基础上,对树脂固定吸附柱床层的吸附过程做出假设:无轴向扩散、无温度变化、液相为扩散活塞流等,建立物料平衡方程、吸附速率方程和吸附等温方程.以一阶迎风差分法为离散方法,Aspen Adsorption软件作为计算工具,研究进料浓度、流量和传质系数对吸附过程的影响,结论如下:

1)运用Aspen Adsorption软件对树脂吸附锌离子过程进行数值模拟,并将模拟数据与实验数据相对比,验证了所建数学模型的可靠性;

2)在吸附初期,吸附柱出料口中含有少量铜、锌离子,随着吸附的进行,由于树脂量是固定的,吸附床层逐渐接近饱和,吸附能力逐渐降低;

3)在进料流量、进料铜、锌离子浓度增加时,穿透曲线均左移,但进料流量的作用更为突出.建议通过提高进料流量以提高树脂吸附铜、锌离子效率;

4)在传质系数大于临界值时,吸附穿透曲线基本重合,其对吸附性能并没有显著的影响.临界KS(Zn2+)=4.3×10-3 s-1, KS(Cu2+)=1.00×10-2 s-1.

参考文献
[1]
邱廷省, 郝志伟, 成先雄. 含氰废水处理技术评述与展望[J]. 江西冶金, 2002, 22(3): 25–29.
[2]
兰新哲, 宋永辉, 廖赞, 等. 含氰尾液综合回收研究[J]. 稀有金属, 2005, 29(4): 493–497. DOI: 10.3969/j.issn.0258-7076.2005.04.027.
[3]
XIE F, DREISINGER D, DOYLE F. A review on recovery of copper and cyanide from waste cyanide solutions[J]. Mineral Processing & Extractive Metallurgy Review, 2013, 34(6): 387–411.
[4]
ZHANG Y L, FANG T, YU X J. Adsorption of zinc and cyanide from cyanide effluents on anionic ion-exchange resin[J]. Chemical Research in Chinese Universities, 2013, 29(1): 144–149. DOI: 10.1007/s40242-013-2070-6.
[5]
刘旸, 刘静欣, 江晓健, 等. 废弃电路板中非金属组分的回收利用[J]. 有色金属科学与工程, 2016, 7(2): 1–7.
[6]
何敏, 兰新哲, 朱国才, 等. 离子交换树脂处理含氰废水进展[J]. 黄金, 2006, 27(1): 45–48.
[7]
黄爱华. 提金含氰废水处理工艺研究现状及发展趋势分析[J]. 黄金科学技术, 2014, 22(2): 83–89.
[8]
王碧侠, 兰新哲, 宋永辉, 等. 用离子交换树脂从含氰溶液中回收有价成分的研究[J]. 有色金属工程, 2002, 54(7): 124–127.
[9]
吴永胜. 金置换率下降原因分析及对策[J]. 有色金属科学与工程, 2011, 2(5): 36–40.
[10]
何敏, 兰新哲, 朱国才. D296树脂对氰化物的吸附和解吸性能研究[J]. 金属矿山, 2005(10): 56–59. DOI: 10.3321/j.issn:1001-1250.2005.10.016.
[11]
张笠, 李慕玲, 迟杰, 等. 树脂循环复用处理镀金废液的研究[J]. 离子交换与吸附, 1997, 13(6): 583–589.
[12]
SHU Z N, XIONG C H, WANG X. Adsorption behavior and mechanism of amino methylene phosphonic acid resin for Ag(I)[J]. Transactions of Nonferrous Metals Society of China, 2006, 16(3): 700–704. DOI: 10.1016/S1003-6326(06)60124-6.
[13]
TORRE M, BACHILLER D, RENDUELES M. Cyanide recovery from gold extraction process waste effluents by ion exchange I. Equilibrium and kinetics[J]. Solvent Extraction & Ion Exchange, 2006, 24(24): 99–117.
[14]
DASH R R, DASH R R, BALOMAJUMDAR C. Treatment of cyanide bearing effluents by adsorption, biodegradation and combined processes:effect of process parameters[J]. Desalination & Water Treatment, 2014, 52(16/17/18): 3355–3366.
[15]
彭松水. 基于Aspen Adsorption的CO2净化工艺动态模拟[J]. 石化技术, 2015, 22(5): 13–16. DOI: 10.3969/j.issn.1006-0235.2015.05.008.
[16]
李建伟, 张三莉, 宣自润, 等. 煤层气单塔吸附过程的模拟[J]. 天然气化工:c1化学与化工, 2016, 41(1): 54–57.
[17]
杨中山, 刘庭耀, 王宝, 等. 活性炭吸附提取氰化亚金离子的数学模拟[J]. 过程工程学报, 2015, 15(4): 594–598.
[18]
刘金昌, 曹俊雅, 刘娟, 等. 基于Aspen Adsorption的活性炭吸附溶液中镍离子的模拟研究[J]. 计算机与应用化学, 2013, 30(4): 439–442.
[19]
杨金垚, 方梦祥, 岑旗钢, 等. 活性炭捕集燃煤烟气中二氧化碳的模拟分析[J]. 浙江大学学报(工学版), 2018, 52(11): 2142–2149. DOI: 10.3785/j.issn.1008-973X.2018.11.013.
[20]
薛才红, 郑娆, 赵鹏翔, 等. 基于Aspen Adsorption的乙醇变压吸附脱水工艺模拟[J]. 石油化工, 2016, 45(9): 1107–1111. DOI: 10.3969/j.issn.1000-8144.2016.09.014.
[21]
竹涛, 刘恩彤, 王艳霞, 等. 煤矿乏风瓦斯变压吸附技术操作条件优化试验研究[J]. 矿业科学学报, 2016, 1(2): 196–202.
[22]
李立清, 孙政, 邢俊冬, 等. 固体床吸附过程传热传质耦合影响数值模拟[J]. 化工学报, 2010, 61(11): 2810–2816.
[23]
马立成.铜离子印迹磁性壳聚糖微球的制备和选择吸附性能研究[D].北京: 北京有色金属研究总院, 2015: 40-65. http://cdmd.cnki.com.cn/Article/CDMD-86403-1015639006.htm