数模空气质量分析预测.doc

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对数据进行统计，并作出折线图表，用线性图像相关法对其进行分析，得出相应结论，并就两种方法的出的结论进行比较分析，就两种方法进行评价，突出更有效方法的科学性、时效性及便民性。首先分析天气所包括的因素，再对不同因素进行分开探究。运用控制变量法将其错综复杂的关系分解，清晰地判别其与环境质量的关系。并运用折线图、饼状图等以图表形式，形象、具体、直观地反映出其相关关系。关键词：spss matlab AQI空气质量指数 SO2 NO2 CO PM2.5 PM10 一 问题重述空气质量的好坏直接反映了空气的污染程度，它是依据空气中污染物浓度的高低来判断的，所以控制污染物的排放是改善空气质量的根本措施。空气污染的污染物主要有二氧化硫（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.17912014/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.80.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-483471162014-49342842014-4-151381041920.8246632014-4-16110721630.947862014-48342732014-47328412014-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作为参考序列 (k)（k=1,230）， （1）其他因素作为比较因素序列(k)（i=1,5;k=1,30）。各因数按照下面的式子进行初值化处理，得各标准化序列（i=1,5;k=1,30） i=1,2,.5 （2）（2） 根据求出差 （3）得出：=0； =1.6664（3）令，=1/5 （i=15） 利用公式 （4）按照关联分析的要求及步骤，很容易计算出各污染物的关联度，如下表所示：表4：各污染物的关联度污染物SO2NO2PM10COPM2.5关联度r0.7455 0.76200.79230.73090.8517由此可知：各因素关联度大小排序：由此可知：环境空气质量的变化规律与PM2.5、PM10等相似。图4.6根据图4.6 AQI数据分析图很容易看出AQI指数随着PM2.5、PM10、NO2、SO2、CO的变化而变化，而且根据spss曲线拟合的效果，AQI指数在13年11月迅速上升，到12月时达到顶峰。而14年开始AQI值降低，并且降低速率也是逐渐减慢，在4月份以后就趋于缓和状态，在120上下徘徊。4.3 空气质量预测 要预测未来七月份的空气质量，我们首先算出每个月份AQI的平均值，然后用MATLAB的CURVE FITING TOOL的高斯分布类型进行数据拟合与预测图如下图4.7： 图4.7拟合分析如下：Fit found when optimization terminated:General model Gauss2:f(x) = a1*exp(-(x-b1)/c1)2) + a2*exp(-(x-b2)/c2)2)Coefficients (with 95% confidence bounds):a1 = 165.5 (-3.331e+004, 3.364e+004)b1 = 2.406 (-25.48, 30.29)c1 = 0.5852 (-85.9, 87.07)a2 = 150.9 (-482.7, 784.5)b2 = 1.866 (-129.8, 133.6)c2 = 11.84 (-154.7, 178.4)Goodness of fit:SSE: 323.5R-square: 0.9797Adjusted R-square: 0.9289RMSE: 12.72 图4.7以上图的横坐标18代表2013/112014/6,通过以上分析可以看出拟合度很好，并带入9数值我们可以预测未来七月AQI平均值为125.0352。 AQI急剧减少是不现实的，因此这个模型只适合预测短时间内的AQI，对于长期的预测会出现较大的误差。对于这个问题，我们只预测未来七月份的AQI值。 最近几个月的AQI均值约为162，约属于二级范围，因此空气质量为中等。这年开始空气质量波动起伏较大，四月中旬后趋于平缓，基本与PM10 的变化趋势一致。其中变化较平缓，且呈下降趋势，由前期的平均80到90降到后来的30到40；变化起伏比较大，但整体变化不大平均在到90之间；PM10 的起伏较大，有过几次大的波动，前期尤为明显，后期趋缓，PM2.5波动较大，且数值也较大。PM10是影响空气质量的主要因素，但同时PM2.5也成为影响空气质量的重要原因，并影响了环境污染程度的评测。五 问题二的解答5.1分析三个月的空气质量值图5.1（1） 由图中的折线图，可以清晰地看出：四月份的AQI值普遍比较低。（2） 四月份以后，气温较高，无论是企业还是家庭，燃料用量会减少，而每年11月初到下一年3月初，气温较低，基本处在供暖时段，对燃料的消耗会增大。另外，季节的不同还会导致大气压、风速以及湿度等的不同，而大气压和风速对气态污染物（，）有显著的影响，风速对PM10有显著影响。温度和湿度对有微弱影响。 得出结论季节的差异会通过其对气温、大气压、风速、湿度等直接或通过影响人们生产生活间接影响空气质量。5.2天气对空气质量的影响小雨多云晴阴良553轻度污染210225中度污染1061重度污染284严重污染341表5.1从图表可以看出，天气晴朗的时候污染程度比较严重，在枣庄地区下雨的时候对空气污染有一定的缓解，AQI呈下降的趋势，说明雨天是空气质量一个影响因素。因为降雨对于空气污染物能够起到清除和冲刷的作用。在风利比较大的季节如秋冬季节，AQI的数值呈上升趋势，且风力的大小与空气质量成负相关的关系，究其原因，我们认为由于枣庄的地理环境所决定的，枣庄空气比较干燥，有沙尘，所以风速使空气质量下降。分析温度与空气质量的关系，我们发现之间并没有明显的关系。综上所诉，季节和天气对空气质量有一定的影响，其中风速和降雨量是重要影响因素，对于枣庄而言，风速与空气质量成负相关，降雨成正相关。5.3风向对空气质量的影响图5.2根据图5.2我们可以知道，当无持续风时轻度污染占了一半,重度污染次之，然后是中度污染，良只占百分之九，所以无持续风向时污染时空气相对好点。图5.3根据图5.3，我们可以知道，南风的时候虽然轻度污染占得比重是最多，但是相比较无持续风向时，污染程度还是低一点。重度污染次之，然而南风时良的时间是最多的，所以当是南风时空气质量最好。 图5.4根据图5.4,我们可以知道，北风的时候轻度污染没有前两个比重大，但是重度污染占得比重是三个中最大的，而且严重污染程度也是比较严重的，所以当是北风时空气质量最差。5.3.1数据分析 图表分别对北风、东南风、南风、无持续风四种情况下的空气质量进行了比例统计，从中可以很清晰地看出，南风条件下空气质量较为好，无持续风时次之，而北风时，空气质量较差。5.3.2得出结论 枣庄市出现较多的北风会对环境质量造成不良影响，而南风情况下环境质量一般较好。5.4温度对空气质量的影响图5.7 11月是冬季温度相对来说比较低的月份，在不考虑其他条件的的影响下，随着温度的升高，AQI的指数会不断升高。这是由于冬季天气比较干燥，温度越高空气中的污染物越多。图5.8 1月份是春季相对来说温度中等，在不考虑其他因素的影响下，AQI也随着温度的增加而增加。 图5.9 4月份是春夏交替的时间，温度普遍适中，温度越低空气中的湿度就越高，空气中的污染物就越多，所以AQI就越高。 5.4.1数据分析 根据以上三个折线图，我们分别分析了11月、1月、4月的温度和AQI