環洱海地區耕地與農業勞動力空間分布特征及其耦合關系

張磊 武友德 李君

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

關鍵詞：耕地；農業勞動力；空間耦合；環洱海地區

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

0 引言

在城鎮化與工業化快速發展的背景下，區域耕地快速流動的同時，農業勞動力能否隨之非農轉移，促進勞耕比例的適應性調節，已成為保障鄉村系統有序發展的重要前提[1]。作為反映區域農業發展狀況的重要指標，耕地與農業勞動力已成為鄉村地域人地關系研究中的重要內容。

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

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

1 研究區域與數據來源

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

本文數據主要由經緯度空間數據和區域發展屬性數據構成，為保障數據連續性與精確性，選取環洱海地區96個行政村為研究對象。基于實地調研分析，采用GPS地理定位法，提取各行政村幾何重心的經緯度，并基于ArcGIS10.1構建行政村坐標體系網作為空間數據；取環洱海地區各行政村1988年以來的年末集體所有耕地總面積和年末匯總勞動力數作為屬性數據，屬性數據主要取自《大理市農村經濟收支情況統計年報》（1988-2013）、《大理市統計年鑒》（1988-2013）、《洱源縣國民經濟和社會發展統計年鑒》（1988-2013）及課題組實地調研資料①。

2 研究方法與理論模型

2.1 重心測度模型

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

3 耕地與農業勞動力空間重心耦合分析

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

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

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

3.2 農業勞動力空間分布重心變動分析

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

3.3 耕地與農業勞動力重心空間耦合特征

基于公式（2）、（3），通過測度兩重心的空間重疊性及變動一致性，對兩重心的空間耦合態勢進行分析。由表2發現，1988年～2012年，環洱海地區耕地與農業勞動力重心的空間耦合呈以下特征：

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

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

3）隨時間演進，環洱海地區耕地與農業勞動力重心的空間重疊性與變動一致性均明顯提升，重心耦合態勢日趨明顯，兩重心空間分布的偏離度明顯減小。

4 耕地與農業勞動力空間分布耦合分析

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

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

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

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

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

4.2 農業勞動力地理集中度分布格局

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

從集中指數看，農業勞動力的分布呈明顯的區域差異。海西中南部的下關鎮北部、大理及銀橋鎮主要由第一、二、三級村莊構成，為農業勞動力集中分布區；隨時間演進，海西中北部喜洲、灣橋兩鎮由農業勞動力相對集中分布區變為集中分布區；海西太邑鄉及海東四鎮為第四、五級村莊的主要分布區，為農業勞動力稀疏分布區；值得注意的是，海西南北兩側的上關及下關兩鎮分別受自然及社會經濟發展的影響，內部南北分異突出，農業勞動力主要在G5611以南及G320以北集聚。

從集中區域看，隨時間演進，農業勞動力整體呈現出由分散轉為集中的布局態勢。25年間，海西地區由以中南部為主的單核式布局變為分別以中南部和中北部為核心的雙核分布；海東地區空間布局變動較小，呈分段集中態勢，主要在環海東路雙廊及挖色鎮段沿線、海東鎮上登工業園區及金梭島旅游景區、鳳儀鎮大理經濟技術開發區鳳儀片區等地域集中。

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

4.3 耕地與農業勞動力空間分布耦合特征

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

4.4 耕地與農業勞動力空間分布耦合類型

為詳細分析環洱海地區各行政村耕地與農業勞動力的空間耦合特征，采用不一致指數對該區耕地與農業勞動力的空間匹配關系進行表征[31]。不一致指數計算公式如下：

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

（1） 耕地優勢集聚型（I>1.4）。此類村莊指的是村域范圍內耕地與農業勞動力的空間分布明顯失調，耕地地理集中度顯著高于農業勞動力。隨著時間演進，此類村莊數量明顯增加，由1988年的13個增至2012年的22個；空間布局的區域差異縮小，由“雙核式”集中于海西中部的銀橋、灣橋及北部的上關三鎮變為“多段式”分布于海西中部的銀橋和灣橋、南部太邑鄉及海東北部的雙廊鎮。

（2）農業勞動力分布滯后型（1.1

（3）耕地-農業勞動力空間協調型（0.95）、均較低（R≤0.5）或均中等（0.5

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

（5）農業勞動力優勢集聚型（I≤0.6）。此類村莊指的是村域范圍內農業勞動力地理集中度顯著高于耕地，耕勞分布明顯失調。隨時間演進，此類村莊數量大幅增加，空間布局由相對集中變為沿交通干線擴散。25年間，此類村莊由集中于北部的喜洲、雙廊兩鎮變為主要分布于北部的上關、雙廊及南部的下關三鎮，呈明顯的南向擴展態勢，海西地區變動幅度明顯大于海東。

5 結論與討論

綜合運用區域重心、地理集中度及不一致指數等模型對環洱海地區耕地與農業勞動力的空間耦合關系進行測度，得出以下結論：

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

（2） 隨時間演進，環洱海地區耕地與農業勞動力重心的空間重疊性與變動一致性均明顯提升，重心耦合態勢日趨明顯，兩重心空間分布的偏離度明顯減小。

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

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

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

參考文獻：

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

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

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

[4] 任平，吳濤，周介銘. 基于GIS和空間自相關模型的耕地空間分布格局及變化特征分析——以成都市龍泉驛區為例[J]. 中國生態農業學報，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]. 農業工程學報，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]. 經濟地理，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] 劉彥彤，張延軍，趙玲. 長春市耕地動態變化及其驅動力分析[J]. 地理科學，2011，31（7）：868-873.

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

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

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

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

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

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

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

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

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

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

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

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

[23] 于慧. 人地關系視角下陜西省耕地與農業勞動力協調研究[D]. 西安：西北大學，2012.

[24] 馬彩虹，趙先貴. 陜西省人口-耕地-糧食系統耦合態勢研究[J].干旱地區農業研究，2005，23（5）：217-221.

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

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

[27] 劉影，肖池偉，李鵬，等.1978-2013年中國糧食主產區“糧-經”關系分析[J].資源科學，2015，37（10）：1891-1901.

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

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

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

[31] 鐘業喜，陸玉麒. 鄱陽湖生態經濟區人口與經濟空間耦合研究[J]. 經濟地理，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