敖寧玉 鮑安紅 胡嫚 張森林

摘要：宜機(jī)化改造是解決丘陵山區(qū)機(jī)械化發(fā)展不平衡的必由之路，在改造之前進(jìn)行適宜性評價能有效提高改造過程中的科學(xué)性和針對性。選取地塊形狀指數(shù)、連片性、農(nóng)村建筑、交通條件、地表粗糙度、地形起伏度6個宜機(jī)化限制因子，以重慶市開州區(qū)一個待改造項目為例，結(jié)合層次分析法（Analytic Hierarchy Process，AHP）和K-means聚類分析方法以及模糊理論，以ArcGIS為平臺構(gòu)建丘陵山地宜機(jī)化改造適宜性綜合評價模型。結(jié)果表明，各影響因子的權(quán)重排序為：地形起伏度>地表粗糙度>連片性>形狀指數(shù)>交通條件>農(nóng)村建筑，為后續(xù)的加權(quán)疊加分析提供基礎(chǔ)依據(jù)；評價目標(biāo)的宜機(jī)化改造適宜性綜合評價結(jié)果為比較適宜，其中適宜、比較適宜、比較不適宜、不適宜的區(qū)域分別占改造區(qū)域面積的41.87%、44.79%、12.8%、0.54%?；贏rcGIS的適宜性評價方法能有效進(jìn)行丘陵耕地適宜性評價，提高宜機(jī)化改造效率，加快丘陵耕地機(jī)械化進(jìn)程。

關(guān)鍵詞：丘陵山地；宜機(jī)化；層次分析法；K-means聚類算法；模糊綜合評價；適宜性評價

中圖分類號：S284? 文獻(xiàn)標(biāo)識碼：A? 文章編號：20955553 （2024） 03-0260-09

Research on the suitable mechanization evaluation method of terraced farmlandin hilly mountainous areas based on ArcGIS

Ao Ningyu， Bao Anhong， Hu Man， Zhang Senlin

（College of Engineering and Technology， ?South West University， Chongqing， 400715， China）

Abstract：

The suitable mechanization transformation ?is the only way to address the imbalance of the mechanization development in hilly and mountainous areas. ?The suitability evaluation before initiating the transformation can effectively increase scientificity and targeted effectiveness throughout the transformation process. In this study， six suitability limiting factors including shape index， connectivity， rural buildings， transportation condition， surface roughness， and terrain variability were selected， taking a project in Kaizhou District of Chongqing as an example， a comprehensive evaluation model for terrace mechanization transformation suitability was constructed using Analytic Hierarchy Process（AHP） method， K-means clustering analysis and fuzzy theory on the platform of ArcGIS. The results indicated that the weight ranking of each influencing factor was as follows： terrain variability > surface roughness > connectivity >shape index > transportation condition> rural buildings. This result provided a basis for subsequent weighted overlay analysis. The comprehensive evaluation result of suitability for the transformation objective was relatively suitable， with suitable， relatively suitable， relatively unsuitable， and unsuitable areas accounting for 41.87%， 44.79%， 12.8%， and 0.54% of the transformation area， respectively. The suitability evaluation method based on ArcGIS can effectively evaluate the suitability of hilly farmland， improve the efficiency of mechanized transformation， accelerate ?the process of mechanization of hilly farmland.

Keywords：hilly and mountainous; suitable mechanization; analytic hierarchy process; K-means clustering algorithm; fuzzy comprehensive evaluation; suitability evaluation

0 引言

農(nóng)業(yè)機(jī)械化是實現(xiàn)鄉(xiāng)村振興的重要手段，隨著綜合國力的提升，我國農(nóng)業(yè)機(jī)械化進(jìn)程顯著，但區(qū)域發(fā)展不平衡問題較為明顯，尤其是丘陵地區(qū)由于地形限制機(jī)械化發(fā)展極為緩慢［1］。重慶是典型的丘陵山地，耕地2800多萬畝，由于地形影響，坡耕地比重較大且破碎化程度較高［2］，導(dǎo)致綜合機(jī)械化水平比平原地區(qū)顯著偏低，嚴(yán)重影響了機(jī)械化進(jìn)程的順利推進(jìn)。因此，如何提高丘陵山地農(nóng)業(yè)機(jī)械化效率成為限制重慶農(nóng)業(yè)機(jī)械化發(fā)展亟待解決的問題。實踐證明發(fā)明適宜丘陵山地的小型“微耕機(jī)”不能從根本上提高丘陵耕地的機(jī)械化效率［3］，有研究指出，改地適機(jī)才是提高丘陵山地機(jī)械化的根本出路［4］，開展農(nóng)田宜機(jī)化改造是破解重慶市農(nóng)業(yè)機(jī)械化和農(nóng)業(yè)現(xiàn)代化發(fā)展最大短板的有效途徑。

日本、韓國的山地丘陵面積占國土面積的70%～80%，通過宜機(jī)化改造了機(jī)械化程度達(dá)到100%［5， 6］，相關(guān)文獻(xiàn)表明國外的土地整理更加重視生態(tài)環(huán)境的保護(hù)［7， 8］。目前國內(nèi)以宜機(jī)化為對象的研究還處在探索階段，已有的研究多側(cè)重于土地整治后的效益評價［9， 10］，牛坡等［11］對宜機(jī)化進(jìn)行了綜合效益評價，是對宜機(jī)化改造的定量評價，開創(chuàng)了宜機(jī)化相關(guān)定量研究的先河，李伶俐等［12］首次提出基于地形復(fù)雜度的宜機(jī)化改造適宜性評價。劉光盛等［13］運用改進(jìn)障礙度模型對宜機(jī)化整治難度進(jìn)行分區(qū)，鐘守琴［14］、胡義萍［15］等對宜機(jī)化改造過程中田塊修筑工程參數(shù)、田間道路設(shè)計等具體操作標(biāo)準(zhǔn)進(jìn)行了探討?，F(xiàn)有研究在推動宜機(jī)化發(fā)展方面發(fā)揮了一定作用，但是缺乏一個完整的評價體系和宜機(jī)化適宜性評價的模型。

基于此，本文以重慶市開州區(qū)一個待改造項目為評價目標(biāo)，選擇地塊形狀指數(shù)、連片性、農(nóng)村建筑、交通條件、地表粗糙度和地形起伏度等為限制因子，以模糊理論為基礎(chǔ)，結(jié)合AHP和K-means聚類分析方法，以ArcGIS為平臺對研究區(qū)域宜機(jī)化改造適宜性進(jìn)行評價，為提高宜機(jī)化改造效率提供科學(xué)依據(jù)，加快推進(jìn)我國農(nóng)業(yè)機(jī)械化發(fā)展和農(nóng)機(jī)裝備產(chǎn)業(yè)轉(zhuǎn)型升級。

1 研究方法

本文所用數(shù)據(jù)主要來源于所研究項目分辨率為30m的數(shù)字高程影像（DEM）。研究區(qū)域中的耕地、道路、農(nóng)村建筑區(qū)域數(shù)據(jù)均根據(jù)數(shù)字高程影像圖進(jìn)行繪制，并使用ArcGIS將測繪所得數(shù)據(jù)轉(zhuǎn)換為柵格數(shù)據(jù)格式進(jìn)行分析。

本文的研究思路分為以下3個步驟：（1）從宜機(jī)化改造要求出發(fā)，從立地條件、基礎(chǔ)設(shè)施、自然特征3個方面構(gòu)建宜機(jī)化改造適宜性評價指標(biāo)體系；（2）運用AHP和K-means聚類分析方法計算各指標(biāo)權(quán)重；（3）以ArcGIS為平臺結(jié)合模糊數(shù)學(xué)的方法構(gòu)建宜機(jī)化改造適宜性綜合評價模型對研究區(qū)域進(jìn)行綜合評價。

1.1 模糊綜合評價法

基于模糊理論的模糊綜合評價法是運用模糊數(shù)學(xué)原理和最大隸屬度原則，考慮與評價對象相關(guān)的各項影響因素，對評價對象做出綜合評價，具有評價過程清晰，結(jié)果明確的優(yōu)點，適于丘陵山地“宜機(jī)化”適宜性評價。

1.1.1 指標(biāo)體系構(gòu)建

丘陵山地宜機(jī)化改造需要考慮多種限制因素的綜合影響，各種影響因素紛繁復(fù)雜，但總體而言，影響丘陵山地宜機(jī)化改造的主要因素分為立地條件μ1、基礎(chǔ)設(shè)施μ2、自然特征μ3三個方面。本文綜合考慮三方面因素影響，結(jié)合建設(shè)高標(biāo)準(zhǔn)農(nóng)田評價指標(biāo)和《丘陵山區(qū)‘宜機(jī)化地塊整理整治技術(shù)規(guī)范》《丘陵山區(qū)坡改梯宜機(jī)化土地整治技術(shù)規(guī)范》，本著指標(biāo)選取的科學(xué)合理性并遵循主導(dǎo)性、分析可行性、全面性、獨立性、差異性等原則，建立適合于丘陵山地宜機(jī)化改造的評價指標(biāo)體系。

為了在簡化整體評價過程的同時能更科學(xué)準(zhǔn)確而全面的評價丘陵山地宜機(jī)化改造適宜性，本文從宜機(jī)化改造的要求和特點中篩選出評價單元地塊形狀指數(shù)μ11和地塊連片性μ12兩個“立地條件”指標(biāo)，“基礎(chǔ)設(shè)施”中選取農(nóng)村建筑μ21和交通條件μ22兩個評價指標(biāo)，從“地塊特征”中選取地表粗糙度μ31和地形起伏度μ32兩個指標(biāo)。

地塊的形狀指數(shù)和連片性主要通過影響地塊的破碎度來影響宜機(jī)化改造，地塊破碎度在一定程度上能反應(yīng)研究單元地塊的完整性和連片性，在土地整治方面起著重要作用［16］。張正峰等［17］研究發(fā)現(xiàn)，對地塊的合理改造能降低地塊的形狀指數(shù)和分形指數(shù)，通過改變地塊規(guī)模，可以提高單個地塊規(guī)模，消除地塊死角，提高地塊的連片性，降低地塊的破碎度，從而提高機(jī)械耕作效率。農(nóng)機(jī)是否能到達(dá)耕地取決于耕地到交通干線以及田間道路的距離，在進(jìn)行宜機(jī)化改造時要避免和農(nóng)村建筑用地沖突，并在規(guī)劃設(shè)計機(jī)耕道時要將農(nóng)村建筑考慮在內(nèi)。地表粗糙度［18］和地形起伏度是反映研究區(qū)域內(nèi)地形起伏的宏觀因子，研究區(qū)域內(nèi)的坡度和高差越大，宜機(jī)化改造難度越大。這些指標(biāo)不僅包括評價單元本身的地塊特征限制因素，而且從地塊的自然特征出發(fā)，綜合考慮了非特征限制因素，能夠科學(xué)有效的反映出丘陵山地宜機(jī)化改造的適宜程度。

根據(jù)評價因子建立因素集U，本文將評價指標(biāo)體系分為U={U1，U2，U3}三個一級指標(biāo)，分別表示立地條件、基礎(chǔ)設(shè)施、自然特征，將三個一級指標(biāo)分別細(xì)分為2個，共6個二級指標(biāo)U1={u11，u12}，U2={u21，u22}，U3={u31，u32}。耕地宜機(jī)化改造適宜性評價的指標(biāo)體系如表1所示。計算公式如式（1）～式（6）所示。

I=L/S（1）

式中：I——形狀指數(shù)；L——地塊周長，m；S——地塊面積。

D=1（radius）2∑ni=13π·popi1-distiradius22（2）

式中：D——核密度；i——耕地圖斑中心點；popi——點i的population字段值；radius——核半徑范圍；disti——點i和（x，y）位置之間的距離。

（x-t1）2+（y-t2）2=r2（3）

式中：t1，t2——農(nóng)村建筑圖斑中心點坐標(biāo)；r——緩沖區(qū)半徑。

Distance=（x2-x1）2+（y2-y1）（4）

式中：Distance——兩點之間的歐式距離；（x1，y1）、（x2，y2）——圖層中兩點的坐標(biāo)。R=1/cosα（5）

式中：R——地形粗糙度；α——地形坡度。

Z=Amax-Amin（6）

式中：Z——特定區(qū)域地形起伏度；Amax——該區(qū)域內(nèi)最大高程值；Amin——該區(qū)域內(nèi)最小高程值。

1.1.2 評價集

不同的影響因素會對綜合評價結(jié)果產(chǎn)生不同的影響，因此需要將各影響因子根據(jù)數(shù)值對應(yīng)一個評語集。根據(jù)相關(guān)的技術(shù)規(guī)定結(jié)合工程實際情況構(gòu)建的各影響因子評價集如表2所示。

1.1.3 權(quán)重集

本文選用AHP法結(jié)合K-means聚類算法進(jìn)行各因子權(quán)重的計算，層次分析法（AHP）最早由美國Satty提出，經(jīng)過多年的發(fā)展已經(jīng)逐漸成熟，被廣泛運用于工程評價領(lǐng)域，但是其具有主觀性較大的缺點，因此本文選用K-means聚類算法增強(qiáng)指標(biāo)權(quán)重計算的客觀性。K-means聚類算法又稱為K均值聚類方法，核心方法是根據(jù)樣本數(shù)據(jù)點之間的歐氏距離構(gòu)造一個空間度量矩陣，距離較近的樣本數(shù)據(jù)點視為一簇，一個簇中的數(shù)據(jù)點具有較大的相似性［19］。

為更加真實地計算重慶市丘陵山區(qū)耕地宜機(jī)化改造適宜性各評價指標(biāo)的權(quán)重值，本文共邀請15位相關(guān)領(lǐng)域?qū)＜覍υu價指標(biāo)進(jìn)行重要性賦值，將15位專家給出的判別矩陣計算得到各評價指標(biāo)權(quán)重值，進(jìn)行一致性檢驗之后作為樣本構(gòu)建一個空間矩陣進(jìn)行K-means聚類分析，根據(jù)K-means聚類分析結(jié)果最終將15個樣本分為了兩簇，選取簇內(nèi)平方誤差和SSE較小的聚類中心作為最終的權(quán)重向量。權(quán)重確定方法主要分為以下3個步驟。

1） 構(gòu)建判別矩陣。

構(gòu)建某一層次指標(biāo)相對于上一層次指標(biāo)的判別矩陣

A=1…a1j…a1nai1…aij…ainan1…anj…1

式中：aij——i因子與j因子相比的重要程度；n——第n個因子。

上述判別矩陣的構(gòu)建根據(jù)AHP評價尺度表如表3所示。

2） 一致性檢驗。

將上訴判別矩陣A按列進(jìn)行歸一化處理得到矩陣C，對矩陣C按行相加后得到矩陣

C-再進(jìn)行歸一化處理得到權(quán)向量W，計算矩陣的最大特征根λmax，如式（7）、式（8）所示。

AW=λmaxW（7）

λmax=∑ni=1（AW）inWi（8）

然后進(jìn)行一致性檢驗，一致性檢驗公式如式（9）、式（10）所示。

CI=λmax-nn-1（9）

CR=CIRI（10）

式中：CI——一致性指標(biāo)；CR——一致性檢驗率；RI——隨機(jī)一致性指標(biāo)，本文n=6，查表得RI=1.24。

當(dāng)CR≤0.1時，權(quán)向量W即為某一層次指標(biāo)對于上一層次相關(guān)指標(biāo)的權(quán)重值，當(dāng)CR＞0.1時，則沒通過一致性檢驗，需要對判別矩陣進(jìn)行修正，直到通過一致性檢驗為止。

3） K-means聚類分析。

K-means聚類分析采用計算樣本間的歐氏距離作為衡量樣本間相似性的標(biāo)準(zhǔn)，將AHP方法計算得到的權(quán)重作為K-means聚類分析時的權(quán)值空間矩陣，將用AHP法計算的各指標(biāo)權(quán)向量作為數(shù)據(jù)樣本集X={x1，x2，x3，…，xm}，m表示第m個權(quán)重向量，將樣本劃分K個簇C={C1，C2，C3，…，CK}，計算簇Ci的質(zhì)心，如式（11）所示。

μi=1|Ci|∑x∈Cix（11）

計算每個簇中所有樣本數(shù)據(jù)的平方誤差和SSE，平方誤差和SSE越小，簇內(nèi)樣本數(shù)據(jù)間的相似度越高，計算公式如式（12）所示。

SSE=∑ki=1∑x∈Ci‖x-μi‖2（12）

通過計算Xi與質(zhì)心之間的歐氏距離將數(shù)據(jù)集中的每個數(shù)據(jù)點分配到距離最近的聚類中心，然后重新計算每個聚類中心的數(shù)據(jù)平均值，并將該均值作為新的聚類中心，不斷重復(fù)上述過程，直到計算所得的新聚類中心與舊聚類中心不再發(fā)生變化或者他們之間的最小距離小于設(shè)置的閾值，聚類結(jié)果不再發(fā)生變化即認(rèn)為算法收斂。K-means聚類算法在VScode軟件中通過Python語言實現(xiàn)，實現(xiàn)流程如圖1所示。

通過上述步驟計算的結(jié)果如表4所示。根據(jù)K-means聚類分析，將15位專家所給出的權(quán)重值分為兩個簇，其中編號為1、2、3、4、5、6、11、12、13、14、15共11個權(quán)重向量分為簇1，編號為7、8、9、10共4個權(quán)重向量分為簇2，簇1和簇2的聚類中心值以及各簇的誤差平方和SSE如表5所示。

簇1包含了11個數(shù)據(jù)樣本，且簇內(nèi)誤差平方和最小，因此選擇簇一聚類中心作為評價指標(biāo)權(quán)重值λ=［0.101，0.195，0.026，0.081，0.240，0.357］

1.1.4 建立模糊隸屬函數(shù)

不同的評價因子對宜機(jī)化改造評價的影響不同，需要建立不同的隸屬度函數(shù)來表示各因子的模糊隸屬度，根據(jù)影響因素值與評價結(jié)果關(guān)系，選用梯形型隸屬度函數(shù)來映射形狀指數(shù)、連片性、交通條件、地表粗糙度以及地形起伏度的模糊隸屬度，選擇矩陣型函數(shù)來映射農(nóng)村建筑的隸屬度函數(shù)。隸屬度函數(shù)表示如表6所示。各項評價指標(biāo)得隸屬度函數(shù)圖整理如圖2所示。

1.2 基于ArcGIS的宜機(jī)化適宜性模糊評價

ArcGIS具有較強(qiáng)的空間數(shù)據(jù)分析與查詢能力，基于ArcGIS的柵格加權(quán)疊加板塊，便可以高效、直觀的得出區(qū)域范圍內(nèi)的評價結(jié)果，具體評價步驟如下：將研究區(qū)域的高程數(shù)據(jù)導(dǎo)入ArcGIS中利用Arctoolbox中的坡度分析工具以高程柵格圖為基礎(chǔ)生成坡度柵格圖像Slope of DEM，隨后使用柵格計算器計算生成地表粗糙度柵格圖像；運用焦點統(tǒng)計分析中的range函數(shù)生成地形起伏度柵格圖層，根據(jù)影像圖在ArcGIS中繪制耕地、道路、農(nóng)村建筑圖斑并將其轉(zhuǎn)換為柵格數(shù)據(jù)格式，隨后分別運用核密度分析、歐式距離分析、緩沖區(qū)分析得到連片性、交通條件、農(nóng)村建筑的分析柵格圖像。完成單因素分析之后，運用柵格計算器根據(jù)不同評價標(biāo)準(zhǔn)下各評價因子的隸屬度函數(shù)生成適宜、比較適宜、比較不適宜、不適宜4項隸屬度柵格圖像。由于ArcGIS中無法進(jìn)行矩陣疊加，因此將柵格隸屬度按照適宜、比較適宜、比較不適宜、不適宜4項評價等級分別進(jìn)行加權(quán)疊加，即可得到6項評價因子所對應(yīng)的4個評價等級下的單項評價結(jié)果，對單項評價結(jié)果進(jìn)行最大隸屬度判別后進(jìn)行柵格重分類得到最終的綜合評價結(jié)果。

2 結(jié)果與分析

本文以重慶市開州區(qū)一宜機(jī)化改造項目為例，驗證宜機(jī)化改造適宜性評價模型。該區(qū)域呈南北向帶狀分布，總面積為101 312.76m2，最高處海拔為319.1m，最低處海拔為389.4m。去除研究區(qū)域內(nèi)的居民建筑等不能進(jìn)行宜機(jī)化改造的區(qū)域，對剩余73 922.56m2的地塊按上述模型進(jìn)行宜機(jī)化改造適宜性評價。

2.1 影響因子分析

2.1.1 立地條件分析

按照研究區(qū)域的衛(wèi)星影像圖中的原始地塊分布將研究區(qū)域內(nèi)可改造地塊分為68個圖斑，在ArcGIS中將數(shù)據(jù)轉(zhuǎn)換成柵格數(shù)據(jù)格式，利用字段計算器和核密度分析工具計算分析研究區(qū)域的地塊連片性和形狀指數(shù)。結(jié)果如圖3所示。

研究范圍內(nèi)的地塊集中分布在3個區(qū)域，3個區(qū)域間也具有一定的相連性，在進(jìn)行宜機(jī)化改造時比較適宜。研究區(qū)域的原始地塊多呈不規(guī)則條狀分布，形狀指數(shù)I≤6的地塊面積為50224.05m2，占待改造面積的67.97%，從立地條件的連片性和形狀指數(shù)兩個指標(biāo)來看，此研究區(qū)域比較適宜進(jìn)行宜機(jī)化改造。

2.1.2 基礎(chǔ)設(shè)施

將研究區(qū)域內(nèi)的居民點作為圖斑進(jìn)行緩沖區(qū)分析，圖中設(shè)置的是5m的緩沖區(qū)域，即居民點圖斑5m內(nèi)的區(qū)域是禁止開發(fā)區(qū)域。將研究區(qū)域內(nèi)現(xiàn)有的主要道路進(jìn)行歐氏距離分析，由圖4分析可知，待改造耕地離現(xiàn)有道路最遠(yuǎn)距離為110m，現(xiàn)有道路30m內(nèi)的耕地面積為54673.56m2，占待改造面積的73.96%，區(qū)域內(nèi)地塊與田間道路、田間道路與外部路網(wǎng)之間連接情況較好，滿足機(jī)械作業(yè)和進(jìn)出通行需求，有利于進(jìn)行宜機(jī)化改造工程的實施。

2.1.3 自然特征

地表粗糙度與地形起伏度分析結(jié)果如圖5所示。

根據(jù)圖5（a）可知，研究區(qū)域內(nèi)地形坡度在（0°，7°］的面積為17949.75m2，地形坡度在（7°，15°］的面積為48713.66m2，根據(jù)宜機(jī)化改造技術(shù)規(guī)范，地形坡度在（0°，15°］，也就是地表粗糙度在（1，1.0353］之間適宜宜機(jī)化改造，研究區(qū)域內(nèi)適宜改造的地塊面積占總面積的90.18%。坡度在（0°，7°］的區(qū)域集中分布在研究區(qū)域的上部與下部，中部的地形坡度主要為（７°，１５°］，可為宜機(jī)化改造設(shè)計提供思路。由圖5（b）所示，研究區(qū)域內(nèi)的地形起伏度較大，起伏度在（０，３］范圍內(nèi)的區(qū)域面積僅占改造區(qū)域的37.21%。中部區(qū)域的起伏度較大，進(jìn)行宜機(jī)化改造時比較困難。

２.2 模糊綜合評價結(jié)果分析

利用ArcGIS中的柵格計算器對各評價因子?xùn)鸥癜凑赵u價標(biāo)準(zhǔn)和相應(yīng)的模糊隸屬函數(shù)進(jìn)行隸屬度賦值，再將四項評價標(biāo)準(zhǔn)隸屬度柵格根據(jù)計算出的權(quán)重值λ=［0.101，0.195，0.026，0.081，0.240，0.357］進(jìn)行加權(quán)疊加分析，根據(jù)權(quán)重計算結(jié)果可知，地表粗糙度與地形起伏度權(quán)重值分別為0.240和0.357，對宜機(jī)化改造適宜性評價結(jié)果影響最大；立地條件下的形狀指數(shù)和連片性對評價結(jié)果影響較大，分別為0.101和0.195，研究區(qū)域內(nèi)的基礎(chǔ)設(shè)施對宜機(jī)化綜合評價結(jié)果影響最小，農(nóng)村建筑和交通條件所占權(quán)重僅為0.026和0.081。加權(quán)疊加分析結(jié)果如圖6所示。

按照最大隸屬度原則對單項評價結(jié)果進(jìn)行最大隸屬度判別，得到整個評價區(qū)域的宜機(jī)化適宜性模糊綜合評價的結(jié)果為比較適宜，對“比較適宜”的隸屬度柵格按照隸屬度0

３ 結(jié)論

本文以重慶市開州區(qū)一個改造項目為目標(biāo)，從宜機(jī)化改造的要求和特點出發(fā)，建立了丘陵耕地宜機(jī)化適宜性評價指標(biāo)體系，并結(jié)合AHP和K-means聚類分析法更加客觀的確定了指標(biāo)權(quán)重，在ArcGIS中對評價因子數(shù)據(jù)采用緩沖區(qū)分析、核密度分析、加權(quán)疊加分析等方法進(jìn)行處理，以ArcGIS為平臺結(jié)合模糊綜合評價方法建立了丘陵山地農(nóng)田宜機(jī)化適宜性評價模型，對研究區(qū)域的宜機(jī)化改造適宜性進(jìn)行評價。

1） 不同的影響因子對于丘陵山地宜機(jī)化改造適宜性的影響力存在顯著差別，本文通過AHP法結(jié)合K-means聚類分析法計算出評價指標(biāo)權(quán)重值，計算結(jié)果顯示：代表地形特征的地表起伏度和地形粗糙度的因子權(quán)重值為0.357和0.240，對丘陵山地宜機(jī)化改造適宜性評價的影響最大；代表立地條件的連片性和形狀指數(shù)的權(quán)重值分別為0.195和0.101，影響相對較小；基礎(chǔ)設(shè)施農(nóng)村建筑和交通條件的權(quán)重值為0.026和0.081，對宜機(jī)化改造的影響最小，各影響因子的排序為：地形起伏度>地表粗糙度>連片性>形狀指數(shù)>交通條件>農(nóng)村建筑。

2） 本文以模糊綜合評價法為基礎(chǔ)，以ArcGIS為平臺搭建了丘陵山地宜機(jī)化改造適宜性綜合評價模型，以重慶市開州區(qū)某宜機(jī)化改造項目為例，生成各影響因素的柵格數(shù)據(jù)，并運用柵格計算器，加權(quán)疊加分析以及最大隸屬度判別得到該項目整體的宜機(jī)化改造適宜性評價結(jié)果為比較適宜。對比較適宜隸屬度柵格進(jìn)行柵格重分類得到研究區(qū)域中各具體地塊的綜合評價結(jié)果，分析結(jié)果顯示，適宜、比較適宜、比較不適宜、不適宜的區(qū)域分別占改造區(qū)域面積的41.87%、44.79%、12.8%、0.54%。該結(jié)果對后續(xù)的工程施工有一定的指導(dǎo)意義，在設(shè)計改造方案時重點考慮不適宜和比較不適宜區(qū)域，通過消除死角，聯(lián)通地塊，并小為大，合理規(guī)劃農(nóng)田布局和農(nóng)田邊界的方法進(jìn)行合理改造，同時，可根據(jù)待改造區(qū)域的宜機(jī)化改造適宜性的差異，合理安排施工順序，從整治難度較大的不適宜，比較不適宜的區(qū)域入手，由點到面建設(shè)為互聯(lián)互通，布局合理的宜機(jī)化耕地。

參 考 文 獻(xiàn)

［1］楊敏麗， 李世武， 惲竹恬. 區(qū)域農(nóng)業(yè)機(jī)械化發(fā)展問題研究［J］. 中國農(nóng)機(jī)化， 2010（1）： 8-13.Yang Minli， Li Shiwu， Hui Zhutian. Study on regional development of agricultural mechanization in China ［J］. Chinese Agricultural Mechanization， 2010（1）： 8-13.

［2］張宗毅. “十四五”期間丘陵山區(qū)農(nóng)田宜機(jī)化改造若干重大問題與舉措［J］. 中國農(nóng)村經(jīng)濟(jì)， 2020（11）： 13-28.Zhang Zongyi. Some important problems and measures of farmland construction suitable for mechanization in hilly and mountainous areas during the 14th five-year plan period ［J］. Chinese Rural Economy， 2020（11）： 13-28.

［3］羅錫文. 對發(fā)展丘陵山區(qū)農(nóng)業(yè)機(jī)械化的思考［J］. 農(nóng)機(jī)科技推廣， 2011（2）： 17-20.Luo Xiwen. Reflections on the development of agricultural mechanization in hilly and mountainous areas ［J］. Agriculture Machinery Technology Extension， 2011（2）： 17-20.

［4］劉小偉. “改地適機(jī)”是丘陵山區(qū)農(nóng)機(jī)化發(fā)展的治本之策［J］. 中國農(nóng)機(jī)化學(xué)報， 2017， 38（7）： 1-2.Liu Xiaowei. Fundamental strategy for developing agricultural mechanization in hilly areas is adapting land for mechanization ［J］. Journal of Chinese Agricultural Mechanization， 2017， 38（7）： 1-2.

［5］陳進(jìn)， 陳敏. 日本、韓國及臺灣地區(qū)丘陵山區(qū)農(nóng)業(yè)機(jī)械化發(fā)展的經(jīng)驗與啟示［C］.中國科學(xué)技術(shù)協(xié)會，重慶市人民政府. 自主創(chuàng)新與持續(xù)增長第十一屆中國科協(xié)年會論文集（2）. 自主創(chuàng)新與持續(xù)增長第十一屆中國科協(xié)年會論文集（2）， 2009： 174-180.

［6］Xiong Yufei. Application of the principles of landscape ecology in land consolidation engineering ［J］. Journal of Physics： Conference Series， 2021， 1865（3）： 1-4.

［7］André Sorensen. Conflict， Consensus or consent： implications of Japanese land readjustment practice or developing countries ［J］. Habitat International， 2000， 24（1）： 51-73.

［8］吳萍， 曹光喬. 韓國農(nóng)田整治對我國丘陵山區(qū)宜機(jī)化改造借鑒與啟示［J］. 中國農(nóng)機(jī)化學(xué)報， 2022， 43（9）： 158-164， 181.Wu Ping， Cao Guangqiao. Reference and enlightenment of farmland construction in South Korea to the farmland mechanized construction in hilly and mountainous areas of China ［J］. Journal of Chinese Agricultural Mechanization， 2022， 43（9）： 158-164， 181.

［9］陳梅， 曹中華. 丘陵山區(qū)宜機(jī)化土地整治項目綜合效益評價研究——基于改進(jìn)AHP—模糊綜合評價法［J］. 中國農(nóng)機(jī)化學(xué)報， 2019， 40（1）： 145-151.Chen Mei， Cao Zhonghua. Research on comprehensive benefits evaluation of land consolidation in hilly and mountainous areas based on method of improved AHP-fuzzy comprehensive evaluation ［J］. Journal of Chinese Agricultural Mechanization， 2019， 40（1）： 145-151.

［10］馬雪瑩， 邵景安， 曹飛. 重慶山區(qū)縣域高標(biāo)準(zhǔn)基本農(nóng)田建設(shè)綜合成效評估——以重慶市墊江縣為例［J］. 自然資源學(xué)報， 2018， 33（12）： 2183-2199.Ma Xueying， Shao Jingan， Cao Fei. Comprehensive performance evaluation of high standard farmland construction in mountainous counties—A case study in Dianjiang， Chongqing ［J］. Journal of Natural Resources， 2018， 33（12）： 2183-2199.

［11］牛坡， 夏飛， 尤顯平， 等. 重慶市宜機(jī)化土地整治效益研究［J］. 農(nóng)業(yè)工程， 2021， 11（8）： 86-89.Niu Po， Xia Fei， You Xianping， et al. Benefits of land consolidation for suitable large-medium farm mechanization of Chongqing ［J］. Agricultural Engineering， 2021， 11（8）： 86-89.

［12］李伶俐， 李江文， 鐘守琴， 等. 基于地形復(fù)雜度的重慶市耕地“宜機(jī)化”改造適宜性評價［J］. 中國生態(tài)農(nóng)業(yè)學(xué)報（中英文）， 2022， 30（2）： 302-313.Li Lingli， Li Jiangwen， Zhong Shouqin， et al. Suitability evaluation of suitable-for-mechanization transformation of cultivated land based on topographic complexity in Chongqing， China ［J］. Chinese Journal of Eco-Agriculture， 2022， 30（2）： 302-313.

［13］劉光盛， 趙樂松， 程迎軒， 等. 基于限制因子的粵北丘陵山區(qū)耕地宜機(jī)化整治分區(qū)［J］. 農(nóng)業(yè)工程學(xué)報， 2021， 37（12）： 262-270.Liu Guangsheng， Zhao Lesong， Cheng Yingxuan， et al. Land consolidation zoning of Northern Guangdong for suitable mechanization transformation in hilly and mountainous areas based on limiting factors ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2021， 37（12）： 262-270.

［14］鐘守琴， 劉涓， 劉衛(wèi)平， 等. 山地丘陵區(qū)耕作田塊修筑工程設(shè)計參數(shù)及田塊特征——以重慶為例［J］. 中國農(nóng)業(yè)科學(xué)， 2017， 50（23）： 4558-4574.Zhong Shouqin， Liu Juan， Liu Weiping， et al. Engineering design parameters of farming plots construction and plot characteristics in hilly area： A case study of Chongqing ［J］. Scientia Agricultura Sinica. 2017， 50（23）： 4558-4574.

［15］胡義萍， 陳超， 曹磊， 等. 改進(jìn)土地整治技術(shù)、推進(jìn)丘陵山區(qū)農(nóng)業(yè)機(jī)械化［J］. 國土資源導(dǎo)刊， 2014， 11（10）：19-22.Hu Yiping， Chen Chao， Cao Lei， et al. Improving land consolidation technology and promoting agricultural mechanization in hilly and mountainous areas ［J］. Land & Resources Herald， 2014， 11（10）： 19-22.

［16］賈興利， 許金良， 楊宏志， 等. 基于GIS的地表破碎指數(shù)計算［J］. 重慶大學(xué)學(xué)報， 2012， 35（11）： 126-130.Jia Xingli， Xu Jingliang， Yang Hongzhi， et al. Calculation of broken index of surface based on GIS ［J］. Journal of Chongqing University， 2012， 35（11）： 126-130.

［17］張正峰， 楊紅， 谷曉坤. 土地整治對平原區(qū)及丘陵區(qū)田塊利用的影響［J］. 農(nóng)業(yè)工程學(xué)報， 2013， 29（3）： 1-8.Zhang Zhengfeng， Yang Hong， Gu Xiaokun. Effects of land consolidation in plains and hills on plots use ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2013， 29（3）： 1-8.

［18］張暉， 王曉峰， 余正軍. 基于ArcGIS的坡面復(fù)雜度因子提取與分析——以黃土高原為例［J］. 華中師范大學(xué)學(xué)報（自然科學(xué)版）， 2009， 43（2）： 323-326.Zhang Hui， Wang Xiaofeng， Yu Zhengjun. Slope surface complexity factor extract and analysis based on ArcGIS ［J］. Journal of Central China Normal University（Natural Sciences）， 2009， 43（2）： 323-326.

［19］李云霞. 基于核K-Means與GIS的縣域地質(zhì)環(huán)境承載力評價——以彭山區(qū)為例［D］. 成都： 成都理工大學(xué)， 2016.Li Yunxia. Evaluation of counties geological environment carrying capacity Based on nuclear K-Means and GIS—A case study Pengshan ［D］. Chengdu： Chengdu University of Technology， 2016.

基金項目：重慶市技術(shù)創(chuàng)新與應(yīng)用發(fā)展專項重點項目（cstc2019jscx-gksbX0118）

第一作者：敖寧玉，女，1997年生，重慶人，碩士研究生；研究方向為丘陵耕地宜機(jī)化適宜性評價。E-mail： aoningyu97@163.com

通訊作者：鮑安紅，女，1969年生，重慶人，博士，教授；研究方向為農(nóng)業(yè)生態(tài)建筑。E-mail： 1193580252@qq.com