• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.350-350
• /
• 2017

저수지의 설계홍수량 산정 시 인근의 기상관측 자료를 활용하고 있으나 인근에 기상관측 자료가 없거나 저수지 배후 유역이 큰 경우에는 단일 기상관측 자료를 이용하기에는 한계가 있다. 따라서 실무적으로 지점별 기상관측소의 자료를 이용하여 설계홍수량을 산정할 때에는 각 관측소 자료를 이용하여 확률강우량을 산정하고 Thiessen 가중평균을 한 후 면적우량환산계수 (ARF)를 곱하여 사용하고 있는데, Thiessen 방법의 경우 방법이 간단하지만 지형 고도 효과는 무시되고 우량계의 지배면적에 의한 우량계의 분포 상태만을 고려하게 된다. 그러므로 설계홍수량 산정시 사용되는 Thiessen 방법은 공간적 불확실성을 내포하고 있고, 특히 소규모 저수지의 설계홍수량을 산정하는 경우에는 저수지 유역의 국소적인 특징을 나타내기 어렵다. 본 연구에서는 설계홍수량 산정 시 저수지 위치에 해당하는 확률강우량의 공간적 불확실성을 평가하기 위하여 SIS(Sequential Indicator Simulation) 방법을 이용하였다. SIS 방법은 Kriging 기법과 마찬가지로 베리오그램으로부터 얻어지는 공간적 상관관계를 기반으로 하고 있는 방법으로 Kriging 기법과 달리 공간분포의 국소적인 특성을 평가할 수 있다는 장점을 가지고 있다.

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

• Tunnel and Underground Space
• /
• v.16 no.2 s.61
• /
• pp.135-145
• /
• 2006

In this study, conditional simulation was conducted to estimate rock mass rating(RMR) in unsurveyed regions. Sequential Gaussian simulation(SGS) and sequential indicator simulation(SIS) were applied for estimating RMR from the bore hole logging data. The uncertainty of SGS and SIS was verified by sample cross validation. A subset composed of 5 bore hole logging data among the original 30 bore hole logging data was set aside as test data. The remainder was training data. The quality of SGS and SIS estimation on the testing data reflects how well it would perform in an unsupervised setting. SGS and SIS were useful stochastic methods to estimate the spatial distribution of rock mass classes correctly and assess the uncertainty of estimation quantitatively. The result of conditional simulation can offer useful information of rock mass classes such as RMR in unsurveyed regions.

• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.257-267
• /
• 2013

3-D ore modeling was performed to understand the configuration of ore bodies by integrated analysis of borehole and geophysical data in iron-mine area. Five representative indices of rocks were designated, which were obtained from geological survey and borehole. The five indices of rocks were geostatistically simulated by Sequential Indicator Simulation method to delineate boundary of the ore bodies. And Ordinary Kriging and Sequential Gaussian Simulation was applied to make secondary information using resistivity data from magnetotellurics and DC resistivity survey, and this information was used for simple kriging with local varying means, one of integrated kriging techniques. From the correlation analysis between each properties, it was found that high grade of ore is characterized by increased density, whereas the electrical resistivity decreases. With the integrated results of geophysical and borehole data, it was also found that the real configuration of ore body was similar to the modeled result and information about ore grade in 3-D space was obtained.

• Journal of the Korean earth science society
• /
• v.23 no.1
• /
• pp.15-29
• /
• 2002

Bayesian inversion for gravity and resistivity data was performed to investigate the cavity structure appearing as a lava tunnel in Cheju Island, Korea. Dipole-dipole DC resistivity data were proposed for a prior information of gravity data and we applied the geostatistical techniques such as kriging and simulation algorithms to provide a prior model information and covariance matrix in data domain. The inverted resistivity section gave the indicator variogram modeling for each threshold and it provided spatial uncertainty to give a prior PDF by sequential indicator simulations. We also presented a more objective way to make data covariance matrix that reflects the state of the achieved field data by geostatistical technique, cross-validation. Then Gaussian approximation was adopted for the inference of characteristics of the marginal distributions of model parameters and Broyden update for simple calculation of sensitivity matrix and SVD was applied. Generally cavity investigation by geophysical exploration is difficult and success is hard to be achieved. However, this exotic multiple interpretations showed remarkable improvement and stability for interpretation when compared to data-fit alone results, and suggested the possibility of diverse application for Bayesian inversion in geophysical inverse problem.

• Journal of Energy Engineering
• /
• v.22 no.2
• /
• pp.156-168
• /
• 2013

Drilling expedition #1 in 2007 and drilling expedition #2 in 2010 were performed for gas hydrate resources evaluation and optimal site selection of pilot test in Ulleung basin, East Sea, Korea. This study presents to build the 3D spatial distribution models using the estimated sedimentary facies, porosity, and gas hydrate saturation derived by well logs and core analysis data from UBGH1-4, UBGH1-9, UBGH1-10, UBGH1-14, UBGH2-2-1, UBGH2-2-2, UBGH2-6, UBGH2-9, UBGH2-10 and UBGH2-11. The objective of 3D spatial distribution modeling is to build a geological representation of the gas hydrate-bearing sediment that honors the heterogeneity in 3D grid scale. The facies modeling is populating sedimentary facies into a geological grid using sequential indicator simulation. The porosity and gas hydrate saturation modeling used sequential Gaussian simulation to populate properties stochastically into grid cells.

• Tunnel and Underground Space
• /
• v.15 no.4 s.57
• /
• pp.264-274
• /
• 2005

The optimum blasting pattern to excavate a quarry efficiently and economically can be determined based on the minimum production cost which is generally estimated according to rock fragmentation. Therefore it is a critical problem to predict fragment size distribution of blasted rocks over an entire quarry. By comparing various prediction models, it can be ascertained that the result obtained from Kuz-Ram model relatively coincides with that of field measurements. Kuz-Ram model uses the concept of rock factor to signify conditions of rock mass such as block size, rock jointing, strength and others. For the evaluation of total production cost, it is imperative to estimate 3-D spatial distribution of rock factor for the entire quarry. In this study, a sequential indicator simulation technique is adopted for estimation of spatial distribution of rock factor due to its higher reproducibility of spatial variability and distribution models than Kriging methods. Further, this can reduce the uncertainty of predictor using distribution information of sample data The entire quarry is classified into three types of rock mass and optimum blasting pattern is proposed for each type based on 3-D spatial distribution of rock factor. In addition, plane maps of rock factor distribution for each ground levels is provided to estimate production costs for each process and to make a plan for an optimum blasting pattern.

• Journal of Soil and Groundwater Environment
• /
• v.19 no.4
• /
• pp.31-44
• /
• 2014

In the present study, a few of recently developed geostatistical models are comparatively studied. The models are two-point statistics based sequential indicator simulation (SISIM) and generalized coupled Markov chain (GCMC), multi-point statistics single normal equation simulation (SNESIM), and object based model of FLUVSIM (fluvial simulation) that predicts structures of target object from the provided geometric information. Out of the models, SNESIM and FLUVSIM require additional information other than conditioning data such as training map and geometry, respectively, which generally claim demanding additional resources. For the comparative studies, three-dimensional fluvial reservoir model is developed considering the genetic information and the samples, as input data for the models, are acquired by mimicking realistic sampling (i.e. random sampling). For SNESIM and FLUVSIM, additional training map and the geometry data are synthesized based on the same information used for the objective model. For the comparisons of the predictabilities of the models, two different measures are employed. In the first measure, the ensemble probability maps of the models are developed from multiple realizations, which are compared in depth to the objective model. In the second measure, the developed realizations are converted to hydrogeologic properties and the groundwater flow simulation results are compared to that of the objective model. From the comparisons, it is found that the predictability of GCMC outperforms the other models in terms of the first measure. On the other hand, in terms of the second measure, the both predictabilities of GCMC and SNESIM are outstanding out of the considered models. The excellences of GCMC model in the comparisons may attribute to the incorporations of directional non-stationarity and the non-linear prediction structure. From the results, it is concluded that the various geostatistical models need to be comprehensively considered and comparatively analyzed for appropriate characterizations.

• 지반과기술
• /
• v.3 no.3
• /
• pp.15-22
• /
• 2006

The optimum blasting pattern to excavate a quarry efficiently and economically can be determined based on the minimum production cost, which is generally estimated according to rock fragmentation. Therefore, it is a critical problem to predict fragment size distribution of blasted rocks over an entire quarry. By comparing various prediction models, it can be ascertained that the result obtained from Kuz-Ram model relatively coincides with that of field measurements. Kuz-Ram model uses the concept of rock factor to signify conditions of rock mass such as block size, rock jointing, strength and others. For the evaluation of total production cost, it is imperative to estimate 3-D spatial distribution of rock factor for the entire quarry. In this study, a sequential indicator simulation technique is adopted for estimation of spatial distribution of rock factor due to its higher reproducibility of spatial variability and distribution models than Kriging methods. Further, this can reduce the uncertainty of predictor using distribution information of sample data. The entire quarry is classified into three types of rock mass and optimum blasting pattern is proposed for each type based on 3-D spatial distribution of rock factor. In addition, plane maps of rock factor distribution for each ground level are provided to estimate production costs for each process and to make a plan for an optimum blasting pattern.

• Journal of Soil and Groundwater Environment
• /
• v.10 no.5
• /
• pp.52-60
• /
• 2005

The conceptual model of under-sampled study area will include a great amount of uncertainty. In this study, we investigate the applicability of Markov chain model in a spatial domain as a tool for minimizing the uncertainty arose from the lack of data. A new formulation is developed to generalize the previous two-dimensional coupled Markov chain model, which has more versatility to fit any computational sequence. Furthermore, the computational algorithm is improved to utilize more conditioning information and reduce the artifacts, such as the artificial parcel inclination, caused by sequential computation. A generalized 20 coupled Markov chain (GCMC) is tested through applying a hypothetical soil map to evaluate the appropriateness as a substituting model for conventional geostatistical models. Comparing to sequential indicator model (SIS), the simulation results from GCMC shows lower entropy at the boundaries of indicators which is closer to real soil maps. For under-sampled indicators, however, GCMC under-estimates the presence of the indicators, which is a common aspect of all other geostatistical models. To improve this under-estimation, further study on data fusion (or assimilation) inclusion in the GCMC is required.