• The Journal of Engineering Geology
• /
• v.10 no.3
• /
• pp.247-261
• /
• 2000

The geostatistical analyses for the chemical components of pH, TS, KMnO4 Demand, Cl, SO$_4$ and NO$_3$-N are carried out to understand the groundwater contamination in Pusan. The average values of each component are 7.2 for pH, 336.4mg/$\ell$ for TS, 2.3mg/$\ell$ for KMnO$_4$ Demand, 44.3mg/$\ell$ for Cl, 36.0mg/$\ell$ for SO$_4$, and 4.6mg/$\ell$ for NO$_3$-N. The ratios over the drinking standard of each component are 0.34% for pH, 2.27% for TS, 1.55% for KMnO$_4$ Demand, 1.59% for Cl, 0.57% for SO$_4$, and 3.7% for NO$_3$-N. The highest ratio of NO$_3$-N results from the municipal sewage and exhaust gas of vehicles. The isopleth maps of 6 chemical components show that the high values of groundwater contamination come from the inland of Pusan, and that some high values appear at the coastal area. The isopleth maps of Cl and SO$_4$ related with seawater intrusion also show that the high values appear only at the particular coastal area, not at the whole area. On the isopleth maps of Cl and SO$_4$, the anomalies of the concentration contours were compared with the directions of two large fault zones, the Ilkwang Fault and the Dongrae Fault. Apparently, they don't have the particular correlation. Therefore, it is concluded that the main source of groundwater contamination in Pusan is not the seawater, but the municipal sewage and other sources such as the exhaust gas of vehicles, the contaminated surface water, the waste water of factories, and the leachate of waste landfills.

• Journal of the Korea Concrete Institute
• /
• v.15 no.4
• /
• pp.574-584
• /
• 2003

The most common deterioration cause of concrete structures in urban environment is carbonation. Recently, the $CO_2$ concentration and temperature at atmosphere is sharply increased with time due to global warming phenomena. In this study, the climate scenario IS92a, which was suggested by the IPCC, is used to consider temperature and atmospheric $CO_2$ concentration change in the model of service life prediction. The modified mathematical solution, which was based on the Fick's 1st law of diffusion, was used to reflect concrete materials properties such as the degree of hydration of concrete with elapsed time, and important parameters, which associated with deterioration rate. The techniques of service life prediction are developed introducing the method of reliability and stochastic concept to consider microclimatic condition in Seoul, South Korea. From the result of service life prediction, concrete containing high W/C ratio is shown fast carbonation rate due to $CO_2$ concentration increase. It is concluded that the deterioration of concrete structures due to carbonation is insignificant problem on the conditions that below W/C 55%, well curing concrete.

• Korean Journal of Remote Sensing
• /
• v.20 no.6
• /
• pp.383-396
• /
• 2004

The application of spatial statistics to obtain the spatial uncertainty distributions in classification of remote sensing images is investigated in this paper. Two quantitative methods are presented for describing two kinds of uncertainty; one related to class assignment and the other related to the connection of reference samples. Three quantitative indices are addressed for the first category of uncertainty. Geostatistical simulation is applied both to integrate the exhaustive classification results with the sparse reference samples and to obtain the spatial uncertainty or accuracy distributions connected to those reference samples. To illustrate the proposed methods and to discuss the operational issues, the experiment was done on a multi-sensor remote sensing data set for supervised land-cover classification. As an experimental result, the two quantitative methods presented in this paper could provide additional information for interpreting and evaluating the classification results and more experiments should be carried out for verifying the presented methods.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.640-640
• /
• 2015

2010년과 2011년 서울에서 발생한 집중호우와 2014년 부산에서 발생한 집중호우의 발생으로 막대한 재산상의 피해와 사상자를 냈다. 2010년 9월 21일에 발생한 집중호우는 1908년 관측시작이래 가장 많은 비가 내린 것으로 기록되었으며 주거지 4,727호, 상가 1,164호, 공장 126동 등이 침수되고 13시를 기준으로 강서지점의 경우 시간당 98.5mm의 기록적인 강우를 기록하였으나, 관악지점은 5.5mm에 그쳐 두 지점간의 시간당 강우량의 편차가 약 200배 가까이 차이가 나는 것으로 나타났다. 이와 같이 최근 도시지역에서 국지성 집중호우가 증가하고 있으며 지역별 강우 편차가 크고 이에 따라 침수피해발생 여부도 지역에 따라 달라진다. 강수의 공간적 분포와 그로 인한 침수해석은 도시돌발홍수 예경보 시스템에 있어 무엇보다도 중요하다. 본 연구의 목적은 도시지역 돌발홍수 예경보 시스템 구축을 위한 정량적 강수추정 QPE(Quantitative Precipitation Estimation)기법에 따른 수리 수문학적 영향을 평가하는 것이다. 정량적 강수추정을 위해 AWS, SKP, 레이더 자료를 활용하여 250m의 해상도를 가지도록 크리깅을 적용하였다: QPE 1은 34개의 AWS의 지점우량을 지구통계학적 기법 중의 하나인 크리깅을 이용하여 산정한 기법, QPE 2는 AWS와 156개의 SKP의 강우데이터를 크리깅을 이용하여 산정한 기법, QPE 3는 광덕산 레이더를 이용한 기법, QPE 4는 AWS, SKP, 광덕산 레이더 자료를 조건부 합성한 기법이다. 월류량을 산정하기 위해 도시유출해석모형인 SWMM을 강남역 일대를 대상으로 구축하고 우수관로 시스템으로 유입되지 못한 노면류(Surface flow)를 함께 고려하였다. 침수해석을 위해서는 DHM모델을 적용하였으며 2013년 7월 기간에 발생한 호우에 대하여 분석을 수행하였다. 비교수행을 위해서 인접한 서초 AWS와 강남 AWS의 지점강우량도 함께 고려하였으며 모의결과를 국가 재난관리 정보 시스템(NMDS)에 침수피해가 확인된 가옥 및 빌딩 정보와 일치여부를 적합도로 산정하였다. 산정된 적합도를 통하여 정량적 강수추정기법에 따른 수리?수문학적 영향을 평가하였다. 실제 침수흔적정보와 비교 결과, QPE 2와 QPE4가 가장 적합도가 높았으며 이에 따라 고밀도의 관측망의 구성이 도시지역 침수해석결과에도 적합할 것으로 판단된다.

• Journal of the Korean GEO-environmental Society
• /
• v.4 no.3
• /
• pp.79-88
• /
• 2003

This paper presents a numerical study of the spatial behavior of a linear absorbing solute in a heterogeneous porous medium. The spatially correlated log-normal hydraulic conductivity field is generated in a given two-dimensional domain by using the geostatistical method (Turning Bands algorithm). The velocity vector field is calculated by applying the two-dimensional saturated groundwater flow equation to the Galerkin finite element method. The simulation of solute transport is carried out by using the random walk particle tracking model with CD(constant displacement) scheme in which the time interval is automatically adjusted. In this study, the spatial behavior of a solute is analyzed by the longitudinal center-of-mass displacement, longitudinal spatial spread moment and longitudinal plume skewness.

• Journal of the Korean Geophysical Society
• /
• v.3 no.3
• /
• pp.161-174
• /
• 2000

Bayesian inversion is a stable approach to infer the subsurface structure with the limited data from geophysical explorations. In geophysical inverse process, due to the finite and discrete characteristics of field data and modeling process, some uncertainties are inherent and therefore probabilistic approach to the geophysical inversion is required. Bayesian framework provides theoretical base for the confidency and uncertainty analysis for the inference. However, most of the Bayesian inversion require the integration process of high dimension, so massive calculations like a Monte Carlo integration is demanded to solve it. This method, though, seemed suitable to apply to the geophysical problems which have the characteristics of highly non-linearity, we are faced to meet the promptness and convenience in field process. In this study, by the Gaussian approximation for the observed data and a priori information, fast Bayesian inversion scheme is developed and applied to the model problem with electric well logging and dipole-dipole resistivity data. Each covariance matrices are induced by geostatistical method and optimization technique resulted in maximum a posteriori information. Especially a priori information is evaluated by the cross-validation technique. And the uncertainty analysis was performed to interpret the resistivity structure by simulation of a posteriori covariance matrix.

• Journal of the Korean earth science society
• /
• v.32 no.6
• /
• pp.525-536
• /
• 2011

In any applications using various types of spatial data, it is very important to account for the scale differences among available data sets and to change the scale to the target one as well. In this paper, we propose to use a geostatistical downscaling approach based on vaiorgram deconvloution and block simulation to generate fine-scale categorical thematic maps from coarse-scale fraction data. First, an iterative variogram deconvolution method is applied to estimate a point-support variogram model from a block-support variogram model. Then, both a direct sequential simulation based on area-to-point kriging and the estimated point-support variogram are applied to produce alternative fine-scale fraction realizations. Finally, a maximum a posteriori decision rule is applied to generate the fine-scale categorical thematic maps. These analytical steps are illustrated through a case study of land-cover mapping only using the block fraction data of thematic classes without point data. Alternative fine-scale fraction maps by the downscaling method presented in this study reproduce the coarse-scale block fraction values. The final fine-scale land-cover realizations can reflect overall spatial patterns of the reference land-cover map, thus providing reasonable inputs for the impact assessment in change of support problems.

• Journal of Soil and Groundwater Environment
• /
• v.7 no.3
• /
• pp.3-17
• /
• 2002

Data analysis of groundwater monitoring wells and geostatistical methods are used to identify the local characteristics of sea water intrusion and the range of sea water intrusion at the southeastern coastal area of Busan, Korea. Rainfall and groundwater level of two monitoring wells show a linear correlation because of the direct groundwater recharge by the precipitation. However, rainfall and electric conductivity have the inverse relationship because of the increase of groundwater. Electric conductivity rapidly increased at 24m depth and exceeded 20,000$\mu\textrm{s}$/cm near 26m depth in the monitoring wells. The variations of groundwater level and electric conductivity show that the interface between sea water and fresh water tends to move upward when groundwater level goes down. In the cross correlation analysis, groundwater level versus rainfall represents the largest cross correlation coefficient in 0 time lag but the cross correlation coefficient of electric conductivity versus rainfall is the largest when the time lag is 9 days. This suggests that the fluctuations of groundwater level respond to rainfall in a short time, but the interface between sea water and fresh water respond very slow to rainfall. Horizontal extents of sea water intrusion are estimated to 14 m from the east of Line 1, and 25 m from the southeast end of Line 2 in the inversion of dipole-dipole profiling data of two survey lines. The data of VES by the Schulumberger array in May and July show lognormal distributions. In the kriged apparent resistivity and earth resistivity distributions, the resistivities of July are increased comparing to those of May. This reflects that the concentration of sea water in aquifer is reduced due to the increased groundwater recharge from the rainfall in June and July. In analyzing the vertical and horizontal apparent resistivities and earth resistivity distributions, the geostatistical methods are very useful to identify the variations of earth resistivity distributions at the coastal area.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.4 no.3
• /
• pp.217-227
• /
• 2002

Geostatistical methods were used for the groundwater flow analysis on the ${\bigcirc}{\bigcirc}$ tunnel area. Linear regression analysis shows that the topographic elevation and ground water level of this area have very high correlation. Groundwater-level contour maps produced by ordinary kriging and cokriging have little differences in mountain areas. But, comparing two maps on the basis of an elevation contour map, a groundwater-level contour map using cokring is more accurate. Analyzing the groundwater flow on two groundwater-level contour maps, the groundwater of study area flows from the north-west mountain areas to near valleys, and from the peak of the mountain to outside areas. In the design steps, the groundwater-level distribution is reasonably considered in the tunnel construction area by cokriging approach. And, geostatistics will provide quantitative information in the unknown groundwatrer-level area.

• Tunnel and Underground Space
• /
• v.17 no.6
• /
• pp.493-502
• /
• 2007

Measurements of hydraulic conductivity are point or interval values, and are highly limited in their number. Meanwhile, results of geophysical prospecting can provide the information of spatial variation of geology, and abundant in number. In this study, it was aimed to develop a data integration tool for constructing a hydraulic conductivity map by integrating geophysical data and hydraulic conductivity measurements. The developed code employed a geostatistical optimization method, simulated annealing (SA), and consists of 4 distinct computation modules by which from exploratory data analysis to postprocessing of the simulation were processed. All these modules are equipped with Graphical User Interface (GUI). Validation of the developed code was evaluated in-situ in characterizing hydraulic characteristics of highly permeable fractured zone.