Effect of Rock Mass Structure and Block Size on the Slope Stability-Physical Modeling and Discrete Element Simulation

This paper studies the stability of jointed rock slopes by using our improved three-dimensional discrete element methods (DEM) and physical modeling. Results show that the DEM can simulate all failure modes of rock slopes with different joint configurations. The stress in each rock block is not homogeneous and blocks rotate in failure development. Failure modes depend on the configuration of joints. Toppling failure is observed for the slope with straight joints and sliding failure is observed for the slope with staged joints. The DEM results are also compared with those of limit equilibrium method (LEM). Without considering the joints in rock masses, the LEM predicts much higher factor of safety than physical modeling and DEM. The failure mode and factor of safety predicted by the DEM are in good agreement with laboratory tests for any jointed rock slope.

节理裂隙岩体爆破试验研究

首先利用超声波探测法对节理裂隙岩体进行了声波测试,对不同的节理裂隙岩体进行了质量评估.然后通过对节理裂隙岩体进行爆破漏斗试验,分析了节理裂隙对岩体中爆炸应力波传播、岩石破碎及爆破漏斗形成的影响,得出了一些有益结论,对节理裂隙岩体理论研究和指导工程实践具有一定的参考价值.

岩体中的节理裂隙对爆破影响的试验研究

利用超声波探测法对节理裂隙岩体进行了声波测试,对不同的节理裂隙岩体进行了质量评估,通过对节理裂隙岩体进行爆破漏斗试验,分析了节理裂隙对岩体中爆炸应力波传播、岩石破碎及爆破漏斗形成的影响.结果表明,节理裂隙等结构弱面的存在严重阻碍了爆炸能量均匀分布,因此在爆破设计和施工中,应尽量使炮孔与结构弱面垂直布置,如果遇到含有夹层情况宜加大装药量;炮孔与结构弱面平行时,应尽量减小炮眼孔距和最小抵抗线,尽量采用单孔螺旋掏槽.

环渤海花岗岩体地下油库建设的工程地质力学研究

As a kind of strategic resource，petroleum play an very important role in current social stability, economic development and state safety. Since 1993 China has turned from a net oil exporter into a net oil importer, the figure of imported oil increased from then on. In 2004 China's total energy consumption exceeded Japan’s, and ranked in the second place, just inferior to America. Today China is the world’s third-largest importing nation, accounting for 6% of world imports and 8% of world consumption. Comparing with other strategic petroleum reserve schemes, underground oil storage possess many advantages, such as security, economy, less pollution, save land, suited for strategic reserve and so on, so it is the most ideal form for strategic petroleum reserve. In the background of China Strategic Petroleum Reserve Program started just now, this paper choose Circum-Bo sea region as a study area, and do some system study on the underground oil storage caverns constructed in inter-large granite rock masses in Circum-Bo sea region. On the foundation of a great amount of information come from both home and abroad, firstly this paper analysed the principle, economy, cavern shape, profile dimension, and gain some cognizances and logos, as follows: ①Hard rock mass such as granite is the major rock, in which underground oil storage are constructed; ②Unlined underground oil storage caverns had been wide spread used as a sort of oil storage form abroad, there already exist a suit of skilled experience and technologies to prevent oil product from leaking; ③Compared with surface tanks, underground oil storage cavern possess predominance in economy clearly. In general, it will be more economical when the storage capacity exceed 50000m3. The quality of rock mass is the most important factor for underground storage cost, however such as hydrogeology, storage capacity, the number of storage galleries, the length, storage product, mechanical equipments, geographic location also influent the cost. In designed depth of the underground storage, the rock mass of Jinzhou mainly belong to class Ⅱ, but parts with dykes, clayization alteration, and dense joints are Ⅲ, Ⅳ; ④Now, there are few underground oil storages span more than 25m in both abroad and home. The examples of some ancient underground works and modern underground excavation with wide span surely give us many precious elicitations to construct more great unlined storage caverns, when the rock mass quality is good, cavern shape and construction method also are proper, it is quite possible to construct underground oil storage cavern with span more than 30m . The main axis orientation of Jinzhou underground oil storage cavern is NW direction, the cavern's elevation locate between -53msl and -76msl. The storage's total volume is about 3×106m3, composed of 8 parallel galleries with 950m length, the pillars between them are 45m, and every two of galleries form one unit, which can deposit 75×104m3 for each unit. The product will be stored are Saudi light and Saudi medium crude oil, the main cavern's section is 411.5m2, with 23m height and 19m width. According to the principle and technique of engineering geomechanics, this study supply a sort of system scientific thinking and method for sitting location of underground oil storage in granite region: ① On the foundation of the earth crust stability sub-zone appraise of Circum-Bo sea region, farther research concerning granite distribution, genesis, geological period and fault structure are conducted in stable areas, generally, this paper select Liaoxi, east shore of Liaotung peninsula and Jiaotung peninsula as target areas for underground oil storage regions, where Mesozoic granite is magnitude; ②After roundly comparison in facts of geologic structure, engineering geology, hydrogeology, topography, transportation and so on of three granite distributed areas, at last, selecting Jingzhou granite zone in Liaoxi out as an ideal construction area; ③ Detailed investigation is conducted in the southeast of Baimashi in Jingzhou development district, the final field. Ultrasonic Borehole Television, as a major way to collect original information of borehole rock mass were used, which is very effective to appraise the quality of deep rock mass; ④ According to the field data of tectonic stress, rock mass quality, the spatial distribution of fracture water, some optimum designs in cross section, axial direction and cavern span have been designed for the underground oil storage cavern layout in Jinzhou. To understand the characteristics of swelling alteration rock in Jinzhou granite mass, collected abundant swelling alteration rock engineering examples in granite, which study them in detail, concluded the swelling alteration rock distribute nearly everywhere in China, intruded medium-basic dykes alteration, along discontinuities and mineral hydrothermal alteration with genesis of granite are three main forms clayization alteration rock in granite rock mass. In Jinzhou field, from macro to micro studied the swelling rock which induced by mid-basic dyke intrusion, with weak swelling. In conclusion, this paper conclude the distribution rule and features of expansion alteration rock in filed, and advise some technical suggestions for excavation at swelling alteration rock part. The main features of this paper: ①In the process of site selection, investigation and design, a suit of technique and method of engineering geomechanics metasynthesis were formed, which is significative to guide the large scale underground oil storage cavern sitting location, investigation and design in granite rock mass; ②The detailed discussion on the engineering geology problems in granite mass, such as weathering crust, faults, dykes and clayization alteration rock, are useful for other projects in aspects of site selection, engineering geology evaluation and stability estimation; ③The summary and integration of the genesis, type, countermeasure relate to swelling alteration rock, also is likely to be used for other underground oil storage caverns constructed in swelling alteration granite. In conclusion, this study is meaningful for guiding the large scale underground oil storage for site selection, investigation and design in granite rock mass.

基于精细结构描述及数值试验的断续节理岩体力学参数智能选取

Evaluating the mechanical properties of rock masses is the base of rock engineering design and construction. It has great influence on the safety and cost of rock project. The recognition is inevitable consequence of new engineering activities in rock, including high-rise building, super bridge, complex underground installations, hydraulic project and etc. During the constructions, lots of engineering accidents happened, which bring great damage to people. According to the investigation, many failures are due to choosing improper mechanical properties. ‘Can’t give the proper properties’ becomes one of big problems for theoretic analysis and numerical simulation. Selecting the properties reasonably and effectively is very significant for the planning, design and construction of rock engineering works. A multiple method based on site investigation, theoretic analysis, model test, numerical test and back analysis by artificial neural network is conducted to determine and optimize the mechanical properties for engineering design. The following outcomes are obtained: (1) Mapping of the rock mass structure Detailed geological investigation is the soul of the fine structure description. Based on statistical window，geological sketch and digital photography，a new method for rock mass fine structure in-situ mapping is developed. It has already been taken into practice and received good comments in Baihetan Hydropower Station. (2) Theoretic analysis of rock mass containing intermittent joints The shear strength mechanisms of joint and rock bridge are analyzed respectively. And the multiple modes of failure on different stress condition are summarized and supplied. Then, through introducing deformation compatibility equation in normal direction, the direct shear strength formulation and compression shear strength formulation for coplanar intermittent joints, as well as compression shear strength formulation for ladderlike intermittent joints are deducted respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. (3) Model test of rock mass containing intermittent joints Model tests are adopted to study the mechanical mechanism of joints to rock masses. The failure modes of rock mass containing intermittent joints are summarized from the model test. Six typical failure modes are found in the test, and brittle failures are the main failure mode. The evolvement processes of shear stress, shear displacement, normal stress and normal displacement are monitored by using rigid servo test machine. And the deformation and failure character during the loading process is analyzed. According to the model test, the failure modes quite depend on the joint distribution, connectivity and stress states. According to the contrastive analysis of complete stress strain curve, different failure developing stages are found in the intact rock, across jointed rock mass and intermittent jointed rock mass. There are four typical stages in the stress strain curve of intact rock, namely shear contraction stage, linear elastic stage, failure stage and residual strength stage. There are three typical stages in the across jointed rock mass, namely linear elastic stage, transition zone and sliding failure stage. Correspondingly, five typical stages are found in the intermittent jointed rock mass, namely linear elastic stage, sliding of joint, steady growth of post-crack, joint coalescence failure, and residual strength. According to strength analysis, the failure envelopes of intact rock and across jointed rock mass are the upper bound and lower bound separately. The strength of intermittent jointed rock mass can be evaluated by reducing the bandwidth of the failure envelope with geo-mechanics analysis. (4) Numerical test of rock mass Two sets of methods, i.e. the distinct element method (DEC) based on in-situ geology mapping and the realistic failure process analysis (RFPA) based on high-definition digital imaging, are developed and introduced. The operation process and analysis results are demonstrated detailedly from the research on parameters of rock mass based on numerical test in the Jinping First Stage Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Then the applicable fields are figured out respectively. (5) Intelligent evaluation based on artificial neural network (ANN) The characters of both ANN and parameter evaluation of rock mass are discussed and summarized. According to the investigations, ANN has a bright application future in the field of parameter evaluation of rock mass. Intelligent evaluation of mechanical parameters in the Jinping First Stage Hydropower Station is taken as an example to demonstrate the analysis process. The problems in five aspects, i. e. sample selection, network design, initial value selection, learning rate and expected error, are discussed detailedly.

大型反倾岩体变形破坏机制及稳定性分析———以金沙江宽谷坝址龙蟠边坡为例

A large number of catastrophic accidents were aroused by the instability and destruction of anti-dip rock masses in the worldwide engineering projects, such as hydropower station, mine, railways and so on. Problems in relation to deformation and failure about anti-dip rock slopes are significant for engineering geology research. This dissertation takes the Longpan slope in the Jinsha River as a case to study the deformation mechanism of large-scale anti-dip rock masses and the slope stability analysis method. The primary conclusions are as follows. The Dale Reach of Jinsha River, from Longpan to the debouchment of Chongjiang tributary, is located in the southeastern margin of the Qinghai-Tibet Plateau. Longpan slope is the right embankment of Dale dam, it is only 26 km to the Shigu and 18 km to Tiger Leaping Gorge. The areal geology tectonic structures here area are complicated and blurry. Base on the information of geophysical exploration (CSAMT and seismology) and engineering geological investigation, the perdue tectonic pattern of Dale Reach is put forward for the first time in this paper. Due to the reverse slip of Longpan fault and normal left-rotation of Baihanchang fault, the old faulted valley came into being. The thick riverbed sediments have layered characters of different components and corresponding causes, which attribute to the sedimentary environments according with the new tectonic movements such as periodic mountain uplifting in middle Pleistocene. Longpan slope consists of anti-dip alternate sandstone and slate stratums, and the deformable volume is 6.5×107m3 approximately. It was taken for an ancient landslide or toppling failure in the past so that Dale dam became a vexed question. Through the latest field surveying, displacement monitoring and rock masses deforming characters analyses, the geological mechanism is actually a deep-seated gravitational bending deformation. And then the discrete element method is used to simulate the deforming evolution process, the conclusion accords very well with the geo-mechanical patterns analyses. In addition strength reduction method based on DEM is introduced to evaluate the factor of safety of anti-dip rock slope, and in accordance with the expansion way of the shear yielding zones, the progressive shear failure mechanism of large-scale anti-dip rock masses is proposed for the first time. As an embankment or a close reservoir bank to the lower dam, the stability of Longpan slope especially whether or not resulting in sliding with high velocity and activating water waves is a key question for engineering design. In fact it is difficult to decide the unified slip surface of anti-dip rock slope for traditional methods. The author takes the shear yielding zones acquired form the discrete element strength reduction calculation as the potential sliding surface and then evaluates the change of excess pore pressure and factor of stability of the slope generated by rapid drawdown of ponded water. At the same time the dynamic response of the slope under seismic loading is simulated through DEM numerical modeling, the following results are obtained. Firstly the effective effect of seismic inertia force is resulting in accumulation of shear stresses. Secondly the discontinuous structures are crucial to wave transmission. Thirdly the ultimate dynamic response of slope system takes place at the initial period of seismic loading. Lastly but essentially the effect of earthquake load to bringing on deformation and failure of rock slope is the coupling effect of shear stresses and excess pore water pressure accumulation. In view of limitations in searching the critical slip surface of rock slope of the existing domestic and international software for limit equilibrium slope stability analyses, this article proposes a new method named GA-Sarma Algorithm for rock slope stability analyses. Just as its name implies, GA-Sarma Algorithm bases on Genetic Algorithm and Sarma method. GA-Sarma Algorithm assumes the morphology of slip surface to be a broken line with traceability to extend along the discontinuous surface structures, and the slice boundaries is consistent with rock mass discontinuities such as rock layers, faults, cracks, and so on. GA-Sarma Algorithm is revolutionary method that is suitable for global optimization of the critical slip surface for rock slopes. The topics and contents including in this dissertation are closely related to the difficulties in practice, the main conclusions have been authorized by the engineering design institute. The research work is very meaningful and useful for the engineering construction of Longpan hydropower station.

溶蚀岩体渗透水力学特性研究

Seepage control in karstic rock masses is one of the most important problems in domestic hydroelectric engineering and mining engineering as well as traffic engineering. At present permeability assessment and leakage analysis of multi-layer karstic rock masses are mainly qualitative, while seldom quantitative. Quantitative analyses of the permeability coefficient and seepage amount are conducted in this report, which will provide a theoretical basis for the study of seepage law and seepage control treatment of karstic rocks. Based on the field measurements in the horizontal grouting galleries of seepage control curtains on the left bank of the Shuibuya Hydropower Project on the Qingjiang river, a hydraulic model is established in this report, and the computation results will provide a scientific basis for optimization of grouting curtain engineering. Following issues are addressed in the report. (1) Based on the in-situ measurements of fissures and karstic cavities in grouting galleries, the characteristics of karstic rock mass is analyzed, and a stochastic structural model of karstic rock masses is set up, which will provide the basis for calculation of the permeability and leakage amount of karstic rock mass. (2) According to the distribution of the measured joints in the grouting galleries and the stochastic results obtained from the stochastic structural model of karstic rock mass between grouting galleries, a formula for computation of permeability tensor of fracturing system is set up, and a computation program is made with Visual Basic language. The computation results will be helpful for zoning of fissured rock masses and calculation of seepage amount as well as optimization of seepage control curtains. (3) Fractal theory is used to describe quantitatively the roughness of conduit walls of karstic systems and the sinuosity of karstic conduits. It is proposed that the roughness coefficient of kastic caves can be expressed by both fractal dimension Ds and Dr that represent respectively the extension sinuosity of karstic caves and the roughness of the conduit walls. The existing formula for calculating the seepage amount of karstic conduits is revised and programmed. The seepage amount of rock masses in the measured grouting galleries is estimated under the condition that no seepage control measures are taken before reservoir impoundment, and the results will be helpful for design and construction optimization of seepage curtains of the Shuibuya hydropower project. This report is one part of the subject "Karstic hydrogeology and the structural model and seepage hydraulics of karstic rock masses", a sub-program of "Study on seepage hydraulics of multi-layer karstic rock masses and its application in seepage control curtain engineering", which is financially supported by the Hubei Provincial key science and technology programme.

Synchronous Multi-Scale Observations on Rock Damage and Rupture

This paper reports a multi-scale study on damage evolution process and rupture of gabbro under uniaxial compression with several experimental techniques, including MTS810 testing machine, white digital speckle correlation method, and acoustic emission technique. In particular, the synchronization of the three experimental systems is realized for the study of relationship of deformation and damage at multiple scales. It is found that there are significant correlation between damage evolution at small and large length scales, and rupture at sample scale, especially it displays critical sensitivity at multiple scales and trans-scale fluctuations.

Acoustic Emission Experiments of Rock Failure under Load Simulating the Hypocenter Condition

A series of acoustic emission (AE) experiments of rock failure have been conducted under cyclic load in tri-axial stress tests. To simulate the hypocenter condition the specimens are loaded by the combined action of a constant stress, intended to simulate

Numerical Study on Rainfall Infiltration in Rock-Soil Blocks

A mathematical model for the rain infiltration in the rock-soil slop has been established and solved by using the finite element method. The unsteady water infiltrating process has been simulated to get water content both in the homogeneous and heterogeneous media. The simulated results show that the rock blocks in the rock-soil slop can cause the wetting front moving fast. If the rain intensity is increased, the saturated region will be formed quickly while other conditions are the same. If the rain intensity keeps a constant, it is possible to accelerate the generation of the saturated region by properly increasing the vertical filtration rate of the rock-soil slop. However, if the vertical filtration rate is so far greater than the rain intensity, it will be difficult to form the saturated region in the rock-soil slop. The numerical method was verified by comparing the calculation results with the field test data.

3D Anisotropic Elastoplastic-Damage Model and its Application in Simulating the Behavior of Rock Materials

A 3D anisotropic elastoplastic-damage model was presented based on continuum damage mechanics theory. In this model, the tensor decomposition technique is employed. Combined with the plastic yield rule and damage evolution, the stress tensor in incremental format is obtained. The derivate eigenmodes in the proposed model are assumed to be related with the uniaxial behavior of the rock material. Each eigenmode has a corresponding damage variable due to the fact that damage is a function of the magnitude of the eigenstrain. Within an eigenmodes, different damage evolution can be used for tensile and compressive loadings. This model was also developed into finite element code in explicit format, and the code was integrated into the well-known computational environment ABAQUS using the ABAQUS/Explicit Solver. Numerical simulation of an uniaxial compressive test for a rock sample is used to examine the performance of the proposed model, and the progressive failure process of the rock sample is unveiled.

Mechanism on Progressive Failure of a Faulted Rock Slope Due to Slip-Weakening

Slip-weakening is one of the characteristics of geological materials under certain loadings. Non-uniform rock structure may exist in the vicinity of the slip surface for a rock slope. Some portion of the slip surface may be penetrated but the other not. For the latter case, the crack or the fault surface will undergo shear deformation before it becomes a successive surface under a certain loading. As the slipped portion advances,slip-weakening occurs over a distance behind the crack tip. In the weakening zone, the shear strength will decrease from its peak value to residual friction level. The stress will redistribute along the surface of crack and in the weakening zone. Thus the changed local stress concentration leads the crack to extend and the ratio of penetration of the slip surface to increase. From the view of large-scale for the whole slip surface, the shear strength will decrease due to the damage of interior rock structure, and the faulted rock behaves as a softening material. Such a kind of mechanism performs in a large number of practical landslides in the zones experienced strong earthquakes. It should be noted that the mechanism mentioned above is different from that of the breakage of structural clay,in which the geological material is regarded as a medium containing structural lumps and structural bands. In this paper, the softening behavior of a faulted rock should be regarded as a comprehensive result of the whole complicated process including slip-weakening, redistribution of stress, extension of crack tip, and the penetration of the slip surface. This process is accompanied by progressive failure and abrupt structural damage. The size of slip-weakening zone is related to the undergoing strain. Once the relative slide is initiated (local or integrated), the effect of slip-weakening will behave in a certain length behind the crack tip until the formation of the whole slip surface.

Experimental study on pore pressure in rock-soil slope during reservoir water level fluctuation

A test system was developed for measuring the pore pressure in porous media, and a new model was devised for the pore pressure testing in both saturated and unsaturated rock-soil. Laboratory experiments were carried out to determine the pore pressure during water level fluctuation. The variations of transient pore pressure vs. time at different locations of the simulated rock-soil system were acquired and processed, and meanwhile the deformation and failure of the model are observed. The experiment results show that whether the porous media are saturated or not, the transient pore pressure is mainly dependent on the water level fluctuation, and coupled with the variation of the stress field.

Numerical simulation of rock failure and earthquake process on mesoscopic scale

On the basis of the lattice model of MORA and PLACE, Discrete Element Method, and Molecular Dynamics approach, another kind of numerical model is developed. The model consists of a 2-D set of particles linked by three kinds of interactions and arranged into triangular lattice. After the fracture criterion and rules of changes between linking states are given, the particle positions, velocities and accelerations at every time step are calculated using a finite-difference scheme, and the configuration of particles can be gained step by step. Using this model, realistic fracture simulations of brittle solid (especially under pressure) and simulation of earthquake dynamics are made.

Distinct Element Analysis on Propagation Characteristics of P-Wave in Rock Pillar with Finite length

以节理岩体等效刚度的概念为基础,讨论了离散元刚性块体模型中节理刚度的选取问题。采用面-面接触模型模拟了纵波在一维岩体中的传播,给出了纵波波形；研究了阻尼比、软弱夹层以及节理间是否可拉对波传播规律的影响。