[论文解读] Statistical considerations on safety analysis

Л. Пал, Mihály Makai|ArXiv.org|Nov 16, 2005

Risk and Safety Analysis参考文献 1被引用 5

一句话总结

本文挑战核安全分析中普遍采用的0.95/0.95方法，揭示其统计缺陷：重复运行可能以不可忽略的概率产生更高的容忍限值，从而增加未被检测到的限值违规风险。本文提出符号检验作为更可靠的替代方法，用于估计安全关键系统中的分位数，尤其适用于非对称分布情况，并纠正了多变量情况下Wilks公式的应用。

ABSTRACT

Alerting experience with a well-acknowledged safety analysis code initiated the authors to pay attention to safety issues of complex systems. Their first concern was the statistical characteristics of such a code. We point out a remarkable weakness of the so called 0.95/0.95 methodology: when repeating the search for the tolerance limit, we get a higher value with non-negligible probability. We propose the sign test as an alternative method. We point out the correct form of Wilks' formula when the number of parameters subjected to limitation is two or more.

研究动机与目标

• 识别并解决核反应堆安全分析中广泛使用的0.95/0.95方法的统计缺陷。
• 证明0.95/0.95方法的重复应用可能以不可忽略的概率产生更高的容忍限值，从而削弱安全保证。
• 提出符号检验作为安全关键系统中分位数估计的统计上可靠替代方法。
• 纠正多个输出变量受限制时Wilks公式的应用。
• 倡导在复杂系统的最优估计建模中整合概率特性，包括混沌分析与输出分布分析。

提出的方法

• 使用顺序统计量基于N次独立运行的最大值估计上容忍限值，应用公式 β = 1 − γ^N 计算覆盖概率。
• 采用贝叶斯方法估计真实接受概率w超过阈值ω的概率，基于N次运行中故障的二项分布。
• 提出符号检验作为0.95/0.95方法的非参数替代方案，尤其适用于非对称或未知分布情况。
• 推导k个独立最大值中的k个超过参考最大值的概率，证明 P_k = 1 − k/(n+1)，该结果支撑符号检验的统计有效性。
• 纠正多变量情况（n ≥ 2）下Wilks公式的应用，确保覆盖概率保持分布无关。
• 使用对数正态分布输出的蒙特卡洛模拟，实证检验0.95/0.95方法在非对称条件下的失效模式。

实验结果

研究问题

• RQ1 为何0.95/0.95方法在非对称分布中会失效？重复运行产生高于原始值的容忍限值的概率是多少？
• RQ2 如何在安全分析中应用符号检验，以提供比0.95/0.95方法更稳健的分位数估计？
• RQ3 当多个输出变量受安全约束时，Wilks公式的正确形式是什么？
• RQ4 最大顺序统计量在重复模拟中的统计行为如何？其对安全裕量估计有何影响？
• RQ5 模型输出的概率特性如何影响复杂系统中最佳估计安全分析的可靠性？

主要发现

• 在非对称分布中，0.95/0.95方法可能失效，因为重复运行可能以不可忽略的概率产生更高的容忍限值，从而增加限值违规风险。
• 符号检验提供了一种分布无关的分位数估计方法，可保证覆盖概率，且不依赖于底层分布。
• 对于N次运行，最大值超过γ分位数的概率为 β = 1 − γ^N，这是一个已确立的结果，用于验证置信水平。
• 当涉及多个变量（n ≥ 2）时，正确应用Wilks公式可确保联合覆盖概率保持分布无关。
• k个独立最大值中有k个超过参考最大值的概率为 P_k = 1 − k/(n+1)，该结果支持符号检验的统计有效性。
• 即使在N次运行中无故障发生，真实接受概率w低于阈值ω的概率仍为非零，凸显了采用贝叶斯或符号检验方法估计的必要性。