APP下载

大中河水库顺层岩质边坡覆盖层开挖整治模拟分析

2015-11-24沈新建

水科学与工程技术 2015年5期
关键词:风化层库岸河水库

沈新建

(惠州市水电建筑工程有限公司,惠州516000)

1 前言

随着国内大量水库的修建,形成了不少库岸边坡,特别是顺层岩质边坡,其稳定性问题对水库的正常运行影响很大[1-3]。大中河水库位于云南省思茅市六顺乡境内,距思茅城50km,是一座以农业灌溉为主,兼顾发电、人畜饮水和防洪的中型水利工程。为进一步合理利用水资源,在保证下游河道生态、本流域农业灌溉和生活用水需求后,引水至信房水库,向思茅坝区调水,供给城镇生活和工业用水,解决普洱市城市规模快速发展对水资源的需求问题,支撑经济社会可持续发展。调水工程完工后,预计每年可从大中河水库调入思茅城区水量4730万m3,可满足预计的用水需求。大中河水库思茅坝区引水工程2013年5月8日开工建设,工程概算总投资5.6亿元,工程直接费用4.2亿元。

大中河水库引水工程涉及大量边坡开挖,库岸边坡基岩大部分裸露,且表层覆盖较厚的第四系残坡积层,加之岩层走向近乎与坡面线平行,属于顺层岩质边坡,导致开挖边坡的稳定性很差。因此,为保障引水工程建设的顺利开展,应对库岸边坡进行整治加固。

基于大中河水库地质勘测报告,建立水库典型边坡断面,采用有限元法[4-6]对边坡覆盖层的开挖和加固过程进行模拟分析,旨在对施工过程中边坡的各项指标变化进行监测,并优化加固参数。

2 工程概况

2.1 大中河水库引水工程

大中河水库多年平均来水量7180万m3,水库死水位1039m,正常蓄水位1075.53m;水库总库容5600万m3,兴利库容4327万m3,死库容210万m3。按2030水平年,P=95%供水保证率下大中河水库可向思茅坝区调水3150万m3,以满足思茅坝区需水要求。

工程总体布置,从大中河水库坝后零级电站压力钢管取水,一级泵站提水后绕大中河水库右岸向思茅坝区方向输水,一级泵站提水至1211.5m高程,随后自流经过大中河1#箱涵、大中河水库1#隧洞、大中河水库2#箱涵、大中河水库2#隧洞等14个建筑物,到达位于曼弯村西北侧900m山脊的二级泵站,经二级泵站提水至山顶水池,水位高程1391m,随后自流经过平掌梁子箱涵、丫口寨渡槽、大地山隧洞、曼歇倒虹吸等26个建筑物,于信房水库西侧副坝右岸进入信房水库,线路全长43.64km。

2.2 二级泵站工程地质条件

二级泵站位于曼湾村西北方向900m处,泵房处于山坡坡脚下,所处区域为构造侵蚀浅切割低中山沟谷斜坡地貌,下伏基岩为中侏罗统花开左组(J2h)砂泥岩,岩层微倾向山内,倾角30°~35°,地形坡度一般30°左右,岩层走向与坡面线近平行,属于典型的顺层岩质边坡。另外,边坡表面覆盖残坡积层(Qedl),分布于区内山坡表层,主要为褐红色砂质黏土、粉土,厚度1.5~3m,胶结差,浸水后易于崩解。厂房开挖边坡为顺向坡,表面覆盖层和强风化岩石抗剪强度低,开挖边坡稳定性较差,需及时采取加固处理措施。

3 边坡模型建立与计算参数选取

根据大中河水库地质勘测报告得到了引水工程二级泵站某典型边坡断面内部岩土体岩层走向及分布,边坡断面地层从上至下可划分为:表面覆盖层、强风化带、中风化带和弱风化带,另外,选取断面内含泥质砂岩藏于中风化带内部。

基于相关勘测结果及地质资料,建立了边坡有限元模型,边坡模型高40m,长80m,坡角30°,如图1。

图1 某开挖边坡断面有限元模型及岩层分布

选取不同风化带土样进行相关土工试验,得到了数值计算所需的力学参数,如表1。

表1 边坡岩体有限元计算力学参数

4 计算结果及分析

4.1 边坡位移场

大中河水库库岸边坡覆盖层开挖后边坡发生的水平位移场,如图2。其中图2(a)为开挖不支护,图2(b)为开挖并支护,规定位移向左为正,向右为负。

图2 边坡坡积层开挖后水平位移

由图2(a)可知,当覆盖层开挖后不采取支护时,位于强风化带的岩体出露,在开挖卸荷作用的影响下,岩体应力得到释放,以前受压的岩体变为受拉,所以边坡表面发生面向临空方向的水平位移,且主要集中在边坡中上部位,最大水平位移4.75mm,而边坡坡脚位置发生位移较小;当边坡覆盖层开挖,并及时采取加固措施后(边坡表面采用网格梁植草护坡,边坡中下部加抗滑桩),边坡表面发生位移和范围均大大减小,最大水平位移也仅有1.7mm。说明边坡覆盖层开挖后,及时采取加固措施,可大大减小边坡表面发生的位移。

为进一步分析大中河水库库岸边坡覆盖层开挖前后边坡稳定性变化情况,采用边坡节点最大位移突变作为指标,计算了3种情况不同折减系数[7-9]下边坡的最大节点位移,得到了3种情况下边坡的安全系数,计算结果如图3。

图3 不同工况边坡最大节点位移随折减系数增加变化曲线

分析发现,边坡最大位移随折减系数的增加而逐渐增大,刚开始增加幅度不大,最后发生突变,急剧增加。计算结果表明:边坡在自然状态下,当折减系数增加到1.4时,位移由11mm突然增加到26mm,可认为自然状态下计算边坡的安全系数为1.4;同理,可以得出边坡开挖后(不加固)安全系数为1.9,说明削坡减压有助于提高边坡稳定性,而当采取加固措施后,安全系数进一步增加到2.6。

4.2 边坡破坏模式

通过有限元法自动搜索得到了3种工况下大中河水库库岸边坡的危险滑裂面,如图4。

图4(a)表明,自然状态下边坡的最危险位置位于距离坡脚一定高度的边坡覆盖层,容易在外界作用下发生覆盖层滑动,为了保障工程安全,应当及时对表层覆盖体进行开挖清除;当覆盖层开挖后,边坡稳定性得到一定提升,滑裂面向内扩展,危险滑裂面位于强风化层和中风化层之间,最终沿坡脚位置剪切滑出,如图4(b);当进行加固后,如图4(c)。网格梁可有效控制边坡表面发生的位移,抗滑桩则插入到不稳定地层以下,将强风化层和中风化层连接起来,通过抗滑桩将受力传递到中风化层的岩体,大大增加了边坡的抗滑力,与上述安全系数分析结果吻合。此时,边坡破坏时的滑裂面则进一步向内扩展,位于中风化层和弱风化层之间,说明采用网格梁和抗滑桩加固边坡效果明显。

图4 不同工况下边坡的危险滑裂面

分析可见,对于顺层岩质边坡,边坡的稳定性往往受岩层走向控制。特别是岩层倾角和边坡坡角相当时,边坡处于最不稳定状态,更应引起重视。

4.3 加固参数对边坡稳定性影响

大中河水库边坡加固主要采用了表面网格梁植草护坡和中下部抗滑桩加固处理。其中网格梁护坡作用通过在模型中输入对应强度参数即可,网格梁截面如图5。

图5 网格梁植草护坡设计

抗滑桩采用板单元模拟,选用弹性材料,主要参数包括轴向刚度EA和抗弯刚度EI,表示抗滑桩单位宽度受力(kN/m)。抗滑桩轴向刚度主要体现其抗压和抗拉能力,而抗弯刚度则反映了抗剪断能力,由于顺层边坡破坏模式多属于失稳剪切破坏,故只分析抗滑桩不同抗弯刚度对边坡的稳定性影响。

计算结果如图6。当抗弯刚度在1600kN/m以下时,随着抗弯刚度的增加,边坡的安全系数增加明显;而当抗弯刚度超过1600kN/m时,随着抗弯刚度的增加,边坡的安全系数增加很小,几乎维持在2.8左右。计算结果表明,本工程中抗滑桩抗弯刚度在1600kN/m为宜,太小难以保证加固效果,过大则会造成浪费。

图6 边坡安全系数随抗滑桩抗弯刚度增加变化

5 结语

基于大中河水库引水工程地质勘测报告相关数据,建立了二级泵站厂房边坡断面,对边坡表面覆盖层的开挖和加固过程进行模拟分析,计算结果表明:

(1)泵站厂房边坡覆盖层开挖后,边坡的稳定性得到一定提升,安全系数由1.4增加到1.9,但开挖将导致强风化带岩层出露,局部位置将发生较大位移,采取加固措施后,位移得到有效控制,边坡稳定性进一步提升。

(2)计算得到了不同工况下边坡的危险滑裂面。自然状态下,边坡最危险位置位于坡脚上部的覆盖层,易出现覆盖层滑移;当边坡覆盖层开挖后,滑裂面位于强风化层和中风化层之间;而采用网格梁和抗滑桩加固后,危险滑裂面向内扩展,位于中风化层和弱风化层之间。

(3)计算了抗滑桩不同抗弯刚度对提高边坡稳定性大小的影响,发现抗弯刚度在1600kN/m时,提升边坡稳定性较大,且最为经济。

[1]李杨鹏,薛新华,杨兴国,等.锦屏一级水电站左岸高边坡施工动态仿真模拟[J].水利水电技术,2014,45(8):121-124.

[2]邓华峰,李建林,易庆林,等.软岩高边坡开挖卸荷变形研究[J].岩土力学,2009,30(6):1731-1734.

[3]皱丽春,王国进,汤献良,等.复杂高边坡整治理论与工程实践[M].北京:中国水利水电出版社,2006.

[4]郑颖人,陈祖煜,王恭先,等.边坡与滑坡工程治理[M].北京:人民交通出版社,2007.

[5]黄梦宏,丁桦.边坡稳定性分析极限平衡法的简化条件[J].岩石力学与工程学报,2006,25(12):2529-2536.

[6]王汉辉.边坡稳定的有限元塑性极限分析法研究[D].武汉:武汉大学,2003.

[7]郑颖人,赵尚毅,张鲁渝.用有限元强度折减法进行边坡稳定分析[J].中国工程科学,2002,4(10):57-62.

[8]赵 川,付成华.基于局部强度折减法的某水电站边坡稳定性分析[J].人民珠江,2015,36(3):66-69.

[9]刘小丽,周德培.有软弱夹层岩体边坡的稳定性评价[J].西南交通大学学报,2002,37(4):382-386.

猜你喜欢

风化层库岸河水库
岩质边坡风化层的关键滑动面及其稳定性研究
近地表速度模型精度影响分析
开阳县三口河水库工程高边坡预裂爆破控制
新疆BEJ山口水库近坝库岸HP1滑坡体稳定性分析
开阳县三口河水库(岩溶地区)溶洞处理
合肥地区古近系红层岩土工程特性分析
万州江南新区密溪沟库岸特征及塌岸风险预测
酒泉市洪水河水库坝址方案比选
恰甫其海水库库岸侵蚀坍塌及其防护措施
湾兔河水库除险加固工程洪水调节计算