逆矩陣解析γ譜

1、逆矩陣法解析譜一、 實(shí)驗(yàn)?zāi)康模海?）了解逆矩陣法解析復(fù)雜譜的基本原理；（2）測(cè)量137Cs、54Mn、60Co混合樣品中各核素的未知活度。二、 實(shí)驗(yàn)儀器：NaI(Tl)閃爍體探測(cè)器一個(gè)、光電倍增管一個(gè)、高壓直流電源一個(gè)、線形放大器1個(gè)、UMS微機(jī)多道采集系統(tǒng)、樣品臺(tái)和屏蔽體1套、同軸電纜若干已經(jīng)活度的137Cs、54Mn、60Co源各一個(gè)、含有137Cs、54Mn、60Co的混合樣品源一個(gè)、無(wú)源塑料片一個(gè)三、 實(shí)驗(yàn)原理：（1） 逆矩陣算法首先測(cè)定第j種核素標(biāo)準(zhǔn)源的能譜，記錄第i道域的計(jì)數(shù)，進(jìn)而得到第j種核素對(duì)第i道域的響應(yīng)系數(shù)aij，其意義是第j中核素的衰變率在i道域上所引起的計(jì)數(shù)率。然后測(cè)量

2、已知核素種類未知活度的混合源樣品，記錄各個(gè)第i道域上的計(jì)數(shù)，根據(jù)標(biāo)準(zhǔn)源的響應(yīng)系數(shù)，得出各核素源對(duì)道域i的計(jì)數(shù)貢獻(xiàn)的線性方程組，令，即有，所以當(dāng)測(cè)得響應(yīng)系數(shù)矩陣A和混合源道域計(jì)數(shù)矩陣M后，對(duì)矩陣A求逆再乘以M即可得混合源各核素放射性活度矩陣X，即（2） 逆矩陣法解譜的條件： 混合樣品中的核素都是一致的，可以選定相應(yīng)的標(biāo)準(zhǔn)樣品和建立標(biāo)準(zhǔn)譜； 測(cè)量標(biāo)準(zhǔn)譜和混合樣品譜時(shí)，譜儀的測(cè)量條件一致，儀器要求穩(wěn)定； 各單一核素放射性線形疊加時(shí)，探測(cè)系統(tǒng)的脈沖幅度不隨計(jì)數(shù)率改變； 混合樣品中各核素具有各自的能譜，能夠選出表征各自能量的特征峰；（3）逆矩陣法解譜得誤差：逆矩陣法對(duì)某核素的分析精度，取決于他的特征道域

3、中該核素計(jì)數(shù)在樣品總計(jì)數(shù)中所占的份額。對(duì)各種成份放射性含量比相近的樣品，解譜誤差可達(dá)2%左右，而對(duì)弱成分，誤差會(huì)顯著增大。被測(cè)樣品在各特征道域上的計(jì)數(shù)和譜儀本底都會(huì)有統(tǒng)計(jì)漲落，并通過(guò)逆矩陣計(jì)算互相傳遞。如果忽略響應(yīng)系數(shù)的誤差，則成分j的活度xj的標(biāo)準(zhǔn)誤差為其中，Bj為本底計(jì)數(shù)統(tǒng)計(jì)漲落對(duì)誤差的貢獻(xiàn)，bi為第i道域本底計(jì)數(shù)率，t為樣品測(cè)量時(shí)間。同時(shí)，結(jié)果的誤差會(huì)隨著混合樣品中核素的數(shù)目增加而增大，因此實(shí)際應(yīng)用中混合樣品中的核素種類不超過(guò)5種。四、 實(shí)驗(yàn)步驟： 檢查儀器連接，預(yù)熱10min，打開微機(jī)，進(jìn)入多道UMS； 將60Co源置于樣品臺(tái)上，調(diào)整NaI（Tl）譜儀工作電壓至600V，調(diào)整放大倍數(shù)，

4、使1.33MeV對(duì)應(yīng)的光電峰位于800道左右，放大倍數(shù)最后選擇1.86； 測(cè)量137Cs、54Mn、60Co標(biāo)準(zhǔn)樣品的能譜，測(cè)量時(shí)間1800秒。在三個(gè)脈沖高度譜圖上分別記錄三個(gè)光電峰位對(duì)應(yīng)道址左右18道范圍內(nèi)的面積計(jì)數(shù)； 分別對(duì)混合樣品和本底采集1800s，記錄第步中所確定的三個(gè)道域中的計(jì)數(shù)； 采用逆矩陣法求出混合樣品中三種源的活度，并求出方差； 以第步測(cè)得的4個(gè)光電峰位和相應(yīng)的能量做能量刻度曲線，并用137Cs能譜確定譜儀能量分辨率； 降下高壓，將放射源送回鉛盒，歸整儀器。五、 NaI閃爍體探測(cè)系統(tǒng)線路框圖圖1、NaI閃爍體探測(cè)系統(tǒng)線路框圖六、實(shí)驗(yàn)結(jié)果和分析（1）系統(tǒng)工作狀態(tài)參數(shù)；表1、系統(tǒng)

5、工作狀態(tài)參數(shù)表高壓600V放大增益因子1.86單次測(cè)量譜采集時(shí)間1200s（2）能量刻度曲線和閃爍體探測(cè)器的能量分辨率：采用137Cs、54Mn、60Co的四個(gè)峰能量進(jìn)行刻度測(cè)量137Cs的峰能量下閃爍體探測(cè)器的能量分辨率表2、能量刻度曲線數(shù)據(jù)能量E/MeV0.6620.8351.171.33道址V380492701804線性擬合結(jié)果：E=B*V+A MeVParameterValueError-A0.059860.00481B0.001587.79135E-6-RSDNP-0.999980.00264<0.0001-圖2、閃爍體探測(cè)器能量刻度曲線表3、137Cs的峰能量下系統(tǒng)分辨率峰位

6、能量半高寬半高寬能量分辨率3800.662MeV44.10.06968MeV11.7%（3）標(biāo)準(zhǔn)源活度：表4、3種標(biāo)準(zhǔn)源活度核素半衰期/d初始活度/Bq制作日期衰變時(shí)間/d當(dāng)前活度/Bq137Cs110122.381x10548931.74985x10554Mn312.27.77 x1052004.1.2318721.2173x10460Co1925.21.597x10547262.9130x104（4）標(biāo)準(zhǔn)源響應(yīng)系數(shù)：表5、3種標(biāo)準(zhǔn)源對(duì)3個(gè)指定道域的響應(yīng)系數(shù)道域362-398474-510786-822本底計(jì)數(shù)769474160137Cs總計(jì)數(shù)3973097366558606凈計(jì)數(shù)39703

7、28361818446響應(yīng)系數(shù)22.68952.0677x10-14.8267x10-254Mn總計(jì)數(shù)49071370998168凈計(jì)數(shù)483023705248響應(yīng)系數(shù)3.968030.43826.5719x10-460Co總計(jì)數(shù)244185290415331773凈計(jì)數(shù)243416289941331613響應(yīng)系數(shù)8.35629.953311.3839（5）逆矩陣法解混合樣品中三種源的活度表6、混合源在3個(gè)道域上的計(jì)數(shù)362-398474-510786-822本底總計(jì)數(shù)凈計(jì)數(shù)本底總計(jì)數(shù)凈計(jì)數(shù)本底總計(jì)數(shù)凈計(jì)數(shù)76928687872868018474432018431544160224226224066由表6得混合源道域計(jì)數(shù)矩陣由表5得3種標(biāo)準(zhǔn)源響應(yīng)系數(shù)矩陣?yán)肕ATLAB求逆矩陣A-1，乘以M得到混合源中3種放射核素的活度矩陣 Bq即混合樣品中，137Cs的活度為Bq，54Mn的活度為Bq， 60Co的活度為B