孫意菲,孫意翔
(1.黑龍江省慶達(dá)水利水電工程有限公司,哈爾濱150080,2.河海大學(xué)水利水電學(xué)院,南京210098)
生態(tài)彈性度是指生態(tài)在自然和人為干擾下表現(xiàn)出的響應(yīng),用來(lái)反映生態(tài)穩(wěn)定性與生態(tài)環(huán)境問(wèn)題的可能性大小。在自然狀況下,各種生態(tài)過(guò)程維持著一種相對(duì)穩(wěn)定的耦合關(guān)系,保證著生態(tài)系統(tǒng)的相對(duì)穩(wěn)定,而當(dāng)外界干擾超過(guò)一定限度時(shí),這種耦合關(guān)系將被打破,適度干擾可增加其異質(zhì)性;而嚴(yán)重干擾往往使異質(zhì)性迅速降低。所以,生態(tài)彈性度評(píng)價(jià)實(shí)質(zhì)上是對(duì)現(xiàn)狀自然環(huán)境背景下潛在的生態(tài)穩(wěn)定進(jìn)行明確的辨識(shí),并將其落實(shí)到具體的空間區(qū)域,而目前將濕地作為研究對(duì)象的生態(tài)彈性度并不多見(jiàn)。
三江平原是國(guó)家重要的商量糧基地,大面積開(kāi)荒、墾建脫節(jié)、重用輕養(yǎng),使人類(lèi)活動(dòng)對(duì)其自然環(huán)境的干擾強(qiáng)度日益增大,導(dǎo)致生態(tài)環(huán)境的退化[1]。因此,開(kāi)展?jié)竦厣鷳B(tài)彈性度的評(píng)價(jià)是非常必要的。本文根據(jù)國(guó)內(nèi)外已有研究成果及三江平原濕地的調(diào)研數(shù)據(jù),嘗試?yán)媚:镌椒▽?duì)其進(jìn)行評(píng)價(jià),并在GIS支持下,與其分析屬性功能結(jié)合,對(duì)濕地生態(tài)彈性度進(jìn)行區(qū)劃,研究生態(tài)穩(wěn)定性的分布規(guī)律及其地域差異,以期為流域生態(tài)環(huán)境建設(shè)和保護(hù)規(guī)劃方案的制定提供依據(jù)。
三江平原位于黑龍江省東北部,其范圍東起烏蘇里江、西接湯旺河及牡丹江分水嶺,南至大興凱湖及綏芬河分水嶺,北界黑龍江,包括黑龍江、松花江與烏蘇里江匯流的三角地帶以及倭肯河與穆棱河流域和興凱湖平原。南北長(zhǎng)500 km,東西寬430 km,全區(qū)總控制面積10.88 km2,其中山區(qū)面積為3.74萬(wàn)km2,占總面積的 34.4%;丘陵面積為0.94萬(wàn)km2,占總面積的8.6%;平原面積為6.2萬(wàn)km2,占總面積的57%。該區(qū)西南高東北低,除西部和西南部邊界的小興安嶺、老爺嶺、張廣才嶺和橫亙中部的完達(dá)山為森林覆蓋的山區(qū)外,廣闊的沖積低平原和河流形成的階地、河漫灘上廣泛發(fā)育著沼澤和沼澤化草甸區(qū)。
三江平原現(xiàn)有河流、湖泊、庫(kù)塘、沼澤和沼澤化草甸類(lèi)型濕地,總面積為151.78萬(wàn)hm2,占三江平原地區(qū)的總面積的13.95%,其中可劃分為蘿北地區(qū)、同撫地區(qū)、撓力河地區(qū)、安邦河地區(qū)、倭肯河地區(qū)、穆棱河地區(qū)見(jiàn)圖1。
2.1.1 模糊物元和復(fù)合模糊物元[2-3]
給定事物的名稱(chēng)N,它關(guān)于特征c的量值為v,以有序三元R=(N,c,v)組作為描述事物的基本元,簡(jiǎn)稱(chēng)物元。事物名稱(chēng)N,特征c和量值v稱(chēng)為物元的三要素。如果其中量值具有模糊性,便稱(chēng)為模糊物元。如果事物N有n個(gè)特征c1,c2,···,cn和相應(yīng)的模糊量值 v1,v2,···,vn,稱(chēng) R 為n 維模糊物元,簡(jiǎn)記為 R=(N,C,V)。
圖1 三江平原濕地劃分區(qū)
如果m個(gè)事物的n維物元組合在一起,便構(gòu)成m個(gè)事物n維復(fù)合物元,記作Rmn。若將式(1)的量值改為模糊物元量值,稱(chēng)為m個(gè)事物n維復(fù)合模糊物元,記作R-mn,見(jiàn)式(2)。
式中:Rmn為m個(gè)事物n維復(fù)合物元;R-mn為m個(gè)事物n維復(fù)合模糊物元;Mi為第 i個(gè)事物,i=1,2,···m;ck為第 k項(xiàng)特征,k=1,2,···n;Xik為 第i個(gè)事物第k項(xiàng)特征對(duì)應(yīng)的量值,i=1,2,···m ,k=1,2,···n;uik為 第 i個(gè)事物第 k項(xiàng)特征對(duì)應(yīng)的模糊量值,i=1,2,···m ,k=1,2,···n 。
本文采用從優(yōu)隸屬度原則對(duì)uik進(jìn)行標(biāo)準(zhǔn)化,即采用以下指標(biāo):
越大越優(yōu)型:
越小越優(yōu)型:
式中:uik為第i個(gè)事物第k項(xiàng)特征對(duì)應(yīng)的模糊量值,i=1,2,···m ,k=1,2,···n;Xik為第 i個(gè)事物第 k項(xiàng)特征對(duì)應(yīng)的量值,i=1,2,···m ,k=1,2,···n;maxXik、minXik各事物中每一項(xiàng)特征所有量值Xik中的最大值和最小值。
2.1.2 標(biāo)準(zhǔn)模糊物元與差平方復(fù)合模糊物元
由式(2)可以構(gòu)成標(biāo)準(zhǔn)方案的n維模糊物元R0n,其中各項(xiàng)由R-mn中各方案從優(yōu)隸屬度中的最大值或最小值確定,則得到
若以Δji(i=1,2,…,n;j=1,2,… ,m)表示標(biāo)準(zhǔn)模糊物元R0n與復(fù)合模糊物元R-mn中各項(xiàng)差的平方,則組成差平方復(fù)合模糊物元RΔ,即
其中Δji=(u0i-uji)2,i=1,2,…,n;j=1,2,…,m。
2.1.3 權(quán)重的確定
本文采用熵權(quán)系數(shù)法來(lái)確定各評(píng)價(jià)指標(biāo)間的相對(duì)重要性次序,從而得到各評(píng)價(jià)指標(biāo)的權(quán)重,具體方法如下:
(1)計(jì)算第i個(gè)評(píng)價(jià)指標(biāo)下第j個(gè)待評(píng)價(jià)監(jiān)測(cè)點(diǎn)評(píng)價(jià)指標(biāo)特征值比重:
(2)計(jì)算第i個(gè)評(píng)價(jià)指標(biāo)的熵:
(3)計(jì)算第i個(gè)評(píng)價(jià)指標(biāo)的權(quán)重:
2.1.4 歐氏貼近度
由于本文具有綜合評(píng)價(jià)的意義,采用M(·,+)算法,即先乘后加運(yùn)算歐氏貼近度ρHj,則
式中:ρHj為第n個(gè)方案與標(biāo)準(zhǔn)方案之間的相互接近程度,其值越大,表示兩者越接近;反之,則相差越大。然后以此來(lái)構(gòu)造歐氏貼近度復(fù)合模糊物元RρH,則
歐氏貼近度是表示各方案與標(biāo)準(zhǔn)方案(最優(yōu)方案)之間的貼近程度,可以根據(jù)歐氏貼近度的大小對(duì)各方案進(jìn)行優(yōu)劣排序,并可進(jìn)行分類(lèi)。
傳統(tǒng)的評(píng)價(jià)方法,通??紤]的是時(shí)間因素,主要從時(shí)間上、區(qū)域經(jīng)濟(jì)發(fā)展的歷史、現(xiàn)狀出發(fā)來(lái)預(yù)測(cè)、規(guī)劃及確定經(jīng)濟(jì)的未來(lái)發(fā)展方向,主要有時(shí)間序列分析法、數(shù)量經(jīng)濟(jì)分析法、線(xiàn)性規(guī)劃及一些經(jīng)濟(jì)統(tǒng)計(jì)分析模型等。事實(shí)上,在經(jīng)濟(jì)增長(zhǎng)及發(fā)展過(guò)程中,人類(lèi)和自然活動(dòng)的空間布局變化總是和時(shí)間演變過(guò)程交織在一起,區(qū)域間差異不可避免地存在。專(zhuān)家已逐漸意識(shí)到在評(píng)價(jià)與決策中,必須考慮區(qū)域評(píng)價(jià)差異的空間特征[4]。GIS提供了一種更好的理解各規(guī)劃區(qū)之間的相互作用,使得區(qū)域規(guī)劃的實(shí)施得到更好保證。借助GIS技術(shù)分析生態(tài)問(wèn)題可解決:圖解各種環(huán)境要素位置;確定環(huán)境要素空間變化趨勢(shì);匹配滿(mǎn)足環(huán)境要素經(jīng)濟(jì)布局條件或要求場(chǎng)所;分析及可視化環(huán)境要素空間變化模式;模擬環(huán)境要素空間變化;制定適合優(yōu)化環(huán)境的經(jīng)濟(jì)發(fā)展規(guī)劃[5]。
將GIS與評(píng)價(jià)模型實(shí)現(xiàn)集成主要包括以下一些工作[6]:①將來(lái)自于GIS的空間數(shù)據(jù)和屬性數(shù)據(jù)用于模型的參數(shù)估計(jì)、率定及模型的檢驗(yàn);②以GIS為主要技術(shù)平臺(tái),實(shí)現(xiàn)模擬結(jié)果的可視化;③采用GIS技術(shù)來(lái)定義表達(dá)模型研究的區(qū)域單元。本研究是在ArcGIS軟件環(huán)境下,針對(duì)濕地生態(tài)問(wèn)題選擇若干影響因子,轉(zhuǎn)換為單因子?xùn)鸥駭?shù)據(jù)GRID。通過(guò)綜合各影響因子進(jìn)行評(píng)價(jià),采用Spatial Analyst模塊地圖代數(shù)分析功能得出各生態(tài)穩(wěn)定的彈性度分布圖。
利用上述模糊物元評(píng)價(jià)模型,模型參數(shù)包括年平均降水量、天然植被面積、調(diào)蓄洪水、濕地面積占本地區(qū)百分比、地表徑流、地下水(可開(kāi)采量)等影響因素,結(jié)合GIS對(duì)三江平原濕地生態(tài)彈性度進(jìn)行評(píng)價(jià),以檢驗(yàn)本文所建模型的實(shí)用性。表1中給出了濕地生態(tài)彈性度評(píng)價(jià)中各地區(qū)及各要素包括的評(píng)價(jià)指標(biāo)。
根據(jù)前述的評(píng)價(jià)方法,按照各地區(qū)情況選用年平均降水量、天然植被面積、調(diào)蓄洪水、濕地面積占本地區(qū)百分比、地表徑流、地下水(可開(kāi)采量)的評(píng)價(jià)指標(biāo)的樣本構(gòu)造復(fù)合物元;分別計(jì)算出各地區(qū)的從優(yōu)隸屬度,指標(biāo)按越大越優(yōu)計(jì)算;利用各評(píng)價(jià)樣本模糊物元與標(biāo)準(zhǔn)樣本模糊物元指標(biāo)之間差的平方確定差平方復(fù)合物元,這里取最大值組成標(biāo)準(zhǔn)樣本模糊物元;再利用熵權(quán)法確定指標(biāo)的權(quán)重值(表2);根據(jù)公式2-10得出歐式貼近度(表3),確定濕地生態(tài)穩(wěn)定度;最后通過(guò)GIS技術(shù)顯示評(píng)價(jià)結(jié)果空間特征(圖2)。
表1 三江平原生態(tài)彈性度定量性數(shù)據(jù)
表2 指標(biāo)權(quán)重值
表3 三江平原濕地生態(tài)穩(wěn)定度評(píng)價(jià)結(jié)果
圖2 濕地生態(tài)穩(wěn)定分級(jí)圖
從表3計(jì)算結(jié)果可以看出,同撫地區(qū)的生態(tài)彈性度為最大,說(shuō)明該地區(qū)生態(tài)系統(tǒng)穩(wěn)定程度較高、植被情況良好,沼澤濕地保護(hù)措施得當(dāng)。生態(tài)系統(tǒng)在抵御外界干擾能力和自我恢復(fù)能力還是比較強(qiáng)的,基本上處于穩(wěn)定狀態(tài),但是對(duì)于安邦河和倭肯河兩個(gè)地區(qū),彈性度比較低,處于較穩(wěn)定狀態(tài),但已接近不穩(wěn)定狀態(tài),尤其是安邦河地區(qū)生態(tài)彈性度只有0.387856,說(shuō)明如果再對(duì)該地區(qū)進(jìn)行開(kāi)發(fā),那么生態(tài)系統(tǒng)將遭到進(jìn)一步破壞,其后果將不堪設(shè)想,應(yīng)該引起足夠重視。
1)本文從生態(tài)彈性度的概念出發(fā),論述生態(tài)穩(wěn)定性及評(píng)價(jià)方法,提出了模糊物元評(píng)價(jià)模型在生態(tài)彈性度評(píng)價(jià)中的應(yīng)用,并對(duì)其方法進(jìn)行介紹。
2)評(píng)價(jià)結(jié)果得出:同撫地區(qū)的生態(tài)彈性度為最大,說(shuō)明該地區(qū)生態(tài)系統(tǒng)穩(wěn)定程度較高、植被情況良好,沼澤濕地保護(hù)措施得當(dāng)。安邦河地區(qū)生態(tài)彈性度最低,說(shuō)明如果再對(duì)該地區(qū)進(jìn)行開(kāi)發(fā),那么生態(tài)系統(tǒng)將遭到進(jìn)一步破壞,其后果將不堪設(shè)想,應(yīng)該引起足夠重視。
3)本文通過(guò)GIS對(duì)空間信息表達(dá),能夠發(fā)現(xiàn)區(qū)域生態(tài)不穩(wěn)定的規(guī)律,更直觀的體現(xiàn)區(qū)域間生態(tài)穩(wěn)定的空間變化特征,可以作為模型應(yīng)用中分析理論的基礎(chǔ)。
[1]李偉業(yè).三江平原沼澤濕地生態(tài)承載力與可持續(xù)調(diào)控模式研究[D].哈爾濱:東北農(nóng)業(yè)大學(xué)碩士論文,2007.
[2]蔡文.物元模型及其應(yīng)用[M].北京:科學(xué)技術(shù)文獻(xiàn)出版社,1994.
[3]張斌,雍歧東,肖芳淳.模糊物元分析[M].北京:石油工業(yè)出版社,1997.
[4]第寶鋒,楊忠,艾南山等.基于RS與GIS的金沙江干熱河谷區(qū)退化生態(tài)系統(tǒng)評(píng)價(jià)——以云南省元謀縣為例[J].地理科學(xué),2005,25(4):484 -489.
[5]DrummondW J.GIS as a visualization tool for economic development[J].Comput,Environ and Urban Systems,1993,17(5):469-479.
[6]侯英姿,陳曉玲,王方雄.基于GIS的水環(huán)境價(jià)值模糊綜合評(píng)價(jià)研究[J].地理科學(xué),2008,28(1):89-93.