產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率的空間效應研究
——來自京津冀城市群的經(jīng)驗證據(jù)

高素英，王 迪，馬曉輝

（河北工業(yè)大學 經(jīng)濟管理學院，天津 300401）

一、引 言

“十四五”規(guī)劃綱要關于加快京津冀協(xié)同發(fā)展的新思路中再次提及綠色發(fā)展問題，指出生態(tài)環(huán)境保護應立體治理、山海統(tǒng)籌。京津冀城市群作為國家戰(zhàn)略的重點推進區(qū)域，憑借其獨特的要素稟賦及資源集聚等優(yōu)勢成為引領我國經(jīng)濟增長的重要引擎。然而，在過去我國粗放型經(jīng)濟增長方式的大背景下，其經(jīng)濟績效的取得很大程度上以犧牲環(huán)境績效為代價，有悖于“綠水青山就是金山銀山”的發(fā)展理念。根據(jù)2020年中國生態(tài)環(huán)境部公布的168個重點城市空氣質(zhì)量狀況，在排名后10位的城市中，京津冀地區(qū)占據(jù)4席。生態(tài)環(huán)境問題制約著經(jīng)濟高質(zhì)量發(fā)展，因此，全面構筑京津冀城市群綠色發(fā)展之路至關重要。2021年5月19日，京津冀數(shù)字經(jīng)濟聯(lián)盟成立，旨在借助數(shù)字經(jīng)濟推動產(chǎn)業(yè)協(xié)同互補、環(huán)保協(xié)同共治，該戰(zhàn)略將為京津冀城市群綠色發(fā)展提供重要的政策支持。

綠色全要素生產(chǎn)率是將資源利用和環(huán)境代價納入傳統(tǒng)全要素生產(chǎn)率的核算框架，準確衡量實際投入產(chǎn)出效率的同時又能體現(xiàn)綠色發(fā)展理念，現(xiàn)已成為評價經(jīng)濟增長質(zhì)量的重要指標［1-2］?，F(xiàn)有研究表明，產(chǎn)業(yè)集聚已成為現(xiàn)代城市經(jīng)濟中的重要現(xiàn)象，對于區(qū)域經(jīng)濟發(fā)展及效率提升具有顯著影響［3］。如陳陽和唐曉華（2018）［4］選取城市等級視角驗證了制造業(yè)集聚對城市綠色全要素生產(chǎn)率的影響效應；余奕杉等（2021）［5］研究表明，生產(chǎn)性服務業(yè)集聚對城市綠色全要素生產(chǎn)率亦存在顯著影響。隨著產(chǎn)業(yè)結構調(diào)整和協(xié)同發(fā)展戰(zhàn)略的不斷推進，“單輪驅(qū)動”已不再符合經(jīng)濟高質(zhì)量發(fā)展要求，區(qū)域產(chǎn)業(yè)布局逐漸由傳統(tǒng)的制造業(yè)集聚向制造業(yè)與生產(chǎn)性服務業(yè)協(xié)同集聚模式轉變。集聚的生產(chǎn)率效應研究也逐步轉向協(xié)同集聚視角，且研究發(fā)現(xiàn)協(xié)同集聚對經(jīng)濟效率的影響作用更強［6］。因此，在新經(jīng)濟時代我國加快推動經(jīng)濟綠色轉型以及“十四五”時期加快推動京津冀協(xié)同發(fā)展戰(zhàn)略的雙重背景下，通過有效引導京津冀協(xié)同發(fā)展及促進區(qū)域內(nèi)產(chǎn)業(yè)間滲透融合的方式，促使資源得以合理集約利用，減少環(huán)境污染物的排放，提高綠色全要素生產(chǎn)率水平，促進各城市經(jīng)濟綠色協(xié)調(diào)發(fā)展，具有重要的現(xiàn)實意義。

二、文獻綜述

Ellison和Glaeser（1997）［7］于1997年首次提出產(chǎn)業(yè)協(xié)同集聚的概念，將其界定為異質(zhì)性關聯(lián)產(chǎn)業(yè)間的空間集聚現(xiàn)象?；仡櫸覈a(chǎn)業(yè)發(fā)展進程，制造業(yè)為提高效率，“服務外包”行為大量涌現(xiàn)，逐漸由“以生產(chǎn)為中心”向“以服務為中心”轉型。生產(chǎn)性服務業(yè)作為從制造業(yè)中分離出來的產(chǎn)業(yè)［8］，以“中間投入”屬性與制造業(yè)協(xié)同，實現(xiàn)了產(chǎn)業(yè)垂直一體化向垂直專業(yè)化轉型。Kelle（2013）［9］、唐曉華等（2018）［10］認為，制造業(yè)與生產(chǎn)性服務業(yè)之間存在著相輔相成的動態(tài)協(xié)同發(fā)展關系，這一互動關聯(lián)則為產(chǎn)業(yè)協(xié)同集聚奠定了基礎。兩業(yè)協(xié)同集聚有助于優(yōu)化經(jīng)濟結構、盤活城市經(jīng)濟，對促進經(jīng)濟高質(zhì)量發(fā)展具有重要的驅(qū)動作用［11-13］。

梳理文獻發(fā)現(xiàn)，選取制造業(yè)及生產(chǎn)性服務業(yè)視角分析產(chǎn)業(yè)集聚對全要素生產(chǎn)率影響的研究較為成熟，可概括為兩個層面。第一，單一產(chǎn)業(yè)集聚對全要素生產(chǎn)率的影響。以Marshall（1890）［14］和Jacobs（1969）［15］為代表的經(jīng)典外部性理論認為，產(chǎn)業(yè)集聚外部性能夠提高地區(qū)競爭力，促進經(jīng)濟增長；程中華和張立柱（2015）［16］通過空間計量研究發(fā)現(xiàn)，制造業(yè)集聚對全要素生產(chǎn)率具有負向作用，且這一影響主要是由Marshall和Porter倡導的外部性主導；而楊浩昌等（2018）則認為，制造業(yè)集聚和生產(chǎn)性服務業(yè)集聚均有助于促進城市全要素生產(chǎn)率提升［17］。第二，兩業(yè)協(xié)同集聚對全要素生產(chǎn)率的影響。相較于單一產(chǎn)業(yè)集聚，產(chǎn)業(yè)協(xié)同集聚作為經(jīng)濟結構轉型發(fā)展的主要方式，其對經(jīng)濟增長效率能夠產(chǎn)生重要影響。國內(nèi)研究中，伍先福（2018）［18］關于制造業(yè)與生產(chǎn)性服務業(yè)協(xié)同集聚對全要素生產(chǎn)率的影響效應研究起步較早，并指出兩者間隨產(chǎn)業(yè)專業(yè)化集聚度變化呈非線性關系；王靜田等（2021）［19］認為，大量生產(chǎn)要素聚集容易導致擠占其他地區(qū)資源，相關企業(yè)會因沉沒成本難以自由進出市場，協(xié)同集聚所帶來的負鎖定效應可能會降低城市全要素生產(chǎn)率。隨著經(jīng)濟高質(zhì)量發(fā)展，對經(jīng)濟增長提出了更高的要求，學術界關于產(chǎn)業(yè)集聚對綠色全要素生產(chǎn)率的研究愈發(fā)豐富。早期研究中，一些學者認為集聚的擁擠效應會對集聚區(qū)域內(nèi)的環(huán)境承載力發(fā)起挑戰(zhàn)，不利于經(jīng)濟可持續(xù)發(fā)展［20-21］；與之相反，任陽軍等（2019）［22］通過構建GMM模型驗證了生產(chǎn)性服務業(yè)及制造業(yè)集聚有利于提高綠色全要素生產(chǎn)率；朱風慧和劉立峰（2021）［23］證實了制造業(yè)集聚對綠色全要素生產(chǎn)率的影響并非呈簡單的線性關系，兩者間關系隨對外開放程度的增加由負轉正。

綜上，首先，因研究方法或研究范圍的不同，學術界關于集聚的生產(chǎn)率效應尚未得到統(tǒng)一定論，還需進一步研究；其次，目前聚焦于制造業(yè)與生產(chǎn)性服務業(yè)協(xié)同集聚能否有效提高綠色全要素生產(chǎn)率的研究尚少，少量相關研究也是以環(huán)境污染或創(chuàng)新效率為落腳點［24-25］?！笆奈濉睍r期是京津冀協(xié)同發(fā)展戰(zhàn)略承上啟下的關鍵時期，各城市間經(jīng)濟綠色協(xié)調(diào)可持續(xù)發(fā)展問題更加值得關注。鑒于此，本文運用基于非徑向、非角度包含非期望產(chǎn)出的DEA-SBM模型同時結合GML生產(chǎn)率指數(shù)，對京津冀13個城市的綠色全要素生產(chǎn)率進行測算，并通過相對差異指數(shù)法測算制造業(yè)與生產(chǎn)性服務業(yè)協(xié)同集聚指數(shù)，就測算結果實證分析京津冀城市群產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率的影響效應。

三、理論分析與研究假設

制造業(yè)與生產(chǎn)性服務業(yè)協(xié)同集聚并不是簡單地將單一產(chǎn)業(yè)集聚效應雙倍擴大，兩個異質(zhì)性產(chǎn)業(yè)間能夠通過互動關聯(lián)形成“1+1＞2”的效應作用于京津冀城市群綠色全要素生產(chǎn)率。

一方面，產(chǎn)業(yè)協(xié)同集聚可通過正外部效應促進城市群綠色全要素生產(chǎn)率提高。一是規(guī)模經(jīng)濟效應。根據(jù)新經(jīng)濟地理理論，制造業(yè)與生產(chǎn)性服務業(yè)在城市內(nèi)集聚有助于加強異質(zhì)性產(chǎn)業(yè)間的聯(lián)系，通過降低運輸成本及共享集聚區(qū)域公共資源等使得交易成本降低［26］，提高資源利用率的同時降低了環(huán)境污染。此外，制造業(yè)發(fā)展在很大程度上依賴于能源消耗，生產(chǎn)性服務業(yè)自帶“清潔”及“高效”的屬性，有助于促進制造業(yè)向低污染、高創(chuàng)新的方向發(fā)展。為疏解非首都功能，北京市制造業(yè)向津冀轉移，規(guī)模經(jīng)濟效應有助于緩解因京津冀城市群自身“核心-外圍”特征所引致的環(huán)境問題。二是知識溢出效應。制造業(yè)與生產(chǎn)性服務業(yè)協(xié)同集聚有助于豐富集聚區(qū)域內(nèi)公共知識池，促使產(chǎn)業(yè)間難以傳遞的隱形知識顯性化［25］，便于京津冀在環(huán)境污染與經(jīng)濟協(xié)調(diào)發(fā)展問題上達成一致。三是技術進步效應。兩大異質(zhì)產(chǎn)業(yè)間既能夠相互借鑒學習同時也存在著競爭關系，有利于淘汰高污染產(chǎn)業(yè)，促使經(jīng)濟結構清潔化。集聚區(qū)內(nèi)企業(yè)通過與高校、科研院所和政府等組織合作還有利于推進產(chǎn)學研一體化，進一步提高京津冀城市群特別是河北省11個城市的綠色技術發(fā)展水平。

另一方面，產(chǎn)業(yè)協(xié)同集聚也會通過負外部效應制約城市群綠色全要素生產(chǎn)率提高。一是環(huán)境擁擠效應。根據(jù)Hoover（1948）［27］提出的產(chǎn)業(yè)集聚最佳規(guī)模論，產(chǎn)業(yè)集聚效應并非始終有效，而是存在一個最佳規(guī)模，當集聚超過最優(yōu)臨界值后，規(guī)模效應將不再顯著。要素的持續(xù)過度投入會導致規(guī)模報酬遞減，帶來要素擁擠問題，最終對集聚區(qū)域的資源環(huán)境承載能力提出挑戰(zhàn)。某城市的環(huán)境承載能力超過上限，則不利于整個城市群綠色協(xié)同發(fā)展。二是路徑鎖定效應。當集聚發(fā)展到一定階段時，集聚區(qū)域內(nèi)便容易形成路徑鎖定，僅限于接收本區(qū)域內(nèi)部的信息而“拒收”區(qū)域外的新信息和新技術，封閉式的學習模式會阻礙本城市群減排新技術的變革。此外，京津冀各城市間的產(chǎn)業(yè)發(fā)展現(xiàn)狀存在著嚴重的梯度化特征，還會容易導致“搭便車”現(xiàn)象發(fā)生。三是能源回彈效應。根據(jù)杰文斯悖論［28］，協(xié)同集聚所帶來的技術進步效應有助于促進經(jīng)濟快速增長，由此會對能源產(chǎn)生新的需求，最終導致因技術進步效應節(jié)約的能源甚至可能被經(jīng)濟快速增長所帶來的新耗能源完全抵消，如此一來，則會導致技術水平相對落后的城市無法有效依托技術進步效應提高經(jīng)濟發(fā)展質(zhì)量。

綜上所述，當產(chǎn)業(yè)協(xié)同集聚水平處于適度規(guī)模時，其可通過正外部效應促進綠色全要素生產(chǎn)率提高；反之，處于過度協(xié)同集聚狀態(tài)時則會通過負外部效應發(fā)揮抑制作用。基于此，繪制作用機制圖，如圖1所示，并提出本文假設1。

圖1 產(chǎn)業(yè)協(xié)同集聚對京津冀城市群綠色全要素生產(chǎn)率的影響機制

H1：協(xié)同集聚的生產(chǎn)率效應存在最佳規(guī)模，產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率存在“倒U”型影響。

此外，京津冀城市群內(nèi)存在的要素流動、資源共享等通常使得各城市間存在空間交互作用，現(xiàn)實中相鄰城市之間也會相互模仿學習［29］。因此，某一城市的產(chǎn)業(yè)協(xié)同集聚水平不僅僅影響本城市內(nèi)的綠色全要素生產(chǎn)率，往往也會對相鄰城市的綠色全要素生產(chǎn)率水平帶來空間溢出效應影響。據(jù)此，本文提出假設2。

H2：產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率的影響存在空間溢出效應。

四、實證設計

（一）空間計量模型構建

根據(jù)前文分析，為檢驗產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率是否存在非線性影響，本文在模型中引入產(chǎn)業(yè)協(xié)同集聚指數(shù)的二次項。此外，本文同時將空間因素納入模型之中，構建空間計量模型探究產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率的空間效應影響?？臻g杜賓模型作為空間誤差模型及空間滯后模型的一般形式，具有能夠同時捕捉地區(qū)空間異質(zhì)性和競爭特點的優(yōu)勢。因此，本文采用空間杜賓模型，設定如下：

其中：GTFPit為被解釋變量，即綠色全要素生產(chǎn)率；COit為核心解釋變量，即產(chǎn)業(yè)協(xié)同集聚指數(shù)；C表示核心解釋變量產(chǎn)業(yè)協(xié)同集聚的二次項；ISit、FDIit、PGDPit、GOINit均為控制變量，分別代表產(chǎn)業(yè)結構、外資依賴程度、經(jīng)濟發(fā)展水平和政府科技投入；i和t分別表示城市和年份；ρ為空間滯后系數(shù)，表示相鄰城市綠色全要素生產(chǎn)率對本城市綠色全要素生產(chǎn)率觀測值的影響方向和程度；Wij為經(jīng)濟距離權重矩陣；β0為常數(shù)項；β1-β6為各解釋變量的回歸系數(shù)；θ1-θ6為各解釋變量相應的空間回歸系數(shù)；μi為地區(qū)固定效應；νt為時間固定效應；εit為隨機擾動項。

此外，空間權重矩陣是反映各單元之間空間關系的重要工具，在設定空間計量模型時還需構建合適的空間權重矩陣?？紤]京津冀城市間存在著經(jīng)濟發(fā)展不平衡的問題，故本文選取經(jīng)濟距離空間權重矩陣進行空間計量實證分析。經(jīng)濟距離空間權重矩陣構建如下：

其中，Y i、Y j分別代表京津冀地區(qū)中i市和j市的人均GDP。

（二）變量選取及測度

1.被解釋變量

本文被解釋變量為綠色全要素生產(chǎn)率（GTFP）。參考Fukuyama和Weber（2009）［30］、Oh（2010）［31］的研究，采用DEA-SBM模型并構造GML指數(shù)對2006—2019年京津冀城市群的綠色全要素生產(chǎn)率進行測度，所需指標包括投入指標、期望產(chǎn)出以及非期望產(chǎn)出，具體見表1所列。

首先，依據(jù)Fukuyama和Weber的研究，構建同時考慮能源投入以及環(huán)境污染產(chǎn)出的SBM模型如式（3）所示。其中：( )xti,y ti,d ti為t時期i城市的投入產(chǎn)出向量表示方向向量表示投入產(chǎn)出的松弛變量。

在此基礎上，本文依據(jù)Oh的研究，采用SBMGML指數(shù)法進一步測算京津冀城市群綠色全要素生產(chǎn)率的動態(tài)變化，測算公式如式（4）所示。由于GML指數(shù)測算的是相鄰年份間的環(huán)比增長，且其具有累積性，因此本文以2005年為基期，取其為1，將測算出的GML指數(shù)通過累積轉換為最終的綠色全要素生產(chǎn)率，以此進行實證分析。

2.核心解釋變量

本文核心解釋變量為產(chǎn)業(yè)協(xié)同集聚（CO）。借鑒陳建軍等（2016）［32］的研究，對2006—2019年京津冀城市群制造業(yè)與生產(chǎn)性服務業(yè)協(xié)同集聚水平加以測算。

首先，利用區(qū)位熵分別測算制造業(yè)或生產(chǎn)性服務業(yè)單一產(chǎn)業(yè)集聚程度，區(qū)位熵公式為：

其中：LQiu為i城市u產(chǎn)業(yè)的集聚指數(shù)；q iu為i城市u產(chǎn)業(yè)的就業(yè)人數(shù)；qi為i城市所有產(chǎn)業(yè)的就業(yè)人數(shù)；qu為全國u產(chǎn)業(yè)的就業(yè)人數(shù)；q為全國總就業(yè)人數(shù)。

其次，采用產(chǎn)業(yè)協(xié)同集聚相對差異指數(shù)法測算產(chǎn)業(yè)協(xié)同集聚指數(shù)，測算公式如下：

其中，LQim、LQip分別代表i城市制造業(yè)及生產(chǎn)性服務業(yè)的區(qū)位熵指數(shù)。測算出的產(chǎn)業(yè)協(xié)同集聚指數(shù)越大，代表制造業(yè)與生產(chǎn)性服務業(yè)協(xié)同集聚水平越高。

3.控制變量

通過對影響綠色全要素生產(chǎn)率的相關文獻加以梳理可知，除本文研究的核心解釋變量產(chǎn)業(yè)協(xié)同集聚外，還有一些其他因素會影響綠色全要素生產(chǎn)率變化。本文結合已有文獻，選取如下控制變量：①產(chǎn)業(yè)結構（IS）。為促進產(chǎn)業(yè)轉型升級，經(jīng)濟資源通常會從低生產(chǎn)率產(chǎn)業(yè)流向高生產(chǎn)率產(chǎn)業(yè)，有利于提高綠色全要素生產(chǎn)率。②外資依賴程度（FDI）。FDI對綠色全要素生產(chǎn)率存在不確定性影響，外資增加一方面有利于引進國外先進技術促進綠色全要素生產(chǎn)率的提高，同時國外污染密集型產(chǎn)業(yè)因規(guī)避本國環(huán)境規(guī)制政策進行產(chǎn)業(yè)轉移也會制約外資流入國的綠色經(jīng)濟發(fā)展。③經(jīng)濟發(fā)展水平（PGDP）。一般而言，經(jīng)濟發(fā)展水平的提高有利于環(huán)境治理，但過快的經(jīng)濟發(fā)展往往也會以犧牲環(huán)境為代價。④政府科技投入（GOIN）?？萍纪度胗欣诖龠M經(jīng)濟增長，也有利于提高相應的污染治理技術、增強環(huán)境污染治理能力。

各變量的指標選取及計算方法見表1所列。

表1 變量選取及計算方法

（三）數(shù)據(jù)來源與描述性統(tǒng)計

以上指標均以京津冀地區(qū)13個城市為研究對象，研究區(qū)間為2006—2019年，涉及的數(shù)據(jù)主要來源于《中國城市統(tǒng)計年鑒》《中國統(tǒng)計年鑒》《北京統(tǒng)計年鑒》《天津統(tǒng)計年鑒》《河北經(jīng)濟年鑒》，少量缺失數(shù)據(jù)通過查閱各市對應年份的國民經(jīng)濟和社會發(fā)展統(tǒng)計公報獲得。此外，為保證數(shù)據(jù)的平穩(wěn)性，對各變量進行取對數(shù)處理。各變量原始數(shù)據(jù)的描述性統(tǒng)計結果見表2所列。

表2 變量的描述性統(tǒng)計結果

五、實證結果分析

（一）空間相關性檢驗

1.全局空間自相關檢驗

根據(jù)表3的莫蘭檢驗指數(shù)結果可以看出，各年份均能通過顯著性檢驗，并且表現(xiàn)為正相關，由此說明京津冀地區(qū)各城市間的綠色全要素生產(chǎn)率存在著較為顯著的空間相關性。因此，在研究產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率的影響時，有必要將空間因素考慮在內(nèi)。

表3 京津冀城市群綠色全要素生產(chǎn)率的全局莫蘭指數(shù)

2.局部空間自相關檢驗

根據(jù)全局莫蘭指數(shù)檢驗結果，京津冀城市間綠色全要素生產(chǎn)率整體呈現(xiàn)正向聚集特征，但是全局Moran′sI無法具體反映各個城市間的空間分布及特點，進而選取2006年、2011年、2015年和2019年的綠色全要素生產(chǎn)率進行局部空間自相關檢驗。從圖2可以看出，2006年、2011年、2015年和2019年京津冀城市群綠色全要素生產(chǎn)率均呈現(xiàn)正向空間自相關，即大部分城市分布在第一象限和第三象限，以“H-H”和“L-L”類型集聚為主。

圖2 綠色全要素生產(chǎn)率局部莫蘭指數(shù)散點圖

（二）空間計量模型檢驗

在進行空間計量分析之前，還需進行相關檢驗以確定合適的空間計量模型，結果見表4所列。首先，參照Anselin（1988）［33］的研究，利用LM方法對SEM模型以及SAR模型進行有效性檢驗。由表4可以看出，兩個模型的LM值均在1%的顯著性水平上顯著，其Robust LM值也均在1%的顯著性水平上顯著，拒絕無空間誤差和空間滯后效應的原假設，即本文選擇SEM模型和SAR模型均較為合適。其次，Hausman檢驗結果表明，固定效應模型更優(yōu)。進一步地，本文需要借助LR及Wald統(tǒng)計量對SDM模型的適用性加以檢驗。結果顯示，空間杜賓模型的LR及Wald聯(lián)合檢驗均能夠通過1%的顯著性檢驗，說明本文建立空間杜賓模型是有效的。

表4 空間計量模型檢驗結果

（三）基本回歸結果分析

1.空間計量結果分析

為提高模型估計結果的有效性和穩(wěn)健性，本文基于經(jīng)濟距離權重矩陣對SEM、SAR和SDM模型的回歸結果進行對比分析，空間計量結果見表5所列。

從表5可以看出，核心解釋變量產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率影響的一次項系數(shù)在三個模型中始終為正，二次項系數(shù)在三個模型中始終為負，并均能夠通過顯著性檢驗，說明產(chǎn)業(yè)協(xié)同集聚對京津冀城市群綠色全要素生產(chǎn)率存在“倒U”型影響，H1成立。這一結果也與王燕和孫超（2020）［34］等學者所得出的研究結果一致，即產(chǎn)業(yè)協(xié)同集聚存在適度規(guī)模，符合Hoover（1948）［27］的產(chǎn)業(yè)集聚最佳規(guī)模論。具體來看，基于SDM的實證結果，求得產(chǎn)業(yè)協(xié)同集聚指數(shù)的拐點值為1.084，將其換算為原始數(shù)據(jù)后對應的結果為2.956。說明當京津冀城市群產(chǎn)業(yè)協(xié)同集聚指數(shù)小于2.956時，正外部效應發(fā)揮作用，其對綠色全要素生產(chǎn)率的影響處于“倒U”型曲線的上升段；而當協(xié)同集聚指數(shù)大于2.956時，產(chǎn)業(yè)協(xié)同集聚則會對綠色全要素生產(chǎn)率產(chǎn)生負外部效應，表現(xiàn)為抑制作用。

表5 空間計量結果

根據(jù)本文測算結果，2019年京津冀城市群中秦皇島市、唐山市的產(chǎn)業(yè)協(xié)同集聚水平已經(jīng)超過了拐點值，處于過度集聚狀態(tài)，其余11個城市的產(chǎn)業(yè)協(xié)同集聚水平均在拐點值以下。秦皇島市臨近京津，地處環(huán)渤海經(jīng)濟圈中心地帶，具有獨特的自然稟賦，近年來，秦皇島市不斷推動產(chǎn)業(yè)轉型升級，其產(chǎn)業(yè)協(xié)同集聚水平較高。唐山市地處渤海灣中心地帶，工業(yè)基礎雄厚，唐山灣“四點一帶”區(qū)域發(fā)展規(guī)劃成效顯著，促進了唐山市經(jīng)濟結構優(yōu)化升級。因此，秦皇島市和唐山市應適當降低產(chǎn)業(yè)協(xié)同集聚程度，使之回歸于拐點值水平；對于其余11個城市而言，則應進一步加大兩業(yè)協(xié)同集聚程度，特別是對于滄州市、廊坊市和保定市，這些城市的綠色全要素生產(chǎn)率尚處于較低水平，需以產(chǎn)業(yè)協(xié)同集聚為抓手趕追其他城市，以期提高整個城市群綠色經(jīng)濟效率。

此外，從控制變量來看，產(chǎn)業(yè)結構具有顯著的正向影響，產(chǎn)業(yè)結構升級有利于提高綠色全要素生產(chǎn)率；外資依賴程度具有正向影響，但不顯著；經(jīng)濟發(fā)展水平具有一定的負向影響，說明城市經(jīng)濟發(fā)展過快時會犧牲環(huán)境績效；政府科技投入對綠色全要素生產(chǎn)率未表現(xiàn)出正向作用。綠色全要素生產(chǎn)率在三個模型中的空間自回歸系數(shù)分別為0.667、0.582、0.382，且均能夠通過顯著性檢驗，由此表明空間因素對綠色全要素生產(chǎn)率具有顯著的正向影響。

2.空間溢出效應分解

Lesage和Pace（2009）［35］認為，當空間計量模型含有內(nèi)生空間交互效應時，無法直接利用點估計結果來判定空間溢出效應是否存在。因此，為進一步驗證本文假設，參考Lesage和Pace的研究，運用偏微分方法深入分解影響效應，空間杜賓模型溢出效應分解結果見表6所列。

表6 空間杜賓模型溢出效應分解結果

根據(jù)表6，從核心解釋變量的直接效應來看，產(chǎn)業(yè)協(xié)同集聚的一次項系數(shù)顯著為正，二次項系數(shù)顯著為負，同樣說明適度規(guī)模的產(chǎn)業(yè)協(xié)同集聚能夠促進京津冀地區(qū)本城市綠色全要素生產(chǎn)率提高，過度集聚則會抑制本城市綠色全要素生產(chǎn)率增加，H1再次得以驗證。從間接效應來看，產(chǎn)業(yè)協(xié)同集聚的一次項系數(shù)為正，二次項系數(shù)為負，且均通過1%的顯著性檢驗，說明適度規(guī)模的產(chǎn)業(yè)協(xié)同集聚同樣能夠帶動相鄰城市綠色全要素生產(chǎn)率提高，驗證了京津冀城市群產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率具有顯著的空間溢出效應影響，H2成立。在京津冀協(xié)同發(fā)展背景下，三地在交通方面共同謀劃多條高速鐵路、公路及城際列車；在港口方面深化改革通關一體化；在生態(tài)方面實施大氣污染聯(lián)防聯(lián)控；在產(chǎn)業(yè)對接方面加大三地間資金流動。京津冀協(xié)同發(fā)展戰(zhàn)略增強了城市間的要素流動，加強了各城市生產(chǎn)單位之間的聯(lián)系，進而表現(xiàn)出相鄰城市間綠色全要素生產(chǎn)率的正向溢出效應。

控制變量方面，產(chǎn)業(yè)結構對綠色全要素生產(chǎn)率的直接效應系數(shù)為0.358，在1%的顯著性水平上顯著，間接效應系數(shù)不顯著，說明產(chǎn)業(yè)結構僅能夠促進京津冀地區(qū)本城市綠色全要素生產(chǎn)率上升；經(jīng)濟發(fā)展水平對綠色全要素生產(chǎn)率的直接效應顯著為負，間接效應顯著為正，這表明當經(jīng)濟發(fā)展水平過高時，會導致環(huán)境洛倫茲曲線處于“倒U”型的下降段，城市污染問題會由于過度追求經(jīng)濟增長速度而加重，但空間外溢效應使其能夠促進相鄰城市的綠色全要素生產(chǎn)率上升；外資依賴程度和政府科技投入對綠色全要素生產(chǎn)率的分解結果不顯著。

（四）進一步分析

1.產(chǎn)業(yè)異質(zhì)性分析

本文參照楊桐彬等（2020）［36］對生產(chǎn)性服務業(yè)的分類方法，將“交通運輸、倉儲和郵政業(yè)”“租賃和商務服務業(yè)”這兩類低附加值產(chǎn)業(yè)歸為低端生產(chǎn)性服務業(yè)，將“信息傳輸、計算機服務和軟件業(yè)”“金融業(yè)”“科學研究和技術服務業(yè)”歸為高端生產(chǎn)性服務業(yè)。制造業(yè)與高端生產(chǎn)性服務業(yè)或低端生產(chǎn)性服務業(yè)協(xié)同集聚的測算方法同前。

基于此，進一步運用SDM模型考察產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率的影響是否存在產(chǎn)業(yè)異質(zhì)性，結果見表7所列。

根據(jù)表7，首先，就直接效應而言，制造業(yè)與高端生產(chǎn)性服務業(yè)協(xié)同集聚的一次項系數(shù)為8.706，二次項系數(shù)為-3.709，制造業(yè)與低端生產(chǎn)性服務業(yè)協(xié)同集聚的一次項系數(shù)為4.655，二次項系數(shù)為-2.279，且均通過了1%的顯著性檢驗。由此說明，無論是制造業(yè)與生產(chǎn)性服務業(yè)協(xié)同集聚，抑或是制造業(yè)與其細分行業(yè)協(xié)同集聚，均會對京津冀城市群綠色全要素生產(chǎn)率產(chǎn)生“倒U”型影響，再次驗證了產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率的影響作用存在最佳規(guī)模。對比分析，制造業(yè)與高端生產(chǎn)性服務業(yè)協(xié)同集聚相比于低端生產(chǎn)性服務業(yè)對京津冀城市群綠色全要素生產(chǎn)率的直接影響效應更大，且具有更高的拐點值，說明制造業(yè)與高端生產(chǎn)性服務業(yè)協(xié)同集聚對京津冀城市群綠色全要素生產(chǎn)率的影響會更晚進入拐點，能夠延遲負外部效應影響。其次，就間接效應來看，無論是制造業(yè)與高端生產(chǎn)性服務業(yè)協(xié)同集聚，抑或是制造業(yè)與低端生產(chǎn)性服務業(yè)協(xié)同集聚，其對綠色全要素生產(chǎn)率均存在空間溢出效應影響，但制造業(yè)與高端生產(chǎn)性服務業(yè)對相鄰城市綠色全要素生產(chǎn)率的“帶動”作用更強。

表7 產(chǎn)業(yè)異質(zhì)性分析回歸結果

2.對綠色全要素生產(chǎn)率分解項的影響分析

上文實證檢驗了制造業(yè)與生產(chǎn)性服務業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率的空間效應影響，包括總效應、直接效應和間接效應。本部分則進一步探究京津冀城市群產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率分解項的空間效應影響，即產(chǎn)業(yè)協(xié)同集聚對綠色技術效率（GEC）和綠色技術進步（GTC）的影響，做進一步回歸，結果見表8所列。

表8 綠色全要素生產(chǎn)率分解項的基本回歸結果

從表8綠色全要素生產(chǎn)率分解項的基本回歸結果可以看出，在SDM模型下，產(chǎn)業(yè)協(xié)同集聚對綠色技術效率和綠色技術進步影響的一次項系數(shù)均顯著為正，二次項系數(shù)均顯著為負，由此表明，在空間計量模型中，京津冀城市群產(chǎn)業(yè)協(xié)同集聚能夠通過綠色技術效率和綠色技術進步“雙引擎”影響綠色全要素生產(chǎn)率。從控制變量來看，產(chǎn)業(yè)結構對綠色技術效率和綠色技術進步具有一定的正向影響；外資依賴程度對綠色技術效率和綠色技術進步未表現(xiàn)出顯著影響；經(jīng)濟發(fā)展水平對綠色技術效率存在顯著的負向影響，而在SAR模型下對綠色技術進步表現(xiàn)出正向影響；政府科技投入對綠色技術效率和綠色技術進步均未表現(xiàn)出正向影響。此外，從空間自回歸系數(shù)ρ/λ來看，綠色全要素生產(chǎn)率分解項的空間自回歸系數(shù)均在1%的顯著性水平上顯著，由此表明京津冀城市間的綠色技術效率和綠色技術進步均具有空間溢出效應，相鄰城市的綠色技術效率和綠色技術進步提高有助于帶動本城市綠色技術效率和綠色技術進步水平的提升。

（五）穩(wěn)健性檢驗

空間權重矩陣的選取對于空間計量模型的結果具有很大影響，因此本文參照現(xiàn)有研究中的主流做法，通過更換空間權重矩陣對實證結果進行穩(wěn)健性檢驗?；凇?-1”鄰接權重矩陣的穩(wěn)健性檢驗結果見表9所列。

表9 穩(wěn)健性檢驗結果

從表9的穩(wěn)健性檢驗結果可以看出，基于“0-1”鄰接權重矩陣與前述經(jīng)濟距離權重矩陣的估計結果基本一致，由此說明本文的實證結果具有穩(wěn)健性。

六、結論與啟示

本文以京津冀城市群為例，構建空間杜賓模型實證檢驗制造業(yè)與生產(chǎn)性服務業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率的直接影響效應以及空間溢出效應。研究結論如下：

第一，產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率存在直接影響。在正、負外部效應的雙重作用下，產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率具有非線性影響，表現(xiàn)為“倒U”型。即只有協(xié)同集聚水平處于適度規(guī)模時，其對城市綠色全要素生產(chǎn)率才能有效發(fā)揮促進作用。

第二，空間溢出效應存在。一方面，在京津冀協(xié)同發(fā)展戰(zhàn)略背景下，各城市間加強了要素流動及聯(lián)系，因此本城市綠色經(jīng)濟效率受相鄰城市綠色全要素生產(chǎn)率的影響；另一方面，產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率存在空間溢出效應，城市間能夠通過發(fā)揮“示范”作用影響鄰近城市綠色經(jīng)濟發(fā)展水平。

第三，進一步分析表明存在產(chǎn)業(yè)異質(zhì)性以及“雙引擎”作用路徑。一方面，產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率的影響存在產(chǎn)業(yè)異質(zhì)性。相比而言，高端生產(chǎn)性服務業(yè)集聚了更高端的人力、技術等要素，因此制造業(yè)與高端生產(chǎn)性服務業(yè)協(xié)同集聚的正外部效應更強。另一方面，就產(chǎn)業(yè)協(xié)同集聚對綠色全要素生產(chǎn)率的作用路徑而言，京津冀城市群產(chǎn)業(yè)協(xié)同集聚能夠通過綠色技術效率和綠色技術進步“雙引擎”影響綠色全要素生產(chǎn)率。

結合上述研究結論，得到以下政策啟示：

第一，把控集聚區(qū)域產(chǎn)業(yè)協(xié)同集聚程度。各城市在鼓勵產(chǎn)業(yè)協(xié)同集聚時，不能一味重視“量”的積累，而應注意“度”的把控，因地制宜地結合本區(qū)域產(chǎn)業(yè)協(xié)同集聚現(xiàn)狀適當調(diào)整產(chǎn)業(yè)結構。以本文測算結果為例，2019年京津冀地區(qū)大多城市的產(chǎn)業(yè)協(xié)同集聚程度均在拐點值以內(nèi)，而唐山市和秦皇島市的產(chǎn)業(yè)協(xié)同集聚程度已經(jīng)超過拐點值。因此，其他城市可繼續(xù)適度增加協(xié)同集聚規(guī)模，唐山和秦皇島兩市則應適度降低協(xié)同集聚規(guī)模，使各城市集聚區(qū)域內(nèi)的協(xié)同集聚正外部效應均能最大化，以期促進整個城市群綠色發(fā)展。

第二，以數(shù)字經(jīng)濟加速城市間要素流動與共享。京津冀三地應落實數(shù)字經(jīng)濟聯(lián)盟戰(zhàn)略，以命運共同體視角為產(chǎn)業(yè)協(xié)同集聚空間效應提供條件，減少在綠色經(jīng)濟建設進程中的地方保護主義。一方面促進先進數(shù)字技術在13個城市間的流動與共享，另一方面借助人工智能和大數(shù)據(jù)等技術實現(xiàn)數(shù)據(jù)共享，促進污染聯(lián)防聯(lián)治，以數(shù)字經(jīng)濟賦能京津冀城市群經(jīng)濟綠色協(xié)同可持續(xù)發(fā)展。

第三，著重引導制造業(yè)與高端生產(chǎn)性服務業(yè)深度融合。高端生產(chǎn)性服務業(yè)因“高智、高效、高資”的自身優(yōu)勢，其與制造業(yè)協(xié)同集聚所產(chǎn)生的正外部效應更強。因此，京津冀城市群應根據(jù)優(yōu)勢產(chǎn)業(yè)所帶來的經(jīng)濟效應，著重引導制造業(yè)與信息傳輸、計算機服務和軟件業(yè)以及科學研究和技術服務業(yè)等高端生產(chǎn)性服務業(yè)深度融合，同時也要注重提高低端生產(chǎn)性服務業(yè)的知識技術水平，發(fā)動綠色技術效率和綠色技術進步“雙引擎”作用，助推京津冀城市群經(jīng)濟綠色發(fā)展。