气体扩散模型

放射气体模型的预估模型 摘要摘要 本文是以日本福岛核电站遭遇自然灾害发生核泄漏的背景而提出的 且结 合了高斯烟羽模型 线性拟合 以及微分方程模型 运用 MATLAB 软件 分 析泄漏源强度 风速 大气稳定度参数等因素对放射性气体扩散的影响 预测 了放射性气体浓度在不同时间 不同地区的浓度变化 并且本文模型中数据可 以根据不同的实际情况而加以改变 因而是本文的应用范围大大增加 可以适 用于具有较强的应用型 对于问题一 讨论在无风的情况下 放射性气体以 m s 的匀速在大气中向四s 周扩散 本问中由于不考虑风力的影响 且扩散出来的气体匀速向四周散开 这样经过任意时刻 t 扩散的气体围成一个半径为 st 的球 且距球心位置不同 的地方浓度值不同 采用列数列的表现方法 设定相同时间段 t 把条件进行整 理 并经过简单计算得出每段时间所预测得到的扩散距离 r 和浓度 C 利用 MATLAB 软件对数据进行线性拟合 采用微分方程模型得到核电站周边放射性 气体在不同地区 不同时间段的浓度变化 得出随着离泄漏源距离的延伸 最 后放射性物质的浓度越来越小 趋近于零 即当 x 趋向无穷时 C x y z t 趋向 于零 当时间趋向于无穷时 C x y z t 也趋于无穷 对于问题二 要探究风速对放射性物质浓度分布的影响 风速的处理是问题的 中心 采用大气污染的经典高斯扩散模型 实现了高斯烟团气体扩散模型的动态 预测 分析计算了气体扩散过程中的各关键参数 对于问题三 本文在问题二的基础上 结合考虑风速和放射性物质扩散速度在 空间中的矢量运算 将在上风和在下风不同情况下与传播速度 之间的比较的s 分析 利用高斯烟羽模型对核电站周边地区的浓度进行预测 然后 利用 MATLAB 软件 将相关数据代入程序 我们得到核电站周边地区的浓度分布的 等高曲线 对于问题四 本文参阅整理大量气象 地理 新闻资料 选择我国东海岸 典型地域 山东半岛作为研究对象 综合考虑对应海域平均风速及风向 地理 距离 海水对放射性物质扩散的部分反射系数等因素 集合核电站周边的浓度 等高线 可 关键词 关键词 放射性气体 扩散 浓度变化 高斯修正模型 预测 1 问题的提出 由于重大的突发性核泄漏紧急灾害事件具有爆发性 空间分布不连续性 对周边地形和气象条件的敏感性的特点 研究核事故所释放的物质的时空分布 需要高度精确的技术 但是在对于更好地保护环境有着极其重要的意义 在有 一座核电站遇自然灾害发生泄漏 浓度为的放射性气体以匀速排出 速度为 0 p 在无风的情况下 匀速在大气中向四周扩散 速度为 m kg ss m s 问题一 若能建立一个描述核电站周边不同距离地区 不同时段放射性物 质浓度的预测模型 这对于研究核污染模式具有重要的意义 问题二 当风速为时 给出核电站周边放射性物质浓度的变化情况 k m s 这对于研究核电站附近浓度的在实际环境下有着重要的作用 问题三 当风速为时 计算出上风和下风公里处的放射性物质浓k m sL 度的预测模型就显得更加急迫的了 问题四 将建立的模型应用于福岛核电站的泄漏 计算出福岛核电站的泄 漏对我国东海岸的影响 这一实际的有意义预测 可以明确我们实际的污染情 况 为我们的核应急决策提供技术支持 2 问题的分析 对于问题一 在无风的情况下 放射性气体 以 m s 的速度 匀速在大气ss 中向四周扩散 在此条件下 探求一个模型来对核电站周边不同距离地区 不 同时段放射性物质浓度进行预测 我们要明确此问题研究的核扩散是点源连续 泄露的扩散问题 虽然只是要求考虑在无风情况下放射性物质浓度分布 但为 了使模型更贴切实际 需考虑地面反射 泄漏源有效高度等因素对浓度分布的 影响 根据 泄露放射性物质质量守恒定律 和 气体泄漏连续性原理 进行分 析 发现要得出核电站周边不同距离地区 不同时段放射性物质浓度的预测模 型 最后对于该方程进行分析求解 对于问题二 为了探究风速对发生核泄漏的核电站周边放射性物质浓度分 布的影响 运用概率学知识 通过图解和数学推导得出 连续点源放射性物质 高斯扩散模型 应在 连续点源放射性物质高斯扩散模型 的基础上经多次 合理修正后得到更好的 优化高斯模型 对于问题三 该问题要求建立泄漏源上风口和下风口处放射性物质浓度的 预测模型 在参考第二问的基础上 主要考虑风速和放射性物质扩散速度在空 间中的矢量运算 在对上风口分析时 要分类讨论风速和自然扩散速度之间的 大小关系 当风速小于自然扩散速度时 放射性物质是无法到达上风口的 对于问题四 应参考大量气象 地理 新闻资料 选择我国东海岸典型地 域 山东半岛 作为研究对象 综合考虑对应海域平均风速及风向 地理距离 海水对放射性物质扩散的部分反射系数等因素 预测出放射性核物质与实际情 况比较 3 模型的假设 考虑到放射性气体扩散的复杂性 为简单起见 在讨论扩散模型时都作了如 下假设 1 瞬时泄漏假定瞬时完成 连续泄漏假定泄漏速率恒定 2 气云在平整 无障碍物的地面上空扩散 3 气云中不发生化学反应 地面对气云无吸收 4 为水平风向 风速和风向不随时间变 5 气体的传播服从扩散定律 即单位时间通过单位发祥面积的流量与他的浓 度梯度成正比 6 气体的扩散看作空中末已连续店员向四周等强度瞬时释放气体 放射性气体在无穷空间适房的过程不发生性质变化 4 符号说明及名词解释 4 1 符号说明 符号符号说明说明 s 放射性气体的传播速度 k风速 单位 m s H泄漏点距有效地面的高度O t 任意扩散时刻 C x y z t 空间任意一点的放射性物质浓度 i ix y z 空间任意一点的放射性物质的扩散系数的方差 i ix y z 空间任意一点的放射性物质的扩散系数 空间域 V空间域其体积 S一规则的球面面积 1 Q 在内通过的流量 t tt 2 Q 内放射性物质的增量 符号符号说明说明 0 Q 从泄漏源泄漏的放射性物质的总量 h 附加高度 s T 核泄漏出口处的温度 0 T 环境温度 设地面反射系数 Q 源强 单位为 kg s xyz 分别为用浓度标准差表示的轴上的扩散参数 x y z Vs 沉降速度 单位为 m s d W 地面干沉积率 冲洗系数 0 5 T 放射性核素的半衰期 4 2 名词解释 烟羽烟羽又称烟云 smoke cloud 烟流 smoke plume 从烟囱中连续排放到 大气中的烟气流 由于烟羽各部分的运动速度不同 因而其外形也千变万化 不同的烟羽形状表示污染物浓度的空间分布不同 它与大气湍流 大气稳定度 地形地物 排放参数等有密切的关系 动力抬升动力抬升 暖气流受锋面 辐合气流的作用被迫上抬 或者在运行中受地 形阻挡产生上升运动 这种空气在运动中由外力 不包括重力和浮力 使一部 分空气被抬上升 湍流扩散 湍流扩散 是指湍流运动导致大气或水体中的污染物质或其他物质与周围 洁 流体的混合 5 模型的建立和求解 5 1 问题一模型的建立与求解 5 1 1 模型一的建立 模型一的建立和求解模型一的建立和求解 核电站源源不断泄漏引起气体扩散船舶可以看做在无穷空间有连续点源导致的 扩散过程 能够由二阶抛物型偏微分方程描述放射性气体扩散过程中浓度变化 规律 本问中由于不考虑风力的影响 且扩散出来的气体匀速向四周散开 这 样经过任意时刻 t 扩散的气体围成一个半径为 st 的球 且距球心位置不同的 地方浓度值不同 讨论在无风的情况下 放射性气体以 m s 的匀速在大气中向四周扩散 核s 电站源源不断泄漏引起气体扩散船舶可以看做在无穷空间有连续点源导致的扩 散过程 能够由二阶抛物型偏微分方程描述放射性气体扩散过程中浓度变化规 律 本问中由于不考虑风力的影响 且扩散出来的气体匀速向四周散开 这样 经过任意时刻 t 扩散的气体围成一个半径为 st 的球 且距球心位置不同的地 方浓度值不同 利用 MATLAB 软件对数据进行线性拟合 采用微分方程模型 得到核电站周边放射性气体在不同地区 不同时间段的浓度变化 假设核电站泄漏点 O 距离有效地面的高度为 H 以核泄漏点正下方的地面 为坐标原点 X 轴指向下风向 Y 轴水平垂直于风向轴 Z 轴为铅直方向 建 立空间直角坐标系 则核泄漏点的空间坐标为 O 0 0 H 图 1 空间坐标系示意图 将气体从泄露时刻记作 t 0 泄漏点选为坐标原点 时刻 t 无穷空间中任一 点 x y z 的气体浓度记为 C x y z t 假设单位时间内通过单位法向量面积的流 量与浓度梯度成正比 则有 1 i qgradC 其中 是扩散系数 表示浓度梯度 其中的负号代表放射性物质的 i ix y zgrad 浓度是由高到低的地方扩散 假设空间域的体积为 V 包围空间域的曲面为一规则的球面 设其表面 面积为 S 外法线向量为 1 xy n zz 则在内流通过空间域的流量可表示为 t tt 2 1 tt t s Qq nddt 空间域所包围的区域内放射性物质的增量可表示为 3 2 V QC x y z ttC x y z t dV 而由泄漏源泄漏出的放射性物质的总量可表示为 4 00 tt t Qp dVdt 根据质量守恒定律和连续性原理 单位时间内通过所选曲面的向外扩散S 的放射性物质与曲面内放射性物质增量之和 等于泄漏源在单位时间内向外S 泄漏的放射性物质 则有 5 012 QQQ 即 0 tt tt t t Vs C x y z ttC x y z t dVq nd dtp dVdt 6 根据曲面积分的 Gauss 公式得 7 sV q nddivqdV 则式子 6 可以转换成 0 tt tt t t VV C x y z ttC x y z t t dVdivqdVdtp dVdt t 8 由于 00 limlim tt t tt kdiv gradC dt CC x y z ttC x y z t ttt 9 故式子 8 即可转换成 10 0 VV C dVtdivqdVtp dVt t 即 11 0 VV C dVdivqdVp t 根据 A Fick 扩散微分方程式中 222 222x xyzyz CCCCVVV txyzxyz uuu 其中 C 为气体浓度 t 为时间 zyx uuu 为 x y z 方向风速 为 x y z x y z 方向上的扩散系数 由于无风时0 zyx uuu 为简化模型便于求解 假设扩散各向同性且扩 散系数为常数 即 则根据式 1 无风条件下扩散方程为 XYZ 222 222x yz CVVV txyz 结合式子 12 可解出式子 11 得到结果为 222 0 1 50 5 exp 4 444 xyzxyz pxyzH C x y z t tttt 13 这结果表明 对于任意时刻 t 烟雾浓度 C 的等值是球面 并且 2222 Rzyx 随着球面半径 R 的增加 C 的只是连续减少的 当 R或 t时 C x y z t 0 此模型仅能预测在无边界空间环境中 且不考虑放射性物质释放出时的初 动量的理想化情况下不同时间里放射性物质的浓度变化情况 为使模型更贴切 实际 下面分别通过考虑有效泄漏源及地方反射这两方面因素对模型进行优化 考虑有效泄漏源对模型一进行优化考虑有效泄漏源对模型一进行优化 图 2 如图2所示 核泄漏口的有效高度为 h 放射性物质从泄漏口排出时由 于受到热力抬升和本身动力抬升 产生的附加高度为 h 因此有效泄漏源 的高度 H h h 放射性物质受热力抬升的高度主要由其排放时本身的初始动量和因其温 度高于环境温度产生的静浮力决定 同时还会受到温度和风速等因素的影 响 根据烟气抬升高度的综合分析公式 0 40 6 0 925 25 sb s V DFh h u 14 其中 h 泄漏源的实际高度 泄漏源出口处的风速 已知为 m s s uk 泄漏源出口的有效直径 D 放射性气体的扩散速度 已知为 m s s Vs 浮力通量 b F 4 ms 根据 Briggs 公式 浮力通量 2 2 sa bs s TTD FV g T 15 其中 核泄漏出口处的温度 环境温度 s TK 0 TK a 在有风 且释放气体温度与环境温度差 35K 时 1 0 s um s 35 sa TTK 抬升高度 16 0 40 6 0 920 792 s s V DQh h u 其中 为大气压强 2 0 0 0 275 s s TT QPD V T PhPa b 小风 时 且温度差 35K 抬升高度为 1 0 s um s 35 sa TTK 17 0 375 0 25 0 5 50 0098 dT hQ dZ 其中 为泄漏源的有效高度处上的环境温度梯度 0 dT dZ K m c 当温度差 35K 时 抬升高度 0 35 s TTK 18 2 1 50 01 s s V DQ h u 其中 当风速时 取 1 0 s um s s u1 0 m s 由本题题目知风速为 k 放射性物质的扩散速度为 s 分别代入式子 16 17 和 18 可得抬升高度 0 0 40 6 0 375 0 25 0 0 0 0 920 792 1 0 35 5 50 0098 1 0 35 2 1 50 01 1 0 1 035 s s s sDQh km sTTK k dT Qkm sTTK dZ sDQ kkkTTK k h 适用于且 适用于且 适用于不限 时 取 其中 H 为泄漏源的有效高度 K 为泄漏源出口处的风速 D 为泄漏源出口处的有效直径 S 为放射性物质的扩散速度 Ta 为泄漏源出口处的温度 K T0 为环境温度 K 2 0 0 0 275 s TT QPD VP s T 为标准大气压 综上所述 泄漏源的有效高度为 Hhh 考虑地面反射对模型一进行优化考虑地面反射对模型一进行优化 地面会对排放出来的放射性物质有一定的反射作用 而放射性物质因受 沉降和重力等因素的影响并不能全部被地面反射回去 因此放射性物质只能 部分返回空气中 由此可知空气中的放射性物质由两部分组成 其中包括由 泄漏口直接排向空气中的和由地面反射回到空气中的 假设地面反射系数为 根据几何物理学知识可得 通过地面反射进入大气的核扩散物质相当于虚 拟泄漏源泄露的核辐射物质在原来空间的一个浓度叠加 如果所设核泄漏点 O 源在距有效地面的高度为的地方 本文求空间任一点的浓度值 OH A x y z 则实际泄漏源对点的影响部分可用来表示 在考虑反射系数后 A 2 exp 2 z zH 虚泄漏源对点的影响部分可用来表示 因此式子 13 所 OA 2 exp 2 z zH 得模型可以优化为 222 0 1 50 52 expexp 4 4442 xyzxyzz pxyzHzH C x y z t tttt 其中 泄漏源的有效抬升高度为 Hhh 5 2 问题二模型的建立与求解问题二模型的建立与求解 5 2 1 模型二的建立 问题二 针对在有风速时 由于地面风速对大气扩散的影响起着至 k m s 关重要的作用 不同的风速其浓度计算方法有较大差异 因此 下面的分析从 有风 大于 1 5 m s 小风 小于1 5m s大于0 5m s 和静风 小于0 5m s 三种情况分别论述 在此模型中 我们同样采用模型一中的空间坐标系 5 2 1 1 无风时 当风速k 1m s时 可以看作是无风时的情况 即此时采用模型一的瞬时点源扩 散模型求解 即此时的浓度场有 222 0 1 50 5 exp 4 444 xyzxyz pxyzH C x y z t tttt 14 再利用部分反射倾斜烟云模式 可得出该情况下的浓度分布 2222 0 3 22222 0 exp exp exp exp 2 2222 T ss xyzxyzz pzHV TzHV TxuTy c x y z HdT 15 上式中假定 T代表积分时间 即静风持续时间 通常按照每小 xy 时排放6个烟团进行积分 即可达到需要的精度 5 2 1 2 有风时 即 当风速大于 1 0m s 时 选取高斯连续点源扩散模型进行分析 这也是在所 有有风速时 放射性物质扩散模型中最常见也最方便的一种泄漏方式 因此时的风速 u 不变 y z 方向上的速度为 0 则此时的扩散微分方程为 22 0 22 yz ccc up xyz 16 其中 模拟连续点源烟羽扩散的定解条件是 当 0c x y z yx ucd dQ 则解是 17 22 0 0 5 exp 4 44 yzyz pyz C x y z t tutt 以上是一三维高斯函数式 表示连续排放时的烟团浓度呈以烟团中心为原 点的正态分布 以表示分布的方差 他们和扩散系数的关系是 222 y xz 18 2 2 xxt 2 2 y yt 2 2 zzt 由于高斯连续点源扩散模型有边界点源 这里设泄漏源有效高度为 取H 其地面投影为坐标原点 轴指向风向 继续考虑地面反射作用 x 图4 边界点源示意图 则此时的 17 式可以改成为 222 0 222 exp exp exp 2222 yzyzz pyzHzH c x y z H k 19 其中 为源强 为污染物质量浓度 分别为Q kg sc 3 kg m xyz 用浓度标准差表示的轴上的扩散参数 为泄漏高度的平均风速 m s x y zk 为泄漏有效高度 Hm 考虑到动力抬升的作用 将有效高度 即有效烟云高度 HeHH 图5 动力抬升示意图 因此 19 式可以写成为 222 0 222 exp exp exp 2222 yzyzz pyzHezHe c x y z H k 20 式 20 是高架连续点源的一般解析式 由 20 式可以推导出各种条件 下的常用大气扩散模型 1 泄露源为近地面 当为地面连续点源时 即He 0时 由 20 得 21 22 0 22 0 exp exp 22 yzyz pyz c x y z k 地面连续点源在地面上任一点产生的浓度为 22 2 0 2 0 0 exp 2 yzy py c x y k 2 泄漏源为高架点时 当为高架连续点源时 地面任意一点浓度模型 此时令 20中的 z 0 即可得到高架连续点源的地面浓度模型 23 22 0 22 0 exp 22 yzyz pyHe c x yHe k 5 2 2模型的计算求解 5 2 2 1 模型参数的确定 为了求解出核电站周边放射性物质浓度的具体的变化情况 则应首先把大 气条件分类 从而确定出每类天气条件下的扩散参数的值 在事故过程 yz 模型中所需参数的选取及确定十分重要 通常情况下气象参数的选取是利用该 地区多年气象资料 采取工业安全与环保统计的方法进行有关参数的确定 而 其他扩散参数是以实际测定为基础的 a 大气稳定度的计算大气稳定度的计算 对于大气条件的分类 根据国家标准 GB T 13201 1991 制定地方大气 污染物排放标准的技术方法的规定 划分大气稳定度的级别 共分为6级A F A为极不稳定 F为极稳定 其具体分类情况如下 表 1 大气稳定度的级别参考表 白天太阳辐射有云的夜晚地面风速 1 m s 强中弱 阴天的白天或夜 间薄云遮天或低云 0 5 云量 0 4 6CDDDDD b 扩散参数的计算 扩散参数的计算 有风时的扩散参数 的确定采用Briggs给出一套扩散参数 则每一 yz 个大气稳定度所对应的扩散参数如表1和表2所示 表2 Briggs 扩散参数 开阔平原田野 大气稳定度 y z A 2 1 0001 01 22 0 xx x2 0 B 2 1 0001 01 16 0 xx x12 0 C 2 1 0001 01 11 0 xx 2 1 0002 0 1 08 0 xx D 2 1 0001 01 08 0 xx 2 1 0015 0 1 06 0 xx E 2 1 0001 01 06 0 xx 0003 0 1 03 0 xx F 2 1 0001 01 04 0 xx 0003 01 016 0 xx 表3 Briggs 扩散参数 城市 大气稳定度 x y A B 2 1 0004 0 1 32 0 xx 2 1 0001 0 1 24 0 xx C 2 1 0004 0 1 22 0 xx x20 0 D 2 1 0004 0 1 16 0 xx 2 1 0003 0 1 14 0 xx E F 2 1 004 0 1 11 0 xx 2 1 0015 0 1 08 0 xx 则欲求解出核电站周边放射性物质浓度的变化情况 只需将各个大气稳定 度所对应的带入计算即可 xy 5 2 2 2 模型的求解思路整合 图6 模型二计算流程图 5 2 3 对高斯浓度计算模型的修正 利用标准的高斯模型 计算大气中放射性核素云团的扩散 目前已有较多 的研究 而实际的核素扩散过程还存在着粒子的重力沉降 雨洗沉积以及核素 衰变等的因素对浓度分布的影响 因此 有必要对高斯浓度扩散模式进行修正 方可较为准确 真实地反映实际的核素扩散规律 a 考虑干沉积时的连续点源扩散 考虑干沉积时的连续点源扩散 粒径大于10 m 的粒子有明显的重力沉降 粒子的沉降速度取决于空气阻力 和重力平衡 可用斯托克斯公式表示 24 2 18 gD Vs 式中 粒子密度 3 kg m 重力加速度 9 806 5 g m s 粒子直径 Dm 空气的动力粘性系数 可取 5 1 8 10 kgm s 沉降速度 含碘放射性核素的干沉积速度为 1 1 cm s 该Vs m sVs 值也可由式 计算得到 因为在扩散过程中同时有重力沉降的位移迭加到羽流中心线上 中心线就 会向下倾斜 所有粒子相当于在下倾的中心线上扩散 该类扩散和沉降的迭加 可认为是羽流运行过程中 实源以的速度向下移动 在 处向下移动的高度 s Vx 为 即源高由降到了 实际上 由于大气湍动及其他动力作 s s V x V t u H s V x H u 用 地面不是全吸收表面 应考虑地面的反射作用 但毕竟存在粒子的沉降作 用 又不是全反射 因此 需对反射项乘以一个反射系数 由于反射 1 项的有效源高度也变成了 故相应的浓度计算公式为 s V x H u 25 22 2 222 exp exp exp 2222 ss yzyzz V xV x zHzH Qy uu c x y z H k 式 20 可得到经过干沉积作用后的浓度分布结果 式 20 中 反射系数需外 部给定 通常对于放射性核素可取0 5 1 b 地面的干沉积量地面的干沉积量 根据扩散理论和动量传递的普朗特理论 可以导出干沉积的地面沉积量为 27 0 ds WVc x y 式中 前面各种情形计算的地面污染物的浓度 0 c x y 干沉积速度 s V m s 地面干沉积率 d W 21 Bq ms c 雨洗作用雨洗作用 降雨对烟羽中的颗粒物及气溶胶具有清洗作用 可溶性气体与蒸汽亦可溶 于雨水中 降雨过程造成的这类湿沉积是导致放射性气溶胶和气体向地面沉积 的另一重要机制 通常以冲洗系数 描述降雨对烟羽中污染物清洗作用 1 S 的大小 与雨强的关系可以表达为 28 b aI 式中 为雨强 mm h I a b为经验系数 式中 按释放物质为含碘 不含碘情况分别取值 对含碘物质 取 对于不含碘物质 取 5 8 10 0 6ab 5 1 2 10 0 5ab 对于湿沉积导致的烟羽耗减 可采用湿沉积耗减因子对源强进行修正 Q 有 29 exp x Q xQ u d 放射性衰变的影响放射性衰变的影响 除了干沉积和湿沉积外 放射性物质的衰变也是影响大气中核素浓度分布的主 要原因 由于放射性物质服从简单的衰变规律 其浓度随时间的变化可由下式 计算 30 0 t cc e 式中 为初始浓度 0 c 衰变常数 经过的时间 t 可由无衰变的浓度公式计算 0 c 参照前面的模式 同样采用衰变耗减因子来对源项进行修正 经过理论推导 可以得出 31 0 5 0 693 exp 3600 x Q xQ T u 其中 放射性核素的半衰期 0 5 T 如以上情况都存在 应同时考虑源项修正 即式 24 和式 26 相乘后代替 5 3问题三模型的建立与求解问题三模型的建立与求解Q 5 3 1模型三的建立和求解模型三的建立和求解 在问题二的基础上 结合风速和放射性物质的扩散速度在空间中的矢量运算 将在上 风和下风情况下的风速与放射性物质的传播速度 s 之间的比较分析 知道必须在满足风速 大于放射性物质的传播速度的情况下 放射性物质才能到达上风口 上风情况下 满足条件k m s s m s 将模型二中的求解结果中的风速 k 用 k s 代替 即可求出在上风情况下放射性物质浓度的变化 222 222 exp exp exp 2 222 yzyzz QyzHezHe c x y z H ks 31 下风情况下 风速的大小对放射性物质能否到达风口均无影响 将模 型二中的风速 k 用 k s 代替 即可求出在下风情况下放射性物质浓度的变化 222 222 exp exp exp 2 222 yzyzz QyzHezHe c x y z H ks 32 模型四的建立和求解模型四的建立和求解 在结合模型二 模型三的条件下 在参阅整理大量气象 地理 新闻资料 选择我国东海岸 山东半岛作为研究对象 综合考虑对应海域平均风速及风 向 地理距离 海水对放射性物质扩散的部分反射系数等因素 并通过 C 编 程模拟计算 预测出放射性核物质将经过 6 5 天到达我国东海岸 山东半岛 且 131I 浓度值为 0 100 与 与实际情况比较吻合 3 mBq m 6 结果分析与检验 由于在前四问题中均无具体的数据 所以对于结果 在检验上存在一些麻 烦 但是在模型三预测上下风的浓度具有较好的实际意义 而在模型四中 对 于预测出浓度的变化和对我国的影响也较为合理 故该模型具有可靠性和参考 性 可以进行推广和其它方面的应用 7 模型的优缺点与改进方向 优点 优点 该方法可计算核事故中连续点源和瞬时点源在不同气象 地形条件 下的浓度分布 该结果在核事故的应急救援过程中 对救援人员划定警戒区和 确定周围居民的疏散范围具有重要意义 并可为制定救援方案和应急决策提供 科学依据 放射性云团在空中迁移和扩散提供浓度的定量描述 对可能发生的 核事件的放射性核素浓度监测及监测时间范围提供相关信息 缺点 缺点 模型在建立时只考虑了两种下垫面情况 即农田和城市 实际的扩 散情况可能因地形条件的复杂化而有所改变 由于扩散是一个非常复杂 影响 因素众多的过程 上述方法还有一定的局限性 有待于进一步检验 另外 该 方法只给出了浓度的计算结果 没有换算成人的吸收剂量 今后将对其计算过 程和参数的选取进行多方面验证 使其进一步完善 参考文献 1 金勇进 统计学 北京 中国人民出版社 2010 2 梅长林 周家良 实用统计方法 北京 科学出版社 2002 3 罗艾民 魏利军 有毒重气泄露安全距离数值方法 J 中国安全科学学报 2005 15 8 98 100 4 周燕等 Matlab 在统计与工程数据分析中的应用 北京 电子工业出版社 2010 5 丁信伟 王淑兰 徐国庆 可燃及毒性气体泄漏扩散研究综述 化学工业与工程 1999 2 16 118 122 6 胡世明 气体释放源的三维瞬态重气扩散研究 劳动保护科学技术 2O02 3 20 28 30 7 蒋军成 潘旭海 描述重气扩散的一种新模型 南京工业大学学报 2OO2 I 24 41 46 附 录 附录一 模型二在上下风附录一 模型二在上下风 L 公里的浓度求解源代码公里的浓度求解源代码 include include include using namespace std define e 2 718281828459 define PI 3 1415926 void qiuxy char wending double x double break case B dy 0 16 x sqrt 1 0 0001 x dz 0 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