通風管道長度變化對鈾礦井下氡濃度影響的模擬研究

蘇曉書，冀 東，劉 穎，李元崗

(中核第四研究設計工程有限公司，河北 石家莊 050021)

鈾礦開采屬于天然放射性物質開采的開放型作業，放射性物質向環境釋放，會對環境、公眾產生影響。因此，必須在鈾礦開采過程中采取有效的防護措施，對井下主要危害因素：長壽命氣溶膠、γ射線、氡及其子體進行防治[1]。其中，氡及氡子體的內照射對井下工作人員輻射危害最大。據統計，礦工職業照射劑量中由氡及氡子體引起的照射劑量占89.8%，氡及其子體的內照射是防治重點[2-5]。因此，掌握氡在鈾礦山井下介質中傳輸規律以及影響氡濃度遷移分布的因素是提高防護水平，控制井下氡濃度不超過國家標準限值的重要前提[6]。

目前，國內外未見針對井下通風管道長度變化對氡濃度影響的研究，但英國、荷蘭、印度等國對土壤、混凝土、大理石等材料的氡析出率有較深入研究，且利用流體動力學計算模塊(CFD)模擬了地下室氡遷移全過程以及室內222Rn、220Rn及其子體的濃度與三維分布。另外，我國研究者也研究了溫濕度、孔隙度、含水量、風壓等氡析出率影響因素，井下防氡措施的種類及效果，井下通風方式、通風量對氡濃度的影響，這對本研究一些定性分析有指導意義[7-12]。

本研究在完成軟件模擬井下風場現狀，篩選空氣流動性較差位置等前期工作的基礎上，有針對性的改變相應位置處通風管道的長度，分析管道長度變化對氡濃度的影響，判斷此措施能否起到降低井下氡濃度，為后續工程實施提供理論參考。

1 軟件簡介

ANSYS是融結構、流體、電場、磁場、聲場分析于一體的大型通用有限元分析軟件，內嵌上述各模塊所需的數學模型。本研究使用的流體模塊，自帶質量、動量、能量守恒方程，以及標準k-ε模型、RNGk-ε模型等湍流數學模型。該軟件由ANSYS公司開發，能與多數計算機輔助設計(CAD)軟件接口，實現數據的共享和交換，是現代產品設計中最為廣泛使用的計算機輔助工程(CAE)工具之一，ANSYS Workbench15.0是ANSYS公司出品的新一代仿真平臺，可滿足多種工程的仿真需要[13]。

本研究利用ANSYS Workbench15.0中流體動力學計算模塊(CFD)對礦井進行建模，然后經網格劃分、參數設置、模型運算后，觀察改變井下通風管長度后對氡運移的影響。

2 鈾礦井下模型建立

2.1 模型建立

本研究以江西某礦井的一個中段為模型進行模擬分析，以開采圖紙為基準，按照1∶1的比例在坐標系中創建模型，繪圖時保證所有巷道均連通。該中段有1個中段車場，1條主巷道，9條主穿，2個采場，詳細參數列于表1，平面布置及關注點位情況示于圖1。

理論及前期國內外研究成果表明[7-12]，風速、風壓、溫度、濕度均會對井下氡析出率產生影響，預計本研究中通風管道的變化會對井下整個風場產生影響，進而影響氡濃度。因此，在工程實施前利用ANSYS Workbench15.0軟件驗證管道變化對氡濃度影響的效果。該中段共有9根通風管，本研究建立的模型調整了其中5根通風管的長度，詳細調整方式列于表2。

表1 某礦井中段車場、巷道、主穿參數Table 1 Middle section of a mine yard, roadway, the main wear parameters

圖1 鈾礦井某中段平面布置及模擬點位圖Fig.1 A mid-section uranium mine layout and simulation of bitmap

名稱長度變化前/mm長度變化后/mm通風管直徑/mm主穿1通風管53 13085 000300主穿2通風管46 37084 000300主穿3通風管62 22096 000300主穿5通風管161 000155 000300主穿7通風管47 78044 000300

2.2 網格劃分

由于該中段尺寸較大，故在模型建立過程中，對巷道內流體流動影響較小的巷道壁面的凸起、凹陷、折邊等進行光滑化處理。

進行巷道內氡濃度擴散模擬計算時，湍流模型采用k-ε模型。對于本模型，采用非結構化網格的計算結果與采用結構化網格的計算結果相差不大，在可接受范圍內，考慮到計算量及計算資源，采用較為簡單的非結構化網格進行計算。由于通風管道與巷道尺寸有一定差距，因此對通風管壁面區域進行加密。通風管道長度變化后的模型中一共有3 738 933個節點，17 541 147個網格元素，模型網格質量最小為0.162 7，最大為1，平均值為0.832 62，滿足網格質量最小不低于0.1，網格質量平均值高于0.80的要求，可進行數值模擬計算，通風管道長度變化后模型網格劃分情況示于圖2。

圖2 通風管道長度變化后模型網格劃分圖Fig.2 Model grid division diagram after ventilation duct length change

2.3 邊界條件

根據現場情況及相關實測數據，結合數學模型和Fluent數值模擬方法，對區域網格進行自適應等調試，數值模擬參數設定如下：湍流模型采用k-ε模型，系統壓力為105 400 Pa；入口邊界為velocity inlet類型(速度入口)，入口速度為1.8 m/s，代表風流以一定速度吹進巷道內；水力直徑(等效直徑)為2 m，湍流強度為3.32%；出口邊界為pressure outlet類型(壓力出口)，代表排風井出口為壓力排出，表壓為0 Pa；采場和巷道為mass flow inlet類型，代表氡以一定的析出速率進入計算模型(氡初始濃度為0)。由于采場壁面和巷道壁面的析出率有較大差異，因此依據圖1中31個點位處分別設置，不同溫濕度下，礦樣及巖樣氡析出率實測數據：采場1處，氡析出率取值0.58 Bq·m-2·s-1(溫度21.5 ℃)，采場2處，氡析出率取值0.55 Bq·m-2·s-1(溫度16.8 ℃)，對于圍巖處(以10點位為例)氡析出率取值0.005 Bq·m-2·s-1(溫度18 ℃)。另外，由于井下壁面較粗糙，采用面積展平法對壁面進行光滑處理，避免氡析出量數值偏小，但該尺度下壁面是否光滑對整體風場影響較小。

混合物材料設置為氡氣在空氣中擴散，物質輸送的問題在物質模型(species model)面板設置，選擇組分傳輸模型。在軟件的材料庫中，自定義氡氣的輸入參數(1 Bq氡相當于1.77×10-16g)。在材料面板選取自定義的氡氣作為計算材料。

3 模擬結果及分析

3.1 模擬結果

根據現場調研資料和監測數據，經過多次調試，最終將通風管入口的邊界條件設置為風扇入口，壓力為200 Pa，出口為自然流出邊界，通風管中的空氣經通風管出口后進入巷道內。

參數設置完全后，通過ANSYS Workbench15.0軟件CFD模塊計算，通風管長度變化后氡濃度變化情況列于表3。巷道內31個點位氡濃度在通風管道長度變化后也發生了變化。其中，氡濃度降低幅度較大的典型點位，如10點位、20點位的計算結果示于圖3。

表3 通風管長度變化后氡濃度變化Table 3 The changes in radon concentration after ventilation

3.2 結果分析

本次研究設置31個預測點位，基本均勻分布于主巷內，有采場的支巷加密預測點位，詳見圖1。由表4中統計結果可知，當通風管道長度發生變化時，有29個點位氡濃度降低，僅有2個點位氡濃度升高，整體降氡效果較好。

a——10點位通風管長度變化前氡質量分數；b——10點位通風管長度變化后氡質量分數；c——20點位通風管長度變化前氡質量分數;d——20點位通風管長度變化后氡質量分數圖3 典型模擬點位氡濃度變化情況a——10 point ventilation pipe length before the change in radon mass fraction; b——10 point ventilation pipe length change;c——20 point ventilation pipe length before the change in radon mass fraction; d——20 point ventilation pipe length changeFig.3 Typical simulated point radon concentration changes

降低幅度<00%～10%10%～20%20%～30%30%～40%40%～50%>50%點位個數2688421

結果表明，更改通風管長度可有效降低井下氡濃度。其理論基礎是井下通風措施優化導致風場發生變化，進而優化了井下通風效果，降低了氡濃度。本研究提供的方法僅是降氡措施的一種，其他如改變井下通風量，局部涂抹防氡涂料等都會對降氡有一定效果，具體降低數值需要建模計算，可在工程優化措施實施之前對其進行模擬論證，確保工程優化措施有效性。

4 結論

首次利用ANSYS Workbench15.0軟件對某鈾礦山一個中段通風管道長度變化前后的氡濃度情況進行了模擬分析，通過中段整體建模、通風管道模型、網格劃分、參數設置、運行計算等，得到了通風管道變化前后的氡濃度模擬結果。經分析，在其他條件不變的情況下，排風管長度變化可有效降低井下氡濃度，該方法可在工程優化措施實施前論證其有效性，可為后期實施提供有效依據。

參考文獻：

[1] 胡鵬華,李先杰. 我國鈾礦通風降氡現狀分析[J]. 輻射防護,2011,31(3):178-183.

Hu Penghua, Li Xianjie. Practice in inspection for low and intermediate level radioactive waste packages[J]. Radiation Protection, 2011, 31(3): 178-183(in Chinese).

[2] 鄧文輝,李先杰,黃銀萍. 鈾礦山井下氡及氡子體濃度管理限值的探討[J]. 輻射防護通訊,2004,24(3):27-30.

Deng Wenhui, Li Xianjie, Huang Yinping. Discussion on concentration control limit of radon and its progeny inunderground uranium mine[J]. Radiation Protection Bulletin,2004, 24(3): 27-30(in Chinese).

[3] Gherghel C, De Souza E. Ventilation requirements for uranium mines[C]∥12th U.S./North American Mine Ventilation Symposium. Reno, NV, USA: Wallace K 2005: 65-69.

[4] Mudd G M. Radon sources and impacts: a review of mining and non-mining issues[J]. Rev Environ Sci Biotechnol, 2008, (7): 325-353.

[5] 李先杰,鄧文輝，薛建新,等. 關于我國鈾礦冶個人劑量約束值的討論[J]. 輻射防護,2008,28(1):13-17.

Li Xianjie, Deng Wenhui, Xue Jianxin, et al. Discussion about individual dose constraint value at uranium mine and mill in china[J]. Radiation Protection, 2008, 28(1): 13-17(in Chinese).

[6] 戴劍勇,石競羽. 基于氡析出機制的鈾礦井巷道通風可靠性分析[J]. 核技術,2015,38(1):1-5.

Dai Jianyong, Shi Jingyu. Reliability analysis of tunnel ventilation based on the mechanism of radon exhalation in uranium mine[J]. Nuclear Techniques, 2015, 38(1): 1-5(in Chinese).

[7] Wang Fan. Radon entry, migration and reduction in houses with cellars[J]. Building and Environment, 2002, 37(11): 1 153-1 165.

[8] Zhuo W, Iida T, Moriizumi J, et al. Simulation of the concentrations and distributions of indoor radon and thoron[J]. Radiation Protection Dosimetry, 2001, 93(4): 357-368.

[9] Merkel B J, Hasche-Berger A. Uranium in the environment mining impact and consequences[M]. Berlin: Springer, 2006: 593-599.

[10] 李萬偉,李曉紅,劉興榮. 密閉環境中溫度和相對濕度對墻磚氡析出率的影響[J]. 環境與健康雜志,2007,24(7):519-521.

Li Wanwei, Li Xiaohong, Liu Xingrong. Effect of temperature and relative humidity on radon exhalation rate of glazed tile in hermetic environment[J].Journal of Environment and Health, 2007, 24(7): 519-521(in Chinese).

[11] 劉小松,劉純魁. 鈾礦石濕度對氡析出率影響的研究[J]. 鈾礦冶,2004,23(3):163-165.

Liu Xiaosong, Liu Chunkui. Effect of uranium ore humidity on radon emanation rate[J]. Uranium Mining and Metallurgy, 2004, 23(3): 163-165(in Chinese).

[12] 葉勇軍,丁德馨,羅潤,等.體積含水率對鈾尾礦氡析出率影響的試驗研究[J]. 安全與環境學報,2012,12(3):124-126.

Ye Yongjun, Ding Dexin, Luo Run, et al. On the effect of volumetric water content on the radon exhalation rate of uranium tailings[J].Journal of Safety and Environment, 2012, 12(3): 124-126(in Chinese).

[13] 劉浩. ANSYS15.0有限元分析從入門到精通[M]. 北京:機械工業出版社,2014.