基于MATLAB的放射性气体扩散模拟及应用

基于MATLAB的放射性气体扩散模拟及应用

发表时间：2020-02-27T22:44:59.213Z 来源：《防护工程》2019年19期作者：杨震何维锋冷宇飞邢佳奇周文平（通讯作者）[导读] 随着核电在世界范围内广泛应用，越来越多的核电站建立起来，但是，运用核能的同时我们也要防范于未然。

杨震何维锋冷宇飞邢佳奇周文平（通讯作者）沈阳工程学院核工程与核技术沈阳 110121 摘要：随着核电在世界范围内广泛应用，越来越多的核电站建立起来，但是，运用核能的同时我们也要防范于未然。因此，模拟核事故放射性气体扩散有着非常重要的意义。本文运用概率动力学的相关知识，结合高斯扩散模型，建立放射性气体扩散模型，模拟了福岛核电站放射性气体扩散对我国东海岸的影响。

关键词：核事故放射性气体扩散；高斯扩散模型；matlab放射性气体扩散模型

1.前言

北京时间2011年3月11日，日本福岛县的福岛第一核电站发生了一起重大核事故，大量的放射性污染气体从事故的核电站泄露进入大气，对大气环境产生了非常严重的污染，短时间内事故等级从四级跃升到最高等级——七级核事故，引起了国际社会的广泛关注。3月15日，专家组分析相关数据得出较低浓度的放射性气体正从核电站向福岛以东地区扩散，并可能在将来几天内到达北美地区，最终到达欧洲地区。同时事故核电站10km范围内的所有居民被日本政府要求紧急撤离。核电站周围的各个监测站检测到碘，氩，钚等多种放射性同位素从核电站泄出，23日，在核电站厂区内检测出中子辐射。随着核电在世界范围内广泛应用，越来越多的核电站建立起来，但是，运用核能的同时我们也要防范于未然。因此，模拟核事故放射性气体扩散有着非常重要的意义。本文运用概率动力学的相关知识，结合高斯扩散模型，建立放射性气体扩散模型，模拟了福岛核电站放射性气体扩散对我国东海岸的影响。

首先，我们考虑到风向和风速对放射性气体浓度分布有一定影响，由于当时日本发生了大地震和大海啸，这使得当地的气象环境十分复杂。本文我们结合高斯烟羽模型，并考虑烟气抬升，地面反射，干湿沉积，放射性气体的衰变等多种因素对模型进行反复修正，得到最终的模型。之后，对于上风向和下风向L公里处的放射性气体浓度，只需在上述的基础上，令x=L或-L，同时，将风速k用（k+s）和（k-s）来代替，y=0，z=0。建立完这个基础模型后，我们就可以以此为基础研究在风速一定的情况下，位处上下风L公里处，放射性气体浓度的估计模型，并用matlab软件运行模型，模拟出在下风向时的浓度分布图。最后，我们参阅整理了气象、地理、新闻等资料，确定大气稳定度的等级为D 级，并且我们还查阅了在风速一定情况下，位于上下风L公里处的放射性气体浓度估计模型中所需要的各种参数。我国东海岸位于日本福岛的上风方向，所以，我们可以利用上风向L公里处的浓度预估模型，来预测放射性气体到达我国东海岸时的浓度。同时使用matlab 计算得出I的浓度为0.14，所得数据与实际测量值0.10基本吻合。

对于一个地区的核事故空气传播影响的分析预测，我们应参考大量气象资料、地理条件、新闻资料，确定大气稳定度等级，气体扩散系数，地表的平均风速、地理距离、当地地表对放射性物质的反射系数、空气黏度、放射性物质的半衰期等各种有用参数，并对其进行严谨的分析设计。事故发生之后我国也是立刻做出反应紧急应对，时刻关注事故对我国造成的影响。在我国东北部地区的黑龙江省某个气象监测站对空气样品进行检验时发现了来自日本福岛核电站的放射性核同位素I-131。统计我国31个省，市的大气放射性浓度记录数据，发现到达我国I-131浓度最大值位于我国东北吉林省，这与我国当时东部沿海地区的持续东风有关。不仅仅是我国，此次福岛核电站的重大事故也引起了国际社会的高度关注，美国，俄罗斯，澳大利亚，新加坡等日本周边国家都对核泄漏造成的空气中放射性气体浓度进行了紧密的持续经跟踪测量，一方面可以第一时间了解辐射污染物的准确数据，令国民提前做好预防工作，确保本国国民的人身安全；另一方面我们都了解核能是目前世界上众所周知的高效清洁能源，核电站的建设为各国带来的巨大的经济效益同时又减少了二氧化碳等温室气体的排放。核电站已经逐渐替代了火电站。但同时也没有人不了解原子弹的恐怖以及核泄漏的危害，核能的使用究竟是利大于弊还是弊大于利值得各国思考，毕竟人类已经有了切尔诺贝利和三里岛核电站事故的可怕阴影。因此我们小组成员以此次福岛核事故为例预测我国东海岸的放射性气体浓度，并与实际测量值进行比较。希望能够为我国对核事故的监测能力献出一份力。

2.模型假设与约定

1、假设连续云团的泄露时间和速度恒定，连续且均匀，初始泄露时刻连续云团内部的浓度、温度呈均匀的分布；

2、假设气体在扩散过程中只发生衰变，不发生分解以及任何化学反应等。

3、污染物的数量和浓度在y、z时间轴上的浓度分布通常是高斯分布（正态高斯分布）的；

4、扩散的处理过程中不能只考虑热量扩散后的云团内部外界温度的剧烈变化，忽略了热量的传递、热对流及其他外部热辐射；

5、泄漏的气体是我们所认为的理想气体，这些气体遵守理想气体的状态方程；

6、假设地面水平，且风向与地面保持水平；

7、在相对垂直和横向水平大气运动不同方向，大气扩散系数通常呈各向不相同性；

3.符号说明及名词解释 3.1符号说明：见下表1

表1 参数说明表

3.2 名词解释

烟羽：指烟气是泛指从一个烟囱中连续地将气排放出再进入输送到其他气体大气中的污染物气体中的一种烟雾或气流。由于大型烟羽各个组成部分的飞行运动和降落速度不同，因而其的整体外形也千变万化。不同的排放烟羽浓度大小和烟囱形状可以表示不同大气中的污染物平均排放烟气浓度的主要来源空间影响因素及其分布不同。它与其周围大气的变化湍流、大气流的稳定度、地形的湿地物、排放环境参数等诸多因素都具有密切的相互关系。

烟气抬升：暖一些暖气流体因受锋地表面、辐射耦合和冷气流的强大空气外力作用被迫进行上升和抬，或者在其高空运行中部分气体空气外力受强大地形的外力阻挡从而产生了强迫上升的飞行运动

干沉积：是主要定义是广泛指一些具有颗粒性的物体在通过之时受到强大重力的干扰作用或与其它具有颗粒性的物体之间发生强烈碰撞后继续运动时所可能发生的堆积沉降。

4.模型建立与求解

4.1风速对放射性物质浓度分布影响的高斯模型的建立与求解

4.1.1模型的建立

该模型运用了概率动力学的知识，通过图解和数学的推导得出连续点源周边的放射性化学物质高斯扩散的模型。

下面进行模型公式推导

由正态喷射气体中的浓度能量分布的公式假设正态气体中的浓度分布可以被通过推导而得出下式为风向任意一点处处当泄漏正态气体