数学建模空气质量.docx
- 文档编号:10607613
- 上传时间:2023-02-21
- 格式:DOCX
- 页数:64
- 大小:1.43MB
数学建模空气质量.docx
《数学建模空气质量.docx》由会员分享,可在线阅读,更多相关《数学建模空气质量.docx(64页珍藏版)》请在冰豆网上搜索。
数学建模空气质量
2011高教社杯全国大学生数学建模竞赛
承诺书
我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括、电子、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。
如有违反竞赛规则的行为,我们将受到严肃处理。
我们参赛选择的题号是(从A/B/C/D中选择一项填写):
A
我们的参赛报名号为(如果赛区设置报名号的话):
所属学校(请填写完整的全名):
中国农业大学
参赛队员(打印并签名):
1.王万能
2.董祥祥
3.靖翔
指导教师或指导教师组负责人(打印并签名):
日期:
2014年9 月8日
赛区评阅编号(由赛区组委会评阅前进行编号):
2011高教社杯全国大学生数学建模竞赛
编号专用页
赛区评阅编号(由赛区组委会评阅前进行编号):
赛区评阅记录(可供赛区评阅时使用):
评
阅
人
评
分
备
注
全国统一编号(由赛区组委会送交全国前编号):
全国评阅编号(由全国组委会评阅前进行编号):
市空气质量状况探究
摘要
空气与我们的生存是息息相关的,它直接参与人体的新代、物质代和体温调节等过程。
随着现代工业和交通的迅猛发展,烟尘和汽车尾气等的排放,超越了大气的自净界限,接踵而至的却是一个十分严峻的问题—大气污染。
作为我国的首都,的空气质量又是什么情况呢?
我们通过数学模型来分析这个问题。
对于第一问“查找相关数据(包括近期连续数据及零星数据)”,我们认为,由于空气质量指数(AQI)是2012年上半年才被出台规定取代原有空气污染指数(API)的,且参与空气质量评价的主要污染物为细颗粒物(PM2.5)、可吸入颗粒物(PM10)、二氧化硫、二氧化氮、臭氧、一氧化碳等六项。
为方便后续几问的研究,我们收集了2013年至2014年8月每日的AQI、空气质量等级、PM2.5、PM10、SO2、NO2、CO这7项指标的数据(即连续型数据)和2008年至2013年PM10、SO2、NO2、CO的年度平均指标数据(即离散数据)。
第二问中,本文从时间和空间两个维度分析污染情况。
从时间角度分析时,为保证充分利用数据,本文分别从污染物年浓度平均变化和持续污染日数和空气质量超标日数年变化两个层次对空气质量进行了分析,从而归纳出空气污染的年度特征。
从空间角度分析时,本文根据市各区县PM10、SO2、NO2、CO的年度平均数据,将各区县各项指标浓度反映在了的区域地图上,绘制得各项指标的浓度空间分布示意图,同时将两年数据求平均值并按照PM10降序排列,结合图、表定量分析并得出结论。
第三文中,本文以各个监测点为圆心,以各监测点所能测得的准确数据的最大距离为半径作圆,根据所有圆是否将市全部覆盖判断空气质量监测点的布局是否合理。
最终本文得出结论:
市空气质量监测点的布局是合理的。
对于第四问“分析空气污染的主要成因和解决办法”,本文以抽样的方式,抽取每月1号、5号、10号、15号、20号、25号、30号,二月取28号的数据为样本,统计2013年9月1日到2014年8月31日抽样日期每天的23个检测站点监测到的当日首要污染物为“颗粒物”、“臭氧”、“二氧化氮”、“二氧化硫”、“一氧化碳”、及“无”的监测点个数,求出这84天各项指标为首要污染物的频率,比较频率大小,得出空气污染的主要气体是颗粒物、臭氧、二氧化氮、一氧化碳,并据此提出解决办法。
第五问中,本文选取主成分神经网络模式识别的方法,对未来一周(9月6日至9月12日)的空气状况进行预测。
我们选择预报下周空气质量的等级,使用的是MATLAB提供的模式识别工具箱。
我们将空气质量看做六个模式:
优,良,轻度污染,中度污染,重度污染,严重污染;并将spss软件求解出来的主成分与对应的模式交给神经网络学习。
训练好后用于预测。
根据主成分分析的方法,对全部160个样本的8个因子进行分析。
将数据导入spss,设置抽取系数为0.5,得到相关矩阵;然后由解释总方差表格和成分矩阵找出spss提取的主成分,将八维数据降成五维数据;并将样本转化,用于神经网络的训练;并把不同的空气质量等级二进制化,得到训练结果后进行预测。
将要预测时间段的平均气温,最高气温,最低气温,大气压,相对湿度,相对风速,水平能见度,降水量数据按照公式变为五维数据并输入网络进行求解,最终得到未来一个星期的空气质量为:
重度污染,重度污染,优,良,良,轻度污染,轻度污染。
第六问,本文假设自己想从事相关商业活动,给出了计划及广告词。
关键词:
空气质量 AQI PM2.5
主成分神经网络 空气污染 PM10
1问题重述
空气与我们的生存是息息相关的,它直接参与人体的新代、物质代和体温调节等过程。
一个人每天呼吸的空气约为1万多升,折合质量约为12.9 kg,约为每天所需食物和饮水量的10倍。
随着现代工业和交通的迅猛发展,烟尘和汽车尾气等的排放,超越了大气的自净界限,接踵而至的却是一个十分严峻的问题—大气污染。
对大气污染影响较大的污染物有:
硫氧化物、氮氧化物、碳氢化合物、碳氧化合物、粉尘等。
大气污染物对我们身体的影响是巨大的,而2013年以来日益频发的雾霾天气已经严重影响到人们的日常生活。
“生态文明”是被列为当前的一大重要议题,“我们将加强生态环境保护,扎实推进资源节约,为人民创造良好生产生活环境,为应对全球富气候变化作出新的贡献。
”,的一系列讲话为坚持节约资源和保护环境基本国策,努力走向生态文明新时代赋予了新容,提出了新要求,彰显出他对于生态文明领域建设的决心与魄力。
本文中,我们
(1)查找相关数据(包括近期连续数据及零星数据);
(2)分析空气污染情况(时间、空间);
(3)分析空气质量监测点的布局是否合理;
(4)分析空气污染的主要成因及解决办法;
(5)对未来一周进行预测;
(6)假设自己想从事相关商业活动,给出计划及广告词。
(不超过一页)
2问题分析
对于第一问“查找相关数据(包括近期连续数据及零星数据)”,我们认为,由于空气质量指数(AQI)是2012年上半年才被出台规定取代原有空气污染指数(API)的,且参与空气质量评价的主要污染物为细颗粒物(PM2.5)、可吸入颗粒物(PM10)、二氧化硫、二氧化氮、臭氧、一氧化碳等六项。
为方便后续几问的研究,我们收集了2013年至2014年8月每日的AQI、空气质量等级、PM2.5、PM10、SO2、NO2、CO这7项指标的数据(即连续型数据)和2008年至2013年PM10、SO2、NO2、CO的年度平均指标数据(即离散数据)。
第二问中,本文从时间和空间两个维度分析污染情况。
从时间角度分析时,为保证充分利用数据,本文分别从污染物年浓度平均变化和持续污染日数和空气质量超标日数年变化两个层次对空气质量进行了分析,从而归纳出空气污染的年度特征。
从空间角度分析时,本文根据市各区县PM10、SO2、NO2、CO的年度平均数据,将各区县各项指标浓度反映在了的区域地图上,绘制得各项指标的浓度空间分布示意图,同时将两年数据求平均值并按照PM10降序排列,结合图、表定量分析并得出结论。
第三文中,本文以各个监测点为圆心,以各监测点所能测得的准确数据的最大距离为半径作圆,根据所有圆是否将市全部覆盖判断空气质量监测点的布局是否合理。
最终本文得出结论:
市空气质量监测点的布局是合理的。
对于第四问“分析空气污染的主要成因和解决办法”,本文以抽样的方式,抽取每月1号、5号、10号、15号、20号、25号、30号,二月取28号的数据为样本,统计2013年9月1日到2014年8月31日抽样日期每天的23个检测站点监测到的当日首要污染物为“颗粒物”、“臭氧”、“二氧化氮”、“二氧化硫”、“一氧化碳”、及“无”的监测点个数,求出这84天各项指标为首要污染物的频率,比较频率大小,得出空气污染的主要气体是颗粒物、臭氧、二氧化氮、一氧化碳,并据此提出解决办法。
第五问中,本文选取主成分神经网络模式识别的方法,对未来一周(9月6日至9月12日)的空气状况进行预测。
我们选择预报下周空气质量的等级,使用的是MATLAB提供的模式识别工具箱。
我们将空气质量看做六个模式:
优,良,轻度污染,中度污染,重度污染,严重污染;并将spss软件求解出来的主成分与对应的模式交给神经网络学习。
训练好后用于预测。
根据主成分分析的方法,对全部160个样本的8个因子进行分析。
将数据导入spss,设置抽取系数为0.5,得到相关矩阵;然后由解释总方差表格和成分矩阵找出spss提取的主成分,将八维数据降成五维数据;并将样本转化,用于神经网络的训练;并把不同的空气质量等级二进制化,得到训练结果后进行预测。
将要预测时间段的平均气温,最高气温,最低气温,大气压,相对湿度,相对风速,水平能见度,降水量数据按照公式变为五维数据并输入网络进行求解,最终得到未来一个星期的空气质量为:
重度污染,重度污染,优,良,良,轻度污染,轻度污染。
第六问,本文假设自己想从事相关商业活动,给出了计划及广告词。
3假设与约定
假设:
(1)假设预测时间段没有发生重大的天气异常;
(2)假设预测时间段没有污染源;
(3)假设以pm2.5为空气质量等级划分的标准为:
ps:
pm2.5标准
0-35 优
35-75良
75-115轻度污染
115-150中度污染
150-250重度污染
250- 严重污染
4符号说明
符号
PM2.5
PM10
AQI
含义
细颗粒物
可吸入颗粒物
空气质量指数
5模型建立与求解
5.1问题一
本题主要收集了市环境污染近期相关连续数据与离散数据。
由于空气质量指数(AQI)是2012年上半年才被出台规定取代原有空气污染指数(API)的,且参与空气质量评价的主要污染物为细颗粒物(PM2.5)、可吸入颗粒物(PM10)、二氧化硫、二氧化氮、臭氧、一氧化碳等六项。
为方便后续几问的研究,我们收集了2013年至2014年8月每日的AQI、空气质量等级、PM2.5、PM10、SO2、NO2、CO这7项指标的数据(即连续型数据)和2008年至2013年PM10、SO2、NO2、CO的年度平均指标数据(即离散数据)。
见附录。
5.2问题二
本问主要从时间和空间两个维度分析空气污染情况。
5.2.1模型建立
(1)本题主要从时间和空间两个维度分析空气污染情况。
(2)从时间角度分析时,为保证充分利用数据,本文分别从污染物年浓度平均变化和持续污染日数和空气质量超标日数年变化两个层次对空气质量进行了分析,从而归纳出空气污染的年度特征。
(3)从空间角度分析时,本文根据市各区县PM10、SO2、NO2、CO的年度平均数据,将各区县各项指标浓度反映在了的区域地图上,绘制得各项指标的浓度空间分布示意图,同时将两年数据求平均值并按照PM10降序排列,结合图、表定量分析并得出结论。
5.2.2模型求解
5.2.2.1 时间特征
(1)按污染物年浓度平均变化
本文分析了2008年至2013年PM10、CO、NO2、SO2的平均浓度变化,变化情况如下表:
时间
PM10
CO
NO2
SO2
2008
122
1.4
49
36
2009
121
1.6
53
34
2010
121
1.5
57
32
2011
114
1.4
55
28
2012
109
1.4
52
28
2013
108.1
56
26.5
变化趋势如下
由上表我们可以得出:
近年来,市NO2浓度处于缓慢上升趋势;CO浓度与往年相比基本保持持平;PM10浓度略有下降;SO2年平均浓度处于下降趋势。
(2)按持续污染日数和空气质量超标日数年变化
统计分析2001-2010年空气质量超标日数和空气持续污染日数年变化图,我们可以看出:
2001年-2009年年空气持续污染日数在年超标日数中所占比例都在50%以上。
可以说一半以上的空气质量超标日都是持续污染日。
近年来的年空气质量超标日数和持续污染日数总体呈下降趋势。
但年持续污染日数在年超标日数中所占的比例一直较大。
综上,空气污染年际变化特征如下:
污染物排放量的大力削减使市空气质量得到明显改善,年空气质量超标日数和持续污染日数总体也呈下降趋势,但是持续污染对全年空气质量超标日数贡献依然较大,空气质量仍有待改善。
5.2.2.2 空间特征
为了分析污染的空间分布情况,我们将2012与2013年的SO2,NO2,以及可吸入颗粒物的各区县浓度反映在了的区域地图上。
(1)2013年:
各区县各项指标年平均数据
区县
PM2.5
SO2
NO2
PM10
东城区
93.6
26.8
58.0
109.6
西城区
91.7
28.8
59.6
112.5
区
91.3
29.7
64.0
112.4
海淀区
98.1
26.9
63.6
115.0
丰台区
96.9
28.1
57.5
118.5
石景山区
92.8
24.9
63.3
116.4
门头沟区
91.1
24.6
51.8
114.8
房山区
106.8
31.2
61.9
131.7
通州区
105.7
38.6
55.8
123.5
顺义区
84.8
20.8
44.8
98.5
大兴区
107.8
33.7
65.7
130.3
昌平区
79.2
25.9
43.5
94.7
平谷区
84.8
20.6
35.0
98.7
怀柔区
76.1
22.3
37.9
95.3
密云县
71.6
21.3
43.6
85.9
延庆县
68.0
19.2
34.4
78.3
亦庄
104.9
33.6
57.5
123.2
空气中可吸入颗粒物浓度空间分布示意图
空气中细颗粒物颗粒物浓度空间分布示意图
空气中二氧化氮浓度空间分布示意图
空气中二氧化硫浓度空间分布示意图
(2)2012年:
各区县各项指标年平均数据
区县
SO2
NO2
PM10
东城区
29
56
113
西城区
32
59
111
区
29
60
114
海淀区
31
54
114
丰台区
28
57
113
石景山区
24
55
124
门头沟区
29
51
109
房山区
36
59
122
通州区
42
56
119
顺义区
23
45
98
大兴区
35
64
124
昌平区
29
40
97
平谷区
24
35
98
怀柔区
22
30
87
密云县
29
40
85
延庆县
26
39
82
亦庄
36
54
126
(3)两年数据平均值(按关键字PM10降序排列)
区县
PM10
SO2
NO2
大兴区
127.15
27.9
57
房山区
126.85
30.4
59.3
亦庄
124.6
29.35
62
通州区
121.25
28.95
58.8
石景山区
120.2
28.05
57.25
丰台区
115.75
24.45
59.15
海淀区
114.5
26.8
51.4
区
113.2
33.6
60.45
门头沟区
111.9
40.3
55.9
西城区
111.75
21.9
44.9
东城区
111.3
34.35
64.85
平谷区
98.35
27.45
41.75
顺义区
98.25
22.3
35
昌平区
95.85
22.15
33.95
怀柔区
91.15
25.15
41.8
密云县
84.45
22.6
36.7
延庆县
80.15
34.8
55.75
(4)结合图表定量分析
由上表易知2013年:
(4.1)
按PM2.5值由大到小,各区县排列为大兴区>房山区>通州区>亦庄>海淀区>丰台区>东城区>石景山区>西城区>区>门头沟区>平谷区>顺义区>昌平区>怀柔区>密云县>延庆县,按此顺序,各区县空气质量愈好;
(4.2)
各区县的SO2、NO2、PM10取值及大小排序;
易知2012年各区县的SO2、NO2、PM10取值及大小排序;
易知2013年和2012年各区县的SO2、NO2、PM10平均取值及大小排序;
结合图,根据各区县的地理位置分布,我们得出结论:
市空气质量南北差异显著。
位于北部的生态涵养发展区空气质量优于其他区域。
5.3问题三
本题主要探讨市环境质量监测点的布局是否合理。
模型建立及求解
5.3.1环境空气质量评价城市点的定义:
以监测城市建成区的空气质量整体状况和变化趋势为目的而设置的监测点,参与城市环境空气质量评价。
其设置的最少数量根据本标准由城市建成区面积和人口数量确定。
每个环境空气质量评价城市点代表围一般为半径500米至4千米,有时也可扩大到半径4千米至几十千米(如对于空气污染物浓度较低,其空间变化较小的地区)的围。
可简称城市点。
5.3.2
这是市环境监测点的分布图:
按照比例尺计算,距离最近的两个监测点的距离约为10公里。
,这样说来就按照最保守的估计每一个监测点的检测半径为5公里,则每一个监测点的检测的面积为78.5平方公里,市区10个监测点总共为785平方公里,正好和市区的面积相仿。
所以,从面积上来说,监测点的分布是合理的。
而对于其他的区来说,位于郊区人口密度明显比市区低许多,因而空气质量的变化不大,一个监测点所能代表的围的半径也有数十公里。
5.4问题四
本题主要探究空气污染的主要成因,并提出解决办法。
5.4.1模型建立
(1)以抽样的方式,抽取每月1号、5号、10号、15号、20号、25号、30号,二月取28号的数据为样本,统计2013年9月1日到2014年8月31日抽样日期每天的23个检测站点监测到的当日首要污染物为“颗粒物”、“臭氧”、“二氧化氮”、“二氧化硫”、“一氧化碳”、及“无”的监测点个数。
(2)求出这84天“颗粒物”、“臭氧”、“二氧化氮”、“二氧化硫”、“一氧化碳”、及“无”为首要污染物的频率。
(3)比较污染成分为首要污染物的出现频率,随即便可得出空气污染的主要成因,并据此提出解决办法。
5.4.2模型求解
(1)统计2013年9月1日到2014年8月31日抽样日期每天的23个检测站点监测到的当日首要污染物为“颗粒物”、“臭氧”、“二氧化氮”、“二氧化硫”、“一氧化碳”、及“无”的监测点个数。
统计结果见下表:
序号
颗粒物
臭氧
二氧化氮
二氧化硫
一氧化碳
无
合计
1
23
0
0
0
0
0
23
2
0
21
1
0
0
1
23
3
0
23
0
0
0
0
23
4
0
23
0
0
0
0
23
5
0
21
0
0
0
0
21
6
1
14
3
0
0
5
23
7
23
0
0
0
0
0
23
8
7
16
0
0
0
0
23
9
22
1
0
0
0
0
23
10
0
23
0
0
0
0
23
11
0
23
0
0
0
0
23
12
1
22
0
0
0
0
23
13
1
22
0
0
0
0
23
14
22
1
0
0
0
0
23
15
0
23
0
0
0
0
23
16
0
23
0
0
0
0
23
17
20
2
0
0
0
1
23
18
4
19
0
0
0
0
23
19
0
23
0
0
0
0
23
20
23
0
0
0
0
0
23
21
0
23
0
0
0
0
23
22
22
1
0
0
0
0
23
23
4
12
8
0
0
1
25
24
16
2
5
0
0
0
23
25
23
0
0
0
0
0
23
26
14
9
0
0
0
0
23
27
23
0
0
0
0
0
23
28
0
22
0
0
0
1
23
29
23
0
0
0
0
0
23
30
19
3
0
0
0
1
23
31
23
0
0
0
0
0
23
32
22
0
1
0
0
0
23
33
20
1
2
0
0
0
23
34
23
0
0
0
0
0
23
35
13
10
0
0
0
0
23
36
23
0
0
0
0
1
24
37
4
0
9
0
0
10
23
38
22
0
0
0
0
1
23
39
22
0
1
0
0
0
23
40
2
5
1
0
0
15
23
41
23
0
0
0
0
0
23
42
9
3
11
0
0
0
23
43
23
0
0
0
0
0
23
44
23
0
0
0
0
0
23
45
3
0
6
1
0
13
23
46
23
0
0
0
0
0
23
47
23
0
0
0
0
0
23
48
23
0
0
0
0
0
23
49
23
0
0
0
0
0
23
50
19
0
4
0
0
0
23
51
23
0
0
0
0
0
23
52
18
0
5
0
0
0
23
53
22
0
1
0
0
0
23
54
2
0
0
0
0
21
23
55
8
0
7
0
1
8
24
56
23
0
0
0
0
0
23
57
21
0
2
0
0
0
23
58
13
0
5
0
0
5
23
59
14
0
0
0
1
8
23
60
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 建模 空气质量