(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210577101.X (22)申请日 2022.05.25 (71)申请人 西安交通大 学 地址 710049 陕西省西安市碑林区咸宁西 路28号 (72)发明人 贺清明　舒瀚林　曹良志　 (74)专利代理 机构 西安智大知识产权代理事务 所 61215 专利代理师 何会侠 (51)Int.Cl. G06F 30/23(2020.01) G06F 30/12(2020.01) G06T 17/20(2006.01) (54)发明名称 基于有限元-蒙特卡罗-点核积分耦合获得 γ剂量率分布的方法 (57)摘要 一种基于有限元 ‑蒙特卡罗 ‑点核积分耦合 获得γ剂量率分布的方法， 通过分析核装置辐射 屏蔽问题中源的源强和空间分布自动生成剖分 面， 基于核装置的三维CAD模型以剖分面为边界 划分非结构网格， 采用有限元方法对剖分面内部 区域进行共轭输运计算， 得到源偏倚和权窗参 数， 提供给蒙特卡罗方法对剖分面内部区域进行 前向输运计算， 得到剖分面的中子面源和光子面 源， 将剖分面的中子面源和光子面源作为源项， 采用点核积分方法对剖分面外的区域进行计算， 得到γ剂量率分布； 本发明方法具有较强的复杂 几何适应性， 效率高于传统的蒙特卡罗方法， 在 保证计算精度的情况下极大提高了计算效率， 为 核装置辐射屏蔽问题的快速、 精确计算提供了可 靠方案。 权利要求书1页 说明书3页 附图2页 CN 114781230 A 2022.07.22 CN 114781230 A 1.一种基于有限元 ‑蒙特卡罗 ‑点核积分耦合获得γ剂量率分布的方法， 其特征在于： 包括如下步骤： 步骤1： 根据核装置辐射屏蔽问题中源的空间分布信息和源强， 计算获得每组源的等效 中心； 步骤2： 对于每组源， 以源的等效中心为中心分别画源的外切球体、 外切长方体及外切 圆柱体， 比较外切球体、 外切圆柱体及外切长方体的体积并选取体积最小的几何体的表面 作为剖分面； 步骤3： 根据核装置辐射屏蔽 问题中的几何、 材料信息在CAD软件中建立核装置三维CAD 模型； 步骤4： 基于步骤3建立的核装置三维CAD模型， 利用网格剖分工具以步骤2所确定的剖 分面为边界生成非结构网格； 步骤5： 基于步骤4生成的非结构网格， 采用有限元方法对剖分面内部区域进行共轭输 运计算， 得到源偏倚和权窗参数； 步骤6： 利用步骤5的得到的源偏倚和权窗参数， 采用蒙特卡罗方法对剖分面内部区域 进行前向输运计算， 得到剖分面的中子面源和光子面源； 步骤7： 将步骤6算出的剖分面的中子面源和光子面源作为源项， 采用点核积分方法对 剖分面外的区域进行计算， 得到γ剂量 率分布。权　利　要　求　书 1/1 页 2 CN 114781230 A 2基于有限元 ‑蒙特卡罗 ‑点核积分耦合获得 γ剂量率分布的 方法 技术领域 [0001]本发明涉及核装置辐射屏蔽计算领域， 具体涉及一种基于有限元 ‑蒙特卡罗 ‑点核 积分耦合获得γ剂量 率分布的方法。 背景技术 [0002]确定辐射场水平是开展辐射防护工作的重要依据之一， 理论计算的手段不受时空 限制， 是确定辐射场水平的重要方法。 目前辐射场水平计算有三种基本方法： (1)蒙特卡罗 方法： 非确定论方法， 可求解复杂几何， 但计算时间长且处理深穿透问题时表现不够好； (2) 确定论方法： 如 有限元方法， 计算速度相对较快， 但在求解复杂几何时不如蒙特卡罗方法精 确； (3)点核积分方法： 半经验方法， 计算速度快， 但求 解复杂几何时计算结果显著偏高。 [0003]三种方法各有优劣， 单独使用难以满足工程实际对核装置辐射场水平进行快速、 精确求解的需求。 以核设施退役辐射场计算为例， 在核设施退役过程中， 随着设备的逐步拆 除， 需对其三维辐 射场的γ剂量率进行快速、 精准的计算与评估作为退役策略和辐 射防护 工作的重要依据。 若单独以蒙特卡罗方法进行辐 射场水平理论计算则耗时过长， 限制退役 工作按时开展， 且在部 分重要区域(如一回路周围厂房)将因闯过厚屏蔽结构(如反射层)而 出现计数率偏低的问题， 无法保证计算结果的可靠性； 若单独以有限元方法进行辐 射场水 平理论计算则难以保证含源项复杂几何构件如堆芯所在区域计算结果的正确性； 若单独以 点核积分方法进行辐射场水平理论计算， 将显著高估γ剂量 率， 增加退役工作的成本 。 [0004]针对目前辐射场计算方法无法兼顾计算精度和计算效率的缺点， 发明一种基于有 限元‑蒙特卡罗 ‑点核积分耦合获得γ剂量 率分布的方法。 发明内容 [0005]蒙特卡罗方法具有计算精度 高、 复杂几何适应能力强的优点， 十分适合对核装置 辐射屏蔽问题进行计算， 但计算耗时长， 难以满足工程应用要求， 如果采用有限元方法对求 解区域进行共轭输运计算提供源偏倚和权窗参数， 蒙特卡罗计算的时间将大大缩短。 点核 积分方法的计算速度快， 适于求解不含源项的真空区域或含源项几乎可忽略的区域。 因此 将计算区域进 行剖分， 以最小相切的方法自动选取剖分面， 利用有限元 ‑蒙特卡罗耦合方法 对含源项、 几何复杂的区域进行计算， 为点核积分提供源项， 利用点核积分对不含源项、 几 何简单的区域计算得到γ剂量 率分布。 [0006]为实现上述目的， 本发明采用以下技 术方案予以实施： [0007]一种基于有限元 ‑蒙特卡罗 ‑点核积分耦合获得γ剂量率分布的方法， 包括如下步 骤： [0008]步骤1： 根据核装置辐射屏蔽问题中源的空间分布信息和源强， 计算获得每 组源的 等效中心； [0009]步骤2： 对于每组源， 以源的等效中心为中心分别画源的外切球体、 外切长方体及说　明　书 1/3 页 3 CN 114781230 A 3