基于MATLAB的放射性气体扩散模拟及应用
Matlab在危险气体扩散模拟分析中的应用
图 1 气体扩散影响区域划分
一般情况下 ,气体扩散后影响区域可划分为 :致死区域 、 重伤区 、致伤区 、吸入反应区和安全区域 。由于致死区和重 伤区都极易导致人口的伤亡 ,因此在分析中将这 2 个分区统 一合并为致死区[7] 。
[ x ,y ] = meshgrid (50 : d :1000 , - 100 : d :100) ; %设置解的 范围
c = (Q. / 3. 14. / u. / by. / bz) . 3 tempy2. 3 tempz2 ; %计算 浓度
·26 ·
cs = input (’cs = ’) ; %设置等高线浓度值 contour (x ,y ,c ,cs) ; a = contour (x ,y ,c ,cs) ; %作等高线 clabel (a) ; shading interp ; colorbar ; grid ; xlabel (’下风向距离’) ; ylabel (’横风向距离’) ; title (’危险气体扩散浓度分布图’) ; 运行程序 ,得出的结果见图 3 、图 4 。
当 H = 0 时 ,得到地面连续点源扩散模式为 :
c ( x , y , z ,0)
=πuσQσy z exp (
-
y2 2σ2x
)
·exp
(
-
z2 2σ2y
)
(2)
由式 (1) 和式 (2) 可以看出 , 气体扩散浓度和源强 、风速
及 x , y , z 轴上的扩散参数都有着一定关系 。在后文研究中
速 ,可以按照气象部门的统计资料来确定 ,也可以按照实际 如表 1 。
气体扩散模拟实验报告
气体扩散模拟实验报告工程力学1001 陈金刚3100104568一．实验背景：要求用matlab编程模拟分子碰撞，演示气体扩散情况。
本实验中的模型采用简化形式，所发生碰撞均为完全弹性碰撞。
壁面压强产生的原因是大量气体分子对容器壁的持续的、无规则撞击产生的。
各个壁面所受力和压强为单位时间内的平均值来代替。
由Ft=mv-（mv）=2mv，四个壁面在经历相同时间t的情况下，所受的平均压力F与2mv成正比。
又四个壁面长度S一样，所以由P=F/S 可知，壁面所受压强P与所受平均压力成正比。
因此，只需统计分析在相同时间内，四个壁面所受碰撞的总的2mv 即可知道各个壁面所受压力和压强的情况。
二．实验基本情况说明：区域尺寸：2维空间，200*200粒子数目：N个，N可变化粒子半径：各不相同，根据实验情况在一定范围内随机产生粒子质量：各不相同，与粒子半径的平方成正比初始位置：随机分布在左半区域初始速度：随机大小三．实验结果分析：（所有运行时间都一致）由表格数据分析可知：对于不同粒子数目，在运行足够长时间下，区域内分子运动会近似稳定，各个壁面所受压力，趋向稳定，在误差允许范围内，各个壁面所受压强可认为近似相同。
四．源代码：N=50;%球数，可变rad=rand(N,1)*2+4;%球半径，[4，6]pos=[rand(N,1)*90+5 rand(N,1)*190+5];%初始位置：左半边区域vel=rand(N,2)*20-10;%各球初始速度color=rand(N,3);%各球颜色，随机产生mass=5*rad.^2;%各球质量，与半径的平方成正比Left=0; %左边界，下同Right=200;Up=200;Down=0;figure;axis manual;axis equal square; %固定坐标axis([0 200 0 200]);hold on;%===============画不重叠的N个小球====for i=2:Nindex=0;while index==0index2=0;for j=1:i-1pp(j,:)=[pos(i,1)-pos(j,1),pos(i,2)-pos(j,2)];rr(j)=rad(i)+rad(j);ppp(j)=norm(pp(j,:));if ppp(j) <rr(j)+sqrt((vel(i,1)-vel(j,1))^2+(vel(i,2)-vel(j,2))^2)*0.1;%修正一开始若球相切时，后面判断可能误认为相撞index2=1;break;endendif ~index2index=1;else pos(i,:)=[rand()*90+5 rand()*190+5];endendend%==================================================================== couleft=0;%左壁面撞击的总mv，以下类似couright=0;couup=0;coudown=0;dt=0;%最小时间qq=0;%用于循环次数控制while qq<200%=========以下各球两两碰撞最小时间计算k=1;for i=1:Nfor j=i:Ndis=[pos(j,1)-pos(i,1),pos(j,2)-pos(i,2)];vv=[vel(j,1)-vel(i,1),vel(j,2)-vel(i,2)];dist=norm(dis);rr=rad(i)+rad(j);cosAlpha =abs( sqrt(1-(rr/dist)^2));cosTheta=(dot(dis,vv)/norm(vv)/dist);if cosTheta>=cosAlpha && cosTheta<1dd=dist*cosTheta-sqrt(rr^2-(dist*sqrt(1-cosTheta^2))^2);time(k)=dd/norm(vv);k=k+1;endendendtball=min(time);%================各球碰墙最小时间计算for i=1:Nif vel(i,1)>0tx(i)=(Right-pos(i,1))/vel(i,1);elsetx(i)=(Left-pos(i,1))/vel(i,1);endif vel(i,2) >0ty(i)=(Up-pos(i,2))/vel(i,2);elsety(i)=(Down-pos(i,2))/vel(i,2);endendtwall=min(tx,ty);%====tBallWall=min(tball,twall);dt=dt+tBallWall;%与墙相撞改变速度for i=1:Nif pos(i,1)-rad(i)<=Left && vel(i,1)~=0pos(i,1)=rad(i);%修正边couleft=couleft+mass(i)*abs(vel(i,1));vel(i,1)=-vel(i,1);elseif pos(i,1)+rad(i)>=Right && vel(i,1)~=0 pos(i,1)=200-rad(i);%修正边couright=couright+mass(i)*abs(vel(i,1));vel(i,1)=-vel(i,1);endif pos(i,2)-rad(i)<=Down && vel(i,2)~=0pos(i,2)=rad(i);coudown=coudown+mass(i)*abs(vel(i,2));vel(i,2)=-vel(i,2);elseif pos(i,2)+rad(i)>=Up && vel(i,2)~=0pos(i,2)=200-rad(i);couup=couup+mass(i)*abs(vel(i,2));vel(i,2)=-vel(i,2);endend%两球碰撞改变速度for i=1:Nfor j=i+1:NtwoBall=[pos(i,1)-pos(j,1),pos(i,2)-pos(j,2)];D=norm(twoBall);if D-rad(i)-rad(j)<sqrt((vel(i,1)-vel(j,1))^2+(vel(i,2)-vel(j,2))^2)*0.1;if D<rad(i)+rad(j) %修正碰撞误差pos(j,1)=pos(j,1)-vel(j,1)*0.05;pos(j,2)=pos(j,2)-vel(j,2)*0.05;endc=((mass(i)-mass(j))*vel(i,1)+2*mass(j)*vel(j,1))/(mass(i)+mass(j));vel(j,1)=((mass(j)-mass(i))*vel(j,1)+2*mass(i)*vel(i,1))/(mass(i)+mass(j)); vel(i,1)=c;c=((mass(i)-mass(j))*vel(i,2)+2*mass(j)*vel(j,2))/(mass(i)+mass(j));vel(j,2)=((mass(j)-mass(i))*vel(j,2)+2*mass(i)*vel(i,2))/(mass(i)+mass(j)); vel(i,2)=c;endendend%===============每隔一段时间画图while dt>0.05cla;for i=1:Nalpha=0:pi/20:2*pi;xx=pos(i,1)+rad(i)*cos(alpha);yy=pos(i,2)+rad(i)*sin(alpha);plot(xx,yy,'k-','LineWidth',2);fill(xx,yy,color(i,:));enddt=0;drawnowendpos=pos+0.1*vel;qq=qq+1;end%输出各个统计量clc;disp(['LeftSide=',num2str(couleft),' ','RightSide=',num2str(couright),' ',...'UpSide=',num2str(couup),' ','DownSide=',num2str(coudown)]);。
MATLAB在模拟煤层气扩散中的应用
ｇ 单位 时 间交换 质 量 ，ｇ （ ・ ） 一 ｋ／ ｍ ｓ ；
Ｄ 扩 散 系数 ， ｓ 一 ｍ／。
２ 模 拟 实列
２ １ 地 质概 况 ．
宁武 煤 田地 处 山西 地 区北 中部 ， 田南北 长 ， 煤 东 西窄 ， 呈北 东 向斜 列 ， 布 于鄂 尔 多 斯 地 台 ， 展 吕梁 地
器 Ｐ Ｅ Ｏ ＬＧ Ｉｐｄｏ１ ， 以形 象直观 的模 拟煤层 中的煤层 气在储 层 中扩散 与分 布 。现 以宁武 煤 Ｄ Ｔ Ｏ Ｕ （ ｅｔｏ） 可
田太原组 ４号为例 ， 模拟结果表明， 煤层气浓度扩散与实际具有较好的吻合特征 。
关键 词 质 量 守恒 ； 克定律 ； 菲 数值 模拟 ； Ｔ Ａ ＭＡ Ｌ Ｂ
块 、 台地 块及 内蒙 地轴 之 间为一 继 承性上 叠 构造 盆 五 地 ， 总体 看来 ， 造 复 杂 程 度 尚属 简单 类 型 。本 区 从 构
含 煤地 层 为二叠 系 下 统 山 西组 、 炭 系上 统 太 原 组 、 石
中统本 溪组 ， 要可 采 煤 层 为 太 原组 ４， ９号 煤 层 ， 主 ６，
式中：
马
胜等： ＴＡ ＭＡ Ｌ Ｂ在 模 拟煤层 气扩散 中的 应 用 动态模 拟 的结果 ， 以在 对话 框 中选 择 ａｉａｉ （ 可 ｎｍ ｔｎ 动 ｏ 态 模拟 ） 最后再 点 击 Ｐｏ选 项 （ ２ 。 ， ｌ ｔ 图 ）
Ｃ 浓度 ，ｇ ｍ ； 一 ｋ ／
￡ 时间 ，； 一 Ｓ
认 为煤 层气 在煤 粒 孔 隙 扩散 遵 从 菲 克 （ ｉ ） 二 定 Ｆｃ 第 ｋ 律 既 非稳 态 扩 散 … ， 稳 态 扩 散 模 式 中 ， 要 有 两 非 主 点 ： 是认 为煤 基质块 内 甲烷浓 度从 中心 到边缘 是 变 一
气体扩散过程的数值模拟与分析
气体扩散过程的数值模拟与分析随着科技的发展和工业的进步，气体扩散过程的数值模拟与分析变得越来越重要。
气体扩散涉及到很多领域，如化工、环保、安全等，因此对于气体扩散过程的研究和掌握，对于我们的生活和工作具有重要意义。
一、气体扩散过程的定义气体扩散是指气体分子在压力差或浓度梯度的作用下，从高浓度区域向低浓度区域的传输过程。
气体扩散是一种无需介质（如流体）的传输方式，因此它与传统的流体传输方式有所不同。
而气体扩散过程的数值模拟与分析可以推算出气体扩散规律，更好地理解气体扩散的机制。
二、气体扩散过程的机制气体扩散的机制主要有两种，一种是分子扩散，另一种是对流扩散。
1. 分子扩散分子扩散是指气体分子在浓度梯度作用下，自由移动并和周围分子相互碰撞传输的过程。
气体分子在高浓度区域与周围分子相互碰撞的次数较多，因此相对速度较小；而在低浓度区域，气体分子之间的碰撞次数较少，相对速度也较大。
分子扩散是一种相当自由度高的传输方式，但也存在一定的局限性，如直接关系到气体分子的性质和大小等。
2. 对流扩散对流扩散是指气体在流体中的扩散过程，即气体随着液体或气体的流动而传播。
在对流扩散中，气体的扩散速率与气体分子的速度、流体的流速和其他因素都有关系。
对流扩散通常发生在工业、环保、药品等领域。
三、数值模拟与分析的重要性数值模拟与分析是研究气体扩散过程的重要方法之一。
通过数值模拟和分析，可以有效地模拟和分析气体扩散的过程和规律，预测气体扩散的变化趋势，提供相关的参考和指导意见。
此外，数值模拟还可以减少实验的成本和时间，并能够更好地控制操作和管理。
四、常用的数值模拟方法目前，常用的数值模拟方法包括有限元法、有限差分法、有限体积法等。
这些数值模拟方法都能够通过计算机程序模拟出气体扩散的过程和规律，为研究和开发气体扩散相关的产品和设备提供可靠的依据。
总的来说，气体扩散过程的数值模拟和分析是一个非常重要的研究领域，通过模拟和分析可以更好地掌握气体扩散规律，为相关领域的研究和开发工作提供有力的支撑。
放射性的气体预测模型
放射性气体扩散的预测摘要 2011年日本近海发生地震并引发海啸，沿海的核电站受到破坏，释放了大量具有放射性的物质，福岛第一核电站的核泄漏引起了来了国际社会的广泛关注。
在本文中，我们将针对核泄漏释放出的放射性物质的浓度的预测建立模型。
针对问题一：我们考虑核电站周围不同距离地区、不同时段放射性物质上网浓度建立模型。
设气体扩散在空气中的密度梯度为（drdp 0）,质量流J （单位时间内流经管道任一截面流体的质量）质量流的表达公式为J=-31V ∙λ∙(dzdp )∙S ∆ (dz dp )∙S ∆(其中Z=Z 0)。
根据菲克扩散定律：J=t M ∆∆=-D ∙(dz dp )∙S ∆(其中Z=Z 0) 通过积分最后可得物质浓度为:P(r →r+∆r)=]2^2^[32]2)^([*2]2^2^[320000）（）（Vr t V P kK V r r r r M V r t V P K -+∆+∆-πππ(r ∆取1米左右)（M 表示摩尔质量;V 表示摩尔体积） 针对问题二：由于要考虑风速的影响，我们仍可用问题一建立的模型来求解。
在上风部分和下风部分的气体浓度不同，在上风处风速的影响下扩散速度为（V-K ），在下风处风速的影响下扩散速度为（V+K ），而此时假设气体向外扩散的整体分布为一个半椭球形，椭球的短半轴为Vt ，长半轴为（V-K ）∙t+Kt=Vt.故可以把此问题的扩散范围理想化为一个半球体，上风处大于下风处的浓度而整个半球在以风速方向过电站的直线为对称轴的浓度左右具有对称性。
针对问题三：目前比较常见的有害气体泄漏与扩散机理有高斯云羽扩散、高斯云团扩散、重气云扩散和非重气云扩散、FEM3模型等，在本题中我们用高斯扩散模型进行模型建立。
可得到高架连续点源泄漏的浓度分布为：C(X 、Y 、Z 、H)=]}2)^(22)^(exp[]2)^(22)^(exp[{]2)^(22^[exp 2x H Z y H Z x y z y u Q σσσσσπ+-+--∙- 在应用高斯模型时最关键的是确定扩散参数，再应用程序计算出各处的质量浓度。
LNG低温储罐泄漏事故危险性分析综述(全文)
第三届中国LNG 论坛论文编号：1260304LNG 低温储罐泄漏扩散危险性分析方法综述张文冬张永信路江华（中国石化青岛液化天然气有限责任公司，山东青岛，266400）摘要:液化天然气（LNG ）低温储罐泄漏可能导致喷射火、闪火、沸腾液体扩展蒸气爆炸和蒸气云爆炸等事故，国内外学者在LNG 泄漏扩散的危险性分析方面取得了一定的理论和实验研究成果。
文章从数值模拟和实地物理模拟两个方面，总结了近年来国内外学者LNG 泄漏事故危险性分析尤其是泄漏扩散方面的成果。
数值模拟方面，介绍了基于不同原理泄漏扩散的数值计算模型，如三维传递模型、高斯模型、经验模型、BM 模型、一维积分模型、浅层模型等的建立及其在LNG 泄漏扩散中的应用，并总结了目前广泛应用于定量风险评估的计算机软件。
物理模拟试验以风洞试验为主，本文归纳了国内外LNG 泄漏试验的研究成果，作为数值模型的验证和进一步补充。
通过对比分析不同的分析研究方法在LNG 储罐泄漏事故中的应用，为进一步研究提供参考，以便更加合理地进行LNG 接收站安全设计。
关键词：LNG 储罐；泄漏；扩散；数值模拟；试验随着社会经济的发展，液化天然气（LNG ）作为一种清洁能源越来越受到人们的普遍关注，其在储存过程中的安全性也备受关注。
液化天然气储罐一旦发生泄漏，将可能导致喷射火、闪火、沸腾液体扩展蒸气爆炸和蒸气云爆炸等事故，会对周围的环境、人员、设备等造成极大危害。
国内外此类事故常有发生，如美国俄亥俄州发生的天然气泄漏事故造成128人死亡和难以计数的财产损失。
许多国家对LNG 的相关设施制定了严格的安全标准，如美国联邦规范将LNG 设备的安全归于管网安全。
国内外学者对于LNG 泄漏扩散危险性分析的理论和试验研究工作已开展多年，并且得到了很多具有实际指导意义的成果。
1.LNG泄漏模型LNG 储罐的泄漏形式包括气相泄漏、液相泄漏和两相流泄漏，不同的泄漏模式对应有不同的泄漏模型。
由于LNG 储罐主要发生液相泄漏和两相流泄漏，这里只介绍这两种泄漏模式的预测模型。
matlab 爆炸力学
matlab 爆炸力学一、什么是爆炸力学爆炸力学是研究爆炸现象的科学，包括爆炸物质的性质、爆炸过程中物质的运动和变化规律、爆炸波的传播规律以及对周围环境的影响等方面。
二、matlab在爆炸力学中的应用matlab是一种高级技术计算语言和交互式环境，具有强大的数值计算和数据可视化功能。
在爆炸力学中，matlab可以用于模拟和分析各种爆炸现象，如冲击波传播、气体扩散等。
1. 冲击波传播模拟冲击波是一种由高压气体形成的压缩波，它是爆炸产生的主要效应之一。
利用matlab可以建立数值模型来模拟冲击波在不同介质中传播时的变化规律。
例如，可以利用有限元法建立三维模型来模拟地下核试验产生的地震波。
2. 爆轰反应分析爆轰反应是指在极端条件下（如高压、高温等）下发生的快速氧化反应。
利用matlab可以建立化学反应动力学模型来分析爆轰反应的机理和过程。
例如，可以利用matlab分析TNT等炸药的燃烧过程。
3. 气体扩散模拟在爆炸过程中，气体会迅速扩散并对周围环境产生影响。
利用matlab 可以建立数值模型来模拟气体扩散的过程和规律。
例如，可以利用Navier-Stokes方程建立流体力学模型来分析气体扩散的特性。
4. 爆炸后果评估爆炸事件会对周围环境和人类产生严重影响，因此需要进行后果评估以制定相应的安全措施。
利用matlab可以建立数值模型来预测爆炸事件的后果，如伤亡人数、房屋损坏等。
三、matlab在爆炸力学中的优势1. 强大的数值计算能力matlab具有强大的数值计算能力，能够快速处理大量数据，并进行高精度计算。
2. 丰富的可视化功能matlab具有丰富的可视化功能，可以将计算结果以图表或动画的形式展示出来，方便用户进行分析和理解。
3. 灵活的编程语言matlab具有灵活的编程语言，可以根据用户需求进行定制化开发，并与其他软件进行集成。
4. 大量的工具箱支持matlab拥有大量的工具箱支持，如信号处理、图像处理、优化等，可以满足不同领域的需求。
基于MATLAB的放射性气体扩散模拟及应用
基于MATLAB的放射性气体扩散模拟及应用摘要：随着核电在世界范围内广泛应用，越来越多的核电站建立起来，但是，运用核能的同时我们也要防范于未然。
因此，模拟核事故放射性气体扩散有着非常重要的意义。
本文运用概率动力学的相关知识，结合高斯扩散模型，建立放射性气体扩散模型，模拟了福岛核电站放射性气体扩散对我国东海岸的影响。
关键词：核事故放射性气体扩散；高斯扩散模型；matlab放射性气体扩散模型1.前言北京时间2011年3月11日，日本福岛县的福岛第一核电站发生了一起重大核事故，大量的放射性污染气体从事故的核电站泄露进入大气，对大气环境产生了非常严重的污染，短时间内事故等级从四级跃升到最高等级——七级核事故，引起了国际社会的广泛关注。
3月15日，专家组分析相关数据得出较低浓度的放射性气体正从核电站向福岛以东地区扩散，并可能在将来几天内到达北美地区，最终到达欧洲地区。
同时事故核电站10km范围内的所有居民被日本政府要求紧急撤离。
核电站周围的各个监测站检测到碘，氩，钚等多种放射性同位素从核电站泄出，23日，在核电站厂区内检测出中子辐射。
随着核电在世界范围内广泛应用，越来越多的核电站建立起来，但是，运用核能的同时我们也要防范于未然。
因此，模拟核事故放射性气体扩散有着非常重要的意义。
本文运用概率动力学的相关知识，结合高斯扩散模型，建立放射性气体扩散模型，模拟了福岛核电站放射性气体扩散对我国东海岸的影响。
首先，我们考虑到风向和风速对放射性气体浓度分布有一定影响，由于当时日本发生了大地震和大海啸，这使得当地的气象环境十分复杂。
本文我们结合高斯烟羽模型，并考虑烟气抬升，地面反射，干湿沉积，放射性气体的衰变等多种因素对模型进行反复修正，得到最终的模型。
之后，对于上风向和下风向L公里处的放射性气体浓度，只需在上述的基础上，令x=L或-L，同时，将风速k用（k+s）和（k-s）来代替，y=0，z=0。
建立完这个基础模型后，我们就可以以此为基础研究在风速一定的情况下，位处上下风L公里处，放射性气体浓度的估计模型，并用matlab软件运行模型，模拟出在下风向时的浓度分布图。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
基于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模型的建立该模型运用了概率动力学的知识，通过图解和数学的推导得出连续点源周边的放射性化学物质高斯扩散的模型。
下面进行模型公式推导由正态喷射气体中的浓度能量分布的公式假设正态气体中的浓度分布可以被通过推导而得出下式为风向任意一点处处当泄漏正态气体时其浓度的乘积函数为其形式为：上式为我们所求的高斯模型的公式，其应用条件为无界限的空间中，但是在我们生活中，由于没有存在地面的泄露气体存在，导致地球上的烟羽无法得到无限的释放和扩散。
所以若根据高斯模型假设，我们可以把存在地面的烟气看做一个反射体，该反射体对无界空间泄漏的气体起全反射的作用，便等于我们可以直接采用像源法对其进行扩散处理。
对于任何一个无界空间p点的浓度都是可以用像源法看做两个组成部分的影响与贡献之和：其中一部分指的是不泄露气体存在于地面时的无限扩散排放所导致的泄漏物扩散排放浓度；另一部分指的是由于无界空间地球表面对无界空间泄露气体的无限全反射作用导致增加的无界空间泄漏物扩散排放浓度。
即我们可以用像源法表示为该处的无界空间泄漏物扩散排放浓度相当于不泄露气体存在地面时位于（0，0，h）的实源和同时位于（0，0，-h）的像源在无界空间p点处所扩散排放造成的增加泄漏物浓度之和。
图2 地面反射示意图其中，实源的贡献为：以上式列表中的一个泄漏源有效值的高度一般也都是用来指一个泄漏源在烟云中由于气体流动形成的高度在泄漏气云基本上已经被转变成一个类似水平状的物体时候。
对于直接影响出口烟云的抬升和出口高度的影响因素主要有很多，这其中主要因素包括：烟云泄漏气体的初始流动速度和其方向、初始的温度、泄漏口的直径、环境中的风速及出口风速岁高度的变化率、环境中的温度及大气的稳定度。
查资料得：a. 由于有风（k>1.5m/s），中性和不稳定的条件，按照以下式计算烟气抬升的高度（1）当两个烟气热能量释放相应温度与烟气功率相应qh的差分数值平均大于或小过等于35k时也就是2100kj/s，且（2）当两个烟气热能量释放相应温度与功率相应的烟气环境热能量释放相应温度的差分数值△t平均大于或小于等于35k时，△h的差分数值可以采用上列下式方法进行精确计算：式中：no----烟气热状况及地表系数，见表1；n1----烟气热释放率指数，见表1；n2----排气筒高度指数，见表1；Qh----烟气热释放率，KJ/s；H----泵与排气筒间的距离离地面的平均几何高度，m，超过高度系数减去240m时，取h=240m； Qv----实际排烟率，m3/s；△t----烟气发生器出口的温度与环境的温度差，k；Ts----烟气出口温度，K；ta----k为环境中在大气中的平均温度，k，如无任何实测值，可能会取决于其所邻近的国家气象台（包括观测气象站）季或月的年平均值；u-排气筒出口处平均输出风速，m/s.表2 n0，n1，n2选取式中：D----排气筒出口直径，m；△H2----按a方法计算，no、n1、n2按表5中Qh值较小的一类选取；Qh、U----与a中的定义相同。
式中，为一个介于烟囱所在几何图形高度以上的平均温度大气压的温度梯度，k/m。
这里我们还需要特别注意的是保证>-0.0098。
4.1.2 模型的修正对于一个标准的高斯模型，我们计算得出的是一个理想的模型。
但是在实际的计算中核素云团扩散的过程还是存在着其他放射性因素的直接影响。
这些因素比如：大气中粒子的重力运动沉降、雨水的冲洗沉积以及大气中核素云团衰变等的放射性因素对高斯浓度的分布等。
因此，我们必须一定要对理想的高斯浓度核素云团扩散的模式进行相应的修正，才使计算能够更准确、有效的反映出实际的扩散规律。
4.1.3 模型求解Matlab模拟下风向的浓度分布图如下：图3 软件模拟图4.2将模型应用于福岛核电站的泄漏对我国东海岸影响结合海域平均风速对人工利用放射性物质气体浓度的分布产生影响的高斯模型与上风向和下风向l公里以内处的放射性物质气体平均浓度的模型，再通过参阅和整理大量的气象资料、地理条件、新闻材料等，确定了大气稳定度标准等级，扩散系数以及相应海域平均风速、地理位置和距离等各种相关的参数，对于我国东北部海岸的人工利用放射性物质气体的浓度分布进行了预测，与实际测量的平均值分布进行了比较。
以我国东北部海岸——山东半岛海域为中心的实例，综合分析考虑海域平均风速风向，地理位置和距离等。