劉 瑩，王立權(quán)，孔 達(dá)，王子龍

(1.黑龍江大學(xué)水利電力學(xué)院，哈爾濱 150080；2.東北農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院，哈爾濱 150030)

0 引 言

生態(tài)系統(tǒng)健康是指在生態(tài)系統(tǒng)基本單元、區(qū)域或景觀以及全球等多個(gè)尺度上生態(tài)系統(tǒng)的自我保持與更新的能力和社會(huì)經(jīng)濟(jì)服務(wù)功能[1]。隨著社會(huì)經(jīng)濟(jì)的快速發(fā)展，人類高強(qiáng)度活動(dòng)導(dǎo)致的生存環(huán)境惡化、水土資源萎縮、生態(tài)功能退化、承載能力降低、生物多樣性縮減等生態(tài)系統(tǒng)問題日益突出，由此給社會(huì)經(jīng)濟(jì)可持續(xù)性發(fā)展帶來了嚴(yán)重影響。生態(tài)系統(tǒng)健康日益受到國際社會(huì)的廣泛關(guān)注，成為了生態(tài)系統(tǒng)管理的重要問題和熱點(diǎn)領(lǐng)域[2]。生態(tài)系統(tǒng)健康評(píng)價(jià)興起于20世紀(jì)80年代，在區(qū)域或景觀尺度上已經(jīng)開展了河流[3]、湖泊[4]、森林[5]、濕地[6]、流域[7]、海島[8]、城市[9]等多種類型區(qū)域的生態(tài)系統(tǒng)健康評(píng)價(jià)，由簡單到復(fù)雜、由定性到定量、由關(guān)注生物指標(biāo)到重視服務(wù)功能，已成為該領(lǐng)域研究的重要趨勢(shì)。

三江平原是三江(松花江、黑龍江和烏蘇里江)匯流、沖積形成的沼澤化低平原，地理位置位于黑龍江省東部區(qū)域，行政區(qū)劃包括佳木斯市、鶴崗市、雙鴨山市、七臺(tái)河市、雞西市、穆棱市、依蘭縣，同時(shí)包括分布其中的幾十個(gè)國有農(nóng)場(chǎng)和林業(yè)局[10]。該區(qū)域經(jīng)過幾十年的開發(fā)，生態(tài)系統(tǒng)受到了一定的影響，特別是以地下水過度利用為代表的水循環(huán)變化引起的生態(tài)環(huán)境效應(yīng)凸顯，有必要定量評(píng)估三江平原的生態(tài)系統(tǒng)健康狀況。

1 評(píng)價(jià)指標(biāo)體系的構(gòu)建

評(píng)價(jià)指標(biāo)體系的構(gòu)建是生態(tài)系統(tǒng)健康綜合評(píng)價(jià)的基礎(chǔ)和核心，不同研究中提出了多種生態(tài)系統(tǒng)健康綜合評(píng)價(jià)指標(biāo)體系的框架或方案[11-15]。針對(duì)三江平原生態(tài)系統(tǒng)特征，借鑒其他學(xué)者相關(guān)研究成果[16-18]，依據(jù)評(píng)價(jià)指標(biāo)體系建立的原則，從社會(huì)、經(jīng)濟(jì)、生態(tài)3個(gè)子系統(tǒng)相互依存和相互作用的關(guān)系入手，注重生態(tài)系統(tǒng)的服務(wù)功能和價(jià)值，構(gòu)建了三江平原生態(tài)系統(tǒng)健康評(píng)價(jià)指標(biāo)體系包含目標(biāo)層、系統(tǒng)層和指標(biāo)層三個(gè)層次，如圖1所示。

第一層為目標(biāo)層(三江平原生態(tài)系統(tǒng)健康A(chǔ))，為本研究所追求的總目標(biāo)；第二層為系統(tǒng)層，是與此研究目標(biāo)緊密聯(lián)系的三個(gè)系統(tǒng)，即生態(tài)系統(tǒng)B1、社會(huì)系統(tǒng)B2、經(jīng)濟(jì)系統(tǒng)B3；第三層為目標(biāo)層，共包括12項(xiàng)指標(biāo)，各指標(biāo)的計(jì)算公式及指標(biāo)表征含義見表1。

圖1 三江平原生態(tài)系統(tǒng)健康評(píng)價(jià)指標(biāo)體系

表1 三江平原生態(tài)系統(tǒng)健康評(píng)價(jià)指標(biāo)計(jì)算公式及其含義

2 生態(tài)系統(tǒng)健康評(píng)價(jià)方法

生態(tài)系統(tǒng)健康評(píng)價(jià)的模型方法所采用的指標(biāo)體系法，通過指標(biāo)層次和指標(biāo)選取體現(xiàn)待評(píng)價(jià)生態(tài)系統(tǒng)的特點(diǎn)與價(jià)值，能夠更為全面地體現(xiàn)復(fù)雜生態(tài)系統(tǒng)健康狀況[19]，本研究采用指標(biāo)體系法，構(gòu)建生態(tài)系統(tǒng)健康評(píng)價(jià)的指標(biāo)體系，運(yùn)用數(shù)學(xué)模型對(duì)三江平原生態(tài)系統(tǒng)健康進(jìn)行綜合評(píng)價(jià)。

2.1 粒子群優(yōu)化算法

粒子群優(yōu)化(Particle Swarm Optimization，簡稱PSO)算法是一種源于對(duì)鳥群捕食行為的研究而發(fā)明的進(jìn)化計(jì)算技術(shù)[20]。粒子群體中的每個(gè)粒子代表解空間中一個(gè)可能的解，群體中每個(gè)粒子在優(yōu)化過程中所找到的最好位置就是該粒子本身所找到的最好解，整個(gè)粒子群體所找到的最好位置就是整個(gè)粒子群體目前找到的最好解。前者稱為個(gè)體極值，后者稱為全局極值，每個(gè)粒子都通過這兩個(gè)極值不斷更新自己位置和速度，從而產(chǎn)生新一代粒子群體，在整個(gè)過程中粒子群體對(duì)解空間進(jìn)行全面搜索尋找最優(yōu)[21]。

2.2 基于粒子群優(yōu)化的投影尋蹤模型

運(yùn)用投影尋蹤模型開展評(píng)價(jià)問題研究的基本步驟是：根據(jù)數(shù)據(jù)結(jié)構(gòu)尋找最佳投影方向，運(yùn)用投影值計(jì)算公式，計(jì)算被評(píng)價(jià)對(duì)象各樣本的最佳投影值，將各樣本的最佳投影值進(jìn)行相對(duì)比較，以確定待評(píng)價(jià)樣本的優(yōu)劣。對(duì)于投影尋蹤模型而言，如何獲得最佳投影方向是評(píng)價(jià)模型構(gòu)建的核心問題，最佳投影方向的尋求可以通過一個(gè)復(fù)雜的非線性優(yōu)化問題來解決[22]。本研究采用操作相對(duì)簡單的粒子群算法優(yōu)化最佳投影方向，建立投影尋蹤評(píng)價(jià)模型[23]。

3 結(jié)果與分析

3.1 三江平原不同時(shí)期生態(tài)系統(tǒng)健康變化分析

采用粒子群優(yōu)化算法對(duì)投影尋蹤評(píng)價(jià)模型中的投影方向進(jìn)行優(yōu)化。利用MATLAB軟件編制基于粒子群優(yōu)化的投影尋蹤評(píng)價(jià)模型程序，其中粒子群優(yōu)化算法的參數(shù)設(shè)置如下：種群規(guī)模N=400，最大迭代次數(shù)Gmax=30，慣性權(quán)重w=0.996 5，學(xué)習(xí)因子c1=c2=2。在MATLAB環(huán)境下運(yùn)行程序，得到投影指標(biāo)函數(shù)最大值Q*(a)、最佳投影方向a*以及粒子群優(yōu)化算法的迭代尋優(yōu)過程，如圖2所示。

圖2 粒子群優(yōu)化算法的尋優(yōu)過程(三江平原)

采用2000-2011期間12年的數(shù)據(jù)，在MATLAB環(huán)境下運(yùn)行程序，得出投影指標(biāo)函數(shù)最大值Q*(a)=0.067 205 862 065 981，最佳投影方向a*=(0.727 930 300 198 792，0.001 433 060 959 954，0.287 577 467 128 917，0.622 426 401 101 767)。將最佳投影方向 代入公式即可得評(píng)價(jià)指標(biāo)的最佳值a*=(1.431 847 331 635 9，1.600 378 867 016 74，1.594 041 503 740 58，1.393 663 791 829 91，1.393 663 537 376 41，1.32711440223315，1.393 667 335 629 2，1.254 983 159 846 58，1.155 477 796 397 66，1.263 885 522 669 75，1.091 869 707 395 63，1.011 533 228 609 15)。

計(jì)算可得三江平原地區(qū)2000-2011年12年的最佳投影值，如圖3所示。由圖3可知，2000-2011年，三江平原生態(tài)健康指數(shù)處于遞減趨勢(shì)，也就是說該期間待評(píng)價(jià)地區(qū)的生態(tài)健康情況處于下滑階段。

圖3 2000-2011年三江平原生態(tài)系統(tǒng)健康指數(shù)

分析三江平原生態(tài)系統(tǒng)健康評(píng)價(jià)的12年結(jié)果可知，地下水資源模數(shù)由2000年的8.81 萬m3/(km2·a)逐步下降至2011年的5.00 萬m3/(km2·a)，而地下水供水比趨于穩(wěn)定，所以地下水資源的壓力越來越大，對(duì)三江平原生態(tài)系統(tǒng)健康有逆向影響。而此階段大力發(fā)展工業(yè)及相關(guān)行業(yè)，使人均GDP得到了迅猛發(fā)展，由2000年的6 440.97 元/人增加到2011年的30 247.17 元/人，說明經(jīng)濟(jì)發(fā)展水平、人民的生活質(zhì)量得到了提高，對(duì)水利的投入提供了保障。但同時(shí)，工業(yè)的發(fā)展也對(duì)三江平原生態(tài)系統(tǒng)的生態(tài)功能、物種的多樣性、水質(zhì)的標(biāo)準(zhǔn)起到了惡化的效果。

3.2 三江平原不同地區(qū)生態(tài)系統(tǒng)健康對(duì)比分析

根據(jù)前述評(píng)價(jià)模型的理論和方法，分別對(duì)2006年和2011年三江平原各地區(qū)生態(tài)健康狀況進(jìn)行評(píng)價(jià)。運(yùn)行程序語言，計(jì)算投影指標(biāo)函數(shù)最大值為Q*(a)=0.252 847 734 301 184，最佳投影方向a*=(0.002 922 068 891 988，0.485 527 633 730 488，0.840 718 624 038 648，0.016 525 437 155 085，0.000 381 150 453 593，0.000 074 440 839 243，0.136 551 761 430 636，0.196 282 829 488 327)，并得到粒子群優(yōu)化算法的迭代尋優(yōu)過程，如圖4所示。

圖4 粒子群優(yōu)化算法的尋優(yōu)過程(三江平原不同地區(qū))

依據(jù)粒子群優(yōu)化算法獲得的最佳投影方向計(jì)算可得2006年和2011年評(píng)價(jià)指標(biāo)數(shù)據(jù)集的最佳投影值，如圖5、圖6所示。

圖5 2006年三江平原不同地區(qū)生態(tài)系統(tǒng)健康狀況

圖6 2011年三江平原不同地區(qū)生態(tài)系統(tǒng)健康狀況

由圖5可知，2006年穆棱市的生態(tài)系統(tǒng)健康狀況處于三江平原平均水平之上，依蘭縣次之，佳木斯市、雞西市、鶴崗市、雙鴨山、七臺(tái)河市六個(gè)地區(qū)生態(tài)健康狀況處于三江平原平均水平，其中，穆棱市生態(tài)系統(tǒng)健康狀況最優(yōu)，雙鴨山市生態(tài)系統(tǒng)健康狀況最差。

對(duì)比分析圖5與圖6，自2006-2011年，三江平原區(qū)域內(nèi)的依蘭縣、穆棱市、佳木斯市、雞西市、鶴崗市、雙鴨山市、七臺(tái)河市生態(tài)系統(tǒng)健康狀況均處于下滑狀態(tài)，即三江平原生態(tài)健康狀況持續(xù)惡化。對(duì)比分析各指標(biāo)數(shù)據(jù)，相對(duì)于2006年，雖然2011年三江平原各地區(qū)的人均GDP等各項(xiàng)正向指標(biāo)有所上升，但是2011年三江平原各地區(qū)的地下水供水比、人類活動(dòng)強(qiáng)度、人口健康狀況等逆向指標(biāo)上升幅度更為明顯，表明人類活動(dòng)對(duì)生態(tài)系統(tǒng)健康的影響很大。

4 結(jié) 語

隨著三江平原各地區(qū)經(jīng)濟(jì)的快速發(fā)展和人口數(shù)量的增多，水資源需求量會(huì)日益增加，對(duì)生態(tài)環(huán)境的污染 日益加劇，濕地面積不斷減小。因此，應(yīng)盡快提高水資源利用率，在加快經(jīng)濟(jì)的發(fā)展的同時(shí)，力爭做到社會(huì)-經(jīng)濟(jì)-生態(tài)的和諧與可持續(xù)發(fā)展是保持三江平原生態(tài)系統(tǒng)健康發(fā)展的必要措施。

[1] Rapport D J, Costanza R, McMichael A J. Assessing ecosystem health [J]. Trends in Ecology & Evolution, 1998,13(10):397-402.

[2] Wang R S. Integrative eco-management for resource, environment and industrial transformation [J]. Systems Engineering-theory & Practice, 2003,(2):125-132.

[3] 朱衛(wèi)紅, 曹光蘭, 李 瑩, 等. 圖們江流域河流生態(tài)系統(tǒng)健康評(píng)價(jià)[J]. 生態(tài)學(xué)報(bào), 2014,34(14):3 969-3 977.

[4] 張艷會(huì), 楊桂山, 萬榮榮. 湖泊水生態(tài)系統(tǒng)健康評(píng)價(jià)指標(biāo)研究[J]. 資源科學(xué), 2014,36(6):1 306-1 315.

[5] 袁 菲, 張星耀, 梁 軍. 基于干擾的汪清林區(qū)森林生態(tài)系統(tǒng)健康評(píng)價(jià)[J]. 生態(tài)學(xué)報(bào), 2013,33(12):3 722-3 731.

[6] 郝文淵, 李文博, 王忠斌, 等. 西藏拉薩河谷拉魯濕地生態(tài)系統(tǒng)健康評(píng)價(jià)[J]. 干旱區(qū)資源與環(huán)境, 2013,27(5):95-99.

[7] 廖靜秋, 曹曉峰, 汪 杰, 等. 基于化學(xué)與生物復(fù)合指標(biāo)的流域水生態(tài)系統(tǒng)健康評(píng)價(jià)----以滇池為例[J]. 環(huán)境科學(xué)學(xué)報(bào), 2014,34(7):1 845-1 852.

[8] 陸麗珍, 詹遠(yuǎn)增, 葉艷妹, 等. 基于土地利用空間格局的區(qū)域生態(tài)系統(tǒng)健康評(píng)價(jià)----以舟山島為例[J]. 生態(tài)學(xué)報(bào), 2010,30(1):245-252.

[9] 趙 帥, 柴立和, 李鵬飛, 等. 城市生態(tài)系統(tǒng)健康評(píng)價(jià)新模型及應(yīng)用----以天津市為例[J]. 環(huán)境科學(xué)學(xué)報(bào), 2013,33(4):1 173-1 179.

[10] 何 璉. 中國三江平原[M]. 哈爾濱: 黑龍江科學(xué)技術(shù)出版社, 2000.

[11] 劉明華, 董貴華. RS和GIS支持下的秦皇島地區(qū)生態(tài)系統(tǒng)健康評(píng)價(jià)[J]. 地理研究, 2006,25(5):930-938.

[12] 崔保山, 楊志峰. 濕地生態(tài)系統(tǒng)健康評(píng)價(jià)指標(biāo)體系Ⅱ.方法與案例[J]. 生態(tài)學(xué)報(bào), 2002,22(8):1 231-1 239.

[13] 陳 銘, 張樹清, 王志強(qiáng), 等. 基于GIS的蛟流河流域濕地生態(tài)系統(tǒng)健康評(píng)價(jià)[J]. 山東農(nóng)業(yè)大學(xué)學(xué)報(bào)自然科學(xué)版, 2007,22(1):165-168.

[14] 戴全厚, 劉國彬, 田均良, 等. 侵蝕環(huán)境小流域生態(tài)經(jīng)濟(jì)系統(tǒng)健康定量評(píng)價(jià)[J]. 生態(tài)學(xué)報(bào), 2006,26(7):2 219-2 228.

[15] 官冬杰, 蘇維詞. 城市生態(tài)系統(tǒng)健康評(píng)價(jià)方法及其應(yīng)用研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 2006,26(10):1 716-1 722.

[16] 王慧亮, 王學(xué)雷, 莫明浩, 等. 基于生態(tài)健康的洪湖濕地恢復(fù)評(píng)價(jià)[J]. 武漢大學(xué)學(xué)報(bào)(理學(xué)版), 2010,56(5):557-563.

[17] 王文杰, 張 哲, 王 維, 等. 流域生態(tài)健康評(píng)價(jià)框架及其評(píng)價(jià)方法體系研究(一)----框架和指標(biāo)體系[J]. 環(huán)境工程技術(shù)學(xué)報(bào), 2012,2(4):271-277.

[18] 李春暉, 鄭小康, 崔 嵬, 等. 衡水湖流域生態(tài)系統(tǒng)健康評(píng)價(jià)[J]. 地理研究, 2008,27(3):565-573.

[19] 馬克明, 孔紅梅, 關(guān)文彬, 等. 生態(tài)系統(tǒng)健康評(píng)價(jià):方法與方向[J]. 生態(tài)學(xué)報(bào), 2001,21(12):2 106-2 116.

[20] Kennedy J, Eberhart R. Particle swarm optimization[C]∥ IEEE International Conference on Neural Networks, 1995. Proceedings. IEEE, 1995:1942-1948 vol.4.

[21] 雷崇民, 任再敏, 高岳林. 約束優(yōu)化問題帶有混沌變異的PSO-DE混合算法[J]. 蘭州理工大學(xué)學(xué)報(bào), 2010,36(6):94-98.

[22] 陳廣洲, 汪家權(quán). 基于模擬退火算法的投影尋蹤方向優(yōu)化[J]. 合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2008,31(8):1 315-1 317.

[23] 蔡怡馨, 王 靜, 劉鵬玲,等. 利用組合模型評(píng)價(jià)云南省水資源利用效率[J]. 人民黃河, 2015,37(5):58-61.