郝天平，張洋，王志龍

(河北工程大學 資源學院，河北 邯鄲056038)

土地節(jié)約集約評價是根據(jù)土地的投入和產(chǎn)出狀況，考慮土地的經(jīng)濟效益、社會效益、生態(tài)效益進行的綜合評價[1－2]。我國土地節(jié)約集約利用的研究開始于20世紀90年代，起步較晚，主要集中在土地節(jié)約集約利用的相關(guān)理論和內(nèi)涵、土地節(jié)約集約利用定量評價、土地節(jié)約集約利用途徑及區(qū)域?qū)嵺`研究等方面[3－5]。在定量評價方面施建剛等[6]利用因子分析和聚類分析法，趙旭等[7]利用聚類分析和文獻資料法對城市土地節(jié)約集約評價進行了研究，從總體上看這些方法并沒有涉及對評價指標的篩選，一些評價僅對結(jié)果進行了文字性說明，表現(xiàn)不直觀。本文運用主成分分析法(Principal Component Analysis，PCA)[8]對土地節(jié)約集約的各評價指標進行相關(guān)性分析，篩選出對評價影響較大的主要指標，在此基礎(chǔ)上利用層次分析法(Analytical Hierarchy Process，AHP)獲取各指標的權(quán)重值，利用GIS對土地節(jié)約集約評價結(jié)果進行制圖表達。

1 模型建立

1.1 評價指標篩選

表1 土地節(jié)約集約評價指標Tab.1 Evaluation index of land saving

表2 層次分析模型層次結(jié)構(gòu)Tab.2 Hierarchy structure of the AHP model

土地節(jié)約集約評價主要包括土地利用強度、土地投入強度、土地的經(jīng)濟效益和生態(tài)效益4個方面，具體指標如表1所示[10]，由于相同指標在不同的地區(qū)發(fā)揮的作用不相同，所以利用主成分分析對評價指標進行篩選，得到影響評價結(jié)果的主要指標。

將研究區(qū)分成若干小的評價區(qū)域，分別計算各區(qū)域的指標值，計算結(jié)果用矩陣X表示。

式中:Xij－第i小區(qū)域的第j指標(變量)的值，0 ＜i＜n，0＜j＜p。

對式(1)中的數(shù)據(jù)進行標準化處理，統(tǒng)一各評價指標的量綱和取值范圍，得到矩陣X'

對式(2)標準化處理得到相關(guān)矩陣R

對矩陣R的特征值λi(0＜i≤p)按大小排序，根據(jù)特征值得到R的特征向量ei。

取Fi達到0.85～0.95的特征值時，第1～ i個主成分所對應的指標便為篩選得到的評價指標。

1.2 權(quán)重確定

利用主成分析法對評價指標進行篩選，采用層次分析法建立層次結(jié)構(gòu)(表2)，根據(jù)Saaty等提出的1～9尺度(表3)，將評價指標量化處理。

表3 定性到定量的轉(zhuǎn)化表Tab.3 Convert qualitative into quantitative

步驟1根據(jù)準則層對目標層的重要性建立判斷矩陣A

式中:aij－準則層Ci與Cj比較對于目標層A的重要程度(貢獻大小)。

步驟2根據(jù)評價指標對準則層的重要程度建立判斷矩陣C(k)

步驟3分別對判斷矩陣A，C(k)處理，得到最大特征值及對應的特征向量 ea，e(k)。

利用式(7)判斷矩陣的一致性，當CR＜0.1時，判斷矩陣具有一致性，否則調(diào)整判斷矩陣。

式中:n－判斷矩陣的維數(shù)，RI通過查平均隨機一致性指標表獲得。

利用式(8)判斷總體一致性，同理當CR＜0.1時，說明判斷矩陣具有一致性，否則調(diào)整各判斷矩陣。

指標的權(quán)重計算公式為

土地節(jié)約集約評價結(jié)果f為

其中P表示評價指標標準化處理后組成的向量。

2 GIS分析

在主成分分析和層次分析建立的數(shù)量關(guān)系的基礎(chǔ)上，利用GIS的空間處理能力進行綜合分析[10－11]。GIS 分析的步驟如下:

步驟1計算指標值。對研究區(qū)按照一定規(guī)則劃分成多個區(qū)域，計算各區(qū)域不同評價指標標準化處理后的值。

步驟2矢量數(shù)據(jù)柵格化。即將多邊形矢量數(shù)據(jù)轉(zhuǎn)化為柵格數(shù)據(jù)，柵格的Value值為標準化處理后的指標值。不同評價指標生成不同柵格圖件。

步驟3代數(shù)運算。將評價指標柵格圖件Value值與對應的權(quán)值相乘作為柵格圖件的新的Value值。

步驟4疊加分析。將評價指標柵格圖件在GIS系統(tǒng)下進行疊加，取全部圖層Value值總和，所得到的圖像即為土地節(jié)約集約評價專題圖，其Value值為評價值，反映了該區(qū)域土地節(jié)約集約的平均情況。

3 案例分析

扎蘭屯市位于內(nèi)蒙古自治區(qū)東部，資源豐富，經(jīng)濟發(fā)展迅速，全市行政區(qū)域總面積1.69萬km2，總?cè)丝?5萬，人地矛盾相對比較緊張。

3.1 評價指標篩選和權(quán)重確定

利用主成分分析法對表1中的評價指標進行分析(表4)，確定該市的最終評價因子為人口密度、地均GDP、地均第二產(chǎn)業(yè)GDP、人均建設(shè)用地面積、地均固定資產(chǎn)投資額、地均第三產(chǎn)業(yè)GDP、人均道路面積等7個因子。

表4 特征值與貢獻率Tab.4 Characteristic values and contribution rate

在確定判斷矩陣的過程中，應用Matlab軟件對3位專家給出的判斷矩陣進行權(quán)重計算和一致性檢驗，結(jié)果均符合要求。通過對結(jié)果算術(shù)平均，得到該市建設(shè)用地節(jié)約集約評價指標矩陣為

各評價指標的權(quán)重為W=[0.198 0.134 0.091 0.096 0.137 0.062 0.060]。

3.2 空間分析與評價結(jié)果

將該市土地劃分為18個小區(qū)域，新增的7個屬性字段，分別對應標準化處理后的評價指標值。在ArcGIS下利用ArcToolBox的Polygon to Raster工具將18個區(qū)域的矢量圖轉(zhuǎn)化為柵格圖，其中Value field選擇評價指標，共計產(chǎn)生7個柵格圖件。用ArcMap打開7個柵格圖件，利用Spatial Analyst工具條下的Raster Calculator工具做代數(shù)運算，使柵格圖件的Value值為原值與該評價指標權(quán)重的乘積。利用同樣的工具對7個評價指標柵格圖求和，生成新的柵格圖層。新的柵格圖即為該市土地節(jié)約集約利用評價專題圖(圖1)，其value值便為該區(qū)域的評價值。

從圖1可知，顏色較深的2、5、6區(qū)域評價值較高，均達到了0.78以上，反應了這部分區(qū)域經(jīng)濟發(fā)展水平較高，土地節(jié)約集約程度較高;顏色越淺的區(qū)域評價值越低，土地節(jié)約集約利用程度越低;該市東北部為山地，土地節(jié)約集約利用程度最低，經(jīng)濟發(fā)展水平也是該市的最低。將評價結(jié)果與該市土地節(jié)約集約普查的結(jié)果進行對照，發(fā)現(xiàn)評價結(jié)果與普查結(jié)果是基本吻合的?？傊C合運用主成分分析和層次分析的評價反應了區(qū)域土地利用節(jié)約集約情況。該實例證明利用主成分分析確定評價指標，運用層次分析法獲取各評價指標權(quán)重值，最后發(fā)揮GIS的統(tǒng)計分析能力進行土地節(jié)約集約評價的方法是合理可行的。

4 結(jié)論

1)使用數(shù)學知識和專家經(jīng)驗模型相結(jié)合確定評價因子及其權(quán)重，不僅避免了數(shù)學方法的機械性，同時也避免了專家的主觀性。

2)在評價過程中引入了GIS技術(shù)，使評價的過程和結(jié)果數(shù)據(jù)都可以直觀的表現(xiàn)，減少了工作時間和工作量，節(jié)約了傳統(tǒng)評價方法需要的后期人工繪圖的人力、財力和物力。

3)通過主成分分析和層次分析，不僅體現(xiàn)了各評價指標的表面信息，更挖掘了評價指標間的關(guān)聯(lián)關(guān)系等深層次信息，使土地節(jié)約集約評價更科學，更具說服力。

[1]肖勁松.堅持節(jié)約集約用地，促進經(jīng)濟可持續(xù)發(fā)展[J].國土資源導刊，2006，3(2):35.

[2]魏國銳.土地節(jié)約集約利用面臨的問題集對策[J].山東國土資源，2006，22(2):57－59，62.

[3]王靜，邵曉梅.土地節(jié)約集約利用技術(shù)方法研究:現(xiàn)狀、問題與趨勢[J].地理科學進展，2008，27(3):68－74.

[4]劉敏.城市土地集約利用及其評價指標體系構(gòu)建[J].標準科學，2009，(9):38－40.

[5]周志剛.開發(fā)區(qū)土地集約利用評價關(guān)鍵技術(shù)集成創(chuàng)新研究[D].杭州:浙江大學，2010.

[6]施建剛，宋英瀟.城市土地節(jié)約集約利用時空評價研究－以上海市為例[J].土地市場，2010，8(356):43－46.

[7]趙旭，劉小平，薛 姝.兩型社會視角下湖南城市土地節(jié)約集約利用研究[J].城市發(fā)展研究，2010，17(5):102－107.

[8]范杰.主成分分析法的數(shù)值實現(xiàn)算法[J].河北工程大學學報:自然科學版，2007，24(4):103－105.

[9]馮 帆，朱 剛，高會軍.晉陜蒙接壤地區(qū)土地利用/覆蓋變化及其驅(qū)動力研究——以神華集團大柳塔礦區(qū)為例[J].中國煤炭地質(zhì)，2011，23(3):24 －26.

[10]楊大兵，李俊付，王 鳳.基于GIS的城鎮(zhèn)土地潛力評價系統(tǒng)研究[J].河北工程大學學報:自然科學版，2009，26(3):89－92.

[11]馮帆，王佟.基于GIS煤炭勘查地測空間信息系統(tǒng)[J].中國煤炭地質(zhì)，2011，23(4):64 －67.