東童童，鄧世成，李軒楓

(1.中共廣東省委黨校 經(jīng)濟(jì)學(xué)教研部，廣東 廣州 510053；2.華東理工大學(xué) 商學(xué)院，上海 200237)

粵港澳大灣區(qū)是我國首個(gè)經(jīng)國家認(rèn)可的灣區(qū)，也是我國開放度最高、經(jīng)濟(jì)活力最強(qiáng)的城市群之一。近年來，粵港澳大灣區(qū)城市間、產(chǎn)業(yè)間協(xié)作更加密切，人口與產(chǎn)業(yè)集聚程度日益加深。與此同時(shí)，日益深化的經(jīng)濟(jì)關(guān)聯(lián)也給大灣區(qū)發(fā)展帶來了經(jīng)濟(jì)負(fù)外部性，最顯著且亟待解決的問題之一便是大氣污染[1]。污染物排放是規(guī)模經(jīng)濟(jì)的附產(chǎn)品之一，生產(chǎn)過程中污染源大量集中，造成污染物通過大氣環(huán)境在區(qū)域內(nèi)各城市間輸送[2]。城市之間的動(dòng)態(tài)關(guān)聯(lián)構(gòu)成了一個(gè)以城市為節(jié)點(diǎn)的復(fù)雜網(wǎng)絡(luò)，加之大氣污染的現(xiàn)實(shí)特性，由此導(dǎo)致污染區(qū)域內(nèi)所有城市都無法獨(dú)善其身，大氣污染問題呈現(xiàn)“一榮俱榮，一損俱損”的特征[3,4]。

社會(huì)網(wǎng)絡(luò)分析法是一種針對(duì)關(guān)系數(shù)據(jù)的跨學(xué)科分析方法，為研究環(huán)境污染問題提供了一種新的分析視角和技術(shù)工具[5]。目前，國內(nèi)外學(xué)者們運(yùn)用社會(huì)網(wǎng)絡(luò)分析法對(duì)大氣污染問題開展了相關(guān)的研究。Diks，C.和Panchenko，V.[6]運(yùn)用非線性Granger因果檢驗(yàn)方法識(shí)別了我國省際霧霾污染與經(jīng)濟(jì)增長(zhǎng)之間的非線性傳導(dǎo)效應(yīng)，并在此基礎(chǔ)上揭示了其聯(lián)動(dòng)網(wǎng)絡(luò)結(jié)構(gòu)特征；劉華軍、劉傳明、孫亞男等[7,8]的研究發(fā)現(xiàn)，空間網(wǎng)絡(luò)結(jié)構(gòu)分析可有效揭示環(huán)境污染所呈現(xiàn)的復(fù)雜的、多線程的空間溢出關(guān)系，并且這種關(guān)系在區(qū)域大氣污染中表現(xiàn)為一種動(dòng)態(tài)關(guān)聯(lián)的、穩(wěn)定性強(qiáng)的多線程復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu)形態(tài)。另有學(xué)者從其他視角對(duì)區(qū)域大氣污染開展了廣泛深入的研究，多數(shù)學(xué)者主要聚焦在大氣污染的影響因素和大氣污染治理兩個(gè)方面。茹少峰、雷振宇[9]認(rèn)為，導(dǎo)致霧霾天氣形成的主要原因是粗放的經(jīng)濟(jì)發(fā)展方式；邵帥、李欣、曹建華等[10]研究發(fā)現(xiàn)，大氣污染與經(jīng)濟(jì)增長(zhǎng)之間存在著顯著的“U”形曲線關(guān)系，對(duì)大氣污染防治的前提在于對(duì)霧霾污染社會(huì)經(jīng)濟(jì)根源的準(zhǔn)確識(shí)別；王立平、陳俊[11]研究發(fā)現(xiàn)，產(chǎn)業(yè)結(jié)構(gòu)、能源消費(fèi)結(jié)構(gòu)、人口規(guī)模等是對(duì)大氣污染空間相關(guān)性和具有抗干擾性的“穩(wěn)健性”影響因素。

目前學(xué)者們從經(jīng)濟(jì)視角對(duì)大氣污染問題開展了廣泛深入的研究，但存在需要改進(jìn)和進(jìn)一步研究的領(lǐng)域。一是部分研究采用傳統(tǒng)空間計(jì)量方法分析大氣污染的空間關(guān)聯(lián)問題，由此導(dǎo)致空間關(guān)聯(lián)局限于經(jīng)濟(jì)地理學(xué)上的“相近”，而區(qū)域范圍的大氣污染存在多線程復(fù)雜網(wǎng)絡(luò)關(guān)系，空間計(jì)量分析方法無法深入揭示這一重要特性。二是多數(shù)學(xué)者對(duì)大氣污染空間關(guān)聯(lián)的研究采用“屬性數(shù)據(jù)”，而這種數(shù)據(jù)分析方法無法刻畫大氣污染的整體空間網(wǎng)絡(luò)關(guān)聯(lián)關(guān)系。為彌補(bǔ)上述研究不足，本文將通過引入修正引力模型，確定城市大氣污染的關(guān)系矩陣與像矩陣進(jìn)而確定空間關(guān)聯(lián)關(guān)系；采用社會(huì)網(wǎng)絡(luò)分析法(SNA)，以城市的關(guān)聯(lián)數(shù)據(jù)為基礎(chǔ)，對(duì)大氣污染整體空間關(guān)聯(lián)特征和內(nèi)部結(jié)構(gòu)關(guān)系效應(yīng)進(jìn)行分析；通過空間關(guān)聯(lián)關(guān)系整體指標(biāo)網(wǎng)絡(luò)密度、網(wǎng)絡(luò)等級(jí)度、網(wǎng)絡(luò)關(guān)聯(lián)度，刻畫大氣污染空間關(guān)聯(lián)網(wǎng)絡(luò)的整體特征及其變化趨勢(shì)；通過個(gè)體指標(biāo)點(diǎn)度中心度、中介中心度、接近中心度，刻畫單個(gè)城市在網(wǎng)絡(luò)中所處地位與作用；通過大氣污染的空間聚類分析與角色劃分，刻畫關(guān)聯(lián)網(wǎng)絡(luò)內(nèi)部結(jié)構(gòu)及其作用關(guān)系。

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

1.1 研究方法

社會(huì)網(wǎng)絡(luò)分析(SNA)是以單個(gè)或多個(gè)社會(huì)行動(dòng)者之間的關(guān)系作為基本量化分析單位，以代數(shù)方法、圖論工具描述網(wǎng)絡(luò)關(guān)系模式，從而研究空間關(guān)系網(wǎng)絡(luò)的整體特征、個(gè)體特征與網(wǎng)絡(luò)結(jié)構(gòu)的跨學(xué)科分析方法，其應(yīng)用范圍已逐漸從社會(huì)學(xué)向經(jīng)濟(jì)學(xué)、管理學(xué)等領(lǐng)域拓展，成為一種新的研究范式[12,13]。

1.2 網(wǎng)絡(luò)確定與數(shù)據(jù)來源

本文將區(qū)域經(jīng)濟(jì)聯(lián)系量化分析的引力模型運(yùn)用到之間粵港澳大灣區(qū)11市大氣污染的空間聯(lián)系上，基于典型的公式進(jìn)行修正，修正后的引力模型為：

(1)

式中，i、j表示粵港澳大灣區(qū)11個(gè)城市；yij表示城市i和城市j之間大氣污染的相互關(guān)系；kij表示城市i在城市j之間大氣污染聯(lián)系中的貢獻(xiàn)率；Pi、Pj表示城市i、j的年末總?cè)丝?；Ei、Ej分別表示城市i、j的大氣污染物排放量，包括PM2.5、PM10、NO2、SO2等4類污染物；Gi、Gj分別表示城市i、j的實(shí)際城市GDP；dij表示城市i與城市j之間的球面距離，gi、gj分別表示城市i、j的人均GDP，用dij與gi-gj之比Dij的平方表示城市之間的“距離”。根據(jù)式(1)，得到Y(jié)=(yij)11×11的引力矩陣，取矩陣中每一列的平均值作為該列的“列閾值”，然后用該列的各值與“列閾值”比較。若大于閾值，記為1，表示該列的城市與該行的城市在大氣污染上具有關(guān)聯(lián)關(guān)系;反之記為0，表示不存在關(guān)聯(lián)關(guān)系。由此構(gòu)建出城市間大氣污染的空間關(guān)聯(lián)網(wǎng)絡(luò)矩陣，由于計(jì)算的引力矩陣X是非對(duì)稱的，所以這是一個(gè)有向空間關(guān)聯(lián)網(wǎng)絡(luò)。

本文以粵港澳大灣區(qū)11個(gè)城市作為網(wǎng)絡(luò)節(jié)點(diǎn)，實(shí)證考察了大灣區(qū)大氣污染的空間關(guān)聯(lián)，樣本的時(shí)期跨度為2013—2016年，大氣污染數(shù)據(jù)(除PM2.5數(shù)據(jù))來源于2014—2017年廣東省《環(huán)境狀況公報(bào)》、香港環(huán)境保護(hù)署、澳門環(huán)境保護(hù)局空氣質(zhì)量在線分析平臺(tái)；PM2.5數(shù)據(jù)來源于巴特爾研究所、哥倫比亞大學(xué)國際地球科學(xué)信息網(wǎng)絡(luò)中心。引力模型計(jì)算數(shù)據(jù)來源于相應(yīng)年份《廣東統(tǒng)計(jì)年鑒》《香港統(tǒng)計(jì)年鑒》和《澳門統(tǒng)計(jì)年鑒》，地區(qū)之間的地理距離以地區(qū)之間的球面距離表示。

2 大氣污染空間關(guān)聯(lián)整體網(wǎng)絡(luò)特征

2.1 整體網(wǎng)絡(luò)特征及演變趨勢(shì)

本文基于修正的引力模型確定粵港澳大灣區(qū)城際大氣污染的空間關(guān)聯(lián)矩陣，并利用UCINET軟件繪制空間關(guān)聯(lián)網(wǎng)絡(luò)圖。由于篇幅限制，本文僅報(bào)告了2013年和2016年粵港澳大灣區(qū)的空間網(wǎng)絡(luò)結(jié)構(gòu)形態(tài)。圖1—4分別展示了2013年和2016年粵港澳大灣區(qū)城際PM2.5、PM10、SO2和NO2空間網(wǎng)絡(luò)結(jié)構(gòu)形態(tài)。從圖中顯示的情況可以看出，粵港澳大灣區(qū)城際大氣染物主要污染物的空間關(guān)聯(lián)均呈現(xiàn)較為典型的網(wǎng)絡(luò)結(jié)構(gòu)特征，同時(shí)發(fā)展趨勢(shì)較為穩(wěn)定，空間關(guān)聯(lián)性趨于緊密。

圖1 2013年、2016年粵港澳大灣區(qū)城際PM2.5空間網(wǎng)絡(luò)結(jié)構(gòu)

圖2 2013年、2016年粵港澳大灣區(qū)城際PM10空間網(wǎng)絡(luò)結(jié)構(gòu)

圖3 2013年、2016年粵港澳大灣區(qū)城際SO2空間網(wǎng)絡(luò)結(jié)構(gòu)

圖4 2013年、2016年粵港澳大灣區(qū)城際NO2空間網(wǎng)絡(luò)結(jié)構(gòu)

2.2 空間關(guān)聯(lián)的網(wǎng)絡(luò)密度關(guān)系

網(wǎng)絡(luò)密度反映了各城市之間大氣污染空間關(guān)聯(lián)網(wǎng)絡(luò)的緊密程度。網(wǎng)路密度越大，則各城市節(jié)點(diǎn)間大氣污染的空間關(guān)聯(lián)關(guān)系越為緊密。大氣污染網(wǎng)絡(luò)密度的計(jì)算公式為：

(2)

2013—2016年粵港澳大灣區(qū)大氣主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)密度趨勢(shì)見圖5。整體來看，各主要污染物空間網(wǎng)絡(luò)密度均保持在0.45水平以上，城際大氣污染存在一定程度上的空間關(guān)聯(lián)性。其中，PM2.5和SO2空間關(guān)聯(lián)網(wǎng)絡(luò)密度波動(dòng)幅度較大，增長(zhǎng)趨勢(shì)明顯。PM2.5空間關(guān)聯(lián)網(wǎng)絡(luò)密度從2013年的0.4909上升到2016年的0.6000，SO2空間關(guān)聯(lián)網(wǎng)絡(luò)密度從2013年的0.5091上升到2016年的0.5273。而PM10和NO2的空間關(guān)聯(lián)網(wǎng)絡(luò)密度均保持在0.4909水平，趨于穩(wěn)定，與之對(duì)應(yīng)的網(wǎng)絡(luò)關(guān)聯(lián)關(guān)系數(shù)也呈現(xiàn)相同的變化趨勢(shì)。

圖5 2013—2016年粵港澳大灣區(qū)大氣主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)密度趨勢(shì)

2.3 空間關(guān)聯(lián)的網(wǎng)絡(luò)關(guān)聯(lián)關(guān)系

網(wǎng)絡(luò)關(guān)聯(lián)度反映了整體大氣污染關(guān)聯(lián)網(wǎng)絡(luò)的穩(wěn)健性強(qiáng)度。大氣的空間網(wǎng)絡(luò)關(guān)聯(lián)度越高，表明城市之間大氣污染存在越強(qiáng)的空間關(guān)聯(lián)關(guān)系。整體網(wǎng)絡(luò)穩(wěn)健性對(duì)城市大氣關(guān)聯(lián)關(guān)系依賴程度越高，大氣污染自身的空間網(wǎng)絡(luò)穩(wěn)健性則越弱。大氣污染網(wǎng)絡(luò)關(guān)聯(lián)度的計(jì)算公式為：

(3)

2013—2016年，粵港澳大灣區(qū)大氣主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)關(guān)聯(lián)關(guān)系數(shù)趨勢(shì)見圖6。

圖6 2013—2016年粵港澳大灣區(qū)大氣主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)關(guān)聯(lián)關(guān)系數(shù)趨勢(shì)

從圖6可見，走勢(shì)與網(wǎng)絡(luò)密度趨于一致，整體在54上下波動(dòng)，關(guān)聯(lián)關(guān)系數(shù)目最大值為66個(gè)(2016年P(guān)M2.5關(guān)系數(shù))與最大可能的空間關(guān)聯(lián)關(guān)系數(shù)目(110)相差較大。分析表明，粵港澳大灣區(qū)大氣污染的城際空間關(guān)聯(lián)強(qiáng)度處于合理水平，網(wǎng)絡(luò)中各城市的聯(lián)系強(qiáng)度逐漸提高，尤其在PM2.5和SO2污染的空間網(wǎng)絡(luò)中，各城市空間關(guān)聯(lián)越來越密切。

2.4 空間關(guān)聯(lián)的網(wǎng)絡(luò)等級(jí)度

網(wǎng)絡(luò)等級(jí)度反映大氣污染城市之間關(guān)聯(lián)關(guān)系的非對(duì)稱可達(dá)性。網(wǎng)絡(luò)等級(jí)越高，表明霧大氣污染關(guān)聯(lián)網(wǎng)絡(luò)中等級(jí)結(jié)構(gòu)層次越分明，越多的城市在大氣污染關(guān)聯(lián)網(wǎng)絡(luò)中處于從屬地位。在城際大氣污染空間網(wǎng)絡(luò)中，網(wǎng)絡(luò)效率越低，說明城市間存在更多的連線，溢出渠道多，網(wǎng)絡(luò)也就更穩(wěn)定。網(wǎng)絡(luò)等級(jí)刻畫了網(wǎng)絡(luò)中各城市非對(duì)稱可達(dá)性，反映了各城市的支配地位和相互影響程度，計(jì)算公式為：

(4)

式中，G為大氣污染網(wǎng)絡(luò)等級(jí)度；S為對(duì)稱可達(dá)城市節(jié)點(diǎn)對(duì)數(shù)；max(S)為關(guān)聯(lián)網(wǎng)絡(luò)中存在大氣污染關(guān)聯(lián)城市節(jié)點(diǎn)總對(duì)數(shù)。

2013—2016年，粵港澳大灣區(qū)大氣主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)等級(jí)度趨勢(shì)見圖7。

圖7 2013—2016年粵港澳大灣區(qū)大氣主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)等級(jí)度趨勢(shì)

整體來看，大灣區(qū)大氣主要污染物的網(wǎng)絡(luò)等級(jí)度波動(dòng)趨勢(shì)明顯，處于下降態(tài)勢(shì)。以2015年為節(jié)點(diǎn)可分為兩個(gè)階段。2015年以前為震蕩下滑期，其中NO2下降幅度最大，網(wǎng)絡(luò)等級(jí)度從2013年0.6116下降到2015年的0.3305，降幅達(dá)46.0%，PM2.5的網(wǎng)絡(luò)等級(jí)度較為平穩(wěn)，穩(wěn)定在0.6154水平；2015年后，PM10與SO2的網(wǎng)絡(luò)等級(jí)度變化趨勢(shì)一致，均表現(xiàn)為增長(zhǎng)趨勢(shì)，2016年網(wǎng)絡(luò)等級(jí)度分別較2015年增長(zhǎng)了10.1%和40.0%，但仍低于2103年的水平，PM2.5與NO2的網(wǎng)絡(luò)等級(jí)度變化趨勢(shì)一致，均表現(xiàn)為下降趨勢(shì)，2016年網(wǎng)絡(luò)等級(jí)度分別較2015年降低了45.2%和22.9%。

粵港澳大灣區(qū)大氣污染空間網(wǎng)絡(luò)具有穩(wěn)定性特征，PM2.5和SO2的空間溢出得到有效控制，而PM10和NO2的空間溢出具有長(zhǎng)期的穩(wěn)定性。大灣區(qū)城際大氣污染存在空間關(guān)聯(lián)性，城際大氣污染的相互影響呈現(xiàn)長(zhǎng)期的穩(wěn)定性。大灣區(qū)各城市在大氣污染的空間網(wǎng)絡(luò)中的聯(lián)系增強(qiáng)，尤其在PM2.5與NO2污染空間網(wǎng)絡(luò)中的聯(lián)系更為緊密。分析結(jié)果表明，區(qū)域大氣污染的防控機(jī)制不再適用，大氣防治效果必然受到空間關(guān)聯(lián)效應(yīng)的制約。同時(shí)，細(xì)微顆粒污染物引起的大氣污染空間效應(yīng)最為顯著，亟需構(gòu)建以細(xì)微顆粒污染物為防治重點(diǎn)的大氣污染聯(lián)防機(jī)制。

3 大氣污染空間關(guān)聯(lián)中心性分析

3.1 空間關(guān)聯(lián)網(wǎng)絡(luò)點(diǎn)度中心度分析

點(diǎn)度中心度用來衡量大氣污染城市在整體空間關(guān)聯(lián)網(wǎng)絡(luò)中所處中心地位程度。某個(gè)城市大氣污染點(diǎn)度中心度越大，表明該城市與其他城市的直接關(guān)聯(lián)關(guān)系越多，越處于關(guān)聯(lián)網(wǎng)絡(luò)的中心地位。點(diǎn)度中心度計(jì)算公式為：

(5)

2016年，粵港澳大灣區(qū)城際大氣主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)點(diǎn)度中心度結(jié)果見表1。

表1 2016年大灣區(qū)大氣主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)點(diǎn)度中心度

從表1可見，各城市大氣主要污染物的點(diǎn)度中心度均值均在70以上。從分項(xiàng)來看，PM2.5的點(diǎn)度中心度均值為87.27，中心度較高的城市包括廣州、珠海、佛山、惠州、中山；PM10的點(diǎn)度中心度均值為76.36，中心度較高的城市包括香港、珠海、深圳、惠州、中山；NO2的點(diǎn)度中心度均值為74.55，中心度較高的城市包括珠海、佛山、江門、澳門；SO2的點(diǎn)度中心度均值為80.00，中心度較高的城市包括深圳、佛山、惠州、中山、江門、肇慶。其中，珠海、佛山、惠州、中山4個(gè)城市在4類污染物的點(diǎn)度中心度均較高，且在PM2.5指標(biāo)上的點(diǎn)度中心度最高，均為100，說明這4個(gè)城市在大氣污染空間關(guān)聯(lián)網(wǎng)絡(luò)中與其他城市的關(guān)系數(shù)較多，尤其是以PM2.5污染為代表的空氣污染表現(xiàn)最為明顯。

從點(diǎn)入度和點(diǎn)出度來看，各城市大氣主要污染物的點(diǎn)入和點(diǎn)出度均值均在5—7之間，廣州、澳門在PM2.5指標(biāo)和PM10指標(biāo)上的點(diǎn)入度大于點(diǎn)出度，在SO2指標(biāo)上的點(diǎn)出度大于點(diǎn)入度，說明廣州和澳門的細(xì)微顆粒污染主要受周邊城市的影響較大，SO2有害氣體污染呈顯著的空間溢出效應(yīng)；深圳在PM2.5指標(biāo)和PM10指標(biāo)上的點(diǎn)出度大于點(diǎn)入度，而在SO2指標(biāo)上的點(diǎn)入度大于點(diǎn)出度，說明深圳的細(xì)微顆粒污染呈明顯的空間溢出效應(yīng)，而SO2有害氣體污染主要受周邊城市的影響較大；珠海在PM10、NO2、SO2等指標(biāo)上的點(diǎn)出度大于點(diǎn)入度，在PM2.5指標(biāo)上的點(diǎn)入度大于點(diǎn)出度，說明珠海的PM10、NO2、SO2污染呈明顯的空間溢出效應(yīng)，而PM2.5污染主要受周邊城市的影響較大；佛山4類污染物的點(diǎn)出度均在4以上，點(diǎn)入度均在5以上，且點(diǎn)出度均小于點(diǎn)入度，在區(qū)域大氣污染中屬于“凈流入”型城市，說明佛山的大氣污染主要受周邊城市的影響，大氣污染的溢出效應(yīng)不顯著；惠州、中山在PM2.5指標(biāo)和SO2指標(biāo)上點(diǎn)入度大于點(diǎn)出度，在PM10指標(biāo)上的點(diǎn)出度大于點(diǎn)入度，說明惠州和中山的PM10污染呈明顯的空間溢出效應(yīng)，而PM2.5、SO2污染主要受周邊城市的影響較大；東莞在PM2.5指標(biāo)和SO2指標(biāo)上點(diǎn)出度大于點(diǎn)入度，在PM10指標(biāo)上的點(diǎn)入度大于點(diǎn)出度，說明東莞的PM2.5、SO2污染呈明顯的空間溢出效應(yīng)，而PM10污染主要受周邊城市的影響較大；江門、肇慶在PM10指標(biāo)和SO2指標(biāo)上點(diǎn)入度大于點(diǎn)出度，而在PM2.5指標(biāo)上的點(diǎn)出度大于點(diǎn)入度，說明江門和肇慶的PM2.5污染呈明顯的空間溢出效應(yīng)，而PM10、SO2污染主要受周邊城市的影響較大；香港在PM10指標(biāo)和SO2指標(biāo)上點(diǎn)出度大于點(diǎn)入度，而在PM2.5指標(biāo)上的點(diǎn)入度大于點(diǎn)出度，說明香港的PM10、SO2污染呈明顯的空間溢出效應(yīng)，而PM2.5污染主要受周邊城市的影響較大。

綜上，粵港澳大灣區(qū)中珠海、佛山、惠州、中山在城際大氣污染空間關(guān)聯(lián)網(wǎng)絡(luò)中處于中心地位，大氣污染與其他城市存在較顯著的空間關(guān)聯(lián)。PM2.5污染溢出城市分別是深圳、東莞、江門、肇慶，PM10污染溢出城市分別是深圳、珠海、惠州、中山、香港，NO2污染溢出城市分別是珠海、澳門，SO2污染溢出城市分別是廣州、珠海、東莞、香港、澳門。

3.2 空間關(guān)聯(lián)網(wǎng)絡(luò)中介中心度分析

中介中心度刻畫了某個(gè)城市在關(guān)聯(lián)網(wǎng)絡(luò)中與所處其他城市之間大氣污染關(guān)聯(lián)的“中介”角色大小。若中介中心度越大，則該城市對(duì)其他城市大氣污染關(guān)聯(lián)的控制能力也越大。中介中心度計(jì)算公式為：

(6)

式中，CABi為大氣污染城市的絕對(duì)中介中心度；CRBi為大氣污染城市相對(duì)中介中心度；gjk(i)為城市i經(jīng)過城市j與城市k的捷徑條數(shù)；gjk為城市j與城市k之間的總捷徑條數(shù)。

2016年粵港澳大灣區(qū)城際大氣主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)中介中心度結(jié)果見表2。

表2 2016年大灣區(qū)大氣污染主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)中介中心度

從表2可見，各個(gè)城市大氣主要污染物的中介中心度均值均在1.40以上。分項(xiàng)來看， PM2.5的中介中心度均值為1.41，中心度較高的城市分別是廣州、珠海、佛山、惠州、中山；PM10的中介中心度均值為2.63，中心度較高的城市分別是深圳、肇慶；NO2的中介中心度均值為2.83，中心度較高的城市分別是江門、肇慶、香港、澳門；SO2的中介中心度均值為2.22，中心度較高的城市分別是東莞、中山。說明在PM2.5污染空間關(guān)聯(lián)網(wǎng)絡(luò)中，廣州、珠海、佛山、惠州、中山等城市處于核心地位并發(fā)揮著“中介”和“橋梁”作用；在PM10污染空間關(guān)聯(lián)網(wǎng)絡(luò)中，深圳和肇慶處于核心地位并發(fā)揮著“中介”作用；在NO2污染空間關(guān)聯(lián)網(wǎng)絡(luò)中，江門、肇慶、香港、澳門等城市處于核心地位并發(fā)揮著“中介”作用；在SO2污染空間關(guān)聯(lián)網(wǎng)絡(luò)中，東莞和中山處于核心地位并發(fā)揮著“中介”作用。

3.3 空間關(guān)聯(lián)網(wǎng)絡(luò)接近中心度分析

接近中心度度量某個(gè)城市多大程度上接近于其他城市產(chǎn)生關(guān)聯(lián)關(guān)系，接近中心度越高，表明該城市與其他城市產(chǎn)生關(guān)聯(lián)的距離非常短，與其他城市存在越多的直接短距離聯(lián)系，該城市越處于網(wǎng)絡(luò)的中心行動(dòng)者地位。接近中心度計(jì)算公式為：

(7)

2016年粵港澳大灣區(qū)城際大氣主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)接近中心度結(jié)果見表3。

表3 2016年大灣區(qū)大氣主要污染物空間關(guān)聯(lián)網(wǎng)絡(luò)接近中心度

從表3可見，各城市大氣主要污染物的接近中心度均值均在80以上。從分項(xiàng)來看，PM2.5的接近中心度均值為89.96，中心度較高的城市分別是廣州、珠海、佛山、惠州、中山；PM10的接近中心度均值為81.21，中心度較高的城市分別是佛山、惠州；NO2的接近中心度均值為80.51，中心度較高的城市分別是惠州、東莞、肇慶、香港；SO2的接近中心度均值為84.05，中心度較高的城市分別是廣州、深圳、惠州、東莞、江門、澳門。其中，惠州在4類污染物的接近中心度均較高，均在90以上，明顯高于其他城市。說明惠州在粵港澳大灣區(qū)城際大氣污染空間關(guān)聯(lián)網(wǎng)絡(luò)中扮演著“中心行動(dòng)者”角色，即在大氣污染空間關(guān)聯(lián)網(wǎng)絡(luò)中能快速與其他城市相關(guān)聯(lián)，對(duì)大氣污染的傳導(dǎo)速度較快。

4 結(jié)論與政策啟示

4.1 結(jié)論

本文選取2013—2016年粵港澳大灣區(qū)11市大氣主要污染物數(shù)據(jù)，基于網(wǎng)絡(luò)結(jié)構(gòu)視角，采用修正的引力模型構(gòu)建地區(qū)間大氣污染的空間關(guān)聯(lián)矩陣，并采用社會(huì)網(wǎng)絡(luò)分析法(SNA)對(duì)大灣區(qū)大氣污染的空間關(guān)聯(lián)性進(jìn)行了實(shí)證研究。結(jié)論如下：①大灣區(qū)城際大氣主要污染物的空間關(guān)聯(lián)呈較典型的網(wǎng)絡(luò)結(jié)構(gòu)特征，城際大氣污染存在明顯的空間關(guān)聯(lián)和空間溢出，尤其在PM2.5和SO2污染的空間網(wǎng)絡(luò)中，各城市空間聯(lián)系越來越密切。②大灣區(qū)大氣污染空間網(wǎng)絡(luò)具有穩(wěn)定性特征，PM2.5和SO2的空間溢出得到有效控制，而PM10和NO2的空間溢出具有長(zhǎng)期穩(wěn)定性。③大灣區(qū)中珠海、佛山、惠州、中山在城際大氣污染空間關(guān)聯(lián)網(wǎng)絡(luò)中處于中心地位，其大氣污染與其他城市存在較為顯著的空間關(guān)聯(lián)。④大灣區(qū)各城市網(wǎng)絡(luò)中心度差異明顯，各城市在大氣污染空間網(wǎng)絡(luò)中存在明顯的非均衡地位，惠州在大氣污染空間關(guān)聯(lián)網(wǎng)絡(luò)中主要扮演著中心行動(dòng)者的角色。

4.2 政策啟示

首先，構(gòu)建粵港澳大灣區(qū)大氣污染聯(lián)防聯(lián)控體系，形成跨區(qū)域治污合力。研究表明，粵港澳大灣區(qū)城際大氣污染存在動(dòng)態(tài)關(guān)聯(lián)關(guān)系且趨于密切，區(qū)域內(nèi)單個(gè)城市的空氣質(zhì)量在大氣污染空間網(wǎng)絡(luò)關(guān)聯(lián)下均不能獨(dú)善其身，聯(lián)防聯(lián)控是區(qū)域大氣污染治理的必由之路，能有效破解“以鄰為壑”的治理困境。廣東已建立珠三角區(qū)域大氣污染防治聯(lián)席會(huì)議機(jī)制，初步建立廣佛肇、深莞惠、珠中江3個(gè)經(jīng)濟(jì)圈聯(lián)防聯(lián)控工作機(jī)制。在此基礎(chǔ)之上，要進(jìn)一步依托粵港合作聯(lián)席會(huì)議制度，建立起粵港澳大珠三角聯(lián)防聯(lián)控工作機(jī)制，將環(huán)保合作由粵港、粵澳、港澳雙邊合作推進(jìn)到粵港澳三邊合作。

其次，制定并推行城市差異化的大氣污染防治政策，實(shí)行大氣污染空間分類治理。研究表明，粵港澳大灣區(qū)各城市在大氣污染空間網(wǎng)絡(luò)結(jié)構(gòu)中的中心地位不盡相同，大氣污染的空間溢出效應(yīng)在城市間存在顯著差異。區(qū)域大氣污染的防治政策制定應(yīng)充分認(rèn)識(shí)各城市在大氣污染網(wǎng)絡(luò)中的關(guān)聯(lián)效應(yīng)，因地制宜地制定城市大氣污染防治政策，實(shí)行大氣污染空間分類管理，降低防治成本，提高防治效率。對(duì)于在城際大氣污染空間關(guān)聯(lián)網(wǎng)絡(luò)中處于中心地位的珠海、佛山、惠州、中山等城市，應(yīng)加快產(chǎn)業(yè)結(jié)構(gòu)調(diào)整升級(jí)進(jìn)程，大力開發(fā)清潔能源，引導(dǎo)發(fā)展新興產(chǎn)業(yè)，提高對(duì)高污染、高能耗產(chǎn)業(yè)的準(zhǔn)入門檻，淘汰升級(jí)落后產(chǎn)能，轉(zhuǎn)變經(jīng)濟(jì)發(fā)展方式，實(shí)現(xiàn)綠色發(fā)展。

第三，協(xié)調(diào)推進(jìn)經(jīng)濟(jì)發(fā)展與環(huán)境保護(hù)，實(shí)現(xiàn)大氣污染的源頭治理。研究表明，粵港澳大灣區(qū)4類大氣主要污染物均呈現(xiàn)明顯的空間關(guān)聯(lián)和空間溢出，其中PM2.5和SO2的空間溢出有所緩解，而PM10和NO2的空間溢出具有長(zhǎng)期穩(wěn)定性。治理大氣污染最根本的途徑仍是轉(zhuǎn)變生產(chǎn)和生活方式，調(diào)整產(chǎn)業(yè)和能源結(jié)構(gòu)，從源頭上減少和控制大氣污染物的產(chǎn)生。一方面，優(yōu)化升級(jí)產(chǎn)業(yè)結(jié)構(gòu)，合理調(diào)整產(chǎn)業(yè)布局，促進(jìn)新興產(chǎn)業(yè)發(fā)展，引導(dǎo)形成合理的產(chǎn)業(yè)布局；另一方面，優(yōu)化能源消費(fèi)結(jié)構(gòu)，大力開發(fā)風(fēng)電、核電、氣電等清潔能源，逐步減少煤炭消費(fèi)量，降低煤炭消費(fèi)比重，從根本上治理大氣污染。