青藏高原地區農村居民點空間演化特征及驅動力研究

李媛媛，李鋒，陳春,2*

（1.重慶交通大學建筑與城市規劃學院，重慶 400074；2.重慶交通大學生態人居與綠色交通研究中心，重慶 400074）

伴隨著城鎮化的推進，城市擴張、產業結構的轉變以及基礎設施的革新等問題的出現，對農村居民點空間形態造成了沖擊，改變了傳統農村居民點空間布局，引發了村鎮用地結構單一、規模小、破碎度高、缺乏規劃管理等問題[1]。在此背景下，農村居民點演化受到了眾多關注。在研究內容上，目前的研究主要包括了農村居民點空間格局演化特征及驅動機制[2-3]、空間布局優化[4]以及農村居民點擴張對生態的影響[5]等；在研究方法上，目前的研究中大多采用空間分析和景觀格局指數相結合的方法研究農村居民點的演化特征；在研究尺度上，在區域、縣域、村域尺度均有研究，而這些研究大多以平原地區[6]、東部經濟發達地區[7]、丘陵地區[8-9]、中部傳統農業區[10]典型的市（縣）為研究區域[11]。

青藏高原因其獨特的地形特征，有著地球“第三極”、世界屋脊之稱。并且由于其是亞洲大江大河的發源地，被稱為“亞洲水塔”。其特殊的地理氣候特征，對我國乃至全球的氣候、碳收支平衡、生物多樣性等生態環境問題具有重要的影響作用，是我國重要的生態安全屏障。青藏高原地區生態系統薄弱，人類活動對青藏高原地區環境系統的影響持續增強[12]。1980年后，隨著以家庭或聯合牧場為經營單元的季節輪牧制度的實施[13]，當地居民的生活發生轉變，由游牧生活轉變為定居，推動了農村居民點用地格局的劇烈變化。同時隨著城鄉一體化的推進，農村經濟進一步發展，青藏高原地區農村居民點用地不斷增加，人類活動日漸頻繁，為促進農村經濟的發展，一味地追求土地利用的最大效益，導致草地退化、生物多樣性減少等一系列的生態環境問題[14]。對青藏高原地區農村居民點的演化特征及驅動因素進行研究，對其農村居民點用地的合理規劃、青藏高原地區的生態系統維持有著重要的現實意義。

由于青藏高原地區發展相對落后，基礎資料獲取困難，且鄉村聚落規模較小，密度低，為農村居民點研究增加了難度，導致目前對青藏高原地區農村居民點演化及驅動因素研究薄弱。目前僅有部分學者針對青藏高原地區的城鎮進行了研究，包括城鎮體系的演變[15]、城鎮化時空演變格局[16-17]以及城鎮化與生態環境之間的關系[18]等，也有學者以縣、鎮、鄉尺度的居民點為研究對象，將城市聚落和鄉村聚落作為一個整體，研究青藏高寒牧區的聚落 演化特征及驅動因素[19]、城鄉建設用地時空轉移格局[20]等。專門以農村居民點為研究對象，分析青藏高原農村居民點的空間演化特征及驅動因素的研究相對缺乏。

基于此，本文以青藏高原地區為研究區域，以農村居民點為研究對象，采用空間分析和景觀格局分析相結合的方法，分析1980—2015年青藏高原地區農村居民點用地轉換、空間分布、形態結構以及規模特征變化，進而闡釋演化驅動機制，以期為青藏高原地區農村居民點合理規劃提供借鑒。

1 研究區概況

青藏高原地處我國的西南部，地理位置在26°00′～39°47′N、73°19′～104°47′E之間，北靠昆侖山、祁連山，南至喜馬拉雅山脈，西起帕米爾高原、喀喇昆侖山，東至橫斷山脈。作為世界范圍內地形最為復雜的地貌單元，其地勢西高東低，內部地形起伏度大，4 000 m以上的高海拔地區占整個青藏高原面積的73.11%[21]。由于青藏高原地區地形及氣候條件嚴苛，與同緯度地區相比，其人口密度較低，城鎮化水平極低，社會經濟也相對落后。在考慮研究具有代表性和數據的可獲得性的基礎上，本文并沒有將整個青藏高原定為研究區，僅以青藏高原的主體部分，即西藏自治區和青海省所轄區域（圖1）作為研究區域，面積約194.94 km2，占中國國土面積的19.8%。據統計，直至2019年青海省和西藏自治區共有鄉鎮1 100個，共有農村人口511萬人，占其總人口的53.28%，其中西藏自治區農村人口為240萬，占西藏總人口的68.38%，農村人口占比位居全國第一[22]。

2 研究方法和數據來源

2.1 研究方法

2.1.1 空間韻律指數 空間韻律指數即景觀格局指數，其特點是高度濃縮景觀格局信息，能夠反映其形態、規模及空間分布等特征[23]。農村居民點的景觀表現為不同的農村居民點所組成的大小、規模以及形狀等不同斑塊。本研究在結合研究區實際情況和各景觀指數含義的基礎上，選取斑塊面積（CA）、斑塊個數（NP）、斑塊密度（PD）、平均斑塊面積（MPS）、最大斑塊指數（LPI）、平均斑塊形狀指數（SHAPE_MN），平均斑塊分維數（FRAC_MN）和聚集度（AI）等指標對農村居民點的規模和形態結構變化進行研究。

2.1.2 核密度分析、平均最近鄰指數 1）核密度分析。通常使用柵格數據密度值反映點狀數據分布的密 集程度[24]。通常以研究區所有樣本點為圓心，通過設置合適的搜索半徑生成圓，距離圓心越近柵格密度值越大，農村居民點分布密度則越大，若沒有樣本點落入柵格數據中，則密度值為零[25]，以此方法測度農村居民點空間分布密度。核密度分析公式如下[26]：

式中：f(x，y）是位于（x，y）位置的核密度估計；n為觀測變量；h為平滑參數；k為核函數；di為位置距第i個觀測位置的距離。

2）平均最近鄰指數。通過將農村居民點之間最鄰近距離的平均值和隨機分布模式下最鄰近距離做比較，用其比值大小做農村居民點空間聚集性判斷依據。其比值ANN小于1，則農村居民點為集聚分布，若ANN大于1，則趨向于隨機分布，公式如下[27]：

2.1.3 空間關聯指數 1）全局聚類檢驗。全局聚類檢驗（Getis-Ord GeneralG）主要用于檢測居民點規模的全局性空間分布模式，即是高值空間集聚或低值空間集聚，公式如下[1]：

式中：wij是通過距離規則定義的空間權重；xi和xj分別是i和j區域的觀測值。

2）熱點探測分析。在式3的基礎上，同時引用ArcGIS10.2平臺的Getis-OrdGi*工具檢驗農村 居民點規模在空間上局部規模冷熱點分布，公式 如下[25]：

式中：n為農村居民點總數，wij、xi和xj的含義與式（3）相同。對Gi*(d)進行標準化為Z(Gi*)，當Z(Gi*)為正，居民點為高值聚集；反之則為低值聚集。

2.2 數據來源

本文采用的農村居民點斑塊數據來源于中科院資源環境科學與數據中心的1980、2000以及2015年土地利用/覆被柵格數據（https://www.resdc.cn/），空間分辨率為30 m，并經過ArcGIS平臺提取農村居民點類型所得。由于景觀格局軟件Fragstats4.2對柵格數據的行列數有限制，在進行景觀格局分析時，利用最近鄰法重采樣所獲取到的農村居民點柵格數據為150 m×150 m，故文中提及的農村居民點并非行政意義上的村鎮，而是多個農村居民點集聚在一起滿足150 m×150 m的面積范圍的農村居民點集合[28]。通過與原數據對比發現，采樣后的數據特征與原數據一致，查閱文獻后，發現已有研究的研究區面積遠小于本研究區的面積，但采用了空間分辨率為1 km的柵格數據進行了相關的研究[29]。因此，本研究將柵格數據重采樣為150 m，其尺度效應滿足研究需求。DEM數據來源于地理空間數據云（http://www.gscloud.cn/），空間分辨率為30 m，高程數據由DEM數據計算所得。社會經濟人口等數據來源于1980年、2000年和2015年《西藏統計年鑒》、《青海統計年鑒》。

3 結果與分析

3.1 農村居民點用地變化特征

ArcGIS10.2平臺的可視化工具能直觀的展現農村居民點的空間分布特征，并對其變化特點加以分析，因此對1980、2000和2015年青藏高原地區土地利用數據中農村居民點用地進行提取，并按屬性選擇出不同時間范圍內轉入農村居民點的面狀圖斑，因研究區平均海拔較高，且地域遼闊，人煙稀少，所提取的面狀圖斑無法滿足可視化條件，故采用要素轉點工具，對農村居民點的時空分布特征進行展示分析。如圖2所示，西寧和拉薩作為青海、西藏兩省的省會及首府城市，是青藏高原地區主要的政治、經濟中心，對其周邊近鄰范圍內的農村居民點具有強烈的帶動和吸引作用，1980年以來，農村居民點主要集中分布于西寧市和拉薩市及其周邊市縣；1980—2000年期間，新增的農村居民點主要集中分布于西寧市、海東市、海南藏族自治州、海西蒙古族藏族自治州及玉樹藏族自治州，該期間農村居民點的斑塊增長分布主要是在原有農村居民點近鄰區域呈現外圍擴張的態勢；2000—2015年間，農村居民點的斑塊增長主要分布在拉薩市南部、西寧市中部、昌都市東部、日喀則市中部及東部、林芝市中部、那曲市東南部和阿里地區西部，新增的農村居民點不再局限于原有近鄰區的范圍式擴張，而是在較偏遠的地區由其他用地類型轉入，以分散式點狀分布或小型集聚組團的形式呈現。

GIS平臺的可視化工具只能反應農村居民點整體的空間分布特點，無法探究其內在結構的變化特點，因此，對農村居民點及不同用地類型間的轉換過程進行定量分析，能夠準確的認識到農村居民點變化的內部原因。使用ArcGIS10.2的相交工具對農村居民點的轉入轉出進行分析，得到如表1所示的青藏高原地區農村居民點面積轉移表。由表1可知，1980—2015年間，耕地、草地及林地向農村居民點的轉入量大，其次是城鎮及其他建設用地、林地以及水域；農村居民點用地轉出量最大的是耕地和城鎮及其他建設用地，其次是草地、林地、水域；綜合來看，研究期內，耕地、未利用地和草地的凈轉入量較大，說明農村居民點用地擴張的首要原因是侵占耕地和草地，但也不乏對未利用地的合理開發利用；城鎮及其他建設用地的轉出量不斷增加，并且其是農村居民點用地凈轉入量為負的唯一用地類型，這與青藏高原地區城鎮化不斷推進，城市建設需求用地不斷加大有密切的關系。

表1 1980—2015年青藏高原地區農村居民點面積轉移情況（km2）Table 1 Transfer of village-town settlement area in Qinghai-Tibet region from 1980 to 2015 (km2)

綜上，農村居民點斑塊擴張態勢明顯，擴張方式多樣，由近鄰區的范圍式擴張到分散式點狀擴張；新增農村居民點主要來源于對耕地、未利用地和草地的侵占，為順應新農村建設要求，應當加強未利用地的使用，保障基本農田和生態紅線范圍內用地類型的穩定。

3.2 空間集聚演化特征分析

基于ArcGIS10.2平臺的Spatial Statistics工具計 算平均最近鄰ANN指數，由表2可知，研究區1980、 2000、2015年的ANN指數值均小于1，表明青藏高原地區農村居民點呈現聚集分布的狀態。

表2 1980—2015年青藏高原地區農村居民點ANN指數Table 2 ANN index of villages and towns in Qinghai-Tibet region from 1980 to 2015

根據農村居民點點狀數據，利用Kernel工具，分別生成研究區1980、2000、2015年間核密度分布圖，因青藏高原地區地域廣闊，人口稀少，經多次實驗，確定搜索半徑為25 km，能夠使得青藏高原地區核密度圖達到較好效果。如圖3所示，青藏高原地區農村居民點分布主要以局部集聚的形式呈現，其中青海省東部和南部及西藏南部地區是農村居民點高密度集聚的區域，而以西寧市、玉樹市和拉薩市為中心的近鄰圈則一直是農村居民點集聚分布的核心。1980年以來，西寧市及近鄰圈的農村居民點主要集聚在西寧市、海東市、海南藏族自治州和海北藏族自治州，而隨著城鎮化不斷推進，城鎮建設用地需求量不斷加大，農村居民點用地也在逐年上升，但高密度集聚區域一直處于不變狀態，主要與未利用地的轉入貢獻有關；拉薩市及近鄰圈的農村居民點高密度集聚區域面積在2015年有所減小，主要因為自2005年之后，拉薩市及近鄰圈達到城鎮建設大發展第二階段，重點加強了交通、水利、能源及市政等基礎設施建設，其中不乏對農村居民點的侵占與改建。

3.3 形態特征分析

農村居民點的形態特征主要從整體的斑塊特征和農村居民點自身形態進行分析。在Fragstats4.2平臺的支持下，從景觀格局指數中選取平均斑塊形狀指數（SHAPE_MN）、平均斑塊分維數（FRAC_MN）和聚合度（AI）三個指標對農村居民點的形態特征進行分析，其中SHAPE_MN和FRAC_MN表征農村居民點自身的形態特征，AI表征農村居民點斑塊整體形態特征。根據表3得知，平均形狀指數在1980—2015年一直處于上升狀態，表明青藏高原地區農村居民點斑塊形狀逐漸趨于不規則狀態；平均斑塊分形指數一般取值在1～2之間，若值為1，則表示形狀為最簡單的規則形狀，其值越接近2，則分形特征越明顯，青藏高原地區平均斑塊分形指數十分接近1，且研究期內變化幅度很小，表明人為因素對農村居民點的影響較大，其自身形態較為規則。聚合度可以表明農村居民點的集聚程度，其值越大，表明農村居民點越具有集聚的特征，研究區聚合度處于逐年增長狀態，表明農村居民點存在聚集發展的態勢。可見，在城鎮化的帶動下，研究期內青藏地區農村居民點斑塊逐漸趨于不規則形態，且存在明顯的集聚特征。

表3 1980—2015年青藏高原地區形態特征景觀指數Table 3 Landscape index table of morphological characteristics in Qinghai-Tibet region from 1980 to 2015

3.4 規模演變特征

3.4.1 規模增減特征 農村居民點的規模增減特征主要從整體用地規模變化和局部斑塊自身規模變化兩方面考量，通常以斑塊面積（CA）、平均斑塊面積（MPS）、斑塊個數（NP）、最大斑塊所占景觀面積比例（LPI）以及斑塊密度（PD）等指標表征農村居民點規模增減變化特征。如表4所示，1980—2015年間，農村居民點斑塊面積增加了22.21%，平均斑塊面積增長了25.40%，表明農村居民點用地規模一直在擴張；除此之外，最大斑塊所占景觀面積比例也在逐年增加，表明農村居民點規模呈集聚發展的趨勢；研究期內，斑塊個數逐年減少，1980—2015年間減少了2.54%，表明農村居民點斑塊破碎度減小。1980—2015年間，青藏高原地區 農村居民點斑塊面積不斷增長，而斑塊個數逐年 減少，可見農村居民點斑塊在研究期內逐漸趨于連續性。

表4 1980—2015年青藏高原地區規模變化景觀指數Table 4 Landscape index of scale change in Qinghai-Tibet region from 1980 to 2015

3.4.2 規模分異特征 青藏高原地區農村居民點斑塊面積標準差的Moran’sI指數1980年為0.090 932（E=0.000 005，Z=42.912 21）、2000年 為0.095 155（E=0.000 004，Z=44.942 93）、2015年 為0.012 254（E= -0.000 186，Z=6.354 95），在研究期內Moran’sI指數顯著為正，由此可見青藏高原地區農村居民點的規模布局有明顯的空間集聚特征。由于其常用來描述屬性值之間的相似程度以及空間上的分布情況，對于究竟是高值還是低值的規模空間集聚不能準確地描述。基于此，本研究使用全局空間聚類檢驗的GeneralG統計量對青藏高原地區高/低規模空間集聚特征進一步分析。如表5所示，研究期內GennralG統計量觀測值均大于期望值，全局空間聚類檢驗顯著。由此可見，青藏高原地區的高值聚集現象更為明顯，農村居民點規模分布出現了多個高值集聚區。

表5 1980—2015年青藏高原地區全局空間聚類檢驗Table 5 Global spatial clustering test of Qinghai-Tibet region from 1980 to 2015

對青藏高原地區農村居民點局部規模分異的可視化分析，有助于直觀的展現農村居民點規模的空間分布，通過ArcGIS10.2平臺的“熱點探測工具”，對青藏高原地區農村居民點面積屬性進行統計得出規模局域空間熱點分布圖（圖4），結果顯示： 1）1980—2015年分布于西寧市、海南藏族自治州及海西蒙古族自治州等城市近鄰區的高值集聚區面積明顯減少，主要轉向了次高值集聚區，農村居民點斑塊由較大規模轉變為適宜農村建設發展的合理規模分布狀態，在鄉村振興戰略的領導下，青海省2010年啟動“千村建設、百村示范”，建設高原美麗宜居村鎮計劃，通過生態保護與建設，村莊環境整治等合理規劃發展村鎮建設的初步成果；拉薩市近鄰區的部分農村居民點自2015年轉為線狀分散式高值集聚，表明農村居民點斑塊規模變大，在人才引進政策的帶動下，西藏自治區形成以拉薩為中心帶動周邊縣市經濟發展的長久計劃，農村居民點規模于2015年形成線狀分散式高值集聚，是充分考慮地形限制和交通便利條件下的綜合結果。2）果洛藏族自治州自2015年由低值轉為次高值集聚，斑塊數量增多，規模變大，玉樹自治州自1980年以來一直有大面積的低值集聚農村居民點，主要在于緊鄰四川西部的阿壩自治州，有G214國道和S101省道貫通，是青海省與四川省經濟文化交流的必經之路，在城鎮化加速推進和鄉村振興戰略的共同作用下，兩地都有大規模的農村居民點分布，但果洛藏族自治州有著距省會城市西寧市距離較近的區位優勢，受西寧市的輻射帶動作用較強，經濟發展優于玉樹自治州，在2000年之后逐漸轉變為次高值集聚分布，農村居民點斑塊規模變大，存在一定的規劃管理。3）隨著政策的逐漸落實，交通逐漸便利，原本偏僻的阿里地區、那曲市等在2015年也形成了農村居民點的點狀次高值集聚，表明農村居民點斑塊規模變大，隨著交通、政策的支持力度增大，促進了偏遠地區的發展，優化了農村居民的生活條件，從而促進了農村居民點的集聚。

4 驅動因素分析

4.1 自然因素

4.1.1 地形地貌 青藏高原作為世界上最復雜的地形單元，其大部分地區海拔都在4000 m以上，農村居民點的布局深受區域高程和坡度的影響。高程和坡度作為宏觀的地理背景，對農村居民點的布局和演化有直接的關系。

高程因素是農村居民點形成的主要原因之一，是影響農村居民點演化的主要自然因素，一方面，合適的高程可構建農村居民點形成的基礎條件；但在另一方面也限制著農村居民點的發展變化。本文將青藏高原地區高程劃分為DEM < 2 000 m為低區，2 000 m < DEM < 4 000 m為中區，DEM > 4 000 m為高區，選取合適的景觀指標進行計算。如表6所示，青藏高原地區農村居民點斑塊面積、斑塊所占景觀面積比例、斑塊個數及斑塊密度的平均最大值主要集中在中區，并隨著高程的不斷升高而增大，但當高程達到4 000 m時，斑塊面積、斑塊個數、斑塊密度等指標數值大幅減小，主要是因為高海拔地區不適宜人類的居住。由于青藏高原的平均海拔在4 000 m以上，高程小于2 000 m的區域值很小，所以在低區并未呈現最大值；聚合度（AI）隨著DEM的升高在高區出現最大值，表明高海拔地區農村居民點處于小規模、高集聚分布狀態。從中區的農村居民點特征來看，1980—2015年斑塊面積、斑塊所占景觀面積比例以及聚合度均有所增長，隨著時間的增長，農村居民點在中區更為集中。由此可見，高程對農村居民點布局的影響顯著，青藏高原地區的農村居民點主要分布在海拔2 000～4 000 m的中區，并且大部分農村居民點布局逐漸向中區集聚，而在高區部分表現出了農村居民點小規模、高集聚分布的狀態。

表6 1980—2015年青藏高原地區不同高程緩沖區農村居民點景觀特征Table 6 Settlement landscape characteristics of villages and towns at different elevations in the buffer zone of Qinghai-Tibet region from 1980 to 2015

除了高程外，坡度也是影響農村居民點布局的重要自然因素。本研究根據地貌坡度劃分標準和研究區的地形特征，將坡度劃分為小于8°、8°～25°、大于25°三個等級，進而探究不同坡度范圍下農村居民點分布差異。如表7所示，青藏高原地區超半數的農村居民點分布在坡度小于8°地勢較為平坦的區域，地勢平坦的地區更適合居住與生活。坡度為8°～25°相對陡峭的區域也有農村居民點的集聚，主要是因為青藏高原本身地形限制，且在地勢稍有坡度的地方更適合畜牧業的發展，部分居民選擇在此生活居住，但由于坡度較大，在這個區域內的農村居民點分布十分零散，且規模較小。坡度大于25°的區域極少有農村居民點的布局，主要是因為坡度過大，不適合居住和生活。從聚合度來看，不同坡度范圍的農村居民點均呈現出了分布零散，規模較小的特點，這也與青藏高原本身的地形起伏大有密切的關系。對于坡度小于8°的區域，1980—2015年農村居民點的斑塊面積、聚合度以及斑塊所占景觀面積比例均呈上升趨勢，可見隨著時間的發展，青藏高原地區的農村居民點不再是小規模、零散的分布，而是在地勢較為平緩的地區逐漸朝著大規模、高集聚的方向發展。

表7 1980—2015年青藏高原地區不同坡度緩沖區農村居民點景觀特征Table 7 landscape characteristics of rural residential areas in buffer zones with different slopes in Qinghai Tibet Plateau from 1980 to 2015

4.1.2 水源 水系作為影響農村居民點布局的重要影響因素，是驅動力研究的重要因素。本研究分別以2 km、4 km和6 km建立水系的緩沖區，并采用景觀格局指數相關指標分析不同緩沖區距離內農村居民點的分異特征（表8）。從整體上來看，研究期內農村居民點沿水系分布明顯，在距離水系2 km以內區域，農村居民點斑塊面積、斑塊個數、斑塊密度、形態特征以及聚合度均較另外兩個緩沖區分組大，由此可見農村居民點主要集中于距離水系0～2 km內區域，且集聚特征明顯。在0～2 km緩沖區內，1980—2015年斑塊面積和聚合度都呈現出增長趨勢，而斑塊個數在逐步減少，由此可見農村居民點逐漸在靠近水系的區域集中，且規模有擴大的趨勢。雅魯藏布江作為水能資源僅次于長江的世界海拔最高的河流，其中游地區憑借其豐富的水資源和形成的平原地形，使其沿岸有大規模的農村居民點集聚。靠近水源，一方面可以增加用水的便捷性，保障居民在生活、農耕等方面的用水；另一方面，河谷地帶地形平緩，適合居民生活、耕作，因此農村居民點在距離水系0～2 km的區域分布密集。

表8 1980—2015年青藏高原地區不同距離水系緩沖區農村居民點景觀特征Table 8 Settlement landscape characteristics of villages and towns with different distance from river system buffer zones in Qinghai-Tibet region from 1980 to 2015

4.2 社會經濟因素

4.2.1 區位因素 區位因素在社會經濟大力發展的今天，對農村居民點的演化起到至關重要的作用，距離交通干線遠近及城鎮的遠近，影響著農村居民點的發展。農村居民點的分布形態對交通干線和交通流量也有反作用，影響著交通干線的走向；城鎮的輻射帶動作用對周邊農村居民點也有一定的帶動 作用。

道路作為連接各個居民點的紐帶，在農村居民點聚集過程中有著舉足輕重的作用。因此本研究以2 km為半徑建立道路的多個緩沖區，從而分析不同道路緩沖區距離內農村居民點的變化特征。由于鐵路除站點以外的沿線與農村居民點演化影響關系并非很大，本研究將公路作為主要驅動因素進行分析，所選公路為青藏高原地區主要國道、省道及縣道。

從表9可以看出，農村居民點在道路兩側呈帶狀分布，且1980—2015年新增農村居民點也在該范圍內形成擴展區。從不同緩沖區來看，農村居民點在距離公路0～2 km區域內大規模集聚，隨著與道路距離的增加，農村居民點斑塊面積、斑塊個數以及斑塊密度都逐漸變小，在道路緩沖區0～2 km范圍內新增農村居民點的面積占0～6 km緩沖區內所有新增農村居民點面積的78%。由此可見，道路對農村居民點的集聚有著重要作用，與此同時，農村居民點也影響著公路的走向與分布。在0～2 km緩沖區內，農村居民點的形態指數和聚合度均呈上升趨勢，可見隨著時間的增長，農村居民點在道路兩旁集聚規模不斷增大，集聚特征更為明顯。2009年國道109正式通車，帶動了地處偏遠的那曲市與外界的聯系，109國道青藏段因其路況較好成為外地人進藏的首選，同時109國道上的客運班車也加強了當地與外地的聯系，109國道通車后，那曲市出現農村居民點的點狀集聚中心。由此可見，道路是促成農村居民點發展的重要驅動因素。

表9 1980—2015年青藏高原地區不同距離道路緩沖區農村居民點景觀特征Table 9 Landscape characteristics of villages and towns with different distance from road buffer zones in Qinghai-Tibet region during 1980-2015

城鎮是區域的經濟發展中心，農村居民點與城鎮之間的距離是農村居民點演化的重要影響因子。因此，研究以2 km為半徑建立城鎮距離緩沖區，各級緩沖區內農村居民點分異特征如表10所示，隨著距城鎮距離的增加，農村居民點斑塊面積和斑塊密度整體趨于減少趨勢。斑塊密度處于增大趨勢，表明城鎮的吸引及輻射作用，使得農村居民點向著城鎮周邊聚集分布，越靠近城鎮，規模越大；而斑塊個數則隨著與城鎮距離的增加而減少，由于城鎮周圍分布的農村居民點規模較大，連續性較高，其斑塊個數則比與城鎮距離較遠的農村居民點個數少。斑塊形狀和聚集度隨距離的增加而減小，表明農村居民點在城鎮周邊處于聚集規模分布狀態，但仍缺乏整體的規劃和管理。整體來看，城鎮對于農村居民點的吸引力隨著距離的增加而衰減。同時由表9可得，1980—2015年研究期內，斑塊面積和聚合度隨時間的推移而增長，但斑塊個數呈衰減趨勢，可見城鎮對農村居民點的吸引力隨時間的增長而增強，且周邊聚落的集聚更加明顯。

表10 1980—2015年青藏高原地區不同距離城鎮緩沖區農村居民點景觀特征Table 10 Landscape characteristics of villages and towns in different distance buffer zones in Qinghai-Tibet region from 1980 to 2015

綜上，區位因素是一種外部不可控因素，對農村居民點有強烈的作用力，是農村居民點選址、規模及形態變化最有活力的驅動力之一。區位因素具有動態性[10]，當城鎮、公路等區位因素建成后，農村居民點向其周邊集聚，但已確定的農村居民點也會對道路這樣的人為區位因素產生作用。區位因素對農村居民點的影響具有隨距離增加而衰減的特征，作用范圍一般也很小，只對鄰近區域有較大的影響力，但在自然區位優勢較明顯的范圍內，農村居民點的聚集往往比較隨意，缺乏規劃。

4.2.2 人文因素 農村居民點發展的人文因素主要體現在人口、經濟及政策三方面，在三者相互作用之下對農村居民點的演化起到明顯的作用。

人口因素作為農村居民點發展演化的內部驅動因素，農村居民點建設首先要滿足人口居住，人口增加的直觀展現就是農村居民點的面積增加。通過計算青藏高原地區人口數量增加值及農村居民點面積增加數量相關關系，相關系數R2=0.834 4，人口和農村居民點斑塊呈顯著正相關，故人口數量是影響農村居民點規模與形態的重要因素。

經濟發展作為農村居民點演變的外部驅動因素，經濟發展好，可以為農村居民點的擴展提供可靠的外部基礎，為人才引進及新農村建設發展提供有力的支持。青藏高原地區作為高海拔地區，農村經濟主要以河谷地帶的農業耕作及高寒草原的畜牧業為主，選取農林牧漁總產值變化與農村居民點面積變化做相關關系計算。相關系數R2=0.551 1，可見農村經濟和農村居民點斑塊呈正相關，但其相關性并不算高，究其原因是青藏高原地區地形特殊，其社會經濟受自然條件的約束明顯，農村經濟對農村居民點有一定的影響，但不是決定性因素。因此，適宜農耕的河谷地帶及適宜放牧的高寒草原地帶也影響著農村居民點的分布和規模。

政策因素是農村居民點演化的綜合控制力。政策制度可以直接對農村居民點的規模及選址起到規范作用。通過啟動“千村建設、百村示范”，建設高原美麗宜居村鎮計劃，進行青藏高原地區環境治理、土地整治及生態修復等措施使得農村居民點內部系統更加合理；鄉村規劃、遷村并鎮及富民安居等一系列政策，也在時刻影響著農村居民點形態、規模等特征。

綜上，人文因素是一種集人口、經濟和政策影響的復合因素，對農村居民點的演變具有時效性，動態性和階段性等特征。其作用范圍往往受到行政區劃的影響，但其對農村居民點的分布、形態及規模往往具有最重要的作用。總體來說，因青藏高原地區地形獨特，其農村居民點布局受自然條件限制和人文因素的影響。

5 結論與對策建議

5.1 結論

本文以地形復雜的青藏高原地區為研究區域，利用ArcGIS的空間分析和景觀格局指數等方法，定量地分析了青藏高原地區農村居民點用地分布特征及規模分異特征，并對其驅動機制進行了探討，研究發現：

1）青藏高原地區農村居民點用地擴張明顯，新增的農村居民點用地主要來源于對耕地、未利用地和草地的侵占；

2）在研究期內，農村居民點斑塊破碎度降低，逐漸趨于連續性，主要以青海東部、南部以及西藏南部為集聚中心，高值集聚明顯；

3）青藏高原地區農村居民點布局受自然因素和區位因素的影響顯著，且自然因素對農村居民點的發展起約束性作用。

5.2 對策建議

通過研究發現，青藏高原地區農村居民點用地擴張方式多樣，且居民點斑塊呈不規則形式，農村居民點缺乏規劃先行的思維，農村居民點主要是以西寧和拉薩等經濟發達地區為集聚中心，新增居民點由中心向外擴張，但在政策的扶持和交通的發展下，原本偏遠的那曲市、阿里市等地區也出現了農村居民點的集聚中心，但還未形成大規模集聚，對交通的依賴性也較強。由此可見，在未來青藏高原地區的農村居民點發展過程中，應注意以下兩點：

1）合理規劃農村居民點。中央第七次西藏工作會議中提出，要將生態保護與民生改善相結合，青藏高原的生態保護是地方發展的重中之重。因此在進行農村居民點建設時，要進行合理規劃，嚴格把握對各類用地的占用，保證占補平衡，嚴格控制生態保護紅線。

2）逐步完善基礎設施，推動鄉村高質量發展。依托各類援藏政策，實現脫貧攻堅和鄉村振興的有效銜接。在鄉村建設中對于已形成的村落，可通過對基礎設施的優化、完善，滿足居民日常生活需求，使居民生活便利化、現代化、多樣化。與此同時，通過大力發展農村經濟及交通網絡，形成強有力的農村居民點集聚中心，帶動周圍居住環境質量較差的居民向生活條件相對優越的地區轉移，實現生態保護和鄉村振興協同發展。

本文在研究過程中還存在以下不足：①在驅動因素研究過程中，由于數據不全，主要采取定性的方式對社會經濟因素進行分析，缺乏定量化的研究，在后續的研究中將繼續收集相關數據，進一步對農村居民點演化的驅動因素進行定量化的研究。②僅利用三期數據對比了青藏高原地區農村居民點的空間演化特征，農村居民點的演化是復雜的，后續的研究中還可以通過更多期數據的對比，進而揭示青藏高原地區農村居民點的演化規律。在后續的研究中將針對以上兩點不足，進一步展開研究，以期為青藏高原地區農村居民點的合理規劃及人居環境的改善提供參考借鑒。