一种核反应堆退役方案的数值仿真方法和系统技术方案

本发明专利技术提出一种核反应堆退役方案的数值仿真方法和系统。该方法包括：建立三维模型；对于每一个part的属性进行编辑；在外部边界内，所有未建立part的区域自动填充part；将整个区域自动划分为密铺的正立方体网格；基于自动划分的网格所继承的渲染颜色对其展示；基于网格进行辐射场计算；计算得到辐射场数据，实现辐射场可视化；根据辐射场、人体模型、行进路径和时间等参数进行辐射剂量计算，得到辐射剂量随时间、空间变化的数据；移除一个part；循环执行以上步骤，直至三维模型被全部移除，形成完整的退役方案。本发明专利技术解决了现有系统中不能反映真实世界中的辐射场变化，以及现有模型对应纠错过程中人力、时间成本过大的问题。时间成本过大的问题。时间成本过大的问题。

【技术实现步骤摘要】
一种核反应堆退役方案的数值仿真方法和系统


[0001]本专利技术属于核反应堆数值仿真领域，涉及一种核反应堆退役方案的数值仿真方法和系统。

技术介绍

[0002]核电机组或核反应堆达到使用年限后，将会进行核退役作业，其中核岛退役拆除作业是最具难度与挑战性的工作。由于核岛内部的核反应堆产生的中子、光子对大量构件进行了长时间的辐照，这些构件都不可避免地产生了活化核素，这导致核岛内部必将产生特定的辐射剂量场。
[0003]现有的核退役虚拟仿真系统，可以通过外部导入的三维辐射分布与三维场景模型进行可视化的比对，然后参照核退役作业流程，通过动画编辑虚拟推演核退役作业过程，在过程中进行碰撞检测、人员实时剂量和累积剂量评估等工作。
[0004]上述核退役虚拟仿真系统中，三维辐射分布(以下简称辐射场)的变化来源于用户预先准备好的数据，而在真实世界中，辐射场是由辐射源产生的中子、光子等微观粒子穿越空间而形成的，该粒子穿越空间的过程遵循第一性原理的粒子输运方程。辐射源来自于核岛内部构件的活化核素，因此当用户在虚拟推演核退役作业过程中，一定会模拟添加、移除、移动、拆分构件(即三维场景模型)，则将导致辐射源的位置、强度、种类发生变化，进一步则将导致辐射场的变化。因此，用户预置的辐射场，无法反映虚拟推演核退役作业过程中所引发的辐射场变化，则更进一步无法正确评估人员所受辐射剂量。在另一方面，核退役场场景庞大且琐碎，用户预置的辐射场也需要进行建模才能实现基于粒子输运方程的求解计算，但该计算建模过程与现有系统的建模过程为两个独立过程，相关数据的对应关系很难仅通过人工静态检查的方式予以保障，同时在产生最终方案的过程中会形成大量中间方案，现有系统所采用的MCNP5数据或EXCEL数据仅可通过人工建模计算予以获得，总之这将耗费大量的人力、时间成本。

技术实现思路

[0005]为了克服上述现有技术中存在的问题，本专利技术提出一种核反应堆退役方案的数值仿真方法和系统，通过对CAD技术、粒子输运方程求解技术和源项求解技术进行结合，一方面实现核退役作业过程中对由于三维场景变化而导致的辐射场的动态跟踪和求解，解决现有系统中不能反映真实世界中的辐射场变化的问题，另一方面实现对核退役过程中真实三维场景的建模以及对应辐射场计算建模的统一对应，解决现有系统与方法在模型对应纠错过程中所需的大量人力、时间成本问题。
[0006]本专利技术解决上述问题的技术方案是：
[0007]第一方面，本专利技术提出一种核反应堆退役方案的数值仿真方法，其特殊之处在于，包括以下步骤：
[0008]1)对核退役场景建立可视化的三维模型，三维模型是由多个互不重叠的独立几何
体(以下简称part)构成，同时建立整个三维模型的外部边界；
[0009]2)对于每一个part的属性进行编辑，属性具有三层的结构化数据结构，将对每个part的属性予以赋值；
[0010]3)在外部边界内，所有未建立part的区域自动填充part；
[0011]4)基于三维模型的外部边界，将整个区域自动划分为密铺的正立方体网格，网格边长可以编辑，每个网格单元含有与part相同的属性数据结构；
[0012]5)基于自动划分的网格继承part的渲染颜色属性，通过该属性进行基于网格的可视化展示；同时，用户可以通过编辑模块合并多个相邻的网格或拆分单个网格；
[0013]6)基于网格进行辐射场计算，利用粒子输运方程求解技术和源项求解技术，其中粒子输运方程求解技术包含点核积分方法或离散纵坐标方法或蒙特卡罗方法；
[0014]7)计算执行完成后，得到辐射场数据，将网格辐射通量值在三维空间渲染展示(以下简称辐射场云图)，辐射场云图和三维模型可共同展示，实现辐射场可视化；
[0015]8)根据人体基本参数建立简化人体模型，对行进路径和时间相关参数进行编辑，利用根据辐射场、人体模型、行进路径和时间等参数进行辐射剂量计算，得到辐射剂量随时间、空间变化的数据；
[0016]9)添加、移除、移动、拆分三维模型中的一个part；
[0017]10)循环执行步骤1)至步骤9)，直至整个三维模型被全部移除，即代表形成了一个完整的退役方案。
[0018]进一步地，上述步骤4)中，自动划分出的网格分为两类，第一类为仅与一个part重叠的网格，该类网格的属性继承与其重叠的part。第二类为与一个以上的part重叠的网格(即网格中出现了多个part的交界面)，该网格的属性继承在该网格中所占体积最大的part的属性。除此之外，网格还有一个额外的辐射通量值属性，该属性的值由后续的辐射场计算模块计算得出。
[0019]进一步地，上述步骤6)中，点核积分的粒子输运方法适合快速计算复杂几何空间辐射场问题，对辐射源项进行空间和能量离散使其成为点源，空间中某一点位置受点源辐射的注量率为：
[0020][0021]其中，S(r',E)为位置r'能量E处的源强，x为光学距离，B(E,x)为光学距离x处能量为E下的累积因子，光学距离的计算方法为：
[0022][0023]其中，N为射线穿过材料的总数目，μ
i
为第i种材料的线性衰减系数，t
i
为射线在第i种材料中穿越的距离；本专利技术所用点核积分求解核心为NECP
‑
MCX
‑
PK。
[0024]进一步地，上述步骤6)中，离散纵标法的粒子输运方法适合较快速较精确的计算辐射场问题，该方法通过把粒子通量的角度分布按照一套求积组对应的离散方向进行离散化，构造离散的代数方程组进行求解，离散方法如下：
[0025][0026]其中，φ
g
(r)为位置r能群g处的标通量，φ
g
(r,Ω)为位置r能群g处方向为Ω的角通量，w
m
为第m个离散方向Ω
m
对应的求积组权重，本专利技术所用离散纵标法求解核心为NECP
‑
Hydra。
[0027]进一步地，上述步骤6)中，蒙特卡罗方法适合复杂几何的精确辐射场计算，该方法模拟粒子在介质中的输运过程，通过统计计数的方法获得辐射场分布，抽样计算方法如下：
[0028][0029]其中S(r,Ω,E)为相空间(r,Ω,E)处的源强，通过源项定义可以进行抽样，K(r',Ω',E'
→
r,Ω,E)为迁移算符，通过抽样粒子的输运和碰撞可以获得，Q(r,Ω,E)为相空间处的发射密度，由发射密度即可获得通量值：
[0030][0031]其中，τ(r,r
‑
lΩ,E)为能量E处，r
‑
lΩ到r处的光学厚度，积分过程也通过抽样统计进行计算，本专利技术中的蒙特卡罗求解核心为NECP
‑
MCX。
[0032]进一步地，上述步骤7)中，对于辐射场云图透明度可以进行设置，可任意设置平面对辐射场云图进行切片，以查看辐射场的某一平面的辐射通量分布，辐射场分布的数据本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种核反应堆退役方案的数值仿真方法，其特征在于，包括以下步骤：1)对核退役场景建立可视化的三维模型，三维模型是由多个互不重叠的独立几何体part构成，同时建立整个三维模型的外部边界；2)对于每一个part的属性进行编辑，属性具有三层的结构化数据结构，将对每个part的属性予以赋值；3)在外部边界内，所有未建立part的区域自动填充part；4)基于三维模型的外部边界，将整个区域自动划分为密铺的正立方体网格，网格边长可以编辑，每个网格单元含有与part相同的属性数据结构；5)基于自动划分的网格继承part的渲染颜色属性，通过该属性进行基于网格的可视化展示；同时，用户可以通过编辑模块合并多个相邻的网格或拆分单个网格；6)基于网格进行辐射场计算，该计算利用粒子输运方程求解技术和源项求解技术，其中粒子输运方程求解技术包含点核积分方法或离散纵坐标方法或蒙特卡罗方法；7)计算执行完成后，得到辐射场数据，将网格辐射通量值在三维空间进行渲染展示，得到辐射场云图，辐射场云图和三维模型可共同展示，实现辐射场可视化；8)根据人体基本参数建立简化人体模型，对行进路径和时间相关参数进行编辑，利用根据辐射场、人体模型、行进路径和时间参数进行辐射剂量计算，得到辐射剂量随时间、空间变化的数据；9)添加、移除、移动、拆分三维模型中的一个part；10)循环执行步骤1)至步骤9)，直至整个三维模型被全部移除，即代表形成了一个完整的退役方案。2.根据权利要求1所述的一种核反应堆退役方案的数值仿真方法，其特征在于：所述步骤4)中，自动划分出的网格分为两类，第一类为仅与一个part重叠的网格，该类网格的属性继承与其重叠的part；第二类为与一个以上的part重叠的网格，即网格中出现了多个part的交界面，该网格的属性继承在该网格中所占体积最大的part的属性；除此之外，网格还有一个额外的辐射通量值属性。3.根据权利要求1或2所述的一种核反应堆退役方案的数值仿真方法，其特征在于：所述步骤6)中，点核积分的粒子输运方法适合快速计算复杂几何空间辐射场问题，对辐射源项进行空间和能量离散使其成为点源，空间中某一点位置受点源辐射的注量率为：其中，S(r',E)为位置r'能量E处的源强，x为光学距离，B(E,x)为光学距离x处能量为E下的累积因子，光学距离的计算方法为：其中，N为射线穿过材料的总数目，μ
i
为第i种材料的线性衰减系数，t
i
为射线在第i种材料中穿越的距离。4.根据权利要求1或2所述的一种核反应堆退役方案的数值仿真方法，其特征在于：所述步骤6)中，离散纵标法的粒子输运方法适合较快速较精确的计算辐射场问题，该
方法通过把粒子通量的角度分布按照一套求积组对应的离散方向进行离散化，构造离散的代数方程组进行求解，离散方法如下：其中，φ
g
(r)为位置r能群g处的标通量，φ
g
(r,Ω)为位置r能群g处方向为Ω的角通量，w
m
为第m个离散方向Ω
m
对应的求积组权重。5.根据权利要求1或2所述的一种核反应堆退役方案的数值仿真方法，其特征在于：步骤6)中，蒙特...

【专利技术属性】
技术研发人员：徐嘉隆，郑琪，梁亮，黄星，谭景烋，巨佳龙，
申请(专利权)人：西安核创能源科技有限公司，
类型：发明
国别省市：