我國空氣質(zhì)量區(qū)域劃分

（蘭州財經(jīng)大學，甘肅 蘭州 730030）

1 研究目的及意義

2020年國務(wù)院政府工作報告指出，接下來一年要打好藍天、碧水、凈土保衛(wèi)戰(zhàn)，實現(xiàn)污染防治攻堅戰(zhàn)階段性目標。提高生態(tài)環(huán)境治理成效，突出依法、科學、精準治污，深化重點地區(qū)大氣污染治理攻堅。隨著我國大氣區(qū)域性復(fù)合型污染特征日益明顯，區(qū)域性大氣灰霾頻繁發(fā)生，區(qū)域整體環(huán)境質(zhì)量惡化，嚴重威脅群眾健康，影響環(huán)境安全。建立大氣污染聯(lián)防聯(lián)控機制，顯得尤為迫切。研究空氣質(zhì)量空間分布規(guī)律和影響因素，對于預(yù)防空氣污染問題具有重要意義，空氣污染問題已經(jīng)從單一城市問題上升到了區(qū)域性問題。建立大氣污染聯(lián)防聯(lián)控機制，合理的區(qū)域劃分顯得尤為重要，但是目前針對空氣質(zhì)量的全國性區(qū)域劃分研究很少，大多數(shù)都是根據(jù)特定行政區(qū)劃分作為研究對象，并沒有真正意義上全國性區(qū)域劃分的研究；目前大多數(shù)研究對空氣質(zhì)量的影響因素忽視區(qū)域因素，但是每座城市會因為其特殊的地理位置，在大氣污染傳送網(wǎng)絡(luò)等多種因素的作用下或多或少的受到周邊城市的影響，在這種影響下，空氣污染問題已經(jīng)從單一城市問題上升到了區(qū)域性問題。研究空氣質(zhì)量的影響因素是一個多維問題，會受到多方面因素的影響，比如不同地區(qū)因自然因素和社會經(jīng)濟等因素的差別，其空氣質(zhì)量狀況的表現(xiàn)和特征也會存在很大的區(qū)別。

因此，本文首提出利用動態(tài)時間規(guī)劃模型和系統(tǒng)聚類相結(jié)合，將我國不同城市空氣質(zhì)量指數(shù)變動趨勢及數(shù)值水平相似性強的地區(qū)聚為一類，即認為類內(nèi)城市間空氣質(zhì)量指數(shù)和其變化趨勢相近，并對我國存在大氣污染聯(lián)動性的城市進行有效分割。一方面為治理方式由城市的“單打獨斗”轉(zhuǎn)變?yōu)閰^(qū)域的“聯(lián)防聯(lián)控”提供區(qū)域劃分的依據(jù)；另一方面，為職能部門制定政策提供相應(yīng)的理論依據(jù)，做到因地制宜。

2 國內(nèi)外文獻綜述

美國地理學家Tobler提出的地理學第一定律，所有的事物在空間上都是具有聯(lián)系的，聯(lián)系程度與距離有關(guān)，距離越近聯(lián)系越緊密。這也是最初的空間自相關(guān)理論的形成，即觀察對象在空間上與相鄰單元的觀察對象存在空間上的聯(lián)系，包括高值-高值聚集，低值-低值聚集。張琪敏等以北京、青島等十座具有代表性的城市作為重點研究對象，利用2004年各城市空氣質(zhì)量日報數(shù)據(jù)研究其空氣污染狀況和污染特征差異，用描述統(tǒng)計的方法研究了大氣污染現(xiàn)狀，最終結(jié)論表明大氣污染輸送網(wǎng)路與地理位置、地形有關(guān)[1]。劉曉紅結(jié)合空間計量經(jīng)濟學模型，對我國城市霧霾污染時空分異、動態(tài)演化與影響機制進行探究后發(fā)現(xiàn)，全國和各區(qū)域霧霾污染時空特征顯著，空間集聚性強霧霾污染存在顯著的空間正自相關(guān)，霧霾污染區(qū)域一體化態(tài)勢明顯[2]。龔鵬鵬通過對北京市空氣污染的空間統(tǒng)計分析，揭示了北京市與周邊各省的空氣污染的自相關(guān)性和北京市空氣污染的自相關(guān)性，相鄰城市之間的空氣質(zhì)量指數(shù)在空間上有明顯的聚集效應(yīng)，北京市受保定市、石家莊市污染比較嚴重[3]。以上研究證明，每座城市會因為其特殊的地理位置，在大氣污染傳送網(wǎng)絡(luò)等多種因素的作用下或多或少的受到周邊城市的影響，在這種影響下，空氣污染問題已經(jīng)從單一城市問題上升到了區(qū)域性問題。

賀冉冉等人基于時間序列模型發(fā)現(xiàn)中國東部城市群存在2個值得重視的高空間自相關(guān)區(qū)域：京津冀和長三角地區(qū)[4]。賀祥和南國衛(wèi)等人基于灰色關(guān)聯(lián)模型發(fā)現(xiàn)江蘇省呈沿海低、內(nèi)陸高，南部高、北部低，陜西省呈“北高南低”的空間分布特征?？追宓热死每諝馕廴局笖?shù)（API）對2001年—2015年我國119個重點城市空氣質(zhì)量進行研究后也證明了這一觀點，他們發(fā)現(xiàn)中國北方地區(qū)多年API日值平均狀況高于南方地區(qū)，東部沿海地區(qū)高于西部內(nèi)陸地區(qū)[5]。劉永偉等人基于GIS技術(shù)探究了全國和珠三角空氣質(zhì)量的時空分布規(guī)律，發(fā)現(xiàn)我國空氣污染具有明顯的季節(jié)變化特征和空間集聚特征[6]。

以上研究表明，空氣質(zhì)量指數(shù)空間聚集性顯著，總體表現(xiàn)出“東重西輕、北重南輕”格局，因而有必要研究空間聚集性的具體表現(xiàn)，根據(jù)空氣污染指數(shù)的種種特征將我國劃分為幾大區(qū)域。

3 研究內(nèi)容

運用動態(tài)時間規(guī)整算法計算任意兩座城市的空氣質(zhì)量指數(shù)時間序列變化趨勢及數(shù)值水平的相似性，以該算法計算結(jié)果作為城市間系統(tǒng)聚類的衡量標準，將空氣質(zhì)量指數(shù)變化相似性強的聚為一類，將我國劃分為幾大區(qū)域，即認為類內(nèi)城市之間區(qū)域流動性對空氣質(zhì)量影響相同。

4 區(qū)域劃分

空氣質(zhì)量指數(shù)（Air Quality Index）是我國空氣質(zhì)量評價標準，將監(jiān)測到的6項污染物濃度用統(tǒng)一的標準呈現(xiàn)，指數(shù)越好代表空氣質(zhì)量越好，并將AQI級別分為六級來表示其嚴重程度。

本文基于城市間大氣污染會相互傳輸?shù)恼J識下，認為AQI趨勢變化相同城市同時地理距離較近的幾個城市大氣污染物相互傳播的可能性較大，將這些城市分為一類。希望通過對AQI的變化趨勢和城市間地理距離聚類研究，為我國大氣污染治理方式由城市的“單打獨斗”轉(zhuǎn)變?yōu)閰^(qū)域的“聯(lián)防聯(lián)控”提供區(qū)域劃分的依據(jù)。

4.1 數(shù)據(jù)來源與研究方法

AQI數(shù)據(jù)是從中國空氣質(zhì)量監(jiān)測中心網(wǎng)站爬取的環(huán)境監(jiān)測總站劃列的169個城市污染物濃度的小時數(shù)據(jù)，統(tǒng)計起止時間為2018年1月1日—2018年12月31日。按照《環(huán)境空氣治理標準》（GB 3095-2012）：日均AQI為24小時AQI檢測值的算術(shù)平均。

經(jīng)緯度數(shù)據(jù)從中國城市經(jīng)緯度查詢網(wǎng)站爬取獲得。

4.1.1 動態(tài)時間規(guī)劃

AQI變化趨勢就在大氣污染傳送網(wǎng)絡(luò)的影響下，可能會出現(xiàn)趨勢非當期相關(guān)，而是表現(xiàn)出滯后的現(xiàn)象。動態(tài)時間規(guī)劃（Dynamic Time Warping，DTW）算法可用于衡量時間序列的相似性，解決變化趨勢非當期相關(guān)的問題。DTW采用滿足一定條件的時間序列輸入樣本和參考的時間序列的對應(yīng)關(guān)系，求解兩者匹配時累計距離最小所對應(yīng)的規(guī)劃函數(shù)。

4.1.2 系統(tǒng)聚類

系統(tǒng)聚類是將距離相近的樣本先聚類一類，樣本距離遠的后聚為一類，如此反復(fù)進行下去，直到每個樣本聚到適合的類別中。假設(shè)有n個樣本，首先將n個樣本各自分為一類；第二步，選取合適的距離度量，并將類間距離最近的兩個樣本合并為一類，其他樣本仍各自一類，此時有n-1個類別；第三步，再將距離最近的兩個樣本合并為一類，此時有n-2個類別；如此反復(fù)進行下去，直到n個類別變成一個大類為止。

在進行系統(tǒng)聚類前，首先要定義類間距離。類間距的初次計算通常采用歐氏距離：

本文在聚類時不僅參考兩城市間空氣質(zhì)量AQI指數(shù)的DTW距離，同時，還將兩城市間的地理距離作為聚類的另一評判標準。地理距離指標的加入使得聚類結(jié)果可行性增強，由于距離的限制，使得聚類所得各類別內(nèi)部的城市間具有大氣污染相互影響可能的概率變大。最終，本文在系統(tǒng)聚類時，將標準化后的任意兩城市間 DTW 距離和地理距離分別按照0.5、0.5的比重進行加權(quán)，將所得數(shù)值作為聚類的依據(jù)進行系統(tǒng)聚類。

4.2 空氣質(zhì)量區(qū)域劃分研究

本文首先計算城市間空氣污染指數(shù)的DTW距離矩陣與城市間的地理距離矩陣，將兩個指標同時作為聚類依據(jù)進行聚類。

表1

根據(jù)聚類的結(jié)果，本文將169個城市分為7個不同的城市群，其中拉薩和烏魯木齊被單列出來，不能很好的代表區(qū)域，故將其剔除。

（1）長江流域地區(qū)，長江流域地區(qū)，降雨量大且風力作用旺盛，有效抑制地面揚塵發(fā)生并有助于大氣污染物的擴散、稀釋和沉降，改善空氣質(zhì)量。

（2）京津冀地區(qū)，值得注意的是，中國環(huán)境監(jiān)測總站將亳州、淮北、宿州、徐州劃分為長江三角區(qū)域，但是根據(jù)本文的聚類結(jié)果，亳州、淮北、宿州和徐州AQI變化趨勢與京津冀地區(qū)更為相近。這四個城市同時也位于秦嶺-淮河線以北區(qū)域，此線是我國重要的地理分界線，此分界線南北城市氣候差異顯著。

（3）川渝地區(qū)，地處丘陵、盆地的相對封閉自然條件，常年多靜風和微風，空氣濕度相對較大，不利于大氣污染物的擴散，且該地區(qū)經(jīng)濟發(fā)展迅速，產(chǎn)業(yè)聚集，人口集中。

（4）珠江三角地區(qū)，三面環(huán)山一面環(huán)海，以及受嶺南獨特的氣候條件影響，形成了相對封閉的區(qū)域。污染程度最輕，城市氣候適宜、地形特殊，部分地區(qū)人類活動相對不強烈，因此它們的大氣環(huán)境優(yōu)良。

4.3 結(jié)論與討論

城市空氣污染表現(xiàn)出“東重西輕，北中南輕”的空間格局區(qū)域城市空氣污染的總程度和分布結(jié)構(gòu)具備明顯的分異特征