一种确定梯级水库群连续溃坝风险分析方法技术

本发明专利技术公开了一种梯级水库群连续溃坝风险分析方法。该方法包括溃口洪水计算、溃坝洪水演进、漫坝风险分析三大模块，从风险量化与预警预防角度，提出了解决梯级水库群连续溃坝洪水形成、演进及其风险分析的技术方法。溃口洪水计算采用宽顶堰公式估算溃口流量过程，采用双曲线模型模拟溃口冲刷侵蚀过程，采用简化的Bishop法计算溃口扩展过程；溃坝洪水演进基于HLL格式的数值差分方法，采用中心格式的有限体积法进行计算；漫坝风险分析采用基于时间微元的方法，逐时段地联立求解水量平衡方程和水库的蓄泄方程。本发明专利技术能够快速分析评估流域梯级水库群中的堰塞坝、土石坝、混凝土坝等连续漫顶或溃坝的风险，为提高应急抢险工作效率，快速制定抢险处置方案提供技术支撑。

【技术实现步骤摘要】
一种确定梯级水库群连续溃坝风险分析方法
本专利技术涉及水利水电工程领域，是一种梯级水库群连续溃坝风险分析方法，是一种流域连续溃坝洪水风险分析和数值计算技术，可为梯级水工建筑物风险设计、梯级水库群安全管理、流域中堰塞湖应急处置、梯级水电站应急管理、防洪度汛方案制定提供定量的技术支撑。
技术介绍
梯级水库群是河流水资源开发的基本形式，更是大规模水能资源开发的必然选择。我国已建成各类水库大坝超过9.8万座，形成一系列水库群。这些水库群在防洪减灾、提供清洁电力、保障生态安全以及节能减排等方面，发挥了不可替代的作用，但如若由于不可抗力导致梯级连溃事故，其后果损失也是极其惨重。目前，单一梯级溃坝的研究成果较多，但分析多座梯级水库群连续溃坝风险的研究成果较少。暴雨洪水、地震、滑坡赌江等灾害事件对流域梯级水库群构成极大的连续溃决风险，如2014年牛栏江流域的红石岩堰塞湖、2018年金沙江上游的白格堰塞湖，均对下游梯级水库群或水电站造成极大威胁。现有溃坝数值计算方法，多采用简单的数学模型模拟溃坝过程，未考虑溃坝的物理机制，无法更切合实际地模拟溃坝洪水过程，且未形成系统的梯级水库群连续溃坝风险分析方法，无法快速分析下游梯级水库漫坝或溃坝风险；或而采用复杂的高性能数值计算方法计算溃坝洪水演进，对地形、河道形态等基础数据精度要求较高，且计算耗时长，无法满足流域应急抢险工作的实际需要。
技术实现思路
本专利技术根据流域梯级水库群应急抢险工作需要和现有技术的不足，提供一种梯级水库群连续溃坝风险分析方法，该方法集成溃口洪水计算、溃坝洪水演进、漫坝风险分析三大模块，从风险量化与预警预防角度，提出了解决梯级水库群连续溃坝风险分析的技术方法，为快速制定抢险处置方案，提高应急抢险工作效率，确保处置科学有效提供了技术支撑。为实现上述目的，本专利技术提供如下技术方案：一种梯级水库群连续溃坝风险分析方法，所述梯级水库群包括至少两座或以上水库大坝，所述梯级水库群连续溃坝风险分析方法步骤如下：1)选定连续溃坝风险分析对象。根据本专利技术需在同一河流自上而下依次选定至少两座或以上梯级水库大坝为连续溃坝风险分析对象。2)基础数据收集与确定。根据本专利技术计算分析要求，梯级水库群连续溃坝风险分析需要收集梯级水库库容—水位关系曲线和数据、坝体几何特征参数、入库流量、起调水库水位、泄洪能力曲线和数据，相邻梯级水库之间的河道坡降、糙率、几何形状参数，确定溃口冲刷侵蚀参数、溃口扩展参数、洪水演进计算时间步长等。3)溃口洪水计算。根据本专利技术提出的溃口水流模型、溃口冲刷侵蚀模型和溃口几何扩展模型，确定溃口洪水流量过程和峰值流量。4)溃坝洪水演进。根据本专利技术提出的基于HLL格式的数值差分方法，采用中心格式的有限体积法，进行溃坝洪水演进计算，确定洪水演进至下游梯级水库的流量过程和相应的时间。5)漫坝风险分析。根据本专利技术集成的调洪分析模型，计算分析是否存在漫顶连续溃坝的可能。上述技术方案中，优选地，3)溃口洪水计算的溃口水流模型是在水量平衡的基础上，采用宽顶堰公式进行估算，即式中B为溃口断面的宽度；H为库水位高程；z为溃口进口处床面高程；C为流量系数；q为水库自然入流量；W为库容量，V是自变量为库水位高H的函数。上述技术方案中，优选地，3)溃口洪水计算的溃口冲刷侵蚀采用双曲线模型，即式中dz/dt为溃口冲刷侵蚀率；τ为剪应力；τc为临界剪应力；k为单位变换因子；a、b为双曲线参数。该模型认为水流冲刷土石料时，土石料抵抗冲刷侵蚀的能力不应是无限制的，而应具有一定的“强度”，即双曲线有一当τ＝τc时的渐进线，即dz/dt的极值1/b。上述技术方案中，优选地，3)溃口洪水计算的溃口几何扩展过程采用简化的简化Bishop法进行计算。上述技术方案中，优选地，4)溃坝洪水演进计算过程中，为提高计算效率，将河道实测断面概化为倒梯形，断面面积A和水面宽度B存在如下关系：A＝h(B0+hm)，B＝B0+2hm，其中B0为河道底宽，h为水流深度，m为河道两侧坡比。上述技术方案中，优选地，4)溃坝洪水演进计算过程中，上游边界条件为上游梯级溃坝的流量过程或上游来水流量过程；下游边界条件根据下游断面水位流量关系，采用曼宁公式即可确定。演进过程中，各断面的起始条件，一般选择天然径流量。具体边界条件如下：(1)初始条件：初始流场Φt＝0＝Φ0(x，y)，其Φ分别为x、y方向的水位(Z)、流速(V)、库容(W)关于时间(t)的函数，即：t＝t0，Z＝Z0，V＝V0，W＝W0；(2)水位边界：Z＝Z(t)；(3)流量边界：Q＝Q(t)；(4)水位流量关系边界：Q＝Q(Z)。在以上(3)和(4)中，根据边界上流速分布，即Q＝∫sVnhds，可得边界节点上法向流速。此外，h采用nΔt时刻的节点水深。上述技术方案中，优选地，5)漫坝风险分析采用基于时间微元的方法，调洪漫坝风险计算的实用方程即为逐时段地联立求解水量平衡方程和水库的蓄泄方程，即式中Q1、Q2分别为计算时段初、末的入库流量，q1、q2分别为计算时段初、末的下泄流量，W1、W2分别为计算时段初、末水库的需水量，ΔW为库容变化量，Δt为计算时段。与现有技术相比，本专利技术的有益效果是：本专利技术提供了一种简便易行的梯级水库群连续溃坝风险分析方法，可为流域梯级水库群安全管理、应急抢险、风险处置决策提供准确、快捷和定量化的技术支撑。附图说明图1梯级水库群连续溃坝风险分析流程图图2梯级水库群平面布置图图3梯级水库群位置示意图图4各梯级库容水位曲线图5C库泄流能力曲线图6A坝溃决并演进至B库流量过程及B库水位变化过程图7B坝溃决及演进至C库流量过程图8C库预先泄洪过程(0:00～11:38)图9C库漫坝风险分析过程(11:38～次日6:00)具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚完整的描述。应该了解，下面所描述的实施例仅仅是本专利技术一部分，而非全部的实施例。基于本专利技术中的实施例，在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。1)选定连续溃坝风险分析对象以某河流干流上自上而下依次规划的A、B、C三座梯级水库群为实施例，平面布置和位置示意见图2和图3。C坝为三座梯级所在河段的控制性梯级。2)基础数据收集与确定A坝距离B坝约85.6km，高差约318.9m，B坝距离C坝约89.3km，高差约380.5m。A、B、C三座梯级工程特性参数见表1，库容水位曲线见图4。B坝水库不具有调洪能力，不考虑其泄洪过程。C坝水库按设计方案具备防洪功能，泄洪设施包括溢洪道、竖井泄洪洞、深孔泄洪洞，全闸泄洪流量约8000m3/s，泄流能力见图5。表1三座梯级水库工程特性表3)A库溃坝洪水计算A库溃坝参数取值见表2。假定A坝于0:00开本文档来自技高网...

【技术保护点】
1.一种确定梯级水库群连续溃坝风险分析方法，所述梯级水库群包括两座以上水库大坝，其特征在于，所述梯级水库群连续溃坝风险分析方法步骤如下：/n1)选定连续溃坝风险分析对象：需在同一河流自上而下依次选定至少两座以上梯级水库大坝为连续溃坝风险分析对象；/n2)基础数据收集与确定：根据计算要求，梯级水库群连续溃坝风险分析需要收集梯级水库库容水位关系曲线和数据、大坝主要特征参数、入库洪水流量、水库起调水位、泄洪能力曲线和数据，相邻梯级水库之间的河道坡降、糙率、几何形状参数，确定溃口冲刷侵蚀参数、溃口扩展参数、洪水演进计算时间步长等；/n3)溃口洪水计算：根据本专利技术提出的溃口水流泄流模型、溃口冲刷侵蚀模型和溃口几何扩展模型，确定溃口洪水流量过程和峰值流量；/n4)溃坝洪水演进：基于HLL格式的数值差分方法，采用中心格式的有限体积法，进行溃坝洪水演进计算，确定洪水演进至下游梯级水库的流量过程和相应的时间；/n5)漫坝风险分析：根据集成的调洪分析模型，计算分析是否存在漫顶连续溃坝的可能。/n

【技术特征摘要】
1.一种确定梯级水库群连续溃坝风险分析方法，所述梯级水库群包括两座以上水库大坝，其特征在于，所述梯级水库群连续溃坝风险分析方法步骤如下：
1)选定连续溃坝风险分析对象：需在同一河流自上而下依次选定至少两座以上梯级水库大坝为连续溃坝风险分析对象；
2)基础数据收集与确定：根据计算要求，梯级水库群连续溃坝风险分析需要收集梯级水库库容水位关系曲线和数据、大坝主要特征参数、入库洪水流量、水库起调水位、泄洪能力曲线和数据，相邻梯级水库之间的河道坡降、糙率、几何形状参数，确定溃口冲刷侵蚀参数、溃口扩展参数、洪水演进计算时间步长等；
3)溃口洪水计算：根据本发明提出的溃口水流泄流模型、溃口冲刷侵蚀模型和溃口几何扩展模型，确定溃口洪水流量过程和峰值流量；
4)溃坝洪水演进：基于HLL格式的数值差分方法，采用中心格式的有限体积法，进行溃坝洪水演进计算，确定洪水演进至下游梯级水库的流量过程和相应的时间；
5)漫坝风险分析：根据集成的调洪分析模型，计算分析是否存在漫顶连续溃坝的可能。


2.根据权利要求1所述的确定梯级水库群连续溃坝风险分析方法，其特征在于，所述溃口洪水计算的溃口水流模型是在水量平衡的基础上，采用如下公式进行估算，即



上式中B为溃口断面的宽度；H为库水位高程；z为溃口进口处床面高程；C为流量系数；q为水库自然入流量；W为库容量，V是自变量为库水位高H的函数。


3.根据权利要求1所述的确定梯级水库群连续溃坝风险分析方法，其特征在于，所述溃口洪水计算的溃口冲刷侵蚀采用双曲线模型，即

式中dz/dt为溃口冲刷侵蚀率；τ为剪应力；τc为临界剪应力；k为单位变换因子；a、b为双曲线参数；该模型...

【专利技术属性】
技术研发人员：周建平，周兴波，陈祖煜，杜效鹄，于沭，王琳，王双敬，高洁，张强，
申请(专利权)人：中国电力建设股份有限公司，水电水利规划设计总院，中国水利水电科学研究院，
类型：发明
国别省市：北京;11