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

강우 발생시 유역에 집수된 물이 하천에 이르는 경로에 따른 유출은 지표유출, 지표하유출, 지하수유출로 구분된다. 정확한 수문순환 과정의 해석을 위해서는 지표 흐름뿐만 아니라 지표하 및 지하수 흐름의 해석이 중요한 실정이나 일반적으로 실무에서 사용되는 강우-유출해석 모형은 지표유출을 해석하기 위한 모형이 대부분이며, 지표하 유출과 침투량을 산정하는데 어려움이 있다. 일반적인 강우-유출해석 모형은 Horton 방법, NRCS 방법, Green-Ampt 방법에 의해 유효우량을 분리하며, 이 과정은 침투량을 직접적으로 모형화 할 수 없으므로 지표 및 지표하, 지하수 흐름을 복합적으로 해석할 수 있는 모형이 질적이나 양적으로 부족한 실정이다. 이러한 지표하 흐름과 침투량을 산정하기 위하여 FE-FLOW, PM, MS-VMS, GMS, GW-VISTAS, ARGUS 및 MODFLOW와 같은 지하수 모형을 사용하고 있다. 본 연구에서는 지표하유출 산정을 위한 침투량의 비교분석을 위해 현재 가장 범용되는 지하수 유동 모델링 프로그램인 Visual Modlfow 모형과 GMS 모형을 이용하여 침투량 산정을 위한 수치 모의를 진행하였다. 각 모형의 입력자료는 2009년 국립방재연구원에서 수행한 침투실험시설 자료를 이용하여 동일한 조건을 부여하고, 두 모형의 비교를 위해 Visual Modflow에서는 MODFLOW의 기본 해석방법인 유한차분법(FDM)을 이용하고, GMS 모형에서는 3차원 유한요소해석이 가능한 GMS-FEMWATER를 이용하였다. 두 모형의 수치모의 조건으로 2009년 국립방재연구원에서 수행한 침투실험방법과 동일하게 공극률에 따른 투수성 보도블럭의 구분과 50mm/hr, 100mm/hr, 150mm/hr, 200mm/hr의 강우강도별 선행함수조건에 따른 수치모의를 진행하였으며, 수치모의된 침투량의 적정성을 판단하기 위하여 국립방재연구원의 침투실험 결과자료와 비교분석하였다. 침투실험 자료와 각각 수치모의된 침투량을 비교분석한 결과 서로 유사한 경향을 보이고 있으나 초기 침투시 상대오차가 비교적 크게 발생하였다. 이는 수치모형의 경우 수리실험과는 다르게 모의시작과 동시에 해당 강우강도의 침투가 시작되므로 초기 유입 유출량 발생시간의 차이가 종료시간까지 누적 침투량에 미치는 것으로 판단되며, 매개변수에 많은 영향을 받는 것으로 판단된다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.39 no.5
• /
• pp.579-586
• /
• 2019

Pilot Point Method (PPM) is one of the popular methods to search hydraulic conductivities in the inverse method using groundwater flow equations. In this study, the Simultaneous Search based Pilot Point Method (SSBM) was applied with diverse information (e.g. hydraulic heads and/or tracer concentration) applications over previously developed sensitivity based Pilot Point Method (e.g. D-optimality based Pilot Point Method: DBM). In the case of DBM, due to the minimized the variance size, tracer concentration can be recognized as a tool to control the searching space of hydraulic conductivities. SSBM reduced the procedure of hydraulic conductivity searching, though it produced more variance for exploring hydraulic conductivities. In addition, SSBM was dependent on the initial hydraulic conductivity values for search finalized hydraulic conductivities. When tracer concentration was applied, searching hydraulic conductivities was more preferable than only when hydraulic head was applied. Applications of various data for searching hydraulic conductivities is recommended as a more efficient way.

• Ocean and Polar Research
• /
• v.41 no.2
• /
• pp.63-77
• /
• 2019

In order to investigate the behavior and seasonal variability of Mn as one of the bio-essential metals in the Amundsen sea, which is known as the most biologically productive coastal area around the Antartica, seawaters were collected using a clean sampling system for 10 stations (96 ea) in 2014 (ANA04B) and for 12 stations (139 ea) in 2016 (ANA06B) surveys of RV ARAON. Dissolved and particulate Mn concentration varied in the range of 0.15-4.43 nmol/kg and <0.01 to 2.42 nM in 2014 and in the range of 0.25-4.15 nmol/kg and 0.01-2.64 nM in 2016, respectively. From the sectional distribution of dissolved and particulate Mn, it might be suggested that dissolved/particulate Mn was provided from iceberg melting and diffusion/resuspension from sediments, respectively. Although this sea is highly productive, there was little evidence regarding the biological origin of dissolved Mn, but particulate Mn only in sea ice and offshore areas could be explained as originating from organic matters, e.g. phytoplanktons. And it could be suggested that the subsurface maximum of dissolved Mn was formed by isopycnal transport of melting materials from ice wall to offshore. Compared to early (2014) summer, temperature, salinity, biomass, dissolved and particulate Mn in late (2016) summer indicated that temporal variations might be resulted from the reduction of ice melting and mCDW flow, which induced a reduction in resuspension. In addition, in the late summer, particles including biomass were reduced, which brought about a reduction in the removal rate of dissolved Mn.

• The Journal of Engineering Geology
• /
• v.22 no.3
• /
• pp.343-351
• /
• 2012

This study proposes a modified equation to calculate the factor of safety for an infinite slope considering the saturation depth ratio as a new variable calculated from rainfall infiltration into unsaturated soil. For the proposed equation, this study introduces the concepts of the saturation depth ratio and subsurface flow depth. Analysis of the factor of safety for an infinite slope is conducted by the sequential calculation of the effective upslope contributing area, subsurface flow depth, and the saturation depth ratio based on quasi-dynamic wetness index theory. The calculation process makes it possible to understand changes in the factor of safety and the infiltration behavior of individual rainfall events. This study analyzes stability changes in an infinite slope, considering the saturation depth ratio of soil, based on the proposed equation and the results of soil column tests performed by Park et al. (2011 a). The analysis results show that changes in the factor of safety are dependent on the saturation depth ratio, which reflects the rainfall infiltration into unsaturated weathered gneiss soil. Under continuous rainfall with intensities of 20 and 50 mm/h, the time taken for the factor of safety to decrease to less than 1.3 was 2.86-5.38 hours and 1.34-2.92 hours, respectively; in the case of repeated rainfall events, the time taken was between 3.27 and 5.61 hours. The results demonstrate that it is possible to understand changes in the factor of safety for an infinite slope dependent on the saturation depth ratio.

• Journal of the Korean Society of Groundwater Environment
• /
• v.5 no.2
• /
• pp.88-100
• /
• 1998

Since fractures may serve as major conduits of groundwater flow in crystalline rocks, characterization of conductive fractures is especially important for interpretation of flow system. In this study, characterization of fractures to investigate hydraulically conductive fractures in gneisses at an abandoned mine area was performed. The orientation, width, length, movement sense, infilling materials, spacing, aperture, roughness of both joints and faults and intersection and connectivity to other joints were measured on outcrops. In addition, characteristics of subsurface fractures were examined by core logging in five boreholes, of which the orientations were acquired by acoustic televiewer logging from three boreholes. The dominant fracture sets were grouped from outcrops; GSet 1: N50-82$^{\circ}$E/55-90$^{\circ}$SE, GSet 2: N2-8$^{\circ}$E/56-86$^{\circ}$SE, GSet 3: N46-72$^{\circ}$W/60-85$^{\circ}$NE, GSet 4:Nl2-38$^{\circ}$W/15-40$^{\circ}$SW and from subsurface; HSet 1: N50-90$^{\circ}$E/55-90$^{\circ}$SE, HSet 2: N10-30$^{\circ}$E/50-70$^{\circ}$SE, HSet 3: N20-60$^{\circ}$W/50-80$^{\circ}$NE, HSet 4: N10-50$^{\circ}$E/$\leq$40$^{\circ}$NW. Among them, GSet 1, GSet 3 and HSet 1, HSet 3 are the most intensely developed fracture sets in the study area. The mean fracture spacings of HSet 1 are 30-47cm and code 1 fractures, such as faults and open fractures, comprise 21.0-42.9 percent of the whole fractures in each borehole. HSet 3 shows the mean fracture spacings of 55-57cm and the ratio of code 1 fractures is 15.4-26.9 percent. In spite of the mean fracture spacing of 239cm, code 1 fractures of HSet 4 have the highest ratio of 54.5 percent. From the fact that faults or open fractures have high hydraulic conductivity, it can be inferred that the three fracture sets of N55-85$^{\circ}$E/50-80$^{\circ}$SE, N20-60$^{\circ}$W/50-75$^{\circ}$NE and N10-30$^{\circ}$E/$\leq$30$^{\circ}$NW from a fracture system of relatively high conductivity. It is indirectly verified with geophysical loggings and constant injection tests performed in the boreholes.

• The Journal of Engineering Geology
• /
• v.14 no.1
• /
• pp.29-46
• /
• 2004

It is important to identify geometries of fracture that act as a conduit of fluid flow for characterization of ground water flow in fractured rock. Fracture geometries control hydraulic conductivity and stream lines in a rock mass. However, we have difficulties to acquire whole geometric data of fractures in a field scale because of discontinuous distribution of outcrops and impossibility of continuous collecting of subsurface data. Therefore, it is needed to develop a method to describe whole feature of a target fracture geometry. This study suggests a new approach to develop a method to characterize on the whole feature of a target fracture geometry based on the Fourier transform. After sampling of specimens along a target fracture from borehole cores, effective frequencies among roughness components were selected by the Fourier transform on each specimen. Then, the selected effective frequencies were averaged on each frequency. Because the averaged spectrum includes all the frequency profiles of each specimen, it shows the representative components of the fracture roughness of the target fracture. The inverse Fourier transform is conducted to reconstruct an averaged whole roughness feature after low pass filtering. The reconstructed roughness feature also shows the representative roughness of the target subsurface fracture including the geometrical characteristics of each specimen. It also means that overall roughness feature by scaling up of a fracture. In order to identify the characteristics of permeability coefficients along the target fracture, fracture models were constructed based on the reconstructed roughness feature. The computation of permeability coefficient was performed by the homogenization analysis that can calculate accurate permeability coefficients with full consideration of fracture geometry. The results show a range between $10^{-4}{\;}and{\;}10^{-3}{\;}cm/sec$, indicating reasonable values of permeability coefficient along a large fracture. This approach will be effectively applied to the analysis of permeability characteristics along a large fracture as well as identification of the whole feature of a fracture in a field scale.

• Journal of the Korean Society of Groundwater Environment
• /
• v.4 no.1
• /
• pp.27-40
• /
• 1997

The hydrogeochemical study on the 15 natural waters was carried out in the vicinity of tunnel excavation site of Donghae largely composed of granite and limestone. The water samples can be classified based on their chemical characteristics into two groups; waters draining in the granitic region(group 1) and the limestone region(group 2). This classification was also confirmed by statistical examination through cluster analysis, and the tunnel seepage waters collected at the same site appear to be included in group 1 and 2 by their sampling period, respectively. According to factor analysis, the waters of group 1 art mainly represented by the weathering of plagioclase to kaolinite and those of group 2 are characterized by the dissolution of calcite. Different properties of the tunnel seepage waters are thought to be resulted from the effective waterproofing processes conducted during the sampling interval to the surface and subsurface leakage zones at the granitic region, which contributed to the change of groundwater flow system. However both the tunnel seepage waters seem to have thermodynamically interacted with rock-forming minerals in their wallrocks. The mixing ratio of the waters from two groups and water-rock interactions are evaluated quantitatively for the tunnel seepage waters through the mass balance approach, and the results are identical with the previous conclusions in this study.

• Journal of Korea Water Resources Association
• /
• v.43 no.12
• /
• pp.1011-1027
• /
• 2010

In hydrologic modeling, prediction uncertainty generally stems from various uncertainty sources associated with model structure, data, and parameters, etc. This study aims to assess the parameter uncertainty effect on hydrologic prediction results. For this objective, a distributed rainfall-sediment yield-runoff model, which consists of rainfall-runoff module for simulation of surface and subsurface flows and sediment yield module based on unit stream power theory, was applied to the mesoscale mountainous area (Cheoncheon catchment; 289.9 $km^2$). For parameter uncertainty evaluation, the model was calibrated by a multi-objective optimization algorithm (MOSCEM) with two different objective functions (RMSE and HMLE) and Pareto optimal solutions of each case were then estimated. In Case I, the rainfall-runoff module was calibrated to investigate the effect of parameter uncertainty on hydrograph reproduction whereas in Case II, sediment yield module was calibrated to show the propagation of parameter uncertainty into sedigraph estimation. Additionally, in Case III, all parameters of both modules were simultaneously calibrated in order to take account of prediction uncertainty in rainfall-sediment yield-runoff modeling. The results showed that hydrograph prediction uncertainty of Case I was observed over the low-flow periods while the sedigraph of high-flow periods was sensitive to uncertainty of the sediment yield module parameters in Case II. In Case III, prediction uncertainty ranges of both hydrograph and sedigraph were larger than the other cases. Furthermore, prediction uncertainty in terms of spatial distribution of erosion and deposition drastically varied with the applied model parameters for all cases.

• Journal of Korean Society of Environmental Engineers
• /
• v.33 no.7
• /
• pp.544-551
• /
• 2011

The objective of this study was to analyze VIRULO model, a probabilistic quantitative model, which had been developed by US Environmental Protection Agency. The model could assess the viral attenuation capacity of soil as hydrogeologic barrier using Monte Carlo simulation. The governing equations used in the model were composed of unsaturated flow equations and viral transport equations. Among the model parameters, those related to water flow for 11 soil types were from UNDODA data, and those related to 5 virus species were from the literatures. The model compared the attenuation factor with threshold of attenuation to determine the probability of failure and presented the exceedances and Monte Carlo runs as output. The analysis indicated that among 11 USDA soil types, the viral attenuation capacity of loamy sand and sand were far lower than those of clay and silt soils. Also, there were differences in the attenuation in soil among 5 viruses with poliovirus showing the highest attenuation. The viral attenuation capacity of soil decreased sharply with increasing soil water content and increased nonlinearly with increasing soil barrier length. This study indicates that VIRULO model could be considered as a useful screening tool for viral risk assessment in subsurface environment.

• 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.