1. 引言
大多数经济发达国家在确定河流的防洪标准时,考虑了其防洪保护对象的类型、级别和重要性[1]。城市的防洪标准一般为50~200年一遇,一些国家规定特别重要城市可达到500~1000年一遇防洪标准。肖义[2]等建议我国把基于风险分析的设计洪水标准作为防洪标准规范的使用参考。洪水灾害风险可定义为不同强度洪水发生的概率及其可能造成的洪水灾害损失。洪水灾害风险分析的主要工作是风险识别和风险估计[1]。目前,研究较多、理论体系相对健全的防洪标准确定方法主要有经济分析方法、风险分析方法和综合评价模型法[3];由于我国山洪灾害防治任务的复杂性和艰巨性,目前山洪灾害防治工作整体上还处于初级阶段,防汛方面存在的风险主要包括五方面:遇超标洪水或标准内洪水,堤防仍有发生漫顶和溃决的可能;中小水库尤其是病险水库失事风险较大[4],本研究通过总结反思,按照“尊重自然,顺应自然规律,与自然和谐相处”的要求[5],提高应对镇朔湖拦洪库超标准洪水致灾风险分析的能力,实现最终防御。
贺兰山是宁夏青铜峡灌区的天然屏障,既可以阻止西部寒流及风沙对灌区的侵袭,又能拦截东麓暖湿空气,并使其抬升形成局地暴雨。降雨(暴雨)形成的地面径流多以洪水形式出现,具有干旱、半干旱区暴雨洪水的一般特性;季节性特点较强;暴雨历时不长,大强度暴雨持续时间短;笼罩面积不大,点面折减快。镇朔湖拦洪库位于宁夏回族自治区石嘴山市平罗县崇岗镇境内,是宁夏贺兰山东麓防洪工程体系的重要组成部分,镇朔湖拦洪库建成于1974年,成功抵御了“2016.8.21”“2018.7.22”等区域暴雨洪水过程,有效发挥滞洪调峰作用,确保了区域防洪安全。2019年实施了镇朔湖拦洪库改造修复工程,工程建成后,有效提升了水库拦蓄洪水、错峰调节的能力。镇朔湖拦洪库防洪标准为20年一遇,总库容1857万立方米,防洪库容1152万立方米。
镇朔湖拦洪库超标准洪水致灾风险分析与防御措施研究分析了镇朔湖拦洪库上游大西伏沟、小西伏沟、桃箕沟、大水沟、和小水沟等山洪沟发生洪水超标准洪水的致灾风险,提出了超标准洪水防御措施,强化了超标准洪水情况下防汛调度决策理论技术支撑,对于保护区域第二农场渠、包兰铁路、银石线等交通、灌溉设施及平罗县崇岗镇、国营前进农场约2.3万人生命财产安全和8万亩农田的防洪安全有积极作用。同时,镇朔湖拦洪库超标准洪水致灾风险分析与防御措施研究,对于镇朔湖拦洪库区域防洪工程规划建设,完善贺兰山东麓防洪工程体系也具有积极作用。
2. 资料和研究方法
2.1. 研究区工程体系布置情况
镇朔湖库防洪体系如图1所示,由导洪工程、行洪沟道工程(调洪库)、镇朔湖拦洪库、泄洪沟(闸)工程组成。贺兰山自大西伏沟至大水沟、小水沟间沟道洪水,出山后经贺兰山坡地部分沟道导洪堤、自然沟道、调洪库下泄至110国道,拦洪再经南大闸、一号、二号、北大闸沟四条行洪沟道归入镇朔湖拦洪库。洪水经拦蓄后,由退水沟、泄水闸进入三二支沟、金色河入瀚泉海,最后经典农河(第三排水沟)排入黄河。
图1. 研究区工程体系布置图
2.2. 研究方法
根据地区水文特性,利用实测资料和《宁夏暴雨洪水图集》对镇朔湖拦洪库坝址以上流域的设计暴雨进行分析,推求镇朔湖拦洪库坝址以上集水区域的设计面暴雨过程,进而利用产流方案推求设计净雨过程,再由设计净雨过程求得设计洪水过程。
汇流计算将采用纳希瞬时单位线法,一般来讲,单位线在物理上具有5个特征,超渗降雨在有效时间内强度恒定;超渗降雨均匀分布在整个汇水区;给定持续时间的超渗降雨造成的直接径流持续时间是恒定的;所有直接径流持续时间随时间增加,每个单位线的直接径流总量成正比增加;在一个汇水区内,给定一个超渗降雨产生的汇水区径流特征不变,上述特征能够反映所关注的汇水区降水径流特征[6]。镇朔湖拦洪库坝址以上流域一般为超渗产流,洪水过程涨落迅速,该地区均采用纳希瞬时单位线法进行汇流推算[7]。
2.2.1. 产汇流计算
拦洪库所在流域洪水多为超渗产流,历时短,洪峰高。根据该区域暴雨特性及下垫面条件采用扣损法计算产流过程,其产流计算公式如式(1)所示。
(1)
式中:
为产流期某时段净雨,mm;
为产流期某时段面雨量,mm;
为产流期某时段损失雨量,mm。
汇流计算是根据推求的设计净雨过程
推求流域出口断面的洪水流量过程。按瞬时单位线原理,就是净雨历时趋于零的极小时段内产生的一个单位净雨(1 mm净雨深)在出口处所形成的流量过程线,其公式为:
(2)
式中:
为时段变量;
为
时瞬时单位线纵坐标;
、
为瞬时单位线形状参数;
为伽马函数。
根据上式可导出时段单位线方程式为:
(3)
式中:
表示面积为F的流域上t时段内净雨量为1 mm的单位线在t时段的纵坐标,m3/s;
为计算时段长,h;
为换算系数;
为纳希瞬时单位的累积曲线,又称
曲线;
为
曲线在时间坐标轴上错后Δt时间,用
表示。
2.2.2. 调洪计算
1) 水库调洪计算原理
水库调洪计算原理为水量平衡[8],计算公式如下:
(4)
式中:Q1、Q2表示计算时段初、末的入库流量,m3/s;
表示计算时段中的平均入库流量,m3/s;q1、q2为计算时段初、末的下泄流量,m3/s;
为计算时段中的平均下泄流量,m3/s;V1、V2为计算时段初、末的库存水量,m3;Δt为计算时段,取Δt = 3600 s。
2) 泄洪闸泄流计算公式
t泄洪闸泄流计算参考《灌溉与排水工程设计规范》(GB50288-99)中的涵洞(隧洞)水力计算方法,计算公式如下:
(5)
式中:Q为涵洞设计流量,m3/s;m为无压流时的流量系数,m = 0.385;m1为半压力流时的流量系数,八字墙
进口,m1 = 0.670;m2为压力流时的流量系数,
;
为局部水头损失系数,进口
为渐变段,取
;C为谢才系数,C = R1/6/n;B为涵洞底宽,m;A为涵洞断面面积,m2;L为涵洞长度,m;R为涵洞水力半径;n为涵洞糙率;i为涵洞底部坡降;
、
位修正系数,根据相关规范,取
,
;D为涵洞高度,m;h1为洞前水深,m;H0为行进流速水头的洞前进口水头,m。
3) 溢洪道泄流计算
水库溢洪道的泄流量按宽顶堰流按下式进行计算:
(6)
式中:ε为侧收缩系数,对溢流堰取ε = 0.90;σ为淹没系数;m为流量系数,取0.385;b为堰顶宽度,m;H0包括行近流速的水头,近似取为上游水头,m。
2.2.3. 导洪堤过流能力和冲深计算
导洪堤过流能力主要根据水力学原理,采用比降面积法,按照洪水调查规范要求,调查计算出所需参数,计算不同水深对应流量的一种方法[9]。具体采用如下公式计算导洪堤不同部位的过流能力(本次仅计算最大过流能力),即:
(7)
(8)
式中:i为导洪堤不同断面位置;V为各断面平均流速,m/s;Q为流量,m3/s;R为水力半径,m;J为水面比降;B为断面平均宽度面积,m2。
导洪堤的冲刷按照《水力计算手册》中岸坡冲深公式计算[10],即:
(9)
式中:
为从沟底算起的局部冲深,m;α为水流流向与岸坡交角,度;m为导洪堤迎水面边坡系数;d为坡脚处泥沙计算粒径,m;
为局部冲沙流速,本次计算采用设计流量值除以过水断面面积,m/s。
3. 现状防洪工程体系行洪能力复核
3.1. 产汇流及调洪计算结果
在计算瞬时单位线,进行洪水演进时,须求得参数n、K和滞时M,将依据宁夏回族自治区水旱灾害防御中心提供的镇朔湖拦洪库相关基础数据进行推求,其计算公式如下:
(10)
式中:i为产流区平均净雨强度,mm/h;L为概化长度,km。据此求得镇朔湖拦洪库设计洪水计算成果见表1和图2。
表1. 镇朔湖拦洪库设计洪水成果表
拦洪库名称 |
类别 |
设计频率P (%) |
1 |
2 |
3.33 |
5 |
10 |
镇朔湖拦洪库 |
洪峰流量(m3/s) |
1350 |
1120 |
920 |
773 |
453 |
洪水总量(m3) |
1470 |
1200 |
974 |
805 |
471 |
图2. 镇朔湖拦洪库设计洪水过程线
根据《平罗县镇朔湖拦洪库调度运行计划》,水库汛期最低水位主要受三二支沟控制,水库段三二支沟现状沟底高程1096.1~1096.00 m,汛期正常水位1098.21~1097.95 m,水库调洪计算时釆用三二支沟汛期平均正常水位1098.2 m作为水库最低水位。
根据水文计算成果,本次镇朔湖拦洪库设计洪水标准20年一遇,校核洪水标准100年一遇。设计20年一遇洪水总量805万m3,100年一遇洪水总量1470万m3。计算镇朔湖拦洪库泥沙入库年淤积量约23.5万m3,本次淤积以30年计,总淤积量约705万m3。受三二支沟水位顶托作用,拦洪库死水位高程为1098.10 m,淤积库容以上形成325万m3,死库容。镇朔湖拦洪库设计洪水标准下库容为1835万m3,与此对应的设计洪水位为1099.15 m;校核洪水标准下库容为2500万m3,与此对应的校核洪水位为1100.05 m。有效防洪库容1470万m3,加上淤积库容与死库容后,设计总库容为2500万m3。
3.2. 行洪能力复核
3.2.1. 上游工程体系
在上游防洪体系中,跨山洪沟道以桥梁为控制性断面计算,涉及到的控制性桥梁断面为五分沟桥、三号沟桥、北大闸沟桥、大水沟一号沟南桥和北桥、大水沟二号沟南桥和北桥、箱涵沟桥和平贺桥,五分沟桥不受沟道洪水直接影响,因此不进行桥位设计洪水计算。汇流计算面积采用Google卫星图量算和航拍勘查综合确定,坡地面积按50%折减,并结合各控制性断面的测绘结果、集水区域面积(分为山地面积和坡地面积)、概化长度以及沟道上游来流情况,以北大闸沟桥为例,进行水文计算和行洪复核。
1) 水文计算
北大闸沟桥选取各沟道沟口位置和桥位所在位置作为设计洪水计算断面,上游的集水区域和沟道相关测量值如表2所示。其断面的具体位置为小水沟导洪工程出山口和原110国道过水路面上游200 m处(见图3)。
表2. 小水沟及其南侧小沟计算面积表
沟道名称 |
山地面积(km2) |
坡地面积(km2) |
计算面积(km2) |
概化长度(km) |
小水沟 |
69.85 |
6.31 |
73 |
15.8 |
小水沟附近小沟 |
8 |
9.51 |
3.7 |
图3. 北大闸沟桥来水流路及控制断面位置选取
小水沟沟口及其南侧小沟的设计洪水计算成果见表3。
表3. 小水沟沟口及其南侧小沟设计洪水计算成果表
沟道 |
计算断面 |
设计洪水 |
洪水频率P (%) |
5 |
2 |
1 |
小水沟 |
沟口1# |
洪峰流量(m3/s) |
257 |
377 |
465 |
洪水总量(万m3) |
180 |
263 |
325 |
小水沟南侧小沟 |
沟口2# |
洪峰流量(m3/s) |
40.7 |
60.4 |
76.8 |
洪水总量(万m3) |
17.1 |
25.4 |
32.2 |
小水沟及南侧小沟坡面 |
原G110国道 |
洪峰流量(m3/s) |
13.5 |
20.0 |
25.4 |
洪水总量(万m3) |
5.4 |
8.0 |
10.2 |
根据宁夏水文中心在小水沟沟口的洪水调查断面,连续洪水调查资料系列为1979~2018年,并有1955年、1975年历史调查资料,据此,采用P-III型曲线进行适线计算,将历年调查洪水进行频率计算,其相应的计算成果如表4所示。
表4. 小水沟沟口调查洪水频率计算成果表
沟道 |
计算断面 |
设计洪水 |
洪水频率P (%) |
5 |
2 |
1 |
小水沟 |
沟口 |
洪峰流量(m3/s) |
310 |
450 |
560 |
两种算法成果进行对比分析,考虑工程的安全,两者釆用较大值,即小水沟沟口设计洪水釆用洪水调查资料频率计算的成果。根据小水沟导洪堤选取的分析断面位置,原110国道过水路面上游200 m洪水按照沟口和小水沟南侧小沟的设计洪水叠加计算,其百年一遇设计频率洪峰流量为637 m3/s。
2) 行洪复核
根据小水沟区域的两个控制断面和水位–流量关系曲线(见图4、图5),以100年重现期的频率洪水开展行洪复核计算,由于原110国道过水路面上游200 m断面处沟道较深,本次不计算其冲深。小水沟导洪堤各不同断面位置的过流能力及冲深计算成果见表5。
图4. 小水沟导洪堤沟口实测大断面图及水位–流量关系
图5. 小水沟导洪堤原110国道上游200 m实测大断面图及水位–流量关系
表5. 小水沟沟口调查洪水频率计算成果表
序号 |
沟道 |
断面名称 |
设计洪水(m3/s) |
漫堤流量(m3/s) |
冲深(m) |
1 |
小水沟 |
沟口断面 |
560 |
671 |
0.92 |
2 |
原110国道过水路面上游200 m |
637 |
1295 |
- |
由表5可得,小水沟发生100年一遇洪水时,沟口断面和原110国道过水路面上游200 m断面的设计流量分别为560 m3/s、637 m3/s,均小于这两个断面的漫堤流量(671 m3/s和1295 m3/s),可通过标准内洪水;小水沟沟口断面处在100年一遇洪水时,冲深为0.92 m,小于导洪堤基础埋深2 m,则小水沟导洪堤现状基本满足100年一遇洪水过流。
经上述分析,小水沟现状基本可满足100年一遇洪水过流,则小水沟洪水出沟口后归拢下泄,不会串流至附近其他沟道。小水沟及其南侧小沟的洪水出沟口沿着导洪堤下泄至北大闸沟桥,在此过程中洪峰流量会消减坦化,本次考虑小水沟沟口与桥梁区间增加的坡积洪水与消减坦化的洪水相抵消,故北大闸沟桥桥位断面设计洪水直接釆用小水沟及其南侧小沟沟口汇合后设计洪水。北大闸沟桥100年一遇设计洪峰流量为637 m3/s,具体见表6。
表6. 北大闸沟桥桥位断面设计洪水计算成果
沟道 |
计算断面 |
洪峰流量(m3/s) |
洪水频率P (%) |
5 |
2 |
1 |
小水沟 |
沟口 |
310 |
450 |
560 |
小水沟南侧小沟 |
40.7 |
60.4 |
76.8 |
- |
北大闸沟桥 |
351 |
510 |
637 |
3.2.2. 下游工程体系
研究区内的下游退水工事主要由三二支沟承担,该支沟为人工修建的泄洪沟道,假设该过程为均匀流,根据现状沟道的流量、比降及断面尺寸,采用谢才公式和曼宁公式可计算出桩号0+000~7+357段不同分段水深,计算结果如表7所示。由表7可知,根据现状沟道沟深在3.0~4.0 m左右,则断面满足镇朔湖泄洪(8 m3/s)和现状农田排水(6.67 m3/s)共14.67 m3/s的过流能力需求。
表7. 三二支沟各沟道段的过流能力计算成果
桩号
(km + m) |
流量(m3/s) |
底宽(m) |
水深(m) |
比降 |
边坡 |
糙率 |
流速(m/s) |
过水面积(km2) |
0+000~2+150 |
6.67 |
12 |
1.16 |
6500 |
2 |
0.03 |
0.403 |
16.539 |
2+150~3+000 |
14.67 |
30 |
1.09 |
6500 |
2 |
0.03 |
0.416 |
35.246 |
3+000~7+150 |
14.67 |
30 |
1.06 |
6000 |
2 |
0.03 |
0.426 |
34.365 |
7+150~7+357 |
14.67 |
30 |
1.30 |
6000 |
2 |
0.03 |
0.479 |
30.60 |
4. 超标准洪水应对措施
通过本文研究分析,镇朔湖拦洪库自身具有较大调蓄空间,标准内洪水通过水库泄洪设施退泄洪水。水库上游沟道普遍发生超标准洪水情况,仅通过水库泄洪设施难以有序退泄洪水,釆取破堤分洪的方式应对超标准洪水。
4.1. 危险区基本情况
根据2016年《平罗县山洪灾害普查危险区划定报告》和2019年《平罗县山洪灾害调查及预警指标复核》,危险区内主要有学校、企业、养殖场等,按行政村及重要场所将平罗县危险区划定为8处,危险区分区见图6,其山洪灾害危险区基本情况见表8。根据各危险区的房屋数和危险区人口,从水库下游进行应对措施分析。
表8. 平罗县山洪灾害危险区基本情况表
序号 |
乡镇 |
危险区地点 |
家庭户数(户) |
危险区人口(人) |
房屋数(座) |
1 |
崇岗镇 |
崇岗村银汝路 |
20 |
70 |
20 |
2 |
崇岗村4队 |
50 |
211 |
50 |
3 |
崇岗村姚汝路 |
1 |
3 |
1 |
4 |
崇岗村2队 |
6 |
21 |
6 |
5 |
崇岗村汝箕沟 |
100 |
237 |
100 |
6 |
崇岗村中耀路 |
28 |
116 |
28 |
7 |
崇富村1队(小水沟) |
13 |
19 |
13 |
8 |
下庙10队 |
46 |
194 |
46 |
9 |
下庙1队 |
20 |
72 |
20 |
10 |
下庙村2队 |
52 |
215 |
52 |
11 |
常青村1队 |
42 |
170 |
42 |
12 |
下庙6、9队(二) |
24 |
106 |
24 |
13 |
下庙5队 |
51 |
220 |
51 |
14 |
下庙3队(一) |
33 |
134 |
52 |
15 |
下庙6、9队(一) |
22 |
81 |
22 |
16 |
下庙4队 |
33 |
136 |
33 |
17 |
下庙3队(二) |
11 |
23 |
11 |
18 |
常青村4队 |
56 |
224 |
56 |
19 |
常青5队(一) |
41 |
165 |
41 |
20 |
常青村6队 |
24 |
55 |
24 |
21 |
常青2、7队(二) |
28 |
112 |
28 |
22 |
常青2、7队(一) |
59 |
238 |
59 |
23 |
常青村8队(一) |
19 |
50 |
19 |
24 |
常青村8队(二) |
33 |
135 |
33 |
25 |
常青3队 |
93 |
356 |
93 |
26 |
常青5队(二) |
23 |
79 |
23 |
27 |
常青2队 |
4 |
6 |
4 |
28 |
红崖子乡 |
王家沟村4队 |
49 |
189 |
49 |
29 |
王家沟村3队 |
55 |
244 |
55 |
30 |
王家沟村2队 |
58 |
245 |
58 |
4.2. 下游应对方案
利用镇朔湖拦洪库泄洪闸分洪,在退水闸和溢洪道处设置紧急开口,超标准洪水时打开开口进行泄洪。在三二支沟和沙湖公路及包银铁路之间的区域形成非常洪水淹没区,面积约为19 km2,分洪区地面高程1098.10~1101.8 m之间,见图6,该区域地形具备分洪条件,可采取相应工程措施,预估库容1000万m3。根据超标准洪水设定分洪和不分洪方案,得出各方案灾损情况如表9所示。
表9. 不同方案灾损情况
是否采取非常措施 |
淹没面积(km2) |
淹没居民地面积
(万m2) |
淹没公路(km) |
受影响人口
(万人) |
淹没农业面积
(公顷) |
受影响GDP
(万元) |
直接经济损失
(万元) |
否 |
3.586 |
14.658 |
6.041 |
0.3231 |
5.32 |
1561.758 |
161.491 |
是 |
1.738 |
5.231 |
4.251 |
0.0148 |
1.27 |
1007.227 |
56.589 |
图6. 分洪区示意图
5. 结论
镇朔湖拦洪库防洪体系通过多年的建设,现已初步形成了“导引、滞蓄、控制排泄”的防洪体系,在防洪调度、保障生命财产安全方面发挥了重大作用。但由于工程体系不完善,防洪标准低,应对超标准洪水和突发事件的能力不足,通过分析研究主要结论如下:
1) 以小水沟为例,开展水文计算和行洪能力复核得出:小水沟现状基本可满足100年一遇洪水过流,在出沟口后可归拢下泄,在经过导洪堤的消减、坦化之后,与小水沟沟口到桥梁之间的区间汇流相互抵消,故北大闸沟桥桥位断面设计洪水直接釆用小水沟及其南侧小沟沟口汇合后设计洪水,则其100年一遇的设计洪峰流量为637 m3/s。
2) 结合现场调研和资料收集工作,对平罗县行政范围内的山洪危险区进行分析,并据此对超标准洪水条件下,为镇朔湖拦洪库上、下游提出“上拦蓄,下分洪”的山洪防御措施,以此来缓解超标准洪水对拦洪库和危险区人民的危害。
基金项目
本研究由国家重点研发计划项目(2023YFC3209305、2023YFC3209303)、黄河水利科学研究院中央级公益性科研院所基本科研业务费专项资金项目(HKY-JBYW-2024-12、HKY-YF-2024-05、HKY-JBYW-2025-08)等资助。
NOTES
*通讯作者简介:杨东旭(1987-),男,宁夏银川人,硕士研究生,工程师,从事水利工程规划设计、建设、运行管理及防汛抗旱工作方面的研究,Email: 31970086@qq.com