侯小宇,何 剛,王曉東,谷雅嫻
(安徽理工大學(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ū)域、評價指標體系、土地生態(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)安全評價指標上來看,張鳳太[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)安全評價指標體系,利用SPSS25對安徽省各個地級市進行主成分分析和聚類分析,最后依據(jù)綜合得分,對安徽省16個地級市土地生態(tài)安全進行排名和等級劃分.
本文以安徽省為研究區(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)安全提出了新的要求.
數(shù)據(jù)來源于《2020年安徽省統(tǒng)計年鑒》《2020年國民經(jīng)濟與社會發(fā)展統(tǒng)計公報》.
2.2.1 主成分分析法
主成分分析的數(shù)學(xué)模型為:
(1)
由式(1)可得
(2)
(3)
(4)
(5)
(6)
4)var依次遞減:
(7)
其中:n為樣本數(shù);p為樣本的觀測指標數(shù);X為隨機變量;μ為隨機變量均值;δ為隨機變量的協(xié)方差矩陣;y1,y2,…,yp為生成的新的綜合指標即主成分.
本文利用SPSS 25對標準化后的指標進行主成分分析,實現(xiàn)有效降維、消除指標間線性相關(guān)性,得到相對獨立、重要的、少量的綜合指標.利用各個主成分得分,計算安徽省16個地級市的綜合指標得分,對其進行土地生態(tài)安全等級劃分和排名.
2.2.2 聚類分析
為了消除主成分分析中指標之間的差異性、更大空間尺度分析指標的不同,本文采用聚類分析方法中的系統(tǒng)聚類法進行進一步的驗證[11].系統(tǒng)聚類分析(Hierarchical Cluster Method)包含以下步驟:
1)構(gòu)造n個類,其中只包含一個樣品.
2)計算每個類之間的距離,本文采用中間距離公式,得到距離矩陣,類之間的距離公式如下
(8)
3)將距離最近的兩類合并成一個新類.
4)計算合并后的新類與其他類之間的距離.
5)畫出聚類圖,并根據(jù)聚類圖決定類的個數(shù).
本文以主成分分析的輸出作為聚類分析的輸入,進一步驗證對安徽省16個地級市土地生態(tài)安全等級劃分及排名的準確性.
研究土地生態(tài)安全關(guān)鍵就是評價指標的選擇,選取的指標要有合理性、科學(xué)性,要能覆蓋各個地區(qū)的土地生態(tài)特征,能充分反映各個地區(qū)土地生態(tài)健康狀況.故從指標體系構(gòu)建的科學(xué)性、系統(tǒng)性原則出發(fā),參照相關(guān)研究[12],立足于安徽省土地生態(tài)安全現(xiàn)狀,以“經(jīng)濟-環(huán)境-社會”(EES)模型為基礎(chǔ)構(gòu)建土地生態(tài)安全評價指標體系.其中經(jīng)濟因素有四個正向指標,即經(jīng)濟密度(萬元/km2)、地區(qū)生產(chǎn)總值(億元)、人均GDP(元/人)、第三產(chǎn)業(yè)占GDP比重(%),它們反映了研究區(qū)經(jīng)濟發(fā)展水平和速度;環(huán)境因素有兩個正向指標、兩個負向指標,即農(nóng)藥施用強度(t/hm2)、人均道路面積(m2/人)、單位面積播種糧食產(chǎn)量(t/hm2)、工業(yè)固體廢物綜合利用率(%),它們反映出對土地的利用程度,與土地生態(tài)健康有著直接關(guān)系;社會因素有兩個正向指標、兩個負向指標,即城鎮(zhèn)化水平(%)、建成區(qū)面積(km2)、農(nóng)業(yè)機械總動力(104kW)、人口自然增長率(%),它們通過城鎮(zhèn)化、建成區(qū)、機械動力和人口的增加反映出社會發(fā)展對土地的影響.共選取了12個多維指標,如表1所示.
表1 安徽省土地生態(tài)評價指標體系
為了消除不同指標在單位、數(shù)量級上的差距,需對指標進行標準化處理[13].本文考慮到指標體系的特性,利用SPSS25統(tǒng)計軟件下的描述統(tǒng)計,對安徽省16個地級市12項指標進行z-score標準化.公式如下:
(9)
其中:Z為標準化后的輸出結(jié)果;C為具體指標數(shù)值;μ為樣本平均數(shù);σ為樣本標準差.標準化輸出的結(jié)果即Z值代表著原始指標數(shù)值與總數(shù)值平均數(shù)之間的距離,固有正負之分.
在原始數(shù)據(jù)Z-score標準化之后,需要對指標之間的相關(guān)性進行檢驗.指標之間的相關(guān)性決定著主成分分析提取的主成分是否有效,故需進行因子的適宜性檢驗[14].KMO與巴特利特球形檢驗用于比較變量指標間的簡單相關(guān)系數(shù)和偏相關(guān)系數(shù),其中KMO數(shù)值越接近1,代表著指標變量之間的相關(guān)性越強,越適合做因子分析;巴特利特球形檢驗值是常用的KMO度量標準,0.9以上表示非常合適,0.5以下表示極不合適.SPSS25輸出的結(jié)果如表2,KMO檢驗統(tǒng)計量數(shù)值為0.648,在相關(guān)系數(shù)矩陣中有64.8%的系數(shù)個數(shù)大于0.3,這表明指標體系之間具有很強的相關(guān)性,適合做因子分析;巴特利特球形度檢驗顯著性為0.000,說明原始數(shù)據(jù)服從正態(tài)分布總體,具有進一步研究的價值.
表2 KMO與巴特利特檢驗
指標特征值與貢獻率結(jié)果如圖1、表3所示.結(jié)果顯示:到第四個主成分時,曲線開始趨于平穩(wěn),特征值達到1,故在所有指標中提取出3個主成分因子最為合適.前三個主成分因子的初始方差貢獻率分別為42.209%、24.376%、11.044%,累積方差貢獻率達到77.629%,即提取出來的3個主成分因子能夠涵蓋絕大部分的具有相關(guān)性的指標體系信息,避免了指標之間具有線性相關(guān)關(guān)系或者指標與研究結(jié)果相關(guān)程度不高,從而導(dǎo)致研究結(jié)果不準確等問題.
圖1 碎石圖
表3 各個主成分的總方差解釋和方差貢獻
因子載荷矩陣是提取的主成分命名的主要依據(jù),由于本研究指標未進行因子旋轉(zhuǎn)的因子載荷矩陣數(shù)值較小,無法準確、直觀的看出提取三個主成分代表的含義,故以四次方最大旋轉(zhuǎn)法(Quartimax)進行因子旋轉(zhuǎn),得出因子的載荷矩陣,如表4.
表4 因子載荷矩陣
F1、F2、F3分別表示從12個相關(guān)指標中提取的第一、第二、第三個主成分.由表4可知F1在經(jīng)濟密度、地區(qū)生產(chǎn)總值、人均GDP上的載荷量比較高,這些指標代表著地區(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ā)展狀態(tài)因子.該評價指標體系綜合12個指標全面的反應(yīng)了F1、F2、F3以及土地生態(tài)安全綜合指數(shù)F三個方面的信息,能夠有效地對安徽省土地生態(tài)安全狀況作出評價.
4.5.1 土地生態(tài)安全狀況綜合得分
在指標信息量保持完整的前提下,將安徽省16個地級市12個不同維度的指標數(shù)據(jù)綜合成能夠囊括所有指標信息的三個主成分,即主成分因子F1、F2、F3.由SPSS25統(tǒng)計軟件輸出的結(jié)果,得到成分得分系數(shù)矩陣,如表5所示,以標準化后的數(shù)據(jù)與成分得分系數(shù)相乘可計算出三個主成分數(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 確定評價標準
目前,國內(nèi)關(guān)于土地生態(tài)安全的統(tǒng)一評判標準尚未達到一致,故根據(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)安全評判標準值
為避免主成分分析各個指標之間相關(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ū)面積等指標均居于全省首位,經(jīng)濟發(fā)展因子、人口發(fā)展狀態(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ā)展狀態(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)都為負值,且低于全省平均水平.在經(jīng)濟發(fā)展方面,池州的2019年經(jīng)濟密度、地區(qū)生產(chǎn)總值、第三產(chǎn)業(yè)占GDP密度均為全省排名倒數(shù)第二,經(jīng)濟發(fā)展嚴重依賴旅游業(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)建的評價指標體系中,經(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ù)遭到破壞.
根據(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個多維指標對安徽省16個地級市土地生態(tài)安全進行評價,最終將其劃分為三個不同等級,即安全狀態(tài)、臨界安全狀態(tài)和極不安全狀態(tài),所得結(jié)果與安徽省土地生態(tài)安全實際情況相吻合,反映出本研究具有科學(xué)合理性,可為提高安徽省各個地級市土地生態(tài)安全水平、科學(xué)制定相關(guān)政策提供決策與參考.