放射性气体扩散的预估模型

基于倾斜烟团模型的放射性气体预估方法

摘要

日本震后核泄露造成世界范围恐慌，通过合理建模可以有效的预测放射性气体的扩散情况，为预防核辐射提供科学指导。为解决这一问题，我们对核事故中放射性云团在大气中的扩散规律进行了研究，在高斯烟羽模型的基础上进行了改进，考虑实际过程中核素粒子的重力沉降、风速变化等因素的影响，得到了改进后的气体扩散模型--倾斜烟团模型，并给出一种迅速估算放射性核素扩散浓度的新方法。该方法可计算核事故中连续点源和瞬时点源在不同气象、地形条件下的浓度分布，与典型的高斯模型比较，其计算结果在10m的范围内总体合理，但在近源处有所偏小。经过验证，该模型具有很好的实用性，在核事故的应急救援过程中可为制定救援方案提供科学依据。

关键字：核泄漏高斯烟羽倾斜烟团

一、问题提出

有一座核电站遇自然灾害发生泄漏，浓度为p的放射性气体以匀速排出，速度为m kg/s，在无风的情况下，匀速在大气中向四周扩散, 速度为s m/s. 现在，就此情况，请搜集数据，建立数学模型，回答以下问题：

1）请你建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。

2）当风速为k m/s时，给出核电站周边放射性物质浓度的变化情况。

3）当风速为k m/s时，分别给出上风和下风L公里处，放射性物质浓度的预测模型。

4）将你建立的模型应用于福岛核电站的泄漏，计算出福岛核电站的泄漏对我国东海岸的影响。

二、问题分析

目前，应用于核事故应急决策系统中的大气扩散模型可分为3类，即高斯烟羽模型，拉格朗日烟团模型和三维数值计算模型。实际应用较多的为高斯烟羽模型，如直线高斯烟羽模型和变轨迹高斯烟羽模型。但是高斯模式适用于比较均匀稳定的流动条件，而且没有考虑其他因素的影响，在实际使用时模拟结果的可靠性往往较低。烟团模型和粒子模型能克服高斯烟羽模型的局限性，在复杂条件下可以较为真实地模拟气载放射性物质在大气传输中的时空分布。由于三维数值计算模型要求较大的计算机花费较长的计算时间，在实时的核事故应急响应中往往不能满足要求。我们在考虑了实际地形、地貌对放射性污染物在大气中输运扩散的影响后，根据实时的各种风场、气候潮湿程度的预报，分析放射性粒子的沉降、放射性气体的空间分布和时间分布，通过倾斜烟团模型来预测几小时内至未来数日内的放射性核素传播的区域和数量，为核泄漏后实施各种应急计划和措施提供了科学的依据。

三、模型假设与符号假设

3.1 符号假设

kg；

,

,

x处的浓度，3m

C——空间点)

y

(z

y

,

)

x

,

,

z

(H

m——泄漏源强，s

kg；

k

——环境风速，s m ；

k ——平均环境风速，s m ；

s ——放射性物质扩散速度，s m ；

s V ——沉降速度，s m ，一般含碘放射性核素的干沉积速度为s V =1.1s cm ；

D ——粒子直径，m ；

g

——重力加速度，9.80652s m ；

ρ——粒子密度，3

m kg ；

μ——空气的动力粘性系数，可取)(108.15

s m kg ∙⨯-；

d

——出口直径，m ；

x σ——下风向扩散系数，m ；

y σ——侧风向扩散系数，m ；

z σ——垂直风向扩散系数，m ；

H

——有效源高，m ，它等于泄漏源高度与抬升高度之和, 即H =0H + ∆H ；

H

——泄漏源高度，m ；

∆H ——泄漏源抬升高度，m 。 3.2 高斯模型基本假设

3.2.1 模型的基本表示及相应假设

云羽中心线

风速k

横风向距离y （m

m ）

垂直距离z 释放高度抬升高度

图1 高斯模型示意图

高斯模型用来描述危险物质泄漏形成的非重气云扩散行为, 或描述重气云在

重力作用消失后的远场扩散行为。为了便于分析, 建立如下坐标系OXYZ: 其中原点O 是泄漏点在地面上的正投影，X 轴沿下风向水平延伸，Y 轴在水平面上垂直于X 轴，Z 轴垂直向上延伸，如图1所示。高斯模型除了须满足平板模型中的一般假设外，还需要满足如下假设: ①气云在平整、无障碍物的地面上空扩散；②气云中不发生化学反应和相变反应，也不发生液滴沉降现象；③气体泄漏速率不随时间变化，即为稳定泄漏；④风向为水平方向，风速和风向不随时间、地点和高度变化，云团中心的移动速度或云羽轴向蔓延速度等于环境风速；⑤气云和环境之间无热量交换。⑥云团内部或云羽横截面上浓度、密度等参数服从高斯分布(即正态分布)。

3.2.2 高斯模型计算公式

根据高斯模型，以泄漏源为原点，风向方向为轴的空间坐标系中一点(x, y, z)处的浓度为:

)2)(2)((22),,,(2

2

222

2z z y y

x H z e H z e y e

k m H z y x C σσπ++-⨯⨯∙=

-

-

-

（1）

(高斯模型的C 语言实现可见附录1)

四、模型建立与求解

4.1 对高斯模型的修正

利用标准的高斯模型，计算大气中放射性核素云团的扩散，目前已有较多的研究，而实际的核素扩散过程还存在着粒子的干沉积(重力沉降)、湿沉积(雨洗)以及核素衰变等的因素对浓度分布的影响。因此，有必要对高斯浓度扩散模式进行修正，方可较为准确、真实地反映实际的核素扩散规律。由于地面风速对大气扩散的影响起着至关重要的作用，不同的风速其浓度计算方法有较大差异，因此，下面的分析将主要考虑风速的影响。 4.2 考虑干沉积时的连续点源扩散

粒径大于10肿的粒子有明显的重力沉降，粒子的沉降速度取决于空气阻力和重力平衡，可用斯托克斯公式表示：μ

ρ182

gD V s =

。

我们采用部分反射的倾斜烟云模式来计算干沉积的影响，其基本的思路是：因为在扩散过程中同时有重力沉降的位移迭加到羽流中心线上，中心线就会向下倾斜，所有粒子相当于在下倾的中心线上扩散。该类扩散和沉降的迭加可认为是