譚君印, 李益, 文傳浩，3

（1.蘭州大學(xué)經(jīng)濟(jì)學(xué)院，甘肅 蘭州 730000；2.重慶工商大學(xué)長江上游經(jīng)濟(jì)研究中心，重慶 南岸 400067；3.云南大學(xué)經(jīng)濟(jì)學(xué)院，云南 昆明 650091）

一、引言

改革開放四十多年來，我國創(chuàng)造了巨大的經(jīng)濟(jì)奇跡，也出現(xiàn)了一些深層次矛盾和問題，隨著資源稟賦發(fā)生改變，以及生產(chǎn)要素投入達(dá)到瓶頸，資源短缺與環(huán)境惡化對可持續(xù)發(fā)展的影響日益凸顯。 黨中央高度重視生態(tài)文明建設(shè)，鮮明提出了創(chuàng)新、協(xié)調(diào)、綠色、開放、共享的新發(fā)展理念，將綠色發(fā)展作為高質(zhì)量發(fā)展的生態(tài)本底。 重慶市地處長江上游和三峽庫區(qū)腹心地帶，在推動內(nèi)陸地區(qū)對外開放和帶動區(qū)域協(xié)調(diào)發(fā)展的大棋盤中，發(fā)揮著承東啟西、牽引南北、通達(dá)江海的重要作用［1］。從“構(gòu)建山清水秀美麗之地”，到“加快建設(shè)山清水秀美麗之地”，再到“推進(jìn)長江經(jīng)濟(jì)帶綠色發(fā)展中發(fā)揮示范作用”，習(xí)總書記對重慶綠色發(fā)展提出了更高的要求、更多的期盼。 近年來，重慶出臺了一系列生態(tài)文明建設(shè)規(guī)劃，提出將生態(tài)優(yōu)勢轉(zhuǎn)化為發(fā)展動能，把“綠色+”融入經(jīng)濟(jì)社會發(fā)展各方面。 各區(qū)縣也堅定不移地探索以生態(tài)優(yōu)先、綠色發(fā)展為導(dǎo)向的高質(zhì)量發(fā)展新路徑，綠色全要素生產(chǎn)率不斷提高，產(chǎn)業(yè)生態(tài)化和生態(tài)產(chǎn)業(yè)化進(jìn)程不斷加快。 正確認(rèn)識重慶綠色發(fā)展的現(xiàn)狀，洞悉各區(qū)縣綠色全要素生產(chǎn)率的變化差異和收斂趨勢，對黨和政府部署重慶綠色產(chǎn)業(yè)發(fā)展工作具有重要的現(xiàn)實意義。

全要素生產(chǎn)率是所有投入要素對產(chǎn)出增長貢獻(xiàn)的一種能力，常被用來衡量經(jīng)濟(jì)體的增長質(zhì)量和增長潛力。 Tinbergen 首次提出較為全面的能夠反映生產(chǎn)率的指標(biāo)——全要素生產(chǎn)率［2］，Davis 進(jìn)一步明確了全要素生產(chǎn)率的內(nèi)涵［3］，指出其測算過程應(yīng)該針對全部的要素投入而非某一種要素，包括在生產(chǎn)過程中使用的勞動、資本、原材料、能源等所有生產(chǎn)要素。 Solow 將增長核算方法與全要素生產(chǎn)率的測速結(jié)合在一起［4］，提出索洛余值法，后來學(xué)者對其進(jìn)行豐富和補充，不斷創(chuàng)新關(guān)于全要素生產(chǎn)率的理論和實踐［5］。隨著資源短缺與環(huán)境惡化對可持續(xù)發(fā)展的約束日益凸顯，越來越多的學(xué)者關(guān)注經(jīng)濟(jì)增長與資源環(huán)境之間的關(guān)系。 Pittman 將資源環(huán)境剛性影響共同納入社會福利與經(jīng)濟(jì)績效評價，提出用綠色全要素生產(chǎn)率（簡稱GTFP 指數(shù)）評估產(chǎn)業(yè)綠色發(fā)展現(xiàn)狀［6］，“波特假說”認(rèn)為環(huán)境規(guī)制會誘導(dǎo)技術(shù)變遷以提高生產(chǎn)效率并改善環(huán)境績效［7-8］，現(xiàn)實中獲取這種雙贏機(jī)會的潛力較大［9］，但也存在不容忽視的“污染避難所效應(yīng)”［10］和“回彈效應(yīng)”［11］。國內(nèi)學(xué)者關(guān)于綠色全要素生產(chǎn)率的早期研究主要是對國外研究成果的推廣性應(yīng)用研究，主要集中在綠色全要素生產(chǎn)率測算、演化差異［12］及其影響因素［13］等方面。隨著研究的逐漸深入，學(xué)者結(jié)合國家制度環(huán)境、政策實踐，在研究方法、研究視角、研究內(nèi)容等方面完善和拓展綠色全要素生產(chǎn)率體系框架。 一方面，學(xué)者從宏觀、中觀和微觀等視角，運用指 數(shù) 思 路［14］、參 數(shù) 思 路［15］和 非 參 數(shù) 思 路［16］等 方 法 測 度 旅 游 業(yè)［16］、智 能 制 造［17］、互 聯(lián) 網(wǎng)經(jīng)濟(jì)［18］、綠色金融［19］等特定部門或新興行業(yè)以及微觀企業(yè)［20］的綠色全要素生產(chǎn)率水平，探討不同維度綠色全要素生產(chǎn)率的時空演進(jìn)特征及收斂趨勢。另一方面，學(xué)者運用各種實證方法［20］，深入研究貿(mào)易開放［21］、環(huán)境規(guī)制［22］、智能化［17］、互聯(lián)網(wǎng)［18］、金融信貸［23］、新型城鎮(zhèn)化［24］等因素對綠色全要素生產(chǎn)率的影響，多視角驗證機(jī)制效應(yīng)的可靠性。

通過文獻(xiàn)梳理，發(fā)現(xiàn)國內(nèi)外學(xué)者對于全要素生產(chǎn)率理論基礎(chǔ)的研究比較完善，關(guān)于綠色全要素生產(chǎn)率也形成了較為豐富的研究體系， 不過研究視角多從省級層面或地級市層面展開，對經(jīng)濟(jì)區(qū)域或城鎮(zhèn)群的研究有待進(jìn)一步完善，關(guān)于重慶市區(qū)縣綠色全要素生產(chǎn)率的研究有待深入拓展?；诖耍疚囊灾貞c市38 個區(qū)縣為研究對象，運用ML 指數(shù)法測度重慶市綠色全要素生產(chǎn)率水平，并分析其增長變動、差異變動及收斂特征。 其研究結(jié)論將有助于揭示重慶市綠色發(fā)展現(xiàn)狀和差距，為重慶市發(fā)揮好綠色發(fā)展示范提供政策參考。

二、研究方法與數(shù)據(jù)處理

（一）研究方法

1.Malmquist-Luenberger生產(chǎn)率指數(shù)法

基于DEA 方法的Malmquist-Luenberger 指數(shù)（簡稱ML 指數(shù)）可用于計算隨時間變化的GTFP 指數(shù)變化速率，被廣泛應(yīng)用于各個領(lǐng)域［25］。 本文將每一個區(qū)縣視為一個獨立的決策單元，每個單元包括N 個輸入x=［x1+x2+…+xN］，x∈，M個預(yù)期輸出y=［y1+y2+…+yM］，y∈，以及I種非期望產(chǎn)出b=［b1+b2+…+bI］,y∈。 在此基礎(chǔ)上，建立t=1，2，…，T時期區(qū)縣k=1，2，…，T的一系列生產(chǎn)可能性函數(shù)：

基于此，預(yù)期產(chǎn)出和非期望產(chǎn)出的SBM 函數(shù)計算公式為：

基于SBM 函數(shù)，構(gòu)建從t到t+1 期的綠色全要素生產(chǎn)率指數(shù)（簡稱GTFP指數(shù)，下同）：

上式說明， 可以用ML指數(shù)代表GTFP指數(shù)變化， 對于一個區(qū)縣來說，ML的測算值大于（小于）1，分別表示在t到t+1 期間GTFP指數(shù)增多（減少）。

2.核密度估計法

核密度估計利用連續(xù)的密度曲線描述隨機(jī)變量的分布形態(tài)，反映其要素在空間分布中的距離衰減效應(yīng)［26］。 本文選用核密度估計法模擬重慶市GTFP指數(shù)的集散特征，公式表示為：

其中，n是觀測值數(shù)，GTFPit是i區(qū)縣在第t年份的綠色全要素生產(chǎn)率指數(shù)，GTFPt為重慶第t年份的均值，k（GTFPit）為高斯核函數(shù)，h為窗寬。 為減少窗寬h對核密度估計結(jié)果的影響，更好地呈現(xiàn)核函數(shù)的典型分布特征，本文將窗寬設(shè)置為1 年。

3.泰爾指數(shù)

泰爾指數(shù)作為衡量某地域各種社會發(fā)展差異性的重要工具被廣泛應(yīng)用，具備良好的可分解性質(zhì)，其數(shù)值越低，說明該區(qū)域社會發(fā)展的差異性越?。?7］。 本文使用泰爾指數(shù)分析重慶市GTFP指數(shù)組內(nèi)外差距。 根據(jù)重慶市“一區(qū)兩群”戰(zhàn)略布局，將重慶38 個區(qū)縣分為主城都市區(qū)、渝東南武陵山區(qū)城鎮(zhèn)群、渝東北三峽庫區(qū)城鎮(zhèn)群3 個群組，每個群組用gk（k=1，2，3）表示，nk代表第k群組的區(qū)縣數(shù)，有nk=38；用gtfpit表示第t年份i區(qū)縣GTFP指數(shù)占總體的份額，用gtfpkt表示第t年份第k組的GTFP指數(shù)占總體的份額， 記Tbt與Twt分別為組間差距和組內(nèi)差距，則可將泰爾指數(shù)分解為：

公式中將Tkt=gtfpkt×，Tkt設(shè)為第t年份第k組組內(nèi)差距。進(jìn)一步可以計算第t年份第k組組內(nèi)差距貢獻(xiàn)率和組間差距的貢獻(xiàn)率，分別用Dkt和Dbt表示：

4.收斂性分析

收斂理論最初用于考察區(qū)際人均收入是否存在收斂或發(fā)散的問題， 近年來的應(yīng)用范圍逐漸擴(kuò)展到產(chǎn)業(yè)效率方面。 在本文中，收斂性分析有助于檢驗重慶市區(qū)際GTFP指數(shù)差距是在擴(kuò)大還是在縮小，以及未來是否會趨于平衡發(fā)展［16］。 收斂性分析包括σ 收斂、β 收斂和俱樂部收斂。

（1）綠色全要素生產(chǎn)率σ 收斂

σ 收斂反映的是重慶市綠色全要素生產(chǎn)率隨時間演變的波動趨勢， 常用的衡量樣本偏離程度的指標(biāo)包括標(biāo)準(zhǔn)差、變異系數(shù)、泰爾指數(shù)、基尼系數(shù)。 其中，標(biāo)準(zhǔn)差是反映各觀測值離散程度的一個絕對指標(biāo)，但容易受隨機(jī)變量量綱影響，而變異系數(shù)是無量綱的，可以消除平均值差異對樣本的干擾，在平均值差別較大時，變異系數(shù)可以近似地評價樣本分散度［28］。 因此，本文采用標(biāo)準(zhǔn)差（σt）和變異系數(shù)（cvt）共同描述重慶區(qū)際GTFP指數(shù)的σ 收斂狀態(tài)，公式表示為：

（2）綠色全要素生產(chǎn)率β 收斂

β 收斂反映的是隨時間推移，落后地區(qū)能夠以較快的增速趕上發(fā)達(dá)地區(qū)，最后保持相同水平的收斂狀態(tài)［29］。β 收斂又分為絕對β 收斂和條件β 收斂，區(qū)別在于是否嚴(yán)格控制一系列影響因素。鑒于經(jīng)典β 收斂側(cè)重于從時間上考察重慶GTFP指數(shù)收斂趨勢，如果不考慮空間因素，可能會導(dǎo)致回歸結(jié)果的有偏性［30］。 因而，本文在經(jīng)典β 收斂分析的基礎(chǔ)上，引入空間因素，檢驗重慶市GTFP指數(shù)收斂的空間β 收斂特征。不過在進(jìn)行空間β 收斂回歸前，需要進(jìn)行全局空間相關(guān)性檢驗，討論GTFP的聚集程度和空間溢出效應(yīng)。 公式表示為：

公式中分母可構(gòu)成樣本方差S2=；w1為空間鄰接矩陣。 具體的模型描述如下。

空間β 收斂：由于地區(qū)間普遍存在不同程度的空間依賴性，因此本文建立了空間杜賓模型（Sdm）、空間誤差模型（Sem）和空間滯后模型（Sar）三種空間計量模型，根據(jù)相關(guān)檢驗判斷GTFP空間β 收斂最優(yōu)模型，分別表示為：

其中，GTFPit+1/GTFPit為i區(qū)縣在第t+1 年份的綠色全要素生產(chǎn)率增長率，GXit為第s個控制變量的估計系數(shù)，其中s=1，2，3，4。 當(dāng)?θs=0 時，上述模型為絕對收斂GXit=GTFPit，反之則為條件收斂，GXit為s個控制變量與GTFPit構(gòu)成的集合。 β 為收斂系數(shù)，當(dāng)β＜0 時，則GTFP指數(shù)趨于收斂，反之則趨于擴(kuò)散。w1為空間權(quán)重矩陣，考慮到區(qū)際協(xié)同創(chuàng)新發(fā)展突破了距離的限制，是經(jīng)濟(jì)轉(zhuǎn)型升級和綠色發(fā)展的支點力量，因此本文采取經(jīng)濟(jì)矩陣，用區(qū)際人均GDP年絕對差額的倒數(shù)表征。

（3）綠色全要素生產(chǎn)率俱樂部收斂

具有相似特征的地區(qū)間趨于收斂，而不具有相似特征的地區(qū)間不存在收斂，這種地區(qū)間收斂被稱為俱樂部收斂［26］。 本文結(jié)合綠色全要素生產(chǎn)率特征與“一區(qū)兩群”協(xié)同戰(zhàn)略布局，選擇對“一區(qū)兩群”GTFP指數(shù)進(jìn)行俱樂部收斂檢驗。

（二）指標(biāo)選取

1.核心解釋變量——綠色全要素生產(chǎn)率（GTFP）

結(jié)合已有文獻(xiàn)在投入產(chǎn)出指標(biāo)上綜合考量多種因素， 基于重慶市區(qū)縣數(shù)據(jù)的可獲取性，在投入指標(biāo)上，考慮了勞動投入、資源要素投入、資本投入、技術(shù)投入四個方面。 資本投入選擇全社會固定資產(chǎn)投資（百億元）表征；資源要素投入包括能源要素和土地要素，受限于數(shù)據(jù)，土地資源投入指標(biāo)采用非住宅類房屋建筑竣工面積（萬平方米）來衡量，能源投入指標(biāo)選擇工業(yè)企業(yè)能源消費總量（萬噸標(biāo)準(zhǔn)煤）表征；技術(shù)投入選用發(fā)明專利授權(quán)量；勞動投入使用城鎮(zhèn)非私營單位職工人數(shù)表征。 用在產(chǎn)出指標(biāo)上，期望產(chǎn)出以GDP（百億元）代表當(dāng)前的產(chǎn)出能力，非期望產(chǎn)出采用PM2.5濃度（μg/m3）和CO2排放量（百萬噸）表示。

2.控制變量

（1）基礎(chǔ)設(shè)施（Inds）：采用公路里程數(shù)與轄區(qū)面積的比重（km/km2）來進(jìn)行表征；（2）數(shù)字金融（Fin）：采用北京大學(xué)數(shù)字金融研究中心公布的縣級數(shù)字普惠金融發(fā)展指數(shù)表征；（3）市場化（Mar）： 采用國有經(jīng)濟(jì)單位職工人數(shù)和集體經(jīng)濟(jì)單位職工人數(shù)占城鎮(zhèn)非私營單位在崗職工人數(shù)表征；（4）城鎮(zhèn)化（Town）：采取城鎮(zhèn)常住人口占總?cè)丝诘陌俜直葋磉M(jìn)行衡量。

（三）數(shù)據(jù)來源

本文選取2011—2019 年重慶市38 個區(qū)縣的面板數(shù)據(jù)，所有數(shù)據(jù)來源于《重慶市統(tǒng)計年鑒（2012—2020）》、重慶市知識產(chǎn)權(quán)局（http://zscqj.cq.gov.cn/）、重慶市生態(tài)環(huán)境局（http://sthjj.cq.gov.cn/）、國家地球物理數(shù)據(jù)中心（https://www.ngdc.noaa.gov/）和北京大學(xué)數(shù)字金融研究中心（https://idf.pku.edu.cn/），全社會固定資產(chǎn)投資額用各區(qū)縣2011 年為基期的固定資產(chǎn)投資價格指數(shù)平減，GDP以各區(qū)縣2011 年為基期的GDP平減指數(shù)進(jìn)行平減， 數(shù)字普惠金融指數(shù)均除以100 表征，部分缺失值使用插值法補充。 表1 列出了各變量描述性統(tǒng)計結(jié)果。

表1 各變量描述性統(tǒng)計結(jié)果

三、研究過程與結(jié)果分析

（一）綠色全要素生產(chǎn)率水平測算

本文測算出2012—2019 年重慶市GTFP指數(shù)水平，并整理得到其空間分布圖（見圖1）。從圖1 顯示的空間分布情況來看， 總體上重慶市產(chǎn)業(yè)生態(tài)化、 生態(tài)產(chǎn)業(yè)化發(fā)展態(tài)勢良好，其GTFP指數(shù)呈現(xiàn)波動上升趨勢，其中，2012—2017 年重慶區(qū)際GTFP指數(shù)呈現(xiàn)高速增長，年均增速高達(dá)5.09%，2017 年后增速逐漸放緩，年均增速僅1.90%。 從“一區(qū)兩群”協(xié)調(diào)發(fā)展戰(zhàn)略布局來看，渝東南武陵山區(qū)城鎮(zhèn)群后發(fā)趕超特征明顯，渝東南武陵山區(qū)城鎮(zhèn)群2012 年GTFP均值（0.93）略高于重慶市平均水平（0.92），到2017 年達(dá)到1.48，成功趕超主城都市區(qū)和渝東北三峽庫區(qū)城鎮(zhèn)群。 主城都市區(qū)和渝東北三峽庫區(qū)城鎮(zhèn)群則呈現(xiàn)交替上升現(xiàn)象。 2013—2017 年渝東北三峽庫區(qū)城鎮(zhèn)群GTFP指數(shù)均值高于主城都市區(qū)， 而2018—2019 年和2012 年期間主城都市區(qū)高于渝東北三峽庫區(qū)城鎮(zhèn)群，不過兩個地區(qū)GTFP指數(shù)差值較為穩(wěn)定，絕對值始終低于0.09。從代表性區(qū)縣來看，豐都、梁平、綦江、秀山和南川5 個區(qū)縣GTFP指數(shù)2012—2019 年均增速都超過10%，綜合排名也大幅上升，如綦江區(qū)2012 年GTFP指數(shù)排名第33 位，到2019年上升到第5 位，長壽、大足、九龍坡、萬州、石柱和武隆6 個區(qū)縣GTFP指數(shù)呈現(xiàn)負(fù)增長，綜合排名大幅下滑，6 年間下滑了14～24 名不等， 如武隆縣GTFP指數(shù)排名從2012 年的第4 位下滑到2019 年的第27 位，萬州區(qū)GTFP指數(shù)排名從2012 年的第13 位下滑到2019 年的第36 位。

圖1 重慶市GTFP 指數(shù)的空間分布

（二）分布動態(tài)

1.時間差異

根據(jù)重慶市GTFP指數(shù)繪制出其核密度曲線，如圖2 所示。 重慶市整體GTFP指數(shù)核密度曲線在波動中緩慢右移，說明重慶市整體GTFP指數(shù)普遍是有效的，具有上升趨勢。 從各年份分布形態(tài)來看，2012—2016 年曲線由寬峰分布逐漸向尖峰分布轉(zhuǎn)變， 主體面積位于橫坐標(biāo)1的左側(cè)，波峰高度不斷上升，右擺尾逐漸縮小，表明在2012—2016 年重慶市整體GTFP指數(shù)呈現(xiàn)低水平集聚趨勢。 2016 年后，曲線由尖峰分布逐漸向?qū)挿宸植嫁D(zhuǎn)變，主體面積位于橫坐標(biāo)1的右側(cè)，波峰高度不斷下降，右擺尾明顯拉長，表明在2016 年以后重慶市GTFP指數(shù)呈現(xiàn)中高水平離散趨勢，說明了整個區(qū)域內(nèi)GTFP指數(shù)區(qū)際不協(xié)調(diào)性正在擴(kuò)大。從“一區(qū)兩群”協(xié)調(diào)發(fā)展來看，三大區(qū)域和重慶市整體保持一致，但也呈現(xiàn)出部分差異。 主城都市區(qū)GTFP指數(shù)核密度曲線在2016 年以后， 右擺尾面積不斷擴(kuò)大， 存在多峰極化趨勢， 顯示主城都市區(qū)內(nèi)部區(qū)縣GTFP指數(shù)在后期中高水平離散分布更加明顯。 渝東北三峽庫區(qū)城鎮(zhèn)群GTFP指數(shù)曲線在2012—2015 年期間，由單寬峰分布逐漸向多尖峰分布轉(zhuǎn)變，尤其是在2013 年出現(xiàn)四個尖峰分布特征，左右擺尾均明顯縮短，說明渝東北三峽庫區(qū)城鎮(zhèn)群GTFP指數(shù)在早期處于低水平離散分布特征。渝東南武陵山區(qū)城鎮(zhèn)群GTFP指數(shù)核密度曲線波峰數(shù)量和高度在2016 年前明顯波動，2016 年以后曲線右移明顯加快，波峰也逐漸消失，顯示出渝東南武陵山區(qū)城鎮(zhèn)群GTFP指數(shù)在早期提升具有不穩(wěn)定性，但在2016 年后逐漸朝著高水平均衡方向發(fā)展。

圖2 重慶市GTFP 指數(shù)的核密度曲線

2.空間差異

根據(jù)重慶市GTFP指數(shù)計算得到泰爾指數(shù)，如表2 所示。 可以發(fā)現(xiàn)，2012—2019 年重慶市GTFP指數(shù)總差距呈現(xiàn)U 型分布，2012—2016 年總體差距在波動中逐漸縮小，之后逐年擴(kuò)大，2019 年總差距已經(jīng)超過2012 年總體差距， 組間差距呈現(xiàn)先擴(kuò)大后縮小的態(tài)勢，2014 年和2017 年組間差距急劇增加，組內(nèi)差距則表現(xiàn)為先減小后增加趨勢，這是總體差距的主要來源。從“一區(qū)兩群”組間差距來看，主城都市區(qū)GTFP指數(shù)的組內(nèi)差距呈下降趨勢，渝東北三峽庫區(qū)城鎮(zhèn)群GTFP指數(shù)的組內(nèi)差距則基本無變化， 且處于低水平均衡， 渝東南武陵山區(qū)城鎮(zhèn)群GTFP指數(shù)的組內(nèi)差距在波動中擴(kuò)大。 “一區(qū)兩群”組內(nèi)差距貢獻(xiàn)率均值分布順序為：主城都市圈（0.525）＞渝東南武陵山區(qū)城鎮(zhèn)群（0.225）＞渝東北三峽庫區(qū)城鎮(zhèn)群（0.190）。

表2 重慶市GTFP 指數(shù)泰爾指數(shù)

（三）收斂性檢驗

1.σ 收斂分析

根據(jù)重慶市GTFP指數(shù)可以繪制出其標(biāo)準(zhǔn)差和變異系數(shù)隨年份的變化趨勢，如圖3 所示。從中可以發(fā)現(xiàn)，重慶市總體和“一區(qū)兩群”GTFP指數(shù)標(biāo)準(zhǔn)差和變異系數(shù)的變化趨勢基本一致，只是渝東北三峽庫區(qū)城鎮(zhèn)群和渝東南武陵山區(qū)城鎮(zhèn)群的波動性更明顯。2012—2016 年綠色全要素生產(chǎn)率標(biāo)準(zhǔn)差和變異系數(shù)呈微弱下降趨勢， 表現(xiàn)出明顯的σ 收斂趨勢，2016 年后標(biāo)準(zhǔn)差和變異系數(shù)快速上升，說明不存在顯著的σ 收斂。 此外，從上文GTFP指數(shù)泰爾指數(shù)重慶市總體差距和“一區(qū)兩群”組內(nèi)差距分析結(jié)果也充分證明了該結(jié)論。

圖3 重慶市GTFP 指數(shù)標(biāo)準(zhǔn)差與變異系數(shù)

2.β 收斂分析

首先，在進(jìn)行空間β 收斂檢驗時，需要采用Moran’s I指數(shù)方法對綠色全要素生產(chǎn)率的空間相關(guān)性進(jìn)行檢驗。 表3 所示的檢驗數(shù)據(jù)顯示，2012—2019 年重慶市綠色全要素生產(chǎn)率Moran’s I指數(shù)均顯著為負(fù)或為正，說明重慶市鄰近區(qū)縣之間綠色全要素生產(chǎn)率的影響存在空間依賴性［31］。

表3 重慶市綠色全要素生產(chǎn)率Moran’s I 指數(shù)

此外，空間收斂模型在進(jìn)行空間β 收斂過程中需要確定并識別空間杜賓模型（Sdm）是否會退化為空間誤差（Sem）模型或者空間滯后（Sar）模型。 如表4 所示的Lr 檢驗、Lm 檢驗和Wald 檢驗可以發(fā)現(xiàn)空間杜賓模型將退化成空間誤差模型。 對于隨機(jī)效應(yīng)和固定效應(yīng)的選擇，通過表5 和表6 中的hausman檢驗結(jié)果，可以發(fā)現(xiàn)前述10 個模型均接受原假設(shè)，選擇隨機(jī)效應(yīng)更為合理。 綜上所述，本文選擇隨機(jī)效應(yīng)的空間誤差模型分析綠色全要素生產(chǎn)率空間收斂性。

表4 空間回歸檢驗結(jié)果

表5 顯示了重慶市GTFP指數(shù)β 收斂的回歸結(jié)果，回歸結(jié)果表現(xiàn)為四個方面。 （1）重慶市GTFP指數(shù)存在明顯收斂趨勢。GTFP指數(shù)系數(shù)均在1%顯著性水平下為負(fù)，說明GTFP指數(shù)與其增長率之間是反向關(guān)系，收斂趨勢顯示重慶市GTFP指數(shù)的差異在縮小，低GTFP指數(shù)地區(qū)存在后發(fā)趕超的特點。 （2）空間因素具有減速效應(yīng)。 加入空間因素后，絕對收斂與條件收斂的回歸系數(shù)分別縮小了17.15%和15.00%，可能的原因在于，由于空間效應(yīng)的存在，GTFP指數(shù)不均衡分布使得經(jīng)濟(jì)水平相近地區(qū)之間的相互影響更加顯著，特別是區(qū)際生產(chǎn)要素和信息的空間交流與互動，進(jìn)一步加劇了GTFP指數(shù)空間溢出與擴(kuò)散效應(yīng)，從而使得地區(qū)間GTFP指數(shù)空間差異隨時間推移呈現(xiàn)出進(jìn)一步放大的趨勢，減緩了收斂速度，進(jìn)而拉長了收斂周期。 （3）條件收斂的收斂速度均顯著大于絕對收斂速度。 加入控制變量后，經(jīng)典和空間模型下的收斂速度提高了5.62%和8.34%，可能是因為區(qū)際生產(chǎn)條件的異質(zhì)性，控制變量加速了GTFP指數(shù)收斂，同時壓縮了收斂周期，進(jìn)而更加準(zhǔn)確可靠地反映模型的收斂特征。 （4）加入控制變量后，條件模型2 中Fin、Inds和Mar系數(shù)顯著大于零，條件模型4 中Fin系數(shù)顯著為正，數(shù)字金融、基礎(chǔ)設(shè)施覆蓋率和市場化促進(jìn)了綠色全要素生產(chǎn)率提升，數(shù)字金融通過解決信息不對稱和降低資本的使用成本等問題在一定程度上解決了融資約束問題，從而提高綠色技術(shù)升級型項目投資成功的概率，并最終作用于GTFP指數(shù)；市場化能夠改善勞動力市場要素扭曲程度和資本市場要素扭曲程度的作用路徑， 提高重慶市GTFP指數(shù)； 良好的基礎(chǔ)設(shè)施能夠促進(jìn)物質(zhì)文化交流，城市規(guī)劃和土地利用更加合理，從而形成產(chǎn)業(yè)集聚效應(yīng)，大幅提升了企業(yè)生產(chǎn)效率和全社會運行效率。 加入空間因素后，Town系數(shù)顯著為負(fù)， 表明城鎮(zhèn)化抑制綠色全要素生產(chǎn)率的提升，各區(qū)縣城鎮(zhèn)人口容量能力有限，大規(guī)模農(nóng)業(yè)人口向城市轉(zhuǎn)移，帶來生產(chǎn)要素流動不均衡，加大了區(qū)縣生態(tài)環(huán)境負(fù)擔(dān)。

表5 β 收斂回歸結(jié)果

3.俱樂部收斂分析

表6 顯示了重慶市GTFP指數(shù)俱樂部收斂的回歸結(jié)果，從中可以發(fā)現(xiàn)“一區(qū)兩群”GTFP指數(shù)條件收斂系數(shù)都通過了顯著性檢驗，且小于0，表明“一區(qū)兩群”GTFP指數(shù)都存在俱樂部收斂，這與重慶市GTFP指數(shù)β 收斂趨勢一致。 隨著城鎮(zhèn)化、基礎(chǔ)設(shè)施、數(shù)字金融、市場化等因素調(diào)整，“一區(qū)兩群”GTFP指數(shù)均將達(dá)到各自的收斂水平。 從系數(shù)絕對值來看，不論是經(jīng)典模型還是空間模型，主城都市區(qū)收斂速度最快，渝東北三峽庫區(qū)城鎮(zhèn)群次之，而渝東南武陵山區(qū)城鎮(zhèn)群收斂速度低于重慶市β 收斂速度，顯示出主城都市區(qū)、渝東北三峽庫區(qū)城鎮(zhèn)群和渝東南武陵山區(qū)城鎮(zhèn)群的GTFP指數(shù)俱樂部收斂依次遞減的分布格局。從空間減速效應(yīng)來看，空間因素對渝東南武陵山區(qū)城鎮(zhèn)群的減速效應(yīng)最強，達(dá)到22.28%，而主城都市區(qū)和渝東北三峽庫區(qū)城鎮(zhèn)群空間減速效應(yīng)相差不大，分別是17.59%和17.27%。 這可能是因為渝東南武陵山區(qū)城鎮(zhèn)群文旅融合發(fā)展示范區(qū)建設(shè)帶動其生態(tài)康養(yǎng)、文化旅游產(chǎn)業(yè)高質(zhì)量協(xié)同發(fā)展，空間溢出提升了區(qū)際GFTP指數(shù)水平，形成對其他區(qū)縣趕超之勢。 從控制變量來看，數(shù)字金融對渝東南武陵山區(qū)城鎮(zhèn)群GTFP指數(shù)促進(jìn)效應(yīng)最強，主城都市區(qū)次之，渝東北三峽庫區(qū)城鎮(zhèn)群最弱；基礎(chǔ)設(shè)施主要作用于對渝東南武陵山區(qū)城鎮(zhèn)群GTFP指數(shù)， 市場化促進(jìn)主城都市圈和渝東南武陵山區(qū)城鎮(zhèn)群GTFP指數(shù)提升，而城鎮(zhèn)化在空間效應(yīng)作用下，抑制主城都市區(qū)和渝東南武陵山區(qū)城鎮(zhèn)群GTFP指數(shù)提升。

表6 俱樂部收斂回歸結(jié)果

四、結(jié)論與啟示

本文測度了2012—2019 年重慶市區(qū)際GTFP指數(shù)水平， 在討論其時空分布動態(tài)的基礎(chǔ)上，將空間因素引入經(jīng)典收斂模型，分析重慶市及其“一區(qū)兩群”GTFP指數(shù)的σ 收斂、β 收斂和俱樂部收斂特征，得到以下研究結(jié)論。

從時空演進(jìn)來看，重慶市GTFP指數(shù)呈現(xiàn)不斷上升趨勢。 其中，重慶市GTFP在早期呈現(xiàn)低水平集聚趨勢，區(qū)際差距逐漸縮小，2016 年后呈現(xiàn)中高水平離散趨勢，區(qū)際差距逐漸擴(kuò)大。“一區(qū)兩群”與重慶市整體保持一致，但也呈現(xiàn)出部分差異，渝東南武陵山區(qū)城鎮(zhèn)群GTFP指數(shù)早期提升具有不穩(wěn)定性，2016 年后逐漸朝著高水平均衡方向發(fā)展， 主城都市區(qū)在后期高水平離散分布更加明顯，渝東北三峽庫區(qū)城鎮(zhèn)群GTFP指數(shù)在早期呈現(xiàn)低水平多極化特征。 豐都、梁平、綦江、秀山和南川5 個區(qū)縣GTFP指數(shù)年均增速和綜合排序大幅上升，而長壽、大足、九龍坡、萬州、石柱和武隆6 個區(qū)縣GTFP指數(shù)呈現(xiàn)負(fù)增長，綜合排名也大幅下滑。

從收斂性來看，重慶GTFP指數(shù)收斂速度正趨于一致。 2012—2016 年綠色全要素生產(chǎn)率標(biāo)準(zhǔn)差和變異系數(shù)呈微弱下降趨勢， 表現(xiàn)出σ 收斂特征，2016 年后標(biāo)準(zhǔn)差和變異系數(shù)快速上升，σ 收斂逐漸消失。 重慶市總體GTFP指數(shù)隨時間推移呈現(xiàn)逐漸縮小的β 收斂趨勢，區(qū)域差異擴(kuò)大的現(xiàn)象將得到遏制， 主城都市區(qū)、 渝東北三峽庫區(qū)城鎮(zhèn)群和渝東南武陵山區(qū)城鎮(zhèn)群GTFP指數(shù)收斂速度呈現(xiàn)依次遞減的分布格局。 空間因素對GTFP指數(shù)收斂具有減速效應(yīng)，并隨時間不斷增強，尤其是對渝東南武陵山區(qū)城鎮(zhèn)群的減速效應(yīng)最明顯。 數(shù)字金融、基礎(chǔ)設(shè)施覆蓋率和市場化促進(jìn)重慶市總體GTFP指數(shù)提升，城鎮(zhèn)化表現(xiàn)為空間抑制作用，四個控制變量對渝東南武陵山區(qū)城鎮(zhèn)群的影響明顯，主城都市區(qū)次之，渝東北三峽庫區(qū)城鎮(zhèn)群GTFP指數(shù)最弱。

本文的研究給我們兩個方面的啟示。（1）立足當(dāng)前經(jīng)濟(jì)，推動產(chǎn)業(yè)生態(tài)化。一是存量升級，加強信息基礎(chǔ)設(shè)施建設(shè)，提升大數(shù)據(jù)在綠色產(chǎn)品加工全產(chǎn)業(yè)鏈中的應(yīng)用水平，對生產(chǎn)設(shè)備和工藝的清潔化、低碳化和循環(huán)化更新，筑牢綠色產(chǎn)業(yè)基礎(chǔ)支撐；二是培育增量，積極培育生態(tài)產(chǎn)業(yè)，尤其是堅持“市場主導(dǎo)+政府引導(dǎo)”模式，推動多元化戰(zhàn)略新興產(chǎn)業(yè)發(fā)展壯大；三是流量更新，積極探索全產(chǎn)業(yè)鏈金融模式，總結(jié)推廣“農(nóng)銀企”產(chǎn)業(yè)共同體創(chuàng)新試點經(jīng)驗，扎實推進(jìn)綠色金融投融資改革和創(chuàng)新，開發(fā)推廣覆蓋特色優(yōu)勢綠色產(chǎn)品的險種，推動商業(yè)銀行將生態(tài)環(huán)保指標(biāo)分析納入貸款流程；四是變量創(chuàng)新，大力推進(jìn)企業(yè)綠色技術(shù)創(chuàng)新，培育科技型綠色產(chǎn)品加工龍頭企業(yè)，進(jìn)一步激活企業(yè)主體的內(nèi)生動力，創(chuàng)造新的經(jīng)濟(jì)增長點。 （2）立足現(xiàn)有資源，推動差異化、特色化及協(xié)同發(fā)展。 黨的十八大以來，重慶市在探索綠色發(fā)展方面取得了重大成就，但市內(nèi)各區(qū)縣由于要素稟賦差異巨大，仍有很長的路要走。 主城都市區(qū)是國家層面的重點開發(fā)區(qū)域，應(yīng)當(dāng)以產(chǎn)業(yè)數(shù)字化和數(shù)字產(chǎn)業(yè)化為轉(zhuǎn)型方向，推動數(shù)字技術(shù)在傳統(tǒng)產(chǎn)業(yè)中的深度應(yīng)用，全面提升制造業(yè)、建筑業(yè)、農(nóng)業(yè)和服務(wù)業(yè)數(shù)字化水平，適當(dāng)控制城鎮(zhèn)化發(fā)展的速度，建設(shè)緊湊低耗高效彈性的都市區(qū)［32］。 渝東南武陵山區(qū)城鎮(zhèn)群屬于國家重點生態(tài)功能區(qū)，應(yīng)當(dāng)按照“面上保護(hù)、點上開花”的原則，立足武陵山區(qū)生態(tài)資源和民族特色，推進(jìn)農(nóng)商文旅體融合發(fā)展，建立健全生態(tài)產(chǎn)品價值實現(xiàn)機(jī)制，帶動生態(tài)保護(hù)地區(qū)和鄉(xiāng)村地區(qū)走出一條生態(tài)優(yōu)先綠色發(fā)展的新路子［33］。 渝東北三峽庫區(qū)城鎮(zhèn)群是重慶市重要的老工業(yè)基地，傳統(tǒng)制造業(yè)轉(zhuǎn)型是重點，應(yīng)當(dāng)以生態(tài)產(chǎn)業(yè)化和產(chǎn)業(yè)生態(tài)化為發(fā)展方向，聚焦綠色食品、綠色建材、特色輕工等細(xì)分領(lǐng)域，高起點承接產(chǎn)業(yè)轉(zhuǎn)移，發(fā)揮數(shù)字金融的支撐作用，統(tǒng)籌構(gòu)建高效、清潔、低碳、循環(huán)的“三峽制造”綠色制造體系［34］。