李勛昭，吳宏春，鄭友琦，周生誠(chéng)，何明濤

（西安交通大學(xué) 核科學(xué)與技術(shù)學(xué)院，陜西 西安 710049）

近年來，加速器驅(qū)動(dòng)次臨界系統(tǒng)（ADS）作為可使核能大規(guī)模、可持續(xù)發(fā)展的創(chuàng)新技術(shù)路線在國(guó)際上形成了一個(gè)研究熱點(diǎn)，研究表明ADS是理想的長(zhǎng)壽命放射性廢物焚燒爐［1］，是嬗變次錒系核素（MA）的有效技術(shù)途徑之一。ADS主要是由質(zhì)子加速器、散裂靶及次臨界堆芯組成，質(zhì)子加速器產(chǎn)生的高能質(zhì)子轟擊重核散裂靶產(chǎn)生高能的散裂中子，從而驅(qū)動(dòng)次臨界堆芯發(fā)生核反應(yīng)，以達(dá)到嬗變核廢料和增殖核燃料的目的。

目前，對(duì)于ADS反應(yīng)堆核設(shè)計(jì)采用的核數(shù)據(jù)庫主要是基于現(xiàn)有的反應(yīng)堆核設(shè)計(jì)的核數(shù)據(jù)庫或改進(jìn)后的核數(shù)據(jù)庫。大多核數(shù)據(jù)庫的能量上限為20 MeV，而ADS 中，散裂中子源的中子能量分布范圍廣、跨度大，可從10-5eV 到GeV?，F(xiàn)有的專用核數(shù)據(jù)庫忽略20 MeV 以上的中子數(shù)據(jù)對(duì)ADS的影響，這必將給ADS的核設(shè)計(jì)帶來一定的偏差。

為分析高能核數(shù)據(jù)對(duì)ADS性能參數(shù)的影響，本文基于JENDL-HE-2007高能中子評(píng)價(jià)庫［2］，使用NJOY程序［3］加工適用于蒙特卡羅程序計(jì)算的ACE格式的ADS專用核數(shù)據(jù)庫NECL-HE／MC，并采用一系列的基準(zhǔn)題實(shí)驗(yàn)例題對(duì)該核數(shù)據(jù)庫進(jìn)行校核以檢驗(yàn)數(shù)據(jù)庫的合理性。最后基于NECL-HE／MC核數(shù)據(jù)庫，分析不同能量段的散裂中子源對(duì)ADS外中子源效率的影響。

1 NECL-HE／MC核數(shù)據(jù)庫的設(shè)計(jì)與制作

1.1 NECL-HE／MC核數(shù)據(jù)庫的NJOY建模

NECL-HE／MC核數(shù)據(jù)庫主要是針對(duì)ADS的蒙特卡羅程序模擬計(jì)算使用。該核數(shù)據(jù)庫的格式為ACE格式，制作流程如圖1所示。評(píng)價(jià)庫來源日本的JENDL-HE-2007。JENDL-HE-2007包含能量高達(dá)3GeV 的106種核素的中子核反應(yīng)數(shù)據(jù)。由于NJOY 程序版本的問題，其中235U 和14N 兩種核素?zé)o法成功加工制作。本文共加工104種核素的高能核數(shù)據(jù)庫。

圖1 NECL-HE／MC核數(shù)據(jù)庫制作流程Fig.1 Processing sequence of NECL-HE／MC nuclear data library

1.2 NECL-HE／MC核數(shù)據(jù)庫制作參數(shù)選擇

本文加工的核素的溫度均取300K。制作過程中主要參數(shù)選擇列于表1。

表1 NECL-HE／MC核數(shù)據(jù)庫制作參數(shù)Table 1 Parameter of NECL-HE／MC nuclear data library

2 NECL-HE／MC核數(shù)據(jù)庫的數(shù)值驗(yàn)證

2.1 散裂靶模型的計(jì)算

使用MCNPX［4］模擬計(jì)算質(zhì)子束和208Pb靶的相互作用。計(jì)算模型參考MCNPX 說明書中附錄E 的例題1，靶直徑為10cm，長(zhǎng)為30cm，質(zhì)子束直徑為7 cm。分別采用LA150N 數(shù)據(jù) 庫、IAEA-ADS截 面 庫［5］網(wǎng) 站 上1GeV 的208Pb數(shù)據(jù)庫以及NECL-HE／MC 數(shù)據(jù)庫進(jìn)行MCNPX 模擬計(jì)算，質(zhì)子束的能量分別取1、1.5、2、2.5、3GeV。散裂中子產(chǎn)額如圖2所示。3種核數(shù)據(jù)庫吻合良好。1GeV 質(zhì)子束轟擊208Pb靶產(chǎn)生的散裂中子能譜如圖3所示。散裂中子能譜分布趨勢(shì)一致。在低能量段，由于評(píng)價(jià)庫來源不同，散裂中子能譜偏差較大，在高能段，散裂中子能譜曲線幾乎重合。

圖2 采用不同核數(shù)據(jù)庫計(jì)算的不同能量質(zhì)子束轟擊208Pb靶的散裂中子產(chǎn)額Fig.2 Spallation neutron yield of 208Pb target with different energy incident protons and nuclear data libraries

圖3 1GeV 質(zhì)子束轟擊208Pb靶產(chǎn)生的散裂中子能譜Fig.3 Spallation neutron spectrum of 208Pb target with energy of 1GeV incident proton

2.2 臨界基準(zhǔn)實(shí)驗(yàn)裝置的驗(yàn)證

利用國(guó)際臨界安全基準(zhǔn)評(píng)價(jià)實(shí)驗(yàn)（ISBEP2006［6］）中 的 高 富 集235U 裝 置、中 富集235U 裝 置、低 富 集235U 裝 置、Pu 裝 置、U-Pu裝 置和233U 裝 置6 類 臨 界 裝 置 對(duì)NECL-HE／MC 數(shù)據(jù)庫進(jìn)行了數(shù)值驗(yàn)證?；贜ECLHE／MC數(shù)據(jù)庫和西安交通大學(xué)NECP實(shí)驗(yàn)室自主開發(fā)的基于ENDF／B-Ⅶ.0［7］的NECL-1.0／MC［8］數(shù)據(jù)庫，使用MCNP5［9］統(tǒng)計(jì)keff，以及在10×10-5～0.625eV、0.625eV～100keV和100keV 以上3個(gè)能量區(qū)間的中子注量率、并將統(tǒng)計(jì)結(jié)果與臨界安全基準(zhǔn)實(shí)驗(yàn)裝置實(shí)驗(yàn)值進(jìn)行比較。各實(shí)驗(yàn)裝置的描述列于表2。具體的幾何及材料描述見國(guó)際臨界安全基準(zhǔn)評(píng)價(jià)實(shí)驗(yàn)手冊(cè)ISBEP2006。keff計(jì)算結(jié)果如圖4所示。NECL-HE／MC和NECL-1.0／MC的結(jié)果基本一致，NECL-HE／MC 大多的計(jì)算結(jié)果與實(shí)驗(yàn)值吻合良好。

臨界基準(zhǔn)題在10×10-5～0.625 eV、0.625eV～100keV 和100keV 以上3個(gè)能量區(qū)間的中子注量率列于表3，NECL-HE／MC和NECL-1.0／MC 的 結(jié) 果 基 本 一 致，NECLHE／MC的計(jì)算結(jié)果更傾向于實(shí)驗(yàn)值。NECLHE／MC與實(shí)驗(yàn)的keff的結(jié)果基本吻合。

3 NECL-HE／MC核數(shù)據(jù)庫的應(yīng)用

3.1 外中子源效率的定義和計(jì)算模型的建立

外中子源效率是評(píng)價(jià)ADS的重要參數(shù)之一。本文基于NECL-HE／MC核數(shù)據(jù)庫，建立簡(jiǎn)化的ADS計(jì)算模型，模擬散裂中子源對(duì)ADS外中子源效率的影響。外中子源效率φ＊為：

將散裂中子源的能量分為4段，分別為0～5、5～20、20～150、150～1 500 MeV。對(duì)于每個(gè)能量段Si內(nèi)的外中子源效率φ＊i 可表達(dá)為：

通過簡(jiǎn)單的推導(dǎo)，可得到總外中子源效率與各能量段外中子源效率的關(guān)系，表達(dá)式如下：

式中，〈〉表示積分算符。

表2 臨界安全基準(zhǔn)實(shí)驗(yàn)裝置描述Table 2 Description of criticality safety benchmark experiment assembly

圖4 臨界安全基準(zhǔn)題keff比較Fig.4 Comparison of criticality safety benchmark keff

質(zhì)子束能量為1.5GeV，燃料為彌散體燃料（TRU-10Zr）-Zr＊。超鈾核素（TRU）中包含Pu向量和241Am。分別設(shè)置keff為0.94、0.96及0.98。

3.2 計(jì)算結(jié)果及分析

不同次臨界度下各能量段散裂中子源對(duì)外中子源效率的影響列于表4。計(jì)算結(jié)果表明，2.5%的高能散裂中子源（高于20 MeV）對(duì)外中子源效率的貢獻(xiàn)在20%左右，同時(shí)在一定的次臨界范圍內(nèi)，這種貢獻(xiàn)值基本保持穩(wěn)定。散裂中子源的能譜分布對(duì)外中子源效率的影響至關(guān)重要，由于本文所設(shè)計(jì)的模型中，散裂中子源主要分布在0～5MeV的能量區(qū)間內(nèi)，但很明顯的是低能段由于散裂中子源能譜相對(duì)高能區(qū)較軟，故相應(yīng)的外中子源效率比較小，隨著散裂中子源能量的增加，外中子源效率急劇增加，甚至高達(dá)20左右，遠(yuǎn)大于20MeV 以下中子能量段所引起的外中子源效率。散裂中子源的能譜對(duì)外中子源效率的影響很大，同時(shí)，如果忽略高能散裂中子源，將會(huì)給ADS設(shè)計(jì)帶來較大的偏差。

表3 臨界安全基準(zhǔn)題中子注量率比較Table 3 Comparison of criticality safety benchmark neutron fluence rate

表4 不同次臨界度下各能量段散裂中子源對(duì)外中子源效率的影響Table 4 Effect of spallation neutron source for different energy segments on external neutron source efficiency under different sub-criticalities

4 討論及結(jié)論

本工作針對(duì)ADS具有中子能量范圍廣、跨度大的特點(diǎn)，基于JENDL-HE-2007中子評(píng)價(jià)庫，使用NJOY 程序成功地加工NECL-HE／MC核數(shù)據(jù)庫。經(jīng)過簡(jiǎn)單的高能散裂靶和國(guó)際臨界安全基準(zhǔn)題實(shí)驗(yàn)裝置的驗(yàn)證，驗(yàn)證了該核數(shù)據(jù)庫的可靠性與合理性。針對(duì)散裂中子源對(duì)ADS外中子源效率的影響，設(shè)計(jì)了一簡(jiǎn)化的ADS模型，分析了不同次臨界度下，0～5、5～20、20～150、150～1 500 MeV 4個(gè)能量段的散裂中子源對(duì)外中子源效率的影響，研究表明，2.5%的高能散裂中子源（高于20 MeV）對(duì)總外中子源效率的貢獻(xiàn)在20%左右，同時(shí)在一定的次臨界范圍內(nèi)，這種貢獻(xiàn)值基本保持穩(wěn)定。高能散裂中子源不容忽略。

由于高能數(shù)據(jù)基準(zhǔn)題（尤其是達(dá)到GeV 以上數(shù)據(jù)）的缺乏，NECL-HE／MC 核數(shù)據(jù)庫的驗(yàn)證還存在不足，有待進(jìn)一步的驗(yàn)證。高能散裂中子源對(duì)ADS的外中子源影響較大。對(duì)ADS設(shè)計(jì)的其他參數(shù)（如結(jié)構(gòu)材料的原子離位損傷）的影響有待進(jìn)一步分析研究。

［1］ 趙志祥，夏海鷗.加速器驅(qū)動(dòng)次臨界系統(tǒng)（ADS）與核能可持續(xù)發(fā)展［J］.中國(guó)核電，2009，2（3）：202-211.ZHAO Zhixiang，XIA Haihong.Study on ADS and the sustainable development of nuclear energy［J］.China Nuclear Power，2009，2（3）：202-211（in Chinese）.

［2］ JENDL-HE-2007［DB／OL］.2014-10-25.http：∥www-nds.iaea.org／public／download-endf／JENDLHE-2007／.

［3］ MacFARLANE R，MUIR D.The NJOY nuclear data processing system，version 91，LA-12740-M［R］.USA：LANL，1994.

［4］ PELOWITZ D B.MCNPX［TM］user's manual，version 2.5.0，LA-CP-05-0369［R］.USA：LANL，2005.

［5］ ADS-HE nuclear data library［DB／OL］.2014-10-25.http：∥www-nds.iaea.org／ads／adshe.html.

［6］ YOLANDA R，ALI N，ENRICO S，et al.International criticality safety benchmark evaluation working group participants［R］.France：Organization for Economic Cooperation and Development Nuclear Energy Agency，2006.

［7］ ENDF／B-Ⅶ.0［DB／OL］.2014-10-25.http：∥wwwnds.iaea.org／public／download-endf／ENDF-B-Ⅶ.0／.

［8］ 劉志彥.中子學(xué)計(jì)算中的連續(xù)能量核數(shù)據(jù)庫和多群核數(shù)據(jù)庫的開發(fā)［D］.西安：西安交通大學(xué)，2011.

［9］ X-5 Monte Carlo Team. MCNP：A general Monte Carlo N-particle transport code，Version 5，LA-UR-03-1987［R］.USA：LANL，2003.