CN107133412B - 一种河道中多组分放射性核素浓度预报的快速计算方法 - Google Patents

一种河道中多组分放射性核素浓度预报的快速计算方法 Download PDF

Info

Publication number
CN107133412B
CN107133412B CN201710341976.9A CN201710341976A CN107133412B CN 107133412 B CN107133412 B CN 107133412B CN 201710341976 A CN201710341976 A CN 201710341976A CN 107133412 B CN107133412 B CN 107133412B
Authority
CN
China
Prior art keywords
river
concentration
radionuclide
multicomponent
equation
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
CN201710341976.9A
Other languages
English (en)
Other versions
CN107133412A (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.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
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 China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN201710341976.9A priority Critical patent/CN107133412B/zh
Publication of CN107133412A publication Critical patent/CN107133412A/zh
Application granted granted Critical
Publication of CN107133412B publication Critical patent/CN107133412B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Molecular Biology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种河道中多组分放射性核素浓度预报的快速计算方法,包括;首先利用数值算法计算不考虑衰变情形下河道中核素的浓度分布,并利用解析法将考虑衰变作用的核素浓度同不考虑衰变作用的核素浓度建立转换关系,从而快速获得河道中多组分放射性核素的浓度分布。该河道中多组分放射性核素浓度预报的快速计算方法计算参数获取精准并且计算速度快的,其能够通过对河道水深、二维向流速及相应各系数的获取构建具有普适性的计算方式及模型,得到河道中多组分放射性核素的浓度分布。本发明的计算速度快、精度高,可适用于不同河道的计算模型构建。

Description

一种河道中多组分放射性核素浓度预报的快速计算方法
技术领域
本发明属于核环境检测技术,具体的涉及一种评价突发事故状态下不同半衰期核素的环境影响,对河道中多组分放射性核素浓度进行预报的快速计算方法。
背景技术
近年来,为保障能源安全、优化能源结构,国家制定了积极发展核电的能源政策,提出要在大型先进压水堆核电站示范工程、高温气冷堆核电机组商业运行方面取得突破。核电站在突发事故期间会向海洋、河流、湖泊、水库等环境水域排放放射性核素。核素在衰变过程中会发出α射线、β射线、γ射线以及X射线,从而对接触核素的人体产生不同程度的危害。在大剂量照射情况下,神经系统、免疫系统、造血系统、生殖系统和消化系统会受到严重损伤,出现细胞蛋白质凝固、新陈代谢终止、甚至大量细胞立即死亡等危及生命安全的严重后果。与放射性核素在海域中输运相比,其在河道中输运时,受流域水量的限制,影响范围更大,持续时间更长,安全性问题更突出。日本福岛核电事故所引发的放射性核素泄漏使得核电站周围80公里范围内的公众受到超过安全阈值的辐射剂量的照射,20万居民被迫从居住地撤离,对生态环境的影响将长达数十年。为保障流域的水安全,最大限度地发挥流域的生态服务功能,必须对突发事故状态下核素的影响范围进行快速预报,以便采取科学合理的应对措施,而河道中多组分放射性核素浓度的快速算法正是相关部门有效开展这些工作的重要技术支撑,对核电水安全评价与应对具有重要意义。
核电站排出的放射性核素种类较多,达数十种,比如氚、铯-137、锶-90、钴-60、碘-131等。不同核素的半衰期差别较大,比如碘-131的半衰期为8天,钴-60的半衰期为5.27年,氚的半衰期12.43年,锶-90的半衰期为28年,铯-137的半衰期30.17年。因此,对于同一河道,不同组分核素的影响程度存在较大差异。为了评价突发事故状态下不同半衰期核素的环境影响,就需要快速预报多组分放射性核素的浓度分布。对于实际河道,由于地形与岸线十分复杂,已有的算法是采用纯数值理论,比如有限差理论、有限元理论、有限体积理论,分别计算每种核素的浓度。这类算法的不足在于需要对每种核素所满足的控制方程分别进行数值离散与求解。如果核素种类较多,计算时间会变得很长,难以满足突发性核事故情况下快速预测的要求。因此,探索一种快速计算多组分放射性核素浓度的算法是核电水安全预报领域亟待解决的问题,对核电水环境风险评价具有十分重要的意义。
发明内容
本发明提供了一种计算参数获取精准并且计算速度快的河道中多组分放射性核素浓度预报的快速计算方法,其能够通过对河道水深、二维向流速及相应各系数的获取构建具有普适性的计算方式及模型,得到河道中多组分放射性核素的浓度分布。本发明的计算速度快、精度高,可适用于不同河道的计算模型构建。
本发明所采用的技术方案如下:
一种河道中多组分放射性核素浓度预报的快速计算方法,其特征在于所述计算方法包括;
首先利用数值算法计算不考虑衰变情形下河道中核素的浓度分布,并利用解析法将考虑衰变作用的核素浓度同不考虑衰变作用的核素浓度建立转换关系,从而快速获得河道中多组分放射性核素的浓度分布。
具体的讲,所述计算方法进一步包括:
进行核素的浓度分布计算时,首先采用数值算法求解连续性方程、二维深度平均的水动力学方程,获得河道的流速分布;所述数字算法包括有限差分法、有限元法或者有限体积法。
另外,所述计算方法进一步包括:
所述连续性方程为:
二维平均的水动力学方程包括:
x方向动量方程:
y方向动量方程:
其中:h为水深,t为时间变量,为梯度算子,Z为相对基准面水位,u1、u2分别为x、y方向的垂向平均流速,u为流速矢量,g为重力加速度,Sx、Sy分别为动量方程中x、y方向的源项,vt为紊动粘性系数。
一实施方式中,河道上游采用流量边界条件,河道下游采用水位边界条件,固壁采用自由滑移边界条件。初始条件采用静流条件即初始速度设置为0起算。
另一实施方式中,所述计算方法进一步包括:
进行核素的浓度分布计算时,其次在不考虑核素衰变作用,计算河道中核素浓度分布;
所述不考虑衰变情形下,核素浓度的控制方程、边界条件、初始条件为:
C0(x,y,0)=Iδ(x-x0,y-y0) (6)
其中:C0为核素浓度,vT为紊动扩散系数,I为突发事故状态下瞬间排放的核素强度,δ为狄拉克δ函数,(x0,y0)为排放点位置。
再一实施方式中,所述计算方法进一步包括:
进行核素的浓度分布计算时,最后通过数值算法求得无衰变核素浓度C0后,通过如下变换公式计算衰变速率为k(量纲为1/[时间量纲])的核素的浓度分布:
Ck(x,y,t)=C0(x,y,t)exp(-kt) (7)
其中K的量纲为1/时间量纲。
另外,所述设定水文条件为恒定流条件,h、u、vT只是空间位置(x,y)的函数。
本发明采用数值算法与解析方法相结合的思路,其预报河道中多组分放射性核素浓度的快速算法原理为:利用数值算法计算不考虑衰变情形下河道中核素的浓度分布,并利用解析法将考虑衰变作用的核素浓度同不考虑衰变作用的核素浓度建立转换关系,从而快速获得河道中多组分放射性核素的浓度分布。
本发明的有效果在于,该河道中多组分放射性核素浓度预报的快速计算方法计算参数获取精准并且计算速度快的,其能够通过对河道水深、二维向流速及相应各系数的获取构建具有普适性的计算方式及模型,得到河道中多组分放射性核素的浓度分布。本发明的计算速度快、精度高,可适用于不同河道的计算模型构建。
下面结合具体实施方式对本发明做进一步的阐述。
具体实施方式
本发明预报河道中多组分放射性核素浓度的快速算法原理为:利用数值算法计算不考虑衰变情形下河道中核素的浓度分布,并利用解析法将考虑衰变作用的核素浓度同不考虑衰变作用的核素浓度建立转换关系,从而快速获得河道中多组分放射性核素的浓度分布。
具体算法过程如下:
首先,采用常规的数值算法,例如有限差分法、有限元法、有限体积法求解连续性方程、二维深度平均的水动力学方程,获得河道的流速分布。控制方程与边界条件如下:
连续性方程:
x方向动量方程:
y方向动量方程:
其中:h为水深,t为时间变量,为梯度算子,Z为相对基准面水位,u1、u2分别为x、y方向的垂向平均流速,u为流速矢量,g为重力加速度,Sx、Sy分别为动量方程中x、y方向的源项,vt为紊动粘性系数。
河道上游采用流量边界条件,河道下游采用水位边界条件,固壁采用自由滑移边界条件。初始条件采用静流条件起算,即初始速度设置为0。
其次,不考虑核素衰变作用,计算河道中核素浓度分布。不考虑衰变情形下,核素浓度的控制方程、边界条件、初始条件为:
C0(x,y,0)=Iδ(x-x0,y-y0) (6)
其中:C0为核素浓度,vT为紊动扩散系数,I为突发事故状态下瞬间排放的核素强度,δ为狄拉克δ函数,(x0,y0)为排放点位置。
第三步,通过数值算法求得无衰变核素浓度C0后,通过如下变换公式计算衰变速率为k(量纲为1/[时间量纲])的核素的浓度分布:
Ck(x,y,t)=C0(x,y,t)exp(-kt) (7)
变换公式(7)的推导过程如下:
对于河道中核素浓度的水环境风险评价,通常考虑的水文条件为恒定流条件,因此h、u、vT只是空间位置(x,y)的函数。
衰变速率为k的核素的浓度控制方程、边界条件以及初始条件为:
Ck(x,y,0)=Iδ(x-x0,y-y0) (10)
令Ck(x,y,t)=Ω(x,y,t)exp(-kt),将其带入方程式(6)-(8)式并利用h、u、vT只是空间位置(x,y)的函数,可得:
Ω(x,y,0)=Iδ(x-x0,y-y0) (13)
方程式(9)-(11)式表明令Ω(x,y,t)即为不考虑衰变的核素浓度C0(x,y,t),可得出方程式(7)式成立。

Claims (3)

1.一种河道中多组分放射性核素浓度预报的快速计算方法,其特征在于所述计算方法包括:
首先,利用有限差分法求解连续性方程、二维深度平均的水动力学方程,获得河道的流速分布,
所述连续性方程为:
所述二维深度平均的水动力学方程包括:
x方向动量方程:
y方向动量方程:
其中:h为水深,t为时间变量,▽为梯度算子,Z为相对基准面水位,u1、u2分别为x、y方向的垂向平均流速,u为流速矢量,g为重力加速度,Sx、Sy分别为动量方程中x、y方向的源项,νt为紊动粘性系数;
其次,利用有限差分法计算不考虑核素衰变情形下河道中放射性核素的浓度C0的分布,
所述不考虑衰变情形下,核素浓度的控制方程、边界条件、初始条件为:
n·▽C0=0 (5)
C0(x,y,0)=Iδ(x-x0,y-y0) (6)
其中:C0为核素浓度,νT为紊动扩散系数,n为计算域边界指向计算域外的单位法线向量,I为突发事故状态下瞬间排放的核素强度,δ为狄拉克δ函数,(x0,y0)为排放点位置;
第三,通过如下变换公式计算考虑核素衰变情形下核素的浓度分布:
Ck(x,y,t)=C0(x,y,t)exp(-kt) (7)
其中:x、y为核素的空间位置坐标,k为衰变速率(1/时间量纲),t为时间变量。
2.根据权利要求1所述的河道中多组分放射性核素浓度预报的快速计算方法,其特征在于河道上游采用流量边界条件,河道下游采用水位边界条件,固壁采用自由滑移边界条件,初始条件采用静流条件即初始速度设置为0起算。
3.根据权利要求1所述的河道中多组分放射性核素浓度预报的快速计算方法,其特征在于设定水文条件为恒定流条件,h、u、νT只是空间位置(x,y)的函数。
CN201710341976.9A 2017-05-15 2017-05-15 一种河道中多组分放射性核素浓度预报的快速计算方法 Active CN107133412B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710341976.9A CN107133412B (zh) 2017-05-15 2017-05-15 一种河道中多组分放射性核素浓度预报的快速计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710341976.9A CN107133412B (zh) 2017-05-15 2017-05-15 一种河道中多组分放射性核素浓度预报的快速计算方法

Publications (2)

Publication Number Publication Date
CN107133412A CN107133412A (zh) 2017-09-05
CN107133412B true CN107133412B (zh) 2018-08-31

Family

ID=59731781

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710341976.9A Active CN107133412B (zh) 2017-05-15 2017-05-15 一种河道中多组分放射性核素浓度预报的快速计算方法

Country Status (1)

Country Link
CN (1) CN107133412B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110111854B (zh) * 2019-04-01 2023-09-15 中国辐射防护研究院 放射性核素从地表水释放途径向水生生物浓集的计算方法
CN116401500B (zh) * 2023-04-06 2023-09-29 中国水利水电科学研究院 周期排放方式下长距离暗管内外水体核素浓度的计算方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6741955B2 (en) * 2002-03-05 2004-05-25 Uson, L.P. System and method for leak rate testing during adiabatic cooling
CN104182600B (zh) * 2013-05-21 2017-04-12 环境保护部核与辐射安全中心 核电厂气载放射性物质源空间活度浓度评价方法
CN104156630B (zh) * 2014-09-05 2017-06-23 西南科技大学 三维核素扩散的计算方法
CN104933294B (zh) * 2015-05-25 2017-10-17 苏州热工研究院有限公司 放射性物质迁移和传播评估模型构建方法
CN105468831B (zh) * 2015-11-19 2019-02-01 厦门大学 一种反应堆工程仿真机辐射值仿真方法

Also Published As

Publication number Publication date
CN107133412A (zh) 2017-09-05

Similar Documents

Publication Publication Date Title
CN106528926B (zh) 一种地下核电站烟羽应急计划区划分方法
CN103413062B (zh) 一种放射性核素扩散的计算方法
Margvelashvily et al. THREETOX-a computer code to simulate three-dimensional dispersion of radionuclides in stratified water bodies
Budyansky et al. Role of mesoscale eddies in transport of Fukushima-derived cesium isotopes in the ocean
CN107133412B (zh) 一种河道中多组分放射性核素浓度预报的快速计算方法
Lyubimova et al. The modeling of the formation of technogenic thermal pollution zones in large reservoirs
Shen et al. Shelter construction for fish at the confluence of a river to avoid the effects of total dissolved gas supersaturation
Zheleznyak et al. Hydrological dispersion module of JRODOS: renewed chain of the emergency response models of radionuclide dispersion through watersheds and rivers
Periáñez Redissolution and long-term transport of radionuclides released from a contaminated sediment: a numerical modelling study
Raza et al. A 3-D Lagrangian stochastic model for the meso-scale atmospheric dispersion applications
Johnson VAHM: A Vertically Averaged Hydrodynamic Model using boundary-fitted coordinates
Salomon et al. Mathematical model of 125Sb transport and dispersion in the Channel
Nie et al. Numerical study on tritium dispersion in coastal waters: The case of Hangzhou Bay, China
de Carvalho Gomes et al. Tritium (3H) as a tracer for monitoring the dispersion of conservative radionuclides discharged by the Angra dos Reis nuclear power plants in the Piraquara de Fora Bay, Brazil
KR101724396B1 (ko) 침몰사고 시 방사성폐기물 영향 평가 시스템 및 방법
Zheleznyak The mathematical modelling of radionuclide transport by surface water flow from the vicinity of the Chernobyl Nuclear Power Plant
Devkota et al. Response characteristics of the Perdido and Wolf Bay system to inflows and sea level rise.
Suh et al. Estimation of the radionuclides emission region using trajectory and atmospheric dispersion models
Deleersnijder et al. A simple model of the tracer flux from the Mururoa lagoon to the Pacific
Xin et al. Reservoir operation schemes for water pollution accidents in Yangtze River
Guan et al. The numerical simulation of caesium-137 transportation in ocean and the assessment of its radioactive impacts after Fukushima NPP release
CN205679771U (zh) 一种走航式水体放射性巡测与展示系统
Margvelashvili et al. 3-D numerical modelling of mud and radionuclide transport in the Chernobyl Cooling Pond and Dnieper-Boog Estuary
Jeong et al. Radiological risk assessment for an urban area: Focusing on a drinking water contamination
Aly et al. Dispersion modeling: methodology and new trends

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