程珍旎，蔣貴國(guó)，王 濤，董威威，鐘孟君

（1.四川師范大學(xué) 地理與資源科學(xué)學(xué)院，四川 成都 610066；2.四川師范大學(xué) 西南土地資源評(píng)價(jià)與監(jiān)測(cè)教育部重點(diǎn)實(shí)驗(yàn)室，四川 成都 610066）

0 引言

國(guó)土空間規(guī)劃是我國(guó)國(guó)土空間資源開發(fā)與保護(hù)的整體謀劃和頂層設(shè)計(jì)，是優(yōu)化空間格 局、統(tǒng)籌資源配置的總體方案，對(duì)實(shí)現(xiàn)區(qū)域協(xié)調(diào)發(fā)展和城鄉(xiāng)融合發(fā)展具有重要意義［1］。位于國(guó)家構(gòu)建“兩橫三縱”空間開發(fā)格局的主軸線上的成渝地區(qū)雙城經(jīng)濟(jì)圈，是促進(jìn)東中西部協(xié)同互動(dòng)發(fā)展的重要力量和輻射帶動(dòng)西部地區(qū)發(fā)展的重要地區(qū)，是我國(guó)面向中亞、西亞、東南亞等國(guó)家和地區(qū)的開放前沿。作為全國(guó)經(jīng)濟(jì)高質(zhì)量發(fā)展的新晉“增長(zhǎng)極”，成渝地區(qū)雙城經(jīng)濟(jì)圈發(fā)展水平雖與粵港澳大灣區(qū)、長(zhǎng)三角城市群、京津冀城市群有較大差距，但其城市產(chǎn)業(yè)發(fā)展各具優(yōu)勢(shì)和梯度發(fā)展態(tài)勢(shì)，為各產(chǎn)業(yè)有序轉(zhuǎn)移和承接、推動(dòng)一體化發(fā)展、提升整體競(jìng)爭(zhēng)力提供了方向指引［2］。中共中央、國(guó)務(wù)院印發(fā)的《成渝地區(qū)雙城經(jīng)濟(jì)圈建設(shè)規(guī)劃綱要》提出應(yīng)把握要素流動(dòng)和產(chǎn)業(yè)分工規(guī)律，圍繞重慶主城和成都培育現(xiàn)代化都市圈，帶動(dòng)中心城市周邊市地和區(qū)縣加快發(fā)展。2023 年，四川省發(fā)布的《成德眉資同城化發(fā)展暨成都都市圈建設(shè)成長(zhǎng)期三年行動(dòng)計(jì)劃（2023—2025 年）》確定了同城化發(fā)展成長(zhǎng)期的總體思路：以成渝地區(qū)雙城經(jīng)濟(jì)圈建設(shè)為總牽引，認(rèn)真落實(shí)《成都都市圈發(fā)展規(guī)劃》持續(xù)深化成德眉資同城化發(fā)展，全面推進(jìn)國(guó)土空間布局優(yōu)化、基礎(chǔ)設(shè)施互聯(lián)互通、創(chuàng)新協(xié)同產(chǎn)業(yè)協(xié)作、對(duì)外開放協(xié)同共進(jìn)、公共服務(wù)便利共享、生態(tài)環(huán)境共保共治。全面優(yōu)化成都都市圈土地利用結(jié)構(gòu)、提升土地利用效率，對(duì)其作為推動(dòng)成渝地區(qū)雙城經(jīng)濟(jì)圈建設(shè)的先手棋、服務(wù)成渝地區(qū)雙城經(jīng)濟(jì)圈建設(shè)的支撐性工程具有重要作用。

隨著對(duì)于土地需求的不斷提高，國(guó)內(nèi)外學(xué)者在土地利用方面做了大量研究。在土地利用結(jié)構(gòu)方面，Kgpk Weerakoon運(yùn)用基于GIS 的景觀格局研究城市土地利用結(jié)構(gòu)變化，并對(duì)科倫坡區(qū)土地利用結(jié)構(gòu)變化模式進(jìn)行識(shí)別［3］；李小康研究了土地利用結(jié)構(gòu)變化對(duì)碳排放變化的影響及機(jī)理，得出建設(shè)用地變化是影響凈碳排放量變化的主要原因［4］。同時(shí)，學(xué)者們也關(guān)注土地效益的提高。在土地利用效率方面，Gianni等對(duì)不同規(guī)模城市的土地利用效率進(jìn)行空間計(jì)量分析，強(qiáng)調(diào)了規(guī)劃政策對(duì)土地利用效率的作用［5］；朱從謀對(duì)長(zhǎng)三角城市群建設(shè)用地時(shí)空關(guān)聯(lián)性及利用效益進(jìn)行研究，提出分區(qū)制定差異化調(diào)控策略以促進(jìn)土地利用效率提升［6］；匡兵探討了土地利用結(jié)構(gòu)變化與土地利用績(jī)效的關(guān)聯(lián)度，認(rèn)為城鎮(zhèn)村及工礦用地與土地利用效率關(guān)聯(lián)性較為顯著［7］。同時(shí)，近年來，隨著國(guó)土空間規(guī)劃的逐步開展，學(xué)者們也逐漸將視角移向都市圈中的空間作用效應(yīng)，原女原等基于空間溢出效應(yīng)視角研究發(fā)現(xiàn)高技術(shù)產(chǎn)業(yè)集聚能推動(dòng)城市群協(xié)同創(chuàng)新［8］。

綜上所述，學(xué)者們關(guān)于土地利用結(jié)構(gòu)的合理性與時(shí)空演變特征進(jìn)行了大量的研究，并將研究落實(shí)到土地利用結(jié)構(gòu)的優(yōu)化方面，從而對(duì)土地資源進(jìn)行合理配置，以此尋求土地利用效率最大的土地利用方式。國(guó)內(nèi)外學(xué)者多從某一類用地或是社會(huì)經(jīng)濟(jì)對(duì)土地利用效率的影響以及從經(jīng)濟(jì)學(xué)、生態(tài)環(huán)境等交叉學(xué)科視野展開研究，而對(duì)土地利用結(jié)構(gòu)與效率二者的直接關(guān)系研究，尤其是土地利用效率對(duì)土地利用結(jié)構(gòu)的影響研究上還需進(jìn)一步增強(qiáng)。本研究選擇我國(guó)西部重要經(jīng)濟(jì)中心——成都都市圈，以其35 個(gè)縣（市、區(qū)）為研究單位，基于信息熵法、非期望產(chǎn)出模型及Malmquist指數(shù)等方法對(duì)其土地利用結(jié)構(gòu)與土地利用效率及變化特征進(jìn)行研究，在此基礎(chǔ)上對(duì)二者的交互影響程度進(jìn)行分析，并以空間計(jì)量模型為基礎(chǔ)研究其空間效應(yīng)，分析土地利用結(jié)構(gòu)與效率之間的交互影響關(guān)系，實(shí)現(xiàn)兩者的相互促進(jìn)。

1 研究區(qū)概況與數(shù)據(jù)來源

1.1 研究區(qū)概況

成都都市圈位于29°24'—31°41'E，102°51'—105°45'N 之間，地處成都平原經(jīng)濟(jì)區(qū)內(nèi)圈，是以成都市為中心，由與成都市聯(lián)系緊密的德陽(yáng)市、眉山市、資陽(yáng)市共同組成的都市圈，是成渝地區(qū)雙城經(jīng)濟(jì)圈的重要組成部分（圖1）。成都都市圈的GDP 總量較大，且在研究的4 個(gè)時(shí)間段內(nèi)呈現(xiàn)快速上升趨勢(shì)，2005—2019 年間由3 194.39 億元提高到了21 506.5億元（圖2）。其中成都市、德陽(yáng)市、眉山市、資陽(yáng)市四市的生產(chǎn)總值差異明顯，成都市占據(jù)了絕大部分比重，成都市的地區(qū)國(guó)民生產(chǎn)總值由2005年的2 294.51 億元上升到了2019 年的17 012.5 億元，其占總體比重在2005—2019 年間由71%提高到了79%，地區(qū)間的社會(huì)經(jīng)濟(jì)差異充分體現(xiàn)了成都市在成都都市圈中的中心地位。

圖1 成都都市圈概況Figure 1 Overview of Chengdu Metropolitan Area

圖2 成都都市圈地區(qū)生產(chǎn)總值Figure 2 Chengdu Metropolitan Area's GDP

成都都市圈擁有十分豐富的旅游資源，地形、氣候、水文、土壤、資源等良好的自然條件，為成都都市圈的發(fā)展奠定了良好的基礎(chǔ)。同時(shí)成都都市圈地處我國(guó)內(nèi)陸腹地，位于亞歐航路中點(diǎn)，是中亞、南亞、東南亞等地區(qū)的重要交匯點(diǎn)，有著廣闊的市場(chǎng)腹地。因而成都都市圈的土地利用結(jié)構(gòu)的改善與其土地利用效率的提高對(duì)于成渝地區(qū)雙城經(jīng)濟(jì)圈、四川省以及全國(guó)土地資源管理有著重要意義。

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

本研究所用的土地利用類型數(shù)據(jù)即耕地、林地、草地、水域、城鄉(xiāng)工礦居民用地以及未利用地六類一級(jí)用地與DEM柵格數(shù)據(jù)均來源于中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心。土地利用效率分析所用的各類社會(huì)經(jīng)濟(jì)指標(biāo)數(shù)據(jù)主要來源于2006 年、2011 年、2016 年及2020 年的《四川統(tǒng)計(jì)年鑒》《成都統(tǒng)計(jì)年鑒》《德陽(yáng)統(tǒng)計(jì)年鑒》《眉山統(tǒng)計(jì)年鑒》《資陽(yáng)統(tǒng)計(jì)年鑒》以及各地區(qū)的《國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展公報(bào)》《政府工作報(bào)告》。相關(guān)矢量數(shù)據(jù)來源于國(guó)家基礎(chǔ)地理信息系統(tǒng)網(wǎng)站。

2 研究方法

2.1 土地利用結(jié)構(gòu)研究

信息熵。在地理研究中，熵是一定區(qū)域內(nèi)土地利用結(jié)構(gòu)有序度的集中體現(xiàn)，信息熵反映土地利用結(jié)構(gòu)的有序度，熵值越大，表明各類用地面積的差距越小，有序度越低，反之則越高。根據(jù)信息熵的定義［9，10］，其表達(dá)式為：

式中：n為區(qū)域內(nèi)的土地利用類型種類數(shù)量，其中i ＝1，2，…，n；Pi為區(qū)域內(nèi)各類型土地利用面積與總面積之比。

泰爾指數(shù)。該指數(shù)以熵概念為基礎(chǔ)，不僅能計(jì)算區(qū)域整體的差異，還能進(jìn)一步分析整體差異中區(qū)域內(nèi)和區(qū)域間的差異情況，并且可以比較不同區(qū)域?qū)φw差異的貢獻(xiàn)程度。計(jì)算公式為：

式中：Tb與Tw分別為區(qū)域間差距和區(qū)域內(nèi)差距；n為研究區(qū)個(gè)體數(shù)量，其指標(biāo)被分為k 個(gè)區(qū)域，則各區(qū)域分別為gk（k ＝1，2，…，k），第k 區(qū)域gk中的個(gè)體數(shù)目為nk，則nk總和為n；yi表示個(gè)體i指標(biāo)占總熵值的比例；yk表示第k 區(qū)域指標(biāo)占總熵值的比例。其中，Tk為第k 組的組內(nèi)差距（k ＝1，2，…，k）。

2.2 土地利用效率指標(biāo)體系構(gòu)建

本研究在對(duì)國(guó)內(nèi)外相關(guān)研究進(jìn)行系統(tǒng)梳理后，借鑒相關(guān)成果［11-15］，并充分考慮成都都市圈的社會(huì)經(jīng)濟(jì)環(huán)境狀況，從投入、期望產(chǎn)出以及非期望產(chǎn)出3個(gè)部分構(gòu)建土地利用效率測(cè)算指標(biāo)體系（表1）。投入指標(biāo)由土地投入、勞動(dòng)力投入以及資金投入構(gòu)成。土地是社會(huì)經(jīng)濟(jì)發(fā)展的基礎(chǔ)，2020 年中華人民共和國(guó)自然資源部提出的《自然資源價(jià)格評(píng)估通則》對(duì)包括土地在內(nèi)的自然資源進(jìn)行價(jià)格估算，表明各種土地利用類型屬于國(guó)家重要資產(chǎn)，具有重要的經(jīng)濟(jì)價(jià)值。因此本研究將土地要素作為評(píng)價(jià)土地利用效率的重要投入指標(biāo)，在綜合考慮成都都市圈城市建設(shè)以及農(nóng)業(yè)發(fā)展的基礎(chǔ)上，選取耕地、林地、草地以及城鄉(xiāng)、工礦、居民用地作為土地投入指標(biāo)；資金投入在土地利用中可以直接影響社會(huì)經(jīng)濟(jì)效應(yīng)產(chǎn)出，選取固定資產(chǎn)投資作為成都都市圈土地利用效率測(cè)評(píng)的資金投入指標(biāo)；勞動(dòng)力是將要素投入轉(zhuǎn)化為各類產(chǎn)出的重要媒介之一，勞動(dòng)力投入則指區(qū)域土地利用過程中投入的勞動(dòng)力數(shù)量，選取成都都市圈全社會(huì)從業(yè)人口數(shù)作為勞動(dòng)力投入指標(biāo)。

表1 土地利用效率指標(biāo)體系Table 1 Land use efficiency index system

本研究選取地區(qū)生產(chǎn)總值與地方財(cái)政收入作為期望產(chǎn)出的衡量指標(biāo)。地區(qū)生產(chǎn)總值是經(jīng)濟(jì)規(guī)模的絕對(duì)量指標(biāo)，可以反映經(jīng)濟(jì)發(fā)展的規(guī)模和水平；地方財(cái)政收入是當(dāng)?shù)卣惺蛊渎毮艿闹匾Ｕ?，也是衡量地區(qū)經(jīng)濟(jì)發(fā)展水平的重要指標(biāo)。

對(duì)于非期望產(chǎn)出指標(biāo)，本研究參考以往文獻(xiàn)中將廢氣、廢水和固體廢棄物作為非期望產(chǎn)出的做法［16］，同時(shí)在綜合考慮數(shù)據(jù)可獲取性的基礎(chǔ)上，選擇碳排放量和年P(guān)M2.5 總量作為非期望產(chǎn)出指標(biāo)，以表征經(jīng)濟(jì)發(fā)展過程中的環(huán)境負(fù)效應(yīng)的影響。

2.3 土地利用效率研究

非期望產(chǎn)出模型。非期望產(chǎn)出SBM—Undesirable超效率模型，解決了傳統(tǒng)模型的不足，克服了使松弛變量問題，并在土地利用效率測(cè)算中加入了非期望產(chǎn)出，綜合考慮了環(huán)境負(fù)效應(yīng)的影響，使研究更加科學(xué)合理［17］。根據(jù)模型原理，非期望產(chǎn)出的SBM—Undesirable模型表達(dá)式如下：

式中：s-、sg、sb分別代表各類投入與產(chǎn)出指標(biāo)的松弛變量；λ為權(quán)重向量；ρ*為目標(biāo)函數(shù)值，且當(dāng)ρ*＝1 時(shí)，即s-＝0、sg＝0、sb＝0 時(shí)決策單元有效，當(dāng)ρ*＜1 時(shí)，即s-、sg、sb三者中至少有1 個(gè)不為零時(shí)，決策單元無效，效率有一定的提升空間［18］。

Malmquist指數(shù)。其原理為，當(dāng)全要素生產(chǎn)率值大于1 時(shí)表示其在某一時(shí)期到下一時(shí)期呈現(xiàn)正增長(zhǎng)，小于1 時(shí)則下降；其分解結(jié)果意為，變化比率大于1 表示該變化比率對(duì)全要素生產(chǎn)率有促進(jìn)作用，變化比率小于1 則會(huì)導(dǎo)致生產(chǎn)率下降。該方法最早于1953 年由Malmquist.S 提出，是在面板數(shù)據(jù)上通過距離函數(shù)比值測(cè)算所確定的全要素生產(chǎn)效率［19］。Fare等將t時(shí)期的技術(shù)當(dāng)作參考值，在t ＋1 時(shí)期和t時(shí)期之間的馬氏全要素生產(chǎn)效率變化指數(shù)可以表示為：

馬氏全要素生產(chǎn)效率指數(shù)一般以技術(shù)變化與效率變化的幾何平均解釋，其表達(dá)式如下：

2.4 交互影響研究

在空間計(jì)量模型中，考慮到被解釋變量與解釋變量的空間聯(lián)系選擇空間杜賓模型。從模型構(gòu)建的角度，空間杜賓模型如下所示：

式中：i表示區(qū)域；t表示時(shí)間；Yit表示被解釋變量向量；Xit表示解釋變量向量；W 表示空間權(quán)重矩陣；WXit反映了鄰居地區(qū)解釋變量對(duì)本地被解釋變量由于空間傳導(dǎo)機(jī)制所產(chǎn)生的影響［20］；β、θ 表示待估參數(shù)；εi表示隨機(jī)誤差項(xiàng)向量；ρ表示空間自回歸系數(shù)。

3 結(jié)果及分析

3.1 土地利用結(jié)構(gòu)結(jié)果分析

為更直觀地體現(xiàn)成都都市圈各縣（市、區(qū)）的土地利用結(jié)構(gòu)信息熵，本研究利用ArcGIS 軟件，將土地利用結(jié)構(gòu)信息熵值數(shù)據(jù)與成都都市圈行政區(qū)劃矢量圖進(jìn)行關(guān)聯(lián)，得出2005—2019 年成都都市圈各縣（市、區(qū)）土地利用結(jié)構(gòu)信息熵的變化特征，結(jié)果如圖3 所示。

圖3 土地利用結(jié)構(gòu)信息熵分布Figure 3 Land use structure information entropy distribution

由圖3 可知，2005—2019 年間成都都市圈土地利用結(jié)構(gòu)信息熵值在空間變化特征上呈現(xiàn)由西部向東部遞減態(tài)勢(shì)。時(shí)序變化特征上具備明顯的圈層結(jié)構(gòu)，在成都市主城區(qū)及其周邊的龍泉驛區(qū)、雙流區(qū)等縣（市、區(qū)）信息熵值呈現(xiàn)逐年下降趨勢(shì)，而周圍城市土地利用結(jié)構(gòu)信息熵呈現(xiàn)上升趨勢(shì)，更外圍地區(qū)土地利用結(jié)構(gòu)信息熵則基本保持平穩(wěn)或有小幅度上升。信息熵值呈現(xiàn)下降趨勢(shì)，是土地利用有序程度增加的體現(xiàn)，說明成都都市圈以成都市為土地利用有序度較高區(qū)，由其向外有序度逐漸降低。進(jìn)一步分析成都都市圈內(nèi)部土地利用信息熵的差異，成都都市圈土地利用結(jié)構(gòu)信息熵泰爾指數(shù)其分解如表2所示。

表2 泰爾指數(shù)及其分解結(jié)果Table 2 Theil index and its decomposition results

由表2 可知，成都市的土地利用結(jié)構(gòu)信息熵在研究時(shí)段內(nèi)逐步降低，而周圍縣（市、區(qū)）土地利用結(jié)構(gòu)信息熵變動(dòng)幅度較小，甚至有所下降。其他三組組內(nèi)差異則呈現(xiàn)逐年下降的趨勢(shì)，表明德陽(yáng)、眉山、資陽(yáng)等地區(qū)域內(nèi)各縣（市、區(qū)）的土地利用結(jié)構(gòu)信息熵差異不斷縮小。上述變化充分體現(xiàn)了成都都市圈在發(fā)展過程中成都市的中心地位，區(qū)域間的差異貢獻(xiàn)度在逐年下降，表明了成都都市圈4 市間的土地利用結(jié)構(gòu)信息熵差距在逐年縮小，其內(nèi)部土地利用結(jié)構(gòu)和諧發(fā)展。

3.2 土地利用效率結(jié)果分析

3.2.1 非期望超效率模型計(jì)算

根據(jù)表1 進(jìn)行成都都市圈土地利用超效率測(cè)算。為進(jìn)一步掌握各區(qū)域土地利用效率的空間分布情況及演變規(guī)律，本研究運(yùn)用ArcGIS軟件構(gòu)建成都都市圈土地利用效率數(shù)據(jù)庫(kù)，通過對(duì)已有成果進(jìn)行參考［21］，并結(jié)合研究區(qū)的實(shí)際土地利用效率情況，采用自然斷點(diǎn)法，以0.6、1、1.2 為劃分標(biāo)準(zhǔn)，將2005、2010、2015、2019 年四期的土地利用效率值劃分為4 個(gè)相對(duì)類型：低效型、中效型、有效型、高效型，并對(duì)其進(jìn)行可視化表達(dá)，其空間分布圖如圖4 所示。

圖4 成都都市圈土地利用超效率空間分布Figure 4 Spatial distribution of land use super-efficiency in Chengdu Metropolitan Area

由圖4 可知，低效型的縣（市、區(qū)）數(shù)量較多，主要分布于研究區(qū)的東南部；中效型分布較為零散，隨年份不同，其分布也有較大差別；有效型主要分布于成都都市圈中部以及西部地區(qū)；2019 年出現(xiàn)的3 個(gè)高效型均位于成都都市圈中心區(qū)域，即成都市主城區(qū)。從土地利用效率各類型空間分布的變化趨勢(shì)來看，有效型與高效型呈現(xiàn)向成都都市圈中心區(qū)域收縮的態(tài)勢(shì)。

3.2.2 土地利用全要素生產(chǎn)率測(cè)算

根據(jù)表1 對(duì)研究區(qū)進(jìn)行不包含非期望產(chǎn)出的土地利用全要素生產(chǎn)率測(cè)算，計(jì)算當(dāng)期對(duì)上一期的生產(chǎn)率變化，得到成都都市圈2005—2010 年、2010—2015年及2015—2019 年3 個(gè)時(shí)段的土地利用全要素生產(chǎn)率（MI）、純技術(shù)效率變化（PEC）、規(guī)模效率變化（SEC）、技術(shù)進(jìn)步指數(shù)（TC）、技術(shù)效率變化指數(shù)（EC）。成都都市圈土地利用全要素生產(chǎn)率指數(shù)測(cè)算結(jié)果如表3 所示。

表3 Malmquist效率指數(shù)及其分解結(jié)果Table 3 Malmquist efficiency index and its decomposition results

由表3 可知，成都都市圈全要素生產(chǎn)率在2005—2010 年間呈現(xiàn)負(fù)成長(zhǎng)，而純技術(shù)效率變化、規(guī)模效率變化以及技術(shù)效率變化指數(shù)均是其負(fù)增長(zhǎng)的根源。而2010—2015 年及2015—2019 年成都都市圈全要素生產(chǎn)率均大于1，其生產(chǎn)率呈現(xiàn)正增長(zhǎng)，且在2015—2019 年間增速有所放緩，表明成都都市圈發(fā)展方式轉(zhuǎn)為高質(zhì)量平穩(wěn)發(fā)展，技術(shù)進(jìn)步指數(shù)正是兩個(gè)時(shí)段全要素生產(chǎn)率正增長(zhǎng)的主要根源。在2010—2015 年更是達(dá)到了1.459。技術(shù)進(jìn)步指數(shù)表示在這一時(shí)期生產(chǎn)系統(tǒng)的生產(chǎn)技術(shù)的創(chuàng)新程度，成都都市圈的技術(shù)進(jìn)步指數(shù)始終大于1，表明在研究時(shí)段內(nèi)土地利用生產(chǎn)技術(shù)的創(chuàng)新始終是土地利用效率提升的關(guān)鍵。

3.3 交互影響結(jié)果分析

3.3.1 空間相關(guān)性分析

本研究采用基于空間鄰接矩陣的全局Moran's I對(duì)土地利用結(jié)構(gòu)信息熵與土地利用效率的空間相關(guān)性進(jìn)行檢驗(yàn)，結(jié)果如表4 所示。土地利用結(jié)構(gòu)信息熵的Moran's I在樣本考察期內(nèi)均顯著為正，且均通過了1%的顯著性水平檢驗(yàn)，同時(shí)其土地利用結(jié)構(gòu)信息熵空間相關(guān)性在逐漸增強(qiáng)；成都都市圈的土地利用效率全局Moran's I指數(shù)均為正，說明成都都市圈土地利用效率存在空間相關(guān)性。以上表明空間計(jì)量模型適用于本研究的分析。

表4 Moran's I檢驗(yàn)結(jié)果Table 4 Moran's I test results

3.3.2 模型識(shí)別與檢驗(yàn)

對(duì)于空間計(jì)量模型的辨別，根據(jù)相關(guān)文獻(xiàn)［22］處理。由表5 可知，本研究選擇時(shí)間、空間雙固定效應(yīng)的空間杜賓模型。

表5 最優(yōu)模型選擇檢驗(yàn)結(jié)果Table 5 Test results for optimal model selection

3.4 模型結(jié)果分析

3.4.1 面板模型計(jì)量初步分析

利用Stata16.0 軟件，代入時(shí)空鄰接權(quán)重矩陣，對(duì)成都都市圈以及都市圈各區(qū)域2005、2010、2015、2019年四期土地利用結(jié)構(gòu)信息熵（x）與土地利用效率（y）進(jìn)行雙固定效應(yīng)的空間杜賓模型分析，結(jié)果如表6 所示。

表6 成都都市圈土地利用結(jié)構(gòu)信息熵與土地利用效率交互影響Table 6 Interactive influence of land use structure information entropy and land use efficiency in Chengdu Metropolitan Area

由表6 可知，成都都市圈土地利用結(jié)構(gòu)信息熵與土地利用效率之間存在顯著的交互影響關(guān)系。從土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率的影響來看，成都都市圈土地利用結(jié)構(gòu)信息熵的空間自回歸系數(shù)的估計(jì)值為0.370，表明研究區(qū)內(nèi)相鄰的縣（市、區(qū)）其土地利用效率會(huì)產(chǎn)生相互影響，即存在正的空間相關(guān)性，本縣（市、區(qū)）的土地利用效率會(huì)受到相鄰縣（市、區(qū)）的土地利用效率變化影響；系數(shù)估計(jì)值為-1.742，表明成都都市圈的土地利用效率會(huì)隨著土地利用結(jié)構(gòu)信息熵的提高而降低；從空間溢出的角度看，其空間溢出效應(yīng)系數(shù)為0.936，表明成都都市圈土地利用結(jié)構(gòu)信息熵存在一定的空間擴(kuò)散效應(yīng)。

從土地利用效率對(duì)土地利用結(jié)構(gòu)信息熵的影響來看，其估計(jì)系數(shù)為-0.163，即土地利用效率每提高1%將會(huì)使土地利用結(jié)構(gòu)信息熵下降0.163%；空間自回歸系數(shù)估計(jì)值為0.566，表明成都都市圈各相鄰縣（市、區(qū)）的土地利用結(jié)構(gòu)信息熵會(huì)相互影響；但其空間溢出效應(yīng)系數(shù)為0.009 且不顯著，表明成都都市圈的土地利用效率不具有明顯的空間擴(kuò)散效應(yīng)。對(duì)比成都都市圈土地利用結(jié)構(gòu)信息熵與土地利用效率之間的相互影響情況，可以看出，土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率的影響要強(qiáng)于土地利用效率對(duì)利用結(jié)構(gòu)信息熵的影響，分別為-1.742 與-0.163，表明二者之間的交互影響以土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率的影響為主。

分析成都都市圈各市計(jì)算結(jié)果（表7—10）可知，成都都市圈各市土地利用結(jié)構(gòu)信息熵與土地利用效率存在較為顯著的交互影響。其中成都市的土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率的空間溢出效應(yīng)系數(shù)為-0.738，表明成都市的土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率的影響在空間上具有一定的收斂性，但并不顯著。土地利用效率對(duì)土地利用結(jié)構(gòu)信息熵的影響系數(shù)為-0.144，同時(shí)其空間溢出效益為正，具有不顯著的空間擴(kuò)散效應(yīng)。成都市二者交互影響空間自回歸系數(shù)均顯著為正，表明成都市的各縣（市、區(qū)）本地區(qū)的土地利用結(jié)構(gòu)信息熵上升或下降會(huì)造成鄰居地區(qū)的土地利用結(jié)構(gòu)信息熵上升或下降，土地利用效率同理。

表7 成都市土地利用結(jié)構(gòu)信息熵與土地利用效率交互影響Table 7 Interactive influence of land use structure information entropy and land use efficiency in Chengdu

表8 德陽(yáng)市土地利用結(jié)構(gòu)信息熵與土地利用效率交互影響Table 8 Interactive influence of land use structure information entropy and land use efficiency in Deyang

表9 眉山市土地利用結(jié)構(gòu)信息熵與土地利用效率交互影響Table 9 Interactive influence of land use structure information entropy and land use efficiency in Meishan

表10 資陽(yáng)市土地利用結(jié)構(gòu)信息熵與土地利用效率交互影響Table 10 Interactive influence of land use structure information entropy and land use efficiency in Ziyang

德陽(yáng)市的土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率的影響系數(shù)與土地利用效率對(duì)土地利用結(jié)構(gòu)信息熵的影響系數(shù)分別為-2.257 與-0.248，同時(shí)，二者空間溢出效應(yīng)系數(shù)分別為0.336、0.105，空間自回歸系數(shù)分別為0.372、0.054，均擁有不顯著的空間擴(kuò)散效應(yīng)與空間相關(guān)性，表明德陽(yáng)市各縣（市、區(qū)）與其鄰近地區(qū)的土地利用結(jié)構(gòu)之間、土地利用效率之間相互影響程度不高。

眉山市與資陽(yáng)市的土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率影響系數(shù)分別為-3.415 與-3.845，土地利用效率對(duì)土地利用結(jié)構(gòu)信息熵的影響系數(shù)分別為-0.15 與-0.126。但資陽(yáng)市、眉山市的土地利用結(jié)構(gòu)信息熵與土地利用效率空間溢出效應(yīng)系數(shù)與空間自回歸系數(shù)均為負(fù)值且不顯著，呈現(xiàn)出一定的空間收斂效應(yīng)以及不顯著的空間負(fù)相關(guān)關(guān)系，表明眉山市、資陽(yáng)市各縣（市、區(qū)）間的關(guān)聯(lián)性不強(qiáng)。

3.4.2 空間效應(yīng)分解結(jié)果分析

由于在空間計(jì)量模型中構(gòu)建了空間權(quán)重矩陣，使得模型具有非線性結(jié)構(gòu)，因此回歸系數(shù)不再反映自變量對(duì)因變量的影響，空間滯后項(xiàng)系數(shù)已經(jīng)不能準(zhǔn)確地反映空間計(jì)量模型的空間溢出效應(yīng)［23］。Le Sage等指出，以求解偏微分的方式給出空間杜賓計(jì)量模型（SDM）的參數(shù)釋義，可以更為有效的檢驗(yàn)變量是否存在空間溢出效應(yīng)［24］。成都都市圈偏微分估計(jì)結(jié)果如表11 所示。

表11 土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率影響分解結(jié)果Table 11 The decomposition results of the influence of land use structure information entropy on land use efficiency

表11 中直接效應(yīng)表示成都都市圈及其各市土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率的直接影響，以及土地利用結(jié)構(gòu)信息熵對(duì)鄰近地區(qū)土地利用效率產(chǎn)生影響后反饋給當(dāng)?shù)赝恋乩眯实挠绊?；間接效應(yīng)表示鄰近地區(qū)的土地利用結(jié)構(gòu)信息熵對(duì)當(dāng)?shù)氐耐恋乩眯十a(chǎn)生的影響；總效應(yīng)表示成都都市圈或各市所有區(qū)域土地利用結(jié)構(gòu)信息熵對(duì)該區(qū)域土地利用效率的影響［25］。從成都都市圈及4 市土地利用結(jié)構(gòu)的信息熵對(duì)土地利用效率影響分解結(jié)果來看，成都都市圈土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率影響的直接效應(yīng)為-1.720，成都市、德陽(yáng)市、眉山市、資陽(yáng)市分別為-1.063、-2.346、-3.409、-3.804，表明成都都市圈與其各市內(nèi)各縣（市、區(qū)）的土地利用結(jié)構(gòu)信息熵和受該地土地利用結(jié)構(gòu)信息熵影響的鄰居地區(qū)土地利用效率的反饋?zhàn)饔?，?duì)該縣（市、區(qū)）的土地利用效率有著顯著的負(fù)向影響。即該縣（市、區(qū)）的土地利用結(jié)構(gòu)信息熵上升，在影響鄰居地區(qū)土地利用效率后，在與其共同作用下將導(dǎo)致該縣（市、區(qū)）的土地利用效率下降。

成都都市圈以及4 個(gè)城市的土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率影響的間接效應(yīng)分別為-0.359、-1.518、-0.908、-1.575、-1.583，其中除了成都市的間接效應(yīng)在10%的顯著性水平下顯著以外，其余地區(qū)均不顯著。表明除成都市以外其余各地區(qū)內(nèi)各縣（市、區(qū)）土地利用效率受鄰居地區(qū)的土地利用結(jié)構(gòu)信息熵的負(fù)向影響，但并不顯著。成都市內(nèi)各縣（市、區(qū)）的土地利用效率，會(huì)受到區(qū)域內(nèi)其他縣（市、區(qū)）土地利用結(jié)構(gòu)信息熵較為顯著的負(fù)向影響。雖然成都都市圈及各市的間接效應(yīng)顯著性不強(qiáng)，但其總效應(yīng)均較為顯著，即成都都市圈與各市各區(qū)域內(nèi)的土地利用結(jié)構(gòu)信息熵對(duì)該區(qū)域的土地利用效率有顯著的負(fù)向影響。成都市相比其余地區(qū)，在間接效應(yīng)上較為顯著，表明成都市內(nèi)各區(qū)域之間土地利用結(jié)構(gòu)與效率存在較為明顯的相互影響關(guān)系。從對(duì)成都市與其余地區(qū)總體的交互影響情況進(jìn)行對(duì)比可以發(fā)現(xiàn)，成都市周圍地區(qū)的土地利用結(jié)構(gòu)變化對(duì)土地利用效率的影響強(qiáng)于成都市，但其各縣（市、區(qū)）的土地利用結(jié)構(gòu)與效率之間的關(guān)聯(lián)性不強(qiáng)，同時(shí)成都市周圍地區(qū)的間接效應(yīng)與總效應(yīng)均較成都市小，且顯著性水平低，表明成都都市圈周圍地區(qū)各區(qū)域間相互影響程度不高。因此，成都都市圈在成都市發(fā)揮其引領(lǐng)作用同時(shí)，周圍地區(qū)應(yīng)加強(qiáng)區(qū)域間的聯(lián)系與交流，強(qiáng)化對(duì)成都市發(fā)展的支撐作用。

4 結(jié)論與政策建議

4.1 結(jié)論

本研究利用成都都市圈4 市35 個(gè)縣（市、區(qū)）2005、2010、2015、2019 年的面板數(shù)據(jù)，在對(duì)土地利用結(jié)構(gòu)信息熵和利用效率定量測(cè)度的基礎(chǔ)上，以成都都市圈整體、各市之間及各市內(nèi)部之間為研究單元，利用空間面板杜賓模型探討了土地利用結(jié)構(gòu)與利用效率的交互影響關(guān)系。主要結(jié)論如下：①成都都市圈土地利用結(jié)構(gòu)信息熵經(jīng)歷了一個(gè)先升后降的變動(dòng)趨勢(shì)，組間差距貢獻(xiàn)度結(jié)果表明各區(qū)域間的差異不斷縮小。其土地利用效率總體呈上升趨勢(shì)，中心區(qū)有效型與高效型土地利用效率逐漸增多，技術(shù)進(jìn)步指數(shù)與環(huán)境負(fù)產(chǎn)出對(duì)土地利用效率影響十分顯著。②成都都市圈土地利用結(jié)構(gòu)信息熵（土地利用效率）存在空間相關(guān)性，其土地利用結(jié)構(gòu)信息熵與土地利用效率存在顯著的交互影響關(guān)系，且土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率的影響要強(qiáng)于土地利用效率對(duì)利用結(jié)構(gòu)信息熵的影響；成都市及各地級(jí)市中發(fā)展水平越低的地區(qū)，土地利用結(jié)構(gòu)信息熵對(duì)土地利用效率的影響越大。③成都都市圈空間溢出效應(yīng)以直接效應(yīng)為主，是變化的內(nèi)因，間接效應(yīng)顯著性水平不高，是變化的外因，表明成都都市圈的被解釋變量受鄰居地區(qū)的解釋變量影響較??；成都都市圈土地利用結(jié)構(gòu)信息熵與土地利用效率的交互影響以本地為主，成都市與周圍地區(qū)交互影響情況差異明顯。

4.2 對(duì)策建議

成都都市圈作為《四川省主體功能區(qū)規(guī)劃》中的城鎮(zhèn)化戰(zhàn)略格局規(guī)劃的重要核心，是四川省經(jīng)濟(jì)的核心區(qū)域和帶動(dòng)西部經(jīng)濟(jì)發(fā)展的重要增長(zhǎng)極，同時(shí)也是成渝地區(qū)雙城經(jīng)濟(jì)圈發(fā)展的動(dòng)力源，有著十分重要的戰(zhàn)略地位。依據(jù)成渝地區(qū)雙城經(jīng)濟(jì)圈以及成都都市圈規(guī)劃，成都都市圈要實(shí)現(xiàn)更高水平發(fā)展、處理好內(nèi)部梯度結(jié)構(gòu)問題，優(yōu)化都市圈內(nèi)土地利用結(jié)構(gòu)與提高土地利用效率是必不可少的一環(huán)。

首先，在土地利用結(jié)構(gòu)方面，在四川省國(guó)土空間規(guī)劃的指導(dǎo)下，成都市周圍土地利用結(jié)構(gòu)信息熵較高地區(qū)以及土地利用結(jié)構(gòu)信息熵不降反增地區(qū)應(yīng)以市場(chǎng)化為主、政府手段為輔的方式來促進(jìn)土地資源的開發(fā)和整理，注重各地類面積的合理分配，減小都市圈內(nèi)部土地利用結(jié)構(gòu)差異。在都市圈土地利用結(jié)構(gòu)有序發(fā)展的同時(shí)，成都市更應(yīng)按照國(guó)土空間規(guī)劃嚴(yán)格用地突出其圈層核心地位，并為周圍地區(qū)帶來相關(guān)受益。

其次，在土地利用效率方面，重視影響土地利用效率因素中的主導(dǎo)因素，即技術(shù)效率變化因素。建議形成以成都市為主體、德陽(yáng)市為次主體、資陽(yáng)市與眉山市為輔的技術(shù)創(chuàng)新區(qū)、資陽(yáng)市與眉山兩市應(yīng)當(dāng)通過加大對(duì)高素質(zhì)人才的培養(yǎng)與引進(jìn)，如采取引進(jìn)高校、提高人才引進(jìn)待遇等措施提高技術(shù)效率，進(jìn)而提高兩市以及成都都市圈整體土地利用效率。

最后，把握土地利用結(jié)構(gòu)與土地利用效率兩者之間的交互影響機(jī)制，利用成都都市圈以土地利用結(jié)構(gòu)對(duì)土地利用效率的影響為主的特點(diǎn)，重點(diǎn)進(jìn)行土地利用結(jié)構(gòu)有序化。在各市進(jìn)行內(nèi)部土地利用結(jié)構(gòu)優(yōu)化的同時(shí)，應(yīng)注意成都都市圈內(nèi)部各市之間的相互影響，構(gòu)建信息共享機(jī)制，通過學(xué)術(shù)、經(jīng)濟(jì)會(huì)議以及社會(huì)工程項(xiàng)目加強(qiáng)區(qū)域間的聯(lián)系與交流。土地資源配置方面，加強(qiáng)與發(fā)展改革、工業(yè)和信息化、工商等部門的協(xié)調(diào)聯(lián)動(dòng)，共同推進(jìn)土地利用結(jié)構(gòu)優(yōu)化，提高土地利用效率，構(gòu)建城鎮(zhèn)發(fā)展帶。眉山市與資陽(yáng)市兩市應(yīng)加強(qiáng)與周圍地區(qū)尤其是與成都市的聯(lián)系與交流，相互促進(jìn)，強(qiáng)化對(duì)成都市的支撐作用，形成“一圈一軸、兩翼三帶”城鎮(zhèn)格局，強(qiáng)化成渝地區(qū)雙城經(jīng)濟(jì)圈的城鎮(zhèn)支撐。

成都都市圈作為成渝地區(qū)雙城經(jīng)濟(jì)圈的重要組成部分，應(yīng)當(dāng)抓住機(jī)遇，進(jìn)行相關(guān)基礎(chǔ)設(shè)施建設(shè)，減少都市圈內(nèi)部要素流動(dòng)路徑差異性，爭(zhēng)取削弱都市圈市域甚至縣域之間的要素流動(dòng)壁壘，加強(qiáng)區(qū)域之間生產(chǎn)要素流動(dòng)與先進(jìn)技術(shù)學(xué)習(xí)，最大程度上發(fā)揮土地利用結(jié)構(gòu)與土地利用效率的優(yōu)化與提高地區(qū)間的溢出與擴(kuò)散效應(yīng)，從而形成成都市引領(lǐng)周圍地區(qū)，周圍地區(qū)支撐成都市，成都都市圈帶動(dòng)成渝地區(qū)雙城經(jīng)濟(jì)圈共同進(jìn)步的發(fā)展格局。