CN109583520B - 一种云模型与遗传算法优化支持向量机的状态评估方法 - Google Patents

一种云模型与遗传算法优化支持向量机的状态评估方法 Download PDF

Info

Publication number
CN109583520B
CN109583520B CN201811608793.XA CN201811608793A CN109583520B CN 109583520 B CN109583520 B CN 109583520B CN 201811608793 A CN201811608793 A CN 201811608793A CN 109583520 B CN109583520 B CN 109583520B
Authority
CN
China
Prior art keywords
sample
state
data
protection device
classification
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.)
Active
Application number
CN201811608793.XA
Other languages
English (en)
Other versions
CN109583520A (zh
Inventor
张学敏
王斌
施迎春
王文林
党军朋
杨永旭
马智鹏
唐一恒
李雷
刘祺
韩宗延
张志强
周洪胜
戴伟康
白建林
杜林强
矣林飞
许鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yuxi Power Supply Bureau of Yunnan Power Grid Co Ltd
Original Assignee
Yuxi Power Supply Bureau of Yunnan Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yuxi Power Supply Bureau of Yunnan Power Grid Co Ltd filed Critical Yuxi Power Supply Bureau of Yunnan Power Grid Co Ltd
Priority to CN201811608793.XA priority Critical patent/CN109583520B/zh
Publication of CN109583520A publication Critical patent/CN109583520A/zh
Application granted granted Critical
Publication of CN109583520B publication Critical patent/CN109583520B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • G06F18/241Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2411Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/12Computing arrangements based on biological models using genetic models
    • G06N3/126Evolutionary algorithms, e.g. genetic algorithms or genetic programming
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Biophysics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Genetics & Genomics (AREA)
  • Biomedical Technology (AREA)
  • Computational Linguistics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种云模型与遗传算法优化支持向量机的状态评估方法。包括:首先根据采集单元采集的运行状态数据经数据预处理后形成样本数据库,选择径向基核函数作为支持向量机的核函数,采用遗传算法优化支持向量机的核函数参数及错误惩罚因子,建立遗传算法优化的支持向量机状态评估模型,将评估问题转化为一个分类问题,对保护装置的运行状态进行评估,利用云模型的随机性和稳定的趋向性来实现健康度与评语域之间的不确定性映射,使得评估结果更加符合实际情况。本发明极大的提升了状态检修工作的效率,使得运维检修人员及时掌握保护装置运行状态,防止设备出现安全事故,保证供电的安全可靠,且评估准确率显著提高。

Description

一种云模型与遗传算法优化支持向量机的状态评估方法
技术领域
本发明属于电力行业并涉及继电保护装置,具体涉及一种云模型与遗传算法优化支持向量机的状态评估方法。
背景技术
目前,继电保护装置通常采用定期检修模式,该方法可能会存在“检修不足、检修过剩”,导致设备状态不确定。我国电力系统发展迅速,社会对供电质量以及可靠性提出了更高的要求,保证智能变电站的继电保护装置安全、稳定运行至关重要。而状态检修以设备状态为依据,通过对设备的连续性观测并且综合其他各种因素来识别设备已有的或潜在的劣化现象,然后对设备的各个状态量进行合理的评估,最后进行预测评估确定最佳检修时间。状态评估是状态检修的基础,只有准确、合理的评估设备运行状态才能制定合理的状态检修策略,及时有效的开展检修工作,为运维检修工作人员提供参考,实现智能变电站安全、可靠运行。
发明内容
为了实现上述目的,本发明提出了一种云模型与遗传算法优化支持向量机的状态评估方法。
本发明的具体技术方案为一种云模型与遗传算法优化支持向量机的状态评估方法,具体包括以下步骤:
步骤1:根据继电保护装置运行状态信息,选取工作电压、CPU温度、绝缘性能、设备故障次数、家族缺陷率、异常告警率、断路器不正确动作次数作为支持向量机的输入特征向量;
步骤2:分别对工作电压、CPU温度、绝缘性能、设备故障次数、家族缺陷率、异常告警率、断路器不正确动作次数进行数据预处理,得到训练样本数据;
步骤3:将训练样本数据通过人工标记方法得到标记后数据预处理样本,将标记后数据预处理样本作为支持向量机的输入,支持向量机的核函数为径向基核函数,使用遗传算法对支持向量机的核函数参数、错误惩罚因子进行参数寻优,得到分类效果最好的参数取值以构建优化后支持向量机;
步骤4:将训练样本数据通过人工标记方法得到标记后数据预处理样本,通过遗传算法优化后支持向量机对标记后数据预处理样本进行分类训练,得到失效状态样本分类边界、最优状态分类面、良好状态样本分类边界;
步骤5:将测试样本数据根据步骤1-4得到测试样本数据的输入向量,计算测试样本数据点到最优分类超平面的距离,根据测试样本数据点到最优分类超平面的距离判断测试样本数据的状态;
步骤6:利用云模型产生云滴的随机性和稳定的趋向性来模拟不同专家对测试样本数据点到良好状态样本分类边界的距离的不同评估值,实现设备的评估值到评语域的不确定转化。
作为优选,步骤1中所述工作电压ai为第i个时间点的继电保护装置运行状态的工作电压;
步骤1中所述CPU温度bi为第i个时间点的继电保护装置运行状态的CPU温度;
步骤1中所述绝缘性能ci为第i个时间点的继电保护装置运行状态的绝缘性能;
步骤1中所述设备故障次数di为第i个时间点的继电保护装置运行状态的设备故障次数;
步骤1中所述家族缺陷率ei为第i个时间点的继电保护装置运行状态的家族缺陷率;
步骤1中所述异常告警率fi为第i个时间点的继电保护装置运行状态的异常告警率;
步骤1中所述断路器不正确动作次数gi为第i个时间点的继电保护装置运行状态的断路器不正确动作次数;
i∈[0,M],M为继电保护装置运行时间;
作为优选,步骤2中所述对工作电压进行数据预处理为:
Figure BDA0001924219460000021
其中,M为继电保护装置运行时间,N为数据预处理后的样本数量,ai为第i个时间点的继电保护装置运行状态的工作电压,ai *为第i个时间点的保护装置运行状态的工作电压安全稳定阈值,aj *为数据预处理后第j个样本的的工作电压,其为0-1之间的一个数值;
步骤2中所述对CPU温度进行数据预处理为:
Figure BDA0001924219460000022
其中,bi为第i个时间点的保护装置运行状态的CPU温度,bi *为第i个时间点的保护装置运行状态的CPU温度安全稳定阈值,bj *为数据预处理后第j个样本的CPU温度,其为0-1之间的一个数值;
步骤2中所述对绝缘性能进行数据预处理为:
Figure BDA0001924219460000023
其中,ci为第i个时间点的保护装置运行状态的绝缘性能,ci *为第i个时间点的保护装置运行状态的绝缘性能安全稳定阈值,cj *为数据预处理后第j个样本的绝缘性能,其为0-1之间的一个数值;
步骤2中所述对设备故障次数进行数据预处理为:
Figure BDA0001924219460000024
其中,di为第i个时间点的保护装置运行状态的设备故障次数,di *为第i个时间点的保护装置运行状态的设备故障次数安全稳定阈值,dj *为数据预处理后第j个样本的设备故障次数,其为0-1之间的一个数值;
步骤2中所述对家族缺陷率进行数据预处理为:
Figure BDA0001924219460000031
其中,ei为第i个时间点的保护装置运行状态的家族缺陷率,ei *为第i个时间点的保护装置运行状态的家族缺陷率安全稳定阈值,ej *为数据预处理后第j个样本的家族缺陷率,其为0-1之间的一个数值;
步骤2中所述对异常告警率进行数据预处理为:
Figure BDA0001924219460000032
其中,fi为第i个时间点的保护装置运行状态的异常告警率,fi *为第i个时间点的保护装置运行状态的异常告警率安全稳定阈值,fj *为数据预处理后第j个样本的异常告警率,其为0-1之间的一个数值;
步骤2中所述对不正确动作次数进行数据预处理为:
Figure BDA0001924219460000033
其中,gi为第i个时间点的保护装置运行状态的不正确动作次数,gi *为第i个时间点的保护装置运行状态的不正确动作次数安全稳定阈值,gj *为数据预处理后第j个样本的不正确动作次数,其为0-1之间的一个数值;
步骤2中所述训练样本数据为:
数据处理后第j个样本的工作电压aj *,数据处理后第j个样本的CPU温度bj *,数据处理后j个样本的绝缘性能cj *,数据处理后第j个样本的设备故障次数dj *,数据处理后第j个样本的家族缺陷率ej *,数据处理后第j个样本的异常告警率fj *,数据处理后第j个样本的断路器不正确动作次数gj *
作为优选,步骤3中所述将训练样本数据通过人工标记方法得到标记后数据预处理样本为:
(xj,yj),j∈[1,N]
xj=(aj *,bj *,cj *,dj *,ej *,fj *,gj *)T
yj∈{-1,1}
其中,N为标记后数据预处理样本的数量,(xj,yj)为第j组标记后数据预处理样本点,xj为第j组标记后数据预处理样本的输入向量;
若通过人工标记判断xj在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,则yj=1,若通过人工标记判断xj在在保护装置运行状态出现故障的情况下获得的样本数据,则yj=-1;
yj∈{-1,1}为第j组标记后数据预处理样本的输出结果,yj=1代表第j个样本的输入向量xj为状态良好的样本数据,状态良好的样本数据是在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,yj=-1代表第j个样本的输入向量xj为状态失效的样本数据,状态失效的样本数据是在保护装置运行状态出现故障的情况下获得的样本数据;
aj *为数据处理后第j个样本的工作电压,bj *为数据处理后第j个样本的CPU温度,cj *为数据处理后第j个样本的绝缘性能,dj *为数据处理后第j个样本的设备故障次数,ej *为数据处理后第j个样本的家族缺陷率,fj *为数据处理后第j个样本的异常告警率,gj *为数据处理后第j个样本的断路器不正确动作次数;
步骤3中所述支持向量机的核函数为径向基核函数,径向基核函数模型为:
Figure BDA0001924219460000041
其中,x为支持向量机的输入向量,x∈{x1,x2,...,xN},xi *为yi=1或yi=-1任意的一个支持向量,支持向量定义为正好落在良好状态分类分类边界上的向量或者失效状态分类分类边界上的向量;
步骤3中所述参数寻优具体为:
遗传算法中的输入数据为支持向量机的核函数参数C、错误惩罚因子δ、以及经过支持向量机的预测分类结果yp *,遗传算法中的优化目标定义为测试样本的平均相对百分比误差;
支持向量机的输入为第j个样本标记后数据预处理样本的输入向量xj,支持向量机的输出为第j个样本标记后数据预处理样本的输出结果yj
支持向量机的核函数类型、核函数参数C、错误惩罚因子δ均会对支持向量机的预测分类结果yp *产生影响;
遗传算法的优化目标模型:
Figure BDA0001924219460000042
其中,Q表示预测数量,Q∈[1,N],yp *表示支持向量机的预测分类输出结果,yp(yp∈{-1,1})表示为第j组标记后数据预处理样本的输出结果;
当遗传算法的优化目标模型函数达到最小值时,即认为此时的分类效果最好,取此时的支持向量机参数为最优参数;
通过核函数模型、最优核函数参数C*、最优错误惩罚因子δ*构建遗传算法优化后支持向量机;
作为优选,步骤4中所述标记后数据预处理样本为:
(xj,yj),j∈[1,N]
xj=(aj *,bj *,cj *,dj *,ej *,fj *,gj *)T
yj∈{-1,1}
其中,N为标记后数据预处理样本的数量,(xj,yj)为第j组标记后数据预处理样本点,xj为第j组标记后数据预处理样本的输入向量;
若通过人工标记判断xj在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,则yj=1,若通过人工标记判断xj在在保护装置运行状态出现故障的情况下获得的样本数据,则yj=-1;
yj∈{-1,1}为第j组标记后数据预处理样本的输出结果,yj=1代表第j个样本的输入向量xj为状态良好的样本数据,状态良好的样本数据是在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,yj=-1代表第j个样本的输入向量xj为状态失效的样本数据,状态失效的样本数据是在保护装置运行状态出现故障的情况下获得的样本数据;
aj *为数据处理后第j个样本的工作电压,bj *为数据处理后第j个样本的CPU温度,cj *为数据处理后第j个样本的绝缘性能,dj *为数据处理后第j个样本的设备故障次数,ej *为数据处理后第j个样本的家族缺陷率,fj *为数据处理后第j个样本的异常告警率,gj *为数据处理后第j个样本的断路器不正确动作次数;
步骤4中所述对标记后数据预处理样本进行分类训练为:
标记后数据预处理样本(xj,yj),j∈[1,N]可以被分类超平面wT·xj+b=0分开,其中,w=(w1,w2,w3,...,wN)T为分类超平面的法向量,b为代表分类超平面到原点之间的距离;
(xj,yj),j∈[1,N]中与分类超平面最近的距离之和称为分类间隔,分类间隔等于2/||w||,分类间隔最大时的超平面为最优分类超平面;
要使分类间隔最大即是使||w||2最小,因此可以转化为如下的约束优化问题:
Figure BDA0001924219460000061
s.t.yj(wT·xj+b)≥1
其中,w=(w1,w2,w3,...,wN)T为分类超平面的法向量,b为代表分类超平面到原点之间的距离,xj为第j组数据预处理样本的输入向量,yj为第j组数据预处理样本的输出结果;
这个约束优化问题可以通过构造拉格朗日函数求解,求解拉格朗日函数的鞍点,引入拉格朗日因子λj≥0,构造拉格朗日函数如下:
Figure BDA0001924219460000062
其中,λj≥0,j∈[1,N],为拉格朗日因子。
依据拉格朗日对偶理论将
Figure BDA0001924219460000063
转化为对偶问题,即:
Figure BDA0001924219460000064
Figure BDA0001924219460000065
应用二次规划方法可求解,求解得到的最优解α*=[λ1 *2 *,...,λN *]T,则可以得到最优的w*,b*
Figure BDA0001924219460000066
Figure BDA0001924219460000067
其中,xr、xs为yr=1,ys=1或yr=-1,ys=-1任意的一对支持向量,支持向量定义为正好落在良好状态分类分类边界上的向量或者失效状态分类分类边界上的向量;
可知w*T·xj+b*=0为最优分类超平面,w*T·xj+b*=+1为良好状态样本分类边界,w*T·xj+b*=-1为失效状态样本分类边界;
作为优选,步骤5中所述测试样本数据到最优分类面的距离为
Figure BDA0001924219460000068
其中,d为测试样本数据到最优分类边界面的距离,x为步骤5中所述测试样本数据的输入向量,其中w*=(w*1,w*2,w*3,...,w*7)T为平面的法向量,b*为一个实数,代表平面到原点之间的距离,w*、b*为步骤四中求得的最优值:
Figure BDA0001924219460000069
Figure BDA00019242194600000610
其中,xr、xs为两个类别中任意的一对支持向量,xi为支持向量机的输入向量,yi为支持向量机的输出结果,λi *为拉格朗日因子;
步骤5中所述根据测试样本数据点到最优分类超平面的距离判断测试样本数据的状态为:
若d>1则测试样本数据属于良好状态;
若0≤d<1则测试样本数据属于注意状态,测试样本数据到良好状态样本分类边界的距离为d'=1-d;
若-1≤d≤0则测试样本数据属于异常状态,测试样本数据到良好状态分类边界面的距离为d'=1-d;
若d<-1则测试样本数据属于失效状态;
作为优选,步骤6中所述云模型为:构建云模型,即期望Er、熵En、超熵Hq
其中,期望Er表示云模型的重心位置,反映了定性概念Q的中心值,熵En表示模糊性与随机性相关联程度的大小,超熵Hq为熵的熵,间接反映了云模型的离散程度;
步骤6中所述模拟不同专家对测试样本数据点到良好状态样本分类边界的距离的不同评估值为:
定义健康度H为:保护装置健康状态的度量,其值越大表示健康状态越好;健康度H用d'对应的云模型的期望值表示,云模型的数学期望表达式为:
H(d')=exp[-(d'-Er)2/2En 2]
其中,d'为测试样本数据点到良好状态分类边界面的距离,Er为云模型期望值,En为云模型熵;
云模型的重心Er=0,En=2/3,超熵Hq依据经验取值Hq=0.1;
测试样本数据点到最优分类超平面的距离为:
Figure BDA0001924219460000071
若d>1,测试样本数据点属于良好状态,可以不经过云模型直接对设备状态进行判断为良好状态;
若d<-1,测试样本数据点属于失效状态,可以不经过云模型直接对设备状态进行判断为失效状态;
若d∈[-1,1],需经过云模型的模糊映射:
将d'作为云模型的期望值,随机产生K个随机数,每个随机数有个相对应的云模型隶属度Vi,Vi∈[0,1],i=1,2,...,K;
通过比较Vi与H(d')的大小,H(d')=exp[-(d'-Er)2/2En 2],得到测试样本处于注意状态的数量和异常状态的数量;
当0<Vi<H(d'),i=1,2,...,K时,统计此时Vi的个数计数为N,N<K,此时Vi的状态为异常状态;
当H(d')<Vi<1,i=1,2,...,K时,统计此时Vi的个数计数为K-N,此时Vi的状态为注意状态;
分别统计健康状态处于注意状态的数量K-N和异常状态的数量N,取数量最大的状态为最终的状态,并给出可信度R如下:
Figure BDA0001924219460000072
由支持向量机分类结果可知d'=1为良好状态和注意状态的分类边界,d'为测试样本数据点到良好状态分类边界面的距离;
步骤6中所述实现设备的评估值到评语域的不确定转化:
由云模型的数学期望表达式确定其值为H(d'),从而确定由健康度H向评语域转化;
将得到的健康度转化为百分制:规定良好状态的样本健康度为100,失效状态的样本健康度为0,
注意状态的样本转化公式如下:
Figure BDA0001924219460000081
异常状态的样本转化公式如下:
Figure BDA0001924219460000082
其中,H为云模型的数学期望值,H'为健康状态转化为百分制后的得分。得到继电保护装置最终状态评估结果;
测试样本数据点到最优分类超平面的距离为:
Figure BDA0001924219460000083
若d>1,测试样本数据点属于良好状态,可以不经过云模型直接对设备状态进行判断为良好状态;
若d<-1,测试样本数据点属于失效状态,可以不经过云模型直接对设备状态进行判断为失效状态;
若d∈[-1,1],需经过云模型的模糊映射:
将d'作为云模型的期望值,随机产生K个随机数,每个随机数有个相对应的云模型隶属度Vi,Vi∈[0,1],i=1,2,...,K;
通过比较Vi与H(d')的大小,H(d')=exp[-(d'-Er)2/2En 2],得到测试样本处于注意状态的数量和异常状态的数量;
当0<Vi<H(d'),i=1,2,...,K时,统计此时Vi的个数计数为N,N<K,此时Vi的状态为异常状态;
当H(d')<Vi<1,i=1,2,...,K时,统计此时Vi的个数计数为K-N,此时Vi的状态为注意状态;
分别统计健康状态处于注意状态的数量K-N和异常状态的数量N,取数量最大的状态为最终的状态。
本发明优点在于:
本发明实现了对继电保护装置的运行状态进行准确状态评估,从而极大的提升了状态检修工作的效率,使得运维检修人员及时掌握保护装置运行状态,防止设备出现安全事故,保证供电的安全可靠。同时,本发明比层次分析法等传统状态评估方法有着明显的优势,评估准确率显著提高。
本发明通过对继电保护装置运行状态数据进行分析,达到对继电保护装置运行状态进行评估的目的,做到“应修必修”,防患于未然,为运维检修人员提供科学检修依据,保证智能变电站安全稳定运行。
附图说明
图1:继电保护装置状态指标评估体系;
图2遗传算法优化支持向量机参数流程图;
图3:支持向量机分类示意图;
图4:本发明的方法流程图。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明,下面结合附图及实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
下面结合图1至图4介绍本发明的实施方式,具体为:
步骤1:根据继电保护装置运行状态信息,选取工作电压、CPU温度、绝缘性能、设备故障次数、家族缺陷率、异常告警率、断路器不正确动作次数作为支持向量机的输入特征向量;
步骤1中所述工作电压ai为第i个时间点的继电保护装置运行状态的工作电压;
步骤1中所述CPU温度bi为第i个时间点的继电保护装置运行状态的CPU温度;
步骤1中所述绝缘性能ci为第i个时间点的继电保护装置运行状态的绝缘性能;
步骤1中所述设备故障次数di为第i个时间点的继电保护装置运行状态的设备故障次数;
步骤1中所述家族缺陷率ei为第i个时间点的继电保护装置运行状态的家族缺陷率;
步骤1中所述异常告警率fi为第i个时间点的继电保护装置运行状态的异常告警率;
步骤1中所述断路器不正确动作次数gi为第i个时间点的继电保护装置运行状态的断路器不正确动作次数;
i∈[0,M],M为继电保护装置运行时间;
步骤2:分别对工作电压、CPU温度、绝缘性能、设备故障次数、家族缺陷率、异常告警率、断路器不正确动作次数进行数据预处理,得到训练样本数据;
步骤2中所述对工作电压进行数据预处理为:
Figure BDA0001924219460000091
其中,M为继电保护装置运行时间,N为数据预处理后的样本数量,ai为第i个时间点的继电保护装置运行状态的工作电压,ai *为第i个时间点的保护装置运行状态的工作电压安全稳定阈值,aj *为数据预处理后第j个样本的的工作电压,其为0-1之间的一个数值;
步骤2中所述对CPU温度进行数据预处理为:
Figure BDA0001924219460000092
其中,bi为第i个时间点的保护装置运行状态的CPU温度,bi *为第i个时间点的保护装置运行状态的CPU温度安全稳定阈值,bj *为数据预处理后第j个样本的CPU温度,其为0-1之间的一个数值;
步骤2中所述对绝缘性能进行数据预处理为:
Figure BDA0001924219460000101
其中,ci为第i个时间点的保护装置运行状态的绝缘性能,ci *为第i个时间点的保护装置运行状态的绝缘性能安全稳定阈值,cj *为数据预处理后第j个样本的绝缘性能,其为0-1之间的一个数值;
步骤2中所述对设备故障次数进行数据预处理为:
Figure BDA0001924219460000102
其中,di为第i个时间点的保护装置运行状态的设备故障次数,di *为第i个时间点的保护装置运行状态的设备故障次数安全稳定阈值,dj *为数据预处理后第j个样本的设备故障次数,其为0-1之间的一个数值;
步骤2中所述对家族缺陷率进行数据预处理为:
Figure BDA0001924219460000103
其中,ei为第i个时间点的保护装置运行状态的家族缺陷率,ei *为第i个时间点的保护装置运行状态的家族缺陷率安全稳定阈值,ej *为数据预处理后第j个样本的家族缺陷率,其为0-1之间的一个数值;
步骤2中所述对异常告警率进行数据预处理为:
Figure BDA0001924219460000104
其中,fi为第i个时间点的保护装置运行状态的异常告警率,fi *为第i个时间点的保护装置运行状态的异常告警率安全稳定阈值,fj *为数据预处理后第j个样本的异常告警率,其为0-1之间的一个数值;
步骤2中所述对不正确动作次数进行数据预处理为:
Figure BDA0001924219460000105
其中,gi为第i个时间点的保护装置运行状态的不正确动作次数,gi *为第i个时间点的保护装置运行状态的不正确动作次数安全稳定阈值,gj *为数据预处理后第j个样本的不正确动作次数,其为0-1之间的一个数值;
步骤2中所述训练样本数据为:
数据处理后第j个样本的工作电压aj *,数据处理后第j个样本的CPU温度bj *,数据处理后j个样本的绝缘性能cj *,数据处理后第j个样本的设备故障次数dj *,数据处理后第j个样本的家族缺陷率ej *,数据处理后第j个样本的异常告警率fj *,数据处理后第j个样本的断路器不正确动作次数gj *
步骤3:将训练样本数据通过人工标记方法得到标记后数据预处理样本,将标记后数据预处理样本作为支持向量机的输入,支持向量机的核函数为径向基核函数,使用遗传算法对支持向量机的核函数参数、错误惩罚因子进行参数寻优,得到分类效果最好的参数取值以构建优化后支持向量机;
步骤3中所述将训练样本数据通过人工标记方法得到标记后数据预处理样本为:
(xj,yj),j∈[1,N]
xj=(aj *,bj *,cj *,dj *,ej *,fj *,gj *)T
yj∈{-1,1}
其中,N为标记后数据预处理样本的数量,(xj,yj)为第j组标记后数据预处理样本点,xj为第j组标记后数据预处理样本的输入向量;
若通过人工标记判断xj在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,则yj=1,若通过人工标记判断xj在在保护装置运行状态出现故障的情况下获得的样本数据,则yj=-1;
yj∈{-1,1}为第j组标记后数据预处理样本的输出结果,yj=1代表第j个样本的输入向量xj为状态良好的样本数据,状态良好的样本数据是在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,yj=-1代表第j个样本的输入向量xj为状态失效的样本数据,状态失效的样本数据是在保护装置运行状态出现故障的情况下获得的样本数据;
aj *为数据处理后第j个样本的工作电压,bj *为数据处理后第j个样本的CPU温度,cj *为数据处理后第j个样本的绝缘性能,dj *为数据处理后第j个样本的设备故障次数,ej *为数据处理后第j个样本的家族缺陷率,fj *为数据处理后第j个样本的异常告警率,gj *为数据处理后第j个样本的断路器不正确动作次数;
步骤3中所述支持向量机的核函数为径向基核函数,径向基核函数模型为:
Figure BDA0001924219460000121
其中,x为支持向量机的输入向量,x∈{x1,x2,...,xN},xi *为yi=1或yi=-1任意的一个支持向量,支持向量定义为正好落在良好状态分类分类边界上的向量或者失效状态分类分类边界上的向量;
步骤3中所述参数寻优具体为:
遗传算法中的输入数据为支持向量机的核函数参数C、错误惩罚因子δ、以及经过支持向量机的预测分类结果yp *,遗传算法中的优化目标定义为测试样本的平均相对百分比误差;
支持向量机的输入为第j个样本标记后数据预处理样本的输入向量xj,支持向量机的输出为第j个样本标记后数据预处理样本的输出结果yj
支持向量机的核函数类型、核函数参数C、错误惩罚因子δ均会对支持向量机的预测分类结果yp *产生影响;
遗传算法的优化目标模型:
Figure BDA0001924219460000122
其中,Q表示预测数量,Q∈[1,N],yp *表示支持向量机的预测分类输出结果,yp(yp∈{-1,1})表示为第j组标记后数据预处理样本的输出结果;
当遗传算法的优化目标模型函数达到最小值时,即认为此时的分类效果最好,取此时的支持向量机参数为最优参数;
通过核函数模型、最优核函数参数C*、最优错误惩罚因子δ*构建遗传算法优化后支持向量机;
步骤4:将训练样本数据通过人工标记方法得到标记后数据预处理样本,通过遗传算法优化后支持向量机对标记后数据预处理样本进行分类训练,得到失效状态样本分类边界、最优状态分类面、良好状态样本分类边界;
步骤4中所述标记后数据预处理样本为:
(xj,yj),j∈[1,N]
xj=(aj *,bj *,cj *,dj *,ej *,fj *,gj *)T
yj∈{-1,1}
其中,N为标记后数据预处理样本的数量,(xj,yj)为第j组标记后数据预处理样本点,xj为第j组标记后数据预处理样本的输入向量;
若通过人工标记判断xj在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,则yj=1,若通过人工标记判断xj在在保护装置运行状态出现故障的情况下获得的样本数据,则yj=-1;
yj∈{-1,1}为第j组标记后数据预处理样本的输出结果,yj=1代表第j个样本的输入向量xj为状态良好的样本数据,状态良好的样本数据是在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,yj=-1代表第j个样本的输入向量xj为状态失效的样本数据,状态失效的样本数据是在保护装置运行状态出现故障的情况下获得的样本数据;
aj *为数据处理后第j个样本的工作电压,bj *为数据处理后第j个样本的CPU温度,cj *为数据处理后第j个样本的绝缘性能,dj *为数据处理后第j个样本的设备故障次数,ej *为数据处理后第j个样本的家族缺陷率,fj *为数据处理后第j个样本的异常告警率,gj *为数据处理后第j个样本的断路器不正确动作次数;
步骤4中所述对标记后数据预处理样本进行分类训练为:
标记后数据预处理样本(xj,yj),j∈[1,N]可以被分类超平面wT·xj+b=0分开,其中,w=(w1,w2,w3,...,wN)T为分类超平面的法向量,b为代表分类超平面到原点之间的距离;
(xj,yj),j∈[1,N]中与分类超平面最近的距离之和称为分类间隔,分类间隔等于2/||w||,分类间隔最大时的超平面为最优分类超平面;
要使分类间隔最大即是使||w||2最小,因此可以转化为如下的约束优化问题:
Figure BDA0001924219460000131
s.t.yj(wT·xj+b)≥1
其中,w=(w1,w2,w3,...,wN)T为分类超平面的法向量,b为代表分类超平面到原点之间的距离,xj为第j组数据预处理样本的输入向量,yj为第j组数据预处理样本的输出结果;
这个约束优化问题可以通过构造拉格朗日函数求解,求解拉格朗日函数的鞍点,引入拉格朗日因子λj≥0,构造拉格朗日函数如下:
Figure BDA0001924219460000132
其中,λj≥0,j∈[1,N],为拉格朗日因子。
依据拉格朗日对偶理论将
Figure BDA0001924219460000133
转化为对偶问题,即:
Figure BDA0001924219460000134
Figure BDA0001924219460000135
应用二次规划方法可求解,求解得到的最优解α*=[λ1 *2 *,...,λN *]T,则可以得到最优的w*,b*
Figure BDA0001924219460000141
Figure BDA0001924219460000142
其中,xr、xs为yr=1,ys=1或yr=-1,ys=-1任意的一对支持向量,支持向量定义为正好落在良好状态分类分类边界上的向量或者失效状态分类分类边界上的向量;
可知w*T·xj+b*=0为最优分类超平面,w*T·xj+b*=+1为良好状态样本分类边界,w*T·xj+b*=-1为失效状态样本分类边界。
步骤5:将测试样本数据根据步骤1-4得到测试样本数据的输入向量,计算测试样本数据点到最优分类超平面的距离,根据测试样本数据点到最优分类超平面的距离判断测试样本数据的状态;
步骤5中所述测试样本数据到最优分类面的距离为
Figure BDA0001924219460000143
其中,d为测试样本数据到最优分类边界面的距离,x为步骤5中所述测试样本数据的输入向量,其中w*=(w*1,w*2,w*3,...,w*7)T为平面的法向量,b*为一个实数,代表平面到原点之间的距离,w*、b*为步骤四中求得的最优值:
Figure BDA0001924219460000144
Figure BDA0001924219460000145
其中,xr、xs为两个类别中任意的一对支持向量,xi为支持向量机的输入向量,yi为支持向量机的输出结果,λi *为拉格朗日因子;
步骤5中所述根据测试样本数据点到最优分类超平面的距离判断测试样本数据的状态为:
若d>1则测试样本数据属于良好状态;
若0≤d<1则测试样本数据属于注意状态,测试样本数据到良好状态样本分类边界的距离为d'=1-d;
若-1≤d≤0则测试样本数据属于异常状态,测试样本数据到良好状态分类边界面的距离为d'=1-d;
若d<-1则测试样本数据属于失效状态。
步骤6:利用云模型产生云滴的随机性和稳定的趋向性来模拟不同专家对测试样本数据点到良好状态样本分类边界的距离的不同评估值,实现设备的评估值到评语域的不确定转化;
步骤6中所述云模型为:构建云模型,即期望Er、熵En、超熵Hq
其中,期望Er表示云模型的重心位置,反映了定性概念Q的中心值,熵En表示模糊性与随机性相关联程度的大小,超熵Hq为熵的熵,间接反映了云模型的离散程度;
步骤6中所述模拟不同专家对测试样本数据点到良好状态样本分类边界的距离的不同评估值为:
定义健康度H为:保护装置健康状态的度量,其值越大表示健康状态越好;健康度H用d'对应的云模型的期望值表示,云模型的数学期望表达式为:
H(d')=exp[-(d'-Er)2/2En 2]
其中,d'为测试样本数据点到良好状态分类边界面的距离,Er为云模型期望值,En为云模型熵;
云模型的重心Er=0,En=2/3,超熵Hq依据经验取值Hq=0.1;
测试样本数据点到最优分类超平面的距离为:
Figure BDA0001924219460000151
若d>1,测试样本数据点属于良好状态,可以不经过云模型直接对设备状态进行判断为良好状态;
若d<-1,测试样本数据点属于失效状态,可以不经过云模型直接对设备状态进行判断为失效状态;
若d∈[-1,1],需经过云模型的模糊映射:
将d'作为云模型的期望值,随机产生K个随机数,每个随机数有个相对应的云模型隶属度Vi,Vi∈[0,1],i=1,2,...,K;
通过比较Vi与H(d')的大小,H(d')=exp[-(d'-Er)2/2En 2],得到测试样本处于注意状态的数量和异常状态的数量;
当0<Vi<H(d'),i=1,2,...,K时,统计此时Vi的个数计数为N,N<K,此时Vi的状态为异常状态;
当H(d')<Vi<1,i=1,2,...,K时,统计此时Vi的个数计数为K-N,此时Vi的状态为注意状态;
分别统计健康状态处于注意状态的数量K-N和异常状态的数量N,取数量最大的状态为最终的状态,并给出可信度R如下:
Figure BDA0001924219460000152
由支持向量机分类结果可知d'=1为良好状态和注意状态的分类边界,d'为测试样本数据点到良好状态分类边界面的距离;
步骤6中所述实现设备的评估值到评语域的不确定转化:
由云模型的数学期望表达式确定其值为H(d'),从而确定由健康度H向评语域转化;
将得到的健康度转化为百分制:规定良好状态的样本健康度为100,失效状态的样本健康度为0,
注意状态的样本转化公式如下:
Figure BDA0001924219460000153
异常状态的样本转化公式如下:
Figure BDA0001924219460000154
其中,H为云模型的数学期望值,H'为健康状态转化为百分制后的得分。得到继电保护装置最终状态评估结果;
测试样本数据点到最优分类超平面的距离为:
Figure BDA0001924219460000155
若d>1,测试样本数据点属于良好状态,可以不经过云模型直接对设备状态进行判断为良好状态;
若d<-1,测试样本数据点属于失效状态,可以不经过云模型直接对设备状态进行判断为失效状态;
若d∈[-1,1],需经过云模型的模糊映射:
将d'作为云模型的期望值,随机产生K个随机数,每个随机数有个相对应的云模型隶属度Vi,Vi∈[0,1],i=1,2,...,K;
通过比较Vi与H(d')的大小,H(d')=exp[-(d'-Er)2/2En 2],得到测试样本处于注意状态的数量和异常状态的数量;
当0<Vi<H(d'),i=1,2,...,K时,统计此时Vi的个数计数为N,N<K,此时Vi的状态为异常状态;
当H(d')<Vi<1,i=1,2,...,K时,统计此时Vi的个数计数为K-N,此时Vi的状态为注意状态;
分别统计健康状态处于注意状态的数量K-N和异常状态的数量N,取数量最大的状态为最终的状态。
为了检验本方法的可靠性和分类效率,在保证实验数据不变的情况下,利用神经网络和标准支持向量机方法对继电保护装置状态进行评估,并且和遗传算法优化支持向量机参数状态评估模型进行比较。本文所提方法能有效评估继电保护装置的运行状态,使得运维检修人员能够及时掌握继电保护装置的运行状态,为检修人员合理安排检修计划提供参考意见,防止设备出现安全事故,保证供电的安全可靠。
本文中所描述的具体实施案例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施案例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (6)

1.一种云模型与遗传算法优化支持向量机的状态评估方法,其特征在于,包括以下步骤:
步骤1:根据继电保护装置运行状态信息,选取工作电压、CPU温度、绝缘性能、设备故障次数、家族缺陷率、异常告警率、断路器不正确动作次数作为支持向量机的输入特征向量;
步骤2:分别对工作电压、CPU温度、绝缘性能、设备故障次数、家族缺陷率、异常告警率、断路器不正确动作次数进行数据预处理,得到训练样本数据;
步骤3:将训练样本数据通过人工标记方法得到标记后数据预处理样本,将标记后数据预处理样本作为支持向量机的输入,支持向量机的核函数为径向基核函数,使用遗传算法对支持向量机的核函数参数、错误惩罚因子进行参数寻优,得到分类效果最好的参数取值以构建优化后支持向量机;
步骤4:将训练样本数据通过人工标记方法得到标记后数据预处理样本,通过遗传算法优化后支持向量机对标记后数据预处理样本进行分类训练,得到失效状态样本分类边界、最优状态分类面、良好状态样本分类边界;
步骤5:将测试样本数据根据步骤1-4得到测试样本数据的输入向量,计算测试样本数据点到最优分类超平面的距离,根据测试样本数据点到最优分类超平面的距离判断测试样本数据的状态;
步骤6:利用云模型产生云滴的随机性和稳定的趋向性来模拟不同专家对测试样本数据点到良好状态样本分类边界的距离的不同评估值,实现设备的评估值到评语域的不确定转化;
步骤6中所述云模型为:构建云模型,即期望Er、熵En、超熵Hq
其中,期望Er表示云模型的重心位置,反映了定性概念Q的中心值,熵En表示模糊性与随机性相关联程度的大小,超熵Hq为熵的熵,间接反映了云模型的离散程度;
步骤6中所述模拟不同专家对测试样本数据点到良好状态样本分类边界的距离的不同评估值为:
定义健康度H为:保护装置健康状态的度量,其值越大表示健康状态越好;健康度H用d'对应的云模型的期望值表示,云模型的数学期望表达式为:
H(d')=exp[-(d'-Er)2/2En 2]
其中,d'为测试样本数据点到良好状态分类边界面的距离,Er为云模型期望值,En为云模型熵;
云模型的重心Er=0,En=2/3,超熵Hq依据经验取值Hq=0.1;
测试样本数据点到最优分类超平面的距离为:
Figure FDA0003810295250000011
若d>1,测试样本数据点属于良好状态,可以不经过云模型直接对设备状态进行判断为良好状态;
若d<-1,测试样本数据点属于失效状态,可以不经过云模型直接对设备状态进行判断为失效状态;
若d∈[-1,1],需经过云模型的模糊映射:
将d'作为云模型的期望值,随机产生K个随机数,每个随机数有个相对应的云模型隶属度Vi,Vi∈[0,1],i=1,2,...,K;
通过比较Vi与H(d')的大小,H(d')=exp[-(d'-Er)2/2En 2],得到测试样本处于注意状态的数量和异常状态的数量;
当0<Vi<H(d'),i=1,2,...,K时,统计此时Vi的个数计数为N,N<K,此时Vi的状态为异常状态;
当H(d')<Vi<1,i=1,2,...,K时,统计此时Vi的个数计数为K-N,此时Vi的状态为注意状态;
分别统计健康状态处于注意状态的数量K-N和异常状态的数量N,取数量最大的状态为最终的状态,并给出可信度R如下:
Figure FDA0003810295250000021
由支持向量机分类结果可知d'=1为良好状态和注意状态的分类边界,d'为测试样本数据点到良好状态分类边界面的距离;
步骤6中所述实现设备的评估值到评语域的不确定转化:
由云模型的数学期望表达式确定其值为H(d'),从而确定由健康度H向评语域转化;
将得到的健康度转化为百分制:规定良好状态的样本健康度为100,失效状态的样本健康度为0,
注意状态的样本转化公式如下:
Figure FDA0003810295250000022
异常状态的样本转化公式如下:
Figure FDA0003810295250000023
其中,H为云模型的数学期望值,H'为健康状态转化为百分制后的得分,得到继电保护装置最终状态评估结果。
2.根据权利要求1所述的云模型与遗传算法优化支持向量机的状态评估方法,其特征在于:
步骤1中所述工作电压ai为第i个时间点的继电保护装置运行状态的工作电压;
步骤1中所述CPU温度bi为第i个时间点的继电保护装置运行状态的CPU温度;
步骤1中所述绝缘性能ci为第i个时间点的继电保护装置运行状态的绝缘性能;
步骤1中所述设备故障次数di为第i个时间点的继电保护装置运行状态的设备故障次数;
步骤1中所述家族缺陷率ei为第i个时间点的继电保护装置运行状态的家族缺陷率;
步骤1中所述异常告警率fi为第i个时间点的继电保护装置运行状态的异常告警率;
步骤1中所述断路器不正确动作次数gi为第i个时间点的继电保护装置运行状态的断路器不正确动作次数;
i∈[0,M],M为继电保护装置运行时间。
3.根据权利要求1所述的云模型与遗传算法优化支持向量机的状态评估方法,其特征在于:
步骤2中所述对工作电压进行数据预处理为:
Figure FDA0003810295250000031
其中,M为继电保护装置运行时间,N为数据预处理后的样本数量,ai为第i个时间点的继电保护装置运行状态的工作电压,ai *为第i个时间点的保护装置运行状态的工作电压安全稳定阈值,aj *为数据预处理后第j个样本的工作电压,其为0-1之间的一个数值;
步骤2中所述对CPU温度进行数据预处理为:
Figure FDA0003810295250000032
其中,bi为第i个时间点的保护装置运行状态的CPU温度,bi *为第i个时间点的保护装置运行状态的CPU温度安全稳定阈值,bj *为数据预处理后第j个样本的CPU温度,其为0-1之间的一个数值;
步骤2中所述对绝缘性能进行数据预处理为:
Figure FDA0003810295250000033
其中,ci为第i个时间点的保护装置运行状态的绝缘性能,ci *为第i个时间点的保护装置运行状态的绝缘性能安全稳定阈值,cj *为数据预处理后第j个样本的绝缘性能,其为0-1之间的一个数值;
步骤2中所述对设备故障次数进行数据预处理为:
Figure FDA0003810295250000034
其中,di为第i个时间点的保护装置运行状态的设备故障次数,di *为第i个时间点的保护装置运行状态的设备故障次数安全稳定阈值,dj *为数据预处理后第j个样本的设备故障次数,其为0-1之间的一个数值;
步骤2中所述对家族缺陷率进行数据预处理为:
Figure FDA0003810295250000035
其中,ei为第i个时间点的保护装置运行状态的家族缺陷率,ei *为第i个时间点的保护装置运行状态的家族缺陷率安全稳定阈值,ej *为数据预处理后第j个样本的家族缺陷率,其为0-1之间的一个数值;
步骤2中所述对异常告警率进行数据预处理为:
Figure FDA0003810295250000041
其中,fi为第i个时间点的保护装置运行状态的异常告警率,fi *为第i个时间点的保护装置运行状态的异常告警率安全稳定阈值,fj *为数据预处理后第j个样本的异常告警率,其为0-1之间的一个数值;
步骤2中所述对不正确动作次数进行数据预处理为:
Figure FDA0003810295250000042
其中,gi为第i个时间点的保护装置运行状态的不正确动作次数,gi *为第i个时间点的保护装置运行状态的不正确动作次数安全稳定阈值,gj *为数据预处理后第j个样本的不正确动作次数,其为0-1之间的一个数值;
步骤2中所述训练样本数据为:
数据处理后第j个样本的工作电压aj *,数据处理后第j个样本的CPU温度bj *,数据处理后j个样本的绝缘性能cj *,数据处理后第j个样本的设备故障次数dj *,数据处理后第j个样本的家族缺陷率ej *,数据处理后第j个样本的异常告警率fj *,数据处理后第j个样本的断路器不正确动作次数gj *
4.根据权利要求1所述的云模型与遗传算法优化支持向量机的状态评估方法,其特征在于:
步骤3中所述将训练样本数据通过人工标记方法得到标记后数据预处理样本为:
(xj,yj),j∈[1,N]
xj=(aj *,bj *,cj *,dj *,ej *,fj *,gj *)T
yj∈{-1,1}
其中,N为标记后数据预处理样本的数量,(xj,yj)为第j组标记后数据预处理样本点,xj为第j组标记后数据预处理样本的输入向量;
若通过人工标记判断xj在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,则yj=1,若通过人工标记判断xj在保护装置运行状态出现故障的情况下获得的样本数据,则yj=-1;
yj∈{-1,1}为第j组标记后数据预处理样本的输出结果,yj=1代表第j个样本的输入向量xj为状态良好的样本数据,状态良好的样本数据是在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,yj=-1代表第j个样本的输入向量xj为状态失效的样本数据,状态失效的样本数据是在保护装置运行状态出现故障的情况下获得的样本数据;
aj*为数据处理后第j个样本的工作电压,bj *为数据处理后第j个样本的CPU温度,cj *为数据处理后第j个样本的绝缘性能,dj *为数据处理后第j个样本的设备故障次数,ej *为数据处理后第j个样本的家族缺陷率,fj *为数据处理后第j个样本的异常告警率,gj *为数据处理后第j个样本的断路器不正确动作次数;
步骤3中所述支持向量机的核函数为径向基核函数,径向基核函数模型为:
Figure FDA0003810295250000051
其中,x为支持向量机的输入向量,x∈{x1,x2,...,xN},xi *为yi=1或yi=-1任意的一个支持向量,支持向量定义为正好落在良好状态分类边界上的向量或者失效状态分类边界上的向量;
步骤3中所述参数寻优具体为:
遗传算法中的输入数据为支持向量机的核函数参数C、错误惩罚因子δ、以及经过支持向量机的预测分类结果yp *,遗传算法中的优化目标定义为测试样本的平均相对百分比误差;
支持向量机的输入为第j个样本标记后数据预处理样本的输入向量xj,支持向量机的输出为第j个样本标记后数据预处理样本的输出结果yj
支持向量机的核函数类型、核函数参数C、错误惩罚因子δ均会对支持向量机的预测分类结果yp *产生影响;
遗传算法的优化目标模型:
Figure FDA0003810295250000052
其中,Q表示预测数量,Q∈[1,N],yp *表示支持向量机的预测分类输出结果,yp(yp∈{-1,1})表示为第j组标记后数据预处理样本的输出结果;
当遗传算法的优化目标模型函数达到最小值时,即认为此时的分类效果最好,取此时的支持向量机参数为最优参数;
通过核函数模型、最优核函数参数C*、最优错误惩罚因子δ*构建遗传算法优化后支持向量机。
5.根据权利要求1所述的云模型与遗传算法优化支持向量机的状态评估方法,其特征在于:
步骤4中所述标记后数据预处理样本为:
(xj,yj),j∈[1,N]
xj=(aj *,bj *,cj *,dj *,ej *,fj *,gj *)T
yj∈{-1,1}
其中,N为标记后数据预处理样本的数量,(xj,yj)为第j组标记后数据预处理样本点,xj为第j组标记后数据预处理样本的输入向量;
若通过人工标记判断xj在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,则yj=1,若通过人工标记判断xj在保护装置运行状态出现故障的情况下获得的样本数据,则yj=-1;
yj∈{-1,1}为第j组标记后数据预处理样本的输出结果,yj=1代表第j个样本的输入向量xj为状态良好的样本数据,状态良好的样本数据是在保护装置刚刚投入运行,且运行状态良好的情况下获得的样本数据,yj=-1代表第j个样本的输入向量xj为状态失效的样本数据,状态失效的样本数据是在保护装置运行状态出现故障的情况下获得的样本数据;
aj *为数据处理后第j个样本的工作电压,bj *为数据处理后第j个样本的CPU温度,cj *为数据处理后第j个样本的绝缘性能,dj *为数据处理后第j个样本的设备故障次数,ej *为数据处理后第j个样本的家族缺陷率,fj *为数据处理后第j个样本的异常告警率,gj *为数据处理后第j个样本的断路器不正确动作次数;
步骤4中所述对标记后数据预处理样本进行分类训练为:
标记后数据预处理样本(xj,yj),j∈[1,N]可以被分类超平面wT·xj+b=0分开,其中,w=(w1,w2,w3,...,wN)T为分类超平面的法向量,b为代表分类超平面到原点之间的距离;
(xj,yj),j∈[1,N]中与分类超平面最近的距离之和称为分类间隔,分类间隔等于2/||w||,分类间隔最大时的超平面为最优分类超平面;
要使分类间隔最大即是使||w||2最小,因此可以转化为如下的约束优化问题:
Figure FDA0003810295250000061
s.t.yj(wT·xj+b)≥1
其中,w=(w1,w2,w3,...,wN)T为分类超平面的法向量,b为代表分类超平面到原点之间的距离,xj为第j组数据预处理样本的输入向量,yj为第j组数据预处理样本的输出结果;
这个约束优化问题可以通过构造拉格朗日函数求解,求解拉格朗日函数的鞍点,引入拉格朗日因子λj≥0,构造拉格朗日函数如下:
Figure FDA0003810295250000071
其中,λj≥0,j∈[1,N],为拉格朗日因子;
依据拉格朗日对偶理论将
Figure FDA0003810295250000072
转化为对偶问题,即:
Figure FDA0003810295250000073
Figure FDA0003810295250000074
应用二次规划方法可求解,求解得到的最优解α*=[λ1 *2 *,...,λN *]T,则可以得到最优的w*,b*
Figure FDA0003810295250000075
Figure FDA0003810295250000076
其中,xr、xs为yr=1,ys=1或yr=-1,ys=-1任意的一对支持向量,支持向量定义为正好落在良好状态分类边界上的向量或者失效状态分类边界上的向量;
可知w*T·xj+b*=0为最优分类超平面,w*T·xj+b*=+1为良好状态样本分类边界,w*T·xj+b*=-1为失效状态样本分类边界。
6.根据权利要求1所述的云模型与遗传算法优化支持向量机的状态评估方法,其特征在于:
步骤5中所述测试样本数据到最优分类面的距离为
Figure FDA0003810295250000077
其中,d为测试样本数据到最优分类边界面的距离,x为步骤5中所述测试样本数据的输入向量,其中w*=(w*1,w*2,w*3,...,w*7)T为平面的法向量,b*为一个实数,代表平面到原点之间的距离,w*、b*为步骤四中求得的最优值:
Figure FDA0003810295250000078
Figure FDA0003810295250000079
其中,xr、xs为两个类别中任意的一对支持向量,xi为支持向量机的输入向量,yi为支持向量机的输出结果,λi *为拉格朗日因子;
步骤5中所述根据测试样本数据点到最优分类超平面的距离判断测试样本数据的状态为:
若d>1则测试样本数据属于良好状态;
若0≤d<1则测试样本数据属于注意状态,测试样本数据到良好状态样本分类边界的距离为d'=1-d;
若-1≤d≤0则测试样本数据属于异常状态,测试样本数据到良好状态分类边界面的距离为d'=1-d;
若d<-1则测试样本数据属于失效状态。
CN201811608793.XA 2018-12-27 2018-12-27 一种云模型与遗传算法优化支持向量机的状态评估方法 Active CN109583520B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811608793.XA CN109583520B (zh) 2018-12-27 2018-12-27 一种云模型与遗传算法优化支持向量机的状态评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811608793.XA CN109583520B (zh) 2018-12-27 2018-12-27 一种云模型与遗传算法优化支持向量机的状态评估方法

Publications (2)

Publication Number Publication Date
CN109583520A CN109583520A (zh) 2019-04-05
CN109583520B true CN109583520B (zh) 2023-04-07

Family

ID=65932037

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811608793.XA Active CN109583520B (zh) 2018-12-27 2018-12-27 一种云模型与遗传算法优化支持向量机的状态评估方法

Country Status (1)

Country Link
CN (1) CN109583520B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110032981B (zh) * 2019-04-19 2022-07-26 电子科技大学 基于改进支持向量机的旋转机械故障识别方法
CN110082136B (zh) * 2019-04-19 2021-03-30 电子科技大学 基于云遗传算法优化支持向量机的旋转机械故障诊断方法
CN110287543B (zh) * 2019-05-30 2022-05-13 武汉大学 一种继电保护装置寿命预测方法
CN110880062B (zh) * 2019-10-31 2022-07-08 广西电网有限责任公司电力科学研究院 一种配电设备状态检修时间的确定方法
CN111239516A (zh) * 2020-01-19 2020-06-05 广东电网有限责任公司计量中心 一种互感器寿命预测方法及装置
CN113780338B (zh) * 2021-07-30 2024-04-09 国家计算机网络与信息安全管理中心 基于支持向量机的大数据分析中置信度评价方法、系统、设备及存储介质
CN114298208A (zh) * 2021-12-24 2022-04-08 东南大学 一种用于变压器故障诊断的支持向量机参数优化算法
CN115994305B (zh) * 2022-12-06 2024-04-26 中国电力科学研究院有限公司 一种确定继电保护设备运行状态的方法及系统
CN116451876B (zh) * 2023-06-15 2023-09-22 国网江西省电力有限公司信息通信分公司 一种基于人工智能的配电网故障预测和主动检修系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011085435A1 (en) * 2010-01-14 2011-07-21 The University Of Sydney Classification process for an extracted object or terrain feature
CN105866725A (zh) * 2016-04-20 2016-08-17 国网上海市电力公司 一种基于聚类分析和云模型的智能电表故障分类方法
CN106053067A (zh) * 2016-05-24 2016-10-26 广东石油化工学院 基于量子遗传算法优化支持向量机的轴承故障诊断方法
CN106096771A (zh) * 2016-06-07 2016-11-09 广州供电局有限公司 支持向量机的参数选取方法和系统与负荷预测方法和系统
CN108051660A (zh) * 2017-10-31 2018-05-18 华北电力大学(保定) 一种变压器故障组合诊断模型建立方法及诊断方法
CN108171432A (zh) * 2018-01-04 2018-06-15 南京大学 基于多维云模型-模糊支持向量机的生态风险评价方法
CN108680814A (zh) * 2018-08-29 2018-10-19 国网河北省电力有限公司电力科学研究院 一种多维度变压器运行状态评价方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011085435A1 (en) * 2010-01-14 2011-07-21 The University Of Sydney Classification process for an extracted object or terrain feature
CN105866725A (zh) * 2016-04-20 2016-08-17 国网上海市电力公司 一种基于聚类分析和云模型的智能电表故障分类方法
CN106053067A (zh) * 2016-05-24 2016-10-26 广东石油化工学院 基于量子遗传算法优化支持向量机的轴承故障诊断方法
CN106096771A (zh) * 2016-06-07 2016-11-09 广州供电局有限公司 支持向量机的参数选取方法和系统与负荷预测方法和系统
CN108051660A (zh) * 2017-10-31 2018-05-18 华北电力大学(保定) 一种变压器故障组合诊断模型建立方法及诊断方法
CN108171432A (zh) * 2018-01-04 2018-06-15 南京大学 基于多维云模型-模糊支持向量机的生态风险评价方法
CN108680814A (zh) * 2018-08-29 2018-10-19 国网河北省电力有限公司电力科学研究院 一种多维度变压器运行状态评价方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
GA-SVM modeling of multiclass seizure detector in epilepsy analysis system using cloud computing;Chia-Ping Shen等;《Soft Computing》;20151114;第21卷;2139–2149 *
Research and Application of an Air Quality Early Warning System Based on a Modified Least Squares Support Vector Machine and a Cloud Model;Jianzhou Wang等;《Public Health》;20170302;第14卷(第03期);1-11 *
基于云理论的智能变电站二次设备状态评估;王月月等;《电力系统保护与控制》;20180109;第46卷(第01期);71-77 *
基于支持向量机和云模型的网络健康状态评估;温祥西等;《北京邮电大学学报》;20111123;第35卷(第01期);10-14 *

Also Published As

Publication number Publication date
CN109583520A (zh) 2019-04-05

Similar Documents

Publication Publication Date Title
CN109583520B (zh) 一种云模型与遗传算法优化支持向量机的状态评估方法
CN105117602B (zh) 一种计量装置运行状态预警方法
CN108320043B (zh) 一种基于电力大数据的配电网设备状态诊断预测方法
CN113156917B (zh) 基于人工智能的电网设备故障诊断方法及系统
CN107358366B (zh) 一种配电变压器故障风险监测方法及系统
CN103793854B (zh) 多重组合优化的架空输电线路运行风险信息化评估方法
CN102496069B (zh) 基于模糊层次分析法的电缆多状态安全运行评估方法
CN106651169A (zh) 基于模糊综合评价的配电自动化终端状态评价方法及系统
CN104914327A (zh) 基于实时监测信息的变压器故障检修预测方法
CN105046591A (zh) 一种电力用户用电能效评估方法
CN111598352A (zh) 一种基于贝叶斯网络的混凝土梁式桥综合评估方法
CN113887846B (zh) 一种电容式电压互感器的超差风险预警方法
CN110287543B (zh) 一种继电保护装置寿命预测方法
CN109492790A (zh) 基于神经网络与数据挖掘的风电机组健康管理方法
CN108764659A (zh) 一种基于模糊综合法的配电网可靠性影响因素分析方法
CN112257914B (zh) 一种基于随机森林的航空安全因果预测方法
CN117560300B (zh) 一种智能物联网流量预测与优化系统
CN106548284A (zh) 一种面向运行调控的自适应模块化电网安全预警评价方法
CN103675518A (zh) 油浸式高压设备检修方法及装置
CN112434887A (zh) 一种结合网络核密度估计和svm的供水管网风险预测方法
CN105741184A (zh) 一种变压器状态评估方法及装置
CN115864644A (zh) 一种继电保护装置状态评价方法、系统、设备及介质
CN115936663A (zh) 一种电力系统的维护方法及装置
CN112507290B (zh) 配电设备故障概率预判方法、装置及存储介质
CN115564170A (zh) 一种配电智能终端可靠性评估方法

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Zhang Xuemin

Inventor after: Liu Qi

Inventor after: Han Zongyan

Inventor after: Zhang Zhiqiang

Inventor after: Zhou Hongsheng

Inventor after: Dai Weikang

Inventor after: Bai Jianlin

Inventor after: Du Linqiang

Inventor after: Yi Lin Fei

Inventor after: Xu Xin

Inventor after: Wang Bin

Inventor after: Shi Yingchun

Inventor after: Wang Wenlin

Inventor after: Dang Junpeng

Inventor after: Yang Yongxu

Inventor after: Ma Zhipeng

Inventor after: Tang Yiheng

Inventor after: Li Lei

Inventor before: Zhang Xuemin

Inventor before: Li Lei

Inventor before: Liu Qi

Inventor before: Han Zongyan

Inventor before: Zhang Zhiqiang

Inventor before: Zhou Hongsheng

Inventor before: Dai Weikang

Inventor before: Bai Jianlin

Inventor before: Du Linqiang

Inventor before: Yi Lin Fei

Inventor before: Xu Xin

Inventor before: Wang Bin

Inventor before: Shi Yingchun

Inventor before: Wang Wenlin

Inventor before: Dang Junpeng

Inventor before: Yang Yongxu

Inventor before: Ma Zhipeng

Inventor before: Tang Yiheng

Inventor before: Chen Haitao

GR01 Patent grant
GR01 Patent grant