CN107132310A - 基于高斯混合模型的变压器设备健康状态判别方法 - Google Patents
基于高斯混合模型的变压器设备健康状态判别方法 Download PDFInfo
- Publication number
- CN107132310A CN107132310A CN201710193925.6A CN201710193925A CN107132310A CN 107132310 A CN107132310 A CN 107132310A CN 201710193925 A CN201710193925 A CN 201710193925A CN 107132310 A CN107132310 A CN 107132310A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msubsup
- msub
- transformer
- transformer 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 63
- 230000003862 health status Effects 0.000 title claims abstract description 58
- 238000012549 training Methods 0.000 claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000036541 health Effects 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000007476 Maximum Likelihood Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 230000000737 periodic effect Effects 0.000 claims 1
- 238000013459 approach Methods 0.000 abstract description 6
- 230000007812 deficiency Effects 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 59
- 230000006870 function Effects 0.000 description 20
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 13
- 238000012360 testing method Methods 0.000 description 9
- OTMSDBZUPAUEDD-UHFFFAOYSA-N Ethane Chemical compound CC OTMSDBZUPAUEDD-UHFFFAOYSA-N 0.000 description 8
- VGGSQFUCUMXWEO-UHFFFAOYSA-N Ethene Chemical compound C=C VGGSQFUCUMXWEO-UHFFFAOYSA-N 0.000 description 8
- HSFWRNGVRCDJHI-UHFFFAOYSA-N alpha-acetylene Natural products C#C HSFWRNGVRCDJHI-UHFFFAOYSA-N 0.000 description 6
- 239000001257 hydrogen Substances 0.000 description 6
- 229910052739 hydrogen Inorganic materials 0.000 description 6
- 230000007257 malfunction Effects 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 5
- 206010020843 Hyperthermia Diseases 0.000 description 5
- 125000002534 ethynyl group Chemical group [H]C#C* 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- 150000002430 hydrocarbons Chemical class 0.000 description 5
- 230000036031 hyperthermia Effects 0.000 description 5
- 238000012423 maintenance Methods 0.000 description 5
- 150000002431 hydrogen Chemical class 0.000 description 4
- 230000004069 differentiation Effects 0.000 description 3
- 238000004868 gas analysis Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010891 electric arc Methods 0.000 description 2
- 238000013021 overheating Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 239000005977 Ethylene Substances 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 150000001335 aliphatic alkanes Chemical class 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000001684 chronic effect Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000003449 preventive effect Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/12—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
- G01R31/1227—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Combustion & Propulsion (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Housings And Mounting Of Transformers (AREA)
Abstract
本发明公开了一种数据驱动的变压器设备健康状态判别方法，该方法利用不同规格型号、故障类型的变压器及其油中溶解气体数据构建训练样本，通过高斯混合模型描述样本中多种特征气体在高维空间内的分布情况，并运用所得模型的参数对变压器设备的健康状态进行判别。该方法弥补了现有比值方法在描述复杂特征方面的不足，较人工经验相比具有较高的严谨性和创新性，能够对变压器设备的运行状态进行更为准确的判断。
Description
技术领域
本发明涉及统计分析领域，具体涉及一种基于高斯混合模型的变压器设备健康状态判别方法。
背景技术
变压器设备发生故障时往往会释放较高的电和热，导致变压器油分解，产生不同的油中溶解气体，包括氢气、甲烷、乙烷、乙烯等。油中溶解气体分析法通过计算、比较变压器在运行过程中产生气体的浓度和组成成分，分析设备存在的潜在风险，对变压器设备的健康状态给予初步的判断。油中溶解气体分析法是一种非常重要的变压器故障诊断方法，该方法对判断变压器慢性局部缺陷十分有效，具有方便在线监测、无需拆解检测的优点，电力工业部1997年颁布的《电力设备预防性试验规程(DL/T 596—1996)》中将油中溶解气体分析置于首要地位。
油中溶解气体是变压器故障的重要判断指标。根据大量的实际经验，相关研究人员已总结出了一些针对不同故障类型体的规律和判断方法，如IEC 60599(三比值)方法、Rogers(四比值)方法和Doernenburg(五比值)方法等。这些方法的思路相同，都是利用变压器绝缘油产生气体的浓度和比值来判断变压器的故障类型，是一种来源于实际经验阈值方法。以IEC 60599(三比值)方法为例，该方法使用“乙炔/乙烯”、“甲烷/氢气”、“乙烯/乙烷”三种比值来判断变压器设备的故障类型。例如，当甲烷与氢气的比值小于0.1且乙烯与乙烷的比值小于0.2时，该方法将变压器设备的状态判别为局部放电故障。IEC 60599(三比值)详细的故障判别标准如表1所示。
目前，我国推行的《变压器油中溶解气体分析和判断导则(国家电力行业标准DL/T722-2014)》中针对变压器设备故障类型判断的气体含量比值法，是在IEC60599(三比值)方法基础上的改进方法。该方法共进行了两方面的改进，一是在气体浓度和增长速率方面重新设置了注意值(见表2)，二是根据国内的实践经验对比值的编码和故障类型进行了细化(见表3、4)。改进后的方法先将待测变压器设备的气体数据与注意值进行对照，当待测数据超过注意值时，相关变压器维护人员应缩短设备的检测周期，并结合气体比值及对应的编码对设备的故障类型进行初步的判断。
表1IEC 60599三比值法解释表
| 故障类型 | C2H2/C2H4 | CH4/H2 | C2H4/C2H6 |
| 局部放电 | NS | <0.1 | <0.2 |
| 低能量放电 | >1 | 0.1～0.5 | >1 |
| 高能量放电 | 0.6～2.5 | 0.1～1 | >2 |
| 低温过热(低于300℃) | NS | >1但NS | <1 |
| 中温过热(300℃～700℃) | <0.1 | >1 | 1～4 |
| 高温过热(高于700℃) | <0.2 | >1 | >4 |
注：表中NS表示不重要。
表2变压器油中溶解气体含量注意值与绝对产气速率注意值
表3比值方法编码原则
表4变压器故障类型判断方法
从现有的变压器设备运行状态判别方法来看，IEC 60599(三比值)方法、Rogers(四比值)方法和Doernenburg(五比值)方法属于早期应用的经典方法，这些方法以气体比值作为设备故障的判断依据，能够在不对变压器设备进行解体的前提下对设备的运行状态进行初步的判断。然而，这些比值只有在设备发生故障后才具有意义，在设备正常的情况下极易造成误判，且应用当中容易产生比值计算结果与故障判别表无法对应的情况，因此在实践方面仍存在一定的问题。我国目前推行的《变压器油中溶解气体分析和判断导则(国家电力行业标准DL/T 722-2014)》在上述比值方法的基础上设置了注意值，减少了设备正常情况下的误判现象，并对比值编码和故障类型进行了细化，排除了比值结果与故障类型无法对应的情况，为变压器设备运行状态的判别提供了可靠的依据和标准。但是，此类比值方法仍存在一定的局限，具体表现在三个方面，第一，变压器设备故障的实际情况十分复杂，仅用简单的比值的形式对其进行描述存在一定的欠缺，比值运算无法体现出气体指标之间的复杂特征。第二，现有技术仅考虑设备的正常和故障两种状态，而当设备已经出现故障状态时，相关维护人员只能对其进行事后的维修工作，成本较高。变压器故障识别方法应体现出设备亚健康状态方面的信息，做到“早发现，早预防”，降低维护成本。第三，比值方法的本质是一种阈值方法(例如上述编码原则以0.1、1、3三种阈值对气体比值进行划分)，而这些阈值大多来自研究人员的实践经验，缺乏一定的严谨性，阈值附近的故障判别往往比较模糊，此类方法在判别的准确度方面还有进一步的提升空间。
发明内容
鉴于上述，本发明提供了一种基于高斯混合模型的变压器设备健康状态判别方法，具体是一种通过高斯混合模型进行聚类分析，并在此基础上完成健康状态判别的方法。该方法弥补了现有比值方法在描述复杂特征方面的不足，较人工经验相比具有较高的严谨性和创新性，能够对变压器设备的运行状态进行更为准确的判断。
一种基于高斯混合模型的变压器设备健康状态判别方法，包括以下步骤：
(1)利用故障变压器设备与正常变压器设备的的油中溶解气体历史浓度数据对高斯混合模型进行训练，得到健康状态模型Mm,s与正常状态模型M，其中，m代表变压器设备的故障类型，s代表变压器设备的规格型号；
(2)利用健康状态模型Mm,s与正常状态模型M，找到与预处理后的待测变压器的油中溶解气体浓度数据相匹配的模型M′，并根据模型M′的参数计算得到变压器当前时刻的健康状态。
所述变压器油中溶解气体包括氢气、甲烷、乙烷、乙烯、乙炔和总烃，其中，总烃是指甲烷、乙烷、乙烯、乙炔四种烃类气体含量总和；所述气体浓度数据为气体的体积浓度，单位为μL/L。
所述故障变压器的油中溶解气体历史浓度数据是指设备从运行良好的状态开始，逐渐向故障状态发展，直到设备出现持续的故障状态为止的过程，相应的气体浓度数据是一个时间序列，时间跨度设为定值；所述正常变压器的油中溶解气体历史浓度数据是指从投运开始至今未出现任何异常的变压器设备，其在正常运行情况下产生的油中溶解气体浓度数据，该数据同样为时间序列，时间跨度与故障变压器设备的气体浓度数据的时间跨度相同。
所述的步骤(1)的具体过程为：
(1-1)采集故障变压器设备与正常变压器设备的油中溶解气体历史浓度数据，并对其进行预处理，得到故障变压器设备数据库与正常变压器设备的数据库；
(1-2)对故障变压器设备数据库中的数据按照数据所属变压器设备的故障类型进行分类，得到针对不同故障类型的数据集Cm，代表变压器设备的故障类型，然后对数据集Cm中的数据按照其所属的规格型号进行分类，得到数据集Cm,s；
(1-3)将数据集Cm,s作为高斯混合模型的训练样本，利用高斯混合模型对数据集Cm,s进行聚类分析，聚类数目设置为3，对应设备正常、亚健康、故障异常三种状态，并以似然函数L收敛为目标，确定高斯混合模型模型的最优参数得到健康状态模型Mm,s；
(1-4)将正常变压器设备的数据库作为高斯混合模型的训练样本，利用高斯混合模型对该数据库进行聚类分析，以最大似然函数L收敛为目标，确定高斯混合模型模型的最优参数得到正常状态模型M。
所述变压器设备故障类型为《变压器油中溶解气体分析和判断导则(DL/T 722-2014)》中气体比值法运用的故障类型，包括低温过热、中温过热、高温过热、局部放电、低能放电、低能放电兼过热、电弧放电和电弧放电兼过热，共8种。
所述的变压器设备的规格型号指变压器的电压等级，包括110kV、220kV、330kV等。
步骤(1-1)中，所述对油中溶解气体历史浓度数据进行预处理采用线性插值方法，具体过程为：以1天为单位对所有样本进行线性插值处理，补全数据中的缺失值，得到时间间隔为1天的气体浓度数据插值序列。再按照事先设定的时间间隔(如1周)对其进行等间隔取样，使得所有训练样本的时间间隔相同。
所述高斯混合模型(Gaussian Mixture Model,GMM)是一种利用多个高斯函数精确量化事物的模型，该模型可以通过聚类分析的方式将训练样本中的数据自发式的分类，每一类对应设备的不同健康状态，模型通过概率分布来反映数据在高维空间的分布情况。
步骤(1-3)中，所述的似然函数为：
其中，n表示第n台变压器设备，N表示变压器设备的总个数，t表示时刻，时间长度为T，为算法迭代前的初始参数，π′i,μ′i,∑′i为算法迭代后的更新参数，更新目标为最大化似然函数值L。表示设备在第i种健康状态时气体数据为xnt的概率，由计算所得。φi为设备在第i种健康状态时气体数据的概率密度函数。
在对模型训练的过程中，高斯混合模型将故障变压器设备的训练样本的气体浓度数据自发的分为三类，气体浓度较低的一类对应变压器设备正常状态，气体浓度较高的一类对应变压器设备故障状态，介于二者之间的一类对应变压器设备亚健康状态，模型相应地产生三种参数。对于正常变压器设备的训练样本，由于样本只包含变压器正常状态的信息，高斯混合模型将所有数据视为一类，对应变压器设备正常状态，模型只产生一种参数。
所述步骤(2)的具体步骤为：
(2-1)采集待测变压器设备当前时刻的油中溶解气体浓度数据，并对其进行线性插值处理，得到待检测样本X＝(X1…X6)；
(2-2)利用待检测样本X＝(X1…X6)、健康状态模型Mm,s及其最优参数正常状态模型M及其最优参数计算似然函数值f，计算得似然函数值序列f1,f2,…,fd×h,fd×h+1，MAX(f1,f2,…,fd×h,fd×h+1)所对应的模型为待检测数据O′匹配的模型M′其中，d为变压器设备故障类型的总个数，h为变压器设备规格型号的总个数，似然函数f为：
其中，φi为设备在第i种健康状态时气体数据的概率密度函数。
(2-3)利用待检测样本X＝(X1…X6)、匹配模型M′及其最优参数计算待检测变压器设备在不同健康状态下的概率，具体计算公式为：
概率最大的健康状态为待检测变压器设备当前时刻的健康状态。
本发明方法通过高斯混合模型描述多种特征气体在高维空间内的分布情况，弥补了比值方法在描述复杂特征方面的局限。同时，与人工经验相比，该方法从数据本身挖掘出具有价值的故障信息，具有较高的严谨性和创新性，能够变压器设备运行的健康状态进行准确的判断。
附图说明
图1是本发明基于高斯混合模型的变压器设备健康状态判别方法的流程图。
具体实施方式
为了更为具体地描述本发明，下面结合附图及具体实施方式对本发明的技术方案进行详细说明。
如图1所示，本发明基于高斯混合模型的变压器设备健康状态判别方法，包括：
S01，采集故障变压器设备与正常变压器设备的油中溶解气体历史浓度数据，并对其进行线性插值处理，得到故障变压器设备数据库与正常变压器设备的数据库；
S02，对故障变压器设备数据库中的数据按照数据所属变压器设备的故障类型进行分类，得到针对不同故障类型的数据集Cm，然后对数据集Cm中的数据按照其所属的规格型号进行分类，得到数据集Cm,s；
S03，将数据集Cm,s作为高斯混合模型的训练样本，利用高斯混合模型对数据集Cm,s进行聚类分析，聚类数目设置为3，对应设备正常、亚健康、故障异常三种状态，并以似然函数L收敛为目标，确定高斯混合模型模型的最优参数得到健康状态模型Mm,s；
S04，将正常变压器设备的数据库作为高斯混合模型的训练样本，利用高斯混合模型对该数据库进行聚类分析，以最大似然函数L收敛为目标，确定高斯混合模型模型的最优参数得到正常状态模型M。
高斯混合模型是一种利用多个高斯函数精确量化事物的模型，该模型可以通过聚类分析的方式将训练样本中的数据自发式的分类，每一类对应设备的不同健康状态，模型通过概率分布来反映数据在高维空间的分布情况。
对于某种类型的训练样本，设训练样本中共有N台变压器，每台变压器的气体浓度数据拥有T个时间点，第n台变压器t时刻的气体数据记为n＝1,2…N，t＝1,2…T，分别对应6种油中溶解气体(氢气、甲烷、乙烷、乙烯、乙炔和总烃)。高斯混合模型通过联合概率密度来描述气体数据在高维空间的分布情况，具体形式为：
式中，i表示样本的健康状态(正常、亚健康、故障共有k种)。μi，∑i，πi为模型的参数，分别对应第i种健康状态气体数据的均值、协方差矩阵以及训练样本中第i种健康状态出现的概率。φi为第i种健康状态气体数据的概率密度函数，计算方式如下：
高斯混合模型通过最大期望算法(EM)算法对训练样本中的数据进行分类，并对模型中的参数进行估计，具体步骤如下：
(1)随机赋予模型中各参数一个初始值，
(2)计算数据xnt属于第i种状态的概率具体计算方式如下：
(3)利用数据xnt与状态概率对模型中的参数进行更新，得到参数的更新结果μi′，∑′i和π′i，具体计算方式如下所示。同时，将所得更新结果μi′，∑′I和π′i代入似然函数L，并计算相应似然函数值L。
(4)将步骤三中得到的参数更新值μi′,∑′i,π′i代替步骤(2)中的参数值不断重复步骤(2)与步骤(3)，直至步骤(3)中所得的似然函数值L收敛，此时的参数值μi′,∑′i,π′i即为模型的最优参数结果，记为
(5)根据步骤(4)中所得的最优参数计算训练样本中每个数据的状态概率以最大概率为原则对数据进行分类。对于某种类型的故障样本，模型将样本数据分为三类，气体浓度较低的一类对应变压器设备的正常状态，气体浓度较高的一类对应变压器设备的故障状态，介于二者之间的一类对应变压器设备的亚健康状态。对于某种类型的正常样本，模型将样本中的数据视为一类，对应变压器设备的正常状态。
S05，采集待测变压器设备当前时刻的油中溶解气体浓度数据，并对其进行线性插值处理，得到待检测样本X＝(X1…X6)；
S06，利用待检测样本X＝(X1…X6)、健康状态模型Mm,s及其最优参数正常状态模型M及其最优参数计算似然函数值f，计算得似然函数值序列f1,f2,…,fd×h,fd×h+1，MAX(f1,f2,…,fd×h,fd×h+1)所对应的模型为待检测数据O′匹配的模型M′其中，d为变压器设备故障类型的总个数，h为变压器设备规格型号的总个数，似然函数f为：
S07，利用待检测样本X＝(X1…X6)、匹配模型M′及其最优参数计算待检测变压器设备在不同健康状态下的概率，具体计算公式为：
概率最大的健康状态为待检测变压器设备当前时刻的健康状态。
实施例1
首先，本实施例收集了6台故障变压器设备(故障类型为高温过热)与32台正常变压器设备的变压器油中溶解气体数据，设备的电压等级为220kV，气体类型包括氢气、甲烷、乙烷、乙烯、乙炔和总烃，时间跨度为1600天，组成了针对故障类型为高温过热的故障设备的训练样本和正常设备的训练样本。同时，按照本发明中提出的数据预处理方法，利用线性插值法以1天为单位对所有训练样本中每台变压器气体浓度数据进行线性插值处理，补全数据中的缺失值，得到时间间隔为1天的气体浓度数据插值序列。再以20天为时间间隔等间隔取样，使得所有训练样本中的变压器气体浓度数据时间间隔均为20天。
其次，对于故障设备的训练样本，利用高斯混合模型对其进行聚类分析，将训练样本中的数据分为三类，对应设备的正常、亚健康、故障三种状态，相应的参数估计结果如下所示：
所得最优参数中，第一类对应的气体浓度均值最低(最小)，对应设备的正常状态，第三类对应的气体浓度均值最高(最大)，对应设备的故障状态，介于二者之间的第二类对应设备的亚健康状态。
对于正常设备的训练样本，由于样本中的数据只包含变压器正常状态的信息，高斯混合模型将所有数据视为一类，对应设备的正常状态，相应的计算结果如下。
最后，本实施例选取了3台220kV的高温过热故障设备作为测试设备，并对其健康状态进行判别。变压器设备健康状态判别方法分为两个部分，一是似然度计算，即利用待测设备的油中溶解气体数据与模型的最优参数计算不同训练样本对应的高斯混合模型中的似然函数值，筛选出与待测设备实际情况最为相似的模型。二是健康状态判别，即利用筛选出的模型及对应的最优参数计算待测设备油中溶解气体数据处于不同健康状态下的概率。三台待测设备的油中溶解气体数据、设备的真实健康状态、似然函数计算结果与设备健康状态判别结果分别如表1、2、3所示。从判别结果来看，三台测试设备真实健康状态为故障的数据范围与判别结果为故障的数据范围完全重合，本发明提供的变压器健康状态判别方法能够准确的识别出设备的故障状态。同时，该方法识别出的亚健康状态往往对应设备发生故障之前的过渡阶段，能够对设备的保养与维护起到一定的警示作用。
表1测试设备1健康状态判别结果
表2测试设备2健康状态判别结果
表3测试设备3健康状态判别结果
以上所述的具体实施方式对本发明的技术方案和有益效果进行了详细说明，应理解的是以上所述仅为本发明的最优选实施例，并不用于限制本发明，凡在本发明的原则范围内所做的任何修改、补充和等同替换等，均应包含在本发明的保护范围之内。
Claims (5)
1.一种基于高斯混合模型的变压器设备健康状态判别方法，包括以下步骤：
(1)利用故障变压器设备与正常变压器设备的的油中溶解气体历史浓度数据对高斯混合模型进行训练，得到健康状态模型Mm,s与正常状态模型M，其中，m代表变压器设备的故障类型，s代表变压器设备的规格型号；
(2)利用健康状态模型Mm,s与正常状态模型M，找到与预处理后的待测变压器的油中溶解气体浓度数据相匹配的模型M′，并根据模型M′的参数计算得到变压器当前时刻的健康状态。
2.如权利要求1所述的基于高斯混合模型的变压器设备健康状态判别方法，其特征在于，所述的步骤(1)的具体过程为：
(1-1)采集故障变压器设备与正常变压器设备的油中溶解气体历史浓度数据，并对其进行预处理，得到故障变压器设备数据库与正常变压器设备的数据库；
(1-2)对故障变压器设备数据库中的数据按照数据所属变压器设备的故障类型进行分类，得到针对不同故障类型的数据集Cm，代表变压器设备的故障类型，然后对数据集Cm中的数据按照其所属的规格型号进行分类，得到数据集Cm,s；
(1-3)将数据集Cm,s作为高斯混合模型的训练样本，利用高斯混合模型对数据集Cm,s进行聚类分析，聚类数目设置为3，对应设备正常、亚健康、故障异常三种状态，并以似然函数L收敛为目标，确定高斯混合模型模型的最优参数得到健康状态模型Mm,s；
(1-4)将正常变压器设备的数据库作为高斯混合模型的训练样本，利用高斯混合模型对该数据库进行聚类分析，以最大似然函数L收敛为目标，确定高斯混合模型模型的最优参数得到正常状态模型M。
3.如权利要求2所述的基于高斯混合模型的变压器设备健康状态判别方法，其特征在于，所述对油中溶解气体历史浓度数据进行预处理采用线性插值方法，具体过程为：
以1天为单位，对所有设备案例库中每台变压器的溶解气体浓度数据进行线性插值处理，补全数据中的缺失值，得到时间间隔为1天的气体浓度数据插值序列。再按照事先设定的时间间隔对气体浓度数据插值序列进行等间隔取样，使得所有设备案例库中的变压器的溶解气体浓度数据的时间间隔相同。
4.如权利要求2所述的基于高斯混合模型的变压器设备健康状态判别方法，其特征在于，所述的似然函数为：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>L</mi>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</msubsup>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</msubsup>
<msub>
<mi>p</mi>
<mrow>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
<msubsup>
<mi>&mu;</mi>
<mi>i</mi>
<mn>0</mn>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&Sigma;</mi>
<mi>i</mi>
<mn>0</mn>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&pi;</mi>
<mi>i</mi>
<mn>0</mn>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mi>log</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&phi;</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
<msup>
<msub>
<mi>&mu;</mi>
<mi>i</mi>
</msub>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<msubsup>
<mi>&Sigma;</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>T</mi>
</msubsup>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</msubsup>
<msub>
<mi>p</mi>
<mrow>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
<msubsup>
<mi>&mu;</mi>
<mi>i</mi>
<mn>0</mn>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&Sigma;</mi>
<mi>i</mi>
<mn>0</mn>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&pi;</mi>
<mi>i</mi>
<mn>0</mn>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mi>log</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&pi;</mi>
<mi>i</mi>
<mo>&prime;</mo>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，n表示第n台变压器设备，N表示变压器设备的总个数，t表示时刻，时间长度为T，为算法迭代前的初始参数，π′i,μ′i,∑′i为算法迭代后的更新参数，更新目标为最大化似然函数值L，表示设备在第i种健康状态时气体数据为xnt的概率，由计算所得，φi为设备在第i种健康状态时气体数据的概率密度函数。
5.如权利要求1所述的基于高斯混合模型的变压器设备健康状态判别方法，其特征在于，所述步骤(2)的具体步骤为：
(2-1)采集待测变压器设备当前时刻的油中溶解气体浓度数据，并对其进行线性插值处理，得到待检测样本X＝(X1…X6)；
(2-2)利用待检测样本X＝(X1…X6)、健康状态模型Mm,s及其最优参数正常状态模型M及其最优参数计算似然函数值f，计算得似然函数值序列f1,f2,…,fd×h,fd×h+1，MAX(f1,f2,…,fd×h,fd×h+1)所对应的模型为待检测数据O′匹配的模型M′其中，d为变压器设备故障类型的总个数，h为变压器设备规格型号的总个数，似然函数f为：
<mrow>
<mi>f</mi>
<mo>=</mo>
<mi>l</mi>
<mi>o</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</msubsup>
<msubsup>
<mi>&pi;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<msub>
<mi>&phi;</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>X</mi>
<mn>1</mn>
</msub>
<mo>...</mo>
<msub>
<mi>X</mi>
<mn>6</mn>
</msub>
<mo>|</mo>
<msubsup>
<mi>&mu;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&Sigma;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
其中，φi为设备在第i种健康状态时气体数据的概率密度函数；
(2-3)利用待检测样本X＝(X1…X6)、匹配模型M′及其最优参数计算待检测变压器设备在不同健康状态下的概率，具体计算公式为：
<mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>X</mi>
<mn>1</mn>
</msub>
<mn>...</mn>
<msub>
<mi>X</mi>
<mn>6</mn>
</msub>
<mo>|</mo>
<msubsup>
<mi>&mu;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&Sigma;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&pi;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&pi;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<msub>
<mi>&phi;</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>X</mi>
<mn>1</mn>
</msub>
<mn>...</mn>
<msub>
<mi>X</mi>
<mn>6</mn>
</msub>
<mo>|</mo>
<msubsup>
<mi>&mu;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&Sigma;</mi>
<mi>i</mi>
<mo>*</mo>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</msubsup>
<msubsup>
<mi>&pi;</mi>
<mi>j</mi>
<mo>*</mo>
</msubsup>
<msub>
<mi>&phi;</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>X</mi>
<mn>1</mn>
</msub>
<mn>...</mn>
<msub>
<mi>X</mi>
<mn>6</mn>
</msub>
<mo>|</mo>
<msubsup>
<mi>&mu;</mi>
<mi>j</mi>
<mo>*</mo>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&Sigma;</mi>
<mi>j</mi>
<mo>*</mo>
</msubsup>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
概率最大的健康状态为待检测变压器设备当前时刻的健康状态。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710193925.6A CN107132310B (zh) | 2017-03-28 | 2017-03-28 | 基于高斯混合模型的变压器设备健康状态判别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710193925.6A CN107132310B (zh) | 2017-03-28 | 2017-03-28 | 基于高斯混合模型的变压器设备健康状态判别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107132310A true CN107132310A (zh) | 2017-09-05 |
| CN107132310B CN107132310B (zh) | 2019-04-30 |
Family
ID=59721684
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710193925.6A Expired - Fee Related CN107132310B (zh) | 2017-03-28 | 2017-03-28 | 基于高斯混合模型的变压器设备健康状态判别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107132310B (zh) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109299844A (zh) * | 2018-07-03 | 2019-02-01 | 国网浙江省电力有限公司电力科学研究院 | 一种电力设备状态静态阈值评估方法 |
| CN109376881A (zh) * | 2018-12-12 | 2019-02-22 | 中国航空工业集团公司上海航空测控技术研究所 | 基于维修成本优化的复杂系统维修决策方法 |
| CN109656202A (zh) * | 2018-12-06 | 2019-04-19 | 南通大学 | 基于局部和结构一致性流形gmm的多模态过程监测方法 |
| CN109738782A (zh) * | 2018-12-07 | 2019-05-10 | 上海交通大学 | 基于测量数据机器学习的电路板故障检测方法和系统 |
| CN110222980A (zh) * | 2019-06-05 | 2019-09-10 | 上海电气集团股份有限公司 | 轨道交通轴承的健康评估方法及系统 |
| CN110929768A (zh) * | 2019-11-14 | 2020-03-27 | 国电大渡河检修安装有限公司 | 一种机器故障的预测方法 |
| CN110988563A (zh) * | 2019-12-23 | 2020-04-10 | 厦门理工学院 | 一种ups故障检测方法、装置、设备和存储介质 |
| CN111597655A (zh) * | 2020-05-19 | 2020-08-28 | 西安航天动力技术研究所 | 一种基于故障发生概率的产品健康判断方法 |
| CN115758260A (zh) * | 2023-01-05 | 2023-03-07 | 青岛明思为科技有限公司 | 基于高斯混合模型的机械设备状态检测的方法 |
| CN115982602A (zh) * | 2023-03-20 | 2023-04-18 | 济宁众达利电气设备有限公司 | 一种光伏变压器电故障检测方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101299183B1 (ko) * | 2012-06-22 | 2013-08-22 | 성균관대학교산학협력단 | 음성 처리 기법을 이용한 전력품질 분류 방법, 전력품질 모니터링 장치 및 전력품질 모니터링 시스템 |
| CN105335759A (zh) * | 2015-11-12 | 2016-02-17 | 南方电网科学研究院有限责任公司 | 一种基于生成概率模型的变压器故障检测方法 |
| CN105930901A (zh) * | 2016-07-18 | 2016-09-07 | 河海大学 | 一种基于rbpnn的变压器故障诊断方法 |
-
2017
- 2017-03-28 CN CN201710193925.6A patent/CN107132310B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101299183B1 (ko) * | 2012-06-22 | 2013-08-22 | 성균관대학교산학협력단 | 음성 처리 기법을 이용한 전력품질 분류 방법, 전력품질 모니터링 장치 및 전력품질 모니터링 시스템 |
| CN105335759A (zh) * | 2015-11-12 | 2016-02-17 | 南方电网科学研究院有限责任公司 | 一种基于生成概率模型的变压器故障检测方法 |
| CN105930901A (zh) * | 2016-07-18 | 2016-09-07 | 河海大学 | 一种基于rbpnn的变压器故障诊断方法 |
Non-Patent Citations (2)
| Title |
|---|
| 刘娜等: "基于组合神经网络模型的电力变压器故障诊断方法", 《电工技术学报》 * |
| 徐永军等: "基于Borda模型的变压器故障诊断系统的研究", 《电工电气》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109299844A (zh) * | 2018-07-03 | 2019-02-01 | 国网浙江省电力有限公司电力科学研究院 | 一种电力设备状态静态阈值评估方法 |
| CN109656202A (zh) * | 2018-12-06 | 2019-04-19 | 南通大学 | 基于局部和结构一致性流形gmm的多模态过程监测方法 |
| CN109738782A (zh) * | 2018-12-07 | 2019-05-10 | 上海交通大学 | 基于测量数据机器学习的电路板故障检测方法和系统 |
| CN109376881B (zh) * | 2018-12-12 | 2022-06-03 | 中国航空工业集团公司上海航空测控技术研究所 | 基于维修成本优化的复杂系统维修决策方法 |
| CN109376881A (zh) * | 2018-12-12 | 2019-02-22 | 中国航空工业集团公司上海航空测控技术研究所 | 基于维修成本优化的复杂系统维修决策方法 |
| CN110222980A (zh) * | 2019-06-05 | 2019-09-10 | 上海电气集团股份有限公司 | 轨道交通轴承的健康评估方法及系统 |
| CN110929768A (zh) * | 2019-11-14 | 2020-03-27 | 国电大渡河检修安装有限公司 | 一种机器故障的预测方法 |
| CN110988563A (zh) * | 2019-12-23 | 2020-04-10 | 厦门理工学院 | 一种ups故障检测方法、装置、设备和存储介质 |
| CN110988563B (zh) * | 2019-12-23 | 2022-04-01 | 厦门理工学院 | 一种ups故障检测方法、装置、设备和存储介质 |
| CN111597655A (zh) * | 2020-05-19 | 2020-08-28 | 西安航天动力技术研究所 | 一种基于故障发生概率的产品健康判断方法 |
| CN111597655B (zh) * | 2020-05-19 | 2023-04-21 | 西安航天动力技术研究所 | 一种基于故障发生概率的产品健康判断方法 |
| CN115758260A (zh) * | 2023-01-05 | 2023-03-07 | 青岛明思为科技有限公司 | 基于高斯混合模型的机械设备状态检测的方法 |
| CN115982602A (zh) * | 2023-03-20 | 2023-04-18 | 济宁众达利电气设备有限公司 | 一种光伏变压器电故障检测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107132310B (zh) | 2019-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107132310A (zh) | 基于高斯混合模型的变压器设备健康状态判别方法 | |
| CN103136587B (zh) | 基于支持向量机的配电网运行状态分类识别方法 | |
| CN106156485A (zh) | 电力变压器故障诊断方法及装置 | |
| CN103245861B (zh) | 一种基于贝叶斯网络的变压器故障诊断方法 | |
| CN101587155B (zh) | 一种油浸式变压器的故障诊断方法 | |
| CN107862114A (zh) | 基于三比值特征量的小波pso‑svm变压器故障诊断方法 | |
| CN109145949A (zh) | 基于集成学习的非侵入式电力负荷监测与分解方法及系统 | |
| CN104020401B (zh) | 基于云模型理论的变压器绝缘热老化状态的评估方法 | |
| CN104535865A (zh) | 基于多参数的电力变压器运行故障综合诊断方法 | |
| CN103745119A (zh) | 一种基于故障概率模型的油浸式变压器故障诊断方法 | |
| CN105242155A (zh) | 一种基于熵权法和灰色关联分析的变压器故障诊断方法 | |
| CN106529124A (zh) | 基于主成分分析与支持向量机的变压器绝缘状态评估方法 | |
| CN110133410A (zh) | 基于模糊c均值聚类算法的变压器故障诊断方法及系统 | |
| CN101738998B (zh) | 一种基于局部判别分析的工业过程监测系统及方法 | |
| CN103198175A (zh) | 基于模糊聚类的变压器故障诊断方法 | |
| CN106770939A (zh) | 一种基于支持向量描述和k近质心近邻的变压器故障诊断方法 | |
| CN107037306A (zh) | 基于隐马尔科夫模型的变压器故障动态预警方法 | |
| CN102289682A (zh) | 基于集成学习Bagging算法的变压器故障诊断方法 | |
| CN113342784A (zh) | 一种电网主变设备风险评估的数据库设计方法 | |
| CN104123678A (zh) | 一种基于状态等级评估模型的电力继电保护状态检修方法 | |
| CN108874984A (zh) | 一种对质量较差电网设备缺陷文本的质量提升方法 | |
| CN111797566A (zh) | 一种表征变压器健康状态的关键特征量确定方法和系统 | |
| CN110689068A (zh) | 一种基于半监督svm的变压器故障类型诊断方法 | |
| CN110503305A (zh) | 一种变压器性能评估方法 | |
| CN110596490A (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190430 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |