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

　　航空γ能谱探测技术因具备降低工作人员所受辐射剂量、强辐射测量效果媲美航线上可探测区域地面原位测量均值等优势被作为核应急常备的技术手段之一，在核事故污染范围圈定、核动力卫星放射性碎片与人工放射源定位中得到了广泛应用[1?4]。 　　因人工放射性航空γ能谱刻度模型缺失，王南萍等[5]在室内环境下发现采用航空γ能谱中04?14与14?30间计数率之比可有效揭示人工放射性是否存在。 　　但等[6]验证发现采用总道计数率变化规律来定位人工放射性热点的方法受天然放射性的分布不均影响。 　　等[7]通过768个无人工放射性污染地区所测得天然放射性核素的航空γ能谱，建立单位含量、或40地层的航空γ谱仪响应谱，并结合三道法计算得到的当前测点、和40的含量，利用剥谱法提取铯响应谱，野外实验发现该方法能很好地揭示人工放射性位置信息。 　　利用γ射线能谱数据估计放射性核素浓度的标准方法是加权最小二乘法拟合[8]，其核心为γ场叠加原理。 　　随着探测器工艺及核电子学水平的提高，能谱仪的能量线性、稳定性等性能提升，利用蒙特卡罗,方法模拟标准源响应谱库的加权最小二乘法拟合法已被广泛应用于能量色散荧光谱[9]、水体原位γ能谱[10]和航空γ能谱[11?14]分析。 　　2016年4月在核工业航测遥感中心举办了提升核应急航空监测能力中欧技术研讨会，会议双方确定将在使用蒙特卡罗方法进行人工放射性核素信息提取方面进行深入的研究。 　　在应用最小二乘法解析航空γ能谱中，以往并未考虑本底宇宙射线、仪器设备自身放射性和大气氡响应谱的影响[11?13]，依据我国《航空γ能谱测量规范》[15]，用测区内大面积水域上同高度实测航空γ能谱替代航空γ能谱仪本底响应谱，经10野外测线实验后验证其效果良好[14]。 　　倘若测区内无类似大面积水域的情况下，如何获取航空γ能谱仪本底响应谱成为决定该方法应用效果的关键因素之一。 　　为此本文拟建立一套航空γ能谱仪本底响应谱估计方法，提高基于最小二乘法航空γ能谱解析的普适性，并将其应用于提取人工放射性信息的实践，验证所建立方案的可靠性。 　　1宇宙射线响应谱校正模型在远离海岸线的深海上空1800及以上高空实测的航空γ能谱可认为仅由宇宙射线和仪器设备自身放射性的响应谱组成[15]。 　　从理论上来说，航空γ谱仪对仪器设备自身放射性的响应谱与探测高度无关，可将其作为航空γ谱仪对宇宙射线响应谱剥离效果的比对依据。 　　1997年[16]认为可利用幂函数来描述探测高度上航空γ谱仪对宇宙射线的响应谱中第道计数,,=??131式中为航空γ能谱第道所对应的能量，刻度方法详见文献[17]；为拟合系数，随探测高度变化。 　　为了避免天然放射性核素的影响，可采用285?303范围内的实测谱分布来拟合值。 　　将相应探测高度上的实测航空γ能谱与式1反演得到的宇宙射线响应谱对应道计数相减，获得航空γ谱仪对仪器设备自身放射性响应谱随探测高度的相对变化规律图1，可以看出湮灭辐射峰后各道计数基本重合，说明此能区内宇宙射线估计准确。 　　但湮灭辐射峰前各道计数不重合，这是因为宇宙射线中的μ+介子等衰变产生了大量正电子并发生湮灭放出0511的γ射线，且该辐射对航空γ能谱低能区的影响不可忽略，其强度随探测高度而增大，说明宇宙射线产生的湮灭辐射[18]随穿透大气层厚度而改变。 　　11宇宙射线中湮灭辐射峰强度校正模型在测区内一平坦区域让搭载航空γ谱仪的飞机平缓从地面起飞盘旋爬升至3000高空，获得不同探测高度的实测航空γ能谱数据。 　　考虑到航空γ能谱测量时间间隔1，随探测高度升高，单位时间实测航空γ能谱受统计涨落影响更大，在数据处理时，以200探测高度间隔统计实测航空γ能谱的平均谱。 　　利用自适应峰形切削法扣除本底[19]、多高斯函数-算法拟合[20]获得湮灭辐射峰净面积随探测高度的变化规律如图2所示。 　　采用谱线比法[21]获得地面-208产生的湮灭辐射对的贡献量，两者相减则为宇宙射线对的贡献量。 　　拟合得到随探测高度的变化规律如下拟合优度2=09742=00121?146352图2湮灭辐射峰净面积组成2由于飞行过程中高度难以稳定，难以运用式2进行野外校正。 　　而宇宙射线越多，其产生湮灭光子也就越多，得到与宇宙射线道计数间的变化关系如图3所示，拟合方程拟合优度2=09672如下=01083?132193图3与间的关系3因单个航空γ能谱测量时间短，亦受到统计涨落的影响，后续研究中将采用测区内高差5范围内的航空γ能谱宇宙射线道计数平均值来表征。 　　12宇宙射线中湮灭辐射响应谱模拟由于宇宙射线中μ+介子等衰变产生的湮灭辐射在空气介质中的质量衰减系数为008610752??1[22]，则在密度为0001293??3的空气中半