Experimental study on characteristics of impact force of tailing flow under dam break of tailing reservoir
-
摘要: 尾矿库溃坝下泄尾砂流冲击力的计算目前还未有成熟的模型或计算方法,而直接采用泥石流的经验公式或现有的水力学计算模型并修改相应的参数不能真实反映尾砂流的冲力特性。采用尾砂流模型槽装置,开展了不同密度的尾砂浆体、不同流速、不同埋深条件下冲击力的演化特征。试验结果表明:冲击力在时间分布上受密度的影响规律为,密度越大最大冲击力出现的时间越早;冲击力在纵向空间上可分为两段,一段为液面以下的埋深范围内,冲击力随深度呈线性分布,另一段为涌高部分,冲击力呈对数型快速衰减;尾砂流冲击力与流速、密度和埋深直接相关,不同密度下流速与冲击力呈幂函数相关,同时密度越大冲击力越大,呈幂函数相关性,埋深与冲击力呈线性关系。并以此建立了涵盖了流速、密度和埋深3因素的冲击力模型,并设定了3个参数,分别为扰流系数k1受密度影响的冲击力调整系数α和受埋深对冲击力的调整系数k2,该模型能综合反映3因素的影响,通过验证表明了模型的合理性及参数物理意义的正确性,可为相关研究和工程应用提供一定的理论参考。Abstract: At present, there is no mature model or method for the calculation of the impact force of tailing flow discharged from the dam break of tailing pond. Instead, the empirical formula for debris flow or the existing hydraulic model are directly used, and the corresponding parameters are modified. This method can not truly reflect the impact characteristics of tailing flow. The evolution characteristics of impact force under different densities of tailing slurry, different velocities and different depths of impact are studied by using the tailing flow model groove devices. The test results show that the impact force is affected by density shadow in time distribution. The higher the density is, the earlier the maximum impact force appears. In the longitudinal space, the impact force can be divided into two parts. One is the buried depth below the liquid level, the impact force is linearly distributed with the depth, the other is the surge height part. The impact force is a logarithmic fast attenuation. The impact force of tailing flow is directly related to velocity, density and buried depth, and the velocity and impact force are power function-related under different densities. At the same time, the greater the density is, the greater the impact force is the power function related. The buried depth and impact force are linear relationship. Based on this, a model for impact force covering three factors of the velocity, density and buried depth is established, and three parameters are set up, namely, the coefficient of turbulence k1, the coefficient of impact force adjustment α affected by the density, and the coefficient of impact force adjustment k2 affected by the buried depth. The model can comprehensively reflect the influences of three factors. The rationality of the model and the correctness of the physical meaning of the parameters are verified, which may provide some theoretical reference for the related researches and engineering applications
-
Keywords:
- tailings reservoir /
- dam break /
- tailing flow /
- impact force /
- model test
-
0. 引言
近年来,国内外学者针对经历循环荷载历史后土体的剪切特性进行了大量研究[1-5]。一般来说,循环荷载会使土体的超孔隙水压力增大,平均有效应力降低,从而导致土体的循环后强度降低。王淑云等[1]针对重塑粉质黏土,在不同围压下进行了一系列静三轴和动-静三轴不排水试验,发现粉质黏土的振后不排水强度衰减程度取决于动载引起的动应变和孔压值;郑刚等[2]认为原状土样的振后不排水抗剪强度显著衰减,而重塑土的抗剪强度变化不甚明显;Moses等[3]发现土体的振后不排水强度随循环荷载幅值的增加而减小;Yasuhara等[5]研究了重塑Ariake黏土的循环后剪切特性,发现循环荷载作用后土体累积孔压可以很好地评价振后不排水剪切强度,并提出了预测土体振后抗剪强度的经验模型。
然而,经历循环荷载作用后的土体往往处于未固结和完全固结的中间状态,该状态可用振后固结度这一概念进行描述。目前针对原状软土在不同振后固结度条件下的剪切特性研究较少。因此,有必要针对不同振后固结度下的软土振后剪切特性进行研究。本文选取珠江入海口原状软土,通过一系列的动-静三轴试验,分析了初始围压,循环应力比以及振后固结度对土体振后剪切特性的影响,以期加深对软土振后静力特性的理解,为工程设计提供试验基础。
1. 土体试样和试验方案
1.1 土样基本物理力学特性
试验选取珠江入海口原状软土为研究对象,土壤呈深灰色,取土深度为0~7 m。依据《土工试验规程(GB/T 50123—2019)》获得土体的基本物理力学指标和土体的粒径分布曲线分别见表1和图1所示。
表 1 原状软土基本物理力学指标值Table 1. Indices of basic physical and mechanical properties of undisturbed soft soils物理力学特性 取值 天然密度ρ/(g·cm-3) 1.54~1.79 含水率w/% 38.3~69.1 孔隙比e 1.12~1.91 相对质量密度Gs 2.71 液限wL/% 38.5~61.8 塑限wP/% 19.4~28.0 塑性指数IP 19.1~33.8 渗透系数K/(10-7cm·s-1) 8.79~10.80 压缩系数av/MPa-1 0.48~1.15 侧压力系数K0 0.44 1.2 试验方案
通过薄壁取土器进行现场取土,并将原状土样保存在恒温恒湿箱中。按照《土工试验方法标准(GB/T 50123—2019)》的要求,将土体制成直径38 mm,高76 mm的原状试样,并采用真空饱和及反压饱和法对试样进行饱和。首先对所有圆柱试样进行真空饱和,然后将试样置于压力室内进行反压饱和,当B值达到0.95以上时,可认为土体已经饱和。随后,对试样施加一定的固结压力,进行等向固结,当试样排水体积速率小于100 mm3/h时,认为土样固结完成。为研究有效固结围压对振后强度的影响,本次试验中选取有效固结围压分别为20,40,60 kPa。
对固结完成后的土体在不排水状态下开展循环三轴试验,并采用Sakai等[6]提出的循环应力比CSR描述循环偏应力大小,即
CSR=q/2p′0=q/2σ′3, (1) 式中,q为循环偏应力幅值,
p′0 为固结完成后的有效围压。为研究CSR对振后土体强度的影响,在围压为60 kPa下,CSR值分别取0.08,0.17,0.25,0.33。循环轴向偏应力采用半正弦波形,频率为0.1 Hz,且循环振次为1000次。对振动后的试样再次进行固结过程,并采用振后固结度这一概念描述振后土体固结程度。试验过程中为得到不同振后固结度的土体,可通过向振后试样施加不同的反压。具体为:向经历循环荷载作用后的土体施加一定的反压pu,此时,振后土体在一定围压作用下进行固结,振动过程中产生的超孔隙水压力将逐渐降低直至与反压相等。因此,振后不同固结度Ur可通过下式计算得到,
Ur=1−pu(Δu)cy, (2) 式中,pu为再固结过程中施加在试样上的反压,(∆u)cy为振动过程中产生的超孔隙水压力。
当试样再次固结完成后,对土体进行三轴固结不排水剪切试验。其中,剪切过程中均采用应变控制,剪切速率为0.1%/min。当应变达到20%时,试验结束。对于经历和未经历循环荷载作用的土体而言,若静力剪切过程中应力应变曲线出现峰值点,其不排水抗剪强度为峰值偏应力的1/2;若未出现峰值点,则土体不排水抗剪强度取20%应变处对应偏应力的1/2。整个试验方案如表2所示,其中ST-20,ST-40和ST-60为未经历循环荷载作用的土体在不同围压条件下的静三轴试验。经历和未经历循环荷载作用的土体强度用Su表示。
表 2 试验方案Table 2. Test schemes试样编号 p′0 /kPaqampl/kPa CSR 循环次数 Ur/% Su/kPa U01 20 10 0.25 1000 0 20.2 U02 20 10 0.25 1000 25 24.4 U03 20 10 0.25 1000 75 22.8 U04 20 10 0.25 1000 100 33.0 U05 40 20 0.25 1000 100 40.4 U06 60 30 0.25 1000 100 39.5 U07 60 10 0.08 1000 100 32.7 U08 60 20 0.17 1000 100 31.5 U09 60 40 0.33 1000 100 42.9 ST-20 20 — — — — 28.5 ST-40 40 — — — — 33.4 ST-60 60 — — — — 33.5 2. 试验结果
2.1 固结围压对振后土体强度的影响
图2展示了循环荷载后完全固结试样在不同围压下的偏应力-应变曲线。其中,CSR=0.25,Ur= 100%。可以看出,不同围压下的q-ε曲线变化趋势相同。当应变较低时,应力在小应变范围内迅速增加。随着应变的增加,应力增速放缓,当轴向应变为20%左右时,偏应力趋于稳定。不同围压下完全固结的振后试样强度大于相同围压无循环荷载作用历史的土体静强度。例如,当围压为60 kPa时,试样的振后强度比未经历循环荷载作用下的强度大6.0 kPa。此外,对比发现,实验条件相同的情况下,振后土体的抗剪强度随着围压的增加基本呈增大趋势。具体为完全固结的振后试样在初始固结围压为20,40,60 kPa时,对应的抗剪强度分别为33.0,40.4,39.5 kPa。
2.2 循环应力比对振后土体强度的影响
图3展示了不同循环应力比循环荷载作用后完全固结试样(Ur=100%)的偏应力-应变曲线。其中,初始固结围压为60 kPa。可以看出,初始应变随着CSR的增加而增加。当CSR值分别为0.08,0.17,0.25和0.33时,对应的初始应变分别为0.24%,0.43%,0.66%,4.27%。当CSR>0.25时,试样的振后剪切强度明显大于静态剪切强度(表2)。相较于未经历循环荷载的土体静强度,当CSR=0.08时,振后土体强度增幅为-0.8 kPa;当CSR=0.33时,振后土体强度增幅为9.4 kPa。同样地,土体振后剪切强度随着CSR的增大而增大。当CSR从0.08增大到0.33时,土体振后剪切强度从32.7 kPa增大到42.9 kPa。
2.3 振后固结度对振后土体强度的影响
图4为不同振后固结度下土体的偏应力-应变曲线。其中,初始固结围压为20 kPa,CSR=0.25。通过与未经历动荷载的试样q-ε曲线对比发现,Ur在75%~100%存在临界值,使得振后土体强度等于相同围压下未经历动荷载作用的土体强度。当固结围压为20 kPa,固结度为0%,25%,75%,100%时,抗剪强度分别为20.2,24.4,22.8,33.0 kPa。试验结果表明土体的振后抗剪强度随Ur的增大而增加。
3. 结论
本文通过对珠江入海口原状软土开展一系列动、静三轴试验,对其振后强度特性进行了研究。
(1)经历循环荷载作用后的土体,不同因素影响下对应的偏应力-应变曲线变化趋势基本一致。不同固结围压影响下,振后完全固结试样的振后抗剪强度大于未经历动荷载土体静强度。另一方面,土体振后抗剪强度随初始固结围压的增大呈增加趋势。
(2)当循环应力比大于0.25时,振后完全固结状态下土体剪切强度明显大于相同围压下未经历循环荷载作用的土体静态剪切强度。另一方面,振后完全固结的土体,其振后剪切强度随着循环应力比的增大而增加。
(3)土体振后抗剪强度随振后固结度的增加而增大,且振后固结度存在某一临界值,使得土体的振后抗剪强度与相同围压条件下未经历循环荷载作用的土体抗剪强度基本一致。
-
表 1 尾矿砂粒径累计含量
Table 1 Accumulated contents of characteristic particle size of tailings sand
(%) 包线类型 颗粒粒径/mm 2 0.5 0.25 0.075 0.005 0.002 综合细包线 100 100 75 24 17 7 综合粗包线 97 84 50 15 10 2 综合平均线 100 90 57 20 12 5 表 2 尾砂流试验控制参数
Table 2 Control parameters of tailings flow tests
密度/(g·cm-3) 水砂比 含砂量/(g·cm-3) 重度/(kN·m-3) 流速/(m·s-1) 最大埋深/cm 1.0 1∶0 0 9.80 0.5~3 20 1.1 5.58∶1 0.17 10.78 0.2~1.5 20 1.2 2.59∶1 0.33 11.76 0.2~1.5 20 1.3 1.59∶1 0.50 12.74 0.2~1.5 20 1.4 1.09∶1 0.67 13.72 0.2~1.5 20 1.5 0.80∶1 0.84 14.70 0.2~1.5 20 表 3 尾砂流冲击特性试验结果
Table 3 Test results of impact characteristics of tailings flow
密度/(g·cm-3) 流速/(m·s-1) 动压力/kPa 5 cm 10 cm 15 cm 20 cm 1.1 0.30 0.33 0.78 1.23 1.82 0.32 0.45 0.86 0.96 1.58 0.39 0.52 1.08 1.07 1.62 0.50 0.61 1.12 1.22 1.68 0.56 0.72 1.27 1.27 1.93 0.70 1.00 1.42 1.52 2.17 0.90 1.51 2.09 2.21 2.71 1.10 2.22 2.75 2.87 3.39 1.15 2.43 2.98 3.09 3.65 1.2 0.23 0.48 0.69 1.13 2.17 0.30 0.64 0.94 1.41 2.15 0.33 0.83 1.13 1.99 2.52 0.50 1.12 1.66 2.10 2.78 0.65 1.44 2.02 2.35 3.10 0.71 1.75 2.25 2.88 3.46 0.90 2.22 2.83 3.16 3.88 1.11 2.97 3.48 3.90 4.72 1.24 3.62 4.75 4.58 5.40 1.3 0.30 0.81 1.09 1.37 2.72 0.50 1.49 1.89 2.20 3.35 0.70 2.11 2.75 3.20 3.68 0.79 2.32 2.94 3.48 4.11 0.91 2.92 3.39 4.18 4.62 1.11 3.77 4.20 5.01 5.58 1.21 4.61 5.20 6.49 6.25 1.4 0.22 0.85 1.24 1.56 1.77 0.30 1.17 1.67 1.94 2.11 0.37 1.50 2.06 2.23 2.67 0.50 2.22 2.65 3.09 3.32 0.69 3.27 3.95 4.54 4.98 1.5 0.30 1.48 2.10 2.69 3.04 0.50 2.76 3.46 4.28 4.54 0.61 3.17 3.91 4.67 6.08 0.71 4.21 4.75 5.43 6.92 0.93 5.97 6.55 7.47 8.58 1.05 6.96 7.58 8.06 9.87 1.12 8.28 8.81 9.50 11.97 表 4 冲击模型拟合结果取值
Table 4 Values of fitting results impact model
密度ρ/(g·cm-3) 埋深h/m k1 α k2 相关系数 1.1 0.05 1.65 7.24 0.69 0.99 1.2 0.10 1.80 5.47 0.69 0.97 1.3 0.15 1.92 4.81 0.66 0.96 1.4 0.20 2.43 4.28 0.69 0.97 1.5 0.20 2.65 4.07 0.67 0.97 表 5 文献实验结果与模型计算结果对比
Table 5 Comparison between model results and experimental data of Reference [12]
序号 试验数据 密度/(g·cm-3) 流速/(m·s-1) 平均埋深/m 冲击力/kPa 1 1.7 3.85 0.09 65.48 2 1.8 3.74 0.09 59.14 3 1.9 3.57 0.09 64.75 4 2.0 3.05 0.09 63.15 5 2.1 2.78 0.09 67.75 序号 模型参数及结果 k1 α k2 冲击力/kPa 1 2.81 2.17 0.68 66.91 2 2.91 1.83 0.69 60.79 3 3.04 1.81 0.68 63.07 4 3.25 2.05 0.68 63.82 5 3.38 2.21 0.68 68.60 -
[1] 张家荣, 刘建林, 朱记伟. 我国尾矿库事故统计分析及对策建议[J]. 武汉理工大学学报(信息与管理工程版), 2016, 38(6): 682-685. https://www.cnki.com.cn/Article/CJFDTOTAL-WHQC201606009.htm ZHANG Jia-rong, LIU Jian-lin, ZHU Ji-wei. Statistical analysis and countermeasures of tailings accident in China[J]. Journal of Wuhan University of Technology, 2016, 38(6): 682-685. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-WHQC201606009.htm
[2] 吴宗之, 梅国栋. 尾矿库事故统计分析及溃坝成因研究[J]. 中国安全科学学报, 2014, 24(9): 70-76. https://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK201409014.htm WU Zong-zhi, MEI Guo-dong. Statistical analysis of tailings pond accidents and cause analysis of dam failure[J]. China Safety Science Journal, 2014, 24(9): 70-76. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZAQK201409014.htm
[3] 闪淳昌, 张振东, 钟开斌, 等. 襄汾“9·8”特别重大尾矿库溃坝事故处置过程回顾与总结[J]. 中国应急管理, 2011(10): 13-18. https://www.cnki.com.cn/Article/CJFDTOTAL-YIGU201110008.htm SHAN Chun-chang, ZHANG Zhen-dong, ZHONG Kai-bin, et al. Review and summary of the disposal process of Xiangfen "9.8" extremely serious tailings dam break accident[J]. China Emergency Management, 2011(10): 13-18. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YIGU201110008.htm
[4] 殷跃平. 山西襄汾县塔山矿区尾矿溃决泥流灾难[J]. 中国地质灾害与防治学报, 2008, 19(4): 70-148. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGDH200804014.htm YIN Yue-ping. Debris flow disaster caused by tailings outburst in Tashan mining area, Xiangfen County, Shanxi Province[J]. The Chinese Journal of Geological Hazard and Control, 2008, 19(4): 70-148. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGDH200804014.htm
[5] HU Kai-heng, CUI Peng, GE Yong-gang. Building destruction patterns by august 8, 2010 Debris flow in Zhouqu, Western China[J]. Journal of Mountain Science, 2012, 30(4): 484-490. doi: 10.3969/j.issn.1008-2786.2012.04.015
[6] ZOU Qiang, GUO Xiao-jun, ZHU Xing-hua, et al. Hazard characteristics and causes of “7·10 ” debris flow along highways in the upper reaches of minjiang riverhazard characteristics and causes of “7·10 ” debris flow along highways in the upper reaches of minjiang river[J]. Journal of Mountain Science, 2014(6): 747-753.
[7] IVERSON R M. Scaling and design of land slide and debris-flowexperiments[J]. Geomorphology, 2015, 244: 9-20.
[8] 尹光志, 敬小非, 魏作安, 等. 尾矿坝溃坝相似模拟试验研究[J]. 岩石力学与工程学报, 2010, 29(增刊2): 3830-3838. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2010S2056.htm YIN Guang-zhi, JING Xiao-fei, WEI Zuo-an, et al. Experimental study of similar simulation of tailings dam-break[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(S2): 3830-3838. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2010S2056.htm
[9] JING F X. Research on the movement mechanism of pulp surging from tailings dam-break[J]. Advanced Materials Research, 2012, 516-517: 997-1000.
[10] WU T, QIN J. Experimental study of a tailings impoundment dam failure due to overtopping[J]. Mine Water & the Environment, 2018, 37(2): 272-280.
[11] LIU L, ZHANG H W, ZHONG D Y, et al. Research on tailings dam break due to overtopping[J]. Journal of Hydraulic Engineering, 2014, 45(6): 675-681.
[12] 刘道川, 游勇, 杜杰, 等. 泥石流冲击力的时空分布特征[J]. 工程科学与技术, 2019, 51(3): 17-25. https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201903003.htm LIU Dao-chuan, YOU Yong, DU JIE, et al. Spatio-temporal distribution of the impact force of debris flow[J]. Advanced Engineering Ences, 2019, 51(3): 17-25. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201903003.htm
-
期刊类型引用(13)
1. 杨东东. 基于孔压静力触探的土层分类研究. 中国新技术新产品. 2024(04): 102-104 . 百度学术
2. 马昌慧,陈义平,郭峰,李栋广. 浅析基坑工程的优化设计方法. 工程勘察. 2024(04): 26-30 . 百度学术
3. 魏军. 大型自升式海工平台桩靴地基土工程特性试验研究. 江苏建筑. 2024(02): 96-100+120 . 百度学术
4. 李伟. 软基处理PST桩复合地基加固机理及CPTU验证. 山西建筑. 2024(11): 76-80+85 . 百度学术
5. 时文峰,丁川,仇涛,童立元,王逸凡,马海洋. 地下工程预制梁钢筋骨架安装与成型技术研究. 建筑技术. 2024(16): 1924-1928 . 百度学术
6. 肖乾,袁颖,汪婷. 基于系统聚类算法的CPT土体分层试验研究. 湖北理工学院学报. 2024(06): 36-41 . 百度学术
7. 卢云龙,戴正彬,汤国毅. 基于静力触探对深基坑突涌区地层扰动的分析. 土工基础. 2023(01): 157-160 . 百度学术
8. 丁剑桥. 基于静力触探的桩基承载力及数值模拟分析. 安徽建筑. 2023(03): 135-136+173 . 百度学术
9. 吴早生,胡春东,白浩东,沙鹏,黄曼,王天佐. 静力触探深度与垂直度实时智能监测系统的应用研究. 中国新技术新产品. 2023(04): 90-93 . 百度学术
10. 曹振平,曾彪,张雨波,孟韬,周盛生,谷军. 基于CPTU的长江河漫滩相软土夹层分布特征与渗透特性评价. 勘察科学技术. 2023(03): 43-48+64 . 百度学术
11. 吴早生,白浩东,胡春东,沙鹏,黄曼,张鑫,王天佐. 基于静力触探数据的绍兴平原地区土层参数分析评价. 岩土工程技术. 2023(06): 692-699 . 百度学术
12. 郭全元,钟仕兴,肖旺. 基于孔压静力触探(CPTU)及扁铲侧胀试验(DMT)南沙软土力学特性研究. 广东土木与建筑. 2022(06): 26-29 . 百度学术
13. 周燕国,王智男,马强,梁甜,陈云敏. ZJU-400离心机机载微型CPTu装置研制与性能测试. 岩土工程学报. 2022(S2): 11-15 . 本站查看
其他类型引用(4)