• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2002.04a
• /
• pp.342-346
• /
• 2002

불포화 층상 해안 대수층 내에서의 밀도 의존적 지하수 유동 및 염분 이동에 대한 연구를 위해 하나의 지하수 유동-용질 이동 연동 수치 모델이 제시되었다. 이 수치 모델은 밀도 의존적 지하수 유동 지배 방정식, 염분 이동 지배 방정식 및 농도와 밀도의 관계식, 그리고 유한 요소법에 기초하여 개발되었다. 서로 다른 두가지 성질의 불포화 대수층이 고려되었다. 하나는 사질토층 위에 점토층이 존재하는 층상 대수층이고, 다른 하나는 사질토층과 점토층이 혼합된 두가지 물질로 구성된 균질화된 대수층이다. 수치모델의 결과는 층상 불균질성 (layered heterogeneity)가 해안 대수층 내에서의 밀도의존적 지하수 유동과 염분 이동에 있어서 매우 중요한 역할을 하고 있음을 보여준다. 그러한 층상 불균질성의 효과는 사질토층과 점토층과의 현저한 수리학적 및 수리역학적 성질의 차이에 기인한다 따라서 실제 해안 대수층 내에서 관찰되는 점토층을 적절히 고려하는 것이 보다 합리적고 타당한 해안 대수층내에서의 밀도 의존적 지하수 유동 및 염분 이동 해석을 가능하게 할 것이다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.62-62
• /
• 2023

강우가 지표면 아래로 침투할 때 초기에는 투수층이 불포화 상태이어서 부압이 작용하면서 침투할 것이다. Richards 식(Richards, 1931)을 써서 불포화 투수층의 침투를 모의할 수 있다. 강우가 지속되는 동안 하상 아래 어느 구간은 포화 상태가 되어 Richards 식을 더 이상 사용할 수 없다. 하지만 현재까지의 연구는 Richards 식을 사용하여 침투를 모의하는 오류를 범하고 있다. 강우에 의한 침투를 예측할 때 지표면에서의 침투율 q_b 가 필요한 데 현존하는 연구에서는 Horton 식(Horton, 1941)을 사용하여 초기 침투율 f_o 와 장시간 후 침투율 f_c 와 시간에 따라 지수함수로 감소하는 계수 k 의 3가지 계수값을 실험이나 현장 관측값에서 찾아서 쓰고 있다. 그런데, 이 계수값은 강우강도 r_i 가 클수록 침투율 q 가 커지는 물리 현상을 반영하지 못하는 한계가 있다. 본 연구에서 먼저 포화 투수층에서의 침투를 모의하는 식을 개발하였다. 지표면 아래에서 불포화 투수층에는 Richards식을 사용하고 포화 투수층에는 개발한 식을 사용하여 침투를 모의하였다. 또한 지표면에서의 침투율 q_b 를 구하는 공식을 개발하였다. 하상에서의 침투율의 최대값은 $q_{bmax}=-{\lambda}{\sqrt{2g(s-b)}}$ 일 것이다. 여기서 λ 는 투수층의 공극율, s 는 유출수면의 위치, b 는 지표면의 위치이다. 지표면에서의 침투율의 최소값 q_bmin 은 지표면 바로 아래 지점에서의 침투율일 것이다. 지표면에서의 침투율 q_b 로 q_bmax 와 q_bmin 사이의 적절한 값을 선택한다. 강우강도를 r_i 라고 하면 지표면 위 유출수의 연속방정식은 다음과 같다: $s-b={\int}(r_i-{\mid}q_b{\mid})dt$. 즉, 유출수면의 위치 s 는 강우강도 r_i 가 클수록 또는 지표면에서의 유출율의 크기 |q_b| 가 작을수록 크다. 또한 지표면에서의 침투율 q_b 와 지표면 아래에서의 침투율 q 는 s - b 가 클수록 크다. 따라서, 강우강도 r_i 가 클수록 침투율 q_b, q 가 큰 현상이 잘 반영되었다. 강우-침투-유출 모형실험을 수행하여 강우강도에 따라 침투율과 유출량이 다른 현상을 관측하여 수치실험 결과와 비교·검증하였다.

• The Journal of Engineering Geology
• /
• v.28 no.4
• /
• pp.727-740
• /
• 2018

The media in the undersaturated zone is defined as the uppermost layer of the water table at which the groundwater is unsaturated or saturated discontinuously. The properties of the unsaturated zone can affect the reduction of contaminants that flow from the lower part of soil to the water table. In recent, there have been problems in evaluating groundwater contaminations vulnerability because weighted value for permeability is given, regardless of anisotropy and heterogeneity in the unsaturated media. Geological media have various ranges of permeability. When applying the weighted value, representative of permeability for grain sizes standardized, to construction of contamination vulnerability, it will produce more exaggerated result than the case that considers unsaturated geological properties. In this study, we performed laboratory column tests considering two sets of the unsaturated layers in order to investigate the permeability in anisotropic unsaturated zone with anisotropy. On the basis of the tests, average permeability coefficients were calculated considering the properties of unsaturated media obtained from drill cores in the field. The final contamination vulnerability map constructed shows that the contamination vulnerability map applying the properties of geological media of the unsaturated zone coincides much better with the results measured in the field, compared to the case of contamination vulnerability considering the weighted value in the unsaturated zone.

• Journal of Soil and Groundwater Environment
• /
• v.10 no.5
• /
• pp.18-24
• /
• 2005

Water quality changes of wastewater effluent in the shallow aquifier condition was investigated using laboratory unsaturated and saturated SAT columns for over five months. Average DOC removal was 31.9% in the unsaturated SAT column whereas no removal occurred in the saturated SAT column. Under the shallow aquifer condition, nitrification was not completed in the unsaturated SAT column, releasing residual ammonium nitrogen into the saturated SAT column. Short retention time (one day) in the shallow unsaturated SAT column rendered DO of about 2 mg/L to the influent of the saturated SAT column. Phosphate was not removed at all in the unsaturated SAT column while complete removal was achieved in the saturated column. Consequently, organic and inorganic compounds were removed under the shallow aquifer condition as effectively as was in deep aquifer, except for the release of ammonium and relatively high DO into the saturated SAT column.

• The Journal of Engineering Geology
• /
• v.26 no.1
• /
• pp.87-99
• /
• 2016

We present the results of two-dimensional numerical simulations predicting the flow of groundwater in a fractured unsaturated zone. We applied the k-field distribution of permeability derived from discrete fracture network (DFN) modeling as the hydraulic properties of a model domain. To model an unsaturated zone, we set the depth from the ground surface to the underground aquifer. The rate of water infiltration into the unsaturated zone was divided into two parts, an artificial structure surface and unsaturated soil zone. The movement of groundwater through the unsaturated zone was simulated with particular emphasis on contaminant transport. It was clearly observed that the contaminants dissolved in groundwater transported vertically from the ground surface to the saturated zone.

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

Air Sparging (IAS, AS) is a ground-water remediation technique, in which organic contaminants are volatilized into air as they rise from saturated to vadose soil zone. This study was conducted to investigate the variation characteristics of TPH, VOCs and $CO_2$ for air sparging of diesel contaminated saturated soil. Initial TPH concentration was 10,000 mg/kg for saturated soil phase and 1,001 mg/L for soil aquifer phase. After 36 days of air sparging, the equilibrium temperature of 2-Dimension experiment system was $24.9{\pm}1.5^{\circ}C$. The saturated soil TPH concentration (in the C10 port close to air diffuser) was reduced to 66.0% of the initial value. The mass amount of $CO_2$ was 3,800 mg and 3,200 mg in air space (C70 port) and in unsaturated soil zone (C50 port), respectively. The VOCs production kinetic parameter was 0.164/day in the air space (C70 port) and 0.182/day in the unsaturated soils (C50 port).

• Journal of the Korean Society of Hazard Mitigation
• /
• v.9 no.3
• /
• pp.9-17
• /
• 2009

Most of the current drainage criteria have been developed on the basis of experimental field results and theoretical analyses of infiltration under saturated conditions. The objective of this study is to extend the understanding of pavement drainage systems by considering unsaturated condition in the sublayers. Analyses of unsaturated flows through pavements was performed by running finite element program(SEEP/W) with a range of pavement materials and drainage parameters. Meanwhile, the widely used DRIP program developed by FHWA is based on assumption of saturated condition of pavements. Differences between saturated and unsaturated condition in the sublayers of the pavements are verified. It is verified that for unsaturated conditions time to drain would take longer time compared to saturated condition.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.3C
• /
• pp.167-178
• /
• 2008

Rainfall-induced slope failures were simulated by seepage and stability analyses for actual slopes of weathered soils. After undisturbed sampling and testing on a specimen of unsaturated conditions, a seepage analysis was performed under actual rainfall and it was found that the pore water pressure increased at the boundary of soil and rock layers. The safety factor of slope stability decreased below 1.0 and the failure of actual slope could be simulated. Under design rainfall intensity, the seepage analysis could not include the effects of the antecedent rainfall and the rainfall duration. Due to these limitations, the safety factor of slope stability resulted in above 1.0, since the hydraulic head of soil layers had not be affected significantly. In the analysis of another slope failure, the parameters of unsaturated conditions were evaluated using artificial neural network (ANN). In the analysis of seepage, the boundary of soil and rock was saturated sufficiently and then the safety factor could be calculated below 1.0. It was found that the failure of actual slope can be simulated by ANN-based estimation.

• Economic and Environmental Geology
• /
• v.47 no.1
• /
• pp.49-59
• /
• 2014

Rainfall-induced landslides are caused by reduction of effective stress and shear strength due to rainfall infiltration. In order to analyze the susceptibility of landslides, the statistical analysis approach has been used widely but this approach has the limitation which cannot take into account of landslide triggering mechanism. Therefore, the physically based model which can consider the process of landslide occurrence was proposed and commonly used. However, the most previous physically based model analyses evaluate and consider the strength characteristics for saturated soil only in the susceptibility analysis. But the strength parameters for unsaturated soil such as matric suction should be considered with the strength parameters for saturated soil since the shear strength in unsaturated soil also plays important role in the stability of slope. Consequently this study suggested the modified physically based slope model which can evaluate strength characteristics for both of saturated and unsaturated soils. In addition, this study evaluated the thickness of saturated part in slope with rainfall intensity and hydraulic characteristics of slope on the basis of physically based model. In order to evaluate the feasibility, the proposed model was applied to practical example in Jinbu area, Gangwon-do, which was experienced large amount of landslides in July 2006. The ROC graph analysis was used to evaluate the validation of the model, and the analysis results were compared with the results of the previous analysis approach.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 1999.10a
• /
• pp.39-43
• /
• 1999

불포화대에서 LNAPL의 이동과 분포를 수치 모의를 통하여 예측하였다. 균질한 매질에서 LNAPL의 이동은 조립질 매질에서 빠르고, 세립질 매질에서 더 많은 면적으로 확산되며, 더 많은 LNAPL이 불포화대에 잔류한다. 조립질 매질내에 세립질층이 존재할 경우, 이 층이 지하수면으로부터 멀수록 LNAPL이 많이 포획된다. 조립질 매질에 세립질 또는 더 조립질인 매질이 렌즈 상으로 존재하는 환경에서는, LNAPL이 이들 렌즈를 통과하지 못한다. 불균질한 렌즈가 존재할 때의 LNAPL 분포를 초기조건으로 이용하여, 지하수 면의 수직 이동과 물의 침투에 따른 LNAPL의 이동을 모의하였다. 두 경우 모두 불포화대에 잔류되어 있던 LNAPL의 수직방향 이동이 증가되었다. 특히, 지하수면의 하강 시 LNAPL이 조립질 렌즈를 통해 이동하나, 세립질 렌즈를 통해서는 이동하지 못한다.