CN108733948B - 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 - Google Patents
一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 Download PDFInfo
- Publication number
- CN108733948B CN108733948B CN201810531520.3A CN201810531520A CN108733948B CN 108733948 B CN108733948 B CN 108733948B CN 201810531520 A CN201810531520 A CN 201810531520A CN 108733948 B CN108733948 B CN 108733948B
- Authority
- CN
- China
- Prior art keywords
- ice
- lake
- moraine
- formula
- risk
- 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
Links
Images
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
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)
- Testing Or Calibration Of Command Recording Devices (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
本发明公开了一种无冰芯的冰碛堤冰湖溃决的危险性判断方法,其特征在于:包括以下步骤:a、判断冰碛堤内无冰芯的冰碛湖;b、通过谷歌地球确定无冰芯的冰碛湖的基本参数;c、计算无冰芯冰碛湖的危险性综合判别因子G;d、通过满足基本条件:α>8°、R>0.1和V>0.75,并且结合无冰芯冰碛湖的危险性综合判别因子G对无冰芯冰碛湖的危险性级别进行划分;G值越大,危险性越大;反之,G值越小,危险性越小。本发明通过分析各影响因子对冰湖溃决发生过程进行内部机理研究,建立了无冰芯冰碛湖溃决危险性计算模型,以定量的方式精确划分无冰芯冰碛湖溃决危险性级别,保障了判断准确性,进而利于提高防治效果。
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最有利于冰湖溃决。除了以上因素对冰湖溃决的影响外,冰川坡向因为日照和辐射的不同而影响冰崩或冰滑坡的产生,从而间接对冰湖溃决产生影响。
公开号为CN105513285A,公开日为2016年04月20日的中国专利文献公开了一种冰湖溃决预警方法,以潜在溃决冰湖为预警监测区域,测量确定潜在溃决冰湖的冰舌坡度α和冰川坡向θ,查看气象资料确定潜在冰湖溃决点前两年的平均正积温Ty2、多年平均正积温Ty,实时监测某个时间之前的正积温T0及该时间前7天平均温度和T7,以所述预警监测区域冰湖溃决临界值Cr为监测值,根据监测值大小划分冰湖溃决的预警级别。
公开号为CN105740616A,公开日为2016年07月06日的中国专利文献公开了一种冰湖溃决的危险性判断方法,其特征在于,包括以下步骤:通过谷歌地球确定溃决冰湖及周边未溃决冰湖的基本地形数据:冰湖后缘冰舌坡度α,冰川的坡向θ;计算冰湖溃决的地形综合判别因子S:由冰舌坡度、冰川坡向:S=tan(α)+0.5sin(θ/2),以冰川的地形综合判别因子S对冰湖溃决的危险性等级进行划分。以上述专利文献为代表的现有技术,采用了冰湖后缘冰舌坡度、冰川的坡向、潜在冰崩体质心至冰湖垂直高度、潜在冰崩体质心至冰湖运动距离、溃决口宽度、溃决口厚度、溃决前冰湖宽度、溃决前冰湖长宽比、溃决前冰湖岸边斜坡和连续性冰川的面积这些参数来识别危险性冰湖与预警冰湖的发生。但均没有全面考虑冰川与冰湖之间的关系以及冰碛堤的特征,特别是潜在冰崩体积与冰湖体积关系、冰湖面积与冰碛堤宽度和长度的比例关系、冰碛堤背水坡坡度这些关键参数,因此对冰碛湖的危险性判断有较大的欠缺。
发明内容
本发明为了克服上述现有技术的缺陷,提供一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用,本发明通过分析各影响因子对冰湖溃决发生过程进行内部机理研究,作出综合判断,建立了无冰芯冰碛湖溃决危险性计算模型,以定量的方式精确划分无冰芯冰碛湖溃决危险性级别,保障了判断准确性,进而利于提高防治效果。
本发明通过下述技术方案实现:
一种无冰芯的冰碛堤冰湖溃决的危险性判断方法,其特征在于,包括以下步骤:
a、判断冰碛堤内无冰芯的冰碛湖;
通过谷歌地球影像来区别,若表面圆滑并有分层的弯曲山脊冰碛堤,则为有冰芯的冰碛湖;若小的冰川前有一个不成比例的大冰碛堤,则为有冰芯的冰碛湖;若有狭窄的尖顶或横断面为三角形的冰碛堤,则为无冰芯的冰碛湖;
b、通过谷歌地球确定无冰芯的冰碛湖的基本参数;
包括冰湖后缘冰舌坡度α、母冰川坡向θ、母冰川面积F、冰舌面积 AS、冰舌质心至冰湖垂直高度差H、冰舌质心至冰湖运动距离L、冰碛湖面积A、冰碛湖溢流口宽度W、冰碛湖溢流口长度C、冰碛堤背水坡坡度β;
c、计算无冰芯冰碛湖的危险性综合判别因子G,依式1计算确定;
G =X+D+S+0.5R+0.17Log(MV) 式1
式中,G 为无冰芯冰碛湖的危险性综合判别因子;X为母冰川坡向因子,由式2计算;
X=sin(θ/2) 式2
其中,θ为母冰川坡向,单位°,适用于北半球,是指母冰川的平面法线在平面投影指向与正北方向夹角,0°<θ≤360°;
D为冰碛堤背水坡坡度因子,由式3计算;
D=tanβ 式3
其中,β为冰碛堤背水坡坡度,单位°;
S为冰舌坡度因子,由式4计算;
S=tanα 式4
其中,α为冰湖后缘冰舌坡度,单位°;
R为冰崩运动因子,由式5计算,若R>1,取R=1;
R=H/L 式5
其中,H为冰舌质心至冰湖垂直高度差,单位m;L为冰舌质心至冰湖运动距离,单位m;
M为冰碛湖因子,由式6计算;
M=A/(WC) 式6
其中,A为冰碛湖面积,单位m2;W为冰碛湖溢流口宽度,单位m;C 为冰碛湖溢流口长度,单位m;
V为体积比因子,由式7计算;
V=VG/VL 式7
其中,VG为冰舌体积,单位m3,由式8计算;VL为冰碛湖体积,单位m3,由式9计算;
VG=ASh 式8
VL=0.104A1.42 式9
式中,AS为冰舌面积,单位m2;h为冰舌厚度,计算获得,单位m;
d、通过满足基本条件:α>8°、R>0.1和V>0.75,并且结合无冰芯冰碛湖的危险性综合判别因子G对无冰芯冰碛湖的危险性级别进行划分;G值越大,危险性越大;反之,G值越小,危险性越小;当G>2.33时,危险性大;当1.76<G≤2.33时,危险性中等;当G≤1.76时,危险性小。
所述步骤a中,有冰芯的冰碛湖通过冰碛堤的渗流进行排水,无冰芯的冰碛湖通过冰碛堤表面的溢流口进行排水。
所述步骤c中,若为海洋性冰川,则冰舌厚度h通过式10计算;
h=5.2+15.4F0.5 式10
若为大陆性冰川,则冰舌厚度h通过式11计算;
h=-11.32+53.21F0.3 式11
若为悬冰川,则冰舌厚度h通过式12计算;
h=34.4F0.45 式12。
本发明,适用于高海拔或高纬度冰冻圈的冰湖溃决的危险性级别判断。
本发明的基本原理如下:
无冰芯的冰碛湖溃决,是因为补给其水源的冰川前部,特别是伸入冰湖内或距冰湖很近的冰舌部位运动速度瞬间发生了突变,由缓慢前进变为瞬间高速运动造成。冰舌坡度是影响冰崩和冰滑坡发生的最主要因素之一,冰舌坡度的陡缓不仅影响着冰体的运动形式,而且影响着冰舌段的汇流条件。随着气候的干湿冷暖波动,在气候转向湿冷或干热年代,随着年平均温度的不断上升,冰川的温度也慢慢提高,冻结粘滞力慢慢变小,冰川的塑性增强,运动速度逐渐加快。特别在冰川前部和冰舌地段,由于海拔较低,温度较高,冰川温度相应较上部要高,塑性更强,运动速度更快。当冰川的活动性提高到一定水平,盛夏或秋天的消融水流沿着冰裂隙、孔隙下渗,缓慢地暖切割冰体,润滑底床,不断地减小冰川向下前进运动的阻力。当这种阻力减小到一定程度时,冰川内部积累的应力瞬间释放,将冰川前部、特别是冰舌地段的冰体一齐挤推入湖。冰舌地段如果坡度很大,则以冰崩形式发生;冰舌地段如果坡度稍小,则以快速运动的冰滑坡形式发生。温度升高可以造成冰崩或冰滑坡,地震也可以造成冰崩或冰滑坡。
由于母冰川的坡向不同,日照和辐射条件各异,导致气温有明显的差异。在我国,一般地,位于北坡的温度最低,因此在地形上,坡向与正北方向夹角越大,越有利于冰川内部裂隙的产生,越有利于冰崩或冰滑坡的产生。冰崩运动因子代表了冰崩可能运动的距离,与冰舌质心与冰湖的高差和冰舌质心与冰湖的距离有关。冰崩运动因子R越大,代表冰崩体进入冰碛湖的可能性越大,造成冰碛湖溃决的可能性越大。冰碛堤的背水坡坡度代表冰碛湖被侵蚀、拉槽并溃决的难易程度。冰碛湖溃决机理是冰崩或冰滑坡进入冰湖,形成涌浪和冰湖水位抬升,侵蚀冰碛湖溢流口的表面粗大块石,因为溢流口表面坡度较小,而背水坡坡度较大,因此侵蚀从溢流口表面与背水坡结合部开始,沿着溢流口从下游往上游逐渐侵蚀,直到侵蚀到迎水坡,完成溃决。因此背水坡的坡度对冰碛湖的溃决影响很大。
以往的研究表明,较大的冰碛湖体积不会发生溃决,这是因为太大的冰碛湖容量太大,因为冰崩造成的水位上涨太小而不能造成足够的冰碛堤侵蚀形成冰湖溃决,因此有冰碛湖面积在100000-1000000m2范围内最有利于冰碛湖溃决的研究结果。但实际上在这个范围之外仍然有冰碛湖溃决的案例。事实上,冰碛湖的体积较大是一个相对值,即相对于一定的冰崩体积而言可能会太大,造成水位上涨太小;如果冰崩体积较大,即使较大的冰碛湖体积仍然可能有较大的水位上涨并形成足够的冰碛堤侵蚀,形成冰碛湖溃决。因此冰崩体体积与冰湖体积的比值才具有评价的意义。因为冰崩体体积很难确定,而冰崩体体积与冰舌体积相关,因此本发明采用冰舌体积与冰湖体积比因子V用于危险冰碛湖的危险性,可以很好地用于冰碛湖的危险性:冰舌体积与冰湖体积比因子V越大,越有利于冰碛湖溃决。如前所述,冰碛湖的面积太小,也不容易发生冰碛湖溃决,这是因为太小的冰碛湖面积,在溢流口有大量的溢流水流流出情况下,水位下降很快,还没有侵蚀到冰碛堤的迎水坡时,水位就下降到无法侵蚀溢流口粗大块石的程度,因此发生冰碛湖溃决也很困难。但该冰碛湖面积也是一个相对值,相对于溢流口溢流水的总量,即溢流流量与溢流时间的乘积:流量,与溢流口的宽度成正比;形成溃决需要侵蚀的时间,由冰碛堤的长度决定。因此冰崩运动因子R用于评估冰碛湖面积与溢流口宽度及长度的比例关系,可以很好地危险性危险冰碛湖:冰崩运动因子R越大,越有利于冰碛湖溃决。
本发明的有益效果主要表现在以下方面:
一、本发明,“a、判断冰碛堤内无冰芯的冰碛湖;b、通过谷歌地球确定无冰芯的冰碛湖的基本参数;c、计算无冰芯冰碛湖的危险性综合判别因子G,依式1计算确定;d、通过满足基本条件:α>8°、R>0.1和V>0.75,并且结合无冰芯冰碛湖的危险性综合判别因子G对无冰芯冰碛湖的危险性级别进行划分;当G>2.33时,危险性大;当1.76<G≤2.33时,危险性中等;当G≤1.76时,危险性小”,通过分析各影响因子对冰湖溃决发生过程进行内部机理研究,作出综合判断,建立了无冰芯冰碛湖溃决危险性计算模型,以定量的方式精确划分无冰芯冰碛湖溃决危险性级别,保障了判断准确性,进而利于提高防治效果。
二、本发明,区别无冰芯冰碛湖和有冰芯冰碛湖的溃决机理,首先判断出无冰芯冰碛湖,再进行冰碛湖危险性判断,使得危险性判断更有目的性,更加准确;危险性判断全面考虑冰川条件,母冰川坡向因子X、冰舌坡度因子S;冰川与冰湖之间的条件,体积比因子V、冰崩运动因子R;冰湖条件,冰碛堤背水坡坡度因子D、冰碛湖因子M的影响因素,能够以定量的方式精确划分无冰芯冰碛湖溃决危险性级别。
三、本发明,对于冰湖后缘冰舌坡度α设定下限值8°,在无冰芯冰碛湖的危险性综合判别因子G中同时考虑到了冰湖后缘冰舌坡度α的作用,使危险性判断更简单,更合理。
四、本发明,对于冰崩运动因子R设定下限值0.1,在无冰芯冰碛湖的危险性综合判别因子G中同时考虑到了冰崩运动因子R的作用,使危险性判断更简洁,更合理;对于冰崩运动因子R设定上限值,避免了在冰舌处于冰碛湖上方,冰舌质心至冰湖运动距离L为0,冰崩运动因子R无穷大而无法判断的情况,使危险性判断更合理,更可靠。
五、本发明,综合考虑冰湖体积大小的影响机理与原因,将冰舌与冰湖体积比因子用于判断冰碛湖危险性,并设定其下限值0.75,避免了单一冰碛湖体积判断的片面性,使得危险性判断更加严谨和准确。
六、本发明,综合考虑冰湖体积大小的影响机理与原因,将冰碛湖面积与冰碛湖溢流口宽度和长度比因子用于判断冰碛湖危险性,避免了单一冰碛湖面积判断的片面性,使得危险性判断更加合理可靠。
七、本发明,考虑了对冰碛堤侵蚀有很大影响的冰碛堤背水坡坡度β,将主要的冰碛堤溃决影响因素完全考虑在内,使得危险性的判断更全面,更可靠;计算无冰芯冰碛湖的危险性综合判别因子G所涉及的母冰川坡向因子X、冰碛堤背水坡坡度因子D、冰舌坡度因子S、冰崩运动因子R、冰碛湖因子M和体积比因子V均为无量纲参数,在各种无冰芯冰碛湖溃决条件下都可以使用,适用性强。
八、本发明,适用于高海拔或高纬度冰冻圈的冰湖溃决的危险性级别判断,能够为这些地区的防灾减灾提供依据,为这些地区的发展和规划提供参考,还能够为冰湖溃决的监测预警提供基本数据。
附图说明
图1是无冰芯的冰碛湖的基本参数的获取示意图;
图2是图1中A-A示意图;
图3是图1中B-B示意图;
其中:1、母冰川面积F,2、母冰川坡向θ,3、冰舌面积AS,4、冰湖后缘冰舌坡度α,5、冰舌质心至冰湖垂直高度差H,6、冰舌质心至冰湖运动距离L,7、冰碛湖面积A,8、冰碛堤背水坡坡度β,9、冰碛湖溢流口长度C,10、冰碛湖溢流口宽度W。
具体实施方式
实施例1
一种无冰芯的冰碛堤冰湖溃决的危险性判断方法,包括以下步骤:
a、判断冰碛堤内无冰芯的冰碛湖;
通过谷歌地球影像来区别,若表面圆滑并有分层的弯曲山脊冰碛堤,则为有冰芯的冰碛湖;若小的冰川前有一个不成比例的大冰碛堤,则为有冰芯的冰碛湖;若有狭窄的尖顶或横断面为三角形的冰碛堤,则为无冰芯的冰碛湖;
b、通过谷歌地球确定无冰芯的冰碛湖的基本参数;
包括冰湖后缘冰舌坡度α、母冰川坡向θ、母冰川面积F、冰舌面积 AS、冰舌质心至冰湖垂直高度差H、冰舌质心至冰湖运动距离L、冰碛湖面积A、冰碛湖溢流口宽度W、冰碛湖溢流口长度C、冰碛堤背水坡坡度β;
c、计算无冰芯冰碛湖的危险性综合判别因子G,依式1计算确定;
G =X+D+S+0.5R+0.17Log(MV) 式1
式中,G 为无冰芯冰碛湖的危险性综合判别因子;X为母冰川坡向因子,由式2计算;
X=sin(θ/2) 式2
其中,θ为母冰川坡向,单位°,适用于北半球,是指母冰川的平面法线在平面投影指向与正北方向夹角,0°<θ≤360°;
D为冰碛堤背水坡坡度因子,由式3计算;
D=tanβ 式3
其中,β为冰碛堤背水坡坡度,单位°;
S为冰舌坡度因子,由式4计算;
S=tanα 式4
其中,α为冰湖后缘冰舌坡度,单位°;
R为冰崩运动因子,由式5计算,若R>1,取R=1;
R=H/L 式5
其中,H为冰舌质心至冰湖垂直高度差,单位m;L为冰舌质心至冰湖运动距离,单位m;
M为冰碛湖因子,由式6计算;
M=A/(WC) 式6
其中,A为冰碛湖面积,单位m2;W为冰碛湖溢流口宽度,单位m;C 为冰碛湖溢流口长度,单位m;
V为体积比因子,由式7计算;
V=VG/VL 式7
其中,VG为冰舌体积,单位m3,由式8计算;VL为冰碛湖体积,单位m3,由式9计算;
VG=ASh 式8
VL=0.104A1.42 式9
式中,AS为冰舌面积,单位m2;h为冰舌厚度,计算获得,单位m;
d、通过满足基本条件:α>8°、R>0.1和V>0.75,并且结合无冰芯冰碛湖的危险性综合判别因子G对无冰芯冰碛湖的危险性级别进行划分;G值越大,危险性越大;反之,G值越小,危险性越小;当G>2.33时,危险性大;当1.76<G≤2.33时,危险性中等;当G≤1.76时,危险性小。
“a、判断冰碛堤内无冰芯的冰碛湖;b、通过谷歌地球确定无冰芯的冰碛湖的基本参数;c、计算无冰芯冰碛湖的危险性综合判别因子G,依式 1计算确定;d、通过满足基本条件:α>8°、R>0.1和V>0.75,并且结合无冰芯冰碛湖的危险性综合判别因子G对无冰芯冰碛湖的危险性级别进行划分;当G>2.33时,危险性大;当1.76<G≤2.33时,危险性中等;当G≤1.76时,危险性小”,通过分析各影响因子对冰湖溃决发生过程进行内部机理研究,作出综合判断,建立了无冰芯冰碛湖溃决危险性计算模型,以定量的方式精确划分无冰芯冰碛湖溃决危险性级别,保障了判断准确性,进而利于提高防治效果。
实施例2
一种无冰芯的冰碛堤冰湖溃决的危险性判断方法,包括以下步骤:
a、判断冰碛堤内无冰芯的冰碛湖;
通过谷歌地球影像来区别,若表面圆滑并有分层的弯曲山脊冰碛堤,则为有冰芯的冰碛湖;若小的冰川前有一个不成比例的大冰碛堤,则为有冰芯的冰碛湖;若有狭窄的尖顶或横断面为三角形的冰碛堤,则为无冰芯的冰碛湖;
b、通过谷歌地球确定无冰芯的冰碛湖的基本参数;
包括冰湖后缘冰舌坡度α、母冰川坡向θ、母冰川面积F、冰舌面积 AS、冰舌质心至冰湖垂直高度差H、冰舌质心至冰湖运动距离L、冰碛湖面积A、冰碛湖溢流口宽度W、冰碛湖溢流口长度C、冰碛堤背水坡坡度β;
c、计算无冰芯冰碛湖的危险性综合判别因子G,依式1计算确定;
G =X+D+S+0.5R+0.17Log(MV) 式1
式中,G 为无冰芯冰碛湖的危险性综合判别因子;X为母冰川坡向因子,由式2计算;
X=sin(θ/2) 式2
其中,θ为母冰川坡向,单位°,适用于北半球,是指母冰川的平面法线在平面投影指向与正北方向夹角,0°<θ≤360°;
D为冰碛堤背水坡坡度因子,由式3计算;
D=tanβ 式3
其中,β为冰碛堤背水坡坡度,单位°;
S为冰舌坡度因子,由式4计算;
S=tanα 式4
其中,α为冰湖后缘冰舌坡度,单位°;
R为冰崩运动因子,由式5计算,若R>1,取R=1;
R=H/L 式5
其中,H为冰舌质心至冰湖垂直高度差,单位m;L为冰舌质心至冰湖运动距离,单位m;
M为冰碛湖因子,由式6计算;
M=A/(WC) 式6
其中,A为冰碛湖面积,单位m2;W为冰碛湖溢流口宽度,单位m;C 为冰碛湖溢流口长度,单位m;
V为体积比因子,由式7计算;
V=VG/VL 式7
其中,VG为冰舌体积,单位m3,由式8计算;VL为冰碛湖体积,单位m3,由式9计算;
VG=ASh 式8
VL=0.104A1.42 式9
式中,AS为冰舌面积,单位m2;h为冰舌厚度,计算获得,单位m;
d、通过满足基本条件:α>8°、R>0.1和V>0.75,并且结合无冰芯冰碛湖的危险性综合判别因子G对无冰芯冰碛湖的危险性级别进行划分;G值越大,危险性越大;反之,G值越小,危险性越小;当G>2.33时,危险性大;当1.76<G≤2.33时,危险性中等;当G≤1.76时,危险性小。
所述步骤a中,有冰芯的冰碛湖通过冰碛堤的渗流进行排水,无冰芯的冰碛湖通过冰碛堤表面的溢流口进行排水。
区别无冰芯冰碛湖和有冰芯冰碛湖的溃决机理,首先判断出无冰芯冰碛湖,再进行冰碛湖危险性判断,使得危险性判断更有目的性,更加准确;危险性判断全面考虑冰川条件,母冰川坡向因子X、冰舌坡度因子S;冰川与冰湖之间的条件,体积比因子V、冰崩运动因子R;冰湖条件,冰碛堤背水坡坡度因子D、冰碛湖因子M的影响因素,能够以定量的方式精确划分无冰芯冰碛湖溃决危险性级别。
对于冰湖后缘冰舌坡度α设定下限值8°,在无冰芯冰碛湖的危险性综合判别因子G中同时考虑到了冰湖后缘冰舌坡度α的作用,使危险性判断更简单,更合理。
对于冰崩运动因子R设定下限值0.1,在无冰芯冰碛湖的危险性综合判别因子G中同时考虑到了冰崩运动因子R的作用,使危险性判断更简洁,更合理;对于冰崩运动因子R设定上限值,避免了在冰舌处于冰碛湖上方,冰舌质心至冰湖运动距离L为0,冰崩运动因子R无穷大而无法判断的情况,使危险性判断更合理,更可靠。
实施例3
一种无冰芯的冰碛堤冰湖溃决的危险性判断方法,包括以下步骤:
a、判断冰碛堤内无冰芯的冰碛湖;
通过谷歌地球影像来区别,若表面圆滑并有分层的弯曲山脊冰碛堤,则为有冰芯的冰碛湖;若小的冰川前有一个不成比例的大冰碛堤,则为有冰芯的冰碛湖;若有狭窄的尖顶或横断面为三角形的冰碛堤,则为无冰芯的冰碛湖;
b、通过谷歌地球确定无冰芯的冰碛湖的基本参数;
包括冰湖后缘冰舌坡度α、母冰川坡向θ、母冰川面积F、冰舌面积 AS、冰舌质心至冰湖垂直高度差H、冰舌质心至冰湖运动距离L、冰碛湖面积A、冰碛湖溢流口宽度W、冰碛湖溢流口长度C、冰碛堤背水坡坡度β;
c、计算无冰芯冰碛湖的危险性综合判别因子G,依式1计算确定;
G =X+D+S+0.5R+0.17Log(MV) 式1
式中,G 为无冰芯冰碛湖的危险性综合判别因子;X为母冰川坡向因子,由式2计算;
X=sin(θ/2) 式2
其中,θ为母冰川坡向,单位°,适用于北半球,是指母冰川的平面法线在平面投影指向与正北方向夹角,0°<θ≤360°;
D为冰碛堤背水坡坡度因子,由式3计算;
D=tanβ 式3
其中,β为冰碛堤背水坡坡度,单位°;
S为冰舌坡度因子,由式4计算;
S=tanα 式4
其中,α为冰湖后缘冰舌坡度,单位°;
R为冰崩运动因子,由式5计算,若R>1,取R=1;
R=H/L 式5
其中,H为冰舌质心至冰湖垂直高度差,单位m;L为冰舌质心至冰湖运动距离,单位m;
M为冰碛湖因子,由式6计算;
M=A/(WC) 式6
其中,A为冰碛湖面积,单位m2;W为冰碛湖溢流口宽度,单位m;C 为冰碛湖溢流口长度,单位m;
V为体积比因子,由式7计算;
V=VG/VL 式7
其中,VG为冰舌体积,单位m3,由式8计算;VL为冰碛湖体积,单位m3,由式9计算;
VG=ASh 式8
VL=0.104A1.42 式9
式中,AS为冰舌面积,单位m2;h为冰舌厚度,计算获得,单位m;
d、通过满足基本条件:α>8°、R>0.1和V>0.75,并且结合无冰芯冰碛湖的危险性综合判别因子G对无冰芯冰碛湖的危险性级别进行划分;G值越大,危险性越大;反之,G值越小,危险性越小;当G>2.33时,危险性大;当1.76<G≤2.33时,危险性中等;当G≤1.76时,危险性小。
所述步骤a中,有冰芯的冰碛湖通过冰碛堤的渗流进行排水,无冰芯的冰碛湖通过冰碛堤表面的溢流口进行排水。
所述步骤c中,若为海洋性冰川,则冰舌厚度h通过式10计算;
h=5.2+15.4F0.5 式10
若为大陆性冰川,则冰舌厚度h通过式11计算;
h=-11.32+53.21F0.3 式11
若为悬冰川,则冰舌厚度h通过式12计算;
h=34.4F0.45 式12。
适用于高海拔或高纬度冰冻圈的冰湖溃决的危险性级别判断。
综合考虑冰湖体积大小的影响机理与原因,将冰舌与冰湖体积比因子用于判断冰碛湖危险性,并设定其下限值0.75,避免了单一冰碛湖体积判断的片面性,使得危险性判断更加严谨和准确。
综合考虑冰湖体积大小的影响机理与原因,将冰碛湖面积与冰碛湖溢流口宽度和长度比因子用于判断冰碛湖危险性,避免了单一冰碛湖面积判断的片面性,使得危险性判断更加合理可靠。
考虑了对冰碛堤侵蚀有很大影响的冰碛堤背水坡坡度β,将主要的冰碛堤溃决影响因素完全考虑在内,使得危险性的判断更全面,更可靠;计算无冰芯冰碛湖的危险性综合判别因子G所涉及的母冰川坡向因子X、冰碛堤背水坡坡度因子D、冰舌坡度因子S、冰崩运动因子R、冰碛湖因子M和体积比因子V均为无量纲参数,在各种无冰芯冰碛湖溃决条件下都可以使用,适用性强。
适用于高海拔或高纬度冰冻圈的冰湖溃决的危险性级别判断,能够为这些地区的防灾减灾提供依据,为这些地区的发展和规划提供参考,还能够为冰湖溃决的监测预警提供基本数据。
下面采用本发明方法对加拿大不列颠哥伦比亚省的无冰芯冰碛湖进行分析:
加拿大不列颠哥伦比亚省广泛分布有冰碛湖,其中有相当部分为无冰芯冰碛湖。在该省的海岸山区区域共发现有20处无冰芯冰碛湖,采用本发明的公式对无冰芯冰碛湖溃决危险性进行计算划分,本发明所用危险性判断计算方法选用的无冰芯冰碛湖的基本参数包括:冰湖后缘冰舌坡度α、母冰川坡向θ、母冰川面积F、冰舌面积AS、冰舌质心至冰湖垂直高度差H、冰舌质心至冰湖运动距离L、冰碛湖面积A、冰碛湖溢流口宽度W、冰碛湖溢流口长度C、冰碛堤背水坡坡度β。
根据式1-12,计算出危险性判断的参数母冰川坡向因子X、冰舌坡度因子S;冰川与冰湖之间的条件,体积比因子V、冰崩运动因子R;冰湖条件,冰碛堤背水坡坡度因子D、冰碛湖因子M,并计算出无冰芯冰碛湖的危险性综合判别因子G,20处无冰芯的冰碛湖参数计算情况、危险性分级以及无冰芯的冰碛湖溃决的实际发生情况如表1所示。
表1
通过满足基本条件:α>8°、R>0.1和V>0.75,并且结合无冰芯冰碛湖的危险性综合判别因子G对无冰芯冰碛湖的危险性级别进行划分;G值越大,危险性越大;反之,G值越小,危险性越小;当G>2.33时,危险性大;当1.76<G≤2.33时,危险性中等;当G≤1.76时,危险性小。
表1中G值计算结果显示:23处潜在无冰芯的冰碛湖冰湖中,危险性大的冰湖3处,危险性中等的冰湖2处,危险性小的冰湖18处。对比实际发生情况,3个危险性大的冰湖(见表1:序号1、2、3)都发生了冰湖溃决;2个危险性中的冰湖,都没有发生冰湖溃决(见表1:序号5、11); 18个危险性小的冰湖(见表1:序号3、4、6-10、12-23)都没有发生冰湖溃决。综上所述,应用本发明所述方法对冰湖溃决危险性划分的准确性较高。
Claims (4)
1.一种无冰芯的冰碛堤冰湖溃决的危险性判断方法,其特征在于,包括以下步骤:
a、判断冰碛堤内无冰芯的冰碛湖;
通过谷歌地球影像来区别,若表面圆滑并有分层的弯曲山脊冰碛堤,则为有冰芯的冰碛湖;若小的冰川前有一个不成比例的大冰碛堤,则为有冰芯的冰碛湖;若有狭窄的尖顶或横断面为三角形的冰碛堤,则为无冰芯的冰碛湖;
b、通过谷歌地球确定无冰芯的冰碛湖的基本参数;
包括冰湖后缘冰舌坡度α、母冰川坡向θ、母冰川面积F、冰舌面积AS、冰舌质心至冰湖垂直高度差H、冰舌质心至冰湖运动距离L、冰碛湖面积A、冰碛湖溢流口宽度W、冰碛湖溢流口长度C、冰碛堤背水坡坡度β;
c、计算无冰芯冰碛湖的危险性综合判别因子G,依式1计算确定;
G =X+D+S+0.5R+0.17Log(MV) 式1
式中,G 为无冰芯冰碛湖的危险性综合判别因子;X为母冰川坡向因子,由式2计算;
X=sin(θ/2) 式2
其中,θ为母冰川坡向,单位°,适用于北半球,是指母冰川的平面法线在平面投影指向与正北方向夹角,0°<θ≤360°;
D为冰碛堤背水坡坡度因子,由式3计算;
D=tanβ 式3
其中,β为冰碛堤背水坡坡度,单位°;
S为冰舌坡度因子,由式4计算;
S=tanα 式4
其中,α为冰湖后缘冰舌坡度,单位°;
R为冰崩运动因子,由式5计算,若R>1,取R=1;
R=H/L 式5
其中,H为冰舌质心至冰湖垂直高度差,单位m;L为冰舌质心至冰湖运动距离,单位m;
M为冰碛湖因子,由式6计算;
M=A/(WC) 式6
其中,A为冰碛湖面积,单位m2;W为冰碛湖溢流口宽度,单位m;C为冰碛湖溢流口长度,单位m;
V为体积比因子,由式7计算;
V=VG/VL 式7
其中,VG为冰舌体积,单位m3,由式8计算;VL为冰碛湖体积,单位m3,由式9计算;
VG=ASh 式8
VL=0.104A1.42 式9
式中,AS为冰舌面积,单位m2;h为冰舌厚度,计算获得,单位m;
d、通过满足基本条件:α>8°、R>0.1和V>0.75,并且结合无冰芯冰碛湖的危险性综合判别因子G对无冰芯冰碛湖的危险性级别进行划分;G值越大,危险性越大;反之,G值越小,危险性越小;当G>2.33时,危险性大;当1.76<G≤2.33时,危险性中等;当G≤1.76时,危险性小。
2.根据权利要求1所述的一种无冰芯的冰碛堤冰湖溃决的危险性判断方法,其特征在于:所述步骤a中,有冰芯的冰碛湖通过冰碛堤的渗流进行排水,无冰芯的冰碛湖通过冰碛堤表面的溢流口进行排水。
3.根据权利要求1所述的一种无冰芯的冰碛堤冰湖溃决的危险性判断方法,其特征在于:所述步骤c中,若为海洋性冰川,则冰舌厚度h通过式10计算;
h=5.2+15.4F0.5 式10
若为大陆性冰川,则冰舌厚度h通过式11计算;
h=-11.32+53.21F0.3 式11
若为悬冰川,则冰舌厚度h通过式12计算;
h=34.4F0.45 式12。
4.根据权利要求1所述的一种无冰芯的冰碛堤冰湖溃决的危险性判断方法的应用,其特征在于:适用于高海拔或高纬度冰冻圈的冰湖溃决的危险性级别判断。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810531520.3A CN108733948B (zh) | 2018-05-29 | 2018-05-29 | 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810531520.3A CN108733948B (zh) | 2018-05-29 | 2018-05-29 | 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108733948A CN108733948A (zh) | 2018-11-02 |
CN108733948B true CN108733948B (zh) | 2021-11-19 |
Family
ID=63936596
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810531520.3A Active CN108733948B (zh) | 2018-05-29 | 2018-05-29 | 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108733948B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109697314B (zh) * | 2018-12-18 | 2022-11-01 | 国网西藏电力有限公司 | 一种冰崩的最大运动距离的计算方法及其应用 |
CN110210173B (zh) * | 2019-06-20 | 2021-09-10 | 中国地质科学院探矿工艺研究所 | 一种冰碛土物理力学参数的快速估算方法 |
CN110820668A (zh) * | 2019-10-25 | 2020-02-21 | 成都理工大学 | 一种冰崩坠落在冰湖中造成的涌浪高度的计算方法 |
CN110909293B (zh) * | 2019-10-25 | 2022-04-26 | 成都理工大学 | 一种冰湖涌浪高度的计算方法 |
CN112735098B (zh) * | 2020-12-31 | 2022-08-09 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种震后冰湖溃决预警方法、装置及存储介质 |
CN114282766B (zh) * | 2021-11-26 | 2023-04-07 | 成都理工大学 | 一种无埋藏冰冰碛堤的冰碛湖溃决的易发性判断方法 |
CN118279752B (zh) * | 2024-05-31 | 2024-08-02 | 中国地质科学院探矿工艺研究所 | 一种冰崩型泥石流灾害链隐患识别方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101547069B1 (ko) * | 2014-12-05 | 2015-08-25 | 연세대학교 산학협력단 | 지반의 침식과 연행 작용을 고려한 토석류 거동 해석 방법 및 장치 |
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 | 成都理工大学 | 一种冰湖溃决的早期识别方法及其应用 |
-
2018
- 2018-05-29 CN CN201810531520.3A patent/CN108733948B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101547069B1 (ko) * | 2014-12-05 | 2015-08-25 | 연세대학교 산학협력단 | 지반의 침식과 연행 작용을 고려한 토석류 거동 해석 방법 및 장치 |
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 | 成都理工大学 | 一种冰湖溃决的早期识别方法及其应用 |
Non-Patent Citations (2)
Title |
---|
Preliminary study of the effect of earthquakes on the rainfall threshold of debris flows;B Yu,et al.;《Engineering Geology》;20141231;第130-135页 * |
藏东南典型冰湖溃决机制及危险性研究;常鸣等;《南水北调与水利科技》;20171231;第15卷(第6期);第115-121页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108733948A (zh) | 2018-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108733948B (zh) | 一种无冰芯的冰碛堤冰湖溃决的危险性判断方法及应用 | |
CN108681651B (zh) | 一种无冰芯的冰碛堤冰湖溃决的早期识别方法 | |
Costa | Physical geomorphology of debris flows | |
Evans | The geomorphology and morphometry of glacial and nival areas | |
Pierson et al. | Flow and deposition of posteruption hot lahars on the east side of Mount Pinatubo, July–October 1991 | |
Duman | The largest landslide dam in Turkey: Tortum landslide | |
Bovis et al. | The July 29, 1998, debris flow and landslide dam at Capricorn Creek, Mount Meager Volcanic Complex, southern Coast Mountains, British Columbia | |
Ge et al. | Case history of the disastrous debris flows of Tianmo Watershed in Bomi County, Tibet, China: Some mitigation suggestions | |
Decaulne et al. | Debris flow triggered by rapid snowmelt: a case study in the glei. arhjalli area, northwestern iceland | |
CN108090670B (zh) | 一种基于水流侵蚀能力的堰塞坝形成评估方法 | |
CN112735098B (zh) | 一种震后冰湖溃决预警方法、装置及存储介质 | |
Chigira et al. | Landslides in weathered granitic rocks in Japan and Malaysia | |
Turcotte et al. | Ice processes in a steep river basin | |
Gaillot et al. | Impact of gravel-mining on stream channel and coastal sediment supply: example of the Calvi Bay in Corsica (France) | |
Addison | Debris flow during intense rainfall in Snowdonia, North Wales: a preliminary survey | |
CN112836280A (zh) | 泥石流堰塞坝溃口下切速率计算方法 | |
PIERSON | TRANSFORMATION OF WATER FLOOD TO | |
Dotseva et al. | Assessment of debris flows-prone watersheds in southern slopes of Stara Planina Mountain by combined raster and morphometric analysis | |
Waskito et al. | Simulation of Saguling dam break using the HEC-RAS software | |
Daoed et al. | Backwater flow analysis due to the construction of the check dam series on Aie Dingin River using HEC-RAS | |
Ma et al. | Research on the debris flow hazards after the Wenchuan earthquake in Bayi Gully, Longchi, Dujiangyan, Sichuan Province, China | |
CN114282766B (zh) | 一种无埋藏冰冰碛堤的冰碛湖溃决的易发性判断方法 | |
Song et al. | Insights into the long-term stability of landslide dams on the eastern margin of the Tibetan Plateau, China–A case study of the Diexi area | |
Meyer et al. | Lahars from the 1992 eruptions of crater peak, Mount Spurr volcano, Alaska | |
Romanescu et al. | Ice dams and backwaters as hydrological risk phenomena–case study: the Bistrita River, upstream of the Izvorul Muntelui Lake (Romania) |
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 |