气体扩散模型论文

1、1 放射气体模型的预估模型 摘要 本文是以日本福岛核电站遭遇自然灾害发生核泄漏的背景而提出的。 对丁问题一，考虑到放射性物质的泄漏是连续不断的。 本文根据“泄漏放射 性物质质量守包定律”和“气体泄漏连续性原理”建立了微积分方程，应用了高 等数学中散度、梯度、流量等数学概念，通过Guass公式、四维二阶偏微分方程， 因而得到了核电站周边不同距离地区、 不同时段放射性物质浓度的预测模型。 同 时为使模型适用范围更广，本文引入了地面反射系数，考虑了由丁放射性物质从 泄漏口喷出时具有初动量而使其泄漏源有效高度被抬高等因素， 进而得到了在无 风环境中适用范围更广的“高斯修正模型”。 对丁问题二，要探究风

2、速对放射性物质浓度分布的影响。本文运用概率学 1 知识，通过图解和数学推导得出“连续点源放射性物质高斯扩散模型”。本文依 次考虑了 “重力沉积”、“雨水沉积”、“核衰变”等因素对浓度分布的影响。 并通过构建“耗减因子”、“衰变因子”等方法将耗减和衰变的放射性物质“投 影”到泄漏源浓度中，得到了经多次合理修正后的“优化高斯模型”，并据此分 析了泄漏源周边地区放射性物质的浓度变化。 针对丁问题三，本文在问题二的基础上，结合考虑风速和放射性物质扩散速 度在空间中的欠量运算。得出在对上风口分析时，要分类讨论风速和自然扩散速 度之间的大小关系，当风速小丁自然扩散速度时，放射性物质是无法到达上风口 的。

3、对丁问题四，本文参阅整理大量气象、地理、新闻资料，选择我国东海岸典 型地域-山东半岛作为研究对象，综合考虑对应海域平均风速及风向、地理距离、 海水对放射性物质扩散的部分反射系数等因素， 并通过C+编程模拟计算，预测 出放射性核物质将经过 6.5天到达我国东海岸，且 131131I I浓度预测值为： 0.100mBq m与实际情况比较吻合 2 关键词：放射性气体 扩散 浓度变化 高斯修正模型 预测 1问题的提出 由于重大的突发性核泄漏紧急灾害事件具有爆发性、 空间分布不连续性、对 周边地形和气象条件的敏感性的特点，研究核事故所释放的物质的时空分布需要 高度精确的技术，但是在对于更好地保护环境有着

4、极其重要的意义。 在有一座核 电站遇自然灾害发生泄漏，浓度为po的放射性气体以匀速排出，速度为m kg / s , 在无风的情况下，匀速在大气中向四周扩散，速度为 sm/s。 问题一，若能建立一个描述核电站周边不同距离地区、 不同时段放射性物质 浓度的预测模型，这对于研究核污染模式具有重要的意义。 问题二，当风速为km/s时，给出核电站周边放射性物质浓度的变化情况， 这对于研究核电站附近浓度的在实际环境下有着重要的作用。 问题三，当风速为k m/s时，计算出上风和下风L公里处的放射性物质浓度 的预测模型就显得更加急迫的了。 问题四，将建立的模型应用于福岛核电站的泄漏，计算出福岛核电站的泄漏 对

5、我国东海岸的影响，这一实际的有意义预测，可以明确我们实际的污染情况， 为我们的核应急决策提供技术支持。 2问题的分析 对于问题一，在无风的情况下，放射性气体s以sm/s的速度，匀速在大气中 向四周扩散。在此条件下，探求一个模型来对核电站周边不同距离地区、不同时 段放射性物质浓度进行预测。我们要明确此问题研究的核扩散是点源连续泄露的 扩散问题。虽然只是要求考虑在无风情况下放射性物质浓度分布，但为了使模型 更贴切实际，需考虑地面反射、泄漏源有效高度等因素对浓度分布的影响。 根据“泄露放射性物质质量守恒定律”和“气体泄漏连续性原理”进行分 析，发现要得出核电站周边不同距离地区、不同时段放射性物质浓度

6、的预测模 型，最后对于该方程进行分析求解。 对于问题二，为了探究风速对发生核泄漏的核电站周边放射性物质浓度分 布的影响，运用概率学知识，通过图解和数学推导得出“连续点源放射性物质 高斯扩散模型”。应在“连续点源放射性物质高斯扩散模型”的基础上经多次合 理修正后得到更好的“优化高斯模型” 对丁问题三，该问题要求建立泄漏源上风口和下风口处放射性物质浓度的 预测模型，在参3 考第二问的基础上，主要考虑风速和放射性物质扩散速度在空 间中的欠量运算。在对上风口分析时，要分类讨论风速和自然扩散速度之间的 大小关系，当风速小丁自然扩散速度时，放射性物质是无法到达上风口的。 对丁问题四，应参考大量气象、地理、

7、新闻资料，选择我国东海岸典型地 域-山东半岛，作为研究对象，综合考虑对应海域平均风速及风向、地理距离、 海水对放射性物质扩散的部分反射系数等因素，预测出放射性核物质与实际情 况比较。 3模型的假设 考虑到放射性气体扩散的复杂性，为简单起见，在讨论扩散模型时都作了如 下假设； (1) 瞬时泄漏假定瞬时完成，连续泄漏假定泄漏速率包定； (2) 气云在平整、无障碍物的地面上空扩散； (3) 气云中不发生化学反应，地面对气云无吸收； (4) 为水平风向，风速和风向不随时间变； (5) 气体的传播服从扩散定律，即单位时间通过单位发祥面积的流量与他的浓 度梯度成正比； (6)气体的扩散看作空中末已连续店员

8、向四周等强度瞬时释放气体， 放射性气体在无穷空间适房的过程不发生性质变化 4符号说明及名词解释 4.1符号说明 符号 说明 s 放射性气体的传播速度 k 风速，单位m/ s H 泄漏点。距有效地面的高度 t 任息扩散时刻 C(x, y,z,t) 空间任意一点的放射性物质浓度 4 、i(i =x,y,z) 空间任意一点的放射性物质的扩散系数 Q 空间域 V 空间域其体积 S 一规则的球面而枳 QI 在(t,t+At)(t,t+At)内通过 Q Q 的流量 Q2 Q内放射性物质的增量 符号 说明 QO 从泄漏源泄漏的放射性物质的总量 .h 附加高度 Ts 核泄漏出口处的温度 TO 环境温度 Ct

9、设地面反射系数 Q 源强，单位为kg / s 二 x, y, ；z 分别为用浓度标准差表小的x, y,z轴上的扩散参数 Vs 沉降速度，单位为m/s Wd 地面干沉积率 W 冲洗系数 T0.5 放射性核素的半哀期 4.2名词解释 烟羽乂称烟云(smoke cloud)、烟流(smoke plume)：从烟囱中连续排放到 大气中的烟气流。由丁烟羽各部分的运动速度不同，因而其外形也千变万化。不 同的烟羽形状表示污染物浓度的空间分布不同。 它与大气湍流、大气稳定度、地 形地物、排放参数等有密切的关系。 动力抬升：暖气流受锋面、辐合气流的作用被迫上抬，或者在运行中受地5 形阻挡产生上升运动，这种空气在

10、运动中由外力(不包括重力和浮力)使一部分 空气被抬上升。 湍流扩散：是指湍流运动导致大气或水体中的污染物质或其他物质与周围 洁流体的混合。 5模型的建立和求解 5.1问题一模型的建立与求解 5.1.1模型一的建立 以核泄漏点正下方的地面为坐标原点(0,0,0)，平均风向为X轴、指向下风 方向，铅直方向为Z轴，水平垂直丁风向轴(X轴)为Y向，建立空间坐标系， 则核电站泄漏点O距有效地面的高度为H，则泄漏点位置坐标为0(0,0, H)。 图1空间坐标系示意图 并记t时刻时，空间任意一点的放射性物质浓度为 C(x,y,z,t)。根据假设设 单位时间通过单位法向面积的流量 与浓度梯度成正比，有： q

11、= q gradC (1) (i=x,y,z)是扩散系数，grad表示梯度，负号表示由浓度高向浓度低的地 方扩散。 先考察空间域。,其体积为V，包围。的曲面为S , S为一规则的球面,S 外法线向量为=(-=-,1)。 z z 则在(t,t (t,t +A+At)t)内通过 Q Q 的流量为： - Q Q1 = t ! ! q ndq nd。dtdt(2) 6 Q内放射性物质的增量为: Q2 = C(x,y,z,t ：t)-C(x,y,z,t)dV V 从泄漏源泄漏的放射性物质的总量为: 上述模型仅是一个最理想化的预测浓度模型5,因为它将环境视为无边界空 问，且没有考虑放射性物质从泄漏口释放出

12、时的初动量。 为使所建立模型更加贴 合实(3) (4) 根据质量守包定律和连续性原理，单位时间内通过所选曲面S的向外扩散的 放射性物质与S曲面内放射性物质增量之和，等丁泄漏源在单位时间内向外泄漏 的放射性物质。有: Qo = Qi Q2 (5) t - t 即， C(x,y,z,t ：t) -C(x,y,z,t)dV .q nd；dt =； podVdt (6) 乂根据曲面积分的Gauss公式3： iiq nd；：= divqdV （其中div是散度记号） (7) C(x,y,z,t ：t) C(x, y,z,t) tdVI iiidivqdVdt= t 111 podVdt ,t 、 . 巴

13、=肺324 pcmy/t) .:t 0 t kdiv(gradC )dt =lim .t)0 由以上两式得: :C IIIdV t divqdV t = podV t ft V 即为: ；C I!dV “ i l l divqdV (8) 根据参考文献 M =、 ;:t x 2 :c (9) 解得：C(x,y,z,t) Po i.5/ . . .、o.5exp (4 t) (xyz) (z-H)2 .4、xt 4、yt 4、zt (10) 7 际，本文先从“有效泄漏源”和“地面反射”这两个方面对所建立模型进行 修正 5.1.2考虑热力抬升作用对模型的修正 如图2所示，H为核泄漏点源的有效源高。

14、它是由两部分构成：一是核泄漏 口的有效高度h；二是在实际核扩散中核泄漏气团从泄漏口排出时，由丁受到热 力抬升和本身动力抬升，进而产生的一个附加高度 Ah。因而H =h + Ah。 图2热力抬升示意图 对丁 Ah,主要由浮升力和泄漏的初始动量决定，同时还要受到泄漏口温度、 大气温度、风速、地形地貌等多种因素的影响。我们直接引用气体污染扩散学中 应用较广范的，有关烟气抬升高度的综合分析公式6: i (0.92VsD 5.25Fb0.4h0.6) h= - (11) Us 式中：h泄漏源的实际高度； us泄漏源出口处的风速，已知为km/s； D泄漏源出口的有效直径； Vs放射性气体的扩散速度，已知为

15、s m/s； Fb浮力通量,m4/s； 由Briggs抬升公式知，浮力通量Fb = fD I Vsg TTa , 2 - Ts 其中Ts一核泄漏出口处的温度K ; To 一环境温度K。 a：在有风(/ (/ 21.0m/s)且释放气体温度与环境温度差 35K(35K(； - -Ta占35K ) )时，抬 升高度： ffi念*:泅K8 5.1.3考虑到地面反射对模型的修正 ,(0. 9/2D 0. 7Q)2h0.6 ) _h = - Us (12) 其中Q =0.275PD2Vs , P为大气压强hPa To b:b:小风（us 35K 35KTs-Ta芝35K ）, ,抬升高度为: ：h=5.

16、5Q .25 dT0 0.0098 dZ _0.375 (13) 其中 匹 为泄漏源的有效高度处上的环境温度梯度（K /m） dZ c：当温度差 35KTsT0 圣35K k idT 矿家5 h = 5.5Q .25. 0 +0.0098 ，适用丁k Ts T 芝 35K | IdZ ） 2（1.5sD + 0.01Q） ebh 才用, * - - ，适用丁k不限（k1.0时，取k=1.0） , Ts -T0 35K 注：H为泄漏源的有效高度；k为泄漏源出口处的风速；D为泄漏源出口的有效直径; 胡放射性气体的扩散速度；Ta核泄漏出口处的温度（K） ; T。环境温度（K）; Q =0.275PD

17、2V TzT0 s T （P为标准大气压） 综上所述，泄漏源的有效高度为： H = h hH = h h (15) 9 考虑到地面会对扩散来放射性气团有一定的反射作用， 同时扩散的核素粒 子受沉降等作用，核扩散物质乂不可能被全部反射回去，因而核辐射物质只能是10 部分反射回大气，为了便丁描述和模拟，设地面反射系数为 a 。这样进入大气 的核辐射物质可以看成是两个部分：一是从泄漏源 O直接扩散到空间A点；二 是从地面反射进入空间A点(见图3) 图3连续点源扩散地面部分示意图 从几何物理学分析可知，通过地面反射进入大气的核扩散物质相当丁虚拟 一.八、一、一 . . . . . . . . . 、

18、. . . . 泄漏源 O O 泄露的核辐射物质在原来空间的一个浓度叠加。如果所设核泄漏点源 。在距有效地面的高度为H的地方，本文求空间任一点 A(x,y,z)的浓度值，则 实际泄漏源对A点的影响部分可用exp (Z H)来表示。在考虑反射系数a后， 2；弓2 虚泄漏源。对A点的影响部分可用口，exp-三斗来表示，丁是(10)式所得模 2七2 型可以修正为： 2 C(x,y,z,t)= 人.5；、不、0.5exp4- ， (Ot)(知 8y5z) I 4知t 、H =h十凯 5.2问题二模型的建立与求解 5.2.1.1模型二的建立 问题二，针对在有风速km/s时，选取高斯连续点源扩散模型进行分

19、析，这 也是在所有有风速时，放射性物质扩散模型中最常见也最方便的一种泄漏方式。 考放射性气体云团在大气中实际迁移和扩散的数值计算基本上可分为二步：(16) 虚源o 高架核泄漏建睇点源扩散地面部分反射图 ：exp 11 第一步，根据大气动力学理论进行所关注区域中风场的计算， 其理论基础是 大气运动方程、连续性方程、状态方程、热力学方程和水汽方程构成的基本方程 组； 第二步，进行已知风场中放射性气体云团迁移和扩散的计算， 可采用类似丁 处理大气污染的方法，假设放射性气体云团不影响大气流体速度和温度， 求解出 放射性气体云团的连续性方程。 高斯连续点源扩散模型有边界点源，这里设泄漏源有效高度为 H

20、,取其地 面投影为坐标原点，x轴指向风向，继续考虑地面反射作用， 图4边界点源示意图 可得到高斯连续点源泄漏的浓度分布7: 其中，Q为源强，kg/s ; c为污染物质量浓度，kg / m3 ； bx,cry,crz分别为 用浓度标准差表示的x, y,z轴上的扩散参数；k为泄漏高度的平均风速，m/ s； H 为泄漏有效高度，m o 考虑到动力抬升的作用，将有效高度（即有效烟云高度） He = H +AH 。c(x,y,z,H) Q 2 二 EL exp( 2 y , )exp 2 二 J y (z-H)2 exp (z H)2 2二2 (17) 12 热力抬升阶段 S 瘠 7 /7/ 7 777

21、 *77/ /，77 / / / 图5动力抬升示意图 因此（17）式可以写成为: 2 2 2 expJ-Xexp _ z - ： +exp一 (18) 2- y 2；一 z 2一 z 5.2.1.2对高斯浓度计算模型的修正 利用标准的高斯模型，计算大气中放射性核素云团的扩散，目前已有较多的 研究，而实际的核素扩散过程还存在着粒子的重力沉降、 雨洗沉积以及核素衰变 等的因素对浓度分布的影响。因此，有必要对高斯浓度扩散模式进行修正， 方可 较为准确、真实地反映实际的核素扩散规律。 由丁地面风速对大气扩散的影响起着至关重要的作用， 不同的风速其浓度计 算方法有较大差异，因此，下面的分析从有风（大丁

22、1.5 m/s）、小风（小丁 1.5m/s 大丁 0. 5m/s）和静风（小丁 0. 5m/ s）,三种情况分别论述。 a：考虑干沉积时的连续点源扩散 粒径大丁 10 m的粒子有明显的重力沉降，粒子的沉降速度取决丁空气阻力 和重力平衡，可用斯托克斯公式表示: Vs* 式中，P :粒子密度，kg/m3 ； g :重力加速度，9. 806 5m/s； D :粒子直径，m ; H：空气的动力粘性系数,可取1.8乂10* kg/（m司；互介阶段 扩散阶段 Q c(x,y,z,H)= 2 二 EyOz (19) He 13 Vs :沉降速度，m/s,含碘放射性核素的干沉积速度为 Vs= 1. 1 cm/

23、 s,该 值也可由式()计算得到。 因为在扩散过程中同时有重力沉降的位移迭加到羽流中心线上， 中心线就会 向下倾斜， 所有粒子相当丁在下倾的中心线上扩散。该类扩散和沉降的迭加可 认为是羽流运行过程中，实源以 也的速度向下移动，在x处向下移动的高度为 ”=告，即源高由H降到了 H岑。实际上，由于大气湍动及其他动力作用, 地面不是全吸收表面，应考虑地面的反射作用。但毕竟存在粒子的沉降作用，乂 不是全反射，因此，需对反射项乘以一个反射系数 a(a1),由丁反射项的有效 源高度也变成了 H -VsX,故相应的浓度计算公式为: u y2 (z-H )2 (z H exp(-土)exp - exp 2 二

24、 k；z 2；y 2了 z 式(20)可得到经过干沉积作用后的浓度分布结果。式 (20)中，反射系数口需外 部给定，a(a 1),通常对丁放射性核素可取0.5。 b:b:小风(小丁 1.5m/1.5m/以丁 0. 5m/s)0. 5m/s)和静风(小丁 0.5m/ s)0.5m/ s)时的计算 参考小风和静风时的浓度计算方法， 再利用部分反射倾斜烟云模式，可得出 该情况下的浓度分布： Q (x-UT)2 y2 (zH 0T)2 】 (z+H-VsT)2 (21) c(x,y,z, H) = - 3/2 - exp(- - )exp( _ J)exp - 2 二exp_ - 2 dT k 1 。

25、(2二)；二；3；二 2馈 2;y 2;z 2. 通常按照每小时排放6个烟团进行积分，即可达到需要的精度。 c:c:地面的干沉积量 根据扩散理论和动量传递的普朗特理论，可以导出干沉积的地面沉积量为: Wd =乂 c(x,y,0) 式中，c(x, y,0)前面各种情形计算的地面污染物的浓度; Vs干沉积速度，m/s Wd地面干沉积率，Bq m s-1。 d:d:雨洗作用 降雨对烟羽中的颗粒物及气溶胶具有活洗作用， 可溶性气体与蒸汽亦可溶丁 雨水中，降雨过程造成的这类湿沉积是导致放射性气溶胶和气体向地面沉积的另VsX)2 _r) (20) 2z Q c(x, y, z, H )= (22) 14

26、一重要机制。通常以冲洗系数中(S)，描述降雨对烟羽中污染物活洗作用的大小。 与雨强的关系可以表达为： 里=3 (23) 式中，I为雨强(mm/ h); a, b为经验系数。 式中，按释放物质为含碘、不含碘情况分别取值。对含碘物质，取 a =8乂10七b= 0.6;对于不含碘物质，取 a = 1.2x103,b = 0.5。 对丁湿沉积导致的烟羽耗减，可采用湿沉积耗减因子对源强 Q进行修正，有 - - x Q(x) =Qexp( = ) (24) u e：放射性衰变的影响 除了干沉积和湿沉积外，放射性物质的衰变也是影响大气中核素浓度分布的主要 原因。由丁放射性物质服从简单的衰变规律，其浓度随时间

27、的变化可由下式计算： c=c0eT (25) 式中，C0为初始浓度； 赤衰变常数； t经过的时间； C0可由无衰变的浓度公式计算。 参照前面的模式，同样采用衰变耗减因子来对源项进行修正。 经过理论推导，可 以得出： 其中，TO.5放射性核素的半衰期。 如以上情况都存在，应同时考虑源项修正，即式 (24)和式(26)相乘后代替Q 5.2.2模型二的求解 5.2.2.1模型参数的确定 在事故过程模型中所需参数的选取及确定十分重要， 通常情况下气象参数的 选取是利用该地区多年气象资料，采取工业安全与环保统计的方法进行有关参数 的确定，而其他扩散参数是以实际测定为基础的。 a:a:大气稳定度的计算 根

28、据国家标准（GB/ T 13201 - 1991）制定地方大气污染物排放标准的技术 Q(x) = Qexp( 0.693x 3600T0.5u (26) 15 方法的规定，划分大气稳定度的级别，共分为6级A F, A为极不稳定;F为极稳 定。首先，根据释放源所在地的经度和纬度以及泄漏的日期和时间计算当时的太 阳高度角h0;然后，由太阳高度角h0和云量查出太阳辐射等级；最后，再根据地 面风速确定当时的大气稳定度，计算细节可参考文献。 表1 大气稳定度的级别参考表 地面风速 白天太阳辐射 阴天的白天或夜 有云的夜晚 （ms） 强 中 弱 间 薄云遮天或彳氐云 M 0.5 0.4 6 C D D D

29、 D D 有风时的扩散参数，（QyQz）的确定采用Briggs给出一套扩散参数，如表1和 表2所示。 表2 Briggs扩散参数（开阔平原田野） 大气稳定度 CT y A 0.22x(1 +0.0001x)12 0.2x B 0.16x(1 +0.0001x)12 0.12x C 0.11x(0.0001x)12 0.08x(1 + 0.0002x)2 D 0.08x(0.0001x)12 0.06x(1 + 0.0015x)12 E 0.06x(0.0001x)12 0.03x(1 0.0003x) F 0.04x(0.0001x)12 0.016x(1 0.0003x) 表3 Briggs

30、扩散参数（城市） 16 大气稳定度 -x -y A-B 1 0.32x(1 0.0004乂)一2 0.24x(1 0.0001刈2 C 0.22x(1 0.0004x)-12 0.20 x D 0.16x(1 0.0004x)-12 0.14x(1 0.0003x)-12 E-F 0.11x(1 0.004x)12 0.08x(1 0.0015x)12 5.2.2.2模型的求解 图6模型二计算流程图 5.3问题三模型的建立与求解 5.3.1模型三的建立 17 问题三，针对在有风速km/s时，选取高斯连续点源扩散模型进行分析在给 出的上风和下风L公里的距离处在结合，在模型二的基础上考虑浓度等效圆

31、的情 况。在模型二里面将式（18）中的风速k在上风和在下风不同情况下与传播速度 s之间的比较的分析。 a:a:在上风的情况下，此时满足（km/s芝sm/s） 此时若km/ssm/s时，即风速小丁自然扩散速度，此时的放射性物质是 无法到达上风口的，则必须满足风速大丁自然扩散速度， 仍然根据模型二可以求 解，只是将式（18）中的k用k=ks来代入。 b:b:在上风的情况下，此时k无较大的影响和要求 此时不论风速小丁自然扩散速度，还是风速大丁自然扩散速度，此时的放射 性物质都可到达下风口的，则仍然根据模型二可以求解，步骤同 5.2.2. 模泄溺模拟图 4Q 时 80 100 120 甘 Q 160

32、13C 网 Q 图7核泄漏对上下风在距离上的影响 5.4问题四的求解 日本丁 2011年3月15日下午，世界气象组织和国际原子能机构北京区域环 境紧急响应中心分析认为：日本中北部区域在中低层大气中的风向由西南风转为 西北风；高空大气主要以偏西风气流为主， 近期由丁降水发生，有利丁核物质沉 降，影响范围缩小。未来三天（16日至18日），日本核电站核泄漏产生的放射 性污染物主要影响区域为日本中部、北部及其以东的北太平洋区域。 在结合模型二、模型三的条件下，在参阅整理大量气象、地理、新闻资料， 选择我国东海岸一一山东半岛作为研究对象，综合考虑对应海域平均风速及风 向、地理距离、海水对放射性物质扩散的

33、部分反射系数等因素， 并通过C+编程 模拟计算，预测出放射性核物质将经过 6.518 天到达我国东海岸一一山东半岛，且 131131I I浓度值为：0.100 mBq，mm与，与实际情况比较吻合。 6结果分析与检验 由丁在前四问题中均无具体的数据，所以对丁结果，在检验上存在一些麻烦， 但是在模型三预测上下风的浓度具有较好的实际意义， 而在模型四中，对丁预测 出浓度的变化和对我国的影响也较为合理， 故该模型具有可靠性和参考性，可以 进行推广和其它方面的应用。 7模型的优缺点与改进方向 优点：该方法可计算核事故中连续点源和瞬时点源在不同气象、 地形条件下 的浓度分布。该结果在核事故的应急救援过程中

34、， 对救援人员划定警戒区和确定 周围居民的疏散范围具有重要意义，并可为制定救援方案和应急决策提供科学依 据。放射性云团在空中迁移和扩散提供浓度的定量描述，对可能发生的核事件的 放射性核素浓度监测及监测时间范围提供相关信息。 缺点：模型在建立时只考虑了两种下垫面情况， 即农田和城市，实际的扩散 情况可能因地形条件的复杂化而有所改变。由丁扩散是一个非常复杂、影响因素 众多的过程，上述方法还有一定的局限性，有待丁进一步检验。另外，该方法只 给出了浓度的计算结果，没有换算成人的吸收剂量，今后将对其计算过程和参数 的选取进行多方面验证，使其进一步完善。 参考文献 1 金勇进，统计学，北京：中国人民出版社

35、，2010。 2 梅长林，周家良,实用统计方法，北京:科学出版社,2002。 3 罗艾民，魏利军.有蠹重气泄露安全距离数值方法 J.中国安全科学学报， 2005,15 (8) : 98-100 4 周燕等,Matlab在统计与工程数据分析中的应用，北京:电子工业出版社，2010。 5 丁信伟,王淑兰，徐国庆.可燃及蠹性气体泄漏扩散研究综述 .化学工业与工 程,1999,2(16):11 122. 6 胡世明.气体释放源的三维瞬态重气扩散研究.劳动保护科学技 术,2002,3(20):2830. 7 蒋军成，潘旭海.描述重气扩散的一种新模型.南京工业大学学报,2002,1(24):41 46.1

36、 附录 附录一：模型二在上下风 L L公里的浓度求解源代码 #include #include #include using namespace std; #define e 2.718281828459 #define PI 3.1415926 void qiuxy(char wending ,double x,double &dy,double &dz) switch(wending) case A: dy=0.22*x/(sqrt(1+0.0001*x), dz=0.2*x; break; case B: dy=0.16*x/(sqrt(1+0.0001*x), dz=0

37、.12*x; break; case C: dy=0.11*x/(sqrt(1+0.0001*x), dz=0.08*x/(sqrt(1+0.0002*x); break; case D: dy=0.08*x/(sqrt(1+0.0001*x), dz=0.06*x/(sqrt(1+0.0015*x); break; case E: dy=0.06*x/(sqrt(1+0.0001*x), dz=0.03*x/(1+0.0003*x); break; case F: dy=0.04*x/(sqrt(1+0.0001*x), dz=0.016*x/(1+0.0003*x); break; void main() double c1,c2,x,z,t,p0,H,dx,dy,dz,a,s,k,I,Ts; char daqi; cout请按以下格式输入: endl; cout”下风口 L 处坐标(米)：x; c