1. 研究区概况
潇水为湘江的主源,发源于湖南省蓝山县子良乡九嶷山脉的野狗岭南麓,蓝山县往下游至永州萍岛河段一般称为潇水,永州萍岛以下称为湘江。潇水干流全长354 km,流域面积12,099 km2,其中湖南境内11,924 km2,占总流域面积的98.55%,平均坡降0.76‰。2022年,潇水流域遭遇10次强降雨、8次流域性洪水过程,江永县永明河、江华县冯河发生超历史洪水过程。涔天河水库出现建库以来最大入库流量4410 m3/s,流域内部分地区遭受了严重洪涝灾害损失。潇水流域位于湘江的上游河段,潇水流域洪水不仅直接影响潇水沿岸人民的生活,同时对湘江流域的防洪影响巨大。研究潇水流域洪水,可为潇水流域水库群防洪联合优化调度提供可靠的数据支持,对减少洪涝灾害,避免潇水洪水与湘江干流洪水遭遇,有着重要现实意义 [1] 。
2. 洪水特性及洪峰P-III型频率曲线
2.1. 洪水特性分析
潇水流域洪水主要由暴雨产生,流域内雨量较充沛,洪水主要来源于暴雨。潇水流域洪水主要发生在4~9月,根据对流域内各水文站建站至今的实测资料统计分析,年最大洪水多发生在4~7月,占全年的85%;以5、6月发生次数最多,5、6月份出现年最大洪水的几率超过55% [2] 。
自五十年代以来,有关设计单位先后多次在潇水流域进行历史洪水调查与考证,调查到的主要洪水年份有:1974年、1877年、1906年、1914年、1915年、1934年。涔天河以上流域以1877年洪水为首大,1915年洪水为次大,1934年洪水位第三位。涔天河以下则以1915年洪水为首大,次大洪水沿程有变化,其中双牌以上1914年为次大,而以下则以1974年为次大,1906年为第三位。
2.2. 汛期及其划分
根据已有的潇水流域的实测降雨和洪水资料,将每年4月到9月划分为汛期,10月到次年3月为非汛期,参考相关研究和潇水流域所处的地理位置,划分汛期与非汛期进一步精准。根据对流域内各水文站建站至今的实测资料统计分析,年最大洪水多发生在4月至7月,占全年的85%;以5、6月发生次数最多。
以图1所示老埠头站部分年份逐日流量为例:2017、2015、2013年来水较大,多数汛期洪水发生在4月份之后,8月份之前,个别年份最大洪水发生在10月之后。
2.3. 洪峰流量配线
部分站点已有洪峰流量设计成果,据此可直接进行P-Ⅲ型曲线 [3] 参数求取,已有设计洪水成果见表1和表2。
Figure 1. Daily flow discharge hydrograph of Laobutou Station from 2013 to 2018
图1. 老埠头站2013~2018年逐日流量过程线
Table 1. Existing design floods results at the Xiaoshui River (Unit: m3/s)
表1. 潇水已有设计洪水成果(单位:m3/s)
Table 2. Each station controlled catchment area in the Xiaoshui River (Unit: km2)
表2. 潇水各站控制集水面积(单位:km2)
无设计成果的站点按照实测径流数据配线。年最大洪峰系列的适线结果如图2~6所示。
3. 设计洪水
3.1. 单站设计洪水
参照《涔天河水库扩建工程任务和规模》,选取1968年6月实测洪水为全潇水流域典型洪水。涔天河入库、双牌入库、老埠头站、道县、江华、大路铺等站点有部分设计洪峰设成果,其余站点无设计成果。对于有设计成果的站点,不同频率的洪峰洪量设计值采用设计成果,对于无设计成果的站点,综合考虑控制站集水面积关系与当地实测资料排频配线成果,进行设计洪峰与设计洪量的推求。
对老埠头、南津渡、双牌入库、道县、江华县、江华站、大路铺七个控制站点的七个频率的单站设计洪水,共49场次,对应频率分别为100年一遇、50年一遇、20年一遇、10年一遇、5年一遇、3年一遇、2年一遇。
对于无设计成果的站点,采用临近有设计成果的站点设计值,根据测站控制面积进行放大或者缩小,获得
Figure 2. P-III distribution curve of annual maximum flood peak Figure 3. P-III distribution curve of annual maximum inflow
at Laobutou Station peak at Shuanpai Reservoir
图2. 老埠头站年最大洪峰P-Ⅲ型曲线 图3. 双牌入库年最大洪峰P-Ⅲ型曲线配线
Figure 4. P-III distribution curve of annual maximum flood peak Figure 5. P-III distribution curve of annual maximum flood peak
at Daoxian Station at Jianghua Station
图4. 道县站年最大洪峰P-Ⅲ型曲线配线 图5. 江华站年最大洪峰P-Ⅲ型曲线配线
Figure 6. P-III distribution curve of annual maximum flood peak at Dalupu Station
图6. 大路铺站年最大洪峰P-Ⅲ型曲线配线
Figure 7. Design flood hydrograph of Nanjindu Station
图7. 南津渡站设计洪水过程线
此类站点的洪峰设计值。设计考虑山区性河流,因此本项目洪水放大统一采用洪峰流量同倍比放大进行。单站设计洪水成果以南津渡站为例展示,设计洪水过程线如表3和图7所示。
设计成果汇总:
Table 3. Summary of design flood results for Nanjindu Station
表3. 南津渡站设计洪水成果汇总
3.2. 区间典型洪水与设计洪水
利用搬家指数法 [4] ,可求得两控制站之间的区间洪峰或者洪量情况。搬家指数法如下定义:搬家指数np即两个测站(控制断面)河道特征值对数之差与控制面积对数之差的比值。具体为:
1测站和2测站之间的区间在1%频率下洪峰流量Q的搬家指数:
(3.1)
1测站和2测站之间的区间在5%频率下最大3 d洪量W3d的搬家指数:
(3.2)
式中:n为搬家指数,下标明确了上下测站、设计频率以及设计指标,F为不同测站的控制集水面积。
利用搬家指数,可实现上下游两站的设计洪水推求区间设计洪水,即可以计算缺少径流资料系列支流区间的设计洪水(假设1、2测站之间的区间来水包括支流来水加区间来水)。
区间或者无设计资料测站设计洪峰的计算按照如下公式:
Figure 8. Typical flood hydrograph of Xiangjiang River
图8. 湘江典型洪水过程线
(3.3)
式中:Q1%为待求设计洪峰,F为待求区间或测站的控制集水面积。
根据以上方法求得各个区间不同频率的设计洪峰,同样考虑山区型河流特性,将各个区间典型洪水过程按照洪峰同倍比放大得到各个频率的区间设计洪水过程 [5] 。各个区间的典型洪水过程则同样使用规定的1968年6月实测洪水作为典型过程。区间典型洪水过程以南津渡–老埠头区间(湘江)区间为例展示,典型洪水如图8所示。
3.3. 流域设计洪水
技术路线:按照设计洪水地区组成的思路,提出流域设计洪水组合方案–推求各区间设计洪水–洪水组合方案比较–得出特定需求下的流域设计洪水 [6] 。
同时考虑河流以山区河流为主,且洪水多暴雨成因,设计中进行如下假设:单站设计洪水时段洪量符合相应频率;设计洪水的历时内包括上游涨水与下游落水的全过程;某段区间来水和此区间内支流来水概化为同一个入口汇入。潇水流域概化示意图见图9。
1) 组合洪水分析是以相邻两控制点与其区间构成的整体为单位进行的,分析时,以下游控制断面的某一设计频率的设计洪水为基础,得到甲乙两种组合洪水。甲种洪水:上下游发生设计频率洪水、区间发生相应洪水,以上游来水为主;乙种洪水:下游与区间发生设计频率洪水,上游发生相应洪水,以区间来水为主。自上游向下游进行上述分析,总共形成25 = 32种洪水组合方案。
2) 根据每种方案下区间洪水的洪量,对区间典型洪水或不同频率设计洪水按照总洪量同倍比缩放,由此得到潇水流域以下四个区域的区间来水:江华县–道县区间、道县–双牌区间、双牌–南津渡区间、南津渡–老埠头区间。这些区间设计洪水成果,应用于步骤3中,当发生区间来水为主时,区间设计频率洪水过程的推求。
3) 结合步骤1)和2),得出不同频率下不同洪水遭遇组合方案的源头与5个区间来水情况。为了模拟效率和实际需求,还要人为对各个洪水组合方案进行筛选,剔除无意义方案,突出对于某些防洪保护对象极为不利的方案。由于当洪水以支流来水为主时,干流的水库调蓄作用被缩小,相应的为了保护防洪对象,对水库的调控能力需要提出更高的要求,因此选择支流来水为主的洪水模式相较于上游干流来水的模式,更具备极端性,是风险更高的洪水情景。
因此按照“下游与区间发生设计频率洪水,上游发生相应洪水”为不利的原则,筛选出具有代表性的对重点保护对象不利的组合洪水方案。
Figure 9. Sketch map of the Xiaoshui Basin
图9. 潇水流域概化示意图
最终选择了五种情景,分别为四种支流来水为主的情况,以及一种上游河源为主的情况,共计五种洪水组合方案,作为案例,进行流域洪水设计。
以1号洪水为具体案例,同时比较1至5号洪水设计成果,具体阐述设计步骤与设计成果:
1) 首先确定潇水流域全流域本场次设计洪水的设计频率为2年一遇。
2) 根据上一章节单站设计洪水成果中的最下游控制站——老埠头站的单站设计洪水成果,确定全流域发生2年一遇洪水时流域出口控制站总洪量。
3) 老埠头主要洪水来源选取确定为上游干流,则南津渡控制断面洪量为2年一遇频率设计洪水洪量,根据水量平衡得到湘江及其区间相应洪量。利用搬家指数方法求得湘江及其区间的各频率设计洪水后,依据湘江的相应洪量对设计成果进行总洪量的同倍比缩放。
4) 根据上一步确定的南津渡洪水洪量,将南津渡洪水主要来源确定为上游干流,那么双牌控制断面洪量为2年一遇,则双牌–南津渡区间洪水洪量按照水量平衡得到。类似于湘江洪水过程的求取可得到双牌–南津渡区间洪水过程。
5) 根据上一步确定的双牌洪水洪量,将双牌洪水主要来源确定为上游干流,那么道县控制断面洪量为2年一遇,则道县–双牌区间洪水洪量按照水量平衡得到。类似的也可以求得道县–双牌区间洪水过程。
6) 根据上一步确定的道县洪水洪量,将道县洪水主要来源确定为江华–道县区间,即江华–道县区间发生2年一遇洪水,江华县即上游干流洪量为水量平衡的相应洪量。类似的,得到江华–道县区间洪水过程。
Figure 10. Design results of No. 1 flood in each tributary and section
图10. 各支流与区间1号洪水设计成果
Figure 11. Flood hydrographs for No. 1 to No. 5
图11. 1至5号洪水过程线
7) 根据上一步确定的江华县洪水洪量,将江华县洪水主要来源确定为萌渚水及其区间。这时,由于江华县洪水洪量频率不是2年一遇,萌渚水在作为洪水主源的情况下洪水频率也不会达到2年一遇。因此假设江华县为2年一遇洪量,类似于前几步可以确定在这一假设条件下萌渚水及其区间的洪水洪量于上游干流的洪水洪量,由此确定了萌渚水作为洪水主源时洪量占江华县总洪量的比率。根据江华县实际洪量与江华县2年一遇洪量的比值,进一步对萌渚水及其区间设计洪水进行缩小,得到萌渚水及其区间的洪水过程,同理得到干流源头即江华站洪水过程。
8) 在完成以上全部步骤后,得到的设计成果包括:干流源头以及流域全部支流和区间的洪水过程,即潇水流域的一场设计洪水过程。将干流源头与支流区间来水过程进行河道洪水演进,即可得到未经调度的流域出口的洪水过程。共设计流域洪水25场次,以1号洪水为例,源头及各个支流与区间来水设计成果如图10所示;1号洪水至5号洪水在未经调度调洪情况下演进至流域出口的洪水过程线如图11所示。
4. 结论与展望
本文以潇水流域为研究对象,分析了流域各测站五十多年的年最大洪峰流量系列,并基于皮尔逊III型曲线利用总偏差最小和目估配线法进行洪水频率计算,通过峰比放大和搬家指数与同频率法进行洪水设计,主要得出以下结论:
1) 潇水流域老埠头等七个测站的历史径流数据满足代表性与一致性要求,具备明显的丰平枯变化与周期性,同时在95%置信程度,满足一致性要求。潇水流域内洪水主要由暴雨产生,根据对流域内各水文站建站至今的实测资料统计分析,年最大洪水多发生在4~7月,占全年的85%;以5、6月发生次数最多,5、6月份出现年最大洪水的几率超过55%。
2) 确定流域干流主要控制断面的典型洪水过程与大部分区间和支流典型洪水过程。在已有流域各站设计成果的基础上,进行了干流7个控制站的单站设计洪水设计,补充了相关控制断面各个频率下的设计洪水成果共计49场次,同时对四个区间和支流也进行了洪水设计。
3) 在同频率流域洪水组合分析的基础上,设计了多场次流域组合洪水情景共计25场次。为潇水干流部分测站的历史径流数据整理工作提供一定参考,对于流域水情水势分析提供部分支撑,更补充了潇水流域的洪水资料数据库,一定程度上满足了流域关键节点防洪保障与全流域联合防洪规划的数据要求。
但本文洪水设计成果也存在一些不足以待改进:
流域的洪水组合情景是人为进行的主观设计和选择,对于特定需求比如防洪模拟调洪演算等适用性较高,但对于流域洪水分析存在局限性。
全流域选取同一场洪水作为典型是依据相关文件进行的,但典型过程是否符合气候变化条件下的现实情况有待考证。
潇水流域一场洪水持续时间一般不超过一周,因此洪水设计历时选择10天长度,可以保证涨水落水阶段的完整性,同时保证了总洪量在频率上的对应性和水量平衡,而小时段内则因考虑河道的蓄水作用和洪水波的演进而未进行过程线修正。
基金项目
国家重点研发计划资助(项目编号:2022YFC3800700);湖南水利科技重大项目(XSKJ2019081-08)。