袁偉，曲延濤，楊宏偉，左莉，陳顯波

（1.第二炮兵裝備研究院，北京 100085；2.中國原子能科學(xué)研究院，北京 102413）

基于部隊(duì)核事故后果分析的系統(tǒng)開發(fā)

袁偉1，曲延濤2，楊宏偉2，左莉1，陳顯波1

（1.第二炮兵裝備研究院，北京 100085；2.中國原子能科學(xué)研究院，北京 102413）

通過分析國內(nèi)核設(shè)施事故后果評價(jià)技術(shù)現(xiàn)狀及相關(guān)系統(tǒng)的功能特征，研究部隊(duì)核事故后果評價(jià)技術(shù)的發(fā)展情況，比較兩者在評價(jià)技術(shù)方法和功能模式上的異同點(diǎn)。在此基礎(chǔ)上，針對前期研制的一種適用于部隊(duì)核事故后果評價(jià)的軟件系統(tǒng)，在功能組成、開發(fā)應(yīng)用和優(yōu)化升級等方面進(jìn)行了探討，為部隊(duì)核武器事故后果評價(jià)技術(shù)發(fā)展提供了啟示。

核事故；后果評價(jià)；系統(tǒng)開發(fā)

當(dāng)核設(shè)施發(fā)生嚴(yán)重的放射性物質(zhì)泄漏事故時(shí)，有必要采取有針對性的行動(dòng)措施減少或避免公眾可能受到的放射性照射。這些年來，國內(nèi)相關(guān)核工業(yè)單位在后果評價(jià)方面開展了一些有效的工作，取得了明顯的成效。同時(shí)，各個(gè)研究中心和科研機(jī)構(gòu)也針對核設(shè)施運(yùn)行狀況以及周邊環(huán)境條件，收集大量的資料信息獲得的實(shí)用數(shù)據(jù)能夠用于完善后果評價(jià)模型，也可對同類核事故應(yīng)急處理提供處置行動(dòng)建議，使人員和環(huán)境得到更好的保護(hù)。

軍隊(duì)武器雖然經(jīng)過長時(shí)間的發(fā)展，其安全性和可靠性得到了極大提升，但是仍無法完全避免事故的發(fā)生。一旦發(fā)生核事故，需要在較短時(shí)間內(nèi)進(jìn)行事故后果評估，為后續(xù)應(yīng)急處置行動(dòng)提供決策建議。目前，國內(nèi)部隊(duì)核事故后果評價(jià)軟件系統(tǒng)的研究工作，還處于起步階段，許多基礎(chǔ)性的研究工作急需開展。本文主要針對前期研制的一種適用于部隊(duì)核事故后果評價(jià)的軟件系統(tǒng)，在功能組成、開發(fā)應(yīng)用和優(yōu)化升級等方面進(jìn)行了探討，為部隊(duì)核事故后果評價(jià)技術(shù)發(fā)展提供了一定的啟示。

1 民用核事故后果評價(jià)技術(shù)現(xiàn)狀

目前，國內(nèi)研究的核事故后果評價(jià)系統(tǒng)可分為兩個(gè)階段。民用最早的核電站（大亞灣和秦山核電站）于20世紀(jì)80年代末研制的核事故后果評價(jià)軟件，可作為第一階段的典型代表。由于當(dāng)時(shí)技術(shù)條件限制和科研水平薄弱，僅采取了比較簡單的分析方法，包括高斯煙羽大氣彌散模式和拉格朗日煙團(tuán)模式。

其與現(xiàn)階段的分析水平相比已落后，并且與實(shí)際條件差別比較大［1］。

隨后，國內(nèi)逐漸出現(xiàn)了以地理信息系統(tǒng)為平臺的核事故后果評價(jià)分析系統(tǒng)，可作為第2階段的代表。這類系統(tǒng)結(jié)合核事故發(fā)生區(qū)域的實(shí)際地理環(huán)境條件，發(fā)展了不同風(fēng)場的預(yù)報(bào)模式，大氣彌散主要采用三維蒙特卡羅粒子模式，劑量模塊包括有預(yù)期劑量、可防止劑量的計(jì)算等。系統(tǒng)分析結(jié)果貼近實(shí)際事故場景，計(jì)算誤差較小，已取代了早期評價(jià)系統(tǒng)。

2 部隊(duì)核事故后果評價(jià)特點(diǎn)

由于軍隊(duì)武器自身的特殊性，因此其發(fā)生核事故與民用核設(shè)施核事故相比，在后果評價(jià)方面具有以下顯著特點(diǎn)：

2.1 事故評價(jià)的源項(xiàng)核素特殊

核設(shè)施使用的核燃料鈾235的豐度僅占3%～5%，而軍隊(duì)武器中所使用的鈾材料，鈾235的豐度達(dá)90%以上。軍隊(duì)武器中還使用了極毒钚材料部件，钚是一種具有強(qiáng)放射性的物質(zhì)，其致死劑量為0.6微克。另外，還含有滲透性很強(qiáng)的氚材料，對防護(hù)要求很高。

2.2 事故后果評價(jià)模式復(fù)雜

核電站、核設(shè)施主要發(fā)生放射性物質(zhì)泄漏、堆芯熔化等事故，釋放源為穩(wěn)定源，可采用近似模型進(jìn)行分析。而軍隊(duì)武器中含有大量的火工品和炸藥部件，一旦發(fā)生爆炸等事故，放射性物質(zhì)隨著爆炸產(chǎn)物瞬間釋放出來，初始狀態(tài)物質(zhì)擴(kuò)散模型不易建立，放射性物質(zhì)運(yùn)動(dòng)規(guī)律和軌跡難以模擬。

2.3 事故后果評估技術(shù)要求高

軍隊(duì)武器是集高、精、尖技術(shù)于一體的武器裝備，內(nèi)部結(jié)構(gòu)精細(xì)，事故誘因多，事故機(jī)理復(fù)雜。一旦發(fā)生事故，對事故后果評估的全面性和準(zhǔn)確性有很高的技術(shù)要求，需要建立準(zhǔn)確的評估技術(shù)方法以指導(dǎo)后續(xù)應(yīng)急處置行動(dòng)。

3 系統(tǒng)主要功能

3.1 系統(tǒng)介紹

多年來結(jié)合部隊(duì)開展的核事故后果評價(jià)研究工作，在“十二五”初期通過與中國原子能科學(xué)研究院合作，開發(fā)研制了一套適用于軍隊(duì)武器事故后果評價(jià)軟件系統(tǒng)，適用于軍隊(duì)武器作戰(zhàn)使用過程中發(fā)生核事故的機(jī)理及模擬技術(shù)研究，能夠?qū)χ饕耸鹿暑愋偷挠绊懸蛩丶爱a(chǎn)生的影響進(jìn)行預(yù)測，為軍隊(duì)武器使用安全研究提供基礎(chǔ)技術(shù)支持（圖1）。

該軟件系統(tǒng)主要用途包括：能夠提供主要事故（地面和近地面化爆、核材料泄露、火災(zāi)等事故模式）計(jì)算模塊及接入接口；能夠提供各種典型源項(xiàng)輸入模塊接口；能夠根據(jù)氣象站的觀測數(shù)據(jù)進(jìn)行風(fēng)場預(yù)報(bào)計(jì)算，給出預(yù)報(bào)溫度場、湍能、風(fēng)場和壓強(qiáng)的具體分布；能夠根據(jù)預(yù)報(bào)的風(fēng)場和氣流特點(diǎn)，分析給出核事故條件下氣溶膠的運(yùn)動(dòng)過程、煙羽濃度梯度分布；可對應(yīng)急處置人員輻照評估和公眾劑量后果評價(jià)；能夠在數(shù)據(jù)地圖上同步展現(xiàn)二維的計(jì)算結(jié)果；能夠在數(shù)據(jù)地圖上同步展現(xiàn)三維空間氣溶膠擴(kuò)散沉降情況；報(bào)告、報(bào)表生成管理功能。

3.2 模塊組成

該系統(tǒng)主要包括源項(xiàng)數(shù)據(jù)模塊、地理信息模塊、氣象診斷模塊、氣象預(yù)報(bào)模塊、污染物擴(kuò)散模塊、后果評價(jià)模塊、對策優(yōu)化模塊和實(shí)時(shí)圖形展示模塊等部分。模塊之間的相互串聯(lián)以Net Framework微軟官方版本為平臺，以VC作為開發(fā)程序語言，以SQL Server作為數(shù)據(jù)管理系統(tǒng)（圖2）。

3.3 模塊功能

3.3.1源項(xiàng)數(shù)據(jù)模塊

針對軍隊(duì)武器不同的事故情景模式（爆炸、火燒、放射性泄漏等），開發(fā)不同的源項(xiàng)數(shù)據(jù)，這是整個(gè)系統(tǒng)評價(jià)的核心。爆炸源項(xiàng)描述地面或高空爆炸事故產(chǎn)生的瞬時(shí)高溫放射性物質(zhì)，主要數(shù)據(jù)包括爆炸當(dāng)量、爆炸相對地面高度、釋放核素分類、放射性核素釋放量等；火燒源項(xiàng)描述火災(zāi)事故產(chǎn)生的連續(xù)高溫放射性物質(zhì)，主要數(shù)據(jù)包括燃燒物質(zhì)種類，燃燒量、燃燒時(shí)間、釋放核素分類、放射性核素釋放量等；放射性泄漏源項(xiàng)主要包括放射性核素釋放高度、釋放功率、釋放速率、釋放核素分類、釋放量、事故時(shí)間等［2］（圖3）。

3.3.2地理信息模塊

基于事故發(fā)生區(qū)域的ArcGIS9.0地理信息（包括人口分布、特殊單位、公鐵路分布、高程參數(shù)）的數(shù)據(jù)管理平臺，既能夠方便儲(chǔ)存、修改、更新和查詢數(shù)據(jù)，也能夠定制應(yīng)急有關(guān)的技術(shù)數(shù)據(jù)，并且為輔助決策以及應(yīng)急評估提供有力的數(shù)據(jù)服務(wù)（圖4）。

3.3.3氣象診斷模塊

無論是定性分析放射性氣溶膠的運(yùn)動(dòng)情況，還是氣溶膠分布模式定量分析，都需要掌握當(dāng)?shù)禺?dāng)時(shí)的氣象數(shù)據(jù)。在地面粗糙度較高情況下，邊界層周圍的風(fēng)場的時(shí)空變化很大。風(fēng)場診斷結(jié)果的準(zhǔn)確性，取決于實(shí)際測量風(fēng)場的精度、密度以及代表性。在中尺度風(fēng)場問題分析中，由于實(shí)際條件所限，在一個(gè)模擬范圍內(nèi)不易獲得數(shù)個(gè)風(fēng)速輪廓的測量數(shù)據(jù)，一般利用地面風(fēng)的觀測數(shù)據(jù)向上推演，以代替風(fēng)速輪廓線［3］。本系統(tǒng)的風(fēng)場診斷模塊可以處理100 km×100 km范圍，高2 000m的區(qū)域風(fēng)場診斷問題，計(jì)算精度取決于計(jì)算區(qū)域的氣象觀測站點(diǎn)的密度和高質(zhì)量的氣象觀測數(shù)據(jù)（圖5）。

3.3.4氣象預(yù)報(bào)模塊

核事故應(yīng)急響應(yīng)、評價(jià)分析中涉及的預(yù)報(bào)風(fēng)場，一般是中等尺度的氣象模擬。中尺度可定義為水平尺度比天氣尺度小，但比單個(gè)積云尺度要大的尺度。水平尺度從幾公里到幾百公里，時(shí)間尺度約為1到12 h，垂直尺度則從幾十米到整個(gè)對流層。本系統(tǒng)采用非靜力、完成可壓動(dòng)力學(xué)模型和交錯(cuò)網(wǎng)格技術(shù)，根據(jù)氣象特征、初始風(fēng)場、邊界條件預(yù)報(bào)24 h內(nèi)的三維風(fēng)場演變情況［4］。預(yù)報(bào)風(fēng)場模塊采用風(fēng)場診斷模塊的結(jié)果作為初始風(fēng)場，氣象特征需要用戶進(jìn)行選擇和輸入，計(jì)算區(qū)域要小于風(fēng)場診斷的區(qū)域，默認(rèn)為80 km× 80 km。

3.3.5污染物擴(kuò)散模塊

采用粒子隨機(jī)游走模式來模擬放射性污染物在大氣中的擴(kuò)散，該模式可以模擬復(fù)雜形狀源項(xiàng)在復(fù)雜地形情況下的擴(kuò)散問題，適用面廣，計(jì)算精度高，并可以跟蹤污染物的遷移軌跡，能夠?yàn)槿S動(dòng)態(tài)展示提供數(shù)據(jù)。系統(tǒng)還通過不同精度計(jì)算網(wǎng)格的給出，反映污染物的不同濃度，包括近地面煙云的即時(shí)濃度、近地面煙云的積累濃度、地面沉積濃度和地面濕沉積濃度等（圖6）。

3.3.6后果評價(jià)模塊

核事故后果評價(jià)的根本目的是評價(jià)核事故后產(chǎn)生的大量放射性污染物對人的影響，在輻射防護(hù)中采用當(dāng)量劑量這個(gè)物理量來評價(jià)不同輻射對人體不同組織器官的影響，對不同組織器官的當(dāng)量劑量進(jìn)行加權(quán)平均就得到有效當(dāng)量劑量。評價(jià)事故后放射性物質(zhì)通過不同途徑對人的有效劑量是事故后果評價(jià)的重要任務(wù)。本系統(tǒng)根據(jù)污染物的擴(kuò)散濃度，沉積濃度，計(jì)算不同途徑的有效劑量，包括空氣浸沒外照射劑量率（各組織器官的劑量）、地面沉積外照射劑量率（干沉積和濕沉積后各組織器官的劑量）、吸入內(nèi)照射劑量（各組織器官的劑量）和食入內(nèi)照射劑量（各組織器官的劑量）（圖7）。

3.3.7對策優(yōu)化模塊

該模塊的實(shí)現(xiàn)主要分3步完成：第1步主要為生成對策，對策的生成將結(jié)合后果評價(jià)的計(jì)算結(jié)果和地理基礎(chǔ)信息來完成，生成每個(gè)網(wǎng)格可能的行動(dòng)方案；第2步是在考慮各種相關(guān)因素的條件下，用各種不同的優(yōu)化評估方法對可能的對策進(jìn)行優(yōu)化評估和排序，得到盡可能優(yōu)的決策方案；第3步對各種優(yōu)化后的決策方案做進(jìn)一步的代價(jià)分析、對比和排序，以供決策者參考。

3.3.8實(shí)時(shí)圖形展示模塊

該模塊能夠?qū)崿F(xiàn)對空間相關(guān)數(shù)據(jù)的表現(xiàn)。二維渲染圖包括地面沉積濃度、近地表時(shí)間積分濃度、空氣浸沒劑量率、地面沉積劑量率和吸入劑量的計(jì)算結(jié)果，渲染圖將用不同顏色標(biāo)注計(jì)算結(jié)果的大小，并用等值線給出不同結(jié)果的分布區(qū)域。三維擴(kuò)散效果圖采用三維模擬對風(fēng)場和放射性污染物擴(kuò)散軌跡進(jìn)行展示，風(fēng)場將在ArcGIS地圖的網(wǎng)格點(diǎn)上用帶箭頭的線段表示，線段的長度表示風(fēng)的大小，箭頭指向風(fēng)的方向，污染物的擴(kuò)散軌跡將用帶顏色的小球模擬污染物粒子的運(yùn)動(dòng)［5-6］。

4 結(jié)論

運(yùn)用計(jì)算機(jī)軟件平臺及有關(guān)程序語言，筆者開發(fā)完成了適用于部隊(duì)使用的核事故后果仿真分析系統(tǒng)。該系統(tǒng)軟件界面友好，操作簡便且功能實(shí)用，充分發(fā)揮了ArcGIS9.0軟件在地理空間數(shù)據(jù)管理和分析方面的優(yōu)勢，易于部隊(duì)核應(yīng)急專業(yè)人員學(xué)習(xí)掌握，為開展核事故的物理模型、發(fā)生模式和影響因素分析，以及大空間尺度事故模擬技術(shù)的研究提供技術(shù)支持。在“十二五”初期該套軟件已通過驗(yàn)收，下階段的主要工作將重點(diǎn)圍繞軍隊(duì)武器不同核事故場景條件下的源項(xiàng)數(shù)據(jù)（如放射性核素釋放高度、釋放功率、釋放速率、釋放核素分類和釋放量等）以及積累局部地區(qū)氣象數(shù)據(jù)，開展進(jìn)一步深入研究和探討。

［1］王醒宇，康凌.核事故后果評價(jià)方法及其新發(fā)展［M］.北京：原子能出版社，2003.

［2］胡二邦.大氣環(huán)境研究若干領(lǐng)域的動(dòng)態(tài)與進(jìn)展［J］.輻射防護(hù)通訊，2001，21（6）:13-21.

［3］曲靜原，曹建主.RODOS4.0與RODOS的未來發(fā)展［J］.輻射防護(hù)，2002，22（4）:2-8.

［4］田鋼，王醒宇.廣東省核電站核事故場外后果預(yù)測評價(jià)系統(tǒng)［J］.輻射防護(hù)，2001，21（4）:17-21.

［5］李紅，徐火根，陳立航，等.核設(shè)施正常工況下氣載放射性排出物環(huán)境影響評價(jià)軟件的設(shè)計(jì)與開發(fā)［J］.核動(dòng)力工程，2003，21（4）:209-212.

［6］樊申.秦山核電廠全廠斷電事故研究和廠外后果分析［D］.北京：中國原子能科學(xué)研究院，2005.

The system development of consequence analysis for arm y nuclear accident

YUANWei，QU Yan-tao，YANG Hong-wei，ZUO Li，CHEN Xian-bo

（1.The Secondary Artillery Equipment Academy，Beijing 100085，China；
2.China Atomic Energy Institute，Beijing 102413，China）

Through the analysis of the assessment technique and the function characteristics of related system，domestic nuclear facilities accident consequence，assessment technique for army nuclear accident consequence was studied and the difference between the two evaluation methods were also compared.On this basis，the function，the application and the optimization of a software system that was developed for army was discussed.All of these will provide suggestion for the nuclear weapon accident consequence assessment.

nuclear accident；consequence assessment；system development

P631.1+2

A

1672-0636（2014）02-0120-05

10.3969/j.issn.1672-0636.2014.02.010

2013-04-19；

2013-04-27

袁偉（1983—），男，江蘇南京人，工程師，研究方向：輻射防護(hù)。

E-mail:wjg.0031@163.com