《工程科学学报》：硫化铜矿粒孔隙模型重构与溶液渗流模拟

工程科学学报 Chinese Journal of Engineering 硫化铜矿粒孔隙模型重构与溶液渗流棋拟 尹升华宋庆陈威陈勋 Pore model reconstruction of copper sulfide ore agglomerate and simulation of solution seepage YIN Sheng-hua.SONG Qing.CHEN Wei.CHEN Xun 引用本文： 尹升华，宋庆，陈威，陈勋.硫化铜矿粒孔隙模型重构与溶液渗流模拟[J.工程科学学报，2021,43(4：495-502.doi: 10.13374j.issn2095-9389.2020.02.27.002 YIN Sheng-hua,SONG Qing.CHEN Wei,CHEN Xun.Pore model reconstruction of copper sulfide ore agglomerate and simulation of solution seepage[J.Chinese Journal of Engineering,2021,43(4):495-502.doi:10.13374/j.issn2095-9389.2020.02.27.002 在线阅读View online::htps:/doi.org/10.13374.issn2095-9389.2020.02.27.002 您可能感兴趣的其他文章 Articles you may be interested in 次生硫化铜矿制粒试验 Agglomeration experiment of secondary copper sulfide ore 工程科学学报.2019,41(9：1127 https::/doi.org10.13374j.issn2095-9389.2019.09.003 次生硫化铜矿微生物浸出实验 Secondary copper sulfide bioleaching experiments 工程科学学报.2017,3910)：1498 https:/1doi.org10.13374.issn2095-9389.2017.10.006 我国铜矿微生物浸出技术的研究进展 Progress of research in copper bioleaching technology in China 工程科学学报.2019,41(2：143htps:/1doi.og/10.13374j.issn2095-9389.2019.02.001 低品位氧硫混合铜矿的酸性制粒及机理 Acid agglomeration and mechanism analysis of a low-grade oxide-sulfide mixed copper ore 工程科学学报.2017,399%：1321 https:1oi.org10.13374.issn2095-9389.2017.09.004 电弧焊接数值模拟中热源模型的研究与发展 Research and development of a heat-source model in numerical simulations for the arc welding process 工程科学学报.2018.40(4：389 https:/doi.org10.13374.issn2095-9389.2018.04.001 连铸坯脱氢退火数值模拟 Numerical simulation of dehydrogenation annealing in bloom 工程科学学报.2020.42(7)：862htps/doi.org/10.13374j.issn2095-9389.2020.03.16.003

硫化铜矿粒孔隙模型重构与溶液渗流模拟 尹升华 宋庆 陈威 陈勋 Pore model reconstruction of copper sulfide ore agglomerate and simulation of solution seepage YIN Sheng-hua, SONG Qing, CHEN Wei, CHEN Xun 引用本文: 尹升华, 宋庆, 陈威, 陈勋. 硫化铜矿粒孔隙模型重构与溶液渗流模拟[J]. 工程科学学报, 2021, 43(4): 495-502. doi: 10.13374/j.issn2095-9389.2020.02.27.002 YIN Sheng-hua, SONG Qing, CHEN Wei, CHEN Xun. Pore model reconstruction of copper sulfide ore agglomerate and simulation of solution seepage[J]. Chinese Journal of Engineering, 2021, 43(4): 495-502. doi: 10.13374/j.issn2095-9389.2020.02.27.002 在线阅读 View online: https://doi.org/10.13374/j.issn2095-9389.2020.02.27.002 您可能感兴趣的其他文章 Articles you may be interested in 次生硫化铜矿制粒试验 Agglomeration experiment of secondary copper sulfide ore 工程科学学报. 2019, 41(9): 1127 https://doi.org/10.13374/j.issn2095-9389.2019.09.003 次生硫化铜矿微生物浸出实验 Secondary copper sulfide bioleaching experiments 工程科学学报. 2017, 39(10): 1498 https://doi.org/10.13374/j.issn2095-9389.2017.10.006 我国铜矿微生物浸出技术的研究进展 Progress of research in copper bioleaching technology in China 工程科学学报. 2019, 41(2): 143 https://doi.org/10.13374/j.issn2095-9389.2019.02.001 低品位氧硫混合铜矿的酸性制粒及机理 Acid agglomeration and mechanism analysis of a low-grade oxide-sulfide mixed copper ore 工程科学学报. 2017, 39(9): 1321 https://doi.org/10.13374/j.issn2095-9389.2017.09.004 电弧焊接数值模拟中热源模型的研究与发展 Research and development of a heat-source model in numerical simulations for the arc welding process 工程科学学报. 2018, 40(4): 389 https://doi.org/10.13374/j.issn2095-9389.2018.04.001 连铸坯脱氢退火数值模拟 Numerical simulation of dehydrogenation annealing in bloom 工程科学学报. 2020, 42(7): 862 https://doi.org/10.13374/j.issn2095-9389.2020.03.16.003

工程科学学报.第43卷，第4期：495-502.2021年4月 Chinese Journal of Engineering,Vol.43,No.4:495-502,April 2021 https://doi.org/10.13374/j.issn2095-9389.2020.02.27.002;http://cje.ustb.edu.cn 硫化铜矿粒孔隙模型重构与溶液渗流模拟 尹升华2)，宋庆2)四，陈威2，陈勋2 1)北京科技大学金属矿山高效开采与安全教育部重点实验室，北京1000832)北京科技大学土木与资源工程学院，北京100083 ☒通信作者，E-mail:SQ1213ueen@163.com 摘要以次生硫化铜矿粉为原料，添加黏结剂、氯化钠制备矿粒，并借助CT扫描技术、图像处理及三维重构方法，开展了 单个矿粒浸出试验，探究了溶浸前后矿粒内部的孔隙变化；运用COMSOL Multiphysics模拟仿真软件，构建了溶液在孔隙通 道中流动的仿真模型.结果表明：经过一周时间的溶浸，矿粒内部孔隙的数目、平均体积、平均表面积及孔隙平均等效直径分 别增长了99%、151%、223%和90%.孔隙率增长了4倍，孔隙连通度增长了近2倍.在孔隙通道较狭窄的区域和底部区域， 溶液的流速、压力急剧增加，对矿粒结构的稳定性产生较大影响. 关键词硫化铜矿旷：制粒技术；CT扫描技术；孔隙演变；数值模拟 分类号TD862 Pore model reconstruction of copper sulfide ore agglomerate and simulation of solution seepage YIN Sheng-hua2),SONG Qing2,CHEN Wei2),CHEN Xun2) 1)Key Laboratory of Ministry of Education for High-Efficient Mining and Safety of Metal,University of Science and Technology Beijing,Beijing 100083,China 2)School of Civil and Resource Engineering,University of Science and Technology Beijing,Beijing 100083,China Corresponding author,E-mail:SQ1213ueen@163.com ABSTRACT Bioleaching technology,a method used in mining,utilizes organisms or microorganisms to dissolve and leach useful components from ore.This technology is applied widely in many countries in mining and production of various metals.Although considerable economic benefits have been obtained with the help of bioleaching technology,this technology faces many challenges.For example,as the leaching process progresses,the permeability of the ore heap is worsened from the presence of fine particles and muddy ore.As such,improving heap permeability is a key focus in the field of leaching mining technology.Agglomeration technology,an effective method for improving the permeability of ore heaps,involves crushing the ore and agglomerating it with binders before leaching,which alters the structure of the ore,especially ore powder.The resulting granular structure provides space for the flow of the solution and greatly improves the permeability of the leaching system.The pore structure is constantly changing as the leaching process progresses,which affects the seepage of the solution and even the effect of leaching the ore.In this study,leaching tests were conducted using single-ore granulation of secondary copper sulfide ore powder,binders,and NaCl.Changes in the pore structure during leaching were investigated using X-ray computerized tomography,digital image processing,and three-dimensional reconstruction methods.In addition,COMSOL Multiphysics simulation software was used to construct a simulation model of the flow of the solution in the pore channels.The results reveal that the number,average volume,average surface area,and average equivalent diameter of pores inside the 收稿日期：2020-02-27 基金项目：国家优秀青年科学基金资助项目(51722401)：国家自然科学基金重点资助项目(51734001)：中央高校基本科研业务费专项资金 资助项目(FRF-TP.18-003C1)

硫化铜矿粒孔隙模型重构与溶液渗流模拟 尹升华1,2)，宋    庆1,2) 苣，陈    威1,2)，陈    勋1,2) 1) 北京科技大学金属矿山高效开采与安全教育部重点实验室，北京 100083    2) 北京科技大学土木与资源工程学院，北京 100083 苣通信作者，E-mail：SQ1213ueen@163.com 摘    要    以次生硫化铜矿粉为原料，添加黏结剂、氯化钠制备矿粒，并借助 CT 扫描技术、图像处理及三维重构方法，开展了 单个矿粒浸出试验，探究了溶浸前后矿粒内部的孔隙变化；运用 COMSOL Multiphysics 模拟仿真软件，构建了溶液在孔隙通 道中流动的仿真模型. 结果表明：经过一周时间的溶浸，矿粒内部孔隙的数目、平均体积、平均表面积及孔隙平均等效直径分 别增长了 99%、151%、223% 和 90%，孔隙率增长了 4 倍，孔隙连通度增长了近 2 倍. 在孔隙通道较狭窄的区域和底部区域， 溶液的流速、压力急剧增加，对矿粒结构的稳定性产生较大影响. 关键词    硫化铜矿；制粒技术；CT 扫描技术；孔隙演变；数值模拟 分类号    TD862 Pore model reconstruction of copper sulfide ore agglomerate and simulation of solution seepage YIN Sheng-hua1,2) ，SONG Qing1,2) 苣 ，CHEN Wei1,2) ，CHEN Xun1,2) 1) Key Laboratory of Ministry of Education for High-Efficient Mining and Safety of Metal, University of Science and Technology Beijing, Beijing 100083, China 2) School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China 苣 Corresponding author, E-mail: SQ1213ueen@163.com ABSTRACT    Bioleaching  technology,  a  method  used  in  mining,  utilizes  organisms  or  microorganisms  to  dissolve  and  leach  useful components  from  ore.  This  technology  is  applied  widely  in  many  countries  in  mining  and  production  of  various  metals.  Although considerable economic benefits have been obtained with the help of bioleaching technology, this technology faces many challenges. For example, as the leaching process progresses, the permeability of the ore heap is worsened from the presence of fine particles and muddy ore.  As  such,  improving  heap  permeability  is  a  key  focus  in  the  field  of  leaching  mining  technology.  Agglomeration  technology,  an effective  method  for  improving  the  permeability  of  ore  heaps,  involves  crushing  the  ore  and  agglomerating  it  with  binders  before leaching, which alters the structure of the ore, especially ore powder. The resulting granular structure provides space for the flow of the solution and greatly improves the permeability of the leaching system. The pore structure is constantly changing as the leaching process progresses, which affects the seepage of the solution and even the effect of leaching the ore. In this study, leaching tests were conducted using single-ore granulation of secondary copper sulfide ore powder, binders, and NaCl. Changes in the pore structure during leaching were  investigated  using  X-ray  computerized  tomography,  digital  image  processing,  and  three-dimensional  reconstruction  methods.  In addition, COMSOL Multiphysics simulation software was used to construct a simulation model of the flow of the solution in the pore channels. The results reveal that the number, average volume, average surface area, and average equivalent diameter of pores inside the 收稿日期: 2020−02−27 基金项目: 国家优秀青年科学基金资助项目（51722401）；国家自然科学基金重点资助项目（51734001）；中央高校基本科研业务费专项资金 资助项目（FRF-TP-18-003C1） 工程科学学报，第 43 卷，第 4 期：495−502，2021 年 4 月 Chinese Journal of Engineering, Vol. 43, No. 4: 495−502, April 2021 https://doi.org/10.13374/j.issn2095-9389.2020.02.27.002; http://cje.ustb.edu.cn

496 工程科学学报，第43卷，第4期 granulated ore increase respectively by 99%,151%,223%,and 90%.In addition,the porosity increases fourfold and the pore connectivity nearly triples.The flow velocity and pressure of the solution was found to increase sharply as the solution flows through the narrow and bottom areas of the pore channels,which has a great impact on the stability of the ore granulation structure. KEY WORDS copper sulfide ore;agglomeration technology;CT technology;pore evolution;numerical simulation 生物浸矿系统是一个固、液、气三相共存的体 的三维数值模型，探究渗流特性的变化，为进一步 系，溶液与矿石接触并发生化学反应，金属铜从固 了解制粒后的孔隙演变、探索矿粒内部溶液的渗 态矿石转移至液态溶液中，浸出液经收集、萃取等 流规律等提供参考 工序达到“采铜”的目的-)生物浸矿技术因其高 1试验概述 效率、低成本、绿色等优点，被广泛应用于多个国 家的金属矿开采-然而，在创造了可观的经济 1.1试验材料 效益的同时，生物浸铜技术面临着重大挑战，如浸 本试验所用矿样取自福建省某铜矿，其主要 出过程中无法保证堆浸系统的渗透性、浸矿菌种 成分以及铜物相分析结果如表1、表2所示.铜的 活性随浸出过程逐渐变差和浸矿体系基础理论不 品位为0.70%，铜的主要存在形态是次生硫化铜， 完善等问题？).浸矿系统中的孔隙为溶液流动提 其质量分数为0.60%，其次为原生硫化铜和自由氧 供了空间与路径，随着反应的发生，孔隙结构也不 化铜，结合铜的含量最少.预先将矿样破碎以方便 断发生变化，影响溶液的渗流特性，并最终影响浸 制粒工序，破碎后矿样粒径均小于5mm,90%以 矿效果0四制粒预处理技术的出现，较好地解决 上矿样粒径小于0.8mm,矿样粒径分布图如图1 了浸矿系统渗透性差的问题，改善了矿堆的渗透 所示 性：此外，开采及破碎过程中产生的粉矿也得到 表1矿样主要元素质量分数 了回收利用，减少了资源的流失与浪费4 Table 1 Mass fractions of major elements in mineral samples % 浸矿过程中，孔隙演变与溶液渗流相互作用， Elements Cu Fe S CaO Mgo AlO SiO, 相互影响.近年来，国内外专家开展了大量有关于 Mass fraction 0.701.671.100.300.045.1991.00 溶液渗流特性的研究，随着CT技术的广泛应用， 浸矿系统中的孔隙演变及溶液渗流规律的研究已 表2铜物相分析结果（质量分数） 取得了较好的成果-刘超研究了酸浸氧化铜 Table 2 Cu phase analysis results of mineral samples 矿石颗粒内部裂隙结构的演化规律，结果表明酸 Phase Mass fraction/% 浸作用下，矿粒孔隙率、裂隙尺寸和数量分别增加 Copper oxide 0.04 了约40%、3.5%和25%；氧化铜矿岩散体线孔隙率 Primary copper sulfide 0.05 也随时间呈现出一定的演化规律P.吴爱祥等以 Secondary copper sulfide 0.60 高泥氧化铜矿为试验材料，开展矿岩散体的柱浸 Combined copper 0.01 试验，结果表明松散颗粒在浸出过程中发生迁移， Total 0.70 浸出后平均孔隙率降低了2.63%，浸出5h后渗透 100 系数降低了37.5%2.L0等利用XMT技术，开展 -The line of mineral powder 了腐泥土(SAP)和针铁矿(G)两种镍矿石制粒后 80 柱浸试验，结果表明两组矿柱的孔隙度变化与质 量损失和坍落度变化一致，但前者有更好的孔隙 率与浸出稳定性四).尽管基于CT扫描技术的运 40 用，浸矿系统内部孔隙演变规律的研究已取得了 较好的成绩，但矿石制粒后矿粒内部孔隙演变的 研究仍较为空白 0.074 0.0980.125020.450.81.2525 本文利用CT扫描技术和数字图像处理技术， Ore particle size/mm 分析溶浸后矿粒内部的孔隙演变，利用COMSOL 图1矿样粒径分布图 Multiphysics模拟仿真软件构建矿粒内部溶液渗流 Fig.I Particle size distribution of ore samples

granulated  ore  increase  respectively  by  99%,  151%,  223%,  and  90%.  In  addition,  the  porosity  increases  fourfold  and  the  pore connectivity nearly triples. The flow velocity and pressure of the solution was found to increase sharply as the solution flows through the narrow and bottom areas of the pore channels, which has a great impact on the stability of the ore granulation structure. KEY WORDS    copper sulfide ore；agglomeration technology；CT technology；pore evolution；numerical simulation 生物浸矿系统是一个固、液、气三相共存的体 系，溶液与矿石接触并发生化学反应，金属铜从固 态矿石转移至液态溶液中，浸出液经收集、萃取等 工序达到“采铜”的目的[1−3] . 生物浸矿技术因其高 效率、低成本、绿色等优点，被广泛应用于多个国 家的金属矿开采[4−6] . 然而，在创造了可观的经济 效益的同时，生物浸铜技术面临着重大挑战，如浸 出过程中无法保证堆浸系统的渗透性、浸矿菌种 活性随浸出过程逐渐变差和浸矿体系基础理论不 完善等问题[7−9] . 浸矿系统中的孔隙为溶液流动提 供了空间与路径，随着反应的发生，孔隙结构也不 断发生变化，影响溶液的渗流特性，并最终影响浸 矿效果[10−12] . 制粒预处理技术的出现，较好地解决 了浸矿系统渗透性差的问题，改善了矿堆的渗透 性[13] ；此外，开采及破碎过程中产生的粉矿也得到 了回收利用，减少了资源的流失与浪费[14−16] . 浸矿过程中，孔隙演变与溶液渗流相互作用， 相互影响. 近年来，国内外专家开展了大量有关于 溶液渗流特性的研究，随着 CT 技术的广泛应用， 浸矿系统中的孔隙演变及溶液渗流规律的研究已 取得了较好的成果[17−19] . 刘超研究了酸浸氧化铜 矿石颗粒内部裂隙结构的演化规律，结果表明酸 浸作用下，矿粒孔隙率、裂隙尺寸和数量分别增加 了约 40%、3.5% 和 25%；氧化铜矿岩散体线孔隙率 也随时间呈现出一定的演化规律[20] . 吴爱祥等以 高泥氧化铜矿为试验材料，开展矿岩散体的柱浸 试验，结果表明松散颗粒在浸出过程中发生迁移， 浸出后平均孔隙率降低了 2.63%，浸出 5 h 后渗透 系数降低了 37.5% [21] . Lo 等利用 XMT 技术，开展 了腐泥土（SAP）和针铁矿（G）两种镍矿石制粒后 柱浸试验，结果表明两组矿柱的孔隙度变化与质 量损失和坍落度变化一致，但前者有更好的孔隙 率与浸出稳定性[22] . 尽管基于 CT 扫描技术的运 用，浸矿系统内部孔隙演变规律的研究已取得了 较好的成绩，但矿石制粒后矿粒内部孔隙演变的 研究仍较为空白. 本文利用 CT 扫描技术和数字图像处理技术， 分析溶浸后矿粒内部的孔隙演变，利用 COMSOL Multiphysics 模拟仿真软件构建矿粒内部溶液渗流 的三维数值模型，探究渗流特性的变化，为进一步 了解制粒后的孔隙演变、探索矿粒内部溶液的渗 流规律等提供参考. 1    试验概述 1.1    试验材料 本试验所用矿样取自福建省某铜矿，其主要 成分以及铜物相分析结果如表 1、表 2 所示. 铜的 品位为 0.70%，铜的主要存在形态是次生硫化铜， 其质量分数为 0.60%，其次为原生硫化铜和自由氧 化铜，结合铜的含量最少. 预先将矿样破碎以方便 制粒工序，破碎后矿样粒径均小于 5 mm，90% 以 上矿样粒径小于 0.8 mm，矿样粒径分布图如图 1 所示. 表 1    矿样主要元素质量分数 Table 1    Mass fractions of major elements in mineral samples % Elements Cu Fe S CaO MgO Al2O3 SiO2 Mass fraction 0.70 1.67 1.10 0.30 0.04 5.19 91.00 表 2    铜物相分析结果（质量分数） Table 2    Cu phase analysis results of mineral samples Phase Mass fraction/% Copper oxide 0.04 Primary copper sulfide 0.05 Secondary copper sulfide 0.60 Combined copper 0.01 Total 0.70 0.074 0.098 0.125 0.2 0.45 0.8 1.25 2 5 0 20 40 60 80 100 Cumulative mass fractions of ore particles vsmaller than certain sizes/ % Ore particle size/mm The line of mineral powder 图 1    矿样粒径分布图 Fig.1    Particle size distribution of ore samples · 496 · 工程科学学报，第 43 卷，第 4 期

尹升华等：硫化铜矿粒孔隙模型重构与溶液渗流模拟 497. 1.2试验原理 本次试验预先在制粒材料中添加适量的 NaCL,拟用NaCI溶解于水的过程模拟浸矿过程， 即金属铜从固态矿粒转移至液态溶液中.在实际 的硫化铜矿石溶浸过程中，往往使用酸性溶剂进 行溶浸，硫化铜溶解在酸性溶液中，而矿石其他成 分无法溶解，经固液分离即可达到提取铜的目的： 在本试验中，将制备好的矿粒溶浸在水中时， 图2试验矿粒.(a)溶浸时：(b)晾干后 NaCI极易溶解在水中，而矿粉等材料是不溶于水 Fig.2 Photographs of mineral agglomerate used in the test:(a)during 的.以上两个“溶解”过程就矿粒结构变化而言是 leaching;(b)after drying 相同的 三维扫描成像系统(nano Voxel3502E)获取矿粒各 基于NaCI的溶解作用，矿粒内部会产生新的 切面扫描图像，试验保证溶浸前后两次CT扫描各 孔隙并发生一系列变化.一方面，附着在NaCl上 参数完全相同.扫描原图如图3所示，其中，黑色 的矿石颗粒随NaCI溶解离开矿粒，掉落下来；另 部分代表矿粒内部的孔隙，白色部分表示矿粒中 一方面，矿粒表面及孔隙周围的矿石颗粒受溶液 的金属成分，灰色部分表示矿粒中的矿石 渗流流动压力的作用发生掉落.在试验溶浸过程 中，由于固态NaCl溶解和矿石颗粒的掉落，矿粒 (a) (b) 内部结构及孔隙结构发生变化，溶液渗流受路径 及场所变化的影响发生变化，用以模拟硫化铜矿 石溶浸过程中发生的变化. 1.3试验方案 选择合适的制粒黏结剂与试验矿粉混合均 匀，并添加质量分数为30%的NaCL,在100g的制 粒材料中加入30g水进行制粒，后经干燥熟化等 图3扫描原图.(a)浸出前：(b)浸出后 工序制好矿粒.选择一个形态较好的矿粒进行第 Fig.3 Original CT scan images of mineral agglomerate:(a)before 一次CT扫描.将试验矿粒放置在塑料圆桶中，添 leaching.(b)after leaching 加去离子水直至没过矿粒，浸泡矿粒，7d后取出 由CT扫描原图对比可得，经过7d时间的溶 矿粒并晾干，将矿粒进行第二次CT扫描.基于此， 浸，矿粒的整体形态发生了变化，溶浸前矿粒呈较 得到矿粒两次CT扫描的图像 为规整的椭圆形，内部孔隙（黑色部分）较少，且分 布较为分散：溶浸后矿粒整体形态极其不规则，孔 2分析与讨论 隙数目增多，且面积增大.这说明经过溶浸，矿粒 2.1矿粒的制备 内部的孔隙明显增多.基于CT扫描得到的数据， 以硫化铜矿石为试验原料，开展制粒试验，探 可得溶浸前后矿粒3个方向上的尺寸，浸矿前为 究最佳黏结剂及最佳制粒条件.通过比较相同试 23.08mm×24.73mm×26.38mm.浸矿后为24.76mm× 验条件下矿粒平均粒径以及酸浸下发生破裂的时 24.76mm×27.41mm,矿粒的三维尺寸均增大了，即 间，综合考虑制粒黏结效果、黏结剂的易得性及成 矿粒发生了膨胀.矿粒发生膨胀是由以下3个方 本等因素，选择SFS-2作为本次试验的制粒黏结 面的变化引起的：(1)矿粒中的固态NaCI溶解在 剂，SFS-2黏结剂中包括：硅酸钠、生石灰、粉煤 水中，矿粒内孔隙增大、增多，导致矿粒体积增大； 灰及水泥：通过比较矿粒湿强度、矿粒抗压强度和 (2)由于NaCI的溶解，紧紧附着其上的矿粉颗粒 酸浸时间，确定制粒最佳条件为：黏结剂添加量为 掉落至水中，矿粒内部结构发生了较大变化，矿石 16g,熟化加酸量为25kgt,制粒喷水量占矿粉质 颗粒之间的排列分布发生了变化，导致整体矿粒 量的30%.基于以上制粒试验结果，制备浸出试验 膨胀：(3)黏结剂中的生石灰与水发生反应，生成 所用矿粒，如图2所示 了Ca(OH)2沉淀附着在矿粒上，导致矿粒体积增 2.2孔隙模型重构 大，整体效果表现为膨胀 本次试验，采用北京科技大学高分辨X射线 首先对CT扫描原图像进行图像预处理，以消

1.2    试验原理 本次试验预先在制粒材料中添加适量 的 NaCl，拟用 NaCl 溶解于水的过程模拟浸矿过程， 即金属铜从固态矿粒转移至液态溶液中. 在实际 的硫化铜矿石溶浸过程中，往往使用酸性溶剂进 行溶浸，硫化铜溶解在酸性溶液中，而矿石其他成 分无法溶解，经固液分离即可达到提取铜的目的； 在本试验中 ，将制备好的矿粒溶浸在水中时 ， NaCl 极易溶解在水中，而矿粉等材料是不溶于水 的. 以上两个“溶解”过程就矿粒结构变化而言是 相同的. 基于 NaCl 的溶解作用，矿粒内部会产生新的 孔隙并发生一系列变化. 一方面，附着在 NaCl 上 的矿石颗粒随 NaCl 溶解离开矿粒，掉落下来；另 一方面，矿粒表面及孔隙周围的矿石颗粒受溶液 渗流流动压力的作用发生掉落. 在试验溶浸过程 中，由于固态 NaCl 溶解和矿石颗粒的掉落，矿粒 内部结构及孔隙结构发生变化，溶液渗流受路径 及场所变化的影响发生变化，用以模拟硫化铜矿 石溶浸过程中发生的变化. 1.3    试验方案 选择合适的制粒黏结剂与试验矿粉混合均 匀，并添加质量分数为 30% 的 NaCl，在 100 g 的制 粒材料中加入 30 g 水进行制粒，后经干燥熟化等 工序制好矿粒. 选择一个形态较好的矿粒进行第 一次 CT 扫描. 将试验矿粒放置在塑料圆桶中，添 加去离子水直至没过矿粒，浸泡矿粒，7 d 后取出 矿粒并晾干，将矿粒进行第二次 CT 扫描. 基于此， 得到矿粒两次 CT 扫描的图像. 2    分析与讨论 2.1    矿粒的制备 以硫化铜矿石为试验原料，开展制粒试验，探 究最佳黏结剂及最佳制粒条件. 通过比较相同试 验条件下矿粒平均粒径以及酸浸下发生破裂的时 间，综合考虑制粒黏结效果、黏结剂的易得性及成 本等因素，选择 SFS–2 作为本次试验的制粒黏结 剂，SFS–2 黏结剂中包括：硅酸钠、生石灰、粉煤 灰及水泥；通过比较矿粒湿强度、矿粒抗压强度和 酸浸时间，确定制粒最佳条件为：黏结剂添加量为 16 g，熟化加酸量为 25 kg·t−1，制粒喷水量占矿粉质 量的 30%. 基于以上制粒试验结果，制备浸出试验 所用矿粒，如图 2 所示. 2.2    孔隙模型重构 本次试验，采用北京科技大学高分辨 X 射线 三维扫描成像系统（nano Voxel 3502E）获取矿粒各 切面扫描图像，试验保证溶浸前后两次 CT 扫描各 参数完全相同. 扫描原图如图 3 所示，其中，黑色 部分代表矿粒内部的孔隙，白色部分表示矿粒中 的金属成分，灰色部分表示矿粒中的矿石. 由 CT 扫描原图对比可得，经过 7 d 时间的溶 浸，矿粒的整体形态发生了变化，溶浸前矿粒呈较 为规整的椭圆形，内部孔隙（黑色部分）较少，且分 布较为分散；溶浸后矿粒整体形态极其不规则，孔 隙数目增多，且面积增大. 这说明经过溶浸，矿粒 内部的孔隙明显增多. 基于 CT 扫描得到的数据， 可得溶浸前后矿粒 3 个方向上的尺寸，浸矿前为 23.08 mm×24.73 mm×26.38 mm，浸矿后为 24.76 mm× 24.76 mm×27.41 mm，矿粒的三维尺寸均增大了，即 矿粒发生了膨胀. 矿粒发生膨胀是由以下 3 个方 面的变化引起的：（1）矿粒中的固态 NaCl 溶解在 水中，矿粒内孔隙增大、增多，导致矿粒体积增大； （2）由于 NaCl 的溶解，紧紧附着其上的矿粉颗粒 掉落至水中，矿粒内部结构发生了较大变化，矿石 颗粒之间的排列分布发生了变化，导致整体矿粒 膨胀；（3）黏结剂中的生石灰与水发生反应，生成 了 Ca(OH)2 沉淀附着在矿粒上，导致矿粒体积增 大，整体效果表现为膨胀. 首先对 CT 扫描原图像进行图像预处理，以消 (a) (b) 图 2    试验矿粒. （a）溶浸时；（b）晾干后 Fig.2    Photographs of mineral agglomerate used in the test: (a) during leaching; (b) after drying (a) (b) 图 3    扫描原图. （a）浸出前；（b）浸出后 Fig.3     Original  CT  scan  images  of  mineral  agglomerate:  (a)  before leaching; (b) after leaching 尹升华等： 硫化铜矿粒孔隙模型重构与溶液渗流模拟 · 497 ·

498 工程科学学报，第43卷，第4期 除因噪点、模糊等对后续实验分析造成的干扰 连通度为贯穿上界面和下界面的孔隙通道的体积 对CT扫描图像进行预处理包括降噪滤波、阈值调 占孔隙总体积的比值).在生物浸矿过程中，孔隙 整、孔隙分割及矿粒切割等，CT图像预处理流程 连通度越大，溶液越容易进入矿粒内部流动，即生 如图4所示，得到便于试验分析的700像素×700 物与矿石的接触越充分.基于矿粒的三维重构结 像素×700像素的立方体矿块图像数据，实际尺寸 构，溶浸前矿块在X、Y、Z三个方向上均没有连通 为1.2cm×1.2cm×1.2cm. 上界面和下界面的孔隙通道，无法直接进行孔隙 连通度的计算.假设矿粒内部最大的孔隙是小孔 隙相互连通贯穿的结果.基于此，通过矿粒中最大 孔隙的体积可确定矿粒的孔隙连通度，如公式(2) 所示.孔隙率及孔隙连通度计算结果如表4所示. Filter and Adjust Cut 3D de-nosing thresholds Segment pores structure 9=7×100% (1) 0=x100% (2) V 其中，为孔隙率；θ为孔隙连通度；'1为矿粒内部 图4图像预处理流程 孔隙总体积，m3;V为矿粒总体积，um3;Vo为矿粒 Fig.4 Workflow of CTscan image preprocessing 中最大的孔隙的体积，um3. 浸前矿粒的孔隙率为3.2%，浸后矿粒的孔隙 由于孔隙形态不规则，无法从数据的角度定 率达到16%，说明溶浸过程中，矿粒内部孔隙增 量对比分析，将孔隙近似等效为球体，基于三维重 多，矿粒内部的结构也发生了明显改变.矿粒内部 构的矿粒结构，得到矿粒内部孔隙参数数据，并计 孔隙增多的原因有：(1)矿粒内部的固态NaCI溶 算其增长率，结果如表3所示.对比浸前与浸后矿 解在水中，导致矿粒结构改变、孔隙增多；(2)随 粒的特征参数，矿粒内部孔隙数目增长了99%，孔 着NaCI的溶解，附着在NaCI上的矿石颗粒随之 隙平均体积增长了151%，平均表面积增长了 掉落至溶液中，亦造成了矿粒结构改变，孔隙增 223%,平均等效直径增长了90%.溶浸过程后矿粒 多，如图5所示. 的孔隙数目、孔隙平均体积、孔隙平均面积以及 溶浸前矿粒的孔隙连通度为29.96%，而溶浸 孔隙等效直径均明显增加 后孔隙连通度达到了88.26%，连通度增加了195%， 孔隙率是指多孔物质中孔隙体积占物质总体 说明在模拟溶浸过程中，矿粒内部的部分孔隙相 积的百分数，是表征矿粒结构及孔隙特征的一个 互连通，形成了孔隙通道.孔隙通道的形成为溶液 重要物理量.浸矿系统的孔隙率与系统渗透性联 在矿粒内部流动提供了空间和通路，基于此，矿粒 系密切：孔隙率越大，系统的渗透性越好.基于以 与溶液接触更加充分，反应也更加充分，从而促进 上数据，根据公式(1)计算矿粒内部孔隙率.孔隙 了新孔隙的生成与连通 表3矿粒孔隙参数变化 Table 3 Variation of pore parameters Parameter Pore number Mean pore volume/(10%um') Average pore surface area/I0°μm) Average equivalent pore diameter/um Before leaching 288 1.71 1.16 304 After leaching 575 4.29 3.75 579 Growth rate/% 99 151 223 90 表4矿粒孔隙率及孔隙连通度变化 Table 4 Evolution of porosity and pore connectivity Parameter Porosity/% Growth rate/% Porosity connectivity/% Growth rate/% Before leaching 3.20 29.96 400 195 After leaching 16.00 88.26

除因噪点、模糊等对后续实验分析造成的干扰. 对 CT 扫描图像进行预处理包括降噪滤波、阈值调 整、孔隙分割及矿粒切割等，CT 图像预处理流程 如图 4 所示，得到便于试验分析的 700 像素×700 像素×700 像素的立方体矿块图像数据,实际尺寸 为 1.2 cm×1.2 cm×1.2 cm. 由于孔隙形态不规则，无法从数据的角度定 量对比分析，将孔隙近似等效为球体，基于三维重 构的矿粒结构，得到矿粒内部孔隙参数数据，并计 算其增长率，结果如表 3 所示. 对比浸前与浸后矿 粒的特征参数，矿粒内部孔隙数目增长了 99%，孔 隙平均体积增长 了 151%，平均表面积增长 了 223%，平均等效直径增长了 90%. 溶浸过程后矿粒 的孔隙数目、孔隙平均体积、孔隙平均面积以及 孔隙等效直径均明显增加. 孔隙率是指多孔物质中孔隙体积占物质总体 积的百分数，是表征矿粒结构及孔隙特征的一个 重要物理量. 浸矿系统的孔隙率与系统渗透性联 系密切；孔隙率越大，系统的渗透性越好. 基于以 上数据，根据公式（1）计算矿粒内部孔隙率. 孔隙 连通度为贯穿上界面和下界面的孔隙通道的体积 占孔隙总体积的比值[23] . 在生物浸矿过程中，孔隙 连通度越大，溶液越容易进入矿粒内部流动，即生 物与矿石的接触越充分. 基于矿粒的三维重构结 构，溶浸前矿块在 X、Y、Z 三个方向上均没有连通 上界面和下界面的孔隙通道，无法直接进行孔隙 连通度的计算. 假设矿粒内部最大的孔隙是小孔 隙相互连通贯穿的结果. 基于此，通过矿粒中最大 孔隙的体积可确定矿粒的孔隙连通度，如公式（2） 所示. 孔隙率及孔隙连通度计算结果如表 4 所示. φ = V1 V ×100% （1） θ = V0 V1 ×100% （2） 其中， φ 为孔隙率；θ 为孔隙连通度；V1 为矿粒内部 孔隙总体积，μm3 ；V 为矿粒总体积，μm3 ；V0 为矿粒 中最大的孔隙的体积，μm3 . 浸前矿粒的孔隙率为 3.2%，浸后矿粒的孔隙 率达到 16%，说明溶浸过程中，矿粒内部孔隙增 多，矿粒内部的结构也发生了明显改变. 矿粒内部 孔隙增多的原因有：（1）矿粒内部的固态 NaCl 溶 解在水中，导致矿粒结构改变、孔隙增多；（2）随 着 NaCl 的溶解，附着在 NaCl 上的矿石颗粒随之 掉落至溶液中，亦造成了矿粒结构改变，孔隙增 多，如图 5 所示. 溶浸前矿粒的孔隙连通度为 29.96%，而溶浸 后孔隙连通度达到了 88.26%，连通度增加了 195%， 说明在模拟溶浸过程中，矿粒内部的部分孔隙相 互连通，形成了孔隙通道. 孔隙通道的形成为溶液 在矿粒内部流动提供了空间和通路，基于此，矿粒 与溶液接触更加充分，反应也更加充分，从而促进 了新孔隙的生成与连通. 表 3 矿粒孔隙参数变化 Table 3 Variation of pore parameters Parameter Pore number Mean pore volume/(108 μm3 ) Average pore surface area/(106 μm2 ) Average equivalent pore diameter/μm Before leaching 288 1.71 1.16 304 After leaching 575 4.29 3.75 579 Growth rate/% 99 151 223 90 表 4 矿粒孔隙率及孔隙连通度变化 Table 4 Evolution of porosity and pore connectivity Parameter Porosity/% Growth rate/% Porosity connectivity/% Growth rate/% Before leaching 3.20 400 29.96 195 After leaching 16.00 88.26 Filter and de-nosing Adjust thresholds Segment pores Cut 3D structure 图 4    图像预处理流程 Fig.4    Workflow of CT scan image preprocessing · 498 · 工程科学学报，第 43 卷，第 4 期

尹升华等：硫化铜矿粒孔隙模型重构与溶液渗流模拟 499. 流动→单向流→层流物理场，进行稳态研究 pu.V)u=pl+μ(u+(Ww】]+F (3) p7.(w)=0 其中，p为溶液的密度，kgm3;u为溶液的流速， ms;p为压强，Pa;I为单位矩阵；μ为溶液的动力 黏度，Pas:F为流体所受外力，N 2.3.3设置参数及边界条件 由数据确定坐标系长度单位为微米，导入网 格数据.添加材料为水，设置密度为1100kgm, 动力黏度为0.9Pas,具体参数设置如表5所示.设 置研究对象壁条件为无滑移，如图7设置入口及 图5矿粉掉落 出口位置.设置边界条件，入口初始压力为0.715Pa, Fig.5 Photograph of ore powder that dropped off the agglomerate 出口压力为0Pa 2.3溶液渗流数值模拟 基于CT扫描的图像及三维重构的结果，在 表5构建模型的关键参数 溶浸后矿粒构建的三维图像中，裁取200像素× Table 5 Key parameters used in the model 200像素×200像素大小的立方体，实际尺寸为 Parameters Symbol Value 3.4mm×3.4mm×3.4mm.基于此数据，运用COMSOL Density/(kg-m) 1100 298.13 Multiphysics仿真软件构建孔隙通道模型，进行溶 Temperature/K T Dynamic viscosity/Pa's 0.9 液在孔隙通道中流动的数值模拟，分析溶液的流 Initial pressure/Pa Po 动特征 0.715 23.1生成网格数据 2.3.4溶液流动特性 基于CT扫描图像，三维重构孔隙结构，以孔 基于构建的数值模型，对孔隙通道中水的流速 隙中的孔隙通道为基础，构建溶液流动模型.基 与压力进行计算与分析，探究孔隙通道中溶液的 于CT扫描三维重构的孔隙形态极其不规整，无法 流动特性及变化规律包括流速与压力两个方面. 在COMSOL Multiphysics软件中直接构建形态相 图8显示的是孔隙通道中水的速度流线图像， 当的几何体，对孔隙通道数据进行一系列处理，包 即水的运动轨迹，图8(a)中流线间隔为0.1μm, 括移除不连通的小孔隙、生成表面以及划分体网 图8(b)中流线间隔为0.01um.由速度流线图可 格，如图6所示，获取网格数据用以构建模型24 得，通道中水的运动轨迹呈现出自上而下贯穿，整 体上先聚集后分散的形态：且部分区域流线分布 密集，部分区域流线分布稀疏.当孔隙通道中某处 的流线分布较为密集时，说明溶液更倾向于经过 Pore structure Pore channel Generate surface Partition grid 此处流动：相反，当某处的流线较为稀疏时，说明 图6生成网格数据流程 孔隙中的溶液不易流经此处，即此处更有可能会 Fig.6 Workflow for generating grid data 形成溶液空区，溶液无法流经此处.溶液在孔隙通 2.3.2基本假设 道中的流动存在优先流与空区的现象，优先流与 基于CT图像构建模型前，需要进行以下假 空区的存在都会导致溶液在孔隙通道中分布不均 设：溶液流动仅发生在孔隙通道中，未渗透到矿石 匀，溶液与矿石接触不充分、不彻底) 结构中：溶液流动过程中，孔隙通道不再发生变 分别在灯、YZ、XZ三个平面上对孔隙通道做 化；孔隙中的溶液是不可压缩、连续的.基于以上 切面处理，设置切面数为5，得到流体速度切面图， 假设，并且网格数据是基于孔隙通道建立的，即网 如图9所示，X、Y、Z三个方向的最大速度分别为 格显示的就是溶液在矿粒内部流动的空间与路 7.35×107、4.95×10-7和4.75×10-7ms、.三个方向 径，溶液的流动可以看作是在孔隙通道中做单向 上的最大速度均出现在孔隙通道底部、较为狭窄 层流运动，运动方程如式(3).在软件中选择流体 的位置，且同一切面图上，由四周向中心，速度逐

2.3    溶液渗流数值模拟 基于 CT 扫描的图像及三维重构的结果，在 溶浸后矿粒构建的三维图像中，裁取 200 像素× 200 像素×200 像素大小的立方体 ，实际尺寸为 3.4 mm×3.4 mm×3.4 mm. 基于此数据，运用COMSOL Multiphysics 仿真软件构建孔隙通道模型，进行溶 液在孔隙通道中流动的数值模拟，分析溶液的流 动特征. 2.3.1    生成网格数据 基于 CT 扫描图像，三维重构孔隙结构，以孔 隙中的孔隙通道为基础，构建溶液流动模型. 基 于 CT 扫描三维重构的孔隙形态极其不规整，无法 在 COMSOL Multiphysics 软件中直接构建形态相 当的几何体，对孔隙通道数据进行一系列处理，包 括移除不连通的小孔隙、生成表面以及划分体网 格，如图 6 所示，获取网格数据用以构建模型[24] . 2.3.2    基本假设 基于 CT 图像构建模型前，需要进行以下假 设：溶液流动仅发生在孔隙通道中，未渗透到矿石 结构中；溶液流动过程中，孔隙通道不再发生变 化；孔隙中的溶液是不可压缩、连续的. 基于以上 假设，并且网格数据是基于孔隙通道建立的，即网 格显示的就是溶液在矿粒内部流动的空间与路 径，溶液的流动可以看作是在孔隙通道中做单向 层流运动，运动方程如式（3）. 在软件中选择流体 流动→单向流→层流物理场，进行稳态研究. ρ(u · ∇)u = ∇ [ −pI+µ ( ∇u+(∇u) T )]+ F ρ∇ ·(u) = 0 （3） 其中， ρ 为溶液的密度， kg·m−3 ； u 为溶液的流速， m·s−1 ；p 为压强，Pa；I 为单位矩阵；μ 为溶液的动力 黏度，Pa·s；F 为流体所受外力，N. 2.3.3    设置参数及边界条件 由数据确定坐标系长度单位为微米，导入网 格数据. 添加材料为水，设置密度为 1100 kg·m−3 ， 动力黏度为 0.9 Pa·s，具体参数设置如表 5 所示. 设 置研究对象壁条件为无滑移，如图 7 设置入口及 出口位置. 设置边界条件，入口初始压力为 0.715 Pa， 出口压力为 0 Pa. 2.3.4    溶液流动特性 基于构建的数值模型，对孔隙通道中水的流速 与压力进行计算与分析，探究孔隙通道中溶液的 流动特性及变化规律,包括流速与压力两个方面. 图 8 显示的是孔隙通道中水的速度流线图像， 即水的运动轨迹，图 8（ a）中流线间隔为 0.1 μm， 图 8（b）中流线间隔为 0.01 μm. 由速度流线图可 得，通道中水的运动轨迹呈现出自上而下贯穿，整 体上先聚集后分散的形态；且部分区域流线分布 密集，部分区域流线分布稀疏. 当孔隙通道中某处 的流线分布较为密集时，说明溶液更倾向于经过 此处流动；相反，当某处的流线较为稀疏时，说明 孔隙中的溶液不易流经此处，即此处更有可能会 形成溶液空区，溶液无法流经此处. 溶液在孔隙通 道中的流动存在优先流与空区的现象，优先流与 空区的存在都会导致溶液在孔隙通道中分布不均 匀，溶液与矿石接触不充分、不彻底[25] . 分别在 XY、YZ、XZ 三个平面上对孔隙通道做 切面处理，设置切面数为 5，得到流体速度切面图， 如图 9 所示，X、Y、Z 三个方向的最大速度分别为 7.35×10−7、4.95×10−7 和 4.75×10−7 m·s−1、. 三个方向 上的最大速度均出现在孔隙通道底部、较为狭窄 的位置，且同一切面图上，由四周向中心，速度逐 表 5    构建模型的关键参数 Table 5    Key parameters used in the model Parameters Symbol Value Density/(kg·m−3) ρ 1100 Temperature/K T 298.13 Dynamic viscosity/Pa·s μ 0.9 Initial pressure/Pa P0 0.715 图 5    矿粉掉落 Fig.5    Photograph of ore powder that dropped off the agglomerate Pore structure Pore channel Generate surface Partition grid 图 6    生成网格数据流程 Fig.6    Workflow for generating grid data 尹升华等： 硫化铜矿粒孔隙模型重构与溶液渗流模拟 · 499 ·

500 工程科学学报，第43卷，第4期 (a) (b) 50300 5005s060060 (c)Xlum Y/um 750 100150 700 700 600 600 500 550 um/Z 550 500 500 600 450 450 650 650 600 600 Xum 550 Xum 500 图7边界条件.(a)导入网格数据：(b)人口：(c)出口 Fig.7 Boundary conditions:(a)imported grid data;(b)entrance;(c)export b 地出现在通道边界的位置上，且较为狭窄；而靠近 通道中心的位置，溶液流动较为复杂，压力无法确 定.随着流体压力的增加，矿粒结构需要承受的压 力也越来越大.基于数值模拟的结果，孔隙通道底 部区域较为狭窄的位置受压力影响较为严重，当 压力超过其承受能力时，矿粒会发生破裂，影响溶 浸效果 图8速度流线图.(a)流线间隔为0.1um:(b)流线间隔为0.01m 3结论 Fig.8 Streamline distribution of the solution:(a)streamline interval of 本文基于CT扫描技术和模拟仿真软件构建 0.1 um;(b)streamline interval of 0.01 um 了次生硫化铜矿粒内部孔隙通道的数值模型，分 渐增大.流速过大会导致溶液与矿石的接触时间 析了孔隙的演变规律，研究了溶浸过程中矿粒内 变短，溶浸过程中化学反应不彻底，无法提取出所 部溶液的渗流作用，研究结果表明： 有的有价组分，造成资源的流失与浪费 (1)利用CT扫描技术、数字图像处理与三维 图10为孔隙通道内溶液所受压力的特征图 重构技术，得到了矿粒孔隙的三维模型，且能准确 像，其中，图10(a)为通道边界上溶液的压力分布 地表征硫化铜矿粒内部的孔隙结构特征 图，从入口到出口，自上而下，压力逐渐增大，最大 (2)浸矿试验中，矿粒内部孔隙度与孔隙连通 压力可达0.8Pa,超过了设置的初始压力(0.715Pa): 度均有大幅度地增长，使用NaCl溶解模拟浸矿过 最大压力出现在孔隙通道底部的位置；图10(b)为 程，可以得到较为理想的孔隙变化，便于孔隙结构 溶液所受压力的等值线图像，速度等值线圈较多 的研究 The arrow points:velocity field Streamline:velocity field Sections:velocity /10-7msJ /10-7ms） Ylum ul(10-7m-s-) ▲7.35×10-7 ▲4.95×10-7 ▲4.75×10- (a) (b)xhum (c) 650 850800750700 650 1 4.5 4.5 X/um 500506000850 750 600 6 700 4.0 GO 4.0 550 35 550 500 500 3.0 50 600 600 550 2.5 600 50 2.0 550 320 2 450 1.0 500 1.0 450 0.5 450 850800 750 700 0 0 Yμm 70 70 图9速度切面图.(a)X方向：(b)Y方向：(c)Z方向 Fig.9 Velocity sections in different directions:(a)Xdirection;(b)Ydirection;(c)Zdirection

渐增大. 流速过大会导致溶液与矿石的接触时间 变短，溶浸过程中化学反应不彻底，无法提取出所 有的有价组分，造成资源的流失与浪费[26] . 图 10 为孔隙通道内溶液所受压力的特征图 像，其中，图 10（a）为通道边界上溶液的压力分布 图，从入口到出口，自上而下，压力逐渐增大，最大 压力可达 0.8 Pa，超过了设置的初始压力（0.715 Pa）； 最大压力出现在孔隙通道底部的位置；图 10（b）为 溶液所受压力的等值线图像，速度等值线圈较多 地出现在通道边界的位置上，且较为狭窄；而靠近 通道中心的位置，溶液流动较为复杂，压力无法确 定. 随着流体压力的增加，矿粒结构需要承受的压 力也越来越大. 基于数值模拟的结果，孔隙通道底 部区域较为狭窄的位置受压力影响较为严重，当 压力超过其承受能力时，矿粒会发生破裂，影响溶 浸效果. 3    结论 本文基于 CT 扫描技术和模拟仿真软件构建 了次生硫化铜矿粒内部孔隙通道的数值模型，分 析了孔隙的演变规律，研究了溶浸过程中矿粒内 部溶液的渗流作用，研究结果表明： （1）利用 CT 扫描技术、数字图像处理与三维 重构技术，得到了矿粒孔隙的三维模型，且能准确 地表征硫化铜矿粒内部的孔隙结构特征. （2）浸矿试验中，矿粒内部孔隙度与孔隙连通 度均有大幅度地增长，使用 NaCl 溶解模拟浸矿过 程，可以得到较为理想的孔隙变化，便于孔隙结构 的研究. (a) 850 800 750 700 600 650 X/μm 600 550 500 500 550 450 Z/μm Y/μm (b) 850 800 750 700 600 650 X/μm 600 550 500 500 550 450 Z/μm Y/μm (c) 750 650 550 500 700 800 850 600 600 550 500 450 Z/μm Y X/μm /μm 图 7    边界条件. （a）导入网格数据；（b）入口；（c）出口 Fig.7    Boundary conditions: (a) imported grid data; (b) entrance; (c) export (a) (b) 图 8    速度流线图. （a）流线间隔为 0.1 μm；（b）流线间隔为 0.01 μm Fig.8    Streamline distribution of the solution: (a) streamline interval of 0.1 μm; (b) streamline interval of 0.01 μm 650 600 600 550 550 450 850 800 750 7 u/(10−7 m·s−1) ▲ 7.35×10−7 ▲ 0 6 5 4 3 2 1 0 700 500 500 (a) X/μm Z/μm Y/μm 650 600 600 550 550 450 850 800 750 u/(10−7 m·s−1) ▲ 4.95×10−7 ▲ 0 700 500 500 (b) X/μm Z/μm Y/μm 4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0 650 600 600 550 550 450 850 800 750 u/(10−7 m·s−1) ▲ 4.75×10−7 ▲ 0 4.5 4.0 3.5 3.0 2.5 2.0 1.5 1.0 0.5 0 700 500 500 (c) X/μm Z/μm Y/μm The arrow points: velocity field Streamline: velocity field Sections: velocity 图 9    速度切面图. （a） X 方向；（b） Y 方向；（c） Z 方向 Fig.9    Velocity sections in different directions: (a) X direction; (b) Y direction; (c) Z direction · 500 · 工程科学学报，第 43 卷，第 4 期

尹升华等：硫化铜矿粒孔隙模型重构与溶液渗流模拟 501· (a) 850 P10-3Pa) (b) 6006s0 850 P/(10-3 Pa) Xμm 800 799.61 800 750 55 750 目691.44 700 500 700 567.33 600 露 525.96 443)1 600 401.84 550 360.47 550 目319.10 277.72 500 236. 500 自94.8 450 841 450 650 29.49 -11.88 600 -53.25 550 9462 500 图10压力分布图.(a)压力分布图：(b)压力等值线图 Fig.10 Pressure profile:(a)pressure distribution;(b)pressure contour map (3)在试验假设与参数设定条件下，溶液渗流 (曾毅君，李建华，李铁球，等.中国含泥轴矿酸法堆浸制粒技术 模拟仿真试验中溶液的流速与压力在孔隙通道底 的应用.轴矿冶，2002,21(4)：182) 部区域和较狭窄的区域急剧增加，对矿粒结构的 [8] Liu XX,Wang L,Xie J P,et al.Experimental study on column 稳定性造成了影响 leach of low grade copper oxide.Min Metall Eng,2016,36(1):83 (刘新星，王龙，谢建平，等.低品位氧化铜矿柱浸试验研究.矿 (4)将矿粉预先制粒再浸矿，矿粒内部孔隙结 冶工程，2016,36(1)：83) 构的变化与溶液的渗流作用相互影响、相互作用， [9] Chen W,Yin S H,Wu A X,et al.Bioleaching of copper sulfides 会影响到浸矿效率与浸出率 using mixed microorganisms and its community structure succession in the presence of seawater.Bioresour Technol,2020. 参考文献 297:122453 [1]Akinci G,Guven D E.Bioleaching of heavy metals contaminated [10]Zhou A,Zhang Q,Bai R N,et al.Characterization of coal micro- sediment by pure and mixed cultures of Acidithiobacillus spp pore structure and simulation on the seepage rules of low-pressure Desalination,2011,268(1-3:221 water based on CT scanning data.Minerals,2016,6(3):78 [2]Zeng J,LiJ,Gou M.et al.Effective strategy for improving sludge [11]Tang B W,Gao S,Wang Y G,et al.Pore structure analysis of treatment rate and microbial mechanisms during chromium electrolytic manganese residue based permeable brick by using bioleaching of tannery sludge.Process Biochem,2019,83:159 industrial CT.Construct Build Mater,2019,208:697 [3]Yin S H,Chen W,Chen X,et al.Bacterial-mediated recovery of [12]Bell S L,Welch G D,Bennett P G.Development of ammoniacal copper from low-grade copper sulfide using acid-processed rice lixiviants for the in-situ leaching of chalcopyrite.Hydromegy. straw.Bioresource Technol,2019,288:121605 1995,39(1-3):11 [4]Yin S H,Wang L M,Wu A X,et al.Progress of research in copper [13]Yin S H,Chen W,Liu J M,et al.Agglomeration experiment of bioleaching technology in China.ChinJEng.019,41(2):143 secondary copper sulfide ore.Chin JEng,019,41(9):1127 (尹升华，王雷鸣，吴爱祥，等.我国铜矿微生物浸出技术的研究 (尹升华，陈威，刘家明，等.次生硫化铜矿制粒试验.工程科学 进展.工程科学学报，2019,41(2)：143) 学报，2019,41(9)：1127) [5]Petersen J.Heap leaching as a key technology for recovery of [14]Nosrati A,Quast K,Xu D F,et al.Agglomeration and column values from low-grade ores-a brief overview.Hydrome. leaching behaviour of nickel laterite ores:effect of ore mineralogy 1 allurgy,2016.165:206 and particle size distribution.Hydromerallurgy,2014,146:29 [6]Chen W.Yin S H.Qi Y,et al.Effect of additives on bioleaching of [15]Yin J S,He R G,Shen K N.Commercial tests of tailing copper sulfide ores.J Central South Univ Sci Technol,2019, briquetting-heap leaching in a certain gold mine.Gold,2007, 50(7):1507 28(2):42 (陈威，尹升华，齐炎，等.添加剂对硫化铜矿生物浸出规律的影 (尹江生，贺锐岗，沈凯宁.某金矿选矿厂尾矿制粒堆浸工业试 响.中南大学学报：自然科学版，2019,50(7)：1507) 验.黄金，2007,28(2)：42) [7]Zeng Y J,Li J H,Li T Q,et al.Technical application of [16]Quaicoe I,Nosrati A,Skinner W,et al.Agglomeration and column agglomerated acidic heap leaching of clay bearing uranium ore in leaching behaviour of goethitic and saprolitic nickel laterite ores. China.Uran Min Metall.2002.21(4):182 Miner Eng,2014,65:1

（3）在试验假设与参数设定条件下，溶液渗流 模拟仿真试验中溶液的流速与压力在孔隙通道底 部区域和较狭窄的区域急剧增加，对矿粒结构的 稳定性造成了影响. （4）将矿粉预先制粒再浸矿，矿粒内部孔隙结 构的变化与溶液的渗流作用相互影响、相互作用， 会影响到浸矿效率与浸出率. 参    考    文    献 Akinci G, Guven D E. Bioleaching of heavy metals contaminated sediment  by  pure  and  mixed  cultures  of  Acidithiobacillus  spp. Desalination, 2011, 268（1-3）: 221 [1] Zeng J, Li J, Gou M, et al. Effective strategy for improving sludge treatment  rate  and  microbial  mechanisms  during  chromium bioleaching of tannery sludge. Process Biochem, 2019, 83: 159 [2] Yin S H, Chen W, Chen X, et al. Bacterial-mediated recovery of copper  from  low-grade  copper  sulfide  using  acid-processed  rice straw. Bioresource Technol, 2019, 288: 121605 [3] Yin S H, Wang L M, Wu A X, et al. Progress of research in copper bioleaching technology in China. Chin J Eng, 2019, 41（2）: 143 （尹升华, 王雷鸣, 吴爱祥, 等. 我国铜矿微生物浸出技术的研究 进展. 工程科学学报, 2019, 41（2）：143） [4] Petersen  J.  Heap  leaching  as  a  key  technology  for  recovery  of values  from  low-grade  ores ——a  brief  overview. Hydrome￾tallurgy, 2016, 165: 206 [5] Chen W, Yin S H, Qi Y, et al. Effect of additives on bioleaching of copper  sulfide  ores. J Central South Univ Sci Technol,  2019, 50（7）: 1507 （陈威, 尹升华, 齐炎, 等. 添加剂对硫化铜矿生物浸出规律的影 响. 中南大学学报:自然科学版, 2019, 50（7）：1507） [6] Zeng  Y  J,  Li  J  H,  Li  T  Q,  et  al.  Technical  application  of agglomerated acidic heap leaching of clay bearing uranium ore in China. Uran Min Metall, 2002, 21（4）: 182 [7] （曾毅君, 李建华, 李铁球, 等. 中国含泥铀矿酸法堆浸制粒技术 的应用. 铀矿冶, 2002, 21（4）：182） Liu X X, Wang L, Xie J P, et al. Experimental study on column leach of low grade copper oxide. Min Metall Eng, 2016, 36（1）: 83 （刘新星, 王龙, 谢建平, 等. 低品位氧化铜矿柱浸试验研究. 矿 冶工程, 2016, 36（1）：83） [8] Chen W, Yin S H, Wu A X, et al. Bioleaching of copper sulfides using  mixed  microorganisms  and  its  community  structure succession in the presence of seawater. Bioresour Technol, 2020, 297: 122453 [9] Zhou A, Zhang Q, Bai R N, et al. Characterization of coal micro￾pore structure and simulation on the seepage rules of low-pressure water based on CT scanning data. Minerals, 2016, 6（3）: 78 [10] Tang  B  W,  Gao  S,  Wang  Y  G,  et  al.  Pore  structure  analysis  of electrolytic  manganese  residue  based  permeable  brick  by  using industrial CT. Construct Build Mater, 2019, 208: 697 [11] Bell S L, Welch G D, Bennett P G. Development of ammoniacal lixiviants for the in-situ leaching of chalcopyrite. Hydrometallurgy, 1995, 39（1-3）: 11 [12] Yin  S  H,  Chen  W,  Liu  J  M,  et  al.  Agglomeration  experiment  of secondary copper sulfide ore. Chin J Eng, 2019, 41（9）: 1127 （尹升华, 陈威, 刘家明, 等. 次生硫化铜矿制粒试验. 工程科学 学报, 2019, 41（9）：1127） [13] Nosrati  A,  Quast  K,  Xu  D  F,  et  al.  Agglomeration  and  column leaching behaviour of nickel laterite ores: effect of ore mineralogy and particle size distribution. Hydrometallurgy, 2014, 146: 29 [14] Yin  J  S,  He  R  G,  Shen  K  N.  Commercial  tests  of  tailing briquetting-heap  leaching  in  a  certain  gold  mine. Gold,  2007, 28（2）: 42 （尹江生, 贺锐岗, 沈凯宁. 某金矿选矿厂尾矿制粒堆浸工业试 验. 黄金, 2007, 28（2）：42） [15] Quaicoe I, Nosrati A, Skinner W, et al. Agglomeration and column leaching behaviour of goethitic and saprolitic nickel laterite ores. Miner Eng, 2014, 65: 1 [16] 650 600 600 550 550 450 850 800 750 799.67 P/(10−3 Pa) 734.85 637.62 507.98 410.75 345.93 281.11 216.29 572.80 443.16 767.26 702.44 605.21 475.57 670.03 540.39 378.34 313.52 248.70 54.25 −6.52 86.66 21.84 183.88 151.47 119.07 700 500 500 (a) X/μm Z/μm Y/μm 650 600 550 850 800 750 700 500 X/μm Y/μm 600 550 450 691.44 P/(10−3 Pa) 567.33 443.21 360.47 277.72 194.98 650.07 484.58 525.96 608.70 401.84 319.10 236.35 −11.88 −94.62 29.49 −53.25 153.61 112.24 70.86 500 (b) Z/μm 图 10    压力分布图. （a）压力分布图；（b）压力等值线图 Fig.10    Pressure profile: (a) pressure distribution; (b) pressure contour map 尹升华等： 硫化铜矿粒孔隙模型重构与溶液渗流模拟 · 501 ·

502 工程科学学报，第43卷，第4期 [17]Yin S H,Wang L M,Chen X,et al.Seepage law of solution inside 律的形响.矿业研究与开发，2011,31(5)：22) ore granular under condition of different heap constructions.J [22]Lo A,Nosrati A,Addai-Mensah J.Particle and pore dynamics Central South Univ Sci Technol,2018,49(4):949 under column leaching of goethitic and saprolitic nickel laterite (尹升华，王雷鸣，陈勋，等.不同堆体结构下矿岩散体内溶液渗 agglomerates.Ady Powder Technol,2016,27(6):2370 流规律.中南大学学报：自然科学版，2018,49(4)：949) [23]Yang B H,Wu A X,Miao XX.3D micropore structure evolution [18]Yin S H,Wu A X,Hu K J,et al.Visualization of flow behavior of ore particles based on image processing.Chin J Eng,2016, during bioleaching of waste rock dumps under saturated and 38(3):328 unsaturated conditions.Hydrometallugy,2013,133:1 (杨保华，吴爱样，缪秀秀.基于图像处理的矿石颗粒三维微观 [19]Xue ZL,Zhang YZ,Liu Z Y,et al.Effects of ore shape on pore 孔隙结构演化.工程科学学报，2016,38(3)：328) structure and seepage velocity field in heap leaching.China Min [24]Bird M B,Butler S L,Hawkes C D,et al.Numerical modeling of Mag,2018,27(12:128 fluid and electrical currents through geometries based on (薛振林，张有志，刘志义，等.矿石形状对浸堆结构及渗流场影 synchrotron X-ray tomographic images of reservoir rocks using 响机制.中国矿业，2018,27(12)：128) Avizo and COMSOL.Comput Geosci,2014,73:6 [20]Liu C.Evolution of the Pore Structure and Flow in Granular Ore [25]Wu A X,Yin S H,Yao B H,et al.Study on preferential flow in Heaps for Copper Oxides Acid Leaching[Dissertation].Beijing: dump leaching of low-grade ores.Hydrometallurgy,2007,87(3- University of Science and Technology Beijing,2017 4):124 (刘超.酸浸条件下氧化铜矿岩散体孔隙结构及渗流演化规律 [26]Lai C S,Yu W Z,Chen S Y.The parameters affecting on heap [学位论文]北京：北京科技大学，2017) leaching of gold ore.Multipurpose Utilization Miner Resour, [21]Wu A X,Wang S Y,Yang B H.Effect of particle structure on 2013(6):9 permeability of leaching dump.Min Res Dev,2011,31(5):22 (赖才书，余文章，陈森煜.彩响金矿堆浸工艺的因素.矿产资源 (吴爱祥，王少勇，杨保华.堆浸散体颗粒结构对溶浸液渗流规 综合利用，2013(6)：9)

Yin S H, Wang L M, Chen X, et al. Seepage law of solution inside ore  granular  under  condition  of  different  heap  constructions. J Central South Univ Sci Technol, 2018, 49（4）: 949 （尹升华, 王雷鸣, 陈勋, 等. 不同堆体结构下矿岩散体内溶液渗 流规律. 中南大学学报: 自然科学版, 2018, 49（4）：949） [17] Yin S H, Wu A X, Hu K J, et al. Visualization of flow behavior during  bioleaching  of  waste  rock  dumps  under  saturated  and unsaturated conditions. Hydrometallurgy, 2013, 133: 1 [18] Xue Z L, Zhang Y Z, Liu Z Y, et al. Effects of ore shape on pore structure  and  seepage  velocity  field  in  heap  leaching. China Min Mag, 2018, 27（12）: 128 （薛振林, 张有志, 刘志义, 等. 矿石形状对浸堆结构及渗流场影 响机制. 中国矿业, 2018, 27（12）：128） [19] Liu C. Evolution of the Pore Structure and Flow in Granular Ore Heaps for Copper Oxides Acid Leaching[Dissertation].  Beijing: University of Science and Technology Beijing, 2017 （ 刘超. 酸浸条件下氧化铜矿岩散体孔隙结构及渗流演化规律 [学位论文]. 北京: 北京科技大学, 2017） [20] Wu  A  X,  Wang  S  Y,  Yang  B  H.  Effect  of  particle  structure  on permeability of leaching dump. Min Res Dev, 2011, 31（5）: 22 （吴爱祥, 王少勇, 杨保华. 堆浸散体颗粒结构对溶浸液渗流规 [21] 律的影响. 矿业研究与开发, 2011, 31（5）：22） Lo  A,  Nosrati  A,  Addai-Mensah  J.  Particle  and  pore  dynamics under  column  leaching  of  goethitic  and  saprolitic  nickel  laterite agglomerates. Adv Powder Technol, 2016, 27（6）: 2370 [22] Yang B H, Wu A X, Miao X X. 3D micropore structure evolution of  ore  particles  based  on  image  processing. Chin J Eng,  2016, 38（3）: 328 （杨保华, 吴爱祥, 缪秀秀. 基于图像处理的矿石颗粒三维微观 孔隙结构演化. 工程科学学报, 2016, 38（3）：328） [23] Bird M B, Butler S L, Hawkes C D, et al. Numerical modeling of fluid  and  electrical  currents  through  geometries  based  on synchrotron  X-ray  tomographic  images  of  reservoir  rocks  using Avizo and COMSOL. Comput Geosci, 2014, 73: 6 [24] Wu A X, Yin S H, Yao B H, et al. Study on preferential flow in dump  leaching  of  low-grade  ores. Hydrometallurgy,  2007,  87（3- 4）: 124 [25] Lai  C  S,  Yu  W  Z,  Chen  S  Y.  The  parameters  affecting  on  heap leaching  of  gold  ore. Multipurpose Utilization Miner Resour, 2013（6）: 9 （赖才书, 余文章, 陈森煜. 影响金矿堆浸工艺的因素. 矿产资源 综合利用, 2013（6）：9） [26] · 502 · 工程科学学报，第 43 卷，第 4 期