011引则采用不同溫度Θι，θ2作为基准恒定溫度Θη时，此直定溫度Θ据行时间t时的h%分 别为的：
[0119] 在相同条件下，变压器的故障概率相同，即1/〇山01，0〇=11%山02，0〇，则有：
[0120]
[012。 将公式13再次变化为exp(丫 3Z3+丫 qZq)的形式，则有：
[0122]
[0123] 步骤(4).参数估计
[0124] 模型的参数估计采用极大似然估计。一般形式的极大似然函数为：
[0125]
[01%] 公式15中参数定义如公式6。
[0127]相应似然函数为：
[0132] 公式17中样本总体个数为n，失效个数为q，ti表示第i个失效个体在基准溫度Θη下 的等效时间，协变量条件的个数为P，s为设备个体在Θη下的等效检测时间点。
[0133] 则油浸式变压器的极大似然函数的对数形式为：
[0134]
[0135] 将公式18参数β、丫求一阶导并置为零，可得一组超越方程，之后可采用相关算法 求解获得貧、7。
[0136] 步骤(5).故障概率求解
[0137] 根据参数估计的结果，代入油浸式变压器等效运行时间Teq、健康状态、天气状况、 实时热点溫度，即可计算变压器实时故障概率。
【主权项】
1. 一种面向风险评估的输变电设备故障概率模型的建模方法，其特征在于包括如下步 骤： 步骤(1).输变电设备风险因素分析与风险指标筛选，具体如下： 1-1.输变电设备风险因素分析 通过深入分析输变电设备的老化失效机理和缺陷劣化机理，基于设备实时状态监测、 微气象及自然灾害预报信息，针对不同设备个体分析全寿命周期信息和非周期随机信息， 进而分析运行风险因素与设备故障停运的内在关联； 所述的全寿命周期信息包括物理表征、老化状态、健康指数； 所述的非周期随机信息包括微气象、自然灾害、负载状况； 1- 2.风险指标筛选 利用大数据平台以及大数据筛选分析方法，将各项风险因素整合为几项风险指标，整 合后的风险指标包括老化失效指标、健康状态指标、天气状况指标以及包括负载状况指标 在内的设备相关运行工况指标； 步骤(2).模型基础分析 2- 1.比例风险模型建模 PHM在电力系统输变电设备的可靠性评估和检修计划中有广泛的应用；根据模型定义 可知故障概率函数为： h(t；Z)=h〇(t)exp( γ Ζ)公式 1 公式1中，t为当前时刻，ho (t)为基准故障概率函数，Ζ为反映设备不同状态的协变量向 量，γ为协变量向量参数，exp( γ Ζ)为协变量连接函数； ΡΗΜ认为状态指标与设备风险之间具有乘法效应;通过该模型能有效建立设备各类状 态指标和故障概率之间的联系，从而获得设备的实时故障概率； 所述的状态指标即为步骤1中整合后的风险指标； 2- 2.温度-负载状况相依模型建模 设备负载状况对故障概率的影响十分显著；同时变压器、输电线路等输变电设备时刻 与外界环境发生热交换，需建立相应温度-负载状况相依模型，如变压器的热点温度计算模 型、输电线路的导体温度模型等； 步骤(3).故障概率模型建模 3- 1.基准故障概率函数h〇(t)的建模 ho(t)作为基准故障概率函数，可用来表征输变电设备的老化失效情况;威布尔分布充 分契合浴盆曲线变化趋势，已广泛应用于输变电设备的失效建模中，因此本模型的基准故 障概率函数采用威布尔分布：公式2中，β为形状参数，τι为比例参数;若忽略浴盆曲线中后期偶然故障的影响，可认为 β>〇，此时故障概率随时间的增长而上升； 3-1-1.比例参数η的推导 η在威布尔分布中表征寿命期望，可通过分析设备老化失效机理推导； 3-1-2.等效运行时间1^的推导 公式2中的t值为恒定基准温度ΘΗ下的等效运行时间Teq;需要将波动的设备温度曲线对 应的役龄时间折算为恒定ΘΗ下的等效运行时间；可按导体温度值将运行时间划分成η个小 区间t，每个小区间内的导体温度θ Ηι恒定，累加可得等效服役时间； 3-1-3.结合η与Teq推导基准故障概率函数h〇( t)表达式； 3-2.协变量连接函数exp( γ Z)的建模 exp( γ Ζ)是协变量连接函数，在ΡΗΜ中反映不同状态协变量对设备故障概率的影响；当 协变量个数为Ρ时，对应的协变量连接函数为：除老化失效外，影响输变电设备故障停运的因素还有很多，归结起来主要有三类:设备 自身的健康状态指标、天气状况指标以及负载状况指标；同时不同的设备可能受除这三类 风险因素外的其他风险因素的影响，所述的负载状况指标是设备相关运行工况指标的一个 因素;因此协变量连接函数中协变量应包含这几方面风险因素影响，则协变量连接函数为：公式4中，Ζι、Ζ2、Ζ3、γ 1、γ 2、γ 3分别表不健康状态、天气状况、负载状况以及相应的协 变量参数;Zk、γ k分别表示其他风险因素以及相应的协变量参数； 3-2-1.设备自身的健康状态建模 将设备状态分为正常、注意、异常、严重四类； 3-2-2.设备所处天气状况建模 仿照健康状态的建模标准对天气状态进行建模，根据国家气象局规定的台风、暴雨、暴 雪、寒潮、大风、高温、雷电、霜冻、大雾、道路结冰、冰雹、干旱、霾等13项天气预警信号中选 取对输变电设备故障概率影响较大的雷电、覆冰、降雨、大风、气温、台风、冰雹、降雪8项天 气风险因素，根据气象信息，将天气状况分为正常天气、恶劣天气、灾变天气三类； 3-2-3.设备负载状况建模 设备负载状况通过设备实时温度落实到故障概率上；设备实时温度能够通过温度-负 载状况相依模型以及微气象信息计算而得;基准故障概率函数以及协变量连接函数均为设 备实时温度的因变量，因此协变量连接函数中设备负载状况无法像健康状态和天气状况一 样进行简单建模，需进行单独分析;根据PHM的定义，状态变量与设备风险之间具有乘法效 应，此外，在一定环境条件下可引入温度影响函数g(9)来代替负载状况协变量的作用函数， 则协变量连接函数为：3-2-4.其他风险因素建模 仿照设备健康状态和天气状况的建模方式以其他风险因素为依据，对不同的设备个体 进行状态分级； 3-2-5.温度影响函数g(0)建模 温度影响函数g(9)是实时设备温度01和设备基准恒定温度ΘΗ共同作用的结果;在相同 健康状态和天气状况下，采用不同温度91，02作为基准恒定温度ΘΗ时，以恒定温度0 h运行时 间t时的故障概率是相同的，以此为依据进行温度影响函数g(0)的计算;求得g(0)后将其再 次变化为exp( γ 3Ζ3+ γ qZq)的形式代入协变量连接函数，则协变量连接函数为：步骤(4).参数估计 模型的参数估计采用极大似然估计;一般形式的极大似然函数为：公式7中样本总体个数为η;是否截尾的状态量为δ,，截尾数据用〇表示，失效数据用1表 示;观测到的寿命时间的样本集为Q，其中失效个数为q，R为可靠度函数，f为故障概率密度 函数;h为故障概率函数;均表示时间； 则相应似然函数为：公 式8 取对数得公式9中样本总体个数为η，失效个数为q，t表示第i个失效个体在基准恒定温度ΘΗ下的 等效时间，协变量条件的个数为P，s为设备个体在ΘΗ下的等效检测时间点； 将式中参数β、γ、η求一阶导并置为零，可得一组超越方程，其中η已知，可采用相关算 法求解获得/L f; 步骤(5).故障概率求解 根据步骤4求解获得的参数|、f，估计结果，代入设备等效运行时间Tm、健康状态、天 气状况、负载状况等风险指标，即可计算设备故障概率。
【专利摘要】本发明公开了一种面向风险评估的输变电设备故障概率模型的建模方法。本发明基于设备实时状态监测、微气象及自然灾害预报信息，针对不同设备个体分析其全寿命周期信息和非周期随机信息，利用大数据筛选分析方法，将各项风险因素整合为几项风险指标作为基于比例风险模型的输变电设备故障概率模型的输入，进而量化各风险因素对故障趋势的影响。本发明通过风险因素分析与风险指标筛选、模型基础分析、模型建立、参数估计、故障概率计算等手段，为电网的风险评估提供理论基础的综合故障概率模型。本发明通过模型计算结果，结合负载能力评估、故障预测等结果，可分析计算设备缺失所引发的系统潜在损失和相关风险评估指标。
【IPC分类】G06Q50/06, G06Q10/06
【公开号】CN105631578
【申请号】CN201510908694
【发明人】何迪, 章禹, 杨祎, 郭创新, 张金江, 周加斌, 黄锐, 张行, 夏伊乔
【申请人】浙江大学, 国网山东省电力公司
【公开日】2016年6月1日
【申请日】2015年12月10日