CN110362925A - 一种包含库区的土石坝漫顶溃决洪水数值模拟方法 - Google Patents
一种包含库区的土石坝漫顶溃决洪水数值模拟方法 Download PDFInfo
- Publication number
- CN110362925A CN110362925A CN201910640947.1A CN201910640947A CN110362925A CN 110362925 A CN110362925 A CN 110362925A CN 201910640947 A CN201910640947 A CN 201910640947A CN 110362925 A CN110362925 A CN 110362925A
- Authority
- CN
- China
- Prior art keywords
- dam
- crevasse
- area
- reservoir area
- unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Abstract
本发明公开了一种包含库区的土石坝漫顶溃决洪水数值模拟方法。采用三角形非结构网格对库区和下游淹没区进行剖分,坝体所在区域不参与剖分。将溃口出流量作为库区的出流边界和下游淹没区的入流边界,分别进行库区和下游淹没区的水动力演进计算,出流边界和入流边界所对应的坝前单元和坝后单元根据溃口的实时宽度动态调整。该方法不仅能够提供下游的淹没形态,而且能够提供库区内的水流运动过程,从而扩大溃坝模型的适用范围。在计算溃口出流时,使用跟溃口直接相连的坝前单元的平均水位,而不是整个库区的平均水位,可显著提高溃口出流过程的准确性,从而完善现有的土石坝溃决洪水计算方法,提高土石坝溃决洪水计算的数值精度。
Description
技术领域
本发明涉及水利工程领域,尤其涉及防洪减灾领域,具体为一种包含库区的土石坝漫顶溃决洪水数值模拟方法。
背景技术
我国是世界上修建水库大坝最多的国家,根据第一次全国水利普查工作统计,我国现有水库大坝98002座,其中土石坝作为应用最广泛的坝型,约占已建坝总数的95%以上。已经修建的土石坝中,绝大多数为小型水库,由于特殊的历史原因,这些小型土石坝在修建时多缺乏严格的设计和施工程序,当遭遇超标准洪水时,存在比较较大的溃决风险。水库大坝一旦溃决失事,将会给下游带来巨大的灾难。在土石坝的溃决失事方式中,最常见的为洪水漫顶溃决,洪水一旦漫顶,会迅速在坝后形成冲刷,直至溃口稳定。另外值得说明的是,土石坝溃决洪峰流量远大于水库的正常泄洪流量,所以在溃坝分析计算时,一般可忽略溢洪道泄洪或底孔泄洪的影响。
对土石坝溃决洪水分析计算来讲,非常重要的一步是计算溃口的出流过。溃口出流计算的方法总体来说可以分为两大类:一类为采用溃坝的历史数据资料得出的回归方程来计算出流过程,这类回归过程一般是建立溃口出流过程与库容和坝高的相关关系的表达式。另一类为采用带有溃口发展物理机制的理论或半理论方法来预测出流过程,在这其中又可分为两小类,一类为参数型溃口模型,另一类为物理-参数型溃口模型。所谓参数型溃口模型是指通过设定溃口尺寸,形状和发展时间,利用堰流公式来计算溃坝洪水出流过程,这些模型通常采用相对简单的参数化经验方法来估算溃口发展过程。所谓物理-参数型溃口模型指的是溃口的发展考虑了溃口冲刷的物理机制,通过泥沙输移公式来计算溃口侵蚀过程。这类模型相对较多,代表了土石坝溃口出流模型的一个发展方向。
现在的土石坝溃决计算实际上人为分成了两个步骤,首先是溃口出流过程的计算,得到完整的溃口出流过程后,再进行下一步洪水演进计算。然而,真实的溃决过程是库区水位降低、溃口发展和下游洪水演进是同时进行的,也是相互影响的,溃口前局部水位的高低直接影响着溃口的发展,溃口的发展决定了溃口出流量的大小,出流量的大小决定了下游淹没形势的不同,因此,将库区包含在模型内,进行库区、溃口和下游淹没区一体化的模拟是非常有必要的,能够提高当前土石坝溃决洪水的模拟精度。
发明内容
本发明的目的在于提供一种包含库区的土石坝漫顶溃决数值模拟方法,在每个时间步长内同时完成库区洪水演进计算、溃口出流计算和下游淹没过程的计算,充分考虑土石坝溃决过程的整体性,提高现有土石坝溃决洪水分析计算的精度。
本发明是通过以下技术方案实现的:
将库区包含在整个计算区域内,土石坝溃口出流过程的模拟与库区和下游淹没区洪水演进模拟统一进行,具体步骤如下:
1)获取所需的基础数据:所述基础数据包括库区水下地形数据(比例尺1:5000及以上);坝体相关数据,包含坝高、坝长、坝顶宽度、迎水面和背水面坡度、坝体材料粘性系数、粒径级配等;溃坝淹没区地形数据(比例尺1:10000及以上);获取整个计算区域DOM地貌数据(分辨率2m及以上);入库洪水过程数据。
2)计算区域网格离散:库区与坝体相连接的边定义为坝前边,下游淹没区与坝体相连接的边称为坝后边,坝前边与坝后边长度相等,相互对应;坝前边和坝后边作为网格剖分的控制边,离散方式相同,离散尺度控制在5m左右;采用三角形非结构网格对库区和下游淹没区分别离散,坝体所在的区域不参与网格离散,库区和下游淹没区网格单元统一编码;离散后,坝前边对应的网格单元称为坝前单元,坝后边对应的网格单元为坝后单元,坝前单元和坝后单元数量相同,一一对应。
3)计算区域赋值:采用库区和下游淹没区地形数据对网格单元型心位置进行高程插值;根据地貌数据对网格单元糙率赋值;对网格单元的水力变量赋初值;将入库洪水过程设置为库区的入流边界;下游淹没区的出流边界设置为自由出流。
4)溃口设置:设定溃口形状为矩形,溃口初始宽度设置为0.5m,初始深度设置为0.2m;不失一般性,溃口初始位置设置在坝前边的中间位置,此时溃口对应的坝前单元编号为U0,以U0单元为中心,顺水流方向左侧坝前单元编号依次为U1,U3,U5……Un,其中n为奇数,,顺水流方向右侧坝前单元编号依次为U2,U4,U6……Um,其中m为偶数;同样的,与溃口初始位置对应的坝后单元编号为D0,以D0单元为中心,顺水流左侧坝后单元编号依次为D1,D3,D5……Dn,其中n为奇数,顺水流右侧坝后单元编号依次为D2,D4,D6……Dm,其中m为偶数;按照各坝前单元(坝后单元)的编号顺序,与各坝前单元(坝后单元)对应的坝前边(坝后边)的长度依次定义为B0,B1,B2,B3,B4,B5,B6……Bi,i为坝前边或坝后边离散后线段的总条数,i为自然数。
5)确定时间步长dt:由各网格单元的水力变量值,获取满足计算稳定CFL条件的时间步长dt。
6)溃口出流过程计算:根据溃口出流公式计算t时刻溃口的出流量,公式如下式所示:
Qb=kbB(Zu-Zd)1.5 (1)
其中:Qb为溃口出流量,Kb为流量系数,B为溃口宽度,Zu为与溃口直接相连的所有坝前单元的平均水位值,Zd为溃口底部高程。
根据de Vries输沙率公式来计算溃口处的泥沙输移量,以确定下一时刻t+dt的溃口宽度和深度,计算溃口发展。
7)洪水演进过程计算:洪水演进计算包括库区洪水演进计算和下游淹没区洪水演进计算,采用守恒型的二维浅水方程组描述洪水运动;使用溃口流量Qb驱动库区和下游淹没区的洪水演进计算,在库区洪水演进计算中Qb作为出流边界处理,在下游淹没区洪水演进计算中Qb作为入流边界处理;
根据溃口在当前时刻的宽度B来确定出流边界所对应的坝前单元和入流边界对应的坝后单元。在选择坝前单元时,,如果B小于B0,出流边界所对应的坝前单元为U0;当B大于B0时,以U0为中心,每一次在原来已选择的坝前单元的基础上再向两侧各选取1个坝前单元,直到所选的所有坝前单元对应的坝前边之和大于等于溃口宽度B时结束,这些选中的坝前单元作为当前时刻出流边界所在单元参与本时段溃口流量Qb的分配;同理,与这些选中的坝前单元对应的坝后单元作为当前时刻入流边界所在单元参与本时段溃口流量Qb的分配。
8)更新t=t+dt,重复步骤5)~7)直至计算结束。
进一步的,步骤6)中溃口发展的计算公式如下:
其中,dZd为时段dt内溃口的冲刷深度,cs和βs为输沙率公式系数,Qb为溃口出流量,Ab为溃口的过水断面面积。设定溃口的冲刷深度与侧向展宽的宽度相同,溃口向两侧对称展宽。
进一步的,步骤7)中,描述洪水运动所采用的方程组为守恒型二维浅水方程组,采用具有良好激波捕捉能力的Godunov格式离散该方程组。采用的守恒型二维浅水方程组如下式所示:
式中:
h为水深,u,v分别为x,y方向的流速,t为时间,分别x,y方向的坡度,Zb为地面高程,g为重力加速度,分别为x,y方向的摩阻项,其中n为Manning糙率系数。
本发明的有益效果:
在进行土石坝漫顶溃决模拟时,通过将溃口出流过程分别作为上下游区域的出流和入流边界条件的方法将库区、溃口和下游淹没区统一起来建立一体化的土石坝溃决模拟模型,该方法不仅能够提供下游的淹没形态,而且能够提供库区内的水流运动过程,在计算溃口出流时,使用溃口前的局部水位,而不是整个库区的平均水位,可显著提高溃口出流过程的准确性,从而整体上完善现有的土石坝溃决洪水计算方法,提高土石坝溃决洪水计算的数值精度。
附图说明
图1为本发明的一种包含库区的土石坝漫顶溃决洪水数值模拟方法流程图;
图2为本发明的一种包含库区的土石坝漫顶溃决洪水数值模拟方法示意图。
具体实施方式
下面结合附图1和附图2对本发明作进一步说明。
本发明为一种包含库区的土石坝漫顶溃决洪水数值模拟方法,将库区包含在整个计算区域内,土石坝溃口出流过程的模拟与库区和下游淹没区洪水演进模拟统一进行,具体步骤如下:
1)获取所需的基础数据:所述基础数据包括库区水下地形数据(比例尺1:5000及以上);坝体相关数据,包含坝高、坝长、坝顶宽度、迎水面和背水面坡度、坝体材料粘性系数、粒径级配等;溃坝淹没区地形数据(比例尺1:10000及以上);获取整个计算区域DOM地貌数据(分辨率2m及以上),入库洪水过程数据。
2)计算区域网格离散:库区与坝体相连接的边定义为坝前边,下游淹没区与坝体相连接的边称为坝后边,坝前边与坝后边长度相等,相互对应;坝前边和坝后边作为网格剖分的控制边,离散方式相同,离散尺度控制在5m左右;采用三角形非结构网格对库区和下游淹没区分别离散,坝体所在的区域不参与网格离散,库区和下游淹没区网格单元统一编码;离散后,坝前边对应的网格单元称为坝前单元,坝后边对应的网格单元为坝后单元,坝前单元和坝后单元数量相同,一一对应。
3)计算区域赋值:采用库区和下游淹没区地形数据对网格单元型心位置进行高程插值;根据地貌数据对网格单元糙率赋值;对网格单元的水力变量赋初值;将入库洪水过程设置为库区的入流边界;下游淹没区的出流边界设置为自由出流。
4)溃口设置:设定溃口形状为矩形,溃口初始宽度设置为0.5m,初始深度设置为0.2m;不失一般性,溃口初始位置设置在坝前边的中间位置,此时溃口对应的坝前单元编号为U0,以U0单元为中心,顺水流方向左侧坝前单元编号依次为U1,U3,U5……Un,其中n为奇数,顺水流方向右侧坝前单元编号依次为U2,U4,U6……Um,其中m为偶数;同样的,与溃口初始位置对应的坝后单元编号为D0,以D0单元为中心,顺水流左侧坝后单元编号依次为D1,D3,D5……Dn,其中n为奇数,,顺水流右侧坝后单元编号依次为D2,D4,D6……Dm,其中m为偶数;按照各坝前单元(坝后单元)的编号顺序,与各坝前单元(坝后单元)对应的坝前边(坝后边)的长度依次定义为B0,B1,B2,B3,B4,B5,B6……Bi,i为坝前边或坝后边离散后线段的总条数,i为自然数。
5)确定时间步长dt:由各网格单元的水力变量值,获取满足计算稳定CFL条件的时间步长dt。
6)溃口出流过程计算:根据溃口出流公式计算t时刻溃口的出流量,公式如下式所示:
Qb=kbB(Zu-Zd)1.5 (1)
其中:Qb为溃口出流量,Kb为流量系数,B为溃口宽度,Zu为与溃口直接相连的所有坝前单元的平均水位值,以图2为例,此时的Zu为U0,U1和U2三个坝前单元水位的平均值,Zd为溃口底部高程。
根据de Vries输沙率公式来计算溃口处的泥沙输移量,以确定下一时刻t+dt的溃口宽度和深度,溃口发展的计算公式如下:
其中,dZd为时段dt内溃口的冲刷深度,cs和βs为输沙率公式系数,Qb为溃口出流量,Ab为溃口的过水断面面积。设定溃口的冲刷深度与侧向展宽的宽度相同,溃口向两侧对称展宽。当达到设定的溃口展宽宽度后,溃口发展达到稳定状态,溃口的形状不在发生变化。
7)洪水演进过程计算:洪水演进过程计算包括库区洪水演进计算和下游淹没区洪水演进计算,采用守恒型的二维浅水方程组描述洪水运动。守恒型二维浅水方程组如下式所示:
式中:h为水深,u,v分别为x,y方向的流速,t为时间,分别x,y方向的坡度,Zb为地面高程,g为重力加速度,分别为x,y方向的摩阻项,其中n为Manning糙率系数。
对以上浅水方程组进行非结构离散,采用Roe格式的近似Riemann解计算通过非结构单元界面的数值通量,对底坡项采用特征分解的方法处理,摩阻项采用半隐的方法处理。详细的数值求解过程可参见如下著作P85-P90页(张大伟.基于Godunov格式的堤坝溃决水流数值模拟,2014,中国水利水电出版社)。
使用溃口流量Qb驱动库区和下游淹没区的洪水演进计算,在库区洪水演进中Qb作为出流边界处理,在下游淹没区洪水演进计算中Qb作为入流边界处理。根据溃口在当前时刻的宽度B来确定出流边界所对应的坝前单元和入流边界对应的坝后单元。在选择坝前单元时,如果B小于B0,出流边界所对应的坝前单元为U0;当B大于B0时,以U0为中心,每一次在原来已选择的坝前单元的基础上再向两侧各选取1个坝前单元,直到所选的所有坝前单元对应的坝前边之和大于等于溃口宽度B时结束,如图2所示,此时的溃口宽度B明显大于B0,则以U0为中心,选取U1和U2后,B0+B1+B2之和已大于B,选取过程结束,选中的坝前单元为U0,U1和U2,这些选中的坝前单元作为当前时刻出流边界所在单元参与本时段溃口流量Qb的分配;同理,与这些选中的坝前单元对应的坝后单元作为当前时刻入流边界所在单元参与本时段溃口流量Qb的分配。详细的Qb分配处理方法可参见如下著作P92页(张大伟.基于Godunov格式的堤坝溃决水流数值模拟,2014,中国水利水电出版社)。
8)更新t=t+dt,重复步骤5)~7)直至计算结束。
上述的实施例仅是本发明的部分体现,并不能涵盖本发明的全部,在上述实施例以及附图的基础上,本领域技术人员在不付出创造性劳动的前提下可获得更多的实施方式,因此这些不付出创造性劳动的前提下获得的实施方式均应包含在本发明的保护范围内。
Claims (3)
1.一种包含库区的土石坝漫顶溃决洪水数值模拟方法,其特征在于:将库区包含在整个计算区域内,土石坝溃口出流过程的模拟与库区和下游淹没区洪水演进模拟统一进行,具体步骤如下:
1)获取所需的基础数据:所述基础数据包括库区水下地形数据、坝体相关数据、溃坝淹没区地形数据、整个计算区域DOM地貌数据、入库洪水过程数据;所述坝体相关数据包括坝高、坝长、坝顶宽度、迎水面和背水面坡度、坝体材料粘性系数、粒径级配;
2)计算区域网格离散:库区与坝体相连接的边定义为坝前边,下游淹没区与坝体相连接的边定义为坝后边,坝前边与坝后边长度相等,相互对应;坝前边和坝后边作为网格剖分的控制边,坝前边与坝后边的离散方式相同,离散尺度控制在5m;采用三角形非结构网格对库区和下游淹没区分别离散,坝体所在的区域不参与网格离散,库区和下游淹没区网格单元统一编码;离散后,坝前边对应的网格单元称为坝前单元,坝后边对应的网格单元称为坝后单元,坝前单元和坝后单元数量相同,一一对应;
3)计算区域赋值:采用库区和下游淹没区地形数据对网格单元型心位置进行高程插值;根据地貌数据对网格单元糙率赋值;对网格单元的水力变量赋初值;将入库洪水过程设置为库区的入流边界;下游淹没区的出流边界设置为自由出流;
4)溃口设置:设定溃口形状为矩形,溃口初始宽度设置为0.5m,初始深度设置为0.2m;不失一般性,溃口初始位置设置在坝前边的中间位置,此时溃口对应的坝前单元编号为U0,以U0单元为中心,顺水流方向左侧坝前单元编号依次为U1,U3,U5……Un,其中n为奇数,顺水流方向右侧坝前单元编号依次为U2,U4,U6……Um,其中m为偶数;与溃口初始位置对应的坝后单元编号为D0,以D0单元为中心,顺水流左侧坝后单元编号依次为D1,D3,D5……Dn,其中n为奇数,顺水流右侧坝后单元编号依次为D2,D4,D6……Dm,其中m为偶数;按照各坝前单元及坝后单元的编号顺序,与各坝前单元及坝后单元分别对应的坝前边及坝后边的长度依次定义为B0,B1,B2,B3,B4,B5,B6……Bi,i为坝前边或坝后边离散后线段的总条数,i为自然数;
5)确定时间步长dt:由各网格单元的水力变量值,获取满足计算稳定CFL条件的时间步长dt;
6)溃口出流过程计算:根据溃口出流公式计算t时刻溃口的出流量,溃口出流公式如式(1)所示:
Qb=kbB(Zu-Zd)1.5 (1)
其中:Qb为溃口出流量,kb为流量系数,B为溃口宽度,Zu为与溃口直接相连的所有坝前单元的平均水位值,Zd为溃口底部高程;
根据de Vries输沙率公式来计算溃口处的泥沙输移量,以确定下一时刻t+dt的溃口宽度和深度,计算溃口发展;
7)洪水演进过程计算:洪水演进计算包括库区洪水演进计算和下游淹没区洪水演进计算,采用守恒型的二维浅水方程组描述洪水运动;使用溃口流量Qb驱动库区和下游淹没区的洪水演进计算,在库区洪水演进计算中Qb作为出流边界处理,在下游淹没区洪水演进计算中Qb作为入流边界处理;
根据溃口在当前时刻的宽度B来确定出流边界所对应的坝前单元和入流边界对应的坝后单元:在选择坝前单元时,如果B小于B0,出流边界所对应的坝前单元为U0;当B大于B0时,以U0为中心,每一次在原来已选择的坝前单元的基础上再向两侧各选取1个坝前单元,直到所选的所有坝前单元对应的坝前边之和大于等于溃口宽度B时结束,选中的坝前单元作为当前时刻出流边界所在单元参与本时段溃口流量Qb的分配;与选中的坝前单元对应的坝后单元作为当前时刻入流边界所在单元参与本时段溃口流量Qb的分配;
8)更新t=t+dt,重复步骤5)~7)直至计算结束。
2.一种包含库区的土石坝漫顶溃决洪水数值模拟方法,其特征在于:步骤6)中溃口发展的计算公式如下:
其中,dZd为时段dt内溃口的冲刷深度,cs和βs为输沙率公式系数,Qb为溃口出流量,Ab为溃口的过水断面面积;设定溃口的冲刷深度与侧向展宽的宽度相同,溃口向两侧对称展宽。
3.一种包含库区的土石坝漫顶溃决洪水数值模拟方法,其特征在于:步骤7)中,描述洪水运动所采用的方程组为守恒型二维浅水方程组,采用具有良好激波捕捉能力的Godunov格式离散该方程组,采用的守恒型二维浅水方程组如下式所示:
式中:
h为水深,u,v分别为x,y方向的流速,t为时间,分别x,y方向的坡度,Zb为地面高程,g为重力加速度,分别为x,y方向的摩阻项,其中n为Manning糙率系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910640947.1A CN110362925B (zh) | 2019-07-16 | 2019-07-16 | 一种包含库区的土石坝漫顶溃决洪水数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910640947.1A CN110362925B (zh) | 2019-07-16 | 2019-07-16 | 一种包含库区的土石坝漫顶溃决洪水数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110362925A true CN110362925A (zh) | 2019-10-22 |
CN110362925B CN110362925B (zh) | 2020-05-19 |
Family
ID=68219591
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910640947.1A Active CN110362925B (zh) | 2019-07-16 | 2019-07-16 | 一种包含库区的土石坝漫顶溃决洪水数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110362925B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110485363A (zh) * | 2019-08-31 | 2019-11-22 | 华北水利水电大学 | 洪水作用下河床动态变化的模拟装置及多点实时监测方法 |
CN111241757A (zh) * | 2020-01-10 | 2020-06-05 | 中核第四研究设计工程有限公司 | 基于计算流体力学的铀尾矿库溃坝三维数值模拟方法 |
CN112697673A (zh) * | 2020-12-14 | 2021-04-23 | 中国水利水电科学研究院 | 一种穿堤无压涵管接触渗流破坏的可视化试验装置与方法 |
CN113111489A (zh) * | 2021-03-11 | 2021-07-13 | 重庆交通大学 | 一种堰塞坝洪水漫顶溃坝过程仿真方法及仿真系统 |
CN113297753A (zh) * | 2021-06-24 | 2021-08-24 | 四川大学 | 山区河流透过性坝沿程水深预测方法 |
CN115455867A (zh) * | 2022-10-31 | 2022-12-09 | 武汉大学 | 基于回归分析的坝区流态的推求方法 |
CN116757125A (zh) * | 2023-08-18 | 2023-09-15 | 水利部交通运输部国家能源局南京水利科学研究院 | 基于sph-dem算法的土石坝漫顶溃决模拟方法 |
CN117034801A (zh) * | 2023-08-07 | 2023-11-10 | 中国水利水电科学研究院 | 一种基于河道与防洪保护区一体化的溃堤洪水模拟方法 |
CN117494478A (zh) * | 2023-12-29 | 2024-02-02 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种模拟心墙坝漫顶溃坝流量及溃口演化过程的计算方法 |
CN117521404A (zh) * | 2023-11-30 | 2024-02-06 | 中国水利水电科学研究院 | 一种平原水库溃坝洪水危险等级划分方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107451372A (zh) * | 2017-08-09 | 2017-12-08 | 中国水利水电科学研究院 | 一种运动波与动力波相结合的山区洪水过程数值模拟方法 |
CN108256193A (zh) * | 2018-01-10 | 2018-07-06 | 中国水利水电科学研究院 | 一种河道溃决洪水一、二维数学模型耦合分析方法 |
CN108446502A (zh) * | 2018-03-22 | 2018-08-24 | 中国水利水电科学研究院 | 一种利用完整二维浅水方程组获得流域单位线的方法 |
KR101906858B1 (ko) * | 2017-06-28 | 2018-10-12 | 홍익대학교 산학협력단 | 수치모형 및 침수선 보간을 이용한 하천 범람 예측 방법 |
CN109284531A (zh) * | 2018-08-03 | 2019-01-29 | 中国水利水电科学研究院 | 一种基于空间拓扑的一二维水动力学耦合方法 |
CN109543275A (zh) * | 2018-11-15 | 2019-03-29 | 中国水利水电科学研究院 | 一种城区地表径流二维数值模拟方法 |
CN109815574A (zh) * | 2019-01-16 | 2019-05-28 | 中国恩菲工程技术有限公司 | 尾矿库漫顶溃坝的洪水条件确定方法及装置 |
CN109918821A (zh) * | 2019-03-15 | 2019-06-21 | 中国水利水电科学研究院 | 一种迎风守恒型河道漫溢出流数值模拟方法 |
CN109948176A (zh) * | 2019-01-16 | 2019-06-28 | 中国恩菲工程技术有限公司 | 尾矿坝安全评价方法、存储介质、电子设备 |
-
2019
- 2019-07-16 CN CN201910640947.1A patent/CN110362925B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101906858B1 (ko) * | 2017-06-28 | 2018-10-12 | 홍익대학교 산학협력단 | 수치모형 및 침수선 보간을 이용한 하천 범람 예측 방법 |
CN107451372A (zh) * | 2017-08-09 | 2017-12-08 | 中国水利水电科学研究院 | 一种运动波与动力波相结合的山区洪水过程数值模拟方法 |
CN108256193A (zh) * | 2018-01-10 | 2018-07-06 | 中国水利水电科学研究院 | 一种河道溃决洪水一、二维数学模型耦合分析方法 |
CN108446502A (zh) * | 2018-03-22 | 2018-08-24 | 中国水利水电科学研究院 | 一种利用完整二维浅水方程组获得流域单位线的方法 |
CN109284531A (zh) * | 2018-08-03 | 2019-01-29 | 中国水利水电科学研究院 | 一种基于空间拓扑的一二维水动力学耦合方法 |
CN109543275A (zh) * | 2018-11-15 | 2019-03-29 | 中国水利水电科学研究院 | 一种城区地表径流二维数值模拟方法 |
CN109815574A (zh) * | 2019-01-16 | 2019-05-28 | 中国恩菲工程技术有限公司 | 尾矿库漫顶溃坝的洪水条件确定方法及装置 |
CN109948176A (zh) * | 2019-01-16 | 2019-06-28 | 中国恩菲工程技术有限公司 | 尾矿坝安全评价方法、存储介质、电子设备 |
CN109918821A (zh) * | 2019-03-15 | 2019-06-21 | 中国水利水电科学研究院 | 一种迎风守恒型河道漫溢出流数值模拟方法 |
Non-Patent Citations (2)
Title |
---|
张大伟 等: "堰塞坝漫顶溃决试验及相关数学模型研究", 《水利学报》 * |
陈生水 等: "土石坝漫顶破坏溃口发展数值模型研究", 《水利水运工程学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110485363A (zh) * | 2019-08-31 | 2019-11-22 | 华北水利水电大学 | 洪水作用下河床动态变化的模拟装置及多点实时监测方法 |
CN111241757A (zh) * | 2020-01-10 | 2020-06-05 | 中核第四研究设计工程有限公司 | 基于计算流体力学的铀尾矿库溃坝三维数值模拟方法 |
CN112697673A (zh) * | 2020-12-14 | 2021-04-23 | 中国水利水电科学研究院 | 一种穿堤无压涵管接触渗流破坏的可视化试验装置与方法 |
CN112697673B (zh) * | 2020-12-14 | 2021-09-17 | 中国水利水电科学研究院 | 一种穿堤无压涵管接触渗流破坏的可视化试验装置与方法 |
CN113111489A (zh) * | 2021-03-11 | 2021-07-13 | 重庆交通大学 | 一种堰塞坝洪水漫顶溃坝过程仿真方法及仿真系统 |
CN113297753A (zh) * | 2021-06-24 | 2021-08-24 | 四川大学 | 山区河流透过性坝沿程水深预测方法 |
CN115455867A (zh) * | 2022-10-31 | 2022-12-09 | 武汉大学 | 基于回归分析的坝区流态的推求方法 |
CN117034801A (zh) * | 2023-08-07 | 2023-11-10 | 中国水利水电科学研究院 | 一种基于河道与防洪保护区一体化的溃堤洪水模拟方法 |
CN117034801B (zh) * | 2023-08-07 | 2024-01-02 | 中国水利水电科学研究院 | 一种基于河道与防洪保护区一体化的溃堤洪水模拟方法 |
CN116757125A (zh) * | 2023-08-18 | 2023-09-15 | 水利部交通运输部国家能源局南京水利科学研究院 | 基于sph-dem算法的土石坝漫顶溃决模拟方法 |
CN117521404A (zh) * | 2023-11-30 | 2024-02-06 | 中国水利水电科学研究院 | 一种平原水库溃坝洪水危险等级划分方法 |
CN117494478A (zh) * | 2023-12-29 | 2024-02-02 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种模拟心墙坝漫顶溃坝流量及溃口演化过程的计算方法 |
CN117494478B (zh) * | 2023-12-29 | 2024-03-08 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种模拟心墙坝漫顶溃坝流量及溃口演化过程的计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110362925B (zh) | 2020-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110362925A (zh) | 一种包含库区的土石坝漫顶溃决洪水数值模拟方法 | |
CN108256193B (zh) | 一种河道溃决洪水一、二维数学模型耦合分析方法 | |
CN110046469B (zh) | 多约束条件下水电站坝前河床冲淤变形的计算方法 | |
Fang et al. | Flood management selections for the Yangtze River midstream after the Three Gorges Project operation | |
Wang et al. | Long-term (1842–2006) morphological change and equilibrium state of the Changjiang (Yangtze) Estuary, China | |
CN109918821A (zh) | 一种迎风守恒型河道漫溢出流数值模拟方法 | |
CN108021780B (zh) | 一种基于无规则非结构网格模型的山洪动态仿真方法 | |
CN112784502B (zh) | 一种水文水力学动态双向耦合的洪水预测方法 | |
CN112464584B (zh) | 自由表面流的水位和流量推求方法 | |
CN104598701B (zh) | 一种三维场景分段迭代的河网污染仿真推演方法 | |
CN106407530A (zh) | 一种梯级水库泥沙冲淤同步联合计算方法 | |
CN104933268A (zh) | 一种基于一维非恒定流数值模型的洪水分析方法 | |
CN104750972A (zh) | 复杂水网地区楔形水体容积的计算方法 | |
Isaac et al. | Sediment management studies of a run-of-the-river hydroelectric project using numerical and physical model simulations | |
CN110147646B (zh) | 一种数值模拟框架下线性挡水构筑物的过流处理方法 | |
Pisaturo et al. | A procedure for human safety assessment during hydropeaking events | |
Ishikawa et al. | Assessment of flood risk management in lowland Tokyo areas in the seventeenth century by numerical flow simulations | |
CN103774605A (zh) | 一种提高环抱式港池水体交换能力的设计方法 | |
Burkham | Effects of changes in an alluvial channel on the timing, magnitude, and transformation of flood waves, southeastern Arizona | |
CN115293037A (zh) | 基于机器学习的河湖复合系统水动力建模方法 | |
Sharma | Mathematical modelling and braid indicators | |
Altinakar et al. | New generation modeling and decision support tools for studying impacts of dam failures | |
Ali | Assessment of vanes effectiveness in controlling erosion and deposition zones at a river confluence using a 2D model | |
Le et al. | F28: A novel coupling strategy for 1D/2D hydraulic models for flood risk assessment of the Mekong Delta | |
CN114459571B (zh) | 一种基于初态水位曲线簇河道初始水位确定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |