地浸采铀的多重耦合动力学模型与三维模拟程序.pdfVIP

中国核科学技术进展报告（第一卷） 铀矿冶分卷 Progress Report on China Nuclear Science Technology （Vol.1 ） 2009 年11 月 地浸采铀的多重耦合动力学模型及三维模拟程序1 谭凯旋，曾 晟，桑 潇，孙 冰 (南华大学 核资源与核燃料工程学院，湖南 衡阳 421001) 摘 要：地浸采铀是一个复杂的非线性动力系统，涉及到多种因素与过程的耦合和反馈作用。为此建立了一个多因素、 多过程耦合的综合动力学模型和模拟方法来研究地浸采铀。该模型主要包括：(1) 岩石结构力学及其演化；(2)岩石变 形的递增应力流变学；(3)三维粘弹性/塑性多重岩石变形作用；(4) 水力断裂作用；(5) 张量断裂和岩石渗透率；(6) 流 体-岩石反应和质量迁移；(7)化学溶解引起的压实作用；(8)多相流体流动。基于Fortran 与C++语言编制了3 维模拟 程序。以新疆某地浸矿山为例，对地浸采铀的酸化、生产、溶浸末期等各阶段进行了数值模拟与分析。该模型和程 序可应用于地浸采铀的理论研究、矿山设计、生产管理和地下水污染物迁移与复原研究。 关键词：多重耦合；动力学模型；数值模拟；地浸采铀 地浸是将溶浸液直接注入地下含矿岩层中，利用矿物与水溶液的化学反应来获取有用金属和化 合物的一种新的采矿理论与采矿方法，目前已广泛应用于砂岩型铀矿的开采[1-3] 。地浸采铀是一个极 其复杂的非线性动力系统, 一方面含矿含水层都是具有初始非均匀性的，不仅孔隙度大小在空间上 是变化的，而且矿物组成特别是有价矿物含量（矿石品位）在空间上也是非均匀的和变化较大的； 另一方面，地浸中涉及到许多因素与过程的耦合和反馈作用，主要的有溶浸液的流动、水溶液组分 的扩散、矿物溶解与沉淀反应、溶解压实与水力断裂等引起的岩石变形作用、孔隙度与渗透率的动 力演化等，有的耦合具正反馈作用，而有的耦合可能是负反馈[4-5] 。因此，建立一个多因素、多过程 耦合的综合动力学模型和模拟方法来研究地浸采铀是有必要的。 1 动力学模型 1.1 原地溶浸采矿过程中的耦合与反馈作用 原地溶浸采矿涉及到许多过程的耦合作用，其中最主要的有溶浸液的流动、水溶液组分的扩散、 矿物溶解、成核、矿物颗粒的生长等。当可反应流体在孔隙岩矿中流动时，矿物的溶解将使孔隙度 增大、渗透率增高，使得流体流动会在渗透率更高的地方汇聚。所有天然岩矿石都是具有初始非均 匀性的，如孔隙度大小、矿石品位在空间上是非均匀的和变化较大的。溶浸采矿过程中，浸矿液在 这种初始非均匀岩矿石中流动会在孔隙度大的局部地方流体量大，流体量大有利于矿物的溶解（浸 出），矿物的溶解使岩矿石孔隙度增大，孔隙度的增大促使流入的流体量增大，流量的增大又促进矿 物的溶解并增大孔隙度，这样就会由于流体流动与化学反应的耦合产生正的反馈作用。 溶浸采矿过程中，除了流体流动与溶解反应发生耦合与正反馈作用外，流体中溶解组分的扩散 还将与溶解反应产生耦合与负反馈作用，此时从指形前锋的两侧壁处溶解的矿物溶质通过扩散迁移 到它的顶部，推动顶部处的流体趋向饱和而使溶解反应停止。这种溶解组分扩散与反应的负反馈作 用将抑制反应前锋的非均衡性和不稳定性，阻止指形前锋的无限制生长。此外某些溶质组分的过饱 和还会引起一些新矿物的沉淀，这将降低岩矿石的孔隙度甚至堵塞流动通道。 作者简介：谭凯旋(1963—)，男，湖南湘潭人，教授，博士生导师，主要从事地球化学和溶浸采铀的研究 基金项目：国家自然科学基金项目 111 1.2 岩石变形与结构动力学演化 地浸采铀过程中的岩石变形作用相对较弱，主要有矿物溶解引起的压实、孔隙弹性变形、粘塑 性变形、以水力断裂作用为主的断裂变形，这些变形将主要影响含矿层的渗透性，因此考虑岩石变 形对地浸采铀是有意义的。变形的模拟采用递增应力流变学方法计算[6-10] ，应变速率计算公式如式(1): ep up sp fr