程詩(shī)思，吳 晗

（1.中國(guó)核動(dòng)力研究設(shè)計(jì)院核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室，成都 610041；2.環(huán)境保護(hù)部核與輻射安全中心，北京 100082）

小型動(dòng)力堆彈棒事故環(huán)境釋放源項(xiàng)分析

程詩(shī)思1，吳 晗2，*

（1.中國(guó)核動(dòng)力研究設(shè)計(jì)院核反應(yīng)堆系統(tǒng)設(shè)計(jì)技術(shù)重點(diǎn)實(shí)驗(yàn)室，成都 610041；2.環(huán)境保護(hù)部核與輻射安全中心，北京 100082）

彈棒事故作為壓水堆設(shè)計(jì)基準(zhǔn)事故，包括兩種放射性向環(huán)境釋放的途徑：安全殼泄漏及汽輪機(jī)、安全閥釋放。本文以小型動(dòng)力堆為研究對(duì)象，參考AP1000的最新事故源項(xiàng)分析思路，并與M 310型核電廠的分析思路進(jìn)行對(duì)比，得出小型動(dòng)力堆彈棒事故的環(huán)境釋放源項(xiàng)分析方法，計(jì)算了彈棒事故造成的環(huán)境釋放源項(xiàng)。

小型動(dòng)力堆；彈棒事故；源項(xiàng)

在放射性核素和裝置失控的情況下，可能對(duì)工作人員、公眾和周圍環(huán)境造成輻射危害［1］。核電廠運(yùn)行狀態(tài)和事故工況下的源項(xiàng)分析是核電廠安全分析的重要組成部分，進(jìn)行源項(xiàng)分析主要是因?yàn)椋海?）確保設(shè)計(jì)是最優(yōu)化的，使源項(xiàng)合理可行盡量低；（2）證明設(shè)計(jì)能夠滿足輻射防護(hù)的要求，包括劑量限制；（3）理論源項(xiàng)作為應(yīng)急所需測(cè)量參數(shù)的早期指示，可為制定應(yīng)急計(jì)劃提供依據(jù)；（4）證明設(shè)計(jì)基準(zhǔn)事故所要求的設(shè)備具有足夠的可靠性，包括工具和氣體處理系統(tǒng)［2］?；谑鹿试错?xiàng)分析，進(jìn)而可以開(kāi)展事故放射性后果評(píng)價(jià)、設(shè)備鑒定、屏蔽設(shè)計(jì)等［3］。本文研究了小型動(dòng)力堆彈棒事故環(huán)境釋放源項(xiàng)分析方法，就計(jì)算方法等與M 310型核電廠進(jìn)行了比較，并給出了小型動(dòng)力堆彈棒事故環(huán)境釋放源項(xiàng)的計(jì)算結(jié)果。

1 計(jì)算方法

彈棒事故是由于控制棒驅(qū)動(dòng)機(jī)構(gòu)耐壓殼機(jī)械損壞，導(dǎo)致控制棒組件和驅(qū)動(dòng)軸彈出堆芯外［4］，為極限事故工況［5］。彈棒事故發(fā)生后，事故前和事故期間因燃料元件包殼破損而從破損燃料元件釋放到反應(yīng)堆冷卻劑中的裂變產(chǎn)物同時(shí)通過(guò)兩種途徑釋放到環(huán)境：途徑一，反應(yīng)堆冷卻劑內(nèi)的裂變產(chǎn)物隨反應(yīng)堆冷卻劑通過(guò)控制棒驅(qū)動(dòng)機(jī)構(gòu)耐壓殼破口釋放到安全殼大氣中，安全殼大氣內(nèi)的裂變產(chǎn)物再通過(guò)安全殼泄漏釋放到環(huán)境；途徑二，反應(yīng)堆冷卻劑內(nèi)的裂變產(chǎn)物隨反應(yīng)堆冷卻劑通過(guò)蒸汽發(fā)生器傳熱管泄漏到蒸汽發(fā)生器二次側(cè)工質(zhì)中，蒸汽發(fā)生器二次側(cè)工質(zhì)內(nèi)的裂變產(chǎn)物再通過(guò)汽輪機(jī)和蒸汽發(fā)生器安全閥釋放到環(huán)境中［6］。

對(duì)于小型動(dòng)力堆，采用分時(shí)間步的計(jì)算方法，根據(jù)事故分析得到的事故序列，由各時(shí)間點(diǎn)的破口流量、安注流量等數(shù)據(jù)計(jì)算事故后各時(shí)間點(diǎn)裂變產(chǎn)物在反應(yīng)堆冷卻劑中的比活度及向安全殼大氣的釋放速率，進(jìn)而得到通過(guò)途徑一造成的環(huán)境釋放源項(xiàng)。

而對(duì)于M 310型核電廠，在計(jì)算通過(guò)釋放途徑一造成的環(huán)境釋放源項(xiàng)時(shí)，采用包絡(luò)的瞬時(shí)釋放模型，不用再單獨(dú)考慮事故序列，現(xiàn)實(shí)模型中假設(shè)的惰性氣體及碘從破損燃料元件氣隙釋放到反應(yīng)堆冷卻劑的釋放份額（除85Kr外的惰性氣體：2%；85K r：30%；碘：3%），就是包絡(luò)的份額。

小型動(dòng)力堆及M 310型核電廠都是采用蒸汽發(fā)生器一次側(cè)向二次側(cè)的泄漏率、蒸汽發(fā)生器向環(huán)境的釋放流量、蒸汽發(fā)生器的給水流量等數(shù)據(jù)計(jì)算蒸汽發(fā)生器二次側(cè)水和蒸汽中放射性核素的比活度，進(jìn)而得到通過(guò)釋放途徑二造成的環(huán)境釋放源項(xiàng)。

此外，對(duì)于釋放到安全殼內(nèi)的碘的去除機(jī)制，小型動(dòng)力堆參考了AP1000的思路：設(shè)計(jì)中不依靠能動(dòng)系統(tǒng)去除事故后釋放到安全殼大氣中的元素碘或氣溶膠。元素碘通過(guò)在安全殼內(nèi)結(jié)構(gòu)件表面的沉積被去除，考慮去污因子；安全殼內(nèi)氣溶膠通過(guò)重力沉降去除，考慮去污因子。

在AP1000彈棒事故源項(xiàng)分析中，元素碘去污因子為200，氣溶膠去污因子為1000［7］。

由于元素碘通過(guò)表面沉積去除，當(dāng)單位體積內(nèi)元素碘的個(gè)數(shù)少到某種程度時(shí)，沉積到壁面的概率很小，不應(yīng)再考慮元素碘的去除。因此小型動(dòng)力堆不直接采用AP1000的去污因子，而是根據(jù)AP1000彈棒事故后元素碘去污因子達(dá)到200時(shí)，安全殼大氣中元素碘的濃度，得到小型動(dòng)力堆彈棒事故后安全殼大氣中元素碘達(dá)到該濃度時(shí)所對(duì)應(yīng)的去污因子［8］。

對(duì)于氣溶膠，AP1000由方程e-λt=1/1 000得到去污因子達(dá)到1 000所需的時(shí)間。式中，λ為氣溶膠去除系數(shù)［9］。計(jì)算去污因子達(dá)到1 000所需的時(shí)間時(shí)，只考慮重力沉降對(duì)安全殼大氣中氣溶膠濃度的影響，不考慮衰變及安全殼泄漏。上述方程是基于瞬時(shí)釋放模型，小型動(dòng)力堆參考上述方法，基于分時(shí)間步釋放模型，得到事故后各時(shí)間步安全殼大氣中氣溶膠的濃度及安全殼大氣中氣溶膠的濃度達(dá)到最大濃度的1/1 000所需時(shí)間，即氣溶膠去除持續(xù)時(shí)間。安全殼大氣中氣溶膠的濃度也只考慮重力沉降的去除，不考慮衰變及安全殼泄漏。

而對(duì)于M 310型核電廠，在安全殼內(nèi)考慮的碘的去除機(jī)制包括重力沉降和安全殼噴淋除碘，噴淋系統(tǒng)在事故后手動(dòng)啟動(dòng)［10］。安全殼噴淋的作用是去除元素碘；對(duì)粒子碘不起作用［11］。

2 計(jì)算程序及模型

計(jì)算程序?yàn)門ACT［12］及GVACT［13］。

（1）TACTIII程序計(jì)算模型如下：

式中，SAnm，i：在第m個(gè)時(shí)間步，核素n在控制容積i（安全殼）凈空間中的放射性濃度；

Vi：控制容積i凈空間體積；

λn：核素n的衰變常數(shù)；

βn，v：核素v衰變成核素n的分支比；

m：時(shí)間步，（m=1，2，3，…，M）；

n：衰變鏈中核素的標(biāo)號(hào)，（n=1，2，3，…，N）；

v：衰變鏈中核素的標(biāo)號(hào)，（1≤v＜n）。

（2）GVACT程序計(jì)算事故后環(huán)境釋放源項(xiàng)的模型如下：

Ai：事故發(fā)生后，核素i的環(huán)境釋放量；

t1：事故發(fā)生0時(shí)刻；

t2：事故發(fā)生后，二回路工質(zhì)停止向環(huán)境釋放時(shí)刻；

C1i：核素i在受影響蒸汽發(fā)生器二次側(cè)蒸汽中的比活度；

C2i：核素i在完好蒸汽發(fā)生器二次側(cè)蒸汽中的比活度；

D2：完好蒸汽發(fā)生器二次側(cè)蒸汽通過(guò)安全閥的蒸汽流量；

C3i：核素i在受影響蒸汽發(fā)生器二次側(cè)水中的比活度；

D3：受影響蒸汽發(fā)生器二次側(cè)滿溢流量；

D1：受影響蒸汽發(fā)生器二次側(cè)到汽輪機(jī)的蒸汽流量；

FP：核素i經(jīng)過(guò)冷凝器的汽水分配系數(shù)；

D4：完好蒸汽發(fā)生器二次側(cè)到汽輪機(jī)的蒸汽流量。

3 計(jì)算參數(shù)

彈棒事故源項(xiàng)分析中用到的參數(shù)包括：裂變產(chǎn)物堆芯積存量、事故造成的燃料元件破損份額、燃料元件氣隙內(nèi)的裂變產(chǎn)物占破損燃料元件積存量的份額、反應(yīng)堆冷卻劑向蒸汽發(fā)生器二次側(cè)的泄漏率等。

（1）小型動(dòng)力堆彈棒事故源項(xiàng)分析所用的參數(shù)包括：

①事故造成的燃料元件破損份額：4%；

②破損燃料元件氣隙內(nèi)的惰性氣體和碘核素占破損燃料元件積存量的份額［14］：10%；

③事故發(fā)生后，反應(yīng)堆冷卻劑向安全殼釋放的釋放速率隨時(shí)間的變化，如圖1所示；

圖1 反應(yīng)堆冷卻劑向安全殼釋放的釋放速率隨時(shí)間的變化Fig.1 Release rate of reactor coolant to containm ent

④事故發(fā)生后，堆芯補(bǔ)水流量隨時(shí)間的變化，如圖2所示；

圖2 堆芯補(bǔ)水流量隨時(shí)間的變化Fig.2 Flow rateof coremakeup

⑤事故發(fā)生后，釋放到安全殼內(nèi)的反應(yīng)堆冷卻劑中的裂變產(chǎn)物釋放到安全殼大氣中的份額［14］：惰性氣體：100%；碘：100%；

⑥釋放到安全殼大氣中的各化學(xué)形態(tài)的碘占總碘的份額［14］：粒子碘：0.95；元素碘：0.048 5；有機(jī)碘：0.0015；

⑦各化學(xué)形態(tài)的碘在安全殼內(nèi)的去除率［13］：粒子碘的去除系數(shù)：0.1 h-1，粒子碘去除的持續(xù)時(shí)間為事故發(fā)生后75.7 h；元素碘的去除系數(shù)：0.11 h-1，元素碘去除的持續(xù)時(shí)間為事故發(fā)生后27.7h；有機(jī)碘無(wú)去除；

⑧事故發(fā)生后安全殼大氣的泄漏率［14］：24 h內(nèi)，0.3%·d-1；24h后，0.15%·d-1；

⑨蒸汽發(fā)生器給水流量隨時(shí)間的變化，如圖3所示；

圖3 蒸汽發(fā)生器給水流量隨時(shí)間的變化Fig.3 Flow rate of steam generator feedwater

⑩事故發(fā)生后39 s，二回路工質(zhì)停止向環(huán)境釋放；二回路工質(zhì)向環(huán)境的釋放流量隨時(shí)間的變化，如圖4所示。

圖4 二回路工質(zhì)向環(huán)境的釋放流量隨時(shí)間的變化Fig.4 Release rate of activity contained in thesecondary side to theenvironment

（2）M 310型核電廠彈棒事故源項(xiàng)分析現(xiàn)實(shí)模型所用的參數(shù)包括：

①假設(shè)燃料包殼破損率為10%，在堆芯中可能熔化的燃料棒估計(jì)為總數(shù)的1%［15］；

②破損燃料元件氣隙內(nèi)的惰性氣體和碘核素占破損燃料元件積存量的份額［15］：

惰性氣體（85Kr除外）：2%；

85K r：30%；

碘：3%；

③對(duì)于堆芯內(nèi)熔化的燃料元件，裂變產(chǎn)物的釋放份額［15］：

惰性氣體：100%；

碘：50%；

④事故發(fā)生后，釋放到安全殼內(nèi)的反應(yīng)堆冷卻劑中的裂變產(chǎn)物釋放到安全殼大氣的份額［15］：惰性氣體：100%；碘：50%；

⑤釋放到安全殼大氣的各化學(xué)形態(tài)的碘占總碘的份額［15］：粒子碘：0.1；元素碘：0.9；

⑥各化學(xué)形態(tài)的碘在安全殼內(nèi)的去除率：事故發(fā)生6 h后手動(dòng)啟動(dòng)安全殼噴淋系統(tǒng)，元素碘的去污因子為1000［16］（由于M 310型核電廠彈棒事故源項(xiàng)分析時(shí)采用的瞬時(shí)釋放模型，對(duì)元素碘考慮瞬時(shí)去除，去污因子為1 000，沒(méi)有去除系數(shù)的概念），粒子碘無(wú)去除；

⑦事故發(fā)生后安全殼大氣的泄漏率［15］：6 h內(nèi)：0.25%·d-1；6h后：0.15%·d-1。

4 計(jì)算結(jié)果

表1給出了小型動(dòng)力堆彈棒事故通過(guò)安全殼泄漏的累積環(huán)境釋放源項(xiàng)；表2給出了以小型動(dòng)力堆彈棒事故發(fā)生到事故發(fā)生后39 s，通過(guò)汽輪機(jī)和蒸汽發(fā)生器安全閥釋放的環(huán)境釋放源項(xiàng)；表3給出了小型動(dòng)力堆彈棒事故后環(huán)境釋放源項(xiàng)。

表1 事故發(fā)生后，通過(guò)安全殼泄漏的累積環(huán)境釋放源項(xiàng)Table 1 The source term release to the environment through containment leakage

表2 通過(guò)二回路釋放途徑的環(huán)境釋放源項(xiàng)Table 2 The source term release to the environment through secondary activity release

表3 事故后環(huán)境釋放源項(xiàng)Table3 Thesource term release to theenvironmentafter accident

5 結(jié)論

綜上所述，小型動(dòng)力堆和M 310型核電廠彈棒事故環(huán)境釋放源項(xiàng)分析方法的主要差異在于：

（1）小型動(dòng)力堆采用分時(shí)間步的方法計(jì)算反應(yīng)堆冷卻劑中放射性核素濃度（考慮核素的衰變、向安全殼釋放造成的消失項(xiàng)以及補(bǔ)水對(duì)冷卻劑中核素濃度的影響）及向安全殼的釋放速率，而M 310型核電廠采用瞬時(shí)釋放模型計(jì)算反應(yīng)堆冷卻劑中放射性核素向安全殼的釋放量；

（2）小型動(dòng)力堆采用的燃料元件破損率為事故分析的結(jié)果，而M 310型核電廠采用的燃料元件破損、熔化份額為RCC-P法規(guī)中的值；

（3）小型動(dòng)力堆采用的以下參數(shù)參考RG 1.183，而M 310型核電廠采用的參數(shù)參考RCC-P法規(guī)：

①破損燃料元件氣隙內(nèi)惰性氣體和碘核素占破損燃料元件積存量的份額；

②釋放到安全殼內(nèi)的反應(yīng)堆冷卻劑中的惰性氣體和碘核素釋放到安全殼大氣的份額；

③釋放到安全殼大氣的各化學(xué)形態(tài)的碘占總碘的份額；

④事故后安全殼大氣的泄漏率。

小型動(dòng)力堆彈棒事故環(huán)境釋放源項(xiàng)分析方法，較M 310型核電廠的瞬時(shí)釋放模型得到更為符合實(shí)際事故進(jìn)程的環(huán)境釋放源項(xiàng)。

［1］楊慶明，汪俊，唐濤.事故運(yùn)行與事故分析的關(guān)系［J］.核安全，2014，13（4）：6-10.

［2］張志剛，范育茂.輻射安全文化的建立與培育［J］.核安全，2010（3）：22-26.

［3］吳斌，賈銘椿，龔軍軍.船用堆核事故狀態(tài)下源項(xiàng)特性及計(jì)算方法研究［J］.海軍工程大學(xué)學(xué)報(bào)，2003，15（5）：87-90.

［4］陳濟(jì)東.大亞灣核電站系統(tǒng)及運(yùn)行［M］.北京：原子能出版社，1994.

［5］俞冀陽(yáng)，俞爾俊.核電廠事故分析［M］.北京：清華大學(xué)出版社，2012.

［6］朱繼洲.核反應(yīng)堆安全分析［M］.西安：西安交通大學(xué)出版社，2000.

［7］林誠(chéng)格.非能動(dòng)安全先進(jìn)壓水堆核電技術(shù)［M］.北京：原子能出版社，2010.

［8］張帆，朱波，邾明亮，等.小型動(dòng)力堆大破口事故下安全殼內(nèi)氣載放射性物質(zhì)濃度計(jì)算［J］.原子能科學(xué)技術(shù)，2010，44：269-273.

［9］J.L.Grover.AP1000 Rod Ejection Accident Doses［R］.CNCRA-02-9.W estinghouse Electric Company LLC，2007.

［10］陳曉秋，岳會(huì)國(guó)，林權(quán)益.雙層安全殼壓水堆LOCA事故的放射性后果分析審查［J］.輻射防護(hù)通訊，2005，25（4）：1-8.

［11］喬清黨，李冰，岳會(huì)國(guó)，等.核事故應(yīng)急氣載放射性物質(zhì)擴(kuò)散的精細(xì)化模擬方法研究［J］.核安全，2014，13（1）：28-33.

［12］A Guide for the TACT IIIComputer Code［R］.NUREG/ CR-3287.Oak RidgeNational Lab，1983.

［13］趙科.GVACT程序使用手冊(cè)［R］.1997，3.

［14］Regulatory Guide 1.183 A lternative Radiological Source Terms for Evaluating Design Basis Accidents at Nuclear Power Reactors s［Z］.U.S.Nuclear Regulatory Comm ission，2000.

［15］RCC-P 900MWe壓水堆核電廠系統(tǒng)設(shè)計(jì)和建造規(guī)則［S］.法國(guó)電力公司和法馬通公司，1991.

［16］IAEA Specific Safety GuideNo.SSG-2［S］.IAEA，2009.

Analysisof theSource Term Released to theEnvironmentin Caseof a Rod Ejection Accident in a SmallNuclear Power Reactor

CHENG Shisi1，WU Han2，*
（1.Scienceand Technologyon Reactor System Design Technology Laboratory，NuclearPower InstituteofChina，Chengdu610041，China；2.Nuclearand Radiation Safety Center，MEP，Beijing100082，China）

Asa design basisaccident（DBA），a rod ejection accident leads to the releaseofactivity to the environment through containment leakage，and the releaseofactivity contained in the secondary side isalso considered.Thispaperaimsatanalyzing the source term released to theenvironmentin caseofa rod ejection accident in a smallpower reactorw ith themethodsof the AP1000.Themethodswere compared w ith thoseof aM 310nuclearpowerplant.Theanalysismethod and the resultof thesourceterm released to theenvironment fora rod ejectionaccidentaregiven in thispaper.

smallpower reactor；rod ejectionaccident；source term

TL732

：A

：1672-5360（2015）04-0058-06

2015-03-27

2015-05-12

程詩(shī)思（1988—），女，四川榮縣人，助理工程師，現(xiàn)主要從事輻射屏蔽與環(huán)境安全分析工作

*通訊作者：吳 晗，E-mail：wuhan@chinansc.cn