褚俊英，嚴登華，周祖昊，吳 迪，劉 琳

（1.中國水利水電科學研究院 流域水循環(huán)模擬與調(diào)控國家重點實驗室，北京 100038；2.中建水務(wù)環(huán)保有限公司，北京 100070）

1 研究背景

隨著城鎮(zhèn)化進程的加快，我國城市河湖水系正面臨著水資源短缺、水質(zhì)惡化、水體萎縮與調(diào)蓄能力降低、河湖水力和生態(tài)聯(lián)系阻隔、濱湖地帶過度開發(fā)、生物多樣性減少等突出問題，河湖的健康生命受到嚴重威脅，并對人類社會的長遠發(fā)展帶來嚴峻挑戰(zhàn)。城市河湖生態(tài)系統(tǒng)是一個完整的有機整體，水是生態(tài)過程的重要驅(qū)動因子，生態(tài)流量是河湖健康維系的重要條件［1］?？茖W給出城市河湖生態(tài)流量閾值尤為重要，是維護城市河湖基本功能，實施最嚴格水資源管理的基石［2］，有利于促進城市生態(tài)系統(tǒng)與社會經(jīng)濟的協(xié)調(diào)發(fā)展。

在城市河湖生態(tài)流量計算方面，國內(nèi)外學者當前提出多種計算方法，如水文學、水力學、棲息地法以及整體法等［3］，河湖生態(tài)流量的研究日益注重水分-生態(tài)的耦合作用機理、人類活動對生態(tài)用水過程影響機制且對資料的要求較高［4］。如在棲息地法中研究者可通過河湖中存在的敏感性、指示性的生物物種調(diào)查，分析確定其生活史過程中（如棲息、產(chǎn)卵等）所需要的水深、流速等指標要求，從而推求河湖生態(tài)流量。Marsili-Libelli等［5］采用鯉魚和歐洲鰱魚作為指示物種，依據(jù)適宜水流和污染物濃度曲線確定生態(tài)流量。徐宗學等［6］依據(jù)魚類洄游的流速范圍，采用生態(tài)水力學法計算生態(tài)流量。由于我國城市河湖缺乏長期的水量、水質(zhì)與水生生態(tài)系統(tǒng)結(jié)構(gòu)與群落的觀測資料，多數(shù)河湖都已渠道化且受到人工調(diào)控影響嚴重，這些方法在實踐中應(yīng)用通常面臨數(shù)據(jù)缺失、適應(yīng)性和代表性有待提升等瓶頸，其實踐應(yīng)用的范圍與效果有限［7］。

城市河湖相對自然河湖在服務(wù)功能方面有明顯差異，通常大部分喪失了傳統(tǒng)上的灌溉、飲水、航運、維持河口鹽分平衡等功能，多保留了排水、景觀、娛樂、防洪、沖沙等功能［8］。為體現(xiàn)城市河湖水系這些特定的生態(tài)服務(wù)功能，面向生態(tài)服務(wù)功能的河湖生態(tài)流量計算方法日益得到發(fā)展。如賈仰文等［9］考慮生態(tài)基流、魚類產(chǎn)卵、輸沙、濕地與景觀、滲漏與蒸發(fā)等多個要素，綜合給出了渭河流域城市河湖生態(tài)流量；楊毅等［10］體現(xiàn)基本生態(tài)需水、凈化、景觀娛樂等功能要求開展了資料有限條件下北京市河流的生態(tài)需水計算。這些方法當前主要針對城市特定的河湖構(gòu)建，而河湖多種復(fù)合生態(tài)服務(wù)功能的合理辨識、多種功能所需生態(tài)流量的量化方法等有待于進一步提升，河流與湖泊、濕地之間的耦合關(guān)系也需要進一步明確。本文針對城市河湖的特點，基于綜合功能辨識構(gòu)建了河湖生態(tài)流量計算與保障程度分析的總體框架與量化模型，選取株洲市區(qū)的河湖水系進行了典型案例分析，可為城市河流生態(tài)流量的保障與調(diào)控提供基礎(chǔ)支撐。

2 基本概念與總體框架

2.1 基本概念 當前，城市河湖生態(tài)流量尚缺乏明確的定義，生態(tài)用水、生態(tài)耗水、生態(tài)缺水、生態(tài)基流、生態(tài)需水、環(huán)境需水及生態(tài)環(huán)境用水等多個概念尚未統(tǒng)一。生態(tài)流量對河湖生態(tài)系統(tǒng)所需物質(zhì)、能量、信息及物種輸移有重要意義，為了保證河湖健康與生態(tài)系統(tǒng)功能的發(fā)揮，要求河湖生態(tài)流量的大小、情勢與水質(zhì)不能超過生態(tài)系統(tǒng)的耐受性范圍，生態(tài)系統(tǒng)仍然能夠順利完成其生活史。

城市河湖生態(tài)流量是指城市河湖健康所需要的、符合水質(zhì)要求且滿足生態(tài)需求節(jié)律的流動水量。河湖生態(tài)流量由消耗分項和非消耗分項構(gòu)成［11］。其中，消耗分項主要是為實現(xiàn)河湖生態(tài)系統(tǒng)蒸散發(fā)耗水和滲漏等功能；非消耗分項主要是滿足河湖生態(tài)健康所需的水動力學過程以維持連通性。在消耗分項和非消耗分項基礎(chǔ)上，河湖生態(tài)流量可進一步細分成存量水和通量水。其中，存量水是指賦存于河道與湖泊且具有一定體積的連續(xù)水分體，通量水是指單位時間流經(jīng)某一單位截面積的水量（如流量或換水量）。存量主要以河湖特征水深體現(xiàn)，可基于湖泊水位-容積關(guān)系曲線和水位-面積關(guān)系曲線，由特征水位確定相應(yīng)的容積和水面面積相除所得，如滿足河湖景觀功能水深等；通量表現(xiàn)為河湖的流動水量，如河流生態(tài)基流、河湖維持自凈流量、河流輸沙功能流量以及河流景觀流量等。消耗分項主要由通量構(gòu)成，非消耗分項則包括存量和通量。同一流域或區(qū)域的水生態(tài)系統(tǒng)分布格局及生態(tài)用水的時空特征不同，產(chǎn)水區(qū)域也不盡相同，需要進行整合計算，以確定城市河湖主要控制斷面或控制點的生態(tài)流量過程。通常，城市河道與湖泊水體以外的農(nóng)田灌溉、綠地澆灌、城市供水等則不計入河湖生態(tài)流量核定的范疇。

圖1 城市河湖生態(tài)流量分析的總體框架

2.2 總體框架 基于綜合功能辨識的城市河湖生態(tài)流量分析的總體框架如圖1所示，具體如下：（1）確定城市河湖的主要控制斷面與控制節(jié)點。針對各城市河湖水系的特點，綜合考慮支流匯入、上下游關(guān)聯(lián)、生態(tài)功能要求等多種因素，選取河湖生態(tài)流量計算的代表性控制斷面與節(jié)點。每個斷面或節(jié)點控制的集水區(qū)都可視為獨立單元，能夠發(fā)揮特定的生態(tài)服務(wù)功能；（2）城市河湖的保護目標與主導功能的識別。逐個識別各控制斷面與控制節(jié)點的主導生態(tài)服務(wù)功能，如敏感生物、景觀、娛樂、防洪、沖沙以及航運等方面，明確生態(tài)保護目標；（3）城市河湖生態(tài)流量過程線的綜合計算與優(yōu)化分析。面向各控制斷面與控制節(jié)點的功能和保護目標，建立上下游的產(chǎn)匯流水力聯(lián)系，借助水文、水量、水質(zhì)、生態(tài)等多方面數(shù)據(jù)分析，進行河湖生態(tài)流量消耗與非消耗分項、存量與通量的計算，綜合優(yōu)化得到河湖生態(tài)流量過程線、流量構(gòu)成（即流量中消耗與非消耗量、存量與增量、各項生態(tài)服務(wù)功能所需流量）與變異性系數(shù)；（4）提出城市河湖生態(tài)流量的保障程度和調(diào)控要求。分析不同頻率的來水過程，計算得出不同來水頻率條件下城市河湖生態(tài)流量的保障程度，提出生態(tài)流量的調(diào)控指標要求以及調(diào)控時機，為城市河湖的健康管理提供決策支持。

3 計算模型

3.1 模型結(jié)構(gòu) 基于綜合功能辨識的城市河湖生態(tài)流量計算與保障程度分析模型（Model for urban ecological streamflow calculation and guarantee degree analysis based on systhesized ecosystem service function identification，簡稱MUSEF模型）的基本結(jié)構(gòu)如圖2所示。該模型面向城市河湖多種生態(tài)服務(wù)功能目標，綜合體現(xiàn)河流斷面與湖泊濕地水力聯(lián)系，權(quán)衡生態(tài)流量消耗分項與非消耗分項、存量與通量、水量與水質(zhì)之間的關(guān)系，量化得出城市河湖生態(tài)流量大小、過程、構(gòu)成、空間分布、滿足程度與調(diào)控要求，為城市河湖生態(tài)流量的保障與健康管理決策提供支持。

圖2 基于綜合功能辨識的城市河流生態(tài)流量計算模型結(jié)構(gòu)（MUSEF模型）

3.2 調(diào)控機制 河段水平衡關(guān)系為上游匯入流量、自身流域匯入量、補給量、降水量與人工補水量之和等于社會經(jīng)濟最大耗水限值與生態(tài)流量（如圖2所示，具體計算如下式所示）。生態(tài)流量指標閾值尤為關(guān)鍵，該指標的確定將對社會經(jīng)濟（如工業(yè)、農(nóng)業(yè)等）用耗水過程進行限制，從而最大限度減少人類水資源開發(fā)利用對河湖生態(tài)系統(tǒng)的破壞［2］。為維系生態(tài)流量，人類可以減少本河段社會經(jīng)濟耗水，也可從其他水源如地表水、地下水或非常規(guī)水源對河段進行補水［12］，具體調(diào)控措施的選取需要進行技術(shù)經(jīng)濟的綜合權(quán)衡。

式中：St為第t時間上游匯入流量，m3/s；ZQt為第t時間自身流域匯入流量，m3/s；Tt為第t時間人工補水流量，m3/s；PQt為第t時間水體降水補給流量，m3/s，KQt為第t時間其他補給流量，m3/s；Ut為第t時間社會經(jīng)濟最大耗水流量限值，m3/s；Wt為第t時間河流生態(tài)流量閾值，m3/s。

3.3 生態(tài)流量過程線計算 城市河流生態(tài)流量過程線由河段生態(tài)流量以及與河段連接的湖泊/濕地生態(tài)流量組合構(gòu)成，包括消耗和非消耗分項兩部分，具體如下：（1）河流生態(tài)流量的消耗分項包括河流蒸散通量、滲漏通量以及其他消耗（如生物消耗）等通量；非消耗分項包括河流生態(tài)系統(tǒng)基流、河流維持自凈流量、河流輸沙功能流量、河流景觀與其他功能流量等通量和存量；（2）湖泊/濕地生態(tài)流量的消耗分項包括湖泊/濕地蒸散量、滲漏量和其他消耗等通量，非消耗分項為自凈與其他功能水量。鑒于非消耗分項實現(xiàn)多項生態(tài)服務(wù)功能具有重疊性，取外包值（即max函數(shù)）計算，如下式所示：

式中：Qt為第t時間河流生態(tài)流量，m3/s；Mt為第t時間與河流連接的湖泊濕地的生態(tài)流量，m3/s；Qet、Qgt和Qnt分別為第t時間河流蒸散通量、河道滲漏通量以及其他消耗通量，m3/s；Qmt、Qat、Qbt、Qct、Qdt分別為第t時間河流非消耗生態(tài)流量、河流生態(tài)系統(tǒng)基流、河流維持自凈流量、河流輸沙功能流量、河流景觀與其他功能流量等通量，m3/s；Met、Mgt、Mnt和Mbt分別為湖泊/濕地蒸散通量、湖泊/濕地滲漏通量、湖泊/濕地其他消耗通量以及湖泊/濕地自凈與其他水量等通量，m3/s。

3.3.1 消耗分項計算

（1）河湖水面蒸發(fā)與濕地植被蒸散發(fā)通量。蒸散發(fā)作為城市河湖水循環(huán)過程中的一個重要環(huán)節(jié)，是水量的損失項，是河湖生態(tài)流量的組成部分和消耗項。通常，濕地植被的蒸散量包括植物同化過程耗水、植物體內(nèi)包含的水分、濕地植株表面蒸發(fā)耗水等部分，其中后兩部分所占比重較大［13］。在數(shù)據(jù)稀疏時濕地植被與水面面積可按照一定的比例計算［9］。河道和湖泊的水面面積以及濕地植被面積通常是水深（Hb和Hq）的函數(shù)，具體計算如下：

式中：Qet、Met分別為河湖水體第t時間的蒸散發(fā)量，m3/s；E1t、E2t分別為河湖水面第t時間水面蒸發(fā)強度、濕地植被蒸散發(fā)強度，mm/s；Aqt、Amt分別為第t時間的河道和湖泊的水面面積，m2；A1bt、A2bt分別為第t時間的河流濕地與湖泊濕地植被面積，m2。

（2）河湖滲漏通量。滲漏是河湖水量自然損失的方式，當河湖水位高于地下水位時，河湖水體通過底部滲漏和岸邊側(cè)滲補給地下水。考慮到城市河湖底部通常會有護坡和防滲處理，滲透系數(shù)非常小。對于未經(jīng)防滲的自然河湖水體，可采用如下公式計算［14］。自然河湖的滲透系數(shù)通常受到土壤類型、河流形態(tài)、河道水位與地下水位差、地質(zhì)條件等多種因素影響，通常采用水量平衡法、地下水動力學方法或經(jīng)驗公式法等方法計算，或通過河道滲漏現(xiàn)場試驗以及動態(tài)模擬等方法進行確定［15］。

式中：Qgt、Mgt分別為第t時間的河道和湖泊滲漏水量，m3/s；Kq、Km分別為河道和湖泊滲透系數(shù)，m/s；Iq、Im為河道和湖泊坡降，‰；Lqt、Lmt分別為第t時間的河流和湖泊過水斷面面積，m2；如果僅考慮垂直滲漏，則滲漏面積即為河湖的水面面積，可采用如下公式計算［11］。

（3）其他消耗通量。城市河湖其他消耗量（Qnt和Mnt）通常較小，可忽略不計。對于特殊的河湖生態(tài)系統(tǒng)，需根據(jù)實際情況采用經(jīng)驗公式或量化方法確定。

3.3.2 非消耗分項計算 城市河湖生態(tài)流量非消耗分項的計算主要包括兩部分，一是存量特征水深的計算；二是體現(xiàn)特征水深的通量計算。其中，存量特征水深是指湖泊生態(tài)系統(tǒng)最低生態(tài)水深、河湖景觀與其他功能水深；體現(xiàn)特征水深的通量包括河流生態(tài)系統(tǒng)基流、河湖維持自凈流量、河流輸沙功能流量、河流景觀娛樂與其他功能流量。

（1）存量的特征水深計算。

（i）湖泊生態(tài)系統(tǒng)基本水深。維持湖泊形態(tài)以及水生生物生存與棲息等生態(tài)服務(wù)功能不喪失，需要一定的湖泊蓄水量（存量）。對于湖泊水生生態(tài)系統(tǒng)而言，枯水期實際可保持的滿足其生存需求的水深尤為重要。計算湖泊維系生態(tài)系統(tǒng)功能的基本生態(tài)水深Hs可以采用湖泊形態(tài)法、生物空間法等方法。其中，湖泊形態(tài)法主要基于歷史資料，量化分析湖泊水面面積變化率與湖泊水深關(guān)系曲線，依據(jù)該曲線尋找枯水期水深附近面積變化率最大值對應(yīng)的數(shù)值作為基本生態(tài)水深，如公式所示。要維持該湖泊基本生態(tài)水深，需要計算湖泊蒸發(fā)、滲漏等消耗分項（如前所述），作為湖泊所需的生態(tài)流量指標。

式中：F為湖面面積，m2；f為函數(shù)關(guān)系；Hs為湖泊基本生態(tài)水深，m；Hmin為湖泊多年平均最低水深，m；a、b分別為和湖泊水位變幅相比較小的正數(shù)，m。

（ii）河湖景觀與其他功能基本水深。城市河湖系統(tǒng)通常受到水閘和橡皮壩的人工擾動，流速和流量受到控制。保持特定的水深通常是城市河湖實現(xiàn)景觀與其他功能的重要途徑［1］。基于數(shù)據(jù)資料，綜合體現(xiàn)景觀、娛樂、航運、灌溉、發(fā)電等多方面功能要求，確定河湖景觀與其他功能基本水深（Hw和Hq）。

（2）通量計算。

（i）河流生態(tài)基流。維持河流基本水生生態(tài)系統(tǒng)的穩(wěn)定，確保一定的河流基流尤為重要。河流生態(tài)基流多采用水文學方法計算，如美國、法國多采用多年平均流量的10%，作為大多數(shù)水生生物生存所推薦的最低流量；我國《水利水電建設(shè)項目水資源論證導則》（SL525-2011）提出原則上按多年平均流量的10%～20%確定。隨著城市河湖的水文、水量、水質(zhì)與水生生態(tài)系統(tǒng)（如浮游植物、浮游動物、底棲動物等）與群落觀測資料的不斷完善，突出水分與生態(tài)系統(tǒng)響應(yīng)機理的生態(tài)流量計算方法將得到更多應(yīng)用［16］。

式中：Qat為第t時間的河流生態(tài)基流，m3/s；Rr為多年平均徑流量，m3/s；Cat為第t時間的生態(tài)基流的比重系數(shù)，無量綱；Prt為多年平均條件第t時間的降雨量，mm/s；Kr為降雨徑流系數(shù)，無量綱；Dq為河段陸域的匯流面積，m2。

（ii）河湖維持自凈流量。河流自凈流量是保護河流達到水質(zhì)目標必須維持的最小流量。我國《地方水污染物排放標準的技術(shù)原則和方法》（GB3839-83），規(guī)定一般河流采用近10年最枯月平均流量或90%保證率最枯月平均流量，用于計算污染物允許排放量?！端蚣{污能力計算規(guī)程》（GB/T25173）提出“計算河流水域納污能力，應(yīng)采用90%保證率最枯月平均流量或近10年最枯月平均流量作為設(shè)計流量。”除了統(tǒng)計學方法外，研究者還構(gòu)建了水質(zhì)模型模擬污染物濃度，結(jié)合本底情況計算得到所需的維持自凈流量。如康然然等［17］利用QUAL2K建立了沈陽城區(qū)南運河水系的豐、枯水期的水質(zhì)模型，以COD和氨氮指標達到水質(zhì)目標為約束，不斷調(diào)整進水閘水量，直到得到維持自凈的最小生態(tài)流量。

對于城市湖泊而言，采用新鮮水源對水體進行稀釋不利于污染物的源頭控制，但可促進河湖水體流動而提高水環(huán)境承載力，也成為當前城市河湖管理舉措，例如北京六海換水頻率為每年0.5～1.25次［18］，渭河流域主要城市湖泊換水周期為60 d［9］，長沙市公園湖泊換水頻率約每年1.7次［19］。換水周期或頻率的大小通常受到城市水資源配置調(diào)度、湖泊水系景觀娛樂等服務(wù)功能要求的影響［20］。隨著湖泊水污染防治與生態(tài)修復(fù)技術(shù)的發(fā)展及應(yīng)用，該水量將逐漸減少至消除，具體計算如下式所示：

式中：Mbt為第t時間的城市湖泊維持環(huán)自凈水量，m3/s；Hb為湖泊最小生態(tài)水深，m；Uh為湖泊換水系數(shù)，通常根據(jù)換水周期、換水頻率或換水量比率推算，次/s；Hw為湖泊景觀與其他功能基本水深，m；Hs為湖泊基本生態(tài)水深，m。

（iii）河流輸沙功能流量。河流輸沙功能通常受到水量和含沙量兩方面影響，主要在汛期完成。為保障河流輸沙功能的實現(xiàn)，生態(tài)流量計算如下式所示：

式中：Qct為第t時間的河流輸沙水量，m3/s；St為第t時間的河流輸沙量，t/s；Cmax為長系列最大河流平均含沙量，t/m3；Ct為第t時間河流含沙量，t/m3。

（iv）河流景觀與其他功能流量。除了上述生態(tài)服務(wù)功能外，健康的河流還需要滿足河流景觀與其他特定功能如航運、發(fā)電等。河流景觀與其他功能流量的計算通?？刹捎媒?jīng)驗法、水力學法等多種方法進行計算。其中，水力學法可采用下式計算：

式中：N為河道天然糙率，通常根據(jù)現(xiàn)場勘測的河床組成、斷面形狀、岸坡及河床中卵石礫石的分布、疏密度等情況確定，無量綱；J為坡降，‰；R為生態(tài)水力學半徑，當河流寬深比較大時，可近似等于河流景觀與其他功能基本水深Hq，m；B為河道寬度，m。

3.4 生態(tài)流量檢驗與優(yōu)化分析 對于城市河湖生態(tài)流量合理性檢驗通常采用多年平均流量的合理百分比確定，枯水期和豐水期比重應(yīng)有所不同，需要根據(jù)各城市河湖水系的實際情況確定。通常，在數(shù)據(jù)缺失條件下，城市河湖的部分生態(tài)服務(wù)功能無法量化辨識，為了最大限度的確保河湖生態(tài)系統(tǒng)的健康，通常取兩者的外包絡(luò)線確定，如下式所示。

式中：Qtt為第t時間的河湖生態(tài)流量檢驗值，m3/s；Ntt第t時間的生態(tài)流量檢驗系數(shù)，無量綱；WTt為第t時間綜合優(yōu)化后的河湖生態(tài)流量，m3/s。

3.5 過程變異性 城市河湖生態(tài)系統(tǒng)具有特定的水文節(jié)律要求，在不同時期河湖生態(tài)流量的大小通常有所差異。這種過程變異性程度可采用過程變異性系數(shù)來體現(xiàn)，以反應(yīng)不同時間生態(tài)流量指標偏移總體樣本均值的大小，具體計算如下式所示：

式中：Kwm，t為第m個控制斷面或節(jié)點第t時間的河湖生態(tài)流量過程變異性系數(shù)，無量綱；WAm為第m個控制斷面或節(jié)點河湖生態(tài)流量從t=1到T時刻的均值，m3/s；KW為河湖生態(tài)流量所有控制斷面或節(jié)點總過程變異性系數(shù)，無量綱；Np為控制斷面或節(jié)點總數(shù)，個。

3.6 不同來水頻率條件下的保障程度 不同來水頻率條件下城市河湖各斷面或控制節(jié)點生態(tài)流量的保障程度計算如下式所示。其中，城市河湖生態(tài)流量總體保障程度等于所有河流控制斷面與節(jié)點保障程度的平均值。基于生態(tài)系統(tǒng)對干擾的抗性能力特征，通常采用月平均流量進行生態(tài)需水保障程度分析［21］。

式中：BWx分別為x來水頻率條件下的總保障程度，%；Bwx，t為x來水頻率條件下的第t時間的保障程度，%；Rx，t為x來水頻率條件下第t時間的來水流量，m3/s；Px，t為x來水頻率條件下第t時間的降雨量，mm/s；x為來水頻率，根據(jù)不同城市河湖水系的特點確定，如多年平均、75%、90%和95%等不同來水頻率條件。

3.7 模型主要參數(shù)與輸入特征 MUSEF模型的輸入輸出與主要參數(shù)如表1所示。河湖生態(tài)流量各分項所選擇的計算方法不同，要求的精細化程度不同，所需的輸入數(shù)據(jù)和參數(shù)也差異較大?？傮w上看，降水、蒸發(fā)、滲漏以及河湖水體的特征參數(shù)（如面積、水深、水質(zhì)以及生態(tài)系統(tǒng)特征）的可得性與可靠性是影響城市河湖生態(tài)流量計算合理性的重要因素。

表1 基于綜合功能辨識的城市河流生態(tài)流量計算模型輸入輸出與主要參數(shù)

4 案例應(yīng)用

4.1 研究區(qū)域與基礎(chǔ)數(shù)據(jù) 株洲市位于湖南省東部，湘江中下游，屬于湘江生態(tài)經(jīng)濟帶，是湖南工業(yè)重鎮(zhèn)、南方交通樞紐，素有“火車拖來的城市”之稱，享有“中國電力機車的搖籃”和“中國電力機車之都”的美譽，是湖南省“一點一線”區(qū)域經(jīng)濟帶的重要城市，是國家重點開發(fā)區(qū)長株潭城市群（“金三角”）的重要組成部分。株洲市地處濕潤多雨丘陵區(qū)，河網(wǎng)縱橫、湖庫密布、水系發(fā)達，湘江東側(cè)（河東）地形破碎復(fù)雜，地勢高低起伏，分布有湘江的4個一級支流，即白石港、建寧港、楓溪港、霞灣港。

基于株洲市區(qū)河湖水系特點，研究確定了河湖生態(tài)流量的31個主要控制斷面和控制節(jié)點。其中，河流生態(tài)流量斷面21個（干流為11個，其余為支流），湖泊控制節(jié)點10個，如圖3所示。株洲市河湖生態(tài)流量計算的主要數(shù)據(jù)來源如下：湖泊、河流長度、面積與斷面形狀等基礎(chǔ)資料來自株洲市水利局的長期監(jiān)測和中國水利水電科學研究院2016年8月現(xiàn)場測定的結(jié)果（表2）。河湖滲透系數(shù)來源于株洲市水利水電規(guī)劃勘測設(shè)計院的勘測結(jié)果，均值為3.0×10-7m/s，其余數(shù)據(jù)見表3。降雨資料來源于國家氣象局株洲站1955—2016年共計62年逐月資料，降雨產(chǎn)流系數(shù)來自《株洲市水資源公報》，蒸發(fā)數(shù)據(jù)來源于《株洲統(tǒng)計年鑒》，如表4所示。

圖3 株洲市河湖生態(tài)流量的控制斷面與控制節(jié)點

表2 株洲市區(qū)河流生態(tài)流量計算斷面的參數(shù)與輸入數(shù)據(jù)

表3 湖泊濕地多年平均面積與最小生態(tài)水深數(shù)據(jù)

表4 株洲市區(qū)河流主要控制斷面生態(tài)流量過程線（單位：mm）

4.2 河流生態(tài)需水過程線、構(gòu)成與變異性 在確定株洲市河湖主要控制斷面和控制節(jié)點多重生態(tài)服務(wù)功能目標基礎(chǔ)上，經(jīng)MUSEF模型進行長系列計算，得到控制斷面與節(jié)點逐月生態(tài)流量過程線，如表5所示。可以看出，河流干流斷面汛期（4—9月）的流量總體大于非汛期，以滿足汛期特定生態(tài)服務(wù)功能的要求。支流入江口斷面（以霞灣港為例）生態(tài)流量的過程線與構(gòu)成如圖4所示?？傮w上，株洲市區(qū)河湖生態(tài)流量主要為非消耗生態(tài)流量，其中最主要的是景觀娛樂需水，其次是環(huán)境自凈流量，輸沙流量、生態(tài)基流所占比重較小。從生態(tài)流量的過程檢驗看，汛期生態(tài)流量過程有所優(yōu)化提升，屬偏安全計算。

表5 株洲市區(qū)河流主要控制斷面生態(tài)流量過程線（單位：m3/s）

圖4 支流入江口斷面生態(tài)流量構(gòu)成（霞灣港，S41斷面）

河湖生態(tài)流量的月變異系數(shù)反應(yīng)了逐月生態(tài)流量偏離全年尺度月均值的程度，計算結(jié)果如圖5所示。這些斷面可分成如下三類：（1）“二值型”，即汛期（4—9月）逐月生態(tài)流量基本相同，月變異系數(shù)較大；非汛期（其余月份）逐月生態(tài)流量基本相同，月變異系數(shù)較小，包括14個斷面（即S1-2，S11，S12，S13，S22，S3-1，S3-4，S3-5，S3-7，S32，S33，S34，S41和 S42）；（2）“一值型”，即汛期、非汛期逐月的生態(tài)流量過程大致相當，基本在月均值附近，如斷面S1-1，S1-3，S3-2，S3-3，S3-6和S31；（3）“分散型”，即逐月的生態(tài)流量過程不同，如斷面S21，主要原因是湖泊濕地生態(tài)用水所占比重較大，受到湖泊濕地蒸散發(fā)耗水、環(huán)境自凈等月過程變化的影響，生態(tài)流量逐月過程有明顯波動，具體如圖6所示。

圖5 株洲市區(qū)河流生態(tài)流量的月變異系數(shù)

圖6 株洲河流連接的所有湖泊生態(tài)水量月過程構(gòu)成

圖7 株洲河流各斷面生態(tài)需水逐月滿足程度

4.3 不同來水頻率條件下河湖生態(tài)流量保障程度與調(diào)控要求 在多年平均、75%、90%和95%來水頻率條件下，株洲市各河流斷面生態(tài)流量的總保障程度分別為98.8%、91.7%、79.0%和65.9%，具體斷面保障程度如圖7所示。通過各種水資源調(diào)控途徑所需增加的河湖生態(tài)流量，在多年平均、75%、90%和95%來水頻率條件下折合年水量分別為15.6萬、333.4萬、3975.8萬和4335.8萬m3。從調(diào)控時機看，河湖生態(tài)流量調(diào)控主要分布在非汛期，在多年平均條件下12月需增加生態(tài)流量；在75%來水頻率條件下2月和12月需增加生態(tài)流量；在90%和95%來水頻率條件下，1—3月、8—10月和12月需增加生態(tài)流量。從空間分布看，河湖生態(tài)流量所需增加最大的為白石港，其次為楓溪港，建寧港和霞灣港所需增加的生態(tài)流量較少。

5 結(jié)論

研究構(gòu)建了基于綜合功能辨識的城市河湖生態(tài)流量計算與保障程度分析模型，可為不同城市的河湖生態(tài)流量的安全保障和健康管理提供科學基礎(chǔ)。以株洲市區(qū)河湖為案例應(yīng)用的分析結(jié)果表明，在多年平均、75%、90%和95%來水頻率條件下，株洲市河湖生態(tài)流量的年均保障程度分別為98.8%、91.7%、79.0%和65.9%。在此基礎(chǔ)上，量化提出了所需增加的生態(tài)流量指標要求和時空分布特征。

城市河湖生態(tài)流量隨著河湖水系特征、生態(tài)系統(tǒng)結(jié)構(gòu)和功能的差異化具有相當?shù)膹?fù)雜性。應(yīng)進一步加強城市河湖生態(tài)系統(tǒng)健康狀況調(diào)查與評估，對河流及其關(guān)聯(lián)的湖泊、濕地的水文、水質(zhì)、水生生物多樣性與群落演替（如浮游植物、浮游動物、底棲動物、魚類等）進行系統(tǒng)的本底調(diào)查與長系列的觀測，為河湖生態(tài)流量的精細化計算提供進一步支撐。應(yīng)加強河湖生態(tài)流量計算的參數(shù)不確定性分析，并結(jié)合城市社會經(jīng)濟發(fā)展耗水過程與趨勢，綜合考慮再生水、雨水、地表水和地下水，開展城市河湖生態(tài)流量的多水源優(yōu)化配置與調(diào)度，為變化環(huán)境下城市河湖生態(tài)流量閾值的定量識別與動態(tài)保障提供決策支持。