滑坡计算参数反演分析

稳定评估指标主要依据边坡宏观变形状况来确定,，建议参考表2并结合实际情况选定。稳定评估指标具有一定的先验性，必须考虑边坡不同发育阶段的变形性质并详细查勘边坡前、后缘变形量和地形变化后才能做出正确选择。
表2 滑坡不同发展阶段的稳定系数
Table 2 Stability coefficient fordifferentstages of landslide
对天然状态下滑带土抗剪强度参数的试验值由小到大进行排序，其饱和强度参数应与天然状态的一一对应。对排序后的参数值进行分段，对各段内的天然强度参数和相应的饱和强度参数求算术平均值，得到相应的关系曲线如图3。
图3 滑带土天然和饱和抗剪强度参数关系曲线
Fig.3 The relation curve ofnatural shearstrengthsparameters and saturation shearstrengthsparameters of sliding zone
3 滑带土抗剪强度参数统计
对三峡库区二期崩塌滑坡治理工程和三期规前勘（调）察中的崩塌滑坡点的勘察试验资料进行分类统计，得到本区滑带土抗剪强度参数值，可以用于验证和优化反演得到的参数。经统计得到适合该滑坡的抗剪强度参数分布函数如表1，图2是滑带抗剪强度参数统计直方图。
表1 T2b1和T2b3滑带土的抗剪强度参数统计表
计算模型如图图4：
对于第i（i=1，2，…，n）个条块，沿平行及垂直条块底面方向建立局部坐标系，由力平衡方程 得到：
图4 条块受力分析示意图
Fig.4 Schematic diagram of forces on slices
（1）
（2）
其中 （3）
式中Ci、 、 分别为第I条块的粘聚力、内摩擦角、底长度
由式（1）～（3）消去 ， ，得到 与 的关系：
（4）
如果 ，则令
其中
（5）
（6）
（7）
式中 为推力传递系数
经试算迭代，当最后一条块剩余推力为0时，所对应的K即为滑坡的稳定性系数。
4.2 确定计算状态和稳定性评估指标
通常将反分析的状态称为临界状态，临界状态是指在确定工况的评估指标下的边坡即时状态，包括坡面形态、地下水位、滑带赋存条件和外荷载等因素。再确定了计算状态后应该使边坡的临界状态各因素符合实际情况。
2 滑坡概况
欧家湾滑坡位于奉节县白帝镇坪上村2、3组，长江支流石马河左岸，属于三峡库区三期专业监测崩塌滑坡灾害点。滑坡无详细的勘察资料，仅在监测设计阶段做了地面调查。
滑坡自然坡角约25~40°，滑坡前缘临近石马河处零星分布石马河一级阶地，滑坡区属低山丘陵剥蚀地貌。滑坡体的主滑方向为5°，平面形态呈箕形，由后缘向前缘逐渐变宽，滑坡东西宽约350~400m，南北向主轴长约420m。后缘高程约325m，前缘高程约170m，左侧以山脊为界，右侧以冲沟为界，总变形规模约507×104m3。滑体主要由第四系碎块石土夹粘性土组成，滑床为巴东组第三段（T2b3）的泥灰岩，岩层产状为280°∠3°，为斜交坡，图1是滑坡的工程地质剖面图。
目前，滑带土抗剪强度参数反演分析的方法分为单参数反演和双参数反演两种。前者假定一个参数已知的前提下，反算另外一个参数，通常选择对滑坡稳定性影响较敏感的作为未知参数。后者在反演中有两个未知的参数，通常选择两个距主滑动面等距的剖面建立极限平衡方程求解。本文以三峡库区太山庙滑坡为例，在C、 值未知的情况下，综合采用经验类比和反演分析方法确定滑带土的抗剪强度参数，分析时兼顾了滑坡的区域相似性和个体特性，所得到的结果更为准确、可靠。
正态分布
饱和
（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
4 反演计算模型的建立
4.1稳定性计算方法
常用于计算滑坡稳定系数的方法有Bishop法、Sarma法、剩余推力法。本文采用剩余推力法（residual thrust method，RTM）建立反演分析所需的极限平衡方程。条块间的作用力，亦即上一条块的剩余下滑力，其方向与上一条块的底面平行，且条块之间传压不传拉；滑坡整体的剩余下滑力，即最后一个条块的剩余下滑力为0。
Table 1 The shear strengthparametersstatistic table of sliding zone of T2b1and T2b3strata
工况
参数
指标
分布区间
样本个数
分布形式
分布函数
天然
（Kpa）
（10.0，43.1）
65
正态分布
（°）
（6.0，23.0）
92
图1 欧家湾庙滑坡工程地质剖面图
Fig.1 The engineering geological profile of Oujiawan landslide
滑坡为老滑坡，滑坡区经过过去的剧烈滑动后，在改变了当时的地形地貌后形成了现今的老滑坡体地形。经对现场的调查踏勘发现，滑坡体上树木歪斜，现仍有滑移变形产生。在滑坡中部多户民房附近，近年每逢雨季都有蠕动滑移。从地表调查和发展趋势上看，目前该滑坡处于不稳定状态。
滑坡计算参数反演分析的优化算法
1 引言
在滑坡稳定性计算和工程设计中，滑带土的粘聚力（C）和内摩擦角（ ）取值正确与否至关重要。目前确定滑带土抗剪强度参数（C、 ）值的方法有试验、工程类比和反演分析3种。滑带土剪切试验分为现场或室内两种，受试样和试验条件的限制，滑带土试验数据通常很离散，需要进行分析计算来确定。工程类比法在确定滑带土的抗剪强度参数时具有很强的主观性，在确定类比指标时又受到类比滑坡客观条件的限制。反演分析是确定滑带土抗剪强度参数的一种有效的方法，根据滑坡的宏观变形状况假设滑坡的稳定性系数，再反算滑带土抗剪强度参数。反算是滑坡稳定性计算的逆过程，得到的参数更符合滑坡的变形情况，参数可以作为试验数据选取的参考，若没有试验数据时，可以直接作为稳定性计算、工程设计的参数。