林支康，殷煜皓，梁國興

（上海交通大學(xué)核科學(xué)與工程學(xué)院，上海 200240）

與傳統(tǒng)的二代核電技術(shù)壓水堆不同，AP1000核能系統(tǒng)采用了第三代核電技術(shù)的非能動冷卻系統(tǒng)，可以在設(shè)計基準事故下提供應(yīng)急堆芯冷卻，堆芯的安全性將大大提高。采用非能動冷卻技術(shù)，一旦核電廠在發(fā)生失水事故（LOCA），系統(tǒng)的堆芯補水箱（CMT）、蓄壓安注箱（ACC）、安全殼內(nèi)置換料水箱（IRWST）以及安全殼再循環(huán)水都可以立即向反應(yīng)堆冷卻系統(tǒng)注水。特別是設(shè)計中采用了自動降壓系統(tǒng)（ADS）和非能動余熱排出系統(tǒng)（PRHRHX），一旦發(fā)生小破口失水事故（SBLOCA），隨著反應(yīng)堆冷卻劑系統(tǒng)從破口和ADS降壓，冷卻劑可以通過堆芯補水箱和蓄壓安注箱添加到反應(yīng)堆壓力容器中；當(dāng)系統(tǒng)降壓到IRWST系統(tǒng)注入壓力時，IRWST中的水將不斷地注入反應(yīng)堆堆芯，讓堆芯保持在長期冷卻狀態(tài)。由于自動降壓系統(tǒng)可以實現(xiàn)反應(yīng)堆冷卻系統(tǒng)的快速降壓，因此IRWST可以依靠自身的重力注水，對堆芯實現(xiàn)長期的冷卻。本文將采用分析程序RELAP5［1］建立AP1000的失水事故分析模式，并用這個模式對AP1000發(fā)生小破口失水事故后系統(tǒng)的特性進行分析，分析破口面積和功率提升對事故嚴重性的影響。

1 模式的建立

本文利用RELAP5最佳估算程式，對AP1000核電廠發(fā)生SBLOCA后的運行建立模型，目的是分析系統(tǒng)中壓力容器（RPV）、蒸汽發(fā)生器（SG）、反應(yīng)堆冷卻劑泵（RCP）、穩(wěn)壓器及主管道等重要部件的核蒸汽供應(yīng)系統(tǒng)（NSSS）；為了方便分析，二次側(cè)與安全殼則進行了必要的簡化。同時，為了模擬SB－LOCA，還根據(jù)US Westinghouse公司的AP1000設(shè)計文件DCD［2］對非能動安全系統(tǒng)，如非能動堆芯冷卻系統(tǒng)（PXS）、非能動余熱導(dǎo)出換熱器（PRHR－HX）及反應(yīng)堆停堆保護系統(tǒng)進行了模擬；重要建模數(shù)據(jù)包括幾何數(shù)據(jù)、控制邏輯、初始條件以及邊界條件等。

堆芯的節(jié)點劃分對包殼峰值溫度影響比較大，因此對堆芯區(qū)域的水力部件和熱構(gòu)件進行了詳細的劃分，堆芯水力部件被分成熱通道（hot channel）和平均通道（average channel）兩部分，并且在兩者之間加橫向通道（cross flow junction），使之能相互流動，其中平均通道和熱通道軸向均劃分為10個節(jié)點。在熱構(gòu)件中，堆芯被分為三部分：最熱棒（hot rod），最熱棒束（hot bundle）和平均棒束（average bundle）。前兩者對應(yīng)的水力部件是熱通道，后者對應(yīng)的水力部件是平均通道。AP1000核電廠RELAP5節(jié)點劃分如圖1所示。

2 模式的校驗

模式的校驗，采用US Westinghouse公司NOTRUMP［3］分析模式的小破口（SBLOCA）事故分析結(jié)果進行對比。為了讓對比結(jié)果更準確，在SBLO－CA發(fā)生前先要進行相應(yīng)的穩(wěn)態(tài)調(diào)試，確保計算條件穩(wěn)定并且已達到期望的初始條件。在計算條件穩(wěn)定后，開始進行動態(tài)分析，對事故現(xiàn)象以及非能動安全系統(tǒng)的行為進行研究。在分析時，電廠初始狀態(tài)的設(shè)定值與西屋公司NOTRUMP程式采用的包絡(luò)電廠運行狀態(tài)設(shè)置一致。分析涉及的重要狀態(tài)參數(shù)，包括主泵冷卻劑流量、反應(yīng)堆冷卻劑平均溫度（Tavg）、功率分布（包括功率峰值因子FQ和焓升熱管因子FΔH）、穩(wěn)壓器壓力和穩(wěn)壓器水位等。經(jīng)過200s穩(wěn)態(tài)運行后，模式進入了預(yù)期的狀態(tài)，運行穩(wěn)定值的校驗結(jié)果見表1；破口發(fā)生后的主要參數(shù)變化與Westinghouse公司NOTRUMP程式分析結(jié)果的對比結(jié)果見圖2至圖7。

表1 穩(wěn)態(tài)運行值與NOTRUMP程序運行值比較表

圖1 AP1000—回路模型的節(jié)點

圖2 SBLOCA包殼峰值溫度

圖3 SBLOCA穩(wěn)壓器的壓力變化

圖4 CMT注水流量

圖5 ACC注水流量

圖6 IRWST注水流量

圖7 系統(tǒng)水和壓力容器內(nèi)水裝量變化

可以看出，本文采用的模式所進行的分析結(jié)果，在時間和數(shù)量上與西屋公司模擬出的結(jié)果相當(dāng)一致；由于DCD中并未有SBLOCA PCT的計算趨勢圖，而僅記有PCT之峰值為1 016.5K，該值和圖2RELAP5所計算之PCT峰值相當(dāng)接近。

3 模式計算結(jié)果的分析

利用模式計算結(jié)果可以分析AP1000發(fā)生SBLOCA后，系統(tǒng)特性受到的影響，也可以對破口大小對事故嚴重性的影響進行定量的分析，還可以對事故發(fā)生后功率提升與包殼的溫度的關(guān)系進行定量分析。

3.1 設(shè)計特性分析

AP1000采用了與傳統(tǒng)PWR不同的特殊設(shè)計。采用這種設(shè)計，一旦發(fā)生SB－LOCA，由于一回路系統(tǒng)的壓力下降相對較慢，特殊設(shè)計的ADS自動降壓系統(tǒng)可以保證系統(tǒng)泄壓和確保安注系統(tǒng)重力補水動作，將小破口失水事故下的一回路系統(tǒng)壓力降到IRWST注水壓力，從而使IRWST能夠依靠自身重力向堆芯提供長期冷卻水。

在AP1000的系統(tǒng)設(shè)計中，自動降壓系統(tǒng)（ADS）共分4級，由連接在穩(wěn)壓器和熱管段上的一系列閥門組成，向反應(yīng)堆冷卻劑系統(tǒng)提供分階段的降壓，其中第1到第3級連接在穩(wěn)壓器頂端，并與IRWST相連接，排氣經(jīng)各自的卸壓噴頭［4］（鼓泡器）被IRWST中的水冷凝。ADS的第4級連接在熱管段上，卸壓閥采用結(jié)構(gòu)十分簡單，控制和動作原理也采用簡單的爆破閥。由于ADS的閥口徑較大，因此可以從RCS壓力邊界直接向安全殼的環(huán)路隔間卸壓。

AP1000采用了創(chuàng)新的自動降壓系統(tǒng)，利用RELAP5LOCA分析模式對自動降壓系統(tǒng)對小破口失水事故的影響進行分析。分析結(jié)果表明：當(dāng)破口直徑為7.62cm時，如果4級ADS都失效不能開啟，那么一回路系統(tǒng)的壓力下降就會變慢，ACC和CMT的注水時間推遲，雖然出現(xiàn)了7.62cm破口，但是這樣大小破口帶來的泄壓并不能使一回路系統(tǒng)壓力下降到IRWST的注水壓力，這樣IRWST不能保證向堆芯提供長期冷卻水。計算結(jié)果還表明，在事故發(fā)生后17 500s內(nèi)，壓力容器中的水位無法有效覆蓋燃料，燃料包殼溫度將迅速升高。ADS正常動作與失效時的包殼峰值溫度比較，結(jié)果如圖8?？梢钥吹?，ADS的失效導(dǎo)致的系統(tǒng)壓力下降明顯要慢，導(dǎo)致ACC注水和CMT注水往后推遲（圖9）。由圖10可以看出，ADS失效后，由于IRWST沒有提供長期的冷卻水，使壓力容器內(nèi)的水位比ADS正常時的水位要低很多，導(dǎo)致堆芯裸露。分析結(jié)果進一步顯示，至少需維持1—3級ADS正常功能，IRWST才能在SBLOCA中提供有效重力補水（圖11）。

圖8 ADS失效與正常時PCT的比較

圖9 ADS失效與正常時穩(wěn)壓器壓力的比較

圖10 壓力容器內(nèi)水位變化圖

圖11 IRWST在ADS維持1—3級的注水流量

3.2 破口大小對事故嚴重性的影響

在AP1000堆型的安全設(shè)計中，面積小于0.093m2破口稱為小破口［5］。在小破口失水事故中，破口的大小會影響一回路系統(tǒng)的泄壓速率和冷卻劑的噴射速率，進而影響到安注系統(tǒng)的動作，特別是影響到ACC和IRWST的注水時間，進而影響到壓力容器內(nèi)的水量，影響到包殼的溫度，甚至燒毀燃料包殼。

包殼的溫度可以作為事故嚴重程度的度量，利用RELAP5模式對破口大小事故的嚴重性進行分析，破口的大小對包殼的溫度（PCT）的影響，分析結(jié)果如圖12所示。

圖12 包殼峰值溫度隨破口面積的變化

從圖12中可以看出，在小破口失水事故工況下，包殼峰值溫度隨著破口變大而升高，并且破口面積每增加122.6cm2時（5%管道面積），PCT溫度上升約36K。

3.3 事故發(fā)生后功率提升與包殼溫度的關(guān)系

AP1000堆型發(fā)生小破口失水事故后，功率提升也會造成PCT上升。假定AP1000堆型發(fā)生破口當(dāng)量直徑為25.4cm的失水事故，利用模式對功率提升與PCT的關(guān)系進行分析，結(jié)果如圖13所示。

圖13 包殼峰值溫度隨功率提升的變化關(guān)系

從圖13中可以看出，SBLOCA發(fā)生后包殼峰值溫度與功率提升，接近線性關(guān)系，功率每增加5%，PCT上升約42.26K。

4 結(jié)論

本文根據(jù)西屋公司的設(shè)計文件建立了AP1000核電廠的RELAP5SB－LOCA分析模式，通過與西屋公司的NOTRUMP分析模式的對比，驗證了采用RELAP5分析模式模擬AP1000核電廠小破口事故的準確性。

利用RELAP5SB－LOCA分析模式對AP1000的系統(tǒng)設(shè)計特性進行分析，結(jié)果發(fā)現(xiàn)，當(dāng)AP1000核電廠發(fā)生SBLOCA，并且破口直徑在7.62cm時，自動降壓系統(tǒng)（ADS）至少要保證1—3級維持正常動作。否則，一回路系統(tǒng)中的壓力下降速率較慢，不能下降到IRWST的注水壓力，在17 500s之內(nèi)堆芯中的水位無法有效覆蓋燃料。

計算結(jié)果還表明，在小破口失水事故后，破口面積越大，包殼的峰值溫度會越高，破口面積每增加5%的管道面積，PCT溫度上升約36K；當(dāng)量直徑為25.4cm時，每增加5%功率，PCT上升約42.26K。

［1］RELAP5／MOD3.3，2001［R］／／Code Manual Volume 1：Code Structure，System Models，and Solution Methods NUREG／CR－5535.

［2］Westinghouse Electric Company LLC.AP1000TMDesign Control Document，Revision 17［R］／／US：Westinghouse Electric Company，2008：Tier 2Material 15.6－29－15.6－263

［3］Lewis D.Thorson.NOTRUMP，An Updated Version of TRUMP［R］／／.U.S.：USNRC UCID－18682.1980－9－29

［4］K.B.Welter and S.M.Bajorek.APEX－AP1000Confirmatory Testing To Support AP1000Design Certification（nonproprietary）［C］／／U.S.Nuclear Regulatory Commission.August 2005

［5］Accident Analysis For Nuclear Power Plants With Pressurized Water Reactors.Safety Reports Series NO.30［R］／／International Atomic Energy Agency，Vienna，2003