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

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;

sigz=0.0757182*x".00770;

%c=@(x,y)Q./(2*pi*sigy*sigz*u+eps).*exp(-0.5*(y.A2)./((sigy+eps).A2)).*(exp(-0.5*(z-H).A2./((sigz+eps)42))+exp(-0.5*(z+H).A

2./((sigz+eps)八2)));

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)八2))); %g=dblquad(c,-51000,51000,-51000,51000,100),

%g=abs(g/51000/51000*1000),

mesh(x,y,c);

xlabel('X'),ylabel('Y'),zlabel(C),

clear all

clc

[x,y]=meshgrid(-51000:100:51000,-51000:100:51000);

Q=1836.7; z=1.5; H=50; u=1.7;

sigy=0.3914238*x.P.865014;

sigz=0.0757182*x".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)八2)));

mesh(x,y,c);

xlabel('X'),ylabel('Y'),zlabel(C),

clear all

[x,y]=meshgrid(-51000:100:51000,-51000:100:51000);

Q=1836.7; z=1.5; H=50; u=1.7;

sigy=0.3914238*x.A0.865014;

sigz=0.0757182*x".00770;

c=@(x,y)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),

g=dblquad(c,-51000,51000,-51000,51000,100),

g=abs(g/51000/51000*1000),

其具体标准如下：当空气污染指数达0—50时为1级；51 —100时为2级；101 —200时为3级；201 —300时为4级；300以上时为5级。其中3级属于轻度污染，4级属于中度污染，5级则属于重度污染

按照国家统一规定，空气质量达到优良标准即达到国家质量二级标准是指空气污染指数小于等于100，如果空气污染指数小于

等于50，说明空气质量为优。空气污染指数大于50且小于等于100时，说明空气质量为良好。具体如何判断，下面就给出空

气质量的评定标准。

空气质量的评定标准：

我国空气质量采用了空气污染指数进行评价。空气污染指数是根据环境空气质量标准和各项污染物对人体健康和生态环境的影

响来确定污染指数的分级及相应的污染物浓度值。我国目前采用的空气污染指数(API)分为五个等级：

API值小于等于50，说明空气质量为优，相当于国家空气质量一级标准，符合自然保护区、风景名胜区和其它需要特殊保护地

区的空气质量要求；

API值大于50且小于等于100，表明空气质量良好，相当于达到国家质量二级标准；

API值大于100且小于等于200，表明空气质量为轻度污染，相当于国家空气质量三级标准；

API值大于200表明空气质量差，称之为中度污染，为国家空气质量四级标准；

API大于300表明空气质量极差，已严重污染。

根据我国空气污染特点和污染防治重点，目前计入空气污染指数的项目暂定为：二氧化硫、氮氧化物和总悬浮颗粒物。随着环境保护工作的深入和监测技术水平的提高，将调整增加其它污染项目，以便更为客观地反映污染状况。

大气污染物扩散模式

第四章大气扩散浓度估算模式第一节湍流扩散的基本理论一湍流 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．地面足够平坦。

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

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;

城市空气质量模式研究进展

城市空气质量模式研究进展大气污染物扩散模式的应用受到物理化学特征、污染源特征等多种因素的制约，综述了不同的大气污染物扩散模式适用的范围，分别阐述了ISC、AERMOD、ISCST、ADMS等适用于中小尺度的和CALPUFF等应用于大尺度的不同大气污染扩散模式的研究进展。标签：大气污染；扩散模式；空气质量 1 引言工业污染源的排放对大气环境质量有直接的影响，同时也直接影响了周围居民的身体健康及生活。但由于高成本和相关实验的难度，对污染物浓度进行准确的动态分布监测不是十分可行，因此大气污染物扩散模式被广泛用来模拟预测污染物的扩散分布情况，评估大气环境质量。大气污染物扩散模式结合污染物浓度和气象资料定量分析污染物在大气中输送、扩散特征。最初，模式的研究理论核心是高斯扩散理论，应用范围是小尺度。随着研究的逐渐深入和计算机的发展，开始利用计算机进行数值计算，突破了高斯扩散理论均匀平稳湍流的限制，可以求解非均匀、非定常的污染物扩散问题，且模式的适用范围向中尺度、大尺度扩展。目前，数值计算已经成为研究的主流方法，研究范围也在逐步扩大。大气污染物扩散模式的应用受到地形、气象、大气污染物的物理化学特征、污染源特征等多种因素的制约，不同的扩散模式都有各自不同的考虑因素和适用范围，选择恰当的扩散模式能够较为准确地模拟污染物的扩散及分布，用于城市环境空气质量预报。目前，己经有许多发展成熟的污染物扩散模式应用于不同尺度的污染物扩散研究，其中不少模式经过验证都获得了较好的结果，国内的相关研究中尤其以高斯类模式的应用最为广泛。 2 大气污染物扩散模式的类型近年来针对工业点源的大气污染物扩散模式的应用，结合模式的理论核心以及应用特征，目前常用的扩散模式以高斯公式、拉格朗日方法、欧拉方法为基础。高斯模式是半经验型扩散模式，假定下风向的污染物浓度符合正态分布，是很多实用模式发展的基础。基于高斯理论的大气污染扩散模式被广泛应用于各种尺度的研究区域，其中适用于中小尺度的有ISC（Industrial Source Complex Model）、AERMOD （AMS/EPA Regulatory Model）、ISCST （Industrial Source Complex Short Term Model）、ADMS （Advanced Dispersion Modeling System）、ADMS-Urban等，应用于大尺度的有CALPUFF等。拉格朗日方法是用来跟随流体移动的粒子来描述污染物浓度及其变化。它是一种描述污染物分布的自然方式。在基于拉格朗日方法的大气污染物扩散模式中，适用于中小尺度的有LS （Lagrangian Stochastic Model）、TAPM （the Air Pollution Model）、LPM

水污染模型

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

大气污染扩散模型

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

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

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

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

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

第五章大气污染扩散

第五章大气污染扩散第一节大气结构与气象有效地防止大气污染的途径，除了采用除尘及废气净化装置等各种工程技术手段外，还需充分利用大气的湍流混合作用对污染物的扩散稀释能力，即大气的自净能力。污染物从污染源排放到大气中的扩散过程及其危害程度，主要决定于气象因素，此外还与污染物的特征和排放特性，以及排放区的地形地貌状况有关。下面简要介绍大气结构以及气象条件的一些基本概念。一、大气的结构气象学中的大气是指地球引力作用下包围地球的空气层，其最外层的界限难以确定。通常把自地面至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高度范围内的对流层称为大气边界层，该层空气流动受地表影响

城市空气污染程度的分析和预测模型

城市空气污染程度的分析和预测摘要本文讨论了有关城市污染程度、污染因素及污染扩散的问题。对于问题一，本文主要从大气污染、噪声污染和水体污染这三个面选取主要污染物，查阅北京、天津、上海、重庆和西安五座城市2007-2012年的年度平均污染数据，采用降维的思想，运用主成分分析法减少变量个数，再借助Matlab 软件计算各主成分的贡献率，分析知可选取前三个主成分作为衡量污染程度的标准，最后根据综合指标得到这五个城市的污染程度从高到低依次为：重庆、上海、北京、天津、西安。通过判断相关系数的大小，确定五个城市影响人们生活的主要污染因素是水污染，其四项指标依次为化学需氧量、总氮、总磷和氨氮。对于问题二，以北京市大气污染为例。首先，利用GPS 记录北京市14个城区观测点的位置，并查阅2013年污染指标2SO 、2CO 、5.2PM 与10PM 的污染数据，绘制出相应的空间浓度分布图，估计这四种污染物的大致污染源位置依次为：)100,110(附近、)83,130(附近，)85,125(附近和)80,132(附近；其次，根据污染扩散原理和方式，建立Cauchy 污染传播模型，根据各地区空气污染物的浓度分布，运用Matlab 软件对数据非线性拟合，得出扩散模型各参数的值，计算得出各项污染指标的污染源位置依次为：)3.97,5.115(，)3.85,2.128(，)8.80,1.129(和)6.87,5.125(；最后，比较污染物位置的计算值与实际值，发现误差相差较小，故模型建立较为合理。对于问题三，分析西安市的主要污染——大气污染。收集西安市2014年4.1-7.31日的空气污染数据，根据时间序列的平稳性特点及AIC 定阶准则选取合适的时间序列模型)11(ARMA ，，利用Matlab 软件对序列模型的各项参数进行估计并检验模型的合理性，并将模型用于数据预报。利用时间序列模型预测西安市未来10天的空气污染状况总体等级为良。对于问题四，基于问题一、二、三对污染因素的分析和污染扩散的特点，主要从减少污染物的产生和治理净化已产生的污染物两方面，针对大气污染、水体污染和噪声污染为相关部门提供合理化防治建议。关键词主成分分析；Cauchy 污染传播模型；时间序列模型；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公里处，放射性物质浓度的预测模型。

2011A数学建模城市污染物

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

编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：

城市表层土壤重金属污染分析摘要随着工农业的不断发展，我国土壤污染问题日益突出，进入土壤的持久性有机污染物、重金属污染等与日俱增，土壤环境面临着前所未有的压力。因此，了解土壤污染物来源是切实有效地控制土壤污染，保障环境安全和农业可持续发展的重要前提。本文针对城市表层土壤重金属污染问题，在合理的假设下，先绘制出重金属元素空间分布图；再根据单因子指数与内梅罗综合污染指数分析方法对各功能区污染程度分析，说明重金属污染的主要原因；接着，在此基础上建立二维污染物空间扩散模型，采用搜索算法在污染较为严重的区域搜索出适量的样本点，再采用多元非线性回归方法对二维扩散模型进行拟合，求解相应的参数；最后通过对第三问模型优缺点的分析，建立时空扩散模型，更加精确反映城市地质环境的演变情况。对于问题1，根据附件数据，调用Matlab 软件中的surf 函数绘制出重金属元素空间分布图，采用三维图像表示城区地貌特征，用灰度描绘该地受污染程度，清晰反映出在不同地貌，不同功能区的重金属污染物空间分布特点。选用单因子指数和内梅罗指数模型评价污染状况，运用Visual C++ 6.0编程计算出不同功能区单因子指数与内梅罗综合指数，由数据反映出工业区和交通区Hg 和Cu 严重超标，各功能区染程度由重到轻依次为：工业区、交通区、生活区、公园绿地区、山区。对于问题2，对各功能区的污染物含量进行分析，分析推测该城市存在规模较大或者数量较多的Hg,Cu,Zn 严重超标的工厂。为进一步揭示表层土壤重金属污染来源，对各个功能区的不同重金属污染元素进行相关性分析，运用spss 软件计算出各相关系数，揭示了不同功能区内土壤表层重金属污染的不同来源，分析可得，交通区污染主要来源于汽车尾气排排放，工业区污染主要来源于工厂的超标排放。对于问题3，根据污染物扩散的特征，建立一维污染物扩散模型：2 2 () 0()x u Q x Q e σ -- =，以此推广建立污染物的二维空间扩散模型：2 22 000()() (,)exp[]x x y y Q x y Q σ -+-=- ；再运用搜索算法在局部极值点附近搜索适量样本点；采用多元非线性拟合方法对二维空间扩散模型进行拟合，得到区域的最大值，即为污染源所在地,以Hg 的分布为例说明。对于问题4，通过对上述模型优缺点的分析，进一步搜集相关数据，依据重金属传播所具有的富集性和不易分解性等特征，在空间分布模型的基础上引进时间变量，突出时间对于污染物扩散的影响，建立时空扩散模型，以更好地研究地质环境的演化模式。本文描述了城市土壤重金属在城市不同功能区的含量分布、不同功能区各元素间分布格局的差异性、给出定量描述城区不同区域重金属的污染程度评价模型，通过对数据的分析以及数据间相关性的分析推测不同污染物的来源，并根据重金属的污染物的传播特征，建立确定污染源位置的模型，为合理规划和利用城市土壤、改造和提高城市环境质量、保障人类健康等提供重要依据。关键词：城市土壤重金属污染单因子指数内梅罗综合指数扩散模型非线性拟

高斯扩散模式在瞬间排放空气污染物模拟中的应用

高斯扩散模式在瞬间排放空气污染物模拟中的应用摘要：在文章中提出应用高斯模式模拟和预测在瞬间排放状况下空气污染等级，用FORTRAN 语言编写的高斯模式程序还可应用于区域污染影响评价中，模式不仅可以从GIS 中输入数据而且还可以应用GIS 格式输出结果。关键词：高斯模式空气污染地理信息系统瞬时污染源浓度场瞬间排放是指工业企业或电厂的事故性污染物排放，如贮油罐或输油管道发生事故等。排放的污染物污染了空气、土壤、地面及地下水，影响植被和影响环境。模拟瞬间空气污染要求得到污染区域面积、污染浓度和等级、污染预测等。本文提出用高斯模式的解析解来模拟和预测瞬间排放空气污染状况。基于烟羽扩散上的解析公式求解的高斯模式非常广泛的应用于评价区域污染状况。高斯数学模式作为一个污染物扩散的基础模式被国际原子能机构广泛推广。从瞬间点源污染源排放的污染物，其转换和扩散可以用以下的扩散方程来表示： t C ??+div(CV )=?(K ?C )+Ri +Q δ(t ?t 0)δ(x ?x 0)δ(y ?y 0)δ(z ?z 0) (1) 式中：C(x, y, z, t)为污染物浓度 V 为风速 K 为扩散系数 R 为污染物光化学转化率 Q 为污染物排放量

x 0, y 0, z 0为污染源相对坐标在一定的风速和扩散系数条件下，公式（1）有其高斯扩散模式的解析解。因此，污染物浓度值C 由点源污染源的高度H 决定。H 在高斯扩散模式中由下述公式计算： C （x,y,z,t ）= )() 2(22 22 2 22 2 2/) 2(2/) (2/)(2/) (2 /3z z y x wt h H z wt h z vt y ut x z y x e e e e Q σσσσσσσπ-++--------+ (2) 式中：t 为时间 Q 为排放量 u ，v ，w 为风速分别在x ，y ，z 方向的分量 σx , σy , σz 分别在x ，y ，z 方向的扩散系数 h 为点源高度 H 为混合层高度高斯模式中，假设X 轴与风向方向一致，Z 轴铅直向上，V=W=0。公式（2）可以转化为以下形式。 C （x,y,z,t ）= )() 2(22 22 2 22 2 2/) (2/) (2/2/) (2 /3z z y x H z H z y ut x z y x e e e e Q σσσσσσσπ+------+ (3) 从公式(3)我们可以看出，每一个烟团需要用不同的坐标系进行计算，当我们计算多源的污染浓度时，我们需要用到几个坐标系，这样计算起来很复杂。因此，公式（3）必须做相应的转化到同一个坐标系中。我们建立一个相对的坐标系，I 表示原点，坐标轴为ξ和η（见1）。并以I 为原点建立第二个坐标系，LX 表示X 轴，其方向与风向

污染物扩散模型

错误！未指定书签。该模块采用突发性水污染扩散模型，利用一维水质模型，通过对河段长度与扩散时间进行微分，后利用四点隐式差分格式进行模型的数值求解。详解如下： 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 时刻的浓度。通过距离量算来确定排污口与测量点的距离，计算测量点的浓度，并得到污染物在河道断面上的扩散规律。

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

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;

大气模型

大气模型发展简史与简介 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（不稳态高斯）对于比较大的区域

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

大气污染物扩散模式

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

城市空气质量模式研究进展

水污染模型

大气污染扩散模型

推荐-基于修正高斯扩散模型的城市表层土壤重金属污染探究 精品

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

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

第五章 大气污染扩散

城市空气污染程度的分析和预测模型

污染空气的扩散模型

2011A数学建模城市污染物

高斯扩散模式在瞬间排放空气污染物模拟中的应用

污染物扩散模型

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

大气模型

推荐-基于修正高斯扩散模型的城市表层土壤重金属污染探究精品

第五章大气污染扩散