白成斐,王 磊,付學峰,王 闖,張松文,李軍德,張 洪
(1.中科華核電技術研究院有限公司,廣東 深圳518026;2.環境保護部核與輻射安全中心,北京 100082;3.中廣核工程有限公司,廣東 深圳518124)
寧德核電站一號機組首次在首循環使用釓作為可燃毒物,燃料管理方案采用18個月換料模式,因此在啟動物理試驗項目設置上按照“系列首堆”進行考慮,同時借鑒大亞灣核電站和嶺澳核電站相關經驗[1],寧德核電站一號機組首次物理啟動試驗基本按原型堆的項目進行。在調試運行之前,根據相關調試大綱、試驗導則和細則以及燃料組件出廠性能參數,計算了反應堆物理試驗所涉及的各項測量參數的理論預計值(包括臨界硼濃度、控制棒價值、反應性系數和堆芯功率分布等)。反應堆物理試驗的理論分析一方面可以對試驗實施進行指導,另一方面通過理論預計值與實測值的比較,可以達到驗證堆芯核設計的正確性以及計算方法和程序的有效性的目的。
寧德核電站一號機組反應堆于2012年11月24日21時28分順利達到首次臨界,其后進行了一系列物理試驗。
本文將介紹該電站一號機組反應堆首循環啟動物理試驗的理論分析結果及其與實測值的比較。
寧德核電站一號機組反應堆熱功率為2 895MW,堆芯由157個AFA 3G燃料組件構成;為了展平堆芯功率分布,首循環堆芯燃料組件按235U富集度分三區裝載,富集度為1.8%,2.4%和3.1%的燃料組件數分別為53,52和52。較低富集度的兩種燃料組件按不完全棋盤格式排列在堆芯內區,最高富集度的燃料組件裝在堆芯外圍。首循環采用Gd2O3作為可燃毒物吸收體材料,并均勻彌散在較低富集度的UO2芯塊內。首循環預計的循環長度為332個等效滿功率天(EFPD),反應堆從第二循環開始快速向18個月換料過渡。反應堆共布置61束控制棒,控制棒布置如圖1所示。
寧德核電站一號機組首循環啟動物理試驗計算分析包括兩大部分內容:① 物理試驗所需參數的計算;② 反應堆物理試驗的預計分析或模擬計算。為了完成上述計算分析內容,根據反應堆首循環堆芯裝載、堆芯各區燃料組件235U富集度的實測值等參數,采用SCIENCE V2核程序包[2],進行了三維堆芯計算。

圖1 堆芯控制棒的布置Fig.1 Location of rod cluster control assemblies in core
SCIENCE V2核程序包主要包括:APOLLO2-F、SMART和ESPADON等程序。APOLLO2-F程序采用碰撞幾率方法進行組件輸運計算。對于一個燃料組件,該程序求解99群輸運方程,并為SMART程序提供兩群均勻化的中子截面參數。SMART是一個三維兩群堆芯擴散—燃耗計算程序,采用先進節塊法技術,可以進行穩態和瞬態工況的計算。ESPADON是一維二群穩態擴散燃耗程序。處理軸向幾何,節點數可達100個。它能處理水溫度、多普勒(Doppler)、氙和釤等非線性反饋效應和控制棒等效毒物截面。
SCIENCE V2核程序包已通過大量的基準例題和實驗測量結果的檢驗,其計算精度滿足工程設計要求。
反應堆首次臨界試驗,在熱態零功率(HZP)狀態下,首先將停堆棒組和功率補償棒組提出堆芯,將溫度調節棒組R置于130步位置,然后進行硼稀釋操作,最后通過提升R棒組使堆芯達到臨界[3]。表1中給出了首次臨界時對應的硼濃度及R棒組的棒位。

表1 首次臨界點Table 1 Initial criticality
在啟動物理試驗分析中,零功率下的臨界硼濃度測量包括所有控制棒全部提出堆芯狀態(ARO)、R棒全部插入堆芯狀態(R in)、R棒和G1棒全部插入堆芯狀態(R+G1in),所有試驗結果都在設計允許范圍之內,表2列出了HZP工況和不同棒位下堆芯臨界硼濃度的理論計算值與實測值的比較[4]。從表2可以看出,計算結果與實測值符合得非常好(驗收準則為±50ppm)。

表2 臨界硼濃度計算值與實測值(ppm)Table 2 Calculated and measured values of critical boron concentration(ppm)
表3給出了HZP狀態下控制棒價值計算值與測量值的比較。控制棒價值驗收準則要求范圍為計算值的±10%。從表中可以看出,寧德核電站一號機組控制棒價值的最大計算偏差為3.50%,試驗結果滿足驗收準則要求。

表3 控制棒價值計算值與測量值 (pcm)Table 3 Calculated and measured values of control rod worth(pcm)
堆芯慢化劑溫度系數直接影響反應堆的自穩性能,而慢化劑溫度系數的測量是通過等溫溫度系數的測量來獲得,所以測得的等溫溫度系數包含了慢化劑溫度系數和多普勒溫度系數的綜合效應。表4給出了HZP工況和不同插棒狀態下的等溫溫度系數計算值與測量值的比較。從表中可看出最大計算偏差為0.7pcm/℃,遠小于該試驗的驗收準則±3.6pcm/℃。

表4 等溫溫度系數(aiso)計算值與實測值 (pcm/℃)Table 4 Calculated and measured values of isothermal temperature coefficient(aiso)(pcm/℃)
表5給出了在滿功率狀態下等溫溫度系數與功率系數之比(aiso/ap)的計算值與實測值。此試驗的目的是為了驗證反應堆的自穩特性。在100%功率平臺完成,在30分鐘內,維持所有控制棒不動,堆芯硼濃度不變的條件下,降低40MW電功率,反應堆功率變化引入的正反應性完全用堆芯平均溫度升高引入的負反應性補償。該試驗的評價標準是等溫溫度系數與功率系數比值的測量值與理論值的相對偏差小于20%,由表可知,試驗結果滿足驗收準則的要求。

表5 aiso/aP計算值與測量值Table 5 aiso/aPcalculated and measured values
落棒試驗是驗證堆芯中子通量密度快速降低能夠正確引發反應堆停堆保護。該試驗在50%功率平臺進行。寧德核電站一號機組采用D12和M12位置的控制棒束突然落下的方法引發堆芯中子通量密度的快速降低,試驗中采用記錄儀和數字化儀控系統(DCS)共同記錄的方法,驗證測量通道能正確引發反應堆停堆保護。該試驗的驗收準則是四個功率量程通道中至少有三個通道觸發閾值保護。
表6給出了落棒試驗的試驗結果,通過表中數據可以看出,四個功率量程通道都觸發了閾值,滿足安全準則要求。

表6 落棒試驗結果Table 6 Results of rod drop test
選取功率補償棒組G1中D08位置的棒束作為HZP工況下模擬彈棒的棒束,模擬彈棒試驗后核焓升因子F△H和熱點因子FQ實測結果及堆芯功率分布測量結果見表7。其中,彈棒后熱點因子FQ需滿足彈棒事故的安全準則,即FQ(測量值)×1.114<15.5(安全分析值)。

表7 模擬彈棒試驗結果Table 7 Results of pseudo rod ejection test
從表中可以看出,模擬彈棒后,對于相對功率P≥0.9的組件,功率測量值與理論計算值最大偏差為5.6%;對于相對功率P<0.9的組件,測量值與理論計算值最大偏差為4.7%,均在驗收準則范圍之內(驗收準則分別為±10%和±15%);彈棒后,堆芯FQ=3.40,滿足彈棒事故的安全準則。
模擬落棒試驗的目的是檢驗單束控制棒組件由于機械或電器故障落入堆芯而引發堆芯焓升因子F△H或徑向功率峰值因子FXY的改變,以作為安全分析或設計改進的試驗數據。模擬落棒試驗是在50%滿功率下,全部功率控制棒提至頂部,選取安全停堆棒組(SB)中一束棒(G13位置)下插到堆芯底部,在下插過程中利用稀釋法維持反應堆臨界。當這束棒下插到堆底后用堆芯儀表系統(RIC)測量堆芯的功率分布,將其結果與參考堆芯功率分布測量結果(下插前)進行比較。
表8給出了理論預計值與實測值的比較。從表中可以看出,模擬落棒后,對于相對功率P≥0.9的組件,功率測量值與理論計算值最大偏差為3.1%;對于相對功率P<0.9的組件,測量值與理論計算值最大偏差為5.7%,均在驗收準則范圍之內(驗收準則分別為±5%和±8%),同時,堆芯F△H變化計算值與實測值也符合良好。

表8 模擬落棒試驗結果Table 8 Result of Pseudo rod drop test
在反應堆啟動調試階段,將進行零功率和各功率平臺下的功率分布測量試驗,這些試驗結果可以檢驗堆芯燃料組件裝載的正確性;并檢驗堆芯核設計的正確性,確保反應堆的安全運行。本文計算了各種堆芯狀態下的三維堆芯功率分布,并將理論計算結果與各個功率平臺下利用堆芯儀表系統測量的堆芯功率分布情況進行比較,從而驗證堆芯在相應的功率水平運行時各項安全指標和設計指標滿足驗收準則的要求。驗證結果表明:熱態零功率狀態下,最大計算偏差為6.4%(相對功率P≥0.9的邊緣組件),滿足驗收準則(±10%);而在反應堆功率運行狀態下,最大計算偏差為4.0%,出現在10%滿功率功率平臺(相對功率P<0.9的邊緣組件),滿足驗收準則(±15%)。當反應堆運行在滿功率狀態下,反應堆功率分布的最大計算偏差僅為2.7%,也滿足驗收準則(±8%)。由此可見,反應堆功率分布理論預計值與實測值符合得非常好。
作為世界上首例在首循環采用釓可燃毒物的大型商用壓水堆核電站,進行了一系列的反應堆物理試驗,測得了大量的表征反應堆堆芯設計性能的參數,如堆芯臨界硼濃度、控制棒價值、反應性系數和堆芯功率分布,同時還進行了模擬落棒試驗、模擬彈棒試驗以及落棒試驗等。
試驗結果表明:所有理論預計值與實測值符合良好,滿足有關驗收準則,由此證明寧德核電站一號機組首循環含釓堆芯設計滿足各項設計準則和安全準則,反應堆堆芯核設計是安全、可靠的。
[1] 廉志坤,芮旻 .大亞灣核電站18個月換料啟動物理試驗改進[J].核動力工程,2003,24(1):12-14.
[2] S.Rauck,B.Barbier.SCIENCE V2Nuclear Code Package -Qualification Report [S].FRAMATOME ANP,2004.
[3] 何濤,史永謙,朱慶福,等 .在臨界實驗中本底計數對倒數外推臨界的影響[J].核動力工程,2003,24(6).
[4] 章宗耀 .大亞灣核電站堆芯換料設計[J].核動力工程,2000,21(1).