CN112504435A - 一种高压断路器基座螺丝松动故障检测方法 - Google Patents
一种高压断路器基座螺丝松动故障检测方法 Download PDFInfo
- Publication number
- CN112504435A CN112504435A CN202011316591.5A CN202011316591A CN112504435A CN 112504435 A CN112504435 A CN 112504435A CN 202011316591 A CN202011316591 A CN 202011316591A CN 112504435 A CN112504435 A CN 112504435A
- Authority
- CN
- China
- Prior art keywords
- vibration
- operator
- fdexc
- circuit breaker
- voltage circuit
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H17/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
-
- 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/327—Testing of circuit interrupters, switches or circuit-breakers
- G01R31/3271—Testing of circuit interrupters, switches or circuit-breakers of high voltage or medium voltage devices
- G01R31/3275—Fault detection or status indication
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种高压断路器基座螺丝松动故障检测方法,引入两个震动传感器,构建了加权分段震动差异算子来衡量两个传感器之间震动差异,进而衡量震动差异数据与高压断路器正常运行状态下的距离来判断高压断路器基座螺丝是否存在松动故障。利用本发明对于两个震动传感器的安装位置和高压断路器连接基座的材质没有严格的要求,降低设备装配条件限制;同时在数据上本发明仅需要收集震动传感器正常运行时的数据,不会与具体的基座螺丝松动样本绑定,提高检测成功率,保障高压断路器安全稳定运行。
Description
技术领域
本发明涉及电力系统中高压断路器故障检测方法,具体涉及一种高压断路器基座螺丝松动故障检测方法,属于高压电力系统安全防护装置技术领域。
背景技术
高压断路器能切断或闭合高压线路,作为电力系统的核心器件之一它起着对线路的控制和保护的关键作用。高压断路器其安全运行、维护及检修与整个电网系统的正常运行密切相关。高压断路器操作机构易于出现机械故障,比如在长时间操作过程中震动导致基座螺丝松动;这些螺丝的松动会进一步增强高压断路器的震动降低其使用寿命,甚至影响高压断路器的切断或闭合操作导致开关拒动,因此及时准确的检查高压断路器基座螺丝是否存在松动情况是十分必要的。
由于高压断路器所处位置存在危险性,所以通常不适合高频次的人工检查;对于高压断路器基座螺丝的松动的自动检测,目前采取的主流方法是利用震动传感器收集基座螺丝的震动数据,并基于震动数据的特征利用人工智能模型进行预测,此类方法在实验室进行模拟实验通常较为成功,但是在实际实施过程中会产生问题:(1)由于震动与高压断路器安装位置方向、海拔、气压和固定基座的材质均相关,而且在高压断路器震动的频率和幅度的变化也与震动传感器位置也相关,所以要使得实验室获得的方法生效需要严格要求设备的装配位置和过程;当影响震动的某一个条件发生改变时,传统的智能模型需要构建实验数据集对其进行适应;而实验室模拟能力有限,不可能穷举所有的安装时出现的可能性,所以传统模式在实际使用过程中往往会失效,出现大量的故障误报、漏报的情况。(2)传统方法需要在实验室构造高压断路器基座螺丝松动和螺丝不松动样本,由于螺丝松动样本不可能穷举所有松动情况,所以这种情况下利用神经网、决策树、支持向量机等浅层人工智能模型可能会导致模型与特定松动样本绑定,在实际应用中这种绑定导致预测能力较低;一些深度神经网模型增加了高层特征抽取功能来防止绑定,但是这类神经网需要较为强大的计算机或服务器支持;而高压断路器所处位置通常难于安放这些服务器,高频的震动信息上传到服务器需要额外的网络和存储设备,这些设备增加了电网的总体维护和建设成本。
发明内容
针对现有技术存在的问题,本发明提供一种高压断路器基座螺丝松动故障检测方法。该方法同时引入两个震动传感器,构建了加权分段震动差异算子来衡量两个传感器之间震动差异,进而衡量震动差异数据与高压断路器正常运行状态下的距离来判断高压断路器基座螺丝是否存在松动故障。
本发明所述的一种高压断路器基座螺丝松动故障检测方法,包括以下步骤:
S1.在高压断路器中安装震动传感器A1和震动传感器B1;利用A1和B1收集高压断路器正常运行状态并执行打开动作时的震动数据,形成震动数据集VSet;输入参考窗口变量VWindowsSize;获取数据条目数变量NumRows;
S101,在高压断路器中安装震动检测器A1和震动检测器B1, A1可以安装在高压断路器的任意位置,B1可以安装在高压断路器的任意位置但需要保证距离A1直线距离在15厘米以上;震动检测器的采集频率为800Hz;
S102,利用A1和B1收集高压断路器正常运行状态并执行打开动作时的数据,形成震动数据集VSet;VSet为包含二个字段的数据表,字段包括:
FieldA1:震动传感器A1获取的数据,为一个800个元素的数组;
FieldB1:震动传感器B1获取的数据,为一个800个元素的数组;
S103,输入参考窗口变量VWindowsSize, VWindowsSize的默认值为5;
S104,数据条目数变量NumRows=VSet的数据条目数;
S105,本步骤处理过程结束;
S2. 建立分段振动值抽取算子FDExc, 该算子的输入包括分段振动数组FDExcInputArray,计算起始位置FDExcStart,计算终止位置FDExcEnd;分段振动值抽取算子输出变量为FDExcOutput;
S201,建立分段振动值抽取算子FDExc, 该算子的输入包括分段振动数组FDExcInputArray,计算起始位置FDExcStart,计算终止位置FDExcEnd;
S202,分段振动值抽取算子暂存数组FDExcArray=取出FDExcInputArray数组中以FDExcStart位置为起点并以FDExcEnd为终点的所有元素;
S203, 分段振动值抽取算子最大位置FDExcMaxPos=FDExcArray中值最大值元素在FDExcArray中所处位置;
S204,分段振动值抽取算子输出变量FDExcOutput=建立5个元素的数组,元素初始值为0;
S205,分段振动值抽取算子计数器FDExcCounter=1;
S206,计算分段振动值抽取算子暂存起始位置值:
FDExcTempStart=FDExcMaxPos+(FDExcCounter-1)×VWindowsSize;
S207,计算分段振动值抽取算子暂存终止位置值:
FDExcTempEnd=FDExcTempStart+VWindowsSize;
S208,分段振动值抽取算子循环暂存数组FDExcItTemp=取出FDExcArray中以FDExcTempStart位置为起始并以FDExcTempEnd为终止的所有元素; S209,FDExcItTemp=ABS(FDExcItTemp), 其中ABS为绝对值函数;
S210,FDExcOutput[FDExcCounter]=计算FDExcItTemp的均值;
S211,FDExcCounter=FDExcCounter+1;
S212,如果FDExcCounter大于5则转到S213,否则转到S206;
S213,将FDExcOutput作为算子计算机结果输出;
S3, 建立加权分段震动差异算子FDDiff, FDDiff的输入为加权分段震动差异算子第一输入数组FDDiffArray1和加权分段震动差异算子第二输入数组FDDiffArray2,输出为加权分段震动差异结果变量FDDiffResult;
S301,建立加权分段震动差异算子FDDiff, FDDiff的输入为加权分段震动差异算子第一输入数组FDDiffArray1和加权分段震动差异算子第二输入数组FDDiffArray2;
S302,建立加权分段震动差异结果变量FDDiffResult=建立15个元素的数组,元素初始值为0;
S303,建立分段震动差异第一暂存变量FDDiffTemp1=0;建立分段震动差异第二暂存变量FDDiffTemp2=0;
S304,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=1,FDExcEnd=150;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S305,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=1,FDExcEnd=150;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S306,设置FDDiffResult第1至5元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×2.5, 其中TANH为双曲正切函数;
S307,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=151,FDExcEnd=570;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S308,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=151,FDExcEnd=570;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S309,设置FDDiffResult第6至10元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×1.2, 其中TANH为双曲正切函数;
S310,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=571,FDExcEnd=800;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S311,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=571,FDExcEnd=800;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S312,设置FDDiffResult第11至15元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×0.5, 其中TANH为双曲正切函数;
S313,将FDDiffResult作为算子计算结果输出;
S4,基于加权分段震动差异算子FDDiff对VSet进行处理,获取震动差异中心矢量CenterVector,获取震动差异最大值变量MaxD;
S401,建立震动差异暂存列表DiffListTemp=空的列表;
S402,处理计数器ProcessCounter=1;
S403, 建立当前差异暂存矢量DiffTempVector=利用FDDiff进行计算,算子输入FDDiffArray1=VSet[ProcessCounter].FieldA1;
FDDiffArray2=VSet[ProcessCounter].FieldA2; FDDiff的计算结果FDDiffResult返回给DiffTempVector;
S404,将DiffTempVector加入到DiffListTemp之中;
S405,ProcessCounter=ProcessCounter+1;
S406,如果ProcessCounter大于NumRows则转到S407,否则转到S403;
S407,建立震动差异中心矢量CenterVector=DiffListTemp的所有矢量的均值;
S408,ProcessCounter=1;建立震动差异最大值变量MaxD=0;
S409,处理暂存距离值TempD=norm1(DiffListTemp[ProcessCounter]-CenterVector), 其中norm1为计算矢量的1-范数;
S410,如果TempD>MaxD则MaxD=TempD;
S411,ProcessCounter=ProcessCounter+1;
S412,如果ProcessCounter大于NumRows则转到S413,否则转向S409;
S413,本步骤处理过程结束;
S5,在高压断路器执行打开动作时,利用A1获取震动数据PredArray1,利用B1获取震动数据PredArray2;利用加权分段震动差异算子FDDiff,震动差异中心矢量CenterVector和震动差异最大值变量MaxD对是否出现基座螺丝松动故障进行判断;
S501,在高压断路器执行打开动作时,利用A1获取震动数据PredArray1,利用B1获取震动数据PredArray2;
S502,判定暂存矢量PredTempVector=利用FDDiff进行计算,算子输入FDDiffArray1= PredArray1, FDDiffArray2= PredArray2; FDDiff的计算结果FDDiffResult返回给PredTempVector;
S503,判断暂存距离PredTempD=norm1(PredTempVector-CenterVector), 其中norm1为计算矢量的1-范数;
S504,如果PredTempD大于MaxD则转到S505,否则转到S506;
S505,输出存在基座螺丝松动故障,转到S507;
S506,输出不存在基座螺丝松动故障;
S507,本步骤处理过程结束。
本发明的有益效果是:
针对现有技术存在的问题,本发明提供一种高压断路器基座螺丝松动故障检测方法。该方法同时引入两个震动传感器,构建了加权分段震动差异算子来衡量两个传感器之间震动差异,进而衡量震动差异数据与高压断路器正常运行状态下的距离来判断高压断路器基座螺丝是否存在松动故障。利用本发明对于两个震动传感器的安装位置和高压断路器连接基座的材质没有严格的要求,降低设备装配条件限制;同时在数据上本发明仅需要收集震动传感器正常运行时的数据,不会与具体的基座螺丝松动样本绑定,提高检测成功率,保障高压断路器安全稳定运行。
具体实施方式
通过以下实施例进一步举例描述本发明,并不以任何方式限制本发明,在不背离本发明的技术解决方案的前提下,对本发明所作的本领域普通技术人员容易实现的任何改动或改变都将落入本发明的权利要求范围之内。
实施例1
S1.在高压断路器中安装震动传感器A1和震动传感器B1;利用A1和B1收集高压断路器正常运行状态并执行打开动作时的震动数据,形成震动数据集VSet;输入参考窗口变量VWindowsSize;获取数据条目数变量NumRows;
S101,在高压断路器中安装震动检测器A1和震动检测器B1, A1可以安装在高压断路器的任意位置,B1可以安装在高压断路器的任意位置但需要保证距离A1直线距离在15厘米以上;震动检测器的采集频率为800Hz;
S102,利用A1和B1收集高压断路器正常运行状态并执行打开动作时的数据,形成震动数据集VSet;VSet为包含二个字段的数据表,字段包括:
FieldA1:震动传感器A1获取的数据,为一个800个元素的数组;
FieldB1:震动传感器B1获取的数据,为一个800个元素的数组;
S103,输入参考窗口变量VWindowsSize, VWindowsSize的默认值为5;
S104,数据条目数变量NumRows=VSet的数据条目数;
S105,本步骤处理过程结束;
S2. 建立分段振动值抽取算子FDExc, 该算子的输入包括分段振动数组FDExcInputArray,计算起始位置FDExcStart,计算终止位置FDExcEnd;分段振动值抽取算子输出变量为FDExcOutput;
S201,建立分段振动值抽取算子FDExc, 该算子的输入包括分段振动数组FDExcInputArray,计算起始位置FDExcStart,计算终止位置FDExcEnd;
S202,分段振动值抽取算子暂存数组FDExcArray=取出FDExcInputArray数组中以FDExcStart位置为起点并以FDExcEnd为终点的所有元素;
S203, 分段振动值抽取算子最大位置FDExcMaxPos=FDExcArray中值最大值元素在FDExcArray中所处位置;
S204,分段振动值抽取算子输出变量FDExcOutput=建立5个元素的数组,元素初始值为0;
S205,分段振动值抽取算子计数器FDExcCounter=1;
S206,计算分段振动值抽取算子暂存起始位置值:
FDExcTempStart=FDExcMaxPos+(FDExcCounter-1)×VWindowsSize;
S207,计算分段振动值抽取算子暂存终止位置值:
FDExcTempEnd=FDExcTempStart+VWindowsSize;
S208,分段振动值抽取算子循环暂存数组FDExcItTemp=取出FDExcArray中以FDExcTempStart位置为起始并以FDExcTempEnd为终止的所有元素; S209,FDExcItTemp=ABS(FDExcItTemp), 其中ABS为绝对值函数;
S210,FDExcOutput[FDExcCounter]=计算FDExcItTemp的均值;
S211,FDExcCounter=FDExcCounter+1;
S212,如果FDExcCounter大于5则转到S213,否则转到S206;
S213,将FDExcOutput作为算子计算机结果输出;
S3, 建立加权分段震动差异算子FDDiff, FDDiff的输入为加权分段震动差异算子第一输入数组FDDiffArray1和加权分段震动差异算子第二输入数组FDDiffArray2,输出为加权分段震动差异结果变量FDDiffResult;
S301,建立加权分段震动差异算子FDDiff, FDDiff的输入为加权分段震动差异算子第一输入数组FDDiffArray1和加权分段震动差异算子第二输入数组FDDiffArray2;
S302,建立加权分段震动差异结果变量FDDiffResult=建立15个元素的数组,元素初始值为0;
S303,建立分段震动差异第一暂存变量FDDiffTemp1=0;建立分段震动差异第二暂存变量FDDiffTemp2=0;
S304,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=1,FDExcEnd=150;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S305,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=1,FDExcEnd=150;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S306,设置FDDiffResult第1至5元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×2.5, 其中TANH为双曲正切函数;
S307,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=151,FDExcEnd=570;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S308,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=151,FDExcEnd=570;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S309,设置FDDiffResult第6至10元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×1.2, 其中TANH为双曲正切函数;
S310,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=571,FDExcEnd=800;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S311,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=571,FDExcEnd=800;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S312,设置FDDiffResult第11至15元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×0.5, 其中TANH为双曲正切函数;
S313,将FDDiffResult作为算子计算结果输出;
S4,基于加权分段震动差异算子FDDiff对VSet进行处理,获取震动差异中心矢量CenterVector,获取震动差异最大值变量MaxD;
S401,建立震动差异暂存列表DiffListTemp=空的列表;
S402,处理计数器ProcessCounter=1;
S403, 建立当前差异暂存矢量DiffTempVector=利用FDDiff进行计算,算子输入FDDiffArray1=VSet[ProcessCounter].FieldA1;
FDDiffArray2=VSet[ProcessCounter].FieldA2; FDDiff的计算结果FDDiffResult返回给DiffTempVector;
S404,将DiffTempVector加入到DiffListTemp之中;
S405,ProcessCounter=ProcessCounter+1;
S406,如果ProcessCounter大于NumRows则转到S407,否则转到S403;
S407,建立震动差异中心矢量CenterVector=DiffListTemp的所有矢量的均值;
S408,ProcessCounter=1;建立震动差异最大值变量MaxD=0;
S409,处理暂存距离值TempD=norm1(DiffListTemp[ProcessCounter]-CenterVector), 其中norm1为计算矢量的1-范数;
S410,如果TempD>MaxD则MaxD=TempD;
S411,ProcessCounter=ProcessCounter+1;
S412,如果ProcessCounter大于NumRows则转到S413,否则转向S409;
S413,本步骤处理过程结束;
S5,在高压断路器执行打开动作时,利用A1获取震动数据PredArray1,利用B1获取震动数据PredArray2;利用加权分段震动差异算子FDDiff,震动差异中心矢量CenterVector和震动差异最大值变量MaxD对是否出现基座螺丝松动故障进行判断;
S501,在高压断路器执行打开动作时,利用A1获取震动数据PredArray1,利用B1获取震动数据PredArray2;
S502,判定暂存矢量PredTempVector=利用FDDiff进行计算,算子输入FDDiffArray1= PredArray1, FDDiffArray2= PredArray2; FDDiff的计算结果FDDiffResult返回给PredTempVector;
S503,判断暂存距离PredTempD=norm1(PredTempVector-CenterVector), 其中norm1为计算矢量的1-范数;
S504,如果PredTempD大于MaxD则转到S505,否则转到S506;
S505,输出存在基座螺丝松动故障,转到S507;
S506,输出不存在基座螺丝松动故障;
S507,本步骤处理过程结束。
实施例2
以北方地区某地安装的高压断路器为实验对象实施本发明技术方案:
S1.在高压断路器中安装震动传感器A1和震动传感器B1;利用A1和B1收集高压断路器正常运行状态并执行打开动作时的震动数据,形成震动数据集VSet;输入参考窗口变量VWindowsSize;获取数据条目数变量NumRows;
S101,在高压断路器中安装震动检测器A1和震动检测器B1, A1可以安装在高压断路器的任意位置,B1可以安装在高压断路器的任意位置但需要保证距离A1直线距离在15厘米以上;震动检测器的采集频率为800Hz;
S102,利用A1和B1收集高压断路器正常运行状态并执行打开动作时的数据,形成震动数据集VSet;VSet为包含二个字段的数据表,字段包括:
FieldA1:震动传感器A1获取的数据,为一个800个元素的数组;
FieldB1:震动传感器B1获取的数据,为一个800个元素的数组;
S103,输入参考窗口变量VWindowsSize, VWindowsSize的默认值为5;
S104,数据条目数变量NumRows=VSet的数据条目数;
S105,本步骤处理过程结束。
对应的震动数据集VSet为一个包含30条记录数据集:
输入参考窗口变量VWindowsSize=5;
数据条目数变量NumRows=30;
S2. 建立分段振动值抽取算子FDExc, 该算子的输入包括分段振动数组FDExcInputArray,计算起始位置FDExcStart,计算终止位置FDExcEnd;分段振动值抽取算子输出变量为FDExcOutput:
S201,建立分段振动值抽取算子FDExc, 该算子的输入包括分段振动数组FDExcInputArray,计算起始位置FDExcStart,计算终止位置FDExcEnd;
S202,分段振动值抽取算子暂存数组FDExcArray=取出FDExcInputArray数组中以FDExcStart位置为起点并以FDExcEnd为终点的所有元素;
S203, 分段振动值抽取算子最大位置FDExcMaxPos=FDExcArray中值最大值元素在FDExcArray中所处位置;
S204,分段振动值抽取算子输出变量FDExcOutput=建立5个元素的数组,元素初始值为0;
S205,分段振动值抽取算子计数器FDExcCounter=1;
S206,计算分段振动值抽取算子暂存起始位置值:
FDExcTempStart=FDExcMaxPos+(FDExcCounter-1)×VWindowsSize;
S207,计算分段振动值抽取算子暂存终止位置值:
FDExcTempEnd=FDExcTempStart+VWindowsSize;
S208,分段振动值抽取算子循环暂存数组FDExcItTemp=取出FDExcArray中以FDExcTempStart位置为起始并以FDExcTempEnd为终止的所有元素;
S209,FDExcItTemp=ABS(FDExcItTemp), 其中ABS为绝对值函数;
S210,FDExcOutput[FDExcCounter]=计算FDExcItTemp的均值;
S211,FDExcCounter=FDExcCounter+1;
S212,如果FDExcCounter大于5则转到S213,否则转到S206;
S213,将FDExcOutput作为算子计算机结果输出;
S3, 建立加权分段震动差异算子FDDiff, FDDiff的输入为加权分段震动差异算子第一输入数组FDDiffArray1和加权分段震动差异算子第二输入数组FDDiffArray2,输出为加权分段震动差异结果变量FDDiffResult;
S301,建立加权分段震动差异算子FDDiff, FDDiff的输入为加权分段震动差异算子第一输入数组FDDiffArray1和加权分段震动差异算子第二输入数组FDDiffArray2;
S302,建立加权分段震动差异结果变量FDDiffResult=建立15个元素的数组,元素初始值为0;
S303,建立分段震动差异第一暂存变量FDDiffTemp1=0;建立分段震动差异第二暂存变量FDDiffTemp2=0;
S304,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=1,FDExcEnd=150;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S305,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=1,FDExcEnd=150;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S306,设置FDDiffResult第1至5元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×2.5, 其中TANH为双曲正切函数;
S307,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=151,FDExcEnd=570;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S308,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=151,FDExcEnd=570;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S309,设置FDDiffResult第6至10元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×1.2, 其中TANH为双曲正切函数;
S310,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=571,FDExcEnd=800;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S311,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=571,FDExcEnd=800;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S312,设置FDDiffResult第11至15元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×0.5, 其中TANH为双曲正切函数;
S313,将FDDiffResult作为算子计算结果输出。
S4,基于加权分段震动差异算子FDDiff对VSet进行处理,获取震动差异中心矢量CenterVector,获取震动差异最大值变量MaxD;
S401,建立震动差异暂存列表DiffListTemp=空的列表:
S402,处理计数器ProcessCounter=1;
S403, 建立当前差异暂存矢量DiffTempVector=利用FDDiff进行计算,算子输入FDDiffArray1=VSet[ProcessCounter].FieldA1, FDDiffArray2=VSet[ProcessCounter].FieldA2; FDDiff的计算结果FDDiffResult返回给DiffTempVector;
S404,将DiffTempVector加入到DiffListTemp之中;
S405,ProcessCounter=ProcessCounter+1;
S406,如果ProcessCounter大于NumRows则转到S407,否则转到S403;
S407,建立震动差异中心矢量CenterVector=DiffListTemp的所有矢量的均值;
S408,ProcessCounter=1;建立震动差异最大值变量MaxD=0;
S409,处理暂存距离值TempD=norm1(DiffListTemp[ProcessCounter]-CenterVector), 其中norm1为计算矢量的1-范数;
S410,如果TempD>MaxD则MaxD=TempD;
S411,ProcessCounter=ProcessCounter+1;
S412,如果ProcessCounter大于NumRows则转到S413,否则转向S409;
S413,本步骤处理过程结束。
获取震动差异中心矢量CenterVector=[1.58, 1.59, 1.57, 1.57, 1.57, 0.94,0.85, 0.92, 0.91, 0.94, 0.75, 0.72, 0.77, 0.75, 0.61];
获取震动差异最大值变量MaxD=6.41;
S5,在高压断路器执行打开动作时,利用A1获取震动数据PredArray1,利用B1获取震动数据PredArray2;利用加权分段震动差异算子FDDiff,震动差异中心矢量CenterVector和震动差异最大值变量MaxD对是否出现基座螺丝松动故障进行判断;
S501,在高压断路器执行打开动作时,利用A1获取震动数据PredArray1,利用B1获取震动数据PredArray2;
S502,判定暂存矢量PredTempVector=利用FDDiff进行计算,算子输入FDDiffArray1= PredArray1, FDDiffArray2= PredArray2; FDDiff的计算结果FDDiffResult返回给PredTempVector;
S503,判断暂存距离PredTempD=norm1(PredTempVector-CenterVector), 其中norm1为计算矢量的1-范数;
S504,如果PredTempD大于MaxD则转到S505,否则转到S506;
S505,输出存在基座螺丝松动故障,转到S507;
S506,输出不存在基座螺丝松动故障;
S507,本步骤处理过程结束。
当该高压断路器正常运行时,利用A1获取震动数据PredArray1,利用B1获取震动数据PredArray2;利用加权分段震动差异算子FDDiff和震动差异中心矢量CenterVector计算获得PredTempD=3.22; 由于PredTempD小于MaxD,所以发明方法输出不存在基座螺丝松动故障,本发明方法做出了正确判断。
在故意松动高压断路的一个基座螺丝之后,利用A1获取震动数据PredArray1,利用B1获取震动数据PredArray2;利用加权分段震动差异算子FDDiff和震动差异中心矢量CenterVector计算获得PredTempD=7.22; 由于PredTempD大于MaxD,所以发明方法输出存在基座螺丝松动故障,本发明方法正确检测到了基座螺丝松动故障。
Claims (1)
1.一种高压断路器基座螺丝松动故障检测方法,包括以下步骤:
S1.在高压断路器中安装震动传感器A1和震动传感器B1;利用A1和B1收集高压断路器正常运行状态并执行打开动作时的震动数据,形成震动数据集VSet;输入参考窗口变量VWindowsSize;获取数据条目数变量NumRows;
S101,在高压断路器中安装震动检测器A1和震动检测器B1, A1可以安装在高压断路器的任意位置,B1可以安装在高压断路器的任意位置但需要保证距离A1直线距离在15厘米以上;震动检测器的采集频率为800Hz;
S102,利用A1和B1收集高压断路器正常运行状态并执行打开动作时的数据,形成震动数据集VSet;VSet为包含二个字段的数据表,字段包括:
FieldA1:震动传感器A1获取的数据,为一个800个元素的数组;
FieldB1:震动传感器B1获取的数据,为一个800个元素的数组;
S103,输入参考窗口变量VWindowsSize, VWindowsSize的默认值为5;
S104,数据条目数变量NumRows=VSet的数据条目数;
S105,本步骤处理过程结束;
S2. 建立分段振动值抽取算子FDExc, 该算子的输入包括分段振动数组FDExcInputArray,计算起始位置FDExcStart,计算终止位置FDExcEnd;分段振动值抽取算子输出变量为FDExcOutput;
S201,建立分段振动值抽取算子FDExc, 该算子的输入包括分段振动数组FDExcInputArray,计算起始位置FDExcStart,计算终止位置FDExcEnd;
S202,分段振动值抽取算子暂存数组FDExcArray=取出FDExcInputArray数组中以FDExcStart位置为起点并以FDExcEnd为终点的所有元素;
S203, 分段振动值抽取算子最大位置FDExcMaxPos=FDExcArray中值最大值元素在FDExcArray中所处位置;
S204,分段振动值抽取算子输出变量FDExcOutput=建立5个元素的数组,元素初始值为0;
S205,分段振动值抽取算子计数器FDExcCounter=1;
S206,计算分段振动值抽取算子暂存起始位置值:
FDExcTempStart=FDExcMaxPos+(FDExcCounter-1)×VWindowsSize;
S207,计算分段振动值抽取算子暂存终止位置值:
FDExcTempEnd=FDExcTempStart+VWindowsSize;
S208,分段振动值抽取算子循环暂存数组FDExcItTemp=取出FDExcArray中以FDExcTempStart位置为起始并以FDExcTempEnd为终止的所有元素; S209,FDExcItTemp=ABS(FDExcItTemp), 其中ABS为绝对值函数;
S210,FDExcOutput[FDExcCounter]=计算FDExcItTemp的均值;
S211,FDExcCounter=FDExcCounter+1;
S212,如果FDExcCounter大于5则转到S213,否则转到S206;
S213,将FDExcOutput作为算子计算机结果输出;
S3, 建立加权分段震动差异算子FDDiff, FDDiff的输入为加权分段震动差异算子第一输入数组FDDiffArray1和加权分段震动差异算子第二输入数组FDDiffArray2,输出为加权分段震动差异结果变量FDDiffResult;
S301,建立加权分段震动差异算子FDDiff, FDDiff的输入为加权分段震动差异算子第一输入数组FDDiffArray1和加权分段震动差异算子第二输入数组FDDiffArray2;
S302,建立加权分段震动差异结果变量FDDiffResult=建立15个元素的数组,元素初始值为0;
S303,建立分段震动差异第一暂存变量FDDiffTemp1=0;建立分段震动差异第二暂存变量FDDiffTemp2=0;
S304,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=1,FDExcEnd=150;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S305,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=1,FDExcEnd=150;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S306,设置FDDiffResult第1至5元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×2.5,其中TANH为双曲正切函数;
S307,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=151,FDExcEnd=570;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S308,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=151,FDExcEnd=570;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S309,设置FDDiffResult第6至10元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×1.2,其中TANH为双曲正切函数;
S310,计算FDDiffTemp1=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray1,FDExcStart=571,FDExcEnd=800;FDExc的计算结果FDExcOutput返回给FDDiffTemp1;
S311,计算FDDiffTemp2=利用FDExc进行计算,FDExc算子的FDExcInputArray=FDDiffArray2,FDExcStart=571,FDExcEnd=800;FDExc的计算结果FDExcOutput返回给FDDiffTemp2;
S312,设置FDDiffResult第11至15元素的值=TANH(FDDiffTemp1-FDDiffTemp2)×0.5,其中TANH为双曲正切函数;
S313,将FDDiffResult作为算子计算结果输出;
S4,基于加权分段震动差异算子FDDiff对VSet进行处理,获取震动差异中心矢量CenterVector,获取震动差异最大值变量MaxD;
S401,建立震动差异暂存列表DiffListTemp=空的列表;
S402,处理计数器ProcessCounter=1;
S403, 建立当前差异暂存矢量DiffTempVector=利用FDDiff进行计算,算子输入FDDiffArray1=VSet[ProcessCounter].FieldA1;
FDDiffArray2=VSet[ProcessCounter].FieldA2; FDDiff的计算结果FDDiffResult返回给DiffTempVector;
S404,将DiffTempVector加入到DiffListTemp之中;
S405,ProcessCounter=ProcessCounter+1;
S406,如果ProcessCounter大于NumRows则转到S407,否则转到S403;
S407,建立震动差异中心矢量CenterVector=DiffListTemp的所有矢量的均值;
S408,ProcessCounter=1;建立震动差异最大值变量MaxD=0;
S409,处理暂存距离值TempD=norm1(DiffListTemp[ProcessCounter]-CenterVector), 其中norm1为计算矢量的1-范数;
S410,如果TempD>MaxD则MaxD=TempD;
S411,ProcessCounter=ProcessCounter+1;
S412,如果ProcessCounter大于NumRows则转到S413,否则转向S409;
S413,本步骤处理过程结束;
S5,在高压断路器执行打开动作时,利用A1获取震动数据PredArray1,利用B1获取震动数据PredArray2;利用加权分段震动差异算子FDDiff,震动差异中心矢量CenterVector和震动差异最大值变量MaxD对是否出现基座螺丝松动故障进行判断;
S501,在高压断路器执行打开动作时,利用A1获取震动数据PredArray1,利用B1获取震动数据PredArray2;
S502,判定暂存矢量PredTempVector=利用FDDiff进行计算,算子输入FDDiffArray1=PredArray1, FDDiffArray2= PredArray2; FDDiff的计算结果FDDiffResult返回给PredTempVector;
S503,判断暂存距离PredTempD=norm1(PredTempVector-CenterVector), 其中norm1为计算矢量的1-范数;
S504,如果PredTempD大于MaxD则转到S505,否则转到S506;
S505,输出存在基座螺丝松动故障,转到S507;
S506,输出不存在基座螺丝松动故障;
S507,本步骤处理过程结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011316591.5A CN112504435B (zh) | 2020-11-23 | 2020-11-23 | 一种高压断路器基座螺丝松动故障检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011316591.5A CN112504435B (zh) | 2020-11-23 | 2020-11-23 | 一种高压断路器基座螺丝松动故障检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112504435A true CN112504435A (zh) | 2021-03-16 |
CN112504435B CN112504435B (zh) | 2022-09-20 |
Family
ID=74958792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011316591.5A Active CN112504435B (zh) | 2020-11-23 | 2020-11-23 | 一种高压断路器基座螺丝松动故障检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112504435B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011103230A (ja) * | 2009-11-11 | 2011-05-26 | Sumitomo Metal Ind Ltd | 開閉器の診断方法及び開閉器の診断装置 |
CN103487749A (zh) * | 2013-09-18 | 2014-01-01 | 国家电网公司 | 高压断路器机械状态在线监测诊断系统及方法 |
CN109100646A (zh) * | 2018-08-17 | 2018-12-28 | 国网江苏省电力有限公司检修分公司 | 一种高压断路器故障诊断方法 |
CN109357846A (zh) * | 2018-08-06 | 2019-02-19 | 长春工程学院 | 一种高压断路器电磁铁卡涩故障检测方法 |
CN110796120A (zh) * | 2019-11-21 | 2020-02-14 | 杭州电力设备制造有限公司 | 基于时域特征的断路器机械故障XGBoost诊断方法 |
CN111175698A (zh) * | 2020-01-18 | 2020-05-19 | 国网山东省电力公司菏泽供电公司 | 一种基于声振结合的变压器噪声源定位方法、系统和装置 |
WO2020107828A1 (en) * | 2018-11-28 | 2020-06-04 | Zhejiang Future Technology Institute (jiaxing) | Method and system for detecting fastening state of fastening structure |
CN111473860A (zh) * | 2020-05-29 | 2020-07-31 | 国网安徽省电力有限公司电力科学研究院 | 一种用于高压交流断路器振动信号特征参数提取方法 |
-
2020
- 2020-11-23 CN CN202011316591.5A patent/CN112504435B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011103230A (ja) * | 2009-11-11 | 2011-05-26 | Sumitomo Metal Ind Ltd | 開閉器の診断方法及び開閉器の診断装置 |
CN103487749A (zh) * | 2013-09-18 | 2014-01-01 | 国家电网公司 | 高压断路器机械状态在线监测诊断系统及方法 |
CN109357846A (zh) * | 2018-08-06 | 2019-02-19 | 长春工程学院 | 一种高压断路器电磁铁卡涩故障检测方法 |
CN109100646A (zh) * | 2018-08-17 | 2018-12-28 | 国网江苏省电力有限公司检修分公司 | 一种高压断路器故障诊断方法 |
WO2020107828A1 (en) * | 2018-11-28 | 2020-06-04 | Zhejiang Future Technology Institute (jiaxing) | Method and system for detecting fastening state of fastening structure |
CN110796120A (zh) * | 2019-11-21 | 2020-02-14 | 杭州电力设备制造有限公司 | 基于时域特征的断路器机械故障XGBoost诊断方法 |
CN111175698A (zh) * | 2020-01-18 | 2020-05-19 | 国网山东省电力公司菏泽供电公司 | 一种基于声振结合的变压器噪声源定位方法、系统和装置 |
CN111473860A (zh) * | 2020-05-29 | 2020-07-31 | 国网安徽省电力有限公司电力科学研究院 | 一种用于高压交流断路器振动信号特征参数提取方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112504435B (zh) | 2022-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110933031A (zh) | 一种基于lstm的智能电网配电终端单元入侵检测方法 | |
CN105703258B (zh) | Gis开关设备动作状态监测系统及其使用方法 | |
CN107219453B (zh) | 一种基于多维混合量的变电站继电保护隐性故障诊断方法 | |
CN117590223B (zh) | 断路器的在线监测系统及其方法 | |
CN108418304B (zh) | 变电站二次回路状态监测方法、装置及系统 | |
CN117078017A (zh) | 一种用于电力电网设备监控智能化决策分析系统 | |
CN115516724A (zh) | 用于智能电气保护系统的数据分析 | |
CN111551352B (zh) | 一种gis设备的断路器的状态的检测方法及系统 | |
CN113379252A (zh) | 基于多尺度排列熵值的贝叶斯网络高压开关柜健康系统评估方法 | |
CN109490715A (zh) | 一种极端环境下的电力系统故障判别方法 | |
WO2022105285A1 (zh) | 一种基于冲击事件来监测异常状态的系统及方法 | |
CN112504435B (zh) | 一种高压断路器基座螺丝松动故障检测方法 | |
CN113161164B (zh) | 一种用于真空断路器触头磨损的测量方法 | |
CN107147143B (zh) | 一种风机连锁脱网故障早期预警模型建立方法 | |
Li et al. | Real-time detecting false data injection attacks based on spatial and temporal correlations | |
CN114967571B (zh) | 基于互联网的充气柜远程监控系统 | |
CN115632379A (zh) | 一种多途径检测的vcs隔离开关箱用电弧光保护系统 | |
CN113125949B (zh) | 一种基于电场的高压断路器绝缘监测方法及系统 | |
US11933834B2 (en) | Method and system for detecting self-clearing, sub-cycle faults | |
CN109357846B (zh) | 一种高压断路器电磁铁卡涩故障检测方法 | |
CN107329018A (zh) | 一种变电站继电保护系统测量环节隐性故障检测方法 | |
Hou et al. | Discriminant approach to the GIS mechanical fault diagnosis based on the KL divergence of vibration signals | |
CN117851910B (zh) | 输电线路新型智能接地装置挂接与拆除自感知方法及系统 | |
CN207424178U (zh) | 监测系统 | |
CN113484573B (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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20220831 Address after: No. 395, KUANPING Road, Chaoyang District, Changchun City, Jilin Province Applicant after: CHANGCHUN INSTITUTE OF TECHNOLOGY Applicant after: LIAOYUAN POWER SUPPLY COMPANY OF STATE GRID JILINSHENG ELECTRIC POWER SUPPLY Co. Address before: No. 395, KUANPING Road, Chaoyang District, Changchun City, Jilin Province Applicant before: CHANGCHUN INSTITUTE OF TECHNOLOGY |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |