• 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.

• Journal of the Korean Society of Groundwater Environment
• /
• v.4 no.3
• /
• pp.111-115
• /
• 1997

Discrete fracture model has become one of the alternatives for the classical continuum model to simulate the irregular aspects of the fluid flow and the solute transport in fractured rocks. It is based on the assumptions that the discharge in a single fracture is proportional to the cube of the aperture and the fractured rock can be represented by the statistical assemblage of such single fractures. This study is intended to evaluate the effect of the fracture junction on the cubic law. Numerical solution of flow in junction system was obtained by using the Boundary-Fitted Coordinate System (BFCS) method. Results with different intersection angles in crossing fractures show that the geometry of the junction affects the discharge pattern under the same simulation conditions. Therefore, strict numerical and experimental examinations on this subject are required.

• Journal of the Korean GEO-environmental Society
• /
• v.2 no.1
• /
• pp.37-56
• /
• 2001

As a preliminary study to establish fracture network model in crystalline rocks, detail investigation on fracture characteristics were performed. Five fracture sets were determined on the basis of regional survey of geological structures and fractures on outcrops. Among the fracture sets, S1 set has the highest density and longest trace length of fractures which was identified on surface in the study area. S4 and S5 sets are composed of foliations and foliation parallel shear joints of gneisses, which are very important sets at the aspect of weighting of fracture length. For characterization of subsurface fractures, detail core logging was performed to identify fractures and fracture zones from five boreholes. Acoustic televiewer logging and borehole geophysical loggings produced images, orientations and geophysical properties of fractures which intersect with boreholes. According to the result of the investigations, subsurface fractures can be grouped as three preferred orientations(B1, B2 and B3), which correspond to S1, S2 and S4/S5 of surface fracture sets, respectively. Actually, B1 set is expected to be intensely developed at subsurface. However, it has low frequency of intersection with boreholes due to its parallel or sub-parallel direction to boreholes. According to the inference of conductive fractures, B1 and B3 sets have possibilities of water flow and their intersection lines are also thought to consist of important conduits of groundwater flow. In particular, faults which are parallel to foliations control major groundwater flow in the study area.

• Proceedings of the Korean Radioactive Waste Society Conference
• /
• 2005.06a
• /
• pp.231-238
• /
• 2005

On the basis of flow resistance theory the conceptual model and related mathematical descriptions is proposed for resistance modeling of groundwater flow in CPM(continuum Porous medium), DFN(discrete fracture network) and fractured-porous medium. The proposed model is developed on the basis of finite volume method assuming steady-state, constant density groundwater flow. The basic approach of the method is to evaluate inter-block flow resistance values for a staggered grid arrangement, i.e. fluxes are stored at cell walls and scalars at cell centers. The balance of forces, i.e. the Darcy law, is utilized for each control volume centered around the point where the velocity component is stored. The transmissivity (or permeability) at the interface is assumed to be the harmonic average of neighboring blocks. Flow resistance theory was utilized to relate the fluxes between the grid blocks with residual pressures. The flow within porous medium is described by three dimensional equations and that within an individual fracture is described by a two dimensional equivalent of the flow equations for a porous medium. Newly proposed models would contribute to develop flow simulation techniques with various matrix characteristics.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.4 no.1
• /
• pp.9-15
• /
• 2006

For the accurate safety assessment of potential radioactive waste disposal site which is located in the crystalline rock it is important to simulate the mass transportation through engineered and natural barrier system precisely, characterized by porous and fractured media respectively. In this work the methods to construct discrete fracture network for the analysis of flow and mass transport through fractured-porous medium are described. The probability density function is adopted in generating fracture properties for the realistic representation of real fractured rock. In order to investigate the intersection between a porous and a fractured medium described by a 2 dimensional rectangular and a cuboid grid respectively, an additional imaginary fracture is adopted at the face of a porous medium intersected by a fracture. In order to construct large scale flow paths an effective method to find interconnected fractures and algorithms of swift detecting connectivities between fractures or porous medium and fractures are proposed. These methods are expected to contribute to the development of numerical program for the simulation of radioactive nuclide transport through fractured-porous medium from radioactive waste disposal site.

• The Journal of Engineering Geology
• /
• v.10 no.1
• /
• pp.13-25
• /
• 2000

A three-dimensional discrete fracture network model based on probabilistic characteristics of fracture geometry and transmissivity was designed to calculate the conductivity tensor and to estimate theanisotropy of conductivity. The conductivities, $K_p$, obtained from the numerical simulation of single-holepacker test corresponded well to those from the field tests. From this, it can be concluded that thefracture network model designed in this study can represent hydraulic characteristics of in-situ fractured rock mass. Block-scale conductivities, $K_b$, estimated from the modelling of steady-state flow through the REV-scale block were ranged between the arithmetic mean and harmonic mean of theconductivity estimates from packer tests. The conductivity along north-south direction was 1.4 timesgreater than that along the east-west direction. It was concluded that the anisotropy of conductivitywas insignificant. It was also found that there was a little correlation between $K_b$ and $K_p$. This would be to that the conductivities from the packer test simulation was strongly dependent on thetransmissivity and the number of fractures within the packer test intervals.

• The Journal of Engineering Geology
• /
• v.31 no.1
• /
• pp.55-66
• /
• 2021

Hydraulic-mechanical (H-M) coupled numerical modeling was used to evaluate the evolution of hydrogeological properties in response to the installation and expansion of a borehole. A domain with a discrete fracture network was adopted for discontinuum modeling to simulate changes in fracture apertures. Comparison with real hydraulic test data shows that the effects of principal stress direction and expansion of borehole diameter were reasonably simulated by H-M coupled numerical modeling. The modeling confirmed that aperture changes depended on the principal stress direction, with an increase in aperture size due to vertical displacement being the dominant effect. A concentration of shear dilation around the borehole had an additional, subsidiary, effect on the hydrogeological evolution. These results show that the permeability of fractured rock can be increased by changing the hydraulic properties of a fracture through stress redistribution caused by the installation and expansion of a borehole.

• The Journal of Engineering Geology
• /
• v.6 no.1
• /
• pp.1-13
• /
• 1996

The fracture transmissivity($T_f$) is the most important parameter of fracture in assessing groundwater flow in fractured rock masses by using the DFN(Discrete Fracture Network) modeling. $T_f$, the most sensitive parameter m DFN modeling, is dependent upon aperture, size and filling characteristics of each fracture set. In the field test, the accuracy of $T_f$ can be increased with Borehole Acoustic Scanning (Televiewer) and Fixed Interval Length(FIL) test in constant head. $T_f$ values measured from FIL test was modified and estimated by each fracture set on the basis of the Cubic Law and the information of aperture and filling characteristics obtained from Televiewer. The modified $T_f$ results in the increase of confidence and reliability of modeling results including the amount of tunnel inflow.And, this approach would reduce the uncertaintity of the assessment for groundwater flow in fractured rock masses using the DFN modeling.

• Proceedings of the Korean Radioactive Waste Society Conference
• /
• 2006.11a
• /
• pp.230-231
• /
• 2006

• Proceedings of the Korean Radioactive Waste Society Conference
• /
• 2015.10a
• /
• pp.169-170
• /
• 2015