进一步对寿光市上口镇三日和24小时降雨过程数据分析。自18日6:00 ~20日2:00,持续降雨45小时,累计降雨量401.5mm,雨峰在三日最后。其中最大24小时降雨发生在19日0时-19日24时,累计降雨量346.1mm;最大2小时降雨发生在8月19日15-16时,累计降雨量100.0mm。峰值雨量发生在降雨过程后期,前期土壤在小降雨下达到饱和状态,后期大暴雨来临产流系数高,属于极不利情形。
弥河中下游三日降雨时程分布图
《寿光市防汛预警预案》定义红色I级预警为:“雨量已经达到100mm,降雨可能持续,且预报降雨区域在过去15天内累积降雨达到200mm以上…”。从三日、最大24小时雨量可以判断,应该是I级红色预警。形势异常紧急。
3.4 819洪水水库调度方式评估
3.4.1 冶源水库调蓄分析
水库库区自2018年8月18日06:00~2018年8月19日 21:00,持续降雨40小时,累计降雨量241mm,最大2小时降雨发生在降雨末期,累计降雨量82mm,雨峰在后,加剧了洪水风险。整场降雨过程中,冶源水库累计入库水量1.11亿m3,累计出库水量0.99亿m3,水库调蓄水量1240万m3。
冶源819暴雨实际水库调度流量过程
据新闻资料,在暴雨到来前,水库水位实现预降至136.2m,相应库容8200万m3。宜水环境以此为初始水位,经调洪演算得出,汛中最高水位138.5m;,汛中最大蓄水量10800万m3,本次峰段洪水利用了调洪库容2600万m3。
冶源819暴雨调洪演算水库水位、库容变化过程
3.4.2 冶源水库防洪库容潜力估计
为了利用更多调洪库容,假设当时调度时冶源水库下泄控制至300m3/s(实际700m3/s)。调洪演算表明冶源水库累计入库水量1.11亿m3,出库水量从现状0.99亿m3,降低到0.71亿m3;整场洪水水库调蓄4000万m3,比实际增加2600万m3。
汛中最高水位从现状138.5m升高到140.23m;汛中最大库容从现状10800万m3升高到13100万m3,实际峰段调洪容积从现状2600万m3升高到4900万m3。
冶源水库校核洪水位为141.95 m,上游移民高程线在139.0 m(摘自2016.4 冶源水库提高汛限水位的可行性分析)。这个也是约束条件,若减少下泄量,水库水位将上升至140.23m,则库区将产生淹没损失。从中可见,冶源水库理论上可以继续利用调洪库容2000万m3左右,但会因此而加重上游的淹没损失。
819入库流量峰值仅发生于雨峰之后2小时左右,峰值逼近2000m3/s,每小时入库700万m3水量。如果降雨继续,几小时后水库将超越校核洪水位,保坝泄流或溃坝风险大大增加。
上述分析表明,819暴雨水库调度,冶源水库具备额外2600万m3相对安全可用的调洪容积。但前提是必须有准确的至少未来6-12小时气象预报和先进的基于实时模型的水库流域洪水预警预报系统。否则决策的风险极大。
冶源819暴雨水库下泄预案对比示意图(下泄700m3/s与300 m3/s对比)
冶源819暴雨调洪演算水库水位变化(下泄700m3/s与300 m3/s对比)
3.4.3 冶源水库防洪库容潜力利用效果分析
根据上述假设,即冶源水库最大下泄流量控制为300m3/s,利用模型模拟,得出淹没范围及风险等级。
根据实际下泄和本方案模拟成果,对洪水风险进行对比分析。结果表明,本方案控制冶源水库下泄流量不超过300m3/s,高风险区域面积减少2.6km2,削减率3.5%;中风险区域面积减少0.21km2,削减率0.5%,低风险区域面积降低0.77km2,削减率0.8%;整体风险区域减少3.6km2,削减率1.7%。整体削减效果并不明显,说明对减轻下游洪水灾情作用甚微。
具体结果如下:
上表对比分析说明,819暴雨条件下,利用冶源水库防洪库容潜力,即下泄流量控制在300m3/s,下游洪水风险并未有显著减小,对减轻下游洪水灾情作用甚微。说明,819暴雨灾区本地洪水影响占相当大作用。
3.4.4 库群调度分析
19日暴雨前,水库水位估计在汛限水位附近(缺乏实测资料)。根据实际发布的819暴雨水库下泄流量(下图)冶源水库自19日中午开始从200-400m3/s 下泄,黑虎山水库19日中午开始从80-160m3/s 下泄,说明根据气象预报,调度管理部门采取了预泄措施。
降雨和水库入流峰段黑虎山水库大约下泄900m3/s, 源水库控制了下泄,并实现了错峰。
没有嵩山水库的水位和出库流量资料。
暴雨水库实际下泄流量过程
新闻资料显示,黑虎山水库流域面积190km2,总库容5600万m3,历史最高水位165.38米,警戒水位168.19米。819暴雨黑虎山水库最高水位达到167.84米。8月24日上午8时黑虎山水库水位为164.22米,仍超过汛末水位1.22米。估计前期水位略有偏高,导致19日晚雨峰到来时迅速优先下泄,缓解大坝险情。
3.4.5 冶源水库汛限水位抬高影响评估
根据2016年4月《山东水利》发表的林立廷等关于《冶源水库提供汛限水位可行性分析》,冶源水库原来的汛限水位是136.5m (相应库容7935 万m3),拟提高到137.72 m (相 应库容9678 万m3),增加约1700万m3兴利库容,或减少1700万m3防洪库容。
以上述对“冶源水库防洪库容潜力利用效果分析”类推,假如调度采用原定的136.5 m汛限水位,增加1700万m3防洪库容,下游洪水淹没影响面积可能可以减少1.7%左右。
在缺乏对抬高汛限水位可研成果认真学习的前提下,简单的认可或否定此可研的结论是对复杂技术工作的不尊重或轻视。我们认为,水库特征水位的确定,需要长期水文气象资料,利益影响范围内社会经济和城市发展规划资料,并利用水资源、防洪、经济等技术分析,综合确定。由于社会发展水资源需求、用水结构特征变化、水文气象资料的不断丰富和计算技术、管理水平的提高,这些特征水位的修正并不少见。
但是特征水位的修正需要考虑配套的防汛预案、保护范围防洪工程建设、预警预报和风险评估调度技术、洪水风险图的发布与应用等相应措施。
3.4.6 水库调度评价
综合上述分析,在19日的水库水情和紧迫雨情条件下,水库调度的技术决策基本合理。
其改进的空间在于气象预报的准确性、实时流域洪水模型和降雨不确定性预报决策技术的应用、流域内水库群优化调度以及应急联合调度规则的研究与制定等。2018年8月19日8:00国家气象台发布了未来24小时降水量预报(见下图),其描述的总量、覆盖范围,与本研究团队获取的GPM雨量特征很相似。819暴雨防汛管理部门依据的气象预报和水情形势评估在技术上值得研究与回顾。
国家气象局2018.8.19 未来24小时暴雨预报
3.5 寿光下游河道断面平面走势特征
浏览Google 2018年影像资料,总体发现弥河寿光市下游开始河势逐渐呈流荡型,弯度增大。这是多年以来河流与自然动态平衡的结果。图中可见,弥河分段筑有橡胶坝,下游有分叉,河滩空间较大。这种河势,一旦发生大水,由于水流的动能和取直惯性,行洪滩地是否保护好非常重要。
2018年《寿光市弥河河道防汛应急预案》根据省防指《关于编制修订全省骨干河道防洪预案的通知》(鲁汛旱字[2004]第37号)要求,寿光市分流口控制断面现状水利工程情况下和弥河行洪能力,确定弥河干流洪水,从而划分出洪水等级,结合实测纵横断面资料成果,进行水面线推算。
洪水等级划分如下:
- “一般洪水”分流口断面流量≤2722m3/s。
- “现状标准洪水”分流口断面流量2722 m3/s~3943 m3/s。
- “超标准洪水”分流口断面流量≥3943 m3/s。
本研究估算得出弥河寿光以下流量不到上述“一般洪水”等级。
弥河寿光市区下游断面形式复杂,且受分叉控制工程以及分叉河道的影响,弥河上下游防洪能力、相应流量估算的技术方法很重要。这在很大程度上将影响工程的合理布局、建设投资和保护区洪水安全。
4、结论
根据降雨量空间分布、时程特征和模型成果分析,819寿光洪水主要成因为:
- 本地特大暴雨是导致寿光洪灾的主要原因。本地24小时大于300mm的极为不利的特大暴雨降雨条件;由于前二日近100mm左右降雨,导致819特大暴雨几乎100%产流;根据寿光市防汛预案,本地的暴雨量已超过I级红色预警,本身险情异常紧急;
- 弥河下游河道排洪能力衰退或尚未配套,导致寿光市至下游段弥河道洪水漫溢。模型分析表明,弥河寿光下游排洪能力可能仅仅在1500m3/s左右,不足寿光市防洪预案II级橙色预警要求。
- 黑虎山水库可能前期水位略有偏高,导致雨峰来临时迅速下泄。冶源水库实现了下泄错峰并拦截了近1000m3/s的入库洪峰,防洪效果显著。尽管仍有2000多万m3可用调洪库容,在综合考虑水库大坝重大设施安全风险和决策紧迫的条件下,认为调度得当。
- 上游一大二中三座水库819洪水的优化调度略有潜力。但可能由于目前流域水文和洪水预报技术的局限,对水情和风险形势判断的准确性和自信心稍有不足。尽管这不是本次寿光洪水至灾主要因素,但为发挥水库综合效益,在水资源利用、防洪风险管理决策、支撑下游社会经济发展及城市安全方面,决策技术的提高已经非常必要,也很迫切。