《大气污染评价与预报分析模型.docx》由会员分享,可在线阅读,更多相关《大气污染评价与预报分析模型.docx(11页珍藏版)》请在优知文库上搜索。
1、大气污染评价与预报分析模型摘要本文对空气质量的评价及污染预报问题进展了分析.运用层次分析法依据处理后的数据对六个城市的空气质量进展了具体细致的排序;对2010年9月15日至9月21日的各项污染物浓度、各气象参数运用一元多项式回归模型进展了预测;就气象参数所属城市问题及污染物浓度与其的关系建设了相关性分析模型和多元线性回归模型;最后.根据建模过程和结果.我们对相关部门提出了几个具体的建议。通过将数据附件所给有效数据.即日污染物浓度.转化为对应的月污染物浓度的均值.根据各城市月均污染浓度做出其随时间的走势折线图.分析了各个城市SO?、NO2、PMlO之间的特点。我们拟根据APl指数值.以二级达标次
2、数为准.对各城市之间的空气质量进展排名.但由于依据APl的区分空气质量等级时灵敏度较低.故采用了层次分析法对空气质量进展排名。由于我们采用了全部数据进展排名.而E、F数据较少.故只对ABCD进展了排名。依据层次分析法得出的排名为:A、B、D、Co为了准确预测各城市短期内的数据.本文选用一元多项式回归模型。对2010年的数据进展分析整理.依据回归模型得出其与时间的关系.得出预测值.并得出其置信度为95舟的置信区间.结果显示模型的预测效果尚能承受.能够对所要预测数据进展预测。但由于F城市数据缺失.根据假设做了合理的定性分析.并未对其进展定量预测。分析空气质量与气象参数之间的关系时.首先根据数据完整
3、性.气象参数应只属于其中一个城市.排除了D、E、F的可能性.再根据相关性分析的方法.确定了气象参数属于A城市。根据污染物与气象参数之间的因果关系.建设了多元线性回归模型.由于季节对污染物的浓度存在影响.分季节得出各污染物与各气象参数之间的相关系数.定性分析该相关系数.得出污染物与气象参数之间的关系。最后对该系数的理论与实际意义做了检验。根据以上分析及结果.确定局部与空气质量控制相关的部门.针对其职能提出了诚恳建议。关键词:API评价模型层次分析一元多项式回归模型问题重述大气是指包围在地球外围的空气层.是地球自然环境的重要组成局部之一。人类生活在大气里.干净大气是人类赖于存在的必要条件。随着地球
4、上人口的急剧增加.人类经济增长的急速增大.地球上的大气污染日趋严重.其影响也日趋深刻。因此.加强大气质量的监测和预报是非常必要。目前对大气质量的监测主要是监测大气中5。2、NO2、悬浮颗粒物主要为PMI0等的浓度.研究说明.城市空气质量好坏与季节及气象条件的关系十分密切。现有城市A、B、C、D、E、F从2003年3月1日至2010年9月14日测量的污染物含量及气象参数的数据。1 .找出各个城市SO2、No2、PM10之间的特点.并将几个城市的空气质量进展排序。2 .对未来一周即2010年9月15日至9月21日各个城市的SO?、No?、PM10以及各气象参数作出预测。3 .分析空气质量与季节、气
5、象参数之间的关系。4 .就空气质量的控制对相关部门提出你的建议。二、问题分析问题一寻找各城市SO2、NO2、PMlo的特点.最直观的方法就是作图.把各城市的三种污染物浓度做到一张图中进展比照分析可较容易的得到其特点.而排序题目中给出的是三种污染物浓度.必须先用一个指标将它们统一起来综合的对城市的空气质量进展评价.用同一个指标进展排序。问题二是依据所给的2010年1月1日至9月14日的数据.预测2010年9月15日至9月21日各个城市的SO2、N02、PM10以及各气象参数.预测的时期较短.数目多.选择时间序列进展预测。同时将数据序列选取为2010年所有测量日的数据。根据给定的数据.利用一元多项
6、式回归.求得回归模型.从而预测出需要的7项数据.并进展预测误差估计。问题三是研究SO2、No2、PMIO的浓度与气象参数之间的关系。首先应对六个城市的So2、No2、PMIO的浓度同气象参数进展相关性检验.以找出气象参数与SQ、No2、PMIo的浓度的对应关系。又由于城市空气质量好坏与季节及气象条件的关系十分密切。故分季节对SO?、NO2、PMIO的浓度与气象参数进展回归分析.并进展检验。问题四则是通过对气象参数.季节等因素的考量.分析得出提高空气质量减小So2、No2、PMlO的浓度的方法.从而给有关部门提出合理性的建议。三、模型假设1 .各城市的检测时间具有随机性;2 .F城市的开展是平衡
7、开展.政府对环境治理干预较小.即F城市的环境不会出现强烈波动;3 .数据附件所给的各个城市的污染物浓度及气象参数不随测量地点的改变而改变.且有效数据都准确可靠。同时不考虑人为因素.检测仪器准确度不同的影响具有统计、预测意义;4 .在对预测期内即2010年9月15日至9月21日时间段内.各个城市不会出现重大空气污染事故.或环境不会变好;5 .城市空气质量好坏与气象条件的关系十分密切.与其他因素关系不大。四、模型的建设和求解(一) 找出各个城市SO2、No2、PMIO之间的特点.并将几个城市的空气质量进展排序数据处理:剔除1.第86项2005年的.与前后年份不符很明显是错误数据;2.第986项气压
8、值的明显错误。大气污染程度与空气中有害气体的含量有关.根据题目要求我们只考虑SO2、NO2、PMlo的影响.由于数据量大我们按月为周期计算出了每个城市每种污染物浓度的月平均值.据此画出折线图从而观察分析其特点。对于排序问题.我们参考相关资料采用了判断大气污染的空气污染指数(API.API的计算依据为表1污染指数污染物浓度毫克/立方米)APIS02C日均值N02(日均值PM10日均值5002.6200.9400.6004002.1000.7500.5003001.6000.5650.4202000.2500.1500.3501000.1500.1000.150500.0500.0500.050有
9、计算公式:其中.I为某污染物的污染指数.c为该污染物的污染浓度。C人与CJ、分别为上表中最贴近C值的两个限值.C大为大于C的限值.C小为小于C的限值.同样./大与/小也是最靠近C值的两个限值。得出每种污染物的月平均污染指数。取三种污染物中API中的最大值作为该市的月平均污染指数。等级判断标准:空气污染指数APl空气质量状况对安康的影响速议采取的措施0-50优可正常活动。51-100良101-150轻微污染易感人群病症有轻度加剧.安康人群出现刺激病症。心脏病和呼吸系统疾病患者应减少体力消耗和户外活动。151-200轻度污染201-250中度污染心脏病和肺病患者病症显著加剧.运动耐受力降低.安康人
10、群中普遍出现病症。老年人和心脏病、肺病患者应当停留在室内.并减少体力活动。251-300中度重污染300重污染安康人运动耐受力降低.有明显强烈病症.提前出现某些疾病。老年人和病人应当留在室内.防止体力消耗.一般人群应防止户外活动得到APl值。但假设单纯的采用APl值进展排序.API值的分类太过粗糙所以我们采用层次分析法对来处理.最终得到方案层对目标层权重.进展排序。建设以各个城市为方案层.空气质量为准侧层.空气质量排名为目标层的层次分析模型,由于E、F市数据过少.这里我们只考虑A、B、C、D市的排名.层次图如下根据两两比照法建设准则层对目标层的判断矩阵A根据等级的19比照尺度定性的两两比照得到
11、判断矩阵A计算出权重向量E建设方案层对准则层的判断矩阵B根据我们求得的各城市每月空气污染指数统计得表3各城市污染指数统计IIlIllIVVA540151OB3312520C11232141D1301841从而再次利用两两比照法分别得到各城市对不同空气等级的判断矩阵Bi.根据特征根法确定权重向量F=w1,w2,w3,w4,w5通过一致性检验后.由怕ETF得到最终方案层对目标层的权重.它表示的是受污染程度的权重.所以权重越大.受污染程度越高.空气质量越差。求解各城市直线图.以以以下图中蓝色的实线为S02浓度折线图.红色的虚线为PM10的浓度折线图.黑色的点线为N02浓度折线图。A城市各污染物月平均
12、浓度分析0.71I110.6-0.5从图中可看出A城市中S02指数除前几个月有明显上升外.有明显的下降趋势.大局部低于PM10的指数.且变化与PM10有一定的一致性。PM10指数有较明显的振荡.但总体变化不是很大。N02指数普遍最低.较平稳。8城市各污染物月平均浓度分析0.35I1111B城市的So2、PMIO变化很是相似.都以较大的幅度振荡且总趋势是下降的。S02的振荡更明显.N02指数均低于前两种污染物指数.较平稳。与A、B一样C市中SO2、PM10的振荡具有一致性.但C市的PM10普遍比S02高.且在2025月左右有十清楚显的上升.S02振荡幅度不是很大.较平稳。城市各污染物月平均浓度分
13、析O45111111D市中S02和PM10又出现大幅度的振荡总体趋于降低的趋势N02指数还是最低最平稳。0.18E城市各污染物月平均浓度分析E币中PM10指数高于另两种指数.在5月和22月左右各出现一次大幅度振荡.且下降趋势不明显.S02指数在PM10下振荡下降.N02出项小幅的振荡但总体平稳。F城市各污染物月平均浓度分析 0.110.10908070605040302Ua0.0.0.0.0.0.月份F市数据太少.仅有的三个月数据看较平稳.但数据太少缺乏信。1根据两两比照的到准则层对目标层的判断矩阵AI1/21/31/41/52I2/32/42/533/2I3/43/544/24/3I4/55
14、5/25/35/4I计算的X(max)=5.00.CI=O.CR1=0,一致性通过.所以得到权向量E=(0.1238,0.1599,0.3857,0.3317)r又根据表3得至U方案层对准则层的判断矩阵B-15/35/115-140/3240/234/33/513/11331/40131/2331/30B1=11/511/3111B2=23/4023/31123/301/51/31/111_3/430/3130/231-115/2515/21151811/21/41425/15125/2125/18B4=211/21/221/1521/25121/18421118/1518/2518/2114211由于第五等级的数据较少.可以定性的判断其权重运用Matlab编程算出其权向量吗.最大特征根乙和一致性指标C,结果列入下表表4123450.2500.32260.19000.098000.1500.25000.31640.196000.5500.18550.26580.