本发明专利技术提出了一种采用溴化镧探测器测量环境中的钾40的方法，属于辐射探测和辐射环境监测方法领域。该方法包括对溴化镧探测器自身放射性进行标定，得到该探测器本底谱响应y；采用溴化镧探测器对环境中的40K进行测量，并拟合得到能谱响应y*；根据溴化镧探测器本底谱响应y和能谱响应y*，计算得到环境中40K活度。本发明专利技术方法实现了对环境中的40K的精确测量。
【技术实现步骤摘要】
本专利技术属于辐射探测和辐射环境监测方法领域，具体涉及采用溴化镧探测器测量环境中的钾40(40K)的方法。
技术介绍
掺铈溴化镧(LaBr3:Ce)具有出色的性能：1.原子序数高(有效原子量为46.9接近碘化钠的50)、密度大(密度5.29g/cm3大于碘化钠的3.67g/cm3)，使其具有较高的本征探测效率；2.能量分辨率好(最好的可达到2.6％@662keV)，其能量分辨率已经接近了闪烁体泊松分布的统计极限，表明光产额的非线性对分辨率的影响要远小于其他的闪烁体；3.时间分辨率好(数百皮秒)；4.输出光谱与常用的光阴极相匹配。目前掺铈溴化镧已经实现商业化，如Saint-Gobain公司的BrilLanCe380型晶体、Ortec公司的digiDART-LF探测器等。溴化镧探测器基本组成结构如图1所示，溴化镧晶体(LaBr3:Ce)作为探测晶体，外围三面设置的反射层，另一面通过光耦合结构与光电倍增管相连，光电倍增管分别与高压分压器、以及前置放大器相连；输入高压通过高压分压器加载在光电倍增管，输出信号依次经过前置放大器、线性放大器和多道分析器的处理形成最终的输出信号。相比于其他的闪烁晶体，溴化镧晶体的主要缺点是其自身的放射性，其放射性的主要来源是晶体中天然存在的镧的放射性同位素138La(天然丰度0.09％)，半衰期1.02×1011年。138La的衰变类型有两种：轨道电子俘获和β-衰变。138La有66.4％的概率发生轨道电子俘获到138Ba的激发态，在这个过程中，伴随着特征X射线和俄歇电子发射，其中有63.7％的概率放出KαX射线，27.5％的概率放出LαX射线，8.8％的概率放出M层X射线。对于KαX射线，一个Ba的L层电子填充K层空穴可以释放出一个31.84keV的KαX射线，然后由于俄歇效应会产生5.6keV的俄歇电子，由于LaBr3:Ce的非正比响应，实际探测到的能量大约为35.5keV，同理LαX射线的能量约为4.5keV。之后再由激发态通过γ跃迁的形式衰变为稳定的138Ba并释放出1436keV的γ射线，此衰变产生的特征X射线和俄歇电子会与γ射线在溴化镧晶体内发生符合相加，从而形成一个1440keV和1472keV两个峰。此外，138La还有33.6％的概率发生β衰变到138Ce的激发态，β的最大能量为255.26keV，再以γ跃迁的形式衰变为稳定的138Ce并释放789keV的γ射线，其中的β粒子与γ射线的沉积能量也会发生符合相加，从而在800～1000keV形成一段连续谱。其自身本底谱如图2所示。除了138La外，在溴化镧晶体(LaBr3:Ce)中还存在着少量钾的放射性同位素钾40(40K)。40K是一种环境中常见的天然放射性核素，其天然丰度为0.0117％，半衰期是1.251×109年。40K的衰变类型有两种：轨道电子俘获和β-衰变，其中轨道电子俘获产生的激发态的40Ar进一步退激会释放出能量为1460.851keV的γ射线，发射强度为10.55％，这也是探测40K的主要指标。K在地壳中分布非常广泛，所以40K在环境中是常见的天然核素本底。例如，海水中的K含量大约为0.40g/kg，40K在海水中的天然活度大约在13Bq/L左右。然而，如果用溴化镧探测器来探测环境中的40K的活度，例如用溴化镧探测器对环境开展就地γ谱测量，由于138La自身的放射性，其在1440keV和1472keV各有一个峰，而40K的γ峰的能量是1461keV，位于138La的两个峰之间，于是在测量得到的溴化镧的能谱图上，这3个峰会互相重叠在一起，这就给测量环境中的40K的活度带来很大的影响，如图3是测量的40K的能谱，可以看出较为明显的多峰重叠现象，此外，溴化镧探测器中自身也存在着一定量的40K，虽然40K的活度不高，但是对于测量环境中的40K还是存在着一定的影响。由于探测器中138La和40K的存在，而环境中的40K活度通常都不高，通过直接扣除自身放射性的方法来得到40K的活度会带来很大的误差。由于这个缘故，通常都认为不能采用溴化镧探测器测量环境中的40K的活度。对于一定量的γ放射性核素，其活度可以通过不同的方式得到：采用理论计算的方法，即如式(1.1)A＝λN (1.1)其中，A为放射性核素的活度，λ是放射性核素的衰变常数，该常数与核素种类有关，为固定常数，N是放射性核素的核素数目，N可通过公式(1.2)计算得到：N＝m/M·NA (1.2)其中，m为核素质量，M为核素的原子质量，NA是阿伏伽德罗常数。除了理论计算之外，另一种确定核素活度的方法是采用γ谱测量的方法，将放射性物质放置在低本底γ谱仪中，通过事先确定的谱仪探测效率，结合测量的γ谱中的计数率，可以得出放射性核素的活度。通常，为了避免周围环境的影响，在γ谱测量过程中需建立较为完善的屏蔽环境。对于探测器自身放射性来说，也可以通过探测器自身的屏蔽测量得出放射性活度。对于γ探测器来说，刻度实验指的是通过测量标准刻度源，根据测量得到的刻度能谱，对探测器的性能进行刻度，主要包括能量刻度和能量分辨率刻度。从测量得到的γ谱中可以得到对应特征峰的能量和能量分辨率，从而建立起能量、能量分辨率和γ谱中道址之间的对应关系，从而实现探测器的能量刻度和能量分辨率刻度。针对某一能量Eγ的γ射线，通常重点关注的区域集中在能量在Eγ附近的区域，由于探测器响应存在一定的涨落，能量为Eγ的γ射线在探测器完全沉积能量后在能谱中显示的能量服从一定的分布形式，而不是单纯的Eγ。在扣除了其他因素的干扰后，这一分布可以用正态分布来描述，即能谱中γ射线的全能峰服从正态分布的形式，其均值为Eγ，标准差由探测器自身的能量分辨率决定。通常可以采用一个高斯函数来表示能谱中的全能峰，可以表示成：H·N(μ,σ)的形式，这里H代表着全能峰的面积，也即全能峰的总计数率，理论上可以表示为：H＝ε·A·Iγ (1.3)其中，ε是溴化镧探测器对γ射线的探测效率，A为目标核素的活度，Ir是目标核素发射γ射线的分支比。N(μ,σ)表示一个均值为μ，标准差为σ的标准正态分布，这里μ＝Eγ，σ由探测器自身的能量分辨率决定，在本专利技术中，统一用H·N(μ,σ)来表示能谱中的全能峰，不同的全能峰之间通过下标的形式进行区分。感兴趣区指的是在能谱处理过程中，需要重点关注或是分析的区域，通常γ谱分析过程中的感兴趣区都集中在全能峰附近，可以表示为：Eγ-k·σ～Eγ+k·σ，这里k可以根据实际情况进行选取，通常取k＝3。
技术实现思路
本专利技术针对掺铈溴化镧探测器高精度测量环境中40K的需求，提出了一种采用溴化镧探测器测量环境中的钾40(40K)的方法，该方法的核心在于通过高斯拟合的方式将测得的环境中的40K全能峰和溴化镧自身放射性造成的本底谱分离开来，得到环境中的40K的全能峰，从而有效的提高环境中的40K的测量精度，实现对环境中的40K的精确测量。本专利技术提出的一种采用溴化镧探测器测量环境中的钾40(40K)的方法，该方法包括以下步骤：1)对溴化镧探测器自身放射性进行标定，得到该探测器本底谱响应y；在138La的1440keV和1472keV全能峰感兴趣区中，溴化镧探测器的自身放射性形成3个γ全能峰，分别是138本文档来自技高网...
【技术保护点】
一种采用溴化镧探测器测量环境中的钾40的方法，其特征在于，该方法包括以下步骤：1)对溴化镧探测器自身放射性进行标定，得到该探测器本底谱响应y；在138La的1440keV和1472keV全能峰感兴趣区中，溴化镧探测器的自身放射性形成3个γ全能峰，分别是138La的1440keV、1472keV全能峰和溴化镧探测器自身含有的40K的1461keV全能峰，用三个高斯函数将感兴趣区内溴化镧探测器自身放射性的本底谱响应表示为：y＝H1·N(μ1，σ1)+H2·N(μ2，σ2)+H3·N(μ3，σ3) (1)其中，y代表溴化镧探测器自身放射性的本底谱响应；H1·N(μ1，σ1)、H3·N(μ3，σ3)为为138La的两个全能峰响应，H2·N(μ2，σ2)为40K的全能峰响应；H1、H2、H3分别为上述三个全能峰的总计数率；N(μ1，σ1)、N(μ2，σ2)、N(μ3，σ3)分别为上述三个全能峰的标准正态分布；σ1、σ2、σ3分别为上述三个全能峰的标准正态分布的标准差，μ1、μ2、μ3为均值；通过下述两种方法的之一，确定上述公式(1)中溴化镧探测器自身放射性谱响应中的各全能峰的高斯函数的参数，实现对探测器的标定，得到该探测器本底谱响应y；第一种方法具体包括：通过下述公式(1.1)、(1.2)计算或采用γ谱测量得到溴化镧探测器中138La的活度：A＝λN (1.1)其中，A为放射性核素的活度，λ是放射性核素的衰变常数，为固定常数，N是放射性核素的核素数目；N＝m/M·NA (1.2)其中，m为核素质量，M为核素的原子质量，NA是阿伏伽德罗常数；采用γ谱测量得到溴化镧探测器中40K的活度；通过摩卡模拟计算分别得到该探测器在138La的1440keV、1472keV全能峰以及40K的1461keV全能峰的γ射线探测效率ε，从核数据库中分别获取138La和40K的γ射线分支比Iγ，根据下述公式(1.3)分别计算获得138La的1440keV、1472keV全能峰以及40K的1461keV全能峰的总计数率H1、H2、H3；H＝ε·A·Iγ (1.3)其中，H为全能峰的总计数率；通过溴化镧探测器的刻度实验获取各高斯函数的标准差σ1、σ2、σ3；各高斯函数中的均值取理论值μ1＝1440，μ2＝1461，μ3＝1472；将得到的H1、H2、H3，σ1、σ2、σ3，μ1、μ2、μ3代入公式(1)中完成对该探测器自身放射性的标定；第二种方法具体包括：通过γ谱测量得到溴化镧探测器中40K的活度，并通过蒙卡模拟计算得到该探测器在40K的1461keV全能峰的γ射线探测效率ε，从核数据库中获取40K的γ射线分支比Iγ，通过上述公式(1.3)计算得到40K的总计数率H2，通过溴化镧探测器的刻度实验得到高斯函数的标准差σ2，得到40K的全能峰响应H2·N(μ2，σ2)；将溴化镧探测器放置在屏蔽环境中直接测量得到本底谱，结合40K的全能峰响应H2·N(μ2，σ2)，通过高斯拟合，确定138La的全能峰响应H1·N(μ1，σ1)和H3·N(μ3，σ3)；将上述得到的结果代入公式(1)中完成对该探测器自身放射性的标定；2)采用溴化镧探测器对环境中的40K进行测量，并拟合得到能谱响应y*；将溴化镧探测器放置在环境中进行测量，并使用三高斯拟合的方法将能谱中的响应表示为三个高斯函数的叠加：y*=H1*·N(μ1*,σ1*)+H2*·N(μ2*,σ2*)+H3*·N(μ3*,σ3*)---(2)]]>其中，y*代表溴化镧探测器在环境中感兴趣区内的能谱响应；为138La的两个全能峰响应，为40K的全能峰响应；分别为上述三个全能峰的总计数率；分别为上述三个全能峰的标准正态分布；分别为上述三个全能峰的标准正态分布的标准差，为均值；实测能谱中138La的贡献与本底谱中138La的贡献是一致的，也即：3)根据步骤1)得到的溴化镧探测器本底谱响应y和步骤2)得到的能谱响应y*，计算得到环境中40K活度：根据拟合得到的实测能谱响应y*和溴化镧探测器本底谱响应y，计算得到环境中40K对探测器总计数率的贡献为：通过蒙卡模拟计算得到溴化镧探测器对环境中40K的γ射线的探测效率εsp，从核数据库中得到40K的γ射线分支比为Iγ，从而根据下述公式得到环境中40K的活度A(K40)：A(K40)=(H2*-H2)/(ϵsp·Iγ)---(3).]]>...
【技术特征摘要】
1.一种采用溴化镧探测器测量环境中的钾40的方法，其特征在于，该方法包括以下步骤：1)对溴化镧探测器自身放射性进行标定，得到该探测器本底谱响应y；在138La的1440keV和1472keV全能峰感兴趣区中，溴化镧探测器的自身放射性形成3个γ全能峰，分别是138La的1440keV、1472keV全能峰和溴化镧探测器自身含有的40K的1461keV全能峰，用三个高斯函数将感兴趣区内溴化镧探测器自身放射性的本底谱响应表示为：y＝H1·N(μ1，σ1)+H2·N(μ2，σ2)+H3·N(μ3，σ3) (1)其中，y代表溴化镧探测器自身放射性的本底谱响应；H1·N(μ1，σ1)、H3·N(μ3，σ3)为为138La的两个全能峰响应，H2·N(μ2，σ2)为40K的全能峰响应；H1、H2、H3分别为上述三个全能峰的总计数率；N(μ1，σ1)、N(μ2，σ2)、N(μ3，σ3)分别为上述三个全能峰的标准正态分布；σ1、σ2、σ3分别为上述三个全能峰的标准正态分布的标准差，μ1、μ2、μ3为均值；通过下述两种方法的之一，确定上述公式(1)中溴化镧探测器自身放射性谱响应中的各全能峰的高斯函数的参数，实现对探测器的标定，得到该探测器本底谱响应y；第一种方法具体包括：通过下述公式(1.1)、(1.2)计算或采用γ谱测量得到溴化镧探测器中138La的活度：A＝λN (1.1)其中，A为放射性核素的活度，λ是放射性核素的衰变常数，为固定常数，N是放射性核素的核素数目；N＝m/M·NA (1.2)其中，m为核素质量，M为核素的原子质量，NA是阿伏伽德罗常数；采用γ谱测量得到溴化镧探测器中40K的活度；通过摩卡模拟计算分别得到该探测器在138La的1440keV、1472keV全能峰以及40K的1461keV全能峰的γ射线探测效率ε，从核数据库中分别获取138La和40K的γ射线分支比Iγ，根据下述公式(1.3)分别计算获得138La的1440keV、1472keV全能峰以及40K的1461keV全能峰的总计数率H1、H2、H3；H＝ε·A·Iγ (1.3)其中，H为全能峰的总计数率；通过溴化镧探测器的刻度实验获取各高斯函数的标准差σ1、σ2、σ3；各高斯函数中的均值取理论值μ1＝1440，μ2＝1461，μ3＝1472；将得到的H1、H2、H3，σ1、σ2、σ3，μ1、μ2、μ3代入公式(1)中完成对该探测器自身放射性的标定；第二种方法具体包括：通过γ谱测量得到溴化镧探测器中40K的活度，并通过蒙卡模拟计算得到该探测器在40K的1461keV全能峰的γ射线探测效率ε，从核数据库中获取40K的γ射线分支比Iγ，通过上述公式(1.3)计算得到40K的总计数率H2，通过溴化镧探测器的刻度实验得到高斯函数的标准差σ2，得到40K的全能峰响应H2·N(μ2，σ2)；将溴化镧探测器放置在屏蔽环境中直接测量得到本底谱，结合40K的全能峰响应H2·N(μ2，σ2)，通过高斯拟合，确定138La的全能峰响应H1·N(μ1，σ1)和H3·N(μ3，σ3)；将上述得到的结果代入公式(1)中完成对该探测器自身放射性的标定；2)采用溴化镧探测器对环境中...
【专利技术属性】
技术研发人员：曾志，潘兴宇，程建平，马豪，曾鸣，李君利，
申请(专利权)人：清华大学，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。