放射气体模型的预估模型摘要本文是以日本福岛核电站遭遇自然灾害发生核泄漏的背景而提出的。且结合 了高斯烟羽模型、线性拟合，以及微分方程模型，运用MATLAB软件，分析泄 漏源强度、风速、大气稳定度参数等因素对放射性气体扩散的影响，预测了放射 性气体浓度在不同时间，不同地区的浓度变化，并且本文模型中数据可以根据不 同的实际情况而加以改变，因而是本文的应用范围大大增加，可以适用于具有较 强的应用型。对于问题一，讨论在无风的情况下，放射性气体以 s m/s 的匀速在大气中向四 周扩散。本问中由于不考虑风力的影响，且扩散出来的气体匀速向四周散开，这 样经过任意时刻t,扩散的气体围成一个半径为st的球，且距球心位置不同的地 方浓度值不同。采用列数列的表现方法，设定相同时间段t,把条件进行整理， 并经过简单计算得出每段时间所预测得到的扩散距离r和浓度Co利用MATLAB 软件对数据进行线性拟合，采用微分方程模型得到核电站周边放射性气体在不同 地区，不同时间段的浓度变化，得出随着离泄漏源距离的延伸，最后放射性物质 的浓度越来越小，趋近于零，即当x趋向无穷时，C(x,y,z,t)趋向于零；当时间趋 向于无穷时，C(x,y,z,t)也趋于无穷。对于问题二，要探究风速对放射性物质浓度分布的影响。风速的处理是问题的中 心，采用大气污染的经典高斯扩散模型,实现了高斯烟团气体扩散模型的动态预 测，分析计算了气体扩散过程中的各关键参数。对于问题三，本文在问题二的基础上，结合考虑风速和放射性物质扩散速度在空 间中的矢量运算，将在上风和在下风不同情况下与传播速度s之间的比较的分 析，利用高斯烟羽模型对核电站周边地区的浓度进行预测，然后，利用 MATLAB 软件，将相关数据代入程序，我们得到核电站周边地区的浓度分布的等高曲线。对于问题四，本文参阅整理大量气象、地理、新闻资料，选择我国东海岸典 型地域-山东半岛作为研究对象，综合考虑对应海域平均风速及风向、地理距离、 海水对放射性物质扩散的部分反射系数等因素，集合核电站周边的浓度等高线， 可。关键词：放射性气体 扩散 浓度变化 高斯修正模型 预测1 问题的提出由于重大的突发性核泄漏紧急灾害事件具有爆发性、空间分布不连续性、对 周边地形和气象条件的敏感性的特点，研究核事故所释放的物质的时空分布需要 高度精确的技术，但是在对于更好地保护环境有着极其重要的意义。在有一座核 电站遇自然灾害发生泄漏，浓度为p0的放射性气体以匀速排出，速度为m kg /s , 在无风的情况下，匀速在大气中向四周扩散，速度为s m/s。问题一，若能建立一个描述核电站周边不同距离地区、不同时段放射性物质 浓度的预测模型，这对于研究核污染模式具有重要的意义。问题二，当风速为km/s时，给出核电站周边放射性物质浓度的变化情况， 这对于研究核电站附近浓度的在实际环境下有着重要的作用。问题三，当风速为km/s时，计算出上风和下风L公里处的放射性物质浓度 的预测模型就显得更加急迫的了。问题四，将建立的模型应用于福岛核电站的泄漏，计算出福岛核电站的泄漏 对我国东海岸的影响，这一实际的有意义预测，可以明确我们实际的污染情况， 为我们的核应急决策提供技术支持。2 问题的分析对于问题一，在无风的情况下，放射性气体s以sm/s的速度，匀速在大气中 向四周扩散。在此条件下，探求一个模型来对核电站周边不同距离地区、不同时 段放射性物质浓度进行预测。我们要明确此问题研究的核扩散是点源连续泄露的 扩散问题。虽然只是要求考虑在无风情况下放射性物质浓度分布，但为了使模型更贴切实际，需考虑地面反射、泄漏源有效高度等因素对浓度分布的影响根据“泄露放射性物质质量守恒定律”和“气体泄漏连续性原理”进行分析，发现要得出核电站周边不同距离地区、不同时段放射性物质浓度的预测模 型，最后对于该方程进行分析求解。对于问题二，为了探究风速对发生核泄漏的核电站周边放射性物质浓度分 布的影响，运用概率学知识，通过图解和数学推导得出“连续点源放射性物质 高斯扩散模型”。应在“连续点源放射性物质高斯扩散模型”的基础上经多次合 理修正后得到更好的“优化高斯模型”。对于问题三，该问题要求建立泄漏源上风口和下风口处放射性物质浓度的预测模型，在参考第二问的基础上，主要考虑风速和放射性物质扩散速度在空间中的矢量运算。在对上风口分析时，要分类讨论风速和自然扩散速度之间的大小关系，当风速小于自然扩散速度时，放射性物质是无法到达上风口的。对于问题四，应参考大量气象、地理、新闻资料，选择我国东海岸典型地 域-山东半岛，作为研究对象，综合考虑对应海域平均风速及风向、地理距离、 海水对放射性物质扩散的部分反射系数等因素，预测出放射性核物质与实际情 况比较。3 模型的假设考虑到放射性气体扩散的复杂性，为简单起见，在讨论扩散模型时都作了如 下假设；（1）瞬时泄漏假定瞬时完成，连续泄漏假定泄漏速率恒定；（2）气云在平整、无障碍物的地面上空扩散；（3）气云中不发生化学反应，地面对气云无吸收；（4）为水平风向，风速和风向不随时间变；5）气体的传播服从扩散定律，即单位时间通过单位发祥面积的流量与他的浓度梯度成正比;(6)气体的扩散看作空中末已连续店员向四周等强度瞬时释放气体 放射性气体在无穷空间适房的过程不发生性质变化。4 符号说明及名词解释4.1 符号说明符号说明s放射性气体的传播速度k风速，单位m / sH泄漏点O距有效地面的高度t任意扩散时刻C(x, y, z,t)空间任意一点的放射性物质浓度8 .(i 二 x, y, z)i空间任意一点的放射性物质的扩散系数的方差k (i = x, y, z) i空间任意一点的放射性物质的扩散系数Q空间域V空间域其体积S一规则的球面面积Q1在(t,t + A)内通过Q的流量Q2Q内放射性物质的增量符号说明Q 0从泄漏源泄漏的放射性物质的总量Ah附加高度Ts核泄漏出口处的温度T0环境温度a设地面反射系数Q源强，单位为kg / sb xQ yQ z分别为用浓度标准差表示的x, y, z轴上的扩散参数Vs沉降速度，单位为m / sWd地面干沉积率申冲洗系数T0.5放射性核素的半衰期4.2 名词解释烟羽又称烟云(smoke cloud)、烟流(smoke plume)：从烟囱中连续排放到 大气中的烟气流。由于烟羽各部分的运动速度不同，因而其外形也千变万化。不 同的烟羽形状表示污染物浓度的空间分布不同。它与大气湍流、大气稳定度、地 形地物、排放参数等有密切的关系。动力抬升：暖气流受锋面、辐合气流的作用被迫上抬，或者在运行中受地 形阻挡产生上升运动，这种空气在运动中由外力(不包括重力和浮力)使一部分 空气被抬上升。湍流扩散：是指湍流运动导致大气或水体中的污染物质或其他物质与周围 洁 流体的混合。5 模型的建立和求解5.1 问题一模型的建立与求解5.1.1 模型一的建立模型一的建立和求解 核电站源源不断泄漏引起气体扩散船舶可以看做在无穷空间有连续点源导致的 扩散过程，能够由二阶抛物型偏微分方程描述放射性气体扩散过程中浓度变化规 律。本问中由于不考虑风力的影响，且扩散出来的气体匀速向四周散开，这样经 过任意时刻t,扩散的气体围成一个半径为st的球，且距球心位置不同的地方浓 度值不同。讨论在无风的情况下，放射性气体以S m/s的匀速在大气中向四周扩散。核 电站源源不断泄漏引起气体扩散船舶可以看做在无穷空间有连续点源导致的扩 散过程，能够由二阶抛物型偏微分方程描述放射性气体扩散过程中浓度变化规 律。本问中由于不考虑风力的影响，且扩散出来的气体匀速向四周散开，这样经 过任意时刻t,扩散的气体围成一个半径为st的球，且距球心位置不同的地方浓 度值不同。利用MATLAB软件对数据进行线性拟合，采用微分方程模型得到核 电站周边放射性气体在不同地区，不同时间段的浓度变化。假设核电站泄漏点O距离有效地面的高度为H，以核泄漏点正下方的地面为坐标原点，X轴指向下风向，Y轴水平垂直于风向轴，Z轴为铅直方向，建立空间直角坐标系，则核泄漏点的空间坐标为 O（0，0，H）将气体从泄露时刻记作t=0,泄漏点选为坐标原点，时刻t无穷空间中任 点（x,y,z）的气体浓度记为C（x,y,z,t）假设单位时间内通过单位法向量面积的流量 与浓度梯度成正比，则有1)q = -k gradCi其中k i（i = x，y，z）是扩散系数严 表示浓度梯度，其中的负号代表放射性物质的浓 度是由高到低的地方扩散。假设空间域Q的体积为V，包围空间域Q的曲面为一规则的球面，设其表面 面积为S,外法线向量为n二（-，-2,1）z z ，则在（t,t + At）内流通过空间域Q的流量可表示为：Q = J +A 1 ff q nd dt1ts （2）空间域Q所包围的区域内放射性物质的增量可表示为： QC(x, y, z, t + At) - C(x, y, z, t)dV2V(3)而由泄漏源泄漏出的放射性物质的总量可表示为：Q =J t+A t BI0tQp dVdt04)根据质量守恒定律和连续性原理，单位时间内通过所选曲面S的向外扩散的放射性物质与S曲面内放射性物质增量之和，等于泄漏源在单位时间内向外泄漏的放射性物质。则有：Q = Q + Q0 1 25)即：fff C (x, y, z, t + At) C (x, y, z, t )dV + A+A JJ q -ndq dt = J t+At JJJ tsV6)p dVdtt0Q根据曲面积分的 Gauss 公式得：JJV7)则式子(6)可以转换成JJJC(xyzt+ At)-C(x，y，z，t) .AtdV + J+AJJJdivqdVdt = J t+aJJJ p dVdt Attt0VVQ(8)由于fAtAt竺-lim C(x，y,z,t + At) C(x，y,z,t) - iimTkd讽gradCCdt AtTOAtAtTO9）故式子（8）即可转换成p dV At010)即：B! aC dV + B! divqdV = pdfo11)VV根据 A.Fick 扩散微分方程式中：QC+ u QC+ u匹+ U QC :Q22-K+ K Q 2+ K Q 2atx Qxy Qyz Qzx Qx 2y Qy 2z Qz 2其中：C为气体浓度；t为时间；u ,u ,u为x,y,z方向风速；k k k为x,y,z x y zx y z方向上的扩散系数。由于无风时u二u =u = 0，为简化模型便于求解，假设扩散各向同性且扩 xyz散系数为常数，即K =K =K =K，则根据式（1），无风条件下扩散方程为 XYZQCd 2 K d 2K d 2=K+Qtx Qx2y Qy 2z Qz2结合式子（12），可解出式子（11）得到结果为：c（八p x2y2（z - H ）2 （4兀