CN110081968B - 一种变压器振动信号的模拟方法和装置 - Google Patents

一种变压器振动信号的模拟方法和装置 Download PDF

Info

Publication number
CN110081968B
CN110081968B CN201910467079.1A CN201910467079A CN110081968B CN 110081968 B CN110081968 B CN 110081968B CN 201910467079 A CN201910467079 A CN 201910467079A CN 110081968 B CN110081968 B CN 110081968B
Authority
CN
China
Prior art keywords
vibration signal
signal sequence
transformer
calculating
absolute value
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.)
Expired - Fee Related
Application number
CN201910467079.1A
Other languages
English (en)
Other versions
CN110081968A (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.)
Guangdong University of Petrochemical Technology
Original Assignee
Guangdong University of Petrochemical Technology
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 Guangdong University of Petrochemical Technology filed Critical Guangdong University of Petrochemical Technology
Priority to CN201910467079.1A priority Critical patent/CN110081968B/zh
Publication of CN110081968A publication Critical patent/CN110081968A/zh
Application granted granted Critical
Publication of CN110081968B publication Critical patent/CN110081968B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Protection Of Transformers (AREA)
  • Testing Relating To Insulation (AREA)

Abstract

本发明的实施例公开一种变压器振动信号的模拟方法和装置,所述方法包括:步骤1,输入实测的振动信号序列P=[P1,P2,…,PN],N为振动信号序列的长度;骤2,计算所述振动信号序列的幅度参数A0:步骤3,计算所述振动信号序列的主频f0;步骤4,计算所述振动信号序列的衰减系数η0;步骤5,根据所述幅度参数、主频、衰减系数,建立所述振动信号的模型。

Description

一种变压器振动信号的模拟方法和装置
技术领域
本发明智能电网领域,尤其涉及一种变压器振动信号的模拟方法和装置。
背景技术
随着智能电网的高速发展,电力设备安全稳定运行显得尤其重要。目前,对超高压及以上电压等级的电力设备开展运行状态检测,尤其是对异常状态的检测显得愈加重要和迫切。电力变压器作为电力系统的重要组成部分,是变电站中最重要的电气设备之一,其可靠运行关系到电网的安全。一般而言,变压器的异常状态可分为铁芯异常与绕组异常。铁芯异常主要表现为铁芯饱和,绕组异常通常包括绕组变形、绕组松动等。
变压器异常状态检测的基本原理是提取变压器运行中的各特征量,分析、辨识并跟踪特征量以此监测变压器的异常运行状态。检测方法按照接触程度可分为侵入式检测和非侵入式检测;按照是否需停机检测可分为带电检测和停电检测;按照检测量类型可以分为电气量法和非电气量法等。相比而言,非侵入式检测可移植性强,安装更方便;带电检测不影响变压器运行;非电气量法与电力系统无电气连接,更为安全。当前变压器运行状态的常用检测方法中,包括检测局部放电的脉冲电流法和超声波检测法、检测绕组变形的频率响应法以及检测机械及电气故障的振动检测法等。这些检测方法主要检测变压器绝缘状况及机械结构状况,其中以变压器振动信号(振声)的检测最为全面,对于大部分变压器故障及异常状态均能有所反应。
变压器在运行过程中,铁芯硅钢片的磁致伸缩与绕组电动力引起的振动会向四周辐射不同幅值和频率的振声信号。变压器正常运行时对外发出的是均匀的低频噪声;如果发出不均匀声音,则属不正常现象。变压器在不同运行状态下会发出有区别性的声音,可通过对其发出声音的检测,掌握变压器的运行状况。值得关注的是,对变压器不同运行状态下发出声音的检测不仅可以检测很多种引起电气量变化的严重故障,还可以检测许多并未危及绝缘的没有引起电气量变化的异常状态,比如变压器内外部零部件松动等。
综上所述,变压器运行状态监测的常用方法是振声检测,此方法利用了变压器发出的振动信号,振动信号很容易受到环境噪声的影响,环境噪声对此检测方法的性能有较大的影响,因此在不同信噪比和不同噪声类型下评估振声检测算法性能非常重要,这就需要对变压器振动信号进行有效模拟,但现在缺少变压器振动信号的模拟方法。
发明内容
有鉴于此,本发明实施例提供一种变压器振动信号的模拟方法和装置,能够有效地模拟变压器振动信号。
一种变压器振动信号的模拟方法,包括:
步骤1,输入实测的振动信号序列P=[P1,P2,…,PN],N为振动信号序列的长度;
步骤2,计算所述振动信号序列的幅度参数A0
步骤3,计算所述振动信号序列的主频f0
步骤4,计算所述振动信号序列的衰减系数η0
步骤5,根据所述幅度参数、主频、衰减系数,建立所述振动信号的模型。
一种变压器振动信号的模拟装置,包括:
输入单元,输入实测的振动信号序列P=[P1,P2,…,PN],N为振动信号序列的长度;
第一计算单元,计算所述振动信号序列的幅度参数A0
第二计算单元,计算所述振动信号序列的主频f0
第三计算单元,计算所述振动信号序列的衰减系数η0
建模单元,根据所述幅度参数、主频、衰减系数,建立所述振动信号的模型。
本发明提出一种变压器振动信号模拟方案,提出一种振动信号模型,并利用实测得到的数据确定振动信号模型中的参数。利用这种方法,可以有效地模拟变压器振动信号,原理简单易行,精度高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本发明的实施例变压器振动信号的模拟方法的流程示意图;
图2为本发明应用场景的变压器振动信号的模拟方法的流程示意图;
图3为本发明的实施例变压器振动信号的模拟装置的连接示意图。
具体实施方式
下面结合附图对本发明实施例进行详细描述。
应当明确,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
为了描述的方便,描述以上装置是以功能分为各种单元/模块分别描述。当然,在实施本发明时可以把各单元/模块的功能在同一个或多个软件和/或硬件中实现。
如图1所示,为本发明所述的一种变压器振动信号的模拟方法,包括:
步骤1,输入实测的振动信号序列P=[P1,P2,…,PN],N为振动信号序列的长度;
步骤2,计算所述振动信号序列的幅度参数A0:所述步骤2包括:
A0=m0+10.5σ;
Figure BDA0002079755930000041
Figure BDA0002079755930000042
步骤3,计算所述振动信号序列的主频f0;所述步骤3包括:
f0=max[FFT(P)]
其中FFT(P)表示对所述振动信号序列进行快速傅立叶变换。
步骤4,计算所述振动信号序列的衰减系数η0;所述步骤4包括:
Figure BDA0002079755930000043
max[|P|]表示对振动信号序列P中的每个元素取绝对值,并求绝对值的最大值;
min[|P|]表示对振动信号序列P中的每个元素取绝对值,并求绝对值的最小值;
nmax表示绝对值的最大值所对应的序号;
nmin表示绝对值的最小值所对应的序号。
步骤5,根据所述幅度参数、主频、衰减系数,建立所述振动信号的模型。
所述步骤5包括:
Figure BDA0002079755930000044
x(t)为模拟的变压器的振动信号;t表示时间,e为自然指数。
本发明提出一种变压器振动信号模拟方法,此模拟方法提出一种振动信号模型,并利用实测得到的数据确定振动信号模型中的参数。利用这种方法,可以有效地模拟变压器振动信号,并可以评估变压器运行状态监测算法在不同信噪比下的检测性能。所提出的方法原理简单易行,精度高。
如图2所示,为本发明所述的变压器振动信号模拟方法的应用场景,所述方法包括:
1.输入数据
输入实测的振动信号序列P=[P1,P2,…,PN],N为振动信号序列的长度。
2.确定幅度参数A0
A0=m0+10.5σ
Figure BDA0002079755930000051
Figure BDA0002079755930000052
3.确定主频f0
f0=max[FFT(P)]
其中FFT(P)表示对振动信号进行快速傅立叶变换。
4.确定衰减系数η0
Figure BDA0002079755930000053
max[|P|]:对振动信号矢量P中的每个元素取绝对值,并求绝对值的最大值
min[|P|]:对振动信号矢量P中的每个元素取绝对值,并求绝对值的最小值
nmax:绝对值的最大值所对应的序号
nmin:绝对值的最小值所对应的序号
5.确定振动信号模型
Figure BDA0002079755930000054
x(t)即为模拟的变压器振动信号
如图3所示,为本发明所述的一种变压器振动信号的模拟装置,包括:
输入单元31,输入实测的振动信号序列P=[P1,P2,…,PN],N为振动信号序列的长度;
第一计算单元32,计算所述振动信号序列的幅度参数A0
第二计算单元33,计算所述振动信号序列的主频f0
第三计算单元34,计算所述振动信号序列的衰减系数η0
建模单元35,根据所述幅度参数、主频、衰减系数,建立所述振动信号的模型。
综上所述,变压器运行状态监测的常用方法是振声检测,此方法利用了变压器发出的振动信号,振动信号很容易受到环境噪声的影响,环境噪声对此检测方法的性能有较大的影响,因此在不同信噪比和不同噪声类型下评估振声检测算法性能非常重要,这就需要对变压器振动信号进行有效模拟,但现在缺少变压器振动信号的模拟方法。
本发明提出一种变压器振动信号模拟方法,此模拟方法提出一种振动信号模型,并利用实测得到的数据确定振动信号模型中的参数。利用这种方法,可以有效地模拟变压器振动信号,并可以评估变压器运行状态监测算法在不同信噪比下的检测性能。所提出的方法原理简单易行,精度高。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random AccessMemory,RAM)等。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。

Claims (2)

1.一种变压器振动信号的模拟方法,其特征在于,包括:
步骤1,输入实测的振动信号序列P=[P1,P2,…,PN],N为振动信号序列的长度;
步骤2,计算所述振动信号序列的幅度参数A0
步骤3,计算所述振动信号序列的主频f0
步骤4,计算所述振动信号序列的衰减系数η0
步骤5,根据所述幅度参数、主频、衰减系数,建立所述振动信号的模型;
其中,所述步骤2包括:
A0=m0+10.5σ;
Figure FDA0002984977330000011
Figure FDA0002984977330000012
其中,所述步骤3包括:
f0=max[FFT(P)]
其中FFT(P)表示对所述振动信号序列进行快速傅立叶变换;
其中,所述步骤4包括:
Figure FDA0002984977330000013
max[|P|]表示对振动信号序列P中的每个元素取绝对值,并求绝对值的最大值;
min[|P|]表示对振动信号序列P中的每个元素取绝对值,并求绝对值的最小值;
nmax表示绝对值的最大值所对应的序号;
nmin表示绝对值的最小值所对应的序号;
其中,所述步骤5包括:
Figure FDA0002984977330000014
x(t)为模拟的变压器的振动信号;t表示时间,e为自然指数。
2.一种变压器振动信号的模拟装置,其特征在于,包括:
输入单元,输入实测的振动信号序列P=[P1,P2,…,PN],N为振动信号序列的长度;
第一计算单元,计算所述振动信号序列的幅度参数A0
第二计算单元,计算所述振动信号序列的主频f0
第三计算单元,计算所述振动信号序列的衰减系数η0
建模单元,根据所述幅度参数、主频、衰减系数,建立所述振动信号的模型;
其中,所述第一计算单元用于:
A0=m0+10.5σ;
Figure FDA0002984977330000021
Figure FDA0002984977330000022
其中,所述第二计算单元用于:
f0=max[FFT(P)]
其中FFT(P)表示对所述振动信号序列进行快速傅立叶变换;
其中,所述第三计算单元用于:
Figure FDA0002984977330000023
max[|P|]表示对振动信号序列P中的每个元素取绝对值,并求绝对值的最大值;
min[|P|]表示对振动信号序列P中的每个元素取绝对值,并求绝对值的最小值;
nmax表示绝对值的最大值所对应的序号;
nmin表示绝对值的最小值所对应的序号;
其中,所述建模单元用于:
Figure FDA0002984977330000024
x(t)为模拟的变压器的振动信号;t表示时间,e为自然指数。
CN201910467079.1A 2019-05-31 2019-05-31 一种变压器振动信号的模拟方法和装置 Expired - Fee Related CN110081968B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910467079.1A CN110081968B (zh) 2019-05-31 2019-05-31 一种变压器振动信号的模拟方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910467079.1A CN110081968B (zh) 2019-05-31 2019-05-31 一种变压器振动信号的模拟方法和装置

Publications (2)

Publication Number Publication Date
CN110081968A CN110081968A (zh) 2019-08-02
CN110081968B true CN110081968B (zh) 2021-06-11

Family

ID=67422826

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910467079.1A Expired - Fee Related CN110081968B (zh) 2019-05-31 2019-05-31 一种变压器振动信号的模拟方法和装置

Country Status (1)

Country Link
CN (1) CN110081968B (zh)

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010271073A (ja) * 2009-05-19 2010-12-02 Nissin Electric Co Ltd 機器の異常診断装置
JP5582063B2 (ja) * 2011-02-21 2014-09-03 Jfeスチール株式会社 電力変換器の故障診断方法、及び故障診断装置
JP5631915B2 (ja) * 2012-03-29 2014-11-26 株式会社東芝 音声合成装置、音声合成方法、音声合成プログラムならびに学習装置
CN105973621A (zh) * 2016-05-02 2016-09-28 国家电网公司 一种基于异常振动分析的gis机械故障诊断方法和系统
CN107101714A (zh) * 2017-05-09 2017-08-29 华北电力大学(保定) 一种基于多测点振动信号特征的变压器健康状况评估方法
CN107271809A (zh) * 2017-05-18 2017-10-20 国家电网公司 一种面向大数据应用的电力设备状态量动态阈值获取方法
CN107894969A (zh) * 2017-09-13 2018-04-10 中国石油大学(华东) 一种基于趋势分析的变压器潜伏性故障预警方法
CN107907778A (zh) * 2017-10-31 2018-04-13 华北电力大学(保定) 一种基于多特征参量的变压器综合故障诊断方法
CN107992804A (zh) * 2017-11-16 2018-05-04 华北电力大学(保定) 一种基于改进多尺度熵的变压器振动信号特征提取方法
CN108664741A (zh) * 2018-05-14 2018-10-16 平顶山学院 基于时间序列模型特征的变电站故障检测方法
CN109359271A (zh) * 2018-12-21 2019-02-19 浙江大学 一种基于逻辑回归的变压器绕组变形程度在线检测方法
CN109443528A (zh) * 2018-11-16 2019-03-08 国网江苏省电力有限公司盐城供电分公司 一种基于振动信号分析的变压器故障诊断系统及其诊断方法
CN109597967A (zh) * 2018-11-20 2019-04-09 江苏云上电力科技有限公司 一种电力系统配变电站负荷数据异常检测与修复方法
CN109708748A (zh) * 2019-02-01 2019-05-03 国网山东省电力公司昌乐县供电公司 变电站gis组合电器振动及噪声异常故障定位方法

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010271073A (ja) * 2009-05-19 2010-12-02 Nissin Electric Co Ltd 機器の異常診断装置
JP5582063B2 (ja) * 2011-02-21 2014-09-03 Jfeスチール株式会社 電力変換器の故障診断方法、及び故障診断装置
JP5631915B2 (ja) * 2012-03-29 2014-11-26 株式会社東芝 音声合成装置、音声合成方法、音声合成プログラムならびに学習装置
CN105973621A (zh) * 2016-05-02 2016-09-28 国家电网公司 一种基于异常振动分析的gis机械故障诊断方法和系统
CN107101714A (zh) * 2017-05-09 2017-08-29 华北电力大学(保定) 一种基于多测点振动信号特征的变压器健康状况评估方法
CN107271809A (zh) * 2017-05-18 2017-10-20 国家电网公司 一种面向大数据应用的电力设备状态量动态阈值获取方法
CN107894969A (zh) * 2017-09-13 2018-04-10 中国石油大学(华东) 一种基于趋势分析的变压器潜伏性故障预警方法
CN107907778A (zh) * 2017-10-31 2018-04-13 华北电力大学(保定) 一种基于多特征参量的变压器综合故障诊断方法
CN107992804A (zh) * 2017-11-16 2018-05-04 华北电力大学(保定) 一种基于改进多尺度熵的变压器振动信号特征提取方法
CN108664741A (zh) * 2018-05-14 2018-10-16 平顶山学院 基于时间序列模型特征的变电站故障检测方法
CN109443528A (zh) * 2018-11-16 2019-03-08 国网江苏省电力有限公司盐城供电分公司 一种基于振动信号分析的变压器故障诊断系统及其诊断方法
CN109597967A (zh) * 2018-11-20 2019-04-09 江苏云上电力科技有限公司 一种电力系统配变电站负荷数据异常检测与修复方法
CN109359271A (zh) * 2018-12-21 2019-02-19 浙江大学 一种基于逻辑回归的变压器绕组变形程度在线检测方法
CN109708748A (zh) * 2019-02-01 2019-05-03 国网山东省电力公司昌乐县供电公司 变电站gis组合电器振动及噪声异常故障定位方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《Evolution of transformer health index in the form of mathematical equation》;A.AzmiJ.Jasni;《Renewable and Sustainable Energy Reviews》;20170930;全文 *
《基于可听声的变压器故障诊断技术综述》;谢荣斌;《宁夏电力》;20170428;全文 *
《基于盲源分离的电力变压器振声自适应提取与异常状态检测方法》;刘晗;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20190115;第9页-第59页 *
刘勇.《 基于动态阈值的变压器异常状态检测》.《电测与仪表》.2017, *

Also Published As

Publication number Publication date
CN110081968A (zh) 2019-08-02

Similar Documents

Publication Publication Date Title
CN110031089B (zh) 一种变压器运行状态振声检测信号的滤波方法和装置
CN110017894B (zh) 变压器运行状态振声检测中随机噪声的滤除方法和装置
CN110703149A (zh) 一种利用字符间距的变压器运行状态振声检测方法和系统
CN110081968B (zh) 一种变压器振动信号的模拟方法和装置
CN110545086A (zh) 一种利用全局优化的变压器振声信号滤波方法和系统
CN111664934A (zh) 一种利用特征选择的变压器状态振声检测信号滤波方法和系统
CN111780868A (zh) 一种利用Jeffery差异量的变压器运行状态振声检测方法和系统
CN110703145A (zh) 一种利用多优化理论的变压器振声信号重构方法和系统
CN110286289B (zh) 一种变压器振声检测信号滤波方法
CN110646691B (zh) 一种利用拉伸变换的变压器振声信号滤波方法和系统
CN110286287B (zh) 一种基于小波变换的变压器运行状态振声检测信号滤波方法和系统
CN110161330B (zh) 基于灰度理论的变压器运行状态的振声检测方法以及装置
CN110514295B (zh) 一种利用svd分解的变压器运行状态振声检测信号滤波方法和系统
CN110161363B (zh) 基于主频表征量的变压器运行状态振声检测方法和系统
CN112304419A (zh) 一种利用泛化稀疏编码的振声检测信号重构方法和系统
CN111751098A (zh) 一种利用高斯预测模型的振声检测信号重构方法和系统
CN111649819A (zh) 一种利用迭代软阈值的变压器状态振声检测信号滤波方法和系统
CN111141384A (zh) 一种利用Frechet正则化的变压器状态振声检测信号重构方法和系统
CN110017892B (zh) 一种变压器的异常状态振声的检测方法和装置
CN112345226B (zh) 一种利用块协调最小化的振声检测信号重构方法和系统
CN112417994B (zh) 一种利用正则化因子的振声检测信号滤波方法和系统
CN110545087A (zh) 一种利用非线性正则化的变压器振声信号滤波方法和系统
CN110657881B (zh) 一种利用稀疏反演的变压器振声信号滤波方法和系统
CN110703144B (zh) 基于离散余弦变换的变压器运行状态检测方法和系统
CN112284520B (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: 20210611

CF01 Termination of patent right due to non-payment of annual fee