堆浸工艺中流动反应变形传质耦合过程数值模拟.pdfVIP

堆浸工艺中流动一反应一变形 —传质耦合过程数值模拟 刘金枝1 吴爱祥2苏永定2 (1．中南大学数学科学与计算技术学院，湖南长沙。410075； 2．北京科技大学土木与环境工程学院，北京，100083) 摘要：本文基于堆浸过程中的孔隙变形，发展了一孔隙发育模型，并建立了堆浸工艺中流动．反应．变形．传质全 耦合模型控制方程组。在喷淋强度不变和常水头的情况下，给出了方程的数值解，数值计算采用有限特征差分 方法。数值结果描述了浸出过程中孔隙率的变化以及堆中目的金属离子浓度的分布。其中孔隙率呈上升趋势， 浸出离子浓度在顶部、中部和底部的分布情形与实际物理行为相一致。 关键词：堆浸孔隙发育传质数值解 1 引言 溶浸采矿是通过注入某些化学溶剂，以溶解矿石或矿床中有用成分的采矿方法，一般采用原地浸出和堆 浸两种形式。其中堆浸技术以其工艺简单、投资少、见效快、管理简单等诸多优点而被广泛应用于从低品位 矿石中提取有用金属[1~31。堆浸过程中，由于目的金属与注入溶浸剂之间的化学或生化反应，以及目的金属的 溶解／沉淀，导致多孑L介质的机械变形【4硼。 近年来，对堆浸工艺的理论研究越来越引起人们的广泛关注[7。10】，值得指出的是已有的研究工作都没有 考虑多孔介质的机械变形。本文将在已有研究成果的基础上，考虑堆浸过程中的变形因素，发展一孔隙发育 模型，并建立堆浸过程中流动一反应一变形一传质全耦合模型控制方程组。在喷淋强度不变和常水头的情况 下，给出了一个一维情形的数值算例，数值结果给出了浸出过程中孔隙率的分布、浸出矿石浓度的分布以及 矿堆中溶浸液饱和度的分布。 2流动一反应一变形一传质全耦合模型控制方程 2．1孔隙发育模型 假设矿石颗粒中的反应锋面随时间而向内部推进，反应区中的有 用金属全部浸出，表征单元体(REV)／111内的矿岩散体介质骨架不会因 浸出而垮塌(见图1)，且化学反应速度很快，浸出过程反应的速度为 扩散控制。由浸出反应收缩核模型，此时孔隙率的改变量可用下式表 示为： 4 ，4， △目=秒一岛=喜6％至二芋=喜Gk(·一号) cl， 式中：N为单元体中颗粒个数；a；为浸出前颗粒的半径；n为未反应 图1反应区示意图 核的半径；G行为单颗粒中目的金属的体积比；0为动态有效孔隙率。 Fi2．1Schematicofreactivezone 若颗粒粒径一致，反应均匀：口，-口，垆，_，2i=2，则方程(1)变形为： 第一作者简介：刘金枝，1971年生，女，副教授，博士，从事偏微分方程理论及应用以及冶金工程数值模拟方向的研究。 E-mail：liu0619@mail．CSU．edu．Cll ·500· △目=口一岛=G(·一事) (2) ^， 式中：q=∑瓯为整个单元体中目的金属所占体积比。 i=I 因矿物反应绝大多数受矿岩表面和溶浸剂浓度控制，因此可定义反应速率为与反应前后颗粒的表面积之 比和溶浸剂的浓度Cl的乘积成正比【12l： 昙(轷一后卧 式甲：k为比例常数。 结合式(2)和式(3)可得孔隙率的发育模型为： 署=慨(z一等卜 ㈣ 令够为所有有用金属全部浸出之后介质的孔隙率，显然，G卢够一oo。方程(4)变形为：