張之鳳, 唐玉萍, 潘洪義, 馬紅菊

（1.四川師范大學(xué)西南土地資源評(píng)價(jià)與監(jiān)測(cè)教育部重點(diǎn)實(shí)驗(yàn)室/地理與資源科學(xué)學(xué)院, 成都 610066；2.四川省農(nóng)業(yè)農(nóng)村廳農(nóng)業(yè)生態(tài)資源保護(hù)中心, 成都 610041）

隨著城鎮(zhèn)化進(jìn)程的推進(jìn), 城鎮(zhèn)規(guī)模不斷擴(kuò)張, 出現(xiàn)耕地持續(xù)被占用的現(xiàn)象, 而城鎮(zhèn)周邊一般耕地質(zhì)量相對(duì)較高, 往往會(huì)引起耕地?cái)?shù)量和質(zhì)量的雙降低。新冠肺炎疫情后糧食安全關(guān)注度增加, 但當(dāng)下區(qū)域經(jīng)濟(jì)高質(zhì)量發(fā)展目標(biāo)實(shí)現(xiàn)缺乏堅(jiān)實(shí)的后備力量, 因此, 城鎮(zhèn)擴(kuò)張占用耕地的效率引起了學(xué)者的關(guān)注[1, 2］。2016 年4 月, 國務(wù)院批準(zhǔn)《成渝城市群發(fā)展規(guī)劃》, 強(qiáng)調(diào)打造以成都為核心具有國際影響力的現(xiàn)代化都市圈, 其中, 成都都市圈成為主體區(qū)域, 成為成渝地區(qū)雙城經(jīng)濟(jì)圈戰(zhàn)略高地[3-5］。隨著該規(guī)劃的落實(shí)必然會(huì)加快相應(yīng)的城鎮(zhèn)化進(jìn)程, 在此期間, 實(shí)現(xiàn)農(nóng)業(yè)和二三產(chǎn)業(yè)的有效銜接, 推進(jìn)現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)的轉(zhuǎn)型升級(jí), 提高城鎮(zhèn)擴(kuò)張占用耕地的效率, 促使投入產(chǎn)出達(dá)到有效均衡, 推動(dòng)成都都市圈產(chǎn)業(yè)協(xié)同發(fā)展。

改革開放后, 城鎮(zhèn)快速擴(kuò)張給當(dāng)時(shí)以農(nóng)業(yè)為主的中國經(jīng)濟(jì)帶來一些負(fù)面影響, 學(xué)者從定性層面上對(duì)城鎮(zhèn)擴(kuò)張占用耕地的問題以及應(yīng)對(duì)政策進(jìn)行探討, 而對(duì)區(qū)域數(shù)量增減和變化率定量分析的研究深度較淺[6-8］。20 世紀(jì)90 年代后國內(nèi)市場(chǎng)經(jīng)濟(jì)快速發(fā)展, 城鎮(zhèn)擴(kuò)張占用耕地力度大小與地方經(jīng)濟(jì)繁榮程度密切相關(guān), 退耦模型、相關(guān)度、脫鉤模型、聯(lián)立方程模型等剖析二者之間關(guān)系的方法得以運(yùn)用, 深入模型定量研究初現(xiàn)端倪[9-12］。隨著遙感技術(shù)、ArcGIS軟件平臺(tái)等的廣泛應(yīng)用, 時(shí)空兩維度變動(dòng)成為當(dāng)前主要研究尺度, 空間相關(guān)性、冷熱點(diǎn)分析、空間加權(quán)矩陣等研究方法成為其時(shí)空變動(dòng)的主要表達(dá)形式[13-15］。同時(shí)為了更有針對(duì)性地解決不同區(qū)域城鎮(zhèn)擴(kuò)張占用耕地產(chǎn)生的問題, 多位學(xué)者采用主成分分析、匯合OLS、地理探測(cè)器等統(tǒng)計(jì)模型方法定量探討驅(qū)動(dòng)城鎮(zhèn)化占用大量耕地的因素[15-17］。近年來, 部分學(xué)者開始利用時(shí)空配置效率和綜合效率分析被占用耕地的利用程度, 以便緩解當(dāng)下經(jīng)濟(jì)發(fā)展與資源保護(hù)的雙重壓力[17, 18］。

在現(xiàn)有城鎮(zhèn)擴(kuò)張占用耕地方向研究中, 采用象限法來分析占用方向變動(dòng)的研究相對(duì)較少, 一般從綜合效率分析其經(jīng)濟(jì)效益, 未能考慮不同區(qū)域產(chǎn)生的規(guī)模和技術(shù)效率差異, 對(duì)其經(jīng)濟(jì)效益進(jìn)行聯(lián)立分析的鮮有報(bào)道。因此, 本研究在刻畫4 時(shí)段成都都市圈城鎮(zhèn)擴(kuò)張占用耕地現(xiàn)狀強(qiáng)度、貢獻(xiàn)度指數(shù)的基礎(chǔ)上, 采取劃分象限區(qū)的等扇形分析法, 有效地揭示了研究區(qū)近20 年城鎮(zhèn)擴(kuò)張占用耕地的方向變動(dòng), 再從多投入角度構(gòu)建BCC 模型, 計(jì)算4 時(shí)段綜合、規(guī)模、技術(shù)效率, 分析被占用耕地的經(jīng)濟(jì)產(chǎn)出效益, 以期在實(shí)現(xiàn)同城化區(qū)域經(jīng)濟(jì)高質(zhì)量發(fā)展的同時(shí), 能夠?yàn)樘岣叱擎?zhèn)擴(kuò)張占用耕地的效率、減少耕地資源的浪費(fèi)、保障糧食安全提供參考依據(jù)。

1 研究區(qū)概況和研究方法

1.1 研究區(qū)概況

成都都市圈地處四川省東南部（102°54′—105°41′E, 29°24′—31°42′N）, 包括成都市、德陽市、眉山市、資陽市下轄的35 個(gè)縣（市、區(qū)）, 隸屬成渝地區(qū)雙城經(jīng)濟(jì)圈西部重要節(jié)點(diǎn), 位于 “一帶一路” 和西部大開發(fā)戰(zhàn)略主要核心區(qū)域。其地勢(shì)西高東低, 兩山夾平原, 平原廣布。位于亞熱帶濕潤(rùn)季風(fēng)氣候區(qū), 雨熱同期, 夏季炎熱多暴雨。地處長(zhǎng)江上游水系, 河網(wǎng)密布, 以岷江、沱江水系為主。紫色土、水稻土覆蓋較廣, 土質(zhì)較好, 風(fēng)化度低。2019 年底, 地區(qū)生產(chǎn)總值達(dá)22 352.04 億元, 一二三產(chǎn)業(yè)結(jié)構(gòu)比為5.5:33.4:61.1。區(qū)域內(nèi)常住總?cè)丝? 564 萬人, 城鎮(zhèn)化率達(dá)65.5%。以漢族為主, 部分藏族、羌族, 西南客家群體多聚集在洛帶片區(qū), 含道教圣山青城山、都江堰水利工程、寬窄巷子等多個(gè)自然與人文景區(qū), 第三產(chǎn)業(yè)發(fā)達(dá)。

1.2 數(shù)據(jù)來源與處理

各期30 m 分辨率的土地利用數(shù)據(jù)和行政邊界矢量來源于中國科學(xué)院資源與環(huán)境數(shù)據(jù)中心；投入產(chǎn)出指標(biāo)基礎(chǔ)統(tǒng)計(jì)數(shù)據(jù)來源于四川省統(tǒng)計(jì)年鑒及各市統(tǒng)計(jì)年鑒。基于近19 年成都都市圈土地利用現(xiàn)狀, 單獨(dú)提取2000—2005 年、2005—2010 年、2010—2015 年、2015—2018 年的城鎮(zhèn)建設(shè)用地及耕地柵格數(shù)據(jù), 運(yùn)用ArcGIS 疊加得到其城鎮(zhèn)擴(kuò)張占用耕地?cái)?shù)據(jù)并進(jìn)行分區(qū)統(tǒng)計(jì)；根據(jù)統(tǒng)計(jì)年鑒中各縣（市、區(qū)）的社會(huì)固定資產(chǎn)投資、二三產(chǎn)業(yè)從業(yè)人數(shù)、社會(huì)消費(fèi)零售品銷售總額、財(cái)政收入指標(biāo)計(jì)算4 時(shí)段增量, 再與新增建設(shè)用地中耕地所占比例相乘后, 采用極值標(biāo)準(zhǔn)化方法構(gòu)建含有土地投入、資本投入等多角度的效率評(píng)價(jià)數(shù)據(jù)集。

1.3 研究方法

1.3.1 城鎮(zhèn)擴(kuò)張占用耕地測(cè)度 運(yùn)用擴(kuò)張速度（V）、貢獻(xiàn)度指數(shù)（I）[19］, 分析各時(shí)段內(nèi)每年城鎮(zhèn)擴(kuò)張占用耕地的面積, 并確定新增城鎮(zhèn)建設(shè)用地中耕地對(duì)城鎮(zhèn)化貢獻(xiàn)的面積比重, 剖析城鎮(zhèn)擴(kuò)張與耕地轉(zhuǎn)換之間關(guān)系, 計(jì)算公式如下：

式中, Δi為時(shí)間差,Si為某個(gè)時(shí)段新增城鎮(zhèn)擴(kuò)張占用的耕地面積,Ci是某個(gè)時(shí)間段內(nèi)新增城鎮(zhèn)建設(shè)用地面積。

1.3.2 等扇形分析法 以研究區(qū)幾何中心為原點(diǎn), 基于ArcGIS 平臺(tái)構(gòu)建象限坐標(biāo), 將研究區(qū)劃分為16 個(gè)象限區(qū), 分別用區(qū)域所在方向的字母表示。將其與研究時(shí)段內(nèi)城鎮(zhèn)擴(kuò)張信息進(jìn)行聯(lián)立, 以獲取各時(shí)段不同方位城鎮(zhèn)擴(kuò)張占用耕地的方向變動(dòng), 從而進(jìn)一步刻畫被占用耕地的空間格局。

1.3.3 BCC 模型 BCC 模型是學(xué)者在考量實(shí)際經(jīng)濟(jì)運(yùn)行的情況下, 由原來DEA 投入產(chǎn)出模型進(jìn)行改進(jìn), 通過在聚合輸入或輸出時(shí)引入一個(gè)常數(shù), 使生產(chǎn)邊界從原點(diǎn)移開[20, 21］, 更加符合經(jīng)濟(jì)發(fā)展的規(guī)律, 計(jì)算公式如下：

式中,s.t.為約束條件, minθB為目標(biāo)函數(shù),θB為投入產(chǎn)出相對(duì)有效程度,X和Y分別為某個(gè)縣（市、區(qū)）的某項(xiàng)投入指標(biāo)值和產(chǎn)出指標(biāo)值,X0、Y0分別為當(dāng)前決策單元的投入值和產(chǎn)出值。

結(jié)合數(shù)據(jù)可得性和代表性[17］, 完成投入產(chǎn)出指標(biāo)體系構(gòu)建（表1）。

表1 投入產(chǎn)出指標(biāo)體系

2 結(jié)果與分析

2.1 城鎮(zhèn)擴(kuò)張占用耕地?cái)U(kuò)張分析

2.1.1 城鎮(zhèn)擴(kuò)張占用耕地時(shí)間變動(dòng)分析 經(jīng)過Arc-GIS 10.2 軟件將2000、2005、2010、2015、2018 年5 年節(jié)點(diǎn)的耕地和城鎮(zhèn)建設(shè)用地進(jìn)行疊加分析, 獲取4時(shí)段城鎮(zhèn)擴(kuò)張占用耕地的數(shù)據(jù), 經(jīng)過式（1）和式（2）計(jì)算得到各時(shí)段擴(kuò)張速度、貢獻(xiàn)度指數(shù)（表2）。

表2 2000—2018 年城鎮(zhèn)擴(kuò)張占用耕地?cái)U(kuò)張變動(dòng)

2000—2018 年, 成都都市圈城鎮(zhèn)擴(kuò)張占用耕地面積呈先增后減態(tài)勢(shì)。其中, 2000—2005、2005—2010 年2 個(gè)時(shí)段占用耕地面積逐步增加, 之后在2010—2015 年迅速減少, 2015—2018 年有小幅增加；擴(kuò)張速度與占用耕地面積走勢(shì)相同；貢獻(xiàn)度指數(shù)整體呈下降趨勢(shì), 前3 個(gè)時(shí)段占用耕地面積小幅度遞減, 2015—2018 年陡降。該走勢(shì)與近年來中國城鎮(zhèn)化在2015 年后開始內(nèi)涵式發(fā)展, 研究區(qū)域在政策支持下城鎮(zhèn)擴(kuò)張面積雖有增加, 但耕地保護(hù)政策越發(fā)嚴(yán)格密切相關(guān)。

2.1.2 城鎮(zhèn)擴(kuò)張占用耕地空間變化分析 為了增加研究區(qū)域近20 年城鎮(zhèn)擴(kuò)張占用耕地的直觀性, 利用GIS 平臺(tái)進(jìn)行分區(qū)統(tǒng)計(jì), 獲取35 個(gè)縣（市、區(qū)）4 時(shí)段城鎮(zhèn)建設(shè)用地?cái)U(kuò)張占用耕地的空間變動(dòng)（圖2）。

圖2 2000—2018 年4 時(shí)段占用耕地雷達(dá)圖

2000—2018 年, 研究區(qū)35 個(gè)縣（市、區(qū)）中金牛區(qū)、武侯區(qū)、龍泉驛區(qū)、新都區(qū)、溫江區(qū)、雙流區(qū)、郫都區(qū)、旌陽區(qū)等城鎮(zhèn)擴(kuò)張占用耕地范圍較大。武侯區(qū)、雙流區(qū)、郫都區(qū)占用耕地面積超過40 km2, 其次, 彭州市、成華區(qū)、青白江區(qū)、都江堰市、彭州市、崇州市、廣漢市占用耕地面積超過10 km2。2000—2005 年, 武侯區(qū)、郫都區(qū)城鎮(zhèn)擴(kuò)張占用耕地面積最多；2005—2010 年以雙流區(qū)最突出, 而這2 個(gè)時(shí)段成都市的武侯區(qū)、成華區(qū)、龍泉驛區(qū)、新都區(qū)、溫江區(qū)、郫都區(qū)城鎮(zhèn)擴(kuò)張迅速, 占用了較多耕地；2010—2015 年城鎮(zhèn)擴(kuò)張占用耕地面積明顯減少, 各縣（市、區(qū)）新增均不超過10 km2。成都市老主城區(qū)占用耕地面積減小, 而青白江區(qū)、雙流區(qū)在2010—2015 年占用耕地面積最多, 2015—2018 年雙流區(qū)、龍泉驛區(qū)為占用耕地主要區(qū)域, 其他縣（市、區(qū)）占用耕地相比前幾個(gè)時(shí)間段明顯減小。

2.2 城鎮(zhèn)擴(kuò)張占用耕地方向分析

利用GIS 平臺(tái)構(gòu)建211 km 為半徑的緩沖區(qū), 將研究區(qū)域等分為16 部分, 展現(xiàn)各方位城鎮(zhèn)擴(kuò)張占用耕地方向變動(dòng)（圖3）。

圖3 2000—2018 年城鎮(zhèn)擴(kuò)張占用耕地方向變動(dòng)

2000—2018 年研究區(qū)整體城鎮(zhèn)擴(kuò)張占用耕地重心集中在成都市主城區(qū)域內(nèi)部遷移, 其中, 前2 個(gè)時(shí)段由西部往北部移動(dòng), 之后由南部向東部遷徙。2000—2005 年, 西部、西北部的縣（市、區(qū)）城鎮(zhèn)擴(kuò)張占用耕地面積較多。成都市在該時(shí)段強(qiáng)調(diào)發(fā)展西部, 大力完善區(qū)域交通體系, 增加西部交通通達(dá)度, 為西部區(qū)域城鎮(zhèn)化提供基礎(chǔ)條件；2005—2010 年集中在北部、東北部、西北部, 該時(shí)段北部區(qū)域工業(yè)化進(jìn)程加快, 以青白江區(qū)的物流業(yè)和新都區(qū)制造業(yè)快速發(fā)展區(qū)為典型。什邡市、綿竹市由于2008 年汶川大地震后受損嚴(yán)重, 為了盡快恢復(fù)區(qū)域經(jīng)濟(jì)發(fā)展, 同時(shí)考慮區(qū)域地質(zhì)條件特殊性, 在耕地富集平原地區(qū)占用耕地, 保障和改善民生居住環(huán)境, 房地產(chǎn)迎來新一波發(fā)展潮；2010—2015 年開拓方向主要集中南部、東南部, 東北偏東方向也有部分耕地被占用, 成都市為分散中心城市壓力, 不斷往南部雙流區(qū)開拓衛(wèi)星城, 占用了大量耕地。眉山市作為成都市南沿城市, 迎來了新的發(fā)展機(jī)遇, 但發(fā)展還較緩慢, 此時(shí)廣漢市、青白江區(qū)運(yùn)輸業(yè)、技術(shù)園區(qū)逐步引進(jìn), 其占用耕地面積小幅增加；2015—2018 年, 城鎮(zhèn)擴(kuò)張占用耕地力度大的區(qū)域集中在東部、東南部、南部, 其中, 龍泉驛區(qū)、雙流區(qū)成為該時(shí)段加快東進(jìn)、南拓的重要節(jié)點(diǎn)區(qū)域[22］, 多類新城策略推動(dòng)成都東南部城鎮(zhèn)擴(kuò)張。

2.3 城鎮(zhèn)擴(kuò)張占用耕地效率分析

2.3.1 不同效率時(shí)序差異分析 利用DEAP 2.1 軟件計(jì)算得到2000—2018 年成都都市圈4 時(shí)段占用耕地綜合效率、技術(shù)效率及規(guī)模效率。

從圖4 可知, 2000—2018 年成都都市圈城鎮(zhèn)擴(kuò)張占用耕地綜合效率、技術(shù)效率、規(guī)模效率整體均呈先降低后上升趨勢(shì)。就綜合效率而言, 2000—2005年比2010—2015 年效率降低了11.9%, 而后回升至0.950。這是因?yàn)?015 年成渝政策加快實(shí)施后, 簡(jiǎn)陽市、樂至縣等部分較落后縣（市、區(qū)）前期荒置的被占用耕地區(qū)在此時(shí)段開始大量運(yùn)轉(zhuǎn)利用, 在政策吸引下城市周圍村莊勞動(dòng)力回流, 土地綜合效益增加。

圖4 2000—2018 年效率時(shí)序變化

技術(shù)效率在前3 個(gè)時(shí)段內(nèi)呈逐段下降后上升的趨勢(shì), 最低達(dá)到0.857。前2 時(shí)段差距甚微, 可能是由于2000—2010 年國內(nèi)房地產(chǎn)為快速發(fā)展的黃金時(shí)期, 2008 年汶川地震后, 地方獲得大量外部投資幫扶, 同時(shí)加上當(dāng)?shù)孛癖娤M(fèi)觀念的改變, 推動(dòng)了地區(qū)經(jīng)濟(jì)快速恢復(fù), 尤其是綿竹市、什邡市、都江堰市等受災(zāi)嚴(yán)重區(qū)域, 吸引了外資和區(qū)域創(chuàng)新發(fā)展要素、技術(shù)的匯入。2010 年后地產(chǎn)發(fā)展受到控制, 東坡區(qū)、綿竹市、樂至縣等部分縣（市、區(qū)）在土地運(yùn)轉(zhuǎn)管理機(jī)制和技術(shù)創(chuàng)新方面投入減少, 技術(shù)效率得以制約。2015 年受政策驅(qū)使, 區(qū)域迎來新發(fā)展機(jī)遇, 吸引外資、技術(shù)和人才涌入, 為地區(qū)技術(shù)效率提升提供后備力量。

規(guī)模效率在2005—2010 年達(dá)到最低, 前2 時(shí)段相差0.067, 而后逐段上升。2008 年地震造成大量城鎮(zhèn)擴(kuò)張占用耕地區(qū)房屋坍塌, 大片城鎮(zhèn)嚴(yán)重受損, 2010 年后為控制地產(chǎn)投機(jī), 政府對(duì)樓地購置的條件進(jìn)行限制, 同時(shí)企業(yè)為降低成本和減少對(duì)新耕地的渴求, 原已占用耕地荒置區(qū)域在該時(shí)段成為重點(diǎn)發(fā)展對(duì)象, 規(guī)模效率提高。2015 年后, 糧食安全與耕地保護(hù)理念與政策逐步落實(shí)在實(shí)踐上, 成渝城市群的政策出臺(tái)又促進(jìn)了地區(qū)的城鎮(zhèn)化, 二者共同作用下致使規(guī)模效率小幅增加。

2.3.2 綜合效率空間差異分析 為了更直觀地感知其利用效率空間變化, 參考已有研究成果[23, 24］, 將城鎮(zhèn)擴(kuò)張占用耕地的綜合效率值劃分為高效區(qū)（TE=1）、較高效區(qū)（0.8≤TE＜1）、低效區(qū)（0.5≤TE＜0.8）, 對(duì)各時(shí)段35 個(gè)縣（市、區(qū)）綜合利用效率進(jìn)行可視化表達(dá)（圖5）。

圖5 2000—2018 年城鎮(zhèn)擴(kuò)張占用耕地綜合效率空間分異

高效區(qū)整體由中西部往東南遷移, 數(shù)量先增后減, 與城鎮(zhèn)擴(kuò)張占用耕地增量變動(dòng)總體趨勢(shì)一致。2000—2005 年, 除錦江區(qū)、青羊區(qū)、金牛區(qū)、武侯區(qū)外, 彭州市、什邡市、綿竹市等7 個(gè)縣（市、區(qū)）均達(dá)到DEA 有效, 其中, 什邡市、綿竹市憑借相對(duì)較少的土地、勞動(dòng)力投入, 在該時(shí)段生產(chǎn)總值、財(cái)政收入、銷售總額產(chǎn)出較高；2005—2010 年, 主城區(qū)內(nèi)青羊區(qū)效率較低, 成華區(qū)、龍泉驛區(qū)、青白江區(qū)等11 個(gè)縣（市、區(qū)）DEA 的值為1。多為受地震影響較弱區(qū)域, 經(jīng)濟(jì)效益比西部、北部震區(qū)高；2010—2015 年, 雁江區(qū)、中江縣等8 個(gè)縣（市、區(qū)）綜合效率最優(yōu)。成都市龍泉驛區(qū)等投入產(chǎn)出比相當(dāng), 仁壽縣、中江縣能在較低勞動(dòng)力、土地投入下獲取較高的資本產(chǎn)出, 而雁江區(qū)高投入高產(chǎn)出, 均達(dá)到預(yù)期效果；2015—2018 年, 資陽市整體達(dá)到DEA 有效。雁江區(qū)在該時(shí)段土地勞動(dòng)力、資金投入產(chǎn)出都高, 錦江區(qū)等4 個(gè)縣（市、區(qū)）在土地和資本投入較小情況下產(chǎn)出較多, 效率高。

較高效區(qū)呈先減少后增加趨勢(shì), 2005—2010 年和2010—2015 年2 時(shí)段內(nèi)集中在東北到西南方向布局。2000—2005 年, 區(qū)域除郫都區(qū)和原高效區(qū)外, 多在較高效的情形下有高產(chǎn)出, 獲取經(jīng)濟(jì)效益較高；但2005—2015 年受自然災(zāi)害影響, 較高效區(qū)逐段減少。該時(shí)段內(nèi)區(qū)域經(jīng)濟(jì)仍在恢復(fù), 2010 年后只有新津區(qū)、羅江區(qū)等6 個(gè)縣（市、區(qū)）因?yàn)橥度氘a(chǎn)出差距不大效率較高, 加之控制地產(chǎn)熱以及保護(hù)生態(tài)環(huán)境, 土地投入逐步減少, 勞動(dòng)力、資本投入較多, 經(jīng)濟(jì)效益也較高, 但仍有規(guī)模報(bào)酬遞減趨勢(shì)制約其達(dá)到最優(yōu)。2015 年后伴隨地區(qū)經(jīng)濟(jì)快速發(fā)展, 城鎮(zhèn)化轉(zhuǎn)向內(nèi)部填充為主, 地區(qū)被占用耕地經(jīng)濟(jì)效益得到有效提升。

低效區(qū)集中在西部區(qū)域, 數(shù)量變化則與較高效區(qū)相反, 先增后減。其中, 2005—2010 年和2010—2015 年2 時(shí)段內(nèi)數(shù)量較多, 是因?yàn)?008 年地震造成西北部綿竹市、什邡市、都江堰市等經(jīng)濟(jì)損失嚴(yán)重, 同時(shí)地方投入大量資金、人力推動(dòng)經(jīng)濟(jì)復(fù)蘇, 但恢復(fù)過程較慢。南部青神縣、蒲江縣、安岳縣等地政府投入力度雖大, 但距離主城較遠(yuǎn), 產(chǎn)出較少；2000—2005 年郫都區(qū)投入較高, 但當(dāng)時(shí)其距離主城較遠(yuǎn), 地理位置優(yōu)越性還不突出, 資本產(chǎn)出相對(duì)較少；2015年后, 成渝城市群聯(lián)動(dòng)發(fā)展, 經(jīng)濟(jì)發(fā)展迎來新機(jī)遇, 區(qū)域效益整體提升, 但洪雅縣由于位置偏且地形險(xiǎn), 近年來瓦屋山等景區(qū)開發(fā)助推該地區(qū)經(jīng)濟(jì)產(chǎn)業(yè)結(jié)構(gòu)更新, 但整體效率仍然較低。

3 小結(jié)與討論

3.1 小結(jié)

本研究基于GIS 平臺(tái)獲取成都都市圈2000—2018 年4 時(shí)段城鎮(zhèn)擴(kuò)張占用耕地?cái)?shù)據(jù), 在此基礎(chǔ)上分析區(qū)域強(qiáng)度和貢獻(xiàn)度指數(shù)、方向及利用效率。

1）研究區(qū)域城鎮(zhèn)擴(kuò)張占用耕地整體強(qiáng)度呈現(xiàn)先增強(qiáng)后減弱的趨勢(shì)。整體占用耕地面積470.528 km2, 4 時(shí)段在2005—2010 年擴(kuò)張速度達(dá)到最大, 為41.383 km2/年, 2010—2015年減至最小, 為4.198 km2/年。貢獻(xiàn)度指數(shù)在2000—2005 年呈下降趨勢(shì)。

2）研究區(qū)域城鎮(zhèn)擴(kuò)張占用耕地重心整體由西部、北部往東部、南部遷移, 后期主要集中在龍泉驛區(qū)、雙流區(qū)等新興縣（市、區(qū)）。

3）在研究區(qū)域效率時(shí)序變化上, 綜合效率、技術(shù)效率、規(guī)模效率整體呈現(xiàn)先下降后上升趨勢(shì)。在空間變化上, 綜合效率高效區(qū)由中西部往東南部集中, 數(shù)量先增后減；而較高效區(qū)在首尾兩時(shí)段覆蓋廣, 中間兩時(shí)段主要在東北到西南方向分布；低效區(qū)多位于西部, 數(shù)量先增后減。

3.2 討論

2000—2018 年成都都市圈城鎮(zhèn)擴(kuò)張占用耕地面積呈現(xiàn)先增長(zhǎng)后減少趨勢(shì), 符合中國2000 年后城鎮(zhèn)化加速和2015 年轉(zhuǎn)向內(nèi)涵式發(fā)展的城鎮(zhèn)化規(guī)律。整體向東部和南部方向演進(jìn), 但其城鎮(zhèn)擴(kuò)張占用耕地高效區(qū)減少, 與已有研究結(jié)果較為一致[17］。主要存在某些區(qū)域由于自身經(jīng)濟(jì)實(shí)力弱、資本和人才投入還較少、部分區(qū)域存在占用地閑置乃至荒廢等問題。為有效提高其經(jīng)濟(jì)效益, 高效區(qū)應(yīng)當(dāng)繼續(xù)保持當(dāng)下資本和人力投入, 同時(shí)利用完善交通體系, 為周圍地區(qū)發(fā)展鋪設(shè)基礎(chǔ)條件[25］；較高效區(qū)要抓住發(fā)展機(jī)遇, 利用資本和政策進(jìn)行技術(shù)更新。同時(shí)保持與中心城市的耕地占補(bǔ)平衡機(jī)制同步, 實(shí)現(xiàn)二者的有效聯(lián)動(dòng)；低效區(qū)應(yīng)該充分挖掘自身產(chǎn)業(yè)特色, 延長(zhǎng)產(chǎn)業(yè)結(jié)構(gòu)鏈, 保障被占用耕地能夠有充足后備條件支持其發(fā)展。

同時(shí)本研究也有不足之處, 在產(chǎn)出指標(biāo)選取時(shí), 由于2010 年后個(gè)別縣（市、區(qū)）農(nóng)業(yè)機(jī)械投入數(shù)據(jù)獲取困難, 城鎮(zhèn)擴(kuò)張占用耕地的技術(shù)投入無法有效顯現(xiàn), 在產(chǎn)出上未能考慮生態(tài)環(huán)境指標(biāo), 具有一定的局限性。因而, 下一步應(yīng)該通過其他途徑或者尋找有效替代數(shù)據(jù), 分析技術(shù)投入對(duì)于當(dāng)下占用耕地經(jīng)濟(jì)效益的作用。為了追尋經(jīng)濟(jì)高質(zhì)量的發(fā)展和優(yōu)先保護(hù)生態(tài)環(huán)境, 計(jì)量占用耕地產(chǎn)生的經(jīng)濟(jì)效益和生態(tài)效益并進(jìn)行二者之間的耦合關(guān)聯(lián)將是主要改進(jìn)方向。