玉馬經(jīng)濟帶農(nóng)用地耕地果園轉換時空特征研究

楊雪萃　張勇榮　徐異龍　童賢倩　李萬錄

摘 ?要：農(nóng)用地中“耕地-果園”轉化作為鄉(xiāng)村振興背景下重要的土地利用轉化方式，對促進區(qū)域高質量可持續(xù)發(fā)展有深刻影響，為揭示“玉馬經(jīng)濟帶”耕地-果園轉換時空變化特征及影響因素，采用2012年與2018年遙感影像和土地利用數(shù)據(jù)，通過轉移矩陣、核密度分析等方法，旨在揭示2012—2018年期間土地利用格局時空演變特征和影響因素。結果表明：（1）2012—2018年地類轉化趨勢明顯，耕地面積顯著減少，果園面積逐年增加且總體呈穩(wěn)定增長態(tài)勢，總體形成集中連片分布模式。（2）“耕地-果園”轉換熱點區(qū)域主要表現(xiàn)為三類，分別是一號地塊的獨立擴展模式、二號地塊的帶狀延伸模式和三號地塊的團狀發(fā)展模式。（3）“耕地-果園”轉換，其影響因素主要受制于宏觀政策、坡度、交通、等綜合因素驅動?！坝?馬經(jīng)濟帶”土地利用結構未來將趨向規(guī)?；?、均衡化發(fā)展，將形成科學高效的土地利用結構。

關鍵詞：玉-馬經(jīng)濟帶 ?耕地-果園轉換 ?遙感解譯 ?時空特征 ?核密度分析

中圖分類號：F323.211??????? ?文獻標識碼：A? ????????

Research on the Spatial-Temporal Characteristics of Farmland-Orchard Transformation in Agricultural Land in the Yushe-Machang Economic Belt

YANG Xuecui? ZHANG Yongrong XU Yilong? TONG Xianqian? LI Wanlu

（Liupanshui Normal University， Liupanshui， Guizhou Province， 553000 China）

Abstract： In agricultural land， as an important land use transformation method in the context of rural revitalization， "farmland-orchard" transformation has a profound impact on promoting the high-quality and sustainable development of the region. In order to reveal the spatial-temporal variation characteristics and influence factors of farmland-orchard transformation in the "Yushe-Machang Economic Belt"， this paper uses remote sensing images and land use data in 2012 and 2018 to reveal the spatial-temporal evolution characteristics and influence factors of land use patterns from 2012 to 2018 through methods such as transition matrix and kernel density analysis.?The results are as follows. （1） From 2012 to 2018， the trend of land type transformation was obvious， cultivated land area decreased significantly， and orchard area increased year by year and showed an overall stable growth trend， which generally formed a centralized contiguous distribution mode. （2） There were three main types of the hot spot regions of "farmland-orchard" conversion， namely， the independent expansion mode presented by Plot 1， the banded extension mode presented by Plot 2， and the cluster development mode presented by Plot 3. （3） The influence factors of "farmland-orchard" conversion were mainly driven by comprehensive factors such as macro policies， slope and traffic. The land use structure in the "Yushe-Machang Economic Belt" will trend to scale and balanced development in the future， and a scientific and efficient land use structure will be formed.

Key Words： Yushe-Machang Economic Belt; Farmland-orchard transformation; Remote sensing interpretation; Spatial-temporal characteristics; Kernel density analysis

全球環(huán)境變化與可持續(xù)發(fā)展越來越注重人類對土地資源的開發(fā)、利用及土地利用類型的轉化，其中農(nóng)業(yè)用地的發(fā)展對生態(tài)環(huán)境與社會經(jīng)濟的影響極為顯著，由此農(nóng)用地變化成了土地利用/土地覆被變化的重要研究內容。在當今快速城市化進程發(fā)展中和鄉(xiāng)村農(nóng)業(yè)振興發(fā)展的背景下，農(nóng)用地類轉化顯著，大量耕地轉化為其他經(jīng)濟類型農(nóng)用地類。目前，對土地利用轉型研究主要在地類轉型動因、土地利用轉型影響因素，還有一些關于特殊地類利用轉型以及轉型研究的探測、討論^[1]。土地利用研究有多種維度、多種方法，系統(tǒng)推進其理論與實際研究的升華，提升研究方法的準確科學性、實地適應性對推進土地利用研究有深遠意義。

土地利用類型變化是一個極具國際性的前沿課題。國外現(xiàn)有研究多集中于農(nóng)業(yè)用地的功能、轉型、未來演變趨向等方面^[2]，以及林地、建設用地、住宅用地等轉型研究。在土地利用類型中，農(nóng)用地土地資源利用方式、變化格局深刻影響到區(qū)域生產(chǎn)結構穩(wěn)定性、合理性問題。耕地作為生產(chǎn)用地的主要類型和基本資源的重要來源，對其優(yōu)化配置可有效提升耕地資源優(yōu)勢，促進區(qū)域生產(chǎn)可持續(xù)發(fā)展?，F(xiàn)如今，大量新技術應用到土地利用變化研究，尤其是以遙感解譯、地理信息系統(tǒng)為主的技術應用到地類變化特征、驅動因素和直接影響等方面的研究。

1 研究區(qū)概要、數(shù)據(jù)來源和研究方法

1.1? 研究區(qū)概要

研究區(qū)位于貴州省六盤水市水城區(qū)，包含玉舍鎮(zhèn)、勺米鎮(zhèn)、米蘿鎮(zhèn)和阿戛鎮(zhèn)等，如圖1所示。經(jīng)緯度位置為104?34′30″～104?57′30″E，26?24′15″～25?31′43″N，研究區(qū)總面積約139 km²，境內山脈延伸呈西北東南走向，地勢西高東低中部低洼，平均海拔約為1 623 m，是典型的北亞熱帶高原喀斯特巖溶景觀區(qū)。

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

文中采用的數(shù)據(jù)源于2012年和2018年的Google?Earth高分辨率遙感影像數(shù)據(jù)，以及空間分辨率為30m的DEM數(shù)據(jù)，坡度根據(jù)其計算得出。使用ENVI和ArcGIS軟件并結合實際情況進行面向對象分類，加以人工修改有誤圖斑得出2012年和2018年的土地利用類型圖，如圖2所示作為基礎數(shù)據(jù)來源。

1.3 ?研究方法

1.3.1 土地實際利用功能分類

根據(jù)2017年11月1日出版的中國土地利用現(xiàn)狀分類（Current land use classification）^[3]并結合研究區(qū)實際概況，將其土地利用類型分為12類，具體如表1所示。

1.3.2 土地利用轉移矩陣分析

利用轉移矩陣，可以系統(tǒng)、直觀地表明不同時期土地利用類型之間的轉變情況^[4]。通過2012年與2018年土地利用矢量數(shù)據(jù)相交并疊加分析，得出反映2012年與2018年研究區(qū)土地利用類型轉化的面積變化表，如圖3所示。

1.3.3 土地利用/覆蓋變化率指數(shù)

為了反映各類型土地資源面積的變化幅度與速度以及研究區(qū)土地利用或覆蓋變化中的類型差異，使用土地利用/覆蓋類型的變化指數(shù)，通過對土地利用程度的分級，量化人類活動對于土地系統(tǒng)的影響程度，定量地描述該地區(qū)土地利用的綜合水平和變化趨勢，是衡量指定區(qū)域內土地利用情況的一個重要指標^[5]。從不同方面反映變化率指數(shù)，主要通過時段變化率和年平均變化率進行分析，表達式分別為：

P = （U_b-U_a）/U_a×100%

R_s=[（U_b-U_a）/U_aT]×100%

式中，Ｐ為某種類型的時段變化率；Ｒ_ｓ為年均變化率；Ｕ_ａ、Ｕ_ｂ分別為研究初期、末期某一類型的面積；Ｔ為研究時段，單位設為年^[6]。

1.3.4 核密度分析

核密度分析是使用核函數(shù)根據(jù)點要素計算每單位面積的量值以將各個折點擬合為光滑錐狀表面，用于計算點要素測量值在指定鄰域范圍內的單位密度。其熱力圖像表達能直觀地反映出離散測量值在連續(xù)區(qū)域內的分布情況，也能更直觀地對比地理事物分布特征與關系^[7]，能較好地反映耕-果轉換的熱點區(qū)域，主要利用spatial analyst工具中的密度分析中的核密度分析進行計算。

先將研究區(qū)2012年與2018年土地利用矢量數(shù)據(jù)進行相交，再將相交的矢量數(shù)據(jù)進行轉化，通過數(shù)據(jù)管理工具中的要素轉點轉化為興趣點即POI^[8]，研究區(qū)POI分布較為密集且分布集中；將興趣點放入ArcMap中的核密度分析，核密度變化特征數(shù)據(jù)，如圖4所示。

2 ?結果分析

2.1? 時間變化特征分析

研究區(qū)土地利用格局的數(shù)量變化特征鮮明^[9]。通過圖5以及土地利用/覆蓋變化率指數(shù)、土地利用轉矩表可知2012—2018年間農(nóng)用地類“耕地-果園”變化明顯，呈上升趨勢。其他地類也有明顯變化，其他用地的比例在不斷增加，且變化趨勢明顯。

草地面積所占比例不斷減少；灌木林地變化比例較小，但整體呈增加趨勢；林地面積由所占比例有所增加，幾種地類轉化說明研究區(qū)進行了大面積的退耕還林和生態(tài)修復，該數(shù)據(jù)符合研究區(qū)域作為重點生態(tài)功能區(qū)的特征，這對于本研究區(qū)的生態(tài)環(huán)境和土地資源綜合利用建設有重要意義。工礦與倉儲用地、住宅用地、交通用地和公共管理與服務用地面積變化較小，總體呈小幅度上升趨勢，水域設施用地面積無變化較為穩(wěn)定，這說明研究區(qū)內城鎮(zhèn)化進程和基礎設施建設在逐步進行中，未來變化趨勢將更加明顯。

2.2? 耕地-果園轉換空間特征分析

研究區(qū)呈西北-東南走向的長條狀，其耕地主要分布于中部河谷、低地地區(qū)，還有少量分布在海拔較高的坡耕地，導致耕地-果園轉化的空為中部地區(qū)。為了更直觀地展現(xiàn)研究區(qū)“耕地-果園”轉換熱點區(qū)域，將研究區(qū)核密度分析與土地利用空間結構耦合分析后，將其分類為3種模式：獨立擴展模式、帶狀延伸模式、團狀發(fā)展模式。具體分析如下。

2.2.1 獨立擴展模式

圖中一號地塊，是玉舍鎮(zhèn)中心城區(qū)，主要為工礦倉儲用地、建筑用地、住宅用地等生活或生產(chǎn)用地。此區(qū)域發(fā)生地類轉變的熱點區(qū)域主要是中部及北部，其中北部地類熱點變化是以建筑用地為中心的獨立輻射擴展，而中部的熱點變化是以工礦倉儲用地、公共管理與服務用地為中心的獨立輻射擴展；一號地塊中農(nóng)用地轉換占比較小，因為農(nóng)用地面積較少，城鎮(zhèn)化水平在研究區(qū)內屬較高水平。在2012—2018年期間此地塊地類變化持續(xù)且變化顯著。

2.2.2 帶狀延伸模式

圖中二號地塊，主要是以農(nóng)用地、林地、草地為主的鄉(xiāng)村區(qū)域，此區(qū)域發(fā)生熱點地類轉換的主要是中部河谷低地及東南部，呈連續(xù)帶狀延伸，其中部河谷周邊區(qū)域耕地轉換為果園的面積最大，“耕地-果園”轉換最突出，是研究區(qū)中農(nóng)用地類轉換的重點區(qū)域，此地區(qū)地勢較為平坦，中部有河流流經(jīng)提供灌溉水域，有主要交通干線經(jīng)過，利于農(nóng)產(chǎn)品的運輸，總體來說區(qū)域優(yōu)勢明顯有利于耕地-果園轉化；二號地塊中其他地類轉換占比較小，主要為一些生態(tài)用地的轉換。2012—2018年間總體為先呈現(xiàn)退耕還林、坡改梯的土地利用變化格局，后呈現(xiàn)耕地轉為果園的變化格局，整體變化顯著且持續(xù)。

2.2.3 團狀發(fā)展模式

圖中三號地塊，此區(qū)域發(fā)生地類轉變的熱點區(qū)域主要是其南部，地勢低平，且有河流和交通線經(jīng)過，以建筑用地和農(nóng)用地為主。其中以耕地-果園轉化最為突出，但其他地類轉換占比較小。此外，還有大量設施農(nóng)用地轉化，新建有大規(guī)模塑料大棚、農(nóng)產(chǎn)品加工處理設施等。2012—2018年期間耕地轉換為果園、設施農(nóng)用地的變化較為頻繁且集中。

3種模式表現(xiàn)了其轉化的空間特征與空間格局，其中三號地塊南部，是整個研究區(qū)中耕地-果園轉換面積最大的地區(qū)；二號地塊耕地-果園轉化也較為突出，但其他地類轉化活躍；一號地塊主要為生活用地的轉化，即建筑用地、公共管理與服務用地、住宅用地等。

2.3? 農(nóng)用地時空變化特征影響因素分析

（1）了解“玉-馬經(jīng)濟帶”農(nóng)用耕地轉化為果園的影響因素，可為后期耕地轉換利用提供指導，也有利于促進當?shù)亟?jīng)濟高質量和可持續(xù)發(fā)展。在研究區(qū)生態(tài)較脆弱的喀斯特地區(qū)中，土地利用方式的空間布局易受自然因素如土壤肥力、水源、地形等因素的影響 ，其中坡度影響最為明顯，如圖7所示，將研究區(qū)地形坡度分為5類，再與土地利用數(shù)據(jù)進行疊合分析后得出：坡度大于25°及以上的山地地區(qū)“耕地-果園”轉換比例較小，坡度小于5°的河谷低地或丘陵、平地“耕地-果園”轉換比例大，且其他地類以及總轉化面積都處于較高數(shù)值。

（2）土地利用方式空間布局除了受自然因素影響之外還易受交通路線分布、地區(qū)開發(fā)狀況、基礎設施建設狀況、經(jīng)濟發(fā)展狀況和居民點分布等社會因素影響，其中交通路線分布對其影響具體表現(xiàn)為交通通達度對于農(nóng)業(yè)日常耕作、經(jīng)營管理、集約化生產(chǎn)和產(chǎn)品運輸有重要影響（見圖8）。

研究區(qū)中“玉馬線”從中間經(jīng)過并與其他交通線共同形成密集的交通網(wǎng)絡，基于交通線路在ArcMap中進行緩沖區(qū)分析，設置四級緩沖區(qū)單位為100 m，結果如圖9所示，將交通緩沖區(qū)與地類數(shù)據(jù)疊加分析后得出，緩沖100 m的地區(qū)耕地-果園轉換面積最大，緩沖距離增加而轉換面積減小，由此說明交通線分布以及距離交通線的距離對其轉化有顯著影響。

（3）2012年之后，受國家生態(tài)保護政策和居民生態(tài)環(huán)保意識的影響，研究區(qū)大量山地退耕還林、坡改梯治理、水土保持等生態(tài)環(huán)保工程有條不紊地進行，顯著表現(xiàn)為草地面積減小，灌木林地、林地面積增加，這在一定程度上影響了“玉-馬經(jīng)濟帶”的土地利用結構變化和生態(tài)可持續(xù)發(fā)展，提升了植被覆蓋率和生態(tài)指數(shù)。

3 討論

在研究區(qū)農(nóng)用地類轉換過程中，果園和設施農(nóng)用地面積逐漸增加，其轉換來源主要是耕地。而由于研究區(qū)地形、坡度的限制，河谷低地、平地面積有限，在未來的果園產(chǎn)業(yè)擴展過程中勢必面臨可轉換土地資源短缺的狀況。由于果園是“玉-馬經(jīng)濟帶”土地轉換利用的主流地類，與之對應的是未來耕地、其他地類面積的減少和建設用地、交通用地和公共管理與服務用地等地類面積的增加，這在一定程度上會對當?shù)赝恋刭Y源的開發(fā)和合理配置產(chǎn)生一定影響^[10]，改變土地利用格局和效能。

農(nóng)用地利用方式空間布局的變化，是由多種影響因素綜合影響而形成的。在“玉-馬經(jīng)濟帶”中耕地-果園轉換主要受到政策、地形和交通線等影響，深刻且顯著地改變了土地資源的利用方式。地形作為喀斯特山區(qū)中自然因素中最重要的因素，其影響巨大。其中坡度可影響土壤的水分分布、養(yǎng)分分布以及作物運輸便捷性等，例如：地形因子地勢較高處土壤有機碳含量高，地勢低洼處含量低，不同坡位土壤有機碳含量上坡>坡頂>中坡>下坡>坡腳>洼地^[11]。此外，交通便捷程度是社會經(jīng)濟因素中相當重要的因素，交通通達度以及距離主要交通干線的距離影響了耕地-果園的轉換。首先果品對運輸時效性有較高要求，其次規(guī)模化經(jīng)營的果園需要大型機械進入采摘運輸?shù)?。在貴州實施“縣縣通高速”的背景下，研究區(qū)也在大力建設小規(guī)模的交通線路方便進入管理、采摘等生產(chǎn)活動。這些交通網(wǎng)絡的建設減少了交通運輸對農(nóng)業(yè)發(fā)展長期的“瓶頸”制約，有效地促進農(nóng)業(yè)高質量可持續(xù)發(fā)展^[12]。

基于此，怎樣推進“玉-馬經(jīng)濟帶”土地資源的合理和高效利用，實現(xiàn)綜合效益共贏，加強研究區(qū)土地資源科學規(guī)劃和可持續(xù)利用將是未來“玉-馬經(jīng)濟帶”發(fā)展規(guī)劃中仍需解決的關鍵問題，因此，未來的研究可引入相關理論和方法，數(shù)字地球、地理信息系統(tǒng)與農(nóng)用地研究相結合^[13]。以促進形成高效的土地利用模式和促進當?shù)剞r(nóng)業(yè)的科學發(fā)展與繁榮興盛，為日后的鄉(xiāng)村振興提供理論支持和指導。

4? 結論

本文通過對“玉-馬經(jīng)濟帶”2012年、2018年高分辨率遙感影像進行遙感解譯、人工修改有誤圖斑和數(shù)據(jù)整理，得到兩期土地利用類型數(shù)據(jù)，并利用ArcMap分析工具和一系列研究方法整理得出2012—2018年6年間農(nóng)用地耕地轉換為果園的時間空間變化特征和影響因素，結論如下。

（1）農(nóng)用地類在六盤水“三變政策”和鄉(xiāng)村振興政策等影響下“耕地-果園”轉化趨勢明顯，耕地面積顯著減少，果園面積逐年增加且呈穩(wěn)定增長態(tài)勢，總體形成集中連片的分布模式。耕地在六年間由45.41 km²減少到27.23 km²，所占比例也由32%減少到20%；果園則由6.08 km²增加到15.34 km²，所占比例也由4%增加到22%；設施農(nóng)用地所占比例則由1%增加到9%，其轉化揭示了研究區(qū)內農(nóng)業(yè)用地的變化速率與變化格局。

（2）“耕地-果園”轉換熱點區(qū)域在中部河谷低地最為明顯，主要表現(xiàn)為二號地塊中部河谷低地及東南部呈現(xiàn)帶狀延伸發(fā)展模式，三號地塊南部呈團狀集中發(fā)展模式。其他地類的轉換熱點區(qū)域在一號地塊最為明顯，主要集中于中部及北部，呈現(xiàn)獨立擴展模式。3種變化模式深刻影響了研究區(qū)土地資源空間配置。

（3）“耕地-果園”轉換不斷發(fā)展演化，主要受制于宏觀政策影響，以及自然因素驅動，主要表現(xiàn)為受地形坡度的影響較大；還受到社會經(jīng)濟因素影響，顯著表現(xiàn)為交通路線和距離交通遠近的影響等。 “玉-馬經(jīng)濟帶”農(nóng)用地轉化未來演變將趨向規(guī)?；?、均衡化發(fā)展，形成科學高效的土地利用格局。

參考文獻

[1] 龍花樓. 論土地利用轉型與土地資源管理[J]. 地理研究，2015，34（9）：1607-1618.

[2] 梁鑫源，李陽兵. 三峽庫區(qū)“耕—果”轉換時空變化特征及其啟示：以草堂溪流域為例[J]. 自然資源學報，2019，34（2）：385-399.

[3] 李森盛. 基于遙感影像處理技術的土地利用變化研究[J]. 測繪與空間地理信息，2023，46（1）：94-96，99.

[4] 鄧佩如，羅力兵，田丹，等. 基于轉移矩陣的六盤水市鐘山區(qū)土地利用研究[J]. 綠色科技， 2021， 23（8）：194-197.

[5] 王子龍，朱彤，姜秋香. 黑龍江省土地利用/覆被變化鐘山區(qū)土地利用研究[J]. 綠色科技， 2021， 23（8）：194-197.

[5] 王子龍，朱彤，姜秋香. 黑龍江省土地利用/覆被變化及生境質量特征變化分析[J]. 東北農(nóng)業(yè)大學學報，2022，53（1）：77-86.

[6] 祖立輝，楊靜，苗世源. 基于遙感影像的臨湘市土地利用／覆蓋時空變化研究[J]. 測繪地理信息，2021，46（2）：30-35.

[7] 王雯雯，谷優(yōu)陽，張馳，等. 基于POI 及核密度分析的咸豐縣城鄉(xiāng)交錯帶分布研究[J]. 資源環(huán)境與工程，2022，36（1）：117-123.

[8] 胡良柏. 基于核密度分析與關聯(lián)規(guī)則的天水盆地黃土滑坡規(guī)律研究[J]. 測繪與空間地理信息，2022，45（1）：30-32，9.

[9] 左岍，周勇，李晴，等. 鄂西南地區(qū)土地利用格局時空變化及軌跡特征分析[J]. 水土保持學報，2022，36（1）：161-169.

[10] 何春陽，張金茜，劉志鋒，等.1990—2018 年土地利用/覆蓋變化研究的特征和進展[J]. 地理學報，2021，76（11）：2730-2748.

[11] 劉紀遠，劉明亮，莊大方，等. 中國近期土地利用變化的空間格局分析[J]. 中國科學（D 輯：地球科學），2002（12）：1031-1040，1058-1060.

[12] 施園園. 大都市區(qū)耕地多功能評價[D]. 北京：中國地質大學，2016

[13] 李會. 喀斯特小流域土壤有機碳空間分布及影響因素研究[D]. 貴陽：貴州大學，2016.