基于局部空間自相關(guān)的耕地質(zhì)量空間集聚性和保護(hù)分區(qū)

郭敏+李淑杰

摘要：空間自相關(guān)統(tǒng)計量是空間某位置的地理數(shù)據(jù)與其他位置地理數(shù)據(jù)相互關(guān)聯(lián)程度的衡量標(biāo)準(zhǔn)。探討局部空間自相關(guān)在耕地研究中的應(yīng)用，將耕地質(zhì)量等指數(shù)作為空間變量，通過局部空間自相關(guān)，以ArcGIS軟件、Open Geoda軟件為平臺，分析了九臺市耕地的自然質(zhì)量、利用質(zhì)量、經(jīng)濟(jì)質(zhì)量的空間集聚性，提出九臺市耕地保護(hù)分區(qū)的最優(yōu)方案。該方法將耕地自然條件、利用條件、經(jīng)濟(jì)條件進(jìn)行分類考慮，科學(xué)合理并具有一定的可行性，為實現(xiàn)耕地的精細(xì)化管理和差異性保護(hù)提供基礎(chǔ)。

關(guān)鍵詞：耕地質(zhì)量等指數(shù)；局部空間自相關(guān)；耕地保護(hù)分區(qū)；九臺市

中圖分類號： F323.211文獻(xiàn)標(biāo)志碼： A

文章編號：1002-1302（2017）03-0206-05

[HJ1.45mm]

收稿日期：2016-05-09

基金項目：國家自然科學(xué)基金（編號：71303006）；中國科學(xué)院知識創(chuàng)新項目（編號：KZCX2-YW-Q1-07）。

作者簡介：郭敏（1990—），女，內(nèi)蒙古包頭人，碩士研究生，主要從事土地經(jīng)濟(jì)與管理研究。E-mail：597981287@qq.com。

通信作者：李淑杰，博士，教授，主要從事土地經(jīng)濟(jì)與管理研究。E-mail：450127875@qq.com。[HJ]

耕地作為人類生存與發(fā)展不可或缺的提供養(yǎng)育功能的土地資源，其數(shù)量和質(zhì)量備受關(guān)注。保護(hù)耕地的一項重要內(nèi)容就是耕地的質(zhì)量管理。中國先后頒布了《土地管理法》《基本農(nóng)田管理條例》《農(nóng)用地質(zhì)量分等規(guī)程》等法律法規(guī)保護(hù)耕地。耕地的空間分布掌握以及質(zhì)量評定是耕地質(zhì)量管理中的重要措施。耕地質(zhì)量差異除了體現(xiàn)在數(shù)量上，還體現(xiàn)在空間分布上，其空間差異性主要表現(xiàn)為一定的空間聚集性和離散分布性。本研究將耕地的自然等指數(shù)、利用等指數(shù)、經(jīng)濟(jì)等指數(shù)作為空間變量，采用Moran散點圖的方法，結(jié)合局部空間自相關(guān)分析，研究吉林省九臺市耕地質(zhì)量的空間聚集性規(guī)律，提出相應(yīng)的耕地保護(hù)分區(qū)方案，旨在為更好地保護(hù)耕地提供依據(jù)。

1研究區(qū)概況

九臺市位于吉林省長春市東北部，下轄9個鎮(zhèn)、2個民族鄉(xiāng)。九臺市境內(nèi)地勢由西南向東逐漸傾斜，西南高、東北低，年平均氣溫4.7 ℃，年平均降水量577 mm。有效積溫 2 880 ℃，無霜期140～155 d。根據(jù)九臺市2011年土地利用變更調(diào)查成果顯示，該市土地面積33.7萬hm2，其中耕地面積達(dá)21.9萬hm2，占土地面積的65%。

2研究方法

2.1局部空間自相關(guān)

空間自相關(guān)分析是地統(tǒng)計學(xué)中最常用的方法之一，主要是根據(jù)要素位置和要素值來度量空間的相關(guān)性。在給定一組要素及相關(guān)屬性的前提下，通過計算Morans I指數(shù)值、z得分、P值對該指數(shù)的顯著性進(jìn)行評估，其中局部空間自相關(guān)指數(shù)是空間上某一單元與鄰接單元相關(guān)程度的衡量標(biāo)準(zhǔn)。本研究以Local Morans I為統(tǒng)計量進(jìn)行局部空間自相關(guān)分析，以其來反映空間單元i與其鄰域單元之間的空間差異性及顯著性。

局部空間自相關(guān)分析中，主要以Moran散點圖反映空間格局的可視化程度，4個象限分別代表4種關(guān)聯(lián)類型。其中，正相關(guān)類型分別用第1象限的HH型和第3象限的LL型表示；負(fù)相關(guān)類型由第2、4象限代表，分別用HL、LH型表示[1]。

3數(shù)據(jù)來源與處理

3.1數(shù)據(jù)來源

本研究區(qū)域為九臺市，涉及數(shù)據(jù)的時間節(jié)點為2014年，具體資料包括：（1）文本資料，九臺市土地利用現(xiàn)狀、地形圖等成果文件及相關(guān)農(nóng)田水利建設(shè)等統(tǒng)計資料，來源于九臺市國土資源局；（2）數(shù)據(jù)庫資料，2014年九臺市土地利用變更調(diào)查數(shù)據(jù)庫、2013年九臺市耕地質(zhì)量等別評價成果和九臺市第二次土壤普查數(shù)據(jù)庫等，來源于九臺市國土資源局。

3.2數(shù)據(jù)處理

3.2.1空間權(quán)重的確定

為了防止出現(xiàn)比例效應(yīng)，在數(shù)據(jù)處理時要求空間權(quán)重頻率的分布遵循正態(tài)分布，根本目的是為了防止因降低數(shù)據(jù)估計精確程度而導(dǎo)致的某些結(jié)構(gòu)特征無法明顯表達(dá)的問題[2]。采用局部空間自相關(guān)分析研究九臺市耕地質(zhì)量過程中，依據(jù)共邊相鄰（Rook）、共點或共邊相鄰（Queen）、閾值權(quán)重（Threshlod）、指定周圍多邊形個數(shù)（K-nearest）等原則，構(gòu)建空間鄰接矩陣，從空間鄰接性頻率直方圖的分布特征來看，只有共邊相鄰（Rook）符合正態(tài)分布，因此本研究中采用共邊相鄰原則確定空間權(quán)重。

3.2.2耕地質(zhì)量指數(shù)的確定

不同學(xué)者對耕地質(zhì)量內(nèi)涵有不同解釋，本研究采用GB/T 28407—2012《農(nóng)用地質(zhì)量分等規(guī)程》中對耕地質(zhì)量的描述，以自然等、利用等、經(jīng)濟(jì)等作為評價耕地質(zhì)量的依據(jù)，自然等、利用等、經(jīng)濟(jì)等分別對應(yīng)自然質(zhì)量等指數(shù)、利用等指數(shù)、經(jīng)濟(jì)等指數(shù)，自然質(zhì)量等指數(shù)是在自然質(zhì)量分的基礎(chǔ)上乘以光溫生產(chǎn)潛力指數(shù)及指定作物產(chǎn)量比系數(shù)得到，反映耕地自然屬性的優(yōu)劣；在自然質(zhì)量等指數(shù)的基礎(chǔ)上乘以土地利用系數(shù)得到利用等指數(shù)；經(jīng)濟(jì)等指數(shù)是在利用等指數(shù)的基礎(chǔ)上乘以土地經(jīng)濟(jì)系數(shù)得到[3-5]。

如果以耕地質(zhì)量評價數(shù)據(jù)庫最小單元（分等單元）為研究對象，會造成局部空間自相關(guān)分析結(jié)果辨識性差，導(dǎo)致統(tǒng)計結(jié)果不夠準(zhǔn)確?？紤]到九臺市研究區(qū)域范圍較廣，在分析九臺市局部空間自相關(guān)時以行政村為單位。因此，本研究需要將地籍庫中的分等單元進(jìn)行合并，合并后的行政村耕地自然等指數(shù)等于各分等單元耕地自然等指數(shù)按圖斑面積比例求取加權(quán)平均數(shù)。行政村的利用等指數(shù)和經(jīng)濟(jì)等指數(shù)的算法同自然等指數(shù)的算法一致。參照《中國耕地質(zhì)量等級調(diào)查與評定（吉林卷）》，在ArcGIS軟件平臺中[6]，以400分為等別間距確定國家級自然質(zhì)量等，以200分為等別間距確定國家級利用等和經(jīng)濟(jì)等。

4結(jié)果與分析

4.1耕地等別空間特征

2014年九臺市村級耕地質(zhì)量等別空間分布情況如圖1所示。從空間分布上來看，飲馬河和松花江沿岸平原區(qū)附近的耕地自然等別較高，主要耕地類型為水田，具有優(yōu)越的耕地自然條件；九臺市東部、南部、西南部的低山丘陵區(qū)集中分布等別較低的耕地，耕地自然條件相對較差，且耕地類型以旱地為主。整體來看，九臺市耕地的自然質(zhì)量呈東西兩側(cè)高于中部的空間分布特征。其次，利用等別較高的耕地主要分布于飲馬河和松花江沿岸的平原地區(qū)，該地區(qū)擁有較高的耕地生產(chǎn)投入和高水平的田間管理技術(shù)，土地經(jīng)營效益高，產(chǎn)投比高；東南部的低山丘陵區(qū)大量分布有低等別的耕地，該地區(qū)土地生產(chǎn)投入和田間管理技術(shù)均較差；同自然等別一致，從整體來看，利用等別呈東西兩側(cè)高于中間的空間分布特征。經(jīng)濟(jì)等別的空間分布與利用等別的空間分布類似，不再贅述，詳見圖1。[FL）]

[FK（W18][TPGM1.tif；S+3mm][FK）]

4.2自相關(guān)指數(shù)分析

本研究采用Open GeoDa軟件對九臺市耕地質(zhì)量指數(shù)進(jìn)行局部空間自相關(guān)分析，得到圖2所示的Moran散點分布圖。其中，Morans I值反映了以行政村為單元的局部空間自相關(guān)的特性，提示異常值或集聚所在的區(qū)域。Morans I指數(shù)的取值范圍為[-1，1]，即越接近1，表示正相關(guān)性越顯著；而越接近-1，則表示負(fù)相關(guān)性越顯著；越接近0，表示不存在空間聚集性[7-8]。

圖2中，自然等指數(shù)的Morans I為0.66，利用等指數(shù)的Morans I為0.61，經(jīng)濟(jì)等指數(shù)的Morans I為0.57，可以看出九臺市耕地質(zhì)量在空間分布上均表現(xiàn)出較為強(qiáng)烈的正相關(guān)性，表明九臺市耕地質(zhì)量自然條件、利用水平、經(jīng)濟(jì)狀況均具有較大的空間關(guān)聯(lián)性。

4.3局部空間自相關(guān)結(jié)果分析

對村級耕地自然等指數(shù)、利用等指數(shù)、經(jīng)濟(jì)等指數(shù)進(jìn)行局部[CM（25]空間自相關(guān)分析，得到九臺市耕地質(zhì)量指數(shù)局部空間自相[CM）]

關(guān)表（表1）。由表1可知，在95%置信度下，空間上表現(xiàn)為NS（非顯著）型的行政村比例為55%左右。由于本研究以行政村為單元進(jìn)行局部分析，這便將研究區(qū)域的地形地貌、利用管理、人類活動、基礎(chǔ)設(shè)施建設(shè)以及農(nóng)村居民點的空間區(qū)位等因素細(xì)碎化，在本研究單元上存在一定差異性，進(jìn)而導(dǎo)致空間質(zhì)量相關(guān)性不顯著，反向觀之，行政村之間會存在一定的差異性，差異即為合理，若大量行政村之間是相同的，在實際中反而不合理，故即使NS型所占比例為55%，對其進(jìn)行空間分析是合理的。圖3為耕地的空間關(guān)聯(lián)聚集圖，與表1數(shù)據(jù)相對應(yīng)，可以直觀反映出相關(guān)性的空間分布狀態(tài)。自相關(guān)類型共有5種：（1）HH型，表示質(zhì)量等指數(shù)高，且被質(zhì)量等指數(shù)高的行政村包圍的區(qū)域；（2）HL型，表示質(zhì)量等指數(shù)高，但被質(zhì)量等指數(shù)低的行政村包圍的區(qū)域；（3）LH型，表示質(zhì)量等指數(shù)低，但被質(zhì)量等指數(shù)高的行政村包圍的區(qū)域；（4）LL型，表示質(zhì)量等指數(shù)低，且被質(zhì)量等指數(shù)低的行政村包圍的區(qū)域；（5）NS型，代表質(zhì)量等指數(shù)不存在顯著空間相關(guān)性且隨機(jī)分布的區(qū)域[9]。

4.3.1自然等指數(shù)自相關(guān)分析

結(jié)合表1、圖3可以看出，自然等指數(shù)下的局部空間自相關(guān)分析中，正相關(guān)類型中HH、LL型的行政村數(shù)量共占44.25%；然而，負(fù)相關(guān)類型中有4個行政村為HL型，沒有LH型的行政村，負(fù)相關(guān)類型的行政村占行政村總數(shù)的1.21%；NS型行政村有180個。空間分布上，HH型主要分布在興隆鎮(zhèn)東部、龍家堡鎮(zhèn)東部、葦子溝鎮(zhèn)西部、九郊辦事處西部、其塔木鎮(zhèn)東部、莽卡滿族鄉(xiāng)北部；LL型主要分布于上河灣鎮(zhèn)西部、城子街鎮(zhèn)東部交匯處，胡家回族鄉(xiāng)南部，沐石河鎮(zhèn)南部、土們嶺鎮(zhèn)、西營城鎮(zhèn)、波泥河鎮(zhèn)交匯處等；HL型主要零星分布于沐石河鎮(zhèn)北部和葦子溝鎮(zhèn)中南部。

實際情況下，興隆鎮(zhèn)東部、龍家堡鎮(zhèn)東部、葦子溝鎮(zhèn)西部、九郊辦事處西部耕地位于飲馬河沿岸的平原地帶，其塔木鎮(zhèn)東部、莽卡滿族鄉(xiāng)北部位于松花江沿岸的平原區(qū)，地勢平坦，土壤肥沃，具有好的灌排條件和優(yōu)越的自然條件，在空間上與HH型相重疊。耕地自然條件較差的區(qū)域主要集中在東南部低山丘陵地帶，以胡家回族鄉(xiāng)南部以及沐石河鎮(zhèn)南部、土們嶺鎮(zhèn)、西營城鎮(zhèn)、波泥河鎮(zhèn)交匯處為典型代表，上述鄉(xiāng)鎮(zhèn)部分區(qū)域地勢不平，且不具有適宜的有效土壤厚度、優(yōu)良的排水條件和灌溉水源條件，在空間分布上與LL型相匹配。綜上，通過采用局部空間自相關(guān)的方法對九臺市耕地的自然等指數(shù)進(jìn)行分析，得出的結(jié)果與九臺市的實際情況相符，表明該方法科學(xué)合理。

4.3.2利用等指數(shù)自相關(guān)分析

結(jié)合表1、圖3可以看出，利用等指數(shù)下的局部空間自相關(guān)分析中，正相關(guān)類型中HH、LL型的行政村數(shù)量共占行政村數(shù)量的43.94%；負(fù)相關(guān)類型中HL型僅有2個行政村，LH型僅有1個行政村，二者分別占行政村數(shù)量的0.61%、0.30%；NS型行政村有182個。從空間分布角度進(jìn)行分析，正相關(guān)類型多集中連片并且呈“團(tuán)狀”形式出現(xiàn)；而負(fù)相關(guān)類型則呈零星離散狀分布，與正相關(guān)類型具有鮮明的對比特征。HH型主要分布于興隆鎮(zhèn)東部、龍家堡鎮(zhèn)東部、葦子溝鎮(zhèn)西部、九郊辦事處西部（“飲馬河沿岸”組團(tuán)）、其塔木鎮(zhèn)東部、莽卡滿族鄉(xiāng)的北部（“松花江沿岸”組團(tuán)）；LL型主要分布于“上河灣-城子街”組團(tuán)、“沐石河-土門嶺-西營城-波泥河”組團(tuán)， 包括上河灣鎮(zhèn)西部、城子街[FL）]

[FK（W17][TPGM3.tif][FK）]

鎮(zhèn)東部交匯處，胡家回族鄉(xiāng)南部、沐石河鎮(zhèn)南部、土們嶺鎮(zhèn)、西營城鎮(zhèn)、波泥河鎮(zhèn)交匯處等。負(fù)相關(guān)類型（HL、LH型）零星分布在靠近“飲馬河沿岸”組團(tuán)的興隆鎮(zhèn)中部，以及沐石河鎮(zhèn)北部的“上河灣-城子街”組團(tuán)周邊。

從實際上來說，興隆鎮(zhèn)東部、龍家堡鎮(zhèn)東部、葦子溝鎮(zhèn)西部、九郊辦事處西部、其塔木鎮(zhèn)東部、莽卡滿族鄉(xiāng)北部的耕地因多分布于市域河流兩側(cè)，使得耕作灌溉水源豐富且充沛，對水資源的利用比較方便，因此土地利用條件相對較好，空間分布上“飲馬河沿岸”組團(tuán)和“松花江沿岸”組團(tuán)在很大程度上與農(nóng)村居民點相鄰，這不僅縮短了農(nóng)戶與耕作區(qū)域的距離，方便農(nóng)戶耕作，而且對耕作區(qū)的管理也大有裨益，在空間上與HH型分布重疊。而“上河灣-城子街”組團(tuán)、“沐石河-土門嶺-西營城-波泥河”組團(tuán)均處于典型的低山丘陵區(qū)域，地勢不平坦，高低起伏，且可耕作的范圍小，可利用條件差，因此不便于耕作管理，該實際情況與對耕地的利用等指數(shù)進(jìn)行局部空間自相關(guān)分析的LL型耕地聚集區(qū)在空間上十分吻合。

4.3.3經(jīng)濟(jì)等指數(shù)自相關(guān)分析

結(jié)合表1、圖3可以看出，經(jīng)濟(jì)等指數(shù)下的局部空間自相關(guān)分析中，正相關(guān)類型中HH、LL型分別有75、71個行政村，二者之和占行政村數(shù)量的 44.25%，所占比例較大；負(fù)相關(guān)類型中HL型有3個行政村，LH型僅有1個行政村，二者分別占行政村數(shù)量的0.91%、0.30%；NS型行政村有180個。從空間分布角度進(jìn)行分析，正相關(guān)類型多集中連片并且呈“團(tuán)狀”形式出現(xiàn)；而負(fù)相關(guān)類型則呈零星離散狀分布，與正相關(guān)類型具有鮮明的對比特征。HH型主要分布在“飲馬河沿岸”組團(tuán)、“松花江沿岸”組團(tuán)處；LL型主要在“上河灣-城子街”組團(tuán)、“沐石河-土門嶺-西營城-波泥河”組團(tuán)處集中分布。負(fù)相關(guān)類型零星分布在興隆鎮(zhèn)中部，附著于“飲馬河沿岸”組團(tuán)以及沐石河鎮(zhèn)北部、城子街鎮(zhèn)南部依附于“上河灣-城子街”組團(tuán)周邊，無明顯聚集性特征。

在實際中，興隆鎮(zhèn)東部、龍家堡鎮(zhèn)東部、葦子溝鎮(zhèn)西部、九郊辦事處西部、其塔木鎮(zhèn)東部、莽卡滿族鄉(xiāng)北部區(qū)域內(nèi)耕地，市域內(nèi)穿插有多條省道，道路通達(dá)性較高，地勢平坦，利于提升耕地的經(jīng)濟(jì)效益。在空間分布上與HH型耕地聚集區(qū)保持重合。LL型的胡家回族鄉(xiāng)南部、沐石河鎮(zhèn)南部、土們嶺鎮(zhèn)、西營城鎮(zhèn)、波泥河鎮(zhèn)屬典型低山丘陵地帶，地勢起伏不定，不利于交通運輸，且經(jīng)濟(jì)欠發(fā)達(dá)，不利于增加耕地的經(jīng)濟(jì)效益。

4.4耕地保護(hù)分區(qū)的劃分

4.4.1耕地保護(hù)分區(qū)結(jié)果

在耕地保護(hù)實踐方案制定過程中，可從耕地的利用水平、經(jīng)濟(jì)狀況、自然水平3方面，設(shè)計恰當(dāng)?shù)姆謪^(qū)方案。其中，若該區(qū)域耕地自然條件的空間聚集性在利用條件和經(jīng)濟(jì)條件中最顯著，則該分區(qū)方案是該區(qū)域耕地保護(hù)分區(qū)劃分的最優(yōu)方案。本研究以九臺市為實證對象，根據(jù)九臺市實際情況，其自然條件的空間聚集性最顯著（Morans I=0.66），本研究優(yōu)先選擇基于耕地自然質(zhì)量的分區(qū)保護(hù)方案。

耕地自然質(zhì)量不同，其可改良的角度便不同，因此可將耕地保護(hù)劃分為4種類型，分別是重點改良區(qū)、適度改良區(qū)、綜合改良區(qū)、限制建設(shè)區(qū)[10]。對于不同自相關(guān)類型來說，HH型區(qū)域內(nèi)自然質(zhì)量等指數(shù)分值高的耕地十分集中，并且局部范圍內(nèi)呈明顯均質(zhì)狀態(tài)，基于這種情況，該類型區(qū)域應(yīng)該禁止非農(nóng)建設(shè)，且劃入高標(biāo)準(zhǔn)基本農(nóng)田保護(hù)區(qū)。LL型的情況與HH型背道而馳，該區(qū)域?qū)儆诰C合改良區(qū)，要增加有效耕地面積，并且提高耕地整體質(zhì)量，以便經(jīng)濟(jì)社會可持續(xù)發(fā)展和后續(xù)土地利用都需要進(jìn)行整體綜合性整治。HL型區(qū)域內(nèi)高自然等指數(shù)分值的耕地集中分布在被低分值耕地包圍的區(qū)域中部，由于空間極化效應(yīng)對區(qū)域的影響，質(zhì)量好的耕地將慢慢被同化為質(zhì)量較差的耕地，HL型將最終演變成LL型，因此該類型需要進(jìn)行適度保護(hù)。LH型耕地分布情況與HL型恰恰相反，在極化效應(yīng)的作用下演變成HH型，在改良過程中應(yīng)將重點放在高質(zhì)量耕地的改良上，在類型上屬于重點改良區(qū)（表2）。

4.4.2耕地分區(qū)保護(hù)措施

針對九臺市耕地保護(hù)分區(qū)結(jié)果，提出具體保護(hù)措施如下：（1）繼續(xù)維護(hù)HH型對應(yīng)的興隆鎮(zhèn)東部、龍家堡鎮(zhèn)東部、葦子溝鎮(zhèn)西部、九郊辦事處西部、其塔木鎮(zhèn)東部、莽卡滿族鄉(xiāng)北部等耕地自然質(zhì)量高值區(qū)域，以高自然等指數(shù)分值單元為中心，加強(qiáng)其擴(kuò)散效應(yīng)對周邊耕地的影響，限制其非農(nóng)業(yè)建設(shè)；（2）沐石河鎮(zhèn)北部、葦子溝鎮(zhèn)中南部區(qū)域?qū)儆贖L型，為減少周邊低自然等指數(shù)分值耕地對中心高自然等指數(shù)分值的影響，應(yīng)加強(qiáng)對中心區(qū)域耕地的保護(hù)并逐步加大保護(hù)范圍；（3）九臺市上河灣鎮(zhèn)西部、城子街鎮(zhèn)東部交匯處，胡家回族鄉(xiāng)南部，沐石河鎮(zhèn)南部、土們嶺鎮(zhèn)、西營城鎮(zhèn)、波泥河鎮(zhèn)交匯處等區(qū)域?qū)儆贚L型，應(yīng)進(jìn)行適度改良。這些區(qū)域的耕地自然條件普遍較差，需要采取有針對性的改良保護(hù)手段，具體保護(hù)實踐過程中要分主次、逐步推進(jìn)綜合改良措施；同時，LL型區(qū)域也是進(jìn)行非農(nóng)建設(shè)的理想?yún)^(qū)域。根據(jù)上述結(jié)果，結(jié)合空間分布的差異性，遵循因地制宜的原則，科學(xué)合理地制定相應(yīng)的九臺市耕地保護(hù)方式，詳見表2。

5結(jié)論

本研究探索了局部空間自相關(guān)方法在耕地質(zhì)量中的實際應(yīng)用，將耕地自然等指數(shù)、經(jīng)濟(jì)等指數(shù)、利用等指數(shù)作為空間變量，通過應(yīng)用局部空間自相關(guān)方法，以ArcGIS軟件、Open Geoda軟件為平臺，分析了九臺市耕地質(zhì)量的空間集聚性，提出針對九臺市實際情況的耕地保護(hù)分區(qū)的最優(yōu)方案。研究結(jié)果如下：（1）在空間分布上，九臺市耕地質(zhì)量從宏觀上呈現(xiàn)出一定規(guī)律的聚集性特點。具體來說，從耕地自然質(zhì)量來看，九臺市東西兩側(cè)地區(qū)明顯高于中部地區(qū)；從耕地利用狀況來看，市域東部、西部利用質(zhì)量普遍高于東南部的低山丘陵區(qū)。耕地經(jīng)濟(jì)質(zhì)量分布規(guī)律同利用質(zhì)量相同，也是市域東西兩側(cè)地區(qū)高于南部的低山丘陵區(qū)。同時，耕地質(zhì)量在空間分布上呈正相關(guān)類型的區(qū)域集聚性較強(qiáng)，且大部分以“組團(tuán)”形式結(jié)團(tuán)出現(xiàn)；而呈負(fù)相關(guān)類型的區(qū)域多零星分布，無集中連片現(xiàn)象。（2）在Morans I值方面，自然等指數(shù)的Morans I為0.66，利用等指數(shù)的Morans I為0.61，經(jīng)濟(jì)等指數(shù)的Morans I為 0.57，可見九臺市耕地質(zhì)量在空間分布上均表現(xiàn)出較強(qiáng)的正相關(guān)性，反映九臺市耕地質(zhì)量自然條件、利用水平、經(jīng)濟(jì)狀況均具有較大的空間關(guān)聯(lián)性。（3）耕地自然質(zhì)量不同，其可改良的角度便不同。本研究采用局部空間自相關(guān)方法對九臺市耕地空間分布特征及其相關(guān)程度進(jìn)行分析，根據(jù)研究結(jié)果提出了九臺市耕地保護(hù)分區(qū)方案，將耕地保護(hù)劃分為4種類型：重點改良區(qū)、適度改良區(qū)、綜合改良區(qū)、限制建設(shè)區(qū)。針對九臺市耕地實際情況，只有適度改良區(qū)、綜合改良區(qū)、限制建設(shè)區(qū)3大類。對于適度改良區(qū)，對其進(jìn)行保護(hù)的措施集中于積極保護(hù)高值地塊耕地的自然條件，加強(qiáng)對其周邊低值地塊的保護(hù)力度，增加農(nóng)業(yè)投資，防止高值區(qū)向低值區(qū)的退化；對于綜合改良區(qū)，需要實施綜合性的改良措施，同時由于該地區(qū)均為低值區(qū)，耕地較為貧瘠，耕作效率不高，因此該區(qū)域也是農(nóng)地非農(nóng)業(yè)化的理想?yún)^(qū)域；對于限制建設(shè)區(qū)，該區(qū)域耕地條件較好，應(yīng)維持現(xiàn)有水平的耕地條件，該區(qū)可劃為高標(biāo)準(zhǔn)基本農(nóng)田保護(hù)區(qū)，保證糧食產(chǎn)量和安全。

綜上，本研究根據(jù)耕地的自然質(zhì)量、利用質(zhì)量、經(jīng)濟(jì)質(zhì)量的空間集聚類型劃定九臺市耕地保護(hù)分區(qū)。該方法將耕地的自然條件、利用條件、經(jīng)濟(jì)條件分類考慮，科學(xué)合理并具有一定的可行性，為實現(xiàn)耕地的精細(xì)化管理和差異性保護(hù)提供基礎(chǔ)。

[HS2*3]參考文獻(xiàn)：[HJ1.85mm]

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

[2]孫英君，王勁峰，柏延臣. 地統(tǒng)計學(xué)方法進(jìn)展研究[J]. 地球科學(xué)進(jìn)展，2004，19（2）：268-274.

[3]陳印軍，肖碧林，方琳娜，等. 中國耕地質(zhì)量狀況分析[J]. 中國農(nóng)業(yè)科學(xué)，2011，44（17）：3557-3564.

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

[5]胡科，石培基. 甘肅省耕地質(zhì)量評價研究[J]. 中國土地科學(xué)，2008，22（11）：38-43.

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

[7]張松林，張昆. 空間自相關(guān)局部指標(biāo)Moran指數(shù)和G系數(shù)研究[J]. 大地測量與地球動力學(xué)，2007，27（3）：31-34.

[8]張松林，張昆. 局部空間自相關(guān)指標(biāo)對比研究[J]. 統(tǒng)計研究，2007，24（7）：65-67.

[9]Anselin L. Local indicators of spatial association-LISA[J]. Geographical Analysis，1995，27（2）：93-115.

[10]熊昌盛，譚榮，岳文澤. 基于局部空間自相關(guān)的高標(biāo)準(zhǔn)基本農(nóng)田建設(shè)分區(qū)[J]. 農(nóng)業(yè)工程學(xué)報，2015，31（22）：276-284.