等离子体图像诊断中带吸收校正及轴对称重建.docVIP

等离子体图像诊断中带吸收校正的轴对称重建 朱宏权1,2，王奎禄2，宋顾周2 1. 清华大学工程物理系，北京，100084 2. 西北核技术研究所，西安，710024 摘要：在轴对称等离子体的图像诊断中，根据等离子体在某个方向的X射线强度图像可以重建出辐射系数的空间分布。由于X射线在等离子体传播过程中存在一定吸收等离子体FBP迭代算法。通过对钠等离子体辐射系数的轴对称重建模拟，发现该算法重建精度高，收敛性好，取得了很好的重建结果。 关键词：等离子体，轴对称重建，滤波反投影，吸收校正 中图分类号：O571.44 TP391.9 文献标识码:A 1 引言 对等离子体辐射的X射线空间分的研究,可以给出有关靶加热过程发生的重要数据。X射线针孔相机(PHC)等离子体沿观测方的积分等离子体具有的性质那么逆Abel变换等离子体内的X射线径向分布。但实际上X射线在等离子体传播过程中存在一定吸收当吸收较大时，光性薄假设在此情况下，的影响。江少恩含有吸收的逆Abel变换。，Morozumi迭代法的基础上，结合X射线在等离子体FBP迭代算法。然后，以温度为100eV的钠等离子体为例，对单峰分布和环状分布两种情况下等离子体图像的轴对称重建模拟表明，该算法重建精度高，逆Abel变换。 2.1 等离子体的吸收系数 对于高温稠密等离子体，在局域热动平衡的条件下，根据细致平衡原理，系统的自发辐射系数与吸收系数之间的关系为 (1) 其中，是光子频率，是普朗克常数，是波尔兹曼常数，为黑体辐射强度，其表达式为 （2） 诱发辐射系数为 (3) 其中为等离子体的光子辐射强度。 通常将诱发辐射系数看作负吸收项，将其与吸收系数合并，得到有效吸收系数为 (5) 根据式(1)，上式可以转换为 (6) 取 (7) 对于温度为100eV的钠等离子体，能量在150eV～450eV范围内的辐射光子[3]，的变化曲线如图1所示。 图1 随光子能量的变化曲线 由图1知，在上述能量范围内可以近似为一个常数，式(6)可以简化为 (8) 其中，，即有效吸收系数正比于自发辐射系数。 2.2 等离子体的吸收衰减 等离子体的辐射输运方程为： (9) 其中，是光子沿传播方向的径迹长度。上式左边第一项表示等离子内辐射强度随时间的增量，左边第二项表示等离子内辐射强度随光子传播方向的增量。由于光子通过等离子体的时间很短，等离子体辐射强度在该时间段内的变化可以忽略。再根据式(8)，上式可以改写为： (10) 解微分方程，得到： (11) 可以看出，在给定的情况下，钠等离子体某个方向的辐射强度只与沿该方向的自发辐射系数有关。而自发辐射系数正比于钠等离子体密度的平方，所以，重建出辐射系数的分布，在一定程度上也可以反映出钠等离子体密度的空间分布情况。 3 带吸收校正的迭代算法 FBP法是目前CT重建问题中一种实时高效的方法，它根据逆Radon变换公式，由0～(方向的投影数据来重建物体的三维空间分布。在满足轴对称的假设条件下，我们可以认为等离子体在各个方向的X射线投影图像都相同，因此，可以采用FBP的方法来重建等离子体区域的自发辐射系数分布。 若不考虑等离子体吸收衰减的影响，根据Radon变换，可以由自发辐射系数j的空间分布，得到任意方向的辐射强度I (如图2)。其公式为： （11） 其中R为等离子体的半径，为s轴与x轴的夹角，为方向的辐射强度数据。 图2 FBP法示意图 反之，由0～(方向的辐射强度数据做反向投影，并在空间中相同的位置处迭加，并与重建滤波器卷积，就可以重建自发辐射系数的空间分布，这就是逆Radon变换(即FBP法)。可以表示为： （12） 其中h为重建滤波器，最简单的为R-L滤波器，其定义为： （13） 其中d为s的取样[4] 如果考虑钠等离子体的吸收衰减，根据式（11），得到带衰减的Radon变换为： （14） 其中，为光学厚度。 由于带衰减的逆Radon变换很难求解，所以，根据带衰减Radon变换，FBP法与Morozumi迭代法[5]的原理，得到FBP迭代算法如下： 用FBPI做逆Radon变换得到自发辐射系数分布。 对辐射系数分布做带衰减的Radon变换得到强度曲线。 采用动态校正因子(DMF)作衰减校正。DMF为无衰减的与有衰减的之比。把乘以 DMF得到新的。 重复)和3)，直到辐射强度误差收敛时，停止迭代。 4 数值计算结果