Hydrogeomechanics for rock engineering: coupling subsurface hydrogeomechanical assessement and hydrogeotechnical mapping on fracturated rock masses

The present work aims to achieve and further develop a hydrogeomechanical approach in Caldas da Cavaca hydromineral system rock mass (Aguiar da Beira, NW Portugal), and contribute to a better understanding of the hydrogeological conceptual site model. A collection of several data, namely geology, hydrogeology, rock and soil geotechnics, borehole hydraulics and hydrogeomechanics, was retrieved from three rock slopes (Lagoa, Amores and Cancela). To accomplish a comprehensive analysis and rock engineering conceptualisation of the site, a multi‐technical approach were used, such as, field and laboratory techniques, hydrogeotechnical mapping, hydrogeomechanical zoning and hydrogeomechanical scheme classifications and indexes. In addition, a hydrogeomechanical data analysis and assessment, such as Hydro‐Potential (HP)‐Value technique, JW Joint Water Reduction index, Hydraulic Classification (HC) System were applied on rock slopes. The hydrogeomechanical zone HGMZ 1 of Lagoa slope achieved higher hydraulic conductivities with poorer rock mass quality results, followed by the hydrogeomechanical zone HGMZ 2 of Lagoa slope, with poor to fair rock mass quality and lower hydraulic parameters. In addition, Amores slope had a fair to good rock mass quality and the lowest hydraulic conductivity. The hydrogeomechanical zone HGMZ 3 of Lagoa slope, and the hydrogeomechanical zones HGMZ 1 and HGMZ 2 of Cancela slope had a fair to poor rock mass quality but were completely dry. Geographical Information Systems (GIS) mapping technologies was used in overall hydrogeological and hydrogeomechanical data integration in order to improve the hydrogeological conceptual site model.

Delineating environmental groundwater vulnerability and protection zones mapping in fractured rock masses

Hard‐rock watersheds commonly exhibit complex geological bedrock and morphological features. Hydromineral resources have relevant economic value for the thermal spas industry. The present study aims to develop a groundwater vulnerability approach in Caldas da Cavaca hydromineral system (Aguiar da Beira, Central Portugal) which has a thermal tradition that dates back to the late 19th century, and contribute to a better understanding of the hydrogeological conceptual site model. In this work different layers were overlaid, generating several thematic maps to arrive at an integrated framework of several key‐sectors in Caldas da Cavaca site. Thus, to accomplish a comprehensive analysis and conceptualization of the site, a multi‐technical approach was used, such as, field and laboratory techniques, where several data was collected, like geotectonics, hydrology and hydrogeology, hydrogeomorphology, hydrogeophysical and hydrogeomechanical zoning aiming the application of the so‐called DISCO method. All these techniques were successfully performed and a groundwater vulnerability to contamination assessment, based on GOD‐S, DRASTIC‐Fm, SINTACS, SI and DISCO indexes methodology, was delineated. Geographical Information Systems (GIS) technology was on the basis to organise and integrate the geodatabases and to produce all the thematic maps. This multi‐technical approach highlights the importance of groundwater vulnerability to contamination mapping as a tool to support hydrogeological conceptualisation, contributing to better decision‐making of water resources management and sustainability.

Comparisons between two empirical yield criteria for rock masses

In general, rock masses are inhomogeneous, discontinuous media composed of rock material and naturally occurring discontinuities such as joints, fractures and bedding planes. Because of these features, the strength of rock masses is notoriously difficult to assess. Nonetheless, many criteria have been proposed for estimating rock mass strength. Based on the finite element upper and lower bound limit analysis methods, this study examined two empirical yield criteria for rock masses, the Hoek-Brown failure criterion (2002) and the Douglas criterion (2002). The comparisons showed that very different results may be obtained using the same input parameters. Therefore, it is interesting to discuss the source of these differences.

Discontinuity spacing analysis in rock masses using 3D point clouds

The complete characterization of rock masses implies the acquisition of information of both, the materials which compose the rock mass and the discontinuities which divide the outcrop. Recent advances in the use of remote sensing techniques – such as Light Detection and Ranging (LiDAR) – allow the accurate and dense acquisition of 3D information that can be used for the characterization of discontinuities. This work presents a novel methodology which allows the calculation of the normal spacing of persistent and non-persistent discontinuity sets using 3D point cloud datasets considering the three dimensional relationships between clusters. This approach requires that the 3D dataset has been previously classified. This implies that discontinuity sets are previously extracted, every single point is labeled with its corresponding discontinuity set and every exposed planar surface is analytically calculated. Then, for each discontinuity set the method calculates the normal spacing between an exposed plane and its nearest one considering 3D space relationship. This link between planes is obtained calculating for every point its nearest point member of the same discontinuity set, which provides its nearest plane. This allows calculating the normal spacing for every plane. Finally, the normal spacing is calculated as the mean value of all the normal spacings for each discontinuity set. The methodology is validated through three cases of study using synthetic data and 3D laser scanning datasets. The first case illustrates the fundamentals and the performance of the proposed methodology. The second and the third cases of study correspond to two rock slopes for which datasets were acquired using a 3D laser scanner. The second case study has shown that results obtained from the traditional and the proposed approaches are reasonably similar. Nevertheless, a discrepancy between both approaches has been found when the exposed planes members of a discontinuity set were hard to identify and when the planes pairing was difficult to establish during the fieldwork campaign. The third case study also has evidenced that when the number of identified exposed planes is high, the calculated normal spacing using the proposed approach is minor than those using the traditional approach.

Datasets and results of the research [Discontinuity spacing analysis in rock masses using 3D point clouds]

.bin files should be opened using CloudCompare

Proposal of a New Parameter for the Weathering Characterization of Carbonate Flysch-Like Rock Masses: The Potential Degradation Index (PDI)

The susceptibility of clay bearing rocks to weathering (erosion and/or differential degradation) is known to influence the stability of heterogeneous slopes. However, not all of these rocks show the same behaviour, as there are considerable differences in the speed and type of weathering observed. As such, it is very important to establish relationships between behaviour quantified in a laboratory environment with that observed in the field. The slake durability test is the laboratory test most commonly used to evaluate the relationship between slaking behaviour and rock durability. However, it has a number of disadvantages; it does not account for changes in shape and size in fragments retained in the 2 mm sieve, nor does its most commonly used index (Id2) accurately reflect weathering behaviour observed in the field. The main aim of this paper is to propose a simple methodology for characterizing the weathering behaviour of carbonate lithologies that outcrop in heterogeneous rock masses (such as Flysch slopes), for use by practitioners. To this end, the Potential Degradation Index (PDI) is proposed. This is calculated using the fragment size distribution curves taken from material retained in the drum after each cycle of the slake durability test. The number of slaking cycles has also been increased to five. Through laboratory testing of 117 samples of carbonate rocks, extracted from strata in selected slopes, 6 different rock types were established based on their slaking behaviour, and corresponding to the different weathering behaviours observed in the field.

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.

Application of neural network to rock slope stability assessments

It is known that rock masses are inhomogeneous, discontinuous media composed of rock material and naturally occurring discontinuities such as joints, fractures and bedding planes. These features make any analysis very difficult using simple theoretical solutions. Generally speaking, back analysis technique can be used to capture some implicit parameters for geotechnical problems. In order to perform back analyses, the procedure of trial and error is generally required. However, it would be time-consuming. This study aims at applying a neural network to do the back analysis for rock slope failures. The neural network tool will be trained by using the solutions of finite element upper and lower bound limit analysis methods. Therefore, the uncertain parameter can be obtained, particularly for rock mass disturbance.

Mathematical development of stress-strain law for rock block contacts

The fracture behavior of rock block contacts has been studied for many years. Unfortunately, up to now, there is not a rigorous formulation or a solid theoretical foundation to support it. A mathematical development to represent the failure mechanism which occurs in the contacts between rock blocks is presented to evaluate the performance of breaking mechanism of such blocks relating it to the morphology of the contact and mechanical parameters of the material. The examined framework includes the evaluation of the surface roughness of first order in the failure mechanism of the granular particles of large size and the development of a theoretical model describing the morphology of the contact between rock blocks.

Characterization of rock joints by fractal analysis

From a physical perspective, a joint experiences fracturing processes that affect the rock at both microscopic and macroscopic levels. The result is a behaviour that follows a fractal structure. In the first place, for saw-tooth roughness profiles, the use of the triadic Koch curve appears to be adequate and by means of known correlations the JRC parameter is obtained from the angle measured on the basis of the height and length of the roughnesses. Therefore, JRC remains related to the geometric pattern that defines roughness by fractal analysis. In the second place, to characterise the geometry of irregularities with softened profiles, consequently, is proposed a characterisation of the fractal dimension of the joints with a circumference arc generator that is dependent on an average contact angle with regard to the mid-plane. The correlation between the JRC and the fractal dimension of the model is established with a defined statistical ratio.

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

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

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

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

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

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.