Copyright © 2006 by The McGraw-Hill Companies, Inc. Click here for terms of use

As noted previously, the accuracy of the model output statistics and probabil­ity distribution (e. g., probability that a specified safety level will be exceeded) obtained from Monte Carlo simulation is a function of the number of simu­lations performed. For models or problems with a large number of uncertain basic variables and for which low probabilities (<0.1) are of interest, tens of thousands of simulations may be required. Rules for determining the number of simulations required for convergence are not available, and thus replication of the Monte Carlo simulation runs for a given number of simulations is the only way to check convergence (Melching, 1995). Cheng et al. (1982) considered the convergence characteristics of Monte Carlo simulation for a simple case of Z = X3X4 — (X1 + X2), where the distributions and statistics of the variables are listed in Table 6.1. They found that failure probabilities (i. e., probability of Z < 0) down to 0.0025 could be estimated reliably with 32,000 simulations, failure probabilities down to 0.015 could be estimated reliably with 8000 sim­ulations, and failure probabilities down to 0.2 could be estimated reliably with 1000 simulations.

Problems involving more complex system functions Z and more basic vari­ables may require more simulations to obtain similar accuracy. For example, Melching (1992) found that 1000 simulations were adequate to estimate the mean, standard deviation, and quantiles above 0.2 for an application of the HEC-1 (U. S. Army Corps of Engineers, 1990) and RORB (Laurenson and Mein, 1985) rainfall-runoff models and that 10,000 simulations were needed to ac­curately estimate quantiles between 0.01 and 0.2. Brown and Barnwell (1987) reported that for the QUAL2E multiple-constituent (dissolved oxygen, nitrogen cycle, algae, etc.) steady-state surface water-quality model, 2000 simulations were required to obtain accurate estimates of the output standard deviation. With the computational speed of today’s computers, making even 10,000 runs is not prohibitive for simpler models. However, increased computational speed has made possible the use of computational fluid dynamics codes in three di­mensions for hydrosystems design work. When such codes are applied, the variance-reduction techniques described in Sec. 6.7 may be preferred to Monte Carlo simulation.

This chapter focuses on the basic principles and applications of Monte Carlo simulations in the reliability analysis of hydrosystems engineering problems. Section 6.2 describes some basic concepts of generating random numbers, fol­lowed by discussions on the classifications of algorithms for a generation of ran­dom variates in Sec. 6.3. Algorithms for generating univariate random numbers

TABLE 6.1 Basic Variable Statistics and Distributions for Evaluation of Monte Carlo Simulation of Convergence

Variable

Mean value

Coefficient of variation

Distribution function

X1

0.5

0.2

Uniform

X 2

1.5

0.4

Uniform

X 3

1.0

0.005

Lognormal

X 4

1.5

0.1

Lognormal

SOURCE: After Cheng et al. (1982).

are described in Sec. 6.4 for several commonly used distribution functions. In Sec. 6.5, attention is given to algorithms that generate multivariate ran­dom numbers. As reliability assessment involves mathematical integration, Sec. 6.6 describes several Monte Carlo simulation techniques for reliability evaluation. Given that Monte Carlo simulations, in essence, are sampling tech­niques, they provide only estimations, which inevitably are subject to certain degrees of errors. To improve the accuracy of the Monte Carlo estimation while reducing excessive computational time, several variance-reduction techniques are discussed in Sec. 6.7. Finally, resampling techniques are described in Sec. 6.8, which allow for assessment of the uncertainty of the quantity of interest based on the available random data without having to make assumptions about the underlying probabilistic structures.

Updated: 19 ноября, 2015 — 7:49 пп