基于有效应力原理的尾矿坝稳定性评价.docxVIP

基于有效应力原理的尾矿坝稳定性评价 在中国，每年有大约3亿t的采矿尾矿，包括大约1500个尾矿储量。这些储备的最大设计水库为260m，超过100m，10个储水库大于108m，80%的库高小于30m。然而，20%的大型和中型采矿水库占总设计水库的80%以上。同时，上海港水库的最大设计库占总设计库的90%以上。同时，上游法政储备水库占总尾矿渣库的90%以上。随着对矿产资源的巨大需求，水库的数量和高度迅速增加。然而，对尾矿渣容的大规模扩建研究提供的技术支持无法满足大规模建设项目发展的规模和速度。尾矿渣的扩张存在一定的风险，中国没有完整的尾矿渣容技术的综合影响。作者研究了辽宁省风水沟尾矿渣容设计后的现状。风水沟尾矿库位于辽阳县风水沟村，北距辽阳县约9.0km，南距马鞍市约11.0km，西距齐山区约7.5m。尾矿库周围群山环抱，防护良好，东西约4.5m，南北约2.5m，总集水区约7.88km。2.尾矿库由主坝和积坝组成。沉积坝采用上游处理法填充原有矿，多管分散开采，平均外滩比为1.4.5m。原设计的最终坝高为85.0m，可以延长服务时间。尾矿坝可增加至150.0m，总坝高95.0m。根据《尾矿库安全技术规划》aq2005-2006和zbj1-90的相关规定，处理当前坝体，提高坝体当前的稳定性，提高坝体当前的稳定性。同时，我们可以评估水库继续设计是否安全。作者评估了设计现状的条件水库高度以及设计条件下的104m高度。在加工期，对水库进行了静力、动态和地表径流的分析和研究。 1 数值模型的建立 为了进行理论上的计算,需要对实体模型进行必要的简化.该坝最危险计算剖面的选取,由现状条件下的钻孔实测资料决定.首先把钻孔所得岩土颗粒的分布状况加以分析,归并粒径相近的尾砂层,较大的尾砂层兼并小的夹层,最后将属性组成相近的夹层连在一起.每一层底部有断续分布的细小夹层亦归并为大层之内.现状未加高的模型总长约为32 0m,约高87 m.到最终设计高程,推得地质结构分层图及计算模型(如图1所示).模型总长约为320 m,约高104 m.按平面应变问题进行分析,单元统一采用四面体单元,并通过网格检验达到网格最优化.最终该模型被划分为8 187个四面体单元,节点2 698个. 2 流量分析和静力稳定 2.1 水面边界面影响 渗流边界条件:视模型基岩为不透水,初期坝为堆石透水坝;按现场实测数据施加水头边界,堆积坝面为自由透水面.静力边界条件:模型左右两侧加水平向约束;模型底部加竖直方向约束.尾矿砂渗流及静力计算所需参数见表1. 2.2 材料属性的设定 土体的静力本构模型选用摩尔—库仑模型,以反应土体的非弹性及剪塑性行为.坝坡安全系数及潜在滑裂面位置的确定采用有限差分强度折减法.其基本原理是:将坡体强度参数黏聚力c和内摩擦角φ值,同时除以一个折减系数Ftrial,得到一组新的ctrial,φtrial值,其数学表达式可写为 ctrial=1Ftrailc?(1)φtrial=arctan(1Ftrialtanφ).(2)ctrial=1Ftrailc?(1)φtrial=arctan(1Ftrialtanφ).(2) 然后作为新的材料属性参数输入,再进行计算.当满足相应的失稳判据时,对应的Ftrail被称为坝坡的最小安全系数,此时的坝坡达到极限状态,发生剪切破坏,同时可得到临界滑动面. 2.3 坝体内部渗流对坝体稳定性的影响 浸润面发生在孔隙水压力为零与不为零的分界面上,所以用孔隙水压力等值线图既可以显示浸润线的位置,又可以显示坝体内部孔隙压力的分布.从预测的浸润线埋深数据分布图(见图2),可以看到浸润线均没有在堆积坝坝坡处出溢,浸润线距离坝坡仍有一段距离,不会对坝坡产生渗流破坏.从剪应力τxz等值线图(见图3)可以发现坡脚处、坡体内部均呈现压应力状态,坝体沉积滩区域呈现部分拉应力状态,由于其拉应力分布主要集中在坝顶及沉积滩表层区域,故不会对整个坝体的稳定性构成显著影响.由图3还可以发现最大剪应力分布在初期坝及坝体与基岩的接触部位,故这些位置应是工程师重点关注的地方. 坝体潜在滑裂面的位置如图4所示.经计算其抗滑稳定安全系数为1.34,大于现行《选矿厂尾矿设施设计规范》中关于“二级尾矿库在正常运行工况下坝坡抗滑稳定最小安全系数不得小于1.25”的规范要求.因此,该坝是静力稳定的. 3 土体自振及液化危险性分析 使用渗流及静力状态结果作为动力分析的初始状态,计算参数见表2.按最不利工况组合,即洪水与地震同时作用条件下,对该尾矿坝进行加高稳定性分析.由于地震时间很短,故假定在地震期间的孔隙水不向外排出,而是封闭在土体骨架中.动力计算的边界条件:下部边界采用静态边界(黏滞边界)条件;水平向通过施加自由边界条件的办法处理弹性波的能量逸散. 在分析坝体的性状中,不考虑地震引