空气中PM25的建模研究

1、承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则 .我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子 邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关 的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其 他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式 在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违 反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从 A/B/C/D中选择一项填写）：我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的

2、全名）：参赛队员（打印并签名）：1._指导教师或指导教师组负责人（打印并签名广日期：年月日编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：空气中PM2.5问题的建模研究摘要空气质量问题始终是政府、环境保护部门和全国人民关注的热点问题，因此， 对空气质量的监测、预报和控制工作显得尤为重要。本文通过建立相关因素分析 模型、PM2.5的时空分布及污染评估模型、基于满意度的最优化模型和多目标非 线性规划模型等对空气质量监测指标一一对人类健康危害极大的细颗粒物 PM2.5的浓度进行深入研究分析。对于问题一

3、，在关于PM2.5的相关因素分析中，众多的因素对PM2.5有无统 计上的联系，直接观察无法得到结论。因此本文构建了Pearson相关系数模型，通过积差方法计算相关系数，以两变量与各自平均值的离差为基础， 通过两个离 差的乘积来反映两变量之间相关程度，并通过 SPS敞件进行相关性分析从而给 出6个基本监测指标的相关性分析。对于问题二，要求描述该地区内PM2.5的时空分布及其规律，并结合环境保 护部新修订的环境空气质量标准分区进行污染评估。用EXCEL绘制西安市13个监测点的PM2.5含量随时间变化图，并选取两组方差最大的地区绘制它们 的PM2.5含量随时间变化图。根据这两个图分析该地区内PM2.

4、5的时空分布及其 规律，并分区进行污染评估。对于问题三，空气质量的控制管理主要包括综合治理和专项治理两个方面， 实际生活中我们采用二者结合的方案逐年达到治理目标。为平衡预定PM2.5减排 计划和投入经费合理化要求，本节建立基于满意度的最优化模型和多目标非线性 规划模型，从而给出五年投入总经费和逐年经费投入的合理预算计划。关键词：Pearson相关因素模型PM2.5的时空分布及污染评估模型基于满意度的最优化模型 多目标非线性规划模型 SPSS软件一、问题重述空气质量问题始终是政府、环境保护部门和全国人民关注的热点问题，因此， 对空气质量的监测、预报和控制工作显得尤为重要，国家和地方政府均制定了相

5、 应的政策、法规和管理办法。启用空气质量指数AQI作为新的空气质量监测指标， 以代替原来的空气质量监测指标空气污染指数API （Air Pollution Index）AQI作为无量纲指数，它的分项监测指标为 6个基本监测指标（二氧化硫SQ、二氧化氮NO2、可吸入颗粒物PM10细颗粒物PM2.5臭氧O3和一氧化碳COI亨6项）。新标准中，首次将产生灰霾的主要因素一一对人类健康危害极大的细颗 粒物PM2.5的浓度指标作为空气质量监测指标2 o新监测标准的发布和实施，将 会对空气质量的监测，改善生存环境起到重要的作用。针对空气质量问题的研究， 本题主要对如下三个问题进行讨论研究：一、关于PM2.5

6、的相关因素分析：对AQI中6个基本监测指标的相关与独立 性进行定量分析，尤其是对其中 PM2.5（含量）与其它5项分指标及其对应污染 物（含量）之间的相关性及其关系进行分析。二、PM2.5的分布与演变及应急处理：描述该地区内 PM2.5的时空分布及其 规律，并结合环境保护部新修订的环境空气质量标准分区进行污染评估。三、空气质量的控制管理根据要求给出合理的PM2.5年平均治理计划，即给出每年的全年年终平均治理指标。为数据 1所在地区设计有效的专项治理计划， 使得既达到预定PM2.5减排计划，同时使经费投入较为合理，要求你给出五年投 入总经费和逐年经费投入预算计划，并论述该方案的合理性。二、问题分

7、析问题一要求对AQI中6个基本监测指标的相关与独立性进行定量分析， 尤其 是对其中PM2.5（含量）与其它5项分指标及其对应污染物（含量）之间的相关 性及其关系进行分析。由于PM2.5的影响因素太多，直接观察无法得到结果，因 此本文建立Pearson相关系数模型，通过两个离差的乘积来反映两变量之间相关程度。问题二要求描述该地区内PM2.5的时空分布及其规律，并结合环境保护部新 修订的环境空气质量标准分区进行污染评估。根据对附件2、附件3中的数据的分析，构建PM2.5的时空分布及污染评估模型，由于高压开关厂和小寨地区 的PM2.5含量数据方差最大，因此画出其 PM2.5随时间的变化规律图，以研究

8、 PM2.5的空间变化问题。用EXCE份别彳出13个监测点PM2.5含量的平均值1, 按照国家环境保护部新修订的环境空气质量标准，对13个监测点进行环境空气功能分区与评估2o问题三要求给出合理的PM2.5年平均治理计划、五年投入总经费以及逐年经 费投入预算计划，通过构建效用函数来表明规划决策的满意度，治理满意度F越高，则说明治理方案越合理。减少PM2.5的排放与投入总经费最少这两个约束 条件之间存在矛盾性，不能够同时满足，因此本题为多目标规划问题，分别给出 对各目标的分析。同时，为简化多目标非线性函数的求解， 我们采用主要目标法 将多目标非线性函数转化为以投入总经费最少为主要目标的单目标规划3

9、,建立模糊多目标非线性规划模型，从而给出五年投入总经费和逐年经费投入的合理预 算计划。三、问题假设1、假设每年各自减少的年平均浓度符合一个递减等比数列，其等比数列的 公比为q。2、假设忽略对PM2.5产生影响的其他因素，仅考虑 PM2.5的专项治理和综 合治理问题。四.符号说明符号解释说明皮尔逊相关系数XY表示第i年减少的年平均浓度U(x)为关于x的效用函数表示为PM2.5减少的年平均浓度为为的前提下综合治理的满意度投入总经费为V表示经专项治理PM2.5减少的年平均浓度注：在此仅给出主要符号说明，具体细分符号论文过程中会给出详细说明。五、模型的建立与求解PM2.5的相关因素分析模型相关性分析对

10、于众多的因素对于PM2.5有无统计上的联系，直接观察是得不到结论的， 所以我们使用了 Pearson相关系数模型，相关系数是按积差方法计算，同样以两 变量与各自平均值的 离差为基础，通过两个离差相乘来反映两变量之间相关程度。_Cov(X,Y)_XY D(X)1D(Y)(X x)(yi y) i 1 nn-2-(X x) *(yi y)i 1i 1n xyix yii 1i 1 i 1nnnnJn Xi2 ( x)2jn y2 ( yi)2, i 1i 1； i 1i 1我们使用了附件一中的数据，为了简化计算，在这里我们使用了SPS欹件进行相关性分析，得到的结论如下表 5-1所示：表5-1 描述

11、性统计量均值 标准差N二氧化硫28.1917.724238二氧化氮63.1124.475238可吸入颗粒物79.1633.129236一氧化碳25.8114.110238臭氧52.0533.795238PM2.5108.9474.696238表5-2相关性二氧化硫二氧化氮可吸入颗粒物一氧化碳PM2.5臭氧二氧化硫Pearson相关性1* .807*.678*.659* .726*-.1*79显著性（双侧）.000.000.000.000.006N238238236238238238二氧化氮Pearson相关性* .807*1*.727*.626* .734*-.063显著性（双侧）.000.0

12、00.000.000.336N238238236238238238可吸入颗粒物Pearson相关性* .678* .727*1*.586* .779*-.069显著性（双侧）.000.000.000.000.295N236236236236236236一氧化碳Pearson相关性* .659* .626*.5861* .822*-.3*81显著性（双侧）.000.000.000.000.000N238238236238238238PM2.5Pearson相关性* .726* .734*.779*.8221-.3*52显著性（双侧）.000.000.000.000.000N23823823623

13、8238238臭氧Pearson相关性-.17*9-.063-.069*-.381-.35*21显著性（双侧）.006.336.295.000.000N238238236238238238*.在.01水平（双侧）上显著相关。红色字体标注的为我们所需要的结果结果分析:由运行结果可知，PM2.5与二氧化硫、二氧化氮、可吸入颗粒物、一氧化碳、 臭氧均显著相关，其 Pearson相关系数分别为0.726、0.734、0.779、0.822、 -0.352 ,则表明PM2.5与二氧化硫、二氧化氮、可吸入颗粒物、一氧化碳呈显著 正相关；与臭氧呈负相关。PM2.5的分布与演变及应急处理总体思路：该问题要求描

14、述该地区内PM2.5的时空分布及其规律，并结合环境保护部新 修订的环境空气质量标准分区进行污染评估。用 excel绘制西安市13个监 测点的PM2.5含量随时间变化图，并选取两组方差最大的地区绘制它们的 PM2.5 含量随时间变化图。根据这两图来分析该地区内PM2.5的时空分布及其规律，并 分区进行污染评估。PM2.5的时空分布及污染评估模型：步骤一 研究PM2.5随时间变化问题：由excel绘制13个监测点的PM2.5含量随 时间变化图如图1:市人民西育新 开体翻区经区A-4- -3102 3T-4- -3TOZ 71-4- -3TOZ or-4- -3TOZ 4-3- -310267-3

15、-3TOZ 2-3- -3102 31-3- -3TOZ 7-ox- -3102 67- -o,- -3TOZ2_231-ox- -3TOZ 7-1- -3102QT-1- -3TOZ C2I -1- -3TOZ 61- -1- -Qwoz 1-1- -3102关厂兴庆小 区 纺织 城 小寨图1 13个监测点的PM2.5含量随时间变化图从下图可以看出13个监测点PM2.5数据变化趋势和幅度基本一致，而且随 着时间的推移在波动中呈逐渐下降趋势。 且可以看出冬季PM2.5的含量明显高于 春季，说明季节因素对PM2.5有很大影响。而季节因素主要影响风向变化，由此 可推断出风向也是影响PM2.5含量的

16、重要因素。步骤二研究PM2.5随空间问变化问题：由于高压开关厂和小寨地区的 PM2.5含 量数据方差最大，因此画出其PM2.5随时间的变化规律图，以研究PM2.5的空间 问题，变化图2如所示。由图可看出，小寨的PM2.5含量比高压开关厂明显滞后。查阅资料得知高压 开关厂位于西安西面，距离市中心约 5kmi小寨位于西安正南方向，距离市中心 约3kmi西安由于地理因素的影响，全年盛行西北风，将正处于小寨西北方向的 高压开关厂的污染物带至小寨，对这里的空气质量造成了较大影响。局压开关J 小寨5004003002001009-4-O-QN 4-4-O-QN 3o_31n0 4J_31n0 1-4-o-

17、QN 53-O-QN n,3302 n4oo_31n0 oz-00-3102 SXJO_31n0 nv3302 62-O-QN 7oo_31n0 2oo_31n0 RXJO_31n0 3JO_31n0 9-lo-QN 4-1-o-QN QU71302 3ol_31n0 1302 4J.1,31n0图2高压开关厂和小寨地区的 PM2.5随时间变化图步骤三用excel分别作出13个监测点PM2.5含量的平均值1 ,如下表所示:表5-3各个监测点PM2.5含量的平均值高压开关厂兴庆小区纺织城小寨市人民体育场局新西区经开区214.893805195.5333189.4636176.1363636206

18、.1414141212.3818197.89长安区阎良区临潼区曲江文化集团广运潭草滩全市平均179.872727191.8455180.2449183.8990826205.4205607211.5405195.0522由上表可得：高压开关厂的PM2.5含量是13个监测站中最大的；市人民体 育场、 广运潭、草滩PM2.5含量也较高；小寨、临潼区的 PM2.5含量较低。 步骤四按照国家环境保护部新修订的环境空气质量标准，对13个监测点进 行环境空气功能分区与评估20则一级区为；曲江文化集团、广运潭、草滩；二 级区为：高压开关厂、兴庆小区、纺织城、小寨、高新西区、市人民体育场、经 开区、长安区、阎

19、良区、临潼区。西安市大气环境中的可吸入颗粒物含量高，13个监测点pm2.5超标现象均不低于70%出现在高压开关厂超标率甚至高达88.6%, 市人民广场超标率达84%可见工业污染、商业区人流、车流多，绿化率低导致 了这一现象的出现。高新技术开发区、旅游度假区的绿化率高，工业企业的排污 量少，可吸入颗粒物的年日均值较少。上述分析可见，研究区域范围内的大气环境质量最差的区域出现在工业区和 商业区，高新技术开发区和旅游度假区的空气质量较好，但全市PM2.5质量普遍不好。5.3空气质量的控制管理1总体思路空气质量的控制管理主要包括综合治理和专项治理两个方面， 实际生活中我 们采用二者结合的方案逐年达到治

20、理目标。为平衡预定PM2.5减排计划和投入经 费合理化要求，本节建立基于满意度的最优化模型和多目标非线性规划模型， 从 而给出五年投入总经费和逐年经费投入的合理预算计划。基于满意度的最优化模型步骤一根据未来五年内逐年减少PM2.5的年平均浓度这一问题要求，本节 假设每年各自减少的年平均浓度符合一个递减等比数列，其等比数列的公比为q本问题属于规划与决策分析问题，可以通过构建效用函数3来表明规划决策的满 意度。PM2.5减少的年平均浓度越高，决策部门越满意，但是同时当年投入的治 理总经费就会相对较高，因此可取效用函数：xU(x) 一(3) xi aX表示第i年减少的年平均浓度，系数a通过“高频滤波

21、系统中下限截止频率”的计算方法确定为0.2 o步骤二 假设治理?f意度为F,满意度越高，则说明治理方案越合理，则未 来五年的最大满意度表示为5U(xi)MaxF (4)5其中，X表示第i年减少的年平均浓度，U(x)为关于x的效用函数，F表F的目标规划示为PM2.5减少的年平均浓度为x的前提下，综合治理的满意度。步骤三 根据上述对最大满意度的目标分析，构建基于满意度函数为:5MaxU(x)i 15st.U(Xi)5Xi i 1XiXiX 0.2,280 35,qxi 1,q 0,(6)99.59%,逐年Xi 0,i 1,2,3,4,5运用Lingo编程（见附件2）可求最优结果为：最大满意度为 减

22、少上年减少量的95.31%,未来五年的年平均浓度治理计划如表所示:表5-4未来五年PM2.5的具体年平均浓度治理计划表年份该年PM2.5年平均减少的浓度（单位为 m/m ）该年年终PM2.5要求达到的年平均浓度（单位为 m/m3）第一年53.8165226.1835第二年51.29264174.89086第三年48.88714126.00375第四年46.5944579.40928第五年44.40928355.3.3模糊多目标非线性规划模型经上述分析可知，减少PM2.5的排放与投入总经费最少这两个约束条件之间 存在矛盾性，不能够同时满足，因此本小题为多目标非线性规划问题4,下面是 对各目标的分

23、析: 目标一投入总经费最少根据题目要求，本题忽略对PM2.5产生影响的其他因素，仅考虑PM2.5的专 项治理和综合治理问题。假设投入总经费为Z ,专项治理和综合治理投入的费用 分别为乙,Z2，则投入的总费用最少可表示为 TOC o 1-5 h z Min Z Z1 Z255其中，ZiXi , Z20.005(Xi yi), Xi表示第i年经综合治理PM2.51i 1减少的年平均浓度，y表示经专项治理PM2.5减少的年平均浓度。目标二PM2.5减排计划满意度最大化经分析可知，影响PM2.5减排计划的主要因素为专项治理和综合治理的投入 费用，构建关于专项治理和综合治理的总满意度函数F,总的满意度越

24、高，说明PM2.5的减排计划越合理，则未来五年内PM2.5减排计划的实施满意度最大可表 示为 TOC o 1-5 h z / 55Max F - U(x)/5(U(xi) U(yi)/5(8)i 1i 1其中，F表示综合治理和专项治理下的满意度，U(Xi)表示综合治理的效用5函数，U(y)表示专项治理下的效用函数，U (为)/5表示未来五年综合治理下i 15的满意度，U(yi)/5表示未来五年专项治理下的满意度。i 1模糊多目标非线性规划模型的建立基于上述的多目标分析，结合相关约束条件，构造多目标非线性规划函数：Min Z Z1 Z255Max F - U(Xi)/5(U(Xi) U(yi)/

25、5i 1i 1st. TOC o 1-5 h z ZiU(x)i 15_2Z20.005(x yji 155Xiyi 280 35Xiyi 1)2i 1 i 120.005(xiyi)qx 1 0.005(xi 1U(Xi)U(yJqX 0,(iy 0,(ixxi 0.2yi0.201,234,5)1,234,5)5.3.4模型的求解为简化多目标非线性函数的求解，采用主要目标法将多目标非线性函数转化为以投入总经费最少为主要目标的单目标规划5,简化后的其他约束不变，仅考虑以下变化、Min Z Z1 Z2(10)1 55F -U(x)/5 (U (xi) U(yJ)/5b2 i 1i 1b的取值影

26、响着满意度 F的大小，本文经综合考虑，取 b=90%运用Lingo 编程（见附件2）可求最优结果为：在保证满意度的前提下，投入治理的总经费 最少为64.03053百万单位，q为1,未来五年具体的投入总经费和逐年经费投入 预算计划表如下：表5-5未来五年综合治理和专项治理投入经费的年预算表年份该年综合治理投入经费（白力单位）该年专项治理投入经费（白力单位）第一年0.821112.005第二年0.821112.005第三年0.821112.005第四年0.821112.005第五年0.821112.005合计4.105560.025六、模型评价基于Pearson相关系数的统计推断对数据的分布类型是很敏感的，不论变量 之间的联合概率密度函数是不是正态的,Pearson相关系数都是用来衡量两个随 机变量之间的线性相关程度的。对于二元正态数据，样本的相关系数是总体相关 系数的