徐浩田,周林飛,成 遣

沈陽農(nóng)業(yè)大學(xué)水利學(xué)院, 沈陽 110866

濕地作為地球上十分重要的自然資源和生態(tài)系統(tǒng),在調(diào)節(jié)徑流、凈化環(huán)境、控制污染、涵養(yǎng)水源、抵御洪水、保護(hù)生物多樣性和維持區(qū)域生態(tài)平衡等方面起著其他類型生態(tài)系統(tǒng)不可替代的作用[1- 2]。然而,隨著人類活動的不斷加劇,人類的需求發(fā)展同自然資源、環(huán)境的相互矛盾日趨尖銳,濕地資源受到了嚴(yán)重的干擾和破壞,同時導(dǎo)致了一系列生態(tài)環(huán)境問題的發(fā)生[3]。對濕地進(jìn)行生態(tài)系統(tǒng)健康評價與預(yù)警研究,不僅可以評估濕地資源所面臨的各種環(huán)境壓力,分析其現(xiàn)在所處的狀態(tài),還能夠預(yù)測濕地的演變規(guī)律,提出保護(hù)修復(fù)措施,從而確保濕地及其資源的可持續(xù)利用,為流域生態(tài)環(huán)境保護(hù)、工農(nóng)業(yè)生產(chǎn)及協(xié)調(diào)區(qū)域經(jīng)濟發(fā)展提供科學(xué)依據(jù),具有重要意義[4- 5]。

1 研究區(qū)概況與基礎(chǔ)資料

1.1 研究區(qū)概況

凌河口濕地(40°45′—41°00′N,121°00′—121°30′E)自然資源豐富,生態(tài)環(huán)境美好,位于渤海遼東灣區(qū)域的北部,包含遼寧省凌海市約69 km的海岸線,是遼河流域主要的濕地保護(hù)對象。凌河口濕地自然保護(hù)區(qū)擁有近8.36萬hm2的面積,其生態(tài)系統(tǒng)屬于典型的沿海濕地混合生態(tài)系統(tǒng)。凌河口受遼東灣北部潮汐、海流、波浪等海洋因素的影響,形成有大面積的灘涂、沼澤濕地。由于生態(tài)環(huán)境適宜,濕地內(nèi)部動植物物種豐富,約有239科,1024種的動植物,每年經(jīng)過凌河口濕地進(jìn)行遷徙的鳥類約有7萬余只,因此這里是國際上非常重要的候鳥遷徙停歇之地,也是目前世界上瀕危鳥類中特有的黑嘴鷗(Larus saundersi)和丹頂鶴(Grus japonensis)的主要棲息地和繁殖場所[6- 7]。

1.2 基礎(chǔ)資料

本研究的主要數(shù)據(jù)源是USGS官網(wǎng)提供的衛(wèi)星遙感影像,它們是1995年、2000年、2005年、2009年和2014年的分辨率為30 m的單幅TM衛(wèi)星遙感影像。研究過程中借助地理信息系統(tǒng)軟件ArcGIS為信息技術(shù)平臺,同時利用ERDAS 8.6軟件對獲取的影像依次進(jìn)行了輻射校正(輻射定標(biāo)和大氣校正)、幾何糾正、信息增強等影像預(yù)處理。其它收集的資料包括凌河口濕地研究區(qū)的水系圖、周邊行政界限圖和大凌河地區(qū)以往的調(diào)查資料以及與研究區(qū)相關(guān)的地質(zhì)地貌、水文氣象和社會經(jīng)濟資料等。

根據(jù)遼寧省土地利用分類系統(tǒng)劃分及對凌河口濕地進(jìn)行的實地調(diào)查資料,將研究區(qū)土地利用景觀類型劃分為八類,它們分別是：居民地、林地、蘆葦沼澤、水田、湖泊-養(yǎng)殖塘、旱地、灘涂及河流。運用目視解譯信息提取方法提取研究區(qū)地物信息,建立研究區(qū)主要解譯標(biāo)識,并通過GPS野外定點驗證解譯標(biāo)識的準(zhǔn)確性,最后運用ArcGIS 10.0軟件中的統(tǒng)計分析工具(Summary statistics)對校正好的影像進(jìn)行分析計算,最后得出研究區(qū)5個時期的遙感信息數(shù)據(jù)(表1)。

2 凌河口濕地生態(tài)系統(tǒng)健康評價

2.1 PSR概念模型的選取

壓力-狀態(tài)-響應(yīng)模型(Pressure-State-Response,簡稱PSR)最初是由Tony Friend和David Rapport提出的,用來分析研究自然環(huán)境所處的壓力、現(xiàn)狀與響應(yīng)之間的關(guān)系。在20世紀(jì)70年代,國際經(jīng)濟合作與發(fā)展組織(OECD)對PSR模型進(jìn)行了進(jìn)一步的修改并將其應(yīng)用于環(huán)境報告研究中；到90年代初期,OECD在對生態(tài)環(huán)境相關(guān)重要指標(biāo)進(jìn)行分析研究時,對PSR模型進(jìn)行了適用性和有效性的分析[8]。PSR模型就“發(fā)生了什么、為什么發(fā)生、我們將如何做”這3個人與自然和諧共處的基本問題給出了具體答案,尤其是其創(chuàng)造的將評價對象自身所面臨的相關(guān)壓力-狀態(tài)-響應(yīng)指標(biāo)和規(guī)范的參照標(biāo)準(zhǔn)進(jìn)行對比的分析模式,得到了國內(nèi)外眾多學(xué)者的肯定,同時被廣泛地應(yīng)用在區(qū)域環(huán)境、水文水資源及天然濕地自然資源等具體環(huán)境指標(biāo)體系的研究中[9]。

表1 凌河口濕地屬性信息

采用PSR模型進(jìn)行分析的優(yōu)點是：研究區(qū)域的各個環(huán)境影響因子相互之間的邏輯關(guān)系明確,外界人為因素壓力的干擾、研究區(qū)內(nèi)部各個環(huán)境指標(biāo)的變化和響應(yīng)的具體措施得到了充分的分析。依據(jù)凌河口濕地的勘察資料和研究區(qū)具體內(nèi)部情況,將PSR模型引入到研究區(qū)生態(tài)系統(tǒng)健康評價中,構(gòu)建凌河口濕地PSR基礎(chǔ)框架模型(圖1)。

圖1 凌河口濕地PSR模型Fig.1 PSR model in Linghekou Wetland

2.2 評價指標(biāo)體系的構(gòu)建

對濕地生態(tài)系統(tǒng)健康進(jìn)行評價,要將生態(tài)環(huán)境、經(jīng)濟影響和社會因素綜合考慮進(jìn)來,同時還應(yīng)顧及天然濕地生態(tài)結(jié)構(gòu)、經(jīng)濟體系以及社會構(gòu)成在不同管理政策條件下所產(chǎn)生的動態(tài)變化,從而有利于維持濕地系統(tǒng)的持續(xù)性[10]。

本研究依據(jù)凌河口天然濕地生態(tài)系統(tǒng)所面臨的生態(tài)問題,從整體性、代表性和可操作性3個方面出發(fā),結(jié)合生態(tài)系統(tǒng)的組成、生態(tài)特征、功能作用、經(jīng)濟社會的影響等方面選取了10種評價指標(biāo),通過PSR評價模型構(gòu)建了凌河口濕地生態(tài)系統(tǒng)的健康評價指標(biāo)體系(表2)。

表2 凌河口濕地生態(tài)系統(tǒng)健康評價指標(biāo)體系

2.2.1 壓力指標(biāo)

一般情況下可依據(jù)人類對自然的干擾程度將天然濕地分為4種類型,即完全自然型、受擾自然型、退化自然型和人工修復(fù)型。當(dāng)前,天然濕地生態(tài)系統(tǒng)所面臨的主要壓力是人類頻繁的對天然濕地進(jìn)行干擾所造成的,以前理想中的濕地已逐漸消失殆盡,絕大部分濕地都遭受到了人為因素的破壞[11]。經(jīng)實地勘察,凌河口濕地的生態(tài)健康主要是受到了人類活動的影響,如亂砍亂伐、大面積墾殖、肆意排放污水、破壞生態(tài)環(huán)境等,這些行為造成了濕地景觀的嚴(yán)重破碎化、生物多樣性降低,濕地生態(tài)健康面臨嚴(yán)重威脅。因此,本文選取了人口密度指數(shù)和人類干擾指數(shù)這兩個指標(biāo)來構(gòu)建凌河口濕地生態(tài)系統(tǒng)健康評價的壓力指標(biāo)。

(1) 人口密度

人口密度是單位面積土地上居住的人口數(shù)。它是表示研究區(qū)內(nèi)人口密集程度的指標(biāo)。通常以每平方千米或每公頃內(nèi)的常住人口為計算單位。

人口密度=研究區(qū)人口數(shù)/濕地總面積；(單位：人/km2)

(2)人類干擾指數(shù)

人類干擾是指由于人類生產(chǎn)、生活和其他社會活動形成的干擾體對環(huán)境和生態(tài)系統(tǒng)施加的各種影響。

人類干擾指數(shù)=(旱地面積+城鎮(zhèn)居民用地面積)/濕地總面積×100%

2.2.2 狀態(tài)指標(biāo)

狀態(tài)指標(biāo)不僅反映了環(huán)境要素的變化還體現(xiàn)了環(huán)境政策的終極目標(biāo),其中主要包括自然環(huán)境的現(xiàn)狀和生態(tài)系統(tǒng)狀況等,指標(biāo)選擇主要考慮環(huán)境或生態(tài)系統(tǒng)的生物、物理化學(xué)特征及生態(tài)功能[12- 13]。本文從景觀、濕地功能、生物及物理化學(xué)等多個角度和方面選取了以下六項指標(biāo),六項指標(biāo)值分別是景觀多樣性指數(shù)、平均斑塊面積、水文調(diào)節(jié)指數(shù)、景觀優(yōu)勢度指數(shù)、均勻度指數(shù)和初級生產(chǎn)力。簡述如下：

(1)景觀多樣性指數(shù)(H)

景觀多樣性指數(shù)是指不同類型的景觀要素或生態(tài)系統(tǒng)在空間結(jié)構(gòu)、功能機理和時間動態(tài)等方面的多樣性和變異性,揭示的是濕地中景觀類型的豐富程度、多樣性程度和復(fù)雜程度。H的大小反映了景觀要素的多少和各景觀要素數(shù)量占總要素比例的變化。H=0時說明景觀是均質(zhì)的并且是由單一要素所構(gòu)成；當(dāng)H趨于最大值時說明此時構(gòu)成景觀要素的元素多于兩個并包括兩個,且每種斑塊類型面積占總面積的比例是一致的。計算公式如下：

式中,Pi為第i類景觀類型所占的面積比例；m為景觀類型的數(shù)目。

(2)平均斑塊面積(MPS)

平均斑塊面積就是某類型斑塊的面積和占該類型斑塊總面積的比值。有研究表明可以用MPS來表示景觀的破碎水平,MPS值的變化能夠反饋諸如濕地地類、植物面積、景觀變化等多種景觀生態(tài)信息,是反映景觀性質(zhì)改變的關(guān)鍵性指標(biāo)。計算公式如下：

式中,Si是指第i種濕地景觀類型的總面積,Ni是指第i種濕地景觀類型的斑塊數(shù)目。

(3)水文調(diào)節(jié)指數(shù)

水文調(diào)節(jié)是濕地生態(tài)系統(tǒng)一個重要的服務(wù)功能.濕地具有巨大的滲透能力和蓄水能力,它能儲蓄降水,減少并滯后了降水進(jìn)入江河,在洪水期,它能削減并滯后了洪峰,從而減少了洪水徑流,而在枯水期,它能夠提供水源,以滿足工農(nóng)業(yè)生產(chǎn)的需要。凌河口濕地主要的服務(wù)功能體現(xiàn)在水文調(diào)節(jié)方面,如調(diào)蓄洪水、提供水源,以及眾多水禽鳥類的棲息地和遷徙的中轉(zhuǎn)站。因此,在本文的評價過程中,將濕地的水文調(diào)節(jié)指數(shù)作為評價濕地生態(tài)系統(tǒng)健康評價的一個指標(biāo)。

水文調(diào)節(jié)指數(shù)=(河流面積+灘地面積)/濕地總面積

(4)景觀優(yōu)勢度指數(shù)(D)

景觀優(yōu)勢度指數(shù)指一種或某一類構(gòu)成景觀空間格局中的景觀類型占其支配地位的程度,用最大可能多樣性指數(shù)(Hmax)的離差來表示。D值由小到大依次表示優(yōu)勢度由低到高；同時也表示土地和覆蓋的植被被利用和受支配程度由低到高。優(yōu)勢度指數(shù)與多樣性指數(shù)為負(fù)相關(guān)且存在一定的相關(guān)系數(shù)。計算公式如下：

式中,Hmax為各類型景觀所占面積比例相等的情況下,景觀的最大多樣性指數(shù),Hmax=log2m；Pi為第i類景觀類型所占的面積比例；m為景觀類型的總數(shù)目。

(5)均勻度指數(shù)(SHEI)

均勻度是用來反映濕地景觀類型中每種斑塊在濕地總面積上散布的不均勻程度,用SHEI來表示,其比值由小到大表示各類景觀在濕地中的分布由不均勻到均勻。計算公式如下：

式中,Pi為第i類景觀類型所占的面積比例；m為濕地景觀類型的總數(shù)目。

(6)初級生產(chǎn)力

植被覆蓋指數(shù)(NDVI)是一種指示因子,指植被的生長狀態(tài)和植被覆蓋度的最佳指示因子。研究表明,植被覆蓋指數(shù)(NDVI)與植被的生產(chǎn)能力兩者為正相關(guān),所以衡量研究區(qū)初級生產(chǎn)力的主要指標(biāo)應(yīng)首選濕地內(nèi)年均NDVI值。計算公式如下：

2.2.3 響應(yīng)指標(biāo)

響應(yīng)指標(biāo)是管理部門為使?jié)竦厣鷳B(tài)系統(tǒng)遵循自身繁衍生息的規(guī)律向著健康、良性的方向發(fā)展而采取的必要管理措施。在人類活動的不斷干擾下,濕地生態(tài)系統(tǒng)的內(nèi)部結(jié)構(gòu)和功能將會發(fā)生改變,造成景觀的嚴(yán)重破碎化,進(jìn)而逐步喪失生物多樣性。當(dāng)人類意識到環(huán)境破壞的嚴(yán)重后果時,會采取相應(yīng)的保護(hù)措施。這種一系列的鏈?zhǔn)椒磻?yīng)即濕地生態(tài)系統(tǒng)對人類干擾的響應(yīng)、人類對濕地生態(tài)系統(tǒng)變化的響應(yīng)[14]。因此,響應(yīng)指標(biāo)選取了濕地退化指數(shù)和斑塊破碎化指數(shù),簡述如下：

(1)濕地退化指數(shù)

濕地退化指數(shù)=濕地減少的面積/原始濕地總面積

(2)斑塊破碎化指數(shù)(PD)

濕地景觀總體的斑塊分化程度和破碎化程度用景觀斑塊密度來表示。PD值大,表明一定面積上斑塊規(guī)模小,破碎化程度高；反之,表明景觀類型保存完好。通過對濕地景觀要素的邊緣密度進(jìn)行分析后就可以說明濕地景觀各要素的動態(tài)特性。計算公式如下：

式中,PD表示班塊密度；∑Ni為研究區(qū)景觀斑塊總數(shù)或某景觀要素斑塊類型的斑塊總數(shù)目；∑Ai為研究區(qū)總面積或某景觀斑塊類型的面積。

2.3 單因子評價

邏輯斯蒂增長曲線關(guān)系模型又被叫做自我抑制性曲線方程,是20世紀(jì)20年代Lotka和Volterra在種群生態(tài)學(xué)中的總?cè)簲?shù)量增長過程的研究中提出的,至今應(yīng)用仍比較廣泛。濕地生態(tài)系統(tǒng)健康評價中各個指標(biāo)的測試值并不能線性反應(yīng)出濕地生態(tài)系統(tǒng)健康各個方面和層次狀態(tài)水平。鑒于此,李永健、林茂昌等使用此模型分別對拉魯濕地、閩江河口區(qū)濕地生態(tài)系統(tǒng)健康進(jìn)行了單因子評價：

式中,P表示單項指標(biāo)的生態(tài)安全評價指標(biāo)評價值,R表示單項指標(biāo)測度值,a、b均為常數(shù)。

其中a、b的確定方法為：當(dāng)R=0.01時,P的值近似取0.001；當(dāng)R=0.99時,P的值近似取0.999,則此時方程中的a和b的值求解分別為4.595和9.19,因此,單項指標(biāo)評價模型最終為：

(1)

(2)

其中,當(dāng)單項指標(biāo)量值增加方向與生態(tài)環(huán)境質(zhì)量的增加方向相同時,采用公式(1)進(jìn)行評價；當(dāng)單項指標(biāo)量值的增加方向與生態(tài)環(huán)境質(zhì)量增加方向相反時,采用公式(2)進(jìn)行評價[15- 16]。

根據(jù)凌河口濕地2005年和2014年遙感影像的解譯信息和各個單項指標(biāo)的公式計算,利用單因子評價法得到兩個時期的單因子指標(biāo)值及指標(biāo)測度值(表3)。

表3 2005年和2014年各個單項指標(biāo)值及指標(biāo)測度值

2.4 綜合評價

根據(jù)各個單項指標(biāo)的評價值,采用AHP確定各個指標(biāo)的權(quán)重,再用加權(quán)平均法求得濕地生態(tài)系統(tǒng)健康的綜合評價指數(shù)(CEI),計算公式如下：

式中,Wi為第i個單項指標(biāo)的權(quán)重值,Pi為第i個單項指標(biāo)的評價值。

2.5 評價標(biāo)準(zhǔn)

根據(jù)一般的評價方法,采用連續(xù)的實數(shù)區(qū)間[0,1]來表示濕地生態(tài)系統(tǒng)健康各個等級的標(biāo)準(zhǔn)值。標(biāo)準(zhǔn)值為1時,對應(yīng)最佳狀態(tài)時的濕地生態(tài)系統(tǒng)健康狀況；當(dāng)標(biāo)準(zhǔn)值為0時,對應(yīng)最差時的濕地生態(tài)系統(tǒng)健康狀況。將濕地的生態(tài)系統(tǒng)健康狀況評價等級分為五級：很健康、健康、亞健康、一般病態(tài)、病態(tài),同時結(jié)合大量參考文獻(xiàn)及凌河口濕地的實際情況,確定了凌河口濕地健康評價分類標(biāo)準(zhǔn)(表4)[17- 19]。

表4 濕地健康等級劃分標(biāo)準(zhǔn)

2.6 評價結(jié)果分析

通過綜合評價法計算出凌河口濕地2014年的生態(tài)系統(tǒng)健康指數(shù)為0.405(表5),根據(jù)凌河口濕地生態(tài)系統(tǒng)健康評價等級標(biāo)準(zhǔn)可知,當(dāng)前凌河口濕地生態(tài)系統(tǒng)健康屬于第Ⅲ等級(亞健康)。

從以上評價結(jié)果可以看出,目前凌河口濕地生態(tài)系統(tǒng)健康狀況不容樂觀,由于人類活動的頻繁干擾,濕地內(nèi)部結(jié)構(gòu)遭受到一定程度的破壞,內(nèi)部生態(tài)系統(tǒng)仍保持穩(wěn)定,但是具有極強的敏感性,已有多處生態(tài)異?，F(xiàn)象產(chǎn)生,主要表現(xiàn)在：天然濕地面積大幅度退化、水禽棲息地遭受破壞、水質(zhì)差且富營養(yǎng)化程度嚴(yán)重、景觀多樣性降低、斑塊破碎化嚴(yán)重、生物多樣性降低、生產(chǎn)能力下降、濕地功能水平嚴(yán)重退化、對外界的干擾響應(yīng)迅速等。

2.6.1 壓力分析

目前,凌河口濕地生態(tài)系環(huán)境的惡化主要是受到了人類頻繁的活動干擾,并且對濕地本身產(chǎn)生的壓力也在逐漸增加。經(jīng)計算,2005年研究區(qū)內(nèi)的人口密度為834.57 人/km2,2014年增長到1026.521 人/km2；期間濕地面積減少了近100 km2。由此可知,人口密度變化和濕地面積變化成負(fù)相關(guān),隨著人類干擾壓力的不斷增強,研究區(qū)原有濕地面積日益減少、濕地功能下降、生物多樣性降低。

2.6.2 狀態(tài)分析

從單項指標(biāo)評價結(jié)果的趨勢可以看出,研究區(qū)內(nèi)平均斑塊面積和景觀多樣性顯著降低,平均斑塊面積由2005年的2.908 km2降到2014年的2.569 km2,景觀多樣性由2005年的2.664降低到2014年的2.631。以上結(jié)果足以表明研究區(qū)內(nèi)濕地的景觀異質(zhì)性在逐漸的增加,濕地斑塊破碎化現(xiàn)象也愈來愈嚴(yán)重,濕地功能水平持續(xù)降低,濕地生態(tài)健康面臨病態(tài)的威脅。

表5 2014年凌河口濕地生態(tài)系統(tǒng)健康綜合評價計算結(jié)果

2.6.3 響應(yīng)分析

濕地生態(tài)系統(tǒng)遭受外界環(huán)境的持續(xù)干擾或面臨一定程度的壓力時,就會發(fā)生與之相對應(yīng)的響應(yīng)。研究區(qū)內(nèi)2005年到2014年的天然濕地面積由368.587 km2減少到276.440 km2,減少了92.147 km2,約占原有濕地總面積的25%；景觀破碎度也逐漸增加。結(jié)果表明,隨著研究區(qū)周邊城鎮(zhèn)、農(nóng)村的不斷發(fā)展建設(shè),大量的濕地被侵占、開墾為居民地或建設(shè)用地,濕地資源隨之消耗殆盡,致使水禽棲息地逐漸萎縮,研究區(qū)內(nèi)珍貴水禽的數(shù)量和種類日益減少,生物多樣性逐年降低。而景觀破碎度的增加也說明了凌河口濕地資源遭到了嚴(yán)重的破壞,濕地功能的穩(wěn)定性降低,濕地保護(hù)水平較低。因此,需要控制凌河口濕地生態(tài)系統(tǒng)健康的壓力和生態(tài)威脅,并對其采取相應(yīng)的保護(hù)和修復(fù)措施,逐步改善當(dāng)前的生態(tài)狀況。

3 凌河口濕地生態(tài)系統(tǒng)健康狀況預(yù)警研究

根據(jù)本文第二部分對凌河口濕地2014年生態(tài)系統(tǒng)健康評價的計算方法,分別計算出1995年、2000年、2005年和2009年與之對應(yīng)的綜合評價指數(shù),依次為0.642、0.617、0.524和0.436。結(jié)果表明,凌河口濕地生態(tài)系統(tǒng)健康狀況有由健康狀態(tài)向病態(tài)發(fā)展的趨勢。利用凌河口濕地1995年到2014年的生態(tài)系統(tǒng)健康狀況綜合評價指數(shù),采用灰色系統(tǒng)理論中的系統(tǒng)預(yù)測方法,通過建立GM模型群,構(gòu)建狀態(tài)方程,進(jìn)而實現(xiàn)對凌河口濕地未來20年生態(tài)系統(tǒng)健康狀況的合理預(yù)測。

3.1 灰色預(yù)測GM(1,1)模型構(gòu)建

數(shù)列預(yù)測的基礎(chǔ)是基于累加生成的數(shù)列的GM(1,1)模型,參考鄧聚龍、朱衛(wèi)紅等人的研究[20- 23],將凌河口濕地1995年、2000年、2005年、2009年和2014年的生態(tài)系統(tǒng)健康綜合評價指數(shù)記為原始數(shù)據(jù)序列X(0)(i),(i=1,2,3,4,5)。依據(jù)X(0)(i)建立GM(1,1)模型,并對該預(yù)測模型進(jìn)行檢驗分析。

3.1.1累加算子序列X(1)(i)

由樣本原始數(shù)據(jù)序列X(0)={x(0)(1),x(0)(2),x(0)(3),x(0)(4),x(0)(5)}=(0.6424,0.6172,0.5244,0.4363,0.4051)

對其作1-AGO,得：X(1)={x(1)(1),x(1)(2),x(1)(3),x(1)(4),x(1)(5)}=(0.6424,1.2596,1.7840,2.2203,2.6254)

3.1.2原始序列X(0)(i)準(zhǔn)光滑性檢驗

3.1.3累加算子序列X(1)(i)準(zhǔn)指數(shù)律檢驗

由序列級比公式m(1)(k)∈[1,1.5],δ=0.5,得：m(1)(3)=1.4163,m(1)(4)=1.2446,m(1)(5)=1.1825所以當(dāng)k>3時,m(1)(k)∈[1,1.5],δ=0.5,準(zhǔn)指數(shù)律成立。

綜上可知,可以對序列X(1)(i)建立GM(1,1)模型。

3.1.4構(gòu)建GM(1,1)模型

X(1)(i)的緊鄰均值生成：Z(1)={z(1)(2),z(1)(3),z(1)(4),z(1)(5)}=(0.9510,1.5218,2.0022,2.4229)

構(gòu)造數(shù)據(jù)矩陣B和數(shù)據(jù)向量Y,設(shè)P=(a,b)T,根據(jù)最小二乘法估計參數(shù),有：

由于GM(1,1)模型得到的是一次累加量,需要將預(yù)測值還原為x(0)(k),即

3.2 預(yù)測模型精度檢驗

3.2.1 相對誤差ε檢驗

表6 預(yù)測模型殘差檢驗

3.2.2灰色絕對關(guān)聯(lián)度O檢驗

所以,灰色絕對關(guān)聯(lián)度

3.2.3 小誤差概率P檢驗

綜上,小誤差概率P=1>0.95；灰色絕對關(guān)聯(lián)度O=0.9978>0.90；后驗差比值C=S2/S1=0.1107<0.35,根據(jù)精度檢驗等級參照表(表7)可知：3個檢驗指標(biāo)的精度檢驗等級均為一級,預(yù)測模型的預(yù)測結(jié)果好,說明本文所建立的GM(1,1)模型是可行的,因此可以對未來凌河口濕地生態(tài)系統(tǒng)健康狀況進(jìn)行預(yù)測分析。

表7 精度檢驗等級參照表

3.3 預(yù)測結(jié)果分析

圖2 凌河口濕地生態(tài)系統(tǒng)健康綜合評價指數(shù)模擬與預(yù)測 Fig.2 The predictive value of Linghekou Wetland ecological health competitive evaluation

4 結(jié)論

本研究依據(jù)凌河口濕地生態(tài)系統(tǒng)實際面臨的生態(tài)問題,采用壓力-狀態(tài)-響應(yīng)模型(PSR)構(gòu)建了凌河口濕地生態(tài)系統(tǒng)健康評價指標(biāo)體系,利用AHP方法確定指標(biāo)體系中的各個單項指標(biāo)的權(quán)重,然后運用邏輯斯蒂增長曲線模型對其進(jìn)行單因子指標(biāo)評價,最后通過測算CEI值的綜合評價法對1995年到2014年凌河口濕地生態(tài)系統(tǒng)健康進(jìn)行綜合評價,由此得出結(jié)論：1995年到2000年濕地生態(tài)系統(tǒng)健康狀況為健康；2005年到2014年濕地生態(tài)系統(tǒng)健康狀況為亞健康,并有向一般病態(tài)發(fā)展的趨勢。

[1] KEDDY P A. Wetland Ecology: Principles and Conservation. Cambridge: Cambridge University Press, 2000: 124- 238.

[2] Holland M M. Wetlands and environmental gradients//Mulamoottil G, Warnser B G, Mcbean E A, eds. Wetlands Environmental Gradients, Boundaries, and Buffers. Boca Raton: CRC Press Inc., 1996: 112- 131.

[3] 林佳, 宋戈, 宋思銘. 景觀結(jié)構(gòu)動態(tài)變化及其土地利用生態(tài)安全——以建三江墾區(qū)為例. 生態(tài)學(xué)報, 2011, 31(20): 5918- 5927.

[4] 崔保山, 楊志峰. 濕地生態(tài)系統(tǒng)健康評價指標(biāo)體系I理論. 生態(tài)學(xué)報, 2002, 22(7): 1005- 1011.

[5] 肖篤寧, 陳文波, 郭福良. 論生態(tài)安全的基本概念和研究內(nèi)容. 應(yīng)用生態(tài)學(xué)報, 2002, 13(3): 354- 358.

[6] 李新香. 凌河口濕地科學(xué)考察. 遼寧林業(yè)科技, 2004, (1): 34- 35.

[7] 周林飛, 徐浩田, 張靜. 凌河口濕地自然保護(hù)區(qū)景觀格局變化及功能區(qū)劃分. 濕地科學(xué), 2016, 14(3): 403- 47.

[8] 麥少芝, 徐頌軍, 潘穎君. PSR模型在濕地生態(tài)系統(tǒng)健康評價中的應(yīng)用. 熱帶地理, 2005, 25(4): 317- 321.

[9] 仝川. 環(huán)境指標(biāo)研究進(jìn)展與分析. 環(huán)境科學(xué)研究, 2000, 13(4): 53- 55.

[10] 張崢, 張建文, 李寅年,王欣,劉泓,張敬國. 濕地生態(tài)評價指標(biāo)體系. 農(nóng)業(yè)環(huán)境保護(hù), 1999, 18(6): 283- 285.

[11] 崔保山, 楊志峰. 濕地生態(tài)系統(tǒng)健康研究進(jìn)展.生態(tài)學(xué)雜志, 2001, 20(3): 31- 36.

[12] 王一涵, 周德民, 孫永華. RS和GIS支持的洪河地區(qū)濕地生態(tài)健康評價. 生態(tài)學(xué)報, 2011, 31(13): 3590- 3600.

[13] 周林飛, 許士國, 孫萬光. 基于壓力-狀態(tài)-響應(yīng)模型的扎龍濕地健康水循環(huán)評價研究. 水科學(xué)進(jìn)展, 2008, 19(2): 205- 213.

[14] 朱衛(wèi)紅, 郭艷麗, 孫鵬, 苗承玉, 曹光蘭. 圖們江下游濕地生態(tài)系統(tǒng)健康評價. 生態(tài)學(xué)報, 2012, 32(21): 6609- 6618.

[15] 李永健. 拉魯濕地生態(tài)環(huán)境質(zhì)量評價的景觀生態(tài)學(xué)方法應(yīng)用研究[D]. 成都: 四川大學(xué), 2002.

[16] 林茂昌. 基于RS和GIS的閩江河口區(qū)濕地生態(tài)環(huán)境質(zhì)量評價[D]. 福州: 福建師范大學(xué), 2005.

[17] 高桂芹. 東平湖濕地生態(tài)系統(tǒng)健康評價研究[D]. 濟南: 山東師范大學(xué), 2006.

[18] 王薇. 黃河三角洲濕地生態(tài)系統(tǒng)健康綜合評價研究——以墾利縣為例[D]. 泰安: 山東農(nóng)業(yè)大學(xué), 2007.

[19] 鄭耀輝. 濱海紅樹林濕地生態(tài)系統(tǒng)健康評價——以珠江口淇澳島紅樹林為例[D]. 廣州: 中山大學(xué), 2009.

[20] 劉思峰, 鄧聚龍. GM(1,1)模型的適用范圍. 系統(tǒng)工程理論與實踐, 2000, 20(5): 121- 124.

[21] 崔立志. 灰色預(yù)測技術(shù)及其應(yīng)用研究[D]. 南京: 南京航空航天大學(xué), 2010.

[22] 賈得海, 曾建初. 灰色系統(tǒng)理論GM(1,1)預(yù)測模型的應(yīng)用. 昆明理工大學(xué)學(xué)報: 自然科學(xué)版, 2013, 38(6): 115- 120.

[23] 朱衛(wèi)紅, 苗成玉, 鄭小軍, 曹光蘭, 王凡凡. 基于3S技術(shù)的圖們江流域濕地生態(tài)安全評價與預(yù)警研究. 生態(tài)學(xué)報, 2014, 34(6): 1379- 1390.