李子超， 周濤， 石順， 秦雪猛

(1.華北電力大學(xué) 核科學(xué)與工程學(xué)院，北京 102206；2.東南大學(xué) 能源與環(huán)境學(xué)院，江蘇 南京 211096)

傳統(tǒng)的化石燃料導(dǎo)致了嚴(yán)重的環(huán)境問(wèn)題，已經(jīng)不能滿足可持續(xù)發(fā)展的需求，發(fā)展核電是解決這一問(wèn)題的有效途徑。然而核電發(fā)展也有其局限性，一旦發(fā)生核泄漏事故，將會(huì)造成極大的影響。日本福島核電站受地震海嘯影響，發(fā)生嚴(yán)重核泄漏事故，造成不可估量的損失[1]。目前我國(guó)核電廠全部分布在沿海地區(qū)，例如2018年10月22日正式投入商運(yùn)的海陽(yáng)核電站AP1000機(jī)組，緊鄰半封閉海域的黃海，一旦發(fā)生核泄漏事故，放射性核素很可能通過(guò)海洋擴(kuò)散污染生態(tài)環(huán)境。研究放射性核素在海洋中的遷移機(jī)制，有效應(yīng)對(duì)海洋核污染已經(jīng)成為一個(gè)重要議題，國(guó)內(nèi)外學(xué)者對(duì)放射性核素在海洋中的遷移進(jìn)行了系列研究。福島核事故后，Buesseler等[2]分析了核電站附近海域放射性核素濃度，發(fā)現(xiàn)事故發(fā)生4個(gè)月后，137Cs的含量比事故前高近萬(wàn)倍。Nakano等[3]用全球輸運(yùn)模式對(duì)福島泄漏核素在海洋中的長(zhǎng)期遷移進(jìn)行了數(shù)值模擬，分析了核素的遷移路徑。王輝等[4]對(duì)福島核泄漏物質(zhì)在海洋中的運(yùn)輸進(jìn)行了10 a的模擬和預(yù)測(cè)。Kanai等[5]指出核泄漏事故后大氣中的放射性核素輸運(yùn)較快，幾個(gè)月之內(nèi)便可因快速擴(kuò)散和沉降迅速減少。此外，文獻(xiàn)[6]估計(jì)了大氣中放射性核素的沉降作用，研究表明大氣中87%的131I和78%的137Cs沉降到了海洋中。喬清黨等[7]分析了放射性核素在海洋中遷移的數(shù)值模擬方法，給出了建立事故后果評(píng)價(jià)系統(tǒng)的技術(shù)路線，指出擬建立的系統(tǒng)只考慮了核素衰變影響，沒(méi)有考慮大氣沉降、懸浮物吸附沉降等因素影響。放射性核素在海洋中的遷移過(guò)程十分復(fù)雜，受到源項(xiàng)、對(duì)流擴(kuò)散、大氣沉降、自身衰變、懸浮物吸附沉降和生物富集等多種因素影響。目前對(duì)核泄漏事故后放射性核素在海洋中遷移的研究已有不少探討，但較少考慮懸浮物吸附沉降、大氣核素入海通量和生物富集等因素對(duì)核素遷移的影響。因此，本文通過(guò)分析放射性核素在海洋中的遷移機(jī)制，建立了相應(yīng)模型，研究了不同因素對(duì)核素?cái)U(kuò)散的影響程度。

1 AP1000機(jī)組環(huán)境

1.1 源項(xiàng)

海陽(yáng)核電站規(guī)劃建設(shè)6臺(tái)百萬(wàn)千瓦級(jí)壓水堆機(jī)組，留有2臺(tái)擴(kuò)建余地，單機(jī)容量1 250 MW，熱功率為3 415 MWt。首臺(tái)機(jī)組2018年8月并網(wǎng)發(fā)電，因此，選擇海陽(yáng)核電站AP1000機(jī)組進(jìn)行分析。根據(jù)美國(guó)反應(yīng)堆安全研究報(bào)告[8]和福島核事故放射性核素在海洋中的泄漏情況[9]，設(shè)定131I大氣泄漏總量為1016Bq，137Cs的大氣泄漏總量為1016Bq，10 d內(nèi)均勻泄漏，泄漏率均為1.16×1010Bq/s。

1.2 大氣環(huán)境

放射性核素的干沉降和濕沉降受多種環(huán)境因素影響，其中顆粒物粒徑對(duì)放射性核素沉降有重要影響，切爾諾貝利核泄漏事故釋放的放射性氣溶膠中，50%到80%的粒徑小于1.1 μm[10]。目前干沉降速率還沒(méi)有確定值，因此根據(jù)德國(guó)輻射防護(hù)委員會(huì)第十七卷出版物(SSK 17號(hào)報(bào)告)內(nèi)容進(jìn)行選取。沖刷比例常數(shù)α值根據(jù)非穩(wěn)態(tài)拉格朗日煙團(tuán)模型系統(tǒng)CALPUFF模式選取，各情景參數(shù)值如表1所示。

表1 干、濕沉降環(huán)境參數(shù)Table 1 Parameters of dry and wet sedimentation

1.3 海洋懸浮物環(huán)境

海洋中懸浮物對(duì)核素的吸附沉降受懸浮物沉積通量、分配系數(shù)、水深和時(shí)間等多種因素影響。黃海平均水深為44 m，海底西部平坦，東部復(fù)雜多變，海陽(yáng)核電站位于南黃海西部，平均水深可設(shè)定為44 m。關(guān)于中國(guó)近海岸懸浮物沉積通量已有不少研究，北黃海的沉積通量為 0.06～1.18 g/(cm2·a)[11]；渤、黃、東海沉積通量0.2 g/(cm2·a)[12]；南黃海的沉積通量為 0.033～0.76 g/(cm2·a)[13]。關(guān)于懸浮物對(duì)核素的吸附能力，銫分配系數(shù)小于500 L/kg，碘幾乎不能被懸浮物吸附，分配系數(shù)很小[14]?；谝陨衔墨I(xiàn)條件，懸浮物吸附沉降參數(shù)設(shè)定如表2所示。

表2 懸浮物吸附沉降參數(shù)Table 2 Parameters of sorption and sedimentation by suspended solids

1.4 海洋生物環(huán)境

關(guān)于海洋生物對(duì)放射性核素的富集作用，銫的分配系數(shù)小于50 L/kg，而碘的分配系數(shù)小于103L/kg。黃海生物量密度[15]在1～103kg/km2，具體參數(shù)設(shè)定如表3所示。

表3 生物富集參數(shù)Table 3 Parameters of bioaccumulation

2 影響程度分析模型

2.1 衰變模型

放射性活度和衰變常數(shù)關(guān)系：

A=A0e-λt

(1)

式中：A0為原來(lái)活度，Bq；t為時(shí)間，s；A為經(jīng)過(guò)t時(shí)間衰變后的活度，Bq；λ為衰變常數(shù)。

2.2 大氣核素沉降模型

2.2.1 大氣核素?cái)U(kuò)散

事故情況下，人們關(guān)心的是放射性煙云的遷移方向和沉降區(qū)域，因此計(jì)算快速的高斯煙羽模式[16]很適合應(yīng)急反應(yīng)。高斯煙羽模式是計(jì)算大氣中氣載污染物下風(fēng)向濃度中應(yīng)用最廣泛的模式，此模式假定煙羽中污染物濃度在水平方向和垂直方向都遵循高斯分布，高斯煙羽模式：

(2)

2.2.2 干沉降

假定氣載放射性物質(zhì)的沉積通量密度正比于近地面處氣載放射性核素濃度，則干沉降通量模型：

Wd=VdC(x,y,0)

(3)

式中：Wd為干沉積通量密度，Bq/(s·m2)；C(x,y,0)為地面處放射性核素密度，Bq/m3；Vd為沉積速度，m/s。

2.2.3 濕沉降

煙羽中的污染物濃度除干沉降外，還會(huì)因降水沖刷減小，污染物的濕沉降通量模型：

(4)

式中：Wd為濕沉積通量密度，Bq/(s·m2)；Λ表示沖洗系數(shù)，1/s；C(x0,y0,z)為(x0,y0)處高度為z的核素濃度，Bq/m3。

2.3 海洋懸浮物吸附沉降模型

海洋中懸浮物來(lái)源很廣，主要來(lái)源于河流攜帶、大氣灰塵、生物顆粒和地殼風(fēng)化物等過(guò)程，不是一個(gè)靜止?fàn)顟B(tài)，物質(zhì)交換過(guò)程十分復(fù)雜。提出了懸浮物對(duì)放射性核素吸附沉降影響的箱式模型，考慮到懸浮物的更新過(guò)程，做出如下假設(shè)：

1)分析范圍內(nèi)核素濃度均勻一致；

2)進(jìn)入水中的懸浮物和沉積的懸浮物數(shù)量相等，處于動(dòng)態(tài)平衡狀態(tài)；

3)懸浮物對(duì)放射性核素的吸收量很短時(shí)間內(nèi)可以達(dá)到穩(wěn)定狀態(tài)，吸收時(shí)間忽略不計(jì)；

4)分析海域內(nèi)水深為一平均值h；

5)Δt時(shí)間內(nèi)水環(huán)境放射性核素濃度不變。

建立了海洋中核素動(dòng)態(tài)吸附的箱式模型，懸浮物作用后水中放射核素活度計(jì)算模型：

(5)

式中：C為懸浮物吸附沉降后水中放射性活度，Bq/L；C0為懸浮物吸附沉降前水環(huán)境放射性活度，Bq/L；φ為懸浮物沉積通量，g/(cm2·a)；Δt為間隔時(shí)間，a，當(dāng)Δt=0.01 a時(shí)，達(dá)到收斂精度要求；h為水深，m；Kd為分配系數(shù)L/kg。

2.4 生物富集模型

海洋中有浮游動(dòng)物、浮游植物和魚(yú)類等生物，會(huì)通過(guò)表面吸附、體內(nèi)吸收等多種方式富集核素，因此生物是富集核素的一個(gè)重要因素，為建立海洋生物富集模型，做出如下假設(shè)：

1)分析范圍內(nèi)核素濃度均勻一致；

2)分析范圍內(nèi)水中生物量分布均勻；

3)海洋生物總量處于動(dòng)態(tài)平衡，且吸收的核素在生物圈內(nèi)循環(huán)。

生物富集后水中放射性核素濃度計(jì)算模型為：

(6)

式中：C0為生物富集前海水中放射性活度，Bq/L；C為生物富集后海水中放射性核素濃度，Bq/L；Fd為分配系數(shù)，L/kg；Wb為水中生物量，kg/L。

3 計(jì)算結(jié)果及分析

3.1 衰變影響

根據(jù)2.1節(jié)中的模型，選擇半衰期為30.17 a的137Cs和半衰期為8.04 d的131I，分別計(jì)算了其放射性活度隨時(shí)間衰變后的剩余份額，如圖1所示。

圖1 衰變后剩余份額Fig.1 Remaining share after decay

從圖1可以看出，131I隨時(shí)間衰變較快，一個(gè)月之內(nèi)就可以衰變?yōu)樵瓉?lái)的十分之一，因此，對(duì)于短期核素?cái)U(kuò)散，需要考慮其衰變影響。137Cs半衰期較長(zhǎng)，對(duì)于短期核素?cái)U(kuò)散，可以忽略衰變影響。

3.2 干沉降影響

核泄漏事故發(fā)生后，氣溶膠粒子、元素碘和化合物碘混合在一起，很難精確測(cè)量各組分所占比例，因此選擇最小干沉降速率5×10-4m/s和最大干沉降速率4×10-3m/s進(jìn)行分析。基于2.2節(jié)中的模型，計(jì)算了干沉降通量的影響，地面高度下風(fēng)向核素濃度如圖2所示，干沉降通量如圖3所示。

圖2 地面高度下風(fēng)向放射性活度Fig.2 Radioactivity in wind direction at ground height

從圖2和圖3可以看出，地面核素濃度隨下風(fēng)向距離的增大而先增大后迅速減小，在下風(fēng)向約300 m達(dá)到最大值。干沉降通量比地面核素濃度低3～4個(gè)數(shù)量級(jí)，但隨著下風(fēng)向距離的增大，作用面積增大，沉降總量會(huì)隨著作用面積和作用時(shí)間的增大而增大。發(fā)生核泄漏事故后，計(jì)算海洋中放射性核素的遷移和擴(kuò)散，是否考慮大氣核素干沉降通量的影響，需要根據(jù)海洋核素濃度和大氣沉降作用時(shí)間確定。針對(duì)核泄漏事故應(yīng)急響應(yīng)的短期模擬，如果泄漏到海洋中的放射性核素很多，核素濃度很大，可以忽略大氣核素干沉降通量的影響；如果海洋中核素濃度較小，需要考慮核素干沉降通量的影響。長(zhǎng)期核素?cái)U(kuò)散模擬，需考慮大氣中核素干沉降影響。干沉降的核素適合以面源作為海洋核素?cái)U(kuò)散的源項(xiàng)。

圖3 干沉降通量Fig.3 Dry sedimentation flux

3.3 濕沉降影響

基于2.2節(jié)中的模型，計(jì)算了濕沉降對(duì)海洋中核素濃度的影響。根據(jù)1.2節(jié)中的大氣環(huán)境，設(shè)定沖刷比例常數(shù)Λ=1.6×10-4/s，則濕沉降通量如圖4所示。

圖4 濕沉降通量Fig.4 Wet sedimentation flux

從圖4可以看出，濕沉降通量隨著下風(fēng)向距離的增大迅速減小。300 m范圍內(nèi)濕沉降通量比干沉降通量大很多，300 m外濕沉降通量比干沉降通量大2個(gè)數(shù)量級(jí)。計(jì)算海洋中放射性核素的遷移和擴(kuò)散，是否考慮大氣核素濕沉降的影響，與干沉降判斷方法相同。小范圍降雨情況下，降雨區(qū)內(nèi)濕沉降通量適合以點(diǎn)源作為海洋核素?cái)U(kuò)散的源項(xiàng)。

3.4 懸浮物吸附沉降影響

3.4.1 分配系數(shù)影響

根據(jù)2.3節(jié)懸浮物對(duì)核素的吸附沉降模型和1.3節(jié)懸浮物的環(huán)境參數(shù)，不考慮核素的衰變等效應(yīng)，只考慮懸浮物作用。設(shè)定水深為44 m，懸浮物的沉積通量0.76 g/(cm2·a)，懸浮物的分配系數(shù)為10～103L/kg，計(jì)算了不同分配系數(shù)時(shí)，海洋中核素剩余放射性活度隨時(shí)間的變化，如圖5所示。

圖5 不同分配系數(shù)下放射性活度隨時(shí)間的變化Fig.5 Radioactivity change with time under different partition coefficient

從圖5可以看出，在理想的假定條件下，分配系數(shù)為10 L/kg時(shí)，核素的放射性活度基本沒(méi)有變化，懸浮物對(duì)核素吸附沉降作用可忽略。分配系數(shù)為100 L/kg時(shí)，懸浮物對(duì)核素吸附沉降作用開(kāi)始變得明顯，30 a后剩余60%左右。分配系數(shù)為500和103L/kg時(shí)，懸浮物的作用變得很大，30 a通過(guò)吸附沉降作用清除了大部分核素，剩余10%以內(nèi)。

3.4.2 沉積通量影響

設(shè)定水深為44 m，懸浮物的分配系數(shù)為500 L/kg，懸浮物的沉積通量為0～5 g/(cm2·a)，計(jì)算了在不同懸浮物沉積通量下，海洋中核素剩余放射性活度比率隨時(shí)間的變化，如圖6所示。

圖6 不同沉積通量下放射性活度隨時(shí)間的變化Fig.6 Radioactivity change with time under different sedimentation flux of suspended solids

從圖6可以看出，懸浮物沉積通量越大，懸浮物吸附沉降作用越明顯。當(dāng)沉積通量為0.01 g/(cm2·a)時(shí)，放射性活度在30 a內(nèi)基本沒(méi)有變化；當(dāng)沉積通量為0.1 g/(cm2·a)時(shí)，放射性活度開(kāi)始降低，30 a后降低為原來(lái)的70%左右；當(dāng)沉積通量為1～5 g/(cm2·a)時(shí)，放射性活度降低很快，懸浮物的影響程度較大，30 a后降低為原來(lái)的5%以下。

實(shí)際情況下懸浮物對(duì)核素有再懸浮和再解析作用，同時(shí)核素的分配系數(shù)和懸浮物的沉積通量也較小，懸浮物的作用會(huì)小很多。

3.5 生物影響

根據(jù)2.4節(jié)懸浮物對(duì)核素的富集模型和1.4節(jié)的參數(shù)，不考慮核素的衰變、懸浮物吸附沉降等效應(yīng)，只考慮生物富集作用，分別計(jì)算了生物量為10-5、10-7和10-9kg/L時(shí)，核素放射性活度隨分配系數(shù)的變化，如圖7所示。

圖7 生物吸收率隨富集系數(shù)變化Fig.7 Bioaccumulation rate with partition coefficient

從圖7可以看出，生物量一定時(shí)，剩余放射性活度的比率隨分配系數(shù)的增加而減小，生物的富集作用越明顯。當(dāng)海洋生物濃度為10-7kg/L和10-9kg/L時(shí)，海洋生物的富集作用非常小，可以忽略其影響。當(dāng)海洋生物濃度為10-5kg/L時(shí)，生物的富集作用開(kāi)始顯現(xiàn)，隨著分配系數(shù)的增加，放射性活度逐漸降低，當(dāng)分配系數(shù)達(dá)到104L/kg時(shí)，在生物富集作用下，放射性活度減小了10%左右。實(shí)際情況下，海洋中的生物量遠(yuǎn)低于10-5kg/L；平均分配系數(shù)也遠(yuǎn)小于104L/kg，海洋生物的富集作用對(duì)海水中放射性活度影響非常小，短期和長(zhǎng)期的核素?cái)U(kuò)散模擬均可以忽略。

3.6 各因素對(duì)比分析

綜合比較各因素的影響程度，研究了短期和長(zhǎng)期核素在海洋中的擴(kuò)散機(jī)制，給出了各因素的影響程度，遷移機(jī)制示意圖如圖8所示。

圖8 核素在海洋中的遷移機(jī)制Fig.8 Radionuclide migration mechanism in the ocean

從圖8可以看出，核泄漏事故開(kāi)始幾個(gè)月內(nèi)海洋中放射性核素濃度主要受源項(xiàng)、對(duì)流擴(kuò)散、大氣沉降和半衰期較短核素衰變的影響，影響程度依次減弱。幾十年后，在海洋環(huán)流作用下，核素在大洋中分布均勻，同時(shí)大氣中大部分核素通過(guò)干、濕沉降均已落入海洋，131I基本全部衰變，對(duì)流擴(kuò)散、大氣沉降和131I衰變影響可以忽略；而半衰期較長(zhǎng)的核素137Cs，衰變影響變大，同時(shí)懸浮物對(duì)放射性核素吸附沉降的影響程度越來(lái)越明顯，生物富集效應(yīng)也開(kāi)始有一定影響，但影響仍很小。

4 結(jié)論

1)核素在海洋中短期擴(kuò)散主要受源項(xiàng)、對(duì)流擴(kuò)散、大氣沉降和短半衰期核素衰變的影響，影響程度依次減弱；核素長(zhǎng)期擴(kuò)散主要受源項(xiàng)、長(zhǎng)半衰期核素衰變和懸浮物影響；海洋生物對(duì)核素濃度影響一直很小。

2)與已有工作相比，不但分析了核素在海洋中遷移的影響因素，而且建立模型量化研究了各因素的影響程度。

3)完善核素在海洋中遷移機(jī)制和遷移模型方面的理論。為核素?cái)U(kuò)散及核事故應(yīng)急處置工作提供了理論支持。

在此基礎(chǔ)上，還需進(jìn)一步研究核素在海洋中遷移的多因素耦合模型，建立核泄漏事故后果評(píng)價(jià)系統(tǒng)，為短期和長(zhǎng)期的事故處理提供技術(shù)支撐。