多球或嵌套(3He、金箔)能谱仪(Bonner Sphere or Multi-foil Spectrometer)是常用的中子能谱测量方法,而且是唯一可覆盖热中子至高能中子的探测方法。然而,多球中子谱仪或多箔中子谱仪法的测量结果需要通过专门的软件进行解谱或“能谱调整”,才能得到适当的中子能谱结果。
中子能谱仪解谱软件NSUP实现了离线或实时在线多球或嵌套中子能谱仪的可视化解谱功能,该软件采用多种解谱算法,其中,探测器读数、预置谱、响应矩阵数据等均可通过SQL数据库进行传递,软件操作简便、使用方便。
其中,SPUNIT算法采用迭代法进行解谱,也支持“1/E谱 + 麦克斯韦谱”搜索。最大似然期望最大化(MLEM)算法也属于迭代算法。
NSUP集成了IAEA发布的NMF-90解谱软件包中的三个解谱程序,即MSANDB、MIEKEB、MSITER,它们采用的方法分别为迭代方法、蒙卡方法和最小二乘法。其中MIEKEB和MSITER解谱算法最好提供协方差数据,这样才能获得更准确的解谱结果。
NSUP可用于BNCT照射装置中子能谱解谱。后续将集成更多解谱算法。
图1 D-T中子源解谱(迭代法,初始谱为D-T聚变谱)
图 2 D-T中子源解谱(谱搜索,不使用初始谱)
图3 重水慢化锎源解谱(谱搜索,不使用初始谱)
图4 重水慢化锎源解谱(迭代法,初始谱为重水慢化锎源,数据库输入解谱文件)
图5 NSUP支持通过文件导入进行解谱(由数据库解谱时可以导出此文件)
图6 重水慢化锎源解谱(迭代法,初始谱为重水慢化锎源能谱,文件载入解谱数据)
已完成最大似然期望最大化(MLEM)算法的集成,后续将集成更多解谱算法。
图7 锎中子源裸源解谱(MLEM算法,数据库输入解谱数据)
已完成MSANDB算法(SAND-II)的集成。
图8 镅铍源解谱(MSANDB算法,H3管嵌套中子谱仪)
已完成MIEKEB算法(蒙卡)、MSITER算法(最小二乘法)的集成。
图9 MIEKEB算法(蒙卡方法)(未提供协方差数据)
图10 MSITER算法(最小二乘法)(未提供协方差数据)
中子能谱解谱软件NSUP目前已集成6种解谱算法,每种算法各有特点,迭代法简单易用,蒙卡方法和最小二乘方法需要配合中子源参数的协方差数据(也可以不提供协方差数据)。当然,解谱程序只是个经验工具,重点还是要对中子源做细致的测算和分析。
图11 NSUP已集成6种中子能谱解谱算法
进一步信息请参阅“NSUP中子能谱解谱程序使用手册 v1.0”。