Summary of the invention
The objective of the invention is to overcome existing power system voltage stabilization weak node and confirm the deficiency of method, provide a kind of have calculate simple, result accurately, be easy to realize and need not the repetition trend calculate and confirm that voltage stablizes the method for weak node.
Realize that the technical scheme that the present invention adopted is, a kind ofly confirm that based on the matrix perturbation theory voltage stablizes the weak node method, it is characterized in that it may further comprise the steps:
1) based on the structure of the Jacobian matrix conditional number of matrix perturbation theory
The equation that is used for electric power system tide calculating is:
JV=W (1)
Wherein: J is used for the Jacobian matrix that system load flow calculates,
V representes the variation column vector [Δ θ Δ U/U] of node phase angle and voltage
T,
W representes the variation column vector [Δ P Δ Q] of node active power and reactive power
T
When system receives extraneous disturbance; System can meritoriously on the basis of present operating point redistribute with reactive power; Make system load flow reach a new operating point; Corresponding variation has also taken place in the element that this moment, total system was used for the Jacobian matrix of trend iteration; Power flow equation reaches balance at new operating point place, and there is following relation in system at new operating point place:
(J+ΔJ)(V+ΔV)=(W+ΔW) (2)
Wherein: Δ J is that system is measured by the change of Jacobian matrix after the disturbance,
Δ V is that system is subjected to disturbance posterior nodal point voltage to change column vector,
Δ W is that system is subjected to disturbance posterior nodal point power to change column vector,
Make J '=J+ Δ J, then J ' is the Jacobian matrix of system at new operating point place power flow equation;
Establishing a capital corresponding to the Jacobian matrix J at two operating point places of trend convergence and J ' is non-singular matrix, and when system configuration did not change, J was identical with J ' dimension, can get according to the disturbance enclosure theorem of matrix inversion and equation group in the matrix perturbation theory:
With
Wherein:
κ=||J||
2*||J
-1||
2 (5)
Then, (5) formula can be defined the conditional number of Jacobian matrix according to the matrix perturbation theory;
Wherein: J is a Jacobian matrix in the power flow equation,
J
-1Be the inverse matrix of J,
|| ||
2The 2-norm of representing matrix or vector;
2) stablize the structure of weak node index based on the voltage of Jacobian matrix conditional number sensitivity
State variable and control variables are respectively X=(x in the note system
1, x
2..., x
n)
T, and Y=(y
1, y
2... y
n)
T, n is the node number in the system,
By the differentiate rule of compound function, conditional number to the sensitivity of state variable does in the formula (5)
Wherein: δ
MaxAnd δ
MinBe the maximum and the minimum singular value of Jacobian matrix,
According to the differentiate rule of implicit function, singular value to the partial derivative of state variable in the system is:
That is:
To sum up, can get the sensitivity degree of conditional number to control variables:
Wherein: V
MaxAnd V
MinBe respectively δ
MaxAnd δ
MinCorresponding left singular vector,
U
MaxAnd U
MinBe respectively δ
MaxAnd δ
MinCorresponding right singular vector;
3) foundation of node load nargin index
The size of Jacobian matrix conditional number sensitivity has reflected the sensitivity that the stability of a system changes node power, and the definite condition number is a node load nargin index to the inverse of node control variable sensitivity, that is:
The present invention is based on the matrix perturbation theory and confirm that it is to utilize the definite voltage of Jacobian matrix conditional number sensitivity realization to stablize the method for weak node that voltage is stablized the weak node method; Confirm that with traditional weak node method compares the advantage applies that has and exist in the large-scale electrical power system static electric voltage stability analysis in modern times: 1. Jacobian matrix conditional number and sensitivity mathematical derivation process thereof are rigorous, the explicit physical meaning of model; 2. under current running status, can confirm the weak node in the system, calculate node load nargin index simultaneously; 3. Jacobian matrix conditional number and Sensitivity calculation thereof are simple, and be consuming time few; The present invention is applied to the matrix perturbation theory in the power system voltage stabilization analysis; Opened up new approach for the research and analysis of power system voltage stabilization property, improved traditional static voltage stable in weak node confirm to repeat the defective that trend is calculated in the continuous tide method with calculated load nargin; It calculates simply, speed is fast, and practical applications is worth high.
Embodiment
Of the present inventionly a kind ofly confirm that based on the matrix perturbation theory voltage stablizes the weak node method, be used for the static voltage stability analysis of electric power system, may further comprise the steps:
1) based on the structure of the Jacobian matrix conditional number of matrix perturbation theory
The equation that is used for electric power system tide calculating is:
JV=W (1)
Wherein: J is used for the Jacobian matrix that system load flow calculates,
V representes the variation column vector [Δ θ Δ U/U] of node phase angle and voltage
T,
W representes the variation column vector [Δ P Δ Q] of node active power and reactive power
T
When system received an extraneous disturbance (like the switching of load, the starting of motor etc.), system can redistribute meritorious and idle distribution on the basis of present operating point, make system load flow reach a new operating point.Corresponding variation has also taken place in the element that this moment, total system was used for the Jacobian matrix of trend iteration, and power flow equation reaches balance at new operating point place.There is following relation in system at new operating point place:
(J+ΔJ)(V+ΔV)=(W+ΔW) (2)
Wherein: Δ J is that system is measured by the change of Jacobian matrix after the disturbance;
Δ V is that system is subjected to disturbance posterior nodal point voltage to change column vector;
Δ W is that system is subjected to disturbance posterior nodal point power to change column vector;
Make J '=J+ Δ J, then J ' is the Jacobian matrix of system at new operating point place power flow equation.
Establishing a capital corresponding to the Jacobian matrix J at two operating point places of trend convergence and J ' is non-singular matrix, and when system configuration did not change, J was identical with J ' dimension, can get according to the disturbance enclosure theorem of matrix inversion and equation group in the matrix perturbation theory:
With
Wherein:
κ=||J||
2*||J
-1||
2 (5)
Then, (5) formula can be defined the conditional number of Jacobian matrix according to the matrix perturbation theory;
Wherein: J is a Jacobian matrix in the power flow equation,
J
-1Be the inverse matrix of J,
|| ||
2The 2-norm of representing matrix or vector;
2) stablize the structure of weak node index based on the voltage of Jacobian matrix conditional number sensitivity
State variable and control variables are respectively X=(x in the note system
1, x
2..., x
n)
T, and Y=(y
1, y
2... y
n)
T, n is the node number in the system,
By the differentiate rule of compound function, conditional number to the sensitivity of state variable does in the formula (5)
Wherein: δ
MaxAnd δ
MinBe the maximum and the minimum singular value of Jacobian matrix,
According to the differentiate rule of implicit function, singular value to the partial derivative of state variable in the system is:
That is:
To sum up, can get the sensitivity degree of conditional number to control variables:
Wherein: V
MaxAnd V
MinBe respectively δ
MaxAnd δ
MinCorresponding left singular vector,
U
MaxAnd U
MinBe respectively δ
MaxAnd δ
MinCorresponding right singular vector;
3) foundation of node load nargin index
The size of Jacobian matrix conditional number sensitivity has reflected the sensitivity that the stability of a system changes node power, and the definite condition number is a node load nargin index to the inverse of node control variable sensitivity, that is:
Can know by (10) formula, weak node in the system, the value of conditional number sensitivity is big more, and its load margin index is more little.When system was in voltage and stablizes critical condition, it is infinitely great that the conditional number sensitivity of each node is tending towards, and the nargin index will approach 0.
Instantiation:
It is that example is calculated and analyzed the validity of Jacobian matrix and conditional number thereof that the present invention comes verification system with the IEEE57 node system.IEEE57 node system structure chart as shown in Figure 1.
Under different trend modes, the conditional number of system is calculated the variation tendency of conditional number when being illustrated in figure 2 as the total system load, the number of times of k for increasing among the figure respectively with ratio 20% increase.
Can find out that in Fig. 2 the conditional number that different load level is corresponding is also inequality, and along with the increase of load level, conditional number increases gradually, when system was put near voltage collapse, the conditional number increasing degree was bigger.Aforementioned calculation result shows that the size of Jacobian matrix conditional number has reflected the voltage maintenance level of system to a certain extent, and the current operating point of its big more illustrative system of value is more near the voltage collapse point.
In the IEEE57 node system, as control variables, its node voltage is as state variable, respectively in the sensitivity of two kinds of load level computing system conditional numbers to the load bus control variables with the reactive power of each load bus.Its result of calculation as shown in Figure 3.The bigger part of nodes of conditional number Sensitirity va1ue has only drawn among Fig. 3.
State 1: the initial load state of system;
State 2: under the initial load state, all load buses increase power 10% with permanent power factor (PF).
In Fig. 3, the difference in size of sensitivity is comparatively remarkable, and it is more obvious that weak node and weakness zone seem.Sensitivity under the state 2 is greater than the sensitivity under the state 1 on the whole, and this is because the load level relative status 1 of state 2 wants high, compares state 1 stability of a system and worsens, and more approaches the critical condition of voltage collapse.But magnitude relationship does not change between each the node sensitivity under the two states.Explanation is under above-mentioned two kinds of trend modes, and the weak node of system and weakness zone do not change.As space is limited, only listed bigger part of nodes data and the ordering of system sensitivity in the table 1.
Weak node ordering under the table 1IEEE57 node system different conditions
Utilize the continuous tide method to calculate the meritorious nargin of each weak node, its result of calculation is as shown in table 2.
The meritorious nargin of each weak node of table 2IEEE57 node system
Table 1 is weak node ranking results and the voltage margin index result of calculation thereof under two kinds of load methods of IEEE57 node system.Can find out that by table 1 data in system, its nargin index of node a little less than the relative thin shows that more near 0 its load margin is more little.Table 2 is two kinds of continuous tide result of calculations under the load method.Contrast table 1 and table 2 result, the nargin index under node 33 two states is respectively 0.0958 and 0.0592, and the continuous tide result of calculation of nargin is 19.38MW and 10.45MW under its two states simultaneously.The checkout result of other nodes has These characteristics equally, and characteristics show the degree of closeness that can cross good reflection system and critical condition based on the nargin index of Jacobian matrix conditional number sensitivity.
Be the characteristic of further checking system weak node, utilize the continuous tide method to calculate the meritorious nargin of each weak node, its result of calculation is as shown in table 3.
The meritorious nargin of each weak node of table 3IEEE57 node system
With the power of each load bus in the growth pattern increase system of constant power factor, change curve such as Fig. 4 of each the weak node load margin index in the table 3.
Can find that from above example the size of Jacobian matrix conditional number can reflect the distance of system apart from voltage collapse point preferably.Through the calculating checking to static load nargin of continuous tide method, declaration condition is counted sensitivity and can be accurately the weak node of system be sorted, and its nargin index has been reacted in the system each node current state accurately apart from the distance of voltage collapse critical point.