, ,,
(環境保護部核與輻射安全中心,北京 100082)
CAP1400是我國自主開發的具有自主知識產權的三代非能動壓水堆核電技術,在AP1000的基礎上提升了功率,并對反應堆冷卻劑系統、專設安全設施、主要核島輔助系統和主設備以及核島廠房布置等進行了重新設計和系統性的優化[1-2]。因此,針對CAP1400主蒸汽管道破裂事故開展分析研究是有必要的。
主蒸汽管道破裂是設計基準事故中反應堆一次側熱輸出增加的重要工況之一[3],破口流量的上升使得堆芯有重返臨界和重返功率的可能性。本文重點關注這方面對事故的影響。
本文采用SNAP程序進行圖形化模型建立及計算結果后處理,SNAP(Symbolic Nuclear Analysis Package)程序是由NRC資助開發的圖形化分析程序包,功能包括創建和編輯工程分析軟件的輸入、提交計算申請、過程監視以及程序之間的相互調用,可以有效提高事故分析的效率。使用RELAP5/MOD 3.3程序進行瞬態計算,Relap 5是輕水堆冷卻系統事故工況的瞬態行為最佳估算程序[4],其功能幾乎覆蓋了核電廠全部熱工水力工況。RELAP5/MOD 3.3版本是由美國核管會(NRC)于2001年,經過大量的實驗數據和運行經驗反饋,對RELAP5系列程序的升級和改進后發布的,已有廣泛應用[5]。采用VIPRE程序進行子通道熱工水力分析。VIPRE-W程序是由美國西屋公司研制用來考慮水力和核對堆芯及熱通道焓升影響的一個三維子通道程序。
根據事故分析需要,采用RELAP5/MOD 3.3程序,建立了CAP1400反應堆系統模型,示意圖如圖1所示。模型包括兩條環路,每條環路根據冷卻劑流動方向,分別包括一根冷卻劑主管道熱管段,一臺蒸汽發生器,兩臺反應堆冷卻劑泵,兩根冷管段。另外,系統還包括穩壓器、波動管、閥門及用于運行控制和安全觸發的儀表。
主蒸汽管道破裂事故假設主蒸汽管道破裂引起蒸汽從破口排放,破口的蒸汽流量在發生破裂后的短暫時間內迅速增加,而后,隨著蒸汽壓力下降而減小。反應堆冷卻劑系統能量的過多移出導致反應堆冷卻劑溫度和壓力下降。由于堆芯具有負的慢化劑溫度系數,溫度下降將引入正的反應性,從而使得反應堆功率上升[6]。
如果在停堆后有一束具有最大反應性價值的控制棒組件卡在它完全抽出的位置上,這就增加了堆芯重返臨界和重返功率的可能性。主蒸汽管道破裂之后重返功率可能是一個潛在的重要問題。因為假設一束具有最大反應性價值的控制棒組件卡在它完全抽出的位置上,會產生一個高的功率峰因子。本節的目的是分析反應堆功率運行時發生主蒸汽管道破裂事故,以關注在反應堆停堆前或停堆期間堆芯保護是否依然有效。分析從初始熱態滿功率開始至PRHR和CMT投入。
主蒸汽管道破裂事故分析的初始工況采用與安分報告相同的初始數值。
根據不同的主蒸汽管道破口尺寸,事故被劃分為Ⅲ類或Ⅳ類工況,保守考慮,本報告采用Ⅱ類工況的驗收準則,具體包括:最小DNBR高于95/95DNBR限值,以保證燃料包殼完整性;超功率限值條件下,不應發生燃料熔化;RCS和主蒸汽系統壓力不應高于設計壓力的110%;在不發生其他單一失效的情況下,不能導致更嚴重的事故工況。該事故的主要假設包括:
(1)與堆芯相關的假設
1)慢化劑溫度系數:假設慢化劑溫度系數為絕對值最大值;
2)多普勒效應:多普勒反應性系數為絕對值最小值;
3)快中子壽命:取0.000 019 8 s;
4)停堆反應性:假設控制棒組件中負反應性價值最大的一組卡死在堆頂。
(2)控制和保護系統
根據破口尺寸,下列停堆信號將給熱態滿功率主蒸汽管道破裂事故提供必要的保護:
1)超溫ΔT;
2)穩壓器低壓力;
3)S信號:①蒸汽管道低壓力信號;②冷段低溫信號。
(3)初因事件與功能假設。
反應堆功率為100%,0 s時刻,主蒸汽管道破裂,保守假設給水流量與蒸汽流量相匹配。分析時認為喪失廠外電是該事故的潛在后果,不需要詳細分析喪失廠外電工況。
本節分析了CAP1400熱態滿功率主蒸汽管道0.009~0.15 m2的破口譜,如圖2所示。結果表明,結果表明,直到0.058 m2的蒸汽管道破口都不會觸發反應堆停堆。反應堆不停堆工況與負荷過增事故相似,反應堆達到一個新的與蒸汽增加相對應的功率水平。對于0.059~0.105 m2的蒸汽管道破口,反應堆由超功率ΔT信號觸發停堆。對于0.106~0.15 m2的蒸汽管道破口,反應堆由蒸汽管道低壓力安注信號觸發停堆。

圖1 CAP1400一回路系統節點模型Fig.1 The node model for primary coolant system of CAP1400

圖2 破口譜Fig.2 The break spectrum
從DNB和燃料中心熔化(kW/m)保護角度考慮,在分析計算中,超功率ΔT觸發停堆時,破口尺寸越大,達到的停堆功率越大,最極限的工況是破口尺寸為超功率觸發停堆的最大破口尺寸0.105 m2。
針對極限工況進行詳細計算,0 s時刻不含穩壓器回路的主蒸汽管道破裂,引起蒸汽流量迅速增加,事故序列如表1所示。

表1 CAP1400 MSLB事故序列Table 1 Sequence of MSLB accident for CAP1400
將計算所得的關鍵參數進行分析,詳細變化曲線如圖3~圖9所示。
如圖9所示,破口發生后,破口的蒸汽流量在發生破裂后的短暫時間內迅速增加,而后,隨著蒸汽壓力下降而減小。蒸汽流量的增加使得二回路帶熱增加,一回路冷卻劑溫度下降,受慢化劑溫度系數和多普勒功率系數共同影響,堆芯引入正反應性,功率和熱流密度不斷上升,如圖3、圖4所示。功率上升至115%超功率停堆整定值后,超功率ΔT信號觸發停堆,反應堆功率迅速下降。
由于存在慢化劑溫度系數和多普勒功率系數,控制棒插入、功率驟降的過程中會引入較大的正反應性,使得功率在停堆后存在略微反彈的現象,如圖3。停堆后功率約降至40%,分析是由于RELAP5為國際通用的最佳估算系統分析程序,破口流量使用Henry-Fauske模型進行計算,上游出現氣泡即認為進入兩相臨界流,計算的破口流量波動較小,使得停堆后得到的核功率較高。
堆芯熱流密度變化趨勢同核功率變化趨勢類似,如圖4所示。在破口發生后逐漸上升,停堆后快速降低,由于存在較大的慢化劑溫度系數,停堆后核功率較大,熱流密度較高。

圖3 瞬態歸一化核功率Fig.3 Nuclear power vs time

圖4 堆芯熱流密度Fig.4 Heat flus vs time
破口發生后,二回路帶熱增加,導致一回路冷卻劑系統壓力不斷下降,如圖5所示。引入正反應性后,反應堆功率上升,產熱增加,至停堆前,壓力下降趨勢逐漸減弱消失,壓力出現微弱回升。停堆后,由于堆芯產熱驟降,壓力也隨之快速下降。圖6為穩壓器和波動管水體積變化的對比曲線。穩壓器和波動管水體積變化趨勢同穩壓器壓力變化趨勢相同。為關注在反應堆停堆前或停堆期間堆芯保護是否依然有效,本文所計算的主蒸汽管道破裂事故從初始熱態滿功率開始,至PRHR和CMT投入結束。其中穩壓器低-2水位觸發CMT動作的極限整定值為26 m3,瞬態計算時間為23.6 s。

圖5 穩壓器壓力Fig.5 RCS pressure vs time

圖6 穩壓器和波動管水體積Fig.6 RCS water volume vs time
圖7、圖8為壓力容器入口溫度和蒸汽發生器壓力變化的曲線。停堆前,由于破口導致的二回路帶熱增加,壓力容器入口溫度和蒸汽發生器壓力不斷降低。停堆后,堆芯產熱和二回路流量都快速下降。二回路流量下降如圖9。二回路流量降低造成的帶熱能力下降同停堆造成的堆芯產熱降低相互作用,但由于二回路流量下降幅度較小,使得停堆后,二回路帶熱能力較強,故壓力容器入口溫度和蒸汽發生器壓力除在剛停堆后出現微弱的回升外,隨后仍繼續下降。

圖7 壓力容器入口溫度Fig.7 RCS temperature vs time

圖8 蒸汽發生器壓力Fig.8 Steam generator pressure vs time

圖9 蒸汽流量Fig.9 Steam mass flow vs time
根據系統程序的計算結果選取極限點,將極限點的狀態參數(如堆芯流量、熱流密度、溫度、壓力等)提供給物理方向專業人員,采用三維堆芯核設計程序ANC計算該極限點對應的不同功率分布。計算分為4個循環,每個循環根據功率向上、向下傾斜情況不同又分為6個工況。針對這24個不同的功率分布,采用同安分報告相同的F對DNBR進行計算,得到的最小DNBR 1.914即為事故對應的最小DNBR。滿功率MSLB事故DNBR的驗收準則為1.5, DNBR滿足驗收準則。
本文使用最佳估算程序RELAP5/mod3.3為CAP1400核電廠熱態滿功率主蒸汽管道破裂事故建立穩態和瞬態模型,參照安分報告中的基本假設,開展破口譜分析,并就關鍵熱工水力參數進行了分析。
本文分析了0.009~0.15m2的蒸汽管道破口譜。結果表明直到0.058 m2的蒸汽管道破口都不會觸發反應堆停堆。對于0.059~0.105 m2的蒸汽管道破口,反應堆由超功率ΔT觸發停堆。對于0.106~0.15 m2的蒸汽管道破口,反應堆由蒸汽管道低壓力安注信號觸發停堆。從燃料中心熔化保護角度考慮,最極限的工況是0.105 m2。分析表明,對于最極限的工況仍滿足壓力和DNBR的驗收準則。結論認為,CAP1400反應堆在主蒸汽管道破裂事故下安全可靠。