CN100573593C - 滑坡灾变过程时空预测的智能方法 - Google Patents
滑坡灾变过程时空预测的智能方法 Download PDFInfo
- Publication number
- CN100573593C CN100573593C CNB2008100692308A CN200810069230A CN100573593C CN 100573593 C CN100573593 C CN 100573593C CN B2008100692308 A CNB2008100692308 A CN B2008100692308A CN 200810069230 A CN200810069230 A CN 200810069230A CN 100573593 C CN100573593 C CN 100573593C
- Authority
- CN
- China
- Prior art keywords
- landslide
- sliding surface
- sliding
- gliding mass
- sin
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Abstract
本发明涉及一种滑坡灾变过程进行时空预测的智能方法。该方法通过计算机运行,包括如下步骤:1)选定待预测的具体滑坡,确定计算参数;2)建立滑坡智能预测的几何模型;3)对滑坡三维稳定性的智能预测;4)滑坡失稳状态的智能预测;5)滑坡活动强度的智能预测:计算滑坡失稳后的变形过程,输出任意时刻的变形图,直至滑坡结束。本发明的工作量较现有方法大大减少,提高了精度,同时简化了滑坡灾变过程时空预测的判断方法,使滑坡灾变过程时空预测更加方便和可靠,增强了滑坡灾变过程时空预测的可操作性。因此本发明是一种简单、直观、经济和高效的滑坡灾变过程时空预测方法。
Description
技术领域
本发明涉及一种滑坡时空预测方法,尤其是一种能对滑坡灾变过程进行时空预测的智能方法。更具体地说是一种岩土力学与现场监测相结合的智能预测方法。
背景技术
我国是滑坡地质灾害多发区,滑坡每年给我国造成巨大的经济损失和重大的人员伤亡。随着社会经济的不断发展,滑坡作为一种地质灾害,其严重的危害性已经给当地居民的生活造成了极大的影响。尽管对滑坡时空预测方法的研究工作已经取得了巨大成就,目前每年仍发生大量的滑坡地质灾害,表明现有的滑坡时空预测方法仍存在缺陷。
目前,滑坡时空预测还没有系统方法,工程实践中常常采用的方法是:首先根据二维极限平衡法判断滑坡稳定不稳定;然后根据短期的地表变形观测,凭经验判断滑坡大致失稳时间,最后决定是否加固。传统的滑坡时空预测方法存在着缺乏理论依据,以及不能预测滑体变形的不足,而且精度无法保证。另外,传统的滑坡时空预测方法是采用二维极限平衡法判别滑坡稳定性,实际上滑坡是一个三维问题,采用二维计算方法明显不能反映滑坡的真实应力状态;再者,传统方法通过对滑坡位移监测数据的简单处理来判断滑坡的失稳时间,这种方法没有考虑滑坡介质的流变特性,因此,物理意义不明确。
总之,传统的滑坡时空预测方法没有反映滑坡岩土体材料的非线性本质,没有反映滑坡的灾变过程,因此有必要发明一种既能反映滑坡岩土体材料非线性本质又能方便工程应用的滑坡时空预测系统方法。
发明内容
本发明的目的在于提供一种简单、直观、经济和高效的滑坡灾变过程时空预测的智能方法,该方法能判别滑坡的三维稳定性和滑坡的失稳时间,同时能够模拟滑坡失稳后的变形过程。
该方法通过计算机运行,包括如下步骤:
1)选定待预测的具体滑坡,确定计算参数;
2)建立滑坡智能预测的几何模型;
3)对滑坡三维稳定性的智能预测:对单滑面滑坡或多滑面滑坡的三维安全系数进行计算,判断安全系数是否大于1;当安全系数大于1时,再进入滑坡失稳状态判断;否则,程序结束;
4)滑坡失稳状态的智能预测:采用分叉集方程判断滑坡是否突滑失稳;如果结果是否,则继续采用饱和状态的参数,计算滑坡三维安全系数;如果结果是是,则预测出滑坡失稳时间,再进入滑坡活动强度预测;
5)滑坡活动强度的智能预测:计算滑坡失稳后的变形过程,输出任意时刻的变形图,直至滑坡结束。
具体地说,本发明是采用以下步骤对滑坡的灾变过程进行时空预测的:
1)选定具体待预测的滑坡,确定计算参数;如,确定三维滑动面的产状、孔隙水压力大小以及滑面上岩土材料在天然和饱和状态下的黏聚力c和内摩擦角对三维滑体进行离散化,根据滑面上剪切力方向与滑动面倾向相反近似确定每个离散块体的滑动方向;根据滑动面贯通率现场测量滑动面剪切段和蠕滑段的长度,确定滑面上材料西原体流变模型参数;
2)建立滑坡智能预测的几何模型:将三维滑体离散为竖直条柱,确定离散条柱的滑动方向,从而确定滑坡时间预测的二维剖面
3)对滑坡三维稳定性的智能预测:
对于单滑面滑坡,由极限平衡法求解二维滑面上正应力分布函数,根据三维滑体3个方向上力的平衡和绕3个坐标轴的力矩平衡得到安全系数Fs;
对于多滑面滑坡,采用极限分析上限法求解安全系数;
4)滑坡失稳时间的智能预测:
由滑坡位移实测数据拟合得到滑坡位移随时间变化多项式为:
u=u(0)+h1t+h2t2+h3t3+h4t4+h5t5
式中,u(0)为t=0时刻的初始水平位移,参数h1,h2,h3,h4,h5为时间变量t不同次方项前的常系数,可以通过拟合实测数据得到;
将步骤1)确定的滑面介质西原体流变参数输入尖点突变理论的分叉集方程4p3+27q2=0,得到滑坡突滑时间,同时根据能量转化原理,计算出滑体突滑初速度;
5)滑坡活动强度的智能预测:将滑体进行垂直条分离散化,得到每一条块的宽度、面积和弹性模量,计算滑落体活动过程任意时刻的变形、速度和滑动距离,直至滑坡结束。
上述预测方法的第2)步所建立的三维滑体几何模型中的竖直条柱底滑面上的滑动方向各不相同,每个条柱的滑动方向定义为条柱底滑面倾角的反向,其滑动方向单位矢量为:
式中,S为滑动方向单位矢量;θ为条柱底滑面相邻两条边的夹角;αx和αy分别为底滑面两边与坐标轴的夹角;α为滑面倾向与条柱底边的夹角。
上述预测方法的第3)步中三维滑体被垂直离散为m行和n列条柱,每个条柱按所在的行号i和列号j定义为(i,j)。
上述预测方法的第4)步中滑坡失稳时间预测以突变理论的分叉集方程作为判别条件,其判别条件为:方程4p3+27q2=0有正实根,式中p和q为控制变量。
上述预测方法第5)步中滑落体运动轨迹预测首先将滑体视为弹性滑块,滑块之间通过弹簧传递下滑力和抗滑力,同时滑块在下滑力和抗滑力等荷载作用下将发生弹性变形,然后根据能量守恒求解滑块的变形和速度。
本发明的工作量较现有方法大大减少,提高了精度,同时简化了滑坡灾变过程时空预测的判断方法,使滑坡灾变过程时空预测更加方便和可靠,增强了滑坡灾变过程时空预测的可操作性。因此本发明是一种简单、直观、经济和高效的滑坡灾变过程时空预测方法。本方法通过计算机实现,不但能判断滑坡的三维稳定性,而且可以根据突变理论的分叉集判断滑坡的突滑失稳时间以及模拟滑坡失稳后的变形过程。解决了现有的滑坡灾变全过程智能预测的难题,为滑坡灾变过程的跟踪与治理提供了一种简单实用的智能预测方法。
附图说明
图1为本发明实施例1的流程图;
图2为本发明实施例2的滑坡剖面图;
图3为本发明实施例3重庆某滑坡平面图;
图4(a)-图4(c)为某滑坡滑动过程变化图。
具体实施方式
实施例1:参见图1,一种滑坡灾变过程时空预测的智能方法,其具体实施步骤如下:
1)选定具体待预测的滑坡,确定计算参数;如,确定三维滑动面的产状、孔隙水压力大小以及滑面上岩土材料在天然和饱和状态下的黏聚力c和内摩擦角对三维滑体进行离散化,根据滑面上剪切力方向与滑动面倾向相反近似确定每个离散块体的滑动方向;根据滑动面贯通率现场测量滑动面剪切段和蠕滑段的长度,确定滑面上材料西原体流变模型参数;
2)建立滑坡智能预测的几何模型:将三维滑体离散为竖直的条柱,根据三维滑面的产状确定每个离散条柱的滑动方向,从而确定滑坡时间预测的二维剖面;
3)对滑坡三维稳定性的智能预测:对单滑面滑坡或多滑面滑坡的滑坡三维安全系数进行计算,判断安全系数是否大于1;当安全系数大于1时,再进入滑坡失稳状态判断;否则,程序结束;
对于单滑面滑坡,由极限平衡法求解二维滑面上正应力分布函数,根据三维滑体3个方向上力的平衡和绕3个坐标轴的力矩平衡得到安全系数Fs;
对于多滑面滑坡,采用极限分析上限法求解安全系数;
4)滑坡失稳时间的智能预测:判断滑坡是否突滑失稳;如果结果是否,则继续采用饱和状态的参数,计算滑坡三维安全系数;如果结果是是,则预测出滑坡失稳时间;再进入滑坡活动强度预测;
由滑坡位移实测数据拟合得到滑坡位移随时间变化多项式为:
u=u(0)+h1t+h2t2+h3t3+h4t4+h5t5
式中,u(0)为t=0时刻的初始水平位移,参数h1,h2,h3,h4,h5为时间变量t不同次方项前的常系数,可以通过拟合实测数据得到;
将步骤1)确定的滑面介质西原体流变参数输入滑坡失稳时间预测模型;根据尖点突变理论的分叉集方程4p3+27q2=0,得到滑坡突滑时间,同时根据能量转化原理,计算出滑体突滑初速度;
5)滑坡活动强度的智能预测:将滑体进行垂直条分离散化,得到每一条块的宽度、面积和弹性模量,计算滑落体活动过程任意时刻的变形、速度和滑动距离,直至滑坡结束。
其中第2)步所建立的三维滑体几何模型中的竖直条柱底滑面上的滑动方向各不相同,每个条柱的滑动方向定义为条柱底滑面倾角的反向,其滑动方向单位矢量为:
式中,S为滑动方向单位矢量;θ为条柱底滑面相邻两条边的夹角;αx和αy分别为底滑面两边与坐标轴的夹角;α为滑面倾向与条柱底边的夹角。
第3)步中三维滑体被垂直离散为m行和n列条柱,每个条柱按所在的行号i和列号j定义为(i,j)。
第4)步中滑坡失稳时间预测以突变理论的分叉集方程作为判别条件,其判别条件为:方程4p3+27q2=0是否有正实根,式中,p和q为控制变量。
第5)步中滑落体运动轨迹预测首先将滑体视为弹性滑块,滑块之间通过弹簧传递下滑力和抗滑力,同时滑块在下滑力和抗滑力等荷载作用下将发生弹性变形,然后根据能量守恒求解滑块的变形和速度。
实施例2:某滑坡三维安全系数确定
(1)主要计算参数
(2)计算流程
①步骤一:确定三维滑面形状及滑体几何尺寸
滑面函数s可以表示为
坡面函数z1为
根据式(1)和式(2)确定滑体的体积、厚度和滑面产状。
②步骤二:将三维滑体离散化
计算所用的离散模型如图2所示,滑体被离散为3017个垂直条柱。
③步骤三:根据滑面产状确定每个条柱的滑动方向
根据滑面上剪切力方向与滑动面倾向相反近似确定每个离散块体的滑动方向。
④步骤四:根据二维极限平衡法,假定滑面正应力分布函数
假设滑体达到极限平衡状态时滑面正应力分布为:
式中,Fs0为初始设定的安全系数;c为滑面岩土材料的黏聚力;为内摩擦角;λ1,λ2,λ3,λ4为待定参数;d1,d2为滑体最大宽度处y方向的坐标值;xa,xb分别为滑体出口和后缘x方向的坐标值;σ(x,y)为法向应力;q为孔隙水压力;ξ(x,y)为修正函数;Δ为滑面面积;σ0(x,y)为初始应力;αx为滑面倾角;w为滑体重力;x和y为条柱中心点坐标。
⑤步骤五:根据静力平衡方程求得安全的显式解
对于三维滑坡整体,需要考虑滑体的6个平衡条件,即x,y和z方向力的平衡,绕x轴、y轴和z轴的力矩平衡。
(3)计算结果分析
计算结果如表1所示,对比分析表明本发明方法的计算结果与二维极限平衡法和二维极限分析法计算结果接近。
表1不同计算方法求得例一的安全系数计算结果对比
可见,本发明所述方法的预测结果与国家规范的方法比较接近,但工作量却大大减少。
实施例3:重庆某滑坡灾变过程时空预测
(1)选定重庆某滑坡,确定计算参数
某滑坡位于渝黔高速公路重庆南坪段,该边坡于1997年12月到1998年3月开挖路基期间,由于受到大雨和人工扰动的影响逐步失稳形成了圆弧形拉裂缝。滑体坡面走向为SN向,滑坡在平面上呈不规则的马蹄形,横向宽约200~360m,纵向长约230m,滑坡后缘陡壁高10~20m,坡度70~80°。滑坡体厚度变化较大,呈周围浅中部深的规律,中部最厚达到40m,平均厚度约19.5m。滑坡区面积约7万m2,体积约140万m3,属大型滑坡,如图3所示。
(2)建立滑坡智能预测的几何模型并对滑坡三维稳定性的智能预测:
三维安全系数求解的具体步骤如下:
①步骤一:建立几何模型
将滑体按照折线形滑面的转折点离散为一系列垂直的条块,视每一条块为刚体,滑带和条块间错动部位被视为塑性体。对上下滑体中的离散块进行编号,同时找出代表主滑方向的中性面;该滑坡实际为三层滑面,由于上层滑体已经经过了一定的治理,上层滑体中的两个滑面较为接近且实测变形基本相近,因此,计算时仍按一层考虑。将上下滑体划分为在oxy平面投影为2m×2m的竖直条柱,条柱总数为8016。
②步骤二:求解机动许可的速度场
确定整个条块系统的速度场的计算过程主要包括以下几步:
a)假定中性面下滑体上的第c列条柱第一行的条块速度为V1,方向与底滑面的夹角为且平行于中性面,根据二维边坡极限分析法构造出中性面各条块机动许可的速度场;
b)按每一行第一列不存在下侧条块的边界条块速度绝对值等于相邻的左侧条块速度假定计算边界条块的速度场,同理可以计算每一行最后一列不存在下侧条块的边界条块;
c)根据相邻条块之间速度矢量闭合关系,从中性面向两侧逐步计算各条块的速度以及条块行与行、列与列之间的相对速度;
d)根据同一位置上下滑体条块的速度计算上下滑块的相对速度。
③步骤三:迭代求解安全系数
根据外部作用荷载和滑体自重所做的外功率与塑性变形区的内部能量耗损相等的平衡条件建立三维滑体虚功率方程为
式中,左边第一项代表列界面上的内能耗散,第二项代表行界面的内能耗散,第三项代表底滑面上的内能耗散,第四项代表上滑面上的内能耗散;右边第一项为滑体重力所做的功,第二项为外荷载所作的功。可以通过求解上式得到安全系数。
该滑坡为多滑面滑坡,为了考虑多滑面之间的相互影响,选用多滑面滑坡三维稳定性分析模型计算。主要计算参数如表2所示。
表2滑面计算参数
分别计算了该滑坡上层滑面和下层滑面的稳定性,得到了考虑滑面间相互作用后的安全系数,具体计算结果见表3。计算结果表明:下层滑面基本稳定,上层滑面不稳定或处于极限平衡状态。其中,上层滑面的计算结果比传递系数法的计算结果小,主要是因为考虑了滑面的相互作用。
表3不同方法安全系数计算结果对比
(3)滑坡失稳时间智能预测
具体实施步骤如下:
①步骤一:根据常规室内蠕变试验确定滑面介质的本构模型参数
计算主要参数为:剪切段剪切模量λ=1.64GPa,峰值位移uc=0.0109m,初始位移u0=0.0091m,软化系数φ=0.78,含水量ζ=0.05,滑体质量M=29328kg,剪切段倾角α1=35°,蠕滑段倾角α2=20°,后缘拉裂段长度hw=0.01m,滑体容重γ=24.92kN/m3,瞬时剪切模量G0=40kPa,长期剪切模量G∞=20kPa,黏弹性系数η1=4.05×1011Pa·s,黏塑性系数η2=2.887×1011Pa·s。
②步骤二:拟合位移监测数据得到位移时间关系
根据该滑坡3-3’剖面(如图3所示)的大地变形监测结果预测滑坡失稳时间,前面三维稳定性分析表明上层滑面处于极限平衡状态,下层滑面处于稳定状态,因此,本节只计算上层滑面的失稳时间。监测结果得到的位移时间关系为:
u=2×10-37t5-8×10-30t4-6×10-23t3+8×10-16t2-4×10-9t+0.0091 (5)
式中,位移单位为米(m),时间单位为秒(s)。
③步骤三:求解突滑时间及突滑速度
根据尖点突变理论的分叉集方程,得到了滑坡突滑时间的计算公式,同时根据能量转化原理,导出了滑体突滑初速度的计算公式。
按加固前的滑坡计算表明:当t=169.6天时,滑坡将要突滑失稳。通过本发明方法的计算表明:如果不采取支挡加固措施,该滑坡可能要突滑失稳造成严重损失。
(4)滑坡活动强度的智能预测
具体实施步骤如下:
①步骤一:将滑坡进行离散化并确定各条块的几何参数
滑落体失稳前如图4(a)所示,将整个滑落体被分为十个条块。其具体计算参数如下:弹性模量E0=15MPa,初速度vi,0=0,摩擦系数f=0.4,条块初始宽度bi,0=8.06m,计算时间步长Δt=0.02s。滑体被分为十个条块,每个条块的宽度均为8.06m,各条块的面积分别为表4所示。
表4各条块的面积
条块号 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
面积(m<sup>2</sup>) | 16.81 | 75.86 | 105.52 | 130.66 | 157.85 | 173.29 | 192.77 | 188.52 | 164.21 | 133.58 |
滑床各段的倾角为
②步骤二:根据弹簧-滑块模型求解各条块的速度
其详细计算方法为:
a)求解各条块的质量;
b)根据运动学公式,计算滑体突滑瞬间的条间力合突滑加速度;
c)选取一个微小的时间步增量Δt,根据滑体在滑动过程保持连续不分离的条件求得t+Δt时刻滑块的速度;
d)根据弹性力学公式求解条块变形,由后往前逐块求解t+Δt时刻滑块的条间力;
e)根据运动学公式求得条块在t+Δt时刻的加速度合摩擦阻力;
f)重复步骤b)~e)直到各条块的速度均为零,得到最后滑动时间和最大滑动距离。
③步骤三:根据每一时刻的速度和变形绘制任意时刻的滑坡整体剖面图。
计算结果表明:当t=8.06s时滑坡停止滑动,滑体前缘最大滑移距离为42.3m,滑坡灾变过程变形如图4(a)-4(c)所示。
可见,本发明所述方法能对多滑面滑坡三维安全系数进行求解,并能对滑坡灾变过程的活动强度定量预测。这些都是现有的滑坡时空预测方法所不能达到的。
Claims (2)
1、一种滑坡灾变过程时空预测的智能方法,其特征在于:
1)选定具体待预测的滑坡,确定三维滑动面的产状、孔隙水压力大小以及滑面上岩土材料在天然和饱和状态下的黏聚力c和内摩擦角;对三维滑体进行离散化,根据滑面上剪切力方向与滑动面倾向相反近似确定每个离散块体的滑动方向;根据滑动面贯通率现场测量滑动面剪切段和蠕滑段的长度,确定滑面上材料西原体流变模型参数;
2)建立滑坡智能预测的几何模型:将三维滑体离散为竖直条柱,确定离散条柱的滑动方向,从而确定滑坡时间预测的二维剖面;
3)对滑坡三维稳定性的智能预测:
对于单滑面滑坡,由极限平衡法求解二维滑面上正应力分布函数,根据三维滑体3个方向上力的平衡和绕3个坐标轴的力矩平衡得到安全系数Fs;
对于多滑面滑坡,采用极限分析上限法求解安全系数;
4)滑坡失稳时间的智能预测:
由滑坡位移实测数据拟合得到滑坡位移随时间变化多项式为:
u=u(0)+h1t+h2t2+h3t3+h4t4+h5t5
式中,u(0)为t=0时刻的初始水平位移,参数h1,h2,h3,h4,h5为时间变量t不同次方项前的常系数,通过拟合实测数据得到;
将步骤1)确定的滑面介质西原体流变参数输入尖点突变理论的分叉集方程4p3+27q2=0,得到滑坡突滑时间,同时根据能量转化原理,计算出滑体突滑初速度;
5)滑坡活动强度的智能预测:将滑体进行垂直条分离散化,得到每一条块的宽度、面积和弹性模量,计算滑落体活动过程任意时刻的变形、速度和滑动距离,直至滑坡结束。
2、如权利要求1所述的滑坡灾变过程时空预测的智能方法,其特征在于:
所述预测方法的第2)步所建立的三维滑体几何模型中的竖直条柱底滑面上的滑动方向各不相同,每个条柱的滑动方向定义为条柱底滑面倾角的反向,其滑动方向单位矢量为:
式中,S为滑动方向单位矢量;θ为条柱底滑面相邻两条边的夹角;αx和αy分别为底滑面两边与坐标轴的夹角;α为滑面倾向与条柱底边的夹角;
所述预测方法的第3)步中三维滑体被垂直离散为m行和n列条柱,每个条柱按所在的行号i和列号j定义为(i,j);
所述预测方法的第4)步中滑坡失稳时间预测以突变理论的分叉集方程作为判别条件,其判别条件为:方程4p3+27q2=0有正实根,式中p和q为控制变量;
所述预测方法第5)步中滑落体运动轨迹预测首先将滑体视为弹性滑块,滑块之间通过弹簧传递下滑力和抗滑力,同时滑块在下滑力和抗滑力等荷载作用下将发生弹性变形,然后根据能量守恒求解滑块的变形和速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2008100692308A CN100573593C (zh) | 2008-01-10 | 2008-01-10 | 滑坡灾变过程时空预测的智能方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2008100692308A CN100573593C (zh) | 2008-01-10 | 2008-01-10 | 滑坡灾变过程时空预测的智能方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101216955A CN101216955A (zh) | 2008-07-09 |
CN100573593C true CN100573593C (zh) | 2009-12-23 |
Family
ID=39623383
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2008100692308A Expired - Fee Related CN100573593C (zh) | 2008-01-10 | 2008-01-10 | 滑坡灾变过程时空预测的智能方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100573593C (zh) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102279134A (zh) * | 2011-04-06 | 2011-12-14 | 中国科学院武汉岩土力学研究所 | 基于现场直剪试验的滑坡滑动面力学参数的计算方法 |
CN102619209A (zh) * | 2012-04-25 | 2012-08-01 | 重庆大学 | 基于监测位移的岩土工程灾变动态实时智能预警方法 |
CN103135128B (zh) * | 2013-01-31 | 2015-07-01 | 重庆大学 | 一种地震荷载作用下三维边坡稳定性预测方法 |
CN103163563B (zh) * | 2013-01-31 | 2015-08-19 | 重庆大学 | 一种三维边坡稳定性预测方法 |
CN103543098B (zh) * | 2013-10-30 | 2015-10-07 | 西南交通大学 | 用于高速远程滑坡裹气流变特性研究的剪切流变装置 |
CN104063616B (zh) * | 2014-07-04 | 2017-05-10 | 重庆大学 | 一种滑坡失稳运动作用下建筑物破坏损失程度评价方法 |
CN104200100B (zh) * | 2014-09-01 | 2017-02-15 | 重庆大学 | 一种基于滑动位移分析的三维边坡稳定性预测方法 |
CN104182646B (zh) * | 2014-09-01 | 2017-02-15 | 重庆大学 | 一种基于滑动位移分析的三维地震边坡发生滑坡的时间预测方法 |
CN106503434A (zh) * | 2016-10-18 | 2017-03-15 | 河南理工大学 | 一种滑坡预报方法及装置 |
CN106776472B (zh) * | 2016-11-29 | 2019-04-12 | 曹杰 | 一种月尖突变可视化呈现方法及系统 |
CN106777585B (zh) * | 2016-12-01 | 2019-10-11 | 桂林理工大学 | 一种区域浅表滑坡时空分异规律的esda分析法 |
CN106710148A (zh) * | 2016-12-28 | 2017-05-24 | 中国科学院深圳先进技术研究院 | 地质灾害预测方法及系统 |
CN107688895B (zh) * | 2017-08-09 | 2018-08-07 | 中国地质环境监测院 | 一种冻融型滑坡安全性分析及运动距离测算方法 |
CN108418873A (zh) * | 2018-02-12 | 2018-08-17 | 福建网即通网络科技有限公司 | 一种智能农业管控方法和系统 |
CN110263968B (zh) * | 2019-05-07 | 2021-04-30 | 深圳大学 | 一种新的时间序列的预测方法、装置、电子设备及计算机可读存储介质 |
CN110189505B (zh) * | 2019-06-14 | 2021-03-23 | 成都理工大学 | 一种岩层层面滑动面岩质滑坡的预警方法及其应用 |
CN112183874B (zh) * | 2020-10-12 | 2021-06-18 | 兰州大学 | 一种滑坡滑距预测方法 |
CN112733412B (zh) * | 2020-12-28 | 2023-03-24 | 河海大学 | 一种水动力作用滑坡运动机制研究的速度等效表征方法 |
CN113065702B (zh) * | 2021-03-31 | 2023-08-08 | 中国地质大学(武汉) | 基于st-seep分段法和时空arma模型的滑坡位移多线性预测方法 |
CN114491737B (zh) * | 2021-12-30 | 2023-07-14 | 北京市政路桥股份有限公司 | 一种路堑砌体挡墙的稳定性力学分析方法 |
-
2008
- 2008-01-10 CN CNB2008100692308A patent/CN100573593C/zh not_active Expired - Fee Related
Non-Patent Citations (8)
Title |
---|
A three-dimensional slope stability analysismethodusingtheupper bound theorem - Part I: theory andmethods. Chen Z., Wang X ET AL.International Journal of Rock Mechanics and Mining Sciences,Vol.38 No.3. 2001 |
A three-dimensional slope stability analysismethodusingtheupper bound theorem - Part I: theory andmethods. Chen Z., Wang X ET AL.International Journal of Rock Mechanics and Mining Sciences,Vol.38 No.3. 2001 * |
二维高速滑坡力学模型. 肖盛燮,周小平,杨海清,陈山林.岩石力学与工程学报,第25卷第3期. 2006 |
二维高速滑坡力学模型. 肖盛燮,周小平,杨海清,陈山林.岩石力学与工程学报,第25卷第3期. 2006 * |
滑坡灾害的时间预测预报研究. 李秀珍.成都理工大学硕士学位论文. 2004 |
滑坡灾害的时间预测预报研究. 李秀珍.成都理工大学硕士学位论文. 2004 * |
非规则滑坡三维稳定性分析. 芮勇勤,吕强. 中外公路,第22卷第5期. 2002 |
非规则滑坡三维稳定性分析. 芮勇勤,吕强. 中外公路,第22卷第5期. 2002 * |
Also Published As
Publication number | Publication date |
---|---|
CN101216955A (zh) | 2008-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100573593C (zh) | 滑坡灾变过程时空预测的智能方法 | |
CN106557608B (zh) | 一种基于混合数值离散的节理岩体的极限分析上限法 | |
CN103485353B (zh) | 基于全局最优化的边坡稳定性分析条分法 | |
Niu et al. | Experimental research on the aerodynamic characteristics of a high-speed train under different turbulence conditions | |
CN101344389B (zh) | 一种利用神经网络预测隧道围岩位移的方法 | |
CN108627137B (zh) | 一种滑坡变形预测计算方法 | |
Zhou et al. | An enhanced equivalent continuum model for layered rock mass incorporating bedding structure and stress dependence | |
CN107330146A (zh) | 一种同时考虑平动和转动效应的节理岩质边坡极限承载力分析上限法 | |
CN103135128A (zh) | 一种地震荷载作用下三维边坡稳定性预测方法 | |
Safaei Pirooz et al. | Comparison of speed-up over hills derived from wind-tunnel experiments, wind-loading standards, and numerical modelling | |
CN103150421A (zh) | 利用位移监测同时确定抗滑桩桩位和临界深度的方法 | |
CN110348125B (zh) | 一种块体稳定性的可视化判别方法 | |
CN111931369B (zh) | 降雨型滑坡稳定性分析及运动距离测算方法、设备及介质 | |
CN110210776A (zh) | 一种基于监控量测的隧道施工动态风险评估方法 | |
CN106250590A (zh) | 一种基于cfd数值模拟的高空露天平台行人风环境评估方法 | |
CN103163563A (zh) | 一种三维边坡稳定性预测方法 | |
CN104200100A (zh) | 一种基于滑动位移分析的三维边坡稳定性预测方法 | |
CN113408184B (zh) | 基于遗传算法和离散元分析法的锚固边坡安全性评价方法 | |
CN103266617A (zh) | 岩质边坡楔形体最优锚固角的计算方法 | |
KR102024195B1 (ko) | 궤도-지반 상호작용을 고려한 야지궤도차량의 구동성능 평가방법 및 평가시스템 | |
CN104778336A (zh) | 一种飞机系留载荷计算方法 | |
Mansour | Approximate probabilistic method of calculating ship longitudinal strength | |
Ren et al. | Comprehensive safety evaluation method of surrounding rock during underground cavern construction | |
CN108170898A (zh) | 一种节理岩质边坡可靠度分析下限法 | |
CN115455759A (zh) | 一种含软弱夹层矿山边坡的动态监测预警方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091223 Termination date: 20130110 |
|
CF01 | Termination of patent right due to non-payment of annual fee |