Summary of the invention
Provided in the embodiment of the present invention SCADA bad data elimination method based on transmission line of electricity current imbalance amount and
System, to solve the problems, such as existing method, the practicability is poor, error is big.
In order to solve the above-mentioned technical problem, the embodiment of the invention discloses following technical solutions:
First aspect present invention provides the SCADA bad data elimination method based on transmission line of electricity current imbalance amount,
The following steps are included:
S1 obtains the SCADA metric data at route both ends to be measured;
The metric data is converted to single-phase metric data by S2;
S3, obtains the susceptance parameter over the ground of route to be measured, and the length of confidence interval is arranged;
S4 calculates the amount of unbalance of line current to be measured according to the single-phase metric data, and determines the confidence interval
Endpoint, by the metric data being distributed in except confidence interval reject.
Further, the acquisition of a metric data is carried out every the fixed period in step S1;Confidence described in step S4
Section is used to detect the metric data obtained in current all periods.
Further, the SCADA metric data includes voltage magnitude, current amplitude and the wattful power at route both ends to be measured
Rate.
Further, the detailed process for calculating line current amount of unbalance to be measured are as follows:
Calculate the power angle of route both ends to be measured metric data;
Based on KCL principle, the transition current at route both ends to be measured is calculated separately;
According to transition current amplitude, the current imbalance amount of present period is calculated.
Further, it is determined that the detailed process of the confidence interval endpoint are as follows:
Calculate the average value of N number of period current imbalance amountWith standard deviation S;
Confidence interval is determined asλ is the length of the confidence interval.
Second aspect of the present invention provides the SCADA bad data eliminating system based on transmission line of electricity current imbalance amount,
The system comprises,
Data acquisition module, for obtaining the SCADA metric data and susceptance parameter over the ground at route both ends to be measured;
Data preprocessing module, for the two sides data to be converted to single-phase metric data;
Accuracy measure module, for the length of confidence interval to be arranged;
Calculation process module calculates the amount of unbalance of line current to be measured according to the single-phase metric data, and determines institute
The endpoint for stating confidence interval rejects the metric data being distributed in except confidence interval.
Further, the system is communicated with regulation center.
Specific embodiment
In order to clarify the technical characteristics of the invention, below by specific embodiment, and its attached drawing is combined, to this hair
It is bright to be described in detail.Following disclosure provides many different embodiments or example is used to realize different knots of the invention
Structure.In order to simplify disclosure of the invention, hereinafter the component of specific examples and setting are described.In addition, the present invention can be with
Repeat reference numerals and/or letter in different examples.This repetition is that for purposes of simplicity and clarity, itself is not indicated
Relationship between various embodiments and/or setting is discussed.It should be noted that illustrated component is not necessarily to scale in the accompanying drawings
It draws.Present invention omits the descriptions to known assemblies and treatment technology and process to avoid the present invention is unnecessarily limiting.
The present invention is acquired using the susceptance parameter of the voltage magnitude at route both ends, current amplitude, active power and route
The current imbalance amount of route, and then the confidence interval of the metric data based on current imbalance amount is acquired, judge metric data
Validity, reject bad metric data.The transition current of route is calculated from two different line nodes respectively.It is based on
It is equal by 2 calculated transition currents of different nodes of same route, when gap is larger between the two, illustrate this section
SCADA data validity is relatively low, then is rejected bad data according to confidence interval.
As shown in Figure 1, the present invention is based on the SCADA bad data elimination methods of transmission line of electricity current imbalance amount, including
Following steps:
S1 obtains the SCADA metric data at route both ends to be measured;
The metric data is converted to single-phase metric data by S2;
S3, obtains the susceptance parameter over the ground of route to be measured, and the length of confidence interval is arranged;
S4 calculates the amount of unbalance of line current to be measured according to the single-phase metric data, and determines the confidence interval
Endpoint, by the metric data being distributed in except confidence interval reject.
In step S1, the metric data of acquisition includes voltage magnitude, current amplitude and the active power at route both ends to be measured,
And each fixed period acquires a metric data, and it is 1 minute that data acquisition intervals, which are such as arranged,.Current direction is as shown in Figure 2.
It include N number of period in the present embodiment, the value of N can be determined according to actual current transformation range, if in certain day first time period
Interior, electricity demand is big, then corresponding current is big;Electricity demand is small in second time period, then corresponding current is small, when can be by first
Between section and second time period separate computations, i.e., respectively carry out step S1~S4 operation, N number of adopt such as is set in first time period
The sample period, by taking the metric data of the first period as an example, the metric data got is 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 table 1, U1(kV)、I1(A) and P1(MW) be respectively 1 side of route to be measured phase voltage, phase current and three-phase wattful power
Rate;U2(kV)、I2(A) and P2(MW) be respectively 2 side of route to be measured phase voltage, phase current and three phases active power.
In step S2, the step S1 data obtained are pre-processed, the SCADA data read in practical power systems
For three-phase data, route both ends three-phase metric data is uniformly converted into single-phase metric data.By active-power P divided by 3
It is single-phase three-phase inversion;The measurement unit of electric current is converted into kA by A simultaneously.By taking the metric data in table 1 as an example, by pre-
Treated, and metric data is as 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 parameter of the susceptance over the ground B of current measurement circuit is obtained1, in the present embodiment, route to be measured it is electric over the ground
It receives as B1=7.377 × 10-5Ω.The length lambda of confidence interval is set simultaneously.The size of λ value determines the precision of result, λ's
Value is smaller, and precision is higher.In this example, length lambda=5 of confidence interval.As long as the operating process of step S3 is in step S4
It carries out before.
In step S4, when practical operation, pass through the analysis to a large amount of route both ends SCADA metric data, route both ends 1,2
The voltage phase angle difference very little of side metric data, almost 0.Therefore, process to simplify the calculation, it is assumed that 1 side voltage phase angle is 0
Degree, 2 side voltage phase angles are 0 degree.
The power angle of route two sides to be measured metric data is calculated first.According to the voltage U, electric current I, active in metric data
The power angle α of power P calculating route both ends metric data.Its calculation formula is:
1 side power angle of route:
2 side power angle of route:
By taking the first period metric data as an example, gained is calculated according to formula (1) (2), the power angle at route both ends is ∠ α1=-
0.391°;∠α2=-2.813 °.
Then the transition current of route is calculated.According to fig. 2, by KCL (Kirchhoff laws, Kirchhoff's law) principle
It is found that
KCL equation is listed to 1. point to obtain,
KCL equation is listed to 2. point to obtain,
Wherein U1、I1Voltage, current amplitude for 1 end of route, α1For the power angle at 1 end of route, β1For 1 end voltage of route
Phase angle;Wherein U2、I2Voltage, current amplitude for 2 end of route, α2The power angle at 2 end of route, β2For 2 end voltage phase angle of route;j
Indicate the imaginary part in plural number.
Assuming that route both end voltage phase angle is 0, then (1) (2) formula can be equivalent to following (5) (6) formula
Wherein,The transition current at 2 ends is referred to as directed toward by 1 end,The transition current at 1 end is referred to as directed toward by 2 ends.With
For the metric data of one period, calculating gained electric current phasor is,
Then the amount of unbalance of line current to be measured is calculated.In the ideal case, current amplitude is equal, i.e. I12=I21, by
The transition current amplitude that 1 end flows to 2 ends is equal to the transition current amplitude that 1 end is flowed to by 2 ends.Due to the presence of error, I12With I21
There are errors, referred to as the amount of unbalance E of electric current.Its calculation formula is:
E=| I12-I21| (7)
With the metric data of the first period, current imbalance amount E=| 0.480-0.469 |=0.011, i.e. the first period
Current imbalance amount E=0.011kA.
Calculate N number of period current imbalance amount one by one, and calculate the average value of current imbalance amount in first time period with
And standard deviation.According to the current imbalance amount E of each period, the average value of N number of period current imbalance amountStandard deviation S.Its
Calculation formula is
Take N=448, the average value of current imbalance amount
Finally, according to the calculation formula of confidence intervalThe endpoint value of confidence interval is calculated, thus
To confidence interval, in the present embodiment, the confidence interval of metric data current imbalance amount is [0.3588A, 18.08A].Such as Fig. 3
Shown, two horizontal lines parallel with horizontal axis indicate that confidence interval has 1 number according to the confidence interval in the present embodiment in figure
Strong point is distributed in except confidence interval, is marked in figure with " O ".The corresponding metric data validity of two data points compared with
Difference rejects it.
For the using effect for verifying the above method, the validity of following Example Verification the method for the present invention is designed.Assuming that attached drawing 2
Middle line-to-ground susceptance parameter B1=7.377 × 10-5Ω。
A certain period route both ends SCADA metric data is as shown in table 3 below,
U1(kV) |
I1(A) |
U2(kV) |
I2(A) |
P1(MW) |
P2(MW) |
304.431 |
437.047 |
303.903 |
439.188 |
133.050 |
-132.911 |
Table 3
In the metric data of table 3, random error is added in a certain period, forms the multiple periods for having different errors in measurement
SCADA simulating measuring data.Error pattern are as follows: the active power of a route both ends a certain period has 5% error.Such as the following table 4
It is shown
Error pattern |
Period where bad data |
Bad data current imbalance amount |
a |
3 |
-0.663 |
Table 4
By the simulation result of Fig. 4 it is found that bad metric data a can be identified effectively, marked with " O ".
Compared with existing bad data elimination method, this method does not need to establish complicated equation, and method is simple, can be fast
Speed rejects bad data;This method is mainly used in SCADA metric data, and being equally applicable to WAMS metric data, (WAMS is measured
Data characteristics is characterized in similar with SCADA metric data, and WAMS metric data ratio SCADA metric data is more advanced.WAMS amount
Measured data clock synchronization containing having time, that is, WAMS metric data is containing angled), it has wide range of applications;By to specific example
Analysis improves the accuracy of metric data it is found that this method can be applied in Practical Project, to carry out state estimation, trend meter
The offers references such as calculation.
As shown in figure 5, including number the present invention is based on the SCADA bad data eliminating system of transmission line of electricity current imbalance amount
According to acquisition module, data preprocessing module, accuracy measure module and calculation process module.Data acquisition module is to be measured for obtaining
The SCADA metric data at route both ends and over the ground susceptance parameter;Data preprocessing module is for being converted to the two sides data
Single-phase metric data;Accuracy measure module is used to be arranged the length of confidence interval;Calculation process module, according to the single-phase measurement
Data calculate the amount of unbalance of line current to be measured, and determine the endpoint of the confidence interval, will be distributed in except confidence interval
Metric data reject.
System of the invention is communicated with regulation center, and the SCADA metric data for weeding out bad data is sent to
Regulation center.
The above is the preferred embodiment of the present invention, for those skilled in the art,
Without departing from the principles of the invention, several improvements and modifications can also be made, these improvements and modifications are also regarded as this hair
Bright protection scope.