肖建能，丁華祥，趙耀龍，溫宥越

1. 廣東省國土資源測繪院，廣州 510663；

2. 自然資源部華南熱帶亞熱帶自然資源監(jiān)測重點實驗室，廣州 510663；

3. 華南師范大學 地理科學學院 廣東省智慧國土工程技術研究中心， 廣州 510631；

4. 生態(tài)環(huán)境部華南環(huán)境科學研究所， 廣州 510535

1 引 言

農(nóng)村宅基地發(fā)展是自然和人文等多種因素綜合作用的結果，宅基地的空間形態(tài)變化反映了社會經(jīng)濟發(fā)展的痕跡，也是推進鄉(xiāng)村振興過程重點研究方向（陳榮清等，2010；楊文杰等，2023）。研究農(nóng)村宅基地時空演變規(guī)律，有助于正確把握和理解其在城市化進程中所展現(xiàn)的特征、發(fā)展趨勢，為農(nóng)村未來發(fā)展預測提供技術支撐，有利于促進農(nóng)村建設用地合理利用。長期以來，眾多學者在農(nóng)村宅基地利用管理及制度建設方面進行了大量研究，其主要研究歐盟地區(qū)最先進和最落后的18 個農(nóng)村地區(qū)揭示文化與農(nóng)村發(fā)展的關系（Terluin，2003）、韓國農(nóng)村農(nóng)業(yè)變革與旅游發(fā)展（Kim 和Jamal，2015）、國內政府經(jīng)濟發(fā)展對農(nóng)村轉型促進作用（Long 等，2011）、晉城市村莊空心化原因及相關治理方案（原野等，2015）、山東禹城市空心村整治的基本準則（劉彥隨等，2009）、河南鄲城縣重點農(nóng)村空心化特征形成機制和調控策略（崔衛(wèi)國等，2011）。

在農(nóng)村宅基地空間變化分析方面，也有學者開展研究，包括全國（王磊和李哲敏，2021）、西南地區(qū)的重慶和云南兩省（市）三個農(nóng)村（蘇康傳，2020）、環(huán)渤海地區(qū)農(nóng)村（Yang 等，2015）、貴州惠水縣（楊智謀，2021）、江西貴溪市（張晗等，2019）等地區(qū)的農(nóng)村宅基地分析。作為改革開放先行省份的廣東，經(jīng)過40 多年的發(fā)展，農(nóng)村的發(fā)展及其宅基地發(fā)生巨大變化。但現(xiàn)有研究對廣東省的農(nóng)村宅基地時空演化特性分析較少，其演變背后的驅動力不清楚，因此，開展典型農(nóng)村宅基地的空間演變規(guī)律研究，對廣東省農(nóng)村發(fā)展規(guī)劃具有重大意義。

農(nóng)村宅基地演變與發(fā)展，受到人文（姜廣輝等，2007；鄭先迪等，2019）、自然（師滿江等，2016；董光龍等，2017；王涵和趙翠薇，2020）、社會（Tian等，2014；劉同山，2018；胡斯斌等，2021）、經(jīng)濟（徐濤等，2022）、政策（陳陽等，2014；孔雪松等，2014；趙路等，2022）等多方面因素的影響。因此，研究農(nóng)村，遵循農(nóng)村發(fā)展的客觀規(guī)律，把握農(nóng)村的系統(tǒng)要素空間分異特征尤為重要（劉彥隨等，2019）。為探究粵西沿海地區(qū)的農(nóng)村宅基演變特征，本文以茂名市電白區(qū)為例，利用2017 年廣東省農(nóng)村地籍調查成果、第一次全國土地調查（以下簡稱“一調”）數(shù)據(jù)、第二次全國土地調查（以下簡稱“二調”）數(shù)據(jù)等已有的農(nóng)村土地調查數(shù)據(jù)，分析農(nóng)村宅基地空間演變特征，進而探究宅基地演變的主要驅動力，以期為粵西沿海地區(qū)的村莊土地利用規(guī)劃提供技術支撐。

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

2.1 研究區(qū)概況

電白區(qū)是茂名市的轄區(qū)，位于粵西沿海地區(qū)，北靠高州市、陽春市，東部交接陽西縣，西接茂南區(qū)，南部瀕臨中國南海。電白區(qū)北部屬中低山地，中部屬沿江平原和低丘陵地，西南部為黃土丘陵，南部屬沿海臺地，海岸線總長約220 km。電白區(qū)自然條件優(yōu)越，農(nóng)業(yè)資源豐富，交通便利，具有獨特農(nóng)產(chǎn)品和生產(chǎn)優(yōu)勢，發(fā)展農(nóng)業(yè)有著得天獨厚的優(yōu)越，電白區(qū)也在發(fā)揮沿海、平原、山區(qū)“三位一體”的地理優(yōu)勢，不斷推進農(nóng)業(yè)產(chǎn)業(yè)化，全區(qū)的農(nóng)業(yè)農(nóng)村經(jīng)濟在粵西沿海地區(qū)較為靠前（張志勇，2017）。因此，本文以電白區(qū)作為粵西沿海地區(qū)典型區(qū)域，開展農(nóng)村宅基地空間演變特征研究。

2.2 數(shù)據(jù)來源

電白區(qū)的2017 年農(nóng)村地籍調查成果和地理國情監(jiān)測數(shù)據(jù)來源于廣東省自然資源廳；電白區(qū)“一調”土地利用現(xiàn)狀圖、“二調”數(shù)據(jù)、1∶1 萬數(shù)字高程模型（digital elevation model，DEM）數(shù)據(jù)來源于廣東省國土資源檔案館。

2.3 研究方法

2.3.1 農(nóng)村宅基地數(shù)據(jù)處理

首先，將“二調”數(shù)據(jù)坐標系轉換為2000 年國家大地坐標系；其次，以“二調”的影像作為底圖，以“二調”地類圖斑中屬性是203 類作為“二調”農(nóng)村的范圍；再次，用2017 年的農(nóng)村地籍調查成果與“二調”的村莊范圍進行套合分析，結合“二調”影像，將超出村莊范圍外，以及不是房屋的2017 年原宅基地進行更新；最后，提取出“二調”時期農(nóng)村宅基地數(shù)據(jù)。

由于“一調”的土地利用現(xiàn)狀圖數(shù)據(jù)僅有1∶1萬比例尺的紙質圖集，首先，用1∶1 萬標準圖幅框，根據(jù)圖幅編號進行匹配，將紙質圖集轉換成具有空間坐標的數(shù)據(jù)集；其次，用“二調”的農(nóng)村宅基地數(shù)據(jù)與矯正后“一調”土地利用數(shù)據(jù)套合，對超出“一調”村莊范圍外的原“二調”宅基地進行更新；最后，提取出“一調”時期的農(nóng)村宅基地數(shù)據(jù)。

2.3.2 驅動力因子數(shù)據(jù)處理

本文選取地形復雜度（Ter）、人口密度（PopD）、路網(wǎng)密度（RodDen）、到行政中心距離（D）、親水性（RIVDen）、聚耕比（Ag）共六個因子，采用相關性分析方法，分析影響電白區(qū)農(nóng)村宅基地時空演變的主要驅動力。相關方法及各因子計算及說明如下。

1）研究方法

通過Pearson 相關性分析方法，分析“一調”至2017 年，電白區(qū)農(nóng)村宅基地占地面積的演變驅動力因子，Pearson 相關系數(shù)計算原理如下（Eggers等，2003）：

式中，r為要素x與y之間的相關系數(shù)，取值為[–1，1]。r>0 表示正相關，r<0 表示負相關；r的絕對值越大，變量之間的相關性越強，即相關系數(shù)越接近1 或—1 相關性越強，反之，相關系數(shù)r越接近0，變量之間相關性就越弱。

2）增幅

選取增幅指標判定農(nóng)村宅基地的變化情況，增幅的計算方式如下：

3）地形復雜度（Ter）

基于DEM 數(shù)據(jù)，利用ArcGIS 分別計算地形起伏度（?）、坡度（Δ）、坡向（θ），將三者進行加權計算，得到地形復雜度（Ter）。這用以分析典型縣區(qū)農(nóng)村宅基地與地形復雜度之間的關系。

在計算出地形起伏度、坡度、坡向后，分別進行歸一化處理，并由下式計算出地形復雜度：

4）路網(wǎng)密度（RodDen）

本文選取電白區(qū)的高速路、國道、省道、縣道、城市道路、匝道、鄉(xiāng)道和二級鄉(xiāng)道，根據(jù)道路等級，研究設計了不同的道路等級權重值，具體如表1 所示，村級的路網(wǎng)密度：

表1 道路等級權重設置Tab.1 Weight setting of road grade

式中，RodDen為路網(wǎng)密度，km/km2；RodL為路網(wǎng)總長度，km；A為村域面積，km2。

5）親水性（RivDen）

親水性是自古以來村鎮(zhèn)選擇區(qū)位時要考慮的重要因素之一，已有研究證明，區(qū)域河流的密度對農(nóng)村宅基地的分布有著顯著的影響（任平等，2014）。本文從2017 年電白區(qū)地理國情監(jiān)測數(shù)據(jù)中，提取出一級至五級的河流，對不同河流級別分別賦予不同的河流等效長度比值（表2）：

表2 河流等效長度對照表Tab.2 Comparison table of river equivalent length

表 3 電白區(qū)宅基地面積增幅與驅動力因子相關性分析Tab.3 Correlation analysis of homestead growth and driving factors in Dianbai District

式中，RivDen為親水性（河流密度），km/km2；RivL為河流等效總長度，km；A為村域面積，km2。

6）聚耕比（Ag）

本文采用聚耕比指數(shù)來衡量研究區(qū)人地關系的協(xié)調程度。聚耕比指數(shù)越大，意味著相應區(qū)域內的農(nóng)村聚落的耕地資源越不足，人地矛盾越突出；聚耕比指數(shù)越小，意味著相應區(qū)域內鄉(xiāng)村聚落的耕地資源越豐富，人地矛盾越調和（劉香云，2019）。聚耕比計算方式：

式中，Ag 為聚耕比；La 為鄉(xiāng)鎮(zhèn)居民占地面積，m2；Fa 為鄉(xiāng)鎮(zhèn)耕地面積，m2。

7）到行政中心的距離（D）

行政中心附近居民生活配套設施較為完善，政策執(zhí)行效率較高，對于周邊的農(nóng)村居民具有較高吸引力。距離縣（區(qū)）級行政中心較近的鄉(xiāng)鎮(zhèn)由于交通等因素較為方便，農(nóng)民的離心力小于遠離縣（區(qū)）行政中心的鄉(xiāng)鎮(zhèn)（原野等，2015）。因此，本文選取電白區(qū)村級行政中心到電白區(qū)鎮(zhèn)級和區(qū)行政中心的距離，作為分析電白區(qū)農(nóng)村宅基地演化的驅動力因子之一。

8）人口密度（PopD）

人口密度越大的村鎮(zhèn)，村莊的內生的自我改造能力越大（原野等，2015）。合適的人口密度可以提高宅基地節(jié)約集約利用水平，從而有利于宅基地的擴張減少空心化的現(xiàn)象。本文所用的人口數(shù)據(jù)來自農(nóng)村地籍調查中的人口統(tǒng)計數(shù)據(jù)。

3 實驗及結果分析

3.1 時間變化特征

由圖1 看出，電白區(qū)從“一調”至2017 年，農(nóng)村宅基地占地面積增幅達79.86%。由時間段分析，從“一調”至“二調”，面積增幅為43.95%；“二調”至2017 年增幅為24.95%。由此可以得出，從“一調”至2017 年，電白區(qū)的農(nóng)村宅基地占地擴張迅速。其中，“一調”至“二調”的宅基地占地面積增長增幅，大于“二調”至2017 年的宅基地占地面積增長的增幅。這可能是因為“一調”至“二調”，電白區(qū)的農(nóng)村受到政策影響較小，且有較多的土地可用于宅基地建設，于是這一時期農(nóng)村宅基地呈現(xiàn)快速增長；“二調”至2017 年，電白區(qū)農(nóng)村逐漸受到相關農(nóng)村建房政策影響，同時，也因“一調”至“二調”占用了大量土地，可用于建設的土地在先前時期開發(fā)后逐漸減少，因此，這段時期的增幅比“一調”至“二調”的增幅要小。

圖1 不同時期宅基地面積增幅情況Fig.1 The increase of homestead area in different periods

由圖2、圖3 看出，“一調”至“二調”，宅基地面積增幅前五名的鎮(zhèn)（街道）分別是博賀鎮(zhèn)、坡心鎮(zhèn)、羊角鎮(zhèn)、南海街道、高地街道?！岸{”至2017 年宅基地面積增幅前五名的鎮(zhèn)（街道）是沙瑯鎮(zhèn)、林頭鎮(zhèn)、電城鎮(zhèn)、水東鎮(zhèn)、馬踏鎮(zhèn)，兩個時期增幅前五名中都沒有出現(xiàn)相同的鎮(zhèn)（街道），且“一調”至“二調”宅基地面積增幅前五名的增長幅度，遠大于“二調”至2017 年這一時間段。

圖2 “一調”至“二調”宅基地面積增幅前五名的鎮(zhèn)（街道）Fig.2 Top five towns (streets) with increased homestead area from the first national land survey to the second national land survey

圖3 “二調”至2017 年宅基地面積增幅前五名的鎮(zhèn)（街道）Fig.3 Top five towns (streets) with increased homestead area from the second national land survey to 2017

3.2 空間變化特征

3.2.1 農(nóng)村宅基地占地面積增長量

“一調”至2017 年，茂名市電白區(qū)農(nóng)村宅基地占地面積增長量的空間變化情況，如圖4 所示?！耙徽{”至2017 年，電白區(qū)農(nóng)村宅基地占地面積增長較多的村落主要集中在羊角鎮(zhèn)、坡心鎮(zhèn)、水東鎮(zhèn)、博賀鎮(zhèn)、高地街道及南海街道。分析增長的分布位置可以發(fā)現(xiàn)，這些村落主要集中在茂南區(qū)的毗鄰鎮(zhèn)區(qū)及沿海地區(qū)。說明茂南區(qū)作為茂名市中心城區(qū)，對于電白區(qū)鄉(xiāng)村規(guī)劃和發(fā)展具有積極影響；沿海區(qū)域是電白區(qū)水產(chǎn)品產(chǎn)業(yè)集群地區(qū)，農(nóng)村宅基地也出現(xiàn)大面積增長?？梢灶A見，隨著電白區(qū)水產(chǎn)品產(chǎn)業(yè)鏈日趨成熟，對于農(nóng)村規(guī)模的需求會不斷增加，沿海地區(qū)的農(nóng)村仍然還有巨大發(fā)展?jié)摿Αｋ姲讌^(qū)農(nóng)村宅基地增量較少的村落則集中在電白區(qū)的東北部，這部分地區(qū)遠離茂名市區(qū)，也缺乏沿海農(nóng)村發(fā)展的天然地域優(yōu)勢，因此，這些區(qū)位的農(nóng)村發(fā)展演變相對較為緩慢。

圖4 電白區(qū)“一調”至2017 年各階段村級宅基地占地增加面積分布Fig.4 Increased area distribution of village-level homestead in Dianbai District from the first national land survey to 2017

分時期分析農(nóng)村宅基地的變化情況，由圖5 可以看出，“一調”至“二調”，農(nóng)村宅基地面積增加較多的區(qū)域主要集中分布在羊角鎮(zhèn)、坡心鎮(zhèn)、高地街道和南海街道等電白區(qū)的西部區(qū)域，這段時期電白區(qū)農(nóng)村發(fā)展集中在靠近茂名市和電白區(qū)中心的地方；“二調”至2017 年，電白區(qū)的宅基地增長分布較為均勻，農(nóng)村宅基地增加較多的區(qū)域主要集中分布在羊角鎮(zhèn)、坡心鎮(zhèn)、林頭鎮(zhèn)、馬踏鎮(zhèn)、觀珠鎮(zhèn)等。

圖5 電白區(qū)不用時期農(nóng)村宅基地面積增量分布Fig.5 Incremental distribution of Dianbai District area rural homestead in unused period

3.2.2 農(nóng)村宅基地占地面積增幅

本文選用“一調”至2017 年村級農(nóng)村宅基地的增幅，來反映這一時期電白區(qū)各村莊宅基地面積變化水平，分析村級宅基地增幅結果，如圖6 所示。電白區(qū)從“一調”至2017 年農(nóng)村宅基地占地面積增長量較多，絕大多數(shù)農(nóng)村宅基地增幅大于0.5%。從空間上看，“一調”至2017 年農(nóng)村宅基地面積增幅大的區(qū)域分散分布，主要分布在博賀鎮(zhèn)、電城鎮(zhèn)、羊角鎮(zhèn)、坡心鎮(zhèn)、羅坑鎮(zhèn)與黃嶺鎮(zhèn)等電白區(qū)外圍的鎮(zhèn)區(qū)；在研究時間段內，電白區(qū)中部鎮(zhèn)區(qū)如關珠鎮(zhèn)、林頭鎮(zhèn)、霞洞鎮(zhèn)與馬踏鎮(zhèn)的村級農(nóng)村宅基地面積增長幅度比較小。結合圖4、圖5，可以發(fā)現(xiàn)，羅坑鎮(zhèn)在“一調”至“二調”和“二調”至2017 年農(nóng)村宅基地增長量都比較小，但在增幅上表現(xiàn)比較高，這是與該鎮(zhèn)區(qū)的“一調”時期的農(nóng)村宅基地占地面積比較小有關。

圖6 電白區(qū)“一調”至2017 年村級宅基地面積增幅分布Fig.6 Increased area distribution of village-level homestead in Dianbai District from the first national land survey to 2017

3.3 驅動力分析

將電白區(qū)的村級農(nóng)村宅基地面積增幅和親水性、到行政中心距離、路網(wǎng)密度、聚耕比、地形復雜度、人口密度進行相關性分析，分析結果如表3所示。電白區(qū)“一調”至2017 年，農(nóng)村宅基地占地面積增幅與親水性、到行政中心距離、路網(wǎng)密度、地形復雜度、人口密度呈顯著相關。其中，增幅與地形復雜度、到行政區(qū)中心距離在0.01 級別上呈顯著負相關，與聚耕比在0.05 級別上呈顯著負相關；增幅與親水性、路網(wǎng)密度、人口密度在0.01 級別上呈顯著正相關。表明電白區(qū)“一調”至2017 年，農(nóng)村宅基地的增長主要集中在親近水源、交通基礎條件好及人口多的村莊。在近水源、交通條件好、人口多的村莊，農(nóng)村宅基地占地面積增長較多；在地形復雜的地區(qū)及距離鎮(zhèn)、區(qū)行政中心較遠的村莊，農(nóng)村宅基地增長較慢；在聚耕比較大，即耕地面積比例較少的村莊，人地矛盾相對緊張，不利于村莊宅基地的擴張。

4 結 論

本文以電白區(qū)農(nóng)村地籍調查成果等數(shù)據(jù)分析了農(nóng)村宅基地的時間和空間變化特征，并進一步分析了影響電白區(qū)農(nóng)村宅基地演變的驅動因素，得出以下結論。

（1）從“一調”至2017 年，電白區(qū)農(nóng)村宅基地占地擴張迅速，農(nóng)村宅基地占地面積增幅達79.86%。分時間段看，從“一調”至“二調”，宅基地占地面積增幅為43.95%；“二調”至2017 年，占地面積增幅為24.95%。從空間上看，從“一調”至2017 年，電白區(qū)宅基地占地面積增加量較大的區(qū)域主要分布在西南部，南部的電城鎮(zhèn)與博賀鎮(zhèn)的面積增幅相較于其他鎮(zhèn)區(qū)更大。宅基地增加量較大的村落主要集中分布在茂名市區(qū)的毗鄰鎮(zhèn)區(qū)及沿海地區(qū)，說明市中心城區(qū)對于電白區(qū)鄉(xiāng)村規(guī)劃具有積極影響；沿海區(qū)域是電白區(qū)水產(chǎn)品產(chǎn)業(yè)集群地區(qū)，也對電白區(qū)農(nóng)村發(fā)展有積極影響。

（2）分析電白區(qū)農(nóng)村宅基地驅動力因素，電白區(qū)的“一調”至2017 年的農(nóng)村宅基地面積增幅與親水性、到行政區(qū)中心距離、路網(wǎng)密度、地形復雜度和人口密度等因子的相關性在0.01 級別上顯著相關。其中，電白區(qū)的面積增幅與親水性、路網(wǎng)密度和人口密度因子呈正向相關，與到行政區(qū)中心距離呈負相關，與聚耕比因子的相關性在0.05 級別上顯著負相關。因此，電白區(qū)“一調”至2017 年，農(nóng)村宅基地面積的增長主要發(fā)生在距離河流近、到行政中心近、路網(wǎng)密度高及地形平坦的區(qū)域。