IEC 61710：幂律模型——拟合优度检验与置信区间

要点提示：IEC 61710 为可修系统的可靠性建模提供了完整的统计框架。该标准基于幂律模型（也称威布尔过程）为可靠性工程师提供了一套系统的方法来评估故障数据、检验模型适用性并计算置信区间。

一、标准范围与数学基础

IEC 61710 规定了将幂律模型拟合到可修系统故障时间数据的程序，以及评估拟合充分性的方法。幂律模型，也称为威布尔过程或具有幂律强度函数的非齐次泊松过程（NHPP），将瞬时故障强度描述为时间的函数：

幂律强度函数：

λ(t) = (β/θ) · (t/θ)^β-1

其中：
λ(t) = 时间 t 处的瞬时故障强度
β = 形状参数（无量纲）
θ = 尺度参数（时间单位）

截至时间 t 的累积故障数为：
E[N(t)] = (t/θ)^β

形状参数 β 决定了系统的可靠性趋势。当 β = 1 时，过程退化为具有恒定故障强度的齐次泊松过程（HPP），表明既无可靠性增长也无退化。当 β < 1 时，故障强度随时间递减，表明可靠性增长——这是在研制试验阶段观察到的典型模式。当 β > 1 时，故障强度随时间递增，表明系统退化或耗损失效。这种简洁而有力的特性使幂律模型成为可修系统可靠性工程中最广泛使用的工具之一。

注意事项：可靠性分析中一个常见错误是将幂律模型应用于不可修系统，或应用于”修复”未能将系统恢复到故障前可比运行状态的数据。该模型假设”最小维修”——恢复功能但不改变整体故障强度过程。

标准涵盖三个主要应用：模型参数的点估计、形状参数 β 和尺度参数 θ 的区间估计（置信区间），以及判定幂律模型是否充分拟合观测故障数据的拟合优度检验。同时描述了极大似然估计（MLE）和图形估计方法，使工程师能够根据数据质量和可用计算工具灵活选择。

二、统计程序与拟合优度检验

IEC 61710 定义了多种拟合优度检验方法。检验的选择取决于故障数据是时间截尾（在预定时间 T 停止试验）还是故障截尾（在预定故障数 n 后停止试验）。

检验方法	检验统计量	截尾类型	零分布
Cox–Lewis（修正）	在 HPP 假设下，转换后的间隔时间服从指数分布	时间截尾	自由度为 2n 的 χ²
军用手册（MIL-HDBK-189）	CM = 2 ∑ ln(T/t_i)	时间截尾	自由度为 2n 的 χ²
Bartlett 检验（修正）	B = 2n [ln(∑t_i/n) − ∑ln(t_i)/n] / [1 + (n+1)/(6n)]	故障截尾	自由度为 n−1 的 χ²
Laplace 趋势检验	U = (∑t_i − nT/2) / (T √(n/12))	时间截尾	标准正态 N(0,1)

实践中应用最广泛的是 MIL-HDBK-189 检验（常称为 Crow 检验）。检验统计量 CM 与卡方分布的临界值进行比较。若检验统计量落在接受域内，则无充分证据拒绝幂律模型。若结果显著，表明幂律模型可能无法充分描述故障过程，应考虑对数线性过程或调制幂律模型等替代方案。

标准还提供了计算置信区间的详细程序。对于形状参数 β，置信区间基于如下事实：2nβ/βˆ（其中 βˆ 是 β 的 MLE）对于故障截尾数据服从自由度为 2(n−1) 的卡方分布。对于时间截尾数据，自由度调整有所不同，反映了停止条件带来的额外不确定性。这些置信区间对于判断观测到的可靠性趋势是否具有统计显著性（而非仅仅随机波动）至关重要。

实践建议：始终计算 β 的点估计和 90% 或 95% 置信区间。若置信区间包含 1.0，则不能判定系统表现出可靠性增长或退化——数据与恒定故障强度（HPP）一致。在此情况下，应考虑延长试验时间或收集更多故障数据。

三、工程设计要点与实际应用

IEC 61710 中的幂律模型在多个行业得到广泛应用。在航空航天和国防领域，它是 MIL-HDBK-189 和 Crow/AMSAA 方法中规定的研制试验可靠性增长跟踪标准模型。在汽车行业，它用于分析车辆修复后返回使用的保修索赔数据。在工业设备维护中，该模型通过早期检测耗损失效（β > 1）来帮助优化预防性维护间隔。

最强大的工程应用之一是使用幂律模型进行可靠性增长规划。通过在对数-对数坐标上绘制累积故障数与累积试验时间，工程师可以估计当前可靠性水平并预测在持续试验下的未来可靠性。对数-对数图的斜率为 β，截距提供初始可靠性水平的信息。典型的可靠性增长计划可能设定 β 目标值为 0.5 至 0.7，表示随试验推进故障强度显著降低。

β 范围	可靠性趋势	工程解释
β < 0.5	强增长	纠正措施积极；早期可能过度试验
0.5 ≤ β < 0.8	中度增长	FRACAS 有效；健康的发展计划
0.8 ≤ β < 1.0	弱增长	系统成熟；试验收益递减
β = 1.0	无趋势	恒定故障强度；适用 HPP
1.0 < β < 1.5	轻度退化	早期耗损；审查预防性维护计划
β ≥ 1.5	严重退化	故障率加速上升；需要重新设计

重要提醒：幂律模型假设故障过程是非齐次泊松过程，要求给定强度函数后故障间隔时间条件独立。若因共因机制（如批次缺陷、环境应力）导致故障聚集，模型可能产生误导性结果。务必结合工程判断和故障模式分析来进行统计检验。

IEC 61710 的实际实施需要重视数据质量。故障时间必须以连续时间（累积试验时间）记录，而非离散的运行间隔。右删失数据（在观察期结束时未发生故障的系统）可以处理，但标准提供了处理中止数据的具体程序。对于多个相同受试单元的系统，标准建议合并数据，假设共同的 β 但可能不同的 θ 值，反映不同的起始可靠性但相同的改进速率。

现代软件工具使 IEC 61710 的计算方面变得简单，但工程解释仍然具有挑战性。一个 β = 0.6 且置信区间狭窄 [0.55, 0.65] 的系统提供了可靠性增长的强有力证据，允许自信地预测未来可靠性。相比之下，来自有限数据的 β = 0.6 且置信区间宽泛 [0.35, 1.05] 应谨慎解释——系统可能根本没有改进。标准强调置信区间宽度与观测到的故障数成反比，为规划试验时间提供了定量依据。

常见问题解答

问题1：应用 IEC 61710 所需的最小故障数是多少？

标准建议至少 5 到 10 个故障才能进行有意义的参数估计。少于 5 个故障时，置信区间变得极宽，拟合优度检验也缺乏检测模型不适用的统计功效。对于研制试验中的可靠性增长跟踪，实际最低要求是 20 个故障以建立可靠趋势。

问题2：幂律模型能否处理具有不同试验持续时间的多个系统？

可以。IEC 61710 提供了组合多个系统数据的程序。标准建议在所有系统中拟合共同的 β，同时允许个体 θ 参数。当所有系统经历类似的纠正措施过程但起始可靠性水平不同时，这种方法是适当的。累积试验时间是各系统试验时间的总和。

问题3：幂律模型与标准威布尔分布有何不同？

标准（两参数）威布尔分布对不可修群体的首次故障时间建模。IEC 61710 中的幂律模型（威布尔过程）对可修系统随时间变化的故障强度建模。两者数学形式相同，但描述的是根本不同的物理情境。对可修系统数据使用威布尔分布是常见但严重的统计错误。

问题4：当幂律模型拟合不充分时，应考虑哪些替代模型？

当拟合优度检验拒绝幂律模型时，可考虑对数线性模型（指数强度函数）、调制幂律模型（也称 Musa-Okumoto 过程）或分段幂律模型（不同试验阶段使用不同 β 值）。选择取决于残差分析和趋势图揭示的偏离幂律假设的性质。