周佳松 夏建國 龍丹

摘要 耕地資源是人類生存和發(fā)展的基礎性資源，科學判定耕地空間分布特征，對促進耕地保護及空間優(yōu)化布局等方面具有重要意義。選用耕地質量的國家利用等別，應用空間自相關分析方法來探究漢源縣耕地利用等別空間分布特征。結果表明：①漢源縣耕地質量狀況較好，耕地利用等別主要分布在7～11等，7～8等為高等質量地，占全縣耕地面積的10.27%;9～11等為中等質量地，占全縣耕地面積的89.73%;高等地主要分布在唐家鄉(xiāng)、富春鄉(xiāng)、九襄鎮(zhèn)等鄉(xiāng)鎮(zhèn)，中等地各個鄉(xiāng)鎮(zhèn)均有分布。②全局自相關結果顯示，全局Morans I指數為0.648 947，漢源縣耕地質量有顯著空間集聚性，局部自相關結果顯示漢源縣中部有低—低的低值聚集區(qū)和四周有高—高的高值聚集區(qū)2個顯著相關性區(qū)域，高—低聚集區(qū)和低—高聚集區(qū)零星分布區(qū)域。該研究結果可為調整和補劃基本農田、評價更新耕地質量、管理土地資源及農業(yè)生產提供科學依據。

關鍵詞 耕地利用等別;空間自相關分析;空間分布特征;漢源縣

Abstract Cultivated lands are among the most fundamental resources for the existence and development of mankind. Scientific analysis of the spatial distribution characteristics of cultivated land utilization grade is highly useful in the protection and resource optimization of cultivated lands. Based on the national cultivated land quality grade， the spatial distribution characteristics of cultivated land utilization grade in Hanyuan County was studied by using spatial autocorrelation analysis method. The results showed that ①The cultivated land utilization grade in Hanyuan County was good and mainly distributed from level 7 to 11. Level 7-8 was the highgrade land， accounting for 10.27% of cultivated land area in Hanyuan County， which was mainly distributed in Tangjia， Fuxiang and Jiuxiang towns. Level 9-11 was mediumgrade land， up to 89.63% of the cultivated land area in Hanyuan County and distributed in all villages and towns. ②The global autocorrelation results show that the global Morans I is about 0.648 947， indicating that the quality of cultivated land in Hanyuan County has significant characteristics of spatial agglomeration. The local autocorrelation results show that there are two significant correlation regions， the lowlow low value accumulation area in the middle and the highhigh high value accumulation area in the surrounding， while the highlow and lowhigh accumulation areas are scattered sporadically. The research results provide scientific evidence to adjust and rearrange basic farmlands， to evaluate and update the quality of cultivated land， to manage land resources and promote agricultural production.

Key words Cultivated land utilization grade;Spatial autocorrelation analysis;Spatial distribution characteristics;Hanyuan County

耕地是人類賴以生存的基本資源和重要的農業(yè)資源。耕地質量可以保證國家糧食安全和農產品質量安全[1]。從分布范圍來看，我國的耕地質量總體偏低，呈現優(yōu)等地分布范圍小、數量少，低等地數量多、分布廣泛，各省域內、省際間的耕地等別存在復雜的分異規(guī)律[2-3]。從空間分布上看，耕地質量存在空間聚集性差異性的特征。當前，我國不斷推進新型城鎮(zhèn)化和工業(yè)化進程，調整農業(yè)產業(yè)結構，占用了大量耕地資源[4]。這致使我國耕地資源數量持續(xù)減少，且減少的大部分屬于優(yōu)質耕地，加上我國后備耕地資源嚴重不足，這些因素加劇了我國經濟社會發(fā)展與耕地資源保護的矛盾，嚴重阻礙我國社會經濟的可持續(xù)發(fā)展。通過開展耕地質量等別監(jiān)測與評價工作，我國形成了較齊全的耕地質量數據。對這些數據處理和分析，掌握耕地資源質量空間分布格局及其動態(tài)變化規(guī)律，把握影響空間演變的主要因素，對優(yōu)化耕地質量空間布局與合理利用土地等具有重要意義，同時也為耕地保護和建設提供依據[5-6]。

眾多學者對耕地質量空間分布作出了不同的貢獻，國外學者的研究側重于較為微觀的方面，主要針對耕地的物理指標、化學指標以及生物指標展開研究[7-8]，國內學者主要研究區(qū)域尺度大耕地質量等別分布變化、質量空間分布格局、影響因素等方面[9-10]?？傮w來看，目前多是對全國或省級等較大尺度范圍的耕地質量空間特性分析研究，而缺少對縣級區(qū)域的研究[11-13]。筆者選用耕地質量等別的國家利用等別，應用空間自相關分析方法探究漢源縣耕地質量等別空間分布格局，為漢源縣調整和補劃基本農田、評價更新耕地質量、管理土地資源及農業(yè)生產提供科學依據。

1 研究區(qū)域概況

漢源縣（102°16′～103°01′E、29°05′～29°43′N）轄于四川省雅安市，位于橫斷山脈北段的東部邊緣，屬川西高原與四川盆地的過渡帶。流沙河貫穿南北，大渡河橫貫東西，形成中部低平，四周環(huán)山的河谷區(qū)，海拔550～4 021 m，縣城海拔低，由縣城向四周海拔逐漸遞增，全縣南北長70.1 km，東西寬71.4 km。東至樂山金口河、眉山洪雅縣，西至瀘定縣和雅安市石棉縣，南至涼山彝族自治州甘洛縣，北至雅安市滎經縣?？h城位于富林鎮(zhèn)，距成都市210 km，距雅安市81 km[14]。

漢源縣耕地質量空間分布特征研究范圍是全縣行政轄區(qū)范圍內的所有耕地，包括水田、旱地。具體范圍包括8個鎮(zhèn)、32個鄉(xiāng)轄區(qū)范圍內的耕地，耕地總面積為28 665.60 hm2。

2 數據來源與方法

2.1 數據來源

數據來源于2017年漢源縣耕地質量等別更新評價數據。

2.2 研究方法

使用GIS軟件完成屬性數據的編輯、提取、管理以及制作耕地質量利用等別圖和離散圖;使用Excel軟件完成等別面積量算統(tǒng)計、制作表格。

2.2.1 GIS數據管理、提取和空間分析。

GIS具有數據提取、管理和空間分析的功能，通過提取圖件中的屬性資料，可以導出dbf文件，完成面積的匯總，同時也能為耕地分等提供技術平臺，生成耕地質量等別分布圖、耕地質量等別離散圖，能分析漢源縣耕地質量分布特征和空間相關性[15-17]。

2.2.2 Excel數據分析。

Excel有函數和計算的功能，可打開GIS空間數據庫中輸出的dbf文件，通過應用其函數和計算功能[18]，計算導出的數據，對數據進行匯總和分析，制作耕地面積匯總表。

2.2.3 空間自相關分析。

空間自相關是指同一分布區(qū)域內某些變量之間的相互依賴關系。空間自相關分析能確定某一變量在空間上是否相關及其相關程度。空間自相關度量空間分布需同時根據要素的空間位置和要素屬性值度量。最常用空間自相關的指標和方法的是Morans I統(tǒng)計，其統(tǒng)計量一般可分為局部空間自相關、全局空間自相關2種[12]。

（1）全局自相關。

全局空間自相關是對整個區(qū)域中地理要素屬性值的空間特征描述。全局Morans I指數可分析該區(qū)域的總體空間相關性和空間差異的平均程度[19-20]。故該文選取全局Morans I指數來檢測全局空間自相關程度，其計算模型如下：

式中：n為樣本量，xi、xj分別為空間位置i和j上的觀測值，wij表示空間位置i和j的鄰近關系，當i和j為鄰近的空間位置時wij=1;反之wij=0。全局Morans I的取值介于[-1，1]。標準化Z值常用檢驗全局Morans I指數的顯著性水平，其計算公式如下：

式中：E（I）、VAR（I）分別表示Morans I的期望值和方差。若Z值為正且顯著時，表明存在正的空間自相關，也就是說相似的觀測值（高值或低值）有空間集聚性;若當Z值為負且呈顯著時，則表明存在負的空間自相關，并且類似的觀測值趨向于分散分布;當Z值為零時，不具有相關性，觀測值為獨立隨機分布。

（2）局部自相關。

全局空間自相關是對整個研究區(qū)域的總體描述，但由于局部空間的環(huán)境、經濟、社會等其他因素的不同，其并不能很好的反映，所以就需要探討研究區(qū)的局部空間自相關。

在隨機分布的假設下，用標準化形式Z檢驗局域Morans I系數的顯著性水平。利用GIS軟件可在Z檢驗的基礎上繪制LISA圖[21]。

3 結果與分析

3.1 耕地質量利用等別數量分布特征

按照《農用地質量分等規(guī)程》（GB/T 28407-2012），利用2017年農用地分等成果及其他相關資料對耕地質量等別評價，取得了漢源縣耕地利用等別成果。利用GIS軟件的數據管理、提取功能，提取圖件中的屬性資料，生成2015年全縣耕地國家利用等別圖（圖1）。

全國耕地共分為15等，按照13～15等劃分為低等地，9～12等為中等地、5～8等為高等地、1～4等為優(yōu)等地。漢源縣耕地利用等別區(qū)間為7～11等，為中、高等地。其中以10等地分布最多，占全縣耕地面積的44.52%，其面積為12 706.01 hm2;其次是11等地和9等地，其面積分別為5 911.77、7 049.56 hm2，分別占全縣耕地面積的20.62%、24.59%;7等地和8等地較少，其面積分別占全縣耕地面積的8.02%和2.25%（表1）。

3.2 耕地質量利用等別區(qū)域分布特征

漢源縣耕地質量利用等別有明顯的差異，鄉(xiāng)鎮(zhèn)耕地利用等別為7～11等，屬于中、高等質量耕地（圖1）。漢源縣耕地質量等別總體呈中間低、四周高的分布特征，即漢源中部耕地質量高，四周耕地質量低。

如表2所示，各鄉(xiāng)鎮(zhèn)耕地數量差別較大，九襄鎮(zhèn)、宜東鎮(zhèn)、桂賢鄉(xiāng)、富莊鎮(zhèn)耕地數量較多，均大于1 000 hm2，其中以九襄鎮(zhèn)耕地數量最多，其余鄉(xiāng)鎮(zhèn)面積在200～1 000 hm2。從耕地質量利用等別空間分布上來看，漢源縣耕地質量9～11等分布廣泛，幾乎各個鄉(xiāng)鎮(zhèn)均有分布，只有萬工鄉(xiāng)無10等地和11等地，西溪鄉(xiāng)無11等地。相反7～8等耕地等別分布較為聚集，僅分布在17個鄉(xiāng)鎮(zhèn)。耕地質量相對較差的區(qū)域主要分布在三交鄉(xiāng)、建黎鄉(xiāng)、片馬彝族鄉(xiāng)、前域鄉(xiāng)、后域鄉(xiāng)、富鄉(xiāng)鄉(xiāng)、坭美彝族鄉(xiāng)、永利彝族鄉(xiāng)、順河彝族鄉(xiāng)等鄉(xiāng)鎮(zhèn)，這些鄉(xiāng)鎮(zhèn)耕地利用等別為9～11等，無7等地、8等地分布。其中三交鄉(xiāng)、建黎鄉(xiāng)、坭美彝族鄉(xiāng)、永利彝族鄉(xiāng)、順河彝族鄉(xiāng)耕地質量等別僅為10等、11等，耕地質量相對最差。這些鄉(xiāng)鎮(zhèn)大多位于高寒山區(qū)，大多是高山和坡地，氣候潮濕，水源稀少，耕地數量少，土地貧瘠，土壤質地干硬，自然條件差，土地利用效率和利用程度低。耕地質量較好的區(qū)域主要分布在富莊鎮(zhèn)、大嶺鄉(xiāng)、富泉鄉(xiāng)、河西鄉(xiāng)、西溪鄉(xiāng)等鄉(xiāng)鎮(zhèn)，這些鄉(xiāng)鎮(zhèn)耕地質量等別雖多為9～11等，但含有少量7等、8等耕地分布。耕地質量相對最好的區(qū)域主要分布在九襄鎮(zhèn)、唐家鄉(xiāng)、富春鄉(xiāng)、雙溪鄉(xiāng)、桂賢鄉(xiāng)、大田鄉(xiāng)等鄉(xiāng)鎮(zhèn)，這些區(qū)域耕地質量利用等別為7等、8等分布數量相對最多，其中以九襄鎮(zhèn)的耕地質量利用等別為7等的數量最多，其耕地質量最好。這些區(qū)域內耕地數量多且較為集中，地形坡度較小，臨近河流，水利設施比較完善，離建制鎮(zhèn)和村莊距離較近，交通便利，土地利用效率和利用程度高。

3.3 耕地利用等別空間分布的相關性分析

3.3.1 耕地利用等別的全局自相關分析。

筆者主要運用Arc GIS中的Spatial Statistics模塊計算完成，用Morans I檢驗漢源縣耕地質量是否存在顯著的空間自相關，并使用近似正態(tài)分布假設來檢驗其顯著性?？臻g權重是在局域空間自相關分析的基礎上，用于度量空間單元之間在位置上的依賴程度[22]。該研究選擇距離權重矩陣參與比較計算。首先作出假設：

H0：漢源縣耕地質量不存在空間自相關;

H1：漢源縣耕地質量存在空間自相關。

運用公式（1）（2）對漢源縣的耕地質量數據分析，結果如圖2所示，Morans I指數值為0.648 947，Z得分為415.118 970，表明空間分布上存在顯著的正相關性，拒絕H0，接受H1，即耕地質量存在顯著空間聚集性，即低質量等別耕地與低質量等別耕地相鄰，反之高質量等別耕地與高質量等別耕地相鄰，由此說明漢源縣的耕地質量存在較強的空間關聯性。

3.3.2 耕地利用等別的局部自相關分析。

局部空間自相關是檢驗具有空間位置某要素的觀測值與鄰近地區(qū)該觀測值的空間差異程度，通常用LISA聚類圖來體現，其目的在于分析全局顯著相關的情況下，分析相鄰空間單元之間空間變量的相關程度，局部Morans I值是對全局 Morans I的分解，反映的是各單元的空間集聚程度。通過計算單元Morans I值，并進行顯著性檢驗，生成局部空間自相關LISA聚類分布（圖3）。圖3中黃色是高—高關聯的聚集點，表示該耕地圖斑及周圍圖斑耕地質量利用等別高于全區(qū)該觀測值的均值，具有較強的空間聚集性;藍色代表低—低關聯聚集點，表示耕地質量利用等別低于全區(qū)該觀測值的均值，具有較強的空間聚集性;紅色和玫紅色則是高—低和低—高的關聯點，說明這些耕地圖斑與周圍耕地圖斑有較強的空間異質性，空間聚集性不強;而未通過局部Morans I顯著性檢驗的呈灰色，代表具有低空間集聚性。

由圖3和表3可知，漢源縣52.60%的耕地呈無顯著空間相關性，這個區(qū)域分布較為分散和廣泛，各個地區(qū)均有分布，多分布在漢源縣四周。漢源縣耕地有47.4%具有空間顯著相關性，顯著相關性主要集中在漢源縣中部的鄉(xiāng)鎮(zhèn)。其中高—高聚集區(qū)的比例最大達23.71%，分布較為分散，主要分布在前域鄉(xiāng)、富鄉(xiāng)鄉(xiāng)、安樂鄉(xiāng)、片馬彝族鄉(xiāng)等鄉(xiāng)鎮(zhèn)，有分布在縣域四周的空間分布特征。這些是由研究區(qū)內中等質量的耕地圖斑及其周邊耕地質量平均水平中等的耕地圖斑組成;其次低—低聚集區(qū)比例為20%，主要分布漢源縣中部即在九襄鎮(zhèn)、富春鎮(zhèn)、大田鄉(xiāng)、萬工鄉(xiāng)、桂賢鄉(xiāng)等鄉(xiāng)鎮(zhèn)，大體趨勢是從中部向兩邊遞減，該區(qū)域耕地質量等別多為7等、8等，耕地質量高;高—低和低—高2個區(qū)域的面積總體較小，分別為2.42%和0.51%，主要是零星地分布在低—低、高—高2種集聚地的附近區(qū)域，其空間相關性較小。從表3中可得到各聚類區(qū)等別面積分布情況，該研究區(qū)內高—高集聚區(qū)為10等、11等地，且以11等較多，利用等別較高，屬于中等質量耕地;低—低則主要分布有 7～9等地，且以7 等、9等居多，利用等別較低，質量較好;高—低集聚區(qū)為10等、11等地;低-高集聚區(qū)為7～9 等地，多為9等地。

4 結論

（1）漢源縣耕地質量利用等別分布在7～11等，耕地質量總體為中、高等，其中10.27% 的耕地處在高等地水平，89.73%的耕地處在中等地水平，總體耕地質量較好。

（2）漢源縣耕地質量利用等別空間分布不平衡，受地形地貌、土地利用水平和水資源分布影響較大，總體呈現中間耕地質量高、兩邊耕地質量低的特點，即高等地主要分布在漢源縣中部的富春鎮(zhèn)、大田鄉(xiāng)、九襄鎮(zhèn)、萬工鄉(xiāng)、桂賢鄉(xiāng)等鄉(xiāng)鎮(zhèn);中等地主要分布在四周的片馬彝族鄉(xiāng)、前域鄉(xiāng)、后域鄉(xiāng)、富鄉(xiāng)鄉(xiāng)、坭美彝族鄉(xiāng)、永利彝族鄉(xiāng)、順河彝族鄉(xiāng)等鄉(xiāng)鎮(zhèn)。研究結果與該地區(qū)的實際情況相吻合，可為土地資源管理和農業(yè)生產提供依據。

（3）全局Morans I值為0.648 947，漢源縣耕地質量具有顯著空間集聚性，局部自相關結果顯示有中部低—低低值聚集區(qū)和四周高—高高值聚集區(qū)兩個顯著相關性區(qū)域，高—低和低—高零星分布區(qū)域。這種空間分布特征可為改造中低產田土、評價更新耕地質量、補劃和調整基本農田等方面提供依據。

參考文獻

[1] 沈仁芳，陳美軍，孔祥斌，等.耕地質量的概念和評價與管理對策[J].土壤學報，2012，49（6）：1210-1217.

[2] 王洪波，程鋒，張中帆，等.中國耕地等別分異特性及其對耕地保護的影響[J].農業(yè)工程學報，2011，27（11）：1-8.

[3] 陳藜藜，宋戈，鄒朝暉，等.基于改進局部空間自相關分析的蘿北縣耕地保護分區(qū)[J].資源科學，2016，38（10）：1871-1882.

[4] 譚克均，夏忠敏，韓峰，等.貴州省重點農業(yè)園區(qū)耕地質量現狀及提升對策[J].耕作與栽培，2017（1）：42-44.

[5] 于杰，寧靜，董芳辰，等.1950-2013年三江平原東北部耕地分布變化特征分析[J].干旱區(qū)資源與環(huán)境，2017，31（12）：79-86.

[6] 趙建軍，張洪巖，王野喬，等.基于AHP和GIS的省級耕地質量評價研究：以吉林省為例[J].土壤通報，2012，43（1）：70-75.

[7] GLOVER J D，REGANOLD J P，ANDREWS P K，et al.Systematic method for rating soil quality of conventional，organic，and integrated apple orchards in Washington state[J].Agriculture ecosystems & environment，2000，80（1）：29-45.

[8] ZORNOZA R，MATAIXSOLERA J，GUERRERO C，et al.Evaluation of soil quality using multiple lineal regression based on physical，chemical and biochemical properties[J].Science of the total environment，2007，378（1/2）：233-237.

[9] 孔祥斌，張青璞.中國西部區(qū)耕地等別空間分布特征[J].農業(yè)工程學報，2012，28（22）：1-7.

[10] 宋戈，李丹，梁海鷗，等.松嫩高平原黑土區(qū)耕地質量特征及其空間分異：以黑龍江省巴彥縣為例[J].經濟地理，2012，32（7）：129-134.

[11] 張貞，魏朝富，尚慧.丘陵山區(qū)耕地質量的空間格局分析[J].長江流域資源與環(huán)境，2010，19（8）：901-907.

[12] 張超，張海鋒，楊建宇，等.耕地自然質量空間分布及其影響因子分析：以北京市大興區(qū)為例[J].測繪學報，2015，44（S1）：75-81.

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

[14] 張景華，張建龍，歐陽淵.遙感技術在漢源縣地質災害調查中的應用[J].遙感信息，2010（5）：87-92.

[15] 韓鑫峰.淺析GIS在土地管理中的應用[J].科技風，2014（18）：78.

[16] 楊永俠，孫婷，張麗紅，等.京津冀地區(qū)耕地質量空間分布分形機制研究[J].農業(yè)機械學報，2017，48（2）：165-171.

[17] 施冰臣，文杰，李建華，等.云南省瀾滄江流域耕地自然質量等空間分布特征研究[J].西南農業(yè)學報，2015，28（2）：696-699.

[18] 陳蘭康，段正松，羅崇連，等.田東縣耕地質量空間分布特征研究[J].安徽農業(yè)科學，2014，42（5）：1432-1434.

[19] 任平，吳濤，周介銘.基于GIS和空間自相關模型的耕地空間分布格局及變化特征分析：以成都市龍泉驛區(qū)為例[J].中國生態(tài)農業(yè)學報，2016，24（3）：325-334.

[20] 李黎，王會豪，王玉貴.四川省縣域耕地空間分布與變化研究[J].樂山師范學院學報，2015，30（8）：70-75.

[21] 錢夢茹，梁棟棟，吳旭.基于空間自相關的耕地質量空間分布特征研究：以蕪湖市4區(qū)為例[J].安徽農業(yè)科學，2015，43（35）：315-317，342.

[22] 周巖.基于局域空間自相關的吉林省耕地質量空間分布分析[J].安徽農業(yè)科學，2017，45（32）：200-203，208.