Method for obtaining symmetric power transfer coefficients under simultaneous change of sources and loads in AC power networks

09774188 · 2017-09-26

Assignee

Inventors

Cpc classification

International classification

Abstract

The present invention belongs to electric power engineering field, relates to a method for obtaining a set of symmetric power transfer coefficients under simultaneous change of sources and loads in AC power networks, which comprises the steps of: firstly establishing a linear function of a lossy branch transferred power in terms of buses voltage angles according to a nonlinear function of a branch transferred power, parameters and operation features of the AC power network; then establishing symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles by combining the linear function of the lossy branch transferred power in terms of buses voltage angles and buses injection powers of power sources and loads, in turn establishing symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads; and finally obtaining the symmetric power transfer coefficients from buses injection powers of power sources and loads to the lossy branch transferred power under simultaneous change of sources and loads by using the linear functions just mentioned above. The obtained set of symmetric power transfer coefficients is unique, follows electric circuit laws, is applicable for the practical situations that power sources and loads change at the same time and transmission losses need considering in the AC power network, and truly reflects the substantive characteristics of the power transfer from bus injection powers to branch transferred powers.

Claims

1. A method for obtaining a set of symmetric power transfer coefficients under simultaneous change of sources and loads in alternating current (AC) power networks, which comprises: establishing a linear function of a lossy branch transferred power in terms of buses voltage angles according to given AC power network parameters, a nonlinear function of a branch transferred power in the AC power network, and operation features of the AC power network; establishing symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles according to the linear function of the lossy branch transferred power in terms of buses voltage angles and given buses injection powers of power sources and loads; establishing symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads according to the symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles; obtaining the symmetric power transfer coefficients from buses injection powers of power sources and loads to the lossy branch transferred power under simultaneous change of sources and loads according to the symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads and the linear function of the lossy branch transferred power in terms of buses voltage angles.

2. The method for obtaining a set of symmetric power transfer coefficients under simultaneous change of sources and loads in AC power networks according to claim 1, wherein the step of establishing a linear function of the lossy branch transferred power in terms of buses voltage angles according to the given AC power network parameters, the nonlinear function of the branch transferred power in the AC power network, and operation features of the AC power network comprises the step of: according to the given AC power network parameters, the nonlinear function of the branch transferred power in the AC power network, and the operation features of the AC power network, establishing a linear function of the lossy branch transferred power in terms of buses voltage angles by the following equation:
P.sub.ij=−b.sub.ij(θ.sub.i−θ.sub.j) where i and j are the numbers of two arbitrary buses in the AC power network respectively, i and j are natural numbers and equal to 1, 2 . . . , n, n is the total number of buses in the AC power network, n is natural number and is the given AC power network parameter; ij is the branch between bus i and bus j; P.sub.ij is the lossy branch transferred power entering branch ij from bus i; θ.sub.i and θ.sub.j are voltage angles at bus i and bus j respectively; b.sub.ij is a constant determined by the formula of b ij = - 0.5 r ij ( θ i - θ j ) + x ij r ij 2 + x ij 2  and is a pseudo branch susceptance of branch ij, where r.sub.ij and x.sub.ij are the resistance and reactance of branch ij respectively and are the given AC power network parameters; θ′.sub.i and θ′.sub.j are the initial buses voltage angles at bus i and bus j respectively and are the given AC power network parameters.

3. The method for obtaining a set of symmetric power transfer coefficients under simultaneous change of sources and loads in AC power networks according to claim 2, wherein the step of establishing symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles according to the linear function of the lossy branch transferred power in terms of buses voltage angles and the given bus injection powers of power sources and loads comprises the steps of: according to the linear function of the lossy branch transferred power in terms of buses voltage angles and the given bus injection powers of power sources and loads, establishing the linear function of particular bus injection power of power sources and loads in terms of buses voltage angles by the following equation: P Gi - P Di = - .Math. k = 1 , k i n b ik ( θ i - θ k ) where P.sub.Gi and −P.sub.Di are the bus injection powers of power sources and loads connected to bus i respectively; k is the number of each bus in the AC power network; k is natural number and k is equal to 1,2 . . . , n; θ.sub.k is the bus voltage angle at bus k; ik is the branch between bus i and bus k; b.sub.ik is the pseudo branch susceptance of branch ik determined by the formula of b ik = - 0.5 r ik ( θ i - θ k ) + x ik r ik 2 + x ik 2 ,  where r.sub.ik and x.sub.ik are the resistance and reactance of branch ik respectively and are the given AC power network parameters; θ′.sub.k is the initial bus voltage angle at bus k and is the given AC power network parameter; and according to the linear function of particular bus injection powers of power sources and loads in terms of buses voltage angles, establishing the symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles by the following matrix equation: [ P G 1 - P D 1 .Math. P Gi - P Di .Math. P Gn - P Dn ] = ( c ij ) [ θ 1 .Math. θ i .Math. θ n ] , c ij = { b ij , for j i - .Math. k = 1 , k i n b ik , for j = i where P.sub.G1 and −P.sub.D1 are the bus injection powers of power sources and loads connected to bus 1 respectively; where P.sub.Gn and −P.sub.Dn are the bus injection powers of power sources and loads connected to bus n respectively; θ.sub.1 and θ.sub.n are the bus voltage angles at bus 1 and bus n respectively; (c.sub.ij) is a pseudo bus susceptance matrix, c.sub.ij is row i and column j element of the pseudo bus susceptance matrix (c.sub.ij).

4. The method for obtaining a set of symmetric power transfer coefficients in AC power network under simultaneous change of sources and loads according to claim 3, wherein the step of establishing symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads according to the symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles comprises the step of: according to the symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles, establishing symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads by the following matrix equation: [ θ 1 .Math. θ i .Math. θ n ] = ( a ij ) [ P G 1 - P D 1 .Math. P Gi - P Di .Math. P Gn - P Dn ] , ( a ij ) = ( c ij ) + where (a.sub.ij) is the generalized pseudo bus reactance matrix of the AC power network and equals to Moore-Penrose pseudoinverse of pseudo bus susceptance matrix (c.sub.ij) of the AC power network, “+” is the operator of the Moore-Penrose pseudoinverse in the algebra.

5. The method for obtaining a set of symmetric power transfer coefficients under simultaneous change of sources and loads in AC power networks according to claim 4, wherein the step of obtaining the symmetric power transfer coefficients from the buses injection powers of power sources and loads to the lossy branch transferred power under simultaneous change of sources and loads according to the symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads and the linear function of the lossy branch transferred power in terms of buses voltage angles comprises the steps of: according to the symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads and the linear function of the lossy branch transferred power in terms of buses voltage angles, establishing the symmetric linear function of the lossy branch transferred power in terms of buses injection powers of power sources and loads by the following equation: P ij = - b ij .Math. k = 1 n ( a ik - a jk ) ( P Gk - P Dk ) where a.sub.ik is row i and column k element of the generalized pseudo bus reactance matrix (a.sub.ij) of the AC power network; a.sub.ik is row j and column k element of the generalized pseudo bus reactance matrix (a.sub.ij) of the AC power network; according to the symmetric linear function of the lossy branch transferred power in terms of buses injection powers of power sources and loads, obtaining the symmetric power transfer coefficients from the buses injection powers of power sources and loads to the lossy branch transferred power under simultaneous change of sources and loads by the following equations:
D′.sub.ij,k=−b.sub.ij(a.sub.ik−a.sub.jk)
D″.sub.ij,k=b.sub.ij(a.sub.ik−a.sub.jk) where D′.sub.ij,k is the symmetric power transfer coefficient from the bus injection power of the power sources connected to bus k to the lossy branch transferred power of branch ij; D″.sub.ij,k is the symmetric power transfer coefficient from the bus absorbing power of the load connected to bus k to the lossy branch transferred power of branch ij.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a flow chart of a method for obtaining a set of symmetric power transfer coefficients under simultaneous change of sources and loads in AC power network according to an embodiment of the present invention; and

(2) FIG. 2 is a structure diagram of a universal model in AC power network according to an embodiment of the present invention.

DESCRIPTION OF EMBODIMENTS

(3) To make the objectives, technical solutions, and advantages of the present invention clearer and more comprehensible, the following further describes the present invention in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiment described herein is merely used to explain the present invention but is not used to limit the present invention to it.

(4) An embodiment of the present invention provides a method for obtaining symmetric power transfer coefficients under simultaneous change of sources and loads in alternating current (AC) power network, which comprises:

(5) establishing a linear function of a lossy branch transferred power in terms of buses voltage angles according to given AC power network parameters, a nonlinear function of a branch transferred power in the AC power network, and operation features of the AC power network;

(6) establishing symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles according to the linear function of the lossy branch transferred power in terms of buses voltage angles and given buses injection powers of power sources and loads;

(7) establishing symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads according to the symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles;

(8) obtaining the symmetric power transfer coefficients from buses injection powers of power sources and loads to the lossy branch transferred power under simultaneous change of sources and loads according to the symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads and the linear function of the lossy branch transferred power in terms of buses voltage angles.

(9) In the embodiment of the present invention, the method comprises the steps of: firstly establishing a linear function of a lossy branch transferred power in terms of buses voltage angles according to a nonlinear function of a branch transferred power, parameters and operation features of the AC power network; then establishing symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles by combining the linear function of the lossy branch transferred power in terms of buses voltage angles and buses injection powers of power sources and loads, in turn establishing symmetric linear functions of buses voltage angles in terns of buses injection powers of power sources and loads; and finally obtaining the symmetric power transfer coefficients from buses injection powers of power sources and loads to the lossy branch transferred power under simultaneous change of sources and loads by using the linear functions just mentioned above. The obtained set of symmetric power transfer coefficients is unique, follows electric circuit laws, is applicable for the practical situations that power sources and loads change at the same time and transmission losses need considering in the AC power network, and truly reflects the substantive characteristics of the power transfer from bus injection powers to branch transferred powers, thus solve the problem that the obtained results of the existing methods for obtaining a set of power transfer coefficients in AC power network not only don't comply with the uniqueness theorem for electromagnetic field in circuits, but also don't follow electric circuit laws, are not applicable to the practical situation that the power network is lossy and bus injection powers of power sources and loads change at the same time, and thus fail to truly reflect the substantive characteristics of the power transfer from bus injection powers to branch transferred powers.

(10) FIG. 1 shows the flow chart of the method for obtaining a set of symmetric power transfer coefficients under simultaneous change of sources and loads in AC power network according to an embodiment of the present invention. For illustration purposes, only the parts related to the embodiment of the present invention are shown, which are described in detail as follow:

(11) In step S1, a linear function of a lossy branch transferred power in terms of buses voltage angles is established according to given AC power network parameters, a nonlinear function of a branch transferred power in the AC power network, and operation features of the AC power network.

(12) The step S1 comprises:

(13) according to the given AC power network parameters, the nonlinear function of the branch transferred power in the AC power network(i.e. the branch transferred power equals the bus voltage times complex conjugation of the branch current, which is the common AC power network knowledge), and operation features of the AC power network (i.e. the bus voltage angle difference across any branch in the AC power network is always quite small), establishing a linear function of the lossy branch transferred power in twins of buses voltage angles by the following equation:
P.sub.ij=−b.sub.ij(θ.sub.i−θ.sub.j)

(14) As shown in FIG. 2, i and j are the numbers of two arbitrary buses in the AC power network respectively, i and j are natural numbers and equal to 1, 2 . . . , n, n is the total number of buses in the AC power network (the buses herein refer to those operating independently in the AC power network), n is natural number and is the given AC power network parameter; ij is the branch between bus i and bus j; P.sub.ij is the lossy branch transferred power entering branch ij from bus i; θ.sub.i and θ.sub.j are voltage angles at bus i and bus j respectively; b.sub.ij is a constant determined by the formula of

(15) b ij = - 0.5 r ij ( θ i - θ j ) + x ij r ij 2 + x ij 2 ,
its property is similar to a branch susceptance, and b.sub.ij is called-pseudo branch susceptance of branch ij, where the r.sub.ij and x.sub.ij are the resistance and reactance of branch ij respectively and are the given AC power network parameters; θ′.sub.i and θ′.sub.j are the initial buses voltage angles at bus i and bus j respectively and are the given AC power network parameters.

(16) The lossy branch transferred power mentioned above refers to the loss-retained branch transferred power, the above function of the branch transferred power P.sub.ij contains inapparently the branch resistance r.sub.ij by the pseudo branch susceptance b.sub.ij, and the implied branch resistance reflects the transmission losses, therefore P.sub.ij is called lossy branch transferred power (i.e. the loss-retained branch transferred power).

(17) In step S2, symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles are established according to the linear function of the lossy branch transferred power in terms of buses voltage angles and given buses injection powers of power sources and loads.

(18) The step S2 comprises:

(19) according to the linear function of the lossy branch transferred power in terms of buses voltage angles and the given buses injection powers of power sources and loads, establishing the linear function of particular bus injection power of power sources and loads in terms of buses voltage angles by the following equation:

(20) P Gi - P Di = - .Math. k = 1 , k i n b ik ( θ i - θ k )

(21) As shown in FIG. 2, P.sub.Gi and −P.sub.Di are the bus injection powers of power sources and loads connected to bus i respectively; k is the number of each bus in the AC power network; k is natural number and k is equal to 1,2, . . . , n; θ.sub.k is the bus voltage angle at bus k; ik is the branch between bus i and bus k; b.sub.ik is the pseudo branch susceptance of branch ik determined by the formula of

(22) b ik = - 0.5 r ik ( θ i - θ k ) + x ik r ik 2 + x ik 2 ,
where the r.sub.ik and x.sub.ik are the resistance and reactance of branch ik respectively and are the given AC power network parameters; θ′.sub.k is the initial bus voltage angle at bus k and is the given AC power network parameter; and

(23) according to the linear function of particular bus injection powers of power sources and loads in terms of buses voltage angles, establishing the symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles by the following matrix equation:

(24) [ P G 1 - P D 1 .Math. P Gi - P Di .Math. P Gn - P Dn ] = ( c ij ) [ θ 1 .Math. θ i .Math. θ n ] , c ij = { b ij , for j i - .Math. k = 1 , k i n b ik , for j = i

(25) As shown in FIG. 2, where P.sub.G1 and −P.sub.D1 are the bus injection powers of power sources and loads connected to bus 1 respectively (the bus injection power of the power source refers to the time-average power of the power source entering bus 1, the bus injection power of the load refers to the minus of time-average power of the load absorbing from bus 1); where P.sub.Gn and −P.sub.Dn are the bus injection powers of power sources and loads connected to bus n respectively (the bus injection power of the power source refers to the time-average power of the power source entering bus n, the bus injection power of the load refers to the minus of time-average power of the load absorbing from bus n); θ.sub.1 and θ.sub.n are the bus voltage angles at bus 1 and bus n respectively; (c.sub.ij) is the AC power network complete bus susceptance matrix built by pseudo branch susceptances, its property is similar to the actual bus susceptance matrix and (c.sub.ij) is thus called pseudo bus susceptance matrix, c.sub.ij is row i and column j element of the pseudo bus susceptance matrix (c.sub.if).

(26) In the above matrix-equation-described linear functions of buses injection powers of power sources and loads in terms of buses voltage angles, the forms of the functions of buses injection powers of power sources and loads at all buses are identical, that is each bus is treated equally, so the linear functions are called symmetric linear functions.

(27) In step S3, symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads are established according to the symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles.

(28) The step S3 comprises:

(29) according to the symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles, establishing symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads by the following matrix equation:

(30) [ θ 1 .Math. θ i .Math. θ n ] = ( a ij ) [ P G 1 - P D 1 .Math. P Gi - P Di .Math. P Gn - P Dn ] , ( a ij ) = ( c ij ) +

(31) where (a.sub.ij) is the generalized pseudo bus reactance matrix of the AC power network and equals to Moore-Penrose pseudoinverse of pseudo bus susceptance matrix (c.sub.ij) of the AC power network, “+” is the operator of the Moore-Penrose pseudoinverse in the algebra.

(32) In the above matrix-equation-described linear functions of buses voltage angles in terms of buses injection powers of power sources and loads, the forms of the functions of buses voltage angles at all buses are identical, that is each bus is treated equally, so the linear functions are called symmetric linear functions.

(33) In step S4, the symmetric power transfer coefficients from buses injection powers of power sources and loads to the lossy branch transferred power under simultaneous change of sources and loads are obtained according to the symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads and the linear function of the lossy branch transferred power in terms of buses voltage angles.

(34) The step S4 comprises:

(35) according to the symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads and the linear function of the lossy branch transferred power in terms of buses voltage angles, establishing the symmetric linear function of the lossy branch transferred power in terms of buses injection powers of power sources and loads by the following equation:

(36) P ij = - b ij .Math. k = 1 n ( a ik - a jk ) ( P Gk - P Dk )

(37) where a.sub.ik is row i and column k element of the generalized pseudo bus reactance matrix (a.sub.ij) of the AC power network; a.sub.jk is row j and column k element of the generalized pseudo bus reactance matrix (a.sub.ij) of the AC power network;

(38) according to the symmetric linear function of the lossy branch transferred power in terms of buses injection powers of power sources and loads, obtaining the symmetric power transfer coefficients from the buses injection powers of power sources and loads to the lossy branch transferred power under simultaneous change of sources and loads by the following equations:
D′.sub.ij,k=−b.sub.ij(a.sub.ik−a.sub.jk)
D″.sub.ij,k=b.sub.ij(a.sub.ik−a.sub.jk)

(39) where D′.sub.ij,k is the symmetric power transfer coefficient from the bus injection power of the power sources connected to bus k to the lossy branch transferred power of branch ij; D″.sub.ij,k is the symmetric power transfer coefficient from the bus absorbing power of the loads connected to bus k to the lossy branch transferred power of branch ij.

(40) D′.sub.ij,k and D″.sub.ij,k are obtained according to the symmetric linear function of the lossy branch transferred power (i.e. the loss-retained branch transferred power) in terms of buses injection powers of power sources and loads, and the symmetric linear function not only considers simultaneous change of sources and loads, but also is identical for all branches, so D′.sub.ij,k and D″.sub.ij,k are called the symmetric power transfer coefficients under simultaneous change of sources and loads.

(41) Using the above sequence of functions-calculates the symmetric power transfer coefficients for each couple of buses and branches then the symmetric power transfer coefficients from all buses injection powers of power sources and loads to each lossy branch transferred power under simultaneous change of sources and loads are obtained. The symmetric power transfer coefficients not only comply with the uniqueness theorem for electromagnetic field in circuits, but also follow electric circuit laws, are applicable to the practical situation that the power network is lossy and bus injection powers of power sources and loads change at the same time, and thus can truly reflect the substantive characteristics of the power transfer from bus injection powers to branch transferred powers.

(42) In the embodiment of the present invention, the method comprises the steps of: firstly establishing a linear function of a lossy branch transferred power in terms of buses voltage angles according to a nonlinear function of a branch transferred power, parameters and operation features of the AC power network; then establishing symmetric linear functions of buses injection powers of power sources and loads in terms of buses voltage angles by combining the linear function of the lossy branch transferred power in terms of buses voltage angles and buses injection powers of power sources and loads, in turn establishing symmetric linear functions of buses voltage angles in terms of buses injection powers of power sources and loads; and finally obtaining the symmetric power transfer coefficients from buses injection powers of power sources and loads to the lossy branch transferred power under simultaneous change of sources and loads by using the linear functions just mentioned above. The obtained set of symmetric power transfer coefficients is unique, follows electric circuit laws, is applicable for the practical situations that power sources and loads change at the same time and transmission losses need considering in the AC power network, and truly reflects the substantive characteristics of the power transfer from bus injection powers to branch transferred powers, thus solve the problem that the obtained results of the existing methods for obtaining a set of power transfer coefficients in AC power network not only don't comply with the uniqueness theorem for electromagnetic field in circuits, but also don't follow electric circuit laws, are not applicable to the practical situation that the power network is lossy and bus injection powers of power sources and loads change at the same time, and thus fail to truly reflect the substantive characteristics of the power transfer from bus injection powers to branch transferred powers.

(43) The foregoing descriptions are merely exemplary embodiment of the present invention, but are not intended to limit the present invention to it. Any modification, equivalent replacement, or improvement made without departing from the spirit and principle of the present invention shall fall within the protection scope of the present invention.