城市污染物扩散模型

大气污染物扩散模式

第四章大气扩散浓度估算模式第一节湍流扩散的基本理论一湍流 1．定义：大气的无规则运动风速的脉动风向的摆动 2．类型：按形成原因热力湍流：温度垂直分布不均（不稳定）引起，取决于大气稳定度机械湍流：垂直方向风速分布不均匀及地面粗糙度引起 3．扩散的要素风：平流输送为主，风大则湍流大湍流：扩散比分子扩散快105～106倍二湍流扩散理论（主要阐述湍流与烟流传播及湍流与物质浓度衰减的关系） 1.梯度输送理论通过与菲克扩散理论类比建立起来的（菲克定律：单位时间内通过单位断面上的物质的数量与浓度梯度呈正比）类比于分子扩散，污染物的扩散速率与负浓度梯度成正比 x C k F ??-= 式中，F — 污染物的输送通量 k — 湍流扩散系数 C — 污染物的浓度 X — 与扩散截面垂直的空间坐标（扩散过程的长度） x C ??— 浓度梯度要求得各种条件下某污染物的时、空分布，由于边界条件往往很复杂，不能求出严格的分析解，只能是在特定的条件下求出近似解，再根据实际情况进行修正。 2.湍流统计理论泰勒首先将统计理论应用在湍流扩散上图4-1显示：从原点O 放出的粒子，在风沿着x 方向吹的湍流大气中扩散。粒子的位置用y 表示，则结论为： ①y 随时间变化，但其变化的平均值为零 ②若从原点放出很多粒子，则在x 轴上粒子的浓度最高，浓席分布以x 轴为对称轴，并符合正态分布。萨顿实用模式：解决污染物在大气中扩散的实用模式高斯模式：应用湍流统计理论得出正态分布假设下的扩散模式 3.相似理论第二节高斯扩散模式一坐标系的建立—右手坐标系

1．原点O ：无界点源或地面源，O 为污染物的排放点高架源，O 为污染物的排放点在地面上的投影点补充：点源高架源连续源固定源线源地面源间歇源流动源面源 2．x 轴：正向为平均风向，烟流中心线与x 轴重合 3．y 轴：垂直于x 轴 4．z 轴：垂直于xoy 平面二高斯模式的有关假定 1．污染物浓度在y 、z 轴上的分布为正态分布； )2exp(21 )(22 y y y y f σπ σ-= )2exp(21 )(22 z z z z f σπ σ-= y σ，z σ— 分别为污染物在y 和z 方向上分布的标准差，m 2．全部高度风速均匀稳定，即风速u 为常数； 3．源强是连续均匀稳定的，源强Q 为定值； 4．扩散中污染物是守恒的，不考虑转化，即烟云在扩散过程中没有沉降、化合、分解及地面吸收、吸附作用发生； 0=??t C 5．在x 方向上，输送作用远远大于扩散作用，即 )(x C k x x C u x ????>>??； 6．地面足够平坦。

水污染模型

基于GIS 的环境污染应急分析系统的开发重点是实现水体污染扩散模拟。目前, 国外在此方面的研究成果很多，已经进行到了三维水体污染扩散模拟，国内的起步则较晚, 至今的研究成果在一维的较多，二维和三维的较少。鉴于目前网络的发展, 有必要将互联网与系统结合起来。一维水体污染扩散数学模型：一维水质模型是水环境模型中相对简单的一种，是河流、河口和湖泊遭受污染时，实际的断面浓度分布与断面浓度的平均值偏差不大时常采用的水污染预测模型。它主要研究污染物浓度分布沿程的变化以及各个断面上污染物浓度随时间的变化，其中河流以一维水质模型最为常见。在突发性河道水源地污染事故发生时。污染物的排放存在两种情况，即一维稳定排放和一维瞬时排放, 二维水体污染扩散数学模型：二维计算模型模拟速度快、实时而精度无需很高, 可忽略基本控制方程中的一些非主要因素,模型结构简单、实用性强。目前最为常用的有限差分数值计算方法对控制方程进行离散, 按物理分步法将二维偏微分方程化简成较简单的一维方程, 应用广为采用的ADI隐式格式联合求解水动力模型与水污染模型。算法具有编程简单、占用计算机内存较小、无条件稳定、可适当增大空间步长、计算效率高、易于实现自动化的实时模拟计算等显著优点, 适合于在应急处置中应用。并且利用GIS 的强大的空间分析、处理和表现功能, 将水力计算与GIS 结合在一起, 实现了污染模拟结果的二维可视化, 为应急处置提供一个形象、直观的表现平台, 能有效地辅助应急决策。三维水体污染扩散数学模型：水污染三维可视化包含两方面的内容：河道地形地貌三维仿真与污染扩散可视化，二者通过地理坐标进行空间叠加形成河道污染扩散可视化展示平台，在此基础上进行各种统计分析功能。

点污染源空气污染扩散模型

8 点、中午12 点、晚上9 点都没有排放气体，该怎么算，是不是需要找到一个关于时间t的函数，来计算多长时间之后污染还剩下多少 c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.^2)./((sigy+eps).^2)).*(exp(-0.5*(z-H).^2./((sigz+eps).^2))+exp(-0.5*(z+H).^2./((sigz+ eps).^2))); 这个函数对吗？该调用什么函数？问题：建立单污染源空气污染扩散模型，描述其对周围空气污染的动态影响规律。现有河北境内某一工厂废气排放烟囱高50m，主要排放物为氮氧化物。早上9 点至下午 3 点期间的排放浓度为406.92mg/m3，排放速度为1200m3 /h；晚上10 点-凌晨4 点期间的排放浓度为1160mg/m3，排放速度为5700m3 /h；通过你的扩散模型求解该工厂方圆51 公里分别在早上浓度8 点、中午12 点、晚上9 点空气污染分布和空气质量等级。源代码 clear all clc [x,y]=meshgrid(0:20:5100,0:20:5100); Q=135.64; z=1.5; H=50; u=1.94; sigy=0.3914238*x.^0.865014; sigz=0.0757182*x.^1.00770; c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.^2)./((sigy+eps).^2)).*(exp(-0.5*(z-H).^2./((sigz+eps).^2))+exp(-0.5*(z+H).^2./((sigz+ eps).^2))); mesh(x,y,c); xlabel('X'),ylabel('Y'),zlabel('C'), clear all clc [x,y]=meshgrid(-5100:20:5100,-5100:20:5100); Q=1836.7; z=1.5; H=50; u=1.7; sigy=0.3914238*x.^0.865014; sigz=0.0757182*x.^1.00770; c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.^2)./((sigy+eps).^2)).*(exp(-0.5*(z-H).^2./((sigz+eps).^2))+exp(-0.5*(z+H).^2./((sigz+ eps).^2))); mesh(x,y,c); xlabel('X'),ylabel('Y'),zlabel('C'), 分享到： 2015-05-29 16:32 提问者采纳 clear all [x,y]=meshgrid(-51000:100:51000,-51000:100:51000); Q=135.64; z=1.5; H=50; u=1.94; sigy=0.3914238*x.^0.865014;

大气污染扩散模型

第一节大气污染物的扩散一、湍流与湍流扩散理论 1. 湍流低层大气中的风向是不断地变化，上下左右出现摆动；同时，风速也是时强时弱，形成迅速的阵风起伏。风的这种强度与方向随时间不规则的变化形成的空气运动称为大气湍流。湍流运动是由无数结构紧密的流体微团——湍涡组成，其特征量的时间与空间分布都具有随机性，但它们的统计平均值仍然遵循一定的规律。大气湍流的流动特征尺度一般取离地面的高度，比流体在管道内流动时要大得多，湍涡的大小及其发展基本不受空间的限制，因此在较小的平均风速下就能有很高的雷诺数，从而达到湍流状态。所以近地层的大气始终处于湍流状态，尤其在大气边界层内，气流受下垫面影响，湍流运动更为剧烈。大气湍流造成流场各部分强烈混合，能使局部的污染气体或微粒迅速扩散。烟团在大气的湍流混合作用下，由湍涡不断把烟气推向周围空气中，同时又将周围的空气卷入烟团，从而形成烟气的快速扩散稀释过程。烟气在大气中的扩散特征取决于是否存在湍流以及湍涡的尺度(直径)，如图5－7所示。图5－7（a）为无湍流时，烟团仅仅依靠分子扩散使烟团长大，烟团的扩散速率非常缓慢，其扩散速率比湍流扩散小5～6个数量级；图5 －7（b）为烟团在远小于其尺度的湍涡中扩散，由于烟团边缘受到小湍涡的扰动，逐渐与周边空气混合而缓慢膨胀，浓度逐渐降低，烟流几乎呈直线向下风运动；图5－7（c）为烟团在与其尺度接近的湍涡中扩散，在湍涡的切入卷出作用下烟团被迅速撕裂，大幅度变形，横截面快速膨胀，因而扩散较快，烟流呈小摆幅曲线向下风运动；图5－7（d）为烟团在远大于其尺度的湍涡中扩散，烟团受大湍涡的卷吸扰动影响较弱，其本身膨胀有限，烟团在大湍涡的夹带下作较大摆幅的蛇形曲线运动。实际上烟云的扩散过程通常不是仅由上述单一情况所完成，因为大气中同时并存的湍涡具有各种不同的尺度。根据湍流的形成与发展趋势，大气湍流可分为机械湍流和热力湍流两种形式。机械湍流是因地面的摩擦力使风在垂直方向产生速度梯度，或者由于地面障碍物(如山丘、树木与建筑物等)导致风向与风速的突然改变而造成的。热力湍流主要是由于地表受热不均匀，或因大气温度层结不稳定，在垂直方向产生温度梯度而造成的。一般近地面的大气湍流总是机械湍流和热力湍流的共同作用，其发展、结构特征及强弱决定于风速的大小、地面障碍物形成的粗糙度和低层大气的温度层结状况。 2. 湍流扩散与正态分布的基本理论气体污染物进入大气后，一面随大气整体飘移，同时由于湍流混合，使污染物从高浓度区向低浓度区扩散稀释，其扩散程度取决于大气湍流的强度。大气污染的形成及其危害程度在于有害物质的浓度及其持续时间，大气扩散理论就是用数理方法来模拟各种大气污染源在

污染物扩散模型-深圳数学建模

赛区评阅编号（由赛区组委会填写）： 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小时运转）。从构建环境动态监控体系、并根据潜在污染风险对周围居民进行合理经济补偿的需求出发，有关部门希望能综合考虑垃圾焚烧厂对周围带来环境污染以及其他危害的多种因素（例如，焚烧炉的污染物排放量、居住点离开垃圾焚烧厂的距离、风力和风向及降雨等气象条件、地形地貌以及建筑物的遮挡程度等等），在进行科学定量分析的基础

大气污染物扩散高斯模型模拟

大气污染物扩散的高斯模型模拟：可视化模拟点源大气污染的扩散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)。同理，高斯烟羽模型的表达式如：技术方法若用高斯模型算出空间每一个点在一个时刻的污染浓度，这个计算量是很大的。因此所设计的系统一般都是采用先进行图层网格化，由高斯模型计算出有限个网格点的上的污染物浓度，在进行空间内插得到面上每一个点的污染物浓度，并由此得到污染物浓度的等值线。整个过程的示意图如图所示

第五章大气污染扩散

第五章大气污染扩散第一节大气结构与气象有效地防止大气污染的途径，除了采用除尘及废气净化装置等各种工程技术手段外，还需充分利用大气的湍流混合作用对污染物的扩散稀释能力，即大气的自净能力。污染物从污染源排放到大气中的扩散过程及其危害程度，主要决定于气象因素，此外还与污染物的特征和排放特性，以及排放区的地形地貌状况有关。下面简要介绍大气结构以及气象条件的一些基本概念。一、大气的结构气象学中的大气是指地球引力作用下包围地球的空气层，其最外层的界限难以确定。通常把自地面至1200 km左右范围内的空气层称做大气圈或大气层，而空气总质量的98．2％集中在距离地球表面30 km以下。超过1200 km的范围，由于空气极其稀薄，一般视为宇宙空间。自然状态的大气由多种气体的混合物、水蒸气和悬浮微粒组成。其中，纯净干空气中的氧气、氮气和氩气三种主要成分的总和占空气体积的99.97％，它们之间的比例从地面直到90km高空基本不变，为大气的恒定的组分；二氧化碳由于燃料燃烧和动物的呼吸，陆地的含量比海上多，臭氧主要集中在55～60km高空，水蒸气含量在4％以下，在极地或沙漠区的体积分数接近于零，这些为大气的可变的组分；而来源于人类社会生产和火山爆发、森林火灾、海啸、地震等暂时性的灾害排放的煤烟、粉尘、氯化氢、硫化氢、硫氧化物、氮氧化物、碳氧化物为大气的不定的组分。大气的结构是指垂直（即竖直）方向上大气的密度、温度及其组成的分布状况。根据大气温度在垂直方向上的分布规律，可将大气划分为四层：对流层、平流层、中间层和暖层，如图5－1所示。 1. 对流层对流层是大气圈最靠近地面的一层，集中了大气质量的75％和几乎全部的水蒸气、微尘杂质。受太阳辐射与大气环流的影响，对流层中空气的湍流运动和垂直方向混合比较强烈，主要的天气现象云雨风雪等都发生在这一层，有可能形成污染物易于扩散的气象条件，也可能生成对环境产生有危害的逆温气象条件。因此，该层对大气污染物的扩散、输送和转化影响最大。大气对流层的厚度不恒定，随地球纬度增高而降低，且与季节的变化有关，赤道附近约为15km,中纬度地区约为10～12 km，两极地区约为8km；同一地区，夏季比冬季厚。一般情况下，对流层中的气温沿垂直高度自下而上递减，约每升高100m平均降低0．65℃。从地面向上至1～1.5 km高度范围内的对流层称为大气边界层，该层空气流动受地表影响

重金属污染

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人(打印并签名)：日期： 2012 年 8 月15 日赛区评阅编号（由赛区组委会评阅前进行编号）：

2010高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：

城市表层土壤重金属污染分析摘要本文运用Matlab、SPSS及Excel软件，综合分析该城区各区重金属污染程度及其重金属污染的主要原因，并根据重金属的传播特征建立数学模型，寻找出污染源，同时更进一步修改模型以更好研究城市地质环境的演变模式。针对问题一，利用Matlab作出各重金属浓度空间分布图，可看出重金属浓度大多呈高浓度向外逐步扩散现象。同时，利用单因素污染指数公式及污染等级规则，通过模糊综合评价模型，得出各类区域重金属污染程度如下：针对问题二，利用SPSS软件对各重金属变量进行R型聚类分析，将重金属污染分为四类：Cr、Ni、Cu为一类，Cd、Pb、Zn为一类，而As、Hg各为一类。并计算出各区域各单金属污染指数进行对照分析，说明该城区重金属污染主要来源于工业活动及其生产品和交通主干道污染。针对问题三，我们根据重金属物质在土壤中扩散的特点，利用有衰减的扩散模型，构造出浓度与坐标的关系，并利用matlab编程，多次试验，合理筛选数据，通过回归分析的方法，求解参数，从而得到不同重金属的不同污染源。针对问题四，由于问题三求解的局限性，我们在模型三的基础上，增加了时间、水文、人类活动等因素，得出一个更为具体并符合实际的模型。有利于计算土壤中重金属的浓度。关键词单因素污染指数模糊综合评价 R型聚类分析有衰减的扩散模型扩散规律

城市表层土壤重金属污染分析数学建模国家二等奖

承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则. 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人(打印并签名)：日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：

编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：

城市表层土壤重金属污染分析摘要本文以某城区为例，对其土壤地质环境重金属含量进行取样分析，找出取样点被污染的主要原因并寻找到取样点的污染源。为此，对于第（1）和第（2）小题我们对取样点的数据统计分析，用MATLAB 软件画出8张元素的分布图，各代表同一重金属元素在不同功能区不同含量范围内的分布点，统计出点的数目，对均值和方差进行对比得出该城区内不同区域重金属的污染程度和主要污染原因。针对第（3）小题，我们建立了2个模型来研究，干尘扩散模型——模型一，高斯水流扩散模型——模型二。模型一是为了求解重金属元素附着在干尘粒上的扩散方式和范围，运用烟雾扩散模型，结论是根据取样点的某金属含量程度来判断该点距污染源的距离，则污染源在以该距离为半径的圆的圆周上，在根据取样点的金属含量高低分布情况，找到污染源。模型二是研究重金属在土壤中的扩散情况，在相间物质交换为平衡的条件下 ,可用阻滞系数来表示其影响，可得到一个微分方程 3t c ρ= ??带入参数，可以找到某金属的污染源坐标。关键词：烟雾扩散模型图例统计污染物运移物质交换 MATLAB

污染空气的扩散模型

放射性气体扩散的预估模型摘要：由于放射性气体泄漏造成惨重损失的报道在国际屡见不鲜，近日日本福岛核电站的放射性气体的泄漏事件更让我们关注放射性气体泄漏时在环境中的浓度问题，为了今后事故发生后提供积极的补救措施, 所以对放射性气体的扩散作深入的研究是很有必要的。本文结合高斯烟羽模型、线性拟合，以及微分方程模型，运用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.模型推导：污染物在全断面混和后，其迁移转化过程可用一维模型来描述，基本控制方程为：S S h A KAC x c E D A x x AUC t AC r x x ++-??+??=??+??])([)()( 其中：C 为污染物质的断面平均浓度，U 为断面平均流速，A 为断面面积，h 为断面平均水深，x D 为湍流扩散系数，K 为污染物降解系数。x E 为纵向扩散系数r S 为河床底泥释放污染物的速率，S 为单位时间内，单位河长上的污染物排放量。实践证明，水的纵向流速是引起污染物浓度变化的主要参数，因此河流各断面的污染物浓度变化主要由这一项引起。因此该模型可以简化。不考虑湍流扩散，河床底泥释放污染物以及沿河其他污染物排放的影响，水污染模型的基本方程为： AKC x C AE x AUC t AC -??=??+??22)()( 2.模型求解：采用有限差分法中的四点隐式差分格式对上式进行数值求解： )(2121121111111j i j i j i j i j i j i j i j i j i C C K x C C C E x C C U t C C -++-++--++-?+-=?-+?- 整理可得 : 其中 2x E a i ?-=；2212K x E t i +?+?=β；2x E i ?-=γ；)2()1(1K x U C x U t C j i j i i -?+?-?=-δ 将上游边界条件带入上式得：将下游边界条件带入，得：从而组成方程组，利用追赶法求解出j i C ; 3:具体实现：本模块通过的含酚污染物污染扩散情况作为实验典型代表来粗略模拟实现扩散过程。系统默认提供河流参数等数据。设置K 为2/d ，U 为流速为10m/s 。x E 为1d km /2。t ?为100s ，x ?为1000m ；根据上述参数计算出方程组的参数。定义二维数组M[i,j]表示在i 断面j 时刻的浓度。通过距离量算来确定排污口与测量点的距离，计算测量点的浓度，并得到污染物在河道断面上的扩散规律。

重金属污染物扩散FICK 模型

FICK 模型： <一>，Fick 第二定律简介： 22x C D T C ??=?? C T C =),0( 0)0,(C x C = 0),(C T C =∞ 符号说明： ),(T x C : T 时间后距离扩散源x 米处气体浓度。 0C : 初始浓度。 C : 气体浓度。 D : 扩散系数。 <二>，模型假设：（1）假设扩散系数为常数。（2）扩散过程中没有特殊外界条件影响（刮风或下雨等）。（3）扩散过程中外界温度恒定。（4）扩散过程中气体没有与外界发生化学反应。 <三>，模型求解：假设气体发生泄漏（T=0）时刻，所以气体泄漏中心周围处浓度C 与时间和距离（r ）的关系可用Fick 第二定律求解，公式如下： 22r C D T C ??=?? （1）令T r =λ代入（1）式可得： T C D r C D 1.222λ??=?? (2) 可得出(1)式为： λλλ .222??=??-C C D (3) 若n=2，D 41=α，代入化简，积分可得：

?+-=λλλ02)4exp(B d D A C (4) 令DT r D 22==λ β，由高斯误差积分公式可求得： ? ∞=-022)exp(πββd （5）根据初始化条件可知，02C A π- =,0C B = 最终求得的模型公式为： )2()exp(2000200Dt r erf C C d C C C -=--=?β ββπ （6）不妨取D=0.00001,发生泄漏时初始气体浓度1000=C 量纲，扩散时间T=1.0*10 8秒，利用Matlab 编程模拟，得到模拟图像如下：模型评价： Fick 第二定律模拟出了理想状态下气体源扩散的过程，从仿真图（地平面点浓度）可以直观的看出，在气体源位置浓度最大，然后均匀向四周浓度逐步减少扩散，大致符合实际气体扩散过程，模型缺点为，条件过于苛刻，现实气体扩散中，往往收到外界条件影响，比如风向等等。

大气模型

大气模型发展简史与简介 By Laiwf | Published: 2010年04月16日 1.1 第一代空气质量模型―高斯模型和拉格朗日烟团轨迹模型第一代空气质量模型主要包括了高斯扩散模型和拉格朗日轨迹模型。这两类模型都是利用风的运动轨迹来模拟近地层大气层中复杂的物理和化学过程。它的物理表述即模拟均匀混合的大气物质沿风向运动的情况。在大气物质从地面向高层运动的过程中，其运动规则受到垂直方向上风速以及温度的不均匀分布的影响而不断的发生变化。具体过程见图。 1 EIAA （典型高斯）适用于<50km的区域 EIAA大气环评助手“是宁波环科院六五软件工作室开发的软件。《 HJ/T2.2-93 环评导则–大气环境》、《JTJ005-96 公路建设项目环评规范-大气部分》，中国环境影响评价培训教材等文献中推荐的模型和计算方法作为主要框架，内容涵盖了导则中的全部要求，并进行了适当地拓展与加深。可以处理点源、面源、体源、线源对于预测计算结果，可以查看 §各接受点地面高程及其等高线图 §各接受点的背景浓度及其分布图 §各污染源的浓度和总的浓度及其分布图 §各污染源的分担率及其分布图 §各污染源或总的浓度的平均评价指数和超标面积

§还可以任意改变各污染源的排放率（排放强度）以观察不同排放率下的浓度变化情况 §也可查看任意一个横截面或竖截面上的浓度变化图广泛应用的版本是EIAA2.5，EIAA2.6。版本中均有bug，大家谨慎使用。 2 aermod（稳态高斯）适用于<50km的区域 AERMOD由美国国家环保局联合美国气象学会组建法规模式改善委员会（AERMIC）开发。 AERMIC的目标是开发一个能完全替代ISC3的法规模型，新的法规模型将采用ISC3的输入与输出结构、应用最新的扩散理论和计算机技术更新ISC3 计算机程序、必须保证能够模拟目前ISC3能模拟的大气过程与排放源。 20世纪90年代中后期，法规模式改善委员会在美国国家环保局的财政支持下，成功开发出AERMOD扩散模型。该系统以扩散统计理论为出发点，假设污染物的浓度分布在一定程度上服从高斯分布。模式系统可用于多种排放源（包括点源、面源和体源、线源）的排放，也适用于乡村环境和城市环境、平坦地形和复杂地形、地面源和高架源等多种排放扩散情形的模拟和预测。 Aermod作为正在开发的模型，模型中还存在bug，但其不在改进，而且模型在小范围预测的精准性是其他模型不能比拟的，模型需要地面气象数据、高空气象数据和地形数据（平坦地形不需要） 3 CALPUFF（不稳态高斯）对于比较大的区域

污染物迁移分析模型POLLUTEv7原理与应用

污染物迁移分析模型POLLUTEv7原理与应用摘要：垃圾填埋场渗沥液中的污染物严重威胁填埋场底部的土体和地下水环境。为了对垃圾渗沥液中各种无机和有机污染物的迁移规律进行研究，需要采用分析功能强大的计算模型进行模拟。污染物迁移模拟模型POLLUTEv7可以用来提供快速、准确、全面的污染物运移分析能力。模型能够根据污染物迁移边界条件建立一维无限空间对流-弥散方程。采用拉普拉斯变换和拉普拉斯逆变换求解控制方程。可以考虑的范围从在天然粘土隔水简单系统到垃圾填埋场的设计模型复合衬垫，多重障碍和多含水层。除了对流-弥散运输， POLLUTEv7可以考虑非线性吸附、放射性衰变和生物降解衰变，运输通过各类垃圾填埋场防渗衬垫，随时间变化的特性，和相位的变化。本文对污染物迁移模拟模型POLLUTEv7原理进行分析，并介绍了其在工程中的应用。关键词：填埋场，污染物，迁移，原理 1污染物迁移分析原理城市生活垃圾填埋场产生的渗沥液将对填埋场周边的土壤和地下水造成严重的污染[1]。为了控制渗沥液对周边环境的污染，需要对填埋场渗沥液中污染物的迁移规律进行准确分析。POLLUTEv7[2]可以用来提供快速、准确、全面的污染物运移分析能力。其基本分析原理如下： 1.1 控制方程 POLLUTEv7采用一维对流-扩散模型模拟污染物沿竖直方向通过防渗衬里，模型的基本假设为：（a ）污染源位于土层顶部，并且污染物在土层中的扩散是一维的，不考虑其他外部污染源；（b ）每层土的物理特性（如扩散系数、流速、孔隙率等）为均匀分布；（c ）土层中土壤为饱和状态，水流为平流，不考虑瞬时流动；（d ）土壤颗粒对污染物的的吸附为线性吸附。在上述基本假设前提下，污染物在土层中一维对流-扩散模型的控制方程为： 22d C C C R D v C t z z λ???=--??? (1) 式（1）中C 为污染物的浓度，D 为污染物通过土层的扩散系数，v 为达西速度，λ为衰变系数， λ=ln2/T 1/2，T 1/2为污染物在土层中的半衰期。R d 为阻滞因子，其表达式为： 1d d K R n ρ=+ (2) 式（2）中，ρ为土层的干密度，K d 为，n 为土层的孔隙率。 1.1 方程求解采用Rowe [3]提出的有限层方法求解控制方程（1），该方法为半解析半数值的方法，其主要步骤为：（a ）采用Laplace 变换简化控制方程（1），然后采用解析方法求解变换后的偏微分方程；（b ）采用数值方法将上一步骤中得到的解析解求Laplace 逆变换，得到控制方程（1）的数值解。 2污染物迁移分析应用 2.1分析模型选取美国规范RCRA 的典型垃圾填埋场。该填埋场由的复合衬垫和主渗滤液收集系统。该复合衬垫是由一个 1.5毫米土工膜与0.9米厚压实粘土衬垫构成，土工膜与压实粘土之间接触良好。根据工程实际情况，假定土工膜上缺陷的面积为0.1平方厘米，缺陷的数量为1个缺陷每英亩（2.5公顷）。采用吉罗等人（1992 ）提出的方法计算污染物通

点污染源空气污染扩散模型

8点、中午12点、晚上9点都没有排放气体，该怎么算，是不是需要找到一个关于时间t的函数，来计算多长时间之后污染还剩下多少 c=Q./(2*pi*sigy.*sigz*u+eps)*exp(-0.5*(y.A2)./((sigy+eps).A2)).*(exp(-0.5*(z-H).A2./((sigz+eps).A2))+exp(-0.5*(z+H).A2./((sigz+ eps)A2))); 这个函数对吗？该调用什么函数？问题：建立单污染源空气污染扩散模型，描述其对周围空气污染的动态影响规律。现有河北境内某一工厂废气排放烟囱高50m，主要排放物为氮氧化物。早上9点至下午 3点期间的排放浓度为406.92mg/m3，排放速度为1200m3 /h;晚上10点-凌晨4点期间的排放浓度为1160mg/m3，排放速度为5700m3 /h;通过你的扩散模型求解该工厂方圆51 公里分别在早上浓度8点、中午12点、晚上9点空气污染分布和空气质量等级。源代码 clear all clc [x,y]=meshgrid(0:20:5100,0:20:5100); Q=135.64; z=1.5; H=50; u=1.94; sigy=0.3914238*x.A0.865014; sigz=0.0757182*x.A1.00770; c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.A2)./((sigy+eps).A2)).*(exp(-0.5*(z-H).A2./((sigz+eps).A2))+exp(-0.5*(z+H).A2./((sigz+ eps)A2))); mesh(x,y,c); xlabel('X'),ylabel('Y'),zlabel(C), clear all clc [x,y]=meshgrid(-5100:20:5100,-5100:20:5100); Q=1836.7; z=1.5; H=50; u=1.7; sigy=0.3914238*x.A0.865014; sigz=0.0757182*x.A1.00770; c=Q./(2*pi*sigy.*sigz*u+eps).*exp(-0.5*(y.A2)./((sigy+eps).A2)).*(exp(-0.5*(z-H).A2./((sigz+eps).A2))+exp(-0.5*(z+H).A2./((sigz+ eps)A2))); mesh(x,y,c); xlabel('X'),ylabel('Y'),zlabel(C), 分享到： 2015-05-29 16:32 提问者采纳 clear all [x,y]=meshgrid(-51000:100:51000,-51000:100:51000); Q=135.64; z=1.5; H=50; u=1.94; sigy=0.3914238*x.A0.865014;

城市污染物扩散模型

大气污染物扩散模式

水污染模型

点污染源空气污染扩散模型

大气污染扩散模型

污染物扩散模型-深圳数学建模

大气污染物扩散高斯模型模拟

第五章 大气污染扩散

重金属污染

城市表层土壤重金属污染分析数学建模国家二等奖

污染空气的扩散模型

最新城市表层土壤重金属污染分析

污染物扩散模型

重金属污染物扩散FICK 模型

大气模型

污染物迁移分析模型POLLUTEv7原理与应用

点污染源空气污染扩散模型

第五章大气污染扩散