環(huán)渤海地區(qū)包容性綠色增長(zhǎng)效率的空間關(guān)聯(lián)網(wǎng)絡(luò)結(jié)構(gòu)及其影響因素

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

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

0 引言

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

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

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

1 研究方法與評(píng)價(jià)指標(biāo)體系

1.1 研究方法

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

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

(1)

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

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

表1 空間關(guān)聯(lián)網(wǎng)絡(luò)特征主要指標(biāo)含義[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)絡(luò)中成員總數(shù)量；L為實(shí)際擁有關(guān)聯(lián)關(guān)系數(shù)；V為網(wǎng)絡(luò)中不可達(dá)成員對(duì)的數(shù)量；K為網(wǎng)絡(luò)中對(duì)稱可達(dá)成員對(duì)的數(shù)量；M為網(wǎng)絡(luò)中多余線條數(shù)量；n為網(wǎng)絡(luò)中與某成員直接關(guān)聯(lián)的其他成員數(shù)量；gjk為成員j、k間捷徑數(shù)量，gjk(i)為成員j、k間捷徑經(jīng)過(guò)成員i的數(shù)量，bjk(i)為成員i位于j、k間捷徑上的概率，j≠k≠i，且j

1.2 評(píng)價(jià)指標(biāo)體系

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

表2 包容性綠色增長(zhǎng)效率評(píng)價(jià)體系Table 2 Evaluation system of inclusive green growth efficiency

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

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

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

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

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

3 結(jié)果分析

3.1 包容性綠色增長(zhǎng)效率測(cè)度結(jié)果

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

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

3.2 空間關(guān)聯(lián)網(wǎng)絡(luò)結(jié)構(gòu)特征

3.2.1 整體結(jié)構(gòu)特征分析

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

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

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

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

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

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

3.3 塊模型分析

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

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

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

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

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

4 空間關(guān)聯(lián)網(wǎng)絡(luò)結(jié)構(gòu)的影響因素分析

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

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

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

(8)

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

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

4.2 QAP相關(guān)性分析

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

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

4.3 QAP回歸分析

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

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

5 結(jié)論與政策啟示

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

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