污染空气的扩散模型

1、放射性气体扩散的预估模型摘要 ：由于放射性气体泄漏造成惨重损失的报道在国际屡见不鲜，近日日本福岛核电站的放射 性气体的泄漏事件更让我们关注放射性气体泄漏时在环境中的浓度问题，为了今后事故发生后提供 积极的补救措施 , 所以对放射性气体的扩散作深入的研究是很有必要的。本文结合高斯烟羽模型、 线性拟合，以及微分方程模型，运用 MA TLAB 软件，分析了泄漏源强度、风速、大气稳定度参数、 地面粗糙度参数和计算精确度等的因素对放射性气体扩散的影响，预测了放射性气体浓度在不同时 间，不同地区的浓度变化，并且本文模型中的数据可以根据不同的实际情况而加以改变，因而使本 文的应用范围大大增加， 可以适用于具

2、有较强的应用性。 文章首先在第一问中利用 MATLAB 软件对 数据进行线性拟合，采用微分方程模型得到核电站周边放射性气体在不同地区，不同时间段的浓度 变化，得出随着离泄漏源距离的延伸，最终放射性物质的浓度越来越小，趋近于零，即当 L 趋向无 穷是，C(x,y,z,t)趋向于零；当时间趋于无穷时 ，C (x,y,z,t)也趋于无穷。问题二，问题三中，建立以 核电站周边不同地区得距离以及风速为因变量，设置各个主要因素的参考数据，同时，利用高斯烟 羽模型对核电站周边地区的浓度进行预测，然后，利用 MATLAB 软件，将相关数据代入程序，我们 得到核电站周边地区的浓度分布的等高曲线。问题四中，通过实

3、际收集数据，集合核电站周边地区 的浓度等高曲线，可以直观的看出日本福岛核电站对我国东海岸以及美国西海岸的影响。一 问题的提出1.1 背景的介绍目前，核电的发展给国家带来了巨大的经济效益和社会效益， 但核电正常运行以 及发生泄露时不可避免的会有气载放射性核素排出， 这样就给周围的环境产生了一 定的影响，因此，正确的测出大气中放射性物质的浓度在环境检测以及安全评估中 具有重要意义 。1.2 需要解决的问题 设有一座核电站遇自然灾害发生泄漏，浓度为 p0 的放射性气体以匀速排出， 速度为 m kg/s ，在无风的情况下，匀速在大气中向四周扩散 , 速度为 s m/s.(1) 请你建立一个描述核电站周

4、边不同距离地区、不同时段放射性物质浓度的 预测模型。(2) 当风速为 k m/s 时，给出核电站周边放射性物质浓度的变化情况。当风速为 k m/s 时，分别给出上风和下风 L 公里处，放射性物质浓度的预 测模型。（4）将你建立的模型应用于福岛核电站的泄漏，计算出福岛核电站的泄漏对我国东海岸，及美国西海岸的影响。计算所用数据可以在网上搜索或根据具体情况自己模拟基本假设1. 气体的扩散看作空中某一连续点源向四周等强度地瞬时释放气体，放射性气体在无穷空间扩散的过程中不发生性质变化2. 气体的传播服从扩散定律，即单位时间通过单位法向面积的流量与它的浓度 梯度成正比3. 定常态，即所有的变量不随时间变化

5、4. 假设释放的气体的密度与空气相差不多（不考虑重力或浮力的作用），且气体 扩散过程中没有发生化学反应5. 扩散气体的性质与空气相同6. 扩散气体达到地面时，完全反射，没有任何吸收7. 假定地面水平8在下风向上的湍流扩散相对于稳流相可忽略不计9.风向与地面水平，且在气体扩散的过程中保持不变符号说明与名词解释t 气体扩散时间，气体由泄露源泄漏时刻t=0x,y,z以泄漏源为坐标原点，空间任意一点的坐标C 空间中任一点的气体浓度k 气体扩散系数Q气体由扩散源扩散时施放的气体总量J -平均风速二y -用浓度标准偏差表示y轴上的扩散参数-Z-用浓度标准偏差表示的z轴上的扩散参数 H气体扩散的有效高度x下

6、风方向到泄漏点源的距离y侧风方向离泄漏源点的距离z垂直向上方向离泄漏源点的距离I-距离泄漏源的距离5 泄漏源的总浓度m-放射性气体排出的速度k 风速s-放射性气体排除后向四周扩散的速度四问题分析4.1问题（1）核电站源源不断泄漏引起的气体扩散传播可以看作在无穷空间由连续点源导 致的扩散过程，能够由二阶抛物型偏微分方程描述放射性气体扩散过程中浓度变 化的规律。本问中由于不考虑风力的影响，且扩散出来的放射性气体匀速向四周 散开，这样经过任意时刻t，扩散的气体围成一个半径为st的球体，且距离球心 位置不同的地方浓度值不同。4.2问题（2）当环境中空气流动时，在均匀湍流场中，扩散参数与下风向距离的关系

7、是明确的， 所以泄放时间较长时，可以认为扩散是定常的。在下风向上的湍流扩散相对于移 流相可忽略不计时，在流动方向上建立 x轴，横向速度为V,在不考虑垂直速度， 并且假设空间中放射性气体云的浓度服从高斯分布的情况下，运用高斯模型可以 较合理的计算出核电站周边地区的放射性气体浓度。4.3问题（3）在考虑风速的情况下，我们之前已经假设风向与水平 x轴正方向一致，由于气体是向 四周扩散，这样在下风处，气体相对与地面的扩散的最大速度为（ k+s）m/s,在上风处, 气体相对与地面扩散的最大速度为（k-s）m/s（令k>s）,然后我们分别代入到2中建立的 高斯烟羽改进模型中，分别用（k+s）和（k-

8、s ）去代替方程中的u，同时分别令x=l,y=z=O 或者x=-l，y=z=0,这样就可以求出当风速一定时，上风和下风l公里处，放射性物质浓 度的估计模型。4.4问题（4）在以上的分析中，我们可知，通过收集福岛核电站和我国东海岸以及美国西海岸之间的 距离，以及和两地之间的风向，风速数据，可以大致判断出福岛核电站的泄漏事故对我 国东海岸以及美国西海岸的影响，然后进行仿真模拟得出比较准确的结果。五.模型的建立5.1问题（1）将气体从泄漏源泄漏时刻记作 t=0,泄漏点选为坐标原点，时刻t无穷空间中任一点 （x,y,z ） d的气体浓度记为C（x,y,z,t ）.根据假设，单位时间通过单位法向面积的流

9、量q =k gradek是扩散系数，grad表示梯度，负号表示由浓度高向浓度低的地方扩散，考察空间域门， 门的体积V，包围门的曲面为S，S的外法线向量为n，则在t,t*t内通过门的流量 为Q!qn d；dts而门内气体的增量为QC（x, y, z,t）-C（x,y, z,t t）dVV由质量守恒定律Qi根据曲面积分的奥氏公式q* nd；二 divqdVs其中div是散度记号。C kdiv(gradC) t(1)V由以上公式再利用积分中值定理不难得到2 22c C c C c C二 k( 222),t 0,：x,y,zx2:y2z2这是无界区域的抛物线型偏微分方程。函数，可记作C(x, y,z,

10、O)二 Q (x,y, z)(2)Q表示泄漏源泄漏施放的气体总量， 方程(1)满足条件(2)的解为Q根据假设1,初始条件为作用在坐标原点的电源(x, y, z)是单位强度的电源函数。x 2 y 2 z24 ktC ( x , y, z, t)3/ 2(4 兀 kt )这个结果表明，对于任意时刻t烟雾浓度球面半径R的增加C的值是连续减少的；当5.2问题(2)气体在泄露过程中，气体从裂口泄漏的速度与流动状态有关，计算泄漏量时首先要判断泄漏时气 体的流动属于声速还是亚音速流动，根据下式判断：rP£Po(r+1/2)r一一(1) 当成立时，气体流动属亚声速流动，则泄漏流量为：C的等值面是球面

11、x y2z R2,并且随着Rr '或 tr :时 C(x, y, Z,t)r 0(1 1)c式中：Q气体泄漏流量，kg/s； Cd 排放系数，通常取1.0； A 泄漏口面积，m2 p 容器内气体压力，Pa;尺一环境压力，Pa;丫 绝热指数，是等压比热容与等容比热容的比值；M 气体的分子量，kg/mol ； R 气体常数，8.314J/(molK)；T 容器内气体温度，K。1p-p0(1/(2 )当2成立时，气体流动属声速流，则泄漏流量为：Q = Cd pAMRT(1 2)式中各变量意义同上本问中，我们将核泄漏视为高架点连续点源扩散模型，即连续源或泄放时间大于或等于扩散时间。 而且本题中

12、气体恒以 m kg/s的速度排出，故泄露流量 Q亘定。(3)设泄漏源有效高度为H,取其在地面投影为坐标原点，x轴指向风向。考虑地面 反射作用，可得到烟羽模型的浓度分布算式为：y2 (z h)2C (x, y ,乙 t)eCFZ(z H )22-1式中：C 污染物浓度, 泄漏高度的平均风速， 轴上的扩散参数；H y、 z 侧风方向、-源强即气体泄漏流量或速度，y二-分别用浓度标准偏差表示的m ； x -下风方向到泄漏源点的距离,kg / mm / s ；-泄漏有效高度，垂直向上方向的离泄漏源点的距离，m okg / sy7轴及（4）扩散参数的确定应用高斯模型的关键是确定扩散参数。根据定义，在均匀

13、湍流场中，扩散参数与下风向，距离的关 系是明确的。确定扩散参数的方法一般有：P-G扩散曲线、P-T法、世界银行推荐扩散系数法、sminth法等。本文采用P-G扩散曲线确定扩散参数法。P-G扩散曲线是帕斯奎尔在大量观测和研究的基础上，于1961年总结提出一套根据常规气象观测资料划分大气稳定度级别和估算扩散参数的方法。为了便于使用计算机计算大气污染物浓度分布，可用幕函数式近X _1- -x _ 2似表示P-G扩散曲线，将 7，'- z表示为下风距离 x的函数， 7 一 1x ，、一 z- 2X在进行计算时，首先从表 1-2中确定大气稳定度，我国在标准GB/T13201-1991中规定，当确

14、定稳定度级别后，实际的扩散参数选择见中国有关环境评价标准采用的系数值-P-G扩散曲线幕函数数据。表1-2帕斯奎尔稳定度级别划分表5.3问题（3）在考虑风速的条件下，由于我们在前面的假设中认为风向与x轴相同，故我们要求计算上风和下风I处的放射性物质浓度，只需在建立的模型下，分别用（k+s）和（k-s ） 去代替方程中的u，同时分别令x=l,y=z=0或者x=-l ，y=z=0,这样就可以求出当风速一 定时，上风和下风l公里处，放射性物质浓度的估计模型。上风l处浓度计算式如下：m_ -H 2ie 2k-s； 2 J下风l公里处浓度计算式如下:-H2eC(x, y,乙 t)二 s(3 1)C(x,

15、y,乙 t)二2叮根据扩散系数的定义可以知道,（3 2） l -2，z 一 2l，代入上两式就可以求出上风和下风I处的浓度5.4问题（4）由于本问中，我们只需要通过检测福岛的核泄漏是否对我国东海岸以及美国西海岸有影 响，故我么只需要在前面建立的模型的基础上，通过网上查询相关的核泄漏在我国东海 岸以及美国西海岸的风向，风速以及之间距离，代入建立的模型中，并用matlab进行仿真模拟即可。六模型求解和结果分析5.1问题一求解和结果分析：我们通过将中广核以前的测试数据（如下图5.1.1 ）进行仿真模拟后，得出实际值与预测值的图形（5.1.2 ）如下：通过图形，我们可以看到，模拟值与实际值拟合的比较好

16、， 从而验证了模型的准确性，同时我们可以看到，随着离泄漏源距离的延伸，最终放射性 物质的浓度越来月小，趋近于零，即当l趋向无穷是，c（x,y,z,t） 趋向于零。l/m100300500700900c(kg/m3)0.0029540.000818260.000361010.000205150.00013352l/m11001300150017001900C(kg/m3)9.656e-57.514e-56.0584e-55.0186e-53.8988e-5距离泄漏源不同距离处的放射1性物质浓度（图5.1.1 ）5.2问题二求解和结果分析：当风速为k m/s时，我们根据上面的高斯烟羽模型，将相对速

17、度代入到式中的u,即可得到核电站周围放射性物质浓度的变化情况。假设风速k=2.1m/s ;泄漏源强度Q=1kg/s；地面粗糙度参数Z0=0.4 ;计算精度d=1m由matlab仿真结果可得如下图所示的结果：图5.2.1为上风向浓度分布，图522为下风向浓度分布。由两图我们可以直观的看出核电 站周边地区的放射性气体的浓度变化规律。该程序中的源强、风速、大气稳定度参数、 地面粗糙度参数和计算精确度等都可根据实际晴况分析需要设置。inl x|£ile £dit 匝艸Xcoli Uzktop Window Help爼 q is ef亦爲m背運)逅応匡斤戸x 10'1(

18、63;<逼叵旻o气体扩歆上风向浓度分布图-r1I fIiI150200250300350400450500x轴冋距离(km)图5.2.1上风向浓度等高分布曲线Eile Edit yi ew Insert Tools Hesktap Liiidow HelpL>005o11001502002603003504004S0500浩由向距离(km)图522下风向浓度等高分布曲线5.3问题三求解和结果分析：我们通过对第二问的分析可以从图中得出在距上风处和下风处L公里时可以对应的得出相应地点的气体浓度，方便直观。比如当L=100km时，通过观察上风和下风浓度等高曲线可以很直观，很方便的得出：在下风处得浓度是千分之四；在上风处时千分之五。而 在实际中，则要根距不同的实际情况，确定好各个参数，更准确的得出相应点的气体浓 度。我们通过查询大亚湾核电站周围放射性物质的浓度数据用matlab进行仿真模拟，得出实际的图形与模拟的图形如下（5.3.1 ）和（5.3.2），通过观察图形我们可以知道， 两图拟合较好，说明以上所建立模型的准确性较高。图 5.3.1图5325.4问题4求解和分析七模型的平价与推广Gaussian模式是在大量实测资料分析的基础上，应用湍流扩散的统计理论得到的正态分 布假设卞的扩散模式、采用正