一种降雨条件下的边坡三维稳定性预测方法
【专利摘要】本发明公开了一种降雨条件下的边坡三维稳定性预测方法，包括以下步骤：1、根据待预测的边坡，确定土体强度参数和土体渗透参数；2、对滑坡体划分土柱，沿竖直方向滑坡体被划分为行和列土柱；3、利用有限元软件GEO/SEEP，建立边坡渗流模型，计算滑坡体内的孔隙水压力分布情况；4、建立降雨条件下边坡三维极限平衡方程组，通过联立方程组解得的边坡稳定系数Fs；5、求得不同降雨强度、降雨持时和降雨型式下边坡的最小稳定系数，将该稳定系数与边坡稳定临界值进行比较。本发明的优点是：考虑了降雨条件对边坡稳定性的影响，计算精度有所提高，预测结果更为可靠。
【专利说明】
一种降雨条件下的边坡三维稳定性预测方法
技术领域
[0001] 本发明属于地质灾害防控技术领域，具体涉及一种针对降雨条件下预测边坡三维 稳定性的方法。
【背景技术】
[0002] 随着全球气候变暖，世界范围内出现极端降雨天气的频率越来越高，由于我国山 地、丘陵地带分布广泛，因而由于降雨引起的滑坡灾害占每年地质灾害的百分之八十以上。 降雨前边坡处于安全状态，降雨过程中或者降雨后一定时间范围内出现滑坡现象，这使得 预测降雨条件下边坡的稳定性问题变得尤为迫切。
[0003] 中国专利文献CN103149340A于2013年6月12日公开了一种用降雨量测定滑坡稳定 性的动力监测方法，包括以下步骤:步骤1、选取滑坡位移监测点与基准点；步骤2、在基准点 和监测点布置安装监测设备；步骤3、通过监测设备能够获得降雨量与位移速率的相关数 据，对数据进行处理;步骤4、确定滑坡动力加载率与位移动力响应率，在此处动力加载率是 某一降雨过程中降雨量的增量变化与初始降雨量的比值，位移动力响应率是相应的滑坡位 移变化速率与初始位移速率的比值;步骤5、确定滑坡稳定性的判据，将滑坡的位移降雨动 力加载率作为滑坡稳定性的判据，在此处将滑坡失稳的动力加载率与位移动力响应率的比 值成为滑坡位移动力加载率;步骤6、进行滑坡的监测与预警预报。该方法利用实时测量的 方法，建立降雨量与位移的关系，利用位移变化速率来预测预报边坡的稳定性。该方法基准 点选择的准确与否与预报效果有很大关系，同时通过降雨量与边坡位移的关系来预测具有 一定的不确定性。
【发明内容】
[0004] 本发明所要解决的技术问题就是提供一种降雨条件下的边坡三维稳定性预测方 法，它在预测降雨条件下边坡三维稳定性时，降雨条件包括降雨强度、降雨持时和降雨型 式，得出不同岩土体边坡在不同降雨条件下的稳定系数，能够提高计算精度，增大预测结果 的可靠性。
[0005] 本发明所要解决的技术问题是通过这样的技术方案实现的，本发明包括以下步 骤：
[0006] 步骤1、根据待预测的边坡，确定土体强度参数和土体渗透参数，选定边坡坡体表 面几何尺寸和三维滑面形状，将滑裂面形状和边坡表面形状用代数方程表示；确定需要进 行渗流分析截面的几何尺寸，设定不同降雨强度、降雨持时和降雨型式的组合条件，进行渗 流分析；
[0007] 步骤2、将边坡进行离散化，滑坡体被垂直离散为m行和η列土柱，每个土柱按所在 的行号i和列号j定义为[i、j];假定行方向的土柱间作用力与水平面夹角为±〇，假定列方 向的土柱间作用力与水平面夹角为±β;
[0008] 步骤3、利用有限元软件GE0/SEEP，建立边坡渗流模型，计算滑坡体内的孔隙水压 力分布情况，得到孔隙水压力分布图，获得各位置点的孔隙水压力；
[0009] 步骤4、建立降雨条件下边坡三维极限平衡方程组，通过联立方程组解得的边坡稳 定系数Fs;
[0010] 步骤5、求得不同降雨强度、降雨持时和降雨型式下边坡的最小稳定系数，将该稳 定系数与边坡稳定临界值进行比较，若该最小稳定系数大于边坡稳定的临界值，则该边坡 处于稳定状态;若该最小稳定系数小于边坡稳定的临界值，则该边坡处于失稳状态。
[0011] 由于本发明考虑了比二维情况更真实的三维情况，同时在离散化过程中考虑了更 多力和力矩的平衡关系，使得计算精度提高，所以本发明的预测结果更准确可靠。另外，所 有建模过程都程序化，便于操作和编程，大大的减少了人为的计算量，由计算机实现边坡的 稳定系数的预测。所以本发明的优点是:考虑了降雨条件对边坡稳定性的影响，计算精度有 所提高，预测结果更为可靠。
【附图说明】
[0012] 本发明的【附图说明】如下：
[0013]图1为本发明实施例的土体饱和渗透系数曲线图；
[0014] 图2为本发明实施例的土体含水率曲线图；
[0015] 图3为本发明实施例的边坡三维几何尺寸；
[0016] 图4为图3的渗流剖截面的几何尺寸；
[0017] 图5为滑坡体滑裂面和边坡表面的剖面图；
[0018] 图6为本发明实施例的孔隙水压力分布图；
[0019]图7为离散化土柱受力图。
【具体实施方式】
[0020] 下面结合附图和实施例对本发明作进一步说明：
[0021] 步骤1、根据待预测的边坡，确定土体强度参数和土体渗透参数，土体强度参数包 括土体重度γ、有效粘聚力C'、孔隙气压力μ α、岩土体的有效内摩擦角炉'和抗剪强度随基质 吸力变化的内摩擦角土体渗透参数包括饱和渗透系数曲线和土体体积含水率曲线;选 定边坡坡体表面几何尺寸和三维滑面形状，将滑裂面形状和边坡表面形状用代数方程表 示;确定需要进行分析的渗流截面几何尺寸，设定不同降雨强度、降雨持时和降雨型式的组 合条件，进行渗流分析；
[0022] 本实施例中，待预测边坡是典型的残坡积土边坡，土体重度γ为18.7kg · ΠΓ3,有效 粘聚力c '为7 · OkPa，有效内摩擦角供'为26 · 0°，取孔隙气压力uaS〇Kpa，抗剪强度随基质吸 力变化的内摩擦角/为定值16°。根据"降雨条件下土质边坡瞬态稳定性分析"，吴长富，岩 土力学报，第29卷第2期，387~389页，2008年2月的记载:土体的饱和渗透系数如图1所示， 饱和渗透系数为0.415m · cf1;土体含水率曲线如图2所示。选取降雨条件为暴雨98mm均强型 工况，降雨持续时间为24h。降雨入渗到土体内部，使得土体内浸润线不断上升，影响土体内 的基质吸力，使其不断减小，同时增大了土体内孔隙水压力，这两个因素引起土体抗剪强度 较低。
[0023] 三维滑坡体如图3所示，边坡表面为倾斜面，长高比为1:2,边坡斜面的水平方向的 投影为20m，在竖直方向上的投影为10m;滑裂面为椭球体滑裂面，设y方向上滑裂面宽度的 半长轴为40m;X、z方向的半轴长度为21m。渗流剖截面为三维滑坡体的横截面，其尺寸如图4 所示。
[0024] 参照中国专利文献CN103163563A于2013年6月19日公开的一种三维边坡稳定性预 测方法的步骤1:将边坡表面和滑裂面用方程表示。如图5所示，边坡表面的参数有边坡斜面 在水平面上的投影长度1，边坡斜面在竖直方向上的投影长度H;滑裂面的参数依据实际几 何形状确定。
[0025]边坡表面的表达式为
[0027]滑裂面的表达式为：
[0029] 步骤2、参照中国专利文献CN103163563A于2013年6月19日公开的一种三维边坡稳 定性预测方法的步骤2:对三维滑坡体离散化，沿竖直方向滑坡体被划分为m行和η列土柱， 每个土柱按所在的行号i和列号j定义为[i，j];假定行方向的土柱间作用力与水平面夹角 为土α，假定列方向的土柱间作用力与水平面夹角为土β;
[0030]本实施例中，三维滑坡体被垂直离散为30行和30列土柱。
[0031]步骤3、利用有限元软件GE0/SEEP，建立边坡渗流模型，如图4所示，该边坡的边界 条件设置为:abed为降雨入渗边界、ah、de、gf为0(m/s)流量边界，hg、ef为初始水头边界6m， 将步骤1中的岩土参数代入，计算即可得到边坡的孔隙水压力分布图，如图6所示。
[0032]步骤4、对每个土柱进行受力分析，土柱的受力如图7所示。根据力的平衡，每个土 柱上可建立如下的三个方向上力的平衡方程：
[0036]式（1)、（2)和(3)中：
[0037] Ni>j为第[i，j] 土柱底面支持力；
[0038] 1^」为第[1，]_]土柱底面孔隙水压力的合力；
[0039] 1^'」为第[丨，幻土柱底面孔隙气压力的合力；
[0040] Si>j为第[i，j] 土柱底面的剪切力；
[0041] Q1，」为第[i，j-l] 土柱与第[i，j] 土柱之间的条间力；
[0042] QU+1为第[i，j] 土柱与第[i，j+l] 土柱之间的条间力；
[0043] G1，」为第[i_l，j] 土柱与第[i，j] 土柱之间的条间力；
[0044] G1+U为第[i，j] 土柱与第[i+l，j] 土柱之间的条间力；
[0045] ,'」为第[丨，幻土柱的重力；
[0046] α行方向的土柱间作用力与水平面夹角；
[0047] β列方向的土柱间作用力与水平面夹角；
[0048] 表示第[i，j ] 土柱底面上正应力的X方向余弦；
[0049] 表示第[i，j]土柱底面上正应力的y方向余弦；
[0050] 表示第[i, j]土柱底面上正应力的z方向余弦；
[0051 ] f表示第[i，j]土柱底面上剪应力的X方向余弦；
[0052] t表示第[i，j]土柱底面上剪应力的y方向余弦；
[0053] 表示第[i, j]土柱底面上剪应力的z方向余弦；
[0054]设该边坡没有加固措施，则在滑体上没有支护力的作用，那么条间力的边界条件 可以表示为
（m为土柱的行总数和η为土柱的列总数）。 因为条间力是相邻两个土柱之间的作用力与反作用力，其绕X轴方向的力矩之和等于0。 [0055]利用非饱和土强度准则，并引入稳定系数表达式，得到下式：
[0057]式(4)中，
[0058] c '为岩土体的有效粘聚力；
[0059] 梦'为岩土体的有效内摩擦角；
[0060] #为抗剪强度随基质吸力变化的内摩擦角；
[0061] Α1，」为[i，j]土柱横截面的面积；
[0062] Fs为稳定系数。
[0063]将前面(1)、（2)、（3)、（4)联立求解，可得每个土柱的支持力和剪切力：
[0066]所有力绕x、y、z轴方向的力矩分别为：
[0070] 式(7)、(8)、(9)中：
[0071] AGi>j是平面ABB'A'和平面CDD'C'之间的力之差，AGi'LGMj-GU;
[0072] &9^」是平面8〇：'8'和平面4001'之间的力之差，&9^」=9^」+1-9^ ;
[0073] χΔ??ν/ΔΡ是力Λ Gi>j或者Δ QU的X坐标，其值等于点Η的X坐标值；
[0074] 是力AG1，域者AQ1，啪y坐标，其值等于点油％坐标值；
[0075] 是力Λ Gi>j或者Δ QU的z坐标，其值等于点Η的z坐标值；
[0076] <是点Η的X坐标值; 是点Η的y坐标值。
[0077] 将式(5)、（6)代入上面式(7)、（8)、（9)，可求得边坡在不同时刻的稳定系数。稳定 系数通过多次迭代求得。
[0078] 步骤5、计算暴雨98mm均强型降雨工况下的边坡稳定系数。
[0079] 表1暴雨98mm均强型降雨工况下的边坡稳定系数
[0080]
[0081]将计算的不同时刻边坡的稳定系数，与边坡稳定性临界值进行比较;若计算出的 稳定系数大于边坡稳定性临界值，则边坡处于稳定状态;若计算出的稳定系数小于边坡稳 定性临界值，则边坡处于失稳状态。本实例计算出的最危险时刻是在降雨结束时，稳定系数 最小，且接近1，该边坡很有可能出现失稳情况。
【主权项】
1. 一种降雨条件下边坡的Ξ维稳定性预测方法，其特征是，包括w下步骤： 步骤1、根据待预测的边坡，确定±体强度参数和±体渗透参数，选定边坡坡体表面几 何尺寸和Ξ维滑面形状，将滑裂面形状和边坡表面形状用代数方程表示;确定需要进行渗 流分析截面的几何尺寸，设定不同降雨强度、降雨持时和降雨型式的组合条件，进行渗流分 析； 步骤2、将边坡进行离散化，滑坡体被垂直离散为m行和11列±柱，每个±柱按所在的行 号i和列号j定义为[i、j];假定行方向的±柱间作用力与水平面夹角为±〇,假定列方向的 ±柱间作用力与水平面夹角为±β; 步骤3、利用有限元软件GEO/SEEP，建立边坡渗流模型，计算滑坡体内的孔隙水压力分 布情况，得到孔隙水压力分布图，获得各位置点的孔隙水压力； 步骤4、建立降雨条件下边坡Ξ维极限平衡方程组，通过联立方程组解得的边坡稳定系 数Fs; 步骤5、求得不同降雨强度、降雨持时和降雨型式下边坡的最小稳定系数，将该稳定系 数与边坡稳定临界值进行比较，若该最小稳定系数大于边坡稳定的临界值，则该边坡处于 稳定状态;若该最小稳定系数小于边坡稳定的临界值，则该边坡处于失稳状态。2. 根据权利要求1所述的降雨条件下边坡Ξ维稳定性预测方法，其特征是，步骤4中，沿 Ξ个轴向方向的力的平衡方程为：式(1)、（2)和(3)中： Nl'^为第[i，j]±柱底面支持力； 1^'^为第[1，引±柱底面孔隙水压力的合力； 山1'^为第[1，引±柱底面孔隙气压力的合力； 81'^为第[1，引±柱底面的剪切力； Qi'j为第[iJ-1]上柱与第[ij]上柱之间的条间力； Qi'w为第[ij]上柱与第[iJ+1]上柱之间的条间力； 为第[i-l，j]±柱与第[i，j]±柱之间的条间力； G1+U为第[i，j]±柱与第[i+l，j]±柱之间的条间力； 胖1'^为第[1，引±柱的重力； α行方向的±柱间作用力与水平面夹角； 0列方向的±柱间作用力与水平面夹角； 表示第[ij]±柱底面上正应力的X方向余弦； /<:-/^表示第[1，引±柱底面上正应力的7方向余弦； 表示第[i，j]±柱底面上正应力的Z方向余弦； 表示第[ij]上柱底面上剪应力的X方向余弦； 吟表示第[i，j]上柱底面上剪应力的y方向余弦； 表示第[i，j]±柱底面上剪应力的Z方向余弦； 每个±柱的支持力和剪切力：式巧）、（6)中： C'为岩±体的有效粘聚力； 批为岩±体的有效内摩擦角； 辞6为抗剪强度随基质吸力变化的内摩擦角； Al'^为[i，j]±柱横截面的面积； Fs为边坡稳定系数； 所有力绕x、y、z轴方向的力矩分别为：式(7)、（8)、（9)中： AGi'j是平面ABB'A'和平面C孤'C'之间的力之差，ΔGi'j = GW'j-Gi'j; AQi'j是平面BCC'B'和平面A孤'A'之间的力之差，ΔQi'j = Qi'j+l-Qi'j; yAG' 'MG'·'是力Δ护域者Δ QIJ的X坐标，其值等于点Η的X坐标值； '是力Δ护域者Δ Q1J的y坐标，其值等于点Η的y坐标值； '.'心―'是力AG1'域者AQ1J的Z坐标，其值等于点Η的Z坐标值； 瑪是点Η的X坐标值;冷是点Η的y坐标值； 将式(5)、（6)代入上面式(7)、（8)、（9)，可求得边坡在不同时刻的稳定系数。
【文档编号】G01N33/24GK106093344SQ201610405311
【公开日】2016年11月9日
【申请日】2016年6月6日
【发明人】周小平, 陈洪, 程浩, 寿云东, 毕靖
【申请人】重庆大学