725 数模 空气污染问题研究 (15 五一 B)

空气污染问题研究

问题

一般认为影响空气质量的主要因素有PM2.5、PM10、二氧化氮、二氧化硫、一氧化碳、臭氧、硫化氢、碳氢化合物和烟尘等，以京津冀地区为研究对象解决以下问题：
（1）参考现有国标和美标，建立衡量空气质量优劣程度等级的数学模型．
（2）查找数据并列出京津冀地区主要污染源及其污染参数，分析影响空气质量的主要污染源的性质和种类．
（3）建立单污染源空气污染扩散模型，描述其对周围空气污染的动态影响规律．现有河北境内某一工厂废气排放烟囱高50m，主要排放物为氮氧化物．早上9点至下午3点期间的排放浓度为406.92mg/m3，排放速度为1200m3/h；晚上10点至凌晨4点期间的排放浓度为1160mg/m3，排放速度为5700m3/h；通过你的扩散模型求解该工厂方圆51公里分别在早上8点、中午12点、晚上9点空气污染浓度分布和空气质量等级．
（4）建立多污染源空气污染扩散模型，并以汽车尾气污染源为例求解分析以下问题：北京在2020年1月15日已经连续三天发生重污染，假设从16日开始北京启动汽车单双号限行交通管制措施，求解北京市二环、四环、六环路在16日早上8点、中午12点、晚上9点时空气污染浓度梯度变化及空气质量等级．
（5）根据你们的模型和求解结果，分析总结影响空气质量的关键参数，为京津冀地区环保部门撰写一份建议报告，给出实现蓝天计划的可行性措施和建议．

分析

1. 参考现有国标和美标，建立衡量空气质量优劣程度等级的数学模型

考虑到国标和美标的合理性，我们只对轻度污染及以下进行了改进。

采用模糊综合评定法

2. 查找数据并列出京津冀地区主要污染源及其污染参数，分析影响空气质量的主要污染源的性质和种类．

收集数据
我们收集了北京、天津、河北承德、河北石家庄四地今年7月以来的PM2.5、PM10、CO、NO2、SO2、O3含量，日均浓度对比如下：

根据国家标准污染物基本项目浓度限值表，四地臭氧和PM2.5均轻微超标。根据这两种污染物的主要污染源。我们收集到了京津冀地区主要污染源及其污染参数。

采用系统合成的聚类算法，通过计算两类数据点间的距离，对最为接近的两类数据点进行组合，并反复迭代这一过程，直到将所有数据点合成一类，并生成聚类谱系图。

选用 欧几里得距离 类平均法 进行计算。

我们定义各个类畸变程度之和为该类重心与其内部成员位置距离的平方和。
根据聚合系数折线图可知，当类别数为2时，折线的下降趋势趋缓，故可将类别数设定为2。
生成的聚类谱系图如下：

同样方法分析石家庄与承德

根据聚类谱系图可知，
北京和天津的空气污染来源主要分为季节性污染和综合性污染。
季节性污染的主要成分是臭氧。由于2021年7月我国东部地区气温偏高，可能出现阶段性高温过程，北京天津以及周边地区气温较常年同期略偏高，降水接近常年到偏少，气象条件有利于臭氧污染的生成。这是一种具有显著显著的区域分布和季节变化特征的污染。综合性污染主要包括PM2.5、PM10、CO、NO2、SO2。这些污染物主要来源于二次硫酸盐、机动车尾气、化石燃料燃烧、生物质燃烧等人类生产活动。

承德和石家庄的空气污染来源主要分为重金属污染和非重金属污染。
重金属污染主要是由PM2.5构成。河北省共有涉重金属企业305家，其中涉铬企业198家、涉铅企业83家，这些企业大都没有污染防治设施、污染防治设施运行不正常或超标排放，导致PM2.5过度，而PM2.5大多含有重金属等有毒物质。非重金属污染主要是由03、PM10、CO、NO2、SO2构成，这些污染物主要来源于机动车尾气、化石燃料燃烧、生物质燃烧、道路扬尘等复杂的综合性人类生产活动。

3. 建立单污染源空气污染扩散模型，描述其对周围空气污染的动态影响规律．现有河北境内某一工厂废气排放烟囱高50m，主要排放物为氮氧化物．早上9点至下午3点期间的排放浓度为406.92mg/m3，排放速度为1200m3/h；晚上10点至凌晨4点期间的排放浓度为1160mg/m3，排放速度为5700m3/h；通过你的扩散模型求解该工厂方圆51公里分别在早上8点、中午12点、晚上9点空气污染浓度分布和空气质量等级．

本题专业性很强，我们大量阅读借鉴地理文献，其中以武汉大学戴万瑞的论文 大气点源污染扩散模拟研究 为主。

采用高斯扩散模型

根据大气污染物的大小形态可将其分为气态污染物和颗粒污染物。由于氮氧化物的粒径为纳米级，远小于15um，故采用气体扩散模型。

（1）稳定状态下模型解析

当大气运动稳定且污染源是源源不断排放时，污染物在空间上的分布状况是稳定的，也就是说污染物在某一空间位置的浓度是不随时间变化的。在大气条件呈无边界情况下，一个稳定点源连续排放的污染物扩散模型的表示方法如下：

上式成立是基于以下四个假定：1.污染物在无限空间内均匀稳定连续排放；2.连续点源排放的污染物在水平和垂直方向上的浓度分布均呈正态分布；3.污染物输送扩散过程中遵守质量守恒定律；4.区域风向，风速稳定

注意到，横向扩散系数和垂直扩散系数都是x的函数假定，这意味着预测点一带的烟气在x，y两个方向上尺度变化不能太大，也就是要控制烟气扩张角度。而考虑到河北省地处温带大陆性季风气候，大型工厂多建于沿海平原处，大气常处在不稳定状态且风速较大的大气条件下，可以将烟气扩张的角度控制在一个比较小的范围内，符合模型假定。

然而得出上式依然需要很强的假定，另外它也并没有把地形条件，其他气象条件，尤其是污染源的空间位置等因素考虑在内，所以该公式在实际中需要考虑上述因素并加以修正。

（2）地面反射作用下的修正

由于我国改革开放以来大力发展工业，京津唐工业基地的不断建设，基本满足其在可计算范围内扩散是无边界的。然而需要考虑地面作用所增加的浓度。当污染物到达地面后，由于污染物的密度很小，地面对污染物的反射作用非常大。为了简单起见，假设地面对污染物反射系数为1，即认为污染物到达地面后全部反射到大气中，采用镜像计算空间任一点污染物的浓度。任一点的污染物浓度值可以看作是两部分浓度的叠加值，一部分为不考虑地面影响时该点所具有的浓度，另一部分为地面反射作用增加的浓度。

（3）参数的确定

扩散参数指的是大气污染物在横向和垂直方向上的扩散参数,表征了湍流扩散运动的剧烈程度,这对污染物的浓度分布影响很大。通过上述对大气扩散模型的了解,可发现任意一点污染物地面浓度是由模型中的各个参数所决定,因此在模拟大气污染物扩散的过程中,要根据特定的地形及气象条件对参数做必要的修正并确定具体值。扩散参数的大小与污染点源到下风向的计算点距离及大气稳定度有关,目前对其的计算仍采取经验公式。

我们考察了河北省大型工业区主要分步，鉴于河北省常年风速达3-4.5m/s，大型工厂大多建于沿海，大气长期处于不稳定的状态，我们假定了稳定剂度数（P.S）处于B等级，因此只考虑有风条件下的参数：

依此我们拟合了题目要求的两个时间段排放在0.5小时以内的风向轴地面浓度图
MATLAB

clc; clear;
%以9:00-15:00为例，其中大气稳定度为B，排放时间在0.5h以内
x=(1:100:50000);
oy=0.281846*((x.*(x>=1 & x<=1000)).^0.914371)+0.396353*((x.*(x>1000)).^0.865014);
oz=0.12719*((x.*(x>=1 & x<=500)).^0.964435)+0.057025*((x.*(x>500).^1.09356));
fm=12*pi*(oy.*oz);
zs=-1250./(oz.^2);
f=(607.5./fm).*exp(zs);

plot(x,f)

呈正态分布，检验是合理的。

然而当排放时间变长，参数还需修正
查阅相关资料可以看到在实际中由于污染物在大气中的扩散是在湍流作用下形成的，而湍流统计量与取样时间的长短有关，当取样时间大于0.5h时，由于湍流速度脉动标准差增大，因而扩散参数也会随之增大，但两者增大程度不同，一般认为只对有增大作用，因此污染物由向下风轴方向扩散为主转为向Y方向扩散为主。因此需对扩散参数修正：

模拟出 NOx 在9:00-15:00时段不同采样时间下的X轴上的浓度值分布比较如图10所示

Y轴：

查阅第一问空气质量优劣程度评测表进行评估。

4. 建立多污染源空气污染扩散模型，并以汽车尾气污染源为例求解分析以下问题：北京在2020年1月15日已经连续三天发生重污染，假设从16日开始北京启动汽车单双号限行交通管制措施，求解北京市二环、四环、六环路在16日早上8点、中午12点、晚上9点时空气污染浓度梯度变化及空气质量等级．

一开始看汽车以为是线源，然而北京二环、四环、六环并不是一般意义上的狭长道路，得按照面源来分析。

MATLAB

% 8:00 与 12:00
clc; clear;
L=[2000 2600 1400;1800 2200 1200];
l=[32.7 65.3 187.6;32.7 65.3 187.6];
u=[3 3 3;2 2 2];
Q=(L.*l.*85.2)./3600;
x=70;
zz=0.2*70;
zw=0.14*70*(1+0.0003*70)^(-1/2);

C=((2/pi)^0.5).*(Q./(u.*zz)).*(exp(-15^2/(2*(zz)^2)));
C

%21:00
clc; clear;
L=[1600 1400 600;8000 7000 3000];
l=[32.7 65.3 187.6];
u=[3.2 3.2 3.2];
Q=(L.*l.*85.2)./3600;
x=70;
zz=0.2*70;
zw=0.14*70*(1+0.0003*70)^(-1/2);

C=((2/pi)^0.5).*(Q./(u.*zw)).*(exp(-15^2/(2*(zw)^2)));
C

（懒着打字了，他妈的，这是数学建模还是地理报告哇）

5. 根据你们的模型和求解结果，分析总结影响空气质量的关键参数，为京津冀地区环保部门撰写一份建议报告，给出实现蓝天计划的可行性措施和建议．

针对京津冀地区的空气污染问题,我们特此向环保部门提出如下几条意见:
（1）重污染企业限产甚至停产:针对重污染企业,进行强制限产停产,减少污染物的排放。
（2）部分施工工地停工,减少扬尘污染:加强对工地施工现场的扬尘管理.加强城市道路清扫保洁和冲洗。
（3）私家车单双号出行:机动车采取单双号行驶,增加公交的数量和班次,鼓励市民乘坐公共交通,绿色出行,限行大排量私家车。
（4）禁止焚烧秸秆、垃圾:枯秆焚烧对加重空气污染的作用明显,环保部门要采取保障措施完全杜绝焚烧秸秆、垃圾的现象。
（5）控制高耗能行业,促进产业升级:相比较其他较发达地区,京津翼地区单位工业增加值燃煤消耗量最大,高耗能行业也达70%以上,加重了对大气环境的污染。一方面,严格限制高耗能行业发展;另一方面,优化工业内部结构,提升高技术制造业份额.
（6）推广清洁能源的使用:引导出租车、公交车转型为电动汽车或混合动力汽车,推广太阳能电池在社区、公司、电力公司的应用,加强研发绿色高效能源的力度。

期待交流斧正