基于ESDA的湖南省耕地質量空間差異分析

唐 暉,黃 翔

(1.湖南省不動產登記中心，中國 長沙 410000；2.湖南萬源評估咨詢有限公司，中國 長沙 410000)

基于ESDA的湖南省耕地質量空間差異分析

唐 暉1,黃 翔2*

(1.湖南省不動產登記中心，中國 長沙 410000；2.湖南萬源評估咨詢有限公司，中國 長沙 410000)

利用湖南省2014年度耕地質量等別縣級年度更新數據，采用面積加權平均法，計算各行政區的耕地平均等別；采用探索性空間分析方法(ESDA)分析湖南省耕地質量等別的空間分布情況.研究結果表明：湖南省耕地質量等別在空間分布上呈現較強的聚集性規律，湖南省東北部的等別普遍低于西部及西北部，即東北部的耕地質量要優于西部及西北部.耕地質量等別空間正相關類型多以“抱團”的形式出現，聚集性較強；而負相關類型則無明顯的集中區域，國家經濟等無負相關區域.最后結合耕地質量等別空間差異分析結果，提出耕地質量保護分區建議.

耕地質量等別；空間差異；ESDA方法；湖南省

耕地作為人類必不可少的生存資料，對保證糧食安全、維持社會經濟可持續發展具有重要的現實意義[1].為了更好地保障國家糧食安全，國土資源部在1999—2009年組織全國范圍內的分省農用地分等調查評價任務，并于2009年完成了全國農用地分等工作，為我國耕地質量研究提供了一套完善的成果體系[2]，對我國耕地質量等別差異性有了初步的了解.但耕地質量的差異不僅僅體現在數量等級上，在空間分布上也存在一定的差異，表現出一定的空間集聚性或者離散性[3]，而空間差異是耕地保護工作需要考慮的一個重要因素，鑒于此，本文以湖南省為例，基于2014年度耕地質量等別更新數據，借助GeoDA 1.10和ArcGIS 10.2空間分析軟件，利用ESDA方法深入剖析湖南省各縣、區、市耕地質量國家自然等、利用等、經濟等分布的空間差異性，以期為湖南省耕地保護分區提供科學依據.

1 研究區及數據概況

1.1 研究區概況

圖1 2014年度湖南省耕地質量等別劃分成果Fig.1 Cultivated land quality grade in Hunan province in 2014

湖南省地處中國中部、長江中游，東經108°47′～114°15′，北以濱湖平原與湖北接壤.東北部是幕連九山脈，東部是羅霄山脈，西部是武陵山脈和雪峰山脈，南部是南嶺，中部是湘中盆地群，北部是兩湖盆地里的洞庭湖平原.湖南省共轄13個地級市、1個自治州；35個市轄區、16個縣級市、64個縣、7個自治縣，共122個縣區級行政單元.湖南年日照時數為1 300～1 800 h，年平均溫度為15～18 ℃.年平均降水量為1 200～1 700 mm.根據湖南省2014年度耕地質量更新評價調查結果，湖南省耕地總面積4 153 222.61 hm2，耕地自然質量等分為15個等別，其中1等地自然質量水平最高，15等地自然質量水平最低.從圖1可以看出，湖南省耕地自然質量等別主要集中分布在1～6等，占全部耕地面積的74.06%，整體質量偏好.

1.2 主要數據來源

本文所使用的數據均來自2014年度湖南省耕地質量等別年度更新數據庫以及湖南省行政區劃圖，涉及1 720 722個圖斑.

2 研究方法

2.1 面積加權平均法

面積加權平均法主要是用來計算區域耕地質量的總體平均水平[4].本文采用面積加權平均法計算湖南省各級行政區域的平均自然等別、平均利用等別和平均經濟等別，用以分析該地區的耕地質量狀況，并為進一步進行空間差異分析奠定基礎.計算公式為：

2.2 ESDA方法

2.2.1 空間權重矩陣 通常，在進行空間自相關分析之前，需要定義一個二元對稱空間權重矩陣W來表現空間對象間相互作用的空間關系.ANSELIN[5]認為空間位置的相鄰關系劃分為三類：鄰接、距離和最近K點.采用Rook contiguity標準鄰接規則定義的權重矩陣.空間權重的定義表達式如下：

2.2.2 全局空間自相關 全局空間自相關是對屬性值在整個區域空間特征的描述，反映整個區域的空間集聚程度[6-7].通過全局空間自相關來衡量區域間整體上的空間關聯與差異[8].目前，應用于全局自相關的最主要方法是Moran’sI指數，本文選用該指數來對湖南省耕地質量等別空間模式整體進行定量描述，其計算公式為：

全局空間自相關中Moran指數I的取值范圍為[-1,1]，當I值接近于1時，說明相同屬性的單元集聚在一起(高值與高值相鄰、低值與低值相鄰)；當值接近-1時，相異屬性的單元集聚(高值與低值相鄰、低值與高值相鄰).若Moran’sI指數值接近于0，代表單元屬性隨機分布，單元之間不存在空間相關性[9].

3 結果與分析

3.1 全局空間自相關分析

利用GeoDA 1.10軟件計算出湖南省耕地質量國家自然等別、國家利用等別以及國家經濟等別的全局Moran’sI值分別為0.69,0.84和0.77(見表1)，均趨近于1，所以認為湖南省縣區級耕地質量等別在全局上都呈現較強的正相關性.為了檢驗各年度Moran’sI是否顯著，采用蒙特卡洛模擬方法進行檢驗，結果顯示各項P-value均為0.001，通過顯著性檢驗.這說明湖南省縣區耕地質量等別相似的地區在空間上分布集中(見圖2)，就相關程度來說，耕地國家利用等>國家經濟等>國家自然等.

表1 湖南省各縣區耕地質量等別Global Moran’s I 估計值

(a)國家自然等 (b)國家利用等 (c)國家經濟等圖2 湖南省耕地質量等別Moran散點圖Fig.2 Moran’s scatter diagrams of cultivated land quality grade in scale of county in Hunan province

3.2 局部空間自相關分析

通過計算全局Moran’s I指數，可以知道湖南省的耕地國家自然等、國家利用等和國家經濟等全局自相關關系顯著，高值區和高值區、低值區與低值區出現了明顯的集聚現象.但并不能了解集聚發生的具體位置，為了進一步了解集聚發生的位置，采用局部自相關分析進行探索.運行GeoDA1.10中Univariate Local Moran’s I命令，得到2014年度國家自然等、國家利用等和國家經濟等Cluster Map(見圖3).

表2 湖南省耕地質量等別局部空間自相關類型與行政單元個數統計匯總

注1：P<0.05.高-高(HH) 型(High-High type)；高-低(HL) 型(High-Low type)；低-高(LH) 型(Low-High type)；低-低(LL) 型(Low-Low type).

注2：湖南省包含122個縣區，但常德市鼎城區的行政范圍分為兩個獨立的區域，為了進行空間分析，所以將鼎城區分為兩個圖斑，但相關屬性保持一致.

(a)國家自然等 (b)國家利用等 (c)國家經濟等圖3 湖南省耕地質量等別LISA聚集圖Fig.3 Local spatial agglomeration of cultivated land quality grade in Hunan province

3.2.1 耕地質量國家自然等局部空間自相關分析 根據表2中耕地質量國家自然等局部空間自相關類型統計可知，同類聚集分布的HH型和LL型縣區個數占總行政單元個數的43.90%；在空間上呈現異常分布的HL型和LH型只有3個縣區，僅占湖南省縣區行政單元總數的2.44%；非顯著型呈現隨機分布的縣區有66個，占湖南省縣區行政單元總數的53.66%.結合圖3湖南省耕地質量國家自然等集聚圖，LL型主要“組團”分布在常德市、益陽市、岳陽市、長沙市以及衡陽市部分區域；而HH型主要集中分布在湘西北的湘西、張家界地區以及湘西南的邵陽市.

結合湖南省實際地形地勢分析，LL型“組團”出現在省域東北區域，該地區屬于環洞庭湖平原地區以及湘南南嶺丘陵槽谷地帶，該區域地勢平坦、土壤肥沃，多為良田沃土，歷來就是“魚米之鄉”，這個結果與局部空間關聯圖基本符合；HH型則分布在湘西與湘南的丘陵山地區，易產生水土流失，多開墾為梯田，一部分為坡耕旱地或林地，耕地質量一般較差，這個結果與局部空間關聯圖基本符合.

3.2.2 耕地質量國家利用等局部空間自相關分析 根據表2中耕地質量國家利用等局部空間自相關類型統計可知，同類聚集分布的HH型和LL型縣區個數為70個，占總行政單元個數的56.91%；在空間上呈現異常分布的HL型和LH型只有1個縣區；非顯著型呈現隨機分布的縣區有52個，占湖南省縣區行政單元總數的42.28%.空間分布上，正相關類型呈現鮮明的“組團”分布特征，聚集性較強，這跟全局空間自相關得到的結果一致.HH型主要分布在湘南、湘西南和湘西北地區，LL型主要分布在湘東北地區，包括岳陽、常德、益陽、長沙等地區.

根據湖南省耕地利用實際情況分析，岳陽—常德—長沙—株洲—益陽這一區域的耕地多分布與京廣線以及交通主干道附近，地勢較為平緩，距離居民點較近，能夠得到較好的管理，耕地利用情況好，因此耕地質量國家利用等表現為低值聚集分布；而湘西、懷化、邵陽、永州和郴州大部分地區位于湘西北山原山地區、湘西山地區、湘南山丘區，地形起伏較大，難以進行大規模機械化作業，對耕地的利用較困難，因此耕地質量國家利用等表現為高值聚集分布.

3.2.3 耕地質量國家經濟等局部空間自相關分析 根據表2中耕地質量國家經濟等局部空間自相關類型統計可知，同類聚集分布的HH型和LL型縣區個數為63個，占總行政單元個數的51.22%；在空間上無異常分布的HL型和LH型縣區；非顯著型呈現隨機分布的縣區有70個，占湖南省縣區行政單元總數的48.78%.空間分布上，正相關類型呈現鮮明的“組團”分布特征，聚集性較強.HH型主要分布在湘南、湘西南和湘西北地區，LL型主要分布在湘東北地區，包括岳陽、常德等地區.

根據湖南省耕地投入產出實際情況分析，岳陽、常德和長沙大部分地區處于環洞庭湖平原地區，土壤肥沃、交通條件良好，因此耕作的投入產出比較高，表現為耕地質量國家經濟等別低值聚集分布；而湘西州、懷化、邵陽和永州大部分地區位于湘西北山地區、湘南山丘區，地形起伏大，土壤貧瘠，交通不便，該區域的耕地投入產出比較低，經濟狀況較差，表現為耕地質量國家經濟等別高值聚集分布.

4 結論與討論

采用探索性空間分析方法(ESDA)，以耕地質量等別為變量，探討湖南省耕地質量的空間集聚性規律及空間差異情況，通過研究發現：湖南省耕地質量等別在空間分布上呈現較強的集聚性規律.就耕地質量國家自然等別而言，湖南省東北部的等別低于西部及西北部，也就是說東北部的耕地質量要由于西部及西北部；國家利用等和經濟等在聚集區域在空間分布上較為一致，高等別區域集中在省域西部及南部，而低等別區域集中在東北部.耕地質量等別空間正相關類型多以“抱團”的形式出現，集聚性較強；而負相關類型則無明顯的集中區域，國家經濟等無負相關區域.

湖南省耕地質量等別空間分布及差異分析是進行永久基本農田劃定、耕地保護分區的基礎和前提，根據前文的耕地質量等別空間分布特征及差異分析結果，提出以下耕地質量保護分區建議：

(1) 空間自相關類型呈現LL型區域的耕地應作為永久基本農田劃定的重點區域.該地區耕地質量等別呈現空間低值聚集分布，即高質量耕地集聚區域，該區域耕地應優先劃定為永久基本農田，維持高質量現狀.

(2) 空間自相關類型呈現HH型區域的耕地應作為非農化理想區域.該地區耕地質量等別呈現空間高值聚集分布，即低質量、難利用、產出低下的耕地集聚區域，對于該區域的耕地可以進行綜合性改良，如果需要進行退耕還林應優先考慮這部分耕地.

[1] 李秀彬. 中國近20年來耕地面積的變化及其政策啟示[J]. 自然資源學報,1999,14(4):329-333.

[2] 孔祥斌,張青璞. 中國西部區耕地等別空間分布特征[J]. 農業工程學報,2012,28(22):1-7.

[3] 韋仕川,熊昌盛,欒喬林,等. 基于耕地質量指數局部空間自相關的耕地保護分區[J]. 農業工程學報,2014,30(18):249-256.

[4] 楊建宇,徐 凡,趙 龍,等. 基于綜合質量空間自相關性分析的耕地保護分區方法[J]. 農業機械學報,2017,48(3):156-163.

[5] ANSELIN L．Local indicators of spatial association-LISA[J]．Geogr Anal，1995，27(2):93-115．

[6] 李勝坤,張 毅,閆 欣,等. 武漢城市圈人口空間分布格局及演變分析[J]. 湖南師范大學自然科學學報,2014,37(4):11-15.

[7] 謝 磊,袁華斌,李景保,等. 洞庭湖生態經濟區經濟差異時空格局[J]. 湖南師范大學自然科學學報,2014,37(2):7-11，27.

[8] 蔣國富,劉長運. 河南省縣域經濟的空間分異[J]. 經濟地理,2008,28(4):636-639.

[9] 趙檐瑾. 山西省縣域經濟發展時空格局演變研究[D].西安：陜西師范大學,2012.

[10] 柯新利,黃 翔,楊柏寒. 基于ESDA的武漢城市圈耕地利用集約度空間差異分析[J]. 水土保持研究,2016,23(2):287-291，2.

(編輯 CXM)

Space Difference Analysis of Cultivated Land Quality Grade in Hunan Province Based on ESDA

TANGHui1,HUANGXiang2*

(1. Hunan Real Estate Registration Center, Changsha 410000, China; 2. Hunan Wanyuan Assessment Advisory Co.Ltd.,Changsha 410000, China)

The mean grade in each county-level city was calculated with annual database of 2014 about the cultivated land quality by using the weighted average method. Analysis space difference of cultivated land quality grade in Hunan province based on ESDA. The results showed that the cultivated land quality, which included the natural condition, utilizing situation and economic level, exhibited significant aggregation of global spatial distribution in Hunan province. Besides, the local distribution of cultivated land quality showed that the level of northeast area less than west and northwestern. HH type and LL type emerged as the shape of cluster, while the HL type and LH type were barely sporadic distributed in Hunan province.At last, kinds of protection zoning were put forward according to the results.

cultivated land quality grade; space difference; ESDA; Hunan province

10.7612/j.issn.1000-2537.2017.04.004

2017-07-16

國家社會科學基金資助項目(10CJY027)

X825

A

1000-2537(2017)04-0021-05

*通訊作者,E-mail：hxh66698@sohu.com