CN107132451A - The winding state detection method and system of transformer - Google Patents
The winding state detection method and system of transformer Download PDFInfo
- Publication number
- CN107132451A CN107132451A CN201710400131.2A CN201710400131A CN107132451A CN 107132451 A CN107132451 A CN 107132451A CN 201710400131 A CN201710400131 A CN 201710400131A CN 107132451 A CN107132451 A CN 107132451A
- Authority
- CN
- China
- Prior art keywords
- mrow
- transformer
- frequency domain
- domain sequence
- matrix
- 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
- 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/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
- G01R31/72—Testing of electric windings
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B17/00—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations
- G01B17/04—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations for measuring the deformation in a solid, e.g. by vibrating string
Abstract
The present invention relates to the winding state detection method and system of a kind of transformer.The winding state detection method of above-mentioned transformer, comprises the following steps:Gather test vibration signal of the transformer in the early stage in test process, the first frequency domain sequence matrix of transformer according to the test vibration signal generation;Gather operational shock signal of the transformer in practical work process, the second frequency domain sequence matrix of transformer according to the operational shock signal generation;The first frequency domain sequence matrix and the second frequency domain sequence matrix are subjected to correlation analysis, the correlation matrix and amplitude coefficient matrix of transformer is obtained;Calculate the ∞ norm value S of the correlation matrix, and the amplitude coefficient matrix ∞ norm values F;The winding state of transformer is detected according to the ∞ norm values S and ∞ norm value F.The data that the detection process of above-mentioned transformer winding state is covered more fully, effectively increase the Detection results of deformation of transformer winding state-detection.
Description
Technical field
The present invention relates to technical field of electric power, the winding state detection method and system of more particularly to a kind of transformer.
Background technology
One of equipment particularly significant and expensive in power system during transformer.Its operation conditions not only influences itself
Safety, and affect the stability and reliability of whole Operation of Electric Systems.For a long time, transformer is safe and reliable
Operation is constantly subjected to operation power and the most attention of administrative department, and this is also that system is safe and stable and economical operation important
Index.With the rapid development of the national economy, people are increasing to the demand of electricity, the effect that transformer is played also increasingly is weighed
Will, and develop towards the bigger direction of voltage class and capacity.
The winding of transformer is when it occurs short trouble by by huge short-circuit force effect, and winding is under this masterpiece use
It can be easy to deform, collapse etc. and damaging, the stable operation to transformer causes potential potential safety hazard.Therefore recent year
The outer diagnostic method to deformation of transformer winding carries out numerous studies, mainly there is short-circuit impedance hair, frequency response method, low pressure at present
Impulse method, dissolved gas method etc., but the winding state detection method of above-mentioned traditional transformer normally only can be used for offline
Detection;And in actual use, transformer generally requires continuous service, stoppage in transit detection is difficult to, and easily influence is directed to
Transformer Winding carries out the effect of deformation state detection.
The content of the invention
Based on this, it is necessary to easily influence to carry out deformation state Detection results for Transformer Winding for traditional scheme
There is provided the winding state detection method and system of a kind of transformer for technical problem.
A kind of winding state detection method of transformer, comprises the following steps:
Test vibration signal of the transformer in the early stage in test process is gathered, according to the test vibration signal generation
First frequency domain sequence matrix of transformer;
Operational shock signal of the transformer in practical work process is gathered, according to the operational shock signal generation
Second frequency domain sequence matrix of transformer;
The first frequency domain sequence matrix and the second frequency domain sequence matrix are subjected to correlation analysis, the phase of transformer is obtained
Relation matrix number and amplitude coefficient matrix;
Calculate the ∞ norm value S of the correlation matrix, and the amplitude coefficient matrix ∞ norm values F;
The winding state of transformer is detected according to the ∞ norm values S and ∞ norm value F.
A kind of winding state detecting system of transformer, including:
First acquisition module, for gathering test vibration signal of the transformer in the early stage in test process, is surveyed according to described
Try the first frequency domain sequence matrix that vibration signal generates the transformer;
Second acquisition module, for gathering operational shock signal of the transformer in practical work process, according to the fortune
Row vibration signal generates the second frequency domain sequence matrix of the transformer;
Analysis module, for the first frequency domain sequence matrix and the second frequency domain sequence matrix to be carried out into correlation analysis,
Obtain the correlation matrix and amplitude coefficient matrix of transformer;
The ∞ of computing module, the ∞ norm value S for calculating the correlation matrix, and the amplitude coefficient matrix
Norm value F;
Detection module, the winding state for detecting transformer according to the ∞ norm values S and ∞ norm value F.
The winding state detection method and system of above-mentioned transformer, can gather survey of the transformer in the early stage in test process
Vibration signal is tried, the first frequency domain sequence matrix of transformer is generated;Gather operational shock of the transformer in practical work process
Signal, generates the second frequency domain sequence matrix of transformer;Again by the first frequency domain sequence matrix and the second frequency domain sequence matrix
Correlation analysis is carried out, the correlation matrix and amplitude coefficient matrix of transformer is obtained;To calculate the coefficient correlation square
Battle array ∞ norm value S, and the amplitude coefficient matrix ∞ norm values F;Examined according to the ∞ norm values S and ∞ norm value F
Survey the winding state of transformer;The winding state of above-mentioned transformer can be according to transformer factory data (the first frequency domain sequence
Matrix) and the course of work in service data (the second frequency domain sequence matrix) accordingly detected, the number that detection process is covered
According to the Detection results for more fully, effectively increasing deformation of transformer winding state-detection.
A kind of computer-readable recording medium, is stored thereon with computer program, and the computer program is executed by processor
The winding state detection method of Shi Shixian transformers as described above.
The computer program stored on above computer readable storage medium storing program for executing, allows the winding state of transformer according to change
Service data in the factory data and the course of work of depressor is accordingly detected, with higher Detection results.
A kind of computer equipment, including memory, processor and storage can be run on a memory and on a processor
Computer program, realizes the winding state detection method of transformer as described above during the computing device described program.
In above computer equipment, factory data (the first frequency domain sequence that the winding state of transformer can be according to transformer
Column matrix) and the course of work in service data (the second frequency domain sequence matrix) accordingly detected, what detection process was covered
Data more fully, make the Detection results for Transformer Winding progress deformation state detection higher.
Brief description of the drawings
Fig. 1 is the winding state detection method flow chart of the transformer of one embodiment;
Fig. 2 is the winding state detecting system structural representation of the transformer of one embodiment;
Fig. 3 is the computer equipment structural representation of one embodiment.
Embodiment
The winding state detection method of transformer and the embodiment of system to the present invention are made below in conjunction with the accompanying drawings
It is described in detail.
With reference to Fig. 1, Fig. 1 show the winding state detection method flow chart of the transformer of one embodiment, including as follows
Step:
S10, the test vibration signal of collection transformer in the early stage in test process, according to the test vibration signal generation
First frequency domain sequence matrix of the transformer;
Above-mentioned test process at initial stage can include Transformer Plant Test process or initial operation stage running.Specifically
Ground, it is possible to use related vibration measurement instrument gathers original of the transformer under unloaded and rated load condition in above-mentioned test process at initial stage
Beginning multiple spot vibration signal, to obtain corresponding test vibration signal.
In one embodiment, above-mentioned test process at initial stage includes delivery test process or initial operation stage running.
The test process at initial stage that the present embodiment is provided includes delivery test process or initial operation stage running, so,
It is original many under unloaded and rated load condition using related vibration measurement instrument collection transformer in above-mentioned test process at initial stage
Point vibration signal, can the correlation of more complete acquisition transformer dispatch from the factory vibration signal, make resulting test vibration signal
More accuracy and integrality.
As one embodiment, the test vibration signal of above-mentioned collection transformer in the early stage in test process, according to described
The process of first frequency domain sequence matrix of transformer described in test vibration signal generation can include:
In the early stage in test process, the multiple spot vibration signal of collection transformer tank surface in the unloaded state obtains original
Begin unloaded vibration signal sequence matrix Xmm, the multiple spot vibration signal of transformer tank surface under rated load condition is gathered, is obtained
To original load vibration signal frequency domain sequence matrix Ymm;
According to the original unloaded vibration signal sequence matrix XmmWith original load vibration signal frequency domain sequence matrix YmmMeter
Calculate the first frequency domain sequence matrix Zmm;Wherein, Zmm=Ymm-Xmm。
In the present embodiment, collection transformer, can be to above-mentioned in the unloaded state after the multiple spot vibration signal of tank surface
The multiple spot vibration signal collected under Light Condition carries out FFT (fast Fourier transform) processing, to determine original unloaded vibration
Signal sequence matrix Xmm;The multiple spot vibration signal of transformer tank surface under rated load condition is gathered, can be to above-mentioned volume
The multiple spot vibration signal collected under fixed load state carries out FFT processing, to determine original load vibration signal frequency domain sequence
Matrix Ymm.In above-mentioned original unloaded vibration signal sequence matrix Xmm, original load vibration signal frequency domain sequence matrix YmmAnd
First frequency domain sequence matrix ZmmIn, the dimension of the corresponding square formation in tank surface vibration acquisition face of subscript m indication transformer.
S20, operational shock signal of the collection transformer in practical work process, according to the operational shock signal generation
Second frequency domain sequence matrix of the transformer;
Specifically, above-mentioned steps can respectively be adopted using related vibration measurement instrument in the practical work process of above-mentioned transformer
Collect to be measured multiple spot vibration signal of the transformer under unloaded and rated load condition, to obtain corresponding operational shock signal.
In one embodiment, operational shock signal of the above-mentioned collection transformer in practical work process, according to described
The process of second frequency domain sequence matrix of transformer described in operational shock signal generation can include:
In practical work process, the multiple spot vibration signal of collection transformer tank surface in the unloaded state is treated
Survey unloaded vibration signal sequence matrix Amm, the multiple spot vibration signal of transformer tank surface under rated load condition is gathered, is obtained
To load vibration signal frequency domain sequence matrix B to be measuredmm;
According to the original unloaded vibration signal sequence matrix AmmWith original load vibration signal frequency domain sequence matrix BmmMeter
Calculate the second frequency domain sequence Matrix Cmm;Wherein, Cmm=Bmm-Amm。
In practical work process, collection transformer is in the unloaded state after the multiple spot vibration signal of tank surface, can be with
FFT processing is carried out to the multiple spot vibration signal collected under above-mentioned Light Condition, to determine unloaded vibration signal sequence to be measured
Matrix Amm;The multiple spot vibration signal of transformer tank surface under rated load condition is gathered, can be to above-mentioned nominal load shape
The multiple spot vibration signal collected under state carries out FFT processing, to determine load vibration signal frequency domain sequence matrix B to be measuredmm.
Above-mentioned unloaded vibration signal sequence matrix A to be measuredmm, it is to be measured load vibration signal frequency domain sequence matrix BmmAnd the second frequency domain sequence
Column matrix CmmIn, the dimension of the corresponding square formation in tank surface vibration acquisition face of subscript m indication transformer.
As one embodiment, it is determined that load vibration signal frequency domain sequence matrix B to be measuredmm, and be difficult to zero load to be measured
Vibration signal sequence matrix Amm, can be by above-mentioned original unloaded vibration signal sequence matrix X when accurately being determinedmmIt is defined as
Unloaded vibration signal sequence matrix A to be measuredmmApproximation, now, above-mentioned second frequency domain sequence Matrix Cmm≈Bmm-Xmm。
S30, carries out correlation analysis by the first frequency domain sequence matrix and the second frequency domain sequence matrix, obtains transformer
Correlation matrix and amplitude coefficient matrix;
Above-mentioned steps can be respectively to the first frequency domain sequence matrix ZmmWith the second frequency domain sequence Matrix CmmCarry out correlation point
The data processings such as analysis, to determine the correlation matrix and amplitude coefficient matrix of transformer.Above-mentioned first frequency domain sequence matrix Zmm
With the second frequency domain sequence Matrix CmmFor line number and columns equal matrix, the first frequency domain sequence matrix ZmmWith the second frequency domain sequence
Column matrix CmmThe constituent element element of position identical one can be to should determine that phase a coefficient correlation k and amplitude coefficient f, according to each constituent element
The corresponding phase coefficient correlation k of element can build correlation matrix, and width can be built according to the corresponding amplitude coefficient f of each group element
Value coefficient matrix.
S40, calculates the ∞ norm value S of the correlation matrix, and the amplitude coefficient matrix ∞ norm values F;
S50, the winding state of transformer is detected according to the ∞ norm values S and ∞ norm value F.
Above-mentioned steps can according to above-mentioned ∞ norm values S and ∞ norm value F respectively residing scope determine transformer around
Group state;Specifically, if above-mentioned ∞ norm values S is more than the first correlation coefficient threshold, and ∞ norm values F is less than first amplitude coefficient
Threshold value, then can be determined that the winding state of transformer is good;If above-mentioned ∞ norm values S is in the second correlation coefficient threshold to the first phase
Close on coefficient threshold this numerical intervals, and ∞ norm values F first amplitude coefficient threshold to the second amplitude coefficient threshold value this
On numerical intervals, then can be determined that the winding state of transformer is slight loosening;If above-mentioned ∞ norm values S is in third phase relation number
On threshold value to the second correlation coefficient threshold this numerical intervals, and ∞ norm values F is in the second amplitude coefficient threshold value to the 3rd amplitude
On this numerical intervals of coefficient threshold, then it can be determined that the winding state of transformer loosens for severe;If above-mentioned ∞ norm values S exists
4th correlation coefficient threshold to third phase is closed on coefficient threshold this numerical intervals, and ∞ norm values F is in the 3rd amplitude coefficient threshold
On value to the 4th amplitude coefficient threshold value this numerical intervals, then the winding state that can be determined that transformer is torsional deformation;If on
∞ norm values S is stated less than the 4th correlation coefficient threshold, and ∞ norm values F is more than the 4th amplitude coefficient threshold value, then can be determined that change
The winding state of depressor is short-circuit condition;Above-mentioned first correlation coefficient threshold is more than the second correlation coefficient threshold, the second phase relation
Number threshold value is more than third phase and closes coefficient threshold, and third phase closes coefficient threshold and is more than the 4th correlation coefficient threshold, the first coefficient correlation
Threshold value, the second correlation coefficient threshold, third phase close coefficient threshold and the 4th correlation coefficient threshold respectively can be according to transformer
The corresponding square formation dimension m in tank surface vibration acquisition face determines that above-mentioned 4th amplitude coefficient threshold value is more than the 3rd amplitude coefficient threshold
Value, the 3rd amplitude coefficient threshold value is more than the second amplitude coefficient threshold value, and the second amplitude coefficient threshold value is more than first amplitude coefficient threshold,
First amplitude coefficient threshold, the second amplitude coefficient threshold value, the 3rd amplitude coefficient threshold value and the 4th amplitude coefficient threshold value respectively can be with
The corresponding square formation dimension m in tank surface vibration acquisition face according to transformer is determined.
Specifically, if the corresponding square formation dimension m=3 in tank surface vibration acquisition face of transformer, above-mentioned first is related
Coefficient threshold, which can take the 2.6, second correlation coefficient threshold that 2.4, third phase can be taken to close coefficient threshold, can take the 2.15, the 4th phase
Close coefficient threshold can take 1.9, first amplitude coefficient threshold can take the 0.26, second amplitude coefficient threshold value can take 0.45, the
Three amplitude coefficient threshold values can take the 0.6, the 4th amplitude coefficient threshold value to take 0.75.Above-mentioned steps S40 can be detected according to table 1
The winding state of transformer.
Table 1
The winding state detection method for the transformer that the present invention is provided, can gather transformer in the early stage in test process
Test vibration signal, generates the first frequency domain sequence matrix of transformer;Operation of the transformer in practical work process is gathered to shake
Dynamic signal, generates the second frequency domain sequence matrix of transformer;Again by the first frequency domain sequence matrix and the second frequency domain sequence square
Battle array carries out correlation analysis, obtains the correlation matrix and amplitude coefficient matrix of transformer;To calculate the coefficient correlation
The ∞ norm value S of matrix, and the amplitude coefficient matrix ∞ norm values F;According to the ∞ norm values S and ∞ norm value F
Detect the winding state of transformer;The winding state of above-mentioned transformer can be according to transformer factory data (the first frequency domain sequence
Column matrix) and the course of work in service data (the second frequency domain sequence matrix) accordingly detected, what detection process was covered
Data more fully, effectively increase the Detection results of deformation of transformer winding state-detection.
In one embodiment, it is above-mentioned that the first frequency domain sequence matrix and the second frequency domain sequence matrix are subjected to correlation
Analysis, the process of the correlation matrix and amplitude coefficient matrix that obtain transformer can include:
According to the arrangement position of each element in the first frequency domain sequence matrix and the second frequency domain sequence matrix, set up described
In first frequency domain sequence matrix and the second frequency domain sequence matrix, the corresponding relation between the identical element of position obtains multigroup right
Answer element;
Phase the coefficient correlation k and amplitude coefficient f between each group corresponding element are calculated respectively;
Correlation matrix is built according to each phase coefficient correlation k, amplitude coefficient square is built according to each amplitude coefficient f
Battle array.
Above-mentioned first frequency domain sequence matrix ZmmFor m row m column matrix, the second frequency domain sequence Matrix Cmm ZmmFor m row m row squares
Battle array, the first frequency domain sequence matrix ZmmWith the second frequency domain sequence Matrix CmmThe constituent element of position identical one element is one group of corresponding element, respectively
Group corresponding element can be to should determine that an a phase coefficient correlation k and amplitude coefficient f.
As one embodiment, the above-mentioned process for calculating the phase coefficient correlation k between each group corresponding element respectively can be wrapped
Include:
In formula, N (i) represents numerical value of i-th of element in the first frequency domain sequence matrix or the second frequency domain sequence matrix
Sequence length, specifically, above-mentioned sequence of values length can be i-th of element in the first frequency domain sequence matrix or the second frequency domain
In sequence matrix data sequence number (i.e. in the first frequency domain sequence matrix or the second frequency domain sequence matrix, the data of preceding i element
Length), p (i) represents i-th of element of the second frequency domain sequence matrix, and q (i) represents i-th yuan of the first frequency domain sequence matrix
Element, covpqFor the covariance between i-th group of corresponding element, DpFor the corresponding variance of second i-th of element of frequency domain sequence matrix, Dq
For the corresponding variance of first i-th of element of frequency domain sequence matrix.
As one embodiment, the above-mentioned process for calculating the amplitude coefficient f between each group corresponding element respectively can include:
In formula, N (i) represents numerical value of i-th of element in the first frequency domain sequence matrix or the second frequency domain sequence matrix
Sequence length, p (i) represents i-th of element of the second frequency domain sequence matrix, and q (i) represents i-th of the first frequency domain sequence matrix
Element.
The winding state detection method for the transformer that the present embodiment is provided, calculates the phase phase between each group corresponding element respectively
Coefficient k and amplitude coefficient f are closed, is the correlation matrix according to constructed by each phase coefficient correlation k, according to each amplitude system
Amplitude coefficient matrix constructed by number f is more accurate, so as to ensure that the accuracy of the winding state of detection transformer.
With reference to shown in Fig. 2, Fig. 2 is the winding state detecting system structural representation of the transformer of one embodiment, including:
First acquisition module 10, for gathering test vibration signal of the transformer in the early stage in test process, according to described
First frequency domain sequence matrix of transformer described in test vibration signal generation;
Second acquisition module 20, for gathering operational shock signal of the transformer in practical work process, according to described
Second frequency domain sequence matrix of transformer described in operational shock signal generation;
Analysis module 30, for the first frequency domain sequence matrix and the second frequency domain sequence matrix to be carried out into correlation point
Analysis, obtains the correlation matrix and amplitude coefficient matrix of transformer;
Computing module 40, the ∞ norm value S for calculating the correlation matrix, and the amplitude coefficient matrix
∞ norm values F;
Detection module 50, the winding state for detecting transformer according to the ∞ norm values S and ∞ norm value F.
The winding state for the transformer that the winding state detecting system for the transformer that the present invention is provided is provided with the present invention is examined
Survey method is corresponded, the technical characteristic illustrated in the embodiment of the winding state detection method of the transformer and its beneficial effect
Really suitable for the embodiment of the winding state detecting system of transformer, hereby give notice that.
Based on example as described above, a kind of computer-readable recording medium is also provided in one embodiment, stored thereon
There is computer program, the computer program realizes the winding state detection side of transformer as described above when being executed by processor
Method.
One of ordinary skill in the art will appreciate that realize all or part of flow in above-described embodiment method, being can be with
The hardware of correlation is instructed to complete by computer program, described computer program can be stored in a non-volatile calculating
In machine read/write memory medium, in such as embodiment of the present invention, the program can be stored in the storage medium of computer system, and by
At least one computing device in the computer system, to realize the flow for including the embodiment such as above-mentioned each method.Wherein,
Described storage medium can be magnetic disc, CD, read-only memory (Read-Only Memory, ROM) or random storage note
Recall body (Random Access Memory, RAM) etc..
Based on example as described above, with reference to shown in Fig. 3, the present invention also provides a kind of computer equipment 60, the computer
Equipment includes memory 61, processor 62 and is stored in the computer program that can be run on memory 62 and on processor 61,
The processor 61 realizes the winding state detection of any one transformer in each embodiment as described above when performing described program
Method.
Above computer equipment 60 can including computer etc. Intelligent treatment equipment.One of ordinary skill in the art will appreciate that
The computer program that memory 61 is stored, it is corresponding with the description in the winding state detection method embodiment of above-mentioned transformer,
Processor 62 can also be used to perform other executable instructions that memory 61 is stored.
Each technical characteristic of embodiment described above can be combined arbitrarily, to make description succinct, not to above-mentioned reality
Apply all possible combination of each technical characteristic in example to be all described, as long as however, the combination of these technical characteristics is not deposited
In contradiction, the scope of this specification record is all considered to be.
Embodiment described above only expresses the several embodiments of the present invention, and it describes more specific and detailed, but simultaneously
Can not therefore it be construed as limiting the scope of the patent.It should be pointed out that coming for one of ordinary skill in the art
Say, without departing from the inventive concept of the premise, various modifications and improvements can be made, these belong to the protection of the present invention
Scope.Therefore, the protection domain of patent of the present invention should be determined by the appended claims.
Claims (10)
1. the winding state detection method of a kind of transformer, it is characterised in that comprise the following steps:
Gather test vibration signal of the transformer in the early stage in test process, the transformation according to the test vibration signal generation
First frequency domain sequence matrix of device;
Gather operational shock signal of the transformer in practical work process, the transformation according to the operational shock signal generation
Second frequency domain sequence matrix of device;
The first frequency domain sequence matrix and the second frequency domain sequence matrix are subjected to correlation analysis, the phase relation of transformer is obtained
Matrix number and amplitude coefficient matrix;
Calculate the ∞ norm value S of the correlation matrix, and the amplitude coefficient matrix ∞ norm values F;
The winding state of transformer is detected according to the ∞ norm values S and ∞ norm value F.
2. the winding state detection method of transformer according to claim 1, it is characterised in that the test process at initial stage
Including delivery test process or initial operation stage running.
3. the winding state detection method of transformer according to claim 2, it is characterised in that the collection transformer exists
Test vibration signal in initial stage test process, the first frequency domain sequence of transformer according to the test vibration signal generation
The process of matrix includes:
In the early stage in test process, the multiple spot vibration signal of collection transformer tank surface in the unloaded state obtains original sky
Carry vibration signal sequence matrix Xmm, the multiple spot vibration signal of transformer tank surface under rated load condition is gathered, original is obtained
Begin load vibration signal frequency domain sequence matrix Ymm;
According to the original unloaded vibration signal sequence matrix XmmWith original load vibration signal frequency domain sequence matrix YmmCalculate the
One frequency domain sequence matrix Zmm;Wherein, Zmm=Ymm-Xmm。
4. the winding state detection method of transformer according to claim 1, it is characterised in that the collection transformer exists
Operational shock signal in practical work process, the second frequency domain sequence of transformer according to the operational shock signal generation
The process of matrix includes:
In practical work process, the multiple spot vibration signal of collection transformer tank surface in the unloaded state obtains sky to be measured
Carry vibration signal sequence matrix Amm, the multiple spot vibration signal of transformer tank surface under rated load condition is gathered, is treated
Survey load vibration signal frequency domain sequence matrix Bmm;
According to the original unloaded vibration signal sequence matrix AmmWith original load vibration signal frequency domain sequence matrix BmmCalculate the
Two frequency domain sequence Matrix Csmm;Wherein, Cmm=Bmm-Amm。
5. the winding state detection method of transformer according to claim 1, it is characterised in that described by the described first frequency
Domain sequence matrix and the second frequency domain sequence matrix carry out correlation analysis, obtain the correlation matrix and amplitude coefficient of transformer
The process of matrix includes:
According to the arrangement position of each element in the first frequency domain sequence matrix and the second frequency domain sequence matrix, described first is set up
In frequency domain sequence matrix and the second frequency domain sequence matrix, the corresponding relation between the identical element of position obtains multigroup corresponding element
Element;
Phase the coefficient correlation k and amplitude coefficient f between each group corresponding element are calculated respectively;
Correlation matrix is built according to each phase coefficient correlation k, amplitude coefficient matrix is built according to each amplitude coefficient f.
6. the winding state detection method of transformer according to claim 5, it is characterised in that described to calculate each group respectively
The process of phase coefficient correlation k between corresponding element includes:
<mrow>
<mi>k</mi>
<mo>=</mo>
<mfrac>
<msub>
<mi>cov</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mrow>
<msqrt>
<msub>
<mi>D</mi>
<mi>p</mi>
</msub>
</msqrt>
<msqrt>
<msub>
<mi>D</mi>
<mi>q</mi>
</msub>
</msqrt>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>cov</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>&lsqb;</mo>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>&lsqb;</mo>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>D</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>D</mi>
<mi>q</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>;</mo>
</mrow>
In formula, N (i) represents sequence of values of i-th of element in the first frequency domain sequence matrix or the second frequency domain sequence matrix
Length, p (i) represents i-th of element of the second frequency domain sequence matrix, and q (i) represents i-th of element of the first frequency domain sequence matrix,
covpqFor the covariance between i-th group of corresponding element, DpFor the corresponding variance of second i-th of element of frequency domain sequence matrix, DqFor
The corresponding variance of first i-th of element of frequency domain sequence matrix.
7. the winding state detection method of transformer according to claim 5, it is characterised in that described to calculate each group respectively
The process of amplitude coefficient f between corresponding element includes:
<mrow>
<mi>f</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mo>&lsqb;</mo>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</msubsup>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
In formula, N (i) represents sequence of values of i-th of element in the first frequency domain sequence matrix or the second frequency domain sequence matrix
Length, p (i) represents i-th of element of the second frequency domain sequence matrix, and q (i) represents i-th of element of the first frequency domain sequence matrix.
8. a kind of winding state detecting system of transformer, it is characterised in that including:
First acquisition module, for gathering test vibration signal of the transformer in the early stage in test process, shakes according to the test
First frequency domain sequence matrix of transformer described in dynamic signal generation;
Second acquisition module, for gathering operational shock signal of the transformer in practical work process, shakes according to the operation
Second frequency domain sequence matrix of transformer described in dynamic signal generation;
Analysis module, for the first frequency domain sequence matrix and the second frequency domain sequence matrix to be carried out into correlation analysis, is obtained
The correlation matrix and amplitude coefficient matrix of transformer;
The ∞ norms of computing module, the ∞ norm value S for calculating the correlation matrix, and the amplitude coefficient matrix
Value F;
Detection module, the winding state for detecting transformer according to the ∞ norm values S and ∞ norm value F.
9. a kind of computer-readable recording medium, is stored thereon with computer program, it is characterised in that the computer program is located
Manage the winding state detection method that the transformer as described in claim 1 to 7 any one is realized when device is performed.
10. a kind of computer equipment, including memory, processor and storage are on a memory and the meter that can run on a processor
Calculation machine program, it is characterised in that realized during the computing device described program as described in claim 1 to 7 any one
The winding state detection method of transformer.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710400131.2A CN107132451A (en) | 2017-05-31 | 2017-05-31 | The winding state detection method and system of transformer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710400131.2A CN107132451A (en) | 2017-05-31 | 2017-05-31 | The winding state detection method and system of transformer |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107132451A true CN107132451A (en) | 2017-09-05 |
Family
ID=59735036
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710400131.2A Pending CN107132451A (en) | 2017-05-31 | 2017-05-31 | The winding state detection method and system of transformer |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107132451A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110286289A (en) * | 2019-06-30 | 2019-09-27 | 广东石油化工学院 | A kind of running state of transformer vibration sound detection signal filtering method and system using low-rank matrix recovery |
CN110702043A (en) * | 2019-10-24 | 2020-01-17 | 长春工程学院 | Power transformer winding deformation fault detection method |
CN111006579A (en) * | 2019-12-27 | 2020-04-14 | 广东电网有限责任公司电力科学研究院 | Transformer online winding deformation diagnosis method, system and equipment |
CN111521259A (en) * | 2020-04-30 | 2020-08-11 | 中国恩菲工程技术有限公司 | Grinding machine detection method, device and equipment |
CN112034433A (en) * | 2020-07-09 | 2020-12-04 | 重庆邮电大学 | Through-wall passive moving target detection method based on interference signal reconstruction |
CN113701684A (en) * | 2021-08-05 | 2021-11-26 | 西安交通大学 | Transformer winding state detection method, device, equipment and storage medium |
CN113739730A (en) * | 2021-08-30 | 2021-12-03 | 西安交通大学 | Transient acoustic signal-based transformer winding deformation detection method and system |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101701995A (en) * | 2009-11-12 | 2010-05-05 | 重庆大学 | Impulse response analytical test apparatus and method for detecting deformation of transformer winding |
CN103513138A (en) * | 2013-10-15 | 2014-01-15 | 国家电网公司 | Transformer winding and iron core fault diagnosis method based on vibration characteristic |
-
2017
- 2017-05-31 CN CN201710400131.2A patent/CN107132451A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101701995A (en) * | 2009-11-12 | 2010-05-05 | 重庆大学 | Impulse response analytical test apparatus and method for detecting deformation of transformer winding |
CN103513138A (en) * | 2013-10-15 | 2014-01-15 | 国家电网公司 | Transformer winding and iron core fault diagnosis method based on vibration characteristic |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110286289A (en) * | 2019-06-30 | 2019-09-27 | 广东石油化工学院 | A kind of running state of transformer vibration sound detection signal filtering method and system using low-rank matrix recovery |
CN110286289B (en) * | 2019-06-30 | 2021-06-11 | 广东石油化工学院 | Filtering method for vibration and sound detection signal of transformer |
CN110702043A (en) * | 2019-10-24 | 2020-01-17 | 长春工程学院 | Power transformer winding deformation fault detection method |
CN110702043B (en) * | 2019-10-24 | 2021-05-11 | 长春工程学院 | Power transformer winding deformation fault detection method |
CN111006579A (en) * | 2019-12-27 | 2020-04-14 | 广东电网有限责任公司电力科学研究院 | Transformer online winding deformation diagnosis method, system and equipment |
CN111521259A (en) * | 2020-04-30 | 2020-08-11 | 中国恩菲工程技术有限公司 | Grinding machine detection method, device and equipment |
CN111521259B (en) * | 2020-04-30 | 2022-02-18 | 中国恩菲工程技术有限公司 | Grinding machine detection method, device and equipment |
CN112034433A (en) * | 2020-07-09 | 2020-12-04 | 重庆邮电大学 | Through-wall passive moving target detection method based on interference signal reconstruction |
CN112034433B (en) * | 2020-07-09 | 2024-01-12 | 深圳市领冠检测技术有限公司 | Through-wall passive moving target detection method based on interference signal reconstruction |
CN113701684A (en) * | 2021-08-05 | 2021-11-26 | 西安交通大学 | Transformer winding state detection method, device, equipment and storage medium |
CN113739730A (en) * | 2021-08-30 | 2021-12-03 | 西安交通大学 | Transient acoustic signal-based transformer winding deformation detection method and system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107132451A (en) | The winding state detection method and system of transformer | |
CN111487046B (en) | Fault diagnosis method for circuit breaker voiceprint and vibration entropy feature fusion | |
CN106443316B (en) | Multi-information detection method and device for deformation state of power transformer winding | |
CN102998545B (en) | Online transformer winding work state monitoring method | |
CN103575523B (en) | The rotary machinery fault diagnosis method of kurtosis-envelope spectrum analysis is composed based on FastICA- | |
CN103439631B (en) | Ground net corrosion condition detection method and system | |
CN105064420A (en) | High-pile wharf foundation pile damage diagnosis method based on structural residual modal force | |
CN103822696B (en) | A kind of diagnostic method of transformer winding state | |
CN104897780B (en) | A kind of method positioned using Acoustic Emission Signal Energy to acoustic emission source | |
CN104849050A (en) | Rolling bearing fault diagnosis method based on composite multi-scale permutation entropies | |
CN106338705A (en) | Error verifying method and device for three-phase electric energy metering device | |
CN110007141B (en) | Resonance point detection method based on voltage and current harmonic similarity | |
CN110647924B (en) | GIS equipment state evaluation method based on support vector description and K-nearest neighbor algorithm | |
CN103968939B (en) | Based on the Transformer Winding looseness fault detection method of average displacement method | |
CN104112239A (en) | Method and device for transformer state assessment based on benchmark state analysis | |
CN109657720B (en) | On-line diagnosis method for turn-to-turn short circuit fault of power transformer | |
CN103645248A (en) | High-temperature alloy grain size evaluation method based on ultrasonic phase velocity | |
US20200393347A1 (en) | Imaging Method of Internal Defects in Longitudinal Sections of Trees | |
CN109932053A (en) | A kind of state monitoring apparatus and method for high-voltage shunt reactor | |
CN112415325A (en) | Single-phase earth fault section identification method based on edge calculation | |
CN117421976A (en) | Wind power plant current collection line fault positioning monitoring method and system based on ensemble learning | |
CN106706285B (en) | A kind of brake disc intrinsic frequency online test method | |
CN110098610B (en) | Real-time identification method and system for oscillation leading mode of power system under fault disturbance | |
CN106908688A (en) | A kind of portable power transformer winding failure real-time diagnosis equipment | |
CN110411721A (en) | A kind of marine riser damage positioning method and system |
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: 20170905 |
|
RJ01 | Rejection of invention patent application after publication |