承诺书 我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。 我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。 我们参赛选择的题号是（从A/B/C/D中选择一项填写）： 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员(打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人(打印并签名)： 日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：
编号专用页 赛区评阅编号（由赛区组委会评阅前进行编号）： 全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：
基于系统综合评价的城市表层土壤重金属污染分析 摘要 本文针对城市表层土壤重金属污染问题，首先对各重金属元素进行分析，然后对各种重金属元素的基本数据进行统计分析及无量纲化处理，再对各金属元素进行相关性分析，最后针对各个问题建立模型并求解。 针对问题一，我们首先利用EXCEL 和 SPSS 统计软件对各金属元素的数据进行处理，再利用Matlab 软件绘制出该城区内8种重金属元素的空间分布图最后通过内梅罗污染 模型：2 /12 max 22?? ? ? ??+=P P P 平均综，其中平均P 为所有单项污染指数的平均值，max P 为土壤环境中 针对问题二,我们首先利用EXCELL 软件画出8种元素在各个区内相对含量的柱状图，由图可以明显地看出各个区内各种元素的污染情况，然后再根据重金属元素污染来源及传播特征进行分析，可以得出工业区及生活区重金属的堆积和迁移是造成污染的主要原因，Cu 、Hg 、Zn 主要在工业区和交通区如公路、铁路等交通设施的两侧富集，随时间的推移，工业区、交通区的土壤重金属具有很强的叠加性，受人类活动的影响较大。同时城市人口密度，土地利用率，机动车密度也是造成重金属污染的原因。 针对问题三，我们从两个方面考虑建模即以点为传染源和以线为传染源。针对以点为传染源我们建立了两个模型：无约束优化模型()[]()[]() 22y i y x i x m D -+-=，得到污染源的位置坐标()6782,5567；有衰减的扩散过程模型得位置坐标（8500,5500），模型为： u k z u c y u b x u a h u 222 2222222-??+??+??=??， 针对以线为传染源我们建立了l c be u Y ?-+=0模型，并通过线性拟合分析线性污染源的位置。 针对问题四，我们在已有信息的基础上，还应收集不同时间内的样点对应的浓度以及各污染源重金属的产生率。根据高斯浓度模型建立高斯修正模型，得到浓度关于时间和空间的表达式ut e C C -?=0。 在本题求解过程中，我们所建立的模型与实际紧密联系，有很好的通用性和推广性。但在求点污染源时，我们假设只有一个污染源，而实际上可能有多个点污染源，从而使得误差增大，或者使污染源的位置够不准确。 关键词 内梅罗污染模型 无量纲化 相关性 回归模型 高斯浓度模型
大气污染扩散 第一节大气结构与气象 有效地防止大气污染的途径，除了采用除尘及废气净化装置等各种工程技术手段外，还需充分利用大气的湍流混合作用对污染物的扩散稀释能力，即大气的自净能力。污染物从污染源排放到大气中的扩散过程及其危害程度，主要决定于气象因素，此外还与污染物的特征和排放特性，以及排放区的地形地貌状况有关。下面简要介绍大气结构以及气象条件的一些基本概念。 一、大气的结构 气象学中的大气是指地球引力作用下包围地球的空气层，其最外层的界限难以确定。通常把自地面至1200 km左右范围内的空气层称做大气圈或大气层，而空气总质量的98．2％集中在距离地球表面30 km以下。超过1200 km的范围，由于空气极其稀薄，一般视为宇宙空间。 自然状态的大气由多种气体的混合物、水蒸气和悬浮微粒组成。其中，纯净干空气中的氧气、氮气和氩气三种主要成分的总和占空气体积的99.97％，它们之间的比例从地面直到90km高空基本不变，为大气的恒定的组分；二氧化碳由于燃料燃烧和动物的呼吸，陆地的含量比海上多，臭氧主要集中在55～60km高空，水蒸气含量在4％以下，在极地或沙漠区的体积分数接近于零，这些为大气的可变的组分；而来源于人类社会生产和火山爆发、森林火灾、海啸、地震等暂时性的灾害排放的煤烟、粉尘、氯化氢、硫化氢、硫氧化物、氮氧化物、碳氧化物为大气的不定的组分。 大气的结构是指垂直（即竖直）方向上大气的 密度、温度及其组成的分布状况。根据大气温度在 垂直方向上的分布规律，可将大气划分为四层：对 流层、平流层、中间层和暖层，如图5－1所示。 1. 对流层 对流层是大气圈最靠近地面的一层，集中了大 气质量的75％和几乎全部的水蒸气、微尘杂质。受 太阳辐射与大气环流的影响，对流层中空气的湍流 运动和垂直方向混合比较强烈，主要的天气现象云 雨风雪等都发生在这一层，有可能形成污染物易于 扩散的气象条件，也可能生成对环境产生有危害的 逆温气象条件。因此，该层对大气污染物的扩散、输送和转化影响最大。 大气对流层的厚度不恒定，随地球纬度增高而降低，且与季节的变化有关，赤道附近约
关于金属汞扩散的问题 引言： 我们都知道，重金属丢弃到土地后会严重污染环境，同时对人体健康造成危害。著名的秦始皇陵墓，据专家在陵墓周围取数据观测，周围的汞含量呈现出外渗的趋势。也就是说，随着外围半径的扩大，汞含量浓度递减，并且随着时间的增加，汞渗透的半径越来越大。这就证明了汞金属在泥土中会发生扩散。因此，我们就提出，能否通过在外部取样的观察数据，建立一个数学模型，来判断陵墓中心处汞的浓度呢？ 模型的提出： 由于汞的扩散快慢跟本身的化学性质,物理性质有关。还有，由于在土堆里头，在各个方向上受到的力不相同和各种因素的影响，因此扩散的速度也会有差异。例如东西方向和南北方向会因为地球的自传而扩散速度会不一样。另一方面，汞在扩散的过程，由于泥土的吸收，化学反应等因数的影响，也会影响到汞的扩散。 为此我们引入一个函数u(x, y, z, t)，它表示t时刻在（x,y,z）处汞的浓度。我们的目标就是利用所观测到的数据，来推断出这个函数的表达式。 模型符号的引入： 为了表示汞在想x,y,z 方向上的扩散速度，我们在此引
入扩散系数： 2 a ：x 方向上的扩散系数 2 b ：y 方向上的扩散系数 2 c ：z 方向上的扩散系数 2 k ：由于泥土吸收，化学反应而引起的衰减系数 M ：扩散源汞的质量 模型假设： 1。假设有一汞扩散源，汞从扩散源沿 x ，y ，z 三个方向向四周扩散。 2。扩散前周围空间此物质的浓度为零。 3。扩散过程中没有人为因素的影响。 模型建立： u(x, y, z, t) 是 t 时刻点 (x, y , z) 处某物质的浓度。任取一个闭曲面 S ，它所围的区域是Ω，由于扩散，从 t 到 t t +? 时刻这段时间内，通过 S 流入Ω的质量为 1 M 2 2 2 1(cos cos cos )d d t t t S u u u M a b c S t x y z αβγ+????= ++???? ?? 其中 2 a ,2 b ,2 c 分别是沿 x ，y ，z 方向的扩散系数。 由高斯公式 ： ? ??? ?+Ω ??+??+??= t t t t z y x z u c y u b x u a M d d d d )(2 2 2 2 2 2 2 2 2 1
第一节大气污染物的扩散 一、湍流与湍流扩散理论 1. 湍流 低层大气中的风向是不断地变化，上下左右出现摆动；同时，风速也是时强时弱，形成迅速的阵风起伏。风的这种强度与方向随时间不规则的变化形成的空气运动称为大气湍流。湍流运动是由无数结构紧密的流体微团——湍涡组成，其特征量的时间与空间分布都具有随机性，但它们的统计平均值仍然遵循一定的规律。大气湍流的流动特征尺度一般取离地面的高度，比流体在管道内流动时要大得多，湍涡的大小及其发展基本不受空间的限制，因此在较小的平均风速下就能有很高的雷诺数，从而达到湍流状态。所以近地层的大气始终处于湍流状态，尤其在大气边界层内，气流受下垫面影响，湍流运动更为剧烈。大气湍流造成流场各部分强烈混合，能使局部的污染气体或微粒迅速扩散。烟团在大气的湍流混合作用下，由湍涡不断把烟气推向周围空气中，同时又将周围的空气卷入烟团，从而形成烟气的快速扩散稀释过程。 烟气在大气中的扩散特征取决于是否存在 湍流以及湍涡的尺度(直径)，如图5－7所示。 图5－7（a）为无湍流时，烟团仅仅依靠分子 扩散使烟团长大，烟团的扩散速率非常缓慢， 其扩散速率比湍流扩散小5～6个数量级；图5 －7（b）为烟团在远小于其尺度的湍涡中扩散， 由于烟团边缘受到小湍涡的扰动，逐渐与周边 空气混合而缓慢膨胀，浓度逐渐降低，烟流几乎呈直线向下风运动；图5－7（c）为烟团在与其尺度接近的湍涡中扩散，在湍涡的切入卷出作用下烟团被迅速撕裂，大幅度变形，横截面快速膨胀，因而扩散较快，烟流呈小摆幅曲线向下风运动；图5－7（d）为烟团在远大于其尺度的湍涡中扩散，烟团受大湍涡的卷吸扰动影响较弱，其本身膨胀有限，烟团在大湍涡的夹带下作较大摆幅的蛇形曲线运动。实际上烟云的扩散过程通常不是仅由上述单一情况所完成，因为大气中同时并存的湍涡具有各种不同的尺度。 根据湍流的形成与发展趋势，大气湍流可分为机械湍流和热力湍流两种形式。机械湍流是因地面的摩擦力使风在垂直方向产生速度梯度，或者由于地面障碍物(如山丘、树木与建筑物等)导致风向与风速的突然改变而造成的。热力湍流主要是由于地表受热不均匀，或因大气温度层结不稳定，在垂直方向产生温度梯度而造成的。一般近地面的大气湍流总是机械湍流和热力湍流的共同作用，其发展、结构特征及强弱决定于风速的大小、地面障碍物形成的粗糙度和低层大气的温度层结状况。 2. 湍流扩散与正态分布的基本理论 气体污染物进入大气后，一面随大气整体飘移，同时由于湍流混合，使污染物从高浓度区向低浓度区扩散稀释，其扩散程度取决于大气湍流的强度。大气污染的形成及其危害程度在于有害物质的浓度及其持续时间，大气扩散理论就是用数理方法来模拟各种大气污染源在
学号： 07412225 常州大学 毕业设计（论文） （2011届） 题目重气泄漏扩散模拟及应急救援 学生薛云龙 学院环境与安全工程学院专业班级安全072班 校内指导教师王新颖专业技术职务讲师 校外指导老师专业技术职务 二○一一年六月
重气泄漏扩散模拟及应急救援 摘要：重气泄漏扩散事故是经常发生且危害较大的一种事故形式，由于重气的密度大于空气，因此重气往往沿地面扩散，泄放物质进入人体将引起中毒事故，若泄放物质被点燃或引爆将引起大规模的燃烧爆炸事故。虽然人们对重气泄漏扩散所造成的危害十分重视，但由于缺乏足够有效的数据来提供人们作风险评估及预防改善措施，因此采用数学模型进行模拟是必要的。应在生产过程中，加强管理，强化生产者的安全生产教育。分析了泄漏扩散事故的七大影响因素，提取并建立了泄漏事故模式，并对各种事故模式的泄漏机理和发生条件进行了研究分析。通过试验研究得出在实际环境中大气主导风的风速，泄漏方向对气体扩散浓度分布有重大的影响，泄漏气体在下风向扩散的最快。静风时，随着时间的增加，空间各点的浓度有升高的趋势；在稳定风流中，空间各点的浓度随时间的变化不明显，可以认为是稳态的。泄漏的气体在下风向扩散的最快，在现场一旦发生天燃气泄漏，应综合考虑泄漏源的方向和该点当时的风向，风速等因素，及时准确预测泄漏气体可能扩散到危险区域，做好应对措施。 关键词：相似理论；泄漏模型；泄漏扩散模式；示踪法；重气；应急救援；
Heavy gas leak dispersion modeling and emergency rescue Abstract ： As it is well-known, many industrial and domestic gases are toxic and flammable are stored in highly-pressurized vessels at liquefied state with ambient temperature. If there is by chance a sudden release, it often forms heavy-than-air vapour. The accident release and dispersion of toxic and flammable heavy gas can present a serious ris k to the public’s safety and to the environment. Disease may be caused when the flammable heavy gases are lit. Although great attention has been paid to the hazard of heavy gas dispersion, effective data of filed experiments are still insufficient to make risk assessment and precaution. Through the statistical analysis, draw a conclusion that chemical system in production, transportation and storage process, should first consideration and control of hazardous chemicals, and summarizes the characteristics of the leak diffusion process performance. Subjective factors, equipment inherent defect caused by leakage on China's chemical system is the main reason of the accident. In the process of production, should be strengthen management, strengthen the education of production safety producer. Analysis of the seven factors affecting diffusion of leakage accident, to extract and established the patterns of the leakage accident, and various and leakage accident modes mechanism and the conditions were studied and analyzed. Through the experimental study on practical environment atmosphere that dominated the wind, the wind of gas leakage direction spread concentration distribution, has enormous influence on the spread of gas leakage next wind fastest. Static, as time flies, the space increased concentration of the each point of the trend. In the stable romantic, space the concentration of each point does not change significantly over time, can be considered a steady. Leak gas diffusion next wind fastest, on the site once produce natural gas leak, should be taken into account in the direction and point source leaking the wind direction, wind speed at factors such as timely and accurate prediction leakage, gas may be spread to dangerous area, completes the countermeasures. Key words：Theory of similarity; Leakage model；Leakage diffusion mode；Trace method； heavy gas；Emergency rescue
关于核电站泄漏放射性气体扩散的预估模型 摘要 由于核泄漏导致放射性气体扩散对经济和人身造成巨大损失的报道在国内外屡见不鲜，本文中日本福岛核泄漏事件更加使我们认识到对放射性气体扩散进行合理性的预估从而为以后类似于此的突发性事件作积极有效的补救措施的重要性。 对于问题一我们运用了点源烟羽扩散模型，用抛物型二阶偏微分方程解出理 想状态下的不同时刻、不同地点的浓度表达式： 222 4 32 (,,,) (4) x y z kt Q C x y z t e kt π ++ - =。 此模型是建立在以泄漏点为圆心的一个无界球形区域内的。为了使模型更符合实际情况，能够被应用于现实生活中，我们在泄漏源有效高度的确定和考虑地面反射与吸收作用下对此模型进行了修正，最终得到问题一浓度的确定公式 (14)(,,,) C x y z t的表达式。 对于问题二，我们采用高位连续点源烟羽扩散模式，其扩散服从正态分布，并根据概率论的相关知识通过数学公式推导，得到理想状态下的高斯模型，由泄漏源有效高度，地面反射等因素的影响对其进行修正，又由于重力干沉积，雨洗湿沉积以及核衰变等因素对源强的影响，对高斯烟羽模型再次进行修正，最终得到泄漏源周边浓度变化情况即公式(32)，在风速为k m/s的条件下浓度为(,,,) C x y z H。 对于问题三，我们在第二问建立的模型的基础上，引入时间变量 r t和t，和 扩散速度变量s，在风速和扩散速度的共同影响下，可分别求出上风向和下风向浓度预估模型即公式(40)和(41)。 对于问题四，本文参阅整理大量气象、地理、新闻资料，选择我国东海岸典型地域---山东半岛和美国西海岸典型地域---加利福尼亚州作为研究对象，综合考虑对应海域平均风速及风向、地理距离、海水对放射性物质扩散的部分反射系数等因素，并通过计算机模拟，预测出放射性核物质将经过6天到达我国东海岸，且131I浓度预测值为：0.1053 mBq m- ?,，经过6.8天到达美国西海岸，且氙-133浓度的预测值几乎为零，与实际情况比较吻合。 关键词点源烟羽扩散模式偏微分方程 P-G曲线高斯修正模型仿真
赛区评阅编号（由赛区组委会填写）： 2015高教社杯全国大学生数学建模竞赛 承诺书 我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。 我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。 我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。 我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。 我们参赛选择的题号（从A/B/C/D中选择一项填写）： C 我们的报名参赛队号（12位数字全国统一编号）： 参赛学校（完整的学校全称，不含院系名）：温州医科大学 参赛队员 (打印并签名) ：1. 章成俊 2. 杨超 3. 谢锦 指导教师或指导教师组负责人 (打印并签名)： 日期：年月日
赛区评阅编号（由赛区组委会填写）： 2015高教社杯全国大学生数学建模竞赛 编号专用页 送全国评阅统一编号（由赛区组委会填写）： 全国评阅随机编号（由全国组委会填写）：
对垃圾处理厂污染的动态监控及居民补偿 摘要 城市垃圾处理问题是一个世界性难题。目前垃圾焚烧正逐步成为中国垃圾处理的主要手段之一。本论文构根据题目设置的垃圾处理厂规模，建立了环境动态监控体系，并根据潜在污染风险对周围居民进行了合理经济补偿的设计。 对于问题（1），为了实现对垃圾焚烧厂烟气排放及相关环境影响状况的动态监控，本论文在高斯烟羽模型的基础上进行改进，引入温度、降雨对污染物扩散的影响，建立了新的污染物扩散模型。本论文创新性的提出了风雨影响指数M,用来衡量风向、降雨对颗粒物扩散的影响。本论文将抽象的污染物含量形象化，利用空气污染指数API描述具体的污染程度及其给周围居民带来的影响。并且从不同角度给出了模型检验，验证了所建模型的准确性。 对于问题（1）具体赔偿方案的制定，在综合考虑了不同方位风向频率、受污染时间、受污染程度的基础上，本论文使用了层次分析法，并且进行了一致性检验，使得赔偿方案具有说服力。通过MATLAB编程，计算出当政府和垃圾处理厂共支付风险赔偿金为N时，得出居住地的每位居民应得的赔偿金额计算公式。对于监测点的设置，经计算共需21个，具体布置情况见后文。 对于问题（2），在题目所述的发生事故的情况下，对污染物的具体含量进行了合理的预测与假设。模拟出酸性物质与颗粒物的影响范围，并根据具体的污染程度设置不同的污染区。对每个污染区的不同情况设置更改监测点的设置，并且在问题（1）的基础上对居民的经济补偿进行合理修改。 关键词：高斯烟羽模型，层次分析法，空气污染指数，烟气抬升公式 一、问题重述 “垃圾围城”是世界性难题，在今天的中国显得尤为突出。数据显示，目前全国三分之二以上的城市面临“垃圾围城”问题，垃圾堆放累计侵占土地75万亩。因此，垃圾焚烧正逐步成为中国垃圾处理的主要手段之一。然而，由于政府监管不力、投资者目光短浅等多方面的原因，致使前些年各地建设的垃圾焚烧电厂在运营中出现了环境污染问题，给垃圾焚烧技术在我国的推广造成了很大阻力，许多城市的新建垃圾焚烧厂选址都出现因居民反对而难以落地的局面。在垃圾焚烧厂运行监管方面，目前主要是在垃圾焚烧厂内进行测量监控，缺少从周边环境视角出发的外围动态监控，因而难以形成为民众所信服的全方位垃圾焚烧厂环境监控体系。 深圳市某地点计划建立一个中型的垃圾焚烧厂，计划处理垃圾量1950吨/天（设置三台可处理垃圾650吨/天的焚烧炉，排烟口高度80米，每天24小时运转）。从构建环境动态监控体系、并根据潜在污染风险对周围居民进行合理经济补偿的需求出发，有关部门希望能综合考虑垃圾焚烧厂对周围带来环境污染以及其他危害的多种因素（例如，焚烧炉的污染物排放量、居住点离开垃圾焚烧厂的距离、风力和风向及降雨等气象条件、地形地貌以及建筑物的遮挡程度等等），在进行科学定量分析的基础
放射性气体扩散浓度预估模型 【摘要】本文是以日本地震引起的福岛核电站的核泄漏为背景，并以给出的数据为基础，研究某一假设核电站的核泄漏问题。我们通过收集相关的资料，并结合题目给出的数据，建立了高斯模型、连续点源高斯扩散模型解决了题目提出的四个问题。 针对问题一：考虑到泄漏源是连续、均匀和稳定的，我们运用散度、梯度、流量等数学概念，通过“泄漏放射性物质质量守恒”、“气体泄漏连续性定理”、 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 云团扩散模型 根据物质泄漏后所形成的气云的物理性质的不同，可以将描述气云扩散的模型分为非重气云模型和重气云模型两种[5-13]。 1.1 非重气云模型 高斯模型是一种常用的非重气扩散模型，高斯烟羽(Plume model)模型又称高架点连续点源扩散模型，适用于连续源的扩散，即连续源或泄放时间大于或等于扩散时间的扩散。 高斯烟团(Puff model)模型适用于短时间泄漏的扩散，即泄放时间相对于扩散时间比较短的情形，如突发性泄放等。若假设气体云内空间上的分布为高斯分布，则地面地处风向的烟团浓度分布算式为 式中， c(x，y，H)——点(x，y，H)处浓度值，mg/m3； Q——源强，即单位时问的排放量，mg/s； u——环境平均风速，m/s； σx，σy，σz——扩散参数； H——源高(烟团高度)，m； x——下方向到泄漏原点的距离，m； y，z——侧风方向、垂直向上方向离泄漏原点的距离，m。 高斯模式的实际应用效果很大程度上依赖于如何给定模式中的一些参数，尤其要注意源强、扩散参数等的确定。 源强与污染物的物理化学属性、扩散方式、释放点的地理环境等有关。扩散参数表征大气边界层内
湍流扩散的强弱，是高斯模式的一项重要数据。高斯扩散模式所描述的扩散过程(实质上也包含了在实际应用中对高斯模式的一些限制)主要有： 1)下垫面平坦、开阔、性质均匀，平均流场稳定，不考虑风场的切变。 2)扩散过程中，污染物本身是被动、保守的，即污染物和空气无相对运动，且扩散过程中污染物无损失、无转化，污染物在地面被反射。 3)扩散在同一温度层结中发生，平均风速大于1.0 m/s。 4)适用范围一般小于10～20 km。 1.2 重气云模型 由于重气本身的特殊性，在重气扩散领域也有大量基于不同理论的模型。鉴于重气扩散与中性或浮性气体扩散有着明显的区别，目前国内外已开发大量的不同复杂程度的重气扩散模型，如箱模型、相似模型、LTA－HGDM模型、CFD模型等。 1.2.1 箱(BOX)模型 箱模型是指假定浓度、温度和其他场，在任何下风横截面处为矩形分布等简单形状，这里的矩形分布是指在某些空间范围内场是均匀的，而在其他地方为零。该类模型预报气云的总体特征，如平均半径、平均高度和平均气云温度，而不考虑其在空间上的细节特征。重气效应消失后其行为表现为被动气体扩散，所以该类模型还包括被动扩散的高斯模型及对它的修正。 1.2.2 层流及湍流大气环境中的重气扩散(LTA－HGDM)模型 LTA－HGDM模型(Heavy Gas Dispersion Model in Lsaminar and Turbulent Atmosphere层流及湍流大气环境中的重气扩散模型)以箱模型为基础，结合虚点源模型，能描述重气泄漏扩散整个过程。模型同三维有限元模型相比，具有形式简单、原始输入数据运算速度快等优点。 LTA－HGDM模型的建立基于以下几点假设： 1)危险性气体初时泄漏时，其外形呈正圆柱形(H=2R)。 2)初始时刻泄漏源即此核电站内部的浓度、温度呈均匀分布。 3)扩散过程不考虑泄漏源即此核电站内部温度的变化，忽略热传递、热对流及热辐射。
大气污染物扩散的高斯模型模拟：可视化模拟点源大气污染的扩散Gaussian Atmospheric Dispersion Model 突发性大气污染事故时有发生，对大气污染扩散进行模拟和分析，有利于减小事故的危害，减轻人员伤亡和财产损失。高斯扩散模型是国际原子能机构(IAEA)推荐使用于重气云扩散模拟的数学模型,该模型在非重气云扩散的应用日益广泛。高斯扩散模型是描述大气对有害气体的输移、扩散和稀释作用的物理或数学模型,是进行灾害预测和救援指挥的有力手段之一。 高斯扩散模型 高斯模型又分为高斯烟团模型和高斯烟羽模型。大气污染物泄漏分为瞬时泄漏和连续泄漏，瞬时泄漏是指污染物泄放的时间相对于污染物扩散的时间较短如突发泄漏等的情形，连续泄漏则是指污染物泄放的时间较长的情形。瞬时泄漏采用高斯烟团模型模拟，而连续泄漏采用高斯模型烟羽模型模拟。高斯模型适用于非重气云气体，包括轻气云和中性气云气体。要求气体在扩散过程中，风速均匀稳定。 在高斯烟团模型中，选择风向建立坐标系统，即取泄漏源为坐标原点，x轴指向风向，y轴表示在水平面内与风向垂直的方向，z轴则指向与水平面垂直的方向，具体公式见式: (mg/s); x、y、z轴上的扩散系数，需根据大气稳定度选择参数计算得到(m)；x、y、z表示x、y、z上的坐标值(m)；u 表示平均风速(m/s)；t表示扩散时间(s)；H 表示泄漏源的高度(m)。 同理，高斯烟羽模型的表达式如： 技术方法 若用高斯模型算出空间每一个点在一个时刻的污染浓度，这个计算量是很大的。因此所设计的系统一般都是采用先进行图层网格化，由高斯模型计算出有限个网格点的上的污染物浓度，在进行空间内插得到面上每一个点的污染物浓度，并由此得到污染物浓度的等值线。整个过程的示意图如图所示
2 扩散模型 2.1 高斯模型 燃气泄漏后会在泄漏源附近形成气团，气团在大气中的扩散计算通常采用高斯模型。高斯模型的基本形式是在如下的假设条件下推导出来的[1、9]：假定燃气在扩散的过程中没有沉降、化合、分解及地面吸收的发生；燃气连续均匀地排放；扩散空间的风速、大气稳定度都均匀、稳定；在水平和垂直方向上都服从正态分布。 泄漏燃气相对密度小于或接近1的连续泄漏采用高斯烟羽模型。以泄漏点为原点，风向方向为x轴的空间坐标系中的某一点(x，y，z)处的质量浓度计算公式如下[9]： 平均风速＞1m/s时： 平均风速=0.5～1m/s时： 平均风速＜0.5m/s时，假设气团围绕泄漏点浓度均匀分布，则距离泄漏点r 处的燃气质量浓度为： 式中ρ d (x，y，z)——扩散燃气在点(x，y，z)处的质量浓度，kg/m3 x、y、z——x、y、z方向上距泄漏点的距离，m u a ——平均风速，m/s δ x 、δ y 、δ z ——x、y、z方向的扩散系数，m
h——泄漏点高度，m ρ (r)——距离泄漏点r处的燃气质量浓度，kg/m3 d r——空间内任意一点到泄漏点的距离，m a、b——扩散系数，m t——静风持续时间，s，取3600的整数倍 扩散系数可查HJ/T 2.2—93《环境影响评价技术导则大气环境》得到。2.2 重气扩散模型 液化石油气密度比空气密度大，属于重气。该类气体泄漏时在重力的作用下会下沉，这时使用高斯模型计算的结果会使泄漏燃气扩散速度偏大，泄漏源附近的浓度偏小。为了解决这个问题，可以引入最早由Van Ulden提出，并由Manju Mohan等发展的箱式模型[1]。箱式模型分为两个阶段：泄漏后的重气扩散阶段和重气效应消失后的被动气体扩散阶段。 重气泄漏后首先是重气扩散阶段。在这个阶段，重气云团由于重力作用逐渐下沉并不断卷吸周围的空气，在卷吸空气的同时，气云受热，最终当重气云团与空气的密度差＜0.001kg/m3时，可认为气云转变成中性状态。 随着重气的继续扩散，气云所受的重力不再是影响扩散的主要因素，而大气湍流扩散逐渐占主要地位，这时便是被动气体扩散阶段，可以应用高斯模型计算泄漏燃气的扩散。 3 结论 使用泄漏模型可以计算出燃气泄漏的理论量，此量为扩散计算提供基础数据，可以依据此量分析泄漏后的扩散范围以及预测评价事故后果。使用扩散模型可以对燃气泄漏后的危险区域进行预测。泄漏模型和扩散模型都有各自的适用条件和范围，应该根据泄漏扩散的具体情况分析选择相应模型。
放射性气体扩散的预估模型 摘要：由于放射性气体泄漏造成惨重损失的报道在国际屡见不鲜，近日日本福岛核电站的放射性气体的泄漏事件更让我们关注放射性气体泄漏时在环境中的浓度问题，为了今后事故发生后提供积极的补救措施, 所以对放射性气体的扩散作深入的研究是很有必要的。本文结合高斯烟羽模型、线性拟合，以及微分方程模型，运用MA TLAB软件，分析了泄漏源强度、风速、大气稳定度参数、地面粗糙度参数和计算精确度等的因素对放射性气体扩散的影响，预测了放射性气体浓度在不同时间，不同地区的浓度变化，并且本文模型中的数据可以根据不同的实际情况而加以改变，因而使本文的应用范围大大增加，可以适用于具有较强的应用性。文章首先在第一问中利用MA TLAB软件对数据进行线性拟合，采用微分方程模型得到核电站周边放射性气体在不同地区，不同时间段的浓度变化，得出随着离泄漏源距离的延伸，最终放射性物质的浓度越来越小，趋近于零，即当L趋向无穷是，C(x,y,z,t)趋向于零；当时间趋于无穷时,C（x,y,z,t）也趋于无穷。问题二，问题三中，建立以核电站周边不同地区得距离以及风速为因变量，设置各个主要因素的参考数据，同时，利用高斯烟羽模型对核电站周边地区的浓度进行预测，然后，利用MATLAB软件，将相关数据代入程序，我们得到核电站周边地区的浓度分布的等高曲线。问题四中，通过实际收集数据，集合核电站周边地区的浓度等高曲线，可以直观的看出日本福岛核电站对我国东海岸以及美国西海岸的影响。 一．问题的提出 1.1背景的介绍 目前，核电的发展给国家带来了巨大的经济效益和社会效益，但核电正常运行以及发生泄露时不可避免的会有气载放射性核素排出，这样就给周围的环境产生了一定的影响，因此，正确的测出大气中放射性物质的浓度在环境检测以及安全评估中具有重要意义。 1.2需要解决的问题 的放射性气体以匀速排出，设有一座核电站遇自然灾害发生泄漏，浓度为p 速度为m kg/s，在无风的情况下，匀速在大气中向四周扩散, 速度为s m/s. （1）请你建立一个描述核电站周边不同距离地区、不同时段放射性物质浓度的预测模型。 （2）当风速为k m/s时，给出核电站周边放射性物质浓度的变化情况。 （3）当风速为k m/s时，分别给出上风和下风L公里处，放射性物质浓度的预测模型。
污染源扩散模拟的实现及应用 突发性环境污染事故是一种威胁人类安全和健康、破坏生态环境、危害性大的污染事故。近年来，随着经济的发展，越来越多的突发环境事故爆发，造成严重的环境污染，不仅给国家人民财产造成了损失，同时还严重的危害了周边人民的健康。因此，加强突发性环境污染事故应急监测，研究其处理技术，对污染物的扩散进行预报是环境监测和环境保护领域中一项非常重要的工作。 1 大气扩散模型研究的意义 突发性环境污染事故主要是由于高压容器、储罐、输送管道节门的破裂等诸多原因引起的，它可导致有毒有害气体外泄。其特点是没有固定的排放方式和排放途径，突然发生、来势凶猛，在短时间内排放大量有毒有害的污染物，有毒气体外泄后，随大气弥散，中心位置浓度最高，向外逐步扩散稀释，下风向形成相应的时空浓度分布。 对于重大突发事故分析，目前国内外普遍采用仿真技术，通过建立数学模型进行分析，而且已建立了很多适用于不同条件的数学模型。当前应用较为广泛的应急大气扩散模型有：SLAB，DEGADIS，ALOHA，ARCHIE，DEMRA和LPDM,其中比较著名的有RADM、ADOM、STEMD等模型、美国Sigma公司于八十年代中期建立的HPDM模型以及英国剑桥研究院开发的ADMS模型。这些模型通过对早期的CRSTER的法规式模式做了一些改进和发展，从而产生新一代扩散模型。 这些算法以扩散统计理论为出发点，假设污染物的浓度分布在一定程度上服从高斯分布。模式系统可用于多种排放源（包括点源、面源和体源）的排放，也适用于乡村环境和城市环境、平坦地形和复杂地形、地面污染物排放模拟、区域环境容量计算与总量控制等多种功能。 这些扩散模型的特点是基于专有平台实现，自成系统。并且在这些系统中大多考虑了扩散模型源排放、平流输送、湍流扩散、干沉积、湿沉积、气象化学等众多因素，系统功能庞大。但同时它们基于专有平台，自成系统，所以很难同环保局具体的业务系统进行整合。考虑到天津某区环保监控与应急指挥系统的实际情况和具体需求，我们基于高斯扩散模型，并根据实际情况加以改造，最终在GIS系统中进行了展现。 2基于高斯的实用性大气扩散模型的算法实现 高斯扩散模型是高斯应用湍流统计理论，在大量实验数据资料分析以及正态分布假设的基础上，得到的污染物在大气中扩散的数学模型。经过多年的研究试验，国内外建立了多种高斯扩散模型，包括高斯点源扩散模式、点源封闭式扩散模型、高斯面源(虚拟点源)扩散模式以及多种特殊气象条件和复杂地形条件下的高斯扩散模式。高斯扩散模型是目前运用得最普遍的大气扩散的数学模型。例如UK一ADMS模型的Urban部分中就直接采用的模型是一个三维高斯模型，以高斯分布公式为主计算污染物浓度，在非稳定条件下的垂直扩散使用了倾斜式的高斯模型烟羽扩散。 3、基于Supermap的扩散模型图形化显示 SuperMap GIS是国内具有完全自主知识产权的大型地理信息系统软件平台。包括组件式GIS开发平台、服务式GIS开发平台等应用开发平台，同时具备相关的空间数据生产、加工和管理工具。其中B/S开发的组件SuperMap IS .NET采用面向Internet的分布式计算技术，支持跨区域、跨网络的复杂大型网络应用系统集成，提供可伸缩、多种层次的WebGIS 解决方案，全面满足网络GIS应用系统建设的需要。
数学建模高斯扩散模 型
§4-2高斯扩散模式 ū —平均风速； Q—源强是指污染物排放速率。与空气中污染物质的浓度成正比，它是研究空气污染问题的基础数据。通常： （ⅰ）瞬时点源的源强以一次释放的总量表示； （ⅱ）连续点源以单位时间的释放量表示； （ⅲ）连续线源以单位时间单位长度的排放量表示； （ⅳ）连续面源以单位时间单位面积的排放量表示。 δy—侧向扩散参数，污染物在y方向分布的标准偏差，是距离y的函数，m； δz—竖向扩散参数，污染物在z方向分布的标准偏差，是距离z的函数，m； 未知量—浓度c、待定函数A(x)、待定系数a、b； 式①、②、③、④组成一方程组，四个方程式有四个未知数，故方程式可解。 二、高斯扩散模式 （一）连续点源的扩散 连续点源一般指排放大量污染物的烟囱、放散管、通风口等。排放口安置在地面的称为地面点源，处于高空位置的称为高架点源。 1. 大空间点源扩散 高斯扩散公式的建立有如下假设：①风的平均流场稳定，风速均匀，风向平直；②污染物的浓度在y、z轴方向符合正态分布；③污染物在输送扩散中质量守恒； ④污染源的源强均匀、连续。 图5－9所示为点源的高斯扩散模式示意图。有效源位于坐标原点o处，平均风向与x轴平行，并与x轴正向同向。假设点源在没有任何障碍物的自由空间扩散，不考虑下垫面的存在。大气中的扩散是具有y与z两个坐标方向的二维正态分布，当两坐
标方向的随机变量独立时，分布密度为每个坐标方向的一维正态分布密度函数的乘积。由正态分布的假设条件②，参照正态分布函数的基本形式式（5－15），取μ＝0，则在点源下风向任一点的浓度分布函数为： （5－16）式中 C—空间点（x，y，z）的污染物的浓度，mg/m3； A（x）—待定函数； σy、σz—分别为水平、垂直方向的标准差，即y、x方向的扩散参数，m。 由守恒和连续假设条件③和④，在任一垂直于x轴的烟流截面上有： （5－17） 式中 q—源强，即单位时间内排放的污染物，μg/s； u—平均风速，m/s。 将式（5－16）代入式（5－17）, 由风速稳定假设条件①，A与y、z无关，考虑到③和④，积分可得待定函数A（x）： （5－18） 将式（5－18）代入式（5－16），得大空间连续点源的高斯扩散模式 （5－19） 式中，扩散系数σy、σz与大气稳定度和水平距离x有关，并随x的增大而增加。当y＝0，z＝0时，A（x）＝C（x，0,0），即A（x）为x轴上的浓度，也是垂直于x轴截面上污染物的最大浓度点C max。当x→∞，σy及σz→∞，则C→0，表明污染物以在大气中得以完全扩散。 2．高架点源扩散
SLAB 用户手册:模拟重气体泄漏的空气扩散模型 中文简要用户使用手册 环境保护部环境工程评估中心 国家环境保护环境影响评价数值模拟重点实验室
手册说明 本用户手册基于《USER’S MANUAL FOR SLAB: AN ATMOSPHERIC DISPERSION MODEL FOR DENSER-THAN-AIR RELEASES》(1990.06)编写，仅对美国EPA网站所提供的模拟重气体泄漏的空气扩散模型SLAB的使用方法提供中文版简要说明，更详细的程序使用说明请查阅相关的软件手册及文档，或采用带图形界面版的商业软件。 本手册由环境保护部环境工程评估中心国家环境保护环境影响评价数值模拟重点实验室负责编写，参与人员包括：易爱华、陈陆霞、胡翠娟、梁昊、杨晔、丁峰等。 本手册版权所有，转载及印刷请与环境保护部环境工程评估中心联系。 本手册所涉及的模型系统及本手册电子版本下载地址：
一、SLAB简介 SLAB是用于模拟重气体泄漏的空气扩散模型。该模型最初基于Zeman于1982年提出的关于重气体云的空气卷吸和重力扩散的理念而开发。SLAB早期相关工作由美国能源署支持。SLAB的进一步开发由USAF工程和服务中心（1986年开始）和美国石油学会（1987开始）共同提供支持。 现行的SLAB版本可以模拟连续的、限时的和瞬时的物质泄漏，泄漏源包括以下4种：地面液池蒸发、高于地面的水平射流、烟囱或高于地面的垂直射流，以及瞬时释放。 SLAB除可以用于模拟重气体的扩散，还可以模拟中性浮力气体的烟云扩散，以及烟云轻于空气时的上升过程。 泄漏时的空气扩散过程可以通过求解质量、动量、能力和物质的守恒方程来计算，如图1所示。为了简化守恒方程的求解过程，方程可以通过将烟云作为稳态烟羽或瞬时烟团在空间上进行平均。连续排放（持续时间非常长的排放源）可以作为稳态烟羽。有限时间的排放采用稳态烟羽模式描述最初烟云的扩散，而且在该排放源持续泄漏的时间段内，可以一直使用稳态烟羽模式。释放一旦终止，烟团被视为瞬时烟团，之后的扩散采用瞬时烟团模式来计算。对于瞬时泄漏的排放源，整个过程都均使用瞬时烟团扩散模式。 二、理论介绍 2.1重气体扩散模型简介 重气体泄漏的空气扩散模型受到多种物理现象的影响，这些物理现象在中性或浮力气体泄漏中可能不会发生或者即便发生也不是很重要。这些物理现象包括：重气体烟云的稳定密度分层导致的湍流衰减；由于重力流和初始排放源动量导致的环境速率场的改变；由于液滴形成和挥发以及在过热或低温液体排放情况下的地面加热对烟云温度、浮力和湍流的热力学效应；此外，我们所关注的某种特定的重气体的浓度可能和典型大气污染物关注的累积浓度差别很大。例入，对于易燃气体，关注的是瞬间浓度；而对于有毒气体，关注的则是几分钟到几小时的浓度，以及累积浓度。因此，为了能够更好的预测出重气体泄漏时有毒浓度区的大小和持续时间，所有重要的物力现象都需要进行考虑，而且预测过程中要使用最合理的浓度平均时间。 为了满足重气体泄漏情形的要求，SLAB模型以质量、动量、能量和物质守恒方程的平均形式为起点，在该理论框架的基础上进行开发（如图1所示）。这些方程用于计算扩散烟团的空间平均性质，并且以两种方式来代表两种不同的扩散模式：稳态烟羽扩散模式和瞬时烟团扩散模式。