汪磊+張覺文

摘要：土地生態(tài)安全關系著區(qū)域生態(tài)環(huán)境可持續(xù)發(fā)展及社會經(jīng)濟系統(tǒng)的良性運行，科學評價土地生態(tài)安全是生態(tài)文明建設的重要內(nèi)容。本研究基于經(jīng)濟、社會、生態(tài)3個維度，構建土地生態(tài)安全評價指標體系，借助主成分聚類組合評價模型，對2015年山東省土地生態(tài)安全狀況進行綜合評價，并結合評價結果進行區(qū)域等級劃分。結果表明，主成分聚類分析結果與實際情況高度吻合；山東省土地生態(tài)安全空間分布極不均衡，風險區(qū)劃之間差異性大，其中土地生態(tài)安全低度風險區(qū)（第Ⅰ類區(qū)域）2個（濟南市、青島市）、中度風險區(qū)（第Ⅱ類區(qū)域）14個、高度風險區(qū)（第Ⅲ類區(qū)域）1個（東營市）。今后，山東省除確保Ⅰ類區(qū)域2座城市的生態(tài)安全水平外，應將土地生態(tài)保護工作的重心放在第Ⅱ類區(qū)域，謹防其滑入第Ⅲ類區(qū)域，同時也要積極投入相應的人力、物力、財力，制定科學有效的政策措施以提高東營市的土地生態(tài)安全水平。

關鍵詞：土地生態(tài)安全；綜合評價；主成分聚類分析；山東??；生態(tài)文明；可持續(xù)發(fā)展

中圖分類號： F301.2文獻標志碼： A文章編號：1002-1302（2017）17-0246-04

收稿日期：2017-05-17

基金項目：國家社會科學基金（編號：17BJY011）。

作者簡介：汪磊（1974—），男，湖北黃岡人，博士，教授，主要從事土地資源管理及生態(tài)治理研究。E-mail：80290257@qq.com。

通信作者：張覺文，碩士研究生，研究方向為土地資源管理。E-mail：1040461115@qq.com。土地資源是人類賴以生存和發(fā)展的重要物質(zhì)基礎。土地生態(tài)安全關系著土地資源的可持續(xù)利用，維系著經(jīng)濟安全以及社會穩(wěn)定，是生態(tài)文明建設的重要內(nèi)容。隨著我國城市化、工業(yè)化進程的不斷加快，土地資源的粗放式開發(fā)引發(fā)了諸如水土流失加劇、耕地質(zhì)量下降、草地退化嚴重、生物多樣性減少等系列問題。近年來，山東省經(jīng)濟發(fā)展勢頭強勁，城鎮(zhèn)化速度迅猛，2015年常住人口城鎮(zhèn)化率已達57.0%，生態(tài)用地逐漸減少，人地矛盾越發(fā)突出，土地生態(tài)功能逐步衰退。統(tǒng)計數(shù)據(jù)表明，截至2015年年底，山東省水土流失面積高達 22 160 km2；國家林業(yè)局抽查結果表明，山東省經(jīng)濟林土壤中均不同程度檢測出重金屬。土地生態(tài)安全問題已成為制約山東省經(jīng)濟社會可持續(xù)發(fā)展的重要因素。

當前，土地生態(tài)安全問題日趨嚴峻，對土地生態(tài)安全的研究具有很強的重要性和緊迫性，引起了眾多學者的關注。國內(nèi)對土地生態(tài)安全的研究自2006年后不斷升溫，雖起步較晚但發(fā)展迅速。其中，學術界對土地生態(tài)安全的研究熱點集中于土地生態(tài)安全評價。從評價指標體系看，目前國內(nèi)對土地生態(tài)安全的研究，大多基于壓力-狀態(tài)-響應（pressure-state-response，簡稱PSR）和環(huán)境-經(jīng)濟-社會（environment-economy-society，簡稱EES）及其修正模型2種基本體系構建。如王磊等基于PSR框架選取了30項指標因子構建土地生態(tài)安全評價指標體系，系統(tǒng)評價了寧夏回族自治區(qū)的土地生態(tài)安全狀況[1]；莊偉等以EES為三維度選取了16項相關指標構建了三峽庫區(qū)土地生態(tài)安全評估指標體系并對其進行評價[2]。從評價方法來看，基于數(shù)學模型法的綜合指數(shù)分析法在土地生態(tài)安全研究中應用范圍最為廣泛。該方法的關鍵在于確定各評價指標的權重。比較成熟的賦權方法有熵值法[3]、層次分析法[4]、主成分分析法[5]、逼近于理想解的排序技術（technique for order preference by similarity to an ideal solution，簡稱TOPSIS法）[1]、德爾菲法[6]等，經(jīng)文獻梳理可知，關于土地生態(tài)安全評價大多采用單一的評價方法，忽略了單一評價方法存在的自身缺陷。如層次分析法、德爾菲法等賦權方法主觀性強；主成分分析法給定的權重與數(shù)據(jù)自身的質(zhì)量高度相關；TOPSIS法則將各個評價指標的權重均等化。為克服單一評價方法存在的不足，本研究通過將主成分分析與聚類分析2種方法進行有效組合，實現(xiàn)二者之間的優(yōu)勢互補、取長補短，極大地提高了土地生態(tài)安全評價的精確性和有效性，即將主成分分析的結果（提煉出的主成分）作為聚類分析的輸入，一方面通過主成分分析實現(xiàn)有效降維、消除提煉出主成分之間的線性相關；另一方面，借助聚類分析可以從更大的空間尺度分析不同類別之間的差異性，有助于提高制定政策的針對性和有效性。從應用的角度來看，根據(jù)實證分析得出的主成分聚類分析結果，一方面可以從宏觀層面分析不同類別之間的空間差異；另一方面，可以從微觀層面分析每個地市在不同主成分得分上的不足，可以為提升山東省土地生態(tài)安全水平，科學制定相關政策措施提供決策參考。

1樣本數(shù)據(jù)來源

本研究所構建的山東省土地生態(tài)安全評價指標體系所需的截面數(shù)據(jù)來源于山東省各政府職能部門公布的官方公報或年鑒：《2016山東省統(tǒng)計年鑒》《2015年度山東省林業(yè)資源概況與發(fā)展目標》《2015年山東省環(huán)境狀況公報》等。

2指標體系構建

土地生態(tài)安全評價包括指標體系構建和評價方法這2個方面的內(nèi)容。其中，構建山東省土地生態(tài)安全評價指標體系是實施評價過程的第一步。當前，既有文獻大多從人類與環(huán)境之間的關系出發(fā)，基于壓力、狀態(tài)、響應3個維度（PSR）構建相應的指標體系。實際上，影響土地生態(tài)安全的因素多、機制復雜。土地生態(tài)系統(tǒng)是以土地為載體，受自然環(huán)境、經(jīng)濟發(fā)展、社會發(fā)展3大因素影響而形成的有機整體，較少的學者從環(huán)境、經(jīng)濟、社會3個方面來構建土地生態(tài)安全的評價指標體系。根據(jù)指標的可比性、數(shù)據(jù)的易獲取性、指標的合理性、體系的簡潔性等原則，從《2016山東省統(tǒng)計年鑒》《2015年度山東省林業(yè)資源概況與發(fā)展目標》《2015年山東省環(huán)境狀況公報》等官方權威資料中提取相應的指標，構建山東省土地生態(tài)安全水平的指標體系（表1），并運用該指標體系對山東省17個市進行實證分析。endprint

該指標體系包括的指標有地均GDP（X1）、第三產(chǎn)業(yè)增加值（X2）、第三產(chǎn)業(yè)占比（X3）、外商直接投資（X4）、單位GDP能耗（X5）、森林覆蓋率（X6）、濕地面積占行政面積比重（X7）、建設用地占比（X8）。其中，地均GDP是一個復合指標，能反映土地的使用效率，該指標值越高，意味著經(jīng)濟發(fā)展空間上更為集中，理論上用于生態(tài)層面的土地則相對較多，生態(tài)安全水平會更高；第三產(chǎn)業(yè)是指不生產(chǎn)物質(zhì)產(chǎn)品的行業(yè)，與第一、第二產(chǎn)業(yè)相比，第三產(chǎn)業(yè)對土地資源依賴程度較低。因此，第三產(chǎn)業(yè)增加值越大、占比越高可以對土地生態(tài)安全帶來積極影響。實踐中，近年來山東省在引入外資項目中逐步實現(xiàn)了由以工業(yè)為主的第二產(chǎn)業(yè)向以金融、投資、教育為主的第三產(chǎn)業(yè)轉變；單位GDP能耗表征一個國家經(jīng)濟活動中對能源的利用程度，經(jīng)濟結構足夠合理、能源利用效率高有利于提高土地集約利用程度；森林覆蓋率和濕地面積占行政面積比重2項指標則從生態(tài)的角度直接體現(xiàn)土地生態(tài)安全的狀況；建設用地占比能更直觀地反映區(qū)域范圍內(nèi)土地利用結構是否合理。表1山東省土地生態(tài)安全評價指標體系

目標層指標層指標含義指標單位指標性質(zhì)指標代碼土地生態(tài)安全地均GDP指每km2土地創(chuàng)造的GDP億元/km2+X1第三產(chǎn)業(yè)增加值指在一定時期內(nèi)第三產(chǎn)業(yè)單位產(chǎn)值的增加值億元+X2第三產(chǎn)業(yè)占比指第三產(chǎn)業(yè)占總產(chǎn)業(yè)產(chǎn)值的比重%+X3外商直接投資指外國企業(yè)和經(jīng)濟組織或個人在我國境內(nèi)開辦外商獨資企業(yè)、與我國境內(nèi)的企業(yè)或經(jīng)濟組織共同舉辦中外合資經(jīng)營企業(yè)、合作經(jīng)營企業(yè)或合作開發(fā)資源的投資，以及經(jīng)政府有關部門批準的項目投資總額內(nèi)企業(yè)從境外借入的資金萬美元+X4單位GDP能耗指一定時期內(nèi)一個國家（地區(qū)）每生產(chǎn)一個單位的國內(nèi)（地區(qū)）生產(chǎn)總值所消耗的能源t標準煤/萬元-X5森林覆蓋率指一個國家或地區(qū)森林面積占土地面積的百分比%+X6濕地面積占行政面積比重指一個國家或地區(qū)濕地面積占行政面積的百分比%+X7建設用地占比指一個國家或地區(qū)建設用地占土地面積的百分比%-X8注：“+”表示效益型指標；“-”表示成本型指標。

3山東省土地生態(tài)安全評價過程

3.1主成分聚類分析

3.1.1主成分分析過程本研究使用SPSS 20.0軟件作為分析工具，考慮不同指標是從不同維度對山東省土地生態(tài)安全進行刻畫，彼此間存在著量綱和數(shù)量級影響，因此事先需要對原始數(shù)據(jù)進行正向化處理后再實施標準化處理（表2），以消除量綱差異和數(shù)量級影響，需要說明的是標準化處理后相關指標值變?yōu)樨摂?shù)，表明該樣本的指標值低于全省該指標的平均水平。同時得到了兩兩指標之間的相關系數(shù)矩陣R（表3）。由相關系數(shù)矩陣可知，指標之間既存在正負相關也存在強弱相關，且部分指標之間相關性很強，如ZX2與ZX4之間的相關系數(shù)為0.88，說明指標之間所攜帶的信息有很大重疊，若直接用于分析，會帶來嚴重的多重共線性問題。需要說明的是，ZX1表示的是原始指標X1經(jīng)過標準化處理后對應的變量，表3中其他各指標的含義如此類推。

由表4可知，指標間存在較強的相關性，可用主成分分析進行降維。由軟件計算可得KMO（Kaiser-Meyer-Olkin）統(tǒng)計值為0.713，表明指標體系的整體相關程度很高，Bartlett球形檢驗顯著性P<0.000 1，表明樣本數(shù)量充足，相關系數(shù)矩陣R為非單位陣，故可以實施主成分分析。矩陣特征值與累計貢獻率如表4所示，提取前3個主成分即提取了樣本84062%的數(shù)據(jù)信息（前3個特征值的累計貢獻率為84062%）。提取的3個主成分的載荷矩陣如表5所示，其中F1、F2、F3分別表示表示第1、2、3個主成分。由表5可知，F(xiàn)1在ZX1（地均GDP）、ZX2（第三產(chǎn)業(yè)增加值）、ZX3（第三產(chǎn)業(yè)占比）、ZX4（外商直接投資）4個變量上的荷載值都很大，這些指標都是從經(jīng)濟層面刻畫經(jīng)濟發(fā)展對土地生態(tài)安全的影響，將其定義為經(jīng)濟發(fā)展質(zhì)量因子；F2在ZX6（森林覆蓋率（%）、ZX7（濕地面積占行政面積比重）2個指標上有較高的載荷，這些指標都是從生態(tài)環(huán)境層面反映生態(tài)環(huán)境因素對土地生態(tài)安全的影響，可以將其定義為生態(tài)環(huán)境發(fā)展因子；F3在ZX5（單位GDP能耗）、ZX8（建設用地占比）2個變量上的荷載值較高，由于當前山東省各地級市的經(jīng)濟發(fā)展主要依賴于土地財政驅(qū)動，從而導致占用耕地現(xiàn)象比較普遍，一定程度上會惡化土地生態(tài)安全，而ZX5（單位GDP能耗）從能耗角度間接地反映了土地集約利用程度，這些指標都反映了土地利用結構對土地生態(tài)安全的影響，定義為土地利用結構因子。

3.1.2聚類分析過程傳統(tǒng)的聚類分析由于無法消除評價指標之間的相關性而導致聚類結果誤差較大。為消除這一缺陷，本研究運用SPSS 20.0中的聚類分析程序模塊，以上述提煉的3個主成分F1、F2、F3的分析結果作為聚類分析的輸入（F1、F2、F3之間線性無關），對山東省17個市的土地生態(tài)安全情況進行聚類分析，得到山東省17個市土地生態(tài)安全聚類分析的譜系圖（圖1）。一方面，從實證分析的結果來看，聚類分析的類別結果與主成分分析的排序結果高度吻合，即在主成分分析結果中綜合得分排在前兩位的是青島和濟南，排在最后面一位的是東營，而在聚類分析中，如果按照3類來進行劃分，則聚在第一類的也是青島和濟南，聚在最后一類的是東營，聚類結果與主成分排序結果完全一致，從實證的角度表明本研究所構建的主成分聚類分析模型具有較高的科學合理性；另一方面，從山東省土地生態(tài)的實際情況來看，主成分聚類分析的結果與近年來山東省土地生態(tài)安全發(fā)展實際高度相符，且與已經(jīng)發(fā)表的關于山東省土地生態(tài)安全評價學術論文中的結果高度一致。

由上述分析可知，將山東省17個城市土地生態(tài)安全評價結果聚成3類比較合適。主成分綜合得分排在第1、第2位的青島和濟南聚為第Ⅰ類（綜合得分接近或高于1分）；東營的綜合得分最低（綜合得分低于-1分），且遠低于其他城市的土地生態(tài)安全綜合得分，劃為第Ⅲ類；其他城市如淄博、濰坊、棗莊、泰安、萊蕪、煙臺、威海、臨沂、德州、聊城、濱州、日照、濟寧、菏澤的綜合得分排在中間，聚為第Ⅱ類（綜合得分在-1～1分之間）。此外，從3個類別的城市數(shù)量分布來看，聚在Ⅰ類的有2個城市，聚在Ⅱ類的有14個城市，聚在Ⅲ類的有1個城市。3類城市聚類數(shù)量基本服從正態(tài)分布，具有較強的邏輯合理性，空間分布如圖2所示。endprint

3.2內(nèi)涵分析

由表6、圖1的對比分析可知，將山東省土地生態(tài)安全劃分為低度風險區(qū)、中度風險區(qū)、高度風險區(qū)3種類型。聚類分析的結果與綜合排序完全一致。其中，土地生態(tài)安全低度風險區(qū)（Ⅰ類）包括青島、濟南2個城市；土地生態(tài)安全中度風險區(qū)（Ⅱ類）包括淄博、濰坊、棗莊、泰安、萊蕪、煙臺、威海、臨沂、德州、聊城、濱州、日照、濟寧、菏澤14個城市；土地生態(tài)安全高度風險區(qū)（Ⅲ類）包括東營1個城市。

第Ⅰ類中，限于篇幅，僅以青島為例展開分析。從主成分得分來看，青島市F1得分（3.029）排在第1位，表明青島經(jīng)濟發(fā)展質(zhì)量很高，對土地生態(tài)環(huán)境的破壞程度低；F2得分大于0，表明青島生態(tài)環(huán)境發(fā)展因子在全省平均水平之上；土地利用結構因子F3得分略低于全省平均水平。實踐中，從經(jīng)濟發(fā)展狀況來看，青島經(jīng)濟總量大，GDP總量連年位列全省第一。獨特的自然環(huán)境、氣候條件和優(yōu)越的港口交通基礎設施使青島旅游觀光業(yè)和海洋產(chǎn)業(yè)對經(jīng)濟增長貢獻大，對土地的依賴程度??；從生態(tài)環(huán)境因素看，目前青島經(jīng)濟發(fā)展模式已經(jīng)實現(xiàn)了由投資驅(qū)動向創(chuàng)新驅(qū)動轉型。統(tǒng)計數(shù)據(jù)顯示，近年來發(fā)明專利申請量和授權量增長明顯，2015年達44 962件。經(jīng)濟發(fā)展過程中土地資源利用的效率明顯提高。此外，2016年青島設立59處省級生態(tài)紅線區(qū)，以改善和提高生態(tài)系統(tǒng)服務功能，減少水土流失；土地利用結構方面，青島市始終著力優(yōu)化空間布局，切實節(jié)約集約利用土地，嚴格保護耕地和基本農(nóng)田，有效保護海洋生態(tài)環(huán)境。因此，土地生態(tài)安全風險低。

第Ⅱ類中，城市數(shù)量較多，限于篇幅，這里僅以臨沂為例展開分析。由表6可知，臨沂市土地生態(tài)安全主成分F1（-0.444）和F3（-0.571）得分均低于平均分，且綜合得分（-0.17）也略低于全省平均值。目前臨沂市已成為省內(nèi)經(jīng)濟增長速度快、發(fā)展質(zhì)量高的城市。2016年綜合實力躍居全國百強城市榜第61位。其中，僅服務業(yè)對GDP的貢獻率已達到57.3%，經(jīng)濟增長對土地的依賴程度低。在發(fā)展經(jīng)濟的同時，臨沂市重視保護土地資源及生態(tài)平衡，造林面積和耕地面積均全省第一。但臨沂人口眾多，人口數(shù)位列山東省首位，2015年末總人口已達1 031.16萬人，導致人均GDP、人均森林面積、人均耕地面積等與土地生態(tài)安全有著直接關系的指標值不高，同時臨沂還是山東省水土流失最為嚴重的地區(qū)。綜合這些因素，臨沂未能進入Ⅰ類區(qū)。

第Ⅲ類屬于土地生態(tài)安全高度風險區(qū)，僅包括東營1個城市。由表6可知，F(xiàn)1、F2、F3的得分均為負值，遠低于全省平均水平。綜合得分F（-1.04）也是全省最低得分。作為新興的油氣工業(yè)城市，東營油田是中國的第二大油田，石化工業(yè)在區(qū)域經(jīng)濟中起主導作用，2015年東營第二產(chǎn)業(yè)生產(chǎn)總值占比達64.7%，增加值達2 230.61億元。然而東營經(jīng)濟快速增長的同時也污染了生態(tài)環(huán)境，大量而頻繁的油氣開采和化工產(chǎn)業(yè)生產(chǎn)也使得水資源、土壤、植被被污染，土壤的自凈化恢復能力也得不到恢復，土地生態(tài)安全面臨著極大威脅。因此，東營屬于土地生態(tài)安全高度風險區(qū)。

土地生態(tài)安全類別劃分有著重要的現(xiàn)實意義，可以為制定山東省土地資源的開發(fā)和利用提供決策參考。由圖2可知，山東省土地生態(tài)安全空間分布極不均衡，風險區(qū)劃之間差異性大，客觀上需要制定具有針對性的對策建議。結合表6中的聚類結果，山東省政府今后應該將土地生態(tài)保護工作的重心放在第Ⅱ類區(qū)域，積極提升第Ⅱ類城市生態(tài)安全水平，降低土地生態(tài)安全風險，謹防其滑入第Ⅲ類區(qū)域。同時要積極投入相應的人力、物力和財力，制定科學有效的政策措施，不

斷提高東營市的土地生態(tài)安全水平。

4結論

已有同類文獻的土地生態(tài)安全評價存在著選取指標數(shù)量多、數(shù)據(jù)收集困難、指標權重易被稀釋等諸多問題，而本研究僅采用少量具有典型性和代表性的指標（8項）進行分析，所得分析結果不僅與同類文獻中所得的結論高度一致，同時還與山東省土地生態(tài)安全的具體實際高度相符，表明本研究構建的指標體系具有科學合理性。此外，傳統(tǒng)土地生態(tài)評價大多使用單一評價方法展開，無法規(guī)避由單一評價方法帶來的自身缺陷。本研究通過構建基于主成分分析和聚類分析的組合模型，充分利用提煉出的主成分之間線性無關的優(yōu)良特性，將該主成分分析結果作為聚類分析的數(shù)據(jù)輸入，在使用較少指標的前提下極大地提高了對土地生態(tài)安全評價的精度。結果表明，本研究使用聚類分析與主成分分析二者結合所得結果高度一致，同時還與其他文獻關于山東省土地生態(tài)安全水平的結論基本相符[7]，客觀上證實了該組合評價方法的科學合理性，可為其他省份開展土地生態(tài)安全評價提供方法上的借鑒。當然，土地生態(tài)安全水平是一個動態(tài)的發(fā)展過程，如何結合時間序列數(shù)據(jù)動態(tài)分析山東省土地生態(tài)安全的演化過程是今后的一個研究方向。

參考文獻：

[1]王磊，郭燦，李慧明. 基于PSR-TOPSIS模型的寧夏回族自治區(qū)土地生態(tài)安全評價[J]. 水土保持研究，2016，23（6）：154-159.

[2]莊偉，廖和平，潘卓，等. 基于變權TOPSIS模型的三峽庫區(qū)土地生態(tài)安全評估——以巫山縣為例[J]. 西南大學學報（自然科學版），2014，36（8）：106-112.

[3]李悅，張合兵，張小虎，等. 基于熵值法和灰色預測模型的城市土地生態(tài)安全評價[J]. 環(huán)境科學與技術，2015，12（38）：242-247.

[4]馮文斌，李升峰. 江蘇省土地生態(tài)安全評價研究[J]. 水土保持通報，2013，33（2）：285-290.

[5]榮聯(lián)偉，師學義，高奇，等. 黃土高原山丘區(qū)土地生態(tài)安全動態(tài)評價及預測[J]. 水土保持研究，2015，22（3）：210-216.

[6]張利，陳影，王樹濤，等. 濱?？焖俪鞘谢貐^(qū)土地生態(tài)安全評價與預警——以曹妃甸新區(qū)為例[J]. 應用生態(tài)學報，2015，26（8）：2445-2454.

[7]潘潤秋，孫煥芳，劉瑛慧. 基于TOPSIS模型的山東省土地生態(tài)安全評價[J]. 湖北農(nóng)業(yè)科學，2016，2（55）：509-514.蔣云峰，錢瑞雪，王祥賀，等. 免耕玉米秸稈覆蓋對大型土壤動物功能類群的影響[J]. 江蘇農(nóng)業(yè)科學，2017，45（17）：250-253.endprint