SIMULATION MODELING AND SAMPLE SIZE EFFECT ANALYSIS FOR GRAVEL SOIL TRIAXIAL TEST
-
摘要: 借助离散单元理论与室内三轴试验, 分析碎石土物理实验中试样尺寸(直径和高度)变化对应力应变、体应变、粘聚力和内摩擦角等力学性能的尺寸效应影响。研究结果表明, 所提出的块体随机生成法则能较好地实现试样中不同形状碎石块的模拟。不同直径和高度的碎石土试样, 初始2%应变范围内的应力-应变曲线变化规律基本一致, 后部曲线变化较大。直径和高度越小, 围压越低, 应力软化现象越明显; 直径大于250 mm、高度大于350 mm后无应变软化, 残余应力恒定。峰值应力随试样直径增大以25%的增长率呈线性增长, 随高度的增大呈非线性增长, 高度小于200 mm时增长率为11.9%, 大于200 mm后为28.9%。体应变破坏峰值则表现为先增大后减小的趋势。同时粘聚力随直径的增大线性增长, 内摩擦角则减小, 而试件高度的变化对其影响规律则相反。Abstract: Based on the size effect in the gravel soil physics experiment, mechanical properties of the stress-strain, the strain, the cohesive force and inner friction angle of gravel soil was analyzed using discrete element theory and indoor triaxial test. The results showed that the proposed block generated law could realize the simulation of the crushed stone much more. The stress strain curves were almost unanimous (2% strain) at different diameter and height, but the rear part of the curve changes greatly, the smaller the diameter and height, the lower confining pressure stress, the more obvious softening phenomenon. The residual stress was constant after diameter 250 mm, height 350 mm. The peak stress increased with a linear growth rate of 25% with the sample diameter, but nonlinear growth along with the increase of height, when less than 200 mm, the high growth rate was 11.9%, and 28.9% after 200 mm. Damage body strain values were first increasing then decreasing, at the same time cohesive force increases of the linear growth with the diameter and in friction angle was reduced, but rule was just the opposite at different height.
-
Key words:
- gravel soil /
- triaxial simulation /
- size effect
-
大宝山钼多金属矿床是粤北地区最重要的多金属矿床之一。从开采至今已半个多世纪，目前矿区内已探明的资源量已经极度短缺，矿山生产面临严重危机，寻找可接替资源成为当前的首要任务。本文以大宝山矿区的岩性、构造、蚀变及矿化为重点，结合它们与成矿作用之间相互关系，对大宝山矿区外围——白面山成矿预测区开展地质找矿工作。通过对白面山地区的野外地质调查，发现该区断裂构造发育，中基性岩脉广泛分布。沿断裂带及其两侧围岩产生多种类型、不同程度的蚀变及矿化。对区内断裂构造、岩浆活动及与矿化、蚀变特征之间相互关系的研究认为，预测区内成矿条件较好，找矿潜力大。因此，将其圈定为大宝山矿田内首个找矿预测区，开展进一步研究。采集围岩(蚀变灰岩)及岩浆岩(石英闪长玢岩、辉绿岩)样品，进行等离子体质谱(ICP-MS)分析。通过对样品中的成矿元素、微量和稀土元素地球化学特征以及与成矿作用之间关系的研究，结合该区成矿地质特征和相关物化探工作的研究成果，开展综合分析，对白面山地区进行找矿预测。
1. 区域地质背景
白面山地区位于广东省韶关市曲江区大坑口乡，隶属大宝山钼多金属矿床矿田，距大宝山矿区西南约7 km。区域构造位置位于华南褶皱系与湘粤桂褶皱带南缘，粤北古生代凹陷带和吴川—四会深大断裂带内，地处北东向北江断裂带和东西向大东山—贵东断裂带的交汇处，曲仁盆地南缘。
研究区出露地层主要为上泥盆统东岗岭组、天子岭组的灰岩、泥灰岩、生物碎屑灰岩、白云岩；下石炭统帽子峰组、大赛坝组和长来组的粉砂质泥岩、泥灰岩等(见图 1)。
图 1 白面山成矿预测区地质简图1—中泥盆统东岗岭组灰岩；2—上泥盆统天子岭组灰岩；3—上泥盆统-下石炭统帽子峰组泥岩；4—下石炭统大赛坝组+长来组粉砂质泥岩；5—第四系覆盖物；6—石英闪长玢岩脉；7—辉绿岩脉；8—铁帽；9—实测断层；10—推测断层；11—地质界线；12—推测地质界线Figure 1. A sketch geological map of the Baimianshan metallogenic forecasting area区内构造以断裂为主，褶皱不发育。断裂构造主要发育北东向、北西向和近东西向3组，以北西向断裂为主。区内发育3条近平行的北西向断裂，倾向南西，倾角50°~80°，具有多期活动的特征；东西向断裂为北东、北西2组断裂的次级断裂，多终止于与北东向或北西向断裂交汇处。受区域性深大断裂带长期活动的影响，中基性岩脉较发育，主要为石英闪长玢岩脉和辉绿岩脉，分别沿北西向和北东向断裂产出。其中石英闪长玢岩脉条数较多，宽度较大，10~50 m不等，延伸较远；辉绿岩脉发育较少，宽度相对较小，一般小于10 m。
研究区内断裂构造发育，伴随岩浆活动，地表沿断裂带发育大面积的蚀变、矿化现象。矿化、蚀变是围岩遭受岩浆热液活动影响最直观的表现形式。蚀变类型包括白云岩化、硅化以及少量方解石化、重晶石化。其中白云岩化、硅化最为发育。大部分蚀变都具有2种或2种以上蚀变类型共生的特征。矿化主要表现为黄铁矿化、褐铁矿化、方铅矿化、闪锌矿化、辉锑矿化以及少数铜蓝矿化、金银矿化等。金属矿物呈细脉状、浸染状。黄铁矿化最多，与其他金属矿物以伴生或共生的方式产出；方铅矿化、闪锌矿化也相对集中，表现为多个矿化点；其他类型矿化极为分散，偶尔可见。矿化与蚀变受热液活动影响共同产出于断裂带内或断裂带附近，呈带状分布。
2. 地球化学特征
本文样品均采自于断裂带内或断裂带附近具有一定程度蚀变、矿化的灰岩及相对新鲜的石英闪长玢岩、半氧化状态的辉绿岩，包括12块蚀变灰岩样品，2块石英闪长玢岩样品，2块辉绿岩样品(其中1块氧化成为铁帽)。测试分析其成矿元素、微量和稀土元素含量，研究其地球化学元素特征。
样品由核工业北京地质研究院分析测试研究中心测试。将测定的样品粉碎、研磨到0.074 mm(200目)以下粉末备用。称取样品0.05 g放入25 mL专用溶样罐中，先用少量水湿润，轻轻震动使样品均匀，加入1 mL氢氟酸，3 mL硝酸，1 mL高氯酸，盖上专用溶样罐盖，在低温电热板上以200 ℃恒温加热溶解，待样品分解后，打开溶样罐，在低温电热板上加热蒸至近干，滴加2滴高氯酸，再次蒸至近干，后加入1:1硝酸3 mL，盖上专用溶样罐盖焖置一段时间。用1%硝酸提取至50 mL容量瓶中，摇匀后采用美国热电公司ELEMENT XR等离子体质谱(ICP-MS)分析仪利用在线内标(Rh)法对样品进行分析测试，测试结果见表 1—表 4。
表 1 白面山预测区蚀变灰岩微量元素含量(×10-6)Table 1. Trace elements of the meta-alteration limestone in Baimianshan forecasting area序号 岩性 Rb Ba Th U Ta Nb La Ce Sr Nd Zr Hf Sm Y Yb Lu Mo Tl Bi Cd Sb W Zr/Hf Nb/Ta B004 中厚层灰黑色灰岩 6.950 17.300 0.565 1.120 0.164 0.540 2.050 4.030 663 1.800 6.130 0.168 0.372 1.870 0.163 0.024 0.166 0.078 0.047 0.206 0.675 1.200 36.488 3.293 B016 厚层状灰黑色灰岩 27.800 333.0 2.360 1.370 0.205 2.330 6.320 12.70 415 5.850 30.70 0.815 1.240 5.620 0.596 0.097 0.418 0.267 0.221 0.317 5.700 0.590 37.669 11.367 B043 厚层状灰白色灰岩 7.690 3488.0 0.762 0.463 0.045 0.514 4.770 9.000 410 3.890 6.820 0.222 0.804 3.730 0.341 0.057 0.248 0.326 1.240 0.791 3.330 0.441 30.721 11.422 B054-2 厚层状灰白色灰岩 17.500 64.2 1.560 0.830 0.122 1.590 3.220 6.630 180 2.960 23.10 0.657 0.634 3.150 0.366 0.061 0.138 0.137 0.095 11.40 3.030 0.670 35.160 13.033 B057-1 灰白色灰岩 14.200 49.5 2.600 1.230 0.162 2.100 8.540 16.70 392 8.200 25.20 0.737 1.710 8.030 0.682 0.107 0.976 0.131 0.134 0.089 2.520 0.303 34.193 12.963 B058 土黄色泥质灰岩 16.100 90.2 2.100 0.513 0.164 1.840 4.800 10.80 250 5.590 13.90 0.395 1.550 8.250 0.644 0.093 0.499 0.133 0.166 0.180 2.890 0.455 35.190 11.220 B077 中厚层灰白色灰岩 8.220 25.1 0.709 1.230 0.039 0.527 3.690 6.470 336 2.520 8.350 0.243 0.474 2.190 0.210 0.039 0.674 0.154 0.076 0.188 1.240 0.271 34.362 13.513 B081 厚层灰白色灰岩 23.400 57.6 1.920 0.997 0.194 2.220 8.960 15.10 488 6.060 27.10 0.794 1.080 4.850 0.496 0.077 1.010 0.216 0.466 0.368 16.20 0.657 34.131 11.443 B0128 中厚层状灰岩 39.900 232.0 2.560 1.010 0.154 2.030 9.970 16.60 151 8.830 26.00 0.783 2.190 10.20 1.060 0.158 0.224 0.269 0.267 0.175 2.810 0.821 33.206 13.182 B0153 厚层状灰白色灰岩 5.210 27.30 0.828 0.705 0.095 1.060 3.550 6.890 531 3.180 18.60 0.483 0.666 3.470 0.305 0.047 0.322 0.093 0.160 27.20 1.840 0.323 38.509 11.158 B0159 中厚层状微晶灰岩 1.780 51.90 0.128 0.273 0.020 0.097 1.540 2.310 466 0.966 1.160 0.036 0.173 0.924 0.062 0.010 0.118 0.047 0.057 2.930 1.080 0.101 32.222 4.850 B0185 方解石化泥晶灰岩 5.920 35.30 0.569 0.942 0.049 0.448 1.720 2.920 527 1.270 5.670 0.157 0.245 1.360 0.130 0.023 0.227 0.105 0.072 1.460 1.140 0.276 36.115 9.143 注：B0153号样品为预测区内石英闪长玢岩岩脉南侧采集的铅锌矿石 表 2 白面山预测区岩浆岩中微量元素含量(×10-6)Table 2. Trace elements of the meta-magmatic rocks in Baimianshan forecasting area序号 岩性 Rb Ba Th U Ta Nb La Ce Sr Nd Zr Hf Sm Y Yb Lu Mo Tl Bi Cd Sb W Zr/Hf Nb/Ta B020 石英闪长玢岩 35.2 100 3.15 0.802 1.74 26.2 64.8 84.5 24.8 90.3 433 11.8 22.4 103 9.09 1.38 0.681 0.34 0.023 0.656 24.10 1.21 36.695 15.057 B032 石英闪长玢岩 29.3 24.8 1.95 0.423 1.33 19.8 16.3 42.5 89.3 32.5 202 5.88 8.49 48.2 4.3 0.599 0.86 0.257 0.036 0.388 33.70 5.12 34.354 14.887 B0165 辉绿岩 1.95 11.2 0.161 0.639 0.012 0.106 1.61 2.96 516 1.2 1.07 0.033 0.202 0.9 0.072 0.013 0.061 0.031 0.084 0.389 0.708 0.08 32.424 8.833 B0176 辉绿岩 10.8 4608 1.51 1.26 0.062 0.673 10.3 25.6 9.7 12.3 16.4 0.518 3.02 8.11 1.03 0.148 1.21 0.428 0.158 0.604 7.75 0.504 31.660 10.855 注：B0165号辉绿岩样品发育褐铁矿化，B0176号样品氧化成为铁帽 表 3 白面山预测区蚀变灰岩稀土元素含量(×10-6)Table 3. REE elements of the meta-alteration limestone in Baimianshan forecasting area序号 岩性 La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu ΣREE LREE HREE LREE/HREE (La/Yb)N δEu δCe B004 中厚层灰黑色灰岩 2.050 4.030 0.456 1.800 0.372 0.072 0.327 0.054 0.324 0.061 0.165 0.030 0.163 0.024 9.928 8.780 1.148 7.648 9.020 0.617 0.980 B016 厚层状灰黑色灰岩 6.320 12.700 1.430 5.850 1.240 0.215 1.010 0.173 1.000 0.195 0.545 0.102 0.596 0.097 31.473 27.755 3.718 7.465 7.610 0.569 0.990 B043 厚层状灰白色灰岩 4.770 9.000 0.968 3.890 0.804 0.284 0.730 0.118 0.727 0.139 0.364 0.065 0.341 0.057 22.257 19.716 2.541 7.759 10.030 1.112 0.970 B054-2 厚层状灰白色灰岩 3.220 6.630 0.760 2.960 0.634 0.124 0.556 0.098 0.607 0.116 0.326 0.063 0.366 0.061 16.521 14.328 2.193 6.533 6.310 0.624 1.000 B057-1 灰白色灰岩 8.540 16.700 2.050 8.200 1.710 0.347 1.440 0.253 1.460 0.286 0.724 0.118 0.682 0.107 42.617 37.547 5.070 7.406 8.980 0.658 0.950 B058 土黄色泥质灰岩 4.800 10.800 1.260 5.590 1.550 0.428 1.290 0.243 1.500 0.272 0.695 0.121 0.644 0.093 29.284 24.428 4.856 5.030 5.350 0.900 1.050 B077 中厚层灰白色灰岩 3.690 6.470 0.669 2.520 0.474 0.090 0.442 0.072 0.411 0.079 0.209 0.037 0.210 0.039 15.412 13.913 1.499 9.282 12.600 0.591 0.940 B081 厚层灰白色灰岩 8.960 15.100 1.650 6.060 1.080 0.234 1.030 0.168 0.942 0.184 0.494 0.088 0.496 0.077 36.563 33.084 3.479 9.510 12.960 0.669 0.890 B0128 中厚层状灰岩 9.970 16.600 2.190 8.830 2.190 0.518 1.860 0.347 2.090 0.395 1.030 0.186 1.060 0.158 47.424 40.298 7.126 5.655 6.750 0.765 0.830 B0153 厚层状灰白色灰岩 3.550 6.890 0.786 3.180 0.666 0.312 0.617 0.106 0.638 0.124 0.317 0.056 0.305 0.047 17.594 15.384 2.210 6.961 8.350 1.463 0.970 B0159 中厚层状微晶灰岩 1.540 2.310 0.244 0.966 0.173 0.041 0.173 0.026 0.147 0.028 0.072 0.012 0.062 0.010 5.804 5.274 0.530 9.951 17.820 0.717 0.830 B0185 方解石化泥晶灰岩 1.720 2.920 0.336 1.270 0.245 0.049 0.251 0.039 0.232 0.046 0.126 0.022 0.130 0.023 7.409 6.540 0.869 7.526 9.490 0.599 0.880 注：B0153号样品为预测区内石英闪长玢岩岩脉南侧采集的铅锌矿石；测试单位为核工业北京地质研究院分析测试研究中心 表 4 白面山预测区岩浆岩中稀土元素含量(×10-6)Table 4. REE elements of the meta-magmatic rocks in Baimianshan forecasting area序号 岩性 La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu ΣREE LREE HREE LREE/HREE (La/Yb)N δEu δCe B020 石英闪长玢岩 64.80 84.50 19.00 90.30 22.40 7.54 19.30 3.81 24.00 4.32 10.70 1.740 9.090 1.380 362.880 288.540 74.340 3.881 5.11 1.0819 0.58 B032 石英闪长玢岩 16.30 42.50 6.56 32.50 8.49 3.01 7.36 1.57 11.00 2.04 4.92 0.841 4.300 0.599 141.990 109.360 32.630 3.352 2.72 1.1368 1.01 B0165 辉绿岩 1.610 2.96 0.33 1.20 0.20 0.05 0.19 0.03 0.16 0.03 0.08 0.013 0.072 0.013 6.927 6.345 0.582 10.902 16.04 0.7251 0.95 B0176 辉绿岩 10.30 25.60 3.16 12.30 3.03 0.65 2.08 0.38 2.13 0.38 1.02 0.174 1.030 0.148 62.385 55.041 7.344 7.495 7.17 0.7501 1.09 注：B0165号辉绿岩样品发育褐铁矿化，B0176号样品氧化成为铁帽 2.1 成矿元素与微量元素特征
从蚀变灰岩及岩浆岩(石英闪长玢岩、辉绿岩)中主要成矿及微量元素测试结果(见表 1、表 2)可以看出，蚀变灰岩(样品B043、B054-2、B0153) 中成矿元素铅、锌的含量明显偏高，比其他样品高出1~2个数量级，说明蚀变灰岩发生不同程度的铅、锌矿化。所有样品中微量元素Ba含量普遍偏高，而B043样品中Ba的含量比其他样品高出1~2个数量级。
对比蚀变灰岩、辉绿岩和石英闪长玢岩的测试结果(见表 1、表 2)发现，岩浆岩中铅、锌等成矿元素平均含量明显高于蚀变灰岩。因此，可以推测岩浆岩(石英闪长玢岩)可能为蚀变灰岩的矿化提供物质来源。
蚀变灰岩微量元素配分曲线总体一致(见图 2)，呈现平缓的右倾型、M型多峰谷模式。蚀变灰岩中除Rb、La、Sr、Sm、U等元素相对富集外，其他微量元素均有亏损。Nb、Ta、Zr、Hf、Th等高场强元素，由于离子半径小、电荷高，且难溶于水，化学性质稳定，为非活动性元素，因此可指示原始物质的组成[1~2]。这些高场强元素均表现为负异常，且Zr/Hf比值相近，Nb/Ta比值除了蚀变灰岩中的B004、B0159(黑色灰岩、微晶灰岩)相对较低以外，其余白色蚀变灰岩中Nb/Ta比值大致相等。因此，白面山预测区的蚀变灰岩具有相同的物质来源，同时个别样品(B0159中厚层状微晶灰岩)Nb更具负异常。大离子亲石元素Ba既有负异常，也存在少数正异常，主要是受热液活动或重晶石化蚀变的影响。
W、Mo、Bi属高温成矿元素，通常与岩浆热液活动有关，Cd、Tl、Sb是典型的亲硫重金属元素，常作为Pb、Zn的伴生元素赋存于Pb、Zn的硫化物矿床中[3]；Sb又是热液活动的指示元素[4]，因此Sb的富集说明灰岩的蚀变、矿化是受岩浆热液作用所引起的。
石英闪长玢岩和辉绿岩中成矿及微量元素变化特征规律性较差(见表 2，图 3)。两种不同岩性表现出明显的差异。总体而言，Rb、La、Sm、U等元素相对富集，其他微量元素相对亏损。2块石英闪长玢岩样品及褐铁矿化的辉绿岩样品中Sr元素表现为明显亏损，与蚀变灰岩中Sr元素的特征恰好相反。其余微量元素变化特征大致相同，只是含量差异较大。
所有石英闪长玢岩、辉绿岩岩石样品中Zr/Hf及Nb/Ta比值非常相近，说明石英闪长玢岩和辉绿岩2种岩脉可能来源于同一岩浆房[5]。
2.2 稀土元素特征
蚀变灰岩(白云岩化、方解石化、弱硅化灰岩)和岩浆岩(石英闪长玢岩、辉绿岩)稀土元素含量采用文献[6]球粒陨石数据标准化。
白面山地区蚀变灰岩稀土总量(ΣREE)较低(见表 3)，ΣREE平均值为23.524×10-6，轻重稀土含量比值LREE/HREE平均为7.561，(La/Yb)N平均值9.606。δEu变化范围0.569~1.463，平均值0.774；δCe变化范围0.83~1.05，平均值0.94，二者均存在个别正异常，大多数表现为负异常。稀土元素配分曲线呈右倾型(见图 4)，表现为轻稀土富集，重稀土相对亏损。同时反映出大多数蚀变灰岩样品的稀土元素配分曲线变化一致，仅表现出不同样品中稀土元素总量存在一定差异。除此之外，大多数蚀变灰岩δEu值介于0.569~0.765之间，接近于平均值，δEu表现为显著的负异常；而2块厚层状白色灰岩样品(B043和B0153) 中δEu值分别为1.112和1.463，明显表现为正异常。此2块灰岩样品同时发育强烈的方铅矿矿化，且采样点靠近石英闪长玢岩脉。由此推断它们之所以出现δEu正异常，有可能是受到了石英闪长玢岩脉中高温流体的影响[7]。石英闪长玢岩脉岩浆热液中的Pb可能为该区成矿元素Pb的重要成矿物质来源。2块中厚层白色灰岩和土黄色泥质灰岩样品(B054-2和B058) 中δCe值分别为1.00和1.05，几乎无异常；而其他蚀变灰岩中δCe值介于0.83~0.99之间，表现出微弱的负异常。
白面山地区的石英闪长玢岩和辉绿岩2种脉岩的稀土元素特征见表 4、图 5。二者稀土元素配分曲线均为平缓的右倾式，表现为轻稀土富集，重稀土相对略微亏损。2个石英闪长玢岩样品的ΣREE值较高，分别为362.88×10-6和141.99×10-6；LREE/HREE比值相近；(La/Yb)N值相差较大；δEu值分别为1.0819和1.1368，均表现出弱的正异常；δCe值分别为0.58和1.01，前者表现为明显的负异常，后者则几乎无异常。2个辉绿岩样品的ΣREE值分别为6.927×10-6和62.385×10-6，LREE/HREE及(La/Yb)N值差异明显；δEu值分别为0.7251和0.7501，均表现为负异常，δCe值分别为0.95和1.09，前者表现为弱负异常，后者表现为弱正异常。
对比两种脉岩的稀土元素特征发现，石英闪长玢岩的ΣREE值远大于辉绿岩，并且2件辉绿岩样品的ΣREE值差别也很大。这是由于2块辉绿岩样品均遭受不同程度风化作用的影响，前者表现为褐铁矿化，后者则被氧化成为铁帽。因此推断其原岩中的一部分稀土元素可能受风化作用影响而损失。但是辉绿岩样品B0176的ΣREE值是B0165样品ΣREE值的近10倍，而δEu和δCe值的差异很小，几乎相等。表 2中B0176样品成矿元素Cu、Ba、Tl含量都远远高于B0165样品，Ba、Tl是反应热液作用特征的重要指示元素[4]，而Cu的矿化也反映出高温成矿作用的特点。因此推断该样品可能受到成矿期热液作用的影响使稀土元素发生二次叠加。
徐晓春等[8]认为岩浆岩稀土元素地球化学特征发生规律性的变化，反映其可能是同源岩浆演化的结果。对比岩浆岩样品稀土元素地球化学特征，推断预测区内2种脉岩可能为同一岩浆源经岩浆分异作用所产生的2种不同产物[9]。
通过对白面山预测区内蚀变灰岩、岩浆岩脉中成矿及微量元素、稀土元素地球化学特征分析，认为该区围岩(灰岩)的蚀变及矿化是岩浆热液沿断裂带活动的结果，且为成矿作用提供物质来源。
3. 找矿标志及找矿预测
3.1 找矿标志
3.1.1 构造标志
白面山地区断裂构造发育，且以北西向断裂为主，在石英闪长玢岩北部接触带附近的该组断裂内铅、锌矿化岩石标本(见图 6a)，化学分析结果：铅1.4%，锌0.84%；在石英闪长玢岩南部接触带附近的北西向断裂中的铅、锌矿化岩石标本(见图 6b)，化学分析结果：铅2.14%，锌12.0%。北东向断裂为压性，延伸方向与区域构造线方向一致；北东东及北西西—北西向次级断裂均为其伴生构造。北西向断裂是良好的储矿(岩)构造，断裂不仅为成矿热液运移提供通道，也为成矿物质的就位提供空间，尤其是断裂的交汇复合部位是矿体最有利的赋存场所。
图 6 B0153铅锌矿石金属矿物组合特征及野外露头Figure 6. B0153 lead-zinc ore metal mineral assemblages and field outcrop3.1.2 岩性标志
研究区岩性有沉积岩和岩浆岩。沉积岩主要为上泥盆统至下石炭统的碳酸盐岩，是良好的容矿围岩；岩浆岩主要为北北西向断裂带中充填的中酸性次英安斑岩、石英闪长玢岩及中性辉绿岩脉，是成矿热液来源。碳酸盐岩带状蚀变及脉状矿化，显示岩浆热液作用的结果；稀土元素地球化学特征证明岩浆岩来源于同一岩浆源区，并利用物探方法推测在其深部可能有中酸性隐伏岩体存在，可以为成矿提供充足的物质来源。
3.1.3 化探标志
基于广东省地矿局1960、1961、1988、1989年在白面山地区多次开展大比例尺化探填图的成果，结合野外地质调查及室内工作，确认白面山地区发育多处铅、锌异常点，主要集中在飞凤形一带的南部，该处具有适合的成矿条件，地表矿化、蚀变较强且集中。指示该区具有理想的成矿地球化学环境，深部可能存在隐伏的铅、锌矿体。
3.1.4 物探标志
广东省物化探院2008年1:10000高精度磁测资料显示，飞凤形一带深部可能存在隐伏岩体。中国地质调查局武汉地质调查中心2012年在该地区榕树下村一带开展激发极化法测量，发现深部具有多个高阻体和高极化体。结合蚀变、矿化分布及岩脉产出特征，推测深部可能存在北西、东西向岩体以及北西和近南北向的矿化体(见图 7)。
图 7 白面山成矿预测区综合物探推测剖面图(a, b为同一测量点两个相互垂直剖面)(据广东物化探院[2008]、中国地质调查局武汉地质调查中心[2012])Figure 7. A comprehensive geophysical suppositional profile map of Baimianshan metallogenic forecasting area3.1.5 其他标志
主要有矿化和蚀变标志。矿化是成矿作用在地表最直接的反应，通过地表矿化可初步推断深部矿床特征。矿化标志主要有黄铁矿化、方铅矿化、闪锌矿化以及少量的辉锑矿、铜蓝矿等矿化。除黄铁矿化普遍发育，可见立方体状的单颗粒黄铁矿呈浸染状产出外，其他类型矿化仅在个别样品中出现，并且大多数样品中很难见到单颗粒金属矿物，仅可通过化学分析所得到的测试数据来推断其具有某种类型的矿化。
3.2 找矿预测
在该区进行野外地质调查过程中发现大面积、多种类的矿化，主要为多个铅、锌异常带，且异常带均发育于石英闪长玢岩脉附近；蚀变灰岩(矿石)呈密集浸染状铅锌矿化(见图 6c—6f)。对石英闪长玢岩岩脉附近的矿化蚀变灰岩进行采样分析测试，测试结果显示出其中Zn含量最高，为12%，Pb含量大于2%，Ag含量高达5.98 g/t。结合该区的物化探研究成果及硅化、白云岩化和弱的方解石化、重晶石化等蚀变，显示所有的铅、锌异常带均与物化探所揭示出的异常区域非常吻合。因此，推断该预测区地表及地表以下深部可能存在矿化体(矿体)，是最有可能形成中—低温矽卡岩型铅、锌、银多金属矿床的地区。
4. 讨论
白面山地区位于大宝山矿田范围内，同属一个基底构造隆起带。从研究成果来看，白面山预测区成矿地质条件与大宝山矿区相似。
大宝山矿区主要矿产有铜、铅、锌、钼、钨等，铜、铅、锌为矽卡岩型矿床，产于泥盆系中，受北西向、北东向及近东西向断裂构造控制，矿体呈透镜状；成矿岩体主要为受北西向逆冲推覆断裂带所控制的次英安斑岩；成矿时代为燕山早期(距今162~168 Ma)[10~11]。白面山远景区沿北西、北东向断裂带发育有较多的铅、锌矿化，可能与北北西向的次英安斑岩岩脉有关。白面山地区岩浆岩稀土元素地球化学特征研究显示，石英闪长玢岩和辉绿岩脉为同一岩浆源区经岩浆分异作用的产物。但是白面山次英安斑岩是否与大宝山次英安斑岩来自同一岩浆源区还需要进一步验证。如果它们来自同一岩浆源，则大宝山矿区内发育一套高—中—低温的矿化-蚀变组合，白面山表现为中—低温矿化-蚀变组合，推断大宝山为成矿区的中心，白面山为成矿区的边缘，这样白面山就很难形成具工业价值铅锌银矿体。但如果它们来自不同岩浆源区，则说明白面山地区深部可能存在单独的成矿岩体为该区成矿提供物质来源。
5. 结论
预测区内蚀变灰岩和岩浆岩中微量元素、稀土元素地球化学特征显示，该区断裂带中的灰岩发育强烈的蚀变及矿化，是岩浆热液活动的结果。
岩浆岩微量及稀土元素地球化学特征证明，白面山预测区的石英闪长玢岩和辉绿岩来自于同一岩浆源区，是岩浆分异作用的产物；石英闪长玢岩可能为该区的成矿岩体，为矿化及成矿物质的来源。
结合地表围岩矿化及蚀变特征分析，预测区中西部地表以下500~800 m之间可能存在一个中—低温矽卡岩型铅锌银多金属矿体。
-
表 1 模拟级配与试样级配对比
Table 1. Comparison of simulated grading and sample grading
三轴试验试样尺寸和级配 数值模拟模型尺寸和级配 试样尺寸/mm 粒径/mm 含量/% 模型尺寸/mm 粒径/mm 含量/% 101×200 20~10 21.72 101×200 20~10 21.72 10~5 20.19 10~5 20.19 5~2 18.09 5.0~2.5 18.09 ＜2 40 2.5 40 
下载: 导出CSV
表 2 参数标定值
Table 2. The calibration parameters
参数含义 程序代码 参数取值 杨氏模量 md_Ec 4.5e6 平行粘结模量 pb_Ec 1e7 密度 md_dens 2600 土颗粒法向刚度 kn 7.2e4. 块体法向刚度 kn 4e5 法向刚度与切向刚度比 kn/ks 0.8 接触粘结法向力平均值 cb_sn_mean 1.35e5 接触粘结切向力平均值 cb_ss_mean 1.35e5 平行粘结法向力平均值 pb_sn_mean 7e4 平行粘结切向力平均值 pb_ss_mean 7e4
下载: 导出CSV
表 3 模拟实验与三轴试验误差分析
Table 3. Error analysis in simulation experiment and triaxial test
实验类型 峰值偏应力/kPa 粘聚力c/kPa 内摩擦角φ/(°) 围压100 kPa 围压200 kPa 围压400 kPa 三轴试验 501.80 677.90 1119.20 79.89 30.67 模拟实验 466.91 708.41 1016.43 91.33 28.38 误差/% 6.95 4.50 9.18 14.30 7.45
下载: 导出CSV
表 4 数值模拟工况组合
Table 4. Numerical simulation case combinations
编号 试样尺寸/mm 直径 高度 试样1 50 200 试样2 70 200 试样3 101 200 试样4 150 200 试样5 200 200 试样6 250 200 试样7 101 100 试样8 101 150 试样9 101 200 试样10 101 250 试样11 101 300 试样12 101 350
下载: 导出CSV
表 5 不同试样直径下峰值体应变
Table 5. Peak body strain under different diameter
试样直径/
mm峰值体应变 100 kPa 200 kPa 400 kPa 50 0.0264 0.0481 0.0174 70 0.0326 0.0477 0.0452 101 0.0568 0.0587 0.0372 150 0.0286 0.02840 0.03640 200 0.0113 0.00983 0.00282 250 -0.0021 -0.01360 -0.03350
下载: 导出CSV
表 6 不同试样高度下峰值体应变
Table 6. Peak body strain under different height
试样高度/
mm峰值体应变 100 kPa 200 kPa 400 kPa 100 0.0279 0.0121 0.0074 150 0.0379 0.0309 0.0151 200 0.0568 0.0587 0.0372 250 0.0557 0.0448 0.0912 300 0.0563 0.0462 0.0396 350 0.0133 0.0137 0.0888
下载: 导出CSV
-
[1] 谈云志, 郑爱, 吴翩, 等.红黏土承载比的土团尺寸效应研究[J].岩土力学, 2013, 34(5):1242~1246 http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201305002.htmTAN Zhi-yun, ZHEN Ai, WU Pian, et al. Effect of aggregate soil size on Californi a bearing ratio values of laterite soil[J]. Rock and soil mechanics, 2013, 34(5): 1242~1246. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201305002.htm [2] 凌华, 殷宗泽, 朱俊高, 等.堆石料强度的缩尺效应试验研究[J].河海大学学报, 2011, 39(5):540~544 http://www.cnki.com.cn/Article/CJFDTOTAL-HHDX201105016.htmLING Hua, YIN Zong-ze, ZHU Jun-gao, et al. Experimental studyof scaleeffect onstrengthof rockfill materials[J]. Journal of HohaiUniversity, 2011, 39(5): 540~544. http://www.cnki.com.cn/Article/CJFDTOTAL-HHDX201105016.htm [3] 徐文杰, 胡瑞林, 岳中琦, 等.土石混合体细观结构及力学特性数值模拟研究[J].岩石力学与工程学报, 2007, 26(2):300~311 http://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200702010.htmXU Wen-jie, HU Rui-lin, YUE Zhong-qi, et al. Mesostructural character and numerical simulation of mechanical properties of soil-rock mixtures[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(2): 300~311. http://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200702010.htm [4] The Influence of Rock Block Contents on Geo-Mechanical Behavior of Soil-Rock Mixture[C] //2012颗粒材料计算力学研究进展. 2012. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-DLLG201209001032.htm [5] 狄圣杰, 汪明元, 邓南沙, 等.土石混合体REV尺度数值方法研究[J].防灾减灾工程学报, 2012, 32(5):523~526 http://www.cnki.com.cn/Article/CJFDTOTAL-DZXK201205002.htmDI Sheng-jie, WANG Yuan-ming, DENG Nan-sha, et al. Numerical Simulation Study on REV's Measure of Soil-rock Mixture[J]. Disaster Prevention and Mitigation Engineering, 2012, 32(5): 523~526. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXK201205002.htm [6] Minneapolis M. PFC2D particle flow code in 2 dimensions theory and background [R]. Itasca Consulting Group Incorporation, 2004. doi: 10.1134/S1062739116041046 [7] 刘宝琛, 张家生, 杜奇中, 等.岩石抗压强度的尺寸效应[J].岩石力学与工程学报, 1998, 17(5):611~614 http://www.cnki.com.cn/Article/CJFDTOTAL-DZLX201702017.htmLIU Bao-chen, ZHANG Jia-sheng, DU Qi-zhong, et al. A study of size effect for compression strength of rock[J]. Chinese Journal of Rock Mechanics and Engineering, 1998, 17(5): 611~614. http://www.cnki.com.cn/Article/CJFDTOTAL-DZLX201702017.htm [8] 潘一山, 魏建明.岩石材料应变软化尺寸效应的试验和理论研究[J].岩石力学与工程学报, 2002, 21(2):215~218 http://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200202013.htmPAN Yi-shan, WEI Jian-ming. Experimental and theoretical study on size effect on strain softening of rock materials[J]. Chinese Journal of Rock Mechanics and Engineering, 2002, 21(2): 215~218. http://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200202013.htm [9] 李翀, 何昌荣, 王琛, 等.粗粒料大型三轴试验的尺寸效应研究[J].岩土力学, 2008, 29(增刊):563~566 http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2008S1113.htmLI Xu, HE Chang-rong, WANG Chen, et al. Study of scale effect of lage-scale triaxial test of coarse-grained meterials[J]. Rock and soil mechanics, 2008, 29(Supplement.): 563~566. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2008S1113.htm [10] 阳云华, 赵晏, 关沛强, 等.膨胀土抗剪强度的尺寸效应研究[J].人民长江, 2007, 38(9):18~22 http://www.cnki.com.cn/Article/CJFDTOTAL-RIVE200709010.htmYANG Yun-hua, ZHAO Yan, GUAN Pei-qiang, et al. A study of size effect for shear strength of the expansive soil[J]. Ren Min Chang Jiang Journal, 2007, 38(9): 18~22. http://www.cnki.com.cn/Article/CJFDTOTAL-RIVE200709010.htm [11] 杨圣奇, 苏承东, 徐卫亚, 等.岩石材料尺寸效应的试验和理论研究[J].工程力学, 2005, 22(4):112~118 http://www.cnki.com.cn/Article/CJFDTOTAL-GCLX200504021.htmYANG Sheng-qi, SU Cheng-dong, XU Wei-ya, et al. Experimental and theoretical study of size effect on of rock material[J]. Engineering mechanics, 2005, 22(4): 112~118. http://www.cnki.com.cn/Article/CJFDTOTAL-GCLX200504021.htm [12] 吕兆兴, 冯增朝, 赵阳升, 等.岩石的非均质性对其材料强度尺寸效应的影响[J].煤炭学报, 2007, 32(9):917~920 http://www.cnki.com.cn/Article/CJFDTOTAL-MTXB200709004.htmLV Zhao-xing, FENG Zeng-chao, ZHAO Yang-shen, et al. Influence of rock in homogeneity on strength2size effect of rock materials[J]. Journal of coal, 2007, 32(9): 917~920. http://www.cnki.com.cn/Article/CJFDTOTAL-MTXB200709004.htm [13] 肖昭然, 胡霞光, 刘玉, 等.沥青混合料细观结构离散元分析[J].公路, 2007, (4):146~148 http://www.cnki.com.cn/Article/CJFDTOTAL-GLGL200704039.htmXIAO Zhao-ran, HU Xia-guang, LIU Yu, et al. Discrete element method analysis for micro structure of asphalt mixtures [J]. High Way, 2007, (4): 146~148. http://www.cnki.com.cn/Article/CJFDTOTAL-GLGL200704039.htm [14] 李耀旭. 颗粒流方法在土石混合体力学特性研究中的应用[D]. 长江科学院, 2009 http://cdmd.cnki.com.cn/Article/CDMD-10611-1016731846.htmLI Yao-xu. The application of particle flow code in studing of mechanical characteristics of soil-rock mixture [D]. Yangtze River Scientific Research Institute, 2009. http://cdmd.cnki.com.cn/Article/CDMD-10611-1016731846.htm [15] 王新. 土石混合体力学特性影响因素及破坏机制研究[D]. 长江科学院, 2010 http://cpfd.cnki.com.cn/Article/CPFDTOTAL-DZDQ201501006007.htmWANG Xin. Research on influence factors of mechanics charateristics and failure mechanism of soil-rock mixture [D]. Yangtze River Scientific Research Institute, 2010. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-DZDQ201501006007.htm -
下载:
下载:
计量
- 文章访问数: 236
- HTML全文浏览量: 151
- PDF下载量: 10
- 被引次数: 0
