[发明专利]大空间深穿透辐射场空气吸收剂量率的计算方法及介质

权利要求书

1.一种大空间深穿透辐射场空气吸收剂量率的计算方法，其特征在于，该方法包括：

根据核材料贮存场景三维模型数据和放射源数据构建辐射场计算模型；

将所述辐射场计算模型的计算模式设置为通量转剂量模式；

将计算方法设定为探测器法与网格计算法相结合的组合方式；以及

根据所述计算模式和所述计算方法计算辐射场的空气吸收剂量率。

2.如权利要求1所述的大空间深穿透辐射场空气吸收剂量率的计算方法，其特征在于，根据所述计算模式和所述计算方法计算辐射场的空气吸收剂量率的步骤之后，该方法还包括：

读取并显示所述辐射场的空气吸收剂量率。

3.如权利要求1所述的大空间深穿透辐射场空气吸收剂量率的计算方法，其特征在于，根据所述计算模式和所述计算方法计算辐射场的空气吸收剂量率的步骤，具体包括：

将探测器设置在辐射场的关键点位；

根据实际应用场景及关注区域对所述辐射场进行网格划分；

利用截断减方差对所述辐射场进行空间截断；

利用指数变换、强迫碰撞联合法引导粒子输运方向；

利用权窗引导粒子向计数单元加速运行；以及

计算所述辐射场的空气吸收剂量率。

4.如权利要求3所述的大空间深穿透辐射场空气吸收剂量率的计算方法，其特征在于，所述关键点位是根据工作人员活动频率、屏蔽物体设置和放射源距离计算得到的。

5.如权利要求3所述的大空间深穿透辐射场空气吸收剂量率的计算方法，其特征在于，所述网格划分的粗细程度是根据预掌握空间分布规律确定的。

6.如权利要求3所述的大空间深穿透辐射场空气吸收剂量率的计算方法，其特征在于，利用指数变换、强迫碰撞联合法引导粒子输运方向的步骤，具体包括：

利用指数变换法调整输运粒子在介质中的微观截面；以及

利用强迫碰撞法提高计数区粒子的碰撞次数。

7.如权利要求6所述的大空间深穿透辐射场空气吸收剂量率的计算方法，其特征在于，所述微观截面的调整公式为：

其中，为调整后的微观截面；为原始微观截面；P为拉伸方向；μ为粒子方向和拉伸方向之间角度的余弦。

8.如权利要求3所述的大空间深穿透辐射场空气吸收剂量率的计算方法，其特征在于，利用权窗引导粒子向计数单元加速运行的步骤，具体包括：

通过权窗产生器和能量权窗产生器联合获得不同网格的能量区间权窗下限分布文件；以及

调用所述能量区间权窗下限分布文件引导粒子向计数单元加速运行。

9.如权利要求1所述的大空间深穿透辐射场空气吸收剂量率的计算方法，其特征在于，所述网格计算法是将整个辐射场划分为多个立方体后计算每个立方体中的粒子数。

10.一种包含计算机执行指令的计算机存储介质，所述计算机执行指令经由数据处理设备处理时，该数据处理设备执行权利要求1～9任一所述的大空间深穿透辐射场空气吸收剂量率的计算方法。