耕地非農(nóng)化遷移路徑研究

孫佩 尹偉 易潔偉 康全國 張馳 袁知洋

摘要 ［目的］通過對2009—2021年湖北省丹江口市耕地遷移路徑的研究，為有效掌握耕地變化情況及精準(zhǔn)管理提供參考。［方法］運用Arcgis空間疊加分析及統(tǒng)計工具，并借鑒利菲弗方向性分布，分析丹江口市各鄉(xiāng)鎮(zhèn)耕地非農(nóng)化過程及空間分布特征。［結(jié)果］耕地非農(nóng)化方向呈現(xiàn)中部—東南—西北—西北的遷移路徑，方位角有一定的擺動性，呈現(xiàn)順時針—逆時針—逆時針特征，遷移距離呈逐年減小趨勢。丹江口市各鄉(xiāng)鎮(zhèn)耕地非農(nóng)化空間存在一定差異性。整體而言，高值區(qū)集中分布于東部及西部地區(qū)。［結(jié)論］揭示了丹江口市耕地非農(nóng)化過程及遷移路徑，為政府部門把控耕地非農(nóng)化機制及完善國土空間結(jié)構(gòu)提供科學(xué)依據(jù)。

關(guān)鍵詞 丹江口市；耕地非農(nóng)化；遷移路徑；利菲弗方向性分布

中圖分類號 F 323.211? 文獻(xiàn)標(biāo)識碼 A? 文章編號 0517-6611（2023）23-0074-04

doi：10.3969/j.issn.0517-6611.2023.23.017

Study on the Migration Path of Cultivated Land Non Agriculturalization—A Case Study of Danjiangkou City

SUN Pei， YIN Wei， YI Jie-wei et al

（Hubei Institute of Geo-sciences，Wuhan， Hubei 430000）

Abstract ［Objective］Through the research on the migration path of cultivated land in Danjiangkou of Hubei Province from 2009 to 2021， it provides reference for effectively grasping the change of cultivated land and accurate management. ［Method］Using Arcgis spatial overlay analysis and statistical tools， and referring to the Lefevere directional distribution， this paper analyzes the process and spatial distribution characteristics of farmland conversion in various villages and towns of Danjiangkou City. ［Result］The results showed that the direction of cultivated land conversion showed a migration path of central-southeast-northwest-northwest， the azimuth angle had a certain swing， showing the characteristics of clockwise-counterclockwise-counterclockwise， and the migration distance decreased year by year. There are some differences in the non-agricultural space of cultivated land in different villages and towns of Danjiangkou City. On the whole， the high value areas are concentrated in the east and west. ［Conclusion］It reveals the process and migration path of farmland conversion in Danjiangkou City， providing scientific basis for government departments to control farmland conversion mechanism and improve land spatial structure.

Key words Danjiangkou City;Farmland conversion;Migration path;Lefevere directivity distribution

基金項目 荊門市耕地質(zhì)量生態(tài)監(jiān)測試點（續(xù)作）（MSDZ202212）； 編制全省耕地保護(hù)規(guī)劃項目（HBECC-ZB-ZC22097）。

作者簡介 孫佩（1990—），女，湖北荊州人，助理工程師，碩士，從事土地生態(tài)與資源以及土地遙感和農(nóng)情監(jiān)測研究。

收稿日期 2022-12-11

我國一直很重視“三農(nóng)”問題，近年來的變更調(diào)查工作已將耕地列為重點考核對象，對于即將開展的新一輪變更調(diào)查工作，《自然資源部辦公廳關(guān)于開展2022年度全國國土變更調(diào)查工作的通知》［1］中明確提出“各級耕地保護(hù)監(jiān)督部門要對耕地變化情況及原因進(jìn)行分析，認(rèn)真排查整改耕地違法違規(guī)轉(zhuǎn)為其他農(nóng)用地的情況，及時將整改結(jié)果反饋調(diào)查部門”，由此可見，耕地問題是土地資源調(diào)查對象的重中之重。

現(xiàn)階段全國各地從不同角度對耕地展開了研究。錢鳳魁等［2］以遼寧省為研究對象，探討了不同地貌區(qū)耕地景觀格局與耕地質(zhì)量之間的關(guān)聯(lián)性及作用機制，揭示了二者之間的正向關(guān)聯(lián)性特征。對于耕地信息的提取，范影影等［3］以大同市云州區(qū)4鄉(xiāng)鎮(zhèn)為例，利用面向?qū)ο蠓指畹姆椒ㄌ崛×烁丶羞B片經(jīng)營分布的區(qū)域。在耕地時空變化方面，眾多學(xué)者從變化特征［4-6］、過程［7］、影響因子［8-9］、驅(qū)動力［5，10］等方面進(jìn)行了研究，剖析了不同區(qū)域耕地在時間與空間上的變化情況及成因。對于耕地的集約利用，杜國明等［11］解析了耕地集約利用轉(zhuǎn)型的理論，胡賢輝等［12］則通過分析耕地集約利用與生態(tài)效率兩者耦合協(xié)調(diào)機理，探究了耦合協(xié)調(diào)度與空間自相關(guān)2種模型的耦合協(xié)調(diào)現(xiàn)狀及時空分異特征，提升了耕地集約利用與生態(tài)效率的耦合協(xié)調(diào)水平。對于耕地發(fā)展的對策研究，郝軍等［13］立足于生態(tài)文明的背景探究了內(nèi)蒙古耕地的發(fā)展對策，王慧等［14］以江蘇省為研究區(qū)，基于耕地多功能的角度，探討了耕地保護(hù)補償標(biāo)準(zhǔn)，梁鑫源等著［15］眼于多情景糧食安全底線約束科學(xué)解析了我國耕地保護(hù)的彈性空間，為我國耕地保護(hù)制度轉(zhuǎn)型等方面提供了有效的借鑒價值。在耕地質(zhì)量評價方法上，何新瑩等［16］利用改進(jìn)灰靶模型，提出了縣域耕地質(zhì)量評價新方法，為新時代耕地多功能質(zhì)量評價提供理論和方法借鑒，王來剛等［17］以河南省安陽縣為研究對象，基于遙感數(shù)據(jù)對縣域耕地質(zhì)量進(jìn)行了評價與分析。近年來也出現(xiàn)了諸多關(guān)于耕地資源價值的核算［18-19］及耕地低碳利用績效測算［20］等方面的研究方法。

綜上所述，現(xiàn)有的研究絕大部分聚焦于耕地信息的提取、耕地時空動態(tài)變化特征、驅(qū)動因子、耕地質(zhì)量與保護(hù)、耕地評價等方面，對于縣域尺度耕地的流出及遷移路徑的研究較少，從較長時間維度分析耕地的轉(zhuǎn)化遷移特征，對于研究區(qū)耕地的有效保護(hù)、可持續(xù)發(fā)展以及國土整治等方面有重大意義。

1 研究區(qū)概況

丹江口市位于湖北省西北部偏東，110°47′53″～110°34′47″E，32°14′10″～32°58′10″N，區(qū)內(nèi)以山地為主，平均海拔約400 m，地勢中間低，南北高。全市總面積約3 121 km2，2021年耕地面積約252.48 km2，占全市總面積的8.09%。以漢江為界，市域劃分為江南與江北兩大自然區(qū)，氣候以亞熱帶半濕潤季風(fēng)氣候區(qū)為主，年平均氣溫15.9 ℃，轄20個鎮(zhèn)（辦、處、區(qū)）194個村，總?cè)丝诩s47萬，是南水北調(diào)中線工程核心水源區(qū)［21］。

2 數(shù)據(jù)與方法

2.1 數(shù)據(jù)來源及預(yù)處理

采用的矢量數(shù)據(jù)來源于自然資源規(guī)劃局提供的5期國土變更調(diào)查數(shù)據(jù)庫（2009、2012、2015、2018、2021年）。由于收集到的數(shù)據(jù)是全域性的土地利用情況，地類繁多，而主要研究不同時間階段的耕地非農(nóng)化的特征，因此對收集到的原始數(shù)據(jù)需先進(jìn)行重分類處理。依據(jù)非農(nóng)化定義［22］以及《土地利用現(xiàn)狀分類》國家標(biāo)準(zhǔn)（GB/T 21010—2017）［23］，將每期數(shù)據(jù)劃分為耕地、建設(shè)用地與其他三大地類，并在Arcgis軟件中對不同時期非耕地耕地面積進(jìn)行統(tǒng)計輸出。

2.2 研究方法

利菲弗方向性分布工具在空間統(tǒng)計中是綜合能力最突出的工具之一，有著廣泛的應(yīng)用，該算法在1926年最早由美國南加州大學(xué)社會學(xué)教授韋爾蒂·利菲弗（D.Welty Lefever）提出［24］， 可以用來度量概括地理要素中心趨勢、離散和方向趨勢的空間特征［25］，生成的結(jié)果為橢圓，即標(biāo)準(zhǔn)差橢圓，通過計算其重心、方位角、長短半軸之比，揭示研究區(qū)耕地非農(nóng)化的擴散路徑［26］。耕地非農(nóng)化的重心、方位角計算公式［27］如下：

X=ni=1PiXi/ni=1Pi（1）

Y=ni=1PiYi/ni=1Pi（2）

θ=arctan［（Yt+m-Yt）/（Xt+m-Xt）］（3）

式中：X、Y表示研究區(qū)耕地非農(nóng)化重心；Xi、Yi表示第i個單元的幾何重心；Pi表示i單元的耕地非 農(nóng)化面積；n表示研究區(qū)的單元數(shù)量；θ表示重心遷移的角度；Xt、Yt與Xt+m、Yt+m分別表示第t年及t+m年的重心。

3 結(jié)果與分析

3.1 耕地非農(nóng)化遷移路徑

根據(jù)公式（1）～（3），運用 ArcGIS 空間分析工具提取研究區(qū)各時段的耕地非農(nóng)化標(biāo)準(zhǔn)差橢圓并計算所對應(yīng)的重心（圖1）、長短半軸、方位角等參數(shù)（表1），將各時間段的重心相連得到2009—2021年耕地非農(nóng)化重心曲線。由圖1、表1可知，丹江口市耕地非農(nóng)化方向呈現(xiàn)中部—東南—西北—西北的遷移路徑。各時段的長短軸之比由2009—2012年的1.697增加到2012—2015年的2.094，隨后逐漸減小，且越來越接近1，說明耕地非農(nóng)化

方向性越來越弱。2009—2021年方位角有一定的擺動性，2009—2012年較2012—2015年方位角轉(zhuǎn)動較明顯，呈順時針方向轉(zhuǎn)動，2012—2015年到2015—2018年方位角由88.789°減至66.849°，方位角轉(zhuǎn)動相較其他時段明顯，由順時針轉(zhuǎn)動轉(zhuǎn)變?yōu)槌誓鏁r針方向轉(zhuǎn)動，2018年之后方位角轉(zhuǎn)動微弱，轉(zhuǎn)動方向為逆時針。以2009—2012年耕地非農(nóng)化重心為參考起始點，丹江口市耕地非農(nóng)化遷移距離在2009—2021年變化情況為10.578 km—9.255 km—0.707 km，呈逐年減小趨勢，表征了研究區(qū)耕地非農(nóng)化空間不均衡性逐年減弱。

3.2 各鄉(xiāng)鎮(zhèn)耕地非農(nóng)化空間變化

為進(jìn)一步分析丹江口市各鄉(xiāng)鎮(zhèn)耕地非農(nóng)化的空間變化情況，采用 ArcGIS 的自然斷點法工具，將2009—2021年丹江口市各鄉(xiāng)鎮(zhèn)非農(nóng)化面積劃分為低值區(qū)、中值區(qū)和高值區(qū) 3 個類型區(qū)，以表征不同時間段研究區(qū)耕地非農(nóng)化的空間特征［27］。圖2反映了不同時段研究區(qū)各鄉(xiāng)鎮(zhèn)耕地非農(nóng)化的空間分布差異性。從圖2可見，丹江口市耕地非農(nóng)化空間差異顯著，整體而言，高值區(qū)集中分布于東部及西部地區(qū)，東部地區(qū)主要包括三官殿辦事處、丹趙路辦事處以及土關(guān)埡鎮(zhèn)，西部地區(qū)主要是六里坪鎮(zhèn)、官山鎮(zhèn)、武當(dāng)山旅游經(jīng)濟特區(qū)，2018—2021年偏北的習(xí)家店鎮(zhèn)非農(nóng)化面積轉(zhuǎn)變成高值區(qū)，2009—2012年耕地非農(nóng)化處于高值區(qū)的鄉(xiāng)鎮(zhèn)僅有六里坪鎮(zhèn)和三官殿辦事處，中值區(qū)包括官山鎮(zhèn)、武當(dāng)山旅游經(jīng)濟特區(qū)、均縣鎮(zhèn)以及浪河鎮(zhèn)4個鄉(xiāng)鎮(zhèn)，其他均為低值區(qū)。2012—2015年在前一時段的基礎(chǔ)上高值區(qū)增加了丹趙路辦事處，中值區(qū)中官山鎮(zhèn)與均縣鎮(zhèn)變?yōu)榈椭祬^(qū)，但增加了丁家營鎮(zhèn)與新港經(jīng)濟開發(fā)管理處，其他地方均為低值區(qū)。2015—2018年耕地非農(nóng)化面積最小，高值區(qū)僅有六里

坪鎮(zhèn)，中值區(qū)有8個鄉(xiāng)鎮(zhèn)。2018—2021年非農(nóng)化面積最大，僅有6個鄉(xiāng)鎮(zhèn)位于低值區(qū)，其他鄉(xiāng)鎮(zhèn)都處于中值區(qū)及高值區(qū)。

4 結(jié)論與討論

4.1 結(jié)論

基于2009、2012、2015、2018、2021年丹江口市國土變更調(diào)查數(shù)據(jù)，運用 GIS 空間分析工具，剖析了丹江口市耕地非農(nóng)化過程，并通過利菲弗方向性分布構(gòu)建標(biāo)準(zhǔn)差橢圓探明了丹江口市耕地非農(nóng)化空間格局及遷移路徑。主要結(jié)論如下：

（1）丹江口市耕地非農(nóng)化方向呈現(xiàn)中部—東南—西北—西北的遷移路徑。各時段的長短軸之比先增加后逐漸減小且越來越接近1，表明耕地非農(nóng)化方向性越來越弱。方位角呈現(xiàn)順時針—逆時針—逆時針特征，具有一定的擺動性。遷移距離由2009—2012年與2012—2015年之間的10.578 km逐年減小至2018—2021年的0.707 km。

（2）丹江口市各鄉(xiāng)鎮(zhèn)耕地非農(nóng)化空間存在一定的差異性。整體而言，高值區(qū)集中分布于東部及西部地區(qū)，東部地區(qū)主要包括三官殿辦事處、丹趙路辦事處以及土關(guān)埡鎮(zhèn)，西部地區(qū)主要是六里坪鎮(zhèn)、官山鎮(zhèn)、武當(dāng)山旅游經(jīng)濟特區(qū)。

4.2 討論

耕地的合理利用與保護(hù)關(guān)乎糧食安全問題，在城鎮(zhèn)化大力推進(jìn)的過程中，區(qū)域耕地非農(nóng)化問題愈來愈嚴(yán)峻。筆者研究了2009—2021年丹江口市耕地非農(nóng)化的遷移路徑，分析了整個縣域非農(nóng)化的遷移特征及各鄉(xiāng)鎮(zhèn)的非農(nóng)化差異，對研究區(qū)耕地的保護(hù)與精準(zhǔn)化管理及國土空間規(guī)劃等方面提供了借鑒與參考，豐富了耕地保護(hù)制度轉(zhuǎn)型的理論。然而，該研究僅對丹江口市縣域尺度耕地非農(nóng)化動態(tài)遷移路徑進(jìn)行研究，而耕地非農(nóng)化受自然、人文、經(jīng)濟、政策等因素的影響，未來需要在耕地非農(nóng)化及其變化規(guī)律方面進(jìn)行深度剖析，明晰耕地非農(nóng)化時空格局演變的有效驅(qū)動機制。

參考文獻(xiàn)

［1］ 自然資源部辦公廳.自然資源部辦公廳關(guān)于開展2022年度全國國土變更調(diào)查工作的通知：自然資辦發(fā)〔2022〕49號［EB/OL］.（2022-11-25）［2022-12-10］.http：∥gi.mnr.gov.cn/202211/t20221128_2769039.html.

［2］ 錢鳳魁，逄然然，于洋，等.遼寧省不同地貌區(qū)耕地景觀格局與耕地質(zhì)量空間關(guān)聯(lián)特征及作用機制［J］.中國生態(tài)農(nóng)業(yè)學(xué)報，2023，31（1）：113-124.

［3］ 范影影，朱洪芬，劉慧芳，等.基于面向?qū)ο蠓指畹母丶羞B片經(jīng)營分布提取及其影響因素研究：以大同市云州區(qū)四鄉(xiāng)鎮(zhèn)為例［J］.農(nóng)業(yè)資源與環(huán)境學(xué)報，2023，40（3）：548-556.

［4］ 黃珺嫦，華軒軻，賈夢瑤，等.2000—2020年河南省耕地利用顯性轉(zhuǎn)型時空特征［J］.水土保持通報，2023，43（1）：206-213.

［5］ 余李敏，李陽兵，羅光杰，等.貴州中部巖溶山區(qū)坡耕地時空演變特征及驅(qū)動因素分析［J］.農(nóng)業(yè)工程學(xué)報，2022，38（11）：289-298.

［6］ 王麗蒙，王加勝，楊昆，等.1990—2019年石漠化區(qū)耕地時空特征及成因分析：以富源縣為例［J］.測繪通報，2022（4）：77-82.

［7］ 呂瑩瑩，王智，夏欣，等.內(nèi)蒙古大興安嶺林草交錯區(qū)耕地景觀演變及空間過程［J］.生態(tài)學(xué)報，2023，43（3）：1209-1218.

［8］ 金一諾，王楓.珠三角耕地利用轉(zhuǎn)型的時空演化特征及其影響因素分析［J］.水土保持研究，2022，29（4）：352-361.

［9］ 侯青青，陳英，裴婷婷，等.近25 a來甘肅省耕地資源時空變化及其影響因子［J］.干旱區(qū)研究，2022，39（3）：955-967.

［10］ 陳松林，趙翔，劉中秋，等.2009—2018年湖南省縣域耕地時空演變及其驅(qū)動力［J］.水土保持通報，2022，42（3）：274-283.

［11］ 杜國明，薛濡壕，于鳳榮.耕地集約利用轉(zhuǎn)型的理論解析［J］.資源科學(xué)，2022，44（3）：425-435.

［12］ 胡賢輝，劉蒙罷，文高輝.中國耕地集約利用與生態(tài)效率耦合協(xié)調(diào)時空分異特征研究［J］.長江流域資源與環(huán)境，2022，31（10）：2282-2294.

［13］ 郝軍，吳曉光，劉佳雨，等.生態(tài)文明背景下內(nèi)蒙古耕地發(fā)展對策研究［J］.干旱區(qū)資源與環(huán)境，2022，36（12）：49-56.

［14］ 王慧，龍開勝.基于耕地多功能的江蘇省耕地保護(hù)補償標(biāo)準(zhǔn)研究［J］.中國農(nóng)業(yè)資源與區(qū)劃，2022，43（11）：101-111.

［15］ 梁鑫源，金曉斌，孫瑞，等.多情景糧食安全底線約束下的中國耕地保護(hù)彈性空間［J］.地理學(xué)報，2022，77（3）：697-713.

［16］ 何新瑩，聶艷，王樸，等.基于改進(jìn)灰靶模型的耕地質(zhì)量評價方法與實證［J］.土壤學(xué)報，2023，60（4）：1007-1016.

［17］ 王來剛，郭燕，賀佳，等.遙感數(shù)據(jù)輔助下縣域耕地質(zhì)量評價與空間分布研究［J］.中國農(nóng)業(yè)資源與區(qū)劃，2022，43（12）：137-146.

［18］ 張友，劉玉.鄉(xiāng)村多功能視角下耕地資源資產(chǎn)價值核算研究［J］.中國農(nóng)業(yè)資源與區(qū)劃，2022，43（4）：129-138.

［19］ 謝立軍，白中科，楊博宇，等.襄垣縣耕地資源價值核算及空間格局分析［J］.中國農(nóng)業(yè)大學(xué)學(xué)報，2022，27（6）：204-214.

［20］ 吳昊玥，孟越，黃瀚蛟，等.中國耕地低碳利用績效測算與時空分異［J］.自然資源學(xué)報，2022，37（5）：1148-1163.

［21］ 毛亮軍，王淵，王新生，等.丹江口水庫消落帶不同作物類型的土壤總氮分析［J］.中國農(nóng)學(xué)通報，2018，34（8）：78-82.

［22］ 李宏鵬，田達(dá)睿，譚靜斌.延安市2000—2020年耕地非農(nóng)化時空格局演變及其影響因素［J］.水土保持通報，2022，42（4）：330-337，372.

［23］ 高于.山東完善耕地保護(hù)獎懲補償機制 給予成效突出的20個縣（市、區(qū)）6000畝新增建設(shè)用地指標(biāo)和3個億的資金獎勵［EB/OL］.（2022-03-15）［2022-11-17］.https：∥www.163.com/dy/article/H2GDTS7F05 14CFC7.html.

［24］ 李旭卉，楊瑾.西安市高陵區(qū)建設(shè)用地擴張時空變化分析［J］.農(nóng)業(yè)與技術(shù)，2020，40（15）：36-41.

［25］ 白永平，紀(jì)學(xué)朋.蘭州市牛肉面館空間分布特征及其影響因素［J］.開發(fā)研究，2017（1）：128-133.

［26］ 吳海中，田曉四，張樂勤.近40年安徽省耕地非農(nóng)化時空分布特征及擴散路徑［J］.中國農(nóng)業(yè)資源與區(qū)劃，2023，44（6）：110-118.

［27］ 劉少坤，王嘉佳，林樹高，等.廣西邊境地區(qū)耕地非農(nóng)化的空間特征與遷移路徑［J］.中國農(nóng)業(yè)資源與區(qū)劃，2022，43（10）：162-173.