基于高斯烟羽模型的放射性气体的扩散
放射性核素大气扩散程序SDUG的开发与初步验证
放射性核素大气扩散程序SDUG的开发与初步验证曹博;杨晔;陈义学【期刊名称】《原子能科学技术》【年(卷),期】2013(000)0z2【摘要】放射性核素大气扩散过程模拟是进行核事故后果评价和应急决策的主要内容之一。
高斯烟羽模式计算方法简单便捷，在核事故放射性核素释放的后果评价中获得了广泛应用。
本工作在高斯烟羽模式的基础上，考虑了核素粒子的重力沉降、雨洗作用以及放射性衰变等因素的影响，开发了放射性核素大气扩散程序SDUG。
该程序可计算正常工况和事故工况下，连续点源和瞬时点源在不同气象和地形条件下的放射性核素浓度分布。
通过与参考程序的数值结果比较发现，SDUG程序计算结果符合较好，能满足早期核应急和后果评价的需要。
【总页数】5页(P472-476)【作者】曹博;杨晔;陈义学【作者单位】华北电力大学核科学与工程学院，北京 102206;华北电力大学核科学与工程学院，北京 102206;华北电力大学核科学与工程学院，北京 102206【正文语种】中文【中图分类】TL732【相关文献】1.三维氢气爆炸分析程序的开发及初步验证 [J], 黄涛;高颖贤;丁书华;钟明君;吴丹;苏光辉;秋穗正2.GASFLOW程序液膜模型开发及初步验证 [J], 王方年;沈峰;程旭;黄兴冠3.先进中子学栅格程序KYLIN-2的开发与初步验证 [J], 柴晓明;吴文斌;涂晓兰;郭凤晨;尹强;黄世恩;芦韡;卢宗健;姚栋;李庆4.MOC/SN耦合三维中子输运程序rKYCORE开发与初步验证 [J], 唐霄;李庆;柴晓明;涂晓兰;王侃5.基于MCNP5和ORIGEN2.1程序的停堆剂量程序开发与初步验证 [J], 马辉强;苑旭东;于涛;陈珍平;谢金森;吴菱艳因版权原因，仅展示原文概要，查看原文内容请购买。
一种改进的放射性气体扩散高斯预估模型算法
（ ｃｏｌｆＴａ ｃ＆ Ｔａｓｏｔｔ ｎ Ｌ ｎｈ ｕＪａｔｎ 西 ｍ ，Ｌ ｎｈｕ７ ０ ７ ，Ｃ ｉａ Ｓｈｏ ｒｆ ｏ ｆ ｉ ｒｎｐｒ ｉ ， ａｚｏ ｉｏ ｇ ａｏ ｏ ｅ ａｚｏ ３ ０ ０ ｈｎ ）
Ａ ｂｓｒ ｃ ｔ ａ ｔ： Ｂｙ ａ ｌ ｚｎ ｈｅ ｃ ａ ａ ｔ ｒｓｉｓ ｏ ａ ｉ ａ ｔ ｅ ｇ ｓｄｉ ｕ ｉ ｎ ｂ ｓ ｄ ｏ ｕｓｉ ｎ ｍｏ ｅ ’ ｉ ｉ ｉｎ，ｈｉ ａ ｒｉ ｔｏ ｎａｙ ｉ ｇ ｔ ｈ ｒ ｃｅ ｔ ｆｒ ｄ ｏ ｃｉ ａ ｆ ｓｏ ａ ｅ ｎ Ｇａ ｓａ ｄ ｌ Ｓｌｍ ｔ ｏ ｔ ｓｐ ｐｅ ｎｒ － ｉ ｃ ｖ ｆ ｍ ｄｕ ｅ ｈ ｎ ｍｉ ａ ｅ ｏ ｅ ｋａ ｅ，ａ ｄ ｓ ｎ ｈ ｓｚ ｄ ｔ ｅ ｗｉｄ ｓ ｅ ｄ ａ ｄ ｉｆ ｉｎ ｖ ｌ ｃｔ ｆｒ ｄ ｏ ｃｉｅ ｇ ｓｕ ｉｇ ａ ｇｅ’ ｃ ｄ ｔ ｅ ｄｙ ａ ｃ ｒｔ ｆｌ ａ ｇ ｎ ｙ ｔ ｅ ｉｅ ｈ ｎ ｐ ｅ ｎ ｄｆｕｓｏ ｅ ｏ ｉ ｏ ａ ｉ ａ ｔｖ ａ ｓｎ ｎ ｌ Ｓ印 － ｙ
ｒ ｄｉａ ｔ ａ ｅｅ ｓ ｓ ａ ｏ ｃｉ ｇ ｓ ｒｌ ａ ｅ ． ｖｅ
Ｋｅ ｒ ｓ ｙ ｗｏ ｄ ：Ｇａ ｓｉｎ ｍｏ ｅ ； ｒ ｄｏ ｃ ｉｅ ｇ ｓ ｄ ｆ ｓｏ ｕ ｓ ｄ ｌ ａ ｉａ ｔ ａ ； ｉ ｕ ｉ ｎ；ｒ ｔ ｆｌａ ａ ｅ；ｓａｅ ｏ ｉ ｏ ａ ｖ ｆ ａｅ ｏ ｋ ｇ ｅ ｔ ｔ ｆａｒ ｗ； ｄｆ ｓｏ ｅｏ ｉ ｌ ｆ ｉ ｕ ｉｎ ｖ ｌｃｔ ｆ ｙ
ｄ ｉ１ ． ９ ９ ｊｉ ｎ １０ —６ ５ ２ １ ． １０ ４ ｏ：０ ３ ６ ／ ． ｓ ．０ １３ ９ ．０ ２ ０ ．３ ｓ
放射性气体扩散的预估模型(参考版)
核能作为可持续能源，在高新科技领域地位愈发重要。
但对核能技术开发的规范却一直没能跟上。
2011 年3 月12 日，日本福岛地震引起的福岛核电站泄漏，最终引致放射性物质在日本海域扩散，并有传言说中国东海及美国北大西洋海域亦将受到影响，为了检验此种说法，环境安全，核技术开发等领域的学者都结合自身知识设计出相应的预测模型或模拟软件。
本队通过查找文献资料，借鉴已有的标准高斯模型，对模型做出变化风速，干沉积率，湿沉积率以及放射性衰变修正以更好地还原福岛核电站的泄漏事件的具体环境。
最后，本队针对模型的适用性做出了适当推广，并对不足部分做出了解释。
关键词：标准高斯模型变化风速干沉积率湿沉积率放射性衰变放射性气体扩散的预估模型w放射性气体扩散的预估模型4放射性气体扩散的预估模型一、问题提出福岛核电站位于北纬37 度25 分14 秒，东经141 度2 分，地处日本福岛工业区。
它是目前世界最大的核电站，由福岛一站、福岛二站组成，共10 台机。
组历史上福岛核电站曾发生过多起事故：1978 年，福岛第一核电站曾经发生临界事故，但是事故一直被隐瞒至2007 年才公之于众；2005 年8 月，里氏7.2 级地震导致福岛县两座核电站中存储核废料的池子中部分池水外溢；2006 年，福岛第一核电站6 号机组曾发生放射性物质泄漏事故。
最近的一起事故发生在今年3 月份，里氏9.0 级地震导致福岛县两座核电站反应堆发生故障，其中第一核电站中一座反应堆震后发生异常导致核蒸汽泄漏，辐射物质进入空气，并向周边扩散。
本次事故引起了日本周边国家的高度重视，部分邻国国民出现了对放射性物质扩散到本国导致人员伤亡而感到恐慌。
为此，本文将就日本福岛核电站爆炸事故所产生的放射性物质通过空气扩散的实际情况进行建模分析，探究其扩散过程和规律，并对其可能对邻国造成的影响进行定量和定性的评价。
首先，我们对实际问题进行抽象，将其简化为以下几个小问题，逐一解决：设有一座核电站遇自然灾害发生泄漏，浓度为p0 的放射性气体以匀速排出，速度为m kg/s，在无风的情况下，匀速在大气中向四周扩散, 速度为s m/s. 1. 建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。
基于一种改进的高斯模型研究危险气体短时泄漏后扩散规律
２ 模 型扩散标 准差 的计算
表１ 用Ｋ ｌ ｕ ｇ方法计算大气扩散标准差时所需参数
注： ｒ ｏ １ ３ １ ＂ ｙ 叮 分另 ｑ 为ｘ ． ｙ 、 ｚ 方 向扩散标准差。且 ｙ （ ｘ ）＝Ｒ 。 ｘ 、 Ｏ ＂ ｘ （ ｘ ） ＝ｏ ｒ ｙ （ ｘ ） 和ｏ ｒ ｔ （ ｘ ）＝Ｒ ｘ 。
ｃ— ｅ 卜
操作失误或其他原因可能导致上述物质发生短时间 泄漏 ， 并对周围居民和环境造成极大 的危害¨ 。日
常生 活 中 ， 由此类 易燃 易爆 物 质 泄漏 造 成 惨 重 损 失
的案例在 国内外屡见不鲜 Ｊ 。因此 ， 为了避免此类 事故的发生 ， 有必要进行相应 的风险评估工作 。
和理论依据 。
关键词 ： 改进 ； 高斯扩散模型 ； 短时间排放
中图分类号 ： ｘ ５ １ １
在天 然气 、 液 化石油 气 等石 油化 工行 业 中 ， 由于
为有效源高（ 排放高度与抬升高度之和） ， ｍ； 、 Ｙ 、 ＝ 为 空间点各方向的坐标 ， ｍ； 其 中， 为气体在时间 ｔ 内 的向下游传输 的距离 ， 标准差 是关于 的函数。 假 定 短 时排放 源 释放 的烟片是 由无数 个 首尾 相 接 的小 烟 团叠 加后 的结 果 ， 则 一段 时 间 内泄 漏 源 释 放气 体 的质量 为 ： ｄ Ｍ＝ ｑ （ ｔ ‘ ） ‘ （ ０ ≤￡ ‘ ≤ ｔ ， ） （ ２ ） 式中： ｔ 为任一烟团释放的时间， ｓ ； ｔ ｒ 为烟团释放的 总时间 ， ｓ ； ｇ （ ￡ ） 为其释放率， ｇ ／ ｓ 。 将（ ２ ） 式代人（ １ ） 式并积分可得某一时刻 ， 污染 物在 扩散 空 间一点 的浓 度为 ：
基于高斯修正模型的放射性气体扩散浓度预测(精)
基于高斯修正模型的放射性气体扩散浓度预测摘要本文是以日本福岛核电站遭遇自然灾害发生核泄漏的背景而提出的。
对于问题一，考虑到放射性物质的泄漏是连续不断的。
本文根据“泄漏放射性物质质量守恒定律”和“气体泄漏连续性原理”建立了微积分方程，应用了高等数学中散度、梯度、流量等数学概念，通过Guass公式、四维二阶偏微分方程，因而得到了核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。
同时为使模型适用范围更广，本文引入了地面反射系数，考虑了由于放射性物质从泄漏口喷出时具有初动量而使其泄漏源有效高度被抬高等因素，进而得到了在无风环境中适用范围更广的“高斯修正模型”。
对于问题二，要探究风速对放射性物质浓度分布的影响。
本文运用概率学[1]知识，通过图解和数学推导得出“连续点源放射性物质高斯扩散模型”。
本文依次考虑了“重力沉积”、“雨水沉积”、“核衰变”等因素对浓度分布的影响。
并通过构建“耗减因子”、“衰变因子”等方法将耗减和衰变的放射性物质“投影”到泄漏源浓度中，得到了经多次合理修正后的“优化高斯模型”，并据此分析了泄漏源周边地区放射性物质的浓度变化。
针对于问题三，本文在问题二的基础上，结合考虑风速和放射性物质扩散速度在空间中的矢量运算。
得出在对上风口分析时，要分类讨论风速和自然扩散速度之间的大小关系，当风速小于自然扩散速度时，放射性物质是无法到达上风口的。
对于问题四，本文参阅整理大量气象、地理、新闻资料，选择我国东海岸典型地域---山东半岛作为研究对象，综合考虑对应海域平均风速及风向、地理距离、海水对放射性物质扩散的部分反射系数等因素，并通过C++编程模拟计算，预测出放射性核物质将经过6.5天到达我国东海岸，且0.100mBq⋅m-3与实际情况比较吻合。
131I浓度预测值为：关键词：放射性气体扩散浓度变化高斯修正模型预测1 问题的提出由于重大的突发性核泄漏紧急灾害事件具有爆发性、空间分布不连续性、对周边地形和气象条件的敏感性的特点，研究核事故所释放的物质的时空分布需要高度精确的技术，但是在对于更好地保护环境有着极其重要的意义。
基于高斯烟羽模型的放射性气体的扩散
关于核电站泄漏放射性气体扩散的预估模型摘要由于核泄漏导致放射性气体扩散对经济和人身造成巨大损失的报道在国内外屡见不鲜，本文中日本福岛核泄漏事件更加使我们认识到对放射性气体扩散进行合理性的预估从而为以后类似于此的突发性事件作积极有效的补救措施的重要性。
对于问题一我们运用了点源烟羽扩散模型，用抛物型二阶偏微分方程解出理想状态下的不同时刻、不同地点的浓度表达式：222432 (,,,)(4)x y zktQC x y z t ektπ++-=。
此模型是建立在以泄漏点为圆心的一个无界球形区域内的。
为了使模型更符合实际情况，能够被应用于现实生活中，我们在泄漏源有效高度的确定和考虑地面反射与吸收作用下对此模型进行了修正，最终得到问题一浓度的确定公式(14)(,,,)C x y z t的表达式。
对于问题二，我们采用高位连续点源烟羽扩散模式，其扩散服从正态分布，并根据概率论的相关知识通过数学公式推导，得到理想状态下的高斯模型，由泄漏源有效高度，地面反射等因素的影响对其进行修正，又由于重力干沉积，雨洗湿沉积以及核衰变等因素对源强的影响，对高斯烟羽模型再次进行修正，最终得到泄漏源周边浓度变化情况即公式(32)，在风速为k m/s的条件下浓度为(,,,)C x y z H。
对于问题三，我们在第二问建立的模型的基础上，引入时间变量rt和t，和扩散速度变量s，在风速和扩散速度的共同影响下，可分别求出上风向和下风向浓度预估模型即公式(40)和(41)。
对于问题四，本文参阅整理大量气象、地理、新闻资料，选择我国东海岸典型地域---山东半岛和美国西海岸典型地域---加利福尼亚州作为研究对象，综合考虑对应海域平均风速及风向、地理距离、海水对放射性物质扩散的部分反射系数等因素，并通过计算机模拟，预测出放射性核物质将经过6天到达我国东海岸，且131I浓度预测值为：0.1053mBq m-⋅,，经过6.8天到达美国西海岸，且氙-133浓度的预测值几乎为零，与实际情况比较吻合。
基于FICK定律和高斯烟羽模型的放射性气体扩散研究
基于FICK定律和高斯烟羽模型的放射性气体扩散研究放射性气体扩散是核能安全领域的重要研究课题之一，研究其扩散规律有助于评估周围环境的辐射水平，制定相应的应对措施。
本文将基于FICK定律和高斯烟羽模型，对放射性气体扩散进行研究。
FICK定律是描述气体、液体或固体中物质扩散的基本定律。
根据FICK定律，物质的扩散速率正比于物质浓度梯度的变化率，与物质的分子扩散系数成正比。
对于放射性气体的扩散，可以用FICK定律表示为：J=-D*(∂C/∂x)其中，J为扩散通量，D为扩散系数，C为浓度，x为扩散距离。
高斯烟羽模型则是用于预测大气中污染物传输的经典模型。
根据该模型，气体扩散呈现高斯分布，其浓度随距离的增加呈指数衰减。
具体来说，高斯烟羽模型可以表示为：C(x, y, z) = (Q/(2πσu)) * exp(-(x-x0)^2/(2σx^2) - (y-y0)^2/(2σy^2) - (z-z0)^2/(2σz^2))其中，C为扩散浓度，Q为释放速率，(x0,y0,z0)为源的位置，σx、σy和σz分别为扩散系数。
基于上述理论，针对放射性气体扩散研究，可以首先确定材料的扩散系数。
放射性气体通常是从核电站、核工厂等活动中释放出来的，因此首先要进行放射性气体浓度的测量，以便计算扩散系数，然后可以利用FICK定律进行扩散速率的估算。
然后，可以利用高斯烟羽模型进一步研究放射性气体的扩散规律。
首先需要确定放射性气体的释放速率和源位置，然后利用高斯烟羽模型计算不同点的浓度。
通过浓度的计算，可以得到放射性气体在空间中的分布情况，以及随着距离的增加浓度的衰减情况。
最后，将根据模型计算出的数据与实际测量数据进行对比，以验证所建模型的准确性和可靠性。
如果模型与实测数据吻合良好，则可以通过该模型来预测放射性气体的扩散情况，为相关工程和环境保护提供科学依据。
总之，基于FICK定律和高斯烟羽模型的放射性气体扩散研究可以提供对该气体扩散规律的理论解释和预测，并为核能安全领域的决策制定提供科学依据。
放射性的气体预测模型
放射性气体扩散的预测摘要 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 σσσσσπ+-+--∙- 在应用高斯模型时最关键的是确定扩散参数，再应用程序计算出各处的质量浓度。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
关于核电站泄漏放射性气体扩散的预估模型摘要由于核泄漏导致放射性气体扩散对经济和人身造成巨大损失的报道在国内外屡见不鲜，本文中日本福岛核泄漏事件更加使我们认识到对放射性气体扩散进行合理性的预估从而为以后类似于此的突发性事件作积极有效的补救措施的重要性。
对于问题一我们运用了点源烟羽扩散模型，用抛物型二阶偏微分方程解出理想状态下的不同时刻、不同地点的浓度表达式：222432 (,,,)(4)x y zktQC x y z t ektπ++-=。
此模型是建立在以泄漏点为圆心的一个无界球形区域内的。
为了使模型更符合实际情况，能够被应用于现实生活中，我们在泄漏源有效高度的确定和考虑地面反射与吸收作用下对此模型进行了修正，最终得到问题一浓度的确定公式(14)(,,,)C x y z t的表达式。
对于问题二，我们采用高位连续点源烟羽扩散模式，其扩散服从正态分布，并根据概率论的相关知识通过数学公式推导，得到理想状态下的高斯模型，由泄漏源有效高度，地面反射等因素的影响对其进行修正，又由于重力干沉积，雨洗湿沉积以及核衰变等因素对源强的影响，对高斯烟羽模型再次进行修正，最终得到泄漏源周边浓度变化情况即公式(32)，在风速为k m/s的条件下浓度为(,,,)C x y z H。
对于问题三，我们在第二问建立的模型的基础上，引入时间变量rt和t，和扩散速度变量s，在风速和扩散速度的共同影响下，可分别求出上风向和下风向浓度预估模型即公式(40)和(41)。
对于问题四，本文参阅整理大量气象、地理、新闻资料，选择我国东海岸典型地域---山东半岛和美国西海岸典型地域---加利福尼亚州作为研究对象，综合考虑对应海域平均风速及风向、地理距离、海水对放射性物质扩散的部分反射系数等因素，并通过计算机模拟，预测出放射性核物质将经过6天到达我国东海岸，且131I浓度预测值为：0.1053mBq m-⋅,，经过6.8天到达美国西海岸，且氙-133浓度的预测值几乎为零，与实际情况比较吻合。
关键词点源烟羽扩散模式偏微分方程 P-G曲线高斯修正模型仿真一、问题重述现有一座核电站遭遇自然灾害发生泄漏，一种浓度为p的放射性气体以速度m kg/s排出，在无风情况下以速度s m/s匀速向大气四周扩散。
问题一：建立一个数学预测模型来描述核电站周边不同距离地区、不同时段放射性物质浓度。
问题二：给出当风速为k m/s时核电站周边放射性物质浓度的变化情况。
问题三：当风速为k m/s时，建立一个放射性气体的预测模型计算出上风和下风L公里处的浓度。
问题四：利用所建立的模型来计算出福岛核电站的泄漏对我国东海岸，及美国西海岸的影响。
二、问题分析对于问题一，在无风的情况下，匀速在大气中向四周扩散,的速度为s m/s，建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。
对于此问题我们采用放射气体扩散是点源连续扩散模型，首先是在理想的状态下，将环境视作了一个没有边界的空间，建立高斯模型，利用“质量守恒定律”和二阶偏微分方程来描述烟雾浓度的变化规律，求得核电站不同距离、不同时段放射性物质浓度的预估模型。
为了更接近实际，我们对高斯模型进行了改进，考虑泄漏源的有效高度，以及地面和海水对放射性气体的反射，使模型更加接近实际情况。
对于问题二，在外界环境有风，风速为k m/s的情况下,求出核电站周边放射性物质浓度的变化情况。
假设放射性气体在x轴的扩散呈正态分布，利用概率和统计的知识，通过数学推导求出高架连续点源烟羽扩散模式的高斯修正模型中的放射性气体物质浓度的变化函数，考虑扩散过程中地面和海水的反射、泄漏源有效高度、干沉积、雨洗湿沉积和放射性物质的衰变，使高斯模型进一步修正。
对于问题三，当风速为k m/s时，利用求得的模型，分别求出上风和下风L公里处，t时间时，放射性物质浓度的预测模型。
采用第二问的结论，求出t时间内的初始浓度，再加入时间t的变量，得出上风向和下风放射性气体在r向的不同位置，不同高度，不同时间的放射性气体的浓度函数。
对于问题四，参考互联网上的相关资料。
查出当天的天气变化，风速，地理环境，以及福岛核电站爆炸时的相关参考资料。
取中国东海岸距离日本最近的山东半岛为例，美国西海岸以加利福尼亚州为例，进行研究。
综合考虑对应海域平均风速及风向、地理距离、海水对放射性物质扩散的部分反射系数等因素，预测出放射性核物质对中国东海岸和美国西海岸的影响。
三、问题假设1.放射性气体初始泄漏时可看作在空中某一点向四周等强度的瞬时释放，在于风向垂直的三维坐标的两个方向上y 轴和z 轴的分布呈正态分布。
2.初始时刻放射性气体云团的内部、温度呈均匀分布。
3.扩散过程中不考虑云团内部温度的变化，忽略热传递、热对流及热辐射。
4.泄漏气体是理想气体，遵守理想气体状态方程。
5.整个扩散过程中风速的大小、方向保持不变。
6.地面对放射性气体不完全反射。
7.放射性气体的传播服从扩散定律，即单位时间通过单位法向面积的流量与它的浓度梯度成正比。
四、符号的定义和说明(一)符号定义(,,,)C x y z t ：时刻t 无穷空间中任一点(,,)x y z 的烟雾浓度q ：单位时间通过单位法向面积的流量Ω：空间域V ：Ω的体积S ：包围Ω的曲面n ：S 的外法线向量σ：扩散系数1Q ：在[,]t t t +∆内通过空间域Ω的流量2Q ：空间域Ω内放射性气体的增量Q ：核泄漏释放的放射性气体总量(,,)x y z δ：单位强度的点源函数(二)名词解释1.烟羽：从工厂烟囱中连续排放出来的烟体，外形呈羽毛状，因而得名。
烟羽可被看作是由无数个时间间隔为无限短暂的、依次排放的烟团所组成；烟团各部分的运动速率不同。
2.热力抬升：排放物初始温度高于大气环境温度所造成的烟云升高。
若不计烟气与空气成分而造成的密度差异，温度差便决定了作用于烟云的净浮力。
实测资料表明，烟云抬升主要是热力因素的作用。
3.地面反射作用下的浓度处理：由于地面的存在，烟羽的扩散是有界的，可以把地面看作一面镜子，对泄漏气体起反射作用，并采用像源法进行处理，任一点的浓度可以看作两部分的贡献之和。
一部分是不存在地面时所造成的泄漏物浓度，一部分是由地面反射作用增加的浓度。
4.干沉积：在扩散过程中同时有重力沉降的位移迭加到羽流中心线上，中心线就会向下倾斜，所有粒子相当于在下倾的中心线上扩散。
5.雨洗作用：降雨对烟羽中的颗粒物及气溶胶具有清洗作用，可溶性气体与蒸汽亦可溶于雨水中，降雨过程造成的这类湿沉积是导致放射性气溶胶和气体向地面沉积的另一重要机制。
五、模型的建立和求解第一部分准备工作1.3月12日22时，世界气象组织和国际原子能机构北京区域环境紧急响应中心(设在中国气象局国家气象中心)综合分析最新气象资料和放射性污染物扩散模式运行结果权威发布：由于我国位于日本国的西部，日本核泄露发生地区近日风向盛行由西向东，其核泄漏放射性污染物未来3天对我国没有影响。
2.3月13日本上空为西风，未来仍以偏西风为主，放射性污染物主要向日本东部的北太平洋区域扩散。
专家表示，中国位于日本国的西面，中间又有日本海、朝鲜半岛、黄海、东海相隔，因此，在以向北太平洋扩散为主的日本核污染物扩散区域距离中国较远，未来3天内对中国没有影响。
3.日本于2011年3月15日下午，世界气象组织和国际原子能机构北京区域环境紧急响应中心分析认为：日本中北部区域在中低层大气中的风向由西南风转为西北风；高空大气主要以偏西风气流为主，近期由于降水发生，有利于核物质沉降，影响范围缩小。
未来三天（16日至18日），日本核电站核泄漏产生的放射性污染物主要影响区域为日本中部、北部及其以东的北太平洋区域。
4.美国环保署和能源部18日联合发表声明说，美国西海岸已发现可能来自日本的放射性物质，不过剂量“极小”，仅是岩石、砖块和太阳等天然辐射源辐射量的十万分之一。
声明说，能源部设在加利福尼亚州萨克拉门托的监测站18日监测到碘、碲以及铯元素的放射性同位素，不过“量极小，对人体健康没有影响”。
此外，能源部设在华盛顿州太平洋西北国家实验室的监测器16日、17日监测到“微量”的氙-133，对人体健康也没有影响。
5.半岛网3月16日消息日本处于青岛的东北方向，根据中央气象台分析，未来三天日本中北部区域中低空风向由西南风转，西北风，高空以偏西风气流为主，预计主要影响区域为日本本州岛的福岛核电站附近地区，以及日本以东的北太平洋区域。
根据最新气象资料分析，未来一周青岛市1000米以上高空盛行偏西风，近地面16日至19日南风4到5级，19日夜间到23日持续受冷空气影响，北风5到6级阵风7级，海上6到7级阵风8级。
第二部分问题一的模型(一)放射性气体浓度的变化规律图1 模型的理想化图形将核泄漏时刻记作0t =，爆炸点选为坐标原点，时刻t 无穷空间中任一点(,,)x y z 的气体浓度记为(,,,)C x y z t 。
根据假设7，单位时间通过单位法向面积的流量为:q gradC σ=-⋅ (1)其中σ是扩散系数，grad 表示梯度，负号表示由浓度高向浓度低的地方扩散。
考察空间域Ω，Ω的体积为V ，包围Ω的曲面为S ，S 的外法线向量为n ，则在[,]t t t +∆内通过Ω的流量为：1t tt SQ qnd dt σ+∆=⎰⎰⎰ (2) 而Ω内放射性气体的增量为：2[(,,,)(,,,)]VQ C x y z t t C x y z t dV =+∆-⎰⎰⎰ (3)由质量守恒定律可知：12Q Q = (4)然后根据曲面积分的奥式公式可以得到：S Vqnd divqdV σ=⎰⎰⎰⎰⎰ (5)其中div 是散度记号。
由(1)-(5)式再利用积分中值定理不难求得下式：222222()(),0,,,C C C C div gradC t x y z t x y zσσ∂∂∂∂==++>-∞<<+∞∂∂∂∂ (6) 这是无界区域的抛物型偏微分方程。
根据假设1，初始条件为作用在坐标原点的点源函数，能够将其记作：(,,,0)(,,)C x y z Q x y z δ= (7)其中Q 即为核泄漏释放的放射性气体总量，(,,)x y z δ式单位强度的点源函数。
在(7)的约束条件下可以解得方程(6)的解为：22232(,,,)exp()(4)4Q x y z C x y z t t tπσσ++=- (8) 这个结果表明，对于任何时刻t 放射性气体浓度C 的等值面是以初始时刻释放点为球心的无限球面上，并且随着球面半径的增加C 的值是连续减少的。
放射性气体在大气中的扩散是服从一定的物理规律的，但是仅从这些规律得到的结果是不能解释实际生活中放射性气体真正的扩散过程，也不能够得到任意时刻、任意地点放射性气体的浓度。