• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.97-106
• /
• 1989

This paper is on the modeling of two-dimensional groundwater flow, which is the first step of the development of Dynamic System Model for groundwater flow and pollutant transport in subsurface porous media. The particular features of the model are its versatility and flexibility to deal with as many real-world problems as possible. Points as well as distributed sources/sinks are included to represent recharges/pumping and rainfall infiltrations. All sources/sinks can be transient or steady state. Prescribed hydraulic head on the Dirichlet boundaries and fluxes on Neumann or Cauchy boundaries can be time-dependent or constant. Sources/sinks strength over each element and node, hydraulic head at each Dirichlet boundary node and flux at each boundary segment can vary independently of each other. Either completely confined or completely unconfined aquifers, or partially confined and partially unconfined aquifers can be dealt with effectively. Discretization of a compound region with very irregular curved boundaries is made easy by including both quadrilateral and triangular elements in the formulation. Large-field problems can be solved efficiently by including a pointwise iterative solution strategy as an optional alternative to the direct elimination solution methed for the matrix equation approximating the partial differential equation of groundwater flow. The model also includes transient flow through confining leaky aquifers lying above and/or below the aquifer of interest. The model is verified against three simple cases to which analytical solutions are available. The groundwater flow model shall be combined with the model of pollutant transport in subsurface porous media. Then the combined model, with the applications of the Eigenvalue technique and the Dynamic system theory, shall be improved to the Dynamic System Model which can simulate the real groundwater flow and the pollutant transport accurately and effectively for the analyses and predictions.

• The Journal of Engineering Geology
• /
• v.4 no.2
• /
• pp.219-229
• /
• 1994

Since groundwater flow in fractured rocks is controlled by the distribution of fractures irregu1arly developed in space, it is not possible to understand the hydraulic characteristics of fractured aquifers using Theis equation which is applicable only to homogeneous isotropic confined aquifer. This study deals with the theoretical background of the fractal groundwater flow model with leakage, the methodology of calculation of the hydraulic parameters, and the application of the developed model to field data. From the result of the application of the fractal model to two field data in Hongcheon and Yusung areas, we obtained a good match between theoretical curves and observed curves, with the same hydraulic parameters at the pumping well and the observation well. In the two pumping test analyses, we have determined 1.9 of the fractal dimension. This means that the dimension of groundwater flow at these two sites is slightly smaller than radial flow.

• Journal of the Korean Society for Geothermal and Hydrothermal Energy
• /
• v.1 no.1
• /
• pp.1-8
• /
• 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 the effect of groundwater movement, understanding of thermohydraulic processes is necessary. The heat transfer phenomena for an aquifer heat storage are simulated by 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.001 are shaped circular, and the center is moved less than 5 m to the direction of groundwater flow 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 east 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.

• Water for future
• /
• v.27 no.1
• /
• pp.79-88
• /
• 1994

This study was conducted to analyse groundwater flow in the Bugok hot spring area using the MODFLOW model which can simulate three dimensional groundwater flow both in confined and unconfined aquifers. Based on this study the following conclusions were obtained: 1) The hydraulic conductivity and the specific storage of the aquifer were 0.0135 m/day and 0.020, respectively, and the model-predicted groundwater elevation agreed well with the observed one. 2) Simulation results showed that the groundwater level declines at the end of the one-year simulation period when the annual recharge rate is small and the annual pumping rate high, which is the worst combination. Except that combination, the groundwater level does not decline at the end of one-year simulation period indication the pumping rates used were allowable. 3) The safe yield depends upon the magnitudes of the recharge and pumping rates. The pumping rate should not produce excess decline of groundwater level around April when the water level is the lowest in a year.

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

Hydrologic responses to the 4September 2010 $M_W$ 7.1 and 22 February 2011 $M_W$ 6.2 Canterbury earthquakes ranged from near instantaneous co-seismic liquefaction and changes in groundwater levels, to more sustained (days to months) changes in river discharge, spring flow and groundwater level. There was some indication of a sustained change in aquifer properties. This paper presents some of the hydrographs from the September and February events, and compares the response to each event, briefly taking into account the location of the bore relative to each earthquake, together with other factors such as borehole depth. Over the months following the September earthquake, a pattern emerged of relatively short-term responses in the shallow aquifers and in the confined aquifer system, close to the coast. A longer term response appears to have occurred in inland, deep bores, where water levels 12 months after the September event were (in some cases) up to 20 metres higher than would have been expected based on simple modelling (see Figure 3). Some examples of these are highlighted.

• Korean Journal of Microbiology
• /
• v.35 no.1
• /
• pp.47-52
• /
• 1999

The study on the bacterial distribution was performed to evaluate the bacteriological safety in 4 deep confined aquifers drilled for drinking. All the investigated aquifers were oligotrophic, howevec bacterial densities and activities and the number of indicator bacteria in C-2 boreholes were much higher than those in other boreholes. Fecal pollution was observed in 3 boreholcs except C-1 borehole. The ratio of gram positive bacteria of the C-2 borehole was much lower than those 01 another boreholes. The numbers of total colifonns, fecal coliforms, and fecal Slreptococcus in the C-2 boreholes were 61, 33, 55 CFU (250 ml)$^{-1}$

• The Journal of Engineering Geology
• /
• v.34 no.2
• /
• pp.217-228
• /
• 2024

During deep drilling, confined aquifers can present various challenges such as the inability to remove cuttings, rapid groundwater influx, and mud loss. Particularly in flowing well conditions, it is essential to apply the suitable mud density since the aquifer can generates an overpressurized condition. This paper proposes a method for determining the suitable mud density while drilling (SMD) through confined aquifers using mud window analysis and applies it to a case study. The minimum mud density at each depth, which represents the lower limit of the mud window, is determined by the equivalent mud density pore pressure gradient (or by adding a trip margin) at that depth. The pore pressure gradient of a confined aquifer can be calculated using the piezometric level or well head pressure of the aquifer. As the borehole reaches the confined aquifer, there is a significant increase in pore pressure gradient, which gradually decreases with increasing depth. The SMD to prevent a kick can be determined as the maximum value among the minimum mud densities in the open hole section. After entering the confined aquifer, SMD is maintained as the minimum mud density at the top of the aquifer during the drilling of the open hole section. Additionally, appropriate casing installation can reduce the SMD, minimizing the risk of mud loss or invasion into the highly permeable aquifer.

• Economic and Environmental Geology
• /
• v.40 no.4
• /
• pp.403-418
• /
• 2007

The objective of this research was to characterize the fate and transport of Cr(VI) contaminated groundwater in the Daejeon industrial area. Five subsidiary monitoring wells were newly installed and two existing wells were utilized for the investigation and the reduction process of Cr(VI) contaminated groundwater of the Daejeon(Mun-pyeong) national groundwater monitoring station. The Cr(VI) concentrations at the shallow aquifer well of the station were in the range of 3.2-4.5 mg/L indicating continuous contamination. However, Cr was not detected at the deep bedrock well and the other monitoring wells except MPH-1 and 3. The Cr(VI) concentrations of MPH-1 and MPH-3 were below the drinking water guideline value (0.05 mg/L). Therefore, the plume of the Cr(VI) contaminated groundwater was predicted to be confined within the narrow boundary around the station. The soluble/exchangeable Cr(VI) concentrations were below the detection limit in all core and slime samples taken from the five newly installed wells. Although the exact source of contamination was not directly detected in the study area, the spatial Cr(VI) distribution in groundwater and characteristics of the core samples indicated that the source and the dispersion range were confined within the 100 m area from the monitoring station. The contamination might be induced from the unlined landfill of industrial wastes which was observed during the installation of an subsidiary monitoring well. For the evaluation of the natural attenuation of Cr(VI), available reduction capacities of Cr(VI) with an initial concentration of 5 mg/L were measured in soil and aquifer materials. Dark-gray clay layer samples have high capacities of Cr(VI) reduction ranging from 58 to 64%, which is obviously related to organic carbon contents of the samples. The analysis of reduction capacities implied that the soil and aquifer materials controlled the dispersion of Cr(VI) contamination in this area. However, some possibilities of dispersion by the preferential flow cannot be excluded due to the limited numbers of monitoring wells. We suggest the removal of Cr(VI) contaminated groundwater by periodical pumping, and the continuous groundwater quality monitoring for evaluation of the Cr(VI) dispersion should be followed in the study area.

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

Aquifer Thermal Energy Storage (ATES) system can be very cost-effective and renewable energy sources, depending on site-specific parameters and load characteristics. In order to develop an ATES system which has certain hydrogeological characteristics, understanding of the thermo hydraulic processes of an aquifer is necessary for a proper design of an aquifer heat storage system under given conditions. The thermo hydraulic transfer for heat storage is simulated using FEFLOW according to two sets of pumping and waste water reinjection scenarios of heat pump operation in a two layered confined aquifer. In the first set of model, the movement of the thermal front and groundwater level are simulated by changing the locations of injection and pumping well in seasonal cycle. However, in the second set of model the simulation is performed in the state of fixing the locations of pumping and injection well. After 365 days simulation period, the temperature distribution is dominated by injected water temperature and the distance from injection well. The small temperature change is appears on the surface compared to other slices of depth because the first layer has very low porosity and the transfer of thermal energy are sensitive at the porosity of each layer. The groundwater levels and temperature changes in injection and pumping wells are monitored to validate the effectiveness of the used heat pump operation method and the thermal interference between wells is analyzed.

• Journal of Soil and Groundwater Environment
• /
• v.21 no.4
• /
• pp.50-59
• /
• 2016

Spatial and temporal changes in nitrification activities and distribution of microbial population of ammonia oxidizing bacteria (AOB) and ammonia oxidizing archaea (AOA) in paddy soils were investigated. Soil samples were collected in March and October 2015 from rice paddy with and without the presence of confined animal feeding operations. Incubation experiments and quantitative polymerase chain reaction showed that AOA's contribution to nitrification kinetics was much higher in locations where organic nitrogen in animal waste is expected to significantly contribute to overall nitrogen budget, and temporal variations in nitrification kinetics were much smaller for AOA than AOB. These differences were interpreted to indicate that different microbial responses of two microbial populations to the types and concentrations of nitrogen substrates were the main determining factors of nitrification processes in the paddy soils. The copy numbers of ammonium monooxygenase gene showed that AOA colonized the paddy soils in higher numbers than AOB with stable distribution while AOB showed variation especially in March. Although small in numbers, AOB population turned out to exert more influence on nitrification potential than AOA, which was attributed to higher fluctuation in AOB cell numbers and nitrification reaction rate per cells.