CN110021966A - A kind of active distribution network Optimization Scheduling considering dynamic network reconfiguration - Google Patents
A kind of active distribution network Optimization Scheduling considering dynamic network reconfiguration Download PDFInfo
- Publication number
- CN110021966A CN110021966A CN201910172835.8A CN201910172835A CN110021966A CN 110021966 A CN110021966 A CN 110021966A CN 201910172835 A CN201910172835 A CN 201910172835A CN 110021966 A CN110021966 A CN 110021966A
- Authority
- CN
- China
- Prior art keywords
- distribution network
- active distribution
- power
- node
- moment
- 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.)
- Pending
Links
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
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/50—Controlling the sharing of the out-of-phase component
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
- Y02A30/60—Planning or developing urban green infrastructure
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of active distribution network Optimization Schedulings for considering dynamic network reconfiguration, comprising: (1) obtains active distribution network parameter, obtain active distribution network power flow equation constraint condition；(2) consider dynamic network reconfiguration, correct active distribution network power flow equation constraint condition；(3) active distribution network irradiation structure constraint condition is obtained according to active distribution network irradiation structure；(4) with the minimum objective function of active distribution network integrated operation cost, active distribution network Optimal Operation Model is constructed according to every constraint condition；(5) constraint condition non-convex in master mould is subjected to convex relaxationization, Nonlinear Constraints are linearized；(6) model after conversion is solved, obtains active distribution network Optimized Operation scheme.The present invention fully considers the controllable resources of active distribution network, improves the performance driving economy of active distribution network, thus solves the problems, such as that large-scale distributed plant-grid connection active distribution network is difficult to dissolve and increases with voltage.
Description
Technical field
The invention belongs to active distribution networks to optimize running technology field, more particularly, to a kind of consideration dynamic network weight
The active distribution network Optimization Scheduling of structure.
Background technique
In recent years, it is reply climate change, solves energy crisis and promote energy-saving and emission-reduction, greatly develops wind-powered electricity generation and the sun
The renewable energy such as energy power generation have become common recognition and important means both at home and abroad.Country will continue to greatly develop wind-powered electricity generation and solar energy
The distributed development and on-site elimination of renewable energy are realized in power generation by many kinds of measures.With renewable energy on a large scale with
The mode of distributed power generation accesses power distribution network, in the partial period of certain localized power distribution nets, or even distributed generation resource infiltration occurs
Saturating rate is more than 100% situation, and some distributed generation resources connect in active distribution network end can bring serious overvoltage to ask
Topic, influences active distribution network to the digestion capability of renewable energy.
In order to solve the problems, such as that large-scale distributed plant-grid connection power distribution network is brought, the concept of " active distribution network " is met the tendency of
And it gives birth to.Have part research institution and scholar and has studied active distribution network coordinated operation from different angles such as economy, safeties
Mechanism, but as the permeability of distributed generation resource is higher and higher, power distribution network needs more flexible controllable resources.As one kind
The scheduling resource of network topology flexibility and changeability, dynamic network reconfiguration can effectively promote the economic operation level of power distribution network and divide
The consumption rate of cloth power supply.China's power distribution network just develops towards intelligentized direction, and " action plan is transformed in distribution network construction
(2015-the year two thousand twenty) " it will clearly accelerate Modern power distribution net construction retrofit comprehensively, to meet the hair of new energy, distributed generation resource
Exhibition and access demand.Wherein, the application of intellectual remote controlling switch provides condition for the implementation of dynamic network reconfiguration.
Since existing research mostly lays particular emphasis on the coordination in " source-lotus-storage ", as the access control of distributed generation resource, energy storage and
Coordination optimization, microgrid and coordinated operation of active distribution network of the resources such as flexible load etc., about active distribution network network structure
The research for adjusting aspect is less, and the potentiality of active distribution network consumption distributed generation resource are not not fully exerted.And active distribution
The research of net network reconfiguration is confined to the long time scale optimization of active distribution network mostly, using network structure regulation as in a few days
The concern that the dynamic network reconfiguration of scheduling resource is subject to is seldom.In addition, existing in terms of the irradiation structure constraint of active distribution network
Document mostly uses the method for exhaustion, does not provide stringent mathematical expression.
Summary of the invention
Aiming at the above defects or improvement requirements of the prior art, the invention proposes a kind of consideration dynamic network reconfigurations
Active distribution network Optimization Scheduling improves active distribution network its object is to fully consider the controllable resources of active distribution network
Performance driving economy, thus solve large-scale distributed plant-grid connection active distribution network and be difficult to dissolve and the raised technology of voltage
Problem.
To achieve the above object, the present invention provides a kind of active distribution network Optimized Operations for considering dynamic network reconfiguration
Method, comprising:
(1) active distribution network parameter is obtained, active distribution network power flow equation constraint condition is constructed；
(2) consider dynamic network reconfiguration, correct active distribution network power flow equation constraint condition, obtain revised active
Distribution power flow equation constraint condition；
(3) active distribution network irradiation structure constraint condition is obtained according to active distribution network irradiation structure；
(4) with the minimum objective function of active distribution network integrated operation cost, according to revised active distribution network trend
Equation constraint condition, active distribution network irradiation structure constraint condition and active distribution network are safely operated constraint condition, and building is actively
Power distribution network Optimal Operation Model；
(5) the non-convex constraint condition in model is subjected to convex relaxationization, the Nonlinear Constraints in model is subjected to line
Property；
(6) model after conversion is solved, obtains active distribution network Optimized Operation scheme.
Further, active distribution network parameter include distributed electrical source dates, distributed generation resource for 24 hours go out force data, master
Dynamic distribution network structure parameter, active distribution network node load data, active distribution network critical point power limit, reactive power compensator
Parameter, switch motion maximum number of times, every cost coefficient；
Wherein distributed electrical source dates include distributed generation resource maximum active power output Pi,max, take power limitation factor model
Enclose the power factor limit value of operation
It is predicted value that distributed generation resource goes out force data for 24 hours
Active distribution network parameter of double--layer grids includes the conductance g of active distribution network branch i-jij, susceptance bijAnd upper current limit
Iij,max, active distribution network initial primary topology and interconnection configuration information, active distribution network node voltage bound Vi,maxWith
Vi,min；
Active distribution network node load data include burden with power P at t moment node iL,i,tWith load or burden without work QL,i,t；
Active distribution network critical point power limit includes the active exchange power upper lower limit value P in active distribution network critical pointT,min、
PT,maxWith reactive power exchange power upper lower limit value QT,minWith QT,max；
Reactive power compensator parameter includes the minimum value Q of SVC reactive compensation at node iSVC,i,minAnd maximum value
QSVC,i,max；
Switch motion maximum number of times is SWmax；
Every cost coefficient includes Web-based exercise coefficient ηL, abandonment/light unit electricity punishment cost coefficient ηRDG, switch dynamic
Make primary cost coefficient ηS。
Further, active distribution network power flow equation constraint condition are as follows:
Wherein, PG,i,tFor the active injection of distributed generation resource or power transmission network at t moment node i, PL,i,tFor t moment
Burden with power at node i, N (i) are the node set being connected with node i, Pij,tFor the wattful power flowed through on t moment branch i-j
Rate, to be positive from node i outflow；QG,i,tFor the idle injection of power supply at t moment node i, QL,i,tFor at t moment node i without
Workload, QSVC,i,tFor SVC reactive compensation power at t moment node i, Qij,tFor the idle function flowed through on t moment branch i-j
Rate, to be positive from node i outflow；
Pij,tAnd Qij,tFormula can be expressed according to branch apparent energy to derive:
Wherein, Sij,tFor t moment branch i-j apparent energy,For the conjugate phasors of t moment branch i-j electric current, gijWith
bijThe respectively conductance and susceptance of branch i-j, Vi,tFor the voltage magnitude of t moment node i, θij,t=θi,t-θj,tFor t moment section
Phase difference of voltage between point i, j；
It is available after real part, imaginary part separation:
Further, revised active distribution network power flow equation constraint condition are as follows:
Wherein, Zij,tState variable is cut-off for t moment branch i-j, when branch i-j is disconnected, Zij,t=0, it is corresponding
Branch power Zij,t·Pij,t=0.
Further, active distribution network irradiation structure constraint condition are as follows:
Wherein, Eij,tAnd Zij,tIt is Boolean variable, Eij,tFor the direction of tide variable of t moment branch i-j, Eij,t=1
Indicate that t moment node j is the parent node of node i, trend flow direction is from node j to node i；Zij,tFor the company of t moment branch i-j
Logical state, Zij,t=1 indicates t moment branch i-j connection.
Further, active distribution network is safely operated constraint condition further include:
The constraint of voltage magnitude bound:
Vi,min≤Vi,t≤Vi,max
Wherein, Vi,maxWith Vi,minThe upper lower limit value of voltage magnitude respectively at node i.
The constraint of the branch current upper limit:
Wherein, Iij,tFor branch current magnitudes, Iij,maxFor branch current magnitudes upper limit value.
Distributed generation resource units limits:
Wherein, PRDG,i,tAnd QRDG,i,tThe active and idle actual power of distributed generation resource respectively at t moment node i,For distributed generation resource power output predicted value at t moment node i；Power limitation factor range is taken to transport distributed generation resource
Row,Cloth electrical source power factor limits value is punished for node i.
Active distribution network critical point power constraint:
Wherein, PT,tThe active power of the same level power distribution network, P are flowed into from power distribution network root node for t momentT,minAnd PT,maxPoint
It Wei not the active exchange power upper lower limit value in active distribution network critical point；QT,tThe same level distribution is flowed into from power distribution network root node for t moment
The reactive power of net, QT,minAnd QT,maxRespectively active distribution network critical point reactive power exchange power upper lower limit value.
Static passive compensation device SVC constraint:
QSVC,i,min≤QSVC,i,t≤QSVC,i,max
Wherein, QSVC,i,tFor t moment node i reactive compensation power size, QSVC,i,minAnd QSVC,i,maxRespectively node i
Locate the minimum value and maximum value of SVC reactive compensation.
Switch motion count constraint:
ΔZij,tIndicate t moment with respect to t-1 moment switch change number, SWmaxFor switch motion maximum number of times.
Further, the expression formula of objective function are as follows:
Wherein, PL,tFor t moment active distribution network network loss, ηLFor Web-based exercise coefficient；For the distribution of t moment i-node
The abandoning amount of formula power supply, ηRDGFor abandonment/light unit electricity punishment cost coefficient；ΔZij,tIt is t moment with respect to last moment network
Switch motion number, ηSFor the primary cost coefficient of switch motion, ΩRDGFor distributed power access point set, ΩEFor power grid branch
Road set.
Further, non-convex constraint condition includes power flow equation constraint condition, and nonlinear restriction includes revised trend
The absolute value expression of the product term, branch current upper limit constraint condition and switch motion count constraint condition of equation constraint condition
Formula,
The convex relaxation of power flow equation constraint condition are as follows:
Its standard second order tapered are as follows:
Wherein, Ui,t、Wij,tAnd Tij,tMeet:
The product term of revised power flow equation constraint condition linearizes:
Wherein,Node voltage variable U respectively relevant to node i, ji,t、Uj,tIt corresponds to and branch
" i-j " relevant branch voltage variable meets following constraint:
The linearisation of branch current upper limit constraint condition are as follows:
For the absolute value expression of switch motion count constraint condition, following expression can be replaced with:
In this way, former Optimized model is converted into a MIXED INTEGER Second-order cone programming model, available good solution.
Wherein, optimization software is existing optimization problem solving software, the tool box yalmip in such as CPLEX, Matlab, Gurobi
Deng.
The invention proposes a kind of active distribution network Optimized Operation strategies for considering dynamic network reconfiguration, it is intended to sufficiently benefit
With the schedulable resource of active distribution network, distributed generation resource consumption is further increased.Accordingly give the constraint of power distribution network irradiation structure
Strict mathematical expression formula and model method for solving.The deficiency of prior art is compensated for well.
Contemplated above technical scheme through the invention can obtain following compared with prior art
The utility model has the advantages that
1, the present invention is with the minimum target of power distribution network synthesis operating cost, meter and power flow equation constraint condition and switch
Switching number constraint condition, while the adjustable resource such as integrate distributed generation resource, reactive power compensator, establish consideration Dynamic Networks
The active distribution network Optimal Operation Model of network reconstruct, the model can accurately describe to consider large-scale distributed plant-grid connection
When active distribution network operation demand；
2, the model conversation method proposed by the invention based on second order cone relaxation and the linearisation of variable product can be to non-
Convex power flow equation and variable product term carry out convexification processing, and then convert master mould to the MIXED INTEGER two for being easy solution
Rank bores problem, greatly reduces solution difficulty；
3, the active distribution network Optimized Operation strategy proposed by the present invention for considering dynamic network reconfiguration can make full use of master
The dynamic schedulable resource of power distribution network, further increases active distribution network distributed generation resource digestion capability；Idle benefit can be passed through simultaneously
The coordination for repaying device and dynamic network reconfiguration better solves on voltage caused by large-scale distributed plant-grid connection power distribution network
The problem of liter, has good economy and applicability.
Detailed description of the invention
Fig. 1 is a kind of distribution of active distribution network Optimization Scheduling for considering dynamic network reconfiguration provided by the invention
The distribution curve for 24 hours of formula power supply and load；
Fig. 2 is a kind of active of active distribution network Optimization Scheduling for considering dynamic network reconfiguration provided by the invention
Distribution network system topology diagram.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments,
The present invention will be described in further detail.It should be appreciated that specific embodiment described herein is only used to explain this hair
It is bright, it is not intended to limit the present invention.
A kind of active distribution network Optimized Operation strategy of consideration dynamic network reconfiguration of the embodiment of the present invention calculates packet
Include following steps:
(1) obtain that all distributed electrical source dates inside active distribution network are as shown in table 1, distributed generation resource and load
Distribution curve is as shown in Figure 1, active distribution network parameter of double--layer grids is as shown in table 2 for 24 hours, and topological structure is as shown in Fig. 2, active distribution
Net node load data as shown in table 3, active distribution network critical point power limit it is as shown in table 4, reactive power compensator parameter be-
100~300kVar, switch motion maximum number of times are 24 times, every cost coefficient is as shown in table 5；
Table 1
| Distributed generation resource | Access node | Rated power/MW | Power factor |
| PV | 6 | 0.5 | 0 |
| Wind | 7 | 0.6 | 0.9 |
| Wind | 13 | 1 | 0.9 |
| Wind | 18 | 1 | 0.9 |
| Wind | 28 | 0.6 | 0.9 |
| PV | 33 | 0.5 | 0 |
Table 2
| First node | Tail node | Resistance/Ω | Reactance/Ω | First node | Tail node | Resistance/Ω | Reactance/Ω |
| 1 | 2 | 0.0922 | 0.047 | 20 | 21 | 0.4095 | 0.4784 |
| 2 | 3 | 0.493 | 0.2511 | 21 | 22 | 0.7089 | 0.9373 |
| 3 | 4 | 0.366 | 0.1864 | 3 | 23 | 0.4512 | 0.3083 |
| 4 | 5 | 0.3811 | 0.1941 | 23 | 24 | 0.898 | 0.7091 |
| 5 | 6 | 0.819 | 0.707 | 24 | 25 | 0.896 | 0.7011 |
| 6 | 7 | 0.1872 | 0.6188 | 6 | 26 | 0.203 | 0.1034 |
| 7 | 8 | 0.7114 | 0.2351 | 26 | 27 | 0.2842 | 0.1447 |
| 8 | 9 | 1.03 | 0.74 | 27 | 28 | 1.059 | 0.9337 |
| 9 | 10 | 1.044 | 0.74 | 28 | 29 | 0.8042 | 0.7006 |
| 10 | 11 | 0.1966 | 0.065 | 29 | 30 | 0.5075 | 0.2585 |
| 11 | 12 | 0.3744 | 0.1238 | 30 | 31 | 0.9744 | 0.963 |
| 12 | 13 | 1.468 | 1.155 | 31 | 32 | 0.3105 | 0.3619 |
| 13 | 14 | 0.5416 | 0.7129 | 32 | 33 | 0.341 | 0.5302 |
| 14 | 15 | 0.591 | 0.526 | 8 | 21 | 2 | 2 |
| 15 | 16 | 0.7463 | 0.545 | 9 | 15 | 2 | 2 |
| 16 | 17 | 1.289 | 1.721 | 12 | 22 | 2 | 2 |
| 17 | 18 | 0.732 | 0.574 | 18 | 33 | 0.5 | 0.5 |
| 2 | 19 | 0.164 | 0.1565 | 25 | 29 | 0.5 | 0.5 |
| 19 | 20 | 1.5042 | 1.3554 | - | - | - | - |
Table 3
Table 4
| Parameter name | Active lower limit/MW | The active upper limit/MW | Idle lower limit/MVar | The idle upper limit/MVar |
| Power limit | -5 | 10 | -2 | 8 |
Table 5
| Parameter name | Relative value |
| Web-based exercise | 50 |
| Abandonment/light punishment cost | 1 |
| Switch motion cost | 2 |
(2) it establishes active distribution network power flow equation joint constraint condition and considers the amendment constraint item of dynamic network reconfiguration
Part:
Wherein, PG,i,tFor power supply (distributed generation resource or power transmission network) active injection at t moment node i, PL,i,tWhen for t
Burden with power at node i is carved, N (i) is the node set being connected with node i, Pij,tIt is active to be flowed through on t moment branch i-j
Power, to be positive from node i outflow；QG,i,tFor the idle injection of power supply at t moment node i, QL,i,tAt t moment node i
Load or burden without work, QSVC,i,tFor SVC reactive compensation power at t moment node i, Qij,tFor the idle function flowed through on t moment branch i-j
Rate, to be positive from node i outflow.Active/idle total outflow function of all routes that t moment is connected with node i is indicated on the right of equation
The sum of rate.In addition, Pij,tAnd Qij,tFormula can be expressed according to branch apparent energy to derive:
The amendment active distribution network power flow equation constraint condition for considering dynamic network reconfiguration are as follows:
Wherein, Zij,tState variable is cut-off for t moment branch i-j, when branch i-j is disconnected, Zij,t=0, it is corresponding
Branch power Zij,t·Pij,t=0.
(3) constraint condition of active distribution network radial structure is established:
Wherein, Eij,tAnd Zij,tIt is Boolean variable.Eij,tFor the direction of tide variable of t moment branch i-j, Eij,t=1
Indicate that t moment node j is the parent node of node i, trend flow direction is from node j to node i；Zij,tFor the company of t moment branch i-j
Logical state, Zij,t=1 indicates t moment branch i-j connection.First item limits E in constraintij,tAnd Eji,tAt most only one is
1 and maintain the amphicheirality of branch i-j trend；Second formula indicate in addition to power transformation tiny node remaining all node all only one
Parent node；Section 3 indicates that power transformation tiny node does not have parent node.
(4) the active distribution network Optimal Operation Model for considering dynamic network reconfiguration is established, including objective function and correlation are about
Beam；
Wherein, consider the objective function of the active distribution network Optimal Operation Model of dynamic network reconfiguration are as follows:
Specifically, active distribution network irradiation structure is safely operated constraint condition further include:
The constraint of voltage magnitude bound
Vi,min≤Vi,t≤Vi,max
Wherein, Vi,maxWith Vi,minThe upper lower limit value of voltage magnitude respectively at node i.
The constraint of the branch current upper limit
In formula, Iij,tFor branch current magnitudes, Iij,maxFor branch current magnitudes upper limit value.
Distributed generation resource units limits
In formula, PRDG,i,tAnd QRDG,i,tThe active and idle actual power of distributed generation resource respectively at t moment node i,For distributed generation resource power output predicted value at t moment node i；Power limitation factor range is taken to transport distributed generation resource
Row,Cloth electrical source power factor limits value is punished for node i.
Active distribution network critical point power constraint
In order to inhibit the influence caused by power transmission network of the power swing of active distribution network, need active distribution network root section
The critical point exchange power constraint of point takes into account, it may be assumed that
In formula, PT,tThe active power of the same level power distribution network, P are flowed into from power distribution network root node for t momentT,minAnd PT,maxPoint
It Wei not the active exchange power upper lower limit value in active distribution network critical point；QT,tThe same level distribution is flowed into from power distribution network root node for t moment
The reactive power of net, QT,minAnd QT,maxRespectively active distribution network critical point reactive power exchange power upper lower limit value.
Static passive compensation device SVC constraint
QSVC,i,min≤QSVC,i,t≤QSVC,i,max
In formula, QSVC,i,tFor t moment node i reactive compensation power size, QSVC,i,minAnd QSVC,i,maxRespectively node i
Locate the minimum value and maximum value of SVC reactive compensation.
Switch motion count constraint
In formula, Δ Zij,tIndicate t moment with respect to t-1 moment switch change number, SWmaxFor switch motion maximum number of times.
Wherein Δ Zij,tIt is expressed as follows:
ΔZij,t=| Zij,t-Zij,t-1|
(5) the non-convex expression formula such as secondary, trigonometric function will be born in master mould and carries out convex relaxation conversion, is translated into mark
Quasi- second order tapered specifically includes power flow equation and constrains convex relaxation and the convex relaxation of upper current limit constraint；
The convex relaxation of power flow equation constraint condition are as follows:
Its standard second order tapered are as follows:
Wherein, Ui,t、Wij,tAnd Tij,tMeet following formula；
(6) it is linearized by variable product and linearizes nonlinear restriction, and using optimization software to the optimization after conversion
Model is solved, and the active distribution network Optimized Operation strategy for considering dynamic network reconfiguration is obtained.
Power flow equation after linearisation are as follows:
Wherein, in formula,Node voltage variable U respectively relevant to node i, ji,t、Uj,tIt corresponds to and props up
Relevant " branch voltage " variable in road " i-j " meets following constraint:
Upper current limit constraint condition substitution are as follows:
For the absolute value expression of switch motion count constraint condition, it is replaced with following expression:
In this way, former Optimized model is converted into a MIXED INTEGER Second-order cone programming model, the present embodiment using CPLEX and
The tool box yalmip in Matlab is solved.
Example is provided with 5 kinds of scenes to analyze the validity of proposed method.
Scene 1: at single time point, only RDG active reactive is controllable in network；
Scene 2: static network reconstruct is considered on the basis of scene 1；
Scene 3: contributing for the load and RDG of variation in 24 hours, considers hour grade DNR；
Scene 4: on the basis of scene 3, further consider switch motion count constraint proposed in this paper and filter out
Key switch；
Scene 5: on the basis of scene 4, influence of the reactive power compensator to DNR is further probed into.
Table 6
Table 6 is the Comparative result of scene 1 and scene 2, is after considering network reconfiguration it can be seen from Comparative result in table 6
System economic operation level is obviously improved, and the RDG consumption amount of active distribution network totality increases 0.6515MWh, active distribution network net
Damage reduces 29.72%, and active distribution network minimum voltage level is promoted.The distribution it can be seen from the comparison of scene 1 and scene 2
Network reconfiguration has huge potentiality to the promotion of active distribution network economic operation level, and numerical results are to study hour grade
Dynamic restructuring provides data supporting.
Table 7
In analytical table 73 partial period of scene result it is found that RDG power output it is little when, with optimal static reconfiguration strategy
It comparing, DNR is unobvious to active distribution network network loss reduction effect, because RDG can be dissolved completely at this time, a small amount of variation of RDG
It is faint to active distribution network influence on tidal flow；When RDG power output is sufficiently large (permeability >=100%), with optimal static policies phase
Than dynamic restructuring reduces significant effect to active distribution network network loss, and RDG consumption amount also obviously increases.As it can be seen that DNR goes out in RDG
The power biggish period can optimize active distribution network network structure, can preferably cope with load by the adjustment of network switching
With the variation of RDG power output.
Table 8
Table 9
Table 8 is scene 4 and 5 partial period switch state of scene as a result, table 9 is 5 partial period voltage of scene 4 and scene
It is horizontal.As can be seen from Table 8, compared with scene 4, switch state is varied in scene 5, while RDG consumption is also increased
Add.As can be seen from Table 9, network ceiling voltage has reached the voltage magnitude upper limit at node 13 in scene 4；And scene 5
The voltage's distribiuting of network is optimized by reactive power compensator, ceiling voltage does not reach the amplitude upper limit.Known to comprehensive analysis:
There is no reactive power compensator in scene 4, the consumption of RDG receives the limitation of active distribution network voltage constraint, and scene 5 can be with
Flexibly call reactive power compensator to improve active distribution network voltage level by dynamic network reconfiguration, thus in scene 5,
Active distribution network can dissolve more RDG.
Numerical results show: the method that this patent proposes can utilize active distribution network network topology structure well
Dynamic restructuring, the distributed generation resource consumption for further increasing power distribution network is horizontal, improves active distribution network operating status.
As it will be easily appreciated by one skilled in the art that the foregoing is merely illustrative of the preferred embodiments of the present invention, not to
The limitation present invention, any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should all wrap
Containing within protection scope of the present invention.
Claims (8)
1. a kind of active distribution network Optimization Scheduling for considering dynamic network reconfiguration, which is characterized in that include the following steps:
(1) active distribution network parameter is obtained, active distribution network power flow equation constraint condition is obtained；
(2) consider dynamic network reconfiguration, correct the active distribution network power flow equation constraint condition, obtain revised active and match
Electric network swim equation constraint condition；
(3) active distribution network irradiation structure constraint condition is obtained according to active distribution network irradiation structure；
(4) with the minimum objective function of active distribution network integrated operation cost, according to the revised active distribution network trend
Equation constraint condition, the active distribution network irradiation structure constraint condition and active distribution network are safely operated constraint condition, building
Active distribution network Optimal Operation Model；
(5) the non-convex constraint condition in the model is subjected to convex relaxationization, by the Nonlinear Constraints in the model into
Row linearisation；
(6) model after conversion is solved, obtains active distribution network Optimized Operation scheme.
2. the method as described in claim 1, which is characterized in that the active distribution network parameter include distributed electrical source dates,
Distributed generation resource goes out force data, active distribution network parameter of double--layer grids, active distribution network node load data, active distribution network for 24 hours
Critical point power limit, reactive power compensator parameter, switch motion maximum number of times, every cost coefficient；
Wherein the distributed electrical source dates include distributed generation resource maximum active power output Pi,max, take power limitation factor range
The power factor limit value of operation
It is predicted value that the distributed generation resource goes out force data for 24 hours
The active distribution network parameter of double--layer grids includes the conductance g of active distribution network branch i-jij, susceptance bijAnd upper current limit
Iij,max, active distribution network initial primary topology and interconnection configuration information, active distribution network node voltage bound Vi,maxWith
Vi,min；
The active distribution network node load data include burden with power P at t moment node iL,i,tWith load or burden without work QL,i,t；
Active distribution network critical point power limit includes the active exchange power upper lower limit value P in active distribution network critical pointT,min、PT,max
With reactive power exchange power upper lower limit value QT,minWith QT,max；
The reactive power compensator parameter includes the minimum value Q of SVC reactive compensation at node iSVC,i,minWith maximum value QSVC,i,max；
The switch motion maximum number of times is SWmax；
The items cost coefficient includes Web-based exercise coefficient ηL, abandonment or abandon light unit electricity punishment cost coefficient ηRDG, open
Pass acts primary cost coefficient ηS。
3. the method as described in claim 1, which is characterized in that the active distribution network power flow equation constraint condition are as follows:
Wherein, PG,i,tFor the active injection of distributed generation resource or power transmission network at t moment node i, PL,i,tFor t moment node i
Locate burden with power, N (i) is the node set being connected with node i, Pij,tFor the active power flowed through on t moment branch i-j, with
It is positive from node i outflow；QG,i,tFor the idle injection of power supply at t moment node i, QL,i,tFor load or burden without work at t moment node i,
QSVC,i,tFor SVC reactive compensation power at t moment node i, Qij,tFor the reactive power flowed through on t moment branch i-j, with from section
Point i outflow is positive；
Pij,tAnd Qij,tFormula can be expressed according to branch apparent energy to derive:
Wherein, Sij,tFor t moment branch i-j apparent energy,For the conjugate phasors of t moment branch i-j electric current, gijAnd bijPoint
Not Wei branch i-j conductance and susceptance, Vi,tFor the voltage magnitude of t moment node i, θij,t=θi,t-θj,tFor t moment node i, j
Between phase difference of voltage；
Real part, imaginary part separation:
P can be obtainedij,tAnd Qij,t。
4. the method as described in claim 1, which is characterized in that the revised active distribution network power flow equation constraint condition
Are as follows:
Wherein, Zij,tState variable is cut-off for t moment branch i-j, when branch i-j is disconnected, Zij,t=0, corresponding branch function
Rate Zij,t·Pij,t=0.
5. the method as described in claim 1, which is characterized in that the active distribution network irradiation structure constraint condition are as follows:
Wherein, Eij,tAnd Zij,tIt is Boolean variable, Eij,tFor the direction of tide variable of t moment branch i-j, Eij,tWhen=1 expression t
The parent node that node j is node i is carved, trend flow direction is from node j to node i；Zij,tFor the connected state of t moment branch i-j,
Zij,t=1 indicates t moment branch i-j connection.
6. the method as described in claim 1, which is characterized in that the active distribution network is safely operated constraint condition and includes:
The constraint of voltage magnitude bound:
Vi,min≤Vi,t≤Vi,max
The constraint of the branch current upper limit:
Distributed generation resource units limits
Active distribution network critical point power constraint
SVC constraint:
QSVC,i,min≤QSVC,i,t≤QSVC,i,max
Switch motion count constraint:
Wherein, Vi,maxWith Vi,minThe upper lower limit value of voltage magnitude, I respectively at node iij,tFor branch current magnitudes, Iij,maxFor
Branch current magnitudes upper limit value, PRDG,i,tAnd QRDG,i,tThe active and idle practical function of distributed generation resource respectively at t moment node i
Rate,For distributed generation resource power output predicted value at t moment node i；Power limitation factor range is taken to transport distributed generation resource
Row,Cloth electrical source power factor limits value, P are punished for node iT,tThis gradation is flowed into from power distribution network root node for t moment
The active power of power grid, PT,minAnd PT,maxThe respectively active exchange power upper lower limit value in active distribution network critical point；QT,tFor t moment
The reactive power of the same level power distribution network, Q are flowed into from power distribution network root nodeT,minAnd QT,maxThe respectively idle friendship in active distribution network critical point
Change power upper lower limit value, QSVC,i,tFor t moment node i reactive compensation power size, QSVC,i,minAnd QSVC,i,maxRespectively node i
Locate the minimum value and maximum value of SVC reactive compensation, Δ Zij,tIndicate t moment with respect to t-1 moment switch change number, SWmaxTo open
Close the action frequency upper limit.
7. the method as described in claim 1, which is characterized in that the expression formula of the objective function are as follows:
Wherein, PL,tFor t moment active distribution network network loss, ηLFor Web-based exercise coefficient；For t moment i-node distributed generation resource
Abandoning amount, ηRDGFor abandonment/light unit electricity punishment cost coefficient；ΔZij,tIt is acted for t moment with respect to last moment network switching
Number, ηSFor the primary cost coefficient of switch motion, ΩRDGFor distributed power access point set, ΩEFor grid branch set.
8. the method as described in claim 1, which is characterized in that the non-convex constraint condition includes power flow equation constraint condition,
The nonlinear restriction includes product term, branch current upper limit constraint condition and the switch of revised power flow equation constraint condition
The absolute value expression of switching number constraint condition,
The convex relaxation of power flow equation constraint condition are as follows:
Its standard second order tapered are as follows:
Wherein, Ui,t、Wij,tAnd Tij,tMeet following constraint:
The product term of the revised power flow equation constraint condition linearizes are as follows:
Wherein,Meet following constraint:
Node voltage variable U respectively relevant to node i, ji,t、Uj,tCorrespond to branch relevant to branch " i-j "
Road voltage quantities；
The branch current upper limit constraint condition linearisation are as follows:
The absolute value expression of the switch motion count constraint condition linearizes are as follows:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910172835.8A CN110021966A (en) | 2019-03-07 | 2019-03-07 | A kind of active distribution network Optimization Scheduling considering dynamic network reconfiguration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910172835.8A CN110021966A (en) | 2019-03-07 | 2019-03-07 | A kind of active distribution network Optimization Scheduling considering dynamic network reconfiguration |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110021966A true CN110021966A (en) | 2019-07-16 |
Family
ID=67189379
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910172835.8A Pending CN110021966A (en) | 2019-03-07 | 2019-03-07 | A kind of active distribution network Optimization Scheduling considering dynamic network reconfiguration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110021966A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110518580A (en) * | 2019-08-15 | 2019-11-29 | 上海电力大学 | A kind of active distribution network running optimizatin method for considering microgrid and actively optimizing |
| CN111934320A (en) * | 2020-08-20 | 2020-11-13 | 三峡大学 | Active power distribution network rapid reconstruction method based on linear power flow equation and improved firework algorithm |
| CN112039067A (en) * | 2020-09-01 | 2020-12-04 | 国网河北省电力有限公司邢台供电分公司 | Power distribution network new energy power generation utilization rate optimization method and terminal equipment |
| CN112909926A (en) * | 2021-01-25 | 2021-06-04 | 浙江大学 | Radial constraint set generation method and device for power distribution system optimization problem |
| CN113517698A (en) * | 2021-07-21 | 2021-10-19 | 华北电力科学研究院有限责任公司 | Optimal power flow convexity control method and device for active power distribution network |
| CN113609699A (en) * | 2021-08-17 | 2021-11-05 | 上海交通大学 | Calculation method and system for alternating current power flow model of radial power distribution network |
| CN113705883A (en) * | 2021-08-26 | 2021-11-26 | 华南理工大学 | Alternating current power distribution network real-time schedulable domain calculation method based on compact convex relaxation |
| CN113937829A (en) * | 2021-11-16 | 2022-01-14 | 华北电力大学 | Active power distribution network multi-target reactive power control method based on D3QN |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105023058A (en) * | 2015-07-07 | 2015-11-04 | 天津大学 | Power distribution network intelligent soft switch operation optimization method with simultaneous consideration of switch motion |
| CN106208160A (en) * | 2016-07-28 | 2016-12-07 | 东南大学 | Dispatching method based on the sale of electricity company region within the jurisdiction power distribution network that second order cone optimizes |
| CN106803677A (en) * | 2017-04-11 | 2017-06-06 | 四川大学 | A kind of active distribution network voltage management-control method and system based on distributed power source |
-
2019
- 2019-03-07 CN CN201910172835.8A patent/CN110021966A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105023058A (en) * | 2015-07-07 | 2015-11-04 | 天津大学 | Power distribution network intelligent soft switch operation optimization method with simultaneous consideration of switch motion |
| CN106208160A (en) * | 2016-07-28 | 2016-12-07 | 东南大学 | Dispatching method based on the sale of electricity company region within the jurisdiction power distribution network that second order cone optimizes |
| CN106803677A (en) * | 2017-04-11 | 2017-06-06 | 四川大学 | A kind of active distribution network voltage management-control method and system based on distributed power source |
Non-Patent Citations (3)
| Title |
|---|
| ZHUANG TIAN ET AL.: "Mixed-integer second-order cone programing model for VAR optimisation and network reconfiguration in active distribution networks", 《IET GENERATION, TRANSMISSION & DISTRIBUTION》 * |
| 李超等: "考虑动态网络重构的主动配电网优化运行策略研究", 《电工技术学报》 * |
| 马丽叶等: "主动配电网中计及网络重构的分布式电源鲁棒规划", 《电力系统自动化》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110518580B (en) * | 2019-08-15 | 2023-04-28 | 上海电力大学 | Active power distribution network operation optimization method considering micro-grid active optimization |
| CN110518580A (en) * | 2019-08-15 | 2019-11-29 | 上海电力大学 | A kind of active distribution network running optimizatin method for considering microgrid and actively optimizing |
| CN111934320A (en) * | 2020-08-20 | 2020-11-13 | 三峡大学 | Active power distribution network rapid reconstruction method based on linear power flow equation and improved firework algorithm |
| CN111934320B (en) * | 2020-08-20 | 2023-08-01 | 三峡大学 | Active power distribution network rapid reconstruction method based on linear tide equation and improved firework algorithm |
| CN112039067A (en) * | 2020-09-01 | 2020-12-04 | 国网河北省电力有限公司邢台供电分公司 | Power distribution network new energy power generation utilization rate optimization method and terminal equipment |
| CN112909926A (en) * | 2021-01-25 | 2021-06-04 | 浙江大学 | Radial constraint set generation method and device for power distribution system optimization problem |
| CN112909926B (en) * | 2021-01-25 | 2022-06-17 | 浙江大学 | Radial constraint set generation method and device for power distribution system optimization problem |
| CN113517698A (en) * | 2021-07-21 | 2021-10-19 | 华北电力科学研究院有限责任公司 | Optimal power flow convexity control method and device for active power distribution network |
| CN113517698B (en) * | 2021-07-21 | 2024-04-19 | 华北电力科学研究院有限责任公司 | Active power distribution network optimal power flow salifying control method and device |
| CN113609699A (en) * | 2021-08-17 | 2021-11-05 | 上海交通大学 | Calculation method and system for alternating current power flow model of radial power distribution network |
| CN113609699B (en) * | 2021-08-17 | 2023-07-21 | 上海交通大学 | Calculation method and system of alternating current power flow model of radial power distribution network |
| CN113705883A (en) * | 2021-08-26 | 2021-11-26 | 华南理工大学 | Alternating current power distribution network real-time schedulable domain calculation method based on compact convex relaxation |
| CN113937829A (en) * | 2021-11-16 | 2022-01-14 | 华北电力大学 | Active power distribution network multi-target reactive power control method based on D3QN |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110021966A (en) | A kind of active distribution network Optimization Scheduling considering dynamic network reconfiguration | |
| CN101281637B (en) | Electric power system optimal power flow method based on hyperplane form safety field periphery | |
| CN103150606B (en) | A kind of distributed power source optimal load flow optimization method | |
| CN105119280B (en) | Alternating current-direct current mixed structure active power distribution network running optimizatin method based on cone optimization | |
| CN105375513B (en) | A kind of 110 kilovolts of wind power field automatic voltage controls based on real-time on-line equivalent | |
| CN107947175B (en) | Micro-grid economic dispatching method based on distributed network control | |
| CN103593711B (en) | A kind of distributed power source Optimal Configuration Method | |
| CN110034587B (en) | Optimized scheduling method | |
| CN109066694A (en) | Multiple target tide optimization method containing the electric system of flow controller between line | |
| CN106777586B (en) | A kind of operation domain method for solving calculating Distributed Generation in Distribution System and microgrid | |
| Cui et al. | Multi-objective optimal allocation of distributed generation in smart grid | |
| CN103490428A (en) | Method and system for allocation of reactive compensation capacity of microgrid | |
| CN109524988A (en) | A kind of wind-powered electricity generation based on total active power trend prediction collects station voltage control method | |
| CN116345466A (en) | Two-stage power flow optimization method of active power distribution network considering distribution network reconstruction | |
| CN102904266A (en) | Method for determining inactive compensation capacity network adaptability of wind power plant | |
| CN104362650A (en) | Electric power system reactive power optimization method considering cost factor | |
| CN108767891A (en) | Feeder line photovoltaic is pressed to dissolve model in medium voltage distribution network single | |
| CN105226649B (en) | One kind predicting improved provincial power network power generation dispatching optimization method based on bus load | |
| CN109617086A (en) | A kind of power grid dynamic passive compensation equipment coordination control system and method | |
| Yang et al. | A multi-period scheduling strategy for ADN considering the reactive power adjustment ability of DES | |
| Wang et al. | Reactive power planning in distribution network considering the consumption capacity of distributed generation | |
| CN109802448A (en) | A kind of renewable energy maximum consumption capacity analysis calculation method | |
| Xia et al. | Assessment of the carrying capacity of distributed generation in distribution network considering soft open point | |
| CN109120011B (en) | distributed power distribution network congestion scheduling method considering distributed power sources | |
| CN106532729A (en) | Province-regional coordination control method for 220kV collection substation high-voltage bus voltage |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190716 |
|
| WD01 | Invention patent application deemed withdrawn after publication |