賈 建 琦,趙 林,2*，高 曉 彤，曹 乃 剛

(1.曲阜師范大學地理與旅游學院，山東 日照 276826;2.遼寧師范大學海洋經(jīng)濟與可持續(xù)發(fā)展研究中心，遼寧 大連 116029)

0 引言

全球可持續(xù)發(fā)展背景下，包容性綠色增長(Inclusive Green Growth，IGG)已成為國際社會關注的熱點問題[1]。伴隨著2030年新可持續(xù)發(fā)展目標的公布，在經(jīng)濟增長、社會包容和環(huán)境保護三維目標驅動下，包容性綠色增長被認為是實現(xiàn)可持續(xù)發(fā)展的重要途徑[2]。此外，中共十九屆五中全會對生態(tài)文明建設、綠色發(fā)展以及增進民生福祉提出了新要求， “包容”和“綠色”的發(fā)展模式逐漸成為我國轉變發(fā)展方式、實現(xiàn)高質量發(fā)展的重要途徑。包容性綠色增長效率(Inclusive Green Growth Efficiency，IGGE)是衡量包容性綠色增長水平的重要指標，提升IGGE是實現(xiàn)包容性綠色增長的客觀要求。環(huán)渤海地區(qū)作為我國經(jīng)濟發(fā)展的三大增長極之一，是京津冀協(xié)同發(fā)展、東北振興、黃河流域生態(tài)保護和高質量發(fā)展等國家戰(zhàn)略的覆蓋區(qū)域，但由于工業(yè)化和城市化的迅速推進，資源約束趨緊、生態(tài)環(huán)境破壞、區(qū)域經(jīng)濟發(fā)展失衡等問題凸顯，嚴重制約了環(huán)渤海地區(qū)的可持續(xù)發(fā)展進程。2015年《環(huán)渤海地區(qū)合作發(fā)展綱要》明確提出了生態(tài)保護、開放經(jīng)濟和協(xié)調(diào)發(fā)展等具體要求[3]，區(qū)域傳統(tǒng)粗放的發(fā)展方式亟須轉型，包容性綠色增長成為該地區(qū)實現(xiàn)高質量發(fā)展的重要選擇。與此同時，受區(qū)域協(xié)調(diào)發(fā)展戰(zhàn)略和市場調(diào)控作用的雙重影響，區(qū)域間經(jīng)濟聯(lián)系呈現(xiàn)復雜網(wǎng)絡形態(tài)[4,5]。因此，研究典型區(qū)域包容性綠色增長的空間關聯(lián)網(wǎng)絡結構及其影響機制，對于實現(xiàn)區(qū)域經(jīng)濟高質量發(fā)展具有重要意義。

目前，包容性綠色增長作為一個綜合性概念，現(xiàn)有定義較為宏觀，并未形成明確范圍。國內(nèi)外學者圍繞可持續(xù)發(fā)展的現(xiàn)實要求對包容性綠色增長的概念與目標進行諸多討論，部分國際組織及學者從發(fā)展經(jīng)濟學的角度對其進行了界定，如世界銀行(WB)認為包容性綠色增長是可持續(xù)發(fā)展的經(jīng)濟[6]，國際貨幣基金組織(IMF)則認為包容性綠色增長是在經(jīng)濟增長的同時消除貧困和環(huán)境損害以實現(xiàn)可持續(xù)發(fā)展的模式[7]；部分學者則從福利經(jīng)濟學的角度進行闡述，認為包容性綠色增長要平衡增長、綠色和包容性之間的關系，同時關注當下和未來人們的福利增長[8,9]?？傊菪跃G色增長逐步形成以社會、經(jīng)濟與環(huán)境協(xié)調(diào)發(fā)展為核心，以經(jīng)濟增長、社會公平、福利普惠、資源高效和環(huán)境保護為基本訴求的概念體系[10]。包容性綠色增長的測度方法主要包括主/客觀賦權法和數(shù)據(jù)包絡分析法，已有學者通過構建評價指標體系并賦予指標權重，分別測度了埃塞俄比亞、世界各國、G20國家、中國少數(shù)民族地區(qū)、中國省際和長江經(jīng)濟帶等不同區(qū)域尺度的包容性綠色增長水平[1，11-17]，其中，王宇昕等[15]和吳武林等[18]運用空間探索性分析方法分別證實了長江經(jīng)濟帶和中國包容性綠色增長存在空間正相關。近年來，部分學者從投入產(chǎn)出角度對包容性綠色增長效率進行測度。例如：李政大等通過DEA前沿分解模型刻畫中國包容性綠色增長水平[19]；趙林等運用Super-SBM(Slack-Based Measure)模型測度中國省際包容性綠色增長效率并分析其時空分異特征[20]；Sun等使用DDF-SBM(Directional Distance Function and Slacks-Based Measure)模型測算中國285個城市的包容性綠色增長水平[21]；Chen等結合Super-EBM(Epsilon-Based Measure)模型和MML(Metafrontier Malmquist-Luenberger)指數(shù)評價長江經(jīng)濟帶的包容性綠色增長水平[22]。

綜上，現(xiàn)有研究多基于屬性數(shù)據(jù)，采用傳統(tǒng)空間計量方法從地理鄰近角度刻畫包容性綠色增長的空間分異規(guī)律，對包容性綠色增長水平的空間關聯(lián)網(wǎng)絡分析不足，未準確識別各區(qū)域在網(wǎng)絡中的角色定位；另一方面，已有研究雖通過空間自相關討論區(qū)域間包容性綠色增長的關聯(lián)特征，但對空間關聯(lián)網(wǎng)絡的溢出效應及關聯(lián)網(wǎng)絡的形成機制分析尚待強化。此外，隨著區(qū)域間要素流動與信息交流日益密切，區(qū)域間空間關聯(lián)關系也逐漸呈現(xiàn)復雜網(wǎng)絡結構形態(tài)，社會網(wǎng)絡分析方法能有效克服屬性數(shù)據(jù)的局限性，依托關系數(shù)據(jù)進行全局分析[23，24]，目前已得到廣泛應用[25]，但在包容性綠色增長方面的實證案例并不多見。基于此，本文在測度環(huán)渤海地區(qū)44個地級市2005-2018年包容性綠色增長效率及其空間關聯(lián)關系的基礎上，利用社會網(wǎng)絡分析與二次指派程序(Quadratic Assignment Procedure，QAP)方法深入剖析環(huán)渤海地區(qū)包容性綠色增長效率的空間關聯(lián)網(wǎng)絡結構及其影響因素，為環(huán)渤海地區(qū)跨區(qū)域協(xié)同發(fā)展提供參考。

1 研究方法與評價指標體系

1.1 研究方法

1.1.1 考慮非期望產(chǎn)出的Super-EBM模型 EBM模型是兼容徑向和非徑向的混合距離函數(shù)測度方法，既能測算出目標值與實際值間的改進比例，又能反映各投入或產(chǎn)出非徑向的差距，有效提高了測度結果的準確性[26]。為對有效決策單元(DMU)的效率表現(xiàn)進行排序和區(qū)分，將Andersen建立的超效率 DEA 模型與EBM模型相結合，構建Super-EBM模型。此外，將非期望產(chǎn)出納入測量系統(tǒng)，定義考慮非期望產(chǎn)出、非導向、規(guī)模報酬不變的Super-EBM模型，具體公式參考文獻[27]。

1.1.2 修正引力模型 本文參考已有研究成果，通過構建修正引力模型確定環(huán)渤海各地區(qū)包容性綠色增長效率的空間關聯(lián)關系[28]，公式為：

(1)

式中：Fab為地區(qū)a、b間包容性綠色增長效率的空間關聯(lián)強度；IGGEa和IGGEb分別為地區(qū)a、b的效率值；Dab為地區(qū)a、b間的球面距離；ga和gb分別為地區(qū)a、b的生產(chǎn)總值；Kab為地區(qū)a對Fab的貢獻率。

1.1.3 社會網(wǎng)絡分析方法 該方法是一種針對關系數(shù)據(jù)運用圖論工具和代數(shù)模型探索關系結構對結構組分或整體影響的跨學科分析方法[28]。本文采用該方法從整體網(wǎng)絡結構特征、個體網(wǎng)絡結構特征與塊模型三方面分析包容性綠色增長效率空間關聯(lián)網(wǎng)絡。其中，通過網(wǎng)絡密度、網(wǎng)絡關聯(lián)度、網(wǎng)絡等級度和網(wǎng)絡效率表征整體網(wǎng)絡結構特征，通過度數(shù)中心度、中介中心度衡量個體網(wǎng)絡結構特征；通過塊模型分析網(wǎng)絡內(nèi)部結構狀態(tài)及各地區(qū)在板塊中的位置和作用，根據(jù)已有研究，將板塊劃分為“凈溢出”“雙向溢出”“凈受益”和“經(jīng)紀人”4種類型。整體網(wǎng)絡結構特征、個體網(wǎng)絡結構特征各指標說明[29]見表1，計算公式如下：

表1 空間關聯(lián)網(wǎng)絡特征主要指標含義[29]Table 1 Explanation of main indicators for characteristics of spatial correlation network

D=L/(N(N-1))

(2)

C=1-2V/(N(N-1))

(3)

H=1-K/max(K)

(4)

E=1-M/max(M)

(5)

CAD=n/(N-1)

(6)

(7)

式中：N為網(wǎng)絡中成員總數(shù)量；L為實際擁有關聯(lián)關系數(shù)；V為網(wǎng)絡中不可達成員對的數(shù)量；K為網(wǎng)絡中對稱可達成員對的數(shù)量；M為網(wǎng)絡中多余線條數(shù)量；n為網(wǎng)絡中與某成員直接關聯(lián)的其他成員數(shù)量；gjk為成員j、k間捷徑數(shù)量，gjk(i)為成員j、k間捷徑經(jīng)過成員i的數(shù)量，bjk(i)為成員i位于j、k間捷徑上的概率，j≠k≠i，且j

1.2 評價指標體系

包容性綠色增長是包容性增長和綠色增長兩種理念的有機結合，借鑒已有概念內(nèi)涵[30]，本文認為包容性綠色增長是以經(jīng)濟增長、社會公平、福利共享、資源集約和環(huán)境友好為目標，進而實現(xiàn)社會、經(jīng)濟和自然三大系統(tǒng)協(xié)調(diào)可持續(xù)發(fā)展的增長模式，最終實現(xiàn)人的全面發(fā)展以及人與自然和諧共生。根據(jù)包容性綠色增長的內(nèi)涵，包容性綠色增長效率是指在生產(chǎn)要素投入不變的情況下，社會—經(jīng)濟—自然綜合系統(tǒng)在生產(chǎn)過程中所獲取的社會—經(jīng)濟—生態(tài)綜合效益與實際投入資源要素的比率，可綜合反映社會、經(jīng)濟、自然3個系統(tǒng)協(xié)調(diào)發(fā)展狀況，進而構建包容性綠色增長效率理論框架(圖1)。依據(jù)包容性綠色增長效率內(nèi)涵框架并參考相關研究[31,32]，根據(jù)指標選取的系統(tǒng)性、科學性和數(shù)據(jù)可獲取性原則，本文將經(jīng)濟、社會、自然要素納入到包容性綠色增長效率測算框架，充分考慮環(huán)境污染與社會不均衡等問題，突出“包容性”和“綠色”理念以構建包容性綠色增長效率評價體系(表2)。

表2 包容性綠色增長效率評價體系Table 2 Evaluation system of inclusive green growth efficiency

圖1 包容性綠色增長效率理論框架Fig.1 Theoretical framework of inclusive green growth efficiency

投入方面：考慮勞動力、資源、資本和技術四方面，分別選取年末就業(yè)人員總數(shù)、全社會用電量、土地面積、全社會固定資產(chǎn)投資總額和科研經(jīng)費支出表示。產(chǎn)出方面：從社會、經(jīng)濟和環(huán)境三方面構建5個一級指標。其中，期望產(chǎn)出包括經(jīng)濟效益、社會正效益和生態(tài)效益，共涵蓋7個二級指標，其中，地區(qū)國內(nèi)生產(chǎn)總值表征經(jīng)濟增長，人均社會消費品零售總額、職工平均工資、普通中小學師生比、萬人床位數(shù)、城鎮(zhèn)醫(yī)療保險覆蓋率分別從居民消費、居民收入、教育服務、醫(yī)療衛(wèi)生、社會保障五方面表征社會公平與福利普惠，建成區(qū)綠化覆蓋率表征生態(tài)保護成果。非期望產(chǎn)出包括社會負效益和環(huán)境污染，其中，社會負效益選取城鄉(xiāng)居民收入比和失業(yè)率表征機會非均等性和分配不合理性，進而體現(xiàn)經(jīng)濟增長過程中的非包容性，環(huán)境污染選取工業(yè)廢水排放量和工業(yè)二氧化硫排放量表征環(huán)境損害程度。為方便計算，對含有多個二級指標的一級指標運用熵值法進行綜合。

2 研究區(qū)與數(shù)據(jù)來源

本文研究區(qū)域包括環(huán)渤海地區(qū)三省(遼寧省、河北省、山東省)兩直轄市(北京市和天津市)共44個研究單元。2018年環(huán)渤海地區(qū)常住人口和生產(chǎn)總值分別占全國的18.46%和20.76%，在人力資本、金融集聚、科技創(chuàng)新和對外開放等方面具有明顯優(yōu)勢。但資源能源消耗趨緊，環(huán)境污染等問題十分嚴峻，其中，電力消耗占全國的 20.29%，工業(yè)煙粉塵排放量占全國的25.07%，傳統(tǒng)粗放的發(fā)展模式亟待改善。

本文研究時段為2005-2018年，包容性綠色增長效率計算所用數(shù)據(jù)來源于2006-2019年《中國城市統(tǒng)計年鑒》《中國城市建設統(tǒng)計年鑒》《中國區(qū)域經(jīng)濟年鑒》《中國環(huán)境統(tǒng)計年鑒》以及各市區(qū)統(tǒng)計年鑒、政府工作報告和國民經(jīng)濟和社會發(fā)展統(tǒng)計公報。由于資本存量數(shù)據(jù)無法直接獲取，借鑒單豪杰[33]的研究運用永續(xù)盤存法測算資本存量；國內(nèi)生產(chǎn)總值、社會消費品零售總額等數(shù)據(jù)根據(jù)各省市年鑒提供的平減指數(shù)調(diào)整為以2005年為基期的不變價格；各地級市之間的球面距離利用MATLAB軟件測算得到。

3 結果分析

3.1 包容性綠色增長效率測度結果

通過MaxDEA 7 Ultra軟件基于考慮非期望產(chǎn)出的Super-EBM模型計算得到2005-2018年環(huán)渤海地區(qū)市域包容性綠色增長效率值(圖2)。整體呈波動上升趨勢，效率均值由研究初期的0.584上升至末期的0.771，提升顯著。各地區(qū)包容性綠色增長效率存在顯著的空間非均衡性，研究初期效率均值排序為北京>遼寧>天津>河北>山東，研究末期為北京>遼寧>天津>山東>河北。

圖2 2005-2018年環(huán)渤海地區(qū)包容性綠色增長效率變化趨勢Fig.2 Change trend of inclusive green growth efficiency in Bohai Rim Region from 2005 to 2018

3.2 空間關聯(lián)網(wǎng)絡結構特征

3.2.1 整體結構特征分析

(1)選取2005年、2012年和2018年截面數(shù)據(jù)，基于修正引力模型，利用Ucinet 6.0可視化工具繪制包容性綠色增長效率空間關聯(lián)網(wǎng)絡圖(圖3)。由圖3可知，環(huán)渤海地區(qū)包容性綠色增長效率呈多線程、多流向的復雜網(wǎng)絡結構形態(tài)，并具有顯著的通達性和非均衡性特征，網(wǎng)絡中無孤立單元，中心—外圍式圈層結構明顯。2005年空間關聯(lián)網(wǎng)絡中北京、天津的關聯(lián)關系數(shù)最多，處于核心位置，大連、沈陽、唐山、煙臺、青島、濟南、淄博和萊蕪處于次級核心位置，低等級關聯(lián)關系數(shù)地區(qū)主要集中在魯西、冀北和遼東等地區(qū)，構成空間網(wǎng)絡中的邊緣區(qū)；2012年核心區(qū)域更加凸顯，大連、青島也進入核心位置；2018年網(wǎng)絡關聯(lián)度略微下降，核心區(qū)域和次級核心區(qū)域數(shù)量均有所減少?？傮w上，研究期內(nèi)空間關聯(lián)網(wǎng)絡穩(wěn)固鎖定在以北京、天津、青島為核心區(qū)，濟南、煙臺、沈陽和大連為次核心區(qū)，魯西、冀北、遼東等地區(qū)為邊緣區(qū)的中心—外圍圈層結構。

圖3 環(huán)渤海地區(qū)包容性綠色增長效率空間關聯(lián)網(wǎng)絡Fig.3 Spatial correlation network of inclusive green growth efficiency in Bohai Rim Region

(2)借助Ucinet 6.0軟件計算得到2005-2018年整體網(wǎng)絡結構特征衡量指標(圖4)。1)2005-2018年關聯(lián)關系數(shù)呈波動下降趨勢，由研究初期的281個降至末期的252個，與此對應，網(wǎng)絡密度值也表現(xiàn)出微弱波動下降特征，由研究初期的0.149降至末期的0.133，兩者下降幅度均較小。值得注意的是，現(xiàn)有關聯(lián)關系數(shù)與地區(qū)間最大可能關聯(lián)關系數(shù)(1 892)相比，空間關聯(lián)的緊密程度仍處較低水平。2)研究期內(nèi)網(wǎng)絡關聯(lián)度為1，進一步佐證了空間關聯(lián)網(wǎng)絡具有顯著通達性和溢出效應。3)2005-2018年網(wǎng)絡等級度由0.934升至0.941，呈略微上升趨勢，表明空間網(wǎng)絡的等級結構較為森嚴；網(wǎng)絡效率由0.734升至0.768，呈小幅上升趨勢，表明空間網(wǎng)絡多重疊加效應有所減弱，網(wǎng)絡穩(wěn)定性有所下降。

圖4 環(huán)渤海地區(qū)包容性綠色增長效率整體網(wǎng)絡特征Fig.4 Overall network features of inclusive green growth efficiency in Bohai Rim Region

3.2.2 個體結構特征分析 通過社會網(wǎng)絡分析方法測算度數(shù)中心度和中介中心度，揭示各地市在包容性綠色增長效率空間網(wǎng)絡中的地位和作用(表3)。1)度數(shù)中心度。2005年、2012年和2018年度數(shù)中心度均值分別為23、23和21，總體呈略微下降趨勢，度數(shù)中心度較大地區(qū)主要為北京、天津、青島、濟南、煙臺、沈陽和大連，這些地區(qū)在關聯(lián)網(wǎng)絡中占據(jù)主導地位，與其他地區(qū)關聯(lián)較密切；度數(shù)中心度較低地區(qū)主要為分布在京津周邊的張家口、滄州、保定、承德、廊坊及沈陽周邊的撫順和遼陽等，在關聯(lián)網(wǎng)絡中處于從屬地位。點出度方面，北京、天津、青島、萊蕪、沈陽和大連等地較高，為溢出主體；點入度方面，青島、煙臺、濟寧、聊城、菏澤、邯鄲和大連較高，是關聯(lián)網(wǎng)絡中的受益主體。進一步分析凈受益關系可知，魯西南、冀南、遼西北大部分地區(qū)受益關系大于溢出關系，在關聯(lián)網(wǎng)絡中處于凈受益地位；而北京、天津、濟南、青島、萊蕪和沈陽在關聯(lián)網(wǎng)絡中溢出較多，對其他地區(qū)的“涓滴效應”較顯著。2)中介中心度。2005年、2012年和2018年中介中心度均值分別為5、4和4，天津、青島、大連和沈陽等城市的中介中心度處于高水平，表明這些地區(qū)在關聯(lián)網(wǎng)絡中能更快速地與其他地區(qū)產(chǎn)生連接，發(fā)揮著“橋梁”中介作用；中介中心度較低地區(qū)有德州、保定、張家口、滄州和承德等，其包容性綠色增長水平較低且多位于省際邊緣位置，故在關聯(lián)網(wǎng)絡中的控制力較弱，處于網(wǎng)絡的邊緣地帶。

表3 2005年、2012年、2018年環(huán)渤海地區(qū)包容性綠色增長效率空間關聯(lián)的網(wǎng)絡中心性分析Table 3 Network centrality analysis of spatial correlation of inclusive green growth efficiency in Bohai Rim Region in 2005,2012 and 2018

3.3 塊模型分析

本文通過塊模型分析包容性綠色增長效率的空間聚類特征。選取2018年數(shù)據(jù)，采用Ucinet 6.0中的CONCOR方法，設置最大切分深度為2、集中標準為0.2，將環(huán)渤海地區(qū)44個城市劃分為4個板塊，據(jù)此繪制板塊關系圖(圖5)：北京、濟南和青島屬于第Ⅰ板塊，淄博、棗莊等19市屬于第Ⅱ板塊，天津、煙臺、唐山、沈陽和大連屬于第Ⅲ板塊，鞍山、撫順等17市屬于第Ⅳ板塊。

圖5 環(huán)渤海地區(qū)包容性綠色增長效率板塊間關聯(lián)關系Fig.5 Relationship between inclusive green growth efficiency sectors in Bohai Rim Region

分析板塊內(nèi)部和板塊間的空間關聯(lián)(表4)可知，2018年關聯(lián)關系總數(shù)為252個，板塊內(nèi)和板塊間關系數(shù)比重分別為19.44%、80.56%，說明包容性綠色增長效率的空間溢出效應以板塊間的溢出為主。具體而言，板塊Ⅰ溢出關系85個，接受其他板塊溢出關系12個，前者明顯多于后者，屬“凈溢出”板塊，該板塊期望內(nèi)部比例顯著大于實際內(nèi)部比例；板塊Ⅱ溢出關系55個，其中板塊內(nèi)部40個，板塊外部15個，該板塊內(nèi)部及向外溢出關系均較多，屬于“雙向溢出”板塊，主要是山東大部分和河北部分地區(qū)；板塊Ⅲ向其他板塊溢出關系71個，接收其他板塊溢出關系42個，期望內(nèi)部關系比例大于實際內(nèi)部關系比例，即該板塊既向外溢出又接收其他板塊溢出，在關聯(lián)網(wǎng)絡中起 “橋梁”和“中介”作用，屬于“經(jīng)紀人”板塊；板塊Ⅳ接收溢出關系數(shù)79，向其他板塊溢出關系34個，即該板塊接收其他板塊的關系數(shù)遠多于向外溢出的關系數(shù)，屬于“凈受益”板塊，主要包括遼寧大部分和河北部分地區(qū)?？傊?，“凈溢出”和“經(jīng)紀人”板塊聚集了京津冀城市群、山東半島城市群和遼中南城市群的經(jīng)濟發(fā)達地區(qū)，而經(jīng)濟水平滯后地區(qū)則集中在“凈受益”和“雙向溢出”板塊。

表4 環(huán)渤海地區(qū)包容性綠色增長效率空間關聯(lián)板塊溢出效應Table 4 Spillover effect of spatial correlation sectors of inclusive green growth efficiency in Bohai Rim Region

從板塊內(nèi)部關系看(圖5)，成員間均存在關聯(lián)關系，“雙向溢出”板塊內(nèi)部關系數(shù)較多，“俱樂部”效應最顯著，進一步表明山東省各地級市間的包容性綠色增長效率的聯(lián)動性較強；其他3個板塊內(nèi)部溢出關系極少，板塊內(nèi)部尚未體現(xiàn)出明顯的“俱樂部”效應。從板塊間關系看，“凈溢出”板塊主要對“凈受益”和“雙向溢出”板塊產(chǎn)生溢出效應，接受兩板塊溢出數(shù)極少，居溢出主體地位，扮演第一發(fā)動機角色；“經(jīng)紀人”板塊主要對“凈受益”和“雙向溢出”板塊產(chǎn)生溢出效應，同時接收兩板塊的溢出也較多，在網(wǎng)絡中具有顯著“中介”作用；“凈受益”和“雙向溢出”板塊間關系為零，“凈溢出”和“經(jīng)紀人”板塊之間關系也較弱?？傮w看，環(huán)渤海地區(qū)板塊溢出效應主要集中在經(jīng)濟發(fā)達板塊對經(jīng)濟落后板塊的溢出效應，梯度溢出特征顯著，經(jīng)濟發(fā)達地區(qū)既扮演著中心行動者角色，也兼具“中介”和“橋梁”作用。

4 空間關聯(lián)網(wǎng)絡結構的影響因素分析

4.1 影響因素選取與模型建立

根據(jù)包容性綠色增長理論的本質要求，并結合當前我國社會經(jīng)濟發(fā)展現(xiàn)實情況，包容性綠色增長效率的空間關聯(lián)關系不僅局限于地理空間距離，還受經(jīng)濟、結構、政府、技術、環(huán)境等因素的影響[34]，這些因素的地區(qū)差異促使各要素在區(qū)域間溢出與流動，加強了地區(qū)間的交流合作，同時也催生了空間關聯(lián)網(wǎng)絡。總之，包容性綠色增長效率的空間關聯(lián)特征是多維要素共同作用的結果。故本文在借鑒已有研究的基礎上，從地理距離、經(jīng)濟發(fā)展、產(chǎn)業(yè)結構、政府干預、工業(yè)化、技術創(chuàng)新、對外開放和環(huán)境規(guī)制八方面選取自變量，以包容性綠色增長效率空間關聯(lián)網(wǎng)絡為因變量構建回歸模型，深入探討環(huán)渤海地區(qū)包容性綠色增長效率空間關聯(lián)網(wǎng)絡及其溢出效應的生成機制，公式如下：

IGGE=F(GE,EC,IS,GOV,INDU,TI,OP,ER)

(8)

式中：IGGE為包容性綠色增長效率空間關聯(lián)矩陣；GE為各地市間的地理鄰近矩陣，用地市間Rook鄰接權重矩陣表征；EC為經(jīng)濟發(fā)展水平差異矩陣，用人均GDP表示；IS為產(chǎn)業(yè)結構差異矩陣，用第三產(chǎn)業(yè)與第二產(chǎn)業(yè)比值表示；GOV為政府行政干預差異矩陣，用財政支出占GDP比重表示；INDU為工業(yè)化水平差異矩陣，用工業(yè)增加值占GDP比重表示；TI為技術創(chuàng)新水平差異矩陣，用專利授權量表示；OP為對外開放差異矩陣，用實際使用外資占GDP的比重表示；ER為環(huán)境規(guī)制差異矩陣，用工業(yè)固體廢棄物綜合利用率表示。

由于回歸模型所用數(shù)據(jù)為關系數(shù)據(jù)，解釋變量間相似性較高，為避免多重共線性帶來的測量誤差，本文采用社會網(wǎng)絡分析中的QAP方法對包容性綠色增長效率空間關聯(lián)網(wǎng)絡的影響因素進行實證考察。

4.2 QAP相關性分析

本文采用 QAP 相關性分析檢驗包容性綠色增長效率空間關聯(lián)網(wǎng)絡與各影響因素間的相關關系，相關系數(shù)的大小反映該變量對包容性綠色增長效率空間關聯(lián)網(wǎng)絡結構的影響程度。選擇5 000次隨機置換得到的包容性綠色增長效率空間關聯(lián)與8個變量進行相關性分析，結果(表5)表明，地理距離、經(jīng)濟發(fā)展、產(chǎn)業(yè)結構、技術創(chuàng)新、對外開放和環(huán)境規(guī)制6個變量對包容性綠色增長效率空間關聯(lián)網(wǎng)絡結構具有重要影響，而政府干預和工業(yè)化對包容性綠色增長效率空間關聯(lián)的作用不明顯。

表5 環(huán)渤海地區(qū)包容性綠色增長效率空間關聯(lián)影響因素的相關性分析Table 5 Correlation analysis of influencing factors on spatial correlation of inclusive green growth efficiency in Bohai Rim Region

4.3 QAP回歸分析

QAP回歸結果如表6所示，選擇5 000次隨機置換，得到調(diào)整后的判定系數(shù)R2為0.293且在1%水平上顯著，說明6個變量可以解釋環(huán)渤海包容性綠色增長效率空間關聯(lián)關系變異的29%。具體而言：1)地理鄰近矩陣的回歸系數(shù)在1%水平上顯著，說明空間鄰近有利于空間關聯(lián)關系的產(chǎn)生，這與李敬等[4]的研究結論一致，區(qū)位鄰近更有利于地區(qū)間各類要素的空間轉移，進而促進區(qū)域間“極化效應”和“涓滴效應”的產(chǎn)生，形成穩(wěn)定的空間關聯(lián)網(wǎng)絡。2)經(jīng)濟發(fā)展水平差異的回歸系數(shù)均顯著為正，整體上表明經(jīng)濟的不均衡發(fā)展影響包容性綠色增長效率的空間關聯(lián)。地區(qū)間經(jīng)濟差異是多生產(chǎn)要素共同作用的結果，差異越大，經(jīng)濟發(fā)達地區(qū)對落后地區(qū)的溢出效應越顯著，從而加速生產(chǎn)要素流動，為包容性綠色增長效率區(qū)域關聯(lián)效應的產(chǎn)生創(chuàng)造客觀條件。3)技術創(chuàng)新水平差異的回歸系數(shù)均顯著為正，表明區(qū)域間技術創(chuàng)新水平的差距對于包容性綠色增長效率的空間關聯(lián)網(wǎng)絡具有促進作用。技術創(chuàng)新水平的差異既有利于高水平區(qū)對低水平區(qū)的輻射帶動作用，也有利于不同技術之間的優(yōu)勢互補，打破技術壁壘，促進地區(qū)間技術關聯(lián)網(wǎng)絡的形成，進而加強區(qū)域間包容性綠色增長效率空間關聯(lián)網(wǎng)絡。4)對外開放差異矩陣的回歸系數(shù)顯著為正，說明地區(qū)間對外開放差異有利于空間關聯(lián)網(wǎng)絡的形成?？赡茉蛟谟?，對外開放水平高的地區(qū)對開放水平低的地區(qū)形成一定的示范帶動作用，進一步促進地區(qū)間在政策、管理和技術等方面加強聯(lián)系，從而增強空間關聯(lián)網(wǎng)絡的復雜性。5)產(chǎn)業(yè)結構差異矩陣的回歸系數(shù)在1%的水平上顯著為負，說明地區(qū)間相似的產(chǎn)業(yè)結構對包容性綠色增長效率空間網(wǎng)絡形成具有顯著影響，產(chǎn)業(yè)結構相似表明地區(qū)間發(fā)展階段相同，同時在一定程度上強化了地區(qū)間在技術和資本上的聯(lián)系，從而有利于空間關聯(lián)網(wǎng)絡的形成。6)環(huán)境規(guī)制差異回歸系數(shù)也顯著為負，但回歸系數(shù)較小，說明相似的環(huán)境規(guī)制水平對包容性綠色增長效率空間網(wǎng)絡形成具有一定的促進作用，相同的環(huán)境規(guī)制水平有利于地區(qū)間協(xié)同環(huán)境治理并精準施策，進一步提升環(huán)境治理效率和能力，促進地區(qū)間包容性綠色增長效率共同提升。

表6 環(huán)渤海地區(qū)包容性綠色增長效率空間關聯(lián)影響因素的回歸結果Table 6 Regression results of influencing factors on spatial correlation of inclusive green growth efficiency in Bohai Rim Region

5 結論與政策啟示

本文通過修正引力模型、社會網(wǎng)絡分析方法對環(huán)渤海地區(qū)包容性綠色增長效率的空間關聯(lián)網(wǎng)絡特征進行研究，并對其影響因素進行實證分析，得到如下結論：1)2005-2018年環(huán)渤海地區(qū)包容性綠色增長效率整體呈波動上升趨勢，存在顯著的空間不均衡特征。整體網(wǎng)絡方面，環(huán)渤海地區(qū)包容性綠色增長效率的空間關聯(lián)呈多線程、多流向的復雜網(wǎng)絡結構形態(tài)，并具有顯著的通達性和非均衡性特征，呈現(xiàn)以北京、天津、青島為核心區(qū)，濟南、煙臺、沈陽和大連為次核心區(qū)，魯西、冀北、遼東等地區(qū)為邊緣區(qū)的中心—外圍圈層結構。地區(qū)間包容性綠色增長效率空間溢出效應顯著，但網(wǎng)絡的穩(wěn)定性有所下降，且存在森嚴的等級結構。2)個體網(wǎng)絡方面，北京、天津、青島、濟南、煙臺、大連和沈陽等地區(qū)與其他地區(qū)之間關聯(lián)關系較多，在關聯(lián)網(wǎng)絡中占據(jù)主導地位。其中，北京、天津、青島等作為溢出主體對其他地區(qū)具有顯著“涓滴效應”，天津、青島、大連等在關聯(lián)網(wǎng)絡中能更快速地與其他地區(qū)產(chǎn)生連接，控制作用較強，中介作用顯著；魯西南、冀南、遼西北大部分地區(qū)受益關系大于溢出關系，在關聯(lián)網(wǎng)絡中處于凈受益地位。3)環(huán)渤海包容性綠色增長效率的空間溢出效應以板塊間溢出為主，梯度溢出特征顯著，板塊內(nèi)部成員具有一定“俱樂部”效應。其中，北京、青島和濟南屬于“凈溢出”板塊，在網(wǎng)絡中處于核心地位，發(fā)揮引領作用；天津、大連、煙臺、沈陽和唐山屬于“經(jīng)紀人”板塊，發(fā)揮中介作用；威海、承德、滄州、廊坊以及遼寧大部分地市屬于“凈受益”板塊，在網(wǎng)絡中居邊緣位置，扮演“跟隨者”角色；石家莊、邢臺、保定、邯鄲、張家口和衡水以及山東大部分地市屬于“雙向溢出”板塊。4)地理距離、經(jīng)濟發(fā)展、產(chǎn)業(yè)結構、技術創(chuàng)新、對外開放對包容性綠色增長效率空間關聯(lián)網(wǎng)絡結構具有重要影響,可增強環(huán)渤海地區(qū)包容性綠色增長效率的空間關聯(lián),相似的產(chǎn)業(yè)結構和環(huán)境規(guī)制水平也有利于包容性綠色增長效率空間網(wǎng)絡的形成。

環(huán)渤海地區(qū)包容性綠色增長空間關聯(lián)網(wǎng)絡的研究結論首先能為區(qū)域協(xié)同發(fā)展提供決策支持，空間關聯(lián)網(wǎng)絡的優(yōu)化提升能有效促進區(qū)域間的關聯(lián)和溢出效應，所以應積極探尋促進包容性綠色增長效率空間關聯(lián)的有效途徑，通過增加空間關聯(lián)關系的強度提高地區(qū)間的緊密程度，同時應弱化整體關聯(lián)網(wǎng)絡的等級結構，提升網(wǎng)絡穩(wěn)定性，進一步實現(xiàn)跨區(qū)域的協(xié)同發(fā)展。其次，由于各地區(qū)在關聯(lián)網(wǎng)絡中扮演的角色不同，各板塊在空間網(wǎng)絡中發(fā)揮的作用具有差異，所以要有針對性、因地制宜地制定區(qū)域包容性綠色增長政策。發(fā)揮中心城市對邊緣城市的帶動作用，同時注重各板塊在關聯(lián)網(wǎng)絡中的聯(lián)動和溢出關系，充分發(fā)揮“凈溢出”板塊在資金、人才、技術等先進要素方面的溢出效應，強化“經(jīng)紀人”板塊的傳遞和中介功能，而“凈受益”板塊應更好地轉化和吸收包容性綠色生產(chǎn)要素，進一步促進地區(qū)間要素流動的良性循環(huán)。第三，全面考慮包容性綠色增長效率空間關聯(lián)網(wǎng)絡的影響因素，區(qū)域協(xié)調(diào)發(fā)展不應單純依賴空間鄰近地區(qū)，應從整體網(wǎng)絡結構視角構建跨區(qū)域協(xié)同發(fā)展機制，利用地區(qū)間經(jīng)濟發(fā)展、技術和對外開往水平的顯著差異擴大地區(qū)間溢出效應，利用產(chǎn)業(yè)結構和環(huán)境規(guī)制的相似性發(fā)揮地區(qū)間的同化效應，進一步縮小差距。