基于岩体精细化描述的围岩分类及力学参数概率分布特征分析

针对目前水电站地下厂房工程中不同围岩分类方法存在评价结果不一致、围岩力学参数存在室内试验值与实际情况不吻合的现象,现推荐采用岩体精细化描述体系对围岩岩体结构进行定量化评价。将常用围岩分类方法(RMR、Q、RMi、GSI、BQ、HC)评价指标予以归纳分组,并通过各组内不同指标对比分析获得围岩分类方法中的基础评价指标。以大岗山水电站主厂房某区段为分析对象,采取现场岩体精细化地质素描与后期数据挖掘、拟合相结合方法,并依据评价指标间的关联关系,获得了基础、非基础评价指标的分布概型及对应参数,实现对该段围岩岩体精细化描述认知;基于精细化描述结果,应用Monte Carlo法生成符合各评价指标分布概型的大量随机数,而后参照各分类方法评价思路与评分流程,得到评价指标在各分类方法对应的大量随机评分值,通过归纳统计获得不同围岩分类方法评价结果的分布概型;基于各围岩分类方法评价结果与力学参数值之间的关联关系实现对力学参数概率特征分析。该分析方法与思路可为类似工程围岩质量及力学参数的精确确定提供一定借鉴,并可为实现围岩支护极限状态设计提供必要的原始参数支持。     

Three-dimensional numerical modelling of falling rock protection barriers

This paper presents a new numerical strategy for the design and verification of flexible falling rock barriers: passive protection measures for risk mitigation of potentially unstable rock slopes. The key point of the proposed approach is that notwithstanding the complexity of the simulated phenomenon, the resulting highly non-linear, dynamic model is simple and produces an accurate prediction of all the relevant parameters for barrier design, such as anchorage forces, net panel elongations and residual heights.The modelling procedure has been assessed using detailed experimental data obtained from a set of full-scale tests on three barrier prototypes with various energy absorption capacities (5000 kJ, 3000 kJ and 500 kJ). By comparison with the experimental results, the numerical model has shown to be reliable in capturing very accurately the barrier response to a block impact. Consequently, this method can be extended to investigate the behaviour of flexible falling rock protection barriers under conditions different from those encountered in full-scale tests. Therefore, the numerical procedure can be regarded as an effective tool used for designing and testing these structures.     

New Methodology to Characterize Shear Behavior of Joints by Combination of Direct Shear Box Testing and Numerical Simulations

A new procedure is presented, which combines big shear box tests on rocks and corresponding numerical simulations with explicit consideration of joint roughness to get deeper insight into the shear behavior of rock joints. The procedure consists of three parts: (1) constant normal load- or CNS-shear box tests with registration of shear- and normal-components of stress and displacements and deduction of basis rock mechanical parameters; (2) high resolution 3D-scanning of joint surface to deduce joint topography; and (3) set-up, run and evaluation of 3-dimensional numerical model with explicit duplication of joint roughness as back-analysis of shear box tests. The numerical back-analysis provides deeper insight into the joint behavior at the micro-scale. Several parameters can be deduced, like micro-slope angle distribution, aperture size distribution, local normal stress distribution and detailed analysis of dilation in relation to shear direction. The potential of the new procedure is illustrated exemplary by shear box tests on slate.     

Prediction of Drillability of Rocks with Strength Properties Using a Hybrid GA-ANN Technique

The purpose of this paper is to provide a proper, practical and convenient drilling rate index (DRI) prediction model based on rock material properties. In order to obtain this purpose, 47 DRI tests were used. In addition, the relevant strength properties i.e. uniaxial compressive strength and Brazilian tensile strength were also used and selected as input parameters to predict DRI. Examined simple regression analysis showed that the relationships between the DRI and predictors are statistically meaningful but not good enough for DRI estimation in practice. Moreover, multiple regression, artificial neural network (ANN) and hybrid genetic algorithm (GA)-ANN models were constructed to estimate DRI. Several performance indices i.e. coefficient of determination (R²), root mean square error and variance account for were used for evaluation of performance prediction the proposed methods. Based on these results and the use of simple ranking procedure, the best models were chosen. It was found that the hybrid GA-ANN technique can performed better in predicting DRI compared to other developed models. This is because of the fact that the proposed hybrid model can update the biases and weights of the network connection to train by ANN.     

Evaluation of spatial variation of peak horizontal acceleration and spectral acceleration for south India: a probabilistic approach

In this work, an attempt has been made to evaluate the spatial variation of peak horizontal acceleration (PHA) and spectral acceleration (SA) values at rock level for south India based on the probabilistic seismic hazard analysis (PSHA). These values were estimated by considering the uncertainties involved in magnitude, hypocentral distance and attenuation of seismic waves. Different models were used for the hazard evaluation, and they were combined together using a logic tree approach. For evaluating the seismic hazard, the study area was divided into small grids of size 0.1° × 0.1°, and the hazard parameters were calculated at the centre of each of these grid cells by considering all the seismic sources within a radius of 300 km. Rock level PHA values and SA at 1 s corresponding to 10% probability of exceedance in 50 years were evaluated for all the grid points. Maps showing the spatial variation of rock level PHA values and SA at 1 s for the entire south India are presented in this paper. To compare the seismic hazard for some of the important cities, the seismic hazard curves and the uniform hazard response spectrum (UHRS) at rock level with 10% probability of exceedance in 50 years are also presented in this work.     

Recent insights into analytical precision and modelling of DDA and NMM for practical problems

This paper discussed the analytical precision on equations of motion in some practical studies that use 2-D DDA, 2D-Coupled analysis of NMM & DDA, and 3D-DDA. We point out the main reason for numerical instability in DDA is loss of the effective digits when superposing the global stiffness matrix under the variable conditions of large and the small blocks. Categories of analyses are distributed in 2D-DDA static masonry structures such as the Great Pyramid and the Pont du Gard, 2D-DDA dynamic rock fall problems and the sensitivity of its analysis on the velocity ratio and the penalty, earthquake response analysis of rock slopes, 2D-NMM, 2D-coupled analysis of DDA and NMM and 3D-DDA rock fall problem. The selected examples on the Great Pyramid, Pont du Gard and Bayon Temple at Angkor Thom are located in Egypt, France and Cambodia respectively as UNESCO world heritages sites. The authors evaluated the applicability of the methods and the reliability of the results by comparing different methods and site observations from the practical problems.     

Physical and Numerical Modeling of Seismic Soil-Structure Interaction in Layered Soils

The structural response of buildings subjected to seismic loads is affected by local site conditions and the interaction between the structure and the supporting soil media. Seismic centrifuge model tests were conducted on two layered clay soil profiles at 80 g field to investigate soil-structure interaction and dynamic response of foundation. Several earthquake-like shaking events were applied to the models using an electro-hydraulic shaking table to simulate linear and nonlinear soil behavior. Results showed that the foundation input motion was significantly amplified in both models, especially for weak earthquake motions. Seismic soil-structure interaction was found to have an important effect on structure response by increasing the amplification of foundation input motion. A 3D finite difference numerical model was also developed to simulate the response of centrifuge model tests and study the parameters that affect the characteristics of earthquake at the base of the structure. The results indicated that the stiffness and stratification of the soil profiles had a significant effect on modifying the foundation input motion.     

Applicability of existing models to predict the behavior of replicas of natural fractures of welded tuff under different boundary conditions

Assessing the shear behavior of intact rock and rock fractures is an important issue in the design of a potential nuclear waste repository at Yucca Mountain, Nevada. Cyclic direct shear experiments were conducted on replicas of three natural fractures and a laboratory-developed tensile fracture of welded tuff. The tests were carried out under constant normal loads or constant normal stiffnesses with different initial normal load levels. Each test consisted of five cycles of forward and reverse shear motion. In this paper, the results of the constant normal load shear experiments are analyzed using several constitutive models proposed in the rock mechanics literature for joint shear strength, dilatancy, and joint surface damage. It is shown that some of the existing models have limitations. New constitutive models are proposed and are included in a mathematical analysis tool that can be used to predict joint behavior under various boundary conditions. © Rapid Science Ltd. 1998     

岩石力学系统运动稳定性问题及其研究现状

本文从岩石工程破坏过程稳定性的实际出发，提出了岩石力学系统及其运动稳定性的概念，分析了岩石力学系统的特点和岩石力学系统运动稳定性所要研究的内容，比较全面地总结了岩石力学系统稳定性研究的现状，提出在今后的岩石力学系统稳定性的研究中，要从系统运动的角度考虑，对各类稳定性问题可采用统一的研究思路，并对岩石力学系统中的各种参数进行研究，进而分析其对系统稳定性的影响。     

岩石物理建模技术在致密砂岩储层预测中的应用——以鄂尔多斯盆地北部H区块为例

岩石物理建模是叠前反演前重要的基础工作。以致密砂岩储层为目标,剖析了岩石物理建模涉及的主要过程。将测井评价与岩石物理建模视为迭代过程,测井评价提供岩石物理建模的输入数据,岩石物理建模控制测井评价参数的质量。在测井参数具有较高质量的情况下,优选岩石物理建模方法,优化关键参数,较为准确地预测了横波速度等弹性参数,为叠前弹性参数反演提供了可靠的基础数据。利用岩石物理解释模板对叠前反演纵横波速度体进行了解释,提高了致密砂岩气层的识别精度,在鄂尔多斯盆地取得了较好的应用效果。     

Full-scale Modelling of Falling Rock Protection Barriers

Full-scale impact tests, carried out to evaluate the behaviour of flexible falling rock protection barriers, are described and relevant results presented and discussed. Falling rock protection barriers, which may be numbered among passive measures against rockfall, are designed to intercept and stop falling rocks by dissipation of impact energies through the elasto-plastic deformation of a system made up of metallic nets and supporting and connecting components. The testing programme involved models of barriers subjected to the impact of free-falling blocks of kinetic energy ranging from 500 to 5,000 kJ. The experimental test site, set-up in Fonzaso (Italy), and the experimental procedure were developed according to the new European testing standards (ETAG 027) on falling rock protection kits. The paper is aimed at presenting an extensive and high quality database, which can be extremely useful for a better understanding of the actual response of such structures and for any subsequent analytical and numerical modelling.     

A deterministic seismic hazard analysis for shallow earthquakes in Greece

The maximum expected ground motion in Greece is estimated for shallow earthquakes using a deterministic seismic hazard analysis (DSHA). In order to accomplish this analysis the input data include an homogeneous catalogue of earthquakes for the period 426 BC–2003, a seismogenic source model with representative focal mechanisms and a set of velocity models. Because of the discrete character of the earthquake catalogue and of errors in location of single seismic events, a smoothing algorithm is applied to the catalogue of the main shocks to get a spatially smoothed distribution of magnitude. Based on the selected input parameters synthetic seismograms for an upper frequency content of 1 Hz are computed on a grid of 0.2° × 0.2°. The resultant horizontal components for displacement, velocity, acceleration and DGA (Design Ground Acceleration) are mapped. The maps which depict these results cannot be compared with previously published maps based on probabilistic methodologies as the latter were compiled for a mean return period of 476 years. Therefore, in order to validate our deterministic analysis, the final results are compared with PGA estimated from the maximum observed macroseismic intensity in Greece during the period 426 BC–2003.Since the results are obtained for point sources, with the frequency content scaled with moment magnitude, some sensitivity tests are performed to assess the influence of the finite extent of fault related to large events. Sensitivity tests are also performed to investigate the changes in the peak ground motion quantities when varying the crustal velocity models in some seismogenic areas. The ratios and the relative differences between the results obtained using different models are mapped and their mean value computed. The results highlight the importance in the deterministic approach of using good and reliable velocity models.     

A Moisture Index Classification System for Rocks (MiC System)

Summary.  Interaction of rocks with water is leading to a reduction in physical and mechanical properties, and quantifying this effect has always been a problem for design purposes in rock engineering projects. Despite the fact that a number of considerable works have been carried out in this field, they are not thorough and general. In this paper, a new rock classification system is presented in order to evaluate wet rock properties. Laboratory studies showed that six parameters may be recognized as major controlling parameters for wet rocks, including: water absorption, porosity, weathering degree, swelling index, number of microfissures, and average grain size. Many tests were performed and a sensitivity analysis was carried out on the data obtained. This analysis characterized the rock sensitivity when changing each parameter and consequently, a rating value was assigned to each parameter. By adding these rating values, a rating between zero and 100 is allocated to each rock as Moisture index (Mi). The closer this index is to 100, the more sensitive rock will be to moisture. On the basis of this index, rocks are divided into five classes including: very good, good, fair, weak and very weak.     

Vacuum testing of the permeability of the excavation damaged zone

Summary The vacuum permeability method provides a new approach to assessment of permeability in the near-field of rock excavations. The method, which resembles leak testing of vacuum lines, is based on the measurement and interpretation of the pneumatic response to evacuation of a rock mass. As the induced drainage area is in order of meters, the method is particularly suitable on the scale of the excavation damaged zone. The data can be evaluated either analytically or numerically. The evaluation is similar to the approach used in the analysis of drill stem tests in gas wells. Experiments on large scale laboratory models corroborated the validity of the method and the evaluation procedure. Field tests were conducted at a wide range of geomechanical conditions. The collected data indicate that there are three characteristic types of permeability distribution in the excavation damaged zone.     

岩体力学试验的理论分析

本文通过对岩石力学试验中岩石试件破坏方式的分析, 在理论上解释了岩体力学试验中经常遇到的一些力学现象并且作出了定量的回答, 同时还建立了岩石的各种破坏模式的判据。     

Probabilistic prediction of the spatial distribution of potential key blocks during tunnel surrounding rock excavation

Jointed network simulations tend to be more random in nature due to the uncertainty of rock mass structures. In this paper, a series of jointed network models can be established in batches using Monte Carlo simulation (MSC) and loop iteration. Taking the joints, tunnel profile and their intersections as the edges E and vertices V of graph G, the jointed network model can serve as an unweighted undigraph. Then, the breadth-first search is introduced to search the closed paths around the tunnel profile, such as the potential key blocks. With batch simulation of network models, the spatial distribution characteristics and probability distribution rules of blocks can be automatically analysed during the search process. For comparison, the Laohushan tunnel of the Jinan Belt Expressway in China has been analysed using the breadth-first search, discontinuous deformation analysis method and procedure of “Finding the Key Blocks-Unrolled Tunnel Joint Trace Maps”. Each simulation starts from the same probabilistic model of geometrical parameters of joints but develops differently with different outcomes. The spatial distribution rule of potential key blocks simulated by the aforementioned batch jointed network models is essentially identical to the actual rockfall during tunnel excavation.

     

Single‐and multi‐objective genetic algorithm optimization for identifying soil parameters

This paper discusses the quality of the procedure employed in identifying soil parameters by inverse analysis. This procedure includes a FEM‐simulation for which two constitutive models—a linear elastic perfectly plastic Mohr–Coulomb model and a strain‐hardening elasto‐plastic model—are successively considered. Two kinds of optimization algorithms have been used: a deterministic simplex method and a stochastic genetic method. The soil data come from the results of two pressuremeter tests, complemented by triaxial and resonant column testing. First, the inverse analysis has been performed separately on each pressuremeter test. The genetic method presents the advantage of providing a collection of satisfactory solutions, among which a geotechnical engineer has to choose the optimal one based on his scientific background and/or additional analyses based on further experimental test results. This advantage is enhanced when all the constitutive parameters sensitive to the considered problem have to be identified without restrictions in the search space. Second, the experimental values of the two pressuremeter tests have been processed simultaneously, so that the inverse analysis becomes a multi‐objective optimization problem. The genetic method allows the user to choose the most suitable parameter set according to the Pareto frontier and to guarantee the coherence between the tests. The sets of optimized parameters obtained from inverse analyses are then used to calculate the response of a spread footing, which is part of a predictive benchmark. The numerical results with respect to both the constitutive models and the inverse analysis procedure are discussed. Copyright © 2011 John Wiley & Sons, Ltd.     

基于现场试验的岩体变形模量尺寸效应研究

岩体力学参数具有随机性和不确定性,合理确定岩体力学参数是岩土工程稳定性分析的基础和关键。以某水电站前期研究的110个钻孔弹模测试数据为研究样本,对其进行了统计特征分析,确定变形模量的概率分布特征更符合Weibull分布。以钻孔弹模仪承压板尺寸为岩体变形模量研究的样本尺度,岩体由该尺度的岩样叠加而成,建立了不同尺寸的立方体岩体数值分析模型,按照Weibull分布随机给出各岩样的变形模量,在此基础上对不同尺寸的岩体分别进行5 000次单轴压缩试验的数值模拟,得到了不同尺度岩体等效变形模量的变化特征,确定了岩体等效变形模量的REV。基于现场实测岩体变形模量与概率统计理论相结合研究岩体变形模量是一个新的尝试,可为该水电站工程岩体力学参数的合理取值提供依据。     

Effects of Rock Properties on Specific Cutting Energy in Linear Cutting of Sandstones by Picks

Summary. Specific cutting energy (SE) has been widely used to assess the rock cuttability for mechanical excavation purposes. Some prediction models were developed for SE through correlating rock properties with SE values. However, some of the textural and compositional rock parameters i.e. texture coefficient and feldspar, mafic, and felsic mineral contents were not considered. The present study is to investigate the effects of previously ignored rock parameters along with engineering rock properties on SE. Mineralogical and petrographic analyses, rock mechanics, and linear rock cutting tests were performed on sandstone samples taken from sites around Ankara, Turkey. Relationships between SE and rock properties were evaluated using bivariate correlation and linear regression analyses. The tests and subsequent analyses revealed that the texture coefficient and feldspar content of sandstones affected rock cuttability, evidenced by significant correlations between these parameters and SE at a 90% confidence level. Felsic and mafic mineral contents of sandstones did not exhibit any statistically significant correlation against SE. Cementation coefficient, effective porosity, and pore volume had good correlations against SE. Poisson’s ratio, Brazilian tensile strength, Shore scleroscope hardness, Schmidt hammer hardness, dry density, and point load strength index showed very strong linear correlations against SE at confidence levels of 95% and above, all of which were also found suitable to be used in predicting SE individually, depending on the results of regression analysis, ANOVA, Student’s t-tests, and R² values. Poisson’s ratio exhibited the highest correlation with SE and seemed to be the most reliable SE prediction tool in sandstones.