一种时域内降雨作用下的分段顺坡滑坡预测方法技术

本发明专利技术公开了一种时域内降雨作用下的分段顺坡滑坡预测方法，涉及到分段顺坡滑坡预测方法技术领域，通过将基于物理概念的边坡模型与土力学，水力学相结合对边坡的稳定性进行分析，进而得到边坡预警方法，降雨对于边坡的作用通常是雨水大量的下渗导致土的含水率增加，土的抗剪强度降低，因此下部土体发生了滑移导致上部土体继而失稳进行滑动，本发明专利技术的有利效果是将滑坡确定性模型与基于物理概念的边坡模型进行了有效的集成，分析在降雨作用下，降雨时间和土体含水率的关系，通过边坡模型的边界条件得到体积含水率与浸润锋深度之间的关系，在此基础上结合土力学的计算得出土块的土压力。压力。压力。

【技术实现步骤摘要】
一种时域内降雨作用下的分段顺坡滑坡预测方法


[0001]本专利技术涉及分段顺坡滑坡预测方法
，特别涉及一种时域内降雨作用下的分段顺坡滑坡预测方法。

技术介绍

[0002]我国是地质灾害频发的国家，其中滑坡作为一种地质灾害，严重的影响人们的生产生活。预判滑动位置，滑坡方式一直是工程界研究的重点之一，其中理论分析虽然越来越丰富但是由于边坡滑动突发性的特征，其运动规律的不确定性导致有关滑坡方面的检测的精度不是很高。降雨对于滑坡的影响十分明显，雨水大量的下渗导致土的含水率增加，雨水在边坡坡角处积聚，土的抗剪强度指标迅速下降，因此边坡下部土体发生了滑移进而牵引上部土体失稳进行滑动。根据下部土体的滑动的自然规律得出滑块的形状和坡脚处的坡角关系，找出不同边坡体最危险滑块的滑角，可以对不同路基边坡的滑动加以预测。
[0003]滑坡稳定性的预测方法大致分为两种：确定性模型以及非确定性模型。确定性模型是将模型通过力学极限平衡的原理对不同类型的滑坡进行稳定性的分析。非确定模型通过数理统计的方法对数据进行分析。基于选用模型的不同，滑坡预测可分为基于经验概念边坡模型和基于物理概念边坡模型这两大类。基于经验概念的边坡模型相对简单，发展较早。基于物理概念的边坡模型模拟的信息丰富，但是求解复杂对计算的要求高，以此基础建立的滑坡预测模型较晚。
[0004]综合国内外研究现状可以看出,目前关于边坡的评价与预测预报理论方法中,包含有极限平衡分析法和位移时序预测法等。极限平衡理论是一种经典的确定性分析方法,建立在失稳机制和受力条件明确的基础上。
[0005]因此，专利技术一种时域内降雨作用下的分段顺坡滑坡预测方法来解决上述问题很有必要。

技术实现思路

[0006]本专利技术的目的在于提供一种时域内降雨作用下的分段顺坡滑坡预测方法，以解决上述
技术介绍
中提出的问题。
[0007]为实现上述目的，本专利技术提供如下技术方案：一种时域内降雨作用下的分段顺坡滑坡预测方法，所述预测方法包括以下步骤S1、根据太沙基一维固结理论与达西定律确定浸润线位置；
[0008]通过边界条件得到降雨时间和雨水入渗深度及土体含水率的关系：
[0009][0010]利用经典的Green
‑
Amp模型得到体含水率与时间的关系：
[0011][0012]S2、计算出土块的土压力，找到剩余土压力关于β角的函数关系式，滑角β为滑动面与坡面的夹角，基于确定性模型的特点改变滑角β以1
°
为步长进行搜索最危险的滑动面；
[0013]T3＝f(β)＝A sin(a
′‑
β)+ψ2T1‑
Bcos(a
′‑
β)
‑
C
[0014]式中：θ为体含水率，θ
I
为初始体积含水率；θ
s
为饱和体积含水率，h
m
为吸力水头，D为扩散系数，ψ
i
为传导系数，t为降雨时间，h
m
为浸润锋吸力水头。
[0015]本专利技术的技术效果和优点：
[0016]1、以往的边坡稳定性计算是在某一个特定时刻对边坡的稳定性安全性的评价，本专利技术在时间发展中的一个时域内来对边坡的稳定性安全性进行评价；
[0017]2、所用的边坡稳定分析模型相较于传统边的坡模型更为简便，本专利技术计算摒弃了复杂的计算过程使得更加的简单直观。
附图说明
[0018]图1为本专利技术三段式滑坡示意图。
[0019]图2为本专利技术降雨时间与入渗深度关系图。
[0020]图3为本专利技术T3与降雨时间关系图。
[0021]图4为本专利技术土块间剩余土压力的传导示意图。
[0022]图5为本专利技术滑坡预测方法流程图。
具体实施方式
[0023]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚；完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0024]本专利技术提供了如图1
‑
图5所示的一种时域内降雨作用下的分段顺坡滑坡预测方法，
[0025]第一步、根据太沙基一维固结理论以及达西定律可得出降雨时间和雨水的入渗深度的关系，搜索到浸润线的位置。
[0026][0027]C
v
为固结系数
[0028][0029]对于一维固结方程即渗流仅在一个方向上(z)
[0030]若将C
v
替换为(为水分扩散系数常量)
[0031][0032]式(4)为一维入渗Richards[2]方程
[0033]分离变量
[0060]本专利技术以典型三段式顺坡滑动坡体为例进行计算下表为边坡的工程概况：
[0061]表1
[0062][0063]本专利技术这种边坡的稳定安全性的分析方法，包括以下步骤：
[0064](1)首先通过太沙基一维固结理论对坡体进行分析，找到土体含水率与雨水入渗深度的关系来确定浸润线的位置
[0065][0066]依据现有的模型分析出边界条件
[0067][0068]将边界条件带入到土体含水率与雨水入渗深度的公式中求解。
[0069]利用经典的Green
‑
Ampt积水入渗模型令地表处位置水头ξ＝0，浸润锋在吸力水头H
m
和重力水头ξ
w
共同作用下向下推进，则浸润锋上的总水头为H
t
＝h
m
‑
ξ
w
。因此，根据达西定律，可得到以下公式：
[0070][0071]由公式可以得到降雨时间与雨水入渗深度的关系,确定浸润线的位置。
[0072]当降雨时间为30min、60min、90min、120min、150min时将降雨时间代入到式(12)中可计算出雨水入渗的深度。
[0073]表2
[0074][0075]入渗深度随着降雨时间的增长呈现出入渗越深的规律，如图2所示。
[0076]第二步、本专利技术是一种多段式顺坡滑动的边坡模型，以滑角处β的变化不断搜索剩余的滑动推力等于零的临界值进而找到最危险的滑动面。
[0077]根据《公路路基设计规范》(JTG D30—2015)[4]计算滑坡推力时应考虑：滑体内地下水产生的静水压力
[0078]滑坡剩余下滑力可采用传递系数法按(13)式计算条块作用力系见图1
[0079][0080][0081][0082]式中γ
sat
为土的饱和重度；为土块的浸润锋深度；为土块的平均宽度；T
i
T
i
‑1为第i和i
‑
1滑块剩余下滑力(KN/m)；F
s
为稳定系数；W
i...

【技术保护点】

【技术特征摘要】
1.一种时域内降雨作用下的分段顺坡滑坡预测方法，其特征在于，所述预测方法包括以下步骤：S1、根据太沙基一维固结理论与达西定律确定浸润线位置；通过边界条件得到降雨时间和雨水入渗深度及土体含水率的关系：利用经典的Green
‑
Amp模型得到体含水率与时间的关系：S2、计算出土块的土压力，找到剩余土压力关于β角的函数关系式，滑角β为滑动面与坡面的夹角，基于确定性模型的特点改变滑角β以1
°
为步长进行搜索最危险的滑动面；
·<...

【专利技术属性】
技术研发人员：马江平，尹文峰，李霞，宋颖洁，万伟，崔宏环，谢直树，蒋明学，朱晨，王跃庚，李续靖，王雅超，訾亚丹，
申请(专利权)人：河北建筑工程学院，
类型：发明
国别省市：