• Journal of Korea Water Resources Association
• /
• v.31 no.5
• /
• pp.553-564
• /
• 1998

For the proper design and management of underground storage cavern, groundwater flow around cavern should be analyzed. Since this flow is influenced by spatial nonuniformity of hudraulic conductivity, the two-dimensional finite element flow model incorporating stochastic concepts was developed to analyze influences due to this nonuniformity. Monte Carlo technique was applied to obtain an approximate solution for two-dimensional, steady flow in a stochastically defined nonuniform medisu. For this purpose, the values of hydraulic conductivity were generated for each element with known mean and standard deviations. The uncertainty in model prediction depends on both the nonuniformity in hydraulic conductivity and the natures of the flow system such as water curtain and boundary condition. Therefore the uncertainties in predicted hydraulic head and gradient are the greatest where the mean hydraulic gradients are relatively large and far from the boundaries. Especially, we relate these uncertainties with well known gas tightness condition.

• Journal of the Korean Geotechnical Society
• /
• v.34 no.11
• /
• pp.93-105
• /
• 2018

Precise investigation and interpretation of the ground subsidence risk factors needed to predict and evaluate the settlement problems of the surrounding ground due to the ground excavation. There are various geophysical exploration methods to investigate the ground subsidence risk factors. However, there are factors that influence the characteristics of the underground medium in these geophysical methods, and the actual soil contains complex factors affecting geophysical exploration. Therefore, it is necessary to analyze the effects on the geophysical methods based on the understanding of the geotechnical properties of soil. In this study, a test bed was constructed to consider various complicated factors in the complex ground and the ground behavior was analyzed by numerical analysis. In addition, we analyzed the limitations on investigating the ground subsidence risk factors through ground penetration radar (GPR) survey. As a result, ground subsidence of Open-cut Type Excavation is caused by various factors. Especially, in the case of soft ground condition, it was found that it was greatly influenced by the flow change of groundwater level. At the center frequency of GPR of 250 MHz, the attenuation of the electromagnetic wave is severely attenuated in the clay with high electrical conductivity, making it difficult to penetrate deeply into the ground (4 m below the surface). As the electromagnetic waves pass through the groundwater level below the groundwater, the attenuation of the electromagnetic waves becomes severe.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2001.04a
• /
• pp.81-84
• /
• 2001

고준위 방사성 폐기물 영구처분 안전성을 평가하기 위하여 입력 자료로 처분장 주변 각 암 반에서의 지하수 유동 속도 및 유동 시간이 요구된다. 이러한 유동 속도와 시간은 대부분의 경우 단일 값이 요구되지만 고준위 방사성 처분의 경우 지하 매질의 불확실성을 고려하기 위하여 확률론적 분석이 요구된다. 지하수 유동 속도 및 시간이 확률밀도함수로 표시되기 위해서는 기존의 방법에서는 수리 해석의 입력 인자 값들을 변화시키면서 반복적인 계산을 수행하는 방법이 사용되었다. 그러나 이러한 방법론의 한계를 극복하기 위해 최근 섭동 이론을 이용한 adjoint 방법론이 사용되고 있는 바 이를 이용하여 가상 처분장에서의 지하수 유동 속도와 시간을 확률론적으로 해석하였다.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.12 no.3
• /
• pp.285-291
• /
• 1987

Seismic waves are attenuated by losses of energy as they propagate through the earth. One way to model this numerically is to make the velocity a complex number, the real part giving the phase velocity and the imaginary part the attenuation. This models wave propagation in a medium for which the logarithmic decrement is independent of frequency(attenuation coefficient is proportional to frequncy). The aim is to modify forward and inverse numerical modeling so that attenuation can be specified as a function of position.

• Korean Journal of Soil Science and Fertilizer
• /
• v.41 no.5
• /
• pp.330-341
• /
• 2008

The purpose of this study was to investigate the vulnerability for groundwater contamination at the some cemeteries in Gyeonggi Province. Twenty-eight out of 43 cemeteries in Gyeonggi province were selected for this study. The DRASTIC model was applied to those cemeteries, and the reliance of the model was assessed using the water quality data of the target areas. The DRASTIC model was used for the assessment of the potential for groundwater contamination using hydrogeological factors. Seven factors including depth of water, net recharge, aquifer media, soil media, topography, impact of the vadose zone, hydraulic conductivity of the aquifer were assessed. The DRASTIC index of the study area ranged from 82 to 126 with an average value of $113.99(\pm11.48)$. The DRASTIC index was relatively greater in the northern Gyeonggi province than that in the southern area. The DRASTIC index was similar for the areas with the similar burial rate and burial density. This study demonstrated that burial rate and burial density should be considered along with the 7 basic factors for the evaluation of groundwater vulnerability of the cemeteries.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.34-44
• /
• 2001

Forward modeling by construction of synthetic data is usually practiced in a horizontal surface and a few subsurface structures. However, in-situ surveys often take place in such topographic changes that the corrupted field data always make it difficult to interpret the right signals. To examine the propagation characteristic of elastic waves on the irregular surface, a general mesh generation code for finite element method was modified to consider the topography. By implementing this algorithm, the time domain modeling was practiced in some models with surface topography such as mound, channel, etc. The synthetic data obtained by receivers placed on surface also agreed with the analytic solution. The snapshots showing the total wave-field revealed the propagation characteristic of the elastic waves through complex subsurface structures and helped to identify the signals on the time traces. The transmission of Rayleigh waves along the surface, compressive waves, and sheer waves was observed. Moreover, it turned out that the Rayleigh waves behave like a new source at the edge.

• Geophysics and Geophysical Exploration
• /
• v.18 no.4
• /
• pp.232-240
• /
• 2015

Magnetotelluric (MT) survey investigates electrical structure of subsurface by measuring natural electromagnetic fields on the earth surface. For the accurate interpretation of MT data, the precise three-dimensional (3-D) modeling algorithm is prerequisite. Since MT responses are affected by electrical anisotropy of medium, the modeling algorithm has to incorporate the electrical anisotropy especially when analyzing time-lapse MT data sets, for monitoring engineered geothermal system (EGS) reservoir, because changes in different-vintage MT-data sets are small. This study developed a MT modeling algorithm for the simulation MT responses in the presence of electrical anisotropy by improving a pre-existing staggered-grid finite-difference MT modeling algorithm. After verifying the developed algorithm, we analyzed the effect of vertical transversely isotropic (VTI) anisotropy on MT responses. In addition, we are planning to extend the applicability of the developed algorithm which can simulate not only the horizontal transversely isotropic (HTI) anisotropy, but also the tiled transversely isotropic (TTI) anisotropy.

• Tunnel and Underground Space
• /
• v.17 no.2 s.67
• /
• pp.83-89
• /
• 2007

Impact-echo test is a non-destructive testing method to determine dynamic properties of a material. This presentation introduces the experimental set-up and procedure of the test for rock specimens. In addition, the test results of domestic rocks collected in 5 different areas, a cement mortar and aluminium alloy are presented. The test results include resonance frequencies of P- and S-wave as well as damping ratios of the described 7 different materials. The differences between dynamic and static values of elastic moduli are about 10%, while the dynamic Poisson's ratios are greater than the static Poisson's ratios by at least 0.07. The damping ratio is dependent on the joint density and degree of weathering of a rock specimen.

• Geophysics and Geophysical Exploration
• /
• v.17 no.4
• /
• pp.231-241
• /
• 2014

Recently, single channel seismic survey for engineering purpose have been used widely taking advantage of simple processing. However it is very difficult to obtain high fidelity subsurface image by single channel seismic due to insufficient fold coverage. Recently, seismic waveform inversion in multi channel seismic survey is utilized for accurate subsurface imaging even in complex terrains. In this paper, we propose the seismic waveform inversion algorithm for velocity model building using a single channel seismic data. We utilize the Gauss-Newton method and assume that subsurface model is 1-Dimensional. Seismic source estimation technique is used and offset effect is also corrected by removing delay time by offset. Proposed algorithm is verified by applying modified Marmousi2 model, and applied to field data set obtained in port of Busan.

• The Journal of Engineering Geology
• /
• v.6 no.3
• /
• pp.137-153
• /
• 1996

The purpose of this study is to investigate the possibilities of identifying and detecting underground cavities using seismic waves recorded by the fixed and mobile stations. During 18 months of field work we recorded chemical explosions near the Bongdarn station. Seismic Stations were installed on the free surface and underground inside the Samba mine. The seismograms at the fixed(lorg-term) seismic station show abrupt change of polarization characteristics which can he associated with the appearance of P-to-S converted phase(PS) at 150 ~ 200 msec after the first P arrival. This result indicates that converted phases are generated very near to the Bongdarn station at a depth of 190m. Shear-wave splitting phenomena have also been observeci The time delay between fast shear(fS) and slow shear(sS) waves ranges between 30 and 60 msec(average is 42 msec). However, exact time delay between the fast and the slow shear waves can not be accurately measured because of the very short time delay and limitation of sampling rate. Chemical explosion experiments were recorded at stations along various paths to contrast the seismic response of areas with and without cavities. The seismograms recorded at the stations installed at cavity areas show an abrupt change of polarization characteristics but not on the other stations. Seismic waves propagating through the cavity are characterized by the attenuation of high frequency waves and predominantly low frequency seismic waves after the S wave arrivals.