一种用于后期冰水泥石流的预警方法及其应用与流程

1.本发明涉及到泥石流防治工程技术领域，尤其涉及一种用于后期冰水泥石流的预警方法及其应用。

背景技术：

2.冰水泥石流是一种发生在高山冰雪分布区的自然现象。冰水泥石流的发生往往是由于当地温度升高，冰雪融化形成沟道径流，起动沟道中的固体物源，形成冰水泥石流。
3.冰水泥石流的水源来自于冰雪融化：最早的融化全部来自于表面的积雪，可能形成积雪消融期的初期泥石流；随后可能形成积雪消融期的后期的冰水泥石流，这两类冰水泥石流的发生机理不同，因为提供激发泥石流的水源：积雪消融产生大量水源的机理不同；积雪完全消融后，积雪覆盖的冰川全部露出来，如果再形成冰水泥石流，其水源全部来源于冰川融化：可能形成冰川融化期的初期泥石流；随后可能形成冰川融化期的后期的冰水泥石流，这两类冰水泥石流的发生机理不同，因为提供激发泥石流的水源：冰川融化产生大量水源的机理不同。
4.目前，国内外对冰水泥石流形成之预警条件研究较少，其中地形条件的研究比较粗糙，主要通过少量案例统计得出大致范围，以及比较集中的影响范围的简单数据统计，没有深入研究其内在机理，如：西藏泥石流与环境，1999：29-40，其涉及到的地形条件包括：全流域面积，含堆积扇，7-90km2，其中15-32km2最有利冰水泥石流形成；现代冰川和积雪面积比，不含堆积扇，0.09-0.36，其中大于0.1最有利冰水泥石流形成；全流域相对高差：2500-3858m，其中大于3000m最有利冰水泥石流形成；冰雪水源区平均纵坡：0.197-0.815，其中大于0.5最有利冰水泥石流形成；全冰雪水源区以下主沟平均纵坡：0.08-0.361，其中0.1-0.22最有利冰水泥石流形成；其他如沟床纵剖面形状为下凹形、瓢形与栎叶的流域形状都有利冰水泥石流形成的定性描述。这些判断指标较多，范围较广，很难对冰水泥石流的危险性作出正确判断。
5.公开号为cn 106355275a，公开日为2017年01月25日的中国专利文献公开了一种冰川分布区泥石流灾害危险性分级评价方法，其特征在于：包括以下步骤：1)确定用于描述影响泥石流发育常规因子的指标，并分级；2)计算泥石流常规影响因子的加权和，得出无冰川影响下泥石流危险性值，并分级分区；3)分析冰川分布和活动对泥石流发育的影响，确定用于描述冰川影响的指标，并分级；4)确定冰川因子相对常规因子的相对重要性值；5)计算常规因子加权和和冰川影响因子乘以其相对重要性的和，并进行分级，完成冰川分布区泥石流危险性评价。
6.该专利文献公开的冰川分布区泥石流灾害危险性分级评价方法，没有涉及冰水泥石流预警，并未对导致冰水泥石流形成的条件及冰川融化进行内在机理的深入研究，影响整个冰水泥石流的预警准确性。

技术实现要素：

7.本发明为了克服上述现有技术的缺陷，提供一种用于后期冰水泥石流的预警方法及其应用，本发明针对导致冰水泥石流形成的条件及冰川融化进行内在机理的深入研究，建立冰川消融期后期冰水泥石流预警模型，提高了整个冰水泥石流的预警准确性。
8.本发明通过下述技术方案实现：
9.一种用于后期冰水泥石流的预警方法，其特征在于，包括以下步骤：
10.s1、获取冰水泥石流的基本地形数据；
11.s2、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年冰川消融期有效正积温ty，其中冰川消融期有效正积温ty的计算时间段为冰川消融期第一次泥石流发生的那一日为积温统计的起始日到再次发生泥石流时为止；
12.s3、通过式1计算冰水泥石流流域的地形因子r；
13.r＝c+0.65lg(a/w2)
ꢀꢀꢀꢀꢀꢀꢀꢀ
式1
14.其中，r为冰水泥石流地形因子；c为冰水泥石流坡向因子，由式2计算；a为冰川面积，m2；w为冰川区以下堆积区以上的形成区的沟床坡度为17度处的沟床平均宽度，m；
15.c＝sin(θ/2)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式2
16.其中，θ为冰水泥石流的冰川坡向，度，0
°
＜θ≤360
°
；
17.s4、通过式3计算诱发冰水泥石流的温度因子t：
18.t＝(t0+0.27t
7-0.48ty)/t
p
ꢀꢀꢀꢀ
式3
19.其中，t为诱发冰水泥石流的温度因子；t0为当日日均温度，度；t7为前7日积温和，度；ty为冰川消融期有效正积温，度；t
p
为当地多年平均积温，度；
20.s5、计算冰水泥石流的发生指标p，判断冰水泥石流的发生可能性；
21.p＝10r+t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式4
22.其中，p为冰水泥石流的发生指标；
23.当p＜29.4时，泥石流发生的可能性小；
24.当30.2＞p≥29.4时，泥石流发生的可能性中等；
25.当p≥30.2时，泥石流发生的可能性大。
26.所述步骤s1中，获取冰水泥石流的基本地形数据具体是指通过谷歌地球或现场测量确定冰川面积a、冰川区以下堆积区以上的形成区的沟床坡度为17度处的沟床平均宽度w和冰水泥石流的冰川坡向θ。
27.本发明适用于冰川及积雪区以下堆积区以上的形成区的沟床坡度有17度以上的沟道的冰水泥石流预警。
28.所述冰川及积雪区以下堆积区以上的形成区具体是指冰川流域中下游沟道。
29.本发明的基本原理如下：
30.冰水泥石流的形成机理是由于温度的升高，造成冰水泥石流流域的上游冰雪融化，起动中下游的松散固体物源，形成冰水泥石流。因此这类冰水泥石流的形成机理与我国西部地区普遍的沟床起动型泥石流不同，也与我国东部地区的浅层滑坡汇集型泥石流不同。沟床起动型泥石流是由于强降雨形成的大流量山洪侵蚀沟床松散固体物源，形成泥石流。浅层滑坡汇集型泥石流是由于较长时间的降雨和较强降雨，先有较多的浅层滑坡进入沟道形成泥石流物源，再被较大流量洪水搬运，形成泥石流。这两种类型的泥石流都需要较
强的降雨强度形成较大流量洪水，侵蚀搬运沟床固体物源。但是由于温度升高造成的上游冰雪融化形成的流域洪水往往持续时间较长，但流量较小，很难像西部一般的泥石流一样，直接由洪水侵蚀搬运沟道固体物源形成泥石流。在我国冰水泥石流的发生地点以西藏和新疆为主，这些地区的冰水泥石流流域的沟道坡度较大。冰水泥石流的形成机理为：上游冰雪融化形成的流域洪水在沟道内局部水位抬升，当水位抬升较高时，造成局部松散固体物源饱水失稳，并造成失稳区域上游侵蚀，下游冲刷，最后形成冰水泥石流。松散固体物源饱水失稳所需要的坡度远小于松散颗粒的休止角。因此沟道坡度过小，不能达到松散固体物源饱水失稳所需要的最小坡度，冰水泥石流很难发生。如果坡度过大，上游冰雪融化形成的洪水流动过快，水位难以抬升，也很难形成冰水泥石流。而形成这种松散固体物源饱水失稳的最小坡度就是17度。在冰川消融期中，由于短时间内冰川冰自身产生的融水量较少，难以形成起动泥石流所需要的较大流量，所以单纯依靠短时间内冰川融水来激发泥石流较为困难。当冰川表面接受太阳辐射产生后，产生的融水沿冰川内的垂直裂隙下渗至冰川内部，冰川内部存在尚未贯通的裂隙，所以上层产生的融水会储存在冰川内部的裂隙中。但是贯通的裂隙无法储存融水，当冰川内部储存的融水达到一定量时，冰川内部未打通的裂隙将会贯通，其中储存的融水将会在短时间内大量释放，为冰水泥石流的起动提供充足的水源。
31.在冰川消融期第一次泥石流发生时，由于冰川内部不存在已经贯通的裂隙，表层冰川产生的融水绝大部分都进入到未贯通的裂隙中；当温度上升到一定条件，储存在冰川未打通裂隙的融水突然释放，提供足够的水流，起动沟床固体物质形成冰水泥石流。随着冰川消融的持续进行，到了第二次及以后泥石流发生时，表层冰川产生的融水大部分进入能够顺利将融水排出到已打通裂隙中，只有少部分进入能够储存融水的未打通裂隙中，所以对于第二次及以后泥石流来说，还需要更高的气温来增加短时融水量，使得融水进入未打通的裂隙中，才有可能使总融水量再次达到能够激发泥石流的水平。所以冰川消融期第一次泥石流与冰川消融期第二次及以后泥石流有所不同。本发明针对消融过程后期的冰川消融过程，仅仅对冰川消融期后期泥石流进行预警。
32.本发明通过对比研究西藏和新疆公路沿线泥石流的发育与分布特征，研究得出以冰水泥石流流域中下游沟道坡度为依据的冰水泥石流沟道地形条件。从冰水泥石流的形成机理出发，考虑了地形因素中的坡度这个重要因素的作用，体现出影响冰水泥石流形成的决定因素及其重要性。同时，考虑不同区域的差别，考虑了当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年冰川消融期有效正积温ty这些因素，共同组成了诱发冰水泥石流的温度因子t；最终，结合冰水泥石流流域的地形因子r与诱发冰水泥石流的温度因子t，组成了冰水泥石流在冰川消融期的后期冰水泥石流预警模型。
33.本发明的有益效果主要表现在以下方面：
34.一、本发明，s1、获取冰水泥石流的基本地形数据；s2、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年冰川消融期有效正积温ty，其中冰川消融期有效正积温ty的计算时间段为冰川消融期第一次泥石流发生的那一日为积温统计的起始日到再次发生泥石流时为止；s3、通过式1计算冰水泥石流流域的地形因子r；s4、通过式3计算诱发冰水泥石流的温度因子t；s5、计算冰水泥石流的发生指标p，判断冰水泥石流的发生可能性，针对导致冰水泥石流形成的条件及冰川融化进行内在机理的深入研究，建立冰川消融期后期冰水泥石流预警模型，提高了整个冰水泥石流的预警准确性。
35.二、本发明，将冰水泥石流的地形条件与温度条件有机结合起来，充分考虑地形与温度在冰水泥石流形成中的作用，使得冰水泥石流的预警判断更加合理，更加准确。
36.三、本发明，从冰水泥石流的形成机理出发，考虑了地形因素中的坡度这个重要因素的作用，体现出影响冰水泥石流形成的决定因素及其重要性，利于提高预警准确性。
37.四、本发明，考虑不同区域的差别，考虑了当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年冰川消融期有效正积温ty这些因素，共同组成了诱发冰水泥石流的温度因子t，再结合冰水泥石流流域的地形因子r与诱发冰水泥石流的温度因子t，组成了冰水泥石流在冰川消融期的后期冰水泥石流预警模型，极大的提高了冰水泥石流的预警准确度。
38.五、本发明，充分考虑了冰川消融期后期泥石流的形成机理，后期泥石流形成的水源特点与初期泥石流形成的水源特点有重要差别，并通过冰川消融期有效正积温体现出来，极大的提高了冰水泥石流的预警准确度。
具体实施方式
39.实施例1
40.一种用于后期冰水泥石流的预警方法，包括以下步骤：
41.s1、获取冰水泥石流的基本地形数据；
42.s2、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年冰川消融期有效正积温ty，其中冰川消融期有效正积温ty的计算时间段为冰川消融期第一次泥石流发生的那一日为积温统计的起始日到再次发生泥石流时为止；
43.s3、通过式1计算冰水泥石流流域的地形因子r；
44.r＝c+0.65lg(a/w2)
ꢀꢀꢀꢀꢀꢀꢀꢀ
式1
45.其中，r为冰水泥石流地形因子；c为冰水泥石流坡向因子，由式2计算；a为冰川面积，m2；w为冰川区以下堆积区以上的形成区的沟床坡度为17度处的沟床平均宽度，m；
46.c＝sin(θ/2)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式2
47.其中，θ为冰水泥石流的冰川坡向，度，0
°
＜θ≤360
°
；
48.s4、通过式3计算诱发冰水泥石流的温度因子t：
49.t＝(t0+0.27t
7-0.48ty)/t
p
ꢀꢀꢀꢀ
式3
50.其中，t为诱发冰水泥石流的温度因子；t0为当日日均温度，度；t7为前7日积温和，度；ty为冰川消融期有效正积温，度；t
p
为当地多年平均积温，度；
51.s5、计算冰水泥石流的发生指标p，判断冰水泥石流的发生可能性；
52.p＝10r+t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式4
53.其中，p为冰水泥石流的发生指标；
54.当p＜29.4时，泥石流发生的可能性小；
55.当30.2＞p≥29.4时，泥石流发生的可能性中等；
56.当p≥30.2时，泥石流发生的可能性大。
57.本实施例为最基本的实施方式，针对导致冰水泥石流形成的条件及冰川融化进行内在机理的深入研究，建立冰川消融期后期冰水泥石流预警模型，提高了整个冰水泥石流的预警准确性。
58.实施例2
59.一种用于后期冰水泥石流的预警方法，包括以下步骤：
60.s1、获取冰水泥石流的基本地形数据；
61.s2、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年冰川消融期有效正积温ty，其中冰川消融期有效正积温ty的计算时间段为冰川消融期第一次泥石流发生的那一日为积温统计的起始日到再次发生泥石流时为止；
62.s3、通过式1计算冰水泥石流流域的地形因子r；
63.r＝c+0.65lg(a/w2)
ꢀꢀꢀꢀꢀꢀꢀꢀ
式1
64.其中，r为冰水泥石流地形因子；c为冰水泥石流坡向因子，由式2计算；a为冰川面积，m2；w为冰川区以下堆积区以上的形成区的沟床坡度为17度处的沟床平均宽度，m；
65.c＝sin(θ/2)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式2
66.其中，θ为冰水泥石流的冰川坡向，度，0
°
＜θ≤360
°
；
67.s4、通过式3计算诱发冰水泥石流的温度因子t：
68.t＝(t0+0.27t
7-0.48ty)/t
p
ꢀꢀꢀꢀ
式3
69.其中，t为诱发冰水泥石流的温度因子；t0为当日日均温度，度；t7为前7日积温和，度；ty为冰川消融期有效正积温，度；t
p
为当地多年平均积温，度；
70.s5、计算冰水泥石流的发生指标p，判断冰水泥石流的发生可能性；
71.p＝10r+t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式4
72.其中，p为冰水泥石流的发生指标；
73.当p＜29.4时，泥石流发生的可能性小；
74.当30.2＞p≥29.4时，泥石流发生的可能性中等；
75.当p≥30.2时，泥石流发生的可能性大。
76.所述步骤s1中，获取冰水泥石流的基本地形数据具体是指通过谷歌地球或现场测量确定冰川面积a、冰川区以下堆积区以上的形成区的沟床坡度为17度处的沟床平均宽度w和冰水泥石流的冰川坡向θ。
77.本实施例为一较佳实施方式，将冰水泥石流的地形条件与温度条件有机结合起来，充分考虑地形与温度在冰水泥石流形成中的作用，使得冰水泥石流的预警判断更加合理，更加准确。
78.实施例3
79.一种用于后期冰水泥石流的预警方法，包括以下步骤：
80.s1、获取冰水泥石流的基本地形数据；
81.s2、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年冰川消融期有效正积温ty，其中冰川消融期有效正积温ty的计算时间段为冰川消融期第一次泥石流发生的那一日为积温统计的起始日到再次发生泥石流时为止；
82.s3、通过式1计算冰水泥石流流域的地形因子r；
83.r＝c+0.65lg(a/w2)
ꢀꢀꢀꢀꢀꢀꢀꢀ
式1
84.其中，r为冰水泥石流地形因子；c为冰水泥石流坡向因子，由式2计算；a为冰川面积，m2；w为冰川区以下堆积区以上的形成区的沟床坡度为17度处的沟床平均宽度，m；
85.c＝sin(θ/2)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式2
86.其中，θ为冰水泥石流的冰川坡向，度，0
°
＜θ≤360
°
；
87.s4、通过式3计算诱发冰水泥石流的温度因子t：
88.t＝(t0+0.27t
7-0.48ty)/t
p
ꢀꢀꢀꢀ
式3
89.其中，t为诱发冰水泥石流的温度因子；t0为当日日均温度，度；t7为前7日积温和，度；ty为冰川消融期有效正积温，度；t
p
为当地多年平均积温，度；
90.s5、计算冰水泥石流的发生指标p，判断冰水泥石流的发生可能性；
91.p＝10r+t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式4
92.其中，p为冰水泥石流的发生指标；
93.当p＜29.4时，泥石流发生的可能性小；
94.当30.2＞p≥29.4时，泥石流发生的可能性中等；
95.当p≥30.2时，泥石流发生的可能性大。
96.所述步骤s1中，获取冰水泥石流的基本地形数据具体是指通过谷歌地球或现场测量确定冰川面积a、冰川区以下堆积区以上的形成区的沟床坡度为17度处的沟床平均宽度w和冰水泥石流的冰川坡向θ。
97.适用于冰川及积雪区以下堆积区以上的形成区的沟床坡度有17度以上的沟道的冰水泥石流预警。
98.本实施例为又一较佳实施方式，从冰水泥石流的形成机理出发，考虑了地形因素中的坡度这个重要因素的作用，体现出影响冰水泥石流形成的决定因素及其重要性，利于提高预警准确性。
99.实施例4
100.一种用于后期冰水泥石流的预警方法，包括以下步骤：
101.s1、获取冰水泥石流的基本地形数据；
102.s2、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年冰川消融期有效正积温ty，其中冰川消融期有效正积温ty的计算时间段为冰川消融期第一次泥石流发生的那一日为积温统计的起始日到再次发生泥石流时为止；
103.s3、通过式1计算冰水泥石流流域的地形因子r；
104.r＝c+0.65lg(a/w2)
ꢀꢀꢀꢀꢀꢀꢀꢀ
式1
105.其中，r为冰水泥石流地形因子；c为冰水泥石流坡向因子，由式2计算；a为冰川面积，m2；w为冰川区以下堆积区以上的形成区的沟床坡度为17度处的沟床平均宽度，m；
106.c＝sin(θ/2)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式2
107.其中，θ为冰水泥石流的冰川坡向，度，0
°
＜θ≤360
°
；
108.s4、通过式3计算诱发冰水泥石流的温度因子t：
109.t＝(t0+0.27t
7-0.48ty)/t
p
ꢀꢀꢀꢀ
式3
110.其中，t为诱发冰水泥石流的温度因子；t0为当日日均温度，度；t7为前7日积温和，度；ty为冰川消融期有效正积温，度；t
p
为当地多年平均积温，度；
111.s5、计算冰水泥石流的发生指标p，判断冰水泥石流的发生可能性；
112.p＝10r+t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式4
113.其中，p为冰水泥石流的发生指标；
114.当p＜29.4时，泥石流发生的可能性小；
115.当30.2＞p≥29.4时，泥石流发生的可能性中等；
116.当p≥30.2时，泥石流发生的可能性大。
117.所述步骤s1中，获取冰水泥石流的基本地形数据具体是指通过谷歌地球或现场测量确定冰川面积a、冰川区以下堆积区以上的形成区的沟床坡度为17度处的沟床平均宽度w和冰水泥石流的冰川坡向θ。
118.适用于冰川及积雪区以下堆积区以上的形成区的沟床坡度有17度以上的沟道的冰水泥石流预警。
119.所述冰川及积雪区以下堆积区以上的形成区具体是指冰川流域中下游沟道。
120.本实施例为最佳实施方式，考虑不同区域的差别，考虑了当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年冰川消融期有效正积温ty这些因素，共同组成了诱发冰水泥石流的温度因子t，再结合冰水泥石流流域的地形因子r与诱发冰水泥石流的温度因子t，组成了冰水泥石流在冰川消融期的后期冰水泥石流预警模型，极大的提高了冰水泥石流的预警准确度。
121.下面以独库公路为例对本发明进行说明：
122.独库公路北起准格尔盆地西南边缘的石油重地——克拉玛依市独山子区，南迄于天山山脉南麓阿克苏市库车县，由北向南横跨整个天山山脉，其2/3左右的里程穿越海拔2000m以上的寒冻山区，全线穿过4次高山区，途径天山5条较大的河流。根据调查，泥石流在天山公路全线均有分布，但主要集中在北段和南段。天山公路全线各类泥石流共有231处，包括坡面泥石流和沟谷泥石流，其中北段有53处，中段有3处，南段有178处，其中，对公路危胁较严重的泥石流有89处，对公路构成毁灭性威胁重大的泥石流共有12处，其中10处位于北段，主要分布于k629-k660之间。
123.在k629-k660段的泥石流沟中，以k636泥石流沟的冰水泥石流灾害最为严重，在20世纪70年代对g217公路沿线进行勘测时，将k636这条最大泥石流沟误认为是不活跃的泥石流沟，未作防灾工程处理。但是随着近年来全球气温的普遍升高，研究区内的冰雪层消融进一步加剧，在夏季高温时，流域内冰雪融水量开始增大，导致g217公路k636冰水泥石流沟爆发频繁，严重影响了公路畅通。从1984年至今，k636沟由于高温引发的冰川消融期初期冰水泥石流共发生3次。
124.以g217公路k636泥石流沟发生冰川消融期后期泥石流的日期，以及k636泥石流沟和周围泥石流中，未发生冰川消融期后期泥石流且气温较高的日期为例，对冰水泥石流进行预警判断。
125.表1为g217公路k636泥石流沟及周围泥石流沟道中，高温条件下冰水泥石流的40个案例的参数及相应的冰水泥石流的发生指标p，以及泥石流的实际发生情况。由于所选泥石流沟道距离较近，所以当地多年平均积温t
p
采用相同值，取14.13℃。
126.表1
127.[0128][0129]
通过表1中冰水泥石流的发生指标p值计算结果结合实际发生情况显示：
[0130]
判断泥石流发生可能性大的有1个，实际发生泥石流；判断泥石流发生可能性中等的有6个，实际2个发生泥石流，4个未发生泥石流；判断泥石流发生可能性小的有33个，实际都未发生泥石流。
[0131]
综上所述，应用本发明所述方法对冰川消融期后期冰水泥石流的预警判断较为准确。

技术特征：

1.一种用于后期冰水泥石流的预警方法，其特征在于，包括以下步骤：s1、获取冰水泥石流的基本地形数据；s2、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年冰川消融期有效正积温t
y
，其中冰川消融期有效正积温t
y
的计算时间段为冰川消融期第一次泥石流发生的那一日为积温统计的起始日到再次发生泥石流时为止；s3、通过式1计算冰水泥石流流域的地形因子r；r＝c+0.65lg(a/w2)
ꢀꢀꢀꢀ
式1其中，r为冰水泥石流地形因子；c为冰水泥石流坡向因子，由式2计算；a为冰川面积，m2；w为冰川区以下堆积区以上的形成区的沟床坡度为17度处的沟床平均宽度，m；c＝sin(θ/2)
ꢀꢀꢀꢀ
式2其中，θ为冰水泥石流的冰川坡向，度，0
°
＜θ≤360
°
；s4、通过式3计算诱发冰水泥石流的温度因子t：t＝(t0+0.27t
7-0.48t
y
)/t
p
ꢀꢀꢀꢀ
式3其中，t为诱发冰水泥石流的温度因子；t0为当日日均温度，度；t7为前7日积温和，度；t
y
为冰川消融期有效正积温，度；t
p
为当地多年平均积温，度；s5、计算冰水泥石流的发生指标p，判断冰水泥石流的发生可能性；p＝10r+t
ꢀꢀꢀꢀ
式4其中，p为冰水泥石流的发生指标；当p＜29.4时，泥石流发生的可能性小；当30.2＞p≥29.4时，泥石流发生的可能性中等；当p≥30.2时，泥石流发生的可能性大。2.根据权利要求1所述的一种用于后期冰水泥石流的预警方法，其特征在于：所述步骤s1中，获取冰水泥石流的基本地形数据具体是指通过谷歌地球或现场测量确定冰川面积a、冰川区以下堆积区以上的形成区的沟床坡度为17度处的沟床平均宽度w和冰水泥石流的冰川坡向θ。3.根据权利要求1所述的一种用于后期冰水泥石流的预警方法，其特征在于：适用于冰川及积雪区以下堆积区以上的形成区的沟床坡度有17度以上的沟道的冰水泥石流预警。4.根据权利要求3所述的一种用于后期冰水泥石流的预警方法，其特征在于：所述冰川及积雪区以下堆积区以上的形成区具体是指冰川流域中下游沟道。

技术总结

本发明公开了一种用于后期冰水泥石流的预警方法，属于泥石流防治工程技术领域，其特征在于，包括以下步骤：S1、获取冰水泥石流的基本地形数据；S2、获取当地多年平均积温T

技术研发人员：

余斌 马啸宇 王有林 李树武 李为乐 杨贤 曹钧恒 包健

受保护的技术使用者：

中国电建集团西北勘测设计研究院有限公司

技术研发日：

2022.06.14

技术公布日：

2022/9/19