馮慧娟 李建玲 王喆 李丹

1 概 況

1.1 工程基本情況

鳳山水庫位于貴州省黔南州福泉市重安江干流上游河段。水庫總庫容1.04 億m3，死水位868 m，正常蓄水位909 m。水庫多年平均供水量為7 310萬m3，為多年調(diào)節(jié)水庫。工程任務以城鄉(xiāng)生活、工業(yè)供水為主，兼顧發(fā)電等綜合利用。

1.2 研究河段水生生態(tài)系統(tǒng)調(diào)查及代表性魚類生態(tài)習性

根據(jù)對鳳山水庫所在重安江干流河段及部分支流水生生態(tài)現(xiàn)狀調(diào)查及歷史資料收集，該區(qū)域分布有34 種魚類，隸屬于4 目8 科27 屬，以鯉形目為主，其次為鲇形目。產(chǎn)黏性卵魚類是調(diào)查水域最主要類群，多在卵石、礫石、砂質(zhì)或巖縫間隙等底質(zhì)中產(chǎn)卵。

調(diào)查水域，鲇、黃顙魚、白甲魚、大眼鱖等魚類具有重要種質(zhì)資源價值和漁業(yè)資源價值。其中，鲇、黃顙魚、白甲魚繁殖季節(jié)為4—6 月，大眼鱖繁殖季節(jié)為3—4 月。

魚類產(chǎn)卵需一定流速刺激，研究河段魚類產(chǎn)卵流速一般不小于0.2 m/s，水深一般不小于0.5 m，水溫不低于18～20 ℃。

2 鳳山水庫下游減水河段一維河道模型構建

采用一維非恒定流水動力模型HEC-RAS，模擬研究河段不同水文過程下流速、水深過程。與研究河段重要魚類種質(zhì)資源產(chǎn)卵期的棲息生境條件耦合，推薦適合魚類產(chǎn)卵期生態(tài)調(diào)度過程，以滿足魚類產(chǎn)卵需求。

2.1 控制方程

鳳山水庫下游減水河段模型為河道一維非恒定水流模擬，主要基于連續(xù)方程與運行方程進行迭代計算，通過河道上游流量、下游水位確定水位變化過程。主要計算公式如下。

式中B——河流橫斷面寬度，m；

z——河流水位， m；

Q——河流流量， m3/s；

t——時間，s；

x——河流斷面縱向間距，m；

q——單位河長流入或流出的流量，m2/s；

u——流速，m/s；

ρ——水的密度，kg/m3；

f——單位質(zhì)量力，m/s2；

p——壓強，Pa；

ν——運動黏滯系數(shù)，m2/s。

2.2 模型構建

本次模擬河段為重安江干流鳳山水庫壩址斷面至羊昌河斷面，河道總長約22.85 km。選取阿里堡河口上游斷面（壩下6 km）、三江口斷面（壩下10 km）、豆腐橋斷面（壩下14 km）、鴨爪壩斷面（壩下15 km）、羊昌河斷面（壩下22 km）5 個典型斷面進行模擬，各典型斷面位置分布如圖1 所示。

圖1 各典型斷面位置分布圖

采用一維非恒定流水動力模型HEC-RAS 計算典型水文過程。各斷面地形剖面由現(xiàn)場勘測獲得，模型上邊界為流量邊界，下邊界為水位邊界。

2.3 模型率定

提取模型計算結果，對比實測數(shù)據(jù)，調(diào)整糙率值，使得模擬值最大限度接近實測值，最終河道糙率系數(shù)取0.036～0.040，鳳山水庫主體設計河道糙率基本一致。對比模擬與實測值，鳳山水庫壩址斷面水位流量曲線關系模擬值與實測值擬合效果較好，如圖2 所示。

圖2 壩址斷面水位流量曲線關系擬合效果

3 生態(tài)流量的確定

根據(jù)生態(tài)環(huán)境部批復的鳳山水庫工程環(huán)境影響報告書，明確不同水期生態(tài)流量泄放要求。

3.1 生態(tài)流量計算方法

3.1.1 水文學法

根據(jù)Tennant 法，水庫壩址斷面最小生態(tài)流量為多年平均流量的10%，即0.681 m3/s。

根據(jù)最枯月平均徑流法，水庫壩址斷面最小生態(tài)流量為90%保證率最枯月平均流量，即1.43 m3/s。

3.1.2 水力學法

通過5 個斷面特征水文參數(shù)變化論證減水河段所需最小生態(tài)流量?？紤]多年平均流量條件下的壩下區(qū)間支流匯入，以壩址下泄多年平均流量的10%、15%、20%、30%、35%、50%、100% 以及90%最枯月平均流量的生態(tài)基流計算水力參數(shù)，結合水位流量關系曲線確定斷面相應的水面寬度、平均流速、平均水深、濕周等，與魚類棲息生境要求對比分析。

3.2 維持河流水環(huán)境功能的用水量確定

3.2.1 90%最枯月流量

枯水期（10 月—翌年3 月）最低下泄流量為壩址斷面90%保證率最枯月平均流量，即1.43 m3/s，豐水期（4—9 月）最低下泄流量為30%壩址多年平均流量，即為 2.04 m3/s。

3.2.2 保障下游水環(huán)境自凈需水流量

采用穩(wěn)態(tài)水質(zhì)模型，將鳳山水庫壩址至鳳山橋邊國控斷面作為整體單元，按規(guī)劃水平年2030 年區(qū)域水資源利用、水污染防治條件下的退水及污染負荷進行計算，得到鳳山水庫保障下游水環(huán)境納污能力的下泄流量為1.82 m3/s。

3.3 下泄流量綜合分析

綜合考慮生態(tài)需水、水環(huán)境自凈需水，取最大值為各時段生態(tài)流量推薦值。鳳山水庫壩下重安江干流主要斷面不同計算方法計算及推薦生態(tài)流量見表1。

表1 重安江下游關鍵斷面生態(tài)基本流量 m3/s

鳳山水庫壩址斷面枯水期（10 月—翌年3 月）最低下泄流量為1.82 m3/s，占多年平均流量的27%，在豐水期（4—9 月）最低下泄流量為30%壩址多年平均流量，即為2.04 m3/s；在4—5 月，根據(jù)下游魚類產(chǎn)卵需求進行人造脈沖過程。

4 生態(tài)調(diào)度方案

4.1 生態(tài)調(diào)度初步方案

4.1.1 典型過程選擇

根據(jù)1958—2014 年間4—5 月最大日流量資料，采用水文排頻法，選取20%頻率最大日流量（P=20%）14.45 m3/s 作為生態(tài)調(diào)度方案峰值流量。

4.1.1.1 3 d 典型過程選擇

生態(tài)調(diào)度初步計算以天然狀態(tài)下典型水文過程為主。首先篩選1958—2014 年間4、5 月份中，洪峰流量大于14.45 m3/s 的小洪水水文過程。選取1974 年 5 月 、 1960 年 5 月 、 1982 年 4 月 連 續(xù) 3 d 流量過程為典型水文過程，分別設置為典型過程1、2、3。見表2。

表2 鳳山水庫生態(tài)調(diào)度典型3 d水文過程

4.1.1.2 其他典型過程選擇

3 d 過程小洪水需水量較少，通過庫區(qū)調(diào)蓄可實現(xiàn)，主要適用于枯水年。對于平水年和豐水年來水較為充沛，可適當延長生態(tài)調(diào)度時間，或加大生態(tài)流量泄放。本次選取1982 年4 月5 日流量過程、2002 年5 月7 日流量過程為其他典型水文過程，分別設置為典型過程4、5，以研究水量充沛年份的生態(tài)調(diào)度過程。

4.1.2 計算結果

通過模擬天然小洪水過程，得出不同水文條件下研究河段流速、水深，分析是否滿足魚類產(chǎn)卵需求。

4.1.2.1 流速

典型水文過程1 各斷面3 d 流速范圍0.24～3.49 m/s，均滿足魚類產(chǎn)卵流速要求。

典型水文過程2 各斷面3 d 流速范圍0.25～3.48 m/s，均滿足魚類產(chǎn)卵流速要求。

典型水文過程3 各斷面3 d 流速范圍0.20～3.49 m/s，均滿足魚類產(chǎn)卵流速要求。

典型水文過程4 各斷面5 d 流速范圍0.17～10.19 m/s。斷面2、斷面5 部分時間不滿足魚類產(chǎn)卵流速要求。

典型水文過程5 各斷面7 d 流速范圍0.13～12.49 m/s。斷面2、斷面5 部分數(shù)據(jù)不滿足魚類產(chǎn)卵流速要求。

4.1.2.2 水深

典型水文過程1 各斷面3 d 水深范圍0.258～3.633 m，除斷面4 外其他斷面滿足魚類產(chǎn)卵水深要求。

典型水文過程2 各斷面3 d 水深范圍0.422～3.788 m，除斷面4 外其他斷面滿足魚類產(chǎn)卵水深要求。

典型水文過程3 各斷面3 d 水深范圍0.169～3.839 m，除斷面4 外其他斷面均滿足魚類產(chǎn)卵水深要求。

典型水文過程4 各斷面5 d 水深范圍0.152～3.841 m，除斷面4 外其他斷面均滿足魚類產(chǎn)卵水深要求。

典型水文過程5 各斷面7 d 水深范圍0.114～4.133 m，除斷面4 外其他斷面均滿足魚類產(chǎn)卵水深要求。

4.2 生態(tài)調(diào)度推薦方案

4.2.1 推薦水文過程

結合天然狀態(tài)典型水文過程分析結果，設計生態(tài)調(diào)度推薦水文過程。對于典型水文過程，滿足魚類習性指標的最小日流量為11.21 m3/s，以此作為生態(tài)調(diào)度水文過程的最小值；典型過程中的4 個峰值流量均滿足魚類習性指標，分別為17.63、16.72、15.48、15.27 m3/s，選取其中偏大的3 個洪峰值以此作為7、5、3 d 的洪峰流量；水文過程中其他流量則以天然水量為主。調(diào)整后的生態(tài)調(diào)度推薦過程（如圖3 所示）具體如下。

圖3 鳳山水庫生態(tài)調(diào)度推薦水文過程曲線

推薦過程1：生態(tài)調(diào)度時長持續(xù)3 d，峰值流量15.48 m3/s，相應水量 327.47 萬 m3，適用于枯水年等來水量較少的情況。

推薦過程2：生態(tài)調(diào)度時長持續(xù)5 d，峰值流量16.72 m3/s，相應水量 543.28 萬 m3，主要適用于平水年等來水量一般的情況。

推薦過程3：生態(tài)調(diào)度時長持續(xù)7 d，峰值流量17.63 m3/s，相應水量 791.16 萬 m3，主要適用于豐水年等來水量較充沛的情況。

4.2.2 推薦過程1

推薦水文過程1 各斷面3 d 流速及水深計算結果見表3。各斷面流速在0.26～3.49 m/s，均滿足流速0.2 m/s 要求。各斷面水深在0.523～3.756 m，均滿足水深0.5 m 要求。見表3。

表3 推薦水文過程1各斷面3 d流速及水深計算結果

4.2.3 推薦過程2

推薦水文過程2 各斷面5 d 流速及水深計算結果見表4。各斷面流速在0.24～3.53 m/s，均滿足流速0.2 m/s 要求。各斷面水深在0.523～4.213 m，均滿足水深0.5 m 要求。

表4 推薦水文過程2各斷面5 d流速計算結果

4.2.4 推薦過程3

推薦水文過程3 各斷面7 d 流速計算結果見表5。各斷面流速在0.24～3.53 m/s，均滿足流速0.2 m/s要求。各斷面水深在0.523～4.463 m，均滿足水深0.5 m 要求。

表5 推薦水文過程3各斷面7 d流速計算結果

5 生態(tài)調(diào)度運行管理

生態(tài)調(diào)度主要是采用人造洪水脈沖的方式來營造近自然的水文過程，保障下游水生生物產(chǎn)卵期水文需求。鳳山水庫應優(yōu)先滿足生活用水和下游河道生態(tài)用水，其次滿足工業(yè)用水，多余水量用于發(fā)電。水庫全年下泄流量不低于1.82 m3/s，4—9月下泄流量不低于2.04 m3/s。

5.1 生態(tài)調(diào)度與庫容水位關系

圖4 為鳳山水庫生態(tài)優(yōu)化后的運行調(diào)度圖。當庫區(qū)水位在Ⅰ級區(qū)域時來水量較少，僅保證生態(tài)基流與生活供水，該工況可通過水庫年調(diào)節(jié)能力避免；當庫區(qū)水位在Ⅱ級區(qū)域時，滿足Ⅰ區(qū)供水前提下，可滿足3 d 生態(tài)調(diào)度用水量；當庫區(qū)水位在Ⅲ級區(qū)域時，可額外滿足發(fā)電用水需求以及一定工業(yè)用水；當庫區(qū)水位在Ⅳ級區(qū)域時，為庫區(qū)正常蓄水位，生態(tài)流量、生態(tài)調(diào)度、發(fā)電用水、工業(yè)用水均能滿足，多余水量下泄可用于加大生態(tài)用水以及其他供水，在4—5 月可根據(jù)水庫水位和來水進行5 d 和7 d 的生態(tài)調(diào)度。

圖4 鳳山水庫生態(tài)優(yōu)化后的供水調(diào)度與庫區(qū)水位關系

5.2 生態(tài)調(diào)度運行管理

本次生態(tài)調(diào)度分為3 個推薦過程，分別如下。

推薦過程1：生態(tài)調(diào)度時長持續(xù)3 d，峰值流量15.48 m3/s，相應水量為 327.47 萬 m3，主要適用于枯水年來水量較少的情況。

推薦過程2：生態(tài)調(diào)度時長持續(xù)5 d，峰值流量16.72 m3/s，相應水量為 543.28 萬 m3，主要適用于平水年來水量一般的情況。

推薦過程3：生態(tài)調(diào)度時長持續(xù)7 d，峰值流量17.63 m3/s，相應水量為 791.16 萬 m3，主要適用于豐水年來水量較充沛的情況。

6 結 語

本文采用數(shù)值模擬方法探索鳳山水庫生態(tài)流量及生態(tài)調(diào)度方案?；诤恿魉畡恿δＰ虷ECRAS，建立鳳山水庫下游減水河段模型，結合壩址斷面水位流量曲線進行參數(shù)率定和模型驗證。結合天然水文過程提出生態(tài)流量調(diào)度初步方案，經(jīng)生境指標比對和方案優(yōu)化，得出鳳山水庫生態(tài)調(diào)度推薦方案。根據(jù)枯水年、平水年、豐水年來水條件，分別實施3、5、7 d 生態(tài)調(diào)度，相應峰值分別為15.48、16.72、17.63 m3/s，相應水量分別為327.47萬、543.28 萬、791.16 萬 m3。