李娟　何克?！煶枯x　王春林　楊海波

摘 要：PM2.5濃度對空氣的質(zhì)量有著重要影響，對PM2.5濃度及其時(shí)空變化規(guī)律的研究，對由于霧霾引起的空氣污染的治理有著重要意義。本文基于華北地區(qū)59個城市2014年10月-2015年3月、2015年10月-2016年3月、2016年10月-2017年3月的實(shí)時(shí)PM2.5濃度數(shù)據(jù)與MODIS氣溶膠數(shù)據(jù)，利用空間插值、相關(guān)分析、城郊緩沖差異分析和時(shí)間序列分析等方法研究華北地區(qū)主要城市PM2.5濃度及時(shí)空變化的規(guī)律。結(jié)果表明：在空間分布上以華北地區(qū)的中東部或中西部為中心向四周擴(kuò)散，空間分布特征明顯。同時(shí)小型、中型和大型城市PM2.5濃度城郊差異各不相同，小型城市郊區(qū)PM2.5濃度大于市區(qū)，中型城市普遍為郊區(qū)PM2.5濃度小于市區(qū)，大型城市北京市郊區(qū)PM2.5濃度遠(yuǎn)大于市區(qū)。在時(shí)間分布上高發(fā)與12月和1月，各市的月插值PM2.5濃度與MODIS月融合PM2.5濃度因月份的不同其相關(guān)性各異，同時(shí)PM2.5站點(diǎn)統(tǒng)計(jì)數(shù)據(jù)以12月份為軸，研究時(shí)間段內(nèi)PM2.5濃度分布近似于正態(tài)分布。

關(guān)鍵詞：PM2.5濃度；MODIS氣溶膠數(shù)據(jù)；時(shí)空變化；估算模型；正態(tài)分布

Abstract：The concentration of PM2.5 has a significant impact on the quality of the air. It has a great significance of the research on the variation of PM2.5 concentration and its temporal， spatial distribution and the control of air pollution caused by smog and haze. This article is based on real-time PM2.5 concentration data and MODIS aerosol data from 59 cities in North of China from 2014 to 2015， 2015 to 2016， 2016 to 2017， from October to December of the previous year and January to March of next year， spatial interpolation， correlation analysis， suburban buffer difference analysis and time series analysis were used to study the variation of PM2.5 concentration and time-space in major cities of North China. The results show： in terms of spatial distribution， the central and eastern parts or the central and western parts of North China as the center and spread to four weeks， spatial distribution is obvious. At the same time， suburban differences in PM2.5 concentrations vary among small， medium and large cities. The concentration of PM2.5 in the suburbs of small cities is greater than that in the urban areas， and the concentration of PM2.5 in the suburbs of the medium-sized cities is generally less than that in the urban areas. The concentrations of PM2.5 in the suburbs of Beijing is much higher than those in the urban areas. In the time distribution， it was high in December and January. The correlation between PM2.5 concentration interpolation data and remote sensing fusion data in different months is different. Meanwhile， PM2.5 station statistics in December as the axis， during the study period， the PM2.5 concentration distribution is similar to the normal distribution.

Key Words： PM2.5 concentration ；MODIS aerosol data ；Time and space changes； Estimate model； Normal distribution

1.引言

在空氣動力學(xué)中，將微粒直徑小于2.5μm的粒子，稱為PM2.5氣溶膠粒子，簡稱為PM2.5。PM2.5對太陽的輻射具有很強(qiáng)的散射作用，當(dāng)其濃度增加時(shí)大氣能見度將急劇下降[1]。近年來，隨著中國經(jīng)濟(jì)的快速發(fā)展，環(huán)境污染問題越來越嚴(yán)重，PM2.5污染尤為突出，成為危害人類身體健康的主要因素之一[2]，且PM2.5具有較強(qiáng)的時(shí)空性，不同地區(qū)不同季節(jié)PM2.5的濃度變化各有不同。因此，越來越多的學(xué)者研究PM2.5對人體健康的影響及其時(shí)空變化。聶晨暉[3]等對杭州市2015年P(guān)M2.5濃度進(jìn)行了時(shí)空變化分析，證實(shí)PM2.5濃度的空間分布與城鎮(zhèn)化格局相似，PM2.5濃度分布的季節(jié)性特征為：冬季>春季>秋季>夏季。張文[4]等對2013年京津冀地面測站PM2.5濃度數(shù)據(jù)和氣溶膠光學(xué)厚度數(shù)據(jù)（AOT）進(jìn)行線性擬合后，發(fā)現(xiàn)利用AOT數(shù)據(jù)的線性擬合模型能夠有效地估算出PM2.5濃度。

本文以華北地區(qū)為研究對象，華北地區(qū)作為北方經(jīng)濟(jì)發(fā)展的重要樞紐，集農(nóng)業(yè)和工業(yè)于一體，城市人口流動較大，城市化速度較快，特別是京津冀城市群和中原城市群。研究選取的時(shí)間段為2014年到2017年，以霧霾特征明顯的月份為研究時(shí)間，即第一年10月-12月份，第二年1月-3月份。所得的實(shí)時(shí)PM2.5濃度數(shù)據(jù)為環(huán)?？傉救掌骄鶖?shù)據(jù)，但各站點(diǎn)的觀測數(shù)據(jù)只能表示其周圍一定半徑范圍內(nèi)的PM2.5濃度的分布情況，不能反映整個區(qū)域內(nèi)PM2.5濃度的分布情況，所以要通過空間插值的方法對連續(xù)空間的PM2.5濃度進(jìn)行模擬[5]。先利用14-15年的實(shí)時(shí)PM2.5濃度插值數(shù)據(jù)，選取代表性城市進(jìn)行多圈緩沖分析，研究其城郊差異變化特征。同時(shí)統(tǒng)計(jì)各省會的3年實(shí)時(shí)PM2.5濃度數(shù)據(jù)的月平均值，分析其時(shí)間特征；分析2014年12月份華北地區(qū)PM2.5濃度的空間分布特征。最后，以北京市和鄭州市為例，在實(shí)時(shí)PM2.5濃度數(shù)據(jù)和MODIS氣溶膠光學(xué)厚度數(shù)據(jù)的基礎(chǔ)上，分月建立線性擬合模型，并進(jìn)行PM2.5濃度估算和精度分析。

2.數(shù)據(jù)與方法

2.1研究區(qū)域

華北地區(qū)一般指秦嶺-淮河線以北，長城以南的中國廣大區(qū)域。包括：北京市、天津市、河北省、山西省和內(nèi)蒙古自治區(qū)中部。華北地區(qū)是中華民族的發(fā)源地之一，中國的政治文化中心。內(nèi)蒙古區(qū)域呈狹長形狀，區(qū)域較大，且霧霾特征不明顯，東西區(qū)霧霾的數(shù)據(jù)差異較大。華北平原為中國霧霾高發(fā)地，所以我們在研究時(shí)去掉了內(nèi)蒙古自治區(qū)，加入了河南省。故本文所研究的華北地區(qū)包括北京市、天津市、河南省、河北省、山東省和山西省。

2.2數(shù)據(jù)來源

本文使用了2014年-2015年、2015年-2016年、2016年-2017年每上一年的10月、11月、12月和下一年1月、2月、3月的實(shí)時(shí)PM2.5濃度數(shù)據(jù)及MODIS氣溶膠數(shù)據(jù)。該研究時(shí)間段內(nèi)霧霾特征明顯，易于挖掘數(shù)據(jù)的潛在信息。實(shí)時(shí)PM2.5數(shù)據(jù)來源于網(wǎng)站（https：//www.aqistudy.cn/historydata）PM2.5歷史數(shù)據(jù)，該P(yáng)M2.5濃度數(shù)據(jù)是根據(jù)當(dāng)天環(huán)?？傉久啃r(shí)數(shù)據(jù)計(jì)算求平均的結(jié)果。選取了華北地區(qū)59個市的PM2.5濃度數(shù)據(jù)。氣溶膠數(shù)據(jù)來源于美國國家航空航天局（NASA）搭載在Terra和aqua衛(wèi)星上的MODIS傳感器遙感影像產(chǎn)品MOD04_3K數(shù)據(jù)，空間分辨率10km，波長550nm。該數(shù)據(jù)具有分辨率高、多光譜等特點(diǎn)。

2.3數(shù)據(jù)處理

2.3.1實(shí)時(shí)PM2.5濃度數(shù)據(jù)處理

首先利用ArcGIS軟件中ArcToolBox工具箱的普通克里金插值（kriging）方法，對2014年10月至12月和2015年1月至3月華北地區(qū)59個市實(shí)時(shí)PM2.5濃度數(shù)據(jù)按天進(jìn)行空間插值?？死锝鸩逯捣椒ǖ慕Y(jié)果往往可以更接近真實(shí)值，所以，在模擬真實(shí)效果的情況下，克里金插值方法是最好的選擇[6]。然后對每張插值后柵格數(shù)據(jù)按市進(jìn)行分區(qū)統(tǒng)計(jì)，最后得出各市的PM2.5濃度平均值。

2.3.2 MODIS氣溶膠數(shù)據(jù)預(yù)處理

首先使用HEG軟件選出MOD04_3K數(shù)據(jù)的Optical_Depth_Land_And_Ocean波段，并進(jìn)行數(shù)據(jù)格式轉(zhuǎn)換和投影轉(zhuǎn)換。其二使用MRT軟件進(jìn)行批量的裁剪處理，剔除灰度值為空的柵格影像。然后在GIS平臺ArcGIS的支持下，將一天內(nèi)不同時(shí)段的MOD04_3K數(shù)據(jù)進(jìn)行融合處理，由于收集的實(shí)時(shí)PM2.5數(shù)據(jù)為當(dāng)天環(huán)?？傉久啃r(shí)數(shù)據(jù)的平均值，因此這里的影像融合取平均值。最后對2014年10月至12月和2015年1月至3月融合后的MOD04_3K數(shù)據(jù)按市進(jìn)行分區(qū)統(tǒng)計(jì)，得出各市范圍內(nèi)柵格像元的平均值。

3.霧霾時(shí)空變化

3.1時(shí)間趨勢分析

統(tǒng)計(jì)華北地區(qū)六個省會城市，即北京市、天津市、鄭州市、石家莊市、濟(jì)南市、太原市，三年的實(shí)時(shí)PM2.5濃度數(shù)據(jù)的月平均值。統(tǒng)計(jì)結(jié)果如下圖：

2014年10月-12月時(shí)間段內(nèi)，北京市、石家莊市、太原市的月平均PM2.5濃度均呈下降趨勢；天津市、鄭州市、濟(jì)南市的月平均PM2.5濃度均先上升后下降；六個城市的月平均PM2.5濃度在次年1月份都出現(xiàn)回升現(xiàn)象。對比分析2015年及2016年1月-3月和10月-12月的數(shù)據(jù)，可發(fā)現(xiàn)，六個城市的月平均PM2.5濃度存在兩個高峰態(tài)。高峰期均為1月和12月。這說明全年中霧霾天氣高發(fā)期為1月和12月，冬季霧霾現(xiàn)象嚴(yán)重。2017年1月-3月時(shí)間段內(nèi)，六個城市的月平均PM2.5濃度均呈下降趨勢。

隨著時(shí)間的推移，六個城市的月平均PM2.5濃度的峰值逐漸增大，同時(shí)說明了霧霾天氣在逐漸加劇。這與社會的經(jīng)濟(jì)發(fā)展需求緊密相關(guān)，霧霾天氣的治理迫在眉睫。石家莊市2016年-2017年的月平均PM2.5濃度峰值其他時(shí)間段增長劇烈，其他城市則為穩(wěn)定增長。這一現(xiàn)象是由當(dāng)時(shí)的氣象條件、大氣環(huán)流的異常和石家莊市特殊的地貌特征共同作用造成的[7-8]。將石家莊市三個時(shí)間段的月平均PM2.5濃度與其他城市對比分析，發(fā)現(xiàn)石家莊市的月平均PM2.5濃度普遍高于其他省會城市的。經(jīng)查閱相關(guān)資料得知，這一現(xiàn)象源于石家莊的冬天天氣干燥，降雨量小，導(dǎo)致大氣污染物沉降難；大城市居住部落的進(jìn)程推動了熱島效應(yīng)；以煤炭為主要能源；植被面積嚴(yán)重減少[8]。

由以上分析可總結(jié)得，霧霾天氣的高發(fā)期為12月和1月，冬季霧霾現(xiàn)象明顯。各城市的霧霾現(xiàn)象隨時(shí)間逐漸加劇，且較其他省會城市，石家莊市月平均PM2.5濃度普遍高。

3.2城郊差異分析

本文從華北地區(qū)59個城市中，按地區(qū)的面積，選取1個大型城市（北京市）、5個中型城市（鄭州市、唐山市、石家莊市、濟(jì)南市、保定市）2個小型城市（新鄉(xiāng)市、衡水市）進(jìn)行城郊差異分析。時(shí)間段選取2014年的12月，屬于霧霾高發(fā)期。基于ArcGIS平臺，將12月份的實(shí)時(shí)PM2.5濃度的插值數(shù)據(jù)進(jìn)行融合處理，然后以市中心為原點(diǎn)進(jìn)行多圈緩沖分析。其中大型城市的緩沖半徑設(shè)為10、20、30、40、50，中型城市的緩沖半徑設(shè)為8、16、24、32、40，小型城市的緩沖半徑設(shè)為6、12、18、24、30，緩沖半徑一個單位為1km。如圖3所示。

各城市不同緩沖半徑內(nèi)PM2.5濃度平均值統(tǒng)計(jì)結(jié)果如上圖所示。所選的小型城市中衡水市和新鄉(xiāng)市的PM2.5濃度值均隨緩沖半徑R的增大而增大，即沿市中心向外擴(kuò)展PM2.5濃度值越來越大，霧霾天氣越來越嚴(yán)重。說明衡水市和新鄉(xiāng)市的郊區(qū)PM2.5濃度大于市區(qū)的PM2.5濃度。新鄉(xiāng)市與鄭州市相鄰，鄭州市的霧霾天氣對新鄉(xiāng)市郊區(qū)有一定的影響。衡水市周邊的幾個市的霧霾天氣均較嚴(yán)重，這使得衡水市郊區(qū)的PM2.5濃度增大。

所選的5個中型城市中，除唐山市以外的四個城市的PM2.5濃度均隨緩沖半徑R的增大而減小，即沿市中心向外擴(kuò)展PM2.5濃度值越來越小。這表明保定市、石家莊市、唐山市和鄭州市的市區(qū)PM2.5濃度大于郊區(qū)的PM2.5濃度。這一現(xiàn)象與市區(qū)和郊區(qū)的經(jīng)濟(jì)發(fā)展程度、人口居住情況和廢氣排放量等影響霧霾天氣的因素有著密切的聯(lián)系[9]。而唐山市PM2.5濃度隨緩沖半徑的變化規(guī)律與小型城市相同。

所選的大型城市北京市的PM2.5濃度值均隨緩沖半徑R的增大而增大，特別是30公里-40公里范圍內(nèi)的PM2.5濃度顯著增大。造成之一現(xiàn)象的原因是，北京市在政策和經(jīng)濟(jì)的共同作用下產(chǎn)生了郊區(qū)化，進(jìn)而人口和居住郊區(qū)化現(xiàn)象明顯，產(chǎn)業(yè)結(jié)構(gòu)也發(fā)生變化。北京市的第二、第三產(chǎn)業(yè)由中心區(qū)和遠(yuǎn)郊區(qū)向近郊區(qū)聚集[10]。而天津市PM2.5濃度值均隨緩沖半徑變化的幅度非常小，說明天津市霧霾現(xiàn)象的城郊差異不明顯。

3.3空間分布分析

將華北地區(qū)2014年12月份的實(shí)時(shí)PM2.5插值數(shù)據(jù)進(jìn)行融合后，得到PM2.5濃度分布圖如下所示。

從圖中可以看出，2014年和2015年12月霧霾天氣的空間分布大致相同。河北省的保定市、石家莊市、衡水市等、天津市、和山東省的臨沂市、棗莊市、萊蕪市等地的PM2.5濃度明顯高于其余地方。保定市區(qū)的PM2.5濃度最高，華北地區(qū)邊界處PM2.5濃度普遍低。以華北地區(qū)的中東部為中心，向外延伸PM2.5濃度逐漸減小。而2016年2月是以華北地區(qū)的中西部為中心，向外延伸PM2.5濃度逐漸減小。

4. PM2.5濃度模型

本文以北京市和鄭州市為例，選取時(shí)間段為2014年-2015年，以天為單位對實(shí)時(shí)PM2.5插值數(shù)據(jù)和MODIS氣溶膠融合數(shù)據(jù)分月建立線性模型。回歸方程的形式為：，其中 代表MODIS氣溶膠融合數(shù)據(jù)值， 代表對應(yīng)的實(shí)時(shí)PM2.5插值數(shù)據(jù)值，線性關(guān)系及參數(shù)如下表格所示。

相關(guān)系數(shù) 表示實(shí)時(shí)PM2.5插值數(shù)據(jù)和MODIS氣溶膠融合數(shù)據(jù)的線性相關(guān)程度， 越接近1，說明數(shù)據(jù)的擬合度越大。由擬合結(jié)果可得到：北京市2014年-2015年實(shí)時(shí)PM2.5插值數(shù)據(jù)和MODIS氣溶膠融合數(shù)據(jù)的相關(guān)系數(shù)在0.0122到0.4541之間，則兩者有較強(qiáng)的相關(guān)性。其中十月份數(shù)據(jù)的擬合度最高， =0.4541，十一月份數(shù)據(jù)為負(fù)相關(guān)，其余月份的數(shù)據(jù)均為正相關(guān)。

鄭州市2014年-2015年實(shí)時(shí)PM2.5插值數(shù)據(jù)和MODIS氣溶膠融合數(shù)據(jù)的相關(guān)系數(shù)在0.0025到0.7398之間，則兩者有較強(qiáng)的相關(guān)性。其中三月份數(shù)據(jù)的擬合度最高，擬合度為0.7398。十月份、十一月份和十二月份數(shù)據(jù)為負(fù)相關(guān)，一月份、二月份、三月份的數(shù)據(jù)均為正相關(guān)。

5.結(jié)論

通過對華北地區(qū)不同年份、不同月份以及不同空間的PM2.5濃度進(jìn)行對比分析，得知華北地區(qū)PM2.5濃度具有較強(qiáng)的時(shí)空變化規(guī)律。

時(shí)間分布規(guī)律：具有較強(qiáng)的季節(jié)性，冬季和春季PM2.5濃度相差較大，每年P(guān)M2.5濃度基本以十二月份為軸，呈近似正態(tài)分布，以及PM2.5濃度插值數(shù)據(jù)與融合數(shù)據(jù)的相關(guān)性因月份不同各異，隨月份的變化呈現(xiàn)出漸變性。

空間分布規(guī)律：整個華北地區(qū)PM2.5濃度的空間分布規(guī)律明顯，以中東部為中心向四周擴(kuò)散，遍及整個華北地區(qū)，且不同城市類型PM2.5濃度城郊差異明顯，小型大型城市郊區(qū)PM2.5濃度大于市區(qū)PM2.5濃度，中型城市郊區(qū)PM2.5濃度小于市區(qū)PM2.5濃度。

華北地區(qū)PM2.5濃度出現(xiàn)以上時(shí)空變化規(guī)律的原因眾多，隨時(shí)間的變化規(guī)律與不同年份不同月份氣候條件等的變化有很大的關(guān)系，隨空間的變化規(guī)律與華北地區(qū)不同城市群的形成先后、城市的發(fā)展速度、地理位置上差異以及產(chǎn)業(yè)結(jié)構(gòu)的演變有著很大的關(guān)系。

參考文獻(xiàn)：

[1]于鳳蓮，劉東賢，胡英.有關(guān)氣溶膠細(xì)粒子對城市能見度影響的研究.氣象科技，2002，30（6）：379-383.

[2]方冬青，魏永杰，黃偉，等.北京市2014年10月中霾污染特征及有機(jī)碳來源解析[J].環(huán)境科學(xué)研究，2016，29（1）：12-19.

[3] 聶晨暉，潘驍駿，金洪芳. 杭州地區(qū)2015年P(guān)M_（2.5）濃度時(shí)空變化特征分析[J]. 測繪通報(bào)，2016（11）：75-79.

[4] 張文，潘竟虎. 2013年大范圍霧霾期間京津冀PM_（2.5）質(zhì)量濃度遙感估算及時(shí)空變化的經(jīng)驗(yàn)正交函數(shù)分析[J]. 蘭州大學(xué)學(xué)報(bào)（自然科學(xué)版），2016，52（03）：350-356.

[5]陳強(qiáng)，梅琨，朱慧敏，蔡賢雷，張明華. 鄭州市PM_（2.5）濃度時(shí)空分布特征及預(yù)測模型研究[J]. 中國環(huán)境監(jiān)測，2015，31（03）：105-112.

[6]翟進(jìn)乾. 克里金（kriging）插值方法在煤層分布檢測中的應(yīng)用研究[D].太原理工大學(xué)，2008.

[7]袁東敏，馬小會. 2016年12月16～21日重度霾過程及大氣環(huán)流異常[J]. 氣候與環(huán)境研究，2017，22（06）：757-764.

[8] 孫梓瀅. 霧霾污染結(jié)構(gòu)特征及影響因素研究：以石家莊為例[J]. 產(chǎn)業(yè)與科技論壇，2017，16（08）：34-36.

[9]秦昊洋. 淺談霧霾成因的有效分析方法與應(yīng)用[J]. 陜西廣播電視大學(xué)學(xué)報(bào)，2014，16（01）：92-96.

[10]劉志軍. 郊區(qū)化背景下的北京市產(chǎn)業(yè)結(jié)構(gòu)演變的研究[D].首都經(jīng)濟(jì)貿(mào)易大學(xué)，2008.

鄭州大學(xué)創(chuàng)新創(chuàng)業(yè)項(xiàng)目（2017cxcy258）