中子半影成像的软件系统.docVIP

中子半影成像的软件系统 赵宗清1,2丁永坤1郝轶聃1 袁永腾1 李志超1 (1.中国工程物理研究院激光聚变研究中心，四川绵阳，621900；2.清华大学工程物理系，北京 100084) ）μm的分辨率。而在即将建成的美国NIF装置和法国LMJ装置上要实现5～10μm的分辨率。 在神光II激光装置上，中子产额较低，不能满足中子成像的要求 [10]。但随着我国ICF事业的发展，在神光III原型装置上DT中子产额有望达到1010/发以上，这为中子成像技术的发展提供了有利的条件。因此，为了提高我国ICF核诊断的水平，我们将大力发展中子半影成像技术。 中子半影成像系统的建设，包括硬件平台建设和软件平台建设两部分。硬件平台的建设包括编码孔的设计加工，图像记录系统研制和瞄准技术研究等。目前我们在编码孔加工技术和瞄准技术方面还有一定的技术瓶颈，需要进行突破。同时，为了考察中子成像技术的可行性，并且为硬件设计提供定量化的参考意见，软件平台的建设必须先于硬件平台建设。 图1 中子半影成像硬件系统 2成像过程的数值模拟 2.1 蒙特卡罗模拟[11] Monte-Carlo方法是一种随机抽样方法。利用已知的中子反应截面数据, 模拟各种微观物理过程, 通过概率抽样对源粒子的行为进行跟踪,决定每次碰撞后次级粒子的运动方向和速率。MCNP作为计算程序。MCNP是一个大型多功能的蒙特卡罗计算程序, 可处理复杂场所三维几何结构的中子光子耦合输运问题。MCNP具有较强的通用性, 在源描述、空间物质的几何分布上具有很大的灵活性,可处理任意三维几何结构问题, 适用面宽。对中子的输运问题,MCNP详细处理了各种微观物理过程,适用于中子成像的要求。利用MCNP模拟中子成像，能使我们更好地理解中子在成像系统中的输运过程，同时能够得到图像重建所需要的点扩散函数，为图像重建算法的研究提供条件。 2.2 解析模拟[12] 通过物理分析和几何建模，我们得到成像传递函数，由两部分的乘积组成，即几何衰减因子 (1) 及吸收衰减因子 (2) 其中为中子束在锥孔材料中的吸收路程，为材料的中子吸收系数。详细内容可参看文献12。 1.3 两种模拟方法的比较 图2 蒙特卡罗方法和解析方法得到的点扩散函数 图2是蒙特卡罗模拟和解析模拟分别得到的点扩散函数。从图2中可以看出，蒙特卡罗模拟得到的图像考虑了散射等物理因素的影响，存在着统计误差，而解析法得到的图像非常光滑。 Monte-Carlo方法的优势在于通过对单中子事例的微观物理过程的计算来实现宏观模拟，因此其结果在包括噪声分布等的各方面均非常接近真实物理过程。而解析方法的优势在于，忽略噪声、对中子束的几何光学近似及采用各种宏观物理现象的唯象理论，使得计算效率极高。 因此，我们的发展方向定位为，采用解析方法对中子成像系统进行定性或半定量估计，并在此基础上采用Monte-Carlo方法做进一步的细致研究，并且采用并行机提高Monte-Carlo的计算效率。 3 图像重建 3.1 成像模型 中子半影成像，主要有两个步骤[5]： 第一步：由中子源发射中子，经过半影孔形成编码像，编码像由亮区、半影区和本底组成。源的亮度分布o(x,y)通过半影孔点扩散函数h(x,y)之后，像的亮度分布i(x,y)为： （3） 其中，**代表二维卷积，n (x,y)为干扰噪声。 第二步：由编码像经过图像反演得到源的亮度分布，即中子通量分布。图像反演的基本思想是：先分别对i(x,y)，h(x,y)和n(x,y)做傅立叶变换得到其频域分布I(kx,ky)，H(kx,ky)和N(kx,ky)，相应地： （4） 通过一定的图像反演算法，就可以得到o(kx,ky)，再进行逆傅立叶变换便可以得到源的中子强度的二维空间分布。 3.2 维纳滤波 维纳滤波就是寻找传输函数为Hw(kx,ky)的滤波器，使得恢复出的目标估计值Ow(x,y)与真实的O(x,y)的均方差达到最小，即： （5） 经过推导，可得 （6） 其中，为编码孔成像系统的传递函数的复共轭，和分别是靶图像和噪声的功率谱密度。 3.3 Richardson-Lucy方法[14,15] 如果成像方程服从泊松分布： （7） O是真实的源分布，是PSF（从物体真实位置j发出的光散射到观察像素i上的份额）。I是无噪声污染的退化图像。观测计数和期望计数的联合概率L为： （8） 经过一系列的推导，可以得到R-L迭代公式为： （9） 从方程（9）可以明显看出，如果Richardson-Lucy迭代收敛[16]，意味着随着迭代进行，修正因子将趋于