ONLINE VOLTAGE CONTROL METHOD FOR COORDINATING MULTI-TYPE REACTIVE POWER RESOURCES

20220385072 · 2022-12-01

    Inventors

    Cpc classification

    International classification

    Abstract

    An online voltage control method for coordinating multi-type reactive power resources is provided. First, a linearized power flow equation of branch reactive power is established, and an online voltage control model of multi-type reactive power resources including an objective function and constraint conditions is constructed. The constrain conditions includes generator reactive power constraints, reactive power compensator constraints, transformer tap position constraints, a nodal reactive power balance constraint, and slack contained nodal voltage constraints. Then, an optimization result of voltage control is obtained by solving the model. The method makes full use of reactive voltage operation characteristics of a power grid, constructs a practical online solution model for reactive voltage control of large power grid of coordinating multiple reactive power resources, and under a condition of acceptable accuracy loss, takes in account safety of power grid operation, economy of reactive power resource actions and high reliability of online operation.

    Claims

    1. An online voltage control method for coordinating multi-type reactive power resources, comprising steps: (1) establishing a variable set Ω of a base-state operation point model of a power system, wherein
    Ω={V.sub.n,V.sub.n.sup.slack.sup.up,V.sub.n.sup.slack.sup.down,{tilde over (V)}.sub.n,V.sub.i.sup.G,Q.sub.ij.sup.b,Q.sub.i.sup.G,Q.sub.n,Q.sub.n.sup.cp,Q.sub.n.sup.un,Q.sub.n.sup.ld,N.sub.unit.sup.cp,N.sub.unit.sup.rc,μ.sub.unit,B.sub.i.sup.cp,B.sub.j.sup.rc,t.sub.i.sup.adj_up,t.sub.i.sup.adj_down} where, V.sub.n is a voltage after adjustment of a node n, V.sub.n.sup.slack_up and V.sub.n.sup.slack_down respectively are an up slack variable and a down slack variable of the voltage after adjustment of the node n, {tilde over (V)}.sub.n is an adjustment voltage value after slack of the node n, V.sub.i.sup.G is a voltage of a generator i, Q.sub.ij.sup.b is reactive power flowing into a branch b with ports being nodes i and j from the port being the node i, Q.sub.i.sup.G is reactive power of the generator i, Q.sub.n is reactive power of the node n, Q.sub.n.sup.cp is reactive power supplied by a reactive power compensator connected to the node n, Q.sub.n.sup.un is reactive power supplied by a generator connected to the node n, Q.sub.n.sup.ld is reactive power absorbed by a load connected to the node n , N.sub.unit.sup.cp is a number of capacitance compensators in a state of being put in operation under a controller unit, N.sub.unit.sup.rc is a number of inductance compensators in a state of being put in operation under the controller unit, μ.sub.unit is a binary variable of the controller unit, B.sub.i.sup.cp is a Boolean variable of a capacitance compensator i to indicate the capacitance compensator numbered with i whether changes its operation state, B.sub.j.sup.rc is a Boolean variable of an inductance compensator j to indicate the inductance compensator numbered with j whether changes its operation state, t.sub.i.sup.adj_up is an upwards adjusted tap position of a transformer i, and t.sub.i.sup.adj_down is a downwards adjusted tap position of the transformer i; (2) establishing a linearized equation of reactive power of branch, comprising: for each the branch b, calculating a Jacobi matrix of reactive power of the branch to obtain the following values: Q i j b V i , Q i j b V j , Q j i b V i , Q j i b V j , Q i j b t , Q j i b t ; where, i and j are node numbers respectively corresponding to two ports of the branch b, Q i j b V i  is a partial derivative of reactive power of the branch b flowing out from the port node j relative to a voltage V.sub.i after adjustment of the port node i, Q i j b V j  is a partial derivative of the reactive power of the branch b flowing out from the port node j relative to a voltage V.sub.j after adjustment of the port node j, Q j i b V i  is a partial derivative of reactive power of the branch b flowing out from the port node i relative to the voltage V.sub.i after adjustment of the port node i, Q j i b V j  is a partial derivative of the reactive power of the branch b flowing out from the port node i relative to the voltage V.sub.j after adjustment of the port node j, Q i j b t  is a partial derivative of the reactive power of the branch b flowing out from the port node j relative to a tap position t of a transformer connected to the branch b, and Q j i b t  is a partial derivative of the reactive power of the branch b flowing out from the port node i relative to the tap position t of the transformer connected to the branch b; an expression of the linearized equation of reactive power of the branch b is as follows: Q i j b = Q i j 0 b + Q j i b V i .Math. ( V i - V i 0 ) + Q j i b V j .Math. ( V j - V j 0 ) + Q j i b t .Math. ( t - t 0 ) where, Q.sub.ij.sup.b is the reactive power flowing into the branch b from the port node i, Q.sub.ij0.sup.b is an initial reactive power flowing into the branch b from the port node i, V.sub.i0 is an original voltage of the node i, V.sub.j0 is an original voltage of the node j, t is a tap position after adjustment of the transformer connected to the branch b, and t.sub.0 is an original tap position of the transformer connected to the branch b; (3) establishing an online voltage control model of multi-type reactive power resources, wherein the online voltage control model comprises an objective function and constraint conditions, specifically comprising the following steps: (3-1) establishing an expression of the objective function of the online voltage control model as follows:
    minC.sub.total,
    wherein, C.sub.total=Σ(Q.sub.n.sup.cp,op.Math.C.sub.n.sup.cp,op)+Σ(Q.sub.n.sup.un,adj_up.Math.C.sub.n.sup.un,adj_up+Q.sub.n.sup.un,adj_down.Math.C.sub.n.sup.un,adj_down)+Σ(t.sub.i.sup.adj_up+t.sub.i.sup.adj_down).Math.C.sub.i.sup.t,adj+Σ(V.sub.n.sup.slack_up.sup.2+V.sub.n.sup.slack_down.sup.2).Math.C.sub.p, where, C.sub.total is a sum of a total cost and a penalty term of regulation, Q.sub.n.sup.cp,op is a 0-1 variable of the reactive power compensator connected to the node n whether is put in operation, C.sub.n.sup.cp,op is a cost of the reactive power compensator connected to the node n in operation, Q.sub.n.sup.un,adj_up is an amount of upward adjustment of output reactive power of the generator connected to node n, C.sub.n.sup.un,adj_up is a cost of upward adjustment of output reactive power of the generator connected to node n, Q.sub.n.sup.un,adj_down is an amount of downward adjustment of output reactive power of the generator connected to node n, C.sub.n.sup.un,adj_down is a cost of downward adjustment of output reactive power of the generator connected to node n, C.sub.i.sup.t,adj is a cost per unit tap position of adjustment of the transformer i, Σ(V.sub.n.sup.slack_up.sup.2+V.sub.n.sup.slack_down.sup.2).Math.C.sub.p is the penalty term, and C.sub.p is a penalty coefficient of slack variable; (3-2) establishing the constraint conditions as follows: (3-2-1) generator reactive power constraints: when the Q.sub.i.sup.G is adjustable, there is a constraint condition:
    Q.sub.i.sup.G_min≤Q.sub.i.sup.G≤Q.sub.i.sup.G_max, ∀.sub.i∈I.sup.G, where, Q.sub.i.sup.G_min is a lower limit of the reactive power of the generator i, Q.sub.i.sup.G_max is an upper limit of the reactive power of the generator i, and I.sup.G is a collection of generators; when the Q.sub.i.sup.G is non-adjustable and the generator i is a PQ node whose real power P and reactive power Q are specified, there is a constraint condition: Q.sub.i.sup.G=Q.sub.i0.sup.G; when the Q.sub.i.sup.G is non-adjustable and the generator i is a PV node whose real power P and a voltage magnitude V are specified, there is a constraint condition: V.sub.i.sup.G=V.sub.i0.sup.G, where, Q.sub.i0.sup.G and V.sub.i0.sup.G respectively are original reactive power and an original voltage of the generator i; (3-2-2) reactive power compensator constraints:
    N.sub.unit.sup.cp=N.sub.unit.sup.cp0+Σ.sub.i∈I.sub.0B.sub.i.sup.cp−Σ.sub.i∈I.sub.1B.sub.i.sup.cp, where, the controller is a minimum control unit of a group of reactive power compensators, the group of reactive power compensators comprise capacitance compensators and inductance compensators, I.sub.0 is a capacitance compensator set in exit state under the controller unit, I.sub.1 is a capacitance compensator set in the state of being put in operation under the controller unit, N.sub.unit.sup.cp0 is a number of the capacitance compensators in the state of being put in operation under the controller unit at an initial stage;
    N.sub.unit.sup.rc=N.sub.unit.sup.rc0+Σ.sub.j∈J.sub.0B.sub.j.sup.rc−Σ.sub.j∈J.sub.1B.sub.j.sup.rc, where, J.sub.0 is an inductance compensator set in exit state under the controller unit, J.sub.1 is an inductance compensator set in the state of being put in operation under the controller unit, N.sub.unit.sup.rc0 is a number of the inductance compensators in the state of being put in operation under the controller unit at the initial stage; using a big-M method to establish a constraint as follows: { 0 N unit cp μ u n i t N u n i t s u m 0 N u n i t r c ( 1 - μ u n i t ) N u n i t s u m μ u n i t = 0 , 1 where, N.sub.unit.sup.sum is a number of all the reactive power compensators under the controller unit; (3-2-3) slack contained nodal voltage constraints: { V ~ n = V n + V n slack_up - V n slack_down V n slack_up , V n slack_down 0 , when a voltage on the node n has upper and lower limit constraints, there is a constraint condition:
    V.sub.n.sup.min≤{tilde over (V)}.sub.n≤V.sub.n.sup.max, when the voltage on the node n has an objective voltage value, there is a constraint condition:
    {tilde over (V)}.sub.n=V.sub.n.sup.obj, where, V.sub.n.sup.min is an allowable minimum value of the voltage on the node n, V.sub.n.sup.max is an allowable maximum value of the voltage on the node n, and V.sub.n.sup.obj is the objective voltage value of the node n; (3-2-4) transformer tap position constraints: when a tap of the transformer i is adjustable, there is a constraint condition: { t i min t i 0 + t i adj_up - t i adj_down t i max 0 t i adj_up , t i adj_down 2 where, t.sub.i.sup.min is a minimum value of the tap position of the transformer i, t.sub.i.sup.max is a maximum value of the tap position of the transformer i, and t.sub.i.sup.0 is an original value of the tap position of the transformer i; (3-2-5) a nodal reactive power balance constraint:
    Q.sub.n=ΣQ.sub.n.sup.cp+ΣQ.sub.n.sup.un−ΣQ.sub.n.sup.ld−Σ.sub.jQ.sub.nj.sup.b=0, where, ΣQ.sub.n.sup.cp is a sum of reactive power flowing in the node n from all the reactive power compensators connected to the node n, ΣQ.sub.n.sup.un is a sum of reactive power flowing in the node n from all generators connected to the node n, ΣQ.sub.n.sup.ld is a sum of reactive power of flowing from the node n into all loads connected to the node n, Σ.sub.jQ.sub.nj.sup.b is a sum of reactive power of flowing from the node n into all branches each with the node n as its port, and Q.sub.nj.sup.b is calculated by the linearized equation of reactive power of branch established in the step (2); (4) solving the online voltage control model of multi-type reactive power resources established in the step (3) to obtain optimal values of respective optimization variables contained in the variable set Ω, thereby completing the voltage control.

    Description

    DETAILED DESCRIPTION OF EMBODIMENT

    [0032] An embodiment of the invention provides an online voltage control method for coordinating multi-type reactive power resources. The online voltage control method may include the following step (1) through step (4).

    [0033] Step (1), establishing a variable set Ω of a base-state operation point model of a power system,


    Ω={V.sub.n,V.sub.n.sup.slack_up,V.sub.n.sup.slack_down,{tilde over (V)}.sub.n,V.sub.i.sup.G,Q.sub.ij.sup.b,Q.sub.i.sup.G,Q.sub.n,Q.sub.n.sup.cp,Q.sub.n.sup.un,Q.sub.n.sup.ld,N.sub.unit.sup.cp,N.sub.unit.sup.rc,μ.sub.unit,B.sub.i.sup.cp,B.sub.j.sup.rc,t.sub.i.sup.adj_up,t.sub.i.sup.adj_down}

    where, V.sub.n is a voltage after adjustment of a node n, V.sub.n.sup.slack_up and V.sub.n.sup.slack_down respectively are an up slack variable and a down slack variable of the voltage after adjustment of the node n, {tilde over (V)}.sub.n is an adjustment voltage value after slack of the node n, V.sub.i.sup.G is a voltage of a generator i, Q.sub.ij.sup.b is reactive power flowing into a branch b with ports being nodes i and j from the port being the node i, Q.sub.i.sup.G is reactive power of the generator i, Q.sub.n is reactive power of the node n, Q.sub.n.sup.cp is reactive power supplied by a reactive power compensator connected to the node n, Q.sub.n.sup.un is reactive power supplied by a generator connected to the node n, Q.sub.n.sup.ld is reactive power absorbed by a load connected to the node n, N.sub.unit.sup.cp is a number of capacitance compensators in a state of being put in operation under a controller unit, N.sub.unit.sup.rc is a number of inductance compensators in a state of being put in operation under the controller unit, μ.sub.unit is a binary variable of the controller unit, B.sub.i.sup.cp is a Boolean variable of a capacitance compensator i to indicate the capacitance compensator numbered with i whether changes its operation state, B.sub.j.sup.rc is a Boolean variable of an inductance compensator j to indicate the inductance compensator numbered with j whether changes its operation state, t.sub.i.sup.adj_up is an upwards adjusted tap position of a transformer i, and t.sub.i.sup.adj_down is a downwards adjusted tap position of the transformer i.

    [0034] Step (2), establishing a linearized equation of reactive power of branch, specifically: [0035] for each the branch b, calculating a Jacobi matrix of reactive power of the branch to obtain the following values:

    [00012] Q ij b V i , Q ij b V j , Q ji b V i , Q ji b V j , Q ij b t , Q ji b t ;

    where, i and j are node numbers respectively corresponding to two ports of the branch b,

    [00013] Q ij b V i

    is a partial derivative of reactive power of the branch b flowing out from the port node j relative to a voltage V.sub.i after adjustment of the port node i,

    [00014] Q ij b V j

    is a partial derivative of the reactive power of the branch b flowing out from the port node j relative to a voltage V.sub.j after adjustment of the port node j,

    [00015] Q ji b V i

    is a partial derivative of reactive power of the branch b flowing out from the port node i relative to the voltage V.sub.i after adjustment of the port node i,

    [00016] Q ji b V j

    is a partial derivative of the reactive power of the branch b flowing out from the port node i relative to the voltage V.sub.j after adjustment of the port node j,

    [00017] Q ij b t

    is a partial derivative of the reactive power of the branch b flowing out from the port node j relative to a tap position t of a transformer connected to the branch b, and

    [00018] Q ji b t

    is a partial derivative of the reactive power of the branch b flowing out from the port node i relative to the tap position t of the transformer connected to the branch b;

    [0036] Because voltage variations of respective nodes in the power grid are relatively small, variations of reactive power caused by the variations of nodal voltages can be approximately linear. Similarly, a change of the transformer tap position generally does not exceed two positions (gears), and a caused variation of reactive power can also be approximately linear. Therefore, a relationship between the reactive power on the branch b and voltages at both ends of the branch and the transformer tap position can be obtained as follows:

    [00019] Q ij b = Q ij 0 b + Q ji b V i .Math. ( V i - V i 0 ) + Q ji b V j .Math. ( V j - V j 0 ) + Q ji b t .Math. ( t - t 0 )

    where, Q.sub.ij.sup.b is the reactive power flowing into the branch b with ports being i and j from the port node i (i.e., the reactive power flowing out of the node i), Q.sub.ij0.sup.b is an initial reactive power flowing into the branch b with the ports being i and j from the port node i, V.sub.i is a voltage value after adjustment of the node i, V.sub.i0 is an original voltage of the port node i, V.sub.j is a voltage value after adjustment of the node j, V.sub.j0 is an original voltage of the port node j,

    [00020] Q i j b t

    is the partial derivative of the reactive power of the branch b flowing out from the port node j relative to the tap position t of the transformer connected to the branch b, t is a tap position after adjustment of the transformer connected to the branch b, and t.sub.0 is an original tap position of the transformer connected to the branch b.

    [0037] Step (3), establishing an online voltage control model of multi-type reactive power resources, the online voltage control model including an objective function and constraint conditions; specifically, it may include the following sub-step (3-1) and sub-step (3-2).

    [0038] Sub-step (3-1), establishing an expression of the objective function of the online voltage control model as follows:


    minC.sub.total,


    wherein, C.sub.total=Σ(Q.sub.n.sup.cp,op.Math.C.sub.n.sup.cp,op)+Σ(Q.sub.n.sup.un,adj_up.Math.C.sub.n.sup.un,adj_up+Q.sub.n.sup.un,adj_down.Math.C.sub.n.sup.un,adj_down)+Σ(t.sub.i.sup.adj_up+t.sub.i.sup.adj_down).Math.C.sub.i.sup.t,adj+Σ(V.sub.n.sup.slack_up.sup.2+V.sub.n.sup.slack_down.sup.2).Math.C.sub.p,

    where, C.sub.total is a sum of a total cost and a penalty term of regulation (adjustment), Q.sub.n.sup.cp,op is a 0-1 variable of the reactive power compensator connected to the node n whether is put in operation, C.sub.n.sup.cp,op is a cost of the reactive power compensator connected to the node n in operation, Q.sub.n.sup.un,adj_up is an amount of upward adjustment of output reactive power of the generator connected to node n, C.sub.n.sup.un,adj_up is a cost of upward adjustment of output reactive power of the generator connected to node n, Q.sub.n.sup.un,adj_down is an amount of downward adjustment of output reactive power of the generator connected to node n, C.sub.n.sup.un,adj_down is a cost of downward adjustment of output reactive power of the generator connected to node n, C.sub.i.sup.t,adj is a cost per unit tap position of adjustment of the transformer i, Σ(V.sub.n.sup.slack_up.sup.2+V.sub.n.sup.slack_down.sup.2).Math.C.sub.p is the penalty term, and C.sub.p is a penalty coefficient of slack variable.

    [0039] The penalty coefficient C.sub.p generally takes a large number such as 10.sub.8, and therefore, if the final optimization objective result has a relatively large order of magnitude, it indicates that the optimization objective has not been achieved, and a feasible solution close to the objective value is obtained

    [0040] Sub-step (3-2), establishing the constraint conditions shown by (3-2-1) through (3-2-5) as follows.

    [0041] (3-2-1) generator reactive power constraints: [0042] when the Q.sub.i.sup.G is adjustable, there is a constraint condition:


    Q.sub.i.sup.G_min≤Q.sub.i.sup.G≤Q.sub.i.sup.G_max, ∀.sub.i∈I.sup.G,

    where, Q.sub.i.sup.G_min is a lower limit of the reactive power of the generator i, Q.sub.i.sup.G_max is an upper limit of the reactive power of the generator i, and I.sup.G is a collection of generators; [0043] when the Q.sub.i.sup.G is non-adjustable and the generator i is a PQ node (i.e., a node with specified rea power P and reactive power Q), there is a constraint condition: Q.sub.i.sup.G=Q.sub.i0.sup.G, [0044] when the Q.sub.i.sup.G is non-adjustable and the generator i is a PV node (i.e., a node with specified real power P and voltage magnitude V), there is a constraint condition: V.sub.i.sup.G=V.sub.i0.sup.G, [0045] where, Q.sub.i0.sup.G and V.sub.i0.sup.G respectively are original reactive power and an original voltage of the generator i.

    [0046] (3-2-2) reactive power compensator constraints:

    [0047] The controller is defined as the minimum control unit of a group of reactive power compensators (including capacitance compensators and inductance compensators). Generally, one plant with reactive power compensators is one controller. The capacitance compensator and the inductance compensator cannot be put into operation at the same time under the same controller, otherwise oscillation will occur in the power grid;


    N.sub.unit.sup.cp=N.sub.unit.sup.cp0+Σ.sub.i∈I.sub.0B.sub.i.sup.cp−Σ.sub.i∈I.sub.1B.sub.i.sup.cp,

    where, I.sub.0 is a capacitance compensator set in exit state under the controller unit, I.sub.1 is a capacitance compensator set in the state of being put in operation under the controller unit, N.sub.unit.sup.cp0 is a number of the capacitance compensators in the state of being put in operation under the controller unit at an initial stage, B.sub.i.sup.cp is the Boolean variable of the capacitance compensator i to indicate the capacitance compensator numbered with i whether changes its operation state, and if the capacitance compensator numbered with i is non-adjustable, B.sub.i.sup.cp=0;


    N.sub.unit.sup.rc=N.sub.unit.sup.rc0+Σ.sub.j∈J.sub.0B.sub.j.sup.rc−Σ.sub.j∈J.sub.1B.sub.j.sup.rc,

    where, J.sub.0 is a reactance compensator set in exit state under the controller unit, J.sub.1 is a reactance compensator set in the state of being put in operation under the controller unit, N.sub.unit.sup.rc0 is a number of the reactance compensators in the state of being put in operation under the controller unit at the initial stage, B.sub.j.sup.rc is the Boolean variable of the inductance compensator j to indicate the inductance compensator numbered with j whether changes its operation state, and if the inductance compensator numbered with j in non-adjustable, B.sub.j.sup.rc=0;

    [0048] In the same controller, the capacitance compensator and the inductance compensator cannot be put into operation at the same time, otherwise oscillation will occur; and therefore using a big-M method to establish a constraint as follows:

    [00021] { 0 N u n i t cp μ u n i t N u n i t s u m 0 N u n i t r c ( 1 - μ u n i t ) N u n i t s u m μ u n i t = 0 , 1

    where, N.sub.unit.sup.sum is a number of all the reactive power compensators under the controller unit, and μ.sub.unit is a binary variable; as a result, it can be realized that at most one of N.sub.unit.sup.cp and N.sub.unit.sup.rc is greater than 0.

    [0049] (3-2-3) slack contained nodal voltage constraints:

    [00022] { V ~ n = V n + V n slack_up - V n slack_down V n slack_up , V n slack_down 0 , [0050] when a voltage on the node n has upper and lower limit constraints, there is a constraint condition: V.sub.n.sup.min≤{tilde over (V)}.sub.n≤V.sub.n.sup.max, [0051] when the voltage on the node n has an objective voltage value, there is a constraint condition: {tilde over (V)}.sub.n=V.sub.n.sup.obj, [0052] where, V.sub.n.sup.min is an allowable minimum value of the voltage on the node n, V.sub.n.sup.max is an allowable maximum value of the voltage on the node n, and V.sub.n.sup.obj is the objective voltage value of the node n; after adding the slack variables, an unreachable optimization objective can still be solved, and a result relatively close to the optimization objective can be obtained.

    [0053] (3-2-4) transformer tap position constraints: [0054] when a tap of the transformer i is adjustable, there are constraint conditions:

    [00023] { t i min t i 0 + t i adj_up - t i adj_down t i max 0 t i adj_up , t i adj_down 2

    where, t.sub.i.sup.min is a minimum value of the tap position of the transformer i, t.sub.i.sup.max is a maximum value of the tap position of the transformer i, t.sub.i.sup.0 is an original value of the tap position of the transformer i, t.sub.i.sup.adj_up is the upwards adjusted tap position of the transformer i, and t.sub.i.sup.adj_down is the downwards adjusted tap position of the transformer i. Since the tap position adjustment cannot adjust multiple gears at one time in actual, both t.sub.i.sup.adj_up and t.sub.i.sup.adj_down are set as cannot exceed 2.

    [0055] (3-2-5) a nodal reactive power balance constraint: [0056] for any one node, input reactive power should be equal to output reactive power, and thus: Q.sub.n=ΣQ.sub.n.sup.cp+ΣQ.sub.n.sup.un−ΣQ.sub.n.sup.ld−Σ.sub.jQ.sub.nj.sup.b=0,
    where, Q.sub.n is the reactive power of the node n, ΣQ.sub.n.sup.cp is a sum of reactive power flowing in the node n from all the reactive power compensators connected to the node n, ΣQ.sub.n.sup.un is a sum of reactive power flowing in the node n from all generators connected to the node n, ΣQ.sub.n.sup.ld is a sum of reactive power of flowing from the node n into all loads connected to the node n, Σ.sub.jQ.sub.nj.sup.b is a sum of reactive power of flowing from the node n into all branches each with the node n as its port, and Q.sub.nj.sup.b is calculated by the linearized equation of reactive power of branch established in the step (2).

    [0057] Step (4), optimizing and solving: [0058] the established above online voltage control model of multi-type reactive power resources can be solved by Cplex, and optimal solutions of respective optimization variables contained in the variable set Ω are finally obtained, and the voltage control is completed.