• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2003.03a
• /
• pp.113-120
• /
• 2003

The estimation of the coeff. of permeability is a major process that influenced in very large portion to the changing of gound water movement and estimation of the outflow’s quantity in the ground water movement modelling. To overcome the error caused by difficulty of the ground investigation for the hydraulic parameter’s estimation, in this paper, we study of correlation between electric resistivity and coeff. of permeability, and analyse of coeff. of permeability from electric resistivity. To accomplish this purpose, we carry out electrical logging, AMT survey as an in-situ test, and measure electric resistivity, porosity, permeability for the intact rock as a laboratory test. Based on the results of this research, it can be applied that estimation of the coeff. of permeability using electrical logging and AMT survey is a effective method for the impossible area of the ground investigation caused by geographical condition and petition.

• Journal of Soil and Groundwater Environment
• /
• v.14 no.1
• /
• pp.1-10
• /
• 2009

Groundwater flow due to intake of horizontal collector well in riverbank filtration site was analyzed by use of numerical groundwater modeling program (FEFLOW 5.1). Drawdowns of groundwater table nearby collector well were evaluated according to variations of several conditions; pumping rate, thickness of aquifer, offset distance from well to shore line of stream, conductance of streambed. It is observed that the drawdowns of groundwater table are clearly changed according to the variations of these conditions. The results of sensitive analysis shows that the thickness of alluvial aquifer and the offset distance are more sensitive than the conductance of streambed in evaluation of drawdown. This result implies that hydrogeological conditions, as like thickness of aquifer and its distribution in the site are important factors in site selection and evaluating the availability of riverbank filtration intake using horizontal collector well system. It is also revealed that numerical modeling using FEFLOW with 1-D discrete element feature can give efficient quantitative evaluation of horizontal collector well and estimation of availability of riverbank filtration site.

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

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

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.540-540
• /
• 2012

강우 발생시 유역에 집수된 물이 하천에 이르는 경로에 따른 유출은 지표유출, 지표하유출, 지하수유출로 구분된다. 정확한 수문순환 과정의 해석을 위해서는 지표 흐름뿐만 아니라 지표하 및 지하수 흐름의 해석이 중요한 실정이나 일반적으로 실무에서 사용되는 강우-유출해석 모형은 지표유출을 해석하기 위한 모형이 대부분이며, 지표하 유출과 침투량을 산정하는데 어려움이 있다. 일반적인 강우-유출해석 모형은 Horton 방법, NRCS 방법, Green-Ampt 방법에 의해 유효우량을 분리하며, 이 과정은 침투량을 직접적으로 모형화 할 수 없으므로 지표 및 지표하, 지하수 흐름을 복합적으로 해석할 수 있는 모형이 질적이나 양적으로 부족한 실정이다. 이러한 지표하 흐름과 침투량을 산정하기 위하여 FE-FLOW, PM, MS-VMS, GMS, GW-VISTAS, ARGUS 및 MODFLOW와 같은 지하수 모형을 사용하고 있다. 본 연구에서는 지표하유출 산정을 위한 침투량의 비교분석을 위해 현재 가장 범용되는 지하수 유동 모델링 프로그램인 Visual Modlfow 모형과 GMS 모형을 이용하여 침투량 산정을 위한 수치 모의를 진행하였다. 각 모형의 입력자료는 2009년 국립방재연구원에서 수행한 침투실험시설 자료를 이용하여 동일한 조건을 부여하고, 두 모형의 비교를 위해 Visual Modflow에서는 MODFLOW의 기본 해석방법인 유한차분법(FDM)을 이용하고, GMS 모형에서는 3차원 유한요소해석이 가능한 GMS-FEMWATER를 이용하였다. 두 모형의 수치모의 조건으로 2009년 국립방재연구원에서 수행한 침투실험방법과 동일하게 공극률에 따른 투수성 보도블럭의 구분과 50mm/hr, 100mm/hr, 150mm/hr, 200mm/hr의 강우강도별 선행함수조건에 따른 수치모의를 진행하였으며, 수치모의된 침투량의 적정성을 판단하기 위하여 국립방재연구원의 침투실험 결과자료와 비교분석하였다. 침투실험 자료와 각각 수치모의된 침투량을 비교분석한 결과 서로 유사한 경향을 보이고 있으나 초기 침투시 상대오차가 비교적 크게 발생하였다. 이는 수치모형의 경우 수리실험과는 다르게 모의시작과 동시에 해당 강우강도의 침투가 시작되므로 초기 유입 유출량 발생시간의 차이가 종료시간까지 누적 침투량에 미치는 것으로 판단되며, 매개변수에 많은 영향을 받는 것으로 판단된다.

• Journal of Soil and Groundwater Environment
• /
• v.9 no.4
• /
• pp.32-41
• /
• 2004

Analysis of water-level fluctuation due to goundwater-surface water interaction in coastal aquifers is carried out by numerical modeling. The conceptual model used in this study has a stream boundary and a tidal boundary that forms a right angle and the stream partially penetrates the aquifer. We analyzed the effect of each boundary and the simultaneous effect of the two boundary conditions. The area of influence caused by the stream boundary increased during the simulation, while the influence zone of the coastal boundary was relatively constant. The groundwater level near the zone where two boundaries meet may rise by the action of combined effect of the two boundaries or may not change by cancelling the effect of each boundary. Thereafter, care must be taken when hydraulic parameters are estimated using sinusoidal oscillations of hydraulic head in coastal aquifers. Sensitivity analysis is employed to develop insight into the controls on groundwater level fluctuations. In this study our analyses focused on the effect of conductance and the stream width to the aquifer nearby.

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

For various kinds of waters including surface water, shallow groundwater (＜70 m deep) and deep groundwater (500∼810 m deep) from the Pungki area, an integrated study based on hydrochemical, multivariate statistical, thermodynamic, environmental isotopic (tritium, oxygen-hydrogen, carbon and sulfur), and mass-balance approaches was attempted to elucidate the hydrogeochemical and hydrologic characteristics of the groundwater system in the gneiss area. Shallow groundwaters are typified as the 'Ca-HCO$_3$'type with higher concentrations of Ca, Mg, SO$_4$and NO$_3$, whereas deep groundwaters are the 'Na-HCO$_3$'type with elevated concentrations of Na, Ba, Li, H$_2$S, F and Cl and are supersaturated with respect to calcite. The waters in the area are largely classified into two groups: 1) surface waters and most of shallow groundwaters, and 2) deep groundwaters and one sample of shallow groundwater. Seasonal compositional variations are recognized for the former. Multivariate statistical analysis indicates that three factors may explain about 86% of the compositional variations observed in deep groundwaters. These are: 1) plagioclase dissolution and calcite precipitation, 2) sulfate reduction, and 3) acid hydrolysis of hydroxyl-bearing minerals(mainly mica). By combining with results of thermodynamic calculation, four appropriate models of water/ rock interaction, each showing the dissolution of plagioclase, kaolinite and micas and the precipitation of calcite, illite, laumontite, chlorite and smectite, are proposed by mass balance modelling in order to explain the water quality of deep groundwaters. Oxygen-hydrogen isotope data indicate that deep groundwaters were originated from a local meteoric water recharged from distant, topograpically high mountainous region and underwent larger degrees of water/rock interaction during the regional deep circulation, whereas the shallow groundwaters were recharged from nearby, topograpically low region. Tritium data show that the recharge time was the pre-thermonuclear age for deep groundwaters (＜0.2 TU) but the post-thermonuclear age for shallow groundwaters (5.66∼7.79 TU). The $\delta$$\^$34/S values of dissolved sulfate indicate that high amounts of dissolved H$_2$S (up to 3.9 mg/1), a characteristic of deep groundwaters in this area, might be derived from the reduction of sulfate. The $\delta$$\^$13/C values of dissolved carbonates are controlled by not only the dissolution of carbonate minerals by dissolved soil CO$_2$(for shallow groundwaters) but also the reprecipitation of calcite (for deep groundwaters). An integrated model of the origin, flow and chemical evolution for the groundwaters in this area is proposed in this study.

• The Journal of Engineering Geology
• /
• v.20 no.2
• /
• pp.169-181
• /
• 2010

This study aimed at simulating several responses to stresses caused by the ground water level variations around the K-1 oil stockpile. For this simulation, we considered the characteristic hydrogeological condition including the special occurrence of long and thick acidic dyke, which is regarded as the main geological structure dominating the ground water flow system at this study area. We activated twenty-four imaginary wells which are located in northern and southern area around central K-1 site. Each neighboring distance is altogether 300 m and whole distance between K-1 site and remote wells is 1,200 m. Through the modeling, we operated the long-term and continuous pumping tests and finally categorized five zones based on maximum pumping rates for the imaginary wells; zone I within 300 meter distance from K-1 site with a pumping rate of 50 $m^3/day$; zone II between 300 to 600 meter distance from K-1 site with a pumping rate of 75 $m^3/day$; zone III between 600 to 900 meter distance from K-1 site with 150 $m^3/day$; zone IV between 900 to 1,200 meter distance from K-1 site with 300 $m^3/day$; and zone V of acidic dyke area. At zone V, especially because of their possibility of high transmissivity for groundwater flow, it is necessary to control and restrict groundwater discharge.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.9 no.1
• /
• pp.13-21
• /
• 2011

To characterize the site specific properties of a study area for high-level radioactive waste disposal research in KAERI, the several geological investigations such as surface geological surveys and borehole drillings were carried out since 1997. Especially, KURT (KAERI Underground Research Tunnel) was constructed to understand the further study of geological environments in 2006. As a result, the first geological model of a study area was constructed by using the results of geological investigation. The objective of this research is to construct a hydrogeological model around KURT area on the basis of the geological model. Hydrogeological data which were obtained from in-situ hydraulic tests in the 9 boreholes were estimated to accomplish the objective. And, the hydrogeological properties of the 4 geological elements in the geological model, which were the subsurface weathering zone, the log angle fracture zone, the fracture zones and the bedrock were suggested. The hydrogeological model suggested in this study will be used as input parameters to carry out the groundwater flow modeling as a next step of the site characterization around KURT area.

• 한국신재생에너지학회:학술대회논문집
• /
• 2005.06a
• /
• pp.701-704
• /
• 2005

Aquifer Thermal Energy Storage (ATES) can be a cost-effective and renewable geothermal energy source, depending on site-specific and thermohydraulic conditions. To design an effective ATES system having influenced by groundwater movement, understanding of thermo hydraulic processes is necessary. The heat transfer phenomena for an aquifer heat storage are simulated using FEFLOW with the scenario of heat pump operation with pumping and waste water reinjection in a two layered confined aquifer model. Temperature distribution of the aquifer model is generated, and hydraulic heads and temperature variations are monitored at the both wells during 365 days. The average groundwater velocities are determined with two hydraulic gradient sets according to boundary conditions, and the effect of groundwater flow are shown at the generated thermal distributions of three different depth slices. The generated temperature contour lines at the hydraulic gradient of 0.00 1 are shaped circular, and the center is moved less than 5m to the groundwater flow direction in 365 days simulation period. However at the hydraulic gradient of 0.01, the contour center of the temperature are moved to the end of boundary at each slice and the largest movement is at bottom slice. By the analysis of thermal interference data between two wells the efficiency of the heat pump system model is validated, and the variation of heads is monitored at injection, pumping and no operation mode.