基于GIS的秦嶺山區縣域耕地地力調查與評價研究

袁 媛，胡曉黎，雷 蕾，張鴻雁，常慶瑞

（1.商洛市氣象局，陜西 商洛 726000；2. 西北農林科技大學資源環境學院，陜西 楊凌 712100）

【研究意義】耕地作為農業生產最重要的資源之一，是人類生活及社會發展所依賴的物質基礎，是否有效利用土地資源，對區域發展、人民生活水平的高低都有直接影響［1-2］。但是在實際生產過程中，人們常常忽略耕地質量退化、產量變低等不利影響的本質原因，這對我國糧食生產安全和人們生活質量的保障構成嚴重威脅［3-4］。因此，十分有必要建立區域內科學的耕地地力水平評價體系，為土地利用的合理規劃和保持生態平衡提供重要理論依據。【前人研究進展】近年來，眾多國內外學者研究關于耕地利用變化、耕地數量變化、耕地潛力等因素與糧食安全之間的關系，隨著科學技術的進步以及“3S”技術在農業生產中應用地不斷加深，對耕地地力調查與評價研究越來越深入，在該領域研究方法上也不斷更新。【本研究切入點】以陜西省商南縣為例，在全國第二次土壤普查［5-6］基礎上，通過對研究區耕地的地貌特征、土壤結構、養分狀況等的調查，利用“3S”等現代科學技術手段，選取耕地地力評價單元，構建耕地資源基礎數據庫及管理系統，并劃定等級，評價耕地地力水平［7］?！緮M解決的關鍵問題】本研究旨在實現不同等級土地的優化布局，進而指導高效農業生產、培肥土壤，提高研究區耕地資源的可持續發展利用［8］。

1 材料與方法

1.1 研究區概況

商南縣座落在陜西省東南部（110°24' ～111°01'E，33°06' ～ 33°44' N），屬秦嶺東段南麓，3省交匯，地聯8縣，全縣地勢西南和北部高、東南和中部低，海拔落差較大，最高處2 057 m，最低處僅216.4 m。境內有丹江、縣河、滔河、武關河等6條河流，均屬漢江水系?？h域內大部分地區屬暖溫帶氣候，只有東南部為北亞熱帶氣候，歷年平均氣溫13.9 ℃，境內雨量充沛，年降水量862.4 mm，年日照時數1 675.5 h。土壤類型有棕壤、黃棕壤、新積土、潮土和紫色土。全縣土地總面積2 307 km2，主要糧食作物為玉米、小麥、和馬鈴薯。

1.2 數據來源及評價數據庫建立

本研究用于評價的基礎資料包括商南縣土地利用現狀圖（1∶50 000）、土壤圖（1∶50 000）、地形圖（1∶50 000）數據，全國第二次土壤普查的土壤養分數據及相關文字記錄資料，縣域氣候、基本農田及水利設施建設、土壤肥力、土地資源調查和農業生產統計資料等。同時，結合商南縣耕地實際情況，收集野外調查資料，編制《采樣點基本情況調查表》《農戶收入支出情況調查表》《農戶施肥調查表》，調查內容主要包括農田基礎設施、耕作制度、土壤類型、土壤性狀、農民種植業、作物品種、產量及投入產出等方面的情況。根據以上調查內容在全縣耕地上共布設大田采樣點1 444個，采樣時利用GPS定位儀在取樣點確定經緯度，按要求填寫調查表格，并根據實地調查走訪情況，對采樣點位置不合理的情況進行更正。

利用ArcMap 10.2軟件對所用的圖件資料進行矢量化處理，對圖片進行編輯、轉換、空間分析、拓撲分析以后，定性及定量提取商南縣耕地分布、面積等數據，制作輸出耕地資源管理單元圖；對整理的各種調查信息及采樣數據、土肥化驗結果等分類錄入Microsoft Office 2010 Access ，統一將這些數據使用dbf.格式導入ArcGIS 10.2系統中，最終建立商南縣耕地地力調查與評價數據庫［9］。

1.3 研究方法

本研究縣域耕地地力評價技術路線見圖1。

圖1 縣域耕地地力評價技術路線圖Fig. 1 Technical roadmap of county-based cultivated land capacity evaluation

1.3.1 確定評價因子 耕地地力評價是對研究區域土壤理化性狀、地形地貌特征等自然條件對農作物生長影響大小程度的評價。本研究經過實地調查，根據可操作性、穩定性、相關性、主導性等原則，聘請省、市、縣三級農業專家，在全國耕地地力評價指標體系要求下，通過專家組投票方式確定最終指標。選擇12個對耕地地力影響較大且適合當地土地資源實際情況的因素作為評價指標，包括與土壤管理相關的要素（耕地基礎設施具體灌溉情況）、土壤（結構、類型、質地情況）、土壤養分含量（全氮、有效磷、速效鉀、有機質含量）以及地形地貌條件（海拔高度、土地坡度、地貌特征）［10］。

1.3.2 劃分評價單元 為科學劃分研究區耕地地力等級，將商南縣土地利用現狀圖、土壤圖、行政區劃圖在Arc map10.2軟件中疊加形成的色斑圖作為最終評價單元，三類圖形的劃分單元分別為二級利用類型、土種及鎮一級，通過疊加，使得同一評價單元內的行政區域劃分、土種類別及使用方式均達到一致，不同評價單元之間的可比性和差異性則同時體現［11-12］。

1.3.3 計算評價指標權重 不同評價指標對耕地地力評價結果反映土地質量的影響差異較大，各個因子之間相互影響、互相促進，因此應結合當地實際情況對所有選中指標賦予相應的權重［13］。為保證評價的客觀性、準確性、真實性，本研究利用層次分析法進行評價，將所選的12個評價因子按照隸屬關系排成目標層K（耕地地力）、準則層M（有相似特征的要素）、指標層N（參評因子）3個層次［14-15］。專家組在層次結構圖中依據上下兩層對應的重要程度賦予評估值，再經過計算分析形成判斷矩陣，根據矩陣可得到各因子的特征向量（即權重）［16］。各指標的組合權重為MiNj，Mi為M層中決定矩陣的特征向量（i=1，2，3，4）；Nj指N層決定矩陣的特征向量（j=1，2，…12）。經計算，矩陣排列整體一致性較好，層次總排序、單排序及一致性檢驗的一致率（CR）均小于0.1。權重匹配結果見表1。

1.3.4 擬合隸屬函數 根據所選的商南縣耕地地力評價因子，結合模糊數學法，確定3種參評指標與耕地生產關系的函數模型，分別為概念型（耕地基礎設施、灌溉水平、土壤結構、類型、質地、地形地貌）、戒上型（全氮、速效鉀、有機質、有效磷）、直線型（海拔高度和坡度），將后兩種模型的因子進行回歸分析，建立擬合函數模型，通過精度檢驗，置信水平達到顯著的可進行隸屬度計算（表2）。同時，通過經驗法給出評價因素相應的指標值與評估值［17-19］，結果見表3 。最后，

根據不同評價因子與耕地地力等級的相關性，對概念型因子直接賦予隸屬度值［20］（表4）。

表1 評價指標組合權重匹配結果Table 1 Matching results of evaluation index combination weight

表2 所選因子類型及其隸屬函數Table 2 Types of selected factors and the membership function

表3 數量型參評因子專家評估值Table 3 Expert evaluation value of quantitative participation factors

表4 專家組評估非數量型參評因子隸屬度值Table 4 Membership value of non-quantitative type factors evaluated by expert group

1.3.5 耕地地力綜合指數（IFI）的計算 本研究利用加法模型計算IFI［21］，對商南縣耕地地力進行評價，公式如下：

式中，Fi為第i個指標的隸屬度，Ci為第i個指標的組合權重。最后，利用耕地資源管理信息系統中的模塊，計算得出商南縣IFI最小值為0.5261703，最大值為0.862222。

1.4 耕地地力綜合指數分級

在商南縣耕地地力評價中，根據IFI累積頻率曲線圖與實際樣點數，將累積曲線波動劇烈的拐點值作為各個等級的分級指標［22］，將商南縣耕地地力分為5個等級：一級地，IFI＞0.75；二級地，0.7＜IFI≤0.75；三級地，0.65 ＜IFI≤0.7；四級地，0.6＜IFI≤0.65；五級地，IFI≤0.6。

2 結果與分析

基于耕地地力評價手段，利用Arcmap10.2軟件的反距離加權法（IDW）對采樣點的IFI進行空間插值，并根據分級標準對從土地利用現狀圖中提取的耕地進行計算，獲取各評價等級耕地的面積及空間分布情況［23］。

商南縣及各鎮不同耕地地力水平面積統計結果見表5。研究區總耕地面積24 848.26 hm2，占全縣土地面積的10.77%，其中一級地2 286.97 hm2、占總耕地面積的12.19%，二級地4 666.79 hm2、占17.49%，三級地8 733.13 hm2、占32.28%，四級地6 320.22 hm2，占25.10%，五級地2 841.16 hm2、占12.93%。

從評價結果（表5）可以看出，商南縣耕地地力等級與坡度、土壤構造、地貌類型均密切相關，具有明顯的地域分布特征。隨著地力水平的變化，地貌類型從較為平坦的坡塬地和在溝底的溝谷地過渡到高程較高的山地和洪積（或新沖積）形成的灘地，不同類型地貌上各地力等級的耕地面積及占比見表6。

2.1 一級地地力調查評價

從表7可以看出，商南縣一級地面積2 286.97 hm2，占研究區耕地面積的12.19%，全縣各鎮均有分布，主要集中在城關鎮和試馬鎮。地貌類型分布特征是主要集中在河谷川塬上，其余則在淺山丘陵上，其中河谷川塬上的一級地有1 760.76 hm2、占一級地總面積的76.99%，淺山丘陵有526.30 hm2、占24.03%；面積較大的土壤類型為黃棕壤，有1 753.31 hm2、占一級地的76.67%；且一級地坡度較小，大多分布在15°以下，僅有少量地塊坡度略超過15°，無明顯土壤侵蝕現象，土壤養分含量較高。一級地地勢平坦、土層深厚、土壤養分含量高、水土流失較輕，多數地區灌溉條件好，有較高生產力，適宜多種作物的生長。但部分地方易受干旱威脅，在生產上存在肥力后勁不足問題，因此應增施有機肥料，提高土壤抗旱續墑和保肥能力，改善土質結構。在耕種過程中，出現重復開墾、不注重養護等問題，建議合理規劃輪作制度、科學使用配方施肥、改善和加強土壤肥力。

2.2 二級地地力調查評價

商南縣二級地面積4 666.79 hm2，占研究區耕地面積的17.49%，各鎮均有分布，主要集中在城關鎮和富水鎮。其中，分布在淺山丘陵上的面積有4 502.19 hm2、占二級地96.47%，河谷川塬有156.63 hm2、占3.36%，中山地貌中的二級地占比最小、只有7.96 hm2，占0.17%；土壤類型主要為黃棕壤，面積為4 340.03 hm2，占二級地的93.00%，土壤養分含量仍較高（表7）。分布在二級地的地勢也比較平緩，坡度落差不大，土體較厚，水土流失較輕，但部分田塊灌溉條件欠佳，生產性能比一級地低。部分地塊的土種（如泥潮土），質地粘重、通透性較差，在干旱時節常出

現龜裂、雨季積水成澇、耕種易起泥條、限制生產力等問題。因此，需要設法排水除澇，摻沙改善土壤構造，加施有機肥，增強團粒結構，進行土壤培肥與改良等［24］。

表5 商南縣不同地力等級的耕地面積及占比Table 5 Cultivated land area and proportion of different land grades in Shangnan county

表6 商南縣不同類型地貌上各地力等級耕地面積及占比Table 6 Cultivated land area and proportion of different land grades in different geomorphic types in Shangnan county

表7 商南縣不同地力等級耕地土壤養分含量Table 7 Soil nutrient content of different land grades in Shangnan County

2.3 三級地地力調查評價

商南縣三級地面積8 733.13 hm2，占研究區耕地面積的32.28%，各鎮均有分布，主要集中在城關鎮、富水鎮和金絲峽鎮。其中，分布在淺山丘陵上的有三級地 8 619.75 hm2、占98.70%，中山地貌93.36 hm2、占1.07%，河谷川塬僅20.02 hm2、占0.23%；土壤類型主要為黃棕壤，面積為7 724.66 hm2、占88.45%，土壤養分含量略低于前一、二級地（表7）。三級地的坡地耕地面積較前一、二級增多，坡度和水土流失增大，土壤腐殖質層、有效土層的厚度有所降低，且無灌溉能力的耕地面積增加，耕地生產力因此下滑，可通過平整土地、改善土體厚度、增強土壤肥力、提高蓄水保墑能力、防止水土流失、增建灌溉設施等措施，逐步建成高產旱作基本農田。

2.4 四級地地力調查評價

商南縣四級地面積6 320.22 hm2，占研究區耕地面積的25.10%，各鎮均有分布，主要集中在試馬鎮、金絲峽鎮和富水鎮。其中，分布在淺山丘陵上的面積最大、有5 399.22 hm2、占四級地85.43%，中山地貌有917.52 hm2、占14.52%，河谷川塬的分布最少、只有3.48 hm2、占0.05%；土壤類型主要是黃棕壤，面積為4 964.54 hm2、占四級地的78.55%（表7）。四級地的生產性能中等，耕地坡度、水土流失情況繼續增強，肥力減弱，土層變薄，大部分耕地無灌溉能力，嚴重影響生產性能。因此，可進行修建梯田、引水防旱、培肥土壤，提高灌溉能力及土壤濕度，逐步建成穩產旱作基本農田。

2.5 五級地地力調查評價

商南縣五級地面積2 841.16 hm2，僅占研究區耕地面積的12.93%，各鎮均有分布，但主要集中在十里鋪鎮。其中，分布在中山地貌上有2 123.62 hm2、占74.75%，淺山丘陵分布的五級地有717.53 hm2、占25.25%；土壤類型主要是黃棕壤，面積為1 742.37 hm2、占61.33%（表7）。五級地所處的坡度多數較大，土壤表層質地疏松，易導致水土流失，土體厚度、腐殖質層厚度均較低，生產性能中等偏下。對于小于15°的坡耕地，可通過糧草輪作、增施有機肥，提高土壤養分，防止水土流失，建設農田基礎設施，提高耕地生產能力，逐步改良為旱作基本農田。對于大于15°的坡耕地，應全部退耕還林還草。

3 討論

本研究在保護秦嶺生態和發展高效農業、助力脫貧攻堅的大背景下，以地處秦嶺東段南麓的陜西省商南縣作為研究區，建立秦嶺山區縣域耕地地力評價體系，通過對評價指標組合權重匹配結果的分析比較，發現灌溉能力、地貌類型、有機質、坡度為影響商南縣耕地地力評價的主要因子。這些自然條件會限制土層結構、土壤蓄水能力，影響農業生產性能，因此，改善這些因素是提高本區耕地生產地力主要途徑。研究表明，商南縣地力較好的耕地主要集中分布在低海拔（800 m）以下較平坦的川道、丘陵地帶，此區土層深厚，土壤養分含量高，水土流失較輕，多數地區灌溉條件好，是商南縣高產田的主要分布區域，也是科管條件較好，人為干預較多的地區。三、四級耕地多分布在海拔800～1 200 m之間的中高山區，此區海拔較高，坡度增大，水分涵養能力較差，土壤肥力脊薄，交通不便，人為管理難度大的區域。后期可通過加強基礎設施建設，加大科技投入，提高科管水平，更換種植方式（間作、輪作、套種等），改善土體結構，增強土壤肥力，增建灌溉設施等措施提高耕地地力生產水平。五級地所處坡度較大，一般海拔在1 200 m以上、交通極為不便、人為影響較小的區域，此區土壤表層板結，土層薄、腐殖質層厚度均較低，無灌溉條件，管理難度較大，只能采取自然耕作的方式，靠天吃飯。大于15°的坡耕地已不適合耕種，應全部退耕還林還草。分析發現縣域耕地地力固然與海拔、地形地貌等自然因素有關，但可以通人工干預改變或提高地力，隨著現代科學管理水平的提高，耕作方式、種植制度等人為因素對耕地地力也存在一定影響，本研究沒有作深入探討，這為以后研究提出了新的方向。

4 結論

本研究選取土壤管理相關的要素、土體構造、土壤養分含量以及地形地貌條件等4個方面中與耕地地力相關性大且可操作性、穩定性和主導性都較好的12個指標作為參評指標，構建耕地地力優劣評價指標體系，對陜南山區縣域耕地地力進行分級。結果表明，商南縣耕地總面積24 848.26 hm2，1～5級耕地面積分別為 2 286.97、4 666.79、8 733.13、6 320.22、2 841.16 hm2。其中一、二級耕地主要集中分布在海拔800 m以下的河谷川塬和淺山丘陵地帶，三、四級耕地多分布在海拔800～1 200 m之間的淺山丘陵和中山地貌上，而五級地位于海拔1 200 m以上的山坡或高山上。建議通過改善影響地力的主要因子來提高地力，為區域內土地合理利用規劃、農業生產經營、生態平衡保護、耕地資源的可持續發展提供保障。