滑坡计算参数反演分析的优化算法1 引言在滑坡稳定性计算和工程设计中，滑带土的粘聚力（C）和内摩擦角（）取值正确与否至关重要目前确定滑带土抗剪强度参数（C、）值的方法有试验、工程类比和反演分析3种滑带土剪切试验分为现场或室内两种，受试样和试验条件的限制，滑带土试验数据通常很离散，需要进行分析计算来确定工程类比法在确定滑带土的抗剪强度参数时具有很强的主观性，在确定类比指标时又受到类比滑坡客观条件的限制反演分析是确定滑带土抗剪强度参数的一种有效的方法，根据滑坡的宏观变形状况假设滑坡的稳定性系数，再反算滑带土抗剪强度参数反算是滑坡稳定性计算的逆过程，得到的参数更符合滑坡的变形情况，参数可以作为试验数据选取的参考，若没有试验数据时，可以直接作为稳定性计算、工程设计的参数目前，滑带土抗剪强度参数反演分析的方法分为单参数反演和双参数反演两种前者假定一个参数已知的前提下，反算另外一个参数，通常选择对滑坡稳定性影响较敏感的作为未知参数后者在反演中有两个未知的参数，通常选择两个距主滑动面等距的剖面建立极限平衡方程求解本文以三峡库区太山庙滑坡为例，在C、值未知的情况下，综合采用经验类比和反演分析方法确定滑带土的抗剪强度参数，分析时兼顾了滑坡的区域相似性和个体特性，所得到的结果更为准确、可靠。
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滑坡为老滑坡，滑坡区经过过去的剧烈滑动后，在改变了当时的地形地貌后形成了现今的老滑坡体地形经对现场的调查踏勘发现，滑坡体上树木歪斜，现仍有滑移变形产生在滑坡中部多户民房附近，近年每逢雨季都有蠕动滑移从地表调查和发展趋势上看，目前该滑坡处于不稳定状态。
3 滑带土抗剪强度参数统计对三峡库区二期崩塌滑坡治理工程和三期规前勘（调）察中的崩塌滑坡点的勘察试验资料进行分类统计，得到本区滑带土抗剪强度参数值，可以用于验证和优化反演得到的参数经统计得到适合该滑坡的抗剪强度参数分布函数如表1，图2是滑带抗剪强度参数统计直方图表1 T2b1和T2b3滑带土的抗剪强度参数统计表Table 1 The shear strength parameters statistic table of sliding zone of T2b1 and T2b3 strata工况参数指标分布区间样本个数分布形式分布函数天然（Kpa）（10.0，43.1） 65正态分布（°）（6.0，23.0）92正态分布饱和（Kpa）（8.7，31.4）70对数正态分布（°）（5.3，19.8）89正态分布 图2 T2b1和T2b3地层滑带土抗剪强度参数统计直方图Fig.2 The statistic histograms of shear strength parameters of sliding zone of T2b1 and T2b3 strata对天然状态下滑带土抗剪强度参数的试验值由小到大进行排序，其饱和强度参数应与天然状态的一一对应。
对排序后的参数值进行分段，对各段内的天然强度参数和相应的饱和强度参数求算术平均值，得到相应的关系曲线如图3 图3 滑带土天然和饱和抗剪强度参数关系曲线Fig.3 The relation curve of natural shear strengths parameters and saturation shear strengths parameters of sliding zone4 反演计算模型的建立4.1稳定性计算方法常用于计算滑坡稳定系数的方法有Bishop法、Sarma法、剩余推力法本文采用剩余推力法（residual thrust method，RTM）建立反演分析所需的极限平衡方程条块间的作用力，亦即上一条块的剩余下滑力，其方向与上一条块的底面平行，且条块之间传压不传拉；滑坡整体的剩余下滑力，即最后一个条块的剩余下滑力为0计算模型如图图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 滑坡不同发展阶段的稳定系数Table 2 Stability coefficient for different stages of landslide发展阶段变形性质稳定性系数K变形状态局部变形蠕滑1.05～1.00前缘或后缘变形微弱，地表出现未连通微裂缝整体变形微滑1.00～0.90局部坡面变形异常，陡坎出出现小型局部滑崩，裂缝发育，逐渐连通剧滑<0.90坡面出现鼓丘、挤压变形和较大的大裂隙稳定固结固结>1.05根据前面对欧家湾滑坡宏观变形状况的描述，判断滑坡目前处于蠕滑的局部变形阶段，确定稳定性系数为1.00～1.054.3 敏感性分析选取粘聚力和内摩擦角的均值作为基准参数，当一个参数在某一范围内变化时另外一个参数不变。
根据统计结果确定参数变化范围，参数变化区间为敏感性分析结果如图5，内摩擦角是影响滑坡稳定性的敏感因素，根据敏感性分析图的曲线斜率，计算粘聚力和内摩擦角在同方差范围内对滑坡稳定性的影响程度比值＝1：3 图5 滑坡敏感性分析Fig.5 Sensitivity analysis of landslide4.4 滑带土抗剪强度参数反分析及优化取值上面已经确定了该滑坡的稳定性评估指标，即已知滑坡的稳定性系数，可以通过反算得到滑带土的抗剪强度参数C、选择滑坡主滑动剖面作为稳定性的计算剖面，由于存在两个未知参数，首先，确定C值验算范围，并在该范围内选择一系列验算参数值，再通过反算得到相应的一系列值，即将反算得到的一系列抗剪强度参数值与统计的结果进行对比分析，选择最优的参数作为滑坡稳定性计算和设计的强度参数用于比较的计算公式如下：得到的一系列比较函数值Si中最小的则为最优的滑带土抗剪强度参数，表3为经反演得到的一系列参数和相应的比较函数值表3 反演参数和比较函数值Table 3 Back analysis parameters and comparison function valuesK＝1.0C(Kpa)2021.823.625.427.22930.832.634.436.2φ(º)15.8015.5515.3015.0414.8014.5414.2814.0313.7713.51Si1.6141.6061.6081.6091.6011.7982.2062.6143.0213.429K＝1.05C2021.823.625.427.22930.832.634.436.2φ(º)16.6916.4316.1915.9415.6915.4315.1914.9214.6714.42Si0.9300.9040.8960.8890.8901.0871.4851.9032.3002.707根据反算参数的比较值可知，在稳定性系数K=1.0时，滑带土的粘聚力C和内摩擦角φ分别为27.2和14.80；K=1.05时，粘聚力C和内摩擦角φ分别为25.4和15.94。
5 结论（1）根据欧家湾滑坡宏观变形特征，确定了天然状态下滑坡稳定性评估指标，稳定性系数为1.0～1.052）对三峡库区已有勘察资料的滑坡进行分类统计，得到了适合本滑坡滑带土抗剪强度参数的分布函数，并通过分段求平均值的方法得到天然状态和饱和状态下滑带土抗剪强度参数的相关函数3）根据敏感性分析，认为内摩擦角是影响滑坡稳定性的敏感因子，在同方差变化范围内，粘聚力和内摩擦角对稳定性的影响程度比值为1：34）在滑带土抗剪强度参数反演分析的基础上，将得到的一系列参数值与统计资料进行对比分析得到了适合本滑坡的最优参数值，优化过程兼顾了滑坡的区域相似性和个体的特性参考文献：[1] 高德军，徐卫亚，郭其达．长江三峡大石板滑坡计算参数反分析［J］．河海大学学报（自然科学版），2006，34（1）：74～78．[2] 彭亚明，彭军还，张彬等．三峡库区某滑坡抗剪强度参数的反演分析［J］．桂林工学院学报，2003，23（3）：279～283．[3] 富凤丽，佴磊．中里滑坡反分析及强度系数取值研究［J］．长春科技大学学报，2000，30（2）：165～169．[4] 林鲁生，蒋刚，白世伟等．土体抗剪强度参数取值的统计分析法［J］．岩土力学，2002，23（5）：570～574．[5] 龚玉锋，周创兵，梁轶等．参数反演在岩质高边坡变形与稳定分析中的应用［J］．岩石力学，2003，24（5）：809～811．[6] 吴刚，夏艳华，陈静曦等．可靠性理论在边坡反分析中的运用［J］．岩石力学，2005，27（5）：390～393．[7] 徐汉斌，王军．反算法中滑坡稳定系数的取值问题［J］．四川地质学报，1999，19（1）：86～89．[8] 周平根. 滑带土强度参数的估算方法[J]. 水文地质工程地质, 1998, (6): 30－32.[9] 王乐华，杨学堂，李建林．水布垭瓦屋场滑坡参数反分析研究［J］．三峡大学学报（自然科学版），2005，27（5）：390～393．[10] 李端有，甘孝清．滑坡体力学参数反分析研究［J］．长江科学院院报，2005，22（6）：44～48．[11] 徐青，陈士军，陈胜宏．滑坡稳定分析剩余推力法的改进研究［J］．岩土力学，2005，26（3）：465～470．[12] Duncan, J.M., Stark, T.D. Soil strength from back analysis of slope failures. Proc. ASCE Geotechnical Conference. Slope and Embankments. Berkely.1992[13] Maio, C.D. The influence of pore fluid composition on the residual shear strength of some natural clayey soil. Proceedings of 7th Inter. Congress of Landslides. 1996[14] Popesku,M.E.Backanalysisofslopefailures—Apossibilityorachallenge.proceddingsof7thInternationalIAEGcongress.1994[15] Saito, M. Reverse Calculation method to obtain C and &a slip surface, Proceedings Inter symposium. Landslide. New Delhi. 1980[16] 田斌，戴会超，王世梅．滑带土结构强度特征及其强度参数取值研究［J］．岩石力学与工程学报，2004，23（17）：2887～2892．[17] Skempton, A. W. Residual strength of clays in landslides. Folded strata and the laboratory. Geotechnique,. 1985, 35(1).[18] 徐建平，胡厚田，张安松等．边坡岩体物理力学参数的统计特征研究［J］．岩石力学与工程学报，1999，18（4）：382～386．[19] 张社荣，贾世军，郭怀志．岩石边坡稳定的可靠度分析［J］．岩土力学，1999，20（2）：57～66．[20] 罗冲，殷坤龙，陈丽霞等．万州区滑坡滑带土抗剪强度参数概率分布拟合及其优化［J］．岩石力学与工程学报，2005，24（9）：1588～1593．[21] 倪恒，刘佑荣，龙治国．正交设计在滑坡敏感性分析中的应用［J］．岩石力学与工程学报，2002，21（7）：989～992．[22] 汪洋，殷坤龙，安关峰．滑坡敏感因子的灰色关联分析［J］．岩土力学，2004，25（1）：91～93．[23] 杨和雄，王良元．敏感性分析及其模糊方法［J］．南京邮电学院学报，1998，18（1）：99～108．[24] 潘家铮．建筑物的抗滑稳定和滑坡分析[M]．北京：水利出版社，1980．[25] A. Burton, T. J. Arkell , J. C. Bathurst. Field variability of landslide model parameters. Environmental Geology,. 1998, 35(2-3),100~114.。