□莊燕杰 朱莉莎
改革開放以來,長(zhǎng)三角地區(qū)經(jīng)濟(jì)發(fā)展成就有目共睹,然而伴隨的環(huán)境問題也逐漸引起人們關(guān)注。近年長(zhǎng)三角地區(qū)頻繁出現(xiàn)霧霾和水體污染,表明低頻、散發(fā)的環(huán)境污染在區(qū)域經(jīng)濟(jì)一體化推動(dòng)下,借助廣泛的經(jīng)濟(jì)聯(lián)系已擴(kuò)大為區(qū)域性環(huán)境問題,使得污染的空間聯(lián)系超越純粹地理學(xué)意義上的相鄰而呈現(xiàn)復(fù)雜網(wǎng)絡(luò)結(jié)構(gòu),給地方環(huán)境治理帶來新挑戰(zhàn)。
十八大以來,黨中央和國(guó)務(wù)院先后出臺(tái)《關(guān)于加快推進(jìn)生態(tài)文明建設(shè)的意見》和《生態(tài)文明體制改革總體方案》,做出“加快建立系統(tǒng)完整的生態(tài)文明制度體系,加速推進(jìn)生態(tài)文明建設(shè)”的戰(zhàn)略決策。在此背景下,長(zhǎng)三角地區(qū)先后簽署《長(zhǎng)三角城市環(huán)境保護(hù)合作(合肥)宣言》《長(zhǎng)三角地區(qū)跨界環(huán)境污染事件應(yīng)急聯(lián)動(dòng)工作方案》《長(zhǎng)三角區(qū)域空氣質(zhì)量改善深化治理方案(2017-2020 年)》和《長(zhǎng)三角區(qū)域水污染防治協(xié)作實(shí)施方案(2018-2020 年)》,以共同應(yīng)對(duì)環(huán)境污染的挑戰(zhàn)。基于此,掌握長(zhǎng)三角地區(qū)環(huán)境污染現(xiàn)狀,深化對(duì)長(zhǎng)三角環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)形成及其特征的研究,對(duì)于促進(jìn)長(zhǎng)三角地區(qū)環(huán)境污染協(xié)同治理意義重大。
長(zhǎng)江三角洲是中國(guó)第一大經(jīng)濟(jì)區(qū),是我國(guó)的經(jīng)濟(jì)中心、國(guó)際門戶、全球先進(jìn)制造業(yè)基地。根據(jù)2019 年國(guó)務(wù)院批準(zhǔn)的《長(zhǎng)江三角洲區(qū)域一體化發(fā)展規(guī)劃綱要》,長(zhǎng)江三角洲包括上海市、江蘇省、浙江省、安徽省,區(qū)域面積35.8 萬(wàn)平方公里。本文研究范圍則涵蓋上述所有地市,以此來客觀分析長(zhǎng)三角地區(qū)環(huán)境污染現(xiàn)狀及其時(shí)空關(guān)聯(lián)特征。
由于大氣污染和水體污染流動(dòng)性較強(qiáng),往往能跨越轄區(qū)邊界給鄰近地區(qū)帶來影響,故本文以大氣污染和水體污染作為主要研究?jī)?nèi)容。基于研究對(duì)象的特點(diǎn)和數(shù)據(jù)可得性,選取2003-2016 年《中國(guó)城市統(tǒng)計(jì)年鑒》中的工業(yè)廢水排放量、工業(yè)SO2排放量和工業(yè)粉塵排放量等作為主要指標(biāo)。
在選定指標(biāo)的基礎(chǔ)上,對(duì)其進(jìn)行標(biāo)準(zhǔn)化處理;然后,采用熵權(quán)法確定各項(xiàng)指標(biāo)權(quán)重,從而有效避免因主觀因素而導(dǎo)致真實(shí)環(huán)境污染綜合指數(shù)的偏誤,以客觀反映長(zhǎng)三角地區(qū)環(huán)境污染現(xiàn)狀;最后,基于線性加權(quán)模型對(duì)工業(yè)廢水排放量、工業(yè)SO2排放量、工業(yè)粉塵排放量等進(jìn)行加權(quán)綜合,從而測(cè)算得到各地環(huán)境污染綜合排放量指標(biāo)。

圖1 長(zhǎng)三角地區(qū)污染綜合排放演變趨勢(shì)圖
在各地市環(huán)境污染綜合排放量基礎(chǔ)上,匯總得到江蘇、浙江和安徽的環(huán)境污染綜合排放量,進(jìn)而繪制長(zhǎng)三角地區(qū)污染綜合排放演變趨勢(shì)圖(見圖1)。由圖1 可知:2003-2016年,長(zhǎng)三角地區(qū)環(huán)境污染綜合排放總體呈下降趨勢(shì)。2005 年長(zhǎng)三角地區(qū)污染綜合排放總量達(dá)到峰值,然后逐年縮減至2009 年的低谷;但2010 年開始驟增,可能與金融危機(jī)后各城市積極引入人才、資金,急于發(fā)展經(jīng)濟(jì)及開發(fā)過度有關(guān);爾后又波動(dòng)中下降,至2016 年達(dá)到最小值,此時(shí)期環(huán)境的改善顯然與長(zhǎng)三角環(huán)境協(xié)同治理的推動(dòng)密不可分。
為分析長(zhǎng)三角地區(qū)環(huán)境污染的時(shí)變特征,利用非參數(shù)核密度估計(jì)方法繪制出2003 年、2007 年、2011 年和2016 年的環(huán)境污染綜合排放指標(biāo)的核密度曲線圖(見圖2)。
由圖2 可知:2003-2016 年,環(huán)境污染分布總體變動(dòng)十分明顯。其中,2003-2011 年核密度曲線逐漸向右移動(dòng),說明長(zhǎng)三角地區(qū)各地污染排放在此時(shí)期有所增加,環(huán)境狀況輕微惡化;在核密度曲線右移之時(shí),其峰度也變得更為陡峭,說明此時(shí)段長(zhǎng)三角各地污染排放有進(jìn)一步集中趨勢(shì)。與2011 年比,2016 年的環(huán)境污染綜合排放核密度曲線變動(dòng)更加顯著,不僅整體左移明顯,其峰度也有所下降、核密度曲線變得更為平緩,表明長(zhǎng)三角地區(qū)各地污染排放在此時(shí)段有顯著減少,各地環(huán)境狀況有極大改善。
進(jìn)而采用空間統(tǒng)計(jì)中的全局Moran’s I 系數(shù)分析長(zhǎng)三角地區(qū)環(huán)境污染空間集聚特征。通過對(duì)2003-2016 年長(zhǎng)三角地區(qū)環(huán)境污染排放進(jìn)行全局自相關(guān)分析,所得結(jié)果(見表1)。
表1 結(jié)果顯示:2003-2016 年綜合環(huán)境指標(biāo)的全局Moran's I 系數(shù)均為正值,其數(shù)值在0.2~0.4 間波動(dòng),且全域空間自相關(guān)性系數(shù)均通過了1%顯著性水平的檢驗(yàn)。因此可以認(rèn)為,2003-2016 年長(zhǎng)三角地區(qū)各地市環(huán)境污染排放指標(biāo)具有顯著的正向空間自相關(guān)性。長(zhǎng)三角地區(qū)環(huán)境污染的空間分布具有集群化特性,即環(huán)境污染高值地區(qū)與高值地區(qū)、低值地區(qū)與低值地區(qū) 集聚在一起的空間特征;而全局Moran's I 系數(shù)也表現(xiàn)為增長(zhǎng)趨勢(shì),預(yù)示著長(zhǎng)三角地區(qū)環(huán)境污染集群化趨勢(shì)正在進(jìn)一步加強(qiáng)。

圖2 環(huán)境污染綜合排放的核密度曲線圖
長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)的形成,一方面依賴于地理位置的鄰近,地理相近地區(qū)經(jīng)濟(jì)聯(lián)系越廣泛,污染物更易伴隨經(jīng)濟(jì)活動(dòng)而擴(kuò)散,導(dǎo)致污染跨境溢出;另一方面,污染物也會(huì)隨著水流、大氣環(huán)流和化學(xué)作用等與周邊地區(qū)進(jìn)行轉(zhuǎn)移交換,從而與周邊地區(qū)產(chǎn)生密切聯(lián)系,導(dǎo)致網(wǎng)絡(luò)關(guān)系的形成(劉華軍等,2015;陳蘊(yùn)恬等,2017);再者,環(huán)境規(guī)制的差異能直接影響污染型企業(yè)的選址和生產(chǎn)決策,對(duì)地區(qū)污染的關(guān)聯(lián)關(guān)系形成與發(fā)展影響最直接,也是長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)形成的重要因素。

表1 歷年Moran’s I 系數(shù)值及檢驗(yàn)結(jié)果
在上述邏輯基礎(chǔ)上,將長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)抽象為各地污染排放相互影響關(guān)系構(gòu)成的集合,采用格蘭杰因果檢驗(yàn)來驗(yàn)證各地污染間的因果關(guān)系,以此構(gòu)建長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)。
利用Python 編程,本研究對(duì)41個(gè)地市的時(shí)間序列數(shù)據(jù)進(jìn)行兩兩格蘭杰因果檢驗(yàn),得到1640 個(gè)結(jié)果。若檢驗(yàn)P 值小于0.1,認(rèn)為兩地區(qū)有顯著格蘭杰因果關(guān)系,將其值取為1;若P 值大于0.1,認(rèn)為兩地區(qū)無格蘭杰因果關(guān)系,則取值為0。依此構(gòu)造出格蘭杰因果關(guān)系矩陣G,進(jìn)而利用軟件Ucinet 6.0 構(gòu)建長(zhǎng)三角環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)。
在格蘭杰因果關(guān)系矩陣G 基礎(chǔ)上,借助Ucinet 6.0 繪制長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)圖(見圖3)。由圖可知:長(zhǎng)三角地區(qū)污染溢出較為普遍,各地污染排放都或多或少與其他地區(qū)存在空間關(guān)聯(lián)關(guān)系。
1.長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)整體特征刻畫
基于軟件Ucinet 6.0,測(cè)算得到長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)的網(wǎng)絡(luò)密度為0.165,說明各地市污染的空間關(guān)聯(lián)關(guān)系較緊密。進(jìn)一步地,由于長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)中不可達(dá)點(diǎn)為0,故該網(wǎng)絡(luò)的關(guān)聯(lián)度為1,表明網(wǎng)絡(luò)連通性較好,各地市關(guān)聯(lián)程度很高。同時(shí),該網(wǎng)絡(luò)的等級(jí)度為0.308,說明41 個(gè)地市在污染溢出方面扮演著不同角色。此外,該網(wǎng)絡(luò)的網(wǎng)絡(luò)效率為0.733,說明各地市污染溢出的渠道較多、網(wǎng)絡(luò)的穩(wěn)定性良好。
2.長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)個(gè)體特征刻畫
(1)點(diǎn)度中心度。運(yùn)用Ucinet 6.0 軟件可計(jì)算各地市的點(diǎn)度中心度(見表2)。

圖3 長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)圖
由表2 可知:41 個(gè)地區(qū)點(diǎn)出度均值為6.59,其中徐州、連云港、鹽城、泰州、湖州、紹興、銅陵、阜陽(yáng)、宣城的點(diǎn)出度均在10 以上,說明這些地市容易影響其他地區(qū)的環(huán)境,原因在于點(diǎn)出度較高地區(qū)的工業(yè)發(fā)達(dá)、產(chǎn)生污染排放較多,在環(huán)境規(guī)制、地理環(huán)境、自然因素的作用下,污染排放更易向其他地區(qū)溢出。上海、南京、淮安、鎮(zhèn)江、杭州、淮北、池州的點(diǎn)出度很小,甚至為0,說明這些地區(qū)污染排放對(duì)其他地區(qū)影響很小,幾無污染溢出。
從點(diǎn)入度角度看,點(diǎn)入度均值也為6.59,無錫、杭州和宣城點(diǎn)入度最大,南京、揚(yáng)州、蚌埠、亳州點(diǎn)入度也較高,說明這些地市受其他地區(qū)污染的影響較大;反之,蘇州、南通、寧波、嘉興、溫州、宿州的點(diǎn)入度較小甚至為0,表明這些地市受其他地區(qū)污染的影響較小。
從相對(duì)點(diǎn)度中心度角度看,相對(duì)點(diǎn)度中心度均值為16.46%,鹽城的相對(duì)點(diǎn)度中心度最高,為36.25%,宣城、連云港、衢州、合肥、等地次之,表明這些地區(qū)的污染與其他地區(qū)關(guān)聯(lián)度較高,環(huán)境污染溢出現(xiàn)象顯著存在。
(2)中間中心度。根據(jù)Ucinet 6.0 軟件計(jì)算各地市中間中心度結(jié)果(見表3)。
表3 結(jié)果顯示:41 個(gè)地區(qū)的中間中心度均值為46.78,平均相對(duì)中間中心度為3.00%。鹽城、宣城、泰州、衢州和安慶等13 個(gè)地區(qū)的相對(duì)中間中心度高于均值,表明這些地市在長(zhǎng)三角環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)中對(duì)其余地市的控制能力較強(qiáng)。其中,鹽城相對(duì)中間中心度為14.89%,遠(yuǎn)高于其他地區(qū),宣城、泰州和衢州的相對(duì)中間中心度均高于10%,說明這幾個(gè)地區(qū)也處于網(wǎng)絡(luò)較為核心的位置,在污染跨境溢出中處于主導(dǎo)和中介地位;淮安、鎮(zhèn)江、溫州、淮北和宿州的中間中心度、相對(duì)中間中心度均很低,說明其在網(wǎng)絡(luò)中的承接作用不顯著。

表2 長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)點(diǎn)度中心度分析表
本文研究三角地區(qū)環(huán)境污染的動(dòng)態(tài)演化及其時(shí)空關(guān)聯(lián)特征,結(jié)果發(fā)現(xiàn):
1.總體來看,2003-2016 年長(zhǎng)三角地區(qū)環(huán)境污染綜合排放在總體上呈下降趨勢(shì)。其中,2003-2011 年各地污染排放有所增加且有集中的趨勢(shì),總體環(huán)境狀況有一定惡化,2011-2016 年各地污染排放有顯著減少,環(huán)境污染狀況得到極大改善。
2.2003 -2016 年,全局Moran's I系數(shù)顯著為正,表明長(zhǎng)三角地區(qū)環(huán)境污染排放具有顯著正向空間自相關(guān),揭示出污染排放在地域分布上具有明顯空間集聚特征——即環(huán)境污染高值地區(qū)與高值地區(qū)、低值地區(qū)與低值地區(qū)相集聚的模式;而Moran's I系數(shù)的增長(zhǎng)趨勢(shì),預(yù)示著污染排放的集群化趨勢(shì)將進(jìn)一步加強(qiáng)。
3.長(zhǎng)三角地區(qū)環(huán)境污染具有明顯的溢出效應(yīng),環(huán)境污染現(xiàn)象廣泛存在。連云港、鹽城和衢州等地點(diǎn)出度很高,容易對(duì)其他地市產(chǎn)生影響;上海、南京和杭州等地的點(diǎn)出度較小,對(duì)其他地區(qū)的影響微弱。無錫、杭州、宣城等地的點(diǎn)入度較高,容易受到其他地市污染的影響;而蘇州、南通、寧波、嘉興、溫州和宿州等地點(diǎn)入度較小,不易受其他地市污染的影響。
首先,加強(qiáng)環(huán)境管治,持續(xù)推進(jìn)生態(tài)文明建設(shè)。長(zhǎng)三角地區(qū)環(huán)境污染綜合排放正處下降通道,各地應(yīng)繼續(xù)加強(qiáng)污染排放管控力度,推進(jìn)清潔能源使用、推廣綠色生產(chǎn),以創(chuàng)新驅(qū)動(dòng)降低污染排放。同時(shí)堅(jiān)決貫徹科學(xué)發(fā)展理念,既要金山銀山,也要綠水青山,持續(xù)推動(dòng)生態(tài)文明建設(shè)。

表3 長(zhǎng)三角地區(qū)環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)中間中心度分析表
其次,發(fā)揮區(qū)域聯(lián)動(dòng)作用,推進(jìn)跨區(qū)協(xié)同治理。污染跨境溢出不容忽視,這就要求那些發(fā)展快、帶動(dòng)力強(qiáng)的地區(qū)在解決自身環(huán)境問題的同時(shí)也應(yīng)具備大局觀念,不僅要注重解決自身環(huán)境問題,還要考慮對(duì)周邊地區(qū)的影響。各級(jí)政府通過強(qiáng)化環(huán)境領(lǐng)域協(xié)調(diào)合作,發(fā)揮區(qū)域聯(lián)動(dòng)作用,有利于推進(jìn)跨境污染協(xié)同治理。
最后,環(huán)境治理因地施策,污染減排抓住重點(diǎn)。在環(huán)境污染空間關(guān)聯(lián)網(wǎng)絡(luò)中,徐州、蘇州、連云港等地對(duì)其他地市環(huán)境影響較大,這些地區(qū)應(yīng)提高環(huán)保投入,在減少本地污染基礎(chǔ)上,進(jìn)一步降低對(duì)周邊地區(qū)的污染溢出。上海、南京和杭州等易受其他地市污染影響的地區(qū)可適當(dāng)提高周邊勞動(dòng)密集型產(chǎn)業(yè)準(zhǔn)入門檻,通過清潔技術(shù)創(chuàng)新和擴(kuò)散,降低外來污染影響。鹽城、宣城和泰州在空間關(guān)聯(lián)網(wǎng)絡(luò)中處于核心中介地位,政府和公眾尤應(yīng)重視這些地區(qū)的環(huán)境狀況,在節(jié)能減排的同時(shí),增強(qiáng)民眾環(huán)境監(jiān)督意識(shí),督促當(dāng)?shù)卣訌?qiáng)環(huán)境監(jiān)管力度。