滑坡计算参数反演分析

K＝1.05
C
20
21.8
23.6
25.4
27.2
29
30.8
32.6
34.4
36.2
φ(º)
16.69
16.43
16.19
15.94
15.69
15.43
15.19
14.92
14.67
14.42
Si
0.930
0.9பைடு நூலகம்4
0.896
0.889
0.890
1.087
1.485
1.903
2.300
2.707
根据反算参数的比较值可知，在稳定性系数K=1.0时，滑带土的粘聚力C和内摩擦角φ分别为27.2和14.80；K=1.05时，粘聚力C和内摩擦角φ分别为25.4和15.94。
Table 1 The shear strengthparametersstatistic table of sliding zone of T2b1and T2b3strata
工况
参数
指标
分布区间
样本个数
分布形式
分布函数
天然
（Kpa）
（10.0，43.1）
65
正态分布
（°）
（6.0，23.0）
92
4反演计算模型的建立
4.1稳定性计算方法
常用于计算滑坡稳定系数的方法有Bishop法、Sarma法、剩余推力法。本文采用剩余推力法（residual thrust method，RTM）建立反演分析所需的极限平衡方程。条块间的作用力，亦即上一条块的剩余下滑力，其方向与上一条块的底面平行，且条块之间传压不传拉；滑坡整体的剩余下滑力，即最后一个条块的剩余下滑力为0。
[13]Maio,C.D.Theinfluenceofporefluidcompositionontheresidualshearstrengthofsomenaturalclayeysoil.Proceedingsof7thInter.CongressofLandslides.1996
[14]Popesku,M.E.Backanalysisofslopefailures—Apossibilityorachallenge.proceddingsof7thInternationalIAEGcongress.1994
3滑带土抗剪强度参数统计
对三峡库区二期崩塌滑坡治理工程和三期规前勘（调）察中的崩塌滑坡点的勘察试验资料进行分类统计，得到本区滑带土抗剪强度参数值，可以用于验证和优化反演得到的参数。经统计得到适合该滑坡的抗剪强度参数分布函数如表1，图2是滑带抗剪强度参数统计直方图。
表1 T2b1和T2b3滑带土的抗剪强度参数统计表
目前，滑带土抗剪强度参数反演分析的方法分为单参数反演和双参数反演两种。前者假定一个参数已知的前提下，反算另外一个参数，通常选择对滑坡稳定性影响较敏感的作为未知参数。后者在反演中有两个未知的参数，通常选择两个距主滑动面等距的剖面建立极限平衡方程求解。本文以三峡库区太山庙滑坡为例，在C、 值未知的情况下，综合采用经验类比和反演分析方法确定滑带土的抗剪强度参数，分析时兼顾了滑坡的区域相似性和个体特性，所得到的结果更为准确、可靠。
K＝1.0
C(Kpa)
20
21.8
23.6
25.4
27.2
29
30.8
32.6
34.4
36.2
φ(º)
15.80
15.55
15.30
15.04
14.80
14.54
14.28
14.03
13.77
13.51
Si
1.614
1.606
1.608
1.609
1.601
1.798
2.206
2.614
3.021
3.429
2滑坡概况
欧家湾滑坡位于奉节县白帝镇坪上村2、3组，长江支流石马河左岸，属于三峡库区三期专业监测崩塌滑坡灾害点。滑坡无详细的勘察资料，仅在监测设计阶段做了地面调查。
滑坡自然坡角约25~40°，滑坡前缘临近石马河处零星分布石马河一级阶地，滑坡区属低山丘陵剥蚀地貌。滑坡体的主滑方向为5°，平面形态呈箕形，由后缘向前缘逐渐变宽，滑坡东西宽约350~400m，南北向主轴长约420m。后缘高程约325m，前缘高程约170m，左侧以山脊为界，右侧以冲沟为界，总变形规模约507×104m3。滑体主要由第四系碎块石土夹粘性土组成，滑床为巴东组第三段（T2b3）的泥灰岩，岩层产状为280°∠3°，为斜交坡，图1是滑坡的工程地质剖面图。
图1欧家湾庙滑坡工程地质剖面图
Fig.1 The engineering geological profile of Oujiawan landslide
滑坡为老滑坡，滑坡区经过过去的剧烈滑动后，在改变了当时的地形地貌后形成了现今的老滑坡体地形。经对现场的调查踏勘发现，滑坡体上树木歪斜，现仍有滑移变形产生。在滑坡中部多户民房附近，近年每逢雨季都有蠕动滑移。从地表调查和发展趋势上看，目前该滑坡处于不稳定状态。
[15]Saito,M.ReverseCalculationmethodtoobtainCand&aslipsurface,ProceedingsIntersymposium.Landslide.NewDelhi.1980
[2]彭亚明，彭军还，张彬等．三峡库区某滑坡抗剪强度参数的反演分析［J］．桂林工学院学报，2003，23（3）：279～283．
[3]富凤丽，佴磊．中里滑坡反分析及强度系数取值研究［J］．长春科技大学学报，2000，30（2）：165～169．
[4]林鲁生，蒋刚，白世伟等．土体抗剪强度参数取值的统计分析法［J］．岩土力学，2002，23（5）：570～574．
将反算得到的一系列抗剪强度参数值与统计的结果进行对比分析，选择最优的参数作为滑坡稳定性计算和设计的强度参数。用于比较的计算公式如下：
得到的一系列比较函数值Si中最小的则为最优的滑带土抗剪强度参数，表3为经反演得到的一系列参数和相应的比较函数值。
表3反演参数和比较函数值
Table 3 Back analysis parameters and comparison function values
图5滑坡敏感性分析
Fig.5 Sensitivity analysis of landslide
4.4滑带土抗剪强度参数反分析及优化取值
上面已经确定了该滑坡的稳定性评估指标，即已知滑坡的稳定性系数，可以通过反算得到滑带土的抗剪强度参数C、 。选择滑坡主滑动剖面作为稳定性的计算剖面，由于存在两个未知参数，首先，确定C值验算范围，并在该范围内选择一系列验算参数值 ，再通过反算得到相应的一系列 值，即 。
对天然状态下滑带土抗剪强度参数的试验值由小到大进行排序，其饱和强度参数应与天然状态的一一对应。对排序后的参数值进行分段，对各段内的天然强度参数和相应的饱和强度参数求算术平均值，得到相应的关系曲线如图3。
图3滑带土天然和饱和抗剪强度参数关系曲线
Fig.3 The relation curve ofnatural shearstrengthsparameters and saturation shearstrengthsparameters of sliding zone
由式（1）～（3）消去 ， ，得到 与 的关系：
（4）
如果 ，则令
其中
（5）
（6）
（7）
式中 为推力传递系数
经试算迭代，当最后一条块剩余推力为0时，所对应的K即为滑坡的稳定性系数。
4.2确定计算状态和稳定性评估指标
通常将反分析的状态称为临界状态，临界状态是指在确定工况的评估指标下的边坡即时状态，包括坡面形态、地下水位、滑带赋存条件和外荷载等因素。再确定了计算状态后应该使边坡的临界状态各因素符合实际情况。
稳定评估指标主要依据边坡宏观变形状况来确定,，建议参考表2并结合实际情况选定。稳定评估指标具有一定的先验性，必须考虑边坡不同发育阶段的变形性质并详细查勘边坡前、后缘变形量和地形变化后才能做出正确选择。
表2滑坡不同发展阶段的稳定系数
Table 2 Stability coefficient fordifferentstages of landslide
计算模型如图图4：
对于第i（i=1，2，…，n）个条块，沿平行及垂直条块底面方向建立局部坐标系，由力平衡方程 得到：
图4条块受力分析示意图
Fig.4 Schematic diagram of forces on slices
（1）
（2）
其中 （3）
式中Ci、 、 分别为第I条块的粘聚力、内摩擦角、底面长度
4.3敏感性分析
选取粘聚力和内摩擦角的均值作为基准参数，当一个参数在某一范围内变化时另外一个参数不变。根据统计结果确定参数变化范围，参数变化区间为 。敏感性分析结果如图5，内摩擦角是影响滑坡稳定性的敏感因素，根据敏感性分析图的曲线斜率，计算粘聚力和内摩擦角在同方差范围内对滑坡稳定性的影响程度比值 ＝1：3。
[5]龚玉锋，周创兵，梁轶等．参数反演在岩质高边坡变形与稳定分析中的应用［J］．岩石力学，2003，24（5）：809～811．
[6]吴刚，夏艳华，陈静曦等．可靠性理论在边坡反分析中的运用［J］．岩石力学，2005，27（5）：390～393．
[7]徐汉斌，王军．反算法中滑坡稳定系数的取值问题［J］．四川地质学报，1999，19（1）：86～89．
[11]徐青，陈士军，陈胜宏．滑坡稳定分析剩余推力法的改进研究［J］．岩土力学，2005，26（3）：465～470．
[12]Duncan,J.M.,Stark,T.D.Soilstrengthfrombackanalysisofslopefailures.Proc.ASCEGeotechnicalConference.SlopeandEmbankments.Berkely.1992
5结论
（1）根据欧家湾滑坡宏观变形特征，确定了天然状态下滑坡稳定性评估指标，稳定性系数为1.0～1.05。
（2）对三峡库区已有勘察资料的滑坡进行分类统计，得到了适合本滑坡滑带土抗剪强度参数的分布函数，并通过分段求平均值的方法得到天然状态和饱和状态下滑带土抗剪强度参数的相关函数。
（3）根据敏感性分析，认为内摩擦角是影响滑坡稳定性的敏感因子，在同方差变化范围内，粘聚力和内摩擦角对稳定性的影响程度比值为1：3。
正态分布
饱和
（Kpa）
（8.7，31.4）
70
对数正态分布
（°）
（5.3，19.8）
89
正态分布
图2 T2b1和T2b3地层滑带土抗剪强度参数统计直方图
Fig.2 The statistic histograms of shear strengthparametersof sliding zone of T2b1and T2b3strata
滑坡计算参数反演分析的优化算法
1引言
在滑坡稳定性计算和工程设计中，滑带土的粘聚力（C）和内摩擦角（ ）取值正确与否至关重要。目前确定滑带土抗剪强度参数（C、 ）值的方法有试验、工程类比和反演分析3种。滑带土剪切试验分为现场或室内两种，受试样和试验条件的限制，滑带土试验数据通常很离散，需要进行分析计算来确定。工程类比法在确定滑带土的抗剪强度参数时具有很强的主观性，在确定类比指标时又受到类比滑坡客观条件的限制。反演分析是确定滑带土抗剪强度参数的一种有效的方法，根据滑坡的宏观变形状况假设滑坡的稳定性系数，再反算滑带土抗剪强度参数。反算是滑坡稳定性计算的逆过程，得到的参数更符合滑坡的变形情况，参数可以作为试验数据选取的参考，若没有试验数据时，可以直接作为稳定性计算、工程设计的参数。
（4）在滑带土抗剪强度参数反演分析的基础上，将得到的一系列参数值与统计资料进行对比分析得到了适合本滑坡的最优参数值，优化过程兼顾了滑坡的区域相似性和个体的特性。
参考文献：
[1]高德军，徐卫亚，郭其达．长江三峡大石板滑坡计算参数反分析［J］．河海大学学报（自然科学版），2006，34（1）：74～78．
[8]周平根.滑带土强度参数的估算方法[J].水文地质工程地质, 1998, (6): 30－32.
[9]王乐华，杨学堂，李建林．水布垭瓦屋场滑坡参数反分析研究［J］．三峡大学学报（自然科学版），2005，27（5）：390～393．
[10]李端有，甘孝清．滑坡体力学参数反分析研究［J］．长江科学院院报，2005，22（6）：44～48．
发展阶段
变形性质
稳定性系数K
变形状态
局部变形
蠕滑
1.05～1.00
前缘或后缘变形微弱，地表出现未连通微裂缝
整体变形
微滑
1.00～0.90
局部坡面变形异常，陡坎出出现小型局部滑崩，裂缝发育，逐渐连通
剧滑
<0.90
坡面出现鼓丘、挤压变形和较大的大裂隙
稳定固结
固结
>1.05
根据前面对欧家湾滑坡宏观变形状况的描述，判断滑坡目前处于蠕滑的局部变形阶段，确定稳定性系数为1.00～1.05。