基于不確定分析的CPR1000核電廠大破口事故研究

楊 江，羅漢炎，劉仲昊，崔旭陽，陳韻茵

（中國廣核集團中廣核研究院有限公司，廣東深圳 518000）

大破口失水事故（Large Break Loss of Cool?ant Accident，LB LCOA）作為RG 1.70、ANS 18.2和EJ/T 312中劃分的Ⅳ類工況（極限事故），是后果最惡劣的設計基準事故之一。LB LCOA直接制約了核電廠可接受的最高線功率，進而很大程度上限制了核電廠經濟性的提升。

LB LOCA的現象機理復雜、輸入參數繁多，其分析結果受分析方法的影響很大。在傳統的LB LOCA事故分析中，采用了保守的分析方法，即采用保守的軟件模型，以及保守的核電廠狀態參數（初始狀態參數和邊界條件參數），從而得到非常保守的分析結果，所以，傳統的LB LO?CA事故分析中安全裕量均很小。在國內CPR1000核電廠當前執照申請中，LB LOCA關鍵驗收準則包殼峰值溫度（Peak Cladding Tem?perature，PCT）的裕量僅為1%~2%［1］。保守的LB LOCA分析方法不但不利于核電廠的高經濟性，并且有可能錯誤的扭曲事故瞬態和現象。

美國核設施聯邦法規10 CFR 50.46，對輕水堆應急堆芯冷卻系統驗收準則和分析方法做出了規定，即規定了LOCA驗收準則和保守的分析方法。1988年美國核管會（Nuclear Regulatory Commission，NRC）對10 CFR 50.46進行修訂，允許采用最佳估算疊加不確定性方法（Best Esti?mate Plus Uncertainty，BEPU），該方法基于最佳估算程序、同時對結果的不確定性進行量化分析。此后，BEPU方法在國外被廣泛研究和運用［2-4］，比如AP1000核電廠在采用BEPU方法以后，其LB LOCA的PCT的安全裕量為15.4%［5］。

2006年中國核安全局發布的核安全導則HAD 102-17《核動力廠安全評價與驗證》，對LOCA分析方法進行了相關規定：“應用最佳估算程序與初始和邊界條件作現實假設相結合的方法也可接受。該方法應該基于在規定的高可信度條件下相應的考慮了核動力廠工況和程序模型的統計組合的不確定性后以規定的高可信度，使計算結果滿足驗收準則。”可見，目前在我國執照申請中，LOCA事故分析采用BEPU方法也是可以接受的。

本文基于BEPU思想，采用最佳估算程序LOCUST，對CPR1000核電廠LB LOCA進行不確定性分析。首先梳理得到了關鍵的不確定性輸入參數，包括程序模型參數和電廠狀態參數兩方面；然后對這些參數進行59組抽樣，得到59個工況；利用LOCUST進行59個工況的計算；最后對59個工況的計算結果進行統計分析，得到95%置信度下95%概率的PCT（簡稱PCT95%/95%）。結果表明，相對傳統保守分析方法，不確定性分析方法能夠更加真實合理地評估事故瞬態和后果，能識別影響事故后果的關鍵參數，最重要的是，所得到安全裕量有顯著提升。

圖1 保守方法和BEPU方法的保守裕量對比Fig.1 Conservative Margins Comparisons between Conservative method and BEPU method

1 程序與模型

1.1 程序

LOCUST是中廣核研究院有限公司開發的核電廠熱工水力計算軟件，經過了軟件驗證，目前在核電廠安全分析中應用，適用于非破口事故和破口事故的模擬分析。

1.2 程序不確定性分析功能

LOCUST具有不確定性分析功能，可實現以下功能：輸入參數抽樣、多個工況自動計算、輸出參數統計分析等功能。

由用戶給定軟件模型參數的不確定性信息和電廠狀態參數（初始狀態參數和邊界條件參數）的不確定性信息，這些信息包括不確定性參數的數量、每個不確定性參數的不確定性分布范圍和分布形式。

程序對所有不確定性輸入參數進行多組抽樣，每組抽樣結果生成一個計算工況，然后進行多個計算工況的模擬計算。

程序對結果進行自動統計分析，得到95%置信度下95%概率值。

1.3 模型

針對CPR1000的系統結構和運行特點，通過控制體的合理劃分，建立其一回路、部分二回路系統和專設安全設施的LOCUST模型，模型節點如圖2所示。

圖2 CPR1000核電廠的LOCUST節點模型Fig.2 LOCUST nodding model for Nuclear Power Plant CPR1000

2 參數及其不確定性

要量化輸出結果的不確定性，首先需要明確輸入參數的不確定性。LB LOCA計算中輸入參數的不確定性主要來自三個方面：程序模型參數、電廠狀態參數（初始狀態參數和邊界條件參數）、用戶效應，本文主要考慮前兩方面。

實際上，無法明確考慮所有輸入參數的不確定性，因此需要甄選出關鍵參數，縮小輸入參數的不確定性的范圍。對甄選關鍵參數，比較成熟的方法是，采用現象識別和排序表（Phe?nomena Identification and Ranking Table，PIRT），即梳理出對輸出結果影響較大的物理現象，然后識別出與這些現象相關的參數，然后依據對輸出結果的影響程度，對這些參數進行重要度排序，進而甄選出關鍵參數。

對于CPR1000關鍵參數的甄選，主要借鑒了兩方面的研究成果：（1）NEA/CSNI的BEMUSE項目第三和第四階段工作中，得到的傳統三環路壓水堆電廠LB LOCA事故的關鍵現象與參數［6，7］，（2）傳統三環路壓水堆的LB LOCA的PIRT［8，9］。針對PIRT表中重要度為“中”和“高”的現象及參數，最終甄選了47個關鍵參數，這些參數可以分為兩類，第一類是程序模型參數，這類參數一般內置于程序源代碼中，對其修改需要改動程序源碼；第二類是電廠狀態參數，這類參數一般列于程序輸入卡中，對其修改僅需要改動輸入卡。對這些47個參數，進一步確定其不確定性范圍及其分布見表1和表2。對于第一類程序模型參數的不確定性范圍及其分布，主要參考了NEA/CSNI的BEMUSE項目第三和第四階段工作成果［6，7］和LOCUST程序的相關驗證工作［10］；對于第二類電廠參數的不確定性范圍及其分布，主要參考CPR1000的相關設計值和運行值。

表1 核電廠狀態參數Table 1 State parameters of Nuclear Power Plant

表2 程序模型參數Table 2 Model parameters in code

3 抽樣、計算和統計

利用LOCUST的不確定性分析功能，針對47個關鍵參數進行59次隨機抽樣。得到59個計算工況，同時增加1個不考慮參數不確定性的基準工況，然后開展總計60個工況的模擬計算。

運用Wilks方法對60個計算結果進行統計分析。Wilks是一種非參數統計方法，能對未知分布的隨機樣本建立容忍置信區間［11］。

當只希望對單個參數（比如，只針對PCT）進行統計，對于單側容忍區間，Wilks公式為：

式中，β為置信度；

γ為概率份額；

N為計算次數。

因此，需進行59次計算得到59個PCT結果值，其中最大值即為95%置信區間下95%概率值。

圖3給出了由60個工況計算得到的PCT隨事故瞬態時間變化曲線，所有工況PCT呈現出相同的隨時間變化趨勢，且均低于限值驗收準則（1477.15 K）。最大PCT、最小PCT和基準工況PCT見圖4。圖5和圖6中，進一步給出了PCT最大工況、PCT最小工況、基準工況的穩壓器壓力和堆芯坍塌水位，結果顯示，當穩壓器壓力下降更快和堆芯坍塌水位更低時，PCT更大，說明計算的瞬態進程和現象的合理性。

圖3 60個工況的PCTFig.3 PCT of sixty operating conditions

圖4 最大、最小和基準工況的PCTFig.4 Maximum PCT、minimum PCT and PCT of design basis condition

圖5 PCT最大工況、PCT最小工況、基準工況的穩壓器壓力Fig.5 Pressurizer pressure of maximum PCT、minimum PCT and design basis condition

圖6 PCT最大工況、PCT最小工況、基準工況的堆芯坍塌液位Fig.6 Core collapse level of maximum PCT、minimum PCT and design basis condition

59個工況計算得到的最大PCT為1361.3 K，即為95%置信區間下95%概率值，相對于驗收準則（1477.15 K）具有的7.8%的裕量，較CPR1000核電廠傳統保守分析方法的裕量有顯著提升。

對所有PCT進行了柱狀圖統計分析，如圖7所示，其最大概率區間在1180~1200℃，進一步統計分析表明所有這些PCT值分布滿足正態分布。

圖7 PCT的柱形分布Fig.7 Columnar distribution of PCT

4 參數敏感性分析

基于60個工況的輸入參數分布，對PCT進行了輸入參數敏感性分析，得到敏感性因子見如表1和表2的最后一列。結果表明，對PCT“正向”影響程度最高的參數為（按影響大小排序）：塞狀流相間傳熱系數（到氣相）、霧狀流相間傳熱系數（到氣相）、氮氣安注箱初始溫度、燃料棒熱容、氮氣安注箱初始壓力、堆芯初始功率及飽和核沸騰。

5 結論

本文利用LOCUST程序的不確定性分析功能，針對CPR1000核電廠LB LOCA進行分析，得到了PCT的95%置信度下95%概率值為1361.3 K，結果表明，采用基于不確定性分析的方法后，事故安全裕量相對傳統的保守方法有顯著提升。本文對輸入參數對結果的影響進行了敏感性分析，得到了影響程度的排序。