CN109472045A - 一种基于随机截尾数据的继电保护缺陷分布评估方法及系统 - Google Patents
一种基于随机截尾数据的继电保护缺陷分布评估方法及系统 Download PDFInfo
- Publication number
- CN109472045A CN109472045A CN201811123982.8A CN201811123982A CN109472045A CN 109472045 A CN109472045 A CN 109472045A CN 201811123982 A CN201811123982 A CN 201811123982A CN 109472045 A CN109472045 A CN 109472045A
- Authority
- CN
- China
- Prior art keywords
- defect
- distribution
- apparatus module
- parameter
- module
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000007547 defect Effects 0.000 title claims abstract description 212
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 40
- 230000001681 protective effect Effects 0.000 claims abstract description 40
- 230000001186 cumulative effect Effects 0.000 claims abstract description 26
- 238000005070 sampling Methods 0.000 claims description 15
- 238000005315 distribution function Methods 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 6
- 230000002045 lasting effect Effects 0.000 claims 2
- 239000004973 liquid crystal related substance Substances 0.000 description 5
- 241000208340 Araliaceae Species 0.000 description 4
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 4
- 235000003140 Panax quinquefolius Nutrition 0.000 description 4
- 235000008434 ginseng Nutrition 0.000 description 4
- 230000002950 deficient Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 241000607479 Yersinia pestis Species 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于随机截尾数据的继电保护缺陷分布评估方法及系统,包括:根据多个同类型继电保护装置的不同模块在观测期间的装置模块缺陷发生情况,将每个装置模块的无缺陷运行时间作为随机变量,未发生缺陷的装置模块的随机变量具有随机截尾数据特征;根据观测到的随机截尾数据,确定所述每类装置模块的两参数威尔分布的概率密度函数;所述两参数为形状参数β和尺度参数η;建立每类装置模块缺陷分布的似然函数,并根据最大似然函数法和观测到的随机截尾数据,计算每类装置模块的形状参数和尺度参数的极大似然估计值;根据每类装置模块的缺陷分布概率密度函数,建立装置整体缺陷的联合分布,并对继电保护装置的整体缺陷的累计分布进行评估。
Description
技术领域
本发明涉及继电保护技术领域,并且更具体地,涉及一种基于随机截尾数据的继电保护缺陷分布评估方法及系统。
背景技术
可靠性是指一个元件、设备或系统在预定时间内、规定的条件下完成规定功能的能力。对于继电保护而言,“规定功能”大致有两个方面:一是当被保护对象发生故障时保护系统要动作于跳闸,以隔离故障设备;二是当被保护对象处于异常运行状态时,保护系统要能发告警信号。保护设备良好的运行状态是确保保护正确动作的必要条件,但保护缺陷的发生会导致保护拒动、误动的发生,从而影响电网的安全稳定性。保护缺陷一般包括保护装置本体缺陷、二次回路缺陷、通信链路缺陷、定值整定错误等。继电保护缺陷发生模式主要分类如下:1)人为因素导致的保护缺陷,主要包括由于维护不良、误操作、误碰、误接线、误整定、接线设计不合理、未执行反事故措施等因素导致的保护缺陷;2)保护装置本体及相关设备故障导致的保护缺陷,主要包括元器件故障、电源损坏、通信通道问题等引起的缺陷。3)其他原因导致的缺陷,如接触不良、二次回路故障等。
随着设备运行年限的增加,保护装置本体组件老化、运维不良等导致的缺陷逐步突出。目前,对缺陷数据的记录主要包括:缺陷发生时间、缺陷部位、缺陷原因、缺陷严重程度、缺陷处理情况等。由于保护装置运行年限长,存在有些装置自投运到退出运行均未发生缺陷或者自投运到观测期未发生缺陷,缺陷数据的记录具有右截尾试验特性,对于准确评估缺陷分布情况,全面掌握缺陷发生规律对指导保护装置的检修具有重要意义。但由于保护装置缺陷失效为小概率事件,且无法对所有样本进行观测,失效数据具有右截尾特性;且各模块之间具有相对独立性,失效模型存在差异,不同模块的失效事件为独立分布。
因此,需要一种继电保护缺陷分布评估方法用于确定装置缺陷分布,为装置可靠性评估提供依据。
发明内容
本发明提出一种基于随机截尾数据的继电保护缺陷分布评估方法及系统,以解决如何对继电保护缺陷的分布进行评估的问题。
为了解决上述问题,根据本发明的一个方面,提供了一种基于随机截尾数据的继电保护缺陷分布评估方法,其特征在于,所述方法包括:
根据多个同类型继电保护装置的不同模块在观测期间的装置模块缺陷发生情况,将每个装置模块的无缺陷运行时间作为随机变量,其中,未发生缺陷的装置模块的随机变量具有随机截尾数据特征;
根据观测到的随机截尾数据,确定所述每类装置模块的两参数威尔分布的概率密度函数;其中,所述两参数为形状参数β和尺度参数η;
建立每类装置模块缺陷分布的似然函数,并根据最大似然函数法和观测到的随机截尾数据,计算每类装置模块的形状参数和尺度参数的极大似然估计值;
根据每类装置模块的缺陷分布概率密度函数,建立装置整体缺陷的联合分布,并对继电保护装置的整体缺陷的累计分布进行评估。
优选地,其中所述装置模块包括:电源模块、CPU模块、交流采样模块和液晶显示模块。
优选地,其中每类装置模块的缺陷概率分布相互独立,均满足两参数威布尔分布模型,确定的每个装置模块的两参数威尔分布的概率密度函数为:
其中,fi(t)为装置模块i的威尔分布的概率密度函数;β为形状参数;η为尺度参数;m为归一化系数。
优选地,其中建立的每类装置模块的缺陷分布的似然函数为:
其中,L(β,η)为每类装置模块的缺陷分布的似然函数;总体样本为n,观测r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;ci为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k;
其中,形状参数β和尺度参数η的极大似然估计通过对似然函数取极值,并求解下述方程得到:
和
优选地,其中利用如下方式计算每个装置模块的形状参数β和尺度参数η的极大似然估计值:
其中,r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;cj为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k。
优选地,其中利用如下方式对继电保护装置的整体缺陷的累计分布进行评估:
F(t)=1-[1-F1(t)]×[1-F2(t)]...[1-Fi-1(t)]×[1-Fi(t)],
其中,F(t)为继电保护装置的整体缺陷的累计分布函数;Fi为装置第i个模块缺陷累计分布。
根据本发明的另一个方面,提供了一种基于随机截尾数据的继电保护缺陷分布评估系统,其特征在于,所述系统包括:
随机截尾数据确定单元,用于根据多个同类型继电保护装置的不同模块在观测期间的装置模块缺陷发生情况,将每个装置模块的无缺陷运行时间作为随机变量,其中,未发生缺陷的装置模块的随机变量具有随机截尾数据特征;
概率密度函数确定单元,用于根据观测到的随机截尾数据,确定所述每类装置模块的两参数威尔分布的概率密度函数;其中,所述两参数为形状参数β和尺度参数η;
两参数的极大似然估计值计算单元,用于建立每类装置模块缺陷分布的似然函数,并根据最大似然函数法和观测到的随机截尾数据,计算每类装置模块的形状参数和尺度参数的极大似然估计值;
继电保护缺陷分布评估单元,用于根据每类装置模块的缺陷分布概率密度函数,建立装置整体缺陷的联合分布,并对继电保护装置的整体缺陷的累计分布进行评估。
优选地,其中所述装置模块包括:电源模块、CPU模块、交流采样模块和液晶显示模块。
优选地,其中每类装置模块的缺陷概率分布相互独立,均满足两参数威布尔分布模型,确定的每个装置模块的两参数威尔分布的概率密度函数为:
其中,fi(t)为装置模块i的威尔分布的概率密度函数;β为形状参数;η为尺度参数;m为归一化系数。
优选地,其中所述两参数为的极大似然估计值计算单元,建立的每类装置模块的缺陷分布的似然函数为:
其中,L(β,η)为每类装置模块的缺陷分布的似然函数;总体样本为n,观测r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;ci为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k;
其中,形状参数β和尺度参数η的极大似然估计通过对似然函数取极值,并求解下述方程得到:
和
优选地,其中在所述两参数为的极大似然估计值计算单元,利用如下方式计算每个装置模块的形状参数β和尺度参数η的极大似然估计值:
其中,r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;cj为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k。
优选地,其中在所述继电保护缺陷分布评估单元,利用如下方式对继电保护装置的整体缺陷的累计分布进行评估:
F(t)=1-[1-F1(t)]×[1-F2(t)]...[1-Fi-1(t)]×[1-Fi(t)],
其中,F(t)为继电保护装置的整体缺陷的累计分布函数;Fi为装置第i个模块缺陷累计分布。
本发明提供了一种基于随机截尾数据的继电保护缺陷分布评估方法及系统,包括:根据多个同类型继电保护装置的不同模块在观测期间的装置模块缺陷发生情况,将每个装置模块的无缺陷运行时间作为随机变量,其中,未发生缺陷的装置模块的随机变量具有随机截尾数据特征;根据观测到的随机截尾数据,确定所述每类装置模块的两参数威尔分布的概率密度函数;建立每类装置模块缺陷分布的似然函数,并根据最大似然函数法和观测到的随机截尾数据,计算每类装置模块的形状参数和尺度参数的极大似然估计值;根据每类装置模块的缺陷分布概率密度函数,建立装置整体缺陷的联合分布,并对继电保护装置的整体缺陷的累计分布进行评估。本发明采用极大似然估计法,基于保护装置的电源模块、CPU模块、交流采样模块及液晶显示模块存在独立分布特性的基础上,建立装置整体缺陷分布的联合分布模型,实现对装置各模块及整体缺陷概率分布模型的参数估计,对指导评估分析装置的可靠性,指导装置的检修周期确定具有重要意义,为电力系统继电保护缺陷分布、缺陷发生预测及缺陷评估提供有效的算法,完善了继保管理水平。
附图说明
通过参考下面的附图,可以更为完整地理解本发明的示例性实施方式:
图1为根据本发明实施方式的基于随机截尾数据的继电保护缺陷分布评估方法100的流程图;
图2为根据本发明实施方式的A型装置电源模块概率密度分布图;
图3为根据本发明实施方式的A型装置CPU模块概率密度分布图;
图4为根据本发明实施方式的A型装置交流采样模块概率密度分布图;
图5为根据本发明实施方式的A型装置液晶显示模块概率密度分布图;
图6为根据本发明实施方式的基于随机截尾数据的继电保护缺陷分布评估系统600的结构示意图。
具体实施方式
现在参考附图介绍本发明的示例性实施方式,然而,本发明可以用许多不同的形式来实施,并且不局限于此处描述的实施例,提供这些实施例是为了详尽地且完全地公开本发明,并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中,相同的单元/元件使用相同的附图标记。
除非另有说明,此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外,可以理解的是,以通常使用的词典限定的术语,应当被理解为与其相关领域的语境具有一致的含义,而不应该被理解为理想化的或过于正式的意义。
图1为根据本发明实施方式的基于随机截尾数据的继电保护缺陷分布评估方法100的流程图。如图1所示,本发明的实施方式提供的基于随机截尾数据的继电保护缺陷分布评估方法,采用极大似然估计法,基于保护装置的电源模块、CPU模块、交流采样模块及液晶显示模块存在独立分布特性的基础上,建立装置整体缺陷分布的联合分布模型,实现对装置各模块及整体缺陷概率分布模型的参数估计,对指导评估分析装置的可靠性,指导装置的检修周期确定具有重要意义,为电力系统继电保护缺陷分布、缺陷发生预测及缺陷评估提供有效的算法,完善了继保管理水平。本发明的实施方式提供的基于随机截尾数据的继电保护缺陷分布评估方法从步骤101处开始,在步骤101根据多个同类型继电保护装置的不同模块在观测期间的装置模块缺陷发生情况,将每个装置模块的无缺陷运行时间作为随机变量,其中,未发生缺陷的装置模块的随机变量具有随机截尾数据特征。
优选地,其中所述装置模块包括:电源模块、CPU模块、交流采样模块和液晶显示模块。
在本发明的实施方式中,选取保护装置的电源模块、CPU模块、交流采样模块、液晶显示模块的缺陷做为评估装置缺陷发生的具体部位。评估装置可以按型号、软件版本、装置厂家、运维单位、电压等级等不同条件选取。然后,将装置自投运至装置的电源模块、CPU模块、交流采样模块或者液晶显示模块发生缺陷的装置持续运行时间T作为随机变量,将自投运以来未发生缺陷装置的运行时间作为随机截尾数据。
优选地,在步骤102根据观测到的随机截尾数据,确定所述每类装置模块的两参数威尔分布的概率密度函数;其中,所述两参数为形状参数β和尺度参数η。
优选地,其中每类装置模块的缺陷概率分布相互独立,均满足两参数威布尔分布模型,确定的每个装置模块的两参数威尔分布的概率密度函数为:
其中,fi(t)为装置模块i的威尔分布的概率密度函数;β为形状参数;η为尺度参数;m为归一化系数。
优选地,在步骤103建立每类装置模块缺陷分布的似然函数,并根据最大似然函数法和观测到的随机截尾数据,计算每类装置模块的形状参数和尺度参数的极大似然估计值。
优选地,其中建立的每类装置模块的缺陷分布的似然函数为:
其中,L(β,η)为每类装置模块的缺陷分布的似然函数;总体样本为n,观测r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;ci为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k;
其中,形状参数β和尺度参数η的极大似然估计通过对似然函数取极值,并求解下述方程得到:
和
优选地,其中利用如下方式计算每个装置模块的形状参数β和尺度参数η的极大似然估计值:
其中,r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;cj为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k。
优选地,在步骤104根据每类装置模块的缺陷分布概率密度函数,建立装置整体缺陷的联合分布,并对继电保护装置的整体缺陷的累计分布进行评估。
优选地,其中利用如下方式对继电保护装置的整体缺陷的累计分布进行评估:
F(t)=1-[1-F1(t)]×[1-F2(t)]...[1-Fi-1(t)]×[1-Fi(t)],
其中,F(t)为继电保护装置的整体缺陷的累计分布函数;Fi为装置第i个模块缺陷累计分布。
以下具体举例说明本发明的实施方式
以下以某一型号保护装置的历史缺陷数据进行实例分析,具体步骤包括:
Step1:选取某型号A保护装置作为分析样本,分析最近15年型号A装置的历史缺陷数据,分析型号A保护装置电源模块、CPU模块、交流采样模块、液晶显示模块的缺失发生情况。表1为以电源模块的缺陷数据为例,对于在观测期间发生缺陷的数据,将自装置投运至缺陷发生的运行时间(装置无故障运行月数)作为随机变量,将观测期间未发生缺陷的装置运行年限作为随机截尾数据;CPU模块、交流采样模块、液晶显示模块的缺陷采用同样的处理方法。
表1A型号装置电源模块部分缺陷数据及截尾数据
Step2:根据观测到的随机截尾数据,确定所述每类装置模块的两参数威尔分布的概率密度函数。四类模块的缺陷概率分布相互独立,均满足两参数威布尔分布模型,四类模块的威布尔分布的概率密度函数为:
其中,i=1,2,3,4分别代表电源模块、CPU模块、采样模块、MMI模块。
Step3:建立四类模块的缺陷分布的似然函数为:
四类模块基于随机截尾数据的威布尔分布参数估计满足如下两方程,使用数值结算可以计算四模块的β,η两参数的极大似然估计值,进而求出四类模块缺陷的威布尔分布。利用如下方式计算每个装置模块的形状参数β和尺度参数η的极大似然估计值:
图2为根据本发明实施方式的A型装置电源模块概率密度分布图。如图2所示,A型保护装置电源模块的缺陷发生情况的概率密度分布f(t:β,η)=Weibull(1.99,1269.5),β=1.99,η=1269.5。
图3为根据本发明实施方式的A型装置CPU模块概率密度分布图。如图3所示,A型保护装置CPU模块的缺陷发生情况的概率密度分布f(t:β,η)=Weibull(1.64,1507.6),β=1.64,η=1507.6。
图4为根据本发明实施方式的A型装置交流采样模块概率密度分布图。如图4所示,A型保护装置交流采用模块的缺陷发生情况的概率密度分布f(t:β,η)=Weibull(1.6,1149.2),β=1.6,η=1149.2。
图5为根据本发明实施方式的A型装置液晶显示模块概率密度分布图。如图5所示,A型保护装置液晶显示模块的缺陷发生情况的概率密度分布f(t:β,η)=Weibull(1.71,2577.6),β=1.71,η=2577.6。
Step4:根据每类装置模块的缺陷分布概率密度函数,建立装置整体缺陷的联合分布,并对继电保护装置的整体缺陷的累计分布进行评估。装置整体缺陷的累计分布函数为F(t),装置整体分布为T=min{T1,T2,T3,T4}的分布,T1,T2,T3,T4分别为四类模块的随机变量。A型保护装置总体的联合分布为F(t)=1-[1-F1(t)][1-F2(t)][1-F3(t)][1-F4(t)]。装置整体缺陷的累计概率分布为:
图6为根据本发明实施方式的基于随机截尾数据的继电保护缺陷分布评估系统600的结构示意图。如图6所示,本发明的实施方式提供的基于随机截尾数据的继电保护缺陷分布评估系统600,包括:随机截尾数据确定单元601、概率密度函数确定单元602、两参数的极大似然估计值计算单元603和继电保护缺陷分布评估单元604。
优选地,所述随机截尾数据确定单元601,用于根据多个同类型继电保护装置的不同模块在观测期间的装置模块缺陷发生情况,将每个装置模块的无缺陷运行时间作为随机变量,其中,未发生缺陷的装置模块的随机变量具有随机截尾数据特征。
优选地,其中所述装置模块包括:电源模块、CPU模块、交流采样模块和液晶显示模块。
优选地,所述概率密度函数确定单元602,用于根据观测到的随机截尾数据,确定所述每类装置模块的两参数威尔分布的概率密度函数;其中,所述两参数为形状参数β和尺度参数η。
优选地,其中每类装置模块的缺陷概率分布相互独立,均满足两参数威布尔分布模型,确定的每个装置模块的两参数威尔分布的概率密度函数为:
其中,fi(t)为装置模块i的威尔分布的概率密度函数;β为形状参数;η为尺度参数;m为归一化系数。
优选地,所述两参数为的极大似然估计值计算单元603,用于建立每类装置模块缺陷分布的似然函数,并根据最大似然函数法和观测到的随机截尾数据,计算每类装置模块的形状参数和尺度参数的极大似然估计值。
优选地,其中所述两参数为的极大似然估计值计算单元603,建立的每类装置模块的缺陷分布的似然函数为:
其中,L(β,η)为每类装置模块的缺陷分布的似然函数;总体样本为n,观测r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;ci为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k;
其中,形状参数β和尺度参数η的极大似然估计通过对似然函数取极值,并求解下述方程得到:
和
优选地,其中在所述两参数为的极大似然估计值计算单元603,利用如下方式计算每个装置模块的形状参数β和尺度参数η的极大似然估计值:
其中,r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;cj为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k。
优选地,所述继电保护缺陷分布评估单元604,用于根据每类装置模块的缺陷分布概率密度函数,建立装置整体缺陷的联合分布,并对继电保护装置的整体缺陷的累计分布进行评估。
优选地,其中在所述继电保护缺陷分布评估单元,利用如下方式对继电保护装置的整体缺陷的累计分布进行评估:
F(t)=1-[1-F1(t)]×[1-F2(t)]...[1-Fi-1(t)]×[1-Fi(t)],
其中,F(t)为继电保护装置的整体缺陷的累计分布函数;Fi为装置第i个模块缺陷累计分布。
本发明的实施例的基于随机截尾数据的继电保护缺陷分布评估系统600与本发明的另一个实施例的基于随机截尾数据的继电保护缺陷分布评估方法100相对应,在此不再赘述。
已经通过参考少量实施方式描述了本发明。然而,本领域技术人员所公知的,正如附带的专利权利要求所限定的,除了本发明以上公开的其他的实施例等同地落在本发明的范围内。
通常地,在权利要求中使用的所有术语都根据他们在技术领域的通常含义被解释,除非在其中被另外明确地定义。所有的参考“一个/所述/该[装置、组件等]”都被开放地解释为所述装置、组件等中的至少一个实例,除非另外明确地说明。这里公开的任何方法的步骤都没必要以公开的准确的顺序运行,除非明确地说明。
Claims (12)
1.一种基于随机截尾数据的继电保护缺陷分布评估方法,其特征在于,所述方法包括:
根据多个同类型继电保护装置的不同模块在观测期间的装置模块缺陷发生情况,将每个装置模块的无缺陷运行时间作为随机变量,其中,未发生缺陷的装置模块的随机变量具有随机截尾数据特征;
根据观测到的随机截尾数据,确定所述每类装置模块的两参数威尔分布的概率密度函数;其中,所述两参数为形状参数β和尺度参数η;
建立每类装置模块缺陷分布的似然函数,并根据最大似然函数法和观测到的随机截尾数据,计算每类装置模块的形状参数和尺度参数的极大似然估计值;
根据每类装置模块的缺陷分布概率密度函数,建立装置整体缺陷的联合分布,并对继电保护装置的整体缺陷的累计分布进行评估。
2.根据权利要求1所述的方法,其特征在于,所述装置模块包括:电源模块、CPU模块、交流采样模块和液晶显示模块。
3.根据权利要求1所述的方法,其特征在于,每类装置模块的缺陷概率分布相互独立,均满足两参数威布尔分布模型,确定的每个装置模块的两参数威尔分布的概率密度函数为:
其中,fi(t)为装置模块i的威尔分布的概率密度函数;β为形状参数;η为尺度参数;m为归一化系数。
4.根据权利要求1所述的方法,其特征在于,建立的每类装置模块的缺陷分布的似然函数为:
其中,L(β,η)为每类装置模块的缺陷分布的似然函数;总体样本为n,观测r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;ci为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k;
其中,形状参数β和尺度参数η的极大似然估计通过对似然函数取极值,并求解下述方程得到:
和
5.根据权利要求1或4所述的方法,其特征在于,利用如下方式计算每个装置模块的形状参数β和尺度参数η的极大似然估计值:
其中,r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;cj为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k。
6.根据权利要求1所述的方法,其特征在于,利用如下方式对继电保护装置的整体缺陷的累计分布进行评估:
F(t)=1-[1-F1(t)]×[1-F2(t)]...[1-Fi-1(t)]×[1-Fi(t)],
其中,F(t)为继电保护装置的整体缺陷的累计分布函数;Fi为装置第i个模块缺陷累计分布。
7.一种基于随机截尾数据的继电保护缺陷分布评估系统,其特征在于,所述系统包括:
随机截尾数据确定单元,用于根据多个同类型继电保护装置的不同模块在观测期间的装置模块缺陷发生情况,将每个装置模块的无缺陷运行时间作为随机变量,其中,未发生缺陷的装置模块的随机变量具有随机截尾数据特征;
概率密度函数确定单元,用于根据观测到的随机截尾数据,确定所述每类装置模块的两参数威尔分布的概率密度函数;其中,所述两参数为形状参数β和尺度参数η;
两参数的极大似然估计值计算单元,用于建立每类装置模块缺陷分布的似然函数,并根据最大似然函数法和观测到的随机截尾数据,计算每类装置模块的形状参数和尺度参数的极大似然估计值;
继电保护缺陷分布评估单元,用于根据每类装置模块的缺陷分布概率密度函数,建立装置整体缺陷的联合分布,并对继电保护装置的整体缺陷的累计分布进行评估。
8.根据权利要求7所述的系统,其特征在于,所述装置模块包括:电源模块、CPU模块、交流采样模块和液晶显示模块。
9.根据权利要求7所述的系统,其特征在于,每类装置模块的缺陷概率分布相互独立,均满足两参数威布尔分布模型,确定的每个装置模块的两参数威尔分布的概率密度函数为:
其中,fi(t)为装置模块i的威尔分布的概率密度函数;β为形状参数;η为尺度参数;m为归一化系数。
10.根据权利要求7所述的系统,其特征在于,所述两参数为的极大似然估计值计算单元,建立的每类装置模块的缺陷分布的似然函数为:
其中,L(β,η)为每类装置模块的缺陷分布的似然函数;总体样本为n,观测r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;ci为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k;
其中,形状参数β和尺度参数η的极大似然估计通过对似然函数取极值,并求解下述方程得到:
和
11.根据权利要求7或10所述的系统,其特征在于,在所述两参数为的极大似然估计值计算单元,利用如下方式计算每个装置模块的形状参数β和尺度参数η的极大似然估计值:
其中,r为观测数据中对应模块发生缺陷的装置数量,ti为装置模块无缺陷运行时间,i=1,2…r,k为观测数据中对应模块未发生缺陷的装置数量;cj为未发生缺陷装置模块自投运至观测截止的持续运行时间,j=1,2…k,n=r+k。
12.根据权利要求7所述的系统,其特征在于,在所述继电保护缺陷分布评估单元,利用如下方式对继电保护装置的整体缺陷的累计分布进行评估:
F(t)=1-[1-F1(t)]×[1-F2(t)]...[1-Fi-1(t)]×[1-Fi(t)],
其中,F(t)为继电保护装置的整体缺陷的累计分布函数;Fi为装置第i个模块缺陷累计分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811123982.8A CN109472045A (zh) | 2018-09-26 | 2018-09-26 | 一种基于随机截尾数据的继电保护缺陷分布评估方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811123982.8A CN109472045A (zh) | 2018-09-26 | 2018-09-26 | 一种基于随机截尾数据的继电保护缺陷分布评估方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109472045A true CN109472045A (zh) | 2019-03-15 |
Family
ID=65663143
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811123982.8A Pending CN109472045A (zh) | 2018-09-26 | 2018-09-26 | 一种基于随机截尾数据的继电保护缺陷分布评估方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109472045A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111222582A (zh) * | 2020-01-13 | 2020-06-02 | 中国电力科学研究院有限公司 | 一种用于提取继电保护装置的缺陷模式的方法及系统 |
CN112630645A (zh) * | 2020-12-15 | 2021-04-09 | 中国电力科学研究院有限公司 | 一种确定继电保护设备的缺陷率的方法及系统 |
CN113221374A (zh) * | 2021-05-28 | 2021-08-06 | 哈尔滨工程大学 | 一种核电设备可靠性分析的样本数据生成方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103971024A (zh) * | 2014-05-26 | 2014-08-06 | 华北电力大学(保定) | 小样本失效数据下继电保护系统可靠性评估方法 |
-
2018
- 2018-09-26 CN CN201811123982.8A patent/CN109472045A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103971024A (zh) * | 2014-05-26 | 2014-08-06 | 华北电力大学(保定) | 小样本失效数据下继电保护系统可靠性评估方法 |
Non-Patent Citations (4)
Title |
---|
徐涛;李明贞;刘毅刚;周文俊;徐研;周承科;: "基于韦布尔分布的无失效数据参数估计的高压电缆寿命评估", 武汉大学学报(工学版), no. 02 * |
徐福明;林燕贞;: "基于广义证据理论线路继电保护的状态评估", 水电能源科学, no. 03 * |
李海波 等: "基于随机截尾数据下Weibull分布的参数极大似然估计与应用", 强度与环境, vol. 36, no. 4, pages 60 - 64 * |
黄烽等: ""继电保护装置运行年限评价的研究"", 《电力系统保护与控制》, vol. 45, no. 22, pages 159 - 161 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111222582A (zh) * | 2020-01-13 | 2020-06-02 | 中国电力科学研究院有限公司 | 一种用于提取继电保护装置的缺陷模式的方法及系统 |
CN111222582B (zh) * | 2020-01-13 | 2023-12-05 | 中国电力科学研究院有限公司 | 一种用于提取继电保护装置的缺陷模式的方法及系统 |
CN112630645A (zh) * | 2020-12-15 | 2021-04-09 | 中国电力科学研究院有限公司 | 一种确定继电保护设备的缺陷率的方法及系统 |
CN112630645B (zh) * | 2020-12-15 | 2024-05-24 | 中国电力科学研究院有限公司 | 一种确定继电保护设备的缺陷率的方法及系统 |
CN113221374A (zh) * | 2021-05-28 | 2021-08-06 | 哈尔滨工程大学 | 一种核电设备可靠性分析的样本数据生成方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109472045A (zh) | 一种基于随机截尾数据的继电保护缺陷分布评估方法及系统 | |
Cebrian et al. | Probabilistic assessment of financial losses in distribution network due to fault-induced process interruptions considering process immunity time | |
CN102945315B (zh) | 考虑软件失效和人为失效的全数字化继电保护可靠性评估方法 | |
CN103995172B (zh) | 变电站gis母线负荷电流在线监测方法 | |
CN103823972A (zh) | 一种继电保护设备状态评价与评估方法 | |
JP5621967B2 (ja) | 異常データ分析システム | |
CN103149475A (zh) | 电力设备故障诊断方法及系统 | |
CN102654539A (zh) | 一种电子式互感器运行状态评价方法 | |
EP4111191B1 (en) | Systems and methods for dissolved gas analysis | |
CN107239896A (zh) | 一种继电保护安全性评估方法 | |
CN107247198A (zh) | 一种配电设备故障预测方法及装置 | |
Sakthimohan et al. | Optimized Transmission Line Flaw Disclosure and Inkling System Based on IOT | |
CN106529763A (zh) | 一种配电系统运行分析方法和装置 | |
CN103869243B (zh) | 低压微机保护测控装置操作板继电器智能校验方法 | |
JP2010193625A (ja) | 入力電気量の異常兆候判定装置 | |
JP2010216958A (ja) | 接地装置の試験方法及び装置 | |
CN107918690A (zh) | 一种基于bp神经网络算法的继电保护设备状态评估方法 | |
WO2016136391A1 (ja) | 事故点標定装置及び方法、電力系統監視システム、設備計画支援システム | |
US10374548B2 (en) | Inspection support apparatus and control method thereof, inspection system, and control program | |
CN105891637A (zh) | 一种智能变电站二次设备测试方法及系统 | |
CN110516960A (zh) | 一种变电站继电保护设备的可靠性指标定量计算方法 | |
Subramani et al. | Simulation technique for voltage stability Analysis and contingency ranking in power systems | |
JP5018474B2 (ja) | 半導体デバイス試験装置及び半導体デバイス試験方法 | |
Wang et al. | An application example for on-line DSA systems | |
Salim et al. | Identifying severe loading condition during the events of cascading outage considering the effects of protection system hidden failure |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190315 |
|
RJ01 | Rejection of invention patent application after publication |