• The Journal of Engineering Geology
• /
• v.3 no.2
• /
• pp.125-148
• /
• 1993

Groundwater flow in fractured rock masses is controlled by combined effects of fracture networks, state of geostafic stresses and crossflow between fractures and rock matrix. Furthermore the scaie dependent, anisotropic properties of hydraulic parameters results mainly from irregular paftems of fracture system, which can not be evaluated properly with the methods available at present. The basic assumpfion of discrete fracture network model is that groundwater flows only along discrete fractures and the flow paths in rock mass are determined by geometric paftems of interconnected fractures. The characteristics of fracture distribution in space and fracture hydraulic parameters are represented as the probability density functions by stochastic simulation. The discrete fracture network modelling was aftempted to characterize the groundwater flow in the vicinity of existing large cavems located in Wonjeong-ri, Poseung-myon, Pyeungtaek-kun. The fracture data of $1\textrm{km}^2$ area were analysed. The result indicates that the fracture sets evaluated from an equal area projection can be grouped into 6 sets and the fracture sizes are distributed in longnormal. The conductive fracture density of set 1 shows the highest density of 0.37. The groundwater inflow into a carvem was calculated as 29ton/day with the fracture transmissivity of $10^{-8}\textrm{m}^2/s$. When the fracture transmissivity increases in an order, the inflow amount estimated increases dramatically as much as fold, i.e 651 ton/day. One of the great advantages of this model is a forward modelling which can provide a thinking tool for site characterization and allow to handle the quantitative data as well as qualitative data.

• The Journal of Engineering Geology
• /
• v.16 no.1 s.47
• /
• pp.69-83
• /
• 2006

This study aims to evaluate a complex groundwater flow system around the underground oil storage caverns using the concept of hydraulic compartment. For the hydrogeological analysis, the hydraulic testing data, the evolution of groundwater levels in 28 surface monitoring boreholes and pressure variation of 95 horizontal and 63 vertical water curtain holes in the caverns were utilized. At the cavern level, the Hydraulic Conductor Domains(fracture zones) are characterized one local major fracture zone(NE-1)and two local fracture zones between the FZ-1 and FZ-2 fracture zones. The Hydraulic Rock Domain(rock mass) is divided into four compartments by the above local fracture zones. Two Hydraulic Rock Domains(A, B) around the FZ-2 zone have a relatively high initial groundwater pressures up to $15kg/cm^2$ and the differences between the upper and lower groundwater levels, measured from the monitoring holes equipped with double completion, are in the range of 10 and 40 m throughout the construction stage, indicating relatively good hydraulic connection between the near surface and bedrock groundwater systems. On the other hand, two Hydraulic Rock Domains(C, D) adjacent to the FZ-1, the groundwater levels in the upper and lower zones are shown a great difference in the maximum of 120 m and the high water levels in the upper groundwater system were not varied during the construction stage. This might be resulted from the very low hydraulic conductivity$(7.2X10^{-10}m/sec)$ in the zone, six times lower than that of Domain C, D. Groundwater recharge rates obtained from the numerical modeling are 2% of the annual mean precipitation(1,356mm/year) for 20 years.

• The Journal of Engineering Geology
• /
• v.26 no.3
• /
• pp.325-330
• /
• 2016

The second phase of low- and intermediate-level waste (LILW) disposal facility is under planned on the sedimentary rock in unsaturated zone. In this study, we created two meshes which were a matrix continuum mesh and a fracture continuum mesh to carry out 2 dimensional numerical modeling for groundwater flow in the unsaturated zone containing fractures focused on the second phase of LILW disposal facility. Two continuum meshes were developed using MINC in meshmaker module of TOUGH2 code. A fracture continuum mesh was included the k-field distribution of the permeability derived from the Discrete Fractured Network (DFN) modeling. To apply the unsaturated zone for the modeling, the gridding steps to generate mesh were developed. Each step to generate a mesh consisted of definition of materials, setting the initial conditions and creating grids using MINC. The methodology development of meshes in this study will be applied for more precise modeling of groundwater flow and mass transport.

• Economic and Environmental Geology
• /
• v.52 no.1
• /
• pp.81-90
• /
• 2019

In this study, the effects of fracture tensor component and first invariant on block hydraulic behaviors are evaluated in the 2-D DFN(discrete fracture network) systems. A series of regression analysis is performed between connected fracture tensor components and block hydraulic conductivities estimated at every $30^{\circ}$ hydraulic gradient directions for a total of 36 DFN systems having various joint density and size distribution. The directional block hydraulic conductivity seems to have strong relation with the fracture tensor component estimated in direction perpendicular to it. It is found that an equivalent continuum approach could be acceptable for the 2-D DFN systems under condition that the first invariant of fracture tensor is more than 2.0~2.5. The first invariant of fracture tensor seems highly correlated with average block hydraulic conductivity and can be used to evaluate hydraulic characteristics of the 2-D DFN systems. Also, a possibility of upscaling using the first invariant of fracture tensor for the DFN system is addressed through this study.

• Journal of the Korean Society of Groundwater Environment
• /
• v.5 no.4
• /
• pp.192-202
• /
• 1998

We investigated the geochemistry and environmental isotopes of granite-bedrock groundwater in the Yeongcheon diversion tunnel which is located about 300 m below the land surface. The hydrochemistry of groundwaters belongs to the Ca-HCO$_3$type, and is controlled by flow systems and water-rock interaction in the flow conduits (fractures). The deuterium and oxygen-18 data are clustered along the meteoric water line, indicating that the groundwater are commonly of meteoric water origin and are not affected by secondary isotope effects such as evaporation and isotope exchange. Tritium data show that the groundwaters were mostly recharged before pre-thermonuclear period and have been mixed with younger surface water flowing down rapidly into the tunnel along fractured zones. Based on the mass balance and reaction simulation approaches, using both the hydrochemistry of groundwater and the secondary mineralogy of fracture-filling materials, we have modeled the low-temperature hydrogeochemical evolution of groundwater in the area. The results of geochemical simulation show that the concentrations of Ca$\^$2＋/, Na$\^$＋/ and HCO$_3$and pH of waters increase progressively owing to the dissolution of reactive minerals in flow paths. The concentrations of Mg$\^$2＋/ and K$\^$＋/ frist increase with the dissolution, but later decrease when montmorillonite and illitic material are precipitated respectively. The continuous adding of reactive minerals, namely the progressively larger degrees of water/rock interaction, causes the formation of secondary minerals with the following sequence: first hematite, then gibbsite, then kaolinite, then montmorillonite, then illtic material, and finally microcline. During the simulation all the gibbsite is consumed, kaolinite precipitates and then the continuous reaction converts the kaolinite to montmorillonite and illitic material. The reaction simulation results agree well with the observed, water chemistry and secondary mineralogy, indicating the successful applicability of this simulation technique to delineate the complex hydrogeochemistry of bedrock groundwaters.

• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.29-46
• /
• 2004

It is important to identify geometries of fracture that act as a conduit of fluid flow for characterization of ground water flow in fractured rock. Fracture geometries control hydraulic conductivity and stream lines in a rock mass. However, we have difficulties to acquire whole geometric data of fractures in a field scale because of discontinuous distribution of outcrops and impossibility of continuous collecting of subsurface data. Therefore, it is needed to develop a method to describe whole feature of a target fracture geometry. This study suggests a new approach to develop a method to characterize on the whole feature of a target fracture geometry based on the Fourier transform. After sampling of specimens along a target fracture from borehole cores, effective frequencies among roughness components were selected by the Fourier transform on each specimen. Then, the selected effective frequencies were averaged on each frequency. Because the averaged spectrum includes all the frequency profiles of each specimen, it shows the representative components of the fracture roughness of the target fracture. The inverse Fourier transform is conducted to reconstruct an averaged whole roughness feature after low pass filtering. The reconstructed roughness feature also shows the representative roughness of the target subsurface fracture including the geometrical characteristics of each specimen. It also means that overall roughness feature by scaling up of a fracture. In order to identify the characteristics of permeability coefficients along the target fracture, fracture models were constructed based on the reconstructed roughness feature. The computation of permeability coefficient was performed by the homogenization analysis that can calculate accurate permeability coefficients with full consideration of fracture geometry. The results show a range between $10^{-4}{\;}and{\;}10^{-3}{\;}cm/sec$, indicating reasonable values of permeability coefficient along a large fracture. This approach will be effectively applied to the analysis of permeability characteristics along a large fracture as well as identification of the whole feature of a fracture in a field scale.