基于溴化鑭探測(cè)器的本底扣除方法研究

侯鑫等

摘 要：溴化鑭探測(cè)器由于其優(yōu)秀的能量分辨率和探測(cè)效率，受到核行業(yè)工作者的高度關(guān)注。文章針對(duì)LaBr3探測(cè)器在天然γ放射性測(cè)量中譜線本底扣除方法進(jìn)行了比較分析，得出結(jié)論：在實(shí)際工程應(yīng)用中，選擇傅里葉本底扣除方法結(jié)合直線法做本底扣除，效果最佳。

關(guān)鍵詞：溴化鑭探測(cè)器；γ能譜；本底扣除

引言

自2000年首次被推出以來(lái)，LaBr3：Ce3+晶體就以其優(yōu)秀的能量分辨率、良好的能量線性響應(yīng)、較強(qiáng)的光輸出能力引起國(guó)內(nèi)外的普遍關(guān)注并迅速發(fā)展起來(lái)。文章針對(duì)LaBr3探測(cè)器在天然γ放射性測(cè)量中譜線本底扣除方法進(jìn)行了比較分析。

1 譜線測(cè)量

文章中數(shù)據(jù)使用成都中福科技有限公司與成都理工大學(xué)聯(lián)合研發(fā)的IED-3000B型數(shù)字化γ能譜儀系統(tǒng)，搭配法國(guó)圣戈班公司生產(chǎn)的BrilLance 380型號(hào)的溴化鑭晶體探頭，在四川省德陽(yáng)市羅江縣四川省核工業(yè)輻射測(cè)試防護(hù)院計(jì)量模型站的標(biāo)準(zhǔn)飽和模型上測(cè)得。其中本底模型測(cè)量400s，其他模型均測(cè)量200s。模型中的物質(zhì)含量如表1所示。

2 試驗(yàn)方法

為了檢驗(yàn)本地扣除方法對(duì)溴化鑭探測(cè)系統(tǒng)的適應(yīng)性，使用在本底模型、U模型、K模型、Th模型上測(cè)得的數(shù)據(jù)標(biāo)定儀器，使用混合模型上測(cè)得的數(shù)據(jù)進(jìn)行誤差檢驗(yàn)，使用《DZ-T 0205-1999地面伽瑪能譜測(cè)量技術(shù)規(guī)范（行業(yè)標(biāo)準(zhǔn)）》附錄B中推薦的逆矩陣法求解核素含量，以誤差大小判斷本底扣除方法的適應(yīng)性。

2.1 SNIP法本底扣除

SNIP法在實(shí)際扣除γ能譜本底的過(guò)程中，一般首先用LLS運(yùn)算對(duì)每道計(jì)數(shù)進(jìn)行變換，其公式如下：

式中i為道址數(shù)；y（i）為道址i對(duì)應(yīng)的計(jì)數(shù)；v（i）為運(yùn)算結(jié)果保存向量。（1）式中的自然對(duì)數(shù)運(yùn)算符可處理高計(jì)數(shù)譜峰， 而求平方根運(yùn)算符對(duì)增強(qiáng)弱峰很敏感。做此變換能夠有效的減少迭代次數(shù)。

SNIP算法是通過(guò)引入剝離能量譜中變化迅速的特征信息，比較第i道道址信息yi和它鄰近道址的信息。比較vi和它的兩個(gè)鄰近道址的平均值mi大小：

當(dāng)vi的值大于平均值mi，用平均值mi替換道址i的信息，將這種變換在所有道址上運(yùn)行一遍，此時(shí)，能量譜其它位置保持不變，但特征峰峰位置的峰幅度會(huì)降低。此方法導(dǎo)致譜的光滑，趨于連接兩個(gè)局部最小值，最終剝離收斂，留下光滑的本底。再將經(jīng)迭代最終得到的數(shù)據(jù)反變換既可得到本底譜。

2.2 傅里葉變換本底扣除法

采用數(shù)字信號(hào)處理中頻譜的處理方法對(duì)離散的譜數(shù)據(jù)進(jìn)行離散傅里葉變換，將頻域中的低頻部分適當(dāng)?shù)目鄢瑥亩_(dá)到背景扣除的效果。

N是熒光儀的總道數(shù)，這里為1024道。

通過(guò)公式（3）和公式（4）可以實(shí)現(xiàn)離散儀器譜數(shù)據(jù)在時(shí)域和頻譜間相互變換。

由于特征峰為高斯脈沖形狀，其頻率中也含有低頻成分，與散射背景在頻譜中不能完全區(qū)分，因此無(wú)法直接使用低通濾波器來(lái)估計(jì)背景，但是可以使用迭代算法實(shí)現(xiàn)，實(shí)現(xiàn)方法如下：

（1）對(duì)第m次迭代得到的譜數(shù)據(jù)xm（n）進(jìn)行快速傅里葉變換，得到頻譜Xm（k）。

（2）采用理想低通濾波器對(duì)Xm（k）進(jìn)行濾波，得到低頻譜Xm'（k）。

（3）對(duì)Xm'（k）進(jìn)行快速傅里葉逆變換，得到譜線xm'（n）。

（4）因?yàn)楸尘白V每道的計(jì)數(shù)應(yīng)該小于該道原始譜計(jì)數(shù)并且大于零，所以可以構(gòu)造第m+1次迭代的譜線xm+1（n），如式（5）所示。

調(diào)整濾波器閾值，重復(fù)執(zhí)行以上步驟，當(dāng)結(jié)果趨向于穩(wěn)定時(shí)結(jié)束，得到的xm+1（n）即為所求的背景。

3 測(cè)量結(jié)果及分析

根據(jù)《DZ-T 0205-1999地面伽瑪能譜測(cè)量技術(shù)規(guī)范（行業(yè)標(biāo)準(zhǔn)）》中推薦的伽瑪能譜測(cè)量譜段，對(duì)在羅江鄢家鎮(zhèn)鈾、釷、鉀、混合四個(gè)標(biāo)定模型上測(cè)得的γ能譜進(jìn)行特征峰面積計(jì)算，結(jié)果見(jiàn)表2。用逆矩陣解譜法求得鈾釷混合模型含量，結(jié)果見(jiàn)表3。

由表1-表3可以看出，對(duì)于K元素來(lái)說(shuō)，未扣除本底和SNIP法扣除本底，誤差均達(dá)到30%以上，直線本底扣除法雖然誤差有所降低，但任有18%，誤差較大；而采用傅里葉本底扣除方法可以使鉀元素含量誤差降低到3.3%。行業(yè)標(biāo)準(zhǔn)中要求鉀含量低于5%時(shí)，一般精度儀器的絕對(duì)誤差小于1.0，高精度儀器的絕對(duì)誤差小于0.5，由表3可看出，對(duì)于K元素的測(cè)量，傅里葉本底扣除遠(yuǎn)遠(yuǎn)滿足行業(yè)要求。對(duì)于U來(lái)說(shuō)，所有方法的誤差均在2.0%～5.3%之間，傅里葉本底扣除的效果最好；對(duì)于Th來(lái)說(shuō)，所有方法的誤差均在1.7%～5.9%之間，直線法扣除本底的效果最好；行業(yè)標(biāo)準(zhǔn)中要求如果他們含量大于15%時(shí)，一般精度儀器的相對(duì)誤差小于20%，高精度儀器的絕對(duì)誤差小于15%，其測(cè)量結(jié)果均滿足行業(yè)要求。

溴化鑭本身的138La元素產(chǎn)生的γ射線會(huì)對(duì)環(huán)境中的鉀峰有很大影響，在本底模型上測(cè)量的鉀峰計(jì)數(shù)有30～60%來(lái)源于138La。SNIP法由于其方法本身的原因，針對(duì)平緩，起伏不是很大的本底扣除效果較好。對(duì)于溴化鑭探測(cè)器來(lái)說(shuō)，譜線中的本底除了環(huán)境中的自然本底意外，還有溴化鑭自身所帶的本底，而且有明顯的譜峰形成，用SNIP法做本底扣除效果就不理想。對(duì)于Th峰來(lái)說(shuō)，由于其能量在自然界中比較高，環(huán)境中其他元素對(duì)其影響很小，可以用直線法扣除本底，效果較好。

總上所述，在實(shí)際工程應(yīng)用中，選擇傅里葉本底扣除方法做LaBr3探測(cè)系統(tǒng)的本底扣除，計(jì)算效果最佳，使K的絕對(duì)誤差達(dá)到了0.08，而U、Th的相對(duì)誤差分別達(dá)到2.0%、5.9%。對(duì)于Th元素含量的計(jì)算，也可以采用直線法做本底扣除。

參考文獻(xiàn)

[1]MAURIZIOCONTI， LARSERIKSSON， CHARLESHAYDEN. Monitoring energy calibration drift using the scintillator background radiation[J]. IEEE Transactions on Nuclear Science， 2011（58）： 687-694.

[2]ANILKUMAR G， MAZUMDAR I， GOTHE D A. Ef？ciency calibration and simulation of a LaBr3（Ce） detector in close-geometry[J].Nuclear Instruments and Methodsin Physics Research， 2009（609）： 183-186.

[3]高鑫，何元金.LaBr3Ce3閃爍晶體研究進(jìn)展[J].核電子學(xué)與探測(cè)技術(shù)，2010，30（1）：5-11.

[4]高峰，張建國(guó)，楊翊放，等.LaBr3：Ce閃爍體探測(cè)器自發(fā)本底譜研究[J].核電子學(xué)與探測(cè)技術(shù)，2010，32（5）：514-517.

作者簡(jiǎn)介：侯鑫（1991，5-），男，山西臨汾人，成都理工大學(xué)在讀碩士研究生，研究方向：核技術(shù)及應(yīng)用。