[发明专利]基于伽马-广义逆威布尔分布的非精确失效模型构建方法

说明书

本发明公开了一种基于伽马‑广义逆威布尔分布的非精确失效模型构建方法，其包括获取元件寿命的样本数据，建立各分布函数族和期望的非精确模型，构建各分布函数族和期望的区间值上下界的计算模型，计算各分布函数族和期望的区间值上下界。本发明根据在电气应力作用下元件寿命X服从参数为a,b,C的广义逆威布尔分布，转变参数C的先验分布服从参数为α,β的伽马分布，利用贝叶斯公式对元件寿命的各分布函数族和期望进行非精确推断，用非精确结果代替精确结果，有效地避免了由精确概率带来的误差而导致预测规划的偏差，为电力系统非精确可靠性评估和预测提供了依据。

技术领域

本发明属于电力系统和可靠性工程领域，具体涉及一种在电气应力作用下电力系统元件非精确失效模型构建方法。

背景技术

目前广泛采用的失效分布有指数分布、正态分布、对数正态分布以及威布尔分布等。除上述分布外，逆威布尔分布(IW)在可靠性领域也有着广泛地应用，其可用于模拟各种失效特征，如使用寿命和损耗周期；也可用于确定以可靠性为中心的维修活动的成本效益和维护周期。三参数广义逆威布尔分布(GIW)是逆威布尔分布的推广形式，广泛地用于生存数据建模、可靠性问题和疲劳寿命的研究。在许多可靠性工程问题中，如果风险函数为非单调函数时，威布尔分布作为寿命分布将不再合适；当风险函数为单峰函数时，那么广义逆威布尔分布可能是一个合适的模型。相对于指数分布和威布尔分布而言，广义逆威布尔分布在描述具有退化趋势的机械部件故障问题时具有更好的拟合效果。同时，在机械或电气元件寿命试验中，广义逆威布尔分布也适用于建模。总之，广义逆威布尔分布模型包含了电气绝缘应力强度的各种特性，可以让IDHR(先上升后下降的故障率)的特性和ALT(寿命加速试验法)所得的数据更为合理。

关于广义逆威布尔分布的参数估计，目前有相当多基于精确概率的方法，如最大似然估计等。但这些估计方法大多建立在大量样本数据的基础上，而电力系统的某些元件是长寿命的昂贵器件，进行一次寿命试验要花费大量的人力、物力，从而很难获得足够多的样本数据。由于样本数据的缺乏导致参数估计的非精确性，如果再利用传统的精确概率理论去描述研究对象的概率信息可能会产生较大的误差。

发明内容

本发明的发明目的是：为了解决现有技术中在样本数据缺乏情况下无法精确估计分布参数等问题，本发明提出了一种在电气应力作用下基于伽马-广义逆威布尔分布的非精确失效模型构建方法。

本发明的技术方案是：一种基于伽马-广义逆威布尔分布的非精确失效模型构建方法，包括以下步骤：

A、在电气应力作用下电力系统元件寿命X服从参数为a,b,C的广义逆威布尔分布GIW(a,b,C)，其中a为尺度参数，b为形状参数，转变参数C的先验分布服从参数为α,β的伽马分布Γ(c；α,β)，获取电气应力作用下电力系统元件寿命的样本数据；

B、对电气应力作用下电力系统元件寿命的累积分布函数、可靠度函数以及期望进行非精确推断，得到元件寿命在样本数据条件下的累积分布函数族、可靠度函数族和期望的非精确模型；

C、利用参数s构造矩形区域，先验参数α和β在区间[0,s]取值，构建电力系统元件寿命的累积分布函数族、可靠度函数族和期望区间值上下界的计算模型；

D、根据样本数据设定区间宽度ε，确定参数s的值，并结合样本容量N和的值分别计算电力系统元件寿命的累积分布函数族、可靠度函数族和期望的区间值上下界。

进一步地，所述步骤B中得到元件寿命在样本数据条件下的累积分布函数族的非精确模型，具体表示为：

其中，x_i为样本数据，N为样本容量。

进一步地，所述步骤B中得到元件寿命在样本数据条件下的可靠度函数族的非精确模型，具体表示为：