Engineering geological studies and tunnel support design at Sulakyurt dam site, Turkey

This paper presents the results of engineering geological investigations and tunnel support design studies, carried out at the Sulakyurt dam site, northeast of Ankara, Turkey. The Sulakyurt dam will be used for flow control and water storage for irrigation projects. Studies were carried out both in the field and the laboratory. Field studies include engineering geological mapping, intensive discontinuity surveying, core drilling and sampling for laboratory testing. The diversion tunnel will be driven in rock mass, consisting of granite and diorite. Empirical, analytical and numerical methods were combined for safe tunnel design. Rock mass rating (RMR), Rock mass quality (Q) and Geological strength index (GSI) systems were used for empirical rock mass quality determination, site characterization and support design. The convergence–confinement method was used as analytical method and software called Phase², a 2D finite element program, was utilized as numerical method. According to the results acquired from the empirical, analytical and numerical methods, tunnel stability problems were expected in both granite and diorite rock masses. The support system, suggested by empirical methods, was applied and the performance of suggested support system was evaluated by means of numerical modelling. It was concluded that the suggested support systems were adequate, since after applying the suggested support system to granite and diorite, tunnel deformation and the yielded elements around the tunnel decreased significantly. Thus, it is suggested that for more reliable support design empirical, numerical and analytical methods should be combined.     

Support capacity estimation of a diversion tunnel in weak rock

This paper presents the results of the support capacity estimation for the diversion tunnel of the Uru dam site in highly weathered tuff and weak zone. Tunneling in weak rock requires some special considerations, since misjudgment in support design results in costly failures. There are several ways of estimating rock support pressure and selecting support. However, all systems suffer from their characteristic limitations in achieving objectives. Thus, it is more useful to use different methods for estimating support pressure and type of support. The support pressure p_i was established by three different methods. These methods are the (1) empirical methods based on rock mass rating (RMR) and rock mass quality index (Q-classification systems), (2) ground support interaction analysis (GSIA) and (3) numerical methods, namely, Phase² finite element (FEM) program. Rock masses were characterized in terms of RSR, RMR, Q-system and GSI. Drill-core samples were tested in the rock mechanics laboratory to determine physico-mechanical properties. Rock mass strength was estimated by empirical methods. Finally, the required support system is proposed and evaluated by different methods in the highly weathered tuff and weak zone of the diversion tunnel.     

Empirical and numerical analyses of support requirements for a diversion tunnel at the Boztepe dam site, eastern Turkey

This paper presents the engineering geological properties and support design of a planned diversion tunnel at the Boztepe dam site that contains units of basalt and tuffites. Empirical, theoretical and numerical approaches were used and compared in this study focusing on tunnel design safety. Rock masses at the site were characterized using three empirical methods, namely rock mass rating (RMR), rock mass quality (Q) and geological strength index (GSI). The RMR, Q and GSI ratings were determined by using field data and the mechanical properties of intact rock samples were evaluated in the laboratory. Support requirements were proposed accordingly in terms of different rock mass classification systems. The convergence–confinement method was used as the theoretical approach. Support systems were also analyzed using a commercial software based on the finite element method (FEM). The parameters calculated by empirical methods were used as input parameters for the FEM analysis. The results from the two methods were compared with each other. This comparison suggests that a more reliable and safe design could be achieved by using a combination of empirical, analytical and numerical approaches.     

超大跨度深埋地下结构围岩压力计算研究

围岩压力是地下结构设计与施工的重要依据。新建京张城际铁路八达岭车站拟定备选方案之一的单拱超大跨设计方案中,最大的跨度近45 m。对于类似超大跨度的深埋地下结构,传统的围岩压力理论与经验计算方法都存在着明显的局限性,已经不再适用。而结合现场勘查资料和室内外试验成果的数值模拟方法,尽管可以得到相对精确的围岩压力计算结果,但存在着计算过程复杂、计算效率低下的不足。因此,在计算分析了多种传统围岩压力计算方法与数值模拟方法获得的围岩压力所存在差异的基础上,对Q系统的围岩压力公式引入跨度项进行了修正,从而获得了能够快速估算超大跨度深埋地下结构围岩压力的估算方法。最后对不同跨度工程的围岩压力进行估算,并与监测值进行比较,对比结果良好。该方法不仅可以为工程设计人员及现场工作人员提供设计计算依据,还可为相关研究人员提供理论参考。     

Proposed support design, Kaletepe tunnel, Turkey

In this paper, preliminary support design of Kaletepe tunnel, located on Bilecik-Istanbul highway, Turkey, was analyzed by empirical and numerical methods. The rock mass rating (RMR) and rock mass quality (Q) systems were employed for empirical rock mass quality determination. Numerical analysis for the stress–strain distribution of the tunnel excavation and support systems was also carried out. The applied support performance was investigated at different sections of the tunnel route. It was seen that empirical and numerical approaches showed similar results. This indicates that when the empirical method is supported by numerical method, the preliminary support design will be more reliable.     

Preliminary support design for Ankara subway extension tunnel

This paper presents the results of preliminary support design of the subway tunnel for Ankara subway project in accordance with some empirical and numerical methods, using the phase 2D finite element method (FEM). The 5 m diameter subway tunnel will advance through slightly to moderately weathered dacite and weak zones. Rock masses were characterized in terms of rock mass rating (RMR), geological strength index (GSI) and Q System. Core samples were tested in the rock mechanics laboratory to determine uniaxial compressive strength, deformability parameters, unit weight, tensile strength and triaxial compressive strength properties. Finally, rock mass strengths were determined by empirical and numerical methods. Required support system was suggested.     

Tunnel reinforcement via topology optimization

Anchoring is a fundamental method for supporting tunnels. It reinforces the original rock mass and reduces the deformation along the tunnel surface. The topological complexity of its layouts renders a design methodology difficult. A numerical approach to reinforce the host ground becomes desirable. The present paper proposes a topology optimization method based on a two‐phase cell model and finite‐element discretization of the host ground. The element consists of the original rock and the reinforcing material. The design issue involves the distribution of the reinforcing materials. The relative ratios of the two phases in various elements will be optimized to reduce the compliance of the tunnel. The method enables the computer‐aided design for the support of underground structures. The capabilities of the method are demonstrated by the designs to support a deep tunnel under various in situ stresses. The results indicate that oriented reinforcement is needed along the direction of the largest absolute value of the principal stress. Copyright © 2000 John Wiley & Sons, Ltd.     

Upscaling of elastic properties for large scale geomechanical simulations

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.     

Engineering geological studies of the diversion tunnel, focusing on stabilization analysis and support design, Iran

In this paper a detailed engineering geological assessment of rock masses and support design studies at Garmi Chay dam site, has been carried out. This project is located in the northwest of Iran and will be used for flow control and water storage. The diversion tunnel of the dam has a diameter of 5.5 m and a length of 420 m and will be driven in slightly to highly weathered micaschist and trachy andesite rock units. The geological studies include field and laboratory investigations that based on the results; for more exact investigation, tunnel alignment was divided into three geotechnical zones. These zones consist mainly of highly weathered gray micaschists, dark red trachy andesites and slightly weathered gray micaschists, respectively. Then, for every zone, support capacity of rock masses was evaluated by means of empirical and numerical methods. The rock mass classification systems (RMR, Q, GSI, RSR, SRC and RMi), the convergence–confinement method and a 2D finite element computer software, Phase² were used for empirical and numerical method, respectively. According to the results acquired from these methods some stability problems were expected in the tunnel especially in highly weathered micaschist zone, so that in practice two big collapses occurred. Because of high weathering, low constants of rock masses and their soil-like behavior, the stability analysis by analytical method does not give illogical results in lightly weathered micaschist zone. The support system, suggested by empirical method, was applied and its performance was evaluated by means of numerical modeling. After installation the support suggested by Phase² program, the thickness of plastic zone and deformations around the tunnel decreased significantly. Consequently the agreement of these methods with each other was resulted and using combination of them was recommended for more reliable support design.     

Hybrid analytical and finite element numerical modeling of mass and heat transport in fractured rocks with matrix diffusion

Quantification of mass and heat transport in fractured porous rocks is important to areas such as contaminant transport, storage and release in fractured rock aquifers, the migration and sorption of radioactive nuclides from waste depositories, and the characterization of engineered heat exchangers in the context of enhanced geothermal systems. The large difference between flow and transport characteristics in fractures and in the surrounding matrix rock means models of such systems are forced to make a number of simplifications. Analytical approaches assume a homogeneous system, numerical approaches address the scale at which a process is operating, but may lose individual important processes due to averaging considerations. Numerical stability criteria limit the contrasts possible in defining material properties. Here, a hybrid analytical–numerical method for transport modeling in fractured media is presented. This method combines a numerical model for flow and transport in a heterogeneous fracture and an analytical solution for matrix diffusion. By linking the two types of model, the advantages of both methods can be combined. The methodology as well as the mathematical background are developed, verified for simple geometries, and applied to fractures representing experimental field conditions in the Grimsel rock laboratory.     

A coupled hydromechanical system defined through rock mass classification schemes

Changes in the hydraulic conductivity field, resulting from the redistribution of stresses in fractured rock masses, are difficult to characterize due to complex nature of the coupled hydromechanical processes. A methodology is developed to predict the distributed hydraulic conductivity field based on the original undisturbed parameters of hydraulic conductivity, Rock Mass Rating (RMR), Rock Quality Designation (RQD), and additionally the induced strains. The most obvious advantage of the methodology is that these required parameters are minimal and are readily available in practice. The incorporation of RMR and RQD, both of which have been applied to design in rock engineering for decades, enables the stress-dependent hydraulic conductivity field to be represented for a whole spectrum of rock masses. Knowledge of the RQD, together with the original hydraulic conductivity, is applied to determine the effective porosity for the fractured media. When RQD approaches zero, the rock mass is highly fractured, and fracture permeability will be relatively high. When RQD approaches 100, the degree of fracturing is minimal, and secondary porosity and secondary permeability will be low. These values bound the possible ranges in hydraulic behaviour of the secondary porosity within the system. RMR may also be applied to determine the scale effect of elastic modulus. As RMR approaches 100, the ‘softening’ effect of fractures is a minimum and results in the smallest strain-induced change in the hydraulic conductivity because the induced strain is uniformly distributed between fractures and matrix. When RMR approaches zero, the laboratory modulus must be reduced significantly in order to represent the rock mass. This results in the largest possible change in the hydraulic conductivity because the induced strain is applied entirely to the fracture system. These values of RMR bound the possible ranges in mechanical behaviour of the system. The mechanical system is coupled with the hydraulic system by two empirical parameters, RQD and RMR. The methodology has been applied to a circular underground excavation and to qualitatively explain the in situ experimental results of the macropermeability test in the drift at Stripa. Copyright © 1999 John Wiley & Sons, Ltd.     

Modelling the behaviour of a retaining wall in residual soils for a cut and cover construction of a deep station in Metro do Porto

This paper presents a re-appreciation of the ground characterisation and the criteria to select the most representative geomechanical parameters to consider in a numerical model to predict the behaviour of a retaining wall of a deep excavation in highly weathered granite rock masses and residual soils. This study was focused in the construction of a deep station of Metro do Porto, which involved a cut and cover solution, with unusual proportions (in plant and in depth), built in the typical Oporto's granite weathered profiles, being the excavation conducted with retaining walls consisting of multi-anchored concrete piles. Specific sections were carefully instrumented, due to the presence of historic buildings in the vicinity. The definition of representative model parameters was based on precise laboratory tests over high quality soil samples, including oedometer and high-precision triaxial tests. Geotechnical and geological characterisation of all the area for the original design, was initially based on in situ tests, such as SPT and rock masses classification, and on the local experience on this type of ground. After this construction, the assumptions of parameterisation, using a constitutive modelling based on new laboratory tests over high quality block samples, allowed a reanalysis of the assumptions on the design phase. A back-analysis of monitored displacements and forces during the construction was made, assuming the designed structural solutions, which were in fact implemented in construction, but considering the new approaches on the definition of the geomechanical parameters for the prevailing weathered rock masses, necessary for the numerical simulation based on the commercial software Plaxis®, using the Mohr-Coulomb and “Hardening-Soil” models. Some specific changes of the constructive sequence during the excavation and activation of supporting system were attained, by looking at the information found during the construction. The results of this parametrical re-approach and analysis of the singularities of highly weathered granite and corresponding residual soils masses for modelling of retaining walls of large excavations are discussed.     

泥质软岩地基水平抗力系数研究

通过室内水平静载试验和数值模拟对泥质软岩地基水平抗力系数的比例系数m进行研究,给出m值参考值,并讨论了影响地基岩土m值的因素,研究结果表明：试验所用的南宁盆地第三系泥岩m值取值范围为17～35 MN/m4,桩径大时取小值,埋深大时取大值;用ANSYS模拟水平静载试验,计算出的m值与室内试验值十分接近,可为今后数值计算m值提供参考。经分析m值的影响因素后认为适当增加桩的埋深、控制岩土含水率及提高混凝土的浇灌质量等可以提高m值     

Deformation and fracturing behaviour of discontinuous rock mass and damage mechanics theory

The mechanical behaviour of a rock mass is strongly affected by discontinuities such as faults and joints. In this paper, a damage mechanics theory is proposed which deals with some sets of discontinuities distributed in a rock mass, for example, joint systems. In this theory, the distributed discontinuities are characterized by a second-order symmetric tensor, called the damage tensor. By introducing the damage concept, the deformation and fracturing behaviour of the rock mass can be reated in a framework of continuum mechanics. A numerical procedure is developed in order to implement the damage mechanics model by using the finite element method. The theory and numerical analysis are applied to several laboratory tests and a practical underground opening problem. Numerical results are compared with measured data.     

Tunnel blasting techniques in difficult ground conditions

Summary The quality of tunnelling can be improved by proper blast design which takes into account the rock mass conditions. The effects of different rock mass properties on tunnel blast performance need to be assessed. The strength of the formation and joint orientation critically affected fragmentation and overbreak in a model study of blasting. Similar effects were noted in situ when the performance of a blast pattern in different rock mass conditions in the Tandsi inclines (Bihar, India) were analysed. Accordingly, the on-going blast pattern was modified for the poor ground conditions prevailing in the rest of the inclines. Improved fragmentation and smooth profile were obtained as a result; the rate of drivage improved considerably and the cost of excavation was reduced. Based on the observations in the model studies and the investigations at Tandsi, some guidelines for optimum blast design in difficult ground conditions are suggested.     

A Combined Three-Dimensional Geological-Geostatistical-Numerical Model of Underground Excavations in Rock

Summary. This paper exploits geological and borehole geotechnical data obtained in the exploratory phase of a tunneling project to investigate in a first place if the kriging interpolation scheme may effectively reproduce the spatial variability of rock mass quality (Rock Mass Rating, RMR) in the vicinity of tunnels. For this purpose a quick solver in Fortran has been developed that performs variography analysis of 3D spatial data, fast kriging estimations of RMR between borehole sampling locations at the centroids of the elements of the numerical model, and model validation. For the purpose of an integrated underground excavation design, a step further is made by incorporating into the 3D mechanical numerical model of the rock mass, the three-dimensional (3D) solid geological model, thus coupling the geology with the ground (geotechnical) model (i.e. each element of the numerical model is assigned a geological material). The mechanical properties of each finite difference cell (or Representative Elementary Volume) of the ground model were then prescribed according to its geological type, the spatial heterogeneity of the rock mass expressed quantitatively with the kriging model, and the upscaling calculations of the mechanical properties of the intact rocks determined in the laboratory, based on the size-effect (strength dependence on size) and Damage Theory. Furthermore, a preliminary numerical simulation of the advance of unsupported tunnels in the model of the heterogeneous rock mass was performed for illustration purposes.     

New techniques in rock mass classification: application to welded tuffs at the Nevada Yucca Mountain

Summary Many rock mass classification systems exist to assist the engineer in assessing the rock support requirements for underground design. On-going research in this area is directed at attempting to utilize the fractal dimension and the acoustic emission response of the tuffs at the Nevada Yucca Mountain to further aid in rock mass classification. Acoustic emission response is shown to be correlated with the porosity of the sample. Engineering behaviour of the rock varies dramatically with porosity; events and peak amplitude offer a means to distinguish between fracture porosity and pore porosity and consequently the engineering behaviour of the rock. Fractal dimension is used to characterize the roughness of fracture surfaces. Two fractal dimension calculation methods, one based on the semi-variogram for the surface and the other based on the use of dividers, are applied for this purpose. The divider method is shown to resolve deviation from a straight line; the semi-variogram method is shown to identify statistical similarity to various types of noise.Nomenclature D fractal dimension - AE acoustic emission - b b-value determined from log(frequency) against log(amplitude) plots - (h) semi-variogram function - h lag distance for semi-variogram function - H an exponent term related to fractal dimension asD=2 –H     

极震区岩体地震动力破坏若干问题探讨

汶川地震时极震区产生了严重的地震地质灾害,其中强烈地震动造成的岩体动力破坏是造成灾害的根本原因。从极震区含义、地震动特征、岩体地震动力破坏概念、地震松动岩体和方法论等方面初步探讨了极震区岩体地震动力破坏问题。极震区是未来地震的潜在震源区,区内的地震属于直下型。极震区地震动具有不同于非极震区的地震动特征,岩体地震动力破坏的复合性特点就是地震动的不确定性造成的。对极震区岩体动力破坏概念的理解应考虑地震动的特点。地震松动岩体是极震区地震动造成的一种特殊破坏类型,是形成震害次生灾害的重要原因。岩体工程地质力学等学科的思想方法和技术手段为研究极震区岩体地震动力破坏这一命题提供了良好的基础,预测和评价极震区因岩体动力破坏造成的工程震害和地质灾害,减轻和预防未来地震时的灾害损失,是极震区岩体地震动力破坏研究的目标和方向。     

“Tunnel Information and Analysis System”: A Geotechnical Database for Tunnels

Existing experience from the design and construction of underground works is of major importance for the improvement of the construction methods and procedures in tunnelling, especially under adverse and complex geological and geotechnical conditions. This experience can be of great value to geotechnical engineers and engineering geologists, if data acquired through the ground investigation, the design and the construction is systematically collected, categorized and stored in a properly structured database that enables a targeted access to it, as well as to proceed to correlations and analysis, based on engineering criteria. Such a database should be carefully designed to “connect” all available data through all the phases of a tunnel project and premises deep knowledge from the geological and geotechnical investigation to the final design and construction. In order to make substantial use of the experience accumulated from the construction of a great number of tunnels, a database named Tunnel Information and Analysis System (TIAS) was developed. The data source for TIAS database was 62 tunnels of the Egnatia Highway in northern Greece, many of which have been constructed under difficult geological conditions in weak rock masses. The data processed by TIAS came from a variety of sources such as geological mapping, boreholes, laboratory and in situ testing, geotechnical classifications, engineering geological behaviour, groundwater, design parameters, information concerning immediate support measures, construction records and cost. The purpose of the system, besides incorporating extended and multi-source data for easy access, is to provide a tool for turning data into usable information for the comparison of anticipated and encountered geological and geotechnical conditions, the evaluation of geotechnical classification and design methods and the relations regarding rock mass conditions and behaviour and immediate support methods and types.     

Discrete numerical modelling of rockfill dams

The aim of this study is to obtain quantitative information on the behaviour of rockfill used in embankment dams, and particularly on the influence of block breakage on the displacement field, from a numerical analysis using the Distinct element method. A methodology is set up to define the resistance of the 2D particles so that the same probability of breaking blocks may be reproduced as in a 3D material. The model uses the discrete element code PFC^2D (Itasca Consulting Group Inc., PFC2D (Particle Flow Code in Two Dimensions), Version 3.0, 2002) and considers breakable clusters of 2D balls. The different parameters are determined from experimental data obtained from laboratory tests performed on rock blocks. The model is validated by comparing the results of the simulation of shearing tests with actual triaxial tests on rockfill material published in the literature. The numerical analysis of block crushing in an actual dam is proposed in the last part of this paper. Copyright © 2006 John Wiley & Sons, Ltd.