京津冀工業(yè)綠色發(fā)展的空間效應(yīng)分析
——基于動(dòng)態(tài)SDM的供給側(cè)因素分析

王韶華，王 菲，何美璇

(1.燕山大學(xué) 經(jīng)濟(jì)管理學(xué)院，河北 秦皇島 066004；2. 燕山大學(xué) 區(qū)域經(jīng)濟(jì)發(fā)展研究中心，河北 秦皇島 066004)

一、研究背景及文獻(xiàn)綜述

加強(qiáng)生態(tài)文明建設(shè)是打造京津冀世界級(jí)城市群的重要突破口，綠色發(fā)展是推進(jìn)京津冀生態(tài)文明建設(shè)的題中之義，產(chǎn)業(yè)綠色化是構(gòu)成京津冀綠色發(fā)展的內(nèi)核。其中，工業(yè)綠色發(fā)展不僅可以根治京津冀生態(tài)環(huán)境問題，破解發(fā)展與資源環(huán)境不協(xié)調(diào)、區(qū)域發(fā)展不平衡等矛盾，而且會(huì)推動(dòng)京津冀區(qū)域在新一輪工業(yè)革命中搶占領(lǐng)先地位，可見推進(jìn)京津冀工業(yè)綠色發(fā)展取得更大進(jìn)展意義重大。但在“三期疊加”的“經(jīng)濟(jì)新常態(tài)”下，京津冀工業(yè)綠色發(fā)展當(dāng)前面臨著要素成本上升、要素生產(chǎn)效率低下、要素流動(dòng)不順暢導(dǎo)致的傳統(tǒng)要素驅(qū)動(dòng)效應(yīng)下降而新動(dòng)能培育不足、資源環(huán)境約束嚴(yán)重、產(chǎn)品質(zhì)量低下、產(chǎn)能過剩、區(qū)域間梯度差異顯著、轉(zhuǎn)型升級(jí)難等困境[1-2]，必須堅(jiān)持供給側(cè)結(jié)構(gòu)性改革，加強(qiáng)要素創(chuàng)新、產(chǎn)業(yè)升級(jí)、制度改革才能走出一條科技含量高、經(jīng)濟(jì)效益好、資源消耗低、環(huán)境污染少、人力資源優(yōu)勢(shì)得到充分發(fā)揮的新型綠色工業(yè)化道路[3]。綠色發(fā)展和供給側(cè)結(jié)構(gòu)性改革是統(tǒng)一的，共同構(gòu)成了新常態(tài)下京津冀人與生態(tài)和諧發(fā)展，經(jīng)濟(jì)高質(zhì)量發(fā)展的完整體系。其中，綠色發(fā)展是供給側(cè)結(jié)構(gòu)性改革的方向，供給側(cè)結(jié)構(gòu)性改革為綠色發(fā)展提供動(dòng)力。

已有研究將綠色發(fā)展的內(nèi)涵從最初的以生態(tài)為核心拓展至經(jīng)濟(jì)—生態(tài)—社會(huì)協(xié)調(diào)發(fā)展，強(qiáng)調(diào)發(fā)展目標(biāo)的多元化[4-6]。關(guān)于綠色發(fā)展的測(cè)度也從單一指標(biāo)，如綠色全要素生產(chǎn)率測(cè)度[7-8]、綠色經(jīng)濟(jì)效率測(cè)算[9-10]，擴(kuò)展至多指標(biāo)綜合測(cè)度，比較權(quán)威的有北京師范大學(xué)科學(xué)發(fā)展觀與經(jīng)濟(jì)可持續(xù)發(fā)展研究基地等[11]構(gòu)建的中國綠色發(fā)展指數(shù)、四部委[12]聯(lián)合發(fā)布的《綠色發(fā)展指標(biāo)體系》等。在此基礎(chǔ)上，主要通過構(gòu)建計(jì)量經(jīng)濟(jì)模型考察綠色發(fā)展與環(huán)境規(guī)制[13]、技術(shù)創(chuàng)新[7，14-15]、勞動(dòng)力[14]、資本[15]、能源[7，16]、產(chǎn)業(yè)結(jié)構(gòu)[17]、FDI[18]、城市化水平[14]、教育投資[18]等的關(guān)系。一般而言，研究關(guān)于資本密集度、人力資本、產(chǎn)業(yè)結(jié)構(gòu)輕型化、FDI、城市化水平、教育投資等對(duì)綠色增長的影響能夠達(dá)成基本共識(shí)，但由于時(shí)空、變量等的選取不同導(dǎo)致環(huán)境規(guī)制、技術(shù)創(chuàng)新等對(duì)綠色全要素生產(chǎn)率或綠色轉(zhuǎn)型的結(jié)論并不唯一，如彭星，李斌[13]經(jīng)過研究發(fā)現(xiàn)不同類型不同地區(qū)的環(huán)境規(guī)制對(duì)工業(yè)綠色轉(zhuǎn)型的影響存在異質(zhì)性，經(jīng)濟(jì)激勵(lì)型和自愿意識(shí)型環(huán)境規(guī)制能夠促進(jìn)東部地區(qū)綠色轉(zhuǎn)型，而中西部地區(qū)對(duì)命令控制型環(huán)境規(guī)制更敏感；萬倫來，朱琴[15]認(rèn)為在低研發(fā)投入行業(yè)，研發(fā)投入的增加不利于綠色技術(shù)進(jìn)步的提升，國內(nèi)技術(shù)轉(zhuǎn)移作用有限，而在高研發(fā)投入行業(yè)國內(nèi)技術(shù)轉(zhuǎn)移利于綠色技術(shù)進(jìn)步。

綜上所述，雖然國內(nèi)外學(xué)者對(duì)綠色發(fā)展的研究已經(jīng)取得了很大進(jìn)展，為我們開展研究提供了參考，但還存在一些不足：首先，雖然已經(jīng)構(gòu)建了較權(quán)威的綠色發(fā)展指標(biāo)體系，但缺乏對(duì)京津冀這一重大國家戰(zhàn)略發(fā)展區(qū)域和工業(yè)這一重要部門的針對(duì)性研究；其次，雖然已經(jīng)分析了個(gè)別供給側(cè)要素對(duì)綠色發(fā)展的影響，但缺乏對(duì)供給側(cè)要素與綠色發(fā)展間關(guān)系的系統(tǒng)研究。有鑒于此，本文從綠色發(fā)展的內(nèi)涵出發(fā)，參考權(quán)威評(píng)價(jià)指標(biāo)體系，兼顧主觀與客觀構(gòu)建工業(yè)綠色發(fā)展測(cè)度指標(biāo)體系，在考慮空間相關(guān)性的基礎(chǔ)上，從要素、產(chǎn)業(yè)和制度三個(gè)層面探討京津冀工業(yè)綠色發(fā)展的供給側(cè)驅(qū)動(dòng)力，以期推動(dòng)供給側(cè)改革與工業(yè)綠色發(fā)展的深度融合，一方面能夠?yàn)榫G色發(fā)展研究提供一些創(chuàng)新視角，為供給側(cè)改革的實(shí)證研究提供思路；另一方面有利于促進(jìn)供給側(cè)改革政策與綠色發(fā)展政策的統(tǒng)一性，為高質(zhì)量新動(dòng)能的培育提供方向和思路，從而促進(jìn)京津冀協(xié)同發(fā)展取得新的更大進(jìn)展。

二、京津冀工業(yè)綠色發(fā)展測(cè)度

(一)指標(biāo)選取

王韶華等[19]在四部委發(fā)布的《綠色發(fā)展指標(biāo)體系》基礎(chǔ)上構(gòu)建了京津冀工業(yè)綠色發(fā)展指標(biāo)體系，并通過實(shí)證分析論證了其合理性。由于通過《中國統(tǒng)計(jì)年鑒》《北京統(tǒng)計(jì)年鑒》《天津統(tǒng)計(jì)年鑒》《河北經(jīng)濟(jì)年鑒》《中國環(huán)境統(tǒng)計(jì)年鑒》《中國城市統(tǒng)計(jì)年鑒》等僅能完整獲取2004—2016年京津冀工業(yè)化學(xué)需氧量、工業(yè)氨氮排放量、工業(yè)危險(xiǎn)廢物處置率等指標(biāo)數(shù)據(jù)，為了進(jìn)一步增強(qiáng)指標(biāo)體系的系統(tǒng)性、科學(xué)性、可操作性，本文遵循科學(xué)性、數(shù)據(jù)可得性、研究相關(guān)性等原則，將工業(yè)環(huán)境質(zhì)量中的工業(yè)化學(xué)需氧量排放總量降低率替換為工業(yè)廢水排放總量降低率，將工業(yè)氨氮排放總量降低率和工業(yè)二氧化硫排放總量降低率替換為工業(yè)廢氣排放總量降低率，并補(bǔ)充工業(yè)固體廢物產(chǎn)生總量降低率；將工業(yè)環(huán)境治理中的工業(yè)危險(xiǎn)廢物處置率刪除；其他指標(biāo)保持不變。最終構(gòu)建了包括工業(yè)資源利用、工業(yè)環(huán)境質(zhì)量、工業(yè)環(huán)境治理和工業(yè)經(jīng)濟(jì)增長等4個(gè)準(zhǔn)則層15個(gè)評(píng)價(jià)指標(biāo)的區(qū)域工業(yè)綠色發(fā)展測(cè)度指標(biāo)體系，如表1所示。

表1 京津冀工業(yè)綠色發(fā)展測(cè)度指標(biāo)體系

各指標(biāo)間的邏輯關(guān)系如下：

(1)工業(yè)資源利用。對(duì)資源的依賴性和資源利用效率是劃分工業(yè)化進(jìn)程的核心指標(biāo)。當(dāng)前經(jīng)濟(jì)發(fā)展面臨著資源短缺以及資源不合理利用導(dǎo)致的環(huán)境污染等難題，因此降低資源利用總量和提高資源利用效率是提高工業(yè)綠色發(fā)展的基礎(chǔ)[10]。總量方面利用工業(yè)能源消費(fèi)總量和工業(yè)用水總量?jī)蓚€(gè)指標(biāo)加以反映，效率方面利用單位工業(yè)增加值能源消耗降低率和單位工業(yè)增加值用水降低率兩個(gè)指標(biāo)加以反映。

(2)工業(yè)環(huán)境質(zhì)量和工業(yè)環(huán)境治理?！熬G色”一詞最初引入經(jīng)濟(jì)學(xué)領(lǐng)域的主要目的是利用經(jīng)濟(jì)手段解決環(huán)境污染問題，因此提高環(huán)境質(zhì)量是工業(yè)綠色發(fā)展的首要目標(biāo)[5]。工業(yè)污染源主要包括廢水、廢氣和固體廢棄物，其中廢水排放利用工業(yè)廢水排放總量降低率指標(biāo)加以反映，廢氣排放利用單位工業(yè)增加值二氧化碳排放降低率和工業(yè)廢氣排放總量降低率兩個(gè)指標(biāo)加以反映，固體廢棄物用工業(yè)固體廢物產(chǎn)生總量降低率指標(biāo)加以反映。工業(yè)環(huán)境治理是提高環(huán)境質(zhì)量的主要手段，工業(yè)環(huán)境治理應(yīng)綜合運(yùn)用源頭防治、過程控制和末端治理相結(jié)合，因此主要利用一般工業(yè)固體廢物綜合利用率、工業(yè)環(huán)境污染治理完成投資占工業(yè)增加值和工業(yè)廢水治理占工業(yè)環(huán)境污染治理完成投資比重三個(gè)指標(biāo)加以反映。

(3)工業(yè)經(jīng)濟(jì)增長。提高增長質(zhì)量，促進(jìn)經(jīng)濟(jì)、資源與環(huán)境的協(xié)調(diào)統(tǒng)一，實(shí)現(xiàn)經(jīng)濟(jì)可持續(xù)發(fā)展是工業(yè)綠色發(fā)展的最終目標(biāo)。經(jīng)濟(jì)的綠色高質(zhì)量發(fā)展體現(xiàn)在數(shù)量和質(zhì)量的統(tǒng)一[20]。數(shù)量方面利用工業(yè)增加值增長率加以反映，質(zhì)量方面利用工業(yè)資產(chǎn)負(fù)債率、工業(yè)科技創(chuàng)新強(qiáng)度和工業(yè)綠色產(chǎn)業(yè)比重三個(gè)指標(biāo)加以反映。

(二)權(quán)重計(jì)算

計(jì)算指標(biāo)權(quán)重的方法主要有主觀賦權(quán)法和客觀賦權(quán)法，本文將綜合運(yùn)用兩類方法，兼顧主觀賦權(quán)法和客觀賦權(quán)法的優(yōu)勢(shì)。其中主觀賦權(quán)法主要在四部委聯(lián)合發(fā)布的《綠色發(fā)展指標(biāo)體系》的基礎(chǔ)上進(jìn)行計(jì)算，客觀賦權(quán)法選擇粗糙集知識(shí)約簡(jiǎn)理論，該方法的優(yōu)勢(shì)在于能夠根據(jù)信息的不確定性約簡(jiǎn)相對(duì)不重要或不必要的指標(biāo)。

(1)主觀賦權(quán)法

① 依據(jù)四部委聯(lián)合發(fā)布的《綠色發(fā)展指標(biāo)體系》，得到各三級(jí)指標(biāo)在四部委指標(biāo)體系中的重要性εXij；

② 對(duì)各三級(jí)指標(biāo)重要性進(jìn)行歸一化處理，得到相應(yīng)指標(biāo)的四部委綜合權(quán)重；

φXij=εXij/∑jεXij

③ 根據(jù)四部委綜合權(quán)重，分別對(duì)各二級(jí)指標(biāo)下的三級(jí)指標(biāo)進(jìn)行歸一化處理，得到相應(yīng)指標(biāo)的四部委單層權(quán)重；

④ 根據(jù)二級(jí)指標(biāo)權(quán)重、三級(jí)指標(biāo)單層權(quán)重和綜合權(quán)重間的數(shù)量關(guān)系，可得到各二級(jí)指標(biāo)的四部委權(quán)重。

(2)粗糙集理論[21]

① 將各指標(biāo)視作等價(jià)關(guān)系(R)，對(duì)所有評(píng)價(jià)對(duì)象進(jìn)行分類，形成等價(jià)類(U/R)；

②將同一層所有等價(jià)類的交集(ind(P))與該層去除某等價(jià)類(指標(biāo))后其它等價(jià)類的交集(U/ind(P-{P}))進(jìn)行對(duì)比，若相同則表示不必要，刪除該指標(biāo)；若不同則表示必要，保留該指標(biāo)，并利用U/ind(P-{P})不能夠準(zhǔn)確劃分到ind(P)的對(duì)象個(gè)數(shù)與所有評(píng)價(jià)對(duì)象的個(gè)數(shù)之比表示該指標(biāo)的重要性(δXij)；

(3)組合權(quán)重

設(shè)wXij為粗糙集理論計(jì)算得到的綜合權(quán)重，φXij為四部委綜合權(quán)重，則最終的組合權(quán)重為：

θXij=αwXij+(1-α)φXij

(三)測(cè)度模型構(gòu)建

將不同測(cè)度對(duì)象各指標(biāo)的標(biāo)準(zhǔn)化數(shù)據(jù)與對(duì)應(yīng)指標(biāo)的組合權(quán)重加權(quán)求和，得到各對(duì)象的測(cè)度值，即

(1)

其中，正向指標(biāo)的數(shù)據(jù)標(biāo)準(zhǔn)化公式為：

逆向指標(biāo)的數(shù)據(jù)標(biāo)準(zhǔn)化公式為：

(四)數(shù)據(jù)收集與處理

基于科學(xué)性、數(shù)據(jù)可得性等原則，對(duì)2003—2017年京津冀工業(yè)綠色發(fā)展進(jìn)行實(shí)證分析。

(1)工業(yè)能源消費(fèi)總量和單位工業(yè)增加值能源消費(fèi)降低率：北京、天津2003—2017年的工業(yè)能源消費(fèi)總量可以分別從歷年《北京統(tǒng)計(jì)年鑒》和《天津統(tǒng)計(jì)年鑒》直接獲取，河北省2003—2017年的工業(yè)能源消費(fèi)總量選取終端能源消費(fèi)總量來表示，可從歷年《河北經(jīng)濟(jì)年鑒》直接獲取；利用2002年不變價(jià)工業(yè)增加值計(jì)算得到單位工業(yè)增加值能源消費(fèi)降低率。

(2)工業(yè)用水總量和單位工業(yè)增加值用水降低率：北京2003—2017年的工業(yè)用水總量可以從歷年《北京統(tǒng)計(jì)年鑒》中直接獲取，天津和河北2003年的工業(yè)用水總量從《中國統(tǒng)計(jì)年鑒》中獲取，2004—2017年的工業(yè)用水總量可以從歷年《中國環(huán)境統(tǒng)計(jì)年鑒》中獲?。焕?002年不變價(jià)工業(yè)增加值計(jì)算得到單位工業(yè)增加值用水降低率。

(3)單位工業(yè)增加值二氧化碳排放降低率：從2004—2018年的《中國能源統(tǒng)計(jì)年鑒》中獲得京、津、冀2003—2017年的工業(yè)能源消耗實(shí)物量，參考《IPCC溫室氣體排放清單指南》中的估算方法可以分別計(jì)算出京、津、冀2003—2017年的工業(yè)二氧化碳排放總量，利用2002年不變價(jià)工業(yè)增加值計(jì)算得到單位工業(yè)增加值二氧化碳排放降低率。

(4)工業(yè)廢水排放總量降低率和工業(yè)廢氣排放總量降低率：京、津、冀2003—2017年工業(yè)廢水排放總量、工業(yè)二氧化硫排放總量和工業(yè)煙塵排放總量從相應(yīng)年份的《中國城市統(tǒng)計(jì)年鑒》中獲取。

(5)工業(yè)固體廢物產(chǎn)生總量降低率：從2004—2018年的《中國統(tǒng)計(jì)年鑒》中獲取京、津、冀一般工業(yè)固體廢物產(chǎn)生量和工業(yè)危險(xiǎn)廢物產(chǎn)生量，相加得到工業(yè)固體廢物產(chǎn)生量。

(6)一般工業(yè)固體廢物綜合利用率：天津2003—2017年的一般工業(yè)固體廢物綜合利用率可以從歷年《天津統(tǒng)計(jì)年鑒》直接得到，北京和河北的一般工業(yè)固體廢物綜合利用量和一般工業(yè)固體廢物產(chǎn)生量可從對(duì)應(yīng)年份的《中國統(tǒng)計(jì)年鑒》中獲得，然后計(jì)算得到一般工業(yè)固體廢物綜合利用率。

(7)工業(yè)環(huán)境污染治理完成投資占工業(yè)增加值比重和工業(yè)廢水治理占工業(yè)環(huán)境污染治理完成投資的比重：京、津、冀2003—2017年的工業(yè)環(huán)境污染治理完成投資、工業(yè)廢水治理完成投資可以從歷年《中國統(tǒng)計(jì)年鑒》中直接得到；工業(yè)增加值均以2002年不變價(jià)處理。

(8)工業(yè)增加值增長率：分別從《北京統(tǒng)計(jì)年鑒》《天津統(tǒng)計(jì)年鑒》和《河北經(jīng)濟(jì)年鑒》中直接獲取京、津、冀2003—2017年的工業(yè)增加值，以2002年不變價(jià)處理。

(9)工業(yè)資產(chǎn)負(fù)債率：北京和天津2003—2017年的工業(yè)資產(chǎn)負(fù)債率可以分別從歷年《北京統(tǒng)計(jì)年鑒》和《天津統(tǒng)計(jì)年鑒》中直接獲得，河北2003—2017年的數(shù)據(jù)需要從歷年《中國工業(yè)統(tǒng)計(jì)年鑒》中獲取工業(yè)企業(yè)的負(fù)債總額和工業(yè)企業(yè)的資產(chǎn)總額，通過計(jì)算得到。

(10)工業(yè)科技創(chuàng)新強(qiáng)度：由規(guī)模以上工業(yè)R&D經(jīng)費(fèi)內(nèi)部支出與規(guī)模以上主營業(yè)務(wù)收入之比得到。其中，京、津、冀2003—2017年的規(guī)模以上工業(yè)R&D經(jīng)費(fèi)內(nèi)部支出可以從歷年《中國科技統(tǒng)計(jì)年鑒》中直接獲得，京、津、冀2003—2017年的規(guī)模以上主營業(yè)務(wù)收入分別從歷年《北京統(tǒng)計(jì)年鑒》《天津統(tǒng)計(jì)年鑒》和《河北經(jīng)濟(jì)年鑒》中直接獲得。

(11)工業(yè)綠色產(chǎn)業(yè)比重：京、津、冀2001—2016年高技術(shù)產(chǎn)業(yè)利潤總額可分別通過歷年《中國高技術(shù)產(chǎn)業(yè)統(tǒng)計(jì)年鑒》和《中國工業(yè)統(tǒng)計(jì)年鑒》直接獲取，北京2017年高技術(shù)產(chǎn)業(yè)利潤可通過《北京統(tǒng)計(jì)年鑒》直接獲取，津、冀2017年數(shù)據(jù)需要將醫(yī)藥制造業(yè)、鐵路、船舶、航空航天和其他運(yùn)輸設(shè)備制造業(yè)、計(jì)算機(jī)、通信和其他電子設(shè)備制造業(yè)以及儀器儀表制造業(yè)加總得到，相關(guān)數(shù)據(jù)分別通過該年的《天津統(tǒng)計(jì)年鑒》和《河北經(jīng)濟(jì)年鑒》直接獲取；京、津、冀2001—2017年規(guī)模以上工業(yè)利潤總額可分別從歷年的《北京統(tǒng)計(jì)年鑒》《天津統(tǒng)計(jì)年鑒》和《河北經(jīng)濟(jì)年鑒》中直接獲取。

(五)計(jì)算結(jié)果

根據(jù)主觀賦權(quán)法步驟得到四部委二級(jí)指標(biāo)權(quán)重、三級(jí)指標(biāo)單層權(quán)重和綜合權(quán)重如表1所示。根據(jù)粗糙集理論步驟，利用SPSS 16.0對(duì)數(shù)據(jù)進(jìn)行K-均值聚類分析，將其分為4類；對(duì)指標(biāo)進(jìn)行約簡(jiǎn)；最后計(jì)算得到各指標(biāo)單層權(quán)重和綜合權(quán)重，如表1所示。

取α=0.5，可得到各指標(biāo)組合權(quán)重，如表1所示。根據(jù)公式(1)，可得到2003—2017年京津冀工業(yè)資源利用、工業(yè)環(huán)境質(zhì)量、工業(yè)環(huán)境治理、工業(yè)經(jīng)濟(jì)增長以及工業(yè)綠色發(fā)展得分，如圖1所示。

(a) 2003—2017年北京市工業(yè)綠色發(fā)展

(b) 2003—2017年天津市工業(yè)綠色發(fā)展

(c) 2003—2017年河北省工業(yè)綠色發(fā)展

三、 供給側(cè)因素驅(qū)動(dòng)京津冀工業(yè)綠色發(fā)展空間計(jì)量分析

供給側(cè)改革作為新常態(tài)下推進(jìn)中國經(jīng)濟(jì)高質(zhì)量發(fā)展的必由之路，必然會(huì)對(duì)工業(yè)綠色發(fā)展起到積極的促進(jìn)作用。

(一)供給側(cè)因素分析框架構(gòu)建

在遵循系統(tǒng)性、科學(xué)性、數(shù)據(jù)可獲性等原則的基礎(chǔ)上，依據(jù)供給側(cè)結(jié)構(gòu)性改革的內(nèi)涵，從要素、產(chǎn)業(yè)和制度三個(gè)維度構(gòu)建包括3個(gè)二級(jí)指標(biāo)10個(gè)三級(jí)指標(biāo)的區(qū)域工業(yè)部門供給側(cè)效率評(píng)價(jià)指標(biāo)體系[22]，如表2所示，并得到2003—2017年京津冀工業(yè)部門供給側(cè)效率及要素效率、產(chǎn)業(yè)效率、制度效率得分。

表2 區(qū)域工業(yè)部門供給側(cè)效率評(píng)價(jià)指標(biāo)體系

(二)工業(yè)綠色發(fā)展空間相關(guān)性檢驗(yàn)

京津冀工業(yè)綠色發(fā)展的空間相關(guān)性直接影響計(jì)量模型的選擇，因此利用莫蘭指數(shù)考察京津冀工業(yè)綠色發(fā)展是否存在空間相關(guān)性。公式為：

(2)

借助Geoda軟件計(jì)算出2003—2017年京津冀工業(yè)綠色發(fā)展的全局莫蘭指數(shù)，計(jì)算結(jié)果見表3。2003—2017年京津冀工業(yè)綠色發(fā)展的莫蘭指數(shù)均小于-0.2，并且全部通過了1%的顯著性檢驗(yàn)，說明京津冀工業(yè)綠色發(fā)展具有顯著的負(fù)空間相關(guān)性和空間差異性，京津冀的工業(yè)綠色發(fā)展此消彼長；京津冀工業(yè)綠色發(fā)展的莫蘭指數(shù)在波動(dòng)中接近于0，京津冀之間的工業(yè)綠色發(fā)展的差距在逐漸縮小。

表3 2003—2017年京津冀工業(yè)綠色發(fā)展莫蘭指數(shù)及檢驗(yàn)結(jié)果

(三)供給側(cè)因素驅(qū)動(dòng)京津冀工業(yè)綠色發(fā)展的動(dòng)態(tài)SDM構(gòu)建

由于京津冀工業(yè)綠色發(fā)展具有顯著的空間相關(guān)性，且隨著京津冀協(xié)同發(fā)展戰(zhàn)略及供給側(cè)改革的深入推進(jìn)，京津冀工業(yè)綠色發(fā)展不僅受到區(qū)域內(nèi)本單元供給側(cè)因素的影響，還會(huì)受到其他單元供給側(cè)因素的影響，同時(shí)為考察工業(yè)綠色發(fā)展是否具有動(dòng)態(tài)性和延續(xù)性[23]，引入被解釋變量的時(shí)間滯后項(xiàng)，并考慮時(shí)間滯后項(xiàng)的空間相關(guān)性，構(gòu)建動(dòng)態(tài)空間杜賓模型，其公式為：

Ynt=kYnt-1+θWYnt-1+ρWYnt+βiXnit+

λWXnit+πWunt+δ

(3)

式中，n和t分別用于區(qū)分不同的研究單元和時(shí)間；Ynt和Ynt-1分別表示被解釋變量和被解釋變量的時(shí)間滯后項(xiàng)，k和θ分別是被解釋變量的時(shí)間滯后系數(shù)和時(shí)間滯后項(xiàng)的空間滯后系數(shù)；ρ表示空間滯后系數(shù)，反映了本單元與其他單元間研究對(duì)象的空間交互作用；Xnit表示解釋變量，βi為相應(yīng)解釋變量的系數(shù)，λi為相應(yīng)解釋變量的空間滯后系數(shù)；unt表示擾動(dòng)項(xiàng)，π為擾動(dòng)項(xiàng)的空間滯后系數(shù)；W為地理距離空間權(quán)重矩陣。

由于空間杜賓模型中的解釋變量的系數(shù)并不能有效反映解釋變量與被解釋變量間的關(guān)系，因此本文通過偏微分矩陣運(yùn)算求得解釋變量對(duì)被解釋變量的空間直接效應(yīng)與空間溢出效應(yīng)。式(3)可轉(zhuǎn)化為：

Ynt=(I-θW)-1kYnt-1+(I-θW)-1(ρWYnt+βiXnit+

λiWXnit)+(I-θW)-1δ

(4)

由式(4)求偏導(dǎo)數(shù)得到短期效應(yīng)和長期效應(yīng)的公式分別為：

(5)

[β1iIn+β2iW]

(6)

利用2003—2017年京津冀工業(yè)綠色發(fā)展測(cè)度值(Ynt)擬合方程，同樣采用k=2的鄰近規(guī)則確定空間權(quán)重矩陣W，則Ynt-1和WYnt-1分別表示工業(yè)綠色發(fā)展的時(shí)間滯后項(xiàng)及其空間滯后項(xiàng)，是反映工業(yè)綠色發(fā)展是否具有動(dòng)態(tài)延續(xù)性的重要指標(biāo)；由于供給側(cè)效率三級(jí)指標(biāo)數(shù)量較多，為保證計(jì)量模型的科學(xué)性與準(zhǔn)確性，在此利用二級(jí)指標(biāo)要素效率(fe)、產(chǎn)業(yè)效率(ie)和制度效率(se)作為解釋變量構(gòu)建模型，wfe、wie和wse分別表示相鄰區(qū)域的要素效率、產(chǎn)業(yè)效率和制度效率對(duì)本區(qū)域工業(yè)綠色發(fā)展的影響。

(四)擬合結(jié)果及分析

考慮不同類型的解釋變量時(shí)間滯后項(xiàng)的影響，利用Stata軟件擬合三種動(dòng)態(tài)杜賓模型，分別為k不為0，θ為0的動(dòng)態(tài)SDM1模型；k為0，θ不為0的動(dòng)態(tài)SDM2模型；k和θ均不為0的動(dòng)態(tài)SDM3模型，結(jié)果見表4。

表4 動(dòng)態(tài)模型結(jié)果

如表4所示，動(dòng)態(tài)SDM1、動(dòng)態(tài)SDM2和動(dòng)態(tài)SDM3中解釋變量的影響方向全部一致，保證了分析結(jié)果的穩(wěn)健性。其中，動(dòng)態(tài)SDM2的最大似然值LOGL明顯高于其他模型，說明考慮時(shí)空滯后項(xiàng)的動(dòng)態(tài)SDM2中解釋變量的顯著性優(yōu)于動(dòng)態(tài)SDM1和動(dòng)態(tài)SDM3，因此運(yùn)用SDM2分析供給側(cè)效率驅(qū)動(dòng)京津冀工業(yè)綠色發(fā)展的空間效應(yīng)更加合理。被解釋變量的空間滯后項(xiàng)系數(shù)為正，且通過了1%顯著性檢驗(yàn)，說明京津冀工業(yè)綠色發(fā)展具有較強(qiáng)的空間依賴性和正的空間溢出效應(yīng)；但時(shí)空滯后項(xiàng)系數(shù)為負(fù)，且通過了顯著性檢驗(yàn)，說明區(qū)域內(nèi)其他單元上一期工業(yè)綠色發(fā)展會(huì)對(duì)本單元當(dāng)期工業(yè)綠色發(fā)展產(chǎn)生抑制作用，可能的原因在于京津冀區(qū)域行政色彩濃厚，資源利用、環(huán)境質(zhì)量等受政府調(diào)控的影響較大，但過去一段時(shí)間津冀均處于工業(yè)化中期階段，工業(yè)發(fā)展伴隨著高消耗、高排放，工業(yè)增長與資源、環(huán)境間關(guān)系很難實(shí)現(xiàn)協(xié)調(diào)，往往是周邊區(qū)域環(huán)境污染嚴(yán)重，工業(yè)綠色發(fā)展下降時(shí)，會(huì)加大本單元政府調(diào)控力度，控制資源消耗，降低污染排放，由于政府調(diào)控具有滯后性，所以會(huì)促進(jìn)本單元下一期工業(yè)綠色發(fā)展。

根據(jù)公式(5)和(6)測(cè)算的長期效應(yīng)和短期效應(yīng)如表5所示。從總效應(yīng)來看，要素效率、產(chǎn)業(yè)效率和制度效率對(duì)京津冀工業(yè)綠色發(fā)展的促進(jìn)作用較顯著，且制度效率的促進(jìn)效應(yīng)最大，說明京津冀工業(yè)綠色發(fā)展對(duì)制度供給的依賴性較強(qiáng)；短期效應(yīng)明顯大于長期效應(yīng)，可見隨著時(shí)間的推移，促進(jìn)作用趨于下降，可能的原因在于本文所構(gòu)建的工業(yè)綠色發(fā)展指標(biāo)體系中多為速度型指標(biāo)，隨著供給側(cè)改革的持續(xù)推進(jìn)，工業(yè)綠色發(fā)展水平的提升速度減緩。從直接效應(yīng)來看，短期內(nèi)，要素效率和制度效率對(duì)本單元工業(yè)綠色發(fā)展的促進(jìn)作用較顯著，產(chǎn)業(yè)效率的促進(jìn)作用微弱，但產(chǎn)業(yè)效率和制度效率在長期會(huì)產(chǎn)生微弱的抑制作用，可見雖然短期內(nèi)政府調(diào)控的作用明顯，但長期內(nèi)必須發(fā)揮市場(chǎng)在資源配置中的主導(dǎo)作用。從溢出效應(yīng)來看，無論是短期還是長期，要素效率、產(chǎn)業(yè)效率和制度效率均具有正向溢出效應(yīng)，說明其他單元供給側(cè)效率的提升會(huì)促進(jìn)本單元工業(yè)綠色發(fā)展的改善，其中要素效率的長期溢出效應(yīng)小于短期溢出效應(yīng)，產(chǎn)業(yè)效率和制度效率的長期溢出效應(yīng)大于短期溢出效應(yīng)，可能的原因在于隨著京津冀協(xié)同發(fā)展戰(zhàn)略的深入推進(jìn)，產(chǎn)業(yè)協(xié)同、政府合作等不斷完善，但由于定位分工不同，對(duì)要素的需求有所差異。

表5 空間杜賓模型的空間直接效應(yīng)和空間溢出效應(yīng)

四、 結(jié)論及建議

京津冀工業(yè)綠色發(fā)展應(yīng)堅(jiān)持供給側(cè)結(jié)構(gòu)性改革為主線，供給側(cè)改革為工業(yè)綠色發(fā)展提供動(dòng)力。為探討供給側(cè)因素驅(qū)動(dòng)京津冀工業(yè)綠色發(fā)展的效應(yīng)，首先基于工業(yè)綠色發(fā)展的內(nèi)涵，參考四部委的《綠色發(fā)展指標(biāo)體系》，兼顧主觀與客觀，綜合運(yùn)用專家賦權(quán)法和粗糙集理論計(jì)算權(quán)重，并對(duì)2003—2017年京津冀工業(yè)綠色發(fā)展水平進(jìn)行測(cè)度；其次運(yùn)用探索性空間數(shù)據(jù)分析方法檢驗(yàn)京津冀間工業(yè)綠色發(fā)展的空間相關(guān)性；然后在搭建供給側(cè)因素分析框架的基礎(chǔ)上，考慮工業(yè)綠色發(fā)展的空間相關(guān)性、動(dòng)態(tài)性等，建立京津冀工業(yè)綠色發(fā)展與供給側(cè)因素的動(dòng)態(tài)空間杜賓模型，并將影響效應(yīng)分解為直接效應(yīng)和溢出效應(yīng)。結(jié)果發(fā)現(xiàn)：

(1)京、津、冀工業(yè)綠色發(fā)展均受環(huán)境質(zhì)量的影響最大；京、津、冀工業(yè)綠色發(fā)展和工業(yè)增長質(zhì)量總體上分別呈先下降后上升、下降、W型趨勢(shì)，工業(yè)資源利用總體上分別呈上升、下降、波動(dòng)頻繁趨勢(shì)，工業(yè)環(huán)境質(zhì)量總體上均呈上升趨勢(shì)，工業(yè)環(huán)境治理總體上分別呈先下降后上升、先上升后下降、先上升后下降趨勢(shì)。

(2)京津冀工業(yè)綠色發(fā)展具有顯著的空間相關(guān)性和動(dòng)態(tài)性特征；要素效率、產(chǎn)業(yè)效率、制度效率均對(duì)京津冀工業(yè)綠色發(fā)展產(chǎn)生不同程度地促進(jìn)作用，其中，對(duì)制度供給的依賴最大。

(3)從直接效應(yīng)來看，要素效率對(duì)京津冀工業(yè)綠色發(fā)展產(chǎn)生促進(jìn)作用；產(chǎn)業(yè)效率和制度效率在短期內(nèi)會(huì)產(chǎn)生促進(jìn)作用，但長期內(nèi)產(chǎn)生抑制作用。從溢出效應(yīng)來看，要素效率無論短期還是長期均具有正向溢出效應(yīng)；產(chǎn)業(yè)效率和制度效率在短期具有負(fù)向溢出效應(yīng)，在長期均轉(zhuǎn)為正向溢出效應(yīng)，但產(chǎn)業(yè)效率的效應(yīng)并不明顯。

因此，為進(jìn)一步促進(jìn)京津冀工業(yè)綠色發(fā)展，應(yīng)著力解決以下問題：

1.提高要素配置效率

強(qiáng)化產(chǎn)學(xué)研結(jié)合，基于產(chǎn)業(yè)需求加強(qiáng)人才培養(yǎng)，并實(shí)時(shí)發(fā)布高技術(shù)產(chǎn)業(yè)人才需求計(jì)劃，引導(dǎo)人力資本流向，提高行業(yè)間人力資本配置效率；加快落后產(chǎn)能淘汰和僵尸企業(yè)退出，引導(dǎo)投資向高技術(shù)產(chǎn)業(yè)和現(xiàn)代化生產(chǎn)性服務(wù)業(yè)傾斜；加大創(chuàng)新支持力度，營造良好創(chuàng)新環(huán)境，完善專利制度，鼓勵(lì)企業(yè)創(chuàng)新瞄準(zhǔn)節(jié)能減排技術(shù)領(lǐng)域，并促進(jìn)技術(shù)擴(kuò)散。

2.促進(jìn)產(chǎn)業(yè)協(xié)同和地區(qū)協(xié)同

一方面加強(qiáng)先進(jìn)制造業(yè)與現(xiàn)代化生產(chǎn)性服務(wù)業(yè)的協(xié)同發(fā)展，促進(jìn)工業(yè)綠色高質(zhì)量發(fā)展；另一方面加強(qiáng)京津冀工業(yè)綠色協(xié)同發(fā)展，在津冀建立高新技術(shù)產(chǎn)業(yè)園區(qū)，發(fā)揮聚集效應(yīng)，增強(qiáng)北京市生產(chǎn)性服務(wù)業(yè)的輻射能力，形成“核心—外圍區(qū)”結(jié)構(gòu)；但產(chǎn)業(yè)協(xié)同和地區(qū)協(xié)同實(shí)現(xiàn)的關(guān)鍵是建立合理有效的利益協(xié)調(diào)機(jī)制。

3.靈活協(xié)調(diào)市場(chǎng)與政府間關(guān)系

一方面綠色產(chǎn)業(yè)大多具有投入大、風(fēng)險(xiǎn)高、公益性程度高等特點(diǎn)，使得國有企業(yè)在早期工業(yè)綠色發(fā)展中扮演著主力軍的角色，但隨著工業(yè)綠色發(fā)展的不斷深入，市場(chǎng)化機(jī)制的引入十分必要；另一方面不同地區(qū)不同工業(yè)部門應(yīng)有針對(duì)性地綜合應(yīng)用命令控制型環(huán)境規(guī)制、經(jīng)濟(jì)激勵(lì)型環(huán)境規(guī)制與自愿意識(shí)型環(huán)境規(guī)制，為企業(yè)減負(fù)，釋放企業(yè)發(fā)展活力。