湿堆尾矿双向排水大变形固结计算模型.pdf

水利水电科技进展,2014,34(5)摇 Tel:[phone]摇 E鄄mail:jz hhu. edu. cn摇 http:/ / kkb. hhu. edu. cn 第 34 卷第 5 期 Vol. 34 No. 5 水 利 水 电 科 技 进 展 Advances in Science and Technology of Water Resources 2014 年 9 月 Sep. 2014 基金项目:国家自然科学基金(51168004) 作者简介:欧孝夺(1970),男,广西来宾人,教授,博士,主要从事环境岩土工程及尾矿库安全研究。 E鄄mail:ouxiaoduo163. com DOI:10. 3880/ j. issn. 10067647. 2014. 05. 006 湿堆尾矿双向排水大变形固结计算模型 欧孝夺1,2,廖有芳1,2,苏摇 建1,2,潘摇 鑫1,2,李结全1,2 (1. 广西大学土木建筑工程学院,广西 南宁摇 530004; 2. 广西大学工程防灾与结构安全教育部重点实验室,广西 南宁摇 530004) 摘要:根据上游法湿堆尾矿的工艺特点,在分析上游法湿堆尾矿渗流固结特性和存在问题的基础 上,结合一种适合于上游法尾矿库的排水固结装置的排水特点,在考虑竖向和径向双向排水、变渗 透系数及变荷载条件下,采用大变形固结理论推导了尾矿固结计算模型,进而得到以超静孔隙水压 力表示的考虑双向排水作用的大变形固结控制方程,并将该方程与 Barron 固结方程进行比较,验 证了计算模型的正确性。 采用有限差分法对该大变形固结控制方程进行求解,得到了固结控制方 程的差分方程,并利用追赶法等数学方法对差分方程进行了求解。 以中铝广西分公司某尾矿库为 计算实例,采用该排水固结装置的大变形固结计算模型及数值计算程序进行计算,结果表明该计算 模型能有效反映尾矿的固结规律,采用排水固结装置尾矿固结时间缩短了 60 %以上。 关键词:上游法尾矿库;湿堆尾矿;固结装置;大变形固结;径向排水;竖向排水;变荷载;变渗透系数 中图分类号:TV649;TD鄄05摇 摇 摇 文献标志码:A摇 摇 摇 文章编号:10067647(2014)05002807 Large鄄strain consolidation calculation model on both horizontal and vertical drainage of wet filling tailings/ / OU Xiaoduo1, 2, LIAO Youfang1, 2, SU Jian1, 2, PAN Xin1, 2, LI Jiequan1,2(1. College of Civil Engineering and Architecture, Guangxi University, Nanning 530004, China; 2. Ministry of Education Key Laboratory of Disaster Prevention and Structural Safety, Guangxi University,Nanning 530004, China) Abstract: According to the wet landfill process of upstream tailings, its seepage consolidation characteristics and existing problems are analyzed and additionally, it is put forward based on a drainage device that is suitable for about upstream tailings. With the establishment of the mathematical model, it was possible to derive the formula that considers vertical and radial drainage, varying permeability coefficient, and variable load condition at the same time. The governing equations of large鄄strain consolidation with both horizontal and vertical drainage are reduced, which is denoted by excess pore water pressure. By comparing study of that governing equations and Barron consolidation formula, we assessed the accuracy of the proposed model. Using the finite difference method to derive the resulting of large strain consolidation of differential equations, thus obtaining the differential equations of consolidation equations an then solving it by using mathematical methods like Pursuit Method. Using the Guangxi Branch of Chalco aluminum tailings as case study, the results show that this calculation model can effectively reflect the consolidation law of the tailings slurry, and the consolidation time of tailings slurry is shortened by more than 60%. Key words: upstream tailings;wet filling tailings; drainage device; large鄄strain consolidation; horizontal drainage; vertical drainage; time dependent loading; variable permeability coefficient 摇 摇截至 2011 年底,我国已建和在建尾矿库达到 11946 座,其中 90% 以上的尾矿库采用上游法筑 坝1。 此类上游法尾矿库存在坝体稳定性差和尾 矿浸润线高、垂向渗透性差、排水固结缓慢等缺点, 容易出现安全隐患2鄄5。 为了提高尾矿库安全稳定 性,许多学者从絮凝剂、排水工艺、固结措施等方面 开展了加快尾矿排水固结的相关研究,并取得了一 定的研究成果6鄄10。 尾矿的排水固结属于大变形固 结范畴,而目前国内外鲜见涉及尾矿大变形固结的 研究成果。 本文基于大变形固结理论开展上游法湿 堆尾矿固结计算模型研究,以期为上游法湿堆尾矿 固结计算提供借鉴。 1摇 上游法湿堆尾矿渗流固结特性 由于泵送吹填入库的需要,初期排入尾矿库中 的矿泥含水量极高、流动性强、孔隙比大、渗透系数 82 水利水电科技进展,2014,34(5)摇 Tel:[phone]摇 E鄄mail:jz hhu. edu. cn摇 http:/ / kkb. hhu. edu. cn 高、变形能力强,具有大变形固结性质。 根据上游法湿堆尾矿的堆填工艺,尾矿颗粒大 小在水平向(径向)及竖向上均呈不均匀分布,而尾 矿粒径分布直接影响其渗透系数,从而导致渗透系 数在径向及竖向呈现交错分布的情况11。 细粒层 渗透系数可能比粗粒层渗透系数小 1 或 2 个数量 级,因此细粒薄层的存在对竖向上的渗流影响很大, 不可忽略。 由于尾矿泥为分层堆填,对于下层尾矿泥而言, 其上覆荷载呈周期性分级增加。 在新一层尾矿泥排 入库前,尾矿库任意深度处矿泥在自重及上覆荷载 作用下已达到了一定的固结度,但孔隙水压力未完 全消散,土体未完全固结。 当新层尾矿吹填入库后, 即在表层增加了新的荷载,库中矿泥在新荷载作用 下会重新发生固结作用。 在停止排放矿泥的时间里,尾矿泥浆在库内流 动过程中在竖向上粗颗粒首先发生下沉,悬浮于表 面的都是颗粒十分细小的悬浮质颗粒,待表面水分 蒸发完全后便在表面留下厚度极薄的细粒层,即在 矿泥表面会出现一层厚度 2 3 mm 的泥皮,其渗透 系数极小,这就大大地降低了尾矿泥在竖向上的渗 透性,也使得尾矿泥的渗流特性更加复杂。 图 1摇 排水装置示意图 综上所述,上游法湿堆尾矿固结具有变渗透系 数、变荷载、大变形等特性。 2摇 上游法湿堆尾矿排水装置设计 基于上游法湿堆尾矿堆填工艺及排水特点,黄 绍铿等12设计了一种适用于上游法湿堆尾矿的排 水装置,如图 1 所示。 该装置主要由上部悬浮装置 和竖向排水带两部分组成,其中竖向排水带可卷曲 于悬浮装置中。 尾矿泥中的孔隙水经水平(径向) 流向排水带,而后由排水带将水排入底部排渗系统, 使得尾矿库中矿泥在原有竖向排水的基础上,增加 了径向的排水路径,从而加快矿泥排水、缩短固结时 间。 由于初期矿泥含水量极高,上部悬浮结构所用 为轻质材料,故在浮力作用下可实现悬浮结构随尾 矿堆积的同步升高。 在尾矿库平面布置若干个该排水装置(图 2),可 大大缩短尾矿泥的水平排水渗径而加快尾矿泥的排 水进程,在很大程度上降低尾矿库的浸润线,提高尾 矿坝的稳定性。 本文结合该装置的特点,开展大变形 固结计算分析,为该装置的推广应用提供理论参考。 图 2摇 排水装置现场布设剖面示意图 3摇 计算模型的建立 3. 1摇 数学模型简化 a. 尾矿泥的堆填具有成层性,当堆填至第 n 层 时,堆填尾矿泥的总厚度为 H,各层的矿泥厚度依次 记为 h1,h2,hi,hn,如图3 所示(图中 Ki(i=1, 2,n)为各层矿泥侧压力系数;zi(i=0,1,n)为 纵向轴线分层标号),尾矿泥初期坝为透水坝,底部 具有排水特性;堆场顶面处忽略蒸发等微弱影响因 素,为不排水边界条件。 图 3摇 尾矿库成层堆填示意图 b. 根据尾矿泥的堆填排放工艺,作用在某层矿 泥上的上覆荷载是分级加载的变化荷载,故外荷载 按线性变化考虑,其与时间的关系曲线如图 4 所示。 外荷载直至尾矿库满库停排才停止增加,上覆荷载 与时间的关系13可表示为 q(t) = t - tu(i-1) tci- tu(i-1) qui+移 i-1 j = 1 quj摇 t 沂 tu(i-1),tci 移 i j = 1 quj摇 t(沂tu(i-1),tc i (1) 92 水利水电科技进展,2014,34(5)摇 Tel:[phone]摇 E鄄mail:jz hhu. edu. cn摇 http:/ / kkb. hhu. edu. cn 式中:qui为第 i 级最终荷载值;tci,tui分别为第 i 级荷 载增加结束时间和稳定荷载结束时间;t 为总时间。 图 4摇 外加荷载与时间关系曲线 c. 排水装置在尾矿库堆场平面按梅花形布置, 每个竖向排水井的影响范围为六边形(图 5),可用 一个当量直径为 dp的等效面积的圆来代替该六边形 区域14。 其渗流为辐射流,平面上由于水流对称,可 将圆周面视为不排水面进行处理,即无水流穿过圆周 面15。 在尾矿库堆场分离出一个计算单元进行分析, 该单元外边界考虑为不排水边界,分析时同时考虑径 向和竖向的排水固结,半径为 rd,排水井半径为 rD。 图 5摇 大变形固结模型 3. 2摇 基本假定及坐标系选择 a. 土质均质和完全饱和,仅考虑主固结效应, 不考虑次固结和蠕变效应。 b. 土颗粒及孔隙水为不可压缩物质,固结中所 引发的土体压缩变形完全是由土体中的孔隙水外排 引起的。 c. 径向及竖向均可发生排水固结,水的渗流服 从达西渗流定律。 d. 土的压缩变形仅在竖向发生,而在径向上不 发生压缩。 e. 在同一深度处土层径向和竖向的渗透系数 相等;竖向不同深度处渗透系数不同。 f. 固结过程中将忽略井阻和涂抹效应,假定竖 向排水带的排水能力非常强,水平渗流一经流至竖 向排水带,渗流水立即被排出。 根据上述假定,径向上仅发生渗流水的运动,因 此径向可采用欧拉坐标系。 排水井的影响范围考虑 为圆柱体,那么在平面上可用极坐标系(籽,兹),如 图 6 所示。 图 6摇 土体固结坐标系及土体微分单元 在竖向上,由于土颗粒相对水流发生竖向的移 动,上边界及侧边排水边界均发生了改变,因此竖向可 选择流动坐标系,如图6 所示,孜 为流动坐标系坐标。 3. 3摇 固结方程推导 3. 3. 1摇 应力平衡方程 假定土体微分单元体外加荷载沿深度方向传递 为定值,则土体径向应力 滓籽与竖向应力 滓孜有 滓籽= K0滓孜16的关系,其中 K0为土体侧压力系数,那么竖 向和径向应力平衡微分方程可分别表示为 鄣滓孜 鄣孜 - e籽wg + 籽sg 1 + e = 0(2) 鄣滓籽 鄣籽 = K0 鄣滓孜 鄣籽 = 0摇 摇即 摇 摇 鄣滓孜 鄣籽 = 0(3) 式中: 籽s为土颗粒的密度;籽w为水的密度;g 为重力 加速度;e 为孔隙比。 3. 3. 2摇 液相连续性方程 设土颗粒流动速度为 vs,水流速度为 vw(vw1为 竖向孔隙水流速,vw2为径向的孔隙水流速),则微分 单元体渗流情况如图 7 所示。 图 7摇 微分单元体渗流示意图 竖向和径向流量的增量分别为 dQ1= n 鄣 鄣孜(vw1 - vs)籽d籽d兹d孜(4) dQ2= nvw2d籽 + 鄣vw2 鄣籽 籽d 籽 d孜d兹(5) 摇 摇 考虑尾矿泥浆颗粒具有不可压缩性,微分单元 体积的减小就是孔隙体积变化所致,即 03 水利水电科技进展,2014,34(5)摇 Tel:[phone]摇 E鄄mail:jz hhu. edu. cn摇 http:/ / kkb. hhu. edu. cn dV = dVv= d(evs) = vsde = (1 - n)籽d籽d兹d孜de (6) 式中 dVv为土中孔隙体积变化速率。 综合式(4) (5)(6),得到液相连续性方程: (1 - n) 鄣e 鄣t + n vw2 籽 + n 鄣vw2 鄣籽 + 鄣 鄣 孜 n(vw1- vs)= 0 (7) 3. 3. 3摇 达西渗流方程 在竖向上,考虑土体变形为大变形,孔隙水相对 于土颗粒运动,故在采用达西定律时,速度采用相对 速度 vw1-vs,达西渗流方程可表示为 n(vw1- vs) = kiz(8) 式中:k 为渗透系数;iz为 z 方向上的水力梯度。 在 径向上,由于只发生孔隙水渗流,而土颗粒不发生移 动,因此可采用径向孔隙水的绝对速度 vw2,达西渗 流方程表达式如下: nvw2= ki籽(9) 式中 i籽为 籽 方向上的水力梯度。 通过任意时刻 z 和 籽 方向孔隙水压力 ue反解 得到 iz和 i籽,代入式(8)(9)得到达西渗透方程: n(vw1- vs) =- k 籽wg 鄣u 鄣孜 (10) nvw2=- k 籽wg 鄣u 鄣籽 (11) 式中 u 为超静孔隙水压力。 3. 3. 4摇 有效应力方程 根据有效应力原理,在任一时刻饱和土体上作 用的总应力 滓 均由孔隙水压力 ue及作用在土骨架 上的有效应力 滓忆共同承担,总应力表达式为 滓 = 滓忆 + ue= 滓忆 + u0+ u(12) 式中 u0为流动坐标系下某微分单元体静水压力,其 表达式为 鄣u0 鄣孜 = 籽wg(13) 则有效应力竖向和径向的微分表达式分别为 鄣滓 鄣孜 = 鄣滓忆 鄣孜 + 鄣u 鄣孜 + 籽wg(14) K0 鄣滓 鄣籽 = K0 鄣滓忆 鄣籽 + 鄣u 鄣籽 (15) 将式(2)代入式(14)、式(3)代入式(15)可得 d滓忆 de 鄣e 鄣孜 + 鄣u 鄣孜 - 籽sg - 籽wg 1 + e = 0(16) 鄣u 鄣籽 =- K0 d滓忆 de 鄣e 鄣籽 (17) 3. 3. 5摇 以超静孔压为变量的控制方程 将达西渗透方程代入液相连续性方程可得 (1 - n) 鄣e 鄣t - k 籽 1 籽wg 鄣u 鄣籽 - 鄣 鄣籽 k 籽wg 鄣u 鄣 () 籽 - 鄣 鄣孜 k 籽wg 鄣u 鄣 () 孜 = 0(18) 摇 摇 把方程(18)化为以超静孔隙水压力表示的考 虑径向和竖向排水作用的大变形固结控制方程: 鄣 鄣孜 k 籽w () g 鄣u 鄣孜 - (1 + e)鄣滓忆 鄣e k 籽wg 鄣2u 鄣孜2 = 鄣u 鄣t - 鄣q 鄣t (19) 将流动坐标系下的式(19)转换为固相坐标系下的 表达式: (1 + e) k 籽 1 籽wg 鄣u 鄣籽 + (1 + e) k 籽wg 鄣2u 鄣籽2 + 1 1 + e 鄣 鄣e k 籽w () g 鄣u 鄣z + 1 1 + e k 籽wg 鄣2u 鄣z2 = 鄣e 鄣滓忆 鄣q 鄣t - 鄣u 鄣 () t 即摇 A籽( 1 籽 鄣u 鄣籽 + 鄣2u 鄣籽2 ) + Bz 鄣u 鄣z + Cz 鄣2u 鄣z2 = 鄣u 鄣t - 鄣q 鄣t (20) 其中A籽=- k(1 + e) 籽wg琢 琢 =- 鄣e 鄣滓忆 Bz=- 1 1 + e 鄣滓忆 鄣e 鄣 鄣z k 籽w () g Cz=- 1 1 + e 鄣滓忆 鄣e k 籽wg 式中:A籽为径向固结系数:Bz、Cz为竖向大变形系数 项。 式(20)即为双向排水大变形固结计算模型建 立的模型控制方程。 3. 4摇 固结方程比较 Barron 固结理论与太沙基一维固结理论的基本 假定相同,均是建立在小变形基础上的。 如按照这 些假定,具体考虑如下:渗透系数 k 不变,压缩系数 琢=- 鄣e 鄣滓忆,且为常数,引入固结系数 c = k(1+e) 籽wg琢 也为 常数;因为所考虑变形为小变形(变形量微小),因 此可认为流动坐标与固相坐标微量相同,即 d孜=dz; 此外,考虑外荷载是一次施加的,那么鄣q 鄣t =0。 结合 以上分析,式(20)简化后的公式形式如下: c 1 籽 鄣u 鄣籽 + 鄣2u 鄣籽2 + 鄣2u 鄣z () 2 = 鄣u 鄣t (21) 摇 摇 由于竖向与径向渗透系数相同,固结系数也相 同,此时式(21)即为 Barron 固结方程。 即本文推导 的考虑双向排水的大变形固结控制方程经条件简化 后可变为 Barron 固结方程,进一步验证了计算模型 公式(式(20)的正确性。 13 水利水电科技进展,2014,34(5)摇 Tel:[phone]摇 E鄄mail:jz hhu. edu. cn摇 http:/ / kkb. hhu. edu. cn 4摇 差分法求解控制方程 式(20)所表示的大变形固结控制方程为高阶 偏微分方程,一般情况下直接求解方程以获得解析 解不太可能,下面采用有限差分法求解。 4. 1摇 边界条件 根据固结计算模型(图 5),固结方程式需要满 足的一些求解条件如下: a. 初始条件。 初始状态下有效应力全部由孔 隙水承担,数学表达式为 t = 0: ui(z,0) = 籽忆gz = (Gs- 1)籽wgz摇 摇 (zi臆 z 臆 zi-1;i = 1,2,n) (22) 式中:籽忆为浮密度;Gs为土粒相对密度。 b. 上边界条件。 由于不考虑大气蒸发、降水等 因素,因此上边界为不排水边界,数学表达式为 z = 0 颐 鄣u(0,t) 鄣z = 0(23) 摇 摇 c. 下边界条件。 由于尾矿初期坝为透水坝,底 部排水设施排水能力很大,水流流至此处便被完全 排走,因此下边界为排水边界,数学表达式为 z = H:u(H,t) = 0(24) 摇 摇 d. 内径边界条件。 由于竖向排水带的排水能力 很大,因此内径边界考虑为排水边界,数学表达式为 籽 = rD:u(rD,0) = 0(25) 摇 摇e. 外径边界条件。 由于排水固结装置在尾矿 库平面位置上呈梅花形布置,渗流为辐射流,在平面 上水流对称,因此可将圆周面视为不排水边界,数学 表达式为 籽 = rp:鄣u(rp,t) 鄣籽 = 0(26) 摇 摇 f. 层间连续性条件。 相邻层的连续性条件主 要考虑层间超静孔隙水压力及径向和竖向渗流应符 合连续性条件,将其转换为数学表达式如下: z = zi摇 摇 (i = 1,2,n - 1) ui= ui+1 k 籽wg 鄣ui 鄣z = k 籽wg 鄣ui+1 鄣z k 籽wg 鄣ui 鄣籽 = k 籽wg 鄣ui+1 鄣 籽 (27) 4. 2摇 大变形固结控制方程的有限差分方程求解 高阶偏微分方程的有限差分求解需将求解区域 进行空间和时间的离散划分,此外外荷载也需要进 行相应的离散。 4. 2. 1摇 网格设计划分 如图 3 所示,将厚度为 H 的尾矿库分为 n 层尾 矿泥,那么每层尾矿泥的厚度即为 H/ n,在此基础 上,每层再细划分为 x 个薄层。 在径向上划分总网 格数量为 L,竖向上划分总网格数量为 J,网格的划 分如图 8 所示,径向和竖向网格距分别为 驻籽 = (rp- rD) / L(28) 驻z = H/ J(29) 图 8摇 空间网格划分 4. 2. 2摇 时间步长处理 对时间 t 进行离散即得到时间步长 驻t,时间 t 与 驻t 二者的关系可表示为 驻tr= tr- tr-1摇 摇 (r = 1,2,m) t = r驻 t (30) 4. 2. 3摇 外荷载的离散 在第 r 时间段内,相应的荷载增量即离散荷载 段为 摇驻qr= q(tr) - q(tr-1)摇 摇 (r = 0,1,m) (31) 式中 q(tr)为时间点 tr处的上覆土压力荷载。 4. 2. 4摇 差分方程求解 采用交替方向隐式差分法17求解,该方法的特 点是在 tr与 tr+驻t 之间假设有一个过渡的中间时刻 tr+驻t/2,计算分两步进行: 第一步:从 tr到 tr+驻t/2 对径向(籽 方向)采用 隐式差分格式求解,在竖向对 z 采用显式差分格式 求解,得到的差分方程如下: alu r+ 1 () 2 l-1,j + blu r+ 1 () 2 l,j + clu r+ 1 () 2 l+1,j = fl摇 摇(l = 1,2, 摇 摇 摇 ,L;j = 1,2,J;r = 1,2,m)(32) 其中al=- A(r) 籽,l,j 1 驻籽2 bl= 2 驻t + A(r) 籽,l,j rD+ l驻籽 1 驻籽 + 2A(r) 籽,l,j 驻籽2 cl=- A(r) 籽,l,j rD+ l驻籽 1 驻籽 - A(r) 籽,l,j 驻籽2 fl= 2u(r) l,j 驻t + 2驻q(r) 驻t + B(r) z,l,j u(r) l,j+1 - u(r) l,j+1 驻z + C(r) z,l,j u(r) l,j+1 - 2u(r) l,j + u(r) l,j-1 驻z2 23 水利水电科技进展,2014,34(5)摇 Tel:[phone]摇 E鄄mail:jz hhu. edu. cn摇 http:/ / kkb. hhu. edu. cn 式中:A(r) 籽,l,j为径向固结系数的差分格式;B (r) z,l,j、C (r) z,l,j为 竖向大变形系数的差分格式。 第二步:从 tr+驻t/2 到 tr+驻t 对竖向(z 向)采用 隐式差分格式求解,在径向对 籽 采用显式差分格式 求解。 从第一步计算得到的中间时刻 tr+驻t/2 的超 静孔隙水压力变量 u r+ 1 () 2 l,j 可以作为第二步计算的已 知值,得到的差分方程如下: dju(r+1) l,j-1 + hju(r+1) l,j + pju(r+1) l,j+1 = gj摇 摇 (l = 1,2,L; j = 1,2,J;r = 2,3,m)(33) 其中dj=- C r+ 1 () 2 z,l,j 1 驻z2 hj= 2 驻t + B r+ 1 () 2 z,l,j 驻z + 2C r+ 1 () 2 z,l,j 驻z2 pj=- B r+ 1 () 2 z,l,j 驻z - C r+ 1 () 2 z,l,j 驻z2 gj= 2驻q 驻t + 2 驻tu r+ 1 () 2 l,j + A r+ 1 () 2 籽,l,j rD+ l驻籽 u r+ 1 () 2 l+1,j - u r+ 1 () 2 l,j 驻籽 + A r+ 1 () 2 籽,l,j u r+ 1 () 2 l+1,j - 2u r+ 1 () 2 l,j + u r+ 1 () 2 l-1,j 驻籽2 摇 摇 将式(32)和(33)转换为系数矩阵为对角占优 的三对角方程组,采用追赶法可快速求解得到 tr到 tr+驻t/2 中间时刻超静孔隙水压力值 u r+ 1 () 2 l,j ;将其作 为 tr+驻t/2 到 tr+驻t 时段的已知值即可求解时段末 超静孔隙水压力 u(r+1) l,j 的值。 至此便完成了一个时 间步 驻t 的计算,依此求解方法交替进行差分方程的 求解,即可求解各个时间段的待求参数。 计算得到各时刻下的超静孔隙水压力 u 后,通 过有效应力、孔隙比与有效应力之间的非线性关系 式经公式换算得到各个时刻的有效应力 滓忆及孔隙 比 e,而后通过计算得到各薄层对应时刻下的应变 固结度、应力固结度、沉降量等。 平均应变固结度和平均应力固结度分别为 Us= st s肄 = 移 J j = 1 移 L l = 1 cchj 1 + e0llg 滓忆 l 滓忆 0 l 移 J j = 1 移 L l = 1 cchj 1 + e0llg 滓忆 f 滓忆 0 l = 移 J j = 1 移 L l = 1 cchj 1 + e0llg 1 + q - ul 滓忆 0 l 移 J j = 1 移 L l = 1 cchj 1 + e0llg 1 + q 滓忆 0 () l (34) Up= 軍滓忆 軍滓忆f = 移 J j = 1 移 L l = 1 滓忆 0l + q - u () l 移 J j = 1 移 L l = 1 滓忆 0l () + q (35) 式中:cc为土的压缩指数;hj为各层土压缩量;q 为 作用在土层上随时间变化的外加荷载;滓忆 l、滓忆0l分别 为径向第 l 个单元格的有效应力和施加外荷载前的 有效应力;滓忆 f 为固结完成时单元格的有效应力。 5摇 工程实例应用计算分析 以位于广西平果县境内的中铝广西分公司某尾 矿库为例,吹填入库的平果铝尾矿每层矿泥厚度为 0郾 5 m 左右,而后待充填层排水固结由初期含水率 70%降至 40%左右后再重新充填新的一层,而这段 固结过程所需要的时间为 40 d 左右。 采用图 1 所示排水固结装置,利用大变形固结 计算模型及数值计算程序进行求解计算。 计算时土 粒相对密度等于 3郾 26,初始渗透系数取值为 2郾 34伊 10-3cm/ s,初始孔隙比为 1郾 83。 以两层矿泥为计算 模型,每层土厚 0郾 5 m。 通过输入初始条件、边界条 件及实验数据拟合得到的渗透系数和孔隙比与有效 应力之间的函数关系,求解得到每一土层的超静孔 隙水压力分布,进而求解得到每层土层固结度的时 空分布情况。 选取土层平均固结度达到 60% 时所 需时间作为评定标准,对排水半径 rp=1郾 0 m 工况进 行计算分析,结果如下: a. 孔隙水压力在竖向上沿着土层中部并不是 对称分布的,最大孔隙水压力出现在土层中部以下 位置。 随着固结作用的进行,最大孔隙水压力作用 位置逐步向土层中部转移,当靠近固结末期时,孔隙 水压力在竖向上沿土层中部基本上呈对称分布。 b. 获得了底层及上部土层在 rp=1郾 0 m 条件下 的平均应力及应变固结度与时间关系曲线,发现在 相同的条件下,任意时刻平均应变固结度都要大于 平均应力固结度(UsUp),即尾矿泥超静孔隙水压 力的消散速度较沉降速度慢。 c. 在 rp=1郾 0 m 条件下,平均应力固结度达到 60%时,上层矿泥需要 19 d 左右,底层需要的时间 为 16 d 左右;平均应变固结度达到 60% 时,上层矿 泥需要 16 d,底层矿泥需要 14 d。 通过上述工程实例分析可见,尾矿土强度增长 初期呈非对称性;以孔隙水压测算固结度预估土体 强度较沉降测算偏保守;排水固结装置能缩短尾矿 固结时间 60%以上,该大变形固结计算模型能有效 反映尾矿固结规律,在尾矿类超软土固结计算中具 有很大的实用性。 6摇 结摇 语 本文结合一种新型尾矿排水装置,将大变形固 结理论引入该装置作用下尾矿固结计算模型,推导 得到以超静孔隙水压力为变量,可同时考虑径向和 33 水利水电科技进展,2014,34(5)摇 Tel:[phone]摇 E鄄mail:jz hhu. edu. cn摇 http:/ / kkb. hhu. edu. cn 竖向双向排水、变荷载、变渗透系数的大变形固结控 制方程,并利用追赶法得到了有限差分方程的求解 式,进而得到该大变形固结计算模型的数值计算 方法。 工程实例计算结果表明,该排水装置能将尾矿 固结时间缩短 60%以上,计算模型能有效测算尾矿 固结状态及规律,可为该排水固结装置的推广应用 及尾矿质超软土的固结计算提供参考。 参考文献: 1 田文旗,谢旭阳. 我国尾矿库现状及安全对策的建议 J. 中国矿山工程,2009,38(6):42鄄43. (TIAN Wenqi, XIE Xuyang. Tailings pond situation in China and safety countermeasuressuggestions J .ChinaMine Engineering,2009,38 (6):42鄄43. (in Chinese) 2 于斯滢,邵龙潭,刘士乙. 基于有限元极限平衡法的尾 矿坝坝体稳定分析J. 岩土力学,2013,34(4):1185鄄 1190. ( YU Siying, SHAO Longtan, LIU Shiyi. Stability analysis of tailings dam based on finite element limit equilibrium methodJ. Rock and Soil Mechanics,2013, 34(4):1185鄄1190. (in Chinese) 3 刘春刚. 细砂尾矿坝静力和动力稳定性分析研究D. 北京:北京交通大学,2012. 4 李强,张力霆,齐清兰,等. 基于流固耦合理论某尾矿坝 失稳特性及稳定性分析J. 岩土力学,2012,33(增刊 2):243鄄250. (LI Qiang,ZHANG Liting,QI Qinglan,et al. Instability characteristics and stability analysis of a tailings dam based on fluid鄄solid coupling theoryJ. Rock and Soil Mechanics, 2012, 33 ( Sup2 ): 243鄄250. ( in Chinese) 5 张世文,王红艳. 影响尾矿坝安全稳定性因素分析及对 策J. 矿业工程,2004,2(2):61鄄62. (ZHANG Shiwen, WANG Hongyan. Analysis of factors influencing stability of tailing dam and measures for improving dam爷s stability J . MiningEngineering, 2004, 2(2): 61鄄62. ( in Chinese) 6 欧孝夺,曹净,周东,等. 广西平果铝尾矿泥浆化学絮凝 处理试验研究J. 岩土工程学报,2003,25(2):201鄄 203. ( OUXiaoduo, CAOJin, ZHOUdong, etal. Experimental study on chemical flocculating action of waste flush fluid J. Chinese Journal of Geotechnical Engineering,2003,25(2):201鄄203. (in Chinese) 7 尹光志,张千贵,魏作安,等. 尾矿细观结构变形演化非 线性特征试验研究J. 岩石力学与工程学报,2011,30 (8):1604鄄1612. (YIN Guangzhi,ZHANG Qiangui,WEI Zuoan,etal.Experimentalstudyofnonlinear characteristics of deformation evolution for meso鄄scopic structure oftailings J . ChineseJournalofRock Mechanics and Engineering,2011,30(8):1604鄄1612. (in Chinese) 8 李志平,彭振斌,肖尊群,等. 含水率影响下尾矿砂内摩 擦角随标贯击数变化研究J. 岩土力学,2013,34 (5): 1340鄄1344. ( LI Zhiping, PENG Zhenbin, XIAO Zunqun,et al. Variation of internal friction angle of tailings with standard penetration number under influence of water content J. Rock and Soil Mechanics,2013,34 (5): 1340鄄1344. (in Chinese). 9 沈楼燕,王敏杰,王磊. 一种新型尾矿坝同步排渗加固 装置:中国,CN200920188635郾 3P. 2009鄄08鄄19. 10 尹光志,张千贵,魏作安,等. 孔隙水运移特性及对尾矿 细观结构作用机制试验研究J. 岩石力学与工程学 报, 2012, 31 ( 1 ): 71鄄79. ( YIN Guangzhi, ZHANG Qiangui, WEIZuoan, etal.Experimentalstudyof migration characteristics of pore water and its effect on meso鄄structure of tailingsJ. Chinese Journal of Rock Mechanics and Engineering,2012, 31 (1):71鄄79. ( in Chinese) 11 颜学军. 上游法尾矿堆筑坝体沉积规律探讨J. 稀有 金属与硬质合金,2008,36(2):54鄄58. (YAN Xuejun. Deposited law