譚志衛(wèi)，董云仙，趙 潤

(云南省環(huán)境科學(xué)研究院，云南昆明650034;云南省高原湖泊流域污染過程與管理重點(diǎn)實(shí)驗(yàn)室，云南昆明650034)

近年來，生態(tài)系統(tǒng)健康成為生態(tài)系統(tǒng)管理領(lǐng)域的一個(gè)熱點(diǎn)問題［1］，客觀、科學(xué)、全面的評(píng)價(jià)是管理的前提和基礎(chǔ)，其作為資源環(huán)境管理的新方法，被眾多的學(xué)者重視，逐漸成為生態(tài)學(xué)研究的一個(gè)重要方向［2］。由于人口增加和工農(nóng)業(yè)生產(chǎn)的發(fā)展，人類賴以生存的生態(tài)環(huán)境日益退化，特別是湖泊水域富營養(yǎng)化問題，已嚴(yán)重影響并威脅到社會(huì)經(jīng)濟(jì)和人類健康［3］。對(duì)湖泊生態(tài)系統(tǒng)進(jìn)行研究，建立湖泊生態(tài)系統(tǒng)健康評(píng)價(jià)體系和指標(biāo)閾值，不但可以為我國湖泊環(huán)境綜合治理提供可操作的依據(jù)，而且還有利于湖泊生態(tài)系統(tǒng)的可持續(xù)性管理和合理利用，實(shí)現(xiàn)生態(tài)、社會(huì)、經(jīng)濟(jì)三方面效益的協(xié)調(diào)。

1 研究區(qū)概況

程海湖又名黑伍海，位于永勝縣中部，屬永勝縣程海鎮(zhèn)管轄。流域面積318.3km2，其中陸地面積 243.2km2，水域面積 75.8km2。湖南北長19.15km，東西最大寬度5.2km，平均寬度4.3km，湖岸線長45.1km，蓄水量19.8億m3，平均水深25.7m，最大水深35.1m。湖泊地處滇西北高原與云貴高原的銜接部位，世界自然遺產(chǎn)“三江并流”金沙江中段，屬中亞熱帶高原季風(fēng)氣候，主要盛行南風(fēng)，年平均氣溫19.1℃，最冷月平均氣溫8～11℃。

2 研究方法

2.1 數(shù)據(jù)來源及處理

程海流域 1∶25萬數(shù)字高程模型 (DEM);1989年、2000年、2012年程海流域土地利用狀況，其數(shù)據(jù)從這三期TM影像解譯獲取，數(shù)據(jù)格式為ArcGIS的Shp文件;1989年、2000年、2012年程海鎮(zhèn)統(tǒng)計(jì)年鑒及經(jīng)濟(jì)年鑒。

用ArcGIS軟件將Shp格式的程海流域土地利用/覆蓋圖處理成Coverage格式，并利用景觀格局分析軟件Fragstas3.3計(jì)算出3個(gè)時(shí)段的各類景觀結(jié)構(gòu)指數(shù)并進(jìn)行分析。

2.2 評(píng)價(jià)方法

二十多年來，程海流域隨著社會(huì)經(jīng)濟(jì)的發(fā)展及各種自然、地域、水資源利用等因素的影響，流域的土地利用發(fā)生了較大變化，從而使得流域的景觀格局發(fā)生變化，這必然影響到了區(qū)域的生態(tài)系統(tǒng)健康。人類活動(dòng)對(duì)生態(tài)系統(tǒng)的健康也尤為重要，人類活動(dòng)會(huì)威脅生態(tài)系統(tǒng)健康，導(dǎo)致生態(tài)系統(tǒng)結(jié)構(gòu)發(fā)生變化，從而影響生態(tài)系統(tǒng)的服務(wù)功能。因此，本次評(píng)估采用系統(tǒng)本身的生態(tài)指標(biāo)和區(qū)域的社會(huì)、經(jīng)濟(jì)、人類活動(dòng)等方面的指標(biāo)進(jìn)行評(píng)價(jià)。生態(tài)學(xué)指標(biāo)與經(jīng)濟(jì)學(xué)指標(biāo)相結(jié)合使指標(biāo)系統(tǒng)更為全面，評(píng)價(jià)結(jié)果更可靠。

2.3 生態(tài)系統(tǒng)評(píng)價(jià)指標(biāo)體系

在綜合考慮了國內(nèi)外有關(guān)流域生態(tài)系統(tǒng)健康評(píng)價(jià)指標(biāo)體系的基礎(chǔ)上，以景觀生態(tài)學(xué)、生態(tài)系統(tǒng)健康理論為基礎(chǔ)，結(jié)合程海流域當(dāng)?shù)氐膶?shí)際情況，從流域景觀格局、自然生態(tài)系統(tǒng)和人類活動(dòng)干擾3個(gè)方面，共選取了10項(xiàng)指標(biāo)來評(píng)價(jià)，如圖1。

(1)多樣性指數(shù)

景觀多樣性指數(shù)是借用生物多樣性提出的用來描述和評(píng)價(jià)景觀異質(zhì)性水平的一個(gè)概念，有多種不同的測度指標(biāo)。本次評(píng)估采用Shannon－Weaner指數(shù)。

式中:H是景觀多樣性指數(shù)，APi是第i類景觀要素面積占景觀總面積的比例，M指景觀要素類型個(gè)數(shù)。各類景觀要素面積相等時(shí)的景觀多樣性指數(shù)最大，用Hmax表示:

(2)景觀均勻度指數(shù)

景觀均勻度指數(shù)指實(shí)際多樣性和最大多樣性之比，是景觀多樣性的相對(duì)值。其計(jì)算公式為:

景觀多樣性可以定量地描述為景觀中景觀要素斑塊的不確定性，可以反映景觀異質(zhì)性和景觀中不同景觀類型分布的均勻化和復(fù)雜程度。由其計(jì)算公式可見，影響景觀多樣性指數(shù)大小的因素一是景觀中景觀要素類型的數(shù)量，這決定研究對(duì)象的生態(tài)學(xué)尺度和空間分辨率;二是各景觀要素類型間的面積分配關(guān)系。

(3)斑塊密度指數(shù)

斑塊密度包括景觀斑塊密度和景觀要素斑塊密度。景觀斑塊密度是指景觀中所有異質(zhì)性景觀要素中斑塊的單位面積斑塊數(shù);景觀要素斑塊面積則是指景觀要素中某類景觀要素的單位面積斑塊數(shù)。斑塊密度反映景觀的破碎化程度與景觀空間異質(zhì)性程度。斑塊密度越大，破碎化程度越大，空間異質(zhì)性程度也愈大，反之亦然。

式中:PD為景觀斑塊密度;m為景觀要素類型總數(shù);A為研究地區(qū)景觀總面積。

(4)土壤墾殖指數(shù)

土壤墾殖指數(shù)即耕地面積占區(qū)域總面積的百分比，表示土地持續(xù)提供人類生存資源的能力。該指數(shù)越大，土地的利用潛力越弱，可開發(fā)余地越小。

(5)土壤侵蝕指數(shù)

土壤侵蝕指數(shù)指輕度以上侵蝕面積占區(qū)域總面積的百分比。該指數(shù)越大，說明被破壞的土地面積越大，土地保水保肥能力差。

(6)森林覆蓋率

森林覆蓋率指林地占區(qū)域總面積的百分比。

3 綜合參數(shù)評(píng)價(jià)模型

為了避免人為主觀判斷的片面性，本評(píng)估選用客觀賦權(quán)法來進(jìn)行評(píng)價(jià)?？陀^賦權(quán)法中的熵值法是在客觀條件下，由評(píng)估指標(biāo)值構(gòu)成的判斷矩陣來確定指標(biāo)權(quán)重的一種方法，它能盡量消除各因素權(quán)重的主觀性，使評(píng)估結(jié)果更符合實(shí)際。

在信息論中，熵是無系統(tǒng)無序程度的度量，它還可以度量數(shù)據(jù)所提供的有效信息。信息熵越小，系統(tǒng)無序度越小，信息的效用值越大;信息熵越大，系統(tǒng)無序度越高，其信息的效用值越小。其計(jì)算步驟如下:

(1)構(gòu)建n個(gè)樣本m個(gè)評(píng)估指標(biāo)的判斷矩陣。

(2)將判斷矩陣歸一化處理，得到歸一化判斷矩陣B。

在實(shí)際決策中，評(píng)估指標(biāo)通常分為越大越優(yōu)型和越小越優(yōu)型兩類，各類指標(biāo)對(duì)優(yōu)的相對(duì)隸屬度計(jì)算公式分別如下:

越大越優(yōu)型:

越小越優(yōu)型:

式中:xmax，xmin—同指標(biāo)下不同樣本中最滿意或最不滿意者 (越小越滿意或越大越滿意)。

根據(jù)式①和式②可將評(píng)估指標(biāo)特征值矩陣轉(zhuǎn)換為其對(duì)優(yōu)的相對(duì)隸屬度矩陣:

R=rij

式中:rij—方案j指標(biāo)i的特征值對(duì)優(yōu)的相對(duì)隸屬度。

(3)根據(jù)熵的定義，對(duì)于n個(gè)樣本m個(gè)評(píng)估指標(biāo)，可確定評(píng)估指標(biāo)的熵為:

式中:0≤Hi≤1，為使LNfij有意義，假定fij=0，fijLNfij=0，i=1，2，3……m，j=1，2，3……n。

(4)利用熵值計(jì)算評(píng)估指數(shù)的熵權(quán):

(5)生態(tài)系統(tǒng)健康綜合指數(shù):

式中:EHCI—生態(tài)系統(tǒng)健康綜合指數(shù)，其值在0與1之間;

Wi—評(píng)估指標(biāo)在綜合評(píng)估指標(biāo)體系中的權(quán)重值，其值在0與1之間;

Ii—評(píng)估指標(biāo)的歸一化值，其值在0與1之間。

4 生態(tài)系統(tǒng)健康等級(jí)劃分

本研究參考了國內(nèi)外相關(guān)研究的有關(guān)標(biāo)準(zhǔn)，結(jié)合程海流域的實(shí)際情況，提出生態(tài)系統(tǒng)健康評(píng)估標(biāo)準(zhǔn)，把生態(tài)系統(tǒng)健康狀態(tài)分為5個(gè)等級(jí)，取值范圍分別為很健康 (80～100)、健康 (60～80)、亞健康 (40～60)、不健康 (20～40)、病態(tài) (0～20)。各級(jí)具體含義見表1。

5 結(jié)果與分析

本評(píng)估選取1989年、2000年、2012年的結(jié)合度指數(shù)、多樣性指數(shù)、均勻度指數(shù)、斑塊密度指數(shù)、人口密度、人均GDP、城鎮(zhèn)化水平、土壤墾殖水平、土壤侵蝕指數(shù)、森林覆蓋率10個(gè)指標(biāo)進(jìn)行客觀評(píng)估。結(jié)合度指數(shù)、多樣性指數(shù)、均勻度指數(shù)、斑塊密度指數(shù)采用景觀生態(tài)學(xué)的方法，在Arcview3.3和Fragstat3.3的結(jié)合下完成了研究區(qū)景觀格局的分析，得出相關(guān)指數(shù)值。人口密度、人均GDP、城鎮(zhèn)化水平、土壤墾殖水平、土壤侵蝕指數(shù)、森林覆蓋率參考研究年年鑒及土地利用現(xiàn)狀數(shù)據(jù)得出。各指標(biāo)具體數(shù)值見表2。

表1 綜合生態(tài)系統(tǒng)健康等級(jí)劃分

表2 程海陸生生態(tài)系統(tǒng)健康綜合評(píng)估指標(biāo)體系

(1)首先構(gòu)建程海陸生生態(tài)系統(tǒng)健康評(píng)估3個(gè)樣本9個(gè)指標(biāo)的判斷矩陣R:

(2)對(duì)判斷矩陣R進(jìn)行歸一化處理，得到歸一化判斷矩陣B:

(3)根據(jù)熵的定義，可確定評(píng)估指標(biāo)的熵:

(4)計(jì)算評(píng)估指標(biāo)的熵權(quán):

(5)計(jì)算程海陸生生態(tài)系統(tǒng)健康綜合指數(shù)

將表3評(píng)估指標(biāo)的歸一化值及熵權(quán)代入式⑤，可得程海陸生生態(tài)系統(tǒng)的健康綜合指數(shù)，見表4。

表3 程海陸生生態(tài)系統(tǒng)健康綜合評(píng)估指標(biāo)歸一化值

表4 程海陸生生態(tài)系統(tǒng)健康綜合指數(shù)

表5 程海陸生生態(tài)系統(tǒng)健康評(píng)估結(jié)果

從評(píng)價(jià)結(jié)果可以看出，1989年、2000年、2012年的陸生生態(tài)系統(tǒng)健康狀態(tài)都為不健康狀態(tài)。陸生生態(tài)系統(tǒng)景觀被破壞，系統(tǒng)結(jié)構(gòu)出現(xiàn)缺陷，活力較低，外界壓力大，不能滿足生態(tài)系統(tǒng)功能正常發(fā)揮的需要。究其原因，主要是由于一直以來程海流域森林覆蓋率偏低，遠(yuǎn)低于云南省森林覆蓋率(52.93%)的平均水平;另一方面，程海流域位于金沙江干熱河谷地帶，屬中亞熱帶高原季風(fēng)氣候，全年受熱帶大陸性氣團(tuán)和海洋東南季風(fēng)影響，降雨不豐富，加上當(dāng)?shù)氐某赡纲|(zhì)主要為玄武巖、石灰?guī)r風(fēng)化土，雨熱條件和土壤結(jié)構(gòu)造成了程海流域植被類型少。流域主要喬木植被為云南松林和華山松林，而廣泛分布的是灌草叢和中山濕性常綠櫟類灌叢，樹種結(jié)構(gòu)不合理，群落結(jié)構(gòu)簡單，生物多樣性降低，保水保土能力差，不利于生態(tài)系統(tǒng)的良好發(fā)展;第三，程海流域長期以來社會(huì)經(jīng)濟(jì)發(fā)展以傳統(tǒng)農(nóng)業(yè)為主，由于當(dāng)?shù)厮O(shè)施不健全，降雨偏少，使靠天吃飯成為當(dāng)?shù)刂饕纳a(chǎn)方式，粗放的生產(chǎn)生活方式給生態(tài)環(huán)境帶來了人為壓力和脅迫;第四，從土地利用景觀格局動(dòng)態(tài)變化來看，20年來，在政府針對(duì)林業(yè)“宜封則封，宜造則造，封山育林與人工造林相結(jié)合”的強(qiáng)有力政策下，森林覆蓋率有所提高，同時(shí)區(qū)域內(nèi)的景觀結(jié)合度越來越好，有利于生物物質(zhì)的交流，但是總體水平不高，影響整個(gè)生態(tài)系統(tǒng)物質(zhì)和能量的交流。雖然程海流域三個(gè)時(shí)段的生態(tài)系統(tǒng)健康狀態(tài)都為不健康，但是2012年的健康指數(shù)比2000年和1989年都高，說明生態(tài)系統(tǒng)健康狀態(tài)有所改善，森林覆蓋率的提高、水土流失面積的減少對(duì)其也起到了印證作用。

6 措施建議

(1)繼續(xù)實(shí)施“宜封則封，宜造則造，封山育林與人工造林相結(jié)合”的政策，加快林分改造，增強(qiáng)水源林的“清水產(chǎn)流”功能。封山育林簡單易行，經(jīng)濟(jì)省事，可為鄉(xiāng)土樹種創(chuàng)造適宜的生境，促使林木生長，進(jìn)而演替為地帶性植被——中山濕性常綠闊葉林。通過封山育林、自然演替和林分改造提高演替速率。

(2)加強(qiáng)對(duì)土地使用方式的管理，保護(hù)生態(tài)用地資源。保護(hù)耕地資源的數(shù)量穩(wěn)定，嚴(yán)守“耕地紅線”。提高土地集約利用水平，避免對(duì)生態(tài)保護(hù)用地和耕地的侵占。

(3)優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)，“一調(diào)二控三促”。調(diào)整第一產(chǎn)業(yè)的產(chǎn)業(yè)結(jié)構(gòu)，改變傳統(tǒng)農(nóng)業(yè)的生產(chǎn)方式，發(fā)展生態(tài)農(nóng)業(yè)和有機(jī)農(nóng)業(yè)，加大農(nóng)業(yè)面源污染治理力度?？刂啤奥菪濉别B(yǎng)殖業(yè)的發(fā)展速度，并且加強(qiáng)“螺旋藻”養(yǎng)殖業(yè)污染的治理。促進(jìn)第三產(chǎn)業(yè)的發(fā)展。

(4)加強(qiáng)水土流失治理。加強(qiáng)程海流域境內(nèi)小流域治理力度，減少水土流失。

7 結(jié)語

本文由于資料有限，僅對(duì)三個(gè)時(shí)段進(jìn)行了健康評(píng)價(jià)，在今后的理論研究和實(shí)踐中應(yīng)補(bǔ)充其它時(shí)間序列。

［1］羅躍初，周忠軒，孫軼，等．流域生態(tài)系統(tǒng)健康評(píng)價(jià)方法［J］．生態(tài)學(xué)報(bào)，2003，23(8):1066－1075．

［2］方慶，董增川，劉晨，等．基于景觀格局的區(qū)域生態(tài)系統(tǒng)健康評(píng)價(jià)—以灤河流域行政區(qū)為例［J］．南水北調(diào)與水利科技，2012，10(6):37－41．

［3］濮培民．健康水生態(tài)系統(tǒng)的退化及其修復(fù)—理論、技術(shù)及應(yīng)用 ［J］．湖泊科學(xué)，2001，13(3):193－203．