• Geophysics and Geophysical Exploration
• /
• v.18 no.4
• /
• pp.223-231
• /
• 2015

Self-potential (SP) is sensitive to groundwater flow and there are many causes to generate SP. Among many mechanisms of SP, pore-fluid flow in porous media can generate potential without any external current source, which is referred to as electrokinetic potential or streaming potential. When calculating SP responses on the surface due to geothermal fluid within an engineered geothermal system (EGS) reservoir, SP anomaly is usually considered to be generated by fluid injection or production within the reservoir. However, SP anomaly can also result from geothermal water fluid within EGS reservoirs experiencing temperature changes between injection and production wells. For more precise simulation of SP responses, we developed an algorithm being able to take account of SP anomalies produced by not only water injection and production but also the fluid of geothermal water, based on three-dimensional finite-element-method employing tetrahedron elements; the developed algorithm can simulate electrical potential responses by both point source and volume source. After verifying the developed algorithm, we assumed a simple geothermal reservoir model and analyzed SP responses caused by geothermal water injection and production. We are going to further analyze SP responses for geothermal water in the presence of water production and injection, considering temperature distribution and geothermal water flow in the following research.

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

The oceanic tides have an effect on groundwater levels in coastal fractured rock aquifers. The observed groundwater table fluctuations caused by the effective stress through an aquifer are shown as sine curves similar with tidal fluctuation. To estimate a hydrogeologic parameter, tidal method is utilized with groundwater level fluctuations of two monitoring wells. Cross correlation function is used to calculate time lags between observed groundwater levels and tide, and the deeper well shows longer time lag. The storage coefficients calculated by using tidal efficiency and time lag show large differences. The storage coefficients obtained by using time lags are close to the result of slug test, and that of the deeper well shows closer value by slug test. The tidal efficiency is unsatisfied to apply in the tidal method because of an effect of phreatic aquifer and the vertical flow of groundwater through fractured confining bed. This tidal method can be an economical and effective way to define the parameter by considering the location of observation well and hydrogeologic characteristics of a coastal aquifer.

• Tunnel and Underground Space
• /
• v.8 no.4
• /
• pp.275-286
• /
• 1998

The effective permeability and the representative element volume(REV) of fracture network model were evaluated based on the parameters such as permeability tensor, principal permeability and the direction of principal permeability. The effective permeability ranges between the harmonic mean and the arithmetic mean of the local permeabilities of subdivided blocks. From the numerical experiments, which were for investigating the influence of model volume on the variation of flux for the cubic models, it was found that the variation of flux became reduced as the model volume approached REV. The variation of groundwater flux into the tunnel in the fracture network model was mainly dependent on the ratio of the tunnel length to model size rather than REV. And it was found that groundwater flux into the tunnel was not completely consistent between the fracture network model and the equivalent porous media model, especially when the ratio of the tunnel length to model size is small.

• The Journal of Engineering Geology
• /
• v.16 no.4 s.50
• /
• pp.411-427
• /
• 2006

This study aimed to recognize characteristics of groundwater flow in fractured bedrocks based on zonal pump-ing tests, slug tests, water quality logs and borehole TV camera logs conducted on two boreholes (NJ-11 and SJ-8) in the city of Naju. Especially, the zonal pumping tests using sin91e Packer were executed to reveal groundwater flow characteristics in the fractured bedrocks with depth. On borehole NJ-11, the zonal pumping tests resulted in a flow dimension of 1.6 with a packer depth of 56.9 meters. It also resulted in lower flow dimensions as moving to shallower packer depths, reaching a flow dimension of 1 at a 24 meter packer depth. This fact indicates that uniform permissive fractures take place in deeper zones at the borehole. On borehole SJ-8, a flow dimension of 1.7 was determined at the deepest packer level (50 m). Next, a dimension of 1.8 was obtained at 32 meters of packer depth, and lastly a dimension of 1.4 at 19 meters of packer depth. The variation of flow dimension with different packer depths is interpreted by the variability of permissive fractures with depth. Zonal pumping tests led to the utilization of the Moench (1984) dual-porosity model because hydraulic characteristics in the test holes were most suitable to the fractured bedrocks. Water quality logs displayed a tendency to increase geothermal temperature, to increase pH and to decrease dissolved oxygen. In addition, there was an increasing tendency towards electrical conductance and a decreasing tendency towards dissolved oxygen at most fracture zones.

• Journal of Korea Water Resources Association
• /
• v.32 no.3
• /
• pp.301-310
• /
• 1999

In this paper a numerical model and two analytical models in the hydraulically connected stream-aquifer system were analyzed to compare the induced infiltration rate curves derived from each model. And we also examined the effects of anisotropy of hydraulic conductivity and the direction of the ambient ground water flow on the quantification of the induced infiltration rate. The induced infiltration rate curve determined by models is very simple and useful for estimating the induced infiltration rate since it contains only four physical variables such as the induced infiltration rate, the pumping rate, the distance between the pumping well and the stream, and the ambient ground water flow rate. Under the conditions tested in this paper the induced infiltration rate curves resulted from the Wilson's analytical model and FEWA numerical model were in good agreement, and the anisotropic ratio of hydraulic conductivity was evaluated as a physical factor which influences the behaviour of the induced infiltration rate curve. The methods and results of the paper might Icad to improve the understanding of the induced infiltration phenomenon and can be applied to the planning and disign of pumping well and the optimal determination of the induced infiltration rate and pumping rate for water quality management of the water supply wells.

• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.1
• /
• pp.23-32
• /
• 1999

A new numerical method using solution-adaptive grids (SAG) is developed to solve the Richards' equation (RE) for unsaturated flow in porous media. Using a grid generation technique, the SAG method automatically redistributes a fixed number of grid points during the flow process, so that more grid points are clustered in regions of large solution gradients. The method uses the coordinate transformation technique to employ a new transformed RE, which is solved with the standard finite difference method. The movement of grid points is incorporated into the transformed RE, and therefore all computation is performed on fixed grid points of the transformed domain without using any interpolation techniques. Thus, numerical difficulties arising from the movement of the wetting front during the infiltration process have been substantially overcome by the new method. Numerical experiments for an one-dimensional infiltration problem are presented to compare the SAG method to the modified Picard method using a fixed grid. Results show that accuracy of a SAG solution using 41 nodes is comparable with the solution of the fixed grid method using 201 nodes, while it requires only 50% of the CPU time. The global mass balance and the convergence of SAG solutions are strongly affected by the time step size (Δt) and the weighting parameter (${\gamma}$) used for generating solution-adaptive grids. Thus, the method requires automated readjustment of Δt and ${\gamma}$ to yield mass-conservative and convergent solutions, although it may increase computational costs. The method can be effective especially for simulating unsaturated flow and other transport problems involving the propagation of a sharp-front.

• The Journal of Engineering Geology
• /
• v.18 no.2
• /
• pp.185-190
• /
• 2008

Nine wells have been developed for uses of thermal waters at the Icheon hot spa area. Drilling depths of those hot spring wells range from 166 to 294 m and their piezometric heads are located at about 50 m below the surface. Using the differences between the surface and bottom temperatures within all boreholes, we can simply estimate geothermal gradient in this area. Thus, we obtained the highest, lowest and average gradient values as $64^{\circ}C/km$ from SB-2 well, $45^{\circ}C/km$ from SB-1 well and approximately $54.28^{\circ}C/km$, respectively. However, observing the MRD-2 well additionally drilled into the depth of 996 m, we found out that this study area has widely experienced the temperature disturbance due to thermal groundwater penetration through the fracture systems within the depth of 720 m. Unlikely this phenomenon, we can conclude that the groundwater flow below the depth of 720 m does not exist. Therefore, using only those temperature data below the 720 m depth, we can estimate reasonable geo-thermal gradient values as $33^{\circ}C/km$ in this study area. Pumping test shows that outflowing temperature is $36^{\circ}C$ corresponding to the temperature logging data at 720 m depth.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.1918-1923
• /
• 2009

결정질 암반에서 지하수의 수리적 특성을 파악하기 위해 수행되고 있는 수리시험 방법은 정률법, 정압법, 순간주입(회복)법 등 세 가지로 구분할 수 있다. 본 연구에서 다루는 정압주입시험 (Constant Head Injection Test, CHIT)은 위의 정압법의 한 종류로 토목공학, 지질공학 분야에서 대상 구간의 투수계수 추정을 위해 널리 쓰이는 수리 시험이며, 이는 단일 패커나 이중패커를 이용하여 시험 구간을 격리하고, 격리된 구간에 일정한 압력으로 물을 주입하여 주입되는 물의 양을 파악함으로써, 시험 구간의 수리전도도(Hydraulic conductivity)를 산출하는 전통적인 수리시험이다. 본 연구에서는 수치실험을 통해 시험 구간 및 주입 압력의 크기 등 인위적인 요인에 의해 도출되는 투수계수가 어떻게 달라지는지에 대해 평가해 보았다. 일반적으로 단열 암반에서 수행한 정압주입시험의 해석에 있어 매질을 균질, 등방성 다공질이라는 가정으로 구간별 투수량계수를 산출하기 때문에, 다공성 매질의 지하수 유동을 모사하는 MODFLOW를 수치모사 코드로서 이용하였다. 시험구간의 크기 및 주입압력에 대한 민감도 분석 결과, 시험구간의 크기에 상관없이 수치모의에서 입력한 수리전도도 값에 비해 낮은 수리전도도 값이 산출되었으며, 주입 압력이 클수록 산출되는 수리전도도 값이 매질의 수리전도도 값과 차이가 났다. 민감도 분석 결과 현장수리시험에서 정압 주입시험에 의한 구간별 수리전도도 산출함에 있어 시험구간의 크기와 주입 압력 값에 대하여 고려해야 한다고 판단된다.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 1998.11a
• /
• pp.225-231
• /
• 1998

Open dump causes groundwater and soil contamination by leachate, air pollution by LFG (Landfill Gas). In this paper, in order to improve landfill researches which have been done about reduction of high leachate level and LFG collection in the Kimpo landfill separately, the effect of simultaneous flowing of leachate and LFG has been Studied. The HYLGS (Hanyang Leachate Gas Simulator) used in this study is a 3D, 2-phase, transient FDM model which can be applied to venting trenches in a landfill. From present numerical analysis it can be concluded that all the pressures of the Kimpo landfill grid system are almost the same and their maximum value in the center grid block of the system is approximately 26 m $H_2O$ (2.52 atm), that because the pressures of venting trench layer situated in the middle of the landfill have the lowest values and equal with air pressure, the venting trenches play an important role in landfill stabilization, that the flow of gas will be more difficult as time goes by owing to the increase of LGR(Leachate and gas ratio).

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.5 no.2
• /
• pp.113-122
• /
• 2007

In this study, uranium migration experiments have been performed using a natural groundwater and a granite core with natural fractures in a glove-box constructed to simulate an appropriate subsurface environment. Groundwater flow experiments using the non-sorbing anionic tracer Br were carried out to analyze the flow properties of groundwater through the fracture of the granite core. The result of the uranium migration experiment showed a breakthrough curve similar to that of the non-sorting Br. This result may imply that uranium migrates as anionic complexes through the rock fracture since uranium can form carbonate complexes at a given groundwater condition. The distribution coefficient $K_d$ of the uranium between the groundwater and the fracture filling material was obtained as low as 2.7 mL/g from a batch sorption experiment. This result agrees well with the result from the migration experiment, showing a faster elution of the uranium through the rock fracture. In order to analyze retardation properties of the uranium through the rock fracture, the retardation factor $R_d({\sim}16.2)$ was obtained by using the $K_d$ obtained from the batch sorption experiment and it was compared with the $R_d({\sim}14.3)$ obtained by using the result from the uranium migration experiment. The values obtained from the both experiments were very similar to each other. This reveals that the retardation of the uranium is mainly occurred by the fracture filling material when the uranium migrates through the fracture of a granite core.