# Study on the Confidence and Reliability of the Mean Seismic Probability Risk Model.

1. IntroductionThe mean seismic probability risk model has widely been used in seismic design and safety evaluation of critical infrastructures, such as nuclear power plants. Seismic probability risk assessment is one of seismic safety evaluation methodologies for nuclear power plants [1,2]. Nowadays, the U.S. Nuclear Regulatory Commission has adopted the mean seismic probability risk model as the basis for risk-informed decision-making. ASCE/SEI 43-05, whose basis was also the mean seismic probability risk model, was a risk-consistent seismic design criterion in the United States [3]. Kennedy [4] provided the technical basis for the performance-goal based approach presented in the American Society of Civil Engineering Standard ASCE/SEI 43-05. In China, the newest modified seismic design code of nuclear power plants [5] suggested using the seismic evaluation methodology of ASCE/SEI 43-05, whose basis was also the mean seismic probability risk model. In some engineering domains other than nuclear engineering, the mean seismic risk model has also been widely studied and implemented. Lu et al. [6] conducted the seismic risk assessment for a reinforced concrete frame designed according to Chinese codes based on the mean seismic probability risk model.

However, the confidence of the mean seismic probability risk model is unknown. In other words, the mean seismic probability risk model could not convey the sense of confidence directly. Ellingwood and Kinali [7] proposed that some decision-makers may not be comfortable with the mean as a point estimation and prefer a more conservative fractile of the risk distribution, especially if consequences are severe. Many studies on the interval model of seismic risk and structural reliability have been conducted. For the 2000 SAC Federal Emergency Management Agency (FEMA) steel moment frame guidelines [8], the format based on quantitative confidence statements regarding the likelihood of the performance objective was provided. Lu and Yu [9] studied the interval model of seismic risk based on the mean hazard curve and capacity and demand models. Katona [10] studied the uncertainty in seismic safety analysis, including seismic probability risk assessment and seismic margin assessment, based on p-box theory. De Leon and Ang [11] conducted confidence bounds analysis on structural reliability estimations for offshore platforms, and they suggested that the estimation of percentile values, instead of mean values, of the calculated risk should be specified to ensure sufficient low risk levels for decisions-making purposes. Ang [12] proposed that the confidence level of seismic reliability model considering epistemic uncertainties should be as high as 90% to 95% for reducing the influence of epistemic uncertainties on reliability analysis results in some critical systems, while the mean safety index (obtained with the total of the combined aleatory and epistemic uncertainty) has a confidence level of only around 50%. In nuclear engineering domain, high confidence of low probability of failure (HCLPF), which was defined as 95% confidence of less than 5% probability of failure and commonly used as the seismic capacity of structures, systems, and components (SSCs) of nuclear power plants, has been an important concept, which showed that SSCs of nuclear power plants should have low failure risk with high confidence facing earthquake disasters.

The mean seismic probability risk model has been accepted as a basis for risk-informed decision-making by U.S. Nuclear Regulatory Commission. However, for nuclear power plants as critical infrastructure, whose accident consequences are severe, the less failure probability risk with higher confidence should be required by decision-makers. In this paper, the confidence and error equations derivation of the approximate mean seismic probability risk are conducted. Meanwhile, the theoretical basis of ASCE/SEI 4305 code based on the approximate mean seismic probability risk model is extended to the approximate interval model. Confidence levels and error values based on these equations are then calculated. It is suggested that the seismic risk model considering high confidence levels instead of the mean probability risk model should be used for seismic design and safety evaluation of critical infrastructures such as nuclear power plants in the future.

2. Seismic Probability Risk Models and Confidence and Error Equations

2.1. The Approximate Point Estimation Model of Seismic Probability Risk. The limit state probability [P.sub.F] of seismic risk can be determined as the convolution of the seismic hazard curve and the fragility curve by either of the following two mathematically equivalent equations [13]:

[P.sub.F] = [[integral].sup.+[infinity].sub.0] H(a) d[F.sub.C](a)/da da (1)

[P.sub.F] = [[integral].sup.+[infinity].sub.0] [F.sub.C](a) dH(a)/da da (2)

where H(a) is the hazard curve (the annual frequency of exceeding amplitude a) and [F.sub.C] (a) is the fragility curve; (1) is equivalent to (2).

Modern seismic risk analysis, beginning with the seminal paper by Cornell [14] and some later studies by Ellingwood [7, 15], shows that the ground motion intensity can be represented by a type II distribution of extreme values, while the fragility function can be modelled by lognormal distribution.

The seismic hazard function approximation in closed form can be expressed as [7]

[mathematical expression not reproducible], (3)

where u is the scale parameter, [K.sub.H] is the shape parameter, and [k.sub.I] is the constant; that is, [mathematical expression not reproducible].

The seismic fragility can be expressed as [16-18]

[F.sub.C] (a)= [PHI] [ln(a/[a.sub.c])/[square root of [[beta].sup.2.sub.r] + [[beta].sup.2.sub.U]] = [PHI] [ln (a/[a.sub.C])/[[beta].sub.c]] (4)

in which [PHI][] represents the standard normal probability integral, [a.sub.c] is the median capacity, [[beta].sub.R] and [[beta].sub.u] represent the logarithmic standard deviation in capacity of aleatory uncertainty and epistemic uncertainty, respectively, pc is the square root of the sum of the squares (SRSS) of [[beta].sub.R] and [[beta].sub.U].

Substituting (3) and (4) into (1), the so-called "risk equation" can be obtained:

[P.sub.F,mean] = H([a.sub.c]) exp [[[K.sub.H][[beta].sub.c]/2].sup.2]] (5)

[mathematical expression not reproducible]. (6)

There are some approximations for (6): (1) the seismic hazard function, which is approximated by (3), is assumed to be linear on a log-log scale; (2) the seismic fragility function is assumed to follow the lognormal distribution.

2.2. The Approximate Interval Estimation Model of Seismic Probability Risk. The interval function of the fragility model is defined as [16-18]

[F.sub.C] (a) = [PHI] [ln(a/[a.sub.c]) + [[beta].sub.U][[PHI].sup.- 1](Q)/[[beta].sub.R]], (7)

where Q is the confidence parameter. Equation (7) can be transformed to

[F.sub.C] (a) = [PHI] [ln(a/[a.sub.c]) + [-[beta].sub.U][[PHI].sup.- 1](Q))))/[[beta].sub.R]], (7)

Equation (8) can be furtherly expressed as

[F.sub.c] (a) = [PHI] [ln(a/[a.sub.c,q])/[[beta].sub.R] (9)

[a.sub.c,Q] = [a.sub.c] exp (-[[beta].sub.u][[PHI].sup.-1] (Q)). (10)

Cornell [19] has pointed out that epistemic uncertainty can be modelled: (1) by assigning a lognormal (epistemic) distribution to H(x) and (2) by an epistemic lognormal distribution on median [a.sub.c]. McGuire et al. [20] suggested that, given the need to express seismic hazard with a single curve, the mean seismic hazard was the preferred single curve. In this paper, the mean seismic hazard curve H(x) is used, and the epistemic uncertainty is modelled only by the lognormal distribution on the median [a.sup.c]. Substituting (3) and (9) into (1), the interval estimation of the risk equation can be obtained:

[P.sub.F] (Q) = H([a.sub.c,Q]) exp [[[K.sub.H][[beta].sub.R]).sup.2]/2]. (11)

Substituting (10) into (11), the risk equation considering the confidence of the fragility model can be obtained:

[mathematical expression not reproducible]. (12)

Equation (12) can be furtherly expressed as

[mathematical expression not reproducible]. (13)

2.3. The Confidence and Error Equations of the Mean Seismic Risk Model. The mean seismic probability risk model has been widely used in seismic design and safety evaluation of some critical infrastructures, such as nuclear power plants. However, the mean seismic risk model has no direct information of confidence level, which the analyst has in the risk assessment. In order to analyze the confidence of the mean probability risk model, a new approach is proposed below in this paper.

Equation (6) represents the mean probability risk model, whose confidence level is unknown. Meanwhile, (13) is an interval risk model considering confidence level of the fragility model. The risk results could be, respectively, obtained from (6) and (13). Equation (6) represents the point estimation value of seismic risk, while (13) represents the interval estimation values of seismic risk with confidence Q. When Q is taken as a value [Q.sub.mean], [P.sub.F]([Q.sub.mean]) obtained from (13) would be equal to [P.sub.Fmean] obtained from (6). In this paper, [Q.sub.mean] is regarded as the confidence of [P.sub.Fmean] obtained from (6). Following this idea, this paper considers that when all parameters except Q of (6) and (13), which is taken as a value [Q.sub.mean], are the same, the value [Q.sub.mean] shall be regarded as the associated confidence level of the mean probability risk model. Ang [12] has discussed the confidence level of mean reliability model and also expressed the same idea with this paper. Ang [12] studied the design of an underground tunnel to resist a strong-motion earthquake and calculated the mean safety index and safety index with different confidence levels, and it was found that when the mean safety index was used in the design of the tunnel lining, the confidence for its safety was only about 50%. In other words, the mean safety index was equal to a safety index with about 50% confidence level of interval safety index with different confidence levels for the given case study.

When [P.sub.pmean] of (6) is equal to [P.sub.F](Q) of (13) and Q is taken as the value [Q.sub.mean], the equation can be expressed as

[mathematical expression not reproducible] (14)

It is assumed that all parameters of (6) and (13) except Q, which is taken as the value [Q.sub.mean], are the same. When neither [K.sub.H] nor [[beta].sub.U]] is equal to 0, the confidence of the mean probability risk model can be obtained as follows:

[Q.sub.mean] = [PHI] ([K.sub.H][[beta].sub.U]/2), (15)

where [K.sub.H], which represents the slope of the seismic hazard curve, can be furtherly expressed as [3]

[K.sub.H] = 1/log ([A.sub.R]), (16)

where [A.sub.R] represents the ratio of spectra acceleration of probability levels 0.1[H.sub.D] and [H.sub.D], in which [H.sub.D] represents the probability of exceedance at uniform hazard response spectrum (UHRS).

From the definitions of the variables [A.sub.R] and [[beta].sub.U], we could know that [A.sub.R] is larger than 1, so [K.sub.H] is always larger than 0, while [[beta].sub.u] is no less than 0. When [[beta].sub.u] is equal to 0, which means that the perfect knowledge is obtained for the model, (15) would be not suitable, and then the confidence equation of the mean probability risk model would be furtherly expressed as

[Q.sub.mean] = 100%. (17)

When there exists the epistemic uncertainty in the analysis or [[beta].sub.u] is larger than 0, the confidence level of the mean probability risk model would be obtained using (15). Equation (15) shows that [Q.sub.mean], which is the confidence of the mean probability risk model, is proportional to [K.sub.H]. From (15), it can be found that the confidence of the mean probability risk model in the regions with steep slopes of mean seismic hazard curves is higher than the regions with shallow slopes of the mean seismic hazard curves. So, for some critical infrastructures, such as nuclear power plants, which are located in the region with shallower slopes of mean seismic hazard curves, the mean probability risk model might be not appropriate. Typical values of [K.sub.H] in the West United States would be in the range of 3 to 6, while in the Central and Eastern United States, [K.sub.H] is typically 2.5 or less [13]. The confidence of the mean seismic risk model for the West United States and Central and Eastern United States is, respectively, shown in Figures 1 and 2. The results show that for the West United States the confidence of the mean seismic hazard model is in the range of nearly 55% to 96%, while for the Central and Eastern United States the confidence of the mean seismic hazard model is in the range of nearly 52% to 78%. So, when seismic risk is conducted in the Central and Eastern United States, the reliability of risk results based on the mean seismic risk model should be carefully examined and validated because of lower confidence levels.

The mean seismic risk model is obtained through calculating the mean of the distribution representing an epistemic uncertainty with the range of confidence levels. From (15) and (17), it is found that the confidence levels of the mean seismic risk model are always larger than 50%. The physical interpretation of this phenomenon could be expressed as follows: the mean seismic risk model is more sensitive to more severe alternative models considered in the analysis; in other words, the mean seismic risk model is more sensitive to the alternative models with higher confidence levels considered in the analysis.

The confidence levels of the mean probability risk model are calculated using the parameters [[beta].sub.u], which is usually in a range from 0.1 to 0.6, and [A.sub.R], which is usually in a range from 1.5 to 6. The results are listed in Table 1. It is found that, for different values of [[beta].sub.U] and [A.sub.R], the values of the confidence Q are different, so for different sites, the seismic risk probability model has different confidence levels, or it can also be said that the confidence levels of the seismic risk probability model are changed for different sites. The confidence levels of the seismic probability risk model are not high, the lowest and the highest values of which are 52.56% and 95.58%, respectively. We could also find that most confidence levels of the mean seismic risk model in the usual ranges of [[beta].sub.U] and [A.sub.R] are less than 90%, which is low for nuclear power plants.

In fact, the influence of epistemic uncertainty on risk results is not only related to confidence level Q. For fully considering the reliability of the mean seismic risk model, the error equation is provided and analyzed. The error equation, which is defined as the relative distance between [P.sub.p]([Q.sub.target]) and [P.sub.pmean], is expressed as follows:

[mathematical expression not reproducible], (18)

where [Q.sub.target] is the target confidence level; Err represents the relative error of [P.sub.p,mean] with respect to [P.sub.p]([Q.sub.target]).

The relative distance between Pp(Qo.gg) (the target confidence level is taken as 0.99 in this paper) and [P.sub.F,mean] is listed in Table 1 for usual range of and [K.sub.H] values. It is found that the error between Pp(Qo.gg) and PFmean is related not only to the confidence level, but also to the product of and KH. For some sites with large KH, even if the confidence level is high, the error is still large. So at sites with high [K.sub.H] the mean seismic risk model should be more carefully examined and validated because of their larger error.

3. Confidence and Error Analysis of ASCE/SEI 43-05 Seismic Risk-Consistent Design Methodology

3.1. Theoretical Basis of ASCE/SEI 43-05 Risk-Consistent Seismic Design Methodology. The theoretical basis of ASCE/SEI 43-05 risk-consistent seismic design is as follows [13].

For implementation purposes of ASCE/SEI 43-05, the parameters in the risk equation are reformulated in terms of design terms which are familiar to engineers such as structural engineers [13]. First, the p-percentile value of the mean fragility function is expressed as the design capacity:

[a.sub.P] = [a.sub.c] exp [-[X.sub.P][[beta].sub.c]], (19)

where [X.sub.P] is the (1 - p) percentile of the standard normal deviate.

Substituting the value [a.sub.c] of (19) into (6), the following is obtained:

[mathematical expression not reproducible]. (20)

It is assumed that aP is equal to the design response spectrum (DRS) times a seismic margin factor [F.sub.P]:

aP = FP x DRS = FP x [D.sub.F] x UHRS, (21)

where FP is the factor of safety representing the conservatism introduced by the seismic design criteria such as NUREG0800, and [D.sub.F] is the seismic design factor that is related to the uniform hazard response spectra (UHRS):

[mathematical expression not reproducible], (22)

in which [H.sub.D] is the annual probability of exceeding the UHRS.

Substituting (21) into (20) and making use of (22), we can obtain

[mathematical expression not reproducible]. (23)

Equation (23) can be transformed to

[mathematical expression not reproducible]. (24)

The design factor equations of ASCE/SEI 43-05 are expressed as follows:

DF = Maximum ([DF.sub.1], [DF.sub.2]) (25)

[DF.sub.2] = 06[([A.sub.R]).sup.[alpha]],

where [DF.sub.1] and [alpha] are defined in Table 2, respectively, for SDCs 3-5.

3.2. Interval Estimation Functions of ASCE/SEI 43-05 Risk-Consistent Seismic Design Methodology. The mean probability risk model is the theoretical basis of ASCE/SEI 43-05 risk-consistent seismic design methodology. For analyzing the confidence of the risk results based on ASCE/SEI 43-05 design methodology, the interval risk estimation functions of ASCE/SEI 43-05 are obtained. The confidence equations of ASCE/SEI 43-05 risk-consistent seismic design are as follows.

Equation (13) represents the risk equation considering the confidence of the fragility model. Substituting the value ac of (19) into (13), we can obtain

[mathematical expression not reproducible]. (26)

Substituting (21) into (26) and making use of (22), we can obtain

[mathematical expression not reproducible]. (27)

[mathematical expression not reproducible]. (28)

Equations (27) and (28) represent the functions of failure probability [P.sub.F] and design factor [D.sub.F] considering confidence, respectively.

Solving (27) for confidence Q of the mean failure probability [P.sub.Pmean], which could be calculated using (23), the following can be obtained:

[mathematical expression not reproducible]. (29)

Substituting (23) into (29), we could furtherly obtain

[mathematical expression not reproducible], (30)

where neither [K.sub.H] nor [[beta].sub.u] is equal to 0. According to the definitions of [K.sub.H] and [[beta].sub.u], [K.sub.H] is positive, while [[beta].sub.u] is no less than 0. When [[beta].sub.u] is equal to 0, (30) would not be suitable, so the confidence equation of [P.sub.F] could be furtherly expressed as follows:

[mathematical expression not reproducible]. (31)

Solving (28) for confidence Q of the mean design factor [D.sub.Fmean], which could be calculated using (24), the following can be obtained:

[mathematical expression not reproducible]. (32)

Substituting (24) into (32), we could furtherly obtain

QD[F.sub.mean] = [PHI]([K.sub.H][[beta].sub.U]/2), (33)

where neither [K.sub.H] nor [[beta].sub.u] is equal to 0. According to the definitions of [K.sub.H] and [[beta].sub.u], [K.sub.H] is positive, while [[beta].sub.u] is no less than 0. When [[beta].sub.u] is equal to 0, (33) would not be suitable, so the confidence equation of [D.sub.F] could be furtherly expressed as

[mathematical expression not reproducible]. (34)

When neither [K.sub.H] nor [[beta].sub.u] is equal to 0, the confidence functions of [P.sub.F] and [D.sub.F] of ASCE/SEI 43-05 could, respectively, be expressed as (30) and (33). It is found that the confidence levels of the two parameters are proportional to [K.sub.H] and [[beta].sub.u], which are the same as the confidence levels of the mean seismic risk model. This conclusion could be easily understood, as the mean seismic probability risk model is the theoretical basis of ASCE/SEI 43-05.

For three SDCs of ASCE/SEI 43-05 listed in Table 2 [3], the exact design factor (DF) and the confidence of [D.sub.F] could be, respectively, calculated using (24) and (33). When the exact design factor (DF) is calculated, [[beta].sub.c] is usually taken as the range of 0.3 to 0.6, as [[beta].sub.c] lies in the range of 0.3 to 0.6 for most structures, systems, and components (SSCs) of nuclear power plants (NPPs) [13]. The [D.sub.F] in ASCE/SEI 43-05 could be calculated using (25). The confidence levels of the exact DF are shown in Figure 3. For analyzing the confidence of ASCE DF, [[beta].sub.c] is still taken as the range of 0.3 to 0.6 and it is assumed that [[beta].sub.u] is equal to [[beta].sub.R] in this paper. So actually [[beta].sub.u] and [[beta].sub.R] are, respectively, taken as the range of 0.3 x [square root of 2]/2 to 0.6 x [square root of 2]/2 in Figure 3.

ASCE/SEI 43-05 [3] and Braverman et al. [13] have shown that ASCE [D.sub.F] approximation values approach the exact values. However, it is shown that the confidence levels of the exact [D.sub.F] are low, which are less than 90% from Figure 3.

For fully considering the reliability of achieved performance probabilities using ASCE/SEI 43-05 risk-consistent seismic design methodology, the error equation is obtained. The error equation, which is defined as the relative distance between point estimation function (see (23)) and interval estimation function for the target confidence level (see (27)) of achieved performance probabilities based on ASCE/SEI 43-05 risk-consistent seismic design methodology, is expressed as follows:

Err = exp [([[PHI].sup.-1] ([Q.sub.target]) - [K.sub.H][[beta].sub.U]/2 [K.sub.H][[beta].sub.u]] -1. (35)

It is found that the error equation of achieved performance probabilities based on ASCE/SEI 43-05 is the same as the mean seismic riskmodel, as the mean seismic probability risk model is the theoretical basis of ASCE/SEI 43-05.

3.3. Confidence and Error Analysis of a Design Example Based on the ASCE/SEI 43-05 Risk-Consistent Seismic Design Methodology. The typical normalized spectral acceleration hazard curves values are listed in Table 3 [3].

In ASCE/SEI 43-05, both of the following criteria need be obtained [3]:

(1) Less than 1% probability of acceptable performance for the DBE ground motion, which is defined as the Design Response Spectra (DRS)

(2) Less than 10% probability of acceptable performance for a ground motion equal to 150% of the DBE ground motion, which is defined as the DRS.

For the two criteria in ASCE/SEI 43-05, we could calculate the achieved performance probabilities [P.sub.F] using (23) and the confidence levels of the achieved performance probabilities PF using (30). For different SDCs, the values of HD are listed in Table 2. [A.sub.R] can be obtained from Table 3, while [K.sub.H] can be calculated using (16). [X.sub.p] means (1- p) percentile of the standard normal deviate (when p is equal to 1%, [X.sub.p] is 2.326, and when p = 10%, [X.sub.p] = L282). The seismic margin factors [F.sub.p] for different [beta] values are calculated by Kennedy [4], as listed in Table 4.

When the achieved performance probability [P.sub.f] is calculated, is [[beta].sub.c] usually taken as the range of 0.3 to 0.6, as [[beta].sub.c] lies in the range of 0.3 to 0.6 for most structures, systems, and components (SSCs) of nuclear power plants (NPPs) [13]. Kennedy [4] compared the results from the approximation risk equation with the precise convolution risk results, and it was found that the use of approximate power law hazard curve could give acceptable risk results. It could also be found that the failure probability designed according to ASCE/SEI 43-05 can meet two criteria [3, 4, 13]. However, it is shown that the confidence levels of failure probability [P.sub.F] are changed for different sites from Figures 4 and 5. The confidence levels of two example sites are low, which are less than 80%. For analyzing the confidence of failure probability PF, pc is still taken as the range of 0.3 to 0.6 and it is assumed that [[beta].sub.U] is equal to [[beta].sub.R] in this paper. So actually [[beta].sub.u] and [[beta].sub.R] are, respectively, taken as the range of 0.3 x [square root of 2]/2 to 0.6 x [square root of 2]/2 in Figures 4 and 5. We could also find that the confidence levels of the achieved performance probabilities in the region of California with steep slopes of the mean seismic hazard curves are larger than the confidence levels of the achieved performance probabilities in the regions of Eastern United States with shallow slopes of the mean seismic hazard curves. In a word, the uniform risk could be obtained using ASCE/SEI 43-05 criteria, but uniform high confidence level of the specified target risk could not be achieved based on ASCE/SEI 43-05 criteria.

The error analysis for achieved performance probabilities for the above example is conducted. The results are shown in Table 5. It could be found that error results between risk probabilities with target confidence level (the target confidence level is taken as 0.99 in this paper) and mean achieved performance probabilities based on ASCE/SEI 4305 are changed for different sites and the error values are large for some sites.

4. Conclusions

In this study, the confidence equation of the mean seismic risk model is derived. It could be found that the confidence levels of the mean probability risk model for different sites are changed, and the confidence levels are also low for most sites.

The error equation between [P.sub.f]([Q.sub.target]) and [P.sub.f,mean] is given and analyzed. It could be found that the error between [P.sub.F]([Q.sub.target]) (the target confidence is taken as 0.99 in the paper) and [P.sub.F,mean] is related not only to the confidence level, but also to the product of [[beta].sub.u] and [K.sub.H]. So for some sites with large [K.sub.H], even if the confidence level is high, the error is still large. So at sites with high [K.sub.H] the mean seismic risk model should be more carefully examined and validated because of their larger error.

The confidence equation derivation of design factor DF and failure probability [P.sub.F] of ASCE/SEI 43-05 code are conducted. It has been found that the confidence levels of DF and Pp of ASCE/SEI 43-05 are changed, while the confidence levels are low. It is also found that the uniform risk could be obtained, but the uniform high confidence level of the specified target risk could not be achieved using ASCE/SEI 43-05 criteria.

The error equation of achieved performance probabilities based on ASCE/SEI 43-05 risk-consistent seismic design methodology is obtained. It could be found that error results between target confidence level (the target confidence is taken as 0.99 in the paper) and achieved performance probabilities based on ASCE/SEI 43-05 are changed for different sites and the error values are large for some sites.

It is suggested that seismic risk model considering high confidence levels instead of the mean probability risk model should be used for seismic design and safety evaluation of critical infrastructures such as nuclear power plants in the future.

https://doi.org/ 10.1155/2017/2480940

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The financial support received from the National Science Foundation of China (Grant nos. 51678209, 51378162, and 91315301), the Research Fund of Ministry of Science and Technology of China (2013BAJ08B01), and the Open Research Fund of State Key Laboratory for Disaster Reduction in Civil Engineering (SLDRCE12-MB-04) is gratefully appreciated.

References

[1] EPRI, "A Methodology for Assessment of Nuclear Power Plant Seismic Margin," EPRI NP-6041SL, Revision 1, EPRI, Palo Alto, Calif, USA, 1991.

[2] R. Kassawara, "Seismic Probabilistic Risk Assessment Implementation Guide," Tech. Rep., Electric Power Research Institute, 2003.

[3] ASCE, "Seismic design criteria for structures, systems, and components in nuclear facilities," ASCE/SEI 43-05, 2005.

[4] R. P. Kennedy, "Performance-goal based (risk informed) approach for establishing the SSE site specific response spectrum for future nuclear power plants," Nuclear Engineering and Design, vol. 241, no. 3, pp. 648-656, 2011.

[5] Institute of Engineering Mechanics, China Earthquake Administration. "The submitted version of the code for seismic design of nuclear power plants GB50267-97." Harbin: Institute of Engineering Mechanics China, 2012 (Chinese).

[6] D. Lu, X. Yu, M. Jia, and G. Wang, "Seismic risk assessment for a reinforced concrete frame designed according to Chinese codes," Structure and Infrastructure Engineering, vol. 10, no. 10, pp. 1295-1310, 2014.

[7] B. R. Ellingwood and K. Kinali, "Quantifying and communicating uncertainty in seismic risk assessment," Structural Safety, vol. 31, no. 2, pp. 179-187, 2009.

[8] C. A. Cornell, F. Jalayer, R. O. Hamburger, and D. A. Foutch, "Probabilistic basis for 2000 SAC federal emergency management agency steel moment frame guidelines," Journal of Structural Engineering, vol. 128, no. 4, pp. 526-533, 2002.

[9] D. Lu and X. Yu, "Theoretical study of probabilistic seismic risk assessment based on analytical functions of seismic fragility," Journal of Building Structures, vol. 34, no. 10, pp. 41-48, 2013 (Chinese).

[10] T. J. Katona, "Options for the treatment of uncertainty in seismic safety assessment of nuclear power plants," Journal of Disaster Research, vol. 8, no. 3, pp. 465-472, 2013.

[11] D. De Leon and A. H.-S. Ang, "Confidence bounds on structural reliability estimations for offshore platforms," Journal of Marine Science and Technology, vol. 13, no. 3, pp. 308-315, 2008.

[12] A. H.-S. Ang, "Target reliability for design of complex systems Role of PDEM," in Proceedings of the Symposium on Reliability of Engineering System, SRES 2015, Hangzhou, China, October 2015.

[13] J. I. Braverman, J. Xu, B. R. Ellingwood, C. J. Costantino, R. J. Morante, and C. H. Hofmayer, "Evaluation of the Seismic Design Criteria in ASCE/SEI Standard 43-05 for Application to Nuclear Power Plants," NUREG/CR-6926, BNL-NUREG77569-2007, U.S. Nuclear Regulatory Commission, 2007.

[14] C. A. Cornell, "Engineering seismic risk analysis," Bulletin of the Seismological Society of America, vol. 58, pp. 1583-1606, 1968.

[15] B. Ellingwood, "Validation studies of seismic PRAs," Nuclear

Engineering and Design, vol. 123, no. 2-3, pp. 189-196, 1990.

[16] J. W. Reed and P. Kennedy R, "Methodology for Developing Seismic Fragilities," EPRI TR-103959, EPRI, Palo Alto, CA, USA, 1994.

[17] R. Campbell, G. Hardy, and K. Merz, "Seismic Fragility Application Guide," EPRI Report 1002988, EPRI, Palo Alto, CA, USA, 2002.

[18] R. Kennedy, G. Hardy, and K. Merz, "Seismic Fragility Applications Guide Update," Tech. Rep., Electric Power Research Institute, Palo Alto, CA, USA, Product ID Number 1019200, 2009.

[19] C. A. Cornell, "Risk-based structural design," in Proceedings of the Symposium on Risk Analysis, A. S. Nowak, Ed., pp. 37-48, University of Michigan, Michigan, MI, USA, 1994.

[20] R. K. McGuire, C. A. Cornell, and G. R. Toro, "The case for using mean seismic hazard," Earthquake Spectra, vol. 21, no. 3, pp. 879-886, 2005.

Xiao-Lei Wang (1,2) and Da-Gang Lu (1,2,3)

(1) School of Civil Engineering, Harbin Institute of Technology, Harbin 150090, China

(2) Key Lab of Structures Dynamic Behavior and Control of China Ministry of Education, Harbin Institute of Technology, Harbin 150090, China

(3) State Key Laboratory for Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China

Correspondence should be addressed to Da-Gang Lu; ludagang@hit.edu.cn

Received 26 November 2016; Revised 26 January 2017; Accepted 30 July 2017; Published 10 October 2017

Academic Editor: Eugenijus Uspuras

Caption: FIGURE 1: Confidence levels of mean seismic probability risk model used in West United States sites.

Caption: FIGURE 2: Confidence levels of mean seismic risk probability model used in Central and Eastern United States sites.

Caption: FIGURE 3: Confidence levels of exact [D.sub.F] for SDCs 3, 4, and 5 for the range of [A.sub.R] and [[beta].sub.c] values.

Caption: FIGURE 4: Confidence levels of achieved performance probabilities for SDCs 3 and 4 based on the first and second criteria of ASCE/SEI 43-05 for representative hazard curves.

Caption: FIGURE 5: Confidence levels of achieved performance probabilities for SDC 5 based on both the first and second criteria of ASCE/SEI 43-05 for representative hazard curves.

TABLE 1: Confidence levels and error values of mean seismic probability risk model for the usual range of [[beta].sub.U] and [K.sub.H] values. Confidence levels and error results between risk model with 99% confidence and mean risk model (in brackets) [A.sub.R] [K.sub.H] [[beta].sub.U] = 0.1 [[beta].sub.U] = 0.2 1.5 5.6789 61.18% (2.1895) 71.49% (6.3685) 1.75 4.1146 58.15% (1.3930) 65.96% (3.8346) 2 3.3219 56.60% (1.0495) 63.01% (2.7618) 2.25 2.8394 55.64% (0.8594) 61.18% (2.1895) 2.5 2.5129 55.00% (0.7385) 59.92% (1.8374) 2.75 2.2762 54.53% (0.6547) 59.00% (1.5998) 3 2.0959 54.17% (0.5930) 58.30% (1.4286) 3.25 1.9536 53.89% (0.5456) 57.74% (1.2993) 3.5 1.8380 53.66% (0.5079) 57.29% (1.1981) 3.75 1.7421 53.47% (0.4771) 56.91% (1.1166) 4 1.6610 53.31% (0.4515) 56.60% (1.0495) 4.25 1.5914 53.17% (0.4298) 56.32% (0.9933) 4.5 1.5309 53.05% (0.4112) 56.08% (0.9453) 4.75 1.4778 52.95% (0.3950) 55.87% (0.9039) 5 1.4307 52.85% (0.3807) 55.69% (0.8677) 5.25 1.3886 52.77% (0.3681) 55.52% (0.8358) 5.5 1.3507 52.69% (0.3568) 55.37% (0.8075) 5.75 1.3164 52.62% (0.3466) 55.24% (0.7821) 6 1.2851 52.56% (0.3374) 55.11% (0.7592) Confidence levels and error results between risk model with 99% confidence and mean risk model (in brackets) [A.sub.R] [[beta].sub.U] = 0.3 [[beta].sub.U] = 0.4 1.5 80.28% (11.3305) 87.20% (13.9460) 1.75 73.14% (7.2461) 79.47% (10.8745) 2 69.09% (5.1830) 74.68% (8.1008) 2.25 66.49% (4.0473) 71.49% (6.3685) 2.5 64.69% (3.3476) 69.24% (5.2539) 2.75 63.36% (2.8784) 67.55% (4.4937) 3 62.34% (2.5433) 66.25% (3.9476) 3.25 61.53% (2.2926) 65.20% (3.5383) 3.5 60.86% (2.0979) 64.34% (3.2210) 3.75 60.31% (1.9424) 63.62% (2.9681) 4 59.84% (1.8152) 63.01% (2.7618) 4.25 59.43% (1.7092) 62.49% (2.5903) 4.5 59.08% (1.6194) 62.03% (2.4455) 4.75 58.77% (1.5423) 61.62% (2.3215) 5 58.50% (1.4753) 61.26% (2.2141) 5.25 58.25% (1.4165) 60.94% (2.1201) 5.5 58.03% (1.3645) 60.65% (2.0372) 5.75 57.83% (1.3180) 60.38% (1.9633) 6 57.64% (1.2763) 60.14% (1.8971) Confidence levels and error results between risk model with 99% confidence and mean risk model (in brackets) [A.sub.R] [[beta].sub.U] = 0.5 [[beta].sub.U] = 0.6 1.5 92.22% (12.1224) 95.58% (7.3452) 1.75 84.82% (13.4364) 89.15% (13.8175) 2 79.69% (10.9961) 84.05% (13.1603) 2.25 76.11% (8.9240) 80.28% (11.3305) 2.5 73.51% (7.4455) 77.45% (9.7071) 2.75 71.53% (6.3888) 75.27% (8.4360) 3 69.99% (5.6114) 73.53% (7.4551) 3.25 68.74% (5.0212) 72.11% (6.6895) 3.5 67.71% (4.5602) 70.93% (6.0810) 3.75 66.84% (4.1913) 69.94% (5.5885) 4 66.10% (3.8897) 69.09% (5.1830) 4.25 65.46% (3.6389) 68.35% (4.8439) 4.5 64.90% (3.4271) 67.70% (4.5566) 4.75 64.41% (3.2458) 67.12% (4.3101) 5 63.97% (3.0889) 66.61% (4.0963) 5.25 63.58% (2.9517) 66.15% (3.9093) 5.5 63.22% (2.8307) 65.73% (3.7442) 5.75 62.90% (2.7231) 65.35% (3.5974) 6 62.60% (2.6268) 65.01% (3.4660) TABLE 2: Summary of earthquake design provisions of ASCE-SEI 43-05 [3]. Seismic Design Target performance Probability Category (SDC) goal ([P.sub.F]) ratio ([R.sub.p]) 3 1 x [10.sup.-4] 4 4 4 x [10.sup.-5] 10 5 1 x [10.sup.-5] 10 Seismic Design Hazard exceedance [DF.sub.1] [alpha] Category (SDC) probability ([H.sub.D]) 3 4 x [10.sup.-4] 0.8 0.40 4 4 x [10.sup.-4] 1.0 0.80 5 1 x [10.sup.-4] 1.0 0.80 TABLE 3: Typical normalized spectral acceleration hazard curve values of ASCE/SEI 43-05 [3]. Eastern US California [H.sub.(SA)] 1Hz 10 Hz 1Hz 10 Hz SA SA SA SA 5 x [10.sup.-2] 0.014 0.018 0.087 0.046 2 x [10.sup.-2] 0.027 0.034 0.13 0.072 1 x [10.sup.-2] 0.045 0.055 0.175 0.100 5 x [10.sup.-3] 0.07 0.089 0.236 0.139 2 x [10.sup.-3] 0.143 0.169 0.351 0.215 1 x [10.sup.-3] 0.235 0.275 0.474 0.334 5 x [10.sup.-4] 0.383 0.424 0.629 0.511 2 x [10.sup.-4] 0.681 0.709 0.814 0.762 1 x [10.sup.-4] 1.00 1.00 1.00 1.00 5 x [10.sup.-5] 1.46 1.41 1.23 1.22 2 x [10.sup.-5] 2.35 2.13 1.61 1.51 1 x [10.sup.-5] 3.27 2.88 1.89 1.76 5 x [10.sup.-6] 4.38 3.65 2.2 2.05 2 x [10.sup.-6] 6.44 4.62 2.68 2.42 1 x [10.sup.-6] 8.59 5.43 3.1 2.72 5 x [10.sup.-7] 10.34 6.38 3.58 3.06 2 x [10.sup.-7] 13.21 7.9 4.24 3.56 1 x [10.sup.-7] 15.9 9.28 4.67 3.84 TABLE 4: Seismic margin factors for different [beta] values calculated by Kennedy [4]. [beta] [F.sub.1%] [F.sub.5%] [F.sub.10%] 0.30 1.10 1.35 1.5 0.40 1 1.31 1.52 0.50 1 1.41 1.69 0.60 1 1.5 1.87 [beta] [F.sub.50%] [F.sub.70%] 0.30 2.2 2.58 0.40 2.54 3.13 0.50 3.2 4.16 0.60 4.04 5.53 TABLE 5: Error results for both the first and second criteria of ASCE-SEI 43-05 for representative hazard curves. Error results between risk model with 99% confidence and mean risk model SDC Hazard curves [[beta].sub.c] = 0.3 [[beta].sub.c] = 0.4 [F.sub.1%] = 1.1 [F.sub.1%] = 1.0 3 EUS 1 Hz 1.2163 1.8022 EUS 10 Hz 1.4066 2.1043 Calif 1Hz 3.1584 4.9483 Calif 10 Hz 2.3304 3.6001 4 EUS 1 Hz 1.2163 1.8022 EUS 10 Hz 1.4066 2.1043 Calif 1Hz 3.1584 4.9483 Calif 10 Hz 2.3304 3.6001 5 EUS 1 Hz 1.3967 2.0885 EUS 10 Hz 1.6317 2.4652 Calif 1Hz 3.4400 5.4025 Calif 10 Hz 4.1385 6.5117 Error results between risk model with 99% confidence and mean risk model SDC Hazard curves [[beta].sub.c] = 0.5 [[beta].sub.c] = 0.6 [F.sub.1%] = 1.0 [F.sub.1%] = 1.0 3 EUS 1 Hz 2.4889 3.2777 EUS 10 Hz 2.9285 3.8777 Calif 1Hz 7.0248 9.2106 Calif 10 Hz 5.1107 6.8066 4 EUS 1 Hz 2.4889 3.2777 EUS 10 Hz 2.9285 3.8777 Calif 1Hz 7.0248 9.2106 Calif 10 Hz 5.1107 6.8066 5 EUS 1 Hz 2.9056 3.8464 EUS 10 Hz 3.4559 4.5957 Calif 1Hz 7.6479 9.9412 Calif 10 Hz 9.1067 11.5157

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Wang, Xiao-Lei; Lu, Da-Gang |

Publication: | Science and Technology of Nuclear Installations |

Article Type: | Report |

Date: | Jan 1, 2017 |

Words: | 6928 |

Previous Article: | New Method to Measure Crack Extension in Nuclear Graphite Based on Digital Image Correlation. |

Next Article: | ARES: A Parallel Discrete Ordinates Transport Code for Radiation Shielding Applications and Reactor Physics Analysis. |

Topics: |