林 樺,王云偉,張往鎖
(國核自儀系統(tǒng)工程有限公司,上海 200241)
壓水堆核電站功率控制系統(tǒng)對象研究及仿真分析
林 樺,王云偉,張往鎖
(國核自儀系統(tǒng)工程有限公司,上海 200241)
核電站反應堆功率控制系統(tǒng)是核電站的關鍵控制系統(tǒng)之一,在開環(huán)測試的環(huán)境下,控制系統(tǒng)無正確、有效的對象特性反饋,無法驗證功率控制系統(tǒng)的設計和功能實現(xiàn)是否正確?;赗EMARK堆芯物理程序,建立了某壓水堆核電站堆芯物理模型程序,為反應堆功率控制系統(tǒng)的設計和驗證提供控制對象。堆芯模型采用三維有限差分網(wǎng)格進行求解計算,兩群時間相關的三維擴散方程能夠對快中子和熱中子注量率行為進行準確模擬。將堆芯物理模型計算的反應性系數(shù),包括硼微分價值參數(shù)、慢化劑溫度系數(shù)、等溫溫度系數(shù)、多普勒功率系數(shù)以及功率調節(jié)棒組和停堆棒組的反應性價值,與核電廠的堆芯設計參考數(shù)據(jù)進行對比。通過對比分析發(fā)現(xiàn),模型的計算結果與電廠堆芯的設計參考數(shù)據(jù)吻合,說明所建立的堆芯模型能夠有效反映電廠反應堆的堆芯物理特性。該堆芯模型可與反應堆功率控制系統(tǒng)構成有效的閉環(huán)測試環(huán)境,為反應堆功率控制系統(tǒng)的設計和驗證提供有效手段。
能源; 核電站; 壓水堆; 反應堆堆芯; 功率控制; 閉環(huán)
核電站反應堆功率控制系統(tǒng)是核電站的關鍵控制系統(tǒng)之一,其通過改變控制棒的位置來實現(xiàn)對反應堆功率調節(jié)和功率分布控制,并且在執(zhí)行功率控制功能時不會觸發(fā)反應堆緊急停堆或引起蒸汽旁排閥門的開啟。
反應堆功率控制的關鍵是對堆芯的反應性控制,但在開環(huán)測試的環(huán)境下,由于控制系統(tǒng)無正確、有效的對象特性反饋,所以無法驗證功率控制系統(tǒng)的設計和功能實現(xiàn)是否正確。若將控制系統(tǒng)作用于虛擬的堆芯對象模型并進行閉環(huán)驗證測試,則可以有效實現(xiàn)對功率控制系統(tǒng)的設計驗證。由于堆芯的物理特性十分復雜,堆芯的慢化劑系數(shù)、多普勒功率系數(shù)、硼的價值以及控制棒的價值對堆芯反應性影響很大,一般的簡化模型無法準確、有效地反映堆芯的物理特性。
本文基于REMARK堆芯物理程序,對某大型壓水堆核電站的堆芯進行物理建模,并對堆芯模型的慢化劑溫度系數(shù)、等溫溫度系數(shù)、多普勒功率系數(shù)和虧損、總功率系數(shù)和虧損、可溶硼微分價值、裂變產物價值、控制棒積分和微分價值進行了仿真分析。仿真計算結果與堆芯設計分析數(shù)據(jù)的對比結果相吻合,證明該堆芯物理模型能夠有效反映參考電廠的堆芯物理特性。該模型可以應用于反應堆功率控制系統(tǒng)的設計和驗證,以及控制系統(tǒng)在反應堆不同運行壽期下的控制參數(shù)調整和優(yōu)化研究[1-4]。
功率控制系統(tǒng)的作用對象為反應堆堆芯。本文所述的控制對象,其堆芯熱功率為4 040 MW,堆芯由193組燃料組件組成,堆芯活性區(qū)高度為4 267.2 mm。每個燃料組件由264根燃料棒、24根控制棒導向管和1根儀表管,按17×17正方形柵格排列構成。每個燃料組件有8個中間格架、4個攪混格架、1個頂部格架、1個底部格架和1個保護格架。
首循環(huán)堆芯采用5區(qū)(A、B、C、D、E)布置,燃料組件數(shù)分別為12個、61個、32個、60個和28個。首循環(huán)采用了720根通水環(huán)狀可燃吸收體(wet annular burnable absorber, WABA)和8 944根整體型燃料可燃吸收體(integral fuel burnable absorber, IFBA)作為可燃毒物。堆芯采用的棒束控制組件由41束調節(jié)棒和48束停堆棒組成,共計89束(其中73束為黑控制棒,16束為灰控制棒)。棒組中的補償棒(M棒)包含6個子組,分別為MA(4束)、MB(4束)、MC(4束)、MD(4束)、M1(4束)、M2(8束)。軸向偏移控制棒(axial offset rod,AO棒)只有1組(13束)。其中,M棒組主要用于補償因燃耗、溫度、功率水平等運行條件的變化而導致的反應性變化,AO棒組主要用于堆芯軸向功率分布的控制。此外,中心1束AO棒還可用作快速降功率系統(tǒng)的備選子組。停堆棒(shutdown rod, SD棒)共有48束,包含6個子組,分別為SD1(8束)、SD2(8束)、SD3(8束)、SD4(8束)、SD5(8束)、SD6(8束),用于確保反應堆在任何功率水平下都有足夠的停堆裕量。在觸發(fā)快速降功率系統(tǒng)時,停堆棒可按照預先設定的子組落入堆芯,快速將堆芯功率降至較低水平,以避免堆芯停堆。調節(jié)棒組中的MA、MB、MC、MD為灰棒控制組件(gray rod cluster assembly,GRCA),其余調節(jié)棒和停堆棒均為黑棒控制組件(rod cluster control assembly,RCCA)[5]。
堆芯建模設計輸入數(shù)據(jù)為核設計計算程序ANC8.1的計算輸出數(shù)據(jù)。設計輸入數(shù)據(jù)的功率分布和燃耗分布可作為堆芯模型基準對比參數(shù)??刂瓢艚孛嬗蓧燮诔?beginning of life,BOL)、壽期中(middle of life,MOL)和壽期末(end of life,EOL)3個燃耗步提供, 熱態(tài)滿功率(hot full power,HFP)工況下的控制棒截面同樣適用于熱態(tài)零功率(hot zero power,HZP)工況。HFP分布文件對應平衡氙(equilibrium xenon,EQXE)工況,HZP水平對應無氙(no xenon,NOXE)工況。燃料組件平均功率及每個燃料組件的軸向功率分布數(shù)據(jù)作為模型初始工況的設計輸入。建模設計輸入所用的功率分布文件均對應控制棒全提(all rods out,ARO)工況。
堆芯物理模型程序REMARK可用于模擬反應堆正常、異常或應急運行時的響應,模型采用三維有限差分網(wǎng)格進行求解計算。采用與時間相關的兩群中子三維擴散方程(忽略對熱泄漏項的估算),能夠對快中子和熱中子注量率行為進行準確模擬。結合堆芯區(qū)域的三維網(wǎng)格結構,能準確地模擬出詳細的局部中子注量率響應。主要模型特性如下[6-12]。
①以兩群中子截面作為輸入數(shù)據(jù),換料分析數(shù)據(jù)作為數(shù)據(jù)源。壽期初到壽期末的每個堆芯壽期都有一組截面數(shù)據(jù)作為模型輸入,在燃料循環(huán)的不同周期內,燃料特性的改變都能夠準確地得到反映。
②外推注量率邊界。堆芯邊界附近的注量率能被準確地計算,可確保堆外探測器的準確讀數(shù)。反射層的性能包括在截面數(shù)據(jù)中。每個節(jié)點計算6組緩發(fā)中子,通過計算6組緩發(fā)中子的詳細分布,以中子通量密度作為中子擴散方程的分布源項的一部分,可以準確地模擬在所有條件下,特別是在停堆時的中子注量分布。
③反應性影響取決于慢化劑的溫度/密度/空泡、反應堆的壓力、燃料組件的溫度(多普勒)、氙毒、釤毒、硼濃度和控制棒。同時,要考慮每個節(jié)點燃料組件的位置、氙毒和釤毒的濃度、控制棒的位置、硼濃度、慢化劑的溫度/密度/空泡效應、燃料組件的溫度。
④3種同位素(U-235、U-238、Pu-239)的衰變熱計算、詳細的衰變熱空間分布能提供余熱排出工況時的準確溫度分布。
⑤堆外探測器的中子注量率讀數(shù)受到探測器位置和堆外水密度的影響。堆內探測器的中子注量率讀數(shù)取決于探測器的位置和局部中子注量率。堆芯功率的不均勻分布取決于不均勻的熱工水力條件和控制棒的方式。在每一個節(jié)點,計算放射性源項。
對首循環(huán)的核設計參數(shù)建立的堆芯物理模型進行壽期初運行工況下的反應性系數(shù)計算,并將計算結果與堆芯設計計算結果進行對比分析。
3.1 硼微分價值
當在HFP及HZP工況下計算硼微分價值(differential boron worth,DBW)時,將硼濃度變化控制在±25×10-6以內。HFP對應氙平衡EQXE工況,HZP對應無氙NOXE工況。在HFP和HZP工況下,微分硼價值曲線如圖1所示。通過與設計參考值對比分析發(fā)現(xiàn),硼價值與參考值在不同硼濃度下均能很好地吻合。
圖1 微分硼價值曲線
3.2 慢化劑溫度系數(shù)
堆芯的慢化劑平均溫度與功率水平關系如表1所示。計算慢化劑溫度系數(shù)(moderator temperature coefficient,MTC)時,將堆芯入口溫度變化控制在±5 K以內,同時將燃料和包殼溫度維持在參考工況。
表1 慢化劑平均溫度與功率水平關系表
MTC均在HZP、NOXE工況下計算,除M棒組外的其余棒組處于全提位置。HFP對應EQXE工況,HZP對應NOXE工況,均取各自工況下的臨界硼濃度。HFP和HZP工況下的慢化劑溫度系數(shù)(MTC)曲線如圖2所示。對比結果表明,兩者吻合效果較好。
圖2 慢化劑溫度系數(shù)曲線
3.3 等溫溫度系數(shù)
在不同燃耗下,等溫溫度系數(shù)(isothermal temperature coefficient,ITC)的大小隨功率水平及硼濃度變化。計算ITC時,將堆芯入口溫度變化控制在±5 K以內,同時將燃料和包殼溫度設為與慢化劑溫度同步變化。等溫溫度系數(shù)曲線如圖3所示。
圖3 等溫溫度系數(shù)曲線
3.4 多普勒功率系數(shù)
多普勒功率系數(shù)(doppler power coefficient,DPC)與功率水平及燃耗關系如表2所示。
表2 DPC與功率水平及燃耗關系表
計算DPC時,將反應堆額定功率(rated thermal power,RTP)變化控制在±5%以內,同時將燃料和包殼溫度設為與功率水平同步變化。計算時維持HFP工況下的EQXE和臨界硼濃度。本文所計算功率系數(shù)為首循環(huán)的壽期初HFP工況和HZP工況對應的多普勒功率系數(shù),分別為-8.442和-16.07,計算結果與表2中滿功率和零功率水平下的設計值吻合。
3.5 M棒組價值
在HFP和HZP工況下,計算M棒組價值時均維持HFP、ARO所對應的臨界硼濃度和平衡氙。M棒組的重疊步為:灰棒棒組(MA~MD)間重疊83步,黑棒棒組(M1和M2)間以及最后一組灰棒(MA或MD)與M1間重疊12步。M棒組價值曲線如圖4所示。通過與設計參考值的對比,對象模型的計算結果與之吻合。
圖4 M棒組價值曲線
3.6 AO棒組價值
在HFP和HZP工況下,計算AO棒組價值時均維持HFP、ARO所對應的臨界硼濃度和平衡氙。AO棒組價值曲線如圖5所示。通過與設計參考值的對比,對象模型的仿真計算結果與之吻合。
圖5 AO棒組價值曲線
3.7 SD棒組價值
在HFP和HZP工況下,計算SD棒組價值時均維持HFP、ARO所對應的臨界硼濃度和平衡氙。SD棒組價值曲線如圖6所示。通過與設計參考值的對比,對象模型的仿真計算結果與之吻合。
圖6 SD棒組價值曲線
3.8 裂變產物價值
主要裂變產物的價值計算考慮停堆后的氙價值計算。計算時,假定初始工況為HFP EQXE、臨界硼濃度,此時將功率水平直接降為HZP,在核素衰變過程中計算氙價值。停堆后氙價值曲線如圖7所示。通過與設計參考值的對比,對象模型的仿真計算結果與之吻合。
圖7 停堆后氙價值曲線
堆芯反應性的變化是影響反應堆功率控制系統(tǒng)設計的關鍵因素。由于壓水堆核電站反應堆堆芯的物理特性十分復雜,堆芯反應性受堆芯慢化劑系數(shù)、多普勒功率系數(shù)、硼價值以及控制棒價值的影響很大,一般的簡化模型無法準確反映實際物理特性。本文基于REMARK堆芯物理程序,建立了某壓水堆核電站的堆芯物理模型。其慢化劑溫度系數(shù)、等溫溫度系數(shù)、多普勒功率系數(shù)和虧損、總功率系數(shù)和虧損、可溶硼微分價值、裂變產物價值、控制棒積分和微分價值的計算結果與堆芯設計分析數(shù)據(jù)吻合,證明該堆芯物理模型正確,能夠有效反映參考電廠的堆芯物理特性。該模型可用于反應堆功率控制系統(tǒng)的設計和驗證工作,還可為控制系統(tǒng)在反應堆不同運行壽期下的控制參數(shù)調整和優(yōu)化研究工作提供研究對象。
[1] 鄭勇.壓水堆核電站堆芯物理/熱工水力耦合特性研究[J].原子能科學技術,2014,48(12):2298-2303.
[2] 周海翔,徐瑋瑛.三代核電機組數(shù)字化儀控系統(tǒng)及其國產化分析[J].自動化儀表,2010,31(8):61-66.
[3] 楊宗偉,黃鐵明,馮光宇,等.核電站仿真技術在反應堆控制系統(tǒng)調試中的應用[J].核動力工程,2009,30(6):49-53.
[4] 張建民.核反應堆控制[M].西安:西安交通大學出版社,2002.
[5] 林誠格.非能動安全先進壓水堆核電技術[M].北京:原子能出版社,2012.
[6] 伍宇忠,陳世和,李罡,等.基于多模型策略的壓水堆堆芯功率控制與仿真[J].核動力工程,2013,34(5):104-107.
[7] 劉道光.基于PAnySimu的PWR核電站堆芯建模與仿真研究[J].核電子學與探測技術,2011,31(9):1022-1028.
[8] 于濤.基于微機仿真壓水堆核電站堆芯物理數(shù)學模型的建立[J].清華大學學報(理工版),2001,15(1):28-34.
[9] 周海翔.反應堆溫度效應的穩(wěn)定性影響分析[J].控制工程,2013,20(3):423-425.
[10]王云偉,周海翔,王黎澤,等.非能動安全型核電站給水調節(jié)控制方案仿真驗證[J].自動化儀表,2013,34(3):62-65.
[11]何麗華,于濤,鄭衛(wèi)平,等.階躍反應性擾動及溫度變化下壓水堆堆芯響應特性的SIMULINK仿真[J].南華大學學報,2015,3(1):16-21.
[12]李罡,趙福宇,劉洋.壓水堆堆芯非線性系統(tǒng)的全局穩(wěn)定性分析[J].原子能科學技術,2013,47(2):271-276.
《自動化儀表》 郵發(fā)代號: 4-304, 2017年定價: 18.00元,全年價: 216.00元; 國外代號: M 721
歡迎賜稿, 歡迎訂閱, 歡迎寶貴建議, 歡迎惠刊各類廣告
Research and Simulation Analysis of Power Control System in Pressurized Water Reactor Nuclear Power Plant
LIN Hua,WANG Yunwei,ZHANG Wangsuo
(State Nuclear Power Automation System Engineering Company,Shanghai 200241,China)
The power control system of reactor core is one of the critical control systems in nuclear power plant,under the environment of open loop test,the control system does not have correct and effective object characteristics feedback,so it is impossible to verify whether the design and functional implementation of the power control system are correct or not.Based on REMARK reactor core physical program,the physical model program of the reactor core of certain PWR nuclear power plant is established,to provide control object for design and verification of the reactor power control system.The three-dimensional finite difference grid is used to solve the calculation,and the two-group of time dependent three-dimensional diffusion equation can accurately simulate the behavior of fast neutron and thermal neutron flux rate.The reactivity coefficients calculated by the model,including the moderator temperature coefficient isothermal temperature coefficient,Doppler power coefficient and Boron differential value,and all reactivity values of power control rod group and shutdown rod group are compared with the reference data in design of reactor core.Through the comparative analysis,it is found that the simulation results are identical with the design reference data,which shows that the model is correct and effective,and can reflect the physical characteristics of the reactor core.Effective closed loop test environment can be composed of this model and the power control system of the reactor core;the model provides effective means for designing and verifying the power control system of the reactor.
Energy; Nuclear power plant; Pressurized water reactor; Reactor core; Power control; Closed loop
國家重大科技專項基金資助項目(2013ZX06005001)
林樺(1983—),男,碩士,工程師,主要從事核電站數(shù)字化儀控系統(tǒng)驗證及核電站全范圍模擬機系統(tǒng)仿真相關工作。 E-mail:linhua@snpas.com.cn。
TH-3;TP391
A
10.16086/j.cnki.issn1000-0380.201705002
修改稿收到日期:2017-01-10