郭 輝,馮快源,顧漢洋,*
(1.上海交通大學(xué) 核科學(xué)與工程學(xué)院,上海 200240;2.法國(guó)原子能與可替代能源委員會(huì),法國(guó) 卡達(dá)拉什 13115)
快中子反應(yīng)堆能有效提高鈾資源利用效率、減少高放廢物并有助于穩(wěn)定燃料循環(huán)[1-2],是第4代核能系統(tǒng)技術(shù)的重要組成部分,也是我國(guó)重點(diǎn)發(fā)展的先進(jìn)堆型之一[3]。截至目前,國(guó)內(nèi)外共建造了近20臺(tái)各型快堆,積累了400余堆年的運(yùn)行經(jīng)驗(yàn),其中主要為使用钚鈾氧化物混合(MOX)燃料的鈉冷快堆[4-5]。小型模塊化反應(yīng)堆堆芯的尺寸較小,可實(shí)現(xiàn)模塊化制造,具有建造周期短、投資風(fēng)險(xiǎn)低、運(yùn)輸靈活、部署方便等優(yōu)勢(shì)[6-7]。隨著核電技術(shù)的發(fā)展,固有安全的多用途小型模塊化快堆正成為國(guó)內(nèi)外研發(fā)的重點(diǎn)[8-10]。
快堆一般使用以碳化硼(B4C)為吸收劑的控制棒作為反應(yīng)性控制系統(tǒng)[11]。大型快堆通常具有較好的增殖特性,燃耗反應(yīng)性損失較小,可供安裝控制棒的空間較大,天然B4C一般能滿足大型快中子堆的反應(yīng)性控制需求。小型快堆中子泄漏率偏大,增殖能力偏弱,小型實(shí)驗(yàn)快堆對(duì)運(yùn)行周期長(zhǎng)度要求較低,但多用途小型模塊化快堆通常要求較長(zhǎng)的換料周期,因此小型模塊化快堆燃耗反應(yīng)性損失較大[12]。同時(shí)為實(shí)現(xiàn)緊湊堆芯的設(shè)計(jì),小型模塊化快堆可供安裝控制棒的空間有限。為實(shí)現(xiàn)小型模塊化快堆的反應(yīng)性控制,通常需要高10B富集度的B4C作為吸收劑。而高10B富集度的B4C功率密度大、輻照損傷嚴(yán)重,因此其安全裕度較小,不能完全滿足小型模塊化快堆長(zhǎng)期穩(wěn)定運(yùn)行的需求。
本文基于已有新型控制棒研究成果[13],針對(duì)320 MWth小型模塊化鈉冷快堆SMSFR-V0[14],提出以B4C為吸收劑的傳統(tǒng)控制棒和以硼化鉿(HfB2)或氫化鉿(HfH1.62)為吸收劑的新型控制棒概念設(shè)計(jì),并從反應(yīng)性價(jià)值、功率分布、反應(yīng)性反饋系數(shù)、溫度裕度及吸收劑燃耗深度等方面具體分析控制棒的中子特性與安全特性。
SMSFR-V0是以MOX為燃料的320 MWth小型模塊化鈉冷快堆概念設(shè)計(jì)方案[14-15],其堆芯布置如圖1所示,其燃料組件設(shè)計(jì)整合了PHENIX與ASTRID的經(jīng)驗(yàn)反饋。燃料活性區(qū)頂部含有25 cm鈉反射層,用以增加失鈉事故中的軸向中子泄漏率。
圖1 燃料組件軸向布局(a)與堆芯徑向布局(b)Fig.1 Axial layout of fuel assembly (a) and radial layout of SMSFR-V0 (b)
反應(yīng)性控制系統(tǒng)包括6組補(bǔ)償調(diào)節(jié)控制棒(CSD)和6組安全控制棒(DSD)。在SMSFR-V0中,這兩種控制棒使用相同的設(shè)計(jì),但它們?cè)诳刂七壿嬌鲜窍嗷ト哂嗯c獨(dú)立的。SMSFR-V0的控制棒與燃料組件比為12.5%,略高于已有快堆設(shè)計(jì)。由于用于控制棒驅(qū)動(dòng)機(jī)構(gòu)的剩余安裝空間非常有限,因此難以增加SMSFR-V0的控制棒數(shù)量。
該堆芯采用與ASTRID CFV-1500相似的控制棒使用策略[8],即同時(shí)使用CSD與DSD進(jìn)行燃耗反應(yīng)性損失補(bǔ)償、功率分布調(diào)節(jié)與功率水平調(diào)節(jié)。在緊急停堆時(shí),CSD與DSD同時(shí)作為安全停堆棒。該策略使用盡可能多的控制棒進(jìn)行反應(yīng)性控制,因此減少了單個(gè)棒中所平衡的剩余反應(yīng)性,從而減少了控制棒意外提出(CRW)事故的影響。
本文使用確定論中子輸運(yùn)軟件APOLLO3進(jìn)行反應(yīng)堆物理數(shù)值計(jì)算。APOLLO3[16]是法國(guó)CEA、EDF與Framatome聯(lián)合開(kāi)發(fā)的適用于多種堆型的高并行度確定論中子輸運(yùn)軟件。APOLLO3針對(duì)快堆的中子控制棒的計(jì)算已經(jīng)過(guò)詳細(xì)對(duì)標(biāo)與驗(yàn)證[17-18]。
SMSFR-V0主要參數(shù)列于表1。燃料換料周期為1 875 EFPD(當(dāng)量滿功率天),共分5批無(wú)倒料燃料。卸料平均燃耗為120 GW·d/t,峰值燃耗為163 GW·d/t。裂變區(qū)平均功率密度為240 W/cm3,對(duì)應(yīng)于燃料棒的平均線功率密度為294 W/cm。壽期初(BOEC)和壽期末(EOEC)的功率峰因子分別為1.53和1.40。燃耗反應(yīng)性損失為-8.7 pcm/EFPD,單個(gè)換料周期燃耗損失為3 257 pcm,通過(guò)控制棒逐步提出來(lái)補(bǔ)償該燃耗反應(yīng)性損失。
表1 SMSFR-V0主要參數(shù)Table 1 Main parameter of SMSFR-V0
SMSFR-V0的反應(yīng)性控制需求示于圖2。除燃耗反應(yīng)性損失外,為實(shí)際運(yùn)行的靈活性與計(jì)算的不確定性需預(yù)留600 pcm剩余反應(yīng)性裕度。因此部分控制棒的插入需具備3 857 pcm的反應(yīng)性價(jià)值以實(shí)現(xiàn)反應(yīng)堆從BOEC到EOEC的臨界運(yùn)行。此外,停堆后239Np將繼續(xù)β衰變?yōu)?39Pu,從而小幅提高堆芯反應(yīng)性??刂瓢粜枰WC足夠的停堆深度,本文參考超鳳凰堆設(shè)計(jì)準(zhǔn)則,保留10 $裕度以確保發(fā)生換料事故時(shí)堆芯的次臨界安全。因此,所有控制棒的插入需具備大于8 607 pcm的反應(yīng)性價(jià)值,以實(shí)現(xiàn)所有工況下的停堆能力,并為事故工況保留裕度。
圖2 SMSFR-V0中的反應(yīng)性控制Fig.2 Reactivity control in SMSFR-V0
一方面,控制棒的插入降低了堆芯頂部的功率,從而增加了堆芯中央的功率峰值。線功率密度峰值隨控制棒插入深度的變化如圖3所示。另一方面,如果采用較小的壽期初臨界插入深度來(lái)平衡堆芯剩余反應(yīng)性,則會(huì)增加對(duì)吸收體的吸收能力需求。若采用90%10B的B4C為吸收劑,則控制棒的安全使用壽命將非常有限,甚至短于單個(gè)燃料循環(huán)長(zhǎng)度[19]。綜合以上兩方面考慮,最終確定約30 cm為壽期初的臨界插入深度,對(duì)應(yīng)平均線功率密度為294 W/cm、線功率密度峰值為455 W/cm、功率峰因子為1.55。
圖3 不同控制棒插入深度下的峰值線功率密度Fig.3 Peak linear power density at different insertion depth of control rods
SMSFR-V0的10B選擇曲線示于圖4a。燃耗反應(yīng)性損失是反應(yīng)性控制所需控制棒價(jià)值的主要變量,圖4a表示為滿足不同燃耗反應(yīng)性損失控制棒所需要的10B富集度。SMSFR-V0的壽期燃耗反應(yīng)性損失為-3 257 pcm,因此所需的10B富集度約為60%。前期研究所得新型吸收材料與不同10B富集度B4C的吸收能力對(duì)應(yīng)曲線如圖4b[13]所示(鉿、銪與釓為天然豐度;Mod.為氫化鋯慢化劑)。圖4b表明,HfH1.62或90%10B富集度的HfB2(簡(jiǎn)稱HfB2-90)的吸收能力與60%10B富集度的B4C(簡(jiǎn)稱B4C-60)的吸收能力相當(dāng),有可能滿足對(duì)反應(yīng)性控制系統(tǒng)的需求。因此,本文針對(duì)SMSFR-V0主要對(duì)比B4C-60、HfB2-90與HfH1.62。
圖4b是基于3 600 MWth大型鈉冷快堆V2B的分析結(jié)果[13]。V2B與SMSFR-V0具有相近能譜,因此該結(jié)果在小堆中初步可用。圖4組合了10B選擇曲線與新型吸收劑吸收能力曲線,可作為一種在不同設(shè)計(jì)下初步搜尋合適吸收劑的方法。
圖4 SMSFR-V0中10B選擇曲線(a)和控制棒優(yōu)化設(shè)計(jì)的吸收能力(b)Fig.4 10B content selection curve in SMSFR-V0 (a) and absorption ability of optimized design of control rod (b)
圖6 控制棒體軸向示意圖(a)和組件軸向示意圖(b)Fig.6 Axial layout of control rod body (a) and control rod assembly (b)
SMSFR-V0中使用的控制棒如圖5、6所示。圖5中未顯示繞線和穩(wěn)流結(jié)構(gòu)。吸收體高度與燃料活性區(qū)相同,均為85 cm??刂瓢趔w為圓形,與六角形套管的間隙為5 mm,以保證自由旋轉(zhuǎn),避免卡棒。該間距是基于ASTRID CFV-1500的初步設(shè)計(jì),具體數(shù)值有待進(jìn)一步調(diào)整??刂瓢粲?9根吸收劑元件棒組成。B4C、HfB2和HfH1.62與鈉均有較好相容性,因此吸收劑棒和吸收劑與包殼間介質(zhì)可為鈉或氦氣。但由于10B主要通過(guò)(n,α)反應(yīng)吸收中子,B4C與HfB2需要設(shè)置氣孔排氣,其介質(zhì)必須為鈉。HfH1.62沒(méi)有氣體釋放的問(wèn)題,但氦氣介質(zhì)會(huì)導(dǎo)致較大溫差,而HfH1.62的溫度裕度較低,因此HfH1.62采用鈉為吸收劑-包殼介質(zhì),以增強(qiáng)導(dǎo)熱性。
圖5 控制棒徑向布置Fig.5 Radial layout of control rod
控制棒由燃料活性區(qū)頂部向下插入堆芯,假設(shè)所有控制棒同時(shí)插入到同一深度,插入深度定義為控制棒底部與活性區(qū)頂部的距離。SMSFR-V0剩余反應(yīng)性變化和臨界插入深度示于圖7??梢?jiàn),堆芯剩余反應(yīng)性隨反應(yīng)堆運(yùn)行不斷下降,壽期初控制棒插入一定深度,運(yùn)行過(guò)程中控制棒逐步提出以補(bǔ)償燃料燃耗反應(yīng)性損失,因此能維持反應(yīng)堆臨界狀態(tài)的臨界插入深度隨時(shí)間不斷減小。在小型快堆中,該控制棒提出過(guò)程對(duì)堆芯功率分布與控制棒中子注量的影響較為明顯,因此燃耗計(jì)算中考慮了該控制棒提出的過(guò)程。但需要指出的是在燃耗過(guò)程中控制棒的吸收能力不斷降低,因此其臨界插入深度不斷加深,本文計(jì)算僅采用壽期初的臨界插入深度,該假設(shè)可能導(dǎo)致控制棒的中子注量計(jì)算值較實(shí)際略低。
圖7 SMSFR-V0剩余反應(yīng)性變化和臨界插入深度Fig.7 Residual reactivity swing and critical position of control rod in SMSFR-V0
燃料在運(yùn)行過(guò)程中每375 EFPD進(jìn)行1次1/5換料,而控制棒不更換直至經(jīng)受5個(gè)周期輻照,即1 875 EFPD。不同控制棒的反應(yīng)性價(jià)值列于表2。該結(jié)果符合圖4預(yù)測(cè)情況,即B4C-60與HfB2-90控制棒反應(yīng)性價(jià)值接近,略高于HfH1.62控制棒。3種吸收劑在輻照后的反應(yīng)性價(jià)值損失相當(dāng),相比于整根控制棒,其臨界插入部分中子注量較高,反應(yīng)性價(jià)值損失較大。
表2 SMSFR-V0控制棒反應(yīng)性價(jià)值Table 2 Reactivity worth of different control rods in SMSFR-V0
控制棒對(duì)堆芯功率分布有重要影響。控制棒是堆芯中的強(qiáng)吸收體,控制棒的插入會(huì)削弱其周?chē)闹凶油棵芏龋纬删植抗β使?。SMSFR-V0不同控制棒的堆芯線功率密度峰值列于表3??煽闯?,由于控制棒的提出,線功率峰值從BOEC到EOEC有所降低。由于HfH1.62具有慢化作用,因此相對(duì)于B4C控制棒,HfH1.62減小了吸收劑對(duì)周?chē)剂瞎β实南魅踝饔?,使控制棒插入時(shí)的堆芯功率分布更加均勻,如圖8所示。因此,使用HfH1.62控制棒的堆芯線功率峰低于使用其他控制棒的堆芯。
表3 SMSFR-V0不同控制棒的堆芯線功率密度峰值Table 3 Peak linear power density in SMSFR-V0 with different control rod designs
圖8 BOEC時(shí)HfH1.62和B4C-60控制棒的SMSFR-V0堆芯線功率密度間差異Fig.8 Relative difference on core power density between cores with HfH1.62 rod and B4C-60 rod at BOEC
表4對(duì)比了使用不同類(lèi)型控制棒堆芯的鈉空效應(yīng)和多普勒效應(yīng)??煽闯?,不同堆芯全鈉空泡反應(yīng)性價(jià)值接近,由于該小型模塊化快堆采用頂部鈉反射層設(shè)計(jì)及其本身中子泄漏率較大,因此在平衡壽期內(nèi)全鈉空泡反應(yīng)性價(jià)值略偏負(fù)。HfH1.62具有慢化作用,因此其多普勒常數(shù)增強(qiáng)較為顯著。
表4 控制棒對(duì)SMSFR-V0反應(yīng)性反饋系數(shù)的影響Table 4 Influence of control rod on feedback coefficient of SMSFR-V0
使用不同吸收體的控制棒最高溫度變化如圖9所示。對(duì)于硼基吸收體,溫度在前2個(gè)循環(huán)中明顯升高,因?yàn)樵谳椪臻_(kāi)始時(shí)導(dǎo)熱系數(shù)的下降非常明顯。溫度峰值出現(xiàn)在每個(gè)循環(huán)的壽期初,壽期內(nèi)控制棒逐步提出,中子通量密度降低,從而導(dǎo)致功率密度及溫度下降。由于缺少HfH1.62熱導(dǎo)率隨其燃耗變化的數(shù)據(jù),因此采用熱導(dǎo)率不變假設(shè)。由于吸收反應(yīng)率逐漸降低,HfH1.62控制棒的溫度隨時(shí)間緩慢下降。
圖9 SMSFR-V0控制棒最高溫度變化Fig.9 Evolution of maximal temperature in control rod of SMSFR-V0
B4C的熔點(diǎn)約為2 350 ℃,但由于包殼的碳化風(fēng)險(xiǎn),若采用鳳凰堆與超鳳凰堆1 200 ℃的保守限值,B4C的安全使用壽命僅為1個(gè)循環(huán)。較小的吸收劑元件尺寸可用于降低中心溫度,但其使用壽命將受到最大燃耗的限制。由于具有較高的熔點(diǎn)和熱導(dǎo)率,HfB2在5個(gè)輻照循環(huán)中溫度裕度始終大于2 200 ℃。HfH1.62在高溫下有氫氣解離風(fēng)險(xiǎn),其溫度裕度僅約為250 ℃,因此需要進(jìn)一步評(píng)估HfH1.62的安全性。
使用不同吸收體的控制棒燃耗如圖10所示。吸收體的燃耗定義為其單位體積的累計(jì)吸收反應(yīng)率。對(duì)于不同類(lèi)型控制棒,有最大燃耗限值以確保輻照后吸收體與包殼的安全。已有實(shí)驗(yàn)經(jīng)驗(yàn)與實(shí)驗(yàn)數(shù)據(jù)表明,B4C的最大燃耗為2.1×1022cm-3[20],而HfB2的輻照膨脹遠(yuǎn)小于B4C,其燃耗限值約為B4C的3倍[21]。HfH1.62的燃耗限值尚待確定。這些控制棒的吸收能力接近,因此其燃耗相近。B4C在第4個(gè)循環(huán)結(jié)束時(shí)達(dá)到限值,而HfB2仍具有較大裕度。
圖10 SMSFR-V0控制棒峰值燃耗變化Fig.10 Evolution of peak burnup in control rod of SMSFR-V0
由于小型快堆的中子泄漏率較大、增殖能力偏弱、燃耗反應(yīng)性損失較大,控制棒的過(guò)度插入將增加堆芯的功率峰因子。同時(shí)快堆通常采用獨(dú)立控制棒組件形式,導(dǎo)致小堆中控制棒的安裝數(shù)量有限。因此,小型快堆需要高10B富集度的B4C進(jìn)行反應(yīng)性控制。但B4C存在吸收劑燃耗深、功率密度高而導(dǎo)熱能力受輻照削弱大等問(wèn)題,傳統(tǒng)吸收劑B4C的安全使用壽命有限。
本文具體分析了320 MWth小型模塊化鈉冷快堆SMSFR-V0的反應(yīng)性控制要求,并結(jié)合前期新型吸收劑研究,初步篩選了90%10B富集度的HfB2與HfH1.62作為吸收劑,在吸收能力上可替代60%10B富集度的B4C。對(duì)含有這3種吸收劑的控制棒進(jìn)行了綜合對(duì)比,包括反應(yīng)性價(jià)值、堆芯功率分布、堆芯反應(yīng)性反饋系數(shù)、控制棒溫度裕度與吸收劑燃耗深度。與傳統(tǒng)以B4C作為吸收劑的控制棒相比,HfB2有更高的安全裕度和更長(zhǎng)的安全使用壽命,HfH1.62控制棒對(duì)功率分布略有改善,但其高溫氫氣解離問(wèn)題有待進(jìn)一步研究。
含鉿控制棒使用與B4C控制棒相同的吸收劑元件尺寸,這有助于直觀對(duì)比吸收劑的性能,減少對(duì)控制棒其他特性的影響,然而吸收劑元件大小有待進(jìn)一步根據(jù)吸收劑材料特性進(jìn)行優(yōu)化,如HfB2的導(dǎo)熱性較好而輻照膨脹較小,可進(jìn)一步擴(kuò)大其元件尺寸,減小吸收劑-包殼間距以增加控制棒內(nèi)吸收劑比例。此外一些新型設(shè)計(jì)有待探索,如局部慢化劑與含鉿吸收劑的結(jié)合等。