李彬瑜，薛萬(wàn)來(lái)，姜新佩，史紅偉

(1.河北工程大學(xué)，河北 邯鄲 056107；2.北京市水科學(xué)技術(shù)研究院，北京 100048；3.中國(guó)水務(wù)投資有限公司，北京 100053)

0 引 言

河流生態(tài)基流是維護(hù)河湖等生態(tài)系統(tǒng)功能不喪失、需要保留的底限流量過(guò)程中的最小值[1]。合理的生態(tài)基流可以最大限度地平衡水生態(tài)系統(tǒng)與水資源利用之間的矛盾。因此，適合的生態(tài)基流計(jì)算方法則顯得至關(guān)重要。目前，國(guó)內(nèi)外有關(guān)生態(tài)基流的計(jì)算方法已近百種，大體可分為水文學(xué)法、水力學(xué)法、棲息地法、整體分析法四大類[2]。由于我國(guó)生態(tài)資料的缺乏致使棲息地法和整體分析法計(jì)算有一定困難[3]；水文學(xué)法和水力學(xué)法則比較常用[4-5]，水力學(xué)法對(duì)河道斷面資料需求較高則操作性較差；而數(shù)據(jù)需求量少、成本較低的水文學(xué)法成為大多數(shù)學(xué)者的首選。張澤聰?shù)萚6]引入季節(jié)修正系數(shù)并利用中位數(shù)代替平均數(shù)改進(jìn)Tennant法計(jì)算了大凌河生態(tài)基流量；葛金金等[7]選取Tennant法、Q90-Q50法、Lyon法、Tessman法、月流量變動(dòng)法評(píng)估了沙穎河周口段的適宜生態(tài)基流，并綜合水量需求、水文節(jié)律和滿足率認(rèn)為L(zhǎng)yon法最為合理；郭利丹等[8]對(duì)比分析了最小月平均流量法、逐月頻率計(jì)算法、最小連續(xù)平均流量法、逐月最小生態(tài)徑流計(jì)算法，認(rèn)為計(jì)算河流生態(tài)基流的最優(yōu)方法為逐月最小生態(tài)徑流法。本文針對(duì)永定河這類北方季節(jié)性河流，以其官?gòu)d下游段為例采用Tennant法、90%保證率最枯月平均流量法、基流比例法，估算生態(tài)基流，旨在為永定河官?gòu)d下游段生態(tài)需水量調(diào)控，提供科學(xué)依據(jù)。

1 流域概況

永定河地跨山西、內(nèi)蒙古、河北、北京、天津等，流域面積約4.7萬(wàn)km2，上游有桑干河、洋河兩大支流，于河北省張家口懷來(lái)縣朱官屯匯合后稱永定河，在官?gòu)d水庫(kù)納媯水河，從官?gòu)d水庫(kù)起穿越八達(dá)嶺高原形成了官?gòu)d山峽，至三家店流入華北平原，三家店為永定河流域山區(qū)、平原分界，其中山區(qū)流域面積4.51萬(wàn)km2，占95.8%，平原流域面積1 953 km2，占4.2%。永定河流域?qū)贉貛Т箨懶约撅L(fēng)氣候，為半濕潤(rùn)、半干旱型氣候過(guò)渡區(qū)多年平均降水量在360～650 mm之間，不同地區(qū)降水量差異頗大，多雨區(qū)和少雨區(qū)相差將近1倍。

2 生態(tài)基流計(jì)算

依據(jù)官?gòu)d水文站1956年～2016年61a的實(shí)測(cè)月徑流資料，分別用Tennant法、90%保證率最枯月平均流量法、基流比例法分析計(jì)算永定河官?gòu)d下游段的生態(tài)基流。

2.1 豐、平、枯水年劃分

文獻(xiàn)[9]將河流豐、平、枯概念劃分為特豐水年、偏豐水年、平水年、偏枯水年、特枯水年，其中特豐水年、偏豐水年統(tǒng)稱為豐水年；特枯水年、偏枯水年統(tǒng)稱為枯水年。根據(jù)官?gòu)d站1956年～2016年61 a的年徑流頻率曲線，在一定保證率(P)的年徑流量劃分年型見(jiàn)表1。

表1 官?gòu)d站年徑流豐、平、枯水年劃分標(biāo)準(zhǔn)

2.2 Tennant法

Tennant法：認(rèn)為年平均流量的百分比即生態(tài)基流[10]。該方法推薦的生態(tài)基流分為汛期與非汛期，并以占平均流量的比例作為劃分標(biāo)準(zhǔn)，由官?gòu)d站1956年～2016年61 a水文資料，利用Tennant法計(jì)算生態(tài)基流(見(jiàn)表2)。汛期和非汛期分別取對(duì)應(yīng)時(shí)段年均流量的30%(8.21 m3/s)和10%(1.41 m3/s)作為生態(tài)基流。

表2 Tennant方法計(jì)算標(biāo)準(zhǔn)與結(jié)果

2.3 90%保證率最枯月平均流量計(jì)算法

此方法根據(jù)官?gòu)d站1956年～2016年61 a的月流量資料，建立一系列月流量時(shí)段曲線，在90%保證率下的逐月最小生態(tài)需水量[11]計(jì)算結(jié)果如表3所示。最小生態(tài)基流范圍1.54～3.24 m3/s，占多年平均流量20.73 m3/s的7.43%～15.63%。最小生態(tài)基流量、最大生態(tài)基流均分別出現(xiàn)在枯水期的2月份與豐水期7月份。

表3 官?gòu)d站下游90%保證率下最小生態(tài)需水量

2.4 基流比例法

2.4.1 計(jì)算方法

以豐水年作為計(jì)算起點(diǎn)[12]，按各年型平均流量之間的比值關(guān)系，確定基流比例。即

Ti+1=αTi

(1)

式中，Ti為監(jiān)測(cè)斷面第i年型的基流比例，%；i=1、2、3、4、5，分別代表特豐水年、偏豐水年、平水年、偏枯水年、特枯水年；α為比例系數(shù)；并且

α=1+(Qi/Qi+1-1)×μ

(2)

式中，Qi為監(jiān)測(cè)斷面第i種年型的平均流量，m3/s；μ取0.4[13]。

各時(shí)段基流比例與相應(yīng)平均流量相乘即為生態(tài)基流量

Si，j=Qi，jTi

(3)

式中，Sij和Qij分別為斷面第i時(shí)段、第j個(gè)月的的生態(tài)基流量及流量，m3/s。

2.4.2 各年型生態(tài)基流計(jì)算

按照前文豐、平、枯水年的劃分標(biāo)準(zhǔn)，將61 a實(shí)測(cè)逐月平均流量數(shù)據(jù)劃分年型。按Tennant法所論述的河道流量同河道生態(tài)健康狀況關(guān)系，假定豐水年生態(tài)基流量占其平均徑流的20%[14]。由前文所述方法計(jì)算各年型的生態(tài)基流及基流比例，結(jié)果見(jiàn)表4。

表4 官?gòu)d下游段各年型生態(tài)基流

由表4可知，永定河流域官?gòu)d下游段特豐水年至特枯水年的生態(tài)基流量為10.9 ～2.05 m3/s，占多年平均流量20.73 m3/s的9.86%～52.44%，較為符合Tennant法所推薦的比例。由多年生態(tài)基流與實(shí)測(cè)多年徑流作圖比較(見(jiàn)圖1)，可得官?gòu)d站實(shí)測(cè)年徑流在2006年之前均能滿足生態(tài)基流需求，但在2006年至2016年已經(jīng)不能滿足生態(tài)基流需求，反映出自2006年起永定河流域官?gòu)d下游段生態(tài)系統(tǒng)急劇退化。

圖1 1956年～2016年生態(tài)基流分析

2.4.3 逐月生態(tài)基流分析

為能夠更好地了解河道生態(tài)基流在年內(nèi)的變化情況，依據(jù)式(3)統(tǒng)計(jì)計(jì)算了1956年～2016年61 a各典型年逐月生態(tài)基流量以及實(shí)測(cè)流量的多年平均值見(jiàn)表5。

從表5可以得出，豐水年生態(tài)基流在5.13～19.18 m3/s之間，占多年平均流量20.73 m3/s的24.75%～92.52%；偏豐水年生態(tài)基流在3.56～14.95 m3/s之間，占多年平均流量的17.17%～72.12%；平水年生態(tài)基流在3.41～10.20 m3/s之間，占多年平均流量的16.45%～49.20%；偏枯水年生態(tài)基流在3.35～5.78 m3/s之間，占多年平均流量的16.16%～27.88%；特枯水年生態(tài)基流在1.25～2.92 m3/s之間，占多年平均流量的6.03%～14.09%；其中，豐水年、偏豐水年、平水年、偏枯水年的月生態(tài)基流極大值均出現(xiàn)在6月份，這與多年平均下的極大值出現(xiàn)月份一致，特枯水年此極大值提前到了5月份。豐水年、偏豐水年、平水年生、枯水年的月生態(tài)基流極小值均出現(xiàn)在2月份，這與多年平均下的極小值出現(xiàn)月份也是一致的，偏枯水年此極小值則出現(xiàn)在了豐水期的7月份。

表5 官?gòu)d下游段各年型生態(tài)基流 m3/s

3 生態(tài)基流計(jì)算結(jié)果比較與現(xiàn)狀年保障程度分析

3.1 結(jié)果比較

經(jīng)上述3種方法計(jì)算的生態(tài)基流結(jié)果存有一定差異見(jiàn)表6。由表6可看出，采用90%保證率最枯月平均流量法的計(jì)算結(jié)果占多年平均流量的2.85%～11.48%，計(jì)算結(jié)果較低?；鞅壤梢暂^好地涵蓋Tennant法與90%保證率最枯月平均流量法的計(jì)算結(jié)果。由生態(tài)基流過(guò)程線(見(jiàn)圖2)可看出，基流比例法、Tennant法均能較好地反映永定河官?gòu)d下游段生態(tài)基流年內(nèi)豐水期與枯水期的變化，但基流比例法能更好地體現(xiàn)生態(tài)基流年際間變化。考慮到永定河官?gòu)d下游段目前生態(tài)基流缺失嚴(yán)重，河道下游已出現(xiàn)斷流情況，而且生態(tài)基流是滿足生態(tài)功能需求的最小值。因此，暫時(shí)采用基流比例法在偏枯水期和特枯水期的計(jì)算結(jié)果作為生態(tài)基流取值范圍。即：2.05～3.93 m3/s，占多年平均徑流量的比例為9.86%～18.91%(見(jiàn)表6)，既能最大限度地平衡水生態(tài)系統(tǒng)與水資源利用之間的矛盾，也符合Tennant法中一般用水期的推薦比例。

表6 官?gòu)d下游段各年型生態(tài)基流

圖2 不同計(jì)算方法的生態(tài)基流過(guò)程線

3.2 現(xiàn)狀年生態(tài)基流保障程度分析

近年來(lái)，永定河官?gòu)d下游段徑流量急劇減少，考慮實(shí)際徑流是否滿足生態(tài)需水尤為重要；因此，利用2016年官?gòu)d站的逐日徑流數(shù)據(jù)與生態(tài)基流量進(jìn)行比較，對(duì)永定河官?gòu)d下游段的生態(tài)基流量進(jìn)行評(píng)估，現(xiàn)狀年永定河官?gòu)d下游段平均流量3.13 m3/s，年徑流量0.97億m3，屬于特枯水年。從表7可得，90%保證率最枯月平均流量法全年保障天數(shù)最多，為211 d，保障率57.81%；基流比例法(特枯水年)全年保障天數(shù)179 d，保障率49.04%，Tennant法全年保障天數(shù)最少，為89 d，保障率24.38%。由3種計(jì)算方法可知，最高保障率未超過(guò)60%，可見(jiàn)永定河官?gòu)d下游段生態(tài)環(huán)境狀況不容樂(lè)觀。從年內(nèi)月情況來(lái)看上述3種計(jì)算方法得出的生態(tài)基流量在現(xiàn)狀年1月～4月份保障天數(shù)為0，12月份保障天數(shù)也只有1 d，對(duì)永定河這類北方季節(jié)性河流在年內(nèi)枯水期要使得生態(tài)基流得到保障難度很大；對(duì)于豐水期而言，由于在特枯水年天然來(lái)水相對(duì)較少且人工調(diào)蓄等原因使得Tennant法計(jì)算得到的保證率相對(duì)較低，也說(shuō)明Tennant法在計(jì)算永定河此類型的北方季節(jié)性河流且上游有較多調(diào)蓄作用的水利工程時(shí)已有明顯的不足。但90%保證率最枯月平均流量法、基流比例法(特枯水年)在6月～9月份甚至10月、11月份計(jì)算得到的生態(tài)基流保障程度都超過(guò)90%甚至到100%。

表7 官?gòu)d下游段各年型生態(tài)基流

4 結(jié) 論

本文采用Tennant法、90%保證率最枯月平均流量法、基流比例法對(duì)永定河官?gòu)d下游段生態(tài)基流進(jìn)行了計(jì)算。討論分析可得出：基流比例法在該區(qū)域表現(xiàn)最好，而90%保證率最枯月平均流量法計(jì)算得出的結(jié)果偏小，Tennant法在計(jì)算永定河此類型的北方季節(jié)性河流且上游有較多調(diào)蓄作用的水利工程時(shí)已有明顯的不足。對(duì)比分析后認(rèn)為，永定河官?gòu)d下游段生態(tài)基流為2.05～3.93 m3/s；通過(guò)現(xiàn)狀年生態(tài)基流保障程度分析可以得出，對(duì)于永定河此類型的北方季節(jié)性河流來(lái)說(shuō)，上述3種方法還不能很好地平衡水生態(tài)系統(tǒng)與水資源利用之間的矛盾。特別是在枯水期，生態(tài)基流保障程度很難得到滿足。