本發明涉及軍控核查技術領域，特別是涉及一種密閉容器核材料質量的計算方法。
背景技術：
在關系世界安全的核裁軍、核安全、核保障三個重要領域里，都面臨著一個共同的科學技術難題，即對貯存在密閉容器內核材料的識別、認證、性能檢測和屬性測量，而質量屬性是其最重要的屬性之一。基于3he探測器的熱中子多重性探測方法(neutronmultiplicitycounting，nmc)是較為經典的方法，但是由于3he材料價格昂貴等原因，逐漸出現了許多新的探測方法。由于近些年3he氣體的產量降低，需求量上升，因此其價格上升較快，目前價格大體為4萬元/升，生產一套nmc探測設備價格大體在百萬元到千萬元級別。但是液體閃爍體材料相對價格便宜，因此此替代方法的經濟適應性較好。而且nmc方法在探測中子時，需要先將中子慢化而后測量，這樣就不能保留中子探測的時間信息。
上述基于3he探測器的熱中子多重性探測方法在實踐中已經得到了一定的應用，但是由于其所用材料3he方面的高成本的缺點以及方法本身慢化中子導致無法保存中子信號的時間信息，因此該方法實用性有待提高。
技術實現要素：
本發明的目的是提供一種密閉容器核材料質量的計算方法，以克服現有基于3he探測器的熱中子多重性探測方法存在的成本高、無法保存中子信號的時間信息的問題，進而提高密閉容器核材料質量的測量精度。
為實現上述目的，本發明提供了一種密閉容器核材料質量的計算方法，所述計算方法包括：
獲取液體閃爍體探測器探測到中子數的中子脈沖時間序列；
利用前景計數器對所述中子脈沖時間序列進行前景重數計數，得到前景重數分布，所述前景重數分布包含與觸發中子信號相關的中子信號的計數和與所述觸發中子信號無關的中子信號的計數；
利用背景計數器對所述中子脈沖時間序列進行背景重數計數，得到背景重數分布，所述背景重數分布僅包含與所述觸發中子信號無關的中子信號的計數；
根據所述前景重數分布和所述背景重數分布確定真符合重數分布，所述真符合重數分布僅包含與觸發中子信號相關的中子信號的計數；
根據樣品參數和所述真符合重數分布確定真符合重數的總體矩，所述總體矩為所述樣品參數與所述真符合重數分布之間的函數關系，所述樣品參數包括自發裂變率、樣品增殖系數和(α，n)反應中子數與自發裂變中子數的比例系數，所述(α，n)反應為重原子核α衰變產生的α粒子碰撞輕原子核發生的反應；
根據所述前景重數分布、所述背景重數分布和所述真符合重數分布確定真符合重數的樣品矩，所述樣品矩為對所述中子脈沖時間序列的計數與所述所述真符合重數分布之間的函數關系；
根據所述總體矩和所述樣品矩確定總中子計數率、二重符合中子計數率和三重符合計數率；
根據所述總中子計數率、所述二重符合中子計數率和所述三重符合計數率求解所述自發裂變率、所述樣品增殖系數和所述比例系數；
根據所述自發裂變率和核材料每秒的自發裂變次數計算核材料的質量。
可選的，所述利用前景計數器對所述中子脈沖時間序列進行前景重數計數，得到前景重數分布，具體包括：
獲取前景時間內所述中子脈沖時間序列中中子信號個數，記為前景重數，所述前景重數對應一前景計數器，所述前景時間為自觸發信號后時刻t+短延時pd至所述時刻t+所述短延時pd+符合計數門寬g的時間段；
每次計數結束所述前景重數對應的前景計數器的計數加1，得到前景計數器組；
根據所述前景計數器組中前景重數的分布規律，確定前景重數分布。
可選的，所述利用背景計數器對所述中子脈沖時間序列進行背景重數計數，得到背景重數分布，具體包括：
獲取背景時間內所述中子脈沖時間序列中中子信號個數，記為背景重數，所述背景重數對應一背景計數器，所述背景時間為自觸發信號后時刻t+長延時ld至所述時刻t+所述長延時ld+所述符合計數門寬g的時間段，所述長延時ld大于所述短延時pd+所述符合計數門寬g；
每次計數結束所述背景重數對應的背景計數器的計數加1，得到背景計數器組；
根據所述背景計數器組中背景重數的分布規律，確定背景重數分布。
可選的，所述根據樣品參數和所述真符合重數分布確定真符合重數的總體矩，具體包括：
根據所述前景重數分布和所述背景重數分布確定樣品發射中子重數的前一階階乘矩h′(1)、二階階乘矩h″(1)和三階階乘矩h″′(1)：
其中，m是增殖系數，α是(α，n)反應中子數與自發裂變中子數的比例系數，vsf,1自發裂變重數的一階階乘矩，vsf,2自發裂變重數的二階階乘矩，vsf,3自發裂變重數的三階階乘矩，vi1誘發裂變重數的一階階乘矩，vi2誘發裂變重數的二階階乘矩，vi3誘發裂變重數的三階階乘矩；
根據所述樣品發射中子重數與所述樣品發射中子重數的前一階階乘矩h′(1)、二階階乘矩h″(1)和三階階乘矩h″′(1)，確定一個初級事件探測中子重數的前三階階乘矩，所述一個初級事件探測中子重數的前三階階乘矩為單次源事件探測重數的前一階階乘矩d′(1)、二階階乘矩d″(1)和三階階乘矩d″′(1)：
d′(1)＝h′(1)e′(1)＝ε(1+κ)h′(1)
d″(1)＝h″(1)(e′(1))2+h′(1)e″(1)＝ε2(1+κ)2h″(1)+2εκh′(1)
d″′(1)＝h″′(1)(e′(1))3+3h″(1)e″(1)e′(1)+h′(1)e″′(1)＝ε3(1+κ)3h″′(1)+6ε2κ(1+κ)h″(1)
其中，e′(1)表示每一個樣品發射中子產生信號數的一階階乘矩，e″(1)表示每一個樣品發射中子產生信號數的二階階乘矩，e″′(1)表示每一個樣品發射中子產生信號數的三階階乘矩，ε表示探測效率，κ表示散射串擾率，所述源事件為自發裂變和(α，n)反應；
根據所述樣品發射中子重數的前一階階乘矩h′(1)、二階階乘矩h″(1)、三階階乘矩h″′(1)和所述初級事件探測中子重數的前三階階乘矩確定真符合重數的總體矩，所述真符合重數的總體矩為：
其中，f(t)表示探測的中子信號時間分布，f(s)也表示探測的中子信號時間分布，s表示時間，fd表示雙重計數率的門寬因子，
可選的，所述根據所述前景重數分布、所述背景重數分布和所述真符合重數分布確定真符合重數的樣品矩，具體包括：
獲取總觸發率s＝(f+fα)d′(1)+sbkg＝(f+fα)ε(1+κ)h′(1)+sbkg；
根據所述前景重數分布獲取前景重數的階乘矩，所述前景重數的階乘矩fk為：
根據所述總觸發率和所述所述前景重數的階乘矩確定真符合重數的樣品矩，所述真符合重數的樣品矩為：
可選的，所述根據所述總體矩和所述樣品矩確定總中子計數率、二重符合中子計數率和三重符合計數率，具體包括：
令所述總體矩等于所述樣品矩，得到總中子計數率s、二重符合中子計數率d和三重符合計數率t分別為：
s＝fε(1+κ)νsf,1(1+α)m
可選的，所述根據所述總中子計數率、所述二重符合中子計數率和所述三重符合計數率求解所述自發裂變率、所述樣品增殖系數和所述比例系數，具體包括：
利用所述所述總中子計數率、所述二重符合中子計數率和所述三重符合計數率計算所述樣品增殖系數m；
利用公式
利用公式
根據本發明提供的具體實施例，本發明公開了以下技術效果：本發明中采用了成本相對更低的有機液閃探測器對觸發中子進行探測，降低了整個技術的成本。并且本發明所依據的探測中子重數數據是未經慢化的中子信號，因此極大地保留了中子信號的時間信息，克服了現有技術中采用基于3he探測器的熱中子多重性探測方法進行密閉容器內核材料質量測量的缺點。本發明提供的基于液閃的快中子多重性分析方法的密閉容器核材料質量的計算方法的準確性、經濟性更好，具有更好的推廣和應用價值。
附圖說明
為了更清楚地說明本發明實施例或現有技術中的技術方案，下面將對實施例中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發明的一些實施例，對于本領域普通技術人員來講，在不付出創造性勞動性的前提下，還可以根據這些附圖獲得其他的附圖。
圖1為本發明實施例提供的密閉容器核材料質量的計算方法的流程圖。
具體實施方式
下面將結合本發明實施例中的附圖，對本發明實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是本發明一部分實施例，而不是全部的實施例。基于本發明中的實施例，本領域普通技術人員在沒有做出創造性勞動前提下所獲得的所有其他實施例，都屬于本發明保護的范圍。
本發明的目的是提供一種密閉容器核材料質量的計算方法，以克服現有基于3he探測器的熱中子多重性探測方法存在的成本高、無法保存中子信號的時間信息的問題，進而提高密閉容器核材料質量的測量精度。
為使本發明的上述目的、特征和優點能夠更加明顯易懂，下面結合附圖和具體實施方式對本發明作進一步詳細的說明。
本發明的基本思想是：首先，以nmc的六條基本假設為依據，即(1)假設所有誘發裂變中子與其對應的源自發裂變或(α，n)反應是同時發生的，這條假設被稱為超裂變假設；(2)假設對于整個樣品空間，探測效率和發生裂變的可能性是均勻的；(3)假設(α，n)反應與自發裂變以及誘發裂變的中子能譜是一樣的，即它們產生的中子有著相同的探測效率和誘發裂變率；(4)假設所有被俘獲的中子都會引起誘發裂變；(5)假設中子多重性分析與裂變事件發出的中子能量無關聯；(6)假設在樣品探測器結構中的中子衰退時間近似服從單一的e指數分布。逐條進行分析，得出前四條假設可以直接套用，后兩條假設進行修改，即(5)假設中子在測量裝置中的衰退時間全部分布于100ns以內。這是因為快中子的探測沒有聚乙烯慢化中子的過程，1mev的中子飛行1m大約需要70ns，而快中子多重性測量(fastneutronmultiplicitycounting，fnmc)測量裝置樣品與探測器距離一般小于30cm，所以可以假設中子衰退時間全部分布于100ns以內，這樣100ns的符合門寬足夠大。(6)假設每個從樣品中出射的中子在第一個探測器中產生信號后都有相同的機率再在第二個探測器中產生信號，但不會引起第三個探測器的響應。然后，通過方程推導得出在fnmc中s、d、t與f、m、α的關系。使用液體閃爍體探測器對未知質量核素進行探測，得到其中子的探測時間序列，由于中子脈沖序列包含的中子數目較多，因此需要對其進行處理。經過處理后得到前景計數器和背景計數器(說明參見資料)的計數重數分布。而后經過處理，得到中子俘獲時間的s、d、t方程，求解出f、m、α，從而求解出有效質量m。由于fnmc方法直接探測未經慢化的快中子，因此保留了中子信號的時間信息。
如圖1所示，本發明提供的密閉容器核材料質量的計算方法，包括：
步驟101：獲取液體閃爍體探測器探測到中子數的中子脈沖時間序列；
步驟102：利用前景計數器對所述中子脈沖時間序列進行前景重數計數，得到前景重數分布，所述前景重數分布包含與觸發中子信號相關的中子信號的計數和與所述觸發中子信號無關的中子信號的計數；
步驟103：利用背景計數器對所述中子脈沖時間序列進行背景重數計數，得到背景重數分布，所述背景重數分布僅包含與所述觸發中子信號無關的中子信號的計數；
步驟104：根據所述前景重數分布和所述背景重數分布確定真符合重數分布，所述真符合重數分布僅包含與觸發中子信號相關的中子信號的計數；
步驟105：根據樣品參數和所述真符合重數分布確定真符合重數的總體矩，所述總體矩為所述樣品參數與所述真符合重數分布之間的函數關系，所述樣品參數包括自發裂變率、樣品增殖系數和(α，n)反應中子數與自發裂變中子數的比例系數，所述(α，n)反應為重原子核α衰變產生的α粒子碰撞輕原子核發生的反應；
步驟106：根據所述前景重數分布、所述背景重數分布和所述真符合重數分布確定真符合重數的樣品矩，所述樣品矩為對所述中子脈沖時間序列的計數與所述所述真符合重數分布之間的函數關系；
步驟107：根據所述總體矩和所述樣品矩確定總中子計數率、二重符合中子計數率和三重符合計數率；
步驟108：根據所述總中子計數率、所述二重符合中子計數率和所述三重符合計數率求解所述自發裂變率、所述樣品增殖系數和所述比例系數；
步驟109：根據所述自發裂變率和核材料每秒的自發裂變次數計算核材料的質量。
下面詳細描述核材料的質量的計算過程：
1.1fnmc方程推導
在fnmc的基本假設下，影響測量結果的未知參數可以簡化為三個：自發裂變率f、樣品增殖系數m和(α，n)反應中子數與自發裂變中子數的比例系數α。本部分將推導這三個樣品參數與重數計數的結果的數學關系，建立fnmc的分析方程。運用fnmc分析方程，可以由測量結果反推三個樣品參數。
由測量結果分析裂變材料的問題可看成參數估計問題。參數估計采用矩估計方法，本實施方式使用階乘矩，分別推出真符合重數的總體矩(樣品參數與真符合重數的關系)和樣本矩(探測重數與真符合重數的關系)的表達式，令總體矩與樣本矩相等，即得到參數估計方程。具體來說就是：由自發裂變中子重數分布和中子在樣品內輸運過程中的中子數變化性質，求得一次自發裂變的樣品發射中子重數分布，再根據探測器探測和計數方式推導出真符合重數的總體矩rk(f,α,m)；依據前景計數分布、背景計數分布與真符合重數分布之間的關系，推導出由探測重數表示的真符合重數樣本矩rk(f,b)。
推導中通常假定所有探測到的中子信號，除本底信號和噪聲外，都可以追溯到一次自發裂變或者是一次由重原子核(heavynuclei)α衰變產生的α粒子碰撞輕原子核發生的(α，n)反應，所以把樣品中發生的自發裂變和(α，n)反應稱作源事件。并假定中子經過一個可能發生增殖的輸運過程后從樣品中發射出去，然后以一定的概率被探測系統探測到。
串擾會導致探測系統給出的中子計數大于實際的中子計數，從而影響對放射源參數的分析。根據钚裂變發射中子重數分布的統計規律，若探測系統的探測效率在10％左右，二重計數大約為一重的十分之一，三重也大約為二重的十分之一。那么若串擾率為1％，就意味著二重計數有10％的偏差，三重計數有20％的偏差(二重計數的中的每個中子都可能產生串擾)。這對于方程求解分析的準確性可能帶來很大影響。所以散射串擾就是一個原來的表達式中沒有考慮，而在fnmc中又無法避免的因素，要把它考慮進分析使用的表達式中。
本文采用概率母函數(probabilitygeneratingfunction，pgf)作為數學工具，使用h(z)表示發射中子重數的概率母函數，對fnmc方程的進行構建。根據概率母函數的性質，x的k階階乘矩e(x(x-1)…(x-k+1))就是x的pgf在1處的k階導數。可用h′(1),h″(1),h″′(1)代替單次裂變中子的“一階矩”、“二階矩”、“三階矩”，即h′(1),h″(1),h″′(1)是發射中子重數的前三階階乘矩。
1.1.1總體矩
即樣品參數與真符合重數的關系。樣品增殖系數m采用階乘矩表達式表示為：
m＝h′(1)(1)
其中，h′(1)表示核材料中每個初級中子所致樣品發射中子數的一階矩。
本部分從物理過程出發，推導f，α，m三個樣品參數與真符合重數之間的函數關系——真符合重數的總體矩rk(f,α,m)。因為由概率母函數(pgf，probabilitygeneratingfunction)求階乘矩只需做求導運算，且母函數處理復合函數的高階導數也特別方便。所以這里將使用概率母函數完成對rk(f,α,m)的推導。
1)源事件的概率母函數
每次發生自發裂變時發射的中子重數是獨立同分布的非負整值隨機變量，用概率母函數pgf表示，表達示為
式中，pi－表示發射i個中子時的概率。
核材料中自發裂變是泊松過程，泊松過程的強度f稱為自發裂變率。這樣，自發裂變產生中子的過程就是一個復合泊松過程。
(α，n)反應與自發裂變無關，也是泊松過程，設其強度為fα，每次(α，n)反應發射中子數為1，用概率母函數pgf表示為
sα(z)＝z(3)
將自發裂變與(α，n)反應統稱為源事件，則源事件可以看作是強度為f+fα的泊松過程，每次發射的中子數的概率母函數pgf為
根據pgf的性質和參數α的定義式，一個源事件發射中子數的一階矩、二階矩、三階矩分別為
由于自發裂變中子階乘矩是可以根據核數據計算的常量，所有公式中將使用νsf,1,νsf,2,νsf,3代替ssf′(1),ssf″(1),ssf″′(1)。
2)樣品發射中子數的概率母函數
若源事件產生1個中子，從樣品中發射出來的中子數的pgf記為h1(u)。那么，如果源事件中產生n個中子，則根據pgf的性質，發射中子數的pgf為
hn(z)＝[h1(z)]n(8)
將樣品中一個中子發生誘發裂變的概率記為pif，誘發裂變產生i個中子的可能性記為pifν(i)，則這i個中子所致樣品發射中子數的pgf也為hi(z)。一次誘發裂變產生中子數的pgf為
假設中子在樣品中只要被吸收就會誘發裂變，忽略中子在樣品中被吸收而沒有誘發裂變的可能性，根據pgf的性質
為書寫方便，后文統一用h(z)代替h1(z)，即
h(z)≡h1(z)＝(1-pif)z+pifu[h(z)](11)
式(11)給出了當源事件發射一個中子時，樣品所發射中子數h(z)、誘發裂變概率pif和誘發裂變發射中子數u(z)之間的關系。考慮源事件發射的中子數分布，根據pgf的性質，一次源事件產生的發射中子數pgf為
h(z)＝s[h(z)](12)
3)源事件產生的發射中子數的階乘矩
根據pgf的性質，由式(12)可得源事件對應的樣品發射中子數的一階階乘矩、二階階乘矩、三階階乘矩為
h′(1)＝s′(1)h′(1)(13)
h″(1)＝s″(1)[h′(1)]2+s′(1)h″(1)(14)
h″′(1)＝s″′(1)[h′(1)]3+3s″(1)h″(1)h′(1)+s′(1)h″′(1)(15)
對式(11)，求z＝1處的導數，可得
考慮到誘發裂變中子階乘矩是可以由核數據計算得到的常量，式中用νi1代替u′(1)，下文中同樣使用νi2,νi3代替u″(1),u″′(1)。
h′(1)是每個源中子在輸運過程中的增殖期望，所以設樣品發射中子數與源中子數之比增殖參數m為
對式(11)，求z＝1處的二階導數，可獲得中子增殖的二階矩：
對式(11)，求z＝1處的三階導數，可得中子增殖的三階矩：
將式(17)、式(18)、式(19)代入式(13)、式(14)、式(15)，得到樣品發射中子數前一階矩、二階矩、三階矩為
4)真符合重數的總體矩
先假定同一個源事件的中子中有n個中子被探測到，在其中隨機選取一個中子作為觸發信號，推導這個中子觸發信號的真符合計數重數的概率分布p(i|n)。
假設源事件發生于t＝0時刻，探測的中子信號時間分布為f(t)。符合計數門寬為g，短延時為pd。
在(t,t+dt)內能夠探測到1個中子的可能性為
p(t|n)＝nf(t)dt(23)
以它作為觸發信號符合計數門為從t+pd時刻到t+pd+g時刻，其余n-1個中子信號中的任意一個處于此符合計數門中(記為事件a)，對此次觸發的計數重數產生貢獻)的可能性為
因而在(t,t+dt)內探測到1個中子，并且余下的n-1個中子在符合計數門中有i個的可能性為
在整個時間范圍對t做積分，即可得到隨機選取n個信號中的一個作為觸發信號，所記重數為i(符合計數門中有i個中子信號)的可能性p(i|n)，
快中子測量的探測重數與熱中子的不同。對于熱中子探測，一個發射中子被探測并產生一個信號的可能為ε，不被探測的可能性是(1-ε)，所以經典nmc測量中每個發射中子產生的中子信號重數的pgf為e(z)＝(1-ε)+εz，而快中子多重性探測還需要考慮散射串擾的影響，根據fnmc第(6)條基本假設，fnmc測量中每一個樣品發射中子產生信號數的pgf為
e(z)＝(1-ε)+ε(1-κ)z+εcz2(27)
根據pgf的性質，一次源事件的探測重數為：
d(z)＝h(e(z))(28)
記一次源事件探測到的中子數為n的可能性為d(n)，則對于每個探測到的中子都會觸發一次計數，則某一個觸發中子屬于探測到n個中子的源事件的可能性為
式中，d(n)―探測重數分布
由公式(28)、(29)，推導出單次源事件探測重數的前三階矩測量方程。即一個初級事件探測中子重數的前三階矩為：
某次計數重數為i的概率為
由式(31)、式(32)、式(33)、式(34)，可求得出真符合重數總體矩r(i)(信號觸發關聯重數)的一、二階矩分別為：
式中，fd―雙重計數率(doubles)門寬因子；ft―三重計數率(triples)的門寬因子。可以由式(36)、(37)計算。
1.1.2樣品矩
即測量計數與真符合重數的關系。探測到的中子既包括源事件直接或間接誘發裂變產生的中子，又有與源事件無關的本底中子。前景計數器以探測到的中子為觸發，總觸發率為：
s＝(f+fα)d′(1)+sbkg＝(f+fα)ε(1+κ)h′(1)+sbkg(38)
前景重數的階乘矩為：
式中bk―本底中子觸發重數，即偶然符合重數，與背景重數分布相同；gk―源事件中子觸發重數，包含真符合重數rj與偶然符合重數bk。
所以源事件中子觸發重數的階乘矩公式為，
對式(39)兩邊同時減bk，再將式(40)代入得
依次k＝1,2將代入式，可以獲得真符合重數的一、二階矩r'1，r'2使用fk,bk的形式表達。對應的真符合重數樣本矩為：
1.1.3參數估計方程組
令真符合重數的總體矩(式(34)、(35))，與樣本矩(式(42)、(43))相等，即可得到兩個矩估計方程，而式反應了探測到總中子數(觸發數)與樣本參數之間的關系。聯立三個方程，并把測量量與物理參數量分別放在等式兩邊，即可得到考慮散射串擾推導出的fnmc測量方程為：
s-sbkg＝(f+fα)ε(1+κ)h′(1)(44)
式(44)、式(45)、式(46)左側就是fnmc中常用到的測量量singles，doubles和triples，簡寫為s，d，t。
將式(20)、式(21)、式(22)代入式(44)、式(45)、式(46)得：
s＝fε(1+κ)νsf,1(1+α)m(47)
1.2fnmc方程求解
通過對公式(47)、(48)、(49)的相互代入消除f和α，得到方程，m的解可以通過求解該一元三次方程獲得。
a+bm+cm2+m3＝0(50)
其中，
求解出m后，自發裂變率f由下式求出
(α，n)反應系數α為
一旦求出裂變率f后，根據每克240pu平均每秒發生473次裂變的物理性質，可計算等效240pu質量。钚樣品的有效質量可以由式給出。然后用有效質量除以同位素組成得到總質量。
下面給出本發明的計算方法在具體實例中的應用。
現以bc501-a液體閃爍體探測器并結合g4模擬仿真工具箱為例對本發明方法進一步說明。
使用g4仿真工具箱模擬bc501-a液體閃爍體探測器對240pu中子源進行探測，模擬測量82.9g金屬钚(90％240pu，10％239pu，19.8g/cm3)。模擬測量中模擬了4,000,000次初級行為(全部為240pu自發裂變)，輸出的數據包含922,178個探測到中子的時間。按時間先后將數據排序后，為保證數據的準確性，取用中間90％的數據(從第46,108個到第875,522個)，對應的測量時間長度為101.9s。按照上述步驟對其探測結果結果進行處理。
根據中子的能譜模擬發射中子，使其發射概率符合中子的能譜。
根據概率母函數，計算單次源事件探測重數的前三階矩。并得到金屬钚模擬測量重數計數分布。(為簡化計算且在誤差允許范圍內，fd和ft簡單取為1)
金屬钚模擬測量重數計數分布表
由fnmc測量方程及相關參數求解s、d、t的值。
d＝s(f1-b1)＝1.03×103
計算得到有效質量
因樣品是82.9g金屬钚(90％240pu，10％239pu)的，考慮240pu的同位素組成為90％，最后可求出钚的質量相對偏差為8.1％，解算效果較好。
本說明書中各個實施例采用遞進的方式描述，每個實施例重點說明的都是與其他實施例的不同之處，各個實施例之間相同相似部分互相參見即可。
本文中應用了具體個例對本發明的原理及實施方式進行了闡述，以上實施例的說明只是用于幫助理解本發明的方法及其核心思想；同時，對于本領域的一般技術人員，依據本發明的思想，在具體實施方式及應用范圍上均會有改變之處。綜上所述，本說明書內容不應理解為對本發明的限制。