CN110657881B - Transformer vibration sound signal filtering method and system by utilizing sparse inversion - Google Patents
Transformer vibration sound signal filtering method and system by utilizing sparse inversion Download PDFInfo
- Publication number
- CN110657881B CN110657881B CN201910867657.0A CN201910867657A CN110657881B CN 110657881 B CN110657881 B CN 110657881B CN 201910867657 A CN201910867657 A CN 201910867657A CN 110657881 B CN110657881 B CN 110657881B
- Authority
- CN
- China
- Prior art keywords
- matrix
- sparse
- signal sequence
- specifically
- transformer
- 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
Links
Images
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The embodiment of the invention discloses a transformer vibration acoustic signal filtering method and a system by utilizing sparse inversion, wherein the method comprises the following steps: step 1, inputting an actually measured vibration sound signal sequence S; step 2, carrying out noise filtering processing on the signal sequence S according to a sparse inversion theory, wherein the signal sequence after noise filtering is SNEW(ii) a In particular to a method for preparing a high-performance nano-silver alloy,
Description
Technical Field
The invention relates to the field of electric power, in particular to a method and a system for filtering a vibration sound signal of a transformer.
Background
With the high-speed development of the smart grid, the safe and stable operation of the power equipment is particularly important. At present, the detection of the operating state of the power equipment with ultrahigh voltage and above voltage grades, especially the detection of the abnormal state, is increasingly important and urgent. As an important component of an electric power system, a power transformer is one of the most important electrical devices in a substation, and its reliable operation is related to the safety of a power grid. Generally, the abnormal state of the transformer can be divided into core abnormality and winding abnormality. The core abnormality is mainly represented by core saturation, and the winding abnormality generally includes winding deformation, winding looseness and the like.
The basic principle of the transformer abnormal state detection is to extract each characteristic quantity in the operation of the transformer, analyze, identify and track the characteristic quantity so as to monitor the abnormal operation state of the transformer. The detection method can be divided into invasive detection and non-invasive detection according to the contact degree; the detection can be divided into live detection and power failure detection according to whether the shutdown detection is needed or not; the method can be classified into an electrical quantity method, a non-electrical quantity method, and the like according to the type of the detected quantity. In comparison, the non-invasive detection has strong transportability and is more convenient to install; the live detection does not affect the operation of the transformer; the non-electric quantity method is not electrically connected with the power system, so that the method is safer. The current common detection methods for the operation state of the transformer include a pulse current method and an ultrasonic detection method for detecting partial discharge, a frequency response method for detecting winding deformation, a vibration detection method for detecting mechanical and electrical faults, and the like. The detection methods mainly detect the insulation condition and the mechanical structure condition of the transformer, wherein the detection of the vibration signal (vibration sound) of the transformer is the most comprehensive, and the fault and the abnormal state of most transformers can be reflected.
In the running process of the transformer, the magnetostriction of the iron core silicon steel sheets and the vibration caused by the winding electrodynamic force can radiate vibration sound signals with different amplitudes and frequencies to the periphery. When the transformer normally operates, uniform low-frequency noise is emitted outwards; if the sound is not uniform, it is not normal. The transformer can make distinctive sounds in different running states, and the running state of the transformer can be mastered by detecting the sounds made by the transformer. It is worth noting that the detection of the sound emitted by the transformer in different operating states not only can detect a plurality of serious faults causing the change of the electrical quantity, but also can detect a plurality of abnormal states which do not endanger the insulation and do not cause the change of the electrical quantity, such as the loosening of internal and external parts of the transformer, and the like.
Because the vibration sound detection method utilizes the vibration signal sent by the transformer, the vibration sound detection method is easily influenced by environmental noise, and therefore, how to effectively identify the vibration sound and the noise is the key for success of the method. The existing common method has insufficient attention to the problem, and no effective measure is taken to solve the problem.
Disclosure of Invention
The invention aims to provide a transformer vibration acoustic signal filtering method and a transformer vibration acoustic signal filtering system by utilizing sparse inversion. The method has better robustness and simpler calculation.
In order to achieve the purpose, the invention provides the following scheme:
a transformer vibro-acoustic signal filtering method utilizing sparse inversion comprises the following steps:
step 1, inputting an actually measured vibration sound signal sequence S;
step 2, according to sparse inversion theory, placeThe signal sequence S is subjected to noise filtering processing, and the signal sequence after noise filtering is SNEW(ii) a In particular to a method for preparing a high-performance nano-silver alloy,wherein mu is a sparsity factor; g is a sparse matrix; r is a sparse reflection vector; kappa is a sparse adjustment factor; d is a Wolk matrix; p is the intermediate parameter matrix.
A transformer vibro-acoustic signal filtering system using sparse inversion, comprising:
the acquisition module inputs an actually measured vibration sound signal sequence S;
the filtering module is used for carrying out noise filtering processing on the signal sequence S according to a sparse inversion theory, and the signal sequence after noise filtering is SNEW(ii) a In particular to a method for preparing a high-performance nano-silver alloy,wherein mu is a sparsity factor; g is a sparse matrix; r is a sparse reflection vector; kappa is a sparse adjustment factor; d is a Wolk matrix; p is the intermediate parameter matrix.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
although the transformer vibration and sound detection method is widely applied to monitoring the running state of the transformer and the technology is relatively mature, the vibration and sound detection method utilizes the vibration signal sent by the transformer and is easily influenced by the environmental noise, so that the method often cannot obtain satisfactory results when being applied in the actual working environment.
The invention aims to provide a transformer vibration acoustic signal filtering method and a transformer vibration acoustic signal filtering system by utilizing sparse inversion. The method has better robustness and simpler calculation.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments will be briefly described below. It is obvious that the drawings in the following description are only some embodiments of the invention, and that for a person skilled in the art, other drawings can be derived from them without inventive effort.
FIG. 1 is a schematic flow diagram of the process of the present invention;
FIG. 2 is a schematic diagram of the system of the present invention;
FIG. 3 is a flow chart illustrating an embodiment of the present invention.
Detailed Description
The technical solution in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. It is to be understood that the described embodiments are merely exemplary of the invention, and not restrictive of the full scope of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
FIG. 1 is a schematic flow chart of a transformer vibro-acoustic signal filtering method using sparse inversion
Fig. 1 is a schematic flow chart of a transformer vibro-acoustic signal filtering method using sparse inversion according to the present invention. As shown in fig. 1, the transformer vibro-acoustic signal filtering method using sparse inversion specifically includes the following steps:
step 1, inputting an actually measured vibration sound signal sequence S;
step 2, carrying out noise filtering processing on the signal sequence S according to a sparse inversion theory, wherein the signal sequence after noise filtering is SNEW(ii) a In particular to a method for preparing a high-performance nano-silver alloy,wherein mu is a sparsity factor; g is diluteA sparse matrix; r is a sparse reflection vector; kappa is a sparse adjustment factor; d is a Wolk matrix; p is the intermediate parameter matrix.
Before the step 2, the method further comprises:
and 3, solving the sparsity factor mu, the sparse matrix G, the sparse reflection vector R, the sparse adjustment factor kappa and the Wolk matrix D.
The step 3 comprises the following steps:
step 301, obtaining a cyclic delay matrix DCThe method specifically comprises the following steps:
wherein:
sn: the nth element [ N ═ 1,2, …, N of the signal sequence S]
N: length of the signal sequence S
Step 302, obtaining the sparsity factor μ, specifically:
wherein:
Step 303, obtaining the sparse matrix G, specifically:
G=kron(DC,I)
wherein:
kron(DCand I): matrix DCKronecker multiplication with I
I: unit matrix
Step 304, solving the sparse reflection vector R, specifically:
wherein:
Step 305, obtaining the sparse adjustment factor κ, specifically:
wherein:
mS: mean value of the signal sequence S
σS: mean square error of the signal sequence S
Step 306, obtaining the wacker matrix D, specifically:
wherein
Δ=max[|s1-s2|,|s2-s3|,…,|sN-s1|]-min[|s1-s2|,|s2-s3|,…,|sN-s1|]
FIG. 2 structural intention of transformer vibro-acoustic signal filtering system using sparse inversion
Fig. 2 is a schematic structural diagram of a transformer vibro-acoustic signal filtering system using sparse inversion according to the present invention. As shown in fig. 2, the transformer vibro-acoustic signal filtering system using sparse inversion includes the following structure:
the acquisition module 401 inputs an actually measured vibration and sound signal sequence S;
a filtering module 402, configured to perform noise filtering processing on the signal sequence S according to a sparse inversion theory, where the signal sequence after noise filtering is SNEW(ii) a In particular to a method for preparing a high-performance nano-silver alloy,wherein mu is a sparsity factor; g is a sparse matrix; r is a sparse reflection vector; kappa is a sparse adjustment factor; d is a Wolk matrix; p is the intermediate parameter matrix.
The system further comprises:
a calculation module 403 for obtaining the sparse inversion C and the optimal prediction vector mOPTShaping matrix B, system matrix L and aliasing matrix
The calculation module 403 includes the following units:
delay unit 4031 for obtaining cyclic delay matrix DCThe method specifically comprises the following steps:
wherein:
sn: the nth element [ N ═ 1,2, …, N of the signal sequence S]
N: length of the signal sequence S
The first calculation unit 4032 calculates the sparsity factor μ, specifically:
wherein:
The second calculation unit 4033, which calculates the sparse matrix G, specifically is:
G=kron(DC,I)
wherein:
kron(DCand I): matrix DCKronecker multiplication with I
I: unit matrix
The third calculation unit 4034, which calculates the sparse reflection vector R, specifically is:
wherein:
The fourth calculating unit 4035, which calculates the sparse adjustment factor κ specifically is:
wherein:
mS: mean value of the signal sequence S
σS: mean square error of the signal sequence S
The fifth calculation unit 4036, which calculates the wacker matrix D, specifically is:
wherein
Δ=max[|s1-s2|,|s2-s3|,…,|sN-s1|]-min[|s1-s2|,|s2-s3|,…,|sN-s1|]
The following provides an embodiment for further illustrating the invention
FIG. 3 is a flow chart illustrating an embodiment of the present invention. As shown in fig. 3, the method specifically includes the following steps:
1. inputting measured vibration sound signal sequence
S=[s1,s2,…,sN-1,sN]
Wherein:
s: measured PLC signal data sequence with length N
siI is 1,2, …, N is measured PLC signal with serial number i
2. Determining a cyclic delay matrix
Wherein:
sn: the nth element [ N ═ 1,2, …, N of the signal sequence S]
N: length of the signal sequence S
3. Calculating a sparsity factor
Wherein:
4. Obtaining sparse matrices
G=kron(DC,I)
Wherein:
kron(DCand I): matrix DCKronecker multiplication with I
I: unit matrix
5. Finding sparse reflection vectors
Wherein:
6. Finding sparse adjustment factors
Wherein:
mS: mean value of the signal sequence S
σS: mean square error of the signal sequence S
7. Obtaining Wolk matrix
Wherein
Δ=max[|s1-s2|,|s2-s3|,…,|sN-s1|]-min[|s1-s2|,|s2-s3|,…,|sN-s1|]
8. Filtering
And carrying out noise filtering processing on the signal sequence S according to a sparse inversion theory, wherein the signal sequence after noise filtering is SNEW(ii) a In particular to a method for preparing a high-performance nano-silver alloy,wherein mu is a sparsity factor; g is a sparse matrix; r is a sparse reflection vector; kappa is a sparse adjustment factor; d is a Wolk matrix; p is the intermediate parameter matrix.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is simple because the system corresponds to the method disclosed by the embodiment, and the relevant part can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.
Claims (1)
1. A transformer vibro-acoustic signal filtering method utilizing sparse inversion is characterized by comprising the following steps:
step 1, inputting an actually measured vibration sound signal sequence S;
step 2, obtaining a cyclic delay matrix DCThe method specifically comprises the following steps:
wherein:
sn: the nth element [ N ═ 1,2, …, N of the signal sequence S];
N: the length of the signal sequence S;
step 3, obtaining a sparsity factor mu, specifically:
wherein:
step 4, solving a sparse matrix G, specifically:
G=kron(DC,I);
wherein:
kron(DCand I): matrix DCA kronecker product operation with I;
i: an identity matrix;
step 5, solving a sparse reflection vector R, specifically:
wherein:
step 6, obtaining a sparse adjustment factor k, specifically:
wherein:
mS: a mean value of the signal sequence S;
σS: the mean square error of the signal sequence S;
step 7, solving a Wolk matrix D, specifically:
wherein:
Δ=max[|s1-s2|,|s2-s3|,…,|sN-s1|]-min[|s1-s2|,|s2-s3|,…,|sN-s1|];
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910867657.0A CN110657881B (en) | 2019-09-14 | 2019-09-14 | Transformer vibration sound signal filtering method and system by utilizing sparse inversion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910867657.0A CN110657881B (en) | 2019-09-14 | 2019-09-14 | Transformer vibration sound signal filtering method and system by utilizing sparse inversion |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110657881A CN110657881A (en) | 2020-01-07 |
CN110657881B true CN110657881B (en) | 2021-04-06 |
Family
ID=69037068
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910867657.0A Expired - Fee Related CN110657881B (en) | 2019-09-14 | 2019-09-14 | Transformer vibration sound signal filtering method and system by utilizing sparse inversion |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110657881B (en) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010271073A (en) * | 2009-05-19 | 2010-12-02 | Nissin Electric Co Ltd | Diagnosis device of abnormality in equipment |
JP5582063B2 (en) * | 2011-02-21 | 2014-09-03 | Jfeスチール株式会社 | Power converter failure diagnosis method and failure diagnosis apparatus |
JP5631915B2 (en) * | 2012-03-29 | 2014-11-26 | 株式会社東芝 | Speech synthesis apparatus, speech synthesis method, speech synthesis program, and learning apparatus |
CN105973621A (en) * | 2016-05-02 | 2016-09-28 | 国家电网公司 | Abnormal vibration analysis-based GIS (gas insulated switchgear) mechanical fault diagnosis method and system |
CN107101714A (en) * | 2017-05-09 | 2017-08-29 | 华北电力大学(保定) | A kind of transformer health Evaluation method based on multi-measuring point vibration signal characteristics |
CN107271809A (en) * | 2017-05-18 | 2017-10-20 | 国家电网公司 | A kind of status of electric power amount dynamic threshold acquisition methods applied towards big data |
CN107894969A (en) * | 2017-09-13 | 2018-04-10 | 中国石油大学(华东) | A kind of latent transformer fault early warning method based on trend analysis |
CN107907778A (en) * | 2017-10-31 | 2018-04-13 | 华北电力大学(保定) | A kind of Synthesized Diagnosis On Transformer Faults method based on multiple features parameter |
CN108664741A (en) * | 2018-05-14 | 2018-10-16 | 平顶山学院 | Substation fault detection method based on time series models feature |
CN109359271A (en) * | 2018-12-21 | 2019-02-19 | 浙江大学 | A kind of deformation of transformer winding degree online test method that logic-based returns |
CN109443528A (en) * | 2018-11-16 | 2019-03-08 | 国网江苏省电力有限公司盐城供电分公司 | A kind of transformer fault diagnosis system and its diagnostic method based on analysis of vibration signal |
CN109597967A (en) * | 2018-11-20 | 2019-04-09 | 江苏云上电力科技有限公司 | A kind of electric system distribution transforming power station load data abnormality detection and restorative procedure |
CN109708748A (en) * | 2019-02-01 | 2019-05-03 | 国网山东省电力公司昌乐县供电公司 | The vibration of substation's GIS combination electric appliance and noise abnormal failure localization method |
CN110031089A (en) * | 2019-05-15 | 2019-07-19 | 广东石油化工学院 | A kind of filtering method and device of running state of transformer vibration sound detection signal |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103558633A (en) * | 2013-10-21 | 2014-02-05 | 中国海洋石油总公司 | Interlayer multiple suppression method based on sparse inversion |
CN104749631B (en) * | 2015-03-11 | 2017-02-08 | 中国科学院地质与地球物理研究所 | Sparse inversion based migration velocity analysis method and device |
-
2019
- 2019-09-14 CN CN201910867657.0A patent/CN110657881B/en not_active Expired - Fee Related
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010271073A (en) * | 2009-05-19 | 2010-12-02 | Nissin Electric Co Ltd | Diagnosis device of abnormality in equipment |
JP5582063B2 (en) * | 2011-02-21 | 2014-09-03 | Jfeスチール株式会社 | Power converter failure diagnosis method and failure diagnosis apparatus |
JP5631915B2 (en) * | 2012-03-29 | 2014-11-26 | 株式会社東芝 | Speech synthesis apparatus, speech synthesis method, speech synthesis program, and learning apparatus |
CN105973621A (en) * | 2016-05-02 | 2016-09-28 | 国家电网公司 | Abnormal vibration analysis-based GIS (gas insulated switchgear) mechanical fault diagnosis method and system |
CN107101714A (en) * | 2017-05-09 | 2017-08-29 | 华北电力大学(保定) | A kind of transformer health Evaluation method based on multi-measuring point vibration signal characteristics |
CN107271809A (en) * | 2017-05-18 | 2017-10-20 | 国家电网公司 | A kind of status of electric power amount dynamic threshold acquisition methods applied towards big data |
CN107894969A (en) * | 2017-09-13 | 2018-04-10 | 中国石油大学(华东) | A kind of latent transformer fault early warning method based on trend analysis |
CN107907778A (en) * | 2017-10-31 | 2018-04-13 | 华北电力大学(保定) | A kind of Synthesized Diagnosis On Transformer Faults method based on multiple features parameter |
CN108664741A (en) * | 2018-05-14 | 2018-10-16 | 平顶山学院 | Substation fault detection method based on time series models feature |
CN109443528A (en) * | 2018-11-16 | 2019-03-08 | 国网江苏省电力有限公司盐城供电分公司 | A kind of transformer fault diagnosis system and its diagnostic method based on analysis of vibration signal |
CN109597967A (en) * | 2018-11-20 | 2019-04-09 | 江苏云上电力科技有限公司 | A kind of electric system distribution transforming power station load data abnormality detection and restorative procedure |
CN109359271A (en) * | 2018-12-21 | 2019-02-19 | 浙江大学 | A kind of deformation of transformer winding degree online test method that logic-based returns |
CN109708748A (en) * | 2019-02-01 | 2019-05-03 | 国网山东省电力公司昌乐县供电公司 | The vibration of substation's GIS combination electric appliance and noise abnormal failure localization method |
CN110031089A (en) * | 2019-05-15 | 2019-07-19 | 广东石油化工学院 | A kind of filtering method and device of running state of transformer vibration sound detection signal |
Non-Patent Citations (4)
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 |
---|---|
CN110657881A (en) | 2020-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110703149B (en) | Method and system for detecting vibration and sound of running state of transformer by utilizing character spacing | |
CN110703145B (en) | Transformer vibration sound signal reconstruction method and system by using multiple optimization theories | |
CN110545086A (en) | Transformer vibration sound signal filtering method and system by utilizing global optimization | |
CN111665405A (en) | Vibration and sound detection signal filtering method and system based on sparsity minimization | |
CN110017894B (en) | Method and device for filtering random noise in vibration and sound detection of transformer in running state | |
CN111780868A (en) | Transformer running state vibration and noise detection method and system by utilizing Jeffery difference | |
CN111664933A (en) | Method and system for filtering vibration and sound detection signal by utilizing static vector optimization | |
CN111664934A (en) | Transformer state vibration and sound detection signal filtering method and system using feature selection | |
CN110646691B (en) | Transformer vibration sound signal filtering method and system by utilizing stretching transformation | |
CN110657881B (en) | Transformer vibration sound signal filtering method and system by utilizing sparse inversion | |
CN110514295B (en) | Transformer running state vibration and sound detection signal filtering method and system by utilizing SVD (singular value decomposition) | |
CN110286289B (en) | Filtering method for vibration and sound detection signal of transformer | |
CN110286287B (en) | Wavelet transform-based method and system for filtering vibration and sound detection signals of running state of transformer | |
CN111751098A (en) | Vibration and sound detection signal reconstruction method and system by using Gaussian prediction model | |
CN110161363B (en) | Transformer running state vibration and sound detection method and system based on main frequency characteristic quantity | |
CN112304419A (en) | Vibration and sound detection signal reconstruction method and system by using generalized sparse coding | |
CN111141384A (en) | Transformer state vibration and sound detection signal reconstruction method and system by utilizing Frechet regularization | |
CN110286291B (en) | Method and system for detecting vibration and sound of running state of transformer by using principal components | |
CN111649819A (en) | Transformer state vibration and sound detection signal filtering method and system using iteration soft threshold | |
CN110837013A (en) | Transformer state vibration and sound detection signal reconstruction method and system represented by sparse dictionary | |
CN111561992A (en) | Method and system for detecting vibration sound of running state of transformer by using B sampling | |
CN110545088A (en) | Transformer vibration sound signal filtering method and system by using dual variables | |
CN112417994B (en) | Vibration and sound detection signal filtering method and system using regularization factor | |
CN111473861A (en) | Transformer state vibration and sound detection signal reconstruction method and system by using sparse errors | |
CN110545087A (en) | Transformer vibration sound signal filtering method and system by utilizing nonlinear regularization |
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: 20210406 Termination date: 20210914 |
|
CF01 | Termination of patent right due to non-payment of annual fee |