VVER-1000反應(yīng)堆控制棒延壽中子物理學(xué)論證分析

蔣朱敏 王晨琳 李載鵬 于穎銳 王丹 劉同先 王金雨 鐘旻霄

摘 ? 要：田灣VVER-1000反應(yīng)堆控制棒的初始設(shè)計(jì)壽命為10年。通過開展必要的論證工作，將控制棒的使用壽命延長至15年，可以在很大程度上降低控制棒的采購成本和提高經(jīng)濟(jì)性。本文從中子物理學(xué)的角度對控制棒延壽至15年的可行性進(jìn)行了分析評價(jià)。分別從控制棒的鈦酸鏑吸收體中子注量、碳化硼吸收體的10B燃耗份額、下部端塞快中子注量以及控制棒組的反應(yīng)性價(jià)值變化這四個(gè)方面，進(jìn)行了詳細(xì)的計(jì)算分析。結(jié)果表明，田灣VVER-1000反應(yīng)堆控制棒組件在反應(yīng)堆上使用15年后，中子物理特性仍然滿足其初始設(shè)計(jì)的相關(guān)技術(shù)要求。

關(guān)鍵詞：控制棒 ?延壽 ?中子物理分析 ?VVER-1000反應(yīng)堆

中圖分類號：TL48 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)識碼：A ? ? ? ? ? ? ? ? ? ? ? ?文章編號：1674-098X（2019）02（a）-0105-04

Abstract： The designed lifetime of the control rods of VVER-1000 nuclear reactor is 10 years. By extending the lifetime of the control rods to 15 years， the cost of purchasing them will be significantly reduced. For this purpose， necessary assessment needs to be done. In this paper， neutron physics analysis has been carried out， on the following 4 aspects： the neutron fluent of the Dy2O3·TiO2 absorber， the burned percentage of 10B in the B4C absorber， the neutron fluent of the lower plug， and the reactivity change of the control rod banks. The analysis results show that， after operating for a time of 15 years， the neutron-physic characteristics of control rods of the VVER-1000 nuclear reactor still meet the original designed limits.

Key Words： Control rod; Lifetime extending; Neutron physics analysis; VVER-1000 nuclear reactor

核電站反應(yīng)堆控制棒起補(bǔ)償和調(diào)節(jié)堆芯反應(yīng)性和緊急落棒停堆的作用。反應(yīng)堆運(yùn)行時(shí)，控制棒持續(xù)不斷地受到來自堆芯活性區(qū)內(nèi)及活性區(qū)上方的中子輻照。一方面控制棒中的有效吸收成分不斷減少，在達(dá)到一定的服役時(shí)間后，控制棒中子吸收能力降低，反應(yīng)性價(jià)值降低;另一方面控制棒吸收體及包殼材料的抗輻照性能逐漸變差，并可能出現(xiàn)輻照腫脹、包殼變形開裂等現(xiàn)象。

田灣核電站控制棒的初始設(shè)計(jì)壽命為10年。本文通過論證，從中子物理學(xué)的角度評價(jià)控制棒延壽至15年的可行性，以達(dá)到延長控制棒壽命，降低成本和提高經(jīng)濟(jì)性的目的。

田灣1、2號機(jī)組均在第10循環(huán)開始前更換了新的控制棒組件，為此本文的論證分析從第10循環(huán)開始。

1 ?VVER-1000控制棒簡介

田灣核電站1&2號機(jī)組為VVER 1000型壓水堆，其堆芯共有103束控制棒，分為10組。反應(yīng)堆正常運(yùn)行時(shí)，第10組棒位于堆芯高度的90%附近，用于穩(wěn)定堆芯功率，稱為工作棒組;其余棒組提出到堆芯外，主要用于反應(yīng)堆的應(yīng)急停堆保護(hù)，稱為停堆棒組。

每一束控制棒的結(jié)構(gòu)和材料相同。徑向上，采用42XHM合金制成的圓筒形包殼，將圓柱形的中子吸收體包裹在其內(nèi)部。吸收體在軸向采用分段設(shè)計(jì)，全長3.5m，其中上部3.2m的吸收體材料為B4C（碳化硼），下端0.3m的吸收體材料為Dy2O3·TiO2（鈦酸鏑）。包殼的兩端采用42XHM合金制成的端塞進(jìn)行封裝。

碳化硼吸收體中的有效中子吸收成分為10B。10B吸收中子后生成的產(chǎn)物中子吸收截面都很小，因此碳化硼吸收體的中子吸收能力隨燃耗下降較快。另外，反應(yīng)生成的氦氣逐漸累積后可能導(dǎo)致吸收體腫脹及包殼變形。

鈦酸鏑吸收體中主要的吸收材料為鏑（Dy），它包含多種同位素：158Dy、160Dy、161Dy、162Dy、163Dy、164Dy。各同位素的吸收截面均小于10B，因此鈦酸鏑吸收體的反應(yīng)性價(jià)值較小[1]。Dy吸收中子后的反應(yīng)性價(jià)值逐步減小，但相對碳化硼吸收體而言，其反應(yīng)性隨燃耗的下降速度更慢。

2 ?計(jì)算程序

計(jì)算主要采用俄羅斯庫爾恰托夫研究院開發(fā)的KASKAD程序包進(jìn)行。具體包括以下程序：

TVS-M程序：準(zhǔn)備并擬合BIPR-7A和PERMAK-A程序所需要的燃料組件、燃料棒、含釓燃料棒、吸收體元件以及可燃毒物棒的少群中子截面參數(shù)，這些少群中子截面參數(shù)是反應(yīng)堆狀態(tài)和燃料燃耗的函數(shù)。

BIPR-7A程序：計(jì)算臨界參數(shù)、反應(yīng)性系數(shù)、反應(yīng)性效應(yīng)、控制棒微積分價(jià)值、堆芯三維功率分布，提供燃耗計(jì)算和換料計(jì)算，以及計(jì)算氙-135和釤-149的瞬態(tài)變化。

PERMAK-A程序：為多層兩維4群或6群（熱中子分3群）的細(xì)網(wǎng)擴(kuò)散程序，用于計(jì)算燃料棒功率、燃耗以及在控制棒組移動和功率變化狀態(tài)下的線功率密度變化信息。

除了上述程序外，還使用堆用蒙卡分析程序RMC（Reactor Monte Carlo code）進(jìn)行控制棒組件燃耗計(jì)算及頂部反射層內(nèi)的注量率分布計(jì)算。RMC是由清華大學(xué)工程物理系核能科學(xué)與工程管理研究所反應(yīng)堆工程計(jì)算分析實(shí)驗(yàn)室自主研發(fā)的、用于反應(yīng)堆堆芯計(jì)算分析的三維輸運(yùn)蒙卡程序。其采用的蒙卡輸運(yùn)計(jì)算方法適用于任意幾何的堆芯，同時(shí)內(nèi)嵌燃耗計(jì)算功能。

3 ?延壽后的技術(shù)要求

控制棒的包殼和端塞的材料可以視作是不可燃耗的，而鈦酸鏑吸收體的反應(yīng)性隨燃耗下降較慢。對于這兩種材料主要考慮其抗輻照性能。通常隨著其受到的累計(jì)中子注量的增加材料的抗輻照性能逐步降低，因此需對這兩種材料所受的累計(jì)中子注量進(jìn)行限制。

對于碳化硼吸收體，其10B的消耗導(dǎo)致反應(yīng)性價(jià)值下降較快，且并導(dǎo)致氦氣的積累，因此需對10B燃耗份額進(jìn)行限制。

參照田灣1&2號機(jī)組控制棒組件初始壽命的論證依據(jù)及控制棒價(jià)值驗(yàn)收準(zhǔn)則，延壽后控制棒組件仍需滿足與初始設(shè)計(jì)相同的技術(shù)要求，即：

（1）鈦酸鏑受到的快中子（能量> 0.1 MeV）注量< 4×1022 n/cm2，熱中子注量< 7×1022 n/cm2;

（2）10B燃耗份額< 45%;

（3）42XHM合金（包殼和下部端塞）受到的快中子（能量> 0.1 MeV）注量< 3.4×1022 n/cm2;

（4）控制棒價(jià)值滿足事故分析確定的工作棒組價(jià)值及應(yīng)急停堆保護(hù)價(jià)值要求，控制棒組價(jià)值理論計(jì)算值與實(shí)測值的偏差不超過±20%。

4 ?中子物理論證

4.1 鈦酸鏑中子注量分析

工作棒組（第10組）在反應(yīng)堆正常運(yùn)行時(shí)長期插入堆芯，因此其受到的中子注量明顯高于其余棒組。為盡量減少單個(gè)控制棒束受到的中子注量，可在換料大修期間更換用作工作棒組的控制棒束。由于工作棒組僅包含6個(gè)控制棒束，而全堆芯共計(jì)103個(gè)控制棒束，因此可以確保每個(gè)控制棒束在15年的使用壽命內(nèi)，僅有一個(gè)燃料循環(huán)用作工作棒組。

以2號機(jī)組為例，參照其燃料管理報(bào)告及長燃料循環(huán)論證確定后續(xù)循環(huán)堆芯裝載方案，挑選注量最高的一個(gè)循環(huán)將某個(gè)控制棒束用作工作棒組，在其余循環(huán)將其用作停堆棒組。

需要說明的是，目前田灣核電站的燃料循環(huán)按照18個(gè)月的換料周期進(jìn)行設(shè)計(jì)——每個(gè)循環(huán)的循環(huán)長度約為480等效滿功率天（EFPD，在滿功率情況下運(yùn)行的天數(shù)），考慮到運(yùn)行結(jié)束之后的換料大修時(shí)間，及節(jié)假日等因素造成的降功率運(yùn)行，每個(gè)燃料循環(huán)的實(shí)際長度約為1.5年。為此，控制棒15年的工作壽命對應(yīng)為10個(gè)燃料循環(huán)。

計(jì)算時(shí)，保守地假定工作棒組位于80%堆芯高度處。使用KASKAD程序進(jìn)行額定功率下的各燃料循環(huán)的堆芯燃耗計(jì)算，可得到各循環(huán)不同燃耗步下的堆芯快群和熱群中子注量率。將其對時(shí)間積分，即可得到工作棒組出的快群和熱群中子注量[2-4]。

停堆棒組正常情況下全部提出堆芯，其下端位于堆芯上方的水反射層內(nèi)。理論分析表明，反射層內(nèi)的中子注量率明顯低于堆芯內(nèi)部，但考慮到控制棒在此處放置時(shí)間最長可達(dá)15年，其受到的累計(jì)注量不可忽略不計(jì)。

KASKAD程序無法直接模擬計(jì)算反射層內(nèi)的中子注量率分布。為此采用以下方式間接模擬：

（1）利用KASKAD程序計(jì)算得到的堆芯三維注量率分布，并保守的統(tǒng)計(jì)出堆芯頂部的最大中子注量;

（2）采用RMC程序進(jìn)行單個(gè)燃料組件的三維輸運(yùn)計(jì)算，得到堆芯頂部區(qū)域（包括堆芯活性區(qū)最上面若干區(qū)域及頂部反射層）的相對注量率分布，并由此得到停堆棒組下端處與堆芯頂部處的中子注量率的比例系數(shù);

（3）利用第1）步計(jì)算的堆芯頂部的最大中子注量，乘以第2步）得到的比例系數(shù)，得到堆棒組下端處的中子注量。

圖1給出了RMC計(jì)算得到的堆芯頂部區(qū)域的快群和熱群相對注量率分布。該分布用于表征頂部反射層內(nèi)的中子注量率與堆芯頂部中子注量率的相對關(guān)系，因此分別對熱群和快群各自乘以歸一化系數(shù)，使得堆芯頂部的中子注量率為1。圖中A點(diǎn)代表鈦酸鏑底部位置，B點(diǎn)和C點(diǎn)分別代表碳化硼吸收體底部和下端塞底部。

由圖可見，停堆棒組鈦酸鏑吸收體處的熱群和快群相對注重注量率分別為0.854和0.135?？紤]適當(dāng)?shù)谋Ｊ卦Ａ?，分別取為1.0和0.2進(jìn)行后續(xù)計(jì)算。

圖2和圖3分別給出了2號機(jī)組后續(xù)10個(gè)燃料循環(huán)的鈦酸鏑吸收體累計(jì)快中子注量和熱中子注量計(jì)算結(jié)果。由圖可見，鈦酸鏑吸收體受到的快中子累計(jì)注量分別為1.86×1022n/（cm2），小于快中子注量限值4×1022n/（cm2）;熱中子累計(jì)注量為0.80×1022n/（cm2），小于熱中子注量限值7×1022n/（cm2）?？?、熱中子注量均滿足技術(shù)要求。

4.2 10B燃耗份額分析

10B燃耗份額的計(jì)算分為如下3步：

（1）采用RMC程序進(jìn)行插入控制棒的單個(gè)組件徑向二維輸運(yùn)、燃耗計(jì)算得到10B燃耗份額隨控制棒柵元中子注量的變化關(guān)系。

（2）采用鈦酸鏑中子注量分析類似的方法，計(jì)算得到B4C吸收體的累計(jì)中子注量。

（3）根據(jù)第（2）步計(jì)算得到的中子注量，從第（1）步計(jì)算結(jié)果中插值得到10B燃耗份額。

由圖1可見，停堆棒組B4C吸收體處（ULS + 30 cm）的熱群和快群相對通量分別為0.008和0.003?？紤]適當(dāng)?shù)谋Ｊ卦Ａ?，熱群和快群均?.1進(jìn)行后續(xù)碳化硼中子注量計(jì)算。

圖4～圖6給出了相關(guān)計(jì)算結(jié)果。由圖6可見，延壽至15年后，10B累計(jì)燃耗份額為34.27%，小于限值45%，滿足技術(shù)要求。

4.3 下部端塞快中子注量分析

控制棒束的包殼和端塞均采用42XHM合金制成。其中下部端塞插入堆芯的距離最深，用作停堆棒組時(shí)距離堆芯活性區(qū)的距離最近。為此，下部端塞收到的快中子注量是所有42XHM合金中最多的。

下部端塞中子注量的計(jì)算方法與鈦酸鏑類似。所不同的是下部端塞位與鈦酸鏑吸收體下方，受到的中子注量大于鈦酸鏑吸收體。

圖7給出了下部端塞受到的快中子注量計(jì)算結(jié)果。由圖可見，延壽至15年后，下部端塞收受到的快中子累計(jì)注量為2.28×1022n/（cm2），小于快中子注量限值3.4×1022n/（cm2），滿足技術(shù)要求。

4.4 控制棒組反應(yīng)性價(jià)值分析

為了估算控制棒燃耗對其反應(yīng)性價(jià)值的影響，對于碳化硼和鈦酸鏑兩種吸收體材料，采用TVS-M程序分別計(jì)算其燃耗和不燃耗這兩種情況下，一束控制棒插入組件中引入的負(fù)反應(yīng)性大小。相關(guān)計(jì)算結(jié)果見表1和表2。

仍然保守地假定工作棒組位于80%堆芯高度處。在此情況下，計(jì)算可得經(jīng)歷1個(gè)平衡循環(huán)后，工作棒組鈦酸鏑所在堆芯高度處組件的平均燃耗約為19 MWd/kgU，碳化硼所在堆芯高度處組件的平均燃耗約為12 MWd/kgU。由表1和表2可知插入組件內(nèi)的吸收體燃耗導(dǎo)致的反應(yīng)性價(jià)值損失分別約為1584 pcm和1398 pcm，約占初始價(jià)值的11%和6%。

對于停堆棒組，由圖18可知，其碳化硼吸收體在9個(gè)循環(huán)累計(jì)的中子注量約為用作工作棒組1個(gè)循環(huán)累計(jì)中子注量的40%，由此估算可得碳化硼吸收體在整個(gè)控制棒壽期內(nèi)的反應(yīng)性損失約為6%×1.4≈8.4%。對于鈦酸鏑吸收體，其用作停堆棒組時(shí)熱中子注量仍然較高，因此其燃耗份額會大于碳化硼吸收體。

上述僅是針對控制棒下端燃耗最大的部分的計(jì)算和分析結(jié)果，在評估其對整個(gè)棒組反應(yīng)性價(jià)值的影響時(shí)，還需考慮以下因素：

（1）鈦酸鏑吸收體僅占控制棒總長度的約8.6%;

（2）碳化硼吸收體僅在用作工作棒組時(shí)有很小一部分插入堆芯（以位于80%堆芯高度處為例，約為44.5cm，占總長度的12.7%），這部分的反應(yīng)性損失如前述分析（損失約8.4%）。

（3）碳化硼吸收體最下端44.5cm以上的部分僅在用作工作棒組時(shí)，于頂部反射層中吸收中子，由圖1可知，這部分吸收體的燃耗明顯小于插入堆芯的部分（保守地假定損失的反應(yīng)性價(jià)值為插入堆芯部分的10%，即0.6%）;

（4）碳化硼吸收體其余部分幾乎沒有燃耗;

（5）棒組積分價(jià)值是指控制棒完全插入堆芯底部時(shí)引入的總反應(yīng)性，此時(shí)棒束燃耗較多的部分靠近堆芯底部，其附近中子價(jià)值較低，使得吸收體燃耗對于整個(gè)棒束積分價(jià)值的影響進(jìn)一步減小。

考慮靠以上因素，做如下保守假設(shè)：

（1）控制棒束各部分的反應(yīng)性價(jià)值對棒組積分價(jià)值的貢獻(xiàn)與其長度成正比，即假設(shè)棒束各部分微分價(jià)值相等，如此假設(shè)增大了控制棒下端的微分價(jià)值，考慮到燃耗主要發(fā)生在下端，是一種保守的假設(shè);

（2）假定鈦酸鏑吸收體燃耗15年后反應(yīng)性價(jià)值全部損失;

（3）假定碳化硼吸收體最下端44.5cm燃耗15年后反應(yīng)性價(jià)值變?yōu)槌跏贾档?0%;

（4）用作工作棒組時(shí)，堆芯反射層中的碳化硼吸收體燃耗1個(gè)循環(huán)（15年壽期內(nèi)只有1個(gè)循環(huán)用作工作棒組）后反應(yīng)性價(jià)值變?yōu)槌跏贾档?00%-0.6%=99.4%，且假定其長度為40cm;

（5）碳化硼吸收體其余部分15年后的反應(yīng)性價(jià)值不變。

根據(jù)上述假設(shè)，可保守地估算得使用15年后，控制棒組件剩余積分價(jià)值占其初始值的比例為：

可見燃耗導(dǎo)致的控制棒價(jià)值變化小于其計(jì)算誤差范圍（±20%）。

鑒于上述計(jì)算結(jié)果引入了較多保守假設(shè)，為了進(jìn)一步評估實(shí)際運(yùn)行過程中，控制棒燃耗對其總體積分價(jià)值的影響，表3給出了1、2號機(jī)組第2～9循環(huán)應(yīng)急保護(hù)價(jià)值和工作棒組積分價(jià)值的實(shí)際測量結(jié)果，表中計(jì)算值中未考慮控制棒燃耗。由表可知，到第9循環(huán)壽期初為止，應(yīng)急保護(hù)價(jià)值均遠(yuǎn)大于安全限值7100pcm的要求;棒價(jià)值測量值與計(jì)算值的相對偏差均在±10.8%以內(nèi)，滿足設(shè)計(jì)準(zhǔn)則±20%的要求;偏差沒有明顯增大的趨勢，說明控制棒燃耗不會明顯增大控制棒積分價(jià)值的計(jì)算誤差。

5 ?結(jié)語

從中子物理學(xué)的角度，計(jì)算分析了田灣核電站1&2號機(jī)組控制棒組件延壽后，其成分及反應(yīng)性價(jià)值隨其在堆芯輻照時(shí)間的變化情況。參照其控制棒組件初始設(shè)計(jì)技術(shù)要求，分別從鈦酸鏑吸收體和下端塞結(jié)構(gòu)材料受到的累計(jì)中子注量、碳化硼吸收體的10B燃耗份額、控制棒組積分價(jià)值變化等方面進(jìn)行了論證分析。結(jié)果表明，田灣VVER-1000反應(yīng)堆控制棒組件在反應(yīng)堆上使用15年后，其中子物理特性仍然滿足其初始設(shè)計(jì)的相關(guān)技術(shù)要求。

參考文獻(xiàn)

[1] V. D. Risovany， etc. Dysprosium titanate as an absorber material for control rods. Joural of Nuclear materials， 2000.

[2] OmidNoori-Kalkhoran， etc. Development of external coupling for calculation of the control rod worth in terms of burn-up for a WWER-1000 nuclear reactor. Nuclear Engineering and Design， 2016.

[3] Fausto Franceschini， etc. ?Development of a control rod depletion methodology for the Westinghouse NEXUS system. Progress in Nuclear Energy， 2013.

[4] Park， etc. Depletion analysis of control absorber in a small research reactor. Annals of Nuclear Energy， 2013.