葉云 趙小娟 胡月明

摘要：目前關(guān)于耕地質(zhì)量空間格局方面的研究相對較少，以農(nóng)業(yè)、農(nóng)村經(jīng)濟(jì)發(fā)展最快以及經(jīng)濟(jì)發(fā)達(dá)的華南丘陵平原區(qū)珠三角為研究區(qū)域，在對珠三角耕地質(zhì)量進(jìn)行科學(xué)評價的基礎(chǔ)上，分析其鎮(zhèn)級耕地質(zhì)量空間分布特征，利用局部Morans I指數(shù)，探索其空間關(guān)聯(lián)性。結(jié)果表明：珠三角耕地各類型質(zhì)量的空間分布格局均有所差異。耕地質(zhì)量存在局部空間自相關(guān)，表現(xiàn)為質(zhì)量高和質(zhì)量低的地方有較顯著空間聚集性。各類型質(zhì)量空間聚集類型有較大差異，大部分鄉(xiāng)鎮(zhèn)在空間上表現(xiàn)為非顯著型。從空間關(guān)聯(lián)格局來看，自然質(zhì)量正、負(fù)相關(guān)類型均呈現(xiàn)零散分布;經(jīng)濟(jì)質(zhì)量無負(fù)相關(guān)類型，HH和LL型都是以組團(tuán)形式出現(xiàn)，集中性較強(qiáng);利用質(zhì)量中屬于正相關(guān)類型分布范圍較廣，負(fù)相關(guān)類型多零星分布無明顯集中區(qū)域，局部空間關(guān)聯(lián)格局同利用質(zhì)量基本一致;生態(tài)質(zhì)量與經(jīng)濟(jì)質(zhì)量的空間布局類似。

關(guān)鍵詞：中國珠三角地區(qū);耕地質(zhì)量;空間格局;局部空間自相關(guān);Morans I指數(shù);結(jié)構(gòu)特征;耕地分區(qū)保護(hù);政策建議

確保足夠數(shù)量和一定質(zhì)量的耕地，是保障國家糧食安全、經(jīng)濟(jì)發(fā)展和社會穩(wěn)定的需要，也是生態(tài)環(huán)境保護(hù)的需要[1]。對耕地質(zhì)量空間格局的深入分析，有助于耕地分區(qū)保護(hù)等政策或措施的制定，為實(shí)現(xiàn)國土資源空間優(yōu)化配置提供空間依據(jù)，對其他類似地區(qū)耕地資源保護(hù)及合理利用具有重要的理論和實(shí)踐意義[2]。耕地質(zhì)量研究仍然是當(dāng)今耕地研究的核心和重點(diǎn)之一，國內(nèi)外學(xué)者從不同角度對耕地質(zhì)量內(nèi)涵[3]、評價指標(biāo)體系[4-5]、評價方法[6-7]以及評價結(jié)果運(yùn)用等方面展開了相關(guān)研究。但是關(guān)于耕地質(zhì)量空間格局差異性特征等方面的研究卻相對較少，近年來才相繼有學(xué)者從不同側(cè)重點(diǎn)對農(nóng)用地質(zhì)量的空間格局進(jìn)行了相關(guān)研究[8]。任奎等學(xué)者分別對不同省份農(nóng)用地質(zhì)量的空間分布規(guī)律及其影響因素進(jìn)行了研究[9-12]。在研究內(nèi)容上，袁天鳳等分析了耕地質(zhì)量分布與經(jīng)濟(jì)發(fā)展的空間耦合關(guān)系[13];張貞等對重慶丘陵山區(qū)耕地自然質(zhì)量向耕地經(jīng)濟(jì)質(zhì)量轉(zhuǎn)變的過程及其空間分布格局進(jìn)行了分析[14];宋戈等分析了松嫩高平原黑土區(qū)的耕地質(zhì)量總體特征及空間與行政區(qū)域分布特征[15]。Kraemer等對哈薩克斯坦地區(qū)農(nóng)業(yè)用地覆蓋變化及耕地擴(kuò)張潛力格局進(jìn)行了分析[16]。在研究方法上，學(xué)者們也采用不同的方法對耕地質(zhì)量空間格局進(jìn)行分析，如張超等采用空間自相關(guān)分析方法對北京市大興區(qū)耕地自然質(zhì)量的空間分布特征進(jìn)行研究[17];而趙丹等運(yùn)用趨勢面和相關(guān)分析相結(jié)合的方法研究了耕地質(zhì)量空間分布規(guī)律及影響因素[18];李濤等則基于農(nóng)戶決策行為分析耕地質(zhì)量空間變異性特征規(guī)律[19]?？傮w來看，目前對耕地質(zhì)量空間格局進(jìn)行空間自相關(guān)分析的研究還比較少。珠三角是廣東經(jīng)濟(jì)最發(fā)達(dá)的地區(qū)，也是農(nóng)業(yè)、農(nóng)村經(jīng)濟(jì)發(fā)展最快的地區(qū)，同時作為華南典型丘陵平原區(qū)，具有很好的研究價值和代表性。因此本研究以珠三角為研究區(qū)域，在對珠三角耕地質(zhì)量進(jìn)行科學(xué)評價的基礎(chǔ)上，開展其耕地質(zhì)量空間格局特征及相關(guān)研究，以耕地質(zhì)量的空間差異為切入點(diǎn)，考慮質(zhì)量及空間分布影響，尤其是從空間格局分布上探討其差異及結(jié)構(gòu)特征，有助于耕地分區(qū)保護(hù)等政策或措施的制定。

1 研究區(qū)域概況

珠三角位于廣東省的中南部，珠江下游，地處111°59′～115°26′E，21°27′～23°56′N。行政轄域包括9個地級市，涉及48個縣（區(qū)），土地面積占廣東省國土面積的 23.2%。屬南亞熱帶海洋季風(fēng)氣候，雨熱充沛，溫暖濕潤，年均氣溫為21～23 ℃，多年平均降水量在1 600～2 300 mm之間，界內(nèi)有數(shù)千條大小河道縱橫交錯。三角洲平原土地肥沃，土壤類型多樣，可分為5個土類，包括水稻土、赤紅壤、石灰土、紫色土、潮土等。2014年珠三角土地總面積為547.54萬hm2，其中耕地面積為61.59萬hm2（不包含可調(diào)整地類），占全省耕地總面積的23.49%。

2 研究數(shù)據(jù)及方法

2.1 數(shù)據(jù)來源及處理

研究數(shù)據(jù)主要涉及氣象、水文、土壤、地貌等自然條件統(tǒng)計資料。包括珠三角區(qū)域2014年行政區(qū)劃圖、交通道路，土地利用現(xiàn)狀調(diào)查及規(guī)劃資料等;珠三角及周邊區(qū)域41個氣象站點(diǎn)2014年的年日照時數(shù)、年降水量等數(shù)據(jù)，均來源于廣東省氣象局;土壤屬性數(shù)據(jù)來源于第2次土壤普查資料，包括土壤志、土種志、土壤普查報告等;地質(zhì)災(zāi)害易發(fā)區(qū)相關(guān)數(shù)據(jù)來源于《廣東省地質(zhì)災(zāi)害防治規(guī)劃》（2011—2020年）。影像數(shù)據(jù)主要是在地理空間數(shù)據(jù)云平臺下載的2014年多幅30 m分辨率TM遙感數(shù)據(jù)。社會經(jīng)濟(jì)數(shù)據(jù)來源于2015年《廣東農(nóng)村統(tǒng)計年鑒》、《廣東統(tǒng)計年鑒》及2014年各地市統(tǒng)計年鑒、國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計公報等，所涉及的部門包括農(nóng)業(yè)、國土、水利、規(guī)劃、環(huán)保、統(tǒng)計等部門。

對TM遙感影像數(shù)據(jù)的預(yù)處理主要包括輻射定標(biāo)、大氣校正以及影像的拼接和裁剪等。遙感數(shù)據(jù)分類參考GB/T 21010—2007《土地利用現(xiàn)狀分類》以及中科院土地利用覆蓋分類體系，結(jié)合珠三角區(qū)域TM遙感影像特點(diǎn)和景觀異質(zhì)性，提取出耕地、林地、草地、水域、建設(shè)用地和未利用地。分類結(jié)果的精度檢驗(yàn)是從遙感圖像上隨機(jī)選擇了100個樣點(diǎn)，保證各地類都有覆蓋，采用目視解譯的方法結(jié)合實(shí)地調(diào)查與第2次土地調(diào)查成果，結(jié)果表明影像分類總體精度達(dá)到87.3%，Kappa系數(shù)大于0.7，分類達(dá)到了使用要求。

2.2 研究方法

2.2.1 耕地質(zhì)量評價體系及權(quán)重 劃分珠三角耕地質(zhì)量評價單元采用網(wǎng)格法，考慮到評價區(qū)域的因素差異及面積，以500 m×500 m網(wǎng)格作為評價單元。采用資料統(tǒng)計和專家咨詢法，根據(jù)珠三角區(qū)域特點(diǎn)，從自然質(zhì)量、經(jīng)濟(jì)質(zhì)量、利用質(zhì)量、生態(tài)質(zhì)量等4個方面構(gòu)建適用于該區(qū)域的評價體系，依據(jù)層次分析法結(jié)合專家經(jīng)驗(yàn)判斷指標(biāo)權(quán)重，各指標(biāo)量化分析參考相關(guān)文獻(xiàn)[20-23]，具體見表1。

質(zhì)量評價采用多因素綜合法，根據(jù)各參評因子的分級得分值及其權(quán)重，利用加權(quán)指數(shù)和模型計算耕地質(zhì)量綜合得分值。并采用[0，100]實(shí)現(xiàn)指標(biāo)屬性分值到耕地質(zhì)量綜合評價分值之間的轉(zhuǎn)換[24]，利用ArcGIS柵格計算器將上述各個因素圖層進(jìn)行加權(quán)疊加，得出耕地單元的平均評價分值。

3 結(jié)果與分析

3.1 珠三角鎮(zhèn)級耕地質(zhì)量空間分布特征

參照上述方法實(shí)現(xiàn)耕地質(zhì)量分由耕地單元向鄉(xiāng)鎮(zhèn)尺度的轉(zhuǎn)換。在ArcGIS平臺中，對于鄉(xiāng)鎮(zhèn)耕地質(zhì)量分值參照自然斷點(diǎn)法進(jìn)行分級，將其劃分為相應(yīng)等級，并輸出得到如圖1所示的鄉(xiāng)鎮(zhèn)耕地自然質(zhì)量分值、經(jīng)濟(jì)質(zhì)量分值、利用質(zhì)量分值、生態(tài)質(zhì)量分值的分布圖，分值分為5級，分值越高質(zhì)量越好的為1級，依次遞減，最差的為5級，0表示此區(qū)域沒有耕地。

由圖1可知，在空間分布特征上，鎮(zhèn)級耕地自然質(zhì)量等級的空間分布較零散，珠三角縣域中的鄉(xiāng)鎮(zhèn)均有不同自然質(zhì)量級別的耕地分布，整體來看，中部珠三角平原區(qū)如廣州市部分鄉(xiāng)鎮(zhèn)，以及深圳市、中山市的大部分鄉(xiāng)鎮(zhèn)耕地自然質(zhì)量較好;鎮(zhèn)級耕地經(jīng)濟(jì)質(zhì)量地域差異較明顯，質(zhì)量等級分布比較集中，中部發(fā)達(dá)地區(qū)耕地經(jīng)濟(jì)質(zhì)量相對較高，西北部山地丘陵區(qū)經(jīng)濟(jì)相對落后，其耕地經(jīng)濟(jì)質(zhì)量等級也較低;鎮(zhèn)級耕地利用質(zhì)量表現(xiàn)為珠三角區(qū)域偏東南整體水平較差，其余地區(qū)耕地利用質(zhì)量則相對較高;鎮(zhèn)級耕地生態(tài)質(zhì)量等級分布也較集中，生態(tài)質(zhì)量較差的耕地主要集中在珠三角北部廣寧縣、四會市、鼎湖區(qū)的部分鄉(xiāng)鎮(zhèn)以及南部中山市、南沙區(qū)、新會區(qū)的部分鄉(xiāng)鎮(zhèn)，生態(tài)質(zhì)量較好的耕地主要集中在西南部臺山市以及中部地區(qū)的部分鄉(xiāng)鎮(zhèn)。

3.2 局部LISA分析及Moran散點(diǎn)圖

將珠三角鄉(xiāng)鎮(zhèn)行政單元先進(jìn)行處理，由于城市中心部分街道、鄉(xiāng)鎮(zhèn)沒有耕地數(shù)據(jù)，因此采用合并的方式進(jìn)行處理，例如廣州市荔灣區(qū)的花地街、石圍塘街、茶滘街、多寶街、沙面街、嶺南街、華林街、昌華街、彩虹街、南源街、西村街、逢源街、龍津街、金花街、站前街等，將這些沒有耕地的鎮(zhèn)街合并成一個空間單元，便于后續(xù)分析處理，其他相類似區(qū)域也進(jìn)行相同處理，合并后珠三角仍有廣州市越秀區(qū)、深圳市福田區(qū)和鹽田區(qū)3個縣級行政區(qū)沒有耕地分布。

對以珠三角鄉(xiāng)鎮(zhèn)為空間單元的耕地質(zhì)量進(jìn)行局部空間自相關(guān)分析，采用Rook關(guān)系取一階鄰接方式來構(gòu)造空間鄰接矩陣。在確定空間權(quán)重的基礎(chǔ)上，采用GeoDa軟件進(jìn)行數(shù)據(jù)處理，得到珠三角各鄉(xiāng)鎮(zhèn)耕地質(zhì)量指數(shù)局部空間自相關(guān)性散點(diǎn)圖，如圖2所示，Morans I值即為反映局部空間自相關(guān)性的特征值，均表現(xiàn)出正相關(guān)性，其空間相關(guān)性呈現(xiàn)逐漸增加的特征，須要注意的是，自然質(zhì)量的指數(shù)表現(xiàn)出較弱的空間正相關(guān)性。

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

對鎮(zhèn)級各類耕地質(zhì)量分值的局部自相關(guān)性類型進(jìn)行統(tǒng)計得到表2，表中各空間自相關(guān)類型統(tǒng)計個數(shù)是在P<0.05，即置信度為95%下的統(tǒng)計結(jié)果。與統(tǒng)計結(jié)果相對應(yīng)的局部空間關(guān)聯(lián)聚集圖見圖3，反映具體的空間分布情況。

由表2與圖3-a可知，在95%的置信度下，珠三角大部分鄉(xiāng)鎮(zhèn)耕地自然質(zhì)量在空間上表現(xiàn)為非顯著型。針對以鄉(xiāng)鎮(zhèn)為研究單元的耕地自然質(zhì)量局部空間自相關(guān)類型結(jié)果而言，屬正相關(guān)類型的HH型與LL型各有51、67個鄉(xiāng)鎮(zhèn)，分別占珠三角鄉(xiāng)鎮(zhèn)總數(shù)的10.58%、13.90%，總計占鄉(xiāng)鎮(zhèn)總數(shù)的 24.48%;而屬負(fù)相關(guān)類型的HL型與LH型的鄉(xiāng)鎮(zhèn)分別有7、12個，合計占鄉(xiāng)鎮(zhèn)總數(shù)的3.94%;非顯著型鄉(xiāng)鎮(zhèn)有345個，占鄉(xiāng)鎮(zhèn)總數(shù)的71.58%。相較而言，正相關(guān)類型鄉(xiāng)鎮(zhèn)個數(shù)相對于負(fù)相關(guān)類型占多數(shù)，正、負(fù)相關(guān)類型對比來看，正相關(guān)類型處于優(yōu)勢主導(dǎo)地位。從空間分布格局來看，正、負(fù)相關(guān)類型均呈現(xiàn)零散分布，其中HH型主要分布珠三角中部如白云區(qū)、荔灣

區(qū)的大部分區(qū)域，三水區(qū)中南部，番禺區(qū)北部以及南沙區(qū)北部、中山市中部部分鄉(xiāng)鎮(zhèn)，斗門區(qū)與新會區(qū)相鄰區(qū)域鄉(xiāng)鎮(zhèn)，寶安區(qū)大部分區(qū)域等。從實(shí)際分布來看，這些地區(qū)地處珠三角平原地帶，地勢平坦，珠江、西江、北江、東江等河流貫穿而過，有良好的灌溉排水條件，土質(zhì)較好、土壤肥沃，耕作自然條件優(yōu)越;LL型則主要分布于在廣寧縣、德慶縣的大部分區(qū)域，高明區(qū)西部，鶴山市、恩平市北部，順德區(qū)中部，東莞市東南部以及惠城區(qū)和惠陽區(qū)的部分區(qū)域，與此相對應(yīng)的區(qū)域也多數(shù)處于半山丘陵地帶，地勢較高，排灌條件較差，自然適宜性不高。HL型和LH型多呈現(xiàn)零星分布，其中HL型散布于廣寧縣東南部鄉(xiāng)鎮(zhèn)，東莞市西南部，博羅縣的北部和南部以及惠城區(qū)的部分鄉(xiāng)鎮(zhèn);LH型主要分布于花都區(qū)東部，高要區(qū)北部，博羅縣西部，廣州市中心城區(qū)，深圳市南部以及中山市中部靠北的區(qū)域。

由表2與圖3-b可知，從耕地經(jīng)濟(jì)質(zhì)量局部空間自相關(guān)分析的結(jié)果來看，屬于正相關(guān)類型的HH型和LL型各有121、107個鄉(xiāng)鎮(zhèn)，分別占珠三角鄉(xiāng)鎮(zhèn)總數(shù)的25.10%、22.20%，合計占47.30%，負(fù)相關(guān)類型HL型和LH型均有0個鄉(xiāng)鎮(zhèn);其余則為非顯著型，共有254個鄉(xiāng)鎮(zhèn)，占珠三角鄉(xiāng)鎮(zhèn)總數(shù)的52.70%。再從空間格局來分析，HH型和LL型都是以組團(tuán)的形式出現(xiàn)，集中性較強(qiáng)，其中HH型主要集中在珠三角的中部區(qū)域，廣州-佛山-中山一線，包括廣州市北部大部分區(qū)域和南部的南沙區(qū)，佛山市的南海、順德區(qū)域，西部的端州區(qū)、鼎湖區(qū)、高要區(qū)的部分鄉(xiāng)鎮(zhèn)，蓬江區(qū)、江海區(qū)的北部，以及中山市的大部分地區(qū)。因?yàn)檫@一帶區(qū)域內(nèi)經(jīng)濟(jì)發(fā)展水平相對較高，交通道路四通八達(dá)使得產(chǎn)品相互流通，便于耕地經(jīng)濟(jì)效益的提升，耕地產(chǎn)出率也相應(yīng)較高，靠近城市中心因此輻射帶動作用也強(qiáng);LL型也是以組團(tuán)的形式分布在珠三角區(qū)域的四周，主要包括珠三角西北部的懷集縣、封開縣和廣寧縣的西部鄉(xiāng)鎮(zhèn)，珠三角西南部的恩平市和開平市、臺山市的大部分鄉(xiāng)鎮(zhèn)，珠三角東北部的從化市、龍門縣，東部的惠東縣和惠城區(qū)部分區(qū)域。LL型的組團(tuán)區(qū)域多屬于丘陵山地，灌排條件較差，經(jīng)濟(jì)發(fā)展水平相對落后，交通通達(dá)度相對較低，城市中心輻射帶動作用也較弱。

根據(jù)表2與圖3-c可知，以鄉(xiāng)鎮(zhèn)為研究單元的利用質(zhì)量局部空間自相關(guān)結(jié)果中，屬正相關(guān)類型的HH型和LL型分別有86、79個鄉(xiāng)鎮(zhèn)，合計占珠三角鄉(xiāng)鎮(zhèn)總數(shù)的34.23%，屬于負(fù)相關(guān)類型的HL型和LH型分別有5、9個鄉(xiāng)鎮(zhèn)，兩者僅占鄉(xiāng)鎮(zhèn)總數(shù)的2.91%，其余非顯著型鄉(xiāng)鎮(zhèn)占鄉(xiāng)鎮(zhèn)總數(shù)的62.86%。對比正、負(fù)相關(guān)類型來看，正相關(guān)類型處于絕對優(yōu)勢地位。就其空間格局而言，正相關(guān)類型分布范圍較廣，負(fù)相關(guān)類型多零星分布，無明顯集中區(qū)域，各類型空間分布格局同耕地利用質(zhì)量局部空間關(guān)聯(lián)格局基本一致。正相關(guān)類型中HH型分布于懷集縣東部、南部部分鄉(xiāng)鎮(zhèn)，廣寧縣、德慶縣交界處區(qū)域，珠三角東北部從化市、龍門縣、增城市、博羅縣的部分鄉(xiāng)鎮(zhèn)，以及珠三角西南部臺山市、開平市、鶴山市一帶。這些區(qū)域耕地利用條件較好，耕地連片性及人均耕地面積都相對較高，且多分布在住宅周邊，便于農(nóng)民耕種;LL型主要分布在東莞市，深圳市部分區(qū)域，海珠區(qū)、番禺區(qū)的小范圍鄉(xiāng)鎮(zhèn)，順德區(qū)大部分區(qū)域，南部的新會區(qū)、江海區(qū)相接區(qū)域，以及中山市、香洲區(qū)中部地區(qū)。這些區(qū)域耕地利用條件相對較差，可利用面積較少，且耕地較破碎，不利于大范圍機(jī)械化操作與耕作管理。負(fù)相關(guān)類型中HL型僅離散分布于天河區(qū)、斗門區(qū)和中山市的幾個鄉(xiāng)鎮(zhèn);LH型則零星分布在鼎湖區(qū)西北角，端州區(qū)南部，博羅縣的幾個鄉(xiāng)鎮(zhèn)，以及惠城區(qū)、惠東縣部分區(qū)域，均無明顯的集中特征。

根據(jù)表2與圖3-d可知，以鄉(xiāng)鎮(zhèn)為研究單元的生態(tài)質(zhì)量局部空間自相關(guān)結(jié)果中，屬正相關(guān)類型的HH型和LL型分別有106、83個鄉(xiāng)鎮(zhèn)，合計占珠三角鄉(xiāng)鎮(zhèn)總數(shù)的39.21%，屬于負(fù)相關(guān)類型的只有LH型，僅有3個鄉(xiāng)鎮(zhèn)，占鄉(xiāng)鎮(zhèn)總數(shù)的 0.62%，其余的非顯著型鄉(xiāng)鎮(zhèn)占鄉(xiāng)鎮(zhèn)總數(shù)的60.17%。對比正、負(fù)相關(guān)類型來看，正相關(guān)類型處于絕對主導(dǎo)地位。從空間格局來看，與經(jīng)濟(jì)質(zhì)量的空間布局類似，HH型和LL型也都是以組團(tuán)的形式分片區(qū)集中分布。其中HH型主要分為4個片區(qū)，包括懷集縣北部，高要-高明-禪城-番禺一帶，臺山市、開平市大部分區(qū)域以及增城-博羅-惠城片區(qū)。LL型集中分在廣寧-四會-德慶-高要-鼎湖片區(qū)、蓬江-江海-新會-中山-南沙片區(qū)，以及龍門縣部分區(qū)域。負(fù)相關(guān)類型LH型主要零散分布在博羅縣和惠城區(qū)的幾個鄉(xiāng)鎮(zhèn)?？傮w而言，基于耕地生態(tài)質(zhì)量的局部空間自相關(guān)分析所得到的空間關(guān)聯(lián)聚集圖也與珠三角實(shí)際耕地生態(tài)質(zhì)量情況相符合。

由表2與圖3-e可知，以鄉(xiāng)鎮(zhèn)為研究單元的綜合質(zhì)量局部空間自相關(guān)結(jié)果中，屬正相關(guān)類型的HH型和LL型分別有113、78個鄉(xiāng)鎮(zhèn)，合計占珠三角鄉(xiāng)鎮(zhèn)總數(shù)的39.63%，屬于負(fù)相關(guān)類型的也是只有LH型，僅有3個鄉(xiāng)鎮(zhèn)，占鄉(xiāng)鎮(zhèn)總數(shù)的 0.62%，其余的非顯著型鄉(xiāng)鎮(zhèn)占鄉(xiāng)鎮(zhèn)總數(shù)的59.75%。從正、負(fù)相關(guān)類型來看，正相關(guān)類型也是處于絕對主導(dǎo)地位。整體來看，綜合質(zhì)量的空間集聚分布與經(jīng)濟(jì)質(zhì)量相似，以組團(tuán)的形式分布并且整體格局也有部分區(qū)域一致。HH型主要集中在高要-禪城-南海-花都-從化-白云-番禺-蘿崗-增 城- 博羅一帶;LL型主要以組團(tuán)的形式分片區(qū)分布，包括封開-德慶-廣寧片區(qū)，恩平-臺山片區(qū)，以及龍門縣、惠東縣部分區(qū)域;LH型零星分布在鼎湖區(qū)、博羅縣的幾個鄉(xiāng)鎮(zhèn)?？臻g關(guān)聯(lián)聚集圖也與珠三角實(shí)際耕地綜合質(zhì)量情況相符。

4 結(jié)論與討論

珠三角耕地各類型質(zhì)量的空間分布格局均有所差異。自然質(zhì)量等級分布主要集中在2級和3級，耕地自然質(zhì)量呈中等偏上水平，自然質(zhì)量高等級的耕地主要分布于珠三角中部、東部和南部的平原地區(qū)，而自然質(zhì)量低等級耕地主要分布于北部和西部的丘陵地區(qū);經(jīng)濟(jì)質(zhì)量等級屬較高水平，耕地經(jīng)濟(jì)質(zhì)量等級高的耕地基本集中在中部地區(qū)，低等級地多位于珠三角的外圍區(qū)域，與珠三角在此區(qū)域的經(jīng)濟(jì)發(fā)展程度想吻合;耕地利用等級呈現(xiàn)較高水平，利用質(zhì)量等級高的耕地在整個區(qū)域分布較為廣泛，散布在整個區(qū)域，利用質(zhì)量等級低的耕地多位于珠三角的南部地區(qū);耕地生態(tài)質(zhì)量等級呈現(xiàn)中等偏上水

平，生態(tài)質(zhì)量高等級的耕地主要分布于珠三角中部、東部和西南部的部分區(qū)域，生態(tài)質(zhì)量等級低的耕地多集中于珠三角北部山地丘陵區(qū)。

珠三角耕地質(zhì)量局部均存在空間自相關(guān)，表現(xiàn)為耕地質(zhì)量高的地方和耕地質(zhì)量低的地方有較顯著的空間聚集性。而各類型質(zhì)量空間聚集類型有較大差異，大部分鄉(xiāng)鎮(zhèn)在空間上表現(xiàn)為非顯著型。自然質(zhì)量中屬于正相關(guān)類型的鄉(xiāng)鎮(zhèn)占鄉(xiāng)鎮(zhèn)總數(shù)的24.48%，從空間分布格局來看，正、負(fù)相關(guān)類型均呈現(xiàn)零散分布;經(jīng)濟(jì)質(zhì)量無負(fù)相關(guān)類型，HH型和LL型都是以組團(tuán)的形式出現(xiàn)，集中性較強(qiáng);利用質(zhì)量中屬于正相關(guān)類型的鄉(xiāng)鎮(zhèn)占鄉(xiāng)鎮(zhèn)總數(shù)的34.23%，正相關(guān)類型分布范圍較廣，負(fù)相關(guān)類型多零星分布，無明顯集中區(qū)域，各類型空間分布格局同耕地利用質(zhì)量局部空間關(guān)聯(lián)格局基本一致;生態(tài)質(zhì)量屬正相關(guān)類型的鄉(xiāng)鎮(zhèn)占鄉(xiāng)鎮(zhèn)總數(shù)的39.21%，與經(jīng)濟(jì)質(zhì)量的空間布局類似，HH型和LL型也都是以組團(tuán)的形式分片區(qū)集中分布。

就目前來看，空間自相關(guān)研究已應(yīng)用于許多領(lǐng)域，但運(yùn)用其在耕地中的研究還較少，而分析耕地質(zhì)量在空間上的自相關(guān)性及其結(jié)構(gòu)性特征的研究文獻(xiàn)更為少見，空間自相關(guān)在耕地質(zhì)量空間格局分析中的應(yīng)用還較少，本研究也只是簡單地運(yùn)用自相關(guān)法對耕地質(zhì)量的空間關(guān)聯(lián)性特征進(jìn)行分析，下一步可進(jìn)行深入研究，依據(jù)珠三角耕地不同類型質(zhì)量空間關(guān)聯(lián)格局特征進(jìn)行分區(qū)分類保護(hù)，以制定不同措施確保耕地受到保護(hù)。

參考文獻(xiàn)：

[1]吳次芳，徐保根. 土地生態(tài)學(xué)[M]. 北京：中國大地出版社，2003.

[2]熊昌盛，韋仕川，欒喬林，等. 基于Morans I分析方法的耕地質(zhì)量空間差異研究——以廣東省廣寧縣為例[J]. 資源科學(xué)，2014，36（10）：2066-2074.

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

[4]袁秀杰，趙庚星，朱雪欣. 平原和丘陵區(qū)耕地地力評價及其指標(biāo)體系銜接研究[J]. 農(nóng)業(yè)工程學(xué)報，2008，24（7）：65-71.

[5]Morari F，Lugato E，Giardini L. Olsen Phosphorus，exchangeable cations and salinity in two long-term experiments of north-eastern Italy and assessment of soil quality evolution[J]. Agriculture Ecosystems & Environment，2008，124（1/2）：85-96.

[6]Masto R E，Chhonkar P K，Singh D，et al. Alternative soil quality indices for evaluating the effect of intensive cropping，fertilisation and manuring for 31 years in the semi-arid soils of India[J]. Environmental Monitoring and Assessment，2008，136（1/2/3）：419-435.

[7]郄瑞卿，關(guān) 俠，鄢旭久，等. 基于自組織神經(jīng)網(wǎng)絡(luò)的耕地自然質(zhì)量評價方法及其應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報，2014，30（23）：298-305.

[8]張桂花，于 婧，羅清泉，等. 武漢市遠(yuǎn)城區(qū)耕地質(zhì)量空間格局分析[J]. 農(nóng)業(yè)與技術(shù)，2014（2）：221-224.

[9]任 奎，周生路，張紅富，等. 江蘇農(nóng)用地資源質(zhì)量空間格局及影響機(jī)制研究[J]. 農(nóng)業(yè)工程學(xué)報，2008，24（4）：127-134，封3.

[10]張紅富，周生路，吳紹華，等. 江蘇省農(nóng)用地質(zhì)量空間格局及其影響因素分析[J]. 資源科學(xué)，2008，30（2）：221-227.

[11]賈樹海，邱志偉，潘錦華. 遼寧省農(nóng)用地質(zhì)量空間分布格局及影響因素研究[J]. 生態(tài)環(huán)境學(xué)報，2010，19（5）：1143-1150.

[12]張曉燕，張 利，陳 影，等. 河北省農(nóng)用地質(zhì)量空間格局的計量地理分析[J]. 水土保持研究，2010，17（1）：101-106.

[13]袁天鳳，邱道持. 重慶市耕地質(zhì)量分布與經(jīng)濟(jì)發(fā)展的耦合[J]. 經(jīng)濟(jì)地理，2008，28（3）：133-146.

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

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

[16]Kraemer R，Prishchepov A V，Mueller D A，et al. Long-term agricultural land-cover change and potential for cropland expansion in the former Virgin Lands area of Kazakhstan[J]. Environmental Research Letters，2015，10（5）：054012.

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

[18]趙 丹，趙華甫，饒 杰，等. 基于趨勢面的耕地質(zhì)量空間分異特征及影響因素[J]. 水土保持研究，2015，22（6）：219-223.

[19]李 濤，孔祥斌，梁 穎，等. 基于農(nóng)戶決策行為的耕地質(zhì)量評價理論與方法構(gòu)建[J]. 中國農(nóng)業(yè)大學(xué)學(xué)報，2010，15（3）：101-107.

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

[21]奉 婷，張鳳榮，李 燦，等. 基于耕地質(zhì)量綜合評價的縣域基本農(nóng)田空間布局[J]. 農(nóng)業(yè)工程學(xué)報，2014，30（1）：200-210.

[22]路 婕，李 玲，吳克寧，等. 基于農(nóng)用地分等和土壤環(huán)境質(zhì)量評價的耕地綜合質(zhì)量評價[J]. 農(nóng)業(yè)工程學(xué)報，2011，27（2）：323-329.

[23]陳 朝，呂昌河. 基于綜合指數(shù)的湖北省耕地質(zhì)量變化分析[J]. 自然資源學(xué)報，2010，25（12）：2018-2029.

[24]趙小娟，葉 云，周晉皓，等. 珠三角丘陵區(qū)耕地質(zhì)量綜合評價及指標(biāo)權(quán)重敏感性分析[J]. 農(nóng)業(yè)工程學(xué)報，2017，33（8）：226-235.

[25]王 磊，李 濤，曹小曙. 基于ESDA-GIS的廣東省城鄉(xiāng)統(tǒng)籌發(fā)展空間分異[J]. 經(jīng)濟(jì)地理，2012，32（9）：44-50.

[26]王 永，沈 毅. 空間自相關(guān)方法及其主要應(yīng)用現(xiàn)狀[J]. 中國衛(wèi)生統(tǒng)計，2008，25（4）：443-445.

[27]謝花林，李秀彬，張燕婷，等. 基于ESDA的京津冀地區(qū)草地變化空間分異[J]. 自然資源學(xué)報，2012，27（7）：1224-1232.張桀滈，李江風(fēng)，譚 旭，等. 基于ArcScene的丘陵地區(qū)土地整治三維可視化方法[J]. 江蘇農(nóng)業(yè)科學(xué)，2019，47（6）：203-207.