CN107784164B - 一种冰湖溃决的早期识别方法及其应用 - Google Patents

一种冰湖溃决的早期识别方法及其应用 Download PDF

Info

Publication number
CN107784164B
CN107784164B CN201710897519.8A CN201710897519A CN107784164B CN 107784164 B CN107784164 B CN 107784164B CN 201710897519 A CN201710897519 A CN 201710897519A CN 107784164 B CN107784164 B CN 107784164B
Authority
CN
China
Prior art keywords
ice
lake
burst
ice lake
glacier
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.)
Active
Application number
CN201710897519.8A
Other languages
English (en)
Other versions
CN107784164A (zh
Inventor
余斌
范冬捷
刘清华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN201710897519.8A priority Critical patent/CN107784164B/zh
Publication of CN107784164A publication Critical patent/CN107784164A/zh
Application granted granted Critical
Publication of CN107784164B publication Critical patent/CN107784164B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

本发明公开了一种冰湖溃决的早期识别方法,包括以下步骤:a、确定溃决冰湖及周边未溃决冰湖的基本冰川数据;b、确定发生冰湖溃决的先决条件,一是潜在冰崩体质心至冰湖垂直高度H与潜在冰崩体质心至冰湖直线距离L之比r≥0.125;二是潜在冰崩体滑入冰湖造成的上涨水位H1大于冰湖溃决临界水位H2;c、计算冰湖溃决的冰川综合判别因子S;d、当S>0.53时,危险性大;当0.38<S≤0.53时,危险性中等;当S≤0.38时,危险性小。本发明通过地形因子对冰湖溃决发生程度进行内部机理研究,建立了冰湖溃决易发计算模型,全面考虑冰湖溃决的地形影响因素,以定量的方式精确划分冰湖溃决危险性区域,适用性强。

Description

一种冰湖溃决的早期识别方法及其应用
技术领域
本发明涉及到水利工程技术领域,尤其涉及一种冰湖溃决的早期识别方法及其应用。
背景技术
冰湖溃决是一种发生在高山冰川分布区的自然现象。冰湖溃决的发生往往是由于冰湖后援的冰川发生冰崩或冰滑坡,冰崩或冰滑坡进入冰湖使冰湖水位上涨,造成冰湖溃决。而产生冰崩或冰滑坡需要具备两个条件:其一,有利于发生冰崩或冰滑坡的冰川条件、冰湖条件和冰川与冰湖之间的关系条件;其二,一定的温度条件导致冰川内部产生裂隙并最终诱发冰崩或冰滑坡的发生,这是冰湖溃决的诱发因素。这些条件综合影响并决定冰川的稳定性。除了温度的诱发因素外,影响冰湖溃决的主要因素包括:冰湖后缘冰舌的坡度、冰川的坡向、冰舌面积、冰湖面积、潜在冰崩体面积与厚度、冰崩体质心至冰湖垂直高度、冰崩体质心至冰湖直线距离、冰湖溃口宽度、溃口厚度和冰湖面积。
目前,国内外对冰湖溃决形成的地形条件的研究主要集中于冰舌坡度的研究,是对其发育分布规律的简单数据统计,没有深入研究其内在机理(西藏泥石流与环境,1999:71-78;水文地质工程地质,2006,3:88-92;冰川冻土,2004,26:397-402;山地学报,2003,21:128-132;冰川冻土,1986,8:61-71;地理学报,1989,44:343-352;冰川冻土,1987,9:23-34;Can Geotech,2004,41:1068-1083;Nat Hazards,2007,41:131-157;QuaternaryInternational,2000,65:31-47.),研究成果很难用于其他区域。对于冰湖面积在冰湖溃决中的作用也有研究(西藏泥石流与环境,1999:71-78;第四纪研究,2003,23:621-628;):研究表明冰湖面积在1.89×105m3-6.3×105m3最有利于冰湖溃决,但这个研究没有考虑冰崩体积和冰湖面积之间的相互关系,这个范围仅表明是已发生的案例的总结,不代表真正的潜在冰湖溃决的可能性:如冰崩体非常巨大,超过这个面积范围,依然可以发生冰湖溃决。另外还有学者研究了冰舌至冰湖距离对冰湖溃决的影响(西藏泥石流与环境,1999:71-78),但没有考虑冰舌与冰湖的垂直高度差,如果高差很大,即使较远的距离,冰崩体仍然可能大量进入冰湖,如2016年7月发生在西藏的阿鲁错附近的冰崩运动距离超过5km,并进入冰湖,远大于以前研究认为的超过800m就不会有冰湖溃决发生,800m作为不会发生冰湖溃决的临界距离显然不准确。
公开号为CN105740616A,公开日为2016年07月06日的中国专利文献公开了一种冰湖溃决的危险性判断方法,其特征在于,包括以下步骤:通过谷歌地球确定溃决冰湖及周边未溃决冰湖的基本地形数据:冰湖后援冰舌坡度α,冰川的坡向θ;计算冰湖溃决的地形综合判别因子S:由冰舌坡度、冰川坡向:S=tan(α)+0.5sin(θ/2),以冰川的地形综合判别因子S对冰湖溃决的危险性等级进行划分。
该专利文献公开的冰湖溃决的危险性判断方法,通过研究地形因子对冰湖溃决发生程度进行内部机理研究,建立了冰湖溃决易发计算模型;但是,由于没有考虑冰川与冰湖的距离及高差,以及冰湖的特征等影响,致使冰湖溃决危险性划分准确度较低,影响防灾效果。
公开号为CN105513285A,公开日为2016年04月20日的中国专利文献公开了一种冰湖溃决预警方法,以潜在溃决冰湖为预警监测区域,测量确定潜在溃决冰湖的冰舌坡度α和冰川坡向θ,查看气象资料确定潜在冰湖溃决点前两年的平均正积温Ty2、多年平均正积温Ty,实时监测某个时间之前的正积温T0及该时间前7天平均温度和T7,以所述预警监测区域冰湖溃决临界值Cr为监测值,根据监测值大小划分冰湖溃决的预警级别。
该专利文献公开的冰湖溃决预警方法,通过研究温度因子和地形因子对冰湖溃决发生程度进行内部机理研究,建立了冰湖溃决预警预报模型;但是,该方法仅考虑了冰川发生冰崩,进而可能引起冰湖溃决,还没有考虑到冰崩体能否到达冰湖;冰崩体积是否足够大这些因素,因而其预警准确度低,适用性差。
发明内容
本发明为了克服上述现有技术的缺陷,提供一种冰湖溃决的早期识别方法及其应用,本发明通过地形因子对冰湖溃决发生程度进行内部机理研究,建立了冰湖溃决易发计算模型,全面考虑冰湖溃决的地形影响因素,以定量的方式精确划分冰湖溃决危险性区域,适用性强。
本发明通过下述技术方案实现:
一种冰湖溃决的早期识别方法,其特征在于,包括以下步骤:
a、通过谷歌地球确定溃决冰湖及周边未溃决冰湖的基本冰川数据,包括冰湖后缘冰舌坡度α,冰川的坡向θ、潜在冰崩体质心至冰湖垂直高度H、潜在冰崩体质心至冰湖直线距离L、溃决口宽度b、溃决口厚度T、溃决前冰湖宽度w0、溃决前冰湖长宽比C、溃决前冰湖岸边斜坡β和连续性冰川的面积F;
b、确定发生冰湖溃决的两个先决条件,其一是潜在冰崩体质心至冰湖垂直高度H与潜在冰崩体质心至冰湖直线距离L之比r≥0.125,若r<0.125,则冰湖溃决的危险性极小;其二是潜在冰崩体滑入冰湖造成的上涨水位H1大于冰湖溃决临界水位H2,若潜在冰崩体滑入冰湖造成的上涨水位H1小于冰湖溃决临界水位H2,则冰湖溃决的危险性极小;冰湖溃决临界水位H2由式1计算;潜在冰崩体滑入冰湖造成的上涨水位H1由式2计算;最大冰崩体体积V由式3计算;
Figure GDA0002788004340000031
其中,H2为冰湖溃决临界水位,m;d为起动泥沙特征粒径,m;T为溃决口厚度,m;b为溃决口宽度,m;w0为溃决前冰湖宽度,m;C为溃决前冰湖长宽比;β为溃决前冰湖岸边斜坡,度;
H1=0.45V/A′ 式2
其中,H1为潜在冰崩体滑入冰湖造成的上涨水位,m;A′为平均冰湖水面面积,m2;V为最大冰崩体体积,m3
V=A1*h 式3
其中,h为冰川的厚度,m;A1为冰舌所能引发最大冰崩体的面积,m2
平均冰湖水面面积A′由式8计算;
Figure GDA0002788004340000032
c、通过式9计算冰湖溃决的冰川综合判别因子S;
S=tan(α)+0.5sin(θ/2) 式9
其中,S为冰湖溃决的冰川综合判别因子;α为冰湖后缘冰舌坡度,α≥3°;θ为冰川的坡向,0°<θ≤360°;
d、在满足两个先决条件下,以冰湖溃决的冰川综合判别因子S对冰湖溃决的危险性等级进行划分,当S>0.53时,危险性大;当0.38<S≤0.53时,危险性中等;当S≤0.38时,危险性小。
所述步骤b中,冰舌所能引发最大冰崩体的面积A1由式7计算;
A1=0.035A2+13000 式7
其中,A2为冰湖面积,m2;A2≥15000m2
所述步骤b中,冰川为海洋冰川时,冰川的厚度h由式4计算;
h=5.2+15.4F0.5 式4
其中,F为连续性冰川的面积,km2
所述步骤b中,冰川为大陆冰川时,冰川的厚度h由式5计算;
h=-11.32+53.21F0.3 式5
其中,F为连续性冰川的面积,km2
所述步骤b中,冰川为悬冰川时,冰川的厚度h由式6计算;
h=34.4F0.45 式6
其中,F为连续性冰川的面积,km2;F<1km2
本发明,适用于冰川发生冰崩且冰崩体进入冰湖引发冰湖溃决的早期识别。
本发明的基本原理如下:
冰湖后缘冰舌坡度α是影响冰崩和冰滑坡发生的最主要因素,冰湖后缘冰舌坡度α的陡缓不仅影响着冰体的运动形式,而且影响着冰舌段的汇流条件,大多数冰湖溃决发生在冰湖后缘冰舌坡度α为3°以上的坡度上。随着气候的干湿冷暖波动,在气候转向湿冷或干热年代,随着年平均温度的不断上升,冰川的温度也慢慢提高,冻结粘滞力慢慢变小,冰川的塑性增强,运动速度逐渐加快。特别在冰川前部和冰舌地段,由于海拔较低,温度较高,冰川温度相应较上部要高,塑性更强,运动速度更快。当冰川的活动性提高到一定水平,盛夏或秋天的消融水流沿着冰裂隙、孔隙下渗,缓慢地暖切割冰体,润滑底床,不断地减小冰川向下前进运动的阻力。冰舌地段如果坡度很大,则以冰崩形式发生;冰舌地段如果坡度稍小,则以快速运动的冰滑坡形式发生。由于冰川的坡向θ不同,日照和辐射条件各异,导致气温有明显的差异。在我国,一般地,位于北方的温度最低,因此在地形上,冰川的坡向θ越靠近南方越有利于冰川内部裂隙的产生;反之,越靠近北方,越不利于冰川内部裂隙的产生,也就不利于冰崩或冰滑坡的产生。潜在冰崩体质心至冰湖垂直高度H大,冰崩体具有大的势能,发生冰崩或冰滑坡时,冰崩体具有足够的能量进入冰湖;潜在冰崩体质心至冰湖直线距离L大,发生冰崩或冰滑坡时,冰崩体在蹦滑的过程中消耗的能量大,冰崩体有可能在未进入冰湖时就已能量耗尽,停止运动,因此潜在冰崩体质心至冰湖垂直高度H与潜在冰崩体质心至冰湖直线距离L之比r≥0.125才能使冰崩体进入冰湖,从而引发冰湖溃决。式9将这几个影响因素完整地结合在一起,综合地考虑了地形因素的作用,体现出各影响因子相互关系和重要性。
冰湖的溃决条件中,气象条件与冰川条件,以及冰川与冰湖之间的关系,还有冰湖条件等条件缺一不可。冰川与冰湖之间的关系,即冰川到冰湖之间的距离与高差,是决定冰崩体能否到达冰湖的关键,如果冰崩体不能到达冰湖,不会引起冰湖水位上涨,也不会发生冰湖溃决。冰湖条件是决定冰湖溃决的另外一个重要的先决条件:冰湖面积过大(与冰崩体体积相比较),冰湖水位上涨太小,不能侵蚀冰湖溢流口的冰碛堤,形成溃决;冰湖面积过小(与溃决口宽度和厚度相比较),冰湖水位上涨后又迅速回落,还没有将溃决口全部侵蚀就回落到原来的水位,也不能形成溃决。因此,除了诱发冰湖溃决的气象条件与冰川条件外,冰川与冰湖之间的关系、冰湖条件都是冰湖溃决的重要先决条件。
在此机理研究基础上,通过大量研究与准确的推导与计算,提出本发明的方法与判断标准。
本发明的有益效果是:
本发明,通过地形因子对冰湖溃决发生程度的内部机理进行研究,建立了冰湖溃决易发计算模型;全面考虑了冰湖后缘冰舌坡度α,冰川的坡向θ、潜在冰崩体质心至冰湖垂直高度H、潜在冰崩体质心至冰湖直线距离L、溃决口宽度b、溃决口厚度T、溃决前冰湖宽度w0、溃决前冰湖长宽比C、溃决前冰湖岸边斜坡β和连续性冰川的面积F这些影响因素,提出两个发生冰湖溃决的先决条件,不满足先决条件的冰湖溃决危险性极小;并以定量的方式精确划分冰湖溃决危险性区域;式9中S、tan(α)、sin(θ/2)均为无量纲参数,在各种冰湖溃决条件下都可以使用,具有相当高的适用性;由于充分考虑了冰川特征、冰湖特征以及冰川与冰湖之间的距离与高差等的关系,完整地引入了各参数对冰湖溃决的影响,因此对冰湖溃决的早期识别更加准确。
具体实施方式
实施例1
一种冰湖溃决的早期识别方法,包括以下步骤:
a、通过谷歌地球确定溃决冰湖及周边未溃决冰湖的基本冰川数据,包括冰湖后缘冰舌坡度α,冰川的坡向θ、潜在冰崩体质心至冰湖垂直高度H、潜在冰崩体质心至冰湖直线距离L、溃决口宽度b、溃决口厚度T、溃决前冰湖宽度w0、溃决前冰湖长宽比C、溃决前冰湖岸边斜坡β和连续性冰川的面积F;
b、确定发生冰湖溃决的两个先决条件,其一是潜在冰崩体质心至冰湖垂直高度H与潜在冰崩体质心至冰湖直线距离L之比r≥0.125,若r<0.125,则冰湖溃决的危险性极小;其二是潜在冰崩体滑入冰湖造成的上涨水位H1大于冰湖溃决临界水位H2,若潜在冰崩体滑入冰湖造成的上涨水位H1小于冰湖溃决临界水位H2,则冰湖溃决的危险性极小;冰湖溃决临界水位H2由式1计算;潜在冰崩体滑入冰湖造成的上涨水位H1由式2计算;最大冰崩体体积V由式3计算;
Figure GDA0002788004340000051
其中,H2为冰湖溃决临界水位,m;d为起动泥沙特征粒径,m;T为溃决口厚度,m;b为溃决口宽度,m;w0为溃决前冰湖宽度,m;C为溃决前冰湖长宽比;β为溃决前冰湖岸边斜坡,度;
H1=0.45V/A′ 式2
其中,H1为潜在冰崩体滑入冰湖造成的上涨水位,m;A′为平均冰湖水面面积,m2;V为最大冰崩体体积,m3
V=A1*h 式3
其中,h为冰川的厚度,m;A1为冰舌所能引发最大冰崩体的面积,m2
平均冰湖水面面积A′由式8计算;
Figure GDA0002788004340000061
c、通过式9计算冰湖溃决的冰川综合判别因子S;
S=tan(α)+0.5sin(θ/2) 式9
其中,S为冰湖溃决的冰川综合判别因子;α为冰湖后缘冰舌坡度,α≥3°;θ为冰川的坡向,0°<θ≤360°;
d、在满足两个先决条件下,以冰湖溃决的冰川综合判别因子S对冰湖溃决的危险性等级进行划分,当S>0.53时,危险性大;当0.38<S≤0.53时,危险性中等;当S≤0.38时,危险性小。
通过地形因子对冰湖溃决发生程度的内部机理进行研究,建立了冰湖溃决易发计算模型;全面考虑了冰湖后缘冰舌坡度α,冰川的坡向θ、潜在冰崩体质心至冰湖垂直高度H、潜在冰崩体质心至冰湖直线距离L、溃决口宽度b、溃决口厚度T、溃决前冰湖宽度w0、溃决前冰湖长宽比C、溃决前冰湖岸边斜坡β和连续性冰川的面积F这些影响因素,提出两个发生冰湖溃决的先决条件,不满足先决条件的冰湖溃决危险性极小;并以定量的方式精确划分冰湖溃决危险性区域;式9中S、tan(α)、sin(θ/2)均为无量纲参数,在各种冰湖溃决条件下都可以使用,具有相当高的适用性;由于充分考虑了冰川特征、冰湖特征以及冰川与冰湖之间的距离与高差等的关系,完整地引入了各参数对冰湖溃决的影响,因此对冰湖溃决的早期识别更加准确。
实施例2
一种冰湖溃决的早期识别方法,包括以下步骤:
a、通过谷歌地球确定溃决冰湖及周边未溃决冰湖的基本冰川数据,包括冰湖后缘冰舌坡度α,冰川的坡向θ、潜在冰崩体质心至冰湖垂直高度H、潜在冰崩体质心至冰湖直线距离L、溃决口宽度b、溃决口厚度T、溃决前冰湖宽度w0、溃决前冰湖长宽比C、溃决前冰湖岸边斜坡β和连续性冰川的面积F;
b、确定发生冰湖溃决的两个先决条件,其一是潜在冰崩体质心至冰湖垂直高度H与潜在冰崩体质心至冰湖直线距离L之比r≥0.125,若r<0.125,则冰湖溃决的危险性极小;其二是潜在冰崩体滑入冰湖造成的上涨水位H1大于冰湖溃决临界水位H2,若潜在冰崩体滑入冰湖造成的上涨水位H1小于冰湖溃决临界水位H2,则冰湖溃决的危险性极小;冰湖溃决临界水位H2由式1计算;潜在冰崩体滑入冰湖造成的上涨水位H1由式2计算;最大冰崩体体积V由式3计算;
Figure GDA0002788004340000071
其中,H2为冰湖溃决临界水位,m;d为起动泥沙特征粒径,m;T为溃决口厚度,m;b为溃决口宽度,m;w0为溃决前冰湖宽度,m;C为溃决前冰湖长宽比;β为溃决前冰湖岸边斜坡,度;
H1=0.45V/A′ 式2
其中,H1为潜在冰崩体滑入冰湖造成的上涨水位,m;A′为平均冰湖水面面积,m2;V为最大冰崩体体积,m3
V=A1*h 式3
其中,h为冰川的厚度,m;A1为冰舌所能引发最大冰崩体的面积,m2
平均冰湖水面面积A′由式8计算;
Figure GDA0002788004340000072
c、通过式9计算冰湖溃决的冰川综合判别因子S;
S=tan(α)+0.5sin(θ/2) 式9
其中,S为冰湖溃决的冰川综合判别因子;α为冰湖后缘冰舌坡度,α≥3°;θ为冰川的坡向,0°<θ≤360°;
d、在满足两个先决条件下,以冰湖溃决的冰川综合判别因子S对冰湖溃决的危险性等级进行划分,当S>0.53时,危险性大;当0.38<S≤0.53时,危险性中等;当S≤0.38时,危险性小。
所述步骤b中,冰舌所能引发最大冰崩体的面积A1由式7计算;
A1=0.035A2+13000 式7
其中,A2为冰湖面积,m2;A2≥15000m2
所述步骤b中,冰川为海洋冰川时,冰川的厚度h由式4计算;
h=5.2+15.4F0.5 式4
其中,F为连续性冰川的面积,km2
所述步骤b中,冰川为大陆冰川时,冰川的厚度h由式5计算;
h=-11.32+53.21F0.3 式5
其中,F为连续性冰川的面积,km2
所述步骤b中,冰川为悬冰川时,冰川的厚度h由式6计算;
h=34.4F0.45 式6
其中,F为连续性冰川的面积,km2;F<1km2
本发明,适用于冰川发生冰崩且冰崩体进入冰湖引发冰湖溃决的早期识别。
下面采用本发明方法对西藏光谢错附近的4处冰湖以及西藏吉莱错附近的5处冰湖溃决危险性进行分析:
分别量取各冰湖的冰湖后缘冰舌坡度α、冰川的坡向θ、潜在冰崩体质心至冰湖垂直高度H、潜在冰崩体质心至冰湖直线距离L、冰舌所能引发最大冰崩体的面积A1、冰湖面积A2、溃决口宽度b、溃决口厚度T、溃决前冰湖宽度w0、溃决前冰湖长宽比C、溃决前冰湖岸边斜坡β和连续性冰川的面积F。
根据式9计算出各冰湖的冰湖溃决的冰川综合判别因子S,根据式2—式8计算出各冰湖的潜在冰崩体滑入冰湖造成的上涨水位H1,根据式1计算出各冰湖的冰湖溃决临界水位H2
9处冰湖地形因素测量值及S、H1、H2、H/L值计算情况、危险性分级以及冰湖溃决的实际发生情况如表1所示。
编号 H<sub>2</sub> H<sub>1</sub> H<sub>1</sub>≥H<sub>2</sub> α θ r r≥0.125 S 危险等级 溃决
光谢错 1.43 11.17 20 330 0.14 0.49 中等
G1 1.39 1.29 7 30 0.28 0.25 极小
G2 1.43 1.40 12 280 0.19 0.53 极小
G3 1.31 7.76 9 40 0.13 0.33
吉莱错 1.33 1.68 7 250 0.23 0.53
J1 1.38 2.76 9 220 0.12 0.63 极小
J2 1.34 9.09 5 30 0.15 0.22
J3 1.48 5.37 7 0 0.11 0.12 极小
J4 1.45 5.26 7 90 0.13 0.48
表1
首先,根据先决条件1,表1编号中J1、J3两个冰湖潜在冰崩体质心至冰湖垂直高度H与潜在冰崩体质心至冰湖直线距离L之比r<0.125,虽然J1冰湖溃决的冰川综合判别因子S很大,但冰崩体很难进入冰湖,故J1、J3冰湖危险性极小。
其次,根据先决条件2,表1编号中G1、G2两个冰湖冰崩造成的潜在冰崩体滑入冰湖造成的上涨水位H1小于冰湖溃决临界水位H2,虽然G2冰湖溃决的冰川综合判别因子S值很大,但很难发生大的冰崩使冰湖产生较大的上涨水位,故G1、G2冰湖危险性极小。
最后,在满足两个先决条件的基础上,根据危险性划分标准进行危险性划分:当S>0.53时,危险性大;当0.38<S≤0.53时,危险性中等;当S≤0.38时,危险性小。
表1中光谢错危险等级中等,发生了冰湖溃决;表1中,吉莱错危险等级大,发生了溃决;表1中G3、J2、J4危险等级小,都没有发生冰湖溃决;表1中G1、G2、J1、J3危险等级极小,都没有发生冰湖溃决。
可见,应用本发明方法对冰湖溃决危险性划分的准确性较高。

Claims (6)

1.一种冰湖溃决的早期识别方法,其特征在于,包括以下步骤:
a、通过谷歌地球确定溃决冰湖及周边未溃决冰湖的基本冰川数据,包括冰湖后缘冰舌坡度α,冰川的坡向θ、潜在冰崩体质心至冰湖垂直高度H、潜在冰崩体质心至冰湖直线距离L、溃决口宽度b、溃决口厚度T、溃决前冰湖宽度w0、溃决前冰湖长宽比C、溃决前冰湖岸边斜坡β和连续性冰川的面积F;
b、确定发生冰湖溃决的两个先决条件,其一是潜在冰崩体质心至冰湖垂直高度H与潜在冰崩体质心至冰湖直线距离L之比r≥0.125,若r<0.125,则冰湖溃决的危险性极小;其二是潜在冰崩体滑入冰湖造成的上涨水位H1大于冰湖溃决临界水位H2,若潜在冰崩体滑入冰湖造成的上涨水位H1小于冰湖溃决临界水位H2,则冰湖溃决的危险性极小;冰湖溃决临界水位H2由式1计算;潜在冰崩体滑入冰湖造成的上涨水位H1由式2计算;最大冰崩体体积V由式3计算;
Figure FDA0002788004330000011
其中,H2为冰湖溃决临界水位,m;d为起动泥沙特征粒径,m;T为溃决口厚度,m;b为溃决口宽度,m;w0为溃决前冰湖宽度,m;C为溃决前冰湖长宽比;β为溃决前冰湖岸边斜坡,度;
H1=0.45V/A′ 式2
其中,H1为潜在冰崩体滑入冰湖造成的上涨水位,m;A′为平均冰湖水面面积,m2;V为最大冰崩体体积,m3
V=A1*h 式3
其中,h为冰川的厚度,m;A1为冰舌所能引发最大冰崩体的面积,m2
平均冰湖水面面积A′由式8计算;
Figure FDA0002788004330000012
c、通过式9计算冰湖溃决的冰川综合判别因子S;
S=tan(α)+0.5sin(θ/2) 式9
其中,S为冰湖溃决的冰川综合判别因子;α为冰湖后缘冰舌坡度,α≥3°;θ为冰川的坡向,0°<θ≤360°;
d、在满足两个先决条件下,以冰湖溃决的冰川综合判别因子S对冰湖溃决的危险性等级进行划分,当S>0.53时,危险性大;当0.38<S≤0.53时,危险性中等;当S≤0.38时,危险性小。
2.根据权利要求1所述的一种冰湖溃决的早期识别方法,其特征在于:所述步骤b中,冰舌所能引发最大冰崩体的面积A1由式7计算;
A1=0.035A2+13000 式7
其中,A2为冰湖面积,m2;A2≥15000m2
3.根据权利要求1所述的一种冰湖溃决的早期识别方法,其特征在于:所述步骤b中,冰川为海洋冰川时,冰川的厚度h由式4计算;
h=5.2+15.4F0.5 式4
其中,F为连续性冰川的面积,km2
4.根据权利要求1所述的一种冰湖溃决的早期识别方法,其特征在于:所述步骤b中,冰川为大陆冰川时,冰川的厚度h由式5计算;
h=-11.32+53.21F0.3 式5
其中,F为连续性冰川的面积,km2
5.根据权利要求1所述的一种冰湖溃决的早期识别方法,其特征在于:所述步骤b中,冰川为悬冰川时,冰川的厚度h由式6计算;
h=34.4F0.45 式6
其中,F为连续性冰川的面积,km2;F<1km2
6.根据权利要求1所述的一种冰湖溃决的早期识别方法的应用,其特征在于:适用于冰川发生冰崩且冰崩体进入冰湖引发冰湖溃决的早期识别。
CN201710897519.8A 2017-09-28 2017-09-28 一种冰湖溃决的早期识别方法及其应用 Active CN107784164B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710897519.8A CN107784164B (zh) 2017-09-28 2017-09-28 一种冰湖溃决的早期识别方法及其应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710897519.8A CN107784164B (zh) 2017-09-28 2017-09-28 一种冰湖溃决的早期识别方法及其应用

Publications (2)

Publication Number Publication Date
CN107784164A CN107784164A (zh) 2018-03-09
CN107784164B true CN107784164B (zh) 2021-02-12

Family

ID=61433952

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710897519.8A Active CN107784164B (zh) 2017-09-28 2017-09-28 一种冰湖溃决的早期识别方法及其应用

Country Status (1)

Country Link
CN (1) CN107784164B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108733948B (zh) * 2018-05-29 2021-11-19 成都理工大学 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用
CN108681651B (zh) * 2018-05-29 2021-10-29 成都理工大学 一种无冰芯的冰碛堤冰湖溃决的早期识别方法
CN109472445B (zh) * 2018-09-29 2022-04-26 成都理工大学 冰水泥石流的危险性判断方法及其应用
CN109697314B (zh) * 2018-12-18 2022-11-01 国网西藏电力有限公司 一种冰崩的最大运动距离的计算方法及其应用
CN110820668A (zh) * 2019-10-25 2020-02-21 成都理工大学 一种冰崩坠落在冰湖中造成的涌浪高度的计算方法
CN110909293B (zh) * 2019-10-25 2022-04-26 成都理工大学 一种冰湖涌浪高度的计算方法
CN110866330B (zh) * 2019-10-28 2023-03-31 成都理工大学 一种考虑冰湖前地形特征的冰崩产生涌浪的计算方法
CN115953705B (zh) * 2023-03-09 2023-05-16 昆明理工大学 一种冰川冰湖溃决型泥石流判识方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1015025A1 (ru) * 1980-04-23 1983-04-30 Sosnin Feliks Способ борьбы с гл циальными сел ми
CN105740616A (zh) * 2016-01-28 2016-07-06 成都理工大学 冰湖溃决的危险性判断方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1015025A1 (ru) * 1980-04-23 1983-04-30 Sosnin Feliks Способ борьбы с гл циальными сел ми
CN105740616A (zh) * 2016-01-28 2016-07-06 成都理工大学 冰湖溃决的危险性判断方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《喜马拉雅山冰湖溃决型泥石流灾害链研究》;陈宇棠;《中国优秀硕士学位论文全文数据库 基础科学辑》;20081015(第10期);摘要、正文第1-82页 *
《西藏终碛堤冰湖溃决冰川、温度条件的初步研究》;刘秧;《中国优秀硕士学位论文全文数据库 基础科学辑》;20170315(第3期);摘要、正文第1-29页 *

Also Published As

Publication number Publication date
CN107784164A (zh) 2018-03-09

Similar Documents

Publication Publication Date Title
CN107784164B (zh) 一种冰湖溃决的早期识别方法及其应用
CN107564245B (zh) 一种考虑降雨的冰湖溃决的预报方法
CN108681651B (zh) 一种无冰芯的冰碛堤冰湖溃决的早期识别方法
CN108733948B (zh) 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用
Glasser et al. Sedimentary facies and landform genesis at a temperate outlet glacier: Soler Glacier, North Patagonian Icefield
CN105513285B (zh) 冰湖溃决预警方法
Chigira et al. Landslides in weathered granitic rocks in Japan and Malaysia
CN108090670B (zh) 一种基于水流侵蚀能力的堰塞坝形成评估方法
CN105740616A (zh) 冰湖溃决的危险性判断方法
CN111680423B (zh) 一种定量稻田水土界面降水溅蚀力的方法与应用
Erskine et al. Role of a rheophyte in bench development on a sand‐bed river in southeast Australia
Ouakhir et al. Fluvial dynamic in Oued El Abid Basin: monitoring and quantification at an upstream river section in bin El Ouidane dam-2016/2017-(central high atlas/Morocco)
CN109920213B (zh) 基于降雨时程分布的临界雨量进行实时山洪预警的方法
Wohl et al. Episodic wood loading in a mountainous neotropical watershed
Glasser et al. Topographic controls on glacier sediment–landform associations around the temperate North Patagonian Icefield
Xia et al. Variations of soil hydraulic properties along granitic slopes in Benggang erosion areas
Brennand et al. Tunnel channel character and evolution in central southern Ontario
CN110866330B (zh) 一种考虑冰湖前地形特征的冰崩产生涌浪的计算方法
CN114282766B (zh) 一种无埋藏冰冰碛堤的冰碛湖溃决的易发性判断方法
Soltani et al. Influence of vegetated banks on turbulent characteristics of non-uniform flow in gravel bed river
CN110909293B (zh) 一种冰湖涌浪高度的计算方法
Meier Glaciers of the gannett peak-fremont peak area, wyoming
Billi et al. Field investigation on step-pool morphology and processes in steep mountain streams.
Abu-Salim Geomorphological Variables that Determine the Magnitude of Stream Floods in Wadi Al-Hasa-South Jordan
Kempema et al. Variations in anchor-ice crystal morphology related to river flow characteristics

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