差值算法在γ能谱峰面积计算中的应用

《差值算法在γ能谱峰面积计算中的应用》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、广东微量元素科学CUANGD0NGWEIUANGYUANSUKEXUE第20卷第12期文章编号：1006—446X(2013)12—0016—04差值算法在能谱峰面积计算中的应用高红赖万昌李思威曹发明(成都理工大学核技术与自动化工程学院，四川成都610059)摘要：野外多道能谱仪在计算能谱峰面积时，左右边界道址取整数往往产生较大的计算误差，引入拉格朗日插值法对非整数道的道址进行插值，能有效减低这种误差。实验表明，对1024道的NaI(哦)能谱仪，计数率最大校正值可达1．08cps，测量误差则可降低10．61％。关键词：能谱；谱处理；误差；拉格朗日插值法中图分类号：T

2、L817．2文献标识码：AInterpolationAlgorithmApplicationtoCalculateGammaSpectralPeakAreaGAOHong，LAIWanchang，LISiwei，CAOFaming(CollegeofAppliedNuclearTechnologyandAutomationEngineering，ChengduUniversityofTechnology，SichuanChengdu，610059，China)Abstract：Anewtreatmentofpeakareacalculationistheintrod

3、uctionoflagrangeinterpolationchannelofnonintegerchannelassignment，andthemeasurementresultscomparedwiththetraditionalmethodoferror．Experimentsshowthat：the1024NaI(T1)gammaspectrometer，countratemaximumcorectionvalueupto1．08cps．themaximumerrorofmeasurementresultscanbereduced10．61％．Keywords

4、：gammaspectrum；spectrumprocessing；eror；lagrangeinterpolationmethod天然放射性元素K、u、Th的能谱可用各种能谱仪测得，仪器标定后，即可将能谱特征峰面积转换成K、u、Th含量。因此，能谱特征峰峰面积的计算精度非常重要，目前峰面积计算法一般有计数法和函数拟合法等⋯。计数法主要包括总峰面积法、科沃尔峰面积法和瓦森峰面积法；函数拟合法包括扣除本底后净峰拟合法以及峰形函数加基底函数拟合的峰面积法J。在峰面积计算中，左右边界往往是以道址确定的，由于道址的整数特性，往往与能谱的能量边界不能完全一致，因此给峰面积的计

5、算带来计算误差。作者利用插值算法对此进行了插值计算，效果较好。1^y能谱法与能谱仪的能区差异实验仪器是IED一3000B型数字化多道伽玛能谱仪，按照《地面伽玛能谱测量技术规范》设收稿日期：2013一I1—05项目基金：国家863计划项目(2012AA06180303)作者简介：高红(1988～)，女，江西南昌人，成都理工大学核技术与自动化工程学院硕士研究生，核资源与核勘查工程专业。·16·广东微量元素科学GUANGDONGWEIHANGYUANSUKEXUE第20卷第12期定的能区分别为：K的1．46MeV特征峰的边界为1．37MeV、1．57MeV，U的1．76M

6、eV特征峰的边界为1．66MeV、1．86MeV，Th的2．62MeV特征峰的边界为2．41MeV、2．81MeV。采用能量线性刻度法，即使用已知两个特征峰的能量进行刻度，IED一3000B型数字化多道伽玛能谱仪刻度公式如下。．一E×(ch—ch)+(1)实际计算时是根据能区边界计算各谱峰的左右边界道址。依据(1)式可得。Ph一hc×(E—E)+c(2)c的计算结果是浮点数，而在^y能谱仪中c^是整数。因此，必然给谱峰面积的计算带来计算误差2拉格朗日线性插值法的原理对于+1个互不相同的插值节点(i，Y(j))，i=0，1，2，⋯，n，L()为不超过2次多项式，对每个

7、插值节点i可作出相应的n次插值基函数，篙专案羔鲁(3)简化得：L()=f())，(。)(4)式(4)为拉格朗日插值多项式。在能谱中Yi表示第i道的计数。3插值法计算峰面积在特征峰峰区左、右边界应用线性插值计算方法对其进行修正，公式如下，mHS=+((Yit(m，)+1一Yit(m，))×(m一int(mL))+Yi(m￡))×(1一mL+int(，n))+((Yi(m)+1一Yint(mH))×(，nH一(m))+Y州(m))×(mH—n(m))一Y(m)(5)其中int是数值取整的运算，m、m分别为根据公式(2)算出的边界道址，Y表示第i道的计数。对峰面积进行