CN110489789B - 一种核设施退役环境中放射性气体扩散评估方法 - Google Patents

一种核设施退役环境中放射性气体扩散评估方法 Download PDF

Info

Publication number
CN110489789B
CN110489789B CN201910617663.0A CN201910617663A CN110489789B CN 110489789 B CN110489789 B CN 110489789B CN 201910617663 A CN201910617663 A CN 201910617663A CN 110489789 B CN110489789 B CN 110489789B
Authority
CN
China
Prior art keywords
model
radioactive gas
smoke
gas
decommissioning
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910617663.0A
Other languages
English (en)
Other versions
CN110489789A (zh
Inventor
刘永阔
杨立群
彭敏俊
晁楠
龙泽宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201910617663.0A priority Critical patent/CN110489789B/zh
Publication of CN110489789A publication Critical patent/CN110489789A/zh
Application granted granted Critical
Publication of CN110489789B publication Critical patent/CN110489789B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种核设施退役环境中放射性气体扩散评估方法,属于核退役仿真领域。包括构建退役场景数学模型;采用均匀球形烟团模型模拟放射性气体的扩散;使用退役场景的数学模型对均匀球形烟团模型进行修正;计算放射性气体的活度分布,实现退役环境中放射性气体的扩散评估。本发明采用均匀球形烟团模型模拟放射性气体的扩散;考虑了实际退役场景对模型的影响,实现了核设施退役环境中放射性气体浓度分布的动态计算。本发明开发出了一种实时、高效的核设施退役环境中放射性气体扩散评估仿真方法,对核设施退役仿真具有重大的实际意义,应用前景广阔。

Description

一种核设施退役环境中放射性气体扩散评估方法
技术领域
本发明涉及一种核设施退役环境中放射性气体扩散评估方法,属于核退役仿真领域。
背景技术
核设施退役活动会产生大量的放射性气体。为了保证工人在退役活动中的安全,减少放射性气体对工人的内照射,必须分析和评估放射性气体的分布情况。
最早开发的放射性气体扩散模型是高斯模型,它主要应用于点源的扩散研究,包括高斯烟羽模型和高斯烟团模型。前者应用于连续源的扩散研究,后者应用于瞬时泄漏的扩散研究。CFD模拟流体流动是解决气体流动的新兴方法。CFD将数值计算方法和数据可视化技术有机结合起来,能够对流动、换热等相关物理现象进行模拟分析。但由于退役场景空间大,设备所占空间多,空气流通缓慢且放射性气体泄露点多,现有的这些气体扩散模型并不适用于退役场景放射性气体的扩散模拟。
综上所述,开发出一种实时、高效的核设施退役环境中放射性气体扩散评估仿真方法对核设施退役仿真具有重大的实际意义。
发明内容
本发明的目的是为了开发出一种实时、高效的核设施退役环境中放射性气体扩散评估仿真方法而提供一种核设施退役环境中放射性气体扩散评估方法。
本发明的目的是这样实现的,一种核设施退役环境中放射性气体扩散评估方法,具体包括以下步骤:
步骤1、根据核设施退役场景的参数构建虚拟退役环境的数学模型;
步骤2、构建放射性气体均匀球形烟团模型;
步骤3、考虑厂房对模型的影响;
步骤4、考虑退役场景对模型的影响;
步骤5、考虑多个泄露点对模型的影响;
步骤6、考虑排风管道对模型的影响。
本发明还包括这样一些结构特征:
1、步骤1所述数学模型包括厂房、设备、墙体,使用数学公式描述。
2、所述步骤2具体为:放射性气体均匀球形模型通过连续释放一系列球形均匀浓度烟团来模拟放射性气体的扩散,某个空间位置的气体浓度,是所有包含该位置的烟团浓度的累加和,在t时刻时,均匀球形烟团模型的形状公式为:
(x-xi,t(j))2+(y-yi,t(j))2+(z-zi,t(j))2≤Ri,t(j)2
式中,xi,t(j)、yi,t(j)、zi,t(j)为第i个泄露点释放的第j个烟团的中心位置,可根据泄露点的坐标与厂房内气体的流通速度计算得到;Ri(j)是球形烟团模型的半径,可根据放射性气体的扩散速度获得,球形烟团模型内放射性气体的浓度为:
Figure GDA0002217389740000021
式中,Ci,t(j)为t时刻第i个泄露点释放的第j个烟团内放射性气体的浓度,Bq/m3;Qi(j)为第i个泄露点释放第j个烟团时泄露点处放射性气体的释放率,Bq/s;△t为泄露点释放相邻气团的时间间隔,s;Vi,t(j)是t时刻第i个泄露点释放的第j个烟团的体积,单位是m3,理想情况下,烟团模型的体积为:
Figure GDA0002217389740000022
3、所述步骤3具体为:考虑放射性气体烟团形状受到厂房空间的限制,将均匀球形烟团模型的形状公式公式修正为:
Figure GDA0002217389740000023
式中,xmin、ymin、zmin、xmax、ymax、zmax为厂房空间的边界,烟团的体积修正后为:
Figure GDA0002217389740000024
其中,Vout为烟团在厂房边界外的体积。
4、所述步骤4具体为:考虑放射性气体烟团模型与场景的数学模型差运算,对烟团形状公式进行修正,修正后烟团的体积为:
Figure GDA0002217389740000025
其中,Vin为烟团与核设施相交部分的体积。
5、所述步骤5具体为:由于核退役场景中放射性气体泄露点多,且每个烟团都是独立的,因此各空间位置处放射性气体的浓度,为所有泄露点产生的包含该位置烟团的浓度的累加和。
6、所述步骤6具体为:当烟团与排风管道相接时,烟团模型形状不变,单位时间内通过排风管道排走的气体活度为烟团的浓度乘以排风管的气体流量,修正后烟团的浓度为:
Figure GDA0002217389740000026
其中,Vf为单位时间内排风管道抽走气体的体积。
与现有技术相比,本发明的有益效果是:本发明开发出一种基于数学建模方法,将放射性气体简化为内部均匀的球形烟团模型,并动态的计算放射性气体空间分布的仿真方法。本发明实现了采用均匀球形烟团模型模拟放射性气体的扩散;本发明考虑了实际退役场景对模型的影响,实现了放射性气体浓度分布的动态计算。
附图说明
图1是(x,y,z)位置处放射性气体的浓度;
图2是放射性气体扩散过程中体积受厂房空间的限制;
图3是放射性气体扩散过程中体积受场景设施的影响;
图4是多泄露点情况下(x,y,z)位置处的放射性气体的浓度;
图5是放射性气体扩散过程中受排风管道的影响。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述。
如附图所示,图1是(x,y,z)位置处放射性气体的浓度;图2是放射性气体扩散过程中体积受厂房空间的限制;图3是放射性气体扩散过程中体积受场景设施的影响;图4是多泄露点情况下(x,y,z)位置处的放射性气体的浓度;图5是放射性气体扩散过程中受排风管道的影响。
本发明提出一种针对核退役环境中泄露的放射性气体,采用均匀球形烟团模型模拟放射性气体的扩散,并动态的计算放射性气体浓度分布的仿真方法。
本发明采用C++语言编程,主要功能为:构建退役场景的数学模型,并采用均匀球形烟团模型模拟放射性气体的扩散,使用退役场景的数学模型对均匀球形烟团模型进行修正,最终计算放射性气体的活度分布。整个发明包括退役环境建模、均匀球形烟团模型建模、放射性气体浓度分布计算三个模块。
本发明采用以下技术方案:
1、根据核设施退役场景的参数构建虚拟退役环境的数学模型。
根据核设施退役场景的确定参数,构建退役环境的数学模型。数学模型包括厂房、设备、墙体等,可以使用数学公式进行描述。
2、构建放射性气体均匀球形烟团模型。
放射性气体均匀球形烟团模型有以下假设条件:
(a)、泄露处每次产生的烟团是有边界的圆球形,且烟团内放射性气体的浓度分布均匀;
(b)、扩散气体与实物接触时完全不吸收;
(c)、忽略浮力及重力作用,气团中不发生化学反应;
(d)、放射性气体扩散速度不随时间变化,即为匀速向四周扩散;
(e)、气团中心的移动速度等于环境风速。
放射性气体均匀球形模型通过连续释放一系列球形均匀浓度烟团来模拟放射性气体的扩散。在扩散过程中,每个烟团都具有独立行为,不断发生空间位移与形状改变。某个空间位置的气体浓度,是所有包含该位置的烟团浓度的累加和。如图1所示,在放射性气体泄漏持续时间内,泄露源释放了3个烟团,则在(x,y,z)位置处放射性气体的浓度等价为包含该位置所有烟团浓度的叠加。
在t时刻时,均匀球形烟团模型的形状公式为:
(x-xi,t(j))2+(y-yi,t(j))2+(z-zi,t(j))2≤Ri,t(j)2
式中,xi,t(j)、yi,t(j)、zi,t(j)为第i个泄露点释放的第j个烟团的中心位置,可根据泄露点的坐标与厂房内的气体的流通速度计算得到;Ri(j)是球形烟团模型的半径,可根据放射性气体的扩散速度获得。
根据假设条件,烟团内放射性气体的浓度是均匀分布的,因此球形烟团模型内放射性气体的浓度为:
Figure GDA0002217389740000041
式中,Ci,t(j)为t时刻第i个泄露点释放的第j个烟团内放射性气体的浓度,Bq/m3;Qi(j)为第i个泄露点释放第j个烟团时泄露点处放射性气体的释放率,Bq/s;△t为泄露点释放相邻气团的时间间隔,s;Vi,t(j)是t时刻第i个泄露点释放的第j个烟团的体积,m3。在理想情况下,均匀球形烟团模型为圆球形,体积为:
Figure GDA0002217389740000042
3、考虑厂房对模型的影响。
放射性气体扩散形状会受到厂房空间的限制,因此需要对烟团形状公式进行修正:
Figure GDA0002217389740000043
式中,xmin、ymin、zmin、xmax、ymax、zmax为厂房空间的边界。如图2所示,烟团的体积修正后为:
Figure GDA0002217389740000044
其中,Vout为烟团在厂房边界外的体积。
4、考虑退役场景对模型的影响。
核退役场景中的设施体积巨大,会影响放射性气体的扩散。因此需要对烟团模型与场景的数学模型进行差运算,对烟团形状公式进行修正。如图3所示,烟团的体积修正后为:
Figure GDA0002217389740000051
其中,Vin为烟团与核设施相交部分的体积。
5、考虑多个泄露点对模型的影响。
核退役场景中,放射性气体泄露点多,且每个泄露点的泄露量也不相同。如图4所示,由于每个烟团都是独立的,因此(x,y,z)位置处的放射性气体的浓度,为所有泄露点产生的包含该位置烟团的浓度的累加和。
6、考虑排风管道对模型的影响。
放射性气体会通过排风管道,经过滤、吸收等操作后排放到外界环境。如图5所示,当烟团与排风管道相接时,假设模型形状不变,则单位时间内通过排风管道排走的气体活度为烟团的浓度乘以排风管的气体流量。修正后烟团的浓度为:
Figure GDA0002217389740000052
其中,Vf为单位时间内排风管道抽走气体的体积。
综上,本发明公开一种核设施退役环境中放射性气体扩散评估方法,具体涉及一种针对核退役环境中泄露的放射性气体,采用均匀球形烟团模型模拟放射性气体的扩散,并动态的计算放射性气体浓度分布的仿真方法。本发明包括:构建退役场景数学模型;采用均匀球形烟团模型模拟放射性气体的扩散;使用退役场景的数学模型对均匀球形烟团模型进行修正;计算放射性气体的活度分布,实现退役环境中放射性气体的扩散评估。本发明包括退役环境建模、均匀球形烟团模型建模、放射性气体浓度分布计算三个模块,实现了核设施退役环境中放射性气体扩散的动态计算。

Claims (3)

1.一种核设施退役环境中放射性气体扩散评估方法,其特征在于,具体包括以下步骤:
步骤1、根据核设施退役场景的参数构建虚拟退役环境的数学模型;
步骤2、构建放射性气体均匀球形烟团模型;
放射性气体均匀球形模型通过连续释放一系列球形均匀浓度烟团来模拟放射性气体的扩散,某个空间位置的气体浓度,是所有包含该位置的烟团浓度的累加和,在t时刻时,均匀球形烟团模型的形状公式为:
(x-xi,t(j))2+(y-yi,t(j))2+(z-zi,t(j))2≤Ri,t(j)2
式中,xi,t(j)、yi,t(j)、zi,t(j)为第i个泄露点释放的第j个烟团的中心位置,可根据泄露点的坐标与厂房内气体的流通速度计算得到;Ri,t(j)是球形烟团模型的半径,可根据放射性气体的扩散速度获得,球形烟团模型内放射性气体的浓度为:
Figure FDA0004037519680000011
式中,Ci,t(j)为t时刻第i个泄露点释放的第j个烟团内放射性气体的浓度,Bq/m3;Qi(j)为第i个泄露点释放第j个烟团时泄露点处放射性气体的释放率,Bq/s;△t为泄露点释放相邻气团的时间间隔,s;Vi,t(j)是t时刻第i个泄露点释放的第j个烟团的体积,单位是m3,理想情况下,烟团模型的体积为:
Figure FDA0004037519680000012
步骤3、考虑厂房对模型的影响;
考虑放射性气体烟团形状受到厂房空间的限制,将均匀球形烟团模型的形状公式修正为:
Figure FDA0004037519680000013
式中,xmin、ymin、zmin、xmax、ymax、zmax为厂房空间的边界,烟团的体积修正后为:
Figure FDA0004037519680000014
其中,Vout为烟团在厂房边界外的体积;
步骤4、考虑退役场景对模型的影响;
考虑放射性气体烟团模型与场景的数学模型差运算,对烟团形状公式进行修正,修正后烟团的体积为:
Figure FDA0004037519680000015
其中,Vin为烟团与核设施相交部分的体积;
步骤5、考虑多个泄露点对模型的影响;
步骤6、考虑排风管道对模型的影响;
当烟团与排风管道相接时,烟团模型形状不变,单位时间内通过排风管道排走的气体活度为烟团的浓度乘以排风管的气体流量,修正后烟团的浓度为:
Figure FDA0004037519680000021
其中,Vf为单位时间内排风管道抽走气体的体积。
2.根据权利要求1所述一种核设施退役环境中放射性气体扩散评估方法,其特征在于,步骤1所述数学模型包括厂房、设备、墙体,使用数学公式描述。
3.根据权利要求2所述一种核设施退役环境中放射性气体扩散评估方法,其特征在于,所述步骤5具体为:由于核退役场景中放射性气体泄露点多,且每个烟团都是独立的,因此各空间位置处放射性气体的浓度,为所有泄露点产生的包含该位置烟团的浓度的累加和。
CN201910617663.0A 2019-07-10 2019-07-10 一种核设施退役环境中放射性气体扩散评估方法 Active CN110489789B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910617663.0A CN110489789B (zh) 2019-07-10 2019-07-10 一种核设施退役环境中放射性气体扩散评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910617663.0A CN110489789B (zh) 2019-07-10 2019-07-10 一种核设施退役环境中放射性气体扩散评估方法

Publications (2)

Publication Number Publication Date
CN110489789A CN110489789A (zh) 2019-11-22
CN110489789B true CN110489789B (zh) 2023-04-18

Family

ID=68546916

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910617663.0A Active CN110489789B (zh) 2019-07-10 2019-07-10 一种核设施退役环境中放射性气体扩散评估方法

Country Status (1)

Country Link
CN (1) CN110489789B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111695762B (zh) * 2020-04-29 2023-05-05 中国核电工程有限公司 核事故扩散结果的修正方法、装置及后果评价方法、***

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107526852A (zh) * 2016-06-21 2017-12-29 中国辐射防护研究院 一种核设施事故场外后果实时在线评价方法及***
CN107526908A (zh) * 2016-06-21 2017-12-29 中国辐射防护研究院 核事故场外后果评价中拉格朗日烟团大气扩散模拟方法
CN108563839A (zh) * 2018-03-23 2018-09-21 哈尔滨工程大学 一种核设施退役模型程式化仿真方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6966236B2 (en) * 2002-04-12 2005-11-22 New Millennium Nuclear Technologies, Llp Selective profiling of bulk building material for characterization using a concrete contamination profiling technology

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107526852A (zh) * 2016-06-21 2017-12-29 中国辐射防护研究院 一种核设施事故场外后果实时在线评价方法及***
CN107526908A (zh) * 2016-06-21 2017-12-29 中国辐射防护研究院 核事故场外后果评价中拉格朗日烟团大气扩散模拟方法
CN108563839A (zh) * 2018-03-23 2018-09-21 哈尔滨工程大学 一种核设施退役模型程式化仿真方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Drivas P J 等.Modeling Indoor Air Exposure from Short‐Term Point Source Releases.《Indoor Air》.2004,第06卷第271-277页. *
Torben Mikkelsen.DESCRIPTION OF THE RIS PUFF DIFFUSION MODEL.《Nuclear Technology》.1984,第67卷(第01期),第56-65页. *
刘中坤 等.核设施退役虚拟仿真***框架研究.《原子能科学技术》.2011,第45卷(第09期),第1080-1086页. *
王鹏飞 等.气旋活动对福岛核污染物扩散影响的模拟研究.《环境科学研究》.2013,第26卷(第01期),第50-56页. *
陈祖刚 等.基于ArcEngine高斯烟团气体扩散模型的模拟实现.《测绘信息与工程》.2011,第36卷(第02期),第10-12,36页. *

Also Published As

Publication number Publication date
CN110489789A (zh) 2019-11-22

Similar Documents

Publication Publication Date Title
CN103370712A (zh) 用于机架冷却分析的***和方法
CN110826188B (zh) 一种基于gpu加速的天然气管网水力参数仿真方法
CN111881630B (zh) 一种核设施近场环境模拟数值风洞的建立***及方法
Teitel et al. Air exchange and ventilation efficiencies of a monospan greenhouse with one inflow and one outflow through longitudinal side openings
CN110489789B (zh) 一种核设施退役环境中放射性气体扩散评估方法
Chen et al. Numerical study on the integrated effects of supplied air velocity and exhaust velocity on particles removal for industrial buildings
CN114330050A (zh) 抽水蓄能电站地下厂房通风廊道的温湿度数值模拟方法
CN105426984A (zh) 基于多粒子群优化算法的给水管网传感器布置优化方法
Jin et al. Numerical investigation of particle transport characteristics in an isolated room with single-sided natural ventilation
KR101788389B1 (ko) 가스 확산 모델링 장치 및 방법
Shao et al. Long-term prediction of dynamic distribution of passive contaminant in complex recirculating ventilation system
Cybyk et al. Simulation of fluid dynamics around complex urban geometries
CN110472279B (zh) 一种基于涡动扩散模型的放射性气体浓度评估方法
CN113657049B (zh) 一种池式钠冷快堆主冷却剂***传热与流动快速仿真方法
Jin et al. Hybrid Strategy Improved Grey Wolf Optimization Algorithm for Plume Tracking and Localization Method in Indoor Weak Wind Environment
Bi et al. Simulation study of effective pumping speed of divertor pumping system for CFETR with COMSOL Multiphysics
CN112634114A (zh) 一种辐射环境管理与评价方法及***
JP2009193551A (ja) 建造物内における火災及び有毒物質拡散のシミュレーション技術
Fuestel et al. The COMIS infiltration model
Cai et al. Rapid identification of single constant contaminant source by considering characteristics of real sensors
Yang et al. Radioactive gas diffusion simulation and inhaled effective dose evaluation during nuclear decommissioning
Abouelhamd et al. Effect of crowd density, wind direction, and air temperature on the formation of individual human breathing zones in a semi-outdoor environment
Yu et al. Assessment of the Revised 3410 building filtered exhaust stack sampling probe location
CN219233465U (zh) 一种实验室多通道导流风道
Wang et al. Construction and Performance Analysis of a Flow Field Optimization Model for Denitrification System

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant