• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.127-127
• /
• 2021

In order to estimate parameter uncertainty of hydrological models, the consideration of the likelihood functions which provide reliable parameters of model is necessary. In this study, the Bayesian Markov Chain Monte Carlo (MCMC) method with informal likelihood functions is used to analyze the uncertainty of parameters of the SURR model for estimating the hourly streamflow of Gunnam station of Imjin basin, Korea. Three events were used to calibrate and one event was used to validate the posterior distributions of parameters. Moreover, the performance of four informal likelihood functions (Nash-Sutcliffe efficiency, Normalized absolute error, Index of agreement, and Chiew-McMahon efficiency) on uncertainty of parameter is assessed. The indicators used to assess the uncertainty of the streamflow simulation were P-factor (percentage of observed streamflow included in the uncertainty interval) and R-factor (the average width of the uncertainty interval). The results showed that the sensitivities of parameters strongly depend on the likelihood functions and vary for different likelihood functions. The uncertainty bounds illustrated the slight differences from various likelihood functions. This study confirms the importance of the likelihood function selection in the application of Bayesian MCMC to the uncertainty assessment of the SURR model.

• Communications for Statistical Applications and Methods
• /
• v.28 no.5
• /
• pp.425-445
• /
• 2021

A generalization of the log-logistic (LL) distribution called exponentiated log-logistic (ELL) distribution on lines of exponentiated Weibull distribution is considered. In this paper, based on progressive type-II censored samples, we have derived the maximum likelihood estimators and Bayes estimators for three parameters, the survival function and hazard function of the ELL distribution. Then, under the balanced squared error loss (BSEL) and the balanced linex loss (BLEL) functions, their corresponding Bayes estimators are obtained using Lindley's approximation (see Jung and Chung, 2018; Lindley, 1980), Tierney-Kadane approximation (see Tierney and Kadane, 1986) and Markov Chain Monte Carlo methods (see Hastings, 1970; Gelfand and Smith, 1990). Here, to check the convergence of MCMC chains, the Gelman and Rubin diagnostic (see Gelman and Rubin, 1992; Brooks and Gelman, 1997) was used. On the basis of their risks, the performances of their Bayes estimators are compared with maximum likelihood estimators in the simulation studies. In this paper, research supports the conclusion that ELL distribution is an efficient distribution to modeling data in the analysis of survival data. On top of that, Bayes estimators under various loss functions are useful for many estimation problems.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.262-271
• /
• 2002

This study presents a practical procedure for the Bayesian inversion of geophysical data. We have applied geostatistical techniques for the acquisition of prior model information, then the Markov Chain Monte Carlo (MCMC) method was adopted to infer the characteristics of the marginal distributions of model parameters. For the Bayesian inversion of dipole-dipole array resistivity data, we have used the indicator kriging and simulation techniques to generate cumulative density functions from Schlumberger array resistivity data and well logging data, and obtained prior information by cokriging and simulations from covariogram models. The indicator approach makes it possible to incorporate non-parametric information into the probabilistic density function. We have also adopted the MCMC approach, based on Gibbs sampling, to examine the characteristics of a posteriori probability density function and the marginal distribution of each parameter.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.33 no.10
• /
• pp.1163-1170
• /
• 2009

Reliability analysis is of great importance in the advanced product design, which is to evaluate reliability due to the associated uncertainties. There are three types of uncertainties: the first is the aleatory uncertainty which is related with inherent physical randomness that is completely described by a suitable probability model. The second is the epistemic uncertainty, which results from the lack of knowledge due to the insufficient data. These two uncertainties are encountered in the input variables such as dimensional tolerances, material properties and loading conditions. The third is the metamodel uncertainty which arises from the approximation of the response function. In this study, an integrated method for the reliability analysis is proposed that can address all these uncertainties in a single Bayesian framework. Markov Chain Monte Carlo (MCMC) method is employed to facilitate the simulation of the posterior distribution. Mathematical and engineering examples are used to demonstrate the proposed method.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2003.11a
• /
• pp.340-343
• /
• 2003

This study presents a practical procedure for the Bayesian inversion of geophysical data by Markov chain Monte Carlo (MCMC) sampling and geostatistics. We have applied geostatistical techniques for the acquisition of prior model information, and then the MCMC method was adopted to infer the characteristics of the marginal distributions of model parameters. For the Bayesian inversion of dipole-dipole array resistivity data, we have used the indicator kriging and simulation techniques to generate cumulative density functions from Schlumberger array resistivity data and well logging data, and obtained prior information by cokriging and simulations from covariogram models. The indicator approach makes it possible to incorporate non-parametric information into the probabilistic density function. We have also adopted the MCMC approach, based on Gibbs sampling, to examine the characteristics of a posteriori probability density function and the marginal distribution of each parameter. This approach provides an effective way to treat Bayesian inversion of geophysical data and reduce the non-uniqueness by incorporating various prior information.

• Smart Structures and Systems
• /
• v.21 no.5
• /
• pp.601-609
• /
• 2018

The estimated probabilistic model of wind data based on the conventional approach may have high discrepancy compared with the true distribution because of the uncertainty caused by the instrument error and limited monitoring data. A sequential quadratic programming (SQP) algorithm-based finite mixture modeling method has been developed in the companion paper and is conducted to formulate the joint probability density function (PDF) of wind speed and direction using the wind monitoring data of the investigated bridge. The established bivariate model of wind speed and direction only represents the features of available wind monitoring data. To characterize the stochastic properties of the wind parameters with the subsequent wind monitoring data, in this study, Bayesian inference approach considering the uncertainty is proposed to update the wind parameters in the bivariate probabilistic model. The slice sampling algorithm of Markov chain Monte Carlo (MCMC) method is applied to establish the multi-dimensional and complex posterior distribution which is analytically intractable. The numerical simulation examples for univariate and bivariate models are carried out to verify the effectiveness of the proposed method. In addition, the proposed Bayesian inference approach is used to update and optimize the parameters in the bivariate model using the wind monitoring data from the investigated bridge. The results indicate that the proposed Bayesian inference approach is feasible and can be employed to predict the bivariate distribution of wind speed and direction with limited monitoring data.

• Journal of Korea Water Resources Association
• /
• v.50 no.4
• /
• pp.277-288
• /
• 2017

A probabilistic risk analysis of levee system estimates the overall level of flood risk associated with the levee system, according to a series of possible flood scenarios. It requires the uncertainty analysis of all the risk components, including hydrological, hydraulic and geotechnical parts computed by employing MCMC (Markov Chain Monte Carlo), MCS (Monte Carlo Simulation) and FOSM (First-Order Second Moment), presents a joint probability combined each probability. The methodology was applied to a 12.5 km reach from upstream to downstream of the Gangjeong-Goryeong weir, including 6 levee reaches, in Nakdong river. Overtopping risks were estimated by computing flood stage corresponding to 100/200 year high quantile (97.5%) design flood causing levee overflow. Geotechnical risks were evaluated by considering seepage, slope stability, and rapid drawdown along the levee reach without overflow. A probability-based compound risk will contribute to rising effect of safety and economic aspects for levee design, then expect to use the index for riverside structure design in the future.

• International Journal of Reliability and Applications
• /
• v.13 no.1
• /
• pp.37-47
• /
• 2012

In this article we attempted reliability analysis of a component under the stress-strength pattern with both classical as well as Bayesian techniques. The main focus is made to develop the theory for dealing the reliability problems in various circumstances for bivariate environmental set up in context of Bayesian paradigm. A stress-strength based model describes the life of a component which has strength (Y) and is subjected to stress(X). We develop the Bayes and moment estimators of reliability of a component for each of the three possible conditions, under the assumption that the two stresses (i.e. $X_1$ and $X_2$) on a component are dependent and follow a Bivariate exponential (BVE) of Marshall-Olkin distribution, the strength of a component (Y) following exponential distribution is independent of the stresses. The simulation study is performed with Markov Chain Monte Carlo technique via Gibbs sampler to obtain the estimates of Bayes estimators of reliability, are compared with moment estimators of reliabilities on the basis of absolute biases.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.2
• /
• pp.125-132
• /
• 2019

This paper presents a seismic reliability analysis method for an offshore wind turbine with a twisted tripod support structure under earthquake loading. A three dimensional dynamic finite element model is proposed to consider the nonlinearity of the ground-pile interactions and the geometrical characteristics of the twisted tripod support structure where out-of-plane displacement occurs even under in-plane lateral loadings. For the evaluation of seismic reliability, the failure probability was calculated for the maximum horizontal displacement of the pile head, which is calculated from time history analysis using artificial earthquakes for the design return periods. The application of the subset simulation method using the Markov Chain Monte Carlo(MCMC) sampling is proposed for efficient reliability analysis considering the limit state equation evaluation by the nonlinear time history analysis. The proposed method can be applied to the reliability evaluation and design criteria development of the offshore wind turbine with twisted tripod support structure in which two dimensional models and static analysis can not produce accurate results.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.6
• /
• pp.394-404
• /
• 2013

In the present study, the characteristics of spectral peakedness parameter $Q_p$, bandwidth parameter ${\varepsilon}$, and spectral width parameter ${\nu}$ were analyzed as a first step to define the swell waves quantitatively. For the analysis, the joint probability density function of significant wave heights and peak periods were newly developed. The MCMC(Markov Chain Monte Carlo) simulations have been performed to generate the significant wave heights and peak periods from the developed probability density functions. Applying the simulated significant wave heights and peak periods to the theoretical wave spectrum models, the spectral shapes parameters were obtained and analyzed. Among the spectral shape parameters, only the spectral peakedness parameter $Q_p$, is shown to be independent with the significant wave height and peak wave period. It also best represents the peakedness of the spectral shape, and henceforth $Q_p$ should be used to define the swell waves with a wave period. For the field verification of the results, wave data obtained from Hupo port and Ulleungdo were analyzed and results showed the same trend with the MCMC simulation results.