数模灰色预测模型的大气质量

天津农学院承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）天津农学院参赛队员(打印并签名)：1.2.3.指导教师或指导教师组负责人(打印并签名)：日期：2012年08月20日赛区评阅编号（由赛区组委会评阅前进行编号）：天津农学院编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）1基于模糊与灰色模型的大气质量评估摘要近年来，能源短缺和环境污染问题成为世界关注的焦点问题，转变传统高能耗、高污染的经济增长方式，大力推进节能减排，发展以低能耗、低排放为标志的低碳经济，实现可持续发展，正在成为世界各国经济发展的共同选择。本文正是在这样的背景下综合全国31个主要城市，建立模糊综合评价模型，对全国空气质量做出综合评价。利用灰色模型对2012年各城市污染状况做出预测比较，并通过假设在未实行节能减排措施下，历年空气质量指数的变化情况，从而对节能减排措施对环境改善作用做出科学评价。经计算所得历年各城市空气质量等级基本处于二级，但部分地区污染状况较为严重，可能与当地经济发展状况和产业结构有密切关系。在对2012年各城市空气质量预测中应用修正后的灰色模型得出各城市AQI值变化情况，从中发现部分污染严重的城市环境质量有所改善，其余多数城市环境质量处于良好状态且保持稳定，这与实际情况基本符合，说明预测具有一定的合理性。针对问题二我们用时间序列模型对2007-2010年空气质量指数做出预测与实际数据进行对比，发现各城市空气质量指数明显上升，环境趋于恶化。针对问题三和四我们通过对2007-2010年全国2SO、烟尘、工业粉尘三个污染物的排放量与实际排放量相对比，评判出节能减排对大气环境所起的作用。首先我们运用灰色关联模型确定各个污染物之间的权重因子，发现SO2、烟尘、工业粉尘实际权重分别为0.616、0.238、0.146，预测权重为0.659、0.212、0.128，对比发现2SO对大气环境起决定性作用，所以对2SO排放量的控制至关重要，由实际权重和预测权重相对比发现由于节能减排措施的实施，2SO的权重有所下降，说明节能减排对大气环境起到很好的改善作用，政府应该加大推广节能减排措施的力度，从而为市民创造出更优良的生活环境。关键词：灰色模型、时间序列、AQI、模糊综合评价法、灰色关联分析法2一、问题提出环境保护是重大民生问题，随着社会对环境保护的日益重视，人们越来越重视环境的改善，工业革命以来，世界各国尤其是西方国家经济的飞速发展是以大量消耗能源资源为代价的，并且造成了生态环境的日益恶化。节约能源资源，保护生态环境，已成为世界人民的广泛共识。我国从2007年8月起，中央财政开始实施节能减排工作，既是对人类社会发展规律认识的不断深化，也是积极应对全球气候变化的迫切需要。《国民经济和社会发展第十一个五年规划纲要》提出了“十一五”期间单位国内生产总值能耗降低20%左右，主要污染物排放总量减少10%的约束性指标。根据这两个指标，如中国GDP年均增长一成，五年内就需要节能六亿吨标准煤，减排二氧化硫六百二十多万吨、化学需氧量五百七十多万吨。试根据我国近年污染物总量减排和大气环境相关数据，并结合经济发展情况，根据附录中的数据，结合你们收集到的相关资料，建立数学模型，完成以下问题：1、建立模型对全国各省会城市的大气环境质量做出定量的综合评价，并对2012年各地区大气的污染状况进行分析比较。2、假如不采取节能减排，依照过去几年的主要统计数据，对我国大气环境的发展趋势做出预测分析，3、建立模型分析讨论节能减排对大气环境质量改善所起作用。4、建立模型对节能减排实施前后各省会城市大气环境质量改善情况进行科学分析。5、对下一步实施节能减排提出建议。二、基本假设1.假设评价空气质量的各指标间相互作用关系忽略不计；2.假设影响空气质量的因素只有2SO、102PMNO、,其他影响因素忽略不计；3.假设从官方获取的各个省份的指标的统计数据信息真实可靠；4.假设这几年中没有出现重大的自然灾害，及经济发展相对稳；三、符号说明AQI空气质量指数IAQI空气质量分指数P污染项目PIAQI污染物项目P的空气质量PC污染物项目P的质量浓度值HiBP与PC相近的污染物浓度限值的高位值3LoBP与PC相近的污染物浓度限值的低位值HiIAQI与HiBP对应的空气质量分指数LoIAQI与LoBP对应的空气质量分指数t各年份对应的编号)()0(tx综合空气质量指数的参考数列)()1(txAGO生成数列四、问题分析针对问题一：建立模型对全国各省会城市的大气环境质量做出定量的综合评价并对2012年各地区大气的污染状况进行分析比较,通过查阅《中国统计年鉴》[1]获得历年2SO、102PMNO、浓度数据（见附表1）,根据附件中提供的空气质量分指数的计算公式[2]，得出2003-2010年各城市AQI值，从而对全国各省会城市的大气环境质量做出定量的综合评价，近而预测出2012年各地区大气环境质量指数。针对问题二：假如不采取节能减排，依照过去几年的主要统计数据，对我国大气环境的发展趋势做出预测分析，通过观察附录中的数据，根据国家规定的从2007年开始节能减排的要求，利用2003-2006年的AQI值进行预测，进而对我国大气环境的发展趋势做出预测分析。针对问题三和问题四：建立模型分析讨论节能减排对大气环境质量改善所起作用和对节能减排实施前后各省会城市大气环境质量改善情况进行科学分析，根据2007-2010年的实际AQI值与预测所得的AQI值通过作图对比，分析得出节能减排对大气环境质量改善所起作用。五、模型的建立与求解5.1各城市大气环境分析5.1.1问题一的分析为对全国各省会城市的大气环境质量做出定量的综合评价，我们采用模糊综合评价法。通过查阅《中国统计年鉴》，获得所需数据（见附表1），确定因素集、权重集、评价集和综合评价模型，根据附件中所提供的空气质量等级划分标准对31个主要城市进行空气质量的评定，还可以在各个城市之间做比较，看出各个城市的大气环境质量的好坏程度，从而对全国各省会城市的大气环境质量做出定量的综合评价，并通过拟合预测和灰色模型预测出2012年各地区大气的污染状况。45.1.2各城市空气质量综合评价及预测根据模糊综合评价模型的求解过程我们对计算得到的数据及它们在空气质量评价中所占的地位确立以下模糊集合，即因素集权重集评价集由此建立综合评价模型R为评价因素和评价结果的模糊关系建立模糊关系矩阵为了建立从评价因素U到评价结果B的模糊关系矩阵R，首先分别在评价因素的变化域上建立对评价集B的单因素隶属函数（i=1,2,3），这里为第i种因素各级标准的平均值。大气环境质量标准及值见表1。显然，隶属函数对因素是单调下降的。将大气各级标准的数值分别代入相应的隶属函数中去，计算值见表2。表1大气环境质量标准及值级别)/(32mmgSO)/(32mmgNO)/(310mmgPM一级0.050.050.15二级0.150.100.30三级0.250.150.500.150.100.32表2大气三项污染因子的隶属函数值污染因子0.750.670.680.500.500.510.380.400.39至此，我们得出：由于各污染因子的污染程度不同，因此对大气质量的影响也不同。我们采取如下方法确定各因子的权重iiiiisxsxa//5结果如表3：表32003-2010年各城市污染物权重平均值城市北京0.560,190,25天津0.490,300,21石家庄0.550,280,17太原0.560,340,10呼和浩特0.480,280,24沈阳0.560,260,17长春0.600,160,23哈尔滨0.520,210,27上海0.460,250,29南京0.530,220,25杭州0.510,230,26合肥0.710,120,17福州0,550,120,33南昌0.500,300,20济南0.620,250,14郑州0,510,280,20武汉0,530,230,24长沙0,510,300,18广州0,410,260,32南宁0,440,320,24海口0,630,140,23重庆0,490,320,19成都0.520,250,23贵阳0,480,400,12昆明0,450,310,24拉萨0,650,070,27西安0,590,230,18兰州0,580,240,17西宁0,680,160,16银川0,560,280,16乌鲁木齐0,450,340,21隶属度函数为分段线性函数，通过隶属度函数可以计算出各评价参数对各级评价标准的隶属度。隶属函数如下：第一级：第二级：6第三级：式中ic为第i个评价参数的实测值，ijv为第i个评价参数在第j级评价标准规定的标准值（mjni,,2,1;,,2,1）。把各监测点位的实测值和各级浓度标准值代入隶属度函数公式中，即可组成单因素评价矩阵：nmnnmmmnijrrrrrrrrrrR212222111211)(以北京为例，把表3中北京的数据记为权重矩阵A=[0.560.190.25]R=[0.540.180.28]，所以，根据表2各空气质量评判标准得出北京空气质量处于二级，其余各城市求解过程以此类推，结果如表4：表42003-2010年各城市平均空气质量级别图城市空气质量级别城市空气质量级别海口一级南京二级拉萨二级合肥二级南宁二级长沙二级福州二级成都二级昆明二级重庆二级贵阳二级武汉二级呼和浩特二级沈阳二级上海二级西宁二级南昌二级石家庄二级长春二级济南二级广州二级西安二级银川二级太原二级天津二级乌鲁木齐二级郑州二级北京二级哈尔滨二级兰州三级杭州二级根据上述表格，我们可以知道全国94.5%的城市空气质量处于二级标准。海口空气质量最好处于一级标准，兰州空气质量最差处于三级标准，需进一步加强环境整治工作，其余城市都属于二级标准，但部分地区如太原、乌鲁木齐、西安空气质量处于二级标准7偏下，应注重环境保护，以免空气质量进一步恶化。5.1.3预测2012年空气质量根据《中国统计年鉴》中2SO、102PMNO、浓度数据和附件所提供的计算公式，通过代入2SO、102PMNO、浓度数据得到如下计算公式：LoSOLoSOPSOLoSOHiSOLoSOHiSOSOIAQIBPCBPBPIAQIIAQIIAQI22222222（1）LoNOLoNOPNOLoNOHiNOLoNOHiNONOIAQIBPCBPBPIAQIIAQIIAQI22222222（2）LoPMLoPMPPMLoPMHiPMLoPMHiPMPMIAQIBPCBPBPIAQIIAQIIAQI1010101010101010（3）10N22PMOSOIAQIIQAIIAQIMaxAQI，，（4）在模型类型和数据量大小的综合考虑之下，我们选择了灰色模型GM(1,1)[3]对2012年各城市的AQI值进行了预测：GM（1，1）模型其简单的微分方程形式是dxaxudt利用常数变易法解得，通解为()atuxtcea若初始条件为00,()txtx，则可得到微分方程的特解为0()()atuuxtxeaa或时间响应函数(1)(1)((1))atuuxtxe