词条 | Subset simulation |
释义 |
Subset simulation[1] is a method used in reliability engineering to compute small (i.e., rare event) failure probabilities encountered in engineering systems. The basic idea is to express a small failure probability as a product of larger conditional probabilities by introducing intermediate failure events. This conceptually converts the original rare event problem into a series of frequent event problems that are easier to solve. In the actual implementation, samples conditional on intermediate failure events are adaptively generated to gradually populate from the frequent to rare event region. These 'conditional samples' provide information for estimating the complementary cumulative distribution function (CCDF) of the quantity of interest (that governs failure), covering the high as well as the low probability regions. They can also be used for investigating the cause and consequence of failure events. The generation of conditional samples is not trivial but can be performed efficiently using Markov Chain Monte Carlo (MCMC). Subset Simulation takes the relationship between the (input) random variables and the (output) response quantity of interest as a 'black-box'. This can be attractive for complex systems where it is difficult to use other variance reduction or rare event sampling techniques that require prior information about the system behaviour. For problems where it is possible to incorporate prior information into the reliability algorithm, it is often more efficient to use other variance reduction techniques such as importance sampling. It has been shown that subset simulation is more efficient than traditional Monte Carlo simulation, but less efficient than Line Sampling, when applied to a fracture mechanics test problem [2]. Basic ideaLet X be a vector of random variables and Y = h(X) be a scalar (output) response quantity of interest for which the failure probability is to be determined. Each evaluation of h(.) is expensive and so it should be avoided if possible. Using direct Monte Carlo methods one can generate i.i.d. (independent and identically distributed) samples of X and then estimate P(F) simply as the fraction of samples with Y>b. However this is not efficient when P(F) is small because most samples will not fail (i.e., with Y≤b) and in many cases an estimate of 0 results. As a rule of thumb for small P(F) one requires 10 failed samples to estimate P(F) with a coefficient of variation of 30% (a moderate requirement). For example, 10000 i.i.d. samples, and hence evaluations of h(.), would be required for such an estimate if P(F)=0.001. Subset Simulation attempts to convert a rare event problem into more frequent ones. Let be an increasing sequence of intermediate threshold levels. From the basic property of conditional probability, The 'raw idea' of subset simulation is to estimate P(F) by estimating and the conditional probabilities for , anticipating efficiency gain when these probabilities are not small. To implement this idea there are two basic issues:
In the standard algorithm of subset simulation the first issue is resolved by using Markov Chain Monte Carlo.[3] The second issue is resolved by choosing the intermediate threshold levels {bi} adaptively using samples from the last simulation level. As a result, subset simulation in fact produces a set of estimates for b that corresponds to different fixed values of p = P(Y>b), rather than estimates of probabilities for fixed threshold values. There are a number of variations of subset simulation used in different contexts in applied probability and stochastic operations research[4] [5]For example, in some variations the simulation effort to estimate each conditional probability P(Y>bi|Y>bi-1) (i=2,...,m) may not be fixed prior to the simulation, but may be random, similar to the splitting method in rare-event probability estimation. [6]. More generic and flexible version of the simulation algorithms not based on Markov Chain Monte Carlo have been recently developed [7]. See also
Notes
References1. ^{{cite journal|last=Au|first=S.K.|author2=Beck, James L. |title=Estimation of small failure probabilities in high dimensions by subset simulation|journal=Probabilistic Engineering Mechanics|date=October 2001|volume=16|issue=4|pages=263–277|doi=10.1016/S0266-8920(01)00019-4|citeseerx=10.1.1.131.1941}} 2. ^{{cite book|last1=Zio|first1=E|title=Reliability, Risk, and Safety|last2=Pedroni|first2=N|chapter=Subset simulation and line sampling for advanced Monte Carlo reliability analysis|year=2009|doi=10.1201/9780203859759.ch94|isbn=978-0-415-55509-8}} 3. ^{{cite journal|last1=Au|first1=Siu-Kui|title=On MCMC algorithm for Subset Simulation|journal=Probabilistic Engineering Mechanics|date=2016|volume=43|pages=117–120|doi=10.1016/j.probengmech.2015.12.003}} 4. ^ {{Cite conference|first1=Manuel|last1=Villén-Altamirano|first2=José|last2=Villén-Altamirano | title=Restart: a straightforward method for fast simulation of rare events | book-title=Proceedings of the 26th Winter simulation conference|year=1994|series=WSC '94|pages=282–289|location=San Diego, CA, USA | publisher=Society for Computer Simulation International|id=acmid 194044|isbn=0-7803-2109-X|publication-place=Orlando, Florida, United States | url=http://dl.acm.org/citation.cfm?id=193201.194044}} 5. ^{{cite journal|last1=Botev|first1=Z. I. |last2=Kroese|first2=D. P.|title=An Efficient Algorithm for Rare-event Probability Estimation, Combinatorial Optimization, and Counting|journal=Methodology and Computing in Applied Probability|date=2008|volume=10|issue=4|pages=471–505|doi=10.1007/s11009-008-9073-7|citeseerx=10.1.1.399.7912 }} 6. ^{{cite journal|last1=Botev|first1=Z. I. |last2=Kroese|first2=D. P.|title=Efficient Monte Carlo simulation via the generalized splitting method|journal=Statistics and Computing|date=2012|volume=22|issue=1|pages=1–16|doi=10.1007/s11222-010-9201-4}} 7. ^{{cite journal|last1=Au|first1=Siu-Kui |last2=Patelli|first2=Edoardo|title=Rare event simulation in finite-infinite dimensional space|journal=Reliability Engineering & System Safety|date=2016|volume=148|pages=67–77|doi=10.1016/j.ress.2015.11.012}} 8. ^{{cite book|first=S.K. |last=Au |author2=Wang, Y. |title=Engineering Risk Assessment with Subset Simulation|year= 2014|publisher=John Wiley & Sons|location=Singapore|isbn=978-1-118-39804-3|url=}} 9. ^{{cite journal|last=Schuëller|first=G.I.|author2=Pradlwarter, H.J. |title=Benchmark study on reliability estimation in higher dimensions of structural systems – An overview|journal=Structural Safety|year=2007|volume=29|issue=3|pages=167–182|doi=10.1016/j.strusafe.2006.07.010}} 10. ^{{cite book|first=K.K. |last=Phoon|title=Reliability-Based Design in Geotechnical Engineering: Computations and Applications|year= 2008|publisher=Taylor & Francis|location=Singapore|isbn=978-0-415-39630-1|url=}} 11. ^{{cite journal|last=Zio|first=E.|author2=Pedroni, N. |title=How to effectively compute the reliability of a thermal–hydraulic nuclear passive system|journal=Nuclear Engineering and Design|year=2011|volume=241|issue=|pages=310–327|doi=10.1016/j.nucengdes.2010.10.029|citeseerx=10.1.1.636.2126}} 2 : Reliability analysis|Variance_reduction |
随便看 |
|
开放百科全书收录14589846条英语、德语、日语等多语种百科知识,基本涵盖了大多数领域的百科知识,是一部内容自由、开放的电子版国际百科全书。