先進非能動壓水堆核電站燃料組件軸向燃耗分布研究

楊海峰,霍小東,于 淼

(中國核電工程有限公司,北京100840)

先進非能動壓水堆核電站燃料組件軸向燃耗分布研究

楊海峰,霍小東,于 淼

(中國核電工程有限公司,北京100840)

燃料組件的軸向燃耗分布以及末端效應是燃耗信任制技術應用中的難點?；谙冗M非能動壓水堆核電站的運行模式及組件設計特點,結合可能的燃料管理策略,統(tǒng)計軸向兩端使用低富集度抑制區(qū)的乏燃料組件,生成軸向燃耗包絡線。以AP1000堆型的乏燃料貯存單元為例,通過分析統(tǒng)計,證明生成的軸向燃耗包絡線用于臨界安全分析是保守的。在此基礎上,詳細研究燃料組件頂部的低富集度抑制區(qū)對末端效應的貢獻,并對燃料組件進行設計改進,減小至消除末端效應,為簡化乏燃料組件相關的臨界安全分析提供了一個方法。相關研究工作及成果,是先進非能動壓水堆核電站乏燃料組件相關的設施設備的臨界安全設計的基礎,可為其他堆型的相應研究提供參考和借鑒。

軸向燃耗分布;末端效應;燃耗信任制;密集貯存

在臨界安全分析中,考慮燃料組件在堆內輻照導致的反應性降低,即為燃耗信任制。相對于簡單的新燃料組件假設,燃耗信任制技術可在保證安全的前提下,顯著提高乏燃料組件貯存、運輸、后處理的經(jīng)濟性,但卻大大增加了臨界安全分析的難度和復雜度。其中,燃料組件在堆內輻照造成的軸向燃耗非均勻分布及相應的末端效應是一個關鍵的難點。

通常情況下,堆芯頂部慢化劑溫度較高、密度較低,燃耗很淺。先進非能動壓水堆核電站,如AP1000[1]等,采用機械補償模式(MSHIM),依靠控制棒來調節(jié)反應性變化,減少硼化和稀釋的頻率從而簡化硼相關系統(tǒng)。在核電站循環(huán)壽期的相當長時間內,控制棒都插入堆芯,使得燃料組件頂部的燃耗更淺。同時,為減少中子泄漏并提高經(jīng)濟性,燃料組件的軸向兩端使用低富集度抑制區(qū),使得乏燃料組件的軸向燃耗分布更為復雜。因此,不能直接套用現(xiàn)有的基于軸向均一富集度的乏燃料組件的軸向燃耗包絡線。

本文基于先進非能動壓水堆核電站的燃料管理策略,統(tǒng)計分析含軸向抑制區(qū)的乏燃料組件,生成軸向燃耗包絡線,并詳細研究軸向抑制區(qū)對末端效應的貢獻和臨界安全分析的影響。所開展的工作,是先進非能動壓水堆核電站乏燃料組件相關的密集貯存格架、運輸容器等設施設備的臨界安全設計的基礎。

1 計算程序

使用從美國引進的APA程序系統(tǒng)完成燃料管理的計算。在APA程序系統(tǒng)中,主要有組件計算程序PARAGON[2、3]、堆芯計算程序ANC[4、5]、接口程序 ALPHA[6]等。PARAGON程序是一個模塊化的、兩維多群輸運理論的組件計算程序,用于生成組件截面參數(shù)。ALPHA程序根據(jù)用戶輸入,建立組件計算模型,調用PARAGON程序執(zhí)行計算,為ANC程序準備截面數(shù)據(jù)庫。ANC程序用于堆芯的建模和燃料管理計算。

使用從英國引進的蒙特卡羅程序MONK-9A[7]進行臨界安全分析。MONK-9A程序通過模擬中子的產(chǎn)生、輸運、死亡來計算系統(tǒng)的有效增殖因子(keff),同時也可計算中子流、反應率、界面流等。

2 堆芯燃料管理方案

先進非能動壓水堆核電站,大多采用低泄漏裝載模式,采用18個月或16/20個月交替的平衡循環(huán)燃料管理策略,以提高核電站運行的經(jīng)濟性。

在兩種燃料管理策略下,從第二循環(huán)開始,燃料組件中均不再配置離散型可燃毒物,依靠整體型可燃毒物(IFBA,燃料芯塊表面的ZrB2涂層)來控制功率分布,補償循環(huán)堆芯壽期初的后備反應性。普通燃料棒(芯塊表面不含ZrB2涂層)和IFBA棒(芯塊表面含Zr B2涂層)的軸向兩端20.32 cm的區(qū)段是富集度為3.20%的軸向抑制區(qū),區(qū)別在于普通燃料棒的軸向抑制區(qū)使用實心芯塊,IFBA棒的軸向抑制區(qū)使用環(huán)狀芯塊。燃料棒的主體部分,富集度較高,且軸向兩端各有10.16 cm的區(qū)段不含ZrB2涂層。在本文中,燃料組件富集度特指燃料棒主體部分的富集度。

在進行堆芯燃料管理計算時,將堆芯軸向分為42個4英寸的節(jié)塊。為了模擬基負荷MSHIM運行模式下控制棒的插入情況,同時簡化計算分析,采用最佳估算MSHIM模型進行堆芯燃耗計算。目前,我國的核電站大部分以滿功率基負荷模式運行,很少進行負荷調整。因此,最佳估算MSHIM模型貼近我國核電站的實際情況,既考慮控制棒的插入效應,又不至于過分保守。

2.1 18個月平衡循環(huán)

18個月平衡循環(huán)方案[1],每次裝入36組富集度4.45%和28組富集度4.95%的新燃料組件,其中4.45%富集度的燃料組件分別有128、156個IFBA棒,4.95%富集度的燃料組件分別有64、128個IFBA棒。平衡循環(huán)的循環(huán)長度為20 200兆瓦天每噸初始鈾(MWd/t U),對應的循環(huán)壽期為500有效滿功率天(EFPD)。平衡循環(huán)堆芯裝載如圖1所示。

圖1 18個月平衡循環(huán)堆芯裝載圖Fig.1 18-Month Equilibrium Cycle Fuel Loading

圖2 16個月循環(huán)堆芯裝載圖Fig.2 16-Month Equilibrium Fuel Loading

2.2 16/20個月交替平衡循環(huán)

16/20個月交替平衡循環(huán)方案[1]中,16個月循環(huán)堆芯每次裝入33組富集度4.27%和24組富集度4.70%的新燃料組件,其中4.27%富集度的燃料組件分別有80、104、128個IFBA棒,4.70%富集度的燃料組件分別有48、104個IFBA棒。20個月循環(huán)堆芯每次裝入36組富集度4.67%和36組富集度4.95%的新燃料組件,其中4.67%富集度的燃料組件分別有156、200個IFBA棒,4.95%富集度的燃料組件分別有104、128個IFBA棒。16個月循環(huán)堆芯的循環(huán)長度為17 750 MWd/t U,對應的循環(huán)壽期為439EFPD。20個月循環(huán)堆芯的循環(huán)長度為22 750 MWd/t U,對應的循環(huán)壽期為563EFPD。16個月循環(huán)的堆芯裝載圖如圖2所示,20個月循環(huán)的堆芯裝載圖如圖3所示。

圖3 20個月循環(huán)堆芯裝載圖Fig.3 20-Month Equilibrium Fuel Loading

3 軸向燃耗包絡線

大部分壓水堆乏燃料組件具有相似的軸向燃耗分布,中間燃耗很深,兩端燃耗較淺。在臨界安全分析中,軸向燃耗分布模型與平均燃耗模型的keff之差即為末端效應。末端效應是燃耗信任制技術應用中的一個關鍵因素。末端效應的考慮,可使用軸向燃耗包絡線來保守地考慮,也可在平均燃耗模型的計算結果上疊加一個Δkeff來保守地考慮。

3.1 軸向燃耗包絡線生成

目前已知的軸向燃耗包絡線,都是基于軸向均一富集度的乏燃料組件的數(shù)據(jù)匯編而來的,不適合于包含軸向低富集度抑制區(qū)的乏燃料組件相關的臨界安全分析。

基于先進非能動壓水堆可能的燃料管理策略,包括18個月平衡循環(huán)和16/20個月交替平衡循環(huán)等,進行堆芯燃料管理計算,統(tǒng)計首循環(huán)、過渡循環(huán)、平衡循環(huán)等堆芯帶軸向抑制區(qū)的乏燃料組件的卸料燃耗及軸向燃耗分布信息。乏燃料組件的統(tǒng)計如圖4所示,組件富集度為3.70%～4.95%,卸料燃耗為30 000～57 000 MWd/tU。

基于統(tǒng)計的乏燃料組件信息,生成軸向燃耗包絡線。首先,對每個乏燃料組件的軸向燃耗分布進行歸一化處理。經(jīng)匯總統(tǒng)計,活性區(qū)底部的5個節(jié)塊和頂部的7個節(jié)塊,平均歸一化燃耗值小于1.0或者在1.0左右,因此,針對此12個節(jié)塊,統(tǒng)計所有乏燃料組件的相應節(jié)塊的最小歸一化燃耗值,作為相應節(jié)塊包絡的歸一化燃耗值。針對中間30個節(jié)塊,統(tǒng)計所有乏燃料組件的相應節(jié)塊的最大歸一化燃耗值,作為相應節(jié)塊的歸一化燃耗值。然后,按公式(1)計算中間30個節(jié)塊的重新歸一化因子F,最后將中間30個節(jié)塊的歸一化燃耗值乘以重新歸一化因子F,作為相應節(jié)塊包絡的歸一化燃耗值。最終得到歸一化的軸向燃耗包絡線,如圖5所示。

式中:Bui——節(jié)塊i的歸一化燃耗。

所統(tǒng)計的乏燃料組件是源自同一堆型的同類型組件,其軸向燃耗分布形狀相似,因此得到的軸向燃耗包絡線是一個典型的、適合于該堆型及其燃料組件的軸向燃耗包絡線。

3.2 包絡線的保守性分析

為了驗證生成的軸向燃耗包絡線,選定兩個典型組件放入AP1000堆型的乏燃料貯存格架中,分析匯總統(tǒng)計的每一個軸向燃耗分布和軸向燃耗包絡線對應的keff。典型組件I,初始富集度5.0%、含128根IFBA棒。典型組件II,初始富集度4.5%、含156根IFBA棒。根據(jù)燃料管理的計算結果,結合貯存格架的裝載曲線,分別選擇42 600 MWd/t U和35 500 MWd/t U作為兩個典型組件的平均燃耗。

組件燃耗計算采用硬化的中子能譜,使得指定燃耗深度下的核素密度用于臨界安全分析是保守的。熱功率為3 400 MWt(對應于157個組件),慢化劑溫度為323.37℃,對應的慢化劑密度為0.670 9 g/cm3,平均硼濃度為1 000 ppm,控制棒插入。

AP1000堆型的乏燃料貯存格架的貯存單元有兩種:單獨的套管貯存單元和組合貯存單元。套管貯存單元為不銹鋼方管,方管外壁附著含B4C的鋁基中子吸收板(硼鋁板),硼鋁板外有不銹鋼包覆。硼鋁板厚度為0.55 cm,寬度為21.50 cm。組合貯存單元是單獨的套管貯存單元對角焊接的結果。貯存單元中心距為25.00 cm。貯存單元的計算模型如圖6所示。在貯存單元計算模型的四周設置周期性邊界條件,相當于平面上有無窮多個貯存單元,頂部和底部設置足夠厚的水反射層。

圖6 AP1000堆型乏燃料貯存單元計算模型示意圖Fig.6 Calculation model of the spent fuel storage cell of AP1000

在進行臨界分析時,采用錒系加裂變產(chǎn)物的置信水平,考慮主要的錒系核素與裂變產(chǎn)物,忽略強吸收、短壽命的裂變產(chǎn)物以及一些非主要的裂變產(chǎn)物。由于忽略強吸收、短壽命的裂變產(chǎn)物的反應性貢獻,不考慮乏燃料組件的冷卻是保守的。因為在冷卻過程中,155Gd的累積以及241Pu衰變成241Am等引入的負反應性大于由239Np衰變成239Pu引入的正反應性。

臨界計算時,設置足夠的抽樣粒子數(shù)與計算代數(shù),保證計算結果的統(tǒng)計偏差不大于0.000 3。

基于典型組件Ⅰ和典型組件Ⅱ的統(tǒng)計分析結果如圖7、圖8所示。在圖7、圖8中,除了給出每個軸向燃耗分布對應的keff外,還給出所有keff的平均值,加上3倍統(tǒng)計偏差以后的值,以及軸向燃耗包絡線對應的keff值。相對于平均keff,軸向燃耗包絡線引入了很大的保守性,遠遠超過3倍統(tǒng)計偏差。相對于每一個的軸向燃耗分布,軸向燃耗包絡線也是足夠保守的,只有極少的燃耗分布對應的keff與軸向燃耗包絡線對應的keff在3倍統(tǒng)計偏差之內。

基于所匯編的軸向燃耗包絡線,研究末端效應隨平均燃耗的變化。計算結果如表1所示,對于上述兩個典型燃料組件,平均燃耗分別高于35 000 MWd/t U、30 000 MWd/t U時,平均燃耗模型不再保守。相對于軸向均一富集度的乏燃料組件,包含軸向低富集度抑制區(qū)的乏燃料組件,其平均燃耗模型保守與不保守的轉換點有所提高。主要是對keff貢獻很大的頂部區(qū)域采用了低富集度抑制區(qū),使得平均燃耗模型在更大的燃耗范圍內是保守的。

圖7 典型組件Ⅰ的軸向燃耗分布對應的k eff統(tǒng)計Fig.7 Summary of k eff of the axial burnup profiles for typical spent fuel assemblyⅠ

圖8 典型組件Ⅱ的軸向燃耗分布對應的keff統(tǒng)計Fig.8 Summary of keff of the axial burnup profiles for typical spent fuel assemblyⅡ

表1 典型組件不同燃耗深度下的keff值Table 1 keff of different burnup of typical spent fuel assemblies

3.3 軸向抑制區(qū)的影響

與基于軸向均一富集度的乏燃料組件的軸向燃耗分布相比,燃料組件軸向兩端低富集度抑制區(qū)的存在,一方面減小了組件軸向兩端,特別是頂部的燃耗,使其從平均燃耗的50%降至30%,增大了末端效應。另一方面,軸向抑制區(qū)初始富集度的降低又減小了末端效應。

根據(jù)圖4所示的乏燃料組件統(tǒng)計,使用42 600 MWd/t U的平均燃耗,采用同樣的軸向低富集度區(qū)域、不同的主體部分富集度,研究末端效應的變化趨勢。結果如表2所示。從計算結果可以看出,在相同的軸向抑制區(qū)下,隨著主體部分富集度的不斷提高,即主體部分與軸向抑制區(qū)富集度差距越來越大,末端效應越來越小。

表2 軸向抑制區(qū)對末端效應的影響Table 2 Impact of axial blanket on end effect

統(tǒng)計上述計算的裂變份額軸向分布,如圖9所示。從圖中可以看出,裂變峰集中于活性區(qū)頂部區(qū)域。在燃料組件主體部分富集度與頂部抑制區(qū)富集度差別比較小的情況下,趨近于軸向均一富集度的情況,頂部幾個節(jié)塊的裂變份額很高。隨著主體部分富集度與頂部抑制區(qū)富集度的差別越來越大,頂部節(jié)塊的裂變份額逐漸減小。另外,初始富集度越低,達到同樣的燃耗深度,燃料組件頂部裂變份額越高,末端效應越大。

燃料組件頂部2個節(jié)塊為低富集度抑制區(qū),且有較高的中子泄漏。第3個節(jié)塊為高富集度主體部分,且距離組件活性區(qū)頂端有一定的距離,泄漏比較低。因此,第3個節(jié)塊的裂變份額很高。為了降低頂部節(jié)塊的裂變份額,減小末端效應,考慮適當減小燃料組件頂部第3個節(jié)塊的初始富集度,將原4.50%～4.95%富集度統(tǒng)一降為4.0%富集度,將原4.00%～4.45%富集度統(tǒng)一降為3.5%富集度。在不改變堆芯燃料裝載方案的基礎上,進行堆芯燃料管理計算。結果表明,平衡循環(huán)長度減小約100 MWd/t U,對應于2.5EFPD,對于焓升因子FΔH、熱點因子FQ等參數(shù)的影響很小。由于只改變整個燃料組件42個節(jié)塊中的一個節(jié)塊,且處于頂部,對于整個燃料管理的影響很小,對于軸向燃耗分布的影響也很小,只有頂部2個節(jié)塊的歸一化燃耗減小約0.01,第3個節(jié)塊的歸一化燃耗減小約0.04。研究第3個節(jié)塊的富集度降低對末端效應的影響,計算結果如表3所示。

圖9 裂變份額的軸向分布Fig.9 The axial distribution of the fission fraction

表3 改進組件不同燃耗深度下的keff計算結果Table 3 keff of different burnup of improved spent fuel assemblies

圖10 改進組件的裂變份額軸向分布Fig.10 The axial distribution of the fission fraction for the modified fuel assemblies

統(tǒng)計其軸向裂變份額的分布,如圖10所示。與圖9對比可知,頂部第3個節(jié)塊降低富集度后,頂部裂變份額明顯降低,使得改進組件在典型卸料燃耗下,末端效應為零或為負值,從而可以在臨界安全分析中使用組件平均燃耗模型。根據(jù)此研究結果,在進行組件設計時,可考慮適當降低燃料組件頂部芯塊的富集度,一方面提高燃料組件的經(jīng)濟性,另一方面可有效簡化乏燃料組件相關的臨界安全分析。不利的方面在于增加了燃料組件的復雜度。

4 結論

基于先進非能動壓水堆的燃料管理策略,包括18個月平衡循環(huán)和16/20個月交替的平衡循環(huán),統(tǒng)計含軸向低富集度抑制區(qū)的乏燃料組件,生成先進非能動壓水堆乏燃料組件的軸向燃耗包絡線。基于AP1000乏燃料貯存單元,經(jīng)統(tǒng)計分析,證明生成的軸向燃耗包絡線具有足夠的保守性。詳細研究了頂部低富集度抑制區(qū)對末端效應的貢獻和臨界計算的影響,并在此基礎上,對燃料組件的設計進行改進,減小甚至消除末端效應,使得簡單的平均燃耗模型可等效復雜的軸向燃耗分布模型,為先進非能動壓水堆核電站乏燃料組件相關的臨界計算的簡化提供一個新的方法。所開展的工作及所取得的研究結果,可應用于先進非能動壓水堆核電站乏燃料組件相關的臨界設計,研究方法及研究結果,對于其他類型的核電站也有一定的參考作用。

[1] 孫漢虹等．第三代核電技術AP1000[M]．北京:中國電力出版社,2012．

[2] M．Ouisloumen,et al． PARAGON:The New Westinghouse Assembly Lattice Code[C]．Salt Lake City,International Meeting on Mathematical Methods for Nuclear Application,2001

[3] M.Ouisloumen,H．Huria,L.T．Mayhue,et al．Qualification of the Two-Dimensional Transport Code PARAGON[R]．[S.l.]:Westinghouse Electric Company,2004

[4] Y.S．Liu,et al．ANC:A Westinghouse Advanced Nodal Computer Code[R]．[S.l.]:Westinghouse Electric Company．2012

[5] T.Q．Nguyen,K.C．Hoskins,M.M.Weber,et al．Qualification of the phoenix-p/anc nuclear design system for pressurized water reactor cores[R]．[S.l.]:Westinghouse Electric Company,1988

[6] Westinghouse Electric Company．ALPHA User Manual[M]．2012

[7] ANSWERS.MONK,A Monte Carlo Program for Nuclear Criticality Safety and Reactor Physics Analyses[M]．[S.l.]:ANSWERS,2006．

Study of the axial burnup profiles for passive advanced PWR

YANG Hai-feng,HUO Xiao-dong,YU Miao
(China Nuclear Power Engineering Co．Ltd．,Beijing 100840,China)

The axial burnup profiles of the fuel assemblies and the associated end effect play an important role in burnup credit(BUC)．Based on the operational strategy and the fuel assemblies'characteristics of the passive advanced PWR,as well as the possible fuel management schemes,the burnup information of the spent fuel assemblies including low enrichment axial blankets is summarized,and the bounding axial burnup profile is generated．The spent fuel storage cell model of AP1000 is used to statistically evaluate the adequacy of the bounding profile,and the results demonstrate that the bounding axial burnup profile is conservative for criticality analyses．And the impact of the axial blankets on the end effect is studied in detail．Based on the results,fuel assemblies are improved to mitigate and eliminate the end effect,which provides a method to simplify the criticality analyses．The results can be used for the criticalityanalyses related to the spent fuel assemblies of the passive advanced PWR,and the study can be referenced by other PWRs．

2015-07-23

楊海峰(1981—),男,河南南陽人,高級工程師,碩士,現(xiàn)從事反應堆物理及臨界安全分析工作

Axial Burnup Profile;End Effect;Burnup Credit;High Density Storage

TL371

A

0258-0918(2016)01-0279-08