三峡水库建成蓄水后,水库周边区间流域产汇流条件发生改变,新形成的入库洪水较天然条件下的坝址洪水具有“洪峰增大,涨水段洪量增大”等特性。本文基于多输入单输出系统模型模拟计算三峡建库前的入库流量系列,定量分析了三峡入库与坝址的年最大洪峰、3、7、15 d洪量之间的关系,整编得到三峡入库流量不连续序列,据此开展三峡入库洪水频率分析,并与坝址计算成果进行比较。结果表明,入库系列与坝址系列的洪峰以及3 d、7 d洪量差异明显,15 d洪量变化不大,入库百年一遇设计洪峰较坝址设计洪峰增大7.89%,入库百年一遇3 d、7 d、15 d洪量较坝址设计洪量分别增大5.49%、2.52%和0.31%,现状条件下采用入库洪水设计值作为三峡水库防洪调度的科学依据更为合理。 The runoff generation and confluence conditions at intervener basin have been changed after the com-pletion of the Three Gorges Reservoir (TGR). Since the reservoir inflow has higher peak with larger flood volume in the rising limb compared to the floods at dam site, the original designed flood results at dam site can no longer satisfy the current requirement of operation and management for TGR. In this paper, a multiple inputs and single output linear system model was built to calculate the historical reservoir inflow series and then linear regressions were used to determine the quantitative relations of flood peaks, 3 d, 7 d and 15 d flood volumes between inflow series and flood series at dam site. Subsequently, the simulated inflow series combined with historical floods were obtained and used for inflow frequency analysis. Comparing the inflow frequency analysis results with the existing design flood at dam site, the flood peaks and 3d, 7 d flood volumes are increased significantly, while the change of 15 d flood volume is not evident. It is shown that 100-year designed inflow peak discharge, 3 d, 7 d and 15 d inflow volumes are 7.89%, 5.49%, 2.52% and 0.31% larger than those of designed floods at dam site, respectively. Therefore, it’s more rational to take the designed inflow floods as scientific bases for TGR reservoir operation and management.
钟逸轩1,郭生练1,刘章君1,郭海晋2,李妍清2
1武汉大学水资源与水电工程科学国家重点实验室,湖北 武汉
2长江水利委员会水文局,湖北 武汉
Email: zhongyx_92@163.com
收稿日期:2015年7月21日;录用日期:2015年8月5日;发布日期:2015年8月19日
三峡水库建成蓄水后,水库周边区间流域产汇流条件发生改变,新形成的入库洪水较天然条件下的坝址洪水具有“洪峰增大,涨水段洪量增大”等特性。本文基于多输入单输出系统模型模拟计算三峡建库前的入库流量系列,定量分析了三峡入库与坝址的年最大洪峰、3、7、15 d洪量之间的关系,整编得到三峡入库流量不连续序列,据此开展三峡入库洪水频率分析,并与坝址计算成果进行比较。结果表明,入库系列与坝址系列的洪峰以及3 d、7 d洪量差异明显,15 d洪量变化不大,入库百年一遇设计洪峰较坝址设计洪峰增大7.89%,入库百年一遇3 d、7 d、15 d洪量较坝址设计洪量分别增大5.49%、2.52%和0.31%,现状条件下采用入库洪水设计值作为三峡水库防洪调度的科学依据更为合理。
关键词 :三峡水库,入库洪水,坝址洪水,频率分析,系统模型
三峡工程是长江治理开发的关键性工程,是长江综合防洪体系的骨干工程,在长江中下游防洪体系中具有重要地位。研究表明,水库建成后,库区被淹没,坝址以上形成了一个水面宽广的人工湖泊,原有的河槽调蓄能力丧失,使得原坝址洪水转变成具有传播时间缩短、洪峰出现时间提前、洪峰增高、峰形变瘦、洪水历时缩短、涨水段洪量增大等特点的入库洪水,径流系列的一致性遭到破坏[
为获得足够长的入库流量系列以开展入库设计洪水研究,有必要推求三峡水库建库前的入库洪水序列。已有的入库洪水计算方法大致可分为以下四种:1) 流量叠加法。分别计算干支流入库断面流量过程与区间洪水过程,并将二者同时刻叠加作为入库流量过程。2) 坝址流量反演法。采用马斯京根法、槽蓄曲线法或汇流曲线反演法中的一种或多种,将宜昌站实测流量反演至入库点,即可得到入库流量过程。3) 水量平衡法。根据水库的泄流及蓄变量变化过程反推入库洪水过程,可分为考虑静库容或动库容两种情况。4) 相应关系法。通过已有入库洪水资料和实际来水资料,建立坝址洪峰与入库洪峰及坝址流量与入库流量之间的经验关系。然而,以上方法存在一些不足之处。流量叠加法需要获得未控区间的入流过程,但实际上区间过程是一个较为模糊的概念,无法直接观测得到,区间过程计算的正确与否也直接影响入库过程的准确性。坝址流量反演法若采用马斯京根反演法,可能出现锯齿状震荡的错误结果。水量平衡法则需要大量区间沿程水位过程数据,无法应用于无水位资料地区。相应关系法虽简单易行,却无法获得完整的入库洪水过程,只能得到一系列特征值,不利于水库调度及防汛工作的开展。
为准确获得三峡水库建库前的入库流量过程,增加已有数据的代表性,本文基于多输入单输出系统模型模拟得到三峡水库1960~2002年(共43年)的汛期入库洪水过程;建立入库系列与坝址系列的定量关系,插补得到三峡8场历史调查入库洪水与1877~1959年入库流量过程;根据得到的不连序入库洪水系列,开展入库洪水频率分析,并与坝址设计成果对比。
三峡水库位于湖北省宜昌市三斗坪镇,水库控制流域面积为100万km2,坝址多年平均流量为14,300 m3/s,多年平均径流量4510亿m3。水库长度为570~650 km,库面平均宽度约为1.1 km,属河道型水库。三峡水利枢纽是具有防汛、抗旱、发电、航运等多项效益的大型水利水电工程,汛期为每年的6~9月。三峡水库区间呈东西向长条状分布,集水面积为55,907 km2。区间流域的水文站和雨量站分布情况如图1所示,由图可知,三峡入库流量主要来源于三部分:上游寸滩、武隆控制站来水和水库区间入流。
研究数据来源于长江水利委员会水文局,包括三峡寸滩站、武隆站1960~2013年6 h流量,宜昌站1877~2002年6 h流量,三峡区间71个雨量站1960~2013年6 h降雨量,三峡水库2003~2013年动库容法计算得到的6 h入库过程。
三峡区间形状狭长,不宜直接采用面雨量作为输入。为考虑降水空间分布的不均匀性与不同地点河网汇流的影响[
多输入单输出模型系统模型在洪水演进问题上应用效果良好[
式中:
为获得三峡建库前的入库流量,本文基于多输入单输出模型建立了三峡入库洪水模拟模型。模型输入量为三峡上游寸滩、武隆站流量以及区间六个子分区的净雨过程。其中寸滩、武隆站已有实测6 h流量数据,可直接作为输入量;采用降雨径流关系法将三峡区间雨量站网1960~2013年6 h数据计算得到六个子分区的净雨过程。模型的计算形式如下:
式中:
若以矩阵形式表示,则上式可简化为:
其中:
图1. 三峡区间流域站点分布示意图
区间名称 | 集水面积(km2) | 站名 |
---|---|---|
寸~清(左) | 8449 | 寸滩、坛同、大洪河、长寿(二)、狮子滩、高硐、梅市 |
寸~清(右) | 7819 | 太平、沿塘、清溪场(三)、武隆 |
清~万 | 9023 | 两河、忠县、丰都、石柱、桥头、黄水、沙子 |
万~奉 | 12828 | 合兴、余家、南门、南雅、临江、正坝、温泉、大进、岩水、沙沱、盐渠、万县、临溪、走马、建南、谋道、龙驹、龙角、云阳 |
奉~巴 | 9731 | 奉节(二)、明水、龙门、塘坊、红池坝、建楼、高楼、中良、土城、万古、巫溪、福田、巫山、大昌、长安、西宁、徐家坝、双阳、寻乐坪、堆子、庙宇、吐祥 |
巴~宜 | 8057 | 巴东、秭归、兴山、九冲、郑家坪、张官店、水月寺、黄陵庙、杨林桥、石板坪、三斗坪、宜昌 |
表1. 三峡区间流域雨量站分区表
由模型的矩阵形式容易看出,求解模型参数的关键在于求解脉冲响应函数值
三峡水库目前具有2003~2013年共11年采用动库容槽蓄曲线法计算得到的汛期入库流量,取2003~2010年作为模型率定期,2011~2013年为模型检验期,并选用Nash-Sutcliffe效率系数
式中:
模型率定期与检验期的性能指标见表2。其中,率定期与检验期
入库过程实际上无固定入流点,通常可认为集中入流处位于水库回水末端,本文假定三峡水库入库点位于清溪场水文站,并提出以下检验模型模拟结果合理性的方法。采用马斯京根法将模拟得到的入库流量系列演算至宜昌,并与宜昌实测坝址流量系列作对比。由图2中部分展示的典型年对比结果可知,模拟入库流量演算至宜昌后,与实测坝址过程线基本重合,其
入库洪水与坝址洪水之间的关系常通过建库前后的峰量关系体现。以坝址系列为横坐标、入库系列为纵坐标,将两系列1960~2013年的年最大日洪峰、年最大3、7、15 d洪量绘制于散点图上(图3),由散点分布可初步确定坝址与入库系列之间符合线性关系。
对图3各子图中的散点进行线性回归,计算得到入库与坝址年最大日洪峰、年最大3、7、15 d洪量的比值依次为1.11、1.08、1.02与1.01,相关系数ρ分别为0.9623,0.9657,0.9747与0.9798,线性关系显著。随着变量统计时段增大,直线的相关系数不断增大,即入库系列与坝址系列之间相关性也随之增大,且回归直线的系数逐渐减小并趋近于1,符合水量平衡原理以及入库洪水“涨水段洪量增大”的特性。
三峡设计规划时为增加资料的代表性,除采用宜昌站坝址系列外,还用到了经现场调查与文献考证结合的方法得到的8场历史洪水资料,由于长时段洪量难以确定,15 d洪量仅采用1788、1860和1970年的历史调查成果[
通过以上工作,可将三峡水库的入库流量系列由原来的2003~2013年共11年,延长为1877~2013年共137年与8年历史洪水入库特征值组成的不连序入库系列,其中年最大洪峰流量与年最大3、7 d洪量在调查考证期内最大值出现年份为1870年,可调查最远年份为1153年,距今约860余年;年最大15 d洪量可调查最远年份为1788年,距今约220余年。可见,对三峡入库洪水系列进行模拟计算后,极大地延长了系列长度,增加了资料的代表性,使入库系列满足开展入库洪水频率分析的数据要求。
参考三峡原计算成果,采用不连序序列统一处理法,计算各年洪峰值对应的经验频率[
式中:
实测洪水经验频率为:
式中:
参考三峡初设成果[
时期 | NSE | RE | QE |
---|---|---|---|
率定期(2003~2010) | 97.88% | 0.13% | 2.86% |
检验期(2011~2013) | 98.58% | −0.74% | 1.30% |
表2. 多输入单输出系统模型模拟效果
历史年份 | 年最大洪峰(m3/s) | 年最大3 d洪量(亿m3) | 年最大7 d洪量(亿m3) | 年最大15 d洪量(亿m3) | ||||
---|---|---|---|---|---|---|---|---|
坝址 | 入库 | 坝址 | 入库 | 坝址 | 入库 | 坝址 | 入库 | |
1870 | 105,000 | 116,600 | 265 | 286 | 537 | 548 | 975 | 985 |
1860 | 92,500 | 102,700 | 232 | 251 | 474 | 483 | 828 | 836 |
1788 | 86,000 | 95,500 | 216 | 233 | 442 | 451 | 785 | 793 |
1153 | 92,800 | 103,000 | 233 | 251 | 475 | 485 | — | — |
1227 | 96,300 | 106,900 | 242 | 261 | 493 | 502 | — | — |
1560 | 93,600 | 103,900 | 235 | 254 | 479 | 489 | — | — |
1796 | 82,200 | 91,200 | 206 | 222 | 423 | 432 | — | — |
1613 | 81,000 | 89,900 | 203 | 219 | 417 | 426 | — | — |
表3. 宜昌站历史洪水信息
统计量 | 统计参数 | 设计值 | ||||||
---|---|---|---|---|---|---|---|---|
Ex | Cv | Cs/Cv | 0.01% | 0.10% | 1% | 5% | ||
年最大洪峰 | 坝址 | 52,000 | 0.21 | 4 | 113,000 | 98,800 | 83,700 | 72,300 |
入库 | 54,800 | 0.22 | 4 | 123,300 | 107,400 | 90,300 | 77,200 | |
变幅 | 5.38% | — | — | 9.12% | 8.70% | 7.89% | 6.78% | |
年最大3 d洪量 | 坝址 | 130 | 0.21 | 4 | 282.1 | 247 | 209.3 | 180.7 |
入库 | 134 | 0.22 | 4 | 301.6 | 262.5 | 220.8 | 188.7 | |
变幅 | 3.08% | — | — | 6.91% | 6.28% | 5.49% | 4.43% | |
年最大7 d洪量 | 坝址 | 275 | 0.19 | 3.5 | 547.2 | 486.8 | 420.8 | 368.5 |
入库 | 270 | 0.21 | 3.5 | 573.5 | 505.4 | 431.4 | 373.5 | |
变幅 | −1.82% | — | — | 4.81% | 3.82% | 2.52% | 1.36% | |
年最大15 d洪量 | 坝址 | 524 | 0.19 | 3 | 1022 | 911.8 | 796.3 | 702.2 |
入库 | 515 | 0.20 | 3 | 1034.9 | 922.4 | 798.8 | 700.1 | |
变幅 | −1.72% | — | — | 1.26% | 1.16% | 0.31% | −0.30% |
表4. 三峡洪水频率分析结果日均洪峰:m3/s洪量:亿m3
*注:变幅 = (入库值 − 坝址值)/坝址值。
散程度[
对P = 0.01%、0.1%、1%和5%的入库设计洪水与坝址设计洪水进行比较,结果表明,除年最大15 d洪量的5%设计值略小于坝址设计值外,入库设计值均大于坝址设计值。以洪峰为例,三峡水库各重现期的入库设计洪峰较坝址设计值增大6.78%~9.12%,其中0.01%入库设计洪峰较坝址设计值增大10,300 m3/s,增幅高达9.12%,
图2. 1981和1982年宜昌站汛期演算流量与坝址流量对比
图3. 三峡入库系列与坝址系列相关图
图4. 三峡入库洪水频率曲线
5%设计洪峰较坝址设计值增大4900 m3/s,增幅为6.78%,入库年最大洪峰设计值增加显著。各重现期年最大入库3 d洪量增幅为4.43%~6.91%,年最大7 d洪量增幅为1.36%~4.81%,洪量增幅均较为明显。年最大15 d洪量变化不大,变幅位于−0.30%~1.26%之间。
由以上分析可知,三峡建库后,入库洪水分布相比原坝址洪水均有着不同程度的改变,其中年最大洪峰与年最大3 d、7 d洪量的设计值变化显著,年最大15 d洪量分布无明显变化。因此,采用入库洪水系列推求三峡建库后的设计洪水,更符合建库后的实际情况,并保证水库的安全高效运行,也更加科学合理。
本文基于多输入单输出模型计算得到三峡建库前的入库流量,并检验了模拟结果的合理性与可靠性;定量统计了入库流量系列与坝址流量系列之间的经验关系,据此得到宜昌站历史坝址洪水对应的入库特征值;对得到的不连序入库系列进行频率分析,将结果与三峡坝址设计成果进行对比。得到的主要结论如下:
1) 三峡入库MISO模型模拟效果良好,率定期与检验期模型效率系数分别达到97.88%与98.58%,水量平衡系数分别为0.13%与−0.74%,洪峰相对误差分别为2.86%与1.30%,模型精度达到水文情报预报规范的甲等精度。采用马斯京根法将模拟入库流量演算至宜昌站并与实测坝址洪水流量对比,二者高度相符,演算系列的确定性系数达93.73%,模型模拟结果很高。
2) 三峡入库洪水与坝址洪水统计量之间存在显著的线性关系,相关系数均高于0.96。据此延长三峡入库洪水资料,得到三峡1877~2013年入库流量系列与8年历史入库特征值组成的入库洪水不连序系列。
3) 对得到的三峡入库洪水不连序系列进行频率分析,并与原三峡坝址设计成果对比。结果表明,建库前后年最大15 d洪量分布变化不大,年最大洪峰与年最大3 d、7 d洪量分布则发生了显著改变,因此采用入库洪水系列推求三峡建库后的设计洪水,更符合建库后的实际情况,并更好地为三峡水库防洪调度提供科学依据。
钟逸轩,郭生练,刘章君,郭海晋,李妍清, (2015) 三峡水库入库洪水分析计算研究Analysis of Inflow Floods of the Three Gorges Reservoir. 水资源研究,04,330-338. doi: 10.12677/JWRR.2015.44040