• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.68-77
• /
• 2008

We have developed an inversion algorithm for loop-loop electromagnetic (EM) data, based on the localised non-linear or extended Born approximation to the solution of the 2.5D integral equation describing an EM scattering problem. Source and receiver configuration may be horizontal co-planar (HCP) or vertical co-planar (VCP). Both multi-frequency and multi-separation data can be incorporated. Our inversion code runs on a PC platform without heavy computational load. For the sake of stable and high-resolution performance of the inversion, we implemented an algorithm determining an optimum spatially varying Lagrangian multiplier as a function of sensitivity distribution, through parameter resolution matrix and Backus-Gilbert spread function analysis. Considering that the different source-receiver orientation characteristics cause inconsistent sensitivities to the resistivity structure in simultaneous inversion of HCP and VCP data, which affects the stability and resolution of the inversion result, we adapted a weighting scheme based on the variances of misfits between the measured and calculated datasets. The accuracy of the modelling code that we have developed has been proven over the frequency, conductivity, and geometric ranges typically used in a loop-loop EM system through comparison with 2.5D finite-element modelling results. We first applied the inversion to synthetic data, from a model with resistive as well as conductive inhomogeneities embedded in a homogeneous half-space, to validate its performance. Applying the inversion to field data and comparing the result with that of dc resistivity data, we conclude that the newly developed algorithm provides a reasonable image of the subsurface.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.51 no.5
• /
• pp.25-34
• /
• 2009

The objective of this study is to evaluate the future potential climate and vegetation canopy change impact on a dam watershed hydrology. A $6,661.5\;km^2$ dam watershed, the part of Han-river basin which has the watershed outlet at Chungju dam was selected. The SWAT model was calibrated and verified using 9 year and another 7 year daily dam inflow data. The Nash-Sutcliffe model efficiency ranged from 0.43 to 0.91. The Canadian Centre for Climate Modelling and Analysis (CCCma) Coupled Global Climate Model3 (CGCM3) data based on Intergovernmental Panel on Climate Change (IPCC) SRES (Special Report Emission Scenarios) B1 scenario was adopted for future climate condition and the data were downscaled by artificial neural network method. The future vegetation canopy condition was predicted by using nonlinear regression between monthly LAI (Leaf Area Index) of each land cover from MODIS satellite image and monthly mean temperature was accomplished. The future watershed mean temperatures of 2100 increased by $2.0^{\circ}C$, and the precipitation increased by 20.4 % based on 2001 data. The vegetation canopy prediction results showed that the 2100 year LAI of deciduous, evergreen and mixed on April increased 57.1 %, 15.5 %, and 62.5% respectively. The 2100 evapotranspiration, dam inflow, soil moisture content and groundwater recharge increased 10.2 %, 38.1 %, 16.6 %, and 118.9 % respectively. The consideration of future vegetation canopy affected up to 3.0%, 1.3%, 4.2%, and 3.6% respectively for each component.

• Geomechanics and Engineering
• /
• v.14 no.5
• /
• pp.467-477
• /
• 2018

Estimation of hydraulic parameters is a critical step during design of foundation dewatering works. When many piles are installed in an aquifer, estimation of the hydraulic conductivity should consider the blocking of groundwater seepage by the piles. Based on field observations during a dewatering project in Shanghai, hydraulic conductivities are back-calculated using a numerical model considering the actual position of each pile. However, it is difficult to apply the aforementioned model directly in field due to requirement to input each pile geometry into the model. To develop a simple numerical model and find the optimal hydraulic conductivity, three scenarios are examined, in which the soil mass containing the piles is considered to be a uniform porous media. In these three scenarios, different sub-regions with different hydraulic conductivities, based on either automatic inverted calculation, or on effective medium theory (EMT), are established. The results indicate that the error, in the case which determines the hydraulic conductivity based on EMT, is less than that determined in the automatic inversion case. With the application of EMT, only the hydraulic conductivity of the soil outside the pit should be inverted. The soil inside the pit with its piles is divided into sub-regions with different hydraulic conductivities, and the hydraulic conductivity is calculated according to the volume ratio of the piles. Thus, the use of EMT in numerical modelling makes it easier to consider the effect of piles installed in an aquifer.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.10 no.1
• /
• pp.37-47
• /
• 2008

Design and construction of tunnels require understanding the influence of groundwater. Particularly, it is essential to know how the drainage conditions at the tunnel boundary affect flow behavior of ground adjacent to the tunnels. In this study flow behavior of a leaking tunnel was investigated using physical model tests for tunnel depths and various hydraulic boundary conditions. Particular concerns were given to flow lines toward tunnels. Test results showed that the boundary conditions hardly influence on flow patterns and time required to reach steady state conditions. It is revealed that with an increase in water depth, flow lines concentrated to the drain holes. The physical tests were numerically simulated. Numerical results showed that the flow behavior was represented appropriately by considering filter-drain hole drainage rather than boundary drainage all over the lining.

• Explosives and Blasting
• /
• v.27 no.1
• /
• pp.21-31
• /
• 2009

An excavation damaged zone (EDZ) is created by fracturing, excavation or stress redistribution of tunnels. In this zone the mechanical and hydraulic properties of rock are changed, which makes additional cracks and serves as a dominant pathway of groundwater flow. In this study, an assessment on an EDZ size was practiced by the measurement of the deformation modulus at the KAERI underground research tunnel (KURT), and the information was applied to the modelling analysis using FLAC2D software. The EDZ at KURT fell into the range of 0.6~1.8m and the deformation moduli of the EDZ generally correspond to about 40% of intact rock mass. With a consideration of the EDZ in numerical analysis, tunnel displacements increased by about 65% and the maximum principal stress decreased to 58% from the case without EDZ. The plastic zone of the tunnel was enlarged to the crown and invert zones of the tunnel within the range of the length of rock bolts. About 2% of the total tunnel displacement with EDZ was suppressed by the KURT support system. It is anticipated that the investigation of an EDZ can be used as an important and fundamental research for validating the overall performance of a high level waste disposal system.

• Geophysics and Geophysical Exploration
• /
• v.3 no.3
• /
• pp.88-93
• /
• 2000

The present study may be summarized as follows. By means of a pole-pole electrical survey, major geological features, i.e, fresh rocks, fractured & weak & saturated rocks, fresh water bearing aquifer, were successfully delineated in the surveyed area of a granite region in Korea. The subsequent Schlumberger sounding and drilling confirmed the existence of the acquifer at 60 m depth as expected. But one more minor acquifer which does not show up in the resistivity depth section was met at 100 m depth. A simple forward modelling leads the authors to believe that any other electrode configuration, e.g., Wenner, Schlumberger, dipole-dipole, pole-dipole, would not detect the deeper aquifer. Under these circumstances, further studies remain to be done in connection with the spatial resolution in the vertical direction.

• Geophysics and Geophysical Exploration
• /
• v.8 no.1
• /
• pp.73-77
• /
• 2005

To correctly invert and interpret Surface Nuclear Magnetic Resonance (SNMR) data collected in conductive terrains, an accurate estimate of subsurface conductivity structure is required. Given such an estimate, it would be useful to determine, before conducting an SNMR sounding, whether or not the conductivity structure would prevent groundwater being detected. Using SNMR forward modelling, we describe a method of determining the depth range from which most of the SNMR signal originates, given a model of subsurface conductivity structure. We use the method to estimate SNMR depth penetration in a range of halfspace models and show that for conductive halfspaces ($<10{\Omega}.m$) the depth of penetration Is less than 50 m. It is also shown that for these halfspaces, increasing coincident loop size does not significantly improve depth penetration. The results can be used with halfspace approximations of more complicated ID conductivity structures to give a reasonable estimate of the depth range over which signal is obtainable in conductive terrains.

• Geomechanics and Engineering
• /
• v.28 no.3
• /
• pp.283-298
• /
• 2022

Landslides triggered by the combination of heavy precipitation and anthropological disturbance in hilly areas cause severe damage to human lives, properties, and infrastructure constructions. A comprehensive investigation of the influencing factors and failure mechanisms of landslides are significant for disaster mitigation and prevention. This paper utilized the combination of detailed geological investigation, physical experimental testing as well as numerical modelling to determine the failure mechanism, and proposed a countermeasures of the Lantian landslide occurred on 2, July 2017. The results reveal that the Lantian landslide is a catastrophic reactivated slide which occurred in an active tectonic region in Southwest China. Because of the unique geological settings, the fully to highly weathered basalts in the study area with well-developed fractures favored the rainwater infiltration, which is the beneficial to slide reactivation. Engineering excavation and heavy precipitation are the main triggering factors to activate the slide motion. Two failure stages have been identified in the landslide. The first phase involves a shallow mass collapse originated at the upper slopes, which extends from the road to platform at rear part, which is triggered by excavation in the landslide region. Subjected to the following prolonged rainfall from 19 June to 2 July, 2017, the pore water pressure of the slope continually increased, and the groundwater table successively rise, resulting in a significant decrease of soil strength which leads to successive large-scale deep slide. Thereinto, the shallow collapse played a significant role in the formation of the deep slide. Based on the formation mechanisms of the landslide, detailed engineering mitigation measures, involving slope cutting, anchor cable frame, shotcrete and anchorage, retaining wall and intercepting ditch were suggested to reduce the future failure risk of the landslide.

• Journal of Korea Water Resources Association
• /
• v.55 no.11
• /
• pp.923-929
• /
• 2022

The Chuncheon Mullori area is an underprivileged area of water welfare where local water supply is not supplied, and it is supplying water to the villages with small water supply facilities using lateral flow and groundwater as water sources. This is an area with poor water supply conditions, such as relying on water trucks due to water shortages during the recent severe drought. Therefore, in order to solve the problem of water shortage during drought and to prepare for the increasing water demand, a sand dam was installed along the valley, and this facility has been operating since May 2022. In this study, repeated simulations were performed according to the hydraulic conductivity of the filler material and the storage coefficient value for the inflow condition for about two years from mid-March 2020 to mid-March 2022. For each case, the amount of discharge through the perforated drain pipe was calculated. Overall, as the hydraulic conductivity increased, the amount of discharge and its ratio increased. However, when the hydraulic conductivity of the second floor was relatively low, the amount of discharge increased and then decreased as the hydraulic conductivity of the third floor increased. This is considered to be due to the fact that the water level was kept low due to the rapid drainage compared to the net inflow into the third floor because the water permeability of the third floor and the drainage coefficient of the drain pipe were large. As a result of simulating the flow of the open channel in the upper part of the sand dam as a hypothetical groundwater layer with very high hydraulic conductivity, the decrease in discharge rate was slower than the increase in the hydraulic conductivity of the hypothetical layer, but it was clearly shown that the discharge volume decreased relatively as the hydraulic conductivity of the virtual layer increased.

• Journal of the Geological Society of Korea
• /
• v.54 no.5
• /
• pp.545-556
• /
• 2018

This study was conducted to estimate the deep percolation using numerical modeling and field observation data based on rainfall in Haean basin. Soil moisture sensors were installed to monitoring at 30, 60 and 90 cm depths in four sites (YHS1-4) and automatic weather station was installed to around YHS3. Soil moisture and meteorological data was observed from March 25, 2017 to March 25, 2018 and May 06, 2016 to May 06, 2018, respectively. Numerical analysis was performed from June to August, 2017 using the HYDRUS-1D. Average soil moisture contents were high to generally in YHS3 for 0.300 to $0.334m^3/m^3$ and lowest in YHS1 for 0.129 to $0.265m^3/m^3$ during the soil moisture monitoring period. The results of soil moisture flow modeling showed that field observations and modeling values were similar but the peak values were larger in the modeling result. Correlation analysis between observation and modeling data showed that r, $r^2$ and RMSE were 0.88, 0.77, and 0.0096, respectively. This show high correlation and low error rate. The total deep percolation was 744.2 mm during the period of modelling at 500 cm depth. This showed that 61.3% of the precipitation amount (1,214 mm) was recharged in 2017. Deep percolation amount was high in the study area. This study is expected to provide basic data for the estimation of groundwater recharge through unsaturated zone.