• Journal of the Korean earth science society
• /
• v.31 no.7
• /
• pp.681-690
• /
• 2010

Inversion result of geophysical data given as a block type was geostatistically simulated with borehole observation given as a point type and was applied to the rock classifying map. The geophysical data generally involved secondary information for the target material and were obtained for overall region. In contrast, borehole data provided direct information for the target material, but tended to be effective only for a narrow range of region and were dealt as a point type. Integrated simulation or kriging interpolation of these two different kinds of information required the covariance for point-point, point-block and block-block. Using the Bssim module included in SGeMS software, integrated result of geophysical data and borehole data were obtained. The results were then compared with the method of geostatistical inversion proposed by authors. Downscaling method used in this study showed relatively more flexible than the geostatistical inversion.

• Journal of Soil and Groundwater Environment
• /
• v.17 no.5
• /
• pp.56-67
• /
• 2012

In the present study, an enhanced subsurface prediction algorithm based on a non-parametric geostatistical model and a history matching technique through Gibbs sampler is developed and the iterative prediction improvement procedure is proposed. The developed model is applied to a simple two-dimensional synthetic case where domain is composed of three different hydrogeologic media with $500m{\times}40m$ scale. In the application, it is assumed that there are 4 independent pumping tests performed at different vertical interval and the history curves are acquired through numerical modeling. With two hypothetical borehole information and pumping test data, the proposed prediction model is applied iteratively and continuous improvements of the predictions with reduced uncertainties of the media distribution are observed. From the results and the qualitative/quantitative analysis, it is concluded that the proposed model is good for the subsurface prediction improvements where the history data is available as a supportive information. Once the proposed model be a matured technique, it is believed that the model can be applied to many groundwater, geothermal, gas and oil problems with conventional fluid flow simulators. However, the overall development is still in its preliminary step and further considerations needs to be incorporated to be a viable and practical prediction technique including multi-dimensional verifications, global optimization, etc. which have not been resolved in the present study.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2000.11a
• /
• pp.15-22
• /
• 2000

Long ranged groundwater-level data often have the missing intervals because of the trouble of monitoring systems at the national groundwater monitoring wells. Geostatistical methods are very useful for the supplement of the missing data. Ordinary kriging was applied for the interpolation of the missing groundwater-level data with a smooth sinusoidal variation. Conditional simulation was used for the reproduction of the missing data with high fluctuations. Two geostatistical methods produced the very accurate estimates at the missing intervals and reproduced their original variations. This fact is proved by the cross validation test and graphical method, respectively.

• Journal of The Geomorphological Association of Korea
• /
• v.20 no.4
• /
• pp.85-99
• /
• 2013

The objective of this paper is to analyze the spatial distribution of errors in the DEM generated using waterlines from multi-temporal remote sensing data and to assess flood vulnerability. Unlike conventional research in which only global statistics of errors have been generated, this paper tries to quantitatively analyze the spatial distribution of errors from a probabilistic viewpoint using geostatistical simulation. The initial DEM in Baramarae tidal flats was generated by corrected tidal level values and waterlines extracted from multi-temporal Landsat data in 2010s. When compared with the ground measurement height data, overall the waterline-based DEM underestimated the actual heights and local variations of the errors were observed. By applying sequential Gaussian simulation based on spatial autocorrelation of DEM errors, multiple alternative error distributions were generated. After correcting errors in the initial DEM with simulated error distributions, probabilities for flood vulnerability were estimated under the sea level rise scenarios of IPCC SERS. The error analysis methodology based on geostatistical simulation could model both uncertainties of the error assessment and error propagation problems in a probabilistic framework. Therefore, it is expected that the error analysis methodology applied in this paper will be effectively used for the probabilistic assessment of errors included in various thematic maps as well as the error assessment of waterline-based DEMs in tidal flats.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.63-67
• /
• 2008

Geostatistical inverse approach using geophysical data was applied to indirectly make the RMR classification at points apart from boreholes. The geostatistical appoach was usually used to find optimized estimation which supports two or more different physical properties at unsampled points. However, in this study, an approach to solve inverse problem was proposed. The primary variable, RMR values obtained at known boreholes, is geostatistically simulated with many realization at pre-defined grid point according to the variogram model. The simulated values are sequentially compared with the physical property resulted from geophysical survey at an arbitrary grid point, and the most similar one is chosen. This process means that the spatial distribution of primary variable, RMR, is conformed well to the original pattern of the borehole observation, and ensure to fit the geophysical survey result to reflect the correlation between different physical properties.

• Korean Journal of Remote Sensing
• /
• v.26 no.1
• /
• pp.29-38
• /
• 2010

It is necessary to develop an integration model which can account for various data acquired at different measurement scales in environmental thematic mapping with high-resolution ground survey data and relatively low-resolution remote sensing data. This paper presents and applies a multi-scale geostatistical methodology for downscaling of thematic maps generated from lowresolution remote sensing data. This methodology extends a traditional ordinary kriging system to a block kriging system which can account for both ground data and remote sensing data which can be regarded as point and block data, respectively. In addition, stochastic simulation based on block kriging is also applied to describe spatial uncertainty attached to the downscaling. Two downscaling experiments including SRTM DEM and MODIS Leaf Area Index (LAI) products were carried out to illustrate the applicability of the geostatistical methodology. Through the experiments, multi-scale geostatistics based on block kriging successfully generated relatively high-resolution thematic maps with reliable accuracy. Especially, it is expected that multiple realizations generated from simulation would be effectively used as input data for investigating the effects of uncertain input data on GIS model outputs.

• Geophysics and Geophysical Exploration
• /
• v.19 no.4
• /
• pp.187-199
• /
• 2016

In order to delineate the soft ground distribution, an integrated geostatistical analysis was performed using the MASW (Multichannel Analysis of Surface Wave) which has the information of overall region and CPTu (Piezo Cone Penetration Test) which provides the direct information of the measuring point of the ground. MASW results were known to have close relationship with the ground stiffness. This correlation was confirmed through the comparison of MASW data obtained from two survey lines to the laboratory test with extracted soil samples. 3D physical property distribution in the study area was acquired by geostatistical integrated analysis with the data of tip resistance ($q_c$) and pore pressure (u) from the CPTu obtained at 6 points within the study area. The integrated analysis was conducted by applying the COSGSIM (Sequential Gaussian Co-Simulation) technology which can carry out the simulation in accordance with the spatial correlation between the MASW results and both tip resistance and pore pressure. Besides the locations of CPTu, borehole investigations were also conducted at two different positions. As a result, the N value of SPT and borehole log could be secured, so these data were used for the analysis of the geotechnical engineering accuracy of the integrated analysis result. For the verification of reliability of the 3D distribution of tip resistance and pore pressure secured through integrated analysis, the geotechnical information gained from the two drilling areas was compared, and the result showed extremely high correlation.

• Proceedings of the Korean Association of Geographic Inforamtion Studies Conference
• /
• 2004.03a
• /
• pp.463-468
• /
• 2004

이 연구에서는 원격탐사 자료를 이용하여 얻어진 분류 결과로부터 분류 정확도의 공간적 불확실성을 추정하고자 하였다. 기존 분류결과로부터 얻어지는 토지 피복별 확률값을 지구통계학적 시뮬레이션 기법을 이용하여 참조자료의 공간적 분포와 통합하였다. 다중센서 화상 통합에 적용한 결과, 각 토지피복 항목별로 공간적인 정확도 분포를 얻을 수 있을 수 있었으며 이러한 자료는 분류결과를 해석하는데 유용하게 사용될 수 있을 것으로 기대된다.

• Economic and Environmental Geology
• /
• v.30 no.5
• /
• pp.493-498
• /
• 1997

The purpose of this study is to develop geostatistical model for evaluating the abundance of deep-sea manganese nodule. The abundance data used in this study were obtained from the KODOS (Korea Deep Ocean Study) area. The variation of nodule abundance was very high within short distance, while sampling methods was very limited. As the distribution of nodule abundance showed non-gaussian, indicator simulation method was used instead of conditional simulation method and/or ordinary kriging. The abundance data were encoded into a series of indicators with 6 cutoff values. They were used to estimate the conditional probability distribution function (cpdf) of the nodule abundance at any unsampled location. The standardized indicator variogram models were obtained according to variogram analysis. This SIS method had the advantage over other traditional techniques such as the turning bands method and ordinary kriging. The estimating values by indicator conditional simulation near high abundance area were more detailed than by ordinary kriging and indicator kriging. They also showed better spatial characteristics of distribution of nodule abundance.

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