徐軍委，劉志華

（1. 中國勞動關(guān)系學(xué)院 勞動關(guān)系與人力資源學(xué)院，北京 100048；2. 北京林業(yè)大學(xué) 經(jīng)濟(jì)管理學(xué)院，北京 100083）

近年來，隨著經(jīng)濟(jì)的快速發(fā)展，資源環(huán)境的約束進(jìn)一步加劇，環(huán)境承載力與經(jīng)濟(jì)增長的矛盾日益凸顯，綠色發(fā)展理念也因此受到越來越多的關(guān)注，世界各國紛紛制定綠色發(fā)展的相關(guān)政策，以推動本國踐行綠色發(fā)展理念，走可持續(xù)發(fā)展之路。我國黨和政府歷來對綠色發(fā)展十分重視，自黨的十八屆五中全會首次確立綠色發(fā)展理念，到黨的十九大明確提出“必須樹立和踐行綠水青山就是金山銀山的理念”“形成綠色發(fā)展方式和生活方式”和“建設(shè)美麗中國”等，可見我國政府踐行綠色發(fā)展之路的決心。京津冀地區(qū)作為我國北方經(jīng)濟(jì)規(guī)模最大、發(fā)展活力最強(qiáng)的區(qū)域，戰(zhàn)略地位十分重要，大力推動京津冀地區(qū)間協(xié)同發(fā)展是未來大勢所趨。作為踐行綠色發(fā)展的關(guān)鍵區(qū)域，如何將“綠色發(fā)展”理念貫穿到京津冀地區(qū)協(xié)同發(fā)展的全過程，是一個亟待解決的關(guān)鍵問題。當(dāng)務(wù)之急是要科學(xué)評價京津冀地區(qū)各城市的綠色發(fā)展水平，明確差距，找出影響各地區(qū)綠色發(fā)展水平的關(guān)鍵因素。因此，構(gòu)建一套科學(xué)合理、可操作性強(qiáng)的綠色發(fā)展評價指標(biāo)體系，科學(xué)測度京津冀各城市綠色發(fā)展水平并分析其內(nèi)在影響機(jī)制，對于構(gòu)建京津冀綠色協(xié)同發(fā)展長效機(jī)制具有重要意義。

綠色發(fā)展已經(jīng)成為國內(nèi)外學(xué)者研究的熱門話題。從研究內(nèi)容來看，國內(nèi)外學(xué)者主要集中在綠色發(fā)展的內(nèi)涵界定、評價體系、水平測度及提升路徑等方面。從內(nèi)涵來看，綠色與發(fā)展之間的關(guān)系是綠色發(fā)展內(nèi)涵的重點，學(xué)者們最初研究綠色發(fā)展的重點在于綠色，隨著時代發(fā)展，后來逐漸聚焦于對生態(tài)環(huán)境與發(fā)展的關(guān)注，認(rèn)為綠色發(fā)展的內(nèi)涵是將生態(tài)環(huán)境作為核心生產(chǎn)力的一種發(fā)展模式[1-2]。此外，還有學(xué)者從綠色發(fā)展相關(guān)理論[3]、交互機(jī)制[4]、如何在實踐中應(yīng)用等方面[5]對其內(nèi)涵進(jìn)行了研究。目前對于綠色發(fā)展水平評價體系的研究成果最為豐富，但評價體系的側(cè)重點是能源消耗[6]、生態(tài)環(huán)境[7]還是宏觀經(jīng)濟(jì)[8]尚未形成統(tǒng)一定論。從水平測度來看，更多側(cè)重于基于不同方法通過構(gòu)建綠色發(fā)展評價指標(biāo)體系來對研究對象的綠色發(fā)展水平進(jìn)行測度。隨著研究的不斷深入，測度方法呈現(xiàn)多樣化，如熵值法[9]、數(shù)據(jù)包絡(luò)分析[10]、神經(jīng)網(wǎng)絡(luò)[11]、能值法[12]、網(wǎng)絡(luò)層次模型[13]、生態(tài)足跡法[14]等都在綠色發(fā)展水平的測度中得到了應(yīng)用。從提升路徑來看，有的學(xué)者傾向于借鑒國外經(jīng)驗來提升綠色發(fā)展水平[15-16]；部分學(xué)者認(rèn)為，重點在于提升科技發(fā)展水平和轉(zhuǎn)變經(jīng)濟(jì)增長方式[17-18]；還有學(xué)者則認(rèn)為，需要通過政府制定相關(guān)政策來提升綠色發(fā)展水平[19-20]。 從研究范圍來看，可以分為國家級、省級與區(qū)域級三方面，如張玉等[21]、張旭等[22]對中國省際綠色發(fā)展水平及影響因素進(jìn)行了研究，Guo 等[23]、孫永勝和佟連軍[24]對山東、吉林等地的綠色發(fā)展水平進(jìn)行了測度，還有學(xué)者重點關(guān)注了黃河流域[25]、長江經(jīng)濟(jì)帶[26]、珠江—西江經(jīng)濟(jì)帶[27]、京津冀[28]等地區(qū)的綠色發(fā)展問題。

綜上，國內(nèi)外對于綠色發(fā)展的內(nèi)涵、測度等研究日趨完善，但仍存在以下問題：一是從研究范圍來看，現(xiàn)有成果更多關(guān)注國家層面與省級層面，缺乏對京津冀地區(qū)綠色發(fā)展的成熟研究，尤其缺乏對京津冀地區(qū)城市層面的相關(guān)研究；二是從研究內(nèi)容來看，現(xiàn)有對京津冀地區(qū)城市綠色發(fā)展水平的相關(guān)研究多集中在水平測度與區(qū)域比較方面，缺乏對其原因的深層次探討；三是在研究綠色發(fā)展水平的影響因素時，忽略了區(qū)域間的空間關(guān)聯(lián)作用。因此，本文以京津冀地區(qū)13 個城市2002—2018 年的面板數(shù)據(jù)為基礎(chǔ)，構(gòu)建空間杜賓模型，分析其綠色發(fā)展水平的空間格局及影響因素，并提出提升京津冀地區(qū)綠色發(fā)展水平的相關(guān)建議，為京津冀地區(qū)綠色協(xié)同發(fā)展提供有益借鑒。

1 京津冀城市綠色發(fā)展水平測度與區(qū)域差異

1.1 京津冀城市綠色發(fā)展綜合評價指標(biāo)體系構(gòu)建

依據(jù)綠色發(fā)展的內(nèi)涵，堅持系統(tǒng)性、科學(xué)性、全面性的評價原則，結(jié)合國家相關(guān)政策及國內(nèi)外學(xué)者[29-32]對綠色發(fā)展評價指標(biāo)的研究成果，本文從綠色生產(chǎn)、綠色生活及綠色生態(tài)三方面構(gòu)建了京津冀城市綠色發(fā)展綜合評價指標(biāo)體系（表1）。指標(biāo)樣本來源于京津冀13 個城市2002—2018 年的相關(guān)統(tǒng)計數(shù)據(jù)，主要涉及《北京市國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計公報》《北京統(tǒng)計年鑒》《天津統(tǒng)計年鑒》《天津年鑒》以及《河北經(jīng)濟(jì)年鑒》等，其中缺失值采用回歸替換法進(jìn)行插補(bǔ)。

表1 京津冀城市綠色發(fā)展綜合評價指標(biāo)體系

1.2 京津冀城市綠色發(fā)展水平測度

由于綠色發(fā)展評價指標(biāo)分為正向指標(biāo)與負(fù)向指標(biāo)，需要先對數(shù)據(jù)進(jìn)行無量綱化處理，然后采用熵權(quán)TOPSIS 法測算13 個城市2002—2018 年的綠色發(fā)展指數(shù)（表2）。從整體增幅來看，京津冀地區(qū)的綠色發(fā)展水平呈逐年遞增的態(tài)勢，由2002 年的0.40 增長到2018 年的0.89，增幅達(dá)到157%。從城市層面來看，13 個城市的綠色發(fā)展水平也都得到了顯著提升。其中，增幅最大的是北京，達(dá)到335%，究其原因在于北京作為政治文化中心，無論是能源消費結(jié)構(gòu)調(diào)整、產(chǎn)業(yè)結(jié)構(gòu)升級還是生態(tài)環(huán)境保護(hù)都走在全國的前列；其次為承德、張家口與石家莊，增幅分別達(dá)到了323%、280%與203%；增幅最小的是邯鄲，僅為39%，主要原因在于該城市產(chǎn)業(yè)結(jié)構(gòu)極不合理，仍以重工業(yè)為主。從2018 年的綠色發(fā)展水平來看，北京的綠色發(fā)展水平最高，其次為天津與石家莊，這3 個城市的綠色發(fā)展水平均高于京津冀地區(qū)整體的綠色發(fā)展水平；其他10 個城市的綠色發(fā)展水平依次為廊坊＞秦皇島＞唐山＞張家口＞邢臺＞衡水＞承德＞保定＞邯鄲＞滄州，且這些地區(qū)的綠色發(fā)展水平均小于京津冀地區(qū)的整體綠色發(fā)展水平。表明京津冀城市綠色發(fā)展水平極不均衡，需要采取措施進(jìn)一步推動京津冀地區(qū)的綠色協(xié)同發(fā)展。

表2 2002—2018年京津冀城市綠色發(fā)展綜合指數(shù)

進(jìn)一步利用ArcGIS 分析京津冀城市的綠色生產(chǎn)、綠色生活與綠色生態(tài)子系統(tǒng)水平（圖1）。不難看出，2002—2018年13 個城市的三個子系統(tǒng)水平都有較大幅度提升，但總體分布特征沒有發(fā)生明顯變化。從綠色生產(chǎn)水平來看，2002—2018 年排名前三的城市一直為北京、天津、石家莊。從綠色生活水平來看，前三名由2002 年的北京、秦皇島、石家莊變?yōu)?018 年的北京、廊坊與石家莊。綠色生態(tài)排名前三的城市則由北京、天津、邢臺變?yōu)楸本?、天津、石家莊。此外，從起始年份來看，京津冀城市的綠色生活與綠色生態(tài)水平均低于綠色生產(chǎn)水平，因此，倡導(dǎo)綠色生活理念、加大生態(tài)環(huán)保力度是該地區(qū)未來綠色協(xié)同發(fā)展的重點。

2 京津冀城市綠色發(fā)展影響因素的模型設(shè)定

2.1 空間自相關(guān)檢驗

在構(gòu)建空間面板模型之前，首先要引入空間權(quán)重矩陣，通過全局空間自相關(guān)和局部空間自相關(guān)來判斷本地區(qū)的觀測值與相鄰地區(qū)觀測值之間的相關(guān)性是否顯著，只有二者顯著相關(guān)才需要構(gòu)建空間面板模型，否則只需要采用傳統(tǒng)的回歸模型即可。最常用的空間自相關(guān)檢驗包含全局空間自相關(guān)檢驗與局部空間自相關(guān)檢驗。

全局空間自相關(guān)檢驗結(jié)果用全局莫蘭指數(shù)（I）表示，具體公式如下：

圖1 2002—2018年京津冀城市綠色發(fā)展子系統(tǒng)水平

局部空間自相關(guān)檢驗結(jié)果用局部莫蘭指數(shù)（I′）表示，見公式（2）。當(dāng)I′＞0 時，表示研究對象在空間區(qū)域上存在正的空間自相關(guān)關(guān)系，且呈高—高聚集或低—低聚集；當(dāng)I′＜0 時，表示研究對象在空間區(qū)域上存在負(fù)的相關(guān)性，且呈高—低聚集或低—高聚集；當(dāng)I′=0 時，表示研究對象在空間上不相關(guān)。

2.2 空間面板計量模型

（1）空間杜賓模型（SDM）。不僅研究對象的自然經(jīng)濟(jì)社會因素會影響本地區(qū)的綠色發(fā)展水平，相鄰地區(qū)的相應(yīng)因素也會對本地區(qū)的綠色發(fā)展水平有重要影響，構(gòu)建相應(yīng)模型如下：

（2）空間滯后模型（SAR）。由于相鄰地區(qū)的綠色發(fā)展水平可能會影響本地區(qū)的綠色發(fā)展水平，因而構(gòu)建如下模型：

（3）空間誤差模型（SEM）。當(dāng)影響京津冀城市綠色發(fā)展水平的因素有遺漏時，需要構(gòu)建如下空間誤差模型：

式（3）～（5）中：Y 和X 分別為被解釋變量與解釋變量，W 為空間權(quán)重矩陣，λ 與θ 分別代表相鄰地區(qū)被解釋變量和解釋變量對本地區(qū)被解釋變量的影響程度，β 為解釋變量對被解釋變量的影響，ui代表個體效應(yīng)，δt代表空間效應(yīng)，εit為空間誤差項。

（4）空間權(quán)重矩陣。

鄰接空間權(quán)重矩陣，由于因變量不僅受到本地區(qū)自變量的影響，還受到相鄰地區(qū)自變量的影響，參考張軍偉等[33]的相關(guān)研究，將矩陣元素設(shè)置如下：

經(jīng)濟(jì)地理空間權(quán)重矩陣，參考伍國勇等[34]的研究，利用各地區(qū)GDP 均值與距離平方倒數(shù)的乘積構(gòu)建矩陣如下：

產(chǎn)業(yè)結(jié)構(gòu)空間權(quán)重矩陣，本文參考郝淑雙和朱喜安[35]的相關(guān)研究，首先分別計算城市i 和城市j 三次細(xì)分產(chǎn)業(yè)產(chǎn)值占GDP 的比重并將其構(gòu)建為兩個向量；然后，計算兩個城市產(chǎn)業(yè)比重向量的歐氏距離Dij，矩陣元素設(shè)置如下：

3 實證檢驗與結(jié)果分析

3.1 變量選擇

鑒于數(shù)據(jù)限制，本文選擇如下解釋變量：對外開放水平（fdi），用歷年實際利用外資占GDP 比重表示；經(jīng)濟(jì)發(fā)展水平（lngdp），用人均GDP 表示；產(chǎn)業(yè)結(jié)構(gòu)（is），用第三產(chǎn)業(yè)占GDP 比重表示；城鎮(zhèn)化（urban），用城鎮(zhèn)人口占總?cè)丝诒戎乇硎?；科技進(jìn)步（lntec），用三種專利授權(quán)數(shù)表示；環(huán)境規(guī)制（er），參照馮嚴(yán)超和王曉紅[36]的研究，用工業(yè)廢水處理量和工業(yè)廢水排放總量的比值來表示。為了緩解異方差帶來的影響，對人均GDP、三種專利授權(quán)數(shù)的絕對值取對數(shù)處理。各變量描述統(tǒng)計如表3 所示，不難看出無論是對外開放水平、經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)、城鎮(zhèn)化、科技進(jìn)步還是環(huán)境規(guī)制，京津冀地區(qū)各城市都存在顯著差異。

表3 變量描述統(tǒng)計表

3.2 空間自相關(guān)檢驗結(jié)果

由于城市發(fā)展模式相近的其產(chǎn)業(yè)結(jié)構(gòu)也類似，因此本文采用產(chǎn)業(yè)結(jié)構(gòu)空間權(quán)重矩陣，全局空間自相關(guān)檢驗結(jié)果如圖2 所示。從圖中可以看出，2002—2018 年的全局莫蘭指數(shù)均為正且通過10%顯著性水平的檢驗，表明京津冀城市綠色發(fā)展水平呈顯著的正相關(guān)。同時，從時間上來看，莫蘭指數(shù)呈波動上升趨勢，表明京津冀城市綠色發(fā)展水平的空間相關(guān)性在逐漸增強(qiáng)，本地區(qū)綠色發(fā)展水平受相鄰地區(qū)綠色發(fā)展水平的影響也在逐漸加深。

圖2 2002—2018年京津冀城市綠色發(fā)展水平全局莫蘭指數(shù)變動

為了更清晰地反映京津冀城市綠色發(fā)展水平的空間關(guān)聯(lián)，本文繪制了2002 與2018 年京津冀城市綠色發(fā)展水平的莫蘭散點圖（圖3）?？梢钥闯觯?002 年京津冀城市綠色發(fā)展水平主要分布在低—低聚集區(qū)，包括廊坊、張家口、承德、滄州、邢臺、石家莊6 個城市，2018 年增加了秦皇島變?yōu)? 個城市；高—高聚集區(qū)則一直為北京與天津兩市，沒有變化，其他城市分布在低—高聚集區(qū)域。對比起始年份來看，京津冀城市綠色發(fā)展水平主要分布在低—低聚集區(qū)，區(qū)域聚集情況比較穩(wěn)定，路徑依賴性較強(qiáng)。

圖3 2002年、2018年京津冀城市綠色發(fā)展水平莫蘭散點圖

3.3 實證結(jié)果分析

（1）傳統(tǒng)面板模型分析。在進(jìn)行空間面板回歸之前，首先進(jìn)行傳統(tǒng)面板回歸來判斷是否要進(jìn)行空間面板分析。結(jié)果（表4）顯示，空間滯后檢驗的兩個LM 檢驗其中之一拒絕了不存在空間依賴性的原假設(shè)，而空間誤差檢驗的兩個LM 檢驗均拒絕了不存在空間依賴性的原假設(shè)。與上一節(jié)空間自相關(guān)檢驗結(jié)果一致，進(jìn)一步證明了在實證分析時需要考慮空間因素的影響。

表4 傳統(tǒng)面板回歸分析結(jié)果

（2）空間面板模型分析。上一節(jié)已經(jīng)表明需要構(gòu)建空間面板模型，為了判定空間杜賓模型（SDM）是否可以簡化成空間滯后模型（SAR）或空間誤差模型（SEM），進(jìn)一步構(gòu)建Wald 與LR 統(tǒng)計量來進(jìn)行檢驗。表5 給出了三種模型的回歸結(jié)果，從擬合優(yōu)度來看，空間杜賓模型要優(yōu)于其他兩種模型，此外，Wald_spatial_lag、Wald_spatial_error、LR_spatial_lag 和LR_spatial_error 統(tǒng)計量均在不同顯著性水平下拒絕了原假設(shè)，因此，SDM 模型不能簡化為SEM 模型或SAR 模型。此外，Hausman 檢驗結(jié)果也表明拒絕隨機(jī)效應(yīng)的原假設(shè)。因此，應(yīng)使用固定效應(yīng)模型，本文采用時空雙固定空間杜賓模型進(jìn)行分析。

從表5 中可以看出：空間杜賓模型中空間自回歸系數(shù)在1%的水平下顯著為負(fù)，因此，京津冀地區(qū)某個城市綠色發(fā)展水平的提高會降低相鄰城市的綠色發(fā)展水平。從具體指標(biāo)來看，經(jīng)濟(jì)發(fā)展水平每提升1 個百分點，京津冀城市綠色發(fā)展水平下降0.111 7 個百分點，根源在于新常態(tài)下京津冀地區(qū)經(jīng)濟(jì)發(fā)展方式仍處于轉(zhuǎn)變期，結(jié)構(gòu)性、體制性、周期性問題相互交織，經(jīng)濟(jì)發(fā)展對生態(tài)環(huán)境造成了巨大壓力，從而抑制了地區(qū)綠色發(fā)展水平；產(chǎn)業(yè)結(jié)構(gòu)表現(xiàn)出的負(fù)向作用意味著隨著第三產(chǎn)業(yè)比重的增加，京津冀城市綠色發(fā)展水平不升反降，究其原因在于第三產(chǎn)業(yè)比重增加的減排效應(yīng)沒有抵消第二產(chǎn)業(yè)比重減少對技術(shù)進(jìn)步效應(yīng)的弱化，從而阻礙了綠色發(fā)展水平的提升；城鎮(zhèn)化對京津冀綠色發(fā)展水平有抑制作用，城鎮(zhèn)化水平每提升1 個百分點，綠色發(fā)展水平下降0.005 3個百分點，這歸因于京津冀城鎮(zhèn)化發(fā)展極不平衡，且城鎮(zhèn)化進(jìn)程中帶來了人口擁擠、產(chǎn)業(yè)失衡、環(huán)境污染等問題，這種負(fù)向作用超過了新型城鎮(zhèn)化帶來的正向作用；科技進(jìn)步對京津冀綠色發(fā)展水平具有促進(jìn)作用，隨著科技水平的不斷進(jìn)步，不僅能增強(qiáng)對污染物的治理水平，還可以充分利用先進(jìn)技術(shù)降低污染物的排放。對外開放水平與環(huán)境規(guī)制對綠色發(fā)展水平的影響不顯著。

表5 空間面板回歸分析結(jié)果

從各變量的空間滯后效應(yīng)來看，對外開放水平與環(huán)境規(guī)制的空間滯后效應(yīng)為負(fù)，城鎮(zhèn)化的空間滯后效應(yīng)為正，但都不顯著，表明其他地區(qū)對外開放水平、城鎮(zhèn)化和環(huán)境規(guī)制的提升不會對本地區(qū)綠色發(fā)展水平造成實質(zhì)影響；經(jīng)濟(jì)發(fā)展水平的空間滯后效應(yīng)顯著為負(fù)，表明隨著其他地區(qū)經(jīng)濟(jì)發(fā)展水平的提升，本地的綠色發(fā)展水平會惡化?？赡苁怯捎诘胤秸嬖诘牟┺母偁庩P(guān)系導(dǎo)致，如地方政府為了趕超其他地區(qū)的經(jīng)濟(jì)發(fā)展水平而忽視了生態(tài)環(huán)境保護(hù)等；產(chǎn)業(yè)結(jié)構(gòu)與科技進(jìn)步的空間滯后效應(yīng)顯著為正，表明其他地區(qū)產(chǎn)業(yè)結(jié)構(gòu)的優(yōu)化升級會提升本地區(qū)的綠色發(fā)展水平。原因在于地方政府在制定政策時存在明顯的相互“模仿”傾向[37]，當(dāng)其他地方政府制定產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化與科技進(jìn)步相關(guān)政策時，本地區(qū)地方政府部門通過“模仿”也會實現(xiàn)產(chǎn)業(yè)結(jié)構(gòu)、科技進(jìn)步和綠色發(fā)展水平的提升。

（3）空間溢出效應(yīng)分析。由于部分影響京津冀城市綠色發(fā)展水平的空間滯后系數(shù)不顯著為零，因此通過空間杜賓模型計算出的回歸系數(shù)不能直接反映其空間效應(yīng)。本文進(jìn)一步對其直接效應(yīng)、間接效應(yīng)和總效應(yīng)進(jìn)行了計算，結(jié)果見表6。其中，直接效應(yīng)代表上述影響因素對本地區(qū)綠色發(fā)展水平的影響，間接效應(yīng)代表各因素對周邊地區(qū)綠色發(fā)展水平的影響，總效應(yīng)則表示各因素對京津冀城市整體綠色發(fā)展水平的影響。

表6 空間杜賓模型的直接效應(yīng)、間接效應(yīng)與總效應(yīng)估計

經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)和科技進(jìn)步的總效應(yīng)顯著為正，是促進(jìn)京津冀城市綠色發(fā)展水平提升的主要驅(qū)動力；城鎮(zhèn)化的總效應(yīng)顯著為負(fù)，對京津冀城市綠色發(fā)展水平產(chǎn)生負(fù)向作用。經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)和科技進(jìn)步的直接效應(yīng)顯著為正，是促進(jìn)本地區(qū)綠色發(fā)展水平提升的主要因素，其中，經(jīng)濟(jì)發(fā)展水平的促進(jìn)作用最強(qiáng)，城鎮(zhèn)化是抑制本地區(qū)綠色發(fā)展水平提升的主要因素。產(chǎn)業(yè)結(jié)構(gòu)與城鎮(zhèn)化的間接效應(yīng)顯著為負(fù)，表明該影響因素具有負(fù)向空間溢出效應(yīng)；科技進(jìn)步的間接效應(yīng)顯著為正，表明該影響因素具有正向空間溢出效應(yīng)。

從經(jīng)濟(jì)發(fā)展水平來看，其直接效應(yīng)與總效應(yīng)顯著為正，表明經(jīng)濟(jì)發(fā)展水平的提升可以有效促進(jìn)京津冀城市綠色發(fā)展水平的提升，同時，經(jīng)濟(jì)發(fā)展水平的提高可以有效改善本地區(qū)的生產(chǎn)效率、環(huán)境治理和綠色發(fā)展水平等，但其間接效應(yīng)不顯著，不能有效提升其他地區(qū)的綠色發(fā)展水平。因此，加強(qiáng)京津冀城市間經(jīng)濟(jì)合作，促進(jìn)區(qū)域協(xié)同發(fā)展，是未來提升其綠色發(fā)展水平的重要手段。

從產(chǎn)業(yè)結(jié)構(gòu)來看，其直接效應(yīng)顯著為正。隨著產(chǎn)業(yè)結(jié)構(gòu)的優(yōu)化升級，低能耗、低污染、高附加值的知識密集型產(chǎn)業(yè)得到迅速發(fā)展，可以有效改善當(dāng)?shù)氐纳鷳B(tài)環(huán)境問題。同時，本地產(chǎn)業(yè)結(jié)構(gòu)的升級對其他城市的綠色發(fā)展水平產(chǎn)生了抑制作用，可能的解釋是本地區(qū)產(chǎn)業(yè)結(jié)構(gòu)升級必然會將重工業(yè)等產(chǎn)業(yè)轉(zhuǎn)移到其他城市，從而抑制了這些地區(qū)綠色發(fā)展水平的提升。需要注意的是，產(chǎn)業(yè)結(jié)構(gòu)升級雖然可以促進(jìn)京津冀城市整體綠色發(fā)展水平的提升，但這種效應(yīng)并不顯著。

從城鎮(zhèn)化來看，其各項效應(yīng)均顯著為負(fù)，表明隨著城鎮(zhèn)化水平的提升，不僅會抑制本地區(qū)的綠色發(fā)展水平，也會給其他城市與京津冀城市整體的綠色發(fā)展水平帶來抑制作用。究其原因，雖然京津冀地區(qū)一直在采取各種措施發(fā)展新型城鎮(zhèn)化，但各地方政府沒有處理好城鎮(zhèn)化與產(chǎn)業(yè)發(fā)展、基礎(chǔ)設(shè)施配套、生態(tài)環(huán)境保護(hù)等問題，從而致使城鎮(zhèn)化產(chǎn)生了負(fù)向空間溢出效應(yīng)。

從科技進(jìn)步來看，其直接效應(yīng)與間接效應(yīng)均顯著為正，表明科技進(jìn)步不僅促進(jìn)了本地區(qū)的綠色發(fā)展水平，通過空間溢出效應(yīng)也對京津冀其他城市整體綠色發(fā)展水平有顯著促進(jìn)作用。科技進(jìn)步可以有效促進(jìn)資源的循環(huán)利用，提高資源利用率，同時清潔生產(chǎn)技術(shù)可以提升生產(chǎn)效率與環(huán)境治理能力，為城市綠色發(fā)展提供動力。此外，科技進(jìn)步的總效應(yīng)不顯著，原因可能是京津冀城市間的研發(fā)機(jī)構(gòu)、科研成果等沒有得到有效共享，從而抑制了京津冀城市整體的綠色發(fā)展水平。

對外開放水平與環(huán)境規(guī)制效應(yīng)不顯著，這與學(xué)者方時姣和肖權(quán)[38]、張新林等[39]的研究結(jié)論相同。比較直觀的觀點是由于對外開放是把“雙刃劍”，雖然對外開放可以引進(jìn)外資與先進(jìn)技術(shù)、管理經(jīng)驗來促進(jìn)地區(qū)經(jīng)濟(jì)的發(fā)展，但不能忽視的是與此同時引進(jìn)了一些高污染、高能耗的企業(yè)，對生態(tài)環(huán)境造成了巨大破壞。環(huán)境規(guī)制的作用同樣具有兩面性，一方面，可以倒逼企業(yè)變革管理思想，革新技術(shù)降低環(huán)境污染；另一方面，由于監(jiān)管不力、政策落實不到位等使其不能發(fā)揮立竿見影的效果，從而致使其影響不顯著。

4 結(jié)論與建議

本文采用空間杜賓模型對2002—2018 年京津冀地區(qū)13 個城市的綠色發(fā)展水平進(jìn)行了測度，主要結(jié)論為：

（1）2002—2018 年，京津冀地區(qū)13 個城市的綠色發(fā)展水平逐年遞增，增幅最大的是北京，最小的是邯鄲。綠色生產(chǎn)、綠色生活與綠色生態(tài)子系統(tǒng)也均呈增長態(tài)勢，但總體分布態(tài)勢沒有發(fā)生較大改變，三系統(tǒng)排名靠前的城市均為北京、天津、石家莊、秦皇島和廊坊等地。

（2）京津冀城市綠色發(fā)展水平具有顯著的空間相關(guān)性，主要分布在低—低聚集區(qū)。說明綠色發(fā)展水平不僅受本地區(qū)經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)、城鎮(zhèn)化和科技進(jìn)步的影響，也受到其他城市相關(guān)因素的影響。

（3）通過構(gòu)建空間杜賓模型分析和空間溢出效應(yīng)分解不難發(fā)現(xiàn)，經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)和科技進(jìn)步是推動京津冀城市綠色發(fā)展水平提升的主要驅(qū)動力，城鎮(zhèn)化則抑制了地區(qū)綠色發(fā)展水平的提升。經(jīng)濟(jì)發(fā)展水平的直接效應(yīng)最高，其次為科技進(jìn)步與產(chǎn)業(yè)結(jié)構(gòu)，城鎮(zhèn)化是抑制本地區(qū)城市綠色發(fā)展水平提升的主要因素。科技進(jìn)步具有正向的空間溢出效應(yīng)，產(chǎn)業(yè)結(jié)構(gòu)與城鎮(zhèn)化具有負(fù)向的空間溢出效應(yīng)。對外開放水平與環(huán)境規(guī)制的效應(yīng)不顯著。

由于京津冀城市綠色發(fā)展水平具有空間溢出效應(yīng)，因此，在制定綠色發(fā)展水平提升策略時，不僅要合理規(guī)劃本地區(qū)的相關(guān)因素，還要考慮其他地區(qū)相關(guān)因素的影響，具體可以從以下方向進(jìn)行考慮：

（1）樹立綠色發(fā)展觀念，大力推動京津冀地區(qū)協(xié)同發(fā)展。北京的經(jīng)濟(jì)發(fā)展水平、科技創(chuàng)新、產(chǎn)業(yè)發(fā)展水平明顯高于其他城市，要充分發(fā)揮其帶動與輻射作用，深入實施京津冀協(xié)同發(fā)展戰(zhàn)略，加快構(gòu)建京津冀協(xié)同創(chuàng)新共同體，推動各項要素在京津冀城市間自由流動。

（2）推動新型城鎮(zhèn)化建設(shè)，打破區(qū)域發(fā)展不平衡。一方面，應(yīng)根據(jù)區(qū)域特色因地制宜發(fā)展特色產(chǎn)業(yè)，推動產(chǎn)業(yè)結(jié)構(gòu)高級化發(fā)展，充分發(fā)揮核心產(chǎn)業(yè)的外溢作用，穩(wěn)步推進(jìn)京津冀新型城鎮(zhèn)化建設(shè)；另一方面，要提升城市管理能力，保證城鎮(zhèn)化與周邊鄉(xiāng)村規(guī)劃在基礎(chǔ)設(shè)施網(wǎng)絡(luò)、產(chǎn)業(yè)布局、生態(tài)空間布局等方面的協(xié)調(diào)。

（3）加強(qiáng)區(qū)域合作，構(gòu)建京津冀協(xié)同治理機(jī)制與利益共享機(jī)制。通過建立健全京津冀城市間環(huán)境應(yīng)急聯(lián)動機(jī)制、跨區(qū)域生態(tài)補(bǔ)償機(jī)制和碳生態(tài)補(bǔ)償機(jī)制等，構(gòu)建京津冀協(xié)同治理機(jī)制，推動京津冀生態(tài)保護(hù)聯(lián)防聯(lián)控，加大對企業(yè)排污治理的監(jiān)管力度；同時，通過稅收、教育、醫(yī)療等方面的合作探索利益共享機(jī)制，充分發(fā)揮對外開放與環(huán)境規(guī)制對綠色發(fā)展水平提升的正向作用。