CN107328467A - A kind of Transformer Winding thrust change detecting method based on recurrence quantification analysis - Google Patents
A kind of Transformer Winding thrust change detecting method based on recurrence quantification analysis Download PDFInfo
- Publication number
- CN107328467A CN107328467A CN201710617391.5A CN201710617391A CN107328467A CN 107328467 A CN107328467 A CN 107328467A CN 201710617391 A CN201710617391 A CN 201710617391A CN 107328467 A CN107328467 A CN 107328467A
- Authority
- CN
- China
- Prior art keywords
- mrow
- rqa
- winding
- thrust
- 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.)
- Pending
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
Abstract
The present invention discloses a kind of Winding in Power Transformer thrust change detecting method based on recurrence quantification analysis, and method gathers the vibration signal of tank surface to the normal transformer of winding thrust first by short circuit experiment method, uses recurrence quantification analysis(RQA)Method obtains referring to RQA parameters.Then short circuit experiment method vasculum body surface vibration signals are used to the transformer of winding to be checked with identical sample frequency and sampling time, and uses recurrence quantification analysis(RQA)Obtain Transformer Winding to be measured and compress the unknown RQA parameters of degree.Again by RQA parameters with being compared with reference to RQA parameters, to determine whether Transformer Winding thrust changes.The present invention not only sets up the RQA measurements of the overall thrust variation tendency of detection winding, and also the vibration signal of single sensor is directly carried out RQA analyses and sets up local thrust variation tendency to measure.It is not only few, simple to operate using equipment, and diagnosis is accurate, being capable of situation that accurately detection winding loosens by thrust change.
Description
Technical field
It is particularly a kind of to be based on recurrence quantification the present invention relates to voltage transformer condition monitoring and fault diagnosis technical field
The Winding in Power Transformer thrust change detecting method of analysis.
Background technology
Power transformer is one of power system most critical equipment, is also that current intelligent grid construction will realize the energy with future
The nucleus equipment of Internet Strategy.Whether generate electricity, transmit electricity or distribution, be dependent on the normal table operation of transformer.Transformation
Device involves great expense, and volume is big, maintenance difficult, once breaking down, maintenance will take considerable time, man power and material, produces simultaneously
Raw great social influence.Therefore, it is necessary to be monitored in real time to transformer.According to statistics, high-power transformer is made because of external short circuit
The first place of Accident of Transformer is had been raised into mechanical breakdowns such as winding loosenings.Therefore, to the compression shape of Transformer Winding in operation
State is monitored, and winding looseness fault hidden danger is found early, is had important practical significance.
The vibration on power transformer body surface and the compression situation of Transformer Winding and iron core are closely related.Transformer sheet
The vibration of body is mainly derived from basket vibration caused by core vibration caused by magnetostriction and leakage field and electromagnetic force.Drawn by leakage field
The core vibration risen is more faint compared with other vibrations, and the influence of leakage field can be neglected.When secondary winding is short-circuit, applied voltage
Seldom, now core vibration caused by magnetostriction can be neglected magnetic flux, then pass through transformer oil and tank wall in very little, iron core
The vibration for being delivered to tank surface is mainly caused by winding.
The failure of mechanical structure or the evolution process of damage are usually nonlinear, and winding is no exception.Basket vibration is certainly
Body has more obvious nonlinear characteristic, along with vibration signal decay, phase caused by various factors during Vibration propagation
The change of shifting, the evolutionary process that its dynamics changes with thrust is complicated and changeable, and existing thrust change is brought
Magnetic fluxleakage distribution change, have the change of transformer rigidity that winding mechanical structure brings and damping again.In this case,
Go to measure this change often due to ignoring some factors therein using traditional linear method (such as Fourier analysis method)
And the effect that can not have been received.
The content of the invention
The technical problem to be solved in the present invention is:Using recurrence plot (RP) quantitative analysis method pair in nonlinear kinetics
Basket vibration signal is analyzed, to detect the change of Transformer Winding thrust.
Recurrence plot provides a width signal autocorrelative overall figure in various possible time scales, than traditional
Linear transformation method can preferably detect the change on signal structure.The characteristic quantity extracted from recurrence plot is dynamic to system architecture
The change of mechanics is more sensitiveer than conventional linear analysis method, therefore it is higher to use it to have the progress detection of winding thrust
Sensitivity.
The technical scheme that the present invention takes is:A kind of Transformer Winding thrust change detection based on recurrence quantification analysis
Method, including step:
S1, multiple vibrating sensor monitoring points, the vibrating sensor output of each monitoring point are set on transformer-cabinet surface
End connects data collecting instrument respectively;
S2, sets sample frequency and the sampling time of data collecting instrument, winding thrust is gathered using short circuit experiment method
The vibration signal of transformer tank surface when normal, be designated as { x (i) } (i=1,2 ..., n);
S3, calculates the Embedded dimensions m and time delay τ of vibration signal, vibration signal is temporally postponed and Embedded dimensions
Carry out phase space reconfiguration;
S4, calculates any two points X in phase space reconstructioniAnd XjThe distance between | | Xi-Xj| |, and chosen distance threshold value r,
Calculate recursion matrix Rij:
Rij=Θ (r- | | Xi-Xj||) (1)
Wherein, Θ () is Heaviside functions, if x >=0, Θ (x)=1, if x<0, then Θ (x)=0;
S5, using i as abscissa, recursion matrix R is drawn by ordinate of jij, recurrence plot is obtained, is then extracted in recurrence plot
RQA characteristic quantities;
RQA characteristic quantities based on single measuring point, set up local winding thrust RQA measurements;
RQA characteristic quantities based on multiple measuring points, set up the RQA measurements of the overall thrust of winding;
S6, with identical sample frequency and sampling time in S2, using short circuit experiment method collection winding thrust treat
The vibration signal of each monitoring point in transformer-cabinet surface of survey, the step of according to S3 to S5, obtains winding entirety and local compaction
The RQA measurements of power;
S7, the RQA measurements of the overall and local thrust of winding for the correspondence transformer to be measured that S6 is obtained, is obtained with S5
The RQA measurements of overall and local thrust are compared, and define RQA metric difference threshold values:
If the difference of the RQA metric values of the two overall thrust is less than RQA metric difference threshold values, winding to be checked is judged
Thrust do not change, winding does not produce loosening.
If the difference of the RQA metric values of the two overall thrust is more than or equal to RQA metric difference threshold values, it is determined as
Winding thrust changes, and winding produces loosening.
The present invention in use, by experience set different RQA metric differences threshold values can be used for judge winding normally, winding
The situation that not exclusively loosening and winding loosen completely.
It is preferred that, in step S7, the RQA metric differences threshold value includes entirety RQA metric difference threshold values, and correspondence is in
Between phase monitoring point local RQA metric differences threshold value and corresponding local RQA discrepancy thresholds away from interphase monitoring point;
In step S7, the RQA measurements of the overall and local thrust of winding for the correspondence transformer to be measured that S6 is obtained, with S5
The RQA measurements of obtained entirety and local thrust are compared, and then by the difference of the RQA metric values of the two overall thrust
It is different to be compared with overall RQA metric differences threshold value.The interphase is B phases.
Further, step S7 of the present invention also includes, and exceedes part RQA metric differences according to local RQA metric values difference
The monitoring location of different threshold value, judges the position that Transformer Winding to be measured loosens.Even correspond to the local RQA measurements of certain monitoring point
Numerical value differs greatly, then the corresponding winding position of the monitoring location may be the position loosened.
Short-circuit test method described in step S2 and step S6 of the present invention is, by the low pressure winding short circuit of transformer, in high pressure
Winding applies voltage and causes the short circuit current flow of low pressure winding to reach rated current.
It is preferred that, in S2, the vibration signal that data collecting instrument is gathered at each monitoring location at least continuous three times.Every time
Vibration signal is intercepted by the sample frequency and sampling time length of setting complete cycle.
In step S3, G-P algorithms, Cao algorithms etc. can be used by calculating Embedded dimensions m, and calculating time delay τ can be used
Autocorrelation Function method, average displacement method and mutual information method etc..
It is preferred that, in S3, using delay coordinate method, to vibration signal { x (i) }, (i=1,2 ..., n) carry out phase space weight
Structure, reconstruction signal is:
X (i)=x (i), x (i+ τ) ..., x (i+ (m-1) τ) } (1)
Wherein, i=1,2 ..., N;N=n- (m-1) τ, i are the mutually points of reconstruction attractor.
In step S4,1- norms, 2- norms and ∞-model can be used by calculating the distance in phase space reconstruction between any two points
Number etc.;Chosen distance threshold value r uses principle of experience, general to press 15% selection for being less than data standard difference.
It is preferred that, in S5, the RQA characteristic quantities of recurrence plot include recurrence rate RR, determine rate DET, average diagonal line length L and
Entropy ENTR:
Recurrence rate RR represents the ratio shared by recursive point sum in recurrence plot, is:
The rate DET of determination is to constitute the ratio counted parallel to the recurrence points of recurrence plot leading diagonal with total recurrence:
In formula, P (r, l) is the diagonal number that length is l in recurrence plot diagonal arrangement, lminIt is the length that diagonal is taken
Spend initial value.It is preferred that, take lmin=2.
Average diagonal line length L is the average value of catercorner length, is:
In formula,For total diagonal hop count;
Entropy ENTR is:
In formula,Represent the probability that length occurs for l diagonal.
It is preferred that, the overall RQA metric differences threshold value is (RR, DET, L, ENTR)=(2.64 × 10-4, 2.87 × 10-2, 5.35 × 10-1, 2.18 × 10-1), correspondence close to interphase monitoring point local RQA metric differences threshold value for (RR, DET, L,
ENTR)=(1.81 × 10-4, 2.41 × 10-2, 2.54 × 10-1, 2.33 × 10-2), part of the correspondence away from interphase monitoring point
RQA discrepancy thresholds are (RR, DET, L, ENTR)=(1.76 × 10-3, 1.33 × 10-1, 2.46,1.35).Winding is local and overall
The RQA metric difference threshold values of thrust are different.
Beneficial effect
The Multi-channel Vibration Signals of synchronized sampling are included analysis by the present invention, set up the overall thrust variation tendency of detection winding
RQA measurement.Meanwhile, also the vibration signal of single sensor is directly carried out RQA analyses and sets up local thrust to change
Trend is measured.Multicomponent signal RQA measurements can characterize the winding thrust change of the overall situation, and whether detection winding occurs on the whole
Loosen.Unitary signal RQA measurements then show the winding local compaction power change related to sensing station, when multicomponent signal
RQA measurements change, when display winding produces loosening, and unitary signal RQA measurements then can provide beneficial for fault location
Use for reference, judge specific winding failure phase.The present invention relates to equipment it is few, simple to operate, method uses non-linear dynamic
Recurrence plot (RP) quantitative analysis method in is analyzed winding vibration signal, to detect the change of Transformer Winding thrust
The method of change, the change detection sensitivity for inside transformer winding impaction state is higher, and testing result is also more accurate.
Brief description of the drawings
Fig. 1 show the flow chart of the present invention.
Fig. 2 show vibration monitoring point layout drawing.
Fig. 3 show the recurrence plot of vibration signal under winding thrust normal condition;
Fig. 4 show the recurrence plot of vibration signal under winding not exclusively loosening state;
Fig. 5 show the recurrence plot of vibration signal under winding loosening state completely.
Embodiment
Further described below in conjunction with the drawings and specific embodiments.
A kind of Transformer Winding thrust change detecting method based on recurrence quantification analysis, including step:
S1, multiple vibrating sensor monitoring points, the vibrating sensor output of each monitoring point are set on transformer-cabinet surface
End connects data collecting instrument respectively;
S2, sets sample frequency and the sampling time of data collecting instrument, winding thrust is gathered using short circuit experiment method
The vibration signal of transformer tank surface when normal, be designated as { x (i) } (i=1,2 ..., n);
S3, calculates the Embedded dimensions m and time delay τ of vibration signal, vibration signal is temporally postponed and Embedded dimensions
Carry out phase space reconfiguration;
S4, calculates any two points X in phase space reconstructioniAnd XjThe distance between | | Xi-Xj| |, and chosen distance threshold value r,
Calculate recursion matrix Rij:
Rij=Θ (r- | | Xi-Xj||) (1)
Wherein, Θ () is Heaviside functions, if x >=0, Θ (x)=1, if x<0, then Θ (x)=0;
S5, using i as abscissa, recursion matrix R is drawn by ordinate of jij, recurrence plot is obtained, is then extracted in recurrence plot
RQA characteristic quantities;
RQA characteristic quantities based on single measuring point, set up local winding thrust RQA measurements;
RQA characteristic quantities based on multiple measuring points, set up the RQA measurements of the overall thrust of winding;
S6, with identical sample frequency and sampling time in S2, using short circuit experiment method collection winding thrust treat
The vibration signal of each monitoring point in transformer-cabinet surface of survey, the step of according to S3 to S5, obtains winding entirety and local compaction
The RQA measurements of power;
S7, the RQA measurements of the overall and local thrust of winding for the correspondence transformer to be measured that S6 is obtained, is obtained with S5
Entirety and the RQA measurements of local thrust be compared, define RQA metric difference threshold values:
If the difference of the RQA metric values of the two overall thrust is less than RQA metric difference threshold values, winding to be checked is judged
Thrust do not change, winding does not produce loosening.
If the difference of the RQA metric values of the two overall thrust is more than or equal to RQA metric difference threshold values, it is determined as
Winding thrust changes, and winding produces loosening.
Embodiment 1
In order to weigh the compression degree that Transformer Winding is overall, the present embodiment includes the Multi-channel Vibration Signals of synchronized sampling
Analysis, sets up the RQA measurements of the overall thrust variation tendency of detection winding.Meanwhile, it is also straight to the vibration signal of single sensor
Row RQA is tapped into analyze and set up local thrust variation tendency measurement.The present invention is not only few, simple to operate using equipment, and
, being capable of situation that accurately detection winding loosens by thrust change and diagnosis is accurate.
Step S7 also includes, the monitoring location according to local RQA metric values difference more than RQA metric difference threshold values,
Judge the position that Transformer Winding to be measured loosens.The local RQA metric values for even corresponding to certain monitoring point differ greatly, then the prison
The corresponding winding position of point position may be the position loosened.
Short-circuit test method described in step S2 and step S6 is, by the low pressure winding short circuit of transformer, to be applied in high pressure winding
Making alive causes the short circuit current flow of low pressure winding to reach rated current.
In S2, the vibration signal that data collecting instrument is gathered at each monitoring location at least continuous three times, every time by setting
Sample frequency and sampling time length intercept vibration signal complete cycle.
In S3, G-P algorithms, Cao algorithms etc. can be used by calculating Embedded dimensions m, and calculating time delay τ can be used from pass
Join function method, average displacement method and mutual information method etc..
In S3, using delay coordinate method, to vibration signal { x (i) }, (i=1,2 ... n) carry out phase space reconfiguration, reconstruct letter
Number it is:
X (i)=x (i), x (i+ τ) ..., x (i+ (m-1) τ) } (1)
Wherein, i=1,2 ..., N;N=n- (m-1) τ, i are the mutually points of reconstruction attractor.
In S4,1- norms, 2- norms and ∞-norm can be used by calculating the distance in phase space reconstruction between any two points
Deng;Chosen distance threshold value r uses principle of experience, general to press 15% selection for being less than data standard difference.
In S5, the RQA characteristic quantities of recurrence plot include recurrence rate RR, determine rate DET, average diagonal line length L and entropy ENTR:
Recurrence rate RR represents the ratio shared by recursive point sum in recurrence plot, is:
The rate DET of determination is to constitute the ratio counted parallel to the recurrence points of recurrence plot leading diagonal with total recurrence:
In formula, P (r, l) is the diagonal number that length is l in recurrence plot diagonal arrangement, lminIt is the length that diagonal is taken
Spend initial value.It is preferred that, take lmin=2.
Average diagonal line length L is the average value of catercorner length, is:
In formula,For total diagonal hop count;
Entropy ENTR is:
In formula,Represent the probability that length occurs for l diagonal.
Overall RQA metric differences threshold value is (RR, DET, L, ENTR)=(2.64 × 10-4, 2.87 × 10-2, 5.35 × 10-1, 2.18 × 10-1), correspondence close to interphase monitoring point local RQA metric differences threshold value for (RR, DET, L, ENTR)=
(1.81×10-4, 2.41 × 10-2, 2.54 × 10-1, 2.33 × 10-2), local RQA difference of the correspondence away from interphase monitoring point
Threshold value is (RR, DET, L, ENTR)=(1.76 × 10-3, 1.33 × 10-1, 2.46,1.35).The local and overall thrust of winding
RQA metric difference threshold values be different.
Embodiment 2
The present embodiment Nanjing is set up one's own business transformer company a model SFZ10-31500/110 oil-immersed type power become
Depressor carries out winding thrust setting, and couple group is marked as YNd11, and low-pressure side rated voltage is 10.5kV, specified electricity
Flow for 1732A.
Vibration signal is carried out using model JF2020 vibration acceleration sensor and Nicolet data collecting instruments to adopt
Collection.3 measuring points are placed on transformer-cabinet surface, position is as shown in Figure 2.In view of the frequency range of transformer vibration signal,
Sample frequency is set to 10kHz during test.
During test, artificially B phase winding thrusts are configured by hydraulic system, are divided into normal (the specified pretension of winding
Power, 28MPa), not exclusively loosen (0.5 times of specified pretightning force, 14MPa) and completely loosening (pretightning force is zero) three kinds of situations, with
Realize the change of different winding machine performances.
First by the short circuit of experimental transformer low-pressure side three-phase windings, applied voltage is adjusted by pressure regulator in high-pressure side, made low
Side short circuit current flow is pressed close to rated current, high current situation during the specified operation of analogue transformer.When short circuit current flow reach it is specified
During electric current, the vibration signal on oil tank of transformer surface is measured.
The present embodiment selection is analyzed vibration signal close to No. 2 measuring points of B phase windings, after noise reduction process, No. 2 surveys
Vibration signal time-frequency figure and recurrence plot of the point under different winding impaction states are as shown in Figure 3.From time-frequency figure, three kinds of situations
Under vibration signal amplitude it is essentially identical, difference less, not exclusively loosen when amplitude it is slightly larger compared with normal condition, completely loosen
When amplitude it is maximum.Therefore, vibration signal is relevant with winding impaction state, but is only difficult to differentiate between winding difference pressure from amplitude variations
Tight state.Based on this, vibration signal is analyzed from recurrence plot.Basket vibration signal is main it can be seen from recurrence plot
Relevant with diagonal structure, there is notable difference in the recurrence plot under three kinds of winding impaction states, with the reduction of winding thrust,
Recursive point is significantly reduced, and diagonal structure fades away, or even diagonal line segment deteriorates to the recursive point of discrete distribution.This show with
The reduction of winding thrust, winding produces loosening, and the dynamic behavior of system is gradually deviated from original normal condition.Therefore,
Can be with the change of diagonal structure in recurrence plot to can intuitively reflect the trend that winding thrust reduces.
Recurrence plot can only carry out qualitative analysis to the dynamics of winding system, for quantitative description, extract recurrence plot
RQA characteristic quantities carry out quantitative analysis.The RQA characteristic quantities for extracting 3 measuring points constitute entirety RQA measurements, and the RQA of each measuring point is special
The amount of levying directly constitutes local RQA measurements.Tables 1 and 2 sets forth overall and part RQA measurements under three kinds of winding impaction states
Value.
Entirety RQA metrics under 1 three kinds of winding impaction states of table
Thrust/MPa | Recurrence rate | Degree of certainty | Average diagonal line length | Entropy |
28 | 0.007645 | 0.6327 | 8.5349 | 0.6781 |
14 | 0.005739 | 0.5043 | 5.6008 | 2.5574 |
0 | 0.003471 | 0.2587 | 2.8637 | 5.1352 |
Part RQA metrics under 2 three kinds of winding impaction states of table
It can be obtained by table 1, with the decline of B phase winding thrusts, the recurrence rates of overall RQA measurements, degree of certainty and average pair
Diagonal length is all gradually reduced, and entropy gradually increases.RQA of 4 of description selection based on recurrence dot density or diagonal structure
Measurement can embody the variation tendency of thrust decline, can be used in the detection of winding thrust change.
It can be obtained from table 2, with the decline of thrust, the recurrence rate of No. 1 and No. 3 measuring point, degree of certainty and average diagonal
Line length slightly has reduction, and entropy slightly increases, but amplitude of variation is little;And the recurrence rates of No. 2 measuring points, degree of certainty and average diagonal
Line length is obviously reduced, and entropy is significantly increased.Because No. 2 measuring points abut B phases, maximum, No. 1 and 3 are influenceed by thrust change
Number measuring point is influenceed small away from B phases by thrust change.This shows that sensitivity of the different measuring points to structure change is different, close to event
The RQA measurement changes for hindering the measuring point of position become apparent from.This can provide beneficial reference for the positioning of failure, and this is also to count respectively
Calculate the meaning of global thrust RQA measurements and local thrust RQA measurements.
Described above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, without departing from the technical principles of the invention, some improvement and deformation can also be made, these improve and deformed
Also it should be regarded as protection scope of the present invention.
Claims (8)
1. a kind of Transformer Winding thrust change detecting method based on recurrence quantification analysis, it is characterized in that, including step:
S1, multiple vibrating sensor monitoring points, the vibrating sensor output end point of each monitoring point are set on transformer-cabinet surface
Data collecting instrument is not connected;
S2, sets sample frequency and the sampling time of data collecting instrument, and it is normal to gather winding thrust using short circuit experiment method
When transformer tank surface vibration signal, be designated as { x (i) } (i=1,2 ..., n);
S3, calculates the Embedded dimensions m and time delay τ of vibration signal, and vibration signal is temporally postponed and Embedded dimensions are carried out
Phase space reconfiguration;
S4, calculates any two points X in phase space reconstructioniAnd XjThe distance between | | Xi-Xj| |, and chosen distance threshold value r, calculating is passed
Return matrix Rij:
Rij=Θ (r- | | Xi-Xj||) (1)
Wherein, Θ () is Heaviside functions, if x >=0, Θ (x)=1, if x<0, then Θ (x)=0;
S5, using i as abscissa, recursion matrix R is drawn by ordinate of jij, recurrence plot is obtained, the RQA in recurrence plot is then extracted
Characteristic quantity;
RQA characteristic quantities based on single measuring point, set up local winding thrust RQA measurements;
RQA characteristic quantities based on multiple measuring points, set up the RQA measurements of the overall thrust of winding;
S6, with identical sample frequency and sampling time in S2, using short circuit experiment method gather winding thrust it is to be measured
The vibration signal of each monitoring point in transformer-cabinet surface, the step of according to S3 to S5, obtains the overall and local thrust of winding
RQA is measured;
S7, the RQA measurements of the overall and local thrust of the winding of the correspondence transformer to be measured that S6 is obtained, with S5 obtain it is whole
The RQA measurements of body and local thrust are compared, and define RQA metric difference threshold values:
If the difference of the RQA metric values of the two overall thrust is less than RQA metric difference threshold values, the pressure of winding to be checked is judged
Clamp force does not change, and winding does not produce loosening;
If the difference of the RQA metric values of the two overall thrust is more than or equal to RQA metric difference threshold values, it is determined as winding
Thrust changes, and winding produces loosening.
2. according to the method described in claim 1, it is characterized in that, in step S7, the RQA metric differences threshold value includes overall
RQA metric difference threshold values, correspondence is monitored close to the local RQA metric differences threshold value of interphase monitoring point and correspondingly away from interphase
The local RQA discrepancy thresholds of point;
In step S7, the RQA measurements of the overall and local thrust of winding for the correspondence transformer to be measured that S6 is obtained are obtained with S5
Entirety and the RQA measurements of local thrust be compared, and then by the difference of the RQA metric values of the two overall thrust with
Overall RQA metric differences threshold value is compared.
3. method according to claim 2, it is characterized in that, step S7 also includes, super according to local RQA metric values difference
The monitoring location of part RQA metric difference threshold values is crossed, the position that Transformer Winding to be measured loosens is judged.
4. according to the method described in claim 1, it is characterized in that, short-circuit test method described in step S2 and step S6 is to become
The low pressure winding short circuit of depressor, applies voltage in high pressure winding and causes the short circuit current flow of low pressure winding to reach rated current.
5. according to the method described in claim 1, it is characterized in that, in S2, data collecting instrument at least each monitoring of continuous three collections
Vibration signal at point position.
6. according to the method described in claim 1, it is characterized in that, in S3, using delay coordinate method to vibration signal { x (i) } (i
=1,2 ..., phase space reconfiguration n) is carried out, reconstruction signal is:
X (i)=x (i), x (i+ τ) ..., x (i+ (m-1) τ) } (1)
Wherein, i=1,2 ..., N;N=n- (m-1) τ, i are the mutually points of reconstruction attractor.
7. method according to claim 6, it is characterized in that, in S5, the RQA characteristic quantities of recurrence plot include recurrence rate RR, really
Determine rate DET, average diagonal line length L and entropy ENTR:
Recurrence rate RR represents the ratio shared by recursive point sum in recurrence plot, is:
<mrow>
<mi>R</mi>
<mi>R</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msup>
<mi>N</mi>
<mn>2</mn>
</msup>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>R</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
The rate DET of determination is to constitute the ratio counted parallel to the recurrence points of recurrence plot leading diagonal with total recurrence:
<mrow>
<mi>D</mi>
<mi>E</mi>
<mi>T</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<msub>
<mi>l</mi>
<mi>min</mi>
</msub>
</mrow>
<mi>N</mi>
</munderover>
<mi>l</mi>
<mi>P</mi>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>l</mi>
<mi>P</mi>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, P (r, l) is the diagonal number that length is l in recurrence plot diagonal arrangement, lminAt the beginning of being the length that diagonal is taken
Value;
Average diagonal line length L is the average value of catercorner length, is:
<mrow>
<mi>L</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<msub>
<mi>l</mi>
<mi>min</mi>
</msub>
</mrow>
<mi>N</mi>
</munderover>
<mi>l</mi>
<mi>P</mi>
<mo>(</mo>
<mrow>
<mi>r</mi>
<mo>,</mo>
<mi>l</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<msub>
<mi>l</mi>
<mi>min</mi>
</msub>
</mrow>
<mi>N</mi>
</munderover>
<mi>P</mi>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula,For total diagonal hop count;
Entropy ENTR is:
<mrow>
<mi>E</mi>
<mi>N</mi>
<mi>T</mi>
<mi>R</mi>
<mo>=</mo>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>l</mi>
<mo>=</mo>
<msub>
<mi>l</mi>
<mi>min</mi>
</msub>
</mrow>
<mi>N</mi>
</munderover>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mi>ln</mi>
<mi> </mi>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula,Represent the probability that length occurs for l diagonal.
8. according to the method described in claim 1, it is characterized in that, the overall RQA metric differences threshold value for (RR, DET, L,
ENTR)=(2.64 × 10-4, 2.87 × 10-2, 5.35 × 10-1, 2.18 × 10-1), part of the correspondence close to interphase monitoring point
RQA metric differences threshold value is (RR, DET, L, ENTR)=(1.81 × 10-4, 2.41 × 10-2, 2.54 × 10-1, 2.33 × 10-2),
Local RQA discrepancy threshold of the correspondence away from interphase monitoring point is (RR, DET, L, ENTR)=(1.76 × 10-3, 1.33 × 10-1, 2.46,1.35).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710617391.5A CN107328467A (en) | 2017-07-26 | 2017-07-26 | A kind of Transformer Winding thrust change detecting method based on recurrence quantification analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710617391.5A CN107328467A (en) | 2017-07-26 | 2017-07-26 | A kind of Transformer Winding thrust change detecting method based on recurrence quantification analysis |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107328467A true CN107328467A (en) | 2017-11-07 |
Family
ID=60200101
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710617391.5A Pending CN107328467A (en) | 2017-07-26 | 2017-07-26 | A kind of Transformer Winding thrust change detecting method based on recurrence quantification analysis |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107328467A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110119690A (en) * | 2019-04-19 | 2019-08-13 | 河海大学 | A kind of shunt reactor vibration sensing regional selection method based on CRP and RQA |
CN110703149A (en) * | 2019-10-02 | 2020-01-17 | 广东石油化工学院 | Method and system for detecting vibration and sound of running state of transformer by utilizing character spacing |
CN112798688A (en) * | 2021-01-08 | 2021-05-14 | 大连理工大学 | Method for improving detection sensitivity of critical refraction longitudinal wave of material damage |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020370A (en) * | 2014-05-22 | 2014-09-03 | 中国南方电网有限责任公司电网技术研究中心 | Transformer internal fault diagnosis method based on virtual parameter change monitoring |
CN104535865A (en) * | 2014-12-30 | 2015-04-22 | 西安工程大学 | Comprehensive diagnosing method for operation troubles of power transformer based on multiple parameters |
-
2017
- 2017-07-26 CN CN201710617391.5A patent/CN107328467A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020370A (en) * | 2014-05-22 | 2014-09-03 | 中国南方电网有限责任公司电网技术研究中心 | Transformer internal fault diagnosis method based on virtual parameter change monitoring |
CN104535865A (en) * | 2014-12-30 | 2015-04-22 | 西安工程大学 | Comprehensive diagnosing method for operation troubles of power transformer based on multiple parameters |
Non-Patent Citations (2)
Title |
---|
宋晓 等: "基于递归定量特征的变压器励磁涌流识别", 《中南大学学报(自然科学版)》 * |
郭洁 等: "交叉递归图在变压器铁芯压紧力变化检测中的应用", 《高电压技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110119690A (en) * | 2019-04-19 | 2019-08-13 | 河海大学 | A kind of shunt reactor vibration sensing regional selection method based on CRP and RQA |
CN110119690B (en) * | 2019-04-19 | 2022-08-16 | 河海大学 | Parallel reactor vibration sensitive area selection method based on CRP and RQA |
CN110703149A (en) * | 2019-10-02 | 2020-01-17 | 广东石油化工学院 | Method and system for detecting vibration and sound of running state of transformer by utilizing character spacing |
CN112798688A (en) * | 2021-01-08 | 2021-05-14 | 大连理工大学 | Method for improving detection sensitivity of critical refraction longitudinal wave of material damage |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9404957B2 (en) | Fault diagnosis and preliminary location system and method for transformer core looseness | |
CN103245832B (en) | Based on harmonic wave time-frequency characteristic method for parameter estimation and the analyser of quick S-transformation | |
CN101782426B (en) | Detection method of looseness fault vibration of power transformer winding | |
CN104483598B (en) | A kind of Cable fault examination analysis method | |
CN105092024B (en) | A kind of Winding in Power Transformer axial compression state on_line monitoring method and system | |
CN206114822U (en) | Many information detection means of power transformer winding deformation state | |
CN103968939B (en) | Based on the Transformer Winding looseness fault detection method of average displacement method | |
CN106443316A (en) | Power transformer winding deformation state multi-information detection method and device | |
CN107101714A (en) | A kind of transformer health Evaluation method based on multi-measuring point vibration signal characteristics | |
CN101201386B (en) | Method for locating parameter type fault of analogue integrated circuit | |
CN107328467A (en) | A kind of Transformer Winding thrust change detecting method based on recurrence quantification analysis | |
CN101907437A (en) | Wavelet difference algorithm-based cable fault localization method | |
CN110333474A (en) | A kind of method and system of on-line checking mutual inductor metering abnormality | |
CN104748838B (en) | Transformer Winding based on finite element analysis, which loosens, judges system and method | |
CN110991481B (en) | High-voltage shunt reactor internal loosening fault diagnosis method based on cross wavelet transformation | |
CN106441547A (en) | Transformer vibration monitoring method and apparatus | |
CN112485747A (en) | Error state evaluation method for capacitor voltage transformer | |
CN109375010A (en) | Winding in Power Transformer method for diagnosing faults based on vibration distribution characteristic model | |
CN102616156B (en) | Device and method for detecting and elevating contact network state based on electrical capacity analysis | |
CN105675994A (en) | Measuring method for equivalent system harmonic impedance of feeder line in power distribution network | |
CN106932184A (en) | A kind of Diagnosis Method of Transformer Faults based on improvement hierarchical clustering | |
CN109932053A (en) | A kind of state monitoring apparatus and method for high-voltage shunt reactor | |
CN103229062A (en) | Method for measurement of total harmonic distortion | |
CN105093059A (en) | Transformer winding work state detection method based on normalization singular spectrum entropies | |
CN110161375A (en) | A kind of HVDC transmission line computation model based on distributed resistance parameter |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171107 |