大气污染预报问题

摘要本题是研究大气污染预报问题，根据科学方法监测出六个城市的SO2、NO2、PM10三种污染物的浓度值以及一组气象参数值，分析得出大气污染物的特点并对大气质量作出预测，根据其特点得到对于大气治理的建议。问题一：利用EXCEL求出各城市的API（空气污染指数）月平均值并画出折线图，分析图表得到三种污染物的特点。再利用层次分析法得到六个城市大气质量的排名，建立层次模型时，根据国家标准对空气质量进行了五个档次的划分作为因素层，要排名的城市作为对象层，计算出因素层的权重比例，质量越差得到权重越大，接着按因素层的分类分别对对象层进行权值计算，最后得出总排名。问题二：由于污染物浓度变化波动比较大，所以采用自适应滤波法进行预测，得到2010-9-15至2010-9-21各城市SO2、NO2、PM10的值，对于各气象参数的值，先画出其随时间的变化图，看出其数值大概在某条线左右变化，所以采用移动平均法对其进行了预测。问题三：利用典型相关性分析，得到所给气象参数数据来源于C城市，采用多元线性回归模型并考虑到季节对于大气质量的影响，分别得到春季、冬季的SO2、NO2、PM10与气象参数之间的关系。问题四：在前几问的基础上，总结得出大气治理的一些方法，给相关部门以建议。关键词：层次分析法典型相关性多元线性回归一、问题重述众所周知，洁净的空气是人类生存的必要条件之一。世界卫生组织和联合国环境组织发表的一份报告说：“空气污染已成为全世界城市居民生活中一个无法逃避的现实。”工业文明和城市发展，在为人类创造巨大财富的同时，也把数十亿吨计的废气和废物排入大气之中，人类赖以生存的大气圈却成了空中垃圾库和毒气库。因此，对于大气污染物的检测预报以及防治是我们目前的重要任务。目前对大气质量的监测主要是监测大气中2SO、2NO、悬浮颗粒物（主要为PM10）等的浓度，在解决大气污染物的防治问题中，我们要考虑一下两点：（1）、三种污染物2SO、2NO、10PM的特点；（2）、空气质量与气象参数以及季节之间的关系。通过已经观测得到的数据，我们可以预测之后的大气，并且通过一定的分析，对于我国的大气防治提出一些建议。二、模型基本假设本题中假设如下：1、题目所给的六个城市的污染物含量及各城市的气象参数等数据都准确可靠；2、API指标真实可靠，所给数据具有参考统计意义；3、API的月平均值能较好的代表该月空气质量，具有比较意义。三、符号说明API空气污染指数I某污染物的污染指数C该污染物的浓度CC大小、在API分区表中最接近C值得两个值II大小、在API分区表中最接近I值得两个值CI层次模型的一致性指标RI随即一致性指标CR一次性比率p大气压力t温度w湿度v风速四、问题分析第一问是找出A至F城市2SO、2NO、10PM之间的特点，并将几个城市的空气质量进行排序。在分析三种污染物之间的特点时，通过查阅资料，找到API（空气质量污染指数标准）,由此计算出每个城市各项指标的API月平均值，再利用EXCEL进行统计分析，得出各个城市三种污染物的特点。在进行几个城市空气质量的排序时，考虑采用层次分析法进行求解。第二问是预测未来一周即2010-9-15至2010-9-21各城市2SO、2NO、10PM以及各象参数的值。考虑采用自适应滤波法进行预测。自适应滤波法与移动平均法、指数平滑法一样，也是以时间序列的历史观测值进行某种加权平均来预测的，但是移动平均和指数平滑法预测的数据列要求具有较为平稳的趋势变化，而利用excel对所给数据进行分析发现数据随着时间的变化呈现波动变化，虽然各年之间数据呈现随年份增加而下降的趋势，但是本文只要预测最后一年7天的数据，由于前几年的历史数据对最后7天的预测没任何意义，因此也就不用考虑数据的趋势项，只需根据最后三个月数据的波动变化来建立模型即可。它要寻找一组“最佳”的权数，其办法是先用一组给定的权数来计算一个预测值，然后计算预测误差，再根据预测误差调整权数以减少误差。这样反复进行，直至找出一组“最佳”权数，使误差减少到最低限度。第三问是分析空气质量与气象参数之间的关系。由于所给气象参数数据只有一组，所以首先应进行典型相关性分析，得出所给数据来源于哪个城市。在研究空气质量与气象参数之间的关系时，由于气象参数有四个，所以考虑采用多元回归线性模型找出各个空气质量与气象参数之间的关系。近一步考虑，空气质量与季节有关，所以分别建立春季和冬季的多元线性回归模型。第四问是就空气质量的控制对相关部门提出建议。该问题是要在前三问的基础上，找出污染物的影响因素，针对其影响因素提出合理的建议。五、模型建立及求解5.1问题一的求解5.1.1三种污染物2SO、2NO、10PM的特点分析在解决A至F城市2SO、2NO、10PM间的特点这部分问题是，通过查阅资料，可以找到API（空气质量污染指数标准），由此计算每个城市各项指标的API月平均值，对各项数值进行比较，从而可分析得出各个城市三项指标的特点。表1空气污染指数对应的污染物浓度限值污染指数污染物浓度（毫克/立方米）APISO2（日均值）NO2（日均值）PM10（日均值）CO（小时均值）O3（小时均值）500.0500.0800.05050.1201000.1500.1200.150100.2002000.8000.2800.350600.4003001.6000.5650.420900.8004002.1000.7500.5001201.0005002.6200.9400.6001501.200设I为某污染物的污染指数，C为该污染物的浓度。则：-()+-IIICCICC大小小小大小（1）由表一及（1）式可计算出各城市各项指标的污染指数，由于数据量过大，我们用各月的API平均值代表该城市该项指标的污染指数，得到图表如下所示：通过对上图的分析，可以得到以下结果：A、B、C、D、E五个城市2SO、2NO、10PM等污染物浓度均呈现出波动性并且有缓慢下降趋势。A城市与B城市的2NO、10PM两种污染物的浓度变化趋势比较相近。但A城市2SO在03年下半年至04年上半年的浓度值非常大，说明这段时间由于某些因素使得2SO排放较多，但在之后的几年其排放量明显下降，说明该城市采取了某些积极措施进行改进。B城市2SO浓度呈现一定的波动性，说明该城市有一些周期性的污染源。C、D两城市三种污染物浓度的变化趋势基本一致，说明两城市的地理位置以及发展模式都比较相近。其中，在06年时两城市的10PM你浓度非常大，说明这段时间的发展模式使得10PM浓度过高，但之后其浓度有迅速下降，说明该城市的管理者能迅速采取措施来控制10PM的浓度。E城市三种污染物的浓度变化比较平稳，说明该城市的环境质量比较好。F城市由于其数据量过于少，无法看出其变化趋势，但从所给数据来看，该城市的污染物的浓度变化比较平稳。考虑到季节对于大气污染物的浓度随着季节变化，所以，利用EXCEL统计出四个城市在04年至09年取暖季以及非取暖季三种污染物浓度的平均值的比值，结果如下表所示：A地B地SO2NO2PM10SO2NO2PM1004年1.7880.8411.2632.1341.1211.16106年1.6661.0480.7623.3701.3650.88107年0.9710.8781.9792.4421.0071.42208年1.6030.7281.5841.8621.4211.34709年1.5171.0561.5352.7771.2171.639C地D地SO2NO2PM10SO2NO2PM1004年1.7521.2361.0322.2341.0861.24006年1.1371.4790.5662.9481.6390.63407年2.0221.1581.3352.3061.6971.67508年1.9460.6391.4181.6501.1501.47209年3.5591.3681.5565.1601.6661.591表中数据反映出的总体结果是四个城市取暖季的各项污染指数均大于非取暖季的，只有个别年份取暖季的某个污染指数小于非取暖季的；B,D两地取暖季SO2污染明显增强，说明这两地冬季取暖主要以煤炭等为主且使用量较大，且随着年份增加比例逐渐增大，说明B,D两地属于普通居住人口较多的城市；对于A地两季的污染物浓度比相对来说较小，说明A地两季污染物产生情况相差不大，可以反映出A地属于重工业型城市。5.1.2A、B、C、D、E、F空气质量排序由于表中所给观测数据中E,F两城市没有交叉项，所以本文分别对A、B、C、D、E和A、B、C、D、F进行排序，对城市A、B、C、D、E排序取其观测数据的交叉项2007-9-1号的数据到2010-9-14号，计算出每天的API值，由于数据量太大，考虑求出这个时间段内API的月平均值来反映对应城市的空气质量。对A、B、C、D、F排序取其数据交叉项2004-9-1号到2004-12-27号的数据，2004-9-1号F城市的观测值全部为0，与实际情况不符，为坏值剔去，计算出2004-9-2到2004-12-27号每天各个城市的API值用以反映对应城市的空气质量。1、将研究目标（Z）、因素（P）、对象（C）按相关关系分成目标层Z、准则层P、对象层C。层次结构图如图所示：2、给出空气质量一级，二级，三级，四级，五级两两成对比较的判断矩阵P污染级别一级二级三级四级五级12345根据上图得出如下例两两成对比较的判断矩阵P一级二级三级四级五级权重一级１１／２１／３１／４１／５0.0667二级２１２／３２／４２／５0.1333三级３３／２１３／４３／５0.2四级４４／２４／３１４／５0.2697五级５５／２５／３５／４１0.3333由表中数据经计算可得：λ(max)=5，CI=0.00，RI=1.12，CR=0.000.1。因为CR=0.000.1，所以此排序有满意的一致性。3、给出对象层对准则层的各个因素的判断矩阵并进行分析。a、对A、B、C、D、E排序我们利用了EXCEL的相关知识求解出了这31个月中各个城市的空气污染指数：A、B、C、D、E五个城市的空气污染指数统计（单位：月）优良轻度污染中度污染重Ａ921100Ｂ1020100Ｃ524200Ｄ521500Ｅ424300根据上表中的数据可以仿照计算空气质量等级的方法，计算出当空气质量为优时五个城市所占的权重大小：ABCDE权重A19/109/59/59/40.2727B10/9110/510/510/40.3030C5/95/10115/40.1515D5/95/10115/40.1515E4/94/104/54/510.1213由表中数据经计算可得:λ(max)=5,CI=0,RI=1.12，CR=0.000.1所以此排序有比较满意的一致性。同理，可以得出空气质量为良和轻度污染时各个城市的不同权重：空气质量为“良”时五个城市对应的权重表城市ABCDE权重0.19090.18180.21820.19090.2182由表中数据经计算可得:λ(max)=5,CI=0,RI=1.12，CR=0.000.1所以此排序有比较满意的一致性。空气质量为“轻度污染”时五个城市对应的权重表城市ABCDE权重0.08330.08330.16670.41670.25由表中数据经计算可得:λ(max)=5,CI=0,RI=1.12，CR=0.000.1所以此排序有比较满意的一致性。由于中度污染和重污染的指数为零，所以可以不用考虑。最后的出的层次总排序的结果如下：准则优良轻度污染中度污染重污染总排序权值准则层权值0.06670.13330.20.26970.3333方案层单排序权值A0.27270.19090.0833000.0603B0.30300.18180.0833000.0610C0.15150.21820.1667000.0725D0.15150.19090.4167000.1189E0.12130.21820.25000.0872根据5个城市的总权重值进行从小到大依次排序，空气污染严重程度由小到大的排名如下：A、B、C、