浅谈数字模拟分析法在滑坡稳定性分析中应用.docVIP

浅谈数字模拟分析法在滑坡稳定性分析中应用 　　摘 要：滑坡稳定性定量计算评价中常采用极限平衡理论方法计算滑坡稳定系数，此方法无法分析稳定破坏的发生和发展过程，无法考虑局部变形对滑坡稳定的影响。因此，本文章采用数字模拟分析法（FLAC法）对该滑坡体进行分析，其计算结算安全可靠，分析结果可为工程设计及施工提供参考。 　　关键词：杨家坝滑坡；FLAC法；滑坡稳定性分析 　　中图分类号：TV22 文献标识码：A 　　1 概述 　　1.1滑坡概况 　　杨家坝滑坡为三峡库区三期地质灾害重点防治项目，为一岩质古滑坡，位于重庆市云阳县普安乡中阳村4-8社境内。地形地貌属于构造剥蚀低山地貌，堆积体岩土种类有粉质粘土夹碎块石（含量变化较大）、泥岩和砂岩组成，滑体整体厚度变化大。杨家坝滑坡属深层大型岩土质混合滑坡，滑坡体长约626m，宽约875m，总体积约855×104 m3，主滑方向约320°。该滑坡前缘分布于磨刀溪的右岸坡脚上，磨刀溪属长江一级支流。 　　该滑坡由于前期构造运动强烈，在全新世前可能发生过多次变形，形成滑体堆积物成分类型较多，滑坡堆积层根据其结构组成分成上下两层，表层约12～25m厚粉质粘土夹碎块石（块石含量2～30%），底部为7～15m厚的滑体岩质层（泥岩与粉砂岩碎裂体）。各滑坡单元后侧至滑坡后缘为近期堆积形成的崩坡积粉质粘土夹碎块石。 　　1.2滑坡基本特征及变形特征 　　杨家坝滑坡区出露地层主要为第四第崩坡积层和侏罗系上统遂宁组（J3s）地层。该段遂宁组地层为一套紫红色、黄灰色细～中粒长石砂岩，局部含钙质紫红色（粉砂质）泥岩。据勘察钻孔资料显示砂岩厚度30m，层位分布稳定，泥岩多为夹层，其厚度小（夹层厚度约0.1～0.5m）多呈薄层状。本滑坡岩层产状315～350°∠5～65°，倾角上陡（60°）下缓（15°），至溪沟底部岩层最缓为5～10°。滑坡岩层倾向与斜坡地形坡向基本一致，属顺向坡（见图1）。 　　滑坡整体地形南东最高，往北西渐降低，滑坡后缘位于地形最高处“灯盏窝”斜坡上，滑壁光滑，其坡度较陡与基岩倾角相近，约40～75°，沿后缘呈不规则的弧形分布，分布有滑壁陡坎，高1～6m不等，坡度33～50°。滑坡两侧为冲沟，坡度15～25°，两侧边缘基本受冲沟切割控制，沿滑坡两侧边缘堆积体呈台阶状分布，坡体中间则可见鼓起状堆积体。 　　据对滑体现有地质情况和调查分析，该滑坡前缘在全新世前出现过强烈变形，在滑动变形过程沿地形向东西两侧较低处分解破坏，现从地貌形态，主滑方向和结构组成上分析，该滑坡体多处特征具有差异性和特殊性，因此，可分成3个滑坡分单元（L1、L2和L3），其中L2滑坡单元土层深度10m左右多见光滑剪切破坏面，说明产生过次级滑动，体积规模小。 　　2 数字模拟分析 　　2.1 模型建立 　　本次对杨家坝滑坡稳定性的数值模拟分析主要采用“快速（显示）拉格朗日分析法”。在本次模型计算过程中，FLAC采用了离散元的动态松驰法。该法一般不需要求解刚度矩阵，其次，相比其它差分分析法，该模型不但能计算滑坡变形量，还可以模拟滑坡沿滑动面产生的滑动变形。同时，FLAC还能针对不同物质特性，采用本构方程来反映滑坡组成物质的动力行为。 　　选取了L1、L2、L3三个滑坡体在正常运行条件下及三峡水库蓄水后的变形破坏状况，并与现场调查成果和极限平衡法计算成果进行比较。 　　本次计算模型采用依据：地形地貌、地质结构、变形特征及滑坡体物理力学特性。通过模型计算和分析，将本次模型的地层L1滑坡单体分为4层，L2和L3滑坡单体分3层。计算模型所用岩土体变形参数以土工试验统计成果为主要依据。滑坡计算模拟模型见图2、图3和图4，为了直观、清楚读懂模型图，、用图中不同颜色代表不同的岩层或岩性：绿色代表粉质粘土夹碎块石土（未区分上下两层），红色代表粉砂岩（强度较低且碎裂状），夹紫红色泥岩，蓝色代表滑坡下覆稳定基岩层，主要为砂岩。 　　2.2 L1滑坡的变形特征及稳定性评价 　　通过数字模拟对滑坡体重点部位的变形状况进行模拟跟踪，判断滑坡现状的稳定性，同时数字模拟搜索潜在破坏面的位置，预测滑坡潜在滑坡变形、破坏因素和稳定性分析。 　　在滑坡体上确定3点变形特征变化较明显部位，根据相关实验数据模拟其变形破坏情况。初始计算时滑坡体选取的点变形量较大，最大位移量达8.2cm，但变形量随着时间推移逐渐趋于稳定，则表示滑坡已处于稳定状态。根据数据模拟结果并对非选定点位的跟踪计算，其计算过程和形成的结论相近。因此，L1滑坡在天然状态下为稳定，模拟结果与滑坡的定性、定量分析基本一致。 　　通过剪应变增量和剪应变速率来搜索最危险滑动面。剪应变增量最大的部位主要沿硬的滑床（基岩）与其上相对较软的滑体两者的接触面分布，这也是对勘查结果的