《工程科学学报》录用稿，htps:/doi.org/10.13374/i,issn2095-9389.2021.11.02.006©北京科技大学2020 工程科学学报，第x卷，第x期，xx-Xx,2021年10月 Chinese Journal of Engineering,Vol.xx,No.xx-xx,October 2021 DOI:XXXXXXXXXXXXXXX 非饱和矿堆溶液渗流迟滞与毛细扩散行为表征* 王雷鸣2,4)，李希雯2区，尹升华2，周根茂》，李辉2)，刘培延，邓博纳) 1)金属矿山高效开采与安全教有部重点实验室，北京，1000832)北京科技大学土木与资源工程学院，北茶 3)核工业北京化工治金研究院，北京，1011494)金属矿山安全与健康国家重点实验室，安徽，马鞍山， 4 5)武汉工程大学绿色化工过程教有部重点实验室，湖北，武汉，430205 ☒通信作者，E-mail:heavenli1@l63.com 摘要为深入理解非饱和矿堆内溶浸液毛细渗流扩散以及渗流迟滞行为 本文构建适于非饱和矿堆的毛细渗流模型，利用 COMSOL Multiphysics有限元数值平台开展毛细渗流可视化模拟研究 并利用时域反射器(Time Domain Reflector,TDR) 实时探测了非饱和堆内持液率变化，探索了基于Design Expert的免细渗流过程多因素响应规律，讨论了非饱和矿堆持液率、 毛细吸力、孔隙率与喷淋强度间的潜在关联机制。研究结果表明：孔隙率对矿堆持液率的影响高于喷淋强度，矿堆持液率随 喷淋时间的增长收敛性增加，且孔隙率小的矿堆需要更长的时间才能达到稳态持液：不考虑溶液喷淋强度影响时，矿堆持液 率与孔隙比、水力传导系数呈正相关：特别是在喷淋初期020s),喷淋强度、水力传导系数和孔隙比对矿堆持液率的影 响更为显著：初步构建了考虑气液两相运移的非饱和矿堆溶液毛细渗流模型：毛细吸力的变化对孔隙率较小的矿堆更敏感： 喷淋强度较大、孔隙比越小时，矿堆底部的毛细吸力趣大，更易达到稳态持液状态。 关键词非饱和矿堆：流体流动行为：渗流迟潢 :毛细吸力：持液率；COMSOL Multiphysics 分类号TD853 Characterization f liquid seepage hysteresis and capillary diffusion behavior in unsaturated ore heap WANG Lei-ming 2,Li Xi-wen YIN Sheng-hua2),ZHOU Gen-mao,LI Hu),LIU Pei-zheng )DENG Bo-na 1) Key Laboratory of Ministry of Education for High-Efficient Mining and Safety of Metal,Beijing,100083,China; 2) School of Civil and Resource Engineering,University of Science and Technology Beijing,Beijing,100083,China 3) Key Laboratory of Green Chemical Engineering Process of Ministry of Education,Wuhan Institute of Technology,Hubei Province,430205,China 为 State Key Laboratory of Safety and Health for Metal Mines,Maanshan,Anhui Province,243000,China “收稿日期：2021-10-26 基金项目：金属矿山安全与健康国家重点实验室开放基金项目(2021-JSKSSYS-01):绿色化工过程教育部重点实验室开 放基金(GCP202108):国家自然科学基金重点项目(52034001,51734001)：山东省重大科技创新工程项目(2019SDZ Y05):国家科技部重点领域创新团队资助项目(2018RA400):高等学校学科创新引智计划资助(B20041):中央高校 基本科研业务费专项资金资助项目(FRF-TP-18-O03C1)
工程科学学报，第 x 卷，第 x 期，xx-xx，2021 年 10 月 Chinese Journal of Engineering, Vol.xx, No. xx-xx, October 2021 DOI：xxxxxxxxxxxxxxx 非饱和矿堆溶液渗流迟滞与毛细扩散行为表征* 王雷鸣 1,2,4,5)，李希雯 2)，尹升华 1,2)，周根茂 3)，李 辉 2)，刘培正 4)，邓博纳 5) 1) 金属矿山高效开采与安全教育部重点实验室，北京，100083 2) 北京科技大学土木与资源工程学院，北京，100083 3) 核工业北京化工冶金研究院，北京，101149 4) 金属矿山安全与健康国家重点实验室，安徽，马鞍山，243000 5) 武汉工程大学绿色化工过程教育部重点实验室，湖北，武汉，430205 通信作者，E-mail: heavenli1@163.com 摘 要 为深入理解非饱和矿堆内溶浸液毛细渗流扩散以及渗流迟滞行为，本文 构建适于非饱和矿堆的毛细渗流模型，利用 COMSOL Multiphysics 有限元数值平台开展毛细渗流可视化模拟研究，并利用时域反射器（Time Domain Reflector，TDR） 实时探测了非饱和堆内持液率变化，探索了基于 Design Expert 的毛细渗流过程多因素响应规律，讨论了非饱和矿堆持液率、 毛细吸力、孔隙率与喷淋强度间的潜在关联机制。研究结果表明：孔隙率对矿堆持液率的影响高于喷淋强度，矿堆持液率随 喷淋时间的增长收敛性增加，且孔隙率小的矿堆需要更长的时间才能达到稳态持液；不考虑溶液喷淋强度影响时，矿堆持液 率与孔隙比、水力传导系数呈正相关；特别是在喷淋初期（0~20s），喷淋强度、水力传导系数和孔隙比对矿堆持液率的影 响更为显著；初步构建了考虑气液两相运移的非饱和矿堆溶液毛细渗流模型；毛细吸力的变化对孔隙率较小的矿堆更敏感； 喷淋强度较大、孔隙比越小时，矿堆底部的毛细吸力越大，更易达到稳态持液状态。 关键词 非饱和矿堆；流体流动行为；渗流迟滞；毛细吸力；持液率；COMSOL Multiphysics 分类号 TD853 Characterization of liquid seepage hysteresis and capillary diffusion behavior in unsaturated ore heap WANG Lei-ming1,2) , LI Xi-wen2) , YIN Sheng-hua1,2) , ZHOU Gen-mao3) , LI Hui2) , LIU Pei-zheng4) , DENG Bo-na5) 1) Key Laboratory of Ministry of Education for High-Efficient Mining and Safety of Metal, Beijing, 100083, China; 2) School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing, 100083, China 3) Key Laboratory of Green Chemical Engineering Process of Ministry of Education, Wuhan Institute of Technology, Hubei Province, 430205, China 4) State Key Laboratory of Safety and Health for Metal Mines, Maanshan, Anhui Province, 243000, China ** 收稿日期：2021-10-26 基金项目：金属矿山安全与健康国家重点实验室开放基金项目（2021-JSKSSYS-01）；绿色化工过程教育部重点实验室开 放基金（GCP202108）；国家自然科学基金重点项目（52034001，51734001）；山东省重大科技创新工程项目（2019SDZ Y05）；国家科技部重点领域创新团队资助项目（2018RA400）；高等学校学科创新引智计划资助（B20041）；中央高校 基本科研业务费专项资金资助项目（FRF-TP-18-003C1） 《工程科学学报》录用稿，https://doi.org/10.13374/j.issn2095-9389.2021.11.02.006 ©北京科技大学 2020 录用稿件，非最终出版稿
5) Beijing Research Institute of Chemical Engineering Metallurgy,Beijing,101149,China Corresponding author,E-mail:heavenlil@163.com ABSTRACT To deeply understand the capillary diffusion and seepage hysteresis behavior of the leaching solution in the unsaturated ore heaps,this paper builds a capillary seepage model suitable for the unsaturated ore heap,uses COMSOL Multiphysics finite element numerical platform to carry out the capillary seepage visual simulation study,and uses the Time Domain Reflector (TDR)to in-situ, real-time detect the liquid holdup changes of unsaturated heap,explores the multi-factor response regulations of capillary seepage process based on Design Expert,and discusses the potential connection mechanism among the liquid holdup,capillary suction, porosity and irrigation rate of unsaturated ore heaps.The research results show that the heap porosity has an obvious impact on the heap liquid holdup than the irrigation intensity,and the increases convergence of liquid holdup improve with the spraying time,and ore heap with a small porosity takes longer time to reach steady status of liquid holdup.When the effect of liquid irrigation is not considered,the heap iquid holdup is positively correlated with the porosity ratio and hydraulic coduEspecially in the initial stage of irrigation period(0~20 s),the effects of irrigation rate,hydraulic conductivity and porosity ratio onore heap liquid holdup is more significant,an unsaturated ore pile solution capillary seepage model considering gas-liquid two-phase migration is preliminarily constructed;The capillary suction is more sensitive in the ore heap with a smaller porosity;The larger the irrigation rate and the smaller the porosity,the greater the capillary suction at the bottom of the ore heap,and it is easie ore heap to reach a steady-state of liquid holdup. KEY WORDS unsaturated ore heap:fluid flow behavior,seepage hysteresis apillary suction;liquid holdup;COMSOL multiphysics 堆浸技术凭借其基建周期短、投入成本低、作业安全性高等优势，被广泛应用于低品位硫化铜矿、砂 岩型金矿等战略矿产资源回采作业，。其中，堆浸是子种汽、固、液多相介质并存、渗流场、温度场、化 学场等多场耦合的复杂反应体系)；该体系中，溶液是溶准氧、Cu/Cu、Fe+Fe3#等有价金属离子等溶 质赋存的重要媒介，直接影响着反应传质过程与浸矿效率。因此，能否有效表征堆内溶液毛细扩散、流 动与分布规律是制约矿物浸出效率的重要因素。作为之种非饱和多孔介质中多相流，溶浸液的滞后现象是 普遍存在且十分复杂的阿。己有研究发现堆内溶液流动主要有两种形式6，刀：重力驱动的优势流和微观力驱 动的毛细扩散流。其中，由毛细管力驱动的毛细扩散过程，对溶液渗流迟滞行为的影响是十分显著的网。 毛细扩散过程不仅是影响矿堆渗流迟滞行为，更重要的会影响传质过程与浸矿效率，需要进一步深刻理解 和认识。 目前，针对非饱和堆内溶液渗流迟滞行为规律，特别是毛细扩散过程，国内外专家学者主要从物理实 验、 无损探测与数值模拟三介层面展开。其中，英国帝国理工大学Ilankoon等9，o开展矿堆溶液渗流迟滞 实验，利用持液率、残余稳念特液率对非饱和矿石颗粒堆的持液行为进行了研究：王雷鸣等立足于制粒 矿堆，基于自主研发的待液猴为原位监测系统实现了对不可动液、可动液的量化表征：毛细扩散过程对非 饱和堆（如破碎矿堆、制粒矿堆等）持液行为的影响不可忽略，毛细入渗率与毛细速率、含水率呈正相 关3，：随着多种无扰动探测设备技术的进步与引入，对溶液渗流的研究逐渐步入了微细观层面(，英 国剑桥大学Faga等利用核磁共振技术(Magnetic Resonance Imaging,.MRI)对堆内溶液扩散和溶液优势 流形成过程进行了研究，薛振林等s利用粒子图像测速技术(Particle Image Velocity,.PIV)对矿堆结构内溶 液渗流场进行了无损探测，有效实现了饱和溶液渗流场的可视化等。然而，单纯的实验研究与无扰动探测 存在实验误差大、瞬态短时、随机性强等劣势，难以实现对溶液渗流过程进行有效的预测。对此，多类数 学表征模型与渗流模拟软件逐渐被提出应用：国外Dixon、Petersen和Bouffard等l920合作研发了堆浸模拟 程序HeapSIM,可以对工业堆浸喷淋、矿物浸出实现过程监测与结果预测：我国W等2，利用COMSOL Multiphysics数值模拟软件和计算机断层扫描技术(Computed Tomography,CT)构建非饱和矿堆网格模型 并对堆内溶液渗流过程模拟，证实了堆内优势流动区与溶液停滞区：此外，马丹等2训聚焦裂隙岩体渗流， 利用Fluent开展了数值模拟研究，实现了破裂砂岩内的溶液流动规律。 综合已有研究不难发现，以往研究通常利用短时、瞬态可视化表征手段、数值模拟软件与相似物理渗
5) Beijing Research Institute of Chemical Engineering Metallurgy, Beijing, 101149, China Corresponding author, E-mail: heavenli1@163.com ABSTRACT To deeply understand the capillary diffusion and seepage hysteresis behavior of the leaching solution in the unsaturated ore heaps, this paper builds a capillary seepage model suitable for the unsaturated ore heap, uses COMSOL Multiphysics finite element numerical platform to carry out the capillary seepage visual simulation study, and uses the Time Domain Reflector (TDR) to in-situ, real-time detect the liquid holdup changes of unsaturated heap, explores the multi-factor response regulations of capillary seepage process based on Design Expert, and discusses the potential connection mechanism among the liquid holdup, capillary suction, porosity and irrigation rate of unsaturated ore heaps. The research results show that the heap porosity has an obvious impact on the heap liquid holdup than the irrigation intensity, and the increases convergence of liquid holdup improve with the spraying time, and ore heap with a small porosity takes longer time to reach steady status of liquid holdup. When the effect of liquid irrigation is not considered, the heap liquid holdup is positively correlated with the porosity ratio and hydraulic conductivity. Especially in the initial stage of irrigation period (0~20 s), the effects of irrigation rate, hydraulic conductivity and porosity ratio on ore heap liquid holdup is more significant; an unsaturated ore pile solution capillary seepage model considering gas-liquid two-phase migration is preliminarily constructed; The capillary suction is more sensitive in the ore heap with a smaller porosity; The larger the irrigation rate and the smaller the porosity, the greater the capillary suction at the bottom of the ore heap, and it is easier for ore heap to reach a steady-state of liquid holdup. KEY WORDS unsaturated ore heap; fluid flow behavior; seepage hysteresis; capillary suction; liquid holdup; COMSOL multiphysics 堆浸技术凭借其基建周期短、投入成本低、作业安全性高等优势，被广泛应用于低品位硫化铜矿、砂 岩型金矿等战略矿产资源回采作业[1, 2]。其中，堆浸是一种气、固、液多相介质并存、渗流场、温度场、化 学场等多场耦合的复杂反应体系[3]；该体系中，溶液是可溶性氧、Cu2+/Cu+、Fe2+/Fe3+等有价金属离子等溶 质赋存的重要媒介，直接影响着反应传质过程与浸矿效率[4]。因此，能否有效表征堆内溶液毛细扩散、流 动与分布规律是制约矿物浸出效率的重要因素。作为一种非饱和多孔介质中多相流，溶浸液的滞后现象是 普遍存在且十分复杂的[5]。已有研究发现堆内溶液流动主要有两种形式[6, 7]：重力驱动的优势流和微观力驱 动的毛细扩散流。其中，由毛细管力驱动的毛细扩散过程，对溶液渗流迟滞行为的影响是十分显著的[8]。 毛细扩散过程不仅是影响矿堆渗流迟滞行为，更重要的会影响传质过程与浸矿效率，需要进一步深刻理解 和认识。 目前，针对非饱和堆内溶液渗流迟滞行为规律，特别是毛细扩散过程，国内外专家学者主要从物理实 验、无损探测与数值模拟三个层面展开。其中，英国帝国理工大学 Ilankoon 等[9, 10]开展矿堆溶液渗流迟滞 实验，利用持液率、残余稳态持液率对非饱和矿石颗粒堆的持液行为进行了研究；王雷鸣等[11]立足于制粒 矿堆，基于自主研发的持液行为原位监测系统实现了对不可动液、可动液的量化表征；毛细扩散过程对非 饱和堆（如破碎矿堆、制粒矿堆等）持液行为的影响不可忽略[12]，毛细入渗率与毛细速率、含水率呈正相 关[13, 14]；随着多种无扰动探测设备技术的进步与引入，对溶液渗流的研究逐渐步入了微细观层面[15, 16]，英 国剑桥大学 Fagan 等[17]利用核磁共振技术（Magnetic Resonance Imaging, MRI）对堆内溶液扩散和溶液优势 流形成过程进行了研究，薛振林等[18]利用粒子图像测速技术（Particle Image Velocity, PIV）对矿堆结构内溶 液渗流场进行了无损探测，有效实现了饱和溶液渗流场的可视化等。然而，单纯的实验研究与无扰动探测 存在实验误差大、瞬态短时、随机性强等劣势，难以实现对溶液渗流过程进行有效的预测。对此，多类数 学表征模型与渗流模拟软件逐渐被提出应用：国外 Dixon、Petersen 和 Bouffard 等[19, 20]合作研发了堆浸模拟 程序 HeapSIM，可以对工业堆浸喷淋、矿物浸出实现过程监测与结果预测；我国 Wu 等[21, 22]利用 COMSOL Multiphysics 数值模拟软件和计算机断层扫描技术（Computed Tomography, CT）构建非饱和矿堆网格模型 并对堆内溶液渗流过程模拟，证实了堆内优势流动区与溶液停滞区；此外，马丹等[23, 24]聚焦裂隙岩体渗流， 利用 Fluent 开展了数值模拟研究，实现了破裂砂岩内的溶液流动规律。 综合已有研究不难发现，以往研究通常利用短时、瞬态可视化表征手段、数值模拟软件与相似物理渗 录用稿件，非最终出版稿
流实验等方式开展研究，多聚焦溶液优势流动、重力流动、饱和条件、横向毛细扩散等方面s2”,对非饱 和矿堆体系的上向毛细渗流扩散、渗流迟滞行为的认识仍不深刻和不清晰，严重制约着浸矿过程中溶液渗 流调控方法与措施，相关研究有待进一步深入开展。对此，加拿大Top即等提出了时域反射方法，该法利 用电磁波在不同介质中的传播速度的差异来测定土壤含水率的一种方法，具有良好的可预测性2，，因此 可以被有效借鉴到本研究当中。 对此，本文聚焦非饱和矿石颗粒堆，开展毛细溶液渗流迟滞与毛细扩散行为表征模拟，构建非饱和矿 堆渗流模型，利用时域反射器(Time Domain Reflector)实时探测并计算非饱和矿堆内持液率，并结合 Design Expert进行响应曲面分析，深入探究非饱和矿堆内溶液渗流过程中各影响因素之间的相互关系。运 用COMSOL Multiphysics开展溶液毛细渗透过程模拟，相关研究结果对调控矿堆内溶液渗流、强化溶液毛 细扩散，进而提高矿物浸出效率具有良好的指导意义。 石程维见食美德修侧宝装木假货使心筒化为矩形 1非饱和堆涪液渗流模型构建 1.1基本假设 化为矩形二维平面进行 分析，建立尺寸为0.05m×0.2m的数学模型。并且，做出如下基本假设 (1)矿石颗粒的骨架形状、颗粒相对位置不随时间发生变化： (2)颗粒间孔隙率大小、形状、连通性不随时间发生变化 (3)渗流过程中流体不可压缩： (4)矿石颗粒堆内孔隙率、渗透率各向同性： (7)模拟环境温度不变。 1.2控制方程 为探明矿堆非饱和浸出渗流规律，本文以界面作用为切入点，分析浸出液渗流毛细扩散状态。在矿 堆渗透过程中，基于达西渗流定律、质量守恒定律和两相流理论叫，引入毛细吸力和体积含水率的变量， 录用 分别推导矿堆非饱和渗流过程中液相利气相的运移方程，如式1、式2所示： -p) +V Kk.-(Vp.+P.gVD) =0 (1) k.z(Vp.+P.gVD) =( (2) C式中：k为饱和渗透率，m2:k为相对渗透率：n为运动粘度：g为重力加速度，kgms2:D为垂 直坐标：Q为持液率，其中，Q的变化范围从残余持液率0增大到饱和持液率0。 目前，矿堆水分特征曲线目前具有代表性的模型较多，运用较为广泛的模型是非饱和渗流van Genuchten-Mualem(VGM)模型B,为求解函数Pc,根据持液率和渗透率的关系，分别定义O、C、Se和 灯，如式(3)~(6)所示： [0n+Se(0n-0.w) H。>0 w= 日.w H≤0 (3)
流实验等方式开展研究，多聚焦溶液优势流动、重力流动、饱和条件、横向毛细扩散等方面[25-27]，对非饱 和矿堆体系的上向毛细渗流扩散、渗流迟滞行为的认识仍不深刻和不清晰，严重制约着浸矿过程中溶液渗 流调控方法与措施，相关研究有待进一步深入开展。对此，加拿大 Topp 等提出了时域反射方法，该法利 用电磁波在不同介质中的传播速度的差异来测定土壤含水率的一种方法，具有良好的可预测性[28, 29]，因此 可以被有效借鉴到本研究当中。 对此，本文聚焦非饱和矿石颗粒堆，开展毛细溶液渗流迟滞与毛细扩散行为表征模拟，构建非饱和矿 堆渗流模型，利用时域反射器（Time Domain Reflector）实时探测并计算非饱和矿堆内持液率，并结合 Design Expert 进行响应曲面分析，深入探究非饱和矿堆内溶液渗流过程中各影响因素之间的相互关系。运 用 COMSOL Multiphysics 开展溶液毛细渗透过程模拟，相关研究结果对调控矿堆内溶液渗流、强化溶液毛 细扩散，进而提高矿物浸出效率具有良好的指导意义。 1 非饱和矿堆溶液渗流模型构建 1.1 基本假设 为研究溶液在浸出矿堆中的渗流规律，本文构建矿石颗粒堆简化模型，将其简化为矩形二维平面进行 分析，建立尺寸为 0.05m×0.2m 的数学模型。并且，做出如下基本假设： （1） 矿石颗粒的骨架形状、颗粒相对位置不随时间发生变化； （2） 颗粒间孔隙率大小、形状、连通性不随时间发生变化； （3） 渗流过程中流体不可压缩； （4） 矿石颗粒堆内孔隙率、渗透率各向同性； （5） 矿石颗粒堆的孔隙空间中存在气液两相流动状态； （6） 气液两相流动均为低速低扰动的达西渗流； （7） 模拟环境温度不变。 1.2 控制方程 为探明矿堆非饱和浸出渗流规律，本文以界面作用为切入点，分析浸出液渗流毛细扩散状态[30]。在矿 堆渗透过程中，基于达西渗流定律、质量守恒定律和两相流理论[31]，引入毛细吸力和体积含水率的变量， 分别推导矿堆非饱和渗流过程中液相和气相的运移方程，如式 1、式 2 所示： （1） （2） Cw式中：κ 为饱和渗透率，m 2；krw为相对渗透率；η 为运动粘度；g 为重力加速度，kg/m·s2；D 为垂 直坐标；θw为持液率，其中，θw的变化范围从残余持液率 θr增大到饱和持液率 θs。 目前，矿堆水分特征曲线目前具有代表性的模型较多，运用较为广泛的模型是非饱和渗流 van Genuchten-Mualem（VGM）模型[32]，为求解函数 Pc，根据持液率和渗透率的关系，分别定义 θ、C、Se 和 kr，如式（3）~（6）所示： （3） 录用稿件，非最终出版稿
k=[1-1-门 (4) am(6.-0）se1-5e）” C=1 H.>0 (5) 0 H.≤0 He Pe/pwg (6) 迭代后，可求得气相的VGM公式： 0a=03-0,w-0ra (7) ks1-(1-) (8) Sea-1-Se (9) C。=-C (10) 1.3物理模型及网格脚分 终出版稿 为探明矿堆非饱和浸出渗流过程中毛细上升规律， 定义横型为经过对称轴的任意截面。模型的宽为5 cm,高为20cm:设置矿堆底存在高3cm、水平方向无天的恒定饱和流动区。 利用COMSOL Multiphysies5.0数值模拟软件，y构建经维网格模型，如图1所示。由图1可见，较矿 堆中部而言，矿堆顶部和底部的边界网格结构更为细密， 其可更好的表征溶液毛细扩散规律。 (a) 绿用稿件 00 3cm Saturated liquid region 0.20 0.20 (c) 0.15 0.15 0.10 0.10 0.05 0.05 0 -0.05 00.05 0.10 -0.10-0.05 00.05 0.100.15
（4） （5） （6） 迭代后，可求得气相的 VGM 公式： （7） （8） （9） （10） 1.3 物理模型及网格划分 为探明矿堆非饱和浸出渗流过程中毛细上升规律，定义模型为经过对称轴的任意截面。模型的宽为 5 cm，高为 20 cm；设置矿堆底存在高 3cm、水平方向无穷大的恒定饱和流动区。 利用 COMSOL Multiphysics 5.0 数值模拟软件，构建二维网格模型，如图 1 所示。由图 1 可见，较矿 堆中部而言，矿堆顶部和底部的边界网格结构更为细密，其可更好的表征溶液毛细扩散规律。 录用稿件，非最终出版稿
圖1柱浸颗粒堆物理模型构建及网格划分 Fig.1 Physical model and mesh of packed heap in column leaching 借鉴工业矿堆渗流的实际情况及相关文献数据，选取模型主要参数，如表1所示。 表1败学模型的关楚散 Table 1 Key parameters of mathematical model Parameter Symbol Value Unit Liquid density A 1×103 kg.m3 Gas density A 1.28 kg.m Liquid viscosity 103 Pa.S Gas viscosity 加 1.81×105 Pa.S Gravitational acceleration 9.82 m's Residual liquid holdup 0.01 1Size parameters of VGM model 1.89 2n4 Size parameters of VGM model 水 2.811 3 Size parameters of VGM model M N 1.4边界条件 建模前对模型边界条件进行设定。初始状态，矿堆中的液是符合静压分布的。而液相从底部向上运 动，气相是在矿堆顶部出现，且不会消失。假定气液两相流动均为低速低扰动的达西渗流，在入口处、表 面、边壁，气液两相满足式11。 n[-5(Vp (11) 在底部，液相和气相的压力为恒压，如式12所示： p=Po(t) (12) 1.5数值模拟方案与赋值 为进一步探究溶液在浸矿纳的溶液渗流迟滞规律，将物理实验结果和数值模拟相结合，利用孔隙 率对渗透系数进行赋值考察不同孔隙率、喷淋强度对非饱和矿堆渗流过程：探索渗流速度、毛细吸力、 持液率、非饱和系数等相关特性。具体模拟和实验方案，如表2和表3所示。 依据初始实验条件差异，各组互为对照，对实验方案进行简要阐释，如下：(1)不同喷淋强度与孔 隙率条件下，共没置6个实验组，喷淋强度分3个阶梯，即：0、10、50Lmh,矿石颗粒堆孔隙结构用 孔隙比进行表征，选择2种，即：1.04、1.17。(2)不同孔隙率条件下，不进行喷淋作业，即仅依赖毛细 吸力和下部饱和溶液区供给颗粒堆，孔隙比介于1.04117，孔隙率介于51~54%，相应地，矿石颗粒堆的 水力传导系数介于0.02~0.08cms'。 表2不同喷洲强度和孔麻率缘件下存液毛细渗流模拟方案 Table 2 Experimental scheme of liquid capillarity seepage Experimental group Irrigation rate/(L.mh) Porosity ratio Hydraulic conductivity/(cm's) 0 1.040816 0.02 A2 1.173913 0.08
图 1 柱浸颗粒堆物理模型构建及网格划分 Fig.1 Physical model and mesh of packed heap in column leaching 借鉴工业矿堆渗流的实际情况及相关文献数据，选取模型主要参数，如表 1 所示。 表 1 数学模型的关键参数 Table 1 Key parameters of mathematical model Parameter Symbol Value Unit Liquid density ρw 1×103 kg·m-3 Gas density ρa 1.28 kg·m-3 Liquid viscosity ηw 10-3 Pa·S Gas viscosity ηa 1.81×10-5 Pa·S Gravitational acceleration g 9.82 m·s-2 Residual liquid holdup θr 0.01 % 1 st Size parameters of VGM model α 1.89 m-1 2 nd Size parameters of VGM model N 2.811 3 rd Size parameters of VGM model M 1-1/N 1.4 边界条件 建模前对模型边界条件进行设定。初始状态，矿堆中的气液是符合静压分布的。而液相从底部向上运 动，气相是在矿堆顶部出现，且不会消失。假定气液两相流动均为低速低扰动的达西渗流，在入口处、表 面、边壁，气液两相满足式 11。 （11） 在底部，液相和气相的压力为恒压，如式 12 所示： （12） 1.5 数值模拟方案与赋值 为进一步探究溶液在浸出矿堆内的溶液渗流迟滞规律，将物理实验结果和数值模拟相结合，利用孔隙 率对渗透系数进行赋值，考察不同孔隙率、喷淋强度对非饱和矿堆渗流过程；探索渗流速度、毛细吸力、 持液率、非饱和系数等相关特性。具体模拟和实验方案，如表 2 和表 3 所示。 依据初始实验条件差异，各组互为对照，对实验方案进行简要阐释，如下：（1）不同喷淋强度与孔 隙率条件下，共设置 6 个实验组，喷淋强度分 3 个阶梯，即：0、10、50 L·m-2h -1，矿石颗粒堆孔隙结构用 孔隙比进行表征，选择 2 种，即：1.04、1.17。（2）不同孔隙率条件下，不进行喷淋作业，即仅依赖毛细 吸力和下部饱和溶液区供给颗粒堆，孔隙比介于 1.04~1.17，孔隙率介于 51~54%，相应地，矿石颗粒堆的 水力传导系数介于 0.02~0.08 cm·s-1。 表 2 不同喷淋强度和孔隙率条件下溶液毛细渗流模拟方案 Table 2 Experimental scheme of liquid capillarity seepage Experimental group Irrigation rate/ (L·m-2h -1) Porosity ratio Hydraulic conductivity/ (cm·s-1) A1 0 1.040816 0.02 A2 0 1.173913 0.08 录用稿件，非最终出版稿