王毅箴,崔夢(mèng)蕾,郭 炯,鄔穎杰,劉保坤,李 富
(清華大學(xué) 核能與新能源技術(shù)研究院 先進(jìn)核能技術(shù)協(xié)同創(chuàng)新中心先進(jìn)反應(yīng)堆工程與安全教育部重點(diǎn)實(shí)驗(yàn)室,北京100084)
裂變產(chǎn)額描述了裂變產(chǎn)物在裂變系統(tǒng)發(fā)生裂變反應(yīng)過(guò)程中產(chǎn)生的份額。根據(jù)裂變反應(yīng)的進(jìn)程不同,裂變產(chǎn)額可分為獨(dú)立裂變產(chǎn)額和累積裂變產(chǎn)額。獨(dú)立裂變產(chǎn)額描述了瞬發(fā)裂變中子釋放后,在裂變產(chǎn)物發(fā)生衰變前各個(gè)產(chǎn)物產(chǎn)生的份額;累積裂變產(chǎn)額則描述了裂變產(chǎn)物發(fā)生長(zhǎng)時(shí)間衰變后裂變產(chǎn)物的份額[1]。由于裂變過(guò)程的時(shí)間尺度小、復(fù)雜性高,所以裂變產(chǎn)額的測(cè)量與模擬具有較大的不確定度。目前,ENDF/B-VII.1數(shù)據(jù)庫(kù)中的235U熱中子裂變產(chǎn)額數(shù)據(jù)來(lái)自于1994年England等的評(píng)價(jià)結(jié)果[2],該結(jié)果僅提供裂變產(chǎn)額各自的方差信息,不包含裂變產(chǎn)額間的協(xié)方差信息,導(dǎo)致采用基于抽樣統(tǒng)計(jì)不確定度分析方法研究裂變產(chǎn)額的不確定度傳遞時(shí),無(wú)法給出合理自洽的裂變產(chǎn)額概率分布[3]。因此,裂變產(chǎn)額間的協(xié)方差估計(jì)與評(píng)價(jià)工作,十分必要[4-5]。
貝葉斯更新方法廣泛應(yīng)用于數(shù)據(jù)同化、數(shù)據(jù)調(diào)整及模型擬合等領(lǐng)域。在對(duì)裂變產(chǎn)額協(xié)方差的研究中, Kawano等首先基于鏈裂變產(chǎn)額估計(jì)了239Pu裂變產(chǎn)額的協(xié)方差[6]。其次, Pigni等將Kawano等的方法應(yīng)用于ENDF/B-VII.1數(shù)據(jù)庫(kù)中,研究了基于累積裂變產(chǎn)額的獨(dú)立裂變產(chǎn)額調(diào)整方法,但未重點(diǎn)研究更新后的獨(dú)立裂變產(chǎn)額的協(xié)方差矩陣特征[7]。隨后, Fiorito等研究評(píng)價(jià)了分別用鏈裂變產(chǎn)額和累積裂變產(chǎn)額估計(jì)獨(dú)立裂變產(chǎn)額協(xié)方差的差異及其對(duì)傳遞的裂變產(chǎn)額不確定度的影響[8-9]。本文基于貝葉斯更新方法,利用ENDF/B-VII.1數(shù)據(jù)庫(kù)中的累積裂變產(chǎn)額數(shù)據(jù)及獨(dú)立裂變產(chǎn)額自身的物理約束條件,對(duì)ENDF/B-VII.1數(shù)據(jù)庫(kù)中獨(dú)立裂變產(chǎn)額的協(xié)方差進(jìn)行估計(jì),并對(duì)估計(jì)結(jié)果進(jìn)行分析和討論。
貝葉斯更新方法(Bayesian updating method)是一種核數(shù)據(jù)調(diào)整方法,它根據(jù)核數(shù)據(jù)的新認(rèn)知(新評(píng)估結(jié)果或物理約束條件)對(duì)核數(shù)據(jù)進(jìn)行調(diào)整,從而保持核數(shù)據(jù)與新認(rèn)知的一致性[10-11]。該調(diào)整過(guò)程是在概率空間定義下進(jìn)行的,因此,核數(shù)據(jù)x先被賦予某種概率密度分布,稱為先驗(yàn)概率密度函數(shù)p(x)。核數(shù)據(jù)的新認(rèn)知信息y被同樣賦予某種概率密度分布,從而計(jì)算似然值p(y|x)。根據(jù)貝葉斯定理,核數(shù)據(jù)的后驗(yàn)概率分布p(x|y)可表示為
(1)
后驗(yàn)概率分布的最大值即為核數(shù)據(jù)的最佳調(diào)整值。由式(1)可知,貝葉斯更新方法對(duì)核數(shù)據(jù)進(jìn)行調(diào)整具有3個(gè)特點(diǎn):
1)允許按序分級(jí)地引入不同的核數(shù)據(jù)新認(rèn)知y1,y2,…,yN,這樣可以避免每當(dāng)有新認(rèn)知產(chǎn)生時(shí),對(duì)以往核數(shù)據(jù)認(rèn)知的重新引入與計(jì)算,也可以方便地分級(jí)觀察和研究不同新認(rèn)知對(duì)核數(shù)據(jù)調(diào)整值的影響。式(1)可改寫(xiě)為
(2)
2)貝葉斯更新給出的核數(shù)據(jù)最佳調(diào)整值依賴于先驗(yàn)概率密度分布和似然函數(shù),它們可以通過(guò)以往人們對(duì)核數(shù)據(jù)概率分布的研究結(jié)果或基于最大信息熵原理給出[10],從而保證核數(shù)據(jù)調(diào)整的合理性。
3)由于貝葉斯更新方法是基于核數(shù)據(jù)后驗(yàn)概率分布進(jìn)行調(diào)整的,因此,基于獲得的后驗(yàn)概率分布,可以對(duì)核數(shù)據(jù)的高階矩信息進(jìn)行模擬或直接計(jì)算,從而可以獲得核數(shù)據(jù)在新認(rèn)知約束下的協(xié)方差信息。
本文基于貝葉斯更新方法對(duì)獨(dú)立裂變產(chǎn)額進(jìn)行調(diào)整,并對(duì)ENDF/B-VII.1核數(shù)據(jù)庫(kù)中缺失的協(xié)方差信息進(jìn)行估計(jì)。由于ENDF/B-VII.1核數(shù)據(jù)庫(kù)中僅給出了獨(dú)立裂變產(chǎn)額的均值和標(biāo)準(zhǔn)差,由最大信息熵原理可知,采用正態(tài)分布描述獨(dú)立裂變產(chǎn)額能夠得到最大信息熵,從而可以充分利用已知的不確定度信息,最小地引入人為的主觀影響。同理,對(duì)累積裂變產(chǎn)額也同樣賦予正態(tài)分布計(jì)算似然值。
(3)
(4)
其中,Yi為獨(dú)立裂變產(chǎn)額向量,Yi∈RP×1;P為總裂變產(chǎn)物數(shù)目;Yc,0,Yi,0分別為ENDF/B-VII.1提供的累積裂變產(chǎn)額與獨(dú)立裂變產(chǎn)額的均值;Vc,Vi,0分別為累積裂變產(chǎn)額和獨(dú)立裂變產(chǎn)額的協(xié)方差矩陣,它們均為對(duì)角陣,對(duì)角元素為ENDF/B-VII.1數(shù)據(jù)庫(kù)中提供的各個(gè)裂變產(chǎn)額的方差,Vc∈RP×P,Vi,0∈RP×P;M為累積裂變產(chǎn)額與獨(dú)立裂變產(chǎn)額的線性映射矩陣。每個(gè)裂變產(chǎn)物核素可由其質(zhì)量數(shù)A、質(zhì)子數(shù)Z及其同核異能態(tài)I確定,記為(A,Z,I)。累積裂變產(chǎn)額與獨(dú)立裂變產(chǎn)額之間的關(guān)系可表示為[1]
Yc(A,Z,I)=Yi(A,Z,I)+
(5)
其中,Yc(A,Z,I),Yi(A,Z,I)分別表示裂變產(chǎn)物核素(A,Z,I)的累積裂變產(chǎn)額和獨(dú)立裂變產(chǎn)額;b(A′,Z′,I′→A,Z,I)為裂變產(chǎn)物核素(A′,Z′,I′)發(fā)生衰變到裂變產(chǎn)物(A,Z,I)的分支比。聯(lián)立所有裂變產(chǎn)物核素的累積裂變產(chǎn)額計(jì)算式,根據(jù)式(5)可得
Yc=Yi+BYc
(6)
Yc=(E-B)-1Yi
(7)
Yc=MYi
(8)
其中,Yc為累積裂變產(chǎn)額向量,Yc∈RP×1;B為裂變產(chǎn)物核素之間的分支比矩陣,B∈RP×P;E為單位矩陣,E∈RP×P。
由式(1)可以得到累積裂變產(chǎn)額更新后的獨(dú)立裂變產(chǎn)額的后驗(yàn)概率分布為
(9)
經(jīng)推導(dǎo)可得,獨(dú)立裂變產(chǎn)額的后驗(yàn)概率分布也為正態(tài)分布,只是分布參數(shù)發(fā)生了變化。因此,獨(dú)立裂變產(chǎn)額的后驗(yàn)概率分布可表示為
(10)
其中,
(11)
Yi,1=Yi,0+Vi,1MTV-1(Yc,0-MYi,0)
(12)
由于正態(tài)分布均值點(diǎn)的鄰域內(nèi)具有最大概率,所以,獨(dú)立裂變產(chǎn)額的最佳調(diào)整值即為后驗(yàn)概率分布的均值Yi,1。后驗(yàn)概率分布的協(xié)方差矩陣Vi,1,可用于研究獨(dú)立裂變產(chǎn)額不確定度的傳遞。
獨(dú)立裂變產(chǎn)額應(yīng)滿足3個(gè)物理約束條件:
1)在忽略三元裂變的情況下,總獨(dú)立裂變產(chǎn)額應(yīng)為2.0;
2)在忽略三元裂變的情況下,各個(gè)裂變產(chǎn)物的質(zhì)量數(shù)以其各自的獨(dú)立裂變產(chǎn)額為權(quán)重的平均值,應(yīng)為裂變系統(tǒng)裂變前復(fù)合核的總質(zhì)量減去釋放的平均瞬發(fā)中子數(shù);
3)各個(gè)裂變產(chǎn)物的電荷數(shù)以其各自的獨(dú)立裂變產(chǎn)額為權(quán)重的平均值,應(yīng)與裂變系統(tǒng)的總電荷數(shù)相等。
由于本文依據(jù)的ENDF/B-VII.1數(shù)據(jù)庫(kù)中未提供H,He, Li, Be, B, C, N, Ne, Ti等輕核素的相關(guān)數(shù)據(jù)[12],所以,本文分析中均忽略三元裂變。
總裂變產(chǎn)額維持二裂變的約束條件,采用正態(tài)分布計(jì)算似然值,則
p(Yi|Ytot,Yc,0)∝p(Ytot|Yi,Yc,0)p(Yi|Yc,0)∝
(13)
經(jīng)推導(dǎo),引入二裂變守恒后的獨(dú)立裂變產(chǎn)額后驗(yàn)概率分布p(Yi|Ytot,Yc,0),可化為正態(tài)分布形式,即
(14)
其中,
(15)
(16)
其中,Ytot為總獨(dú)立裂變產(chǎn)額;σ2為總裂變產(chǎn)額的求和精度,本文取σ2=10-5;U為單位向量,U∈RP×1。
同理,引入裂變系統(tǒng)總質(zhì)量數(shù)守恒約束條件,則
p(Yi│Atot,Ytot,Yc,0)∝p(Atot│Yi,Ytot,Yc,0)·
(17)
p(Yi│Atot,Ytot,Yc,0)∝
(18)
(19)
(20)
其中,Atot為裂變系統(tǒng)總質(zhì)量數(shù),該值可通過(guò)各個(gè)裂變核素的質(zhì)量數(shù)以其ENDF/B-VII.1數(shù)據(jù)庫(kù)中提供的獨(dú)立裂變產(chǎn)額加權(quán)平均給出;A為各個(gè)裂變核素的質(zhì)量數(shù)向量,A∈RP×1。
引入裂變系統(tǒng)的總電荷數(shù)守恒約束條件,則
p(Yi│Ztot,Atot,Ytot,Yc,0)∝p(Ztot│Yi,Atot,Ytot,Yc,0)·
(21)
p(Yi│Ztot,Atot,Ytot,Yc,0)∝
(22)
(23)
(24)
其中,Ztot為裂變系統(tǒng)總電荷數(shù),該值可通過(guò)各個(gè)裂變核素的電荷數(shù)以其ENDF/B-VII.1數(shù)據(jù)庫(kù)中提供的獨(dú)立裂變裂變產(chǎn)額加權(quán)平均給出;Z為各個(gè)裂變核素的電荷數(shù)向量,Z∈RP×1。
最終得到的獨(dú)立裂變產(chǎn)額的后驗(yàn)概率分布p(Yi│Ztot,Atot,Ytot,Yc,0),能夠給出獨(dú)立裂變產(chǎn)額在累積裂變產(chǎn)額更新及自身的物理約束下的最佳調(diào)整值Yi,4及合理估計(jì)的協(xié)方差矩陣Vi,4。利用Yi,4和Vi,4進(jìn)行獨(dú)立裂變產(chǎn)額抽樣,能夠獲得符合獨(dú)立裂變物理約束條件且累積裂變產(chǎn)額一致的獨(dú)立裂變產(chǎn)額樣本數(shù)據(jù),為后續(xù)開(kāi)展不確定度傳遞研究提供基礎(chǔ)。
獨(dú)立裂變產(chǎn)額的調(diào)整結(jié)果,如圖1所示。其中,藍(lán)色點(diǎn)代表ENDF/B-VII.1數(shù)據(jù)庫(kù)中提供的235U熱中子裂變的獨(dú)立裂變產(chǎn)額,即先驗(yàn)獨(dú)立裂變產(chǎn)額;橙色點(diǎn)代表經(jīng)過(guò)貝葉斯更新后的獨(dú)立裂變產(chǎn)額最佳調(diào)整值Yi,4,即后驗(yàn)獨(dú)立裂變產(chǎn)額的均值。為方便比較分析,縱軸給出的是獨(dú)立裂變質(zhì)量產(chǎn)額Yi,m,即按質(zhì)量數(shù)對(duì)裂變核素分組,并對(duì)每個(gè)質(zhì)量組內(nèi)核素的獨(dú)立裂變產(chǎn)額求和的結(jié)果。由圖1可見(jiàn),獨(dú)立裂變質(zhì)量產(chǎn)額調(diào)整較大的核素在曲線的兩峰處,對(duì)應(yīng)質(zhì)量數(shù)分別為80~110和130~150。這是由于兩峰處裂變產(chǎn)物累積裂變產(chǎn)額的測(cè)量不確定度較小,累積裂變產(chǎn)額較其他位置核素的累積裂變產(chǎn)額更為準(zhǔn)確,通過(guò)貝葉斯更新方法對(duì)兩峰處裂變產(chǎn)物的獨(dú)立裂變產(chǎn)額調(diào)整時(shí),調(diào)整效果較為明顯。
圖1 獨(dú)立裂變產(chǎn)額的調(diào)整結(jié)果Fig.1 Independent fission yields adjustment results
為了分析貝葉斯更新調(diào)整后的獨(dú)立裂變產(chǎn)額Yi,4與ENDF/B-VII.1數(shù)據(jù)庫(kù)提供的累積裂變產(chǎn)額評(píng)價(jià)值Yc,0的一致性,根據(jù)式(8),分別利用先驗(yàn)獨(dú)立裂變產(chǎn)額Yi,0與后驗(yàn)獨(dú)立裂變產(chǎn)額Yi,4計(jì)算對(duì)應(yīng)的累積裂變產(chǎn)額,并記為先驗(yàn)累積裂變產(chǎn)額計(jì)算值與后驗(yàn)累積裂變產(chǎn)額計(jì)算值。為方便分析,將累積裂變產(chǎn)額按照質(zhì)量數(shù)進(jìn)行分組求和,得到先驗(yàn)累積裂變質(zhì)量產(chǎn)額與后驗(yàn)累積裂變質(zhì)量產(chǎn)額計(jì)算值。然后分別計(jì)算它們與ENDF/B-VII.1數(shù)據(jù)庫(kù)中提供的累積裂變產(chǎn)額評(píng)價(jià)值之比r(Yc,m),得到r(Yc,m)隨質(zhì)量數(shù)A的分布曲線,如圖2所示。
由圖2可見(jiàn),用后驗(yàn)獨(dú)立裂變產(chǎn)額計(jì)算的累積裂變產(chǎn)額的估計(jì)不確定度在20%以內(nèi),與用先驗(yàn)獨(dú)立裂變產(chǎn)額計(jì)算的累積裂變產(chǎn)額的估計(jì)不確定度相比減少約30%,提高了ENDF/B-VII.1獨(dú)立裂變產(chǎn)額與累積裂變產(chǎn)額的一致性。對(duì)質(zhì)量數(shù)為110~130的核素,更新后的估計(jì)不確定度有所放大,這是由于這些核素的累積裂變產(chǎn)額的先驗(yàn)計(jì)算結(jié)果與ENDF/B-VII.1數(shù)據(jù)庫(kù)提供的評(píng)價(jià)值之間的偏差較大造成的。為表征先驗(yàn)累積裂變產(chǎn)額計(jì)算值與ENDF/B-VII.1數(shù)據(jù)庫(kù)提供的累積裂變產(chǎn)額評(píng)價(jià)值的一致性,定義一致性因子H為
(25)
其中,NA為質(zhì)量鏈A中裂變產(chǎn)物核素的數(shù)目;Yc,n為質(zhì)量鏈A中裂變產(chǎn)物n先驗(yàn)累積裂變產(chǎn)額計(jì)算值;Yc,0,n為質(zhì)量鏈A中ENDF/B-VII.1提供的裂變產(chǎn)物n累積裂變產(chǎn)額評(píng)價(jià)值;σc,0,n為ENDF/B-VII.1提供的裂變產(chǎn)物n累積裂變產(chǎn)額的標(biāo)準(zhǔn)偏差[13]。H越大,表示計(jì)算值與評(píng)價(jià)值之間的一致性越弱。圖3對(duì)比了先驗(yàn)累積裂變質(zhì)量產(chǎn)額計(jì)算值與ENDF/B-VII.1數(shù)據(jù)庫(kù)中評(píng)價(jià)值的一致性。由圖3可見(jiàn),質(zhì)量數(shù)為110~130和80~90核素的一致性因子H較其周圍核素的H值大,具有較弱的一致性,這種較弱的一致性影響了貝葉斯更新的獨(dú)立裂變產(chǎn)額調(diào)整結(jié)果。由于這些核素的獨(dú)立裂變產(chǎn)額較小(位于產(chǎn)額分布曲線上兩峰之間的低谷處),因此,這些核素的獨(dú)立裂變產(chǎn)額可能需要在未來(lái)的評(píng)價(jià)數(shù)據(jù)庫(kù)中重新評(píng)估。
圖3 先驗(yàn)累積裂變質(zhì)量產(chǎn)額計(jì)算值與ENDF/B-VII.1數(shù)據(jù)庫(kù)中評(píng)價(jià)值的一致性對(duì)比Fig.3 Prior calculated cumulative fission yields consistencycomparison with ENDF/B-VII.1 evaluations
圖4給出了后驗(yàn)與先驗(yàn)獨(dú)立裂變產(chǎn)額標(biāo)準(zhǔn)偏差的比值r(σ)隨質(zhì)量數(shù)A及電荷數(shù)Z的分布。由圖4可見(jiàn),多數(shù)裂變產(chǎn)物獨(dú)立裂變產(chǎn)額的不確定度經(jīng)過(guò)更新后有所減小,約為更新前的70%,而且位于兩峰處的核素,獨(dú)立裂變產(chǎn)額不確定度減小的幅度較大,約為更新前的20%。這是由于隨著累積裂變產(chǎn)額評(píng)價(jià)值及其引入的不確定度的減小,且經(jīng)過(guò)獨(dú)立裂變產(chǎn)額物理約束條件的限制和似然函數(shù)的加權(quán)調(diào)整,獨(dú)立裂變產(chǎn)額的后驗(yàn)概率分布會(huì)在最佳調(diào)整值附近聚集更多,從而使后驗(yàn)概率分布函數(shù)更加尖銳,標(biāo)準(zhǔn)偏差更小。
圖4 后驗(yàn)與先驗(yàn)獨(dú)立裂變產(chǎn)額標(biāo)準(zhǔn)偏差的比值r(σ)隨質(zhì)量數(shù)A及電荷數(shù)Z的分布Fig.4 Independent fission yields posterior to priorstandard deviation ratios r(σ) distribution againstmass number A and charge number Z
隨著累積裂變產(chǎn)額測(cè)量值及獨(dú)立裂變產(chǎn)額物理約束條件的引入,獨(dú)立裂變產(chǎn)額后驗(yàn)協(xié)方差矩陣會(huì)不斷發(fā)生變化。圖5給出了不同更新步驟下獨(dú)立裂變產(chǎn)額相關(guān)系數(shù)矩陣的變化過(guò)程。本文對(duì)裂變產(chǎn)物按照質(zhì)量數(shù)大小進(jìn)行分組并降序排列,在每個(gè)質(zhì)量數(shù)分組內(nèi)按照裂變產(chǎn)物的電荷數(shù)大小升序排列。
由圖5可見(jiàn),相關(guān)系數(shù)呈團(tuán)簇狀出現(xiàn)在主對(duì)角線附近,這是由于累積裂變產(chǎn)額和獨(dú)立裂變產(chǎn)額之間是通過(guò)裂變產(chǎn)物的放射性衰變相聯(lián)系的,相關(guān)性主要出現(xiàn)在有放射性衰變聯(lián)系的核素之間。隨著二裂變約束的引入,可以看到圖中Vi,2在相距較遠(yuǎn)的核素位置處產(chǎn)生了相關(guān)性。而隨著裂變系統(tǒng)的質(zhì)量數(shù)和電荷數(shù)守恒的引入,更多的相關(guān)性被引入到獨(dú)立裂變產(chǎn)額的協(xié)方差矩陣中。
(a)Vi,1:p(Yi|Yc,0)
(b)Vi,2:p(Yi|Ytot,Yc,0)
(c)Vi,3:p(Yi|Atot,Ytot,Yc,0)
(d)Vi,4:p(Yi|Ztot,Atot,Ytot,Yc,0)
圖6為圖5中的Vi,4的局部放大圖。其中,紅色點(diǎn)代表正相關(guān)性系數(shù),黑色點(diǎn)代表負(fù)相關(guān)性系數(shù)。由圖6可見(jiàn),在主對(duì)角線附近呈團(tuán)簇狀出現(xiàn)的多為負(fù)相關(guān)性系數(shù),即由裂變產(chǎn)物衰變引起的核素間的相關(guān)性多為負(fù)相關(guān);而在主對(duì)角線以外由獨(dú)立裂變產(chǎn)額的物理約束引起的相關(guān)性則多為正相關(guān)。圖7為圖6左側(cè)矩陣的下三角元素的相關(guān)系數(shù)散點(diǎn)圖。由圖7可見(jiàn),貝葉斯更新方法得到的獨(dú)立裂變產(chǎn)額之間的相關(guān)性系數(shù)多為負(fù)相關(guān),且多數(shù)的相關(guān)性較弱,相關(guān)系數(shù)小于0.25。少數(shù)核素之間存在較強(qiáng)的負(fù)相關(guān)性,極少數(shù)核素之間存在較強(qiáng)的正相關(guān)性。
圖6 相關(guān)系數(shù)矩陣Vi,4的局部放大圖Fig.6 Correlation matrix of independentfission yields Vi,4 with local zoom-in view
圖7 相關(guān)系數(shù)矩陣Vi,4 下三角元素散點(diǎn)圖Fig.7 Scatter plot of lower triangle elements in correlation matrix Vi,4
近年來(lái),國(guó)外雖然已經(jīng)開(kāi)展了很多基于實(shí)驗(yàn)鏈產(chǎn)額和累積裂變產(chǎn)額數(shù)據(jù)更新獨(dú)立裂變產(chǎn)額和協(xié)方差數(shù)據(jù)的研究工作,但是,目前ENDF/B-VII.1數(shù)據(jù)庫(kù)還未提供獨(dú)立裂變產(chǎn)額的協(xié)方差數(shù)據(jù)。文獻(xiàn)[14]中利用本文研究的獨(dú)立裂變產(chǎn)額的協(xié)方差數(shù)據(jù),研究了球床式高溫氣冷堆裂變產(chǎn)額的不確定度傳遞及獨(dú)立裂變產(chǎn)額不確定度對(duì)堆芯平衡態(tài)下有效增殖因子和核素濃度不確定度的貢獻(xiàn),結(jié)果表明,有效增殖因子的相對(duì)不確定度約為2.44×10-4(燃耗為90 GW·d·t-1),該結(jié)果和REBUS基準(zhǔn)題中單一燃料棒柵元的相應(yīng)評(píng)價(jià)結(jié)果6.0×10-4(燃耗為54 GW·d·t-1)較為接近[9],說(shuō)明本文的評(píng)估工作是合理的。
本文基于貝葉斯更新方法,分級(jí)逐步將ENDF/B-VII.1數(shù)據(jù)庫(kù)提供的累積裂變產(chǎn)額測(cè)量數(shù)據(jù)及獨(dú)立裂變產(chǎn)額服從的二裂變守恒、裂變系統(tǒng)質(zhì)量數(shù)守恒和電荷數(shù)守恒等條件引入?yún)f(xié)方差估計(jì)過(guò)程中,對(duì)ENDF/B-VII.1數(shù)據(jù)庫(kù)中235U熱中子裂變的獨(dú)立裂變產(chǎn)額協(xié)方差進(jìn)行了估計(jì)。結(jié)果表明,與貝葉斯更新前相比,更新后多數(shù)核素獨(dú)立裂變產(chǎn)額的不確定度均出現(xiàn)了不同程度的減小,不確定度約減少了30%。多數(shù)核素的獨(dú)立裂變產(chǎn)額之間為負(fù)相關(guān),且相關(guān)性較弱,極少數(shù)核素的獨(dú)立裂變產(chǎn)額之間存在較強(qiáng)的正相關(guān)。
本文研究是進(jìn)一步開(kāi)展獨(dú)立裂變產(chǎn)額不確定度傳遞研究工作的基礎(chǔ),未來(lái)將基于本文估計(jì)的獨(dú)立裂變產(chǎn)額協(xié)方差矩陣,研究裂變產(chǎn)額不確定度對(duì)反應(yīng)堆燃耗及衰變熱不確定度的貢獻(xiàn)。