亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于貝葉斯更新方法的235U熱中子獨(dú)立裂變產(chǎn)額協(xié)方差估計(jì)

        2021-04-19 04:56:58王毅箴崔夢(mèng)蕾鄔穎杰劉保坤
        現(xiàn)代應(yīng)用物理 2021年1期
        關(guān)鍵詞:質(zhì)量數(shù)產(chǎn)額后驗(yàn)

        王毅箴,崔夢(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)行分析和討論。

        1 獨(dú)立裂變產(chǎn)額的協(xié)方差估計(jì)方法

        1.1 貝葉斯更新方法

        貝葉斯更新方法(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)額不確定度的傳遞。

        1.2 獨(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ǔ)。

        2 估計(jì)結(jié)果

        2.1 獨(dú)立裂變產(chǎn)額調(diào)整結(jié)果與討論

        獨(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

        2.2 獨(dú)立裂變產(chǎn)額相關(guān)性估計(jì)結(jié)果與討論

        隨著累積裂變產(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)估工作是合理的。

        3 結(jié)論

        本文基于貝葉斯更新方法,分級(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)。

        猜你喜歡
        質(zhì)量數(shù)產(chǎn)額后驗(yàn)
        一個(gè)可靠和準(zhǔn)確的光電產(chǎn)額譜模型及應(yīng)用*
        針對(duì)裂變產(chǎn)額和半衰期的燃耗計(jì)算靈敏度和不確定度分析方法
        基于對(duì)偶理論的橢圓變分不等式的后驗(yàn)誤差分析(英)
        貝葉斯統(tǒng)計(jì)中單參數(shù)后驗(yàn)分布的精確計(jì)算方法
        物質(zhì)結(jié)構(gòu)與元素周期律
        一種基于最大后驗(yàn)框架的聚類分析多基線干涉SAR高度重建算法
        一次進(jìn)樣分析氧氣中碳、硫化物和質(zhì)量數(shù)豐度值與比值
        分析儀器(2016年1期)2016-03-04 11:17:00
        裂變產(chǎn)物活度計(jì)算通用程序開(kāi)發(fā)
        一次進(jìn)樣分析氬氣中碳、硫化物和質(zhì)量數(shù)豐度值與比值
        分析儀器(2015年3期)2015-05-09 02:11:01
        Changes in feed intake,nutrient digestion,plasma metabolites,and oxidative stress parameters in dairy cows with subacute ruminal acidosis and its regulation with pelleted beet pulp
        亚洲av综合日韩| 久久中文字幕一区二区| 国产视频在线一区二区三区四区| 国产思思久99久精品| 女同av免费在线播放| 一二三区亚洲av偷拍| 日韩 无码 偷拍 中文字幕| 亚洲精品无码成人a片| 国模无码视频专区一区| 在线视频播放观看免费| 日韩精品熟妇一区二区三区| 欧美人和黑人牲交网站上线| 91亚洲国产成人aⅴ毛片大全| 精品视频一区二区在线观看| 亚洲精品98中文字幕| 无码人妻精品一区二区三区9厂 | 亚洲AV无码成人品爱| 国产精品久久久看三级| 欧美黑人巨大videos精品| 亚洲伊人成综合网| 欧美日韩在线免费看| 日本高清中文一区二区三区| 最新国产不卡在线视频| 成年免费视频黄网站zxgk| 99久久精品国产一区二区蜜芽| 亚洲人成人一区二区三区| av免费网站在线免费观看| 五月天中文字幕日韩在线| 亚洲av无码专区在线播放中文| 91精品国产福利尤物免费| 亚洲天堂av在线免费播放| 亚洲av无码一区二区一二区| 国产精品久久婷婷六月丁香| 国产激情一区二区三区在线蜜臀| 成人久久黑人中出内射青草| 人人妻人人爽人人澡人人| 亚洲欧洲久久久精品| 91国内偷拍精品对白| 四川丰满妇女毛片四川话| 久久亚洲精品无码va大香大香| 国产福利小视频在线观看|