Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses

We present finite element simulations of temperature gradient driven rock alteration and mineralization in fluid saturated porous rock masses. In particular, we explore the significance of production/annihilation terms in the mass balance equations and the dependence of the spatial patterns of rock alteration upon the ratio of the roll over time of large scale convection cells to the relaxation time of the chemical reactions. Special concepts such as the gradient reaction criterion or rock alteration index (RAI) are discussed in light of the present, more general theory. In order to validate the finite element simulation, we derive an analytical solution for the rock alteration index of a benchmark problem on a two-dimensional rectangular domain. Since the geometry and boundary conditions of the benchmark problem can be easily and exactly modelled, the analytical solution is also useful for validating other numerical methods, such as the finite difference method and the boundary element method, when they are used to dear with this kind of problem. Finally, the potential of the theory is illustrated by means of finite element studies related to coupled flow problems in materially homogeneous and inhomogeneous porous rock masses. (C) 1998 Elsevier Science S.A. All rights reserved.

Geological and geotechnical surveys in tectonically active rock masses

The economic development of a region depends on the speed that people and goods can travel. The reduction of people and goods travel time can be achieved by planning smooth road layouts, which are obtained by crossing natural obstacles such as hills, by tunneling at great depths, and allowing the reduction of the road alignment length. The stress state in rock masses at such depths, either because of the overburden or due to the tectonic conditions of the rock mass induces high convergences of the tunnel walls. These high convergence values are incompatible with the supports structural performance installed in the excavation stabilization. In this article it is intended to evaluate and analyze some of the solutions already implemented in several similar geological and geotechnical situations, in order to establish a methodological principle for the design of the tunnels included in a highway section under construction in the region influenced by the Himalayas, in the state of Himachal Pradesh (India) and referenced by "four laning of Kiratpur to Ner Chowk section".

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.

Fracture imaging within a granitic rock aquifer using multiple-offset single-hole and cross-hole GPR reflection data

The sparsely spaced highly permeable fractures of the granitic rock aquifer at Stang-er-Brune (Brittany, France) form a well-connected fracture network of high permeability but unknown geometry. Previous work based on optical and acoustic logging together with single-hole and cross-hole flowmeter data acquired in 3 neighbouring boreholes (70-100 m deep) has identified the most important permeable fractures crossing the boreholes and their hydraulic connections. To constrain possible flow paths by estimating the geometries of known and previously unknown fractures, we have acquired, processed and interpreted multifold, single- and cross-hole GPR data using 100 and 250 MHz antennas. The GPR data processing scheme consisting of timezero corrections, scaling, bandpass filtering and F-X deconvolution, eigenvector filtering, muting, pre-stack Kirchhoff depth migration and stacking was used to differentiate fluid-filled fracture reflections from source generated noise. The final stacked and pre-stack depth-migrated GPR sections provide high-resolution images of individual fractures (dipping 30-90°) in the surroundings (2-20 m for the 100 MHz antennas; 2-12 m for the 250 MHz antennas) of each borehole in a 2D plane projection that are of superior quality to those obtained from single-offset sections. Most fractures previously identified from hydraulic testing can be correlated to reflections in the single-hole data. Several previously unknown major near vertical fractures have also been identified away from the boreholes.

Structural analysis of the vegetation on a highland granitic rock outcrop in Southeast Brazil

Granitic rock outcrops of the Brazilian southeast are either coastal or inland. The latter can often have high altitudes, such as in the summits of "Serra do Mar" and "Serra da Mantiqueira", where they are known as "Campos de Altitude". The landscape on these high altitude plateaux is often a mosaic of shrubs and treelets within a bunchgrass matrix, with sparse pteridophytes and other herbs, interspersed with variable extensions of rock outcrops. Despite the pervasiveness of rock outcrops in the Brazilian landscape, studies on the structural analysis of the vegetation on such formations are scarce. This study aimed to analyze the structure of the vegetation on a highland granitic rock outcrop in "Parque Estadual da Serra do Brigadeiro", state of Minas Gerais, Southeast Brazil (42°20' - 42°40' S and 20°20' - 21°00' W, 1,722 m of elevation ). Quantitative parameters of absolute and relative frequency and dominance (cover) were estimated. The group analysis used the Jaccard similarity index. Trilepis lhotzkiana, Panicum sp. 1, and Vellozia variegata presented the highest relative frequencies, relative dominances and importance values. These three species, along with Dyckia bracteata, Rhynchospora emaciate, and Tibouchina cf. manicata, represented 98.3% of the relative dominance. The remaining 1.7% referred to 22 remaining species. The distinction among quadrats within formed groups by cluster analysis was due to the great number of low frequency species.

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.

Numerical modelling of rock masses interacting with dams: the case of Ridracoli

A three-dimensional Direct Finite Element procedure is here presented which takes into account most of the factors affecting the interaction problem of the dam-water-foundation system, whilst keeping the computational cost at a reasonable level by introducing some simplified hypotheses. A truncated domain is defined, and the dynamic behaviour of the system is treated as a wave-scattering problem where the presence of the dam perturbs an original free-field system. The rock foundation truncated boundaries are enclosed by a set of free-field one-dimensional and two-dimensional systems which transmit the effective forces to the main model and apply adsorbing viscous boundaries to ensure radiation damping. The water domain is treated as an added mass moving with the dam. A strategy is proposed to keep the viscous dampers at the boundaries unloaded during the initial phases of analysis, when the static loads are initialised, and thus avoid spurious displacements. A focus is given to the nonlinear behaviour of the rock foundation, with concentrated plasticity along the natural discontinuities of the rock mass, immersed in an otherwise linear elastic medium with Rayleigh damping. The entire procedure is implemented in the commercial software Abaqus®, whose base code is enriched with specific user subroutines when needed. All the extra coding is attached to the Thesis and tested against analytical results and simple examples. Possible rock wedge instabilities induced by intense ground motion, which are not easily investigated within a comprehensive model of the dam-water-foundation system, are treated separately with a simplified decoupled dynamic approach derived from the classical Newmark method, integrated with FE calculation of dam thrust on the wedges during the earthquake. Both the described approaches are applied to the case study of the Ridracoli arch-gravity dam (Italy) in order to investigate its seismic response to the Maximum Credible Earthquake (MCE) in a full reservoir condition.

Geotecnia dos túneis do maciço de Paranhos: avaliação hidrogeomecânica exploratória

As classificações geomecânicas são uma das abordagens mais reconhecidas para estimar a qualidade do maciço rochoso, face à sua simplicidade e competência para gerir incertezas. As incertezas geológicas e geotécnicas podem ser avaliadas de forma eficaz usando classificações adequadas. Este estudo pretende enfatizar a importância das classificações geomecânicas e índices geomecânicos, tais como a Rock Mass Rating (RMR), a Rock Tunnelling Quality Index (Qsystem), o Geological Strength Index (GSI) e o Hydro‐Potential (HP) Value, para ajuizar a qualidade do maciço rochoso granítico das galerias subterrâneas de Paranhos (setor de Carvalhido ‐ Burgães; área urbana do Porto). Em particular, o valor hidro‐potencial (HP‐value) é uma classificação semi‐quantitativa aplicada a maciços rochosos que permite estimar as infiltrações de água subterrânea em escavações de terrenos rochosos. Para esta avaliação foi compilada e integrada uma extensa base de dados geológico‐geotécnica e geomecânica, apoiada na técnica de amostragem linear de superfícies expostas descontinuidades. Para refinar o zonamento geotécnico do maciço rochoso granítico, previamente realizado em 2010, foram coletadas amostras de rocha em pontos‐chave, com o objetivo de avaliar a sua resistência através do Ensaio de Carga Pontual (PLT). A aplicação das classificações geomecânicas foi realizada de uma forma equilibrada, estabelecendo diferentes cenários e tendo sempre em conta o conhecimento das características do maciço in situ. Apresenta‐se uma proposta de zonamento hidrogeomecânico com o objetivo de compreender melhor a circulação geo‐hidráulica do maciço rochoso granítico. Pretende‐se com esta metodologia contribuir para aprofundar o conhecimento do substrato rochoso do Porto, nomeadamente no que diz respeito ao seu comportamento geomecânico.

Evaluation of rock mechanical behaviors under uniaxial compression with reference to assessed weathering grades

A weathering classification for granitic rock materials from southeastern Brazil was framed based on core characteristics. The classification was substantiated by a detailed petrographic study. Indirect assessment of weathering grades by density, ultrasonic and Schmidt hammer index tests was performed. Rebound values due to Schmidt hammer multiple impacts at one representative point were more efficient in predicting weathering grades than averaged single impact rebound values, P-wave velocities and densities. Uniaxial compression tests revealed that a large range of uniaxial compressive strength (214-153 MPa) exists in Grade I category where weathering does not seem to have played any role. It was concluded that variability in occurrences of quartz intragranular cracks and in biotite percentage, distribution and orientation might have played a key role in accelerating or decelerating the failure processes of the Grade I specimens. Deterioration of uniaxial compressive strength and elastic modulus and increase in Poisson`s ratio with increasing weathering intensity could be attributed to alteration of minerals, disruption of rock skeleton and microcrack augmentation. A crude relation between failure modes and weathering grades also emerged.

Strain-dependent fluid flow defined through rock mass classification schemes

Strain-dependent hydraulic conductivities are uniquely defined by an environmental factor, representing applied normal and shear strains, combined with intrinsic material parameters representing mass and component deformation moduli, initial conductivities, and mass structure. The components representing mass moduli and structure are defined in terms of RQD (rock quality designation) and RMR (rock mass rating) to represent the response of a whole spectrum of rock masses, varying from highly fractured (crushed) rock to intact rock. These two empirical parameters determine the hydraulic response of a fractured medium to the induced-deformations The constitutive relations are verified against available published data and applied to study one-dimensional, strain-dependent fluid flow. Analytical results indicate that both normal and shear strains exert a significant influence on the processes of fluid flow and that the magnitude of this influence is regulated by the values of RQD and RMR.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.