Soil particle size distribution characteristics of ionic rare earth：A case study in rare earth mine of Jiangwozi in Gan County
-
摘要: 颗粒粒度分布对土壤水分运动和溶质迁移，以及水土流失有重要影响. 利用水洗分筛和激光粒度分析技术研究了江西赣县姜窝子稀土矿土壤颗粒粒度并探讨了土壤颗粒粒度在垂直方向的分布规律. 研究表明，姜窝子稀土矿山土壤颗粒粒度分布受风化程度、颗粒垂向迁移因素的影响，呈如下分布:①离子型稀土矿土壤颗粒按直径分为粗、细、粉、黏粒结构，且土壤粗颗粒以石英为主，稀土矿土壤颗粒粒度呈驼峰式分布，即粒径在粗颗粒（ > 2 mm） 和黏粒（ < 0.075 mm）含量较高，粒径在细、粉粒（0.075~2 mm）含量相对较低;②受风化程度和颗粒垂向迁移的影响，稀土矿土壤细、黏粒组分的含量随深度增大呈先升后降的曲线形态；③黏粒（ < 0.075 mm）激光粒度分析表明，颗粒分布随粒径的大小变化与粒径大于0.075 mm的颗粒在深度上的分布相似.Abstract: The particle size distribution could have impact on water movement in soil and solute migration as well as water loss and soil erosion. In order to study the distribution characteristics of soil particle size in the vertical direction, water sieve analysis and laser particle size analysis of soil particle size of rare earth ore in Jiangwozi in Gan county, Jiangxi were conducted. The result shows that the distribution of particle size is affected by weathering degree and vertical migration of particles. The distribution is as follows: (1)According to its diameter, soil particles can be divided into coarse, fine, powder, clay structure and the coarse particles of the soil are mainly quartz, in addition, the distribution of soil particles in the rare earth mine can be characterized as hump type. The contents of particle size in the coarse particles ( > 2mm) and clay ( < 0.075mm) are high, while the content of the particle size in fine powder (0.075-2mm) is relatively low.(2) Because of the influence of weathering degree and vertical migration of particles, the content of fine and clay granules in the rare earth mineral increases first and then decreases in tracing pattern.(3) Laser particle size analysis of clay ( < 0.075) shows that particle size distribution changes with the size of the particle size which is similar to the distribution in depth of particles whose size is greater than 0.075 millimeter.
-
离子型稀土富含世界上短缺的中、重稀土元素，具有很高的经济与战略价值，被我国政府列为保护型开采的矿种之一[1-2].含稀土离子的母岩经风化、溶解，产生稀土阳离子，阳离子经黏土矿物吸附富集而形成离子吸附型稀土矿[3-6].离子型稀土矿一般赋存于花岗岩等火成岩的风化壳中，矿体呈层状分布于被风化、剥蚀、冲刷而形成的准丘陵地带中[7-10].该矿种埋藏较浅，所在山体由浅往深部一般是表土层、全风化层、半风化层、微风化层和基岩，稀土离子主要存在全风化层和半风化层中[11].针对离子型稀土矿，目前主要采用地表原地浅井浸出[12].浸矿施工时，由于向山体风化矿层注入大量浸出药剂，使风化矿层含液量增加，同时浸出药剂与稀土矿的置换过程会减小风化矿层的强度参数，使风化矿层稳定性变差，进而发生滑坡[13-14].基于此，本文通过FLAC3D的流固耦合的实例分析，研究某稀土矿山持续注液对山体边坡的孔压场、位移场和安全系数的影响情况.研究结果为离子型稀土矿原地浸矿采场的施工提供理论依据.
1 某稀土矿工程概况
该稀土矿山体浅层为厚约1~10 m的表土层、全风化矿层，为离子型稀土赋存层，下部基岩为完整性好的花岗岩，是山体含水层下隔水体.山体边坡长约50 m，宽约30 m，高约30 m，山体坡度25°~40°；矿区植被多处中幼龄；矿区地表水主要是山间沟谷溪流，平水期水流量约0.5~100 L/s，受季节影响较大.矿山平均每天向山体内部注入约100 m3浸矿药剂，有研究结果表明，当液固比为0.8时能有效控制母液量和母液浓度[15].图 1为该矿山原地浸矿边坡剖面图.
2 数值计算模型及边界条件
2.1 几何模型
建立数值模型时网格的划分影响到计算的效率. FLAC3D的流固耦合分析模型包含流、固体2种模型，FLAC3D在流固耦合问题上需要加倍的计算步，网格划分的疏密程度极大地影响计算的速度，因此在建模时将数值模型含矿层网格划分较密，基岩强度较高，形变小，网格划分相对较疏.本次数值模型建立过程如下：将原地浸矿边坡的左下侧角点定义为原点O（0，0），边坡最右侧边界的X=50 m；数值模型的底部边界的Z=-10 m，顶部边界纵坐标Z=30 m，边坡总面积约为800 m2.为简化计算，将三维边坡稳定分析简化为二维的平面问题，平面应变计算模型如图 2所示，从上往下分别是全风化矿层、半风化层、基岩，模型单元数为10 600，节点数13 131.
![]()
图 2 原地浸矿边坡数值模型平面模型Figure 2. Local dissolving leaching mining numerical model of slope plane model2.2 边界条件及参数
依据数值模型及原地浸矿施工工艺的特点，确定流固体分析的力学边界条件时固定X=-10 m、X=50 m面的水平位移，固定Y轴方向的前后边界的水平位移，固定Z=-10 m面的垂直位移.做流体分析时设置基岩为不渗透材料，设置坡面与坡顶面为透水自由面.
对岩土工程数值分析影响较大的是材料的物理力学参数，其取值的准确性与合理性直接关系到计算结果的准确性.参考类似矿山及岩土质材料强度参数，本次数值分析材料物理力学参数见表 1所示.
表 1 模型材料物理力学参数Table 1. Material physical and mechanical parameters
3 注液边坡稳定性数值分析
根据基岩在山脚出露与否，集液时采用明沟和暗沟来集液.明沟集液的集液明沟设置在山脚沟谷部位，紧挨矿体一侧，浸出液能沿基岩斜面向山体坡面渗出，进入集液明沟.暗沟网集液的暗集液沟是由地表的坡脚处向山体内打高1.6 m、宽0.6 m的倾斜梯形暗沟，浸出液沿暗沟流出地表.本研究的稀土矿实例中基岩在山体坡脚处出露，浸出液能沿基岩渗出到地表，因此采用明沟集液的原地浅井浸出法采矿.
本次数值分析计算设置为岩土的干密度，由FLAC3D软件自动计算每个单元的饱和重度.流固耦合分析计算过程为：① 计算含矿山体边坡的原始应力；② FLAC3D流固耦合排水分析中，孔压场不与应力场耦合，可以加快运算速率.计算中先单渗流计算，后单力学计算.
3.1 孔隙水分析
浸矿注液使山体坡含液量迅速增加，溶浸液沿基岩斜面向下渗流，并在基岩斜面最先达到饱和，形成有一定厚度的饱和体.注液施工持续时，前期注液量大于渗出的集液量，饱和体厚度增大，液面逐渐接近山体坡面.中期注液量与集液量相当，饱和体厚度将不变，后期注液量减小或停止注液时，饱和体厚度减小.注液量为100 m3/d的不同时刻坡体孔隙水压力云图如图 3所示.因FLAC3D渗流计算中默认未饱和体的孔隙水压力为0，因此认为图 3中蓝色区域均未饱和.由图 3（a）~图 3（d）可以看出，饱和体厚度迅速增大，在坡脚位置增大最快，说明注液5~20 d内的溶浸液注入速度大于浸出液渗出速度.由图 3（d）~图 3（f）可知，饱和体厚度基本不变，说明注液20~30 d内溶浸液注入速度与浸出液渗出速度持平，且注液20 d时饱和矿层与近地表风化层（注液孔区域先受注液渗流浸润作用）都能受到浸矿液浸润. 图 4是孔隙水压力随注液时间变化的折线图.由图 4可以看出，注液5 d、10 d、15 d、20 d、25 d、30 d的孔隙水压力为37.75 kPa、53.11 kPa、69.57 kPa、93.82 kPa、96.51 kPa、97.36 kPa.这些数据可以说明，注液5~20 d的孔隙水压力每5 d最少增大15.36 kPa，20 d后每5 d最大增加2.69 kPa，说明注液20 d后，浸矿边坡每天的注液量比收液量略大.
![]()
图 3 注液期每隔5 d的孔隙水压力云图Figure 3. Pore water pressure graphics of liquid injection every 5 d![]()
图 4 最大孔隙水压力随注液时间变化图Figure 4. Graphics of maximum pore water pressure with liquid injection time variation3.2 位移分析
原地浅井浸出法的施工过程中，浸矿药剂从注液孔注入，使风化矿层岩土的黏结力降低.文献[15]的研究表明，原地浸矿使风化矿层的黏结力在15 d内会减小50 %，因此在模拟分析中需对边坡材料参数进行折减.本次模拟分析中的参数折减方案是：在注液中的第5 d、10 d、15 d，将风化矿层的内聚力折减为注液前的0.85、0.7、0.5倍.
在FLAC3D渗流模式下的原地浅井浸出分析表明，随着注液施工的进行，使原本稳定的边坡发生了较大位移.图 5为注液达20 d的边坡垂直位移云图，图 6为边坡最大垂直位移随注液时间累积折线图.从图 5~图 6可以看出，注液后边坡产生较大位移，边坡注液5 d、10 d、15 d、20 d的最大位移分别是-2.79 mm、-4.66 mm、-10.33 mm、-24.44 mm（文中位移为负，表示位移方向与坐标正轴方向相反，），注液期前期山体内部溶液相对较少，溶浸液与稀土离子的交换反映较缓慢，对岩土边坡的总体强度影响较小，边坡是稳定的.当注液时间达15~20 d时，坡体含液量大，大部分风化矿层都受到溶浸液的浸润作用，离子交换反映剧烈，加上溶液自重的影响，边坡稳定性变差，发生较大位移.当注液时间超过20 d时，计算不收敛，坡脚位置全风化层全部为剪应力屈服状态，表明已经发生滑坡，滑面为全风化矿层沿半风化层滑移.分析表明：连续注液超过20 d时，边坡越来越趋近滑坡临界点，应采取维护措施防止滑坡.
![]()
图 5 注液20 d边坡垂直位移云图Figure 5. Graphics of adding liquid 20 d vertical slope displacement![]()
图 6 边坡最大位移随时间累积折线图Figure 6. Line chart with largest slope displacement accumulated over time3.3 安全系数计算
使用FLAC3D对边坡稳定性安全系数进行求解，得到边坡持续注液20 d内不同时间点的安全系数. 图 7为安全系数变化折线图.由图 7可以看出：未浸矿施工时的安全系数是1.74；施工5 d后边坡安全系数是1.71，与未施工时的安全系数相比变化较小，边坡较稳定；施工10 d、15 d对应的安全系数是1.64、1.56，可以看出浸矿药剂开始发挥作用，减小了坡体岩土强度，对边坡稳定性影响较大；当施工时间达20 d时对应边坡安全系数是1.31，是因为浸矿药剂使边坡的黏结力大幅度减小，而大量药剂存于坡体内增大坡体自重.根据图 7折线曲率变化，以及注液量依旧偏大，当注液达到20~25 d，边坡安全系数将越来越趋近于1.上述分析表明：随着浸矿注液量及浸矿时间的增加，土体含液量大幅增加，黏结力不断减小，使得边坡从稳定状态趋于不稳定状态.
3.4 暗沟集液
根据以上对边坡浸矿的数值分析可知，边坡在注液20 d后的每天浸出母液量比注液量略小，孔隙水压力虽缓慢增大，但位移迅速增大，说明边坡处于不稳定状态.认为20 d后100 m3/d的注液量过大，应适当减小.
为使浸矿施工能持续进行，且不影响浸矿效率，注液量的减小应适当，另可通过在坡脚处开挖暗沟，增加浸出液渗出通道，让坡脚处及附近区域的浸出母液快速流出，以减小坡体的含液量，维持边坡稳定.参考类似暗沟开挖工艺，认为暗沟的布置为高1.6 m，宽0.6 m左右的倾斜梯形暗沟，暗沟间距8 m，每米沟的渗液面积达到4 m2，可以达到很好的收液效果.
为不影响前期注液，暗沟从注液的第20 d开始向外导出浸出液为最佳，则暗沟的开挖时间应该在注液第20 d左右，考虑施工的安全，可根据现场实际情况，将暗沟用密集的倒流管代替.注液第20 d将注液量减小到80 m3/d.图 8~图 9为注液第30 d的孔隙水压力云图和位移云图.图 8中黑色网格为暗沟位置，暗沟延伸至基岩面，计算中设置暗沟为自由透水边界.从图 8可以看出暗沟起到了截流的作用.计算结果表明，注液25 d、30 d的最大孔隙水压力分别为85.33 kPa、85.61 kPa，说明25~30 d内的孔隙水压力变化很小.和开挖暗沟前相比较，第20 d的孔隙水压力明显更小.由图 8可知边坡坡脚处溶液大部分释放出来，边坡最大位移为-26.16 mm. 图 10为30 d内边坡最大位移累积折线图，由图 10可知开挖暗沟后边坡位移变化很小.继续对边坡安全系数求解可得第25 d、30 d的安全系数均为1.27，也说明暗沟对维持边坡稳定起到重要作用.
![]()
图 9 注液30 d边坡垂直位移云图Figure 9. Graphics of vertical slope displacement after 30 d adding liquid![]()
图 10 边坡最大位移随时间累积折线图Figure 10. Line chart with largest slope displacement accumulated over time4 结论
运用FLAC3D软件可以有效模拟计算离子型稀土原地浸矿边坡孔隙水压力、位移、安全系数随注液时间的变化.模拟分析结果表明：
1）注液20 d后，坡体内部饱和液面能覆盖大部分风化矿层.注液持续到5 d、10 d、15 d、20 d时，边坡最大孔隙水压力分别为37.75 kPa、53.11 kPa、69.57 kPa、93.82 kPa，变化较大，20 d、25 d时边坡最大孔压为96.51 kPa、97.36 kPa，变化较小，说明注液0~20 d平均每天注液量远大于集液量，开挖暗沟后5 d内平均每天注液量极为接近集液量.
2）边坡在注液0~20 d内位移不断增大，第20 d时最大位移达-24.44 mm.开挖集液暗沟后，第5 d、10 d（即注液第25 d、30 d）边坡最大位移值为-25.26 mm、-26.16 mm，相比第20 d的位移值变化较小，边坡是稳定的.
3）随着注液施工的进行，边坡安全系数在注液的第0 d、5 d、10 d、15 d、20 d分别为1.74、1.71、1.64、1.56、1.31，开挖暗沟后5 d内边坡的安全系数稳定在1.27，表明边坡能维持稳定，集液暗沟起到重要作用.
-
![]()
图 2 姜窝子稀土矿土壤颗粒粒度比重分布柱状图
Fig 2. Jiang wobi proportion of rare earth soil particle size distribution histogram
![]()
图 3 赣县姜窝子稀土矿土壤粒径小于0.075 mm颗粒级配曲线
Fig 3. Grain -size distribution curve of grain size < 0.075 mm component of REE ore from Jiangwozi.
表 1 姜窝子稀土矿土壤颗粒水洗筛分法分析结果/g
Table 1 Results of particle size analysis of soil samples from Jiang wozi REE ore deposit by washing sieving method /g
样品编号 > 5 mm > 2 mm > 1 mm > 0.5 mm > 0.25 mm > 0.075 mm < 0.075 mm GX-01 55.69 171.02 56.84 63.25 30.90 19.81 165.36 GX--02 36.09 166.56 61.99 77.39 40.71 33.00 142.57 GX-03-2 10.73 153.75 86.61 123.79 58.51 42.04 196.17 GX-03-1 19.10 166.76 78.15 101.79 47.79 36.67 218.07 GX-04-5 43.56 207.57 83.80 98.79 52.82 35.26 171.84 GX-04-4 41.90 242.54 97.79 109.26 59.11 33.74 180.35 GX-04-3 19.30 202.93 92.85 117.18 64.02 46.14 200.23 GX-04-2 33.00 196.14 82.03 90.64 56.36 38.65 181.10 GX-04-1 15.77 175.14 98.16 109.70 51.94 57.70 179.22 GX-05 35.93 208.18 103.71 119.91 61.89 47.18 144.90 
下载: 导出CSV
表 2 粒径特征参数/μm
Table 2 Characteristic parameters of particle size /μm
样品编号 D10 D25 D50 D75 D90 GX-01 0.18 2.06 4.78 11.56 19.08 GX-02 0.18 2.32 7.69 15.6 23.21 GX-03-2 0.21 2.61 8.35 15.19 22.83 GX-03-1 0.17 1.82 5.55 14.76 22.92 GX-04-5 0.17 1.88 5.83 14.42 22.17 GX-04-4 0.18 2.36 8.14 15.07 23.08 GX-04-3 0.2 2.97 9.78 19.82 29.18 GX-04-2 1.83 3.75 9.07 15.26 22.68 GX-04-1 0.2 3.35 10.52 21.3 33.4 GX-05 0.18 2.92 10.06 20.27 32.13 D10：颗粒累积分布为10%的粒径，即小于此粒径颗粒含量占所测颗粒的10%.
下载: 导出CSV
表 3 姜窝子稀土矿土壤颗粒组成与不同类型土壤颗粒组成的对比
Table 3 Ginger snap mine soil particle composition and contrast of different kinds of soil from Jiangwozi
土壤类型 样本度量 粒级/% 参考文献 > 50 μm 10~50 μm 5~10 μm < 5 μm 姜窝子稀土矿 样本最大值 81.15 16.98 5.10 14.99 文中实验数据 样本最小值 69.19 10.47 2.88 7.08 [20] 黄土 样本最大值 59.17 79.76 19.91 31.41 [21] 样本最小值 10.03 31.60 3.78 0.17 [22] 河流沉积物 样本最大值 91.60 45.40 17.70 7.00 [23] 样本最小值 30.80 2.30 2.30 0 [24] 红土 样本最大值 3.64 44.50 29.49 37.58 [19] 样本最小值 0.35 32.46 25.28 28.31 [18]
下载: 导出CSV
-
[1] HUANG G H, ZHANG R D.Evaluation of soil water retention curve with the pore-solid fractal model[J].Geoderma, 2005, 127(2):52-61. http://www.sciencedirect.com/science/article/pii/S001670610400299X
[2] WANG D, FU B J, ZHAO W W, et al.Multifractal characteristics of soil particle size distribution under diferent land-use types on the Loess Plateau, China[J].Catena, 2008, 72(4): 29-36. http://www.irgrid.ac.cn/handle/1471x/666742
[3] 吴爱祥, 尹升华, 李建锋.离子型稀土矿原地溶浸溶浸液渗流规律的影响因素[J].中南大学学报(自然科学版), 2005, 36(3):506-510. http://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD200503031.htm [4] 罗嗣海, 袁磊, 王观石, 胡世丽, 王小玲.浸矿对离子型稀土矿强度影响的试验研究[J].有色金属科学与工程, 2013, 4(3):58-61. http://ysjskx.paperopen.com/oa/DArticle.aspx?type=view&id=201303012 [5] TIAN J.Leaching hydrodynamics of weathered elution-deposited rare earth ore[J].Trans.Nonferrous Met.Soc.China, 2001, 11(3): 434-437. http://www.cqvip.com/QK/85276X/200103/1001449337.html
[6] TIAN J, Chi R, Yin J.Leaching process of rare earthsfrom weathered crust lution-deposited rare earth ore [J].Trans.Nonferrous Met.Soc. China, 2010, 20(5):892-896. doi: 10.1016/S1003-6326(09)60232-6
[7] 金解放, 陶伟, 邱灿, 郭钟群.离子型稀土一维垂直入渗规律及最大粒径的影响试验研究[J].有色金属科学与工程, 2015, (6):125-131. http://ysjskx.paperopen.com/oa/DArticle.aspx?type=view&id=201506023 [8] 汤洵忠, 李茂楠, 杨殿.离子型稀土矿分类之浅见[J].湖南有色金属, 1998, 14(6):1-4. http://www.cnki.com.cn/Article/CJFDTOTAL-HNYJ199806000.htm [9] 张恋, 吴开兴, 陈陵康, 等.离子型稀土矿石颗粒粒度分布及变化规律研究[J].中国稀土学报, 2016, 37(3):67-74. http://www.cnki.com.cn/Article/CJFDTOTAL-XTZZ201603012.htm [10] JOHNSON M J, ANDERSON M A, GREEN S, et al.Physical and hydraulic properties of weathered granitic rock in southern California[J].Soil Science, 1994, 158(5):375-380. doi: 10.1097/00010694-199411000-00009
[11] 尚彦军, 王思敬, 岳中琦, 等.全风化花岗岩孔径分布-颗粒组成-矿物成分变化特征及指标相关性分析[J].岩土力学, 2004, 25(10): 1545-1550. doi: 10.3969/j.issn.1000-7598.2004.10.007 [12] 涂新斌, 王思敬, 岳中琦.风化岩石的破碎分形及其工程地质意义[J].岩石力学与工程学报, 2005, 24(4) : 587-595. http://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200504008.htm [13] 王德, 傅伯杰, 陈利顶, 赵文武, 汪亚峰.不同土地利用类型下土壤粒径分形分析—以黄土丘陵沟壑区为例[J].生态学报, 2007, 27(7):3081-3089. http://www.cnki.com.cn/Article/CJFDTOTAL-STXB200707052.htm [14] 吴景贵.土壤颗粒的功能研究进展[J].吉林农业大学学报, 2008, 30(4):529-537. http://www.cnki.com.cn/Article/CJFDTOTAL-JLNY200804029.htm [15] 张俊华, 李国栋, 南忠仁, 肖洪浪.黑河绿洲区耕作影响下的土壤粒径分布及其与有机碳的关系[J].地理研究, 2012, 31(4):608-618. http://www.cnki.com.cn/Article/CJFDTOTAL-DLYJ201204004.htm [16] 桂东伟, 雷加强, 曾凡江, 穆桂金, 潘燕芳, 李开封.绿洲农田不同深度土壤粒径分布特性及其影响因素—以策勒绿洲为例[J].干旱区研究, 2011, 28(4):622-629. http://www.cnki.com.cn/Article/CJFDTOTAL-GHQJ201104011.htm [17] 吴开兴, 朱平, 张恋, 等.江西龙南乡标稀土矿石干筛法和水洗筛分法颗粒粒度分析[J].中国稀土学报, 2016, 33(5):632-640. http://www.cnki.com.cn/Article/CJFDTOTAL-XTXB201605018.htm [18] 方红, 陈秀玲, 李志忠, 等.福建北部闽江二级阶地第四纪红土粒度特征及其成因研究[J].福建师范大学学报(自然科学版), 2016, 32(5):80-85. http://www.cnki.com.cn/Article/CJFDTOTAL-FJSZ201605014.htm [19] 朱丽东, 谷喜吉, 叶玮, 等.洞庭湖周边地区第四纪红土粒度特征及环境意义[J].地理科学进展, 2014, 33(1):13-22. doi: 10.11820/dlkxjz.2014.01.002 [20] 魏传义, 李长安, 康春国, 常国瑞.哈尔滨黄山黄土粒度特征及其对成因的指示[J].地球科学(中国地质大学学报), 2015, 40(12):1945-1954. http://www.cnki.com.cn/Article/CJFDTOTAL-DQKX201512001.htm [21] 张威, 郭善莉, 李永化, 李云艳.辽东半岛黄土粒度分维特征及其环境意义[J].地理科学进展, 2010, 29(1):79-86. doi: 10.11820/dlkxjz.2010.01.011 [22] 刘庆宏, 贾亮.大连滨海七顶山黄土沉积学特点分析[J].辽宁师专学报(自然科学版), 2010, 12(2):50-51+89. http://www.cnki.com.cn/Article/CJFDTOTAL-LAON201002034.htm [23] 李长安, 张玉芬, 袁胜元, 等."巫山黄土"粒度特征及其对成因的指示[J].地球科学(中国地质大学学报), 2010, 35(5):879-884. http://www.cnki.com.cn/Article/CJFDTOTAL-DQKX201005018.htm [24] 王随继, 闫云霞, 颜明, 等.张家界甘溪砾石沉积物粒度的空间变化及其原因[J].地理科学进展, 2014, 33(1):34-41. doi: 10.11820/dlkxjz.2014.01.004
计量
- 文章访问数: 39
- HTML全文浏览量: 33
- PDF下载量: 9
