无资料小流域水库设计洪水计算方法分析

水利水电快报EWRHI第41卷第5期无资料小流域水库设计洪水计算方法分析

李增国1，张楠2，张振杰1，郭玲梅1

（1.云南省水文水资源局临沧分局，云南临沧2.云南省水文水资源局玉溪分局，云南玉溪

653100）677000；

摘要：水库设计洪水成果对于水库规模的确定和泄洪设施的布置至关重要。针对无资料小流域地区传统设计洪水计算方法的结果与流域实际暴雨洪水特性出入较大的问题，提出了实测暴雨移置法，在充分分析区域暴雨洪水特性基础上，利用邻近暴雨实测资料推求水库设计洪水。以云南省临沧市某水库工程为例，使用暴雨移置法推求得该水库66a一遇洪水为96.5m3/s，与洪水调查结果吻合良好。

关键词：洪水计算；邻近暴雨实测资料；暴雨移置法；小流域；湿润地区中图法分类号：TV122.3

文献标志码：A

DOI：10.15974/j.cnki.slsdkb.2020.05.004

水库设计洪水的推求是水库勘测设计阶段的一项重要工作，关系到水库规模的确定、泄洪设施的设计布置、水库调度运行方式的选择等。大多小（1）型和中型水库工程设计时，流域内无实测暴雨洪水资料或只有短期的实测资料，无法满足计算要求。对于西南湿润无资料地区，通常使用经审批的地区暴雨统计参数等值线图集获取水库流域中心处的暴雨统计参数，由设计暴雨推求设计洪水，另外也可用推理公式法或地区经验公式法推求设计洪水。但对于部分小流域，暴雨图集法查得的暴雨统计参数与流域暴雨特性出入较大，所推求的设计洪水偏差较大。本文以云南省临沧市某水库工程为例，分别采用云南省暴雨统计参数等值线图集（以下简称“暴雨图集法”）和移置邻近实测暴雨资料（以下简称“暴雨移置法”）两种方法推求设计暴雨，通过对比分析确定水库设计暴雨成果，然后推求水库设计洪水。

上游流域位于东经99°50′～99°52′、北纬23°′～23°57′。设计流域内为无人居住区，水库库区附近有少部分耕地，面积为0.6km2，其余大部分为森林面积159km2，天然落差1877m。从河源至水库大坝下游约7km处的河段称为芒佑河，往下称白木河，至小塘坝村附近河段称塘坝河，再往下至勐段村附近河段称勐段河，最后在河底岗附近流入河底岗河。芒佑河流域水系见图1。

覆盖，覆盖率达95%。芒佑河干流全长23km，径流

2暴雨洪水特性

水库流域洪水由暴雨形成，芒佑河属河底岗河

右岸支流，流域地处哀牢山脉西侧，流域暴雨的天气系统主要为冷锋低槽型、冷锋切变型。流域位于孟加拉湾暖湿气流的迎风面，河谷深切，河流走向与水汽入流方向基本一致，故形成流域内暴雨次数多、暴雨量级较大的特点。邻近的河底岗和大兴雨量站与水库同属河底岗河流域，且距离较近，地形地貌、下垫面条件和气候条件基本一致，暴雨洪水特点相似。根据邻近河底岗雨量站39a实测资料统计，历年最大日降水量超过100mm的大暴雨有3场，其中实测最大日降水量168.5mm，发生于2014年7月21日；实测第二大日降水量为143.3mm，发

1流域概况

水库位于云南省临沧市耿马县勐永镇境内的

芒佑河上游，控制径流面积12.2km2，河长4.98km，河道平均比降128‰。芒佑河属于怒江流域南汀河水系，为南汀河左岸加入的一条二级支流。芒佑河

收稿日期：2019-10-10作者简介：李增国，男，工程师，主要从事水文测验整编、水文预测预报、水情水资源分析等方面的工作。E-mail：412808166@qq.com

·24·

李增国等无资料小流域水库设计洪水计算方法分析

图1芒估岗河流域水系

生108.8于19年9月15日；实测第三大日降水量为

日降水量为mm，发生于70mm1986。与暴雨对应，年7月25日。多年平均最大

90%的年最大洪水集中在70%。该流域具有山区河流洪水的特性，7～9月，7～8月最为突出，发生概率超过1～3h，次洪历时一般不超过

洪水陡涨24陡落，h，多为单峰型洪水过程。

起涨历时一般在3

设计暴雨推求

3.1

云南省暴雨图集法

根据SL44-2006《水利水电工程设计洪水计算

规范》［1］规定，对于无资料地区，可从通过审批的暴

雨统计参数等值线图上查算工程所需历时的设计点暴雨量［2-3］2000年，于2007。通过云南省暴雨图集［4］年通过云南省水利厅审定（资料截止库流域暴雨统计参数，其中，最大1，6，24h）暴雨均查得水值分别为35，48mm和77mm，Cv值分别为0.43，0.49

和0.48，Cs=3.5Cv。

根据水库位置的暴雨区划，由云南省分区综合时-面-α，β的关系表查得水库1~24h设计面暴雨量，由云南省1d暴雨分区综合雨型表得水库24h设计面暴雨时程分配过程，即表1中采用暴雨图集法推求计算的结果。

3.2暴雨移置法

3.2.1水库流域下游的河底岗雨量站有实测暴雨资料

2005～2017

年短历时暴雨资料和1979～2017年最大1d暴雨资料系列，邻近流域的大兴雨量站有1979～2017年长系列短历时暴雨资料。河底岗、大兴雨量站与水库流域中心直线距离分别为16km和16.5km，根据水库邻近雨量站点分布及其实测暴雨资料情况，选取

表1

水库24h设计面暴雨时程分配过程

暴雨图集法/mm

暴雨移置法/mm

T/h3.33%P=5%

P=3.610%

P=10.33%P=0.33%P=P=26..04.22..73.33%P=5%

6.010%

P=37..714.2

3.05.049.1

5.74.7

6.53.3511.87.4

5.13.310.06.6

4.05.3112.674.79.8

9.02.57.6617.611.262.453.4

8.111.467.97.256.3

6.17

96.187.867.65.34.916.611.47.7.14.94..112.38.510.5

7.88.9106...13.85.23.53.26.6116.14.23.83.55.83.93.62.7125.73.93.63.20.60.40.43.0135.33.73.43.00.60.40.40.3145.03.53.22.93.22.22.00.31.83.33.02.73.92.62.41.71.63.22.92.61.30.90.82.0174.43.02.82.52.61.71.60.7184.22.92.72.40.60.40.41.3194.12.82.62.33.22.22.00.3203.92.72.52.21.30.90.81.7213.82.62.42.19.76.56.00.75.0223.72.52.32.07.14.8233.62.52.32.08.45.74.43.724

3.52.42.21.95.23.55.23.42.32.1

1.91.93.9

1.33.24.31.8

2.61.22.72.4

1.02.0

这两站的实测暴雨资料对水库及邻近区域暴雨特性进行分析。

对大兴雨量站已有的实测最大1，6，24h暴雨

量与日暴雨建立相关关系，0.481相关系数达、K=0.7290.91、，K另外两个时段的相关系数也在其中得出综合折算系数24h暴雨与日暴雨K1624=1.090，=0.6

以上，说明该区域暴雨时程分配相对规律，用日暴雨资料推求短历时暴雨可行。

根据大兴雨量站1，6，24h暴雨量与日暴雨相关系数，将河底岗站实测日暴雨系列插补延长缺测年份的39a1，6，24h暴雨，得到该站1979～2017年（n=主，）人工校测，1，6，24h暴雨资料系列长度已超过短历时暴雨系列。大兴站以自记为30a，资料可靠且系列具有一定代表性。虽然河底岗站只有10a的实测1，6，24h短历时暴雨系列，但最大日暴雨系列已达32a，且短历时暴雨系列最大值和日暴雨系列最大值均发生在2014年，故通过展延后的短历时暴雨系列资料基本可靠。

依据河底岗、大兴雨量站1979～2017·年25（·n=

2020年5月水利水电快报EWRHI第41卷第5期39a频率，）1，矩法初估统计参数均值和6，24h暴雨系列，用数学期望公式计算经验Cv值，取Cs=3.5Cv，以P-Ⅲ型曲线为线型，用目估适线法确定统计参数，并与河底岗、大兴雨量站1，6，24h查图值进行

比较，见表2。

表2

水库暴雨统计参数

站名方法P1h/mmP6h/mmP24h/mm

Cv1Cv6Cv24

Cs

河底岗查图

实测

36.048.0大兴

查图35.5.677.00.430.520.490.503.5实测

35.036.7

0.430.500.483.5Cv48.078.155.6

78.083.2

0.43

0.480.40

3.5Cv0.45Cv3.2.20.423.5Cv

由表暴雨统计参数移置

2可知，河底岗雨量站1，6，24h实测暴雨

均值与Cv值均比查图值大；大兴雨量站1，6，24h实测暴雨均值比查图值大、Cv值比查图值小。水库流域与大兴和河底岗雨量站直线距离均较近，但河底岗雨量站与水库为同一河谷区，属上下游关系，故河底岗站暴雨特性更能代表水库流域，2014遇，为有实测记录以来的最大值，年最大24h暴雨量达178.5mm且河底岗站实测暴雨资料系，约为66a一列样本较大兴雨量站更具代表性。故移置河底岗站实测暴雨资料推求水库设计暴雨更合理。

分别用暴雨图集查出河底岗雨量站和水库流域中心的1，6，24h暴雨及其参数，根据两个位置查图值的差异计算修正系数，对河底岗站暴雨系列资料进行修正后移置到水库流域，暴雨参数移置成果见表3。

表3

临沧市水库流域1，6，24h暴雨参数移置成果

站名方法P1h/mmPCv1

Cv6

Cv24

Cs

河底岗查图

水库

查图

36.06h48.0/mmP24h/mm

35.00.430.490.48河底岗实测35.8.077.00.483.53.5Cv0.430.49Cv水库移置成果

34.5

.677.0其余修正系数为.6

78.178.1

0.520.501。

0.500.500.503.53.5Cv0.52Cv

注：P1h的修正系数为0.972，在河底岗站实测暴雨资料系列中挑选降水量较大、雨型对工程最不利的暴雨过程为典型暴雨，进行设计暴雨推算，“20140721”实测暴雨过程见图2。

采用1，6，24h移置暴雨成果控制作同频率放大典型暴雨过程，其同频率放大系数按式（1）～（3）计算。

αα1h=H1t/H1d

（α6h=（H6t-H1t）（/H6d-（1）·26·

24h=（H24t-H6t）（/H24d-HH1d6d）

）

（23））

12345671011时间1213/h

01415161718192021222324

105

mm/量15水20降253030图2河底岗雨量站“20140721”实测暴雨过程

式中，αH1h，α6h，α24h分别为1，6，24h放大系数；H1t，H为1，6，24h，典型暴雨量。

6，24h设计点暴雨量；H6t24t分别为11d、H6d、H24d分别，《云南省暴雨径流查算实用手册》［5］

查得水库流

域的点、面换算系数，将水库设计点暴雨成果乘点面系数得设计面暴雨过程，即表1中采用暴雨移置法推求计算的结果。

3.3设计暴雨成果选用

从表3可知，移置邻近实测暴雨资料得到的暴

雨统计参数比采用查图法得到的暴雨统计参数更大，主要是由于云南省暴雨统计参数等值线图集资料截止2000年，而实测资料暴雨资料系列现已增加了17a，特别2014年水库流域发生了较大暴雨，水库下游的河底岗雨量站日暴雨达168.5mm，为39a实测最大值，邻近的大兴雨量站日雨量达107mm，为39a实测第二大值。根据SL44-2006《水利水电

工程设计洪水计算规范》［1］

的规定，从通过审批暴雨

统计参数等值线图上查算工程所需历时的设计点暴雨量。当本地区或邻近地区发生大暴雨，或依据不同年代图集查算成果差别较大时，应对查算图表进行合理性检查，必要时作适当调整。因此，水库设计暴雨成果采用暴雨移置法推求的成果。

4

设计洪水推求

4.1

由设计暴雨推求设计洪水

通过查《云南省暴雨径流查算实用手册》相关

图集得水库产汇流分区，然后根据水库产流分区和设计暴雨过程进行产流计算，采用初损后损法推算得到水库设计净雨过程。根据水库设计净雨过程、汇流分区及流域特征参数，用纳西瞬时单位线法推求出水库断面的设计洪水过程，水库设计洪峰和最大24h洪量见表4。

李增国等

表4

标准P/%0.330.5125

无资料小流域水库设计洪水计算方法分析

水库设计洪水成果

3

Qm.p/(m·s-1)

图推求。本文分析结果表明，移置河底岗雨量站暴

W24.p（/万m3）209.1194.0165.8139.5121.4107.1

雨资料的方法推求设计暴雨成果更符合流域实际，推求的水库设计洪水也更合理。

对于西南湿润地区的无资料区域，中小型水库设计洪水以暴雨推求洪水的方法为主，其设计暴雨可通过查暴雨统计参数等值线图、移置邻近实测暴雨资料、建立区域暴雨高程关系等方法推求，但应在充分分析流域暴雨洪水特性的基础上确定合理的设计暴雨推求方法。

参考文献：

［1］SL44-2006水利水电工程设计洪水计算规范［S］.［2］梁忠民，钟平安，华家鹏.水文水利计算［M］.北京：中

国水利水电出版社，2008.

［3］金光炎.水文统计理论与实践［M］.南京：东南大学出

版社，2012.

［4］云南省水利厅.云南省暴雨统计参数图集［R］.昆明：

云南省水利厅，2007.

［5］云南省水利厅.云南省暴雨洪水查算实用手册［R］.昆

明：云南省水利厅，1992.

（编辑：李慧）

131.0122.0106.092.281.7

3.33

71.9

4.2设计洪水成果合理性分析

2014年该水库流域发生了大洪水，根据当地水

文部门2014年7月25日对水库流域的洪水调查成

果，其洪峰流量为88m3/s，重现期为66a一遇。而用洪调成果相对偏差为8.8%，说明暴雨移置法推算的设计洪水成果合理。

暴雨移置法推求得水库66a一遇洪水为96.5m3/s，与

5结语

水库流域属于无资料地区，仅水库流域下游的河底岗有暴雨资料。根据资料条件，采用暴雨图集法计算设计洪水，而设计暴雨分别采用移置河底岗雨量站暴雨资料和查云南省暴雨统计参数等值线

Discussionondesignfloodcalculatingmethodforreservoirsinsmallungauged

watershed

(1.LincangHydrologyandWaterResourcesBureauofYunnanProvince，Lincang677000,China；

2.YuxiHydrologyandWaterResourcesBureauofYunnanProvince，Yuxi653100,China)

LIZengguo1,ZHANGNan2,ZHANGZhenjie1,GUOLingmei1

Abstract:Reservoirdesignfloodcandeterminethereservoirscaleandthelayoutoftheflooddischargefacilities.

Forsmallwatershedareaswithouthydrologicalrecord,theresultsobtainedbytraditionalfloodcalculatingmethodstransplantingmethodbasedontheanalysisofthecharacteristicsoftheregionalrainstormandfloodisproposed.consistentwiththefloodinvestigation.method;smallwatershed;humidarea

Theproposedmethodusesthemeasuredrainstormdataofthenearbyareatocalculatethereservoirdesignflood.By

maynotbeconsistentwiththeactualcharacteristicsofrainstormandflood.Therefore,ameasuredrainstormusingthismethod,the66-yearfrequencyfloodofareservoirinLinyiCityiscalculated(96.5m3/s)andtheresultisKeywords:floodcalculation;measuredrainstormdataofthenearbyarea;ameasuredrainstormtransplanting

·27·