CN108470239B - Active power distribution network multi-target layered planning method considering demand side management and energy storage - Google Patents
Active power distribution network multi-target layered planning method considering demand side management and energy storage Download PDFInfo
- Publication number
- CN108470239B CN108470239B CN201810170942.2A CN201810170942A CN108470239B CN 108470239 B CN108470239 B CN 108470239B CN 201810170942 A CN201810170942 A CN 201810170942A CN 108470239 B CN108470239 B CN 108470239B
- Authority
- CN
- China
- Prior art keywords
- energy storage
- distribution network
- power supply
- layer
- reliability
- 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.)
- Active
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 71
- 238000004146 energy storage Methods 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000007726 management method Methods 0.000 title claims abstract description 29
- 230000008569 process Effects 0.000 claims abstract description 13
- 230000004044 response Effects 0.000 claims description 41
- 238000010248 power generation Methods 0.000 claims description 32
- 230000005611 electricity Effects 0.000 claims description 14
- 238000012423 maintenance Methods 0.000 claims description 11
- 238000005520 cutting process Methods 0.000 claims description 8
- 230000002829 reductive effect Effects 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 239000002245 particle Substances 0.000 claims description 4
- 238000002271 resection Methods 0.000 claims description 3
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 230000006870 function Effects 0.000 description 12
- 238000009434 installation Methods 0.000 description 11
- 238000005457 optimization Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000007599 discharging Methods 0.000 description 5
- 230000035699 permeability Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 229910052744 lithium Inorganic materials 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Data Mining & Analysis (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Pure & Applied Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- Tourism & Hospitality (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Primary Health Care (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Probability & Statistics with Applications (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention relates to a multi-target hierarchical planning method for an active power distribution network, which takes demand side management and energy storage into consideration. A double-layer model is established: the upper-layer planning model takes various economic costs of a power distribution company and the reliability and stability of a power distribution network as upper-layer targets in a planning process, considers investment and benefits of demand side management and energy storage equipment, and determines a target function and constraint conditions; the lower layer planning model takes the minimum cut-off amount of the distributed power supply in the planning process as an objective function, considers the maximum effective utilization of the distributed power supply, reduces the output of the distributed power supply through measures of wind abandoning, light abandoning and the like, and adjusts a transformer tap, reactive power compensation equipment and the like to perform active control. The invention can not only improve the reliability of power supply and the stability of distribution network operation, but also delay the upgrade of the power grid.
Description
Technical Field
The invention relates to a multi-target hierarchical planning method for an active power distribution network, which takes demand side management and energy storage into consideration.
Background
With the continuous improvement of the ratio of a multi-type Distributed Generation (DG) to be connected to a power distribution network and the continuous increase of the application of controllable equipment such as flexible loads, the traditional power distribution network planning cannot meet the requirement of high-efficiency utilization of the DG with high permeability under the low-carbon economic background. Therefore, research on active power distribution network planning considering demand side management and energy storage is urgently needed, active consumption and multi-stage coordinated utilization of various renewable energy sources are achieved, uncertain factors such as flexible loads are reasonably regulated and controlled, and network topology is changed to improve system reliability and stability.
At present, in the research of an active power distribution network, random uncertainty of distributed power supplies such as fans and photovoltaic power supplies is generally considered, and then constant volume and site selection are carried out on the distributed power supplies. The ADN research also generally adopts double-layer model planning, and the upper layer comprehensively considers the realization of multiple targets from the investment level; the lower layer considers the comprehensive application of DG from the aspect of operation simulation. The demand side management is a new concept that the power of the United states in the last 80 th century is saved on the demand side, and energy is reasonably utilized as replaceable resources on the supply side, so that the consumption of the energy is reduced and the utilization rate is improved by saving the power of the terminal, and the aim of relieving the power shortage of users is fulfilled. The Energy Storage Battery (Battery Energy Storage System, BESS) is also an important interactive resource under the ADN framework, and can optimize the Energy utilization mode and improve the power supply reliability.
Disclosure of Invention
The invention aims to provide an active power distribution network multi-target layered planning method considering demand side management and energy storage, which can improve the power supply reliability and the stability of distribution network operation and delay the upgrading of a power grid.
In order to achieve the purpose, the technical scheme of the invention is as follows: a multi-target layered planning method for an active power distribution network considering demand side management and energy storage comprises the following steps,
s1, determining an objective function and a constraint condition by taking various economic costs of a power distribution company and the reliability and stability of a power distribution network as upper-layer targets and considering investment and benefits of demand side management and energy storage equipment in the planning process, and establishing an upper-layer planning model;
step S2, taking the minimum expected value of annual cutting amount of the distributed power supply as a target, considering the effective utilization of the distributed power supply, and establishing a lower-layer planning model;
and step S3, solving the upper layer model and the lower layer model based on the improved particle swarm optimization.
In an embodiment of the present invention, the step S1 is specifically implemented as follows:
firstly, modeling is carried out on the power requirement of the electric automobile, firstly, a probability model of the charging starting time and the daily mileage is established, and the probability density function is as follows:
in the formula, σtRepresents the variance value, mu, of the first trip timetThe average value of the first trip time is represented and is two parameters of normal distribution;
further establishing a probability model of the charging power, and assuming that the charging period is TrIf the charging power of the user is Δ E, the charging load can be obtained by the following equation:
Socs=Socr-MdEave/(100Bc)
ΔE=Esoc-Socs
Tr=ΔE/(Prηc)
in the formula, PrIndicating rated charging power of electric vehicle, ηcRepresents the charging efficiency of the electric vehicle, MdFor daily mileage, EaveFor electric energy consumed per stroke, BcFor the capacity of a lithium battery, Socr is the initial electric quantity before the electric automobile goes out, Socs is the initial charging electric quantity after the electric automobile goes back, and Esoc is the expected value of full electric quantity;
then, establishing an upper-layer planning model, taking various economic costs of a power distribution company and the reliability and stability of a power distribution network as upper-layer targets in a planning process, considering investment and benefits of demand side management and energy storage equipment, determining a target function and constraint conditions, and establishing the upper-layer planning model:
wherein f is1Representing an economic objective; f. of2Expressed as reliability target, λreliabilityRepresenting the reliability of the power supply of the distribution network, f3The voltage stability index of the system is represented; cLine、CDG、Csub、Closs、CMA、CS,DSR、CIL、CbThe method comprises the following steps of respectively setting investment and operation maintenance cost of a distribution line, investment and operation maintenance cost of a distributed power supply, operation electricity purchasing cost of a transformer substation, total system loss cost, total annual energy storage battery investment cost, annual energy storage battery operation cost, comprehensive load demand response project management cost, annual electricity selling income reduced after demand response is implemented, compensation cost for interruptible load, and government subsidy cost of unit wind power generation and photovoltaic power generation; lambda [ alpha ]reliabilityReliability of power supply to the distribution network; FVSIkFor resolvability of branch voltage, ΩBIs a set of all branches.
In an embodiment of the present invention, the step S2 is specifically implemented as follows:
establishing a lower-layer planning model, taking the minimum cut-off amount of the distributed power supply in the planning process as an objective function, and considering the maximum effective utilization of the distributed power supply, the lower-layer objective function is expressed as follows:
in the formula, omega1、Ω2Is respectively the set of newly built photovoltaic power generation and newly built wind power generation, omegaSFor a set of scenes, Δ tsFor the cumulative running time of the distribution network year under the scene s,the active power cutting amount of the jth photovoltaic power generation and wind power generation under the scene s is obtained;
the lower model reduces the output of the distributed power supply by measures including wind abandoning and light abandoning, and active control is performed by adjusting a transformer tap and reactive compensation equipment, and the constraint conditions are as follows:
(1) DG output resection constraint
In the formula:respectively representing the upper limit and the lower limit of the output cutting of photovoltaic power generation and wind power generation;
(2) reactive compensation equipment switching amount constraint
In the formula: qCiThe switching quantity of the reactive compensation equipment at the ith node is represented,andrespectively represents the lower limit and the upper limit, omega, of the switching quantity of the reactive compensation equipmentCInstalling a node set for reactive compensation;
(3) transformer tap adjustment range constraints
In the formula: t iskIndicating the position of the tap of the transformer,andrespectively representing the lower and upper limits of the adjustment range of the transformer tap.
Compared with the prior art, the invention has the following beneficial effects: the invention comprehensively considers the economy, reliability AND stability of the power distribution network AND considers the optimization of demand side management AND the energy storage battery on AND. The stability and the reliability of the power distribution network are considered while the maximization of the economic benefits of operators is guaranteed.
Drawings
FIG. 1 is a diagram of a two-level programming relationship.
FIG. 2 is a flow chart of model solution.
Fig. 3 is a modified 33-node diagram.
Fig. 4 is a four season industrial load demand response, where (a) is a spring industrial load demand response where (b) is a summer industrial load demand response, where (c) is a fall industrial load demand response where (d) is a winter industrial load demand response.
Fig. 5 is a four-season resident load demand response in which (a) is a spring resident load demand response in which (b) is a summer resident load demand response in which (c) is an autumn resident load demand response in which (d) is a winter resident load demand response.
Fig. 6 is a four-season commercial load demand response, where (a) is a spring commercial load demand response, where (b) is a summer commercial load demand response, where (c) is an autumn commercial load demand response, where (d) is a winter commercial load demand response.
Fig. 7 is a Pareto result without regard to demand response and energy storage planning.
Fig. 8 is Pareto results considering demand response and energy storage planning.
Fig. 9 is a three-dimensional Pareto result considering demand response and energy storage planning.
Fig. 10 is a wind energy storage optimization diagram of the node 6.
Fig. 11 is a wind energy storage optimization diagram of the node 27.
Fig. 12 is a graph of node 8 optical energy storage optimization.
Detailed Description
The technical scheme of the invention is specifically explained in the following by combining the attached drawings 1-12.
The invention relates to a multi-target layered planning method for an active power distribution network, which takes demand side management and energy storage into consideration and comprises the following steps,
s1, determining an objective function and a constraint condition by taking various economic costs of a power distribution company and the reliability and stability of a power distribution network as upper-layer targets and considering investment and benefits of demand side management and energy storage equipment in the planning process, and establishing an upper-layer planning model;
step S2, taking the minimum expected value of annual cutting amount of the distributed power supply as a target, considering the effective utilization of the distributed power supply, and establishing a lower-layer planning model;
and step S3, solving the upper layer model and the lower layer model based on the improved particle swarm optimization.
The following is a specific implementation of the present invention.
As shown in fig. 1 and 2, the multi-objective hierarchical planning method for an active power distribution network, which takes demand side management and energy storage into consideration, is specifically implemented as follows:
(1) and performing time sequence load modeling on the charging and discharging of the electric automobile, and further predicting and obtaining the charging load distribution of the electric automobile 24 hours a day by predicting the probability distribution of the charging starting time and the daily mileage of the electric automobile.
(2) Considering the influence of demand side Management on active power distribution network planning, the invention mainly considers the demand side Interruptible Load Management (ILM) to model the ILM.
(3) BESS has the functions of peak clipping and valley filling and voltage quality improvement, and the BESS charge-discharge model is established in the aspects of residual charge level (SOC) and the like.
(4) And determining a target function and constraint conditions by taking various economic costs of a power distribution company and the reliability and stability of a power distribution network as upper-layer targets and considering investment and benefits of demand side management and energy storage equipment in the planning process, and establishing an upper-layer planning model.
(5) And establishing a lower-layer planning model by taking the minimum expected value of the annual cutting-off quantity of the distributed power supply as a target and considering the effective utilization of the distributed power supply.
(6) And solving the upper layer model and the lower layer model based on an improved particle swarm algorithm.
The active power distribution network multi-target hierarchical planning considering demand side management and energy storage is described as follows:
firstly, modeling is carried out on the power requirement of the electric automobile, firstly, a probability model of the charging starting time and the daily mileage is established, and the probability density function is as follows:
in the formula, σtRepresents the variance value, mu, of the first trip timetThe average value of the first trip time is represented and is two parameters of normal distribution;
further establishing a probability model of the charging power, and assuming that the charging period is TrIf the charging power of the user is Δ E, the charging load can be obtained by the following equation:
Socs=Socr-MdEave/(100Bc)
ΔE=Esoc-Socs
Tr=ΔE/(Prηc)
in the formula, PrIndicating the rated charging power of the EV (electric vehicle), ηcM represents the charging efficiency of this EVdFor daily mileage, EaveFor electric energy consumed per stroke, BcFor the capacity of a lithium battery, Socr is the initial electric quantity before the EV travels, Socs is the initial charging electric quantity after the EV travels back, and Esoc is the expected value of full electric quantity;
then, establishing an upper-layer planning model, taking various economic costs of a power distribution company and the reliability and stability of a power distribution network as upper-layer targets in a planning process, considering investment and benefits of demand side management and energy storage equipment, determining a target function and constraint conditions, and establishing the upper-layer planning model:
wherein f is1Representing an economic objective; f. of2Of the representationFor reliability purposes, λreliabilityRepresenting the reliability of the power supply of the distribution network, f3The voltage stability index of the system is represented; cLine、CDG、Csub、Closs、CMA、CS,DSR、CIL、CbThe method comprises the following steps of respectively setting investment and operation maintenance cost of a distribution line, investment and operation maintenance cost of a distributed power supply, operation electricity purchasing cost of a transformer substation, total system loss cost, total annual energy storage battery investment cost, annual energy storage battery operation cost, comprehensive load demand response project management cost, annual electricity selling income reduced after demand response is implemented, compensation cost for interruptible load, and government subsidy cost of unit wind power generation and photovoltaic power generation; lambda [ alpha ]reliabilityReliability of power supply to the distribution network; FVSIkFor resolvability of branch voltage, ΩBIs a set of all branches.
(1) Economic objective f1
1) Distribution line investment and operation maintenance cost CLine
In the formula:investment cost for distribution lines;the maintenance cost is equal to the annual operation cost of the line; omegaL1、ΩL2And ΩLRespectively setting a newly-built circuit, an upgraded circuit and a set of all circuits;the costs required for newly building a line, upgrading the line and the unit length of all the lines are respectively; lijIs the line ij length; r is the discount rate; t isLineThe investment recovery period of the distribution line is set; gamma rayMMaintaining the cost rate for the line operation;
2) investment and operation maintenance cost C of distributed power supplyDG
In the formula:equal annual value investment cost for a distributed power supply of a power distribution network;operating and maintaining costs for the distributed power supply at an annual value; cPVGAnd CWTGInstallation costs per unit capacity for photovoltaic power generation (PVG) and wind power generation (WTG), respectively; omega1、Ω2The method comprises the steps of respectively integrating newly built photovoltaic power generation and newly built wind power generation;respectively the installation capacity of each photovoltaic and fan of the jth; t isDGA recovery period for the investment of the distributed power supply;andrespectively unit operation and maintenance costs of photovoltaic power generation and wind power generation; omegaSIs a collection of scenes; Δ tsAccumulating the running time of the distribution network in a scene s;the active output of the jth photovoltaic power generation under the scene s is obtained;the active power output of the jth wind power generation under the scene s is obtained;
3) cost C for purchasing electricity during operation of transformer substationsub
In the formula: ceThe unit electricity purchasing cost; n is the total number of the load nodes of the power distribution network;the active load power of the distribution network node j under the scene s is shown. Due to the distributed power supply in the system, the cost required by the transformer substation for supplying power to the system is that all load power and the power which is subtracted from the distributed power supply are supplied.
4) Total loss cost of system Closs
In the formula: omegaSubA substation node set in the system is obtained; i isij,sAnd Ij,sThe magnitudes of the currents flowing through line ij and node j in the s-th scenario, respectively. ZSubIs the equivalent resistance of the substation, rijIs the resistance of line ij.
5) Investment cost and profit for energy storage equipment
The upper typeIn order to realize the annual investment total cost of the energy storage battery,the unit investment cost of the wind energy storage equipment of the distributed wind power node,unit investment cost, omega, of optical energy storage devices for distributed photovoltaic nodesWTGAnd ΩPVGRespectively a distributed wind power generator node and a photovoltaic node,investing decision variables at this point for energy storage whenThe energy storage device is invested in the distributed power node,it means that no energy storage equipment is invested at this point.
The annual running cost of the energy storage battery is saved;andthe unit cost of wind storage equipment and light storage charging in unit time respectively;andrespectively obtaining unit benefits of wind storage equipment and light storage and discharge electricity in unit time; t isi chargeAnd Ti unchargeRespectively representing the charging time and the discharging time of the node i.
6) Demand response interrupt load cost
Above formula CMAManagement of costs for integrated load demand response projects, CS,DSRReduced annual electricity sales revenue after implementing demand response for the s-th scenario; ccurA compensation fee for interruptible loads; cmaThe cost is managed for the demand response item per unit power. Cp,s,tThe electricity selling price provided by the power distribution network company in the tth scene time is provided; the original total load of the network of the power grid company in the t-th period of the s scene is large; pCur,s,tActive power of active interrupt load for a scene s time period t, PIncre,s,tAnd active power actively increasing the load in the second period in the s-th scene.
7) Subsidy cost C for new energy power generationb
In the formula: c. CbThe method is the government subsidy cost of wind power generation and photovoltaic power generation.
(2) Reliability target f2
f2The middle ASAI represents the average power supply availability of the power distribution network, the larger the ASAI is, the higher the reliability of the power distribution network is, and the ASAI is converted into f2The smaller the reliability target, the higher the reliability of the distribution network.
The specific calculation steps are as follows:
1) inputting various parameters of the power distribution network;
2) dividing areas by taking the switch element as a boundary to form element groups and load areas, and numbering each area for each parameter;
3) the component group is equivalent upwards from the maximum level to the level 1, and the final equivalent is the simplest radiation network;
4) calculating the reliability of the element group on the main feeder;
5) the element group is equivalent downwards, and the reliability parameters of the maximum level load area can be calculated in sequence;
6) calculating a load point reliability parameter;
7) and then calculates the system reliability index ASAI.
The specific calculation formula is as follows:
ASAI=(N×8760-T)/(N×8760)
in the formula: n is the total number of power supply users, and T is the total power failure time of the users.
(3) Stability goal f3
f3Medium FVSI refers to a rapid voltage stabilization index, and the smaller the FVSI, the more stable the power distribution network is. FVSIijThe solvability of the branch voltage is represented, the value of which may represent the voltage stability of the branch. Firstly, the FVSI of each branch is calculatedijThen, the largest of them is FVSI, and the calculation formula is as follows:
FSVIij=4Z2Qj/Vi 2X
FSVI=max(FSVIij)
in the formula: z and X represent the impedance and reactance of the branch, respectively; qjRepresenting reactive power flow of the tail end branch; viRepresenting the branch head voltage magnitude.
In order to accelerate the calculation speed, a weight method is introduced to process a reliability target and a stability target, three targets are converted into double targets to be processed, and the specific calculation formula is as follows:
F'=βF2+(1-β)F3
in the formula: f' represents the comprehensive index after the weight is introduced; β is a weight, and is taken herein to be 0.98.
Considering that the distributed power supply is connected into a system and then impacts on the original electrical characteristics of the power distribution network, such as power flow distribution, voltage level, short-circuit capacity and the like, the permeability is required to meet a certain range. The specific constraints are expressed as follows:
(1) distribution network flow constraint
In the formula: x is the number ofijIs the reactance of line ij, ΩBThe node sets are all node sets of the system, (j) is a branch end node set taking j as a head end node, and pi (j) is a branch head end node set taking j as a tail end node; qjk,s、Qij,sReactive power for branches jk and ij; respectively providing reactive power for a transformer substation node, a wind power generation node, a photovoltaic node and a load node; u shapej,sIs the node voltage.
(2) Probability constraint of upper and lower limits of node voltage
P{Umin≤U≤Umax}=k1/N≥β1
In the formula: u shapeminAnd UmaxThe lower limit and the upper limit of the node voltage are respectively; k is a radical of1The number of scenes satisfying the voltage upper and lower limit constraint in all scenes β1Is the confidence level of the node voltage constraint.
(3) Branch power probability constraint
P{Pij≤Pij,max}=k2/N≥β2
In the formula: pij,maxThe upper power limit allowed for branch ij; k is a radical of2For the number of scenarios satisfying the branch power constraint among all scenarios β2Confidence level for branch power constraints。
(4) Inhibiting reverse power probabilistic constraints
P{PΣDG≤PΣL}=kB/N≥βB
In the formula: pΣDGAnd PΣLRespectively the total active power of the distributed power supply and the total active power of the load; k is a radical ofBFor the number of scenarios satisfying the power back-off prohibition constraint among all scenarios βBTo inhibit confidence levels of the reverse power constraint.
(5) Distributed power source installation capacity constraints
In the formula:andPVG and WTG total installation capacity respectively; sigma is the maximum permeability of new energy access;the sum of the maximum active load of the distribution network;andrespectively obtaining the maximum installation capacity of PVG and WTG of the grid-connected node i to be selected; pi PVGAnd Pi WTGThe installation capacities of PVG and WTG, respectively, for the candidate installation node i.
(6) Energy storage charge and discharge restraint
Above formula Ti chargeAnd Ti unchargeRespectively representing the charging time and the discharging time of the node i;andfor the charging power and the discharging power of the energy storage device per unit time,andis the maximum charging power per unit time of the energy storage device,the total charge capacity of the energy storage device.
(7) Demand response management constraint
PTL,s,t=PTLO,s,t-PTLI,s,t
In the formula: pTLO,s,tIs PTLI,s,tLoad transfer-out and load transfer-in for the tth scene period respectively, respectively transferring the upper limit value and the lower limit value of the load transferring response coefficient and the load transferring response coefficient for the t-th time period of the s-th scene;and (4) loading interruption power upper and lower limit values for the t-th time period of the s-th scene.
(8) Radial confinement
Firstly, generating an undirected graph by using a minimum spanning tree algorithm, and then generating a directed graph based on a Kruskal idea.
(9) Connectivity constraints
And (4) solving an adjacency matrix and a reachability matrix of the graph, and then judging the connectivity of the graph through the reachability matrix of the graph.
Then, a lower-layer planning model is established, the minimum cut-off amount of the distributed power supply in the planning process is taken as an objective function, and the maximum effective utilization lower-layer objective function of the distributed power supply is considered to be expressed as follows:
The lower layer model reduces the output of the distributed power supply by measures such as abandoning wind and abandoning light, and carries out active control by adjusting a transformer tap, reactive compensation equipment and the like, and the constraint conditions are as follows:
(1) DG output resection constraint
In the formula:respectively representing the upper limit and the lower limit of the output cutting of photovoltaic power generation and wind power generation;
(2) reactive compensation equipment switching amount constraint
In the formula: qCiThe switching quantity of the reactive compensation equipment at the ith node is represented,andrespectively represents the lower limit and the upper limit, omega, of the switching quantity of the reactive compensation equipmentCInstalling a node set for reactive compensation;
(3) transformer tap adjustment range constraints
In the formula: t iskIndicating the position of the tap of the transformer,andrespectively representing the lower and upper limits of the adjustment range of the transformer tap.
The following is illustrated by specific examples:
the simulation system of the invention adopts an improved 33-node system, and the topological diagram of the system is shown in figure 3.
For example, the improved 33-node system has 61 lines and 39 nodes, where the nodes 34-39 are newly added load nodes, and the branches 38-61 are lines to be newly added. The installation reference capacity of the fan and the photovoltaic power supply in the system is 100kW, and the maximum allowable permeability is 50%. Setting the nodes 3, 6, 16 and 27 as the nodes for accessing the fan, wherein the upper limit of the installation number of the nodes is respectively 20, 18, 28 and 18; the number of the nodes to be connected into the photovoltaic power supply is 8, 10, 28 and 30, and the upper limits of the installation number of the nodes are 10, 20, 8 and 10 respectively. The line operation and maintenance cost rate and the discount rate are respectively set to be 3% and 0.1. As shown in tables 1 and 2, the fixed investment recovery periods for the line and DG were 20 years and 10 years, respectively.
TABLE 1 line parameters
TABLE 2 DG parameters
In this example, the nodes 4, 7, 10, 13, 16, 19, 22, 25, 28, 31, 34, and 37 are assumed to be the residential load nodes; nodes 2, 5, 8, 11, 14, 17, 20, 23, 26, 29, 32, 35, 38 are commercial load nodes; nodes 3, 6, 9, 12, 15, 18, 21, 24, 27, 30, 33, 36, 39 are industrial loads.
The improved PSO algorithm parameters are set to: and (4) iterating for 50 times, wherein the population size is 80, the initial value and the final value of the inertia weight are respectively 0.8 and 0.4, the initial value and the final value of one learning factor are respectively 2.5 and 0.5, and the initial value and the final value of the other learning factor are respectively 0.5 and 2.5.
Through the solution simulation of the model, demand response time sequence curves before and after optimization can be obtained, and the demand response time sequence curves are shown in fig. 4, 5 and 6. Four-season demand response curves for industrial, residential and commercial loads, respectively.
It can be seen from the curve results of the above graph that, through the demand response, after the user responds to the electricity price strategy, the interruption or reduction of the electricity consumption is allowed at some high electricity price moments, and the overall power balance of the system is improved; and the power consumption requirement is increased at certain low-price stages, the resource distribution of a power distribution system can be improved, and the phenomena of wind abandonment and power abandonment are prevented. If the price incentive strategy is adjusted, the demand response distribution curve changes along with the strategy, and reasonable control of part of loads in the active power distribution network can be realized.
The same system is respectively planned, and planning results without considering the demand response and the energy storage equipment and with considering the demand response and the energy storage equipment are compared. Fig. 7 and fig. 8 can be obtained through simulation, and it can be seen that the Pareto optimal boundary considering demand response and energy storage is obviously better than the Pareto optimal boundary not considering demand response and energy storage, and all solutions are non-dominant solutions relative to the situation not considered, that is, all the former solutions answer economic cost and comprehensive indexes smaller than the latter.
Meanwhile, fig. 9 also shows an optimal solution set in a three-dimensional space, i.e., a comprehensive index is decomposed into reliability and stability indexes.
The pareto optimal frontier shown in the upper graph proves the effectiveness of the solution of the method of the invention, the set of machines for obtaining the optimal solution for planning can be shown in table 3,
TABLE 3 planning scenario results
The results in the table above show that scheme 5, although there are more upgraded lines and newly built lines, leads to lower overall cost through energy storage and demand response management.
To further illustrate the optimization effect of energy storage on the distributed power source, fig. 10 to 12 show the comprehensive curves of the energy storage device power in the planning scheme 5. The dotted line is the output power of the distributed power supply after the energy storage device is added, the star curve is the output power before the energy storage is not added, and the solid line is the charge and discharge power of the energy storage device. Through the charging and discharging of the energy storage equipment, the output of the distributed power supply is improved to a certain extent, the input cost of the newly-built distributed power supply is saved, the system loss cost is reduced in a specific period, and the planning is influenced significantly.
The above are preferred embodiments of the present invention, and all changes made according to the technical scheme of the present invention that produce functional effects do not exceed the scope of the technical scheme of the present invention belong to the protection scope of the present invention.
Claims (1)
1. A multi-target layered planning method for an active power distribution network considering demand side management and energy storage is characterized by comprising the following steps,
s1, determining an objective function and a constraint condition by taking various economic costs of a power distribution company and the reliability and stability of a power distribution network as upper-layer targets and considering investment and benefits of demand side management and energy storage equipment in the planning process, and establishing an upper-layer planning model;
step S2, taking the minimum expected value of annual cutting amount of the distributed power supply as a target, considering the effective utilization of the distributed power supply, and establishing a lower-layer planning model;
s3, solving an upper layer model and a lower layer model based on an improved particle swarm algorithm;
the step S1 is specifically implemented as follows:
establishing an upper-layer planning model, taking various economic costs of a power distribution company and the reliability and stability of a power distribution network into comprehensive consideration in the planning process as upper-layer targets, considering investment and benefits of demand side management and energy storage equipment, determining a target function and constraint conditions, and establishing the upper-layer planning model:
wherein f is1Representing an economic objective; f. of2Expressed as reliability target, λreliabilityRepresenting the reliability of the power supply of the distribution network, f3The voltage stability index of the system is represented; cLine、CDG、Csub、Closs、CMA、CS,DSR、CIL、CbThe method comprises the following steps of respectively setting investment and operation maintenance cost of a distribution line, investment and operation maintenance cost of a distributed power supply, operation electricity purchasing cost of a transformer substation, total system loss cost, total annual energy storage battery investment cost, annual energy storage battery operation cost, comprehensive load demand response project management cost, annual electricity selling income reduced after demand response is implemented, compensation cost for interruptible load, and government subsidy cost of unit wind power generation and photovoltaic power generation; lambda [ alpha ]reliabilityReliability of power supply to the distribution network; FVSIkFor resolvability of branch voltage, ΩBA set composed of all branches;
the step S2 is specifically implemented as follows:
establishing a lower-layer planning model, taking the minimum cut-off amount of the distributed power supply in the planning process as an objective function, and considering the maximum effective utilization of the distributed power supply, the lower-layer objective function is expressed as follows:
in the formula, omega1、Ω2Is respectively the set of newly built photovoltaic power generation and newly built wind power generation, omegaSFor a set of scenes, Δ tsFor the cumulative running time of the distribution network year under the scene s,the active power cutting amount of the jth photovoltaic power generation and wind power generation under the scene s is obtained;
the lower model reduces the output of the distributed power supply by measures including wind abandoning and light abandoning, and active control is performed by adjusting a transformer tap and reactive compensation equipment, and the constraint conditions are as follows:
(1) DG output resection constraint
In the formula:respectively representing the upper limit and the lower limit of the output cutting of photovoltaic power generation and wind power generation;
(2) reactive compensation equipment switching amount constraint
In the formula: qCiThe switching quantity of the reactive compensation equipment at the ith node is represented,andrespectively represents the lower limit and the upper limit, omega, of the switching quantity of the reactive compensation equipmentCInstalling a node set for reactive compensation;
(3) transformer tap adjustment range constraints
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810170942.2A CN108470239B (en) | 2018-03-01 | 2018-03-01 | Active power distribution network multi-target layered planning method considering demand side management and energy storage |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810170942.2A CN108470239B (en) | 2018-03-01 | 2018-03-01 | Active power distribution network multi-target layered planning method considering demand side management and energy storage |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108470239A CN108470239A (en) | 2018-08-31 |
CN108470239B true CN108470239B (en) | 2020-09-04 |
Family
ID=63264090
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810170942.2A Active CN108470239B (en) | 2018-03-01 | 2018-03-01 | Active power distribution network multi-target layered planning method considering demand side management and energy storage |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108470239B (en) |
Families Citing this family (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109214597B (en) * | 2018-10-25 | 2021-08-03 | 浙江工业大学 | Optimization planning method for micro-grid power and operation reserve capacity |
CN109390940A (en) * | 2018-11-30 | 2019-02-26 | 国家电网公司西南分部 | A kind of sending end electric network source planing method considering demand response and comprehensive benefit |
CN109558984A (en) * | 2018-12-05 | 2019-04-02 | 中南大学 | A kind of active distribution network distributed energy storage system Bi-objective hybrid optimization planing method |
CN109711706B (en) * | 2018-12-21 | 2023-04-07 | 国网江西省电力有限公司经济技术研究院 | Active power distribution network transformer substation planning method considering distributed power sources and demand response |
CN109871981B (en) * | 2019-01-10 | 2021-07-13 | 国家电网有限公司 | Load characteristic prediction method considering distributed power supply and electric vehicle influence |
CN109871989A (en) * | 2019-01-29 | 2019-06-11 | 国网山西省电力公司吕梁供电公司 | A kind of power distribution network hierarchical reconfiguration planning method containing distributed generation resource |
CN109934487A (en) * | 2019-03-11 | 2019-06-25 | 国网福建省电力有限公司 | A kind of active distribution network coordinated planning method considering multiagent interest game |
CN109948954B (en) * | 2019-04-04 | 2021-07-13 | 华北电力大学 | Power distribution network bidirectional blocking scheduling method for distributed resources of power system |
CN110119886B (en) * | 2019-04-18 | 2022-11-25 | 深圳供电局有限公司 | Active distribution network dynamic planning method |
CN110210659B (en) * | 2019-05-24 | 2021-04-02 | 清华大学 | Power distribution network planning method considering reliability constraint |
CN110504708A (en) * | 2019-08-09 | 2019-11-26 | 国家电网有限公司 | The power distribution network multiple target collaborative planning method of meter and charging station and distributed generation resource |
CN110502814B (en) * | 2019-08-09 | 2023-07-21 | 国家电网有限公司 | Active power distribution network multi-target planning method considering energy storage and load management technology |
CN110429653B (en) * | 2019-08-28 | 2020-11-17 | 国网河北省电力有限公司邢台供电分公司 | Rural power grid distributed photovoltaic absorption method considering energy storage and DR (digital radiography) and terminal equipment |
CN111079972A (en) * | 2019-11-04 | 2020-04-28 | 深圳供电局有限公司 | Method, device and medium for planning reliability of active power distribution network |
CN110895638B (en) * | 2019-11-22 | 2022-12-06 | 国网福建省电力有限公司 | Active power distribution network model establishment method considering electric vehicle charging station site selection and volume fixing |
CN111125887B (en) * | 2019-12-04 | 2022-03-22 | 广东电网有限责任公司 | Resource optimization configuration model establishment method for emergency recovery of power distribution network |
CN111490554B (en) * | 2020-04-16 | 2023-07-04 | 国网江苏省电力有限公司淮安供电分公司 | Multi-objective optimal configuration method for distributed photovoltaic-energy storage system |
CN111404206B (en) * | 2020-04-21 | 2021-07-09 | 清华大学 | Wind-solar energy storage power generation system capacity double-layer planning method considering investment return constraint |
CN111652411B (en) * | 2020-05-15 | 2022-07-15 | 三峡大学 | Distributed power supply double-layer planning method |
CN112039069B (en) * | 2020-09-04 | 2022-02-22 | 国网山东省电力公司济宁供电公司 | Double-layer collaborative planning method and system for power distribution network energy storage and flexible switch |
CN112182907B (en) * | 2020-10-19 | 2022-12-27 | 贵州电网有限责任公司 | Reliability constraint-based planning method for energy storage device of electric-gas coupling system |
CN112380694B (en) * | 2020-11-13 | 2022-10-18 | 华北电力大学(保定) | Power distribution network optimization planning method based on differential reliability requirements |
CN112448403B (en) * | 2020-11-16 | 2022-05-17 | 国网四川省电力公司经济技术研究院 | Decoupling configuration method for energy storage of power distribution network |
CN112821451B (en) * | 2021-01-11 | 2023-05-09 | 国网福建省电力有限公司泉州供电公司 | Urban power distribution network photovoltaic access response method based on demand side management and energy storage |
CN113572843B (en) * | 2021-07-26 | 2023-06-02 | 成都秦川物联网科技股份有限公司 | Gas data off-network operation method and system based on energy metering |
US11898895B2 (en) | 2021-07-26 | 2024-02-13 | Chengdu Qinchuan Iot Technology Co., Ltd. | Methods and systems for natural gas data computation outside gas internet of things based on energy measuring |
CN113988392B (en) * | 2021-10-19 | 2024-05-28 | 华北电力大学(保定) | Micro-grid optimization planning method considering reliability demand response |
CN115333950B (en) * | 2022-05-12 | 2024-05-03 | 四川大学 | Coordinated planning method for information physical active power distribution system considering hybrid communication networking |
CN114881377B (en) * | 2022-07-11 | 2024-01-30 | 国网江苏省电力有限公司苏州供电分公司 | Communication base station energy storage optimization method, system and device considering 5G Internet of vehicles |
CN115147016B (en) * | 2022-09-01 | 2023-01-24 | 国网浙江省电力有限公司宁波供电公司 | Distributed power supply and demand side resource combined planning method considering low-carbon benefits |
CN116780529A (en) * | 2023-06-30 | 2023-09-19 | 国网北京市电力公司 | Power distribution network fault recovery method, device, equipment and medium |
CN117993740A (en) * | 2024-04-03 | 2024-05-07 | 国网山西省电力公司营销服务中心 | Multi-element power distribution network configuration method considering N-1 fault load loss cost |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106815657A (en) * | 2017-01-05 | 2017-06-09 | 国网福建省电力有限公司 | A kind of power distribution network bi-level programming method for considering timing and reliability |
CN107301470A (en) * | 2017-05-24 | 2017-10-27 | 天津大学 | A kind of power distribution network Expansion Planning stores up the dual blank-holder of addressing constant volume with light |
CN107565573A (en) * | 2017-09-28 | 2018-01-09 | 华中科技大学 | A kind of active distribution network operational plan optimization method for considering Demand-side Respondence to the Price of Electric Power |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8862279B2 (en) * | 2011-09-28 | 2014-10-14 | Causam Energy, Inc. | Systems and methods for optimizing microgrid power generation and management with predictive modeling |
-
2018
- 2018-03-01 CN CN201810170942.2A patent/CN108470239B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106815657A (en) * | 2017-01-05 | 2017-06-09 | 国网福建省电力有限公司 | A kind of power distribution network bi-level programming method for considering timing and reliability |
CN107301470A (en) * | 2017-05-24 | 2017-10-27 | 天津大学 | A kind of power distribution network Expansion Planning stores up the dual blank-holder of addressing constant volume with light |
CN107565573A (en) * | 2017-09-28 | 2018-01-09 | 华中科技大学 | A kind of active distribution network operational plan optimization method for considering Demand-side Respondence to the Price of Electric Power |
Non-Patent Citations (1)
Title |
---|
考虑不同类型DG和负荷建模的主动配电网协同规划;高红均 等;《中国电机工程学报》;20160920;第36卷(第18期);4911页-4920页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108470239A (en) | 2018-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108470239B (en) | Active power distribution network multi-target layered planning method considering demand side management and energy storage | |
CN110071505B (en) | Power transmission network extension and energy storage configuration combined planning method with large-scale wind power access | |
CN110502814B (en) | Active power distribution network multi-target planning method considering energy storage and load management technology | |
CN107887903B (en) | Micro-grid robust optimization scheduling method considering element frequency characteristics | |
Yang et al. | Optimal sizing and placement of energy storage system in power grids: A state-of-the-art one-stop handbook | |
Elsied et al. | Optimal economic and environment operation of micro-grid power systems | |
CN109325608B (en) | Distributed power supply optimal configuration method considering energy storage and considering photovoltaic randomness | |
CN103986190B (en) | Based on the wind-solar-storage joint electricity generation system smooth control method of generated output curve | |
CN108446796A (en) | Consider net-source-lotus coordinated planning method of electric automobile load demand response | |
Lotfi et al. | Optimal design of a hybrid solar-wind-diesel power system for rural electrification using imperialist competitive algorithm | |
CN107994595A (en) | A kind of system of peak load shifting control method and system and the application control method | |
CN109948849B (en) | Power distribution network frame planning method considering energy storage access | |
CN110581571A (en) | dynamic optimization scheduling method for active power distribution network | |
CN109830976B (en) | Elastic operation regulation and control method for alternating current/direct current hybrid power distribution network | |
Hossain et al. | Design a novel controller for stability analysis of microgrid by managing controllable load using load shaving and load shifting techniques; and optimizing cost analysis for energy storage system | |
Zaree et al. | An MILP formulation for centralized energy management strategy of microgrids | |
CN112072711A (en) | Power distribution network flexibility optimization scheduling method based on dynamic priority | |
Bruno et al. | Predictive dispatch across time of hybrid isolated power systems | |
CN111009914A (en) | Active power distribution network-oriented energy storage device location and volume determination method | |
Pippi et al. | Assessing the provision of ancillary services considering BES capacity degradation | |
Chen et al. | Cost-benefit analysis of distributed energy storage in distribution grids with renewables | |
CN104836242A (en) | Reactive power optimization method for regional power grid | |
Kasturi et al. | Optimal PV & BES units integration to enhance power distribution network Performance | |
Koutroulis et al. | Optimal design and economic evaluation of a battery energy storage system for the maximization of the energy generated by wind farms in isolated electric grids | |
CN114398777A (en) | Power system flexibility resource allocation method based on Bashi game theory |
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 |