基于CFD的鈾尾礦庫分層覆土中氡的運移規律研究

吳慧 鄒樹梁

摘要：參照低滲透的鈾尾礦庫的分層覆土設計圖，建立分層覆土模型。根據多孔介質模型修正基本控制方程，并基于CFD建立湍流模型和多組分輸運模型，結合多層覆土孔隙度、滲透系數、氡衰變系數等因素，對鈾尾礦庫氡在多層覆土中的運移規律進行模擬，模擬結果表明采用合理有效的多層覆土的治理方式較之均勻覆土能夠有效減少氡運移速度。在氡入口通量為24Bq/m2下，分層覆土可有效降低灘面0.02307mol/m3的氡析出濃度。并采用tecplot后處理軟件展示可視化流場模擬結果，探討氡在分層覆土中的運移規律，將對運行中或退役后的鈾尾礦庫的治理提供必要的技術指導和理論支持。

關鍵詞：CFD;多孔介質;分層覆土;湍流

中圖分類號：X13 文獻標識碼：A 文章編號：2095-672X（2020）03-0-04

DOI：10.16647/j.cnki.cn15-1369/X.2020.03.074

Study on radon migration in layered overburden of uranium tailings pond based on CFD

Wu Hui1，Zou Shuliang2

（1.School of Environment and Safety Engineering，University of South China，Hengyang Hunan 421001，China;

2.Hunan Provincial Key Laboratory of Emergency Saftey Technology and Equipment for Nuclear Facilities，Hengyang Hunan 421001，China）

Abstract：In this paper， a layered covering model is established by referring to the layered covering design of low permeability uranium tailings reservoir.The basic governing equation is modified according to the porous medium model， and the turbulence model and multi-component transport model are established based on CFD.The simulation results show that the radon migration rate can be effectively reduced by adopting a reasonable and effective multi-layer covering soil treatment method compared with uniform covering soil.When the entrance flux of radon is 24Bq/m3， the radon concentration of 0.02307mol/m3 can be effectively reduced by layered soil covering. In addition，Tecplot post-processing software is used to show the simulation results of visualized flow field， and to explore the migration rule of radon in layered overburden， which will provide necessary technical guidance and theoretical support for the treatment of uranium tailings pond in operation or after decommissioning.

Key words：CFD;Porous media;Layered soil;Turbulent flow

鈾尾礦庫中堆放大量開采冶煉過程中產生的廢石和尾礦，存在放射性危害，其中氣態放射性核素氡及子體屬于易于擴散的惰性氣體，其對生態及環境具有嚴重危害。鈾尾礦庫通常采用分層覆土覆蓋的治理方式以減少氡析出率，使灘面氡析出率不超過管理限值0.74 Bq/（m2·s）[1]。

氡在覆土中的運移是典型的多孔介質流動行為，多孔介質表面的氡析出則是氡運移到介質表面后析出到外界環境的過程[2-3]。國內外學者對氡的運移過程大多是針對土壤、鈾尾礦等常見的多孔介質內的氣相氡遷移[4-6]。2003年，Kozak[7]等建立了一個描述尾礦庫灘面土壤中氣液兩相流的一維模型，描述了氡隨著這兩種流體的遷移行為。2010年，Yakovleva S[8]通過數值和實驗方法獲得的氡場時間序列，并給出了多層地質介質中氡輸運的靜態和非平穩擴散-平流方程的解。國內最早對尾礦庫析出方面進行系統研究是在1995年劉慶成[9]利用試驗的方法建立了氡運移的一維模型模擬了氡在尾礦壩沙土中的運移過程;2014年，武劍波[10]基于氡運移二維方程和邊界條件理論建立了內部含裂系、裂隙系-斷層帶和非均勻三種覆蓋層物理模型，模擬了覆蓋層中氡運移過程。

氡在鈾尾礦庫分層覆土中的運移受孔隙度、滲透系數、氡衰變系數等多因素的影響，多孔介質流動行為更加復雜。近年來，基于CFD數值模擬方法能夠解決更為復雜的多孔介質流動問題[11-13]，被廣泛應用于復雜多孔介質流動問題的研究中，可視化的流場模擬使實驗效果擺脫了僅僅依靠理論分析的層面。

因此，本文參照低滲透的鈾尾礦庫的分層覆土設計理念建立多層覆土模型，依據多孔介質模型修正基本控制方程，并基于CFD建立湍流模型和多組分輸運模型，結合多層覆土孔隙度、滲透系數、氡衰變系數等因素，對鈾尾礦庫氡在多層覆土中的運移規律進行模擬。進一步采用tecplot后處理軟件展示可視化流場模擬結果，從定性和定量上探討氡在覆土中運移規律，對鈾尾礦庫的治理具有重要意義。

1 運移擴散模型建立

鈾尾礦庫中含有大量放射性核素，采用多層覆土覆蓋，每一層覆土簡化為各向同性、孔隙均勻的多孔介質模型，氡在多孔介質的運移為不可壓流動，采用三維定常不可壓縮方程描述氡在多孔介質中運移的主控方程包括以下質量守恒方程、動量守恒方程、氡濃度輸運方程? [14-15]。

質量守恒方程：

動量方程：

其中，為新的源項，多孔介質動量方程增加動量源項，包含兩部分，一部分是粘性阻力，另一部分是慣性阻力。氡濃度輸運方程[15]：

式中μ表示氡氣在多孔介質中的動力粘度，N s/m2;ρ表示流體的密度，kg/m3;α表示多孔介質的滲透率，m2;D表示氡氣在介質中運移的擴散系數，m2/s;λ表示氡的衰變系數;Va表示氡的體積濃度，Bq/m3。

2 基于CFD數值模擬

2.1 物理模型和網格劃分

鈾尾礦庫灘面覆蓋層的設計可以最大限度的減少鈾尾砂堆積層的異常氡析出對環境造成的有害影響。通常鈾尾礦庫的灘面控氡設計一般是由多層覆蓋組成的，表1為低滲透的鈾尾礦庫的分層覆土示意圖[17]，大體可分為臨時覆蓋層、排水層、隔離沙土層、儲存層、蓄水植被覆蓋層。其中隔離沙土層、儲存層、蓄水植被覆蓋層滲透率相同。

所以本文將覆土模型簡化為4層（1m1m4m），設置4個body。假設每層覆土孔隙均勻，土壤條件各向同性。從下至上第一層為1.5m，孔隙度為0.35，擴散系數為110-6m2/s;第二層為0.5m密封層，孔隙度為0.1，擴散系數為0.510-6m2/s;第三層為0.5m排水層，孔隙度為0.75，擴散系數為310-6m2/s;第四層為1.5m，孔隙度為0.5，擴散系數為210-6m2/s，采用非結構化網格，覆土交界面網格進行加密處理，第一層網格數為103518，第二層網格數為155460，第三層網格數為155231，第四層網格數為166331，如圖1所示。

2.2 邊界條件和離散方式

根據文獻資料，我國裸露鈾尾礦灘面氡析出的通量為1.65-26.5Bq/m2s，1Bq的氡的質量為1.7510-19kg[18]，所以放射性核素氡底面釋放入口邊界條件設置為質量流入口，定義為Mass flow inlet 類型，土壤四面邊定義為symmetry 對稱面條件類型，頂部地表面出口邊界定義為 pressure-outlet 出口條件類型。

氡的衰變系數λ為2.09710-6;氡氣在覆土的運動粘性系數取μ為1.810-5m2/s;慣性阻力設置為為0.01396426，粘性阻力設置為67390。選用基于壓力基的求解方式（Pressure-based）進行穩態計算和瞬態計算。算法采用SIMPLIC算法，中心項采用二階迎風格式。

3 結果與討論

3.1 分層覆土對運移速度的影響

在鈾尾礦庫分層覆土中（假設每層覆土各向同性），孔隙度對流體運移速度有一定的影響。由圖2、圖3可知， 在第一層和第二層覆中，孔隙度較小，運移速度相應較小;在第三層覆土中，孔隙度最大，運移速度相應變大。在每層覆土的中心處運移速度增大，在不同材料覆土的交界面處，速率變化明顯降低，且隨著運移高度增加，在進入到上層覆土中擴散速率加快，達到最大值。入口處氡通量影響氡在覆土中的運移速度，由圖3可知，入口氡通量越大，氡運移速度增大。對比鈾尾礦庫分層覆土治理方式和均勻一層覆土（空隙度0.5）治理方式中氡的運移速度，在均勻空隙度0.5的覆土中，氡的運移速度先逐步降低，后稍有抬升。在多層覆土中，每層覆土交界面處和灘面處的運移速度均小于均勻一層覆土中的速度。

由此可見，對于氡入口通量大的覆土中，氡的運移速度大。且在不同的孔隙度的土壤中，運移速度不同，孔隙度越大，速度越大。與均勻覆土覆蓋治理方式相比分層覆土邊界層處速度減小，采用分層覆土能夠有效降低氡的運移速度。

3.2 分層覆土下氡濃度變化

從模擬結果可知，如圖4所示，在分層覆土中入口處氡通量越大，氡在土壤運移過程中的質量分數和摩爾濃度越大。隨著運移的過程的發生，氡在覆土中的濃度隨著高度增加不斷減小。為進一步分析覆土表層氡濃度分布，選取入口中心點（0.5，0.5，0）和出口中心點（0.5，0.5，4）連線處，選取10個參考點考察氡濃度分布，由圖5可知，當氡入口通量為8Bqm2，16Bq/m2，24Bq/m2時，覆土表層面中心點處氡摩爾濃度分別為0.00045mol/m3， 0.00052mol/m3， 0.00068mol/m3，氡入口處的通量越大，覆土表面氡析出量越大。在0～2m的覆土高度中，氡濃度降低的速度較快，在2～4m的覆土中減低的速度較為平緩。由此可見，孔隙度小，擴散系數小的覆土能夠有效減少氡的析出。

圖6為氡入口通量為24Bq/m3情況下，不分層覆土與分層覆土中氡濃度比較圖。與均勻一層覆蓋治理方式相比，多層覆土表面中心點（0.5，0.5，4）處氡濃度降低了0.02307Bq/m3。尤其在第一層第二層孔隙度低的覆土中，抑制氡濃度的效果更加明顯。

因而，鈾尾礦庫治理采用多層覆土治理方式，并重點關注入口處的氡通量，減小入口處氡通量。對于入口處氡通量較大大的鈾尾礦庫，加蓋孔隙度低和擴散系數小的第一層和第二層覆土，可有效的減少氡的析出。

氡在分層覆土中運移過程中，氡濃度隨著時間變化而變化，圖7 為入口通量為24Bq/m3時，運移時間為15min，30min和60min時的氡濃度分布切面圖，濃度分別為0.00492 Bq /m3，0.011Bq/m3，0.0128 Bq /m3。隨著時間的推移，覆土中氡濃度含量增大。由圖8可知，根據摩爾濃度的分布，將覆土范圍分為高濃度（紅色區域）、中濃度（綠色區域）區域和低濃度（藍色區域）三個區域范圍，高濃度區域隨著時間推移逐漸增大。

基于CFD數值模擬方法，利用tecplot后處理軟件，采用可視化的方式展示氡運移規律。重點關注在不同時刻下上層覆土中氡的運移過程和濃度變化規律，并進一步通過速度矢量和濃度分布展示可視化模擬結果（如圖8所示），為運行中或退役后的鈾尾礦庫的治理提供必要的技術指導和理論支持。

4 結論與展望

本文參照低滲透的鈾尾礦庫的分層覆土設計理念，建立分層覆土模型，依據多孔介質模型修正后的基本控制方程，并基于CFD理論建立湍流模型和多組分輸運模型對鈾尾礦庫氡在多層覆土中的運移規律進行模擬，得到以下結論：

（1）針對鈾尾礦庫治理采用多層覆土治理方式，劃分臨時覆蓋層、排水層、隔離沙土層、儲存層、蓄水植被覆蓋層，相比均勻單層覆土治理形式能夠有效減少降低氡的運移速度，減低氡的灘面的析出濃度。（2）鈾尾礦庫治理采用多層覆土治理方式，并重點關注入口處的氡通量，減小入口處氡通量。對于入口處氡通量較大大的鈾尾礦庫，加蓋孔隙度低和擴散系數小的第一層和第二層覆土，可有效的減少氡的析出。（3）基于CFD數值模擬方法，采用可視化的方式展示氡運移規律。根據摩爾濃度的分布，將覆土范圍分為高濃度、中濃度和低濃度三個區域范圍，高濃度區域隨著時間推移逐漸增大。重點關注在不同時刻下上層覆土中氡的運移過程和濃度變化規律，并進一步通過速度矢量和濃度分布展示，為運行中或退役后的鈾尾礦庫的治理提供必要的技術指導和理論支持。

本論文考慮了分層覆土中孔隙度、擴散系數以及氡衰變系數對氡運移規律的影響，還存在諸多不足，下一步將對影響因素進行深度研究，深入探究土壤含水率、氣象條件等諸多綜合因素影響下的氡運移規律，為鈾尾礦庫治理提供具體客觀的技術支持和理論指導。

參考文獻

[1]IAEA，2006.Joint Convention on the Safety of Spent Fuel Management and on the Safety of Radioactive Waste Management.IAEA International Law Series No.1.IAEA.International Atomic Energy Agency（IAEA），Vienna.

[2]范楠彬，葉勇軍，李志，等.鈾尾礦充填體氡析出規律理論研究[J].安全與環境學報，2015，15（2）：236-239.

[3]王立恒，葉勇軍，丁德馨，等.破碎射氣介質內氡滲流運移和析出規律的實驗研究[J].原子能科學技術，2014，48（9）：1718-1724.

[4]戴劍勇， 石競羽.基于氡析出機制的鈾礦井巷道通風可靠性分析[J].核技術，2015，38（1）：34-38.

[5]葉勇軍，蘇杭，丁德馨，雷波，馮勝洋，梁濤，張運鋒.地下鈾礦山留礦法采場受限空間內氡遷移的數值模擬[J].安全與環境學報，2018，18（01）：217-223.

[6]H.A. Abdel Ghany， Aassy I E E ， Ibrahim E M ， et al. White Sand potentially Suppresses Radon Emission From Uranium Tailings[J]. Radiation Physics & Chemistry， 2018， 144：100-105.

[7]Kozak J A，Reeves H W， Lewis B A. Modeling radium and radon transport through soil and vegetation[J]. Journal of contaminant hydrology，2003，66（3）：179-200.

[8]Yakovleva V S .Solution of diffusion-advection equation of radon transport in many-layered geological media[J]. Nukleonika -Original Edition-， 2010， 55（4）：601-606.

[9]劉慶成，程業勛，章曄，等.介質中氡運移的模擬[J].東華理工大學學報（自然科學版），1995（4）：366-370.

[10]伍劍波，張慧，蘇鶴軍.斷層氣氡在不同類型覆蓋層中遷移規律的數值模擬[J].地震學報，2014，36（01）：118-128+159.

[11]何鴻展，宋翀芳，潘武軒，雷勇剛.基于CFD的防風抑塵網非均勻孔隙率的優化研究[J].中國環境科學，2016，36（06）：1697-1704.

[12]張輝，徐樂昌，李先杰.鈾礦山廢石堆與尾礦庫覆蓋抑氡實驗效果對比研究[J].輻射防護，2017，37（05）：387-392.

[13]葉巖，王岳.大氣土壤耦合的埋地燃氣管道泄漏擴散數值分析[J].中國安全生產科學技術，2018，v.14（4）：109-115.

[14]樂仁昌，賈文懿，吳允平.理想條件下氡及其子體運移新理論及其運移方程[J].物理學報，2003，52（10）：2457-2461.

[15]黃帥，彭小勇，張欣，郭蘭霞，柴旭昇.鈾尾礦庫灘面析出氡在覆土中運移的數值模擬[J].安全與環境學報，2014，14（03）：176-179.

[16]谷云慶，牟介剛，代東順，等.基于氣體射流的氣液兩相流動減阻特性[J].推進技術，2015，36（11）：1640-1647.

[17]馬盼軍，王哲，易發成，牛韶軍，周海田.某鈾尾礦庫周邊土壤中鈾元素的空間分布與污染評價[J].原子能科學技術，2017，51（05）：956-960.

[18]李志強，肖德濤，趙桂芝，等.鈾尾礦庫氡析出率多點快速測量方法研究與應用[J].四川大學學報（自然科學版），2018，55（06）：144-148.

收稿日期：2020-01-13

作者簡介：吳慧（1993-），女，南華大學碩士研究生，研究方向為核安全。

通訊作者：鄒樹梁（1956-），男，南華大學核設施應急安全作業技術與裝備湖南省重點實驗室教授，博士生導師。