環(huán)洱海地區(qū)耕地與農(nóng)業(yè)勞動力空間分布特征及其耦合關(guān)系

張磊 武友德 李君

摘 要：運用區(qū)域重心、地理集中度及不一致指數(shù)等模型，對環(huán)洱海地區(qū)1988年～2012年耕地與農(nóng)業(yè)勞動力的空間分布與耦合關(guān)系進行分析。得出以下結(jié)論：隨時間演進，耕地重心呈向西北為主、東南為輔的變動態(tài)勢，而農(nóng)業(yè)勞動力重心的階段性變動明顯，整體向西南偏移；25年間，耕地與農(nóng)業(yè)勞動力重心的空間重疊性與變動一致性均大幅提高，兩重心耦合態(tài)勢日趨明顯。海西地區(qū)耕地與農(nóng)業(yè)勞動力地理集中度均遠高于海東，洱海西岸、S221及G214沿線分別為耕地與農(nóng)業(yè)勞動力的主次要集聚軸；自1988年以來，耕地向海西中北部集聚，而農(nóng)業(yè)勞動力呈現(xiàn)出由分散轉(zhuǎn)為集中的布局態(tài)勢，整體上兩者呈高度正向相關(guān)性；隨時間演進，村域耕地與農(nóng)業(yè)勞動力空間分布的協(xié)調(diào)度明顯下降，依據(jù)不一致指數(shù)，可分為耕地優(yōu)勢集聚型、農(nóng)業(yè)勞動力分布滯后型、耕地—農(nóng)業(yè)勞動力空間協(xié)調(diào)型、耕地分布滯后型及農(nóng)業(yè)勞動力優(yōu)勢集聚型五類村莊。

關(guān)鍵詞：耕地；農(nóng)業(yè)勞動力；空間耦合；環(huán)洱海地區(qū)

中圖分類號：F301.24 文獻標識碼：A

0 引言

在城鎮(zhèn)化與工業(yè)化快速發(fā)展的背景下，區(qū)域耕地快速流動的同時，農(nóng)業(yè)勞動力能否隨之非農(nóng)轉(zhuǎn)移，促進勞耕比例的適應性調(diào)節(jié)，已成為保障鄉(xiāng)村系統(tǒng)有序發(fā)展的重要前提[1]。作為反映區(qū)域農(nóng)業(yè)發(fā)展狀況的重要指標，耕地與農(nóng)業(yè)勞動力已成為鄉(xiāng)村地域人地關(guān)系研究中的重要內(nèi)容。

近年來，國內(nèi)學者對耕地與農(nóng)業(yè)勞動力空間分布及變化的研究進展明顯。高凱等以高分辨率土地覆被影像為基礎，運用Ripley K函數(shù)[2]、 Logistic 回歸模型[3]及ESDA[4]等數(shù)學與空間統(tǒng)計分析方法，對耕地時空變化的特征[5，6]、格局[7，8]、景觀形態(tài)[2]、驅(qū)動機制[9，10]與未來預測[11，12]等方面進行了探討；蔡昉、李小建等對農(nóng)業(yè)勞動力特別是農(nóng)村剩余勞動力轉(zhuǎn)移或供給的理論政策[13，14]、影響因素[15-17]及時空格局[18，19]等方面進行了分析；姜曉清等分別基于勞耕彈性系數(shù)、超越對數(shù)隨機前沿生產(chǎn)函數(shù)等模型對耕地與農(nóng)業(yè)勞動力的時空耦合[20]、農(nóng)業(yè)勞動力的流動[21]或年齡結(jié)構(gòu)[22]對耕地處置方式或利用效率的影響等方面進行了探討；于慧[23]、馬彩虹[24]及湯進華[25]等分別從省域視角出發(fā)，對陜西、湖北兩省耕地與農(nóng)業(yè)勞動力的耦合與協(xié)調(diào)進行了探討?？傮w來看，現(xiàn)有研究多側(cè)重于從國家、區(qū)域、省市的宏觀及中觀層面對耕地或農(nóng)業(yè)勞動力分布及變化的單要素分析，而從縣級及以下的微觀視角出發(fā)，對耕地與農(nóng)業(yè)勞動力空間分布的關(guān)聯(lián)性及耦合關(guān)系的研究較為欠缺。

環(huán)洱海地區(qū)是典型的高原湖泊平壩區(qū)，人地關(guān)系特殊。一方面，隨著區(qū)域內(nèi)非農(nóng)化進程的加速，耕地面積大幅變動，耕地“紅線”保護嚴格，土地空間組織與形式發(fā)生顯著改變；另一方面，旅游業(yè)的迅速發(fā)展及高度商業(yè)化的經(jīng)濟使該區(qū)農(nóng)業(yè)勞動力流動頻繁，呈獨特的“鎮(zhèn)外縣內(nèi)”流動態(tài)勢。在此背景下，探討環(huán)洱海地區(qū)耕地與農(nóng)業(yè)勞動力的空間耦合狀況，對合理評價該區(qū)勞耕平衡比例，促進區(qū)域人地關(guān)系穩(wěn)定，有重要的實踐意義。

1 研究區(qū)域與數(shù)據(jù)來源

環(huán)洱海地區(qū)位于云南西部，是我國白族的重要聚居地，也是云南省同時具有高原湖泊生態(tài)脆弱區(qū)、民族文化多元融合區(qū)、鄉(xiāng)村經(jīng)濟發(fā)展活躍區(qū)等多重疊合特征的典型區(qū)域[26]。隨云南省“城鎮(zhèn)上山”項目的啟動和大理市海東新區(qū)建設的推進，該區(qū)旅游業(yè)和市區(qū)城鎮(zhèn)化快速發(fā)展，農(nóng)村居民的流動性明顯增強，生產(chǎn)方式日趨多樣，新區(qū)修建、舊村改造使耕地與農(nóng)業(yè)勞動力的空間“脫鉤”分布日趨明顯。

本文數(shù)據(jù)主要由經(jīng)緯度空間數(shù)據(jù)和區(qū)域發(fā)展屬性數(shù)據(jù)構(gòu)成，為保障數(shù)據(jù)連續(xù)性與精確性，選取環(huán)洱海地區(qū)96個行政村為研究對象?；趯嵉卣{(diào)研分析，采用GPS地理定位法，提取各行政村幾何重心的經(jīng)緯度，并基于ArcGIS10.1構(gòu)建行政村坐標體系網(wǎng)作為空間數(shù)據(jù)；取環(huán)洱海地區(qū)各行政村1988年以來的年末集體所有耕地總面積和年末匯總勞動力數(shù)作為屬性數(shù)據(jù)，屬性數(shù)據(jù)主要取自《大理市農(nóng)村經(jīng)濟收支情況統(tǒng)計年報》（1988-2013）、《大理市統(tǒng)計年鑒》（1988-2013）、《洱源縣國民經(jīng)濟和社會發(fā)展統(tǒng)計年鑒》（1988-2013）及課題組實地調(diào)研資料①。

2 研究方法與理論模型

2.1 重心測度模型

區(qū)域重心是衡量某屬性在區(qū)域總體分布狀況的重要指標，可表示其分布的總趨勢和中心區(qū)位[27]。本文引入重心測度模型從整體上測度耕地與農(nóng)業(yè)勞動力的空間分布及變動狀況。假設某區(qū)域由個平面單元構(gòu)成，其中第個平面單元的重心坐標為（xi，yi），Mi是該單元的屬性值，則這一區(qū)域的重心坐標為[28]：

3 耕地與農(nóng)業(yè)勞動力空間重心耦合分析

以環(huán)洱海地區(qū)96個行政村年末集體所有耕地總面積和年末匯總農(nóng)業(yè)勞動力數(shù)為指標，以各村幾何重心地理坐標為著力點，基于公式（1）測算該區(qū)1988年～2012年的耕地與農(nóng)業(yè)勞動力的空間分布重心坐標和變動軌跡（表1、圖1）。

3.1 耕地空間分布重心變動分析

環(huán)洱海地區(qū)耕地重心位于幾何重心（100.2185E，25.7388N）的西北方向，平均偏移約4570.3米。隨時間演進，耕地重心呈兩個明顯的階段性變動特征：1988年～2006年，耕地重心整體向西北傾斜，海西②北部的喜洲及灣橋兩鎮(zhèn)成為耕地開發(fā)重點，洱海東西兩側(cè)耕地分布差距擴大。2006年～2012年，耕地重心轉(zhuǎn)向東南方向，受“海東開發(fā)”政策驅(qū)動，開發(fā)條件相對優(yōu)越的海東③南部鳳儀及海東兩鎮(zhèn)逐步成為耕地開發(fā)重點，洱海東西兩側(cè)耕地分布差距有所減小。

3.2 農(nóng)業(yè)勞動力空間分布重心變動分析

環(huán)洱海地區(qū)農(nóng)業(yè)勞動力重心位于幾何重心西北方向，平均偏移約4421.2米。隨時間演進，農(nóng)業(yè)勞動力分布變動的區(qū)域差異明顯，重心的階段性變動比耕地更復雜，呈現(xiàn)出四個發(fā)展階段：1）1988年～1995年，農(nóng)業(yè)勞動力重心向西南方向移動超過770米，表明此階段下受主城區(qū)吸引，農(nóng)業(yè)勞動力向該區(qū)西南部集聚，區(qū)域間農(nóng)業(yè)勞動力分布差異擴大。2）1995年～2000年，農(nóng)業(yè)勞動力重心向東南方向移動1110米，表明此階段受大理市經(jīng)濟技術(shù)開發(fā)區(qū)驅(qū)動，農(nóng)業(yè)勞動力向東南部集聚態(tài)勢明顯，區(qū)域間農(nóng)業(yè)勞動力分布差異有所減小。3）2000年～2006年，農(nóng)業(yè)勞動力重心向西北方向移動超過760米，表明此階段海西中北部的大理及喜洲等鎮(zhèn)旅游資源的開發(fā)，特別是大理古城景區(qū)的發(fā)展與大理大學的新擴促使農(nóng)業(yè)勞動力向西北部流動，區(qū)域間農(nóng)業(yè)勞動力分布差異再次擴大。4）2006年～2012年，農(nóng)業(yè)勞動力僅向正南偏西方向移動223.3米，表明在此期間海西南部下關(guān)主城區(qū)的收入、教育及社保等優(yōu)勢再次成為吸引農(nóng)業(yè)勞動力的重要驅(qū)動力，農(nóng)業(yè)勞動力再次呈現(xiàn)出向西南部集聚的態(tài)勢。

3.3 耕地與農(nóng)業(yè)勞動力重心空間耦合特征

基于公式（2）、（3），通過測度兩重心的空間重疊性及變動一致性，對兩重心的空間耦合態(tài)勢進行分析。由表2發(fā)現(xiàn)，1988年～2012年，環(huán)洱海地區(qū)耕地與農(nóng)業(yè)勞動力重心的空間耦合呈以下特征：

1）耕地與農(nóng)業(yè)勞動力重心的空間重疊性經(jīng)歷了“升高～降低～升高”的三個變動階段，整體明顯提高。1988年～1995年，兩重心間距離由1108m降至373.7m，空間重疊性顯著提高，兩重心在分布上趨于集聚；1995年～2006年，兩重心間距離由373.7m升至779.9m，空間重疊性有所降低，在分布上趨于離散；2006年～2012年，兩重心間距離再次縮小，由779.9m降至415.4m，兩重心在分布上再次趨于集聚。

2）耕地與農(nóng)業(yè)勞動力重心的變動一致性指數(shù)整體呈上升態(tài)勢，兩重心間的變動一致性持續(xù)增強。1988～1995年，兩重心間的變動一致性指數(shù)小于0，且絕對值較低，兩重心變動方向與軌跡的偏差較大，空間變動的一致性較弱；1995～2000年，兩重心間的變動一致性指數(shù)大幅提升，由前期的-0.258增至0.299，兩重心變動方向與軌跡的偏差逐步縮小，并開始呈現(xiàn)出同向變動態(tài)勢，空間變動的一致性明顯增強；2000～2006年，兩重心間的變動一致性指數(shù)持續(xù)提升，兩重心變動方向與軌跡的同向態(tài)勢更加明顯，空間變動的一致性進一步增強；同前期相比，2006～2012年兩重心間的變動一致性指數(shù)雖小幅減小，兩重心變動方向與軌跡仍保持明顯的同向態(tài)勢，空間變動的一致性仍較強。

3）隨時間演進，環(huán)洱海地區(qū)耕地與農(nóng)業(yè)勞動力重心的空間重疊性與變動一致性均明顯提升，重心耦合態(tài)勢日趨明顯，兩重心空間分布的偏離度明顯減小。

4 耕地與農(nóng)業(yè)勞動力空間分布耦合分析

重心分析可整體測度耕地與農(nóng)業(yè)勞動力的空間分布及變動特征，但對于區(qū)域內(nèi)部各行政村間的空間分布狀況測度不足。鑒于此，本文引入地理集中度模型，依據(jù)公式（4），分別測算環(huán)洱海地區(qū)96個行政村的耕地與農(nóng)業(yè)勞動力地理集中指數(shù)，并基于ArcGIS10.1的自然間斷分級法，將地理集中指數(shù)由高到低依次分成五個等級（圖2、圖3）。

4.1 耕地地理集中度分布格局

環(huán)洱海地區(qū)耕地地理集中度呈明顯的東西分異，海西地區(qū)整體遠高于海東，耕地主要沿洱海西岸及S221集中分布，耕地的次要集聚軸為G214沿線。

從集中指數(shù)看，耕地地理分布的區(qū)域差異明顯。海西中部的大理、銀橋及灣橋三鎮(zhèn)主要由第一、二級村莊構(gòu)成，為耕地高度集中分布區(qū)；海西北部的喜洲及上關(guān)兩鎮(zhèn)第二、三級村莊所占比重較大，為耕地相對集中分布區(qū)；海西南部的太邑鄉(xiāng)及海東的鳳儀、海東、挖色及雙廊四鎮(zhèn)90%以上為第四、五級村莊，為耕地稀疏分布區(qū)；需要指出的是，海西南部下關(guān)鎮(zhèn)各村間耕地地理集中指數(shù)差異明顯，以G320為界，北部村莊遠高于南部。

從集中區(qū)域看，隨時間推移，耕地向海西中北部地區(qū)集聚態(tài)勢明顯。受大麗線開通及海西開發(fā)重心北移的驅(qū)動，海西中北部的銀橋、灣橋及喜洲三鎮(zhèn)耕地地理集中指數(shù)占環(huán)洱海地區(qū)的比重明顯上升；受大理古城旅游及附屬用地開發(fā)、下關(guān)鎮(zhèn)城鎮(zhèn)及工業(yè)用地擴展的影響，海西中南部大理及下關(guān)兩鎮(zhèn)耕地地理集中度中小幅下降；受海東鎮(zhèn)上登工業(yè)園區(qū)、大理技師學院的修建及雙廊鎮(zhèn)南詔風情等旅游項目的影響，海東地區(qū)耕地地理集中度整體下降。這表明，地形平坦的海西中北部壩區(qū)依托優(yōu)越的自然條件，率先成為環(huán)洱海地區(qū)耕地開發(fā)的重點；隨城鎮(zhèn)化及工業(yè)化的發(fā)展，城鎮(zhèn)、旅游、工商及公共基礎設施用地的需求對耕地布局及數(shù)量的影響開始顯現(xiàn)，這在作為新興城區(qū)的海東南部及海西中南部地區(qū)尤為突出。

4.2 農(nóng)業(yè)勞動力地理集中度分布格局

環(huán)洱海地區(qū)農(nóng)業(yè)勞動力空間布局的區(qū)域差異明顯，海西地區(qū)地理集中度雖遠高于海東，但東西差異度略低于耕地。農(nóng)業(yè)勞動力密集區(qū)主要沿洱海西岸、S221及G214下關(guān)鎮(zhèn)北部段至上關(guān)鎮(zhèn)南部段呈帶狀集聚，隨著與洱海間距離的增加，農(nóng)業(yè)勞動力地理集中度呈圈層遞減。

從集中指數(shù)看，農(nóng)業(yè)勞動力的分布呈明顯的區(qū)域差異。海西中南部的下關(guān)鎮(zhèn)北部、大理及銀橋鎮(zhèn)主要由第一、二、三級村莊構(gòu)成，為農(nóng)業(yè)勞動力集中分布區(qū)；隨時間演進，海西中北部喜洲、灣橋兩鎮(zhèn)由農(nóng)業(yè)勞動力相對集中分布區(qū)變?yōu)榧蟹植紖^(qū)；海西太邑鄉(xiāng)及海東四鎮(zhèn)為第四、五級村莊的主要分布區(qū)，為農(nóng)業(yè)勞動力稀疏分布區(qū)；值得注意的是，海西南北兩側(cè)的上關(guān)及下關(guān)兩鎮(zhèn)分別受自然及社會經(jīng)濟發(fā)展的影響，內(nèi)部南北分異突出，農(nóng)業(yè)勞動力主要在G5611以南及G320以北集聚。

從集中區(qū)域看，隨時間演進，農(nóng)業(yè)勞動力整體呈現(xiàn)出由分散轉(zhuǎn)為集中的布局態(tài)勢。25年間，海西地區(qū)由以中南部為主的單核式布局變?yōu)榉謩e以中南部和中北部為核心的雙核分布；海東地區(qū)空間布局變動較小，呈分段集中態(tài)勢，主要在環(huán)海東路雙廊及挖色鎮(zhèn)段沿線、海東鎮(zhèn)上登工業(yè)園區(qū)及金梭島旅游景區(qū)、鳳儀鎮(zhèn)大理經(jīng)濟技術(shù)開發(fā)區(qū)鳳儀片區(qū)等地域集中。

交通通達度對農(nóng)業(yè)勞動力地理集中度的影響日趨凸顯。隨時間演進，海西G214沿線的沙坪村、上關(guān)村、慶洞村、云峰村及海東環(huán)海東路沿線的向陽村均由農(nóng)業(yè)勞動力稀疏的四級村莊（0.5

4.3 耕地與農(nóng)業(yè)勞動力空間分布耦合特征

為深入了解該區(qū)耕地與農(nóng)業(yè)勞動力空間分布的辯證關(guān)系，采用相關(guān)分析和曲線擬合分析對二者的地理集中分布關(guān)系進行測度。由表3發(fā)現(xiàn)，25年間，耕地與農(nóng)業(yè)勞動力地理集中度間的相關(guān)系數(shù)均在0.75以上，呈顯著的正相關(guān)，表明二者主體變動一致；曲線擬合顯示耕地與農(nóng)業(yè)勞動力地理集中度呈三次曲線關(guān)系，R2大于0.79而小于0.93，農(nóng)業(yè)勞動力地理集中度對耕地的解釋在79%～93%之間，隨時間演進，R2呈升降交替的“W”變動，表明環(huán)洱海地區(qū)農(nóng)業(yè)勞動力的空間集散對耕地的地理分布影響明顯但不能完全決定，從側(cè)面印證了該區(qū)耕地與農(nóng)業(yè)勞動力分布重心高頻區(qū)的“相交不重疊”。

4.4 耕地與農(nóng)業(yè)勞動力空間分布耦合類型

為詳細分析環(huán)洱海地區(qū)各行政村耕地與農(nóng)業(yè)勞動力的空間耦合特征，采用不一致指數(shù)對該區(qū)耕地與農(nóng)業(yè)勞動力的空間匹配關(guān)系進行表征[31]。不一致指數(shù)計算公式如下：

基于以上公式分別計算環(huán)洱海地區(qū)1988年和2012年的不一致指數(shù)，并依據(jù)不一致指數(shù)的大小將該區(qū)96個行政村分為5種耦合類型（圖4）。

（1） 耕地優(yōu)勢集聚型（I>1.4）。此類村莊指的是村域范圍內(nèi)耕地與農(nóng)業(yè)勞動力的空間分布明顯失調(diào)，耕地地理集中度顯著高于農(nóng)業(yè)勞動力。隨著時間演進，此類村莊數(shù)量明顯增加，由1988年的13個增至2012年的22個；空間布局的區(qū)域差異縮小，由“雙核式”集中于海西中部的銀橋、灣橋及北部的上關(guān)三鎮(zhèn)變?yōu)椤岸喽问健狈植加诤Ｎ髦胁康你y橋和灣橋、南部太邑鄉(xiāng)及海東北部的雙廊鎮(zhèn)。

（2）農(nóng)業(yè)勞動力分布滯后型（1.1

（3）耕地-農(nóng)業(yè)勞動力空間協(xié)調(diào)型（0.95）、均較低（R≤0.5）或均中等（0.5

（4）耕地分布滯后型（0.6

（5）農(nóng)業(yè)勞動力優(yōu)勢集聚型（I≤0.6）。此類村莊指的是村域范圍內(nèi)農(nóng)業(yè)勞動力地理集中度顯著高于耕地，耕勞分布明顯失調(diào)。隨時間演進，此類村莊數(shù)量大幅增加，空間布局由相對集中變?yōu)檠亟煌ǜ删€擴散。25年間，此類村莊由集中于北部的喜洲、雙廊兩鎮(zhèn)變?yōu)橹饕植加诒辈康纳详P(guān)、雙廊及南部的下關(guān)三鎮(zhèn)，呈明顯的南向擴展態(tài)勢，海西地區(qū)變動幅度明顯大于海東。

5 結(jié)論與討論

綜合運用區(qū)域重心、地理集中度及不一致指數(shù)等模型對環(huán)洱海地區(qū)耕地與農(nóng)業(yè)勞動力的空間耦合關(guān)系進行測度，得出以下結(jié)論：

（1） 環(huán)洱海地區(qū)耕地與農(nóng)業(yè)勞動力重心均位于幾何重心西北部，二者空間分布失衡，海西地區(qū)耕地與農(nóng)業(yè)勞動力規(guī)模明顯高于海東。隨時間演進，耕地分布重心整體呈向西北為主、東南為輔的變動態(tài)勢而農(nóng)業(yè)勞動力分布重心的階段性變動明顯，整體向西南偏移。

（2） 隨時間演進，環(huán)洱海地區(qū)耕地與農(nóng)業(yè)勞動力重心的空間重疊性與變動一致性均明顯提升，重心耦合態(tài)勢日趨明顯，兩重心空間分布的偏離度明顯減小。

（3） 環(huán)洱海地區(qū)耕地與農(nóng)業(yè)勞動力地理集中度均呈明顯的東西分異，海西整體遠高于海東，洱海西岸、S221及G214沿線分別為耕地與農(nóng)業(yè)勞動力的主次要集聚軸。隨時間推移，耕地向海西中北部集聚態(tài)勢明顯，海西南部及海東耕地地理集中度有所下降；農(nóng)業(yè)勞動力整體呈現(xiàn)由分散轉(zhuǎn)為集中的布局態(tài)勢，海西地區(qū)向中北部及大理鎮(zhèn)集聚而海東地區(qū)的交通通達度特別是高速公路及環(huán)海觀景公路對農(nóng)業(yè)勞動力地理集中度的影響日趨凸顯。

（4） 環(huán)洱海地區(qū)耕地與農(nóng)業(yè)勞動力的地理集中度呈顯著正相關(guān)，農(nóng)業(yè)勞動力的空間集散對耕地的地理分布影響明顯，各行政村間耕地與農(nóng)業(yè)勞動力空間分布的耦合存在明顯差異。海西中部一直為不一致指數(shù)高值村的主要分布區(qū)，隨時間演進，該類村莊向海西南北兩側(cè)及海東主要交通干線與工業(yè)園區(qū)附近擴散；空間協(xié)調(diào)型村莊所占比重由33.3%降為13.5%，村域間耕地與農(nóng)業(yè)勞動力空間分布的協(xié)調(diào)度整體下降；不一致指數(shù)低值村空間分布的南北差異縮小，海東地區(qū)數(shù)量明顯提升而海西地區(qū)數(shù)量變動幅度較大。

在城鎮(zhèn)化與非農(nóng)化發(fā)展背景下，作為城鄉(xiāng)轉(zhuǎn)型發(fā)展過程中的重要組成部分，耕地、農(nóng)業(yè)勞動力等生產(chǎn)要素的流動特征、格局及驅(qū)動力必將是未來鄉(xiāng)村地理研究的熱點。由于環(huán)洱海地區(qū)耕地紅線嚴格、生態(tài)保護嚴峻、旅游商業(yè)化程度高且多民族構(gòu)成的區(qū)域環(huán)境，綜合考慮農(nóng)業(yè)勞動力的兼業(yè)化與就業(yè)彈性、耕地與生態(tài)保護、工商發(fā)展與民族扶持等因素，探討不同區(qū)域適宜的勞耕彈性系數(shù)和最小勞均耕地面積將是本文進一步研究的重點。

參考文獻：

[1] 劉彥隨，李裕瑞. 中國縣域耕地與農(nóng)業(yè)勞動力變化的時空耦合關(guān)系[J]. 地理學報，2010，65（12）：1602-1612.

[2] 高凱，周志翔，楊玉萍，等. 基于Ripley K函數(shù)的武漢市景觀格局特征及其變化[J]. 應用生態(tài)學報，2010，21（10）：2621-2626.

[3] 曾凌云，王鈞，王紅亞. 基于GIS和Logistic回歸模型的北京山區(qū)耕地變化分析與模擬[J]. 北京大學學報：自然科學版，2009，45（1）：165-170.

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

[5] Ramankutty N， Foley J A. Estimating historical changes in global land cover： Croplands from 1700 to 1990[J]. Global Biogeochemical Cycles，1999，13（4）：997-1027.

[6] 趙曉麗，張增祥，汪瀟，等. 中國近30a耕地變化時空特征及其主要原因分析[J]. 農(nóng)業(yè)工程學報，2014，30（3）： 1-11.

[7] Guan D J， Gao W J， Watari K， et al. Land use change of Kitakyushu based on landscape ecology and Markov model[J]. Journal of Geographical Sciences，2008，18（4）：455-468.

[8] 曹雋雋，周勇，葉青清，等. 江漢平原耕地資源空間格局變化分析[J]. 經(jīng)濟地理，2013，33（11）：130-135.

[9] Skinner M W， Kuhn R G， Joseph A E. Agricultural land protection in China： A case study of local governance in Zhejiang Province[J]. Land Use Policy，2001，18（4）：329-340.

[10] 劉彥彤，張延軍，趙玲. 長春市耕地動態(tài)變化及其驅(qū)動力分析[J]. 地理科學，2011，31（7）：868-873.

[11] 車明亮，聶宜民，劉登民，等. 區(qū)域耕地數(shù)量變化預測方法的對比研究[J].中國土地科學，2010，24（5）：13-18.

[12] 張豪，羅亦泳，張立亭，等. 基于遺傳算法最小二乘支持向量機的耕地變化預測[J]. 農(nóng)業(yè)工程學報，2009，25（7）：226-231.

[13] 蔡昉. 中國的二元經(jīng)濟與勞動力轉(zhuǎn)移：理論分析與政策建議[M]. 北京：中國人民大學出版社，1990.

[14] 李小建. 農(nóng)戶地理論[M]. 北京：科學出版社，2009.

[15] 劉焱序，任志遠.基于區(qū)域地形起伏度模型的陜西農(nóng)村勞動力時空格局[J].山地學報，2012，30（4）：431-438.

[16] 黃敦平. 農(nóng)村勞動力流動微觀決策分析[J]. 人口學刊，2016，38（5）：54-59.

[17] 劉暢，許菁. 農(nóng)村勞動力供給的影響因素分析[J]. 中國人口·資源與環(huán)境，2016，26（5）：389-392.

[18] 馬草原. 非農(nóng)收入、農(nóng)業(yè)效率與農(nóng)業(yè)投資—對我國農(nóng)村勞動力轉(zhuǎn)移格局的反思[J].經(jīng)濟問題，2009（7）： 66-69.

[19] 龍翠紅，易承志.中國農(nóng)村勞動力流動的格局與現(xiàn)代新型農(nóng)民培養(yǎng)[J]. 經(jīng)濟問題探索，2011（1）：104-110.

[20] 姜曉清，張 果，李叢穎，等.成都市耕地與農(nóng)業(yè)勞動力的時空耦合關(guān)系[J]. 湖北農(nóng)業(yè)科學，2016，55（14）：3790-3794.

[21] 鄧楚雄，謝炳庚，潘艷平，等. 典型農(nóng)業(yè)縣域轉(zhuǎn)出勞動力對耕地處置方式的影響因素分析——以湖南省新化縣為例[J]. 經(jīng)濟地理，2015，35（9）：134-140.

[22] 楊俊，楊鋼橋，胡賢輝. 農(nóng)業(yè)勞動力年齡對農(nóng)戶耕地利用效率的影響—來自不同經(jīng)濟發(fā)展水平地區(qū)的實證[J]. 資源科學，2011，33（9）：1691-1698.

[23] 于慧. 人地關(guān)系視角下陜西省耕地與農(nóng)業(yè)勞動力協(xié)調(diào)研究[D]. 西安：西北大學，2012.

[24] 馬彩虹，趙先貴. 陜西省人口-耕地-糧食系統(tǒng)耦合態(tài)勢研究[J].干旱地區(qū)農(nóng)業(yè)研究，2005，23（5）：217-221.

[25] 湯進華，陳志，郝漢舟，等.湖北省耕地與農(nóng)業(yè)勞動力變化的耦合關(guān)系探究[J].湖北農(nóng)業(yè)科學，2015，54（21）：5481-5488.

[26] 張磊，武友德，李君，等.云南環(huán)洱海地區(qū)農(nóng)村聚落空間分布變化及特征分析[J]. 水土保持研究，2016，23（6）：316-321.

[27] 劉影，肖池偉，李鵬，等.1978-2013年中國糧食主產(chǎn)區(qū)“糧-經(jīng)”關(guān)系分析[J].資源科學，2015，37（10）：1891-1901.

[28] 張磊，武友德，李君. 泛珠江三角洲經(jīng)濟圈人口與經(jīng)濟的空間分布關(guān)系研究[J].西北人口，2015，36（5）：43-48.

[29] 樊杰，陶岸君，呂晨.中國經(jīng)濟與人口重心的耦合態(tài)勢及其對區(qū)域發(fā)展的影響[J].地理科學進展，2010，29（1）：87-95.

[30] 鐘業(yè)喜，劉運偉，賴格英. 紅三角經(jīng)濟圈人口與經(jīng)濟空間分布關(guān)系研究[J]. 華南師范大學學報：自然科學版，2012，44（1）： 118-123.

[31] 鐘業(yè)喜，陸玉麒. 鄱陽湖生態(tài)經(jīng)濟區(qū)人口與經(jīng)濟空間耦合研究[J]. 經(jīng)濟地理，2011，31（2）：195-200.

Abstract： Based on the models of regional gravity center， geographic concentration index and inconsistency index， from 1988 to 2012， the spatial distribution and coupling relationship between arable land and agricultural labor force in Erhai lake Area are analyzed. The results indicate that： with time going on， the evolutionary trend of the gravity center of arable land was mainly to the northwest， while southeast as a secondary； the phased change of the gravity center of agricultural labor force is obviously， overall shifted to the southwest； In 25 years， the space-overlaps and changes in consistency of arable land and agricultural labor force gravity were both significantly increased， the coupling trend of the two gravity centers are becoming more and more obvious. The degree of geography concentration of arable land and agricultural labor force in Haixi area are both far above those in Haidong， the west coast of Erhai Lake， S221 and G214 are the main and minor agglomeration axis of arable land and agricultural labor force respectively； Since 1988， the arable land gathered to the central and northern regions of Haixi， while the agricultural labor force showing a change from decentralized to centralized layout， as a whole， there was a highly positive correlation between them； With time going on， the coordination degree of spatial distribution of arable land and agricultural labor force in the village area decreased obviously， according to the inconsistency index， the Erhai lake Area can be divided into five types of villages which are arable land advantage agglomeration type， agricultural labor force lag distribution type， arable land-agricultural labor force spatial coordination type， arable land lag distribution type and agricultural labor force advantage agglomeration type.

Key words： arable land； agricultural labor force； spatial coupling； Erhai lake Area