胡嘯峰，陳 鵬，曾昭龍

基于WRF的放射性物質(zhì)大氣擴(kuò)散模擬的環(huán)境輻射劑量估算方法

胡嘯峰，陳 鵬，曾昭龍

(中國(guó)人民公安大學(xué)警務(wù)信息工程學(xué)院，北京 100076)

針對(duì)核事故下放射性污染物質(zhì)輻射，提出了一種基于中尺度氣象模型 WRF 的放射性物質(zhì)大氣擴(kuò)散模擬的環(huán)境輻射劑量估算方法。該方法使用放射性物質(zhì)大氣擴(kuò)散數(shù)值模型對(duì)放射性污染物質(zhì)的空氣濃度和地表沉降濃度進(jìn)行計(jì)算，進(jìn)而結(jié)合有效劑量轉(zhuǎn)換因子得到環(huán)境輻射劑量，并基于輻射劑量標(biāo)準(zhǔn)劃分應(yīng)急處置區(qū)域，可以較好地刻畫(huà)瞬時(shí)排放及瞬時(shí)氣象條件導(dǎo)致的濃度變化造成的輻射劑量隨時(shí)間和空間的變化，從而實(shí)現(xiàn)對(duì)放射性污染物質(zhì)環(huán)境輻射劑量的實(shí)時(shí)計(jì)算。以2011年3月的日本福島核事故為案例，以放射性物質(zhì)131I 和137Cs為研究對(duì)象，使用該方法進(jìn)行了環(huán)境輻射劑量的估算。估算結(jié)果表明：131I 與137Cs 的空氣浸沒(méi)外照射有效劑量比較接近；成年人與嬰兒的吸入內(nèi)照射有效劑量分布也比較接近；地表沉積外照射有效劑量在總劑量中所占比例遠(yuǎn)遠(yuǎn)大于空氣浸沒(méi)外照射有效劑量和吸入內(nèi)照射有效劑量。

核事故；放射性物質(zhì)；環(huán)境輻射劑量；WRF 模型；數(shù)值模擬

核能是我國(guó)能源戰(zhàn)略中的重要組成部分，核能的安全利用是核電發(fā)展的前提，一旦發(fā)生事故，造成的后果對(duì)環(huán)境和人員將帶來(lái)巨大危害。自然災(zāi)害對(duì)核電的安全利用具有一定的威脅性。2011 年 3 月，距離日本本州島東海岸線約 120 km 的太平洋區(qū)域發(fā)生了里氏 9.0 級(jí)地震，間接導(dǎo)致了福島第一核電站中的反應(yīng)堆發(fā)生了核蒸汽泄漏，并于 3 月 14 日發(fā)生爆炸，大量的放射性物質(zhì)釋放至大氣中，通過(guò)大氣的擴(kuò)散作用輸送至周邊區(qū)域，并通過(guò)沉降作用對(duì)日本本土乃至周邊范圍的陸地和海洋造成了核污染，從而對(duì)環(huán)境安全和人員健康造成了極大的危害?？梢?jiàn)，自然災(zāi)害的次生效應(yīng)往往可以導(dǎo)致較之原生災(zāi)害更為嚴(yán)重的后果。因此，對(duì)核事故下自然環(huán)境的污染程度進(jìn)行及時(shí)、準(zhǔn)確的評(píng)估，已成為環(huán)境保護(hù)、核應(yīng)急等工作中的重要環(huán)節(jié)，而放射性污染物質(zhì)的輻射劑量是衡量其危害性的一個(gè)重要指標(biāo)。對(duì)此，本文提出一種基于中尺度氣象模型WRF的放射性物質(zhì)大氣擴(kuò)散模擬的環(huán)境輻射劑量估算方法。

1 研究方法

1.1 環(huán)境輻射劑量估算方法

放射性物質(zhì)的輻射途徑主要包括空氣浸沒(méi)外照射、地表沉積外照射、呼吸道吸入內(nèi)照射和消化道攝入內(nèi)照射等多種途徑，而對(duì)上述途徑的輻射劑量進(jìn)行評(píng)估則均需要以空氣、地表及水體等處的放射性物質(zhì)濃度作為依據(jù)，進(jìn)而通過(guò)劑量轉(zhuǎn)換因子得到環(huán)境輻射劑量值。美國(guó)的國(guó)家大氣排放資訊中心(NARAC)[1]開(kāi)展了一系列的輻射劑量與健康效應(yīng)的分析，其中，輻射劑量使用 Oak Ridge 國(guó)家實(shí)驗(yàn)室提供的劑量轉(zhuǎn)換因子數(shù)據(jù)庫(kù)，通過(guò)空氣和地面污染計(jì)算模型來(lái)計(jì)算，轉(zhuǎn)換因子數(shù)據(jù)則由美國(guó)環(huán)境署(EPA)[2]出版的相關(guān)報(bào)告提供。在以往的研究中，使用高斯模型計(jì)算空氣中的放射性物質(zhì)濃度，并結(jié)合劑量轉(zhuǎn)換因子估算輻射劑量的方法得到了廣泛的應(yīng)用。如李智等[3]使用高斯煙羽模型和食物鏈模型相結(jié)合，修改了美國(guó)橡樹(shù)嶺國(guó)家實(shí)驗(yàn)室編制的 AIRDOS-EPA 程序，進(jìn)行了基于實(shí)際觀測(cè)數(shù)據(jù)的實(shí)例研究，計(jì)算了大亞灣核電站常規(guī)排放條件下的放射性物質(zhì)對(duì)核電站周邊80 km范圍內(nèi)的公眾造成的輻射劑量；Jeong等[4]利用高斯擴(kuò)散模型模擬得出各類(lèi)不同的放射性物質(zhì)(如60Co、137Cs、241Am)釋放造成的環(huán)境輻射劑量，并利用EPA提供的由吸入放射性物質(zhì)引起的癌癥發(fā)病率和死亡率相關(guān)系數(shù)，計(jì)算得到了因環(huán)境放射性污染引起的致癌率和死亡率；王海洋等[5]為預(yù)測(cè)“臟彈” 恐怖襲擊的破壞效應(yīng)，建立了“臟彈”恐怖襲擊劑量評(píng)價(jià)模式，用 Gauss 煙團(tuán)模型對(duì)3種典型放射性物質(zhì)在大氣中的擴(kuò)散過(guò)程進(jìn)行計(jì)算，模擬給出有效劑量隨時(shí)間變化的分布；于京燕等[6]針對(duì)核事故環(huán)境應(yīng)急響應(yīng)的關(guān)鍵技術(shù)進(jìn)行了研究，并對(duì)其數(shù)值模擬計(jì)算中出現(xiàn)的難點(diǎn)進(jìn)行了分析。

通常情況下，空氣中的放射性物質(zhì)濃度CA(Bq/m3)可根據(jù)以下公式計(jì)算[7]：

(1)

式中：Pp為一年中由排放源吹向p點(diǎn)的風(fēng)向?qū)?yīng)的時(shí)間所占的比例；F為高斯擴(kuò)散系數(shù)(m2)；Qi為放射性物質(zhì)i的年平均排放率(Bq/s);ua為年平均風(fēng)速(m/s)。

而地表沉降的放射性物質(zhì)濃度Cgr(Bq/m2)則可由沉降率的計(jì)算得到[7]：

(2)

式中:di為放射性物質(zhì)的沉降率(Bq/m2·s)，即干沉降率與濕沉降率的總和;tb為源項(xiàng)排放的持續(xù)時(shí)間(s);λt為總的衰減率(s-1)，即放射性衰變率λ與土壤清除率λs之和。

式(1)對(duì)空氣中放射性物質(zhì)濃度的計(jì)算方法基于高斯擴(kuò)散系數(shù)，計(jì)算結(jié)果的準(zhǔn)確性有限，且各參數(shù)均采用年平均數(shù)值，不適用于事故狀態(tài)下短時(shí)間內(nèi)釋放大量物質(zhì)情況下的環(huán)境輻射劑量評(píng)估。而實(shí)際情況下，式(2)中的沉降率也不是一個(gè)常數(shù)，而是隨時(shí)間和空間變化的值。

從目前的研究進(jìn)展來(lái)看，事故狀態(tài)下釋放出的放射性物質(zhì)，其在空氣中的濃度的瞬時(shí)值對(duì)風(fēng)場(chǎng)、降水及源項(xiàng)均非常敏感[8-14]，則對(duì)濃度的瞬時(shí)值的估計(jì)準(zhǔn)確與否將在很大程度上取決于是否能夠合理估算出環(huán)境輻射劑量，因此利用能夠更好地處理復(fù)雜邊界條件、大氣湍流運(yùn)動(dòng)的大氣擴(kuò)散模型來(lái)進(jìn)行放射性物質(zhì)濃度及輻射劑量預(yù)測(cè)的方法具有一定的研究意義。為此，本文提出基于放射性物質(zhì)大氣擴(kuò)散模擬的方法來(lái)計(jì)算空氣中放射性物質(zhì)濃度，其優(yōu)點(diǎn)是：①可以刻畫(huà)瞬時(shí)排放及瞬時(shí)氣象條件導(dǎo)致的放射性物質(zhì)濃度變化；②)計(jì)算精度高，可以計(jì)算受復(fù)雜氣象條件和復(fù)雜地形條件影響下的放射性物質(zhì)濃度分布；③可以刻畫(huà)輻射劑量隨時(shí)間和空間的變化情況。由于基于放射性物質(zhì)大氣擴(kuò)散模擬估算出的環(huán)境輻射劑量不是一個(gè)固定值，而是隨時(shí)間、空間而變化的值，因此更符合實(shí)際應(yīng)用中對(duì)劑量估算的需求。在模型的選擇上，本文采用Hu等[15]基于WRF的放射性物質(zhì)大氣擴(kuò)散模型，該模型集成了用于計(jì)算放射源排放、放射性衰變過(guò)程、干沉降過(guò)程及濕沉降過(guò)程的模型，且氣象單元與物質(zhì)擴(kuò)散單元的計(jì)算完全同步，不需要由氣象場(chǎng)到擴(kuò)散模型的插值和驅(qū)動(dòng)，適合于對(duì)放射性物質(zhì)大氣擴(kuò)散的實(shí)時(shí)模擬，有利于得到實(shí)時(shí)的環(huán)境輻射劑量分布。

1.2 環(huán)境輻射劑量估算流程

基于放射性物質(zhì)大氣擴(kuò)散模擬的環(huán)境輻射有效劑量估算流程如圖1所示。

本文中對(duì)環(huán)境輻射有效劑量的計(jì)算基于放射性物質(zhì)大氣擴(kuò)散模型的輸出結(jié)果。按照?qǐng)D1的流程，首先由模型計(jì)算得到放射性物質(zhì)的空氣濃度和地表沉降濃度(本文中暫不涉及水體及食物中的濃度)，在此基礎(chǔ)上，計(jì)算出各種途徑以及各種核素的外照射有效劑量和內(nèi)照射有效劑量，從而得到總的有效劑量，最后根據(jù)估算的劑量值確定相應(yīng)的應(yīng)急措施。

圖1 基于放射性物質(zhì)大氣擴(kuò)散模擬的環(huán)境 輻射劑量估算流程Fig.1 Process of the assessment of radiation dose based on the atmospheric dispersion model of radionuclides

核事故中釋放出的放射性物質(zhì)經(jīng)過(guò)大氣的擴(kuò)散輸運(yùn)作用，分布于近地層的空氣中，對(duì)位于空氣與地表界面處的人員形成β或γ射線的外照射為空氣浸沒(méi)外照射?？諝饨](méi)外照射的年有效劑量可由下式計(jì)算得到[16]：

Eim=CA·DFim·Qf

(3)式中:Eim為空氣浸沒(méi)外照射的年有效劑量(Sv/a);CA為空氣中放射性物質(zhì)的濃度(Bq/m3);DFim為空氣浸沒(méi)外照射有效劑量轉(zhuǎn)換因子[(Sv·a-1)/(Bq·m-3)];Qf為考慮公眾在含有放射性物質(zhì)的空氣中暴露的時(shí)間份額及建筑物屏蔽作用導(dǎo)致的劑量降低因子。

本文中空氣中的放射性物質(zhì)濃度CA為由放射性物質(zhì)大氣擴(kuò)散模型計(jì)算得到的隨時(shí)間和空間變化的數(shù)值,在WRF/Chem模型中，該濃度由以下的擴(kuò)散方程描述：

(4)

式中:u為風(fēng)矢量;K為湍流擴(kuò)散系數(shù);Λs為濕沉降清除系數(shù)(s-1)；λ為衰變率(s-1)；Eemis為放射性物質(zhì)排放源項(xiàng)。

大氣中的放射性物質(zhì)經(jīng)過(guò)干沉降與濕沉降作用，沉積于地表，進(jìn)而對(duì)人員形成β或γ射線的外照射為地表沉積外照射。地表沉積外照射的年有效劑量可由下式計(jì)算[16]:

Egr=Cgr·DFgr·Qf

(5)

式中:Egr為地表沉積外照射的年有效劑量(Sv/a);Cgr為地表沉積的放射性物質(zhì)濃度(Bq/m2);DFgr為地表沉積外照射有效劑量轉(zhuǎn)換因子[(Sv·a-1)/(Bq·m-2)];Qf為考慮公眾在含有放射性物質(zhì)的空氣中暴露的時(shí)間份額及建筑物屏蔽作用導(dǎo)致的劑量降低因子。

本文在基于放射性物質(zhì)大氣擴(kuò)散模擬的估算中，地表沉積的放射性物質(zhì)的濃度Cgr可由下面公式計(jì)算得到：

(6)

(7)

空氣中的放射性氣體或氣溶膠粒子經(jīng)呼吸道進(jìn)入人體后，呼吸系統(tǒng)及體內(nèi)其他各組織和器官將會(huì)受到放射性物質(zhì)的內(nèi)照射，為吸入內(nèi)照射。吸入內(nèi)照射的年有效劑量可由下式計(jì)算[16]：

Einh=CA·Rinh·DFinh

(8)

式中:Einh為吸入內(nèi)照射的年有效劑量(Sv/a);CA為空氣中放射性物質(zhì)的濃度(Bq/m3)；DFinh為吸入內(nèi)照射有效劑量轉(zhuǎn)換因子(Sv/Bq);Rinh為人的呼吸率(m3/a)，不同的年齡組，具有不同的有效劑量轉(zhuǎn)換因子及呼吸率，根據(jù)國(guó)際原子能機(jī)構(gòu)(IAEA)的標(biāo)準(zhǔn)[7]，1歲嬰兒的呼吸率為1 400 m3/a，成年人的呼吸率為8 400 m3/a。

總有效劑量為所有核素的所有照射途徑的有效劑量之和，即

(9)

式中:m為核素的種類(lèi)數(shù)；n為照射途徑的總數(shù)。

2 算例分析

2.1 算例設(shè)置

本文選擇發(fā)生于 2011 年 3 月的日本福島核泄漏事故為案例，采用Hu等[15]改進(jìn)的WRF/Chem 模型模擬計(jì)算自UTC 2011/03/11/00至UTC 2011/03/31/00福島核電站所釋放出的放射性物質(zhì)的大氣擴(kuò)散濃度和地表沉降濃度。水平方向采用雙重網(wǎng)格嵌套區(qū)域(見(jiàn)圖2)，網(wǎng)格點(diǎn)數(shù)均為160×160，對(duì)應(yīng)的水平分辨率分別為9 km和3 km，網(wǎng)格中心點(diǎn)位于 37.5°N、141.0°E。模擬區(qū)域的垂直方向共分為 27 層，層頂氣壓為 100 hPa。模擬的初始與邊界條件由全球預(yù)報(bào)系統(tǒng)(Global Forecasting System，GFS)提供的水平分辨率為 0.5°×0.5°的分析數(shù)據(jù)提供，輸出結(jié)果的時(shí)間間隔為 6 h。

圖2 雙重網(wǎng)格嵌套區(qū)域示意圖Fig.2 Double-nested domains of the study area

物理參數(shù)設(shè)置方面，長(zhǎng)波輻射采用RRTM(Rapid Radiative Transfer Model)方案；短波輻射采用Dudhia方案；微物理過(guò)程采用WSM6(WRF single-moment 6-class)方案；積云對(duì)流參數(shù)化在d01中采用New Grell方案，在d02中不采用積云對(duì)流參數(shù)化方案；邊界層采用 MYJ (Mellor-Yamada-Janjic (Eta) TKE)方案；陸面采用Noah Land Surface方案。干沉降的參數(shù)化采用阻力模型[15](Resistance Model)，通過(guò)計(jì)算干沉降速度實(shí)現(xiàn)；濕沉降的參數(shù)化基于降水，通過(guò)計(jì)算清除系數(shù)[15](Scavenging Coefficient)實(shí)現(xiàn)。

本例中以放射性物質(zhì)131I和137Cs為研究對(duì)象，計(jì)算131I和137Cs的大氣擴(kuò)散濃度和地表沉降濃度分布，并據(jù)此計(jì)算其環(huán)境輻射有效劑量。輻射源的排放率采用 IAEA的估計(jì)結(jié)果，見(jiàn)圖3;131I和137Cs的衰變率、地表清除系數(shù)和有效劑量轉(zhuǎn)換因子的取值[7]詳見(jiàn)表1。

圖3 輻射源的排放率Fig.3 Emission rates of the radiative source

放射性物質(zhì)λλsDFimDFgrDFinh嬰DFinh成131I9.98×10-705.8×1071.2×1097.2×1085.4×109137Cs7.33×10101.62×1098.7×1071.8×1087.4×1094.6×109

注：λ為衰變率(s-1);λs為地表清除系數(shù)(s-1)；DFim為空氣浸沒(méi)外照射有效劑量轉(zhuǎn)換因子[(Sv·a-1)/(Bq·m-3)]；DFgr為地表沉積外照射有效劑量轉(zhuǎn)換因子[(Sv·a-1)/(Bq·m-2)]；DFinh嬰為嬰兒的吸入內(nèi)照射有效劑量轉(zhuǎn)換因子(Sv·Bq-1)；DFinh成為成年人的吸入內(nèi)照射有效劑量轉(zhuǎn)換因子(Sv·Bq-1)。

2.2 算例結(jié)果與分析

圖4為2011年3月31日0時(shí)131I和137Cs的空氣浸沒(méi)外照射年有效劑量分布圖。由圖4可見(jiàn)，131I與137Cs的空氣浸沒(méi)外照射年有效劑量分布比較接近，在2011年3月31日0時(shí)這一瞬時(shí)時(shí)刻，131I和137Cs的空氣浸沒(méi)外照射年有效劑量分布呈條帶狀，與空氣濃度的分布形狀相類(lèi)似，也反映了空氣濃度與空氣浸沒(méi)外照射劑量之間的關(guān)系。

圖5為2011年3月31日0時(shí)131I和137Cs的地表沉積外照射年有效劑量分布圖。由圖5可見(jiàn)，地表沉積外照射年有效劑量的分布范圍和中心區(qū)域的劑量值遠(yuǎn)遠(yuǎn)大于空氣浸沒(méi)外照射年有效劑量。其主要原因是:空氣浸沒(méi)外照射劑量的部分是由空氣中的放射性物質(zhì)造成的，而由于風(fēng)場(chǎng)的流動(dòng)性、放射性物質(zhì)的干濕沉降等作用，空氣中的放射性物質(zhì)隨時(shí)間的變化較大，且同一個(gè)空間位置的累積量較小，因此瞬時(shí)特征較為明顯；而地表沉積外照射是由沉降于地表的放射性物質(zhì)造成的，地表沉降放射性物質(zhì)的分布相對(duì)固定，而干濕沉降伴隨著擴(kuò)散輸運(yùn)的全過(guò)程，且土壤清除率較低，因此容易形成大面積的劑量分布。

圖4 2011年3月31日0時(shí)131I和137Cs的空氣浸沒(méi)外照射年有效劑量分布圖Fig.4 Distributions of the annual effective dose of external exposure from airborne radionuclides at 00 UTC on 31 March,2011

圖5 2011年3月31日0時(shí)131I和137Cs的地表沉積外照射年有效劑量分布圖Fig.5 Distributions of the annual effective dose of external exposure from ground deposited radionuclides at 00 UTC on 31 March,2011

圖6為2011年3月31日0時(shí)131I和137Cs的吸入內(nèi)照射年有效劑量分布圖。由圖6可見(jiàn)，在2011年3月31日0時(shí)這一瞬時(shí)時(shí)刻，131I和137Cs的吸入內(nèi)照射年有效劑量分布呈條帶狀，與空氣濃度的分布形狀相類(lèi)似，也反映了空氣濃度與吸入內(nèi)照射有效劑量之間的關(guān)系；此外，成年人與嬰兒的吸入內(nèi)照射年有效劑量分布也比較接近，這是因?yàn)閶雰旱暮粑蔬h(yuǎn)小于成年人的呼吸率，但是對(duì)于同樣放射性活度的放射性物質(zhì)，嬰兒所承受的劑量要遠(yuǎn)大于成年人，也就是嬰兒的有效劑量轉(zhuǎn)換因子較高，因此在呼吸率和劑量轉(zhuǎn)換因子的共同影響下，導(dǎo)致了成年人與嬰兒的吸入內(nèi)照射有效劑量的分布比較接近。

總的有效劑量為131I和137Cs的空氣浸沒(méi)外照射年有效劑量、地表沉積外照射年有效劑量及吸入內(nèi)照射年有效劑量的總和。圖7為2011年3月31日0時(shí)131I和137Cs的總有效劑量分布圖，對(duì)比圖4至圖6中的各種照射途徑的輻射劑量分布可以看出，地表沉積外照射有效劑量中總有效劑量的比例最大。

在實(shí)際的風(fēng)險(xiǎn)評(píng)估中，通常根據(jù)總有效劑量對(duì)暴露人群的健康風(fēng)險(xiǎn)進(jìn)行評(píng)估。例如普通公眾的年有效劑量限值為1 mSv/a，由圖7可見(jiàn)，超過(guò)該劑量限值區(qū)域的陸地面積較大，東西方向超過(guò)100 km，南北方向超過(guò)200 km；從事防災(zāi)減災(zāi)工作人員(如警察、消防員)的年有效劑量限值為50 mSv/a，從事應(yīng)急工作人員的年有效劑量限值為250 mSv/a，由圖7可見(jiàn)，這兩個(gè)限值對(duì)應(yīng)的陸地區(qū)域面積較小，半徑均在50 km之內(nèi)。

圖6 2011年3月31日0時(shí)131I和137Cs的吸入內(nèi)照射年有效劑量分布圖Fig.6 Distributions of the annual effective dose of internal exposure from inhalation at 00 UTC on 31 March,2011

圖7 2011年3月31日0時(shí)131I和137Cs的總有效 劑量分布圖Fig.7 Distributions of the annual total effective doses at 00 UTC on 31 March,2011

在核事故下，根據(jù)一定的應(yīng)急處置劑量標(biāo)準(zhǔn)，可以劃分相應(yīng)的應(yīng)急處置區(qū)域，每一類(lèi)區(qū)域?qū)?yīng)一種或幾種應(yīng)急處置策略。使用放射性物質(zhì)大氣擴(kuò)散模擬結(jié)果的情況下，應(yīng)急處置區(qū)域的劃分可以更為精確，而不是僅以疏散半徑進(jìn)行粗略劃分。此外，由于WRF/Chem模型的氣象場(chǎng)和擴(kuò)散場(chǎng)可以同步計(jì)算，因此可以實(shí)現(xiàn)對(duì)大氣擴(kuò)散和地表沉降濃度的實(shí)時(shí)計(jì)算，從而也可以實(shí)時(shí)計(jì)算應(yīng)急處置區(qū)域的變化范圍和變化趨勢(shì)。表2為國(guó)際原子能機(jī)構(gòu)(IAEA)和日本原子能規(guī)制委員會(huì)(NRA)分別提供的核事故下的應(yīng)急處置劑量標(biāo)準(zhǔn)。由表2可見(jiàn)，NRA的標(biāo)準(zhǔn)比IAEA的標(biāo)準(zhǔn)更為保守，也就是IAEA的每一條應(yīng)急措施對(duì)應(yīng)的劑量下限要高于NRA。

表2 核事故下的應(yīng)急處置劑量標(biāo)準(zhǔn)

圖8和圖9分別為IAEA和NRA標(biāo)準(zhǔn)下的應(yīng)急處置區(qū)域分布圖。由圖8和圖9可見(jiàn)，兩種標(biāo)準(zhǔn)下均不存在需要數(shù)小時(shí)內(nèi)緊急撤離的區(qū)域；NRA標(biāo)準(zhǔn)對(duì)應(yīng)的一周內(nèi)實(shí)施臨時(shí)疏散的區(qū)域和限制食用當(dāng)?shù)厣a(chǎn)食品的區(qū)域均小于IAEA標(biāo)準(zhǔn)對(duì)應(yīng)的區(qū)域。

圖8 IAEA標(biāo)準(zhǔn)下的應(yīng)急處置區(qū)域分布圖Fig.8 Distribution of the emergency disposal region under the standard from IAEA

圖9 NRA標(biāo)準(zhǔn)下的應(yīng)急處置區(qū)域分布圖Fig.9 Distribution of the emergency disposal area under the standard from NRA

表3為IAEA和NRA標(biāo)準(zhǔn)下三種應(yīng)急措施對(duì)應(yīng)的應(yīng)急處置區(qū)域面積和等效半徑。

表3 IAEA和NRA標(biāo)準(zhǔn)下三種應(yīng)急措施對(duì)應(yīng)的

3 結(jié)論與建議

(1) 本文提出了基于放射性物質(zhì)大氣擴(kuò)散模擬的環(huán)境輻射劑量估算方法，并結(jié)合福島核事故案例，以放射性物質(zhì)131I和137Cs為研究對(duì)象，對(duì)核電站周邊區(qū)域的131I和137Cs空氣浸沒(méi)外照射年有效劑量、地表沉積外照射年有效劑量、吸入內(nèi)照射年有效劑量和總有效劑量進(jìn)行了估算。估算結(jié)果表明：131I與137Cs的空氣浸沒(méi)外照射年有效劑量比較接近；成年人與嬰兒的吸入內(nèi)照射年有效劑量分布也比較接近；地表沉積外照射年有效劑量在總有效劑量中所占比例遠(yuǎn)遠(yuǎn)大于空氣浸沒(méi)外照射年有效劑量和吸入內(nèi)照射年有效劑量。可見(jiàn)，與一般大氣污染物不同的是，沉降作用并未減小放射性物質(zhì)污染對(duì)人體的危害，而是使放射性物質(zhì)大量沉積于地表，并通過(guò)外照射途徑作用于人體。

(2) 根據(jù)國(guó)際原子能機(jī)構(gòu)(IAEA)和日本原子能規(guī)制委員會(huì)(NRA)的劑量標(biāo)準(zhǔn)，本文提出了基于放射性物質(zhì)大氣擴(kuò)散模擬的應(yīng)急處置區(qū)域的劃分方法，福島案例中，NRA標(biāo)準(zhǔn)對(duì)應(yīng)的一周內(nèi)實(shí)施臨時(shí)疏散的區(qū)域和限制食用當(dāng)?shù)厣a(chǎn)食品的區(qū)域均小于IAEA標(biāo)準(zhǔn)對(duì)應(yīng)的區(qū)域。

(3) 基于放射性物質(zhì)大氣擴(kuò)散模擬的環(huán)境輻射劑量估算方法較之傳統(tǒng)的高斯模型精確度高、適用范圍廣，且可以得到環(huán)境輻射劑量的實(shí)時(shí)分布，能夠更好地服務(wù)于核應(yīng)急工作。在未來(lái)的研究中，可考慮結(jié)合水體及食物中的放射性物質(zhì)濃度，計(jì)算消化道攝入內(nèi)照射有效劑量。

致謝：本論文受公共安全協(xié)同創(chuàng)新中心支持。

[1] Bradley M M.NARAC:An emergency response resource for predicting the atmospheric dispersion and assessing the consequences of airborne radionuclides[J].JournalofEnvironmentalRadioactivity,2007,96(1/2/3):116-121.

[2] EPA.LimitingValuesofRadionuclideIntakeandAirConcentrationandDoseConversionFactorsforInhalation,Submersion,andIngestion[R].Washington,DC:U.S.Environmental Protection Agency,1988.

[3] 李智,李華.核電站氣載放射性流出物環(huán)境輻射劑量計(jì)算[J].環(huán)境科學(xué)與技術(shù),2010(6):119-123.

[4] Jeong H,Park M,Jeong H,et al.Radiological risk assessment caused by RDD terrorism in an urban area[J].AppliedRadiationandIsotopes,2013,79(5):1-4.

[5] 王海洋，曾志，劉立業(yè),等.“臟彈” 恐怖襲擊煙團(tuán)劑量計(jì)算方法[J].清華大學(xué)學(xué)報(bào) ( 自然科學(xué)版),2009,49(2):168-171.

[6] 于京燕，張大川，譚季青.核事故環(huán)境緊急響應(yīng)的風(fēng)險(xiǎn)研究[J].安全與環(huán)境工程,2014，21(4):149-153.

[7] IAEA.GenericModelsforUseinAssessingtheImpactofDischargesofRadioactiveSubstancestotheEnvironment[R].Vienna:IAEA,2001.

[8] Korsakissok I,Mathieu A,Didier D.Atmospheric dispersion and ground deposition induced by the Fukushima nuclear power plant accident:A local-scale simulation and sensitivity study[J].AtmosphericEnvironment,2013,70:267-279.

[9] Morino Y,Ohara T,Nishizawa M.Atmospheric behavior,deposition,and budget of radioactive materials from the Fukushima Daiichi nuclear power plant in March 2011[J].GeophysicalResearchLetters,2011,38(7):136-147.

[10]Morino Y,Ohara T,Watanabe M,et al.Episode analysis of deposition of radiocesium from the Fukushima Daiichi nuclear power plant accident[J].EnvironmentalScienceandTechnology,2013,47(5):2314-2322.

[11]Basit A,Espinosa F,Avila R,et al.Simulation of atmospheric dispersion of radionuclides using an Eulerian-Lagrangian modelling system[J].JournalofRadiologicalProtection,2008,28(4):539-561.

[12]Mathieu A,Korsakissok I,Quelo D,et al.Atmospheric dispersion and deposition of radionuclides from the Fukushima Daiichi nuclear power plant accident[J].Elements,2012,8(3):195-200.

[13]Ten Hoeve J E,Jacobson M Z.Worldwide health effects of the Fukushima Daiichi nuclear accident[J].Energy&EnvironmentalScience,2012,5(9):8743-8757.

[14] Takemura T,Nakamura H,Takigawa M,et al.A numerical simulation of global transport of atmospheric particles emitted from the Fukushima Daiichi nuclear power plant[J].Sola,2011,7:101-104.

[15]Hu X,Li D,Huang H,et al.Modeling and sensitivity analysis of transport and deposition of radionuclides from the Fukushima Daichi accident[J].AtmosphericChemistryandPhysics,2014,14:11065-11092.

[16]宋妙發(fā)，強(qiáng)亦忠.核環(huán)境學(xué)基礎(chǔ)[M].北京：原子能出版社，1999.

Method of Assessment of Environmental Radiation Dose Based on Simulation of Atmospheric Dispersion of Radionuclides by Using WRF

HU Xiaofeng,CHEN Peng,ZENG Zhaolong

(SchoolofPoliceInformationEngineering,People’sPublicSecurityUniversityofChina,Beijing100076,China)

Aiming at the assement of the environmental radiation dose in nuclear accidents,this paper proposes a method of assessing environmental radiation dose based on simulation of atmospheric dispersion of radionuclides by using the mesoscale meteorological model WRF.The method calculates the air concentrations and ground deposition concentrations of the radioactive contaminants by the atmospheric dispersion numerical model of radionuclides,and further obtains the radiation dose of the surrounding environment combined with the dose conversion coefficients,and divides the region for emergency disposal based on the dose limitation standard.The method can accurately describe the spatial and temporal distribution of the radiation dose with changing air concentration of radionuclides caused by the instantaneous emission and meteorological conditions in order to realize the real-time calculation of the environmental radiation dose of the radioactive contaminants.The paper conducts a case study of the Fukushima accident (happened in March 2011,Japan) to assess the environmental radiation of radionuclides131I and137Cs with the method.The results indicate that the effective dose of external exposure from airborne radionuclides131I and137Cs are close to each other;the effective dose distribution of internal exposure from inhalation by infants and adults are also similar;for131I and137Cs,the proportion of the effective dose to total dose of external exposure from ground deposition is much larger than that of both external exposure from airborne radionuclides and internal exposure from inhalation.Key words:nuclear accident;radionuclides;environmental radiation dose;WRF model;numerical simulation

1671-1556(2016)02-0135-08

2015-08-07

2016-01-08

公安部技術(shù)研究計(jì)劃項(xiàng)目(2015JSYJC 50)

胡嘯峰(1986—)，男，博士，講師，主要從事公共安全與應(yīng)急管理方面的研究。E-mail: huxiaofeng3@163.com

X946；R144

A

10.13578/j.cnki.issn.1671-1556.2016.02.025