桂林市農(nóng)村居民點(diǎn)變化特征及驅(qū)動(dòng)力分析

李麗清

（廣西國(guó)土資源規(guī)劃設(shè)計(jì)集團(tuán)有限公司，廣西 桂林 541000）

農(nóng)村居民點(diǎn)作為農(nóng)村人口的“歸宿地”，是農(nóng)村地區(qū)重要的生活、生產(chǎn)和生態(tài)空間[1]，一直以來(lái)，我國(guó)農(nóng)村居民點(diǎn)建房多是農(nóng)民自主選擇，盲目自建、缺乏統(tǒng)一規(guī)劃是普遍存在的問(wèn)題。隨著新型城鎮(zhèn)化、工業(yè)化進(jìn)程的加快，農(nóng)村居民點(diǎn)用地散亂、無(wú)序、空心化問(wèn)題愈發(fā)突出，極大地阻礙了城鄉(xiāng)一體化的推進(jìn)。鑒于此，合理規(guī)劃農(nóng)村居民點(diǎn)的用地規(guī)模和空間布局顯得尤為重要，而把握農(nóng)村居民點(diǎn)長(zhǎng)期演化規(guī)律是合理規(guī)劃村莊布局的基礎(chǔ)。

目前，針對(duì)農(nóng)村居民點(diǎn)空間布局特征及演變的研究較為廣泛，研究結(jié)果主要集中在空間分布特征[2-4]、布局優(yōu)化及重構(gòu)[5-7]、整治潛力測(cè)算[8-10]、時(shí)空演化及驅(qū)動(dòng)機(jī)制[11-13]、聚落景觀與類(lèi)型劃分[14-16]；研究方法多采用核密度分析[17-18]、重心遷移[19]、空間自相關(guān)[20]、熱點(diǎn)分析[21-22]、景觀格局指數(shù)[23-25]；研究區(qū)域涵蓋了平原、山區(qū)、丘陵等；研究尺度從全國(guó)、省的宏觀尺度到鎮(zhèn)、村的微觀尺度均有涉及。但總體來(lái)看，已有研究數(shù)據(jù)多源于遙感影像提取，利用地面實(shí)際調(diào)查數(shù)據(jù)的研究較少。本文借助桂林市10年來(lái)土地利用現(xiàn)狀數(shù)據(jù)，應(yīng)用統(tǒng)計(jì)分析法、空間分析法、核密度計(jì)算法等研究方法，從圖斑角度分析2009—2018年農(nóng)村居民點(diǎn)變化特征及演變規(guī)律，并借助地理探測(cè)器探測(cè)影響農(nóng)村居民點(diǎn)規(guī)模和空間格局變化的驅(qū)動(dòng)因素，為桂林市農(nóng)村居民點(diǎn)動(dòng)態(tài)演變監(jiān)測(cè)和空間布局優(yōu)化提供理論依據(jù)，促進(jìn)城鄉(xiāng)土地資源優(yōu)化配置。

1 材料與研究方法

1.1 研究區(qū)概況

桂林市位于廣西東北部，全市現(xiàn)轄6區(qū)10縣1市。市內(nèi)地形復(fù)雜，以山地地貌為主，丘陵、盆地、巖溶地貌兼而有之，巖溶（喀斯特）地貌發(fā)育完美，景觀奇特；江河縱橫、溪流遍布，地表水和地下水資源豐富，主要河流有湘江、資江、桂江（漓江）、洛清江、尋江（古宜河）等5條主干河；氣候?qū)僦衼啛釒Ъ撅L(fēng)氣候，氣候溫和，光照充足，無(wú)霜期長(zhǎng)，水熱條件優(yōu)越。桂林屬以旅游和農(nóng)林為重點(diǎn)的廣西桂北經(jīng)濟(jì)區(qū)。

1.2 數(shù)據(jù)來(lái)源

2009年、2012年、2015年、2018年四期土地利用現(xiàn)狀數(shù)據(jù)來(lái)自各年度土地利用變更調(diào)查成果，DEM數(shù)據(jù)來(lái)源于地理數(shù)據(jù)云網(wǎng)站，分辨率30 m×30 m。人口、GDP等數(shù)據(jù)來(lái)源于桂林市歷年國(guó)民經(jīng)濟(jì)與社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)。

1.3 研究方法

1.3.1 空間分析法

借助ArcGIS10.3軟件，統(tǒng)計(jì)研究時(shí)段內(nèi)居民點(diǎn)變化規(guī)模和格局特征，利用緩存區(qū)分析和相交分析得到不同距離范圍內(nèi)居民點(diǎn)規(guī)模特點(diǎn)。

1.3.2 核密度計(jì)算法

核密度計(jì)算法是通過(guò)樣本數(shù)據(jù)來(lái)研究數(shù)據(jù)整體分布特征的一種方法，常用來(lái)反映點(diǎn)狀要素的空間分布特征，是熱點(diǎn)區(qū)和冷區(qū)識(shí)別和分析的一種探測(cè)性方法[26]。其原理是通過(guò)離散點(diǎn)或線數(shù)據(jù)進(jìn)行內(nèi)插，落入搜索區(qū)的點(diǎn)具有不同的權(quán)重，越靠近中心的點(diǎn)其權(quán)重越大，反之，權(quán)重較小。本文借助ArcGIS10.3進(jìn)行核密度計(jì)算，公式如下[27]：

式中：f（x，y）為位于（x，y）位置的密度估計(jì)；h＞0為帶寬；n為觀測(cè)數(shù)量；K為核函數(shù)；di為位置距第i個(gè)觀測(cè)位置的距離。

1.3.3 利用地理探測(cè)器

地理探測(cè)器是一種新的統(tǒng)計(jì)學(xué)工具，可以探測(cè)空間分異性并揭示其背后的驅(qū)動(dòng)因子[28]。本文利用地理探測(cè)器對(duì)農(nóng)村居民點(diǎn)規(guī)模和空間布局演化模式的影響機(jī)理進(jìn)行分析，計(jì)算公式如下：

式中，PD，U為農(nóng)村居民點(diǎn)規(guī)模影響因素D探測(cè)力值；nD，i為次一級(jí)區(qū)域樣本數(shù)；n為整個(gè)區(qū)域樣本數(shù)；m為次級(jí)區(qū)域個(gè)數(shù)；為整個(gè)區(qū)域農(nóng)村居民點(diǎn)規(guī)模的方差，為次一級(jí)區(qū)域的方差。假設(shè)模型成立。PD，U的取值范圍為 [0，1]，PD，U值越大，說(shuō)明D因素對(duì)縣域農(nóng)村居民點(diǎn)用地面積的影響程度越高。

2 結(jié)果與分析

2.1 農(nóng)村居民點(diǎn)規(guī)模和空間格局演變特征分析

2.1.1 規(guī)模及變化特征分析

2009—2018 年，桂林市居民點(diǎn)數(shù)量和用地面積變化表現(xiàn)出逐年遞增的一致性（表1），但平均居民點(diǎn)規(guī)模逐漸降低，從0.91hm2降至0.79hm2，反映了桂林市農(nóng)村居民點(diǎn)呈破碎化增長(zhǎng)趨勢(shì)。

表1 桂林市2009—2018年農(nóng)村居民點(diǎn)用地?cái)?shù)量和面積

從桂林市農(nóng)村居民點(diǎn)用地規(guī)模頻率分布來(lái)看（圖1），頻率呈偏態(tài)分布，數(shù)據(jù)嚴(yán)重偏離中心性。2009—2018年偏度均超過(guò)7，中值和均值均逐年降低，且中值均大于均值，說(shuō)明桂林市農(nóng)村居民點(diǎn)規(guī)模存在大的低值聚集分布和高值離散分布，且10年來(lái)農(nóng)村居民點(diǎn)離散度沒(méi)有改善。

圖1 桂林市農(nóng)村居民點(diǎn)用地規(guī)模頻率分布圖

2.1.2 空間分布及演變特征分析

分析農(nóng)村居民點(diǎn)斑塊的核密度圖（圖2）可以看出，2009—2018年核密度圖分布格局相似，為由東北向南、南東方向的帶狀分布，大致與湘桂鐵路走向一致。2018年農(nóng)村居民點(diǎn)核密度圖與2009年相比，2個(gè)區(qū)域有較為明顯的擴(kuò)大，分別為桂林老城區(qū)和臨桂新區(qū)，表明10年間此兩個(gè)區(qū)域居民點(diǎn)數(shù)量在單位面積內(nèi)增加明顯；灌陽(yáng)縣洞井瑤族鄉(xiāng)和觀音閣鄉(xiāng)核密度有較為明顯的降低，原因?yàn)?015年左右實(shí)施較大面積增減掛鉤拆舊項(xiàng)目。

圖2 2009年和2018年桂林市農(nóng)村居民點(diǎn)核密度分布圖

2.2 基于地理探測(cè)器的空間格局變化驅(qū)動(dòng)因素分析

2.2.1 農(nóng)村居民點(diǎn)格局變化的驅(qū)動(dòng)因素分析

農(nóng)村居民點(diǎn)格局演變受自然環(huán)境本底及經(jīng)濟(jì)社會(huì)發(fā)展等多重因素的影響，其本質(zhì)是農(nóng)村生產(chǎn)、生活與周?chē)h(huán)境因素趨向統(tǒng)一。為探索桂林市農(nóng)村居民點(diǎn)空間格局演變的驅(qū)動(dòng)機(jī)制，綜合考慮指標(biāo)代表性、數(shù)據(jù)可獲取性等因素，選取了海拔（X1）、坡度（X2）、交通可達(dá)性（X3）、鄉(xiāng)鎮(zhèn)中心可達(dá)性（X4）、耕作便捷度（X5）、縣城中心可達(dá)性（X6）、工礦企業(yè)可達(dá)性（X7）、距水源距離（X8）等8項(xiàng)指標(biāo)作為驅(qū)動(dòng)因子，探測(cè)影響桂林市農(nóng)村居民點(diǎn)空間格局變化的驅(qū)動(dòng)因素，q值越大，驅(qū)動(dòng)程度越強(qiáng)，反之亦然（表2）。

表2 桂林市農(nóng)村居民點(diǎn)空間格局驅(qū)動(dòng)因素探測(cè)結(jié)果

從計(jì)算結(jié)果來(lái)看，影響農(nóng)村居民點(diǎn)空間格局變化的核心因子是工礦企業(yè)可達(dá)性、縣城中心可達(dá)性及海拔。從10年來(lái)新增農(nóng)村居民點(diǎn)位置來(lái)看，農(nóng)村居民點(diǎn)的布局變化受工礦企業(yè)分布、縣城位置及道路走向影響顯著，具有明顯的指向性。從減少農(nóng)村居民點(diǎn)分布來(lái)看，農(nóng)村居民點(diǎn)布局變化受縣城位置、海拔影響較大。

2.2.2 道路對(duì)農(nóng)村居民點(diǎn)分布的指向性

近年來(lái)，越來(lái)越多的農(nóng)村居民點(diǎn)向道路靠攏，具體表現(xiàn)在：①?gòu)淖兓娣e上看，距離道路越近，居民點(diǎn)面積凈增加越大；距離道路越遠(yuǎn)，居民點(diǎn)面積凈減少越大。②從占比上看，0～1000m范圍內(nèi)是居民點(diǎn)增加的集中區(qū)域，占增加居民點(diǎn)總數(shù)的88.75%；同時(shí)也是居民點(diǎn)減少的主要區(qū)域，占減少居民點(diǎn)總面積的73.40%（表3）。隨著與道路距離的增加，居民點(diǎn)變化活躍度急劇減少。

表3 道路各級(jí)緩存區(qū)內(nèi)農(nóng)村居民點(diǎn)變化

2.2.3 縣城對(duì)農(nóng)村居民點(diǎn)分布的輻射性

以縣城為中心向周?chē)椛?，隨著與縣城距離的增加，居民點(diǎn)變化的活躍度有所降低，增加居民點(diǎn)面積趨于減少（表4），說(shuō)明農(nóng)村居民點(diǎn)分布與縣城的距離有較大關(guān)系。8km范圍內(nèi)是居民點(diǎn)減少的主要區(qū)域，其中5km范圍內(nèi)減少的居民點(diǎn)中有87.38%變?yōu)槌擎?zhèn)住宅，12.51%變?yōu)榻煌ㄋ捌渌擎?zhèn)建設(shè)用地，反映了縣城周邊的城鎮(zhèn)化現(xiàn)象。

表4 縣城各級(jí)緩存區(qū)內(nèi)農(nóng)村居民點(diǎn)變化

2.2.4 工礦企業(yè)對(duì)農(nóng)村居民點(diǎn)分布的輻射性

隨著與工礦企業(yè)用地距離的增加，增加和減少居民點(diǎn)面積均急劇減少（表5）。距工礦企業(yè)2000m范圍內(nèi)是農(nóng)村居民點(diǎn)變化的主要區(qū)域，1000m范圍內(nèi)增加居民點(diǎn)面積占總增加面積的77.70%，表明工礦企業(yè)對(duì)居民點(diǎn)分布有較大的吸引力。2000m范圍內(nèi)減少居民點(diǎn)面積占總減少面積的72.74%，其中40%的居民點(diǎn)變更為城鎮(zhèn)住宅，21.09%的居民點(diǎn)變更為交通水利及其他城鎮(zhèn)建設(shè)用地，在一定程度上反映出工礦企業(yè)拉動(dòng)了城鎮(zhèn)化進(jìn)程。

表5 工礦企業(yè)各級(jí)緩存區(qū)內(nèi)農(nóng)村居民點(diǎn)變化

2.2.5 海拔對(duì)農(nóng)村居民點(diǎn)分布的約束性

由表6可知，10年來(lái)增加的農(nóng)村居民點(diǎn)中有86.77%位于海拔小于或等于350m范圍內(nèi)，且隨著海拔升高，增加的農(nóng)村居民點(diǎn)急劇減少。從居民點(diǎn)規(guī)模變化來(lái)看，10年來(lái)海拔在350m以下區(qū)域內(nèi)，增加農(nóng)村居民點(diǎn)規(guī)模是減少居民點(diǎn)規(guī)模的3.3倍；海拔大于350m區(qū)域內(nèi)，減少農(nóng)村居民點(diǎn)規(guī)模是增加居民點(diǎn)規(guī)模的1.7倍，說(shuō)明農(nóng)村居民點(diǎn)有向低海拔區(qū)集聚的趨勢(shì)。

表6 不同海拔緩存區(qū)內(nèi)農(nóng)村居民點(diǎn)變化

2.3 基于地理探測(cè)器的用地規(guī)模變化驅(qū)動(dòng)因素分析

為更好地了解影響農(nóng)村居民點(diǎn)用地規(guī)模變化的因素，綜合考慮指標(biāo)代表性、數(shù)據(jù)可獲取性等因素，選取GDP（X1）、全社會(huì)固定資產(chǎn)投資（X2）、第一產(chǎn)業(yè)增加值（X3）、農(nóng)村居民人均純收入（X4）、產(chǎn)業(yè)非農(nóng)業(yè)化率（X5）、人均糧食產(chǎn)量（X6）、單位面積糧食產(chǎn)量（X7）、耕地面積（X8）、園地面積（X9）、城鎮(zhèn)化率（X10）、人口（X11）等11項(xiàng)指標(biāo)作為驅(qū)動(dòng)因子，探測(cè)影響桂林市農(nóng)村居民點(diǎn)用地規(guī)模變化的驅(qū)動(dòng)因素。探測(cè)結(jié)果表明，農(nóng)村居民點(diǎn)用地面積規(guī)模主要受到第一產(chǎn)業(yè)增加值、城鎮(zhèn)化率、農(nóng)村居民人均純收入、GDP因子的影響；其次為全社會(huì)固定資產(chǎn)投資、人口及園地面積；耕地面積、人均糧食產(chǎn)量、單位面積糧食產(chǎn)量及產(chǎn)業(yè)非農(nóng)業(yè)化率對(duì)農(nóng)村居民點(diǎn)用地規(guī)模變化影響均較小（表7）。

表7 桂林市農(nóng)村居民點(diǎn)用地規(guī)模擴(kuò)張強(qiáng)度影響因子探測(cè)結(jié)果表

農(nóng)村居民點(diǎn)面積變化與GDP、人口、農(nóng)村居民人均純收入、全社會(huì)固定資產(chǎn)投資、第一產(chǎn)業(yè)增加值、城鎮(zhèn)化率呈正相關(guān)性，與耕地?zé)o明顯的正相關(guān)性，與園地有較為明顯的負(fù)相關(guān)性（圖3）。反映了隨著農(nóng)村居民收入途徑的多樣化，務(wù)農(nóng)收入在農(nóng)民收入中不再占據(jù)主導(dǎo)地位。隨著農(nóng)村居民人均純收入水平不斷提高，農(nóng)村居民對(duì)建房的需求變得更加迫切，而這一過(guò)程又加大了對(duì)固定資產(chǎn)的投入。

圖3 農(nóng)村居民點(diǎn)面積變化與各影響因素之間的相關(guān)性

3 結(jié)論

（1）從整體上看，2009—2018年，桂林市農(nóng)村居民點(diǎn)用地規(guī)模擴(kuò)張明顯，2018年農(nóng)村居民點(diǎn)用地規(guī)模比2009年增長(zhǎng)了1.51%，平均居民點(diǎn)規(guī)模從0.91hm2降至0.79hm2，表明桂林市農(nóng)村居民點(diǎn)聚集度不夠，且無(wú)明顯改善。

（2）農(nóng)村居民點(diǎn)規(guī)模頻率分布圖表明，桂林市農(nóng)村居民點(diǎn)分布呈現(xiàn)大規(guī)模低密度和小規(guī)模高密度的空間關(guān)聯(lián)特征。2009年和2018年的核密度分布特征基本一致，10年來(lái)，除桂林市城區(qū)的核密度有較為明顯增加，灌陽(yáng)縣洞井瑤族鄉(xiāng)和觀音閣鄉(xiāng)核密度有較為明顯的降低外，其余地區(qū)農(nóng)村居民點(diǎn)分布特征變化并不明顯。

（3）桂林市農(nóng)村居民點(diǎn)空間格局的變化受多重因素的影響和制約。就單因素而言，地勢(shì)低且平緩的區(qū)域是農(nóng)村居民點(diǎn)選址的自然基礎(chǔ)，道路兩側(cè)1km范圍內(nèi)、工礦企業(yè)2km范圍內(nèi)對(duì)農(nóng)村居民點(diǎn)具有強(qiáng)力的吸引作用和衰減特征，縣城對(duì)農(nóng)村居民點(diǎn)分布的影響作用存在顯著的輻射局限性。各因素作用強(qiáng)度依次為工礦企業(yè)可達(dá)性、縣城中心可達(dá)性、海拔、交通可達(dá)性、距水源距離、鄉(xiāng)鎮(zhèn)中心可達(dá)性、耕作便捷度、坡度等。

（4）桂林市農(nóng)村居民點(diǎn)用地規(guī)模變化受多種因素的影響，可歸納為農(nóng)業(yè)發(fā)展和農(nóng)民生活水平、縣域經(jīng)濟(jì)發(fā)展和投資帶動(dòng)、人口增長(zhǎng)和城鎮(zhèn)化率等3類(lèi)驅(qū)動(dòng)力，其中農(nóng)業(yè)發(fā)展和農(nóng)民生活水平、人口增長(zhǎng)和城鎮(zhèn)化率是農(nóng)村居民點(diǎn)規(guī)模增加的內(nèi)在驅(qū)動(dòng)力，是首要驅(qū)動(dòng)力，也是直接驅(qū)動(dòng)因素。縣域經(jīng)濟(jì)發(fā)展和投資帶動(dòng)是推動(dòng)農(nóng)村居民點(diǎn)擴(kuò)張與收縮的外在動(dòng)力，是重要驅(qū)動(dòng)因素，也是間接推動(dòng)力。各因素作用強(qiáng)度依次為第一產(chǎn)業(yè)增加值、城鎮(zhèn)化率、農(nóng)村居民人均純收入、GDP、全社會(huì)固定資產(chǎn)投資、人口及園地面積、耕地面積、人均糧食產(chǎn)量、單位面積糧食產(chǎn)量、產(chǎn)業(yè)非農(nóng)業(yè)化率。

4 啟示

農(nóng)村居民點(diǎn)是農(nóng)村地域重要的生產(chǎn)、生活、生態(tài)空間，其形態(tài)格局是自然、經(jīng)濟(jì)、社會(huì)多種因素影響下的產(chǎn)物。在快速城鎮(zhèn)化進(jìn)程中，農(nóng)村居民點(diǎn)布局發(fā)生了劇烈的變革，在促進(jìn)城鄉(xiāng)融合的同時(shí)也帶來(lái)了村莊建設(shè)粗放無(wú)序、農(nóng)村空心化等一系列問(wèn)題。本文通過(guò)分析10年來(lái)農(nóng)村居民點(diǎn)規(guī)模及空間格局演變特征，探測(cè)影響農(nóng)村居民點(diǎn)變化的驅(qū)動(dòng)因素，揭示了農(nóng)村居民點(diǎn)發(fā)展演變趨勢(shì)，在一定程度上為村莊規(guī)劃指引方向。

在鄉(xiāng)村振興戰(zhàn)略背景下，農(nóng)村居民點(diǎn)布局散亂無(wú)序、用地粗放低效、亂占耕地的現(xiàn)象亟須轉(zhuǎn)變?；诒疚姆治龊徒Y(jié)論，得出以下啟示。

（1）劃定村莊擴(kuò)張邊界，合理引導(dǎo)村莊建房。理論上，隨著城鎮(zhèn)化帶動(dòng)農(nóng)村人口非農(nóng)化速度加快，農(nóng)村居民點(diǎn)用地需求會(huì)相對(duì)減少。然而多年來(lái)，農(nóng)村居民點(diǎn)用地規(guī)模不降反升。一方面是農(nóng)村居民進(jìn)入城市購(gòu)置房屋以后，受傳統(tǒng)思想的影響，原村莊用房依舊保留；另一方面，外出務(wù)工者多數(shù)無(wú)力在工作地購(gòu)房，務(wù)工收入回流農(nóng)村，進(jìn)一步助推了建房熱潮，造成“一戶(hù)多宅”“亂占耕地”現(xiàn)象并存。長(zhǎng)期以來(lái)，我國(guó)建設(shè)規(guī)劃“重城市輕農(nóng)村”，村莊規(guī)劃缺位或規(guī)劃不到位是形成當(dāng)前狀況的根本原因。因此，為防止村莊建房無(wú)序蔓延，需結(jié)合村莊地理?xiàng)l件，把握村莊發(fā)展演變規(guī)律，合理確定村莊發(fā)展方向及規(guī)模，有序引導(dǎo)村民建房選址。

（2）完善政策法規(guī)，堅(jiān)決遏制農(nóng)村亂占耕地建房。制定符合實(shí)際情況的整治策略，為開(kāi)展農(nóng)村亂占耕地建房問(wèn)題專(zhuān)項(xiàng)整治工作提供參考依據(jù)，確保有效解決農(nóng)村亂占耕地建房難題。

（3）盤(pán)活農(nóng)村閑置土地，提高土地集約節(jié)約利用。在鄉(xiāng)村振興戰(zhàn)略發(fā)展背景下，城鎮(zhèn)用地規(guī)模有限和城鎮(zhèn)人口增加、農(nóng)村人口減少與村莊用地增加之間的矛盾越顯突出，多元化利用好農(nóng)村閑置宅基地對(duì)于緩解用地矛盾具有重要意義，也是融合發(fā)展一、二、三產(chǎn)業(yè)，升級(jí)農(nóng)村經(jīng)濟(jì)的有力抓手。宅基地盤(pán)活要因地制宜，要充分結(jié)合地理區(qū)位條件，做到宜耕則耕、宜農(nóng)則農(nóng)、宜游則游、宜水則水、宜業(yè)則業(yè)。