侯小宇，何 剛，王曉東，谷雅嫻

(安徽理工大學(xué) 經(jīng)濟與管理學(xué)院，安徽 淮南 232001)

生態(tài)安全就是指一個國家或地區(qū)生態(tài)環(huán)境資源狀況不受或少受來自于資源和生態(tài)環(huán)境的制約與威脅的狀態(tài)[1].隨著經(jīng)濟的飛速發(fā)展，人類對生態(tài)資源的利用越來越廣泛、強度越來越大，尤其是土地生態(tài)資源.土地生態(tài)安全綜合評價能夠用于評估人類活動對土地生態(tài)環(huán)境的影響及土地生態(tài)環(huán)境保障人類正常生產(chǎn)生活的底線，并對土地生態(tài)環(huán)境功能結(jié)構(gòu)健康狀況進行判斷[2]，因此對省市域的土地生態(tài)安全進行評價是非常有意義的.

不同的研究區(qū)域、評價指標(biāo)體系、土地生態(tài)的測度方法，造成了土地生態(tài)安全領(lǐng)域研究的“百花齊放”.從土地生態(tài)安全研究尺度來看，可將其分為時間研究尺度和空間研究尺度[3].時間研究尺度有靜態(tài)和動態(tài)兩種，單個時間的土地生態(tài)研究是靜態(tài)的，如王一山[4]等人以烏魯木齊為研究區(qū)，研究其2017年土地生態(tài)安全現(xiàn)狀；動態(tài)土地生態(tài)安全研究按時間序列變化展開分析，如何如海[5]等人研究2008～2017年的淮河經(jīng)濟帶各個地級市的土地生態(tài)安全等級及其變化.從土地生態(tài)安全評價指標(biāo)上來看，張鳳太[6]以物元分析-DPSIR模型構(gòu)建重慶市土地生態(tài)安全評價模型；鄭華偉[7]基PSR模型與集對分析對耕地生態(tài)進行安全診斷；除此之外，還有自然-經(jīng)濟-社會(NES)、經(jīng)濟-環(huán)境-社會(EES)[8]等模型應(yīng)用到土地生態(tài)安全評價.從測度方法上來看，陳慧[9]利用系統(tǒng)動力學(xué)分析區(qū)域社會、經(jīng)濟、自然之間的關(guān)系，并進行仿真模擬研究土地生態(tài)安全；李子君[10]等人運用熵權(quán)物元方法構(gòu)建土地生態(tài)安全評判模型，對沂蒙山區(qū)土地生態(tài)安全的時空變化及其影響因素進行實證研究.

本文以EES模型構(gòu)建土地生態(tài)安全評價指標(biāo)體系，利用SPSS25對安徽省各個地級市進行主成分分析和聚類分析，最后依據(jù)綜合得分，對安徽省16個地級市土地生態(tài)安全進行排名和等級劃分.

1 研究區(qū)概況

本文以安徽省為研究區(qū).安徽省是中國內(nèi)陸城市，處于長江、淮河中下游.地跨長江、淮河南北，故安徽省南北地區(qū)土地作物差距較大.其與蘇、浙、鄂、豫、魯五省接壤，土地面積約14.01×104km2，為全國總面積的1.45%.安徽省共16個地級市，不同地級市文化差異、氣候差異較大，常以皖南、皖中和皖北地區(qū)劃分.2019年末，安徽省常住人口6 365.9萬人，常住人口城鎮(zhèn)化率55.81%，GDP達到37 113.98億元，較2018年增長8.36%.安徽省處于全國經(jīng)濟發(fā)展的戰(zhàn)略要沖，擁有兩大經(jīng)濟帶，即淮河生態(tài)經(jīng)濟帶、長江經(jīng)濟帶，部分市區(qū)享受了“西氣東輸”戰(zhàn)略的福利，故安徽省是國家重大戰(zhàn)略發(fā)展區(qū)域.然而由于工業(yè)化、新型城鎮(zhèn)化的持續(xù)推進，對安徽省土地利用與生態(tài)安全提出了新的要求.

2 數(shù)據(jù)來源與評價方法

2.1 數(shù)據(jù)來源

數(shù)據(jù)來源于《2020年安徽省統(tǒng)計年鑒》《2020年國民經(jīng)濟與社會發(fā)展統(tǒng)計公報》.

2.2 評價方法

2.2.1 主成分分析法

主成分分析的數(shù)學(xué)模型為：

(1)

由式(1)可得

(2)

(3)

(4)

(5)

(6)

4)var依次遞減：

(7)

其中：n為樣本數(shù)；p為樣本的觀測指標(biāo)數(shù)；X為隨機變量；μ為隨機變量均值；δ為隨機變量的協(xié)方差矩陣；y1,y2,…,yp為生成的新的綜合指標(biāo)即主成分.

本文利用SPSS 25對標(biāo)準(zhǔn)化后的指標(biāo)進行主成分分析，實現(xiàn)有效降維、消除指標(biāo)間線性相關(guān)性，得到相對獨立、重要的、少量的綜合指標(biāo).利用各個主成分得分，計算安徽省16個地級市的綜合指標(biāo)得分，對其進行土地生態(tài)安全等級劃分和排名.

2.2.2 聚類分析

為了消除主成分分析中指標(biāo)之間的差異性、更大空間尺度分析指標(biāo)的不同，本文采用聚類分析方法中的系統(tǒng)聚類法進行進一步的驗證[11].系統(tǒng)聚類分析(Hierarchical Cluster Method)包含以下步驟：

1)構(gòu)造n個類，其中只包含一個樣品.

2)計算每個類之間的距離，本文采用中間距離公式，得到距離矩陣，類之間的距離公式如下

(8)

3)將距離最近的兩類合并成一個新類.

4)計算合并后的新類與其他類之間的距離.

5)畫出聚類圖，并根據(jù)聚類圖決定類的個數(shù).

本文以主成分分析的輸出作為聚類分析的輸入，進一步驗證對安徽省16個地級市土地生態(tài)安全等級劃分及排名的準(zhǔn)確性.

3 評價指標(biāo)體系構(gòu)建

研究土地生態(tài)安全關(guān)鍵就是評價指標(biāo)的選擇，選取的指標(biāo)要有合理性、科學(xué)性，要能覆蓋各個地區(qū)的土地生態(tài)特征，能充分反映各個地區(qū)土地生態(tài)健康狀況.故從指標(biāo)體系構(gòu)建的科學(xué)性、系統(tǒng)性原則出發(fā)，參照相關(guān)研究[12]，立足于安徽省土地生態(tài)安全現(xiàn)狀，以“經(jīng)濟-環(huán)境-社會”(EES)模型為基礎(chǔ)構(gòu)建土地生態(tài)安全評價指標(biāo)體系.其中經(jīng)濟因素有四個正向指標(biāo)，即經(jīng)濟密度(萬元/km2)、地區(qū)生產(chǎn)總值(億元)、人均GDP(元/人)、第三產(chǎn)業(yè)占GDP比重(%)，它們反映了研究區(qū)經(jīng)濟發(fā)展水平和速度；環(huán)境因素有兩個正向指標(biāo)、兩個負(fù)向指標(biāo)，即農(nóng)藥施用強度(t/hm2)、人均道路面積(m2/人)、單位面積播種糧食產(chǎn)量(t/hm2)、工業(yè)固體廢物綜合利用率(%)，它們反映出對土地的利用程度，與土地生態(tài)健康有著直接關(guān)系；社會因素有兩個正向指標(biāo)、兩個負(fù)向指標(biāo)，即城鎮(zhèn)化水平(%)、建成區(qū)面積(km2)、農(nóng)業(yè)機械總動力(104kW)、人口自然增長率(%)，它們通過城鎮(zhèn)化、建成區(qū)、機械動力和人口的增加反映出社會發(fā)展對土地的影響.共選取了12個多維指標(biāo)，如表1所示.

表1 安徽省土地生態(tài)評價指標(biāo)體系

4 安徽省土地生態(tài)安全評價過程

4.1 原始指標(biāo)標(biāo)準(zhǔn)化

為了消除不同指標(biāo)在單位、數(shù)量級上的差距，需對指標(biāo)進行標(biāo)準(zhǔn)化處理[13].本文考慮到指標(biāo)體系的特性，利用SPSS25統(tǒng)計軟件下的描述統(tǒng)計，對安徽省16個地級市12項指標(biāo)進行z-score標(biāo)準(zhǔn)化.公式如下：

(9)

其中:Z為標(biāo)準(zhǔn)化后的輸出結(jié)果；C為具體指標(biāo)數(shù)值；μ為樣本平均數(shù)；σ為樣本標(biāo)準(zhǔn)差.標(biāo)準(zhǔn)化輸出的結(jié)果即Z值代表著原始指標(biāo)數(shù)值與總數(shù)值平均數(shù)之間的距離，固有正負(fù)之分.

4.2 適宜性檢驗

在原始數(shù)據(jù)Z-score標(biāo)準(zhǔn)化之后，需要對指標(biāo)之間的相關(guān)性進行檢驗.指標(biāo)之間的相關(guān)性決定著主成分分析提取的主成分是否有效，故需進行因子的適宜性檢驗[14].KMO與巴特利特球形檢驗用于比較變量指標(biāo)間的簡單相關(guān)系數(shù)和偏相關(guān)系數(shù)，其中KMO數(shù)值越接近1，代表著指標(biāo)變量之間的相關(guān)性越強，越適合做因子分析；巴特利特球形檢驗值是常用的KMO度量標(biāo)準(zhǔn)，0.9以上表示非常合適，0.5以下表示極不合適.SPSS25輸出的結(jié)果如表2，KMO檢驗統(tǒng)計量數(shù)值為0.648，在相關(guān)系數(shù)矩陣中有64.8%的系數(shù)個數(shù)大于0.3，這表明指標(biāo)體系之間具有很強的相關(guān)性，適合做因子分析；巴特利特球形度檢驗顯著性為0.000，說明原始數(shù)據(jù)服從正態(tài)分布總體，具有進一步研究的價值.

表2 KMO與巴特利特檢驗

4.3 指標(biāo)特征值與貢獻率的確定

指標(biāo)特征值與貢獻率結(jié)果如圖1、表3所示.結(jié)果顯示:到第四個主成分時，曲線開始趨于平穩(wěn)，特征值達到1，故在所有指標(biāo)中提取出3個主成分因子最為合適.前三個主成分因子的初始方差貢獻率分別為42.209%、24.376%、11.044%，累積方差貢獻率達到77.629%，即提取出來的3個主成分因子能夠涵蓋絕大部分的具有相關(guān)性的指標(biāo)體系信息，避免了指標(biāo)之間具有線性相關(guān)關(guān)系或者指標(biāo)與研究結(jié)果相關(guān)程度不高，從而導(dǎo)致研究結(jié)果不準(zhǔn)確等問題.

圖1 碎石圖

表3 各個主成分的總方差解釋和方差貢獻

4.4 主成分因子的確定

因子載荷矩陣是提取的主成分命名的主要依據(jù)，由于本研究指標(biāo)未進行因子旋轉(zhuǎn)的因子載荷矩陣數(shù)值較小，無法準(zhǔn)確、直觀的看出提取三個主成分代表的含義，故以四次方最大旋轉(zhuǎn)法(Quartimax)進行因子旋轉(zhuǎn)，得出因子的載荷矩陣，如表4.

表4 因子載荷矩陣

F1、F2、F3分別表示從12個相關(guān)指標(biāo)中提取的第一、第二、第三個主成分.由表4可知F1在經(jīng)濟密度、地區(qū)生產(chǎn)總值、人均GDP上的載荷量比較高，這些指標(biāo)代表著地區(qū)的經(jīng)濟發(fā)展水平，故將F1稱為經(jīng)濟發(fā)展質(zhì)量因子；F2在農(nóng)業(yè)機械總動力、人口自然增長率、人均道路面積上的載荷量比較大，其中農(nóng)業(yè)機械總動力占據(jù)了明顯的優(yōu)勢，側(cè)面反映了土地集約利用程度，故將F2稱為土地利用結(jié)構(gòu)因子；F3在第三產(chǎn)業(yè)占GDP比重、建成區(qū)面積、人口自然增長率上的載荷量比較高，既有經(jīng)濟因素又有社會因素，故將F3稱為人口發(fā)展?fàn)顟B(tài)因子.該評價指標(biāo)體系綜合12個指標(biāo)全面的反應(yīng)了F1、F2、F3以及土地生態(tài)安全綜合指數(shù)F三個方面的信息，能夠有效地對安徽省土地生態(tài)安全狀況作出評價.

4.5 主成分結(jié)果分析

4.5.1 土地生態(tài)安全狀況綜合得分

在指標(biāo)信息量保持完整的前提下，將安徽省16個地級市12個不同維度的指標(biāo)數(shù)據(jù)綜合成能夠囊括所有指標(biāo)信息的三個主成分，即主成分因子F1、F2、F3.由SPSS25統(tǒng)計軟件輸出的結(jié)果，得到成分得分系數(shù)矩陣，如表5所示，以標(biāo)準(zhǔn)化后的數(shù)據(jù)與成分得分系數(shù)相乘可計算出三個主成分?jǐn)?shù)值結(jié)果，計算公式如下：

表5 成分得分系數(shù)矩陣

F1=0.178zx1+0.148zx2+0.219zx3-0.094zx4-

0.184zx5-+0.018zx6+0.216zx7-0.155zx8+

0.179zx9+0.126zx10-0.031zx11-0.023zx12

(10)

F2=0.025zx1+0.119zx2-0.037zx3-0.108zx4-

0.311zx5+0.257zx6-0.077zx7-0.035zx8-

0.133zx9+0.062zx10+0.333zx11+0.232zx12

(11)

F3=0.046zx1+0.094zx2-0.108zx3+0.52zx4+

0.273zx5-0.291zx6-0.329zx7+0.065zx8-

0.02zx9+0.177zx10-0.017zx11+0.22zx12

(12)

根據(jù)上述公式計算出F1、F2、F3得分后，計算出F1、F2、F3方差貢獻率占總方差貢獻率的比重，依據(jù)方差貢獻率比重進行加權(quán)，計算出各個市的土地生態(tài)安全綜合得分值F.根據(jù)F的得分大小，對安徽省16個地級市生態(tài)安全狀況進行等級劃分和排名，分析生態(tài)狀況原因，從而提出有效建議.F的計算公式如下：

F=(42.209/77.629)*F1+

(24.376/77.629)*F2+(11.044/77.629)*F3

(13)

表6列舉出了2019年安徽省各個地級市土地生態(tài)安全狀況綜合指數(shù)F值和等級排名.

4.5.2 確定評價標(biāo)準(zhǔn)

目前，國內(nèi)關(guān)于土地生態(tài)安全的統(tǒng)一評判標(biāo)準(zhǔn)尚未達到一致，故根據(jù)安徽省土地資源現(xiàn)狀及特點，參考相關(guān)研究成果[15]，將安徽省土地生態(tài)安全分為三個等級，如表7所示.F值大于1的地區(qū)處于Ⅰ等級，為安全狀態(tài)；F值在0～1之間的地區(qū)處于Ⅱ等級，為臨界安全；F值小于0之間的地區(qū)處于Ⅲ等級，為較不安全狀態(tài).具體來看.2019年合肥市的F值為1.854 5，達到安全等級；蕪湖、滁州、蚌埠、阜陽、亳州以及馬鞍山市F值在0～1之間，處于進階安全狀態(tài)；淮北、六安、宣城、安慶、淮南、宿州、銅陵、黃山、池州的F值小于0，土地生態(tài)安全狀況是極不安全狀態(tài).

表6 2019年各個市綜合得分及排序

表7 安徽省土地生態(tài)安全評判標(biāo)準(zhǔn)值

5 聚類分析

為避免主成分分析各個指標(biāo)之間相關(guān)性問題，以主成分F1、F2、F3為為聚類分析數(shù)據(jù)，利用SPSS 25對安徽省16個地級市土地生態(tài)安全進行聚類分析，得到安徽省16個地級市土地生態(tài)安全聚類分析的譜系圖，如圖2所示.

圖2 2019年安徽省土地生態(tài)安全譜系圖

從2019年安徽省生態(tài)安全譜系圖上來看，聚類分析得到的分類結(jié)果與上述主成分分析得到的排序結(jié)果高度吻合.把聚類分析結(jié)果劃分為三類：第一類地區(qū)僅有合肥；第二類地區(qū)有蕪湖、馬鞍山、宣城和滁州；第三類有宿州、阜陽、亳州、淮北、安慶、六安、蚌埠、黃山、銅陵、淮南和池州.合肥市在主成分分析結(jié)果中屬于第一類城市，蕪湖、馬鞍山屬于第二類城市，淮南、銅陵、池州等屬于第三類城市，與聚類分析結(jié)果高度相同.一方面，從實證研究角度反映出本研究所構(gòu)建的主成分分析模型具有較高的合理性和科學(xué)性；另一方面，從安徽省土地生態(tài)環(huán)境的實際情況來看，主成分聚類分析的結(jié)果與近幾年來安徽省土地生態(tài)安全的實際發(fā)展情況高度吻合，且與已經(jīng)發(fā)表的有關(guān)安徽省土地生態(tài)安全評價的研究成果中的結(jié)果高度一致.

第Ⅰ類為合肥市.合肥市是安徽省的省會，全省經(jīng)濟資源、教育資源和科技資源向合肥傾斜，故近幾年合肥經(jīng)濟發(fā)展態(tài)勢迅猛，在經(jīng)濟發(fā)展的同時也在提高生態(tài)文明建設(shè).合肥市的經(jīng)濟密度、人均GDP、城鎮(zhèn)化水平、建成區(qū)面積等指標(biāo)均居于全省首位，經(jīng)濟發(fā)展因子、人口發(fā)展?fàn)顟B(tài)因子數(shù)值較高，對比之下土地結(jié)構(gòu)利用因子較低.雖然合肥市土地生態(tài)健康處于安全狀態(tài)，但是也要重視土地結(jié)構(gòu)的合理運用，優(yōu)化土地結(jié)構(gòu)，經(jīng)濟和生態(tài)共同發(fā)展.

在第Ⅱ類中共有6個地級市，分別為蕪湖、滁州、蚌埠、阜陽、亳州和馬鞍山.限于篇幅，僅以馬鞍山為例進行分析.馬鞍山地處南京都市圈、皖江城市帶承接產(chǎn)業(yè)轉(zhuǎn)移示范區(qū)，還是中國十大鋼鐵基地之一馬鋼所在地，馬鞍山港是長江十大港口之一，皖江第一港口，是中國重要的鋼鐵流通基地.馬鞍山靠近南京，地理位置優(yōu)越，但是經(jīng)濟成分單一.從主成分分析結(jié)果來看，馬鞍山F1(0.798 1)排名全省第三，F(xiàn)2(-0.937 3)和F3(-0.577 4)排名全省倒數(shù)第四，這說明馬鞍山的經(jīng)濟發(fā)展因子水平較高，但土地結(jié)構(gòu)因子和人口發(fā)展?fàn)顟B(tài)因子發(fā)展水平較差，馬鞍山在經(jīng)濟發(fā)展的同時未重視土地結(jié)構(gòu)的合理規(guī)劃和人口發(fā)展對土地生態(tài)的影響.2019年馬鞍山經(jīng)濟運行保持平穩(wěn)，但與合肥、蕪湖相比經(jīng)濟水平還是有一定的差距.在土地利用方面，農(nóng)藥施用強度較高、單位面積播種糧食產(chǎn)量較高，對土地資源的依賴性程度高.在全省范圍內(nèi)，馬鞍山的土地結(jié)構(gòu)利用水平不突出，其土地生態(tài)安全存在一定的風(fēng)險.

第Ⅲ類包括9個城市，分別為淮北、六安、宣城、安慶、淮南、宿州、銅陵、黃山和池州，選取排在最后一位的池州市展開分析.由主成分得分情況看池州市F1(-0.934 1)、F2(-1.387 5)和F3(-0.122 1)都為負(fù)值，且低于全省平均水平.在經(jīng)濟發(fā)展方面，池州的2019年經(jīng)濟密度、地區(qū)生產(chǎn)總值、第三產(chǎn)業(yè)占GDP密度均為全省排名倒數(shù)第二，經(jīng)濟發(fā)展嚴(yán)重依賴旅游業(yè)，經(jīng)濟結(jié)構(gòu)不合理.從土地生態(tài)方面來看，池州2019年農(nóng)藥施用強度為0.027 t/104m2，為全省農(nóng)藥使用強度最大的地級市. 隨著生態(tài)問題越來越突出，對區(qū)域土地生態(tài)安全進行綜合得分計算，能為安徽省土地資源的利用和開發(fā)提供借鑒和參考.由上述主成分聚類分析結(jié)果可知，處于生態(tài)安全狀態(tài)的地級市是安徽省省會城市合肥，蕪湖的生態(tài)安全綜合數(shù)值也接近1，在以EES構(gòu)建的評價指標(biāo)體系中，經(jīng)濟與土地生態(tài)安全水平存在正相關(guān)，經(jīng)濟社會高質(zhì)量發(fā)展可以和生態(tài)安全同向發(fā)展，這給由于經(jīng)濟快速發(fā)展而使環(huán)境遭到破壞的地區(qū)提供發(fā)展參考.根據(jù)上述聚類分析結(jié)果，安徽省應(yīng)格外重視臨界安全和較不安全狀態(tài)下的市級區(qū)域，加大政府財務(wù)撥款，支持科技創(chuàng)新完成土地生態(tài)安全水平的高質(zhì)量發(fā)展，避免土地生態(tài)繼續(xù)遭到破壞.

6 結(jié) 語

根據(jù)評價結(jié)果分析可知，安徽省16個地級市土地生態(tài)安全類別與其經(jīng)濟水平同向發(fā)展，不同地區(qū)應(yīng)采取差異化措施提高自身生態(tài)安全發(fā)展水平.對于Ⅰ類地區(qū)，應(yīng)充分利用優(yōu)越的地理位置和自然環(huán)境資源，發(fā)揮好科技推動經(jīng)濟、生態(tài)高質(zhì)量發(fā)展的優(yōu)勢，在經(jīng)濟社會條件飛速發(fā)展的同時，土地生態(tài)水平不走下坡路.對于Ⅱ類地區(qū)，如蕪湖、馬鞍山，應(yīng)把經(jīng)濟和環(huán)境發(fā)展協(xié)調(diào)起來，重視土地結(jié)構(gòu)和土地承載能力，設(shè)定合理閾值、處理好總量貢獻與發(fā)展速度之間的關(guān)系，穩(wěn)步提高土地生態(tài)安全水平.對于第Ⅲ類地區(qū)，應(yīng)是安徽省提高土地生態(tài)安全質(zhì)量重點關(guān)心的對象，由于經(jīng)濟結(jié)構(gòu)不完整、土地不合理使用造成了其土地生態(tài)的破壞，重視土地生態(tài)安全，根據(jù)不同市域地區(qū)制定不同的投資計劃，把資金真正用到能改善土地生態(tài)安全上來，避免因為經(jīng)濟發(fā)展給土地生態(tài)造成更大的破壞.

本文通過主成分聚類分析，選擇12個多維指標(biāo)對安徽省16個地級市土地生態(tài)安全進行評價，最終將其劃分為三個不同等級，即安全狀態(tài)、臨界安全狀態(tài)和極不安全狀態(tài)，所得結(jié)果與安徽省土地生態(tài)安全實際情況相吻合，反映出本研究具有科學(xué)合理性，可為提高安徽省各個地級市土地生態(tài)安全水平、科學(xué)制定相關(guān)政策提供決策與參考.