摘要
排放后的尾矿沉积规律影响尾矿库坝体的稳定性。为研究不同排浆浓度条件下尾矿沉积特性,以云南某尾矿库为研究对象,通过现场试验和室内缩尺的堆坝模型试验,研究了不同排浆浓度下的尾矿沉积规律,分析了单颗尾砂颗粒的临界状态,推导了尾砂颗粒运动—沉积临界速度关系和不同粒径颗粒在干滩上的运移距离公式。研究结果表明:沉积滩面的坝前位置粗尾砂沉积较多,随着距离子坝越远,细粒尾砂沉积越多;在干滩面上尾砂流速与矿浆浓度呈负相关,与尾砂浆体压力呈正相关;较粘稠高浓度尾砂浆体的浓度越大,流动力就越弱,致使沉积形成的坡度越大;排浆浓度的逐渐增大,矿浆沉积长度也随之变长;尾砂颗粒在干滩面的沉积距离与排浆初始速度成正比,与干滩坡度、排浆浓度及颗粒粒径成反比。研究成果可对尾矿库坝体稳定性分析提供一定参考。
尾矿库一旦溃坝,将会对人民的生命和财产安全造成巨大威胁,而不同浓度筑坝浆体特性直接影响着尾砂颗粒沉积规律、沉积层倾角以及坝体的沉积结构特征。因此,探索尾矿库不同浓度浆体的沉积规律对尾矿库稳定性影响显得尤为重要。针对尾矿库不同浓度筑坝浆体沉积规律研究,国内外学者通过不同的试验方法开展了大量研究。尹光志
笔者以单颗粒尾砂受力分析为基础,通过研究颗粒受力极限平衡方程,提出尾矿浆体入库过程中固体颗粒状态判据,推导尾砂颗粒在干滩面上的沉积距离公式,并通过模型试验对所提出判据进行验证,探究不同排浆浓度对尾砂沉积规律的影响,对更好地了解尾矿库内部结构及渗流特性具有现实意义。
如

图1 尾矿库放矿堆积图
Fig. 1 Tailings pond drawing accumulation
假设尾砂颗粒为球状,根据文献[
。 | (1) |
式中:为摩擦力;为拖拽力;为上举力;为等效重力;、、、分别为摩擦力臂、拖拽力臂、上举力臂和重力力臂。

图2 尾矿颗粒受力分析
Fig. 2 Force analysis of tailings particles
令
, | (2) |
有
, | (3) |
, | (4) |
, | (5) |
。 | (6) |
式中:为初始流速;、、为摩擦力、拖曳力、上举力系数;为颗粒面积;r为颗粒半径;为颗粒等效重力;为重力加速度;为流体粘度;,分别为颗粒密度和浆体密度。
。 | (7) |
由上式可以得出,当尾砂浆的粘度越大时,尾砂的分离流速越小。尾砂浆的粘度系数大小由物质种类、浓度、温度等因素决定。尾砂浆浓度越大时浆体的粘度越大,流速越小。设尾砂浆体的流速为,则当时,尾砂颗粒处于运动状态;当时,尾砂颗粒处于静止状态。
从1.2对尾砂颗粒的受力分析可知,尾砂颗粒的分离流速与浆体粘度、颗粒密度、颗粒粒径、干滩坡度等有关。尾砂颗粒在干滩面的流动过程中,由于动能的消耗损失,颗粒的流动速度不断减小,当尾砂颗粒的流动速度小于分离流速时,该粒径的颗粒与水流分离,发生沉积,在干滩面上形成坝前粗库尾细的分布规律。引入矿浆流动的水流恒定总流能量方
1)尾矿浆体在干滩面为恒定流,总流流量在流动过程中不发生变化;
2)在矿浆上的质量作用力只有重力;
3)干滩面任意过流断面为渐变流;
4)干滩过流面之间的流量保持不变。
即:
, | (8) |
。 | (9) |
式中:为干滩坡度;, 为动能修正因数;为尾砂颗粒到最远沉积距离时的流速;为初始流速;为定值参数;为干滩两端总水头差;为干滩面之间的总流水头损失。计算可得
。 | (10) |
联合
, | (11) |
式中,为该粒径颗粒的最远沉积距离,可简化为
。 | (12) |
式中:为固体颗粒容重;为矿浆容重;、、、均为定值参数。
不同浓度尾砂筑坝模型试验是指在实验室条件下,按照实体原型,用不同比例尺(缩小、放大尺寸)模型,对工程中可能出现的问题或现象进行研究。通过堆坝模型试验,可以揭示和分析尾矿坝堆筑过程中颗粒沉积规律和坝体结构特性,以验证理论并解决工程实际问题。因此,采用模型试验模拟本文尾矿库的堆筑过程是比较理想的方法,试验的主要目的是研究不同浓度尾砂在筑坝过程中水平向和竖直向颗粒的沉积规律。
试验材料为云南某尾矿库入库铁尾矿。现场勘查结果显示:筑坝尾砂分为尾中砂、尾细砂、尾粉砂和尾粉质粘土4种。为获得堆积尾矿物理力学性质,试验前对尾砂取样进行分析,并进行了土工试验及颗粒级配筛分试验,具体各项土工试验结果见
统计项目 | 容重γs/(g∙c | 饱和度Sr /% | 塑性指数I/% | 含水率W/% |
---|---|---|---|---|
尾中砂 | 1.99 | 94.10 | 12.20 | |
尾细砂 | 2.02 | 81.90 | 14.10 | |
尾粉砂 | 2.04 | 95.50 | 9.63 | 18.90 |
尾粉质粘土 | 1.89 | 98.61 | 14.89 | 23.70 |
为了满足原型与模型的主要参数相似,确保试验结果的相对可靠性,此模型材料均选用现场尾矿库内尾砂,现场尾砂粒径级配如

图3 尾砂粒径级配累积曲线
Fig. 3 Accumulation curve of tailings particle size gradation
不同浓度尾矿沉积模型试验装置主要由以下几部分组成:试验沟槽、流态监测系统、排渗管、水位监测系统、搅拌系统、压力泵排放系统、沉积试验沟槽等,可以模拟现场尾矿排放的多种工况。为了实现尾矿排浆过程与现场相近,搅拌系统和压力泵排放系统包括3个2

图4 模型试验装置
Fig. 4 Model test unit
本次试验中设计了7个不同浓度试验组,分别为20%、25%、30%、35%、40%、45%和50%尾矿浆体(指尾砂浆体中尾砂质量占比),分别从低浓度依次试验,每个浓度试验组沉积时间为24 h(确保每组试验能够沉积足够尾矿满足取样),共7组试验,如
试验组数 | 尾矿砂浆体中尾矿质量百分比/% | 沉积时间/h | 浓度划分 |
---|---|---|---|
1 | 20 | 24 | 低浓度 |
2 | 25 | 24 | 低浓度 |
3 | 30 | 24 | 低浓度 |
4 | 35 | 24 | 低浓度 |
5 | 40 | 24 | 高浓度 |
6 | 45 | 24 | 高浓度 |
7 | 50 | 24 | 高浓度 |
模型试验观测到的浆体流动特性与实际尾矿库基本一致,采用模型堆坝试验对不同浓度尾砂筑坝沉积规律进行对比模拟研究,展现了不同浓度尾砂在整个堆积过程中的尾砂粒径效应、坡比、含水率以及干滩的长度变化。试验成功地模拟了尾矿库筑坝过程普遍存在的尾砂沉积现象,如

图5 不同浓度尾砂堆坝形成的干滩面
Fig. 5 Dry beach surface formed by tailing dam with different concentration
尾矿浓度是影响沉积规律的重要因素之一,如

图6 不同浓度尾砂坝前冲击效果
Fig. 6 Front impact effect of different concentration tailings dam
在每组排浆浓度沉积模型试验中,对放矿结束后尾砂堆积体形态、干滩坡度、干滩面长度和水面长度及不同流距干滩含水率进行测量,由图

图7 不同浓度浆体形成干滩坡度
Fig. 7 Different concentrations of slurry form dry beach slope

图8 不同浓度干滩面和水面长度变化曲线
Fig. 8 Curve of length of dry beach surface with different concentration

图9 不同流距干滩含水率
Fig. 9 Dry beach moisture content at different flow distances
根据试验数据得到7种不同浓度下中值粒径d50位置的变化情况,如

图10 不同浓度尾砂平均粒径沿干滩面的变化规律曲线
Fig. 10 Variation of average particle size of tailings with different concentrations along the dry beach surface
尾砂颗粒在竖直向的沉积特性对干滩面固结强度和尾矿库的稳定性评价有着重要的指导意义。

图11 不同浓度尾矿筑坝剖面图
Fig. 11 Profile of tailings dam construction with different concentrations
不同范围颗粒粒径在不同深度的含量变化趋势如

图12 不同粒径在不同深度的含量变化趋势
Fig. 12 The content variation trend of different particle sizes at different depths
干滩面尾砂颗粒分布规律测量通过现场取样,并用激光颗粒分析仪进行粒径分析。本次测量从子坝开始沿干滩面往库区内选4个典型放矿支管进行纵向取样,每条纵向取样线长约45 m,每个尾砂样间隔1.5 m,纵向共取样120个。在距离坝前30 m处进行横向取样,取样的起始点涵盖4个放矿支管,横向取样采取每隔3 m进行,横向共取样35个尾矿试样,如

图13 尾矿库干滩面取样
Fig. 13 Tailing pond dry beach surface sampling
堆积尾砂沉积规律受尾矿性质、粒度、矿浆浓度和排放形式控制。根据勘察结果和各项控制因素的分析及对室内土工试验的统计分析,对堆积尾砂的沉积规律认识如下。
当尾矿浆液由坝顶向库尾流动时,根据水动力规律,随着浆液的流动,水力梯度逐渐消散,浆液所能携带尾矿的能力逐渐降低,从而尾矿逐渐按照颗粒由粗到细依次沉淀下来,从而在干滩面上,形成了坝前粗库尾细的沉积特点。根据本次勘察结果,对位于同一高程面土试样的平均粒径d50进行统计分析,并绘制成平均粒径d50的变化曲线,从

图14 四条取样线平均粒径变化曲线
Fig. 14 Average particle size change curve of four sampling lines
当尾矿浆由坝顶向库内流动时,随着水力梯度的消散,每层平面上尾矿逐渐按照颗粒由粗到细沉淀下来,随着堆积坝的升高,坝顶不断向上游推移,对于每一固定的纵断面来说,其浆液流动的里程在逐段减小,浆液经过该断面时所能携带的尾矿颗粒的粒径在逐渐增大。因此,随坝体加高,子坝向上游的推移,该断面上即呈现出上粗下细的沉积特性。根据现场勘察成果,选取部分钻孔分别对其不同深度的土试样的平均粒径d50进行统计分析,并绘制成平均粒径d50的变化曲线,从

图15 不同钻孔土试样平均粒径d50随孔深变化曲线
Fig. 15 The average particle size d50 of different drilled soil samples varies with hole depth
1)不同排浆浓度尾砂流动沉积过程中,干滩面的沉积尾砂在纵向上表现为前粗后细的特点,而在竖向上具有上粗下细的沉积特性,高浓度尾砂在筑坝过程中由于浓度较高自身拥有的粘滞性较大与沉积滩面之间产生较大摩擦阻力,很难被冲刷重新启动,从而出现明显的分层现象,且低浓度尾砂矿浆排放沉积的交错夹层结构现象更为显著。
2)不同排浆浓度尾砂流动沉积过程中,高浓度尾矿矿浆形成的干滩面坡度明显大于低浓度矿浆形成的干滩面坡度,高浓度矿浆排放形成的坡度一般在2.00%~3.50%,而低浓度尾砂矿浆排放形成的坡度较平缓,为0.15%~2.00%。
3)不同排浆浓度尾砂流动沉积过程中,不同浓度尾砂沉积粒径总体表现为排浆浓度越小粒径沿沉积距离变化越明显,分选性越好,干滩面长度越短,高浓度筑坝干滩面长度较长,尾水面较短。
4)分析了单颗尾砂颗粒的临界状态,推导了尾砂颗粒运动—沉积临界速度关系和不同粒径颗粒在干滩上的运移距离公式。推导出不同浓度尾砂筑坝在干滩面的运移距离,尾矿颗粒在干滩面上的沉积位置与排浆的初始流速成正比,和干滩坡度、排浆浓度及颗粒粒径等成反比。
参考文献
尹光志, 敬小非, 魏作安, 等. 粗、细尾砂筑坝渗流特性模型试验及现场实测研究[J]. 岩石力学与工程学报, 2010, 29(S2): 3710-3718. [百度学术]
Yin G Z, Jing X F, Wei Z A, et al. Study of model test of seepage characteristics and field measurement of coarse and fine tailings dam[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(S2): 3710-3718. (in Chinese) [百度学术]
敬小非, 尹光志, 魏作安, 等. 模型试验与数值模拟对尾矿坝稳定性综合预测[J]. 重庆大学学报(自然科学版), 2009, 32(3): 308-313. [百度学术]
Jing X F, Yin G Z, Wei Z A, et al. Model test and numerical simulation of tailing dam safety forecasting[J]. Journal of Chongqing University (Natural Science Edition), 2009, 32(3): 308-313. (in Chinese) [百度学术]
巫尚蔚, 杨春和, 张超, 等. 干滩表层沉积尾矿的细观几何特征研究[J]. 岩石力学与工程学报, 2016, 35(4): 768-777. [百度学术]
Wu S W, Yang C H, Zhang C, et al. Microscopic geometric characteristics of surface sedimentary tailings[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(4): 768-777. (in Chinese) [百度学术]
史采星, 郭利杰, 杨超, 等. 某铜镍矿尾矿流变参数测试及管道输送阻力计算[J]. 中国矿业, 2018, 27(B10): 138-141. [百度学术]
Shi C X, Guo L J, Yang C, et al. Experimental study of the rheological parameters of a copper nickel mine tailings and calculation of resistance in pipeline transportation[J]. China Mining Magazine, 2018, 27(B10): 138-141. (in Chinese) [百度学术]
李宗楠, 郭利杰, 余斌, 等. 基于宾汉姆体的高浓度尾砂浆剪切变稀规律研究[J]. 黄金科学技术, 2017, 25(4): 33-38. [百度学术]
Li Z N, Guo L J, Yu B, et al. Shearing thinning behavior of high concentration slurry based on Bingham model[J]. Gold Science and Technology, 2017, 25(4): 33-38. (in Chinese) [百度学术]
张东明, 郑彬彬, 尹光志, 等. 高浓缩分级尾矿上游法堆坝及模型试验研究[J]. 岩土力学, 2016, 37(7): 1832-1838, 1867. [百度学术]
Zhang D M, Zheng B B, Yin G Z, et al. Model tests on upstream dam-building method using concentrated and classified tailings[J]. Rock and Soil Mechanics, 2016, 37(7): 1832-1838, 1867. (in Chinese) [百度学术]
唐永俊, 王光进, 许志发, 等. 尾矿库干滩尾矿粒径的淤积特性研究[J]. 泥沙研究, 2018, 43(3): 50-56. [百度学术]
Tang Y J, Wang G J, Xu Z F, et al. Study on deposition characteristics of tailings particle size in dry beach of tailings pond[J]. Journal of Sediment Research, 2018, 43(3): 50-56. (in Chinese) [百度学术]
赵怀刚, 王光进, 许志发, 等. 多因素影响下尾矿库坝体沉积特性的试验研究[J]. 中国安全生产科学技术, 2018, 14(5): 95-101. [百度学术]
Zhao H G, Wang G J, Xu Z F, et al. Experimental study on deposition characteristics of tailings dam body influenced by multiple factors[J]. Journal of Safety Science and Technology, 2018, 14(5): 95-101. (in Chinese) [百度学术]
Jopony M, Usup G, Mohamed M. Particle size distribution of copper mine tailings from lohan ranau Sabah and its relationship with heavy metal content[J]. Pertanika, 1987, 10(1): 37-40. [百度学术]
Fitton T G, Bhattacharya S N, Chryss A G. Three-dimensional modeling of tailings beach shape[J]. Computer-Aided Civil and Infrastructure Engineering, 2008, 23(1): 31-44. [百度学术]
Kwak M, James D F, Klein K A. Flow behaviour of tailings paste for surface disposal[J]. International Journal of Mineral Processing, 2005, 77(3): 139-153. [百度学术]
Henriquez J, Simms P. Dynamic imaging and modelling of multilayer deposition of gold paste tailings[J]. Minerals Engineering, 2009, 22(2): 128-139. [百度学术]
Babaoglu Y, Simms P H. Simulating deposition of high density tailings using smoothed particle hydrodynamics[J]. Korea-Australia Rheology Journal, 2017, 29(3): 229-237. [百度学术]
余绍维, 孙华, 王兆昌, 等. 上游法尾矿坝堆积尾矿的沉积规律与尾矿堆积坝管理的相关关系[J]. 中国新技术新产品, 2013(6): 175-177. [百度学术]
Yu S W, Sun H, Wang Z C, et al. Correlation between the deposition law of tailings accumulated in upstream tailings dam and the management of tailings accumulation dam[J]. China New Technologies and New Products, 2013(6): 175-177. (in Chinese) [百度学术]
梁冰, 吕志强, 金佳旭, 等. 排浆速度对尾矿沉积影响的模型试验研究[J]. 实验力学, 2017, 32(6): 880-887. [百度学术]
Liang B, Lv Z Q, Jin J X, et al. Model experimental study of effect of slurry draining speed on tailings deposition[J]. Journal of Experimental Mechanics, 2017, 32(6): 880-887. (in Chinese) [百度学术]
Huang S H, Li W, Cheng L J. On equation of discrete solid particles' motion in arbitrary flow field and its properties[J].Applied Mathematics and Mechanics, 2000, 21(3): 297-310. [百度学术]
郑艳娜, 朱永英. 水力学[M]. 南京: 东南大学出版社, 2017. [百度学术]
Zheng Y N, Zhu Y Y. Hydraulics[M]. Nanjing: Southeast University Press, 2017. (in Chinese) [百度学术]