[发明专利]一种非均匀源项分布的γ辐射场数据修正计算方法及系统

权利要求书

1.一种非均匀源项分布的γ辐射场数据修正计算方法，其特征在于，包括：

S1、根据核设施运行数据、源项的空间分布信息和源项特征，对所述核设施内部的源项进行多种区域划分，得到每种区域划分下的多个区域源项；

S2、根据所述核设施的部分实测数据包括的每个测量位置的位置信息计算得到每种区域划分下的每个区域源项在每个测量位置处的剂量率贡献，同时结合对应测量位置处的剂量率实测数据，计算得到每个区域源项的活度值，根据获得的活度值计算得到对应测量位置处的剂量率计算数据；

S3、根据每个测量位置处的剂量率实测数据与对应的剂量率计算数据之间的偏差值，计算得到每个区域源项的活度值的不确定度；

S4、根据每种区域划分下的每个区域源项在每个测量位置处的剂量率贡献的变化程度，以及每个区域源项的活度值的不确定度，确定每种区域划分下的总不确定度；

S5、根据每种区域划分下的总不确定度，选出最优化的区域划分，基于所述最优化的区域划分计算得到最终的γ辐射场数据；

步骤S2包括：

S21、根据所述核设施的部分实测数据中包括的每个测量位置的位置信息，对每种区域划分下的每个区域源项在每个测量位置处的剂量率贡献进行计算，得到响应矩阵；

S22、通过矩阵转置及相乘运算，得到所述响应矩阵对应的法方程组的系数矩阵和常数项；

S23、结合每个测量位置处的剂量率实测数据，通过期望最大化算法对所述法方程组进行迭代求解，得到每个区域源项的活度值，根据获得的源活度值计算得到对应测量位置处的剂量率计算数据。

2.根据权利要求1所述的一种非均匀源项分布的γ辐射场数据修正计算方法，其特征在于，步骤S21中，所述部分实测数据中包括的测量位置的个数大于当前区域划分下的区域源项的个数。

3.根据权利要求1所述的一种非均匀源项分布的γ辐射场数据修正计算方法，其特征在于，步骤S23中，在迭代过程中，设置最大迭代次数，当迭代次数超过所述最大迭代次数且仍不收敛时，替换或增加实测数据。

4.一种非均匀源项分布的γ辐射场数据修正计算系统，其特征在于，包括：

区域划分模块，用于根据核设施运行数据、源项的空间分布信息和源项特征，对所述核设施内部的源项进行多种区域划分，得到每种区域划分下的多个区域源项；

第一计算模块，用于根据所述核设施的部分实测数据包括的每个测量位置的位置信息计算得到每种区域划分下的每个区域源项在每个测量位置处的剂量率贡献，同时结合对应测量位置处的剂量率实测数据，计算得到每个区域源项的活度值，根据获得的活度值计算得到对应测量位置处的剂量率计算数据；

第二计算模块，用于根据每个测量位置处的剂量率实测数据与对应的剂量率计算数据之间的偏差值，计算得到每个区域源项的活度值的不确定度；

确定模块，用于根据每种区域划分下的每个区域源项在每个测量位置处的剂量率贡献的变化程度，以及每个区域源项的活度值的不确定度，确定每种区域划分下的总不确定度；

第三计算模块，用于根据每种区域划分下的总不确定度，选出最优化的区域划分，基于所述最优化的区域划分计算得到最终的γ辐射场数据；

所述第一计算模块包括：

第一计算子模块，用于根据所述核设施的部分实测数据中包括的每个测量位置的位置信息，对每种区域划分下的每个区域源项在每个测量位置处的剂量率贡献进行计算，得到响应矩阵；

运算子模块，用于通过矩阵转置及相乘运算，得到所述响应矩阵对应的法方程组的系数矩阵和常数项；

第二计算子模块，用于结合每个测量位置处的剂量率实测数据，通过期望最大化算法对所述法方程组进行迭代求解，得到每个区域源项的活度值，根据获得的源活度值计算得到对应测量位置处的剂量率计算数据。

5.根据权利要求4所述的一种非均匀源项分布的γ辐射场数据修正计算系统，其特征在于，所述第一计算子模块中，所述部分实测数据中包括的测量位置的个数大于当前区域划分下的区域源项的个数。

6.根据权利要求4所述的一种非均匀源项分布的γ辐射场数据修正计算系统，其特征在于，所述第二计算子模块中，在迭代过程中，设置最大迭代次数，当迭代次数超过所述最大迭代次数且仍不收敛时，替换或增加实测数据。