基于SNIP算法扣除_能谱本底的探讨及应用.docVIP

基于 SN IP算法扣除 γ能谱本底的探讨及应用彬 3尹旺明 ,刘宏章 ,汤(东华理工大学核工程技术学院 基于 SN IP算法扣除 γ能谱本底的探讨及应用 彬 3 尹旺明 , 刘宏章 , 汤 (东华理工大学核工程技术学院 ,江西 抚州 344000 ) 摘 要 :在自然 γ能谱全谱分析中 ,为了分离无用信息及提高解释精度 ,须进行 γ全谱连续本底的自动扣除。介绍了 γ能 谱本底的来源及特征 ,探讨了 SN IP算法扣除γ能谱本底的数学原理及过程 ,重点开展了 SN IP方法在 γ能谱本底扣除中的 初步应用。通过对三个不同能量的全能峰的比较分析 ,指出应根据 γ能谱中谱峰宽度的差异 ,选择合适的变换宽度 ,以便 获得最佳的全谱本底扣除效果。实验证明 ,该技术可应用于自然 γ能谱测井的软件处理平台中。 关键词 :γ能谱 ;本底 ; SN IP算法 ;峰面积 中图分类号 : P631. 6 文献标识码 : A 文章编号 : 1674 23504 ( 2009 ) 03 20245 204 核工业的迅速发展 ,特别是核电的异军突起 , 带动了核资源勘查和环境监测技术的发展 ,表现在 核仪器硬件的不断改进和软件中谱分析方法的不 断更新上 (李必红等 , 2008 ) 。在进行自然 γ能谱全 谱分析前需对全谱做本底扣除 (庞巨丰 , 1990; 腰善 丛 , 2002 ) ,正确的扣除本底方法应满足三个特征 (Mo rhac e t a l. , 1997 ) : 首先是简单 、快速且特征参 数的个数少 ;其次是易实现计算机自动处理 ; 第三 , 能剔除康普顿散射射线对谱峰计数的贡献。在 γ 能谱的全谱本底评估方法中 , SN IP 算法已得到国 外核数据计算系统的应用。因此 , 本文引用 SN IP 算法开展了 γ全谱本底扣除的应用实验 ,并对实验 结果进行了分析。 测得的全能峰将出现低能拖尾现象。第三 ,天然本 底 ,主要来自仪器的电子学噪声、宇宙射线 、大气中 氡子体及探头材料中的放射性核素产生的辐射 。 SN IP算法的原理及扣除本底的过程 R yan等 ( 1988 ) 最早分析了 SN IP 算法扣除 γ 能谱本底的原理 , Mo rhac 等 ( 1997 ) 对该算法进行 了改进。在实际扣除 γ能谱本底的过程中 ,一般首 先用 LL S运算符 ( H amp ton e t a l. , 1994 ) 对每道计 数值进行变换 ,该变换公式如下 : 2 v ( i) = ln [ ln ( y ( i) + 1 + 1 ) + 1 ] ( 1 ) 式中 i为道址数 ; y ( i) 为道址 i对应的计数; v为运算 结果保存向量。 ( 1 ) 式中的自然对数运算符可处理高计数谱 峰 , 而求平方根运算符对增强弱峰很敏感 。 1 γ能谱本底的来源及其特征 γ能谱的全谱本底主要来自三个方面 (庞巨 然后根据 SN IP算法依次迭代计算 v i , v ( i) 1 ( ) 2 丰 , 1990 ) : 第一 , 康普顿效应引起的连续谱分布 。 对单能 γ能谱 ,康普顿分布在低能区是随道址缓慢 变化的平滑谱。对于具有多种能量 γ峰的复杂 γ 能谱 ,所有 γ峰的康普顿分布均叠加在一起 。因 此 ,本底的总特征将是一个随能量而平滑下降的函 数。第二 ,γ射线在探测器灵敏体积内的小角度散 射 。此过程使 γ光子仅损失很小的能量 ,因而 ,所 直到 v i , m 为人为设定的参数。即第 p次迭代中 m ( ) 第 i道的 vp ( i) 值确定方式为通过比较 vp - 1 ( i) 跟 vp - 1 ( i + p) , vp - 1 ( i - p) 的平均值的大小 , 并取两者 中的小值 。其数学原理如下 : 1 vp ( i) = m in vp - 1 ( i) , 2 [ vp - 1 ( i + p) + vp - 1 ( i - p) ] ( 2 ) ( 2 ) 式中的 p值有两种变化方式 : 一是从 1 开 始每次迭代运算后加 1直到等于给定的 m 值 (变换 窗宽逐增法 ) ; 二是从给定的 m 值开始每次迭代运 算后减 1直到等于 1 (变换窗宽逐减法 ) 。 当 vm ( i) 计算完后 , 再进行反 LL S运算即可求 收稿日期 : 2009202 226 基金项目 :国家“863 ”项目 ( 2007AA 06 Z111 ) 作者简介 : 尹旺明 ( 1982—) ,男 ,硕士生 ,研究方向为核探测技术与智能 仪器 。 3 通讯作者 : 汤 彬 ( 1963 —) ,男 ,博士 ,教授 ,研究 生导师 ,主要从事放射性地球物理及