于成波，李 敏，廖 路，劉旭東，唐 豪，徐 濤，高燊甫

(中國核動力研究設(shè)計院,四川 成都 610041)

對于商用反應(yīng)堆而言，由于沒有特別的要求，完全可以避開碘坑內(nèi)啟動，但是對于研究堆，受輻照任務(wù)限制以及一些復(fù)雜的工況要求，會存在碘坑內(nèi)啟動的情況。因此研究并估算HFETR碘坑后沿啟動反應(yīng)堆的臨界棒位對提高反應(yīng)堆的安全性和經(jīng)濟性有重要的理論意義和應(yīng)用價值。

在碘坑后沿啟動反應(yīng)堆時,由于中子通量密度突然增加，氙吸收中子后被大量消耗，堆內(nèi)迅速釋放出正反應(yīng)性。這時，自動棒雖會自動跟蹤下降，但跟蹤范圍有限，尤其是當(dāng)在碘坑底部啟動核反應(yīng)堆時,會在極短時間內(nèi)產(chǎn)生很大的正反應(yīng)性，自動控制棒跟蹤完全有可能補償不了，有可能造成正反應(yīng)性引入事故。碘坑內(nèi)啟動是影響反應(yīng)堆運行安全的一個重要問題。

目前，國內(nèi)外對碘坑曲線及臨界棒位的研究鮮見報道。曾道桂等人對西安脈沖堆裂變產(chǎn)物中毒進行了深入地分析研究。利用脈沖堆物理計算程序計算了西安脈沖堆在不同功率下停堆前后氙(135Xe)和釤(149Sm)反應(yīng)性變化，得到氙毒和碘坑的計算值，并與實測值進行比較[1]。張帆等人進行了船用堆停堆后碘坑情況的仿真，并和點堆模型計算數(shù)據(jù)進行了對比，結(jié)果表明用三維雙群堆芯物理計算模型對碘坑下啟動進行實時仿真，比采用點堆模型更能滿足工程仿真的精度要求，適用于碘坑仿真計算[2]。王皓等人使用數(shù)值方法并通過程序的改進完善，在高通量工程試驗堆(HFETR)某爐次裝載方案的基礎(chǔ)上進行不同工況下的氙毒效應(yīng)計算研究，得到HFETR以不同功率運行相同積分功率、以不同功率運行不同積分功率后緊急停堆的氙毒效應(yīng)曲線[3]。黎浩峰研究了反應(yīng)堆碘坑內(nèi)的啟動問題,導(dǎo)出了熱堆在碘坑內(nèi)啟動時氙毒反應(yīng)性、碘和氙濃度隨時間變化的解析表達式。通過對所得解析解的計算分析，對碘坑內(nèi)啟動問題有了更全面、本質(zhì)了解[4]。曹寅等人根據(jù)低濃化后各種堆芯參數(shù)的改變，分析計算了氙(135Xe)和釤(149Sm)這2種反應(yīng)堆中最重要裂變產(chǎn)物的變化對堆芯反應(yīng)性產(chǎn)生的影響，并給出額定工況運行與停堆后2種情況下的反應(yīng)性變化曲線，使反應(yīng)堆操縱人員了解毒物反應(yīng)性的變化規(guī)律，以便于在實際運行中應(yīng)對工況變化，保證反應(yīng)堆正常安全運行[5]。本文以HFETR某爐段進行研究，首先推導(dǎo)出了反應(yīng)堆在碘坑內(nèi)啟堆時碘和氙濃度及氙毒引入的反應(yīng)性隨時間變化的解析表達式；再結(jié)合控制棒價值曲線、溫度效應(yīng)和燃耗效應(yīng)估算碘坑后沿啟堆的臨界棒位。

1 物理模型

HFETR內(nèi)135I和135Xe濃度隨時間變化的方程為[6]：

(1)

式中：NI(t)，NXe(t)——核反應(yīng)堆中135I，135Xe的平均濃度，個/cm3；

γI，γXe——135I，135X的直接裂變產(chǎn)額；

λI，λXe——135I，135X的衰變常數(shù)，s-1；

Σf——燃料的宏觀裂變截面，cm-1；

φ——為堆內(nèi)平均熱中子平通量密度cm-2s-1；

設(shè)反應(yīng)堆在穩(wěn)定功率運行時間超過48 h(即135I和135Xe濃度均接近或達到它們的平衡濃度后)，突然停閉反應(yīng)堆，停堆后堆內(nèi)135I和135Xe的濃度應(yīng)為：

(2)

式中：NI()，NXe()——停堆前堆內(nèi)的135I和135Xe平衡濃度；

根據(jù)反應(yīng)性定義，可以導(dǎo)出毒物所引起反應(yīng)性變化為：

(3)

有效裂變中子數(shù)與燃料宏觀吸收截面的關(guān)系：

(4)

式中：η為有效裂變中子數(shù)，Σf為燃料宏觀裂變截面，ν為每次裂變產(chǎn)生中子數(shù)；

根據(jù)式(2)、式(3)和式(4)可以推出停堆后氙毒反應(yīng)性隨時間變化的計算式：

(5)

式中:φ0為反應(yīng)堆停堆前穩(wěn)定運行時的熱中子通量密度；其中式(5)中的相關(guān)參數(shù)見表1。

表1 相關(guān)參數(shù)Table 1 Related parameters

2 臨界棒位估算過程

為了便于估算，不失一般性選取HFETR某一爐段參數(shù)進行估算，查找運行日志得到這一爐段的首次啟堆的臨界棒位、控制棒價值以及溫度效應(yīng)參數(shù)。

HFETR某爐段臨界棒位：

1，2AB到頂；1，2ZB 500 mm；9，10，11，12，13，14SB到頂；3，6SB 214 mm； 4, 7SB 525 mm；1，2ZB到底；5, 8SB到底。

HFETR某爐段全部控制棒最大價值見表2，另外1、2ZB，1、2SB，3、6SB，4、7SB，5、8SB控制棒相對價值曲線見圖1。

表2 各控制棒價值Table 2 Value of each control rod

圖1 控制棒價值曲線Fig.1 Control rod value curve

HFETR某爐段的溫度效應(yīng)：

由溫度引起的反應(yīng)性溫度變化-2.32×10-4(ΔK/K)/℃

HFETR某爐段的燃耗效應(yīng)：

燃耗引起的反應(yīng)性燃耗系數(shù)為-2×10-4(ΔK/K)/MWd[7]。

臨界棒位估算方法

當(dāng)反應(yīng)堆在碘坑后沿啟動時，反應(yīng)堆處于消毒過程，由于中子通量密度突然增加，氙(135Xe)吸收中子后被大量消耗，反應(yīng)堆內(nèi)會迅速引入正反應(yīng)性，因此在啟動反應(yīng)堆之前給操縱員操提供一個臨界棒位是至關(guān)重要的。

估算碘坑后沿啟動反應(yīng)堆的臨界棒位有兩種方法，第一種是參考第一次啟動反應(yīng)時的臨界棒位；第二種方法是參考停堆時的棒位。

根據(jù)反應(yīng)堆第一次臨界時的臨界棒位可以計算出反應(yīng)堆的初始剩余反應(yīng)性，依據(jù)前面給出的臨界棒位，再結(jié)合控制棒價值曲線計算得到反應(yīng)堆的初始剩余反應(yīng)性為P=14.55βeff。

選取反應(yīng)堆在滿功率運行48 h后由于某些情況導(dǎo)致意外停堆，然后在停堆48h后準(zhǔn)備啟動反應(yīng)堆。使用第一種方法計算臨界棒位，如果參考第一次開堆的臨界棒位，只需考慮燃耗效應(yīng)就可推導(dǎo)出臨界棒位，具體過程分析如圖2。滿功率運行48 h后，燃耗引起的反應(yīng)性變化約為P1=-0.53βeff，氙毒引入的反應(yīng)性變化P2=-4.039βeff，則此時的剩余反應(yīng)性Pi=P+P1+P2=9.980βeff。

圖2 碘坑后沿開堆臨界棒位估算分析過程Fig.2 Estimation and analysis process of critical rod position in the back edge of iodine pit

使用第二種方法計算臨界棒位，參考停堆時的棒位，并需要考慮溫度效應(yīng)和燃耗效應(yīng)推導(dǎo)臨界棒位，具體過程分析如圖2。滿功率運行48 h后，此時各控制棒棒位為：1ZB：364， 2ZB：500，1SB、2SB：375，5SB、8SB到底，其余棒到頂。經(jīng)計算得到反應(yīng)堆此時的剩余反應(yīng)性P3=8.421βeff，氙平衡引入的反應(yīng)性為P4=-4.699βeff。在停堆48 h后準(zhǔn)備啟動反應(yīng)堆后，溫度效應(yīng)引入的反應(yīng)性P5=-1.9βeff，則此時的剩余反應(yīng)性P9=P3-P4+P2-P5=9.982βeff。

最后根據(jù)控制棒價值曲線得到臨界棒位：1，2ZB 500 mm；3，6SB 500 mm；4，7SB 405；1，2SB到底；其余棒到頂。

3 估算程序

基于MATLAB軟件中的GUI工具，結(jié)合第三節(jié)碘坑臨界棒位估算方法開發(fā)程序。圖3和圖4的程序界面為為碘坑后沿啟動反應(yīng)堆臨界棒位的估算程序。

圖3 碘坑后沿啟動反應(yīng)堆臨界棒位的估算程序(一)Fig.3 Estimation program of critical rod position for starting reactor at back edge of iodine pit(1)

圖4 碘坑后沿啟動反應(yīng)堆臨界棒位的估算程序(二)Fig.4 Estimation program of critical rod position for starting reactor at back edge of iodine pit(2)

查閱HFETR運行歷史后，選取存在碘坑后沿啟動反應(yīng)堆的三爐作為研究對象。這里將碘坑后沿啟動反應(yīng)堆的臨界棒位與估算程序得到臨界棒位進行比較,結(jié)果如表3所示。

通過表3可以看出，通過開發(fā)的估算程序得到的臨界棒位與實際棒位的剩余反應(yīng)性差值分別為0.073βeff、0.11βeff、-0.13βeff、0.09βeff。這四組偏差都在0.2βeff之內(nèi)，可以達到用于估算臨界棒位的目的。因此此估算程序?qū)FETR碘坑后沿啟動反應(yīng)堆有一定參考價值。

表3 臨界棒位對比Table 3 Comparison of critical rod positions

4 結(jié)論

本文主要研究了HFETR在碘坑內(nèi)啟動時臨界棒位的估算問題，首先推導(dǎo)出了氙毒引入的反應(yīng)性隨時間變化的解析表達式;接著結(jié)合相關(guān)參數(shù)給出碘坑后沿啟動反應(yīng)堆臨界棒位的估算方法；最后使用MATLAB軟件的GUI工具開發(fā)用于估算碘坑后沿啟動反應(yīng)堆臨界棒位的程序。選取HFETR存在碘坑后沿啟動反應(yīng)堆的三爐作為研究對象，經(jīng)計算得到該程序后備反應(yīng)性的估算偏差在0.2βeff之內(nèi)，可以達到用于估算臨界棒位的目的。因此該估算程序?qū)FETR碘坑后沿啟動反應(yīng)堆有一定參考價值。