CN109524955A - A kind of active distribution network optimization method of consideration source lotus voltage characteristic - Google Patents
A kind of active distribution network optimization method of consideration source lotus voltage characteristic Download PDFInfo
- Publication number
- CN109524955A CN109524955A CN201811196575.XA CN201811196575A CN109524955A CN 109524955 A CN109524955 A CN 109524955A CN 201811196575 A CN201811196575 A CN 201811196575A CN 109524955 A CN109524955 A CN 109524955A
- Authority
- CN
- China
- Prior art keywords
- power
- node
- reactive
- active
- distribution network
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000005457 optimization Methods 0.000 title claims abstract description 50
- 240000002853 Nelumbo nucifera Species 0.000 title abstract 6
- 235000006508 Nelumbo nucifera Nutrition 0.000 title abstract 6
- 235000006510 Nelumbo pentapetala Nutrition 0.000 title abstract 6
- 230000005540 biological transmission Effects 0.000 claims abstract description 42
- 230000001360 synchronised effect Effects 0.000 claims description 57
- 238000010248 power generation Methods 0.000 claims description 30
- 230000005611 electricity Effects 0.000 claims description 13
- 230000005284 excitation Effects 0.000 claims description 12
- 150000001875 compounds Chemical class 0.000 claims description 3
- 230000001939 inductive effect Effects 0.000 claims description 3
- 230000008092 positive effect Effects 0.000 abstract description 2
- 238000005553 drilling Methods 0.000 abstract 1
- 230000003068 static effect Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009429 electrical wiring Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/008—Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The present invention relates to a kind of active distribution network optimization methods of consideration source lotus voltage characteristic, belong to power system active power distribution network technical field.With reactive compensation capacity, distributed generation resource active and reactive power is decision content, minimize active distribution network purchases strategies, in the case where meeting power distribution network safe operation technical requirements, the active distribution network optimal operation model of consideration source lotus voltage characteristic is constructed, and selects CONOPT solver to solve mentioned Optimized model based on GAMS Optimization Platform.The invention proposes the active distribution network optimal operation models for considering source lotus voltage characteristic, digging system itself source lotus voltage characteristic potential quality to a greater degree, avoid the confinement problems of conventional electrical distribution net running optimizatin decision, illustrate that consideration source lotus voltage characteristic more meets distribution system practical operation situation, and the invention exchanges power with positive effect with higher level's power transmission network to reduction distribution system operation totle drilling cost, reduction, is conducive to promote distribution system on-road efficiency.
Description
Technical Field
The invention relates to an active power distribution network optimization method considering source charge voltage characteristics, and belongs to the technical field of active power distribution networks of power systems.
Background
The voltage level is the main indicator of the electrical support performance of the grid given the power generation and load power balance. The power system is an artificial closed-loop control power system, and the power generation and load power in the power system is a function taking the voltage level as a parameter, so that the optimal operation of the power grid is actually a source-load power balance mode within a safe voltage level range meeting the technical requirements of normal operation of the power grid under an expected target.
The voltage characteristics of the load actually influence the computational load on the nodes, thereby also changing the power flow distribution in the electrical network to some extent. The static voltage characteristic of the load is considered in the process of optimizing and reconstructing the power distribution network, the static voltage model in the load power function form is merged into the power distribution network reconstruction model, and the network reconstruction difference between the load in the static voltage characteristic form and the load in the constant power form is contrastively analyzed, so that the static voltage model is more respectful to reality and can reflect the real power flow distribution, node voltage change and system network loss conditions. The voltage regulation effect of source load power is considered in the economic dispatching of the power system, and the effect can promote the power generation consumption of renewable energy sources and improve the economical efficiency of the operation of a power grid, but the method is mainly researched on the level of a power transmission network, the processing of the voltage characteristic of the source load is only limited to the voltage static characteristic of a synchronous generator set with difference voltage regulation characteristic and load, and the voltage characteristics of wind power, photovoltaic and the like are not discussed. Therefore, the influence of the voltage response characteristic on the decision result is analyzed in the optimization decision, so that the promotion of the power grid operation level improvement is considered to a greater extent, and the method is worthy of further research.
Disclosure of Invention
The invention aims to solve the technical problem of providing an active power distribution network optimization method considering source charge voltage characteristics, and reducing the dependence on voltage regulation and frequency modulation resources of a superior transmission network by calling active flexible resources in a power distribution network.
The technical scheme adopted by the invention is as follows: an active power distribution network optimization method considering source charge voltage characteristics comprises the following steps:
step 1: the reactive compensation capacity, the active power and the reactive power of the distributed power supply are taken as decision quantities, and the minimum electricity purchasing cost of the active power distribution network is taken as a target;
step 2: the method comprises the steps of taking electrical physical constraints and safety technical requirements which need to be met by the operation of an active power distribution system as constraint conditions, and constructing an active power distribution network operation optimization model considering source charge voltage characteristics;
and step 3: and (3) solving the optimization model in the step (2) by using a CONOPT solver based on a GAMS optimization platform.
Specifically, the minimum objective function of the electricity purchasing cost of the active power distribution network in the step 1 is as follows:
in the formula, JGThe method is a set formed by all schedulable synchronous units (small hydroelectric generating sets and gas turbines) in the active power distribution network; pex,0Active power from a superordinate transmission network for an active distribution network, Cex,0Representing the marginal price of electricity of the root node; pgThe active power output by the schedulable synchronous unit g; cg() The method is a generating cost characteristic function of the schedulable synchronous unit g.
Specifically, in the step 2, the electrical physical constraint and the safety technical requirement which need to be met by the operation of the power distribution system are taken as constraint conditions, and the constraint conditions are divided into equality constraint and inequality constraint, which are specifically as follows:
① equality constraint
(1) Node active power and reactive power balance constraints (power flow equation):
wherein, PlAnd QlThe active power and the reactive power transmitted on the power transmission element l can be represented as formula (4) and formula (5); j. the design is a squareS,iThe node i is a set formed by all power transmission elements with the node i as a first node; j. the design is a squareE,iTo take node i as the end nodeA set of all power transmission elements of a point; j. the design is a squareG,iThe method is a set formed by all schedulable synchronous units (small hydroelectric generating sets and gas turbines) on a node i; j. the design is a squareW,iThe method comprises the steps that a set is formed by all doubly-fed wind power on a node i; j. the design is a squareV,iThe photovoltaic power generation system is a set formed by all photovoltaic power generation systems on a node i; j. the design is a squareD,iIs a set formed by all power loads on the node i; j. the design is a squareC,iThe method comprises the following steps of (1) forming a set by all reactive compensation equipment on a node i; j. the design is a squareMA set formed by all nodes in the power distribution system; pgAnd QgRespectively the active power and the reactive power output by the schedulable synchronous unit g; pwAnd QwRespectively the active power and the reactive power output by the double-fed wind power w; pvAnd QvRespectively the active power and the reactive power output by the photovoltaic power generation system v; qCThe inductive reactive power output by the reactive compensation equipment c; pdAnd QdThe active power and the reactive power of the power load d are respectively required;
in the formula, thetaijThe phase angle difference of voltage phasors of the node i and the node j is represented; vliRepresenting the voltage amplitude of a first node i node of a power transmission element l; vljRepresenting the voltage amplitude of the node j at the tail node of the power transmission element l; glAnd blRespectively the conductance value and the susceptance value of the power transmission branch l;
(2) related equation constraint of schedulable synchronous units (small hydroelectric generating units and gas turbines):
wherein,the g stator reactive current of the schedulable synchronous unit is obtained; viRepresenting the voltage amplitude of a node i where the synchronous unit g is located;indicating the no-load voltage set by the excitation control of the synchronous unit g; kgRepresenting the voltage difference adjustment coefficient of the synchronous unit g; egRepresenting the internal potential of the synchronous unit g; deltagRepresenting the power angle value of the synchronous unit g;a direct-axis reactor of a synchronous unit g; j. the design is a squareGThe method is a set formed by all schedulable synchronous units (small hydroelectric generating units and gas turbines);
(3) load voltage characteristics:
in the formula,andrespectively the active power and the reactive power of the power load d under the rated voltage levelA power requirement;andthe constant impedance active power and the reactive power of the power load d are respectively;andthe constant current active power and the reactive power of the power load d are respectively part;andrespectively are a constant power active power part and a reactive power part of the power load d; v0Rated voltage of the system; j. the design is a squareDA set of all electrical loads;
② inequality constraint
(1) Output power range constraint of schedulable synchronous units (small hydroelectric generating units and gas turbines):
in the formula,andrespectively an upper limit and a lower limit of active power allowed by the synchronous unit g;andrespectively an upper limit and a lower limit of an excitation potential of the synchronous unit g;the maximum value of the g stator current of the synchronous unit;
(2) and (3) restricting the operation range of the doubly-fed wind turbine generator:
in the formula, JWThe method comprises the steps of (1) forming a set for all double-fed wind power; vwThe voltage is the stator side machine end voltage of the doubly-fed wind turbine generator w;representing the maximum active power which can be output by the doubly-fed wind turbine generator w under the limitation of natural conditions;the reactance is the stator reactance of the doubly-fed wind turbine generator w;the excitation reactance is the excitation reactance of the double-fed wind turbine generator w;the maximum current of the w rotor side of the doubly-fed wind turbine generator set is obtained; j. the design is a squareWThe method comprises the steps of (1) forming a set by all the double-fed wind power;
(3) and (3) restricting the operation range of the photovoltaic power generation system:
wherein:representing the maximum active power that the photovoltaic power generation system v can output due to the limitation of natural conditions; vvRepresenting the node voltage of the grid-connected side of the photovoltaic power generation system v;representing the maximum load current of the photovoltaic power generation system v inverter; j. the design is a squareVThe photovoltaic power generation system is a set formed by all photovoltaic power generation systems;
(4) and (3) limiting the upper limit and the lower limit of the node voltage amplitude:
in the formula, Vi maxAnd Vi minRespectively representing the upper limit and the lower limit of the voltage amplitude of the node i;
(5) the allowable thermal current range constraint of the power transmission element:
in the formula (19), the compound represented by the formula (I),represents the maximum current of the transmission element l; i isl,ijThe current amplitude of the power transmission element/is represented as:
|Il,ij|＝|Yl|*[(Vli)2+(Vlj)2-2VliVljcosθij]1/2(20)
in the formula (20), YlRepresents the admittance modulus value of the power transmission element l;
(6) and (3) limiting the upper and lower limits of the capacity of the reactive compensation equipment:
in the formula (21), QcRepresenting the reactive power compensated by the reactive compensation device c;andrespectively obtaining maximum compensation reactive power and minimum compensation reactive power of the reactive compensation equipment c; j. the design is a squareCThe method is a set formed by all reactive compensation equipment of the power distribution system.
The invention has the beneficial effects that:
(1) aiming at the defects of the traditional power distribution network optimization decision-making method, the invention takes the electrical physical constraint and the safety technical requirement which are required to be met by the operation of a power distribution system as constraint conditions, pursues the minimum electricity purchasing cost of the power distribution network, takes the reactive compensation capacity, the active power and the reactive power of the distributed power supply as decision-making quantities, fully considers the influence of the source charge voltage characteristic on a power balance mode and power flow distribution, constructs an active power distribution network operation optimization model considering the source charge voltage characteristic, and realizes the advanced optimization decision-making of the active power balance mode and the reactive voltage support mode of the power distribution system;
(2) the regulation characteristic of the source charge voltage is considered to be more respected to the reality, so that the power demand of the power equipment on the node is not in a constant power form in the decision, but shows a positive correlation change rule along with the voltage change, thereby increasing the feasible domain of the optimization decision to a certain extent and overcoming the limitation and the conservation of the traditional optimization decision;
(3) under the guidance of the goal of minimizing the electricity purchasing cost of the whole power distribution system, the active power distribution network utilizes the power supply of the system as much as possible to reduce the electric quantity demand on the interconnected superior power transmission systems, and the aim of maximizing the utilization of renewable energy sources in the power distribution network is to generate electricity is pursued.
Drawings
FIG. 1 is an overall flow chart of the present invention;
fig. 2 is an electrical wiring diagram of the 41-node power distribution system of the present invention.
Detailed Description
For the purpose of illustrating the objects, technical solutions and advantages of the present invention, the present invention will be further described in detail with reference to the accompanying drawings and examples.
Example 1: an active power distribution network optimization method considering source charge voltage characteristics comprises the following steps:
step 1: the method comprises the following steps of taking reactive compensation capacity, active power and reactive power of a distributed power supply as decision quantities and minimizing the electricity purchasing cost of an active power distribution network as a target;
step 2: the method comprises the steps of taking electrical physical constraints and safety technical requirements which need to be met by the operation of an active power distribution system as constraint conditions, and constructing an active power distribution network operation optimization model considering source charge voltage characteristics;
and step 3: and (3) solving the optimization model in the step (2) by using a CONOPT solver based on a GAMS optimization platform.
Specifically, the active power distribution network optimization model considering the source charge voltage characteristic takes reactive compensation capacity, active power and reactive power of a distributed power supply as decision variables, and the minimum electricity purchasing cost of a power distribution system is an optimization target, and the minimum objective function of the electricity purchasing cost of the active power distribution network in the step 1 is as follows:
in the formula, JGThe method is a set formed by all schedulable synchronous units (small hydroelectric generating sets and gas turbines) in the active power distribution network; pex,0Active power from a superordinate transmission network for an active distribution network, Cex,0Representing the marginal price of electricity of the root node; pgThe active power output by the schedulable synchronous unit g; cg() The method is a generating cost characteristic function of the schedulable synchronous unit g.
Furthermore, the optimization model fully considers the influence of the source charge voltage characteristic on a source charge power balance mode and a power flow distribution mode in the power distribution network, realizes advanced optimization decision on an active power balance mode and a reactive voltage support mode of the power distribution system, and utilizes the power supply of the system as much as possible to reduce the electric quantity demand on an interconnected superior power transmission system.
Specifically, the equality constraint and the inequality constraint in step 2 are specifically as follows:
① equality constraint
(1) Node active power and reactive power balance constraints (power flow equation):
wherein, PlAnd QlThe active power and the reactive power transmitted on the power transmission element l can be represented as formula (4) and formula (5); j. the design is a squareS,iThe node i is a set formed by all power transmission elements with the node i as a first node; j. the design is a squareE,iThe node i is a set formed by all power transmission elements with the node i as a final node; j. the design is a squareG,iThe method is a set formed by all schedulable synchronous units (small hydroelectric generating sets and gas turbines) on a node i; j. the design is a squareW,iThe method comprises the steps that a set is formed by all doubly-fed wind power on a node i; j. the design is a squareV,iThe photovoltaic power generation system is a set formed by all photovoltaic power generation systems on a node i; j. the design is a squareD,iIs a set formed by all power loads on the node i; j. the design is a squareC,iThe method comprises the following steps of (1) forming a set by all reactive compensation equipment on a node i; j. the design is a squareMA set formed by all nodes in the power distribution system; pgAnd QgRespectively the active power and the reactive power output by the schedulable synchronous unit g; pwAnd QwRespectively the active power and the reactive power output by the double-fed wind power w; pvAnd QvRespectively the active power and the reactive power output by the photovoltaic power generation system v; qCThe inductive reactive power output by the reactive compensation equipment c; pdAnd QdThe active power and the reactive power of the power load d are respectively required;
in the formula, thetaijThe phase angle difference of voltage phasors of the node i and the node j is represented; vliRepresenting the voltage amplitude of a first node i node of a power transmission element l; vljRepresenting the voltage amplitude of the node j at the tail node of the power transmission element l; glAnd blRespectively the conductance value and the susceptance value of the power transmission branch l;
(2) related equation constraint of schedulable synchronous units (small hydroelectric generating units and gas turbines):
wherein,the g stator reactive current of the schedulable synchronous unit is obtained; viRepresenting the voltage amplitude of a node i where the synchronous unit g is located;indicating the no-load voltage set by the excitation control of the synchronous unit g; kgRepresenting the voltage difference adjustment coefficient of the synchronous unit g; egRepresenting the internal potential of the synchronous unit g; deltagRepresenting the power angle value of the synchronous unit g;a direct-axis reactor of a synchronous unit g; j. the design is a squareGThe method is a set formed by all schedulable synchronous units (small hydroelectric generating units and gas turbines);
(3) load voltage characteristics:
in the formula,andrespectively the active power and reactive power requirements of the power load d under the rated voltage level;andthe constant impedance active power and the reactive power of the power load d are respectively;andthe constant current active power and the reactive power of the power load d are respectively part;andrespectively are a constant power active power part and a reactive power part of the power load d; v0Rated voltage of the system; j. the design is a squareDA set of all electrical loads;
② inequality constraint
(1) Output power range constraint of schedulable synchronous units (small hydroelectric generating units and gas turbines):
in the formula,andrespectively an upper limit and a lower limit of active power allowed by the synchronous unit g;andrespectively an upper limit and a lower limit of an excitation potential of the synchronous unit g;the maximum value of the g stator current of the synchronous unit;
(2) and (3) restricting the operation range of the doubly-fed wind turbine generator:
in the formula, JWThe method comprises the steps of (1) forming a set for all double-fed wind power; vwThe voltage is the stator side machine end voltage of the doubly-fed wind turbine generator w;representing the maximum active power which can be output by the doubly-fed wind turbine generator w under the limitation of natural conditions;the reactance is the stator reactance of the doubly-fed wind turbine generator w;the excitation reactance is the excitation reactance of the double-fed wind turbine generator w;the maximum current of the w rotor side of the doubly-fed wind turbine generator set is obtained; j. the design is a squareWThe method comprises the steps of (1) forming a set by all the double-fed wind power;
(3) and (3) restricting the operation range of the photovoltaic power generation system:
wherein:representing the maximum active power that the photovoltaic power generation system v can output due to the limitation of natural conditions; vvRepresenting the node voltage of the grid-connected side of the photovoltaic power generation system v;representing the maximum load current of the photovoltaic power generation system v inverter; j. the design is a squareVThe photovoltaic power generation system is a set formed by all photovoltaic power generation systems;
(4) and (3) limiting the upper limit and the lower limit of the node voltage amplitude:
in the formula, Vi maxAnd Vi minRespectively representing the upper limit and the lower limit of the voltage amplitude of the node i;
(5) the allowable thermal current range constraint of the power transmission element:
in the formula (19), the compound represented by the formula (I),represents the maximum current of the transmission element l; i isl,ijThe current amplitude of the power transmission element/is represented as:
|Il,ij|＝|Yl|*[(Vli)2+(Vlj)2-2VliVljcosθij]1/2(20)
in the formula (20), YlRepresents the admittance modulus value of the power transmission element l;
(6) and (3) limiting the upper and lower limits of the capacity of the reactive compensation equipment:
in the formula (21), QcRepresenting the reactive power compensated by the reactive compensation device c;andrespectively obtaining maximum compensation reactive power and minimum compensation reactive power of the reactive compensation equipment c; j. the design is a squareCThe method is a set formed by all reactive compensation equipment of the power distribution system.
The present invention will be further described with reference to the following specific embodiments.
The invention uses a 41-node radial distribution system as an example to verify the effectiveness of the invention. An electrical wiring diagram for a 41-node radial power distribution system is shown in figure 2. In example simulation analysis, the voltage level of the power distribution system is 27.6kV, the feeder capacity is 14.3MVA, 4 doubly-fed wind turbine generators, 2 doubly-fed wind turbine generators and 5 doubly-fed wind turbine generators are respectively arranged on a node 19, a node 28 and a node 40 in the power distribution system, a small hydroelectric turbine generator is arranged on the node 4, a small gas turbine generator is arranged on the node 9 and a small gas turbine generator is arranged on the node 39, and the node 1 is a boundary node connected with a superior power transmission network. The power reference of the power distribution system is selected to be 10MVA, the allowable node voltage fluctuation range is set to be 1 +/-6%, and the optimization period is selected to be 15 min. Through simulation calculation, the operation optimization results of the synchronous generator set, the reactive power compensation device and the wind turbine generator set are respectively shown in tables 1 to 3, and the target function values, the exchange power from the superior transmission network and the network loss value are shown in table 4.
TABLE 1 optimization results of synchronous unit operation
TABLE 2 reactive power compensator operation optimization results
TABLE 3 doubly-fed wind turbine generator system operation optimization results
Table 4 run optimization comparison results
Compared with the power distribution network optimization method without the power distribution network optimization method, the power distribution network optimization method has the advantages that the total optimization operation cost of decision is low, and the exchange power of the power distribution network optimization operation cost and the exchange power of a superior power transmission network are low. The influence of the reactive-voltage characteristics of the wind turbine generator on the operation optimization result of the power distribution network is further analyzed, and the optimization decision result is shown in a table 5.
TABLE 5 run optimization comparison results
It can be known from the comparison of table 5 that, considering the reactive-voltage characteristic of the wind turbine generator set in the operation optimization of the power distribution network has positive effects on reducing the total operation cost of the power distribution system and reducing the exchange power with the superior power transmission network, considering the reactive-voltage characteristic of the wind turbine generator set can better respect the actual operation of the power distribution network, and is beneficial to improving the operation benefit of the power distribution system.
The invention provides an active power distribution network operation optimization model considering source charge voltage characteristics, which shows that the active power distribution network operation optimization model considering the source charge voltage characteristics is more practical, the feasible region of optimization decision can be enlarged to a certain extent, the limitation and the conservation of the traditional optimization decision are overcome, and the advanced optimization decision-making of an active power balance mode and a reactive voltage support mode of a power distribution system is realized
While the present invention has been described in detail with reference to the embodiments shown in the drawings, the present invention is not limited to the embodiments, and various changes can be made without departing from the spirit and scope of the present invention.
Claims (3)
1. An active power distribution network operation optimization method considering source charge voltage characteristics is characterized by comprising the following steps: the method comprises the following steps:
step 1: the reactive compensation capacity and the active and reactive power of the distributed power supply are taken as decision quantities, so that the electricity purchasing cost of the active power distribution network is minimum;
step 2: the method comprises the steps of taking electrical physical constraints and safety technical requirements which need to be met by the operation of an active power distribution system as constraint conditions, and constructing an active power distribution network operation optimization model considering source charge voltage characteristics;
and step 3: and (3) solving the optimization model in the step (2) by using a CONOPT solver based on a GAMS optimization platform.
2. The active power distribution network optimization method considering source charge voltage characteristics according to claim 1, wherein the method comprises the following steps: in the step 1, the objective function for minimizing the electricity purchasing cost of the active power distribution network is as follows:
in the formula, JGThe method is a set formed by all schedulable synchronous units (small hydroelectric generating sets and gas turbines) in the active power distribution network; pex,0Active power from a superordinate transmission network for an active distribution network, Cex,0Representing the marginal price of electricity of the root node; pgThe active power output by the schedulable synchronous unit g; cg() The method is a generating cost characteristic function of the schedulable synchronous unit g.
3. The active power distribution network optimization method considering source charge voltage characteristics according to claim 1, wherein the method comprises the following steps: step 2, the electrical physical constraint and the safety technical requirement which are required to be met by the operation of the power distribution system are taken as constraint conditions, and the constraint conditions are divided into equality constraint and inequality constraint, and the method specifically comprises the following steps:
① equality constraint
(1) Node active power and reactive power balance constraint:
wherein, PlAnd QlThe active power and the reactive power transmitted on the power transmission element l can be represented as formula (4) and formula (5); j. the design is a squareS,iThe node i is a set formed by all power transmission elements with the node i as a first node; j. the design is a squareE,iThe node i is a set formed by all power transmission elements with the node i as a final node; j. the design is a squareG,iThe method is a set formed by all schedulable synchronous units (small hydroelectric generating sets and gas turbines) on a node i; j. the design is a squareW,iThe method comprises the steps that a set is formed by all doubly-fed wind power on a node i; j. the design is a squareV,iThe photovoltaic power generation system is a set formed by all photovoltaic power generation systems on a node i; j. the design is a squareD,iIs a set formed by all power loads on the node i; j. the design is a squareC,iThe method comprises the following steps of (1) forming a set by all reactive compensation equipment on a node i; j. the design is a squareMA set formed by all nodes in the power distribution system; pgAnd QgRespectively the active power and the reactive power output by the schedulable synchronous unit g; pwAnd QwRespectively the active power and the reactive power output by the double-fed wind power w; pvAnd QvRespectively the active power and the reactive power output by the photovoltaic power generation system v; qCThe inductive reactive power output by the reactive compensation equipment c; pdAnd QdThe active power and the reactive power of the power load d are respectively required;
in the formula, thetaijThe phase angle difference of voltage phasors of the node i and the node j is represented; vliRepresenting the voltage amplitude of a first node i node of a power transmission element l; vljRepresenting the voltage amplitude of the node j at the tail node of the power transmission element l; glAnd blRespectively the conductance value and the susceptance value of the power transmission branch l;
(2) the related equation constraint of the schedulable synchronous unit is as follows:
wherein,the g stator reactive current of the schedulable synchronous unit is obtained; viRepresenting the voltage amplitude of a node i where the synchronous unit g is located;indicating the no-load voltage set by the excitation control of the synchronous unit g; kgRepresenting the voltage difference adjustment coefficient of the synchronous unit g; egRepresenting the internal potential of the synchronous unit g; deltagRepresenting the power angle value of the synchronous unit g;a direct-axis reactor of a synchronous unit g; j. the design is a squareGThe method is a set formed by all schedulable synchronous units (small hydroelectric generating units and gas turbines);
(3) load voltage characteristics:
in the formula,andrespectively the active power and reactive power requirements of the power load d under the rated voltage level;andthe constant impedance active power and the reactive power of the power load d are respectively;andthe constant current active power and the reactive power of the power load d are respectively part;andrespectively are a constant power active power part and a reactive power part of the power load d; v0Rated voltage of the system; j. the design is a squareDA set of all electrical loads;
② inequality constraint
(1) And (3) output power range constraint of the schedulable synchronous unit:
in the formula,andare respectively synchronous machine set gUpper and lower limits of allowable active power;andrespectively an upper limit and a lower limit of an excitation potential of the synchronous unit g;the maximum value of the g stator current of the synchronous unit;
(2) and (3) restricting the operation range of the doubly-fed wind turbine generator:
in the formula, JWThe method comprises the steps of (1) forming a set for all double-fed wind power; vwThe voltage is the stator side machine end voltage of the doubly-fed wind turbine generator w;representing the maximum active power which can be output by the doubly-fed wind turbine generator w under the limitation of natural conditions;the reactance is the stator reactance of the doubly-fed wind turbine generator w;the excitation reactance is the excitation reactance of the double-fed wind turbine generator w;the maximum current of the w rotor side of the doubly-fed wind turbine generator set is obtained; j. the design is a squareWThe method comprises the steps of (1) forming a set by all the double-fed wind power;
(3) and (3) restricting the operation range of the photovoltaic power generation system:
wherein:representing the maximum active power that the photovoltaic power generation system v can output due to the limitation of natural conditions; vvRepresenting the node voltage of the grid-connected side of the photovoltaic power generation system v;representing the maximum load current of the photovoltaic power generation system v inverter; j. the design is a squareVThe photovoltaic power generation system is a set formed by all photovoltaic power generation systems;
(4) and (3) limiting the upper limit and the lower limit of the node voltage amplitude:
in the formula, Vi maxAnd Vi minRespectively representing the upper limit and the lower limit of the voltage amplitude of the node i;
(5) the allowable thermal current range constraint of the power transmission element:
in the formula (19), the compound represented by the formula (I),represents the maximum current of the transmission element l; i isl,ijThe current amplitude of the power transmission element/is represented as:
|Il,ij|＝|Yl|*[(Vli)2+(Vlj)2-2VliVljcosθij]1/2(20)
in the formula (20), YlRepresents the admittance modulus value of the power transmission element l;
(6) and (3) limiting the upper and lower limits of the capacity of the reactive compensation equipment:
in the formula (21), QcRepresenting the reactive power compensated by the reactive compensation device c;andrespectively obtaining maximum compensation reactive power and minimum compensation reactive power of the reactive compensation equipment c; j. the design is a squareCThe method is a set formed by all reactive compensation equipment of the power distribution system.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811196575.XA CN109524955B (en) | 2018-10-15 | 2018-10-15 | Active power distribution network optimization method considering source charge voltage characteristics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811196575.XA CN109524955B (en) | 2018-10-15 | 2018-10-15 | Active power distribution network optimization method considering source charge voltage characteristics |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109524955A true CN109524955A (en) | 2019-03-26 |
| CN109524955B CN109524955B (en) | 2021-11-30 |
Family
ID=65772600
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811196575.XA Active CN109524955B (en) | 2018-10-15 | 2018-10-15 | Active power distribution network optimization method considering source charge voltage characteristics |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109524955B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111130137A (en) * | 2019-12-24 | 2020-05-08 | 上海电力大学 | Micro-grid source-load dynamic frequency control method considering frequency modulation benefit |
| CN112165086A (en) * | 2020-08-12 | 2021-01-01 | 华南理工大学 | Online optimization system of active power distribution network |
| CN112290604A (en) * | 2020-10-15 | 2021-01-29 | 珠海博威电气股份有限公司 | Power distribution network coordination optimization method and device considering load characteristics and storage medium |
| CN116365506A (en) * | 2023-02-03 | 2023-06-30 | 国网山东省电力公司潍坊供电公司 | Energy-saving and loss-reducing optimization method and system for active power distribution network containing distributed photovoltaic |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104377826A (en) * | 2014-11-19 | 2015-02-25 | 云南电网公司电力科学研究院 | Active power distribution network control strategy and method |
| WO2015196745A1 (en) * | 2014-06-25 | 2015-12-30 | 国家电网公司 | Reconstruction method and device for active power distribution network |
| CN107196310A (en) * | 2016-03-15 | 2017-09-22 | 中国电力科学研究院 | The active distribution network Multiple Time Scales optimization method of consideration source net lotus coordination and interaction |
| CN107634518A (en) * | 2017-09-21 | 2018-01-26 | 国网福建省电力有限公司 | The active distribution network economic load dispatching method that a kind of " source net lotus " mutually coordinates |
| CN107688879A (en) * | 2017-10-20 | 2018-02-13 | 云南电网有限责任公司 | A kind of active distribution network distributed power source planing method of consideration source lotus matching degree |
-
2018
- 2018-10-15 CN CN201811196575.XA patent/CN109524955B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015196745A1 (en) * | 2014-06-25 | 2015-12-30 | 国家电网公司 | Reconstruction method and device for active power distribution network |
| CN104377826A (en) * | 2014-11-19 | 2015-02-25 | 云南电网公司电力科学研究院 | Active power distribution network control strategy and method |
| CN107196310A (en) * | 2016-03-15 | 2017-09-22 | 中国电力科学研究院 | The active distribution network Multiple Time Scales optimization method of consideration source net lotus coordination and interaction |
| CN107634518A (en) * | 2017-09-21 | 2018-01-26 | 国网福建省电力有限公司 | The active distribution network economic load dispatching method that a kind of " source net lotus " mutually coordinates |
| CN107688879A (en) * | 2017-10-20 | 2018-02-13 | 云南电网有限责任公司 | A kind of active distribution network distributed power source planing method of consideration source lotus matching degree |
Non-Patent Citations (1)
| Title |
|---|
| 廖剑波: ""源-网-荷"相协调的主动配电网经济调度", 《电力系统及其自动化学报》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111130137A (en) * | 2019-12-24 | 2020-05-08 | 上海电力大学 | Micro-grid source-load dynamic frequency control method considering frequency modulation benefit |
| CN111130137B (en) * | 2019-12-24 | 2021-07-20 | 上海电力大学 | Micro-grid source-load dynamic frequency control method considering frequency modulation benefit |
| CN112165086A (en) * | 2020-08-12 | 2021-01-01 | 华南理工大学 | Online optimization system of active power distribution network |
| CN112290604A (en) * | 2020-10-15 | 2021-01-29 | 珠海博威电气股份有限公司 | Power distribution network coordination optimization method and device considering load characteristics and storage medium |
| CN116365506A (en) * | 2023-02-03 | 2023-06-30 | 国网山东省电力公司潍坊供电公司 | Energy-saving and loss-reducing optimization method and system for active power distribution network containing distributed photovoltaic |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109524955B (en) | 2021-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109524955B (en) | Active power distribution network optimization method considering source charge voltage characteristics | |
| CN111799813B (en) | Marine wind farm reactive power optimization configuration method considering reactive power regulation of wind turbine generator | |
| CN104779638A (en) | Dispatching method and dispatching device for optimizing units in wind power station | |
| CN102255337B (en) | Calculating method for wind power field receiving capability of power grid | |
| Liao et al. | Model predictive control based coordinated voltage control for offshore radial DC-connected wind farms | |
| CN106786766B (en) | A method of the raising wind-powered electricity generation maximum grid connection capacity based on P2G technology | |
| Sharma et al. | Optimistic use of battery energy storage system to mitigate grid disturbances in the hybrid power system. | |
| Tandjaoui et al. | The impact of wind power implantation in transmission systems | |
| Guo et al. | The impact of integrating distributed generations on the losses in the smart grid | |
| CN104037806B (en) | A kind of electric power system tide computational methods based on wind turbine generator basic model | |
| Amita et al. | Power quality comparison of grid connected wind energy system with STATCOM and UPQC | |
| Kamel et al. | Performance Enhancement of an Islanded Microgrid with the Support of Electrical Vehicle and STATCOM Systems. Energies 2023, 16, 1577 | |
| Zhu et al. | An improved newton load flow for distributed generation based on different control strategies | |
| Li et al. | A low voltage ride-through strategy for grid-connected PV converters based on variable power point tracking method | |
| Ni et al. | Optimal control strategy of reactive power and voltage for wind farm based on LinWPSO algorithm | |
| Li et al. | Collaborative Planning and Economic Analyses of High Renewable Penetrated Energy Systems Considering Concentrating Solar Power | |
| CN205051360U (en) | Power transmission network structure that contains distributed generator | |
| CN203175757U (en) | Wind electricity and energy storage combination system | |
| Du et al. | Power distribution strategy considering active power loss for DFIGs wind farm | |
| Li et al. | Analysis of acceptable capacity of microgrid connected to the main power grid | |
| Zanjani et al. | Analysis and simulation of the static synchronous compensation effect in the distribution system with wind farms based on SCIG | |
| Zhang et al. | Typical scenario reactive power optimization of active distribution network based on limit theory | |
| Mhetre et al. | Optimal placement and sizing of distributed generators in unbalanced distribution networks using lion algorithm | |
| Ashglaf et al. | HYBRID WIND-TIDAL ENERGY SYSTEMS–LITERATURE REVIEW | |
| Zhang et al. | Active and Reactive Power Collaborative Optimization by Considering the Deep Peak Regulation of Thermal Power Units |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240822 Address after: No. 101, Building 10, Sports Center, Jinlinwan Garden Community, Qilin District, Qujing City, Yunnan Province 655000 Patentee after: Zhai Haihan Country or region after: China Address before: 650093 No. 253, Xuefu Road, Wuhua District, Yunnan, Kunming Patentee before: Kunming University of Science and Technology Country or region before: China |