山東省綠色發(fā)展效率時(shí)空演化及影響因素分析

許麗夢(mèng)，郭付友，王慧

(曲阜師范大學(xué)地理與旅游學(xué)院，山東 日照 276826)

0 引言

改革開(kāi)放以來(lái)，中國(guó)經(jīng)濟(jì)發(fā)展取得了舉世矚目的成就，被譽(yù)為“中國(guó)奇跡”，但我國(guó)傳統(tǒng)粗放型發(fā)展模式面臨著不可持續(xù)的現(xiàn)實(shí)危機(jī)，資源約束和生態(tài)頑疾已成為制約我國(guó)經(jīng)濟(jì)社會(huì)健康可持續(xù)發(fā)展的難題，經(jīng)濟(jì)增長(zhǎng)與環(huán)境保護(hù)協(xié)調(diào)統(tǒng)一的綠色發(fā)展模式是突破我國(guó)資源環(huán)境瓶頸制約、緩解經(jīng)濟(jì)增長(zhǎng)和資源環(huán)境之間矛盾的關(guān)鍵所在，因此新形勢(shì)下我國(guó)經(jīng)濟(jì)發(fā)展迫切需要走綠色發(fā)展之路。

綠色發(fā)展是傳統(tǒng)發(fā)展模式基礎(chǔ)上的一種創(chuàng)新，是建立在生態(tài)環(huán)境容量和資源承載力的約束條件下，將環(huán)境保護(hù)作為實(shí)現(xiàn)可持續(xù)發(fā)展重要支柱的一種新型發(fā)展模式。綠色發(fā)展也是一種生態(tài)發(fā)展理念，以實(shí)現(xiàn)經(jīng)濟(jì)、社會(huì)和環(huán)境的可持續(xù)發(fā)展為目標(biāo)[1]。綠色發(fā)展的研究一直是國(guó)內(nèi)外學(xué)者研究的熱點(diǎn)領(lǐng)域。國(guó)外學(xué)者最早開(kāi)始綠色發(fā)展實(shí)證研究，主要運(yùn)用數(shù)據(jù)包絡(luò)分析及計(jì)量經(jīng)濟(jì)學(xué)模型來(lái)研究綠色發(fā)展效率，且大多從GDP、能源消耗和碳排放之間的關(guān)系來(lái)研究，其中Ramanathan采用DEA模型對(duì)發(fā)展效率進(jìn)行實(shí)證分析，以GDP和碳排放量作為產(chǎn)出指標(biāo)，非化石能源的消耗作為投入指標(biāo)，驗(yàn)證了能源消費(fèi)與碳排放關(guān)系[2]。Puliafito等[3]分析了人口結(jié)構(gòu)、GDP和能源消耗的相關(guān)關(guān)系，采用碳排放模型和洛特卡-沃爾泰拉模型，結(jié)果表明人口結(jié)構(gòu)與綠色發(fā)展息息相關(guān)。Soytas[4]基于VAR模型研究了能源消耗和碳排放量與GDP之間的關(guān)系，結(jié)果證明能源消耗是碳排放的原因。國(guó)內(nèi)關(guān)于綠色發(fā)展的實(shí)證研究日漸增多，僅從綠色發(fā)展效率的研究來(lái)看，從研究方法選擇上，主要采用DEA模型及其衍生模型，王軍[5]、穆學(xué)英[6]等采用DEA模型對(duì)省際綠色發(fā)展效率進(jìn)行評(píng)價(jià)，聶玉立[7]、張文博[8]、趙領(lǐng)娣[9]、張澤義[10]等采用SBM模型測(cè)度省際或城市綠色發(fā)展效率；從研究?jī)?nèi)容上，部分學(xué)者探究綠色發(fā)展效率時(shí)空格局或探討其影響機(jī)制，楊志江[11]、趙林[12]、林曉[13]、車?yán)赱14]、劉楊[15]等從時(shí)空角度分析其研究區(qū)綠色發(fā)展特征，錢爭(zhēng)鳴[16]、王曉云[17]、鐘水映[18]、葉仁道[19]等探究其研究區(qū)域綠色發(fā)展效率的影響因素。

縱觀國(guó)內(nèi)外研究發(fā)現(xiàn)，雖然關(guān)于綠色發(fā)展效率的研究成果日漸豐富，但研究方法的選取多基于傳統(tǒng)DEA、隨機(jī)前沿函數(shù)以及曼奎斯特指數(shù)，無(wú)法考慮投入產(chǎn)出的松弛性問(wèn)題，亟需尋求新的研究方法。研究區(qū)域多集中在省域?qū)用婊驏|部沿海發(fā)達(dá)地區(qū)，對(duì)山東省綠色發(fā)展效率的研究不多。

進(jìn)入21世紀(jì)以來(lái)，山東省經(jīng)濟(jì)穩(wěn)步發(fā)展，全省生產(chǎn)總值由2001年的9438.3億元增長(zhǎng)至2017年的72678.18億元，年均增長(zhǎng)率達(dá)11.4%，經(jīng)濟(jì)總量在全國(guó)處于領(lǐng)先地位，但長(zhǎng)期以來(lái)高消耗、高污染的經(jīng)濟(jì)增長(zhǎng)方式導(dǎo)致資源環(huán)境承載力接近飽和[20-21]，因此必須改變?cè)械拇址判徒?jīng)濟(jì)增長(zhǎng)方式，推進(jìn)節(jié)能降耗和生態(tài)建設(shè)，促進(jìn)資源的集約使用和生態(tài)環(huán)境的動(dòng)態(tài)平衡，走綠色發(fā)展之路。在山東省全面推進(jìn)綠色發(fā)展的背景下，綠色發(fā)展的質(zhì)量效率成為人們關(guān)注的熱點(diǎn)。該文構(gòu)建綠色發(fā)展效率評(píng)價(jià)指標(biāo)體系，運(yùn)用基于非期望產(chǎn)出的SBM模型測(cè)算2007—2017年山東省17地市綠色發(fā)展效率，分析其時(shí)空演化特征及影響因素，力圖探求山東省綠色發(fā)展時(shí)空演化規(guī)律，并為山東省綠色發(fā)展提升提供參考借鑒。

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

1.1 研究方法

對(duì)于綠色發(fā)展效率評(píng)價(jià)，多數(shù)學(xué)者采用數(shù)據(jù)包絡(luò)分析模型或其衍生模型。數(shù)據(jù)包絡(luò)分析DEA(Data Envelopment Analysis)是評(píng)價(jià)一類具有多輸入、多輸出單元效率的有效方法，傳統(tǒng)的DEA模型存在徑向性等多重嚴(yán)格假定，效率評(píng)價(jià)結(jié)果存在偏差，未將非期望產(chǎn)出納入考慮范圍。在DEA模型的基礎(chǔ)上，2000年Tone Kaoru提出了基于松弛變量的DEA模型，即SBM(Slack Based Measure)模型，該模型有效解決了投入產(chǎn)出的松弛性問(wèn)題，并兼顧非期望產(chǎn)出的影響。因此該文選取基于非期望產(chǎn)出的SBM模型測(cè)度山東省綠色發(fā)展效率[22-25]。

(1)

式中：ρ為綠色發(fā)展效率；m，s1，s2分別為投入、期望和非期望產(chǎn)出因素的個(gè)數(shù)；s-，sg，sb表示投入、期望和非期望產(chǎn)出的松弛量；x,yg,yb分別為投入、期望和非期望產(chǎn)出值；λ是權(quán)重向量；X,Yg,Yb分別為投入、期望和非期望產(chǎn)出構(gòu)成的矩陣。

1.2 數(shù)據(jù)選取

在借鑒前人經(jīng)驗(yàn)的基礎(chǔ)上[26-28]，綜合考慮山東省發(fā)展現(xiàn)狀，建立綠色發(fā)展效率評(píng)價(jià)指標(biāo)體系(表1)。投入指標(biāo)從資本、勞動(dòng)、資源、科技方面來(lái)選擇，由于DEA模型測(cè)算相對(duì)效率，該文選取固定資產(chǎn)投資額來(lái)衡量資本投入，勞動(dòng)投入選取年末就業(yè)人口來(lái)衡量，資源投入以全社會(huì)用電量來(lái)測(cè)度，科技投入以R&D支出來(lái)衡量；期望產(chǎn)出選用地區(qū)生產(chǎn)總值(GDP)及社會(huì)消費(fèi)品零售額來(lái)衡量；非期望產(chǎn)出多為在經(jīng)濟(jì)發(fā)展過(guò)程中產(chǎn)生的環(huán)境問(wèn)題，選取綜合污染物來(lái)衡量，綜合考慮各種污染物的排放量及危害程度，該文將廢水排放量、化學(xué)需氧量排放量、氨氮排放量、二氧化硫排放量、煙塵排放量5個(gè)指標(biāo)作為綜合污染物的子指標(biāo)，采用熵權(quán)法[29]可得各子指標(biāo)的權(quán)重分別為0.14，0.28，0.22，0.12，0.24。

表1 綠色發(fā)展效率評(píng)價(jià)指標(biāo)體系

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

選取山東省17地市作為研究區(qū)域，鑒于數(shù)據(jù)的統(tǒng)計(jì)口徑及可獲得性，仍采用2018年以前的行政區(qū)劃方式，將萊蕪市作為地級(jí)市進(jìn)行分析。研究基礎(chǔ)統(tǒng)計(jì)數(shù)據(jù)主要源自2008—2018年《山東統(tǒng)計(jì)年鑒》及山東省各地市統(tǒng)計(jì)年鑒，部分?jǐn)?shù)據(jù)通過(guò)計(jì)算獲得。

2 山東省綠色發(fā)展效率時(shí)空演變

2.1 山東省綠色發(fā)展效率的時(shí)序演變

2.1.1 省域總體演化特征

根據(jù)綠色發(fā)展效率評(píng)價(jià)指標(biāo)體系，借助于DEA-Solver Pro5軟件，可得2007—2017年基于SBM模型的考慮非期望產(chǎn)出與不考慮非期望產(chǎn)出的山東省17地市綠色發(fā)展效率值(表2)。

表2 考慮(不考慮)非期望產(chǎn)出的綠色發(fā)展效率值

注：表中僅列出部分年份的結(jié)果；括號(hào)內(nèi)表示不考慮非期望產(chǎn)出的綠色發(fā)展效率值。

參考相關(guān)研究成果[12-13]，將綠色發(fā)展效率值ρ劃分為4類：0.9≤ρ≤1為效率最優(yōu)；0.8≤ρ<0.9為效率良好；0.6≤ρ<0.8為效率中等；ρ<0.6為效率無(wú)效。從山東省綠色發(fā)展效率均值來(lái)看，2010年考慮與不考慮非期望產(chǎn)出的綠色發(fā)展效率分別為0.691，0.762，2014年兩值分別為0.707，0.755?？紤]非期望產(chǎn)出的綠色發(fā)展效率要低于不考慮非期望產(chǎn)出的效率值，說(shuō)明不考慮非期望產(chǎn)出導(dǎo)致高估其綠色發(fā)展效率，而考慮非期望產(chǎn)出的綠色發(fā)展效率更加貼合實(shí)際，因此該文選取考慮非期望產(chǎn)出的效率值進(jìn)行分析。由考慮非期望產(chǎn)出的山東省綠色發(fā)展效率均值可知，研究期限內(nèi)山東省總體綠色發(fā)展效率呈下降趨勢(shì)，綠色發(fā)展效率均值由2007年的0.86降至2017年的0.76。盡管山東省強(qiáng)調(diào)推進(jìn)經(jīng)濟(jì)增長(zhǎng)方式由粗放型向集約型轉(zhuǎn)變，大力促進(jìn)“資源節(jié)約型環(huán)境友好型”國(guó)民經(jīng)濟(jì)體系建設(shè)，但發(fā)展新興產(chǎn)業(yè)，形成新的發(fā)展模式，加快新舊動(dòng)力轉(zhuǎn)換需要一定過(guò)程，傳統(tǒng)的經(jīng)濟(jì)發(fā)展方式帶來(lái)的問(wèn)題開(kāi)始顯現(xiàn)，傳統(tǒng)優(yōu)勢(shì)正在減弱，環(huán)境問(wèn)題加重。

2.1.2 各地市時(shí)序演化特征

從變化趨勢(shì)來(lái)看，2007—2017年各地市綠色發(fā)展總體趨勢(shì)良好，部分地區(qū)效果不理想(圖1)。

圖1 17地市綠色發(fā)展效率分類比例變化

研究期限內(nèi)山東省17地市效率最優(yōu)占比先減后增，最終穩(wěn)定在0.3左右；效率良好，占比波動(dòng)較小，最終穩(wěn)定在0.1左右；效率中等占比呈增長(zhǎng)趨勢(shì)，最終穩(wěn)定在0.4左右；效率無(wú)效占比逐漸增大，最終穩(wěn)定在0.2左右，雖然山東省綠色發(fā)展效率良好及以上地市數(shù)量逐漸增加，但綠色發(fā)展效率中等及無(wú)效區(qū)域仍占60%，說(shuō)明山東省仍有多個(gè)城市綠色發(fā)展效率不理想，效率提升之路任重而道遠(yuǎn)。

從時(shí)序變化來(lái)看，研究期限內(nèi)山東省各地市綠色發(fā)展效率隨著時(shí)間波動(dòng)而變化(圖2)，各地市綠色發(fā)展效率變化趨勢(shì)呈現(xiàn)地域性差異，存在部分地市綠色發(fā)展效率上升與部分地市下降并存的趨勢(shì)特征，可將其劃分為增長(zhǎng)型、穩(wěn)定型及衰減型3種類型。

(1)增長(zhǎng)型：省會(huì)城市及沿海較發(fā)達(dá)地市，包括濟(jì)南市、青島市及煙臺(tái)市，濟(jì)南市由效率無(wú)效增長(zhǎng)至效率最優(yōu)，青島及煙臺(tái)兩市均由效率中等提升至效率最優(yōu)。

圖2 2007—2017年山東省各地市綠色發(fā)展效率變化曲線圖

(2)穩(wěn)定型：沿海地區(qū)及部分中西部地區(qū)，包括威海市、萊蕪市、聊城市，其中威海市穩(wěn)定在效率最優(yōu)，萊蕪市穩(wěn)定在效率中等，聊城市穩(wěn)定在效率無(wú)效。

(3)衰減型：中西部地區(qū)以及鋼鐵、煤炭等重工業(yè)占據(jù)主導(dǎo)的城市，包含淄博市、棗莊市、東營(yíng)市、濰坊市、濟(jì)寧市、泰安市、日照市、臨沂市、德州市、濱州市及菏澤市，其中淄博、棗莊、日照、德州、臨沂五市由效率最優(yōu)降至效率中等，濟(jì)寧、泰安兩市由效率良好降至效率中等，濰坊及濱州兩市由效率中等降至效率無(wú)效，東營(yíng)市及菏澤市由效率最優(yōu)降至效率良好，并呈現(xiàn)先減后增的變化趨勢(shì)。

總體而言，沿海地區(qū)及省會(huì)等經(jīng)濟(jì)較發(fā)達(dá)的地區(qū)，其綠色發(fā)展呈現(xiàn)穩(wěn)定或提升趨勢(shì)，其中位于山東半島藍(lán)色經(jīng)濟(jì)區(qū)的青島市，經(jīng)濟(jì)轉(zhuǎn)型較早，海洋產(chǎn)業(yè)及高新技術(shù)產(chǎn)業(yè)繁榮發(fā)展，故其綠色發(fā)展效率呈增長(zhǎng)趨勢(shì)；中西部大部分地區(qū)經(jīng)濟(jì)發(fā)展起步晚，產(chǎn)業(yè)轉(zhuǎn)型慢，承接優(yōu)先發(fā)展地區(qū)的夕陽(yáng)產(chǎn)業(yè)的轉(zhuǎn)移，其中煤炭城市棗莊及鋼鐵產(chǎn)業(yè)逐步興起的日照等市，因重工業(yè)發(fā)展導(dǎo)致其效率呈下降趨勢(shì)；而一些城市因其資源稟賦條件或歷史原因，綠色發(fā)展效率呈現(xiàn)出不同的變化趨勢(shì)，其中依托于石油資源而興起的東營(yíng)市，在發(fā)展石油工業(yè)的同時(shí)積極探索可持續(xù)的發(fā)展方式，堅(jiān)持資源高效利用和生態(tài)環(huán)境改善，故其效率處于呈先減后增的變化趨勢(shì)，而萊蕪市以鋼鐵產(chǎn)業(yè)而興起，傳統(tǒng)產(chǎn)業(yè)轉(zhuǎn)型升級(jí)力度小，新興產(chǎn)業(yè)發(fā)展速度慢，故其綠色發(fā)展處于低效率穩(wěn)定。

2.2 山東省綠色發(fā)展效率的空間演變

2.2.1 區(qū)域綠色發(fā)展差異

標(biāo)準(zhǔn)差反映區(qū)域的絕對(duì)差異，變異系數(shù)反映區(qū)域的相對(duì)差異，計(jì)算可得2007—2017年山東省綠色發(fā)展效率的標(biāo)準(zhǔn)差和變異系數(shù)(圖3)。由標(biāo)準(zhǔn)差分析可知，山東省綠色發(fā)展效率的絕對(duì)差異呈遞增的變化趨勢(shì)，但波動(dòng)幅度較小，由2007年的0.16升至2017年的0.20，區(qū)域絕對(duì)差異逐漸擴(kuò)大。由變異系數(shù)分析可知，變異系數(shù)同樣呈現(xiàn)遞增的變化趨勢(shì)，由2007年的0.18升至2017年的0.26，區(qū)域相對(duì)差異逐漸擴(kuò)大，但變化幅度較小。通過(guò)以上分析可知，山東省綠色發(fā)展效率差異隨時(shí)間變化呈遞增趨勢(shì)，但變化幅度較小，原因在于隨著生態(tài)文明等政策的實(shí)施，各地同時(shí)高速發(fā)展，粗放型經(jīng)濟(jì)發(fā)展方式逐步改變，但長(zhǎng)期的差異積累，加之部分地區(qū)原有經(jīng)濟(jì)發(fā)展方式并未完全改變，致使區(qū)域差異呈現(xiàn)出變化幅度較小的遞增趨勢(shì)。

圖3 2007—2017年山東省綠色發(fā)展效率標(biāo)準(zhǔn)差和變異系數(shù)

2.2.2 綠色發(fā)展效率空間格局演化

借助于Geoda軟件可得2007—2017年山東省各地市綠色發(fā)展效率的全局Moran'sI指數(shù)，結(jié)果發(fā)現(xiàn)研究期限內(nèi)山東省17地市綠色發(fā)展的全局Moran'sI指數(shù)為負(fù)值，且不顯著，說(shuō)明山東省17地市綠色發(fā)展效率無(wú)明顯的趨異性，綠色發(fā)展效率高的城市并無(wú)顯著趨近于綠色發(fā)展效率低的區(qū)域。

為了能夠直觀的顯示山東省綠色發(fā)展效率的空間演化，選取2007年、2012年及2017年截面數(shù)據(jù)，借助于ArcGIS軟件繪制山東省綠色發(fā)展效率空間格局圖(圖4)，借鑒前人研究成果[30],按其地理方位，將山東省劃分為魯東地區(qū)(1)魯東地區(qū)指位于山東省東部沿海的地區(qū)，包括煙臺(tái)市、威海市、青島市、日照市。、魯西北地區(qū)(2)魯西北地區(qū)指山東省西北部的地區(qū)，包括濟(jì)南市、德州市、聊城市、濱州市以及東營(yíng)市。、魯中地區(qū)(3)魯中地區(qū)指位于山東省中部的地區(qū)，包括濰坊市、泰安市、萊蕪市、淄博市、臨沂市。、魯西南地區(qū)(4)魯西南地區(qū)指山東省西南部的濟(jì)寧市、菏澤市、棗莊市。。

(1)從空間演化上來(lái)看，研究期限內(nèi)山東省綠色發(fā)展效率存在顯著的地域性差異，其演變趨勢(shì)大體由西高東低演化為東高西低，高效率區(qū)向東發(fā)展，低效率區(qū)轉(zhuǎn)向中西部。

(2)2007年山東省綠色發(fā)展效率自西向東大體呈“高—低—高—低—高”的分布格局，效率最優(yōu)區(qū)域零散分布，效率良好區(qū)域分布于魯東地區(qū)，效率中等區(qū)域主要分布于魯西南及魯中地區(qū)，效率無(wú)效區(qū)域主要位于魯西北地區(qū)。

(3)2012年綠色發(fā)展效率呈現(xiàn)低效率均衡，效率最優(yōu)區(qū)域及效率良好區(qū)域主要包括魯東、魯西南及魯西北的部分地區(qū)，效率中等區(qū)域主要位于魯中地區(qū)，效率無(wú)效區(qū)域主要分布于魯西北地區(qū)。

(4)2017年山東省綠色發(fā)展效率呈中間低兩側(cè)高的分布格局，效率最優(yōu)區(qū)域主要位于魯東地區(qū)，效率良好區(qū)域分布于魯中地區(qū)，效率無(wú)效及效率中等區(qū)域主要位于魯西北及魯中地區(qū)。

圖4 山東省綠色經(jīng)濟(jì)發(fā)展時(shí)空格局

2.2.3 山東省綠色發(fā)展效率呈現(xiàn)上述空間演化特征的原因

(1)山東省經(jīng)濟(jì)發(fā)展的區(qū)域差異導(dǎo)致其綠色發(fā)展的實(shí)施效果參差不齊，東高西低的經(jīng)濟(jì)發(fā)展格局與綠色發(fā)展高效率區(qū)向東部轉(zhuǎn)移、低效率區(qū)向中西部轉(zhuǎn)移的演化格局密切相關(guān)。沿海地市區(qū)位條件好，經(jīng)濟(jì)發(fā)展水平高，基礎(chǔ)設(shè)施完備，發(fā)展方式轉(zhuǎn)型較早，綠色發(fā)展效率高；中西部?jī)?nèi)陸早期主要依托于資源稟賦發(fā)展工農(nóng)業(yè)，投入少且污染少，但隨著傳統(tǒng)工業(yè)興起，環(huán)境污染加重，綠色發(fā)展相對(duì)滯后。

(2)隨著山東省產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整，各地市三產(chǎn)比重不同，帶動(dòng)綠色發(fā)展效率呈不同空間格局。魯東地區(qū)第三產(chǎn)業(yè)比重顯著提升，魯中及魯西南地區(qū)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整后第二產(chǎn)業(yè)比重仍高于全省均值，魯西北地區(qū)大部分地區(qū)第三產(chǎn)業(yè)占比小于全省均值，其中濟(jì)南市三產(chǎn)比重由2007年6∶45∶49優(yōu)化為2017年4∶36∶60，其綠色發(fā)展呈現(xiàn)明顯的穩(wěn)定提升趨勢(shì)，菏澤市三產(chǎn)占比由27∶48∶25調(diào)整為10∶52∶38，初期低投入低產(chǎn)出，效率高，后期重化工業(yè)發(fā)展，綠色發(fā)展效率降低。

(3)隨著國(guó)家級(jí)發(fā)展戰(zhàn)略的實(shí)施，帶動(dòng)相關(guān)地區(qū)經(jīng)濟(jì)高質(zhì)量發(fā)展。山東半島藍(lán)色經(jīng)濟(jì)區(qū)發(fā)展戰(zhàn)略的實(shí)施，加快了海洋產(chǎn)業(yè)的發(fā)展，帶動(dòng)青島、煙臺(tái)、威海等市的現(xiàn)代海洋經(jīng)濟(jì)的發(fā)展，經(jīng)濟(jì)發(fā)展質(zhì)量提升；黃河三角洲高效生態(tài)經(jīng)濟(jì)區(qū)發(fā)展戰(zhàn)略推動(dòng)高效生態(tài)產(chǎn)業(yè)的發(fā)展，東營(yíng)市等加大經(jīng)濟(jì)結(jié)構(gòu)調(diào)整力度，形成生態(tài)產(chǎn)業(yè)，綠色發(fā)展效率顯著改善。

3 山東省綠色發(fā)展效率的影響因素

3.1 因素選擇與模型設(shè)定

區(qū)域綠色發(fā)展效率差異的形成是多種因素共同作用的結(jié)果，且綠色發(fā)展效率值是[0,1]之間的受限因變量，因此可用面板Tobit模型來(lái)分析綠色發(fā)展效率的影響因素。Tobit模型又稱受限因變量模型，是因變量滿足某種約束條件下取值的模型，使用Tobit模型同時(shí)可以有效避免傳統(tǒng)的普通最小二乘回歸與離散模型在估計(jì)結(jié)果上的偏差[31]。借鑒前人研究成果[16-19]，結(jié)合研究區(qū)現(xiàn)狀，該文選取經(jīng)濟(jì)發(fā)展水平、城鎮(zhèn)化水平、產(chǎn)業(yè)結(jié)構(gòu)、對(duì)外開(kāi)放程度、科技創(chuàng)新等5個(gè)因素進(jìn)行定量分析。

(1)經(jīng)濟(jì)發(fā)展水平(x1)以人均GDP取對(duì)數(shù)(lnGDP)表示，并將人均GDP對(duì)數(shù)的平方納入模型。

(2)城鎮(zhèn)化水平(x2)以城市人口占總?cè)丝诘谋戎乇硎尽?/p>

(3)產(chǎn)業(yè)結(jié)構(gòu)(x3)以第三產(chǎn)業(yè)產(chǎn)值占當(dāng)年GDP的比重表示。

(4)對(duì)外開(kāi)放程度(x4)用可直接利用外商投資占當(dāng)年GDP比重表征。

(5)科技創(chuàng)新(x5)以財(cái)政用于科研的支出占當(dāng)年財(cái)政支出的比重表示，構(gòu)建回歸方程如下：

式中：因變量GDE為綠色發(fā)展效率；i為地市；β0為常數(shù)項(xiàng)，β1，β2，β3，β4，β5，β6分別為各自變量的待估參數(shù)；μ為隨機(jī)擾動(dòng)項(xiàng)。為了降低回歸模型異方差影響，該文對(duì)所有變量均進(jìn)行取對(duì)數(shù)處理，所有基礎(chǔ)數(shù)據(jù)均來(lái)源于《山東省統(tǒng)計(jì)年鑒》(2008—2018)。

3.2 回歸結(jié)果分析

運(yùn)用Stata軟件進(jìn)行Tobit分析，可得山東省綠色發(fā)展效率影響因素回歸結(jié)果如表3所示。

表3 山東省綠色發(fā)展效率影響因素回歸結(jié)果整理

(1)經(jīng)濟(jì)發(fā)展水平(x1)與綠色發(fā)展效率呈負(fù)相關(guān)關(guān)系，人均GDP對(duì)數(shù)的平方系數(shù)為正，表明山東省綠色發(fā)展效率與經(jīng)濟(jì)發(fā)展水平之間存在U型曲線關(guān)系。根據(jù)環(huán)境庫(kù)茲涅茨假說(shuō)，一個(gè)地區(qū)的經(jīng)濟(jì)發(fā)展水平與環(huán)境有著較強(qiáng)的關(guān)聯(lián)性，在經(jīng)濟(jì)發(fā)展的初期階段，隨著經(jīng)濟(jì)的增長(zhǎng)，環(huán)境污染程度加深，意味著綠色發(fā)展效率隨之降低，但隨著經(jīng)濟(jì)實(shí)力不斷增強(qiáng)，達(dá)到一定臨界點(diǎn)后，經(jīng)濟(jì)水平的提高會(huì)使環(huán)境好轉(zhuǎn)，即綠色發(fā)展效率不斷提升。從回歸結(jié)果中可以看出，現(xiàn)階段山東省經(jīng)濟(jì)發(fā)展仍處于庫(kù)茲涅茨曲線拐點(diǎn)左側(cè)，隨著經(jīng)濟(jì)發(fā)展的提升，綠色發(fā)展效率值不斷降低，因此改變現(xiàn)有粗放式經(jīng)濟(jì)發(fā)展方式成為山東省綠色發(fā)展提升的重要途徑。

(2)城鎮(zhèn)化水平(x2)與綠色發(fā)展效率呈顯著負(fù)相關(guān)關(guān)系，回歸結(jié)果表明城鎮(zhèn)化水平每上升1個(gè)百分比，綠色發(fā)展效率將下降0.23個(gè)單位，研究期限內(nèi)山東省城鎮(zhèn)化水平對(duì)于綠色發(fā)展效率具有負(fù)向作用，說(shuō)明現(xiàn)階段山東省城鎮(zhèn)化仍然注重經(jīng)濟(jì)要素的空間集聚，忽視了經(jīng)濟(jì)發(fā)展模式的提升，傳統(tǒng)的規(guī)模擴(kuò)張與總量增長(zhǎng)的發(fā)展對(duì)區(qū)域綠色發(fā)展產(chǎn)生了嚴(yán)重干擾作用，且城鎮(zhèn)居民環(huán)保意識(shí)仍較弱。

(3)產(chǎn)業(yè)結(jié)構(gòu)(x3)，即第三產(chǎn)業(yè)占比與綠色發(fā)展呈顯著正相關(guān)關(guān)系，回歸結(jié)果顯示第三產(chǎn)業(yè)占比每上升一個(gè)百分比，綠色發(fā)展效率值將提升0.13個(gè)單位，第三產(chǎn)業(yè)的加快發(fā)展是生產(chǎn)力提高和社會(huì)進(jìn)步的必然結(jié)果，第三產(chǎn)業(yè)相對(duì)第二產(chǎn)業(yè)需要消耗更少的能源、資源，對(duì)環(huán)境壓力小，因此調(diào)整產(chǎn)業(yè)結(jié)構(gòu)對(duì)綠色發(fā)展效率具有顯著的促進(jìn)作用，同時(shí)說(shuō)明山東省加快產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級(jí)，同時(shí)把轉(zhuǎn)方式調(diào)結(jié)構(gòu)的成效體現(xiàn)到提高發(fā)展的質(zhì)量效益上來(lái)，能夠有效的促進(jìn)綠色發(fā)展效率的提升。

(4)對(duì)外開(kāi)放程度(x4)與綠色發(fā)展呈顯著正相關(guān)關(guān)系，回歸結(jié)果顯示外商直接投資每上升一個(gè)單位，綠色發(fā)展效率值上升0.02個(gè)單位，現(xiàn)階段引進(jìn)外資對(duì)山東省綠色發(fā)展具有顯著的正向驅(qū)動(dòng)作用，說(shuō)明山東省著力提高外資質(zhì)量取得成效，利用有效外資為經(jīng)濟(jì)發(fā)展帶來(lái)新的活力，緩解本地經(jīng)濟(jì)發(fā)展的困境，因此山東省應(yīng)結(jié)合“一帶一路”等戰(zhàn)略，鼓勵(lì)外資投向高新技術(shù)產(chǎn)業(yè)、裝備制造業(yè)、現(xiàn)代服務(wù)業(yè)、現(xiàn)代農(nóng)業(yè)和研發(fā)領(lǐng)域，推動(dòng)開(kāi)放向更高層次發(fā)展，提升經(jīng)濟(jì)發(fā)展的質(zhì)量。

(5)科技創(chuàng)新(x5)與綠色發(fā)展效率具有較強(qiáng)的正相關(guān)關(guān)系，回歸結(jié)果顯示科技投入每上升一個(gè)單位，綠色發(fā)展效率上升0.01個(gè)單位，現(xiàn)階段科技對(duì)綠色發(fā)展效率具有重要的正向驅(qū)動(dòng)作用，說(shuō)明科技創(chuàng)新能夠促進(jìn)新興產(chǎn)業(yè)的發(fā)展，改變山東省高污染、高投入及低產(chǎn)出的發(fā)展現(xiàn)狀，促進(jìn)產(chǎn)業(yè)優(yōu)化升級(jí)，促進(jìn)資源利用效率及清潔生產(chǎn)效率的提升，因此山東省應(yīng)以科技為支撐，以創(chuàng)新驅(qū)動(dòng)發(fā)展，推動(dòng)經(jīng)濟(jì)發(fā)展質(zhì)量變革。

4 結(jié)論與討論

基于綠色發(fā)展內(nèi)涵，綜合構(gòu)建了綠色發(fā)展效率評(píng)價(jià)指標(biāo)體系，采用基于非期望產(chǎn)出的SBM模型測(cè)算山東省綠色發(fā)展效率，并利用Tobit模型對(duì)其影響因素進(jìn)行研究，結(jié)論如下：

(1)2007—2017年山東省綠色發(fā)展整體發(fā)展趨勢(shì)良好，部分地區(qū)效果不理想，各地市綠色發(fā)展效率提升趨勢(shì)的地區(qū)性差異較明顯，東部地區(qū)效率提升的趨勢(shì)快于中西部地區(qū)，可劃分為增長(zhǎng)型、穩(wěn)定型與衰減型三種類型，且各地市綠色發(fā)展效率的時(shí)間演變與山東省經(jīng)濟(jì)發(fā)展水平及經(jīng)濟(jì)發(fā)展方式的轉(zhuǎn)變息息相關(guān)。

(2)研究期限內(nèi)山東省綠色發(fā)展效率存在明顯的空間分異性，并且存在高效率區(qū)向東發(fā)展，低效率區(qū)轉(zhuǎn)向中西部的演化特征，綠色發(fā)展效率空間格局演化與經(jīng)濟(jì)發(fā)展格局具有相對(duì)一致性，隨山東省地區(qū)產(chǎn)業(yè)結(jié)構(gòu)性變革而演變，且國(guó)家級(jí)發(fā)展戰(zhàn)略的實(shí)施推動(dòng)了相關(guān)地區(qū)綠色發(fā)展效率的提升。

(3)產(chǎn)業(yè)結(jié)構(gòu)、對(duì)外開(kāi)放程度以及科技創(chuàng)新是山東省綠色發(fā)展效率優(yōu)化改善的正向驅(qū)動(dòng)因素，而現(xiàn)階段經(jīng)濟(jì)發(fā)展水平及城鎮(zhèn)化水平對(duì)綠色發(fā)展效率具有顯著的負(fù)向影響，顯示了山東省經(jīng)濟(jì)發(fā)展轉(zhuǎn)型不徹底，結(jié)構(gòu)性制約凸顯，新型城鎮(zhèn)化尚未完全落實(shí)等問(wèn)題突出。

從整體上看，山東省綠色發(fā)展效率存在著不均衡、不合理等問(wèn)題，而轉(zhuǎn)變粗放式經(jīng)濟(jì)發(fā)展方式是改善山東省綠色發(fā)展問(wèn)題根本途徑，產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級(jí)是綠色發(fā)展效率提升的重要手段，科技創(chuàng)新是綠色發(fā)展效率提升的保障，更高層次的對(duì)外開(kāi)放能推動(dòng)山東省經(jīng)濟(jì)發(fā)展質(zhì)量的提升。綠色發(fā)展與經(jīng)濟(jì)發(fā)展模式密切相關(guān)，研究發(fā)現(xiàn)山東省西部地區(qū)綠色發(fā)展效率趨于下降，表明落后地區(qū)的經(jīng)濟(jì)發(fā)展應(yīng)適應(yīng)新的工業(yè)化發(fā)展，探尋高質(zhì)量的經(jīng)濟(jì)發(fā)展模式，而非沿用傳統(tǒng)的發(fā)展模式。與以往研究相比，本研究的貢獻(xiàn)在于基于非期望產(chǎn)出的角度分析了山東省17地市的綠色發(fā)展效率的時(shí)空演化特征，并就其影響因素進(jìn)行定量分析，為山東省綠色發(fā)展優(yōu)化改善提供參考。綠色發(fā)展效率實(shí)證研究多重視邏輯推理和計(jì)量分析，而研究尺度多從大區(qū)域如省、地市的角度來(lái)研究，縣級(jí)及以下城市研究相對(duì)較少，未來(lái)需要加強(qiáng)相關(guān)研究，另外綠色發(fā)展效率時(shí)空格局演化的內(nèi)在驅(qū)動(dòng)機(jī)制比較復(fù)雜，本研究在分析綠色發(fā)展效率影響因素時(shí)，僅選取的關(guān)注度較高的5個(gè)因素，有待進(jìn)一步細(xì)化，因此未來(lái)應(yīng)側(cè)重于多尺度，多層次的綠色發(fā)展效率動(dòng)態(tài)監(jiān)測(cè)與評(píng)價(jià)。