• Journal of Korea Water Resources Association
• /
• v.49 no.1
• /
• pp.61-72
• /
• 2016

In Korea, the aquifers at the coastal areas are mostly shallow alluvial unconfined aquifers. To simulate the flow and dispersion in unconfined aquifer, a FDM model has been developed to solve the nonlinear Boussinesq equation. Related analysis and verification have been executed. The iteration method is used to solve the nonlinearity, and the model shows 3-D shape because it is a 2-D y model that consider the undulation of water table and bottom. For the verification of the model, the output of flow module is compared to the 1-D analytic solution of Lee (1989) which have the drawdown or uplift boundary condition, and the two results show almost the same value. and the mass balance of dispersion module shows about 10% error. The developed model can be used for the analysis and design of the flow and dispersion in the unconfined aquifers. The model has been applied to the estuary area of Ssangcheon watershed, and the parameters have been deduced as a result : hydraulic conductivity is 90 m/day, and longitudinal dispersivity is 15 m. And the analysis with these parameters shows that the wells are situated in the influence circle of each others except for No. 7 well. Groundwater discharge to sea is $3700m^3/day$. And the chlorine ion ($cl^-$) concentration at the pumping wells increase at least 1000 mg/L if groundwater dam is not exist, so the groundwater dam plays an important role for the prevention of sea water intrusion.

• Journal of Soil and Groundwater Environment
• /
• v.11 no.6
• /
• pp.53-60
• /
• 2006

Single well injection withdrawal tracer tests with bromide were carried out at two wells developed in a horizontally heterogeneous fractured rock. The hydraulic conductivity of TW-1 well was 5 times larger than TW-2 well, and the average linear velocity of TW-2 well was 1.8 times faster than TW-1 well. The difference of hydrodynamic dispersions of two wells in the fractured rock was studied with the analysis of concentration breakthrough curves and cumulative mass recovery curves of bromide with withdrawal time, and the estimation of average travel distance, pore velocity, longitudinal dispersivity and longitudinal dispersion coefficient. The average travel distances of bromide were estimated to be 3.00 m in TW-1 well and 5.62 m in TW-2 well. The average pore velocities for the injection/withdrawal phase were estimated to be $4.31\;{\times}\;10^{-4}\;m/sec$ in TW-1 well and $8.08\;{\times}\;10^{-4}\;m/sec$ in TW-2 well. Average travel distance and pore velocity were higher in TW-2 well because of small effective porosity. Longitudinal dispersivities were estimated to be 28.73 cm in TW-1 well and 18.49 cm in TW-2 well, and bromide transport was 1.55 times faster in TW-1 well. Longitudinal dispersion coefficients were estimated to be $5.14\;{\times}\;10^{-6}\;m^2/sec$ in TW-1 well and $6.06\;{\times}\;10^{-6}\;m^2/sec$ in TW-2 well, and diffusion area was 1.18 times larger in TW-2 well.

• Journal of Soil and Groundwater Environment
• /
• v.11 no.6
• /
• pp.28-34
• /
• 2006

Laboratory column tests using $Cl^-$ tracer were conducted to study the correlation of soil particle distribution and hydrodynamic dispersion mechanism with three kinds of ununiformed soil samples, in which the ratio of gravel and sand versus silt and clay is 24.5 for S-1 soil, 4.48 for S-2 soil, and 0.4 for S-3 soil. Chloride breakthrough curves with time were fitted with gaussian functions. The relative concentrations of chloride were converged to 1.0 after 0.7 hours for S-1, 6.3 hours for S-2, and 389 hours for S-3. Average linear velocity, longitudinal dispersion coefficient, and longitudinal dispersivity were calculated by chloride breakthrough curves. Longitudinal dispersion coefficients were $1.20{\times}10^{-4}\;m^2/sec$ for S-1, $8.87{\times}10^{-7}\;m^2/sec$ for S-2, and $1.94{\times}10^{-9}\;m^2/sec$ for S-3. Peclet numbers calculated by the molecular diffusion coefficient of chloride and the mean grain diameters of soils were $2.59{\times}10^2$ for S-1, $6.27{\times}10^0$ for S-2, and $1.35{\times}10^{-4}$ for S-3. Mechanical dispersion was dominant for the hydrodynamic dispersion mechanism of S-1. Both mechanical dispersion and molecular diffusion were dominant for the hydrodynamic dispersion mechanism of S-2, but mechanical dispersion was ascendant over molecular diffusion. Hydrodynamic dispersion in S-3 was occurred mainly by molecular diffusion. When plotting three soils on the graph of $D_L/D_m$ versus Peclet number produced by Bijeljic and Blunt (2006), the values of $D_L/D_m$ for S-1 and S-2 were more than 2.0 order compared to their graph. S-3 was not plotted on their graph because the Peclet number was as small as $1.35{\times}10^{-4}$.

• Journal of Soil and Groundwater Environment
• /
• v.15 no.2
• /
• pp.1-9
• /
• 2010

Estimations of porosity and bulk density, particle size analysis of soli samples, tracer test and slug test were performed in a petroleum contaminated area of Kangwon for understanding characteristics of the aquifer. Porosities of the samples were estimated 0.158~0.257, and bulk densities were estimated as $1.73\sim2.10\;g/cm^3$. Majority proportion of the soil samples was 0.5~1.0 mm size. In the soil texture triangle, all samples were distributed at sand area. Uniformity coefficients were estimated as 7.71~10.39, and thus all samples were poorly-sorted. In the tracer test, Darcy velocity was estimated to $4.8\times10^{-6}$ cm/day, effective porosity was 0.175, and longitudinal dispersivity was 0.1 m. According to the slug test, hydraulic conductivities of the test wells were estimated as $2.243\times10^{-2}\sim1.634\times10^{-2}$ cm/sec. These hydrogeologic parameters can be used for efficient remediation design of the petroleum contaminated area.

• The Journal of Engineering Geology
• /
• v.18 no.4
• /
• pp.381-391
• /
• 2008

The eigenvalue technique is introduced to overcome the problem of truncation errors caused by temporal discretization of numerical analysis. The eigenvalue technique is different from simulation in that only the space is discretized. The spatially discretized equation is diagonized and the linear dynamic system is then decoupled. The time integration can be done independently and continuously for any nodal point at any time. The results of eigenvalue technique are compared with the exact solution and FEM numerical solution. The eigenvalue technique is more efficient than the FEM to the computation time and the computer storage in the same conditions. This technique is applied to the solute transport analysis in nonuniform flow fields around underground storage caverns. This method can be very useful for time consuming simulations. So, a sensitivity analysis is carried out by using this method to analyze the safety of caverns from nearly located contaminant sources. According to the simulations, the reaching time from source to the nearest cavern may takes 50 years with longitudinal dispersivity of 50 m and transversal dispersivity of 5 m, respectively.

• Journal of Soil and Groundwater Environment
• /
• v.15 no.6
• /
• pp.17-28
• /
• 2010

As a preliminary survey to improve efficiency of well-based permeable reactive barrier system for groundwater remediation, this site-scale study was carried to identify the flowpaths and controlling factors of plume at a remediation site in Suwon City, Korea. A total of 22 monitoring wells were installed as a grid system in the $4m{\times}4m$ square area by 1-m interval. For the groundwater characterization, various tests were performed including water-level monitoring, water sampling & analysis, pumping and slug tests, and tracer tests. The aquifer appeared to be unconfined with hydraulic conductivities (K) ranging from $2.6{\times}10^{-4}cm/s$ to $9.5{\times}10^{-3}cm/s$. The average linear velocity of groundwater was estimated to be $2.94{\times}10^{-6}m/s$, and the longitudinal dispersivity of a conservative tracer to be $5.94{\times}10^{-7}m^2/s$. Groundwater plume moves preferentially through the high-K zones, and the relatively high ion concentrations along the low-K zones implying deterred groundwater flow. Consequently, the spatial variation of hydraulic conductivity caused by aquifer heterogeneity and anisotropy appears to be the most important factor to maximize the effect of plume treatment system for application of in-situ groundwater remediation techniques.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.1899-1903
• /
• 2008

대수층의 저유량이 풍부한 강변여과수 개발 예정지역의 충적층(지표면하 $25{\sim}35\;m$ 구간)에서 수리전도도와 종분산지수의 규모종속효과를 규명하기 위해 양수시험과 수렴흐름 추적자시험이 수행 되었다. 양수시험과 추적자시험의 규모는 2 m 와 5 m 이었으며 양수시험은 5개 공, 추적자시험은 3개 공을 이용하여 수행되었다. 양수시험은 일정한 양수율($2,500\;m^3/day$)로 수행되었으며, 양수 시작 후 경과시간에 따른 수위변화 자료를 AQTESOLV 3.5 프로그램에 입력하여 해석하였다. 시험 대수층의 수리전도도는 양수정에서 $1.745{\times}10^{-3}\;m/sec$, 양수정에서 이격거리가 2 m 구간에서는 $2.161{\times}10^{-3}\;m/sec$와 $2.270{\times}10^{-3}\;m/sec$ 이었으며, 이격거리가 5 m 구간에서는 $2.452{\times}10^{-3}\;m/sec$와 $2.591{\time}10^{-3}m/sec$로 산정되었다. 그리고, 양수정에서 회복시험 시 Theis(Recovery) 방법에 의해 해석된 수리전도도는 $1.603{\times}10^{-3}\;m/sec$이었다. 양수정에서 관측정의 이격거리(d)에 따른 수리전도도(K) 증가함수는 log K=0.0693 log d-2.671와 log K=0.0817 log d-2.655로 추정되었으며, 결정 계수는 각각 0.965와 0.979로서 매우 높게 나타났다. 따라서 양수정에서의 이격거리가 멀수록 수리전도도가 증가하는 규모종속을 확인하였으며, 또한 시험대수층의 수리전도도가 방사상으로 유사하게 분포하고 있음을 알 수 있었다. 수렴흐름 추적자시험의 양수율은 $2,500\;m^3/day$ 이었으며, 2개의 주입정에 염소이온 5 kg을 순간 주입하였다. 염소이온의 농도이력곡선을 작성하여 초기도달시간과 최고농도의 차이를 분석하였으며, 누적질량회수곡선을 통해 양수 후 경과시간에 따른 염소이온의 질량회수율을 분석하였다. 그리고, 염소이온농도 대 누적질량회수율의 이력그래프를 작성하여 누적질량회수율에 따른 염소이온농도의 증가와 감소 변화를 분석하였다. 또한, 염소이온농도의 증가/감소 구간에 대한 선형회귀분석을 수행하여 농도 증가율과 감소율의 변화를 파악하였다. 양수정에서 관측된 경과시간별 염소이온농도 자료를 CATTI 코드의 "Converging Radial Flow With Instantaneous Injection" 해석법에 적용하여 종분산지수를 추정하였다. 양수정에서 이격거리가 2 m인 경우의 종분산지수는 0.4152 m, 이격거리가 5 m인 경우의 종분산지수는 3.2665 m이었다. 따라서 양수정에서 이격거리가 멀수록 종분산지수가 증가하는 규모종속효과를 확인하였으며, 또한 이격거리에 대한 종분산지수의 비는 각각 0.21과 0.65 정도로서 증가하였다.