Detailed Description
The following describes the embodiments of the present invention in further detail with reference to the drawings.
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
In order to fully exert the adjustment and complementary characteristics of controllable resources such as a grid frame, distributed energy sources, energy storage and the like, the invention provides a power supply capacity-considered power distribution network optimization scheduling method, which establishes an optimal operation scheduling model taking the maximum renewable energy source utilization rate and the minimum line switching action number as multiple targets and a checking model capable of guaranteeing the power supply capacity of a system, and realizes optimal scheduling of the controllable resources in the system while guaranteeing the power supply reliability of the system through the coordinated use of the optimal operation scheduling model and the checking model, as shown in a figure 1, the method comprises the following steps:
Step 101, determining the optimal output of each power supply node in the power distribution network and the optimal switching state of each circuit in the power distribution network according to the load demand of each load node in the power distribution network;
Step 102, respectively controlling the output of each power supply node in the power distribution network and the switching state of each circuit in the power distribution network to be the optimal output and the optimal switching state;
The power supply nodes are new energy nodes or distributed power supply nodes.
Specifically, the optimal output includes an optimal active output and an optimal reactive output, and the step 101 includes:
Substituting the load demands of all load nodes in the power distribution network into a pre-constructed optimal operation scheduling model, solving the pre-constructed optimal operation scheduling model, and obtaining the optimal output of all power supply nodes in the power distribution network and the optimal switching state of all lines in the power distribution network;
Amplifying the load demand of each load node in the power distribution network, substituting the amplified load demand of each load node in the power distribution network, the optimal active output of each new energy node in the power supply node in the power distribution network and the optimal switching state of each line in the power distribution network into a pre-established test model, solving the pre-established test model, and obtaining an objective function value output by the test model;
And c, judging whether the objective function value output by the test model is 0, if so, outputting the optimal active power output and the optimal reactive power output of each power supply node in the power distribution network and the optimal switching state of each circuit in the power distribution network, otherwise, adding test constraint conditions in the pre-built optimal operation scheduling model, and returning to the step a.
Further, an objective function of the optimal operation scheduling model is determined according to the following formula:
minF=w1·f1+w2·f2
Wherein F is an objective function value of an optimal operation scheduling model, w 1 is a weight corresponding to the new energy utilization rate, F 1 is the reciprocal of the new energy utilization rate of the power distribution network in a control period, w 2 is a weight corresponding to the line switch action times of the power distribution network, and F 2 is the line switch action times of the power distribution network in the control period;
The new energy utilization rate f 1 of the power distribution network in the control period is determined according to the following formula:
In the formula, For the output of the xth new energy node in the power nodes of the power distribution network consumed at the t-th moment of the control period, P NE,x (t) is the actual output of the xth new energy node in the power nodes of the power distribution network at the t-th moment of the control period (t epsilon (1-N T),NT is the total time of the control period, x epsilon (1-N J),NJ is the total number of new energy nodes in the power nodes of the power distribution network);
the line switch action times f 2 of the power distribution network in the control period are determined according to the following steps:
Wherein mu ij,t is the switching state of the line ij in the power distribution network at the t-th moment of the control period, mu ij,t-1 is the switching state of the line ij in the power distribution network at the t-1 th moment of the control period, ij epsilon (1-N B),NB is the total number of lines in the power distribution network;
If the switching state of the line ij in the power distribution network at the t-th moment of the control period is on, μ ij,t =1, otherwise μ ij,t =0;
If the switching state of the line ij in the power distribution network at the t-1 time of the control period is on, μ ij,t-1 =1, otherwise μ ij,t-1 =0.
Further, the interactive output constraint condition of the active power distribution network and the large power grid of the objective function of the optimal operation scheduling model is determined according to the following formula:
In the formula, To control the equivalent current of the active power output provided by the large grid to node i in the distribution network at time t of the cycle,For a preset minimum value of the equivalent current of the active output provided by the large grid to node i in the distribution network,For a preset maximum value of the equivalent current of the active output provided by the large grid to node i in the distribution network,For controlling the equivalent current of reactive power output provided by the large grid to node i in the distribution network at time t of the period,For a preset minimum value of the equivalent current of reactive power output supplied by the large grid to node i in the distribution network,The method comprises the steps that the preset maximum value of the equivalent current of reactive power output provided by a large power grid to a node i in a power distribution network is set;
Determining the output linear power flow constraint condition of the objective function of the optimal operation scheduling model according to the following steps:
Wherein G ij is the conductance of a line ij in the power distribution network, B ij is the susceptance of the line ij in the power distribution network, beta ij,t is the product of the switching state of the line ij in the power distribution network at the t moment of a control period and the voltage imaginary part of a node j in the power distribution network, beta ji,t is the product of the switching state of the line ij in the power distribution network at the t moment of the control period and the voltage imaginary part of the node i in the power distribution network, alpha ij,t is the product of the switching state of the line ij in the power distribution network at the t moment of the control period and the voltage real part of the node j in the power distribution network, alpha ji,t is the product of the switching state of the line ij in the power distribution network at the t moment of the control period and the voltage real part of the node i in the power distribution network, To control the equivalent current of the active power output provided by the large grid to node i in the distribution network at time t of the cycle,For controlling the equivalent current of the active output of the distributed power supply input by the node i in the power distribution network at the t-th moment of the period,To control the real part of the voltage at node i in the distribution network at time t of the cycle,For controlling the equivalent current of the active power output of the new energy input by the node i in the power distribution network at the t-th moment of the period,To control the voltage imaginary part of node i in the distribution network at time t of the cycle,For controlling the equivalent current of reactive power output provided by the large grid to node i in the distribution network at time t of the period,For controlling the equivalent current of reactive power output of the distributed power supply input by the node i in the power distribution network at the t-th moment of the period,For controlling the equivalent current of the reactive power output of the new energy input by the node i in the power distribution network at the t-th moment of the period,For the equivalent current of the active load demand of node i in the power distribution network at the t-th moment of the control period, G i,t is the equivalent conductance of the load demand of node i in the power distribution network at the t-th moment of the control period, B i,t is the equivalent susceptance of the load demand of node i in the power distribution network at the t-th moment of the control period,For the equivalent current required by reactive load of a node i in the power distribution network at the t-th moment of the control period, j epsilon (1-n), wherein n is the total number of other nodes except the node i in the power distribution network;
Wherein, the To control the real part of the voltage at node j in the distribution network at time t of the cycle,A maximum value is preset for the real part of the voltage at node j in the distribution network,A minimum value is preset for the real part of the voltage at node j in the distribution network,To control the voltage imaginary part of node j in the distribution network at time t of the cycle,A maximum value is preset for the voltage imaginary part of the node j in the distribution network,Presetting a minimum value for the voltage imaginary part of a node j in the power distribution network;
determining branch current constraint conditions of an objective function of an optimal operation scheduling model according to the following steps:
In the formula, To control the real current of the line ij in the distribution network at the t-th moment of the cycle,To control the imaginary current of the line ij in the distribution network at the t-th moment of the period,Maximum current allowed on line ij in the distribution network;
Determining a voltage constraint condition of an objective function of an optimal operation scheduling model of the power distribution network according to the following steps:
In the formula, For the real part of the voltage of the node i in the power distribution network at the t-th moment of the control period, V i max is the maximum value of the voltage of the node i in the power distribution network, and V i min is the minimum value of the voltage of the node i in the power distribution network;
Determining a distributed power constraint condition of an objective function of an optimal operation scheduling model according to the following steps:
In the formula, For controlling the equivalent current of the active output of the distributed power supply input by the node i in the power distribution network at the t-th moment of the period,The minimum value of the equivalent current of the active output of the distributed power source input to the node i in the power distribution network,The maximum value of equivalent current of active output of the distributed power supply input to a node i in the power distribution network;
determining a network radial constraint condition of an objective function of the optimal operation scheduling model according to the following steps:
Wherein τ ij is an output flow direction state variable of a line ij in the power distribution network, τ ji is an output flow direction state variable of a line ji in the power distribution network, mu ij is a switching state of the line ij in the power distribution network, N i is a set of parent nodes taking a node i in the power distribution network as a child node, and N 1 is a set of line head nodes in the power distribution network;
Wherein, when node i in the power distribution network is the parent node of node j, τ ij =1, otherwise τ ij =0, and when node j in the power distribution network is the parent node of node i, τ ji =1, otherwise τ ji =0.
In the optimal embodiment of the invention, a nonlinear power flow model is adopted by a traditional optimal scheduling model, and compared with a linear model, the optimal scheduling model is difficult to ensure the globally solved and has larger solving difficulty.
Further, the equivalent current of the active power output provided by the large power grid to node i in the power distribution network at the t-th moment of the control period is determined according to the following formula
Determining an equivalent current of reactive power output provided by the large power grid to node i in the power distribution network at a t-th moment of the control period as follows
Determining equivalent current of active output of distributed power source input by node i in power distribution network at t-th moment of control period according to the following method
Determining equivalent current of reactive power output of distributed power supply input by node i in power distribution network at t-th moment of control period according to the following method
Determining equivalent current of active power output of new energy input by node u in power distribution network at t-th moment of control period according to the following method
Determining equivalent current of reactive power output of new energy input by node i in power distribution network at t moment of control period according to the following method
Determining the equivalent current of the active load demand of node i in the distribution network at the t-th moment of the control period as follows
Determining the equivalent current of reactive load demand of node i in the distribution network at the t-th moment of the control period as follows
The equivalent conductance G i,t of the load demand of the node i in the power distribution network at the t-th moment of the control period is determined as follows:
The equivalent susceptance B i,t of the load demand of node i in the distribution network at the t-th moment of the control period is determined as follows:
Wherein V i 0 is the reference voltage of a node i in the power distribution network, To control the active power provided by the large grid to node i in the distribution network at time t of the cycle,To control the reactive power output provided by the large grid to node i in the distribution network at time t of the cycle,For the active output of the distributed power source input by the node i in the power distribution network at the t-th moment of the control period, Z i,t is the reciprocal of the voltage of the node u in the power distribution network at the t-th moment of the control period,For reactive power output of the distributed power supply input by the node i in the power distribution network at the t-th moment of the control period,To control the active output of the new energy input by the node u in the power distribution network at the t-th moment of the period,For the reactive power output of the new energy input by the node i in the power distribution network at the t-th moment of the control period,For controlling the active load demand of node i in the power distribution network at the t-th moment of the period, C I is the proportionality coefficient between the reactive load demand of node i in the power distribution network and the quadratic term of the voltage of node i in the power distribution network,For controlling the reactive load demand of node i in the power distribution network at the t-th moment of the period, C 'I is the proportionality coefficient between the reactive load demand of node i in the power distribution network and the primary term of the voltage of node i in the power distribution network, C' Z is the proportionality coefficient between the active load demand of node i in the power distribution network and the primary term of the voltage of node i in the power distribution network, C Z is the proportionality coefficient between the active load demand of node i in the power distribution network and the secondary term of the voltage of node i in the power distribution network, Is a preset maximum value of the inverse of the voltage at node u in the distribution network,The preset minimum value is the reciprocal of the voltage of the node u in the power distribution network.
In the preferred embodiment of the present invention, the actual operating conditions for the load are mostly related to the actual voltage level of the node at which it is located, so that the impact of fluctuations in voltage level on node load should be considered in optimizing the scheduling process.
The current load model is a ZI model, a ZIP model, an index model and the like, and the ZI model is adopted in the invention.
Wherein P (V i) and Q (V i) are respectively the actual active demand and reactive demand of the node i, P 0 and Q 0 are respectively the initial active demand and reactive demand of the load node, and V i and V 0 are respectively the actual operating voltage and initial reference voltage of the node; And The current phasor, the load phasor and the voltage phasor of the node i are respectively, and V i re and V i im are respectively the real part and the imaginary part of the voltage of the node i. The load is mainly the required power in the present invention. This current flows out of the node in a direction from the node to the load:
combining the 2 formulas can obtain an equivalent calculation formula of the current:
In a power distribution network, the phase angle of the node voltage is generally small, which means that the imaginary part V i im of the node voltage is small. In practice, the real part of the node voltage tends to be hundreds times more than the imaginary part. Mathematically, this means that the imaginary part of the node voltage is negligible.
Decomposing the current phasor into a real part and an imaginary part, a calculation formula of the current can be obtained,AndRepresenting the real and imaginary parts of the current at node i, respectively.
The decomposition of the current phasor into real and imaginary parts can be expressed in turn as:
and (3) performing equivalent substitution according to the corresponding equality of the real part and the imaginary part, so as to obtain:
wherein P i L and The active and reactive demands of the node i are respectively; And The equivalent branch conductance and the equivalent branch susceptance of the node i are respectively; And The voltage value of the reference node of the node i is V i 0, and is generally set as rated voltage.
The equivalent model of the output injected into the active power distribution network node by the large power grid is as follows:
Wherein P i G and Active and reactive power output injected into a node i of the power distribution network by a large power grid are respectively; And Real and imaginary parts of equivalent current injected into a power distribution network node i by a large power grid are respectively:
Considering that the imaginary part approximation of the node voltage is negligible and that the substation node is selected as a balance node in the distribution network, the approximation of the substation node voltage and the balance node voltage can result in the injection current approximation of the power supply equivalent current source as follows:
The node voltage of the substation may be approximately the same as the voltage of the balance node, and the balance node is typically set to a voltage value of 1.0 pu. Other nodes will not meet this setting, except that the node voltage of the transformer can be approximated.
The following transformation is performed when a distributed power source, a new energy source or other power sources are connected to some nodes, namely:
In the formula, For the real part of the equivalent injection current of the distributed power/new energy output at node i,For the imaginary part of the equivalent injection current of the distributed power supply/new energy output at the node i, P i is the active output of the distributed power supply/new energy input at the node i, Q i is the reactive output of the distributed power supply/new energy input at the node i, Z i is the derivative of the node i voltage,Is the minimum value of the reciprocal of the voltage at node i,Is the maximum value of the reciprocal of the voltage at node i.
The current dispatching can not fully take a large amount of external uncertain factors in the operation of the power distribution network, and the factors can cause potential safety hazards in the operation of the system, so that the power supply capacity of the system is also considered on the premise of considering the maximum new energy utilization rate and the minimum tie switch action of the dispatching.
The power supply capability is mainly represented as the capability of meeting the load, different power supply capability of the distribution network is determined by different load levels, and when the load level is lower, the power supply capability of the distribution network is stronger and the irregularity is smaller. The power supply capacity of the power distribution network is measured, the change of the load is to be planned, and the evaluation results of different load growth modes are different, including but not limited to two growth modes, namely that the current total actual load is increased in proportion, the load of a local area is increased in proportion, and other loads are kept unchanged. In order to ensure the reliability of power supply, the invention designs a checking model, wherein the objective function of the checking model is determined according to the following formula:
further, an objective function of the verification model is determined as follows:
In the formula, t epsilon (1-N T),NT is the total time of the control period, Z t is the variation of the sum of actual injection currents of all nodes in the power distribution network compared with the previous value after the load demand of all nodes in the power distribution network is amplified at the t-th time of the control period;
The variation Z t of the sum of the actual injection currents of all nodes in the power distribution network compared with the previous value is amplified at the t-th moment of the following control period:
In the formula, Amplifying the real part increment of the actual injection current of the node i in the power distribution network after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the real part reduction amount of the actual injection current of the node i in the power distribution network after the load demand of each node in the power distribution network at the t-th moment of the control period,The imaginary part increment of the actual injection current of the node i in the power distribution network after the load demand of each node in the power distribution network is amplified for the t-th moment of the control period,In order to amplify the imaginary part reduction amount of the actual injection current of the node i in the power distribution network after the load demand of each node in the power distribution network at the t-th moment of the control period, i epsilon (1-N), N is the total number of nodes in the power distribution network,AndAre not less than 0;
the interactive output constraint conditions of the active power distribution network and the large power network of the objective function of the inspection model are determined as follows:
In the formula, Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the equivalent current of the active output of the node i by the large power grid,In order to amplify the load demand of each node in the distribution network at the t-th moment of the control period and then provide the equivalent current of the active output of the node i by the large power network with the preset minimum value,Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the equivalent current of the active output of the node i by the large power network to the preset maximum value,Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the reactive power output equivalent current of the node i by the large power network,In order to amplify the load demand of each node in the distribution network at the t-th moment of the control period and then provide the reactive output equivalent current of the node i by the large power network with the preset minimum value,Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the preset maximum value of the equivalent current of the reactive power output of the node i by the large power network;
determining an output linear power flow constraint condition of an objective function of the test model according to the following formula:
Wherein G ij is the conductance of the line ij in the power distribution network, B ij is the susceptance of the line ij in the power distribution network, χ ij,t is the switching state of the line ij in the power distribution network at the t-th moment of the control period in the inspection model, Amplifying the real part of the voltage of the node j after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the real part of the voltage of the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the imaginary part of the node j voltage after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the imaginary part of the i-voltage of the node after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the equivalent current of the active output of the node i by the large power grid,Amplifying the equivalent current of the active output of the distributed power supply input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,In order to verify the equivalent current corresponding to the active power output of the new energy input by the node i in the power distribution network at the t-th moment of the control period in the model,Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the reactive power output equivalent current of the node i by the large power network,In order to amplify the equivalent current of the reactive power of the distributed power supply input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the equivalent current of the reactive power output of the new energy input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the equivalent current of the active load demand of the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,To amplify the equivalent conductance of the load demand of node i after the load demand of each node in the distribution network at the t-th moment of the control period,In order to amplify the equivalent susceptance of the load demand of the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the equivalent current of reactive power demand of a node i after load demand of each node in the power distribution network at the t-th moment of a control period, j epsilon (1-n), n being the total number of other nodes except the node i in the power distribution network, For controlling the optimal switching state of the line ij in the distribution network at the t-th moment of the cycle,The method comprises the steps of controlling the active power output of the optimal new energy source of a node i in the power distribution network at the t-th moment of a period;
determining a branch current constraint of an objective function of the test model according to the following formula:
In the formula, To amplify the real part of the current on the line ij in the distribution network after the load demand of each node in the distribution network at the t-th moment of the control period,To amplify the imaginary part of the current on the line ij in the distribution network after the load demand of each node in the distribution network at the t-th moment of the control period,Maximum current allowed on line ij in the distribution network;
determining a voltage constraint of an objective function of the inspection model according to the following formula:
In the formula, Amplifying the real part of the voltage of a node i after the load demand of each node in the power distribution network at the t-th moment of the control period,To amplify the voltage maximum at node i after the load demand at each node in the distribution network,The voltage minimum value of the node i after the load demand of each node in the power distribution network is amplified;
determining a distributed power constraint condition of an objective function of a test model of the power distribution network according to the following formula:
In the formula, Amplifying the equivalent current of the active output of the distributed power supply input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,To amplify the minimum value of the equivalent current of the active output of the distributed power supply of the input node i after the load demand of each node in the power distribution network,The maximum value of the equivalent current of the active output of the distributed power supply of the node i after the load demand of each node in the power distribution network is amplified;
determining new energy reactive power output constraint conditions of an objective function of an inspection model of the power distribution network according to the following steps:
In the formula, Amplifying the minimum value of the equivalent current of the reactive power output of the new energy input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,And amplifying the maximum value of the equivalent current of the reactive power output of the new energy input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period.
In the specific example of the invention, if the objective function value solved by the inspection model is 0, that is, the inspection model can converge, it is indicated that the system can still stably run under the condition that the load requirement of each node is increased in the same proportion (for example, the load requirement of each node is amplified by K times as the original load requirement), otherwise, the power supply capability cannot meet the requirement of the power supply capability, a group of inspection constraint conditions needs to be added to the optimal operation scheduling model and solved again, and because the fluctuation of the load is in a certain range in the actual power supply process, the fluctuation amount of the load is not more than 30% of the load, so K can be 1.3.
Further, the verification constraint is determined as follows:
wherein lambda ij,t is a formula Mu ij,t is the switching state of the line ij in the distribution network at the t-th moment of the control period, and rho ij,t isIs used for the dual-pair variable of (c),For controlling active power output of new energy of a node i in the power distribution network at the t-th moment of the period, ij epsilon (1-N B),NB is the total number of lines in the power distribution network).
The invention provides a power distribution network optimal scheduling system considering power supply capacity, as shown in fig. 2, the system comprises:
The determining module is used for determining the optimal output of each power supply node in the power distribution network and the optimal switching state of each circuit in the power distribution network according to the load demand of each load node in the power distribution network;
The control module is used for respectively controlling the output of each power supply node in the power distribution network and the switching state of each circuit in the power distribution network to be the optimal output and the optimal switching state;
The power supply nodes are new energy nodes or distributed power supply nodes.
Specifically, the determining module includes:
The first generation unit is used for substituting the load demands of all load nodes in the power distribution network into a pre-built optimal operation scheduling model, solving the pre-built optimal operation scheduling model, and obtaining the optimal output of all power supply nodes in the power distribution network and the optimal switching state of all circuits in the power distribution network;
The second substituting unit is used for amplifying the load demands of all load nodes in the power distribution network, substituting the amplified load demands of all load nodes in the power distribution network, the optimal active power output of all new energy nodes in the power supply nodes in the power distribution network and the optimal switching state of all lines in the power distribution network into a pre-established test model, solving the pre-established test model, and obtaining an objective function value output by the test model;
And (c) judging whether the objective function value output by the test model is 0, if so, outputting the optimal active power output and the optimal reactive power output of each power supply node in the power distribution network and the optimal switching state of each circuit in the power distribution network, otherwise, adding test constraint conditions in the pre-built optimal operation scheduling model, and returning to the step (a).
Further, an objective function of the optimal operation scheduling model is determined according to the following formula:
minF=w1·f1+w2·f2
Wherein F is an objective function value of an optimal operation scheduling model, w 1 is a weight corresponding to the new energy utilization rate, F 1 is the reciprocal of the new energy utilization rate of the power distribution network in a control period, w 2 is a weight corresponding to the line switch action times of the power distribution network, and F 2 is the line switch action times of the power distribution network in the control period;
The new energy utilization rate f 1 of the power distribution network in the control period is determined according to the following formula:
In the formula, For the output of the xth new energy node in the power nodes of the power distribution network consumed at the t-th moment of the control period, P NE,x (t) is the actual output of the xth new energy node in the power nodes of the power distribution network at the t-th moment of the control period (t epsilon (1-N T),NT is the total time of the control period, x epsilon (1-N J),NJ is the total number of new energy nodes in the power nodes of the power distribution network);
the line switch action times f 2 of the power distribution network in the control period are determined according to the following steps:
Wherein mu ij,t is the switching state of the line ij in the power distribution network at the t-th moment of the control period, mu ij,t-1 is the switching state of the line ij in the power distribution network at the t-1 th moment of the control period, ij epsilon (1-N B),NB is the total number of lines in the power distribution network;
If the switching state of the line ij in the power distribution network at the t-th moment of the control period is on, μ ij,t =1, otherwise μ ij,t =0;
If the switching state of the line ij in the power distribution network at the t-1 time of the control period is on, μ ij,t-1 =1, otherwise μ ij,t-1 =0.
Still further, the interactive output constraint condition of the active power distribution network and the large power network of the objective function of the optimal operation scheduling model is determined according to the following formula:
In the formula, To control the equivalent current of the active power output provided by the large grid to node i in the distribution network at time t of the cycle,For a preset minimum value of the equivalent current of the active output provided by the large grid to node i in the distribution network,For a preset maximum value of the equivalent current of the active output provided by the large grid to node i in the distribution network,For controlling the equivalent current of reactive power output provided by the large grid to node i in the distribution network at time t of the period,For a preset minimum value of the equivalent current of reactive power output supplied by the large grid to node i in the distribution network,The method comprises the steps that the preset maximum value of the equivalent current of reactive power output provided by a large power grid to a node i in a power distribution network is set;
Determining the output linear power flow constraint condition of the objective function of the optimal operation scheduling model according to the following steps:
Wherein G ij is the conductance of a line ij in the power distribution network, B ij is the susceptance of the line ij in the power distribution network, beta ij,t is the product of the switching state of the line ij in the power distribution network at the t moment of a control period and the voltage imaginary part of a node j in the power distribution network, beta ji,t is the product of the switching state of the line ij in the power distribution network at the t moment of the control period and the voltage imaginary part of the node i in the power distribution network, alpha ij,t is the product of the switching state of the line ij in the power distribution network at the t moment of the control period and the voltage real part of the node j in the power distribution network, alpha ji,t is the product of the switching state of the line ij in the power distribution network at the t moment of the control period and the voltage real part of the node i in the power distribution network, To control the equivalent current of the active power output provided by the large grid to node i in the distribution network at time t of the cycle,For controlling the equivalent current of the active output of the distributed power supply input by the node i in the power distribution network at the t-th moment of the period,To control the real part of the voltage at node i in the distribution network at time t of the cycle,For controlling the equivalent current of the active power output of the new energy input by the node i in the power distribution network at the t-th moment of the period,To control the voltage imaginary part of node i in the distribution network at time t of the cycle,For controlling the equivalent current of reactive power output provided by the large grid to node i in the distribution network at time t of the period,For controlling the equivalent current of reactive power output of the distributed power supply input by the node i in the power distribution network at the t-th moment of the period,For controlling the equivalent current of the reactive power output of the new energy input by the node i in the power distribution network at the t-th moment of the period,For the equivalent current of the active load demand of node i in the power distribution network at the t-th moment of the control period, G i,t is the equivalent conductance of the load demand of node i in the power distribution network at the t-th moment of the control period, B i,t is the equivalent susceptance of the load demand of node i in the power distribution network at the t-th moment of the control period,For the equivalent current required by reactive load of a node i in the power distribution network at the t-th moment of the control period, j epsilon (1-n), wherein n is the total number of other nodes except the node i in the power distribution network;
Wherein, the To control the real part of the voltage at node j in the distribution network at time t of the cycle,A maximum value is preset for the real part of the voltage at node j in the distribution network,A minimum value is preset for the real part of the voltage at node j in the distribution network,To control the voltage imaginary part of node j in the distribution network at time t of the cycle,A maximum value is preset for the voltage imaginary part of the node j in the distribution network,Presetting a minimum value for the voltage imaginary part of a node j in the power distribution network;
determining branch current constraint conditions of an objective function of an optimal operation scheduling model according to the following steps:
In the formula, To control the real current of the line ij in the distribution network at the t-th moment of the cycle,To control the imaginary current of the line ij in the distribution network at the t-th moment of the period,Maximum current allowed on line ij in the distribution network;
Determining a voltage constraint condition of an objective function of an optimal operation scheduling model of the power distribution network according to the following steps:
In the formula, For the real part of the voltage of the node i in the power distribution network at the t-th moment of the control period, V i max is the maximum value of the voltage of the node i in the power distribution network, and V i min is the minimum value of the voltage of the node i in the power distribution network;
Determining a distributed power constraint condition of an objective function of an optimal operation scheduling model according to the following steps:
In the formula, For controlling the equivalent current of the active output of the distributed power supply input by the node i in the power distribution network at the t-th moment of the period,The minimum value of the equivalent current of the active output of the distributed power source input to the node i in the power distribution network,The maximum value of equivalent current of active output of the distributed power supply input to a node i in the power distribution network;
determining a network radial constraint condition of an objective function of the optimal operation scheduling model according to the following steps:
Wherein τ ij is an output flow direction state variable of a line ij in the power distribution network, τ ji is an output flow direction state variable of a line ji in the power distribution network, mu ij is a switching state of the line ij in the power distribution network, N i is a set of parent nodes taking a node i in the power distribution network as a child node, and N 1 is a set of line head nodes in the power distribution network;
Wherein, when node i in the power distribution network is the parent node of node j, τ ij =1, otherwise τ ij =0, and when node j in the power distribution network is the parent node of node i, τ ji =1, otherwise τ ji =0.
Still further, an equivalent current of the active power provided by the large grid to node i in the distribution network at time t of the control period is determined as follows
Determining an equivalent current of reactive power output provided by the large power grid to node i in the power distribution network at a t-th moment of the control period as follows
Determining equivalent current of active output of distributed power source input by node i in power distribution network at t-th moment of control period according to the following method
Determining equivalent current of reactive power output of distributed power supply input by node i in power distribution network at t-th moment of control period according to the following method
Determining equivalent current of active power output of new energy input by node u in power distribution network at t-th moment of control period according to the following method
Determining equivalent current of reactive power output of new energy input by node i in power distribution network at t moment of control period according to the following method
Determining the equivalent current of the active load demand of node i in the distribution network at the t-th moment of the control period as follows
Determining the equivalent current of reactive load demand of node i in the distribution network at the t-th moment of the control period as follows
The equivalent conductance G i,t of the load demand of the node i in the power distribution network at the t-th moment of the control period is determined as follows:
determining an equivalent susceptance B of the load demand of a node i in the distribution network at the t-th moment of the control period as follows i,t:
Wherein V i 0 is the reference voltage of a node i in the power distribution network,To control the active power provided by the large grid to node i in the distribution network at time t of the cycle,To control the reactive power output provided by the large grid to node i in the distribution network at time t of the cycle,For the active output of the distributed power source input by the node i in the power distribution network at the t-th moment of the control period, Z i,t is the reciprocal of the voltage of the node u in the power distribution network at the t-th moment of the control period,For reactive power output of the distributed power supply input by the node i in the power distribution network at the t-th moment of the control period,To control the active output of the new energy input by the node u in the power distribution network at the t-th moment of the period,For the reactive power output of the new energy input by the node i in the power distribution network at the t-th moment of the control period,For controlling the active load demand of node i in the power distribution network at the t-th moment of the period, C I is the proportionality coefficient between the reactive load demand of node i in the power distribution network and the quadratic term of the voltage of node i in the power distribution network,For controlling the reactive load demand of node i in the power distribution network at the t-th moment of the period, C 'I is the proportionality coefficient between the reactive load demand of node i in the power distribution network and the primary term of the voltage of node i in the power distribution network, C' Z is the proportionality coefficient between the active load demand of node i in the power distribution network and the primary term of the voltage of node i in the power distribution network, C Z is the proportionality coefficient between the active load demand of node i in the power distribution network and the secondary term of the voltage of node i in the power distribution network, Is a preset maximum value of the inverse of the voltage at node u in the distribution network,The preset minimum value is the reciprocal of the voltage of the node u in the power distribution network.
Further, an objective function of the verification model is determined as follows:
In the formula, t epsilon (1-N T),NT is the total time of the control period, Z t is the variation of the sum of actual injection currents of all nodes in the power distribution network compared with the previous value after the load demand of all nodes in the power distribution network is amplified at the t-th time of the control period;
The variation Z t of the sum of the actual injection currents of all nodes in the power distribution network compared with the previous value is amplified at the t-th moment of the following control period:
In the formula, Amplifying the real part increment of the actual injection current of the node i in the power distribution network after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the real part reduction amount of the actual injection current of the node i in the power distribution network after the load demand of each node in the power distribution network at the t-th moment of the control period,The imaginary part increment of the actual injection current of the node i in the power distribution network after the load demand of each node in the power distribution network is amplified for the t-th moment of the control period,In order to amplify the imaginary part reduction amount of the actual injection current of the node i in the power distribution network after the load demand of each node in the power distribution network at the t-th moment of the control period, i epsilon (1-N), N is the total number of nodes in the power distribution network,AndAre not less than 0;
the interactive output constraint conditions of the active power distribution network and the large power network of the objective function of the inspection model are determined as follows:
In the formula, Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the equivalent current of the active output of the node i by the large power grid,In order to amplify the load demand of each node in the distribution network at the t-th moment of the control period and then provide the equivalent current of the active output of the node i by the large power network with the preset minimum value,Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the equivalent current of the active output of the node i by the large power network to the preset maximum value,Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the reactive power output equivalent current of the node i by the large power network,In order to amplify the load demand of each node in the distribution network at the t-th moment of the control period and then provide the reactive output equivalent current of the node i by the large power network with the preset minimum value,Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the preset maximum value of the equivalent current of the reactive power output of the node i by the large power network;
determining an output linear power flow constraint condition of an objective function of the test model according to the following formula:
Wherein G ij is the conductance of the line ij in the power distribution network, B ij is the susceptance of the line ij in the power distribution network, χ ij,t is the switching state of the line ij in the power distribution network at the t-th moment of the control period in the inspection model, Amplifying the real part of the voltage of the node j after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the real part of the voltage of the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the imaginary part of the node j voltage after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the imaginary part of the i-voltage of the node after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the equivalent current of the active output of the node i by the large power grid,Amplifying the equivalent current of the active output of the distributed power supply input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,In order to verify the equivalent current corresponding to the active power output of the new energy input by the node i in the power distribution network at the t-th moment of the control period in the model,Amplifying the load demand of each node in the power distribution network at the t-th moment of the control period, and providing the reactive power output equivalent current of the node i by the large power network,In order to amplify the equivalent current of the reactive power of the distributed power supply input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the equivalent current of the reactive power output of the new energy input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the equivalent current of the active load demand of the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,To amplify the equivalent conductance of the load demand of node i after the load demand of each node in the distribution network at the t-th moment of the control period,In order to amplify the equivalent susceptance of the load demand of the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,Amplifying the equivalent current of reactive power demand of a node i after load demand of each node in the power distribution network at the t-th moment of a control period, j epsilon (1-n), n being the total number of other nodes except the node i in the power distribution network, For controlling the optimal switching state of the line ij in the distribution network at the t-th moment of the cycle,The method comprises the steps of controlling the active power output of the optimal new energy source of a node i in the power distribution network at the t-th moment of a period;
determining a branch current constraint of an objective function of the test model according to the following formula:
In the formula, To amplify the real part of the current on the line ij in the distribution network after the load demand of each node in the distribution network at the t-th moment of the control period,To amplify the imaginary part of the current on the line ij in the distribution network after the load demand of each node in the distribution network at the t-th moment of the control period,Maximum current allowed on line ij in the distribution network;
determining a voltage constraint of an objective function of the inspection model according to the following formula:
In the formula, Amplifying the real part of the voltage of a node i after the load demand of each node in the power distribution network at the t-th moment of the control period,To amplify the voltage maximum at node i after the load demand at each node in the distribution network,The voltage minimum value of the node i after the load demand of each node in the power distribution network is amplified;
determining a distributed power constraint condition of an objective function of a test model of the power distribution network according to the following formula:
In the formula, Amplifying the equivalent current of the active output of the distributed power supply input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,To amplify the minimum value of the equivalent current of the active output of the distributed power supply of the input node i after the load demand of each node in the power distribution network,The maximum value of the equivalent current of the active output of the distributed power supply of the node i after the load demand of each node in the power distribution network is amplified;
determining new energy reactive power output constraint conditions of an objective function of an inspection model of the power distribution network according to the following steps:
In the formula, Amplifying the minimum value of the equivalent current of the reactive power output of the new energy input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period,And amplifying the maximum value of the equivalent current of the reactive power output of the new energy input by the node i after the load demand of each node in the power distribution network at the t-th moment of the control period.
Still further, the verification constraint is determined as follows:
wherein lambda ij,t is a formula Mu ij,t is the switching state of the line ij in the distribution network at the t-th moment of the control period, and rho ij,t isIs used for the dual-pair variable of (c),For controlling active power output of new energy of a node i in the power distribution network at the t-th moment of the period, ij epsilon (1-N B),NB is the total number of lines in the power distribution network).
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that the above embodiments are only for illustrating the technical solution of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the above embodiments, it should be understood by those skilled in the art that modifications and equivalents may be made to the specific embodiments of the present invention without departing from the spirit and scope of the present invention, and any modifications and equivalents are intended to be included in the scope of the claims of the present invention.