CN112363012A - Power grid fault early warning device and method - Google Patents

Power grid fault early warning device and method Download PDF

Info

Publication number
CN112363012A
CN112363012A CN202011178835.8A CN202011178835A CN112363012A CN 112363012 A CN112363012 A CN 112363012A CN 202011178835 A CN202011178835 A CN 202011178835A CN 112363012 A CN112363012 A CN 112363012A
Authority
CN
China
Prior art keywords
bus
signal
power grid
early warning
characteristic
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
Application number
CN202011178835.8A
Other languages
Chinese (zh)
Inventor
杜洁
赵轶
王浩
赵旭伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Power Grid Hebei Electric Power Co ltd Xinle City Power Supply Branch
State Grid Corp of China SGCC
State Grid Hebei Electric Power Co Ltd
Original Assignee
State Power Grid Hebei Electric Power Co ltd Xinle City Power Supply Branch
State Grid Corp of China SGCC
State Grid Hebei Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Power Grid Hebei Electric Power Co ltd Xinle City Power Supply Branch, State Grid Corp of China SGCC, State Grid Hebei Electric Power Co Ltd filed Critical State Power Grid Hebei Electric Power Co ltd Xinle City Power Supply Branch
Priority to CN202011178835.8A priority Critical patent/CN112363012A/en
Publication of CN112363012A publication Critical patent/CN112363012A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/086Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/088Aspects of digital computing

Abstract

The invention relates to the technical field of bus fault diagnosis, along with the penetration of artificial intelligence in various industrial aspects, a plurality of intelligent algorithms are developed for fault diagnosis of a bus of a power grid, mainly a fault identification model of a power grid bus group is established by a neural network and the like, but the algorithm needs to be provided with a huge bus temperature and temperature signal library, so that the difficulty in the process of establishing a complete database is large, and the method is in an experimental stage. The invention discloses a power grid fault early warning device and a method, which comprises the following steps: establishing a data storage unit; extracting a signal X to be detected and a reference signal Y, and establishing a signal characteristic matrix A to be detected and a reference signal characteristic matrix B for the signal X to be detected and the reference signal Y respectively; respectively calculating the numerical values of the feature vectors at the same position in the A and the B based on an entropy weight method; and constructing an early warning index a and a threshold value matrix alfa, judging whether the bus of the power grid is normal or not according to the early warning threshold value and the early warning index, and establishing a state table Tab.

Description

Power grid fault early warning device and method
Technical Field
The invention relates to the technical field of power distribution, in particular to a power grid fault early warning device and method. The core idea is to calculate an early warning index through a signal to be detected and a reference signal, construct a threshold value matrix of the early warning index, and finally comprehensively judge whether the bus of the power grid is in fault or not through comparing and analyzing the early warning index and the threshold value thereof.
Background
In recent years, the research and construction of smart grids has risen globally. The intelligent power distribution network is used as an important component for connecting a main network and supplying power to users in the intelligent power grid, and whether the running state of the intelligent power distribution network is normal or not directly influences the power supply of thousands of households. Meanwhile, with the access of distributed power sources, the popularization of electric automobiles and the increase of user interaction electric power, the dynamic behavior of the power distribution network becomes complex, the operation risk is greatly increased, and once a power failure accident of the power distribution network occurs, huge influence and loss can be caused to social life. Therefore, a need exists for further research on fault early warning of the intelligent power distribution network, and guidance and help are provided for relevant management personnel to maintain the intelligent power distribution network.
At present, scholars at home and abroad propose various solutions from different angles aiming at the fault early warning of the intelligent power distribution network. Research shows that most faults of the power distribution network enter pathological operation before destructive faults occur, and the power distribution network has trend and cumulative effects. However, most of the existing solutions utilize local parameters of the power distribution network, such as harmonic current and short-circuit current, or local components, such as a transformer, or are associated with external factors, such as thunderstorm weather, to achieve the effect of early warning the fault of the intelligent power distribution network. The solutions fail to grasp the operation state of the intelligent power distribution network on the whole, and only consider a certain current factor to judge whether the early warning is needed or not in the early warning scheme, and the trend and the cumulative effect of the power distribution network faults cannot be fully utilized, so that the accuracy of fault early warning is slightly deficient in the face of more and more complex intelligent power distribution networks. Therefore, there is an urgent need for an early warning method that can realize the overall control of the operating state of the smart distribution network and can jointly determine whether a failure is about to occur from the current and past operating states.
Disclosure of Invention
The invention aims to provide a power grid fault early warning device and a power grid fault early warning method, which take a series of practical problems into consideration, adopt signals when a bus of a power grid is normal, and establish threshold value matrixes of maximum early warning indexes of different voltage levels and different power grids when the bus is normal. When the bus is diagnosed and pre-warned, the pre-warning index value of the bus is calculated, and whether the bus is pre-warned or not is indicated by comparing and analyzing the threshold value of the bus under the same voltage level. Further, an energy characteristic matrix is extracted from a normal temperature signal of the bus of the power grid by utilizing wavelet transformation, a threshold value matrix of a fault early warning index is established based on an entropy weight method, the early warning index of the bus to be detected obtained through calculation is compared and analyzed with a corresponding threshold value, the steps are repeated by changing the number of layers of wavelet decomposition, and finally, whether the bus to be detected carries out fault early warning or not is comprehensively analyzed.
The technical scheme provided by the invention is that the method comprises the following steps,
a power grid fault early warning device and method comprises the following specific processes:
step 1, establishing a data storage unit, wherein the storage unit periodically stores and updates all bus temperature signals of different voltage levels in a power grid bus group, which are extracted under the fault-free condition; the temperature signal comprises the actual working temperature of all buses under different voltage levels in the power grid bus group; extracting the temperature signal as a standby reference signal for early warning of faults of all buses under different voltage levels in the power grid bus group;
step 2, extracting a signal to be detected X and a reference signal Y from the bus of the power grid, respectively performing wavelet decomposition on the signal to be detected X and the reference signal Y, extracting respective energy coefficient characteristic attribute values, wherein each energy coefficient characteristic attribute value forms an energy coefficient characteristic attribute characteristic dimension, and further respectively constructing a signal to be detected characteristic matrix A and a reference signal characteristic matrix B;
the step 3, respectively extracting the energy coefficient characteristic attribute values at the same position of the nth column of the to-be-detected signal characteristic matrix a and the reference signal characteristic matrix B to form an m × 2 characteristic matrix Q, where n is {1,2,3, … … K }, K is the number of layers of the wavelet decomposition, and m is the number of segments of the signal division; respectively calculating the characteristic weight values of the first dimension and the second dimension of the energy coefficient characteristic dimension in the K characteristic matrixes Q by using an entropy weight method, wherein the characteristic weight values are respectively bn1And bn2
And 4, constructing an index a of the power grid bus group fault early warning to quantitatively describe the operating temperature of the bus, wherein the formula is as follows:
Figure BDA0002749488720000021
constructing a threshold value matrix alfa by the bus fault early warning index a of the power grid;
and 5, traversing the threshold value matrix alfa according to the serial number and the temperature of the bus to be detected, and determining the early warning threshold value a of the bus to be detectedijAnd the bus fault early warning index a of the power grid, and establishing a state information table Tab of the power grid bus group;
and 6, adjusting the number K of wavelet decomposition layers and the change times thereof, recording the calculation result of each time in a bus state information table Tab in the power grid bus group, and judging the operation conditions of all buses in the power grid bus group under different voltage levels.
Further, in step 2, a certain number of sampling points are set for the signal to be measured X and the reference signal Y as sample lengths, which are defined as wlen, and then the signal to be measured X and the reference signal Y are respectively divided into m segments of signals, that is, the signal to be measured X and the reference signal Y are divided into m segments of signals
Figure BDA0002749488720000031
Wherein m is an integer part of a result in the formula, the number of sampling points of wlen at least comprises temperature sampling points of each voltage section of the bus of the power grid, and X, Y represents the signal to be measured and the reference signal.
Further, in the wavelet decomposition in the step 2, a wavelet energy formula is used to calculate energy coefficients after the decomposition of the reference signal Y and the signal to be detected X, that is, En=∑|xn|2Wherein xnEach decomposed signal segment corresponds to a reference value, n is {1,2,3.. K }, and K is the number of layers of wavelet decomposition of the signal; the signal X to be detected and the reference signal Y are respectively divided into m sections of signals, each section of signal is respectively subjected to K-layer wavelet decomposition, the signal energy coefficient after each section of wavelet decomposition is an energy coefficient characteristic attribute value, the signal energy coefficient characteristic attribute after each section of wavelet decomposition is composed of energy coefficient characteristic dimensions, and the m sections of energy coefficient characteristic dimensions are composed of a signal characteristic matrix A to be detected and a reference signal characteristic matrix B.
Further, the step 3 is to calculate entropy weight methods adopted for feature weight values of the feature attributes of the first dimension and the second dimension in the K feature matrices Q, respectively, that is, firstly, normalization is performed according to each feature attribute value in the feature matrices Q, that is, normalization is performed
Figure BDA0002749488720000032
Where i ═ 1,2,3.. ·. m }, j ═ 1,2}, min (X)i) Is the minimum value, max (X), of the energy coefficient characteristic attribute value corresponding to the temperature signal in the ith sectioni) Is the maximum value, X, of the energy coefficient characteristic attribute value corresponding to the temperature signal in the ith sectionijAnd obtaining a normalized matrix Q' of the characteristic matrix Q for the energy coefficient characteristic attribute value corresponding to the temperature signal of the ith dimension of the ith section.
Further, in the step 4, the bus fault early warning index a is constructed
Figure BDA0002749488720000033
Wherein b isn1And bn2And respectively representing the characteristic weight values of the energy coefficient characteristic attributes in the first dimension and the second dimension in the characteristic matrix Q.
Further, the threshold value matrix alfa in step 4 extracts temperature signals of the same voltage class and the same bus from the temperature signal storage units in three time periods, and the bus fault early warning indicators a of the power grid are calculated according to steps 2 to 4, and are respectively calculated by using a1、a2、a3If the voltage level of the bus is represented, the fault early warning index a of the bus at the voltage level is the maximum value of the three values; and (3) performing the same treatment on other buses in the power grid bus group to obtain a threshold value matrix alfa of the early warning index a:
Figure BDA0002749488720000041
wherein u is the total number of the bus numbers in the power grid bus group, v is the number of the buses in the power grid bus group at different temperatures, and the value range of i is [1, u]Is an integer between, j has a value in the range of [1, v]Is an integer of (1).
Further, in the step 5, the early warning threshold value a of the bus to be testedijComparing the fault early warning index a of the bus of the power grid with the fault early warning index a of the bus of the power grid, and when a is smaller than aijIf so, the temperature of the bus to be tested is normal, otherwise, the bus early warning index of the power grid to be tested exceeds the bus early warning threshold value of the power grid, a corresponding bus state information table Tab of the power grid is established, and a diagnosis result is recorded.
Further, according to the energy coefficient in the normalization matrix Q', a calculation formula is used
Figure BDA0002749488720000042
Figure BDA0002749488720000043
Wherein
Figure BDA0002749488720000044
j={1,2},min(Xi) The most energy coefficient characteristic attribute value corresponding to the temperature signal in the ith sectionSmall value, max (X)i) Is the maximum value, X, of the energy coefficient characteristic attribute value corresponding to the temperature signal in the ith sectionijThe characteristic attribute value of the energy coefficient corresponding to the temperature signal of the ith dimension is the ith dimension, if pijIs defined when 0
Figure BDA0002749488720000045
And obtaining the characteristic weight value in the normalization matrix Q'.
Further, the signal X to be detected in step 2 is a section of temperature signal generated by a certain bus in the power grid bus group at a fixed time by the bus of the power grid; the reference signal Y is a temperature signal of the bus corresponding to the signal to be detected X under the same voltage level, and the retention time and the sample length of the signal to be detected X and the sample length of the reference signal Y are consistent.
Further, in the step 2, K layers of wavelet decomposition are respectively performed on the same temperature signal in the two signals, i.e., the signal to be detected X and the reference signal Y, where K is a constant set artificially, and a wavelet energy formula is used to extract respective energy coefficients after decomposition, so as to respectively generate a corresponding signal to be detected feature dimension X1 and a corresponding reference signal feature dimension Y1; and performing the same processing on the temperature signals in the residual sections, and respectively generating a characteristic matrix A of the signal to be detected and a characteristic matrix B of the reference signal, wherein the sizes of the characteristic matrix A and the characteristic matrix B are both mxK, m is the number of rows of the characteristic matrix, namely the number of sections of a bus, and K is the number of columns of the characteristic matrix, namely the number of layers of wavelet decomposition.
Furthermore, the power grid fault early warning device and method further comprise a fault early warning unit, wherein the fault early warning unit is electrically connected with a central processing unit, the central processing unit is electrically connected with a server, the server is provided with a service storage card, and the storage card is electrically connected with the data storage unit.
One aspect of the technical effect brought by the technical scheme of the invention is that the difficulty of establishing a complete data set in the process of establishing the whole power grid fault identification model by adopting a neural network and the like is avoided, the accuracy and the stability of the bus fault early warning of the power grid are improved, and the method is closer to practical application.
Drawings
FIG. 1 is a flow chart of a power grid fault early warning apparatus and method of the present invention;
Detailed Description
Example 1
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The specific implementation method is as follows:
as shown in fig. 1, a flow chart of a power grid fault early warning apparatus and method of the present invention is provided, which includes the steps of:
step 1, aiming at a power grid, establishing a data storage unit of a bus temperature signal of the power grid, wherein a collected bus comprises temperature signals of all power grids under different voltage grade scales, the temperature signals are used as standby reference signals when the bus fault of the power grid is early warned, the storage unit needs to be completely updated at regular intervals, and a time period is measured by month;
step 2, at a fixed moment, when fault diagnosis is performed on the bus of a power grid, acquiring a section of temperature signal sent by the bus temperature at the moment as a signal to be measured X, setting a certain number of sampling points as sample lengths, and defining the sampling points as wlen, wherein the signal to be measured can be divided into m sections of temperature signals, and a specific calculation formula is as follows:
Figure BDA0002749488720000051
wherein m is an integer part of a result in the formula, the number of sampling points of wlen at least comprises temperature sampling points of each voltage section of the bus of the power grid, and X, Y represents the signal to be measured and the reference signal.
Searching a bus which is the same as the signal to be detected and a temperature signal under the same voltage level in a reference signal storage unit, wherein the time and the sample length of the reference signal are consistent with those of the signal to be detected;
performing K-layer wavelet decomposition on each temperature signal of the two signals, wherein K is a constant set artificially, and extracting an energy coefficient after decomposition by using a wavelet energy formula, wherein the calculation formula is as follows:
En=∑|xn|2 (3)
wherein xnFor each decomposed signal, n ═ 1,2,3.. K }, where K is the number of layers in the wavelet decomposition and K ═ 1,2,3, … … K }.
Thus, each section of temperature signal can be represented by a vector generated by an energy coefficient, and energy characteristic matrixes A and B capable of representing the essence of two signals can be extracted from the signal to be measured and the reference signal through the transformation, wherein the large size of the two matrixes is m multiplied by K, namely the row number of the matrixes is the number of sections of a bus, and the column of the matrixes is the number of layers of wavelet decomposition;
step 3, respectively extracting the energy coefficient characteristic attribute values of the same position of the nth column of the to-be-detected signal characteristic matrix A and the reference signal characteristic matrix B to form an m × 2 characteristic matrix Q, wherein the total number of the characteristic matrices is K, n is {1,2,3, … … K }, K is the number of layers of the wavelet decomposition, and m is the number of segments of the signal division; respectively calculating the characteristic weight values of the first dimension and the second dimension of the energy coefficient characteristic dimension in the K characteristic matrixes Q by using an entropy weight method, wherein the characteristic weight values are respectively bn1And bn2
Said step 3 calculates said entropy weight method used for the first and second dimension characteristic weight values in K matrices Q, i.e. normalizes according to each said energy coefficient value in said characteristic matrix Q, i.e. normalizes
Figure BDA0002749488720000061
Where i ═ 1,2,3.. ·. m }, j ═ 1,2}, min (X)i) Energy corresponding to the temperature signal in the ith sectionMinimum value of the quantity coefficient, max (X)i) Is the maximum value, X, of the energy coefficient corresponding to the temperature signal in the ith sectionijAnd normalizing the energy coefficient corresponding to the temperature signal of the ith dimension to obtain a data normalization table so as to obtain a normalization matrix Q' of the characteristic matrix Q. According to the energy coefficient in the normalization matrix Q', a calculation formula is used
Figure BDA0002749488720000062
Wherein
Figure BDA0002749488720000063
If p isijIs defined when 0
Figure BDA0002749488720000071
And obtaining the characteristic weight value in the normalization matrix Q'.
And 4, constructing an index a of the bus fault early warning of the power grid bus group for quantitatively describing the operating temperature of the bus, wherein the formula is as follows:
Figure BDA0002749488720000072
and constructing a threshold value matrix alfa by the bus fault early warning index a of the power grid.
Step 5, constructing a threshold value matrix alfa of the early warning index a by using a storage unit of the temperature signal,
and establishing storage units of temperature signals in three time periods, respectively taking out the temperature signals of the same bus under the same voltage level from the three storage units, calculating a indexes among the temperature signals according to the steps 2 to 4, respectively representing the values of the three indexes which can be calculated under the fixed working condition of the bus by a1, a2 and a3, and then taking the maximum value of the three indexes as the fault early warning index of the bus under the voltage level. The threshold value matrix alfa of the early warning index a can be obtained by performing the same processing on different buses:
Figure BDA0002749488720000073
in the formula (5), u is the total number of the serial numbers of the buses of a certain power grid, v is the number of different voltage grade scales of the power grid, the value range of i is an integer between [1 and u ], the value range of j is an integer between [1 and v ],
step 6, determining the early warning threshold value of the bus from alfa according to the serial number and the temperature condition of the bus to be detected, and when the calculated a is less than aijIf so, the bus temperature is normal, otherwise, early warning is carried out, a bus equipment manager is reminded that the bus possibly has a fault in the power grid, and the obtained result is recorded in a bus temperature state information table Tab;
and 7, changing the number K of wavelet decomposition layers, wherein the number of wavelet decomposition layers involved in the steps 2 to 6 must be kept consistent, the number of times of circulation of the whole process can be realized by setting the change times of the K values, each result is recorded in a bus temperature state information table Tab, if the early warning times in the table are larger than the normal times, the bus early warns the bus fault, otherwise, the bus temperature is displayed to be normal, and the fault early warning of the bus of the power grid is realized.
Figure BDA0002749488720000074
Figure BDA0002749488720000081
Example 2
The embodiment further describes an implementation process and cautions of the power grid fault early warning device and method by combining with an example of a field.
At present, all buses in 5 power grid bus groups need to be established with an entropy weight method to form a power grid bus fault early warning system for managing on-site bus equipment. The specific method comprises the following steps:
1. extracting a temperature signal Y to be measured and a reference temperature signal X under a certain voltage level, and dividing the temperature signal Y to be measured and the reference temperature signal X into 4 sections for 3-layer wavelet decomposition, so that the measured reference temperature signal Y and the measured temperature signal X are as follows:
the reference temperature signal to be measured
Figure BDA0002749488720000082
The reference temperature signal
Figure BDA0002749488720000091
2. Respectively extracting the energy coefficients of the nth columns of the signal characteristic matrix A to be detected and the reference signal characteristic matrix B, and combining the energy coefficients into an m × 2 matrix Q, wherein n is {1,2,3, … … K }, and K matrices Q are counted, where m is 4 and K is 3, so that the characteristic matrices Q are respectively:
Figure BDA0002749488720000092
feature weights for the first and second dimensions, bn1 and bn2 respectively, were calculated using entropy weight methods for Q1, Q2 and Q3 respectively, where n ═ 1,2,3,4 }. Firstly according to a normalization formula
Figure BDA0002749488720000093
Obtaining the respective normalized Q11, Q21, Q31 and Q41, namely:
Figure BDA0002749488720000094
according to a calculation formula
Figure BDA0002749488720000095
Wherein
Figure BDA0002749488720000096
If p isijIs defined when 0
Figure BDA0002749488720000097
And obtaining characteristic weight values, namely bn1 and bn2, in the normalization matrix Q'.
Figure BDA0002749488720000098
Constructing an index a of the bus fault early warning for quantitatively describing the operating temperature of the bus, wherein a calculation formula is as follows:
Figure BDA0002749488720000099
and a | b11-b12| + | b21-b22| + | b31-b32| 0.012+0.278+0.16 | -0.45, and the bus fault early warning index of the power grid is 0.45.
3. And constructing a threshold value matrix alfa of the early warning index a. And establishing storage units of temperature signals in three time periods, adopting the temperature signals of the same bus under the same voltage level, calculating an early warning index a according to the steps, and calculating values of the three indexes under the fixed working condition of the bus, wherein the values are respectively represented by a1, a2 and a3, so that the fault early warning index of the bus under the voltage level is the maximum value of the three indexes. The specific process is as follows, assuming that the bus is a bus No. 1, in three time slots, namely a time slot 1, a time slot 2, a time slot 3 and the like, a temperature signal is extracted respectively under the condition that the bus has no fault and under the actual operation state, and the temperature signal is divided into 4 segments to carry out 3-layer wavelet decomposition, and a series of operations are carried out to obtain corresponding early warning indexes a1, a2 and a 3. The measured reference temperature signal Y and the measured temperature signal X are as follows:
Figure BDA0002749488720000101
the Q matrix corresponding to each group of temperature signals extracts the same row energy coefficient of the reference temperature signal X and the temperature signal Y to be measured respectively, and combines the Q matrix as follows:
Figure BDA0002749488720000102
Figure BDA0002749488720000111
carrying out entropy weight normalization processing on each Q matrix to obtain each normalization matrix, wherein the normalization matrixes are as follows:
Figure BDA0002749488720000112
according to a calculation formula
Figure BDA0002749488720000113
Wherein
Figure BDA0002749488720000114
If p isijIs defined when 0
Figure BDA0002749488720000115
The characteristic weight values in the normalization matrix Q' are derived, i.e. bi1, bi 2. The details are set forth in the following table:
Figure BDA0002749488720000116
the fault early warning index a is used for quantitatively describing the operating temperature of the bus, and the calculation formula is as follows:
Figure BDA0002749488720000121
wherein b isn1And bn2The characteristic weight values respectively represent the characteristic attributes of the energy coefficients in the first dimension and the second dimension of the characteristic matrix Q, so a1, a2, a3 are respectively: 0.12, 0.66, 0.5. The fault early warning index of the No. 1 bus at the voltage level is the maximum value of the three indexes, namely 0.66.
The operation is also carried out on other buses, and the bus fault early warning system is established for the field bus equipment by utilizing an entropy weight method for 5 buses in the field. Extracting a reference temperature signal Y and a temperature signal X to be detected under a certain voltage level, dividing the reference temperature signal Y and the temperature signal X to be detected into 4 sections for 3-layer wavelet decomposition, and sampling a bus in three speed sections, wherein v is 3, u is 5, and the specific steps are as follows:
Figure BDA0002749488720000122
the specific data are presented in the following table, as follows:
power grid bus group Voltage class 1 Voltage class 2 Voltage class 3
No. 1 bus 0.66 0.45 0.42
No. 2 bus 0.64 0.42 0.27
No. 3 bus 0.78 0.38 0.56
No. 4 bus 0.45 0.31 0.69
Number 5 busThread 0.38 0.39 0.36
Calculating a fault early warning index 0.45 of the bus of the power grid bus group power grid according to the step 2 in the embodiment, and comparing the fault early warning actual values of the buses of the power grid bus group in the upper table under different voltage levels, namely determining an early warning threshold value of the bus from alfa according to the serial number and the temperature condition of the bus to be detected, and when a is smaller than aijAnd if not, giving an early warning to remind a bus equipment manager that the bus possibly has a fault and recording the obtained result in a bus temperature state information table Tab, wherein the bus temperature is normal, and the result is specifically as follows:
Figure BDA0002749488720000123
Figure BDA0002749488720000131
finally, the adjustment of the K value and the number of times, namely the number of layers of wavelet decomposition, is tried, and each result is recorded in a bus temperature state information table Tab, as shown in the following table, so as to be used for early warning of the bus fault of the power grid bus group.
Figure BDA0002749488720000132
Figure BDA0002749488720000141
Example 3
The difference between this embodiment and embodiment 2 is that the characteristic attribute weight in the normalization matrix Q' is subjected to secondary entropy weighting by using an entropy weighting method, and the difference and the relation between the final result and embodiment 1 are compared. The following table shows the characteristic weights in the normalized matrix Q' obtained in example 2:
Figure BDA0002749488720000142
according to the weight calculation formula of each characteristic dimension of the characteristic attribute of the bus sample
Figure BDA0002749488720000143
The weights of the feature dimensions in the feature matrix Y can be obtained as shown in the following table:
Figure BDA0002749488720000144
formula is used for index a of fault early warning at the moment
Figure BDA0002749488720000151
The calculated a1, a2 and a3 are respectively: 0.05, 0.08 and 0.03, and the fault early warning index of the No. 1 bus at the voltage level is the maximum value of the three indexes, namely 0.08. The operation processing is also carried out on the bus, and the bus fault early warning system is established for the field bus equipment by utilizing an entropy weight method for 5 buses in the field. Extracting a reference temperature signal Y and a temperature signal X to be detected under a certain voltage level, dividing the reference temperature signal Y and the temperature signal X to be detected into 4 sections for 3-layer wavelet decomposition, and sampling a bus in three speed sections, wherein v is 3, u is 5, and the specific steps are as follows:
Figure BDA0002749488720000152
the specific data are presented in the following table, as follows:
power grid bus group Voltage class 1 Voltage class 2 Voltage class 3
No. 1 bus 0.35 0.02 0.06
No. 2 bus 0.78 0.05 0.03
No. 3 bus 0.09 0.02 0.08
No. 4 bus 0.02 0.03 0.35
No. 5 bus 0.01 0.06 0.08
According to the embodiment, the fault early warning index of the bus of the power grid bus group power grid is 0.08, and the reference tableComparing fault early warning actual values of all buses of the medium power grid bus group under different voltage levels, namely determining an early warning threshold value of each bus from alfa according to the serial number and the temperature condition of the bus to be detected, and when calculated a is less than aijAnd if not, giving an early warning to remind a bus equipment manager that the bus possibly has a fault and recording the obtained result in a bus temperature state information table Tab, wherein the bus temperature is normal, and the result is specifically as follows:
Figure BDA0002749488720000153
Figure BDA0002749488720000161
finally, the adjustment of the K value and the number of times, namely the number of layers of wavelet decomposition, is tried, and each result is recorded in a bus temperature state information table Tab, as shown in the following table, so as to be used for early warning of the bus fault of the power grid bus group.
Figure BDA0002749488720000162
Figure BDA0002749488720000171
The power grid fault early warning device and method further comprise a fault early warning unit, wherein the fault early warning unit is electrically connected with a central processing unit, the central processing unit is electrically connected with a server, the server is provided with a service storage card, and the storage card is electrically connected with the data storage unit.
According to the calculation result of the embodiment, the secondary weight of the energy coefficient matrix characteristic dimension of all power grids of the power grid bus group extracted based on the entropy weight method is consistent with the result of the embodiment 2, namely, the power grid bus group fault can be early-warned and judged by solving the characteristic attribute weight of the normalized matrix or the secondary weight of the characteristic dimension of the normalized matrix by adopting the entropy weight method.
Example 4
In this embodiment, whether the sequence of the reference temperature signal Y or the reference temperature signal X corresponding to the first characteristic dimension and the second characteristic dimension affects the final judgment result will be further described. Similarly, assuming that all power grids of a certain power grid bus group in a certain bus field total 5, performing 3-layer wavelet decomposition processing on the feature matrix extracted by the bus of all power grids of 5 buses, and obtaining the following result:
the temperature signal to be measured
Figure BDA0002749488720000172
The reference temperature signal
Figure BDA0002749488720000173
Respectively extracting the energy coefficients of the nth columns of the signal characteristic matrix to be detected and the reference signal characteristic matrix, and combining the energy coefficients into an m × 2 matrix Q, wherein n is {1,2,3, … … K }, K matrices Q are counted, where m is 4 and K is 3, so that the characteristic matrices Q are respectively:
Figure BDA0002749488720000181
at this time, the characteristic dimension of the characteristic matrix of the signal to be detected is a first characteristic dimension of the characteristic matrix Q, and the characteristic dimension of the characteristic matrix of the reference signal is a second characteristic dimension of the characteristic matrix Q. Actually, an index a of fault early warning is obtained to quantitatively describe the operating temperature of the bus, and a calculation formula is as follows:
Figure BDA0002749488720000182
wherein b isn1And bn2Respectively representing the characteristic weight values of the energy coefficient characteristic attributes of the first dimension and the second dimension in the characteristic matrix Q, and knowing the first characteristic dimension and the second characteristic dimension of the characteristic matrix QThe second characteristic dimension respectively corresponds to the sequence of the reference temperature signal Y or the temperature signal X to be measured, and has no influence on the bus fault early warning index a of the power grid.

Claims (10)

1. The power grid fault early warning device and the method are characterized by comprising the following steps of 1, establishing a data storage unit, wherein the data storage unit is used for periodically storing and updating all bus temperature signals extracted by a power grid bus group under the fault-free condition and used as standby reference signals when all bus faults are early warned; step 2, constructing a characteristic matrix A of the signal to be detected and a characteristic matrix B of the reference signal; step 3, respectively extracting the energy coefficient characteristic attribute values of the same position of the nth column of the to-be-detected signal characteristic matrix A and the reference signal characteristic matrix B to form an m × 2 characteristic matrix Q, wherein the total number of the characteristic matrices is K, n is {1,2,3, … … K }, K is the number of layers of the wavelet decomposition, and m is the number of segments of the signal division; step 4, constructing an index a of the power grid bus group fault early warning, and constructing a threshold value matrix alfa; step 5, traversing the threshold value matrix alfa, and determining the early warning threshold value a of the bus to be testedijAnd the bus fault early warning index a of the power grid; and 6, adjusting the size of the wavelet decomposition layer number K and the change times thereof, and judging the running state of the bus.
2. The power grid fault early warning device and method according to claim 1, wherein: in the step 2, a certain number of sampling points are set for the signal to be measured X and the reference signal Y as sample lengths, defined as wlen, and then the signal to be measured X and the reference signal Y are respectively divided into m segments of signals, that is, the signals are
Figure FDA0002749488710000011
Wherein m is an integer part of a result in the formula, the number of sampling points of wlen at least comprises temperature sampling points of each voltage section of the bus of the power grid, and X, Y represents the signal to be measured and the reference signal.
3. A power grid fault pre-warning system as claimed in claim 1The device and the method are characterized in that: in the step 2, the wavelet decomposition is carried out, and a wavelet energy formula is utilized to calculate energy coefficients, namely E, of the decomposed reference signal Y and the decomposed energy coefficients of the signal X to be detectedn=∑|xn|2Wherein xnEach decomposed signal segment corresponds to a reference value, n is {1,2,3.. K }, and K is the number of layers of wavelet decomposition of the signal; the signal X to be detected and the reference signal Y are respectively divided into m sections of signals, each section of signal is respectively subjected to K-layer wavelet decomposition, the signal energy coefficient after each section of wavelet decomposition is an energy coefficient characteristic attribute value, the signal energy coefficient characteristic attribute after each section of wavelet decomposition is composed of energy coefficient characteristic dimensions, and the m sections of energy coefficient characteristic dimensions are composed of a signal characteristic matrix A to be detected and a reference signal characteristic matrix B.
4. The power grid fault early warning device and method according to claim 1, wherein: step 3 is to calculate the entropy weight method adopted by the feature weight values of the first dimension and the second dimension of the feature attributes in the K feature matrices Q respectively, namely, firstly, normalization is carried out according to each feature attribute value in the feature matrices Q, namely
Figure FDA0002749488710000021
Where i ═ 1,2,3.. ·. m }, j ═ 1,2}, min (X)i) Is the minimum value, max (X), of the energy coefficient characteristic attribute value corresponding to the temperature signal in the ith sectioni) Is the maximum value, X, of the energy coefficient characteristic attribute value corresponding to the temperature signal in the ith sectionijAnd obtaining a normalized matrix Q' of the characteristic matrix Q for the energy coefficient characteristic attribute value corresponding to the temperature signal of the ith dimension of the ith section.
5. The power grid fault early warning device and method according to claim 1, wherein: in the step 4, the bus fault early warning index a is constructed
Figure FDA0002749488710000022
Wherein b isn1And bn2And respectively representing the characteristic weight values of the energy coefficient characteristic attributes in the first dimension and the second dimension in the characteristic matrix Q.
6. The power grid fault early warning device and method according to claim 1, wherein: and step 4, extracting temperature signals of the same voltage grade of the bus from signal storage units of three time periods respectively by using the threshold value matrix alfa, calculating the bus fault early warning index a of the power grid by using the temperature signals of the single time period, and respectively using the a1、a2、a3If the bus is represented, the fault early warning index a of the bus at the temperature is the maximum value of the three values; and (3) performing the same treatment on other buses in the power grid bus group to obtain a threshold value matrix alfa of the early warning index a:
Figure FDA0002749488710000023
wherein u is the total number of the bus numbers in the power grid bus group, v is the number of the buses in the power grid bus group at different temperatures, and the value range of i is [1, u]Is an integer between, j has a value in the range of [1, v]Is an integer of (1).
7. The power grid fault early warning device and method according to claim 1, wherein: the early warning threshold value a of the bus to be tested in the step 5ijComparing the fault early warning index a of the bus of the power grid with the fault early warning index a of the bus of the power grid, and when a is smaller than aijIf so, the temperature of the bus to be tested is normal, otherwise, the bus early warning index of the power grid to be tested exceeds the bus early warning threshold value of the power grid, a corresponding bus state information table Tab of the power grid is established, and a diagnosis result is recorded.
8. The power grid fault early warning device and method according to claim 4, wherein: according to the energy coefficient in the normalization matrix Q', a calculation formula is used
Figure FDA0002749488710000024
Wherein
Figure FDA0002749488710000025
j={1,2},min(Xi) Is the minimum value, max (X), of the energy coefficient characteristic attribute value corresponding to the temperature signal in the ith sectioni) Is the maximum value, X, of the energy coefficient characteristic attribute value corresponding to the temperature signal in the ith sectionijThe characteristic attribute value of the energy coefficient corresponding to the temperature signal of the ith dimension is the ith dimension, if pijIs defined when 0
Figure FDA0002749488710000031
And obtaining the characteristic attribute characteristic weight value in the normalized matrix Q'.
9. The power grid fault early warning device and method according to claim 1, wherein: in the step 2, the signal X to be detected is a section of temperature signal sent by a certain bus of the power grid bus group at a fixed moment by the bus of the power grid; the reference signal Y is the standby reference signal in the traversal data storage unit, a temperature signal of the bus corresponding to the signal to be detected X at the same voltage level is extracted, and the time and the sample length of the two signals of the signal to be detected X and the reference signal Y are consistent; step 2, performing K-layer wavelet decomposition on the same temperature signal in the two signals of the signal to be detected X and the reference signal Y respectively, wherein K is a constant set manually, and extracting respective energy coefficients after decomposition by using a wavelet energy formula to generate corresponding characteristic dimensions of the signal to be detected and the reference signal; and the temperature signals in the residual section are processed, the characteristic matrix A of the signals to be detected and the characteristic matrix B of the reference signals are respectively generated, the sizes of the characteristic matrix A and the characteristic matrix B are both mxK, wherein m is the row number of the characteristic matrix, namely the number of sections of a bus, and K is the number of columns of the characteristic matrix, namely the number of layers of wavelet decomposition.
10. The power grid fault early warning device and method as claimed in claim 1, further comprising a fault early warning unit, wherein the fault early warning unit is electrically connected with a central processing unit, the central processing unit is electrically connected with a server, the server is provided with a service memory card, and the memory card is electrically connected with the data storage unit.
CN202011178835.8A 2020-10-29 2020-10-29 Power grid fault early warning device and method Pending CN112363012A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011178835.8A CN112363012A (en) 2020-10-29 2020-10-29 Power grid fault early warning device and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011178835.8A CN112363012A (en) 2020-10-29 2020-10-29 Power grid fault early warning device and method

Publications (1)

Publication Number Publication Date
CN112363012A true CN112363012A (en) 2021-02-12

Family

ID=74512365

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011178835.8A Pending CN112363012A (en) 2020-10-29 2020-10-29 Power grid fault early warning device and method

Country Status (1)

Country Link
CN (1) CN112363012A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112415338A (en) * 2020-12-17 2021-02-26 石家庄嘉诚联信科技开发有限公司 Power grid fault detection method and device and storage medium
CN112994248A (en) * 2021-04-07 2021-06-18 李春娥 Power distribution network bus fault early warning device and method

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103941162A (en) * 2014-05-12 2014-07-23 福州大学 Resonant earthed system fault line selection method utilizing waveform time domain feature clustering
CN104502795A (en) * 2014-11-26 2015-04-08 国家电网公司 Intelligent fault diagnosis method suitable for microgrid
CN104698343A (en) * 2015-03-26 2015-06-10 广东电网有限责任公司电力调度控制中心 Method and system for judging power grid faults based on historical recording data
CN106249101A (en) * 2016-06-30 2016-12-21 湖南大学 A kind of intelligent distribution network fault identification method
CN106405339A (en) * 2016-11-11 2017-02-15 中国南方电网有限责任公司 Power transmission line fault reason identification method based on high and low frequency wavelet feature association
CN108562828A (en) * 2018-03-14 2018-09-21 哈尔滨理工大学 The method for improving electrical network low voltage ride-through capability based on Wavelet Detection
CN108663600A (en) * 2018-05-09 2018-10-16 广东工业大学 A kind of method for diagnosing faults, device and storage medium based on power transmission network
CN110488152A (en) * 2019-09-27 2019-11-22 国网河南省电力公司电力科学研究院 A kind of distribution network fault line selection method based on Adaptive Neuro-fuzzy Inference
CN110598281A (en) * 2019-08-28 2019-12-20 桂林理工大学 Entropy weight method based normal cloud model karst collapse prediction analysis method

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103941162A (en) * 2014-05-12 2014-07-23 福州大学 Resonant earthed system fault line selection method utilizing waveform time domain feature clustering
CN104502795A (en) * 2014-11-26 2015-04-08 国家电网公司 Intelligent fault diagnosis method suitable for microgrid
CN104698343A (en) * 2015-03-26 2015-06-10 广东电网有限责任公司电力调度控制中心 Method and system for judging power grid faults based on historical recording data
CN106249101A (en) * 2016-06-30 2016-12-21 湖南大学 A kind of intelligent distribution network fault identification method
CN106405339A (en) * 2016-11-11 2017-02-15 中国南方电网有限责任公司 Power transmission line fault reason identification method based on high and low frequency wavelet feature association
CN108562828A (en) * 2018-03-14 2018-09-21 哈尔滨理工大学 The method for improving electrical network low voltage ride-through capability based on Wavelet Detection
CN108663600A (en) * 2018-05-09 2018-10-16 广东工业大学 A kind of method for diagnosing faults, device and storage medium based on power transmission network
CN110598281A (en) * 2019-08-28 2019-12-20 桂林理工大学 Entropy weight method based normal cloud model karst collapse prediction analysis method
CN110488152A (en) * 2019-09-27 2019-11-22 国网河南省电力公司电力科学研究院 A kind of distribution network fault line selection method based on Adaptive Neuro-fuzzy Inference

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112415338A (en) * 2020-12-17 2021-02-26 石家庄嘉诚联信科技开发有限公司 Power grid fault detection method and device and storage medium
CN112994248A (en) * 2021-04-07 2021-06-18 李春娥 Power distribution network bus fault early warning device and method
CN112994248B (en) * 2021-04-07 2023-05-26 斯普屹科技(北京)有限公司 Power distribution network bus fault early warning device and method

Similar Documents

Publication Publication Date Title
CN112699913B (en) Method and device for diagnosing abnormal relationship of household transformer in transformer area
CN110097297B (en) Multi-dimensional electricity stealing situation intelligent sensing method, system, equipment and medium
CN109871976B (en) Clustering and neural network-based power quality prediction method for power distribution network with distributed power supply
Bashkari et al. Outage cause detection in power distribution systems based on data mining
CN112149873B (en) Low-voltage station line loss reasonable interval prediction method based on deep learning
CN104966161B (en) A kind of power quality recorder data calculation and analysis methods based on gauss hybrid models
CN111950585A (en) XGboost-based underground comprehensive pipe gallery safety condition assessment method
CN110570012B (en) Storm-based power plant production equipment fault early warning method and system
CN112734128A (en) 7-day power load peak value prediction method based on optimized RBF
CN116125361B (en) Voltage transformer error evaluation method, system, electronic equipment and storage medium
CN112363012A (en) Power grid fault early warning device and method
CN110826228B (en) Regional power grid operation quality limit evaluation method
CN114519514B (en) Low-voltage transformer area reasonable line loss value measuring and calculating method, system and computer equipment
CN111723839A (en) Method for predicting line loss rate of distribution room based on edge calculation
CN110210670A (en) A kind of prediction technique based on power-system short-term load
CN111709668A (en) Power grid equipment parameter risk identification method and device based on data mining technology
CN112949201B (en) Wind speed prediction method and device, electronic equipment and storage medium
CN117114161A (en) Method for predicting wind deflection flashover risk of power transmission line based on meta-learning
CN112415338A (en) Power grid fault detection method and device and storage medium
CN115146715A (en) Power utilization potential safety hazard diagnosis method, device, equipment and storage medium
CN114066068A (en) Short-term power load prediction method, device, equipment and storage medium
CN114091344A (en) Power transmission line risk assessment model training method and device based on data coupling
CN114446019A (en) Alarm information processing method, device, equipment, storage medium and product
CN115693692A (en) Voltage qualification rate improving method based on power distribution network voltage data analysis
CN214036158U (en) Water pump fault early warning device based on entropy weight method

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: 20210212