岩村滑坡稳定性评价(作业)

岩村滑坡稳定性评价一、目的学会滑坡机理分析、稳定性定价和定量计算的基本方法，了解滑带土抗剪强度指标选择的基本途径，掌握滑坡防治工程要点。二、滑坡概况l、自然地理岩村滑坡位于四川盆地某城市市中区，地处长江和佳江的交汇地带，呈半岛状，土地资源十分紧张。在经济建设迅速发展的80年代，市中区斜坡土地得到了大量的利用，交通线路不断改进，高层建筑逐渐增多。但与此同时滑坡灾害事件也日趋严重，岩村滑坡就是灾害之一。该地区属于亚热带气候，温暖潮湿，雨量充沛，多年平均降雨量在1200mm以上，并常有暴雨出现。长江和佳江是市中区两大地表水系，水位年平均变化幅度达20m以上，平均低水位158m，高水位181m，1981年为百年一遇的特大洪水，水位达193m，正在筹建中的三峡工程，按175m蓄水方案修建大坝，该地区最高拱水位将达205m左右。2、地质概况滑坡区基岩地质构造属川东隔档式褶皱中的一复向斜内部，岩层产状平缓，倾角10°以下，倾向在SW200°～270°范围变化。无明显的断裂构造，优势节理产状：75°∠82°；346°∠81°，263°∠85°。基岩地层为侏罗系（J25）泥岩砂岩互层，为内陆河潮沉积，呈紫红色。相对坚硬的砂岩组成了滑坡区的上部平台状地形，泥岩及崩积物则组成斜坡主体。崩积物(Q4col)主要由砂岩块石及泥岩风化粘土组成，厚度分布特点是斜坡上部薄，中前部相对较厚。人工堆石为近期在砂岩体中开挖地下洞室而堆弃于斜坡后部的基岩大块石。滑坡区属河流侵蚀、剥蚀的低山丘陵地貌，斜坡顷部为平台，河谷岸坡的坡度由上至下逐渐变缓，在纵剖面上呈内凹的地形。下伏基岩相对不透水，为弱含水层。据洞室调查，基岩洞室绝大多数为干洞，偶见裂隙有渗水现象。斜坡地带入渗的地表水则汇集于基岩顷面，形成崩积层中的上层滞水。该地区新构造运动不强烈，属受活断裂包围的稳定地块，地震基本烈度为Ⅵ度。3、滑坡特征滑坡主滑方向为NW方向，后缘有一系列NE-SW方向的拉张裂缝，居民建筑物受到严重影响。据调查，人工洞室开挖于1970-1980年之间，地面裂缝最早发现在1981年。1981年四川盆地普降暴雨，江河水位达百年一遇特大水位。目前，滑坡的活动已严重威胁经由滑坡区的主干公路的正常通车。滑坡现处于蠕滑阶段，且在每年的雨季，位移明显增大。下图为岩村滑坡平面图：岩村滑坡平面图三、稳定性计算3．1计算剖面的选取和计算模型的建立计算剖面的选取和计算模型的建立依据勘察过程中钻孔的布置情况,选取上图的A--A′剖面作为计算剖面。该剖面可以利用最多的钻孔资料,同时也大致反映了主滑线的特征。本文用剩余推力法计算滑坡稳定性系数。将主滑体进行条块划分,建立稳定性计算的模型。在条块划分的过程中,条块的多少和条块界线的位置都会影响最终计算结果的准确性和计算的繁杂性。所以在该剖面的条块划分过程中,考虑了地形、地下水和滑动面的形态因素,做到尽量准确又保证计算的简易性。在分析过程中,考虑了以下七种工况:，见下表1所示表1工况1、旱季（175m库水位）2、饱和地下水（暴雨，175m库水位）3、旱季205m库水位4、205m库水位+饱和地下水（暴雨）5、旱季205m库水位+地震烈度VI度6、205m库水位+饱和地下水（暴雨）+地震烈度VI度7、205m库水位下降至175m库水位分析以上七种工况，最基本的要素为：1、旱季正常库水位。2、雨季正常库水位。3、旱季205m库水位。4、雨季205m库水位。即对应表1前四种工况。后两种工况（5、6）只要在稳定性计算时加上水平地震力的影响即可，岩村地震基本烈度为Ⅵ度，故按地震烈度确定的水平地震影响系数为0.064。第7种工况为？？？？？？？？下图2至图5为基本工况下的岩村滑坡A--A′主剖面稳定性计算条分图。岩村滑坡共划分了13块,在坡面线、滑动面的明显拐点的地方以及地下水进入和退出滑坡体的地点都有分界线。（黑色线为地质界线，蓝色线为水面，红色线为条块划分线。）1：1350（cm）图2旱季正常库水位1：1350（cm）图3雨季正常库水位1：1350（cm）图4旱季205m库水位1：1350（cm）图5雨季205m库水位3.2岩体力学参数的选取3.2.1滑坡体各滑块容重的确定组成滑坡体各滑块的岩性不同且不均一，需要根据钻孔资料确定岩性分布后计算出每个滑块的容重，以下表2为各个钻孔滑坡体部分的岩性含量表和各钻孔滑坡体部分的容重平均值，可见，容重平均值相对接近，故可取其容重平均值的平均值代表滑坡体各滑块的容重。经计算，订滑坡体各滑块容重为23.26kN/m3。岩性k1-1k1-2k1-3k1-4k1-5k1-6R(kN/m3)粘土(相对含量)0.0250.0880.6000.4000.400020.188泥岩(相对含量)00000024.794砂岩(相对含量)0.9750.9120.4000.6000.600024.598R平均24.4924.2121.9522.8322.830-------表2各个钻孔滑坡体部分的岩性含量表和容重平均值3.2.2滑坡体各滑块抗剪强度指标的确定由于滑坡岩土体的力学参数的获取过程中存在着诸如系统误差、人为误差以及试样获取的随机性,使其参数存在着离散性和随机性。因此,需要进行统计和分析,以获得相对合理的参数。在统计分析的过程中,主要是对滑带土的抗剪试验结果抗剪强度指标进行统计,得到滑动面岩土体内聚力C和内摩擦系数tanΦ的均值和方差。由于滑体已经有大的变形,处于蠕滑阶段,而且从钻孔资料中也可以看出滑带土中已经有滑动擦痕出现,所以这里选取的参数为残余强度参数。统计结果见表3。残余强度序号123456tanφr(°)0.2221520.128990.1250.4142140.2207790.230868Cr(KPa)12.847.0510.712.077.654.9序号78910111213tanφr(°)0.3057310.2493280.2235260.2276530.2962130.2125570.19438Cr(KPa)12.520.18.8214.018.827.92.8序号1415161718均值方差tanφr(°)0.2493280.1763270.2679490.2034520.236390.2324910.004042Cr(KPa)4.68.17.33.44.18.20388919.7459表三滑带土抗剪强度指标统计表由于残余强度参数不全，而峰值强度参数很全。建议采用如下方法确定缺少的残余强度：内摩擦系数计算值在峰值内摩擦系数的60%~90%之间选取，内聚力计算值在峰值内聚力的10%~30%之间选取。表三数据由此而来。在滑坡体稳定性计算时，可取内聚力C和内摩擦系数tanΦ的均值为各滑块的抗剪强度指标。3.3各工况稳定性计算自编C语言程序，分工况计算岩村滑坡的稳定性，得稳定性结果如下表四：工况K现今天然状态0.990饱和地下水（暴雨）0.962205m库水位0.889205m库水位+饱和地下水（暴雨）0.870205m库水位+地震烈度VI度0.801205m库水位+饱和地下水（暴雨）+地震烈度VI度0.785205m库水位下降至175m库水位0.867附、滑坡稳定性计算C语言程序#includeiostream.h#includestdio.h#includestdlib.h#includemath.h#definesum13//#definelieduxishu0.04#definelieduxishu0#definerongzhong23.26#defineshuirongzhong9.8#defineccc15.98527778#defineyyy0.274822381#definebili1voidmain(){structpoint{floaty,x;};structthreepoint{structpointone,two,three;}unit[sum+1]={0};floatg[20]={0},a[20]={0},f[20]={0},q[20]={0},u[20]={0},v[20]={0},c[20]={0},y[20]={0},l[20]={0},next[20]={0},k={0},kk={0},fr[20]={0},fs[20]={0},frsum={0},fssum={0};inti;chartt;FILE*fp;if((fp=fopen(根据情况改变数据源.txt,r))==NULL){printf(cannotopenthisfile\n);exit(0);}for(i=0;isum+1;i++){fscanf(fp,%f\t%f\t\n,&unit[i].one.y,&unit[i].one.x);fscanf(fp,%f\t%f\t\n,&unit[i].two.y,&unit[i].two.x);fscanf(fp,%f\t%f\t\n,&unit[i].three.y,&unit[i].three.x);}fclose(fp);for(i=0;isum;i++){g[i]=(unit[i].one.y-unit[i].three.y+unit[i+1].one.y-unit[i+1].three.y)*(unit[i+1].one.x-unit[i].one.x)*0.5*rongzhong*bili*bili;a[i]=atan((unit[i].three.y-unit[i+1].three.y)/(unit[i+1].three.x-unit[i].three.x));q[i]=lieduxishu*g[i];l[i]=bili*sqrt((unit[i+1].three.x-unit[i].three.x)*(unit[i+1].three.x-unit[i].three.x)+(unit[i+1].three.y-unit[i].three.y)*(unit[i+1].three.y-unit[i].three.y));v[i]=((unit[i].two.y-unit[i].three.y)*(unit[i].two.y-unit[i].three.y)-(unit[i+1].two.y-unit[i+1].three.y)*(unit[i+1].two.y-unit[i+1].three.y))*shuirongzhong*0.5*bili*bili;u[i]=((unit[i].two.y-unit[i].three.y)+(unit[i+1].two.y-unit[i+1].three.y))*shuirongzhong*l[i]*0.5*bili;}for(i=0;isum;i++){c[i]=ccc;y[i]=yyy;}k=1;do{for(i=0;isum;i++){if(i0)next[i-1]=cos(a[i-1]-a[i])-(y[i]/k)*sin(a[i-1]-a[i]);if(i==0)f[i]=(g[i]*sin(a[i])+(q[i]+v[i])*cos(a[i]))-((c[i]*l[i]/k+(g[i]*cos(a[i])-u[i]-(q[i]+v[i])*sin(a[i]))*y[i]/k));if(i0)f[i]=(g[i]*sin(a[i])+(q[i]+v[i])*cos(a[i]))-((c[i]*l[i]/k+(g[i]*cos(a[i])-u[i]-(q[i]+v[i])*sin(a[i]))*y[i]/k))+f[i-1]*next[i-1];if(i==(sum-1))f[i]=0;if(f[i]0)f[i]=0;}for(i=0;isum;i++){fr[i]=c[i]*l[i]+(g[i]*cos(a[i])-u[i]-(q[i]+v[i])*sin(a[i]))*y[i];if(i==0)fs[i]=(g[i]*sin(a[i])+(q[i]+v[i])*cos(a[i]))-f[i];if(i