陳經(jīng)明， 周澤超， 陳茜茜， 李尋， 羅躍

（1.廣東華南環(huán)保產(chǎn)業(yè)技術(shù)研究院有限公司，廣州 511466；2.核工業(yè)二九〇研究所，廣東 韶關(guān) 512029；3.江西省勘察設(shè)計研究院，南昌 330095；4.東華理工大學，南昌，330013）

開發(fā)鈾資源，需要有高效采鈾工藝技術(shù)。經(jīng)過多年實踐證明，原地浸出采鈾（簡稱地浸采鈾）技術(shù)具有顯著優(yōu)勢，該技術(shù)不需井巷工程、不需要建立礦石、廢石場，生產(chǎn)成本低、對地表生態(tài)環(huán)境破壞小、能適應(yīng)于低品位鈾礦開采，是國內(nèi)外鈾礦采冶的主流技術(shù)之一[1-5]。由于地浸采鈾技術(shù)優(yōu)勢明顯，不適合地浸的南方硬巖型鈾礦，已相繼政策性停產(chǎn)。新疆、內(nèi)蒙古等地相繼發(fā)現(xiàn)多個適合地浸的砂巖型鈾礦后[6]，我國鈾資源采冶已基本采用地浸方式采鈾。因此地浸采鈾已成為我國提高鈾產(chǎn)量的重要技術(shù)，甚至是唯一技術(shù)。

地浸采鈾技術(shù)基本原理是：首先通過注液井向含礦含水層中注入溶浸液，讓其與礦石充分接觸，使鈾礦溶解于地下水中，然后通過抽液井將富含溶解態(tài)鈾的地下水抽出，最后將抽出的地下水（稱為浸出液）輸送至浸出液處理廠提取鈾，尾液將作為溶浸液又注入含水層中，形成地下水抽注循環(huán)（圖1和圖2）。按照溶浸液使用的試劑不同將地浸分為酸法、堿法和中性，其中酸法使用硫酸作為溶浸液（圖2），使用較普遍。

圖1 地浸采鈾原理示意Fig.1 Principle of in-situ leaching uranium

圖2 酸法地浸采鈾過程中地下水抽注循環(huán)示意Fig.2 Schematic diagram of groundwater extraction cycle during acid in-situ leaching uranium

溶浸范圍內(nèi)抽注井網(wǎng)布設(shè)與鉆井生產(chǎn)能力、礦層連續(xù)性及滲透性等因素有關(guān)，目前各個鈾礦抽注孔數(shù)量主要根據(jù)單孔抽注液量及水文地質(zhì)試驗結(jié)果來確定[7-8]。地浸井場鉆孔布設(shè)類型主要有網(wǎng)格式和行列式2種，網(wǎng)格式主要有五點型（四注一抽）和七點型（六注一抽）2種，行列式抽注井各自成行平行排列[9]。

目前國內(nèi)外關(guān)于地浸采鈾模擬的研究較少，大部分學者進行的是現(xiàn)場或者室內(nèi)試驗，在實驗的基礎(chǔ)上建立模型研究。趙春虎等較早提出了將數(shù)值模擬應(yīng)用于地浸采鈾中，利用數(shù)值模擬軟件Visual MODFLOW模擬研究了我國某地浸礦山中水動力場與溶質(zhì)運移過程[10-11]。周義朋等通過結(jié)合現(xiàn)場試驗利用數(shù)值模擬研究了某砂巖型鈾礦床地浸過程中NO3-的運移情況[12]。田亮等模擬研究了某礦山地浸試驗流場SO42-的運移規(guī)律[13]。何智等通過成熟的地下水數(shù)值模擬軟件GMS模擬研究了某鈾礦床中地下水溶質(zhì)運移規(guī)律，同時預(yù)測了未來一段時間內(nèi)該采區(qū)浸出鈾濃度變化[14]。焦友軍等對比模擬研究了地下水含水層中六價鈾吸附反應(yīng)模型[15]。陳茜茜等通過PHT3D對酸法地浸過程中在低含氧情況下鈾水解與遷移的時空演化規(guī)律進行了模擬研究[16]。王兵等利用PHT3D模擬分析了酸法地浸采鈾中不同酸度對鈾礦浸出效率影響，及雜質(zhì)礦物對溶浸液的消耗降低鈾礦浸出率[17]。邱文杰等利用TOUGHREACT模擬研究了地浸采鈾過程中溶浸液與礦層水-巖相互作用及反應(yīng)動力學過程[18]。

大部分學者對于地下水的模擬研究更多的著重點是地下水滲透系數(shù)的模擬。桂勇等模擬研究了單孔原地浸礦的滲流量計算模型，結(jié)果表明計算模型穩(wěn)定有效[19]。張將偉等利用不確定分析構(gòu)建了地下水與地表水耦合模型，確定了滲透系數(shù)等參數(shù)在模型中的靈敏度較高[20]。趙勇勝等模擬研究了低滲透系數(shù)非均質(zhì)含水層中污染物遷移規(guī)律[21]。金磊等通過開發(fā)多孔介質(zhì)滲流耦合模擬平臺，揭示了介質(zhì)滲透性降低的本質(zhì)是滲流與介質(zhì)顆粒骨架耦合作用結(jié)果[22]。馬鵬飛等構(gòu)建的近場多孔介質(zhì)滲透模型，可應(yīng)用于多孔介質(zhì)滲流過程模擬[23]。

巴彥烏拉鈾礦是內(nèi)蒙古開發(fā)的第一個酸法地浸砂巖型鈾礦山，自2015年8月正式投產(chǎn)運營以來，礦層滲透性顯著降低，導致鈾浸出量下降，浸鈾成本增高，資源回收率降低，嚴重阻礙了該礦床的鈾資源開發(fā)。本文以巴彥烏拉鈾礦床礦層堵塞為背景，采用數(shù)值模擬法，計算地浸采鈾過程中礦層滲透系數(shù)的時空演化，預(yù)測堵塞發(fā)生的時間和位置,為礦區(qū)開采過程中解決堵塞等問題提供科學依據(jù)。

1 研究方法

1.1 模擬軟件

TOUGH（Trransport of Unsaturated Groundwater and Heat）是由美國勞倫斯-伯克利國家實驗室基于MULKOM模擬程序開發(fā)的數(shù)值模擬程序；是一個可以模擬裂隙或孔隙介質(zhì)中一維到三維的多組分、多相流及非等溫的水流和熱運移的數(shù)值模擬程序，不同程序及所適用的領(lǐng)域如表1所列。本文所用到的數(shù)值模擬軟件TOUGHREACT是由許天福等基于TOUGH基礎(chǔ)上開發(fā)的可應(yīng)用于裂隙或孔隙介質(zhì)中一維到三維的多組分、多相流及非等溫的水流和熱反應(yīng)運移的數(shù)值模擬程序[25-26]，不同程序及所適用的領(lǐng)域如表2所列。TOUGH系列軟件自開發(fā)以來被應(yīng)用于飽和-非飽和帶地下水科學、環(huán)境影響評價和修復(fù)、地熱開發(fā)工程、CO2地質(zhì)封存、核廢物處置、白云石化、頁巖氣及石油開采等領(lǐng)域[27-34]。

表1 TOUGH家族代碼的發(fā)展歷史（自施小清（2009）修[24]）Table 1 Development history of TOUGH family codes(repaired from Shi(2009)[24])

表2 TOUGHREACT模塊功能（自XU（2004）修[35]）Tab le 2 TOUGHREACT module function(repaired from XU(2004)[35])

1.2 控制方程

本文相應(yīng)控制方程均引TOUGHREACT模擬軟件使用說明[25]。在TOUGHREACT中的質(zhì)量守恒示意圖如圖3。

圖3 質(zhì)量（能量）守恒示意Fig.3 Conservation of mass(energy)

質(zhì)量守恒方程可以寫成以下形式：

式（1）中：Vn為 研究區(qū)內(nèi)的任何積分子區(qū)域；Гn為積分區(qū)域表面,m2；M為單位體積的能量或質(zhì)量；k表示不同的組分，如：水、氣、溶質(zhì)等；F為質(zhì)量(能量)對流量,m3/s；q表示源匯相；n表示通過表面積為dΓn的法向量。

式（2）中：k表示組分；Φ表示孔隙度；Sβ表示β相的飽和度；ρβ表示β相的密度,kg/m3；Xkβ表示在相β中組分k的質(zhì)量分數(shù)；β表示液相。

式（5）中：μβ表示β相滲流速度，m/s；k表示絕對滲透率，m/s；krβ表示相對滲透率,m/s；μβ表示液體黏滯系數(shù),Pa；g表示重力加速度常量，g/s2。

式（6）中：Pβ表示液相壓強，Pa；P為參考系壓強，Pa；Pcβ為毛管壓強，Pa。

所采用的反應(yīng)速率表達式來自Lasaga等的研究成果[36]。

式（7）中：n表示礦物;rn表示反應(yīng)速率，（mol/（L·s-1）），正表示溶解，負表示沉淀;kn表示隨溫度變化的速率常數(shù)（mol/（m2·s-1））；An表示每千克水的礦物比表面積（cm2/g）；Ωn表示礦物飽和比率；θ與η由實驗得出，通常取1。

式（8）中：k25表示25℃的速率常數(shù)，mol/（m2·s）；Ea為活化能，kJ/mol；R為氣體常數(shù)，kJ/（mol·K），通常取8.314 51；T表示絕對溫度，K。

滲透率k（Permeability）表征巖層對不同流體的固有滲透性能，滲透率僅取決于巖石的空隙性質(zhì)，與滲流的液體性質(zhì)無關(guān)。滲透率與滲透系數(shù)K的關(guān)系為[37]：

式（9）中:ρ為液體密度，kg/m3;g為重力加速度,m/s2;μ為液體動力黏滯系數(shù),N·s/m2;k為滲透率,常用單位為μm2（即達西）或cm2。

在進行多場耦合建模時需要給定不同物理參數(shù)之間的關(guān)系，如孔隙度-滲透率之間的關(guān)系[38]?？紫抖?滲透率之間復(fù)雜的關(guān)系取決于孔徑分布、孔隙形狀及孔隙間的連通性等因素的相互作用，無法僅僅通過簡單的立方根定理和Konzeny-Carman孔隙度-滲透率方程反映[39]。有室內(nèi)實驗表明，產(chǎn)生適量的沉淀物會導致孔隙度減少且滲透系數(shù)會大幅降低[40]?，F(xiàn)場實驗表明介質(zhì)滲透率對于孔喉變化較為敏感[41]，這與在天然情況下滲流通道的收斂-發(fā)散性一致，由于孔喉被沉淀物堵塞，導致介質(zhì)中存在不連續(xù)的孔隙空間[39]。滲透率降低不僅取決于總孔隙度的減少，同時也取決于孔隙空間幾何形狀和沉淀物的空間分布有關(guān)[42]。模型中采用VERNA等改進的孔隙度-滲透率關(guān)系式[39]。

Verman—Pruess模型：

式（10）中:ki和φi分別表示初始條件的滲透率和孔隙度；φc表示臨界孔隙度，當孔隙度達到這個值時滲透率為0；n表示冪指數(shù)；φc和n的取值由介質(zhì)的性質(zhì)決定。

根據(jù)式（9）、式（10）可得滲透系數(shù)與孔隙度之間的關(guān)系，如式（11）：

本文中孔隙度初始值取0.2；前人的研究認為當多孔介質(zhì)中孔隙度約20%時，滲透系數(shù)會降低幾個數(shù)量級[43]，因此φc臨界孔隙度取初始值的80%即0.16；滲透率為8.03×10-12m2；由此可確定模型中滲透系數(shù)（K，m/d）與孔隙度的關(guān)系式為：

滲透系數(shù)與孔隙度關(guān)系曲線見圖4。

圖4 滲透系數(shù)—孔隙度關(guān)系Fig.4 Relationship between permeability coefficient and porosity

2 數(shù)值模擬

2.1 模型概化

基于內(nèi)蒙古某砂巖型鈾礦的多個抽注系統(tǒng)，模擬酸法地浸采鈾過程中出現(xiàn)滲透系數(shù)減小的現(xiàn)象。該場地的地下水動力場以及地浸采鈾的堵塞物質(zhì)已有詳細研究[44-45]。模型共有6個注液孔（Z1-Z6），2個抽液孔（C1、C2）（圖5），構(gòu)成兩組五點型二維地浸抽注系統(tǒng)；在確定模擬區(qū)域時分別模擬了邊長為2 000、1 000、500、200 m正方形區(qū)域，模擬結(jié)果顯示溶浸液的滲流范圍較?。▓D6模擬結(jié)果流場圖），為了減少模擬時所需的計算量最終確定將邊界大小為200 m較為合適。

圖5 二維模型概念示意Fig.5 Schematic diagram of the 2-D model

圖6 二維流場矢量示意Fig.6 Vector diagram of two-bit flow field

2.2 模型條件

本研究的模擬條件來自鄭和秋野等對巴彥烏拉鈾礦現(xiàn)場地浸的條件試驗和地浸場水動力模擬。模型邊界大小確定為200 m，含礦含水層厚度為58.3 m，每個注液孔與抽液孔的距離為30 m，每個注液孔的注入流量為4 m3/h，抽液孔抽出流量為12 m3/h；含礦層厚度58.3 m，含礦層的初始孔隙度假設(shè)為0.2，初始滲透系數(shù)為7.6 m/d（滲透率為8.03×10-12m2）,總模擬時間為690 d。在地浸采鈾中將尾液中的有用元素（鈾）沉淀之后，再次作為溶浸液注入注液孔，這種工藝也稱為循環(huán)抽注系統(tǒng)。這種工藝容易造成礦層滲透系數(shù)減?。ㄉ踔炼氯蚴腔毓嗟降V層中的地浸尾液礦化度較高。通過對區(qū)域地浸尾液室內(nèi)試驗數(shù)據(jù)的分析得出：研究區(qū)域地浸開采100 d左右后，浸出液中的主要離子濃度趨于穩(wěn)定（尤其是Ca2+、Mg2+），因此模型中的邊界條件各類離子濃度采用100 d后的平均值。初始條件與邊界條件的各離子濃度列于表3。根據(jù)已有的資料整理得到模擬區(qū)域雜質(zhì)礦物含量列于表4，用于計算礦物動力學速率常數(shù)的參數(shù)列于表5。

表3 模型初始/邊界條件Table 3 Model initial/boundary conditions

表4 模擬區(qū)域雜質(zhì)礦物含量Tab l e 4 Impurity mineral content in the simulated area

表5 用于計算礦物動力學速率常數(shù)的參數(shù)Tab le 5 Parameters for calculating the kinetic rate contents of minerals

3 模擬結(jié)果與討論

基于模型概化，選取模擬結(jié)果中各因子變化較明顯時刻，進行討論分析；對模擬結(jié)果中滲透系數(shù)（K）的時空演化規(guī)律進行討論分析。模擬結(jié)果中各因子變化較明顯的時間分別為：6.9、41.4、200、283、469、690 d；同時選取剖面A、B對滲透系數(shù)的空間演化結(jié)果進行討論分析。

3.1 不同時刻滲透系數(shù)演化結(jié)果

在t=6.9 d時，模擬區(qū)域內(nèi)的滲透系數(shù)均大于初始值（7.6 m/d），滲透系數(shù)介于7.6~7.69 m/d之間，變化量較小，相對初始值最大僅增大1.2%。滲透系數(shù)最大值均位于注液孔附近，與注液孔距離約4 m。在抽液孔附近的滲透系數(shù)幾乎不發(fā)生變化。從圖7中可以看出滲透系數(shù)發(fā)生變化的區(qū)域也主要集中在注液孔附近，最大變化半徑約為10 m，且變化方向朝著注液孔一側(cè)。

圖7 不同時刻滲透系數(shù)演化結(jié)果Fig.7 Evolution results of permeability coefficient at different times

在溶浸液（H2SO4）反應(yīng)前鋒作用下與注液孔距離越大，滲透系數(shù)越接近初始值，不同注液孔周圍的滲透系數(shù)演化情況相似，沒有井群干擾現(xiàn)象。

與上一個時刻（6.9 d）對比，在t=41.1 d時模擬區(qū)域內(nèi)滲透系數(shù)發(fā)生變化的面積有所增大，且該時刻開始出現(xiàn)滲透系數(shù)小于初始值的網(wǎng)格點，以注液孔為圓心滲透系數(shù)發(fā)生變化區(qū)域最大半徑為22.63 m。

根據(jù)圖7可知，t=41.1 d時模擬區(qū)域滲透系數(shù)介于6.65~7.75 m/d之間；滲透系數(shù)最大值和最小值較初始值分別增大2.0%和-12.5%。

說明隨著時間的推移模擬區(qū)域內(nèi)的滲透系數(shù)區(qū)間會不斷擴大，并且主要的變化是滲透系數(shù)最小值會不斷減小，從圖7中可以看出，在滲透系數(shù)改變的區(qū)域內(nèi)，主要表現(xiàn)為滲透系數(shù)不斷減小，僅有少量緊鄰注液孔的位置滲透系數(shù)增大。

根據(jù)圖7給出的t=200 d模擬區(qū)域滲透系數(shù)演化結(jié)果可知，在該時刻模擬區(qū)域內(nèi)滲透系數(shù)介于5.99~9.08 m/d之間，相對初始值最小值減少21.14%，最大值增大19.54%。

從圖7中可以看出，模擬區(qū)域內(nèi)滲透系數(shù)增大的范圍僅出現(xiàn)在靠近注液孔周圍的極小區(qū)域內(nèi)（半徑約為4 m），沿溶浸液滲流范圍內(nèi)的其他區(qū)域均出現(xiàn)了滲透系數(shù)減小現(xiàn)象。通過與之前的幾個時刻（6.9 d、41.4 d）對比，說明在模擬區(qū)域內(nèi)靠近注液孔處的礦層隨著時間的延長不斷與溶浸液發(fā)生溶解反應(yīng)，且在滲流作用和反應(yīng)動力學共同作用下，在滲流路徑中不斷會有沉淀物生成，同時也有溶解反應(yīng)發(fā)生。

根據(jù)模擬結(jié)果顯示，該時刻滲透系數(shù)小于初始值的大部分礦層內(nèi)滲透系數(shù)均介于5.99~6.77 m/d之間，少部分區(qū)域滲透系數(shù)介于7.15~7.45 m/d之間。

從模擬結(jié)果可知在抽液孔附近（半徑小于2 m）滲透系數(shù)已經(jīng)開始發(fā)生變化。說明模擬到200 d時，溶浸液與雜質(zhì)礦物發(fā)生的沉淀-溶解反應(yīng)已經(jīng)開始影響抽液周圍的滲透系數(shù)。且與注液孔距離越遠滲透系數(shù)減少值則越小，說明在滲流通道上沉淀-溶解量不斷發(fā)生變化；與注液孔距離近的區(qū)域沉淀物的量較多，溶解量較少；靠近抽液孔處沉淀量相對減少，溶解量也相對增加。

圖7給出了t=283 d模擬區(qū)域滲透系數(shù)演化結(jié)果，t=283d時模擬區(qū)域內(nèi)滲透系數(shù)區(qū)間為5.88~9.44 m/d。最大值相對初始值增大了30.34%,最小值相對初始值減小了22.61%。對比前文中討論的幾個時刻，該時刻最大的區(qū)別是模擬區(qū)域內(nèi)滲透系數(shù)的最大值變化量比最小值變化量大。

滲透系數(shù)最大值（9.94 m/d）均出現(xiàn)在注液孔附近，且滲透系數(shù)發(fā)生變化的范圍幾乎沒有變化；從模擬結(jié)果圖可以看出該區(qū)域內(nèi)滲透系數(shù)均大于7.9 m/d。說明在注液孔附近滲透系數(shù)增大的區(qū)域內(nèi)，不斷發(fā)生溶解反應(yīng)，產(chǎn)生的離子隨滲流通斷不斷運移；在該區(qū)域內(nèi)礦物溶解量遠大于沉淀量。

從圖7可知該時刻滲透系數(shù)發(fā)生改變的區(qū)域內(nèi)只有少部分區(qū)域滲透系數(shù)區(qū)間為[7.4 m/d,7.9 m/d]，說明在模擬達到283 d時，礦層內(nèi)發(fā)生沉淀和溶解的區(qū)域分界線較明顯。

在滲透系數(shù)發(fā)生改變的區(qū)域中滲透系數(shù)介于區(qū)間為[5.88 m/d,6.9 m/d]的區(qū)域占比較大。說明在這些區(qū)域內(nèi)，反應(yīng)生成的沉淀物量大于溶解量；且該區(qū)域內(nèi)大部分位置滲透系數(shù)均小于6.39 m/d，只有少量靠近抽液孔處區(qū)域的滲透系數(shù)為6.9 m/d。

隨著模擬時間不斷增加，滲流場-化學場耦合作用對模擬區(qū)域內(nèi)的滲透系數(shù)以及其他因子影響越明顯。

圖7給出496 d時的滲透系數(shù)演化結(jié)果，該時刻礦層內(nèi)滲透系數(shù)最大值為12.01 m/d，相對初始值增大了57.99%；滲透系數(shù)最小值為5.65 m/d，相對初始值減小了25.71%。

在前文中討論到283 d時模擬區(qū)域滲透系數(shù)在兩組五點型抽注系統(tǒng)的影響范圍有公共區(qū)域，從模擬結(jié)果看出，該時刻（496 d）滲透系數(shù)在模擬區(qū)域內(nèi)分布較均一。在滲透系數(shù)減小的區(qū)域內(nèi)，不同抽注液孔之間幾乎沒有差異；兩組抽注系統(tǒng)的滲透系數(shù)空間分布關(guān)于直線X=100對稱，這是之前幾個時刻沒有出現(xiàn)的情況；說明經(jīng)過496 d的反應(yīng)運移，滲流場-化學場相互影響已經(jīng)開始趨于穩(wěn)定狀態(tài)，至少模擬區(qū)域內(nèi)的沉淀物空間分布差異性開始減小。

對于區(qū)域內(nèi)滲透系數(shù)增大和減小的區(qū)域，與前文的分析相似。區(qū)別在于該時刻礦層內(nèi)滲透系數(shù)減小的范圍有所增大，這是由于不斷注入的溶浸液與伴生礦物反應(yīng)增加了沉淀物的量，同時也由于反應(yīng)動力學的作用，沉淀物會在滲流路徑中生成，同時可以解釋，礦層內(nèi)滲透系數(shù)小于初始值的區(qū)域內(nèi)滲透系數(shù)差異性較小。

根據(jù)圖7給出的模擬結(jié)束時（690 d）的滲透系數(shù)演化結(jié)果，該時刻礦層內(nèi)滲透系數(shù)最大值為14.68 m/d，相對初始值增大了93.18%；滲透系數(shù)最小值為5.37 m/d，相對初始值減小了29.31%。

隨著時間的推移，礦層內(nèi)滲透系數(shù)改變的區(qū)域不斷增大。從模擬結(jié)果可以看出，模擬區(qū)域內(nèi)滲透系數(shù)發(fā)生改變的區(qū)域主要集中在抽注液孔之間，在注液孔外側(cè)也有部分區(qū)域滲透系數(shù)減小。從圖7可知，注液孔Z1、Z3、Z4、Z6演化情況一致，Z2、Z5演化情況一致。

滲透系數(shù)區(qū)間值不斷增大，模擬結(jié)束時（690 d）滲透系數(shù)最大值（14.68 m/d）是最小值（5.37 m/d）的2.73倍，且滲透系數(shù)增大的區(qū)域除了集中在注液孔附近，在外側(cè)也有少量位置滲透系數(shù)所有增大；滲透系數(shù)減小的區(qū)域不斷增大。由于不斷注入的溶浸液與雜質(zhì)礦物發(fā)生反應(yīng)，在注液孔處大量的雜質(zhì)礦物發(fā)生溶解反應(yīng)，但是在反應(yīng)動力學的作用下沉淀反應(yīng)不會馬上發(fā)生，因此沉淀物會分布在滲流通道上。對于滲透系數(shù)減小的區(qū)域內(nèi)，也會有溶解反應(yīng)生成，但是模擬結(jié)果顯示生成的沉淀物量大于溶解量。

對比前文中討論的幾個時刻，模擬結(jié)束時最大的特點是在滲透系數(shù)減少的區(qū)域外側(cè)，有少量區(qū)域的滲透系數(shù)相對初始值有所增大（8.86 m/d）。分析可能的原因是這些區(qū)域也是滲流場影響的范圍內(nèi)，運移到這些網(wǎng)格的溶浸液內(nèi)還有少量的H+未反應(yīng)，與礦層發(fā)生反應(yīng)導致滲透系數(shù)增大。

圖7給出了不同時刻（6.9、41.4、200、690 d）二維模型滲透系數(shù)演化結(jié)果。從演化結(jié)果來看，滲透系數(shù)發(fā)生變化的區(qū)域隨時間不斷增大，模擬區(qū)域內(nèi)滲透系數(shù)范圍隨時間不斷增大。

在滲流作用和化學反應(yīng)前鋒不斷向抽液孔推進，井群干擾現(xiàn)象越明顯，6個注液孔周圍的滲透系數(shù)演化情況略有差異。Z1、Z3、Z4、Z6 4個注液孔演化情況相似，從圖中可以看出滲透系數(shù)變化主要方向不斷向C1或C2抽液孔推進；Z2、Z5 2個注液孔演化情況相似，從圖7中可以看出滲透系數(shù)變化主要方向不斷同時向C1和C2抽液孔推進。出現(xiàn)這種情況的原因是C1、C2共同影響Z2、Z5注液孔周圍的流場分布，其他4個注液孔周圍的流場僅受C1或者C2單獨影響。

3.2 不同空間滲透系數(shù)演化結(jié)果

圖8給出了兩個剖面（A、B）滲透系數(shù)在t=6.9 d、t=41.4 d、t=200 d、t=283 d、t=469 d及t=690 d時的演化結(jié)果。0 m表示注液孔位置，30 m表示抽液孔位置。

從圖8中可以看出在A、B兩個剖面中，從注液孔到抽液孔之間的滲透系數(shù)演化情況相同之處是：滲透系數(shù)隨距離的總體趨勢均是不斷減?。辉赿=17.23 m到注液孔處兩組抽注孔之間的滲透系數(shù)演化結(jié)果相同。

圖8 不同空間位置滲透系數(shù)演化結(jié)果Fig.8 Evolution results of permeability coefficients at different spatial locations

不同之處：A和B兩個剖面上滲透系數(shù)在該時刻演化結(jié)果有所差異。A剖面上滲透系數(shù)區(qū)間為[7.60 m/d，7.68 m/d]，B剖面上滲透系數(shù)區(qū)間為[7.60 m/d,7.69 m/d]。

在d=17.23 m之前兩組抽注孔之間的滲透系數(shù)演化雖然是沿滲流路徑不斷減小，但是相同距離下的滲透系數(shù)均不相同；在d=17.23 m之前的滲流路徑內(nèi)大部分（占抽注液孔間距的81.7%）位置中Z2-C1之間的滲透系數(shù)大于Z4-C1之間的滲透系數(shù)。

根據(jù)圖8給出的兩個剖面（A、B）滲透系數(shù)在t=41.4 d時演化結(jié)果。對比2個剖面滲透系數(shù)演化結(jié)果，相同點是：兩個剖面上滲透系數(shù)沿滲流路徑的演化趨勢相同，均是先減小后增大；兩個剖面上該時刻滲透系數(shù)最小值出現(xiàn)的位置相同（均是于注液孔相距5.58 m處）；兩個剖面在14.82 m之后的空間位置上滲透系數(shù)差異較小。

在A、B兩個剖面上從注液孔到滲透系數(shù)最小值處（5.63 m處）之間區(qū)域內(nèi)也有所差異，A剖面上滲透系數(shù)沿滲流路徑不斷減小，B剖面從注液孔（7.50 m/d）到2.83 m（7.53 m/d）之間有所增大且均小于初始值。

兩個剖面上滲透系數(shù)演化的不同點：兩個剖面上滲透系數(shù)最小值有所不同，其中A剖面上滲透系數(shù)最小值為7.02 m/d，B剖面上滲透系數(shù)最小值為6.99 m/d。

A剖面在靠近注液孔處的滲透系數(shù)大于初始值（7.75 m/d），在B剖面上靠近注液孔處的滲透系數(shù)小于初始值（7.50 m/d）。

兩個剖面從注液孔到5.88 m處兩個剖面上滲透系數(shù)差異性較大，主要體現(xiàn)在B剖面不同位置滲透系數(shù)演化較復(fù)雜，從注液孔處出發(fā)有少量位置的滲透系數(shù)相對注液孔出有所增大；從5.58 m至抽液孔處滲透系數(shù)演化也有所差異，B剖面大部分空間位置的滲透系數(shù)均小于A剖面相同位置。

根據(jù)圖8給出的兩個剖面（A、B）滲透系數(shù)在t=200 d時演化結(jié)果。A、B兩個剖面滲透系數(shù)最小值均出現(xiàn)在5.66 m處，該時刻兩個剖面上滲透系數(shù)最小值差異較小（A剖面滲透系數(shù)最小值為5.99 m/d，B剖面滲透系數(shù)最小值為6.05 m/d，相差1%）；兩個剖面上滲透系數(shù)演化結(jié)果均出現(xiàn)了大于初始值和小于初始值的空間位置。

A和B兩個剖面上滲透系數(shù)在該時刻演化結(jié)果有所差異。A剖面上滲透系數(shù)區(qū)間為[5.99 m/d，9.08 m/d]，B剖面上滲透系數(shù)區(qū)間為[6.05 m/d，8.41 m/d]，B剖面上出現(xiàn)了滲透系數(shù)明顯大于初始值的位置，A剖面上滲透系數(shù)區(qū)間明顯大于B剖面上滲透系數(shù)區(qū)間，說明B剖面上的滲透系數(shù)較A剖面滲透系數(shù)演化較慢。

模擬結(jié)果看出在該時刻A、B剖面大部分位置處，B剖面上的滲透系數(shù)率大于A剖面上的滲透系數(shù)（注液孔除外）；14.15~30 m之間的區(qū)域內(nèi)，B剖面上滲透系數(shù)值曲線一直在A剖面滲透系數(shù)值曲線上方；在5.63~14.15 m之間存在少部分區(qū)域（約2.81 m）B剖面上滲透系數(shù)值曲線一直在A剖面滲透系數(shù)值曲線下方。

圖8給出的兩個剖面（A、B）滲透系數(shù)在t=283 d時演化結(jié)果。A、B兩個剖面滲透系數(shù)最小值均出現(xiàn)在5.66 m處，該時刻兩個剖面上滲透系數(shù)最小值差異較?。ˋ剖面滲透系數(shù)最小值為5.88 m/d，B剖面滲透系數(shù)最小值為6.94 m/d，相差1%）。

A和B兩個剖面上滲透系數(shù)在該時刻演化結(jié)果有所差異。A剖面上滲透系數(shù)區(qū)間為[5.88 m/d,9.94 m/d]，B剖面上滲透系數(shù)區(qū)間為[5.94 m/d,9.17 m/d]，B剖面上出現(xiàn)了滲透系數(shù)明顯大于初始值的位置，A剖面上滲透系數(shù)區(qū)間明顯大于B剖面上滲透系數(shù)區(qū)間，說明B剖面上的滲透系數(shù)較A剖面滲透系數(shù)演化較慢；在該時刻，兩個剖面在滲透系數(shù)較小的區(qū)域內(nèi)（5.66~16 m）之間的演化情況一致。

兩個剖面滲透系數(shù)演化的差異性體現(xiàn)在靠近注液孔處，以及16~30 m這兩段區(qū)域內(nèi)；在16~30 m之間，B剖面滲透系數(shù)均大于A剖面上相同位置處的滲透系數(shù)，注液孔到5.66 m之間兩個剖面上滲透系數(shù)的差異與前文中討論一致，就不贅述。

圖8給出的兩個剖面（A、B）滲透系數(shù)在t=469 d時演化結(jié)果。A、B兩個剖面滲透系數(shù)最小值均出現(xiàn)在5.66 m處，該時刻兩個剖面上滲透系數(shù)最小值差異較?。ˋ剖面滲透系數(shù)最小值為5.88 m/d，B剖面滲透系數(shù)最小值為6.94 m/d，相差0.8%）；在該時刻，兩個剖面從滲透系數(shù)最小值到抽液孔處（5.66~30 m）之間的演化情況一致。

A和B兩個剖面上滲透系數(shù)在該時刻演化結(jié)果有所差異。A剖面上滲透系數(shù)區(qū)間為[5.65 m/d,12 m/d]，B剖面上滲透系數(shù)區(qū)間為[5.70 m/d,10.94 m/d]，B剖面上出現(xiàn)了滲透系數(shù)明顯大于初始值的位置，A剖面上滲透系數(shù)區(qū)間明顯大于B剖面上滲透系數(shù)區(qū)間，說明B剖面上的滲透系數(shù)較A剖面滲透系數(shù)演化較慢。

通過兩個具有代表性的剖面滲透系數(shù)演化結(jié)果可知，在該時刻，模擬區(qū)域內(nèi)大部分空間位置的滲透系數(shù)均小于初始值。

根據(jù)圖8給出的兩個剖面（A、B）滲透系數(shù)在t=690 d（模擬結(jié)束）時的演化結(jié)果。A、B兩個剖面滲透系數(shù)最小值均出現(xiàn)在5.66 m處，該時刻兩個剖面上滲透系數(shù)最小值差異較小（A剖面滲透系數(shù)最小值為5.37 m/d，B剖面滲透系數(shù)最小值為5.42 m/d，相差1.11%）；在該時刻，兩個剖面在滲透系數(shù)較小的從滲透系數(shù)最小值到抽液孔處（5.66~30 m）之間的演化情況一致。

A和B兩個剖面上滲透系數(shù)在該時刻演化結(jié)果有所差異。A剖面上滲透系數(shù)區(qū)間為[5.37 m/d，14.68 m/d]，B剖面上滲透系數(shù)區(qū)間為[5.42 m/d，13.20 m/d]，B剖面上出現(xiàn)了滲透系數(shù)明顯大于初始值的位置，A剖面上滲透系數(shù)區(qū)間明顯大于B剖面上滲透系數(shù)區(qū)間，說明B剖面上的滲透系數(shù)較A剖面滲透系數(shù)演化較慢。

在該時刻，模擬區(qū)域內(nèi)僅有注液孔出的滲透系數(shù)大于初始值，且從模擬結(jié)果看出，在溶液滲流通道內(nèi)大部分區(qū)域的滲透系數(shù)均在6 m/d左右。

通過對比模擬區(qū)域內(nèi)兩個剖面內(nèi)滲透系數(shù)演化的差異可知，在多組抽注系統(tǒng)共同作用下，不同抽注井之間滲透系數(shù)演化差異性隨時間推移不斷減小。在文中給的抽注系統(tǒng)中，滲透系數(shù)最小值出現(xiàn)在與注液孔相距大約5.66 m處；受單個抽液孔影響的注液孔周圍滲透系數(shù)在滲流路徑上演化較為單一，受多個（文中僅有兩個）抽液孔影響的注液孔附近滲透系數(shù)在模擬初期演化較為多樣性；

隨著模擬時間的推移，溶浸液滲流通內(nèi)滲透系數(shù)會不斷較小，且區(qū)域內(nèi)滲透系數(shù)會區(qū)域一致。

從選取的兩個特征剖面來看，不同抽注系統(tǒng)內(nèi)滲透系數(shù)存在一定的相互影響的現(xiàn)象，隨著時間的推移影響越明顯，主要體現(xiàn)在模擬末期兩個剖面上的滲透系數(shù)演化趨于一致，差異性較小。

4 結(jié)論與展望

本文基于前人的研究基礎(chǔ)，運用數(shù)值模擬的方法，對某砂巖型鈾礦在原位地浸采鈾過程中滲透系數(shù)的時空分布規(guī)律進行了模擬研究。通過討論分析反應(yīng)動力學參數(shù)、地下水滲流場對滲透系數(shù)及溶質(zhì)運移的影響，得到了以下主要結(jié)論：

模擬區(qū)域內(nèi)滲透系數(shù)區(qū)間隨時間推移增大，模擬結(jié)束時滲透系數(shù)最大值為14.7 m/d，最小值為5.37 m/d。從滲透系數(shù)空間演化規(guī)律可知滲透系數(shù)最大值出現(xiàn)在注液孔處，滲透系數(shù)最小值出現(xiàn)在與注液孔中心相距約5.66 m處的圓周上；模擬結(jié)果表明，在模擬初期不同抽注系統(tǒng)之間溶質(zhì)運移受井群干擾作用較明顯。

本文基于參考文獻[44-45]的研究基礎(chǔ)，建立多場耦合模型對內(nèi)蒙古某砂巖型鈾礦開采中滲透系數(shù)進行模擬研究，對提升鈾礦開采效率具有一定的科學意義。

由于水文地質(zhì)單元概化及模擬程序限制等因素的影響，本文在酸法地浸采鈾過程的模擬研究中未對介質(zhì)的非均一性、微生物過程、溶浸液中不同離子相互作用及溶浸液循環(huán)抽注等因素對結(jié)果影響。