CN107016185A - 一种冰碛湖溃决洪水峰值流量预测的计算方法 - Google Patents
一种冰碛湖溃决洪水峰值流量预测的计算方法 Download PDFInfo
- Publication number
- CN107016185A CN107016185A CN201710207771.1A CN201710207771A CN107016185A CN 107016185 A CN107016185 A CN 107016185A CN 201710207771 A CN201710207771 A CN 201710207771A CN 107016185 A CN107016185 A CN 107016185A
- Authority
- CN
- China
- Prior art keywords
- dam
- formula
- lake
- drift
- peak flow
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Computer Hardware Design (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Separation Of Suspended Particles By Flocculating Agents (AREA)
Abstract
本发明公开了一种冰碛湖溃决洪水峰值流量预测的计算方法,包括以下步骤:测量初始涌浪高度、冰碛湖水深、冰碛湖的长度、冰碛坝坝体鞍部沿水流方向长度和冰碛坝坝体组成物质的分选系数;将各因子进行组合并消除量纲,与峰值流量进行多元回归方程分析,确定各因子式的影响系数;得到冰碛湖溃决洪水峰值流量Qm的预测计算公式:式中,Qm为冰碛湖溃决洪水峰值流量,m3/s;A0为初始涌浪高度,m;Ld为冰碛坝坝体鞍部沿水流方向长度,m;H0为冰碛湖水深,m;L为冰碛湖长度,m;Cu为分选系数,表征冰碛坝坝体组成的均匀程度,无量纲。本发明极大提高了冰碛湖溃决洪水峰值流量预测计算依据的科学性和准确性。
Description
技术领域
本发明属于水利工程领域,具体地说,涉及一种冰碛湖溃决洪水峰值流量预测的计算方法。
背景技术
近半个世纪以来,随着全球气候变化,世界各高山地区冰碛湖溃决事件频繁发生,造成严重的经济损失和人员伤亡,如喜马拉雅山、天山、阿尔卑斯山脉、高加索山脉、科迪勒拉山系和斯堪的纳维亚山脉等地区。中国是冰湖分布较多的国家之一,冰碛湖溃决主要发生在西藏的喜马拉雅山脉,所形成的洪水和泥石流在该地区造成严重的危害。同时,前人的科学研究发现,冰崩入湖造成湖水涌浪过流漫顶诱发冰碛坝的失效破坏,是冰碛湖溃决的重要成因机制。在我国西藏地区,历史上先后发生的12次冰碛湖溃决事件,涌浪是导致冰湖溃决最直接的诱发因素。
目前,国内外可以直接应用于冰碛湖溃决流量计算的公式较少,且均未考虑涌浪条件的特征参数。现有洪峰流量计算一般借鉴简单的回归公式或已有的土石坝溃决模型,计算方法主要通过已发生溃决的各种坝体(混凝土和土石水工坝,滑坡坝,泥石流坝,冰川侧碛、冰碛坝等)的坝高和库容得出溃决峰值流量的回归公式。如Kirkpatrick(Kirkpatrick,GA.1977.Evaluation guidelines for spillway adequacy.In theEvaluation of dam safety,Proc.Eng.Found.Conf.,New York,Am.Soc.Civ.Eng.,pp:395-414.)根据21个历史上和假设的溃决事件中的坝体高度,得出溃决峰值流量的计算公式:Qm=2.297(H+1)2.5,H为坝高,美国水土保持学会(刘宁,程尊兰,崔鹏和陈宁生等.2013.堰塞湖及其风险控制.北京:科学出版社.pp:141-142.)通过另外的13次溃决将上述公式修订为:Qm=65H1.85,又通过31次坝高在1.8-84m范围的溃决事件将公式修订为:Qm=48H1.63。实践证明,现有此类峰值流量计算方法很难应用于冰碛湖溃决的案例预测中,归咎于不同的学者依据不同的资料得出不同经验公式,且计算误差较大,适用性较低。
发明内容
有鉴于此,本发明针对涌浪条件对冰碛湖溃决形成机制的影响,提高冰碛湖溃决洪水峰值流量预测计算依据的科学性和准确性,提供了一种冰碛湖溃决洪水峰值流量预测的计算方法,该计算方法根据野外案例调查和室内模型试验对冰碛湖溃决形成机制影响的研究结果,计算参数中包括初始涌浪高度、冰碛湖水深、冰碛湖的长度、冰碛坝坝体鞍部沿水流方向长度和冰碛坝坝体组成物质的分选系数5个因子,提高冰碛湖溃决洪水峰值流量预测的科学性和准确性。
为了解决上述技术问题,本发明公开了一种冰碛湖溃决洪水峰值流量预测的计算方法,包括以下步骤:
步骤1、测量初始涌浪高度A0、冰碛湖水深H0、冰碛湖的长度L、冰碛坝坝体鞍部沿水流方向长度Ld和冰碛坝坝体组成物质的分选系数Cu;
步骤2、将各因子进行组合并消除量纲,与峰值流量进行多元回归方程分析,确定各因子式的影响系数;
步骤3、根据各因子式的影响系数,得到冰碛湖溃决洪水峰值流量Qm的预测计算公式:
式中,Qm为冰碛湖溃决洪水峰值流量,m3/s;A0为初始涌浪高度,m;Ld为冰碛坝坝体鞍部沿水流方向长度,m;H0为冰碛湖水深,m;L为冰碛湖长度,m;Cu为分选系数,表征冰碛坝坝体组成的均匀程度,无量纲。
进一步地,步骤2中的将各因子进行组合并消除量纲,与峰值流量进行多元回归方程分析,确定各因子式的影响系数具体为:
步骤2.1、建立峰值流量与诸因素之间的基础关系式,根据基础关系式建立初始峰值流量的预测公式;
步骤2.2、利用多元线性回归分析对初始峰值流量的预测公式进行求解。
进一步地,步骤2.1中的建立峰值流量与诸因素之间的基础关系式具体为:
步骤2.1.1、此式表示以初始涌浪波幅为表征的涌浪条件与坝体鞍部沿水流方向长度的对比关系;组合式分母采用坝体鞍部沿水流方向长度Ld,单位与涌浪波幅一致以消去量纲;
步骤2.1.2、此式表示作为湖水冲刷能力主要影响因素的湖区水深H0与作为坝体抗冲能力主要影响因素的鞍部横截面积E之间的比例关系;分母采用横截面积E的平方根形成,从而分子分母单位一致以消去量纲;
步骤2.1.3、此式表示作为湖水冲刷能力主要影响因素的湖区长度L与作为坝体抗冲能力主要影响因素的鞍部横截面积E之间的比例关系;分母采用横截面积E的平方根形成,从而分子分母单位一致以消去量纲;
步骤2.1.4、Cu:此因子表征坝体组成的均匀程度,拣选系数Cu是过筛质量占75%的粒径D75与过筛质量占25%的粒径D25比值的平方根,其表达式为Cu=(D75/D25)1/2,为无量纲参数;
步骤2.1.5、将上述4个基础关系式,按照其之间的相关性,列出构成涌浪溢流溃坝峰值流量预测公式的一般表达式(1):
将公式(1)进行转换得到公式(2):
其中,Qm为冰碛湖溃决洪水峰值流量,m3/s;A0为初始涌浪高度,m;Ld为冰碛坝坝体鞍部沿水流方向长度,m;H0为冰碛湖水深,m;L为冰碛湖长度,m;Cu为分选系数,表征冰碛坝坝体组成的均匀程度,无量纲;a,α,β,γ和δ为经验公式未知系数,无量纲;
将公式(2)中各组合式取对数构成直线方程,得式变换式(3)如下:
设Y=log Qm,b=log a,X4=logCu,代入式(3)得到多元直线方程转换式(4)如下:
Y=b+αX1+βX2+γX2+δX4 (4)。
进一步地,步骤2.2中利用多元线性回归分析对建立峰值流量的预测公式进行求解计算得出b=2.766,α=0.244,β=1.307,γ=0.538和δ=0.592。
进一步地,步骤3中根据各因子式的影响系数,得到冰碛湖溃决洪水峰值流量Qm的预测计算公式具体为:根据b=2.766,α=0.244,β=1.307,γ=0.538和δ=0.592,从而得到方程的表达式(5)如下:
Y=2.766+0.244X1+1.307X2+0.538X3+0.592X4 (5)
将Y=log Qm,b=log a,X4=logCu代入公式(5)得到转换式(6)如下:
将公式(6)等式两边的各组合式去对数,得到方程转换式(7)如下:
将公式(7)等式两边参数因子合并简化,得到方程转换式(8)如下:
公式(8)即为涌浪溢流溃坝的峰值流量预测公式;其中,Qm是无量纲峰值流量,其单位与因子式单位一致;
将Ld=2E/H0代入公式(8),将指数小数点保留2位,得到计算公式(9):
式中,Qm为冰碛湖溃决洪水峰值流量,m3/s;A0为初始涌浪高度,m;Ld为冰碛坝坝体鞍部沿水流方向长度,m;H0为冰碛湖水深,m;L为冰碛湖长度,m;Cu为分选系数,表征冰碛坝坝体组成的均匀程度,无量纲。
与现有技术相比,本发明可以获得包括以下技术效果:
1)本发明通过野外案例调查和室内模型试验对冰碛湖溃决形成机制影响的研究发现,涌浪是冰碛湖溃决形成的直接原因:当初始涌浪高度越大时,溃决发展速率快,产生的峰值流量偏大;当初始涌浪高度越小时,溃决发展速率慢,产生的峰值流量偏小。基于上述研究分析,本发明针对涌浪条件对冰碛湖溃决形成机制的影响,计算参数中包括初始涌浪高度、冰碛湖水深、冰碛湖的长度、冰碛坝坝体鞍部沿水流方向长度和冰碛坝坝体组成物质的分选系数5个因子,提出一种冰碛湖溃决洪水峰值流量预测的计算方法。
2)本发明极大提高了冰碛湖溃决洪水峰值流量预测计算依据的科学性和准确性。
当然,实施本发明的任一产品并不一定需要同时达到以上所述的所有技术效果。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本发明的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明冰碛湖溃决洪水峰值流量预测计算参数因子简化模型图;其中,H0:冰碛湖水深,A0:初始涌浪高度,Cu:冰碛坝组成物质分选系数,L:冰碛湖长度,Ld:冰碛坝坝体鞍部沿水流方向长度;
图2是本发明预测值Qm *与实测值Qm对比。
具体实施方式
以下将配合实施例来详细说明本发明的实施方式,藉此对本发明如何应用技术手段来解决技术问题并达成技术功效的实现过程能充分理解并据以实施。
本发明公开了一种冰碛湖溃决洪水峰值流量预测的计算方法,包括以下步骤:
步骤1、通过室内模型试验结果中各因子对冰碛湖溃决形成机制的研究分析,确定影响溃决洪水峰值流量Qm的5个主要因子,分别是:初始涌浪高度A0、冰碛湖水深H0、冰碛湖的长度L、冰碛坝坝体鞍部沿水流方向长度Ld和冰碛坝坝体组成物质的分选系数Cu;
其中,初始涌浪高度A0:表征冰崩入湖形成涌浪且传播至冰碛坝坝顶首个波浪的高度,根据冰碛湖湖盆边缘波痕实地调查获取。
冰碛坝坝体鞍部沿水流方向长度Ld:根据冰碛坝实地调查测量获取;Ld=2E/H0,。
冰碛湖水深H0:根据冰碛坝实地测量获取得到的坝体高度,冰碛湖水深取值冰碛坝坝高值。
冰碛湖长度L:根据冰碛湖实地调查或地形图测量获取。
冰碛坝组成物质分选系数Cu:表征坝体组成的均匀程度,根据冰碛坝实地调查和冰碛物取样筛分获取,分选系数Cu是过筛质量占75%的粒径D75与过筛质量占25%的粒径D25比值的平方根,其表达式为Cu=(D75/D25)1/2,为无量纲参数。
步骤2、将各因子进行组合并消除量纲,与峰值流量进行多元回归方程分析,确定各因子式的影响系数;
步骤2.1、建立峰值流量与诸因素之间的基础关系式,建立峰值流量的预测公式;
按照各因素对峰值流量的影响机制分析,需要对各因子进行组合,并达成无量纲,以达到经验公式应用到实践的可能。
(1)此式表示以初始涌浪波幅为表征的涌浪条件与坝体鞍部沿水流方向长度的对比关系。根据实验分析,涌浪条件是涌浪溢流溃坝溃口雏形形成的直接因素,亦与溃口侵蚀速率和峰值流量呈正比关系。因此,与常见的普通溢流溃坝峰值流量经验公式不同,考虑到涌浪溢流溃坝形成机制的特殊性,本发明将创新性地在关系式中引入以初始涌浪波幅A0为表征参数的涌浪条件。组合式分母采用坝体鞍部沿水流方向长度Ld(Ld是坝体鞍部面积E的组成参数之一),单位与涌浪波幅一致以消去量纲。同时,根据“坝体特征对溃坝形成机制的影响”章节的实验分析,以沿水流方向坝体长度为主要特征参数的鞍部横截面面积越大,则溃口侵蚀速率和峰值流量相反越小。
(2)此式表示作为湖水冲刷能力主要影响因素的湖区水深H0与作为坝体抗冲能力主要影响因素的鞍部横截面积E之间的比例关系。根据实验分析,湖区水深H0与湖区长度L虽都是影响湖区库容的主要因子,但与溃口侵蚀速率和峰值流量的关系却稍有不同。因此,考虑到这种差异影响机制,本发明中不采用湖区库容W与坝体规模的直接对比关系,而是采用分别表征的方法。为达到量纲一致,分母采用横截面积E的平方根形成。此式比值越大,即湖水冲刷能力相对大,坝体抵抗湖水的抗冲刷能力相对小,从而溃决速度快,峰值流量大。
(3)此式表示作为湖水冲刷能力主要影响因素的湖区长度L与作为坝体抗冲能力主要影响因素的鞍部横截面积E之间的比例关系。如组合式(2)所述,湖区长度L与湖区水深H0均为影响湖区库容的主要因子,考虑两者对溃坝形成机制的差异影响,因此分别进行讨论。类似地,为达到量纲一致,分母采用横截面积E的平方根形成。
(4)Cu:此因子表征坝体组成的均匀程度,拣选系数Cu是过筛质量占75%的粒径D75与过筛质量占25%的粒径D25比值的平方根,其表达式为Cu=(D75/D25)1/2,为无量纲参数。根据实验分析,拣选系数Cu是影响坝体的抗冲刷能力主要因子之一,与溃口侵蚀速率和峰值流量关系密切。其值越大,溃口侵蚀速率越快,峰值流量越大。
将上述四个组合式(或因子),按照其之间的相关性,列出构成涌浪溢流溃坝峰值流量预测公式的一般表达式(1),并把结果整理成幂次关系,经转换得到表达式(2):
其中,Qm为冰碛湖溃决洪水峰值流量,m3/s;A0为初始涌浪高度,m;Ld为冰碛坝坝体鞍部沿水流方向长度,m;H0为冰碛湖水深,m;L为冰碛湖长度,m;Cu为分选系数,表征冰碛坝坝体组成的均匀程度,无量纲;a,α,β,γ和δ为经验公式未知系数,无量纲;
为了使预测公式适用于涌浪溢流溃坝,采用模拟实验中成功溃决的56组涌浪案例的实验参数,通过统计分析,建立无量纲峰值流量Qm的预测公式。
将公式(2)中各组合式取对数构成直线方程,得式变换式(3)如下:
设Y=log Qm,b=log a,X4=logCu,代入式(3)得到多元直线方程转换式(4)如下:
Y=b+αX1+βX2+γX3+δX4 (4)
下一步工作即是将56组实验数据代入公式(4),解得多元回归方程的未知系数b,α,β,γ和δ。
其中,56组试验,试验中改变试验变量中的A0、H0、Ld、L、E、Cu6个参数,从而得到不同的Qm。a,α,β,γ和δ是根据这些数据得到的一个回归系数。
步骤2.2、多元线性回归分析求解
本发明中,多元线性回归方程的建立采用SPSS17.0进行分析求解。将56实验数据分别转化为Y,X1,X2,X3和X4后输入软件,得到定义模型的输出结果如下:
(1)模型概述
列出了模型的R,R2,调整R2及估计标准误。R2越大所反映的自变量与因变量的共变量比率越高,模拟与数据的拟合程度越好,其输出结果如表1所示。
表1模型概述
a预测变量:(常量),X1,X2,X3和X4
研究定义模型确定系数的平方根R为0.784,确定系数R2为0.614,调整后的确定系数R2为0.580,估计标准误为0.069。
(2)方差分析
列出了变异源、自由度、均方,F值及对F值的显著性检验。显著性水平是判别回归方程的是否有效的重要标准,输出结果如表2所示。
表2方差分析(ANOVAb)
a预测变量:(常量),X1,X2,X3和X4
b因变量:Y
研究定义模型的回归平方和为0.341,残差平方和为0.214,总平方和味0.556,F统计量的值为17.9263,Sig.<0.05,可认为所建立的回归方程有效。
(3)回归系数
列出了常数及回归系数的值及标准化的值,同时对其进行显著性检验,输出结果如表3所示。
表3回归系数(Coefficientsa)
a因变量:Y
研究定义模型中因变量Y对四个自变量X1,X2,X3和X4的回归非标准化回归系数分别是0.244,1.307,0.538和0.592;对应显著性检验的t值分别是2.995,3.293,3.440和4.401;对应四个回归系数B的显著性水平Sig.分别是0.004,0.002,0.001和0.000,Sig.<0.05,可以认为自变量X1,X2,X3和X4对因变量Y均有显著影响。
根据上述分析结果,解得公式(4)方程中4个未知系数,即b=2.766,α=0.244,β=1.307,γ=0.538和δ=0.592。
步骤3、建立预测公式
根据b=2.766,α=0.244,β=1.307,γ=0.538和δ=0.592,从而得到方程的表达式(5)如下:
Y=2.766+0.244X1+1.307X2+0.538X3+0.592X4 (5)
将Y=log Qm,b=log a,X4=logCu代入公式(5)得到转换式(6)如下:
将公式(6)等式两边的各组合式去对数,得到方程转换式(7)如下:
将公式(7)等式两边参数因子合并简化,得到方程转换式(8)如下:
公式(8)即为涌浪溢流溃坝的峰值流量预测公式。其中,Qm是无量纲峰值流量,其单位与因子式单位一致。将56组实验参数代入预测公式进行计算,得到涌浪溢流溃坝产生的峰值流量预测值,图2是实测值与预测值对比结果。
将Ld=2E/H0代入公式(8),然后指数小数点保留2位简化即可,得到本发明的计算公式(9):
式中,Qm为冰碛湖溃决洪水峰值流量,m3/s;A0为初始涌浪高度,m;Ld为冰碛坝坝体鞍部沿水流方向长度,m;H0为冰碛湖水深,m;L为冰碛湖长度,m;Cu为分选系数,表征冰碛坝坝体组成的均匀程度,无量纲。
实施例1
如图1所示,藏东南某冰碛湖溃决前水深H0=15m,湖区长度为L=950m;溃决前冰碛坝坝体鞍部沿水流方向长度L=100m。1988年7月,冰碛湖发生溃决,形成规模巨大的泥石流,席卷沟内村落及农田,堵塞主河,冲毁川藏公路,波及下游县城,造成断道堵车一年之久的危害。
采用本发明提出的计算方法进行事件还原和验证。除已知冰碛湖和冰碛坝的尺寸参数外,通过实地测量和牧民调访,到达冰碛坝的初始涌浪高度A0=3m;根据现场冰碛坝取样和筛分,得到冰碛坝坝体组成的分选系数Cu=3.53。将5个计算参数因子代入发明提出的计算方法,得到此次冰碛湖溃决洪水的峰值流量Qm=1483m3/s,相对绝大部分经验公式和半经验理论公式的值偏小,虽略高但接近的实测值(1270m3/s),误差为+17%。
为了检验预测公式的合理性,将溃决前的冰湖和冰碛坝特征代入已有的经验公式进行计算对比,计算结果如表4所示。
表4光谢错溃决峰值流量经验公式计算结果
通过表1的计算结果对比可以看出,本发明涉及的预测方法所计算得到的峰值流量值(1529m3/s)相对绝大部分经验公式和半经验理论公式的值偏小,虽略高但接近的实测值(1270m3/s),误差仅为+17%。
实施例2
如图1所示。藏东南某冰碛湖呈长方形,湖区长度1.94km,湖区最大宽度0.63km,平均宽度0.50km。湖面面积0.976km2,平均湖区水深36m左右,据估算湖水体积达到35.71×106m3。一旦冰碛湖发生溃决,将会给川藏线和沿江城镇、村庄的生命财产造成极其严重的灾难。
根据室内地形图测量、现场调查测量和冰碛坝取样分析,得到5个计算因子分别为:冰碛湖水深H0=36m,湖区长度为L=1940m,冰碛坝坝体鞍部沿水流方向长度L=182m,到达冰碛坝的初始涌浪高度A0=5m,冰碛坝坝体组成的分选系数Cu=3.53。将5个计算参数因子代入发明提出的计算方法,预测得到该冰碛湖溃决洪水的峰值流量Qm=1806m3/s。
上述说明示出并描述了发明的若干优选实施例,但如前所述,应当理解发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述发明构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离发明的精神和范围,则都应在发明所附权利要求的保护范围内。
Claims (5)
1.一种冰碛湖溃决洪水峰值流量预测的计算方法,其特征在于,包括以下步骤:
步骤1、测量初始涌浪高度A0、冰碛湖水深H0、冰碛湖的长度L、冰碛坝坝体鞍部沿水流方向长度Ld和冰碛坝坝体组成物质的分选系数Cu;
步骤2、将各因子进行组合并消除量纲,与峰值流量进行多元回归方程分析,确定各因子式的影响系数;
步骤3、根据各因子式的影响系数,得到冰碛湖溃决洪水峰值流量Qm的预测计算公式:
式中,Qm为冰碛湖溃决洪水峰值流量,m3/s;A0为初始涌浪高度,m;Ld为冰碛坝坝体鞍部沿水流方向长度,m;H0为冰碛湖水深,m;L为冰碛湖长度,m;Cu为分选系数,表征冰碛坝坝体组成的均匀程度,无量纲。
2.根据权利要求1所述的冰碛湖溃决洪水峰值流量预测的计算方法,其特征在于,步骤2中的将各因子进行组合并消除量纲,与峰值流量进行多元回归方程分析,确定各因子式的影响系数具体为:
步骤2.1、建立峰值流量与诸因素之间的基础关系式,根据基础关系式建立初始峰值流量的预测公式;
步骤2.2、利用多元线性回归分析对初始峰值流量的预测公式进行求解。
3.根据权利要求2所述的冰碛湖溃决洪水峰值流量预测的计算方法,其特征在于,所述步骤2.1中的建立峰值流量与诸因素之间的基础关系式具体为:
步骤2.1.1、此式表示以初始涌浪波幅为表征的涌浪条件与坝体鞍部沿水流方向长度的对比关系;组合式分母采用坝体鞍部沿水流方向长度Ld,单位与涌浪波幅一致以消去量纲;
步骤2.1.2、此式表示作为湖水冲刷能力主要影响因素的湖区水深H0与作为坝体抗冲能力主要影响因素的鞍部横截面积E之间的比例关系;分母采用横截面积E的平方根形成,从而分子分母单位一致以消去量纲;
步骤2.1.3、此式表示作为湖水冲刷能力主要影响因素的湖区长度L与作为坝体抗冲能力主要影响因素的鞍部横截面积E之间的比例关系;分母采用横截面积E的平方根形成,从而分子分母单位一致以消去量纲;
步骤2.1.4、Cu:此因子表征坝体组成的均匀程度,拣选系数Cu是过筛质量占75%的粒径D75与过筛质量占25%的粒径D25比值的平方根,其表达式为Cu=(D75/D25)1/2,为无量纲参数;
步骤2.1.5、将上述4个基础关系式,按照其之间的相关性,列出构成涌浪溢流溃坝峰值流量预测公式的一般表达式(1):
将公式(1)进行转换得到公式(2):
其中,Qm为冰碛湖溃决洪水峰值流量,m3/s;A0为初始涌浪高度,m;Ld为冰碛坝坝体鞍部沿水流方向长度,m;H0为冰碛湖水深,m;L为冰碛湖长度,m;Cu为分选系数,表征冰碛坝坝体组成的均匀程度,无量纲;a,α,β,γ和δ为经验公式未知系数,无量纲;
将公式(2)中各组合式取对数构成直线方程,得式变换式(3)
如下:
设Y=log Qm,b=log α,X4=log Cu,代入式(3)得到多元直线方程转换式(4)如下:
Y=b+αX1+βX2+γX3+δX4 (4)。
4.根据权利要求3所述的冰碛湖溃决洪水峰值流量预测的计算方法,其特征在于,所述步骤2.2中利用多元线性回归分析对建立峰值流量的预测公式进行求解计算得出b=2.766,α=0.244,β=1.307,γ=0.538和δ=0.592。
5.根据权利要求3所述的冰碛湖溃决洪水峰值流量预测的计算方法,其特征在于,步骤3中根据各因子式的影响系数,得到冰碛湖溃决洪水峰值流量Qm的预测计算公式具体为:根据b=2.766,α=0.244,β=1.307,γ=0.538和δ=0.592,从而得到方程的表达式(5)如下:
Y=2.766+0.244X1+1.307X2+0.538X3+0.592X4 (5)
将Y=log Qm,b=log a,X4=log Cu代入公式(5)得到转换式(6)如下:
将公式(6)等式两边的各组合式去对数,得到方程转换式(7)如下:
将公式(7)等式两边参数因子合并简化,得到方程转换式(8)如下:
公式(8)即为涌浪溢流溃坝的峰值流量预测公式;其中,Qm是无量纲峰值流量,其单位与因子式单位一致;
将Ld=2E/H0代入公式(8),将指数小数点保留2位,得到计算公式(9):
式中,Qm为冰碛湖溃决洪水峰值流量,m3/s;A0为初始涌浪高度,m;Ld为冰碛坝坝体鞍部沿水流方向长度,m;H0为冰碛湖水深,m;L为冰碛湖长度,m;Cu为分选系数,表征冰碛坝坝体组成的均匀程度,无量纲。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710207771.1A CN107016185B (zh) | 2017-03-31 | 2017-03-31 | 一种冰碛湖溃决洪水峰值流量预测的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710207771.1A CN107016185B (zh) | 2017-03-31 | 2017-03-31 | 一种冰碛湖溃决洪水峰值流量预测的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107016185A true CN107016185A (zh) | 2017-08-04 |
CN107016185B CN107016185B (zh) | 2020-04-14 |
Family
ID=59445763
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710207771.1A Active CN107016185B (zh) | 2017-03-31 | 2017-03-31 | 一种冰碛湖溃决洪水峰值流量预测的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107016185B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107564245A (zh) * | 2017-09-28 | 2018-01-09 | 成都理工大学 | 一种考虑降雨的冰湖溃决的预报方法 |
CN108038273A (zh) * | 2017-11-24 | 2018-05-15 | 西安理工大学 | 一种瞬间溃坝坝址流量过程线的快速预测方法 |
CN108681651A (zh) * | 2018-05-29 | 2018-10-19 | 成都理工大学 | 一种无冰芯的冰碛堤冰湖溃决的早期识别方法 |
CN108733948A (zh) * | 2018-05-29 | 2018-11-02 | 成都理工大学 | 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 |
CN108959681A (zh) * | 2018-04-03 | 2018-12-07 | 四川大学 | 一种基于点云数据的堰塞湖引流道设计方法 |
CN110210173A (zh) * | 2019-06-20 | 2019-09-06 | 中国地质科学院探矿工艺研究所 | 一种冰碛土物理力学参数的快速估算方法 |
CN110263453A (zh) * | 2019-06-25 | 2019-09-20 | 长江水利委员会水文局 | 一种自适应的堰塞湖逐渐溃耦合预测模型 |
CN110472183A (zh) * | 2019-08-21 | 2019-11-19 | 中国地质科学院探矿工艺研究所 | 一种基于灾害时空演化的泥石流容重计算方法 |
CN110779499A (zh) * | 2019-10-25 | 2020-02-11 | 成都理工大学 | 一种考虑冰湖特征的冰崩形成涌浪高度的计算方法 |
CN110836663A (zh) * | 2019-10-25 | 2020-02-25 | 成都理工大学 | 一种冰湖涌浪高度随传播距离的计算方法及应用 |
CN110866330A (zh) * | 2019-10-28 | 2020-03-06 | 成都理工大学 | 一种考虑冰湖前地形特征的冰崩产生涌浪的计算方法 |
CN110909293A (zh) * | 2019-10-25 | 2020-03-24 | 成都理工大学 | 一种冰湖涌浪高度的计算方法及其应用 |
CN113405769A (zh) * | 2021-06-18 | 2021-09-17 | 中国科学院、水利部成都山地灾害与环境研究所 | 模拟冰碛湖漫顶溃决灾害全过程的实验装置及实验方法 |
CN113420453A (zh) * | 2021-06-30 | 2021-09-21 | 中国科学院、水利部成都山地灾害与环境研究所 | 泥石流堰塞坝溃决洪峰流量计算方法 |
CN113552628A (zh) * | 2021-07-20 | 2021-10-26 | 黄河勘测规划设计研究院有限公司 | 地震涌浪高度计算方法 |
CN113779491A (zh) * | 2021-09-08 | 2021-12-10 | 北京市地质研究所 | 泥石流隐患易发性分析方法 |
CN115637673A (zh) * | 2022-12-23 | 2023-01-24 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种冰崩型冰湖溃决涌浪生成及高度测量装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102086635A (zh) * | 2010-12-31 | 2011-06-08 | 中国科学院水利部成都山地灾害与环境研究所 | 一种主河输移控制型泥石流防治方法 |
CN103306242A (zh) * | 2013-07-02 | 2013-09-18 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种调节洪峰流量的泥石流控流坝及其设计方法和应用 |
JP2014196656A (ja) * | 2013-03-05 | 2014-10-16 | 株式会社山辰組 | サイフォン送水装置及びサイフォン送水方法 |
CN104217424A (zh) * | 2014-07-08 | 2014-12-17 | 西藏自治区水利电力规划勘测设计研究院 | 冰湖终碛垄颗粒的分析方法 |
CN105862653A (zh) * | 2016-04-08 | 2016-08-17 | 华电电力科学研究院 | 一种基于数理统计模型的水库防洪效益动态计算方法 |
-
2017
- 2017-03-31 CN CN201710207771.1A patent/CN107016185B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102086635A (zh) * | 2010-12-31 | 2011-06-08 | 中国科学院水利部成都山地灾害与环境研究所 | 一种主河输移控制型泥石流防治方法 |
JP2014196656A (ja) * | 2013-03-05 | 2014-10-16 | 株式会社山辰組 | サイフォン送水装置及びサイフォン送水方法 |
CN103306242A (zh) * | 2013-07-02 | 2013-09-18 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种调节洪峰流量的泥石流控流坝及其设计方法和应用 |
CN104217424A (zh) * | 2014-07-08 | 2014-12-17 | 西藏自治区水利电力规划勘测设计研究院 | 冰湖终碛垄颗粒的分析方法 |
CN105862653A (zh) * | 2016-04-08 | 2016-08-17 | 华电电力科学研究院 | 一种基于数理统计模型的水库防洪效益动态计算方法 |
Non-Patent Citations (4)
Title |
---|
刘建康 等: ""云南鲁甸红石岩堰塞湖溃坝风险及其影响"", 《山地学报》 * |
刘建康 等: ""藏东南典型冰湖溃决危险性分析"", 《灾害学》 * |
王世金 等: ""冰湖溃决灾害风险研究进展及其展望"", 《水科学进展》 * |
黄金辉 等: ""涌浪规模对冰碛湖溃决的影响实验"", 《山地学报》 * |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107564245B (zh) * | 2017-09-28 | 2019-07-12 | 成都理工大学 | 一种考虑降雨的冰湖溃决的预报方法 |
CN107564245A (zh) * | 2017-09-28 | 2018-01-09 | 成都理工大学 | 一种考虑降雨的冰湖溃决的预报方法 |
CN108038273A (zh) * | 2017-11-24 | 2018-05-15 | 西安理工大学 | 一种瞬间溃坝坝址流量过程线的快速预测方法 |
CN108959681B (zh) * | 2018-04-03 | 2020-08-04 | 四川大学 | 一种基于点云数据的堰塞湖引流道设计方法 |
CN108959681A (zh) * | 2018-04-03 | 2018-12-07 | 四川大学 | 一种基于点云数据的堰塞湖引流道设计方法 |
CN108681651A (zh) * | 2018-05-29 | 2018-10-19 | 成都理工大学 | 一种无冰芯的冰碛堤冰湖溃决的早期识别方法 |
CN108733948A (zh) * | 2018-05-29 | 2018-11-02 | 成都理工大学 | 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 |
CN108733948B (zh) * | 2018-05-29 | 2021-11-19 | 成都理工大学 | 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 |
CN108681651B (zh) * | 2018-05-29 | 2021-10-29 | 成都理工大学 | 一种无冰芯的冰碛堤冰湖溃决的早期识别方法 |
CN110210173A (zh) * | 2019-06-20 | 2019-09-06 | 中国地质科学院探矿工艺研究所 | 一种冰碛土物理力学参数的快速估算方法 |
CN110210173B (zh) * | 2019-06-20 | 2021-09-10 | 中国地质科学院探矿工艺研究所 | 一种冰碛土物理力学参数的快速估算方法 |
CN110263453A (zh) * | 2019-06-25 | 2019-09-20 | 长江水利委员会水文局 | 一种自适应的堰塞湖逐渐溃耦合预测模型 |
CN110472183B (zh) * | 2019-08-21 | 2022-04-05 | 中国地质科学院探矿工艺研究所 | 一种基于灾害时空演化的泥石流容重计算方法 |
CN110472183A (zh) * | 2019-08-21 | 2019-11-19 | 中国地质科学院探矿工艺研究所 | 一种基于灾害时空演化的泥石流容重计算方法 |
CN110836663A (zh) * | 2019-10-25 | 2020-02-25 | 成都理工大学 | 一种冰湖涌浪高度随传播距离的计算方法及应用 |
CN110909293A (zh) * | 2019-10-25 | 2020-03-24 | 成都理工大学 | 一种冰湖涌浪高度的计算方法及其应用 |
CN110779499A (zh) * | 2019-10-25 | 2020-02-11 | 成都理工大学 | 一种考虑冰湖特征的冰崩形成涌浪高度的计算方法 |
CN110779499B (zh) * | 2019-10-25 | 2021-10-29 | 成都理工大学 | 一种考虑冰湖特征的冰崩形成涌浪高度的计算方法 |
CN110909293B (zh) * | 2019-10-25 | 2022-04-26 | 成都理工大学 | 一种冰湖涌浪高度的计算方法 |
CN110866330A (zh) * | 2019-10-28 | 2020-03-06 | 成都理工大学 | 一种考虑冰湖前地形特征的冰崩产生涌浪的计算方法 |
CN110866330B (zh) * | 2019-10-28 | 2023-03-31 | 成都理工大学 | 一种考虑冰湖前地形特征的冰崩产生涌浪的计算方法 |
CN113405769A (zh) * | 2021-06-18 | 2021-09-17 | 中国科学院、水利部成都山地灾害与环境研究所 | 模拟冰碛湖漫顶溃决灾害全过程的实验装置及实验方法 |
CN113420453A (zh) * | 2021-06-30 | 2021-09-21 | 中国科学院、水利部成都山地灾害与环境研究所 | 泥石流堰塞坝溃决洪峰流量计算方法 |
CN113420453B (zh) * | 2021-06-30 | 2023-04-25 | 中国科学院、水利部成都山地灾害与环境研究所 | 泥石流堰塞坝溃决洪峰流量计算方法 |
CN113552628A (zh) * | 2021-07-20 | 2021-10-26 | 黄河勘测规划设计研究院有限公司 | 地震涌浪高度计算方法 |
CN113552628B (zh) * | 2021-07-20 | 2023-08-15 | 黄河勘测规划设计研究院有限公司 | 地震涌浪高度计算方法 |
CN113779491A (zh) * | 2021-09-08 | 2021-12-10 | 北京市地质研究所 | 泥石流隐患易发性分析方法 |
CN113779491B (zh) * | 2021-09-08 | 2023-07-21 | 北京市地质研究所 | 泥石流隐患易发性分析方法 |
CN115637673A (zh) * | 2022-12-23 | 2023-01-24 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种冰崩型冰湖溃决涌浪生成及高度测量装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107016185B (zh) | 2020-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107016185A (zh) | 一种冰碛湖溃决洪水峰值流量预测的计算方法 | |
Bosa et al. | Shallow water numerical model of the wave generated by the Vajont landslide | |
Aksoy et al. | Hydrological and hydraulic models for determination of flood-prone and flood inundation areas | |
CN106971085A (zh) | 一种泥石流沟道侵蚀量的计算方法 | |
Sharifan et al. | Uncertainty and sensitivity analysis of SWMM model in computation of manhole water depth and subcatchment peak flood | |
Herrera-Granados et al. | Numerical and physical modeling of water flow over the ogee weir of the new Niedów barrage | |
Liang et al. | A meso-scale model toward concrete water permeability regarding aggregate permeability | |
Mohammed-Ali et al. | Assessing the Feasibility of an Explicit Numerical Model for Simulating Water Surface Profiles over Weirs. | |
Ji et al. | Progressive geomorphic evolution of reservoir bank in coarse-grained soil in East China–Insights from long-term observations and physical model test | |
Orozco et al. | Comparison between HEC-RAS and TELEMAC-2D hydrodynamic models of the Loire River, integrating levee breaches | |
Zhao et al. | Experimental study of the debris flow slurry impact and distribution | |
Tahershamsi et al. | Reservoir's geometry impact of three dimensions on peak discharge of dam-failure flash flood | |
Andersen et al. | Overtopping on breakwaters with a permeable crest | |
JP2009123026A (ja) | 流量変動予測プログラム | |
Siregar | Numerical Analysis of Flood modeling of upper Citarum River under Extreme Flood Condition | |
Oveici et al. | Flow Field in a Sloped Channel with Damaged and Undamaged Piers: Numerical and Experimental Studies | |
CN117744535B (zh) | 一种松散岩土体堰塞坝泥石流流量峰值预测方法及系统 | |
Reistad et al. | Setting up rating curves using HEC-RAS | |
CN113569444B (zh) | 一种堤身材料渗透系数随机序贯反演方法 | |
Wang et al. | Research on Prediction Method of Foundation Settlement | |
Li et al. | Study on the Wind Erosion Resistance Ability and Slope Stability of a Wind-Eroded Desert Roadbed | |
Chaves et al. | Envoltórias regionais de cheias milenares e decamilenares para o estado do Ceará: Ferramenta para verificação da segurança hidrológica de barragens | |
Mohammed et al. | Seepage through homogenous and non-homogenous earth dams: Comparison between observation and simulation | |
Estación et al. | Predictive numerical modeling of the behavior of asphalt-core rockfill dam | |
Ye et al. | Study on the effect of dry-wet cycle on the dam slope stability |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221110 Address after: No.1, Section 2, North 1st ring road, Chengdu, Sichuan 610000 Patentee after: Chengdu Huajian Geological Engineering Technology Co.,Ltd. Address before: 610081 No.1, North 2nd section of 1st ring road, Jinniu District, Chengdu City, Sichuan Province Patentee before: INSTITUTE OF EXPLORATION TECHNOLOGY OF CHINESE ACADEMY OF GEOLOGICAL SCIENCES |