吳宜燦，鄒 俊，郝麗娟，王明煌，楊 琪，宋 婧，汪 進，尚雷明，龍鵬程，王 芳，胡麗琴，何 桃，F(xiàn)DS團隊

(中國科學院核能安全技術(shù)研究所，中子輸運理論與輻射安全重點實驗室，安徽合肥230031)

?

混合評價核數(shù)據(jù)庫系統(tǒng)HENDL3.0研發(fā)及其在先進核能系統(tǒng)設(shè)計中應(yīng)用

吳宜燦，鄒 俊，郝麗娟，王明煌，楊 琪，宋 婧，汪 進，尚雷明，龍鵬程，王 芳，胡麗琴，何 桃，F(xiàn)DS團隊

(中國科學院核能安全技術(shù)研究所，中子輸運理論與輻射安全重點實驗室，安徽合肥230031)

FDS團隊結(jié)合聚變堆、先進裂變堆、聚變驅(qū)動次臨界堆、加速器驅(qū)動次臨界堆等核能系統(tǒng)對應(yīng)用核數(shù)據(jù)的需求，設(shè)計開發(fā)出混合評價核數(shù)據(jù)庫HENDL。HENDL包括輸運核數(shù)據(jù)庫、嬗變與活化核數(shù)據(jù)庫、輻照損傷核數(shù)據(jù)庫等，并針對先進核能系統(tǒng)的物理特點，從能量自屏效應(yīng)、熱散射效應(yīng)、溫度多普勒效應(yīng)等方面進行了精確的截面修正。HENDL已通過國際基準模型與實驗的驗證與確認，并在ADS-CLEAR 、FDS-SFB為代表的重大核工程與核能研究中得到了廣泛應(yīng)用。

先進核能系統(tǒng)；多群數(shù)據(jù)庫；基準檢驗

在進行核反應(yīng)堆中子學計算時，需要知道不同能量的粒子與各種物質(zhì)(包括核燃料、冷卻劑、結(jié)構(gòu)材料等)相互作用的核反應(yīng)截面，以及衰變常數(shù)、衰變分支比和裂變產(chǎn)物份額等有關(guān)參數(shù)，這些數(shù)據(jù)統(tǒng)稱為核數(shù)據(jù)。核數(shù)據(jù)是核科學技術(shù)研究和核工程設(shè)計所必需的基本數(shù)據(jù)，也是反應(yīng)堆中子學計算的出發(fā)點和依據(jù)。

在實驗核數(shù)據(jù)方面，國際原子能機構(gòu)核數(shù)據(jù)科(International Atomic Energy Agency/Nuclear Data Section，IAEA/NDS)發(fā)布的實驗核數(shù)據(jù)庫EXFOR[1]匯總了全世界的核數(shù)據(jù)測量的主要實驗室的測量數(shù)據(jù)。在評價核數(shù)據(jù)方面，主要有美國的ENDF[2]、歐洲的JEFF[3]、俄羅斯的BROND[4]、日本的JENDL[5]、中國的CENDL[6]等。在應(yīng)用核數(shù)據(jù)方面，主要包括應(yīng)用于裂變堆的WIMSD[7]多群庫、應(yīng)用于聚變裝置的FENDL[8]數(shù)據(jù)庫、應(yīng)用于加速器驅(qū)動次臨界堆的ADS.lib[9]數(shù)據(jù)庫等。

傳統(tǒng)裂變堆設(shè)計與安全分析所使用的核數(shù)據(jù)庫相對成熟與完善，但針對尚處在研究階段的反應(yīng)堆，如加速器驅(qū)動次臨界堆、聚變驅(qū)動次臨界堆等先進核能系統(tǒng)，尚無很成熟的應(yīng)用核數(shù)據(jù)庫。國際上的ADS核數(shù)據(jù)庫主要有IAEA/NDS發(fā)布的ADS-lib、美國LANL發(fā)布的La150[10]、日本原子能機構(gòu)發(fā)布的JENDL-He[11]，但這些核數(shù)據(jù)存在能量范圍無法覆蓋ADS能譜范圍、核素種類過少、物理效應(yīng)修正不充分等問題。另外，國際上相繼開發(fā)了應(yīng)用于聚變堆、裂變堆應(yīng)用數(shù)據(jù)庫，但缺乏專門應(yīng)用于聚變驅(qū)動次臨界堆的數(shù)據(jù)庫。

中國科學院核能安全技術(shù)研究所·FDS團隊結(jié)合聚變反應(yīng)堆、先進裂變堆、聚變驅(qū)動次臨界堆、加速器驅(qū)動次臨界堆等先進核能系統(tǒng)對應(yīng)用核數(shù)據(jù)的需求，以最新發(fā)布的國際主流評價數(shù)據(jù)源為基礎(chǔ)，依據(jù)國際基準臨界安全實驗、中子屏蔽積分實驗?zāi)M結(jié)果對評價數(shù)據(jù)進行了篩選，結(jié)合對先進核能系統(tǒng)的物理分析，設(shè)計開發(fā)出混合評價核數(shù)據(jù)庫系統(tǒng)HENDL[12-19]。本論文主要從該數(shù)據(jù)庫的體系架構(gòu)與特點、關(guān)鍵技術(shù)研究、核數(shù)據(jù)驗證與應(yīng)用展開闡述。

1 體系架構(gòu)與特色

混合評價核數(shù)據(jù)庫系統(tǒng)HENDL總體目標是為先進核能系統(tǒng)的核設(shè)計及安全分析提供應(yīng)用核數(shù)據(jù)庫支持，目前最新版本為3.0。HENDL主要包括粒子輸運數(shù)據(jù)庫、嬗變/活化數(shù)據(jù)庫、輻照損傷庫等功能子庫。

HENDL同時具備支持多種能群結(jié)構(gòu)，涵蓋多種堆型的能譜，包括多種物理計算功能子庫等特點，服務(wù)于多種核能系統(tǒng)尤其是先進核能系統(tǒng)的核設(shè)計與安全分析。

在能群結(jié)構(gòu)方面，HENDL包括連續(xù)點截面數(shù)據(jù)庫HENDL/MC[13]以及群狀數(shù)據(jù)庫，主要有366群中子與42群光子能群結(jié)構(gòu)的HENDL-ADS/MG[14]數(shù)據(jù)庫、315群中子與42群光子能群結(jié)構(gòu)的HENDL/FG[12]細群庫、175群中子與42群光子能群結(jié)構(gòu)的HENDL/MG[9-10]多群庫、27群中子與21群光子能群結(jié)構(gòu)的HENDL/CG[11]粗群數(shù)據(jù)庫。在堆型能譜支持方面，HENDL的能譜包括標準裂變譜、聚變譜、聚變-裂變混合譜以及加速器驅(qū)動次臨界堆能譜，滿足不同類型反應(yīng)堆物理計算的需求。在多種物理計算功能子庫方面，HENDL以基本的群狀主庫為基礎(chǔ)，通過加工處理得到了各個輸運庫對應(yīng)的其他功能子庫，包括HENDL-A(Activation)嬗變/活化庫、HENDL-RF(Response-Function)輻照損傷庫等，這些子庫涵蓋MG/FG/CG能群結(jié)構(gòu)下的數(shù)據(jù)。

在應(yīng)用領(lǐng)域方面，HENDL/MG主要服務(wù)于裂變堆、聚變堆、聚變裝置的核分析；HENDL/FG主要服務(wù)于聚變驅(qū)動次臨界堆、裂變堆的核分析；HENDL-ADS/MG主要服務(wù)于加速器驅(qū)動次臨界堆的核分析。HENDL/MC可以為反應(yīng)堆輸運計算提供連續(xù)點截面支持。

2 關(guān)鍵技術(shù)

FDS團隊針對聚變驅(qū)動次臨界堆、加速器驅(qū)動次臨界堆等先進核能系統(tǒng)的物理特點，從能群結(jié)構(gòu)、權(quán)重函數(shù)的設(shè)計，以及能量共振自屏、熱中子散射、多普勒等物理效應(yīng)的精確修正展開的設(shè)計研究，并發(fā)展了支持先進核能系統(tǒng)設(shè)計與分析的混合評價核數(shù)據(jù)庫系統(tǒng)HENDL。

2.1 能群結(jié)構(gòu)設(shè)計

在多群數(shù)據(jù)庫的設(shè)計中，影響群截面設(shè)計精度的一個重要因素就是能群結(jié)構(gòu)的合理設(shè)計，理論上能群劃分越精細計算越精確，但是能群過細也將帶來計算耗時等問題。因此根據(jù)核素核反應(yīng)截面隨能量變化的特性進行合理的群截面設(shè)計，可以達到快速精確計算的目的。能群結(jié)構(gòu)的設(shè)計需考慮因素包括：各種閾能反應(yīng)，非彈性散射、某些核素的裂變反應(yīng)、(n，2n)和(n，3n)反應(yīng)等的閾能應(yīng)當作為群的間隔點；截面發(fā)生明顯變化的中子能量即拐點處應(yīng)該作為能群的邊界。

針對先進核能系統(tǒng)20 MeV以上中高能中子與次臨界堆包層材料相互作用的理論模型(光學散射、復合核、直接反應(yīng)及預蒸發(fā)模型)分析基礎(chǔ)上，獲取先進核能系統(tǒng)關(guān)鍵核素的中高能新的反應(yīng)道、高能核反應(yīng)閾以及截面明顯變化點，在此基礎(chǔ)上綜合分析中子能群結(jié)構(gòu)的設(shè)計方法。

2.2 權(quán)重函數(shù)設(shè)計

在能群結(jié)構(gòu)選取一定的情況下，權(quán)重函數(shù)越精確對應(yīng)的群截面的精度也越高。經(jīng)過幾十年的研究，中子物理學家給出了一些適用于不同的堆型近似譜作為權(quán)重譜，但是對于先進核能系統(tǒng)中子權(quán)重函數(shù)并未給出合適的分布。在國內(nèi)外目前的設(shè)計中，對于先進核能系統(tǒng)高能區(qū)的權(quán)重函數(shù)大多采用平滑的譜分布，如EASY自帶的燃耗庫Vitamin-J+以及TRPOLI+均采用平滑權(quán)重譜函數(shù)(權(quán)重譜為常數(shù))。

在HENDL核數(shù)據(jù)庫的權(quán)重函數(shù)設(shè)計中，應(yīng)用典型的先進核能系統(tǒng)設(shè)計方案，對其中子能譜展開研究與分析，重點探討20 MeV以上中子譜分布特點，并采用最小二乘法對中高能中子譜的分布進行研究，從而獲取合適的先進核能系統(tǒng)權(quán)重函數(shù)分布。

2.3 能量共振自屏效應(yīng)修正

由于核素的吸收截面在共振能區(qū)(錒系與次錒系主要集中在4 eV～9.118 keV；中等質(zhì)量核主要集中在0.1 keV～1 MeV)急劇增大，形成強烈的共振峰，這使得反應(yīng)堆內(nèi)共振能區(qū)中子通量密度明顯下降，這個現(xiàn)象稱為能量自屏效應(yīng)[19]。在聚變驅(qū)動次臨界堆、加速器驅(qū)動次臨界堆等先進核能系統(tǒng)中，裂變包層中子主要分布于錒系與次錒系的共振能區(qū)，因此存在強烈的自屏效應(yīng)，如果中子學計算中忽略共振自屏效應(yīng)的影響，可能會低估系統(tǒng)的中子有效增殖因數(shù)keff。

實際應(yīng)用中，通常采用基于窄共振近似的Bondarenko方法[20]來進行共振自屏效應(yīng)修正，但這種方法不能有效處理共振干涉效應(yīng)，所以無法精確地對燃料類型復雜的先進核能系統(tǒng)包層進行能量自屏效應(yīng)修正。研究發(fā)現(xiàn)，當核素的參考本底截面選取的與實際堆型環(huán)境下的本底越接近，自屏效應(yīng)修正的結(jié)果就越精確，可以大幅度減少共振干涉效應(yīng)對自屏修正精度的影響[21]。因此通過對參考本底截面的優(yōu)化設(shè)計來實現(xiàn)能量自屏效應(yīng)的精確修正。

參考本底截面優(yōu)化法首先使用普適方法確認該核素的本底截面的大致范圍，再由典型反應(yīng)堆的設(shè)計方案計算出核素的本底截面，然后考慮反應(yīng)堆相關(guān)區(qū)域的不同情況下核素密度發(fā)生變化時對應(yīng)的本底截面范圍，最后綜合確定關(guān)鍵共振核素的本底截面優(yōu)化方案。

2.4 熱散射效應(yīng)修正

在先進核能系統(tǒng)設(shè)計研究中，水冷次臨界堆的能譜復雜，既有高能中子，又有大量熱中子，當中子能量低于一定范圍時(例如4 eV)，由于中子能量與靶核的熱能相當，中子在與靶核發(fā)生碰撞中可能獲得能量而使其能量升高，從而發(fā)生熱中子上散射效應(yīng)。已有研究表明，熱散射效應(yīng)對于超臨界水堆HMT1模型多群截面keff的影響可以達到14%。

因此需要專門針對熱散射修正問題開展研究，使HENDL可以適用于水冷次臨界堆的設(shè)計。熱中子散射截面通?？煞譃椋?) 相干彈性散射，如石墨、鈹?shù)龋?) 非彈性散射，包括相干散射和非相干散射，對所有材料都存在，用散射率S(α，β)表示；3) 非相干彈性散射，對含氫的固態(tài)物質(zhì)重要，如ZrH、固態(tài)輕水等。

熱散射數(shù)據(jù)的處理如下：以熱散射評價文件為原始輸入文件；在核數(shù)據(jù)處理程序NJOY[20]輸入文件中添加中子熱化處理卡片，實現(xiàn)熱能區(qū)材料中的熱散射矩陣計算，然后結(jié)合MATXSR模塊[20]的多群化處理，生成熱散射修正的多群數(shù)據(jù)文件。

2.5 溫度多普勒效應(yīng)修正

由于靶核的熱運動隨溫度的升高而增加，此時共振峰的寬度會隨著溫度上升而增加，同時峰值截面數(shù)值逐漸減小，這個現(xiàn)象稱為溫度多普勒效應(yīng)[20]。在核數(shù)據(jù)庫制作中，根據(jù)多普勒效應(yīng)對不同溫度的截面進行修正，稱為多普勒展寬。

先進核能系統(tǒng)內(nèi)材料的溫度分布范圍和梯度大，如裂變包層溫度高達數(shù)百攝氏度，而超導線圈和磁體的溫度則低到零下二百多攝氏度，因此HENDL核數(shù)據(jù)在設(shè)計中考慮多普勒效應(yīng)，從而保證核數(shù)據(jù)的精度，以適用于不同設(shè)計溫度下的包層計算分析。

針對核數(shù)據(jù)的多普勒展寬目前最為精確的方法是基于Kernal Broading[20]的多普勒展開方法，該方法從反應(yīng)率守恒和靶核熱運動速率的Maxwell-Boltzmann分布出發(fā)，推導出溫度多普勒展寬公式，并通過解析積分計算公式中的積分項，得到指定溫度下的經(jīng)過多普勒效應(yīng)修正的核反應(yīng)截面。

3 核數(shù)據(jù)驗證

FDS團隊針對混合評價核數(shù)據(jù)庫系統(tǒng)HENDL的驗證，匯總建立了一套系統(tǒng)化的核數(shù)據(jù)驗證例題庫。該例題庫包含2 000余道實驗與理論例題，在驗證核素方面包括錒系、次錒系、結(jié)構(gòu)材料、冷卻劑、功能材料等，驗證能量方面包括慢中子、快中子、高能中子等，較為系統(tǒng)地保證了測試的全面性。

通過熱譜/快譜臨界實驗、聚變中子源屏蔽積分實驗以及ITER-Benchmark等例題的測試，結(jié)果顯示HENDL的模擬計算數(shù)值與實驗值的符合較好[12-18]，驗證了HENDL在聚變級聚變驅(qū)動次臨界堆能譜下核數(shù)據(jù)的可靠性與準確性。通過快譜臨界實驗、高能中子屏蔽積分實驗以及BFS、KUCA、YALINA、BN600、IAEA/OECD-ADS等例題測試，結(jié)果顯示HENDL的模擬計算數(shù)值與實驗值的符合較好，驗證了HENDL在先進裂變堆、加速器驅(qū)動次臨界堆能譜下核數(shù)據(jù)的可靠性與準確性。本文選取聚變裝置停堆劑量率實驗例題ITER-T426與日本原子能機構(gòu)800 MW ADS的中子增殖因子測試案例進行介紹。

3.1 ITER-T426停堆劑量率實驗

在意大利ENEA的中子發(fā)生器FNG(Frascati Neutron Generator)上開展的ITER-T426[22]實驗，是最為典型的聚變裝置中子學計算程序與數(shù)據(jù)庫校驗測試實驗。該試驗采用FNG產(chǎn)生的D-T中子來照射實驗裝置，然后使用探測器測量不同冷卻時間該裝置中的劑量率，將實驗結(jié)果與程序計算結(jié)果進行比較分析，以校核相關(guān)計算程序、方法及核數(shù)據(jù)庫的可靠性[22]。ITER-T426模型由兩部分組成，中子產(chǎn)生器模型及立體屏蔽層模型，該裝置的模型，如圖1所示。

圖1 ITER-T426模型Fig.1 ITER-T426 model

應(yīng)用HENDL/MC輸運核數(shù)據(jù)庫與HENDL/MG嬗變/活化核數(shù)據(jù)庫，結(jié)合SuperMC計算程序?qū)TER-T426模型的停堆劑量率進行了計算，為了進一步驗證HENDL核數(shù)據(jù)庫的可靠性，并將計算結(jié)果與MCR2S、MCFISP等的計算結(jié)果[22]進行了對比。圖2給出了網(wǎng)格計數(shù)停機劑量率程序計算結(jié)果與實驗測量結(jié)果隨冷卻時間的變化。從該圖中可以看出，采用HENDL數(shù)據(jù)庫計算結(jié)果與實驗測量結(jié)果以及其他程序計算結(jié)果趨勢相吻合，除了初始時刻，HENDL計算結(jié)果與實驗值偏差都在5%左右，對比其他模擬結(jié)果，HENDL模擬計算結(jié)果更接近于實驗值。

圖2 計算值與實驗值的比值Fig.2 C/E calculation of dose rate in the cavity center

3.2 日本原子能機構(gòu)800 MW ADS測試例題

日本原子能機構(gòu)的800 MW加速器驅(qū)動次臨界堆ADS測試例題[23-24]，是隸屬于國際原子能機構(gòu)的核能聯(lián)合研究計劃，該計劃為“加速器驅(qū)動次臨界堆實驗與模擬測試”，其主要目的是為了研究ADS模擬工具的計算精度問題。該ADS的入射質(zhì)子能量為1.0 GeV，靶材為鉛鉍合金，燃料為次錒系钚與MA組成的混合物，如圖3所示。

圖3 RZ計算模型和燃料成分Fig.3 RZ calculation model and fuel composition

應(yīng)用HENDL-ADS/MG/MC核數(shù)據(jù)庫結(jié)合SuperMC計算程序?qū)υ摾}進行了模擬計算，并將計算結(jié)果與文獻[24]中JENDL-4.0、ENDF/B-VII.1、JEFF-3.1的計算結(jié)果進行比較，見表1所示。

表1 不同核數(shù)據(jù)庫模擬計算keff結(jié)果

* 相對偏差=(keff-平均值)/平均值×100

從計算結(jié)果來看，HENDL-ADS/MG/MC計算keff結(jié)果與平均值吻合較好，相對偏差小與0.05%。這些測試結(jié)果表明，HENDL核數(shù)據(jù)庫在計算加速器驅(qū)動次臨界堆ADS的次臨界度具備一定的精度與可靠性。

4 核數(shù)據(jù)庫應(yīng)用

目前混合評價核數(shù)據(jù)庫系統(tǒng)HENDL已經(jīng)成功應(yīng)用于國際熱核聚變實驗堆ITER、中國全超導托克馬克EAST、中國聚變示范堆CFETR、FDS系列反應(yīng)堆(FDS-I/II/III /SFB/MFX)以及加速器驅(qū)動次臨界堆等核能系統(tǒng)的設(shè)計與分析中。本文以聚變驅(qū)動次臨界堆FDS-SFB[25-30]的燃耗計算分析與ADS-CLEAR[31-32]的屏蔽分析為例介紹HENDL的應(yīng)用情況，計算所應(yīng)用的程序為超級蒙特卡羅核模擬軟件SuperMC[33-39]。

4.1 聚變驅(qū)動次臨界堆核設(shè)計分析

聚變驅(qū)動次臨界乏燃料焚燒堆FDS-SFB (Fusion Driven subcritical System for Spent Fuel Burning)是以焚燒裂變電站乏燃料為目的的聚變裂變混合堆概念，可實現(xiàn)嬗變核廢料、增殖核燃料、生產(chǎn)能量等功能。其聚變堆芯為常規(guī)托卡馬克裝置，大半徑為4 m，小半徑為1 m，聚變功率為150 MW。包層采用高壓氦氣和液態(tài)鉛鋰合金共同冷卻的雙冷嬗變包層，按照功能可以分為乏燃料區(qū)、氚增殖區(qū)與貧鈾區(qū)，其中乏燃料區(qū)主要用于產(chǎn)生能量和嬗變核廢料，氚增殖區(qū)主要用于增殖氚，貧鈾區(qū)主要用于增殖核燃料。

采用HENDL/FG輸運數(shù)據(jù)庫、嬗變/活化庫和材料輻照損傷庫結(jié)合SuperMC軟件，對FDS-SFB開展了優(yōu)化設(shè)計與分析。由于FDS-SFB的包層中子主要分布于錒系與次錒系的共振能區(qū)，因此存在強烈的自屏效應(yīng)，HENDL/FG采用參考本底截面的優(yōu)化[21]法以實現(xiàn)自屏效應(yīng)的精確修正。圖4顯示了采用HENDL/FG核數(shù)據(jù)庫計算的次錒系核素質(zhì)量隨燃耗的變化情況[40]。經(jīng)過分析與對比，超鈾含量為15%的包層燃料方案可以使得FDS-SFB獲得較高的燃料增殖與廢料嬗變性能。

圖4 次錒系核素質(zhì)量隨燃耗的變化曲線Fig.4 Mass of minor actinide as a function of burnup

4.2 中國鉛基研究實驗堆核設(shè)計分析

中國鉛基研究實驗堆CLEAR-I(China Lead-based Research Reactor)是中國科學院戰(zhàn)略性先導專項“未來先進核裂變能-ADS嬗變系統(tǒng)”的反應(yīng)堆參考方案。CLEAR-I具有臨界和次臨界雙模式運行能力，可在同一個裝置上開展第四代鉛冷快堆和外源驅(qū)動次臨界系統(tǒng)耦合技術(shù)研究。CLEAR-I選用富集度為19.75%的UO2為燃料，選擇鉛鉍合金作為冷卻劑，堆芯額定熱功率為10 MW。

基于HENDL-ADS/MG核數(shù)據(jù)庫結(jié)合SuperMC開展了CLEAR-I的堆芯核設(shè)計與屏蔽設(shè)計[41-42]，確定了堆芯基本布置與堆本體屏蔽方案。由于CLEAR-I在次臨界運行模式下系統(tǒng)中子能量跨度大、能譜復雜，經(jīng)過迭代優(yōu)化設(shè)計開發(fā)的HENDL-ADS/MG可以滿足大跨度能譜下物理計算需求。采用HENDL-ADS/MG計算了滿功率運行下堆芯功率分布、堆本體中子通量密度分布(見圖5)及關(guān)鍵部件(如包殼、圍桶、安全容器與主容器等)的材料輻照損傷，同時評估分析了堆頂包容小室的輻射來源及典型部件維修過程中的工作人員輻射劑量。計算結(jié)果為反應(yīng)堆設(shè)計方案優(yōu)化提供指導，最終堆芯與屏蔽設(shè)計方案得到了國家專家評審組的認可。

圖5 滿功率運行下CLEAR-I堆本體中子通量密度Fig.5 Neutron flux of CLEAR with full power operation

5 總結(jié)與展望

FDS團隊結(jié)合聚變堆、先進裂變堆、聚變驅(qū)動次臨界堆、加速器驅(qū)動次臨界堆等核能系統(tǒng)對應(yīng)用核數(shù)據(jù)的需求，設(shè)計開發(fā)出混合評價核數(shù)據(jù)庫HENDL。HENDL包括輸運核數(shù)據(jù)庫、嬗變與活化核數(shù)據(jù)庫、輻照損傷核數(shù)據(jù)庫等，并針對先進核能系統(tǒng)的物理特點，從能量自屏效應(yīng)、熱散射效應(yīng)、溫度多普勒效應(yīng)等方面進行了精確的截面修正。

HENDL已通過臨界安全實驗、積分屏蔽實驗、裂變堆基準實驗、聚變實驗堆基準模型等2 000余個國際基準模型與實驗的驗證與確認，并在ADS-CLEAR 、FDS-SFB為代表的 20 余個重大核工程與核能研究中得到了廣泛應(yīng)用。

下一步將繼續(xù)深入研究HENDL系列核數(shù)據(jù)庫在聚變堆、先進裂變堆、ADS等先進核能系統(tǒng)核設(shè)計與安全分析方面的應(yīng)用。同時，針對先進核能系統(tǒng)的特點，有計劃地開展針對核數(shù)據(jù)的微觀與宏觀實驗，建立HENDL的實驗數(shù)據(jù)庫，從根本上提高核數(shù)據(jù)的精度，結(jié)合數(shù)據(jù)的理論評價工作，全面充實與提高HENDL核數(shù)據(jù)庫的數(shù)據(jù)內(nèi)容與質(zhì)量。

[1] Otuka N, Dupont E, Semkova V, et al. Towards a more complete and accurate experimental nuclear reaction data library (EXFOR) [J]. Nuclear Data Sheets, 2014, 120: 272-276.

[2] Chadwick M B, Herman M, et al. ENDF/B-VII. 0: Next generation evaluated nuclear data library for nuclear science and technology[J]. Nuclear data sheets, 2006, 107(12): 2931-3060.

[3] Santamarina A, Bernard D, Blaise P, et al. The JEFF-3.1.1 nuclear data library[R]. JEFF report, 2009.

[4] Lemmel H D, McLaughlin P K. BROND-2.2 Russian evaluated neutron reaction data library[R]. International Atomic Energy Agency, 1994.

[5] Shibata K, Kawano T, Nakagawa T, et al. Japanese evaluated nuclear data library version 3 revision-3:JENDL-3.3[J]. Journal of Nuclear Science and Technology, 2002, 39: 1125-1136.

[6] Ge Z G, Yu H W, Zhuang Y X, et al. The updated version of the Chinese evaluated nuclear data Library (CENDL-3.1) and China nuclear data evaluation activities[J]. EDP Sciences, 2007, 753-757.

[7] D. López Aldama, F. Leszczynski, A. Trkov. WIMSD-D Library Update[R]. Document INDC(NDS)-467. International Atomic Energy Agency, 2003.

[8] D. López Aldama, et al. FENDL2.1—Update of an evaluated nuclear data library for fusion applications[R]. Document INDC(NDS)-467. International Atomic Energy Agency, 2004.

[9] D. López Aldama, et al. ADS-2.0——A Test Library for Accelerator Driven Systems and New Reactor Designs[R]. Document INDC(NDS)-0545. International Atomic Energy Agency, 2008.

[10] M.B. Chadwick, et al., LA150 of Cross Sections, Heating, and Damage: A (Incident Neutrons) and Part B(Incident Protons), Nuclear Science and Engineering, 293, 131, 1999.

[11] Y. Watanabe, K. Kosako,S. Kunieda, et al, Status of JENDL High Energy File, Journal of the Korean Physical Society,59:1040-1045, 2011.

[12] Y Wu, U. Fischer. Integral Data Tests of the FENDL-1, EFF-2, EFF-3 and JENDL-FF Fusion Nuclear Data Libraries[R]. Forschungszentrum Karlsruhe, Technik and Umwelt Wissenschaftliche Berichte FZKA 5953, Februar 1998.

[13] Yican Wu,Yixue Chen, U. Fischer,et al. Integral Data Tests of FENDL-2 Fusion Nuclear Data Library with Neutronic Integral Experiments,Journal of Nuclear Science and Technology,supplement 1(March 2000).

[14] Gao Chunjing, Xu Dezheng, Li Jingjing, et al. Integral Data Test of HENDL1.0/MG and VisualBUS with Neutronics Shielding Experiments (I)[J]. Wu Yican, Deng Tieru. Plasma Science and Technology, 2004, 6(5): 2507-2513.

[15] Gao Chunjing, Xu Dezheng, Li Jingjing, et al.Integral Data Test of HENDLl.0/MG with Neutronics Shielding Experiments (II)[J]. Plasma Science and Technology, 2004, 6(6): 2596-2600.

[16] Dezheng Xu, Zhaozhong He, Jun Zou, et al. Production and testing of HENDL-2.1/CG coarse-group cross-section library based on ENDF/B-VII.0[J], Fusion Engineering and Design, 85 (2010) : 2105-2110.

[17] Jun Zou, Qin Zeng, Dezheng Xu, et al. Design and Production of Fine-group Cross Section Library HENDL3.0/FG Subcritical System[J], International Conference on the Physics of Reactors 2012, PHYSR 2012: Advances in Reactor Physics,v5,p 5483-5490,2012.

[18] Jun Zou, Qin Zeng,Chong Chen,et al, Development and Application of HENDL-ADS/MG Cross Section Library for ADS System[J],Transactions—American Nuclear Society;108;883-886(2013).

[19] Hanshaw Heath L, et al. Estimation of the Effects of Self-Shielding on Multigroup Reactor Vessel Fluence Calculation [J], Trans.Am.Nucl. Soc, 71 :386-388, 1994.

[20] R.E.MacFarlane, et al. The NJOY Nuclear Data Processing System Version 91[R]. LA-12740-M, Los Alamos National Laboratory, 1994.

[21] 鄒俊，用于次臨界堆設(shè)計分析的核數(shù)據(jù)庫系統(tǒng)研發(fā)與應(yīng)用，博士，中國科學院等離子體物理研究所，2010.

[22] P. Batistoni, L. Petrizzi, Task T426 - Neutronics Experiments, Experimental Validation of Shut Down Dose Rates, EFF-Doc-726, March 2000.

[23] K. Nishihara, T. Sugawara, H. Iwamoto, et al. Investigation of Nuclear Data Accuracy for the Accelerator-Driven System with Minor Actinide Fuel, OECD/NEA 11th Information Exchange Meeting, San Francisco (2010) 315.

[24] Sugawara T, Zolbadral T, Nishihara K, et al. Accelerator-Driven System Analysis by Using Different Nuclear Data Libraries[J]. Atom Indonesia, 2013, 38(2).

[25] Y. Wu, FDS Team. Conceptual Design Activities of FDS Series Fusion Power Plants in China. Fusion Engineering and Design, 2006, 81(23-24): 2713-2718.

[26] L. Qiu, Y. Wu, B. Xiao, et al. A Low Aspect Ratio Tokamak Transmutation System. Nuclear Fusion, 2000, 40: 629-633.

[27] Y. Wu, J. Qian, J. Yu.The Fusion-Driven Hybrid System and Its Material Selection. Journal of Nuclear Materials, 2002, 307-311:1629-1636.

[28] Y. Wu, J. Jiang, M. Wang, et al. A Fusion-Driven Subcritical System Concept Based on Viable Technologies. Nuclear Fusion, 2011, 51(10):103036.

[29] Y. Wu, FDS Team.Conceptual Design of the China Fusion Power Plant FDS-II. Fusion Engineering and Design, 2008, 83(10-12): 1683-1689.

[30] Y. Wu, FDS Team.Fusion-Based Hydrogen Production Reactor and Its Material Selection. Journal of Nuclear Materials, 2009, 386-388:122-126.

[31] Y.Wu, Y.Bai, W.Wang, et.al, “Overview of China Lead Alloy cooled Reactor development and ADS Program in China”, NUTHOS-9, Kaohsiung, Taiwan, September 9-13,2012.

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

[33] Y. Wu, J. Li, Y. Li, et al., An integrated multi-functional neutronics calculation and analysis code system: VisualBUS, Chinese Journal of Nuclear Science and Engineering, 27 (2007) 365-373.

[34] Yican Wu, X.George Xu, The Need for Further Development of CAD/MCNP Interface Codes, Invited presentation at the American Nuclear Society Annual Meeting, Transactions of American Nuclear Society, 2007, 96: TRANSAO 961-882, ISSN: 0003-018X, June 24-28, 2007, Boston, Massachusetts.

[35] Y. Wu, FDS Team, CAD-based Interface Programs for Fusion Neutron Transport Simulation. Fusion Engineering and Design, 2009, 84:1987-1992.

[36] Y. Li, L. Lu, A. Ding, et al. Benchmarking of MCAM 4.0 with the ITER 3D model, Fusion Engineering and Design 82 (2007) 2861-2866. Also presented at the 24th Symposium on Fusion Technology (SOFT-24), Warsaw, Sep.11-15, 2006.

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

[38] Y.Wu, J. Song, H.Q. Zheng, et al. CAD-based Monte Carlo Program for Integrated Simulation of Nuclear System SuperMC. Annals of Nuclear Energy, DOI: 10.1016/j.anucene.2014.08.058.

[39] J. Song, G.Y. Sun, Z.P. Chen, et al. Benchmarking of CAD-based SuperMC with ITER Benchmark Model. Fusion Engineering and Design, DOI: 10.1016/j.fusengdes.2014.05.003.

[40] 王明煌，次臨界乏燃料焚燒堆中子學設(shè)計研究，博士，中國科學技術(shù)大學，2012.

[41] Qi Yang, Bin Li, Chao Chen, et al. Radiation Shielding Design and Analysis for CLEAR-I [J], Paper No. ICONE21-15995, pp. V001T04A022; 5pages.

[42] Bin Li, Qi Yang, Jou Zou, et al. Shielding design and analysis for the proton tube of China LEAD-ALLOY Cooled Reasearch Reactor(CLEAR-I), Paper No. ICONE21-15763, pp. V001T04A020; 5pages.

Development and Application of Hybrid Evaluated Nuclear Data Library HENDL3.0 for Advanced Nuclear System

WU Yi-can, ZOU Jun, HAO Li-juan, WANG Ming-huang, YANG Qi, SONG Jing, WANG Jin, SHANG Lei-ming, LONG Peng-cheng, WANG Fang, HU Li-qin, HE Tao, FDS Team

(Key Laboratory of Neutronics and Radiation Safety, Institute of Nuclear Energy Safety Technology, Chinese Academy of Sciences, Hefei, Anhui Prov. 230031, China)

To meet the need of nuclear analysis for fusion reactors, fission reactors, the fusion-driven sub-critical and accelerator-driven sub-critical systems, HENDL (Hybrid Evaluated Nuclear Data Library) has been developed by FDS team. Corrections of resonance self-shielding, Doppler and thermal scattering effects for advanced nuclear system were accomplished in the design for HENDL. HENDL include different energy structures and weight functions for different reactors. To validate and qualify the reliability of the HENDL, lots of benchmarks were performed. Now HENDL has been applied in fission reactors, fusion reactors, FDS serial reactors and ADS-CLEAR systems.

Advanced nuclear systems; Multi-group data library; Benchmark

2017-02-11

國家自然科學基金(11405204、11305205、10675123)，中國科學院戰(zhàn)略先導專項(XDA03040000)，國家ITER 973計劃(2014GB112001)，國家自然科學基金(11405204、11305205、10675123)，中國科學院信息化專項(XXH12504-1-09)，產(chǎn)業(yè)化基金

吳宜燦(1964—)，男，安徽人，研究員，博導，主要從事先進核能系統(tǒng)研究

TL32

A

0258-0918(2017)02-0242-08