基于高斯模型的放射性物质扩散模型
- 格式:doc
- 大小:1.08 MB
- 文档页数:18
核能作为可持续能源,在高新科技领域地位愈发重要。
但对核能技术开发的规范却一直没能跟上。
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. 建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。
放射性气体扩散的预估模型摘要日本福岛核电遭遇日本近海9.0级地震并引发了大海啸的破坏,发生了核泄漏,核事故发生后放射性物质经由大气的扩散过程, 可在短时间内对环境产生大范围的影响【1】。
因此, 在进行核事故的后果评估和应急决策时,大气扩散过程的模拟是非常重要的。
本文中我们运用大气扩散随机游走模型、高斯烟羽优化模型以及多元线性回归分析模型,结合MATLAB、SPSS软件进行计算,最终建立了理想中的模型。
问题一中,我们建立了大气扩散随机游走模型,对扩散物的示踪粒子进行跟踪,形成大量的随机游走轨道。
对穿过单位体积的轨道数和粒子滞留时间进行统计计算即可获得不同距离地区、不同时段放射性物质的浓度。
问题二中,我们基于大气污染的经典高斯烟羽模型,以期实现在风速k sm/的动态预测,但由于高斯烟羽模型的限制条件太过于理想,我们就充分考虑气体扩散过程中的风速、地面反射、地面粗糙程度等因素对传统的高斯烟羽模型进行优化,据此分析了核电站周边放射性物质浓度的变化情况。
问题三中,根据江苏工业学院油气储运安全综合实验模型平台,利用示踪技术,模拟了有害气体瞬时扩散的整个过程的实验,基于此我们又结合多元线性回归的基本方法,建立放射性气体扩散的多元线性回归模型。
问题四中,我们搜索整理了大量气象、地理、新闻资料,分析了日本福岛发生核泄漏时风向、风速等重要因素,并选取了我国东海岸山东青岛和美国西海岸洛杉矶作为研究对象,综合考虑两者距离扩散源的距离、主要地形等因素进行分析,运用已建立的高斯烟羽优化模型得出最终结果。
关键词:随机游走大气扩散模型线性回归分析模型高斯烟羽优化模型一、问题重述本文是以日本福岛核电站遭遇自然灾害发生核泄漏的背景而提出的。
为了便于为事故提供积极的补救措施,对放射性核污染物大气扩散过程的模拟非常重要。
设有一座核电站在遭遇自然灾害后浓度为0 的气体匀速泄漏,速度为s mkg /,在无风的条件下,气体以s s m /的速度在大气中扩散。
放射性污染物的迁移与扩散模拟近年来,随着人们生活水平的提高和工业化的不断发展,放射性污染也成为了人们关注的重要问题之一。
放射性污染物和核事故的发生不仅会对人类和环境造成极大的危害和影响,更会引起全球范围内的关注和警惕。
如何准确地模拟放射性污染物的迁移与扩散,成为了科学家们研究的重大课题。
一、放射性物质的污染来源放射性物质主要来自两个方面:自然放射性和人工放射性。
自然放射性是指自然界中存在着的放射性物质,比如岩石、土壤、气体等,而人工放射性则是指人类活动造成的放射性污染,如核电站、核试验等。
二、放射性物质的迁移与扩散放射性物质的迁移和扩散是指放射性物质在环境中移动和扩散的过程。
这个过程是很复杂的,受到许多因素的影响,如地形、气象、底物等,而这些因素又相互交织,形成了复杂的生态条件。
在放射性物质的迁移和扩散的过程中,主要受到以下几种因素的影响:1.气象条件:气象条件对放射性物质的迁移和扩散起着极其重要的作用。
风向、气温、湿度等因素都会影响放射性物质的扩散方式和速度。
2.地形地貌:不同的地形和地貌特征也是影响放射性物质迁移和扩散的重要因素之一。
山脉、丘陵和平原等地域区分之间有很大的差异。
3.土壤性质:土壤不仅影响放射性物质的迁移和扩散,而且也是放射性物质富集和释放的重要媒介。
土壤质地、声速、通透性等因素都会影响放射性物质的传播和释放。
三、放射性物质的模拟方法针对放射性物质的复杂迁移和扩散过程,科学家们根据其具体情况,制定了不同的模拟方法。
其中最常用的是数值模拟方法及推理模型等。
1.数值模拟方法数值模拟方法是用计算机来模拟放射性物质的迁移和扩散过程。
该方法主要依据一定的物理方程和相应的数学方法,利用计算机模拟流体在多相介质中的运动和传输过程,具有清晰、可视化等优点,能够提供具有空间和时域信息的放射性物质分布和状况等方面的数据。
2.推理模型推理模型是指在放射性物质的迁移和扩散过程中,根据一定的评价标准,找出放射性物质的扩散规律,并通过统计分析得到相关参数。
基于高斯修正模型的放射性气体扩散浓度预测摘要本文是以日本福岛核电站遭遇自然灾害发生核泄漏的背景而提出的。
对于问题一,考虑到放射性物质的泄漏是连续不断的。
本文根据“泄漏放射性物质质量守恒定律”和“气体泄漏连续性原理”建立了微积分方程,应用了高等数学中散度、梯度、流量等数学概念,通过Guass公式、四维二阶偏微分方程,因而得到了核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。
同时为使模型适用范围更广,本文引入了地面反射系数,考虑了由于放射性物质从泄漏口喷出时具有初动量而使其泄漏源有效高度被抬高等因素,进而得到了在无风环境中适用范围更广的“高斯修正模型”。
对于问题二,要探究风速对放射性物质浓度分布的影响。
本文运用概率学[1]知识,通过图解和数学推导得出“连续点源放射性物质高斯扩散模型”。
本文依次考虑了“重力沉积”、“雨水沉积”、“核衰变”等因素对浓度分布的影响。
并通过构建“耗减因子”、“衰变因子”等方法将耗减和衰变的放射性物质“投影”到泄漏源浓度中,得到了经多次合理修正后的“优化高斯模型”,并据此分析了泄漏源周边地区放射性物质的浓度变化。
针对于问题三,本文在问题二的基础上,结合考虑风速和放射性物质扩散速度在空间中的矢量运算。
得出在对上风口分析时,要分类讨论风速和自然扩散速度之间的大小关系,当风速小于自然扩散速度时,放射性物质是无法到达上风口的。
对于问题四,本文参阅整理大量气象、地理、新闻资料,选择我国东海岸典型地域---山东半岛作为研究对象,综合考虑对应海域平均风速及风向、地理距离、海水对放射性物质扩散的部分反射系数等因素,并通过C++编程模拟计算,预测出放射性核物质将经过6.5天到达我国东海岸,且0.100mBq⋅m-3与实际情况比较吻合。
131I浓度预测值为:关键词:放射性气体扩散浓度变化高斯修正模型预测1 问题的提出由于重大的突发性核泄漏紧急灾害事件具有爆发性、空间分布不连续性、对周边地形和气象条件的敏感性的特点,研究核事故所释放的物质的时空分布需要高度精确的技术,但是在对于更好地保护环境有着极其重要的意义。
基于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定律和高斯烟羽模型的放射性气体扩散研究可以提供对该气体扩散规律的理论解释和预测,并为核能安全领域的决策制定提供科学依据。
基于核电站泄漏放射性气体扩散的预估模型摘要由于核泄漏导致放射性气体扩散对经济和人身造成巨大损失的报道在国内外屡见不鲜,本文中日本福岛核泄漏事件更加使我们认识到对放射性气体扩散进行合理性的预估从而为以后类似于此的突发性事件作积极有效的补救措施的重要性。
对于问题一我们运用了点源烟羽扩散模型,用抛物型二阶偏微分方程解出理想状态下的不同时刻、不同地点的浓度表达式: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 定律求出扩散系数D ,进而得出扩散的粒子流量与其浓度梯度的正比关系,得出描述扩散情况的偏微分方程。
随后,用傅里叶变换求解得到扩散方程的解,并在MATLAB 中绘制此微分方程的图形,发现预测图形与东京市测得的实际数据的图形基本吻合,即离泄漏源越远浓度越低。
对于问题二,要探究风速对放射性物质浓度分布的影响。
风速的处理是此问题的核心,因此我们采用大气污染的经典高斯扩散模型。
通过查阅相关资料,我们发现连续点源的平均气体流,其浓度分布符合正态分布规律,因此污染物浓度在y 、z 轴上的分布为正态分布。
取烟云轴线为x 方向(平均风向),得出无界情况下下风向空间某一点的浓度函数(,,,)C x y z H 。
同时考虑到泄漏点的高度,对模型进行了修正,得出的浓度分布情况与实际情况相似 。
对于问题三,由问题二得到的浓度分布函数(,,,)C x y z H ,可以分上风和下风两种情况,根据当时的实际情况,假设自然风速大于泄露的自身扩散速度,则可将第二问中的风速u 替换成k 与s 的线性组合即可,即下风向的速度为k s +,上风向的速度为k s -。
将平均风速u 分别代入浓度函数,得出了上风向和下风向L 处浓度分布函数。
结果显示两个地点气体浓度变化情况与第二问得出的结果一致。
对于问题四,本文参考了大量地理、气象、专题报道等资料。
假设风向不随时间变化而变化,即我国东海岸一直处于上风,美国一直处于下风,并且取大气稳定度均在D 时的扩散参数。
第43卷第1期2021年2月Vol.43No.1Feb.2021武汉理工大学学报(信息与管理工程版)JOURNAL OF WUT(INFORMATION&MANAGEMENT ENGINEERING)文章编号:2095-3852(2021)01-0001-08文献标志码:A基于Lagrangian-Eulerian耦合的放射性物质大气扩散模型研究邰扬,申世飞(清华大学工程物理系,北京100084)摘要:针对中尺度扩散下拉格朗日模型和欧拉模型在计算效率和准确度方面的矛盾,建立了适用于模拟核事故下放射性物质大气扩散的基于拉格朗日模型与欧拉模型耦合的预测模型,以国内某核电站为例进行了放射性物质扩散的虚拟算例模拟,对模型进行验证和分析。
结果表明:拉格朗日-欧拉耦合模型可以较好地模拟放射性物质在大气中的扩散,不仅可以提升中尺度和大尺度下的计算效率,还可以有效减小欧拉模型初始条件所带来的系统误差。
关键词:核事故;大气扩散模型;拉格朗日-欧拉耦合模型;1311;核电站中图分类号:X946DOO:10.3963/j.issn.2095-3852.2021.01.001近年来,我国核电建设持续增长,核电在建规模长期保持全球第一,按目前已形成的核电发展条件预测,2021—2035年每年将核准建设6-8台核电机组,2035年我国核电装机规模将达到1.5-1.8亿kW,核电发电量占比将提升至10%以上,核能在国民经济中将发挥越来越重要的作用⑴。
保证核能的安全使用是开发和利用核能的前提,对核电来说,除了通过合理选址进行前期安全规划和建立并维护行之有效的安全运营体系之外,对核事故发生后放射性物质的扩散过程进行模拟并根据结果制定安全预案,也是保障核电安全性的重要方面。
以2011年福岛核事故为例,大量的放射性物质持续性地排放到大气中,并在事故发生后的一个多月内随大气运动而大面积扩散,将周围大面积区域都拖入放射性污染的威胁中。
《CALPUFF模式用于放射性核素不同尺度的迁移扩散研究》篇一一、引言随着核能利用的普及,放射性核素迁移扩散的问题越来越受到关注。
在放射性核素的研究中,其不同尺度的迁移扩散是一个关键的科学问题。
近年来,计算流体动力学模型如CALPUFF在放射性核素迁移扩散的研究中得到了广泛的应用。
本文将详细探讨CALPUFF模式在放射性核素不同尺度的迁移扩散研究中的应用。
二、CALPUFF模式简介CALPUFF是一种基于高斯烟羽模型的气象模型,适用于模拟大气中污染物的长距离传输和扩散。
它利用大气环境中的物理、化学和气象条件来模拟放射性核素的迁移和扩散过程,并可根据不同尺度的研究需求调整模型参数。
三、CALPUFF模式在放射性核素迁移扩散研究中的应用1. 大尺度迁移扩散研究在大尺度迁移扩散研究中,CALPUFF模式能够有效地模拟放射性核素在大气中的长距离传输和扩散过程。
通过调整模型参数,可以更好地反映不同气象条件对放射性核素迁移扩散的影响。
此外,CALPUFF模式还可以与其他地理信息系统(GIS)相结合,实现空间数据的可视化处理和分析。
2. 中尺度迁移扩散研究在中尺度迁移扩散研究中,CALPUFF模式能够更精细地模拟放射性核素在局部地区的扩散过程。
通过建立详细的地理和环境信息数据库,可以更准确地模拟放射性核素在不同地形、气象条件下的迁移和扩散过程。
这对于评估局部地区放射性核素污染风险和制定相应的防护措施具有重要意义。
3. 小尺度迁移扩散研究在小尺度迁移扩散研究中,CALPUFF模式可以用于模拟放射性核素在微观尺度上的迁移和扩散过程。
这有助于研究放射性核素在土壤、水体等介质中的迁移和转化过程,为评估核事故等紧急情况下的放射性核素污染风险提供科学依据。
四、案例分析以某核事故为例,我们利用CALPUFF模式对该事故中放射性核素的迁移和扩散过程进行了模拟。
通过调整模型参数,我们成功地模拟了放射性核素在不同时间、不同空间尺度的迁移和扩散过程。
放射性气体扩散浓度预估模型【摘要】本文是以日本地震引起的福岛核电站的核泄漏为背景,并以给出的数据为基础,研究某一假设核电站的核泄漏问题。
我们通过收集相关的资料,并结合题目给出的数据,建立了高斯模型、连续点源高斯扩散模型解决了题目提出的四个问题。
针对问题一:考虑到泄漏源是连续、均匀和稳定的,我们运用散度、梯度、流量等数学概念,通过“泄漏放射性物质质量守恒”、“气体泄漏连续性定理”、 Guass 公式及积分中值定理得到了无界区域的抛物线型偏微分方程,然后再通过电源函数解出空间任意一点的放射性物质浓度的表达式,把此表达式定为模型一的前身。
鉴于放射性物质的扩散受到诸多因素的影响,如:泄漏源的实际高度、地面反射等。
我们以泄漏口为坐标原点建立三维坐标系,通过“像源法”处理地面反射对放射性物质浓度的影响,并由此对模型一的前身进行修正完善,得到模型一:高斯模型,即放射性物质浓度的预测模型。
最后我们模拟了放射性物质无风扩散仿真图。
针对问题二:当风速为k m/s 时,我们根据放射性核素云团在大气中迁移和扩散的数值计算的基本方法和步骤,并以泄漏点源在地面的投影点为坐标原点,以风向方向为x 轴,铅直方向为z 轴,与x 轴水平面垂直方向为y 轴建立三维坐标系,地面的反射作用同样利用“像源法”进行处理,得到连续点源高斯扩散模型。
考虑到地面反射、烟云抬升、放射性物质自身的沉降及雨水的吸附等对浓度的影响,我们对连续点源高斯扩散模型进行了修正,建立了修正的连续点源高斯扩散模型。
最后利用大气稳定度确定了扩散参数,进而求解了模型。
针对问题三:经分析,问题三的提出是以问题二为基础的,模型三的建立只需要将模型二加以调整即可。
我们以风速方向为x 轴正方向,将风速与放射性物质的扩散速度进行矢量运算,此问题则转化为求(,0,)L z 和(,0,)L z -两点处的放射性物质浓度,由此建立模型三,即上风和下风L 公里处放射性物质浓度浓度的预测模型。
针对问题四:首先,我们通过网络收集了相关数据,然后,我们结合模型二、模型三对数据进行整理代入,算出了日本福岛核电站泄漏的放射性物质扩散到中国东海岸和美国西海岸的浓度分别为334.242910/g m -⨯、432.385410/g m -⨯。
关键词:高斯模型 连续点源高斯扩散模型 核泄漏一问题的重述1.1问题背景目前,核电站的发展能带来巨大的经济效益和社会效益,但核电站一旦发生核泄漏,将会给人们的生命健康和周边环境带来巨大的危害性影响。
2011年3月日本的福岛核电站的放射性气体的核泄漏事件更让我们关注放射性气体泄漏时的浓度问题。
因此,正确的测出大气中放射性物质的浓度在环境监测和安全评估中具有重要意义1.2问题提出p的放射性气体以速度m kg/s 有一座核电站遇自然灾害发生泄漏,浓度为匀速排出,在无风的情况下,以速度s m/s匀速在大气中向四周扩散.1)在无风的情况下,建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。
2)当风速为k m/s时,给出核电站周边放射性物质浓度的变化情况。
3)当风速为k m/s时,分别给出上风和下风L公里处,放射性物质浓度的预测模型。
二符号说明三模型假设1、扩散过程中浓度在y、z轴上的变化分布是高斯分布。
2、放射性物质的扩散看作是空间某一连续点源向四周等强度地瞬时释放放射性物质,放射性物质在无穷空间扩散过程中不发生性质变化,且不计地形影响。
3、放射性物质扩散服从扩散定律,即单位时间通过单位法向面积的流量与它的浓度梯度成正比。
4、放射性物质在穿过降雨区域时,其强度由于雨水的吸收而减少,减少比率为常数。
5、假设地面对放射性气体起全反射作用,地面和海面对放射性气体没有吸附,将海面视为平原地区6、假设风向为水平风向,且风向风速不随时间变化。
7、扩散过程中不考虑泄漏点内部温度的变化对气体扩散的影响。
四问题的分析4.1问题(1)的分析p放射性气体以速度m kg/s匀速一座核电站遇自然灾害发生核泄漏,浓度为排除,这近似于放射性物质源是连续均匀稳定的。
在无风情况下,放射性气体以速度s m/s匀速在大气中向四周扩散,放射性气体的扩散服从扩散定律,即单位时间通过单位法向面积的流量与它的浓度梯度成正比。
在这些条件下,我们明确了要研究的问题是点源连续泄漏的扩散问题,为了使建立的模型更加贴近实际,需考虑地面反射、核泄漏源的实际高度、降雨等因素对浓度分布的影响。
由“扩散定律”“放射性物质质量守恒定律”“气体泄漏连续性定理”可得出无界区域的抛物线型偏微分方程。
再通过假设条件建立未考虑地面反射、核泄漏源的实际高度、降雨等浓度影响因素的初步模型,然后从这些影响因素对模型进行完善,最终得出核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。
4.2 问题(2)的分析本问是探究风速为/km s时,核电站周边放射性物质浓度的变化情况。
当环境中空气流动产生风力时,在均匀湍流场中,扩散参数与下风向距离的关系是明确的,核泄漏时间较长时,可认为扩散是稳定的。
在下风向的湍流扩散相对于风力引起的移流相可忽略不计,在流动方向建立x轴,不考虑横向速度和垂直速度。
根据假设,空间中放射性物质的浓度服从高斯分布,可利用连续点源放射性物质的高斯扩散模型。
放射性物质在大气中扩散受诸多因素影响,考虑泄漏源有效高度、放射性物质自身重力产生的重力沉降、雨水的吸附等因素对放射性物质浓度的影响是必要的,通过这些影响因素对高斯模型进行修正,然后利用修正后的高斯模型探究核电站周边放射性物质浓度的变化情况。
4.3 问题(3)的分析本问是要求当风速为/km s时,建立上风和下风L公里处的放射性物质浓度的预测模型。
经分析,此问是问题(2)的延伸,我们只需建立合适的坐标系,将此问题转化为求具体两处的放射性物质的浓度,便能得出上风和下风L公里处的放射性物质浓度预测模型。
五模型的建立与求解5.1模型一的建立与求解5.1.1模型一的初步建立以核泄漏点正下方的地面为坐标原点(0,0,0),平均风向为X轴、指向下风方向,铅直方向为Z轴,水平垂直于风向轴(X轴)为Y向,建立空间坐标系,则核电站泄漏点O距有效地面的高度为H,则泄漏点位置坐标为(0,0,)O H。
图一:空间坐标系示意图将气体从泄漏源泄漏时刻记作t=0,时刻t 无穷空间中任意一点坐标为(x,y,z )的浓度记为C (x,y,z,t ),根据假设2,单位时间通过单位法向面积的流量与浓度梯度成正比,则:i q gradC δ=-⋅ (1)(,,)i i x y z δ=是扩散系数,grad 表示梯度,负号表示由浓度高向浓度低的地方扩散。
考察空间域Ω,其体积为V ,包围Ω的曲面为S ,S 为一规则的球面,S 外法线向量为(-,-,1)x y n z z=。
则在(,)t t t +∆内通过Ω的流量为: 1t t t sQ q nd dt σ+∆=⋅⎰⎰⎰ (2)Ω内放射性物质的增量为:2[(,,,)(,,,)]V Q C x y z t t C x y z t dV =+∆-⎰⎰⎰ (3)从泄漏源泄漏的放射性物质的总量为:00t tt Q p dVdt +∆Ω=⎰⎰⎰⎰(4)根据“质量守恒定律”和“气体泄漏连续性原理”,单位时间内通过所选曲面S 的向外扩散的放射性物质的量与S 曲面内放射性物质增量之和,等于泄漏源在单位时间内向外泄漏的放射性物质。
则:012Q Q Q =+ (5)即:0[(,,,)(,,,)]t t t t t t V s C x y z t t C x y z t dV q nd dt p dVdt σ+∆+∆Ω+∆-+⋅=⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰(6)又根据曲面积分的Gauss 公式:s V q nd divqdV σ⋅=⎰⎰⎰⎰⎰ (其中div 是散度记号) (7) 0(,,,)(,,,)[]t t t t t t V V C x y z t t C x y z t t dV divqdVdt p dVdt t +∆+∆Ω+∆-∴⋅∆+=∆⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰00()(,,,)(,,,)lim lim t t tt t kdiv gradC dt C C x y z t t C x y z t t t t +∆∆→∆→∂+∆-∴==∂∆∆⎰由以上两式得:0[]V V C dV t divqdV t p dV t t Ω∂⋅∆+⋅∆=⋅∆∂⎰⎰⎰⎰⎰⎰⎰⎰⎰即为:0[]V V VC dV divqdV p dV t ∂+=∂⎰⎰⎰⎰⎰⎰⎰⎰⎰ (8)由以上公式并利用积分中值定理得:222222()(),0,,,i i y z C C C C div gradC t x y z t x y zδδδδ∂∂∂∂==+⋅+⋅>-∞<<∞∂∂∂∂ (9) 这是无界区域的抛物线型偏微分方程,根据假设1,初始条件为作用在坐标原点的电源函数,记作(,,,0)(,,)C x y z Q x y z δ= (10) 0(,,)Q x y z δ表示泄漏源漏泄释放的放射性物质总量,是单位强度的电源函数。
方程(9)满足方程(10)的解为:222403/2(,,,),(,,)(4)i x y z t i Q C x y z t e i x y z t δπδ++-==此模型只是在不考虑风速的情况下建立的,但为了使模型具有更加的实用性,下面我们将考虑泄漏源的实际高度、地面反射、降雨等因素对浓度的影响,完善模型。
5.1.2 模型一的最终建立1. 地面反射对模型的完善泄漏源有一定的高度,且泄漏点源是连续点源,则泄漏点源可视为高架连续点源,考虑到地面对扩散来的放射性气体有反射作用,根据假设4,地面对到达地面的扩散气体完全反射。
这儿可认为地面就像镜子一样,对放射性气体起全放射作用,可用“像源法”处理,如图3,建立三个坐标系,一是以泄漏源(实源)为坐标原点;二是以泄漏源在地面的投影点为原点,p 点是空间的任意一点,坐标为(,,)x y z ;三是以泄漏源关于地面的像对称源(像源)为原点。
把p 点放射性气体浓度看成两部分(实源与像源)作用之和。
图二: 高架连续点源扩散示意图从以上分析知,p 点放射性气体的浓度为实源和像源的放射性气体扩散至此点浓度的叠加。
则实际泄漏源(实源)对p 点的浓度贡献部分可用2()2z z H eσ--来表示;因为地面对扩散物质完全反射,则像对称源(像源)对p 点的浓度贡献部分可用2()2z z H e σ+-来表示。
于是对(11)式所修正完善的模型为: 2222()()()44403/21/2(,,,),(,,)(4)x y z x y z H z H t t t i mp C x y z t e i x y z t δδδπδ-++-+-== (12)我们自己模拟一组数据,利用matlab 进行仿真模拟,可实现该模型的模拟图像。