張苗 劉璇 彭山桂 張玉臻 陳銀蓉 文蘭嬌

摘要：研究目的：探究土地利用碳排放效率的空間關(guān)聯(lián)網(wǎng)絡特征及其形成機制，為我國省域?qū)用嫱恋乩锰寂欧判实膮f(xié)同提升提供參考依據(jù)。研究方法：非參數(shù)SBM效率測算法和社會網(wǎng)絡分析法。研究結(jié)果：（1）2002—2019年中國省域土地利用碳排放效率波浪式下降，呈現(xiàn)“東高西低”的空間分布特征；（2）樣本研究期內(nèi)形成了以京、滬、蘇、浙為核心點的復雜但相對穩(wěn)定的空間關(guān)聯(lián)網(wǎng)絡結(jié)構(gòu)，東部發(fā)達省份在信息資源、權(quán)力、聲望及影響方面優(yōu)勢明顯，省域之間等級結(jié)構(gòu)呈現(xiàn)強化趨勢；（3）凈受益板塊和經(jīng)紀人板塊成員為北京、長三角和珠三角地區(qū)的發(fā)達省份，雙向溢出板塊成員主要位于東北、津冀、黃河中下游地區(qū)，凈溢出板塊成員主要為長江中下游、西北和西南地區(qū)的中西部欠發(fā)達省份；（4）空間關(guān)聯(lián)網(wǎng)絡的形成是地理鄰近、經(jīng)濟發(fā)展、土地利用強度、土地利用結(jié)構(gòu)和城市化等因素差異綜合作用的結(jié)果，資源稟賦差異、市場機制調(diào)節(jié)、政府宏觀調(diào)控和信息技術(shù)進步為主要驅(qū)動機制。研究結(jié)論：優(yōu)化土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡，對于協(xié)同提升土地利用碳排放效率以促進高質(zhì)量發(fā)展和踐行“雙碳”戰(zhàn)略具有重要意義。

關(guān)鍵詞：土地利用；碳排放效率；社會網(wǎng)絡分析法；空間關(guān)聯(lián)網(wǎng)絡

中圖分類號：F301.2 文獻標志碼：A 文章編號：1001-8158（2023）10-0091-11

基金項目：國家自然科學基金項目（42001252，42101272，42271270）；山東省社科規(guī)劃基金青年項目（23DGLJ24）；教育部人文社科基金規(guī)劃項目（22YJA790065）。

中共二十大報告提出“優(yōu)化高質(zhì)量發(fā)展的區(qū)域經(jīng)濟布局和國土空間體系”及“積極穩(wěn)妥推進碳達峰和碳中和”。當前，以化石能源效率提升和清潔利用為標志的技術(shù)性減排潛力不斷縮小[1]，基于系統(tǒng)視角增加減排空間被廣泛提及[2]。土地作為不可流動的要素資源，對資本、勞動力和信息技術(shù)等流動性生產(chǎn)要素產(chǎn)生“引力”作用，帶動了能源的傳遞和流通，在土地利用類型轉(zhuǎn)變及作為經(jīng)濟發(fā)展載體過程中帶來了大量的碳排放[3]。提升碳排放約束下的土地利用效率成為促進高質(zhì)量發(fā)展和踐行“雙碳”戰(zhàn)略新時代生態(tài)要義的重要路徑。

研究表明，不考慮碳排放的土地利用效率測算會高估土地利用的實際生產(chǎn)率[4]。近年來，關(guān)于土地利用效率的測算研究多將碳排放納入考慮，實現(xiàn)了將碳排放作為投入指標[5]到非期望產(chǎn)出指標[6]的轉(zhuǎn)換。在低碳經(jīng)濟發(fā)展和“雙碳”戰(zhàn)略時代背景之下，學者們將只包括碳排放作為非期望產(chǎn)出的“碳排放效率”從包含其他非期望產(chǎn)出的“生態(tài)效率”研究中單獨剝離出來，土地利用碳排放效率也由此作為學術(shù)術(shù)語被采納和使用[6]。已有文獻對土地利用碳排放效率的區(qū)域差異[9-10]、動態(tài)演進[11]及其影響因素[7，9]等展開了系統(tǒng)性研究，論證了經(jīng)濟發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)、城鎮(zhèn)化率等典型因素對土地利用碳排放效率的影響。進一步地，基于土地利用和碳排放都存在較強的空間互動性的事實[11]，學者們采用空間計量模型對不同空間尺度的土地利用碳排放效率進行空間量化與表達，以此探究土地利用碳排放效率的空間集聚特征、溢出效應及其影響因素等[12-14]，并尋求空間層面土地利用碳排放效率的協(xié)同提升對策。

傳統(tǒng)計量經(jīng)濟學視角下的空間溢出效應探究的是地理上“相近”或“相鄰”地區(qū)的空間關(guān)系，對空間權(quán)重矩陣具有較強的依賴性[15]。隨著信息技術(shù)不斷發(fā)展和交通基礎設施不斷完善，地區(qū)間的聯(lián)系更多地突破了地理距離限制。由此得出的政策含義往往限于局部，易導致區(qū)域政策制定靶向性錯誤，難以從整體上把握區(qū)域間土地利用碳排放效率的空間關(guān)聯(lián)特征，不利于土地調(diào)控碳減排目標的全局推進。相較于空間計量總體層面的分析，社會網(wǎng)絡分析法以點對點“關(guān)系數(shù)據(jù)”構(gòu)成的網(wǎng)絡為切入點，刻畫某一研究對象的空間網(wǎng)絡特征[16]，可更為精準地揭示土地利用碳排放效率地區(qū)間的互動規(guī)律及其形成機制。盡管已有研究基于“關(guān)系數(shù)據(jù)”對土地利用碳排放[15，17]和碳排放效率[2，18]的空間關(guān)聯(lián)網(wǎng)絡結(jié)構(gòu)進行了分析，但仍缺乏土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡特征性事實的深入刻畫，特別是空間關(guān)聯(lián)網(wǎng)絡的動態(tài)演進及形成機制分析，難以為區(qū)域土地利用碳排放效率的協(xié)同提升提供決策支持。

綜上，本文采用社會網(wǎng)絡分析法構(gòu)建引力模型刻畫土地利用碳排放效率的空間關(guān)聯(lián)網(wǎng)絡特征，采用QAP回歸分析法分析該特征的影響因素，并建立土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡形成機制的一般分析框架。采用該方法能夠有效掌握處于土地利用碳排放效率網(wǎng)絡中的核心點和邊緣點，明確要素的流出地和流向地，診斷網(wǎng)絡形成的障礙因素，為協(xié)同提升碳排放約束下我國省域?qū)用嫱恋乩眯侍峁┛茖W依據(jù)。本文的邊際貢獻在于：提供了基于社會網(wǎng)絡分析法研究土地利用碳排放效率空間互動關(guān)系的新思路，一方面采用“關(guān)系數(shù)據(jù)”克服了傳統(tǒng)計量分析方法采用“屬性數(shù)據(jù)”的缺陷，另一方面基于長時間序列數(shù)據(jù)刻畫了中國省域土地利用碳排放效率的空間關(guān)聯(lián)網(wǎng)絡演進特征與形成機制。

1 研究方法與數(shù)據(jù)來源

1.1 測度模型與方法

1.1.1 非期望產(chǎn)出模型

SBM（Slack Based Measure）模型在碳排放作為非期望產(chǎn)出的土地利用效率測算中得到了廣泛運用[6]。主要優(yōu)點在于解決了徑向DEA模型對無效率的測量無法包含松弛變量的問題。本文構(gòu)建的土地利用碳排放效率測算的投入指標與產(chǎn)出指標之間主要呈現(xiàn)為非徑向關(guān)系，SBM模型作為非導向模型，同時能將碳排放作為非期望產(chǎn)出納入模型，能夠滿足本文獲得土地利用碳排放效率值的分析目的。故采用TNOE定義的包含非期望產(chǎn)出的SBM模型測算中國省域土地利用碳排放效率，具體公式參見文獻[19]和[20]。

1.1.2 社會網(wǎng)絡分析法

土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡實質(zhì)是在碳排放約束下，“資本”和“勞動力”等“流動性要素資源”在“非流動性要素資源”土地上基于“外界作用力”而導致的要素集散，進而導致地區(qū)間土地利用碳排放效率的空間差異和傳導，最終形成點、線、面相結(jié)合的復雜空間網(wǎng)絡組織。社會網(wǎng)絡分析（Social Network Analysis， SNA）采用圖論和矩陣方法描述關(guān)系模式并探究這些關(guān)系模式對結(jié)構(gòu)中成員或整體的影響[16]，能夠刻畫出“資本”和“勞動力”等流動性稀缺資源在非流動性稀缺資源“土地”上配置效率的空間網(wǎng)絡結(jié)構(gòu)特征，同時又能探究造成這種空間網(wǎng)絡結(jié)構(gòu)特征的影響因素。運算過程為：采用修正的引力模型構(gòu)建二值矩陣，運用中心性分析描述空間網(wǎng)絡特征，通過二次指派程序QAP探究形成空間關(guān)聯(lián)網(wǎng)絡特征的影響因素，并解析形成機制。

（2）網(wǎng)絡特征刻畫指標。采用網(wǎng)絡關(guān)系數(shù)、網(wǎng)絡密度、關(guān)聯(lián)度、網(wǎng)絡效率和網(wǎng)絡等級來刻畫土地利用碳排放效率的整體網(wǎng)絡特征；采用點入度、點出度、度數(shù)中心度、中介中心度和接近中心度來刻畫其個體網(wǎng)絡特征；采用塊模型分析進行分類。具體指標含義參考文獻[15]和[16]。

（3）QAP回歸分析。QAP（Quardratic Assignment Procedure）是一種用來檢驗關(guān)系矩陣之間關(guān)系的非參數(shù)檢驗方法，避免了常規(guī)參數(shù)檢驗法檢驗屬性變量時存在的內(nèi)生性問題。采用QAP回歸分析研究自變量矩陣對土地利用碳排放效率空間關(guān)聯(lián)矩陣（因變量矩陣）的回歸關(guān)系，解讀自變量差異矩陣的顯著性，以此作為構(gòu)建土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡特征形成機制框架的依據(jù)。其基本計量模型為G = f（X），其中X代表自變量矩陣，G為因變量矩陣。

1.2 數(shù)據(jù)來源與指標體系

1.2.1 數(shù)據(jù)來源

本文數(shù)據(jù)主要包括土地利用數(shù)據(jù)、碳排放數(shù)據(jù)和經(jīng)濟社會數(shù)據(jù)三大類。其中：（1）土地利用數(shù)據(jù)為全國土地利用調(diào)查數(shù)據(jù)及在其基礎上的變更數(shù)據(jù)。“一調(diào)”數(shù)據(jù)及其變更數(shù)據(jù)來源于《中國國土資源年鑒》（2003—2009年）；“二調(diào)”數(shù)據(jù)及其變更數(shù)據(jù)來源于《中國國土資源統(tǒng)計年鑒》（2010—2018年），其中2018年的數(shù)據(jù)未公布，故2018年不列入研究時段；“三調(diào)”數(shù)據(jù)及其變更數(shù)據(jù)來源于自然資源部國土調(diào)查成果共享應用服務平臺，更新至2021年。（2）碳排放數(shù)據(jù)采用中國碳核算數(shù)據(jù)庫（CEADs）的省級排放清單的表觀二氧化碳排放量，由于西藏缺乏能源消費數(shù)據(jù)導致無法估算碳排放量，同時港澳臺地區(qū)未被列入核算數(shù)據(jù)庫，由此確定本文的研究對象為西藏和港澳臺地區(qū)除外的30個省級行政區(qū)；同時，因省級碳排放數(shù)據(jù)暫時更新至2019年，由此確定本文的研究時間段上限為2019年。（3）經(jīng)濟社會發(fā)展數(shù)據(jù)主要來源于《中國統(tǒng)計年鑒》（2003—2020年）、《中國環(huán)境統(tǒng)計年鑒》（2005—2020年）。

1.2.2 指標體系

（1）效率測算指標。以獲得碳排放約束下的土地利用效率為導向，本文參考現(xiàn)有文獻[6-7，9-10]構(gòu)建土地利用碳排放效率的投入產(chǎn)出指標體系（表1）。其中，固定資本存量采用張軍的方法[22]，以2000年為基期計算，為保持數(shù)據(jù)的一致性，將第一、二和三產(chǎn)業(yè)的GDP增加值同樣以2000年為基期進行修正。

（2）QAP自變量矩陣指標。區(qū)域差異是影響空間關(guān)聯(lián)緊密程度的重要原因，本文基于土地利用碳排放效率的影響因素差異關(guān)系作為自變量矩陣，檢驗其影響土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡結(jié)構(gòu)的顯著性。自變量矩陣選取如下：①現(xiàn)有文獻中關(guān)于碳排放效率、碳排放的空間網(wǎng)絡分析都將地理鄰近和經(jīng)濟發(fā)展水平差異作為首要影響因素[18]，結(jié)合土地利用碳排放效率空間關(guān)聯(lián)結(jié)構(gòu)的個體網(wǎng)絡特征和塊模型分析，將地理鄰近和經(jīng)濟發(fā)展水平差異作為檢驗的影響因素，分別用省會城市間的地理距離和人均GDP表示；②土地要素作為不可流動的“稀缺資源”，土地資源稟賦和配置結(jié)構(gòu)差異可能是導致地區(qū)間土地利用碳排放效率差異的基礎因素，而資本、勞動力等“流動要素”在土地上的投入差異可能是影響土地利用碳排放效率差異的動態(tài)可控因素，因此將土地利用結(jié)構(gòu)和土地利用強度差異納入分析，分別用建設用地面積占比和地均固定資本存量投入表示[23]；③在對非期望產(chǎn)出約束的條件下，城市化、產(chǎn)業(yè)結(jié)構(gòu)和環(huán)境規(guī)制等作為傳統(tǒng)而被廣泛證實的碳排放影響因素，應被考慮為影響土地利用碳排放效率差異的因素，分別用非農(nóng)人口比率、二三產(chǎn)業(yè)產(chǎn)值比值和單位GDP工業(yè)污染物排放量[24]表示。

2 結(jié)果分析

2.1 土地利用碳排放效率

基于SBM模型采用STATA軟件測算了土地利用碳排放效率值，均值a為某年30個地區(qū)的效率均值，均值b為i地區(qū)2002—2019年的效率均值（圖1），研究期內(nèi)30個地區(qū)的效率均值為0.630。從均值a來看，中國省域土地利用碳排放效率呈現(xiàn)波浪式下降，由2002年的0.723降低到2019年的0.472，意味著在研究期內(nèi)，碳排放作為非期望產(chǎn)出約束和要素投入邊際報酬遞減的事實，沖抵了技術(shù)進步帶來的效率提升，反映出提高碳排放約束下的土地利用效率存在較大的阻力。從均值b來看，上海、江蘇、福建、廣東、北京、浙江、湖南、海南、湖北、山東、河南、四川、安徽、河北和廣西等省份具有較高的土地利用碳排放效率（以大于均值0.630為標準），主要集中于東部地區(qū)經(jīng)濟發(fā)達省份；中西部地區(qū)省份的土地利用碳排放效率值大部分低于平均值，反映出我國省域土地利用碳排放效率空間差異和空間集聚特征并存。

2.2 空間關(guān)聯(lián)網(wǎng)絡特征

2.2.1 整體網(wǎng)絡特征

基于修正的引力模型構(gòu)建中國省域土地利用碳排放效率的空間關(guān)聯(lián)矩陣，采用UCINET 6.740軟件計算得到網(wǎng)絡關(guān)系數(shù)、關(guān)聯(lián)度、網(wǎng)絡密度、網(wǎng)絡等級度和網(wǎng)絡效率等指標（圖2），篇幅所限并考慮到2009年和2019年分別為“二調(diào)”和“三調(diào)”數(shù)據(jù)截至點，以2009年和2019年作為代表年份繪制空間關(guān)聯(lián)網(wǎng)絡圖（圖3）。此外，2002—2019年中國省域土地利用碳排放效率的網(wǎng)絡關(guān)聯(lián)度均為1，意味著空間關(guān)聯(lián)網(wǎng)絡結(jié)構(gòu)具有較好的聯(lián)通性和穩(wěn)健性，所有省市自治區(qū)均處于空間關(guān)聯(lián)網(wǎng)絡當中，相互之間存在明顯的空間關(guān)聯(lián)和溢出效應。

圖2顯示，網(wǎng)絡關(guān)系數(shù)呈現(xiàn)三個階段變化：2002—2007年保持相對低值平穩(wěn)，2008—2014年上升為相對高值波動，2015—2019年逐漸下降為低值波動。在研究期間內(nèi)，網(wǎng)絡密度值最大為2013年的0.189，最小為2019年的0.160，表明中國省域之間存在一定的土地利用碳排放效率關(guān)聯(lián)緊密度；但網(wǎng)絡密度較低且呈現(xiàn)階段動態(tài)變化特征，意味著存在著較大的協(xié)同提升和優(yōu)化空間。網(wǎng)絡等級度在2004年表現(xiàn)為明顯高值，表明在2004年出現(xiàn)了短暫的等級森嚴特征；2015—2019年表現(xiàn)為小幅上升趨勢，意味著省域之間等級結(jié)構(gòu)變得明顯，會阻礙要素在不同地域間的合理流動。網(wǎng)絡效率相對平穩(wěn)，均在0.75之上，反映出空間關(guān)聯(lián)網(wǎng)絡具有較強的穩(wěn)定性，存在較多的溢出渠道。

由圖3來看，中國省域土地利用碳排放效率形成了較為復雜的網(wǎng)絡結(jié)構(gòu)圖，非鄰近省份突破傳統(tǒng)的地理空間限制而產(chǎn)生跨區(qū)聯(lián)動效應，形成以北京、上海、江蘇、浙江等為核心點的核心—邊緣空間結(jié)構(gòu)特征。相比于2009年，2019年網(wǎng)絡關(guān)聯(lián)增加幅度小于減少幅度，意味著省域之間的網(wǎng)絡關(guān)聯(lián)有所減少?？傮w來看，空間關(guān)聯(lián)網(wǎng)絡呈現(xiàn)相對穩(wěn)定性和動態(tài)變化性。

2.2.2 個體網(wǎng)絡特征

采用UCINET 6.740軟件計算得到2002—2019年中國省域土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡中心度數(shù)值，考慮到個體網(wǎng)絡特征的相對穩(wěn)定性，以下分析仍基于2009年和2019年。

（1）度數(shù)中心度。通過局部中心度來區(qū)分受益主體和溢出主體，通過相對中心度來比較大小變化，整理結(jié)果見表2。由局部中心度來看，其一，受益主體分布較為均衡，東中西省份均有涉及，而溢出主體則主要集中在北京、上海、江蘇、浙江和廣東等東部沿海發(fā)達省份，表明在空間關(guān)聯(lián)網(wǎng)絡中經(jīng)濟發(fā)展水平越高越有利于發(fā)揮正向溢出效應，帶動其他省份土地利用碳排放效率的提高；其二，從受益主體和溢出主體的省份數(shù)量比較來看，受益主體省份數(shù)量明顯多于溢出主體數(shù)量，同時凈受益關(guān)系體現(xiàn)出北京、上海、江蘇、浙江等省份點入度明顯高于點出度，這意味著在土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡中多呈現(xiàn)為“虹吸效應”，各省份都在競爭勞動力、資本等流動要素。由相對中心度來看，2019年新晉升為核心主導地位的省份為福建和湖南，天津由2009年排名第三的核心主導地位變成2019年排名倒數(shù)第一的邊緣位置，不再具備核心特征。同時，遼寧、青海、河南和山東4個省份新增為邊緣位置，以西部、中部省份為主的邊緣位置省份逐漸向中部、東部擴散，意味著空間關(guān)聯(lián)網(wǎng)絡逐漸打破了地理鄰近的經(jīng)濟集聚趨勢。相對度數(shù)中心勢指數(shù)反映了2009年和2019年的空間關(guān)聯(lián)網(wǎng)絡向某點集中趨勢的程度，其中2009年為0.431，2019年為0.439，上升2.33%，意味著2019年空間關(guān)聯(lián)網(wǎng)絡圖的整體中心性更高，效率高地呈現(xiàn)向某些省份集中趨勢。

（2）中介中心度。中介中心度數(shù)值越高意味著該點可以更多地控制其他行動者，處于網(wǎng)絡的核心，擁有較大的權(quán)力。表3整理了2009年和2019年高于均值的省份，其中，2019年浙江、天津和重慶不再為網(wǎng)絡核心，新增江蘇為網(wǎng)絡核心，網(wǎng)絡核心越來越集聚于發(fā)達省份。

（3）接近中心度。一個點的入接近中心度越大，說明其他點到這個點越容易；出接近中心度越大，說明這個點到其他點越容易，該點不是網(wǎng)絡的核心點。表3整理了2009年和2019年高于均值的省份，可知入接近中心度高的省份主要集中在東部發(fā)達省份，而出接近中心度高的省份正好相反，主要集中于西部欠發(fā)達省份。

通過上述結(jié)果，有效判斷出了處于空間關(guān)聯(lián)網(wǎng)絡中的核心點和非核心點，其中東部發(fā)達省份在信息資源、權(quán)力、聲望及影響方面較強，而西部欠發(fā)達省份處于弱勢，在網(wǎng)絡核心關(guān)系上呈現(xiàn)東部>中部>西部。時間趨勢來看，該特征在2002—2019年動態(tài)變化相對穩(wěn)定。

2.2.3 塊模型分析

塊模型分析是對社會角色的描述性代數(shù)分析[16]。采用UCINET 6.740軟件，運用CONCOR方法，選擇最大分割密度為2，收斂標準為0.2，將30個省份劃分為4個板塊。參照文獻劃分標準[21]，將各板塊接受關(guān)系比例和內(nèi)部關(guān)系比例劃分為雙向溢出、凈受益、凈溢出和經(jīng)紀人四大板塊，根據(jù)發(fā)出關(guān)系數(shù)和接受關(guān)系數(shù)可判斷板塊與板塊之間的互動關(guān)系。其中，雙向溢出板塊成員與外界成員聯(lián)系相對較少；凈受益板塊成員既接受來自外部成員的關(guān)系，也有來自自身成員的關(guān)系；凈溢出板塊成員與其他位置成員之間的關(guān)系比自己成員之間的關(guān)系多；經(jīng)紀人板塊成員既發(fā)送也接受外部關(guān)系，內(nèi)部成員之間的聯(lián)系比較少[16]。以2009年和2019年為例解讀中國省域土地利用碳排放效率形成的板塊間的空間關(guān)聯(lián)（圖4）。

由圖4來看，2019年，在成員空間分布上，凈受益板塊和經(jīng)紀人板塊成員為北京、長三角和珠三角地區(qū)的發(fā)達省份，雙向溢出板塊成員主要位于東北、津冀、黃河中下游地區(qū)，凈溢出板塊成員主要為長江中下游、西北和西南地區(qū)的中西部欠發(fā)達省份；在板塊間互動上，板塊Ⅰ與其他三個板塊均有互動，接受板塊Ⅳ關(guān)系數(shù)占絕對比重，外部發(fā)出關(guān)系較少，說明北京、上海、江蘇和浙江成為資金、勞動力以及信息技術(shù)等要素的流入地，接受了大量其他板塊特別是中西部地區(qū)省份的要素溢出，成為空間關(guān)聯(lián)網(wǎng)絡中的受益方，也因此稱之為凈受益板塊；對比來看，凈溢出板塊Ⅳ發(fā)出關(guān)系總數(shù)為74，而接受關(guān)系總數(shù)僅為14，意味著大量的中西部欠發(fā)達省份為要素資源流出地，成為空間關(guān)聯(lián)網(wǎng)絡中的虧損方；經(jīng)紀人板塊Ⅱ成員數(shù)量較少且發(fā)送關(guān)系數(shù)和接受關(guān)系數(shù)比重較低，與板塊Ⅰ、Ⅳ存在互動關(guān)系，意味著溢出效應更多的是直接效應而非間接效應；板塊Ⅲ僅與板塊Ⅰ發(fā)生互動，意味著東北和華北地區(qū)省份發(fā)揮溢出效應的范圍相對有限，局限于北京和長三角地區(qū)，并僅接受了該地區(qū)省份的溢出效應。對比2009年，2019年各板塊成員保持相對穩(wěn)定，成員變化集中于凈受益和經(jīng)紀人板塊，尤其是東部地區(qū)省份。其中，天津不再為凈受益板塊成員，浙江和江蘇由經(jīng)紀人板塊調(diào)整為凈受益板塊，福建由凈溢出板塊調(diào)整為經(jīng)紀人板塊，說明東部地區(qū)省份活躍度更高，要素資源的吸引力增加。從板塊間互動聯(lián)系來看，主要表現(xiàn)為凈溢出板塊對經(jīng)紀人板塊發(fā)送關(guān)系數(shù)減少，對凈受益板塊成員發(fā)送關(guān)系數(shù)增加。

綜合來看，在板塊成員數(shù)量上，凈溢出板塊和雙向溢出板塊成員數(shù)量占絕對比重，意味著空間關(guān)聯(lián)網(wǎng)絡中多數(shù)成員以要素輸出為主，而土地利用碳排放效率較高省份發(fā)揮“溢出效應”十分有限，存在大量的流動要素資源向少數(shù)省份集聚的特征；在板塊內(nèi)部關(guān)系上，存在板塊內(nèi)部成員之間缺乏聯(lián)系的特征，如經(jīng)紀人板塊內(nèi)部關(guān)系為0。

2.3 土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡形成機制分析

2.3.1 QAP回歸分析

借助UCINET 6.740軟件采用QAP回歸分析法識別中國省域土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡的影響因素。考慮到該空間網(wǎng)絡關(guān)聯(lián)的相對穩(wěn)定性，及與上述分析中樣本年份的對應性，選取2002年、2009年、 2015年、2019年的因變量矩陣和自變量矩陣進行回歸，設置隨機置換次數(shù)為5 000，得到QAP回歸分析結(jié)果（表4）。由表4可知，調(diào)整后的R2在0.218～0.231之間較小幅度浮動，表明構(gòu)建的7個自變量差異矩陣因素可以穩(wěn)定地解釋中國省域土地利用碳排放效率空間關(guān)聯(lián)效應20%以上。

由表4可知，地理距離和經(jīng)濟發(fā)展水平差異因素在研究期內(nèi)表現(xiàn)出穩(wěn)定的顯著負向影響，其他因素的顯著性呈現(xiàn)階段性變化，具體來看：

（1）地理距離（DL）的回歸系數(shù)顯著為負，意味著地理距離加大則會阻礙土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡的形成。地理距離越近的省份空間傳導和擴散更頻繁，容易建立空間網(wǎng)絡關(guān)聯(lián)關(guān)系，這與現(xiàn)有基于“關(guān)系數(shù)據(jù)”的實證研究結(jié)果相吻合[18]。同時，該結(jié)果印證了塊模型分析中板塊成員以地理臨近抱團為主，地理距離增加會阻礙流動要素資源的傳遞和流通。

（2）經(jīng)濟發(fā)展水平（GDP）的回歸系數(shù)顯著為負，說明省域之間經(jīng)濟發(fā)展水平差異越小越有利于土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡的形成。在市場機制的作用下，資本和勞動力等流動要素更容易發(fā)生在經(jīng)濟水平相近的地區(qū)間。從標準化系數(shù)大小來看，經(jīng)濟發(fā)展水平差異的影響由2002年的0.629降低到2019年的0.510，該降低趨勢反映出中國省域間經(jīng)濟發(fā)展水平差異縮小的事實，且越來越有利于促進土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡的形成。

（3）土地利用結(jié)構(gòu)（LS）和土地利用強度（LI）呈現(xiàn)出交替反向影響，前者表現(xiàn)為“不顯著—顯著負向—不顯著”，后者表現(xiàn)為“顯著為正—不顯著—顯著為正”的階段性變化。在土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡的形成初期，資本和勞動力等流動要素在某些省份快速“集聚”，省域間土地利用強度差異增大，促進了土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡的形成。根據(jù)土地報酬遞增遞減規(guī)律，流動要素資源由“集聚飽和”轉(zhuǎn)為“分散均衡”，先發(fā)獲得流動要素資源的省份逐漸發(fā)揮“溢出效應”；根據(jù)結(jié)構(gòu)決定功能理論，資本和勞動力首先流向土地利用結(jié)構(gòu)差異小的省份。由此，此階段土地利用結(jié)構(gòu)差異發(fā)揮顯著影響，差異越小越有利于土地利用碳排放效率網(wǎng)絡的形成。當“溢出效應”充分發(fā)揮，土地利用結(jié)構(gòu)差異不再具有顯著影響。然而，隨著信息技術(shù)的不斷發(fā)展，新型的流動要素資源再次在特定省份“集聚”，土地利用強度和土地利用結(jié)構(gòu)差異再次交替反向影響。QAP回歸結(jié)果顯示2015年為出現(xiàn)交替的時間點，這與2015年空間關(guān)聯(lián)網(wǎng)絡變化相對應。

（4）城鎮(zhèn)化（CR）進程在研究前期表現(xiàn)為顯著的正向影響，后期不再顯著。研究期內(nèi)，省域之間的城鎮(zhèn)化差異先增大后縮小，在差異增大時期，資本和勞動力快速流向城鎮(zhèn)化進程較高的省份，促進了土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡的形成；在差異縮小時期，要素資源流動速度減緩甚至停滯，對空間關(guān)聯(lián)網(wǎng)絡的形成不再具有顯著作用。

（5）產(chǎn)業(yè)結(jié)構(gòu)（IS）和環(huán)境規(guī)制（ER）因素在研究期內(nèi)未表現(xiàn)出顯著影響。在實證檢驗中，分別采用第三、二產(chǎn)業(yè)GDP增加值比值和第二三產(chǎn)業(yè)GDP增加值之和占比的差異矩陣作為產(chǎn)業(yè)結(jié)構(gòu)的自變量矩陣，但顯著性和系數(shù)大小均沒有明顯改變；環(huán)境規(guī)制因素的加入也沒有明顯改變其他自變量的顯著性和大小。這意味著研究期內(nèi)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整和環(huán)境規(guī)制并未對碳排放這個非期望產(chǎn)出產(chǎn)生明顯的約束作用。事實上，我們期望省域間具有較小的產(chǎn)業(yè)結(jié)構(gòu)差異和環(huán)境規(guī)制差異，以此避免碳排放在省域間的轉(zhuǎn)移；通過產(chǎn)業(yè)結(jié)構(gòu)調(diào)整帶動能源結(jié)構(gòu)調(diào)整、環(huán)境規(guī)制倒逼技術(shù)創(chuàng)新等手段帶來碳排放的實際減少。

2.3.2 形成機制解析

由QAP回歸結(jié)果可知，中國省域土地利用碳排放效率的空間關(guān)聯(lián)受地理距離、經(jīng)濟發(fā)展水平差異等穩(wěn)定性因素，以及土地利用結(jié)構(gòu)、土地利用強度及城鎮(zhèn)化差異等階段性因素的綜合影響，這導致了省域間土地利用碳排放效率存在空間差異。在厘清土地利用碳排放效率內(nèi)涵[7，25]和梳理空間相互作用理論基礎上，借鑒相關(guān)文獻[26-27]將前文所述的“外界”作用力總結(jié)為土地資源稟賦差異、市場機制調(diào)節(jié)、政府宏觀調(diào)控和信息技術(shù)進步。省域間土地資源稟賦的差異是空間網(wǎng)絡形成的“源動力”，在市場機制作用和經(jīng)濟效益最大化的驅(qū)動下，各類資源要素不斷流入土地開發(fā)利用邊際效益高的省份，形成具備率先發(fā)展的優(yōu)勢核心地區(qū)，在空間網(wǎng)絡中表現(xiàn)為土地利用碳排放效率較高的節(jié)點，同時也會剝離出邊緣區(qū)。土地政策是政府宏觀調(diào)控的有效手段，作為土地所有權(quán)的行使方，政府在國家經(jīng)濟高質(zhì)量發(fā)展和“雙碳”戰(zhàn)略的時代背景下，通過控制性分配建設用地指標以及導向性干預土地出讓，對流動要素的流向產(chǎn)生有向和有約束的調(diào)節(jié)控制作用，由此促進了區(qū)域發(fā)展的均衡性，形成多個節(jié)點和眾多連線組成的土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡。加之信息時代發(fā)展和交通基礎設施的不斷完善，為流動要素資源在城市間的聚合與擴散提供了通道，打破了傳統(tǒng)“地理鄰近”區(qū)域聯(lián)系緊密的局面，非“地理鄰近”省域的互動更加頻繁，帶來了土地利用碳排放效率的空間傳導，最終形成了錯綜復雜的土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡。將上述中國省域土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡的形成機制總結(jié)為圖5所示的理論框架。

3 結(jié)論與啟示

本文構(gòu)建了土地利用碳排放效率評價指標體系，采用非期望產(chǎn)出SBM模型測算了2002—2019年中國省域土地利用碳排放效率值，基于社會網(wǎng)絡分析法探討了中國省域土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡的演變特征及其形成機制。主要研究結(jié)論如下：（1）從時序變化來看，2002—2019年中國省域土地利用碳排放效率呈波浪式下降，提高碳排放約束下土地利用效率的緊迫性和重要性并存；從空間分布來看，東部地區(qū)具有較高的土地利用碳排放效率值，中西部地區(qū)效率值較低，空間差異與空間集聚特征并存。（2）從整體網(wǎng)絡特征來看，形成了以北京、上海、江蘇和浙江等為核心點的“核心—邊緣”空間結(jié)構(gòu)特征，網(wǎng)絡密度呈現(xiàn)階段動態(tài)變化，網(wǎng)絡效率相對平穩(wěn)，2005—2019年等級結(jié)構(gòu)更加明顯，存在著較大的協(xié)同提升和優(yōu)化空間。從個體網(wǎng)絡特征來看，2002—2019年動態(tài)變化相對穩(wěn)定，受益主體分布均衡且廣泛，溢出主體集中于北上廣和江浙等東部發(fā)達省份，同時凈受益關(guān)系明顯，空間關(guān)聯(lián)網(wǎng)絡中主要呈現(xiàn)“虹吸效應”；東部發(fā)達省份在信息資源、權(quán)力、聲望及影響方面優(yōu)勢明顯，在網(wǎng)絡核心關(guān)系呈現(xiàn)上東部>中部>西部。（3）塊模型結(jié)果顯示，凈受益板塊和經(jīng)紀人板塊成員為北京、長三角和珠三角地區(qū)的發(fā)達省份，雙向溢出板塊成員主要位于東北、津冀、黃河中下游地區(qū)，凈溢出板塊成員主要為長江中下游、西北和西南地區(qū)的中西部欠發(fā)達省份。中國省域土地利用碳排放效率的空間關(guān)聯(lián)網(wǎng)絡中，多數(shù)成員以要素輸出為主，表現(xiàn)為大量的流動要素資源向少數(shù)省份集聚的特征；省際間溢出效應遠小于區(qū)域間溢出效應，且主要表現(xiàn)為直接效應，東北和華北地區(qū)省份發(fā)揮溢出效應的范圍最小。（4）從形成機制來看，中國省域土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡的形成和演變受資源稟賦差異、市場機制調(diào)節(jié)、政府宏觀調(diào)控和信息技術(shù)進步4種機制的影響。QAP回歸結(jié)果表明，地理距離和經(jīng)濟發(fā)展水平差異越大越不利于空間關(guān)聯(lián)網(wǎng)絡的形成，土地利用結(jié)構(gòu)和土地利用強度差異在空間關(guān)聯(lián)網(wǎng)絡形成中發(fā)揮交替反向影響，城鎮(zhèn)化差異呈現(xiàn)為階段性顯著正向影響，產(chǎn)業(yè)結(jié)構(gòu)和環(huán)境規(guī)制差異影響不明顯。

上述研究結(jié)論為中國省域土地利用碳排放效率協(xié)同提升帶來啟示如下：（1）持續(xù)提升并穩(wěn)定東部省份的土地利用碳排放效率，促進其積極發(fā)揮示范引領作用，強化“溢出高地”功能。經(jīng)濟發(fā)展戰(zhàn)略的變化伴隨著生產(chǎn)要素的集聚和擴散，在碳排放約束下土地利用效率持續(xù)降低的情境下，必須充分發(fā)揮政府宏觀調(diào)控和價格、供求、競爭等市場機制的作用，保證勞動力、資本、技術(shù)、信息等流動性要素在北上廣和江浙等發(fā)達省份上的充足投入，提升并穩(wěn)定上述核心節(jié)點省份的土地利用碳排放效率，以此促進輻射效應和示范效應的形成。（2）打通省際間要素流動通道、打破區(qū)域間條隔，充分發(fā)揮市場機制推動要素改革和合理流動的作用。省域之間要素流動困難是土地利用碳排放效率協(xié)同提升的主要阻礙，而空間關(guān)聯(lián)網(wǎng)絡的緊密性和穩(wěn)定性是其重要保障。對于所屬同一板塊的省份，重點在于建立高質(zhì)量流動要素的共享通道，促進單項功能的轉(zhuǎn)變，如推動以受益為主的北京、上海、江蘇、浙江板塊和廣東、福建的經(jīng)紀人板塊向雙向溢出板塊轉(zhuǎn)型，體現(xiàn)發(fā)達省份的吸引力和輻射力；促進中部地區(qū)部分省份向經(jīng)紀人板塊轉(zhuǎn)變，特別是處于長江經(jīng)濟帶的省份如湖北、湖南和安徽等，發(fā)揮“橋梁”和“中介”作用，傳遞東部對中西部地區(qū)的間接溢出效應；充分利用板塊間所形成的“勢能差”，促進“勢能高”省份的資本、技術(shù)和人才的有序流出，同時加強“勢能低”省份土地的高質(zhì)量供給，并建立產(chǎn)業(yè)轉(zhuǎn)移門檻以減少碳排放的省際間轉(zhuǎn)移。（3）以國家“雙碳”戰(zhàn)略和高質(zhì)量發(fā)展為指引，促成區(qū)域一體化協(xié)同穩(wěn)定發(fā)展，以土地利用碳排放效率空間關(guān)聯(lián)網(wǎng)絡的影響因素為依據(jù)，重點關(guān)注地理距離鄰近、經(jīng)濟發(fā)展水平差異小、處于土地利用結(jié)構(gòu)和強度差異同頻次的省份間的土地利用碳排放效率的協(xié)作提升，規(guī)劃布局上述省份間的交通基礎設施，打通各要素流動通道，以此推動土地利用低碳利用的區(qū)域均衡與協(xié)調(diào)發(fā)展。

本文嘗試將社會網(wǎng)絡分析法應用在土地利用碳排放效率的空間互動關(guān)系研究上，但仍存進一步拓展空間：一是土地利用碳排放效率值應該追求更加精確的排序，如采用超效率SBM模型以克服效率值截尾的不足；二是社會網(wǎng)絡分析法中的整體網(wǎng)絡特征、個體網(wǎng)絡特征、塊模型等已在各領域得到了廣泛運用，需探索該方法的深入應用，如“凝聚子群”“小世界”等，以此探究事物更加全面的空間互動關(guān)系。

參考文獻（References）：

[1] 趙榮欽，黃賢金，鄖文聚，等.碳達峰碳中和目標下自然資源管理領域的關(guān)鍵問題[J] .自然資源學報，2022，37（5）：1123 - 1136.

[2] SHEN W R， LIANG H W， DONG L， et al. Synergistic CO2 reduction effects in Chinese urban agglomerations： perspectives from social network analysis[J] . Science of the Total Environment， 2021， 798. doi： 10.1016/ j.scitotenv.2021.149352.

[3] CHUAI X W， HUANG X J， QI X X， et al. A preliminary study of the carbon emissions reduction effects of land use control[J] . Scientific Reports， 2016， 6. doi： 10.1038/srep36901.

[4] YANG B， CHEN X， WANG Z Q， et al. Analyzing land use structure efficiency with carbon emissions： a case study in the middle reaches of the Yangtze River， China[J] . Journal of Cleaner Production， 2020， 274. doi： 10.1016/ j.jclepro.2020.123076.

[5] 游和遠，吳次芳.土地利用的碳排放效率及其低碳優(yōu)化：基于能源消耗的視角[J] .自然資源學報，2010， 25（11）：1875 - 1886.

[6] 張苗，甘臣林，陳銀蓉，等.中國城市建設用地開發(fā)強度的碳排放效率分析與低碳優(yōu)化[J] .資源科學，2016， 38（2）：265 - 275.

[7] 李國煜，王嘉怡，曹宇，等.碳排放約束下的福建省城鎮(zhèn)建設用地利用效率動態(tài)變化與影響因素[J] .中國土地科學，2020，34（4）：69 - 77.

[8] 羅谷松，李濤.碳排放影響下的中國省域土地利用效率差異動態(tài)變化與影響因素[J] .生態(tài)學報，2019，39（13）：4751 - 4760.

[9] 張詩嘉，劉曉英，杜書云.中原城市群城市土地利用效率時空差異及影響因素研究[J] .長江流域資源與環(huán)境，2021，30（10）：2417 - 2429.

[10] 陳丹玲，盧新海，匡兵.長江中游城市群城市土地利用效率的動態(tài)演進及空間收斂[J] .中國人口·資源與環(huán)境，2018，28（12）：106 - 114.

[11] ZHANG M， LIU X， PENG S G. Effects of urban land intensive use on carbon emissions in China： spatial interaction and multi-mediating effect perspective[J] . Environmental Science and Pollution Research， 2023（30）： 7270 - 7287.

[12] 陳真玲，李金鎧，李靜.中國省域城鎮(zhèn)土地利用效率的影響因素及空間溢出效應[J] .經(jīng)濟經(jīng)緯，2017，34（4）：25 - 30.

[13] 盧新海，唐一峰，匡兵.長江中游城市群城市土地利用效率空間溢出效應研究[J] .長江流域資源與環(huán)境，2018，27（2）：252 - 261.

[14] 于斌斌，蘇宜梅.土地財政如何影響土地利用效率？——基于規(guī)模與技術(shù)視角的動態(tài)空間杜賓模型檢驗[J] .地理研究，2022，41（2）：527 - 545.

[15] YU Z Q， CHEN L Q， TONG H X， et al. Spatial correlations of land-use carbon emissions in the Yangtze River Delta region： a perspective from social network analysis[J] . Ecological Indicators， 2022， 142. doi： 10.1016/j.ecolind.2022.109147.

[16] 劉軍.整體網(wǎng)分析——UCINET軟件使用指南[M] . 3版.上海：格致出版社&上海人民出版社，2019：5 - 385.

[17] 魏燕茹，陳松林.福建省土地利用碳排放空間關(guān)聯(lián)性與碳平衡分區(qū)[J] .生態(tài)學報，2021，41（14）：1 - 11.

[18] 邵海琴，王兆峰.中國交通碳排放效率的空間關(guān)聯(lián)網(wǎng)絡結(jié)構(gòu)及其影響因素[J] .中國人口·資源與環(huán)境，2021，31（4）：32 - 41.

[19] TONE K. A slacks-based measure of efficiency in data envelopment analysis[J] . European Journal of Operational Research， 2001， 130（3）： 498 - 509.

[20] 成剛.數(shù)據(jù)包絡分析方法與MaxDEA軟件[M] .北京：知識產(chǎn)權(quán)出版社，2014：62 - 148.

[21] 趙林，曹乃剛，韓增林，等.中國綠色經(jīng)濟效率空間關(guān)聯(lián)網(wǎng)絡演變特征及影響因素[J] .資源科學，2021，43（10）：1933 - 1946.

[22] 張軍，章元.對中國資本存量K的再估計[J] .經(jīng)濟研究，2003（7）：35 - 43，90.

[23] 張苗，陳銀蓉，程道平，等.土地利用結(jié)構(gòu)和強度變化對碳排放影響分析[J] .資源開發(fā)與市場，2018，34（5）：624 - 628，675.

[24] 王康，李志學，周嘉.環(huán)境規(guī)制對碳排放時空格局演變的作用路徑研究——基于東北三省地級市實證分析[J] .自然資源學報，2020，35（2）：343 - 357.

[25] 張苗.中國土地利用碳排放效率與收斂性研究[D] .武漢：華中農(nóng)業(yè)大學，2017：24 - 26.

[26] 趙林，高曉彤，劉焱序，等.中國包容性綠色效率空間關(guān)聯(lián)網(wǎng)絡結(jié)構(gòu)演變特征分析[J] .經(jīng)濟地理，2021，41（9）：69 - 78，90.

[27] 吉雪強，劉慧敏，張躍松.中國省際土地利用碳排放空間關(guān)聯(lián)網(wǎng)絡結(jié)構(gòu)演化及驅(qū)動因素[J] .經(jīng)濟地理，2023，43（2）：190 - 200.

Evolution Characteristics and Formation Mechanism of Spatial Correlation Network of Provincial Land Use Carbon Emission Efficiency in China

ZHANG Miao1， LIU Xuan2， PENG Shangui1， ZHANG Yuzhen3，4， CHEN Yinrong5， WEN Lanjiao5

（1. School of Economics and Management， Shandong Agricultural University， Taian 271018， China； 2. College of Management and Economics， Tianjin University， Tianjin 300072， China； 3. College of Urban and Environmental Sciences， Peking University， Beijing 100871， China； 4. Peking University-Lincoln Institute Center for Urban Development and Land Policy， Beijing 100871， China； 5. School of Public Administration， Huazhong Agricultural University， Wuhan 430070， China）

Abstract： The purpose of this study is to investigate the evolution characteristics and formation mechanism of the spatial correlation network of land use carbon emission efficiency， to provide a reference for the collaborative improvement in land use carbon emission efficiency. The research methods include non-parametric SBM efficiency model and social network analysis. The results show that： 1） provincial land use carbon emission efficiency in China declined in a wavy pattern from 2002 to 2019. It presented the spatial distribution characteristics of “high in the east and low in the west”. 2） During the research period， it formed a complex but relatively stable spatial correlation network structure with Beijing， Shanghai， Jiangsu and Zhejiang as the core points. The developed eastern provinces had obvious advantages in information resources， power， prestige and influence. The hierarchical structure between provinces exhibited a strengthening trend. 3） The results of block model analysis indicated that Beijing and the developed provinces in the Yangtze River Delta and the Pearl River Delta belonged to the net benefit plate and broker plate. The Northeast region， Tianjin and Hebei region， the middle and lower reaches of the Yellow River region constituted two-way overflow plates. The central and western underdeveloped provinces located in the middle and lower reaches of the Yangtze River， the northwestern and southwestern regions were net overflow plates. 4） The formation of spatial correlation network stemmed from the comprehensive effect of geographical proximity， economic development， land use intensity， land use structure， urbanization， and other factors. Resource endowment differences， market mechanism adjustment， government macrocontrol and information technology progress were the main driving mechanisms. In conclusion， the optimization of spatial correlation network is of great significance for the collaborative improvement in land use carbon emission efficiency. It can promote high-quality development and fulfill the “dual-carbon” goals strategy.

Key words： land use； carbon emission efficiency； social network analysis； spatial correlation network

（本文責編：陳美景）