喬清黨，李 冰，岳會(huì)國，姚仁太，盛 黎，王韶偉，*

（1．環(huán)境保護(hù)部核與輻射安全中心，北京 100082；2．中國輻射防護(hù)研究院，太原 030006；3．國家氣象中心，北京 100081）

2011 年3月的日本福島核電廠事故［1］、2013年2月的朝鮮第三次核試驗(yàn)［2］，使得人們對(duì)于核事故、核試驗(yàn)情況下氣載放射性物質(zhì)對(duì)環(huán)境和公眾的輻射后果的關(guān)注度越來越高。我國及周邊國家和地區(qū)的許多核設(shè)施一旦發(fā)生核事故，則會(huì)對(duì)我國的環(huán)境和公眾造成輻射危害。表1給出了我國（不含我國臺(tái)灣地區(qū)）及周邊國家的主要核電機(jī)組數(shù)量（截止到2014年1月）。

核應(yīng)急是核安全縱深防御體系的最后一道屏障，是有效處置核與輻射突發(fā)事件，減緩事件后果，保障公共安全的重要舉措。核與輻射事故后果評(píng)價(jià)系統(tǒng)是核應(yīng)急響應(yīng)過程中決策的重要技術(shù)支持手段，要不斷增加事故后果評(píng)價(jià)手段、強(qiáng)化事故分析手段［3］。放射性污染實(shí)時(shí)預(yù)測(cè)與評(píng)價(jià)因其對(duì)突發(fā)事件的快速響應(yīng)而得到高度重視，已經(jīng)成為各國應(yīng)對(duì)核事故及核恐怖事件的重要技術(shù)手段之一。美國的NARAC、歐共體的RODOS都是多尺度范圍放射性后果評(píng)價(jià)與決策技術(shù)支持軟件系統(tǒng)，均具有精確預(yù)測(cè)本國及其周邊國家及地區(qū)不同位置核與輻射事故的中小尺度后果評(píng)價(jià)功能。中小尺度后果評(píng)價(jià)軟件在核應(yīng)急中的作用如圖1所示。

表1 我國及周邊國家主要核電機(jī)組數(shù)量Table 1 Numbers of nuclear power units in China and neighboring countries

圖1 中小尺度后果評(píng)價(jià)軟件在核應(yīng)急中的作用Fig.1 Role of mesoscale consequence assessment software in nuclear emergency

通過十多年的研究與開發(fā)，目前，我國已經(jīng)開發(fā)了多個(gè)用于核事故應(yīng)急決策的放射性后果評(píng)價(jià)系統(tǒng)［4］，但發(fā)展速度和技術(shù)水平與國外相比還存在一定的差距。李宗明在《從日本福島核事故審視核安全的政府、法律和監(jiān)管框架》一文中曾指出，“就核安全監(jiān)管的技術(shù)力量而言，在核安全理念的深入了解、專業(yè)覆蓋面、分析和計(jì)算程序、質(zhì)量保證、事故診斷和處理的判定、應(yīng)急響應(yīng)能力等方向需要提高的空間很多”［5］。在精細(xì)預(yù)報(bào)我國和周邊國家（地區(qū)）不同位置的核與輻射事故氣載釋放造成的后果影響，特別是核與輻射事故后果的回顧性評(píng)價(jià)研究方面，目前尚缺乏相應(yīng)的軟件系統(tǒng)，難以快速、準(zhǔn)確地做出精細(xì)化預(yù)報(bào)和評(píng)價(jià)。

本文擬建的具有預(yù)報(bào)和回顧性評(píng)價(jià)功能的精細(xì)化核事故氣載放射性物質(zhì)擴(kuò)散的數(shù)值模擬軟件的研發(fā)方法是對(duì)我國現(xiàn)有核事故后果評(píng)價(jià)軟件的補(bǔ)充，可實(shí)現(xiàn)對(duì)預(yù)報(bào)結(jié)果精度的分析，有效應(yīng)對(duì)核與輻射突發(fā)事件的各種復(fù)雜景象，提升應(yīng)對(duì)我國及鄰近地區(qū)核與輻射事故的應(yīng)急響應(yīng)能力。

1 國內(nèi)外研究現(xiàn)狀及進(jìn)展

1.1 國外研究現(xiàn)狀

在核事故后果評(píng)價(jià)系統(tǒng)應(yīng)用和發(fā)展方面，如美國ARAC/NARAC系統(tǒng)、歐共體RODOS系統(tǒng)、日本SPEEDI/WSPEEDI等國際先進(jìn)的核事故后果評(píng)價(jià)與決策支持系統(tǒng)，都非常關(guān)注中小尺度風(fēng)場(chǎng)預(yù)報(bào)、復(fù)雜景象大氣擴(kuò)散數(shù)值模擬技術(shù)的發(fā)展。隨著計(jì)算機(jī)和計(jì)算理論的發(fā)展，數(shù)值模式正在向全球化、精細(xì)化發(fā)展，以適應(yīng)多尺度、多目的應(yīng)用的要求。這些系統(tǒng)的結(jié)構(gòu)和功能不盡相同，但其共同特點(diǎn)是將中尺度氣象預(yù)報(bào)模式、風(fēng)場(chǎng)診斷模式和拉格朗日或歐拉類大氣擴(kuò)散模式集于一體，根據(jù)源項(xiàng)信息、在線氣象數(shù)據(jù)、地形地理以及劑量數(shù)據(jù)庫等，及時(shí)完成計(jì)算區(qū)域內(nèi)任意位置放射性物質(zhì)釋放的擴(kuò)散和劑量的精細(xì)計(jì)算，并可利用地理信息系統(tǒng)（GIS）動(dòng)態(tài)顯示計(jì)算結(jié)果和相關(guān)的應(yīng)急響應(yīng)信息。

在中小尺度核事故后果評(píng)價(jià)方面，ARAC/NARAC系統(tǒng)、RODOS系統(tǒng)、SPEEDI/WSPEEDI系統(tǒng)采用的主要模式見表2。

NARAC系統(tǒng)［4，5］由一套相互關(guān)聯(lián)的氣象和擴(kuò)散模式構(gòu)成。氣象預(yù)報(bào)模式包括美國國家大氣研究中心/國家環(huán)境預(yù)報(bào)中心（NCAR/NCEP）等聯(lián)合開發(fā)的新一代中尺度預(yù)報(bào)模式和同化系統(tǒng)WRF模式［6］，以及美國海軍艦隊(duì)數(shù)值氣象和海洋中心（FNMOC）的耦合海洋/大氣中尺度預(yù)報(bào)系統(tǒng)COAMPS模式。風(fēng)場(chǎng)診斷模式采用質(zhì)量守恒模式ADAPT。大氣擴(kuò)散采用Lagrangian隨機(jī)粒子擴(kuò)散模式LODI。

表2 用于中小尺度后果評(píng)價(jià)的主要模式Table 2 Models of mesoscale consequence assessment

RODOS系統(tǒng)［7］的氣象預(yù)報(bào)模式ADREA是一個(gè)非靜力動(dòng)力學(xué)風(fēng)場(chǎng)模式。風(fēng)場(chǎng)診斷模式應(yīng)用質(zhì)量守恒模式MCF和線性流動(dòng)模式LINCOM。大氣擴(kuò)散模式為具有輻射劑量計(jì)算功能的Lagrangian煙團(tuán)模型RIMPUFF和分段煙羽擴(kuò)散模型ATSTEP。

日本SPEEDI/WSPEEDI系統(tǒng)［8］采用的氣象預(yù)報(bào)模式PHYSIC是一個(gè)流體靜力學(xué)氣象預(yù)報(bào)模式。在風(fēng)場(chǎng)診斷和大氣擴(kuò)散模式方面，對(duì)日本國內(nèi)和國外核事故分別應(yīng)用不同的模式，風(fēng)場(chǎng)診斷為質(zhì)量守恒模式WIND21（國內(nèi)）和WSYNOP（國外）；大氣擴(kuò)散模式為Lagrangian擴(kuò)散模型PRWDA21（國內(nèi)）和GERAN（國外）。

近年來數(shù)值預(yù)報(bào)進(jìn)入了一個(gè)高速發(fā)展時(shí)期，隨著精細(xì)化要素預(yù)報(bào)在天氣預(yù)報(bào)業(yè)務(wù)中重要性的提升，高分辨率、高時(shí)效性的區(qū)域中尺度數(shù)值預(yù)報(bào)產(chǎn)品越來越受到關(guān)注。在精細(xì)化數(shù)值預(yù)報(bào)方面，各國都在積極研發(fā)高分辨率數(shù)值預(yù)報(bào)模式。中尺度數(shù)值預(yù)報(bào)系統(tǒng)的尺度分辨率已延伸至公里級(jí)。目前，一些發(fā)達(dá)國家的中尺度氣象預(yù)報(bào)模式系統(tǒng)已經(jīng)進(jìn)入實(shí)時(shí)應(yīng)用。其中，美國在該領(lǐng)域研發(fā)了多個(gè)模式，并在國際上得到廣泛應(yīng)用，最具代表性的是WRF模式，該模式的尺度分辨率為1 ～ 10 km。WRF模式采用高度模塊化、并行化和分層的設(shè)計(jì)技術(shù)，集成了迄今為止在中尺度方面最新的研究成果，模式中包括輻射過程、邊界層參數(shù)化過程、對(duì)流參數(shù)化過程、次網(wǎng)格湍流擴(kuò)散過程以及微物理過程等，是一個(gè)完全可壓的非靜力模式。另外，還有美國賓州大學(xué)/國家大氣研究中心（PSU/NCAR）聯(lián)合開發(fā)的中尺度數(shù)值預(yù)報(bào)模式第五代產(chǎn)品MM5模式，美國科羅拉多州立大學(xué)（CSU）等研發(fā)的區(qū)域大氣模擬系統(tǒng)RAMS模式等。

1.2 國內(nèi)研究進(jìn)展

在我國，經(jīng)過二十多年的發(fā)展，總體上可將開發(fā)的核事故后果評(píng)價(jià)系統(tǒng)分成兩代。20世紀(jì)90年代初中國輻射防護(hù)研究院為秦山核電廠開發(fā)的后果評(píng)價(jià)系統(tǒng)可作為第一代的代表，沒有氣象預(yù)報(bào)模式，大氣彌散采用高斯煙羽或拉格朗日煙團(tuán)等較為簡(jiǎn)單的模式。21世紀(jì)初開發(fā)的田灣核電廠核事故場(chǎng)外后果評(píng)價(jià)系統(tǒng)及其隨后的升級(jí)版本可作為第二代的代表，其標(biāo)志之一是RODOS系統(tǒng)的引進(jìn)與開發(fā)。

將中尺度氣象預(yù)報(bào)模式、風(fēng)場(chǎng)診斷模式和拉格朗日類軌跡大氣擴(kuò)散模式（M-C模式、Lagrangian煙團(tuán)模式）集于一體的技術(shù)路線已經(jīng)成為我國核事故后果評(píng)價(jià)軟件系統(tǒng)的發(fā)展趨勢(shì)。但是，與國外的中小尺度后果評(píng)價(jià)軟件相比，由于這些軟件開發(fā)的出發(fā)點(diǎn)基本都是針對(duì)某個(gè)固定的核設(shè)施，因此，不能預(yù)測(cè)模擬范圍內(nèi)任意位置發(fā)生的核與輻射事故后果。

氣象預(yù)報(bào)模式大多是對(duì)美國的MM5模式、RAMS模式和WRF等模式的引進(jìn)吸收與自主研發(fā)相結(jié)合的技術(shù)路線。國家氣象中心作為我國的數(shù)值氣象預(yù)報(bào)中心，承擔(dān)著國家級(jí)數(shù)值天氣預(yù)報(bào)科研及相關(guān)系統(tǒng)的業(yè)務(wù)運(yùn)行、維護(hù)和改進(jìn)工作。我國的國家級(jí)區(qū)域中尺度數(shù)值預(yù)報(bào)系統(tǒng)經(jīng)過多年的發(fā)展，從20世紀(jì)90年代的LAFS模式、HLAFS模式，直到2004年業(yè)務(wù)運(yùn)行的GRAPES區(qū)域中尺度模式，逐步從引進(jìn)吸收與自主研發(fā)并重轉(zhuǎn)入了自主研發(fā)、持續(xù)發(fā)展的新格局。

總體上，與國外的中小尺度后果評(píng)價(jià)軟件相比，雖然在區(qū)域氣象數(shù)值模式、大氣擴(kuò)散模式等方面有了一定的基礎(chǔ)，但是，由于軟件基本都是針對(duì)某個(gè)固定的核設(shè)施而開發(fā)的，目前尚缺乏精細(xì)預(yù)測(cè)我國和周邊國家及地區(qū)不同位置發(fā)生的核與輻射事故后果的軟件系統(tǒng)，不能滿足核事故應(yīng)急響應(yīng)的實(shí)際需求。

2 軟件研發(fā)的總體方案

采用成熟物理模型、移植、二次開發(fā)以及自主開發(fā)相結(jié)合的技術(shù)路線，針對(duì)我國和周邊國家及地區(qū)的地形和氣候特征，開發(fā)精細(xì)化核與輻射突發(fā)事件后果評(píng)價(jià)軟件系統(tǒng)。該系統(tǒng)以服務(wù)器、PC機(jī)為硬件平臺(tái)，主系統(tǒng)設(shè)計(jì)為一個(gè)基于三維GIS的，連接多子系統(tǒng)、多種模型和數(shù)據(jù)庫的綜合系統(tǒng)，用以預(yù)測(cè)、評(píng)價(jià)和動(dòng)態(tài)顯示我國和周邊國家及地區(qū)不同位置的輻射影響后果。軟件研發(fā)的技術(shù)路線如圖2所示，總體框架如圖3所示。

圖2 軟件研發(fā)的技術(shù)路線圖Fig.2 Technology program of the software development

圖3 軟件的總體框架圖Fig.3 General framework of the software

3 物理模塊的建立

3.1 中尺度氣象預(yù)報(bào)模式

（1）中尺度數(shù)值天氣預(yù)報(bào)模式GRAPES-MESO［9-13］

GRAPES-MESO是國家氣象中心在借鑒國際上先進(jìn)的通用模式的基礎(chǔ)上，自主研發(fā)的新一代中尺度氣象預(yù)報(bào)模式，經(jīng)過近十年的運(yùn)行、評(píng)估和完善，已實(shí)現(xiàn)業(yè)務(wù)化應(yīng)用，未來將成為我國在該領(lǐng)域的主流模式。其預(yù)報(bào)產(chǎn)品具有分辨率高、時(shí)效快、預(yù)報(bào)精度高等特點(diǎn)。

GRAPES-MESO模式采用三重嵌套，第一重嵌套區(qū)域覆蓋整個(gè)中國及附近海域，水平分辨率為27l km；第二重區(qū)域覆蓋華北大部分地區(qū)，水平分辨率為9 km；第三重區(qū)域覆蓋北京地區(qū)，水平分辨率為3 km。

GRAPES-MESO自2004年5月開始在國家氣象中心業(yè)務(wù)試運(yùn)行，2006年7月作為中國氣象局的業(yè)務(wù)系統(tǒng)投入運(yùn)行，并于同年8月初正式將GRAPES-30 km的產(chǎn)品下發(fā)到全國各級(jí)氣象臺(tái)站。2007年國家氣象中心根據(jù)實(shí)際業(yè)務(wù)需要，以及運(yùn)行效果評(píng)估報(bào)告對(duì)GRAPES-MESO模式進(jìn)行完善改進(jìn)，并將其水平分辨率由原來的30 km升級(jí)為15 km。針對(duì)系統(tǒng)業(yè)務(wù)化運(yùn)行中存在的問題，2009年再次對(duì)系統(tǒng)進(jìn)行了升級(jí)改進(jìn)。2010年3月新版本模式正式業(yè)務(wù)化，新版本模式對(duì)溫度場(chǎng)、風(fēng)場(chǎng)和降水的改進(jìn)效果比較顯著。

（2）GRAPES-MESO模式適用性研究

針對(duì)擬研發(fā)軟件的適用區(qū)域（ E60～150°、N5 ～ 55°），根據(jù)我國和周邊國家及地區(qū)的地形和氣候特征，通過選取不同地區(qū)模式物理參數(shù)化方案（如積云對(duì)流參數(shù)化方案、邊界層參數(shù)化方案、陸面過程等）和敏感性實(shí)驗(yàn)方法，開展氣象預(yù)報(bào)模式在不同區(qū)域的適用性研究，提高GRAPESMESO模式對(duì)不同區(qū)域的模擬能力（尤其是對(duì)風(fēng)場(chǎng)和降水的預(yù)報(bào)）。

（3）數(shù)據(jù)預(yù)處理和后處理程序接口

針對(duì)擬研發(fā)軟件的預(yù)報(bào)和回顧性評(píng)價(jià)兩大功能，開發(fā)支持多背景場(chǎng)的氣象數(shù)據(jù)接口。針對(duì)美國NCEP-GFS全球預(yù)報(bào)數(shù)據(jù)，國家氣象中心的T639全球預(yù)報(bào)數(shù)據(jù)以及美國NCEP-FNL全球再分析數(shù)據(jù)等，開發(fā)相應(yīng)的GRAPES-MESO預(yù)處理程序；針對(duì)以氣象預(yù)報(bào)模擬結(jié)果作為輸入的風(fēng)場(chǎng)診斷模型，開發(fā)相應(yīng)的GRAPES-MESO后處理程序。

（4）GRAPES-MESO模式運(yùn)行流程

首先，模式標(biāo)準(zhǔn)初始化系統(tǒng)對(duì)GFS/T639以及靜態(tài)資料等進(jìn)行預(yù)處理，為模式運(yùn)行提供背景場(chǎng)和側(cè)邊界條件；然后利用以動(dòng)力框架和物理過程參數(shù)化為主的模式系統(tǒng)進(jìn)行積分預(yù)報(bào)，最后得到圖形產(chǎn)品。其運(yùn)行流程如圖4所示。

圖4 GRAPES-MESO模式運(yùn)行流程Fig.4 Operational procedure of GRAPES-MESO

3.2 風(fēng)場(chǎng)診斷模式

由氣象預(yù)報(bào)模式生成的分辨率較粗的風(fēng)場(chǎng)信息，經(jīng)過氣象預(yù)報(bào)模式與風(fēng)場(chǎng)診斷模式的接口模塊，作為下一步風(fēng)場(chǎng)診斷模式的輸入資料，最終生成分辨率較細(xì)的風(fēng)場(chǎng)資料。氣象預(yù)報(bào)模式運(yùn)行結(jié)束后，將調(diào)用風(fēng)場(chǎng)診斷模式。圖5給出了風(fēng)場(chǎng)診斷模式的流程。擬基于CALMET氣象模式進(jìn)行風(fēng)場(chǎng)診斷模式移植與二次開發(fā)。

圖5 風(fēng)場(chǎng)診斷模式流程Fig.5 Procedure of wind field diagnostic model

CALMET氣象模式［14］是美國國家環(huán)保局（USEPA）長期支持開發(fā)的法規(guī)導(dǎo)則模型——三維非穩(wěn)態(tài)拉格朗日擴(kuò)散模式系統(tǒng)CALPUFF中的核心模塊之一。CALMET由風(fēng)場(chǎng)診斷模式和邊界層模式構(gòu)成。其風(fēng)場(chǎng)診斷分兩步進(jìn)行，第一步是風(fēng)場(chǎng)初始化，即利用氣象預(yù)報(bào)數(shù)據(jù)內(nèi)插外推出初始風(fēng)場(chǎng)。第二步是調(diào)整初始風(fēng)場(chǎng)，使其符合質(zhì)量守恒（即運(yùn)動(dòng)的連續(xù)性）條件的約束，同時(shí)使調(diào)整后的風(fēng)場(chǎng)與初始風(fēng)場(chǎng)的總體偏差為最小。

當(dāng)模擬區(qū)域很大時(shí)，可以用蘭伯特圓錐投影坐標(biāo)系統(tǒng)來校正地球的曲率。診斷風(fēng)場(chǎng)的計(jì)算內(nèi)容包括斜坡流、地形運(yùn)動(dòng)學(xué)矯正、地形阻擋效應(yīng)和差異最小化擬合。經(jīng)過以上計(jì)算獲得三維風(fēng)場(chǎng)。

3.3 非穩(wěn)態(tài)、多源拉格朗日煙團(tuán)模式（LPM）

大氣擴(kuò)散模型采用非穩(wěn)態(tài)、多源拉格朗日煙團(tuán)模型（LPM），該模型是在丹麥Riso實(shí)驗(yàn)室開發(fā)的RIMPUFF模式［15］的基礎(chǔ)上進(jìn)行移植與二次開發(fā)合并完成的，該模型能夠很好地處理非靜態(tài)和非均勻的氣象條件，考慮復(fù)雜地形上風(fēng)切變的影響，模型的核心是對(duì)煙團(tuán)的平流、擴(kuò)散和沉積采用局地氣象數(shù)據(jù)進(jìn)行模擬，并能在一定程度上處理復(fù)雜地形下煙流的分裂，可以應(yīng)用于平坦的或不十分復(fù)雜的地形條件。

在本軟件的研發(fā)過程中，擬根據(jù)核電廠事故特點(diǎn)以及日本福島事故的經(jīng)驗(yàn)，對(duì)模式在可模擬的釋放情景、靈活性等方面做較大的改進(jìn)，以有效應(yīng)用于復(fù)雜情景下的后果評(píng)價(jià)和回顧性評(píng)價(jià)。改進(jìn)后的LPM模式主要特點(diǎn)如下：

（1）可響應(yīng)我國和周邊國家及地區(qū)不同位置的氣載放射性物質(zhì)釋放。

（2）可同時(shí)響應(yīng)多源事故釋放情形。

（3）可處理不均勻（非穩(wěn)態(tài)）釋放源項(xiàng)。

（4）事故后果評(píng)價(jià)的模擬時(shí)間可變，只要在可利用的氣象數(shù)據(jù)（包括實(shí)測(cè)、預(yù)報(bào)數(shù)據(jù)）許可范圍內(nèi)即可。

（5）還可利用不同的氣象數(shù)據(jù)開展多種源項(xiàng)特征下的后果評(píng)價(jià)，以適應(yīng)不同情形下的響應(yīng)：①利用氣象預(yù)報(bào)數(shù)據(jù)，適合預(yù)測(cè)事故后果；②利用再分析數(shù)據(jù)，適合事故后開展回顧性評(píng)價(jià)；③通過預(yù)留實(shí)測(cè)氣象數(shù)據(jù)的接口，提高計(jì)算精度。

3.4 劑量估算模式

根據(jù)放射性核素?cái)U(kuò)散濃度分布的模擬結(jié)果，計(jì)算出事故情況下氣載放射性物質(zhì)在不同時(shí)刻可能引起的環(huán)境潛在輻射劑量水平。輻射劑量包括煙云浸沒外照射劑量、地面沉積外照射劑量和煙云吸入內(nèi)照射劑量。在劑量估算模式中考慮氣載放射性物質(zhì)遷移擴(kuò)散過程中的干、濕沉積，源項(xiàng)損耗等因素。為使軟件適應(yīng)不同事故情形下的氣載放射性物質(zhì)釋放，模式中考慮的放射性核素包括核爆、核設(shè)施突發(fā)事件可能產(chǎn)生的核素。

4 展望

快速、準(zhǔn)確地預(yù)測(cè)發(fā)生在我國及鄰近地區(qū)不同位置的核與輻射事故后果，是及時(shí)制定保護(hù)我國公眾和環(huán)境的應(yīng)急對(duì)策、適時(shí)采取減輕核事故后果的防護(hù)行動(dòng)的重要依據(jù)。因此，研發(fā)精細(xì)化核事故氣載放射性物質(zhì)擴(kuò)散的數(shù)值模擬軟件，是核事故應(yīng)急決策中不可或缺的技術(shù)手段。

本文針對(duì)擬建的精細(xì)化核事故應(yīng)急氣載放射性物質(zhì)擴(kuò)散的數(shù)值模擬軟件，給出了軟件研發(fā)的技術(shù)路線、主要物理模塊的建立方法等，為軟件的順利研發(fā)奠定了堅(jiān)實(shí)的技術(shù)基礎(chǔ)。結(jié)合核與輻射應(yīng)急響應(yīng)工作的實(shí)際需求，在未來的工作中，還需在以下方面開展深入的研究：

擬建軟件具有預(yù)報(bào)和回顧性評(píng)價(jià)兩大功能，如何利用回顧性評(píng)價(jià)結(jié)果，提高預(yù)報(bào)結(jié)果的精度，進(jìn)而提高軟件的整體性能。

目前擬建軟件的適用區(qū)域主要是我國及鄰近地區(qū)。為進(jìn)一步提高我國應(yīng)對(duì)全球其他國家及地區(qū)的核應(yīng)急響應(yīng)能力，需逐步地、持續(xù)地?cái)U(kuò)展軟件的適用范圍。

［1］IAEA. IAEA International Fact Finding Expert Mission Of The Fukushima Dai-Ichi NPP Accident Following The Great East Japan Earthquake and Tsunami［R］. Tokyo：The Great East Japan Earthquake Expert Mission，2011.

［2］Mary B N. North Korea’s Nuclear Weapons：Technical Issues［R］. Washington D.C.：Congressional Research Service，2013.

［3］楊玲，岳會(huì)國，林權(quán)益，等．對(duì)日本福島核事故的應(yīng)急響應(yīng)及對(duì)中心應(yīng)急準(zhǔn)備工作的改進(jìn)建議［J］．核安全，2011（3）：33-38.

［4］陳曉秋，李冰，林權(quán)益．核事故應(yīng)急放射性后果評(píng)價(jià)模式的有效性［J］. 核安全，2008（3）：41-43，62.

［5］李宗明．從日本福島核事故審視核安全的政府、法律和監(jiān)管框架［J］．核安全，2012（2）：1-8.

［6］Sullivan T J，Ellis J S，F(xiàn)oster C S，et al. Atmospheric release advisory capability：Real-time modelling airborne hazardous materials［J］. Bulletin of the American Meteorological Society，1993，74（12）：23-43.

［7］Nasstrom J S，Sugiyama G，Baskett R，et a1. The national atmospheric release advisory center（NARAC） modeling and decision Support System for radiological and nuclear emergency preparedness and response［R］．Livermore：NARAC，2005.

［8］Skamarock W C，Klemp J B，Dudhia J，et a1．A description of the advanced research WRF［M］．3rd ed．Colorado：National Center for Atmospheric Research，2008．

［9］Ehrhardt J，Weis A. Development of RODOS，A comprehensive real-time on-line decision support for nuclear emergencies management in Europe［R］. Karlsruhe：Karlsruhe Institute of Technology （KIT），1996.

［10］Chino M，Ishikawa H，Yamazawa H. SPEEDI and WSPEEDI Japanese emergency response system to predict radiological impacts in local and workplace areas due to a nuclear accident［J］. Radiation Protection Dosimetry，1993，50：2.

［11］王雨，李莉．GRAPES-Meso V3.0模式預(yù)報(bào)效果檢驗(yàn)［J］．應(yīng)用氣象學(xué)報(bào)，2010，21（5）：524-534．

［12］陳德輝，薛紀(jì)善，楊學(xué)勝，等．GRAPES新一代全球／區(qū)域多尺度統(tǒng)一數(shù)值預(yù)報(bào)模式總體設(shè)計(jì)研究［J］.科學(xué)通報(bào)，2008，53（20）：2396-2407.

［13］葉成志，歐陽里程，李象玉，等．GRAPES中尺度模式對(duì)2005年長江流域重大災(zāi)害性天氣過程預(yù)報(bào)性能的檢驗(yàn)分析［J］．熱帶氣象學(xué)報(bào)，2006，26（4）：393-399.

［14］朱振鐸，端義宏，陳德輝．GRAPES TCM業(yè)務(wù)試驗(yàn)結(jié)果分析［J］．氣象，2007，33（7）：44-54．

［15］李勇，王雨．2007年夏季 GRAPES-Meso 15及30 km模式對(duì)比檢驗(yàn)［J］．氣象，2008，34（10）：81-89.

［16］Scire J S，Robe F R，F(xiàn)ernau M E，et a1．A user’s guide for the CALMET meteorological model［R］. 5th ed．San Francisco：Earth Technology Incorporated，2000.

［17］Nielsen S T，Deme S，Mikkelsen T. Description of the Atmospheric Dispersion Module RIMPUFF ［R］.Roskilde：Ris? National Laboratory， 1999.