基于貝葉斯更新方法的235U熱中子獨立裂變產額協方差估計

王毅箴，崔夢蕾，郭 炯，鄔穎杰，劉保坤，李 富

(清華大學 核能與新能源技術研究院 先進核能技術協同創新中心先進反應堆工程與安全教育部重點實驗室，北京100084)

裂變產額描述了裂變產物在裂變系統發生裂變反應過程中產生的份額。根據裂變反應的進程不同，裂變產額可分為獨立裂變產額和累積裂變產額。獨立裂變產額描述了瞬發裂變中子釋放后，在裂變產物發生衰變前各個產物產生的份額；累積裂變產額則描述了裂變產物發生長時間衰變后裂變產物的份額[1]。由于裂變過程的時間尺度小、復雜性高，所以裂變產額的測量與模擬具有較大的不確定度。目前，ENDF/B-VII.1數據庫中的235U熱中子裂變產額數據來自于1994年England等的評價結果[2]，該結果僅提供裂變產額各自的方差信息，不包含裂變產額間的協方差信息，導致采用基于抽樣統計不確定度分析方法研究裂變產額的不確定度傳遞時，無法給出合理自洽的裂變產額概率分布[3]。因此，裂變產額間的協方差估計與評價工作，十分必要[4-5]。

貝葉斯更新方法廣泛應用于數據同化、數據調整及模型擬合等領域。在對裂變產額協方差的研究中， Kawano等首先基于鏈裂變產額估計了239Pu裂變產額的協方差[6]。其次， Pigni等將Kawano等的方法應用于ENDF/B-VII.1數據庫中，研究了基于累積裂變產額的獨立裂變產額調整方法，但未重點研究更新后的獨立裂變產額的協方差矩陣特征[7]。隨后， Fiorito等研究評價了分別用鏈裂變產額和累積裂變產額估計獨立裂變產額協方差的差異及其對傳遞的裂變產額不確定度的影響[8-9]。……