基于文獻(xiàn)分析篩選的武漢城市圈土地生態(tài)安全研究

景鈺涵 周天宇 王雄 黃威

摘要：為了促進(jìn)武漢城市圈整體的土地可持續(xù)利用，均衡系統(tǒng)內(nèi)部各市域土地生態(tài)安全水平，探究城市個(gè)體與城市圈整體的狀態(tài)和差異。以各市域和整體區(qū)域?yàn)檠芯繉?duì)象，以改進(jìn)的文獻(xiàn)參考法為基礎(chǔ)，結(jié)合DPSIR模型構(gòu)建武漢城市圈土地生態(tài)安全評(píng)價(jià)指標(biāo)體系，使用主成分分析法計(jì)算2015年各城市土地生態(tài)安全水平綜合值，并比較了整體與各市域水平間的差異。結(jié)果表明，武漢城市圈整體土地生態(tài)安全水平良好，但城市圈內(nèi)部各市域發(fā)展較不均衡；總體表現(xiàn)為南北縱向水平差異優(yōu)于東西橫向差異，并呈現(xiàn)出從中心至外環(huán)逐漸惡化的趨勢(shì)；城市規(guī)模、城市化水平與土地生態(tài)安全水平?jīng)]有直接的線性關(guān)系；根據(jù)地生態(tài)安全值，將組成城市圈的9個(gè)獨(dú)立市域劃分為4個(gè)階梯并給出相應(yīng)的措施進(jìn)行改善。

關(guān)鍵詞：生態(tài)安全；武漢城市圈；評(píng)價(jià)；主成分分析

中圖分類號(hào)：F301.2；F062.1 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：0439-8114（2018）11-0019-06

DOI：10.14088/j.cnki.issn0439-8114.2018.11.004

Abstract： In order to promote the sustainable use of the land in Wuhan Metropolitan Area，to balance the level of land ecological security in all municipalities within the system，explore the status and differences of urban individual and the whole area. Therefore，based on the improved literature reference method and the DPSIR model，the evaluation index system of land ecological security of Wuhan Metropolitan Area was constructed，using principal component analysis to calculate the comprehensive value of urban land ecological security in 2015，then，compared the difference between the whole and the city level. The results show that：the overall ecological security level of Wuhan Metropolitan Area is good，but the development of cities within the group is relatively uneven. There is no direct linear relationship between urban scale，urbanization level and land ecological security level. According to the ecological security value， the nine independent cities that make up the Wuhan Metropolitan Area are divided into four ranks and the corresponding measures are given.

Key words： ecological security； Wuhan Metropolitan area； evaluation； principal component analysis

中國正處在工業(yè)化快速擴(kuò)張和新型城鎮(zhèn)化進(jìn)展階段，各種不合理的用地方式和經(jīng)濟(jì)建設(shè)給土地的生態(tài)安全帶來了嚴(yán)峻的挑戰(zhàn)。基于快速發(fā)展的經(jīng)濟(jì)與用地矛盾的突出、“五位一體”發(fā)展戰(zhàn)略和可持續(xù)發(fā)展的要求，以及城市圈一體化建設(shè)“兩型社會(huì)”進(jìn)程的不斷推進(jìn)，土地生態(tài)安全問題愈來愈受關(guān)注[1，2]。土地資源生態(tài)安全是區(qū)域可持續(xù)發(fā)展的核心，也一直是學(xué)科內(nèi)關(guān)注的熱點(diǎn)。

以Costanza開展的生態(tài)系統(tǒng)健康基本理論和評(píng)估研究為代表，國外土地生態(tài)安全評(píng)價(jià)側(cè)重于土地生態(tài)系統(tǒng)健康、生態(tài)風(fēng)險(xiǎn)、可持續(xù)評(píng)價(jià)[3]，其他學(xué)者也關(guān)注土地生態(tài)退化評(píng)價(jià)[4]和生態(tài)脆弱性評(píng)價(jià)[4-6]。國內(nèi)學(xué)者們主要以土地生態(tài)綜合效益的評(píng)價(jià)、土地生態(tài)安全預(yù)警（監(jiān)測(cè)）和生態(tài)價(jià)值測(cè)算以及補(bǔ)償為研究方向。對(duì)土地生態(tài)安全水平的研究也主要是以評(píng)價(jià)為工具[1-3]，較為成熟的評(píng)價(jià)方法有生態(tài)足跡、物元模糊、神經(jīng)網(wǎng)絡(luò)、TOPSIS和主成分分析法[1-7]，研究對(duì)象也主要以時(shí)間序列的單一市域?yàn)橹黧w，較少有文章通過比較分析城市圈系統(tǒng)與其內(nèi)部組成個(gè)體的差異來研究區(qū)域的土地生態(tài)安全。然而，與單個(gè)城市相比，城市圈的土地利用是一個(gè)更為開放的復(fù)雜系統(tǒng)，其時(shí)空演變過程不僅受到區(qū)域內(nèi)外自然、社會(huì)、經(jīng)濟(jì)和其他各種不確定因素的綜合影響，表現(xiàn)出顯著的復(fù)雜非線性特征[7]，內(nèi)部差異也會(huì)對(duì)系統(tǒng)整體水平產(chǎn)生直接的影響。以武漢市為首的“8+1”城市圈正處于社會(huì)經(jīng)濟(jì)發(fā)展的高速增長階段，社會(huì)生產(chǎn)總值和人均收入增長迅速，對(duì)各類用地?cái)?shù)量和質(zhì)量的需求都將日益增大。因此，結(jié)合系統(tǒng)內(nèi)外部，差異性比較探究城市圈土地生態(tài)安全水平是促進(jìn)武漢城市圈整體發(fā)展的必要途徑和根本手段，同時(shí)也有助于進(jìn)一步了解各城市生態(tài)安全水平與差異，促進(jìn)各城市的土地生態(tài)良性發(fā)展，使生態(tài)安全地區(qū)帶動(dòng)生態(tài)脆弱地區(qū)，引導(dǎo)區(qū)域利用優(yōu)勢(shì)、規(guī)避劣勢(shì)、提高綜合水平并實(shí)現(xiàn)土地的可持續(xù)發(fā)展。

土地生態(tài)安全是多因素影響與相互作用的復(fù)雜動(dòng)態(tài)系統(tǒng)，合理的評(píng)價(jià)指標(biāo)體系是科學(xué)準(zhǔn)確地評(píng)價(jià)土地生態(tài)安全水平的關(guān)鍵[7]。現(xiàn)有的研究主要以自然-經(jīng)濟(jì)-社會(huì)、PSR、DPSIR為模型，為了促進(jìn)武漢城市圈整體的土地可持續(xù)利用，均衡系統(tǒng)內(nèi)部各市域土地生態(tài)安全水平，探究城市個(gè)體與城市圈整體的狀態(tài)和差異。以各市域和整體區(qū)域?yàn)檠芯繉?duì)象，以改進(jìn)的文獻(xiàn)參考法為基礎(chǔ)，結(jié)合DPSIR模型構(gòu)建武漢城市圈土地生態(tài)安全評(píng)價(jià)指標(biāo)體系，使用主成分分析法計(jì)算2015年各城市土地生態(tài)安全水平綜合值，并比較了整體與各市域水平間的差異，以期為這方面的研究提供參考。

1 研究區(qū)概況

武漢城市圈又稱“8+1”城市圈，是指以湖北省省會(huì)武漢市為圓心，包括咸寧、天門、黃石、潛江、仙桃、黃岡、孝感、鄂州8個(gè)城市的城市圈[8]。

武漢城市圈地處江漢平原地區(qū)、氣候溫暖濕潤，水熱條件優(yōu)越，地貌類型多樣，以平原為主，丘陵、山地為輔，大致呈現(xiàn)出“兩山、兩崗、三丘、三原”的格局；平均海拔不足100 m，地勢(shì)平坦、河湖眾多；區(qū)域土地利用形態(tài)由內(nèi)向外呈現(xiàn)出多層次環(huán)狀分布，其中耕地主要分布于中部和東部地區(qū)；林地主要分布于東北部的黃岡和南部的咸寧，較少的草地面積主要集中于黃岡，建設(shè)用地主要集中在武漢市（圖1）。面積不到全省1/3的武漢城市圈，集中了湖北省一半的人口、60%以上的GDP總量，不僅是湖北經(jīng)濟(jì)發(fā)展的核心區(qū)域，也是中部崛起的重要戰(zhàn)略支點(diǎn)[8，9]。人口、經(jīng)濟(jì)與用地現(xiàn)狀決定了武漢城市圈土地生態(tài)安全水平不容樂觀，研究城市圈與各組成地域差異并分析原因以實(shí)現(xiàn)區(qū)域均衡發(fā)展刻不容緩。

2 方法與數(shù)據(jù)

2.1 研究方法

2.1.1 主成分分析法原理 主成分分析方法是將多個(gè)具有相關(guān)性的要素轉(zhuǎn)化成幾個(gè)不相關(guān)的綜合指標(biāo)的分析與統(tǒng)計(jì)方法。綜合指標(biāo)有可能包含眾多相互重復(fù)的信息，主成分分析在保證信息最少丟失原則下，對(duì)原來指標(biāo)進(jìn)行降維處理，把一些不相關(guān)的指標(biāo)省去，將原來較多的指標(biāo)轉(zhuǎn)換成能反映研究現(xiàn)象的較少的綜合指標(biāo)，這樣能夠簡(jiǎn)化復(fù)雜的研究，在保證研究精確度的前提下提高研究效率[10]。本研究采用SPSS17.0軟件對(duì)所選指標(biāo)進(jìn)行主成分分析，以達(dá)到在多指標(biāo)中篩選重要的綜合指標(biāo)進(jìn)行生態(tài)評(píng)價(jià)。

2.1.2 改進(jìn)的文獻(xiàn)參考法與指標(biāo)體系的構(gòu)建 在中國知網(wǎng)上分別以“土地生態(tài)”“生態(tài)安全水平”“生態(tài)評(píng)價(jià)”為檢索詞條，分別按照“被引量”“下載次數(shù)”和“主題排序”為依據(jù)，各搜集15～20篇相關(guān)文獻(xiàn)。將搜集的文獻(xiàn)篩選、去除重復(fù)并摘錄出各文獻(xiàn)中評(píng)價(jià)指標(biāo)體系中的指標(biāo)層，統(tǒng)計(jì)不同名稱指標(biāo)出現(xiàn)的頻率，同時(shí)將具有相同含義的指標(biāo)合并為同一指標(biāo)、合并頻次（如“單位耕地化肥施用量”和“全年化肥用量”“城市人口與用地增長彈性系數(shù)”和“人口承載力”以及“人口密度”）。最終從相關(guān)文獻(xiàn)[11-22]中摘錄出31個(gè)有效的評(píng)價(jià)指標(biāo)體系，合并315個(gè)指標(biāo)為150個(gè)，并從中篩選出頻次大于5的38個(gè)指標(biāo)。

DPSIR模型有較強(qiáng)的系統(tǒng)性，以該模型為基礎(chǔ)，考慮武漢城市圈的實(shí)際情況根據(jù)客觀性、可操作性和實(shí)用性原則，刪減無效指標(biāo)，構(gòu)建含26個(gè)指標(biāo)的區(qū)域土地生態(tài)安全評(píng)價(jià)指標(biāo)體系（表1）。

2.2 數(shù)據(jù)處理

2.2.1 數(shù)據(jù)來源 數(shù)據(jù)基礎(chǔ)來源于2016年武漢、咸寧、天門、黃石、潛江、仙桃、黃岡、孝感、鄂州統(tǒng)計(jì)年鑒與國民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)、2015年鄂州統(tǒng)計(jì)年鑒、2016年武漢市綠化狀況公報(bào)全省第四次森林資源二類調(diào)查、2016年咸寧市政府工作報(bào)告（咸寧市人民政府/信息公開/工作報(bào)告）、湖北省人民政府網(wǎng)站、潛江熱線網(wǎng)、2016年武漢市綠化狀況公報(bào)、中共天門市委辦公室網(wǎng)站和鄂州市環(huán)境保護(hù)局2015年收支預(yù)算總表（鄂州市環(huán)境保護(hù)局網(wǎng)站）。由于缺乏土地利用結(jié)構(gòu)多樣性數(shù)據(jù)、土壤退化與自然災(zāi)害數(shù)據(jù)，研究結(jié)果與現(xiàn)實(shí)情況可能存在微小的偏差。

2.2.2 數(shù)據(jù)預(yù)處理 要得到標(biāo)準(zhǔn)化評(píng)價(jià)矩陣，可以采用歸一化方法對(duì)原始數(shù)據(jù)進(jìn)行處理，對(duì)于正向指標(biāo)，其公式如下。

對(duì)于負(fù)向指標(biāo)，其公式為：

標(biāo)準(zhǔn)化得到關(guān)于生態(tài)效益評(píng)價(jià)體系矩陣R：

其中，xij為第i個(gè)城市第j個(gè)指標(biāo)的原始數(shù)值；其中，rij為第i個(gè)城市第j個(gè)指標(biāo)的標(biāo)準(zhǔn)化值；i=1，2...，m，m為城市個(gè)數(shù)；j=1，2...，n，n為評(píng)價(jià)指標(biāo)數(shù)。

2.2.3 主成分分值計(jì)算 通過文獻(xiàn)篩選的指標(biāo)數(shù)量較多，存在相關(guān)性較大的問題，主成分分析法在一定程度上能夠減少指標(biāo)數(shù)量，因此選擇主成分分析進(jìn)行綜合評(píng)價(jià)。將已經(jīng)標(biāo)準(zhǔn)化處理的數(shù)據(jù)在SPSS17.0軟件中進(jìn)行處理，通過分析指標(biāo)間關(guān)系，提取了7個(gè)主成分。這7個(gè)成分的累計(jì)值已經(jīng)超過97.359%（表2），因此通過這7項(xiàng)可以計(jì)算出各因子對(duì)原始指標(biāo)的載荷狀況。

根據(jù)成分得分系數(shù)矩陣和主成分特征值計(jì)算因子對(duì)原始指標(biāo)載荷狀況（表3），該矩陣反應(yīng)的主成分系數(shù)表明：第一主成分主要反映了人均GDP、經(jīng)濟(jì)密度、固定資產(chǎn)投入、第三產(chǎn)業(yè)產(chǎn)值以及農(nóng)民人均可支配收入信息；第二主成分主要反映了環(huán)保投資比重、耕地面積比重、森林覆蓋以及城市人均公共綠地面積信息；第三主成分主要反映了排污、機(jī)耕信息；第四主成分主要反映了糧食單產(chǎn)、固體廢棄物的利用和農(nóng)藥使用情況信息；第五主成分主要反映化肥使用情況和建成區(qū)綠化覆蓋率信息；第六主成分主要反映單位GDP能耗、污水排放、環(huán)保投資和建成區(qū)綠地率信息；第七主成分主要反映了城鎮(zhèn)生活污水處理率、人均公共綠地等信息，其中污水處理率占據(jù)絕對(duì)影響地位。

綜合來說，主要貢獻(xiàn)指標(biāo)因子為單位GDP污水排放、建成區(qū)綠化覆蓋、環(huán)保投資占比、農(nóng)民人均可支配收入、農(nóng)業(yè)機(jī)械總動(dòng)力和旱澇保收率，進(jìn)而得到7個(gè)主成分和各指標(biāo)因子之間的關(guān)系模型：

其中，Zpi為第p個(gè)主成分第i個(gè)城市的主成分得分值，Cpj為第p個(gè)主成分第j個(gè)指標(biāo)的因子對(duì)指標(biāo)載荷系數(shù)，Xij為第i個(gè)城市第j個(gè)指標(biāo)的標(biāo)準(zhǔn)化數(shù)值（表3）。

2.2.4 生態(tài)安全評(píng)價(jià)綜合得分與排序 根據(jù)各指標(biāo)特征值的特征向量、標(biāo)準(zhǔn)化數(shù)據(jù)，分別以7個(gè)主成分的特征值為依據(jù)權(quán)重計(jì)算主成分關(guān)系系數(shù)，并構(gòu)建主成分綜合評(píng)價(jià)模型：

其中，F(xiàn)i為第i個(gè)城市的土地生態(tài)安全水平綜合得分值，Wp為第p個(gè)主成分的權(quán)重值，Zpi見式（3）。

計(jì)算各研究區(qū)域在7個(gè)主成分下各子系統(tǒng)安全分值和綜合安全分值（表4）。分值為正則表明該主成分在平均發(fā)展水平1之上，且得分越高說明生態(tài)安全越好。

3 評(píng)價(jià)結(jié)果分析

3.1 土地生態(tài)安全子系統(tǒng)與整體系統(tǒng)水平分析

以9個(gè)城市個(gè)體和整體的總體生態(tài)水平變化參考排序，從武漢市到咸寧市、城市圈整體、仙桃市、黃石市、孝感市、潛江市、黃岡市、天門市和鄂州市土地生態(tài)安全水平逐漸下降，即安全穩(wěn)定性由強(qiáng)到弱。為了更清楚地看出子系統(tǒng)安全性分值與總系統(tǒng)得分值的差異，結(jié)合圖2進(jìn)行比較說明。

從城市來看，城市圈整體土地生態(tài)水平較高，僅低于武漢市和咸寧市，仙桃市、黃石市和孝感市緊追其后位于第二生態(tài)安全梯隊(duì)，潛江市和黃岡市位于土地生態(tài)安全水平第三梯隊(duì)，天門市和鄂州市位于最后，尤其是鄂州市的土地生態(tài)安全水平最不容樂觀。武漢市的土地生態(tài)安全性最好，雖然市區(qū)的人口密度、城鎮(zhèn)化率相對(duì)較高，人類活動(dòng)更加明顯，但是市區(qū)第三產(chǎn)業(yè)占據(jù)絕對(duì)地位，在很大程度上避免了以農(nóng)業(yè)為主的第一二產(chǎn)業(yè)給生態(tài)環(huán)境帶來的種種破壞。武漢市的驅(qū)動(dòng)力子系統(tǒng)為負(fù)值，低于平均水平，分析各主成分分值發(fā)現(xiàn)其第三、四、五、七主成分分值都為負(fù)值且第七主成分占絕對(duì)比重，說明以城鎮(zhèn)生活污水處理率、城市公共人均綠地面積和單位產(chǎn)值GDP排放信息為反映的指標(biāo)水平較差。

從系統(tǒng)來看，隨著土地生態(tài)總系統(tǒng)水平的變化趨勢(shì)，子系統(tǒng)變化呈現(xiàn)平穩(wěn)降低的趨勢(shì)。其中狀態(tài)與響應(yīng)子系統(tǒng)在區(qū)域?qū)哟紊嫌胁▌?dòng)；驅(qū)動(dòng)力系統(tǒng)低于平均水平，影響驅(qū)動(dòng)力系統(tǒng)的指標(biāo)主要是森林覆蓋、人口與城市及經(jīng)濟(jì)水平，這說明對(duì)獨(dú)立的城市來說都存在城市發(fā)展與經(jīng)濟(jì)水平和綠化覆蓋等不協(xié)調(diào)的狀態(tài)影響了土地生態(tài)安全性。咸寧市土地生態(tài)安全值在驅(qū)動(dòng)力子系統(tǒng)中處于最高水平，這是由于咸寧市較高的森林覆蓋率水平、較少的化肥投入、較高的單位GDP中，環(huán)保投資額和合理的耕地面積占據(jù)絕對(duì)優(yōu)勢(shì)。咸寧市是湖北省首批低碳經(jīng)濟(jì)試點(diǎn)市，首批全國旅游標(biāo)準(zhǔn)化城市并榮獲“國家森林城市”的榮譽(yù)稱號(hào)，這一信息也與概述中描述的“林業(yè)集中地主要位于咸寧市”相吻合。

3.2 區(qū)域土地生態(tài)安全均衡發(fā)展分析

為了均衡發(fā)展城市圈，以城市圈整體水平為基準(zhǔn)數(shù)據(jù)，與其他9個(gè)城市的子系統(tǒng)的分值分別作差得到每個(gè)城市與整體水平的差異值，并通過ArcGIS10.2等間距分級(jí)使其可視化（圖3），再進(jìn)行城市與城市、城市個(gè)體與城市圈整體之間的比較。

1）通過區(qū)域差異分析可以發(fā)現(xiàn)，就安全性而言，無論按5個(gè)獨(dú)立子系統(tǒng)還是總體系統(tǒng)安全分值來看，各城市與城市圈整體土地生態(tài)安全值差值皆為負(fù)值，且城市圈整體生態(tài)安全水平較高，說明居于城市圈整體水平之上的城市和之下的城市水平差異顯著，優(yōu)劣兩極分化，城市圈內(nèi)各市域土地安全性并不均衡。

2）天門市和鄂州市的壓力系統(tǒng)差異較為失衡，這也是天門市和鄂州市總體土地生態(tài)安全水平位列第四梯隊(duì)的直接原因。壓力系統(tǒng)反映了人口自然增長、森林覆蓋、城鎮(zhèn)化率和人均生產(chǎn)總值指標(biāo)信息，結(jié)合主成分反應(yīng)各指標(biāo)信息情況得出。為了縮小區(qū)域差距，天門市和鄂州市的人口增長和經(jīng)濟(jì)發(fā)展的突出矛盾是導(dǎo)致土地生態(tài)系統(tǒng)脆弱的主導(dǎo)因素。

3）黃岡市的響應(yīng)系統(tǒng)差異失衡，響應(yīng)系統(tǒng)反映固定資產(chǎn)、環(huán)保投資、經(jīng)濟(jì)密度、單位能耗和廢棄物處理利用情況，表明黃岡市的經(jīng)濟(jì)和環(huán)保治理狀況使得土地生態(tài)安全水平落后于其他城市，應(yīng)該在這方面加以改進(jìn)。

4）武漢市和咸寧市總體水平良好，但因經(jīng)濟(jì)發(fā)展要求和城市發(fā)展定位，這2個(gè)城市在耕地面積、治理方面無法改進(jìn)，因此這兩個(gè)城市保持平穩(wěn)現(xiàn)狀即可。

5）武漢城市圈土地生態(tài)安全總體呈現(xiàn)中心至外環(huán)逐漸惡化的趨勢(shì)。總體來講，區(qū)域南北差異小于東西差異，東北部水平偏低。此外，安全水平較低的城市總體表現(xiàn)為區(qū)域小面積，因此推斷城市規(guī)模也并不一定意味著土地生態(tài)的惡化，或兩者之間沒有之間直接的線性關(guān)系。

4 小結(jié)

通過分析武漢城市圈與個(gè)體、總系統(tǒng)與子系統(tǒng)的土地生態(tài)安全綜合值表明：

1）武漢城市圈整體土地生態(tài)安全水平良好，但城市圈內(nèi)部各市域?qū)哟伟l(fā)展較不均衡；總體表現(xiàn)為南北縱向水平差異優(yōu)于東西橫向差異，并呈現(xiàn)出從中心至外環(huán)逐漸惡化的趨勢(shì)。

2）根據(jù)分析可以將城市圈內(nèi)部劃分為4個(gè)區(qū)域，即：第一梯隊(duì)為武漢市和咸寧市，不需要改善；第二梯隊(duì)為仙桃市、黃石市和孝感市，可以改善；第三梯隊(duì)為潛江市和黃岡市，需要改善；第四梯隊(duì)為天門市和鄂州市，必須改善；并提出了相應(yīng)的改善方向。

3）城市化水平或者建設(shè)用地聚集度與土地生態(tài)安全水平?jīng)]有直接關(guān)系。除了土地的農(nóng)業(yè)和生態(tài)價(jià)值，對(duì)于城市化水平比較高的城市，第三產(chǎn)業(yè)所創(chuàng)造的生態(tài)服務(wù)價(jià)值也可以規(guī)避人口密度等對(duì)土地帶來的生態(tài)破壞。因此，合理的產(chǎn)業(yè)結(jié)構(gòu)可以改善土地生態(tài)安全狀況。此外，土地生態(tài)安全不僅體現(xiàn)在其自然功能上，經(jīng)濟(jì)功能更是重要的因素之一。

4）城市規(guī)模擴(kuò)張并不代表土地生態(tài)系統(tǒng)的惡化，或者說兩者之間沒有直接的線性關(guān)系。

參考文獻(xiàn)：

[1] 余 健，房 莉，倉定幫，等.熵權(quán)模糊物元模型在土地生態(tài)安全評(píng)價(jià)中的應(yīng)用[J].農(nóng)業(yè)工程學(xué)報(bào)，2012（5）：260-266.

[2] 李德勝，王占岐，藍(lán) 希.城市土地生態(tài)安全評(píng)價(jià)及障礙因子研究——以武漢市為例[J].中國國土資源經(jīng)濟(jì)，2017，30（8）：40-44，73.

[3] 李春燕，南 靈.陜西省土地生態(tài)安全動(dòng)態(tài)評(píng)價(jià)及障礙因子診斷[J].中國土地科學(xué)，2015，29（4）：72-81.

[4] 戴 靚，姚新春，周生路，等.長三角經(jīng)濟(jì)發(fā)達(dá)區(qū)金壇市土地生態(tài)狀況評(píng)價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào)，2013，29（8）：249-257.

[5] WANG S Y，LIU J S，YANG C J.Eco-environmental vulnerability evaluation in the Yellow River Basin，China[J].Pedosphere，2008，18（2）：171-182.

[6] TRAN L T，KNIGHT C G，O'NEILL R V，et al. Fuzzy decision analysis for integrated environmental vulnerability assessment of the mid-atlantic region[J].Environmental Assessment，2002，29（6）：845-859.

[7] 徐 美，朱 翔，李靜芝.基于DPSIR-TOPSIS模型的湖南省土地生態(tài)安全評(píng)價(jià)[J].冰川凍土，2012，34（5）：1265-1272.

[8] 雷 征，董 捷.武漢城市圈土地利用結(jié)構(gòu)演變規(guī)律及驅(qū)動(dòng)因素分析[J].農(nóng)業(yè)現(xiàn)代化研究，2010，31（2）：147-151.

[9] 向蕓蕓，蒙吉軍.基于生態(tài)效益的武漢城市圈土地利用結(jié)構(gòu)優(yōu)化[J].長江流域資源與環(huán)境，2013，22（10）：1297-1304.

[10] 王 鵬，況福民，鄧育武，等.基于主成分分析的衡陽市土地生態(tài)安全評(píng)價(jià)[J].經(jīng)濟(jì)地理，2015，35（1）：168-172.

[11] 郭宇倫，師學(xué)義，璩路路.基于PSR-CPM模型的市域土地生態(tài)安全評(píng)價(jià)[J].水土保持研究，2017，24（4）：108-112.

[12] 李進(jìn)濤，李貽學(xué)，鞏玉玲，等.基于IDRISI的山東省土地生態(tài)安全評(píng)價(jià)研究[J].中國人口·資源與環(huán)境，2014，24（S3）：190-193.

[13] 孫奇奇，宋 戈，齊美玲.基于主成分分析的哈爾濱市土地生態(tài)安全評(píng)價(jià)[J].水土保持研究，2012，19（1）：234-238.

[14] 黃輝玲，羅文斌，吳次芳.基于物元分析的土地生態(tài)安全評(píng)價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào)，2010，26（3）：316-322.

[15] 吳澤斌，劉衛(wèi)東，汪友結(jié)，等.城市土地生態(tài)利用水平測(cè)度及障礙因子診斷[J].中國人口·資源與環(huán)境，2010，20（2）：24-29.

[16] 李玉平，蔡運(yùn)龍.河北省土地生態(tài)安全評(píng)價(jià)[J].北京大學(xué)學(xué)報(bào)（自然科學(xué)版），2007（6）：784-789.

[17] 鮑 艷，胡振琪，柏 玉，等.主成分聚類分析在土地利用生態(tài)安全評(píng)價(jià)中的應(yīng)用[J].農(nóng)業(yè)工程學(xué)報(bào)，2006（8）：87-90.

[18] 張虹波，劉黎明，張軍連，等.區(qū)域土地資源生態(tài)安全評(píng)價(jià)的物元模型構(gòu)建及應(yīng)用[J].浙江大學(xué)學(xué)報(bào)（農(nóng)業(yè)與生命科學(xué)版），2007（2）：222-229.

[19] 李明月，賴笑娟.基于BP神經(jīng)網(wǎng)絡(luò)方法的城市土地生態(tài)安全評(píng)價(jià)——以廣州市為例[J].經(jīng)濟(jì)地理，2011，31（2）：289-293.

[20] 劉 慶，陳利根，舒幫榮，等.長株潭城市群土地生態(tài)安全動(dòng)態(tài)評(píng)價(jià)研究[J].長江流域資源與環(huán)境，2010，19（10）：1192-1197.

[21] 杜忠潮，韓申山.基于主成分分析的土地生態(tài)安全評(píng)價(jià)實(shí)證研究——以陜西省10個(gè)省轄市為例[J].水土保持通報(bào)，2009，29（6）：198-202，207.

[22] 高桂芹，韓 美.區(qū)域土地資源生態(tài)安全評(píng)價(jià)——以山東省棗莊市中區(qū)為例[J].水土保持研究，2005（5）：275-277.