CN107016489A - A kind of electric power system robust state estimation method and device - Google Patents
A kind of electric power system robust state estimation method and device Download PDFInfo
- Publication number
- CN107016489A CN107016489A CN201710137890.4A CN201710137890A CN107016489A CN 107016489 A CN107016489 A CN 107016489A CN 201710137890 A CN201710137890 A CN 201710137890A CN 107016489 A CN107016489 A CN 107016489A
- Authority
- CN
- China
- Prior art keywords
- node
- measurement
- power
- state variable
- power system
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000005259 measurement Methods 0.000 claims abstract description 226
- 239000011159 matrix material Substances 0.000 claims abstract description 97
- 238000002347 injection Methods 0.000 claims abstract description 44
- 239000007924 injection Substances 0.000 claims abstract description 44
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000001514 detection method Methods 0.000 claims abstract description 12
- 230000008859 change Effects 0.000 claims abstract description 5
- 238000012937 correction Methods 0.000 claims description 22
- 230000017105 transposition Effects 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 30
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 4
- 230000005611 electricity Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005304 joining Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- JEIPFZHSYJVQDO-UHFFFAOYSA-N ferric oxide Chemical compound O=[Fe]O[Fe]=O JEIPFZHSYJVQDO-UHFFFAOYSA-N 0.000 description 1
- 230000009931 harmful effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Power Engineering (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The embodiment of the present invention provides a kind of electric power system robust state estimation method and device based on piecewise nonlinear weight function, and method includes building the bus admittance matrix of power system, measurement measurement equation and zero injection equality constraint;The Jacobian matrix and measurement residual error of power system are calculated, the weight of each measurement is obtained;State variable is updated, and judges whether iterative algorithm reaches convergence threshold, if then obtaining the node voltage amplitude and phase of final power system.The present invention carries out bad data detection and identification and state estimation simultaneously, while state estimation is calculated, can change measurement weight, bad data weight is constantly reduced in iteration, suppress the influence that bad data is brought, effectively reduce amount of calculation, improve operating efficiency;And directly bring the residual error of measurement into segmentation nonlinear weight function, it is to avoid and the calculating of standardized residual, calculating speed can be greatly promoted, the calculating time is saved, is very suitable for engineer applied.
Description
Technical field
The present invention relates to Power system state estimation technology, and in particular to a kind of electric power based on piecewise nonlinear weight function
System robust state estimation method and device.
Background technology
What traditional Power system state estimation was most widely used is weighted least square algorithm, the letter of its model
Single, computational methods are convenient, and it is minimum and unbiased to obtain variance in the case where error in measurement obeys the hypothesis of preferable normal distribution
Estimated result.
But power system measurement in, in addition to random error and systematic error, also inevitably exist some with
The bad data that true value differs greatly, these bad datas can bring serious harmful effect to estimated result.It is mixed in measurement number
Bad data in can make the estimated result of weighted least-squares method lose reliability and accuracy.Estimate to improve presence
The accuracy and practicality of result are counted, is typically the detection that bad data is added after state estimation procedure in actual applications
With identification link, go to examine bad data using methods such as residual tests.However, this mode needs state estimation and umber of defectives
According to the multiple circular test of two links of detection and identification, its workload is huge, and the calculating time is long, and due to residual test method
Itself intrinsic defect, identification reliability is not high enough.
Existing robust state estimation method includes weighting least absolute value (Weighted least absolute
Absolute, WLAV) and weighting least absolute value robust state estimation based on substitution of variable interior point method etc., its solution procedure is answered
It is miscellaneous, it is computationally intensive, it is impossible to be applied in the real-time estimation of large-scale power system.
The content of the invention
In order to overcome the above-mentioned deficiencies of the prior art, the present invention provides a kind of power train based on piecewise nonlinear weight function
Robust state estimation method of uniting and device, gather network parameter and SCADA system parameter first, and according to the network parameter of collection
Bus admittance matrix, measurement measurement equation and zero injection equality constraint are determined with SCADA system parameter;Then with node electricity
Pressure amplitude value and phase are that state variable calculates Jacobian matrix and measurement residual error, and are obtained respectively by piecewise nonlinear weight function
The weight of measurement;State variable is updated finally by Newton iteration method, and carries out convergence judgement, and then obtains state estimation knot
Really.The present invention can meet computational efficiency requirement, but also with good robustness.
In order to realize foregoing invention purpose, the present invention is adopted the following technical scheme that:
The present invention provides a kind of electric power system robust state estimation method based on piecewise nonlinear weight function, methods described
Including:
According to the network parameter of data acquisition and supervisor control the SCADA power system collected and node and branch
Road operational factor, builds the bus admittance matrix, measurement measurement equation and zero injection equality constraint of power system;
Node voltage amplitude and phase using power system is as state variable, and the node of the power system according to structure
The Jacobian matrix and measurement that admittance matrix, measurement measurement equation and zero injection equality constraint calculate power system are residual
Difference, and obtain by piecewise nonlinear weight function the weight of each measurement;
State variable is updated by iterative algorithm, and judges whether iterative algorithm reaches receipts according to the correction of state variable
Threshold value is held back, if then obtaining the node voltage amplitude and phase of final power system.
The network parameter and SCADA system parameter according to collection determines bus admittance matrix, measurement equation and zero
Injection equality constraint includes:
The bus admittance matrix is expressed as:
Wherein, YiiRepresent the diagonal entry of bus admittance matrix, YijRepresent the non-diagonal line element of bus admittance matrix
Element, yijRepresent the branch impedance z between node i and node jijInverse.
The network parameter and SCADA system parameter according to collection determines bus admittance matrix, measurement measurement equation
And zero injection equality constraint include:
The measurement measurement equation is expressed as:
Z=h (x)+r
Wherein, z represents measurement measuring value, and r represents measurement residual error, and x represents state variable, and h (x) calculates for measurement
Equation, it includes node power measurement equation, circuit branch road power measurement equation and transformer branch power measurement equation;
Wherein, node power measurement equation is expressed as:
Wherein, PiThe active power injected for node i, QiThe reactive power injected for node i, viFor the voltage amplitude of node i
Value, vjFor node j voltage magnitude, θiFor the voltage phase angle for node i, θjFor node j voltage phase angle, θijTo represent node i
With node j phase difference of voltage, GijFor the transconductance in bus admittance matrix between corresponding node i and j, BijFor node admittance
Mutual susceptance in matrix between corresponding node i and j, N is node total number;
The circuit branch road power measurement equation is expressed as:
Wherein, PijFor the active power of circuit branch road node side, QijFor the reactive power of circuit branch road node side, ycFor line
Road branch road susceptance over the ground;
The transformer branch power measurement equation is expressed as:
Wherein, Pij' be transformer branch node i side active power, Qij' be transformer branch node i side idle work(
Rate, Pji' be transformer branch node j sides active power, Qji' be transformer branch node j sides reactive power, K is transformation
The non-standard no-load voltage ratio of device, bTFor the susceptance of transformer standard side, node i side is standard side, and no-load voltage ratio is 1;Node j sides are non-standard
Side, no-load voltage ratio is K.
The network parameter and SCADA system parameter according to collection determines bus admittance matrix, measurement equation and zero
Injection equality constraint includes:
The zero injection equality constraint is expressed as:
C (x)=0
Wherein, c (x) is the measurement equation of zero injection node, represents that node injection active and reactive power is zero.
The amplitude and phase using data acquisition and the node voltage of supervisor control as state variable, and according to
Bus admittance matrix, measurement measurement equation and the zero injection equality constraint of the power system of structure calculate the refined of power system
Include than matrix and measurement residual error:
Jacobian matrix is comprising m measurement accounting equation to the local derviation of n state variable, and it is expressed as:
Wherein, H (x) represents Jacobian matrix, and it is m × n matrix;hm(x) m-th of measurement accounting equation, x are representedn
Represent n-th of state variable;
For node voltage amplitude, its Jacobi's element representation is:
Active power and reactive power are injected for node, its Jacobi's element representation is:
Wherein, GiiFor the self-conductance of node i, BiiIt is node i from susceptance;
For circuit branch road active power and reactive power and transformer branch active power and reactive power, its it is refined can
It is than element representation:
Represent that measurement is indexed with l, l=1,2 ..., m, then l-th of measurement residual error rlIt is expressed as:
rl=zl-hl(xk)
Wherein, zlRepresent l-th of measurement measuring value;xkRepresent k-th of state variable, k=1,2 ..., n;hl(xk) table
Show l-th of measurement accounting equation.
The piecewise nonlinear weight function is expressed as:
Wherein, w (rl) represent l-th of measurement piecewise nonlinear weight function, a be bad data detection threshold value;
Pass through w (rl) obtain the weight R of l-th of measurementl, it is expressed as:
Wherein,For the initial weight of l-th of measurement, andσlFor the standard deviation of l-th of measurement.
It is described that state variable is updated by iterative algorithm, and judge whether iterative algorithm reaches according to the correction of state variable
To convergence threshold, include if then obtaining the node voltage amplitude and phase of final power system:
The correction of state variable is obtained using Newton iteration method, is had:
Wherein, t represents iterations, Δ x(t)The correction of state variable during for the t times iteration, T represents transposition, x(t)For
State variable during the t times iteration, x(t+1)State variable during for the t+1 times iteration, h (x(t)) be the t times iteration when state
The corresponding measurement accounting equation of variable, H (x(t)) be the t time iteration when the corresponding Jacobian matrix of state variable, R for measurement
The weight matrix of amount, and R=[R1,R2,…,Rl,…,Rm]T;
According to obtained Δ x(t)Update state variable x;
Carry out convergence judgement, if meet max (| Δ x(t)|) < εx, then show convergence, terminate iterative calculation, output state
Assessment result;Otherwise iteration count plus one, re-start iterative calculation, wherein εxRepresent convergence threshold.
The present invention also provides a kind of power system robust condition estimating device based on piecewise nonlinear weight function, the dress
Put including:
Module is built, for joining according to the network of data acquisition and supervisor control the SCADA power system collected
Number and node and branch road operational factor, build the bus admittance matrix, measurement measurement equation and zero injection of power system
Equality constraint;
Computing module, for using the node voltage amplitude and phase of power system as state variable, and according to structure
Bus admittance matrix, measurement measurement equation and the zero injection equality constraint of power system calculate the Jacobean matrix of power system
Battle array and measurement residual error, and obtain by piecewise nonlinear weight function the weight of each measurement;
Judge module, for updating state variable by iterative algorithm, and judges iteration according to the correction of state variable
Whether algorithm reaches convergence threshold, if then obtaining the node voltage amplitude and phase of final power system.
It is described structure module specifically for:
The bus admittance matrix is expressed as:
Wherein, YiiRepresent the diagonal entry of bus admittance matrix, YijRepresent the non-diagonal line element of bus admittance matrix
Element, yijRepresent the branch impedance z between node i and node jijInverse.
It is described structure module specifically for:
The measurement measurement equation is expressed as:
Z=h (x)+r
Wherein, z represents measurement measuring value, and r represents measurement residual error, and x represents state variable, and h (x) calculates for measurement
Equation, it includes node power measurement equation, circuit branch road power measurement equation and transformer branch power measurement equation;
Wherein, node power measurement equation is expressed as:
Wherein, PiThe active power injected for node i, QiThe reactive power injected for node i, viFor the voltage amplitude of node i
Value, vjFor node j voltage magnitude, θiFor the voltage phase angle for node i, θjFor node j voltage phase angle, θijTo represent node i
With node j phase difference of voltage, GijFor the transconductance in bus admittance matrix between corresponding node i and j, BijFor node admittance
Mutual susceptance in matrix between corresponding node i and j, N is node total number;
The circuit branch road power measurement equation is expressed as:
Wherein, PijFor the active power of circuit branch road node side, QijFor the reactive power of circuit branch road node side, ycFor line
Road branch road susceptance over the ground;
The transformer branch power measurement equation is expressed as:
Wherein, Pij' be transformer branch node i side active power, Qij' be transformer branch node i side idle work(
Rate, Pji' be transformer branch node j sides active power, Qji' be transformer branch node j sides reactive power, K is transformation
The non-standard no-load voltage ratio of device, bTFor the susceptance of transformer standard side, node i side is standard side, and no-load voltage ratio is 1;Node j sides are non-standard
Side, no-load voltage ratio is K.
It is described structure module specifically for:
The zero injection equality constraint is expressed as:
C (x)=0
Wherein, c (x) is the measurement equation of zero injection node, represents that node injection active and reactive power is zero.
The computing module specifically for:
Jacobian matrix is comprising m measurement accounting equation to the local derviation of n state variable, and it is expressed as:
Wherein, H (x) represents Jacobian matrix, and it is m × n matrix;hm(x) m-th of measurement accounting equation, x are representedn
Represent n-th of state variable;
For node voltage amplitude, its Jacobi's element representation is:
Active power and reactive power are injected for node, its Jacobi's element representation is:
Wherein, GiiFor the self-conductance of node i, BiiIt is node i from susceptance;
For circuit branch road active power and reactive power and transformer branch active power and reactive power, its it is refined can
It is than element representation:
Represent that measurement is indexed with l, l=1,2 ..., m, then l-th of measurement residual error rlIt is expressed as:
rl=zl-hl(xk)
Wherein, zlRepresent l-th of measurement measuring value;xkRepresent k-th of state variable, k=1,2 ..., n;hl(xk) table
Show l-th of measurement accounting equation;
The piecewise nonlinear weight function is expressed as:
Wherein, w (rl) represent l-th of measurement piecewise nonlinear weight function, a be bad data detection threshold value;
Pass through w (rl) obtain the weight R of l-th of measurementl, it is expressed as:
Wherein,For the initial weight of l-th of measurement, andσlFor the standard deviation of l-th of measurement.
The judge module specifically for:
The correction of state variable is obtained using Newton iteration method, is had:
Wherein, t represents iterations, Δ x(t)The correction of state variable during for the t times iteration, T represents transposition, x(t)For
State variable during the t times iteration, x(t+1)State variable during for the t+1 times iteration, h (x(t)) be the t times iteration when shape
The corresponding measurement accounting equation of state variable, H (x(t)) be the t time iteration when the corresponding Jacobian matrix of state variable, R for measure
The weight matrix of measurement, and R=[R1,R2,…,Rl,…,Rm]T;
According to obtained Δ x(t)Update state variable x;
Carry out convergence judgement, if meet max (| Δ x(t)|) < εx, then show convergence, terminate iterative calculation, output state
Assessment result;Otherwise iteration count plus one, re-start iterative calculation, wherein εxRepresent convergence threshold.
Compared with immediate prior art, the technical scheme that the present invention is provided has the advantages that:
The present invention carries out bad data detection and identification and state estimation simultaneously, while state estimation is calculated, energy
Enough change measurement weight, bad data weight is constantly reduced in iteration, suppress the influence that bad data is brought, can be effective
Reduction amount of calculation, improve operating efficiency;
Compared with the weight function assignment of prior art, the present invention directly brings the residual error of measurement into segmentation nonlinear weight letter
Number, it is to avoid the calculating of standardized residual, can greatly promote calculating speed, save the calculating time, be very suitable for engineer applied;
The present invention uses piecewise nonlinear weight function, for measurement of the residual error within detection threshold value, maintains it to weigh
Weight, and the measurement for residual error beyond detection threshold value, reduce its weight, closer to zero by Exponential Type Weights function
And be not equal to zero, can effectively exclude bad data it is excessive when, correspondence neutral element is excessive and nonsingular matrix occur in matrix
Situation so that state estimation can effectively restrain, and greatly improve numerical stability when bad data and error in measurement coexist
Property.
Brief description of the drawings
Fig. 1 is the electric power system robust state estimation method flow based on piecewise nonlinear weight function in the embodiment of the present invention
Figure.
Embodiment
The present invention is described in further detail below in conjunction with the accompanying drawings.
The embodiment of the present invention provides a kind of electric power system robust state estimation method based on piecewise nonlinear weight function, such as
Shown in Fig. 1, this method detailed process is as follows:
S101:According to the network parameter of data acquisition and supervisor control the SCADA power system collected and section
(node includes node voltage amplitude, node injecting power and Branch Power Flow with branch road operational factor to be believed point with branch road operational factor
Breath), build the bus admittance matrix, measurement measurement equation and zero injection equality constraint of power system;
S102:The node voltage amplitude and phase of power system be as state variable using in S101, and according to the electricity of structure
Bus admittance matrix, measurement measurement equation and the zero injection equality constraint of Force system calculate the Jacobian matrix of power system
With measurement residual error, and the weight of each measurement is obtained by piecewise nonlinear weight function;
S103:State variable in S102 is updated by iterative algorithm, and judges that iteration is calculated according to the correction of state variable
Whether method reaches convergence threshold, if then obtaining the node voltage amplitude and phase of final power system.
In above-mentioned S101, according to the network parameter of data acquisition and supervisor control the SCADA power system collected
The node of sum and branch road operational factor, build bus admittance matrix, measurement measurement equation and zero injection of power system etc.
Formula constraint is specific as follows:
1) bus admittance matrix is expressed as:
Wherein, YiiRepresent the diagonal entry of bus admittance matrix, YijRepresent the non-diagonal line element of bus admittance matrix
Element, yijRepresent the branch impedance z between node i and node jijInverse.
2) according to the network parameter of collection and SCADA system parameter determine bus admittance matrix, measurement measurement equation with
And zero injection equality constraint include:
Measurement measurement equation is expressed as:
Z=h (x)+r
Wherein, z represents measurement measuring value, and r represents measurement residual error, and x represents state variable, and h (x) calculates for measurement
Equation, it includes node power measurement equation, circuit branch road power measurement equation and transformer branch power measurement equation;
Wherein, node power measurement equation is expressed as:
Wherein, PiThe active power injected for node i, QiThe reactive power injected for node i, viFor the voltage amplitude of node i
Value, vjFor node j voltage magnitude, θiFor the voltage phase angle for node i, θjFor node j voltage phase angle, θijTo represent node i
With node j phase difference of voltage, GijFor the transconductance in bus admittance matrix between corresponding node i and j, BijFor node admittance
Mutual susceptance in matrix between corresponding node i and j, N is node total number;
Circuit branch road power measurement equation is expressed as:
Wherein, PijFor the active power of circuit branch road node side, QijFor the reactive power of circuit branch road node side, ycFor line
Road branch road susceptance over the ground;
Transformer branch power measurement equation is expressed as:
Wherein, Pij' be transformer branch node i side active power, Qij' be transformer branch node i side idle work(
Rate, Pji' be transformer branch node j sides active power, Qji' be transformer branch node j sides reactive power, K is transformation
The non-standard no-load voltage ratio of device, bTFor the susceptance of transformer standard side, node i side is standard side, and no-load voltage ratio is 1;Node j sides are non-standard
Side, no-load voltage ratio is K.
3) zero injection equality constraint is expressed as:
C (x)=0
Wherein, c (x) is the measurement equation of zero injection node, represents that node injection active and reactive power is zero.
State variable includes node voltage amplitude and node voltage phase place, and measurement includes node voltage amplitude, node and noted
Enter active power, node injection reactive power, circuit branch road active power, circuit branch road reactive power, transformer branch active
Power and transformer branch reactive power.
In above-mentioned S102, the node voltage amplitude and phase of power system are as state variable using in S101, and according to structure
Bus admittance matrix, measurement measurement equation and the zero injection equality constraint for the power system built calculate the refined of power system can
Than matrix and measurement residual error, and the weight detailed process for obtaining each measurement by piecewise nonlinear weight function is as follows:
(1) Jacobian matrix is comprising m measurement accounting equation to the local derviation of n state variable, and it is expressed as:
Wherein, H (x) represents Jacobian matrix, and it is m × n matrix;hm(x) m-th of measurement accounting equation, x are representedn
Represent n-th of state variable;
1) for node voltage amplitude, its Jacobi's element representation is:
2) active power and reactive power are injected for node, its Jacobi's element representation is:
Wherein, GiiFor the self-conductance of node i, BiiIt is node i from susceptance;
3) for circuit branch road active power and reactive power and transformer branch active power and reactive power, its is refined
It is than element representation:
(2) represent that measurement is indexed with l, l=1,2 ..., m, then l-th of measurement residual error rlIt is expressed as:
rl=zl-hl(xk)
Wherein, zlRepresent l-th of measurement measuring value;xkRepresent k-th of state variable, k=1,2 ..., n;hl(xk) table
Show l-th of measurement accounting equation;
(3) piecewise nonlinear weight function is expressed as:
Wherein, w (rl) represent l-th of measurement piecewise nonlinear weight function, a be bad data detection threshold value;
(4) w (r are passed throughl) obtain the weight R of l-th of measurementl, it is expressed as:
Wherein,For the initial weight of l-th of measurement, andσlFor the standard deviation of l-th of measurement.
In above-mentioned S103, the state variable in S102 is updated by iterative algorithm, and sentence according to the correction of state variable
Whether disconnected iterative algorithm reaches convergence threshold, if the node voltage amplitude and phase that then obtain final power system are specifically wrapped
Include process as follows:
(1) correction of state variable is obtained using Newton iteration method, is had:
Wherein, t represents iterations, and the correction of state variable when Δ x (t) is the t times iteration, T represents transposition, x(t)
State variable during for the t times iteration, x(t+1)State variable during for the t+1 times iteration, h (x(t)) be the t times iteration when shape
The corresponding measurement accounting equation of state variable, H (x(t)) be the t time iteration when the corresponding Jacobian matrix of state variable, R for measure
The weight matrix of measurement, and R=[R1,R2,…,Rl,…,Rm]T;
(2) according to obtained Δ x(t)Update state variable x;
(3) carry out convergence judgement, if meet max (| Δ x(t)|) < εx, then show convergence, terminate iterative calculation, export shape
State assessment result;Otherwise iteration count plus one, re-start iterative calculation, wherein εxRepresent convergence threshold.
The embodiment of the present invention also provides a kind of power system robust state estimation based on piecewise nonlinear weight function simultaneously
Device, the device includes building module, computing module and judge module, and concrete function is as follows:
Module is built, for joining according to the network of data acquisition and supervisor control the SCADA power system collected
Number and node and branch road operational factor, build the bus admittance matrix, measurement measurement equation and zero injection of power system
Equality constraint;
Computing module, for using the node voltage amplitude and phase of power system as state variable, and according to structure
Bus admittance matrix, measurement measurement equation and the zero injection equality constraint of power system calculate the Jacobean matrix of power system
Battle array and measurement residual error, and obtain by piecewise nonlinear weight function the weight of each measurement;
Judge module, for updating state variable by iterative algorithm, and judges iteration according to the correction of state variable
Whether algorithm reaches convergence threshold, if then obtaining the node voltage amplitude and phase of final power system.
The bus admittance matrix that above-mentioned structure module is determined is expressed as:
Wherein, YiiRepresent the diagonal entry of bus admittance matrix, YijRepresent the non-diagonal line element of bus admittance matrix
Element, yijRepresent the branch impedance z between node i and node jijInverse.
The measurement measurement equation that above-mentioned structure module is determined is expressed as:
Z=h (x)+r
Wherein, z represents measurement measuring value, and r represents measurement residual error, and x represents state variable, and h (x) calculates for measurement
Equation, it includes node power measurement equation, circuit branch road power measurement equation and transformer branch power measurement equation;
Wherein, node power measurement equation is expressed as:
Wherein, PiThe active power injected for node i, QiThe reactive power injected for node i, viFor the voltage amplitude of node i
Value, vjFor node j voltage magnitude, θiFor the voltage phase angle for node i, θjFor node j voltage phase angle, θijTo represent node i
With node j phase difference of voltage, GijFor the transconductance in bus admittance matrix between corresponding node i and j, BijFor node admittance
Mutual susceptance in matrix between corresponding node i and j, N is node total number;
Wherein, circuit branch road power measurement equation is expressed as:
Wherein, PijFor the active power of circuit branch road node side, QijFor the reactive power of circuit branch road node side, ycFor line
Road branch road susceptance over the ground;
Wherein, transformer branch power measurement equation is expressed as:
Wherein, Pij' be transformer branch node i side active power, Qij' be transformer branch node i side idle work(
Rate, Pji' be transformer branch node j sides active power, Qji' be transformer branch node j sides reactive power, K is transformation
The non-standard no-load voltage ratio of device, bTFor the susceptance of transformer standard side, node i side is standard side, and no-load voltage ratio is 1;Node j sides are non-standard
Side, no-load voltage ratio is K.
The zero injection equality constraint that above-mentioned structure module is determined is expressed as:
C (x)=0
Wherein, c (x) is the measurement equation of zero injection node, represents that node injection active and reactive power is zero.
State variable includes node voltage amplitude and node voltage phase place, and the measurement includes node voltage amplitude, section
Point injection active power, node injection reactive power, circuit branch road active power, circuit branch road reactive power, transformer branch
Active power and transformer branch reactive power.
Above-mentioned computing module is using the node voltage amplitude and phase of power system as state variable, and according to structure
Bus admittance matrix, measurement measurement equation and the zero injection equality constraint of power system calculate the Jacobean matrix of power system
Battle array and measurement residual error, and the weight detailed process for obtaining each measurement by piecewise nonlinear weight function is as follows:
Jacobian matrix is comprising m measurement accounting equation to the local derviation of n state variable, and it is expressed as:
Wherein, H (x) represents Jacobian matrix, and it is m × n matrix;hm(x) m-th of measurement accounting equation, x are representedn
Represent n-th of state variable;
For node voltage amplitude, its Jacobi's element representation is:
Active power and reactive power are injected for node, its Jacobi's element representation is:
Wherein, GiiFor the self-conductance of node i, BiiIt is node i from susceptance;
For circuit branch road active power and reactive power and transformer branch active power and reactive power, its it is refined can
It is than element representation:
Represent that measurement is indexed with l, l=1,2 ..., m, then l-th of measurement residual error rlIt is expressed as:
rl=zl-hl(xk)
Wherein, zlRepresent l-th of measurement measuring value;xkRepresent k-th of state variable, k=1,2 ..., n;hl(xk) table
Show l-th of measurement accounting equation;
Piecewise nonlinear weight function is expressed as:
Wherein, w (rl) represent l-th of measurement piecewise nonlinear weight function, a be bad data detection threshold value;
Pass through w (rl) obtain the weight R of l-th of measurementl, it is expressed as:
Wherein,For the initial weight of l-th of measurement, andσlFor the standard deviation of l-th of measurement.
Above-mentioned judge module updates state variable by iterative algorithm, and judges iteration according to the correction of state variable
Whether algorithm reaches convergence threshold, if then obtaining the node voltage amplitude and phase detailed process of final power system such as
Under:
(1) correction of state variable is obtained using Newton iteration method, is had:
Wherein, t represents iterations, Δ x(t)The correction of state variable during for the t times iteration, T represents transposition, x(t)For
State variable during the t times iteration, x(t+1)State variable during for the t+1 times iteration, h (x(t)) be the t times iteration when state
The corresponding measurement accounting equation of variable, H (x(t)) be the t time iteration when the corresponding Jacobian matrix of state variable, R for measurement
The weight matrix of amount, and R=[R1,R2,…,Rl,…,Rm]T;
(2) according to obtained Δ x(t)Update state variable x;
(3) carry out convergence judgement, if meet max (| Δ x(t)|) < εx, then show convergence, terminate iterative calculation, export shape
State assessment result;Otherwise iteration count plus one, re-start iterative calculation, wherein εxRepresent convergence threshold.
The embodiment of the present invention is first when iterative calculation is proceeded by, and each measurement weights use the initial value of setting, profit
Initial value is tried to achieve with weighted least-squares method;Secondly weight is updated with piecewise nonlinear weight function amount of calculation according to measurement residuals, entered
Row measures the weights that measurement is changed before the iteration that weights change with residual error, each iteration, by calculating iterative equation more new state
Amount, and made comparisons with convergence criterion, until convergence, obtains state estimation result.
It should be understood by those skilled in the art that, embodiments herein can be provided as method, system or computer program
Product.Therefore, the application can be using the reality in terms of complete hardware embodiment, complete software embodiment or combination software and hardware
Apply the form of example.Moreover, the application can be used in one or more computers for wherein including computer usable program code
The computer program production that usable storage medium is implemented on (including but is not limited to magnetic disk storage, CD-ROM, optical memory etc.)
The form of product.
The application is the flow with reference to method, equipment (system) and computer program product according to the embodiment of the present application
Figure and/or block diagram are described.It should be understood that can be by every first-class in computer program instructions implementation process figure and/or block diagram
Journey and/or the flow in square frame and flow chart and/or block diagram and/or the combination of square frame.These computer programs can be provided
The processor of all-purpose computer, special-purpose computer, Embedded Processor or other programmable data processing devices is instructed to produce
A raw machine so that produced by the instruction of computer or the computing device of other programmable data processing devices for real
The device for the function of being specified in present one flow of flow chart or one square frame of multiple flows and/or block diagram or multiple square frames.
These computer program instructions, which may be alternatively stored in, can guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works so that the instruction being stored in the computer-readable memory, which is produced, to be included referring to
Make the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one square frame of block diagram or
The function of being specified in multiple square frames.
These computer program instructions can be also loaded into computer or other programmable data processing devices so that in meter
Series of operation steps is performed on calculation machine or other programmable devices to produce computer implemented processing, thus in computer or
The instruction performed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one
The step of function of being specified in individual square frame or multiple square frames.
Claims (14)
1. a kind of electric power system robust state estimation method, it is characterised in that methods described includes:
According to the network parameter of data acquisition and supervisor control the SCADA power system collected and node transported with branch road
Row parameter, builds the bus admittance matrix, measurement measurement equation and zero injection equality constraint of power system;
Node voltage amplitude and phase using power system is as state variable, and the node admittance of the power system according to structure
Matrix, measurement measurement equation and zero injection equality constraint calculate the Jacobian matrix and measurement residual error of power system, and
The weight of each measurement is obtained by piecewise nonlinear weight function;
State variable is updated by iterative algorithm, and judges whether iterative algorithm reaches convergence threshold according to the correction of state variable
Value, if then obtaining the node voltage amplitude and phase of final power system.
2. electric power system robust state estimation method according to claim 1, it is characterised in that the bus admittance matrix
It is expressed as:
Wherein, YiiRepresent the diagonal entry of bus admittance matrix, YijRepresent the off diagonal element of bus admittance matrix, yij
Represent the branch impedance z between node i and node jijInverse.
3. electric power system robust state estimation method according to claim 1, it is characterised in that the measurement amount measurement side
Journey is expressed as:
Z=h (x)+r
Wherein, z represents measurement measuring value, and r represents measurement residual error, and x represents state variable, and h (x) is measurement calculating side
Journey, it includes node power measurement equation, circuit branch road power measurement equation and transformer branch power measurement equation;
Wherein, node power measurement equation is expressed as:
Wherein, PiThe active power injected for node i, QiThe reactive power injected for node i, viFor the voltage magnitude of node i, vj
For node j voltage magnitude, θiFor the voltage phase angle for node i, θjFor node j voltage phase angle, θijTo represent node i and section
Point j phase difference of voltage, GijFor the transconductance in bus admittance matrix between corresponding node i and j, BijFor bus admittance matrix
Mutual susceptance between middle corresponding node i and j, N is node total number;
The circuit branch road power measurement equation is expressed as:
Wherein, PijFor the active power of circuit branch road node side, QijFor the reactive power of circuit branch road node side, ycFor circuit branch
Road susceptance over the ground;
The transformer branch power measurement equation is expressed as:
Wherein, Pij' be transformer branch node i side active power, Qij' be transformer branch node i side reactive power,
Pji' be transformer branch node j sides active power, Qji' be transformer branch node j sides reactive power, K is transformer
Non-standard no-load voltage ratio, bTFor the susceptance of transformer standard side, node i side is standard side, and no-load voltage ratio is 1;Node j sides are non-standard side,
No-load voltage ratio is K.
4. electric power system robust state estimation method according to claim 1, it is characterised in that the zero injection equation is about
Beam is expressed as:
C (x)=0
Wherein, c (x) is the measurement equation of zero injection node, represents that node injection active and reactive power is zero.
5. electric power system robust state estimation method according to claim 3, it is characterised in that it is described with data acquisition with
The amplitude and phase of the node voltage of supervisor control are as state variable, and the node admittance of the power system according to structure
Matrix, measurement measurement equation and zero injection equality constraint calculate the Jacobian matrix and measurement residual error bag of power system
Include:
Jacobian matrix is comprising m measurement accounting equation to the local derviation of n state variable, and it is expressed as:
Wherein, H (x) represents Jacobian matrix, and it is m × n matrix;hm(x) m-th of measurement accounting equation, x are representednRepresent
N-th of state variable;
For node voltage amplitude, its Jacobi's element representation is:
Active power and reactive power are injected for node, its Jacobi's element representation is:
Wherein, GiiFor the self-conductance of node i, BiiIt is node i from susceptance;
For circuit branch road active power and reactive power and transformer branch active power and reactive power, its Jacobi member
Element is expressed as:
Represent that measurement is indexed with l, l=1,2 ..., m, then l-th of measurement residual error rlIt is expressed as:
rl=zl-hl(xk)
Wherein, zlRepresent l-th of measurement measuring value;xkRepresent k-th of state variable, k=1,2 ..., n;hl(xk) represent l
Individual measurement accounting equation.
6. electric power system robust state estimation method according to claim 5, it is characterised in that described by being segmented non-thread
The weight that property weight function obtains each measurement includes:
The piecewise nonlinear weight function is expressed as:
Wherein, w (rl) represent l-th of measurement piecewise nonlinear weight function, a be bad data detection threshold value;
Pass through w (rl) obtain the weight R of l-th of measurementl, it is expressed as:
Wherein,For the initial weight of l-th of measurement, andσlFor the standard deviation of l-th of measurement.
7. electric power system robust state estimation method according to claim 6, it is characterised in that described to pass through iterative algorithm
State variable is updated, and judges whether iterative algorithm reaches convergence threshold according to the correction of state variable, if then obtaining most
The node voltage amplitude and phase of whole power system include:
The correction of state variable is obtained using Newton iteration method, is had:
Wherein, t represents iterations, Δ x(t)The correction of state variable during for the t times iteration, T represents transposition, x(t)For t
State variable during secondary iteration, x(t+1)State variable during for the t+1 times iteration, h (x(t)) be the t time iteration when state change
Measure corresponding measurement accounting equation, H (x(t)) be the t times iteration when the corresponding Jacobian matrix of state variable, R is measurement
Weight matrix, and R=[R1,R2,…,Rl,…,Rm]T;
According to obtained Δ x(t)Update state variable x;
Carry out convergence judgement, if meet max (| Δ x(t)|) < εx, then show convergence, terminate iterative calculation, output state assesses knot
Really;Otherwise iteration count plus one, re-start iterative calculation, wherein εxRepresent convergence threshold.
8. a kind of power system robust condition estimating device, it is characterised in that described device includes:
Build module, for the network parameter according to data acquisition and supervisor control the SCADA power system collected and
Node and branch road operational factor, build the bus admittance matrix, measurement measurement equation and zero injection equation of power system
Constraint;
Computing module, for using the node voltage amplitude and phase of power system as state variable, and according to the electric power of structure
The bus admittance matrix of system, measurement measurement equation and zero injection equality constraint calculate power system Jacobian matrix and
Measurement residual error, and obtain by piecewise nonlinear weight function the weight of each measurement;
Judge module, for updating state variable by iterative algorithm, and judges iterative algorithm according to the correction of state variable
Whether convergence threshold is reached, if then obtaining the node voltage amplitude and phase of final power system.
9. power system robust condition estimating device according to claim 8, it is characterised in that the structure module is specific
For:
The bus admittance matrix is expressed as:
Wherein, YiiRepresent the diagonal entry of bus admittance matrix, YijRepresent the off diagonal element of bus admittance matrix, yij
Represent the branch impedance z between node i and node jijInverse.
10. power system robust condition estimating device according to claim 8, it is characterised in that the structure module tool
Body is used for:
The measurement measurement equation is expressed as:
Z=h (x)+r
Wherein, z represents measurement measuring value, and r represents measurement residual error, and x represents state variable, and h (x) is measurement calculating side
Journey, it includes node power measurement equation, circuit branch road power measurement equation and transformer branch power measurement equation;
Wherein, node power measurement equation is expressed as:
Wherein, PiThe active power injected for node i, QiThe reactive power injected for node i, viFor the voltage magnitude of node i, vj
For node j voltage magnitude, θiFor the voltage phase angle for node i, θjFor node j voltage phase angle, θijTo represent node i and section
Point j phase difference of voltage, GijFor the transconductance in bus admittance matrix between corresponding node i and j, BijFor bus admittance matrix
Mutual susceptance between middle corresponding node i and j, N is node total number;
The circuit branch road power measurement equation is expressed as:
Wherein, PijFor the active power of circuit branch road node side, QijFor the reactive power of circuit branch road node side, ycFor circuit branch
Road susceptance over the ground;
The transformer branch power measurement equation is expressed as:
Wherein, Pij' be transformer branch node i side active power, Qij' be transformer branch node i side reactive power,
Pji' be transformer branch node j sides active power, Qji' be transformer branch node j sides reactive power, K is transformer
Non-standard no-load voltage ratio, bTFor the susceptance of transformer standard side, node i side is standard side, and no-load voltage ratio is 1;Node j sides are non-standard side,
No-load voltage ratio is K.
11. power system robust condition estimating device according to claim 8, it is characterised in that the structure module tool
Body is used for:
The zero injection equality constraint is expressed as:
C (x)=0
Wherein, c (x) is the measurement equation of zero injection node, represents that node injection active and reactive power is zero.
12. power system robust condition estimating device according to claim 10, it is characterised in that the computing module tool
Body is used for:
Jacobian matrix is comprising m measurement accounting equation to the local derviation of n state variable, and it is expressed as:
Wherein, H (x) represents Jacobian matrix, and it is m × n matrix;hm(x) m-th of measurement accounting equation, x are representednRepresent
N-th of state variable;
For node voltage amplitude, its Jacobi's element representation is:
Active power and reactive power are injected for node, its Jacobi's element representation is:
Wherein, GiiFor the self-conductance of node i, BiiIt is node i from susceptance;
For circuit branch road active power and reactive power and transformer branch active power and reactive power, its Jacobi member
Element is expressed as:
Represent that measurement is indexed with l, l=1,2 ..., m, then l-th of measurement residual error rlIt is expressed as:
rl=zl-hl(xk)
Wherein, zlRepresent l-th of measurement measuring value;xkRepresent k-th of state variable, k=1,2 ..., n;hl(xk) represent l
Individual measurement accounting equation.
13. power system robust condition estimating device according to claim 12, it is characterised in that the piecewise nonlinear
Weight function is expressed as:
Wherein, w (rl) represent l-th of measurement piecewise nonlinear weight function, a be bad data detection threshold value;
Pass through w (rl) obtain the weight R of l-th of measurementl, it is expressed as:
Wherein,For the initial weight of l-th of measurement, andσlFor the standard deviation of l-th of measurement.
14. power system robust condition estimating device according to claim 13, it is characterised in that the judge module tool
Body is used for:
The correction of state variable is obtained using Newton iteration method, is had:
Wherein, t represents iterations, Δ x(t)The correction of state variable during for the t times iteration, T represents transposition, x(t)For t
State variable during secondary iteration, x(t+1)State variable during for the t+1 times iteration, h (x(t)) be the t time iteration when state change
Measure corresponding measurement accounting equation, H (x(t)) be the t times iteration when the corresponding Jacobian matrix of state variable, R is measurement
Weight matrix, and R=[R1,R2,…,Rl,…,Rm]T;
According to obtained Δ x(t)Update state variable x;
Carry out convergence judgement, if meet max (| Δ x(t)|) < εx, then show convergence, terminate iterative calculation, output state assesses knot
Really;Otherwise iteration count plus one, re-start iterative calculation, wherein εxRepresent convergence threshold.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710137890.4A CN107016489A (en) | 2017-03-09 | 2017-03-09 | A kind of electric power system robust state estimation method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710137890.4A CN107016489A (en) | 2017-03-09 | 2017-03-09 | A kind of electric power system robust state estimation method and device |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107016489A true CN107016489A (en) | 2017-08-04 |
Family
ID=59439714
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710137890.4A Pending CN107016489A (en) | 2017-03-09 | 2017-03-09 | A kind of electric power system robust state estimation method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107016489A (en) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107994567A (en) * | 2017-11-22 | 2018-05-04 | 华北电力大学 | A kind of broad sense Fast decoupled state estimation method |
CN108281967A (en) * | 2018-01-04 | 2018-07-13 | 国网浙江省电力有限公司电力科学研究院 | A kind of algorithm startup method improving electric power system tide computational convergence |
CN108879695A (en) * | 2018-07-13 | 2018-11-23 | 海南电网有限责任公司电力科学研究院 | A kind of foundation and calculation method of DC flow model |
CN109888773A (en) * | 2019-02-25 | 2019-06-14 | 武汉大学 | A kind of electric system multizone distributions appraisal procedure |
CN110232199A (en) * | 2018-11-30 | 2019-09-13 | 南京南瑞继保电气有限公司 | A kind of bilinearity WLAV robust state estimation method considering that matrix size is flexible |
CN110417009A (en) * | 2019-07-29 | 2019-11-05 | 天津大学 | Power distribution network based on Different sampling period data mixes robust state estimation method |
CN110619105A (en) * | 2019-10-14 | 2019-12-27 | 国网山东省电力公司潍坊供电公司 | Power transmission line temperature estimation method based on quantity measurement and heat balance equation |
CN110781450A (en) * | 2019-09-17 | 2020-02-11 | 广西电网有限责任公司电力科学研究院 | Bad data detection system and method for distribution feeder voltage measurement data |
CN111711186A (en) * | 2020-05-28 | 2020-09-25 | 中国南方电网有限责任公司 | Power system PQ decomposition state estimation method based on GPU parallel computing |
CN112383048A (en) * | 2020-10-20 | 2021-02-19 | 重庆大学 | Medium-voltage feeder line power grid three-phase state estimation method considering measurement characteristics of public and special distribution transformer |
CN112398117A (en) * | 2020-09-24 | 2021-02-23 | 北京航空航天大学 | False data injection attack construction and defense method causing line load overload |
CN112636325A (en) * | 2020-11-11 | 2021-04-09 | 广东电网有限责任公司广州供电局 | Power distribution network robust state calculation method and device, terminal and storage medium |
CN112865096A (en) * | 2021-03-18 | 2021-05-28 | 国网上海市电力公司 | Power distribution network state estimation method and system considering PMU (phasor measurement Unit) measurement phase angle deviation |
CN112993989A (en) * | 2021-03-05 | 2021-06-18 | 广东电网有限责任公司广州供电局 | Robust state estimation data processing method for active power distribution system |
CN113032718A (en) * | 2021-03-29 | 2021-06-25 | 北京华大九天科技股份有限公司 | Method for solving Newton iterative algorithm dead loop in circuit simulation |
CN113139288A (en) * | 2021-04-21 | 2021-07-20 | 广东电网有限责任公司广州供电局 | Power distribution network distributed robust state estimation method and device based on GPU (graphics processing Unit) secondary acceleration |
CN113255959A (en) * | 2021-04-09 | 2021-08-13 | 广东电网有限责任公司电力调度控制中心 | Power system dynamic state estimation method and system |
CN113313339A (en) * | 2021-03-26 | 2021-08-27 | 贵州大学 | Power distribution network dynamic state estimation algorithm |
CN115470736A (en) * | 2022-09-29 | 2022-12-13 | 长沙学院 | Power system dynamic behavior modeling method adaptive to variable working condition operation of energy storage power station |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101219545B1 (en) * | 2011-09-14 | 2013-01-09 | 주식회사 파워이십일 | Optimized parameter estimation method for power system |
CN104765962A (en) * | 2015-04-08 | 2015-07-08 | 河海大学 | Temperature variation considering type state estimation method of power system |
CN105514977A (en) * | 2015-11-27 | 2016-04-20 | 华北电力大学 | Hyperbolic cosine type robust state estimation method for power system state |
CN105512502A (en) * | 2016-01-13 | 2016-04-20 | 重庆大学 | Weight function least square state estimation method based on residual normalization |
WO2016078477A1 (en) * | 2014-11-18 | 2016-05-26 | 国电南瑞科技股份有限公司 | Transformer substation three-phase linear generalized state estimation method |
-
2017
- 2017-03-09 CN CN201710137890.4A patent/CN107016489A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101219545B1 (en) * | 2011-09-14 | 2013-01-09 | 주식회사 파워이십일 | Optimized parameter estimation method for power system |
WO2016078477A1 (en) * | 2014-11-18 | 2016-05-26 | 国电南瑞科技股份有限公司 | Transformer substation three-phase linear generalized state estimation method |
CN104765962A (en) * | 2015-04-08 | 2015-07-08 | 河海大学 | Temperature variation considering type state estimation method of power system |
CN105514977A (en) * | 2015-11-27 | 2016-04-20 | 华北电力大学 | Hyperbolic cosine type robust state estimation method for power system state |
CN105512502A (en) * | 2016-01-13 | 2016-04-20 | 重庆大学 | Weight function least square state estimation method based on residual normalization |
Non-Patent Citations (1)
Title |
---|
闫丽梅等: "状态估计与不良数据检测方法及其在大庆油田电网中的应用", 《化工自动化及仪表》 * |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107994567A (en) * | 2017-11-22 | 2018-05-04 | 华北电力大学 | A kind of broad sense Fast decoupled state estimation method |
CN107994567B (en) * | 2017-11-22 | 2020-09-01 | 华北电力大学 | Generalized fast decomposition state estimation method |
CN108281967A (en) * | 2018-01-04 | 2018-07-13 | 国网浙江省电力有限公司电力科学研究院 | A kind of algorithm startup method improving electric power system tide computational convergence |
CN108879695A (en) * | 2018-07-13 | 2018-11-23 | 海南电网有限责任公司电力科学研究院 | A kind of foundation and calculation method of DC flow model |
CN110232199A (en) * | 2018-11-30 | 2019-09-13 | 南京南瑞继保电气有限公司 | A kind of bilinearity WLAV robust state estimation method considering that matrix size is flexible |
CN110232199B (en) * | 2018-11-30 | 2023-12-19 | 南京南瑞继保电气有限公司 | Bilinear WLAV robust state estimation method considering matrix scale expansion |
CN109888773B (en) * | 2019-02-25 | 2022-08-05 | 武汉大学 | Multi-region distributed state evaluation method for power system |
CN109888773A (en) * | 2019-02-25 | 2019-06-14 | 武汉大学 | A kind of electric system multizone distributions appraisal procedure |
CN110417009A (en) * | 2019-07-29 | 2019-11-05 | 天津大学 | Power distribution network based on Different sampling period data mixes robust state estimation method |
CN110781450A (en) * | 2019-09-17 | 2020-02-11 | 广西电网有限责任公司电力科学研究院 | Bad data detection system and method for distribution feeder voltage measurement data |
CN110781450B (en) * | 2019-09-17 | 2022-12-16 | 广西电网有限责任公司电力科学研究院 | Bad data detection method and system for distribution feeder voltage measurement data |
CN110619105A (en) * | 2019-10-14 | 2019-12-27 | 国网山东省电力公司潍坊供电公司 | Power transmission line temperature estimation method based on quantity measurement and heat balance equation |
CN111711186B (en) * | 2020-05-28 | 2023-05-02 | 中国南方电网有限责任公司 | Power system PQ decomposition state estimation method based on GPU parallel computing |
CN111711186A (en) * | 2020-05-28 | 2020-09-25 | 中国南方电网有限责任公司 | Power system PQ decomposition state estimation method based on GPU parallel computing |
CN112398117A (en) * | 2020-09-24 | 2021-02-23 | 北京航空航天大学 | False data injection attack construction and defense method causing line load overload |
CN112383048A (en) * | 2020-10-20 | 2021-02-19 | 重庆大学 | Medium-voltage feeder line power grid three-phase state estimation method considering measurement characteristics of public and special distribution transformer |
CN112636325A (en) * | 2020-11-11 | 2021-04-09 | 广东电网有限责任公司广州供电局 | Power distribution network robust state calculation method and device, terminal and storage medium |
CN112993989B (en) * | 2021-03-05 | 2022-12-16 | 广东电网有限责任公司广州供电局 | Robust state estimation data processing method for active power distribution system |
CN112993989A (en) * | 2021-03-05 | 2021-06-18 | 广东电网有限责任公司广州供电局 | Robust state estimation data processing method for active power distribution system |
CN112865096A (en) * | 2021-03-18 | 2021-05-28 | 国网上海市电力公司 | Power distribution network state estimation method and system considering PMU (phasor measurement Unit) measurement phase angle deviation |
CN113313339A (en) * | 2021-03-26 | 2021-08-27 | 贵州大学 | Power distribution network dynamic state estimation algorithm |
CN113032718B (en) * | 2021-03-29 | 2022-05-24 | 北京华大九天科技股份有限公司 | Method and device for solving Newton iteration algorithm dead loop in circuit simulation |
CN113032718A (en) * | 2021-03-29 | 2021-06-25 | 北京华大九天科技股份有限公司 | Method for solving Newton iterative algorithm dead loop in circuit simulation |
CN113255959A (en) * | 2021-04-09 | 2021-08-13 | 广东电网有限责任公司电力调度控制中心 | Power system dynamic state estimation method and system |
CN113139288A (en) * | 2021-04-21 | 2021-07-20 | 广东电网有限责任公司广州供电局 | Power distribution network distributed robust state estimation method and device based on GPU (graphics processing Unit) secondary acceleration |
CN115470736A (en) * | 2022-09-29 | 2022-12-13 | 长沙学院 | Power system dynamic behavior modeling method adaptive to variable working condition operation of energy storage power station |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107016489A (en) | A kind of electric power system robust state estimation method and device | |
CN104134999B (en) | Distribution network based on multi-data source measures the practical method of calculation of efficiency analysis | |
CN104836223B (en) | Power grid parameter error and bad data coordinated identification and estimation method | |
CN103944165B (en) | A kind of bulk power grid parameter identification method of estimation | |
CN103532137B (en) | A kind of method for estimating state of three-phase and four-line low-voltage network | |
CN104184144B (en) | A kind of robust state estimation method for multi-voltage grade electric network model | |
CN107577870A (en) | The distribution network voltage power sensitivity robust estimation method measured based on synchronized phasor | |
CN103972884A (en) | Electric system state estimation method | |
CN103886193B (en) | Fuzzy self-adaptation robust estimation method of electric power system | |
CN110417011A (en) | A kind of online dynamic secure estimation method based on mutual information Yu iteration random forest | |
CN102403720A (en) | Hyper-real-time setting method for superposed time sequences based on transient voltage safety margin | |
CN106228459A (en) | Equivalent reliability estimation method based on Monte Carlo | |
CN105512502A (en) | Weight function least square state estimation method based on residual normalization | |
CN103473602A (en) | Theoretical line loss data prediction system and prediction method of power grid | |
CN108565852A (en) | A kind of Contingency screening and ranking method of the bulk power grid Voltage Stability Evaluation of three progress step by step | |
CN103324858A (en) | Three-phase load flow state estimation method of power distribution network | |
CN106356840A (en) | Method and system for estimating states of regional electric power systems on basis of synchronous phasor measurement | |
CN104182644A (en) | Power distribution network state estimation method fusing distributed generation characteristics | |
CN109829246B (en) | Line parameter identification method based on parameter comprehensive suspicion | |
CN104866714A (en) | Self-adaptive nuclear density robust state estimation method for power system | |
CN104617575A (en) | Method for using static security analysis to assess flow insolubility of anticipated faults | |
CN101788608A (en) | Method for evaluating reactance parameters of independent three-winding transformer | |
CN111756031B (en) | Power grid operation trend estimation method and system | |
CN105977963B (en) | Main and distribution network integrated distributed state estimation method based on two-stage method | |
CN103942416B (en) | Voltage estimation method based on weighted node spanning tree |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
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: 20170804 |