尾矿坝稳定性分析.docxVIP

尾矿坝稳定性分析 0 尾矿库容积设计概述 在中国，每年从黑色、商品、黄金、化工、原子能工业和建筑材料行业产生约3亿t的尾矿，并在约1500个尾矿库中产生约1500个。这些库中最大设计坝高260m，超过100m的有26座，有10座库容大于108m3，坝高小于30m的小库占80%左右。但是20%的大、中型库的库容占总设计库容的80%。同时上游法尾矿坝占全部尾矿堆积坝的90%以上，其数量及坝体高度随国家对矿产资源的巨大需求而快速增加，而相关尾矿库扩容研究所提供的技术支撑却无法与大型工程建设发展的规模和速度相适应，尾矿库扩容存在一定风险，国内尚无尾矿扩容的全面技术规范。扩容加高后的渗流稳定、动力稳定及各种因素综合影响是扩容设计的关键所在。 辽宁某尾矿坝堆积高度将达到设计高度的50%，据《尾矿库安全技术规程》AQ 2005—2006及《选矿厂尾矿设施设计规范》ZBJ1—90的有关规定，当尾矿坝堆积高度达1/2—2/3最终设计坝高时，应对坝体进行一次全面勘查与评价，以校核坝体现状的稳定性；尾矿库的扩容加高存在与原有设施衔接问题，排洪系统、回水系统、排渗系统、尾矿输送系统、增加副坝工程等改造问题。针对该工程现状条件，坝高118m及加高到170m的工况进行评价，对加高条件下的坝体进行了流固耦合、动力液化研究。计算地质剖面、计算参数由勘查单位提出，加高部分地质结构依据现有的钻孔资料和坝坡比推断给出，洪水位标高以库区水边线距离干滩长度100m确定。 1 计算模型及网格划分 为了进行理论上的计算，需要对实体模型进行必要的简化，该坝最危险计算剖面的选取，来源于现状条件下的钻孔实测资料。坝体需要经过合理的简化和延伸，其基本过程如下：分析钻孔所得到的颗粒分布状况，归并粒径相近的尾矿层，较大的尾矿层兼并较小的夹层，最后将属性组成相近的夹层连在一起。底部有断续分布的细小夹层亦归并为大层之内，底部基岩按钻孔位置处连线。现状未加高的模型总长约545m，高约118m。到最终设计高程，推得地质结构分层图及计算模型如图1所示，模型总长约760m，高约170m。按平面应变问题进行分析，单元都统一采用四面体单元，并通过网格检验达到网格最优化。最终，该模型被划分为26689个四面体单元，节点8940个。 2 流量分析和静力稳定 2.1 模型表面约束 渗流边界条件：模型基岩视为不透水；初期坝视为堆石透水坝；按现场实测水压施加水头边界，堆积坝面为自由透水面。静力边界条件：模型左右两侧施加水平向约束；模型底部施加竖直方向约束。尾矿砂渗流及静力计算所需参数如表1所示。 2.2 材料物理性能 土体的静力本构模型选用摩尔-库仑（Mohr-Coulomb）模型，以反映土体的非线弹性及剪塑性力学行为。坝坡安全系数及潜在滑裂面位移的确定采用有限差分强度折减法。其基本原理是：将坡体强度参数黏聚力c和内摩擦角φ值，同时除以折减系数Ftrial，得到一组新的ctrial，φtrial值，表达式为 然后作为新的材料属性参数输入进行计算。当满足相应的失稳判据时，对应的Ftrial称为坝坡的最小安全系数F，此时的坝坡达到极限平衡状态，发生剪切破坏，同时可得到临界滑动面。 2.3 积坝坝体结构及坝体内部渗流分析 浸润面发生在孔隙水压力为0与不为0的分界面上，所以用孔隙水压力等值线图既可以显示浸润线的位置，又可以显示坝体内部孔隙压力的分布。从预测的浸润线埋深数据分布图（图2）可以看出，浸润线均没有在堆积坝坝坡处出溢，浸润线距离坝坡表面仍有一段安全距离，浸润线埋深均较深，不会对坝坡产生渗流破坏。从剪应力τxz等值线图（图3）可以发现坡脚处、坡体内部均呈现压应力状态，坝体沉积滩区域部分呈现拉应力状态，由于其拉应力分布主要集中在坝顶及沉积滩表层区域，故不会对整个坝体的稳定性构成显著影响。由图3还可以发现最大剪应力分布在初期坝及坝体与基岩的接触部位。 经过计算，该坝潜在的滑裂面的位置如图4所示，其抗滑稳定安全系数为1.49，符合现行的《选矿厂尾矿设施设计规范》中关于“二级尾矿库在正常运行工况下坝坡抗滑稳定最小安全系数不得小于1.25”的规范要求。因此，该坝是静力稳定的。 3 坝体结构安全响应 使用渗流及静力状态结果作为动力分析的初始条件，计算参数见表2。按最不利工况组合，即洪水加地震作用条件下，对该尾矿坝进行加高稳定性分析，由于地震时间很短，故假定在地震期间的孔隙水不向外排出，而是封闭在土体骨架中。动力计算的边界条件：下部边界采用静态边界（黏滞边界）条件；水平向施加自由场边界条件。通过施加自由场边界条件的办法处理弹性波的能量逸散。 在分析坝体的性状中，不考虑地震引起的垂直运动影响，因为理论分析和实验观察均发现，动剪应力和动力破坏影响主要取决于地震水平运动，故本文只按剪切波输入计算。按要求，假定基本地震加速度峰值为