不同噴淋強(qiáng)度對(duì)核素Sr-90在土壤中遷移的影響*

朱 君 鄧安嫦 石云峰 陳 超 劉團(tuán)團(tuán) 張艾明

（中國(guó)輻射防護(hù)研究院，太原 030006）

不同噴淋強(qiáng)度對(duì)核素Sr-90在土壤中遷移的影響*

朱 君 鄧安嫦 石云峰 陳 超 劉團(tuán)團(tuán) 張艾明?

（中國(guó)輻射防護(hù)研究院，太原 030006）

基于室內(nèi)土柱法，建立120 ml d-1、60 ml d-1、30 ml d-1三種不同噴淋強(qiáng)度下放射性核素Sr-90在砂土介質(zhì)中的動(dòng)態(tài)土柱遷移試驗(yàn)，試驗(yàn)時(shí)間分別為51 d、102 d、204 d，研究不同噴淋強(qiáng)度對(duì)Sr-90遷移規(guī)律的影響。結(jié)果表明，三種不同噴淋強(qiáng)度下Sr-90的濃度峰在土柱體中的遷移距離分別為42.3 cm、46.2 cm、44.4 cm，遷移距離基本一致，與噴淋強(qiáng)度和時(shí)間成正比。此外，Sr-90在土柱體剖面的濃度分布曲線(xiàn)均存在明顯的不對(duì)稱(chēng)性和拖尾，即Sr-90從源層上洗脫下來(lái)后，受到砂土介質(zhì)的吸附作用，待濃度峰值通過(guò)后，更多的Sr-90從之前的砂土介質(zhì)中解吸出來(lái)，使得很長(zhǎng)時(shí)間段內(nèi)土柱體中保持相對(duì)較高的濃度。應(yīng)用Hydrus-1d軟件建立平衡吸附、非平衡吸附兩種模式下的核素遷移數(shù)值模型，發(fā)現(xiàn)在非平衡吸附模式中考慮一級(jí)速率系數(shù)β能更好地描述濃度分布曲線(xiàn)的拖尾，通過(guò)與實(shí)測(cè)濃度值擬合得到Sr-90的分配系數(shù)為0.73 ml g-1、0.70 ml g-1、0.72 ml g-1，噴淋強(qiáng)度對(duì)分配系數(shù)影響不大。

噴淋強(qiáng)度；土柱試驗(yàn)；核素遷移；平衡吸附；非平衡吸附

核素遷移能力是放射性廢物處置場(chǎng)選址、設(shè)計(jì)、建造和安全評(píng)估的一個(gè)重要參考指標(biāo)，分配系數(shù)作為衡量核素遷移能力的一個(gè)關(guān)鍵參數(shù)，現(xiàn)階段其測(cè)量方法分為靜態(tài)批式法和動(dòng)態(tài)土柱法。靜態(tài)法忽略了固相介質(zhì)的孔隙度、水流速度、固液比等因素對(duì)分配系數(shù)的影響，尤其是“固體濃度效應(yīng)”［1-2］。動(dòng)態(tài)土柱法則較切合實(shí)際更具優(yōu)勢(shì)，是研究核素遷移行為及遷移參數(shù)的主要手段。但是由于延遲作用［3-4］，動(dòng)態(tài)土柱法往往需要較長(zhǎng)的試驗(yàn)周期才能得到結(jié)果，通常采用加速?lài)娏艿霓k法。放射性廢物處置場(chǎng)多位于西北干旱、半干旱地區(qū)，實(shí)際年降雨量遠(yuǎn)遠(yuǎn)小于加速?lài)娏艿乃浚铀賴(lài)娏艿脑囼?yàn)結(jié)果能不能反應(yīng)實(shí)際條件下的核素遷移規(guī)律，一直是爭(zhēng)論的焦點(diǎn)問(wèn)題。關(guān)于上述問(wèn)題的研究，國(guó)內(nèi)報(bào)道并不多見(jiàn)，且有關(guān)內(nèi)容主要集中在利用濃度穿透曲線(xiàn)的峰值遷移速率［5-6］，反映噴淋強(qiáng)度對(duì)核素遷移行為的影響。

放射性廢物處置場(chǎng)中的關(guān)鍵核素鍶-90（以下簡(jiǎn)稱(chēng)Sr-90），因?yàn)榘胨テ陂L(zhǎng)、毒性大，參與鈣代謝，對(duì)骨骼有親合力，引起白血病和骨癌而被關(guān)注。于靜［5］和王志明［6］等采用室內(nèi)土柱試驗(yàn)，通過(guò)底部流出液，分別得到了不同噴淋強(qiáng)度下Sr-90在沙土和黃土中的穿透曲線(xiàn)，發(fā)現(xiàn)噴淋強(qiáng)度與Sr-90濃度峰的遷移速率基本呈線(xiàn)性關(guān)系，但是沒(méi)有從遷移理論和模型更加深入地分析現(xiàn)象原因。Li等［7］通過(guò)土柱試驗(yàn)的穿透曲線(xiàn)，發(fā)現(xiàn)Sr-90在石英砂中的遷移是一個(gè)非平衡過(guò)程，應(yīng)用非平衡吸附模式［8-12］計(jì)算得到的結(jié)果較平衡吸附模式［13-15］能更好地與穿透曲線(xiàn)吻合。但利用濃度穿透曲線(xiàn)研究核素遷移模型和噴淋強(qiáng)度對(duì)核素遷移行為、參數(shù)的影響，只考慮了液相中的核素濃度，忽略了被滯留在土柱體固相介質(zhì)中的濃度，無(wú)法全面地反映核素在土壤剖面的遷移過(guò)程和參數(shù)。在探討噴淋強(qiáng)度對(duì)核素遷移規(guī)律的影響時(shí)，應(yīng)考慮滯留在土壤介質(zhì)中的那部分濃度，并結(jié)合遷移模型分析進(jìn)一步開(kāi)展研究。本文采用室內(nèi)土柱法，設(shè)計(jì)三種不同的噴淋強(qiáng)度，充分考慮固相滯留作用，根據(jù)Sr-90在土柱體垂直剖面的濃度分布曲線(xiàn)，并結(jié)合核素遷移數(shù)值模型，研究加速?lài)娏軐?duì)Sr-90在砂土介質(zhì)中遷移的影響，為干旱區(qū)放射性廢物處置場(chǎng)的安全評(píng)估提供參數(shù)支持。

1 材料與方法

1.1 供試材料

土壤樣品取自甘肅嘉峪關(guān)，為棕漠土（正常干旱土），土壤pH8.34，有機(jī)質(zhì)含量1.7 g kg-1，全氮、全磷、全鉀含量分別為94 mg kg-1、315 mg kg-1、11.6 g kg-1，容重1.56 cm3g-1。土壤質(zhì)地為砂土，其中砂粒含量98.44%、粉粒含量1.01%、黏粒含量0.55%。土壤礦物、化學(xué)成分見(jiàn)表1。

1.2 試驗(yàn)設(shè)計(jì)

采用室內(nèi)土柱法，根據(jù)噴淋強(qiáng)度不同設(shè)置3 個(gè)試驗(yàn)土柱，編號(hào)為ST-1、ST-2、ST-3，對(duì)應(yīng)噴淋強(qiáng)度分別為120 ml d-1、60 ml d-1、30 ml d-1。試驗(yàn)裝置如圖1，由三部分組成：（1）噴淋裝置，包括供水裝置和用于控制流量的蠕動(dòng)泵；（2）土柱裝置；（3）底部流出物收集裝置。試驗(yàn)用水為蒸餾水，噴淋裝置由蠕動(dòng)泵控制注入流量。土柱體為有機(jī)玻璃材質(zhì)，厚度0.5 cm，高85 cm，橫截面內(nèi)徑為11 cm的圓形。管的頂部、底部各留5 cm的空間，便于頂部噴淋水的下滲和防止底部流出液的淤積。

表1 供試土樣礦物、化學(xué)成分分析Table 1 Mineral and chemical composition of test soil（%）

圖1 土柱試驗(yàn)裝置Fig. 1 Soil column experiment device

土柱填裝時(shí)，在距管底部5 cm高度處的多孔板上放置一張紗網(wǎng)，將土樣分10 層裝入土柱，每層裝填高度6 cm，總裝填高度60 cm，裝填量由土樣的容重和土柱的體積計(jì)算得到，在每?jī)蓪油林g，打毛土壤界面。蠕動(dòng)泵控制蒸餾水流量，分別對(duì)三個(gè)試驗(yàn)土柱自上而下進(jìn)行噴淋，待土柱底部流出量達(dá)到穩(wěn)定后鋪放示蹤源層。源層用放射性核素Sr-90（總量1.24×106Bq）與80 目的石英砂（約75 g）混合均勻后平鋪在土柱體上端面，示蹤源層厚約3 mm。最后在示蹤源層上鋪放14.7 cm的聚乙烯細(xì)顆粒，目的是使噴淋水能均勻下滲，且不吸附Sr-90。

ST-1、ST-2、ST-3土柱噴淋時(shí)間分別為51 d、102 d、204 d，試驗(yàn)過(guò)程中每天在底部收集流出液并測(cè)量體積和Sr-90的濃度。試驗(yàn)周期的確定原則上是保證Sr-90遷移曲線(xiàn)的濃度峰不穿透土柱體，而在土柱中有完整的濃度分布曲線(xiàn)。試驗(yàn)結(jié)束后，以源層為零點(diǎn)，向下為正，將土柱解體為每層3 mm的樣品，并測(cè)量Sr-90的濃度。具體為先用6 mol L-1的鹽酸浸取砂土，草酸沉淀鈣、鎂、鐵等金屬離子，再用3 mol L-1的硝酸溶解沉淀，過(guò)Srspec樹(shù)脂柱，最后用蒸餾水將Sr-90洗脫，草酸沉淀Sr2+，質(zhì)量法計(jì)算Sr的化學(xué)回收率，溶解沉淀并用低本底液體閃爍譜儀測(cè)量Sr-90的濃度。另外，取6、12、18、24、30、36、42、48、54、60 cm處的樣品，用烘干法測(cè)量含水率，得到土柱體垂向上的含水率分布情況，作為實(shí)測(cè)資料用于核素遷移水分模型的驗(yàn)證數(shù)據(jù)。

1.3 土壤水力參數(shù)測(cè)定

土壤水力參數(shù)包括土壤水分特征曲線(xiàn)與飽和導(dǎo)水率。土壤水分特征曲線(xiàn)采用砂性漏斗法和機(jī)械張力計(jì)法相結(jié)合的方式測(cè)定，土壤吸力較小時(shí)，采用砂性漏斗法；土壤吸力較大時(shí)，采用張力計(jì)法。通過(guò)試驗(yàn)得到了8 個(gè)不同體積含水率下的土壤吸力值。采用Van Genuchten（1980）模型對(duì)飽和含水率θs、殘余含水率θr、進(jìn)氣值倒數(shù)α和擬合參數(shù)n進(jìn)行了求解。參數(shù)結(jié)果如表2。

此外由于砂土介質(zhì)的滲透性較好，飽和導(dǎo)水率測(cè)定采用定水頭法。從土柱底端通水，使整個(gè)土柱體飽和，上部安裝定水頭裝置（如圖2），調(diào)整入水口供水量，使試驗(yàn)期間水頭差H保持不變，記錄時(shí)間Δt和底部出口流量Q。

表2 土壤水力參數(shù)Table 2 Soil hydraulic parameters

2 結(jié) 果

2.1 不同噴淋強(qiáng)度土柱含水率

三個(gè)土柱體中含水率分布從上到下逐漸增大，越靠近底部越接近飽和含水率。ST-1土柱含水率為0.167～0.324 cm3cm-3、ST-2土柱含水率為0.172～0.334 cm3cm-3、ST-3土柱含水率為0.170～0.332 cm3cm-3（表3）。分析原因是，砂土介質(zhì)的滲透系數(shù)遠(yuǎn)大于噴淋強(qiáng)度，水迅速下滲并在土柱底部累積，當(dāng)達(dá)到飽和含水率后在重力作用下排出。雖然噴淋量強(qiáng)度不同，但是當(dāng)含水率達(dá)到穩(wěn)定后，三個(gè)土柱體中垂向含水率的分布基本相同，區(qū)別在于噴淋量越大，含水率越快達(dá)到穩(wěn)定。

圖2 飽和導(dǎo)水率測(cè)定試驗(yàn)裝置Fig. 2 Test device for measuring saturated water conductivity

2.2 不同噴淋強(qiáng)度下Sr-90在土柱中的遷移

ST-1土柱噴淋51 d后，Sr-90濃度峰垂直向下遷移約42.3 cm，峰值濃度為415.8 Bq cm-3（圖3）。ST-2土柱噴淋102 d后，Sr-90濃度峰向下遷移約46.2 cm，峰值濃度為289.6 Bq cm-3。ST-3土柱噴淋204 d后，Sr-90濃度峰向下遷移約44.4 cm，峰值濃度為323.9 Bq cm-3。

三種噴淋強(qiáng)度下Sr-90的濃度峰在砂土介質(zhì)中的遷移距離基本一致，遷移距離與噴淋強(qiáng)度、時(shí)間成正比。但是三個(gè)土柱上部0～40 cm深度范圍內(nèi)Sr-90的濃度并沒(méi)有很快遞減為零，依然保持約100～200 Bq cm-3，濃度分布曲線(xiàn)存在拖尾的特征，說(shuō)明Sr-90從源層上被洗脫下來(lái)后，被吸附并滯留至砂土介質(zhì)上，而濃度峰值通過(guò)后，更多的Sr-90從砂土介質(zhì)中解吸出來(lái)，使得很長(zhǎng)時(shí)間段內(nèi)土柱體上部保持相對(duì)較高的濃度。ST-2、ST-3土柱中Sr-90的濃度分布曲線(xiàn)和峰值濃度基本一致，而ST-1土柱明顯大于前兩者，意味著噴淋強(qiáng)度越大，Sr-90更容易從源層上洗脫下來(lái)，而當(dāng)噴淋強(qiáng)度降低到60 ml d-1時(shí)，影響變小。從污染場(chǎng)地的修復(fù)治理和劑量估算的角度考慮，當(dāng)Sr-90進(jìn)入砂土介質(zhì)后，高濃度污染將會(huì)是一個(gè)持續(xù)過(guò)程。

表3 實(shí)測(cè)土柱體含水率分布Table 3 Measured moisture content in soil column

圖3 土柱體剖面Sr-90實(shí)測(cè)濃度分布曲線(xiàn)Fig. 3 Concentration distribution curve of Sr-90 in soil column profile

3 核素遷移的數(shù)值模擬

3.1 水流模型的建立及含水率驗(yàn)證

為了更加深入地分析上述現(xiàn)象，首先應(yīng)用Hydrus-1d軟件建立垂向上的一維飽和—非飽和水流模型，然后分別加入平衡吸附遷移模式和非平衡吸附遷移模式。從遷移模式和遷移參數(shù)兩個(gè)方面探討不同噴淋強(qiáng)度對(duì)核素Sr-90在砂土中遷移的影響。

首先，根據(jù)Richard方程建立垂向上的一維飽和—非飽和水流模型，空間上將模型離散為3 mm的單元格，時(shí)間步長(zhǎng)為90 s。將試驗(yàn)得到的土壤物理參數(shù)、土壤水力參數(shù)、噴淋強(qiáng)度、底部流出量代入模型，計(jì)算得到土柱體垂向上含水率分布值，同樣選取6、12、18、24、30、36、42、48、54、60 cm處的計(jì)算值。目的是與試驗(yàn)得到的實(shí)測(cè)含水率值對(duì)比，驗(yàn)證水流模型的可靠性。

式中，θ為體積含水率（cm3cm-3），t為水分運(yùn)移時(shí)間（d），h為非飽和帶壓力水頭（cm），z為空間坐標(biāo)，K（h）為非飽和導(dǎo)水率（cm d-1），θ0為初始含水率（cm3cm-3），V1為入滲通量（cm d-1），V2為底部流出量（cm d-1），H為土柱體高度（cm）。

從含水率計(jì)算結(jié)果來(lái)看（表4），模型中含水率分布也是由上至下逐漸增大，靠近底部接近飽和含水率。ST-1土柱含水率為0.187～0.317 cm3cm-3，與實(shí)測(cè)值的偏差為-5.58%～11.98%。ST-2土柱含水率為0.156～0.298 cm3cm-3，與實(shí)測(cè)值的偏差為-12.80%～8.81%。ST-3土柱含水率為0.148～0.290 cm3cm-3，與實(shí)測(cè)值的偏差為-12.94%～8.90%。對(duì)比結(jié)果發(fā)現(xiàn)兩者基本吻合，所建水流模型精度基本達(dá)到要求，能反映實(shí)際水流運(yùn)動(dòng)情況。

表4 模型含水率計(jì)算結(jié)果及驗(yàn)證Table 4 Calculation results and verification of moisture content by model

3.2 平衡與非平衡吸附模式的對(duì)比及遷移參數(shù)擬合

目前，核素在土壤中的遷移模式分為平衡吸附模式和非平衡吸附模式。平衡吸附模式，意味著核素在固相介質(zhì)中能很快達(dá)到吸附、解吸平衡或者是有足夠長(zhǎng)的時(shí)間讓吸附相與溶解相之間達(dá)到平衡，但實(shí)際上對(duì)于吸附能力弱或者沒(méi)有足夠長(zhǎng)時(shí)間與固相介質(zhì)接觸的核素，其濃度分布曲線(xiàn)往往存在不對(duì)稱(chēng)性和拖尾。而非平衡吸附模式能更好地描述這種現(xiàn)象。

20世紀(jì)50年代，對(duì)流—彌散理論開(kāi)始用于描述溶質(zhì)的遷移運(yùn)動(dòng)，其遷移行為是物理、化學(xué)等一系列復(fù)雜過(guò)程綜合作用的結(jié)果［16-20］。Yeh和Tripathi［21］以及Mangold和Tsang［22］將對(duì)流—彌散和多組分化學(xué)反應(yīng)耦合建立了溶質(zhì)的非平衡遷移方程，它有效地融入了化學(xué)非平衡和物理非平衡兩種過(guò)程?；瘜W(xué)非平衡過(guò)程考慮了不同形式的吸附作用，物理非平衡過(guò)程考慮了不同形式的水流運(yùn)動(dòng)，盡管基于理論的不同，但是在一定的假設(shè)前提下，化學(xué)非平衡過(guò)程和物理非平衡過(guò)程可以耦合成統(tǒng)一的數(shù)學(xué)形式，為非平衡吸附模式的發(fā)展和應(yīng)用創(chuàng)造了條件。

平衡吸附模型［13-15］認(rèn)為核素在固相介質(zhì)中能很快達(dá)到吸附平衡，此時(shí)固相中核素濃度是分配系數(shù)與液相中核素濃度的乘積，數(shù)學(xué)模型如下：

非平衡吸附模型［8-12］（單點(diǎn)吸附模型）與平衡吸附模式的區(qū)別是，遷移方程中引入了一級(jí)速率系數(shù)β，認(rèn)為核素在遷移過(guò)程中被固相介質(zhì)吸附不會(huì)瞬時(shí)完成，而是存在吸附速率的問(wèn)題。數(shù)學(xué)模型如下：

初始條件和邊界條件如下：

式中，θ為體積含水率（cm3cm-3），t為核素運(yùn)移時(shí)間（d），z為空間坐標(biāo)，D為彌散系數(shù)（cm2d-1），ρ為容重（g cm-3），C為溶解相濃度（Bq cm-3），Sk為吸附相濃度（Bq cm-3），Kd為平衡吸附時(shí)的分配系數(shù)（cm3g-1），β為溶解相與吸附相之間的一級(jí)速率系數(shù)（d-1），H為土柱體高度（cm），L為源層厚度（cm），C0為初始濃度（Bq cm-3），C1為流出液濃度（Bq cm-3）。

分別建立平衡吸附、非平衡吸附（單點(diǎn)吸附模型）兩種模式下的核素遷移一維垂直入滲數(shù)值模型。計(jì)算得到兩種模式下Sr-90在土柱體中的計(jì)算濃度分布曲線(xiàn)，將實(shí)測(cè)分布曲線(xiàn)與計(jì)算分布曲線(xiàn)進(jìn)行對(duì)比分析，根據(jù)擬合程度，求得Sr-90在砂土中的分配系數(shù)Kd、彌散度D、一級(jí)速率系數(shù)β。

從濃度計(jì)算結(jié)果來(lái)看（表5），ST-1土柱平衡吸附模式的計(jì)算峰值約634.3 Bq cm-3，非平衡吸附模式的計(jì)算峰值約396.4 Bq cm-3，分配系數(shù)為0.73 ml g-1，β為1.46 d-1。ST-2土柱平衡吸附模式的計(jì)算峰值約523.2 Bq cm-3，非平衡吸附模式的計(jì)算峰值約253.3 Bq cm-3，分配系數(shù)為0.70 ml g-1，β為0.56 d-1。ST-3土柱平衡吸附模式的計(jì)算峰值約695.2 Bq cm-3，非平衡吸附模式的計(jì)算峰值約317.5 Bq cm-3，分配系數(shù)為0.72 ml g-1，β為0.36 d-1。

平衡吸附模式計(jì)算的濃度分布曲線(xiàn)呈明顯的對(duì)稱(chēng)性，高于實(shí)測(cè)濃度，且在濃度峰值通過(guò)后，很快遞減為零，無(wú)法描述實(shí)測(cè)濃度分布曲線(xiàn)的拖尾特征。非平衡吸附模式考慮了溶解相與吸附相之間的一級(jí)速率系數(shù)β，計(jì)算的濃度分布曲線(xiàn)與實(shí)測(cè)濃度點(diǎn)吻合程度較好（圖4）。3 種噴淋強(qiáng)度下，Sr-90的濃度峰在土柱體中的遷移距離基本一致，而非平衡吸附模式中所定義的分配系數(shù)為線(xiàn)性平衡吸附時(shí)的分配系數(shù)，按照理論只要介質(zhì)相同，其值應(yīng)該一致。本次試驗(yàn)擬合得到的分配系數(shù)為0.70～ 0.73 ml g-1，基本上變化不大，符合理論，同時(shí)說(shuō)明Sr-90在砂土介質(zhì)中的遷移距離與噴淋強(qiáng)度和時(shí)間成比例。核素在土壤介質(zhì)中的分配系數(shù)受pH和有機(jī)質(zhì)含量的影響較大，但是本次三個(gè)試驗(yàn)土柱的土樣來(lái)自同一種土壤介質(zhì)，經(jīng)過(guò)測(cè)量土壤中pH、有機(jī)質(zhì)含量、氮、磷、鉀含量基本上相同。ST-2土柱砂土有機(jī)質(zhì)含量較ST-1、ST-3土柱少約16%，因此該土柱Sr-90濃度峰遷移的距離較其他2個(gè)土柱稍遠(yuǎn)，分配系數(shù)稍小。

三種噴淋強(qiáng)度下，三個(gè)土柱體的彌散度為0.08～0.24 cm。這是由于本次試驗(yàn)采用的是擾動(dòng)土，裝填存在非均質(zhì)性導(dǎo)致。

表5 Sr-90遷移參數(shù)Table 5 Migration parameters of Sr-90

4 結(jié) 論

在實(shí)驗(yàn)室內(nèi)建立了120 ml d-1、60 ml d-1、30 ml d-1三種不同噴淋強(qiáng)度下Sr-90在砂土土柱中的遷移試驗(yàn)（對(duì)應(yīng)試驗(yàn)時(shí)間分別為51 d、102 d、204 d），對(duì)比分析了不同噴淋強(qiáng)度對(duì)核素Sr-90在砂土中遷移的影響。得出以下主要結(jié)論：（1）三種噴淋強(qiáng)度下，Sr-90濃度峰的遷移距離基本一致，遷移距離與噴淋強(qiáng)度和時(shí)間成正比。（2）三組試驗(yàn)的噴淋強(qiáng)度和試驗(yàn)時(shí)間均不同，但是噴淋的總水量一致，說(shuō)明淋洗水量對(duì)Sr-90的遷移更為重要。（3）三種噴淋強(qiáng)度下，Sr-90的濃度分布曲線(xiàn)均明顯存在拖尾，濃度峰值過(guò)后，更多的Sr-90會(huì)從砂土介質(zhì)中被解吸出來(lái)，使得很長(zhǎng)時(shí)間段內(nèi)保持較高的濃度。（4）應(yīng)用平衡吸附模式已經(jīng)無(wú)法對(duì)Sr-90濃度分布曲線(xiàn)的拖尾進(jìn)行準(zhǔn)確描述，而采用非平衡吸附模式，引入一級(jí)速率系數(shù)β后，能更好地描述該現(xiàn)象。

圖4 土柱Sr-90濃度擬合計(jì)算Fig. 4 Fitting of Sr-90 concentration

［1］ 李靜，張美一，潘綱，等. 不同吸附模式對(duì)As（V）在TiO2顆粒上亞穩(wěn)態(tài)吸附的影響. 物理化學(xué)學(xué)報(bào)，2013，29（7）：1541—1549

Li J，Zhang M Y，Pan G，et al. Influence of adsorption model on metastable equilibrium adsorption of As（V）on TiO2particles（In Chinese）. Acta Physico-Chimica Sinica，2013，29（7）：1541—1549

［2］ 徐叢，李薇，潘綱. Zn（II）/α-FeOOH吸附體系的固體濃度效應(yīng). 物理化學(xué)學(xué)報(bào)，2009，25（9）：1737—1742

Xu C，Li W，Pan G. Particle concentration effect on Zn（II）adsorption at water-goethite interfaces（In Chinese）. Acta Physico-Chimica Sinica，2009，25 （9）：1737—1742

［3］ Cho W J，Lee J O，Choi H J. Radionuclide migration through an unsaturated clay buffer under thermal and hydraulic gradients for a nuclear waste repository. Annals of Nuclear Energy，2012，50（4）：71—81

［4］ Merk R. Numerical modeling of the radionuclide water pathway with HYDRUS and comparison with the IAEA model of SR 44. Journal of Environmental Radioactivity，2012，105（2）：60—69

［5］ 于靜，司高華，劉旭東，等. 噴淋強(qiáng)度對(duì)Sr在沙土中遷移的影響研究. 甘肅科技，2011，27（9）：47—49

Yu J，Si G H，Liu X D，et al. Study on water flow rate effect on migration of strontium in sand（In Chinese）. Gansu Science and Technology，2011，27（9）：47—49

［6］ 王志明，郭謙，李書(shū)紳，等. 水流速度對(duì)鍶遷移影響的實(shí)驗(yàn)研究. 原子能科學(xué)技術(shù)，2003，37（4）：320—324

Wang Z M，Guo Q Li S S，et al. Experiment on effect of water flow velocity on migration of strontium（In Chinese）. Atomic Energy Science and Technology，2003，37（4）：320—324

［7］ Li Y F，Tian S H，Qian T W. Transport and retention of strontium in surface-modified quartz sand with different wettability. Journal of Radioanalytical and Nuclear Chemistry，2011，289（2）：337—343

［8］ Abulaban A，Nieber J L. Modeling the effects of nonlinear equilibrium sorption of solute plumes in saturated heterogeneous porous media. Advances in Water Resources，2000，23（8）：893—905

［9］ Leij F J，Toride N，van Genuchten Th. Analytical solutions for non-equilibrium solute transport in threedimensional porous media. Journal of Hydrology，1993，151（2/4）：193—228

［10］ 王金生，李書(shū)紳，王志明，等. 核素遷移的二維非平衡吸附數(shù)值模式研究. 應(yīng)用生態(tài)學(xué)報(bào)，1996，7（1）：89—93

Wang J S，Li S S，Wang Z M，et al. Two-dimensional non-equilibrium sorption numerical model of radionuclide migration（In Chinese）. Chinese Journal of Applied Ecology，1996，7（1）：89—93

［11］ 陳家軍，戴杰，郭擇德. 非平衡吸附放射性核素遷移數(shù)值模型研究. 水文地質(zhì)工程地質(zhì)，1994，21（2）：37—40

Chen J J，Dai J，Guo Z D. Non-equilibrium adsorption numerical model for radionuclide migration research（In Chinese）. Hydrogeology and Engineering Geology，1994，21（2）：37—40

［12］ 陳家軍，王金生，李合蓮. 非平衡吸附對(duì)核素遷移的影響及平衡速率參數(shù)的確定. 中國(guó)環(huán)境科學(xué)，2000，20 （1）：73—76

Chen J J，Wang J S，Li H L. Influence on radionuclide migration by non-equilibrium sorption and identification of rate coefficient（In Chinese）. China Environmental Science，2000，20（1）：73—76

［13］ 王志明，李書(shū)紳，楊月娥. 非飽和黃土介質(zhì)中含水量對(duì)85Sr遷移的影響. 原子能科學(xué)技術(shù)，2001，35（4）：320—324

Wang Z M，Li S S，Yang Y E. Effect of water content on migration of85Sr in unsaturated loess（In Chinese）. Atomic Energy Science and Technology，2001，35 （4）：320—324

［14］ 王志明，安永鋒.85Sr在非飽和黃土中的遷移特征. 原子能科學(xué)技術(shù)，2004，38（1）：29—34

Wang Z M，An Y F. Migration characteristics of85Sr in unsaturated loess（In Chinese）. Atomic Energy Science and Technology，2004，38（1）：29—34

［15］ 王志明，李書(shū)紳，郭擇德，等.85Sr在非飽和黃土中的遷移. 輻射防護(hù)，2000，20（1/2）：32—35

Wang Z M，Li S S，Guo Z D，et al. Migration of85Sr in unsaturated Loess（In Chinese）. Radiation Protection，2000，20（1/2）：32—35

［16］ 王康，張仁鐸，周祖昊，等. 土壤非均勻水流運(yùn)動(dòng)和溶質(zhì)遷移顯色示蹤方法研究. 土壤學(xué)報(bào)，2009，46 （4）：603—610

Wang K，Zhang R D，Zhou Z H，et al. Evaluation of method for visualization of heterogeneous soil water flow and solute transport using IODINE-STARCH staining （In Chinese）. Acta Pedologica Sinica，2009，46 （4）：603—610

［17］ 陶亞奇，蔣新，卞永榮，等. 二維虛擬土壤中溶質(zhì)遷移行為的數(shù)值模擬研究. 土壤學(xué)報(bào)，2009，46（6）：991—997

Tao Y Q，Jiang X，Bian Y R，et al. Numerical simulation of solute transports in two dimensional virtual soil（In Chinese）. Acta Pedologica Sinica，2009，46 （6）：991—997

［18］ 馬美景，王軍光，郭忠錄，等. 放水沖刷對(duì)紅壤坡面侵蝕過(guò)程及溶質(zhì)遷移特征的影響. 土壤學(xué)報(bào)，2016，53 （2）：365—374

Ma M J，Wang J G，Guo Z L，et al. Research on sediment and solute transport on red soil slope under simultaneous influence of scouring flow（In Chinese）. Acta Pedologica Sinica，2016，53（2）：365—374

［19］ 鄭紀(jì)勇，王全九，邵明安. 利用溶質(zhì)鋒信息估計(jì)溶質(zhì)遷移參數(shù)方法分析. 土壤學(xué)報(bào)，2007，44（2）：221—226

Zheng J Y，Wang Q J，Shao M A. Estimation of solute transport parameters using information of solute front （In Chinese）. Acta Pedologica Sinica，2007，44 （2）：221—226

［20］ 任理，秦耀東，王濟(jì). 非均質(zhì)飽和土壤鹽分優(yōu)先運(yùn)移的隨機(jī)模擬. 土壤學(xué)報(bào)，2001，38（1）：104—113 Ren L，Qin Y D，Wang J. Stochastic modeling ofchloride travel in heterogeneous saturated soils under the condition of preferential flow（In Chinese）. Acta Pedologica Sinica，2001，38（1）：104—113

［21］ Yeh G T，Tripathi V S. A critical evaluation of recent developments in hydrogeochemical transport models of reactive multichemical components. Water Resource Research，1989，25（1）：93—108

［22］ Mangold D C，Tsang C F. A summary of subsurface hydrological and hydrochemical models. Reviews of Geophysics，1991，29（1）：51—79

Effect of Spraying Intensity on Migration of Sr-90 in Soil

ZHU Jun DENG Anchang SHI Yunfeng CHEN Chao LIU Tuantuan ZHANG Aiming?
（China Institute for Radiation Protection，Taiyuan 030006，China）

Based on the dynamic soil column method in laboratory，experiment was conducted on studying the migration of radionuclide Sr-90 in Sand relative to different spraying intensity（spraying intensity was designed to 120 ml d-1，60 ml d-1，30 ml d-1，and the corresponding experimental time were 51，102 204 day respectively）.Results showed that migration distance of Sr-90 concentration peak in the soil columns more or less the same，were 42.3 cm，46.2 cm，44.4 cm，Which was positively related to intensity and duration of the spraying. In addition，the concentration distribution curves of Sr-90 in soil column section all exhibited asymmetry and“tailing”phenomena under the three different spray intensity，i.e. After Sr-90 was leached down from the source layer，it was adsorbed to the Sand，thus lagging behind and getting desorbed after the peak concentration passing the observation point，hence Sr-90 concentration remained to be quite high in soil columns for a long time. With the aid of HYDRUS-1D software，two numerical models were established for radionuclide migration，one based on equilibrium adsorption and the other on nonequilibrium adsorption. It was found that the non-equilibrium sorption model，which took into account the first-order rate coefficients β，could be used to better explain the asymmetry and tailing phenomena. Based on fitting the measurements of Sr-90 concentrations in the soil columns，distribution coefficient of Sr-90 in the Sand was figured out to be 0.73 ml g-1，0.70 ml g-1and 0.72 ml g-1relative to intensity of the spraying，which shows that intensity of the spraying did not have much effect on distribution coefficient of Sr-90.

Spraying intensity；Soil column experiment；Radionuclide migration；Equilibrium adsorption；Non-equilibrium adsorption

X53

A

（責(zé)任編輯：盧 萍）

* 中國(guó)輻射防護(hù)研究院院基金項(xiàng)目 Supported by the Foundation of China Institute for Radiation Protection

? 通訊作者 Corresponding author，E-mial：13834508698@163.com

朱 君（1985—），男，湖南衡陽(yáng)人，碩士，研究實(shí)習(xí)員，主要從事核素在地表水環(huán)境、地下水環(huán)境、土壤環(huán)境中的吸附和遷移研究。E-mial：405400881@qq.com

2016-08-16；

2016-12-18；優(yōu)先數(shù)字出版日期（www.cnki.net）：2017-03-03

10.11766/trxb201608160166