阅读说明：本技术 一种电力变压器绕组变形故障检测方法 (Power transformer winding deformation fault detection method ) 是由 孙宏斌 谢蓓敏 张轶珠 祝晓宏 宋丹 刘春� 韩秀峰 周会林 黄伟 宋威 张晋菁 于 2019-10-24 设计创作，主要内容包括：本发明提供一种电力变压器绕组变形故障新的检测方法；该方法建立映射算子可以将一定范围的传感器振动数据映射为一个矢量,并基于正常运行的数据构建矢量变化的容忍范围,并进一步实现电力变压器绕组变形故障检测。由于本发明描述的方法经过了非线性函数的映射,并基于的一定变化范围,所以本发明可以做到不与具体的数值绑定,能够达到更加稳定的预测效果。(The invention provides a new detection method for deformation faults of a power transformer winding; the method establishes a mapping operator, can map sensor vibration data in a certain range into a vector, establishes a tolerance range of vector change based on data in normal operation, and further realizes the detection of the deformation fault of the power transformer winding. Because the method described by the invention is mapped by the nonlinear function and is based on a certain change range, the method can not be bound with a specific numerical value, and can achieve a more stable prediction effect.)

1. A method for detecting deformation faults of a power transformer winding comprises the following steps:

s1, respectively installing three vibration sensors ZA, ZB and ZC on the phase A winding, the phase B winding and the phase C winding of the power transformer; the three sensors are used for collecting an A-phase vibration array VA, a B-phase vibration array VB and a C-phase vibration array VC of the power transformer in a normal operation state; obtaining an array Length variable Length;

s101, respectively installing three vibration sensors ZA, ZB and ZC on an A-phase winding, a B-phase winding and a C-phase winding of a power transformer;

s102, vibration signals collected by the three vibration sensors ZA, ZB and ZC are transmitted to a digital collection card, the collection card is provided with 10Hz frequency collection data, and the collected data are respectively stored in an A-phase vibration array VA, a B-phase vibration array VB and a C-phase vibration array VC; VA, VB and VC are arrays, and each element corresponds to a vibration signal acquired once; the collection process is continued for 20 days;

s103, an array Length variable Length is equal to the number of elements of an array VA;

s2, constructing a mapping operator KOperator, wherein the input is an A-phase mapping operator input array TA, a B-phase mapping operator input array TB and a C-phase mapping operator input array TC, each array comprises 100 elements, and the output is a 6-element mapping operator output array TO;

s201, constructing a mapping operator KOperator, wherein the input of the KOperator is an A-phase mapping operator input array TA, a B-phase mapping operator input array TB and a C-phase mapping operator input array TC, and each array comprises 100 elements;

s202, calculating a mean value of TA for the a-phase mean variable AVGTA, and calculating a standard deviation of TA for the a-phase standard deviation variable STDTA; calculating the average value of TB by using a phase B average value variable AVGTB, and calculating the standard deviation of TB by using a phase B standard deviation variable STDTB; calculating the mean value of TC for the C-phase mean variable AVGTC, and calculating the standard deviation of TC for the C-phase standard deviation variable STDTA;

s203, calculating a temporary storage variable TA of the phase A, namely (TA-AVGTA)/2 × STDTA; phase B temporary storage variable TB ═ (TB-AVGTB)/2 × STDTB; the temporary storage variable TC of C phase is (TC-AVGTC)/2 × STDTC;

s204, counting the number of elements with the phase A counting variable NTA being greater than 0 in TA; the number of elements of the B-phase counting variable NTB which is greater than 0 in TB; the number of elements with the phase C counting variable NTC ═ TC larger than 0;

s205, calculating a mapping operator first dimension variable M1 ═ tanh (Σ (TA)/100); mapping operator second dimension variable M2 ═ tanh (Σ (TB)/100); the mapping operator third-dimensional variable M3 ═ tanh (Σ (TC)/100); wherein tanh is the hyperbolic tangent function

S206, calculating a mapping operator fourth-dimensional variable M4 ═ tanh (Abs (M1-M2));

s207, calculating a fifth dimension variable M5 ═ tanh ((M1+ M2+ M3));

s208, calculating a sixth dimension variable M6 ═ tanh (NTC/(NTA + NTB)) × (M4+ M5) of the mapping operator;

s209, constructing a mapping operator output array TO ═ M1, M2, M3, M4, M5 and M6;

s210, taking the TO as the output of the KOperator;

s3, calculating VA, VB and VC by using KOperator to obtain a vibration mapping center vector TCenter, a maximum mapping distance YDist and a maximum offset distance PDist;

s301, setting initial values of TCenter, [0,0,0,0, 0], ydexit, [ 0] and PDist;

s302, initializing a center list variable TCenterList [ ];

s303, calculating a center list counter variable TCenterCount ═ 0;

s304, obtaining a random integer by the position counter POS, wherein the range of the random integer is 1-Length-100;

s305, TA is an element from POS to POS +100 in VA, TB is an element from POS to POS +100 in VB, and TC is an element from POS to POS +100 in VC;

s306, calculating mapping operator output variable MM by using mapping operator KOperator (TA, TB, TC);

s307, adding MM into tcentercist, TCenterCount ═ TCenterCount + 1;

s308, if TCentERCount <10000, go to S304, otherwise go to S309;

s309, for all elements of tcentelist, performing mean statistics on each dimension of the element;

s310, POS is equal to 0, and a central temporary storage variable PO is equal to TCenter;

s311, TA is an element from POS to POS +100 in VA, TB is an element from POS to POS +100 in VB, and TC is an element from POS to POS +100 in VC;

s312, the computation mapping operator recalculates the output variable PP, which is calculated by the mapping operator KOperator (TA, TB, TC);

s313, the first distance variable dd1 ═ PP-TCenter |, the second distance variable dd2 ═ PP-PO |;

where | | | represents l2norm (l2 norm) of the calculation vector;

s314, YDist is dd1 if dd1> YDist; if dd2> PDist then PDist is dd 2;

S315,PO＝PP；

S316,POS＝POS+200；

s317, if the POS is less than the Length-100, turning to S311, otherwise, turning to S318;

s318, outputting a vibration mapping center vector TCenter, a maximum mapping distance YDist and a maximum offset distance PDist;

s4, continuously acquiring 200 vibration data for the power transformer by using vibration sensors ZA, ZB and ZC, respectively storing the acquired data into current A-phase data CurrentZA, current B-phase data CurrentZB and current C-phase data CurrentZC, and outputting winding deformation fault detection results;

s401, continuously acquiring 200 vibration data of the power transformer by using vibration sensors ZA, ZB and ZC, and respectively storing the acquired data into CurrentZA, CurrentZB and CurrentZC;

s402, TA is an element cut from 1 to 100 in CurrentZA, TB is an element cut from 1 to 100 in CurrentZB, and TC is an element cut from 1 to 100 in CurrentZC;

s403, calculating a first variable CurrentP1 of the current operator output result by using a mapping operator kopersonator (TA, TB, TC);

s404, currently temporarily storing a first variable TempD1 ═ CurrentP1-TCenter |;

where | | | represents l2norm of the calculation vector;

s405, TA being an element from 101 to 200 in CurrentZA, TB being an element from 101 to 200 in CurrentZB, and TC being an element from 101 to 200 in CurrentZC;

s406, calculating a second variable CurrentP2 of the current operator output result by using a mapping operator kopersonator (TA, TB, TC);

s407, currently temporarily storing the second variable TempD2 ═ CurrentP1-CurrentP2 |;

where | | | represents l2norm of the calculation vector;

s408, distance variation range index decision is 0.5 × (YDist-TempD1)/YDist +0.5 × (PDist-TempD 2)/PDist;

s409, if decision >0, the winding deformation fault is not generated, and the operation is switched to S411, otherwise, the operation is switched to S410;

s410, outputting: if the winding deformation fault occurs, the step goes to S412;

s411, outputting: if no winding deformation fault occurs, go to S412;

and S412, finishing the winding deformation fault judgment process.

Technical Field

The invention relates to a method for detecting a winding deformation fault of a power transformer, provides a method for detecting a fault of a new power transformer for a power system, and belongs to the technical field of safety control management of power transformers.

Background

The transformer is an important element in a power system and plays an important role in safe and reliable operation of a power grid. Because the short-circuit resistance of the transformer is directly influenced by the winding characteristics of the transformer, the fault of winding deformation is found in advance, safety accidents can be effectively prevented from happening, and the maintenance cost of a power grid management enterprise is obviously reduced.

Two measures are generally adopted for quite deformation of a winding of a power transformer, one measure is that when the power transformer has obvious faults and operation problems, the winding is manually disassembled and checked whether the winding is deformed, and the method belongs to a passive mode and has high cost and time consumption. In another mode, a vibration sensor is arranged on a winding, deformation and non-deformation vibration data sets are constructed, and data are learned by using intelligent models such as a neural network and a decision tree, so that the capability of automatically predicting the deformation fault of the winding is obtained; the intelligent model and the specific vibration numerical value characteristics are bound with each other, and due to the fact that the transmission of vibration is related to factors in multiple aspects such as the accurate installation position of the sensor equipment, the temperature, the installation position of the transformer, the use mode and the like, when the factors change within a certain range, the intelligent model can generate prediction errors, so that the method has high requirements on the assembly and the operation of the equipment, the method is unstable in practical application, and a good application effect is difficult to obtain.

Therefore, it is necessary to provide a new method for detecting the deformation fault of the power transformer winding more stably and reliably.

Disclosure of Invention

Aiming at the problems in the prior art, the invention provides a method for detecting the deformation fault of a power transformer winding; the mapping operator established by the method can map the sensor vibration data in a certain range into a vector, and the tolerance range of vector change is established based on the data in normal operation, so that the deformation fault detection of the power transformer winding is further realized.

The invention relates to a method for detecting the deformation fault of a power transformer winding, which adopts the technical scheme as follows:

s1, respectively installing three vibration sensors ZA, ZB and ZC on the phase A winding, the phase B winding and the phase C winding of the power transformer; the three sensors are used for collecting an A-phase vibration array VA, a B-phase vibration array VB and a C-phase vibration array VC of the power transformer in a normal operation state; obtaining an array Length variable Length;

s101, respectively installing three vibration sensors ZA, ZB and ZC on an A-phase winding, a B-phase winding and a C-phase winding of a power transformer;

s102, vibration signals collected by the three vibration sensors ZA, ZB and ZC are transmitted to a digital collection card, the collection card is provided with 10Hz frequency collection data, and the collected data are respectively stored in an A-phase vibration array VA, a B-phase vibration array VB and a C-phase vibration array VC; VA, VB and VC are arrays, and each element corresponds to a vibration signal acquired once; the collection process is continued for 20 days;

s103, an array Length variable Length is equal to the number of elements of an array VA;

s2, constructing a mapping operator KOperator, wherein the input is an A-phase mapping operator input array TA, a B-phase mapping operator input array TB and a C-phase mapping operator input array TC, each array comprises 100 elements, and the output is a 6-element mapping operator output array TO;

s201, constructing a mapping operator KOperator, wherein the input of the KOperator is an A-phase mapping operator input array TA, a B-phase mapping operator input array TB and a C-phase mapping operator input array TC, and each array comprises 100 elements;

s202, calculating a mean value of TA for the a-phase mean variable AVGTA, and calculating a standard deviation of TA for the a-phase standard deviation variable STDTA; calculating the average value of TB by using a phase B average value variable AVGTB, and calculating the standard deviation of TB by using a phase B standard deviation variable STDTB; calculating the mean value of TC for the C-phase mean variable AVGTC, and calculating the standard deviation of TC for the C-phase standard deviation variable STDTA;

s203, calculating a temporary storage variable TA of the phase A, namely (TA-AVGTA)/2 × STDTA; phase B temporary storage variable TB ═ (TB-AVGTB)/2 × STDTB; the temporary storage variable TC of C phase is (TC-AVGTC)/2 × STDTC;

s204, counting the number of elements with the phase A counting variable NTA being greater than 0 in TA; the number of elements of the B-phase counting variable NTB which is greater than 0 in TB; the number of elements with the phase C counting variable NTC ═ TC larger than 0;

s205, calculating a mapping operator first dimension variable M1 ═ tanh (Σ (TA)/100); mapping operator second dimension variable M2 ═ tanh (Σ (TB)/100); the mapping operator third-dimensional variable M3 ═ tanh (Σ (TC)/100); wherein tanh is the hyperbolic tangent function

S206, calculating a mapping operator fourth-dimensional variable M4 ═ tanh (Abs (M1-M2));

s207, calculating a fifth dimension variable M5 ═ tanh ((M1+ M2+ M3));

s208, calculating a sixth dimension variable M6 ═ tanh (NTC/(NTA + NTB)) × (M4+ M5) of the mapping operator;

s209, constructing a mapping operator output array TO ═ M1, M2, M3, M4, M5 and M6;

s210, taking the TO as the output of the KOperator;

s3, calculating VA, VB and VC by using KOperator to obtain a vibration mapping center vector TCenter, a maximum mapping distance YDist and a maximum offset distance PDist;

s301, setting initial values of TCenter, [0,0,0,0, 0], ydexit, [ 0] and PDist;

s302, initializing a center list variable TCenterList [ ];

s303, calculating a center list counter variable TCenterCount ═ 0;

s304, obtaining a random integer by the position counter POS, wherein the range of the random integer is 1-Length-100;

s305, TA is an element from POS to POS +100 in VA, TB is an element from POS to POS +100 in VB, and TC is an element from POS to POS +100 in VC;

s306, calculating mapping operator output variable MM by using mapping operator KOperator (TA, TB, TC);

s307, adding MM into tcentercist, TCenterCount ═ TCenterCount + 1;

s308, if TCentERCount <10000, go to S304, otherwise go to S309;

s309, for all elements of tcentlist, performing mean statistics on each dimension of the element;

s310, POS is equal to 0, and a central temporary storage variable PO is equal to TCenter;

s311, TA is an element from POS to POS +100 in VA, TB is an element from POS to POS +100 in VB, and TC is an element from POS to POS +100 in VC;

s312, the computation mapping operator recalculates the output variable PP, which is calculated by the mapping operator KOperator (TA, TB, TC);

s313, the first distance variable dd1 ═ PP-TCenter |, the second distance variable dd2 ═ PP-PO |;

where | | | represents l2norm (l2 norm) of the calculation vector;

s314, YDist is dd1 if dd1> YDist; if dd2> PDist then PDist is dd 2;

S315,PO＝PP；

S316,POS＝POS+200；

s317, if the POS is less than the Length-100, turning to S311, otherwise, turning to S318;

s318, outputting a vibration mapping center vector TCenter, a maximum mapping distance YDist and a maximum offset distance PDist;

s4, continuously acquiring 200 vibration data for the power transformer by using vibration sensors ZA, ZB and ZC, respectively storing the acquired data into current A-phase data CurrentZA, current B-phase data CurrentZB and current C-phase data CurrentZC, and outputting winding deformation fault detection results;

s401, continuously acquiring 200 vibration data of the power transformer by using vibration sensors ZA, ZB and ZC, and respectively storing the acquired data into CurrentZA, CurrentZB and CurrentZC;

s402, TA is an element cut from 1 to 100 in CurrentZA, TB is an element cut from 1 to 100 in CurrentZB, and TC is an element cut from 1 to 100 in CurrentZC;

s403, calculating a first variable CurrentP1 of the current operator output result by using a mapping operator kopersonator (TA, TB, TC);

s404, currently temporarily storing a first variable TempD1 ═ CurrentP1-TCenter |;

where | | | represents l2norm of the calculation vector;

s405, TA being an element from 101 to 200 in CurrentZA, TB being an element from 101 to 200 in CurrentZB, and TC being an element from 101 to 200 in CurrentZC;

s406, calculating a second variable CurrentP2 of the current operator output result by using a mapping operator kopersonator (TA, TB, TC);

s407, currently temporarily storing the second variable TempD2 ═ CurrentP1-CurrentP2 |;

where | | | represents l2norm of the calculation vector;

s408, distance variation range index decision is 0.5 × (YDist-TempD1)/YDist +0.5 × (PDist-TempD 2)/PDist;

s409, if decision >0, the winding deformation fault is not generated, and the operation is switched to S411, otherwise, the operation is switched to S410;

s410, outputting: if the winding deformation fault occurs, the step goes to S412;

s411, outputting: if no winding deformation fault occurs, go to S412;

and S412, finishing the winding deformation fault judgment process.

The invention has the beneficial effects that:

providing a new detection method for the deformation fault of the power transformer winding; the method establishes a mapping operator, can map sensor vibration data in a certain range into a vector, establishes a tolerance range of vector change based on data in normal operation, and further realizes the detection of the deformation fault of the power transformer winding. Because the method described by the invention is mapped by the nonlinear function and is based on a certain change range, the method can not be bound with a specific numerical value, and can achieve a more stable prediction effect.

Detailed Description

The following examples are provided to further illustrate specific embodiments of the present invention, but it will be understood by those skilled in the art that these are merely examples and the scope of the present invention is defined by the appended claims, and those skilled in the art can make various changes or modifications to these embodiments without departing from the principle and spirit of the present invention, and these changes and modifications fall within the scope of the present invention.