山東省PM2.5濃度時空變化規(guī)律及其影響因素分析

范開放，單寶艷

（山東建筑大學測繪地理信息學院，山東濟南250101）

山東省PM2.5濃度時空變化規(guī)律及其影響因素分析

范開放，單寶艷*

（山東建筑大學測繪地理信息學院，山東濟南250101）

細顆粒物（PM2.5）是空氣質(zhì)量監(jiān)測的重要對象，研究PM2.5的時空變化規(guī)律對于治理霧霾污染具有重要意義。文章基于山東省17個城市在2014～2015年共728 d的空氣質(zhì)量監(jiān)測數(shù)據(jù)和相關(guān)氣象數(shù)據(jù)，采用時間序列分析、相關(guān)分析、空間自相關(guān)分析和空間插值等方法研究了時空變化規(guī)律和氣象因素對山東省PM2.5濃度的影響。結(jié)果表明：山東省PM2.5濃度呈U型變化，變化的時間序列自相性較強，一階自相關(guān)系數(shù)達0.6896，PM2.5濃度空間分布規(guī)律明顯；風向?qū)ι綎|省不同區(qū)域PM2.5濃度的聚集和擴散具有一定作用，與氣象條件均呈現(xiàn)一定相關(guān)性，不同月份其相關(guān)性存在明顯差異；PM2.5還與SO2、CO、NO2、PM10等污染物存在極強的正相關(guān)，與O3存在負相關(guān)。

PM2.5；時空變化；空間自相關(guān)；空間插值

Key words：PM2.5；spatial-temporal characteristics；spatial autocorrelation；spatial interpolation

0 引言

近年來，我國很多地區(qū)遭受著霧霾天氣的困擾。細顆粒物PM2.5是我國城市大氣的首要污染物，是空氣質(zhì)量狀況的重要監(jiān)測對象。PM2.5不僅會導致大氣能見度下降，還會增加呼吸道系統(tǒng)疾病發(fā)病率和死亡率，已引起人們的廣泛關(guān)注［1-2］。對PM2.5研究涉及環(huán)境科學、醫(yī)學、大氣科學、地理科學等學科，其中環(huán)境科學側(cè)重于對PM2.5的組成成分、質(zhì)量濃度、物理化學特性、來源、監(jiān)測技術(shù)等方面研究［3-4］；醫(yī)學則側(cè)重于分析PM2.5對人類健康的危害［5-6］；大氣科學注重分析PM2.5的成分特征和來源，并從氣象學的角度研究PM2.5濃度與氣象條件之間的耦合關(guān)系［7-8］；地理科學分析PM2.5濃度的時間變化、空間分布、濃度預測、濃度高低的地理條件及其與產(chǎn)業(yè)、交通、人口等因素的相關(guān)性等方面，各學科之間對PM2.5的研究界線不明顯，互有交叉融合［9-12］。我國城市的PM2.5濃度具有較明顯的時空變化規(guī)律［11-12］，深入研究這些規(guī)律，對分析、預測和治理PM2.5具有積極意義。

文章采用時間序列分析、相關(guān)分析、空間自相關(guān)和空間插值等方法，分析了山東省2014、2015兩年間PM2.5濃度的時間規(guī)律和空間格局，分析了氣象因素對PM2.5濃度的影響。在此研究中，污染物監(jiān)測數(shù)據(jù)來源于中國環(huán)境監(jiān)測總站的城市空氣質(zhì)量實時發(fā)布平臺，匯總了山東省17個地級市于2014～2015年共728 d的空氣污染物監(jiān)測值；氣象數(shù)據(jù)由地面氣象站通過氣象數(shù)據(jù)國際自由交換系統(tǒng)提供，其中，德州市2014年4月20、21兩天的O3數(shù)據(jù)缺失，在進行PM2.5與SO2、CO、NO2、O3、PM10的相關(guān)分析時舍去了這兩天的污染物監(jiān)測數(shù)據(jù)。

1 山東省PM2.5濃度的時間變化規(guī)律

1.1 山東省PM2.5濃度的時間變化規(guī)律

1.1.1 PM2.5濃度的月度變化規(guī)律

2014年山東省PM2.5平均濃度為81.0μg/m3，明顯高于環(huán)境空氣質(zhì)量標準（GB 3095—2012）規(guī)定的二級濃度限值75μg/m3［13］。在2014年高于75μg/m3的月份有1、2、3、4、10、11、12等7個月份，其余5個月份低于75μg/m3，其中1月份的PM2.5濃度最高，達到了137.1μg/m3，而10月份的最低，為57.8μg/m3。2015年山東省PM2.5平均濃度為76.2μg/m3，略高于二級濃度限值75μg/m3。2015年高于75μg/m3的月份有1、2、3、10、11、12等6個月份，其余6個月份低于75μg/m3，其中12月份PM2.5濃度最高，為137.0μg/m3。這兩年山東省PM2.5濃度變化相對不大，2015年比2014年略有下降。冬季PM2.5濃度相對較高，夏季則相對較低，呈U型曲線，兩年中PM2.5日平均濃度最高的是2015年12月24日，達到了313.6μg/m3。主要原因是冬季燃煤取暖量增加，向大氣中排放的顆粒物增多。從氣象條件來看，冬季對流層大氣層結(jié)相對穩(wěn)定，在沒有冷空氣到來的情況下，加之城市建設(shè)高樓林立，污染物擴散速度較慢；而夏季，氣溫高、空氣對流旺盛，降水較多，促進了顆粒物的沉降。具體結(jié)果如圖1所示。

1.1.2 PM2.5濃度的時間序列自相關(guān)性

自相關(guān)分析是常用的一種時間序列分析方法［14］。自相關(guān)是指時間序列前后數(shù)據(jù)之間的相關(guān)關(guān)系，度量這種相關(guān)程度的指標為自相關(guān)系數(shù)［15］。自相關(guān)函數(shù)是理論意義上的，是針對時間序列總體而言的，在具體處理時經(jīng)常使用的是基于樣本路徑的樣本自相關(guān)函數(shù)［14］（auto-correlation function，ACF）。時間序列的樣本ACF的計算式（1）為

式中：rk為時間序列自相關(guān)系數(shù)；t為時序；k為時滯；n為樣本長度；xt為第t天的PM2.5濃度值；是xt的平均值。

采用MATLAB R2014a統(tǒng)計工具箱對山東省2014、2015年共728 d的PM2.5平均濃度進行自相關(guān)分析，揭示山東省PM2.5濃度變化的時間序列特征，自相關(guān)系數(shù)結(jié)果如圖2所示。查相關(guān)系數(shù)的臨界值表，在顯著水平α＝0.001下，自由度為700，相關(guān)系數(shù)檢驗的臨界值為r0.001＝0.124，一階自相關(guān)系數(shù)為0.6896，1～20階自相關(guān)系數(shù)均大于該臨界值，可見山東省PM2.5濃度隨時間變化的自相關(guān)性非常強，時間越近相關(guān)性越大，山東省每天PM2.5的平均濃度與其以前1～20 d的濃度相關(guān)。

圖1 山東省PM2.5濃度—時間變化規(guī)律圖

圖2 2014～2015年間山東省PM2.5濃度變化的自相關(guān)系數(shù)圖

1.2 山東省PM2.5濃度的空間分布規(guī)律

1.2.1 空間自相關(guān)

空間自相關(guān)的主要思想源于Tobler提出的地理學第一定律［16］，由時間序列自相關(guān)發(fā)展而來，是將一維自相關(guān)發(fā)展到空間二維自相關(guān)［14］。空間自相關(guān)是判斷某一地理要素的屬性值（指PM2.5濃度）與其相鄰空間的屬性值是否相關(guān)的指標，空間自相關(guān)系數(shù)分為正相關(guān)和負相關(guān)兩類，正相關(guān)表明某地理單元的屬性值（PM2.5濃度）變化與其相鄰空間單元的屬性值具有相同變化趨勢，負相關(guān)則表明兩者變化趨勢相反［15］。

全局空間自相關(guān)一般用Moran指數(shù)和Geary系數(shù)測度［15］，采用全局Moran指數(shù)I，其計算式（2）為

式中：xi、xj分別為位置（區(qū)域）i、j的觀測值；wij表示空間權(quán)重矩陣，文章采用queen鄰接方法確定空間權(quán)重矩陣，當i與j相鄰時取1，不相鄰時則取0是xi的平均值。

全局Moran指數(shù)I的取值范圍為［-1，1］，取值［0，1］表示正相關(guān)，取值為0表示不相關(guān)，取值［-1，0］表示負相關(guān)，其絕對值越接近于1，相關(guān)性越強［15，17-18］。對于Moran指數(shù)I，通常用標準化統(tǒng)計量Z檢驗n個區(qū)域的某一屬性值（PM2.5濃度值）是否存在空間自相關(guān)性，Z值計算式（3）［17］表示為

式中：E（I）為均值；VAR（I）為方差，兩者都是理論上的均值和方差［15］。其顯著性水平可以由標準化Z值的P值檢驗來確定［17］。通過計算Z值的P值，并將其與顯著性水平α作比較，決定拒絕或接收零假設(shè)［15］。若P值小于給定的顯著性水平α，則拒絕零假設(shè)，判定計算的區(qū)域?qū)傩灾荡嬖诳臻g自相關(guān)性；否則接受零假設(shè)［15，17-18］，判定屬性值的空間相關(guān)性不顯著。在實際問題分析中，通常將顯著性水平α定為0.05。

全局空間自相關(guān)只能從總體上反映地理要素之間是否存在空間聯(lián)系［15，17］。若要進一步判斷哪些地區(qū)的PM2.5濃度值存在空間集聚特征，需要采用局部空間自相關(guān)方法。常用的局部空間自相關(guān)分析方法有3種為空間聯(lián)系的局部指標LISA（Local Indicators of Spatial Association）、Moran散點圖和G統(tǒng)計［15］。文章采用局部Moran指數(shù)Ii，局部Moran指數(shù)Ii也要以標準化統(tǒng)計量Z值檢驗其顯著性。局部Moran Ii指數(shù)的計算式（4）為

1.2.2 PM2.5濃度的空間自相關(guān)

文章運用ArcGIS 10.2軟件對山東省2014、2015兩年中17個設(shè)區(qū)城市的PM2.5年平均濃度、4個季節(jié)平均濃度分別進行了全局空間自相關(guān)分析，結(jié)果見表1。

由表1可見，山東省2015年全年平均和春、夏、秋季節(jié)平均的PM2.5濃度的全局自相關(guān)系數(shù)Moran的I比2014年有所提高，冬季則稍有下降。各個時間段PM2.5濃度的全局自相關(guān)系數(shù)的Z值均大于2.58，對應(yīng)的P值均小于0.01，因此，其自相關(guān)性均達到了0.01的顯著水平，這表明山東省區(qū)域PM2.5的濃度變化存在明顯空間集聚特征。

采用GeoDa軟件對17個設(shè)區(qū)城市PM2.5濃度分布進行局部空間自相關(guān)分析，結(jié)果如圖3所示。

由圖3可見，煙臺和青島在各個時間段都是低低集聚，在夏季則增加了濰坊；高集聚主要分布在濟南，夏季則增加了濱州和德州2市，2014年冬季高高集聚由濟南變?yōu)闈鷮帲坏透呒壑怀霈F(xiàn)在魯中的泰安；其他設(shè)區(qū)城市PM2.5濃度值的局部空間自相關(guān)性不顯著。

表1 山東省PM2.5濃度的全局空間自相關(guān)系數(shù)

圖3 2014、2015年中山東省PM2.5濃度的局部空間自相關(guān)性圖

1.2.3 山東省PM2.5濃度的空間格局

空間插值是地統(tǒng)計學的主要研究內(nèi)容之一。城市是人口、經(jīng)濟、交通、建筑的集聚地，往往成為PM2.5質(zhì)量濃度的高值區(qū)，大多數(shù)的空氣質(zhì)量監(jiān)測站集中分布在城市，而廣大鄉(xiāng)村地區(qū)則少有分布。有研究指出，區(qū)域?qū)用娌逯档木缺冗b感反演要更加精確［11，19］。PM2.5濃度的空間分布是連續(xù)的，既具有隨機性又具有結(jié)構(gòu)性，各個城市的PM2.5濃度值大小看似是隨機變量，其實其濃度值存在一定的空間變化規(guī)律，上述的空間自相關(guān)性分析表明山東也有此規(guī)律。不同地區(qū)PM2.5濃度的大小是各個地區(qū)自然地理（地形、氣候等）、人文經(jīng)濟地理（城市規(guī)模、產(chǎn)業(yè)結(jié)構(gòu)及布局、人口等）條件的綜合反映，存在空間相互作用。傳統(tǒng)統(tǒng)計學是基于樣本之間的獨立性假設(shè)進行分析［15］，不考慮樣本間的空間相互作用，不適合應(yīng)用于PM2.5濃度空間分布研究?？死锝鸩逯捣椒ㄟm用于研究區(qū)域化變量存在空間相關(guān)性的問題，其實質(zhì)是利用區(qū)域化變量的實測數(shù)據(jù)和變異函數(shù)的結(jié)構(gòu)特點，對研究區(qū)域內(nèi)未知點的屬性值進行無偏、最優(yōu)估計［11，15］。采用ArcGIS 10.2軟件的地統(tǒng)計分析工具箱進行插值分析，經(jīng)驗證和交叉驗證多次實驗，根據(jù)誤差最小原則，應(yīng)用普通克里金法對山東省2014、2015年的PM2.5濃度進行空間插值分析，具體結(jié)果如圖4所示。

圖4 山東省2014、2015年中PM2.5濃度空間插值圖/（μg·m-3）

由圖4可見，山東省PM2.5濃度的空間分布格局總體上是西高東低，表現(xiàn)出從沿海到內(nèi)陸逐漸升高的趨勢，梯度變化明顯；從盛行風向分析，冬季山東省以西北風為主導，在自身排放細顆粒物的同時，從上風向的京津冀地區(qū)也攜帶了大量的細顆粒物過來，加之魯中有山地的阻擋，一定程度上加劇了山東省西北部地區(qū)細顆粒物濃度的累積；到了夏季，盛行風向變成了東南季風，風從海洋吹來，空氣中污染物很少，在一定程度上稀釋了山東省自身排放的細顆粒物濃度；泰安周圍區(qū)域為PM2.5濃度分布的谷地，原因是泰安地區(qū)地形為山地，較周邊城市高差顯著，對來自西北方向的細顆粒物有阻擋作用。東部的威海市與西北部的德州市PM2.5濃度值相差近3倍；在季節(jié)分布來看，2015年春季和夏季PM2.5濃度值均較2014年有明顯改善，夏季全省PM2.5濃度值基本上低于80μg/m3，而冬季僅有膠東半島濃度值較低。

2 PM2.5濃度與其他污染物的相關(guān)性

空氣質(zhì)量監(jiān)測指標（SO2、CO、NO2、O3、PM10、PM2.5）之間相互影響，文中通過計算空氣質(zhì)量監(jiān)測指標之間的相關(guān)系數(shù)測定其密切程度，相關(guān)系數(shù)由式（5）計算為

式中：rxy為要素x與y之間的相關(guān)系數(shù)和分別表示兩個要素樣本值的平均值。rxy值介于［-1，1］，rxy＞0，表示正相關(guān)，越接近于1相關(guān)性越強；rxy＜0，表示負相關(guān)，越接近于 -1表示負相關(guān)性越強；rxy值接近于0，表示兩要素相關(guān)性越弱［15］。

應(yīng)用MATLAB R2014a軟件對2014、2015兩年共728 d的PM2.5質(zhì)量濃度值與PM10、SO2、CO2、O3、NO2的質(zhì)量濃度值進行相關(guān)系數(shù)分析。表2中上三角為相關(guān)系數(shù)，下三角為相關(guān)系數(shù)對應(yīng)的P值。

由表2中可見，各種污染物相關(guān)系數(shù)的P值非常小，幾乎接近于0，這表明各種污染物的兩兩相關(guān)性均達到了極高的顯著性水平。其中，PM2.5濃度與PM10、SO2、CO、NO2等污染物濃度存在極高的正相關(guān)性，PM2.5與PM10的相關(guān)系數(shù)最高，與O3的濃度值存在負相關(guān)性。PM2.5主要由人為產(chǎn)生，在污染物排放量大和不利于擴散的氣象條件下，PM2.5易成首要污染物［20］。O3以自然來源為主，在高溫和強光照條件下，濃度易升高［20］。O3與PM2.5濃度上升的條件相反，當PM2.5濃度高時，O3污染輕；當O3濃度高時，PM2.5污染低［20］。從空間上來看，城鄉(xiāng)差別明顯，城區(qū)O3濃度高于郊區(qū)和農(nóng)村的O3濃度［21］。

表2 PM2.5濃度與其他污染物的相關(guān)系數(shù)

3 PM2.5濃度與氣象因子的相關(guān)分析

3.1 PM2.5濃度與風向的關(guān)系

選取濃度值最高的2015年12月24日及其前后兩天的數(shù)據(jù)分析PM2.5濃度與風向的關(guān)系。首先對山東省88個監(jiān)測站的PM2.5濃度數(shù)據(jù)進行空間插值分析，結(jié)果如圖5所示。

查閱這3 d的氣象資料（見表3）可知，23日全省以西北風和東北風為主，風力為3～4級微風，PM2.5濃度分布西北高、東南低，梯度變化明顯。24日，風向為北風（東北風、西北風）和南風（東南風、西南風）各占一半，在23日PM2.5濃度分布的基礎(chǔ)上，PM2.5細顆粒物繼續(xù)由西北向東南擴散，加之部分地區(qū)為東南風，PM2.5濃度分布梯度變化與23日比不明顯，且出現(xiàn)局部高值區(qū)（菏澤、德州等）。25日，全省為西南風和東南風各占一半，PM2.5濃度分布梯度變化不明顯，污染明顯減輕，并呈現(xiàn)西北、東南方向的條帶狀分布。

3.2 PM2.5濃度與風速、溫度、氣壓的關(guān)系

PM2.5粒徑小，可長期懸浮于大氣中，其濃度變化必然受到各種氣象條件的影響。考慮到數(shù)據(jù)的可獲取性，選取濟南市2015年1、4、7、10等4個月份的溫度、風速、大氣壓、相對濕度和露點溫度等氣象數(shù)據(jù)（采樣數(shù)據(jù)均為距地面2 m高位置）。為便于比較分析，計算的各月份各指標的平均值見表4。采用式（6）對PM2.5與各氣象因子進行相關(guān)分析，計算的相關(guān)系數(shù)見表5。

圖5 山東省2015年12月23～25日PM2.5濃度的空間格局圖

表3 23～25日山東省各市主要風向

表4 濟南市2015年1、4、7、10月份的溫度、風速、大氣壓、相對濕度和露點溫度的平均值

表5 PM2.5濃度與氣象因子的相關(guān)系數(shù)

由表4和5可見，1月份PM2.5濃度最高，10 月份較高，4和7月相對較低。各月份平均風速差別較小，PM2.5濃度與風速呈現(xiàn)負相關(guān)性，當PM2.5濃度較高時（1和10月份），PM2.5濃度與風速的相關(guān)性不顯著；當PM2.5濃度較低時（4和7月份），PM2.5濃度與風速的相關(guān)性分別達到了0.01和0.001的顯著性水平。PM2.5濃度與溫度的相關(guān)性在1、4和7月份均不顯著，在10月份兩者的相關(guān)性達到了0.05的顯著性水平。當大氣壓最高時（1月份），PM2.5濃度與大氣壓呈負相關(guān)，且達到了0.05的顯著性水平；當大氣壓最低時（7月份），PM2.5濃度與大氣壓呈正相關(guān)，也達到了0.05的顯著性水平。當相對濕度較低時（1和10月份），PM2.5濃度與相對濕度呈正相關(guān)，達到了0.05的顯著性水平，當相對濕度較高時（4和7月份），PM2.5濃度與相對濕度的相關(guān)性不顯著。當露點溫度較高時（7月份），PM2.5濃度與露點溫度的相關(guān)性不顯著；當露點溫度較低時（1、4、10月份），PM2.5濃度與露點溫度的呈正相關(guān)性，且分別達到了0.01和0.001的顯著性水平，兩者的相關(guān)性極為顯著。

4 結(jié)論

通過上述研究可知：

（1）山東省PM2.5濃度時間變化規(guī)律明顯。2014～2015年山東省的PM2.5濃度日均變化曲線呈現(xiàn)冬季高、夏季低的U型狀態(tài)。PM2.5濃度時間序列自相關(guān)性較強，一階自相關(guān)系數(shù)為0.6896，1～20階自相關(guān)系數(shù)均大于臨界值0.124。山東省PM2.5濃度空間分布規(guī)律明顯，PM2.5濃度呈現(xiàn)由東南沿海向西北內(nèi)陸遞增趨勢。

（2）PM2.5的濃度與風速、風向、氣溫、大氣壓、相對濕度和露點溫度等氣象因素存在一定關(guān)系。在不同的月份其相關(guān)性也存在一定差異。

（3）PM2.5的濃度與NO2、SO2、CO、PM10、O3等的濃度相關(guān)性較強，對應(yīng)相關(guān)系數(shù)分別為0.7846、0.6735、0.9043、0.9214、-0.3122。在治理霧霾時，應(yīng)對空氣質(zhì)量監(jiān)測的6項指標綜合治理。

［1］ Tian S.，Pan Y.，Liu Z.，et al..Size-resolved aerosol chemical analysis of extreme haze pollution events during early 2013 in urban Beijing，China［J］.Journal of Hazardous Materials，2014（279）：452-460.

［2］ 朱增銀，李冰，趙秋月，等.對國內(nèi)外PM2.5研究及控制對策的研究與展望［J］.環(huán)境科技，2013（261）：70-74.

［3］ 杜金花，黃曉鋒，何凌燕，等.深圳市大氣細粒子（PM2.5）中汞的污染特征［J］.環(huán)境科學研究，2010，23（6）：667-673.

［4］ 馮茜丹，黨志，黃偉林.廣州市秋季PM2.5中重金屬的污染水平與化學形態(tài)分析［J］.環(huán)境科學，2008，29（3）：569-575.

［5］ 孟紫強，張全喜.沙塵暴PM2.5對大鼠肺細胞DNA損傷效應(yīng)［J］.中國公共衛(wèi)生，2006，22（12）：1458-1459.

［6］ 陳威，郭新彪，鄧芙蓉，等.大氣細顆粒物對A549細胞炎性因子分泌影響［J］.中國公共衛(wèi)生，2007，23（9）：1080-1081.

［7］ 楊洪斌，鄒旭東，汪宏宇，等.大氣環(huán)境中PM2.5的研究進展與展望［J］.氣象與環(huán)境學報，2012，28（3）：77-82.

［8］ 徐敬，丁國安，顏鵬，等.北京地區(qū)PM2.5的成分特征及來源分析［J］.應(yīng)用氣象學報，2007，18（5）：645-654.

［9］ 李松，羅緒強，李戀，等.基于GIS的中國PM2.5濃度的空間分布及影響因素分析［J］.水土保持通報，2015，35（4）：202-205，212.

［10］林巧鶯，陳永山.中國城市冬季PM2.5空間特征及其人為影響因子［J］.生態(tài)與農(nóng)村環(huán)境學報，2015，31（4）：460-465.

［11］王振山，方創(chuàng)琳，許光，等.2014中國城市PM2.5濃度的時空變化規(guī)律［J］.地理學報，2015，70（11）：1720-1734.

［12］王占山，李云婷，陳添，等.2013年北京市PM2.5的時空分布［J］.地理學報，2015，70（1）：110-120.

［13］GB 3095—2012.環(huán)境空氣質(zhì)量標準［S］.北京：中國環(huán)境科學出版社，2012.

［14］陳彥光.基于Matlab的地理數(shù)據(jù)分析［M］.北京：高等教育出版社，2012.

［15］徐建華.計量地理學［M］.北京：高等教育出版社，2006.

［16］Tobler W..A computer movie simulating urban growth in the Detroit region［J］.Economic Geography，1970，46（2）：240-243.

［17］曹蕾.區(qū)域生態(tài)文明建設(shè)評價指標體系及建模研究［D］.上海：華東師范大學，2014.

［18］閔曙輝，張郴，潘穎.基于ESDA-GIS的江蘇省13轄市旅游空間關(guān)聯(lián)分析［J］.旅游研究，2015，7（2）：52-57.

［19］Lee S.J.，Serre M.L.，Van Donkelaar A.，etal..Comparison of geostatistical interpolation and remote sensing techniques for estimating long-term exposure to ambient PM2.5 concentrations across the continental United States［J］.Environment Health Perspect，2012，120（12）：1727-1732.

［20］羅岳平，劉孟佳，甘杰，等.長株潭城市環(huán)境空氣中PM2.5和O3質(zhì)量濃度的相關(guān)性研究［J］.安全與環(huán)境學報，2015，15（4）：313-317.

［21］劉潔，張小玲，徐曉峰，等.北京地區(qū)SO2、NOx、O3和PM2.5變化特征的城郊對比分析［J］.環(huán)境科學，2008，29（4）：1059-1065.

Tem poral and spatial variation and its influencing factors of PM 2.5 concentration in Shandong Province

Fan Kaifang，Shan Baoyan*
（School of Surveying and Geoinformatics，Shandong Jianzhu University，Jinan 250101，China）

PM 2.5 is an important target of air quality monitoring，and understanding the temporal and spatial variation of PM2.5 is of great significance for governing gray haze.Based on 728 days'air quality monitoring data from 2014 to 2015 and relevantmeteorological data of 17 cities in Shandong province，using time series analysis，correlation analysis，spatial autocorrelation and spatial interpolation，we studied the spatial and temporal variation of the concentration of PM2.5 in Shandong Province and analyzed the influence of meteorological factors on PM2.5 concentration.The results show that the wintermonths presentwith a higher concentration，relatively lower in summermonths，showing a U-shape changes.The time series autocorrelation of PM2.5 concentration changes was strong，and first order autocorrelation coefficient reaches 0.6896.The concentration distribution of PM2.5 shows obvious spatial agglomeration characteristics.The direction ofwind has an obvious effect on the aggregation and diffusion of PM2.5 concentration；PM2.5 concentration and meteorological conditions show a certain correlation；in different months the correlation is significantly different. Meanwhile，the concentration of PM2.5 and SO2，CO，NO2，PM10 exists a highly positive correlation coefficientwhile a negative correlation with O3.In the process of governing the haze，we should carry on the comprehensivemanagement of various pollutants.

X513

A

1673-7644（2017）01-0039-08

2017-01-19

教育部人文社會科學研究規(guī)劃基金項目（12YJA790019）；住房和城鄉(xiāng)建設(shè)部研究開發(fā)項目（2012-K8-22）

范開放（1991-），男，在讀碩士，主要從事3S技術(shù)在城市建設(shè)及其管理中的應(yīng)用等方面的研究.E-mail：fankaifang_dili＠163.com

*：單寶艷（1970-），男，教授，碩士，主要從事GIS與區(qū)域規(guī)劃等方面的研究.E-mail：shan7066＠sdjzu.edu.cn