作者简介: 徐立鹏, 1987年生, 四川省辐射环境管理监测中心站工程师 e-mail: 836197886@qq.com
为了实现车载γ谱仪巡测系统对辐射剂量率的准确测定, 提出基于快速傅里叶变换(FFT)本底扣除法的改进型BP神经网络模型(FFT-BP神经网络模型)。 实验采用γ射线能谱分析法, 对不同间距处的137Cs放射源进行车载γ能谱测试, 将得到的谱数据通过快速傅里叶变换(FFT)扣除本底, 获得新的谱线数据。 应用FFT-BP神经网络模型对未知剂量的车载γ谱线作辐射剂量率的定量预测, 将预测结果同3个函数模型的拟合结果比较, 验证FFT-BP神经网络模型的预测效果。 结果表明, FFT扣除法能较好的削弱散射本底对γ谱线的影响, 能有效的降低谱线本底。 通过新谱线获得的特征峰面积和净谱线面积与辐射剂量率的相关系数均为0.99( p<0.05), 相关性显著。 模型拟合分析过程中, FFT-BP神经网络模型表现出较强的学习泛化能力, 预测较理想, 相对误差和累计误差分别低于0.6%和9%, 效果明显优于数学模型和γ能谱全能峰法, 可显著降低γ能谱分析辐射剂量率的误差, 且能有效提升工作效率。 因此, FFT-BP神经网络模型适用于γ能谱辐射剂量的预测分析, 为车载γ谱仪巡测系统测量辐射剂量提供了一种新型有效的分析方法。
In order to measure the radiation dose rate accurately with carborne γ spectrometer patrol system, proposed a modified back-propagation network(BP network) model basised on fast fourier transform background deduct method(FFT-BP network model). Using γ-ray energy spectrum analysis method to test the carborne γ-ray energe spectrum of Cs-137 of different spacings, adopting FFT method to deduct the background of spectrum data then get new spectrum data. The modified B-P network model is applied to qualitatively predict the radiation dose rate of unknow dose carborne γ spectrum, by comparing the predicted results with fitting results of 3 function models to verify the effect of FFT-BP network model. The results show the FFT deduct method can weaken the influence of the scattering background on γ spectrum and reduce spectrum background effectively. The correlation coefficients between characteristic peak area and net area getting from new spectrum are 0.99 ( p<0.05), which shows a remarkable correlation. In the process of model fitting, FFT-BP network model shows strong ability of learning and generalization, the prediction of experimental results is ideal, relative error and accumulative error are below 0.6% and 9% respectively, it has better effect than mathematical methods and gamma spectra method and it also can reduce the error of radiation dose rate analized by γ-spectra analysis method, improve the work efficiency effectively. There fore, FFT-BP network model can apply to predictive analysis of γ-ray energy spectrum radiation dose, which provide a new and efficient method for carborne γ spectrometer patrol system to measure radiation dose.
引 言
车载γ 能谱巡测仪因具有快速巡测、 搜寻面广、 灵敏度高等特点, 已在环境监测的例行巡测、 寻找放射源和核应急等领域崭露头角[1, 2]。 目前, 国外市场的主流NaI(Tl)车载γ 能谱仪众多, 国内已有数家单位在研发相关产品[3, 4]。 但利用NaI(Tl)车载γ 能谱仪对环境和放射性物质进行辐射剂量的准确监测至今仍是一个难题, 我国也尚未发布车载γ 能谱仪剂量标定的相关标准和技术规范。
由于入射光子能量的较大差异, 采用核脉冲计数率的辐射剂量监测设备均存在严重的能量响应问题。 为此, 国内外学者提出采用γ 能谱法换算辐射剂量的方法, 包括总能谱法或总窗法、 贝壳公式法[5]、 HASL标定法和G(E)函数法[6]等分析法, 但因存在噪声影响、 刻度繁琐、 限制因素多和误差较大等缺点, 难以准确分析谱线代表的辐射剂量率。
利用快速傅里叶变换(FFT)扣除法削弱谱线散射本底的优势, 结合BP神经网络在预测预报及评价方面的优点, 研究FFT-BP神经网络模型分析预测车载γ 能谱测量辐射剂量率的效果。 为实现通过γ 能谱对辐射剂量率的精准监测提供重要的技术支持。
使用的车载γ 谱仪巡测系统(简称车载系统)是由成都理工大学研制的基于数字多道技术的大体积NaI(Tl)γ 能谱仪。 系统内的数字梯形通道[7]可准确获取高脉冲峰, 避免高计数率下的数据堆积。
为降低环境对实验结果的影响, 要求在测试环境不变的情况下, 分别使用FH40G剂量仪和车载系统对不同间距相同点位处的137Cs放射点源进行测试, 获得对应间距处的辐射剂量率和γ 能谱数据。 其中FH40G和放射性豁免点源137Cs均通过中国测试研究科学院计量检定, 均在使用有效期内。
对于核事故应急的特殊要求, 车载NaI(Tl) γ 能谱巡测系统需实现快速监测。 实验设计在测量周期为1 s的条件下对137Cs源进行测试。 应用快速傅里叶变换法(FFT)[8]对车载γ 谱线作本底扣除处理, 利用改进型的BP神经网络模型学习训练出输入变量和辐射剂量率之间的非线性关系, 通过对预测样本的准确分析, 验证车载谱仪巡测系统可实现对辐射剂量率的准确测量, 分析FFT-BP神经网络模型的应用效果。
FFT-BP神经网络模型是快速傅里叶变换(FFT)扣除法与误差反向传播网络相结合的一种数据融合模型。 该模型是对车载系统快速巡测的谱线数据作FFT本底扣除处理, 在无需描述映射关系数学函数的情况下, 将原始谱线、 净谱线和特征峰面积等作为输入变量, 利用改进型BP神经网络, 建立以设定精度逼近任何非线性连续函数。
在实际的监测过程中, 为获知核素的活度需对γ 谱线扣除本底[9, 10, 11]。 为了更准确地研究车载系统的探测效果, 实验采用快速傅里叶变换法(FFT)对系统所测谱线进行本底扣除处理。 离散傅里叶变换的基本公式[12]为
车载系统所测谱线包括高频的全能峰和低频的散射本底。 傅里叶变换本底扣除法[13]可实现利用理想滤波器扣除高频部分, 留下低频部分。
BP神经网络(back-propagation network, BP) 是一种将误差逆向传播并训练的多层前馈网络算法。 该神经网络的基本结构包括输入层、 隐含层和输出层, 同层间无关联, 数据信息流前向传递, 通过反向传播误差信息不断调整权重系数与偏置量, 直到误差满足要求, 网络训练结束, 网络结构如图1所示。 BP神经网络通过不断学习、 训练和归纳, 能实现不可见非线性映射关系的存贮, 可以处理车载巡测γ 谱数据与辐射剂量之间的复杂关系。 实验利用Matlab软件编写程序, 建立改进型Levenberg-Marquardt BP神经网络。
FFT扣除法是一种基于快速傅里叶变换法(FFT)扣除谱线本底, 使用面积积分法计算峰面积, 再应用拉依达准则剔除异常数据, 得到净谱线或特征峰净面积的分析方法。 图2是利用FFT扣除法处理137Cs源和60Co源共同作用的谱线效果图。
为了比较FFT扣除法的效果, 实验利用本底谱扣除法, 即借助无源本底谱线对整个谱线数据本底或特征峰本底进行扣除[5]的分析方法。 由图3可见, 本底谱扣除法仅仅对低能区的本底进行了有限的扣除, 对特征峰的作用效果并不明显。 而图2中的FFT扣除法对谱线本底的扣除相较于本底谱扣除法更明显、 更彻底。 对比发现FFT扣除法能更好的扣除低能区和特征峰的本底, 特别是对特征峰散射本底的扣除效果更佳。
研究表明, 利用γ 能谱全能峰法测量辐射剂量率的效果显著, 其中测试137Cs, 60Co, 133Ba, 152Eu, 241Am和226Ra等核素的相对误差均小于7%[14]。 为验证FFT-BP神经网络模型分析预测辐射剂量率的效果, 仅设计车载系统和FH40G对137Cs源进行不同间距处的同步测试实验, 获得123组实验数据。 由于实验所使用的FH40G是通过137Cs源校准的, 故其测量137Cs实验源较精准, 结果见表1所示。 表中661 keV峰和净谱线分别表示通过FFT扣除法处理得到的137Cs特征峰和新谱线计数均值, 测量值为FH40G测量得到的辐射剂量率均值。 分析发现, 特征峰面积和净谱线与辐射剂量率均存在正相关性, 相关系数均为0.99 (p< 0.05), 相关性显著。 间距与辐射剂量率呈负相关性, 相关系数为0.86 (p< 0.05), 相关性显著。
应用FFT扣除法对车载系统测试获取137Cs源的谱线进行处理, 得到谱线的净计数和661 keV特征峰的净峰面积计数, 从而将复杂的γ 谱线转换为定性且定量的输入变量。 因此, 实验设定间距、 谱线原始计数、 谱线净计数和特征峰净面积等相关项为改进型BP神经网络的输入变量, 输入节点为4; 输出变量为辐射剂量率, 即输出节点为1。 在网络训练前, 设置允许误差为0.001, 最大迭代次数为1 000次, 实现以目标精度逼近任何非线性连续函数。
FFT-BP神经网络模型对108组实验数据进行网络学习, 再将15组预测样本代入训练好的预测模型, 计算预测结果。 图4为137Cs源的神经网络回归图, 其归一化系数R2几乎均为1.00, 说明该网络预测模型的建立以及对137Cs源辐射剂量率的预测均比较成功。
图5反映了BP神经网络对137Cs源辐射剂量率的预测效果, 横坐标表示15个样本数, 纵坐标分别为辐射剂量率和相对偏差。 由图5(a)可知BP output(预测值)点和Real value(真值)点几乎重合, 且二者的相对误差均在± 0.6%范围以内, 如图5(b)所示。
实验结果表明该模型的学习泛化能力很强, 能比较真实地体现输入变量与辐射剂量率之间的复杂关系, 适用于车载系统辐射剂量值的预测分析。
为了验证FFT-BP神经网络模型在车载系统中预测辐射剂量率的应用效果, 使用最小二乘法对影响辐射剂量率的变量进行函数拟合, 建立相关数学模型。 分析发现, 通过FFT扣除法处理得到的137Cs源的特征峰净面积与净谱线的相关系数为0.996 (p< 0.05), 线性拟合度为0.99, 表明二者的线性度很高。 为简化函数, 实验选取间距(x1)或净谱线(x2)与辐射剂量率分别建立平方倒数数学模型(y1)和线性数学模型(y2), 函数关系分别见式(3)和式(4)。
由于辐射剂量率受间距和放射性核素活度的共同影响, 而核素活度可通过γ 谱线间接表示。 为了较准确的分析辐射剂量率与间距和谱线之间的关系, 使用Matlab软件建立辐射剂量率随间距(x1)和净谱线(x2)变化的数学模型, 即间距-净谱-剂量率模型(y3), 函数关系见式(5)。
对三种数学模型的函数进行拟合, 结果见表2。 比较发现间距-净谱-剂量率模型的拟合系数最高, 为0.997, 说明该模型的拟合效果优于其他两种数学模型。
将相同的15组γ 谱线测验样本分别代入4中模型中, 比较表2中的结果发现, FFT-BP神经网络模型的相对误差值最小, 仅为0.58%, 优于全能峰法的相对误差1.12%[14]; 累计误差亦最小, 为8.58%。 说明该FFT-BP神经网络模型对车载γ 谱仪辐射剂量率的分析预测准确度相当高, 可实现对车载γ 谱线与辐射剂量率之间复杂关系的校正。 同时, 由于FFT-BP神经网络模型可快速、 准确获知γ 谱线表示的辐射剂量率信息, 相较于拟合函数, 缩简了分析复杂谱线的工作量, 节省了建模时间, 有效提升工作效率。
γ 谱线既能定性核素类型, 又能通过复杂的能谱分析定量核素活度, 但是难以利用简单的数学模型表示谱线和辐射剂量之间的复杂关系, 获得比较准确的分析结果。 结合FFT扣除法和改进型BP神经网络的特点提出了FFT-BP神经网络模型, 该模型能较理想的解决上述难题。 实验结果表明, FFT-BP神经网络模型对137Cs源辐射剂量率的分析预测准确度较高, 相对误差小于0.6%, 累计误差低于9%, 分析预测效果明显优于数学模型和γ 能谱全能峰法。 说明该神经网络模型分析预测车载系统辐射剂量率的准确度较好, 基本克服了复杂车载谱线难以准确预测辐射剂量的难题, 实现了对车载γ 谱线与辐射剂量之间非线性关系的校正, 有效提升了工作效率。 所以FFT-BP神经网络模型适用于车载γ 谱仪巡测系统测量辐射剂量的分析预测, 可为NaI(Tl)车载γ 能谱仪的剂量监测工作提供重要的技术支持。
The authors have declared that no competing interests exist.