一种管涌型堰塞体稳定性判断方法

1.本发明涉及到泥石流防治工程技术领域，尤其涉及一种管涌型堰塞体稳定性判断方法。


背景技术：

2.山区沿江两岸的大型斜坡体或高边坡，在地震、暴雨荷载下极易发生破坏失稳，造成堵江并形成堰塞坝。不同于人工土石坝，堰塞坝坝体结构松散，颗粒级配不均匀，在较高水头作用下坝体可能发生管涌破坏而导致溃坝。因此对堰塞体的稳定性和动力学特征进行研究，预测堰塞体的溃决演化过程，为堰塞体应急抢险和地质处置提供技术支持，具有重要理论价值和现实意义。
3.公开号为cn112418697a，公开日为2021年02月26日的中国专利文献公开了一种堰塞湖稳定性评价方法，其特征在于：所述堰塞湖稳定性评价方法包括以下步骤：
4.s1：先采用无量纲堆积体指数法dbi进行堰塞体稳定性初步评判；
5.s2：针对步骤s1中得到的结果进行分析：dbi＜2.75，则堰塞体稳定；dbi＞3.08，则堰塞体不稳定；2.75＜dbi＜3.08，则堰塞体介于稳定与不稳定之间；
6.s3：评价堰塞体的稳定性要结合地质条件，合理选取堰塞体物质的物理力学指标和渗流场参数；
7.s4：可以采用基于极限平衡理论体系的简化毕肖普法，用简化毕肖谱法计算时考虑了土条间的作用力，可得到堰塞体的稳定性评价结论；
8.s5：对堰塞湖进行严密监测，根据专业人员的意见对堰塞湖的威胁进行正确科学的评价，及时对下游人民群众进行疏散；
9.s6：然后，对于高危堰塞湖进行疏导和分流，一般是用炸药和人工挖开个小口，疏导水流。
10.该专利文献公开的堰塞湖稳定性评价方法，用简化毕肖谱法计算时考虑了土条间的作用力，可得到堰塞体的稳定性评价结论；采用无量纲堆积体指数法dbi进行堰塞体稳定性初步评判，使得在后续的判断过程中，有了一个可以参考的初步依据，从而简化了操作步骤，减轻了工作人员的劳动负担。但是，由于没有考虑堰塞体在渗流作用下的管涌破坏的动力学特征，因此准确度较低，普适性较差。


技术实现要素：

11.本发明为了克服上述现有技术的缺陷，提供一种管涌型堰塞体稳定性判断方法，本发明考虑了在渗流作用下骨架孔隙中的可动颗粒启动破坏的动力学特征，选择可动颗粒进行力学平衡分析，提高了计算准确度，能够为堰塞体灾害的应急抢险及危险性评估提供更合理准确的依据，具有良好的普适性。
12.本发明通过下述技术方案实现：
13.一种管涌型堰塞体稳定性判断方法，其特征在于，包括以下步骤：
14.a、通过野外调查获取堰塞体的基本资料，包括堰塞体几何尺寸、堰塞体级配、堰塞体的重度和堰塞体的内摩擦角，建立堰塞体内部颗粒在渗流力作用下的受力模型；
15.b、计算堰塞体可动颗粒的阻抗力；
16.c、计算堰塞体可动颗粒所受的渗流力；
17.d、将堰塞体可动颗粒所受的渗流力与堰塞体可动颗粒的阻抗力的比值定义为堰塞体的稳定性系数，当渗流力大于摩擦力时，判断为堰塞体失稳。
18.所述步骤b中，堰塞体可动颗粒的阻抗力是指摩擦力，通过式1进行计算；
[0019][0020]
其中：
[0021]
f为摩擦力，kn；
[0022]
γs为堰塞体的重度，kn/m3；
[0023]
γw为水的重度，kn/m3；
[0024]
为堰塞体的内摩擦角，
°
；
[0025]
d为可动颗粒的粒径，m；取太沙基控制粒径d1为可动颗粒的最大粒径。
[0026]
所述步骤d中，渗流力f
p
通过式2进行计算；
[0027][0028]
其中：
[0029]fp
为渗流力，kn；
[0030]
γw为水的重度，kn/m3；
[0031]
i为水力梯度，即渗透水流沿单位渗透距离水头的损失量，无量纲，按式3求取；
[0032]
s为单位体积堰塞体通道中可动颗粒的表面积，m2，按式4求取；
[0033][0034]
其中：
[0035]hw
为坝后堰塞湖水深，m；
[0036]
i为水力梯度，即渗透水流沿单位渗透距离水头的损失量，无量纲；
[0037]
l为堰塞体长度，m；
[0038]
η为堰塞体沟床坡度，
°
；
[0039][0040]
其中：
[0041]
s为单位体积堰塞体通道中可动颗粒的表面积，m2，
[0042]
α为颗粒形状系数，无量纲，对于砂粒取1.16-1.17，锐角颗粒取1.50-1.67；
[0043]
n为孔隙率，无量纲；
[0044]
k为堰塞体的渗透系数，m/s；
[0045]
g为重力加速度，取9.81m/s2。
[0046]
所述步骤d中，堰塞体的稳定性系数fs通过式5进行计算；
[0047][0048]
其中：
[0049]
fs为堰塞体的稳定性系数；
[0050]
f为摩擦力，kn；
[0051]fp
为渗流力，kn。
[0052]
本发明的基本原理如下：
[0053]
参见图2-图5，通过室内试验研究了堰塞体的管涌溃决过程的动力学特征。管涌破坏模式可概括为管涌—小型崩滑—溢流侵蚀—坝体稳定四个阶段。随着上游来流的不断涌入，坝后水位逐渐升高，导致坝体内孔隙水压力增加，细小的可动颗粒在渗流力的作用下在较粗的骨架颗粒形成的通道之间移动，并被逐渐冲出堰塞体，形成管涌；而后，堰塞体管涌出口处发生小型崩滑，降低了堰塞体的整体稳定性；坝后水位上升至坝顶后出现漫顶溢流，堰塞体在坡面水流的冲刷下形成溃口，同时伴随着强烈的侧向展宽和下切侵蚀，溃口流量快速增加并到达峰值；随着堰塞体后水量逐渐减小，断面流量到达峰值之后，流量迅速减小直至趋于稳定。
[0054]
首先分析了堰塞体在渗流作用下的启动机理，骨架孔隙中的可动颗粒受到的力有：渗流水流的拖曳力、水下重力、与孔隙壁之间的摩擦力和颗粒之间的相互作用力。对于堰塞体内部发生的管涌而言，只有渗透水对可动颗粒的拖曳力为渗流力，其它力都为阻力，且可动颗粒在孔隙壁中移动时，摩擦力是最大的阻力，而颗粒之间的碰撞力和电场力都非常小且十分复杂。假设渗流通道水，用重力沿渗流通道向下的渗流力抵消颗粒之间的碰撞力和电场力这些阻力，阻止颗粒起动以管道之间的摩擦力为主。将可动颗粒所受渗流力与摩擦力的比值定义为堰塞体的稳定性系数，当渗流力大于摩擦力时，判断为堰塞体失稳。
[0055]
本发明的有益效果主要表现在以下方面：
[0056]
1、本发明，a、通过野外调查获取堰塞体的基本资料，包括堰塞体几何尺寸、堰塞体级配、堰塞体的重度和堰塞体的内摩擦角，建立堰塞体内部颗粒在渗流力作用下的受力模型；b、计算堰塞体可动颗粒的阻抗力；c、计算堰塞体可动颗粒所受的渗流力；d、将堰塞体可动颗粒所受的渗流力与堰塞体可动颗粒的阻抗力的比值定义为堰塞体的稳定性系数，当渗流力大于摩擦力时，判断为堰塞体失稳；较现有技术而言，考虑了在渗流作用下骨架孔隙中的可动颗粒启动破坏的动力学特征，选择可动颗粒进行力学平衡分析，提高了计算准确度，能够为堰塞体灾害的应急抢险及危险性评估提供更合理准确的依据，具有良好的普适性。
[0057]
2、本发明，是基于力学平衡分析的理论推导而来，经过了室内物理试验的验证，并与野外实测数据作对比，对堰塞体稳定性确定具有较高的科学参考价值。
[0058]
3、本发明，计算公式简单，计算所需参数少，且参数易于获取，计算准确度高；能够适用于在缺少前期大量勘察工作情况下的堰塞体渗流稳定性快速判别。
附图说明
[0059]
下面将结合说明书附图和具体实施方式对本发明作进一步的具体说明：
[0060]
图1为本发明的流程框图；
[0061]
图2为堰塞体形成管涌示意图；
[0062]
图3为堰塞体发生小型崩滑示意图；
[0063]
图4为堰塞体漫顶示意图；
[0064]
图5为堰塞体管涌破坏示意图。
具体实施方式
[0065]
实施例1
[0066]
一种管涌型堰塞体稳定性判断方法，包括以下步骤：
[0067]
a、通过野外调查获取堰塞体的基本资料，包括堰塞体几何尺寸、堰塞体级配、堰塞体的重度和堰塞体的内摩擦角，建立堰塞体内部颗粒在渗流力作用下的受力模型；
[0068]
b、计算堰塞体可动颗粒的阻抗力；
[0069]
c、计算堰塞体可动颗粒所受的渗流力；
[0070]
d、将堰塞体可动颗粒所受的渗流力与堰塞体可动颗粒的阻抗力的比值定义为堰塞体的稳定性系数，当渗流力大于摩擦力时，判断为堰塞体失稳。
[0071]
本实施例为最基本的实施方式，考虑了在渗流作用下骨架孔隙中的可动颗粒启动破坏的动力学特征，选择可动颗粒进行力学平衡分析，提高了计算准确度，能够为堰塞体灾害的应急抢险及危险性评估提供更合理准确的依据，具有良好的普适性。
[0072]
实施例2
[0073]
一种管涌型堰塞体稳定性判断方法，包括以下步骤：
[0074]
a、通过野外调查获取堰塞体的基本资料，包括堰塞体几何尺寸、堰塞体级配、堰塞体的重度和堰塞体的内摩擦角，建立堰塞体内部颗粒在渗流力作用下的受力模型；
[0075]
b、计算堰塞体可动颗粒的阻抗力；
[0076]
c、计算堰塞体可动颗粒所受的渗流力；
[0077]
d、将堰塞体可动颗粒所受的渗流力与堰塞体可动颗粒的阻抗力的比值定义为堰塞体的稳定性系数，当渗流力大于摩擦力时，判断为堰塞体失稳。
[0078]
所述步骤b中，堰塞体可动颗粒的阻抗力是指摩擦力，通过式1进行计算；
[0079][0080]
其中：
[0081]
f为摩擦力，kn；
[0082]
γs为堰塞体的重度，kn/m3；
[0083]
γw为水的重度，kn/m3；
[0084]
为堰塞体的内摩擦角，
°
；
[0085]
d为可动颗粒的粒径，m；取太沙基控制粒径d1为可动颗粒的最大粒径。
[0086]
所述步骤d中，渗流力f
p
通过式2进行计算；
[0087]
[0088]
其中：
[0089]fp
为渗流力，kn；
[0090]
γw为水的重度，kn/m3；
[0091]
i为水力梯度，即渗透水流沿单位渗透距离水头的损失量，无量纲，按式3求取；
[0092]
s为单位体积堰塞体通道中可动颗粒的表面积，m2，按式4求取；
[0093][0094]
其中：
[0095]hw
为坝后堰塞湖水深，m；
[0096]
i为水力梯度，即渗透水流沿单位渗透距离水头的损失量，无量纲；
[0097]
l为堰塞体长度，m；
[0098]
η为堰塞体沟床坡度，
°
；
[0099][0100]
其中：
[0101]
s为单位体积堰塞体通道中可动颗粒的表面积，m2，
[0102]
α为颗粒形状系数，无量纲，对于砂粒取1.16-1.17，锐角颗粒取1.50-1.67；
[0103]
n为孔隙率，无量纲；
[0104]
k为堰塞体的渗透系数，m/s；
[0105]
g为重力加速度，取9.81m/s2。
[0106]
本实施例为较佳实施方式，基于力学平衡分析的理论推导，经过室内物理试验的验证，并与野外实测数据作对比，对堰塞体稳定性确定具有较高的科学参考价值。
[0107]
实施例3
[0108]
一种管涌型堰塞体稳定性判断方法，包括以下步骤：
[0109]
a、通过野外调查获取堰塞体的基本资料，包括堰塞体几何尺寸、堰塞体级配、堰塞体的重度和堰塞体的内摩擦角，建立堰塞体内部颗粒在渗流力作用下的受力模型；
[0110]
b、计算堰塞体可动颗粒的阻抗力；
[0111]
c、计算堰塞体可动颗粒所受的渗流力；
[0112]
d、将堰塞体可动颗粒所受的渗流力与堰塞体可动颗粒的阻抗力的比值定义为堰塞体的稳定性系数，当渗流力大于摩擦力时，判断为堰塞体失稳。
[0113]
所述步骤b中，堰塞体可动颗粒的阻抗力是指摩擦力，通过式1进行计算；
[0114][0115]
其中：
[0116]
f为摩擦力，kn；
[0117]
γs为堰塞体的重度，kn/m3；
[0118]
γw为水的重度，kn/m3；
[0119]
为堰塞体的内摩擦角，
°
；
[0120]
d为可动颗粒的粒径，m；取太沙基控制粒径d1为可动颗粒的最大粒径。
[0121]
所述步骤d中，渗流力f
p
通过式2进行计算；
[0122][0123]
其中：
[0124]fp
为渗流力，kn；
[0125]
γw为水的重度，kn/m3；
[0126]
i为水力梯度，即渗透水流沿单位渗透距离水头的损失量，无量纲，按式3求取；
[0127]
s为单位体积堰塞体通道中可动颗粒的表面积，m2，按式4求取；
[0128][0129]
其中：
[0130]hw
为坝后堰塞湖水深，m；
[0131]
i为水力梯度，即渗透水流沿单位渗透距离水头的损失量，无量纲；
[0132]
l为堰塞体长度，m；
[0133]
η为堰塞体沟床坡度，
°
；
[0134][0135]
其中：
[0136]
s为单位体积堰塞体通道中可动颗粒的表面积，m2，
[0137]
α为颗粒形状系数，无量纲，对于砂粒取1.16-1.17，锐角颗粒取1.50-1.67；
[0138]
n为孔隙率，无量纲；
[0139]
k为堰塞体的渗透系数，m/s；
[0140]
g为重力加速度，取9.81m/s2。
[0141]
所述步骤d中，堰塞体的稳定性系数fs通过式5进行计算；
[0142][0143]
其中：
[0144]
fs为堰塞体的稳定性系数；
[0145]
f为摩擦力，kn；
[0146]fp
为渗流力，kn。
[0147]
本实施例为最佳实施方式，计算公式简单，计算所需参数少，且参数易于获取，计算准确度高；能够适用于在缺少前期大量勘察工作情况下的堰塞体渗流稳定性快速判别。
[0148]
下面结合具体实例对本发明的实施方式进行详细说明：
[0149]
四川省青川县红石河堰塞湖是汶川大地震造成的大型堰塞湖之一，由东河口滑坡堵塞红石河形成。红石河堰塞湖坝体高约50m、坝长200-300m、底宽400m，堰塞体总体呈右高
左低的走势。其上下游水位差30m，蓄水量约200万m3，最大库容可达400万m3，集雨面积70km2。堰塞体材料为第四系全新统崩塌堆积物，结构为块石夹黏土，以块石为主，块碎石粒径5-30cm，个别粒径较大，可达8-9m。因红石河堰塞坝体材料的物源及颗粒物质组成非常复杂，局部可存在颗粒支撑型材料，即渗透性较高的区域。
[0150]
对该堰塞体渗流稳定性进行判别，具体步骤如下：
[0151]
第一步，获取堰塞体的基本资料，建立堰塞体内部颗粒在渗流力作用下的受力模型，红石河堰塞体的基本参数见表1；
[0152]
表1
[0153][0154][0155]
第二步，通过式1计算摩擦力；
[0156][0157]
其中：
[0158]
f为摩擦力，kn；
[0159]
γs为堰塞体的重度，kn/m3；
[0160]
γw为水的重度，kn/m3；
[0161]
为堰塞体的内摩擦角，
°
；
[0162]
d为可动颗粒的粒径，m；取太沙基控制粒径d1为可动颗粒的最大粒径；
[0163]
第三步，通过式2计算渗流力f
p
；
[0164][0165]
其中：
[0166]fp
为渗流力，kn；
[0167]
γw为水的重度，kn/m3；
[0168]
计算结果见表2，表2为红石河堰塞体稳定性计算结果。
[0169]
表2
[0170]
摩擦力f(
×
10-7
kn)渗流力f
p
(
×
10-7
kn)fs实际情况4.3942.6981.629未溃决
[0171]
可见，采用本发明稳定性判断方法，判断准确度高，具有良好的普适性。