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

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;

sigz=0.0757182*x.^1.00770;

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

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)));

%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.^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

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

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

sigy=0.3914238*x.^0.865014;

sigz=0.0757182*x.^1.00770;

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

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表明空气质量极差，已严重污染。

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