航空放射性信息提取方法初探论文.docVIP

航空放射性信息提取方法初探 航空γ能谱探测技术因具备降低工作人员所受辐射剂量、强辐射测量效果媲美航线上可探测区域地面原位测量均值等优势被作为核应急常备的技术手段之一，在核事故污染范围圈定、核动力卫星放射性碎片与人工放射源定位中得到了广泛应用[1#8722;4]。因人工放射性航空γ能谱刻度模型缺失，王南萍等[5]在室内环境下发现采用航空γ能谱中0.4#8722;1.4MeV与1.4#8722;3.0MeV间计数率之比VGC可有效揭示人工放射性是否存在。但Kock等[6]验证发现采用总道计数率变化规律来定位人工放射性热点的方法受天然放射性的分布不均影响。Grasty等[7]通过768个无人工放射性污染地区所测得天然放射性核素的航空γ能谱，建立单位含量eU、eTh或40K地层的航空γ谱仪响应谱，并结合三道法计算得到的当前测点eU、eTh和40K的含量，利用剥谱法提取铯响应谱，野外实验发现该方法能很好地揭示人工放射性位置信息。利用γ射线能谱数据估计放射性核素浓度的标准方法是加权最小二乘法拟合[8]，其核心为γ场叠加原理。随着探测器工艺及核电子学水平的提高，能谱仪的能量线性、稳定性等性能提升，利用蒙特卡罗(MonteCarlo,MC)方法模拟标准源响应谱库的加权最小二乘法拟合法已被广泛应用于能量色散X荧光谱[9]、水体原位γ能谱[10]和航空γ能谱[11#8722;14]分析。2016年4月在核工业航测遥感中心举办了“提升核应急航空监测能力”中欧技术研讨会，会议双方确定将在使用蒙特卡罗方法进行人工放射性核素信息提取方面进行深入的研究。在应用最小二乘法解析航空γ能谱中，以往并未考虑本底（宇宙射线、仪器设备自身放射性和大气氡）响应谱的影响[11#8722;13]，依据我国《航空γ能谱测量规范》[15]，用测区内大面积水域上同高度实测航空γ能谱替代航空γ能谱仪本底响应谱，经10km野外测线实验后验证其效果良好[14]。倘若测区内无类似大面积水域的情况下，如何获取航空γ能谱仪本底响应谱成为决定该方法应用效果的关键因素之一。为此本文拟建立一套航空γ能谱仪本底响应谱估计方法，提高基于最小二乘法航空γ能谱解析的普适性，并将其应用于提取人工放射性信息的实践，验证所建立方案的可靠性。 1宇宙射线响应谱校正模型 在远离海岸线的深海上空1800m及以上高空实测的航空γ能谱可认为仅由宇宙射线和仪器设备自身放射性的响应谱组成[15]。从理论上来说，航空γ谱仪对仪器设备自身放射性的响应谱与探测高度无关，可将其作为航空γ谱仪对宇宙射线响应谱剥离效果的比对依据。1997年Minty[16]认为可利用幂函数来描述H探测高度上航空γ谱仪对宇宙射线的响应谱中第x道计数CH,x：CH,x=BBH#8226;Ex#8722;1.3(1)式中：Ex为航空γ能谱第x道所对应的能量，刻度方法详见文献[17]；BBH为拟合系数，随探测高度H变化。为了避免天然放射性核素的影响，可采用2.85#8722;3.03MeV范围内的实测谱分布来拟合BBH值。将相应探测高度上的实测航空γ能谱与式(1)反演得到的宇宙射线响应谱对应道计数相减，获得航空γ谱仪对仪器设备自身放射性响应谱随探测高度的相对变化规律（图1），可以看出湮灭辐射峰后各道计数基本重合，说明此能区内宇宙射线估计准确。但湮灭辐射峰前各道计数不重合，这是因为宇宙射线中的micro;+介子等衰变产生了大量正电子并发生湮灭放出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#8722;1.4635(2)图2湮灭辐射峰净面积组成Fig.2Compositionofthenetcountrateofannihilationradiation.由于飞行过程中高度难以稳定，难以运用式(2)进行野外校正。而宇宙射线越多，其产生湮灭光子也就越多，得到NHC与宇宙射线道计数NHL