环境管理_城市空气质量评估与预测

2010年西南交通大学“新秀杯”数学建模竞赛参赛论文论文题目：城市空气质量评估及预测 参加2011年建模创新团队选拔： 是 （“是” 或 “否”）2010年11月城市空气质量评估及预测 摘要：本文对我国的成都，杭州，北京，上海，广州，拉萨，乌鲁木齐，郑州，武汉，西安等10座城市的每日空气质量详细列表进行科学分析，利用层次分析法和指数平滑法等数学建模方法对其空气质量进行研究，综合考虑各种因素建立如下数学模型。1、十个城市空气污染严重程度的科学排名。本文采用层次分析法来对10个城市的环境污染情况进行比较。用excel统计出各个城市不同级别的污染天数。在采用层次分析法，首先建立层次结构模型然后构造成对比较阵，用matlab软件计算权向量并做一致性检验，计算组合权向量并做组合一致性检验。最后根据权重的大小就可以比较出10个城市空气污染严重程度的科学排名。2、成都市11月的空气质量状况预测。本文采用指数平滑的方法来建立数学模型。根据前面的统计数据，考虑到污染级别不同和首要污染物的种类两个因素来预测成都11月的空气质量状况。3、分析影响城市空气污染程度的主要因素。先通过excel统计出每个城市的各种首要污染物所占总天数的比例，然后再综合考虑各方面因素找出造成该污染物超标的以原因，以找出影响城市空气污染程度的主要因素。最后本文就10个城市的空气污染严重程度的排名给出了相关的分析以及应对策略。就成都11月份的空气质量状况预测给出出行和生活方面的建议。并结合当下倡导建设环境友好型和资源节约型社会出发，就如何兼顾经济发展与环境保护给出指导性建议。关键词：空气质量污染等级 层次分析法 指数平滑法 首要空气污染物一 问题的提出1.1 背景介绍随着科技的发展，工业的进步和全球人口急剧增多的因素的影响，人们赖以生存的环境遭到了很大的破坏，很多地区相继出现了酸雨、物种灭绝、土地沙化等环境问题，环境问题已经成为当今世界各国普遍关注的问题之一，也是21 世纪人类面临的重大挑战。我国是一个人口大国，城市众多，人口密集。但由于工业的发展，我们的很多城市都受到了不同程度的污染，尤其是空气的污染，直接对我们造成伤害。空气中的污染物主要是可吸入颗粒、二氧化硫、二氧化氮等物质。1.2 问题的解决我们通过统计不同城市的每天空气质量详细列表，运用层次分析法可以对不同城市的空气质量进行统计和比较。在预测成都11月的空气质量状况时采用指数平滑的方法来建立数学模型。二 基本假设1.根据空气污染指数（）来划分为个等级：API值小于等于50，空气质量为优，相当于国家空气质量一级标准；API值大于50且小于等于100，表明空气质量良好，相当于达到国家质量二级标准；API值大于100且小于等于200，表明空气质量为轻度污染，相当于国家空气质量三级标准；API值大于200表明空气质量差，称之为中度污染，为国家空气质量四级标准；API大于300表明空气质量极差，已严重污染。2数据取自国家统计局，具有使用价值。三问题的分析1. 不同城市的空气质量统计和比较(1) 最大特征值(max) 的Matlab 计算方法：V,D=eig(A),其中A为待计算特征值的矩阵，D 为对角矩阵，其对角元素为A的特征值，最大的即为(max) 。(2)一致性指标CI 计算方法： CI=(-n)/(n-1);（其中为矩阵A的最大特征值，n为矩阵的阶数）(3)随即一致性指标RI 的计算方法：RI 与n 有如下关系，如表n12345678910RI000.580.901.121.241.321.411.451.49(5)权重计算方法计算矩阵A 的特征根及特征向量，将所求的特征向量单位化后得到的就是权重值。2. 成都11月份空气质量进行预测 指数平滑法（Exponential Smoothing，ES）是布朗(Robert G.Brown)所提出，布朗、认为时间序列的态势具有稳定性或规则性，所以时间序列可被合理地顺势推延；他认为最近的过去态势，在某种程度上会持续的未来，所以将较大的权数放在最近的资料。指数平滑法是在移动平均法基础上发展起来的一种时间序列分析预测法，它是通过计算指数平滑值，配合一定的时间序列预测模型对现象的未来进行预测。其原理是任一期的指数平滑值都是本期实际观察值与前一期指数平滑值的加权平均。基本模型如下：Ft+1=aYt + (1-a)FtFt+1为第t+1时期的时间序列预测值Yt 为第t时间的时间序列的实际值Ft 为第t时间的时间序列的预测值a 为平滑系数, 其取值范围为0,1；四 模型的建立与求解4.1建立数学模型给出十个城市空气污染严重程度的科学排名 1. 将研究目标（Z）、因素（P）、对象（C）按相关关系分成目标层Z、准则层P、对象层C。层次结构图如图所示： 空气质量一级空气质量二级空气质量三级空气质量四级空气质量五级空气质量 成都 杭州 北京 上海 2.给出空气质量一级，二级，三级，四级和五级两两成对比较的判断矩阵P 污染级别一级 二级三级四级五级12345根据上图得出如下例两两成对比较的判断矩阵P权重五级四级三级二级一级五级1/31 5/45/35/25四级4/154/514/324三级1/53/53/413/23二级2/152/51/22/314一级1/151/51/41/31/21由表中数据， 计算可知： (max) = 5.00 ， CI = 0.00 ， RI = 1.12 ，CR = 0.00 0.1。因为CR = 0.00 0.1，所以此排序有满意的一致性。3.给出对象层对准则层的各个因素的判断矩阵并进行分析。由于各个城市只存在污染程度的不同，所以它们两者之间各因素之间的关系。2010空气污染指数统计优良轻度污染中度污染严重污染成都492243300杭州502203002北京352015921上海12715614701广州951931000拉萨21291301乌鲁木齐192166651郑州152493611武汉522054802西安142444301成都杭州 北京上海广州拉萨乌鲁木齐郑州武汉西安权重成都149/507/549/12749/9549/21249/1949/1549/527/20.073杭州50/49110/750/12725/9525/10650/1910/325/2625/70.074北京5/77/10135/12735/9535/21235/197/335/525/20.052上海127/49127/50127/351127/95127/212127/19127/15127/52127/140.190广州95/4919/1019/795/127195/212519/395/5295/140.142拉萨212/49106/25212/35212/127212/951212/19212/1553/13106/70.317乌鲁木齐19/4919/5019/3519/1271/519/212119/1519/5219/140.028郑州15/493/103/715/1273/1915/21215/19115/5215/140.022武汉52/4926/2552/3552/12752/9513/5352/1952/15126/70.077西安2/77/252/514/12714/957/10614/1914/157/2610.021根据表中数据，类比（2）中方法，计算出各种不同污染等级对不同城市的权重空气质量一级对10个城市的不同权重表由表中数据， 计算可知： (max) = 10.00 ， CI = 0.00 ， RI = 1.49 ，CR = 0.00 0.1 同理，可以计算出其余空气质量等级对10个城市的不同权重。计算方法类似，用matlab软件的计算过程详见本文附表。空气质量二级对10个城市的不同权重表城市成都杭州北京上海广州拉萨乌鲁木齐郑州武汉西安权重0.1120.1100.1000.0780.0960.0450.1080.1240.1030.122由表中数据， 计算可知： (max) = 10.00 ， CI = 0.00 ， RI = 1.49 ，CR = 0.00 0.1 空气质量三级对10个城市的不同权重表城市成都杭州北京上海广州拉萨乌鲁木齐郑州武汉西安权重0.0690.0630.1240.3090.0210.0060.1390.0760.1010.091由表中数据， 计算可知： (max) = 10.00 ， CI = 0.00 ， RI = 1.49 ，CR = 0.00 0.1 空气质量四级对10个城市的不同权重表城市成都杭州北京上海广州拉萨乌鲁木齐郑州武汉西安权重000.250000.6250.12500由表中数据， 计算可知： (max) = 10.00 ， CI = 0.00 ， RI = 1.49 ，CR = 0.00 0.1 空气质量五级对10个城市的不同权重表城市成都杭州北京上海广州拉萨乌鲁木齐郑州武汉西安权重00.20.10.100.10.10.10.20.1由表中数据， 计算可知： (max) = 10.00 ， CI = 0.00 ， RI = 1.49 ，CR = 0.00 0.1 4. 模型求解进行层次总排序，方法： 将上面5个空气质量等级对10个城市的不同权重表单位化后作为列向量构成105矩阵，和空气质量一级，二级，三级，四级和五级两两成对比较的判断矩阵P相乘，结果便得到10个城市的权重值。根据上述三问题的分析中的假设可知，权重值越大，表明空气污染情况越严重。因此，将10个城市的权重值按照从小到大依次排序，得出的结果便是10个城市的空气污染严重程度的排名。最终结果如下表所示：一级(0.067)二级(0.133)三级(0.200)四级(0.267)五级(0.333)总权重成都0.0730.1120.069000.杭州0.0740.1100.06300.20.北京0.0520.1000.1240.250.10.上海0.1900.0780.30900.10.广州0.1420.0960.021000.拉萨0.3170.0450.00600.10.乌鲁木齐0.0280.1080.1390.6250.10.郑州0.0220.1240.0760.1250.10.武汉0.0770.1030.10100.20.西安0.0210.1220.09100.10.总的一致性检验CR = 0.067 0 + 0.133 0 + 0.200 0 +0.2670+0.3330= 0 0.1此结果说明排序结有非常满意的一致性。5.结果分析 根据10个城市的总权重值进行从小到大依次排序，空气污染严重程度的排名如下：(1) 广州(2) 成都(3) 拉萨(4) 西安(5) 杭州(6) 郑州(7) 武汉(8) 上海(9) 北京(10) 乌鲁木齐结论显示广州的空气质量状况最好，而上海，北京，乌鲁木齐的空气质量状况较差。广州今年在筹备2010广州亚运会，因此在城市环境保护方面下大力气进行整治，空气质量状况得到明显改善。而北京上海作为我国的一线城市，城市人口密度大，机动车保有量大，交通拥堵，工业化程度高，因此城市空气污染情况严重。乌鲁木齐由于地理位置的因素，处在大西北，经常遭受沙尘暴等恶劣天气影响，因此空气污染情况严重，特别是可吸入颗粒明显超标。4.2 建立模型对成都市11月的空气质量状况进行预测预测成都市11月的空气质量状况主要采用指数平滑的方法来建立数学模型。1.用excel统计出成都市2005年到2010年的空气污染指数，如下表所示：成都20052010空气污染指数统计优良轻微污染轻度污染中度污染中度重污染重污染2005212666015004200635251631042020073127456220020082529143700020093427838111002010492242850002.根据上表，做出统计表，然后做出分布图。如下所示：年份占全年百分比/%空气污染指数0100100200大于200200578.356220.54791.0959200678.356220.00001.6438200783.561615.89040.2479200886.575313.69860200986.187813.535902762201089.215710.78430分布图如下所示：3.分析2005年到2010 年的空气污染指数由这个分布图我们可以看出每年成都空气污染指数的大体情况，结合2010年每月的空气污染指数，我们对成都11月份的空气污染指数进行预测 。运用指数平滑法进行计算。分析: 为了预测成都11月份的空气污染指数，不仅要知道前面几月的实际空气污染指数，还要知道前面几个月的预测值。例如：我们已经知道了成都11月份前4天的空气污染指数。(这只是一个举例)F1=Y1，由于空气污染指数有波动，但长期趋势变化不大时，可选稍大的值，常在0.10.4之间取值；我们取a=0.3进行计算，则F2=0.3Y1+0.7F1,F3=0.3Y2+0.7F2日期空气污染指数移动平均法预测值预测偏差偏方平方值2010/11/1762010/11/284768642010/11/36078.4-18.4338.62010/11/46872.88-4.8823.8合计426.4解得2010/11/5的空气污染指数为F5=0.3Y4+0.7F4 =0.368+0.772.88 71第5日的预测偏差的平方值为=426.4/3=142.1不同的取值对第5日的预测值和预测偏差的平方值的影响，公式写为Ft+1=a .Yt+(1-a)Ft=a .Yt+Ft-aFt=Ft+a(Yt-Ft)选取不同的值获得结果如下表所示:a第5日预测值第5日预测偏差平方估计值0761280.1741210.2731360.371142.10.470600.569560.6681680.7671700.8671880.9662101.068235当取0.5时，用指数平滑法求得的第5日的预测值为69，预测偏差的平方估计值最小为56。因此，0.5为最适合的取值，69为第5日最精确的预测值。其它日子的空气污染指数计算方法如上。4.最终对成都市11月的空气质量状况的预测情况如下： 成都市11月的空气质量状况预测详细列表日期污染指数首要污染物空气质量级别空气质量状况2010-11-3074可吸入颗粒物良2010-11-2972可吸入颗粒物良2010-11-2851可吸入颗粒物良2010-11-2768可吸入颗粒物良2010-11-2644优2010-11-2581可吸入颗粒物良2010-11-2477可吸入颗粒物良2010-11-2347优2010-11-2249优2010-11-2187可吸入颗粒物良2010-11-2077可吸入颗粒物良2010-11-1969可吸入颗粒物良2010-11-1855可吸入颗粒物良2010-11-1767可吸入颗粒物良2010-11-1650优2010-11-1547优2010-11-1453可吸入颗粒物良2010-11-1391可吸入颗粒物良2010-11-1298可吸入颗粒物良2010-11-11117可吸入颗粒物1轻微污染2010-11-10106可吸入颗粒物1轻微污染2010-11-9115可吸入颗粒物1轻微污染2010-11-8105可吸入颗粒物1轻微污染2010-11-794可吸入颗粒物良2010-11-685可吸入颗粒物良2010-11-569可吸入颗粒物良2010-11-468可吸入颗粒物良2010-11-360可吸入颗粒物良2010-11-284可吸入颗粒物良2010-11-176可吸入颗粒物良4.3 建立模型分析影响城市空气污染程度的主要因素是什么 1.通过对10个城市的空气污染情况分析，用excel统计出污染物种类与污染的天数，结果如下表所示：污染种类可吸入颗粒物二氧化硫成都2530杭州2470北京2593上海1672广州1983拉萨930乌鲁木齐24244郑州26224武汉2540西安2880统计结果显示，污染物分为可吸入颗粒和二氧化硫两类。而且对比每个城市的数据，结果显示可吸入颗粒物为首要污染物的天数占绝大多数，远远大于二氧化硫。2. 影响城市空气污染程度的主要因素 首要污染物的来源 可吸入颗粒物的来源有以下五个方面： a.工业污染（燃煤排放和生产中产生的尘及有机物）； b.机动车尾气排放； c.来自道路、施工、裸露地面的扬尘； d.采暖小煤炉及采暖锅炉的民用燃煤排放； e.餐饮、干洗等其他排放。 二氧化硫的来源主要是煤炭的燃烧 通过调查分析，我们认为在一定时期内，每天的工业污染和餐饮干洗等其他排放，即来源a和e可以看作是不变的。我们只需要根据每天污染源数值的变化量，即民用采暖、汽车尾气、扬尘这三项指标进行空气质量的预报监测。3 各个因素的具体分析(1) 民用采暖，分为采暖季和非采暖季两类，得到如下表： 每日平均污染指数采 暖 季非 采 暖 季SO2可吸入颗粒SO2可吸入颗粒63.8104.528.399.2采暖季二氧化硫（SO2）每日平均污染指数是非采暖季的63.8/28.3=225.4%；多出35.5/28.3*100%=125.4%；多出63.8-28.3=35.5个指数；采暖季多出的35.5个指数所对应的就是来自采暖小煤炉及采暖锅炉等民用燃煤排放；这35.5个指数所对应的是占KLW来源的一部分：采暖小煤炉及采暖锅炉的民用燃煤排放，也是占采暖季二氧化硫（SO2）每日污染指数总数的35.5/63.8*100%=55.6%的部分。因此在预报某日空气质量时，应当分采暖季和非采暖季两个不同时间段。在同等或相似的天气条件下，如果设非采暖季二氧化硫（SO2）污染指数为X，可吸入颗粒物（KLW）污染指数为Y，则有：时间段SO2污染指数KLW污染指数非采暖季XY采暖季225.4%X110.8%Y(2)汽车尾气 北京上海等国际大都市的机动车保有量大，机动车尾气排放污染已成为城市空气质量的主要来源之一。 我们从北京市交管局指挥中心得到了一部分反映北京日平均车流量的几组典型数据：2005年线 路工作日日均流量周六日日均流量4月二环路三环路四环路6月二环路三环路四环路9月二环路三环路四环路10月二环路三环路四环路总 流 量因此，北京地区周六日日均汽车总流量/工作日日均汽车总流量为：/*100%=92.1%即周六日每天的汽车尾气排放量是工作日每天的92.1%. 按照月份划分的尾气排放从访问北京市公安交通管理局指挥中心和新闻中心获得的资料来看，北京市每年的3月份、4月份和9月份的车流量较其他几个月份要大，这三个月之间相比，车流量差异不明显；一年中其余的月份之间差异也不明显。 下图的曲线说明2004、2005年北京6-11月份日均车流量的趋势是大体相同的。市公安交通管理局 介绍说，从每年数据统计结果来看，尽管同一个月份日均车流量不同，但是月和月之间的比值，每年基本上是一样的。 2004、2005年6-11月北京二环路日平均车流量及车速曲线 2004、2005年6-11月北京二环路日平均车流量及车速曲线 (该曲线数据来源：北京市公安交通管理局新闻中心)北京市2005年全月份日均车流量数据表(部分月份)2005年二环路三环路四环路总计1月3月4月5月6月8月9月10月(上表数字来源：北京市公安交通管理局指挥中心)2005年1、5、6、8、10月五个月份的日平均车流量为： （+）/5 = .62005年3、4、9月三个月份的日平均车流量为： （+）/3 = .3 由上述两个算式的结果进一步得出：.6 / .3 100% = 93.7% 市公安交通管理局介绍，北京市每年2