Disclosure of Invention
The embodiment of the invention provides a method and a system for eliminating SCADA bad data based on current unbalance, and aims to solve the problems of poor practicability and large error of the existing method.
In order to solve the technical problem, the embodiment of the invention discloses the following technical scheme:
the invention provides a method for eliminating SCADA bad data based on the current unbalance of a power transmission line, which comprises the following steps:
s1, obtaining SCADA measurement data at two ends of the line to be measured;
s2, converting the measured data into single-phase measured data;
s3, acquiring the ground-to-ground electric parameter of the line to be tested, and setting the length of a confidence interval;
and S4, calculating the unbalance amount of the line current to be measured according to the single-phase measurement data, determining the end point of the confidence interval, and eliminating the measurement data distributed outside the confidence interval.
Further, in step S1, measurement data is acquired every fixed time interval; the confidence interval in step S4 is used to detect the measured data acquired in all current time periods.
Further, the SCADA measurement data comprises a voltage amplitude, a current amplitude and active power at two ends of the line to be measured.
Further, the specific process of calculating the current unbalance amount of the line to be measured includes:
calculating power angles of measured data at two ends of the line to be measured;
based on a KCL principle, respectively calculating transition currents at two ends of a line to be detected;
and calculating the current unbalance amount of the current period according to the transition current amplitude.
Further, the specific process of determining the confidence interval endpoint is as follows:
calculating the average value of the current unbalance amount in N time periods
And the standard deviation S;
the confidence interval is determined as
And λ is the length of the confidence interval.
In a second aspect, the invention provides a system for removing SCADA bad data based on the amount of current unbalance of a power transmission line, the system comprises,
the data acquisition module is used for acquiring SCADA measurement data and ground-to-ground susceptance parameters at two ends of the line to be measured;
the data preprocessing module is used for converting the data on the two sides into single-phase measurement data;
the precision measurement module is used for setting the length of the confidence interval;
and the operation processing module is used for calculating the unbalance amount of the current of the line to be measured according to the single-phase measurement data, determining the end point of the confidence interval and rejecting the measurement data distributed outside the confidence interval.
Further, the system is in communication with a regulatory center.
The effect provided in the summary of the invention is only the effect of the embodiment, not all the effects of the invention, and one of the above technical solutions has the following advantages or beneficial effects:
by taking the current unbalance of the power transmission line as an evaluation index, bad data in the measured data can be accurately eliminated. The calculation process is simple, the accuracy is high, and the method can be applied to actual engineering. The accuracy of the measured data is improved, and the accuracy of programming and effect evaluation of schemes such as ordered power utilization, demand side response and the like is further improved. The method is also suitable for the Measurement data of a WAMS (Wide Area Measurement System) System aiming at the SCADA Measurement data, and has Wide application range.
Detailed Description
In order to clearly explain the technical features of the present invention, the following detailed description of the present invention is provided with reference to the accompanying drawings. The following disclosure provides many different embodiments, or examples, for implementing different features of the invention. To simplify the disclosure of the present invention, the components and arrangements of specific examples are described below. Furthermore, the present invention may repeat reference numerals and/or letters in the various examples. This repetition is for the purpose of simplicity and clarity and does not in itself dictate a relationship between the various embodiments and/or configurations discussed. It should be noted that the components illustrated in the figures are not necessarily drawn to scale. Descriptions of well-known components and processing techniques and procedures are omitted so as to not unnecessarily limit the invention.
The invention utilizes the voltage amplitude, the current amplitude, the active power and the susceptance parameter of the line at two ends of the line to obtain the current unbalance of the line, further obtain the confidence interval of the measured data based on the current unbalance, judge the effectiveness of the measured data and reject the bad measured data. The transition currents of the lines are calculated from two different line nodes. Based on that transition currents calculated by 2 different nodes on the same line are equal, when the difference between the transition currents and the difference is large, the effectiveness of the SCADA data at the section is low, and bad data are removed according to a confidence interval.
As shown in fig. 1, the method for removing the SCADA bad data based on the amount of current unbalance of the power transmission line of the present invention includes the following steps:
s1, obtaining SCADA measurement data at two ends of the line to be measured;
s2, converting the measured data into single-phase measured data;
s3, acquiring the ground-to-ground electric parameter of the line to be tested, and setting the length of a confidence interval;
and S4, calculating the unbalance amount of the line current to be measured according to the single-phase measurement data, determining the end point of the confidence interval, and eliminating the measurement data distributed outside the confidence interval.
In step S1, the acquired measurement data includes the voltage amplitude, the current amplitude, and the active power at the two ends of the line to be measured, and the measurement data is acquired once every fixed time period, for example, the data acquisition interval is set to 1 minute. The current flow is shown in figure 2. The embodiment comprises N time periods, the value of N can be determined according to the actual current transformation range, and if the electricity consumption demand is large in the first time period of a certain day, the corresponding current is large; if the power demand is small in the second time period, the corresponding current is small, the first time period and the second time period can be calculated separately, i.e. the operations of steps S1-S4 are performed separately, for example, N sampling periods are set in the first time period, taking the measurement data of the first time period as an example, the obtained measurement data are as shown in table 1 below,
U1(kV)
|
I1(A)
|
P1(MW)
|
U2 (kV)
|
I2(A)
|
P2 (MW)
|
304.414
|
479.743
|
405.096
|
303.894
|
468.756
|
-404.5 |
TABLE 1
In the context of table 1, the following,U 1(kV)、I 1(A) andP 1(MW) are phase voltage, phase current and three-phase active power of the side of the line 1 to be tested respectively;U 2(kV)、I 2(A) andP 2and (MW) are phase voltage, phase current and three-phase active power of the side of the line 2 to be tested respectively.
In step S2, the data obtained in step S1 is preprocessed, the SCADA data read from the actual power system is three-phase data, and the three-phase measurement data at both ends of the line are collectively converted into single-phase measurement data. The active power P is divided by 3 to convert the three phases into a single phase; and converting the measurement unit of the current from A to kA. Taking the measurement data in table 1 as an example, the measurement data after being preprocessed is shown in table 2.
U1(kV)
|
I1(kA)
|
P1(MW)
|
U2 (kV)
|
I2(kA)
|
P2 (MW)
|
304.414
|
0.480
|
135.032
|
303.894
|
0.469
|
-134.833 |
TABLE 2
In step S3, the ground-to-ground parameter of the current test line is obtainedB 1In this embodiment, the ground susceptance of the line to be tested isB 1 =7.377×10 -5Omega. While setting the length lambda of the confidence interval. The accuracy of the result is determined by the value of lambda, and the smaller the value of lambda is, the higher the accuracy is. In the present example, the length of the confidence interval λ = 5. The operation of step S3 may be performed before step S4.
In step S4, during actual operation, the voltage phase angle difference between the measured data at the two ends 1 and 2 of the line is very small, and is almost 0, by analyzing a large amount of SCADA measured data at the two ends of the line. Therefore, to simplify the calculation process, it is assumed that the 1-side voltage phase angle is 0 degrees and the 2-side voltage phase angle is 0 degrees.
Firstly, calculating the power angle of the measured data at two sides of the line to be measured. And calculating the power angle alpha of the measured data at the two ends of the line according to the voltage U, the current I and the active power P in the measured data. The calculation formula is as follows:
line 1 side power angle:
line 2-side power angle:
taking the first time period measurement data as an example, the power angle at two ends of the line is calculated according to the formulas (1) and (2)
。
The transition current of the line is then calculated. According to FIG. 2, as can be seen from the KCL (Kirchhoff laws, Kirchhoff's Law) principle,
the KCL equation is listed for the points of (i),
the KCL equation is listed for the point II,
whereinU 1、I 1The voltage and current amplitude, alpha, of the line 1 terminal1Is the power angle, beta, of the 1 end of the line1Is the terminal voltage phase angle of the line 1; whereinU 2、I 2The voltage and current amplitude, alpha, of the line 2 terminal2Power angle, beta, of line 2 terminal2Is the terminal voltage phase angle of line 2;jrepresenting the imaginary part of the complex number.
Assuming that the phase angle of the voltages at both ends of the line is 0, the equations (1) and (2) can be equivalent to the following equations (5) and (6),
wherein the content of the first and second substances,
is called as formed byA transition current with the 1 terminal directed to the 2 terminal,
referred to as a transition current from
terminal 2 to
terminal 1. Taking the measured data of the first time interval as an example, the current phasor obtained by calculation is,
and then calculating the unbalance amount of the current of the line to be measured. In the ideal case, the currents are equal in magnitude, i.e.I 12=I 21And the amplitude of the transition current flowing from the 1 terminal to the 2 terminal is equal to the amplitude of the transition current flowing from the 2 terminal to the 1 terminal. Due to the presence of the error(s),I 12andI 21with errors, called unbalance of currentE. The calculation formula is as follows:
using the measured data of the first time interval as the current unbalance amount
I.e. the amount of current imbalance for the first period of time
E=0.011kA。
Compute one by one
NThe current unbalance amount in each time period is calculated, and the average value and the standard deviation of the current unbalance amount in the first time period are calculated. According to the amount of current unbalance of each time period
E,
NAverage value of current unbalance amount in each period
Standard deviation S. The calculation formula is as follows,
get
N=448, current unbalanceAverage value of quantity
。
Finally, a calculation formula according to the confidence interval
Calculating the end point value of the confidence interval to obtain the confidence interval, wherein the confidence interval for measuring the current unbalance amount of the data is [0.3588A, 18.08A ] in the embodiment]. As shown in fig. 3, two horizontal lines parallel to the horizontal axis in the figure represent confidence intervals according to which 1 data point in this example is distributed outside the confidence interval, marked with an "O" in the figure. And the measured data corresponding to the two data points has poor validity and is removed.
In order to verify the using effect of the method, the following calculation example is designed to verify the effectiveness of the method. Assume that the line-to-ground capacitance parameter B1=7.377 × 10-5 Ω in fig. 2. The SCADA measurement data at two ends of a certain period of time is shown in the following table 3,
TABLE 3
In the measurement data of table 3, random errors are added at a certain time interval to form SCADA simulation measurement data with different measurement errors at multiple time intervals. The error types are: the active power of a certain time period at two ends of the line a has an error of 5%. As shown in the following table 4,
type of error
|
Period of bad data
|
Bad data current unbalance amount
|
a
|
3
|
-0.663 |
TABLE 4
As can be seen from the simulation results shown in fig. 4, the bad measurement data a can be effectively identified and marked with "O".
Compared with the existing bad data eliminating method, the method does not need to establish a complex equation, is simple and can quickly eliminate the bad data; the method is mainly applied to SCADA (supervisory control and data acquisition) measurement data, is also applicable to WAMS measurement data (the characteristics of the WAMS measurement data are similar to those of the SCADA measurement data, and the WAMS measurement data are higher than the SCADA measurement data, and the WAMS measurement data contain time synchronization, namely the WAMS measurement data contain angles), and has wide application range; through the analysis of the specific examples, the method can be applied to actual engineering, the accuracy of the measured data is improved, and references are provided for state estimation, load flow calculation and the like.
As shown in fig. 5, the SCADA bad data eliminating system based on the amount of current unbalance of the power transmission line of the present invention includes a data acquisition module, a data preprocessing module, an accuracy measurement module, and an operation processing module. The data acquisition module is used for acquiring SCADA measurement data and ground-to-ground sodium-capacity parameters at two ends of a line to be measured; the data preprocessing module is used for converting the data on the two sides into single-phase measurement data; the precision measurement module is used for setting the length of the confidence interval; and the operation processing module is used for calculating the unbalance amount of the current of the line to be measured according to the single-phase measurement data, determining the end point of the confidence interval and rejecting the measurement data distributed outside the confidence interval.
The system of the invention is communicated with the regulation and control center, and sends the SCADA measurement data with the bad data removed to the regulation and control center.
The foregoing is only a preferred embodiment of the present invention, and it will be apparent to those skilled in the art that various modifications and improvements can be made without departing from the principle of the invention, and such modifications and improvements are also considered to be within the scope of the invention.