Experimental study on infiltration of typical ion-adsorption rare-earth ore in southern Jiangxi
-
摘要: 选择原地浸矿过程中表层滑坡频发矿区,现场测试2个矿块的稀土矿层入渗过程,采用4种常用入渗模型拟合入渗过程.结果表明:整个入渗过程可分为3个阶段:0~45 min为入渗率急剧变化阶段,45~125 min为入渗率缓慢变化阶段和125 min之后逐渐达到稳渗阶段;用Kostiakov模型拟合的效果好于其他3种模型,矿块1初始入渗率平均为0.113 cm/min,稳定入渗率为0.063 cm/min,矿块2初始入渗率平均为0.162 cm/min,稳定入渗率为0.101 cm/min;入渗率随砂砾含量增大而增大,随黏粒含量增大而减小,与矿层含水量呈负相关关系;试验矿区渗透系数随深度增加快速减小,注液强度过大而远超渗透系数时,造成表层矿层溶液无法及时下渗,边坡安全系数减小,易引发表层滑坡.Abstract: Selecting the mining area with frequent surface landslide in the process of in-situ leaching, the infiltration rate of rare earth from two ore blocks in filed was tested based on four common models fitting the water infiltration characteristic. The results show that the whole infiltration process can be divided into three phases: the infiltration rate of 0~45 min is the rapid change stage, the infiltration rate of 45~125 min is slowly change phase and after 125 min the infiltration rate gradually achieves the stability. The fitting effect of Kostiakov is better than the other three models with the average ore initial infiltration rate being 0.113 cm/min, the steady infiltration rate of 0.063 cm/min in the first sample region, the average ore initial infiltration rate of 0.162 cm/min and the steady infiltration rate of 0.101 cm/min in second sample region. The infiltration rate increases with the rise of sand and gravel content and decreases with the rise of clay content. At the same time, the infiltration rate is in the negative correlation with the ore moisture content. Hydraulic conductivity in vertical direction shows a rapid decreasing trend with the increase of depth. When liquid injection intensity is too excessive and far more than the hydraulic conductivity, surface seam liquid can not infiltrate in time, which will lead to the decrease of slope safety factor and the occurrence of surface landslide.
-
Keywords:
- surface landslide /
- infiltration model /
- hydraulic conductivity /
- in-situ leaching
-
0 引言
许多国内外学者针对控制着工程岩体稳定性的岩体结构面力学特性这一热点问题进行了持续性研究.众所周知,结构面的存在对岩体力学性质有着决定性影响,因此,研究岩体结构面力学性质是非常有意义的.
孙广忠教授在文献[1]中归纳了结构面的4 种起伏形态分别为平直结构面、台阶状结构面、锯齿状结构面、 波状结构面,4 种结构面的力学效应是不同的. 除了锯齿状结构面,其他3 种结构面力学特性的研究已较多,一些学者对岩体结构面锚固特性采用室内试验或理论分析进行了相关研究[2-4].但是对结构面注浆加固前后力学特性的研究还比较少,尤其是采用数值分析方法的研究更少.因此,本文通过3 种类型和2 种台阶高度结构面剪切模拟试验,对结构面注浆前后抗剪特性进行了分析研究.该研究对揭示其加固机制具有重要的意义.
1 试验与计算模型
通过ANSYS 有限元分析软件建立数值分析模型,再通过转换导人FLAC3D 软件中进行分析和计算[5-12].为尽可能得到岩石结构面剪切特性的规律,模型尺寸:长×宽×高为100 mm×100 mm×100 mm.为了模拟结构面,在模型中间设置厚度为2mm 的结构面模拟其注浆前后效果. 设置的结构面类型包括1个台阶、2个台阶以及组合(1个台阶与1 个锯齿)3 种类型,台阶宽度为20 mm,每种类型包括2 种台阶高度,分别为10 mm、6 mm,则3 种类型结构面模型高宽比分别为0.5 和0.3. 如图1-图3 所示分别为3 种结构面模型示意图.
边界条件设置为底部固定约束,下半部分模型四周约束水平方向位移,上半部分顶部为自由边界.
计算采用Mohr-Coulomb准则,在模型顶部施加法向荷载,上半部分岩石施加水平位移荷载,荷载大小2×l0-4 mm/时步. 模拟试验中只考虑法向荷载,因为法向荷载远远大于岩石自重,所以计算略去岩石自重[13-15].如表1、表2 所示,岩石和结构面参数是根据RMT-150C 型岩石力学试验系统对岩石试块进行室内单轴和剪切试验的结果选取.
表 1 岩石参数选取表 2 结构面参数选取2 破坏规律分析
通过数值计算分析讨论台阶状结构面剪切破坏形式以及规律,探讨台阶数量及结构面形态与抗剪强度关系、施加法向应力大小与抗剪强度关系.由于试验结果图太多,此处只列出图4~图6,分别为试验中高宽比0.5 和0.3 时,结构面最大、最小主应力云图以及组合剖面塑性云图. 本文中所有图示除非特别说明,剪切方向均为从左往右.
通过分析图4、图5、图6 等1 个台阶、2 个台阶、 组合结构面模型剖面的最大主应力、最小主应力以及塑性区云图得到下述研究结果.
2.1 台阶高宽比为0.5 情况分析
1) 拉应力区域皆形成于3 种结构面的台阶及锯齿下部,拉应力最大处位于底部.
2) 结构面上半部分岩石在剪切作用下有往上抬起效应,导致台阶右边下部岩石受到明显压应力作用,台阶右边底角转折处形成压应力集中区域.
3) 组合类型岩石在剪切作用下锯齿下方和锯齿顶部均表现为拉应力.而压应力较集中区域为锯齿正下方与台阶底部,且在锯齿右面底部达到最大压应力,说明剪切过程上部岩石有往上抬起效应.
4) 由塑性云图6 可知,因岩石抗拉强度远远小于岩石抗压强度,导致台阶底部在剪切应力和法向应力的作用下最先表现为拉伸破坏.组合类型模型中锯齿和台阶的底部都已开始出现拉伸破坏,且最大拉应力区域与破坏区域一致,台阶最终因张剪滑移而破坏.
2.2 台阶高宽比为0.3 情况分析
为了对比分析台阶高宽比不同的情况下3 种结构面剪切破坏规律,进行了台阶高宽比为0.3 时结构面剪切试验.
通过分析台阶高宽比为0.3 的1 个台阶结构面模型剖面的最大主应力、最小主应力以及塑性区云图可知:
1) 台阶高度降低后,台阶高宽比为0.3 时结构面主应力云图与台阶高宽比为0.5 时的应力云图规律是一样的.台阶底部受到拉应力作用,拉应力集中区域出现在台阶右下角.
2) 沿顺剪方向,台阶只在左下部分出现斜向上的拉伸破坏痕迹,表明当台阶高宽比较小时,台阶一般因剪切作用产生压切破坏,且对于2 个台阶结构面模型压切破坏同时发生于2 个台阶.因此,当台阶高宽比较小时,台阶因剪切作用而产生压切破坏.破坏区域分布与文献[1]分析一致.
3) 组合结构面模型中拉应力集中区域位于台阶和锯齿的底部,台阶底部受弯矩作用明显,由此可见:高宽比越大,台阶状结构面张性破坏效应越明显.锯齿底部的拉应力明显小于台阶底部的拉应力;锯齿高度降低后,锯齿下部的拉应力集中效应降低了,由此可见:起伏角α 越大,结构面张性破坏效应越明显.
4) 2 个台阶结构面与1 个台阶结构面模型类似,沿受剪方向,出现于台阶左下角的拉伸破坏区逐步向台阶顶端扩展,最终导致台阶形成压切破坏.锯齿下方拉伸破坏区在左下与右下部位明显.
3 法向应力与抗剪强度分析
通过对3 种类型结构面注浆前后的试验过程监测得到了法向应力为1 MPa、2 MPa、3 MPa、4 MPa、5 MPa下峰值剪切应力的大小,如表3 所示.
表 3 不同法向应力对应峰值剪切应力/MPa度的关系特性.通过数值计算结果可知:抗剪强度和法向应力拟合相关系数近似为1.0. 结构面抗剪强度随法向应力增大而呈线性增长.
4 注浆与未注浆抗剪强度对比
依据文献4凹的研究可知,数值模拟中等效参数方法可以模拟结构面注浆效果,因此,通过改变中间层参数模拟注浆效果,结果如图7~图9 所示.
图7 所示为1 个台阶结构面注浆与未注浆结构面抗剪强度对比条形图,注浆后结构面抗剪强度在5 个等级法向应力下提高幅度明显,分别达到14.11 %、10.49 %、8.60 %、8.64 %、8.70 %.注浆后强度平均提高了10.11 %.
图8 所示为2 个台阶结构面注浆与未注浆结构面抗剪强度对比条形图,注浆后结构面抗剪强度在5 个等级法向应力下提高幅度明显,分别为11.78 %、8.33 %、 13.24 % 、6.68 % 、5.93 %. 注浆后强度平均提高了9.19%.
图9 所示为组合结构面注浆与未注浆抗剪强度对比条形图,注浆后结构面抗剪强度在5 个等级法向应力下提高幅度明显,分别为24.09 %、7.03 %、18.80 %、 16.20 %、15.70 %.注浆后强度平均提高了16.35 %.
通过对图7~图9 的分析,进一步确定岩石结构面注浆效果的等效参数模拟方法是可行的.
5 结构面强度与几何形态关系
表4 为5 个等级法向应力下1 个台阶模型和组合结构面模型抗剪强度大小关系.
表 4 不同法向应力2 种类型结构面抗剪强度比较/MPa由表4 可以看出,组合结构面模型相对于1 个台阶结构面在增加一个锯齿的情况下,结构面抗剪强度并没有明显的提高,说明大爬坡角对结构面抗剪强度大小起控制作用,定义其为控制性爬坡角;小爬坡角对结构面抗剪强度大小不起控制性作用,定义其为非控制性爬坡角.
由表4 中结构面抗剪强度增长率可以看出,组合情况下台阶(即爬坡角为900)结构面模型对结构面抗剪强度起决定性作用,锯齿对抗剪强度影响较小,平均增长率仅为2.57 %.
由此,结合前文分析,可以提出并得到组合情况下岩石结构面抗剪强度修正公式:
(1) 式(1)中: ki 为非控制性爬坡角因子,以本文中锯齿状结构面模型爬坡角为30°情况下,其非控制性爬坡角因子ki=0.025 7.
6 结论
在开展室内单轴、剪切试验的基础上,开展了结构面数值模拟试验.通过数值计算,分析了3 种类型结构面抗剪力学特性,得出结论如下:
1) 在剪切试验中,高宽比较大时,结构面一般发生剪断破坏; 高宽比较小时,结构面一般发生压切破坏.
2) 结构面抗剪强度随着法向应力的增加近似呈线性增长.
3) 通过设置中间层,采用等效参数模拟方法进行结构面注浆模拟试验是可行的. 数值模拟结果表明,注浆后结构面抗剪强度得到显著提高,从而在微观上揭示了采用注浆加固岩石结构面的可行性.
4) 以数值试验为依据提出了非控制性爬坡角因子ki 以及组合结构面抗剪强度修正公式.计算出本文中组合结构面的非控制性爬坡角因子为0.025 7.因此,可以推断出同一岩体中结构面抗剪强度大小主要由其中较大的爬坡角控制.
-
表 1 矿块1入渗过程的拟合结果
Table 1 Fitting results of infiltration process in ore block one
台阶深度/m 考斯加科夫模型 通用经验公式模型 霍顿模型 菲利普模型 均方根误差 相对差值/% 相关系数 均方根误差 相对差值/% 相关系数 均方根误差 相对差值/% 相关系数 均方根误差 相对差值/% 相关系数 1.1 7.4×10-4 9.7 0.81 3.4×10-3 18.28 0.02 2.3×10-4 6.4 0.93 4.4×10-4 14.84 0.88 1.9 8.6×10-7 2.8 0.71 1.6×10-6 3.22 0.39 2.4×10-6 4.42 0.09 1.1×10-6 2.94 0.63 2.7 9.6×10-6 4.26 0.88 4.6×10-6 2.82 0.93 1.5×10-5 4.42 0.8 2.2×10-6 2.08 0.97 3.5 1.4×10-6 1.52 0.97 1.3×10-5 3.57 0.65 1.8×10-6 1.32 0.95 3.1×10-6 1.95 0.93 4.3 1.8×10-5 6.75 0.78 1.1×10-3 32.91 / 8.7×10-5 10.11 / 2.8×10-5 8.48 0.65 5.1 8.6×10-6 5.04 0.88 2.7×10-4 22.29 / 4.2×10-5 11.12 0.31 2.2×10-5 8.39 0.69 5.9 3.3×10-6 2.59 0.9 6.9×10-6 4.44 0.75 4.1×10-6 2.75 0.85 6.7×10-6 4.29 0.79 6.7 3.9×10-6 3.23 0.93 3.9×10-5 5.96 0.18 1.6×10-5 5.71 0.67 1.5×10-5 6.34 0.73 7.5 4.1×10-7 2.9 0.87 6.9×10-7 3.51 0.75 3.6×10-6 8.17 / 4.2×10-7 2.81 0.87 平均 8.7×10-5 4.31 0.86 5.3×10-4 10.78 0.34 4.4×10-5 6.05 0.51 5.7×10-5 5.79 0.79 注:“/” 表示相关系数为负,该模型在此处与实际入渗过程不符. 表 2 矿块2入渗过程的拟合结果
Table 2 Fitting results of infiltration process in ore block two
台阶深度/m 考斯加科夫模型 通用经验公式模型 霍顿模型 菲利普模型 均方根误差 相对差值/% 相关系数 均方根误差 相对差值/% 相关系数 均方根误差 相对差值/% 相关系数 均方根误差 相对差值/% 相关系数 0.8 4.4×10-5 2.65 0.95 3.2×10-4 4.80 0.60 8.8×10-5 3.72 0.89 1.9×10-4 5.21 0.79 1.6 3.5×10-6 1.50 0.99 6.2×10-5 3.30 0.80 4.2×10-5 3.93 0.86 4.5×10-5 4.47 0.87 2.4 7.3×10-5 9.50 0.70 3.6×10-4 11.90 / 4.5×10-5 6.52 0.79 1.4×10-4 13.46 0.44 3.2 2.0×10-5 2.97 0.83 2.1×10-5 3.03 0.80 1.6×10-5 2.50 0.86 1.8×10-5 2.70 0.85 4 1.2×10-5 4.08 0.84 6.9×10-5 6.63 / 1.3×10-5 4.06 0.81 2.7×10-5 6.35 0.64 4.8 1.6×10-5 4.44 0.91 6.2×10-5 6.89 0.62 1.6×10-4 15.50 0.01 4.0×10-5 7.51 0.78 5.6 9.4×10-7 2.67 0.93 2.0×10-6 3.62 0.83 2.1×10-6 3.25 0.82 1.8×10-6 3.44 0.86 6.4 2.1×10-5 6.44 0.75 9.5×10-5 7.70 / 1.7×10-5 5.81 0.77 4.1×10-5 9.01 0.51 7.2 4.0×10-6 2.90 0.84 7.9×10-6 3.62 0.65 3.0×10-5 7.18 / 6.3×10-6 3.03 0.75 平均 2.1×10-5 4.08 0.87 1.1×10-4 5.77 0.49 4.9×10-5 5.93 0.65 5.4×10-4 6.05 0.74 注:“/” 表示相关系数为负,该模型在此处与实际入渗过程不符. 表 3 矿块1和矿块2各台阶的入渗速率
Table 3 Infiltration rate of every testing point in ore block one and two
矿块 1 矿块 2 台阶深度/m 初始入渗率 /(cm•min-1) 稳定入渗率 /(cm•min-1) 台阶深度/m 初始入渗率 /(cm•min-1) 稳定入渗率 /(cm•min-1) 1.1 0.341 0.108 0.8 0.387 0.239 1.9 0.047 0.038 1.6 0.223 0.132 2.7 0.105 0.07 2.4 0.157 0.08 3.5 0.118 0.089 3.2 0.169 0.129 4.3 0.103 0.068 4 0.127 0.083 5.1 0.093 0.052 4.8 0.13 0.067 5.9 0.086 0.06 5.6 0.054 0.038 6.7 0.093 0.055 6.4 0.113 0.064 7.5 0.032 0.024 7.2 0.099 0.073 平均 0.113 0.063 平均 0.162 0.101 表 4 矿块1各台阶矿层的物理参数
Table 4 Physical parameters of every testing point in ore block one
台阶深度/m 粒径分布/% 平均容重/(g•cm-3) 平均孔隙比 平均含水量 /% 细砾 砂粒 粉粒 黏粒 1.1 57.43 0.01 20.55 22.01 1.415 1.165 20.82 1.9 30.19 11.92 41.58 16.31 1.285 1.107 19.61 2.7 38.05 12.22 37.57 12.16 1.409 0.915 18.72 3.5 42.19 13.54 36.14 8.13 1.381 0.999 18.94 4.3 37.15 11.67 39.72 11.46 1.308 1.038 20.45 5.1 34.39 8.95 32.12 24.54 1.376 0.851 21.07 5.9 32.31 16.9 36.65 14.14 1.255 0.926 22.02 6.7 40.23 6.07 34.72 18.98 1.311 1.005 20.41 7.5 31.95 8.25 42.15 17.65 1.27 1.069 21.55 平均 38.21 9.95 35.69 16.15 1.334 1.008 20.4 表 5 矿块2各台阶矿层的物理参数
Table 5 Physical parameters of every testing point in ore block two
台阶深度/m 粒径分布/% 平均容重/(g•cm-3) 平均孔隙比 平均含水量 /% 细砾 砂粒 粉粒 黏粒 0.8 52.23 11.01 25.62 11.14 1.423 0.912 16.11 1.6 53.07 9.19 26.42 11.32 1.448 0.885 13.75 2.4 48.92 12.75 27.65 10.68 1.48 0.808 12.7 3.2 52.99 8.59 29.85 8.57 1.333 1.027 17.79 4 39.93 14.16 33.05 12.86 1.351 0.903 16.74 4.8 45.17 7.29 32.55 14.99 1.397 0.957 17.08 5.6 38.41 8.89 34.98 17.72 1.422 0.922 19.14 6.4 40.5 10.72 35.43 13.35 1.327 1.062 20.44 7.2 46.17 11.59 31.04 11.2 1.391 0.907 18.74 平均 46.38 10.47 30.73 12.42 1.397 0.931 16.94 表 6 矿块1和矿块2 考斯加科夫模型参数及渗透系数
Table 6 Kostiakov model parameters and hydraulic conductivity of ore block one and two
矿块 1 矿块 2 台阶深度/m a1 b1 k/(cm•min-1) 台阶深度/m a1 b1 k/(cm•min-1) 1.1 0.824 0.409 0.162 0.8 0.463 0.133 0.358 1.9 0.049 0.049 0.058 1.6 0.280 0.152 0.198 2.7 0.143 0.143 0.106 2.4 0.176 0.158 0.120 3.5 0.137 0.087 0.133 3.2 0.207 0.095 0.194 4.3 0.131 0.131 0.102 4.0 0.141 0.107 0.124 5.1 0.113 0.156 0.078 4.8 0.175 0.192 0.101 5.9 0.101 0.103 0.091 5.6 0.064 0.107 0.057 6.7 0.109 0.137 0.082 6.4 0.122 0.129 0.096 7.5 0.036 0.085 0.036 7.2 0.108 0.078 0.110 -
[1] 金解放,陶伟,邱灿,等. 离子型稀土一维垂直入渗规律及最大粒径的影响试验研究[J]. 有色金属科学与工程, 2015, 6(6): 125-131. http://ysjskx.paperopen.com/oa/DArticle.aspx?type=view&id=201506023 [2] 朱冰冰,张平仓,丁文峰,等. 长江中上游地区土壤入渗规律研究 [J].水土保持通报, 2008, 28(4): 43-47. http://www.cnki.com.cn/Article/CJFDTOTAL-STTB200804009.htm [3] KOSTIAKOV A N. On the dynamics of the coefficient of water percolation in soils and on the necessity of studying it from a dynamic point of view for purposes of amelioration[J]. Soil Science,1932, 97(1):17-21. http://www.oalib.com/references/7959606
[4] 夏江宝,许景伟,李传荣,等. 黄河三角洲退化刺槐林地的土壤水分生态特征[J]. 水土保持通报, 2010, 30(6): 75-80. http://www.cnki.com.cn/Article/CJFDTOTAL-STTB201006017.htm [5] HORTON R E. An approach to ward a physical interpretation of infiltration-capacity[J]. Soil Science Society of America Journal, 1940, 5(3): 399-417. http://www.sciepub.com/reference/149074
[6] PHILIP J R. The theory of infiltration about sorptivity and algebraic infiltration equations[J]. Soil Science, 1957, 84(4): 257-264.
[7] GREEN W H.AMPT G A. Studies on soil physics partl-the flow of air and water through soils[J]. Journal of Agriculture Science, 1911, 76(4): 1-24.
[8] 蒋定生,黄国俊. 黄土高原土壤入渗速率的研究[J]. 土壤学报, 1986, 23(4): 299-305. http://www.cnki.com.cn/Article/CJFDTOTAL-TRXB198604001.htm [9] 李广文,冯起,张福平,等. 祁连山八宝河流域典型草地土壤入渗特征[J]. 干旱地区农业研究, 2014, 32(1): 60-65. http://www.cnki.com.cn/Article/CJFDTOTAL-GHDQ201401011.htm [10] 刘目兴,聂艳,于婧. 不同初始含水率下粘质土壤的入渗过程[J]. 生态学报, 2012, 32(3): 871-878. http://www.cnki.com.cn/Article/CJFDTOTAL-STXB201203024.htm [11] 刘汗,雷廷武,赵军. 土壤初始含水率和降雨强度对黏黄土入渗性能的影响[J]. 中国水土保持科学, 2009, 7(2): 1-6. http://www.cnki.com.cn/Article/CJFDTOTAL-STBC200902000.htm [12] SL345—2007, 水利水电注水试验规程[S]. [13] 赵西宁,吴发启. 土壤水分入渗的研究进展和评述[J]. 西北林学院学报, 2004, 19(1): 42-45. http://www.cnki.com.cn/Article/CJFDTOTAL-XBLX200401011.htm [14] 赵勇钢,赵世伟,曹丽花,等. 半干旱典型草原区退耕地土壤结构特征及其对入渗的影响[J]. 农业工程学报, 2008, 24(6): 14-20. http://www.cnki.com.cn/Article/CJFDTOTAL-NYGU200806005.htm [15] CERDA A. Effects of rock fragment cover on soil infiltration, interrill runoff and erosion[J]. European Journal of Soil Science 2001, 52(1): 59-68. doi: 10.1046/j.1365-2389.2001.00354.x
[16] 李卓,吴普特,冯浩,等. 不同粘粒含量土壤水分入渗能力模拟试验研究[J]. 干旱地区农业研究, 2009, 27(3): 71-77. http://www.cnki.com.cn/Article/CJFDTOTAL-GHDQ200903016.htm [17] 罗嗣海,黄群群,王观石,等. 离子型稀土浸矿过程中渗透性变化规律的试验研究[J]. 有色金属科学与工程, 2014, 5(2): 95-99. http://ysjskx.paperopen.com/oa/DArticle.aspx?type=view&id=201402017