基于改進SPA法的耕地占補平衡生態(tài)安全評價

施開放，刁承泰，* ，孫秀鋒，2，左太安，3，蔡 朕，禹陽春

(1.西南大學地理科學學院，重慶 400715;2.西南大學園藝園林學院，重慶 400715;3.畢節(jié)學院環(huán)境與生命科學系，畢節(jié) 551700)

耕地作為土地的精華對我國農(nóng)業(yè)和國民經(jīng)濟發(fā)展起著不可替代的基礎(chǔ)性作用［1-2］。隨著我國城市化和工業(yè)化的快速提升，耕地數(shù)量減少趨勢不可避免，耕地保護問題越發(fā)嚴峻。為解決我國的糧食安全和日益短缺的耕地資源問題，國家實行了嚴厲的耕地保護制度，其中耕地占補平衡就是其中之一［3］。所謂耕地占補平衡是指一定時間、區(qū)域范圍內(nèi)，耕地在局部、微觀不斷被占用、補充及其變化的情況下，在對其數(shù)量、質(zhì)量和開發(fā)利用的可持續(xù)調(diào)控基礎(chǔ)上，為滿足區(qū)域?qū)Z食等基本農(nóng)產(chǎn)品的需求，而做到的耕地生產(chǎn)力供求平衡［1］。當然，耕地占補平衡不僅僅強調(diào)耕地數(shù)量和質(zhì)量的平衡，耕地占補平衡的生態(tài)安全維持也屬其中的重要內(nèi)容［4］。所謂耕地占補平衡生態(tài)安全是指在一定時間、區(qū)域范圍內(nèi)，耕地在微觀、局部不斷被占用、補充及其變化情況下，土地生態(tài)系統(tǒng)能夠保持其功能與結(jié)構(gòu)不受或少受威脅的狀態(tài)，同時，土地生態(tài)系統(tǒng)為人類提供服務(wù)的質(zhì)量和數(shù)量能夠持續(xù)滿足人類生存和發(fā)展的需求［5］，從而達到土地環(huán)境、社會和經(jīng)濟復合體的長期協(xié)調(diào)發(fā)展。就目前的文獻來看，我國對生態(tài)安全的研究主要包括區(qū)域生態(tài)安全評價［6］，生態(tài)風險性評價［7］，土地生態(tài)安全評價［8］，但是對耕地占補平衡生態(tài)安全評價研究較少。因此，對耕地占補平衡進行生態(tài)安全評價，準確掌握耕地占補平衡生態(tài)安全狀態(tài)成為區(qū)域可持續(xù)發(fā)展測度的一個重要內(nèi)容，也是進行土地利用規(guī)劃與建設(shè)和土地生態(tài)安全預(yù)測和預(yù)警的重要依據(jù)［9］。

耕地占補平衡的生態(tài)安全評價實際上是對其占補平衡生態(tài)安全復合系統(tǒng)的全面診斷，包括耕地資源的占補、社會經(jīng)濟的支持以及生態(tài)環(huán)境的維持等多方面的內(nèi)容。目前，綜合評價的方法很多，如模糊數(shù)學法、層次分析法、投影尋蹤法、神經(jīng)網(wǎng)絡(luò)法［10-11］等，每種方法都有自身的優(yōu)點和不足，眾多學者針對各種方法的不足作了相應(yīng)的改進，取得了一定成果。但是由于耕地占補平衡的生態(tài)安全評價影響因素的不確定性，評價指標與生態(tài)安全等級之間存在復雜的非線性關(guān)系，所以至今沒有一個統(tǒng)一的評價模型來進行耕地占補平衡的生態(tài)安全評價。因此，為了克服現(xiàn)有評價方法的不足，本文以永川區(qū)為例，運用改進SPA法對其進行耕地占補平衡生態(tài)安全評價。

1 改進集對分析

1.1 改進SPA法評價原理及模型

集對分析是趙克勤在1989年提出的一種新的系統(tǒng)分析理論方法［12-13］，但該評價方法存在一點不足，是b、c的細化問題，即同異反評語細化問題［14］。原創(chuàng)聯(lián)系度(u(A-B)=a+bi+cj)雖然可以對研究對象所處狀態(tài)空間進行“一分為三”的刻畫，但是將研究對象所處的空間狀態(tài)簡單的“一分為三”略顯粗糙［15］，因此，必須根據(jù)原創(chuàng)聯(lián)系度的可展性原理進行深層次的細化，將原來的表達式改寫成式(1):

(1)細化后的b、c的意義

假設(shè)有n個評價對象，用m個指標描述評價對象屬性，xij為對象i的第j項指標(j=1，2，…，m;i=1，2，…，n)，將耕地占補平衡的生態(tài)安全等級劃分為K級(K=1，2，…，s)，文本中令s=4。每個年份的耕地占補平衡的生態(tài)安全評價指標值系列組成一個集合，將耕地占補平衡的生態(tài)安全評價指標的各個等級標準值組成另一個集合，這兩個集合于是構(gòu)成一個集對(Ai，Bk)［16］。若指標xij處于某一等級k的范圍內(nèi)時，則認為是同一，其值為a;若指標xij處于k級別相鄰級別優(yōu)越一邊，則認為是優(yōu)異，記為b1;若指標xij處于k級別相鄰級別劣差一邊，則認為是劣異，記為b2;若處于k級別相隔級別優(yōu)越一邊，則認為是優(yōu)反，記為c1;若處于k級別相隔級別劣差一邊，則認為是劣反，記為c2［17］。

(2)考慮優(yōu)劣識別

在耕地占補平衡的生態(tài)安全評價中，為了方便生態(tài)安全等級的優(yōu)劣識別，可以將式(1)簡記為:

進一步分析耕地占補平衡的生態(tài)安全評價指標分級標準和評價指標值之間的數(shù)量關(guān)系，可以看出，即使不同年份的生態(tài)安全評價指標值處于同一級別中，也會因為指標數(shù)值的差異而有所不同，所以，需要構(gòu)造聯(lián)系度表達式更細致的描述同一、差異和對立的定量關(guān)系。若指標值處于評價級別中，則a=1，其他系數(shù)為0;若指標值處于相鄰級別優(yōu)越的一邊，越靠近相鄰的評價標準a越大b越小，反之a(chǎn)越小b越大;若在評價等級劣差的一邊，越靠近評價標準a越大c越小，反之a(chǎn)越小c越大;若指標值處于相隔級別優(yōu)越的一邊，越靠近評價標準則a、b越大，d越小;若在評價等級劣差的一邊，越靠近評價標準則a、c越大，e越?。?7］。依據(jù)以上原理建立指標評價模型。對于正向性指標，相對于一級的聯(lián)系度:

二級聯(lián)系度:

三級聯(lián)系度:

四級聯(lián)系度:

1.2 基于改進SPA法的評價結(jié)果判定

根據(jù)式(3)得出各個年份對于不同等級的聯(lián)系度表達式后，本文擬用集對勢進行評價結(jié)果判定。

1.2.1 集對勢計算

集對勢的定義為同一度和對立度的比值，集對勢所反映兩個集合具有的某種趨同程度或聯(lián)系趨勢。在耕地占補平衡生態(tài)安全評價中計算年份i對于不同等級的集對勢，可得到年份i對于不同等級的趨同程度。當集對勢越大時，表示其同一度越強，所以取集對勢中的最大值所對應(yīng)的等級為研究對象的判定結(jié)果。改進的SPA法集對勢計算公式如下:

1.2.2 樂觀勢與悲觀勢計算

樂觀勢就是從樂觀的角度出發(fā)［18］，將所有的差異度全部轉(zhuǎn)化成同一度，在系統(tǒng)中，用同一度和對立度的比值來量化系統(tǒng)的態(tài)勢，樂觀勢的表達式為:

悲觀勢就是從悲觀的角度出發(fā)，將所有的差異度均轉(zhuǎn)化成對立度，在系統(tǒng)中，用同一度和對立度的比值來量化系統(tǒng)的態(tài)勢，悲觀勢的表達式為:

1.2.3 集對勢定位

由樂觀勢和悲觀勢所確定的區(qū)間范圍就是集對勢的取值范圍，依據(jù)式(7)作集對勢定位計算，計算集對勢在集對勢的取值范圍相對于樂觀勢的貼近度。

模型中skj為第j項指標評價標準的k級和k+1級的界限值，根據(jù)上述模型對每個年份的各個指標計算聯(lián)系度系數(shù) a、b、c、d、e的值，同理，可擬構(gòu)建負向型指標相對于不同等級的聯(lián)系度公式［17］。

設(shè)權(quán)重向量為W=(wij)n×m，聯(lián)系度矩陣運用模糊算子對W和R進行合成運算，得到評價年份i相對于K級的綜合關(guān)聯(lián)度Cik:

2 實例研究

2.1 研究區(qū)概況

永川區(qū)位于重慶市西南部，東與璧山縣、江津區(qū)為鄰，西以榮昌縣、大足縣為界，南與四川省合江縣、瀘縣相連，北與銅梁縣接壤。地跨28°56'― 29°34'N，105°38'―106°05'E。2009 年全區(qū)總?cè)丝?103438 人，其中農(nóng)業(yè)人口799507人，非農(nóng)業(yè)人口303931人，城鎮(zhèn)化率為27.54%，實現(xiàn)社會生產(chǎn)總值300.04億元，完成社會固定資產(chǎn)投資300.7億元。2009年，農(nóng)用地面積為130425hm2，占土地總面積的82.62%，其中，耕地67459hm2，占農(nóng)用地面積的51.72%;園地17954hm2，占農(nóng)用地面積的13.77%;林地28397hm2，占農(nóng)用地面積的2.77%。永川區(qū)屬于亞熱帶季風性濕潤氣候，平均氣溫18.20℃，，年平均降雨量1042.20mm，平均日照1298.50h，年平均無霜期317d。境內(nèi)地形復雜，丘陵起伏，其中相對高差小于50m的緩丘平壩面積為714.04hm2，占土地總面積的45.23%;丘陵面積577.33hm2，占其土地總面積的36.57%;相對高差大于200m的低山面積為287.24hm2，占其土地總面積18.20%。

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

人口、DGP、化肥施用量、農(nóng)藥使用量、糧食生產(chǎn)總量以及其他的自然地理狀況數(shù)據(jù)等來源于歷年《永川統(tǒng)計年鑒》;2005、2009年耕地資源相關(guān)數(shù)據(jù)均來源于永川區(qū)國土資源與房屋管理局;2015年和2020年永川區(qū)相關(guān)耕地資源數(shù)據(jù)則參考《2006-2020年重慶市永川區(qū)土地利用總體規(guī)劃》相關(guān)成果。

2.3 耕地占補平衡的生態(tài)安全評價指標體系及分級標準

在本研究界定的耕地占補平衡生態(tài)安全概念的基礎(chǔ)上，基于永川區(qū)自身的自然條件、社會經(jīng)濟發(fā)展特點，依據(jù)指標數(shù)據(jù)的可得性和方法的可操作性［19］，充分考慮耕地占補平衡生態(tài)安全各評價因子的復雜關(guān)系，借鑒國內(nèi)最新的相關(guān)研究成果［1，5，8］，結(jié)合國家環(huán)保總局制定的《生態(tài)縣、生態(tài)市、生態(tài)省建設(shè)指標(試行)》關(guān)于生態(tài)安全的目標，重點選取與耕地占補平衡生態(tài)安全密切相關(guān)的人均耕地面積、耕地減少面積與耕地面積比例、耕地增加面積與耕地面積比例、單位面積耕地化肥負荷以及單位面積耕地農(nóng)藥負荷等18個屬性特征構(gòu)建評價指標體系，并結(jié)合永川區(qū)歷年各指標量值，利用信息熵法計算出各指標權(quán)重，結(jié)果見表1。

表1 耕地占補平衡生態(tài)安全指標體系及安全等級劃分標準Table 1 Evaluation index system and different safe degrees of eco-security of cultivated land requisition－compensation balance

2.4 永川區(qū)耕地占補平衡的生態(tài)安全評價指標數(shù)據(jù)

永川區(qū)2015、2020年的自然、社會經(jīng)濟數(shù)據(jù)依據(jù)現(xiàn)有的數(shù)據(jù)利用線性回歸模型進行預(yù)測;2015年和2020年耕地相關(guān)數(shù)據(jù)則主要參考《2006—2020年重慶市永川區(qū)土地利用總體規(guī)劃》相關(guān)成果。2005、2009、2015和2020年耕地占補平衡生態(tài)安全評價指標值見表2。

2.5 權(quán)重確定

耕地占補平衡生態(tài)安全評價離不開指標權(quán)重的確定，權(quán)重確定方法的客觀與否直接關(guān)系到評價結(jié)果的好壞。目前確定指標權(quán)重的方法主要有層次分析法、特爾菲法、信息熵法。因此，為了提高指標權(quán)重確定的客觀性，本文采用信息熵法確定各項指標權(quán)重。按照熵的思想，人們能夠根據(jù)決策中獲取信息的數(shù)量和質(zhì)量，提高決策的精度和可靠性，熵在應(yīng)用于不同決策過程中的評價或案例的效果評價時是一個很理想的尺度［8］。其基本原理就是某項指標的值變異程度越大，信息熵越小，該指標提供的信息量就越大［20］，相應(yīng)指標權(quán)重就越大，反之則越小。具體計算過程如下:

表2 永川區(qū)2005、2009、2015和2020年耕地占補平衡生態(tài)安全評價指標值Table 2 Evaluation index value of eco-security of cultivated land requisition-compensation balance in Yongchuan District from 2005，2009，2015 and 2020

(1)評價指標值矩陣標準化［21］

根據(jù)表1實測值得到原始指標值矩陣X=(xij)8×7，通過對xij標準化處理可得R=(pij)mn，其中，pij為第i個評價單元第j評價指標值標準化后所得值，pij∈［0，1］。

(2)計算參評指標熵值

式中，H(xj)為第j項參評指標的熵值，為計算方便，通常令調(diào)劑系數(shù)k=1/ln＞0。(3)計算指標差異系數(shù)(hj)

第j項參評指標差異系數(shù)定義為:

(4)確定參評指標權(quán)重系數(shù)(Wj)

根據(jù)以上信息熵法的計算，得出永川區(qū)耕地占補平衡生態(tài)安全評價各指標權(quán)重分別為:W1=0.0038、W2=0.1500、W3=0.1980、W4=0.0061、W5=0.0093、W6=0.0798、W7=0.0052、W8=0.0023、W9=0.0429、W10=0.0022、W11=0.0305、W12=0.0021、W13=0.0726、W14=0.0371、W15=0.0618、W16=0.0574、W17=0.0429、W18=0.1960。

2.6 改進SPA法綜合評價計算

(1)SPA法聯(lián)系度矩陣R計算

這里以永川區(qū)2005年為例，計算其SPA法聯(lián)系度矩陣，結(jié)果如下:

同理可以計算出永川區(qū)2009、2015以及2020年的SPA法聯(lián)系度。

(2)進行合成運算

運用Excel中的sumproduct函數(shù)對權(quán)重向量W和永川區(qū)2005年聯(lián)系度矩陣R2005進行合成運算，得出永川區(qū)2005年對于不同等級的綜合聯(lián)系度。

根據(jù)式(4)—(7)計算2005年永川區(qū)相對于不同等級的集對勢、樂觀勢和悲觀勢等，利用集對勢進行耕地占補平衡生態(tài)安全等級判定，取集對勢三參數(shù)間特元［17］最大值者對應(yīng)的等級為最終評價結(jié)果，如下表3所示。

表3 永川區(qū)耕地占補平衡生態(tài)安全改進SPA法綜合評價結(jié)果(2005年)Table 3 Improved SPA comprehensive evaluation results of eco-security of cultivated land requisition-compensation balance in Yongchuan district(2005)

同理可以計算出永川區(qū)2009、2015和2020年耕地占補平衡生態(tài)安全改進SPA法綜合評價結(jié)果。為了驗證改進SPA法評價結(jié)果的合理性，本文與改進模糊綜合評價、多指標綜合評價結(jié)果進行比較(表4)。

從表4可以看出，利用改進SPA法得到永川區(qū)耕地占補平衡生態(tài)安全評價結(jié)果為:2005年耕地占補平衡生態(tài)安全級別為Ⅱ級，即“臨界安全”;2009年生態(tài)安全級別為Ⅲ級，即“較安全”;2015年生態(tài)安全級別為，即

“安全”;2020年生態(tài)安全級別同樣為Ⅳ級，即“安全”，整體上呈現(xiàn)出由Ⅱ級向Ⅳ級上升的趨勢。

表4 永川區(qū)耕地占補平衡生態(tài)安全綜合評價結(jié)果Table 4 Comprehensive evaluation results of eco-security of cultivated land requisition-compensation balance in Yongchuan district

將改進SPA法綜合評價結(jié)果與改進模糊綜合評價結(jié)果進行比較可知:2005年、2009年、2020年兩種方法評價結(jié)果完全一致;只有2015年改進SPA法綜合評價結(jié)果為Ⅳ級，即“安全”，而改進模糊綜合評價結(jié)果為Ⅲ級，即“較安全”，但是它們判定的耕地占補平衡生態(tài)安全等級相差不過一個等級，未出現(xiàn)越級現(xiàn)象。改進SPA法綜合評價結(jié)果與多指標綜合評價結(jié)果比較得知，雖然多指標綜合評價沒有進行評價分級，但是其計算的多指標綜合評價總分值從2005年至2020年是不斷增長的，這與改進SPA法綜合評價得到的耕地占補平衡生態(tài)安全等級在一定程度上正好呈對正應(yīng)關(guān)系。通過綜合分析可知:改進SPA法綜合評價結(jié)果與改進模糊綜合評價結(jié)果以及多指標綜合評價結(jié)果基本一致，因為改進模糊綜合評價法和多指標綜合評價法的實用性和可行性已得到普遍認可［17］，所以利用集對勢進行耕地占補平衡生態(tài)安全評價具有一定的可行性和可行性，評價結(jié)果比較貼近實際情況。

3 結(jié)論與討論

集對分析法通過注重信息處理過程中的模糊性和相對性比較好的解決了問題的不確定性，它是一種辯證思維的綜合評價模型。本文利用原創(chuàng)聯(lián)想度的可展性原理建立了耕地占補平衡生態(tài)安全改進SPA法評價模型，更為細致的描述了評價對象的同一性、差異性、對立性特性。此外，利用信息熵法表征各指標在耕地占補平衡生態(tài)安全評價中的重要程度，避免了評價指標權(quán)重確定略顯主觀性和隨意性，較為準確反應(yīng)耕地占補平衡生態(tài)安全等級狀態(tài)。研究結(jié)果表明:利用改進SPA法可以判定出，永川區(qū)2005年耕地占補平衡生態(tài)安全級別為Ⅱ級，2009年為Ⅲ級，2015年生態(tài)安全級別為Ⅳ級，2020年生態(tài)安全級別同樣為Ⅳ級，整體上呈現(xiàn)出由Ⅱ級向Ⅳ級上升的趨勢;改進SPA法綜合評價結(jié)果與改進模糊綜合評價結(jié)果以及多指標綜合評價評價結(jié)果基本一致，充分說明其在生態(tài)安全等級界定上具有一定的可行性和合理性。

在耕地占補平衡生態(tài)安全評價過程中，由于必須設(shè)定評價指標的等級狀況，而目前很多相關(guān)準則和標準中只有標準等級值，而不是一個域，有的甚至還沒有等級標準。此外，由于當前耕地占補平衡生態(tài)安全評價沒有公認的評價指標，而且還受到豐富多樣的物質(zhì)環(huán)境、文化環(huán)境以及自然生態(tài)環(huán)境的影響，文章尚無法對這些相關(guān)因素進行定量分析。因此，如何合理的界定評價指標的等級狀況以及建立一套更加全面的耕地占補平衡生態(tài)安全評價指標體系以及如何量化各指標的影響度等問題，仍需深入研究。

［1］He T B，Jin L，Deng D D.Early warning of cultivated land requisition-compensation balance in Karst mountainous county area.Transactions of the Chinese Society of Agricultural Engineering，2012，28(1):238-243.

［2］Li P，Shao J A，Zhang Z，Wei C F，Qiu D C.Establishment of Balance of CultivatedLand in Chongqing.Journal of Natural Resources，2011，26(6):919-931.

［3］Tan，Y Z，Wu C F，Wang Q R，Zhou L Q，Yan D.The Change of Cultivated Land and Ecological Environment Effects Driven by the Policy of Dynamic Equilibrium of the Total Cultivated Land.Journal of Natural Resources，2005，20(9):727-734.

［4］Solomon Barry.Farmland protection:A case of quality not quantity.Land Use Policy，1984，1(4):357-366.

［5］Huang H L，Luo W B，Wu C F，Li D M.Evaluation of land eco-security based on matter element analysis.Transactions of the Chinese Society of Agricultural Engineering，2010，26(3):316-322.

［6］Zhang J Y，Sun W C，Zhang F T.Regional land ecological security evaluation in the case of Chongqing Three Gorges Reservoir ecological economy area based on the PSR model.China Environmental Science，2011，31(6):1039-1044.

［7］Zeng Y.The regional ecological risk assessment of Hohhot City.Acta Ecologica Sinica，2010，30(3):668-673.

［8］Yu J，F(xiàn)ang L，Cang D B，Zhu L，Bian Z F.Evaluation of land eco-security in Wanjiang district base on entropy weight and matter element model.Transactions of the Chinese Society of Agricultural Engineering，2012，28(5):260-266.

［9］Li Z，Yang Z S.Research progress of land ecological security in China.China Safety Science Journal，2007，17(12):5-13.

［10］Sheng S K，She Y M.Approach to information systems security risk assessment based on fuzzy-BP neural network.Computer Simulation，2011，28(10):91-94.

［11］Wang S J，Yang Z F，Ding J.Projection pursuit method of comprehensive evaluation on groundwater carrying capacity in Guanzhong Plainl.Resources Science，2006，26(6):104-110.

［12］Wan X，Ding J，Zhang X L.Set pair analysis for comprehensive evaluation on regional groundwater resources carrying capacityl.Urban Environment＆ Urban Ecology，2006，19(2):8-10.

［13］Meng X M，Hu H P.Application of set pair analysis model based on entropy weight to comprehensive evaluation of water qualityl.Journal of Hydraulic Engineering，2009，40(3):258-262.

［14］Zhao K Q，Xuan A l.Set pair theory:A new theory method of non-define and its applicationsl.System Engineering，1996，14(1):18-23，72.

［15］Jiang Y L，Xu C F.Advances in set pair analysis theory and its applicationl.Computer Science，2006，33(1):205-209.

［16］Wei M H，Zheng Z H，Huang Q，Qiu L.Fuzzy comprehensive evaluation of groundwater environment based on improved set pair analysis.Journal of Hydraulic Engineering，2009，40(10):1204-1209.

［17］Nan C Y，Su X l.Comprehensive evaluation for soil-water resources carrying capacity in Guanzhong area Based on improved Set Pair Analysis.Journal of Natural Resources，2012，27(12):104-114.

［18］Li D S，Xu K L，Zhang X J，Dai，X S.Research and application of set pair situation in the set pair analysis.Industrial Safety and Environmental Protection，2009，35(9):8-9.

［19］Ye Y P，Liu L J.A preliminary study on assessment indicator system of provincial eco-environmental quality in China.Research o f Environmental Sciences，2002，13(13):33-36.

［20］Shi K F，Diao C T，Sun Y A，Ma J，Zhang B.Estimation of realistic potential of land consolidation in rural residential areas by combined weight on principle of relative entropy.Journal of Agricultural Mechanization Research，2012，34(7):14-18.

［21］Shi K F，Diao C T，Sun Y A，Sun X F，Zuo T A.Comparative on evaluation and classification of rural residentia1 land consolidation potential based on entropy weight method and AHP:A case study of Yongchuan District，Chongqing City.Journal of Agricultural Mechanization Research，2012，34(6):236-241.

［22］Shi K F，Diao C T.Calculation of realistic circulation potential of the rural collective construction land in Yongchuan district，Chongqing.Tropical Geography，2012，32(3):269-273.

參考文獻:

［1］何騰兵，金蕾，鄧冬冬.喀斯特山區(qū)縣域耕地占補平衡預(yù)警.農(nóng)業(yè)工程學報，2012，28(1):2438-243.

［2］李萍，邵景安，張貞，魏朝富，邱道持.重慶市耕地占補平衡體系構(gòu)建.自然資源學報，2011，26(6):919-931.

［3］譚永忠，吳次芳，王慶日，周煉清，嚴棟.“耕地總量動態(tài)平衡”政策驅(qū)動下中國的耕地變化及其生態(tài)環(huán)境效應(yīng).自然資源學報，2005，20(9):727-734.

［5］黃輝玲，羅文斌，吳次芳，李冬梅.基于物元分析的土地生態(tài)安全評價.農(nóng)業(yè)工程學報，2010，26(3):316-322.

［6］張軍以，蘇維詞，張鳳太.基于PSR模型的三峽庫區(qū)生態(tài)經(jīng)濟區(qū)土地生態(tài)安全評價.中國環(huán)境科學，2011，31(6):1039-1044.

［7］曾勇.區(qū)域生態(tài)風險評價:以呼和浩特市區(qū)為例.生態(tài)學報，2010，30(3):668-673.

［8］余健，房莉，倉定幫，朱琳，卞正富.熵權(quán)模糊物元模型在土地生態(tài)安全評價中的應(yīng)用.農(nóng)業(yè)工程學報，2012，28(5):260-266.

［9］李智國，楊子生.中國土地生態(tài)安全研究進展.中國安全科學學報，2007，17(12):5-13.

［10］申時凱，佘玉梅.模糊神經(jīng)網(wǎng)絡(luò)在信息安全風險評估中的應(yīng)用.計算機仿真，2011，28(10):91-94.

［11］王順久，楊志峰，丁晶.關(guān)中平原地下水資源承載力綜合評價的投影尋蹤方法.資源科學，2006，26(6):104-110.

［12］萬星，丁晶，張曉麗.區(qū)域地下水資源承載力綜合評價的集對分析方法.城市環(huán)境與城市生態(tài)，2006，19(2):8-10.

［13］孟憲萌，胡和平.基于熵權(quán)的集對分析模型在水質(zhì)綜合評價中的應(yīng)用.水利學報，2009，40(3):258-262.

［14］趙克勤，宣愛理.集對論:一種新的不確定理論方法與應(yīng)用.系統(tǒng)工程，1996，14(1):18-23，72.

［15］蔣云良，徐從富.集對分析理論及其應(yīng)用研究進展.計算機科學，2006，33(1):205-209.

［16］魏明華，鄭志宏，黃強，邱林.基于改進SPA法的地下水環(huán)境模糊綜合評判.水利學報，2009，40(10):1204-1209.

［17］南彩燕，粟曉玲.基于改進SPA的關(guān)中地區(qū)水土資源承載力綜合評價.自然資源學報，2012，27(12):104-114.

［18］李德順，許開立，張喜嘉，戴雪松.集對分析集對勢的研究及其應(yīng)用.工業(yè)安全與環(huán)保，2009，35(9):8-9.

［19］葉亞平，劉魯君.中國省域生態(tài)環(huán)境質(zhì)量評價指標體系研究.環(huán)境科學研究，2002，13(13):33-36.

［20］施開放，刁承泰，孫永愛，馬驥，張彪.基于熵組合權(quán)重的農(nóng)村居民點復墾現(xiàn)實潛力估算.農(nóng)機化研究，2012，34(7):14-18.

［21］施開放，刁承泰，孫永愛，孫秀鋒，左太安.不同確權(quán)法在農(nóng)村居民點復墾潛力評價中的應(yīng)用:以重慶市永川區(qū)為例.農(nóng)機化研究，2012，34(6):236-241.

［22］施開放，刁承泰.重慶市永川區(qū)農(nóng)村集體建設(shè)用地流轉(zhuǎn)現(xiàn)實潛力估算.熱帶地理，2012，32(3):269-273.