王慧，王宏衛(wèi)，楊勝天,2,呂夢(mèng)婷，王媛媛，劉香云

(1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院，綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室，新疆 烏魯木齊 830046;2.北京師范大學(xué)地理科學(xué)學(xué)部，北京 100875)

鄉(xiāng)村聚落又稱農(nóng)村居民點(diǎn)，是農(nóng)村人口生活的場(chǎng)所，其形成和發(fā)展受自然、社會(huì)、經(jīng)濟(jì)因素的影響[1].黨的十九大報(bào)告提出高度重視“三農(nóng)”工作，堅(jiān)持農(nóng)業(yè)農(nóng)村優(yōu)先發(fā)展，實(shí)施鄉(xiāng)村振興戰(zhàn)略.揭示多種因素驅(qū)動(dòng)下的鄉(xiāng)村聚落空間格局特征，對(duì)當(dāng)前鄉(xiāng)村聚落布局優(yōu)化和美麗鄉(xiāng)村建設(shè)具有重要意義.

早在19世紀(jì)末，國(guó)外學(xué)者就已經(jīng)對(duì)鄉(xiāng)村聚落的形成、現(xiàn)狀和發(fā)展進(jìn)行了描述[2]。20世紀(jì)后趨向定性與定量相結(jié)合來研究聚落空間結(jié)構(gòu)、土地利用模式和聚落未來發(fā)展[3-4].國(guó)內(nèi)關(guān)于鄉(xiāng)村聚落的研究起步較晚，早期內(nèi)容多以農(nóng)村聚落地理、區(qū)劃、聚落形態(tài)分異和地域類型的研究為主[1,5-6].近年來，內(nèi)容上更加關(guān)注鄉(xiāng)村聚落景觀格局與驅(qū)動(dòng)因素[7-10]，聚落空間演變及其LUCC響應(yīng)分析[1,11-13]，鄉(xiāng)村聚落空間重構(gòu)與結(jié)構(gòu)優(yōu)化[14-15]，鄉(xiāng)村聚落空心化測(cè)度[16-17].研究區(qū)域和方法上，朱彬[6,18]，鄂施璇[10]，李紅波等[19]用景觀分析方法、GIS空間分析方法對(duì)平原地區(qū)的鄉(xiāng)村聚落格局及影響因素進(jìn)行了研究；郭曉東[20-21]，譚雪蘭等[22]更多采用GIS分析及定性分析法分析山地、丘陵地區(qū)鄉(xiāng)村聚落格局的人文社會(huì)影響因素；岳邦瑞等[23]定性地探討了基于水資源約束下的綠洲鄉(xiāng)土聚落形態(tài)特征.總體來說，有關(guān)鄉(xiāng)村聚落格局及影響因素的研究多在內(nèi)地平原、山地、丘陵區(qū)域，對(duì)于西北干旱區(qū)尤其是極端干旱區(qū)新疆綠洲鄉(xiāng)村聚落的研究較少，且定性的研究居多，定量化的分析仍顯不足.運(yùn)用多種方法定量分析西北干旱區(qū)鄉(xiāng)村聚落分布格局特征及影響因素至關(guān)重要.

烏蘇市是奎屯河的主要流經(jīng)區(qū)域，地貌類型復(fù)雜多樣，鄉(xiāng)村聚落分布格局深受內(nèi)部海拔、坡度、河流水系等因素的影響，鑒于此，以烏蘇市為研究對(duì)象，運(yùn)用最近鄰點(diǎn)統(tǒng)計(jì)量、核密度估算、空間關(guān)聯(lián)測(cè)度模型、景觀形狀指數(shù)等方法研究烏蘇市鄉(xiāng)村聚落的空間格局特征并把握其發(fā)展規(guī)律，在此基礎(chǔ)上利用多元Logistic回歸模型定量地分析自然和人文社會(huì)因素對(duì)鄉(xiāng)村聚落空間分布格局的影響程度，為促進(jìn)西北干旱區(qū)城鄉(xiāng)一體化發(fā)展、鄉(xiāng)村聚落的布局優(yōu)化和美麗鄉(xiāng)村建設(shè)提供理論參考.

1 材料與方法

1.1 研究區(qū)概況

本文選取地處奎屯河流域、準(zhǔn)噶爾盆地西南部的烏蘇市為研究區(qū)域(圖1)，所屬新疆塔城地區(qū)，總面積為2.07萬km2，下轄10鎮(zhèn)、7鄉(xiāng)、3個(gè)牧場(chǎng)、1個(gè)林場(chǎng)和5個(gè)街道，此外還包括新疆生產(chǎn)建設(shè)兵團(tuán)所屬農(nóng)七師的123～128團(tuán)6個(gè)團(tuán)場(chǎng)，地勢(shì)南高北低，北部是準(zhǔn)噶爾盆地，有很大部分的半沙漠，中部是盆地—山地過渡帶，南部是北天山山地，境內(nèi)海拔在208～4 898 m之間，為典型的溫帶大陸性氣候，年內(nèi)溫差變幅為42.8 ℃，光照充足，實(shí)際日照時(shí)長(zhǎng)數(shù)2 600～2 800 h，主要靠降水和冰雪融水補(bǔ)給水源，以種植小麥、玉米、棉花為主，2016年底全市國(guó)民生產(chǎn)總值達(dá)160.1億元.

1.2 數(shù)據(jù)來源及處理

1.2.1 數(shù)據(jù)來源 基于Google earth遙感影像，利用Arcgis10.2軟件，經(jīng)過幾何校正、地理坐標(biāo)配準(zhǔn)、人工目視解譯等獲得烏蘇市鄉(xiāng)村聚落圖斑和鄉(xiāng)鎮(zhèn)行政區(qū)劃等矢量數(shù)據(jù)，其中鄉(xiāng)村聚落圖斑為建制鎮(zhèn)和村莊；從烏蘇市2016年Landsat8遙感影像(分辨率30 m)中提取道路、河流等矢量數(shù)據(jù)；烏蘇市DEM高程數(shù)據(jù)(分辨率30 m)和坡度數(shù)據(jù)(分辨率90 m)，來源于地理空間數(shù)據(jù)云.

圖1 研究區(qū)示意圖Figure 1 Location of the study area

基于數(shù)據(jù)的可獲取性和科學(xué)性，選取了表征空間分異的斑塊面積、斑塊密度、形狀指數(shù)3個(gè)因變量和坡度、高程、距國(guó)道距離、距省道距離、距縣道距離、距主要河流距離、距鄉(xiāng)鎮(zhèn)中心距離7個(gè)自變量.

1.2.2 數(shù)據(jù)處理 計(jì)算3個(gè)因變量均值的差值，對(duì)高程和坡度進(jìn)行重分類，對(duì)坡度0～2°、>2～4°和>4°分別賦予1、2、3，其余5個(gè)因變量利用Arcgis10.2中的Feature To Point模塊提取斑塊中心點(diǎn)，再用空間分析模塊計(jì)算斑塊到各因子的空間距離；將所有矢量數(shù)據(jù)轉(zhuǎn)成分辨率50 m×50 m的柵格數(shù)據(jù).

1.3 研究方法

1.3.1 最近鄰點(diǎn)統(tǒng)計(jì)量 判斷點(diǎn)要素的空間分布屬于隨機(jī)、集聚還是分散型.公式如下：

(1)

(2)

Z=(r0-re)/SEr

(3)

式中，r0表示平均距離觀測(cè)值；re表示平均距離期望值；di為聚落i的最近鄰距離；n為聚落總數(shù)；A為總面積；SEr為標(biāo)準(zhǔn)誤差[1].R<1表示聚落分布呈集聚型，R>1為分散型.然后進(jìn)行Z檢驗(yàn)，最終判斷點(diǎn)要素的分布模式是否顯著.置信度在95%時(shí)，如果Z>1.96或Z<-1.96且R<1，表明在α=0.05顯著水平下聚落集聚分布態(tài)勢(shì)顯著，如果Z>1.96或Z<-1.96且R>1，表明在α=0.05顯著水平下聚落分散分布態(tài)勢(shì)顯著.

1.3.2 核密度估算 (KDE) 核密度用來表示點(diǎn)要素的空間分布特征，其模型如下：

(4)

式中，f(x,y)表示位置 (x,y)的密度估計(jì)；n表示觀測(cè)數(shù)；h表示參數(shù)；k表示核函數(shù)，di表示位置(x,y)與第i個(gè)位置的距離[1].

1.3.3 空間關(guān)聯(lián)測(cè)度模型 全局空間聚類檢驗(yàn)(Getis-Ord General G)用來檢驗(yàn)鄉(xiāng)村聚落全局空間規(guī)模呈現(xiàn)高值集聚還是低值集聚的特征[24].計(jì)算公式如下：

(5)

空間“熱點(diǎn)”探測(cè)(Getis-OrdGi*)用來檢驗(yàn)鄉(xiāng)村聚落局部地區(qū)存在統(tǒng)計(jì)顯著的高值還是低值特征，公式如下：

(6)

式中，參數(shù)表示與前面相同.對(duì)Gi*(d)進(jìn)行標(biāo)準(zhǔn)化處理得Z(Gi*) =(Gi*-E(Gi*))/var(Gi*)，如果Z(Gi*)>0，且統(tǒng)計(jì)顯著，表明位置i周圍值較高，屬高值集聚的熱點(diǎn)區(qū)，如果Z(Gi*)<0，且統(tǒng)計(jì)顯著，表明位置i周圍的值較低，屬低值集聚的冷點(diǎn)區(qū).

1.3.4 多元Logistic回歸模型 在研究區(qū)分層隨機(jī)抽樣4 000個(gè)樣本點(diǎn)，再從總樣本中抽取3 000個(gè)樣本點(diǎn)參與計(jì)算，剩余樣本點(diǎn)進(jìn)行對(duì)比分析.因變量在區(qū)間[0,1]內(nèi)大體相等，對(duì)多分類型自變量坡度進(jìn)行啞變量設(shè)置，坡度Ⅰ(0～2°)和Ⅱ(>2～4°)為啞變量，坡度Ⅲ(>4°)為參照對(duì)象，利用SPSS 17.0軟件計(jì)算回歸模型的相關(guān)系數(shù)，并通過回歸系數(shù)來反映自變量與因變量之間的相關(guān)關(guān)系.

2 結(jié)果與分析

2.1 鄉(xiāng)村聚落空間分布格局特征

2.1.1 鄉(xiāng)村聚落空間分布特征 將聚落斑塊面要素轉(zhuǎn)為點(diǎn)要素，計(jì)算最鄰近點(diǎn)統(tǒng)計(jì)量(R)，結(jié)果如表1所示：

由表1可知，R統(tǒng)計(jì)量為0.505<1，表明烏蘇市鄉(xiāng)村聚落分布比隨機(jī)模式集聚，標(biāo)準(zhǔn)化Z值為-38.875，小于-1.96，烏蘇市鄉(xiāng)村聚落整體集聚態(tài)勢(shì)較為顯著，主要集中在中部盆地與盆地-低山過渡地區(qū).

表1 最近鄰點(diǎn)統(tǒng)計(jì)量

由圖2可知，烏蘇市鄉(xiāng)村聚落密度最高區(qū)域?yàn)?.502～4.795個(gè)/km2，密度最低區(qū)域聚落個(gè)數(shù)不足1個(gè)甚至沒有.烏蘇市鄉(xiāng)村聚落呈現(xiàn)低密度-多核心的空間分異特征，密度核心區(qū)主要分布在中部頭臺(tái)鄉(xiāng)、九間樓鄉(xiāng)、八十四戶鄉(xiāng)、西湖鎮(zhèn)以及古爾圖鎮(zhèn).密度次核心區(qū)北部主要分布在126團(tuán)、123團(tuán)、125團(tuán)團(tuán)場(chǎng)各連隊(duì)、石橋鄉(xiāng)，南部密度次核心有四棵樹鎮(zhèn)、哈圖布呼鎮(zhèn)、百泉鎮(zhèn)這幾個(gè)鄉(xiāng)鎮(zhèn)連片分布區(qū)以及西大溝鎮(zhèn).密度核心和次核心區(qū)多分布在河流干渠，公路、省道、縣道附近.

圖2 鄉(xiāng)村聚落密度圖Figure 2 The spatial density of rural settlements

2.1.2 鄉(xiāng)村聚落規(guī)模特征 從表2結(jié)果分析得出，鄉(xiāng)村聚落觀測(cè)值G(d)(0.000 015)<期望值E(d)(0.000 019)，Z(d)為-2.813 363，統(tǒng)計(jì)顯著，表明烏蘇市鄉(xiāng)村聚落規(guī)模分布具有顯著的低值集聚特征，即村莊呈集群分布且規(guī)模普遍偏小.

表2 鄉(xiāng)村聚落規(guī)模General G估計(jì)值

運(yùn)用自然斷點(diǎn)法將Gi*分為7級(jí)，得到鄉(xiāng)村聚落規(guī)模熱點(diǎn)分布圖(圖3).紅色區(qū)域?yàn)闊狳c(diǎn)區(qū)即高值簇群，藍(lán)色區(qū)域?yàn)闊狳c(diǎn)區(qū)即低值簇群，由圖3可知，熱點(diǎn)區(qū)多為兵團(tuán)團(tuán)場(chǎng)連隊(duì)，即兵團(tuán)團(tuán)場(chǎng)連隊(duì)為較大規(guī)模的鄉(xiāng)村聚落集群，冷點(diǎn)區(qū)多為地方鄉(xiāng)鎮(zhèn)村莊，即地方鄉(xiāng)鎮(zhèn)村莊為較小規(guī)模的鄉(xiāng)村聚落集群.“熱點(diǎn)區(qū)”即大規(guī)模集聚分布的區(qū)域主要有北部奎屯河下游的123團(tuán)、125團(tuán)、126團(tuán)、127團(tuán)、128團(tuán)團(tuán)場(chǎng)各連隊(duì)以及石橋鄉(xiāng)車牌子鎮(zhèn)一小部分；“冷點(diǎn)區(qū)”主要分布在古爾圖鎮(zhèn)、頭臺(tái)鄉(xiāng)、西湖鎮(zhèn)、甘河子鎮(zhèn)、八十四戶鄉(xiāng)、西大溝鎮(zhèn)以及甘家湖林場(chǎng)的西北小部分區(qū)域.結(jié)合烏蘇市鄉(xiāng)村聚落密度分布圖，發(fā)現(xiàn)它們?cè)诳臻g上呈現(xiàn)低密度-大規(guī)模和高密度-小規(guī)模的分布特征.

2.1.3 鄉(xiāng)村聚落形態(tài)特征 運(yùn)用景觀形狀指數(shù)來分析鄉(xiāng)村聚落形狀的復(fù)雜程度，將鄉(xiāng)村聚落斑塊要素與烏蘇市格網(wǎng)圖層(2 km×2 km)進(jìn)行空間疊置，以格網(wǎng)化后的鄉(xiāng)村聚落斑塊為變量，得到聚落形狀指數(shù)空間分布圖(圖4).由圖4可知，烏蘇市鄉(xiāng)村聚落大部分分布在綠洲內(nèi)部，形狀指數(shù)值總體不高，波動(dòng)較小，烏蘇市鄉(xiāng)村聚落形狀較為規(guī)整，穩(wěn)定性較好，形狀指數(shù)高值區(qū)多靠近國(guó)道、省道、縣道，形狀指數(shù)低值區(qū)域多為距離道路較遠(yuǎn)，鄉(xiāng)村聚落以單戶或規(guī)模較小的戶數(shù)分布為主，聚落結(jié)構(gòu)較為簡(jiǎn)單，因此形狀指數(shù)較低.

圖3 烏蘇市鄉(xiāng)村聚落規(guī)模分布圖Figure 3 Distribution map of rural settlements in Wusu City

2.2 鄉(xiāng)村聚落空間格局影響因素

2.2.1 鄉(xiāng)村聚落分布差異的Logistic回歸模型分析 自變量容忍度在0.714～0.818之間，可以用來做回歸模型分析，根據(jù)表3模型分析的結(jié)果得知，坡度對(duì)聚落分布具有一定影響，高程沒有進(jìn)入模型分析結(jié)果，說明烏蘇市聚落分布受高程影響不大.對(duì)聚落斑塊密度重要的影響變量是到最近縣道距離、到最近主要河流距離和到最近鄉(xiāng)鎮(zhèn)中心距離.距縣道距離的回歸系數(shù)為-0.465，主要原因是縣道作為鄉(xiāng)村人流、物流交往的媒介，對(duì)鄉(xiāng)村聚落靠近具有很強(qiáng)的吸引力.所以距縣道的距離越近聚落分布越密集；距主要河流距離的回歸系數(shù)為-0.346，主要原因是烏蘇市境內(nèi)雖然有奎屯河、四棵樹河以及古爾圖河以及一系列支流，但地處于西北干旱區(qū)，鄉(xiāng)村農(nóng)業(yè)生產(chǎn)灌溉和生活具有很強(qiáng)的“親水性”，所以隨著河流距離的增加，聚落密度減少；距鄉(xiāng)鎮(zhèn)中心距離的回歸系數(shù)為-0.012，主要是由于工業(yè)化和交通運(yùn)輸業(yè)的發(fā)展吸引著鄉(xiāng)村各要素向鄉(xiāng)鎮(zhèn)靠近，鄉(xiāng)鎮(zhèn)周邊村莊成為人口的聚集地，故距離鄉(xiāng)鎮(zhèn)中心越近聚落越密集.

圖4 烏蘇市鄉(xiāng)村聚落形狀指數(shù)空間分布圖Figure 4 The spatial shape index of rural settlements in Wusu City

2.2.2 鄉(xiāng)村聚落規(guī)模分異的Logistic回歸模型分析 為了變量之間出現(xiàn)多重共線性問題，對(duì)7個(gè)自變量進(jìn)行分析，容忍度顯示在0.534～0.657之間，可以用來做回歸模型分析，根據(jù)顯著性水平(P<0.05)和Waldx2卡方統(tǒng)計(jì)量得知(表3)，對(duì)規(guī)模差異具有重要的影響變量有高程、坡度.高程(Waldx2為99.245)的回歸系數(shù)為-1.032，說明高程每增加100 m聚落斑塊面積減少的概率為1.587倍，主要原因是研究區(qū)內(nèi)海拔在208～4 898 m之間，地貌類型復(fù)雜多樣，北部主要是沖積平原和半沙漠構(gòu)成的盆地，中部是沖積洪積平原和低山構(gòu)成的盆地—山地過渡帶，海拔較低，南部是侵蝕中山、冰緣、冰川作用高山構(gòu)成的天山山地，海拔較高，由北向南隨著海拔高度的不斷增加，聚落個(gè)數(shù)和聚落面積呈減少的趨勢(shì)，從而影響規(guī)模分異；坡度(Waldx2為64.134)的回歸系數(shù)為-0.359，說明坡度每增加2°聚落斑塊面積減少的概率為1.755倍，主要是因?yàn)闉跆K市自北向南由盆地(相當(dāng)部分沙漠)、盆地—低山過渡帶、山地組成，氣候地貌的垂直地帶性分異明顯，由北向南依次為荒漠氣候、平原氣候、山麓氣候和山地氣候，氣候地貌使得烏蘇市鄉(xiāng)村聚落多分布在北部和中部的盆地、盆地—低山過渡帶這些坡度較小的區(qū)域，而南部北天山山地坡度大的地方以山麓、山地氣候?yàn)橹鳎瑲夂蚝溆辛阈欠植嫉哪翗I(yè)隊(duì)和林場(chǎng).

2.2.3 鄉(xiāng)村聚落形狀分異的Logistic回歸模型分析 自變量容忍度在0.964～8.564之間，可以用來做回歸模型分析，由表3得知，高程、坡度、距最近鄉(xiāng)鎮(zhèn)中心沒有進(jìn)入模型分析結(jié)果，說明對(duì)聚落形態(tài)分異的影響不大，對(duì)聚落形態(tài)分異主要的影響變量是距縣道距離、距省道距離.距縣道距離的回歸系數(shù)為-0.267，即距離每增加1 km形狀指數(shù)減少1.651倍，形狀越規(guī)整，破碎化程度越低，主要是因?yàn)榭h道作為鄉(xiāng)村各要素交流的媒介，鄉(xiāng)村聚落為了降低交通和運(yùn)輸成本而靠近縣道，進(jìn)而影響鄉(xiāng)村聚落空間的規(guī)整化.距省道距離的回歸系數(shù)為-0.148，主要是因?yàn)槭〉缹?duì)于鄉(xiāng)村聚落面狀要素的切割作用大于聯(lián)通，進(jìn)而形成鄉(xiāng)村聚落空間的破碎化.

表3 鄉(xiāng)村聚落空間分布格局的Logistic回歸模型相關(guān)系數(shù)

3 結(jié)論

本研究以烏蘇市為例，運(yùn)用最近鄰點(diǎn)統(tǒng)計(jì)量、核密度估算、空間關(guān)聯(lián)測(cè)度模型、景觀形狀指數(shù)等方法研究西北干旱區(qū)鄉(xiāng)村聚落景觀的空間分布特征并把握其發(fā)展規(guī)律，并用多元Logistic回歸模型分析鄉(xiāng)村聚落空間分布格局的影響因素，為當(dāng)前鄉(xiāng)村聚落景觀的布局優(yōu)化和美麗鄉(xiāng)村建設(shè)提供理論參考.

鄉(xiāng)村聚落分布整體呈集聚分布模式，密度分布呈現(xiàn)低密度-多核心的空間分異特征，密度核心和次核心區(qū)多分布在河流干渠，公路、省道、縣道附近.

烏蘇市鄉(xiāng)村聚落規(guī)模分布具有顯著低值集聚特征，即村莊呈集群分布且規(guī)模普遍偏?。粺狳c(diǎn)區(qū)多為兵團(tuán)團(tuán)場(chǎng)連隊(duì)，即兵團(tuán)團(tuán)場(chǎng)連隊(duì)為較大規(guī)模的鄉(xiāng)村聚落集群，冷點(diǎn)區(qū)多為地方鄉(xiāng)鎮(zhèn)村莊，即地方鄉(xiāng)鎮(zhèn)村莊為較小規(guī)模的鄉(xiāng)村聚落集群；通過將烏蘇市鄉(xiāng)村聚落規(guī)模與密度分布圖相比，發(fā)現(xiàn)它們?cè)诳臻g上有較明顯的局部負(fù)相關(guān)關(guān)系，即呈現(xiàn)低密度-大規(guī)模和高密度-小規(guī)模的分布特征.

烏蘇市鄉(xiāng)村聚落大部分分布在綠洲內(nèi)部，形狀指數(shù)值總體不高，波動(dòng)較小，表明烏蘇市鄉(xiāng)村聚落形狀較為規(guī)整，穩(wěn)定性較好，形狀指數(shù)高值區(qū)是公路、國(guó)道、省道、縣道以及由河流干渠經(jīng)過或附近區(qū)域.

聚落分布差異的主要影響因素是距縣道距離、距主要河流距離和距鄉(xiāng)鎮(zhèn)中心距離；聚落規(guī)模分異的主要影響因素是坡度、高程；聚落形態(tài)分異的主要影響因素是距縣道距離、距省道距離.

4 討論

明確不同影響因素對(duì)鄉(xiāng)村聚落格局的影響程度至關(guān)重要，在研究西北干旱區(qū)烏蘇市鄉(xiāng)村聚落分布格局的影響因素過程中，地形地勢(shì)對(duì)聚落規(guī)模具有較大影響，此外，受縣道、鄉(xiāng)鎮(zhèn)中心和河流的引力以及省道的切割作用，使得聚落分布和形態(tài)發(fā)生改變，最終導(dǎo)致烏蘇市鄉(xiāng)村聚落格局呈現(xiàn)分布稀疏、規(guī)模較小、形態(tài)較為規(guī)整的空間特征.鄉(xiāng)村聚落作為人口分布的載體，為更好地響應(yīng)中共十九大提出的鄉(xiāng)村全面振興戰(zhàn)略，各級(jí)政府加快發(fā)展城鄉(xiāng)之間互動(dòng)關(guān)系，推動(dòng)鄉(xiāng)村聚落的可持續(xù)發(fā)展.

此外，文中運(yùn)用多元Logistic回歸模型對(duì)影響鄉(xiāng)村聚落格局分異的自然和人文社會(huì)因素進(jìn)行定量的分析，但研究程度有待深入；囿于數(shù)據(jù)的可獲取性，研究局限于單一時(shí)間點(diǎn)對(duì)鄉(xiāng)村聚落格局和影響因素進(jìn)行分析，沒有時(shí)空序列組合來研究鄉(xiāng)村聚落空間分異的復(fù)雜過程.