馬冰然，馬安青，陳 爽，李正炎，王志凱

（中國海洋大學環(huán)境科學與工程學院，山東 青島266100）

水生態(tài)系統(tǒng)對地球環(huán)境具有重要的調節(jié)功能，為生物提供棲息環(huán)境，同時與人類生存、發(fā)展關系密切。目前，流域水生態(tài)系統(tǒng)主要以各種分區(qū)為單元進行管理，并且將自然資源按照生態(tài)系統(tǒng)功能進行分區(qū)規(guī)劃，制定保護和發(fā)展措施，已經(jīng)在科學上被廣泛接受。1967年，生態(tài)區(qū)的概念第一次出現(xiàn)［1］，這意味著研究從傳統(tǒng)的地理學分區(qū)邁入了生態(tài)學領域。諸多學者都曾在此領域開展過深入研究［2－4］，目前水生態(tài)系統(tǒng)分區(qū)體系大體可以分為三種類型，第一類型是基于景觀要素的分區(qū)體系，如美國水生態(tài)分區(qū)［4］；第二類型是基于水生生物要素的分區(qū)體系，如世界淡水生態(tài)分區(qū)［5］；第三類型是基于價值判斷和功能的分區(qū)體系，如中國水功能區(qū)劃。美國生態(tài)學家Bailey于1976年提出了美國生態(tài)區(qū)域的等級概念，并繪制出一套美國生態(tài)區(qū)域圖［2］。McMahon等提出了美國范圍內生態(tài)區(qū)劃的通用方案［6］。澳大利亞的學者用氣候（直接表現(xiàn)為降水量大小及其季節(jié)性）、地貌以及植被類型3個景觀要素指標來反映水生態(tài)系統(tǒng)的自然差異性，從而將澳洲分成了多種生態(tài)區(qū)［7］。2000年，歐盟委員會提出了歐洲水框架指令，第一次全方面闡述了歐洲淡水生態(tài)系統(tǒng)的管理理念，并對不同水域給予了評價［8］。

在我國，各省根據(jù)生態(tài)系統(tǒng)特征、生態(tài)服務功能的重要性、生態(tài)敏感性以及區(qū)域面臨的生態(tài)環(huán)境問題，將某一個特定的研究區(qū)域劃分為自然生態(tài)區(qū)、生態(tài)亞區(qū)和生態(tài)功能區(qū)3個等級單元，并且該項工作已經(jīng)初步完成。很多城市也對生態(tài)功能區(qū)劃進行了研究。在此基礎上，研究尺度開始拓展到流域層面，如趙銳峰等對塔里木河中下游河流廊道的分區(qū)研究，劃分了3大景觀生態(tài)功能類型區(qū)［9］。燕乃玲等對我國長江源區(qū)以流域生態(tài)系統(tǒng)單元為基礎進行了生態(tài)功能區(qū)劃，共劃分了5個生態(tài)功能分區(qū)，為長江源區(qū)生態(tài)系統(tǒng)管理和保護提供了基礎框架［10］。但是，目前我國大部分重點流域還沒有進行過此方面的研究，因此本文針對大遼河流域進行的水生態(tài)功能區(qū)劃研究具有重要意義。

1 研究區(qū)域概況及環(huán)境現(xiàn)狀

研究區(qū)域總體位于：40°40′N～41°10′N，121°50′E～122°40′E。渾河和太子河在海城三岔河匯合，之后稱為大遼河，自營口流入渤海遼東灣。大遼河流經(jīng)海城、盤山、大石橋、大洼、營口等縣市，穿行于遼河中下游區(qū)的近海地帶，沿岸屬于濱海與堆積平原。大遼河沿岸工農(nóng)業(yè)發(fā)達，城市的工業(yè)污水、生活污水和農(nóng)業(yè)的面源污染嚴重污染了河流和流域的土壤。作為河口型河流，2006年大遼河排入遼東灣污水1.35×109t［11］，渤海水域污染令人堪憂。近幾年來，雖然大遼河污染受到國內外學者的重視并有轉好趨勢，但環(huán)境形勢仍十分嚴峻，需要政府和有關部門進行長期的規(guī)劃與治理。

2 技術路線

本文根據(jù)流域濕地景觀變化分析、生態(tài)環(huán)境敏感性評價、水生態(tài)系統(tǒng)服務功能評價，確定功能分區(qū)和分區(qū)依據(jù)，在此基礎上，調查收集數(shù)據(jù)，識別主導水體類型、人類干擾特征，并結合水生態(tài)系統(tǒng)功能評估內容，進行大遼河流域的水生態(tài)系統(tǒng)功能區(qū)劃。主要技術路線如圖1所示。

圖1 大遼河流域水生態(tài)功能區(qū)劃研究技術路線圖Fig.1 The technology roadmap on function regionalization of water ecosystem for the liaohe drainage area

3 評價分析

3.1 大遼河流域濕地景觀格局變化分析

景觀空間格局是指大小和形狀不同的景觀斑塊在空間上的排列，景觀格局隨時間變化可以反映景觀生態(tài)過程［12］。大遼河流域濕地景觀類型面積隨時間的變化特征主要表現(xiàn)在：庫塘面積變化最顯著，擴大的比例最大；赤堿蓬面積縮小最為明顯；蘆葦田面積變化最不穩(wěn)定；水稻田、養(yǎng)殖區(qū)、建筑用地和灘涂的面積則呈現(xiàn)出逐年增加的趨勢。濕地景觀類型的空間變化特征主要表現(xiàn)在自然濕地大量轉化為人工景觀或人工濕地，具體表現(xiàn)在轉化為建筑用地、水稻田、養(yǎng)殖區(qū)。斑塊的破碎度增大，分布分散，景觀的破碎化程度增強，說明人類活動對濕地的干擾強度增強［13］。

3.2 大遼河流域生態(tài)環(huán)境敏感性評價

生態(tài)環(huán)境敏感性是指生態(tài)系統(tǒng)對自然環(huán)境變化和人類活動干擾的反映程度，在GIS的支持下，運用層次分析法（AHP）［14－15］和模糊綜合評判［16］相結合的方法開發(fā)出GIS系統(tǒng)，得到大遼河流域生態(tài)環(huán)境敏感性評價結果，發(fā)現(xiàn)整個地區(qū)敏感等級為中度敏感和高度敏感。其中，高度敏感是主要的敏感類型，占總面積的69.62%，在全區(qū)除盤山縣外，5個縣市均有分布，其中盤錦市的東部、西南少數(shù)地區(qū)以及大石橋市的中部和南部少數(shù)地區(qū)的敏感度值最高；中度敏感的地區(qū)占總地區(qū)土地面積的份額相對較小，占總土地面積的30.38%，主要分布在盤山縣區(qū)，蓋州東南部分地區(qū)＊。

＊陳爽，王其翔，馬安青，等.基于GIS的大遼河口生態(tài)環(huán)境敏感性評價［J］.中國海洋大學學報：自然科學版（in press）

3.3 大遼河流域生態(tài)系統(tǒng)服務功能重要性評價

生態(tài)系統(tǒng)服務是指生態(tài)系統(tǒng)與生態(tài)過程所形成和維持的人類賴以生存的自然環(huán)境條件與效用。通過指標體系評價［17］的方法，其評價結果為：重要性較高的區(qū)域主要分布在大遼河上游的牛莊鎮(zhèn)地區(qū)，中游的大洼縣、溝沿鎮(zhèn)、水源鎮(zhèn)、老邊區(qū)和大石橋市，這些地區(qū)的城鎮(zhèn)建筑用地較多，養(yǎng)殖業(yè)發(fā)展較好，農(nóng)業(yè)發(fā)展程度比較高；功能等級最高的區(qū)域分布在大遼河的上游和下游入?？诘貐^(qū)，這些地區(qū)分布有蘆葦、灘涂和沼澤，入?？诘貐^(qū)更是分布有大量的養(yǎng)殖區(qū)并且靠近城市建筑區(qū)，多樣性維持功能、水環(huán)境及水文支持功能等都比較高。重要性較低的區(qū)域主要分布在大遼河上游的西北部沙嶺、壩墻子等地以及盤錦的西部和營口東部的感王鎮(zhèn)、虎莊鎮(zhèn)等地，這些地區(qū)旱地較多，土地利用類型較為簡單；功能等級為一般的區(qū)域主要分布在大遼河上游盤錦北部的新開鎮(zhèn)等地、中游地段的大洼縣的西部地區(qū)、營口市老邊區(qū)的北部以及營口市東部西柳鎮(zhèn)等。將水生生物多樣性維持功能、野生生物棲息地維持功能、濱岸帶支持功能、水環(huán)境支持功能和水文支持功能各個單項功能的指標應用SPSS軟件進行了主成分分析，得到前3個單項功能指標的解釋方差累積貢獻率大于85%，對整個大遼河流域水生態(tài)系統(tǒng)功能重要性影響最大，使整個研究區(qū)域的綜合評價結果出現(xiàn)了較為明顯的差異［18］。

4 大遼河流域水生態(tài)系統(tǒng)功能區(qū)劃方法和技術

4.1 研究方法概述

生態(tài)功能區(qū)劃的最終結果是生態(tài)系統(tǒng)單元圖，由于區(qū)劃目的不同，采用的區(qū)劃方法也不盡相同，主要有以下幾種［19］：地理相關法、空間疊置法、景觀制圖法、定量分析法等。大遼河流域水生態(tài)功能區(qū)劃的研究采取了定性分區(qū)與定量分區(qū)相結合的方法，利用計算機圖形空間疊置法、相關分析法和專家集成的方法，依照生態(tài)功能分區(qū)的等級體系，運用自下而上的劃分方法進行分區(qū)劃界。自下而上的劃分方法是將具有相似性水生態(tài)系統(tǒng)功能的小流域聚合成為大的分區(qū)單元，通過水生生物、生境類型、水化學等河道內指標來識別水生態(tài)功能區(qū)。本文研究根據(jù)區(qū)域環(huán)境要素和水生態(tài)系統(tǒng)異質性特征，選取基于生態(tài)系統(tǒng)結構的分類指標，并根據(jù)各指標反映的水生態(tài)功能信息，剔除一些信息量重復或相關度較大的指標，將原始數(shù)據(jù)標準化處理后，在GIS支持下分別提取和計算各指標的指數(shù)值，最后通過聚類分析獲得區(qū)劃結果。

4.2 水體類型識別

對小流域內的主要水體類型進行識別主要從3方面進行判斷：鹽度值、河流彎曲度和河道寬度。鹽度值 根據(jù)鹽度值在樣區(qū)的分布特點以及咸淡水鹽度值的分異規(guī)律，表現(xiàn)為從三岔河向河口，鹽度呈遞增趨勢，在河口地區(qū)已經(jīng)接近為海水。

河流彎曲度 即河流的分維數(shù)＝河流總長度／河流的流域面積，具體計算步驟為：在ArcGIS9.3中將河流線狀圖層與整個研究區(qū)格網(wǎng)面狀矢量圖進行疊加，在矢量格網(wǎng)圖層屬性信息中獲得每個格網(wǎng)分維信息，最后進行歸并與綜合。

河道寬度 對照大遼河流域水系圖，以所有采樣點所在的橫截面長度的平均值作為大遼河干流的河道寬度，一級支流以起點、終點和中間點3點所在橫截面長度的平均值作為河流寬度，二級支流取每條河中間點橫截面長度作為河流寬度，三級以下河流寬度則以隨機采樣（樣點設20余條河流）平均值作為代表。以此獲得整個流域的河道寬度。利用ArcGIS9.3的圖層疊置功能，將河流寬度信息輸入到對應位置格網(wǎng)中。

根據(jù)這3個方面的分布特點，將小流域單元進行水體類型分析（見圖2）。

圖2 小流域水體類型圖Fig.2 The water body types of watershed units

4.3 人類干擾特征識別

根據(jù)2007年TM影像解譯得到大遼河流域土地利用圖，各土地利用類型對生態(tài)環(huán)境施加的消極影響程度，計算各個評價單元的土地利用指數(shù)［21］，具體計算方法為：

式中：Lj為區(qū)域土地利用程度的綜合指數(shù)；Ai為研究區(qū)內第i種土地利用類型的分級指數(shù)，其中，各個類型土地利用的分級指數(shù)根據(jù)表1確定；Ci為研究區(qū)內第i種土地利用類型的面積百分比，n為土地利用類型數(shù)量。

表1 各個類型土地利用的分級指數(shù)Table 1 The classification index of various types of land use

根據(jù)所計算的土地利用指數(shù)，劃分不同區(qū)域的土地利用的影響強度大小，如表2所示。

表2 土地利用指數(shù)的影響強度分級表Table 2 Impact strength grading of land use index

按照上述方法計算各個小流域單元的土地利用指數(shù)，得到人類干擾強度圖3所示。

圖3 人類干擾強度分布圖Fig.3 Distribution of human disturbance

4.4 水生態(tài)系統(tǒng)功能評估

水生態(tài)功能重要性現(xiàn)狀評估采用指標體系評價的方法進行，在選取指標時，遵循主導性、差異性、綜合性和實用性，以及所需指標的可獲取性等原則，建立水生態(tài)功能評估指標體系（見表3）。建立評價指標體系后，在對各單項生態(tài)功能進行評價時，廣泛征詢專家意見并且對大遼河口區(qū)域進行調查分析。根據(jù)各個指標在評價中的影響力差異以及該地區(qū)的實際情況，主要分三類人員進行賦值：第一類為科研專家；第二類為研究區(qū)內的管理人員；第三類為調查人員，最后用加權平均法計算各評價指標的分值。

表3 大遼河流域水生態(tài)系統(tǒng)功能評價指標體系Table 3 Evaluation index system of water ecosystem functions in Liaohe valley

上述評估方法是建立在流域水生態(tài)系統(tǒng)特征現(xiàn)狀的基礎上，在評估過程中往往忽略了一些重要功能區(qū)以及人類對一些重要功能的期望目標。因此，需要參照流域歷史狀況及功能重要性特征，對以現(xiàn)狀指標為基礎的功能評估進行適當修正，以此得出能夠確定水生態(tài)系統(tǒng)保護目標的水生態(tài)功能評價結果，有利于確定流域內重要生態(tài)功能區(qū)，并實現(xiàn)相應的保護及恢復措施，并以此確定流域水生態(tài)功能保護的目標及制定正確的管理措施，修正技術路線見圖4。表3中5種功能評價指標經(jīng)過加權賦分后，統(tǒng)一為介于0～10之間的值，在它們的基礎之上，采用求和的方法，計算水生態(tài)系統(tǒng)綜合指數(shù)。根據(jù)大遼河流域整個區(qū)域的實際生態(tài)特點，經(jīng)過多次試驗，最終制定了表4中的等級劃分標準。利用ArcGIS9.3軟件中的空間疊加分析功能，將各單項功能矢量圖進行綜合疊加，如圖5所示，得到大遼河流域水生態(tài)系統(tǒng)功能重要性評價圖。

表4 水生態(tài)系統(tǒng)功能綜合評估分級表Table 4 Grade assignment of the comprehensive evaluation of water ecosystem functions

圖4 大遼河流域水生態(tài)系統(tǒng)功能重要性評估修正技術路線Fig.4 The modified technology roadmap of water ecosystem functions regionalization for Liaohe drainage area

圖5 大遼河流域水生態(tài)系統(tǒng)功能評價圖Fig.5 Evaluation of water ecosystem functions in Liaohe

5 大遼河流域水生態(tài)功能區(qū)劃分

分區(qū)命名要能夠準確體現(xiàn)研究區(qū)域內各分區(qū)的主要特征，即能夠體現(xiàn)水體的主要類型；能夠體現(xiàn)研究區(qū)域的生態(tài)功能類型；能夠反映人類干擾對生態(tài)環(huán)境的影響程度；同一級別生態(tài)功能區(qū)的名稱要相互對應。

功能區(qū)的命名方式主要按照“水體類型＋干擾類型＋功能類型”進行。首先根據(jù)水體類型，將其劃分為河流、水庫、濕地以及河口等水體類型；其次，根據(jù)人類干擾程度，分為城鎮(zhèn)地區(qū)、農(nóng)村地區(qū)等不同類型的生境地區(qū)；最后，根據(jù)水生生物多樣性維持功能、野生生物棲息地維持功能、濱岸帶支持功能、水環(huán)境支持功能、水文支持功能5個單項功能的評估，確定主導功能，以此將大遼河流域劃分成了5個一級區(qū)和13個生態(tài)亞區(qū)，結果如圖6和表5所示。

圖6 大遼河流域水生態(tài)系統(tǒng)功能分區(qū)圖Fig.6 Water ecosystem functions zoning for Liaohe drainage area

表5 大遼河流域水生態(tài)功能區(qū)劃表Table 5 Water ecosystem functions zoning for Liaohe drainage area

6 區(qū)劃結果分析

本次研究將大遼河流域劃分成了4類重要生態(tài)功能區(qū)和1類低生態(tài)功能區(qū)。

重要生態(tài)功能區(qū)生態(tài)敏感性較強，穩(wěn)定性較差，易受外界因素干擾。但是這類功能區(qū)對維持生態(tài)平衡、促進社會與經(jīng)濟持續(xù)發(fā)展等生態(tài)系統(tǒng)服務方面發(fā)揮了重要的作用。在這些生態(tài)區(qū)域內，人類應該禁止所有可能對生態(tài)功能造成破壞的開發(fā)活動和人為破壞活動，對那些己經(jīng)受到破壞的生態(tài)環(huán)境要加強生態(tài)修復工作；限制人口增長；改變傳統(tǒng)的生產(chǎn)經(jīng)營方式，走生態(tài)經(jīng)濟型可持續(xù)發(fā)展的道路；可以有針對性地建立自然保護區(qū)等。

低生態(tài)功能區(qū)的生態(tài)環(huán)境穩(wěn)定性比較好，它即可以承受一定的人類活動影響，又有一定的生態(tài)服務功能，這樣的功能區(qū)需要人們來引導它的開發(fā)方向，對開發(fā)利用加以一定限制，否則，就可能產(chǎn)生生態(tài)災害。

對于Ⅰ水生生物多樣性維持功能區(qū)和Ⅱ野生生物棲息地維持功能區(qū)，要加強河口濕地的保護力度，盡快建立濕地生態(tài)自然保護區(qū)，維護濕地生物多樣性；對已有的自然保護區(qū)要繼續(xù)加強保護，提高生物多樣性；加大退耕還林的力度，保護森林資源，提高保持水土、涵養(yǎng)水源的功能；加強陸域入海污染源的治理工作，使區(qū)域排放確保達標，改善濕地水環(huán)境質量。發(fā)揮區(qū)域資源優(yōu)勢，發(fā)展好生態(tài)農(nóng)業(yè)，并積極主動地發(fā)展生態(tài)旅游業(yè)。

對于Ⅲ河岸帶支持功能區(qū)可以采取小流域綜合治理的辦法，減少水土流失，關鍵是預防人為干擾活動造成新的水土流失；生態(tài)公益林建設要加強，建立一定的水源涵養(yǎng)保護區(qū)，以此來提高森林覆蓋率和水源涵養(yǎng)能力。

對于Ⅳ水環(huán)境支持功能區(qū)要杜絕污染型企業(yè)的建設，禁止對自然資源的掠奪性開發(fā)，把生態(tài)破壞與環(huán)境污染的源頭給徹底切斷；進一步加強環(huán)境的綜合整治工作，加強鞏固環(huán)境基礎設施的建設，以此減輕對水體的污染。

7 結語

本次研究為有關部門對大遼河流域的生態(tài)管理提供了一定的科學依據(jù)，在以后的研究中，還需繼續(xù)強化分區(qū)技術，使流域水生態(tài)系統(tǒng)分區(qū)體系更加完善。同時，水生態(tài)系統(tǒng)的特征是聯(lián)系的、動態(tài)的、復雜的、不確定的，因此應該將生態(tài)系統(tǒng)結構和功能相結合，對兩者的互相聯(lián)系與互相制約的機理進行辨析、維護和調整研究并付諸于實踐，進一步探索適合大遼河流域的生態(tài)系統(tǒng)管理模式。

［1］ Crowly J M.Biogeography［J］.Cannadian Geographer，1967，11（4）：312－326.

［2］ Bailey R G.Ecoregions of the United States（map）.1∶7500000［M］.Ogden，Utah：USDA Forest Service Intermountain Region，1976.

［3］ McNab W H，Avers P E.Ecological subregions of the United States：section descriptions［R］.Washington DC：USDA Forest Service，1994.

［4］ Omernik J M.Ecoregions of the Conterminous United States（Map Supplement）［J］.Annals of the Association of American Geographers，1987，77（1）：118－125.

［5］ Abell R，Thicme M L，Revenga C，et al.Freshwater ecoregions of the world A new map of biogeographic units for freshwater biodiversity conservation［J］.Bioscience，2008，58（5）：403－414.

［6］ 于海燕.錢塘江流域生態(tài)功能區(qū)劃研究［D］.浙江：浙江大學，2008.

［7］ Wells F，Peter N.An examination of an aquatic ecoregion protocol for Australia［R］.Canberra：Australian and New Zealand Environment and Conservation Council（ANZECC），1997.

［8］ 譚偉.《歐盟水框架指令》及其啟示［J］.法學雜志，2010，6：118－120.

［9］ 趙銳峰，周華榮，肖篤寧，等.塔里木河中下游地區(qū)濕地景觀格局變化［J］.生態(tài)學報，2006，26（10）：3470－3478.

［10］ 燕乃玲，趙秀華，虞孝感，等.長江源區(qū)生態(tài)功能區(qū)劃與生態(tài)系統(tǒng)管理［J］.長江流域資源與環(huán)境，2006，15（5）：598－602.

［11］ 于麗霞.大遼河口生態(tài)環(huán)境綜合評價［D］.青島：中國海洋大學，2011.

［12］ 王仰麟.格局與過程—景觀生態(tài)學的理論前沿［M］.北京：中國科學技術出版社，1995：437－411.

［13］ 陳爽，馬安青，李正炎，等.基于RS／GIS的大遼河口濕地景觀格局時空變化研究［J］.中國環(huán)境監(jiān)測，2011，4（2）：34－40.

［14］ 王以彭，李結松，劉立元.層次分析法在確定評價指標權重系數(shù)中的應用［J］.第一軍醫(yī)大學學報，1999，19（4）：377－379.

［15］ 李福建，馬安青，丁原東，等.基于RS與GIS技術的地質災害危險性評價－以青島市嶗山區(qū)為例［J］.中國海洋大學學報：自然科學版，2010，6（40）：47－52.

［16］ 梅昌艮.基于GDM－AHP的模糊綜合評價法在生態(tài)敏感性評價中的應用研究［D］.四川：西南交通大學，2008：39－41.

［17］ 劉堯讓，于法展，李淑芬.蘇北山丘區(qū)生態(tài)環(huán)境現(xiàn)狀評價［J］.海南師范大學學報：自然科學版，2010，23（2）：214－218.

［18］ 陳爽，王其翔，馬安青，等.基于GIS的大遼河口水生態(tài)系統(tǒng)服務功能重要性評價［J］.中國海洋大學學報：自然科學版，2012，42（9）：84－89.

［19］ 張洪軍.生態(tài)規(guī)劃一尺度，空間布局與可持續(xù)發(fā)展［M］.北京：化學工業(yè)出版社，2007.

［20］ 王秀蘭，包玉海.土地利用動態(tài)變化研究方法探討［J］.地理科學進展，1999，18（1）：81－87.