核電廠運行瞬態(tài)分析軟件PANTO的研發(fā)

劉余 任春明 湯琪芬 鄧堅

?

核電廠運行瞬態(tài)分析軟件PANTO的研發(fā)

劉余1,2任春明2湯琪芬2鄧堅1,2

1（中國核動力研究設(shè)計院核反應(yīng)堆系統(tǒng)設(shè)計技術(shù)重點實驗室 成都610213） 2（中國核動力研究設(shè)計院 成都610213）

針對核電廠運行瞬態(tài)分析的功能需求，中國核動力研究設(shè)計院研發(fā)了PANTO (Program for Analysis of Normal Transient and Overpressure)軟件。該軟件基于成熟可靠的系統(tǒng)分析模型和特殊部件模型，采用模塊化的軟件設(shè)計理念，應(yīng)用面向?qū)ο蟮腃++語言和java語言，具有完全自主知識產(chǎn)權(quán)。PANTO軟件通過了單元測試、集成測試和系統(tǒng)測試，基本消除了所有的代碼缺陷。針對秦山二期核電廠階躍負荷增大10%與額定功率下全部甩負荷瞬態(tài)試驗進行了驗證計算。結(jié)果表明，PANTO軟件能夠較好地模擬瞬態(tài)中關(guān)鍵參數(shù)的變化情況，計算精度滿足工程應(yīng)用要求，適用于壓水堆核電廠運行瞬態(tài)分析。

運行瞬態(tài)，軟件研發(fā)，PANTO

運行瞬態(tài)主要是指核電廠在調(diào)試或功率運行期間經(jīng)?；蚨ㄆ诎l(fā)生的瞬態(tài)事件，比如核電廠升溫升壓、階躍負荷變化、甩負荷等，它屬于典型的I類工況，在瞬態(tài)過程中不允許觸發(fā)反應(yīng)堆保護或?qū)ＴO(shè)安全設(shè)施。設(shè)計必須事先開展運行瞬態(tài)分析，以驗證核電廠控制系統(tǒng)能夠進行有效的調(diào)節(jié)，并滿足運行要求。而對于運行瞬態(tài)分析，通常依賴于專門的分析軟件，比如用于法國M310核電廠的CATIA2[1]軟件、用于美國AP1000核電廠的CENTS[2]軟件等。

為了全面實現(xiàn)核電設(shè)計自主化，中核集團于2011年啟動核電工程軟件自主化研發(fā)項目[3]，其中包括運行瞬態(tài)分析軟件PANTO (Program for Analysis of Normal Transient and Overpressure)的研發(fā)。該軟件針對壓水堆核電廠，模擬的范圍包括反應(yīng)堆堆芯、一回路系統(tǒng)、部分二回路系統(tǒng)、控制及保護系統(tǒng)等，它可用于I類和II類工況瞬態(tài)分析、控制系統(tǒng)優(yōu)化、穩(wěn)壓器及蒸汽發(fā)生器安全閥的設(shè)計評價等。

本文主要介紹PANTO軟件的研發(fā)概況，包括物理模型、軟件設(shè)計、軟件編碼、軟件測試與軟件驗證等。

1 物理模型

作為系統(tǒng)級的設(shè)計與分析軟件，PANTO涉及的物理模型較多。研發(fā)中一方面借鑒了同類軟件中成熟可靠的模型，即采用集總參數(shù)方法，假設(shè)流動和傳熱為一維過程；另一方面，對于存在兩相流的穩(wěn)壓器和蒸汽發(fā)生器，分別建立了部件模型。

1.1 一回路冷卻劑流量計算模型

假設(shè)在給定的時間步長內(nèi)冷卻劑密度為常數(shù)，并且系統(tǒng)壓力保持不變，對一回路建立一維動量方程：

式中：、、、分別為控制體內(nèi)的壓力、質(zhì)量、流通面積和自由容積；為質(zhì)量流量；為流體速度；c為重力加速度。

對于封閉的一回路有：

式中：H為泵壓頭；H為提升壓降；H為摩擦壓降；H為局部壓降。

將一回路劃分為段控制體，總的摩擦壓降為：

式中：K為第段的阻力系數(shù)；W為質(zhì)量流量；ρ為流體密度。

由此，便得到一回路流量計算的控制方程：

除了上述計算模型外，PANTO軟件還包括系統(tǒng)流量計算的經(jīng)驗?zāi)Ｐ?，或用戶直接輸入隨時間變化的流量曲線。

1.2 中子動力學模型

對于堆芯核功率的計算，采用點堆中子動力學方程，考慮6組緩發(fā)中子[4]：

(6)

(7)

式中：()為中子通量；為反應(yīng)性；為總的緩發(fā)中子份額；β為第組緩發(fā)中子份額；λ為第組緩發(fā)中子先驅(qū)核衰變時間常數(shù)；C為第組緩發(fā)中子先驅(qū)核濃度；為中子源項；為中子代時間。

1.3 堆芯傳熱計算模型

假設(shè)燃料在徑向只有一個節(jié)點，堆芯換熱方程如下[5]：

式中：f為燃料質(zhì)量；pf為燃料的比熱容；f為燃料的溫度；為總的核功率；為燃料產(chǎn)生核功率的系數(shù)，=0.974；A為等效熱導(dǎo)；c,w為冷卻劑溫度。

采用有限差分方法，在時間步長時的燃料溫度f可以寫成c,w的函數(shù)，假設(shè)c,w取前一時間步長的值，在第個軸向節(jié)點上有：

(10)

1.4 穩(wěn)壓器模型

穩(wěn)壓器的模擬采用兩區(qū)模型[6]，并假設(shè)氣區(qū)和液區(qū)都是均勻混合的。

1) 壓力方程

式中：pr為穩(wěn)壓器壓力；為一回路系統(tǒng)壓力；為波動管摩擦阻力系數(shù)；SU為波動管流通面積；SU為波動管體積流量，流向穩(wěn)壓器為正方向；SU為波動管質(zhì)量流量。

2) 質(zhì)量守恒方程

每個時間步長結(jié)束時穩(wěn)壓器內(nèi)水和蒸汽的質(zhì)量分別為：

(13)

式中：0、0為?Δ時刻水和蒸汽質(zhì)量；、為時刻的水和蒸汽質(zhì)量。

液相質(zhì)量變化為：

氣相質(zhì)量變化為：

(15)

3) 能量守恒方程

每個時間步長結(jié)束時水和蒸汽的比焓分別為：

(17)

1.5 蒸汽發(fā)生器模型

蒸汽發(fā)生器模型在一次側(cè)包括多個傳熱管區(qū)，二次側(cè)采用點模型，假設(shè)里面充滿飽和汽液混合物。一次側(cè)每個節(jié)點上的換熱量為：

式中：Q為每個節(jié)點的熱流量；P,j為一次側(cè)每個節(jié)點的溫度；S為二次側(cè)蒸汽飽和溫度。

假定二次側(cè)處于飽和狀態(tài)。若二次側(cè)流體溫度已知，則流體壓力、比焓與比體積等可通過水的物性計算得到。根據(jù)質(zhì)量守恒方程和能量守恒方程推導(dǎo)二次側(cè)飽和溫度變化率，由溫度變化率就可以求得每個時刻的二次側(cè)飽和溫度。

質(zhì)量守恒方程：

能量守恒方程：

(20)

式中：為蒸汽發(fā)生器流體質(zhì)量；SG為進入蒸汽發(fā)生器的熱流量；fw、s為給水流量和蒸汽流量；fw、s為給水比焓和蒸汽比焓；為蒸汽發(fā)生器流體平均內(nèi)能。

聯(lián)合質(zhì)量、能量方程可推導(dǎo)得到二次側(cè)飽和溫度變化率：

式中：為蒸汽發(fā)生器二次側(cè)總體積；f和g分別為飽和水與飽和蒸汽比焓；f與g分別為飽和水比容與飽和蒸汽比容；為蒸汽壓力。

對于式(21)中與物性相關(guān)的兩個變量：

(23)

可根據(jù)水物性參數(shù)，建立與飽和溫度s相關(guān)的擬合多項式。

1.6 控制系統(tǒng)模型

PANTO軟件開發(fā)考慮了核電廠一、二回路主要的控制系統(tǒng)，包括[7]：1) 反應(yīng)堆功率控制系統(tǒng)；2) 反應(yīng)堆溫度控制系統(tǒng)；3) 穩(wěn)壓器壓力與液位控制系統(tǒng)；4) 蒸汽發(fā)生器液位控制系統(tǒng)；5) 蒸汽排放控制系統(tǒng)；6) 給水流量控制系統(tǒng)。

這些控制系統(tǒng)在運行瞬態(tài)中根據(jù)設(shè)定邏輯調(diào)節(jié)相關(guān)參數(shù)使得核電廠趨于穩(wěn)定運行。目前，PANTO中內(nèi)置了針對國內(nèi)M310電廠的控制系統(tǒng)邏輯，瞬態(tài)模擬時只需要輸入相關(guān)的控制參數(shù)，比如時間常數(shù)、增益系數(shù)和目標值等。

2 軟件設(shè)計

PANTO軟件總體流程設(shè)計如圖1所示。采用模塊化設(shè)計理念，軟件設(shè)計上考慮了三個層次。第一層次分為4大模塊：計算輸入（讀入）、數(shù)據(jù)初始化、核心計算和結(jié)果輸出；第二層次主要針對核心計算模塊進行了細分，包括：一回路熱工水力計算、蒸汽發(fā)生器相關(guān)計算、堆芯中子動力學計算、給水設(shè)備熱工水力計算、蒸汽排放系統(tǒng)計算、控制保護系統(tǒng)及輔助系統(tǒng)計算；第三層次則是對第二層次各計算模塊的再分解，比如對于一回路熱工水力計算模塊分為：硼濃度計算、冷卻劑流量計算、穩(wěn)壓器計算、熱段與波動管流量計算、堆芯溫度和熱流密度計算等。

另一方面，為了方便用戶使用，設(shè)計了圖形化的輸入界面，如圖2所示。采用列表框和標簽頁相結(jié)合的模式，包括8個分頁面：1) 名義工況與模型選擇；2) 熱工水力參數(shù)；3) 蒸汽發(fā)生器參數(shù)；4) 堆芯物理和中子學參數(shù)；5) 保護和控制系統(tǒng)參數(shù)；6) 蒸汽管線模型；7) 給水設(shè)備特性參數(shù)；8) 初始工況和瞬態(tài)特性數(shù)據(jù)。

圖1 PANTO軟件總體流程圖

圖2 PANTO軟件圖形化輸入界面

3 軟件編碼

考慮到計算部分和圖形化輸入界面部分的設(shè)計目標不同，PANTO軟件計算部分采用了C++語言編寫，開發(fā)環(huán)境為Oracle Solaris Studio。對于輸入界面，為了同其它自主化軟件保持統(tǒng)一的風格，同時有利于相似模塊的直接繼承使用，該部分采用Java語言編寫，開發(fā)環(huán)境為NetBeans IDE 7.1.2。PANTO軟件代碼行總數(shù)約為91900行，其中計算部分代碼約為52500行，輸入界面部分代碼約為39400行。

軟件編譯后運行在中國核動力研究設(shè)計院的核反應(yīng)堆設(shè)計計算軟件平臺之上，操作系統(tǒng)為SUN Solaris 5.10，支撐軟件為Java Runtime Environment 6u10或更高版本，硬件環(huán)境為SUN SPARC- Enterprise。

4 軟件測試

根據(jù)軟件開發(fā)規(guī)范，必須開展代碼級的軟件測試，以驗證軟件實現(xiàn)的功能是否滿足需求說明書和設(shè)計說明書的要求，代碼編寫是否滿足編碼規(guī)范。為此，分三個階段進行了PANTO軟件的測試。第一階段為單元測試，測試各單元代碼對相應(yīng)編碼規(guī)范的遵循情況；第二階段為集成測試，驗證軟件計算結(jié)果的正確性和人機界面的正確性；第三階段為系統(tǒng)測試，對軟件進行總體測試。

PANTO軟件測試過程中共發(fā)現(xiàn)了12個缺陷。其中單元測試過程中發(fā)現(xiàn)兩個，系統(tǒng)測試過程中發(fā)現(xiàn)10個。這些缺陷經(jīng)過開發(fā)人員修復(fù)后再進行回歸測試，實現(xiàn)遺留缺陷數(shù)目為零的目標。最終測試的結(jié)果表明，PANTO軟件源代碼結(jié)構(gòu)清晰，可讀性強，編碼格式滿足規(guī)范要求；軟件計算結(jié)果準確，對異常計算能進行相應(yīng)錯誤提示并終止計算過程；軟件圖形界面友好，具備良好的人機交互能力；軟件功能完整，實現(xiàn)了軟件需求分析中所列功能點。

5 軟件驗證

以秦山二期核電廠為對象，選取了調(diào)試試驗中的階躍負荷增大和額定功率全部甩負荷瞬態(tài)進行了驗證計算。秦山二期核電廠為兩環(huán)路設(shè)計，額定工況下的主要參數(shù)如表1[8]所示。

表1 秦山二期核電廠主參數(shù)

兩個瞬態(tài)計算所需的邊界條件取自調(diào)試試驗結(jié)果，包括一回路流量、汽機第一級壓力、汽機進氣流量、給水流量和給水溫度，如圖3?7所示。對于堆芯中子學參數(shù)（如反應(yīng)性反饋系數(shù)、調(diào)節(jié)棒價值等），由于試驗中沒有實測結(jié)果，驗證計算中采用了設(shè)計值。

階躍負荷增大瞬態(tài)的計算結(jié)果對比如圖8?10所示。前200 s為穩(wěn)態(tài)計算，假設(shè)瞬態(tài)從200 s開始（對應(yīng)于邊界條件的零時刻），瞬態(tài)計算360 s。可以看出，由于二回路汽機負荷突然增大，反應(yīng)堆冷卻劑平均溫度和穩(wěn)壓器壓力略有降低，但同時由于功率控制系統(tǒng)發(fā)出調(diào)節(jié)信號，功率調(diào)節(jié)棒快速提升，導(dǎo)致堆芯核功率快速升高，反應(yīng)堆冷卻劑平均溫度逐漸上升，并達到滿功率運行設(shè)定值，反應(yīng)堆最終穩(wěn)定運行在100%滿功率水平。

圖3 一回路流量

圖4 汽機第一級壓力

圖8 堆芯核功率

圖9 反應(yīng)堆冷卻劑平均溫度

圖10 穩(wěn)壓器壓力

額定功率全部甩負荷瞬態(tài)的計算結(jié)果對比如圖11?13所示。同樣，前200 s為穩(wěn)態(tài)計算，瞬態(tài)計算360 s。瞬態(tài)開始，由于汽機負荷迅速下降，蒸汽發(fā)生器排熱量急劇減少，反應(yīng)堆冷卻劑平均溫度略有升高，其與設(shè)定值間的溫差信號將觸發(fā)蒸汽旁排系統(tǒng)開啟，以緩解一回路、二回路間的能量失配。瞬態(tài)參數(shù)最終將穩(wěn)定在廠用負荷水平。在320s后，相對于試驗結(jié)果，PANTO計算得到的堆芯和功率偏大，主要原因是由于在該時刻后功率控制的調(diào)節(jié)量較小，而計算又采用的是理論的中子學參數(shù)，特別是調(diào)節(jié)棒的價值跟實際值存在一定的差異（比實際值小）。堆芯核功率偏大也導(dǎo)致了反應(yīng)堆冷卻劑平均溫度和穩(wěn)壓器壓力偏大。

圖11 堆芯核功率

圖12 反應(yīng)堆冷卻劑平均溫度

圖13 穩(wěn)壓器壓力

從上述兩個瞬態(tài)計算結(jié)果跟試驗數(shù)據(jù)的對比中可以看出，關(guān)鍵參數(shù)的變化趨勢符合較好，該軟件準確地反映瞬態(tài)特征，在參數(shù)值上與試驗結(jié)果的偏差較小，滿足工程設(shè)計應(yīng)用的要求。由此驗證PANTO軟件開發(fā)的正確性。

6 結(jié)語

PANTO軟件的物理模型完備，它采用模塊化的軟件設(shè)計理念，應(yīng)用了面向?qū)ο蟮腃++語言和java語言，開發(fā)過程規(guī)范，并通過軟件測試。針對秦山二期核電廠階躍負荷增大和額定功率全部甩負荷瞬態(tài)，完成軟件驗證分析。計算結(jié)果表明，PANTO軟件能夠較好地模擬瞬態(tài)中反應(yīng)堆堆芯核功率、冷卻劑平均溫度和穩(wěn)壓器壓力的變化，計算精度滿足工程應(yīng)用要求，該軟件適用于壓水堆核電廠正常運行瞬態(tài)分析。

1 張英, 陳智, 周祖鑒, 等. 用于核電廠數(shù)字化儀表控制系統(tǒng)優(yōu)化設(shè)計的CATIA2程序的改進[J]. 核動力工程, 2008, 29(1): 19?23. ZHANG Ying, CHEN Zhi, ZHOU Zujian,. Improvement of CATIA2 code used in optimization design for digital instrument and control system of NPP[J]. Nuclear Power Engineering, 2008, 29(1): 19?23.

2 劉立欣, 鄭利民, 周全福. AP1000核電廠典型的運行瞬態(tài)分析[J]. 核技術(shù), 2012, 35(11): 869?876. LIU Lixin, ZHENG Limin, ZHOU Quanfu. Preliminary study on operational transient analysis for AP1000[J]. Nuclear Techniques, 2012, 35(11): 869?876.

3 中國核動力院. 中核集團召開核電軟件自主化項目啟動會[EB/OL]. 2011-03-11[2016-05-13]. http://www.cnnc. com.cn/publish/portal0/tab426/ info52484.htm. Nuclear Power Institute of China. CNNC held the kick-off meeting of self-reliance nuclear power plant design software project[EB/OL]. 2011-03-11[2016-05-13]. http:// www.cnnc.com.cn/publish/portal0/ tab426/info52484.htm.

4 謝仲生. 核反應(yīng)堆物理分析[M]. 北京: 原子能出版社, 1980: 310?311. XIE Zhongsheng. Nuclear reactor physics analysis[M]. Beijing: Atomic Energy Press, 1981: 310?311.

5 朱繼洲. 核反應(yīng)堆安全分析[M]. 西安: 西安交通大學出版社, 2004: 42?43. ZHU Jizhou. Nuclear reactor safety analysis[M]. Xi’an: Xi’an Jiaotong University Press, 2004: 42?43.

6 Abdallah M A, Ahmed H M, Mohammad A R. Pressurizer transients dynamic model[J]. Nuclear Engineering and Design, 1982, 73(3): 447?453.

7 張森如, 王建渝. 核電廠控制調(diào)節(jié)系統(tǒng)在運行瞬態(tài)分析中的數(shù)值計算[J]. 核科學與工程, 1988, 8(4): 324?329. ZHANG Senru, WANG Jianyu. Digital calculation of control and adjustment system in operational transient analysis for nuclear power plant [J]. Chinese Journal of Nuclear Science and Engineering, 1988, 8(4): 324?329.

8 李經(jīng)緯. 秦山核電二期工程反應(yīng)堆熱工水力設(shè)計[J]. 核動力工程, 1999, 20(4): 308?312. LI Jingwei. Reactor thermal hydraulic design for Qinshan Phase II nuclear power project[J]. Nuclear Power Engineering, 1999, 20(4): 308?312.

Development of operational transient analysis code PANTO for nuclear power plant

LIU Yu1,2REN Chunming2TANG Qifen2DENG Jian1,2

1(Science and Technology on Reactor System Design Technology Laboratory, Nuclear Power Institute of China, Chengdu 610213, China) 2(Nuclear Power Institute of China, Chengdu 610213, China)

Background: The operational transient of a nuclear power plant mainly refers to transients or events that occur frequently or periodically during commissioning or power operation. The transient process involves the response of primary loop, secondary loop and the control system. Usually, particular software is used for transient analysis. Purpose: This study aims to develop the program for analysis of normal transient and overpressure (PANTO) software that fulfils the requirements of operational transient analysis. Methods:Based on the mature and reliable system analysis model and special component model, object-oriented C++ language and java language were applied to design the modular software with completely independent intellectual property rights for PANTO. Through the unit testing, integration testing and system testing, all detected code defects were eliminated. Validation calculations were carried out for the 10% load step increase and complete load ejection at full power transient test which performed in the Qinshan Phase II nuclear power plant. Results:The results show that the PANTO can simulate the key parameters variation in the transients, and the calculation accuracy can meet the requirements of the engineering application. Conclusion: The PANTO is applicable foroperational transient analysis of PWR nuclear power plant.

Operational transient, Software development, PANTO

TL333，TL364，TP311.52

10.11889/j.0253-3219.2017.hjs.40.050604

劉余，男，1983年出生，2010年于清華大學獲博士學位，主要從事反應(yīng)堆熱工水力和安全分析研究

2017-02-13，

2017-03-02

Supported by Special Technical Project of China National Nuclear Corporation (No.LTJH-NESTORII-1405001)

中核集團公司重點科技專項（中核科發(fā)[2011]29號，No.LTJH-NESTORII-1405001）資助

LIU Yu, male, born in 1983, graduated from Qsinghua University with a doctoral degree in 2010, focusing on reactor thermal hydraulic andsafety analysis

2017-02-13, accepted date: 2017-03-02