鄭麗琳，李旭輝，戴 煒

（1.安徽財(cái)經(jīng)大學(xué) 管理科學(xué)與工程學(xué)院，安徽 蚌埠 233030；2.安徽省統(tǒng)計(jì)局 貿(mào)易外經(jīng)統(tǒng)計(jì)處，合肥 340100）

0 引言

生態(tài)文明下的綠色發(fā)展，應(yīng)重視經(jīng)濟(jì)增長(zhǎng)與資源管理、污染治理、環(huán)境保護(hù)之間的協(xié)調(diào)雙贏。然而現(xiàn)實(shí)一些地區(qū)通過消耗大量的物質(zhì)資源和犧牲一定的生態(tài)環(huán)境來(lái)促進(jìn)經(jīng)濟(jì)增長(zhǎng)，這種不可持續(xù)的經(jīng)濟(jì)發(fā)展模式使得生態(tài)環(huán)境壓力不斷增大。因此，宏觀經(jīng)濟(jì)綠色增長(zhǎng)研究就不能回避生態(tài)環(huán)境承載及資源稟賦限制的影響。如何對(duì)綠色經(jīng)濟(jì)增長(zhǎng)及生態(tài)環(huán)境壓力進(jìn)行科學(xué)的量化測(cè)度，如何描述其內(nèi)在聯(lián)系，對(duì)實(shí)現(xiàn)社會(huì)經(jīng)濟(jì)規(guī)模與自然資源承載能力的協(xié)調(diào)發(fā)展具有重要意義。

近些年來(lái)受益于中部崛起戰(zhàn)略，安徽省經(jīng)濟(jì)發(fā)展取得了顯著成就，與此同時(shí)，生態(tài)環(huán)境也面臨巨大壓力，如何測(cè)度安徽省各地市綠色GDP規(guī)模及生態(tài)環(huán)境壓力大小，分析二者間相關(guān)關(guān)系，對(duì)其實(shí)現(xiàn)可持續(xù)發(fā)展具有重要現(xiàn)實(shí)意義。因此，本文基于能值分析研究安徽省16個(gè)地級(jí)市綠色GDP和生態(tài)環(huán)境壓力之間空間效應(yīng)，一方面可以在統(tǒng)一的能值范疇下計(jì)算實(shí)證數(shù)據(jù)，另一方面利用空間面板模型分析時(shí)空變化下地區(qū)生態(tài)環(huán)境壓力對(duì)綠色經(jīng)濟(jì)增長(zhǎng)影響效應(yīng)，具有較強(qiáng)的科學(xué)性及可操作性。

1 理論方法與數(shù)據(jù)來(lái)源

1.1 能值

能值是指某種流動(dòng)或貯存的能量中包含的另一種類別能量（太陽(yáng)能）的數(shù)值大小，它反映了各種資源、產(chǎn)品或者勞務(wù)服務(wù)在其形成過程中都直接或間接地源自于太陽(yáng)能。地區(qū)綠色發(fā)展不僅包括經(jīng)濟(jì)系統(tǒng)，同樣包括生態(tài)系統(tǒng)、社會(huì)系統(tǒng)等，Odum等人根據(jù)地球生態(tài)經(jīng)濟(jì)系統(tǒng)特征，研究給出了各類資源產(chǎn)品等對(duì)應(yīng)的太陽(yáng)能值轉(zhuǎn)換率?；诖耍疚膶⒉煌貐^(qū)復(fù)合生態(tài)系統(tǒng)中的紛繁生產(chǎn)過程和諸多產(chǎn)品統(tǒng)一折算為太陽(yáng)能值（即太陽(yáng)能焦耳sej）。能值計(jì)算公式為：

其中，Ei表示某種能量或物質(zhì)（J或g），τi表示能量折算系數(shù)或能值轉(zhuǎn)換率 （sej/J或sej/g），數(shù)值出自O(shè)dum（1996）和藍(lán)盛芳（2001）等，Em為總能值和。

1.2 基于能值分析的綠色GDP核算方法

依據(jù)張虹等（2010）的研究，綠色GDP計(jì)算公式為：

其中，∑A表示系統(tǒng)中各種不可再生自然資源的能值-貨幣價(jià)值總和，∑B表示系統(tǒng)中各種使用、消耗的不可再生能源的能值-貨幣價(jià)值總和，∑C表示系統(tǒng)中污染排放造成的生態(tài)環(huán)境損耗能值-貨幣價(jià)值總和。為了表述能值-貨幣價(jià)值，需計(jì)算式（3）能值/貨幣比率，以便于實(shí)現(xiàn)那些無(wú)法用市場(chǎng)價(jià)值衡量的各類資源及環(huán)境損耗核算，獲得綠色GDP貨幣價(jià)值。

1.3 基于能值分析的生態(tài)環(huán)境壓力指數(shù)計(jì)算方法

生態(tài)足跡是指能夠?yàn)槿祟愄峁┵Y源或消化廢棄物的生物生產(chǎn)性土地面積；生態(tài)承載力是指自然界通過自我維持、自我調(diào)節(jié)，所能支撐最大社會(huì)經(jīng)濟(jì)活動(dòng)下的各種類型土地供給加權(quán)之和；生態(tài)環(huán)境壓力指數(shù)則是指一定范圍內(nèi)單位生態(tài)承載面積上的生態(tài)足跡，這個(gè)指標(biāo)表征地區(qū)環(huán)境所承受壓力大小。國(guó)內(nèi)很多學(xué)者利用能值-生態(tài)足跡模型，克服傳統(tǒng)生態(tài)足跡模型中土地生產(chǎn)能力的局限，忽略數(shù)值上具有爭(zhēng)議的產(chǎn)量因子和均衡因子，轉(zhuǎn)而引入全球能值密度或區(qū)域能值密度等概念，計(jì)算生態(tài)足跡和生態(tài)承載力，再延伸獲得生態(tài)盈余、生態(tài)環(huán)境壓力指數(shù)等。本文依據(jù)楊燦等（2016）、譚德明等（2016）的研究，生態(tài)足跡、生態(tài)承載力、生態(tài)環(huán)境壓力指數(shù)計(jì)算公式如下：

其中，ef為人均生態(tài)足跡，i指代化石能源、林地、耕地、草地、水域、建筑用地，ai則分別對(duì)應(yīng)不同類型土地使用面積，ci表示對(duì)應(yīng)項(xiàng)目的能值大小。化石能源項(xiàng)目包括原煤、焦炭、汽油、原油、柴油、燃料油、天然氣（液態(tài)、氣態(tài)）等能值之和；林地項(xiàng)目包括水果及林產(chǎn)品能值之和；耕地項(xiàng)目包括糧食、蔬菜、茶葉、糖料、油料、蛋類等能值之和；草地項(xiàng)目包括奶類、肉類能值之和；水域項(xiàng)目包括水產(chǎn)品能值；建筑用地項(xiàng)目包括社會(huì)電力能值。ec為人均生態(tài)承載力，e為可再生資源的人均能值。為了使安徽省16個(gè)地級(jí)市之間EF、EC、EPI具有可比性，因此沒有采用可再生資源能值密度/區(qū)域能值密度的方法核算，而是采用全球平均能值密度p0統(tǒng)一計(jì)算。全球平均能值密度是全球總能值與全球面積之比，采用楊燦等（2016）一文數(shù)據(jù)，取值為3.104×1014sej/ha。

1.4 數(shù)據(jù)來(lái)源

考慮到2010年、2015年安徽省行政區(qū)劃有重大調(diào)整，從數(shù)據(jù)可比性、可得性出發(fā)選擇2011—2016年安徽省下轄16個(gè)地級(jí)市為研究對(duì)象，主要數(shù)據(jù)來(lái)源于2012—2017年《安徽省統(tǒng)計(jì)年鑒》及部分地級(jí)市統(tǒng)計(jì)年鑒。

2 基于能值分析的安徽省地級(jí)市綠色GDP核算

依據(jù)能值理論，使用Odum和藍(lán)盛芳等的能量折算系數(shù)和能值轉(zhuǎn)換率，對(duì)安徽省16個(gè)地級(jí)市社會(huì)-經(jīng)濟(jì)-自然復(fù)合生態(tài)系統(tǒng)的總能值數(shù)據(jù)進(jìn)行計(jì)算。同時(shí)歸總自然資源和能源產(chǎn)品的消耗能值，廢棄物排放消耗能值等，以便完成∑A、∑B、∑C等能值貨幣價(jià)值量的測(cè)算。利用式（2），獲得各地綠色GDP數(shù)值，并與傳統(tǒng)GDP數(shù)值進(jìn)行比較，計(jì)算綠色GDP占傳統(tǒng)GDP比值。經(jīng)人口數(shù)調(diào)整后得人均值，具體數(shù)據(jù)見表1。

表1 2011—2016年安徽省地級(jí)市人均綠色GDP與占傳統(tǒng)GDP比值核算

2.1 空間動(dòng)態(tài)分析

利用GEODA軟件對(duì)表1中2016年各地人均綠色GDP及其占傳統(tǒng)GDP比值進(jìn)行四等分定級(jí)（見下頁(yè)圖1），顯示兩個(gè)指標(biāo)在空間上對(duì)應(yīng)關(guān)系不強(qiáng)。圖1（a）顯示地區(qū)人均綠色GDP具有空間集聚的特征，具體表現(xiàn)為中部凸起、南高于北、東高于西。第一陣營(yíng)主要出現(xiàn)在皖中皖南的合肥、馬鞍山、蕪湖、銅陵四個(gè)地市；第二陣營(yíng)則主要出現(xiàn)在皖南的池州、宣城、黃山以及皖北的蚌埠市；第三陣營(yíng)較為分散，主要包括滁州、淮北、宿州、安慶；第四陣營(yíng)主要集中在皖北皖中的阜陽(yáng)、亳州、淮南、六安等地市。圖1（b）顯示安徽省16個(gè)地級(jí)市綠色GDP占傳統(tǒng)GDP比值分布狀況。第一陣營(yíng)比值最低，主要包括淮南、淮北、銅陵、馬鞍山四個(gè)傳統(tǒng)的資源型城市；第二陣營(yíng)集中在皖中皖南的安慶、蕪湖、池州、宣城等經(jīng)濟(jì)文化旅游城市；第三陣營(yíng)集中安徽東北部發(fā)展速度靠前的宿州、蚌埠、滁州以及中部的合肥；第四陣營(yíng)比值最高，包括安徽西北部的亳州、阜陽(yáng)、六安等農(nóng)業(yè)型城市以及皖南的旅游型城市黃山。值得一提的是第一陣營(yíng)與第四陣營(yíng)城市發(fā)展模式有較大差別，第一陣營(yíng)四座資源型城市該比值顯著偏低（表1顯示均低于81%），淮南淮北具有百年煤炭開采及工業(yè)發(fā)展史、銅陵是新中國(guó)第一個(gè)銅工業(yè)基地、馬鞍山市因鋼設(shè)市立市，這四座城市都存在資源枯竭、產(chǎn)能過剩、城市轉(zhuǎn)型等現(xiàn)實(shí)問題，因此考慮不可再生資源消耗及環(huán)境排放下的綠色GDP所占比重勢(shì)必最低。第四陣營(yíng)四座城市均為農(nóng)業(yè)大市和旅游強(qiáng)市，而農(nóng)業(yè)及旅游業(yè)則主要涉及可再生更新資源消耗利用問題，因此綠色GDP所占傳統(tǒng)GDP比重最高。

2.2 時(shí)間動(dòng)態(tài)分析

圖1 2016年安徽省地級(jí)市人均綠色GDP（a）及其占傳統(tǒng)GDP比值（b）空間四分位圖

表1反映了近六年來(lái)安徽各地級(jí)市人均綠色GDP規(guī)模擴(kuò)大、占傳統(tǒng)GDP比值總體降低的特征，其中多數(shù)（11個(gè)）地級(jí)市比值呈現(xiàn)下降趨勢(shì)，比值下降最大地區(qū)為安慶市，按幾何平均數(shù)計(jì)算年均下降1.21個(gè)百分點(diǎn)；少數(shù)地級(jí)（馬鞍山、宿州、銅陵、淮南）比值先降后增，變動(dòng)上升最大地區(qū)為銅陵市，比值年均增長(zhǎng)5.1個(gè)百分點(diǎn)。這類比值上升的地區(qū)多為資源型城市。近年來(lái)通過積極深化供給側(cè)結(jié)構(gòu)性改革，有效降低產(chǎn)能、控制污染，積極開拓城市轉(zhuǎn)型發(fā)展道路，資源型城市節(jié)能減排效果初現(xiàn)，綠色經(jīng)濟(jì)發(fā)展勢(shì)頭良好。值得注意的是淮南、銅陵二市2014—2015年數(shù)值跳躍明顯，主要源于2015年其行政區(qū)劃發(fā)生重大變動(dòng)，轄區(qū)內(nèi)土地面積、產(chǎn)值、人口等自然經(jīng)濟(jì)資源變化引發(fā)了相關(guān)數(shù)值調(diào)整。

綠色GDP占傳統(tǒng)GDP比值動(dòng)態(tài)變化與不斷增加的不可再生資源（特別是化石能源）消耗及污染排放量密切關(guān)聯(lián)。安徽省能源消費(fèi)總量由2011年10251.86萬(wàn)噸標(biāo)準(zhǔn)煤上升至2016年的12694.96萬(wàn)噸標(biāo)準(zhǔn)煤，年均增長(zhǎng)率為3.63%，且地區(qū)消耗發(fā)展不平衡。全省工業(yè)三廢排放量除工業(yè)廢水以外，其他污染物排放呈現(xiàn)上升趨勢(shì)，工業(yè)廢水、廢氣、固廢排放量由2010年的70971萬(wàn)噸、17849億立方米、9158萬(wàn)噸變動(dòng)至2016年的49625萬(wàn)噸、25367億標(biāo)立方米、12653萬(wàn)噸。安徽省各地經(jīng)濟(jì)規(guī)模的提升仍以大規(guī)模的自然資源利用為前提。

3 安徽省生態(tài)環(huán)境壓力指數(shù)測(cè)度

基于能值理論，歸納總結(jié)化石能源、林地、耕地、草地、水域、建筑用地能值總和，利用式（5）計(jì)算生態(tài)足跡指標(biāo)；考慮到安徽省處于中部?jī)?nèi)陸，能值生態(tài)承載力的計(jì)算中只考慮了5種可再生資源，包括太陽(yáng)能、風(fēng)能、雨水化學(xué)能、雨水勢(shì)能和地球旋轉(zhuǎn)能。為避免重復(fù)計(jì)算，同一性質(zhì)的能量一般只選取其中的最大值，故可再生資源主要為雨水化學(xué)能和地球旋轉(zhuǎn)能之和代表，利用式（6）計(jì)算生態(tài)承載力指標(biāo)。二者比值即為生態(tài)環(huán)境壓力指數(shù)，數(shù)值越大表明生態(tài)環(huán)境越惡劣，地區(qū)環(huán)境所承受壓力越大。依據(jù)楊青等（2016）對(duì)生態(tài)環(huán)境壓力指數(shù)的劃分，分為安全（0～1）、較安全（1～10）、輕度不安全（10～18）、中度不安全（18～24）、重度不安全（24～30）、極度不安全（大于30）六個(gè)等級(jí)。

表2 2011—2016年安徽省地級(jí)市生態(tài)環(huán)境壓力指數(shù)

3.1 空間動(dòng)態(tài)分析

利用GEODA軟件對(duì)表2中2016年各地人均生態(tài)足跡ef、人均生態(tài)承載力ec及生態(tài)環(huán)境壓力指數(shù)EPI進(jìn)行四等分定級(jí)。下頁(yè)圖2（a）城市生態(tài)足跡分級(jí)顯示，第一陣營(yíng)淮北、淮南、銅陵、馬鞍山四個(gè)資源型城市生態(tài)足跡顯著偏高，究其原因是化石能源消耗量大，2016年僅化石能源生態(tài)足跡面積就分別為57.54、38.51、20.24、20.46hm2，遠(yuǎn)高于同省其他城市。第二陣營(yíng)包括池州、蕪湖、安慶、滁州，作為皖江城市帶“一軸”、“兩翼”重要城市，具有承接產(chǎn)業(yè)轉(zhuǎn)移的重要地位，經(jīng)濟(jì)的快速發(fā)展及產(chǎn)業(yè)轉(zhuǎn)移的加速推進(jìn)，使得化石能源及建筑用地生態(tài)足跡面積居高不下。第三陣營(yíng)包括合肥、蚌埠、六安、宣城，多集中于皖中地區(qū)，合肥作為省會(huì)城市，具有較強(qiáng)的輻射帶動(dòng)作用，城市的快速擴(kuò)張加速了對(duì)用地需求，主要體現(xiàn)在化石能源、建筑用地、耕地等方面需求較大。第四陣營(yíng)包括阜陽(yáng)、亳州、宿州、黃山等地，集中于皖北及南部，這些城市多為農(nóng)業(yè)大市及旅游強(qiáng)市，對(duì)耕地及草地需求反而高于化石能源等，生態(tài)足跡面積相對(duì)較低。圖2（b）城市生態(tài)承載力分級(jí)顯示，安徽省生態(tài)承載力具有南高于北的特征，第一陣營(yíng)為皖南的安慶、池州、黃山、宣城，從自然條件看地形地貌多為丘陵山地，日照降水等相對(duì)豐盈，因此太陽(yáng)能、風(fēng)能、雨水化學(xué)能等可再生資源能值最大；第四陣營(yíng)為皖北的阜陽(yáng)、亳州、淮北、宿州地區(qū)，多為平原地區(qū)，日照降水等相對(duì)偏低，可再生資源能值最??；第二陣營(yíng)為滁州、馬鞍山、六安、銅陵地區(qū)又要比第三陣營(yíng)皖中的合肥、蕪湖、蚌埠、淮南等地在地形地貌上有更多起伏。圖2（c）城市生態(tài)環(huán)境壓力指數(shù)分級(jí)顯示，第一陣營(yíng)仍為淮南、淮北、銅陵、馬鞍山，其生態(tài)足跡偏高、生態(tài)承載力不足的現(xiàn)狀，勢(shì)必決定這類地區(qū)生態(tài)壓力極大，生態(tài)環(huán)境極不安全；第四陣營(yíng)為池州、黃山、宣城、六安，良好的自然環(huán)境和氣候條件、較少的資源消耗污染排放、寬松的人口狀況使得此類地區(qū)高生態(tài)承載力對(duì)應(yīng)低生態(tài)足跡，生態(tài)壓力較小，生態(tài)環(huán)境較安全；第二第三陣營(yíng)其他城市生態(tài)環(huán)境處于輕度不安全、中度不安全狀態(tài)。

3.2 時(shí)間動(dòng)態(tài)分析

表2顯示，2011—2016年各地區(qū)生態(tài)足跡ef數(shù)值總體呈現(xiàn)上升趨勢(shì)，一方面由于收入水平的提高，人們?cè)黾恿思Z食資源需求和肉類需求，擴(kuò)大了對(duì)耕地及草地等生態(tài)性用地的需求，另一方面城市化率的提高、工業(yè)化進(jìn)程的加快發(fā)展使得建筑用地等生態(tài)用地需求持續(xù)增加。但淮南、銅陵兩地2015年生態(tài)足跡數(shù)值發(fā)生明顯下降，究其原因是2015年其行政區(qū)劃變動(dòng)導(dǎo)致轄區(qū)面積、人口等重要變量發(fā)生較大幅度變化，從而導(dǎo)致生態(tài)足跡計(jì)算結(jié)果差異較大。2011—2016年各地區(qū)生態(tài)承載力ec數(shù)值略有起伏，但差異不大，六年間各地波動(dòng)的日照降雨數(shù)量影響可更新環(huán)境資源能值大小，由此改變生態(tài)用地的供給狀況。比較安徽省各地生態(tài)足跡與生態(tài)承載力大小，發(fā)現(xiàn)ef均大于ec，都出現(xiàn)了生態(tài)赤字，安徽省整體生態(tài)發(fā)展處于不可持續(xù)狀態(tài)。同時(shí)，生態(tài)足跡與生態(tài)承載力之間不平衡變化，使得生態(tài)壓力不斷升級(jí)，生態(tài)系統(tǒng)退化、人地關(guān)系緊張。

依據(jù)楊青等（2016）提及的生態(tài)環(huán)境壓力指數(shù)劃分標(biāo)準(zhǔn)，對(duì)四類地區(qū)時(shí)間變化進(jìn)行分析?？傮w來(lái)看，各類地區(qū)2011—2016年生態(tài)環(huán)境壓力指數(shù)雖有波動(dòng)，但變動(dòng)幅度有限；絕大多數(shù)地區(qū)波動(dòng)式上升，正向最大變化地區(qū)為滁州市。其中較安全地區(qū)（黃山、池州、宣城、六安）年均變化為3.9%、2.1%、2.0%、3.1%；較不安全地區(qū)（安慶、滁州、亳州）年均變化為5.7%、9.4%、1.3%；中度不安全地區(qū)（宿州、蚌埠、合肥、蕪湖、阜陽(yáng)）年均變化為4.6%、2.8%、0.4%、3.8%、2.7%；極度不安全地區(qū)（淮南、淮北、銅陵、馬鞍山）年均變化為-6.6%、2.8%、-11.5%、1.3%?；茨稀~陵生態(tài)壓力指數(shù)在2015下降明顯也是由于行政區(qū)劃調(diào)整造成的，所以并不改安徽省地級(jí)市總體生態(tài)環(huán)境壓力時(shí)間變動(dòng)趨勢(shì)。安徽省多數(shù)地區(qū)資源消耗仍在繼續(xù)增長(zhǎng)，對(duì)土地的索取持續(xù)增加，而土地供給相對(duì)有限，生態(tài)赤字較為嚴(yán)重，生態(tài)環(huán)境壓力持續(xù)增大，且有惡化的變動(dòng)趨勢(shì)。

4 綠色GDP比值與生態(tài)環(huán)境壓力指數(shù)空間效應(yīng)分析

4.1 空間面板數(shù)據(jù)模型

王火根等（2007）總結(jié)空間效應(yīng)包括空間依賴性和空間異質(zhì)性。前者是指某地樣本值與其他地區(qū)觀測(cè)值相關(guān)，后者是指由于空間單位異質(zhì)性產(chǎn)生的空間效應(yīng)在區(qū)域?qū)用嫔系姆蔷恍浴？臻g計(jì)量經(jīng)濟(jì)的三種基本模型包括空間滯后模型（SLM）、空間誤差模型（SEM）和空間德賓模型（SDM），利用誤差項(xiàng)和因變量的滯后項(xiàng)反映空間自相關(guān)性。

其中，y為被解釋變量，X為解釋變量向量，β為系數(shù)，ρ、λ分別為空間自回歸系數(shù)和空間自相關(guān)系數(shù)，ε、μ為隨機(jī)誤差項(xiàng)。W為n×n的空間權(quán)重矩陣，權(quán)數(shù)系數(shù)往往根據(jù)實(shí)際情況決定。采用空間鄰接權(quán)重矩陣，遵循Rook相鄰規(guī)則將兩個(gè)地區(qū)擁有共同邊界則視為相鄰，若地區(qū)i和 j邊界相鄰，W 中的元素Wij的值為1，否則即為0；主對(duì)角線上元素為0。權(quán)重矩陣經(jīng)過行標(biāo)準(zhǔn)化處理，用每個(gè)元素除以所在行元素之和，保證每行元素之和為1。

截面空間模型中常用的莫蘭指數(shù)Moran’s I、拉格朗日誤差檢驗(yàn)LMerr、拉格朗日滯后檢驗(yàn)LMlag、穩(wěn)健性檢驗(yàn)Robust LMerr和Robust LMlag等都是針對(duì)截面回歸模型提出的，不能直接用于面板數(shù)據(jù)模型。本文采用STATA軟件中處理方法，用分塊對(duì)角矩陣C代替空間自相關(guān)統(tǒng)計(jì)量計(jì)算公式中的空間權(quán)重矩陣，由此將上述檢驗(yàn)擴(kuò)展至空間面板數(shù)據(jù)。其中，LMerr檢驗(yàn)和LMlag檢驗(yàn)及其穩(wěn)健形式的空間相關(guān)性檢驗(yàn)，不僅可以用來(lái)檢驗(yàn)空間相關(guān)性，還可以為選擇SLM、SEM和SDM模型提供依據(jù)。常用做法及步驟包括：一是做無(wú)約束的OLS估計(jì)，若二者均不顯著，無(wú)需進(jìn)行空間計(jì)量分析；二是檢驗(yàn)SDM能否簡(jiǎn)化為SLM或SEM，若似然比檢驗(yàn)（LR）不可拒絕空間誤差模型零假設(shè)H0∶γ+ρβ=0 ，則SDM可被簡(jiǎn)化為SEM，若似然比（LR）檢驗(yàn)不可拒絕空間滯后模型的零假設(shè)H0∶γ=0，則SDM可被簡(jiǎn)化為SLM，若兩個(gè)假設(shè)都被拒絕，則SDM模型為最優(yōu)模型；三是判斷SLM或SEM優(yōu)劣，若LMlag檢驗(yàn)和LMerr檢驗(yàn)只有一個(gè)顯著，設(shè)定為與顯著統(tǒng)計(jì)量對(duì)應(yīng)的模型；若兩個(gè)均顯著，比較穩(wěn)健性LMlag檢驗(yàn)和穩(wěn)健性LMerr檢驗(yàn)，選擇顯著（相對(duì)顯著）的統(tǒng)計(jì)量對(duì)應(yīng)的空間計(jì)量模型。

4.2 實(shí)證過程及結(jié)果分析

綠色GDP占傳統(tǒng)GDP的比值反映了某國(guó)或地區(qū)剔除資源耗減成本、環(huán)境退化成本、生態(tài)破壞成本以及污染治理成本后最終產(chǎn)品和勞務(wù)的價(jià)值總和占經(jīng)濟(jì)中所生產(chǎn)創(chuàng)造的全部最終產(chǎn)品和勞務(wù)的價(jià)值之比。生態(tài)環(huán)境壓力反映生物生態(tài)性土地需求供給狀況，較高的生態(tài)環(huán)境壓力指數(shù)表明人類對(duì)生態(tài)的需求超過了自然生態(tài)的供給能力。從理論上說二者存在很強(qiáng)的關(guān)聯(lián)關(guān)系，過度的生態(tài)需求將帶來(lái)資源的過度消耗、生活的過度消費(fèi)、污染的過度排放等，而不足的生態(tài)供給則限制了自然生態(tài)對(duì)污染的自凈能力，降低了生態(tài)系統(tǒng)與環(huán)境交換物質(zhì)能量的能力，從而造成環(huán)境退化、生態(tài)破壞，這些問題勢(shì)必會(huì)降低綠色GDP占傳統(tǒng)GDP比值。從數(shù)值分析上看，圖1（b）四分位圖顯示安徽省地區(qū)綠色GDP占比呈現(xiàn)明顯的空間集聚性，表明各地綠色發(fā)展可能存在空間交互效應(yīng)。因此，利用2011—2016年安徽省16個(gè)地級(jí)市綠色GDP比值與生態(tài)環(huán)境壓力指數(shù)，建立空間計(jì)量模型反映二者間存在空間效應(yīng)。

4.2.1 模型選擇

利用STATA13軟件包含時(shí)間固定效應(yīng)的混合OLS面板模型、空間滯后面板模型（SLM）、空間誤差面板模型（SEM）的估計(jì)結(jié)果，具體見表3。

表3 三種面板模型估計(jì)結(jié)果

從回歸結(jié)果看，混合OLS模型Moran’s I指數(shù)、LMlag檢驗(yàn)和LMerr檢驗(yàn)值最大且統(tǒng)計(jì)顯著，表明需進(jìn)行空間計(jì)量分析；空間滯后面板模型（SLM）、空間誤差面板模型（SEM）主要檢驗(yàn)值數(shù)值大幅下降且變得不顯著，說明空間計(jì)量模型能有效降低空間相關(guān)性。但從SLM、SEM模型的似然檢驗(yàn)（LR）看，其分析結(jié)果與混合OLS回歸沒有顯著差異，表明這兩個(gè)空間面板模型雖能緩解空間相關(guān)性，但并不適用。因此進(jìn)一步建立空間德賓面板模型（SDM），具體見表4。根據(jù)SDM模型估計(jì)結(jié)果討論兩個(gè)似然比LR檢驗(yàn)值，在10%的顯著性水平下，可以拒絕空間滯后模型的零假設(shè) H0∶δ=0，也可以拒絕空間誤差模型的零假設(shè)H0∶δ+ρβ=0，即SDM模型不能簡(jiǎn)化為SLM、SEM模型。此外，SDM模型的Moran’s I指數(shù)、LMlag檢驗(yàn)和LMerr檢驗(yàn)值也變得不顯著了，說明空間德賓面板模型能有效描述安徽省綠色經(jīng)濟(jì)增長(zhǎng)的空間相關(guān)性問題。幾類模型可決系數(shù)均較高。因此，本文選擇包含時(shí)間固定效應(yīng)的SDM模型解釋生態(tài)環(huán)境壓力對(duì)地區(qū)綠色GDP比值時(shí)空格局影響。4.2.2 實(shí)證結(jié)果分析

表4顯示在10%的顯著性水平下，地區(qū)綠色GDP比值滯后項(xiàng)系數(shù)ρ為0.022是統(tǒng)計(jì)顯著的，表明安徽省地級(jí)市綠色GDP比值對(duì)臨近地市綠色經(jīng)濟(jì)增長(zhǎng)產(chǎn)生明顯的正向空間溢出效應(yīng)，即存在正外部性。由于存在空間滯后項(xiàng)，SDM模型的邊際效應(yīng)并不能直接利用解釋變量的邊際系數(shù)表示，還需通過直接效應(yīng)、間接效應(yīng)來(lái)反映。

表4 SDM模型估計(jì)結(jié)果

（1）從直接效應(yīng)看，安徽省地區(qū)生態(tài)環(huán)境壓力顯著限制了地區(qū)綠色GDP比值；當(dāng)生態(tài)環(huán)境壓力指數(shù)增長(zhǎng)1時(shí)，導(dǎo)致綠色GDP比值變化-0.0033。生態(tài)環(huán)境壓力指數(shù)上升，意味著人類對(duì)自然生態(tài)的需求更高，即為了維持經(jīng)濟(jì)快速發(fā)展所需投入的能值持續(xù)增長(zhǎng)造成了不少的環(huán)境負(fù)荷壓力，快速發(fā)展的城鎮(zhèn)化進(jìn)程也向地區(qū)生態(tài)系統(tǒng)的恢復(fù)能力提出了挑戰(zhàn)，出口項(xiàng)目的增長(zhǎng)帶來(lái)了社會(huì)經(jīng)濟(jì)系統(tǒng)能值輸出的起伏變化，這些問題均不利于地區(qū)綠色經(jīng)濟(jì)發(fā)展；反之，當(dāng)生態(tài)環(huán)境壓力指數(shù)下降1時(shí)，對(duì)綠色GDP比值產(chǎn)生顯著的正向影響，生態(tài)環(huán)境壓力指數(shù)下降意味著人類對(duì)生態(tài)的需求與供給配置趨于合理，有利于地區(qū)城市生態(tài)健康及可持續(xù)發(fā)展，使得各地經(jīng)濟(jì)在從外界獲取物質(zhì)和能量、輸出產(chǎn)品及排放廢物的同時(shí)，能在社會(huì)系統(tǒng)內(nèi)部保留部分能量，完成最終產(chǎn)品及服務(wù)的生產(chǎn)及提供，以維持健康綠色發(fā)展。

（2）從間接效應(yīng)看，生態(tài)環(huán)境壓力存在負(fù)的空間擴(kuò)散效應(yīng)，即某地生態(tài)環(huán)境壓力變化對(duì)相鄰地市綠色經(jīng)濟(jì)增長(zhǎng)產(chǎn)生消極影響，當(dāng)本地生態(tài)環(huán)境壓力指數(shù)增長(zhǎng)1時(shí)，導(dǎo)致相鄰地市的綠色GDP比值變化-0.0002。一方面安徽省地理環(huán)境、氣象條件存在相似性，皖北皖中多為平原、皖南皖西多為丘陵山地，而日照北多南少、降水北少南多，使得相鄰地市間的生態(tài)承載力互有關(guān)聯(lián)；另一方面化石能源、農(nóng)林牧漁產(chǎn)品的貿(mào)易流通和污染物的擴(kuò)散性，使得生態(tài)足跡不可能只由單獨(dú)某個(gè)地區(qū)承擔(dān)。

（3）EPI的邊際系數(shù)與直接效應(yīng)之差反映反饋效應(yīng)大小，即某地通過影響相鄰地市綠色增長(zhǎng)再反過來(lái)影響本地的綠色增長(zhǎng)，EPI的反饋效應(yīng)為-0.00029，這種反饋效應(yīng)一方面來(lái)自于解釋變量空間交互項(xiàng)（W*EPI），另一方面來(lái)自于間接效應(yīng)。反饋效應(yīng)與間接效應(yīng)差值十分微弱，大小僅為9.37E-05且不顯著，由此反映生態(tài)環(huán)境壓力空間交互項(xiàng)無(wú)反饋效應(yīng)體現(xiàn)。

5 結(jié)論

本文基于能值分析對(duì)安徽省16個(gè)地級(jí)市2011—2016年間綠色GDP及生態(tài)環(huán)境壓力狀況進(jìn)行測(cè)算分析，利用空間面板數(shù)據(jù)模型確定二者間數(shù)量關(guān)系，結(jié)論如下：

（1）安徽省各地人均綠色GDP總量提高，但占傳統(tǒng)GDP比值多數(shù)城市呈現(xiàn)下降趨勢(shì)，且空間發(fā)展不平衡。人均綠色GDP從規(guī)模上逐年提高，地區(qū)分布具有中部凸起、南高于北、東高于西的特征。就綠色GDP占傳統(tǒng)GDP比值分類看，傳統(tǒng)資源型城市淮南、淮北、銅陵、馬鞍山比值最低；農(nóng)業(yè)及旅游強(qiáng)市亳州、阜陽(yáng)、六安、黃山比值最高，數(shù)值相對(duì)穩(wěn)定。安徽省城市發(fā)展應(yīng)因地制宜，資源型城市應(yīng)積極深化供給側(cè)結(jié)構(gòu)性改革，提質(zhì)升級(jí)發(fā)展，有效降低產(chǎn)能，實(shí)現(xiàn)城市轉(zhuǎn)型發(fā)展；農(nóng)業(yè)大市及旅游強(qiáng)市應(yīng)重視開展綠色農(nóng)業(yè)及生態(tài)旅游，穩(wěn)妥推進(jìn)綠色城鎮(zhèn)化；科教文化城市、創(chuàng)新性城市應(yīng)積極開拓創(chuàng)新思路，引導(dǎo)綠色消費(fèi)，推動(dòng)工業(yè)綠色升級(jí)，加快現(xiàn)代服務(wù)業(yè)發(fā)展。

（2）安徽省總體生態(tài)環(huán)境壓力大，整體處于不可持續(xù)狀態(tài)，除行政區(qū)劃調(diào)整的銅陵、淮南以外，其他地區(qū)壓力上行趨勢(shì)明顯，且空間發(fā)展不平衡。所有地區(qū)在2011—2016年間都擴(kuò)大了生態(tài)性土地的需求，生態(tài)足跡上升；生態(tài)承載力具有南高于北的特征，日照降雨等氣象指標(biāo)逐年變化且南強(qiáng)北弱。所有地區(qū)均出現(xiàn)生態(tài)赤字，特別是資源型城市生態(tài)環(huán)境壓力巨大，安全等級(jí)為極度不安全。鑒于生態(tài)承載力多與可再生資源相關(guān)，因此安徽省各地應(yīng)重視耕地、化石能源、建筑用地等主要影響生態(tài)足跡的變化趨勢(shì)，執(zhí)行耕地保護(hù)紅線、生態(tài)保護(hù)紅線、控制城市開發(fā)邊界等措施，提高各城市生態(tài)承載力；注重循環(huán)發(fā)展，嚴(yán)格控制資源總量消耗，減少單位能耗，加強(qiáng)垃圾回收和資源循環(huán)再生利用等。

（3）安徽省地級(jí)市綠色GDP比值與生態(tài)環(huán)境壓力指數(shù)之間存在著顯著空間效應(yīng)?；诳臻g德賓面板模型估計(jì)結(jié)果，安徽省地級(jí)市綠色GDP比值對(duì)臨近地市綠色經(jīng)濟(jì)增長(zhǎng)產(chǎn)生明顯的正向空間溢出效應(yīng)，即存在正外部性；生態(tài)環(huán)境壓力對(duì)地區(qū)綠色增長(zhǎng)產(chǎn)生直接限制效應(yīng)，當(dāng)生態(tài)環(huán)境壓力指數(shù)增長(zhǎng)1時(shí)，直接影響本地綠色GDP比值下降0.33個(gè)百分點(diǎn)；生態(tài)環(huán)境壓力對(duì)相鄰地市綠色經(jīng)濟(jì)增長(zhǎng)產(chǎn)生間接擴(kuò)散效應(yīng)，當(dāng)本地生態(tài)環(huán)境壓力指數(shù)增長(zhǎng)1時(shí)，導(dǎo)致相鄰地市的綠色GDP比值下降0.02個(gè)百分點(diǎn)；但生態(tài)環(huán)境壓力的空間反饋效應(yīng)不顯著。安徽省地處中部，城市的快速發(fā)展使得對(duì)自然生態(tài)的需求維持在較高水平，此間消耗了大量的物質(zhì)資源、產(chǎn)生了較多的污染排放，生態(tài)環(huán)境壓力持續(xù)高位，直接限制了地區(qū)綠色經(jīng)濟(jì)增長(zhǎng)。而地理環(huán)境、氣象條件復(fù)雜，在空間分布上具有一定相似性，且在化石能源、農(nóng)林牧漁產(chǎn)品的貿(mào)易流通和污染物的擴(kuò)散性，使得生態(tài)環(huán)境壓力指數(shù)間接擴(kuò)散負(fù)向影響臨近地區(qū)綠色增長(zhǎng)。因此，應(yīng)進(jìn)一步轉(zhuǎn)變經(jīng)濟(jì)增長(zhǎng)方式，增加資源利用效率，提升可持續(xù)發(fā)展能力，推進(jìn)安徽省走上生產(chǎn)發(fā)展、生活富裕、生態(tài)良好的綠色發(fā)展道路。