核電廠氣載流出物大氣彌散因子不同計(jì)算方法比較研究

張 瓊，郭瑞萍，王 博

（環(huán)境保護(hù)部核與輻射安全中心，北京 100082）

?

核電廠氣載流出物大氣彌散因子不同計(jì)算方法比較研究

張 瓊，郭瑞萍，王 博*

（環(huán)境保護(hù)部核與輻射安全中心，北京 100082）

摘要：當(dāng)前，基于核電廠氣載放射性流出物環(huán)境影響的評價(jià)模式種類較多，然而這些模式及其開發(fā)的程序多用于評價(jià)過程，專門針對環(huán)境后果影響審評的應(yīng)用較少。本文從分析核電廠氣載放射性流出物大氣彌散的相關(guān)導(dǎo)則標(biāo)準(zhǔn)以及基于此開發(fā)的相應(yīng)程序入手，針對典型核素，采用C-Airdos，ROULEA-2. 0和Fluidyn-PANACHE等程序，以我國自主設(shè)計(jì)的CAP1400壓水堆重大專項(xiàng)示范工程石島灣核電廠所在地環(huán)境數(shù)據(jù)為基準(zhǔn)，開展正常運(yùn)行工況下大氣彌散因子模擬計(jì)算，與石島灣核電廠環(huán)境影響評價(jià)報(bào)告所采用的Y3001程序計(jì)算結(jié)果進(jìn)行比較分析，推薦適用于核電廠正常運(yùn)行工況下氣載流出物環(huán)境影響的審評模式及程序，為我國核電廠環(huán)境影響審評提供有益的技術(shù)參考。

關(guān)鍵詞：高斯模式；氣載流出物；大氣彌散因子；擴(kuò)散模式

隨著我國核電事業(yè)的快速發(fā)展，核電廠建設(shè)已成為核電中長期發(fā)展非常重要的規(guī)劃目標(biāo)之一。根據(jù)《中國核電中長期發(fā)展規(guī)劃》，到2015年核電規(guī)劃裝機(jī)容量達(dá)到4000萬千瓦，在建1800萬千瓦；2020年規(guī)劃裝機(jī)容量達(dá)到5800萬千瓦，在建3000萬千瓦［1，2］。環(huán)境保護(hù)部發(fā)布的《核動(dòng)力廠環(huán)境輻射防護(hù)規(guī)定》（GB 6249 -2011）［3］指出，在核電廠環(huán)境影響報(bào)告書中（選址、建造及運(yùn)行等階段）必須考慮核電廠氣載流出物對環(huán)境、生態(tài)和公眾的影響。在核電廠環(huán)境影響審評過程中，應(yīng)針對上述要求，對核電廠正常運(yùn)行和事故工況下氣載放射性核素對周圍環(huán)境以及公眾的影響進(jìn)行定量審評，同時(shí)，核電廠應(yīng)急計(jì)劃區(qū)劃分的基本依據(jù)也來自核素?cái)U(kuò)散后在不同距離段區(qū)域內(nèi)的濃度分布。因此，研究放射性核素在大氣中的擴(kuò)散特征成為環(huán)境影響審評的一個(gè)關(guān)鍵因素。

關(guān)于核電廠大氣彌散模式種類較多，如高斯模式，拉格朗日模式，歐拉模式，CFD模式等，基于此開發(fā)的程序也較多，如C-AIRDOS，ROULEA2，CEIRA，CALPUFF和Fluidyn-PANACHE等［4，5］。上述模式和程序被不同的設(shè)計(jì)單位用來開展核電廠《安全分析報(bào)告》和《環(huán)境影響報(bào)告》中大氣彌散因子及環(huán)境影響的分析與評價(jià)。審核計(jì)算的特點(diǎn)是既具有代表性，簡易可行，又不可過于保守。目前，還沒有專門針對審評計(jì)算的相關(guān)模式和軟件，本研究旨在通過比較分析各個(gè)模式和軟件的差異，在此基礎(chǔ)上探索推薦一套適用于常規(guī)核電廠氣態(tài)流出物環(huán)境影響評價(jià)的審評模式和程序，為我國核電審評提供有益的技術(shù)參考。

1 核電廠大氣擴(kuò)散模式

關(guān)于核電廠氣態(tài)流出物環(huán)境擴(kuò)散模式的選取，國際原子能機(jī)構(gòu)（International Atomic Energy Agency，簡稱IAEA）與美國核管會(huì)（Nuclear Regulatory Commission，簡稱NRC）先后頒布了相關(guān)導(dǎo)則，主要有：國際原子能機(jī)構(gòu)于1980年頒布的安全指南（A Safety Guide，No. 50 - SG -G3）“核電廠選址有關(guān)的大氣彌散”［5］，以及于1986年頒布的技術(shù)文件（TECDOC -379）“應(yīng)用于放射性核素釋放的大氣彌散模式”［6］；美國核管會(huì)曾就核電廠正常及事故工況下的大氣擴(kuò)散發(fā)布過相關(guān)導(dǎo)則，主要包括導(dǎo)則RG1. 111《輕水堆正常工況下排放的氣態(tài)排出物的大氣輸運(yùn)與彌散的計(jì)算方法》（USNRC Regulatory Guide1. 111）［7］和導(dǎo)則RG1. 145《用于核電廠潛在事故后果評價(jià)的大氣彌散模式》（USNRC Regulatory Guide 1. 145，1979）等［8］。其中，導(dǎo)則RG1. 111給出了用于評價(jià)陸上輕水反應(yīng)堆正常運(yùn)行氣載流出物大氣輸送和彌散的計(jì)算模式及其假設(shè)，XOQDOQ是實(shí)現(xiàn)RG1. 111的計(jì)算程序，該程序的基本原理為高斯直線軌跡煙羽模式，模式中考慮了干沉降和放射性衰變的煙羽損耗，以及對非直線軌跡的修正等，但不考慮濕沉降影響。XOQDOQ采用的數(shù)據(jù)為風(fēng)向、風(fēng)速、穩(wěn)定度的三維聯(lián)合頻率，對于平坦地形可單獨(dú)使用XOQDOQ，對于復(fù)雜地形常需聯(lián)合使用MESODIF（分段煙羽軌跡模式）。

關(guān)于我國核電廠擴(kuò)散模式的選取，《核電廠廠址選擇的大氣彌散問題》（HAD 101/102）等有關(guān)的法規(guī)和導(dǎo)則對核電廠釋放的放射性物質(zhì)大氣彌散所涉及的氣象特征、估算模式及有關(guān)參數(shù)等方面做出了具體要求［9］，推薦采用雙高斯方法進(jìn)行廠址大氣彌散計(jì)算，源的類型被劃分為地面點(diǎn)源和高架點(diǎn)源，HAD 101/102同時(shí)給出了放射性核素衰變修正、干沉積、濕沉積、煙羽抬升高度、靜風(fēng)修正計(jì)算方法等。

目前，在我國核電廠環(huán)境影響評價(jià)過程中，有關(guān)正常運(yùn)行下的大氣擴(kuò)散因子估算是以美國核管會(huì)管理導(dǎo)則RG1. 111推薦的方法為基礎(chǔ)。表1為我國環(huán)評單位通常采用的部分大氣擴(kuò)散模型和輻射環(huán)境影響評價(jià)程序。這些評價(jià)程序的主要差異表現(xiàn)在氣象資料的預(yù)處理方法和應(yīng)用形式以及擴(kuò)散模型參數(shù)的選取等方面。我國核電廠環(huán)境影響評價(jià)和審評過程中常用的模型基本為高斯煙羽模型。對長期大氣彌散因子計(jì)算主要為Y3001/ C-AIRDOS、ROULEA和CEIRA等，短期的主要為PAVAN，內(nèi)陸核電廠考慮到地形的復(fù)雜性，環(huán)評報(bào)告中用CFD進(jìn)行模擬計(jì)算。

表1 核電廠正常運(yùn)行氣載放射性污染物釋放后果評價(jià)程序Table 1 Assessment procedures on atmospheric radioactive effluxents under normal operation of nuclear power plants in China

2 CAP1400石島灣核電廠環(huán)境影響評價(jià)分析

核電廠址按其地理位置可分為濱海廠址和內(nèi)陸廠址，目前我國在建和運(yùn)行的均為濱海廠址，因此選取CAP1400石島灣濱海核電廠址為例進(jìn)行分析具有代表性［10，11］。CAP1400為我國自行研發(fā)的裝機(jī)容量為140萬千瓦的先進(jìn)非能動(dòng)核電機(jī)組，屬于第3代核電技術(shù)。石島灣核電廠址瀕臨黃海，地貌類型為濱海緩坡地，陸地大部分自然地面高程在0 m -30 m之間，近海地形較為平緩。該核電廠址具有比較完整的氣象監(jiān)測數(shù)據(jù)和多年的環(huán)境監(jiān)測數(shù)據(jù)，選用2009年-2010年的氣象數(shù)據(jù)進(jìn)行比對計(jì)算。石島灣廠址屬暖溫帶大陸性季風(fēng)氣候，受海洋性氣候影響明顯，四季分明、夏季雨量充沛、秋冬雨量較少、日照充足、冬寒夏熱，溫差較大。此區(qū)域中緯度天氣系統(tǒng)活動(dòng)頻繁，大風(fēng)、暴雨較多；同時(shí)，由于海洋性氣候影響明顯，還具有大陸性氣候向海洋性氣候過渡型的氣候特征。同時(shí)，核電廠的氣載流出物主要通過煙囪排放到大氣，煙囪高度為81 m，內(nèi)徑3 m，安全殼上水箱高度暫定為76 m。

正常工況下核電廠環(huán)境影響評價(jià)的氣態(tài)流出物根據(jù)氣象參數(shù)模擬計(jì)算出大氣彌散因子，結(jié)合食譜和人口數(shù)據(jù)，計(jì)算出個(gè)人和集體劑量；液態(tài)流出物則根據(jù)稀釋因子進(jìn)行計(jì)算。本文選用的C -AIRDOS可以計(jì)算出大氣彌散因子和氣態(tài)流出物導(dǎo)致的個(gè)人和集體劑量；ROULEA -2. 0可以計(jì)算氣態(tài)流出物和液態(tài)流出物的個(gè)人和集體劑量；Fluidyn -PANACHE程序主要用于計(jì)算大氣彌散因子，結(jié)合自行編制的劑量模式進(jìn)行劑量計(jì)算。本文重點(diǎn)針對石島灣CAP1400兩臺(tái)機(jī)組正常運(yùn)行氣態(tài)流出物大氣彌散特征進(jìn)行計(jì)算，并選取核素131I，與環(huán)評報(bào)告采用Y3001程序的計(jì)算結(jié)果進(jìn)行比對分析，在此基礎(chǔ)上推薦一套適用于常規(guī)核電廠氣態(tài)流出物環(huán)境影響審評的模式和程序，為我國的核電審評提供技術(shù)支持。

2. 1 C -AIRDOS模擬計(jì)算

C - AIRDOS采用美國RG1. 111《估計(jì)輕水堆正常釋放的氣態(tài)物質(zhì)的大氣遷移和擴(kuò)散方法》以及我國核安全導(dǎo)則《核電廠廠址選擇的大氣彌散問題》推薦的模式方法，考慮了混合釋放、煙云抬升，以及干濕沉積損耗，核素衰變損耗等因素的修正。對于核素連續(xù)釋放所致下風(fēng)向某點(diǎn)大氣彌散因子采用高斯煙羽直線軌跡模式，公式如下［12］。

式中：

C（x，y，z）為下風(fēng)向某點(diǎn)（x，y，z）處的空氣中污染物濃度（Bq·m-3）；

x為下風(fēng)向距離（m）；y為橫截風(fēng)向距離（m）；z為距地面高度（m）；

Q為氣載核素源強(qiáng)，即釋放率（Bq·S-1）；

σy，σz分別為水平方向和垂直方向擴(kuò)散參數(shù)（m），它們是下風(fēng)向距離x及大氣穩(wěn)定度的函數(shù)；

He為有效排放高度（m）；

u-為有效源高度處的平均風(fēng)速（m·S-1）。

通常將擴(kuò)散因子定義為釋放每單位氣載放射性核素在下風(fēng)向某處造成的氣載放射性核素濃度，其單位為s·m-3。對于年均濃度，考慮長期擴(kuò)散的風(fēng)擺效應(yīng)，通常取22. 5度扇形方位的平均值，每個(gè)扇區(qū)內(nèi)平均地面空氣濃度為：

以CAP1400石島灣廠址主導(dǎo)風(fēng)向SSW和次主導(dǎo)風(fēng)向NW的下風(fēng)向NNE及SE方位為例進(jìn)行大氣彌散因子的對比分析。模擬計(jì)算過程中相關(guān)輸入?yún)?shù)參考《國核壓水堆示范工程環(huán)境影響報(bào)告書》［11］。圖1為分別利用C-AIRDOS和Y3001計(jì)算131I大氣彌散因子的結(jié)果比較圖。由圖可知，環(huán)評報(bào)告（Y3001）（以后綴1表示，下同）的計(jì)算結(jié)果和C-AIRDOS（以后綴2表示，下同）相近，在兩個(gè)主導(dǎo)風(fēng)向下風(fēng)向的大氣彌散因子均是隨著距離的增加而減小，且減小趨勢基本保持一致。圖2為二者的相對誤差圖。相對誤差在釋放點(diǎn)附近相差最大，約為40%，到1 km處下降為0，到2 km處為最小，約為-50%，隨后到廠址邊界5 km處變化不大。之后隨著距離的增大，相對誤差緩慢增大到0，到30 km -40 km處又變?yōu)檎怠?/p>

相對誤差計(jì)算公式為：

相對誤差%=（SE2 -SE1）*100/ SE1

或：相對誤差%=（NNE2 - NNE1）*100/ NNE1。其他核素與此類似。

圖1131I大氣彌散因子在NNE和SE風(fēng)向比較分析圖（C-AIRDOS&Y3001）Fig. 1 Comparative analysis of dispersion factor of131I in NNE and SE directions between C-AIRDOS and Y3001

圖2131I大氣彌散因子在NNE和SE風(fēng)向相對誤差示例圖（C-AIRDOS&Y3001）Fig. 2 Relative error comparison of atmospheric dispersion factor of131I in the NNE and SE directions between C-AIRDOS and Y3001

2. 2 ROULEA模擬計(jì)算

ROULEA -2. 0同樣是依據(jù)美國RG1. 111及HAD101/02中推薦的模式方法，考慮了混合釋放、煙云抬升、干濕沉積損耗，以及核素衰變損耗等因素的修正。地面排放的濃度計(jì)算公式為［13］：

x為下風(fēng)向距離，m；

i為第i風(fēng)速級；

j為第j大氣穩(wěn)定度；

Ui為第i風(fēng)速級的中值；

σzj（x）為距離x處穩(wěn)定度j的垂直煙羽擴(kuò)散系數(shù)，m；

fij（K）為第i風(fēng)速級，第j穩(wěn)定度和第K方位的聯(lián)合頻率；

DECi（x）為x距離第i風(fēng)速級由于放射性衰變的損耗因子；

DEPLij（x，K）為x距離i風(fēng)速級，第j穩(wěn)定度和第K方位由于煙羽損耗的因子；

RF（x，K）為x距離和第K方位的環(huán)流和滯流的修正因子；

Dz為建筑物高度，m。

以上兩個(gè)公式表示了由于建筑物尾流的最大附加擴(kuò)散。比較這兩個(gè)公式的計(jì)算結(jié)果，并取保守的較大的X/ Q值。

圖3為利用ROULEA - 2. 0和Y3001計(jì)算131I大氣彌散因子的比較結(jié)果。Y3001的計(jì)算結(jié)果和ROULEA -2. 0的趨勢相似，在兩個(gè)主導(dǎo)風(fēng)向下風(fēng)向的大氣彌散因子均是隨著距離的增加而減小，且減小趨勢基本保持一致。圖4為二者的相對誤差圖。在釋放點(diǎn)附近為正，約為45%，隨著距離的增加逐漸減小，到1km附近為0，到2km處下降減緩，約為-25%，隨后隨著距離增大到10km處均緩慢下降，10km處為拐點(diǎn)，相對誤差約為-10%，之后隨著距離的增大逐漸下降，最大值則在80km處，超過了-75%。

2. 3 Fluidyn-PANACHE模擬計(jì)算

Fluidyn-PANACHE是一個(gè)用于模擬與污染和危險(xiǎn)相關(guān)大氣過程的三維計(jì)算流體力學(xué)（Computational Fluid Dynamics，簡稱CFD）診斷計(jì)算模型。它使用三維有限體積方法求解模擬空氣運(yùn)動(dòng)的Navier - Stokes方程（動(dòng)量方程），同時(shí)求解每個(gè)網(wǎng)格中物種的濃度、質(zhì)量和能量守恒方程。對于湍流模式，PANACHE求解上述方程的雷諾平均形式。

圖3131I大氣彌散因子在NNE和SE風(fēng)向比較分析圖（ROULEA -2. 0&Y3001）Fig. 3 Comparative analysis of dispersion factor of131I NNE and SE direction between ROULEA -2. 0 and Y3001

圖4131I大氣彌散因子在NNE和SE風(fēng)向相對誤差示例圖（ROULEA -2. 0&Y3001）Fig. 4 Relative error comparison of atmospheric dispersion factor of131I in the NNE and SE directions between ROULEA -2. 0 and Y3001

本案例的地形信息可直接導(dǎo)入軟件并形成等高線文件，可在其上手動(dòng)構(gòu)建廠址構(gòu)筑物、釋放源及監(jiān)測點(diǎn)等；由于本次模擬主要計(jì)算分析廠址構(gòu)筑物和地形對氣載流出物擴(kuò)散沉降的影響，故以2#煙囪單獨(dú)排放為例，選取釋放源周邊15 km開展三維計(jì)算流體力學(xué)模擬；本案例采用4層嵌套且逐層加密方式，總共劃分了46萬個(gè)結(jié)構(gòu)化網(wǎng)格，可平滑模擬周圍丘陵的地形起伏，經(jīng)驗(yàn)證，排除了網(wǎng)格依賴性，計(jì)算網(wǎng)格滿足HJ2. 2 -2008《環(huán)境影響評價(jià)技術(shù)導(dǎo)則大氣環(huán)境》的相關(guān)要求。其他有關(guān)三維計(jì)算流體力學(xué)求解器的設(shè)置，包括湍流方案、數(shù)值方案及松弛因子等具體設(shè)置請參見用戶手冊［14］。一旦上述設(shè)置完成即可開始初始化，包括處理所有的輸入，生成初始風(fēng)場等。

為模擬氣載流出物的年均遷移擴(kuò)散特征，基于廠址風(fēng)向、風(fēng)速和穩(wěn)定度的三維聯(lián)合頻率分布，提取發(fā)生頻率最高的前70類氣象條件（占總氣象條件85%以上）為典型氣象條件，分別計(jì)算該70類氣象條件下的穩(wěn)定風(fēng)場和濃度分布；通過PANACHE的Import Average功能將擴(kuò)散結(jié)果按照相應(yīng)概率計(jì)算其平均值；通過讀入平均后的最終結(jié)果，導(dǎo)出地面污染物濃度數(shù)據(jù)，并最終計(jì)算年均大氣彌散因子。

以典型核素131I為例計(jì)算得到該廠址半徑7. 5km區(qū)域范圍內(nèi)年均擴(kuò)散和沉積狀況，如圖5所示，廠址構(gòu)筑物對2號機(jī)組煙囪釋放氣載流出物的擴(kuò)散及富集效果，結(jié)合年風(fēng)玫瑰圖可知，廠址區(qū)域的主導(dǎo)風(fēng)向?yàn)镾SW，坐落于反應(yīng)堆廠房NNE方位的三廢處理廠房恰好位于主導(dǎo)風(fēng)向的下風(fēng)向，之所以在圖中顯示為正上方，是由于此前為了建模已預(yù)先將地形底圖向逆時(shí)針方向旋轉(zhuǎn)34度，因此，受三廢處理廠房影響，氣載流出物不僅易于在三廢廠房面向反應(yīng)堆廠房的區(qū)域形成滯留和累積，最大可達(dá)1. 13E -05 s·m- 3，還會(huì)在三廢廠房下風(fēng)向尾流區(qū)形成相對高濃度區(qū)域（8. 66E -06 s·m- 3），但對廠址其他區(qū)域，特別是西南側(cè)家屬生活區(qū)的影響可忽略不計(jì)。為便于直觀比較，將整個(gè)廠址區(qū)域16個(gè)方位及不同距離段的大氣彌散因子匯總于表2。

圖5 石島灣正常運(yùn)行煙囪釋放年均大氣擴(kuò)散特征三維計(jì)算流體力學(xué)模擬圖（ROULEA -2. 0&Y3001）Fig. 5 CFD simulation of the annual atmospheric dispersion characteristics of the chimney of the normal operation of the Shi Dao Bay

表2 廠址周邊7. 5 km范圍內(nèi)131I年均大氣彌散因子Table 2 The dispersion factor of131I around Shi Dao Bay around in 7. 5 km單位：s·m- 3

表2中各距離段各方位年均大氣彌散因子最大值以粗體特別標(biāo)注，從中可以看出，在釋放點(diǎn)500 m范圍內(nèi)，受下風(fēng)向構(gòu)筑物影響，年均大氣彌散因子最大值會(huì)略向東偏轉(zhuǎn)，隨著距離增加，地表構(gòu)筑物影響逐步減弱，在800 m后，年均大氣彌散因子最大值保持在主導(dǎo)風(fēng)向下風(fēng)向即NNE風(fēng)向，且與高斯模式的計(jì)算結(jié)果在同一數(shù)量級，相差較小，這表明三維計(jì)算流體力學(xué)數(shù)值模擬軟件Fluidyn -PANACHE可更加細(xì)致地模擬廠區(qū)構(gòu)筑物對大氣擴(kuò)散的影響，結(jié)果較準(zhǔn)確可靠。

圖6為核素131I的PANACH和Y3001在兩個(gè)主導(dǎo)風(fēng)向下風(fēng)向的大氣彌散因子相對誤差圖。由圖可見，相對誤差在釋放點(diǎn)處為最大，在NNE方位幾乎達(dá)到500%，在SE方位也接近200%。隨著距離的增加，誤差急劇減小，到1 km -2 km后變化趨于平緩。

圖6131I大氣彌散因子在NNE和SE風(fēng)向相對誤差示例圖（PANACHE& Y3001）Fig. 6 Relative error comparison of atmospheric dispersion factor of131I in the NNE and SE directions between PANACHE and Y3001

3 討論與結(jié)論

由計(jì)算結(jié)果可知，4個(gè)程序計(jì)算的大氣彌散因子結(jié)果趨勢基本一致，在兩個(gè)主導(dǎo)風(fēng)向下風(fēng)向的大氣彌散因子均是隨著距離的增加而迅速減小，且減小趨勢基本保持一致。ROULEA -2. 0的計(jì)算結(jié)果彌散因子最大，C-AIRDOS次之，表明ROULEA - 2. 0的計(jì)算結(jié)果過于保守，應(yīng)用于環(huán)評審評時(shí)需加以注意；Fluidyn-PANACHE在計(jì)算過程中耗時(shí)較多（歷時(shí)兩個(gè)月），計(jì)算范圍有限（半徑7. 5 km），說明用于正常工況下大氣彌散因子的計(jì)算有一定的局限性。CAIRDOS的彌散因子計(jì)算結(jié)果較Y3001稍大，由于其是Y3001的改進(jìn)版，對相關(guān)參數(shù)進(jìn)行了修正，因此，綜合分析比較已有的4種分析軟件，C-AIRDOS更適用于核電廠氣載流出物環(huán)境影響的審評。

對于核電廠正常運(yùn)行工況下，基于RG1. 111核心體系構(gòu)建了多個(gè)程序，由于在煙羽抬升、干濕沉降等方面修正的差異，以及氣象參數(shù)輸入方式的不同，如實(shí)時(shí)氣象參數(shù)和聯(lián)合頻率等，導(dǎo)致了計(jì)算結(jié)果有一定的差別。對采用Y3001、ROULEA -2. 0和C-AIRDOS所得大氣彌散因子對比分析表明，即使是同一類程序，由于修正模式的不同，結(jié)果也會(huì)出現(xiàn)較大差異。這部分的工作還需要深入開展，如進(jìn)行參數(shù)敏感性分析，可靠性分析等。

參考文獻(xiàn)

［1］國務(wù)院. R大氣污染防治行動(dòng)計(jì)劃［R］.北京：國務(wù)院，2013.

［2］潘自強(qiáng).堅(jiān)定地繼續(xù)發(fā)展核電是解決我國能源可持續(xù)發(fā)展的重要途徑［J］.科技導(dǎo)報(bào)，2012（31）：1 -3.

［3］環(huán)境保護(hù)部.核動(dòng)力廠環(huán)境輻射防護(hù)規(guī)定［S］.北京：環(huán)境保護(hù)部，2011.

［4］Zhang Q，Guo RP，Wang B. Radioactive airborne effluents and the environmental impact assessment of CAP1400 nuclear power plant under normal operation［J］. Nuclear Engineering and Design，2014（280）：579 -585.

［5］朱好，熊章輝，李鳳菊，等.典型內(nèi)陸廠址近場年均大氣彌散和沉積因子不同計(jì)算方法的比較［J］.輻射防護(hù)，2015（35）：83 -89.

［6］IAEA. Atmospheric Dispersion in Relation to Nuclear Power Plant Siting：A Safety Guide［R］. Vienna：IAEA，1980.

［7］IAEA. Atmospheric Dispersion Models for Application in Relation to Radionuclide Releases（TECDOC -379）［R］. Vienna：IAEA，1986.

［8］NRC，Regulatory Guide，1. 111. Methods for Estimating Atmospheric Transport and Dispersion of Gaseous Effluents in Routine Releases form Light-Water-Cooled Reactors［R］. Washington DC：NRC，1977.

［9］NRC，Regulatory Guide 1. 145. Atmospheric dispersion models for potential accident consequence assessments at nuclear power plants［R］. Washington DC：NRC，1983.

［10］國家核安全局.核安全導(dǎo)則（上冊）［M］.北京：中國法制出版社，1998.

［11］北京大學(xué)科技開發(fā)部.華能山東石島灣核電廠高溫氣冷堆核電示范工程大氣擴(kuò)散試驗(yàn)研究報(bào)告［R］.北京：北京大學(xué)，2008.

［12］國核示范電站有限公司.國核壓水堆示范工程環(huán)境影響報(bào)告書（選址階段）［R］.上海：上海核工程研究院，2013.

［13］陳曉秋. C-AIRDOS的模式與參數(shù)［M］.北京：核與輻射安全中心，2003.

［14］中國輻射防護(hù)研究院.核設(shè)施正常工況環(huán)境輻射質(zhì)量評價(jià)軟件說明書［R］.太原：中國輻射防護(hù)研究院，2014.

［15］Fluidyn. Fluidyn-PANACHE User Manual［R］. Paris：Fluidyn，2010.

The Comparative Study on Airborne Effluent Atmospheric Dispersion Factor of Nuclear Power Plant

ZHANG Qiong，GUO Qiong，WANG Bo*
（Nuclear and Radiation Safety Center，MEP，Beijing 100082，China）

Abstract：There are several kinds of models for the environmental impact assessment of airborne radioactive effluent from nuclear power plant. However，most of them and relevant softwares are used for the impact assessment；few of them are specifically developed for the review purpose. This paper analyzes the standards，guidelines and corresponding program related to the airborne radioactive effluent from NPP. Then the programs，such as ROULEA -2. 0，C-AIRDOS，F(xiàn)luidyn-PANACHE，were applied for numerical simulation of airborne effluent atmospheric dispersion factor for typical nuclide such as131I with the local environmental parameters of the Shidaowan NPP（CAP1400 PWR major demonstrative project）innormal operation. And the result was compared with the calculation result by the program Y3001 in the environmental impact report of Shidaowan NPP. Finally a set of suitable models and programs were recommended for reviewing the environmental impact of the airborne radioactive effluents of NPP in normal operation which could be the reference for the review on environmental impact of NPP in China.

Key words：Gaussian mode；airborne effluent；atmospheric dispersion factor；dispersion model

中圖分類號：X591

文章標(biāo)志碼：A

文章編號：1672-5360（2016）02-0052-06

收稿日期：2016-03-21 修回日期：2016-05-04

基金項(xiàng)目：大型先進(jìn)壓水堆及高溫氣冷堆核電站國家科技重大專項(xiàng)《CAP1400安全審評關(guān)鍵技術(shù)研究》，項(xiàng)目編號 2013ZX06002001

作者簡介：張瓊（1975—），女，寧夏固原人，工學(xué)博士，高級工程師，現(xiàn)主要從事輻射防護(hù)、大氣環(huán)境影響評價(jià)等相關(guān)工作

*通訊作者：王 博，E-mail：wangbo@ chinansc. cn