一种基于金刚石探测器测量快中子能谱的方法技术

本发明专利技术公布了一种基于金刚石探测器测量快中子能谱的方法。由于不同能量快中子与碳原子核发生的各类核反应的种类不同、截面及角分布不同、核反应产物的能量不同，金刚石探测器对于不同能量快中子的响应函数不同，该方法首先通过蒙特卡洛模拟金刚石探测器对不同能量快中子的响应函数，得到金刚石探测器的中子响应矩阵，然后测量金刚石探测器对待测中子能谱的响应函数，将实验测量得到的响应函数通过中子响应矩阵反解，得到待测中子能谱。

【技术实现步骤摘要】
一种基于金刚石探测器测量快中子能谱的方法
本专利技术属于中子测量
，具体涉及一种基于金刚石探测器测量快中子能谱的方法。
技术介绍
中子能谱的测量对于核物理的研究具有重要意义，同时也是实现中子技术应用的基础。中子能谱是中子诱发的各种核反应截面的准确测量的基础，这对于核数据、核结构、核反应机制以及核天体等具有重要意义。在应用领域，准确的中子能谱对于反应堆的控制、核反应动力装置以及核聚变系统的诊断等领域具有独特价值。尤其是针对聚变研究，由于聚变装置中的中子注量率非常大，因此要求所使用的中子探测器具有极强的抗辐照性能，而金刚石探测器恰好满足这一特点。
技术实现思路
本专利技术提供了一种基于金刚石探测器测量快中子能谱的方法，以解决现有技术在测量大中子注量率中子能谱的测量问题。具有不同能量的快中子与碳原子核发生核反应的种类不同、各类核反应的截面及角分布不同、核反应产物的能量不同，因此金刚石探测器对于不同能量快中子的响应函数(即金刚石探测器的沉积能量-计数谱)不同。本专利技术基于快中子与碳原子核发生的各类核反应，将中子能谱信息转化为金刚石探测器对中子的响应函数。通过蒙特卡洛模拟的方法得到金刚石探测器的中子响应矩阵，将金刚石探测器对未知中子能谱的响应函数基于中子响应矩阵进行反解，即可得到中子能谱信息。具体的，本专利技术的快中子能谱测量方法包括如下步骤：1)利用蒙特卡洛方法模拟金刚石探测器对不同能量快中子的响应函数，得到金刚石探测器的中子响应矩阵；2)实验测量金刚石探测器对待测中子能谱的响应函数；3)基于步骤1)得到的中子响应矩阵对步骤2)得到的金刚石探测器的响应函数进行反解，得到待测中子能谱。利用蒙特卡洛方法模拟金刚石探测器对不同能量快中子的响应函数的基本思路为：假设某一能量的一个中子进入金刚石探测器，设定一个中子前进步长，中子沿其方向前进一个步长，根据评价核数据库中的中子与碳核各类核反应的截面信息，蒙特卡洛抽样判断中子是否会与碳原子核发生某类核反应：如果不发生核反应，则中子继续前进一个步长；如果发生某类核反应，根据评价核数据库中此类核反应的角分布信息，蒙特卡洛抽样决定某一个核反应产物的出射角度，并根据核反应运动学计算核反应产物的出射能量、出射方向，这次核反应产生的带电粒子，如α粒子、反冲核等，认为其能量沉积在金刚石探测器中并产生一个计数，记录此沉积能量，带电粒子沉积在金刚石中能量的计算需要考虑壁效应的问题，即带电粒子穿出金刚石探测器而只能将部分能量沉积到金刚石探测器中，如果核反应过程中没有中子产生，则认为这个中子“死亡”，如果发生的核反应有中子产生，则新产生中子的能量已通过核反应运动学计算得到，此中子沿其方向继续前进一个步长，重复上述过程，直至中子“死亡”或者其穿出金刚石探测器。重复上述过程多次，即模拟多个某一能量的中子入射，当模拟的中子数目较大时，可以得到统计性较好的金刚石探测器对于此中子能量的沉积能量-计数谱，即金刚石探测器对此中子能量的响应函数，对于不同的中子能量均采用此种模拟方法，可以得到金刚石探测器对不同中子能量的响应函数，即金刚石探测器的中子响应矩阵。不同的金刚石探测器具有不同的能量分辨率，因此模拟金刚石探测器的中子响应矩阵时需要考虑能量分辨率的影响。上述步骤1)中利用蒙特卡洛方法模拟金刚石探测器对不同能量快中子的响应函数的基本流程如图1所示，假设某一能量的一个中子从某一方向、位置入射金刚石探测器，首先设定一个中子前进步长，中子在金刚石探测器中沿其方向前进一个步长，根据评价核数据库中的中子与碳核各类核反应的截面信息，蒙特卡洛抽样判断中子是否会与碳原子核发生某类核反应：如果不发生核反应，则中子继续前进一个步长；如果发生某类核反应，计算核反应出射带电粒子在金刚石探测器中的沉积能量，如果核反应过程中没有中子产生，则认为这个中子“死亡”，如果有中子产生，则计算新产生中子的出射方向和能量，此中子沿其方向继续前进一个步长；重复此过程，直至中子“死亡”或者其穿出金刚石探测器。模拟多个某一能量的中子入射，得到金刚石探测器对于此中子能量的沉积能量-计数谱，即金刚石探测器对此中子能量的响应函数。对于不同的中子能量均采用此种模拟方法，得到金刚石探测器对不同中子能量的响应函数，即金刚石探测器的中子响应矩阵。在上述模拟过程中，蒙特卡洛抽样判断中子是否会与碳原子核发生某类核反应时，如果发生某类核反应，根据评价核数据库中此类核反应的角分布信息，蒙特卡洛抽样决定某一个核反应产物的出射角度，并根据核反应运动学计算核反应产物的出射能量、出射方向，这次核反应产生的带电粒子，如α粒子、反冲核等，认为其能量沉积在金刚石探测器中并产生一个计数，记录此沉积能量。计算带电粒子在金刚石探测器中的沉积能量时需要考虑壁效应的问题和金刚石探测器的能量分辨能力，当金刚石探测器中沉积了能量E0，探测器输出值E并不完全等于E0，而是一个以E0为平均值的高斯分布，这个高斯分布的展宽由具体金刚石探测器的能量分辨能力决定。如图2所示，某一能量、方向、位置的带电粒子在金刚石探测器中沉积一个能量步长，根据金刚石对此种带电粒子的阻止本领求得沉积此能量步长时带电粒子的移动距离，判断该带电粒子是否穿出金刚石探测器，如果没有，则该带电粒子继续在金刚石探测器中沉积一个能量步长，如果穿出了金刚石探测器，则对此带电粒子的追踪结束并得到计算的沉积能量。此时需要考虑金刚石探测器的能量分辨能力，构建一个以上述计算得到的沉积能量为平均值、以实际探测器能量分辨能力为展宽的高斯分布，在此高斯分布中抽取得到此带电粒子在金刚石探测器中的实际总沉积能量。上述步骤2)实验测量金刚石探测器对待测中子能谱的响应函数时，使用电荷灵敏型前置放大器，保证金刚石探测器具有较好的信号质量，噪声较小，测量较长时间以获得足够的统计性。上述步骤3)基于中子响应矩阵对金刚石探测器对待测中子能谱的响应函数进行反解，反解的方法为迭代法。中子能谱φ与金刚石探测器对待测中子能谱的响应函数N之间通过中子响应矩阵R相关联：其中，i为响应函数道址序号，而j为中子能谱道址序号。利用GRAVEL迭代算法反解中子能谱，迭代公式为：其中，k为迭代次数，为第k次迭代获得的中子能谱，Ni为金刚石探测器测量的响应函数，ρi为Ni的标准偏差(如果Ni>0，否则，ρi＝1)，Ri,j为金刚石探测器的中子响应矩阵，为权重因子，其定义如式3：迭代思路如图3所示，从某一初始中子能谱出发进行迭代反解，如果迭代次数达到设定次数(如1000次)，或者当利用中子能谱推出的响应函数与实验测量的响应函数Ni之间的相对偏差降低至预设值(如0.05％)时，迭代停止。步骤3)迭代结束后即可得到中子能谱结果。根据本专利技术的方法，利用蒙特卡洛模拟得到金刚石探测器的快中子响应矩阵，将金刚石探测器对待测中子能谱的响应函数基于中子响应矩阵进行反解，即可得到待测中子能谱。其有益效果简述如下：1、本专利技术利用金刚石探测器测量中子能谱，金本文档来自技高网...

【技术保护点】
1.一种基于金刚石探测器测量快中子能谱的方法，包括如下步骤：/n1)利用蒙特卡洛方法模拟金刚石探测器对不同能量快中子的响应函数，得到金刚石探测器的中子响应矩阵；/n2)实验测量金刚石探测器对待测中子能谱的响应函数；/n3)基于步骤1)得到的中子响应矩阵对步骤2)得到的金刚石探测器对待测中子能谱的响应函数进行反解，得到待测中子能谱。/n

【技术特征摘要】
1.一种基于金刚石探测器测量快中子能谱的方法，包括如下步骤：
1)利用蒙特卡洛方法模拟金刚石探测器对不同能量快中子的响应函数，得到金刚石探测器的中子响应矩阵；
2)实验测量金刚石探测器对待测中子能谱的响应函数；
3)基于步骤1)得到的中子响应矩阵对步骤2)得到的金刚石探测器对待测中子能谱的响应函数进行反解，得到待测中子能谱。


2.如权利要求1所述的方法，其特征在于，步骤1)利用蒙特卡洛方法模拟金刚石探测器对不同能量快中子的响应函数的过程如下：假设某一能量的一个中子从某一方向、位置入射金刚石探测器，首先设定一个中子前进步长，中子在金刚石探测器中沿其方向前进一个步长，根据评价核数据库中的中子与碳核各类核反应的截面信息，蒙特卡洛抽样判断中子是否会与碳原子核发生某类核反应：如果不发生核反应，则中子继续前进一个步长；如果发生某类核反应，计算核反应出射带电粒子在金刚石探测器中的沉积能量，如果核反应过程中没有中子产生，则认为这个中子“死亡”，如果有中子产生，则计算新产生中子的出射方向和能量，此中子沿其方向继续前进一个步长；重复此过程，直至中子“死亡”或者其穿出金刚石探测器；模拟多个某一能量的中子入射，得到金刚石探测器对于此中子能量的沉积能量-计数谱，即金刚石探测器对此中子能量的响应函数；对于不同的中子能量均采用此种模拟方法，得到金刚石探测器对不同中子能量的响应函数，即金刚石探测器的中子响应矩阵。


3.如权利要求2所述的方法，其特征在于，在步骤1)的模拟过程中，蒙特卡洛抽样判断中子是否会与碳原子核发生某类核反应时，如果发生某类核反应，根据评价核数据库中此类核反应的角分布信息，蒙特卡洛抽样决定某一个核反应产物的出射角度，并根据核反应运动学计算核反应产物的出射能量、出射方向，对于这次核反应产生的带电粒子，认为其能量沉积在金刚石探测器中并产生一个计数，记录此沉积能量。


4.如权利要求2所述的方...

【专利技术属性】
技术研发人员：刘杰，张国辉，江浩雨，崔增琪，胡益伟，
申请(专利权)人：北京大学，
类型：发明
国别省市：北京;11