关于北京市气候变化与大气污染的研究
中央民族大学 马越龙、姚琨、朱爱丽
摘 要
本文建立了一个关于分析北京市气候变化趋势和大气污染情况的统计模型。本文主要研究以下三个问题：北京市气候变化的趋势，验证北京市是否存在气候变暖现象；大气污染情况及变化趋势；大气污染和气候变化之间的关系。
首先就问题一，根据主成分分析思想，我们选用降水量、气温和风速来表征北京市气候变化，并从《北京市统计年鉴2010》获得了1978——2009年的相关数据，利用Matlab、
Excel对这些数据进行时间序列分析并根据最小二乘法原理将散点图拟合成趋势曲线，又
根据问题实际需要和判定系数等对其进行筛选，得到我们需要的最优回归方程及曲线。据此验证北京市存在气候变暖现象，并选取气温的气候倾向率（气温线性回归方程的一次项系数）来度量气候变暖的程度。
其次就问题二，根据主成分分析思想，我们选用SO2浓度、NO2浓度、可吸入颗粒物含量来表征北京市大气污染状况，并从《北京市统计年鉴2010》获得了2000——2009年的相关数据。利用内地空气污染指数的计算方法，计算出大气污染的评价度量指标API。并利用问题一中的方法分析各数据，得到SO2浓度、NO2浓度、可吸入颗粒物含量及API的回归模型，从而分析度量大气污染情况的变化趋势。
然后就问题三，利用SPSS软件计算分析API与温度、降水量、风速的相关性。通过分析比较相关系数得出大气污染与气候变化之间的关系。
最后我们利用F检验来判断以上所建立回归模型的回归效果。
综上，我们利用SPSS、Excel 、Matlab等数学软件很好地实现了数据统计分析，并结合主成分分析及时间序列分析思想给出了北京市气候变化和大气污染变化趋势的分析。将数据与图形相结合，既充分发挥了图形的直观功能，又根据数据进行了相关分析，排除了主观经验的干扰，使所建立的数学模型能较好的解决上述问题。
关键词：气候变化 大气污染 变化趋势 API 回归模型 F检验
一、 问题的提出与分析
1、问题的提出
目前，气候变暖和大气污染成为当今人类社会亟待解决的两大问题，是人类必须面对的严峻挑战。北京作为中国的首都，因其特殊的政治、经济地位，在国家乃至世界的发展中都占有特殊地位，因此研究北京的气候变化和大气污染及其之间的关系有着特别重要的意义。
气候是地球上某一地区多年时段大气的一般状态，是该时段各种天气过程的综合表现。气象要素(温度、降水、风速等)的各种统计量(均值、极值 、概率等)是表征气候的
基本依据。大气污染是指自然或人为原因使大气中某些成分超过正常含量或排入有毒有害的物质，对人类、生物和物体造成危害的现象[1]。按中国《大气环境质量标准》规定的常规分析指标有总悬浮微粒、二氧化硫、氮氧化物、一氧化碳和光化学氧化剂。
现在我们所关心的问题就是：
（1）对获得的大量数据怎样进行合理有效的筛选及分析，进而建立适当的数学模型，来分析北京市气候变化的趋势，验证北京市是否存在气候变暖现象并选择一个指标来度量气候变暖的程度。
（2）导致北京市大气污染的因素众多，如何选取适当的指标评价和度量大气污染情况及变化趋势。
（3）如何依据现有数据分析大气污染和气候变化之间的关系。 （4）如何选取适当的方法在模型建立并求解之后进行模型的检验。 2、问题的分析
解决上述问题需要大量的数据，结合主成分分析思想筛选数据、掌握数据所代表的实际意义及变化趋势和数据间的联系至关重要。为此我们需要结合时间序列分析思想对表征北京市气候和大气污染的要素做个别分析及联合分析。
根据问题的要求，我们分三个问题进行分析。
问题一：对气候变化趋势的分析，可以先搜集历年的气候数据，选取适当的指标作散点图，采用一元回归的的方法建立模型，观察其波动性及变化趋势。
问题二：对大气污染的度量和变化趋势的分析，可以先搜集历年相关污染物数据，并计算出适当的指标来度量大气污染情况，然后同问题一的方法进行变化情况分析。
问题三：可以研究度量大气污染情况的指标与表征气候变化的要素之间关系，并用其来代表大气污染与气候变化之间的关系[2]。
二、基本假设
1.温度、降水量、风力可以完全表征北京市气候
2.假设北京市每个区的气候差异不大，北京市气象站的所测得的数据可以代表北京市整体的气候状况
3.在回归模型中，拟合的二次多项式可以较好模拟出数据的变化趋势
4.在变暖程度分析中，拟合的一次多项式系数可以作为气温的气候倾向率来进行变暖程度的度量
5.可吸入颗粒物、SO2和NO2可完全代表北京市大气污染的污染物 6.API指标能正确评价和度量大气污染情况
7.假设气温和API都是正态分布，或接近正态的单峰分布，从而可以用皮尔逊相关系数表示二变量间的相关关系
三、符号说明
1.API:空气污染指数，具体分级标准见附录五 2.I：某污染物的实际污染指数 3.C：该污染物在大气中的实际的浓度
4.C大与C小：在API分级限值表（附录五）中最贴近C值的两个值，C大为大 于C的限值，C小为小于C的限值
5.I大与I小：在API分级限值表（附录三）中最贴近I值的两个值，I大为大于 I的值，I小为小于I的值
6.API＝max(I1,I2，…，In)：取所有污染物污染指数的最大值作为API的实际值，这里只取n=3
7.R2：回归模型的判定系数 8.x:气候特征值 9.a:气候倾向率
四、模型的建立与求解
4.1利用实际数据分析北京市气候变化趋势，验证北京市存在气候变暖现象，并给出变暖程度的度量
4.1.1数据总体分析，并画出散点图
经查阅《北京市统计年鉴2010》，我们得到了1978——2009年的气候数据（见附录一），由模型假定，我们选取温度、降水量、风速表征气候情况，并给出了对应的描述统计量（见表1）及散点图（见附录二）。
1978年-2009年气候数据描述统计量
统计量 降水量（mm） 平均气温 最高气温 最低气温 平均风速（m/s） 平均 标准误差 中位数 众数 标准差 方差 峰度 偏度 区域 最小值 最大值 求和 观测数 最大(1) 最小(1) 置信度(95.0%) 537.615625 12.753125 37.4625 -13.44375 2.334375 0.051118414 2.4 2.4 0.289169417 0.083618952 10.26617233 -2.748726014 1.5 1.1 2.6 74.7 32 2.6 1.1 26.17769 0.12668099 0.2981715 0.3489598 498.3 12.9 12.8 37.3 37.5 -14 -15.4 148.0833769 0.71661587 1.6867128 1.9740147 21928.68652 0.51353831 2.845 3.8967339 -1.1541088 0.57331017 0.4220035 0.0681798 0.089242538 -0.935646 0.7219566 0.7196512 546.3 266.9 813.2 17203.7 32 813.2 266.9 3 11 14 408.1 32 14 11 6.9 35 41.9 1198.8 32 41.9 35 8.3 -17 -8.7 -430.2 32 -8.7 -17 53.38975055 0.25836757 0.6081248 0.7117082 表1
0.104256692 4.1.2气候变化回归模型的建立及气候变化趋势的分析
为分析北京市气候变化趋势，我们根据最小二乘法原理，使用Excel作了气候各表征量的各种回归分析并作出了相应的趋势变化图，包括线性回归、对数回归、多项式回归、乘幂回归、指数回归。根据判定系数R2的大小，我们发现多项式的拟合效果最佳且与阶次呈正相关。又由于实际分析趋势的需求，我们不需要阶次>=3的曲线拟合，最终选取了二阶多项式非线性回归作为数学模型来近似表征气候变化趋势[3]。相应变化趋势图及二阶多项式、R2见下图。
图表一900800700600500400300y = -0.3465x2 + 1376x - 1E+06200R2 = 0.15411000197519801985199019952000200520102015降水量多项式 (降水量)
图1