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

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

Info

Publication number
CN110489789A
CN110489789A CN201910617663.0A CN201910617663A CN110489789A CN 110489789 A CN110489789 A CN 110489789A CN 201910617663 A CN201910617663 A CN 201910617663A CN 110489789 A CN110489789 A CN 110489789A
Authority
CN
China
Prior art keywords
radgas
cigarette
model
elimination
nuclear facilities
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.)
Granted
Application number
CN201910617663.0A
Other languages
English (en)
Other versions
CN110489789B (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

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)是球形烟团模型的半径,可根据放射性气体 的扩散速度获得,球形烟团模型内放射性气体的浓度为:
式中,Ci,t(j)为t时刻第i个泄露点释放的第j个烟团内放射性气体的浓度,Bq/m3;Qi(j) 为第i个泄露点释放第j个烟团时泄露点处放射性气体的释放率,Bq/s;△t为泄露点释放相邻 气团的时间间隔,s;Vi,t(j)是t时刻第i个泄露点释放的第j个烟团的体积,单位是m3,理想 情况下,烟团模型的体积为:
3、所述步骤3具体为:考虑放射性气体烟团形状受到厂房空间的限制,将均匀球形烟团 模型的形状公式公式修正为:
式中,xmin、ymin、zmin、xmax、ymax、zmax为厂房空间的边界,烟团的体积修正后为:
其中,Vout为烟团在厂房边界外的体积。
4、所述步骤4具体为:考虑放射性气体烟团模型与场景的数学模型差运算,对烟团形状 公式进行修正,修正后烟团的体积为:
其中,Vin为烟团与核设施相交部分的体积。
5、所述步骤5具体为:由于核退役场景中放射性气体泄露点多,且每个烟团都是独立的, 因此各空间位置处放射性气体的浓度,为所有泄露点产生的包含该位置烟团的浓度的累加和。
6、所述步骤6具体为:当烟团与排风管道相接时,烟团模型形状不变,单位时间内通过 排风管道排走的气体活度为烟团的浓度乘以排风管的气体流量,修正后烟团的浓度为:
其中,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)是球形烟团模型的半径,可根据放射性气体的 扩散速度获得。
根据假设条件,烟团内放射性气体的浓度是均匀分布的,因此球形烟团模型内放射性气 体的浓度为:
式中,Ci,t(j)为t时刻第i个泄露点释放的第j个烟团内放射性气体的浓度,Bq/m3;Qi(j)为第 i个泄露点释放第j个烟团时泄露点处放射性气体的释放率,Bq/s;△t为泄露点释放相邻气团 的时间间隔,s;Vi,t(j)是t时刻第i个泄露点释放的第j个烟团的体积,m3。在理想情况下, 均匀球形烟团模型为圆球形,体积为:
3、考虑厂房对模型的影响。
放射性气体扩散形状会受到厂房空间的限制,因此需要对烟团形状公式进行修正:
式中,xmin、ymin、zmin、xmax、ymax、zmax为厂房空间的边界。如图2所示,烟团的体积修正后 为:
其中,Vout为烟团在厂房边界外的体积。
4、考虑退役场景对模型的影响。
核退役场景中的设施体积巨大,会影响放射性气体的扩散。因此需要对烟团模型与场景 的数学模型进行差运算,对烟团形状公式进行修正。如图3所示,烟团的体积修正后为:
其中,Vin为烟团与核设施相交部分的体积。
5、考虑多个泄露点对模型的影响。
核退役场景中,放射性气体泄露点多,且每个泄露点的泄露量也不相同。如图4所示, 由于每个烟团都是独立的,因此(x,y,z)位置处的放射性气体的浓度,为所有泄露点产生的包 含该位置烟团的浓度的累加和。
6、考虑排风管道对模型的影响。
放射性气体会通过排风管道,经过滤、吸收等操作后排放到外界环境。如图5所示,当 烟团与排风管道相接时,假设模型形状不变,则单位时间内通过排风管道排走的气体活度为 烟团的浓度乘以排风管的气体流量。修正后烟团的浓度为:
其中,Vf为单位时间内排风管道抽走气体的体积。
综上,本发明公开一种核设施退役环境中放射性气体扩散评估方法,具体涉及一种针对 核退役环境中泄露的放射性气体,采用均匀球形烟团模型模拟放射性气体的扩散,并动态的 计算放射性气体浓度分布的仿真方法。本发明包括:构建退役场景数学模型;采用均匀球形 烟团模型模拟放射性气体的扩散;使用退役场景的数学模型对均匀球形烟团模型进行修正; 计算放射性气体的活度分布,实现退役环境中放射性气体的扩散评估。本发明包括退役环境 建模、均匀球形烟团模型建模、放射性气体浓度分布计算三个模块,实现了核设施退役环境 中放射性气体扩散的动态计算。

Claims (7)

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

Cited By (1)

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

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030213315A1 (en) * 2002-04-12 2003-11-20 New Millennium Nuclear Technologies, Llc Selective profiling of bulk building material for characterization using a concrete contamination profiling technology
CN107526908A (zh) * 2016-06-21 2017-12-29 中国辐射防护研究院 核事故场外后果评价中拉格朗日烟团大气扩散模拟方法
CN107526852A (zh) * 2016-06-21 2017-12-29 中国辐射防护研究院 一种核设施事故场外后果实时在线评价方法及***
CN108563839A (zh) * 2018-03-23 2018-09-21 哈尔滨工程大学 一种核设施退役模型程式化仿真方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030213315A1 (en) * 2002-04-12 2003-11-20 New Millennium Nuclear Technologies, Llc Selective profiling of bulk building material for characterization using a concrete contamination profiling technology
CN107526908A (zh) * 2016-06-21 2017-12-29 中国辐射防护研究院 核事故场外后果评价中拉格朗日烟团大气扩散模拟方法
CN107526852A (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》, vol. 06, 30 April 2004 (2004-04-30), pages 271 - 277 *
TORBEN MIKKELSEN: "DESCRIPTION OF THE RIS PUFF DIFFUSION MODEL", 《NUCLEAR TECHNOLOGY》, vol. 67, no. 01, 31 October 1984 (1984-10-31), pages 56 - 65 *
刘中坤 等: "核设施退役虚拟仿真***框架研究", 《原子能科学技术》, vol. 45, no. 09, 20 September 2011 (2011-09-20), pages 1080 - 1086 *
王鹏飞 等: "气旋活动对福岛核污染物扩散影响的模拟研究", 《环境科学研究》, vol. 26, no. 01, 15 January 2013 (2013-01-15), pages 50 - 56 *
陈祖刚 等: "基于ArcEngine高斯烟团气体扩散模型的模拟实现", 《测绘信息与工程》, vol. 36, no. 02, 5 April 2011 (2011-04-05), pages 10 - 12 *

Cited By (2)

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

Also Published As

Publication number Publication date
CN110489789B (zh) 2023-04-18

Similar Documents

Publication Publication Date Title
Hong et al. Validation of an open source CFD code to simulate natural ventilation for agricultural buildings
CN106021817B (zh) 一种海相气田高含硫天然气集输装置泄漏快速模拟方法
Xia et al. CFD analysis of pre-cooling water spray system in natural draft dry cooling towers
Feng et al. An improved particle swarm optimization method for locating time-varying indoor particle sources
CN102967689B (zh) 一种污染源监测布点方法
CN103914622A (zh) 一种化学品泄漏快速预测预警应急响应决策方法
CN111881630B (zh) 一种核设施近场环境模拟数值风洞的建立***及方法
CN104361157B (zh) 一种建筑物间风环境的评价方法
CN107526908A (zh) 核事故场外后果评价中拉格朗日烟团大气扩散模拟方法
CN107145699A (zh) 气载放射性核素长距离迁移拉格朗日粒子扩散计算方法
Ruiz et al. Experimental study of drift deposition from mechanical draft cooling towers in urban environments
CN110489789A (zh) 一种核设施退役环境中放射性气体扩散评估方法
CN111611698B (zh) 一种超薄体传热仿真方法
Michioka et al. Large-eddy simulation of the effects of wind-direction fluctuations on turbulent flow and gas dispersion within a cubical canopy
CN108629056A (zh) 一种风电功率预测方法及***
Patnaik et al. Large scale urban simulations with the MILES approach
Oaks et al. Eulerian numerical modeling of contaminant transport in Lower Manhattan, New York City, from a point-source release under the dominant wind condition: Insights gained via LES
Cybyk et al. Simulation of fluid dynamics around complex urban geometries
CN105808812A (zh) 一种地表水水龄二维介观数值模拟方法
Yang et al. Numerical study on the smoke flow characterization and phenomenon of plug-holing under lateral smoke exhaust in tunnel fire
Patnaik et al. Large scale urban simulations with FCT
Michioka et al. Effects of fetch on turbulent flow and pollutant dispersion within a cubical canopy
CN108846163A (zh) 一种用于确定安全壳试验前气相初始状态的方法
CN110472279B (zh) 一种基于涡动扩散模型的放射性气体浓度评估方法
Li et al. A new method for simulating aerosols Brownian coagulation based on finite active samples assumption

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