徐士坤,于 濤,謝金森,*,李滿倉,王連杰
(1.南華大學 核科學技術學院,湖南 衡陽 421001;2.南華大學 湖南省數字化反應堆工程技術研究中心,湖南 衡陽 421001;3.中國核動力研究設計院 核反應堆系統設計技術重點實驗室,四川 成都 610213)
核反應堆的功率水平與堆芯壽期是重要的性能指標。理想的核動力反應堆希望既具有較大功率又具有較長的堆芯壽期。板狀燃料組件傳熱面積大,燃料中心與冷卻劑通道距離短,傳熱效率高[1],在同等功率水平和冷卻劑流量下獲得更低的燃料芯體與表面溫度,或在燃料芯體與表面溫度不超出限值的情況下,允許核反應堆以更高的功率運行,同時板狀燃料組件結構緊湊,更利于堆芯的小型化[2-3]。板狀燃料組件是高通量研究堆、緊湊型高功率密度核動力反應堆常用的燃料幾何構型[4-7],采用板狀燃料組件的長壽期反應堆在對功率和壽期均有要求的應用場景具有吸引力。
堆芯壽期的增長可通過增加易裂變核素裝載量或提高堆芯可裂變核素轉換比來實現。在壓水堆中,由于轉換比遠小于1,單純提高轉換比的方法難以實現超長的堆芯壽期。通過增加易裂變核素裝載量增長堆芯壽期,也會遇到剩余反應性控制的難題,簡單地增加控制棒數量或硼酸濃度無法實現壽期內剩余反應性的平緩釋放,同時還面臨一回路壓力邊界泄漏概率增加、堆芯功率分布更加不均勻以及正的慢化劑溫度系數等問題。可燃毒物具有較好的中子吸收能力,吸收中子后的產物可以是對中子“透明”的核素,也可以是能繼續吸收中子的核素,因而具有靈活的剩余反應性控制能力[8]。可燃毒物設計是實現長壽期壓水堆堆芯方案剩余反應性控制的一種重要手段。國內外對壓水堆可燃毒物選型進行了相關研究[8-10],但目前大部分研究對象為傳統棒狀低富集度燃料組件,對于富集度較高的板狀燃料組件的可燃毒物選型,較優可燃毒物中子學特性研究較少。長壽期板狀壓水堆燃料組件具有更高的富集度,如果采用傳統壓水堆可燃毒物,無法滿足反應性平緩釋放等情況,因此需對可燃毒物進行重新選型研究,使得可燃毒物在壽期初不僅能抑制住更大的初始反應性,壽期內反應性的平緩釋放,并且在壽期內燃燒充分,降低壽期末反應性懲罰。
本文以長壽期壓水堆板狀燃料組件為研究對象,對采用不同可燃毒物方案的板狀燃料組件開展計算,通過組件剩余反應性-燃耗時間曲線分析不同可燃毒物的性能,獲得較優的可燃毒物,為長壽期板狀燃料壓水堆可燃毒物選型和設計提供參考。
DRAGON程序是由加拿大蒙特利爾大學開發的確定論反應堆組件-堆芯計算程序[11],DRAGON程序包含多種數值計算方法(碰撞概率法、球諧函數法、特征線法等),可用來求解中子輸運的一維(1D)、二維(2D)、三維(3D)問題。DRAGON程序對幾何的適應性和靈活性較高,可用來處理不同堆型下的結構幾何[12]。DRAGON程序對板狀元件計算的可行性已經得到了驗證[13-16],并且可對具有復雜燃耗的可燃毒物(如Gd等)進行燃耗計算且保證一定的計算精度[17-18]。本文程序版本選用DRAGON4.1.0進行計算,其中中子輸運計算方法采用EXCEL模塊,通量求解采用TYPE K模塊進行處理。
板狀燃料組件示意圖如圖1所示,組件由13塊燃料板、兩塊支撐板和水隙構成,其中燃料選擇為UO2彌散型燃料,包殼、支撐板和基體材料均為Zr-4合金。燃料芯體厚度為3 mm,燃料包殼厚度為0.4 mm,兩塊燃料板之間水隙厚度為2.3 mm,支撐板厚度為3 mm。燃料選擇為60%高富集度彌散型燃料,可燃毒物裝載形式選擇為與燃料均勻混合的彌散型可燃毒物。

圖1 組件示意圖Fig.1 Scheme of assembly
不同于傳統低富集度棒狀燃料組件壓水堆,因長壽期壓水堆堆芯富集度較高,組件初始反應性較大。從中子學角度,無論是天然核素還是富集同位素作為可燃毒物,選型應遵循以下3個原則:1) 可燃毒物應在壽期初可抑制住較大的初始反應性;2) 可燃毒物在壽期中需保證反應性的平緩釋放;3) 含可燃毒物組件在壽期末反應性懲罰較小或無懲罰可使反應堆達到較長的滿功率運行天數[19]。本文從壽期初抑制初始反應性、壽期內釋放的反應性和壽期末的反應性懲罰三方面來評價候選可燃毒物的相關物理性能,篩選出適用于高富集度長壽期壓水堆板狀燃料組件的可燃毒物。由于可燃毒物具有不同的中子吸收截面,因此通過調整燃料板中可燃毒物和燃料的含量,使得組件初始無限介質增殖因數kinf=1.20左右。
本文選擇非增殖可燃毒物、錒系核素可燃毒物、先進聚合物可燃毒物作為研究對象開展計算分析。
1) 非增殖可燃毒物
非增殖可燃毒物中各元素采用天然豐度,主要有B4C、Dy2O3、Er2O3、Eu2O3、Gd2O3、Sm2O3等,這些可燃毒物在組件中僅起到吸收中子的作用,在壽期末不會因子代核素的產生而延長壽期,即不會出現子代核素增殖成易裂變核素的情況。
2) 錒系核素可燃毒物
錒系核素是周期系ⅢB族中原子序數為89~103的15種化學元素的統稱,錒系核素中的一部分在吸收中子后可轉換為易裂變核素,以這類核素做可燃毒物時,可通過吸收中子在壽期初抑制過剩反應性,還能利用轉換出的易裂變核素在組件壽期中、末補充反應性,延長組件壽期,并實現剩余反應性的平緩變化。本文選擇錒系核素231Pa、241Am、237Np、238Pu、240Pu進行比較。
3) 先進聚合物可燃毒物
聚碳硼烷-硅氧烷-乙炔基(PACS)作為先進聚合物可燃毒物展現出良好的毒物特性[19]。PACS化學結構式如圖2所示。本文選擇目前較優的兩種先進聚合物可燃毒物組合:PACS-J和PACS-L,對先進聚合物可燃毒物的命名取自文獻[20],PACS中各元素的原子數/分子組成列于表1。

圖2 PACS化學結構式Fig.2 PACS chemical structure

表1 PACS組成Table 1 PACS composition
候選可燃毒物燃耗曲線如圖3所示。由圖3a可見,高燃料富集度情況下,錒系可燃毒物僅有231Pa在kinf=1.0時表現出增殖效應,其他錒系可燃毒物均出現了不同大小的反應性懲罰。由圖3b可見,kinf=1.0時,兩種先進聚合物可燃毒物在壽期末因自身特性延長了組件壽期。PACS-L的燃耗曲線因在整個組件壽期內反應性波動較大,反應性較難控制,不滿足反應堆安全運行要求。PACS-J的燃耗曲線在初始kinf=1.04時在整個壽期內的反應性波動未超過0.2,但反應性在壽期中后期出現較大釋放,雖然延長了組件壽期,但反應性控制存在困難。由圖3c可見,天然配比的非增殖可燃毒物在kinf=1.0時,除含可燃毒物B4C的燃耗曲線所能達到的燃耗深度與無可燃毒物燃耗曲線重合,其余非增殖可燃毒物在壽期末均造成了不同程度的反應性懲罰。這是因為采用天然配比的可燃毒物存在小吸收截面核素,且產生的子代核素具有大小不等的吸收截面和半衰期,減少了組件壽期。對于高富集度板狀燃料組件,采用天然配比的非增殖可燃毒物不是最優選擇。
因此,本文從中子學角度分析選擇天然配比的非增殖可燃毒物Dy2O3、Er2O3、Eu2O3、Gd2O3和Sm2O3對應的富集同位素氧化物164Dy2O3、167Er2O3、151Eu2O3、157Gd2O3和149Sm2O3,并進行組件燃耗計算,計算結果如3d所示。由圖3d可見,在高燃料富集度下,采用富集同位素氧化物作為可燃毒物時,富集同位素消耗得更快,在壽期末造成的反應性懲罰更小,組件的EFPD較采用天然配比可燃毒物所能達到的更大。采用富集同位素作為可燃毒物后,減少了燃耗鏈中子代核素的產生,使得一些具有吸收截面的子代核素減少,降低了壽期末的反應性懲罰,提高了燃料利用率。
從以上研究可篩選出性能較優的可燃毒物,然后分別進行中子學研究。231Pa、B、157Gd、167Er和PACS-J表現出了較優的中子學性能。在整個燃耗周期內,B4C和167Er作為可燃毒物時反應性釋放平緩,在整個壽期內無明顯波動,在壽期末幾乎未造成反應性懲罰;157Gd作為可燃毒物時,在壽期末幾乎沒有造成反應性懲罰,但反應性在壽期中釋放較快;231Pa和PACS-J作為可燃毒物時,由于231Pa自身燃耗鏈的關系和PACS-J自身特性使得組件壽期得到了延長,超過無可燃毒物燃料組件,提高了燃料利用率。但當PACS-J作為單一可燃毒物時,壽期內的反應性控制存在困難。

a——錒系可燃毒物;b——PACS可燃毒物;c——非增殖可燃毒物;d——富集同位素可燃毒物圖3 可燃毒物燃耗曲線Fig.3 Burn-up curve of burnable poison
為定量分析較優可燃毒物相關物理性能,不同可燃毒物分析方案列于表2,不同數量的可燃毒物組件布置示意圖如圖4所示。

表2 不同可燃毒物分析方案Table 2 Analysis solution of different burnable poisons

可燃毒物:a——157Gd2O3;b——B4C;c——167Er2O3;d——PACS-J圖4 不同可燃毒物組件布置示意圖Fig.4 Layout scheme of different burnable poison assemblies
157Gd的中子吸收截面是所選可燃毒物中最大的。當157Gd作為可燃毒物時,在壽期初期僅加入5.20%的157Gd2O3即可抑制住初期剩余反應性,但大的吸收截面使得可燃毒物消耗較快,在壽期的中后期出現較大的反應性釋放,產生反應性波動,同時157Gd吸收中子后產生的子代核素存在一定的中子吸收截面,在壽期末會造成少量反應性懲罰。
10B的中子吸收截面在兩種同位素中最大,B的燃耗鏈非常簡單,吸收中子后直接轉化為本質上對中子透明的氦和鋰,壽期末基本不存在反應性懲罰,但作為彌散型可燃毒物需考慮氦氣釋放導致的燃料元件腫脹問題。
167Er的中子吸收截面相比較157Gd而言較小,但具有一定的中子吸收截面,相較于157Gd而言燃燒較慢,減少單塊板中可燃毒物含量,使得167Er2O3在壽期末盡可能的消耗完全。同時167Er吸收中子后產生的子代核素存在一定的中子吸收截面,在壽期末會造成少量反應性懲罰。
231Pa的中子吸收截面相較于157Gd和167Er而言相對較小,因此燃燒速度明顯較慢,可對反應性實現較好的控制。231Pa的燃耗鏈與其他核素燃耗鏈不同[22],231Pa作為此燃耗鏈的母核,燃耗鏈中直接存在可裂變的232U和233U,231Pa俘獲中子后生成232Pa,232Pa的半衰期僅為1.3 d,衰變形成232U,232U的俘獲和裂變截面如圖5所示,232U的裂變截面和俘獲截面非常相近,這樣使得232U既能俘獲中子生成易裂變核素233U,也可直接裂變補充反應性。231Pa這種特殊的燃耗鏈為含231Pa組件實現高燃耗提供了可能性。

圖5 232U微觀俘獲和裂變截面Fig.5 Microscopic capture and fission cross sections of 232U
PACS-J具有較高的含氫量與高的含硼量,因自身具有較高的含氫量,使得PACS-J作為可燃毒物時自身具有慢化效應,在壽期末可更好地慢化中子,提高熱中子利用率,并促進裂變核素Pu的消耗。因此PACS-J作為可燃毒物時在壽期末可延長壽期,提高燃料利用率。
當PACS-J作為單一可燃毒物進行反應控制時,因為燃料富集度較高使得壽期初始反應性較大,雖然增加可燃毒物在單塊板中的含量來增加可燃毒物自屏效應,但組件反應性波動依舊是所選可燃毒物中最大的,因此在高燃料富集度情況下可考慮PACS-J與慢燃耗可燃毒物搭配組合。本文對可燃毒物組合進行了初步研究,不改變單塊板中可燃毒物含量,僅調整組件中含PACS-J、167Er2O3和231Pa2O3可燃毒物板數量。其中PACS-J與167Er2O3組合的組件中,含PACS-J和167Er2O3的可燃毒物板數量分別為1塊和4塊,其余為純燃料板元件;PACS-J與231Pa2O3組合的組件中,含PACS-J和231Pa2O3的可燃毒物板數量分別為1塊和10塊,其余為純燃料板元件,組件布置示意圖如圖6所示。

a——PACS-J與231Pa2O3組合;b——PACS-J與167Er2O3組合圖6 組件含不同可燃毒物搭配示意圖Fig.6 Layout scheme of assembly with different numbers of burnable poison
圖7為PACS-J與167Er2O3和231Pa2O3可燃毒物組合后的反應性燃耗曲線,可看出PACS-J與中子吸收截面較小的可燃毒物組合后,組件初始kinf為1.10左右,在整個壽期內組件反應性無較大波動,反應性在壽期內平緩釋放。在組件壽期末沒有造成反應性懲罰,且延長了組件壽期,其中PACS-Er組合在壽期末延長了6 EFPD,PACS-Pa組合在壽期末延長了21 EFPD。兩種可燃毒物組合方案的反應性波動在整個壽期內均為0.1左右。

圖7 可燃毒物組合燃耗曲線Fig.7 Burn-up curve of burnable poison combination
PACS-J作為可燃毒物具有較大的中子吸收截面,同時因組件具有較大的初始反應性,如果采用中子吸收截面較大的單一可燃毒物進行反應性控制,在壽期初可以抑制住初始剩余反應性,但隨著燃耗加深,可燃毒物釋放反應性速率較快,組件出現較大的反應性波動,當快燃耗可燃毒物與慢燃耗可燃毒物組合可以彌補高富集度下組件反應性控制存在困難的缺點,從而獲得更優的反應性燃耗曲線。
本文基于高富集度長壽期壓水堆板狀燃料組件,開展可燃毒物選型和中子學特性研究,對不同可燃毒物進行輸運-燃耗計算,篩選出中子學性能較優的可燃毒物,具體結論如下。
1) 對于板狀燃料組件,在高燃料富集度下,可采用富集同位素中157Gd、167Er的氧化物作為可燃毒物來減小組件末期的反應性懲罰;B4C作為可燃毒物時,在壽期末幾乎無反應性懲罰;同時231Pa和PACS-J作為可燃毒物時,在壽期末均延長了組件壽期,提高了燃料利用率。在較為經濟地獲取157Gd、167Er、231Pa和PACS-J時可作為長壽期壓水堆的候選可燃毒物。
2) 在高燃料富集度下,部分快燃耗可燃毒物不適合在組件中作為單一可燃毒物進行反應性控制,但可通過與慢燃耗可燃毒物的合理組合,獲得更優的反應性燃耗曲線。如先進聚合物可燃毒物PACS可與慢燃耗可燃毒物Er和Pa進行組合來增加對反應性的控制能力。