孫佩 尹偉 易潔偉 康全國(guó) 張馳 袁知洋

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

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

中圖分類號(hào) F 323.211? 文獻(xiàn)標(biāo)識(shí)碼 A? 文章編號(hào) 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

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

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

收稿日期 2022-12-11

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

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

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

1 研究區(qū)概況

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

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

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

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

2.2 研究方法

利菲弗方向性分布工具在空間統(tǒng)計(jì)中是綜合能力最突出的工具之一，有著廣泛的應(yīng)用，該算法在1926年最早由美國(guó)南加州大學(xué)社會(huì)學(xué)教授韋爾蒂·利菲弗（D.Welty Lefever）提出［24］， 可以用來(lái)度量概括地理要素中心趨勢(shì)、離散和方向趨勢(shì)的空間特征［25］，生成的結(jié)果為橢圓，即標(biāo)準(zhǔn)差橢圓，通過計(jì)算其重心、方位角、長(zhǎng)短半軸之比，揭示研究區(qū)耕地非農(nóng)化的擴(kuò)散路徑［26］。耕地非農(nóng)化的重心、方位角計(jì)算公式［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個(gè)單元的幾何重心；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），運(yùn)用 ArcGIS 空間分析工具提取研究區(qū)各時(shí)段的耕地非農(nóng)化標(biāo)準(zhǔn)差橢圓并計(jì)算所對(duì)應(yīng)的重心（圖1）、長(zhǎng)短半軸、方位角等參數(shù)（表1），將各時(shí)間段的重心相連得到2009—2021年耕地非農(nóng)化重心曲線。由圖1、表1可知，丹江口市耕地非農(nóng)化方向呈現(xiàn)中部—東南—西北—西北的遷移路徑。各時(shí)段的長(zhǎng)短軸之比由2009—2012年的1.697增加到2012—2015年的2.094，隨后逐漸減小，且越來(lái)越接近1，說明耕地非農(nóng)化

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

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

為進(jìn)一步分析丹江口市各鄉(xiāng)鎮(zhèn)耕地非農(nóng)化的空間變化情況，采用 ArcGIS 的自然斷點(diǎn)法工具，將2009—2021年丹江口市各鄉(xiāng)鎮(zhèn)非農(nóng)化面積劃分為低值區(qū)、中值區(qū)和高值區(qū) 3 個(gè)類型區(qū)，以表征不同時(shí)間段研究區(qū)耕地非農(nóng)化的空間特征［27］。圖2反映了不同時(shí)段研究區(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)濟(jì)特區(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)濟(jì)特區(qū)、均縣鎮(zhèn)以及浪河鎮(zhèn)4個(gè)鄉(xiāng)鎮(zhèn)，其他均為低值區(qū)。2012—2015年在前一時(shí)段的基礎(chǔ)上高值區(qū)增加了丹趙路辦事處，中值區(qū)中官山鎮(zhèn)與均縣鎮(zhèn)變?yōu)榈椭祬^(qū)，但增加了丁家營(yíng)鎮(zhèn)與新港經(jīng)濟(jì)開發(fā)管理處，其他地方均為低值區(qū)。2015—2018年耕地非農(nóng)化面積最小，高值區(qū)僅有六里

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

4 結(jié)論與討論

4.1 結(jié)論

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

（1）丹江口市耕地非農(nóng)化方向呈現(xiàn)中部—東南—西北—西北的遷移路徑。各時(shí)段的長(zhǎng)短軸之比先增加后逐漸減小且越來(lái)越接近1，表明耕地非農(nóng)化方向性越來(lái)越弱。方位角呈現(xiàn)順時(shí)針—逆時(shí)針—逆時(shí)針特征，具有一定的擺動(dòng)性。遷移距離由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)濟(jì)特區(qū)。

4.2 討論

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

參考文獻(xiàn)

［1］ 自然資源部辦公廳.自然資源部辦公廳關(guān)于開展2022年度全國(guó)國(guó)土變更調(diào)查工作的通知：自然資辦發(fā)〔2022〕49號(hào)［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ī)制［J］.中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2023，31（1）：113-124.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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