• Tunnel and Underground Space
• /
• v.31 no.1
• /
• pp.66-81
• /
• 2021

The effects of directional fracture tensor components and first invariant of fracture tensor on deformation moduli and shear moduli of fractured rock masses is analyzed based on regression analysis performed between 3-D fracture tensor parameters and deformability of DFN blocks. Using one or two deterministic joint sets, a total of 224 3-D discrete fracture network (DFN) cube blocks were generated with various configurations of deterministic density and probabilistic size distribution. The fracture tensor parameters were calculated for each generated DFN systems. Also, deformability moduli with respect to three perpendicular direction of the DFN cube blocks were estimated based on distinct element method. The larger the first invariant of fracture tensor, the smaller the values for the deformability moduli of the DFN blocks. These deformability properties present an asymptotic pattern above the certain threshold. It is found that power-law function describes the relationship between the directional deformability moduli and the corresponding fracture tensor components estimated in same direction.

• Tunnel and Underground Space
• /
• v.30 no.1
• /
• pp.15-28
• /
• 2020

In this study, a procedure was introduced to estimate strength and deformation modulus of the 3-D discrete fracture network(DFN) systems using the distinct element method(DEM). Fracture entities were treated as non-persistent square planes in the DFN systems. Systematically generated fictitious fractures having similar mechanical characteristics of intact rock were combined with non-persistent real fractures to create polyhedral blocks in the analysis domain. Strength and deformation modulus for 10 m cube domain of various deterministic and stochastic 3-D DFN systems were estimated using the DEM to explore the applicability of suggested method and to examine the effect of fracture geometry on strength and deformability of DFN systems. The suggested procedures were found to effective in estimating anisotropic strength and deformability of the 3-D DFN systems.

• Tunnel and Underground Space
• /
• v.25 no.4
• /
• pp.341-358
• /
• 2015

This study analyzes the relationship between parameters of the discontinuity in Discrete Fracture Network model such as fracture intensity, fracture orientation, fracture size, fracture shape etc. In this paper, FracMan code was used to model and analyze 3D DFN. A sensitivity analysis was performed in order to analyze the relationship between linear fracture intensity measure ($P_{10}$) and parameters of the discontinuity in $100m{\times}100m{\times}100m$ model area. As a result the sensitivity analysis showed that key parameters affecting fracture intensity are fracture orientation (Trend / Plunge). Conversion factor($C_{13}$) for $P_{10}$, to calculate volumetric fracture intensity measure ($P_{32}$), is derived in case of vertical well and horizontal well when trend is $10^{\circ}$, $30^{\circ}$, $60^{\circ}$, $90^{\circ}$, $120^{\circ}$, $150^{\circ}$, $180^{\circ}$ (7cases) and plunge is $5^{\circ}$, $15^{\circ}$, $30^{\circ}$, $45^{\circ}$, $60^{\circ}$, $75^{\circ}$, $85^{\circ}$ (7cases). It is expected that this paper can be used effectively for modeling and understanding DFN model.

• The Journal of Engineering Geology
• /
• v.32 no.1
• /
• pp.41-57
• /
• 2022

A workflow is presented to estimate the size of a representative elementary volume and 3-D hydraulic conductivity tensor based on fluid flow analysis for a discrete fracture network (DFN). A case study is considered for a Cretaceous granitic rock mass at Gijang in Busan, Korea. The intensity and size of joints were calibrated using the first invariant of the fracture tensor for the 2-D DFN of the study area. Effective hydraulic apertures were obtained by analyzing the results of field packer tests. The representative elementary volume of the 2-D DFN was determined to be 20 m square by investigating the variations in the directional hydraulic conductivity for blocks of different sizes. The directional hydraulic conductivities calculated from the 2-D DFN exhibited strong anisotropy related to the hydraulic behavior of the study area. The 3-D hydraulic conductivity tensor for the fractured rock mass of the study area was estimated from the directional block conductivities of the 2-D DFN blocks generated for various directions in 3-D. The orientations of the principal components of the 3-D hydraulic conductivity tensor were found to be identical to those of delineated joint sets in the study area.

• Economic and Environmental Geology
• /
• v.53 no.3
• /
• pp.271-285
• /
• 2020

In this study, a numerical experiment related to the stress-strain analysis was performed on 3-D discrete fracture network(DFN) systems based on the distinct element method to evaluate the effect of joint geometry on deformability of jointed rock masses. Using one or two joint sets with deterministic orientation, a total of 12 3-D DFN blocks having 10m cube domain were generated with different joint density and size distribution. Directional deformation modulus of the DFN cube blocks were estimated along the axis directions of 3-D cartesian coordinate. In addition, deviatoric stress directions were chosen at every 30° of trend and plunge in 3-D for some DFN blocks to examine the variability of directional deformation modulus with respect to joint geometry. The directional deformation modulus of the DFN block were found to reduce with the increase of joint size distribution. The increase in joint density was less likely to have a significant effect on directional deformation modulus of the DFN block in case of the effect of rock bridges was relatively large because of short joint size distribution. It, however, was evaluated that the longer the joint size, the increase in the joint density had a more significant effect on the anisotropic deformation modulus of the DFN block. The variation of the anisotropic deformation modulus according to the variations in joint density and size distribution was highly dependent on the number of joint sets and their orientation in the DFN block. Finally, this study addressed a numerical procedure for stress-strain analysis of jointed rock masses considering joint geometry and discussed a methodology for practical application at the field scale.

• Tunnel and Underground Space
• /
• v.22 no.2
• /
• pp.93-105
• /
• 2012

This paper compares the two- and three-dimensional (2D and 3D) approaches for the numerical determination of the equivalent mechanical properties of fractured rock masses. Both orthogonally-fractured model and discrete fracture networks (DFN) were used for the geometry and 2D models were cut in various directions from 3D model to compare their mechanical properties. Geological data were loosely based on the data available from Sellafield, UK. Analytical method based on compliance tensor transformation was used for investigation in orthogonally fractured rock and numerical experiments were conducted on fractured rock mass with DFN geometry. It is shown that 2D approach always overestimates the elastic modulus of fractured rock masses by a factor of up to around two because fractures are assumed to be perpendicular to the model plane in 2D problems. Poisson ratios tend to have larger values in 2D analysis while there is opposite trend in some sections. The study quantitatively demonstrates the limitation of the 2D approach that uses the simplified model from true 3D geometry.

• The Journal of Engineering Geology
• /
• v.32 no.1
• /
• pp.59-72
• /
• 2022

This study employed a 3-D numerical analysis based on the distinct element method to estimate the strength and deformability of a Cretaceous biotite granitic rock mass at Gijang, Busan, Korea. A workflow was proposed to evaluate the scale effect and the representative elementary volume (REV) of mechanical properties for fractured rock masses. Directional strength and deformability parameters such as block strength, deformation modulus, shear modulus, and bulk modulus were estimated for a discrete fracture network (DFN) in a cubic block the size of the REV. The size of the mechanical REV for fractured rock masses in the study area was determined to be a 15 m cube. The mean block strength and mean deformation modulus of the DFN cube block were found to be 52.8% and 57.7% of the intact rock's strength and Young's modulus, respectively. A constitutive model was derived for the study area that describes the linear-elastic and orthotropic mechanical behavior of the rock mass. The model is expected to help evaluate the stability of tunnels and underground spaces through equivalent continuum analysis.

• Tunnel and Underground Space
• /
• v.13 no.6
• /
• pp.465-475
• /
• 2003

Groundwater inflow into the caverns constructed in fractured rock mass was simulated by numerical modeling, NAPSAC (DFN, discrete fracture network model) and NAMMU (CPM, continuous porous media model), a finite-element software package for groundwater flow in 3D fractured media developed by AEA Technology, UK. The input parameters for modeling were determined on surface fracture survey, core logging and single hole hydraulic test data. In order to predict the groundwater inflow more accurately, the anisotropic hydraulic conductivity was considered. The anisotropic hydraulic conductivities were calculated from the fracture network properties. With a minor adjustment during model calibration, the numerical modeling is able to reproduce reasonably groundwater inflows into cavern and the travel length and times to the ground surface along the flow paths in the normal, dry and rainy seasons.

• Geomechanics and Engineering
• /
• v.16 no.1
• /
• pp.21-34
• /
• 2018

Discontinuities considerably affect the mechanical and hydraulic properties of rock mass. These properties of the rock mass are influenced by the geometry of the discontinuities to a great extent. This paper aims to render an account of the geometrical parameters of several discontinuity sets related to the surrounding rock mass of Rudbar Lorestan Pumped Storage Power Plant powerhouse cavern making use of the linear and areal (circular and rectangular) sampling methods. Taking into consideration quite a large quantity of scanline and the window samplings used in this research, it was realized that the areal sampling methods are more time consuming and cost-effective than the linear methods. Having corrected the biases of the geometrical properties of the discontinuities, density (areal and volumetric) as well as the linear, areal and volumetric intensity accompanied by the other properties related to four sets of discontinuities were computed. There is an acceptable difference among the mean trace lengths measured using two linear and areal methods for the two joint sets. A 3D discrete fracture network generation code (3DFAM) has been developed to model the fracture network based on the mapped data. The code has been validated on the basis of numerous geometrical characteristics computed by use of the linear, areal sampling methods and volumetric method. Results of the linear sampling method have significant variations. So, the areal and volumetric methods are more efficient than the linear method and they are more appropriate for validation of 3D DFN (Discrete Fracture Network) codes.

• Tunnel and Underground Space
• /
• v.17 no.6
• /
• pp.438-447
• /
• 2007

This technical report is to introduce the research on SRM (Synthetic Rock Mass) which was presented in 2007 ISRM Congress at Lisbon by Prof, Fairhurst who speak with emphasis on its importance and potential in rock engineering. The Synthetic Rock Mass approach to jointed rock mass characterization (Pierce et al. 2007) is reviewed relative to existing empirical approaches and current understanding of jointed rock mass behaviour. The review illustrates how the key factors affecting the mechanical behaviour of jointed rock masses may be considered and demonstrates that the SRM approach constitutes a significant step forward in this field. This technique, based on two well-established methods, Bonded Particle Modelling in PFC-3D (Potyondy and Cundall, 2004) and Discrete Fracture Network simulation, employs a new sliding joint model that allows for large rock volumes containing thousands of pre-existing joints to be subjected to any non-trivial stress path. Output from SRM testing includes rock mass brittleness and strength, evolution of the full compliance matrix and primary fragmentation.