黃 雄，魏 巍，馬國(guó)揚(yáng)，謝明亮，冉曉隆，魏 興

（1. 中核武漢核電運(yùn)行技術(shù)股份有限公司，湖北 武漢 430223；2. 福建福清核電有限公司，福建 福清 350318）

2011 年日本福島核電廠數(shù)臺(tái)機(jī)組發(fā)生氫氣爆炸，地震和海嘯導(dǎo)致堆芯和乏燃料池的冷卻功能喪失，燃料組件不斷升溫，髙溫下鋯合金包殼和其他結(jié)構(gòu)材料與水蒸氣反應(yīng)產(chǎn)生大量氫氣，最后發(fā)生了氫氣爆炸[1]。目前，三代非能動(dòng)核電廠主要采用非能動(dòng)氫氣復(fù)合器和點(diǎn)火器共同作用的方式來消除安全殼內(nèi)的氫氣[2]，從而緩解安全殼內(nèi)的氫氣風(fēng)險(xiǎn)。

為降低氫氣燃燒造成對(duì)安全殼完整性的威脅，研究人員提出了嚴(yán)重事故情況下事故惰化氫氣風(fēng)險(xiǎn)緩解措施[3]，它是通過向大型安全殼內(nèi)注入惰性氣體（氮?dú)饣蚨趸迹?，控制安全殼?nèi)混合易燃?xì)怏w成分，避免達(dá)到可燃濃度，如事故預(yù)惰化、事故后惰化、事故后稀釋措施。

本文基于一體化嚴(yán)重事故分析MAAP5 程序，建立三代非能動(dòng)核電廠安全殼模型，對(duì)嚴(yán)重事故下事故后惰化氫氣緩解措施的有效性進(jìn)行分析，同時(shí)對(duì)事故后惰化的注入位置進(jìn)行分析，為大型安全殼事故后惰化氫氣緩解措施在核電廠氫氣風(fēng)險(xiǎn)控制領(lǐng)域的推廣提供一定的參考。

1 仿真建模

本文將安全殼、主控室、輔助廠房、大氣環(huán)境等劃分為不同的區(qū)域進(jìn)行模擬，如圖1 所示。安全殼內(nèi)劃分為12 個(gè)節(jié)點(diǎn)（環(huán)境編號(hào)13），其中，9 個(gè)節(jié)點(diǎn)用于模擬安全殼內(nèi)壓力邊界，3個(gè)節(jié)點(diǎn)用于模擬殼外即鋼制安全殼與屏蔽構(gòu)筑物之間環(huán)域。同時(shí)，設(shè)置48 個(gè)分布式熱阱，15個(gè)集中式熱阱，30 個(gè)流動(dòng)通道。

1 號(hào)節(jié)點(diǎn)和2 號(hào)節(jié)點(diǎn)分別模擬關(guān)鍵設(shè)備蒸汽發(fā)生器的兩個(gè)隔間，6～9 號(hào)節(jié)點(diǎn)分別用于模擬安全殼上部、堆芯補(bǔ)水箱（CMT）隔間、化學(xué)和容積控制系統(tǒng)隔間及兩個(gè)非能動(dòng)注入系統(tǒng)；4 號(hào)和5 號(hào)節(jié)點(diǎn)分別模擬堆腔室及IRWST 水箱隔間。

非能動(dòng)安全殼冷卻PCS 系統(tǒng)，由10-13 節(jié)點(diǎn)通過流道的連接來模擬，形成安全殼的非能動(dòng)空冷通道；同時(shí)，安全殼頂部噴淋則是通過一組用戶自定義事件來實(shí)現(xiàn)噴淋流量的控制。

2 計(jì)算分析

核電廠典型嚴(yán)重事故導(dǎo)致堆芯損壞的始發(fā)事件主要有：全廠斷電（SBO），冷卻劑喪失事故（LOCA），蒸汽發(fā)生器傳熱管破裂事故（SGTR），未能緊急停堆的預(yù)計(jì)瞬變（ATWS），根據(jù)NUREG-1150 和文獻(xiàn)［4］的分析結(jié)果顯示LOCA 對(duì)堆芯損壞貢獻(xiàn)較高。

本文選取三代非能動(dòng)核電廠一回路主管道發(fā)生雙端斷裂疊加非能動(dòng)余熱排出（PRHR）失效、安全殼內(nèi)換料水箱（IRWST）重力注射失效作為始發(fā)事件。

2.1 事故進(jìn)程分析

事故進(jìn)程如表1 所示，事故發(fā)生初始時(shí)刻100 s，一回路主冷卻劑通過管道破口快速釋放到安全殼中。隨之堆芯水位快速下降，主系統(tǒng)壓力也快速下降。在101.4 s 時(shí)刻，壓力下降到穩(wěn)壓器停堆信號(hào)的整定值，反應(yīng)堆停堆，105.6 s時(shí)刻主泵停運(yùn)。當(dāng)主主系統(tǒng)壓為下降到堆芯補(bǔ)水箱CMT 開始投入的整定值，CMT 啟動(dòng)。當(dāng)CMT 的液位達(dá)到“低 1”時(shí)，觸發(fā)第一級(jí)自動(dòng)泄壓系統(tǒng) ADS1 開啟，第二級(jí)自動(dòng)泄壓系統(tǒng)ADS2 和第三級(jí)自動(dòng)泄壓系統(tǒng)ADS3 則根據(jù)時(shí)間延遲信號(hào)自動(dòng)投入；CMT 和安注箱ACC 的全部注入無法緩解事故的發(fā)展。隨著主回路冷卻劑不斷喪失，堆芯出現(xiàn)裸露，燃料元件溫度不斷上升，鋯合金包殼被氧化，在6 250.0 s 時(shí)下封頭開始形成熔池，持續(xù)1 869 s 后堆芯熔融物大量聚集在壓力容器下封頭，整個(gè)事故進(jìn)程中鋯合金與水發(fā)生劇烈的氧化反應(yīng)，短期內(nèi)產(chǎn)生大量的氫氣。

表1 事故進(jìn)程Table 1 The accident process

2.2 氫氣源項(xiàng)分析

圖2 和圖3 分別為破口處的氫氣產(chǎn)量及釋放速率。從圖可以看出，整個(gè)事故進(jìn)程中隨著堆芯冷卻劑的減少，溫度迅速升高，金屬（Zr、Fe、Cr）與水蒸氣在高溫下的反應(yīng)、水的福射分解和侵蝕反應(yīng)等產(chǎn)生的大量氫氣從破口釋放出來，堆芯鋯合金與水反應(yīng)氫氣累計(jì)產(chǎn)量為266 kg，釋放速率主要集中在4 000～5 000 s 時(shí)間段內(nèi)，最高速率可達(dá)到0.55 kg/s。此后，由于鋯包殼外側(cè)形成的致密氧化鋯膜抑制了鋯合金與水反應(yīng)，或由于堆芯熔融物堵塞水或蒸汽流道等原因，減少了鋯合金與水的接觸面積，致使氫氣的產(chǎn)生速率急劇下降，隨著冷卻劑喪失最終氫氣產(chǎn)生速率趨近于0.0。

2.3 安全殼內(nèi)隔間氫氣濃度分析

圖4 為安全殼內(nèi)不同隔間的氫氣濃度隨時(shí)間變化曲線。從圖可看出，破口釋放出來的氫氣起初大量聚集在1 號(hào)SG 源項(xiàng)隔間內(nèi)部，4 885 s時(shí)刻氫氣濃度峰值為0.072 16，隨后在浮升力等的作用下快速擴(kuò)散至安全殼上部空間及周圍隔間，最終大量的氫氣會(huì)聚集在安全殼上部隔間。本文選取安全殼內(nèi)氫氣風(fēng)險(xiǎn)相對(duì)較高的隔間為1 號(hào)SG 源項(xiàng)隔間和安全殼上部隔間作為研究對(duì)象，對(duì)事故后惰化緩解措施進(jìn)行分析。

3 惰化氣體選擇

核電廠安全殼內(nèi)事故后惰化氣體主要是考慮N2、CO2、Ne 和He 等惰性氣體雖化學(xué)穩(wěn)定性比較好，但難以液化且成本較高，故在核電廠不適合做惰化氣體。其他常用的滅火氣體（比如Halon130）會(huì)產(chǎn)生難處理的腐蝕和放射性物質(zhì)，造成環(huán)境污染，不適合做惰化氣體。

同時(shí)，事故后若從常溫常壓下開始惰化，達(dá)到同樣的惰化效果時(shí)，采用CO2的質(zhì)量小于采用的N2質(zhì)量；同時(shí)，采用CO2惰化后的壓力也小于采用的N2惰化后的壓力，相比之下，CO2的使用降低了對(duì)安全殼通風(fēng)的要求。N2和CO2的物理性質(zhì)于表2 所示。

表2 CO2 和N2 的物理特性Table 2 Physical properties of CO2 and N2

事故后惰化氣體注入安全殼方式有液態(tài)注入和氣態(tài)注入。核電廠在滿足向安全殼內(nèi)注入惰化氣體的條件下不存在技術(shù)方面的問題，也不會(huì)影響核電廠的正常運(yùn)行。若采用液態(tài)注入安全殼，在防止管道滲透方面存在著技術(shù)方面的問題，且氣體結(jié)晶可能對(duì)安全殼內(nèi)的設(shè)備部件等安全性造成影響。

因此，本文選用CO2氣體作為惰化劑，其注入時(shí)刻為5 000 s，溫度300 K，注入量350 t，CO2以氣態(tài)的方式注入安全殼內(nèi)。整個(gè)事故進(jìn)程期間不考慮氫氣點(diǎn)火器和非能動(dòng)氫氣復(fù)合器的消氫作用，主要假想三種工況進(jìn)行對(duì)比分析：工況 0——不采取是事故惰化氫氣緩解措施，工況 1——事故后期惰化氣體CO2從安全殼上部隔間注入，工況 2——事故后期惰化氣體CO2從1 號(hào)SG 源項(xiàng)隔間注入。

4 結(jié)果分析

4.1 安全殼壓力分析

圖 5 為安全殼壓力隨時(shí)間變化曲線。從圖可以看出，在5 000 s 時(shí)刻，向安全殼內(nèi)注入惰化氣體CO2，安全殼的壓力基本處于線性增長(zhǎng)的狀態(tài)，計(jì)算結(jié)束前安全殼壓力達(dá)到了5.047 × 105Pa，超過了安全殼的設(shè)計(jì)壓力4.07 × 105Pa。事故后大量的惰化氣體注入安全殼會(huì)嚴(yán)重威脅到核電廠的最后一道安全屏障，則在采取事故后惰化措施來緩解安全殼內(nèi)的氫氣風(fēng)險(xiǎn)時(shí)，安全殼抽氣系統(tǒng)在確保安全殼完整性方面是一個(gè)必不可少安全設(shè)施。

4.2 安全殼內(nèi)氫氣濃度分析

圖6 和圖7 分別示出了1 號(hào)SG 源項(xiàng)隔間和安全殼上部隔間的氫氣濃度隨時(shí)間的變化圖。從圖可以看出，在5 000 s 時(shí)刻向安全殼內(nèi)注入惰化氣體CO2后，1 號(hào)SG 源項(xiàng)隔間和安全殼上部隔間中，工況 0 的平均氫氣濃度明顯高于工況 1 和工況 2。同時(shí)，在1 號(hào)SG 源項(xiàng)隔間，工況 3 的平均氫氣濃度低于工況 2 的平均氫氣濃度，而在安全殼上部隔間，工況 3的平均氫氣濃度與工況 2 的平均氫氣濃度基本一致，這是由于CO2的注入位置在安全殼上部隔間，一定時(shí)間內(nèi)，CO2在安全殼上層空間大量聚集后，在重力、擴(kuò)散力等作用下迅速向安全殼下部空間遷移，這無疑抑制了 1號(hào) SG 源項(xiàng)隔間氫氣擴(kuò)散到安全殼上部隔間的速度。因此，對(duì)于1 號(hào)SG 源項(xiàng)隔間，惰化氣體CO2在源項(xiàng)隔間注入，惰化效果明顯比安全殼上部隔間注入好。

4.3 氫氣風(fēng)險(xiǎn)分析

圖8 和圖9 分別為1 號(hào)SG 源項(xiàng)隔間和安全殼上部空間的氫氣風(fēng)險(xiǎn)圖，用夏皮羅圖進(jìn)行表征，夏皮羅圖根據(jù)氫氣、空氣和水蒸氣的比例關(guān)系將安全殼內(nèi)的狀態(tài)點(diǎn)劃分為不可燃區(qū)、燃燒區(qū)、燃爆轉(zhuǎn)化區(qū)和爆炸區(qū)。在該工況下，1號(hào)SG 源項(xiàng)隔間和安全殼上部隔間，工況 0 中大部分的狀態(tài)點(diǎn)都處于氫氣燃燒區(qū)，均存在氫氣燃燒的風(fēng)險(xiǎn)。采取事故后惰化注入CO2氫氣緩解措施，1 號(hào)SG 源項(xiàng)隔間和安全殼上部隔間大部分的狀態(tài)點(diǎn)向不可燃區(qū)轉(zhuǎn)移，氫氣燃燒風(fēng)險(xiǎn)的概率下降。同時(shí)，在源項(xiàng)隔間注入惰化氣體CO2，降低源項(xiàng)隔間氫氣燃燒風(fēng)險(xiǎn)的效果更明顯。

5 結(jié)論

本文基于一體化嚴(yán)重事故分析MAAP5 程序，建立三代非能動(dòng)核電廠的安全殼模型，選取核電廠一回路主管道發(fā)生雙端斷裂疊加PRHR 失效、IRWST 重力注射失效作為始發(fā)事件，對(duì)事故進(jìn)程、氫氣源項(xiàng)、事故后惰化氫氣緩解措施有效性和事故后惰化注入位置進(jìn)行研究，主要結(jié)論如下：

（1） 整個(gè)事故進(jìn)程中，堆芯鋯合金與水反應(yīng)氫氣累計(jì)產(chǎn)量為266 kg，氫氣釋放速率峰值可達(dá)到0.55 kg/s；

（2） 事故后惰化注入 CO2氫氣緩解措施，可以有效降低安全殼內(nèi)氫氣燃燒風(fēng)險(xiǎn)，且對(duì)1 號(hào)SG 源項(xiàng)隔間，源項(xiàng)隔間注入惰化氣體CO2比安全殼上部注入等量CO2的效果更明顯。

（3） 采取事故后惰化措施來緩解安全殼內(nèi)的氫氣風(fēng)險(xiǎn)時(shí)，安全殼抽氣系統(tǒng)在確保安全殼完整性方面是一個(gè)必不可少安全設(shè)施。