基于概率迭代的NDP反演方法.pdfVIP

第29卷第6期 计 算 物 理 V01．29，No．6 2012年11月 CHINESEJOURNALOFCOMPUTATIONAL PHYSICS NOV，2012 文章编号：1001—246X(2012)06_0891．10 基于概率迭代的NDP反演方法 杨 鑫1’2， 李润东1， 刘汉刚1， 王冠博1”， 王 侃2 (1．核物理与化学研究所，四川绵阳621900；2．清华大学工程物理系，北京100084) 摘 要：针对中子深度分析(NDP)技术，基于概率迭代思想，推导一种迭代反演方法，对其和线性正则化方法在 NDP反演问题中的应用进行比较．未引入计数随机误差时，两者都能得到较好的反演结果，迭代法的结果是非负 的，线性正则化方法在源强度发生阶跃处得到的结果更好．在考虑随机误差存在的情况下，迭代法的效果更佳．同 时，研究了迭代矩阵对结果的影响，矩阵元的指数放大可以实现迭代过程的加速． 关键词：中子深度分析(NDP)；反演方法；迭代法；线性正则化方法 中图分类号：TL99 文献标识码：A O 引言 中子深度分析(NDP)技术是中子活法分析的一种，属于物质近表面分析技术，其原理是利用中子与靶 时会与其发生电离等相互作用，在此过程中损失能量，能损与深度的关系… d昱／s(E)， (1) x=』E…o 其中，x为出射粒子穿越样品物质的长度(深度)；E。为出射粒子的初始能量；E(x)为出射粒子离开样品表面 时的能量．S(E)为样品物质对出射粒子的阻止本领． 测量穿出靶物质的粒子能谱，通过反演得到粒子源沿靶物质表面深度的分布，适用的目标元素如Li， 塞、薄膜和浸洗效应研究等¨-． 对能谱的反演是NDP技术中关键的一步，通过反演可以得到目标元素浓度随样品表面(几微米)深度变 化的分布，因此反演方法的选择在NDP技术中意义重大，不同算法的反演结果之间也会有差别．本文主要对 基于概率迭代推导得到的迭代算法和线性正则化方法的反演效果进行比较研究． 1 物理模型 采用多道分析器对NDP粒子能谱进行测量时，设测量到的粒子能量分布在J7、r个道内．同时将样品沿深 度方向上进行离散，分为M层，将各层通过核反应产生的粒子考虑为一个强度为s，的单能粒子源，其强度与 待测元素的含量、中子注量率和反应截面成正比，在测量系统和参数不变时，该强度的分布即是待测元素含 量的分布． 假设第歹层的一个粒子在测量系统的第i道产生计数的概率为P(C；l jsf)，记为PP。与粒子源性质，基 体材料性质，测量几何及探测器系统性能等相关．P。构成Ⅳ×M的响应矩阵P．设第i道的总计数为Ci，则有 _|l， Ci=：PiSi． (2) 收稿日期：2012—01—24；修回日期：2012—06—08 基金项目：中国工程物理研究院科学技术发展基金(201080103009)资助项目 作者简介：杨鑫(1988一)，男，贵州，博士生，主要从事反应堆物理及应用研究。E-mail：yangx05@126．corn 万方数据 892 计 算 物 理 第29卷 写成矩阵形式为PS=C，反演的任务即是根据响应矩阵P及 测量到的能谱c求解得到粒子源强度分布s．对于方程(2)， 当NM时为不定方程，NM时为超定方程，N=M时有唯 一解，在实际的NDP问题中，根据测量的结果及需求，可对M 与Ⅳ的值进行调节． X 典型的NDP得到的P(NM，NM)矩阵的形式如 图1．根据P