葉勇軍，趙婭利，丁德馨，王立恒，曹雅琴，范楠彬

（1.南華大學(xué)環(huán)境保護與安全工程學(xué)院，湖南衡陽 421001；

2.南華大學(xué)鈾礦冶生物技術(shù)國防重點學(xué)科實驗室，湖南衡陽 421001）

堆浸鈾礦堆氡析出規(guī)律的實驗研究

葉勇軍1，趙婭利1，丁德馨2，＊，王立恒1，曹雅琴1，范楠彬1

（1.南華大學(xué)環(huán)境保護與安全工程學(xué)院，湖南衡陽 421001；

2.南華大學(xué)鈾礦冶生物技術(shù)國防重點學(xué)科實驗室，湖南衡陽 421001）

堆浸鈾礦堆是鈾礦山大氣氡的主要來源之一。為了揭示堆浸過程中礦堆表面的氡析出規(guī)律，以我國南方某鈾礦山新上堆的礦石為實驗對象，利用自主設(shè)計并制造的改進型一維堆浸實驗裝置，采用局部靜態(tài)法測定了不同噴淋狀態(tài)下和排水過程中礦堆表面的氡析出率。結(jié)果表明：在布液和改變噴淋強度的瞬間，礦堆表面的氡析出率突然增大，之后隨著噴淋的進行而逐漸減小并趨于穩(wěn)定，且噴淋強度越大，穩(wěn)定期間的氡析出率反而越?。辉谂潘^程中，氡析出率先迅速減小，而后逐漸增大趨于穩(wěn)定。

堆浸鈾礦堆；局部靜態(tài)法；氡析出率；噴淋強度

核能是人類使用的重要能源，已成為電力工業(yè)的重要組成部分。核電因不向空氣中排放碳類污染物，被稱為清潔能源。在人們越來越重視地球溫室效應(yīng)、氣候變化的形勢下，積極推進核電建設(shè)，對保障能源供應(yīng)與安全及保護環(huán)境等都具有重要的意義［1］。當前，我國正積極將核電納入優(yōu)先發(fā)展的能源產(chǎn)業(yè)。鈾作為核電原料，其需求量將隨核電產(chǎn)業(yè)的發(fā)展而增加。目前，鈾礦地表堆浸是我國天然鈾生產(chǎn)的主體技術(shù)之一，在堆浸過程中，礦石顆粒內(nèi)的鐳衰變形成的自由氡一部分溶解進入液相，另一部分溶入氣相，通過滲流和擴散遷移到礦堆表面進而釋放到大氣中，因此，堆浸鈾礦堆也是礦區(qū)大氣氡的主要來源之一［2］。氡的半衰期為3.8d，衰變形成218Po、214Pb和214Bi等短壽命子體，氡及其子體均具有放射性危害，研究表明，長期吸入高濃度的氡及其子體容易誘發(fā)肺癌［3－5］。為加強地表堆浸場區(qū)的氡防護水平，開展堆浸鈾礦堆的氡運移和析出規(guī)律是十分必要的。

氡析出是一個復(fù)雜的物理過程，是介質(zhì)內(nèi)部氡運移的邊界現(xiàn)象。對于堆浸體來說，其內(nèi)部的氡運移同時要受到其內(nèi)部的液體和氣體運移的影響［6］，表面的氡析出還受到當?shù)貧庀髤?shù)的影響，因此，堆浸鈾礦堆內(nèi)的氡遷移與析出過程更為復(fù)雜和特殊。目前，國內(nèi)外學(xué)者對氡析出規(guī)律的研究主要是基于擴散理論和滲流理論［78］，多是針對土壤、建筑材料以及鈾尾礦等常見的多孔介質(zhì)內(nèi)的氣相氡遷移［8－9］，對孔隙介質(zhì)內(nèi)液相的影響多考慮靜態(tài)水，集中在含水率對氡析出的影響規(guī)律的研究方面［10－14］。目前，針對像堆浸鈾礦堆等存在氣液兩相流的破碎射氣介質(zhì)表面氡析出規(guī)律的研究不多，且在氣液兩相流動作用下的氡運移機理逐漸成為國外學(xué)者新的研究熱點［15－17］。因此，深入開展堆浸鈾礦堆的氡析出規(guī)律研究，不僅有助于了解堆浸過程中堆浸鈾礦堆氡析出的機理，為今后合理確定堆浸場周圍的輻射照射劑量和防護措施提供依據(jù)，而且對進一步豐富和發(fā)展氣液兩相流動作用下的氡運移理論和實驗基礎(chǔ)具有重要的意義。為此，本文利用自主設(shè)計并研制的一維堆浸實驗裝置，通過室內(nèi)實驗研究特定實驗條件下堆浸鈾礦堆的氡析出規(guī)律。

1 實驗

1.1 實驗原理

1）堆浸鈾礦堆內(nèi)氡的運移理論

標準狀態(tài)下，222Rn是無色無味的透明氣體，密度為9.727kg／m3；能溶解于水等液體中，其在水中的溶解度系數(shù)隨著水溫的升高而降低，當水溫為20℃時，氡在水中的溶解度系數(shù)為0.252。在堆浸鈾礦堆內(nèi)，由于存在氣液兩相流，依據(jù)流動狀態(tài)和氣液兩相氡的濃度分布，將存在溶解平衡和非平衡兩種狀態(tài)。在一定溫度下，若氣液兩相中的氡濃度未達到溶解平衡，則在兩相之間存在氡的非等量遷移現(xiàn)象，此時，堆浸鈾礦堆內(nèi)一維氡非穩(wěn)定運移的控制方程［15］為：

式中：vg為礦堆中氣相的滲流速度，m／s；vw為礦堆中液相的滲流速度，m／s；cg為氣相中氡的濃度，Bq／m3；cw為液相中氡的濃度，Bq／m3；Dg為氣相中的氡在礦堆中的擴散系數(shù)，m2／s；Dw為液相中的氡在礦堆中的擴散系數(shù)，m2／s；θg為氣體體積分數(shù)；θw為水分體積分數(shù)；λ為氡的衰變常量，s－1；Ag為氡析出到礦堆氣相中的能力，Bq·m－3·s－1；Aw為氡析出到礦堆液相中的能力，Bq·m－3·s－1；Yg為氡從氣相到液相的遷移能力，Bq·m－3·s－1；Yw為氡從液相到氣相的遷移能力，Bq·m－3·s－1；t為時間，s；Q為單位體積礦堆的222Rn產(chǎn)生率，Bq·m－3·s－1。

若氣液兩相中的氡濃度達到溶解平衡，則堆浸鈾礦堆內(nèi)一維氡穩(wěn)定運移的控制方程［18］為：

式中：Deff為氣液混合相的有效擴散系數(shù)［19－20］，m2／s；x為距表面的距離（取向下的方向為正方向），m；K為水／氣分配系數(shù)；v為達西滲透速率，m／s；m為礦堆內(nèi)的液體飽和度；n為礦堆孔隙度。

由上述公式可知，當外界環(huán)境條件相對穩(wěn)定時，針對特定的堆浸鈾礦堆，其表面的氡析出率受堆浸鈾礦堆內(nèi)氣液兩相的流動情況和氡氣在兩相中濃度分布的影響。

2）局部靜態(tài)法測氡析出率

局部靜態(tài)法又稱積累法［21］，長期以來被廣泛用于測量各種介質(zhì)表面上的氡析出率。本實驗采用局部靜態(tài)法測定特定實驗條件下礦石表面的氡析出率，具體方法為：（1）將氡集氣罩罩在待測礦石表面；（2）用夾有橡膠密封墊圈的法蘭將氡集氣罩和裝礦筒體進行密封，防止氡氣泄漏；（3）10min后，通過氡集氣罩中部的取樣孔取樣，然后將所取樣品注入準備好的閃爍室中。按下式計算礦樣表面氡的析出率J：

式中：Δc為氡集氣罩內(nèi)氡濃度c與大氣中氡濃度ca之差，Bq／m3；t為集氣時間，s；H為氡集氣罩的高度，m；S為被測礦樣的表面積，m2。

本實驗采用閃爍室法測量氡濃度，相關(guān)設(shè)備為：0.5L的閃爍室、50mL的玻璃注射器、FD-125型氡釷分析儀、FH463B型智能定標器及真空泵。由于實驗時室內(nèi)通風條件好，大氣中的氡濃度約為27～35Bq／m3，遠小于所測氡濃度，故此處將大氣中的氡濃度ca取為零。

氡濃度c按下式計算：

式中：Vq為取樣體積，mL（本實驗為50mL）；Vs為閃爍室的體積，mL；nc為本底加樣品計數(shù)率的平均值，min－1；nb為本底計數(shù)率的平均值，min－1；Ks為刻度系數(shù)，Bq·m－3·min。

3）噴淋強度的計算

按下式計算噴淋強度：

式中：I為噴淋強度，mL·cm－2·min－1；q為液體流量，mL·min－1；A為實驗筒體的橫截面積，cm2。

1.2 實驗材料與裝置

選用我國南方某鈾礦山新上堆浸場的鈾礦石，將礦石樣按0～0.5、0.5～1、1～2、2～3、3～4、4～5、5～6、6～7、7～8、8～9mm的粒徑進行篩分，并將其按最大粒徑尺寸為9mm、粒徑分布分形維數(shù)D＝1.8進行混合，試樣的粒徑級配列于表1，其孔隙率為42.2%。

表1 D＝1.8時試樣的粒徑級配Table 1 Size grading of samples under D＝1.8

自主設(shè)計并制造的改進型一維堆浸實驗裝置如圖1所示。該裝置由實驗柱和布液裝置組成，實驗柱是由底部集水容器，上、下部裝礦筒體和頂部氡集氣罩等4節(jié)有機玻璃筒體組成，4節(jié)筒體間用夾有橡膠密封墊圈的法蘭串聯(lián)而成，下部裝礦筒體與底部集水容器之間的法蘭夾有厚度為1cm的孔板。每節(jié)裝礦筒體外直徑為22cm，內(nèi)直徑為20cm。裝置總高度約為150cm，其中中間兩節(jié)裝礦筒體的高度為100cm。頂部氡集氣罩的高度為15cm，其上部中心處為氡氣取樣孔，氣壓平衡孔位于距中心10cm處。布液裝置主要由帶溢流孔的高位水槽、L型溢流管、流量計和布液滴管構(gòu)成。為了能讓礦樣表面的布液均勻，在礦樣表面中心放置一扁平鵝卵石；噴淋時，通過氡集氣罩中心的氡氣取樣孔，將布液裝置的布液滴管置入礦體柱內(nèi)部，此時仍通過該取樣孔取氡樣；為防止礦樣經(jīng)過孔板滲漏，在孔板上放有孔徑為0.5mm的砂網(wǎng)。

圖1 改進型一維堆浸實驗裝置示意圖Fig.1 Sketch of improved one-dimensional heap-leaching experimental device

1.3 實驗過程

1）準備：按照粒徑分布分形維數(shù)的要求，配好足夠質(zhì)量的礦石樣，并測閃爍室本底。

2）裝樣：取配好的試樣稱重；將試樣分層裝入下部裝礦筒體，用木制擊錘輕輕搗實。當試樣面距下部裝礦筒體頂邊緣50mm左右時，用夾有橡膠密封墊圈的法蘭將裝礦筒體的上下部分相連，并繼續(xù)采用同樣的方法裝填試樣，直至試樣面與上部裝礦筒體上邊緣齊平；接著將準備好的鵝卵石平放在水平礦樣中心位置。將礦體柱置于通風環(huán)境良好的房間，靜置一段時間。

3）取氣樣：每次測定累積氡濃度前，用夾有橡膠密封墊圈的法蘭將集氣罩與裝礦筒體連接，累積10min后，用50mL的玻璃注射器通過集氣罩中心的氡氣取樣孔進行取氣樣，并將取得的氣樣迅速注入閃爍室內(nèi)；之后，將集氣罩取下并驅(qū)趕內(nèi)部積聚的氡氣，為下個測點作準備。

4）集氣罩內(nèi)累積氡濃度的測量：調(diào)節(jié)并控制流量計的讀數(shù)，采用過程3的方法取樣，測量未噴淋、噴淋過程中及排水期間集氣罩內(nèi)的累積氡濃度。噴淋強度分別為0.031 8、0.063 7、0.127 4、0.191 1mL·cm－2·min－1（由式（8）可知對應(yīng)的流量分別為10、20、40、60mL·min－1）。

2 氡濃度及析出率的測量結(jié)果

未噴淋階段集氣罩內(nèi)累積的氡濃度和礦樣表面的氡析出率測量結(jié)果列于表2，噴淋過程中、噴淋階段及排水階段礦樣表面的氡析出率示于圖2～4。

表2 未噴淋階段氡濃度和氡析出率的測量結(jié)果Table 2 Measured results of radon concentration and radon exhalation rate before spraying

圖2 噴淋過程中氡析出率的變化Fig.2 Variation of radon exhalation rate during spraying

3 破碎射氣介質(zhì)氡運移和析出規(guī)律

由表2可知，當外界環(huán)境條件相對穩(wěn)定時，放置足夠長時間的鈾礦石，其析出率基本保持不變，平均為4.05Bq·m－2·s－1。由圖2可知，噴淋強度一定時，礦石堆表面的氡析出率在噴淋初始突然增加，達到最大值后又逐漸減小，減小到一定值后趨于穩(wěn)定。結(jié)合堆浸鈾礦堆內(nèi)氡的運移理論，可解釋為，剛開始噴淋時，水首先進入礦堆上部礦石的孔隙內(nèi)，由于水的密度大，其將孔隙中的部分氣體擠壓出來，被擠壓出的氣體攜帶氡氣向上運移，此時，由于水沿礦石向下滲流的速度較快且量有限，因此只有少部分氡氣溶解到水中，而大部分氡氣向上運移到集氣罩的集氣空間中。隨著噴淋的持續(xù)進行，在重力等因素的影響下，礦柱內(nèi)濕礦石的量增多且含水量也較大，孔隙中已有的氡氣和固體礦石顆粒新產(chǎn)生的氡氣有足夠的時間溶解于水，這樣使越來越多的氡氣隨著水向下運移，從而向上運移到集氣空間中的氡氣就逐漸減少，當柱內(nèi)液相的流動穩(wěn)定后氣相氡濃度分布也重新獲得穩(wěn)定。

圖3 噴淋階段的氡析出率Fig.3 Radon exhalation rate in stage of spraying

圖4 排水階段氡析出率的變化Fig.4 Variation of radon exhalation rate in stage of drainage

由圖3可看出，同一噴淋強度下，滲流非穩(wěn)定期間的氡析出率高于穩(wěn)定期間的氡析出率。這是由于在噴淋未穩(wěn)定狀態(tài)，實驗柱內(nèi)的水量增加的同時，氡氣溶解于水的速度不及水侵入孔隙擠壓出氡氣的速度，致使較多的氡氣向上擴散至集氣空間。噴淋強度越大，非穩(wěn)定和穩(wěn)定時期氡析出率越小。表明隨著噴淋強度的增大，溶解在水中的氡量增加，流動水帶走的氡量也越多，所以礦石表面的氡析出率隨噴淋強度的增大而減小。

由圖4可看出，在排水過程中，氡的析出率先降低然后逐漸升高，最終趨于穩(wěn)定，穩(wěn)定時的氡析出率小于未噴淋階段礦石表面的氡析出率。這是由于排水初期，在重力作用下，礦柱失水速度較大，使礦石表面的空氣流向礦石內(nèi)，礦石表層空氣滲流的方向與氡擴散運移方向相反，滲流對氡析出起抑制作用，導(dǎo)致氡析出率降低；隨著排水時間的延長，排水速率越來越小，故滲流速度越來越小，從而氡的析出率逐漸增大；當滲流速度足夠小時，氡析出率趨于穩(wěn)定。由于排水穩(wěn)定期間礦柱內(nèi)氡的遷移主要靠擴散，依據(jù)式（5），濕度越大的礦石擴散系數(shù)越小，因此排水穩(wěn)定時濕礦石的氡析出率小于未噴淋階段礦石表面的氡析出率。

4 結(jié)論

1）噴淋強度越大，氡析出率越小，并且在同一噴淋強度下滲流穩(wěn)定時的氡析出率小于不穩(wěn)定時的氡析出率；在排水過程中，氡析出率先降低，之后逐漸升高，到一定值后保持穩(wěn)定，且穩(wěn)定時濕礦石的氡析出率小于未噴淋階段礦石表面的氡析出率。

2）結(jié)合堆浸鈾礦堆內(nèi)氡的運移理論分析可知，本實驗條件下噴淋過程中影響氡析出的主要因素為礦體柱內(nèi)水進入礦石孔隙的速度和氡氣溶解于水的總量。水進入礦石孔隙的速度越快，在一定時間內(nèi)從孔隙中析出的氡氣越多，溶解于水的總量越少，礦石表面的氡析出率越大。

3）在開始布液噴淋或增大噴淋強度時，鈾礦堆表面的氡析出率會有一急劇升高的過程，因此，在實際堆浸噴淋操作時，開始噴淋或增加噴淋強度時工作人員應(yīng)盡可能遠離礦堆，或做好必要的輻射防護措施。

［1］ 中華人民共和國國家發(fā)展和改革委員會.核電中長期發(fā)展規(guī)劃（2005—2020年）［EB／OL］.［2007-11-02］.http：∥www.ndrc.gov.cn／gzdt／200711／t20071102＿170163.html.

［2］ 畢忠偉，丁德馨，段仲沅.鈾礦開采對環(huán)境的影響及治理的特殊性［J］.安全與環(huán)境工程，2004，11（4）：40-42.

BI Zhongwei，DING Dexin，DUAN Zhongyuan.The specificity of impact and harnessing induced by uranium mining［J］.Safety and Environmental Engineering，2004，11（4）：40-42（in Chinese）.

［3］ KREWSKI D，LUBIN J，ZIELINSKI J，et al.Residential radon and risk of lung cancer：A combined analysis of 7North American case-control studies［J］.Epidemiology，2005，16（2）：137-145.

［4］ DARBY S，HILL D，AUVINEN A，et al.Radon in homes and risk of lung cancer：Collaborative analysis of individual data from 13European case-control studies［J］.British Medical Journal，2005，330：223-226.

［5］ TOMá?EK L，DARBY S C，F(xiàn)EARN T.Patterns of lung cancer mortality among uranium miners in west bohemia with varying rates of exposure to radon and its progeny［J］.Radiation Research，1994，137（2）：251-261.

［6］ GIBBONS D A，KALIN R M，JOHNSTON H.Modelling of radon-222as a proxy for vapour phase transport in the unsaturated zone［J］.Calibration and Reliability in Groundwater Modeling，2000，265：2 000.

［7］ SAHOO B K，MAYYA Y S，SAPRA B K.Radon exhalation studies in an Indian uranium tailings pile［J］.Radiation Measurements，2010，45：237-241.

［8］ SAHOO B K，SAPRA B K，GAWARE J J.A model to predict radon exhalation from walls to indoor air based on the exhalation from building material samples［J］.Science of the Total Environment，2011，409：2 635-2 641.

［9］ SAHOO B K，NATHWANI D，EAPPEN K P，et al.Estimation of radon emanation factor in Indian building materials［J］.Radiation Measurements，2007，42：1 422-1 425.

［10］葉勇軍，丁德馨，羅潤，等.體積含水率對鈾尾礦氡析出率影響的試驗研究［J］.安全與環(huán)境學(xué)報，2012，12（3）：124-126.

YE Yongjun，DING Dexin，LUO Run，et al.On the effect of volumetric water content on theradon exhalation rate of uranium tailings［J］.Journal of Safety and Environment，2012，12（3）：124-126（in Chinese）.

［11］劉小松，劉純魁.鈾礦石濕度對氡析出率影響的研究［J］.鈾礦冶，2004，23（3）：162-164.LIU Xiaosong，LIU Chunkui.Effect of uranium ore humidity on radon emanation rate［J］.Uranium Mining and Metallurgy，2004，23（3）：162-164（in Chinese）.

［12］SUN H B，F(xiàn)URBISH D J.Moisture content effect on radon emanation in porous media［J］.Journal of Contaminant Hydrology，1995，18：239-255.

［13］STRANDEN E，KOLSTAD A K，LIND B.The influence of moisture and temperature on radon exhalation［J］.Radiation Protection Dosimetry，1984，7（1-4）：55-58.

［14］SAKODA A，ISHIMORI Y，HANAMOTO K，et al.Experimental and modeling studies of grain size and moisture content effects on radon emanation［J］.Radiation Measurements，2010，45：204-210.

［15］李廣悅.松散破碎射氣介質(zhì)中氡運移的氣液兩相耦合模型及其工程應(yīng)用［D］.長沙：中南大學(xué)，2011.

［16］CALUGARU D G，CROLET J M.Identification of radon transfer velocity coefficient between liquid and gaseous phases［J］.Comptes Rendus Mecanique，2002，330：377-382.

［17］LOISY C，CEREPI A.Radon-222as a tracer of water-air dynamics in the unsaturated zone of a geological carbonate formation：Example of an underground quarry（Oligocene Aquitain Limestone，F(xiàn)rance）［J］.Chemical Geology，2012，296-297：39-49.

［18］HO C K.Analytical risk-based model of gaseous and liquid-phase radon transport in landfills with radium sources［J］.Environmental Modelling ＆Software，2008，23（9）：1 163-1 170.

［19］ROGERS V C，NIELSON K K，KALKWARF D R.Radon attenuation handbook for uranium mill tailing cover design，NUREG／CR-3533［R］.US：Nuclear Regulatory Commission，1984.

［20］ROGERS V C，NIELSON K K.Correlations for predicting222Rn diffusion coefficients of soils［J］.Health Physics，1991，61（2）：225-230.

［21］張哲.地下工程與人居環(huán)境氡防護技術(shù)［M］.北京：原子能出版社，2010.

Experimental Study on Radon Exhalation Rule of Heap-leaching Uranium Ore Heap

YE Yong－jun1，ZHAO Ya-li1，DING De-xin2，＊，WANG Li-heng1，CAO Ya-qin1，F(xiàn)AN Nan-bin1
（1.School of Environmental Protection and Safety Engineering，University of South China，Hengyang421001，China；2.Key Discipline Laboratory for National Defense for Biotechnology in Uranium Mining and Hydrometallurgy，University of South China，Hengyang421001，China）

The heap-leaching uranium ore heap is one of the main sources of atmospheric radon of uranium mine.In order to reveal the radon exhalation rule on the surface of the heap-leaching uranium ore heap，taking the new ores from certain uranium mine in southern China as the experimental samples，the improved one-dimensional heap-leaching experimental device which was designed and made autonomously was used to measure the radon exhalation rate on the surface of the ore heap under different spray conditions and during the drainage with the local static method.The results show that in the moment of beginning spray and changing spray intensity，the radon exhalation rate of the ore heap surface increases suddenly，then gradually decreases and finally reaches stable along with the spray continuing.It is also found that the radon exhalation rate ofthe ore heap surface during the stable period of the seepage in the ore heap decreases with the increase of spray intensity.In the process of drainage，the radon exhalation rate decreases rapidly，then gradually increase and stabilizes finally.

heap-leaching uranium ore heap；local static method；radon exhalation rate；spray intensity

TD72；TL75.1

碼：A

：1000-6931（2015）01-0187-06

10.7538／yzk.2015.49.01.0187

2013-10-23；

2014-01-18

國家自然科學(xué)基金資助項目（11105069）

葉勇軍（1979—），男，湖北浠水人，副教授，博士研究生，從事鈾礦開采安全理論與技術(shù)研究

＊通信作者：丁德馨，E-mail：dingdxzzz＠163.com