吳海成

(中國原子能科學(xué)研究院 核數(shù)據(jù)重點實驗室，北京 102413)

核數(shù)據(jù)協(xié)方差是評估核裝置設(shè)計不確定度的重要輸入，直接影響核裝置設(shè)計積分參數(shù)不確定度評估的質(zhì)量。核數(shù)據(jù)協(xié)方差也是核數(shù)據(jù)調(diào)整的重要輸入，直接影響核數(shù)據(jù)調(diào)整的結(jié)果[1-2]。雖然很難精確評價核數(shù)據(jù)的協(xié)方差，但仍需給出對不確定度的合理估計，以提高設(shè)計不確定度評估結(jié)果和核數(shù)據(jù)調(diào)整結(jié)果的可靠性。

核數(shù)據(jù)重點實驗室開展了大量核數(shù)據(jù)協(xié)方差評價方法研究工作[3-4]，掌握了基于線性最小二乘法(LLS)的模型無關(guān)協(xié)方差評價方法和基于誤差傳遞的模型依賴協(xié)方差評價方法，建立了基于不確定度源項分析(ASEU)的非模型依賴協(xié)方差評價方法，完成了235U核反應(yīng)截面的協(xié)方差數(shù)據(jù)的初步評價[5]，并與即將發(fā)布的現(xiàn)有的235U評價數(shù)據(jù)合并，形成了235U cov版評價數(shù)據(jù)。新的協(xié)方差評價數(shù)據(jù)是否合理有必要通過適當(dāng)?shù)臏y試進(jìn)行評估。

與核數(shù)據(jù)本身的準(zhǔn)確度可通過宏觀檢驗來進(jìn)行驗證不同，核數(shù)據(jù)協(xié)方差是否合理不能直接通過積分實驗進(jìn)行評估。為評估上述協(xié)方差評價結(jié)果的合理性，本工作利用經(jīng)濟(jì)合作組織核能局(OECD/NEA)國際核數(shù)據(jù)評價合作工作組第33研究小組(WPEC/SG33)提出的核數(shù)據(jù)調(diào)整例題[6]對235U cov版235U協(xié)方差評價數(shù)據(jù)進(jìn)行測試。通過與JENDL-4.0庫[7]調(diào)整結(jié)果的比對分析，對新評價的協(xié)方差評價數(shù)據(jù)的合理性進(jìn)行評估。

1 235U協(xié)方差數(shù)據(jù)測試

1.1 測試方法

本工作以WPEC/SG33提出的包含20個積分參數(shù)的核數(shù)據(jù)調(diào)整基準(zhǔn)例題為基礎(chǔ)，補(bǔ)充了國際核臨界安全手冊ICSBEP2006[8]中臨界基準(zhǔn)實驗HMF001(Godiva模型)的有效增殖因數(shù)keff作為第21個積分量，開展了基于JENDL-4.0庫(J40組)以及在JENDL-4.0庫基礎(chǔ)上更新235U截面和協(xié)方差數(shù)據(jù)(C組)的核數(shù)據(jù)調(diào)整比對研究。通過對比兩組數(shù)據(jù)驗前和驗后積分量及其不確定度、截面或裂變譜及其不確定度的變化來分析本次測試的235U協(xié)方差數(shù)據(jù)是否合理。

1.2 基準(zhǔn)例題

SG33基準(zhǔn)例題是針對快堆的核數(shù)據(jù)調(diào)整基準(zhǔn)例題，考慮了7個快譜臨界基準(zhǔn)實驗，涉及11個核素、8個主要核反應(yīng)和20個積分參數(shù)，具體信息列于表1?；鶞?zhǔn)實驗的定義詳見ICSBEP2006，所選積分參數(shù)既包含keff，也包含反應(yīng)率比?？紤]到被測對象是235U數(shù)據(jù)，額外增加Godiva實驗[8]作為補(bǔ)充，積分量增加該實驗的keff。Godiva實驗是ICSBEP手冊中HEU-MET-FAST-001實驗的例題高濃鈾裸球簡化模型，其對235U核反應(yīng)數(shù)據(jù)靈敏。

表1 SG33核數(shù)據(jù)調(diào)整基準(zhǔn)例題輸入Table 1 Input for SG33 nuclear data adjustment benchmark

1.3 測試準(zhǔn)備

1) 多群截面及協(xié)方差數(shù)據(jù)

測試使用的截面及協(xié)方差數(shù)據(jù)以WPEC/SG33使用的33群JENDL-4.0截面和協(xié)方差數(shù)據(jù)為基礎(chǔ)。J40組使用的數(shù)據(jù)與WPEC/SG33發(fā)布的數(shù)據(jù)完全一致。C組的截面和協(xié)方差數(shù)據(jù)是在J40組的基礎(chǔ)上以235U cov版235U數(shù)據(jù)更新截面和多群協(xié)方差數(shù)據(jù)得到的。235U cov版235U評價數(shù)據(jù)雖包含完整的核反應(yīng)截面、角分布、能譜和雙微分截面文檔，但不包含完整的協(xié)方差數(shù)據(jù)文檔，僅在10 keV以上的平滑區(qū)給出了截面的協(xié)方差文檔。其中，裂變截面和俘獲截面的協(xié)方差是利用LLS法基于實驗數(shù)據(jù)評價得到的，非彈性散射反應(yīng)截面的協(xié)方差是通過模型參數(shù)誤差傳遞得到的。235U cov版的評價數(shù)據(jù)既不包含熱能區(qū)和共振區(qū)的協(xié)方差數(shù)據(jù)，也不包含裂變平均中子數(shù)(nubar)、裂變譜(chi)以及散射角分布的不確定度。因此，C組中僅利用235U cov版數(shù)據(jù)對235U的截面以及非彈性散射(n,inl)、裂變、俘獲反應(yīng)10 keV以上能區(qū)的協(xié)方差數(shù)據(jù)進(jìn)行了更新。nubar以及多群裂變譜的協(xié)方差數(shù)據(jù)則采用ENDF/B-Ⅷ.0庫[9]數(shù)據(jù)制作。除了更新截面自身的協(xié)方差，裂變和俘獲截面之間的協(xié)方差也進(jìn)行了更新。

新制作的多群截面和協(xié)方差數(shù)據(jù)采用NJOY99程序[10]加工得到。制作中能群結(jié)構(gòu)采用了WPEC/SG33推薦的33群，權(quán)重譜選用了快堆譜，溫度為常溫(293.6 K)，本底截面為10-18m2。

2) 靈敏度數(shù)據(jù)

考慮到國際主要評價數(shù)據(jù)庫計算的靈敏度差異較小，表1中20個積分量的靈敏度數(shù)據(jù)直接使用了WPEC/SG33中以JENDL-4.0庫為基礎(chǔ)計算的靈敏度數(shù)據(jù)。對于額外增加的Godiva實驗，keff及其對核反應(yīng)截面的靈敏度采用SCALE-5.1程序[11-12]基于CENDL-3.1庫[13]計算，主要反應(yīng)道的靈敏度計算結(jié)果如圖1所示。

3) 核數(shù)據(jù)調(diào)整程序

本次測試的核數(shù)據(jù)調(diào)整程序采用自主研發(fā)的NDAC程序[14]。該程序是基于LLS法的核數(shù)據(jù)調(diào)整程序，輸入量包括驗前的多群截面、多群截面不確定度及相關(guān)系數(shù)矩陣、積分量實驗值及不確定度、積分量計算值及不確定度、積分量相關(guān)系數(shù)矩陣以及積分量靈敏度，輸出量包括驗后多群截面、多群截面不確定度及相關(guān)系數(shù)矩陣、驗后積分量計算值及不確定度、驗后積分量相關(guān)系數(shù)矩陣等。上述驗后量的計算公式詳見WPEC/SG33報告[6]。

圖1 Godiva實驗keff對235U截面和裂變譜的靈敏度Fig.1 Sensitivity of keff to 235U reaction cross section and fission spectrum for Godiva experiment

2 測試結(jié)果

2.1 積分量及不確定度調(diào)整結(jié)果

圖2、3示出了調(diào)整前后積分量及不確定度的比較。除圖2中Godiva的keff計算結(jié)果及驗前不確定度存在約300 pcm的差異，其他驗前和驗后積分量及不確定結(jié)果基本一致，表明新評價的235U協(xié)方差數(shù)據(jù)引起的Godiva例題keff計算值不確定度與JENDL-4.0庫的基本一致。

圖2 調(diào)整前后keff及其不確定度的比較Fig.2 Comparison of priori and posteriori keff and uncertainties

圖3 調(diào)整前后反應(yīng)率比及其不確定度的比較Fig.3 Comparison of priori and posteriori reaction rate ratios and uncertainties

2.2 微分量及協(xié)方差調(diào)整結(jié)果

1) 彈性散射截面

圖4示出了調(diào)整前后235U(n,el)截面及不確定度的比較結(jié)果。測試結(jié)果顯示出兩組數(shù)據(jù)的彈性散射截面調(diào)整量(中間藍(lán)色實線和虛線)存在明顯差異，但調(diào)整后的不確定度(綠色和紅色虛線)基本一致。這主要是因為雖然兩組數(shù)據(jù)中235U(n,el)截面不同但協(xié)方差數(shù)據(jù)相同。

圖4 235U(n,el)截面調(diào)整趨勢及驗前與驗后相對不確定度的比較Fig.4 Comparison of adjustment trend of 235U(n,el) cross section and priori and posteriori uncertainties

2) 非彈性散射截面

測試顯示235U cov版非彈性散射反應(yīng)截面的協(xié)方差評價不合理。圖5所示的測試結(jié)果顯示出兩組數(shù)據(jù)的非彈性散射截面調(diào)整量和不確定度均存在極大差異。J40非彈散射截面幾乎無調(diào)整，而C組的后驗截面在9～15 keV附近被調(diào)低了15.4倍之多，其原因在于235U cov版在接近2～20 keV能區(qū)出現(xiàn)高達(dá)231的相對不確定度(遠(yuǎn)大于ENDF/B-Ⅶ.1[15]的評價值)，在調(diào)整中給出了過高的權(quán)重，如圖6所示。

圖5 235U(n,inl)截面調(diào)整趨勢及驗前與驗后相對不確定度的比較Fig.5 Comparison of adjustment trend of 235U(n,inl) cross section and priori and posteriori uncertainties

圖6 ENDF/B-Ⅶ.1和235U cov版的(n,inl)反應(yīng)截面的相對不確定度比較Fig.6 Comparison of relative uncertainty of 235U(n,inl) cross section from ENDF/B-Ⅶ.1 and 235U cov

3) 裂變截面

235U cov版裂變截面協(xié)方差也不合理，且并非文獻(xiàn)[5]中所述來自ENDF/B-Ⅷ.0庫。裂變截面及不確定度調(diào)整前后的比較示于圖7。結(jié)合圖8中ENDF/B-Ⅷ.0、ENDF/B-Ⅶ.1、JENDL-4.0和235U cov裂變截面不確定度的比較結(jié)果，在圖7中9.1 keV～20 MeV能區(qū)，235U cov版評價數(shù)據(jù)給出的驗前和驗后235U(n,f)多群截面的不確定度非常小，整體上小于0.5%，甚至在9～41 keV能區(qū)范圍內(nèi)小于0.1%；而JENDL-4.0庫的驗前和驗后不確定度均接近0.7%。從調(diào)整前后不確定變化來看，即使測試中已包含了對MeV能區(qū)235U(n,f)截面靈敏的Godiva實驗，測試也不支持9.1 keV～20 MeV能區(qū)的相對不確定度從0.7%下降到0.1%的水平。因此，235U cov版的(n,f)截面協(xié)方差存在顯著低估的問題。

圖7 235U(n,f) 截面調(diào)整趨勢及驗前與驗后相對不確定度的比較Fig.7 Comparison of adjustment trend of 235U(n,f) cross section and priori and posteriori uncertainties

圖8 ENDF/B、JENDL-4.0庫與235U cov的(n,f)反應(yīng)截面的相對不確定度比較Fig.8 Comparison of relative uncertaintiesfor (n,f) cross section from ENDF/B andJENDL-4.0 library and 235U cov

4) 俘獲截面

對于235U cov中基于ASEU方法評價的俘獲截面協(xié)方差，測試結(jié)果也顯示出不確定度偏小的問題。俘獲截面及不確定度調(diào)整前后的比較示于圖9。在800 keV以上能區(qū)，235U cov數(shù)據(jù)的驗前與驗后俘獲截面及其不確定度沒有明顯變化，其中該能區(qū)俘獲截面的最大不確定度小于20%；而J40組在該能區(qū)的驗后最大不確定度約70%且均高于20%。這表明當(dāng)前選擇的基準(zhǔn)實驗集，特別是Godiva，并不能支持235U俘獲截面在該能區(qū)獲得小于20%水平的不確定度約束。235U cov版的俘獲截面協(xié)方差評價值顯然偏低。

圖9 235U(n,γ)截面調(diào)整趨勢及驗前與驗后相對不確定度的比較Fig.9 Comparison of adjustment trend of 235U(n,γ) cross section and priori and posteriori uncertainties

3 總結(jié)

本工作利用開展核數(shù)據(jù)調(diào)整測試的方法，對235U cov版235U裂變、俘獲截面及非彈性散射截面的協(xié)方差數(shù)據(jù)進(jìn)行了測試分析，評估了協(xié)方差數(shù)據(jù)的合理性。測試分析結(jié)果顯示，235U cov版235U數(shù)據(jù)的非彈性散射截面協(xié)方差在閾能附近存在估值過大的問題；裂變和俘獲截面的不確定度在10 keV以上能區(qū)則整體偏小。俘獲截面的不確定度甚至小于基于JENDL-4.0庫的驗后截面不確定度。