基于主成分聚類分析的山東省土地生態(tài)安全評價及其影響因素分析

汪磊+張覺文

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

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

中圖分類號： F301.2文獻(xiàn)標(biāo)志碼： A文章編號：1002-1302（2017）17-0246-04

收稿日期：2017-05-17

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

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

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

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

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

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

2指標(biāo)體系構(gòu)建

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

該指標(biāo)體系包括的指標(biāo)有地均GDP（X1）、第三產(chǎn)業(yè)增加值（X2）、第三產(chǎn)業(yè)占比（X3）、外商直接投資（X4）、單位GDP能耗（X5）、森林覆蓋率（X6）、濕地面積占行政面積比重（X7）、建設(shè)用地占比（X8）。其中，地均GDP是一個復(fù)合指標(biāo)，能反映土地的使用效率，該指標(biāo)值越高，意味著經(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è)轉(zhuǎn)變；單位GDP能耗表征一個國家經(jīng)濟活動中對能源的利用程度，經(jīng)濟結(jié)構(gòu)足夠合理、能源利用效率高有利于提高土地集約利用程度；森林覆蓋率和濕地面積占行政面積比重2項指標(biāo)則從生態(tài)的角度直接體現(xiàn)土地生態(tài)安全的狀況；建設(shè)用地占比能更直觀地反映區(qū)域范圍內(nèi)土地利用結(jié)構(gòu)是否合理。表1山東省土地生態(tài)安全評價指標(biāo)體系

目標(biāo)層指標(biāo)層指標(biāo)含義指標(biāo)單位指標(biāo)性質(zhì)指標(biāo)代碼土地生態(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)政府有關(guān)部門批準(zhǔn)的項目投資總額內(nèi)企業(yè)從境外借入的資金萬美元+X4單位GDP能耗指一定時期內(nèi)一個國家（地區(qū)）每生產(chǎn)一個單位的國內(nèi)（地區(qū)）生產(chǎn)總值所消耗的能源t標(biāo)準(zhǔn)煤/萬元-X5森林覆蓋率指一個國家或地區(qū)森林面積占土地面積的百分比%+X6濕地面積占行政面積比重指一個國家或地區(qū)濕地面積占行政面積的百分比%+X7建設(shè)用地占比指一個國家或地區(qū)建設(shè)用地占土地面積的百分比%-X8注：“+”表示效益型指標(biāo)；“-”表示成本型指標(biāo)。

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

3.1主成分聚類分析

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

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

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

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

3.2內(nèi)涵分析

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

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

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

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

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

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

4結(jié)論

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

參考文獻(xiàn)：

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

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

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

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

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

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

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