CN110909293A - 一种冰湖涌浪高度的计算方法及其应用 - Google Patents
一种冰湖涌浪高度的计算方法及其应用 Download PDFInfo
- Publication number
- CN110909293A CN110909293A CN201911023405.6A CN201911023405A CN110909293A CN 110909293 A CN110909293 A CN 110909293A CN 201911023405 A CN201911023405 A CN 201911023405A CN 110909293 A CN110909293 A CN 110909293A
- Authority
- CN
- China
- Prior art keywords
- ice
- lake
- collapse
- surge
- ice lake
- 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
- 238000000034 method Methods 0.000 title claims abstract description 41
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 57
- 238000013507 mapping Methods 0.000 claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- 238000011144 upstream manufacturing Methods 0.000 claims description 22
- 230000005484 gravity Effects 0.000 claims description 19
- 230000000694 effects Effects 0.000 claims description 12
- 238000005381 potential energy Methods 0.000 claims description 11
- 239000013078 crystal Substances 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 230000002265 prevention Effects 0.000 abstract description 17
- 238000002474 experimental method Methods 0.000 description 4
- 238000011835 investigation Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 241000894006 Bacteria Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Tourism & Hospitality (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种冰湖涌浪高度的计算方法,属于灾害防治工程领域,其特征在于,包括以下步骤:a、通过冰崩体的位置,区别冰川的属性;通过调查测绘确定冰崩体的基本数据;b、确定冰崩在冰湖中形成涌浪高度H的参数,包括冰崩体与冰湖湖面的高程差H0、冰崩体与冰湖末端之间的距离J、冰湖末端宽度b、冰湖面积Ag、冰湖平均水深h和冰湖末端前地表坡度P;c、冰崩在冰湖中形成涌浪高度H通过式七计算获得。本发明建立了更精确的冰崩在冰湖中形成的涌浪计算方法,为判断冰湖是否溃决提供依据,极大的提高了主动防灾适用性。
Description
技术领域
本发明涉及到高寒山区灾害防治工程技术领域,尤其涉及一种冰湖涌浪高度的计算方法及其应用。
背景技术
冰崩是一种发生在高寒山区的自然现象,冰崩的下游往往还有冰湖存在。冰崩发生后,冰块运动到冰湖中,往往形成巨大的涌浪,抬高冰湖水位,引发冰湖溃决,形成冰湖溃决泥石流,再冲入下游河道,堵塞河道,再溃决,形成次生灾害链。因此,冰崩不仅可能危害其附近山坡下的范围,还有可能影响到河道上下游。
冰崩的发生往往由于冰川的前部上缘的拉张裂隙,在持续的张力作用下,由于气温的突变,使得冰川的前部完全脱离冰川体,形成冰崩。降雨也会加剧这类冰崩的发生。
滑坡、崩塌、泥石流、雪崩或冰川运动入水冲击海洋、海湾、河湖及水库区水体时都能产生涌浪波。该类型波是非周期性长波,并且有强烈的非线性,介于中等水波至浅水波之间。目前,国内外学者对冰崩在冰湖中形成的涌浪高度的研究非常少,主要是研究滑坡入水造成的涌浪高度计算。黄波林等(黄波林,陈小婷,殷跃平,张衡,滑坡崩塌涌浪计算方法研究,工程地质学报,2012,20(6):909-915)总结了滑坡涌浪高度及传播距离的计算方法,并给出了危岩体入水的涌浪高度计算方法。Zweifel(Zweifel A,Zuccala D,GattiD.Comparison between computed and experimentally generated impulsewaves.Journal of Hydraulic Engin eering,2007,133(2):208-216.)研究了不同滑坡体入水滑体的密度变化对涌浪高度的影响,给出了滑坡入水的涌浪高度计算方法;Panizzo等(Panizzo A.,Girolamo P.De,Petaccia A.Forecasting impulse waves generated bysubaerial landslides.Journal of Geophysical Research.2005,110(C12025),doi:10.1029/2004JC002778);Evers和Hager(Evers F.M.,Hager W.H.Spatial impulsewaves:wave height decay experiments at laboratory scale.Landslides,2016,13:1395-1403.)都计算了涌浪随传播距离衰减的涌浪高度。这些研究对于冰崩进入冰湖形成的涌浪计算主要存在的问题有:
1、滑坡体与冰崩体在质量上有很大的不同:滑坡体容重大,直接进入水体并下沉;冰崩体容重小,小于水的比重;冰崩体进入冰湖后,先下沉再上浮,因此形成的涌浪高度差别很大。
2、涌浪高度计算中没有考虑滑入角度的影响。
3、冰崩体一般体积较大,不仅仅体现在长度长,厚度厚,其宽度也较宽,尤其是对于长条形,宽度比长度小很多的冰湖,冰崩体的宽度的影响不可忽略,对涌浪高度有影响;而一般滑坡进入开阔水面,湖面宽度往往远大于滑坡宽度时,滑坡宽度的影响可以忽略。
4、冰崩体可能在冰湖上方垂直落入冰湖,而滑坡和崩塌很少有这种现象,即使有,其体积往往很小。
5、滑坡的传播过程中的衰减计算在传播开始时的短距离,运动距离与冰湖水深之比小于1时,会出现非常大的误差,甚至于无穷大,如刚开始时,运动距离约等于0时,计算方法有缺陷。
申请号为CN201610630522.9的中国专利文献公开了一种冰湖溃决型泥石流的防治方法及其应用,该方法结合流域的地形地貌特征参数、降雨特征参数、以及流域内物源特征参数和主河的输移能力等,确定溃决洪水泥石流的洪峰流量阈值;然后基于对冰碛坝体的现场调查、取样、室内测试获得的基本参数,判断每个计算区域内的冰碛坝体的稳定性并进行危险性分类;在此基础上确定冰碛坝体中桩柱结构体的体型参数及其分布方式。该专利文献公开的冰湖溃决型泥石流的防治方法及其应用,针对冰碛坝内埋藏冰消融、冰滑坡和冰崩入湖导致的冰湖突然溃决易形成特大型泥石流灾害,能有效控制坝体下切速度,延长溃决洪水下泄时间,削弱溃决洪水峰值流量,减小和调控溃决洪水形成泥石流的规模。但是,由于不能计算冰崩在冰湖中形成的涌浪,因而不能从根本上判断冰湖是否溃决,不能提供相应依据,致使主动防灾性较差。
发明内容
本发明为了克服上述现有技术的缺陷,提供一种冰湖涌浪高度的计算方法及其应用,本发明建立了更精确的冰崩在冰湖中形成的涌浪计算方法,为判断冰湖是否溃决提供依据,极大的提高了主动防灾适用性。
本发明通过下述技术方案实现:
一种冰湖涌浪高度的计算方法,其特征在于,包括以下步骤:
a、通过冰崩体的位置,区别冰川的属性;大陆性冰川或海洋性冰川气候干燥、降雪量少与负温低而雪线又高,收入少支出也少,活动性弱,冰舌短,冰川地质地貌作用弱,这属于大陆性冰川;反之,属于海洋性冰川;若悬贴于山坡上而不下降到山麓,冰川规模小,冰体厚度薄,冰川面积不足1km2,则为悬冰川;通过调查测绘确定冰崩体的基本数据,包括冰崩后端冰裂缝的位置、冰崩的前端位置、冰崩的平均长度L和冰崩的平均宽度W,计算冰崩面积A和冰崩体积V;测量冰川的长度与冰川的宽度,计算冰川面积,并由冰川面积计算潜在冰崩体的厚度D;
对于海洋性冰川
D=5.2+15.4F0.5 式一
对于大陆性冰川
D=-11.32+53.21F0.3 式二
对于悬冰川
D=34.4F0.45 式三
A=LW 式四
V=AD 式五
式中,D为潜在冰崩体的厚度,单位m;F为冰川面积,单位km2;A为冰崩面积,单位m2;L为冰崩的平均长度,单位m;W为冰崩的平均宽度,单位m;V为冰崩体积,单位m3;
b、确定冰崩在冰湖中形成涌浪高度H的参数,包括冰崩体与冰湖湖面的高程差H0、冰崩体与冰湖末端之间的距离J、冰湖末端宽度b、冰湖面积Ag、冰湖平均水深h和冰湖末端前地表坡度P;冰崩体与冰湖湖面的高程差H0是指冰崩体的质心和冰湖湖面的高程之差,通过调查测绘确定冰崩体及冰崩体质心位置与高程及冰湖湖面高程后获得;冰崩体与冰湖末端之间的距离J是指冰崩体从冰崩体质心位置运动进入冰湖之前的距离,通过调查测绘确定冰崩体及冰崩体质心位置从上到下顺冰崩运动到冰湖之前的路程后获得;冰湖末端前地表坡度P是指冰湖上游临近冰湖段地表坡度,通过调查测绘确定;冰湖末端宽度b是指冰湖最靠近上游端的湖面宽度,通过调查测绘测量;冰湖面积Ag通过调查测绘测量;冰湖平均水深h通过式六计算;
h=0.104Ag0.42 式六
式中,H为涌浪高度,单位m;H0为冰崩体与冰湖湖面的高程差,单位m;J为冰崩体与冰湖末端之间的距离,单位km;b为冰湖末端宽度,单位m;Ag为冰湖面积Ag,单位m2;h为冰湖平均水深,单位m;P为冰湖末端前地表坡度,单位rad;
c、冰崩在冰湖中形成涌浪高度H通过式七计算获得;
H=c0c1E0.4W0.5[1-0.15(x/h)0.3]/[(γg)0.4b0.5P0.55h0.6] 式七
式中:
c0为相对水深系数,当h/D>0.7时,c0=0.16;当h/D<0.7时,c0=0.12;
c1为冰崩体与冰湖距离系数,通过式八计算;
c1=1-0.2J0.5 式八
E为冰崩体势能,单位J,通过式九计算;
E=ρVgH0 式九
ρ为冰的比重,ρ=900kg/m3,g为重力加速度,g=9.81m/s2;
X为涌浪传播距离,是指从冰湖上游端开始往下游的距离,单位m;
γ为水的比重,γ=1000kg/m3。
本发明适用于冰崩体进入冰湖的滑入角度坡度0<P≤π/2的冰崩在冰湖中形成涌浪高度的计算。
进一步,本发明适用于冰崩体宽度与冰湖末端上游宽度之比W/b在0.05-1之间的冰崩在冰湖中形成涌浪高度的计算。
本发明的基本原理如下:
冰崩在冰湖中形成涌浪由于冰崩体巨大的动能,由冰崩体势能转化而来,巨大的入水面积推动冰湖中的水体形成涌浪,可以形成较高的涌浪高度并传播很远的距离。冰崩的主体是冰,因此其比重是固定的ρ=900kg/m3。冰崩体积越大,冰崩与冰湖的高差越大,其势能越大,在冰湖中形成的涌浪高度越大;冰湖上游地形越缓,坡度越小,进入冰湖的冰崩体沿水平方向运动速度越大,形成的涌浪高度越大;冰湖水深越大,冰崩体在冰湖内的影响范围越小,冰崩体在冰湖中形成的涌浪高度越小;冰湖的宽度越小,冰崩体在冰湖内所影响的范围越大,冰崩体进入冰湖推动和挤压湖水体越严重,冰崩体在冰湖中形成的涌浪高度越大;冰崩体在冰湖中形成的涌浪传播距离越远,其能量衰减越多,涌浪高度越小;当冰湖水深与冰崩体的厚度之比较小,在0.7以下时,冰崩体的巨大厚度作用在冰湖上并产生涌浪的作用有限,涌浪高度较小;冰崩体距离冰湖越远,进入冰湖形成涌浪的冰体越少,涌浪高度越小。
本发明在充分考虑冰崩体的尺度特征和冰崩体与冰湖之间的高差情况下,综合考虑冰湖及其周围地形特征,最终获得冰崩引起的冰湖涌浪高度计算方法。
本发明的有益效果主要表现在以下方面:
一、本发明,“a、通过冰崩体的位置,区别冰川的属性;通过调查测绘确定冰崩体的基本数据,包括冰崩后端冰裂缝的位置、冰崩的前端位置、冰崩的平均长度L和冰崩的平均宽度W,计算冰崩面积A和冰崩体积V;测量冰川的长度与冰川的宽度,计算冰川面积,并由冰川面积计算潜在冰崩体的厚度D;b、确定冰崩在冰湖中形成涌浪高度H的参数,包括冰崩体与冰湖湖面的高程差H0、冰崩体与冰湖末端之间的距离J、冰湖末端宽度b、冰湖面积Ag、冰湖平均水深h和冰湖末端前地表坡度P;c、冰崩在冰湖中形成涌浪高度H通过式七计算获得”,通过调查测绘确定冰崩体的长度、宽度和面积,由冰川面积计算冰崩体的厚度及体积,通过确定冰崩体的质心的高程,冰湖水面高程得到冰崩体与冰湖的高程差,能够计算出冰崩体进入冰湖时的势能,以及转化为进入冰湖的动能,较现有技术而言,通过考虑冰崩体的体积以及与冰湖的高差,冰崩体与冰湖之间的位置关系,冰湖后缘地形特征及冰湖形态特征,建立了更精确的冰崩在冰湖中形成的涌浪计算方法,为判断冰湖是否溃决提供依据,极大的提高了主动防灾适用性。
二、本发明,考虑到冰湖水深对涌浪高度的影响,通过全面考虑冰崩形成涌浪的冰湖水深影响因素,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
三、本发明,考虑到冰湖上游地形坡度对入水冰崩体运动和进入冰湖的角度的影响,通过全面考虑冰崩形成涌浪的地形坡度影响因素,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
四、本发明,考虑到冰湖上端冰湖宽度对冰崩体形成涌浪的影响,通过全面考虑冰崩形成涌浪的冰湖宽度影响因素,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
五、本发明,考虑到冰湖水深和冰崩体厚度的关系及其影响,通过全面考虑冰崩形成涌浪的冰崩体厚度及冰湖水深影响因素,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
六、本发明,考虑到冰崩体与冰湖之间的距离以及冰崩体运动中部分冰体会停留在运动通道上,通过全面考虑冰崩体的运动过程以及对涌浪的影响,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
七、本发明,考虑到涌浪在冰湖上传播时,传播距离对涌浪的影响,通过全面考虑冰崩形成涌浪的传播距离影响因素,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
具体实施方式
实施例1
一种冰湖涌浪高度的计算方法,包括以下步骤:
a、通过冰崩体的位置,区别冰川的属性;大陆性冰川或海洋性冰川气候干燥、降雪量少与负温低而雪线又高,收入少支出也少,活动性弱,冰舌短,冰川地质地貌作用弱,这属于大陆性冰川;反之,属于海洋性冰川;若悬贴于山坡上而不下降到山麓,冰川规模小,冰体厚度薄,冰川面积不足1km2,则为悬冰川;通过调查测绘确定冰崩体的基本数据,包括冰崩后端冰裂缝的位置、冰崩的前端位置、冰崩的平均长度L和冰崩的平均宽度W,计算冰崩面积A和冰崩体积V;测量冰川的长度与冰川的宽度,计算冰川面积,并由冰川面积计算潜在冰崩体的厚度D;
对于海洋性冰川
D=5.2+15.4F0.5 式一
对于大陆性冰川
D=-11.32+53.21F0.3 式二
对于悬冰川
D=34.4F0.45 式三
A=LW 式四
V=AD 式五
式中,D为潜在冰崩体的厚度,单位m;F为冰川面积,单位km2;A为冰崩面积,单位m2;L为冰崩的平均长度,单位m;W为冰崩的平均宽度,单位m;V为冰崩体积,单位m3;
b、确定冰崩在冰湖中形成涌浪高度H的参数,包括冰崩体与冰湖湖面的高程差H0、冰崩体与冰湖末端之间的距离J、冰湖末端宽度b、冰湖面积Ag、冰湖平均水深h和冰湖末端前地表坡度P;冰崩体与冰湖湖面的高程差H0是指冰崩体的质心和冰湖湖面的高程之差,通过调查测绘确定冰崩体及冰崩体质心位置与高程及冰湖湖面高程后获得;冰崩体与冰湖末端之间的距离J是指冰崩体从冰崩体质心位置运动进入冰湖之前的距离,通过调查测绘确定冰崩体及冰崩体质心位置从上到下顺冰崩运动到冰湖之前的路程后获得;冰湖末端前地表坡度P是指冰湖上游临近冰湖段地表坡度,通过调查测绘确定;冰湖末端宽度b是指冰湖最靠近上游端的湖面宽度,通过调查测绘测量;冰湖面积Ag通过调查测绘测量;冰湖平均水深h通过式六计算;
h=0.104Ag0.42 式六
式中,H为涌浪高度,单位m;H0为冰崩体与冰湖湖面的高程差,单位m;J为冰崩体与冰湖末端之间的距离,单位km;b为冰湖末端宽度,单位m;Ag为冰湖面积Ag,单位m2;h为冰湖平均水深,单位m;P为冰湖末端前地表坡度,单位rad;
c、冰崩在冰湖中形成涌浪高度H通过式七计算获得;
H=c0c1E0.4W0.5[1-0.15(x/h)0.3]/[(γg)0.4b0.5P0.55h0.6] 式七
式中:
c0为相对水深系数,当h/D>0.7时,c0=0.16;当h/D<0.7时,c0=0.12;
c1为冰崩体与冰湖距离系数,通过式八计算;
c1=1-0.2J0.5 式八
E为冰崩体势能,单位J,通过式九计算;
E=ρVgH0 式九
ρ为冰的比重,ρ=900kg/m3,g为重力加速度,g=9.81m/s2;
X为涌浪传播距离,是指从冰湖上游端开始往下游的距离,单位m;
γ为水的比重,γ=1000kg/m3。
“a、通过冰崩体的位置,区别冰川的属性;通过调查测绘确定冰崩体的基本数据,包括冰崩后端冰裂缝的位置、冰崩的前端位置、冰崩的平均长度L和冰崩的平均宽度W,计算冰崩面积A和冰崩体积V;测量冰川的长度与冰川的宽度,计算冰川面积,并由冰川面积计算潜在冰崩体的厚度D;b、确定冰崩在冰湖中形成涌浪高度H的参数,包括冰崩体与冰湖湖面的高程差H0、冰崩体与冰湖末端之间的距离J、冰湖末端宽度b、冰湖面积Ag、冰湖平均水深h和冰湖末端前地表坡度P;c、冰崩在冰湖中形成涌浪高度H通过式七计算获得”,通过调查测绘确定冰崩体的长度、宽度和面积,由冰川面积计算冰崩体的厚度及体积,通过确定冰崩体的质心的高程,冰湖水面高程得到冰崩体与冰湖的高程差,能够计算出冰崩体进入冰湖时的势能,以及转化为进入冰湖的动能,较现有技术而言,通过考虑冰崩体的体积以及与冰湖的高差,冰崩体与冰湖之间的位置关系,冰湖后缘地形特征及冰湖形态特征,建立了更精确的冰崩在冰湖中形成的涌浪计算方法,为判断冰湖是否溃决提供依据,极大的提高了主动防灾适用性。
实施例2
一种冰湖涌浪高度的计算方法,包括以下步骤:
a、通过冰崩体的位置,区别冰川的属性;大陆性冰川或海洋性冰川气候干燥、降雪量少与负温低而雪线又高,收入少支出也少,活动性弱,冰舌短,冰川地质地貌作用弱,这属于大陆性冰川;反之,属于海洋性冰川;若悬贴于山坡上而不下降到山麓,冰川规模小,冰体厚度薄,冰川面积不足1km2,则为悬冰川;通过调查测绘确定冰崩体的基本数据,包括冰崩后端冰裂缝的位置、冰崩的前端位置、冰崩的平均长度L和冰崩的平均宽度W,计算冰崩面积A和冰崩体积V;测量冰川的长度与冰川的宽度,计算冰川面积,并由冰川面积计算潜在冰崩体的厚度D;
对于海洋性冰川
D=5.2+15.4F0.5 式一
对于大陆性冰川
D=-11.32+53.21F0.3 式二
对于悬冰川
D=34.4F0.45 式三
A=LW 式四
V=AD 式五
式中,D为潜在冰崩体的厚度,单位m;F为冰川面积,单位km2;A为冰崩面积,单位m2;L为冰崩的平均长度,单位m;W为冰崩的平均宽度,单位m;V为冰崩体积,单位m3;
b、确定冰崩在冰湖中形成涌浪高度H的参数,包括冰崩体与冰湖湖面的高程差H0、冰崩体与冰湖末端之间的距离J、冰湖末端宽度b、冰湖面积Ag、冰湖平均水深h和冰湖末端前地表坡度P;冰崩体与冰湖湖面的高程差H0是指冰崩体的质心和冰湖湖面的高程之差,通过调查测绘确定冰崩体及冰崩体质心位置与高程及冰湖湖面高程后获得;冰崩体与冰湖末端之间的距离J是指冰崩体从冰崩体质心位置运动进入冰湖之前的距离,通过调查测绘确定冰崩体及冰崩体质心位置从上到下顺冰崩运动到冰湖之前的路程后获得;冰湖末端前地表坡度P是指冰湖上游临近冰湖段地表坡度,通过调查测绘确定;冰湖末端宽度b是指冰湖最靠近上游端的湖面宽度,通过调查测绘测量;冰湖面积Ag通过调查测绘测量;冰湖平均水深h通过式六计算;
h=0.104Ag0.42 式六
式中,H为涌浪高度,单位m;H0为冰崩体与冰湖湖面的高程差,单位m;J为冰崩体与冰湖末端之间的距离,单位km;b为冰湖末端宽度,单位m;Ag为冰湖面积Ag,单位m2;h为冰湖平均水深,单位m;P为冰湖末端前地表坡度,单位rad;
c、冰崩在冰湖中形成涌浪高度H通过式七计算获得;
H=c0c1E0.4W0.5[1-0.15(x/h)0.3]/[(γg)0.4b0.5P0.55h0.6] 式七
式中:
c0为相对水深系数,当h/D>0.7时,c0=0.16;当h/D<0.7时,c0=0.12;
c1为冰崩体与冰湖距离系数,通过式八计算;
c1=1-0.2J0.5 式八
E为冰崩体势能,单位J,通过式九计算;
E=ρVgH0 式九
ρ为冰的比重,ρ=900kg/m3,g为重力加速度,g=9.81m/s2;
X为涌浪传播距离,是指从冰湖上游端开始往下游的距离,单位m;
γ为水的比重,γ=1000kg/m3。
本发明适用于冰崩体进入冰湖的滑入角度坡度0<P≤π/2的冰崩在冰湖中形成涌浪高度的计算。
实施例3
一种冰湖涌浪高度的计算方法,包括以下步骤:
a、通过冰崩体的位置,区别冰川的属性;大陆性冰川或海洋性冰川气候干燥、降雪量少与负温低而雪线又高,收入少支出也少,活动性弱,冰舌短,冰川地质地貌作用弱,这属于大陆性冰川;反之,属于海洋性冰川;若悬贴于山坡上而不下降到山麓,冰川规模小,冰体厚度薄,冰川面积不足1km2,则为悬冰川;通过调查测绘确定冰崩体的基本数据,包括冰崩后端冰裂缝的位置、冰崩的前端位置、冰崩的平均长度L和冰崩的平均宽度W,计算冰崩面积A和冰崩体积V;测量冰川的长度与冰川的宽度,计算冰川面积,并由冰川面积计算潜在冰崩体的厚度D;
对于海洋性冰川
D=5.2+15.4F0.5 式一
对于大陆性冰川
D=-11.32+53.21F0.3 式二
对于悬冰川
D=34.4F0.45 式三
A=LW 式四
V=AD 式五
式中,D为潜在冰崩体的厚度,单位m;F为冰川面积,单位km2;A为冰崩面积,单位m2;L为冰崩的平均长度,单位m;W为冰崩的平均宽度,单位m;V为冰崩体积,单位m3;
b、确定冰崩在冰湖中形成涌浪高度H的参数,包括冰崩体与冰湖湖面的高程差H0、冰崩体与冰湖末端之间的距离J、冰湖末端宽度b、冰湖面积Ag、冰湖平均水深h和冰湖末端前地表坡度P;冰崩体与冰湖湖面的高程差H0是指冰崩体的质心和冰湖湖面的高程之差,通过调查测绘确定冰崩体及冰崩体质心位置与高程及冰湖湖面高程后获得;冰崩体与冰湖末端之间的距离J是指冰崩体从冰崩体质心位置运动进入冰湖之前的距离,通过调查测绘确定冰崩体及冰崩体质心位置从上到下顺冰崩运动到冰湖之前的路程后获得;冰湖末端前地表坡度P是指冰湖上游临近冰湖段地表坡度,通过调查测绘确定;冰湖末端宽度b是指冰湖最靠近上游端的湖面宽度,通过调查测绘测量;冰湖面积Ag通过调查测绘测量;冰湖平均水深h通过式六计算;
h=0.104Ag0.42 式六
式中,H为涌浪高度,单位m;H0为冰崩体与冰湖湖面的高程差,单位m;J为冰崩体与冰湖末端之间的距离,单位km;b为冰湖末端宽度,单位m;Ag为冰湖面积Ag,单位m2;h为冰湖平均水深,单位m;P为冰湖末端前地表坡度,单位rad;
c、冰崩在冰湖中形成涌浪高度H通过式七计算获得;
H=c0c1E0.4W0.5[1-0.15(x/h)0.3]/[(γg)0.4b0.5P0.55h0.6] 式七
式中:
c0为相对水深系数,当h/D>0.7时,c0=0.16;当h/D<0.7时,c0=0.12;
c1为冰崩体与冰湖距离系数,通过式八计算;
c1=1-0.2J0.5 式八
E为冰崩体势能,单位J,通过式九计算;
E=ρVgH0 式九
ρ为冰的比重,ρ=900kg/m3,g为重力加速度,g=9.81m/s2;
X为涌浪传播距离,是指从冰湖上游端开始往下游的距离,单位m;
γ为水的比重,γ=1000kg/m3。
进一步,本发明适用于冰崩体宽度与冰湖末端上游宽度之比W/b在0.05-1之间的冰崩在冰湖中形成涌浪高度的计算。
考虑到冰湖水深对涌浪高度的影响,通过全面考虑冰崩形成涌浪的冰湖水深影响因素,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
考虑到冰湖上游地形坡度对入水冰崩体运动和进入冰湖的角度的影响,通过全面考虑冰崩形成涌浪的地形坡度影响因素,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
考虑到冰湖上端冰湖宽度对冰崩体形成涌浪的影响,通过全面考虑冰崩形成涌浪的冰湖宽度影响因素,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
考虑到冰湖水深和冰崩体厚度的关系及其影响,通过全面考虑冰崩形成涌浪的冰崩体厚度及冰湖水深影响因素,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
考虑到冰崩体与冰湖之间的距离以及冰崩体运动中部分冰体会停留在运动通道上,通过全面考虑冰崩体的运动过程以及对涌浪的影响,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
考虑到涌浪在冰湖上传播时,传播距离对涌浪的影响,通过全面考虑冰崩形成涌浪的传播距离影响因素,以定量的方式精确计算出涌浪的高度,极大的保障了防灾适用性。
以室内模拟方法,采用本发明对冰崩在冰湖中产生的涌浪高度计算方法进行验证。
室内实验水池长8m,高0.5m,宽度有0.26m和0.5m两种。实验前在水池中蓄水水深h=0.15m、0.2m、0.25m、0.35m。采用接近长方体的塑料桶装水的方法,使其整体比重为900kg/m3,即与冰的比重相同。实验采用的塑料桶冰体有3种规格:冰体1:体积V=0.0058m3,宽度W=0.195m,厚度D=0.115m;冰体2:体积V=0.0113m3,宽度W=0.24m,厚度D=0.15m;冰体3:体积V=0.0237m3,宽度W=0.31m,厚度D=0.185m。实验时冰体在水池水面上方高差0.11m-0.14m范围,分别以15度和30度,弧度分别为0.262和0.524的坡度滑入水池中,测量冰体滑入水池中产生的涌浪高度以及沿传播方向1m,2m,3m位置的涌浪高度H’(m)。因为涌浪在水池远端的阻挡作用,涌浪高度在靠近水池远端会出现增加的现象,因此靠近远端的数据不能用于本发明验证,仅仅3m范围内数据可以用于验证。
表1为实验中的参数以及采用本发明的计算公式计算结果及误差表。
表1实验参数及计算结果和误差表
实验测试结果与本发明计算结果对比,最大误差约30%,有16个的误差在10%以内。可见本发明可以得到较准确的冰崩在冰湖内造成涌浪高度计算结果,能够为判断冰湖是否溃决提供依据,极大的提高了主动防灾适用性。
Claims (3)
1.一种冰湖涌浪高度的计算方法,其特征在于,包括以下步骤:
a、通过冰崩体的位置,区别冰川的属性;大陆性冰川或海洋性冰川气候干燥、降雪量少与负温低而雪线又高,收入少支出也少,活动性弱,冰舌短,冰川地质地貌作用弱,这属于大陆性冰川;反之,属于海洋性冰川;若悬贴于山坡上而不下降到山麓,冰川规模小,冰体厚度薄,冰川面积不足1km2,则为悬冰川;通过调查测绘确定冰崩体的基本数据,包括冰崩后端冰裂缝的位置、冰崩的前端位置、冰崩的平均长度L和冰崩的平均宽度W,计算冰崩面积A和冰崩体积V;测量冰川的长度与冰川的宽度,计算冰川面积,并由冰川面积计算潜在冰崩体的厚度D;
对于海洋性冰川
D=5.2+15.4F0.5 式一
对于大陆性冰川
D=-11.32+53.21F0.3 式二
对于悬冰川
D=34.4F0.45 式三
A=LW 式四
V=AD 式五
式中,D为潜在冰崩体的厚度,单位m;F为冰川面积,单位km2;A为冰崩面积,单位m2;L为冰崩的平均长度,单位m;W为冰崩的平均宽度,单位m;V为冰崩体积,单位m3;
b、确定冰崩在冰湖中形成涌浪高度H的参数,包括冰崩体与冰湖湖面的高程差H0、冰崩体与冰湖末端之间的距离J、冰湖末端宽度b、冰湖面积Ag、冰湖平均水深h和冰湖末端前地表坡度P;冰崩体与冰湖湖面的高程差H0是指冰崩体的质心和冰湖湖面的高程之差,通过调查测绘确定冰崩体及冰崩体质心位置与高程及冰湖湖面高程后获得;冰崩体与冰湖末端之间的距离J是指冰崩体从冰崩体质心位置运动进入冰湖之前的距离,通过调查测绘确定冰崩体及冰崩体质心位置从上到下顺冰崩运动到冰湖之前的路程后获得;冰湖末端前地表坡度P是指冰湖上游临近冰湖段地表坡度,通过调查测绘确定;冰湖末端宽度b是指冰湖最靠近上游端的湖面宽度,通过调查测绘测量;冰湖面积Ag通过调查测绘测量;冰湖平均水深h通过式六计算;
h=0.104Ag0.42 式六
式中,H为涌浪高度,单位m;H0为冰崩体与冰湖湖面的高程差,单位m;J为冰崩体与冰湖末端之间的距离,单位km;b为冰湖末端宽度,单位m;Ag为冰湖面积Ag,单位m2;h为冰湖平均水深,单位m;P为冰湖末端前地表坡度,单位rad;
c、冰崩在冰湖中形成涌浪高度H通过式七计算获得;
H=c0c1E0.4W0.5[1-0.15(x/h)0.3]/[(γg)0.4b0.5P0.55h0.6] 式七
式中:
c0为相对水深系数,当h/D>0.7时,c0=0.16;当h/D<0.7时,c0=0.12;
c1为冰崩体与冰湖距离系数,通过式八计算;
c1=1-0.2J0.5 式八
E为冰崩体势能,单位J,通过式九计算;
E=ρVgH0 式九
ρ为冰的比重,ρ=900kg/m3,g为重力加速度,g=9.81m/s2;
X为涌浪传播距离,是指从冰湖上游端开始往下游的距离,单位m;
γ为水的比重,γ=1000kg/m3。
2.根据权利要求1所述的一种冰湖涌浪高度的计算方法的应用,其特征在于:适用于冰崩体进入冰湖的滑入角度坡度0<P≤π/2的冰崩在冰湖中形成涌浪高度的计算。
3.根据权利要求1所述的一种冰湖涌浪高度的计算方法的应用,其特征在于:适用于冰崩体宽度与冰湖末端上游宽度之比W/b在0.05-1之间的冰崩在冰湖中形成涌浪高度的计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911023405.6A CN110909293B (zh) | 2019-10-25 | 2019-10-25 | 一种冰湖涌浪高度的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911023405.6A CN110909293B (zh) | 2019-10-25 | 2019-10-25 | 一种冰湖涌浪高度的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110909293A true CN110909293A (zh) | 2020-03-24 |
CN110909293B CN110909293B (zh) | 2022-04-26 |
Family
ID=69815646
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911023405.6A Active CN110909293B (zh) | 2019-10-25 | 2019-10-25 | 一种冰湖涌浪高度的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110909293B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118467892A (zh) * | 2024-07-11 | 2024-08-09 | 中国地质科学院探矿工艺研究所 | 一种沟槽型冰崩碎屑流运动距离及危险范围预测方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106250635A (zh) * | 2016-08-02 | 2016-12-21 | 中国科学院水利部成都山地灾害与环境研究所 | 一种冰湖溃决型泥石流的防治方法及其应用 |
RU2618494C1 (ru) * | 2015-11-30 | 2017-05-03 | Федеральное государственное бюджетное учреждение "ВЫСОКОГОРНЫЙ ГЕОФИЗИЧЕСКИЙ ИНСТИТУТ" (ФГБУ "ВГИ") | Способ оценки единовременного максимально возможного объема твердых селевых выносов в селевое русло реки при сходе селя |
CN107016185A (zh) * | 2017-03-31 | 2017-08-04 | 中国地质科学院探矿工艺研究所 | 一种冰碛湖溃决洪水峰值流量预测的计算方法 |
CN107564245A (zh) * | 2017-09-28 | 2018-01-09 | 成都理工大学 | 一种考虑降雨的冰湖溃决的预报方法 |
CN107749144A (zh) * | 2017-09-28 | 2018-03-02 | 成都理工大学 | 一种冰湖溃决的溢流水位预警方法及其应用 |
CN107784164A (zh) * | 2017-09-28 | 2018-03-09 | 成都理工大学 | 一种冰湖溃决的早期识别方法及其应用 |
CN108021764A (zh) * | 2017-12-15 | 2018-05-11 | 四川建筑职业技术学院 | 一种冰川泥石流起动流深的计算方法 |
CN108681651A (zh) * | 2018-05-29 | 2018-10-19 | 成都理工大学 | 一种无冰芯的冰碛堤冰湖溃决的早期识别方法 |
CN108733948A (zh) * | 2018-05-29 | 2018-11-02 | 成都理工大学 | 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 |
CN109697314A (zh) * | 2018-12-18 | 2019-04-30 | 国网西藏电力有限公司 | 一种冰崩的最大运动距离的计算方法及其应用 |
-
2019
- 2019-10-25 CN CN201911023405.6A patent/CN110909293B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2618494C1 (ru) * | 2015-11-30 | 2017-05-03 | Федеральное государственное бюджетное учреждение "ВЫСОКОГОРНЫЙ ГЕОФИЗИЧЕСКИЙ ИНСТИТУТ" (ФГБУ "ВГИ") | Способ оценки единовременного максимально возможного объема твердых селевых выносов в селевое русло реки при сходе селя |
CN106250635A (zh) * | 2016-08-02 | 2016-12-21 | 中国科学院水利部成都山地灾害与环境研究所 | 一种冰湖溃决型泥石流的防治方法及其应用 |
CN107016185A (zh) * | 2017-03-31 | 2017-08-04 | 中国地质科学院探矿工艺研究所 | 一种冰碛湖溃决洪水峰值流量预测的计算方法 |
CN107564245A (zh) * | 2017-09-28 | 2018-01-09 | 成都理工大学 | 一种考虑降雨的冰湖溃决的预报方法 |
CN107749144A (zh) * | 2017-09-28 | 2018-03-02 | 成都理工大学 | 一种冰湖溃决的溢流水位预警方法及其应用 |
CN107784164A (zh) * | 2017-09-28 | 2018-03-09 | 成都理工大学 | 一种冰湖溃决的早期识别方法及其应用 |
CN108021764A (zh) * | 2017-12-15 | 2018-05-11 | 四川建筑职业技术学院 | 一种冰川泥石流起动流深的计算方法 |
CN108681651A (zh) * | 2018-05-29 | 2018-10-19 | 成都理工大学 | 一种无冰芯的冰碛堤冰湖溃决的早期识别方法 |
CN108733948A (zh) * | 2018-05-29 | 2018-11-02 | 成都理工大学 | 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 |
CN109697314A (zh) * | 2018-12-18 | 2019-04-30 | 国网西藏电力有限公司 | 一种冰崩的最大运动距离的计算方法及其应用 |
Non-Patent Citations (1)
Title |
---|
阿如汉 等: "基于Fluent数值模拟的冰湖涌浪分析——以西藏聂拉木县嘉龙错为例", 《冰川冻土》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118467892A (zh) * | 2024-07-11 | 2024-08-09 | 中国地质科学院探矿工艺研究所 | 一种沟槽型冰崩碎屑流运动距离及危险范围预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110909293B (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109697314B (zh) | 一种冰崩的最大运动距离的计算方法及其应用 | |
CN107784164B (zh) | 一种冰湖溃决的早期识别方法及其应用 | |
CN108733948B (zh) | 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 | |
CN108681651B (zh) | 一种无冰芯的冰碛堤冰湖溃决的早期识别方法 | |
CN106202770B (zh) | 一种泥石流拦砂坝后回淤坡度的计算方法及应用 | |
Cui et al. | A catastrophic natural disaster chain of typhoon-rainstorm-landslide-barrier lake-flooding in Zhejiang Province, China | |
CN110909293B (zh) | 一种冰湖涌浪高度的计算方法 | |
CN110866330B (zh) | 一种考虑冰湖前地形特征的冰崩产生涌浪的计算方法 | |
Oliveira et al. | Morphological controls and statistical modelling of boulder transport by extreme storms | |
Tinkler | Conditions for the Entrainment of Cuboid Boulders in Bedrock Streams: An Historical Review of Literature with Respect to Recent Investigations Paul Carling Institute of Environmental and Natural Sciences, Lancaster University, United Kingdom | |
CN108221842B (zh) | 泥石流格子坝过坝流速、流量计算方法 | |
Smiraglia et al. | Epiglacial morphology | |
CN110779499B (zh) | 一种考虑冰湖特征的冰崩形成涌浪高度的计算方法 | |
CN110820668A (zh) | 一种冰崩坠落在冰湖中造成的涌浪高度的计算方法 | |
Gudmestad | Technical and economic challenges for Arctic Coastal settlements due to melting of ice and permafrost in the Arctic | |
Shim et al. | Analysis on sediment transport system in the east coast of Korea | |
CN109631727B (zh) | 一种波浪作用下防波堤块石沉陷深度预测方法 | |
Crapoulet et al. | Using LiDAR topographic data for identifying coastal areas of northern France vulnerable to sea-level rise | |
Yang et al. | Assessment on typhoon prevention capability of fishing port based on numerical model calculations | |
CN110836663A (zh) | 一种冰湖涌浪高度随传播距离的计算方法及应用 | |
CN107288093A (zh) | 一种泥石流弯道的最大超高计算方法及应用 | |
Oktaviariyadi et al. | Study Planning of Check Dam as Sediment Controller at Sumbersari UB Forest Area | |
Sanderson et al. | Sediment–laden ice measurements and observations, and implications for potential interactions of ice and large woody debris with tidal turbines in Minas Passage | |
Dewez et al. | Laser survey and mechanical modeling of chalky sea cliff collapse in Normandy, France | |
Chen et al. | Calculation of discharge of debris flows caused by moraine-dam failure at Midui Gully, Tibet, China |
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 |