杜宇靜+李佩媛

摘 要 本文主要以2016年吉林市冬季和春季每日的PM2.5值與平均溫度、風(fēng)速和相對(duì)濕度的數(shù)據(jù)分別建立回歸模型，并對(duì)模型進(jìn)行多重共線性檢驗(yàn)、自相關(guān)性檢驗(yàn)、異方差性檢驗(yàn)等，得到的結(jié)果為風(fēng)速與PM2.5負(fù)相關(guān)，相對(duì)濕度與PM2.5正相關(guān)，并對(duì)冬季和春季兩種情況作對(duì)比分析，最后根據(jù)得到的模型做預(yù)測(cè)分析，預(yù)測(cè)的結(jié)果較好。

關(guān)鍵詞 PM2.5 回歸分析 參數(shù)估計(jì) 顯著性檢驗(yàn) 預(yù)測(cè)

中圖分類號(hào)：X513 文獻(xiàn)標(biāo)識(shí)碼：A

0引言

PM2.5又稱細(xì)顆粒物、細(xì)粒，對(duì)空氣質(zhì)量和能見(jiàn)度有非常大的影響，且能長(zhǎng)時(shí)間懸浮于空氣中。細(xì)顆粒物以自身直徑小的優(yōu)勢(shì)，能進(jìn)入呼吸道更深的部位，對(duì)人體健康的危害更大。很多人以為，PM2.5只是空氣污染，其實(shí)它對(duì)整體的氣候條件影響更糟糕，這就是我們看到藍(lán)天白云的機(jī)會(huì)越來(lái)越少的原因。各種現(xiàn)象和研究說(shuō)明，PM2.5同時(shí)影響著空氣質(zhì)量和天氣氣候；反而言之，在范圍較大的氣候條件下，各氣象因素也是可以影響空氣中PM2.5的濃度的。近幾年，人們已經(jīng)對(duì)PM2.5的關(guān)注增大，對(duì)其環(huán)境污染方面的研究也逐漸增多。周兆媛研究京、津、冀三市空氣污染與氣象要素的相關(guān)關(guān)系。王琪等研究了北京地區(qū)PM2.5質(zhì)量濃度分布及其氣象條件影響關(guān)系，鄭煜，藺澤一等以哈爾濱冬季為對(duì)象進(jìn)行相關(guān)性分析。

吉林市在我國(guó)東北地區(qū)，屬于溫帶大陸性季風(fēng)氣候，四季分明。由于近年來(lái)全球氣候變化加劇，吉林市的經(jīng)濟(jì)水平也在不斷提高，吉林市的空氣質(zhì)量波動(dòng)也在逐年增強(qiáng)，吉林市的污染情況在冬季尤為嚴(yán)重。本文分別研究了吉林市2016年冬季（2016.11.07-2017.02.02）和2017年春季（2017.02.03-2017.05.04）的空氣污染中影響PM2.5濃度與氣候因素的相關(guān)關(guān)系。

1模型的建立

本文選取吉林市2016年11月7日到2017年5月4日共計(jì)179天的數(shù)據(jù)作為原始數(shù)據(jù)。設(shè)因變量為PM2.5值（冬季：y1，春季：y2），自變量平均溫度（冬季：x11，春季：x12）， 風(fēng)速（冬季：x21，春季：x22）， 相對(duì)濕度（冬季：x31，春季：x32）。

1.1模型的建立

首先繪制散點(diǎn)圖，然后判斷因變量與自變量是否存在線性關(guān)系。下面圖1和圖2分別給出y1與x11、y2與x32散點(diǎn)圖：

由散點(diǎn)圖可以看出，平均溫度和PM2.5值存在線性相關(guān)關(guān)系，同樣風(fēng)速和相對(duì)濕度與PM2.5值也存在線性相關(guān)關(guān)系。

多重共線性判定主要是檢驗(yàn)自變量之間是否近似線性相關(guān)，這里用方差擴(kuò)大因子法，計(jì)算共線性統(tǒng)計(jì)量VIF。冬季模型（1）中VIF值分別為1.021，1.035，1.055， 春季模型（2）中VIF值分別為1.938，1.266，1.8.13，VIF值均小于10，說(shuō)明冬季模型和春季模型中的平均溫度、風(fēng)速和相對(duì)濕度之間不存在多重共線性。

建立模型時(shí)，某些因素會(huì)導(dǎo)致隨機(jī)誤差項(xiàng)產(chǎn)生不同方差，稱為異方差性。利用等級(jí)相關(guān)系數(shù)檢驗(yàn)法來(lái)判斷回歸模型的異方差性，此假設(shè)檢驗(yàn)問(wèn)題睥原假設(shè)為隨機(jī)誤差方差相等，通過(guò)計(jì)算得到冬季模型（1）中殘差的絕對(duì)值與三個(gè)自變量之間的等級(jí)相關(guān)系數(shù)分別為0.114，0.072，-0.048，對(duì)應(yīng)的P值分別為0.091，0.058，0.046。春季模型（2）中殘差的絕對(duì)值與三個(gè)自變量之間的等級(jí)相關(guān)系數(shù)分別為0.000，-0.030，-0.005，對(duì)應(yīng)的P值分別為0.099，0.051，0.035。因此可以判斷不存在異方差性。

若模型中隨機(jī)誤差項(xiàng)之間不滿足cov（ i， j）=0，i≠j這個(gè)條件，則稱隨機(jī)誤差項(xiàng)之間存在自相關(guān)現(xiàn)象。

經(jīng)研究得出，隨機(jī)誤差項(xiàng)目有自相關(guān)性。另外，我們用DW檢驗(yàn)法，對(duì)冬季回歸模型（1）觀測(cè)值的數(shù)目n=88，解釋變量的數(shù)目k=4，查閱DW檢驗(yàn)上下界表可知：dL=1.59，dU=1.73。計(jì)算得出DW值為1.019，滿足0≤DW≤dL，可知誤差項(xiàng)存在自相關(guān)。自相關(guān)系數(shù) 1=10.5DW=10.51.019=0.4905。

同理，可以判斷春季模型（2）中的隨機(jī)誤差項(xiàng)之間存在自相關(guān)。我們可以用迭代法消除自相關(guān)性：

并對(duì)模型進(jìn)行顯著性檢驗(yàn)，F(xiàn)檢驗(yàn)法得到的結(jié)果如下表1：

由表1可知，P值均比0.05大，說(shuō)明模型的結(jié)果并不好。用t檢驗(yàn)法對(duì)回歸系數(shù)進(jìn)行檢驗(yàn)，即檢驗(yàn)回歸模型中每個(gè)變量的系數(shù)是否為零，檢驗(yàn)結(jié)果為對(duì)于冬季模型（3）中的變量x11，x21，x31 t統(tǒng)計(jì)量值分別為1.748，1.398，2.013，對(duì)應(yīng)的P值分別為0.084，0.066，0.047，對(duì)于春季模型（4）的變量x12，x22，x32 t統(tǒng)計(jì)量值分別為1.579，1.301，1.765，對(duì)應(yīng)的P值分別為0.118，0.197，0.081，變量x11、x12的P值相對(duì)較大，所以需要對(duì)模型進(jìn)行改進(jìn)。

將變量x11、x12從模型（5）（6）中去掉，分別得到改進(jìn)的冬季的回歸模型：

春季的回歸模型：

經(jīng)檢驗(yàn)F值分別為3.772，3.014對(duì)應(yīng)的P值分別為0.027，0.041，回歸方程顯著。

3預(yù)測(cè)和結(jié)論

預(yù)測(cè)是指在掌握現(xiàn)有信息基礎(chǔ)上，依照一定的方法和規(guī)律對(duì)未來(lái)的事情進(jìn)行預(yù)測(cè)，以預(yù)先了解事情發(fā)展的過(guò)程與結(jié)果。下面根據(jù)回歸方程進(jìn)行預(yù)測(cè)。經(jīng)過(guò)計(jì)算2017年5月5日的預(yù)測(cè)值為52.1，而真實(shí)值為54，偏差為1.9，預(yù)測(cè)的結(jié)果相對(duì)較好，由預(yù)測(cè)結(jié)果可以看出建立的模型擬合度好。

從結(jié)果可以看出，吉林市冬季和春季PM2.5質(zhì)量濃度有明顯的季節(jié)區(qū)別，冬季較高，春季較低，這和吉林市冬季大范圍取暖有必要的關(guān)系。此外，這兩個(gè)季節(jié)PM2.5質(zhì)量濃度的變化也受到了各季節(jié)不同的天氣氣候的影響。

吉林市冬季和春季PM2.5質(zhì)量濃度的日變化與風(fēng)速負(fù)相關(guān)，和相對(duì)濕度有正相關(guān)關(guān)系，這與鄭煜等對(duì)哈爾濱冬季氣象要素對(duì)空氣質(zhì)量指數(shù)的影響分析的結(jié)果相同。

兩個(gè)季節(jié)中，平均溫度、風(fēng)速和相對(duì)濕度對(duì)PM2.5質(zhì)量濃度的影響作用有所不同，風(fēng)速的影響作用最大，而相對(duì)濕度的影響作用最小。綜合作用下的效果：冬季模型和春季模型中，如果風(fēng)速變大、相對(duì)濕度變小可以使PM2.5的質(zhì)量濃度降低。

如果想通過(guò)改變一些氣象因素的影響來(lái)降低PM2.5的質(zhì)量濃度，可以采取以下個(gè)方法：減少降水量降低空氣的相對(duì)濕度，相對(duì)濕度的降低可以使PM2.5質(zhì)量濃度慢慢降低；可以采取一些倚靠大自然的方法來(lái)增加城市的空氣流通情況，即風(fēng)速對(duì)PM2.5質(zhì)量濃度的影響，空氣流通加快，城市的PM2.5質(zhì)量濃度會(huì)降低；當(dāng)然，也可以通過(guò)多個(gè)氣象因素的互相作用來(lái)改變PM2.5的質(zhì)量濃度。

參考文獻(xiàn)

[1] 周兆媛，張時(shí)煌，高慶先，李文杰，趙凌美，馮永恒，徐明潔，施蕾蕾.京津冀地區(qū)氣象要素對(duì)空氣質(zhì)量的影響及未來(lái)變化趨勢(shì)分析[J].資源科學(xué)，2014（01）.

[2] 王琪，孫巍，張新宇.北京地區(qū)PM_（2.5）質(zhì)量濃度分布及其與氣象條件影響關(guān)系分析[J].計(jì)算機(jī)與應(yīng)用化學(xué)，2014（10）：1193-1196.

[3] 鄭煜，藺澤一，張燕杰，張國(guó)財(cái).哈爾濱冬季氣象要素對(duì)空氣質(zhì)量指數(shù)的影響分析[J].黑龍江大學(xué)自然科學(xué)學(xué)報(bào)，2015（06）：805-809.

[4] 何曉群，劉文卿.應(yīng)用回歸分析[M].中國(guó)人民大學(xué)出版社，2011.endprint