石天琦，蔣 濤

（江蘇師范大學商學院，江蘇南京 221116）

一、引言與文獻綜述

近年來，快速的工業(yè)化和城市化進程使城市面臨持續(xù)的經(jīng)濟增長和生活環(huán)境質(zhì)量下降的困境。中國共產(chǎn)黨第十九次全國代表大會的正式報告首次明確表示要實施最嚴格的生態(tài)系統(tǒng)管理條例。以此實現(xiàn)綠色發(fā)展，并建立和完善綠色，低碳，可持續(xù)發(fā)展的經(jīng)濟體系。因此，協(xié)調(diào)好綠色發(fā)展與科技創(chuàng)新，提高綠色創(chuàng)新效率，能實現(xiàn)保持經(jīng)濟中高速增長，滿足人民日益增長的美好生活需要的目標。

長三角地區(qū)是我國制造業(yè)集中度最高的地區(qū)之一，在高強度開發(fā)階段，它面臨著非常嚴重的生態(tài)系統(tǒng)破壞和嚴重的環(huán)境污染。發(fā)行于2019 年的《長三角一體化模式進一步發(fā)展總體規(guī)劃方案文件》指出，我們需要深化改革，實現(xiàn)創(chuàng)新驅(qū)動發(fā)展，促進環(huán)境的協(xié)調(diào)管理，建設(shè)長三角綠色發(fā)展高級示范區(qū)。因此，通過研究長三角地區(qū)不同城市綠色創(chuàng)新效率的整體特征以及相關(guān)因素的影響效果，可以為制定相關(guān)政策提供更多理論建議。

通過梳理現(xiàn)有文獻，發(fā)現(xiàn)關(guān)于綠色創(chuàng)新效率的研究主要是兩個角度：一是有關(guān)衡量效率水平的研究。李金滟等（2016）[1]對長江中游28 個城市綠色創(chuàng)新效率進行了科學研究；曾冰（2018）[2]進行了中國省級綠色創(chuàng)新效率的研究；黃磊和吳傳清（2019）[3]相關(guān)研究表明，長江經(jīng)濟帶省際工業(yè)綠色創(chuàng)新發(fā)展具有協(xié)同效應(yīng)。在研究方法方面，肖黎明等（2018）[4]基于改進的隨機前沿模型，對2001—2015 年中國30個省的綠色創(chuàng)新效率進行了測度；滕堂偉等（2019）[5]采用超級SBM 模型對長三角26 個城市進行了綠色創(chuàng)新效率的測度。二是對影響綠色創(chuàng)新效率的諸多因素進行了相關(guān)研究。楊樹旺等（2018）[6]研究了長江經(jīng)濟帶各省市，認為市場開放程度、企業(yè)污染成本和技術(shù)市場成熟度對綠色創(chuàng)新效率有正向影響，產(chǎn)業(yè)結(jié)構(gòu)對綠色創(chuàng)新效率的影響不明顯；此外，李健和馬曉芳（2019）[7]探索了京津冀各城市，認為經(jīng)濟發(fā)展水平、城市文化水平、政府資助力度、產(chǎn)業(yè)結(jié)構(gòu)和城市信息化水平對提高綠色創(chuàng)新效率水平有著正向影響，但對外開放水平對效率水平的影響不明顯。

以往對綠色發(fā)展與創(chuàng)新活動的研究大多集中在省級層面，對于都市空間尺度的實證研究還有待進一步增加和豐富。事實上，城市是突破創(chuàng)新生態(tài)結(jié)構(gòu)的最主要載體，是地方政府進行環(huán)境治理的最有效單位。也可以更有效地充分調(diào)動和控制各個方面。因此，在創(chuàng)新成為國內(nèi)重要區(qū)域經(jīng)濟持續(xù)發(fā)展重要動力的背景下，對綠色發(fā)展創(chuàng)新效率和質(zhì)量的研究需要進一步豐富。因此，本文認為從省級角度對綠色創(chuàng)新的效率和質(zhì)量進行測度，不利于考慮空間擴散效應(yīng)，而城市是最適合深入研究的方向。所以，本文以2006—2017 年長三角地區(qū)27 個城市的數(shù)據(jù)為樣本，建立綜合評價體系測度了綠色創(chuàng)新效率，并使用空間計量模型來研究相關(guān)影響因素，為制定相關(guān)政策提供更多理論建議。

二、變量與數(shù)據(jù)

（一）綠色創(chuàng)新效率測度

本文基于市場投入與產(chǎn)出效率的視角，構(gòu)建了一個綜合指標體系，并結(jié)合SBM 超效率方法計算了長三角地區(qū)城市的綠色創(chuàng)新效率。指標值的選擇要考慮實際情況和數(shù)據(jù)的可得性，因此在投入指標值方面，選擇了在崗職工平均人數(shù)作為勞動投入，選擇了固定資產(chǎn)投資作為地區(qū)資本投入，選擇了工業(yè)用電量作為能源投入。從產(chǎn)出指標值來看，期望產(chǎn)出包括經(jīng)濟產(chǎn)出和創(chuàng)新能力產(chǎn)出。制造業(yè)增加值和專利技術(shù)授權(quán)量作為衡量期望產(chǎn)出的相關(guān)指標，工業(yè)三廢排放量作為非期望產(chǎn)出的相關(guān)指標，在很大程度上能夠真實反映城市環(huán)境治理水平。具體評價體系（見表1）。

表1 城市綠色創(chuàng)新效率指標體系

通常，DMU 的數(shù)量應(yīng)大于投入和產(chǎn)出主要指標總數(shù)的乘積，并且數(shù)量不應(yīng)超過投入和產(chǎn)出相關(guān)指標數(shù)量的三倍。因此，用于評價長三角地區(qū)綠色創(chuàng)新效率水平的投入產(chǎn)出關(guān)鍵指標數(shù)量的最大值是9。本文選擇了8 個相關(guān)指標，符合本文中使用的方法的基本規(guī)則。

（二）影響因素設(shè)定

城市創(chuàng)新效率水平受到許多因素的影響，例如都市經(jīng)濟的水平，創(chuàng)新技術(shù)的快速發(fā)展水平，相關(guān)的政策體系以及人們思想上的創(chuàng)新。根據(jù)現(xiàn)有研究結(jié)果和收集到的數(shù)據(jù)，結(jié)合長三角地區(qū)現(xiàn)狀，從經(jīng)濟發(fā)展的水平、污染控制的效率、外資的利用水平、產(chǎn)業(yè)結(jié)構(gòu)、政府資助的水平和城市教育的水平這幾個方面分析了對于綠色創(chuàng)新效率水平的影響效果。表2 顯示了影響因素的定義。

表2 變量定義與說明

（三）樣本選擇與數(shù)據(jù)來源

本文研究的27 個城市基于2019 年12 月1 日由中共中央和國務(wù)院發(fā)布的《長江三角洲區(qū)域一體化發(fā)展規(guī)劃綱要》。選取2006—2017 年的面板數(shù)據(jù)進行綠色創(chuàng)新效率水平的衡量和相關(guān)因素的影響效果分析。所有數(shù)據(jù)主要來源于《中國城市統(tǒng)計年鑒》以及各市的統(tǒng)計年鑒。對于個別年份數(shù)據(jù)的缺失，則采用插值法補齊。

三、模型構(gòu)建

本文從投入產(chǎn)出的角度出發(fā)，構(gòu)建了一個綜合的指標體系，并運用SBM 超效率模型計算了長三角各都市綠色創(chuàng)新的效率值。其次，通過ArcGIS、Geoda 等應(yīng)用軟件，利用自然斷裂法，在圖上表示出各城市綠色創(chuàng)新效率的空間相關(guān)性。

（一）SBM 超效率模型[8]

現(xiàn)有的研究結(jié)果主要是用數(shù)據(jù)包絡(luò)模型（DEA）和隨機前沿分析（SFA）兩個基本模型來評價效率。其中，SFA 建模中設(shè)置的產(chǎn)出項總數(shù)只能配置為1，不適合多產(chǎn)出問題的詳細分析。對多投入多產(chǎn)出的效率進行綜合分析的研究大多采用傳統(tǒng)形式的DEA 模型。然而傳統(tǒng)的模型假設(shè)所有產(chǎn)出都是期望產(chǎn)出，其中一些產(chǎn)出指標會降低效率水平，這種產(chǎn)出稱為非期望產(chǎn)出。對于非預(yù)期的產(chǎn)出，Tone（2001）問題提出了SBM 模型（Slack Based Measure）。該模型的優(yōu)點是解決了徑向模型對無效率的測量沒有包含松弛變量的問題。

在DEA 模型的分析結(jié)果中，許多決策單元會被認為是有效的。尤其是當投入和產(chǎn)出指標數(shù)量較多時，有效DMU 數(shù)量也會較多。DEA 模型得出的效率值最大為1，有效DMU 效率值相同，這些有效DMU的效率高低無法進一步區(qū)分。因此，在Li 等（2013）方法的基礎(chǔ)上，采用SBM 超效率模型，來解決在各種資源限制和存在非預(yù)期產(chǎn)出的情況下，對各個有效DMU 進行區(qū)分。

考慮到非期望產(chǎn)出的SBM 超效率模型示例：

（二）全局莫蘭指數(shù)

為了進一步分析城市之間綠色持續(xù)創(chuàng)新效率的演化過程，使用全局Moran's I 指數(shù)來測試27 個長三角城市之間是否存在空間相關(guān)性，全局Moran's I 指數(shù)表示如下：

在公式中，i 和j 代表每個都市，N 是城市的數(shù)量，xi和xj為綠色創(chuàng)新效率水平，wij是空間權(quán)重流量矩陣。莫蘭指數(shù)的取值一般介于-1～1 之間，大于0，則表示正自相關(guān)，即高值與高值相鄰、低值與低值相鄰；小于0，表示負自相關(guān)，即高值與低值相鄰。如果莫蘭指數(shù)接近于0，則表示空間分布是隨機的，不存在空間自相關(guān)。

（三）空間計量模型

空間相關(guān)關(guān)系的識別是空間計量分析的前提。而根據(jù)空間權(quán)重矩陣，對空間相關(guān)關(guān)系進行判定之后，空間數(shù)據(jù)分析的重點則轉(zhuǎn)向了對空間計量模型的選擇和估算上。空間計量模型是對納入不同空間交互效應(yīng)的普通回歸模型的拓展，主要包含以下三類：

1.空間滯后模型（SLM）

其中i 代表相應(yīng)區(qū)域，t 代表相應(yīng)年份；GIEit是綠色發(fā)展創(chuàng)新水平；是空間自回歸比例因子；Wij是空間權(quán)重矩陣；Xit為經(jīng)對數(shù)化處理后的影響因素；是空間滯后的影響因子；it代表個體固定效應(yīng)；代表隨機誤差項。此模型包含了因變量之間的內(nèi)生交互項，表明一個地區(qū)綠色創(chuàng)新效率值不僅取決于這個地區(qū)一些變量的影響，還受到鄰近地區(qū)綠色創(chuàng)新效率值的影響。

2.空間誤差模型（SEM）

SEM 模型包含了誤差項間的交互效應(yīng)，表明綠色創(chuàng)新效率值的空間相關(guān)性受到一些難以觀察到的誤差擾動項的影響。其基本表達式為：

式（4）、式（5）中uit為空間誤差自相關(guān)，為區(qū)域間溢出效應(yīng)。

3.空間杜賓模型（SDM）

SDM 模型同時包含了內(nèi)生及外生交互效應(yīng)，表示一個地區(qū)的綠色創(chuàng)新效率值不僅受到相鄰地區(qū)的綠色創(chuàng)新效率值的影響，還會受到相鄰地區(qū)的影響因素的影響。SDM 模型為：

四、實證結(jié)果分析

（一）綠色創(chuàng)新效率測算

本文利用MaxDEA 應(yīng)用軟件對投入和產(chǎn)出相關(guān)指標進行再加工，通過SBM 超效率模型得到2006—2017 年長三角27 個城市的綠色創(chuàng)新效率水平，最終結(jié)果如表3 所示。

27 個城市的綠色創(chuàng)新效率平均值最小為0.81，最大為1.09，總體而言，它正在逐年增加。從全局來看，從2006—2012 年，大部分城市的綠色創(chuàng)新效率的平均水平不斷提高，2012 年效率值達到0.97；從2013—2017 年，城市綠色創(chuàng)新效率的整體平均水平發(fā)生了較大變化。

從城市層面來看，2006—2017 年，上海、無錫、蘇州、寧波和溫州的14 個城市的效率值略高于平均標準水平，表明這些城市投入的人力和財力獲得了回報。池州、金華等城市的綠色創(chuàng)新效率高，這與其工業(yè)存在后發(fā)優(yōu)勢密切相關(guān)。在這些城市中，第一產(chǎn)業(yè)仍占很大比重，制造業(yè)比重相對較低，發(fā)展負擔相對較小。它能及時對不斷發(fā)展的綠色創(chuàng)新要求作出回應(yīng)。上海、蘇州等大城市的創(chuàng)新能力讓它們受益匪淺，近年來，隨著我國高校和科研機構(gòu)的增多，創(chuàng)新能力已位居全國各城市前列，隨著環(huán)境相關(guān)政策的實施，綠色創(chuàng)新效率得到了提高。

表3 2006—2017 年長三角地區(qū)綠色創(chuàng)新效率的實證結(jié)果

（二）空間分布及空間相關(guān)性分析

1.綠色創(chuàng)新效率空間分布。根據(jù)綠色創(chuàng)新效率的演化特征，選擇了2006 年、2011 年和2017 年三個特征控制節(jié)點，并使用了使用ArcGIS 應(yīng)用軟件的自然斷裂法[9]，把綠色創(chuàng)新效率水平分為四種類型：高水平、較高水平、較低水平和低水平（見圖1）。

如圖1 所示，長三角區(qū)域綠色創(chuàng)新效率整體空間分布不平衡，但失衡程度正在逐步減小，特別是東部區(qū)域的綠色創(chuàng)新效率水平明顯增加?？傮w而言，從2006—2017 年，城市效率水平為高水平和較高水平的數(shù)量大幅提升，2017 年達到16 個，尤其是合肥、南京、鎮(zhèn)江、寧波等的綠色創(chuàng)新效率水平進步不少，這表明在這些城市的發(fā)展過程中，考慮到了既要發(fā)展經(jīng)濟又要提高生活環(huán)境質(zhì)量。相反，另外兩個類型的城市數(shù)量迅速下降。除了鹽城市的較低水平外，馬鞍山、安慶等城市的綠色創(chuàng)新效率也都有不同程度的退步。

圖1 長三角地區(qū)綠色創(chuàng)新效率空間分布情況

2.空間相關(guān)性分析。利用Geoda 軟件測度了2006—2017 年長三角地區(qū)綠色創(chuàng)新效率的全局Moran 指數(shù)，評價城市間綠色創(chuàng)新效率是否存在空間相關(guān)性，結(jié)果如表4 所示。

表4 2006—2017 年長三角地區(qū)綠色創(chuàng)新效率的全局莫蘭指數(shù)

從表4 可以清楚地看出，2006—2017 年27 個城市的綠色創(chuàng)新效率的Moran 指數(shù)值大于0，并且大多數(shù)年份的指數(shù)值都經(jīng)過了10%的顯著性檢驗，拒絕了不存在空間相關(guān)性的假設(shè)。從時間角度來看，長江三角洲地區(qū)綠色創(chuàng)新效率水平在2006—2017 年之間波動很大，即從2006—2011 年有所下降后在2012 年有所上升，接著從2013—2017 年，效率值又急劇減小，但從總體上看，Moran 指數(shù)趨于0，長三角區(qū)域城市之間的空間相關(guān)性也有所減弱。另外，結(jié)合表4 的分析，可以看出，經(jīng)過連續(xù)十二年的發(fā)展，長三角區(qū)域的綠色創(chuàng)新朝著合作（融合）的方向進一步發(fā)展。

（三）影響因素分析

Moran 指數(shù)的實際結(jié)果表明，綠色創(chuàng)新效率水平有著空間相關(guān)性。因此，使用空間計量經(jīng)濟模型來分析綠色創(chuàng)新效率水平的影響因素。

在進入正式的分析之前，使用Matlab2018a 軟件，首先設(shè)定非空間計量模型對其進行估計，然后檢驗滯后或者殘差空間相關(guān)性，并考察空間與時間特定效應(yīng)的聯(lián)合顯著性以及隨機效應(yīng)是否可以替代它們[10]。

在進入正式的分析之前，首先設(shè)定非空間計量模型對其進行估計，然后檢驗滯后或殘差空間相關(guān)性。表5 報告了采用非空間面板數(shù)據(jù)模型的估計結(jié)果，并進行進一步的檢驗來決定采用空間滯后模型還是空間誤差模型來擬合數(shù)據(jù)。利用時間空間雙固定效應(yīng)下的普通面板OLS 估計得到的殘差序列進行LM 檢驗，得到的LMlag 值為63.73、R-LMlag 值為3.81，均通過了5%顯著性檢驗；而R-LMerror 值為2.31，未通過顯著性檢驗，因此選擇SLM 模型要比SEM 模型更為合適。為了檢查時間固定效應(yīng)聯(lián)合檢驗的原假設(shè)是否具有顯著性，執(zhí)行了似然比（LR）檢驗。檢驗的結(jié)果（統(tǒng)計量83.44）通過了1%顯著性檢驗，拒絕了原假設(shè)。但是，空間固定效應(yīng)聯(lián)合檢驗不顯著的原假設(shè)（統(tǒng)計量21.09）沒有被拒絕。以上檢驗指出應(yīng)該采用具有時間固定效應(yīng)的空間滯后模型。

表5 無空間交互效應(yīng)的面板數(shù)據(jù)模型估計

由于LM 檢驗并沒有考慮到SDM 模型的適用性，因此，為了檢驗空間杜賓模型能否簡化成簡單的空間誤差模型或空間滯后模型的原假設(shè)，即H0∶茲=0 和H0∶茲+=0，可以執(zhí)行Wald 或者LR 檢驗。因此，使用Matlab2018a 軟件，在時間相對固定的影響下，得到時間固定效應(yīng)下SDM 模型中Wald 和LR檢驗的檢驗值，據(jù)此判斷SDM 模型是否可以簡化，具體結(jié)果如表6 所示。

表6 時間固定SDM 模型的Wald 和LR 檢驗結(jié)果表

由表6 可知，空間誤差的Wald 檢驗和LR 檢驗，通過了1%顯著性檢驗，拒絕了原假設(shè)，即SDM模型不能簡化；而空間滯后的Wald 檢驗和LR 檢驗，沒有通過10%顯著性檢驗，不能拒絕原假設(shè)。利用極大似然法（ML）估計不同交互效應(yīng)下的SLM 模型，具體結(jié)果如表7 所示。

表7 不同交互效應(yīng)下SLM 模型的ML 估計結(jié)果表

空間計量模型下，log-likelihood 值的大小表明了擬合數(shù)據(jù)效果的優(yōu)良。從表4 可以看出，時間固定效應(yīng)的SLM 模型的對數(shù)似然值最高，且其擬合優(yōu)度R2與離散度滓2結(jié)果相對較好，再次證實了本研究中時間固定SLM 模型的優(yōu)越性。

由表7 可知，經(jīng)濟發(fā)展水平對綠色創(chuàng)新效率的貢獻度顯著為0.202，這意味著人均生產(chǎn)總值每增加1 個百分點，城市綠色創(chuàng)新效率就增加0.202 個百分點。這表明，人均GDP 的提高，將使許多人對環(huán)境有更嚴格的要求，對綠色低碳產(chǎn)品的日常消費需求也將增加；一般工業(yè)固定廢物綜合利用率每增加1%，綠色創(chuàng)新效率同向變動0.197%且通過了5%的顯著性檢驗；外資利用水平的系數(shù)為0.079 且在1%的顯著性水平上顯著，表明當年實際利用外資金額每增加1%，綠色創(chuàng)新效率同向變動0.079%；產(chǎn)業(yè)結(jié)構(gòu)選取的是第二產(chǎn)業(yè)增加值占GDP 比重，該指標在5%的水平下通過了顯著性檢驗，且影響系數(shù)為-0.257，對綠色創(chuàng)新效率有顯著的負向影響，說明第二產(chǎn)業(yè)的繼續(xù)發(fā)展不利于綠色創(chuàng)新活動的產(chǎn)生和綠色創(chuàng)新效率的提高；政府資助力度對綠色創(chuàng)新效率的影響較弱，但仍具有促進作用；城市教育水平選取的是高等學校在校學生數(shù)，對綠色創(chuàng)新效率的影響因素在5%顯著性水平上顯著為0.050，可知城市教育水平可以有力地推動城市綠色創(chuàng)新效率水平的提高。

五、結(jié)論與建議

本文基于長三角地區(qū)27 個城市2006—2017 年的面板數(shù)據(jù)，通過建立空間計量模型得出了相關(guān)影響因素作用效果的結(jié)論，實證結(jié)果及相關(guān)建議如下：

首先，長三角區(qū)域各城市的綠色創(chuàng)新效率正在朝著相互協(xié)調(diào)（整合）的方向發(fā)展。從城市角度來看，27 個城市的綠色創(chuàng)新效率水平差異正在逐步縮小?？傮w上，可以看出江蘇省可持續(xù)綠色創(chuàng)新效率得到了快速提高，浙江省效率水平穩(wěn)步提高且更加均衡，安徽省綠色創(chuàng)新效率則是先降后升。在這種情況下，上海應(yīng)該繼續(xù)發(fā)揮帶頭作用，蘇浙皖要相互學習，揚長避短，加強合作，打破壁壘，推動區(qū)域融合發(fā)展。

其次，經(jīng)濟發(fā)展水平對城市綠色創(chuàng)新效率有著顯著促進作用。然后，第二產(chǎn)業(yè)增加值占GDP 比重對綠色創(chuàng)新效率有著負向驅(qū)動作用，一般工業(yè)固定廢物綜合利用率對綠色創(chuàng)新效率有著正向影響，外資利用水平對綠色創(chuàng)新效率有較強的正向作用，政府資助水平雖然沒有通過顯著性檢驗但影響作用為正，萬人在校大學生人數(shù)正向促進綠色發(fā)展創(chuàng)新效率的提高。隨著我國經(jīng)濟進入中高速增長階段以及美國對于我國發(fā)展的遏制，以長三角地區(qū)為首的各城市有著較高的經(jīng)濟發(fā)展水平和很好的科研水平，政府要加大對研發(fā)的投入，實現(xiàn)產(chǎn)學研一體化發(fā)展，加快創(chuàng)新成果轉(zhuǎn)化，加強城市間合作，加快推動創(chuàng)新人才引進與整合，加大基礎(chǔ)研究和高新技術(shù)的研究，提高綠色發(fā)展水平，更好地滿足人民日益增長的美好生活需要。