COSINE最佳估算大LOCA評價模型評估矩陣開發

傅孝良，劉麗芳，于楠，杜爭，何斯琪，董博，梁國興，楊燕華

(1.國家核電技術有限公司北京軟件技術中心，北京　102209; 2.上海交通大學核科學與工程學院，上海　200240)

核電技術

COSINE最佳估算大LOCA評價模型評估矩陣開發

傅孝良1，劉麗芳1，于楠1，杜爭1，何斯琪1，董博1，梁國興2，楊燕華1

(1.國家核電技術有限公司北京軟件技術中心，北京102209; 2.上海交通大學核科學與工程學院，上海200240)

依據美國NRC最新的EMDAP方法，基于壓水堆大破口LOCA事故發展特征，識別各發展階段的重要現象和過程，并以此為基礎，結合分析國內外已有實驗數據，開發了3個與大破口LOCA事故發展階段相對應的最佳估算評價模型評估矩陣。

最佳估算;大破口失水事故;評估矩陣

美國10 CFR 50.46［1］要求:用于核電廠安全分析的評價模型必須有足夠的數據支持以證明分析結果真實地描述了核反應堆系統在LOCA事故中的行為。必須將評價模型計算結果與可用的實驗數據相比較，且需要識別并評估安全分析方法和數據參數的不確定性，以確保可以評估計算結果的不確定度。否則，分析應急堆芯冷卻系統(ECCS)行為的評價模型必須滿足10 CFR 50附錄K［2］(簡稱附錄K)的要求和接受準則。對于評價模型的開發，美國核管會(NRC)先后發布了RG1.157［3］和RG1.203［4］兩個導則，并在RG1.203中提出了評價模型開發與評估過程方法——EMDAP(Evaluation Model Development and Assessment Process)。根據該方法，開發評價模型的首要工作是識別評價模型的應用目的，開發與重要現象和過程相對應的模型評估矩陣。

針對我國首個自主開發的堆芯設計和系統分析一體化軟件包——COSINE軟件包中的系統分析程序最佳估算模型，選定其應用范圍為AP1000/CAP1400系列非能動壓水堆(同時兼顧傳統二代3/4環路壓水堆)的大破口LOCA事故分析［5］，識別出事故序列中的重要現象和過程，并結合國內外已有核電實驗數據的研究，開發與大破口LOCA事故進程相對應的三個模型評估矩陣。

1　評價模型開發與評估方法

美國NRC 2005年發布的RG1.203［4］中為評價模型的開發與評估推薦了一種最新的分解方法——EMDAP方法(見圖1)。

EMDAP方法主要包含了以下4項任務(共20個步驟):

(1)建立評價模型能力需求。

定義評價模型的分析能力范圍，如針對哪類具體的堆型和哪類具體的事故序列，據此定義相關的性能指標;并對該堆型對應事故序列下的重要物理過程和現象進行識別和排序，得出一個在評價模型開發與評估過程中發揮重要作用的現象識別與排序表(PIRT表)。

(2)開發評估數據庫。

根據任務(1)所建立的PIRT表，收集或開展與排序高分的現象相關的實驗，并評價這些實驗數據的失真度和不確定度。

(3)開發評價模型。

開發用于分析對應堆型對應事故序列的評價模型(包括輸入輸出處理程序、初始條件邊界條件假設和模型選項選擇等信息)。

(4)評估評價模型的適宜性。

利用任務(2)所收集的實驗數據，對任務(3)所開發的評價模型進行適宜性評估，包括利用分離效應實驗數據進行的閉合關系式評估，和利用整體效應實驗數據進行的評價模型整體評估，并根據評估結果確定評價模型計算結果的偏差和不確定度。

因此，應用于核電廠安全分析的評價模型必須經過大量的實驗數據確認與評估。根據最新的評價模型開發與評估方法，開發與應用目的相對應的現象識別與排序表及模型評估矩陣，是整個評價模型開發與評估工作的基礎。

2　COSINE最佳估算大LOCA評價模型能力需求分析

根據EMDAP方法，評價模型開發與評估的第一項任務是確定評價模型的需求，并識別需要評估的數學求解方法、部件、現象、物理過程和參數。

2.1定義分析目的、瞬態類別和核電廠類型

COSINE最佳估算大破口LOCA評價模型應用的核電廠類型為AP1000/CAP1400系列的非能動壓水堆，兼顧考慮傳統二代3/4環路壓水堆。

2.2定義性能指標

由于針對的是大破口LOCA事故，其可接受的量化的性能指標主要為:

(1)包殼峰值溫度:計算的包殼峰值溫度不能超過2 200℉(約1 204.4℃)［1］。

(2)最大包殼氧化份額:計算的包殼局部最大氧化量不能超過包殼初始厚度的17%［1］。

(3)最大氫氣產生量:計算的包殼與水/蒸汽反應所產生的氫氣總量不能超過全部包殼金屬參與反應所產生的氫氣總量的1%［1］。

2.3現象識別與排序

通過對AP1000/CAP1400系列非能動反應堆大破口LOCA事故的研究，并兼顧考慮傳統二代3/4環路壓水堆的特性，結合國內外多名相關領域專家的專家判斷，識別并排序了COSINE最佳估算大破口LOCA評價模型的重要過程和現象［5］。

3　最佳估算大LOCA評估矩開陣發

基于所開發的最佳估算大破口LOCA評價模型PIRT表［5］，通過調研和分析世界范圍內可用的實驗數據，并參考國際上對評價模型確認與評估的研究結果［6］，按照國際上對大破口LOCA事故發展階段的劃分，對PIRT表中的每一個重要現象制定了相應的實驗確認和評估計劃，開發了最佳估算大破口LOCA評價模型評估矩陣，分別見表1～表3，其中標注“×”表示計劃采用此實驗對對應的現象進行評估。

表1　噴放階段的評估矩陣

表2　再灌水階段的評估矩陣

表3　再淹沒階段的評估矩陣

4　結語

我國首個自主開發的COSINE軟件包針對反應堆LOCA事故開發了最佳估算評價模型，將最先開展最佳估算大破口LOCA評價模型的評估。根據國際最新的EMDAP方法開展了最佳估算大破口LOCA評價模型能力需求分析，完成了最佳估算大破口LOCA評價模型PIRT表的開發。調研并評價國內外過去幾十年發展的實驗研究，按照大破口LOCA事故的發展階段開發了3個對應的評估矩陣，為COSINE最佳估算大破口LOCA評價模型的評估和不確定性分析奠定了基礎。

根據所開發的3個評估矩陣，最佳估算大破口LOCA評價模型的評估需要用到約40個不同臺架的實驗數據，針對如此龐大的數據需求，需要在國內外開展實驗數據收集;同時，還需要針對那些無法獲得而又必需的實驗數據開展補充實驗研究，以滿足COSINE最佳估算大破口LOCA評價模型的評估需求。

［1］USNRC.10 CFR 50.46 Acceptance criteria for emergency core cooling systems for light-water nuclear power reactors［S］.Washington DC:USNRC，2007.

［2］USNRC.10 CFR 50.46 Appendix K ECCS evaluation models［S］.Washington DC:USNRC，2000.

［3］USNRC.RG1.157 Best-estimate calculations of emergency core coolingsystemperformance［S］.WashingtonDC: USNRC，1989.

［4］USNRC.RG1.203 Transient and accident analysis methods［S］.Washington DC:USNRC，2005.

［5］Fu X L，Liu L F，Yu N，et al.Validation requirement analysis for COSINE system analysis code［C］.Pisa，Italy:The 15th International Topical Meeting on Nuclear Reactor Thermal-Hydraulics，2013.

［6］Aksan N，Auria F D，Glaeser H，et al.Separate effects test matrix for thermal-hydraulic code validation，vol.1:phenomena characterisation and selection of facilities and tests［R］.Paris: OECD NEA，1993.

Assessment Matrixes Development of COSINE Best-estimate LB-LOCA Evaluation Model

Fu Xiaoliang1，Liu Lifang1，Yu Nan1，Du Zheng1，He Siqi1，Dong Bo1，Liang Guoxing2，Yang Yanhua1
(1.State Nuclear Power Software Development Center，Beijing 102209，China; 2.School of Nuclear Science and Engineering，Shanghai Jiao Tong University，Shanghai 200240，China)

Major phenomena and processes in different stages of large break loss of coolant accident(LBLOCA)in pressurized water reactor are identified and ranked according to the latest EMDAP method proposed by NRC，based on which three assessment matrixes for best-estimate evaluation models corresponding to the propagation stages of LB-LOCA are developed in combination with analysis on existing experimental data.

best estimate;LB-LOCA;assessment matrix

TM623

A

1671-086X(2016)01-0031-04

2015-07-22

核電關鍵設計軟件自主化技術研究重大專項(2011ZX06004-024)

傅孝良(1985—)，男，工程師，研究方向為核能工程。

E-mail:fuxiaoliang@snptc.com.cn