尤 笛， 蔣金豹， 郭海強， 尤 倩

[1.中國礦業(yè)大學(xué)(北京)地球科學(xué)與測繪工程學(xué)院，北京 100083； 2.沈陽工業(yè)大學(xué)理學(xué)院，遼寧沈陽 110870]

在國家“一帶一路”戰(zhàn)略構(gòu)想下，國家發(fā)展和改革委員會等聯(lián)合發(fā)布《推動共建絲綢之路經(jīng)濟帶和21世紀(jì)海上絲綢之路的愿景與行動》，打造以福建省為“海上絲路”核心的經(jīng)濟樞紐。福建省有著“八山一水一分田”的地理格局，境內(nèi)山地丘陵約占90%，建設(shè)用地嚴(yán)重不足，無法承載21世紀(jì)海上絲綢之路核心區(qū)的建設(shè)。解決這一問題的關(guān)鍵就是加快土地后備資源的開發(fā)，即低丘緩坡的建設(shè)用地空間拓展。為此，2014年4月在福建省國土資源廳下發(fā)《關(guān)于創(chuàng)建節(jié)地型、生態(tài)型低丘緩坡土地開發(fā)利用示范(片)區(qū)的通知》后，將邵武市確定為低丘緩坡地示范區(qū)。而在開發(fā)低丘緩坡地之前，須要對其建設(shè)用地的空間拓展進行適宜性評價，以期提供科學(xué)合理的參考。低丘緩坡指地形起伏適中，坡度小于25°且具備一定成片開發(fā)利用規(guī)模的山坡地和荒灘區(qū)域[1]，廣泛分布在我國東南沿海，浙閩交界的武夷山脈區(qū)域。近年來，低丘緩坡土地資源的適宜性評價越來越受到學(xué)者的關(guān)注。羅鵬等基于“3S”技術(shù)，結(jié)合特爾斐法、層次分析法對浙江金衢地區(qū)低丘紅壤土地的開發(fā)利用潛力進行綜合評價[2]；朱曉蕓以浙江省余姚市為例，在GIS和RS支持下，充分考慮地形、植被覆蓋、土壤及區(qū)位等因素，對研究區(qū)低丘緩坡耕地開發(fā)利用潛力進行評價[3]；徐萍等通過對浙江省低丘緩坡開發(fā)利用現(xiàn)狀的評價，提出低丘緩坡農(nóng)業(yè)資源高效集約利用的對策措施[4]。上述研究主要針對農(nóng)用地適宜性進行評價，但對低丘緩坡建設(shè)用地適宜性評價的研究較少，且已有研究常用特爾斐法(delphi method)[5]、層次分析法(analytic hierarchy process，AHP)確定權(quán)重[6]，然后將多個指標(biāo)疊加得出評價結(jié)果分值，有一定的主觀片面性。此外，在評價指標(biāo)體系的構(gòu)建上也主要從地形地貌、地質(zhì)災(zāi)害角度來考慮[7]，指標(biāo)體系不夠全面。本研究在傳統(tǒng)土地適宜性評價的基礎(chǔ)上，多角度考慮空間區(qū)位、自然地理、社會經(jīng)濟、生態(tài)環(huán)境等因素，采用主成分分析法(principal component analysis，PCA)對福建省邵武市低丘緩坡建設(shè)用地的適宜性進行綜合評價[8]，并通過K-均值(K-means)聚類分析將評價結(jié)果分為4個適宜性等級[9]。該方法不僅能消除指標(biāo)之間的相關(guān)性，還能避免確定權(quán)重的主觀性。

1 研究區(qū)及數(shù)據(jù)處理

1.1 研究區(qū)概況

邵武素有“鐵城”之稱，史稱南武夷，隸屬于福建省南平市，位于福建省西北部，東北鄰南平市的建陽區(qū)，東南連順昌縣，南接三明市的將樂縣、泰寧縣，西與江西省黎川縣毗鄰，西北與光澤縣交界。境內(nèi)東西寬115 km，南北長120 km，占地面積2 859.4 km2。全境大部分由山地、丘陵構(gòu)成，地形由西南和東北向中部河谷傾斜，四界接壤處地勢均較高，1 000 m以上的山峰有130座，森林覆蓋率達75.7%。其中中山(海拔800 m以上)和低山(海拔500～800 m)約占全市總面積的39.7%，高丘陵(海拔250～500 m)和丘陵約占41.6%，低丘緩坡資源非常豐富(圖1)。

1.2 數(shù)據(jù)處理及評價流程

根據(jù)福建省城鄉(xiāng)規(guī)劃設(shè)計研究院提供的邵武市低丘緩坡范圍剔除了禁止開發(fā)區(qū)(3個剛性因子區(qū)域)，以30 m×30 m柵格為評價單元，對所有彈性因子數(shù)據(jù)進行標(biāo)準(zhǔn)化處理，以消除量綱不同和量級差異的影響。為了研究方便，將適度指標(biāo)“坡向”轉(zhuǎn)化為正向指標(biāo)[10]，這樣所有指標(biāo)只有正向指標(biāo)和負(fù)向指標(biāo)。轉(zhuǎn)化方法如下：將坡向定性分為陰坡、半陰坡、半陽坡、陽坡4類(圖2)，并分別賦予分值1、3、5、7分，這樣坡向指標(biāo)的相關(guān)性就轉(zhuǎn)化為正相關(guān)，即正向指標(biāo)。

常用的標(biāo)準(zhǔn)化處理方法有均值化法、等測度法、極差正規(guī)化法等。本研究采用極差正規(guī)化法進行標(biāo)準(zhǔn)化[11]，公式如下：當(dāng)指標(biāo)是正向指標(biāo)(即正相關(guān))時，

(1)

當(dāng)指標(biāo)是負(fù)向指標(biāo)(即負(fù)相關(guān))時，

(2)

2 主成分聚類分析評價

2.1 選取評價指標(biāo)

選取評價指標(biāo)對評價結(jié)果的合理性影響很大。本著科學(xué)性、完整性、層次性、可比性等原則，結(jié)合邵武市獨特的地形地貌，從空間區(qū)位、自然地理、社會經(jīng)濟、生態(tài)環(huán)境4個方面選取12個彈性因子和3個剛性因子構(gòu)建了評價指標(biāo)體系[12](表1)。

2.2 主成分綜合評價

PCA是通過正交變換將多個指標(biāo)轉(zhuǎn)化為保留原始變量絕大部分信息的少數(shù)指標(biāo)， 以達到削弱噪聲和數(shù)據(jù)降維的一種統(tǒng)計方法[13]。由于評價指標(biāo)很多且彼此之間有一定的相關(guān)性，從而反映的信息有所重疊，且在確定權(quán)重時也存在主觀片面性，而主成分分析恰好能解決這些問題[14]。

表1 邵武市低丘緩坡地評價指標(biāo)體系

注：“+”表示正相關(guān)，“-”表示負(fù)相關(guān)。剛性因子，即風(fēng)景名勝及特殊用地、基本農(nóng)田保護區(qū)、森林禁止開發(fā)區(qū)，是明文規(guī)定的禁止開發(fā)區(qū)，所以將這3個因子直接劃為建設(shè)用地禁止開發(fā)區(qū)，即不適宜建設(shè)區(qū)。實際參與主成分評價的只有12個彈性因子。

評價原理如下[15]：(1)設(shè)評價單元數(shù)為n，每個評價單元p個指標(biāo)。當(dāng)前k(k

F1=a11x1+a21x1+…+ap1xp；

(3)

F2=a21x1+a22x2+…+ap2xp；

(4)

Fk=a1kx1+a2kx2+…+apkxp。

(5)

式中：aij表示第i個評價單元在第j個主因子上的權(quán)重系數(shù)。

(2)把不同主分量特征值的方差貢獻率βi(i=1，2，…，k)作為加權(quán)系數(shù)，利用綜合評價函數(shù)F=β1F1+β2F2+…+βkFk計算各評價單元的分值。

根據(jù)上述原理，將12個彈性指標(biāo)合成1個多波段數(shù)據(jù)后，再進行PCA處理，得到各主成分的特征值和貢獻率(表2)。

表2 各主成分的特征值及貢獻率

注：Z1～Z12分別表示前12個主成分。

由表2可知，前6個主成分方差累積貢獻率87.08%>85%，基本保留了原始變量的所有信息，所以選取前6個主成分代表所有變量信息。根據(jù)主成分綜合評價公式F=β1F1+β2F2+…+βkFk，建立綜合評價模型：

F=0.230 8F1+0.218 0F2+0.168 5F3+0.102 2F4+0.093 8F5+0.057 5F6。

得到的主成分綜合評價結(jié)果見圖4。

2.3 K-means聚類分析

聚類分析是一種探索性的分析，它將數(shù)據(jù)劃分到不同的類，類內(nèi)的對象有很大相似性，而類間的對象有很大差異性。其中K-means使用較廣泛，其運算簡潔快速，在分類的過程中它能從樣本數(shù)據(jù)出發(fā)自動進行分類，而不需要樣本先驗信息[16]。聚類步驟如下：(1)聚類中心初始化；(2)計算待分類對象與聚類中心的歐式距離，依據(jù)就近原則將待分類對象歸類；(3)計算各初始類別的均值作為新的聚類中心；(4)按照新聚類中心重新聚類，如此反復(fù)循環(huán)，直至聚類中心不再變化或小于給定閾值時，終止循環(huán)。

根據(jù)K-means聚類算法，將主成分綜合評價結(jié)果進行K-means 聚類，分為最適宜、較適宜、一般適宜、不適宜4類。再將主成分評價前剔除的風(fēng)景名勝及特殊用地、基本農(nóng)田保護區(qū)、森林禁止開發(fā)區(qū)這些劃分為不適宜區(qū)進行合并，最終形成分級評價結(jié)果(圖5)。

3 評價結(jié)果與分析

根據(jù)ArcGIS區(qū)域分析功能，得到邵武市各鄉(xiāng)(鎮(zhèn))適宜性建設(shè)評價結(jié)果(表3)。

3.1 定性分析

由圖5可知，適宜建設(shè)的低丘緩坡地主要集中在中心城區(qū)，此區(qū)域交通便利，是鐵路、公路交通樞紐所在地；具有人口密集、基礎(chǔ)設(shè)施齊全、城鎮(zhèn)化水平較高、財政收入較高、地勢平坦等優(yōu)點，不管是從社會經(jīng)濟角度考慮還是從施工難度來看，都較適宜開發(fā)建設(shè)用地。根據(jù)現(xiàn)有資料，此地帶地質(zhì)災(zāi)害易發(fā)性較小，生態(tài)敏感性較弱，開發(fā)過程中不易破壞生態(tài)系統(tǒng)，這樣既可以實現(xiàn)城鎮(zhèn)的可持續(xù)發(fā)展，又節(jié)約了為恢復(fù)生態(tài)建設(shè)所需的成本。不適宜區(qū)主要分布在邵武市的周邊，這些地區(qū)海拔較高、坡度較大、森林覆蓋度較大、社會經(jīng)濟水平較低，開發(fā)起來不僅對施工技術(shù)要求較高、成本較高，還會破壞生態(tài)環(huán)境，引發(fā)地質(zhì)災(zāi)害，導(dǎo)致環(huán)境的惡性循環(huán)。

3.2 定量分析

由表3可知，邵武低丘緩坡面積為119 512.71 hm2，其中最適宜建設(shè)面積為18 983.88 hm2，占低丘緩坡地的15.88%，占邵武市總面積6.64%；較適宜建設(shè)面積為25 545.69 hm2，占低丘緩坡地的21.37%，占邵武市總面積8.93%；一般適宜建設(shè)面積為27 415.53 hm2，占低丘緩坡地的22.94%，占邵武市總面積9.59%；不適宜建設(shè)面積為47 567.61 hm2，占低丘緩坡地的39.80%，占邵武市總面積16.64%。水北鎮(zhèn)、拿口鎮(zhèn)、沿山鄉(xiāng)低丘緩坡資源最豐富，分別占總低丘緩坡面積的15.79%、13.86%、9.23%，合計達38.88%。最適宜建設(shè)面積最大的是城郊鎮(zhèn)，占最適宜建設(shè)總面積的20.69%；較適宜建設(shè)面積最大的是水北鎮(zhèn)，占較適宜建設(shè)總面積的21.80%；一般適宜建設(shè)面積最大的也是水北鎮(zhèn)，占一般適宜建設(shè)總面積的18.21%；不適宜建設(shè)面積最大的是拿口鎮(zhèn)，占不適宜建設(shè)總面積的13.97%。

3.3 適宜性區(qū)域分析及建議

最適宜建設(shè)面積較大的鄉(xiāng)(鎮(zhèn))主要有城郊鎮(zhèn)、吳家塘鎮(zhèn)、水北鎮(zhèn)、曬口街道。其中水北鎮(zhèn)是邵武市占地面積最大的鄉(xiāng)鎮(zhèn)，位于邵武市北部，富屯溪、鷹廈鐵路、316國道、205省道貫穿境內(nèi)，交通便利；城郊鎮(zhèn)位于邵武市中部，距邵武市人民政府很近，平均海拔300 m，境內(nèi)溪流交錯，京福高速公路引線橫穿全境，是邵武市低丘緩坡開發(fā)項目片區(qū)；吳家塘鎮(zhèn)是綜合性的企業(yè)鄉(xiāng)鎮(zhèn)，已經(jīng)初步形成邵武鐵路線上一座新興的集鎮(zhèn)；曬口街道位于邵武市中心城區(qū)，地勢平坦，基礎(chǔ)設(shè)施完備，常住人口較多。這些區(qū)域最適宜開發(fā)建設(shè)用地，應(yīng)進一步加強土地節(jié)約集約利用，充分挖掘土地利用潛力，優(yōu)化產(chǎn)業(yè)空間格局，如靠近城區(qū)的低丘緩坡地適宜開發(fā)為商住用地或生產(chǎn)用地，臨近工業(yè)園區(qū)的低丘緩坡地應(yīng)因地制宜開發(fā)為工業(yè)基地，為人們提供更多的空間發(fā)展。

表3 邵武市各鄉(xiāng)鎮(zhèn)低丘緩坡地適宜建設(shè)統(tǒng)計分析

較適宜建設(shè)區(qū)主要分布在邵武北部的水北鎮(zhèn)，東南部的拿口鎮(zhèn)、吳家塘鎮(zhèn)，城區(qū)的城郊鎮(zhèn)，西部的沿山鎮(zhèn)，西南部的大埠崗鄉(xiāng)、肖家坊鎮(zhèn)、和平鎮(zhèn)。其中，水北鎮(zhèn)北部海拔較高，且森林覆蓋率較大；拿口鎮(zhèn)人口密集，距市區(qū)較近，距地質(zhì)災(zāi)害易發(fā)區(qū)較遠(yuǎn)，且不易發(fā)生地震；大埠崗鎮(zhèn)地勢由東向西逐漸傾斜，有名勝古跡位于其中；肖家坊鎮(zhèn)地處武夷山脈南麓，距離市區(qū)較遠(yuǎn)，境內(nèi)有將石自然保護區(qū)、雞公山、天成巖，可以就地取材、因地制宜，發(fā)展旅游業(yè)；和平鎮(zhèn)全境屬于丘陵地帶，境內(nèi)山區(qū)、半山區(qū)、平原各占1/3，平均海拔320 m，距離市區(qū)較遠(yuǎn)，交通便利，是一個古鎮(zhèn)。這些區(qū)域可以作為建設(shè)用地的備選資源，進行低強度、低密度開發(fā)，可以因地制宜發(fā)展旅游業(yè)，比如適當(dāng)開發(fā)些農(nóng)家樂用地，景區(qū)內(nèi)部和附近的低丘緩坡地可以建設(shè)些瀏覽和配套服務(wù)設(shè)施，提高景區(qū)的服務(wù)質(zhì)量，這樣不僅優(yōu)化城市結(jié)構(gòu)，還會拉動鄉(xiāng)鎮(zhèn)居民經(jīng)濟收入，提高居民的生活水平。

一般適宜建設(shè)區(qū)零星分布在邵武北部的水北鎮(zhèn)，東南部的拿口鎮(zhèn)，西南部的沿山鎮(zhèn)，南部的和平鎮(zhèn)、肖家坊鎮(zhèn)、大埠崗鄉(xiāng)、洪墩鎮(zhèn)以及衛(wèi)閩鎮(zhèn)，這些鄉(xiāng)鎮(zhèn)人口密度、鄉(xiāng)鎮(zhèn)收入水平相對較低，海拔較高、坡度較大、森林覆蓋率高達70%，開發(fā)建設(shè)用地需要很高的技術(shù)和成本，應(yīng)該以生態(tài)保護為原則，加強生態(tài)景觀建設(shè)同時防止地質(zhì)災(zāi)害的基礎(chǔ)設(shè)施建設(shè)，盡可能少地破壞原始地形。

不適宜建設(shè)區(qū)主要分布在拿口鎮(zhèn)、沿山鎮(zhèn)、水北鎮(zhèn)、洪墩鎮(zhèn)、肖家坊鎮(zhèn)、和平鎮(zhèn)、大埠崗鄉(xiāng)、大竹鎮(zhèn)等，其中沿山鎮(zhèn)、洪墩鎮(zhèn)森林資源豐富，覆蓋率達76%，是邵武市林業(yè)重鎮(zhèn)；大竹鎮(zhèn)有生態(tài)公益林2 800 hm2，竹山面積 1 466.67 hm2，屬于重點保護林業(yè)帶。這些鄉(xiāng)鎮(zhèn)如果進行開發(fā)建設(shè)不僅會導(dǎo)致地質(zhì)災(zāi)害的發(fā)生，還會破壞生態(tài)環(huán)境。建議禁止在該區(qū)域進行建設(shè)用地的開發(fā)建設(shè)，應(yīng)該以生態(tài)景觀恢復(fù)和生物多樣性保護為主，盡可能保持原始風(fēng)貌，實現(xiàn)自然與社會和諧、友好發(fā)展。

4 結(jié)論

本研究結(jié)合邵武市獨特的地形地貌，從空間區(qū)位、自然地理、社會經(jīng)濟、生態(tài)環(huán)境4個方面全面構(gòu)建評價指標(biāo)體系，采用主成分分析和均值聚類初步評價邵武市低丘緩坡地空間拓展的適宜性。該方法不僅能最大比例地保留原始信息，還能客觀確定權(quán)重，避免主觀隨意性。研究結(jié)果表明適宜建設(shè)區(qū)主要分布在主城區(qū)附近，涉及的鄉(xiāng)(鎮(zhèn))主要有城郊鎮(zhèn)、曬口街道、水北鎮(zhèn)、吳家塘鎮(zhèn)，占低丘緩坡地的15.88%，占邵武市總面積的6.64%，與規(guī)劃文本相吻合。該評價結(jié)果不但可以為邵武市低丘緩坡地空間拓展的開發(fā)潛力和方向提供宏觀把握，還可以為各鄉(xiāng)(鎮(zhèn))區(qū)域發(fā)展政策的制定提供可行性建議。

[1]古 艷，余萬軍. 淺析GIS技術(shù)在低丘緩坡建設(shè)適宜性評價中的應(yīng)用[J]. 商業(yè)經(jīng)濟，2017(5)：111-114.

[2]羅 鵬，俞志強，殷國慶，等. 基于3S技術(shù)的金衢地區(qū)低丘紅壤開發(fā)利用潛力評價[J]. 浙江農(nóng)業(yè)學(xué)報，2007，19(2)：101-105.

[3]朱曉蕓. 低丘緩坡土地資源開發(fā)利用評價研究[D]. 杭州：浙江大學(xué)，2008.

[4]徐 萍，衛(wèi) 新，王美青，等. 浙江省低丘緩坡農(nóng)業(yè)資源高效集約利用的路徑與對策研究[J]. 中國農(nóng)業(yè)資源與區(qū)劃，2013，34(3)：73-77.

[5]Keeney S，Hasson F，McKenna H P.A critical review of the delphi technique as a research methodology for nursing[J]. International Journal of Nursing Studies，2012，38(2)：195-200.

[6]Lee Y，Kozar K A.Investigating the effect of website quality on e-business success：an analytic hierarchy process (AHP) approach[J]. Decision Support Systems，2006，42(3)：1383-1401.

[7]袁 泉，黃義忠，謝小棠，等. 云南省羅平縣基于GIS低丘緩坡建設(shè)用地適宜性評價的研究[J]. 浙江農(nóng)業(yè)科學(xué)，2014(3)：408-411.

[8]Stacklies W，Redestig H，Scholz M，et al.pcaMethods-a bioconductor package providing PCA methods for incomplete data[J]. Bioinformatics，2007，23(9)：1164-1167.

[9]Hartigan J A，Wong M A. AK-means clustering algorithm[J]. Applied Statistics，1979，28(1)：100-108.

[10]葉宗裕. 關(guān)于多指標(biāo)綜合評價中指標(biāo)正向化和無量綱化法的選擇[J]. 浙江統(tǒng)計，2003(4)：24-25.

[11]汪 麗. 白銀市低丘緩坡土地建設(shè)開發(fā)適宜性評價研究[D]. 蘭州：甘肅農(nóng)業(yè)大學(xué)，2013.

[12]孫曉莉. 基于GIS的低丘緩坡建設(shè)用地適宜性評價[J]. 貴州大學(xué)學(xué)報(自然科學(xué)版)，2013，30(2)：131-134.

[13]鮑 艷，胡振琪，柏 玉，等. 主成分聚類分析在土地利用生態(tài)安全評價中的應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報，2006，22(8)：87-90.

[14]叢明珠，歐向軍，趙 清，等. 基于主成分分析法的江蘇省土地利用綜合分區(qū)研究[J]. 地理研究，2008，27(3)：574-582.

[15]張雪波，肖世青，杜先鋒，等. 基于主成分分析法的安溪鐵觀音香氣質(zhì)量評價模型的構(gòu)建[J]. 食品科學(xué)，2012，33(22)：225-230.

[16]公麗艷，孟憲軍，劉乃僑，等. 基于主成分與聚類分析的蘋果加工品質(zhì)評價[J]. 農(nóng)業(yè)工程學(xué)報，2014，30(13)：276-285.