基于高斯模型的放射性物质扩散模型

2011重庆邮电大学数学建模第二次模拟

承诺书

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则 •

我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮 件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问 题。

我们知道，抄袭别人的成果是违反竞赛规则的，如果引用别人的成果或其他 公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正 文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反 竞赛规则的行为，我们将受到严肃处理。

我们参赛选择的题号是（从 D/E/F中选择一项填写）： ____________________________ 我们的参赛报名号为（如果赛区设置报名号的话） ： ______________________________ 所属学校（请填写完整的全名）： _________________________________

参赛队员（打印并签名）：1. ____________________________________

2. ____________________________________ 3. __________________________________

指导教师或指导教师组负责人 （打印并签名）： ___________________________

日期：年 月曰

赛区评阅编号（由赛区组委会评阅前进行编号）：

2011重庆邮电大学数学建模第二次模拟

编号专用页

赛区评阅编号（由赛区组委会评阅前进行编号）：

赛区评阅记录（可供赛区评阅时使用）：

评阅人 评分 备注 全国统一编号（由赛区组委会送交全国前编号）：

全国评阅编号（由全国组委会评阅前进行编号）：

放射性气体扩散浓度预估模型

【摘要】本文是以日本地震引起的福岛核电站的核泄漏为背景， 并以给出的数据

为基础，研究某一假设核电站的核泄漏问题。我们通过收集相关的资料，并结合 题目给出的数据，建立了高斯模型、连续点源高斯扩散模型解决了题目提出的四 个问题。

针对问题一：考虑到泄漏源是连续、均匀和稳定的，我们运用散度、梯度、 流量等数学概念，通过“泄漏放射性物质质量守恒”、“气体泄漏连续性定理”、 Guass公式及积分中值定理得到了无界区域的抛物线型偏微分方程，然后再通过 电源函数解出空间任意一点的放射性物质浓度的表达式，

把此表达式定为模型一

的前身。鉴于放射性物质的扩散受到诸多因素的影响，如：泄漏源的实际高度、 地面反射等。我们以泄漏口为坐标原点建立三维坐标系，通过“像源法”处理地 面反射对放射性物质浓度的影响，并由此对模型一的前身进行修正完善， 得到模 型一：高斯模型，即放射性物质浓度的预测模型。 最后我们模拟了放射性物质无 风扩散仿真图。

针对问题二：当风速为k m/s时，我们根据放射性核素云团在大气中迁移和 扩散的数值计算的基本方法和步骤，并以泄漏点源在地面的投影点为坐标原点， 以风向方向为x轴，铅直方向为z轴，与x轴水平面垂直方向为y轴建立三维坐 标系，地面的反射作用同样利用“像源法”进行处理，得到连续点源高斯扩散模 型。考虑到地面反射、烟云抬升、放射性物质自身的沉降及雨水的吸附等对浓度 的影响，我们对连续点源高斯扩散模型进行了修正， 建立了修正的连续点源高斯

扩散模型。最后利用大气稳定度确定了扩散参数，进而求解了模型。

针对问题三：经分析，问题三的提出是以问题二为基础的，模型三的建立只 需要将模型二加以调整即可。我们以风速方向为 x轴正方向，将风速与放射性物 质的扩散速度进行矢量运算，此问题则转化为求

（L,0, z）和（-L,0, z）两点处的放

射性物质浓度，由此建立模型三，即上风和下风L公里处放射性物质浓度浓度的 预测模型。

针对问题四：首先，我们通过网络收集了相关数据，然后，我们结合模型二、 模型三对数据进行整理代入，算出了日本福岛核电站泄漏的放射性物质扩散到中

国东海岸和美国西海岸的浓度分别为 4.2429 10‘g/m3、2.3854 10‘g/m3。

关键词：高斯模型连续点源高斯扩散模型核泄漏

一问题的重述

1.1问题背景

目前，核电站的发展能带来巨大的经济效益和社会效益，但核电站一旦发生

核泄漏，将会给人们的生命健康和周边环境带来巨大的危害性影响。

2011年3

月日本的福岛核电站的放射性气体的核泄漏事件更让我们关注放射性气体泄漏 时的浓度问题。因此，正确的测出大气中放射性物质的浓度在环境监测和安全评 估中具有重要意义

1.2问题提出

有一座核电站遇自然灾害发生泄漏，浓度为 po的放射性气体以速度 m kg/s 匀速排出，在无风的情况下，以速度 s m/s匀速在大气中向四周扩散.

1） 在无风的情况下，建立一个描述核电站周边不同距离地区、

质浓度的预测模型。

不同时段放射性 物

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

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

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

西海岸的影响。

计算所用数据可以在网上搜索或根据具体情况自己模拟。

二符号说明

§(i =x,y,z) C(x,y, z,t) 空间任意一点的放射性物质的扩散系数 空间任意一点的放射性物质浓度 放射性气体的扩散速度,m/ s 泄漏源泄漏的放射性物质总量 空间域 空间域其体积 一规则的球面面积 在（lt+銚）内通过0的流量 s Qo Q V S Q1 Q2 Q内放射性物质的增量

t H 也H 任意扩散时刻 泄漏源距地面的实际高度 烟云抬升高度 泄漏源有效高度 泄漏源咼度处的平均风速，m/ s 源强，kg / s 用浓度标准差表示的x, y,z轴上的扩散参数 放射性气体出口流速,m/s 泄漏源出口的有效内径 泄漏源的热排放率，kw 泄漏源出口处温度，K 环境大气平均温度，K 沉降速度，m/s 空气的动力粘性系数 雨水吸附系数 太阳咼度角 h u Q crx, cry, crz Vs D QH Ts Ta 乂 a p ho 三模型假设

1、 扩散过程中浓度在y、z轴上的变化分布是高斯分布。

2、 放射性物质的扩散看作是空间某一连续点源向四周等强度地瞬时释放放射性 物质，放射性

物质在无穷空间扩散过程中不发生性质变化，且不计地形影响。

3、 放射性物质扩散服从扩散定律，即单位时间通过单位法向面积的流量与它的 浓度梯度成正

比。

4、 放射性物质在穿过降雨区域时，其强度由于雨水的吸收而减少，减少比率为 常数。 5、 假设地面对放射性气体起全反射作用，地面和海面对放射性气体没有吸附，

将海面视为平原地区

6、假设风向为水平风向，且风向风速不随时间变化。

7、扩散过程中不考虑泄漏点内部温度的变化对气体扩散的影响

四 问题的分析

4.1 问题（ 1）的分析

一座核电站遇自然灾害发生核泄漏， 浓度为 p0 放射性气体以速度 m kg/s 匀速 排除，这近似于放射性物质源是连续均匀稳定的。 在无风情况下， 放射性气体以 速度s m/s匀速在大气中向四周扩散，放射性气体的扩散服从扩散定律，即单位 时间通过单位法向面积的流量与它的浓度梯度成正比。 在这些条件下， 我们明确 了要研究的问题是 点源连续泄漏的扩散问题， 题给条件中明确要求不考虑风力的 影响，但为了使建立的模型更加贴近实际， 需考虑地面反射、 核泄漏源的实际高 度、降雨等因素对浓度分布的影响。由“扩散定律” “放射性物质质量守恒定律” “气体泄漏连续性定理” 可得出无界区域的抛物线型偏微分方程。 再通过假设条 件建立未考虑地面反射、核泄漏源的实际高度、 降雨等浓度影响因素的初步模型， 然后从这些影响因素对模型进行完善， 最终得出 核电站周边不同距离地区、 不同 时段放射性物质浓度的预测模型。

4.2 问题（ 2）的分析

本问是探究风速为km/s时，核电站周边放射性物质浓度的变化情况。当环 境中空气流动产生风力时， 在均匀湍流场中， 扩散参数与下风向距离的关系是明 确的，核泄漏泄漏时间较长时， 可认为扩散是稳定的。 在下风向的湍流扩散相对 于风力引起的移流相可忽略不计， 在流动方向建立 x 轴，不考虑横向速度和垂直 速度。根据假设， 空间中放射性物质的浓度服从高斯分布， 可利用连续点源放射 性物质的高斯扩散模型。 放射性物质在大气中扩散受诸多因素影响， 考虑泄漏源 有效高度、放射性物质自身重力产生的重力沉降、 雨水的吸附等因素对放射性物 质浓度的影响是必要的， 通过这些影响因素对高斯模型进行修正， 然后利用修正 后的高斯模型探究核电站周边放射性物质浓度的变化情况。

4.3 问题（ 3）的分析

本问是要求当风速为km/s时，建立上风和下风L公里处的放射性物质浓度 的预测模型。经分析，此问是问题（ 2）的延伸，我们只需建立合适的坐标系， 将此问题转化为求具体两处的放射性物质的浓度， 便能得出上风和下风 L 公里处 的放射性物质浓度预测模型。

4.4 问题（ 4）的分析 本问要求将之前建立的模型应用于日本福岛核电站的泄泄漏， 计算出福

岛核 电站的泄漏中国东海岸和美国西海岸的影响， 此问实际上是将模型二、 模型三具 体化。通过收集求解模型的相关数据， 利用模型和数据模拟日本福岛核泄露对中 国东海岸和美国西海岸的影响。

五模型的建立与求解

5.1模型一的建立与求解 5.1.1模型一的初步建立

以核泄漏点正下方的地面为坐标原点(0,0,0)，平均风向为X轴、指向下风

方向，铅直方向为Z轴，水平垂直于风向轴(X轴)为丫向，建立空间坐标系, 则核电站泄漏点

O距有效地面的高度为H，则泄漏点位置坐标为0(0,0, H)。

图一：空间坐标系示意图

将气体从泄漏源泄漏时刻记作t=0，时刻t无穷空间中任意一点坐标为(x,y,z) 的浓度记为

C(x,y,z,t),根据假设2,单位时间通过单位法向面积的流量与浓度 梯度成正比，贝U：

q - -「grade

! (i二x,y,z)是扩散系数，grad表示梯度，负号表示由浓度高向浓度低的地

方扩散。

考察空间域0，其体积为V，包围0的曲面为S， S为一规则的球面，S外 法线向量为

(1)

n =(-x,-Y,1)。则在(t,r ■ t)内通过i■■的流量为：

z z 出 --

Q1 二 . q nd；「dt

t

(2)

s

11

内放射性物质的增量为：

Q2 ： ! I i[C(x,y,z,t ：t)-C(x, y, z,t)]dV

V

(3)

从泄漏源泄漏的放射性物质的总量为：

Q = Jt站川 p°dVdt

(4)

根据“质量守恒定律”和“气体泄漏连续性原理”

，单位时间内通过所选曲

面S的向外扩散的放射性物质的量与 S曲面内放射性物质增量之和，等于泄漏源 在单位时间内向外泄漏的放射性物质。贝

Qo = Qi ■ Q

即：

t

(5)

划 —■

s

1

4/t

podVdt

'J

(

/ c\\

iii[C(x,y,z,t 迸)-C(x,y, z,t)]dV t q nd二dt t '

V

6)

又根据曲面积分的Gauss公式:

(芮= UJdivqdV (其中div是散度记号)

s

V

(7)

C(x,y,z,t

：t) —C(x,

y,z,t) At

t ■ t

：t]dV 亠 I iiidivqdVdt 二 t podVdt

•t

丁 .—

t -1

.C,. C(x,y, z, t

lim .:t I—

0

：t)-C(x,y, z,t)

.：t

gradC )dt

lim --------------------------- 「— .：t

t kdiv(

o

由以上两式得：!!!PC]dV t亠iiidivqdV t

V t

；

podV t

--

V

即为：

[可]dV 亠 iiidivqdV 二

：

,1

PodV

V

(8)

V t V

由以上公式并利用积分中值定理得:

■

C

二、idiv(gradC)二2 dx

◎2c

；

:2C

'孑)，t 0,」：：

(9)

x,y,z ，初始条件为作用在坐标 -:t 这是无界区域的抛物线型偏微分方程，根据假设 ：：：1

原点的电源函数，记作

C( x, y,乙 o)oQ

Qo表示泄漏源漏泄释放的放射性物质总量，

方程(9)满足方程(10)的解为:

(< ,y z )

(

io)

(x, y, z)是单位强度的电源函数

Q

z

C(x,y*启2 一厂\"(x,y,z) (11)

此模型只是在不考虑风速的情况下建立的，但为了使模型具有更加的实用 性，下面我们将考虑泄漏源的实际高度、地面反射、降雨等因素对浓度的影响， 完善模型。

5.1.2模型一的最终建立 1.地面反射对模型的完善

泄漏源有一定的高度，且泄漏点源是连续点源，则泄漏点源可视为高架连续

点源，考虑到地面对扩散来的放射性气体有反射作用，根据假设

4，地面对到达

地面的扩散气体完全反射。这儿可认为地面就像镜子一样，对放射性气体起全放 射作用，可用“像源法”处理，如图3,建立三个坐标系，一是以泄漏源（实源） 为坐标原点；二是以泄漏源在地面的投影点为原点，

p点是空间的任意一点，坐

标为（x, y, z）;三是以泄漏源关于地面的像对称源（像源）为原点。把 p点放射 性气体浓度看成两部分（实源与像源）作用之和。

图二：高架连续点源扩散示意图

从以上分析知，p点放射性气体的浓度为实源和像源的放射性气体扩散至此

^z-H ） 点浓度的叠加。则实际泄漏源（实源）对P点的浓度贡献部分可用来表示; 因为地面对扩散物质完全反射，则像对称源（像源）对P点的浓度贡献部分可用

(z H)

e込来表示。于是对（11）式所修正完善的模型为:

(x2

y2)(z_H)2 (z H)2

C(x,y, z,t)=

mpo (4\") t

3/2 x1/2

e」

4

(12)

=(x, y, z)

我们自己模拟一组数据，利用matlab进行仿真模拟，可实现该模型的模拟图 像。假设扩散系数

=0.00001，放射性物质的初始浓度 C =100，扩散时间

t=1000000,放射源总量Q =1000000。（程序见附录一）

核泄漏无风影响扩散

图三：核泄漏无风扩散

5.2模型二的建立与求解 5.2.1模型二的建立

放射性核素云团在大气中迁移和扩散的数值计算基本上可分为二步。

第一步

根据大气动力学理论进行所关心区域中风场的计算，其理论基础是大气运动方 程、连续性方程、状态方程、热力学方程和水汽方程构成的基本方程组。在大气 科学研究领域中，已有多个实用的大气环流模式。第二步进行已知风场中放射性 核素云团迁移和扩散的计算，可采用类似于处理大气污染的方法， 假设放射性核 素云团不影响大气流体速度和温度，求射性核素的连续性方程。

当风速为km/s时，利用连续点源高斯扩散模型分析核电站周边放射性物质 浓度的变化情况。此泄漏点源是有边界点源，泄漏点源的实际高度为 H。以泄 漏点源在地面的投影点为坐标原点，以风向方向为x轴，铅直方向为z轴，与x轴 水平面垂直方向为y轴建立三维坐标系，由于扩散过程中浓度在 y、z轴上的变 化分布符合高斯分布，所以下风向的任意一点

C（x,y,z）的浓度函数为:

C(x, y,z)2

eJ3z2

二 A(x)e」y

根据概率统计我们可以得出方差的表达式为:

Q00 y C(x,y,z)dy

z2C(x, y, z)dz

2

Q0 QQ

二

0

C(x,y,z)dy

f0 C(x, y,z)dz

进而源强的积分公式可以根据假设得出:

Q 二 uC(x, y, z)dydz

把（13）式代入（14）积分可以得出:

1 2 y2

1

将（13）式和（16）式代入（15）式可以得出:

最后再将把（16）、（ 17）式代入（13）式可以得出:

C(x,y,z)二 (2 2 \\]

Q

exp

I y + z

2 u y；- z -

l2G 26丿 一

(13)

(14)

(15)

(16)

(17)(18)

上式为无界空间连续点源扩散的高斯模型， 然而在实际中，由于地面的影响， 烟羽是有界的。根据假设可以把地面看做一镜面， 对泄漏的气体起反射作用，同 样我们可以利用“像源法”进行处理，原理和示意图在模型一的修正中提到，因 此我们得出：

实源的贡献为：

2 2

C1(x, y,z)

像源的贡献为:

-exp( 1 y、 / 1(Z-H)、 z L) 2「)exp(-2

2 2

Q

(19)

C2(x,y,z)二

2 -

则该处的实际浓度为:

1 y2 exp( S77)exp(—2二 y；「z

1 (z H)2

(20)

C(x, y,z^C1(x,y,z) C2(x, y,z)

综合上面的公式得到连续点源高斯烟羽扩散模型：

2 2

2

(21)

2 片 z

522模型二的修正

exp

「^)

{expzH)

['

「

巧]

exp

卜冒^]} (22)

连续点源高斯扩散模型虽然能分析风力对浓度的影响情况， 物质自身的沉降作用和雨水吸附作用对放射性物质浓度的影响， 对连续点源高斯扩散模型进行修正。

但为了能更准确

进而

的探究风力对核电站周边放射性物质浓度的影响情况， 我们将考虑泄漏源有效高 度、放射性

5.2.2.1泄漏源有效高度对模型的修正

如图2所示，泄漏源的有效高度h是由两部分组成，一是核泄漏口距有效地 面的高度H ;

二是在实际核扩散中核泄漏气团从泄漏口排出时，由于受到热力 抬升和本身动力抬升，进而产生的一个附加高度

H。因而h = H —旧。

图四：烟云抬升示意图

附加高度•汨，主要由核泄漏处泄漏气体的气流具有一初始动量（使他们继 续垂直上升）和气流温度高于环境温度产生的静浮力决定， 这两种动力引起的烟 云浮力运动称烟云抬升，附加高度

H即烟云抬升高度，烟云抬升有利于降低地

面的污染物浓度。而且厶H还受到风速、地形地貌等多种因素的影响。

A :当大气稳定度为中性时，计算烟气抬升高度利用 Holla nd公式

Ts u

：H =聖(1.5 2.7Ts 一& D)二］(1.5vD 9.6 10‘Q) u

s

H

(23)

式中：U：泄漏源出口处的平均风速，m/s；

Vs :放射性气体出口流速，已知为 mm/s；

D :泄漏源出口的有效内径;

QH :泄漏源的热排放率,

kw

；

Ts:泄漏源出口处温度,

Ta :环境大气平均温度,

K ,取当地近五年的平均值; Briggs公式计算烟气抬升高度

B:当大气条件为不稳定时，利用

当QH

21000kw 时

x ：：10H x 10H

•旧=0. 3 6Q1/3 x2/3U

)1

H =1. 5 5H1 /3 H 2/®U )

当QH

：：21000kw 时

x 3x

H =0. 3 6Q

*

1 / 3

1

)x 犷

1 / 3-..

x 3x

H =0. 3 3Q3/5 H 2/

x* =0.3曲/5 H3/5 (犷

综上所述, 泄漏源的有效高度为:

(24)

522.2考虑放射性物质自身的沉降作用对模型的修正

放射性物质的沉降速度取决于空气阻力和自身重力， 沉降速度：

「gd2

V

利用斯托克斯公式表示

s

18：

':放射性物质粒子密度，kg/m3 ；

(25)

g :重力加速度，9. 806 5m/s ;

d :放射性物质粒子直径，m ;

:空气的动力粘性系数，可取1.8 10&kg/(ms);

V :沉降速度，m/s，含碘放射性核素的沉降速度为 Vs= 1. 1 cm/ s；

在扩散过程中重力沉降的位移叠加在羽流中心线上，

使中心线向下倾斜，放射性

物质粒子则相当于在下倾的中心线上扩散， 放射性物质的扩散与沉降的叠加可认 为是放射源以vs的速度向下移动。在x处向下移动的垂直距离为Vst^^，即泄

OC

漏源的有效高度h下降了 Vst = 丫泌，泄漏源的有效高度成为h -Vst二h -空，考

a a

虑到地面的全反射作用，反射项的有效高度也变成了 的连续点源高斯扩散模型为：

h-Vst二h-VsX。则修正后

Ct

(z_h 当2

C(x, y,z,h)

(z h±)2

exp(

){exp[ 2

2二

-^^] exp[-^^]}(26)

522.3考虑雨水的吸附作用对模型的修正

降雨对放射性物质的浓度有一定影响，

即雨水对放射性物质有一定的吸附作

用。以吸附系数［来表示雨水对放射性物质吸附作用的大小，［与降雨强度的关 系为：［二

alb，式中I为降雨强度，a,b为经验系数。如果放射性物质含碘，则 a =8 10*,b=0.6;反之，a=1.2 10*,b = 0.5。

雨水的吸附作用导致的放射性物质浓度的减小，可对源强进行修正：

Q(x)=Qexp(

则进一步修正的连续点源高斯扩散模型为:

Px

x

) (27)

C(x, y,z,h)

2 兀 ZyJ

^Q^exp(孚){exp[ - J 2by 2J

尹 exp[一

丄

Vsx 2

(z h

-亠)

-^1—]}

综上所述：修正的连续点源高斯扩散模型为:

(28)

Q()

2

(z — h+密)2 (z + h—^X)2

Q(x) y ' 仃' ' J

C(x, y,z,h) = ------ exp(_化){exp[— ---------- — -------- ]+ exp[— ------ — --------- ]}

2兀吨压 2CTy 2巧 2bz Q(x) =Qexp( V且

s

：x - b

),.：二 al

Ct

s

18:

(29)

523模型二的求解

模型所需参数的选取对模型的求解至关重要，通常情况下气象参数的选取是 利用该地区多年气象资料，采取工业安全与环保统计的方法进行有关参数的确 定，而其他扩散参数则以实际测定为准。

A:大气稳定度的计算

根据我国标准(GB/ T 13201 — 91))制订地方大气污染物排放标准的技术方法 的规定，大气稳定度分为6级，分别为A —极不稳定、B —不稳定、C—弱不稳定、 D—中性、E—弱稳定、F—稳定。该方法的技术路线是：根据核泄漏源所在地的 经度和纬度以及泄漏的日期和时间计算当时的太阳高度角 he，利用天气条件确定

辐射等级，然后利用辐射等级和风速确定大气稳定度， 最后查扩散参数幕函数表， 确定扩散参数。

首先，然后，由太阳高度角ho和云量查出太阳辐射等级；最后，再根据地面 风速确定当时的大气稳定度。

总云量/ 低云量 太阳咼度角h0 夜间 ho 兰15° -2 -1 -1 0 0 -1 0 0 0 0 15= h0 兰35° 35= h0 兰65° h0 沪 65° +1 +1 0 0 0 +2 +2 +1 0 0 +3 +3 +1 +1 0 Y4/ 兰4 5~ 7/ 兰4 >8^4 >7/5~ 7 >8/ >8

表一：太阳辐射等级

地面风 速 太阳辐射等级 (m / s) +3 +2 +1 0 -1 -2

<1.9 A A-B B D E E

2~2.9 3~ 4.9 5~ 5.9 >6 A-B B C C B B-C C-D C C D D D D E D D D F E D D D D D 注：地面风 1速系指距地面10m高度处10min平均风速

表二：大气稳定度的等级

B:扩散参数的确定

扩散参数 n 的确定，采用Briggs给出一套扩散参数幕函数表，如表三和

表四：

大气稳定度 6 A B C D E F

0.22x(1 +0.0001X)% 0.16x(1 +0.0001x)“ 0.11x(^0.0001x^2 0.08x(1 +0.0001X)% 0.06x(1 +0.0001x)12 0.04x(1 +0.0001x)'2 0.2x 0.12x i 0.08x(1 +0.0002x) 2 0.06x(1 +0.0015x)'2 0.03x(1 +0.0003x) 0.016x(1+0.0003x) 表三：Briggs扩散参数（平原地区和城市远郊区）

大气稳定度 ▽x J A-B C D E-F 0.32x(1+0.0004x)」2 0.22x(1+0.0004x)—J 0.16x(1+0.0004x)—2 0.11x(1 +0.004x) 2 y 0.24x(1 +0.0001x) 2 0.20x 0.14x(1 +0.0003x)'2 丄 0.08x(1 +0.0015x) 2

轨 表四：Briggs扩散参数（工业区和城市中心区）

5.3模型三的建立与求解

当风速为km/s时，建立上风和下风L公里处放射性物质浓度的预测模型。

此问题可以利用模型二进行求解，将模型二中的

u替换成k - S,其中k为风速的

向量表示，s为泄漏点放射性气体自身扩散速度的向量表示。 根据假设，风速k沿 x轴正方向，恒为正。对上风向 L公里处放射性气体浓度计算时，s方向沿x轴 负方向，为负；对下风向L公里处放射性气体浓度计算时，s方向沿x轴正方向， 为正。因此此问题转化为求（L,O,z）和

（-L,O,z）两处的放射性气体浓度，即风速为

km/s时，上风和下风L公里处放射性物质浓度的预测模型为：

/ -丄VsX 2

v2

—exp(_〉){exp[——

_ 2jr(k _s)Q(x)

/ 丄—VsX 2

——) 2oz

]+exp[——

(z+h———)

]},上风向 L处

2iz (z + h-^)2

( 30)

- Q()

2

(z」+Vsx)2

2crz

——exp(#){exp[—— ^(k+s)^^^ 2cry

]+exp[-—

2crz

]}，下风向 L处

我们根据模型分别模拟下风向和上风向的放射性物质扩散仿真图，如图 五、图六：（程序分别见附录二、三）

核泄漏模拟图一

60

門门r [Jn -IF -0.4 -0.35

J0

20

--0.3 --0.25 --0.2 0

-20

-40

-60

100

150 2Q0 250 300

图五：核泄漏下风向扩散图

核泄漏模拟图二

so

0.3 07

40

20

0.5

0

0.4

-20

0 3 0.2 0 1

-40

-60 -300

-250

*200 门別 -100

-50

图六：核泄漏上风向扩散图

5.4问题(4)的求解

此问题可结合模型二和模型三求解，我们通过网络收集了相关数据，将这些 数据代入模型进行模拟。

数据显示，日本在核泄漏时期的风向主要为偏西风，风速为

2-3级,辐射源强

度Q= 770000000kg /s。建立如图七(东海岸扩散图)中所示的坐标轴，则中国 东海岸的坐标大概可取为(-1500km,1900km)。建立如图八(西海岸扩散图)中 所示的坐标轴，则美国西海岸坐标大概取为(4300,7500)。

根据模型二的求解方法可确定

1 1

6 =0.16 x(1

0.0004 x)l；二=0.14 x(1

0.0003 x)P，风速为 4m/s， 分别求解福岛

核电站核泄漏泄漏的放射性物质扩散到中国东海岸和美国西海岸的浓度：

中国东海岸：c(1500,1900,0) =4.2429 10”g/m3 美国西海岸：c(4300,7500,0) =2.3854 10‘g/m3

图七：东海岸扩散图

图八：西海岸扩散图

六模型的评价及改进

6.1模型的评价 6.1.1模型优点

第一：模型一，我们考虑了泄漏源高度、地面反射等因素对放射性物质浓度

的影响，利用这些影响因素对浓度预测模型进行了修正完善， 型具有实用性、高精度性、合理性。

建立的浓度预测模

第二：模型二，我们建立了连续点源高斯扩散模型， 并考虑了烟云抬升高度， 求得泄漏源的有效高度，通过泄漏源的有效高度、放射性物质自身的沉降作用、 雨水的吸附作用对模型进行修正， 得到了更优化的连续点源高斯扩散模型。 此模 型能较好的求出核电站周边放射性物质浓度，并与模型一求出的浓度进行比较， 从而给出核电站周边放射性物质浓度的变化情况。

第三：模型三是在模型二的基础上加以修改转化的，因此模型三继承了模型 二的优点，使放射性物质浓度的预测更加精确。

第四：我们建立的以上模型在核事故的应急救援过程中， 对救援人员划定警 戒区和确定周围居民的疏散范围具有重要意义， 并可为制定救援方案和应急决策 提供科学依据。 对放射性云团在空中迁移和扩散提供浓度的定量描述， 也对可能 发生的核事件的放射性核素浓度监测及监测时间范围提供相关信息。

6.1.2 模型缺点 第一：模型只考虑了平原和城市两种情况， 不利于扩散参数的确定。 但地形

复杂多变，不能找到充分的数据对其进行分析。

第二：模型缺乏检验。没有相关数据和标准用于模型的检验。

6.2 模型的改进

日本是多山国家， 考虑山林对放射性物质扩散的影响能更好的预测和估算放 射性物质的浓度； 放射性物质的扩散， 特别是向其他国家地区扩散， 放射性物质 的扩散大多要跨洋越海， 考虑海水对放射性物质的吸收显得至关重要， 将此因素 考虑进去，能使模型更加贴近实际。

七 参考文献

[1] 卓金武，魏永生，秦建，李必文． Matlab 在数学建模中的应用 [M] ．北京： 北京航空航

天大学出版社， 2011．

[2] 程勇，于林，姚安林．采用高斯模型分析输气管道泄漏后气体的扩散 [J] ．内 蒙古石油化

工， 2010 ，14： 49-51．

[3] 张斌才，赵军．大气污染扩散的高斯烟羽模型及其 GIS 集成研究 [J] ．环境 监测管理与技

术， 2008， 20（5）：17-19．

八 附录

附录一：

%绘制当没有风影响时放射性物质的扩散图形 function []=expand(k,Q,C,t,r,ws,wa,wd) R=sqrt((-4)*k*t*log((C/((1-r)*Q)).*(4*pi*k*tF(1.5))) for i=0:24:t/3600 if (wd*(i/24))<=R

R=R-wd; R0=R; t=0:0.1:0.5*pi; p=0:0.1:2*pi; wsx=ws*cos(wa)*i; wsy=ws*sin(wa)*i;

[theta,phi]=meshgrid(t,p); x=R*sin(theta).*cos(phi)+wsx; y=R*sin(theta).*sin(phi)+wsy; z=R*cos(theta); hold on surfc(x,y,z) end end R0=R0 %绘制其它相关信息

view(20,70) %观察角度 title(' 核泄漏无风影响扩散 ');%图形标题 xlabel('x'); %x 轴标记

ylabel('y'); %y 轴标记 zlabel('z'); %z 轴标记 axis([-10,15,-10,20,0,12]) %xyz 轴显示范围余量 x0=ws*cos(wa)*1000000/3600 y0=ws*sin(wa)*1000000/3600 R仁 sqrt((x0)2+(y0)2) End

附录二：

Q=2; %输入源强

AA

k=2; %输入风速 s=0.5; u=k+s; d=1; %步长 Z=0.4;

x=10:d:300; %下风向距离

y=-70:d:70; %横风向距离 [x,y]=meshgrid(x,y); by0=0.08*x.*(1+0.0001*x).A(-1/2);

bz0=0.06*x.*(1+0.0015*x).A(-1/2); tempy1=-y.*y./by0./by0./2;

tempy2=2.718282.A(tempy1); c=Q/pi/u*((by0.*bz0).A(-1)).*tempy2 Cs=2000; %输入求解的条数 contour(x,y,c,Cs); shading interp; colorbar; grid;

title(' 核泄漏模拟图一 ') 附录三： Q=2; %输入源强 k=2; %输入风速 s=0.5; u=k-s; d=1; %步长 Z=0.4;

x=-300:d:-10; %下风向距离 y=-70:d:70; %横风向距离 [x,y]=meshgrid(x,y); by0=0.08*x.*(1+0.0001*x)A(-1/2); bz0=0.06*x.*(1+0.0015*x)A(-1/2);

tempy1=-y.*y./by0./by0./2; tempy2=2.7182824(tempy1); c=Q/pi/u*((by0.*bz0)4(-1)).*tempy2

Cs=2000; %输入求解的条数 contour(x,y,c,Cs); shading interp; colorbar; grid; title(' 核泄漏模拟图二 ')