PM2.5污染物的影響因素識別

楊帥帥 倪小雨 李放歌 王璽

摘要：有研究認為，AQI監(jiān)測指標中的二氧化硫（SO2），二氧化氮（NO2），一氧化碳（CO）是在一定環(huán)境條件下形成PM2.5前的主要氣態(tài)物體。西安市環(huán)境監(jiān)測站對這幾種污染物進行了長期監(jiān)測，并每天公布檢測數據。本文以西安市PM2.5監(jiān)測記錄為依據，運用多元回歸方法研究PM2.5污染濃度指標的影響因素，并建立回歸方程。研究表明，西安市的PM2.5污染物濃度與二氧化氮和PM10的濃度有關。

關鍵詞：空氣污染? 多元回歸? 模型修正

近年來，PM2.5污染物受到人們的普遍關注。PM，英文全稱為fine particulate matter，即細顆粒物，PM2.5是指大氣中直徑小于或等于2.5μm的顆粒物，又稱為可入肺顆粒物。粒徑在10μm以上的顆粒物，會被擋在人的鼻子外面;粒徑在2.5μm至10μm之間的顆粒物，能夠進入人體上呼吸道，但部分可通過痰液排出體外;而粒徑在2.5μm以下的亞微顆粒物，會被吸入人體肺部，進入支氣管，干擾肺部的氣體交換。雖然PM2.5只是地球大氣成分中含量很少的組分，但它對空氣質量和人類健康等有重要的影響。

對于空氣污染物的監(jiān)測，我們往往著重研究空氣中SO2、NO2、PM10、CO、O31小時、O38小時以及PM2.5的濃度。PM2.5的濃度可能受到其它污染物濃度的影響，如果確定哪些污染物濃度影響PM2.5濃度，就可以為治理PM2.5污染提供有效的突破口。

本文利用空氣污染物的監(jiān)測數據，從統(tǒng)計角度研究 PM2.5濃度的影響因素。

1 數據采集

本文數據采自西安市環(huán)境監(jiān)測站網站（http：//www.xianemc.gov.cn/）公布的空氣質量日報。西安市共有高壓開關廠等12個監(jiān)測點，每個監(jiān)測點每天對西安市空氣污染物進行檢測。本文取2012、2013年共24個月的平均值。

空氣中SO2濃度、NO2濃度、PM10濃度、CO濃度、O31小時、O38小時與PM2.5有著密切的關系，也從一定的角度反映了一個城市的空氣質量狀況。

表1.1? 各監(jiān)測點6項檢測數值24月平均（部分數據）

■

2 建立回歸模型

以PM2.5值（y）作為因變量，SO2濃度（x1：毫克/立方米）、NO2濃度（x2：毫克/立方米）、PM10濃度（x3：毫克/立方米）、CO濃度（x4：毫克/立方米）、O31小時（x5：毫克/立方米）、O38小時（x6：毫克/立方米）作為自變量，采用強制回歸的方法，進行多元回歸分析，分析結果如下：

表2.1? Model Summary

■

表2.2? ANOVAb

■

表2.3? Coefficientsa

■

回歸方程為

■=-62.855-0.453x■+1.143x■+0.820x■+0.228x■-

0.223x■■+0.461x■

由表2.1可以看出，樣本決定系數R=0.994>0.8，即相關程度很好。由表2.2看出F=744.623，Sig.=0.000，可以認為6個自變量對于因變量影響效果顯著。由表2.3看出t檢驗中，變量x1，x2，x3的Sig.值較小，分別為0.000<0.05，0.010<0.05，0.000<0.005，所以x1，x2，x3自變量對于因變量相關程度較高。

但是，通過方程可以看到SO2與O31小時的系數為負數，即當他們增多時PM2.5值減小，故該兩項數據與 PM2.5值無關，舍去這兩項數據重新做一次回歸。

以PM2.5值（y）作為因變量，NO2濃度（x1：毫克/立方米）、PM10濃度（x2：毫克/立方米）、CO濃度（x3：毫克/立方米）、O38小時（x4：毫克/立方米）作為自變量采用強制回歸的方法，再次進行多元回歸分析（計算表省略），得到如下回歸方程：

■=-61.287+0.922x1+0.838x2

樣本決定系數R=0.994>0.8，即相關程度很好。F=1007.267，Sig.=0.000，可以認為4個自變量對于因變量影響效果顯著。經t檢驗，變量x1，x2的Sig.值較小，分別為0.000<0.05，0.010<0.05，所以x1，x2自變量對于因變量相關程度較高。因此，刪去變量x3（即CO濃度），再次以x1，x2為自變量進行回歸分析，得到回歸方程如下：

■=-56.714+0.981x1+0.836x2

經逐步線性回歸看出2個自變量PM10、NO2與因變量PM2.5之間相關性高。

3 多重共線性診斷與模型修正

當自變量存在多重共線性時，利用最小二乘估計得到的回歸參數估計值很不穩(wěn)定，回歸系數的方差隨著多重共線性強度的增加而加速增長，會造成回歸方程高度顯著的情況下，有些回歸系數通不過顯著性檢驗，甚至出現回歸系數的正負號得不到合理的經濟解釋。需要對上文得到的兩個自變量NO2濃度和PM10濃度進行多重共線性診斷?？梢岳梅讲顢U大因子法以及對原始數據（逐步回歸后的數據）進行診斷。

以PM2.5值（y）作為因變量，NO2濃度（x1：毫克/立方米）、PM10濃度（x2：毫克/立方米）為自變量進行共線性診斷，分析結果如下：

綜合表3.1可以看出，PM10、NO2的方差擴大因子（VIF）比較小，且均小于10，說明建立的回歸模型原始數據中PM10、NO2不存在著明顯的多重共線性。因此，有效的回歸方程仍然為

=-56.714+0.981x1+0.836x2

4 結論

根據多種污染物濃度的數據進行回歸分析，結果表明PM2.5濃度值與PM10和NO2濃度值密切相關，并得到有效的回歸方程。這個方法可以對影響PM2.5進行初步篩選。但是不同的地區(qū)應該有不同的結果，所以這個回歸方程不具有普遍意義。

參考文獻：

[1]GB3095-2012環(huán)境空氣質量標準[S].中華人民共和國環(huán)境保護部.

[2]湯偉.霧霾治理研究與國外城市對策[J].城市管理與科技，2014.1：76-79.

[3]林清泉.計量經濟學[M].北京：中國人民大學出版社，2012.

[4]郭曉澤，單思行.針對PM2.5的綜述[J].能源與節(jié)能，2012，86（11）：58-60.

通訊作者：王璽