一种基于变权重离散的γ剂量率计算方法及系统技术方案

本发明专利技术公开了一种基于变权重离散的γ剂量率计算方法及系统，方法包括：S1、以体源或面源的几何中心为原点，分别建立X、Y、Z坐标轴；S2、根据待计算点到各坐标轴的距离和体源或面源在各坐标轴方向上的最大尺寸设置离散权重参数；S3、将体源或面源离散为若干点源；S4、计算每个点源的γ剂量率贡献；S5、将每个点源的γ剂量率贡献进行累积，得到体源或面源在待计算点处的γ剂量率；S6、将每个体源和每个面源在待计算点处的γ剂量率进行累积，得到待计算点的最终γ剂量率。本发明专利技术所提供的方法及系统，根据点到源之间距离不同，采用不同的权重，使得体源和面源相应的离散数目不同，不仅保证了γ外照射剂量率计算的精度，还能够大大提高计算的效率。

A Method and System for Calculating Gamma Dose Rate Based on Variable Weight Discrete

The invention discloses a method and system for calculating gamma dose rate based on variable weight discretization, which includes: S1, establishing X, Y and Z coordinate axes respectively from the geometric center of volume source or surface source; S2, setting discrete weight parameters according to the distance between the point to be calculated and the maximum size of volume source or surface source in each coordinate axis direction; S3, discretizing volume source or surface source as if. Dry point source; S4. Calculate the gamma dose rate contribution of each point source; S5. Cumulate the gamma dose rate contribution of each point source to get the gamma dose rate of the body source or the non-point source at the point to be calculated; S6. Cumulate the gamma dose rate of each body source and each non-point source at the point to be calculated to get the final gamma dose rate of the point to be calculated. The method and system provided by the invention adopts different weights according to the different distances between points and sources, so that the corresponding discrete number of body source and non-point source is different, which not only ensures the accuracy of calculation of dose rate of gamma external irradiation, but also greatly improves the efficiency of calculation.

【技术实现步骤摘要】
一种基于变权重离散的γ剂量率计算方法及系统
本专利技术涉及辐射防护
，具体涉及一种基于变权重离散的γ剂量率计算方法及系统。
技术介绍
点核积分方法是计算γ外照射剂量的常用方法，被广泛应用在屏蔽计算和外照射剂量计算中，与蒙特卡罗方法相比，它不仅不受空间尺寸和屏蔽体厚度的限制，而且机时量小，具有较快的计算速度，缺点在于该方法不能对粒子散射问题进行准确的考虑，需要引入累积因子对粒子散射进行修正，计算结果往往偏于保守。该方法以各向同性点源的计算为基础，对于体源或者面源采用离散的方法，转化为若干个各向同性点源，将每个点源的计算结果进行累加，从而得到整个体源或者面源的计算结果。将体源或者面源离散成点源的数目，不仅与计算结果的精度有关，还与所需的计算时间相关。一方面目前传统的点核积分算法没有对计算效率进行优化，计算速度已然不能满足当前的需求，另一方面随着点核积分方法在γ辐射场计算中的优势显现，同时随着辐射场可视化显示的需要，需要计算的代表三维辐射场的点阵数量越来越庞大，这对计算的速度提出了严峻的挑战。
技术实现思路
针对现有技术中存在的缺陷，本专利技术的目的在于提供一种基于变权重离散的γ剂量率计本文档来自技高网...

【技术保护点】
1.一种基于变权重离散的γ剂量率计算方法，其特征在于，包括：S1、在一体源或一面源内部建立局部坐标系，以所述体源或所述面源的几何中心为原点，分别建立X、Y、Z坐标轴；S2、计算一待计算点到各坐标轴的距离，根据所述待计算点到各坐标轴的距离和所述体源或所述面源在各坐标轴方向上的最大尺寸设置离散权重参数；S3、根据所述待计算点到各坐标轴的距离、所述最大尺寸以及所述离散权重参数，将所述体源或所述面源离散为若干点源；S4、根据每个点源到所述待计算点的距离计算每个点源的γ剂量率贡献；S5、将每个点源的γ剂量率贡献进行累积，得到所述体源或所述面源在所述待计算点处的γ剂量率；S6、将每个体源和每个面源在所述待...

【技术特征摘要】
1.一种基于变权重离散的γ剂量率计算方法，其特征在于，包括：S1、在一体源或一面源内部建立局部坐标系，以所述体源或所述面源的几何中心为原点，分别建立X、Y、Z坐标轴；S2、计算一待计算点到各坐标轴的距离，根据所述待计算点到各坐标轴的距离和所述体源或所述面源在各坐标轴方向上的最大尺寸设置离散权重参数；S3、根据所述待计算点到各坐标轴的距离、所述最大尺寸以及所述离散权重参数，将所述体源或所述面源离散为若干点源；S4、根据每个点源到所述待计算点的距离计算每个点源的γ剂量率贡献；S5、将每个点源的γ剂量率贡献进行累积，得到所述体源或所述面源在所述待计算点处的γ剂量率；S6、将每个体源和每个面源在所述待计算点处的γ剂量率进行累积，得到所述待计算点的最终γ剂量率。2.根据权利要求1所述的一种基于变权重离散的γ剂量率计算方法，其特征在于，步骤S3包括：S31、根据所述待计算点到各坐标轴的距离、所述最大尺寸和所述离散权重参数，计算得到所述待计算点在各坐标轴方向上的离散数；S32、从预设的各坐标轴方向上的最大离散数和计算得到所述待计算点在各坐标轴方向上的离散数中取较小者，作为所述待计算点在各坐标轴方向上的最终离散数；S33、将所述体源或所述面源离散为若干点源，点源个数为所述待计算点在各坐标轴方向上的最终离散数的乘积。3.根据权利要求2所述的一种基于变权重离散的γ剂量率计算方法，其特征在于，步骤S31具体包括：如果Ri≠0，则根据下式计算并向上取整后得到所述待计算点在各坐标轴方向上的离散数Ni；其中，n为所述离散权重参数，Di为所述体源或所述面源在各坐标轴方向上的最大尺寸，Ri为所述待计算点到各坐标轴的距离，Ni为所述待计算点在i坐标轴方向上的离散数，Ni为正整数，i＝{i|x,y,z}；如果Ri＝0，即所述待计算点位于坐标轴上，则根据下式计算并向上取整后得到所述待计算点在该坐标轴方向上的离散数N；其中，n为所述离散权重参数，R为所述待计算点到原点的距离，D为所述体源或所述面源在该坐标轴方向上的最大尺寸，N为正整数。4.根据权利要求3所述的一种基于变权重离散的γ剂量率计算方法，其特征在于，步骤S2中，根据所述待计算点到各坐标轴的距离和所述体源或所述面源在各坐标轴方向上的最大尺寸设置离散权重参数n，具体包括：如果则n的值设置为如果则n的值根据实际情况设置为10-20之间的任意数；如果则n的值根据实际情况设置为1-10之间的任意数；其中，Di为所述体源或所述面源在各坐标轴方向上的最大尺寸，Ri为所述待计算点到各坐标轴的距离，i＝{i|x,y,z}。5.根据权利要求1-4任一项所述的一种基于变权重离散的γ剂量率计算方法，其特征在于，步骤S4具体包括：根据下式计算每个点源的γ剂量率贡献；其中，F(E)为转换函数，E为γ光子的能量，Ai为第i个点源的活度，B为累积因子，e-μt为指数衰减项，μ为材料对γ光子的线衰减系数，t为γ光子在该材料中穿过的距离，r为第i个点源到所述待计算点的距离，Hi为第i个点源的γ剂量率贡献，i为正整数。...

【专利技术属性】
技术研发人员：刘立业，李华，赵原，曹勤剑，肖运实，汪屿，夏三强，赵日，卫晓峰，潘红娟，
申请(专利权)人：中国辐射防护研究院，
类型：发明
国别省市：山西,14