牛江++黃禮淵++代啟東

摘 要：該試驗(yàn)基于臨界堆芯的充硼試驗(yàn)，充硼試驗(yàn)指的是向臨界堆芯中逐漸添加硼酸溶液，并通過提升某組控制棒來補(bǔ)償硼酸所引入的負(fù)反應(yīng)性，以使堆芯再次返回臨界狀態(tài)。在對(duì)該堆芯進(jìn)行充硼試驗(yàn)的過程中，針對(duì)多個(gè)臨界堆芯狀態(tài)的測(cè)點(diǎn)進(jìn)行試驗(yàn)，使堆芯硼濃度從零逐漸增加，并通過提升控制棒而維持堆芯在不同硼濃度情況下的臨界狀態(tài)，試驗(yàn)時(shí)利用數(shù)字反應(yīng)性儀判定堆芯臨界狀態(tài)，記錄各臨界狀態(tài)下的硼濃度及臨界棒位，并將試驗(yàn)結(jié)果與MCNP的計(jì)算結(jié)果相比較。

關(guān)鍵詞：臨界堆芯 充硼試驗(yàn) 數(shù)字反應(yīng)性儀 MCNP

中圖分類號(hào)：TL48 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1674-098X（2017）03（c）-0045-04

Critical Experiment Research of Boron-filling Core

Niu Jiang Huang Liyuan Dai Qidong

（Reactor Engineering Research Division， Nuclear Institute of China， Cheng du Sichuan， 610041， China）

Abstract：The experiment is based on boron-filling experiment. Boron-filling experiment means adding boron solution into critical core， and elevating control rod to offset the negative reactivity which brought by boron， then making the core return to critical state. In the process of boron-filling experiment， research was carried out for multiple critical core. With increase of boron concentration， control rod was pulled out to make the core critical in different boron state. During the experiment， digital reactivity meter was used to judge whether the core was critical. Boron concentration and control rods critical position were recorded， then the experiment results were compared to MCNPs results.

Key Words：Critical core； Boron-filling experiment； Digital reactivity meter； MCNP

該試驗(yàn)所使用的數(shù)字反應(yīng)性儀基于逆動(dòng)態(tài)方法跟蹤堆芯在臨界狀態(tài)附近的動(dòng)態(tài)反應(yīng)性，可通過在控制臺(tái)上微調(diào)控制棒，監(jiān)測(cè)反應(yīng)性儀可視化界面的反應(yīng)性曲線變化，最終確定堆芯處于臨界狀態(tài)（反應(yīng)性ρ=0）。數(shù)字反應(yīng)性儀采用γ補(bǔ)償電離室將中子計(jì)數(shù)轉(zhuǎn)換成電流信號(hào)進(jìn)行記錄并處理。該試驗(yàn)的目的是完成充硼試驗(yàn)，在逐漸增加堆芯硼濃度的過程中，可利用數(shù)字反應(yīng)性儀給出不同硼濃度狀態(tài)下堆芯的臨界棒位。

1 方法原理

1.1 逆動(dòng)態(tài)方法測(cè)量反應(yīng)性

該試驗(yàn)主要利用數(shù)字反應(yīng)性儀對(duì)堆芯進(jìn)行反應(yīng)性測(cè)量，而數(shù)字反應(yīng)性儀基于逆動(dòng)態(tài)方法，逆動(dòng)態(tài)方法的原理如下：

由點(diǎn)堆動(dòng)力學(xué)方程（1）：

（1）

其中：N（t）為t時(shí)刻堆內(nèi)的平均中子密度，ρ（t）為t時(shí)刻的反應(yīng)性；β為緩發(fā)中子的總份額，為第i組緩發(fā)中子的份額；為中子代時(shí)間；為第i組衰變常數(shù)；Ci（t）為第i組緩發(fā)中子先驅(qū)核的濃度。

經(jīng)過一系列推導(dǎo)，可得出反應(yīng)堆剛啟動(dòng)或已停止運(yùn)行很長(zhǎng)時(shí)間后的反應(yīng)性求解公式（2）：

（2）

公式（2）中的λi及βi為常數(shù)，時(shí)間t及τ為可控變量，中子密度N（t）為待測(cè)參數(shù)，因此，只要控制好時(shí)間變量，并測(cè)量中子密度N（t）隨時(shí)間的變化，即可得出單位為元（1個(gè)β）的反應(yīng)性。

該試驗(yàn)中的數(shù)字反應(yīng)性儀基于逆動(dòng)態(tài)方法來判定堆芯的臨界狀態(tài)，數(shù)字反應(yīng)性儀的探測(cè)器采集到中子計(jì)數(shù)時(shí)，將產(chǎn)生與中子計(jì)數(shù)成正比的電流信號(hào)，當(dāng)產(chǎn)生的電流信號(hào)不隨時(shí)間發(fā)生變化時(shí)，即可確定堆芯此時(shí)已處于臨界狀態(tài)。

1.2 硼效應(yīng)

該試驗(yàn)向堆芯中注入的硼溶液含10B，10B原子具有很大的熱中子反應(yīng)截面（如圖1所示），由于主要靠熱中子與裂變材料發(fā)生裂變反應(yīng)而產(chǎn)生能量，故向堆芯引入的10B吸收大量熱中子時(shí)，會(huì)降低熱中子與裂變材料的反應(yīng)幾率而引入負(fù)反應(yīng)性，即引入負(fù)反應(yīng)性效應(yīng)。

2 試驗(yàn)實(shí)施

2.1 試驗(yàn)裝置

該試驗(yàn)裝置為零功率物理試驗(yàn)堆芯，其由堆桶、核測(cè)系統(tǒng)、控制保護(hù)系統(tǒng)、棒控系統(tǒng)、劑量監(jiān)測(cè)系統(tǒng)、慢化劑回路系統(tǒng)及操縱臺(tái)等構(gòu)成。燃料組件的芯塊材料UO2，控制棒材料為鉿。如圖2所示，試驗(yàn)堆芯包含九盒燃料組件，組件內(nèi)的控制棒包括S1、S2安全棒和C1、C2手動(dòng)棒，組件由2 cm后的堆芯圍板包圍，圍板與堆桶之間充滿了水慢化劑。

2.2 試驗(yàn)過程

該試驗(yàn)過程的步驟是：（1）堆芯不含硼時(shí)，利用數(shù)字反應(yīng)性儀及堆上儀表監(jiān)測(cè)堆芯狀態(tài)，將S1、S2、C1均提升至堆頂，然后提升C2棒而使反應(yīng)堆達(dá)臨界；（2）向堆芯注入硼酸溶液，注入一定量硼酸溶液后，可引入負(fù)反應(yīng)性效應(yīng)而使堆芯處于次臨界狀態(tài)，此時(shí)可分段提升C2棒，使反應(yīng)堆再次返回臨界；（3）重復(fù)步驟（2），直至將C2棒提升至堆頂附近。

3 試驗(yàn)結(jié)果分析

3.1 堆芯建模

MCNP為美國(guó)Los Alamos國(guó)家實(shí)驗(yàn)室（LANL）開發(fā)的一種三維多粒子輸運(yùn)程序，其可用于模擬中子、光子、電子等輸運(yùn)模型，尤其適用于反應(yīng)堆的臨界計(jì)算模型。MCNP程序可對(duì)燃料組件、堆內(nèi)結(jié)構(gòu)及反射層等進(jìn)行精確建模，MCNP在反應(yīng)堆臨界計(jì)算中的應(yīng)用經(jīng)受了國(guó)內(nèi)外多次試驗(yàn)驗(yàn)證，其給出的有效增殖系數(shù)Keff得到了廣泛認(rèn)可。

在利用MCNP對(duì)該試驗(yàn)堆芯（如圖2所示）進(jìn)行建模時(shí)，利用MCNP的柵元卡及曲面卡依次按照燃料柵元、燃料組件、堆芯、反射層、圍板、堆桶的順序進(jìn)行精確描述，并設(shè)置控制棒的提升高度可調(diào)節(jié)。將試驗(yàn)過程中每次充硼后的試驗(yàn)臨界棒位及硼濃度作為輸入，寫入利用MCNP5所建立的模型之中，進(jìn)行臨界計(jì)算，得出當(dāng)前棒位下的有效增殖系數(shù)。

3.2 數(shù)據(jù)分析

表1給出了根據(jù)各測(cè)點(diǎn)的實(shí)際臨界棒位和根據(jù)該棒位所計(jì)算出的增殖系數(shù)及計(jì)算偏差，可以看出，隨著硼濃度的增加，需要逐步提升C2棒來維持反應(yīng)堆臨界。圖3為各測(cè)點(diǎn)的計(jì)算增殖系數(shù)（Keff）變化圖，可看出整個(gè)充硼過程的增殖系數(shù)誤差控制在2‰（約200 pcm）以內(nèi)，圖4為各測(cè)點(diǎn)計(jì)算增殖系數(shù)（Keff）的計(jì)算偏差變化圖，最大計(jì)算偏差為1‰。由于數(shù)字反應(yīng)性儀基于逆動(dòng)態(tài)方法判定堆芯達(dá)臨界時(shí)，控制臺(tái)上的周期表指向無窮大且功率表處于穩(wěn)定狀態(tài)，此時(shí)，利用試驗(yàn)臨界棒位及硼濃度作為輸入的蒙特卡洛程序計(jì)算結(jié)果也呈臨界狀態(tài)，故堆芯建模是準(zhǔn)確的。

4 結(jié)語(yǔ)

在反應(yīng)堆堆芯充硼臨界試驗(yàn)過程中，隨著硼濃度的增加，堆芯由臨界狀態(tài)轉(zhuǎn)入次臨界狀態(tài)，故硼具有負(fù)反應(yīng)性效應(yīng)，需要通過提升控制棒而使堆芯再次達(dá)臨界。通過對(duì)堆芯進(jìn)行多次充硼，最終將所有控制棒提出堆外，使反應(yīng)堆僅在硼溶液的作用下達(dá)臨界，此時(shí)，可排除控制棒棒間干涉對(duì)反應(yīng)性測(cè)量的影響。利用數(shù)字反應(yīng)性儀基于逆動(dòng)態(tài)跟蹤方法判定反應(yīng)堆的臨界狀態(tài)，其臨界判定偏差很小，方法直接，能夠有效提高試驗(yàn)效率，且對(duì)試驗(yàn)過程中硼微積分價(jià)值、控制棒微積分價(jià)值等試驗(yàn)數(shù)據(jù)測(cè)量精度的提高有很大幫助。通過將試驗(yàn)結(jié)果與計(jì)算結(jié)果相比較，可以看出，數(shù)字反應(yīng)性儀適用于不同硼濃度的臨界堆芯。

參考文獻(xiàn)

[1] 史永謙.核反應(yīng)堆中子學(xué)實(shí)驗(yàn)技術(shù)[M].北京：中國(guó)原子能出版社，2011.

[2] 謝仲生.核反應(yīng)堆物理分析[M].西安：西安交通大學(xué)出版社，2004.

[3] 許淑艷.蒙特卡羅方法在實(shí)驗(yàn)核物理中的應(yīng)用[M].北京：中國(guó)原子能出版社，1996.