李原野，王明煌，廉 超，蔣潔瓊

(1.中國科學(xué)技術(shù)大學(xué)，安徽 合肥 230027；2.中國科學(xué)院核能安全技術(shù)研究所，安徽 合肥 230031)

加速器驅(qū)動次臨界釷焚燒堆初步中子學(xué)設(shè)計與分析

李原野1，2，王明煌2，廉 超2，蔣潔瓊2

(1.中國科學(xué)技術(shù)大學(xué)，安徽 合肥 230027；2.中國科學(xué)院核能安全技術(shù)研究所，安徽 合肥 230031)

本文針對加速器驅(qū)動次臨界系統(tǒng)進行釷資源利用的優(yōu)勢，提出了鉛基加速器驅(qū)動次臨界釷焚燒堆(CLEAR-Th)概念。該概念采用釷钚混合氧化物燃料，冷卻劑采用液態(tài)鉛鉍，T91鋼作為包殼材料和結(jié)構(gòu)材料，初始有效增殖因數(shù)keff設(shè)計為0.98，1.5 GeV能量的質(zhì)子流強工作在10mA以內(nèi)，使用大型集成中子學(xué)計算與分析系統(tǒng)VisualBUS和混合評價核數(shù)據(jù)庫HENDL進行計算分析。初步結(jié)果表明：CLEAR-Th設(shè)計實現(xiàn)了穩(wěn)定焚燒釷產(chǎn)能的設(shè)計目標(biāo)，有望實現(xiàn)穩(wěn)定的閉式釷鈾燃料循環(huán)，并具有長壽命次錒系物質(zhì)生成量低、固有安全性良好等特點。

釷；加速器驅(qū)動次臨界系統(tǒng)；中子學(xué)

核能是安全清潔的新能源。核能發(fā)展離不開對核燃料循環(huán)體系的研究和改進。當(dāng)前核燃料循環(huán)體系的主要形式是鈾燃料熱堆一次通過式循環(huán)，這種體系使得核能發(fā)展不得不面對兩大主要問題：鈾資源問題和環(huán)境安全問題[1]。釷燃料的利用可以為核能可持續(xù)發(fā)展開辟一條新途徑。

相比常規(guī)鈾資源，釷燃料具有以下主要特點：① 釷資源儲量豐富；②232Th是很好的可轉(zhuǎn)換材料；③ 更少的超鈾物質(zhì)生成，核廢料處置難度更低；④ 核不擴散特征。故釷燃料擴展了核能工業(yè)可利用資源，也在核不擴散、環(huán)境安全方面相比鈾燃料循環(huán)存在一定優(yōu)勢。加速器驅(qū)動次臨界系統(tǒng)(ADS)由于依靠外源中子驅(qū)動、運行在次臨界狀態(tài)，具備中子能譜硬、中子余額多等特點，適合裝載釷燃料運行。

本文工作針對釷燃料與ADS結(jié)合的優(yōu)勢，給出一個現(xiàn)實可行性較強的鉛基加速器驅(qū)動次臨界釷焚燒堆CLEAR-Th中子學(xué)設(shè)計方案。該概念采用釷钚混合氧化物燃料，冷卻劑采用液態(tài)鉛鉍[2-7]，T91鋼作為包殼和結(jié)構(gòu)材料，初始有效增殖因數(shù)keff為0.98，質(zhì)子能量1.5GeV，流強在10mA以內(nèi)。中子學(xué)分析結(jié)果表明，該設(shè)計實現(xiàn)了深燃耗穩(wěn)定焚燒釷產(chǎn)能的設(shè)計目標(biāo)，具備形成閉式釷燃料循環(huán)的潛力，并有長壽命次錒系物質(zhì)生成量低、固有安全性良好等特點。

1 設(shè)計目標(biāo)、程序與數(shù)據(jù)庫

CLEAR-Th的首要設(shè)計目標(biāo)為在堆運行期間內(nèi)深燃耗穩(wěn)定焚燒釷產(chǎn)能。表1列出了設(shè)計原則與目標(biāo)。本文計算使用了FDS團隊自主研發(fā)的大型集成中子學(xué)計算分析系統(tǒng)VisualBUS4.2及混合評價核數(shù)據(jù)庫HENDL3.0[8-11]。輸運-燃耗計算所用程序和數(shù)據(jù)庫，已經(jīng)過IAEA-ADS、OECD-ADS及BN600等國際基準(zhǔn)例題[12]的校驗測試并符合較好。

表1 設(shè)計原則與目標(biāo)

2 中子學(xué)設(shè)計與分析

2.1 中子學(xué)設(shè)計與優(yōu)化

本文基于在調(diào)研分析國內(nèi)外ADS示范堆及釷燃料反應(yīng)堆設(shè)計基礎(chǔ)上得到的初步方案，經(jīng)過優(yōu)化得到了CLEAR-Th中子學(xué)設(shè)計方案。為簡化計算、突出規(guī)律，本文的堆芯中子學(xué)設(shè)計計算所使用的模型為均勻材料的R-Z模型。

CLEAR-Th燃料初步設(shè)計為釷钚混合氧化物陶瓷燃料。氧化物燃料深燃耗服役性能好，且主體成分ThO2是目前研究最深入、反應(yīng)堆運行經(jīng)驗最為豐富的釷燃料。堆芯啟動需要在釷中加入易裂變材料；钚可以從動力堆乏燃料中獲得，其中子學(xué)性能良好，也為解決乏燃料钚積存問題[12]提供了一種高效途徑。所選用的乏燃料钚同位素含量參照標(biāo)準(zhǔn)壓水堆钚廢料(表2[13])。包殼和結(jié)構(gòu)材料采用T91鋼，機械性能、物理性能及輻照性能均較佳，成本較低，已在工業(yè)中應(yīng)用。冷卻劑采用中子學(xué)和熱工水力性能優(yōu)良的液態(tài)鉛鉍。反射材料選用液態(tài)鉛鉍與T91鋼；屏蔽材料選用B4C。堆芯材料成分以及結(jié)構(gòu)的設(shè)計綜合平衡了中子學(xué)、熱工水力學(xué)、材料結(jié)構(gòu)等多方面的要求[12]，堆體由內(nèi)而外設(shè)計分為：散裂靶、燃料區(qū)、反射層、屏蔽層。燃料區(qū)高度與半徑比H/R初始設(shè)定為1，反射層厚以90cm作為優(yōu)化基礎(chǔ)。散裂靶位于堆芯中心，內(nèi)半徑為10cm、外半徑為24.25cm，1.5GeV高能質(zhì)子與液態(tài)鉛鉍反應(yīng)產(chǎn)生散裂中子。

表2 3000MW壓水堆每年卸料冷卻10年長壽命廢料量參考

2.1.1 易裂變物質(zhì)百分含量優(yōu)化

釷钚混合氧化物燃料中易裂變物質(zhì)(239Pu、241Pu)在全部重核中的百分含量，初步選取了10%、11%、12%、13%、14%、15%六組方案。考慮到設(shè)計目標(biāo)與約束條件，為獲得最高的能量增益設(shè)定循環(huán)初keff≈0.98。基于設(shè)計初步方案首先計算六組方案的在循環(huán)初1000MWt功率下所需的加速器質(zhì)子流強。從表3計算結(jié)果可以看出10%與11%方案的質(zhì)子流強值超過了設(shè)計限值，難以實現(xiàn)設(shè)計目標(biāo)。對12%、13%、14%、15%四個方案進行計算分析比較(見圖1)。14%、15%兩個方案，由于易裂變物質(zhì)含量太高使得初始轉(zhuǎn)換比太小，堆芯中易裂變物質(zhì)新產(chǎn)生量難以彌補消耗量導(dǎo)致keff減小過快，在達到100GWd/t燃耗深度目標(biāo)時已經(jīng)小于循環(huán)初值。12%、13%兩個方案keff隨燃耗變化相對穩(wěn)定，加速器質(zhì)子流強等方面都滿足設(shè)計要求。但計算中兩種方案的keff都超過了0.98的限值，因此在進一步的設(shè)計中需要引入反應(yīng)性控制措施。比較兩者可知：① 12%方案需要付出更多反應(yīng)性控制代價才能實現(xiàn)次臨界安全運行(keff≤0.98)；② 對于超過keff=0.98的“多余”反應(yīng)性，余值越高則越增加了發(fā)生反應(yīng)性失控造成堆芯超臨界事故的潛在風(fēng)險。故從經(jīng)濟性和安全性的角度13%百分含量的方案更符合加速器驅(qū)動次臨界釷焚燒堆設(shè)計的需求。綜上所述，選擇易裂變物質(zhì)百分含量為13%的方案作為CLEAR-Th的核燃料設(shè)計最終確定方案。

表3 六組易裂變物質(zhì)百分含量方案壽期初的質(zhì)子流強

圖1 四組易裂變物質(zhì)百分含量方案的keff隨燃耗變化Fig.1 Variation of keff with Burnup of different fissile percentage composition

2.1.2 燃料區(qū)高度與半徑比優(yōu)化

優(yōu)化堆芯的H/R值可以降低中子泄漏率，減少燃料裝量，獲得較優(yōu)的中子經(jīng)濟性。在上一輪優(yōu)化基礎(chǔ)上，假定堆芯能量輸出為固定值，調(diào)整燃料區(qū)域不同的H/R，得到單位裝料量最大的能量沉積，最大限度提高裝料量的效率。從表4給出的優(yōu)化計算結(jié)果，在維持包層功率基本一致的條件下，單位裝料量的能量沉積值隨著H/R值的增加而增加；在達到一定比例后，單位裝料量的能量沉積值開始下降。最優(yōu)結(jié)果為堆芯燃料區(qū)尺寸取H=230cm，R=124cm。

表4 不同H/R條件下中子學(xué)性能

2.1.3 徑向反射層厚度優(yōu)化

反射層可減少堆芯中子泄漏提高中子經(jīng)濟性；兼顧設(shè)計成本反射層不能過厚。在上文兩輪優(yōu)化基礎(chǔ)上計算了徑向反射層厚度對keff

的影響(見圖2)。由圖2可知，當(dāng)徑向反射層比較薄時，堆芯keff值也相應(yīng)比較??；隨著反射層增厚，keff值也相應(yīng)地增加，堆芯中子的泄漏幾率減小中子經(jīng)濟性提高；當(dāng)反射層厚度達到60cm以后，keff值趨于穩(wěn)定在0.981附近，徑向反射層的加厚也未引起keff顯著的增加。綜合上述，反射層厚度取60cm較為合理。

圖2 徑向反射層厚度對keff的影響Fig.2 The effect of the thickness of radial reflector on keff

基于以上優(yōu)化工作的結(jié)論，本文最終確定了CLEAR-Th的設(shè)計方案。詳細參數(shù)見表5以及圖3。

表5 CLEAR-Th材料成分和尺寸

圖3 CLEAR-Th反應(yīng)堆堆芯R-Z模型圖Fig.3 Radial-axial CLEAR-Th reactor core configuration

2.2 燃耗分析

針對CLEAR-Th方案，固定熱功率1000MWt，以堆芯燃料平均燃耗目標(biāo)100GWd/t進行燃耗計算。表6給出了壽期初/末主要參數(shù)的計算結(jié)果，可知：① 雖然堆芯的平均燃耗達到100GWd/t但keff并沒有減??；② 兩項反應(yīng)性系數(shù)負(fù)反饋加強，固有安全性得到強化；③ 轉(zhuǎn)換比從1.12減小到堆芯壽期末的0.93，這是由在燃料循環(huán)內(nèi)232Th不斷消耗減少以及新生成的233U具有較大的裂變截面造成的；④ 在循環(huán)內(nèi)，釷被焚燒消耗掉13.2%；⑤ 核燃料中主要的易裂變核素已經(jīng)由循環(huán)初的239Pu轉(zhuǎn)變?yōu)樾略鲋钞a(chǎn)生的233U。

表6 堆芯物理性能表

圖4給出keff在100GWd/t燃耗深度內(nèi)變化情況。在設(shè)定100GWd/t燃耗深度目標(biāo)內(nèi)，keff經(jīng)歷了逐漸增大而后逐漸減小的變化。循環(huán)初232Th、240Pu含量較高轉(zhuǎn)換比相對高，易裂變物質(zhì)生成速率高于消耗速率，keff逐漸增大。隨著燃耗加深堆芯轉(zhuǎn)換比持續(xù)變小，新生成的易裂變核素已經(jīng)無法彌補反應(yīng)性的損失，于是在后期keff逐漸減小。

圖4 keff隨燃耗深度的變化Fig.4 Evolution of core raw keff over Burnup

表7選取壽期初和壽期末兩個時刻，對比考察以下5種重要組分的含量：232Th，鈾，钚，232U，長壽命次錒系物質(zhì)(LLMA，包括237Np241Am243Am244Cm)。① 釷利用率至壽期末達到了13.2%，約2.47 t232Th被焚燒；② 到壽期末鈾凈增殖產(chǎn)生了1.383t，有效補充了核燃料中钚的消耗；③ 钚在循環(huán)中被消耗掉40.21%；④ 強放射性的232U至壽期末生成約1.64kg，在鈾中濃度為0.11%，在這個濃度值下分離提取易裂變物質(zhì)需要遠距離加工[14]，有利于防止核擴散；⑤ LLMA中237Np僅生成0.38kg，相比3000MW標(biāo)準(zhǔn)壓水堆年卸料14.5kg[12]生成量而言很少，遠遠低于另外三種LLMA的生成量，這得益于釷的特性。

表7 壽期初/末核燃料5種重要組分質(zhì)量

2.3 反應(yīng)性系數(shù)初步計算分析

反應(yīng)性系數(shù)計算對于反應(yīng)堆設(shè)計瞬態(tài)特性分析和安全運行有重要意義。表6中已經(jīng)給出CLEAR-Th中子學(xué)設(shè)計壽期初的兩項反應(yīng)性系數(shù)，其余計算結(jié)果由表8給出。燃料多普勒系數(shù)、冷卻劑溫度系數(shù)呈現(xiàn)出負(fù)反饋效應(yīng)，這有助于在發(fā)生嚴(yán)重事故反應(yīng)堆溫度升高時減緩或限制事故的擴大。負(fù)膨脹系數(shù)的反饋作用有利于維持堆芯形狀穩(wěn)定。顯著的負(fù)空泡反應(yīng)性系數(shù)得益于釷燃料的使用，當(dāng)出現(xiàn)冷卻劑沸騰時，反應(yīng)堆有較強的負(fù)反饋調(diào)節(jié)能力。而鉛鉍冷卻劑的沸點高達1943K，甚至高于結(jié)構(gòu)材料T91鋼的熔點，發(fā)生冷卻劑沸騰的可能性極小。同時加之ADS系統(tǒng)的次臨界深度，CLEAR-Th具有良好的固有安全性。

表8 CLEAR-Th壽期初反應(yīng)性系數(shù)計算結(jié)果

2.4 燃料循環(huán)管理分析

基于CLEAR-Th的設(shè)計，假設(shè)采用簡單的閉式循環(huán)的燃料方式，即固定1000MWt功率及100GWd/t燃耗，每循環(huán)結(jié)束后去除核燃料中裂變產(chǎn)物并補充釷消耗，隨后制成新燃料返回反應(yīng)堆繼續(xù)運行。圖5給出了6輪連續(xù)的燃料循環(huán)的核燃料主要物質(zhì)成分計算結(jié)果，每個循環(huán)時長約6年。在第3循環(huán)后，核燃料中主要易裂變物質(zhì)233U的成分含量已經(jīng)接近平衡；在第5循環(huán)結(jié)束后，核燃料中鈾的總量基本穩(wěn)定，钚的消耗速率也逐漸降低。第5循環(huán)與第6循環(huán)結(jié)束時，钚的嬗變消耗率分別達到了84.1%和87.8%，堆芯初始啟動用的钚被有效嬗變掉。這種趨于穩(wěn)定的發(fā)展態(tài)勢表明，僅采用上述簡單閉式循環(huán)方法，基于CLEAR-Th初步中子學(xué)概念設(shè)計有望實現(xiàn)閉式釷鈾燃料循環(huán)，從而實現(xiàn)對釷燃料的高效利用。

圖5 在6個燃料循環(huán)中，钚的消耗和鈾的產(chǎn)生情況Fig.5 Depletion of plutonium and creation of uranium in 6 cycles

3 總結(jié)

本文針對核能開發(fā)背景下核燃料循環(huán)體系的研究，提出了將ADS與釷鈾燃料循環(huán)優(yōu)勢結(jié)合的鉛基加速器驅(qū)動次臨界釷焚燒堆CLEAR-Th的初步中子學(xué)概念設(shè)計。該系統(tǒng)以釷钚混合氧化物作為啟動燃料，能夠?qū)崿F(xiàn)深燃耗穩(wěn)定運行。基于該系統(tǒng)有望實現(xiàn)閉式釷鈾燃料循環(huán)。該系統(tǒng)具有長壽命次錒系核素生成量少、固有安全性良好等特點。后續(xù)將在熱工和結(jié)構(gòu)設(shè)計上開展方案的進一步深入研究。

致謝

本文得到了國家自然科學(xué)基金(91026004和51406216)的資助，同時非常感謝FDS團隊其他成員對本工作的大力支持。

[1] 顧忠茂. 核能與先進核燃料循環(huán)技術(shù)發(fā)展動向[J]. 現(xiàn)代電力，2006，23(5)：89-94.

[2] Y. Wu，F(xiàn)DS Team. Conceptual Design and Testing Strategy of a Dual Functional Lithium-Lead Test Blanket Module in ITER and EAST[J]. Nuclear Fusion，2007，47(11)：1533-1539.

[3] Y. Wu，Y. Bai，Y. Song，et al. Development Strategy and Design Options of China Lead-based Research Reactor for ADS Transmutation System[J]. Energy Conversion and Management.

[4] Q. Huang，S. Gao，Z. Zhu，et al. Progress in Compatibility Experiments on Lithium-Lead with Candidate Structural Materials for Fusion in China[J]. Fusion Engineering and Design，2009，84：242-246.

[5] Y. Wu. Design Status and Development Strategy of China Liquid Lithium-Lead Blankets and Related Material Technology[J]. Journal of Nuclear Materials，2007，367-370：1410-1415.

[6] 吳宜燦，柏云清，宋勇，等. 中國鉛基研究反應(yīng)堆概念設(shè)計研究[J]. 核科學(xué)與工程，2014(2).

[7] 吳宜燦，黃群英，朱志強，等. 中國系列液態(tài)鋰鉛實驗回路設(shè)計與研發(fā)進展[J]. 核科學(xué)與工程，2009，29(2)：161-169.

[8] Yican Wu，F(xiàn)DS Team. CAD-based interface programs for fusion neutron transport simulation[J]. Fusion Eng. Des，2009，84：1987-1992.

[9] 吳宜燦，李靜驚，李瑩，等. 大型集成多功能中子學(xué)計算與分析系統(tǒng)VisualBUS的研究與發(fā)展[J]. 核科學(xué)與工程，2007，27(4)：365-373.

[10] Y. Wu，Z. Xie，U. Fischer. A Discrete Ordinates Nodal Method for One-Dimensional Neutron Transport Calculation in Curvilinear Geometries[J]. Nuclear Science and Engineering，1999，133(3)：350-357.

[11] 吳宜燦，胡麗琴，龍鵬程，等. 先進核能軟件發(fā)展與核信息學(xué)實踐[M]. 中國科研信息化藍皮書，2013，北京：科學(xué)出版社，ISBN 978-7-03-039323-4，2013.12：232-24.

[12] 陳忠. 加速器驅(qū)動核廢料嬗變次臨界堆中子學(xué)設(shè)計研究[D]. 中國科學(xué)院研究生院，2012.

[13] 鄭善良，吳宜燦，高純靜，等. 聚變驅(qū)動次臨界堆雙冷嬗變包層中子學(xué)設(shè)計與分析[J]. 核科學(xué)與工程，2004，24(2)：164-170.

[14] 張家驊，包伯榮，夏源賢，等. 高通量堆輻照二氧化釷中233U，232U，233Pa及裂變產(chǎn)物的積累與中子積分通量的關(guān)系[J]. 核科學(xué)與工程，1986，6(2)：192-196.

PreliminaryNeutronicsDesignandAnalysisforAcceleratorDrivenSubcriticalThoriumBurningReactor

LIYuan-ye1，2，WANGMing-huang2，LIANChao2，JIANGJie-qiong2

(1.School of Nuclear Science and Technology，University of Science and Technology of China，Hefei of Anhui Prov. 230027，China；2.Institute of Nuclear Energy Safety Technology，Chinese Academy of Sciences，Hefei of Anhui Prov. 230031，China)

In view of combining the advantages of thorium with accelerator driven subcritical system，the concept of Accelerator Driven Subcritical Thorium Burning Reactor(CLEAR-Th)was put forward in this paper. CLEAR-Th is fueled with thorium-plutonium mixes oxide and cooled by LBE. T91 steel is used as cladding and structure material. Thekeffis set as 0.98 at the beginning of fuel cycle(BOC). The accelerator proton(1.5 GeV)beam current is limited to a maximum of 10 mA. Preliminary neutronics design analyses were performed on it by using multi-functional 4D neutronics simulation system VisualBUS and hybrid evaluated nuclear data library HENDL. The preliminary results showed that the primary design goal of stable energy production though burning thorium has been achieved. Based on this design，stable closed Th-U fuel cycle can be achievable. CLEAR-Th is also featured with low production of long-lived minor actinides and good inherent safety.

Thorium；Accelerator driven subcritical system；Neutronics

2017-09-21

國家自然科學(xué)基金項目“加速器驅(qū)動次臨界堆瞬態(tài)安全過程與影響機理研究”(91026004)、“乏燃料焚燒堆新概念次臨界包層長壽命焚燒特性研究”(51406216)

李原野(1990—)，山東日照人，碩士研究生，現(xiàn)從事反應(yīng)堆物理方面研究工作

王明煌：minghuang.wang@fds.org.cn

TL32

A

0258-0918(2017)06-0895-07