• The Journal of Engineering Geology
• /
• v.13 no.1
• /
• pp.51-65
• /
• 2003

A series of numerical simulations were carried out using STOMP( Subsurface Transport over Multiple Phase) simulator. The flow and distribution of LNAPL were analyzed in homogeneous fine and coarse sand. Vertical movement of LNAPL is faster in the coarser sand. But the total volume of LNAPL retained in the unsaturated zone is larger in the finer sand. A fine layer in the coarse sand domain is also simulated. The results showed that the retained LNAPL volume and shape are highly influenced by the Position of the fine layer. Flow and distributions of LNAPL were simulated when there were heterogeneous lenses in the sand domain. Water table fluctuation was also considered. In these cases, it was found that the heterogeneous lens was a barrier to LNAPL flow, and water table fluctuation stimulated the downward movement of retained LNAPL. The LNAPL flow and distribution observed in these numerical experiments show that in the subsurface environment, the behaviors of LNAPL highly depend on heterogeneities of unsaturated zone and the dynamic hydrogeologic condition such as water table fluctuation. These results can explain some of the complexity of LNAPL flow and distribution Patterns in LNAPL contaminated field sites.

• 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이 조립질 렌즈를 통해 이동하나, 세립질 렌즈를 통해서는 이동하지 못한다.

• 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 Korea Water Resources Association
• /
• v.36 no.6
• /
• pp.1035-1046
• /
• 2003

In order to evaluate the TOPMODEL's prediction ability for spatial distribution of water table depths, two major assumptions and governing equation of water table depth are tested. For the test, data of hydrological observations are used and a soil survey is made in the steep hillslope with thin soils. Responses of water table and hydraulic properties of soil are coincident with two major assumptions of the TOPMODEL's such as water table gradient parallel to the local topographical slope and exponential decline in transmissivity with depths. Soil texture and the decline rate of transmissivity(f) we homogeneous in space at the 0∼0.3m depths of the soil of the hillslope, but they are heterogeneous in space below its 0.3m depths due to the vertical change of soil texture and the ‘f’. It is shown that the TOPMODEL's equation can be used for simulating distribution of water table depth at the depths with uniform values of the 'f'.

• Geophysics and Geophysical Exploration
• /
• v.8 no.1
• /
• pp.7-17
• /
• 2005

The determination of seismic velocities in refractors for near-surface seismic refraction investigations is an ill-posed problem. Small variations in the computed time parameters can result in quite large lateral variations in the derived velocities, which are often artefacts of the inversion algorithms. Such artefacts are usually not recognized or corrected with forward modelling. Therefore, if detailed refractor models are sought with model based inversion, then detailed starting models are required. The usual source of artefacts in seismic velocities is irregular refractors. Under most circumstances, the variable migration of the generalized reciprocal method (GRM) is able to accommodate irregular interfaces and generate detailed starting models of the refractor. However, where the very-near-surface environment of the Earth is also irregular, the efficacy of the GRM is reduced, and weathering corrections can be necessary. Standard methods for correcting for surface irregularities are usually not practical where the very-near-surface irregularities are of limited lateral extent. In such circumstances, the GRM smoothing statics method (SSM) is a simple and robust approach, which can facilitate more-accurate estimates of refractor velocities. The GRM SSM generates a smoothing 'statics' correction by subtracting an average of the time-depths computed with a range of XY values from the time-depths computed with a zero XY value (where the XY value is the separation between the receivers used to compute the time-depth). The time-depths to the deeper target refractors do not vary greatly with varying XY values, and therefore an average is much the same as the optimum value. However, the time-depths for the very-near-surface irregularities migrate laterally with increasing XY values and they are substantially reduced with the averaging process. As a result, the time-depth profile averaged over a range of XY values is effectively corrected for the near-surface irregularities. In addition, the time-depths computed with a Bero XY value are the sum of both the near-surface effects and the time-depths to the target refractor. Therefore, their subtraction generates an approximate 'statics' correction, which in turn, is subtracted from the traveltimes The GRM SSM is essentially a smoothing procedure, rather than a deterministic weathering correction approach, and it is most effective with near-surface irregularities of quite limited lateral extent. Model and case studies demonstrate that the GRM SSM substantially improves the reliability in determining detailed seismic velocities in irregular refractors.

• The Journal of Engineering Geology
• /
• v.7 no.2
• /
• pp.113-126
• /
• 1997

Electrical resistivity methods of dipole - dipole array profiling and Schiumberger array sounding were tested on a segment of the Woraksan granitic batholith for the research into the imaging of irregular attitudes of fracture zones in the crystaaline rock in terms of processing and interpretation schemes. By the dipole - dipole array method, inhomogeneities such as small scale of fracture zones were properly delineated down at some depth even within hard rock environment. Fracture zones were interpreted to be at the boundaries between the high amplitude zone and very low amplitude zone in the resistivity plot and they were also successfully outlined in two - dimensional layer and pseudo - three - dimensional volume constructed by the incorporation of vertical sounding data. The surface location of the fracture zones was correlated by the zero - crossing point in the VLF(very low frequency) electromagnetic data. Pseudo - three - dimensional attitudes of fracture zones were efficiently illuminated by optimum projection angle. The mean of bulk resistivity for the Woraksan granite and the near fracture zones is estimated to be approximately of 4,000 ohm - m which is much higher than the value of 700 ohm - m for the Rwachunri limesilicate environment. This difference is due to both the rock type, i.e., biotite granite vs limesilicate, and the occurrence of secondary openings of fold and fault associated with the intrusion of granite. In this study statistical analyses on the resistivity color plot were performed in terms of three representative statistical moments, i.e., standard deviation, skewness, and kurtosis. The fracture zones in the standard deviation plot were characterized by the higher value, compared to the value of homogeneous portion. The upper boundary of the high resistivity zone was also successfully delineated in the skewness and kurtosis plots.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.17 no.2
• /
• pp.127-140
• /
• 2015

A fault is one of the critical factors that may lead to a possible ground collapse occurring in construction site. A fault core, however, possibly acting as a failure plane in whole fault zone, is composed of fractured rock and gouge nonuniformly distributed and thus can be characterized by its wide range of shear strength which is generally acquired by experimental method for stability analysis. In this study, we performed direct shear test and grain size distribution analysis for 62 fault core samples cropped from 12 different spots located in the vicinity of Kyongju and Ulsan, Korea. As a result, the range of shear strength representing the characteristics of fault cores in the study regions is determined with regard to vertical stress using a regression analysis for experiment data. The weight ratio of gravels in the samples is proportional to the shear strength and that of silt and clay is in inverse proportion to the shear strength. For most samples, the coefficient of determination is over 0.7 despite of inhomogeneity of them and consequently we determined the lower limit and upper limit of the shear strength with regard to the weight ratio by setting the confidence interval of 95%.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.11 no.3
• /
• pp.179-192
• /
• 2013

This study is characterized the groundwater flow pattern near crystalline site of Yangbook-Myeon, Gyeong-ju City, South Korea. From the time series analyses, groundwater level could be classified into 4 types reflecting the hydrogeological characteristics and rainfall pattern. The type I (DB1-1, DB1-2) may be directly influenced by rainfall pattern. The type II (DB1-3, DB1-7, KB-1, KB-2, KB-3, KB-7, KB-14, KB-15) may be influenced by rainfall event as well as groundwater flow through water-conducting features. The type III (DB-5, DB1-6, DB2-2, KB-10, KB-11, KB-13) may be predominantly happens in the crystaline rock mass, groundwater in this type flows through the minor fracture networks rather than direct effect of rainfall event. The type IV (DB1-8, KB-9) may be influenced by irregular variation of the groundwater level due to anisotropy and heterogeneity of crystalline rock.

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