• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.4
• /
• pp.103-112
• /
• 1991

In this study, a numerical model was established and applied to simulate the steady-state groundwater and heat flow in an isotropic, heterogeneous, three dimensional aquifer system with uniform thermal properties and no change of state. This model was developed as an aid in screening large groundwater-flow systems as prospects for underground waste storage. Driving forces on the system are external hydrologic conditions of recharge from precipitation and fixed hydraulic head boundaries. Heat flux includes geothermal heat-flow, conduction to the land surface, advection from recharge, and advection to or from fixed-head boundaries. The model uses an iterative procedure that alternately solves the groundwater-flow and heat-flow equations, updating advective flux after solution of the groundwater-flow equation, and updating hydraulic conductivity after solution of the heat-flow equation. Dierect solution is used for each equation. Travel time is determined by particle tracking through the modeled space. Velocities within blocks are linear interpolations of velocities at block faces. Applying this model to the groundwater-flow system located in Jigyung-ri. Songla-myun, Youngil-gun. Kyungsangbuk-do, the groundwater-flow system including distribution of head, temperature and travel time and flow line, is analyzed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.4
• /
• pp.23-32
• /
• 1988

The rates of infiltration contributed to the flow fo water in an unconfined aquifer under the partially penetrated stream at an ungaged station and the corresponding base flow in channel are coupled by using the hydraulic and/or hydrologic characteristics obtained from the geomorphologic and soil maps. For the determination of groundwater flow, the linearized model which is originally Boussinesq's nonlinear equation is applied in this study. Also, a stream flow routing model for base flow in channel is based on a simplification of the Saint-venant. The distributed runoff model with piecewise spatial uniformity is presented for obtaining its solution based on a finite difference technique of the kinematic wave equations. The method developed in this study was tested to the Bocheong watershed(area : $475.5km^2$) of the natural stream basin which is one of tributaries in Geum River basin in Korea. As a result, it is suggested that the rationality of hydro-graph separation according to a wide variability in hydrogeologic properties be worked out as developing the physically based subsurface model. The results of the present model are shown to be possible to simulate a base flow due to an arbitrary rate of infiltration for ungaged basins.

• Economic and Environmental Geology
• /
• v.38 no.2 s.171
• /
• pp.177-185
• /
• 2005

Data from 122 pumping tests were obtained from 100 boreholes in granites, volcanic rocks, metamorphic rocks, and Cretaceous and Tertiary sedimentary rocks, and then were analyzed using AQTESOLV. Results from 86 of the 122 tests ($71\%$) have an analytical solution corresponding to Theis (1935), Cooper-Jacob (1946), Papadopulos-Cooper (1967), Hantush (1962), Moench (1985), or Hantush-Jacob (1955), whereas the remaining 36 results ($29.5\%$) do not correspond to any of the analytical methods. Of the 86 results, only 17 match the Theis and Cooper-Jacob methods, indicating that the basic methods fer pumping test analysis are useful far only $14\%$ of the total data. This suggests that analytical solutions derived using leaky boundary conditions are appropriate for the analysis of pumping test data in fractured aquifers in this study. Furthermore, the results show the importance of carefully selecting an appropriate model for the analysis of pumping test data. Results from the 122 pumping tests were also analyzed using the GRF model. Using the Barker method, the results show that 77 of the 122 tests ($63\%$) have dimensions ranging between 1.1-2.9. Of these 77 solutions, ($39(44.2{\%})$) have a fractional dimension of 1.1-1.9, ($26(6.5{\%})$) show 2-dimensional radial flow also applicable to the Theis method, and ($38(49.3{\%})$) have dimensions of 2.1-2.9. The results show that groundwater flows according to a fractional flow dimension in fractured aquifers.

• Journal of Soil and Groundwater Environment
• /
• v.12 no.5
• /
• pp.19-32
• /
• 2007

In Korea, there have been various methods of estimating groundwater recharge which generally can be subdivided into three types: baseflow separation method by means of groundwater recession curve, water budget analysis based on lumped conceptual model in watershed, and water table fluctuation method (WTF) by using the data from groundwater monitoring wells. However, groundwater recharge rate shows the spatial-temporal variability due to climatic condition, land use and hydrogeological heterogeneity, so these methods have various limits to deal with these characteristics. To overcome these limitations, we present a new method of estimating recharge based on water balance components from the SWAT-MODFLOW which is an integrated surface-ground water model. Groundwater levels in the interest area close to the stream have dynamics similar to stream flow, whereas levels further upslope respond to precipitation with a delay. As these behaviours are related to the physical process of recharge, it is needed to account for the time delay in aquifer recharge once the water exits the soil profile to represent these features. In SWAT, a single linear reservoir storage module with an exponential decay weighting function is used to compute the recharge from soil to aquifer on a given day. However, this module has some limitations expressing recharge variation when the delay time is too long and transient recharge trend does not match to the groundwater table time series, the multi-reservoir storage routing module which represents more realistic time delay through vadose zone is newly suggested in this study. In this module, the parameter related to the delay time should be optimized by checking the correlation between simulated recharge and observed groundwater levels. The final step of this procedure is to compare simulated groundwater table with observed one as well as to compare simulated watershed runoff with observed one. This method is applied to Mihocheon watershed in Korea for the purpose of testing the procedure of proper estimation of spatio-temporal groundwater recharge distribution. As the newly suggested method of estimating recharge has the advantages of effectiveness of watershed model as well as the accuracy of WTF method, the estimated daily recharge rate would be an advanced quantity reflecting the heterogeneity of hydrogeology, climatic condition, land use as well as physical behaviour of water in soil layers and aquifers.

• The Journal of Engineering Geology
• /
• v.5 no.3
• /
• pp.289-300
• /
• 1995

At the initial stage of the underground reservoir design one should thoroughly consider surface and subsurface hydrology, hydrogeologic characteristics of aquifer system, and the function of cut - off wall because it is linked to the effective management. In this study, three dimensional finite difference model was applied to analyse the function of Ian underground reservoir at Kyungbuk Province. The steady and unsteady state conditions after construction of the underground dam were simulated through the model, and from these results the groundwater budget and the safe yield were determined. The model simulation indicates the infiltration of irrigation water to be one of the major factors of seasonal fluctuation of groundwater level. The recharge rates of irrigation water were estimated as 4.3mm/d during May and June, and 1.7mm/d during July and Agust. Groundwater recharge from the watershed area estimated to about $0.04m^3/s$, almost consistent through the year. In 1984, groundwater discharge through the transverse section of the dam was $0.002m^3/s$ and the optimum yield for two momths(July and Aguest)was $254000m^3$, however, the discharge became $0.013m^3/s$ in1993, implying the failure of cut -off function. without appropaiate of the cut - off wall, optiumum yield during the irrigaton period would be $93, 000m^3$.

• Journal of the Korean GEO-environmental Society
• /
• v.16 no.5
• /
• pp.35-42
• /
• 2015

This is conducted to observe underground temperature and to analyze its change affected by climate condition and soil infiltration in the mountainous area, Yesan region, Chungcheong-namdo province. Additionally, underground temperature change is also simulated using air temperature and soil thermal properties with a numerical model. Soil temperature monitoring data acquired from each depth, 20 cm, 50 cm, and 100 cm, indicates that the data within 50 cm in depth shows peak-shaped big fluctuation directly affected by air temperature and it at 100 cm has open-shaped small fluctuation. Underground temperature variation, a difference between high and low values, during monitoring period is weakly proportional to hydraulic conductivity of the sediment and it is assumed that water plays a part in delivering air temperature in soil. The underground temperature estimated by a numerical model is very similar to the observed data with an average value of 0.99 cross-correlation coefficient. From the result of this study, the aquifer unsaturated hydraulic conductivity of the soil and the groundwater recharge is likely to be able to estimate with underground temperature profile calculated using a numerical model.

• Journal of Soil and Groundwater Environment
• /
• v.21 no.6
• /
• pp.22-35
• /
• 2016

Global climate change could have an impact on hydrological process of a watershed and result in problems with future water supply by influencing the recharge process into the aquifer. This study aims to assess the change of groundwater recharge rate by climate change and to predict the sustainability of groundwater resource in Pyoseon watershed, Jeju Island. For the prediction, the groundwater recharge rate of the study area was estimated based on two future climate scenarios (RCP 4.5, RCP 8.5) by using the Soil Water Balance (SWB) computer code. The calculated groundwater recharge rate was used for groundwater flow simulation and the change of groundwater level according to the climate change was predicted using a numerical simulation program (FEFLOW 6.1). The average recharge rate from 2020 to 2100 was predicted to decrease by 10~12% compared to the current situation (1990~2015) while the evapotranspiration and the direct runoff rate would increase at both climate scenarios. The decrease in groundwater recharge rate due to the climate change results in the decline of groundwater level. In some monitoring wells, the predicted mean groundwater level at the year of the lowest water level was estimated to be lower by 60~70 m than the current situation. The model also predicted that temporal fluctuation of groundwater recharge, runoff and evapotranspiration would become more severe as a result of climate change, making the sustainable management of water resource more challenging in the future. Our study results demonstrate that the future availability of water resources highly depends on climate change. Thus, intensive studies on climate changes and water resources should be performed based on the sufficient data, advanced climate change scenarios, and improved modeling methodology.

• The Journal of Engineering Geology
• /
• v.27 no.3
• /
• pp.305-312
• /
• 2017

The most effective way to distinguish subsurface interfaces that produce various geophysical responses is through the integration of multiple geophysical methods, with each method detecting both a complementary and unique set of distinct physical properties relating to the subsurface. In this study, shallow seismic reflection (SSR) and ground penetrating radar (GPR) surveys were conducted at the Cheongju-Gadeok site of the Korea National Groundwater Monitoring Network to map the water table, which was measured at 12 m depth during the geophysical surveys. The water table proved to be a good target reflector in both datasets, as the abrupt transition from the overlying unsaturated weathered rock to the underlying saturated weathered rock yielded large acoustic impedance and dielectric constant contrasts. The two datasets were depth converted and integrated into a single section, with the SSR and GPR surveys conducted to ensure subsurface imaging at approximately the same wavelength. The GPR data provided detailed information on the upper ~15 m of the section, whereas the SSR data imaged structures at depths of 10-45 m. The integrated section thus captured the full depth coverage of the sandy clay, water table, weathered rock, soft rock, and hard rock structures, which correlated well with local drillcore and water table observations. Incorporation of these two geophysical datasets yielded a synthetic section that resembled a simplified aquifer model, with the best-fitting seismic velocity, dielectric constant, and porosity of the saturated weathered layer being $v_{seismic}=1000m/s$, ${\varepsilon}_r=16$, and ${\phi}=0.32$, respectively.

• Journal of Soil and Groundwater Environment
• /
• v.11 no.4
• /
• pp.10-17
• /
• 2006

Laboratory-scale two-dimensional aquifer physical model studies were conducted to assess the effect of air injection rate and air injection pattern on the removal of disel contaminated soil and groundwater by air/bio-sparging. The experimental results were represented that the optimal conditions in this experiment were as air injection rate of 1,000 ml/min and pulsed air injection pattern(15 min on/off). The results of the TPH reduction, DO consumption and $CO_2$ production indicate the effective biodegradation evidence of diesel. Based on our results, The minimal $O_2$ supply and pulsed air injection pattern could effectively enhance the diesel removal and the pulsing air injection had effect on oxygenation in this system. Thus, the cost of operating air/bio-sparging system will be reduced if optimal air injection rate and pulsed air injection pattern are applied to remediate contaminants.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2004.09a
• /
• pp.29-32
• /
• 2004

Gaseous partitioning tracer test has been used for determining the volume and spatial distribution of residual non-aqueous phase liquid (NAPL) in the unsaturated soils. In this study, an experimental method for measuring the content of gas-exposed NAPL as well as that of total NAPL in a sand during air sparging was developed. Two different gaseous phase NAPL-partitioning tracers were used; n-pentane, with very low water solubility, was used as the tracer that partitions into NAPL that is only in contact with the mobile gas, and chloroform, with fairly good water solubility, was selected for detecting total NAPL content in the sand. Helium and difluromethanewere used as the non- reactive tracer and water-partitioning tracers, respectively. Using n-decane as a model NAPL (NAPL saturation of 0.018), 25.6% of total NAPL was detected by n-pentane at the water saturation of 0.68. Oniy 9.1% of total NAPL was detected by n-pentane at the water saturation of 0.84. This result implies that the quantity of gas-exposed NAPL increased about three times when the water saturation decreased from 0.84 to 0.68. At the water saturation of 0.68, more than 90% of total NAPL was detected by chloroform while 65.8% of total NAPL was detected by chloroform at the water saturation of 0.84. Considering that the removal rate of NAPL during air sparging for NAPL-contaminated aquifer is expected to be greatly dependent upon the spatial arrangement of NAPL phase with respect to the mobile gas, this new approach may provide useful information for investigating the mass transfer process during air-driven remedial processes fer NAPL-contaminated subsurface environment.