航空放射性信息提取方法初探

航空γ能谱探测技术因具备降低工作人员所受辐射剂量、强辐射测量效果媲美航线上可探测区域地面原位测量均值等优势被作为核应急常备的技术手段之一，在核事故污染范围圈定、核动力卫星放射性碎片与人工放射源定位中得到了广泛应用[1−4]。因人工放射性航空γ能谱刻度模型缺失，王南萍等[5]在室内环境下发现采用航空γ能谱中0.4−1.4MeV与1.4−3.0MeV间计数率之比VGC可有效揭示人工放射性是否存在。但Kock等[6]验证发现采用总道计数率变化规律来定位人工放射性热点的方法受天然放射性的分布不均影响。Grasty等[7]通过768个无人工放射性污染地区所测得天然放射性核素的航空γ能谱，建立单位含量eU、eTh或40K地层的航空γ谱仪响应谱，并结合三道法计算得到的当前测点eU、eTh和40K的含量，利用剥谱法提取铯响应谱，野外实验发现该方法能很好地揭示人工放射性位置信息。利用γ射线能谱数据估计放射性核素浓度的标准方法是加权最小二乘法拟合[8]，其核心为γ场叠加原理。随着探测器工艺及核电子学水平的提高，能谱仪的能量线性、稳定性等性能提升，利用蒙特卡罗(MonteCarlo,MC)方法模拟标准源响应谱库的加权最小二乘法拟合法已被广泛应用于能量色散X荧光谱[9]、水体原位γ能谱[10]和航空γ能谱[11−14]分析。2016年4月在核工业航测遥感中心举办了“提升核应急航空监测能力”中欧技术研讨会，会议双方确定将在使用蒙特卡罗方法进行人工放射性核素信息提取方面进行深入的研究。在应用最小二乘法解析航空γ能谱中，以往并未考虑本底（宇宙射线、仪器设备自身放射性和大气氡）响应谱的影响[11−13]，依据我国《航空γ能谱测量规范》[15]，用测区内大面积水域上同高度实测航空γ能谱替代航空γ能谱仪本底响应谱，经10km野外测线实验后验证其效果良好[14]。倘若测区内无类似大面积水域的情况下，如何获取航空γ能谱仪本底响应谱成为决定该方法应用效果的关键因素之一。为此本文拟建立一套航空γ能谱仪本底响应谱估计方法，提高基于最小二乘法航空γ能谱解析的普适性，并将其应用于提取人工放射性信息的实践，验证所建立方案的可靠性。

1宇宙射线响应谱校正模型

在远离海岸线的深海上空1800m及以上高空实测的航空γ能谱可认为仅由宇宙射线和仪器设备自身放射性的响应谱组成[15]。从理论上来说，航空γ谱仪对仪器设备自身放射性的响应谱与探测高度无关，可将其作为航空γ谱仪对宇宙射线响应谱剥离效果的比对依据。1997年Minty[16]认为可利用幂函数来描述H探测高度上航空γ谱仪对宇宙射线的响应谱中第x道计数CH,x：CH,x=BBH•Ex−1.3(1)式中：Ex为航空γ能谱第x道所对应的能量，刻度方法详见文献[17]；BBH为拟合系数，随探测高度H变化。为了避免天然放射性核素的影响，可采用2.85−3.03MeV范围内的实测谱分布来拟合BBH值。将相应探测高度上的实测航空γ能谱与式(1)反演得到的宇宙射线响应谱对应道计数相减，获得航空γ谱仪对仪器设备自身放射性响应谱随探测高度的相对变化规律（图1），可以看出湮灭辐射峰后各道计数基本重合，说明此能区内宇宙射线估计准确。但湮灭辐射峰前各道计数不重合，这是因为宇宙射线中的µ+介子等衰变产生了大量正电子并发生湮灭放出0.511MeV的γ射线，且该辐射对航空γ能谱低能区的影响不可忽略，其强度随探测高度而增大，说明宇宙射线产生的湮灭辐射[18]随穿透大气层厚度而改变。1.1宇宙射线中湮灭辐射峰强度校正模型在测区内一平坦区域让搭载航空γ谱仪的飞机平缓从地面起飞盘旋爬升至3000m高空，获得不同探测高度的实测航空γ能谱数据。考虑到航空γ能谱测量时间间隔1s，随探测高度升高，单位时间实测航空γ能谱受统计涨落影响更大，在数据处理时，以200m探测高度间隔统计实测航空γ能谱的平均谱。利用自适应峰形切削法扣除本底[19]、多高斯函数Levenberg-Marquardt算法拟合[20]获得湮灭辐射峰净面积NHT随探测高度H的变化规律如图2所示。采用谱线比法[21]获得地面Tl-208产生的湮灭辐射对NHT的贡献量NHG，两者相减则为宇宙射线对NHT的贡献量NHC。拟合得到NHC随探测高度H的变化规律如下（拟合优度R2=0.9742）：NHC=0.0121H−1.4635(2)图2湮灭辐射峰净面积组成Fig.2Compositionofthenetcountrateofannihilationradiation.由于飞行过程中高度难以稳定，难以运用式(2)进行野外校正。而宇宙射线越多，其产生湮灭光子也就越多，得到NHC与宇宙射线道计数NHL间的变化关系如图3所示，拟合方程（拟合优度R2=0.9672）如下：NHC=0.1083NHL−13.219(3)图3NHC与NHL间的关系Fig.3RelationsbetweenNHCandNHL.因单个航空γ能谱测量时间短，NHL亦受到统计涨落的影响，后续研究中将采用测区内高差5m范围内的航空γ能谱宇宙射线道计数平均值来表征。1.2宇宙射线中湮灭辐射响应谱MC模拟由于宇宙射线中µ+介子等衰变产生的湮灭辐射在空气介质中的质量衰减系数为0.0861075cm2•g−1[22]，则在密度为0.001293g•cm−3的空气中半衰减厚度为62.26m，说明有650m的空气可将湮灭光子几乎完全衰减掉。在采用MC模拟航空γ谱仪对宇宙射线中湮灭光子的响应谱时，应将航空γ谱仪放置在圆柱体（直径与高均为650m）空气介质中。为减小模拟空间体积，通过介质互换原理，将空气密度提升100倍，此时圆柱体尺寸可减小至原来的1/100，模型图如图4所示。采用MC模拟软件GEANT4编写上述模型，源粒子强度分布按式(2)设置，抽样总数设置为4×1011个，模拟结果如图5所示，不确定度为0.78%。图4宇宙射线中湮灭辐射响应谱MC模拟模型Fig.4MCsimulationmodelofresponsespectrumirradiatedbyannihilationradiationincosmicrays.图5宇宙射线中湮灭辐射的MC模拟响应谱Fig.5MCsimulationresponsespectrumirradiatedbyannihilationradiationincosmicrays.1.3校正效果分析运用式(3)及图5响应谱反演获得图1中6个探测高度上航空γ谱仪对宇宙射线响应谱，并从图1中去除后结果如图6所示。可以看出，此时不同探测高度上航空γ谱仪对仪器设备自带放射性的响应谱基本吻合，与理论规律相符，证实了上述方法的有效性。后续研究中将图6中所有谱线的平均谱作为航空γ谱仪对仪器设备自带放射性的响应谱。图6不同探测高度仪器设备自带放射性本底谱再估计Fig.6Estimatebackgroundofgammaradiationfrominstrumentationbythemethodinthispaperatdifferentheights.

2大气氡子体响应谱的近似替代

在航空γ谱仪对平衡天然铀系地层的响应谱MC模拟时，输入源项为平衡天然铀系中每百次衰变产生量大于1的特征γ射线[13−14,23]，源自234Th、226Ra、214Pb、214Bi和210Pb这5种放射性核素。前两种放射性核素产生的特征γ射线最大能量为186.211keV，说明仅对航空γ谱仪响应谱中低能谱段有贡献，影响人工放射性如214Am的定量精度；仅占模拟特征γ射线源粒子数不到7%，说明影响量可近似忽略。后三者则为大气氡子体，说明在MC模拟时源抽样粒子能谱分布近似相同。以下详细探讨采用“航空γ谱仪对平衡天然铀系地层的响应谱”替代“航空γ谱仪对大气氡子体的响应谱”的可行性。结合上述分析，可近似采用式(4)表征内陆大面积湖泊H探测高度上实测航空γ能谱第x道计数yH,x：yH,x=Sx+CH,x+NHC•Dx+cU•UH,x+cTh•TH,x+cK•KH,x+εx(4)式中：Sx为仪器设备自带放射性对航空γ能谱第x道计数的贡献；CH,x+NHC•Dx为宇宙射线响应谱中第x道计数；UH,x、TH,x和KH,x分别表示MC模拟得到仅含平衡235U&238U系、平衡钍系和40K的地层上H探测高度上航空γ谱仪响应谱（特征峰区最大计数归一化）第x道的计数（详见文献[14]）；cU、cTh和cK为待拟合参数；εx为航空γ能谱第x道实测计数与上述各组成总计数率间的差值；其余符号同前所述。拟合代码采用MINUIT软件包[24]。对水库内陆上空8个探测高度（60m、90m、120m、150m、180m、210m、240m、270m和300m）上实测航空γ能谱（合计60s累积测量谱）进行拟合，得到拟合参数cU、cTh和cK的结果如图7所示，典型全谱拟合结果如图8所示，8个探测高度全谱拟合相对偏差均在±3.63%以内。图7内陆水库上拟合参数值随高度变化规律Fig.7Verticaldistributionoffittingvaluesabovethewaterofareservoirs.图8内陆水库上210m高空航空γ能谱拟合效果Fig.8Fittingeffectofairbornegamma-rayspectrumabovethewaterofareservoirswhenH=210m.从拟合结果可以看出：1)水域上空40K的贡献为0，说明40K发射的特征γ射线被空气衰减几乎殆尽，影响可基本忽略。2)各道计数的主要贡献来源于大气氡子体所发射的γ射线，数据变化规律与文献[21]的对数增长规律类似，说明这部分计数主要为大气氡子体的贡献。3)结果中显示存在少量钍系特征γ射线的贡献，且cTh拟合值成微弱增大趋势，可能源自水域周围地层中的钍系特征γ射线（能量较高的如2.62MeV，穿透能力更强）、空气中220Rn子体等的贡献。综上所述，采用“航空γ谱仪对平衡天然铀系地层的响应谱”替代“航空γ谱仪对大气氡子体的响应谱”是可行的。

3人工辐射环境下实验验证及效果分析

3.1活度反演将137Cs（活度2.588×109Bq）和60Co（活度1×108Bq）点源摆放在单箱晶体底面中心，用氢气艇将单箱晶体升至不同高度（范围5−85m、间隔10m）测得不同探测高度下的实测航空γ能谱，并用式(5)拟合：yH,x=Sx+CH,x+NHC•Dx+cU•UH,x+cTh•TH,x+cK•KH,x+∑cj•IjH,x+εx(5)式中：,jHxI为第j种人工放射性核素面源上H探测高度处实测航空γ能谱中第x道计数；cj为对应待拟合参数；N为人工放射性核素的种类。从图9的拟合效果来看，0.2−2.285MeV（17−192道）间各道计数拟合效果俱佳。利用自适应峰形切削法扣除本底[19]、多高斯函数Levenberg-Marquardt算法拟合[20]获得的特征峰净面积计算得到不同探测高度H下的探测效率值，结合cj拟合值，反推出各探测高度H下地面点源活度与计量标准值间的相对偏差如表1所示，可以看出在探测限[23]范围内吻合度在±15%以内，说明拟合方案正确可靠。3.2辐射热点定位地面查证确定有一废弃137Cs源被埋在地面以下10cm处，对离源中心垂直距离为20m处已进行的580s航空γ能谱测量（飞行速度为70km•h−1、飞行高度为75m，单谱测量时间为1s，137Cs源横坐标为298s）数据进行分析，得到总道计数率(0.4−3MeV)IT分布曲线如图10上部实线所示；为减小统计涨落的影响，采用聚类分析的NASVD方法[25]对航空γ能谱平滑（平滑后航空γ能谱第x道计数为xHy,），平滑能谱的总道计数率TI如图10上部点状线所示；利用xHy,代替式(5)中yH,x，得到拟合谱的总道计数率TI′变化规律。可以看出，从上述三者的变化规律中无法获取人工源的位置信息。平滑后航空γ能谱中VGC的变化规律如图10中部虚线所示。计算发现VGC的平均值为4.002、方差为0.270，可以看出大于3倍标准差的异常值横坐标范围分别为34−36s和296−301s，说明亦受天然放射性影响难以有效定位人工放射性热点的位置。将拟合谱中宇宙射线、仪器设备自身放射性、大气氡子体、40K、铀系和钍系的响应谱剔除，剩余谱中总道计数率TI′变化如图10下部实线所示。计算得到平均值为68.151、方差为60.701，大于3倍标准差异常值横坐标范围为297−300s，中心位于298s处，说明本文方法能很好地定位人工放射性热点的位置。

4结语

本文依据实验数据，获得航空γ能谱中宇宙射线内湮灭辐射强度与宇宙射线道计数间的线性规律，并MC模拟获得航空γ谱仪对宇宙射线中湮灭辐射的响应谱，结合Minty的幂函数估计方法，扣除远离海岸线深海面上2100−3600m的实测航空γ能谱中宇宙射线响应谱，显示航空γ谱仪对仪器设备自身放射性的响应谱与探测高度无关，与理论规律一致，说明上述方法是可靠的。同时利用基于最小二乘的航空γ能谱解谱方法对内陆水库上8个探测高度上实测航空γ能谱进行全谱拟合，发现全谱拟合值与实测值间相对偏差在±3.63%以内符合，cU值随高度变化规律与文献[21]氡剥离规律中相同，证实可运用“航空γ谱仪对平衡天然铀系地层的响应谱”替代“航空γ谱仪对大气氡子体的响应谱”。最终将上述结论应用于点源上空不同高度的实测航空γ能谱解析，实验发现在探测限范围内反演活度与计量标称值在±15%以内吻合；并将其应用于人工放射性137Cs源定位，验证发现该方法能有效剔除天然辐射成分的影响、精确定位人工辐射热点位置。致谢感谢格拉斯哥大学AllysonJD博士提供的航空γ能谱测量数据及在研究过程中给予的指导。

作者:杨悦 吴和喜 刘玉娟 刘义保 孟凡松 张思颖 单位:1.东华理工大学 2.东华理工大学