一种积雪消融期冰水泥石流预警方法与流程

1.本发明涉及到泥石流防治工程技术领域，尤其涉及一种积雪消融期冰水泥石流预警方法。


背景技术：

2.冰水泥石流的发生往往是由于当地温度升高，冰雪融化形成沟道径流，起动沟道中的固体物源，形成冰水泥石流。形成冰水泥石流需要具备两个条件：一是有利于发生冰水泥石流的地形条件；二是一定的温度条件导致冰雪融化并最终诱发冰水泥石流的发生。这些条件综合影响并决定冰水泥石流发生的可能性。其中，地形条件对冰水泥石流的影响主要由沟道坡度、积雪或冰川面积决定；温度是冰水泥石流的激发条件；地质条件也有影响，如堆积在沟道内的固体物质由区域的岩性决定；因此，冰水泥石流的预警，主要通过地形、地质与温度条件判断来实现。
3.冰水泥石流的水源来自于冰雪融化：最早的融化全部来自于表面的积雪，可能形成积雪消融期的第一次冰水泥石流；随后可能形成积雪消融期的第二次及以后的冰水泥石流，这两类冰水泥石流的发生机理不同，因为提供激发冰水泥石流的积雪消融产生的水源机理不同。
4.目前，国内外对冰水泥石流形成之预警条件研究较少，其中地形条件的研究比较粗糙，主要通过少量案例统计得出大致范围，以及比较集中的影响范围的简单数据统计，没有深入研究其内在机理，很难对冰水泥石流的危险性作出准确判断。
5.公开号为cn106355275a，公开日为2017年01月25日的中国专利文献公开了一种冰川分布区泥石流灾害危险性分级评价方法，其特征在于：包括以下步骤：
6.1)确定用于描述影响泥石流发育常规因子的指标，并分级；
7.2)计算泥石流常规影响因子的加权和，得出无冰川影响下泥石流危险性值，并分级分区；
8.3)分析冰川分布和活动对泥石流发育的影响，确定用于描述冰川影响的指标，并分级；
9.4)确定冰川因子相对常规因子的相对重要性值；
10.5)计算常规因子加权和和冰川影响因子乘以其相对重要性的和，并进行分级，完成冰川分布区泥石流危险性评价。
11.公开号为cn109472445a，公开日为2019年03月15日的中国专利文献公开了冰水泥石流的危险性判断方法及其应用，其特征在于，包括以下步骤：a、确定冰川及积雪区面积a0、冰川及积雪区以下堆积区以上的形成区的沟床坡度α和沟床宽度w；b、对不同坡度段的沟床坡度α和沟床宽度w分别进行测量，以沟道长度≥10m为标准测量沟道坡度；c、由冰川及积雪区相对面积a＝a0/w2和沟床坡度α对不同坡度段的冰水泥石流的危险性等级从高到低进行划分；d、以危险性等级最高的判断结果为准。
12.以上述专利文献为代表的现有技术，由于并未对导致冰水泥石流形成的地形、地
质和温度条件，以及积雪融化进行内在机理的深入研究，难于对冰水泥石流的预警作出准确判断。


技术实现要素：

13.本发明为了克服上述现有技术的缺陷，提供一种积雪消融期冰水泥石流预警方法，本发明针对导致冰水泥石流形成的地形条件及温度条件进行内在机理的深入研究，并作出综合判断，建立积雪消融期冰水泥石流预警模型，确定冰水泥石流的预警划分级别，能够准确的对冰水泥石流的发生作出判断。
14.本发明通过下述技术方案实现：
15.一种积雪消融期冰水泥石流预警方法，其特征在于，包括以下步骤：
16.a、通过现场测量确定冰水泥石流的基本地形数据，包括积雪区面积a、冰水泥石流形成区纵坡降为17度处的沟床平均宽度w和积雪区坡向θ；
17.b、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年积雪消融期有效正积温ty；
18.c、通过式1计算冰水泥石流流域地形因子r；
19.r＝x+0.65lg(a/w2)
ꢀꢀ
式1
20.式中：
21.r——冰水泥石流流域地形因子；
22.x——冰水泥石流坡向因子，由式2计算；
23.a——积雪区面积，m2；
24.w——冰水泥石流形成区纵坡降为17度处的沟床平均宽度，m；
25.x＝sin(θ/2)
ꢀꢀꢀꢀꢀ
式2
26.式中：
27.θ——积雪区坡向，
°
；0
°
＜θ≤360
°
，在北半球，指母冰川的平面法线在平面投影指向与正北方向夹角；在南半球，指母冰川的平面法线在平面投影指向与正南方向夹角；
28.d、通过式3计算诱发冰水泥石流的温度因子t；
29.t＝(t0+0.27t7+1.48ty)/t
p
ꢀꢀꢀꢀꢀꢀ
式3
30.式中：
31.t——诱发冰水泥石流的温度因子；
32.t0——当日日均温度，℃；
33.t7——前7日积温和，℃；
34.ty——当年积雪消融期有效正积温，℃；
35.t
p
——当地多年平均积温，℃；
36.e、通过式4计算冰水泥石流的发生指标p；
37.p＝rt
0.28
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式4
38.式中：p——冰水泥石流的发生指标；
39.f、判断积雪消融期冰水泥石流的发生；
40.当p＜3.6时，冰水泥石流发生可能性小；当3.6≤p＜3.7时，冰水泥石流发生可能性中等；当p≥3.7时，冰水泥石流发生可能性大。
41.还包括判断积雪消融期后期冰水泥石流的发生，具体是指重新计算当年积雪消融期有效正积温ty，计算时间段为积雪消融期初期泥石流发生的那一日为积温统计的起始日，到预警时间为止；再通过式5计算诱发积雪消融期后期冰水泥石流的温度因子t1；最后通过式6计算积雪消融期后期冰水泥石流的发生指标p1；
42.t1＝(t0+0.27t
7-0.48t
y1
)/t
p
ꢀꢀꢀꢀ
式5
43.式中：
44.t1——诱发积雪消融期后期冰水泥石流的温度因子；
45.t0——当日日均温度，℃；
46.t7——前7日积温和，℃；
47.t
y1
——积雪消融期冰水泥石流发生后有效正积温，℃；
48.t
p
——当地多年平均积温，℃；
49.p1＝t1+8r
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式6
50.式中：
51.p1——积雪消融期后期冰水泥石流的发生指标；
52.t1——诱发积雪消融期后期冰水泥石流的温度因子；
53.r——冰水泥石流流域地形因子；
54.当p1＜23.6时，积雪消融期后期冰水泥石流的发生可能性小；
55.当23.6≤p1＜24.3时，积雪消融期后期冰水泥石流的发生可能性中等；
56.当p1≥24.3时，积雪消融期后期冰水泥石流的发生可能性大。
57.所述步骤b中，当年积雪消融期有效正积温ty的计算时间段为日温稳定在0℃以上的那一日为积温统计的起始日，到发生冰水泥石流时为止。
58.本发明适用于冰川及积雪区以下堆积区以上的冰水泥石流形成区的沟床坡度大于17度的沟道冰水泥石流预警。
59.本发明所述冰川及积雪区以下堆积区以上的冰水泥石流形成区是指流域中下游沟道。
60.本发明的基本原理如下：
61.冰水泥石流的形成机理是由于温度的升高，造成泥石流流域的上游冰雪融化，起动中下游的松散固体物源，形成冰水泥石流。因此这类泥石流的形成机理与我国西部地区普遍的沟床起动型泥石流不同，也与我国东部地区的浅层滑坡汇集型泥石流不同。沟床起动型泥石流是由于强降雨形成的大流量山洪侵蚀沟床松散固体物源，形成泥石流。浅层滑坡汇集型泥石流是由于较长时间的降雨和较强降雨，先有较多的浅层滑坡进入沟道形成泥石流物源，再被较大流量洪水搬运，形成泥石流。这两种类型的泥石流都需要较强的降雨强度形成较大流量洪水，侵蚀搬运沟床固体物源。但是由于温度升高造成的上游冰雪融化形成的流域洪水往往持续时间较长，但流量较小，很难像西部一般的泥石流一样，直接由洪水侵蚀搬运沟道固体物源形成泥石流。
62.在我国冰水泥石流的发生地点以西藏和新疆为主，这些地区的泥石流流域的沟道坡度较大。冰水泥石流的形成机理为：上游冰雪融化形成的流域洪水在沟道内局部水位抬升，当水位抬升较高时，造成局部松散固体物源饱水失稳，并造成失稳区域上游侵蚀，下游冲刷，最后形成泥石流。松散固体物源饱水失稳所需要的坡度远小于松散颗粒的休止角。因
此沟道坡度过小，不能达到松散固体物源饱水失稳所需要的最小坡度，冰水泥石流很难发生。如果坡度过大，上游冰雪融化形成的洪水流动过快，水位难以抬升，也很难形成泥石流。而形成这种松散固体物源饱水失稳的最小坡度就是17度。
63.在积雪消融期中，由于短时间内积雪自身产生的融水量较少，难以形成起动泥石流所需要的较大流量，所以单纯依靠短时间内积雪融水来激发泥石流较为困难。当积雪表面接受太阳辐射产生后，产生的融水沿积雪颗粒之间的空隙下渗至积雪层内部，积雪层内部由于积雪颗粒对产生的融水具有毛细力作用，所以上层产生的融水会储存在积雪层内的小空隙中，积雪层的这种特性称之为持水能力，但是积雪层中较大的空隙无法为融水提供持水能力，参见图2。当积雪层中储存的融水达到积雪层最大持水能力后，表层继续产生的融水将会打破积雪层内的持水平衡，积雪层内储存的融水将会在短时间内大量释放，为冰水泥石流的起动提供充足的水源。
64.在积雪消融期初期泥石流发生时，由于积雪层内部不存在难以提供持水能力的较大空隙，上层积雪产生的融水绝大部分进入具有持水能力的积雪层中；当温度上升达到一定条件，储存在积雪中的融水突然释放，提供足够的水流起动沟床固体物质形成冰水泥石流。随着积雪消融的持续进行，到了后期泥石流发生时，上层产生的融水大部分进入没有持水能力的较大空隙中，只有少部分进入有持水能力的积雪层中，所以对于后期冰水泥石流来说，还需要更高的气温来增加短时融水量，使得融水进入持水层，才有可能使总融水量再次达到能够激发冰水泥石流的水平。所以积雪消融期初期泥石流与积雪消融期后期泥石流有所不同。
65.本发明通过对比研究西藏和新疆公路沿线泥石流的发育与分布特征，研究得出以冰水泥石流流域中下游沟道坡度为依据的冰水泥石流沟道地形条件。本发明从冰水泥石流的形成机理出发，考虑了地形因素中的坡度这个重要因素的作用，体现出影响冰水泥石流形成的决定因素及其重要性。同时，考虑不同区域的差别，考虑了当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年积雪消融期有效正积温ty多个因素，共同组成了冰水泥石流的温度因子。最终，结合冰水泥石流的地形因子与温度因子，组成了冰水泥石流在积雪消融期的泥石流预警模型。
66.本发明的有益效果主要表现在以下方面：
67.1、本发明，a、通过现场测量确定冰水泥石流的基本地形数据，包括积雪区面积a、冰水泥石流形成区纵坡降为17度处的沟床平均宽度w和积雪区坡向θ；b、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年积雪消融期有效正积温ty；c、通过式1计算冰水泥石流流域地形因子r；d、通过式3计算诱发冰水泥石流的温度因子t；e、通过式4计算冰水泥石流的发生指标p；f、判断积雪消融期冰水泥石流的发生；当p＜3.6时，冰水泥石流发生可能性小；当3.6≤p＜3.7时，冰水泥石流发生可能性中等；当p≥3.7时，冰水泥石流发生可能性大，作为一个完整的技术方案，较现有技术而言，针对导致冰水泥石流形成的地形条件及温度条件进行内在机理的深入研究，并作出综合判断，建立积雪消融期冰水泥石流预警模型，确定冰水泥石流的预警划分级别，能够准确的对冰水泥石流的发生作出判断。
68.2、本发明，充分考虑积雪消融期积雪融化特征，给出了积雪消融期初期以及后期冰水泥石流发生的判断模型，适用性强。
69.3、本发明，以积雪区坡向θ、冰水泥石流形成区纵坡降为17度处的沟床平均宽度w和积雪区面积a，计算冰水泥石流流域地形因子，为冰水泥石流发生的地形条件提供了基本特征，保证了冰水泥石流预警的准确性。
70.4、本发明，以当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年积雪消融期有效正积温ty多个因素，计算诱发冰水泥石流的温度因子，为冰水泥石流发生的温度条件提供了基本特征，进一步提高了冰水泥石流预警的准确性。
71.5、本发明，通过将冰水泥石流的地形条件与温度条件有机结合起来，充分考虑地形与温度在冰水泥石流形成中的作用，使得预警判断更加合理，更加准确。
72.6、本发明，适用于冰川及积雪区以下堆积区以上的冰水泥石流形成区的沟床坡度大于17度的沟道冰水泥石流预警判断，冰水泥石流预警方法简单，易于执行。
附图说明
73.下面将结合说明书附图和具体实施方式对本发明作进一步的具体说明：
74.图1为本发明北半球积雪区坡向示意图；
75.图2为本发明积雪层持水能力示意图；
76.其中：θ、积雪区坡向，a、失去持水能力的较大空洞，b、尚具有持水能力的积雪层，c、融水优先排出方向。
具体实施方式
77.实施例1
78.参见图1和图2，一种积雪消融期冰水泥石流预警方法，包括以下步骤：
79.a、通过现场测量确定冰水泥石流的基本地形数据，包括积雪区面积a、冰水泥石流形成区纵坡降为17度处的沟床平均宽度w和积雪区坡向θ；
80.b、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年积雪消融期有效正积温ty；
81.c、通过式1计算冰水泥石流流域地形因子r；
82.r＝x+0.65lg(a/w2)
ꢀꢀ
式1
83.式中：
84.r——冰水泥石流流域地形因子；
85.x——冰水泥石流坡向因子，由式2计算；
86.a——积雪区面积，m2；
87.w——冰水泥石流形成区纵坡降为17度处的沟床平均宽度，m；
88.x＝sin(θ/2)
ꢀꢀꢀꢀꢀ
式2
89.式中：
90.θ——积雪区坡向，
°
；0
°
＜θ≤360
°
，在北半球，指母冰川的平面法线在平面投影指向与正北方向夹角；在南半球，指母冰川的平面法线在平面投影指向与正南方向夹角；
91.d、通过式3计算诱发冰水泥石流的温度因子t；
92.t＝(t0+0.27t7+1.48ty)/t
p
ꢀꢀꢀꢀꢀꢀ
式3
93.式中：
94.t——诱发冰水泥石流的温度因子；
95.t0——当日日均温度，℃；
96.t7——前7日积温和，℃；
97.ty——当年积雪消融期有效正积温，℃；
98.t
p
——当地多年平均积温，℃；
99.e、通过式4计算冰水泥石流的发生指标p；
100.p＝rt
0.28
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式4
101.式中：p——冰水泥石流的发生指标；
102.f、判断积雪消融期冰水泥石流的发生；
103.当p＜3.6时，冰水泥石流发生可能性小；当3.6≤p＜3.7时，冰水泥石流发生可能性中等；当p≥3.7时，冰水泥石流发生可能性大。
104.本实施例为最基本的实施方式，针对导致冰水泥石流形成的地形条件及温度条件进行内在机理的深入研究，并作出综合判断，建立积雪消融期冰水泥石流预警模型，确定冰水泥石流的预警划分级别，能够准确的对冰水泥石流的发生作出判断。
105.实施例2
106.参见图1和图2，一种积雪消融期冰水泥石流预警方法，包括以下步骤：
107.a、通过现场测量确定冰水泥石流的基本地形数据，包括积雪区面积a、冰水泥石流形成区纵坡降为17度处的沟床平均宽度w和积雪区坡向θ；
108.b、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年积雪消融期有效正积温ty；
109.c、通过式1计算冰水泥石流流域地形因子r；
110.r＝x+0.65lg(a/w2)
ꢀꢀ
式1
111.式中：
112.r——冰水泥石流流域地形因子；
113.x——冰水泥石流坡向因子，由式2计算；
114.a——积雪区面积，m2；
115.w——冰水泥石流形成区纵坡降为17度处的沟床平均宽度，m；
116.x＝sin(θ/2)
ꢀꢀꢀꢀꢀ
式2
117.式中：
118.θ——积雪区坡向，
°
；0
°
＜θ≤360
°
，在北半球，指母冰川的平面法线在平面投影指向与正北方向夹角；在南半球，指母冰川的平面法线在平面投影指向与正南方向夹角；
119.d、通过式3计算诱发冰水泥石流的温度因子t；
120.t＝(t0+0.27t7+1.48ty)/t
p
ꢀꢀꢀꢀꢀꢀ
式3
121.式中：
122.t——诱发冰水泥石流的温度因子；
123.t0——当日日均温度，℃；
124.t7——前7日积温和，℃；
125.ty——当年积雪消融期有效正积温，℃；
126.t
p
——当地多年平均积温，℃；
127.e、通过式4计算冰水泥石流的发生指标p；
128.p＝rt
0.28
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式4
129.式中：p——冰水泥石流的发生指标；
130.f、判断积雪消融期冰水泥石流的发生；
131.当p＜3.6时，冰水泥石流发生可能性小；当3.6≤p＜3.7时，冰水泥石流发生可能性中等；当p≥3.7时，冰水泥石流发生可能性大。
132.进一步的，还包括判断积雪消融期后期冰水泥石流的发生，具体是指重新计算当年积雪消融期有效正积温ty，计算时间段为积雪消融期初期泥石流发生的那一日为积温统计的起始日，到预警时间为止；再通过式5计算诱发积雪消融期后期冰水泥石流的温度因子t1；最后通过式6计算积雪消融期后期冰水泥石流的发生指标p1；
133.t1＝(t0+0.27t
7-0.48t
y1
)/t
p
ꢀꢀꢀꢀ
式5
134.式中：
135.t1——诱发积雪消融期后期冰水泥石流的温度因子；
136.t0——当日日均温度，℃；
137.t7——前7日积温和，℃；
138.t
y1
——积雪消融期冰水泥石流发生后有效正积温，℃；
139.t
p
——当地多年平均积温，℃；
140.p1＝t1+8r
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式6
141.式中：
142.p1——积雪消融期后期冰水泥石流的发生指标；
143.t1——诱发积雪消融期后期冰水泥石流的温度因子；
144.r——冰水泥石流流域地形因子；
145.当p1＜23.6时，积雪消融期后期冰水泥石流的发生可能性小；
146.当23.6≤p1＜24.3时，积雪消融期后期冰水泥石流的发生可能性中等；
147.当p1≥24.3时，积雪消融期后期冰水泥石流的发生可能性大。
148.本实施例为一较佳实施方式，充分考虑积雪消融期积雪融化特征，给出了积雪消融期初期以及后期冰水泥石流发生的判断模型，适用性强。
149.实施例3
150.参见图1和图2，一种积雪消融期冰水泥石流预警方法，包括以下步骤：
151.a、通过现场测量确定冰水泥石流的基本地形数据，包括积雪区面积a、冰水泥石流形成区纵坡降为17度处的沟床平均宽度w和积雪区坡向θ；
152.b、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年积雪消融期有效正积温ty；
153.c、通过式1计算冰水泥石流流域地形因子r；
154.r＝x+0.65lg(a/w2)
ꢀꢀ
式1
155.式中：
156.r——冰水泥石流流域地形因子；
157.x——冰水泥石流坡向因子，由式2计算；
158.a——积雪区面积，m2；
159.w——冰水泥石流形成区纵坡降为17度处的沟床平均宽度，m；
160.x＝sin(θ/2)
ꢀꢀꢀꢀꢀ
式2
161.式中：
162.θ——积雪区坡向，
°
；0
°
＜θ≤360
°
，在北半球，指母冰川的平面法线在平面投影指向与正北方向夹角；在南半球，指母冰川的平面法线在平面投影指向与正南方向夹角；
163.d、通过式3计算诱发冰水泥石流的温度因子t；
164.t＝(t0+0.27t7+1.48ty)/t
p
ꢀꢀꢀꢀꢀꢀ
式3
165.式中：
166.t——诱发冰水泥石流的温度因子；
167.t0——当日日均温度，℃；
168.t7——前7日积温和，℃；
169.ty——当年积雪消融期有效正积温，℃；
170.t
p
——当地多年平均积温，℃；
171.e、通过式4计算冰水泥石流的发生指标p；
172.p＝rt
0.28
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式4
173.式中：p——冰水泥石流的发生指标；
174.f、判断积雪消融期冰水泥石流的发生；
175.当p＜3.6时，冰水泥石流发生可能性小；当3.6≤p＜3.7时，冰水泥石流发生可能性中等；当p≥3.7时，冰水泥石流发生可能性大。
176.还包括判断积雪消融期后期冰水泥石流的发生，具体是指重新计算当年积雪消融期有效正积温ty，计算时间段为积雪消融期初期泥石流发生的那一日为积温统计的起始日，到预警时间为止；再通过式5计算诱发积雪消融期后期冰水泥石流的温度因子t1；最后通过式6计算积雪消融期后期冰水泥石流的发生指标p1；
177.t1＝(t0+0.27t
7-0.48t
y1
)/t
p
ꢀꢀꢀꢀ
式5
178.式中：
179.t1——诱发积雪消融期后期冰水泥石流的温度因子；
180.t0——当日日均温度，℃；
181.t7——前7日积温和，℃；
182.t
y1
——积雪消融期冰水泥石流发生后有效正积温，℃；
183.t
p
——当地多年平均积温，℃；
184.p1＝t1+8r
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式6
185.式中：
186.p1——积雪消融期后期冰水泥石流的发生指标；
187.t1——诱发积雪消融期后期冰水泥石流的温度因子；
188.r——冰水泥石流流域地形因子；
189.当p1＜23.6时，积雪消融期后期冰水泥石流的发生可能性小；
190.当23.6≤p1＜24.3时，积雪消融期后期冰水泥石流的发生可能性中等；
191.当p1≥24.3时，积雪消融期后期冰水泥石流的发生可能性大。
192.进一步的，所述步骤b中，当年积雪消融期有效正积温ty的计算时间段为日温稳定
在0℃以上的那一日为积温统计的起始日，到发生冰水泥石流时为止。
193.本实施例为又一较佳实施方式，以积雪区坡向θ、冰水泥石流形成区纵坡降为17度处的沟床平均宽度w和积雪区面积a，计算冰水泥石流流域地形因子，为冰水泥石流发生的地形条件提供了基本特征，保证了冰水泥石流预警的准确性。
194.实施例4
195.参见图1和图2，一种积雪消融期冰水泥石流预警方法，包括以下步骤：
196.a、通过现场测量确定冰水泥石流的基本地形数据，包括积雪区面积a、冰水泥石流形成区纵坡降为17度处的沟床平均宽度w和积雪区坡向θ；
197.b、由当地气象站或监测站点获取当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年积雪消融期有效正积温ty；
198.c、通过式1计算冰水泥石流流域地形因子r；
199.r＝x+0.65lg(a/w2)
ꢀꢀ
式1
200.式中：
201.r——冰水泥石流流域地形因子；
202.x——冰水泥石流坡向因子，由式2计算；
203.a——积雪区面积，m2；
204.w——冰水泥石流形成区纵坡降为17度处的沟床平均宽度，m；
205.x＝sin(θ/2)
ꢀꢀꢀꢀꢀ
式2
206.式中：
207.θ——积雪区坡向，
°
；0
°
＜θ≤360
°
，在北半球，指母冰川的平面法线在平面投影指向与正北方向夹角；在南半球，指母冰川的平面法线在平面投影指向与正南方向夹角；
208.d、通过式3计算诱发冰水泥石流的温度因子t；
209.t＝(t0+0.27t7+1.48ty)/t
p
ꢀꢀꢀꢀꢀꢀ
式3
210.式中：
211.t——诱发冰水泥石流的温度因子；
212.t0——当日日均温度，℃；
213.t7——前7日积温和，℃；
214.ty——当年积雪消融期有效正积温，℃；
215.t
p
——当地多年平均积温，℃；
216.e、通过式4计算冰水泥石流的发生指标p；
217.p＝rt
0.28
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式4
218.式中：p——冰水泥石流的发生指标；
219.f、判断积雪消融期冰水泥石流的发生；
220.当p＜3.6时，冰水泥石流发生可能性小；当3.6≤p＜3.7时，冰水泥石流发生可能性中等；当p≥3.7时，冰水泥石流发生可能性大。
221.还包括判断积雪消融期后期冰水泥石流的发生，具体是指重新计算当年积雪消融期有效正积温ty，计算时间段为积雪消融期初期泥石流发生的那一日为积温统计的起始日，到预警时间为止；再通过式5计算诱发积雪消融期后期冰水泥石流的温度因子t1；最后通过式6计算积雪消融期后期冰水泥石流的发生指标p1；
222.t1＝(t0+0.27t
7-0.48t
y1
)/t
p
ꢀꢀꢀꢀ
式5
223.式中：
224.t1——诱发积雪消融期后期冰水泥石流的温度因子；
225.t0——当日日均温度，℃；
226.t7——前7日积温和，℃；
227.t
y1
——积雪消融期冰水泥石流发生后有效正积温，℃；
228.t
p
——当地多年平均积温，℃；
229.p1＝t1+8r
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
式6
230.式中：
231.p1——积雪消融期后期冰水泥石流的发生指标；
232.t1——诱发积雪消融期后期冰水泥石流的温度因子；
233.r——冰水泥石流流域地形因子；
234.当p1＜23.6时，积雪消融期后期冰水泥石流的发生可能性小；
235.当23.6≤p1＜24.3时，积雪消融期后期冰水泥石流的发生可能性中等；
236.当p1≥24.3时，积雪消融期后期冰水泥石流的发生可能性大。
237.所述步骤b中，当年积雪消融期有效正积温ty的计算时间段为日温稳定在0℃以上的那一日为积温统计的起始日，到发生冰水泥石流时为止。
238.本发明适用于冰川及积雪区以下堆积区以上的冰水泥石流形成区的沟床坡度大于17度的沟道冰水泥石流预警。
239.本发明所述冰川及积雪区以下堆积区以上的冰水泥石流形成区是指流域中下游沟道。
240.本实施例为最佳实施方式，以当地多年平均积温t
p
、前7日积温和t7、当日日均温度t0和当年积雪消融期有效正积温ty多个因素，计算诱发冰水泥石流的温度因子，为冰水泥石流发生的温度条件提供了基本特征，进一步提高了冰水泥石流预警的准确性。
241.通过将冰水泥石流的地形条件与温度条件有机结合起来，充分考虑地形与温度在冰水泥石流形成中的作用，使得预警判断更加合理，更加准确。
242.适用于冰川及积雪区以下堆积区以上的冰水泥石流形成区的沟床坡度大于17度的沟道冰水泥石流预警判断，冰水泥石流预警方法简单，易于执行。
243.下面以具体实例对本发明进行说明：
244.新疆天山公路g217位于我国西北部的天山山脉和昆仑山脉：g217“独山子——库车”，亦称独库公路段是一条南北疆的重要交通要道，独库公路北起准格尔盆地西南边缘的石油重地——克拉玛依市独山子区，南迄于天山山脉南麓阿克苏市库车县，是由北向南横跨整个天山山脉，其2/3左右的里程穿越海拔2000m以上的寒冻山区，全线穿过4次高山区，途径天山5条较大的河流。根据调查，泥石流在天山公路全线均有分布，但主要集中在北段和南段。天山公路全线各类泥石流共有231处，包括坡面泥石流和沟谷泥石流，其中北段有53处，中段有3处，南段有178处，其中，对公路影响较严重的泥石流有89处、对公路构重大影响的泥石流共有12处，其中10处位于北段，主要分布于k629-k660之间。
245.在g217公路k629-k660段的泥石流沟中，以k636泥石流沟的冰水泥石流灾害最为严重，在20世纪70年代对g217公路沿线进行勘测时，将k636这条最大泥石流沟误认为是不
活跃的泥石流沟，未作防灾工程处理。但是随着近年来全球气温的普遍升高，研究区内的冰雪层消融进一步加剧，在夏季高温时，流域内冰雪融水量开始增大，导致g217公路k636冰水泥石流沟爆发频繁，严重影响了公路畅通。从1984年至今，k636沟由于高温引发的积雪消融期初期冰水泥石流共发生2次。
246.以g217公路k636泥石流沟发生积雪消融期初期泥石流的日期，以及k636泥石流沟和周围泥石流中，未发生积雪消融期初期泥石流且气温较高的日期为例，对冰水泥石流进行预警判断。表1为g217公路k636泥石流沟及周围泥石流沟道中，高温条件下冰水泥石流的29个案例的参数及相应的冰水泥石流的发生指标p，以及泥石流的实际发生情况。由于所选泥石流沟道距离较近，所以当地多年平均积温t
p
采用相同值，取14.13℃。
247.表1
248.[0249][0250]
判断积雪消融期冰水泥石流的发生：
[0251]
当p＜3.6时，冰水泥石流发生可能性小；当3.6≤p＜3.7时，冰水泥石流发生可能性中等；当p≥3.7时，冰水泥石流发生可能性大。
[0252]
通过表1中p值计算结果结合实际发生情况显示：判断冰水泥石流发生可能性大的有1个，实际发生泥石流；判断冰水泥石流发生可能性中等的有2个，实际一个发生泥石流，一个未发生泥石流；判断冰水泥石流发生可能性小的有26个，实际都未发生泥石流。
[0253]
以g217公路k636泥石流沟发生积雪消融期后期泥石流的日期，以及k636泥石流沟和周围泥石流中，未发生积雪消融期后期泥石流且气温较高的日期为例，对冰水泥石流进行预警判断。表2为g217公路k636泥石流沟及周围泥石流沟道中，高温条件下冰水泥石流的33个案例的参数及相应的泥石流发生指标p1，以及泥石流的实际发生情况。从1984年至今，k636沟由于高温引发的积雪消融期后期冰水泥石流共发生2次。由于所选泥石流沟道距离较近，所以当地多年平均积温t
p
采用相同值，取14.13℃。
[0254]
表2
[0255]
[0256][0257]
判断积雪消融期后期冰水泥石流的发生：
[0258]
当p1＜23.6时，积雪消融期后期冰水泥石流的发生可能性小；
[0259]
当23.6≤p1＜24.3时，积雪消融期后期冰水泥石流的发生可能性中等；
[0260]
当p1≥24.3时，积雪消融期后期冰水泥石流的发生可能性大。
[0261]
通过表2中p1值计算结果结合实际发生情况显示：判断积雪消融期后期冰水泥石流的发生可能性大的有1个，实际发生泥石流；判断积雪消融期后期冰水泥石流的发生可能性中等的有3个，实际一个发生泥石流，2个未发生泥石流；判断积雪消融期后期冰水泥石流的发生可能性小的有29个，实际都未发生泥石流。
[0262]
综上，应用本发明预警方法对积雪消融期初期和后期泥石流的预警判断都较为准确。