• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.309-315
• /
• 2009

Integrated interpretation of multi-geophysical data has been continuously used in terms that it has provided more confident information than the result from single-geophysical data. Especially, geostatistical integration has its own superiority that it is possible to deal with spatial characteristics as well as physical properties of survey data and the process of integration is clear. This paper further extends the previous work of geostatistical inversion for integrated interpretation. In this paper, we propose a new way of dealing with the case that the multi-geophysical data do not share the measurement range. According to the geostatistical kriging, the closer between the measurement points, the smaller kriging variance we get, and vice versa. We used this spatial properties as a weighting value to the process of geostatistical inversion for the geophysical data integration. An objective way to integrate different kinds of geophysical data measured at different ranges is provided with this algorithm.

• Korean Journal of Remote Sensing
• /
• v.24 no.5
• /
• pp.453-462
• /
• 2008

The objective of this paper is to compare the prediction performances of different landslide hazard maps based on topographic data stemming from different sources of elevation. The geostatistical framework of kriging, which can properly integrate spatial data with different accuracy, is applied for generating more reliable elevation estimates from both sparse elevation spot heights and exhaustive ASTER-based elevation values. A case study from Boeun, Korea illustrates that the integration of elevation and slope maps derived from different data yielded different prediction performances for landslide hazard mapping. The landslide hazard map constructed by using the elevation and the associated slope maps based on geostatistical integration of spot heights and ASTER-based elevation resulted in the best prediction performance. Landslide hazard mapping using elevation and slope maps derived from the interpolation of only sparse spot heights showed the worst prediction performance.

• Geomechanics and Engineering
• /
• v.34 no.1
• /
• pp.29-41
• /
• 2023

The 3D geospatial modeling of geotechnical information can aid in understanding the geotechnical characteristic values of the continuous subsurface at construction sites. In this study, a geostatistical optimization model for the three-dimensional (3D) mapping of subsurface stratification and the SPT-N value based on a trial-and-error rule was developed and applied to a dam emergency spillway site in South Korea. Geospatial database development for a geotechnical investigation, reconstitution of the target grid volume, and detection of outliers in the borehole dataset were implemented prior to the 3D modeling. For the site-specific subsurface stratification of the engineering geo-layer, we developed an integration method for the borehole and geophysical survey datasets based on the geostatistical optimization procedure of ordinary kriging and sequential Gaussian simulation (SGS) by comparing their cross-validation-based prediction residuals. We also developed an optimization technique based on SGS for estimating the 3D geometry of the SPT-N value. This method involves quantitatively testing the reliability of SGS and selecting the realizations with a high estimation accuracy. Boring tests were performed for validation, and the proposed method yielded more accurate prediction results and reproduced the spatial distribution of geotechnical information more effectively than the conventional geostatistical approach.

• Proceedings of the KSRS Conference
• /
• v.1
• /
• pp.406-408
• /
• 2006

Various geological thematic maps such as grain size or ground water level maps have been generated by interpolating sparsely sampled ground survey data. When there are sampled data at a limited number of locations, to use secondary information which is correlated to primary variable can help us to estimate the attribute values of the primary variable at unsampled locations. This paper applies two multivariate geostatistical algorithms to integrate remote sensing imagery with sparsely sampled ground survey data for spatial estimation of grain size: simple kriging with local means and kriging with an external drift. High-resolution IKONOS imagery which is well correlated with the grain size is used as secondary information. The algorithms are evaluated from a case study with grain size observations measured at 53 locations in the Baramarae beach of Anmyeondo, Korea. Cross validation based on a one-leave-out approach is used to compare the estimation performance of the two multivariate geostatistical algorithms with that of traditional ordinary kriging.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2007.06a
• /
• pp.42-47
• /
• 2007

We have studied feasibility of the geostatistical approach to enhance the result of analysis of the sparsely obtained MT(Magnetotelluric) data by combining with gravity data. We have attempted to use geostatistics for integrating the MT data along with gravity data. To evaluate the feasibility of this approach, we have studied about interrelation between geological boundary and density distribution, and corrected density distribution for conversion to more sensitive to geological boundary by minimization of difference between z-directional variogram values of resistivity distribution obtained MT inversion and density distributions. Then, this method has been tested on model and field data. In model test, the results obtained were good agreement with real model. And in a real field data, the result of analysis demonstrate convincingly that our geostatistical approach is effective.

• Korean Journal of Remote Sensing
• /
• v.29 no.1
• /
• pp.69-79
• /
• 2013

This paper presents a two-stage geostatistical integration approach that aims at downscaling of coarse scale remote sensing data. First, downscaling of the coarse scale sedoncary data is implemented using area-to-point kriging, and this result will be used as trend components on the next integration stage. Then simple kriging with local varying means that integrates sparse precise observation data with the downscaled data is applied to generate thematic information at a finer scale. The presented approach can not only account for the statistical relationships between precise observation and secondary data acquired at the different scales, but also to calibrate the errors in the secondary data through the integration with precise observation data. An experiment for precipitation mapping with weather station data and TRMM (Tropical Rainfall Measuring Mission) data acquired at a coarse scale is carried out to illustrate the applicability of the presented approach. From the experiment, the geostatistical downscaling approach applied in this paper could generate detailed thematic information at various finer target scales that reproduced the original TRMM precipitation values when upscaled. And the integration of the downscaled secondary information with precise observation data showed better prediction capability than that of a conventional univariate kriging algorithm. Thus, it is expected that the presented approach would be effectively used for downscaling of coarse scale data with various data acquired at different scales.

• Geophysics and Geophysical Exploration
• /
• v.7 no.2
• /
• pp.121-129
• /
• 2004

The geostatistical approach was applied to integrate MT (Magneto-telluric) resistivity data and borehole information for the spatial RMR (Rock Mass Rating) evaluation. Generally, resistivity of the subsurface is believed to be positively related to the RMR, thus the resistivity and borehole RMR information was combined in a geostatistical approach. To relate the two different sets of data, we take the MT resistivity data as secondary information and estimate the RMR mean values at unsampled points by identification of the resistivity to the borehole data. Two types of approach are performed for the estimation of RMR mean values. Then the residuals of the RMR values around the borehole sites are geostatistically modeled to infer the spatial structure of difference between real RMR values and estimated mean values. Finally, this geostatistical estimation is added to the previous means. The result applied to a real situation shows prominent improvements to reflect the subsurface structure and spatial resolution of RMR information.

• Korean Journal of Remote Sensing
• /
• v.22 no.6
• /
• pp.581-593
• /
• 2006

Various geological thematic maps have been generated by interpolating sparsely sampled ground survey data and geostatistical kriging that can consider spatial correlation between neighboring data has widely been used. This paper applies multi-variate geostatistical algorithms to integrate secondary information with sparsely sampled ground survey data for geological thematic mapping. Simple kriging with local means and kriging with an external drift are applied among several multi-variate geostatistical algorithms. Two case studies for spatial mapping of groundwater level and grain size have been carried out to illustrate the effectiveness of multi-variate geostatistical algorithms. A digital elevation model and IKONOS remote sensing imagery were used as secondary information in two case studies. Two multi-variate geostatistical algorithms, which can account for both spatial correlation of neighboring data and secondary data, showed smaller prediction errors and more local variations than those of ordinary kriging and linear regression. The benefit of applying the multi-variate geostatistical algorithms, however, depends on sampling density, magnitudes of correlation between primary and secondary data, and spatial correlation of primary data. As a result, the experiment for spatial mapping of grain size in which the effects of those factors were dominant showed that the effect of using the secondary data was relatively small than the experiment for spatial mapping of groundwater level.

• Geophysics and Geophysical Exploration
• /
• v.15 no.4
• /
• pp.177-189
• /
• 2012

Evaluation of three-dimensional ore body modeling has been performed by applying the geostatistical integration technique to multiple geophysical (electrical resistivity, MT) and geological (borehole data, physical properties of core) information. It was available to analyze the resistivity range in borehole and other area through multiple geophysical data. A correlation between resistivity and density from physical properties test of core was also analyzed. In the case study results, the resistivity value of ore body is decreased contrast to increase of the density, which seems to be related to a reason that the ore body (magnetite) includes heavy conductive component (Fe) in itself. Based on the lab test of physical properties in iron mine region, various geophysical, geological and borehole data were used to provide ore body modeling, that is electrical resistivity, MT, physical properties data, borehole data and grade data obtained from borehole data. Of the various geostatistical techniques for the integrated data analysis, in this study, the SGS (sequential Gaussian simulation) method was applied to describe the varying non-homogeneity depending on region through the realization that maintains the mean and variance. With the geostatistical simulation results of geophysical, geological and grade data, the location of residual ore body and ore body which is previously reported was confirmed. In addition, another highly probable region of iron ore bodies was estimated deeper depth in study area through integrated modeling.

• Journal of the Korean earth science society
• /
• v.26 no.5
• /
• pp.429-435
• /
• 2005

3D RMR (Rock Mass Rating) analysis has been performed by applying the Geostatistical integration technique for geophysical and borehole data. Of the various geostatistical techniques for the integrated data analysis, in this study, we applied the SKlvm (Simple Kriging with local varying means) method that substitutes the values of the interpreted geophysical result with the mean values of the RMR at the location to be inferred. The substitution is performed by the indicator transform between the result of geophysical interpretation and the observed RMR values at borehole sites. The used geophysical data are the electrical resistivity and MT result, and 10 borehole sites are investigated to obtain the RMR values. This integrated analysis makes the interpretation to be more practical for identifying the realistic RMR distribution that supports the regional geological situation.