• Economic and Environmental Geology
• /
• v.52 no.4
• /
• pp.291-297
• /
• 2019

Fluid flow through rock fractures has been quantified using equations such as Stokes equations, Reynolds equation (or local cubic law), cubic law, etc. derived from the Navier-Stokes equations under the assumption that linear flow prevails. Therefore, these simplified equations are limited to linear flow regime, and cause errors in nonlinear flow regime. In this study, causal mechanism of nonlinear flow and critical Reynolds number were presented by carrying out fluid flow modeling with both the Navier-Stokes equations and the Stokes equations for a three-dimensional rough-walled rock fracture. This study showed that flow regimes changed from linear to nonlinear at the Reynolds number greater than 10. This is because the inertial forces, proportional to the square of the fluid velocity, increased enough to overwhelm the viscous forces. This tendency was also shown for the unmated (slightly sheared) rock fracture. It was found that nonlinear flow was caused by the rapid increase in the inertial forces with increasing fluid velocity, not by the growing eddies that have been ascribed to nonlinear flow.

• Economic and Environmental Geology
• /
• v.34 no.5
• /
• pp.485-498
• /
• 2001

This study was objected to identify the temporal and spatial variations oi groundwater quality and contamination using monthly groundwater monitoring data collected from 30 wells in Daejung watershed, Jeju Island. Water samples were analyzed for major cations. anions, and environmental isotopes including $^{18}0-H_2$O,$^{2}H-H_2$O and$^{15}$ N-NO$_3$The groundwater represented mostly Na(Mg)-HCO$_3$ type, with local change toward Ca-HCO$_3$ type and Na-Cl type. Groundwater quality depends upon various factors such as the local groundwater flowpaths, input of surface contaminants with recharge events, and sea-water intrusion along the coastal area. Nitrate contamination changed temporally according to recharge events and spatially. $\delta$$^{18}$ O-$\delta$$_2$H data for monthly sampled groundwaters showed distinctive clusters, implying that groundwater was originated from independent precipitation, and subsequently recharged very fash. Using$\delta$$^{15}$ N-NO$_3$ data, major sources of nitrate and its areas of influence could be identified. The areas under influence of livestock farms showed relatively high NO$_3$-N concentrations and$\delta$$^{15}$ N values higher than 5$\textperthousand$. The agricultural areas in southeastern part showed very high concentrations of NO$_3$-N with the $\delta$$^{15}$ N values of lower than 5$\textperthousand$.

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

• Tunnel and Underground Space
• /
• v.27 no.5
• /
• pp.303-311
• /
• 2017

The correlation analysis between the results obtained from DFN flow model and equivalent continuum flow model were conducted on total of 72 DFN blocks having various fracture geometry and domain size. A strong linear relation seems to exist between the two approaches under condition that normalized relative error for continuum behavior (ER) is less than 0.2, and the results from both methods are found to almost identical. To explore the field applicability of equivalent continuum flow model in DFN media, a total of 48 numerical schemes related to inflow of underground circular openings were implemented under various DFN configurations. The equivalent continuum flow model in DFN media with a constant hydraulic aperture was evaluated as valid. However, as the anisotropy increases due to variation of the hydraulic aperture, the results are likely to be overestimated compare to the DFN flow model.

• Journal of Korean Society of Environmental Engineers
• /
• v.31 no.6
• /
• pp.460-466
• /
• 2009

The mixing method of water and chemicals is significant in a small-scale water supply system because drinking water should be supplied with a certain quantity of remaining chemicals maintained. In the present study, the concentration distribution and the mixing index were obtained from four models, which were to find out the optimal mixing method of water and chemicals. The two models brought the good mixing effects out of the four, one for providing chemicals from the center of water supply pipe and the other for setting up the semicircle block at the downstream of the chemicals-providing pipe. As a result, the mixing effect was found out to be increased due to the diffusion and the disturbance of flows. In conclusion, these numerical results are expected to contribute to designing the optimal mixing system.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.129-137
• /
• 2000

This paper presents groundwater flow characteristics in discontinuous rock mass using fracture network program(NAPSAC) by statistical approach. Equivalent permeability coefficients are estimated from borehole data around Mabuk test tunnel site and fracture map on the arch of the tunnel. The reliability of fracture network model is obtained from determination of input data for statistical fracture network analysis from the real data(data of fracture network, data of hydraulic tests). The variation of permeability and mean anisotropic permeability coefficients are calculated from the realized model by increasing the size. As a result of analysis, a strong anisotropy of permeability is observed according to the direction of the fracture sets around the test tunnel.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.378-386
• /
• 2000

This paper presents groundwater flow characteristics in discontinuous rock mass using fracture network program(NAPSAC) by statistical approach. Equivalent permeability coefficients are estimated from borehole data around Mabuk test tunnel site and fracture map on the arch of the tunnel. The reliability of fracture network model is obtained from determination of input data for statistical fracture network analysis from the real data(data of fracture network, data of hydraulic tests). The variation of permeability and mean anisotropic permeability coefficients are calculated from the realized model by increasing the size. As a result of analysis, a strong anisotropy of permeability is observed according to the direction of the fracture sets around the test tunnel.

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

• Economic and Environmental Geology
• /
• v.34 no.3
• /
• pp.271-281
• /
• 2001

The characteristics of the heavy metal contamination in the soils affected by the tailings of the Palbong mine have been studied. The soils in the studied area consist mostly of loam by the particle size analysis, but a little of it, located near the stream, consist of loamy sand to sandy loam, finally to loam downward. The organic contents of soils are significantly low aoom 2 percent and the pH is in acidic ranging 6.0 $\pm$ 0.1. The samples of the parent rocks, the normal soils, the tailings and the channel deposits from the studied area were chemically analysed. From the result, the heavy metal concenlration of the soils is a little Jow compared with that of the parent rocks, shows the hydrologic process of the surface and the groundwater. The contamination of the tailings from the ore mining are high in lead, copper and arsenic. In the channel deposits the concenlrations of lead and copper are abnormally high but that of arsenic is uniquely low. And most of heavy metal contamination are decreased with the distance from the mine. It is caused by the properties of the surface and the ground water during the process of the heavy metal migration. The correlation-coefficient between sand and silt contents and the concentrations of Cd, Cu and Pb are significant but the amounts of As and Hg are increased with the clay contents. The dispersion of the heavy metals with the distance shows that the concentrations of them in the soils sampled at distance of 100 m to 200 m along the stream started near the Palbong mine are extremely high compared with those from other distances. These discrepancies are significant in Cd, Cu, Pb and Hg, but low in As. All the samples contain below detection limit of Cr+6 In the present stream water the concentrations of the heavy metals are not detected. So, it is interpreted that the concentrations in the soils are caused by the activities of the mining during the operation and have been continued by the dispersion from the tailings since after the closure of the mining, especially by the surface and ground water. The concentrations are diminished with the distance from the mining site, but in the interval of 800-2000 m increases abruptly. In the soil samples counted on the dispersion direction by wind, the lowering of the concentration is relatively uniform with the distance from the mining site. So, the rapid increase of the heavy metal concentrations is presumed to have been caused by the ground-water movement. In the migration of the heavy metals, the groundwater conditions, such as pH, Eh, the contents of colloidal particles, and Mn and Fe oxides are closely involveo. Integrating with these factors, it is interpreted that the groundwater conditions which have caused the heavy metal contamination of the studied area are those that the pH is about 3 in oxidized conditions, the contents of the colloidal particles are low, and Mn and Fe oxides are not involved in the migration of the heavy metals. Meanwhile, the vegetables growing on the soils in the studied area are not affected by the contamination of the heavy metals.

• Journal of the Korean Geotechnical Society
• /
• v.28 no.4
• /
• pp.101-113
• /
• 2012

The degradation of a tunnel drainage system leads to increases in pore water pressure around the tunnel and the lining stress, which results in affecting the tunnel stability. In the present study of the Namsan 3th tunnel, more than 30 year old tunnel, the effects of the drainage performance reduction due to drain hole clogging on the tunnel lining stability were investigated by examining pore water pressure distribution around the tunnel and the lining stresses through numerical analysis. Groundwater flow modeling on the Mt. Namsan region was done first and 3D seepage and coupled stress-pore water pressure finite element analysis were performed on the tunnel using the results of the groundwater flow modeling. The pore water pressure distribution and the tunnel lining stresses could be predicted using a drain hole outflow data measured in the tunnel site. This analysis method may be used to evaluate the current stability of old tunnels for which in most cases field investigations and related information are not readily available.