CN108710756A - 基于谱聚类分析下多特征信息加权融合的故障诊断方法 - Google Patents
基于谱聚类分析下多特征信息加权融合的故障诊断方法 Download PDFInfo
- Publication number
- CN108710756A CN108710756A CN201810479820.1A CN201810479820A CN108710756A CN 108710756 A CN108710756 A CN 108710756A CN 201810479820 A CN201810479820 A CN 201810479820A CN 108710756 A CN108710756 A CN 108710756A
- Authority
- CN
- China
- Prior art keywords
- svm
- class
- evidence
- spectral clustering
- faulty equipment
- 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
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
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/257—Belief theory, e.g. Dempster-Shafer
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Test And Diagnosis Of Digital Computers (AREA)
Abstract
本发明涉及一种基于谱聚类分析下多特征信息加权融合的故障诊断方法,该方法首先对故障设备进行谱聚类分析,其次通过获取各SVM局部诊断证据对各故障模式的可靠度,同时,由各SVM局部诊断硬输出判决矩阵构造出基本概率分配,并对基本概率分配进行加权处理,获取可信度和不确定度,最后通过设定的诊断规则,结合可信度和不确定度进行诊断。与现有技术相比,本发明考虑了不同来源的证据对辨识框架中各命题的识别具有不同的可靠性,降低了各SVM局部诊断间的冲突,实现了SVM和改进证据理论的有效结合,进而解决了识别的不可靠性造成的合成结果不能反映客观事实的缺点。
Description
技术领域
本发明涉及电厂设备故障诊断技术领域,尤其是涉及一种基于谱聚类分析下多特征信息加权融合的故障诊断方法。
背景技术
当电厂设备发生故障时,监控***会在较短时间内将大量的信息提供给运行人员,其中也包含了大量不必要上传的无用信息,给故障的及时处理带来了严重的阻碍。另外,SCADA/EMS所能提供的信息有限,这些信息并不能完全满足运行人员对故障进行全面分析的需求,保护及开关的误动、拒动,以及因通信信道干扰所造成的信息缺失也均会使基于单一信息源的故障分析的准确性受到严重影响。而随着电厂自动化程度的不断提高以及各种智能电子设备的使用,运行人员可以获得更多的设备信息,且电力通信网的快速发展也使电网故障分析时利用多源信息成为可能。
因此,在智能化、规模化的背景下,有效地精简用于故障分析的信息,通过冗余、异构的多信息源数据融合,有效地减小因保护及开关的误动、拒动,以及因信道干扰造成的信息缺失等的影响,将大大有助于运行人员对故障的正确分析,进而保证及时对故障进行处理,从而保障电网的安全、稳定运行,具有现实的社会效益及其明显的经济效益。将SVM与证据理论相结合在多信息融合的故障诊断中进行应用具有较好的应用前景,然而,由于证据合成时将每个证据体都视为同等重要,没有考虑到不同来源的证据对辨识框架中各命题的识别具有不同的可靠性这一事实,这就造成了合成结果不能反映客观事实的缺点。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于谱聚类分析下多特征信息加权融合的故障诊断方法。
本发明的目的可以通过以下技术方案来实现:
基于谱聚类分析下多特征信息加权融合的故障诊断方法,该方法包括以下步骤:
S1:对各个电厂故障设备进行判断,若某两个故障设备情况集中,则将二者置于同一分区,否则,采用谱聚类方法对故障设备进行划分,获取多个分区域;具体包括以下步骤:
11)将多个故障设备划分为同一样本集合:
假设有n个故障设备,则需要被划分的样本集合为:
X={X1,X2,...,Xi,...,Xn}∈Rm×n
式中,Rm×n为实数集,m为矩阵行数,Xi为第i台故障设备中通过对该机组数据Xi(t)进行采样所构成的参数集合,Xi的表达式为:
Xi={xi(1),xi(2),...,xi(t),,...,xi(m)}T
12)获取样本集合的Markov转移矩阵:
将各个故障设备对应于高斯权重图的各个顶点,获取两点间的相似矩阵Aij,其表达式为:
式中,|xi-xj|2为两点之间的欧几里得距离,σ为尺度参数;
则样本集合的Markov转移矩阵Pij为:
13)对Markov转移矩阵进行谱分析,获取故障设备的划分区域数量:
假设某样本集合的Markov转移矩阵为P,对其进行分解:
式中,λz为矩阵P的第z个特征值,为λz的左特征向量,ψz为λz的右特征向量。
若矩阵P有q<n个主要特征值,其中λq<<λq-1,λ1,λ2,...,λq数值相近且趋于数值1,则P的相似矩阵Pq可由前q个特征值和对应特征向量计算;若不符合前述情况,但矩阵P相邻特征值之间的差值骤降,即λq-λq+1>>λk-λk+1,q+1≤k≤n,λq+1值小于0.001,则依旧按照前述情况处理。
对相似矩阵Pq获取谱分解式:
即q为所有故障设备中待划分区域数。
14)对样本集合中各个点进行扩散距离计算,确定划分区域:
定义映射ψ:
其中,ψ为R→R,设任意两故障设备xi与xj之间的扩散距离为:
式中,e1、...、en为x1、...、xn对应的单位向量;
根据λ1,λ2,...,λq及其特征向量,将xi与xj之间的扩散距离简化:
考虑各点间的扩散距离以及分区域数,给出定限值η>0,当D2(xi,xj)≤η时,确定其对应的故障设备划分在同一分区中。
S2:对各个分区域中的故障设备获取相关系数,并进行正序排序,将绝对值大于设定阈值的相关系数所对应的故障设备作为待诊断目标设备,具体内容包括:
21)假设存在某一个含有故障设备数为a的分区域,对其中故障设备分别进行编号:1,2,3,...,a,则第i个故障设备的相关系数ri为:
式中,b为采样数,xik为第i个分区域中,第k个采样点实测参数,yk为包含点k所在分区的实测参数;为第i个分区平均参数,为故障设备b个时刻功率的参数平均值;
22)对步骤21)获取的ri的绝对值按正序排序,将绝对值大且大于设定的阈值的故障设备,设为该分区域中的标准故障设备,具体条件如下:
式中,S为前S个相关性最高故障设备,rall为所有故障设备的相关系数。
S3:结合SVM局部诊断及改进的证据理论方法,对待诊断目标设备的故障模式计算基本概率分配并进行加权组合,获取可信度和不确定度,具体内容包括:
31)对待诊断目标设备的故障模式获取特征向量,结合证据理论与SVM局部诊断,获取各故障模式的证据体的隶属度;
假设辨识框架Θ中共有l类故障模式,B为故障模式,则利用两两分类方法建立多分类SVM的输出判决矩阵SM为:
式中,F(Bg,Bh)为故障模式是属于Bi类还是属于Bh类的标准SVM判决的硬输出,若属于Bi类,则F(Bg,Bh)=+1,若属于Bh类,则F(Bg,Bh)=-1,判决矩阵SM第一行表示第B1类与其他类两两分类的SVM硬输出判决结果,判决矩阵第l行表示为第Bl类与其他类两两分类的SVM硬输出判决结果;
对判决矩阵SM的第h行向量可知,第Bh类与其它类进行两两分类的总分类次数为(l-1)次,同时第h行判决硬输出为1的次数之和,即第Bh类与其它类进行两两分类时判决属于第Bj类的次数,将判决属于第Bh类的分类次数与第Bh类参与的总分类次数之比定义为第Bh类的隶属度uh,其表达式为:
将不能确定属于任何一类的隶属度,定义为qΘ:
uΘ=1-max{u1,u2,…,ul}
则基本概率分配m(Bh)的表达式为:
式中,M为SVM多分类的规模;
32)根据证据体的隶属度获取SVM的基本概率分配,并将各证据体按基本概率分配进行加权处理,获取多个证据加权融合后的最终诊断结果:
对于同一辨识框架Θ,假设SVM证据具体对辨识框架中c个命题的可靠度为Re(B)→[0,1],则辨识框架Θ上的可靠度加权系数V(B)为:
进行多证据加权融合,获取可信度W(B)和不确定度W(Θ):
S4:根据可信度和不确定度对辨识框架中所有故障模式进行诊断:
设定诊断规则,根据诊断规则对诊断结果进行判断,假设诊断结果为Fc,所述的诊断规则包括:
规则一:
规则二:W(Fc)-W(Fj)>ε;W(Fc)>W(Θ)
规则三:W(Θ)<γ
规则一为诊断的基本条件,即所判定故障模式具有最大的可信度;规则二表明所判定故障模式的可信度必须比其它故障模式的可信度大,其中ε为设定的判定阈值;规则三表明不确定度必须小于阈值γ,保证故障样本是充分可判断的。
与现有技术相比,本发明具有以下优点:
(1)本发明方法先对电厂故障设备情况等进行现场考察,若两故障设备情况相对集中,各参数数据相差不多则放在同一个分区;反之,则采用谱聚类分析法划分该区域,在保证其准确性的同时,也提高了划分效率;
(2)本发明方法首先对故障设备进行谱聚类分析,以考虑不同来源的证据对辨识框架中各命题的识别具有不同的可靠性,其次通过获取各SVM局部诊断证据对各故障模式的可靠度,同时由各SVM局部诊断硬输出判决矩阵构造出基本概率分配,并对基本概率分配进行加权处理,降低了各SVM局部诊断间的冲突,实现了SVM和改进证据理论的有效结合,进而解决了识别的不可靠性造成的合成结果不能反映客观事实的缺点;
(3)本发明方法采用设定判断规则进行故障诊断,判定阈值可根据实际情况来选择,本发明方法可允许用户按照自己的需要建立融合多种数据的判断规则,并可自动积累,进而提高设备和部件分析和早期主动诊断的智能化程度。
附图说明
图1为本发明方法的流程图;
图2为高斯权重图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
如图1所示,本发明涉及一种基于谱聚类分析下多特征信息加权融合的故障诊断方法,包括以下步骤:
S1:对各个电厂故障设备进行判断,若某两个故障设备情况集中,则将二者置于同一分区,否则,采用谱聚类方法对故障设备进行划分,获取多个分区域;具体内容包括:
11)将多个故障设备划分为同一样本集合。
假设有n个故障设备,则需要被划分的样本集合为:
X={X1,X2,...,Xi,...,Xn}∈Rm×n
式中,Rm×n为实数集,m为矩阵行数,Xi为第i台故障设备中通过对该机组数据Xi(t)进行采样所构成的参数集合,Xi的表达式为:
Xi={xi(1),xi(2),...,xi(t),,...,xi(m)}T
12)获取样本集合的Markov转移矩阵。
将各个故障设备对应于高斯权重图的各个顶点,如图2所示,获取两点间的相似矩阵Aij,其表达式为:
式中,|xi-xj|2为两点之间的欧几里得距离,σ为尺度参数;
则样本集合的Markov转移矩阵Pij为:
13)对Markov转移矩阵进行谱分析,获取故障设备的划分区域数量。
假设某样本集合的Markov转移矩阵为P,对其进行分解:
式中,λz为矩阵P的第z个特征值,为λz的左特征向量,ψz为λz的右特征向量;
若矩阵P有q<n个主要特征值,其中λq<<λq-1,λ1,λ2,...,λq数值相近且趋于数值1,则P的相似矩阵Pq可由前q个特征值和对应特征向量计算;若不符合该情况,但矩阵P相邻特征值之间的差值骤降,即λq-λq+1>>λk-λk+1,q+1≤k≤n,λq+1值小于0.001,则依旧按照由前q个特征值和对应特征向量计算进行处理;
对相似矩阵Pq获取谱分解式可得:
即q为所有故障设备中待划分区域数。
14)对样本集合中各个点进行扩散距离计算,确定划分区域。
定义映射ψ:
其中,ψ为R→R,设任意两故障设备xi与xj之间的扩散距离为:
式中,e1、...、en为x1、...、xn对应的单位向量;
根据λ1,λ2,...,λq及其特征向量,将xi与xj之间的扩散距离简化:
考虑各点间的扩散距离以及分区域数,给出定限值η>0,当D2(xi,xj)≤η时,确定其对应的故障设备划分在同一分区中。
S2:对各个分区域中的故障设备获取相关系数,并进行正序排序,选取绝对值大于设定阈值的相关系数所对应的故障设备作为待诊断目标设备;
S3:结合SVM局部诊断及改进的证据理论方法,对待诊断目标设备的故障模式计算基本概率分配并进行加权组合,获取可信度和不确定度;具体包括以下步骤:
31)对待诊断目标设备的故障模式获取特征向量,结合证据理论与SVM局部诊断,获取各故障模式的证据体的隶属度;
假设辨识框架Θ中共有l类故障模式,B为故障模式,则利用两两分类方法建立多分类SVM的输出判决矩阵SM为:
式中,F(Bg,Bh)为故障模式是属于Bi类还是属于Bh类的标准SVM判决的硬输出,若属于Bi类,则F(Bg,Bh)=+1,若属于Bh类,则F(Bg,Bh)=-1,判决矩阵SM第一行表示第B1类与其他类两两分类的SVM硬输出判决结果,判决矩阵第l行表示为第Bl类与其他类两两分类的SVM硬输出判决结果;
对判决矩阵SM的第h行向量可知,第Bh类与其它类进行两两分类的总分类次数为(l-1)次,同时第h行判决硬输出为1的次数之和,即第Bh类与其它类进行两两分类时判决属于第Bj类的次数,将判决属于第Bh类的分类次数与第Bh类参与的总分类次数之比定义为第Bh类的隶属度uh,其表达式为:
将不能确定属于任何一类的隶属度,定义为qΘ:
uΘ=1-max{u1,u2,...,ul}
则基本概率分配m(Bh)的表达式为:
式中,M为SVM多分类的规模;
32)根据证据体的隶属度获取SVM的基本概率分配,并将各证据体按基本概率分配进行加权处理,获取多个证据加权融合后的最终诊断结果:
对于同一辨识框架Θ,假设SVM证据具体对辨识框架中c个命题的可靠度为Re(B)→[0,1],则辨识框架Θ上的可靠度加权系数V(B)为:
进行多证据加权融合,获取可信度W(B)和不确定度W(Θ):
S4:根据获取的可信度和不确定度,对辨识框架中所有的故障模式进行诊断;
设定诊断规则,根据诊断规则对诊断结果进行判断,假设诊断结果为Fc,本发明设定的诊断规则包括:
规则一:
规则二:W(Fc)-W(Fj)>ε;W(Fc)>W(Θ)
规则三:W(Θ)<γ
规则一为诊断的基本条件,即所判定故障模式具有最大的可信度;规则二表明所判定故障模式的可信度必须比其它故障模式的可信度大,其中ε为设定的判定阈值;规则三表明不确定度必须小于阈值γ,保证故障样本是充分可判断的。
本实施例以发电机转子的正常、不平衡、不对中、径向碰磨和油膜涡动五种单一工况故障模式的诊断识别为例来证明本发明方法的有效性及应用型。本实施例中的故障设备共有23台。
建立概率转移矩阵Pp、Pq:
Markov概率转移矩阵谱分析:
计算Pp和Pq特征值得:
λp=[λ1,λ2,λ3,...,λ23]
=[1,0.904,0.798,0.4187,...,0.1135];λq=[λ1,λ2,λ3,...,λ23]
=[1,0.812,0.735,0.554,0.52,...,0.158];
由上可知,Pp、Pq均为前三个特征值间数值差较大,之后的特征值间数值差则相对较小。因此,应该将其分为三个区域。
分区域划分
设定限值η为0.4。对于故障设备13号、6号的划分仍有争议,两故障设备与各个分区域的有功无功距离计算结果分别为:0.4128、0.3583、0.3675、0.4174、0.3897、0.4415;0.4344、0.3609、0.3713、0.4132、0.3941、0.4064。
由算术平均值法可计算出13号故障设备、6号故障设备与各个分区域间的相似距离为:0.3854、0.3925、0.4156;0.3976、0.3922、0.4002。因此13号故障设备应该纳入1号分区域;6号故障设备应该纳入到2号分区域。因此、最后的分配结果如表1所示:
表1分区域划分结果
本实施例以2号分区域作为分析对象,其他分区域的分析方法可同理。由计算可得6、10、11、16、18、23号机组相关系数是0.010、0.38、0.17、0.200、1.028、0.828。分析可知,18、23号相关系数最高。因此,该分区域的标准故障设备为18号、23号,即该分区域的标准故障设备群由18号、23号组成。
***辨识框架为Θ={B1,B2,B3,B4,B5},其中B1为正常状态,B2为不平衡状态,B3为不对中状态,B4为径向碰磨状态,B5为油膜涡动状态。这里的数据主要是布置在轴承座水平和垂直方向上测点得到的振动信号。
分析辨识框架中各故障模式的发生机理,分别从转轴振动信号中提取频域小波能量分析特征向量(λ1,λ2,…,λ6)、时域AR模型自回归参数特征向量(α1,α2,…,α15)和轴心轨迹不变矩特征向量(β1,β2,…,β10)作为三个独立特征向量子空间,从不同侧面对设备进行局部诊断。
基于频域小波能量分析特征向量搭建SVM进行局部诊断作为证据体一,并选取各故障模式下典型样本各100组进行测试,得到该证据体一的混淆矩阵CM1:
基于时域AR模型自回归参数特征向量搭建SVM进行局部诊断作为证据体二,同样选取各故障模式下典型样本各100组进行测试,得到该证据体二的混淆矩阵CM2:
基于轴心轨迹不变矩特征向量搭建SVM进行局部诊断作为证据体三,同样选取各故障模式下典型样本各100组进行测试,得到该证据体三的混淆矩阵CM3:
对混淆矩阵CM1、CM2、CM3进行计算得到三个证据体对转子五种故障模式识别的可靠度如表2所示;取故障模式类别l=5得到各证据体对各状态模式的加权系数如表3所示。
表2各证据体识别的可靠度
表3各证据体的加权系数
现取该转子的一组不平衡状态(B2)信号样本,分别提取其频域、时域和轴心轨迹特征向量,通过各单一特征SVM局部诊断获得各证据体对各故障模式的“一对一”多分类SVM判决(投票)结果统计如表4所示;根据SVM局部诊断获取的证据体的隶属度、基本概率分配分别如表5、如表6所示;加权后的概率分配如表7所示;进行多证据的合成运算后的结果及采用传统S证据理论的计算对比结果如表8所示。
表4“一对一”多分类SVM判决(投票)结果
表5各证据体的隶属度
表6各证据体的概率分配
表7各证据体的加权概率分配
表8多证据体证据融合诊断对比
由表4、表5和表6可以看出,三个SVM单一特征局部诊断证据的判决结论冲突严重,不利于多证据的融合诊断。但从表6、表7的对比可看出,当各证据体被加权处理后,证据G1对B2正确识别的可信度从0.4511提升到0.6499;证据G3对B4错误识别的可信度从0.4511下降到0.2640。由此可见,各SVM单一特征局部诊断证据体通过加权处理后正确识别的可信度能够提升,错误识别的可信度能够下降,各证据体间的冲突明显减小,有利于提高多证据融合诊断的准确率。
从表8中可以看出,传统的D-S证据理论对G1、G2、G3融合后对B2正确识别的可信度和对B4错误识别的可信度都是0.3793,无法对故障模式进行正确诊断;但采用改进后的加权证据理论对多证据体进行加权合成,G1、G2、G3融合后对B2正确识别的可信度提升到了0.6174,对B4错误识别的可信度下降到了0.1596,能够正确诊断该故障模式为B2。由此可见,多证据体通过加权合成处理后能够使得正确识别的结论相互加强,错误识别的结论相互削弱,融合诊断后的结果可信度具有更好的峰值性和可分性,从而能够提高故障诊断***对故障模式的诊断识别能力。
为了评估该融合诊断***的性能,选取不同故障模式下典型样本各100组分别进行单一特征的SVM故障诊断和多特征融合的SVM-DS故障诊断测试,在相同的决策规则下(ε=0.40,γ=0.20),诊断识别结果对比如表9所示。
表9测试样本故障诊断识别结果对比
识别率越高,则诊断越准确,从表9的对比结果可以看出,加权证据融合SVM-DS的识别率达到0.88,远远高于其他方法,可以证明本发明提出的方法确实能够有效提高故障诊断的准确率。本发明提出的谱聚类分析下,结合SVM和改进证据理论的多数据融合故障诊断方法具有很好的应用前景,在实际使用时需要充分结合电厂实际需求,确定需要接入的数据源,根据电厂实际情况完成基于SVM和改进证据理论的多数据融合故障诊断方法的实施。充分发挥电厂工程师在实际设备方面的技术优势,将这些经验转化为***的规则,在使用过程中可以不断扩充规则内容和使用方法,使得整个故障诊断***更智能,更人性,更高效。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的工作人员在本发明揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (10)
1.基于谱聚类分析下多特征信息加权融合的故障诊断方法,其特征在于,该方法包括以下步骤:
1)对各个电厂故障设备进行判断,若某两个故障设备情况集中,则将二者置于同一分区,否则,采用谱聚类方法对故障设备进行划分,获取多个分区域;
2)对各个分区域中的故障设备获取相关系数,并进行正序排序,将绝对值大于设定阈值的相关系数所对应的故障设备作为待诊断目标设备;
3)结合SVM局部诊断及改进的证据理论方法,对待诊断目标设备的故障模式计算基本概率分配并进行加权组合,获取可信度和不确定度;
4)根据可信度和不确定度对辨识框架中所有故障模式进行诊断。
2.根据权利要求1所述的基于谱聚类分析下多特征信息加权融合的故障诊断方法,其特征在于,步骤1)具体包括以下步骤:
11)将多个故障设备划分为同一样本集合;
12)获取样本集合的Markov转移矩阵;
13)对Markov转移矩阵进行谱分析,获取故障设备的划分区域数量;
14)对样本集合中各个点进行扩散距离计算,确定划分区域。
3.根据权利要求2所述的基于谱聚类分析下多特征信息加权融合的故障诊断方法,其特征在于,步骤11)的具体内容为:
假设有n个故障设备,则需要被划分的样本集合为:
X={X1,X2,...,Xi,...,Xn}∈Rm×n
式中,Rm×n为实数集,m为矩阵行数,Xi为第i台故障设备中通过对该机组数据Xi(t)进行采样所构成的参数集合,Xi的表达式为:
Xi={xi(1),xi(2),...,xi(t),,...,xi(m)}T
4.根据权利要求3所述的基于谱聚类分析下多特征信息加权融合的故障诊断方法,其特征在于,步骤12)的具体内容为:
将各个故障设备对应于高斯权重图的各个顶点,获取两点间的相似矩阵Aij,其表达式为:
式中,|xi-xj|2为两点之间的欧几里得距离,σ为尺度参数;
则样本集合的Markov转移矩阵Pij为:
5.根据权利要求4所述的基于谱聚类分析下多特征信息加权融合的故障诊断方法,其特征在于,步骤13)具体包括以下步骤:
131)假设某样本集合的Markov转移矩阵为P,对其进行分解:
式中,λz为矩阵P的第z个特征值,为λz的左特征向量,ψz为λz的右特征向量;
132)若矩阵P有q<n个主要特征值,其中λq<<λq-1,λ1,λ2,...,λq数值相近且趋于数值1,则P的相似矩阵Pq可由前q个特征值和对应特征向量计算;
133)若不符合步骤132)情况,但矩阵P相邻特征值之间的差值骤降,即λq-λq+1>>λk-λk+1,q+1≤k≤n,λq+1值小于0.001,则依旧按照步骤132)处理;
134)对相似矩阵Pq获取谱分解式:
即q为所有故障设备中待划分区域数。
6.根据权利要求5所述的基于谱聚类分析下多特征信息加权融合的故障诊断方法,其特征在于,步骤14)的具体内容为:
定义映射ψ:
其中,ψ为R→R,设任意两故障设备xi与xj之间的扩散距离为:
式中,e1、...、en为x1、...、xn对应的单位向量;
根据λ1,λ2,...,λq及其特征向量,将xi与xj之间的扩散距离简化:
考虑各点间的扩散距离以及分区域数,给出定限值η>0,当D2(xi,xj)≤η时,确定其对应的故障设备划分在同一分区中。
7.根据权利要求6所述的基于谱聚类分析下多特征信息加权融合的故障诊断方法,其特征在于,步骤2)的具体内容为:
21)假设存在某一个含有故障设备数为a的分区域,对其中故障设备分别进行编号:1,2,3,...,a,则第i个故障设备的相关系数ri为:
式中,b为采样数,xik为第i个分区域中,第k个采样点实测参数,yk为包含点k所在分区的实测参数;为第i个分区平均参数,为故障设备b个时刻功率的参数平均值;
22)对步骤21)获取的ri的绝对值按正序排序,将绝对值大且大于设定的阈值的故障设备,设为该分区域中的标准故障设备,具体条件如下:
式中,S为前S个相关性最高故障设备,rall为所有故障设备的相关系数。
8.根据权利要求7所述的基于谱聚类分析下多特征信息加权融合的故障诊断方法,其特征在于,步骤3)的具体内容包括:
31)对待诊断目标设备的故障模式获取特征向量,结合证据理论与SVM局部诊断,获取各故障模式的证据体的隶属度;
假设辨识框架Θ中共有l类故障模式,B为故障模式,则利用两两分类方法建立多分类SVM的输出判决矩阵SM为:
式中,F(Bg,Bh)为故障模式是属于Bi类还是属于Bh类的标准SVM判决的硬输出,若属于Bi类,则F(Bg,Bh)=+1,若属于Bh类,则F(Bg,Bh)=-1,判决矩阵SM第一行表示第B1类与其他类两两分类的SVM硬输出判决结果,判决矩阵第l行表示为第Bl类与其他类两两分类的SVM硬输出判决结果;
对判决矩阵SM的第h行向量可知,第Bh类与其它类进行两两分类的总分类次数为(l-1)次,同时第h行判决硬输出为1的次数之和,即第Bh类与其它类进行两两分类时判决属于第Bj类的次数,将判决属于第Bh类的分类次数与第Bh类参与的总分类次数之比定义为第Bh类的隶属度uh,其表达式为:
将不能确定属于任何一类的隶属度,定义为qΘ:
uΘ=1-max{u1,u2,...,ul}
则基本概率分配m(Bh)的表达式为:
式中,M为SVM多分类的规模;
32)根据证据体的隶属度获取SVM的基本概率分配,并将各证据体按基本概率分配进行加权处理,获取多个证据加权融合后的最终诊断结果。
9.根据权利要求8所述的基于谱聚类分析下多特征信息加权融合的故障诊断方法,其特征在于,步骤32)的具体内容为:
对于同一辨识框架Θ,假设SVM证据具体对辨识框架中c个命题的可靠度为Re(B)→[0,1],则辨识框架Θ上的可靠度加权系数V(B)为:
进行多证据加权融合,获取可信度W(B)和不确定度W(Θ):
10.根据权利要求9所述的基于谱聚类分析下多特征信息加权融合的故障诊断方法,其特征在于,步骤4)的具体内容为:
设定诊断规则,根据诊断规则对诊断结果进行判断,假设诊断结果为Fc,所述的诊断规则包括:
规则一:
规则二:W(Fc)-W(Fj)>ε;W(Fc)>W(Θ)
规则三:W(Θ)<γ
规则一为诊断的基本条件,即所判定故障模式具有最大的可信度;规则二表明所判定故障模式的可信度必须比其它故障模式的可信度大,其中ε为设定的判定阈值;规则三表明不确定度必须小于阈值γ,保证故障样本是充分可判断的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810479820.1A CN108710756A (zh) | 2018-05-18 | 2018-05-18 | 基于谱聚类分析下多特征信息加权融合的故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810479820.1A CN108710756A (zh) | 2018-05-18 | 2018-05-18 | 基于谱聚类分析下多特征信息加权融合的故障诊断方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108710756A true CN108710756A (zh) | 2018-10-26 |
Family
ID=63869010
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810479820.1A Pending CN108710756A (zh) | 2018-05-18 | 2018-05-18 | 基于谱聚类分析下多特征信息加权融合的故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108710756A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109389806A (zh) * | 2018-11-08 | 2019-02-26 | 山东大学 | 基于多信息融合的疲劳驾驶检测预警方法、***及介质 |
CN109888746A (zh) * | 2019-03-06 | 2019-06-14 | 上海电力学院 | 一种基于继电保护信息的电网故障区域识别方法 |
CN110532512A (zh) * | 2019-08-08 | 2019-12-03 | 合肥通用机械研究院有限公司 | 基于大数据分析的设备失效模式诊断特征参量分析方法 |
CN110849419A (zh) * | 2019-11-26 | 2020-02-28 | 国网冀北电力有限公司廊坊供电公司 | 一种开关柜状态监测诊断方法及*** |
CN111781435A (zh) * | 2019-04-04 | 2020-10-16 | 中车唐山机车车辆有限公司 | 四象限整流器的故障检测方法及装置 |
CN112487709A (zh) * | 2020-11-25 | 2021-03-12 | 江苏科技大学 | 一种基于敏感性分析的船用柴油机故障溯源方法 |
CN112613186A (zh) * | 2020-12-30 | 2021-04-06 | 中国航空发动机研究院 | 一种基于统计分布特征的航空发动机气路故障融合诊断方法 |
CN112710486A (zh) * | 2019-10-24 | 2021-04-27 | 广东美的白色家电技术创新中心有限公司 | 设备故障检测方法、设备故障检测装置及计算机存储介质 |
CN113884805A (zh) * | 2021-09-17 | 2022-01-04 | 国网山东省电力公司德州供电公司 | 基于d-s证据理论的单相接地故障绝缘劣化识别方法及*** |
CN114441698A (zh) * | 2020-11-04 | 2022-05-06 | 中国科学院沈阳自动化研究所 | 一种基于贝叶斯网的气相色谱仪故障诊断方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510729A (zh) * | 2014-10-11 | 2016-04-20 | 国家电网公司 | 一种变压器过热故障诊断方法 |
CN105574848A (zh) * | 2014-11-04 | 2016-05-11 | 诺基亚技术有限公司 | 用于对象的自动分割的方法和装置 |
US20170083507A1 (en) * | 2015-09-22 | 2017-03-23 | International Business Machines Corporation | Analyzing Concepts Over Time |
CN107273924A (zh) * | 2017-06-06 | 2017-10-20 | 上海电力学院 | 基于模糊聚类分析的多数据融合的电厂故障诊断方法 |
CN107274011A (zh) * | 2017-06-05 | 2017-10-20 | 上海电力学院 | 综合马尔科夫模型和概率网络的设备状态识别方法 |
CN107368854A (zh) * | 2017-07-20 | 2017-11-21 | 华北电力大学(保定) | 一种基于改进证据理论的断路器故障诊断方法 |
-
2018
- 2018-05-18 CN CN201810479820.1A patent/CN108710756A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105510729A (zh) * | 2014-10-11 | 2016-04-20 | 国家电网公司 | 一种变压器过热故障诊断方法 |
CN105574848A (zh) * | 2014-11-04 | 2016-05-11 | 诺基亚技术有限公司 | 用于对象的自动分割的方法和装置 |
US20170083507A1 (en) * | 2015-09-22 | 2017-03-23 | International Business Machines Corporation | Analyzing Concepts Over Time |
CN107274011A (zh) * | 2017-06-05 | 2017-10-20 | 上海电力学院 | 综合马尔科夫模型和概率网络的设备状态识别方法 |
CN107273924A (zh) * | 2017-06-06 | 2017-10-20 | 上海电力学院 | 基于模糊聚类分析的多数据融合的电厂故障诊断方法 |
CN107368854A (zh) * | 2017-07-20 | 2017-11-21 | 华北电力大学(保定) | 一种基于改进证据理论的断路器故障诊断方法 |
Non-Patent Citations (3)
Title |
---|
向阳辉等: "结合SVM和改进证据理论的多信息融合故障诊断", 《振动与冲击》 * |
谢义苗等: "考虑多因素的大规模设备分区检修决策", 《2017年江西省电机工程学会年会论文集》 * |
赵帆: "基于区域聚类分析的改进神经网络短期风电功率预测研究", 《中国优秀博硕士学位论文全文数据库(硕士)-工程科技Ⅱ辑》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109389806A (zh) * | 2018-11-08 | 2019-02-26 | 山东大学 | 基于多信息融合的疲劳驾驶检测预警方法、***及介质 |
CN109888746A (zh) * | 2019-03-06 | 2019-06-14 | 上海电力学院 | 一种基于继电保护信息的电网故障区域识别方法 |
CN111781435A (zh) * | 2019-04-04 | 2020-10-16 | 中车唐山机车车辆有限公司 | 四象限整流器的故障检测方法及装置 |
CN111781435B (zh) * | 2019-04-04 | 2023-05-23 | 中车唐山机车车辆有限公司 | 四象限整流器的故障检测方法及装置 |
CN110532512A (zh) * | 2019-08-08 | 2019-12-03 | 合肥通用机械研究院有限公司 | 基于大数据分析的设备失效模式诊断特征参量分析方法 |
CN110532512B (zh) * | 2019-08-08 | 2022-07-19 | 合肥通用机械研究院有限公司 | 基于大数据分析的设备失效模式诊断特征参量分析方法 |
CN112710486B (zh) * | 2019-10-24 | 2022-01-25 | 广东美的白色家电技术创新中心有限公司 | 设备故障检测方法、设备故障检测装置及计算机存储介质 |
CN112710486A (zh) * | 2019-10-24 | 2021-04-27 | 广东美的白色家电技术创新中心有限公司 | 设备故障检测方法、设备故障检测装置及计算机存储介质 |
CN110849419A (zh) * | 2019-11-26 | 2020-02-28 | 国网冀北电力有限公司廊坊供电公司 | 一种开关柜状态监测诊断方法及*** |
CN114441698B (zh) * | 2020-11-04 | 2023-03-24 | 中国科学院沈阳自动化研究所 | 一种基于贝叶斯网的气相色谱仪故障诊断方法 |
CN114441698A (zh) * | 2020-11-04 | 2022-05-06 | 中国科学院沈阳自动化研究所 | 一种基于贝叶斯网的气相色谱仪故障诊断方法 |
CN112487709A (zh) * | 2020-11-25 | 2021-03-12 | 江苏科技大学 | 一种基于敏感性分析的船用柴油机故障溯源方法 |
CN112487709B (zh) * | 2020-11-25 | 2023-11-21 | 江苏科技大学 | 一种基于敏感性分析的船用柴油机故障溯源方法 |
CN112613186A (zh) * | 2020-12-30 | 2021-04-06 | 中国航空发动机研究院 | 一种基于统计分布特征的航空发动机气路故障融合诊断方法 |
CN113884805A (zh) * | 2021-09-17 | 2022-01-04 | 国网山东省电力公司德州供电公司 | 基于d-s证据理论的单相接地故障绝缘劣化识别方法及*** |
CN113884805B (zh) * | 2021-09-17 | 2023-10-20 | 国网山东省电力公司德州供电公司 | 基于d-s证据理论的单相接地故障绝缘劣化识别方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108710756A (zh) | 基于谱聚类分析下多特征信息加权融合的故障诊断方法 | |
CN109635928A (zh) | 一种基于深度学习模型融合的电压暂降原因识别方法 | |
CN109543210A (zh) | 一种基于机器学***台的风电机组故障预测*** | |
CN109932627A (zh) | 一种gis耐压试验局部放电的定位方法及*** | |
CN107817404B (zh) | 一种便携式计量自动化终端故障诊断装置及其诊断方法 | |
CN110333414A (zh) | 电力变压器多层次状态评估方法 | |
CN104299115B (zh) | 基于模糊c均值聚类算法的智能变电站二次***状态分析方法 | |
CN110163075A (zh) | 一种基于权值训练的多信息融合故障诊断方法 | |
CN105512448A (zh) | 一种配电网健康指数的评估方法 | |
CN108133225A (zh) | 一种基于支持向量机的覆冰闪络故障预警方法 | |
Di et al. | Ensemble deep transfer learning driven by multisensor signals for the fault diagnosis of bevel-gear cross-operation conditions | |
CN107463993A (zh) | 基于互信息‑核主成分分析‑Elman网络的中长期径流预报方法 | |
CN102818948A (zh) | 基于模糊故障诊断和相关性模型诊断的dr合成诊断方法 | |
CN107257351A (zh) | 一种基于灰色lof流量异常检测***及其检测方法 | |
CN115392333A (zh) | 一种基于改进端到端ResNet-BiLSTM双通道模型的设备故障诊断方法 | |
CN107037278A (zh) | 一种iec61850标准下的智能变电站故障诊断方法 | |
CN106199332A (zh) | 一种改进的轨道电路红光带故障诊断方法 | |
Li et al. | A novel remaining useful life transfer prediction method of rolling bearings based on working conditions common benchmark | |
CN112600659B (zh) | 基于神经网络的安全芯片侧信道泄漏检测方法及*** | |
Min et al. | Fault prediction for distribution network based on CNN and LightGBM algorithm | |
CN117036103A (zh) | 基于lstm循环神经网络的电力***运行态势预测方法 | |
Hassan et al. | Fault classification of power plants using artificial neural network | |
CN114167837B (zh) | 一种铁路信号***的智能故障诊断方法及*** | |
CN115879775A (zh) | 一种三维变电站涉鸟故障风险等级评估方法及*** | |
Dehghani et al. | Distribution feeder classification based on self organized maps (case study: Lorestan province, Iran) |
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: 20181026 |
|
RJ01 | Rejection of invention patent application after publication |