基于WIMS-D庫燃耗數(shù)據(jù)的更新與基準(zhǔn)驗(yàn)證

肖 向 吳 軍 龔 游 楊仝瑞 董正云 陳義學(xué),

1（華北電力大學(xué) 核科學(xué)與工程學(xué)院 北京 102206）

2（國家核電技術(shù)有限公司 北京軟件技術(shù)中心 北京 100029）

3（中國科學(xué)院核能安全技術(shù)研究所 合肥 230000）

基于WIMS-D庫燃耗數(shù)據(jù)的更新與基準(zhǔn)驗(yàn)證

肖 向1吳 軍1龔 游1楊仝瑞2董正云3陳義學(xué)1,2

1（華北電力大學(xué) 核科學(xué)與工程學(xué)院 北京 102206）

2（國家核電技術(shù)有限公司 北京軟件技術(shù)中心 北京 100029）

3（中國科學(xué)院核能安全技術(shù)研究所 合肥 230000）

隨著AP1000等新一代壓水堆的發(fā)展，燃耗深度在不斷提高，平均卸料燃耗深度提高到60 GW d·t?1。然而，傳統(tǒng)使用的WIMS69群和XMAS172群WIMS-D格式多群常數(shù)庫，其能群結(jié)構(gòu)存在共振峰重疊，核素種類較少，裂變產(chǎn)物產(chǎn)額的偏差較大，并且偽裂變產(chǎn)物中包含的核素種類較多而導(dǎo)致152Gd、160Gd、159Tb、160Tb等重要核素?zé)o法得到精確處理等問題。因此，本文主要針對(duì)AP1000等新一代反應(yīng)堆的設(shè)計(jì)以及運(yùn)行特點(diǎn)，基于ENDF/B-VII.1庫，并且在現(xiàn)有基礎(chǔ)上針對(duì)WIMS-D庫中的偽裂變產(chǎn)物、裂變產(chǎn)物燃耗鏈以及裂變產(chǎn)物產(chǎn)額等燃耗數(shù)據(jù)進(jìn)行更新，再通過NJOY程序開發(fā)了SHEM281群WIMS-D格式多群常數(shù)庫。通過DRAGON程序掛載該WIMSD281庫，對(duì)其進(jìn)行臨界和燃耗兩方面基準(zhǔn)驗(yàn)證。計(jì)算結(jié)果表明，該數(shù)據(jù)庫的計(jì)算結(jié)果與基準(zhǔn)值符合較好，精度較高，結(jié)果可靠，可初步用于壓水堆的相關(guān)計(jì)算。

SHEM281，偽裂變產(chǎn)物，裂變產(chǎn)物鏈，裂變產(chǎn)物產(chǎn)額

目前廣泛使用的堆芯物理庫是國際原子能機(jī)構(gòu)(International Atomic Energy Agency, IAEA)在2007年發(fā)布，基于 WIMS69群和 XMAS172群的WIMS-D庫，其數(shù)據(jù)來源主要是 ENDF/B-VI.7、JENDL-3.2、CENDL-3、JEFF-2.2等。XMAS172能群結(jié)構(gòu)克服了 WIMS69[1]的一些缺點(diǎn)，更適用于MOX (Mixed Oxide Fuel)燃料和深燃耗問題等計(jì)算。但是XMAS172沒有考慮共振峰之間的相互影響，限制了計(jì)算精度的進(jìn)一步提高。因此，為了克服上述兩種能群結(jié)構(gòu)的缺點(diǎn)，法國原子能機(jī)構(gòu)以及加拿大蒙特利爾大學(xué)提出了SHEM281[2]的能群結(jié)構(gòu)，它能精確描述錒系核素及其主要裂變產(chǎn)物在23 eV以下的共振峰，減小由于共振峰重疊而引入的誤差。

隨著燃耗深度的提高，裂變產(chǎn)物積累較多，各產(chǎn)物之間的轉(zhuǎn)化過程變得更為復(fù)雜，IAEA在2003年發(fā)布的WIMS-D Library Update[3]中提出，燃耗數(shù)據(jù)制作方法會(huì)引入較大誤差，并且該燃耗數(shù)據(jù)的評(píng)價(jià)庫來源較為陳舊，限制了計(jì)算精度的進(jìn)一步提高。

本文主要基于 ENDF/B-VII.1庫，采用SHEM281能群結(jié)構(gòu)，結(jié)合AP1000等新一代反應(yīng)堆的特點(diǎn)，通過截面處理程序NJOY[4]，制作WIMS-D格式的多群常數(shù)庫。目前，已制作好的WIMSD281庫已經(jīng)進(jìn)行了金屬鈾、氧化鈾和MOX燃料臨界基準(zhǔn)題的驗(yàn)證[5]，計(jì)算精度較高，結(jié)果可靠。本文主要針對(duì)該數(shù)據(jù)庫中的燃耗數(shù)據(jù)進(jìn)行更新，最后通過DRAGON程序[6]對(duì)該數(shù)據(jù)庫進(jìn)行進(jìn)一步臨界和燃耗驗(yàn)證，基準(zhǔn)題的選取主要針對(duì)AP1000反應(yīng)堆。

1 燃耗數(shù)據(jù)更新

在WIMS庫中，燃耗數(shù)據(jù)是每種核素在發(fā)生核反應(yīng)時(shí)的變化，主要包含：1) 錒系核素的裂變能；2) 錒系核素、裂變產(chǎn)物以及可燃毒物的燃耗鏈；3)裂變產(chǎn)物和偽裂變產(chǎn)物的ID及其產(chǎn)額。在IAEA提供的WIMS69群庫中，包含燃耗數(shù)據(jù)的核素一共有：57個(gè)裂變產(chǎn)物、3個(gè)偽裂變產(chǎn)物（1個(gè)平均的偽裂變產(chǎn)物、1個(gè)232U的偽裂變產(chǎn)物和1個(gè)237U的偽裂變產(chǎn)物）、7個(gè)可燃毒物、24個(gè)錒系核素。本文主要針對(duì)偽裂變產(chǎn)物的重構(gòu)、裂變產(chǎn)物燃耗鏈的優(yōu)化及裂變產(chǎn)額的計(jì)算進(jìn)行著重闡述。

1.1 偽裂變產(chǎn)物重構(gòu)

錒系核素在發(fā)生裂變反應(yīng)時(shí)往往會(huì)生成上千種裂變產(chǎn)物，其中一部分核素對(duì)反應(yīng)性有顯著影響。另一部分核素單獨(dú)作用對(duì)反應(yīng)性的影響不大，但是其整體對(duì)反應(yīng)性的影響卻是不可忽略的。在工程計(jì)算中，習(xí)慣將后者作為一種“裂變產(chǎn)物”來處理，一般將其稱為假想集總裂變產(chǎn)物，俗稱偽裂變產(chǎn)物[7]。

隨著AP1000等新一代反應(yīng)堆的發(fā)展，新型低價(jià)值控制棒中子吸收體材料主要采用鋱(Tb)和鏑(Dy)這兩種核素制成[8]，以及采用釓(Gd)作為可燃毒物棒進(jìn)行反應(yīng)性控制[9]。然而，IAEA提供的偽裂變產(chǎn)物中包含79個(gè)核素，但是這些核素卻包含著目前AP1000等新一代反應(yīng)堆需要著重考慮的核素，如159Tb、160Tb、152Gd、160Gd這4個(gè)核素。因此，需要針對(duì)現(xiàn)階段反應(yīng)堆的運(yùn)行特點(diǎn)，重新制作一版?zhèn)瘟炎儺a(chǎn)物。

圖1是采用SHEM281能群結(jié)構(gòu)，按照IAEA提供的核素而制作的235U、238U、239Pu下，對(duì)應(yīng)偽裂變產(chǎn)物之間截面對(duì)比圖。通過分析可以看出，238U的偽裂變產(chǎn)物在1?100 eV內(nèi)有很大的偏差，而通過這樣的偽裂變產(chǎn)物制作出來的平均偽裂變產(chǎn)物，會(huì)對(duì)計(jì)算結(jié)果引入較大的誤差。

圖1 不同偽裂變產(chǎn)物核素輸運(yùn)截面(a)和吸收截面(b)對(duì)比Fig.1 Transport cross section (a) and absorption cross section (b) of the different old pseudo fission products.

圖2是采用SHEM281能群結(jié)果，在IAEA的基礎(chǔ)上去掉159Tb、160Tb、152Gd、160Gd這4個(gè)核素后，制作的235U、238U、239Pu下對(duì)應(yīng)偽裂變產(chǎn)物之間截面對(duì)比圖。通過分析可以看出，對(duì)比IAEA制作出來的偽裂變產(chǎn)物來說，去除這4個(gè)核素后制作出來的偽裂變產(chǎn)物之間的截面偏差較小。因此，去除159Tb、160Tb、152Gd、160Gd這4個(gè)核素后制作出來的平均偽裂變產(chǎn)物，能在一定程度上減小誤差。

圖2 新不同偽裂變核素輸運(yùn)截面(a)和吸收截面(b)對(duì)比Fig.2 Transport cross section (a) and absorption cross section (b) of the different new pseudo fission products.

1.2 裂變產(chǎn)物燃耗鏈優(yōu)化

在核反應(yīng)堆堆芯中發(fā)生核裂變將產(chǎn)生大量的裂變產(chǎn)物，部分核素的裂變產(chǎn)物有上千種。如果將每一種裂變產(chǎn)物都進(jìn)行精確計(jì)算，則會(huì)花費(fèi)大量時(shí)間，這在工程計(jì)算中是不可行的。因此，在工程計(jì)算中只選擇對(duì)反應(yīng)性影響較大的裂變產(chǎn)物做精確計(jì)算。然而，由于在之前的偽裂變產(chǎn)物中減少了部分核素，所以需要針對(duì)這幾個(gè)額外添加的裂變產(chǎn)物制作出一個(gè)新的裂變產(chǎn)物鏈。由于額外增加的裂變產(chǎn)物為：159Tb、160Tb、152Gd、160Gd，因此，只要制作與這4個(gè)核素有關(guān)的裂變產(chǎn)物鏈即可。目前，這4個(gè)核素相關(guān)的完整燃耗鏈如圖3所示。由于ENDF/B-VII.1庫缺少159Gd、161Gd、161Tb相關(guān)的評(píng)價(jià)核數(shù)據(jù)，所以本文在完整燃耗鏈的基礎(chǔ)上通過簡化省略得到一個(gè)新裂變產(chǎn)物燃耗鏈，如圖4所示。

圖3 完整燃耗鏈Fig.3 The complete burnup chain.

圖4 簡化燃耗鏈Fig.4 The simplified burnup chain.

1.3 裂變產(chǎn)物產(chǎn)額計(jì)算

裂變產(chǎn)額與入射中子能量有關(guān)，隨著入射中子能量變化，裂變產(chǎn)額會(huì)隨之改變。WIMS庫采用的是平均裂變產(chǎn)額，通過 AVRFPY程序[10]提取。AVRFPY程序在計(jì)算平均產(chǎn)額過程中，將中子能群劃分為三群，其能量范圍分別為：10?5?5.5 eV、5.5eV?0.1 MeV、0.1?10 MeV。再根據(jù)裂變反應(yīng)率進(jìn)行加權(quán)求解，具體的裂變產(chǎn)物提取方法見圖5。

然而IAEA69和IAEA172兩個(gè)庫的數(shù)據(jù)來源為ENDF/B-VI.8庫，使得其計(jì)算出來的裂變產(chǎn)額數(shù)據(jù)過于陳舊，會(huì)對(duì)最終計(jì)算結(jié)果引入一定的偏差。隨著 ENDF/B-VII.1庫的發(fā)布，需要針對(duì)新的評(píng)價(jià)核數(shù)據(jù)庫計(jì)算出一版新的裂變產(chǎn)額數(shù)據(jù)。圖6和7主要表示通過兩種不同評(píng)價(jià)核數(shù)據(jù)庫計(jì)算出來的235U裂變產(chǎn)額以及相對(duì)偏差，通過分析可以知道，不同數(shù)據(jù)庫計(jì)算出來的裂變產(chǎn)額偏差較大，部分偏差甚至達(dá)到了100%。因此，在新的ENDF/B-VII.1庫基礎(chǔ)上，通過AVRFPY程序計(jì)算出來的新裂變產(chǎn)額替換IAEA提供的舊裂變產(chǎn)額是非常有必要的。

圖5 裂變產(chǎn)物產(chǎn)額的提取流程Fig.5 The process of extracting the fission product yield.

圖6 不同評(píng)價(jià)庫下235U的裂變產(chǎn)物產(chǎn)額Fig.6 The 235U fission product yield of the different nuclear evaluation libraries.

圖7 235U裂變產(chǎn)物產(chǎn)額的相對(duì)偏差Fig.7 The relative deviation of 235U fission product yield.

2 基準(zhǔn)驗(yàn)證

群常數(shù)庫是在核評(píng)價(jià)庫的基礎(chǔ)上，經(jīng)過群常數(shù)制作程序處理得到的。然而，此時(shí)的群常數(shù)庫是未經(jīng)驗(yàn)證以及確認(rèn)的，尚不能用于工程計(jì)算與科學(xué)研究。因此，只有經(jīng)過驗(yàn)證和確認(rèn)的群常數(shù)才能用于工程設(shè)計(jì)與堆物理計(jì)算，本文主要從臨界計(jì)算和燃耗計(jì)算兩方面進(jìn)行驗(yàn)證。

由于WIMSD281庫主要是為AP1000系列提供設(shè)計(jì)和研究用的多群常數(shù)庫，所以基準(zhǔn)題的驗(yàn)證主要針對(duì)AP1000核反應(yīng)堆。目前AP1000使用的燃料組件[11]主要有兩種：一種是由264根燃料棒、24根導(dǎo)向管以及1根中子測量管組成；另一種則是由IFBA (Intergral Fuel Burnable Absorber)棒、燃料棒以及WABA (Wet Annular Burnable Absorber)棒組成。因此，本文主要針對(duì)AP1000中的典型燃料棒柵元以及IFBA柵元為代表，進(jìn)行基準(zhǔn)驗(yàn)證。

2.1 臨界基準(zhǔn)驗(yàn)證

臨界基準(zhǔn)驗(yàn)證主要反映著反應(yīng)堆距離臨界的程度，選取的兩種柵元結(jié)構(gòu)均采用UO2燃料，幾何模型如圖8、9所示，富集度分別為：0.74%、1.58%、3.2%、3.4%、3.8%、4.0%、4.2%、4.4%、4.8%。通過 DRAGON 程序掛載 WIMSD281、IAEA69、IAEA172多群常數(shù)庫進(jìn)行計(jì)算，并且將其計(jì)算結(jié)果與 MCNP程序[12]進(jìn)行對(duì)比。由表 1可知，在算AP1000典型柵元的kinf時(shí)，WIMSD281、IAEA69、IAEA172數(shù)據(jù)庫的平均偏差分別為?0.00654、?0.00786、?0.00469；由表2可知，在算 IFBA柵元的kinf時(shí)，WIMSD281、IAEA69、IAEA172數(shù)據(jù)庫的平均偏差分別為?0.0047、?0.00352、?0.00138。

臨界計(jì)算的過程中不涉及燃耗數(shù)據(jù)，僅和多群常數(shù)庫中的截面信息有關(guān)，并且這兩種柵元結(jié)構(gòu)介質(zhì)成分簡單，由于共振峰干涉效應(yīng)所引入的誤差相對(duì)較小，所以偏差主要來自于235U和238U的截面數(shù)據(jù)。通過分析可知，三種數(shù)據(jù)庫的計(jì)算結(jié)果相對(duì)MCNP整體偏小，誤差可能來自兩個(gè)方面：1)DRAGON程序采用 CPM (Collision Probability Method)計(jì)算過程中引入了一些簡化近似導(dǎo)致，如柵元的等效、平通量密度近似等[13]；2) MCNP程序采用的評(píng)價(jià)核數(shù)據(jù)庫是 ENDF/B-VII.0庫，而WIMSD281庫采用的評(píng)價(jià)核數(shù)據(jù)庫是ENDF/B-VII.1庫，IAEA69和IAEA172庫則是基于ENDF/B-VI.8庫制作而成，這可能是不同評(píng)價(jià)核數(shù)據(jù)庫之間的偏差所導(dǎo)致。

另一方面，IAEA172的精度要比 WIMSD281高0.002左右，這可能是由于IAEA采用NRSC程序[14]針對(duì)WIMS69和IAEA172能群結(jié)構(gòu)修正了高能區(qū)和共振能區(qū)的有效共振積分，從而導(dǎo)致其計(jì)算精度要高于能群結(jié)構(gòu)更加優(yōu)化的SHEM281。因此，針對(duì)WIMSD281庫的截面數(shù)據(jù)，還需要進(jìn)行進(jìn)一步研究改進(jìn)。

圖8 AP1000典型柵元(a)與IFBA柵元(b)幾何模型Fig.8 Geometry model of AP1000 typical pin (a)and IFBA pin (b).

表1 AP1000典型柵元的kinfTable 1 The kinf of AP1000 typical pin.

表2 IFBA柵元的kinfTable 2 The kinf of IFBA pin.

2.2 燃耗基準(zhǔn)驗(yàn)證

燃耗基準(zhǔn)驗(yàn)證則主要是計(jì)算在不同燃耗點(diǎn)下的kinf以及柵元內(nèi)重要核素隨燃耗深度的變化。本文同樣針對(duì)AP1000典型柵元和IFBA柵元進(jìn)行燃耗基準(zhǔn)驗(yàn)證，幾何模型如圖8所示，僅選擇富集度為3.2%和4.8%兩種情況進(jìn)行考慮。通過DRAGON程序掛載 WIMSD281、IAEA69、IAEA172多群常數(shù)庫進(jìn)行計(jì)算，并將其計(jì)算結(jié)果與SCALE程序[15]進(jìn)行對(duì)比。SCALE程序是由美國橡樹嶺國家實(shí)驗(yàn)室開發(fā)的一個(gè)模塊化程序系統(tǒng)，包含多種不同計(jì)算任務(wù)的程序模塊，主要用于核臨界安全分析、反應(yīng)堆物理計(jì)算和輻射屏蔽計(jì)算等。采用 SCALE程序系統(tǒng)中TRITON控制模塊，選取基于ENDF/B-VI的238群截面庫，調(diào)用KENO輸運(yùn)模塊進(jìn)行輸運(yùn)計(jì)算，再調(diào)用ORIGEN-S燃耗模塊針對(duì)柵元、組件等進(jìn)行燃耗計(jì)算。

圖9和10主要展示了AP1000典型柵元和IFBA柵元在不同富集度中kinf隨燃耗的變化，AP1000典型柵元在燃耗壽期初時(shí)，WIMSD281的計(jì)算結(jié)果與SCALE符合得更好，而IAEA69與IAEA172兩個(gè)數(shù)據(jù)庫的計(jì)算結(jié)果與SCALE偏差較大；IFBA柵元在整個(gè)燃耗壽期內(nèi)，這三種數(shù)據(jù)庫的計(jì)算結(jié)果差別不大，并且隨著燃耗深度的增加，結(jié)果與 SCALE偏差越來越大。這主要是因?yàn)樵谏钊己膮^(qū)內(nèi)：1) 裂變產(chǎn)物較多，可能需要考慮裂變產(chǎn)物的共振吸收；2) 在燃耗過程中，中子能譜是發(fā)生變化的，采用怎樣的能譜制作燃耗截面也對(duì)燃耗結(jié)果有較大影響；3) 中子通量密度的間接影響等[16]。這是在深燃耗區(qū)計(jì)算時(shí)，所共同面臨的主要問題，本文暫不對(duì)此進(jìn)行考慮。

圖9 AP1000典型柵元在燃料富集度為3.2% (a)和4.8% (b)時(shí)kinf隨燃耗的變化Fig.9 The kinf result of AP1000 typical pin between the fuel enrichment of 3.2% (a) and 4.8% (b).

圖10 IFBA柵元在燃料富集度為3.2% (a)和4.8% (b)時(shí)kinf隨燃耗的變化Fig.10 The kinf result of IFBA typical pin between the fuel enrichment of 3.2% (a) and 4.8% (b).

圖11 和12主要展示了AP1000典型柵元在不同富集度中主要錒系核素核子密度隨燃耗的變化，在計(jì)算235U和238U的核子密度時(shí)，各數(shù)據(jù)庫的計(jì)算結(jié)果與SCALE符合較好。在計(jì)算239Pu的核子密度時(shí)，WIMSD281的計(jì)算結(jié)果與SCALE符合得更好，但是在計(jì)算240Pu的核子密度時(shí)，各數(shù)據(jù)庫的結(jié)果在壽期初時(shí)與 SCALE符合較好，但隨著燃耗深度的增加，各結(jié)果相對(duì) SCALE有些偏小。這同樣是深燃耗區(qū)計(jì)算所共同面臨的主要問題導(dǎo)致。

圖13和14主要展示了IFBA柵元在不同富集度中主要錒系核素核子密度隨燃耗的變化，在計(jì)算235U、238U、240Pu的核子密度時(shí)，三種數(shù)據(jù)庫的結(jié)果與SCALE符合較好；在計(jì)算239Pu的核子密度時(shí)，WIMSD281的計(jì)算結(jié)果與SCALE符合得更好。通過對(duì)比分析可知，WIMSD281庫在計(jì)算239Pu和240Pu這兩種核素時(shí)，其精度要比 IAEA69和IAEA172高。然而在深燃耗區(qū)內(nèi)，IAEA69、IAEA172和WIMSD281的計(jì)算結(jié)果與SCALE偏差較大。由此可見，在深燃耗區(qū)發(fā)生的核反應(yīng)過程多而復(fù)雜，采用IAEA在文獻(xiàn)[3]中提出的燃耗數(shù)據(jù)會(huì)產(chǎn)生較多的問題。如果用于核素種類更多的MOX燃料計(jì)算中，可能會(huì)產(chǎn)生更大的偏差。因此，目前開發(fā)的WIMSD281多群截面庫在深燃耗區(qū)計(jì)算仍具有一定的限制，需要進(jìn)一步改進(jìn)。

圖11 AP1000典型柵元在燃料富集度為3.2%時(shí)，235U (a)、238U (b)、239Pu (c)、240Pu (d)核子密度隨燃耗的變化Fig.11 The number density of 235U (a), 238U (b), 239Pu (c), 240Pu (d) with the fuel enrichment of 3.2% in AP1000 typical pin.

圖12 AP1000典型柵元在富集度為4.8%時(shí)，235U (a)、238U (b)、239Pu (c)、240Pu (d)核子密度隨燃耗的變化Fig.12 The number density of 235U (a), 238U (b), 239Pu (c), 240Pu (d) with the fuel enrichment of 4.8% in AP1000 typical pin.

圖13 IFBA柵元在富集度為3.2%時(shí)，235U (a)、238U (b)、239Pu (c)、240Pu (d)核子密度隨燃耗的變化Fig.13 The number density of 235U (a), 238U (b), 239Pu (c), 240Pu (d) with the fuel enrichment of 3.2% in IFBA pin.

圖14 IFBA柵元在富集度為4.8%時(shí)，235U (a)、238U (b)、239Pu (c)、240Pu (d)核子密度隨燃耗的變化Fig.14 The number density of 235U (a), 238U (b), 239Pu (c), 240Pu (d) with the fuel enrichment of 4.8% in IFBA pin.

3 結(jié)語

本文主要針對(duì)AP1000等新一代反應(yīng)堆的運(yùn)行以及設(shè)計(jì)特點(diǎn)，在 IAEA提供的 WIMS69群和XMAS172群庫基礎(chǔ)上，采用SHEM281能群制作出來的一版 WIMS-D格式的多群常數(shù)庫并且還將其燃耗數(shù)據(jù)進(jìn)行更新，主要針對(duì)IAEA在文獻(xiàn)[3]中提出的偽裂變產(chǎn)物進(jìn)行重構(gòu)，把之前的裂變產(chǎn)物燃耗鏈進(jìn)行優(yōu)化，在 ENDF/B-VII.1庫的基礎(chǔ)上重新計(jì)算裂變產(chǎn)物產(chǎn)額。本文還針對(duì)AP1000典型柵元以及 IFBA柵元，通過 DRAGON程序掛載WIMSD281、IAEA69、IAEA172這三個(gè)多群常數(shù)庫進(jìn)行臨界和燃耗基準(zhǔn)驗(yàn)證，并將其計(jì)算結(jié)果與MCNP、SCALE進(jìn)行對(duì)比。通過臨界基準(zhǔn)驗(yàn)證可知，DRAGON程序掛載 WIMSD281庫計(jì)算出來的 kinf相比MCNP程序有些偏小，這可能是由于評(píng)價(jià)核數(shù)據(jù)庫中238U的吸收截面偏大所導(dǎo)致。通過燃耗基準(zhǔn)驗(yàn)證可知，除了深燃耗區(qū)，DRAGON程序掛載WIMSD281庫的結(jié)果與SCALE程序相比，計(jì)算精度較高，結(jié)果可靠，能適用于AP1000等新一代反應(yīng)堆的研究和設(shè)計(jì)。

1 DAN G C, Robert C B, Yaron D, et al. Handbook of nuclear engineering: Vol.1: nuclear engineering fundamentals; Vol.2: reactor design; Vol.3: reactor analysis; Vol.4: reactors of generations III and IV; Vol.5:fuel cycles, decommissioning, waste disposal and safeguards[M]. Springer, 2010.

2 Santamarina A, Hfaiedh N, Letellier R, et al. Advanced neutronics tools for BWR design calculations[J]. Nuclear Engineering and Design, 2008, 238(8): 1965?1974.

3 Aldama D L, Leszczynski F, Trkov A. WIMS-D library update[R]. Final Report of a Co-ordinate Research Project,2003.

4 Macfarlane R E. A code system for producing pointwise and multigroup neutron and photon cross sections from ENDF/B evaluated nuclear data[R]. Los Alamos, New Mexico: Los Alamos National Laboratory, USA Report PSR-171, 1987.

5 Dong Z, Wu J, Ma X, et al. Development and verification of a 281-group WIMS-D library based on ENDF/B-VII.1[J]. Annals of Nuclear Energy, 2016, 91:189?194.

6 Marleau G, Hébert A, Roy R. A user guide for DRAGON 3.06[R]. Report IGE-174 Rev, 2008: 7.

7 Trkov A, Aldama D L. Definition of pseudo fission product data for reactor calculations[J]. Nuclear Energy in Central Europe, 1999, 99: 57?62.

8 楊偉焱, 畢光文, 楊波, 等. 低價(jià)值控制棒中子吸收體材料燃耗相關(guān)數(shù)據(jù)的制作及驗(yàn)證研究[J]. 核科學(xué)與工程, 2014, (3): 369?376.YANG Weiyan, BI Guangwen, YANG Bo, et al. Study on production and validation of low worth control rod material burnup calculation parameters[J]. Nuclear Science and Engineering, 2014, (3): 369?376.

9 胡也, 陳義學(xué), 楊壽海, 等. 基于Monte Carlo方法的壓水堆相關(guān)組件內(nèi)熱源計(jì)算與分析[J]. 核科學(xué)與工程,2013, 33(2): 168?174.HU Ye, CHEN Yixue, YANG Shouhai, et al. Calculation and analysis of heat sour of PWR assemblies based on Monte Carlo method[J]. Nuclear Science and Engineering,2013, 33(2): 168?174.

10 Trkov A. Processing of the fission product yield and decay data[J]. Nuclear Energy in Central Europe, 1999,99: 57?62.

11 曹博, 陳義學(xué), 李曉靜, 等. AP1000典型燃料組件物理特性初步分析[J]. 原子能科學(xué)技術(shù), 2013, 47(增 2):599?602. DOI: 10.7538/yzk.2013.47.S1.0599.CAO Bo, CHEN Yixue, LI Xiaojing, et al. Physical characteristics preliminary analysis on AP1000 fuel assembly[J]. Atomic Energy Science and Technology,2013, 47(Suppl 2): 599?602. DOI: 10.7538/yzk.2013.47.S1.0599.

12 Brown F B. MCNP - a general Monte Carlo N - particle transport code[M]. Version 5. Oak Ridge, TN: Los Alamos National Laboratory, 2003.

13 謝仲生, 尹邦華, 潘國品. 核反應(yīng)堆物理分析[M]. 北京: 原子能出版社, 1996.XIE Zhongsheng, YIN Banghua, PAN Guopin. The analysis of nuclear reactor physics[M]. Beijing: Atomic Energy Science and Technology, 1996.

14 Leszczynski F. NRSC system, user’s manual[R]. Internal Report CAB-CNEA, 1999.

15 Bucholz J A. SCALE: a modular code system for performing standardized computer analyses for licensing evaluation[R]. TN, USA: Oak Ridge National Laboratory,1982.

16 馬續(xù)波, 魯凡, 陳義學(xué), 等. 基于 cosRMC的新一代PWR柵元和組件燃耗計(jì)算[J]. 原子能科學(xué)技術(shù), 2016,50(11): 2010?2017. DOI: 10.7538/yzk.2016.50.11.2010.MA Xubo, LU Fan, CHEN Yixue, et al. Burnup calculation of next generation PWR cell and assembly based on cosRMC[J]. Atomic Energy Science and Technology, 2016, 50(11): 2010?2017. DOI: 10.7538/yzk.2016.50.11.2010.

Burnup information updating and verification based on WIMS-D library

XIAO Xiang1WU Jun1GONG You1YANG Tongrui2DONG Zhengyun3CHEN Yixue1,2
1(School of Nuclear Science and Engineering, North China Electric Power University, Beijing 102206, China)
2(National Nuclear Power Technology Co., Ltd, Beijing Software Technology Center, Beijing 100029, China)
3(Institute of Nuclear Safety Technology, Chinese Academy of Sciences, Hefei 230000, China)

Background: With the development of AP1000, the burnup has been raised constantly and the average burnup is expected to reach 60 GW d·t?1. However, the traditional WIMS-D library which is based on WIMS69 and XMAS172 energy group structure, has too many flaws such as the resonance overlap effect, lack of some nuclides,unreliable value of the fission product yield and thus can’t explicitly treat the nuclide of152Gd,160Gd,159Tb and160Tb in the pseudo fission product. Purpose: A new version of WIMS-D library based on ENDF/B-VII.1 with SHEM281 group structure is developed to improve the accuracy of cell calculation using the NJOY program. Methods: The new library is verified through the analysis of criticality and burnup benchmarks by using DRAGON program. Results:The results show that the calculation values are in good agreement with the benchmark. Conclusion: The new library has high precision and the result is reliable, which can be used for the PWR’s calculation.

XIAO Xiang, male, born in 1993, graduated from Hubei University of Science and Technology in 2015, master student, focusing on the development of nuclear data library and cross section processing

WU Jun, E-mail: wujun2007@ncepu.edu.cn

date: 2017-03-28, accepted date: 2017-05-12

SHEM281, Pseudo fission product, Fission product chain, Fission product yield

TL99

10.11889/j.0253-3219.2017.hjs.40.110602

國家自然科學(xué)基金(No.11275070)、國家科技重大專項(xiàng)(No.2011ZX06004-024-06-03-00)、中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)(No.2014MS51)資助第一作者：肖向，男，1993年出生，2015年畢業(yè)于湖北科技學(xué)院，現(xiàn)為碩士研究生，研究領(lǐng)域?yàn)楹藬?shù)據(jù)庫的開發(fā)與截面處理

吳軍，E-mail: wujun2007@ncepu.edu.cn

2017-03-28，

2017-05-12

Supported by National Natural Science Foundation of China (No.11275070), National Research and Development Projects

(No.2011ZX06004-024-06-03-00), Fundamental Research Funds for the Central Universities (No.2014MS51)