1 天津市空气质量评价与预测
摘要
本文对天津市区的空气质量进行了评价，并选出了主要的污染物进行研究分析，运用综合指数评价法和回归分析等方法对其空气质量进行分析，综合各种因素我们建立了如下模型。
1、本文对2001-2010年的空气污染指数和空气质量状况进行了分析，采用了指标评价法和综合指数评价法分析了对空气污染最主要的物质， 对每每年每种物质用(uuiuouoiICCCCIII)( )这个式子计算它们的污染指数，那么计算得到的最大的值的那种物质即是天津的主要污染物，我们发现对天津市空气质影响最大的物质是10PM；
2、运用spss软件我们对天津过去十年的主要污染物进行多种模型拟合分析，发现三次曲线模型的相关系数最接近1，曲线模型拟合度较高，因此我们运用三次曲线模型进行预测分析，运用下面式子Cxbxbxby12233我们得到了未来五年主要污染物浓度预测值，再根据这些值来分析与评价未来五年的空气质量；
3、研究发现我国各个城市的空气质量有类似的，也有差别很大的，文中我们利用了聚类分析法来分析我国各个城市的空气质量。
[关键词] 综合指数评价 回归分析 主要污染物 三次曲线模型 聚类分析
2
一、问题提出
空气是地球上的生物赖以生存的物质，是必不可少的一种物质。随着人类文明和经济的发展，空气污染越来越严重，尤其是工业城市，如何改善空气质量、合理进行大气环境质量预测预警、寻求有效的控制措施是当前环境科学研究的重要内容。
空气质量的好坏直接反映了空气的污染程度，它是依据空气中污染物浓度的高低来判断的，所以控制污染物的排放是改善空气质量的根本措施。空气污染的污染物主要有二氧化硫（SO2）、二氧化氮（NO2）、可吸入悬浮颗粒物\浮尘（PM10）等等。
目前，城市空气质量污染指数的分级标准是根据空气污染指数（API）的取值界定的，空气污染指数指常规监测的几种空气污染物浓度简化成为单一的概念性指数值形式，并分级表征空气污染程度和空气质量状况。
天津是中国重要的能源与工业城市，其工业化与城市化的快速发展对城市环境产生了重要影响。近年来，市委市政府立足资源型城市可持续发展战略，努力改善全市空气质量。
（1）根据你们所掌握的数据资料、选出影响空气质量的主要污染物、建立适当的数学模型对天津市空气质量进行综合评价；
（2）对天津市未来几年空气质量及未来某段时间的主要污染物浓度进行预测；
（3）定量分析全国主要城市空气质量并进行分类；
（4）根据所建立的数学模型，分析主要污染物的来源及控制措施，向有关当局写份建议性报告。
二、问题分析
空气是地球上的生物赖以生存的物质，空气中含有很多物质，这些物质对生物本身是有益而无害的，但是随着我国第二产业的迅速发展，我国一些工业城市的污染越来越严重，例如天津就是以工业为主的城市，本文也对天津的空气污染指数进行了分析，我们知道物质过量的出现在空气中会给人类甚至各种生物带来一定的危害，我们需要采取一定的措施来改善空气环境。
问题一：根据《中国人民共和国环境保护法》和《中华人民共和国大气污染防治法》规定，空气质量的好坏反映了空气污染程度，一般是依据空气中污染物浓度的高低来判断。以天津的实际具体情况分析，被计入控制污染指数(API)的污染物项目为：SO2、NO2、PM10。因此我们通过综合指数评价法和时间序列法等来对天津这个工业城市进行空气指数的模型建立与求解。
问题二：预测天津未来几年的空气质量首先要对主要污染物的浓度进行分析，再根据计算得到的未来几年污染物的浓度来分析天津未来几年的空气质量，针对这两个方面我们采用了回归分析法来计算天津未来五年的空气质量；
问题三：我国空气质量情况在各个城市都有一定的区别，那么我们可以将情况类似的城市归为一类，这样也比较方便研究问题，所以我们聚类分析法进行研究计算。
3 三、模型假设
3.1 忽略空气中在指标控制范围内的其他污染物；
3.2 假设所给数据均有效；
3.3 假设在较近一段时间内，不发生重大工业事故；
3.4 假设在未来一段时间内，城市自然环境稳定，不发生一些较大的自然灾害，例如：地震、洪灾、海啸等；
3.5 天津未来一段时间内，政府没有出台关于大规模工业的迁入迁出政策；
3.6 假设在较近的一段时间内，国家政治、自然环境稳定。
四、符号说明与名词解释
4.1 符号说明
符号 符号意义
Ii 第i种污染物的污染指数
Ci 第i种污染物的浓度(mg/m3)
Co 标准限度的上限
Cu 标准限度的下限
Io 标准污染指数的上限(mg/m3)
Iu 标准污染指数的下限(mg/m3)
A 常数
4.2 名词解释
1、空气污染指数(Air pollution index ,简称API)[3]，就是将常规监测的几种空气污染物浓度简化成为单一的概念性指数值形式，并分级表征孔子去染程度和空气质量状况，适合于表示城市的短期空气质量状况和变化趋势。中国计入空气污染指数的项目暂定为：二氧化硫、氮氧化物、总悬浮颗粒物。
2、拟合度：拟合度检验是对已制作好的预测模型进行检验，比较它们的预测结果与实际发生情况的吻合程度。
4 五、模型的建立与求解
5.1 综合评价
首先，我们对数据进行预处理，根据数据的顺序分别对每一年的各种污染物进行污染指数计算：
uuiuouoiICCCCIII)( （5-1-1）
且同年的空气污染指数API=max（I1，I2，I3）。 （5-1-2）
表5.1.1.空气污染指数对应的污染物浓度限值
污染指数 污染物浓度（mg/m3）
API SO2 NO2 PM10
（均值） （均值） （均值）
50 0.050 0.080 0.050
100 0.150 0.120 0.150
200 0.800 0.280 0.350
300 1.600 0.565 0.420
400 2.100 0.750 0.500
500 2.620 0.940 0.600
据表5.1.1的标准值，找出每一年的每一种污染物介于的标准限度的上限Co与下限Cu，及标准污染指数的上、下限（Io、Iu）；再将数值带入（5-1-1）和（5-1-2）中进行计算；
计算可得到如表5.1.2结果：
表5.1.2.各种污染物的污染指数
年份 SO2 NO2 PM10 API
2001年 52 30 108.5 108.5
2002年 59.5 28.75 94 94
2003年 62 32.5 91.5 91.5
2004年 61.5 32.5 80.5 80.5
2005年 63 29.375 78 78
2006年 58.5 30 82 82
2007年 56 26.875 72 72
2008年 55.5 25.625 69 69
2009年 52.5 25 75 75
2010年 52 28.125 73 73
5 从表5.1.2中可分析得出，从2001年到2010年，污染物指数最大的都为PM10，及可吸入悬浮颗粒物\浮尘，则根据空气污染指数的定义可确认PM10是最主要的污染物。
为了更好的评价天津市空气质量，我们查得空气污染指数范围及相应的空气质量类别[3]，另外运用了EXCEL对2001-2010年的API画图进行了分析，分别为表5.1.3和图5.2.11.4。
表5.1.3. 空气污染指数范围及相应的空气质量类别
空气污
染指数
API 空气质
量状况 对健康的影响 建议采取的措施
0～50 优 可正常活动 51～100 良
101～150 轻微污染 易感人群症状有轻加剧，健康人群出现刺激症状 心脏病和呼吸系统疾病患者应减少体力消耗和户外活动 151～200 轻度污染
201～250 中度污染 心脏病和肺病患者症状显著加剧，运动耐受力降低，健康人群中普遍出现症状 老年人和心脏病、肺病患者应在停留在室内，并减少体力活动 251～300 中度重污染
>300 重污染 健康人运动耐受力降低，有明显强烈症状，提前出现某些疾病 老年人和病人应当留在室内，避免体力消耗，一般人群应避免户外活动
图5.2.11.4.
观察表5.1.3和图5.2.11.4不难发现在过去的十年间天津空气质量徘徊于50～100之间，对应标准表发现，天津空气质量良好，且呈缓慢下降的趋势，说明在过去十年间对天津环境的治理得到了一定的成果，比较适合人们生活。在国家环境保护部的数据中心[4] 我们可查到，在过去的几年间，天津空气质量好于二级的天数基本都占有效监测天数的百分比在持续上升，说明空气质量正在朝好的方向不停发展。符合我们所建立模型得出的结论，因此可证明该模型在评价空气质量时是可靠的。 2001~2010年API指数
108.5 94 91.5 78 82 72 69 75 73 80.5
0 50 100 150
2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 年份
6 5.2 预测未来年空气质量与主要污染物
首先根据问题一所得到的结果，我们发现天津的主要污染物是PM10，因此我们先建立对于PM10未来几年进行预测的数学模型，对数据进行预处理，画出散点图，通过观测图图形我们可以发现PM10的发展趋势大约为回归模型，运用spss软件对过去十年的PM10浓度数据进行多种模型拟合分析，得出图2.和表4.分别如下：
图5.2.1
表5.2.2
模型汇总
因变量: PM10
方程 模型汇总
R2 F df1 df2 Sig.
线性 .771 26.886 1 8 .001
对数 .922 94.203 1 8 .000
倒数 .882 60.009 1 8 .000
二次 .931 46.976 2 7 .000
三次 .935 28.725 3 6 .001
复合 .794 30.775 1 8 .001
幂 .905 76.300 1 8 .000
S .821 36.688 1 8 .000
增长 .794 30.775 1 8 .001
指数 .794 30.775 1 8 .001