CN103886193A - Fuzzy self-adaptation robust estimation method of electric power system - Google Patents
Fuzzy self-adaptation robust estimation method of electric power system Download PDFInfo
- Publication number
- CN103886193A CN103886193A CN201410093428.5A CN201410093428A CN103886193A CN 103886193 A CN103886193 A CN 103886193A CN 201410093428 A CN201410093428 A CN 201410093428A CN 103886193 A CN103886193 A CN 103886193A
- Authority
- CN
- China
- Prior art keywords
- partiald
- sigma
- lambda
- measuring point
- measuring
- 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.)
- Granted
Links
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a fuzzy self-adaptation robust estimation method of an electric power system, and provides a fuzzy self-adaptation robust estimation algorithm with uncertainty metering and measuring weight of the electric power system. According to the algorithm, a measuring point inferior fuzzy membership conception is put forward, the standard deviation of a measuring point is corrected in an on-line mode according to the membership of the measuring point, and gross error self adaptation is achieved. A modern interior point method is used for solving, convergence is good, results cannot be affected by an initial value easily, the solving large-scale optimization problem is solved, bad data are effectively recognized, residual error contamination and residual error flooding are avoided, excellent robust estimation performance is achieved to the quantity of state and quantity measuring, and a good engineering application background is achieved.
Description
Technical field
The present invention relates to a kind of electric system fuzzy self-adaption robust method of estimation, belong to power system monitoring, analysis and control technical field.
Background technology
State estimation is the basis of Power System Analysis and control, it is energy management system key foundation module, its main task is the real-time information providing according to data acquisition and monitoring system (SCADA), estimate the each busbar voltage of electrical network (amplitude, phase angle) and power, and comprise bad data detection and identification function.
Traditional state estimation algorithm is founded in 1970 by people such as Schweppe.According to the difference that solves objective function, state estimator comprises: weighted least-squares (WLS), Non quadratic criteria (QC, QL etc.), weighting least absolute value (WLAV), minimum median square (LMS), minimum truncated side (LTS) etc.At present most widely used is WLS estimation criterion, advantage is that model is simple, calculated amount is little, to desirable normal distribution measurement amount, can reach optimal estimation result, but easily being measured bad data, estimated result affects, need special bad data detection and identification program (as based on maximum weighted discrepancy principle) to reduce the impact of bad data on state estimation result, in the time there is extensive bad data in measurement, iterations is many, and may not necessarily effectively pick out all bad datas.The estimators such as QC, QL, WLAV, LMS, LTS all belong to robust and estimate category, its advantage is that state estimation algorithm is without extra bad data detection and identification program, but also exist need subjective to set that weighting factor, calculated amount are large, the defect such as robustness weakness under some particular case.
Domestic and international most state estimation algorithms are all that hypothesis measures variances sigma at present
2known, and in Practical Project, in large-scale electrical power system, there are numerous measuring instruments, measuring instrument is checked again and must pay very large economic cost, thereby be difficult to systematically again check measuring instrument; In addition along with the moment aging, running environment of equipment changes, the measurement accuracy of instrument is difficult to continue to keep stable.That is to say, it is time dependent measuring variance, is difficult to accurately estimate, measures weight and exists uncertain.
Summary of the invention
Goal of the invention: the present invention proposes a kind of electric system fuzzy self-adaption robust method of estimation, effectively identification bad data, avoided residual contamination and residual error to flood, have robust estimated performance capitally.
Technical scheme: the technical solution used in the present invention is a kind of electric system fuzzy self-adaption robust method of estimation, comprises the following steps:
1) fuzzy membership function of i measuring point poor quality of definition is:
In formula: r
ifor the residual error of measuring point i, σ
ifor the measuring standard of measuring point i poor, v
i(| r
i|, σ
i) be the fuzzy membership of measuring point i, a, b is greater than 0 constant;
2) taking the Weighted Fuzzy degree of membership sum that minimizes measuring point poor quality as optimization aim, following Optimized model is proposed:
3) consider to measure weight uncertainty, i.e. the poor uncertainty of measuring standard, poor based on measuring point fuzzy membership correction measuring standard, to realize the self-adaptation to measuring rough error, for the k+1 time iteration, order:
4) introduce non-negative relaxation factor l, u, step 2) in Optimized model can be equivalent to:
5) equality constraint in step 3) is made as to barrier function, can be able to lower Lagrangian function:
In formula: λ, α, β are m dimension Lagrange multiplier, i.e. dual variable; L, u is former variable;
6) solution procedure 4) in the KKT condition of Lagrangian function obtain:
In above formula, ▽ h (x) is the Jacobian matrix of h (x), L=diag (l
1, l
2..., l
m), U=diag (u
1, u
2, u
m); A=diag (α
1, α
2..., α
m), B=diag (β
1, β
2..., β
m), e=[1 ... 1]
t, μ is the disturbance factor, and meets
8) will in step 5), for KKT condition, after newton-La Fusenfa linearization, obtain following update equation:
In formula: ▽
2h (x) is the gloomy matrix in sea of h (x);
8), according to the iterative step of former-dual interior point, upgrade quantity of state, until convergence.
Beneficial effect: the present invention proposes one and takes into account measurement weight probabilistic electric system fuzzy self-adaption robust algorithm for estimating (FARE), this algorithm has proposed the fuzzy membership concept of measuring point poor quality, and revise online the standard deviation of measuring point according to the degree of membership of measuring point, realize the self-adaptation to rough error.Solve with modern interior-point method, convergence is good, and result is not subject to initial value affecting, is applicable to solve Large-scale Optimization Problems.Effectively identification bad data, avoided residual contamination and residual error to flood, quantity of state and measurement amount are had to robust estimated performance capitally, there is engineering application background well.
Brief description of the drawings
Fig. 1 is workflow diagram of the present invention;
Fig. 2 is taking IEEE14 as standard testing node, the test design sketch of the inventive method FARE and WLAV, QC, WLS algorithm measurement amount robustness, and the weighted residual average of measurement | r/ σ | value distribution plan;
Fig. 3 is taking IEEE118 as standard testing node, the test design sketch of the inventive method FARE and WLAV, QC, WLS algorithm measurement amount robustness, and the weighted residual average of measurement | r/ σ | value distribution plan.
Embodiment
Below in conjunction with the drawings and specific embodiments, further illustrate the present invention, should understand these embodiment is only not used in and limits the scope of the invention for the present invention is described, after having read the present invention, those skilled in the art all fall within the application's claims limited range to the amendment of various equivalents of the present invention.
1. set up electric system fuzzy self-adaption robust estimation model
Power system state estimation measurement equation is:
z=h(x)+e
In formula: x is state vector, the system nodes of setting up departments is n, amplitude and phase angle that x comprises node voltage, and its dimension is 2n-1; Z is that m dimension measures vector; H (x) is measurement function vectors; E is error in measurement vector.
The residual equation of state estimation is:
For i measuring point, if it estimates residual error r
ivery little, can think that this measuring point is high-quality measuring point; Otherwise, if it estimates residual error r
ivery large, can think that this measuring point is measuring point inferior.Because the measurement accuracy of each measuring point is not quite similar, definition σ
ifor the measuring standard of measuring point i poor, thereby with the weighted residual r of measuring point
i/ σ
ithe good and bad measuring point of size discrimination more reasonable.
The quality that definition event y is measuring point, if set Y has comprised all event y, element y and the available fundamental function of relation of gathering Y---membership function v (y) represents, for classical data acquisition theory, has so:
But the quality of actual measuring point is relative concept, there are not absolute superiority and inferiority, than classical data acquisition theory, fuzzy set allows degree of membership to get [0,1] any value on, choose continuously differentiable bell membership function herein, for i measuring point, the fuzzy membership function of its poor quality is:
In formula: a, b are greater than 0 fuzzy membership characteristic parameter.
Taking weighted least-squares (WLS) as main Length Factor Method in Power System State is taking minimum weight residual sum of squares (RSS) as optimization aim, similarly, fuzzy membership based on measuring point poor quality herein, taking the Weighted Fuzzy degree of membership sum that minimizes measuring point poor quality as optimization aim, following Optimized model is proposed:
Domestic and international most state estimation algorithms are all that hypothesis measures variances sigma at present
2known, and in Practical Project, in large-scale electrical power system, there are numerous measuring instruments, measuring instrument is checked again and must pay very large economic cost, thereby be difficult to systematically again check measuring instrument; In addition due to the increase of service time, the external environment condition that the moment changes, the measurement accuracy of instrument is difficult to continue to keep stable.That is to say, measure variance and be difficult to accurately estimate, measure weight and exist uncertain.The general subjectivity of state estimation of domestic real system arranges the weight factor of each measurement, debugging and safeguard complicatedly, and the subjective weight factor of setting may not measure variance with reality and coincide.
For this reason, poor based on measuring point fuzzy membership correction measuring standard herein, for the k+1 time iteration in solving-optimizing object procedure, order:
The poor correction function f of measuring standard (v) should follow 2 principles: 1. as measuring point fuzzy membership v
iapproach at 0 o'clock, now this measuring point is high-quality measuring point, should keep σ
iapproximate constant; 2. as measuring point fuzzy membership v
iapproach at 1 o'clock, now this measuring point is measuring point inferior, should increase σ
i, while occurring measuring rough error, should reduce the weight of this measuring point in iteration, reduce its impact on state estimation result.
Based on above-mentioned criterion, can choose correction function f (v)=[1-(1/v)]
1/b.
2. solving of the FARE model based on former dual interior point
With reference to Fig. 1, first the present invention obtains network parameter, topological structure, the measurement parameters of electric system, network parameter comprises: bus numbering, the reactance of node shunt capacitance, branch road head end numbering and end numbering, the reactance of branch road resistance, reactance, charging capacitor over the ground, transformer resistance, reactance, standard no-load voltage ratio; Topological structure mainly comprises the on off state connecting between two electrical equipments; Measurement parameters comprises: node injects meritorious, idle measurement, and branch road head end and end are meritorious, idle measurement, and busbar voltage amplitude measures.
Obtaining after the data of state estimation program needs, carry out program initialization, comprising: quantity of state voltage magnitude and phase angle, Lagrange multiplier, former dual relaxation variable, former antithesis Center Parameter, maximum iteration time and convergence precision.And form bus admittance matrix.
Introduce the non-negative relaxation factor l of m dimension, u, optimization aim can be equivalent to:
Equality constraint in above formula is made as to barrier function, can be able to lower Lagrangian function:
In formula: λ, α, β are m dimension Lagrange multiplier, i.e. dual variable; L, u is former variable.
Solving its KKT condition can obtain:
L
x=▽h(x)λ=0
L
λ=z-h(x)+l-u=0
In above formula, ▽ h (x) is the Jacobian matrix of h (x), L=diag (l
1, l
2..., l
m), U=diag (u
1, u
2, u
m); A=diag (α
1, α
2..., α
m), B=diag (β
1, β
2..., β
m), e=[1 ... 1]
t, μ is the disturbance factor, and meets
KKT condition is Nonlinear System of Equations, and available newton-La Fusenfa solves, and will after its linearization, obtain update equation group:
▽
2h(x)λdx+▽h(x)dλ=-L
x
-▽
Th(x)dx+dl-du=-L
λ
In formula: ▽
2h (x) is the gloomy matrix in sea of h (x).Being write as matrix form can obtain:
In formula:
And can calculate [dl du dx] by above formula
t, [d α d β d λ]
t.
Calculate duality gap Gap=α
tl+ β
tu and disturbance factor mu=δ Gap/2m.In addition, permanent in non-negative in order to ensure relaxation factor, the iteration step length (θ of, dual variable former according to following formula
r, θ
d):
Revise former, dual variable:
Correction measuring standard is poor: σ
(k+1)=σ
(k)(1/v)
1/b.
Repeat above-mentioned makeover process until duality gap Gap reaches the maximum iteration time that convergence precision or iterations exceed setting.
3. sample calculation analysis
Test example of the present invention comprises IEEE14, IEEE30, IEEE57, IEEE118 node, Poland 2383,2746 nodes are (with WP-2383, WP-2746 represents), metric data adds random error in measurement by strict trend result and obtains, and considers the poor uncertainty of measuring standard, the poor obedience [0.002 of voltage measuring standard, being uniformly distributed 0.005], power measurement standard deviation is obeyed being uniformly distributed of [0.004,0.01].Measure redundance
between 3.5~4.5.It is poor that bad data is added and subtracted [5,20] maximum measuring standard doubly on the basis of metric data, and add at random 10% bad data.
The superiority of FARE algorithm proposing for checking the present invention, by FARE with proposed in the past state estimation algorithm, comprise WLS, WLAV, QC, carry out comparative study from aspects such as bad data identification, robustness, constringency performances.The voltage of WLS, WLAV, QC algorithm measures weight and gets 4, and meritorious, idle measurement weight gets 1.For the poor initial value σ of measuring standard of FARE algorithm
(0), voltage measures 0.001, and meritorious, idle measurement gets 0.002; Fuzzy membership parameter is got a=2.5, b=3.
3.1 measurement amount robustness
For the bad data identification of more various algorithms and resist the performance that residual contamination floods, measuring true value
the poor true value of measuring standard
in known situation, definition:
In formula:
be the estimated value of i measuring point, index τ has weighed measurement estimated value and has departed from the degree that measures true value.
Taking IEEE14, IEEE118 as test example, under 4 kinds of state estimation algorithms, the percent profile of τ is shown in Figure of description 2, Fig. 3.
From Fig. 2, Fig. 3, the number percent of τ > 3 under the estimation of WLS algorithm is more than 20%, be only 10% and measure bad data ratio, this explanation has occurred that very serious residual contamination floods phenomenon, thereby QC algorithm based on WLS is also difficult to effective identification bad data and resists residual contamination flood.Than WLAV, QC, WLS algorithm, FARE algorithm estimates that the distribution of lower τ mainly concentrates between 0~1, and the ratio of τ > 3 is obviously less than all the other 3 kinds of algorithms, thereby measurement estimated value under FARE algorithm is closer to measuring true value, has bad data identification well and resists residual contamination and flood performance.
3.2 quantity of state robustness
In order to compare the quantity of state robustness of 4 kinds of algorithms, adopt following 2 indexs to evaluate the estimated performance of different conditions algorithm for estimating to quantity of state:
1) square error of voltage magnitude:
In formula: V
ifor the true value of node i voltage magnitude,
for the estimated value of node i voltage magnitude.
2) square error of voltage phase angle:
In formula: θ
ifor the true value of node i voltage phase angle,
for the estimated value of node i phase angle, acquiescence node 1 is balance node.
At 6 service systems, all containing 10% bad data in the situation that, table 1 has provided the MSE1 of 4 kinds of state estimation algorithms, the index of MSE2.
The size of MSE1 and MSE2 in comparison sheet 1, the estimated performance of known each algorithm to quantity of state: FARE>WLAV>QC>W LS, and FARE is obviously better than all the other 3 kinds of algorithms.In addition, FARE algorithm is better than the estimated performance to voltage magnitude to the estimated performance of voltage phase angle.
The robustness of comprehensive measurement amount and quantity of state is known, and in 4 kinds of state estimation algorithms choosing, the FARE that the present invention proposes is top quality Length Factor Method in Power System State.
Table 1: containing MSE1 and the MSE2 index performance of different algorithm for estimating under 10% bad data
Claims (1)
1. an electric system fuzzy self-adaption robust method of estimation, is characterized in that, comprises the following steps:
1) fuzzy membership function of i measuring point poor quality of definition is:
In formula: r
ifor the residual error of measuring point i, σ
ifor the measuring standard of measuring point i poor, v
i(| r
i|, σ
i) be the fuzzy membership of measuring point i, a, b is greater than 0 constant;
2) taking the Weighted Fuzzy degree of membership sum that minimizes measuring point poor quality as optimization aim, following Optimized model is proposed:
3) consider to measure weight uncertainty, i.e. the poor uncertainty of measuring standard, poor based on measuring point fuzzy membership correction measuring standard, to realize the self-adaptation to measuring rough error, for the k+1 time iteration, order:
4) introduce non-negative relaxation factor l, u, step 2) in Optimized model can be equivalent to:
5) equality constraint in step 3) is made as to barrier function, can be able to lower Lagrangian function:
In formula: λ, α, β are m dimension Lagrange multiplier, i.e. dual variable; L, u is former variable;
6) solution procedure 4) in the KKT condition of Lagrangian function obtain:
In above formula, ▽ h (x) is the Jacobian matrix of h (x), L=diag (l
1, l
2..., l
m), U=diag (u
1, u
2, u
m); A=diag (α
1, α
2..., α
m), B=diag (β
1, β
2..., β
m), e=[1 ... 1]
t, μ is the disturbance factor, and meets
7) will in step 5), for KKT condition, after newton-La Fusenfa linearization, obtain following update equation:
In formula: ▽
2h (x) is the gloomy matrix in sea of h (x);
8), according to the iterative step of former-dual interior point, upgrade quantity of state, until convergence.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410093428.5A CN103886193B (en) | 2014-03-13 | 2014-03-13 | Fuzzy self-adaptation robust estimation method of electric power system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410093428.5A CN103886193B (en) | 2014-03-13 | 2014-03-13 | Fuzzy self-adaptation robust estimation method of electric power system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103886193A true CN103886193A (en) | 2014-06-25 |
CN103886193B CN103886193B (en) | 2017-05-24 |
Family
ID=50955083
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410093428.5A Active CN103886193B (en) | 2014-03-13 | 2014-03-13 | Fuzzy self-adaptation robust estimation method of electric power system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103886193B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104795811A (en) * | 2015-03-17 | 2015-07-22 | 河海大学 | Power system interval state estimation method |
CN105187334A (en) * | 2015-06-25 | 2015-12-23 | 桂林电子科技大学 | Channel estimation method for power line communication affected by impact noise |
CN105914738A (en) * | 2016-05-26 | 2016-08-31 | 国网山东省电力公司潍坊供电公司 | Power distribution network bad data detection and identification method based on uncertainty of measurement |
CN107947160A (en) * | 2017-11-28 | 2018-04-20 | 国家电网公司 | A kind of bilinearity WLAV method for estimating state based on matrixing |
CN109255541A (en) * | 2018-09-20 | 2019-01-22 | 华北电力大学 | A kind of power distribution network robust state estimation method the sum of multiplied based on least square and one |
CN110428185A (en) * | 2019-08-08 | 2019-11-08 | 河海大学 | Electric-thermal based on pseudo- measurement model interconnects integrated energy system robust state estimation method |
CN110739687A (en) * | 2019-10-24 | 2020-01-31 | 福州大学 | electric power system distribution robust scheduling method considering wind power high-order uncertainty |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109087014A (en) * | 2018-08-13 | 2018-12-25 | 河海大学 | A kind of electric system Robust filter method |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101599643B (en) * | 2009-04-23 | 2011-04-27 | 清华大学 | Robust state estimation method in electric power system based on exponential type objective function |
CN103279676B (en) * | 2013-06-07 | 2016-08-31 | 河海大学 | A kind of power system WLAV Robust filter method based on substitution of variable |
CN103336904A (en) * | 2013-07-08 | 2013-10-02 | 国家电网公司 | Robust state estimation method based on piecewise linearity weight factor function |
-
2014
- 2014-03-13 CN CN201410093428.5A patent/CN103886193B/en active Active
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104795811A (en) * | 2015-03-17 | 2015-07-22 | 河海大学 | Power system interval state estimation method |
CN104795811B (en) * | 2015-03-17 | 2017-07-18 | 河海大学 | A kind of power system state of section method of estimation |
CN105187334A (en) * | 2015-06-25 | 2015-12-23 | 桂林电子科技大学 | Channel estimation method for power line communication affected by impact noise |
CN105187334B (en) * | 2015-06-25 | 2018-04-17 | 桂林电子科技大学 | A kind of power line communication channel method of estimation for the noise jamming that is hit |
CN105914738A (en) * | 2016-05-26 | 2016-08-31 | 国网山东省电力公司潍坊供电公司 | Power distribution network bad data detection and identification method based on uncertainty of measurement |
CN107947160A (en) * | 2017-11-28 | 2018-04-20 | 国家电网公司 | A kind of bilinearity WLAV method for estimating state based on matrixing |
CN109255541A (en) * | 2018-09-20 | 2019-01-22 | 华北电力大学 | A kind of power distribution network robust state estimation method the sum of multiplied based on least square and one |
CN110428185A (en) * | 2019-08-08 | 2019-11-08 | 河海大学 | Electric-thermal based on pseudo- measurement model interconnects integrated energy system robust state estimation method |
CN110739687A (en) * | 2019-10-24 | 2020-01-31 | 福州大学 | electric power system distribution robust scheduling method considering wind power high-order uncertainty |
CN110739687B (en) * | 2019-10-24 | 2021-06-01 | 福州大学 | Power system distribution robust scheduling method considering wind power high-order uncertainty |
Also Published As
Publication number | Publication date |
---|---|
CN103886193B (en) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103886193A (en) | Fuzzy self-adaptation robust estimation method of electric power system | |
Schenato et al. | Bayesian linear state estimation using smart meters and pmus measurements in distribution grids | |
CN104134999B (en) | Distribution network based on multi-data source measures the practical method of calculation of efficiency analysis | |
CN103840452B (en) | A kind of bulk power grid method for estimating state introducing PMU measurement information | |
CN107016236B (en) | Power grid false data injection attack detection method based on nonlinear measurement equation | |
CN103972884B (en) | A kind of power system state estimation method | |
Clements | The impact of pseudo-measurements on state estimator accuracy | |
CN107016489A (en) | A kind of electric power system robust state estimation method and device | |
CN101969198B (en) | Method for estimating electrical power system state with consideration of load static property | |
Sun et al. | Power system observability and dynamic state estimation for stability monitoring using synchrophasor measurements | |
CN106707061A (en) | Hybrid measurement based power distribution network dynamic state estimation method | |
CN101599643A (en) | A kind of anti-difference of electric power system method for estimating state based on the exponential type target function | |
US20140244189A1 (en) | System For Achieving Real-Time Monitoring and State Estimation in Power Distribution Networks | |
CN104836223A (en) | Power grid parameter error and bad data coordinated identification and estimation method | |
CN103020726B (en) | Robust state estimation method for full PMU measurement | |
CN103944165A (en) | Large power grid parameter recognition and estimation method | |
CN104992010A (en) | Topologic partition based multi-section joint parameter estimation method | |
CN106026086B (en) | The method for dynamic estimation of operating status in a kind of power grid | |
CN107453484A (en) | A kind of SCADA data calibration method based on WAMS information | |
Mínguez et al. | State estimation sensitivity analysis | |
CN103034787A (en) | Method for estimating state of microgrid | |
CN114204560A (en) | Medium voltage distribution network line parameter identification method | |
CN104866715A (en) | Electrical power system robust state estimation method based on self-adaptive kernel density estimation | |
CN102280877B (en) | Method for identifying parameter of poor branch of power system through a plurality of measured sections | |
CN106372440B (en) | A kind of adaptive robust state estimation method of the power distribution network of parallel computation and device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |