数模空气质量分析预测

2010高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：B我们的参赛报名号为（如果赛区设置报名号的话）：031007所属学校（请填写完整的全名）：枣庄学院参赛队员(打印并签名)：1.邱立箭2.崔浩3.张微微指导教师或指导教师组负责人(打印并签名)：日期：2014年6月21日赛区评阅编号（由赛区组委会评阅前进行编号）：2010高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：空气质量评价与预测摘要本文对枣庄市区的空气质量进行了评价，并选出了主要的污染物进行研究分析，运用综合指数评价法和回归分析等方法对其空气质量进行分析，综合各种因素我们建立了如下模型。对于问题一，我们用spss画出污染物SO2、NO2、CO、PM2.5、PM10的时间分布图，然后通过这些图，我们分别对各污染物随时间的变化特征进行分析总结，我们建立模型，进行各污染物的变化规律对环境的关联度分析。然后用matlat拟合工具箱的高斯分布类型画出每月AQI的平均分布图，并根据拟合的高斯公式带入数据预测七月份的AQI平均值。对于问题二，我们用excel、spss和matlab对数据进行统计，并作出折线图表，用线性图像相关法对其进行分析，得出相应结论，并就两种方法的出的结论进行比较分析，就两种方法进行评价，突出更有效方法的科学性、时效性及便民性。首先分析天气所包括的因素，再对不同因素进行分开探究。运用控制变量法将其错综复杂的关系分解，清晰地判别其与环境质量的关系。并运用折线图、饼状图等以图表形式，形象、具体、直观地反映出其相关关系。关键词：spssmatlabAQI空气质量指数SO2NO2COPM2.5PM10一问题重述空气质量的好坏直接反映了空气的污染程度，它是依据空气中污染物浓度的高低来判断的，所以控制污染物的排放是改善空气质量的根本措施。空气污染的污染物主要有二氧化硫（SO2）、二氧化氮（NO2）、可吸入悬浮颗粒物\浮尘（PM10）一氧化碳（CO），大气中直径小于或等于2.5微米的颗粒物（PM2.5）等等。目前，城市空气质量污染指数的分级标准是根据空气质量指数（AQI）的取值界定的，空气污染指数指常规监测的几种空气污染物浓度简化成为单一的概念性指数值形式，并分级表征空气污染程度和空气质量状况。二模型假设1．忽略空气中在指标控制范围内的其他污染物；2．假设所给数据均有效；3．假设在较近一段时间内，不发生重大工业事故；4．枣庄未来一段时间内，政府没有出台关于大规模工业的迁入迁出政策；5.表格中已有的数据具有权威性，值得相信，具有使用价值。6.空气质量相同等级的污染程度相同。7.假设该市各种影响空气质量的软因素（如工业发展，人口数量）保持平稳变化。8.不考虑突发事件即人为因素（如工业事故）造成的空气质量突变。9.假设各种因素对环境的影响最终主要表现在可吸入颗粒、二氧化硫、一氧化碳、二氧化氮上，不考虑其他随机因素的影响。三符号说明：关联分析中的比较因素：初值化处理后的标准化序列W:权向量四问题一的解答4.1各污染物变化特征根据所给的表一，我们经过计算得出空气质量及污染指数均值，用MATLAB对各种污染物的指数进行分析，得出相应的走势，并根据表和图对相应的污染物进行特征分析，最后对未来7月份的均值进行近似的预测。表4.1空气质量指数对应的污染物浓度限值空气质量分指数（AQI）污染物项目浓度限值二氧化硫（SO2）24小时平均/（μg/m3）二氧化氮（NO2）24小时平均/（μg/m3）一氧化碳（CO）24小时平均/（mg/m3）颗粒物（PM10）24小时平均/（μg/m3）颗粒物（PM2.5）24小时平均/（μg/m3）5050402.05035100150804.0150751504751801425011520080028024350150300160056536420250400210075048500350500262094060600500表4.2空气质量指数范围及相应的空气质量类别空气质量指数空气质量指数级别（状况）及表示颜色对健康影响情况建议采取的措施0~50一级（优）空气质量令人满意，基本无空气污染各类人群可正常活动51~100二级（良）空气质量可接受，但某些污染物可能对极少数异常敏感人群健康有较弱影响极少数异常敏感人群应减少户外活动101~150三级（轻度污染）易感人群症状有轻度加剧，健康人群出现刺激症状儿童、老年人及心脏病、呼吸系统疾病患者应减少长时间、高强度的户外锻炼151~200四级（中度污染）进一步加剧易感人群症状，可能对健康人群心脏、呼吸系统有影响儿童、老年人及心脏病、呼吸系统疾病患者避免长时间、高强度的户外锻炼，一般人群适量减少户外运动201~300五级（重度污染）心脏病和肺病患者症状显著加剧，运动耐受力降低，健康人群普遍出现症状儿童、老年人及心脏病、肺病患者应停留在室内，停止户外运动，一般人群减少户外运动300六级（严重污染）健康人群运动耐受力降低，有明显强烈症状，提前出现某些疾病儿童、老年人和病人应停留在室内，避免体力消耗，一般人群避免户外活动表4.3空气质量及污染指数均值日期AQI指数均值PM2.5均值PM10均值CO均值NO2均值SO2均值2013/11150.1112.8191.21.2160.47108.92013/12253.3201.1326.01.7573.7159.92014/1208.5158279.51.43457.591.032014/2143.25106.81781.0136.670.12014/3153.6105.1212.90.9743.179．12014/4123.984.3174.50.8840.9712014/5119.970.11760.7738.663.82014/6120.782.9145.10.8135.538.1由表4.3空气质量及污染指数均值与上面的标准对比可以看出，各月的空气质量都处于污染中。十一、二月、四月、五月、六月属于轻度污染，而十二月和一月达到了重度污染。从表中我们也可以看出，在这些污染物中主要的污染物是PM2.5。图4.1PM2.5：主要来源是汽车尾气的排放、集中供暖产生的废气。判断雾霾天气的主要指标。由于可吸入颗粒物的浓度与绿化植被覆盖率、最近降雨量相关，所以在4月几场大的降雨量后，PM2.5值降低了25%。预计，每年的PM2.5最小值与当年的降雨量相对应。由图4.1不难看出PM2.5从2013年十一月到2014年一月上升的，2014年1月达到这几个月的最大值，然后一直到2014年五月都在下降，五月的时候PM2.5达到最小值，但是到六月份的时候又开始有上升的趋势，但是整体逐渐成下降趋势。PM10：十一月份到十二月份是增高的，也达到了这几个月的最大值，到二月都是呈下降的虽然在三月份的时候上升了，但是在之后都是下降趋势的，而且在六月达到最低，说明空气质量在逐渐变好，整体稳定在100到200之间。PM10主要来源是汽车尾气的排放、不合格烟尘排放。本年1月份至4月份，PM10污染浓度逐渐降低。由于可吸入颗粒物的浓度与绿化植被覆盖率、最近降雨量相关，所以在几场大的降雨量后，PM10值降低了许多。预计，每年的PM10最小值与当年的降雨量相对应。NO2：主要来源是汽车尾气的排放。分析NO2的趋势线可知，2013年11月份至2014年6月份，相对于趋于平缓趋势。但在2014年1月份至6月份NO2污染浓度逐渐降低。其随着枣庄市车辆密度的增加而增加，呈正相关。SO2：主要来源是集中供暖产生的废气。分析SO2的趋势线可知，本年11月份至12月份，SO2污染浓度逐渐升高，次年1月份开始so2污染浓度逐渐下降，枣庄四月份的SO2污染浓度下降到了90左右。这与采暖期污染源增加和非采暖期污染源减少相对应。预计，每年的SO2污染浓度最大值与当年的最冷月相对应。图4.2CO：主要来源是汽车尾气的排放、集中供暖产生的废气。根据图4.2CO数据分析图，我们看到CO随着时间的推移指数有显著地下降，每一阶段下降的速度也不一样，开始的时候，下降速度最快，也表明了，枣庄市环保局对枣庄市做出了一定的努力，随着后期阶段，CO的指数下降速度也逐渐变慢，一直到4-5月，指数逐渐趋于缓和状态，指数稳定在0.8~0.9之间。4.2空气质量变化规律关联度分析关联分析法确定影响环境空气质量的主要污染物指标,为便于数据处理，我们选取近30天作为样本，数据记录如下表：日期AQI指数PM2.5PM10CoNo2So22014-4-1124881740.8337612014-4-2127951800.7935562014-4-31691172490.8944892014-4-4130752080.8447732014-4-5108401650.5430422014-4-695471380.5935572014-4-7115701780.856742014-4-81421052171.05591012014-4-91531092261.1547862014-4-101861342821.28561072014-4-111561182111.0847602014-4-1280551000.631472014-4-13124921800.83471162014-4-14110791650.9342842014-4-151381041920.8246632014-4-16110721630.947862014-4-17116851600.8342732014-4-18105661290.7328412014-4-1987621150.7137472014-4-201811362561.3259892014-4-212241753171.59691182014-4-222181703021.2143772014-4-23122781600.8333552014-4-2499521170.6427372014-4-2582481030.6222412014-4-265427530.6618272014-4-2784561050.7125562014-4-2784561050.7125562014-4-2879481060.728892014-4-29114741620.91511112014-4-30126791911.245682（1）将AQI作为参考序列20minmax010max**iikirw0x(k)（k=1,2···30），（1）其他因素作为比较因素序列ix(k)（i=1,···5;k=1,··30）。各因数按照下面的式子进行初值化处理，得各标准化序列iyk（i=1,···5;k=1,··30）301,....211,iiiiiXXXXkyi=1,2,.......5（2）（2）根据iyk求出差00iikykyk（3）得出：min=0；max=1.6664（3）令0.5，iw=1/5（i=1···5）利用公式301maxmaxminkoiioiwr（4）按照关联分析的要求及步骤，很容易计算出各污染物的关联度，如下表所示：表4：各