Distinct element modelling of fracture plan control in continuum and jointed rock mass in presplitting method of surface mining

Controlled blasting techniques are used to control overbreak and to aid in the stability of the remaining rock formation. Presplitting is one of the most common methods which is used in many open pit mining and surface blast design. The purpose of presplitting is to form a fracture plane across which the radial cracks from the production blast cannot travel. The purpose of this study is to investigate of effect of presplitting on the generation of a smooth wall in continuum and jointed rock mass. The 2D distinct element code was used to simulate the presplitting in a rock slope. The blast load history as a function of time was applied to the inner wall of each blasthole. Important parameters that were considered in the analysis were stress tensor and fracturing pattern. The blast loading magnitude and blasthole spacing and jointing pattern were found to be very significant in the final results.

Mass Balance of Major Elements in Relation to Weathering in Soils Developed on Igneous Rocks in a Semiarid Region,Northwestern Iran

This study was conducted to evaluate the weathering intensity of the major soils developed on igneous rocks in semiarid region of northwestern Iran.Eight parent materials were selected including monzodiorite,alkali granite,granodiorite,syenite,pyroxene diorite,hornblende andesite,pyroxene andesite,and dacite.Representative soil profiles were described and soil samples were collected and analyzed for selected chemical and physical properties and total concentrations of major elements and Zr,V,Ti and Y.Bulk densities as well as Ti,Zr and V concentrations were used to estimate the strain factors and mass balance equations were used to quantify the net result of pedogenic weathering,i.e.elemental loss and gain.The results of clay content and pedogenic iron variability as well as index of compositional variability(ICV),chemical index of alteration(CIA) and,A-CN-K and MFW ternary plots showed that the soils developed on volcanic rocks(hornblende andesite> pyroxene andesite> dacite) were more weathered than those on the plutonic parent rocks(alkali granite,granodiorite,monzodiorite,syenite,pyroxene diorite).The results of mass balance calculations based on the strain factors revealed that the Ca and Na depleted during weathering progress mostly from plagioclase grains.In the semiarid regions Ca is precipitated as pedogenic calcite in the soil horizons.K and Mg depletion is less than Ca and Na especially in the profiles on the hornblende andesite with the highest clay and LOI content.The results of this study clearly suggest that the behavior of K and Mg during the weathering cannot only be explained by the disintegration of the primary minerals,since they are fixed on the secondary clay minerals.Iron did not change in the soils compared to the parent material and was precipitated as the pedogenic iron and conserved in the soil horizons.Overall,the results on the weathering indicators and major elements mass balance enrichment/depletion in the study area confirmed that the soil profiles developed on volcanic rocks are more weathered than those on the plutonic igneous rocks.

Characteristics of Element Migration in the Process of Wall-Rock Alteration in the Shibangou Gold Deposit,Western Henan

Abstract The Shibangou gold deposit in western Henan is associated with irregular quartz veinlets occurring in altered shear zones dissecting a dioritic intrusion. The altered shear zones are characterized by silicification, pyritization, sericilization, chloritization and K-feldspar alteration. Zoning of altered rocks adjacent to the Au-bearing quartz veins is obviously exhibited. Fine-grained sulphides and quartz veinlets of different ages and small-scale fissures are widely distributed in the central part of the altered zones. Major mineralization types in this gold deposit are Au-bearing quartz veinlets and altered rocks in the shear zones. Samples were collected from drilling cores according to the alteration zoning and mineralization type and all samples were analyzed for major and trace elements. Mass balance, volume change (fv=97.3–71.9%) and major element variation sequences are studied in terms of major elements. The changes of mobile components (SiO2, K2O, Fe2O3) and CaO) and immobile component (Al2O3) in the wall-rock alteration are discussed. The gold mineralization is associated with the enrichment of As, Ag, Hg and Pb and depletion of Cu and Zn. The study of compositional variation of altered rocks proves to be a very efficient method for defining the extent of wall-rock alteration, fluid activity and mineralization and enrichment.

Reinforcing a Dangerous Rock Mass Using the Flexible Network Method

Because the main failure type of a dangerous rock mass is collapse, the treatment of such a mass should focus on controlling collapse failure. When treating dangerous rock masses, disturbing the mass （e. g. by blasting） needs to be avoided, as this new damage could cause collapse. So the self-bearing capacity of the mountain mass must be used to treat the dangerous rock mass. This article is based on a practical example of the control of a dangerous rock mass at Banyan Mountain, Huangshi, Hubei Province. On the basis of an analysis of damage mechanism and the stability of the dangerous rock mass, a flexible network reinforcement method was designed to prevent the collapse of the rock mass. The deformations of section Ⅱ w of the dangerous rock mass before and after the flexible network reinforcement were calculated using the two-dimensional finite element method. The results show that the maximum deformation reduced by 55 % after the application of the flexible network reinforcement, from 45.99 to 20.75 ram, which demonstrates that the flexible network method is effective, and can provide some scientific basis for the treatment of dangerous rock masses.

Effects of fracture distribution and length scale on the equivalent continuum elastic compliance of fractured rock masses

Fracture systems have strong influence on the overall mechanical behavior of fractured rock masses dueto their relatively lower stiffness and shear strength than those of the rock matrix. Understanding theeffects of fracture geometrical distribution, such as length, spacing, persistence and orientation, isimportant for quantifying the mechanical behavior of fractured rock masses. The relation betweenfracture geometry and the mechanical characteristics of the fractured rock mass is complicated due tothe fact that the fracture geometry and mechanical behaviors of fractured rock mass are stronglydependent on the length scale. In this paper, a comprehensive study was conducted to determine theeffects of fracture distribution on the equivalent continuum elastic compliance of fractured rock massesover a wide range of fracture lengths. To account for the stochastic nature of fracture distributions, threedifferent simulation techniques involving Oda＇s elastic compliance tensor, Monte Carlo simulation （MCS）,and suitable probability density functions （PDFs） were employed to represent the elastic compliance offractured rock masses. To yield geologically realistic results, parameters for defining fracture distributionswere obtained from different geological fields. The influence of the key fracture parameters andtheir relations to the overall elastic behavior of the fractured rock mass were studied and discussed. Adetailed study was also carried out to investigate the validity of the use of a representative elementvolume （REV） in the equivalent continuum representation of fractured rock masses. A criterion was alsoproposed to determine the appropriate REV given the fracture distribution of the rock mass.

Anisotropy of strength and deformability of fractured rocks

Anisotropy of the strength and deformation behaviors of fractured rock masses is a crucial issue for design and stability assessments of rock engineering structures, due mainly to the non-uniform and non- regular geometries of the fracture systems. However, no adequate efforts have been made to study this issue due to the current practical impossibility of laboratory tests with samples of large volumes con- taining many fractures, and the difficulty for controlling reliable initial and boundary conditions for large-scale in situ tests. Therefore, a reliable numerical predicting approach for evaluating anisotropy of fractured rock masses is needed. The objective of this study is to systematically investigate anisotropy of strength and deformability of fractured rocks, which has not been conducted in the past, using a nu- merical modeling method. A series of realistic two-dimensional （2D） discrete fracture network （DFN） models were established based on site investigation data, which were then loaded in different directions, using the code UDEC of discrete element method （DEM）, with changing confining pressures. Numerical results show that strength envelopes and elastic deformability parameters of tested numerical models are significantly anisotropic, and vary with changing axial loading and confining pressures. The results indicate that for design and safety assessments of rock engineering projects, the directional variations of strength and deformability of the fractured rock mass concerned must be treated properly with respect to the directions of in situ stresses. Traditional practice for simply positioning axial orientation of tunnels in association with principal stress directions only may not be adequate for safety requirements. Outstanding issues of the present study and su^zestions for future study are also oresented.

GLOBAL COMPOSITE ELEMENT ITERATION FOR ANALYSIS OF SEEPAGE FREE SURFACE

As one of the most difficult topics in rock mass hydromechanics, seepage free surface plays an important part in slope stability researches. Based on analysis of numerical methods to solve seepage free surface, global composite element iteration (GCEI) is presented in this paper. FEM program is made by using GCEI. Calculation shows that not only the program is simple using GCEI, but also the tolerance is higher after 5 iterations.

块裂隧道围岩等效岩体变形参数研究

采用Voronoi离散单元方法建立块裂岩体隧道的随机DEM数值模型,通过大量随机数值计算结果验证采用岩体等效变形参数预测隧道收敛变形的可行性。在规则裂隙岩体等效变形参数计算方法的基础上,针对数值模拟揭示的深埋隧道变形特征提出改进的修正计算公式及相应的特征参数。通过与既有解答以及大量随机数值计算结果的对比和统计分析验证修正解的合理性,修正解较为合理改善既有规则裂隙岩体等效变形参数计算方法在反映随机块裂岩体在隧道卸荷条件下变形估算的较大误差。采用修正计算方法得到的岩体等效变形参数可以很好地预测随机块裂岩体中隧道开挖后收敛变形的数学期望值。提出的修正计算方法对当前深部岩体变形参数研究和隧道洞壁收敛位移的估算具有一定的参考价值。

隧道穿越岩堆体围岩剪切特性及分级

为解决隧道穿越岩堆体围岩分级不准确的问题,通过相似材料直剪试验分析岩堆体的直剪特性,将含水率和填充土含量两个因素带入熵权可拓物元法,并结合地质指标进行分级研究。结果表明岩堆体的抗剪强度随着含水率、填充土含量的增大,先增大后减小,最终趋于平稳,过高或过低的含水率和填充土都会“弱化”岩堆体的抗剪强度。根据岩堆体的特点,选取岩石单轴饱和抗压强度、节理裂隙密度、岩体完整性系数、含水率和填充土含量作为分级指标,结合熵权可拓物元法,将所选指标和熵权可拓物元法结合并在工程实例进行验证,分级结果比原方法准确性更好,证明研究提出的围岩分级方法是合理可行的。研究结果能够为岩堆体隧道的围岩分级提供指导,为施工设计和支护方式提供依据。

基于优化组合赋权的可拓学磷矿山岩体质量评价

为解决岩体质量评价相邻等级之间指标参数模糊不确定性问题,引入可拓学理论,提出基于优化组合赋权的可拓学岩体质量分级模型,对昆阳二矿磷矿山地下开采的岩石进行质量评价。首先,根据矿山地质特征选取岩石单轴饱和抗压强度(Rc)、岩石质量指标(RQD)、节理间距(Jd)、结构面条件(Jf)、地下水状态(W)和地应力影响系数(Z)6个影响指标,将每个指标划分为5个等级;然后通过改进的层次分析法和熵权法分别确定指标的主观和客观权重,引入矩估计法对主客观权重进行优化组合;最后应用修正的RMR法和Q系统法对岩体质量进行评价,将分级结果与基于优化组合赋权的可拓学岩体质量分级结果进行比较。研究表明:5个待评价岩体中,除上矿体N3的模型分级结果与修正RMR法存在差异之外,其余岩体的分级结果与修正RMR法相同,使用Q系统法的评价结果整体偏低。此外,该模型可得出岩体质量等级偏向于相邻等级的程度,分级结果符合昆阳二矿磷矿山生产勘探地质报告岩体质量范围,比修正RMR法和Q系统法更加准确,验证了本文方法的可靠性,为地下矿山开采提供了技术支撑。

基于关键块体的井巷围岩稳定性分析及支护研究

为了进一步研究关键块体对巷道围岩稳定性的影响及通过关键块体分析进行支护措施的设计。以蒙库铁矿地下开采某一巷道为研究对象,通过工程地质调查与现场测试工作,对围岩节理裂隙统计与分段,并进行岩块强度试验和岩体的波速测试,获得了岩体完整性参数,依据此进行了详细的岩体工程分级,给出不同岩体的分级力学参数,再通过Unwedge进行关键块体分析,对不稳定岩块进行划定,确定了锚喷支护的参数及重点支护位置,然后采用有限元软件模拟开挖过程并分析了开挖过程中围岩及支护的应力应变变化情况,最后进行现场采动应力的监测与分析。研究结果表明:经过关键块体分析,更准确地辨别出可能导致巷道失稳的不稳定岩块,针对关键块体分析结果提出稳定性支护措施,并应用于现场,进行采动应力监测与分析,取得了良好的效果。研究结果在一定程度上为围岩稳定性分析及支护设计提供了可靠依据。

节理裂隙岩体超大断面隧道开挖工法研究

为研究不同开挖工法对围岩中节理发育的超大断面隧道稳定性影响,以在建大罗山公路隧道为工程背景,运用离散元数值分析方法,研究了双侧壁导坑法、CRD法、CD法和三台阶法四种开挖工法下节理面剪切位移、隧道开挖轮廓位移、地表沉降、围岩位移、支护结构内力的变化规律.研究结果表明:双侧壁导坑法在控制节理裂隙岩体各位移变化及隧道支护结构内力方面效果最好,CRD法和CD法次之,不推荐将三台阶法作为节理发育地层中大断面公路隧道的开挖工法.

露天矿高陡边坡岩体参数获取与稳定性评价方法研究

针对露天矿高陡边坡岩体参数难获取和稳定性易变化的问题,以攀枝花露天矿高陡边坡为研究对象,采用非接触式无人机摄影测量技术获取研究区域的三维模型,统计结构面的分布情况,结合室内岩石力学试验数据,采用《工程岩体分级标准》(GB/T 50218—2014)进行岩体质量分级,并给出边坡岩体强度力学参数建议值。在此基础上,通过有限单元法分析了矿区在开采条件下边坡的稳定性动态演化特征,并对开挖过程边坡的稳定性做出了评价。结果表明:随着矿区的开采,边坡在开挖面附近由于卸荷作用发生位移回弹,最大位移增量为33.2 mm;塑性区主要集中在断层和矿岩界面附近,以剪切屈服为主;边坡稳定性随着矿区的开采逐渐降低,直至达到临界稳定状态。采用锚索对边坡进行加固,并采用塑性区面积占比对加固前后的稳定性进行评价,评价系数α由加固前的0.4降低至加固后的0.05,证明了加固措施的有效性。研究成果可为露天矿高陡边坡岩体参数取值和稳定性评价提供参考。

某高速公路高陡危岩体失稳对拟建桥梁的影响分析

以西南山区某高速公路高陡崖边坡对坡脚拟建桥梁影响分析为例,采用工程地质调查分析结合有限元强度折减法分析可知,该陡崖边坡整体处于欠稳定状态,变形模式为错落式滑坡变形。采取有限元-离散元耦合分析技术动态模拟研究边坡失稳对桥梁安全的影响。由模拟分析结果可知,陡崖边坡产生错落变形后,错落岩体最终与桥梁墩台底部发生撞击,桥墩从根部断裂,桥梁整体倒塌。基于分析结果,提出预应力锚索加固陡崖边坡、坡脚设置拦截挡墙、边坡失稳范围内避免设置桥墩3种方案应对高陡边坡对拟建桥梁的破坏性影响。

COMPOSITE ELEMENT MODEL FOR ROCK MASS SEEPAGE FLOW

In this paper, the algorithm developed on the concept of the Composite Element Method （CEM） is proposed for the seepage flow in discontinuous rock masses, The basic idea is that if an element contains one or more discontinuity segments, it will be denoted as a composite element having a set of hydraulic potential values defined at its nodes. By the variation principle, the governing equation can be established to solve the unknown hydraulic potential values at the composite element nodes. The main advantage of the proposed model is that the generation of calculation mesh will not be restricted strongly by the number, the position and the orientation of the discontinuities, A numerical example is presented at the end of the paper, to show the validity of the proposed model.

海南琼中岩-土体系重金属迁移特征及表土环境质量评价

地球关键带和土壤重金属污染是当今研究的前沿与热点,岩-土体系是地球关键带的重要组成部分。目前海南岛针对不同成壤母质岩-土体系的详细研究较少,并且重金属亟需基于新标准进行评价。为了查清岩-土体系中重金属元素迁移特征及污染状况,本文以海南岛琼中3个土壤垂向剖面及7115件表层土壤样品为研究对象,采用电感耦合等离子体质谱/发射光谱法(ICP-MS/OES)和原子荧光光谱法(AFS)等分析方法,测定岩石和土壤中Cu、Pb、Zn、Cr、Ni、Cd、As和Hg共8种重金属元素含量,利用土壤环境质量农用地土壤污染风险管控标准、内梅罗指数法和潜在生态危害指数法定量评价了表层土壤样品重金属污染状况及生态风险。结果表明,①海南岛琼中二长花岗岩岩-土体系中,不同重金属元素表现出明显的分异特征,Cu、Pb、Zn、Cr、Ni、Cd元素以亏损为主,As元素以富集为主,Hg元素则表现出不同的富集或亏损特征,As、Hg元素尤其富集在表层。②琼中绝大多数表层土壤中重金属元素含量低于土壤环境质量农用地土壤污染风险筛选值,土壤污染风险低。③以琼中和海南岛背景值为评价标准时,应关注Cr、Ni和As污染以及Hg、As和Cd的生态危害。以土壤环境质量为评价标准时,研究区表层土壤的环境质量总体清洁,生态风险较弱。

基于LCEM-GFEM方法爆炸载荷作用下含缺陷岩体损伤机理研究

在爆破作业过程中,岩体缺陷对裂纹扩展具有显著影响。本研究基于全局有限元方法的局部黏结单元(LCEM-GFEM)方法构建数值模拟模型,研究爆破作用下含缺陷岩体的损伤演化和裂纹扩展模式。通过引入能量传递系数,定量分析了缺陷形态对应力波传播和衰减的影响。结果表明,缺陷形态对岩体的损伤特性和裂纹扩展具有显著影响。随着平行缺陷夹角的增加,爆破诱发的裂缝与衍生裂缝的合并路径从缺陷末端转移到了缺陷中间,能量传递系数增加,分形维度减小。随着平行缺陷之间的水平距离增加,平行缺陷之间的裂缝数目和能量传递系数相应减少。随着垂直缺陷之间的垂直距离增加,穿过缺陷C的主水平裂缝长度增加,能量传递系数和分形维度也增加。研究结果可为受到爆破荷载作用下含缺陷岩体的损伤特征提供指导。

超大跨度泄水隧洞出口洞段及边坡的风控分析与设计

基于双江口水电站洞式溢洪道出口洞段及边坡的复杂工程地质条件、施工条件、环水保制约因素,通过施工勘察、危岩体风险分析与评价,借助边坡稳定分析、有限元计算等手段,开展治理方案技术经济比较。计算分析成果及实施情况表明,洞室晚出洞,边坡“少开挖、强支护”,危岩体大面防护治理结合关键部位锁口、拦腰、固脚措施,出口洞段与边坡危岩体“分而治之”的系统治理方案是合理可行的,满足工程进度要求的情况下降低了施工难度及风险,节约了工程投资成本,可供类似工程借鉴。

基于物元可拓分析的岩体质量分级方法研究与应用

影响岩体质量的因素众多,岩体质量分级中往往面临因素考虑不全面、主观性强以及不同因素之间难以进行对比分析等问题,对准确性产生影响。物元可拓分析法的最大优势是可以对不同因素进行综合分析,从定性和定量角度对目标进行确定。针对某矿岩体质量评价工作,基于物元可拓分析法建立岩体质量评价体系,通过确定经典域、节域、待评物元以及权重,完成对岩体的质量评价。结果表明,利用该法对岩体质量的评价结果与工程岩体分级标准相近,证明了方法的合理性和准确性。同时,熵权法确定权重更加客观,物元可拓分析法降低了在分析中考虑不全面和专家主观性的影响。

基于科氏质量流量计的微米级煤尘和岩尘颗粒分类方法

【目的】岩尘的粒径比煤尘的小得多,更容易被吸入到肺中,危害较大。为了解决煤尘和岩尘混合颗粒的分类问题,提出一种微米级煤尘和岩尘混合颗粒的测量方法。【方法】首先采用环形静电传感器进行前端流体流量的测定,根据科氏质量流量计的测量原理,研究测量管两侧的时间差和流体的质量流量的关系;其次采用ANSYS有限元软件进行双向流固耦合,验证对煤尘和岩尘混合颗粒进行分类的可行性;最后对测量管进行静力学分析以及数值模拟分析,实现科氏质量流量计对微米级煤尘和岩尘混合颗粒的分类判别。【结果】在科氏质量流量计的同一个测量管道入口处,颗粒的流入速度与测量时间差正相关;谐振式U型管的激振频率应为二阶振型对应的固有频率113.11 Hz;当煤尘和岩尘颗粒混合粒径不同时,静力学参数和时间差随着煤尘的体积分数的增大而减小,并且全岩尘和全煤尘的时间差明显与混合颗粒的不同;当煤尘和岩尘混合颗粒粒径相同时,时间差随着煤尘的体积流量增大而增大,关系曲线的陡度也随着增大。当煤尘体积分数<50%时,煤尘和岩尘的混合颗粒的粒径越大时间差越大;当煤尘体积分数≥50%时,煤尘和岩尘的混合颗粒的粒径越大时间差越小。【结论】煤尘颗粒和岩尘颗粒具有鲜明且不同的分类判别特征,微米级煤尘和岩尘混合颗粒分类具有可行性,可实现科氏质量流量计对煤尘和岩尘混合颗粒的实时在线精确分类。