CN106998072A - A kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network - Google Patents
A kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network Download PDFInfo
- Publication number
- CN106998072A CN106998072A CN201710337499.9A CN201710337499A CN106998072A CN 106998072 A CN106998072 A CN 106998072A CN 201710337499 A CN201710337499 A CN 201710337499A CN 106998072 A CN106998072 A CN 106998072A
- Authority
- CN
- China
- Prior art keywords
- power
- energy storage
- energy
- distribution network
- type energy
- 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
- 238000004146 energy storage Methods 0.000 title claims abstract description 238
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims description 25
- 238000006243 chemical reaction Methods 0.000 claims description 20
- 230000007704 transition Effects 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 8
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 claims description 7
- 229910001416 lithium ion Inorganic materials 0.000 claims description 7
- 230000005611 electricity Effects 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 238000011217 control strategy Methods 0.000 claims description 3
- 238000012423 maintenance Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000007599 discharging Methods 0.000 claims 2
- 239000006185 dispersion Substances 0.000 claims 1
- 230000009191 jumping Effects 0.000 claims 1
- 238000005457 optimization Methods 0.000 abstract description 11
- 238000007796 conventional method Methods 0.000 abstract description 3
- 230000029087 digestion Effects 0.000 abstract description 3
- 241000196324 Embryophyta Species 0.000 description 10
- 239000000243 solution Substances 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 238000000418 atomic force spectrum Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011234 economic evaluation Methods 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 229910052744 lithium Inorganic materials 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
-
- H02J3/383—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E70/00—Other energy conversion or management systems reducing GHG emissions
- Y02E70/30—Systems combining energy storage with energy generation of non-fossil origin
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network, first obtain light-preserved system participate in distribution network operation it is optimal exert oneself with photovoltaic is actual exert oneself between imbalance power, as the reference power of hybrid energy-storing, imbalance power is decomposed into by a series of intrinsic mode functions using set empirical mode decomposition, the instantaneous frequency time graph of each intrinsic mode function is obtained by recurrence Hilbert transform;According to power-type and the characteristic of energy type energy storage, crossover frequency is at least determined for principle with instantaneous frequency time graph aliasing, the high fdrequency component and low frequency component using power-type energy storage and energy type energy storage respectively to imbalance power are stabilized;Finally, it is considered to efficiency for charge-discharge, state-of-charge and overall life cycle cost, power-type and energy type energy storage rated power and rated capacity economic optimum allocation plan are proposed.The present invention can not only compensate new energy and exert oneself relative to conventional method, improve the digestion capability to new energy, and can be actively engaged in the optimization operation of power distribution network.
Description
Technical field
The invention belongs to participate in the technical field of energy storage of power distribution network optimization operation, and in particular to a kind of towards power distribution network optimization
The mixed energy storage system capacity configuration optimizing method of operation.
Background technology
Clean energy resource generation technology by representative of photovoltaic, is just changing the electric power and energy knot in worldwide deeply
Structure.However, photovoltaic generation there are problems that randomness, high permeability it is grid-connected can to distribution power flow distribution, electricity
Energy quality etc., which is brought, to be had a strong impact on.Therefore, the random wave of storage symmetrical photovoltaic plant power output need to be configured in power distribution network
It is dynamic.
Energy-storage system is engaged with photovoltaic, on the one hand can be effectively reduced intermittence and randomness that photovoltaic is exerted oneself, be made light
The power output of volt tends to be steady, and reduces the impact to power network;On the other hand can to energy carry out across when dispatch, participate in distribution
The optimization operation of net.Energy storage can be divided into energy type energy storage and the major class of power-type energy storage two by its power producing characteristics.With lead-acid accumulator,
Its energy density is big for the energy type energy storage of representative for lithium battery etc., but power density is small, the response time is long, it is adaptable to stabilize high energy
The low-frequency power fluctuations of amount.With ultracapacitor, super conductive magnetic storage energy etc. for the power-type energy storage of representative its power density is big, response
Time is short, but energy density is small, it is adaptable to stabilize the high frequency power fluctuation of low energy.Using having energy type energy storage and power-type concurrently
The mixed energy storage system compensation new energy power output of energy storage advantage turns into current research and the focus of application.
Conventional method stabilizes new energy using mixed energy storage system (Hybrid energy storage system, HESS)
The fluctuation of source power output, does not consider that HESS participates in power distribution network optimization operation.The application of energy-storage system can not only compensate new energy
Source is exerted oneself, and improves the digestion capability to new energy, and can be actively engaged in the optimization operation of power distribution network.
The content of the invention
In view of the above-mentioned problems, the present invention proposes that a kind of mixed energy storage system capacity optimization for optimizing operation towards power distribution network is matched somebody with somebody
Method is put, the rated power and rated capacity of energy-storage system can be significantly reduced, is that HESS capacity is distributed rationally there is provided a kind of new
Method, with certain reference value;Meanwhile, HESS application can not only stabilize the fluctuation that photovoltaic is exerted oneself, and promote photovoltaic energy
Source is dissolved, and can participate in the optimization operation of power distribution network, can effectively improve the economy and reliability of distribution network operation.
Above-mentioned technical purpose is realized, above-mentioned technique effect is reached, the present invention is achieved through the following technical solutions:
A kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network, comprises the following steps:
Set up light-preserved system and participate in power distribution network optimal operation model, solution obtains light-preserved system and participates in distribution network operation most
It is excellent to exert oneself, and using light-preserved system participate in distribution network operation it is optimal exert oneself with photovoltaic is actual exert oneself between imbalance power as
The reference power of mixed energy storage system;
By improving, Hilbert-Huang transform decomposition light-preserved system participation distribution network operation is optimal to exert oneself and photovoltaic plant reality
Border exert oneself between imbalance power, imbalance power is divided into high and low frequency two parts wave component, and be utilized respectively
Power-type energy storage and energy type energy storage compensation high-frequency fluctuation component and low-frequency fluctuation component;
Obtain the mixed energy storage system capacity for optimizing operation towards power distribution network and distribute strategy rationally.
Further, the light-preserved system, which participates in power distribution network optimal operation model, includes object function and constraints,
The object function is:
PPESS.loss(t)=(1- ηch)Pch(t)+(1-ηdis)Pdis(t)+ξ(t)·ΔPESS.loss
In formula:T is the cycle of operation;Δ t is operation step-length;N is system node number, Pi(t) it is the injection of t period node is
Active power;ξ (t) is that energy-storage system changes dimension from the t-1 periods to the charging and discharging state of t periods, and value is 0 or 1;ηch、ηdis
For the efficiency for charge-discharge of energy-storage system;Pch(t)、Pdis(t) it is the charge-discharge electric power of t period energy-storage systems;ΔPESS.lossDuring for t
The loss that section energy-storage system is produced when charging and discharging state switches.
Further, the constraints is:
PPESS(t)=PPV(t)+PESS(t)
PPV(t)=PPV.MPPT(t)
-Pch.max≤PESS(t)≤Pdis.max
S (t+1)=(1- σ) S (t)+PESS(t)Δt
Smin≤S(t)≤Smax
S (0)=S (T)
In formula:I, j are system node number;Ui(t)、Uj(t) it is t periods node i and node j point voltage magnitudes;Gij、BijPoint
Transconductance and mutual susceptance that Wei be between node i and node j;δij(t) it is the phase difference between t periods node i and node j;PKi
(t)、QKi(t) it is respectively that the outlet of t periods feeder line is active and idle;PPESSi(t)、QPESSi(t) it is light-preserved system at t period node is
Active and reactive power;PDi(t)、QDi(t) it is the active and reactive power of load at t period node is;Pij(t) it is t
Line power between node i and node j;The subscript "-" and subscript " _ " of variable represent the upper and lower bound of variable respectively;
PPESSAnd Q (t)PESS(t) be respectively t periods light-preserved system output active and reactive power;PPV(t) it is t period photovoltaic generations
Unit active power of output;PMPPTIt is active according to maximum output during MPPT control strategies for photovoltaic DC/DC current transformers;PESS(t)
For the active output of t period energy-storage systems, it is divided into electric discharge, free time, three kinds of states of charging;PchAnd P (t)dis(t) it is respectively the t periods
The charge-discharge electric power of energy-storage system;ηchAnd ηdisRespectively energy-storage system efficiency for charge-discharge;Pch.maxAnd Pdis.maxRespectively charge and discharge
Electrical power bound;S (t) is the dump energy of t period energy-storage systems;σ is the self-discharge rate of energy-storage system;SminAnd SmaxRespectively
For the bound of energy-storage system dump energy.
Further, it is described solution obtain light-preserved system participate in the optimal of distribution network operation exert oneself, be specially:
Power distribution network optimal operation model is participated in using two-dimension dynamic programming Algorithm for Solving light-preserved system, show that light-preserved system is joined
Optimal with distribution network operation is exerted oneself;Specifically include following steps:
Step 2.1:Light-preserved system participation power distribution network optimal operation model is converted into the operable model of two-dimension dynamic programming,
Including stage t, state vector S (sPt,sQt), decision vector U (uPt,uQt), strategy, state transition equation and target function;
Step 2.2:Initialization, light-preserved system and distribution network system primary data needed for input;
Step 2.3:Determine the permission state set of t-1 period light-preserved systems;
Step 2.4:It is determined that meeting the permission state set of t period light-preserved systems under constraints;
Step 2.5:Calculate correspondence decision variable and target function value;
Step 2.6:Whether judge index function is optimal, if so, then skipping to step 2.7;Otherwise, step 2.8 is skipped to;
Step 2.7:Preserve current state variable, decision variable and object function;
Step 2.8:Judge whether the permission state set of traversal t-1 period light-preserved systems, if so, skipping to step 2.9;It is no
Then, step 2.3 is skipped to;
Step 2.9:Judge whether the permission state set of traversal t period light-preserved systems, if so, skipping to step 2.10;If
It is no, skip to step 2.4;
Step 2.10:Judge whether t=T, finished if so, then calculating, output result;If it is not, then t=t+1, and skip to step
Rapid 2.3;
It is described by light-preserved system participate in distribution network operation it is optimal exert oneself with photovoltaic is actual exert oneself between imbalance power
Calculation formula be:
Pnet(t)=POpt(t)-PPV(t)
In formula:Pnet(t) it is imbalance power;POpt(t) participate in that distribution network operation is optimal exerts oneself for light-preserved system;PPV(t)
For photovoltaic generation unit active power of output.
Further, the stage t is:One full schedule cycle T is divided into several periods, the single time is remembered
Section is the stage, and stage sequence number is labeled as t, t ∈ { 1,2,3 ... T }, and the adjacent phases time difference is Δ t;
State S (the sPt,sQt) be:The dump energy S for choosing energy-storage system in light-preserved system is used as state SPt, and will
Electricity difference between its discretization, adjacent states is Δ S;Choose the remaining reactive capability of light-preserved system inverter and be used as state SQt, and
By its discretization, the reactive capability difference between adjacent states is Δ Q;
Decision-making U (the uPt,uQt) be:By the P in light-preserved system each periodPESSAnd Q (t)PESS(t) become as decision-making
Measure uPtAnd uQt, it must is fulfilled for light-preserved system operation constraint;
The strategy is:The sequence of the decision variable composition in each stage;
The state transition equation includes:SPtState transition equation and SQtState transition equation;The SPtState transfer
Equation is S (t+1)=(1- σ) S (t)+PESS(t)Δt;The SQtAbsolute transfer relationship is not present between adjacent two benches, it is permitted
Perhaps state set is by current state SPtDetermined with inverter capacity, its state transition equation is:
The target function is:It regard t phase targets function as stage target function Vt(S(sPt,sQt),U(uPt,
uQt)), then the target function in t stages is:
In formula:D(uPt,uQt) represent t stage conditions permission decision-making set, VtExpression stage target function.
Further, it is described by improve Hilbert-Huang transform decompose light-preserved system participate in distribution network operation it is optimal go out
Power and photovoltaic plant is actual exert oneself between imbalance power, by imbalance power be divided into high and low frequency two parts fluctuation point
Measure, detailed process is:
Using set empirical mode decomposition EEMD to signal Pnet(t) decomposed, after set empirical mode decomposition
Pnet(t) it is represented by:
In formula:Pnet(t) it is imbalance power signal;ci,0(t) it is i-th of intrinsic mode function IMF;R (t) is remaining point
Amount;For all intrinsic mode function averages;For all residual components averages;
To ci,0(t) recurrence Hilbert transform, signal c are carried outi,0(t) it is represented by:
In formula:ci,0(t) it is i-th of intrinsic mode function IMF;ai,j(t) it is intrinsic mode function ci,j(t) amplitude letter
Number;cosφi,m(t) it is recursive calculation to magnitude function ai,j(t) when being intended to 1, intrinsic mode function ci,m(t) instantaneous phase
Cosine value;M is decomposition number of times;
The instantaneous frequency for obtaining intrinsic mode function IMF is:
In formula:ωi(t) it is signal ci,0(t) instantaneous frequency;φi,m(t) it is recursive calculation to magnitude function ai,j(t) become
To in 1 when, intrinsic mode function ci,m(t) instantaneous phase, m is decomposition number of times;
Determine crossover frequency fgIf, crossover frequency fgBetween g and g+1 IMF, then c1(t),c2(t),…,cg(t) regard
For high-frequency fluctuation component, this portion of energy is stabilized by power-type energy storage;cg+1(t),cg+2(t),…,cn(t) it is considered as low-frequency fluctuation
Component, this portion of energy is stabilized by energy type energy storage, then the reference power of energy type energy storage and power-type energy storage is:
In formula:Pg,ref(t)、Pn,ref(t) be respectively power-type and energy type energy storage reference power;ci(t) it is natural mode
State function, its value substituted into is all intrinsic mode function averagesR (t) is residual components.
Further, it is described to be utilized respectively power-type energy storage and energy type energy storage compensation high-frequency fluctuation component and low-frequency fluctuation
The detailed process of component is:
Set up energy type energy storage rated power configuration calculation formula be:
In formula:Pn,rateFor the rated power of energy type energy storage;tn,0For energy type energy storage discharge and recharge initial time, TnTo grind
Study carefully duration;ηDC-DCAnd ηDC-ACFor the energy conversion efficiency of transverter;ηncAnd ηndFor the efficiency for charge-discharge of energy type energy storage;
Set up energy type energy storage rated capacity configuration calculation formula be:
In formula:En,rateFor energy type energy storage rated capacity;Pn(t) it is consideration transverter energy conversion efficiency and discharge and recharge
The charge-discharge electric power of energy type energy storage after efficiency;ΔTnFor energy type energy storage discharge and recharge time interval;SOCn,0For energy type energy storage
Initial state-of-charge;SOCn,max、SOCn,minFor the upper lower limit value of energy type energy storage charge state;
Work as SOCn,0Meet equation, and En,rateWhen taking equal sign, the rated capacity of energy type energy storage is minimum;
Set up power-type energy storage rated power configuration calculation formula be:
In formula:Pg,rateFor the rated power of power-type energy storage;tg,0For power-type energy storage discharge and recharge initial time;ηDC-DCWith
ηDC-ACFor the energy conversion efficiency of transverter;ηgcAnd ηgdFor the efficiency for charge-discharge of power-type energy storage;
Set up power-type energy storage rated capacity configuration calculation formula be:
In formula:Eg,rateFor power-type energy storage rated capacity;Pg(t) it is consideration transverter energy conversion efficiency and discharge and recharge
The charge-discharge electric power of power-type energy storage after efficiency;ΔTgFor power-type energy storage discharge and recharge time interval;SOCg,0For power-type energy storage
Initial state-of-charge;SOCg,max、SOCg,minFor the upper lower limit value of power-type energy storage charge state;
Work as SOCg,0Meet equation, and Eg,rateWhen taking equal sign, the rated capacity of power-type energy storage is minimum.
Further, the power-type energy storage includes:Ultracapacitor, flywheel energy storage, super conductive magnetic storage energy;The energy type
Energy storage includes lithium ion battery.
Further, described a kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network,
The mixed energy storage system capacity configuration scheme for proposing to optimize operation towards power distribution network, including:Set up hybrid energy-storing complete
Life cycle cost model, for evaluating the economy under different allocation plans, is held with obtaining most economical mixed energy storage system
Configuration scheme is measured, detailed process is:
Min LCC=Cinv+Crep+Com+Cscr-Cres
Cinv=CpinvPrate+CeinvErate
Cscr=(CpscrPrate+CescrErate)(n+1)(P/F,i,T)
Cres=σres(Cinv+Crep)(P/F,i,T)
In formula:Cinv、Crep、Com、Cscr、CresRespectively HESS initial outlay cost, renewal displacement cost, operation maintenance
Cost, scrap processing cost and reclaim residual value;Cpinv、Ceinv、Cprep、Ceinv、Cpom、Ceom、Cpscr、CescrRespectively unit power
Cost coefficient, unit capacity cost coefficient, unit power update cost coefficient, unit capacity and update cost coefficient, unit power
O&M cost coefficient, unit capacity O&M cost coefficient, unit power scrap processing cost coefficient and unit capacity scraps processing
Cost coefficient;Prate、ErateRespectively HESS rated power and HESS rated capacity;(P/F, i, t)=(1+i)-t;Wess
(t) it is year discharge and recharge;σresTo reclaim salvage value rate.
Further, the σresTake 3%~5%.
Beneficial effects of the present invention:
The present invention proposes a kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network, can be notable
The rated power and rated capacity of energy-storage system are reduced, is that HESS capacity is distributed rationally there is provided a kind of new method, with certain
Reference value;Meanwhile, HESS application can not only stabilize the fluctuation that photovoltaic is exerted oneself, and promote dissolving for photovoltaic energy, and
The optimization operation of power distribution network can be participated in, the economy and reliability of distribution network operation can be effectively improved.
Brief description of the drawings
Fig. 1 is the structure chart of light-preserved system in an embodiment of the present invention;
Fig. 2 is the photovoltaic plant schematic diagram containing HESS;
Fig. 3 is typical day operation curve map;
Fig. 4 is imbalance power EEMD decomposition result figures;
Fig. 5 is IMF instantaneous frequencys-time plot.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to embodiments, to the present invention
It is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not used to
Limit the present invention.
The application principle of the present invention is explained in detail below in conjunction with the accompanying drawings.
A kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network of the present invention, is obtained first
Light-preserved system participate in distribution network operation it is optimal exert oneself with photovoltaic is actual exert oneself between imbalance power, as mixing store up
The reference power of energy, is decomposed into a series of intrinsic mode functions, by passing using set empirical mode decomposition by imbalance power
Hilbert transform is returned to obtain instantaneous frequency-time graph of each intrinsic mode function;According to power-type and energy type energy storage
Characteristic, at least determines crossover frequency, using power-type energy storage and energy type energy storage with instantaneous frequency-time graph aliasing for principle
The high fdrequency component and low frequency component to imbalance power are stabilized respectively;Finally, it is considered to efficiency for charge-discharge, state-of-charge and complete
Life cycle cost, proposes power-type and energy type energy storage rated power and rated capacity economic optimum allocation plan.The present invention
Relative to conventional method, it can not only compensate new energy and exert oneself, improve the digestion capability to new energy, and can be actively engaged in
The optimization operation of power distribution network.
Embodiment one
Light-preserved system in the embodiment of the present invention is as shown in figure 1, C1For unidirectional DC/DC converters, main realization boosting work(
Can be with MPPT maximum power point tracking (MPPT);C2For two-way DC/DC converters, the controllable discharge and recharge of power is realized, makes system grid connection
The active optimal active power to participate in power distribution network optimization operation;C3For DC/AC inverters (i.e. photovoltaic DC-to-AC converter), for stable
DC bus-bar voltage, is realized active and reactive optimal grid-connected;Photovoltaic generation unit (photovoltaic battery array in such as Fig. 1) and energy storage
Unit is collected to common DC bus by respective DC/DC current transformers, then is connected to power distribution system through public DC/AC inverters
System;It is assumed that light-preserved system is using to power network power output as positive direction.
A kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network in the embodiment of the present invention,
Comprise the following steps:
1st, set up light-preserved system and participate in power distribution network optimal operation model, solve and obtain light-preserved system participation distribution network operation
It is optimal to exert oneself, and by light-preserved system participate in distribution network operation it is optimal exert oneself with photovoltaic is actual exert oneself between imbalance power make
For the reference power of mixed energy storage system;
Foregoing light-preserved system, which participates in power distribution network optimal operation model, includes object function and constraints;
The embodiment of the present invention chooses the whole network loss minimization as object function, is stored up for the light proposed in the embodiment of the present invention
System, the whole network network loss is except the loss P of transmission line of electricityL.loss(t) (be represented by each node of system active injection power it
With), in addition it is also necessary to consider the running wastage P of light-preserved systemPESS.loss(t) (the predominantly efficiency for charge-discharge of energy-storage system and its charge and discharge
Energy loss when electricity condition is changed), therefore:The whole network loss minimization object function is:
PPESS.loss(t)=(1- ηch)Pch(t)+(1-ηdis)Pdis(t)+ξ(t)·ΔPESS.loss
In formula:T is the cycle of operation;Δ t is operation step-length;N is system node number, Pi(t) it is the injection of t period node is
Active power;ξ (t) is that energy-storage system changes dimension from the t-1 periods to the charging and discharging state of t periods, and value is 0 or 1;ηch、ηdis
For the efficiency for charge-discharge of energy-storage system;Pch(t)、Pdis(t) it is the charge-discharge electric power of t period energy-storage systems;ΔPESS.lossDuring for t
The loss that section energy-storage system is produced when charging and discharging state switches, preferably takes the 0.5% of energy-storage system rated capacity.
The constraints bag in the constraints of power distribution network optimization operation, the embodiment of the present invention is participated in view of light-preserved system
Include:System load flow constraint, the constraint of feeder line outlet power, node voltage constraint, line power constraint and light-preserved system operation are about
Beam;
The grid trend constraint is:
The substation feeder outlet power is constrained to:
The node voltage is constrained to:
The line power is constrained to:
The light-preserved system operation is constrained to:
PPESS(t)=PPV(t)+PESS(t)
PPV(t)=PPV.MPPT(t)
-Pch.max≤PESS(t)≤Pdis.max
Because energy-storage system state of charge has absolute continuity in sequential, it is in strict accordance with time sequencing according to filling
Electric discharge
Watt level carries out cumulative calculation, therefore light-preserved system operation constraint also includes:
S (t+1)=(1- σ) S (t)+PESS(t)Δt
Smin≤S(t)≤Smax
S (0)=S (T)
In above-mentioned formula:I, j are system node number;Ui(t)、Uj(t) it is t periods node i and node j point voltage magnitudes;
Gij、BijTransconductance and mutual susceptance respectively between node i and node j;δij(t) it is the phase between t periods node i and node j
Potential difference;PKi(t)、QKi(t) it is respectively that t periods feeder line exports active and reactive power;PPESSi(t)、QPESSi(t) it is t period nodes
The active reactive of light-preserved system (i.e. inverter C3) is exerted oneself at i;PDi(t)、QDi(t) for load active at t period node is and
Reactive power;Pij(t) it is the line power between t node i and node j;Subscript "-" and subscript " _ " difference table of variable
Show the upper and lower bound of variable;PPESSAnd Q (t)PESS(t) be respectively t periods inverter output active and reactive power;PPV
(t) it is t period photovoltaic generation unit active power of output;PMPPTIt is DC/DC current transformers according to maximum during MPPT control strategies
Output is active;PESS(t) it is the active output of t period energy-storage systems, is divided into electric discharge, free time, three kinds of states of charging;Pch(t) and
Pdis(t) be respectively t period energy-storage systems charge-discharge electric power;ηchAnd ηdisRespectively energy-storage system efficiency for charge-discharge;Pch.maxWith
Pdis.maxRespectively charge-discharge electric power bound;S (t) is the dump energy of t period energy-storage systems;σ is put certainly for energy-storage system
Electric rate;SminAnd SmaxThe respectively bound of energy-storage system dump energy.
The solution obtains light-preserved system participation the optimal of distribution network operation and exerted oneself, and is specially:Using two-dimension dynamic programming
Algorithm for Solving light-preserved system participates in power distribution network optimal operation model, show that light-preserved system participates in the optimal of distribution network operation and exerted oneself;
More specifically comprise the following steps:
Step 2.1:Light-preserved system participation power distribution network optimal operation model is converted into the operable model of two-dimension dynamic programming,
Including stage t, state vector S (sPt,sQt), decision vector U (uPt,uQt), strategy, state transition equation and target function;
Step 2.2:Initialization, light-preserved system and distribution network system primary data needed for input;
Step 2.3:Determine the permission state set of t-1 period light-preserved systems;
Step 2.4:It is determined that meeting the permission state set of t period light-preserved systems under constraints;
Step 2.5:Calculate correspondence decision variable and target function value;
Step 2.6:Whether judge index function is optimal, if so, then skipping to step 2.7;Otherwise, step 2.8 is skipped to;
Step 2.7:Preserve current state variable, decision variable and object function;
Step 2.8:Judge whether the permission state set of traversal t-1 period light-preserved systems, if so, skipping to step 2.9;It is no
Then, step 2.3 is skipped to;
Step 2.9:Judge whether the permission state set of traversal t period light-preserved systems, if so, skipping to step 2.10;If
It is no, skip to step 2.4;
Step 2.10:Judge whether t=T, finished if so, then calculating, output result;If it is not, then t=t+1, and skip to step
Rapid 2.3.
Further, the stage t is:One full schedule cycle T is divided into several periods, remembers that the period is
Stage, stage sequence number is labeled as t, t ∈ { 1,2,3 ... T }, and the adjacent phases time difference is Δ t;
State S (the sPt,sQt) be:The dump energy S for choosing energy-storage system in light-preserved system is used as state SPt, and will
Electricity difference between its discretization, adjacent states is Δ S;Choose the remaining reactive capability of light-preserved system inverter and be used as state SQt, and
By its discretization, the reactive capability difference between adjacent states is Δ Q;
Decision-making U (the uPt,uQt) be:By the P in light-preserved system each periodPESSAnd Q (t)PESS(t) become as decision-making
Measure uPtAnd uQt, it must is fulfilled for light-preserved system operation constraint;
The strategy is:The sequence of the decision variable composition in each stage;
The state transition equation includes:SPtState transition equation and SQtState transition equation;The SPtState
Equation of transfer is S (t+1)=(1- σ) S (t)+PESS(t) Δ t, (i.e. energy-storage system electricity and charge-discharge electric power relation equation);Institute
State SQtIntersegmental in the absence of absolute transfer relationship in adjacent time, it allows state set by current state SPtHold with inverter
Amount determines that its state transition equation is
The target function is:It regard t phase targets function as stage target function Vt(S(sPt,sQt),U(uPt,
uQt)), then the in-process metrics function in t stages is:
In formula:D(uPt,uQt) represent t stage conditions permission decision-making set, VtExpression stage target function.
It is described by light-preserved system participate in distribution network operation it is optimal exert oneself with photovoltaic is actual exert oneself between imbalance power
Calculation formula be:
Pnet(t)=POpt(t)-PPV(t)
In formula:Pnet(t) it is imbalance power;POpt(t) participate in that distribution network operation is optimal exerts oneself for light-preserved system;PPV(t)
For photovoltaic generation unit active power of output.
2nd, decomposing light-preserved system by improving Hilbert-Huang transform participates in that distribution network operation is optimal exerts oneself and photovoltaic plant
Imbalance power between actually exerting oneself, high and low frequency two parts wave component, and profit respectively are divided into by imbalance power
With power-type energy storage and energy type energy storage compensation high-frequency fluctuation component and low-frequency fluctuation component;
In the embodiment of the present invention, participate in that distribution network operation is optimal exerts oneself and light using mixed energy storage system light-preserved system of dissolving
Overhead utility is actual exert oneself between imbalance power, its structure is as shown in Figure 2.Photovoltaic battery array by DC/DC converters with
Dc bus is connected, the mixed energy storage system being made up of energy type energy storage and power-type energy storage by two-way DC/DC converters with
Dc bus is connected, absorbed by hybrid energy-storing, delivered power to stabilize photovoltaic plant real output, meet light-preserved system ginseng
Optimal with distribution network operation is exerted oneself.
Imbalance power is obtained a series of intrinsic mode functions by the embodiment of the present invention by gathering empirical mode decomposition
(IMF) all IMF instantaneous frequency-time graph, and using recurrence Hilbert transform is obtained.According to two adjacent IMF it
Between frequency alias minimum principle determine the crossover frequency of imbalance power, imbalance power is resolved into high and low frequency two parts
Wave component.In figure:Pn,refFor energy type energy storage power output reference value, i.e., all low frequency component sums;Pg,refFor power-type
Energy storage power output reference value, i.e., all high fdrequency component sums;Pn,rate、Pg,rateRespectively energy type energy storage and power-type energy storage
Specified charge-discharge electric power;SOCn,max、SOCn,minThe respectively maximum and minimum value of energy type energy storage charge state;SOCn
(t) it is the state-of-charge of energy type energy storage t;SOCg,max、SOCg,minThe respectively maximum of power-type energy storage charge state
And minimum value;SOCg(t) it is the state-of-charge of power-type energy storage t.
Improving Hilbert-Huang transform includes set empirical mode decomposition (EEMD) and recurrence Hilbert transform two parts.
EEMD essence is into a series of data sequences with different characteristic time scale, i.e. intrinsic mode function by signal decomposition
(IMF), each IMF is imbalance power P in factnet(t) a kind of oscillation mode, is obtained by recurrence Hilbert transform
Each IMF instantaneous frequency-time graph, so that the degree of accuracy of its instantaneous frequency of IMF is improved, and EEMD decomposable processes determine ci
(t) instantaneous frequency is generally higher than ci+1(t) instantaneous frequency, therefore need to determine optimal in many instantaneous frequency-time graphs
Crossover frequency so that frequency alias is minimum between two adjacent IMF, and imbalance power is divided into high and low frequency two parts ripple
Dynamic component, and it is utilized respectively power-type energy storage and energy type energy storage compensation high-frequency fluctuation component and low-frequency fluctuation component.Below will
The detailed process for being improved Hilbert-Huang transform decomposition is illustrated.
1) EEMD decomposes imbalance power
Light-preserved system participate in distribution network operation it is optimal exert oneself and photovoltaic plant is actual exert oneself between imbalance power be:
Pnet(t)=POpt(t)-PPV(t)
In formula:POpt(t) participate in that distribution network operation is optimal exerts oneself for light-preserved system;PPV(t) it is the reality output of photovoltaic plant
Power.
Because imbalance power is nonlinear change, the embodiment of the present invention is decomposed using EEMD to it, EEMD points
Solution adds white noise acoustic disturbance signal to be decomposed, and carries out multiple empirical mode decomposition (EMD), obtained by repeatedly decomposing
IMF population means be considered as final IMF.This method can be prevented effectively from modal overlap, improve EMD discomposing effects, specific steps
It is as follows:
If primary signal is x (t), addition white Gaussian noise disturbance w (x) is obtained on x (t):
sj(t)=x (t)+wj(x), j=1,2 ..., m
In formula:sj(t) the imbalance power signal added for jth time after white Gaussian noise disturbance;wj(x) to be separate
White noise acoustic disturbance;M is that EMD carries out number of times.
sj(t) signal resolves into n IMF, i.e. c by EMD1j(t),c2j(t),…,cnj, and a remainder r (t)j
(t) it is
sj(t)=c1j(t)+c2j(t)+…+cnj(t)+rj(t)
In formula:sj(t) the imbalance power signal added for jth time after white Gaussian noise disturbance;cnj(t) it is to pass through EMD
N-th of the intrinsic mode function decomposed;rj(t) remainder after being decomposed for EMD.
Decomposed by m EMD, obtain m group decomposition results:
[c11(t),c21(t),…ci1(t),…cn1(t),r1(t)],
[c1j(t),c2j(t),…cij(t),…cnj(t),rj(t)],
[c1m(t),c2m(t),…cim(t),…cnm(t),rm(t)]。
In formula:cij(t) i-th obtained of IMF is decomposed for jth time;rj(t) obtained remainder is decomposed for jth time.
Obtained all IMF are decomposed by m EMD and residual components its averages is:
In formula:To decompose the average of obtained all IMF and residual components by m EMD;M is EMD points
Solve number of times.
Therefore, the P after EEMD is decomposednet(t) signal is represented by:
In formula:Pnet(t) it is imbalance power signal;For by m EMD decompose obtained all IMF with
The average of residual components.
2) IMF recurrence Hilbert transform
Primary signalIf ci,0(t) it is H [c after Hilbert transformi,0(t)], then:
Analytic signal zi,0(t) it is represented by:
zi,0(t)=ci,0(t)+jH[ci,0(t)]=ai,0(t)exp[jφi,0(t)]
In formula:ai,0(t) it is ci,0(t) instantaneous amplitude,φi,0(t) it is ci,0(t)
Instantaneous phase,
By ci,0(t) it is expressed as instantaneous amplitude ai,0(t) with pure FM signal cos φi,0(t) product, then:
ci,0(t)=ai,0(t)cosφi,0(t)
In formula:ci,0(t) it is instantaneous amplitude ai,0(t) with pure FM signal cos φi,0(t) product.
With pure FM Function ci,1(t)=cos φi,0(t) as new signal, and Hilbert transform is carried out to it, obtained
To new instantaneous amplitude and pure FM signal.Recursive calculation is constantly carried out, its recurrence formula is:
Often corresponding magnitude function and phase function are just can obtain by a recursive calculation:
Recursive calculation is constantly carried out until magnitude function ai,j(t) it is intended to 1, then ci,m+1(t)=cos φi,m(t) it is, comprehensive
Above-mentioned recursive computing steps, signal ci,0(t) it is represented by:
Instantaneous frequency is represented by:
In formula:ωi(t) it is instantaneous frequency;φi,m(t) it is instantaneous phase.
3) crossover frequency is determined
IMF obtains c by recurrence Hilbert transformi(t) instantaneous frequency-time graph, due to i increase, ci
(t) waveform is more smooth, therefore fi(t) instantaneous frequency is also lower.Adjacent f two-by-twoiAnd f (t)i+1(t) always there is frequency to mix
Folded phenomenon, therefore need to determine that optimal crossover frequency causes ciAnd c (t)i+1(t) energy aliasing is minimum.
If crossover frequency fgBetween g and g+1 IMF, then c1(t),c2(t),…,cg(t) it is considered as high-frequency fluctuation point
Amount, this portion of energy is stabilized by power-type energy storage;cg+1(t),cg+2(t),…,cn(t) it is considered as low-frequency fluctuation component, this part energy
Amount is stabilized by energy type energy storage, then the reference power of energy type energy storage and power-type energy storage is:
In formula:Pg,ref(t)、Pn,ref(t) it is power-type and the reference power of energy type energy storage;ci(t) it is natural mode of vibration letter
Number, its value substituted into is all intrinsic mode function averagesTo be every in the m group decomposition result matrixes that are previously obtained
The average of one row);R (t) is residual components.
The charge-discharge characteristic of energy type energy storage and power-type energy storage is taken into full account, the use longevity of energy storage device is considerably increased
Life, can effectively reduce system operation cost.
3rd, obtain the mixed energy storage system capacity for optimizing operation towards power distribution network and distribute strategy rationally.
It is described propose be towards the detailed process of mixed energy storage system capacity configuration scheme that power distribution network optimizes operation:
The rated power of energy type energy storage is set as Pn,rate, it is considered to transverter energy conversion efficiency and energy type energy storage
Efficiency for charge-discharge, the rated power of energy type energy storage configures calculation formula and is:
In formula:Pn,rateFor the rated power of energy type energy storage;tn,0For energy type energy storage discharge and recharge initial time;ηDC-DCWith
ηDC-ACFor the energy conversion efficiency of transverter;ηncAnd ηndFor the efficiency for charge-discharge of energy type energy storage;Pn,ref(t) stored up for energy type
The reference power of energy.
The rated capacity configuration of energy type energy storage is calculated:
If the initial state-of-charge of energy type energy storage is SOCn,0, the state-of-charge at k moment is:
In formula:ΔTnFor energy type energy storage discharge and recharge time interval;SOCn,kFor the state-of-charge of k moment energy type energy storage;
SOCn,min、SOCn,maxFor the bound of energy type energy storage charge state;Pn(t) it is consideration transverter energy conversion efficiency and charge and discharge
The charge-discharge electric power of energy type energy storage after electrical efficiency, its calculation formula is:
In formula:Pn,ref(t) it is the reference power of energy type energy storage;ηDC-DCAnd ηDC-ACFor the energy conversion efficiency of transverter;
ηncAnd ηndFor the efficiency for charge-discharge of energy type energy storage.
A length of T when researchnWhen, energy type energy storage rated capacity En,rateCalculation formula be
In formula:En,rateFor energy type energy storage rated capacity;Pn(t) it is consideration transverter energy conversion efficiency and discharge and recharge
The charge-discharge electric power of energy type energy storage after efficiency;ΔTnFor energy type energy storage discharge and recharge time interval;SOCn,0For energy type energy storage
Initial state-of-charge;SOCn,min、SOCn,maxFor the bound of energy type energy storage charge state
Work as SOC0Meet equation, and En,rateWhen taking equal sign, the rated capacity of energy type energy storage is minimum.
In formula:Pn(t) it is the charge and discharge electric work of energy type energy storage after consideration transverter energy conversion efficiency and efficiency for charge-discharge
Rate;ΔTnFor energy type energy storage discharge and recharge time interval;SOCn,0For the initial state-of-charge of energy type energy storage;SOCn,max、
SOCn,minFor the upper lower limit value of energy type energy storage charge state.
The rated power configuration of power-type energy storage is calculated:
The rated power of power-type energy storage is set as Pg,rate, it is considered to transverter energy conversion efficiency and power-type energy storage
Efficiency for charge-discharge, the rated power of power-type energy storage configures calculation formula and is:
In formula:Pg,rateFor the rated power of power-type energy storage;tg,0For power-type energy storage discharge and recharge initial time;ηDC-DCWith
ηDC-ACFor the energy conversion efficiency of transverter;ηgcAnd ηgdFor the efficiency for charge-discharge of power-type energy storage;Pg,ref(t) stored up for power-type
The reference power of energy.
The rated capacity configuration of power-type energy storage is calculated:
If the initial state-of-charge of power-type energy storage is SOCg,0, the state-of-charge at k moment is:
In formula:ΔTgFor power-type energy storage discharge and recharge time interval;SOCg,kFor the state-of-charge of k moment power-type energy storage;
SOCg,min、SOCg,maxFor the bound of power-type energy storage charge state;Pg(t) it is consideration transverter energy conversion efficiency and charge and discharge
The charge-discharge electric power of power-type energy storage after electrical efficiency, its calculation formula is:
In formula:Pg,ref(t) it is the reference power of power-type energy storage;ηDC-DCAnd ηDC-ACFor the energy conversion efficiency of transverter;
ηgcAnd ηgdFor the efficiency for charge-discharge of power-type energy storage.
A length of T when researchgWhen, power-type energy storage rated capacity Eg,rateCalculation formula be
In formula:Eg,rateFor power-type energy storage rated capacity;Pg(t) it is consideration transverter energy conversion efficiency and discharge and recharge
The charge-discharge electric power of power-type energy storage after efficiency;ΔTgFor power-type energy storage discharge and recharge time interval;SOCg,0For power-type energy storage
Initial state-of-charge;SOCg,max、SOCg,minFor the upper lower limit value of power-type energy storage charge state;
Work as SOCg,0Meet equation, and Eg,rateWhen taking equal sign, the rated capacity of power-type energy storage is minimum.
In formula:Eg,rateFor power-type energy storage rated capacity;Pg(t) it is consideration transverter energy conversion efficiency and discharge and recharge
The charge-discharge electric power of power-type energy storage after efficiency;ΔTgFor power-type energy storage discharge and recharge time interval;SOCg,0For power-type energy storage
Initial state-of-charge;SOCg,max、SOCg,minFor the upper lower limit value of power-type energy storage charge state.
4th, hybrid energy-storing economic evaluation
The embodiment of the present invention considers the input cost of hybrid energy-storing equipment, sets up the complete of mixed energy storage system (HESS)
Life cycle cost model,.Overall life cycle cost LCC refers to all produced in whole life cycle directly, indirectly
Etc. all expenses, it is assumed that the HESS life-spans are T, discount rate is i, and HESS elements update displacement number of times for n, then life cycle management into
This model is:
Min LCC=Cinv+Crep+Com+Cscr-Cres
Cinv=CpinvPrate+CeinvErate
Cscr=(CpscrPrate+CescrErate)(n+1)(P/F,i,T)
Cres=σres(Cinv+Crep)(P/F,i,T)
In formula:Cinv、Crep、Com、Cscr、CresRespectively HESS initial outlay cost, renewal displacement cost, operation maintenance
Cost, scrap processing cost and reclaim residual value;Cpinv、Ceinv、Cprep、Ceinv、Cpom、Ceom、Cpscr、CescrRespectively unit power
Cost coefficient, unit capacity cost coefficient, unit power update cost coefficient, unit capacity and update cost coefficient, unit power
O&M cost coefficient, unit capacity O&M cost coefficient, unit power scrap processing cost coefficient and unit capacity scraps processing
Cost coefficient;Prate、ErateRespectively HESS rated power and HESS rated capacity (calculate energy type energy storage when, PrateJust
Represent Pn,rate;When calculating power-type energy storage, PrateJust represent Pg,rate, similarly push away to obtain Erate);(P/F, i, t)=(1+i)-t;
Wess(t) it is year discharge and recharge;σresTo reclaim salvage value rate, 3%~5% is taken.
Embodiment two
(adopted using the real output data of Jiangsu Province's photovoltaic plant certain typical day in 2013 in the embodiment of the present invention
The sample time is 1min), simulation analysis are carried out in Matlab.One day excellent is carried out to light-preserved system using 33 Node power distribution systems
Change, be utilized respectively photovoltaic power Forecasting Methodology and load forecasting method obtains photovoltaic plant sunrise force curve and load day operation is bent
Line, as shown in figure 3, the light-preserved system participation active distribution network optimal operation model that the present invention is carried is built, using two-dimentional dynamic
Planning algorithm solving model, obtains light-preserved system participation the optimal of distribution network operation and exerts oneself.The typical daylight storage system is participated in
Distribution network operation is optimal to exert oneself with the actual difference exerted oneself of photovoltaic as imbalance power data, and Fig. 4 gives imbalance power
EEMD analysis results.As seen from the figure, intrinsic mode function ci(t) amplitude and frequency reduces with i increase, and zero
Value is oscillated around.Instantaneous frequency-the time graph obtained with reference to recurrence Hilbert transform, as shown in Figure 5.
Knowable to analysis, c4And c (t)5(t) aliasing is minimum, therefore takes i=4 as frequency partition separation, by uneven work(
Rate is divided into high-low frequency weight, then power-type energy storage and energy type energy storage power output reference value are Pg,ref(t)=c1(t)+c2
(t)+c3(t)+c4(t);Pn,ref(t)=c5(t)+c6(t)+r(t)。
With reference to the characteristic of power-type energy storage and energy type energy storage, the embodiment of the present invention is using the different configuration side of following four
Case.
Scheme 1 is lithium ion battery;Scheme 2 is lithium ion battery and ultracapacitor;Scheme 3 is lithium ion battery and flown
Take turns energy storage;Scheme 4 is lithium ion battery and super conductive magnetic storage energy.Each energy storage parameter as shown in table 1, is carried out different according to the above method
The configuration of rated power and rated capacity under scheme, and estimate mixed energy storage system (HESS) life-span, it is assumed that HESS life-cycle
Cycle is 10 years, then can calculate the cost of investment of HESS under each allocation plan, and this calculating process is prior art, is not gone to live in the household of one's in-laws on getting married
State, specific result of calculation is as shown in table 2.
Table 1
Table 2
By analyzing table 2, the rated power and rated capacity of traditional Hilbert-Huang transform configuration are substantially high
In the configuration result of the inventive method, such as scheme 2, the rated power and specified appearance of traditional Hilbert-Huang transform configuration
Amount is respectively 2.18MW/0.37MWh, 3.60MW/0.19MWh, the rated power and rated capacity of the inventive method configuration
Respectively 1.39MW/0.29MWh, 2.64MW/0.17MWh.It is far excellent using HESS schemes economy from the point of view of cost of investment
In single energy storage scheme, because HESS schemes can greatly increase the service life of lithium ion battery, so as to reduce the investment of energy storage
Cost, the wherein economy of scheme 2 are optimal, therefore can optimize the lower hybrid energy-storing capacity configuration of operation application as towards power distribution network
Suggested design.
The general principle and principal character and advantages of the present invention of the present invention has been shown and described above.The technology of the industry
Personnel are it should be appreciated that the present invention is not limited to the above embodiments, and the simply explanation described in above-described embodiment and specification is originally
The principle of invention, without departing from the spirit and scope of the present invention, various changes and modifications of the present invention are possible, these changes
Change and improvement all fall within the protetion scope of the claimed invention.The claimed scope of the invention by appended claims and its
Equivalent thereof.
Claims (10)
1. it is a kind of towards power distribution network optimize operation mixed energy storage system capacity configuration optimizing method, it is characterised in that including with
Lower step:
Set up light-preserved system participate in power distribution network optimal operation model, solve obtain light-preserved system participate in distribution network operation it is optimal go out
Power, and using light-preserved system participate in distribution network operation it is optimal exert oneself with photovoltaic is actual exert oneself between imbalance power as mixing
The reference power of energy-storage system;
By improve Hilbert-Huang transform decompose light-preserved system participate in distribution network operation it is optimal exert oneself with photovoltaic plant it is actual go out
Imbalance power between power, is divided into high and low frequency two parts wave component, and be utilized respectively power by imbalance power
Type energy storage and energy type energy storage compensation high-frequency fluctuation component and low-frequency fluctuation component;
Obtain the mixed energy storage system capacity for optimizing operation towards power distribution network and distribute strategy rationally.
2. a kind of mixed energy storage system capacity side of distributing rationally for optimizing operation towards power distribution network according to claim 1
Method, it is characterised in that:The light-preserved system, which participates in power distribution network optimal operation model, includes object function and constraints,
The object function is:
PPESS.loss(t)=(1- ηch)Pch(t)+(1-ηdis)Pdis(t)+ξ(t)·ΔPESS.loss
In formula:T is the cycle of operation;Δ t is operation step-length;N is distribution system nodes, Pi(t) being injected with for t period node is
Work(power;ξ (t) is that energy-storage system changes dimension from the t-1 periods to the charging and discharging state of t periods, and value is 0 or 1;ηch、ηdisFor
The efficiency for charge-discharge of energy-storage system;Pch(t)、Pdis(t) it is the charge-discharge electric power of t period energy-storage systems;ΔPESS.lossFor the t periods
The loss that energy-storage system is produced when charging and discharging state switches.
3. a kind of mixed energy storage system capacity side of distributing rationally for optimizing operation towards power distribution network according to claim 2
Method, it is characterised in that:The constraints is:
PPESS(t)=PPV(t)+PESS(t)
PPV(t)=PPV.MPPT(t)
-Pch.max≤PESS(t)≤Pdis.max
S (t+1)=(1- σ) S (t)+PESS(t)Δt
Smin≤S(t)≤Smax
S (0)=S (T)
In formula:I, j are distribution network systems node number;Ui(t)、Uj(t) it is t periods node i and node j point voltage magnitudes;Gij、BijPoint
Transconductance and mutual susceptance that Wei be between node i and node j;δij(t) it is the phase difference between t periods node i and node j;PKi
(t)、QKi(t) it is respectively that the outlet of t periods feeder line is active and idle;PPESSi(t)、QPESSi(t) it is light-preserved system at t period node is
Active and reactive power;PDi(t)、QDi(t) it is the active and reactive power of load at t period node is;Pij(t) it is t
Line power between node i and node j;The subscript "-" and subscript " _ " of variable represent the upper and lower bound of variable respectively;
PPESSAnd Q (t)PESS(t) be respectively t periods light-preserved system output active and reactive power;SPESS.maxFor the specified appearance of inverter
Amount;PPV(t) it is t period photovoltaic generation unit active power of output;PMPPTIt is photovoltaic DC/DC current transformers according to MPPT control strategies
When maximum output it is active;PESS(t) it is the active output of t period energy-storage systems, is divided into electric discharge, free time, three kinds of states of charging;
PchAnd P (t)dis(t) be respectively t period energy-storage systems charge-discharge electric power;ηchAnd ηdisRespectively energy-storage system discharge and recharge is imitated
Rate;Pch.maxAnd Pdis.maxRespectively charge-discharge electric power bound;S (t) is the dump energy of t period energy-storage systems;σ is energy storage
The self-discharge rate of system;SminAnd SmaxThe respectively bound of energy-storage system dump energy.
4. a kind of mixed energy storage system capacity side of distributing rationally for optimizing operation towards power distribution network according to claim 1
Method, it is characterised in that:The solution obtains light-preserved system participation the optimal of distribution network operation and exerted oneself, and is specially:
Power distribution network optimal operation model is participated in using two-dimension dynamic programming Algorithm for Solving light-preserved system, show that light-preserved system participates in matching somebody with somebody
The optimal of operation of power networks is exerted oneself;Specifically include following steps:
Step 2.1:Light-preserved system participation power distribution network optimal operation model is converted into the operable model of two-dimension dynamic programming, including
Stage t, state vector S (sPt,sQt), decision vector U (uPt,uQt), strategy, state transition equation and target function;
Step 2.2:Initialization, light-preserved system and distribution network system primary data needed for input;
Step 2.3:Determine the permission state set of t-1 period light-preserved systems;
Step 2.4:It is determined that meeting the permission state set of t period light-preserved systems under constraints;
Step 2.5:Calculate correspondence decision variable and target function value;
Step 2.6:Whether judge index function is optimal, if so, then skipping to step 2.7;Otherwise, step 2.8 is skipped to;
Step 2.7:Preserve current state variable, decision variable and object function;
Step 2.8:Judge whether the permission state set of traversal t-1 period light-preserved systems, if so, skipping to step 2.9;Otherwise,
Skip to step 2.3;
Step 2.9:Judge whether the permission state set of traversal t period light-preserved systems, if so, skipping to step 2.10;If it is not, jumping
To step 2.4;
Step 2.10:Judge whether t=T, finished if so, then calculating, output result;If it is not, then t=t+1, and skip to step
2.3;
It is described by light-preserved system participate in distribution network operation it is optimal exert oneself with photovoltaic is actual exert oneself between imbalance power meter
Calculating formula is:
Pnet(t)=POpt(t)-PPV(t)
In formula:Pnet(t) it is imbalance power;POpt(t) participate in that distribution network operation is optimal exerts oneself for light-preserved system;PPV(t) it is light
Lie prostrate generator unit active power of output.
5. a kind of mixed energy storage system capacity side of distributing rationally for optimizing operation towards power distribution network according to claim 4
Method, it is characterised in that:
The stage t is:One full schedule cycle T is divided into several periods, remembers that the single period is stage, rank
Section sequence number is labeled as t, t ∈ { 1,2,3 ... T }, and the adjacent phases time difference is Δ t;
State S (the sPt,sQt) be:The dump energy S for choosing energy-storage system in light-preserved system is used as state SPt, and by its from
Electricity difference between dispersion, adjacent states is Δ S;Choose the remaining reactive capability of light-preserved system inverter and be used as state SQt, and by its
Reactive capability difference between discretization, adjacent states is Δ Q;
Decision-making U (the uPt,uQt) be:By the P in light-preserved system each periodPESSAnd Q (t)PESS(t) as decision variable uPt
And uQt, it must is fulfilled for light-preserved system operation constraint;
The strategy is:The sequence of the decision variable composition in each stage;
The state transition equation includes:SPtState transition equation and SQtState transition equation;The SPtState transfer
Equation is S (t+1)=(1- σ) S (t)+PESS(t)Δt;The SQtAbsolute transfer relationship is not present between adjacent two benches, its
Allow state set by current state SPtDetermined with inverter capacity, its state transition equation is:
The target function is:It regard t phase targets function as stage target function Vt(S(sPt,sQt),U(uPt,uQt)), then
The target function in t stages is:
In formula:Represent the permission decision-making set of t stage conditions, VtExpression stage target function.
6. a kind of mixed energy storage system capacity side of distributing rationally for optimizing operation towards power distribution network according to claim 1
Method, it is characterised in that:It is described by improve Hilbert-Huang transform decompose light-preserved system participate in distribution network operation it is optimal exert oneself and
Photovoltaic plant is actual exert oneself between imbalance power, imbalance power is divided into high and low frequency two parts wave component, its
Detailed process is:
Using set empirical mode decomposition EEMD to signal Pnet(t) decomposed, the P after set empirical mode decompositionnet
(t) it is represented by:
In formula:Pnet(t) it is imbalance power signal;ci,0(t) it is i-th of intrinsic mode function IMF;R (t) is residual components;
For all intrinsic mode function averages;For all residual components averages;
To ci,0(t) recurrence Hilbert transform, signal c are carried outi,0(t) it is represented by:
In formula:ci,0(t) it is i-th of intrinsic mode function IMF;ai,j(t) it is intrinsic mode function ci,j(t) magnitude function;
cosφi,m(t) it is recursive calculation to magnitude function ai,j(t) when being intended to 1, intrinsic mode function ci,m(t) instantaneous phase
Position cosine value;M is decomposition number of times;
The instantaneous frequency for obtaining intrinsic mode function IMF is:
In formula:ωi(t) it is signal ci,0(t) instantaneous frequency;φi,m(t) it is recursive calculation to magnitude function ai,j(t) it is intended to 1
When,
Intrinsic mode function ci,m(t) instantaneous phase, m is decomposition number of times;
Determine crossover frequency fgIf, crossover frequency fgBetween g and g+1 IMF, then c1(t),c2(t),…,cg(t) it is considered as height
Frequency wave component, this portion of energy is stabilized by power-type energy storage;cg+1(t),cg+2(t),…,cn(t) it is considered as low-frequency fluctuation component,
This portion of energy is stabilized by energy type energy storage, then the reference power of energy type energy storage and power-type energy storage is:
In formula:Pg,ref(t)、Pn,ref(t) be respectively power-type and energy type energy storage reference power;ci(t) it is natural mode of vibration letter
Number, its value substituted into is all intrinsic mode function averagesR (t) is residual components.
7. a kind of mixed energy storage system capacity side of distributing rationally for optimizing operation towards power distribution network according to claim 6
Method, it is characterised in that:It is described to be utilized respectively power-type energy storage and energy type energy storage compensation high-frequency fluctuation component and low-frequency fluctuation point
The detailed process of amount is:
Set up energy type energy storage rated power configuration calculation formula be:
In formula:Pn,rateFor the rated power of energy type energy storage;tn,0For energy type energy storage discharge and recharge initial time, TnDuring for research
It is long;ηDC-DCAnd ηDC-ACFor the energy conversion efficiency of transverter;ηncAnd ηndFor the efficiency for charge-discharge of energy type energy storage;
Set up energy type energy storage rated capacity configuration calculation formula be:
In formula:En,rateFor energy type energy storage rated capacity;Pn(t) after for consideration transverter energy conversion efficiency and efficiency for charge-discharge
The charge-discharge electric power of energy type energy storage;ΔTnFor energy type energy storage discharge and recharge time interval;SOCn,0For the initial of energy type energy storage
State-of-charge;SOCn,max、SOCn,minFor the upper lower limit value of energy type energy storage charge state;
Work as SOCn,0Meet equation, and En,rateWhen taking equal sign, the rated capacity of energy type energy storage is minimum;
Set up power-type energy storage rated power configuration calculation formula be:
In formula:Pg,rateFor the rated power of power-type energy storage;tg,0For power-type energy storage discharge and recharge initial time;TgDuring for research
It is long;ηDC-DCAnd ηDC-ACFor the energy conversion efficiency of transverter;ηgcAnd ηgdFor the efficiency for charge-discharge of power-type energy storage;
Set up power-type energy storage rated capacity configuration calculation formula be:
In formula:Eg,rateFor power-type energy storage rated capacity;Pg(t) after for consideration transverter energy conversion efficiency and efficiency for charge-discharge
The charge-discharge electric power of power-type energy storage;ΔTgFor power-type energy storage discharge and recharge time interval;SOCg,0For the initial of power-type energy storage
State-of-charge;SOCg,max、SOCg,minFor the upper lower limit value of power-type energy storage charge state;
Work as SOCg,0Meet equation, and Eg,rateWhen taking equal sign, the rated capacity of power-type energy storage is minimum.
8. a kind of mixed energy storage system capacity side of distributing rationally for optimizing operation towards power distribution network according to claim 7
Method, it is characterised in that:The power-type energy storage includes:Ultracapacitor, flywheel energy storage, super conductive magnetic storage energy;The energy type storage
Lithium ion battery can be included.
9. a kind of mixed energy storage system capacity side of distributing rationally for optimizing operation towards power distribution network according to claim 7
Method, it is characterised in that the proposition optimizes the mixed energy storage system capacity configuration scheme of operation, specific bag towards power distribution network
Include:Hybrid energy-storing overall life cycle cost model is set up, it is most economical to obtain for evaluating the economy under different allocation plans
Mixed energy storage system capacity configuration scheme, detailed process is:
Min LCC=Cinv+Crep+Com+Cscr-Cres
Cinv=CpinvPrate+CeinvErate
Cscr=(CpscrPrate+CescrErate)(n+1)(P/F,i,T)
Cres=σres(Cinv+Crep)(P/F,i,T)
In formula:Cinv、Crep、Com、Cscr、CresRespectively HESS initial outlay cost, update displacement cost, operation maintenance into
Originally, scrap processing cost and reclaim residual value;Cpinv、Ceinv、Cprep、Ceinv、Cpom、Ceom、Cpscr、CescrRespectively unit power into
This coefficient, unit capacity cost coefficient, unit power update cost coefficient, unit capacity and update cost coefficient, unit power fortune
Tie up cost coefficient, unit capacity O&M cost coefficient, unit power scraps processing cost coefficient and unit capacity is scrapped and is processed into
This coefficient;Prate、ErateRespectively HESS rated power and HESS rated capacity;(P/F, i, t)=(1+i)-t;Wess(t)
For year discharge and recharge;σresTo reclaim salvage value rate.
10. a kind of mixed energy storage system capacity side of distributing rationally for optimizing operation towards power distribution network according to claim 9
Method, it is characterised in that:The σresTake 3%~5%.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710337499.9A CN106998072A (en) | 2017-05-15 | 2017-05-15 | A kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710337499.9A CN106998072A (en) | 2017-05-15 | 2017-05-15 | A kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106998072A true CN106998072A (en) | 2017-08-01 |
Family
ID=59434594
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710337499.9A Pending CN106998072A (en) | 2017-05-15 | 2017-05-15 | A kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106998072A (en) |
Cited By (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107492901A (en) * | 2017-08-29 | 2017-12-19 | 广东电网有限责任公司电力科学研究院 | A kind of distributed energy storage system real-time optimization method and device |
CN107632239A (en) * | 2017-08-25 | 2018-01-26 | 南京理工大学 | A kind of photovoltaic based on IMF Energy-Entropies sends out line fault phase-selecting method |
CN107947231A (en) * | 2017-12-01 | 2018-04-20 | 国网江苏省电力有限公司电力科学研究院 | A kind of mixed energy storage system control method towards power distribution network optimization operation |
CN107994593A (en) * | 2017-12-08 | 2018-05-04 | 囯网河北省电力有限公司电力科学研究院 | Composite energy storage power system capacity collocation method and terminal device |
CN108110780A (en) * | 2018-01-29 | 2018-06-01 | 广东电网有限责任公司电力科学研究院 | A kind of isolated micro-capacitance sensor stored energy capacitance Optimal Configuration Method and device |
CN108521133A (en) * | 2018-04-18 | 2018-09-11 | 国网山西省电力公司电力科学研究院 | Fire-storage combined dispatching frequency modulation power distribution method based on set empirical mode decomposition |
CN109031946A (en) * | 2018-05-07 | 2018-12-18 | 中车工业研究院有限公司 | A kind of control method and device of mixed energy storage system |
CN109088421A (en) * | 2018-08-28 | 2018-12-25 | 河海大学 | Mixed energy storage system capacity configuration optimizing method based on FDM |
CN109167391A (en) * | 2018-10-11 | 2019-01-08 | 珠海吉瓦科技有限公司 | A kind of echelon battery energy storage power station energy management method and system based on set empirical mode decomposition |
CN109347120A (en) * | 2018-09-25 | 2019-02-15 | 中国电力科学研究院有限公司 | A kind of decision-making technique and system towards mains frequency response |
CN109412192A (en) * | 2018-12-07 | 2019-03-01 | 国网江苏省电力有限公司常州供电分公司 | From the back-to-back soft straight device operation method of energy storage multiterminal |
CN109698557A (en) * | 2017-10-23 | 2019-04-30 | 中国南方电网有限责任公司 | Battery energy storage system configuration method and device, computer equipment and storage medium |
CN110086181A (en) * | 2019-05-31 | 2019-08-02 | 浙江大学 | Power adaptive cooperative control method in electric-thermal alliance micro-capacitance sensor when a kind of off-grid operation |
CN110311396A (en) * | 2019-07-30 | 2019-10-08 | 太原理工大学 | A kind of alternating current-direct current mixing micro-capacitance sensor hybrid energy-storing capacity configuration optimizing method |
CN110535187A (en) * | 2019-08-21 | 2019-12-03 | 国网山东省电力公司济宁供电公司 | A kind of the energy dispatching method and system of the composite energy storage capacity of active distribution network |
CN110752614A (en) * | 2019-11-07 | 2020-02-04 | 山东大学 | Energy storage system control method and system |
CN110797896A (en) * | 2019-11-28 | 2020-02-14 | 国网江苏电力设计咨询有限公司 | Hybrid energy storage capacity optimal configuration method in micro-grid wind and solar energy storage system |
CN111146793A (en) * | 2019-12-18 | 2020-05-12 | 济南大学 | Photovoltaic-energy storage system capacity optimization design method and system based on power feature extraction |
CN111404199A (en) * | 2020-04-22 | 2020-07-10 | 江苏东源电器集团股份有限公司 | Micro-grid hybrid energy storage system based on flywheel energy storage and control method |
CN111654028A (en) * | 2020-06-15 | 2020-09-11 | 湘潭大学 | Final benefit hybrid energy storage optimized capacity allocation ratio based on optimization algorithm |
CN111754361A (en) * | 2020-06-29 | 2020-10-09 | 国网山西省电力公司电力科学研究院 | Energy storage capacity optimal configuration method and computing device of wind-storage combined frequency modulation system |
CN111900712A (en) * | 2020-08-31 | 2020-11-06 | 北方民族大学 | Direct-current micro-grid energy balance control method and system based on hybrid energy storage |
CN112182975A (en) * | 2020-10-10 | 2021-01-05 | 国网江苏省电力有限公司扬州供电分公司 | Power frequency domain decomposition optimization method based on empirical mode decomposition |
CN112260296A (en) * | 2020-10-10 | 2021-01-22 | 国网江苏省电力有限公司扬州供电分公司 | Hybrid energy storage system optimization method and device for stabilizing intermittent load |
CN112636367A (en) * | 2020-12-09 | 2021-04-09 | 南京工程学院 | Hybrid energy storage capacity optimal configuration method for stabilizing wind power fluctuation |
CN113629736A (en) * | 2021-08-12 | 2021-11-09 | 国网江苏省电力有限公司常州供电分公司 | Intraday rolling optimization method based on power distribution network hydrogen energy storage system |
CN113690941A (en) * | 2021-03-24 | 2021-11-23 | 清华大学 | Method, device, equipment and storage medium for configuring capacity of optical storage and charging microgrid |
CN113762626A (en) * | 2021-09-10 | 2021-12-07 | 国网江苏省电力有限公司扬州供电分公司 | Hybrid energy storage system capacity optimization method considering economy |
CN113822707A (en) * | 2021-09-10 | 2021-12-21 | 国网冀北电力有限公司电力科学研究院 | Output decision method and device for power market, computer equipment and storage medium |
CN113904352A (en) * | 2021-10-11 | 2022-01-07 | 国网福建省电力有限公司 | Power distribution optimization method and terminal for hybrid energy storage system |
CN114142473A (en) * | 2021-12-06 | 2022-03-04 | 浙江华云电力工程设计咨询有限公司 | Micro-grid configuration method and system under condition of unknown wind and solar capacity |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103580291A (en) * | 2013-11-27 | 2014-02-12 | 上海空间电源研究所 | Battery energy storage system integrating energy type and power type and control method thereof |
CN105226688A (en) * | 2015-10-12 | 2016-01-06 | 中国电力科学研究院 | Based on the polymorphic type energy storage system capacity configuration optimizing method of Chance-constrained Model |
CN105260797A (en) * | 2015-10-22 | 2016-01-20 | 华北电力大学 | Microgrid energy storage power station program control method based on economical evaluation |
CN106099964A (en) * | 2016-06-16 | 2016-11-09 | 南京工程学院 | A kind of energy-storage system participates in active distribution network runing adjustment computational methods |
CN106230012A (en) * | 2016-09-19 | 2016-12-14 | 华北电力大学(保定) | Ultracapacitor and the Optimal Configuration Method of accumulator capacity in grid-connected photovoltaic system |
CN106451504A (en) * | 2016-10-19 | 2017-02-22 | 上海电力设计院有限公司 | Control method and device for configuration cost of hybrid energy storage system |
-
2017
- 2017-05-15 CN CN201710337499.9A patent/CN106998072A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103580291A (en) * | 2013-11-27 | 2014-02-12 | 上海空间电源研究所 | Battery energy storage system integrating energy type and power type and control method thereof |
CN105226688A (en) * | 2015-10-12 | 2016-01-06 | 中国电力科学研究院 | Based on the polymorphic type energy storage system capacity configuration optimizing method of Chance-constrained Model |
CN105260797A (en) * | 2015-10-22 | 2016-01-20 | 华北电力大学 | Microgrid energy storage power station program control method based on economical evaluation |
CN106099964A (en) * | 2016-06-16 | 2016-11-09 | 南京工程学院 | A kind of energy-storage system participates in active distribution network runing adjustment computational methods |
CN106230012A (en) * | 2016-09-19 | 2016-12-14 | 华北电力大学(保定) | Ultracapacitor and the Optimal Configuration Method of accumulator capacity in grid-connected photovoltaic system |
CN106451504A (en) * | 2016-10-19 | 2017-02-22 | 上海电力设计院有限公司 | Control method and device for configuration cost of hybrid energy storage system |
Non-Patent Citations (3)
Title |
---|
张晴等: "净效益最大的平抑风电功率波动的混合储能容量配置方法", 《电工技术学报》 * |
王仕奎著: "《随机信号分析理论与实践》", 31 August 2016 * |
胡志祥等: "基于递归希尔伯特变换的振动信号解调和瞬时频率计算方法", 《振动与冲击》 * |
Cited By (44)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107632239A (en) * | 2017-08-25 | 2018-01-26 | 南京理工大学 | A kind of photovoltaic based on IMF Energy-Entropies sends out line fault phase-selecting method |
CN107632239B (en) * | 2017-08-25 | 2020-06-05 | 南京理工大学 | Photovoltaic outgoing line fault phase selection method based on IMF energy entropy |
CN107492901B (en) * | 2017-08-29 | 2020-04-07 | 广东电网有限责任公司电力科学研究院 | Distributed energy storage system real-time optimization method and device |
CN107492901A (en) * | 2017-08-29 | 2017-12-19 | 广东电网有限责任公司电力科学研究院 | A kind of distributed energy storage system real-time optimization method and device |
CN109698557A (en) * | 2017-10-23 | 2019-04-30 | 中国南方电网有限责任公司 | Battery energy storage system configuration method and device, computer equipment and storage medium |
CN107947231A (en) * | 2017-12-01 | 2018-04-20 | 国网江苏省电力有限公司电力科学研究院 | A kind of mixed energy storage system control method towards power distribution network optimization operation |
CN107947231B (en) * | 2017-12-01 | 2020-04-21 | 国网江苏省电力有限公司电力科学研究院 | Hybrid energy storage system control method for optimized operation of power distribution network |
CN107994593A (en) * | 2017-12-08 | 2018-05-04 | 囯网河北省电力有限公司电力科学研究院 | Composite energy storage power system capacity collocation method and terminal device |
CN108110780A (en) * | 2018-01-29 | 2018-06-01 | 广东电网有限责任公司电力科学研究院 | A kind of isolated micro-capacitance sensor stored energy capacitance Optimal Configuration Method and device |
CN108521133A (en) * | 2018-04-18 | 2018-09-11 | 国网山西省电力公司电力科学研究院 | Fire-storage combined dispatching frequency modulation power distribution method based on set empirical mode decomposition |
CN109031946A (en) * | 2018-05-07 | 2018-12-18 | 中车工业研究院有限公司 | A kind of control method and device of mixed energy storage system |
CN109088421A (en) * | 2018-08-28 | 2018-12-25 | 河海大学 | Mixed energy storage system capacity configuration optimizing method based on FDM |
CN109347120B (en) * | 2018-09-25 | 2022-11-25 | 中国电力科学研究院有限公司 | Power grid frequency response oriented decision method and system |
CN109347120A (en) * | 2018-09-25 | 2019-02-15 | 中国电力科学研究院有限公司 | A kind of decision-making technique and system towards mains frequency response |
CN109167391A (en) * | 2018-10-11 | 2019-01-08 | 珠海吉瓦科技有限公司 | A kind of echelon battery energy storage power station energy management method and system based on set empirical mode decomposition |
CN109412192B (en) * | 2018-12-07 | 2020-11-13 | 国网江苏省电力有限公司常州供电分公司 | Self-energy-storage multi-end back-to-back flexible straight device operation method |
CN109412192A (en) * | 2018-12-07 | 2019-03-01 | 国网江苏省电力有限公司常州供电分公司 | From the back-to-back soft straight device operation method of energy storage multiterminal |
CN110086181B (en) * | 2019-05-31 | 2020-08-11 | 浙江大学 | Power self-adaptive cooperative control method in electricity-heat combined supply micro-grid during off-grid operation |
CN110086181A (en) * | 2019-05-31 | 2019-08-02 | 浙江大学 | Power adaptive cooperative control method in electric-thermal alliance micro-capacitance sensor when a kind of off-grid operation |
CN110311396B (en) * | 2019-07-30 | 2024-03-19 | 太原理工大学 | Hybrid energy storage capacity optimization configuration method for AC/DC hybrid micro-grid |
CN110311396A (en) * | 2019-07-30 | 2019-10-08 | 太原理工大学 | A kind of alternating current-direct current mixing micro-capacitance sensor hybrid energy-storing capacity configuration optimizing method |
CN110535187A (en) * | 2019-08-21 | 2019-12-03 | 国网山东省电力公司济宁供电公司 | A kind of the energy dispatching method and system of the composite energy storage capacity of active distribution network |
CN110752614B (en) * | 2019-11-07 | 2021-05-07 | 山东大学 | Energy storage system control method and system |
CN110752614A (en) * | 2019-11-07 | 2020-02-04 | 山东大学 | Energy storage system control method and system |
CN110797896A (en) * | 2019-11-28 | 2020-02-14 | 国网江苏电力设计咨询有限公司 | Hybrid energy storage capacity optimal configuration method in micro-grid wind and solar energy storage system |
CN111146793A (en) * | 2019-12-18 | 2020-05-12 | 济南大学 | Photovoltaic-energy storage system capacity optimization design method and system based on power feature extraction |
CN111146793B (en) * | 2019-12-18 | 2021-11-26 | 济南大学 | Photovoltaic-energy storage system capacity optimization design method and system based on power feature extraction |
CN111404199A (en) * | 2020-04-22 | 2020-07-10 | 江苏东源电器集团股份有限公司 | Micro-grid hybrid energy storage system based on flywheel energy storage and control method |
CN111654028A (en) * | 2020-06-15 | 2020-09-11 | 湘潭大学 | Final benefit hybrid energy storage optimized capacity allocation ratio based on optimization algorithm |
CN111754361A (en) * | 2020-06-29 | 2020-10-09 | 国网山西省电力公司电力科学研究院 | Energy storage capacity optimal configuration method and computing device of wind-storage combined frequency modulation system |
CN111754361B (en) * | 2020-06-29 | 2022-05-03 | 国网山西省电力公司电力科学研究院 | Energy storage capacity optimal configuration method and computing device of wind-storage combined frequency modulation system |
CN111900712A (en) * | 2020-08-31 | 2020-11-06 | 北方民族大学 | Direct-current micro-grid energy balance control method and system based on hybrid energy storage |
CN112182975A (en) * | 2020-10-10 | 2021-01-05 | 国网江苏省电力有限公司扬州供电分公司 | Power frequency domain decomposition optimization method based on empirical mode decomposition |
CN112260296A (en) * | 2020-10-10 | 2021-01-22 | 国网江苏省电力有限公司扬州供电分公司 | Hybrid energy storage system optimization method and device for stabilizing intermittent load |
CN112636367A (en) * | 2020-12-09 | 2021-04-09 | 南京工程学院 | Hybrid energy storage capacity optimal configuration method for stabilizing wind power fluctuation |
CN113690941A (en) * | 2021-03-24 | 2021-11-23 | 清华大学 | Method, device, equipment and storage medium for configuring capacity of optical storage and charging microgrid |
CN113690941B (en) * | 2021-03-24 | 2024-07-16 | 清华大学 | Method, device, equipment and storage medium for configuring capacity of optical storage and filling micro-grid |
CN113629736A (en) * | 2021-08-12 | 2021-11-09 | 国网江苏省电力有限公司常州供电分公司 | Intraday rolling optimization method based on power distribution network hydrogen energy storage system |
CN113762626A (en) * | 2021-09-10 | 2021-12-07 | 国网江苏省电力有限公司扬州供电分公司 | Hybrid energy storage system capacity optimization method considering economy |
CN113822707A (en) * | 2021-09-10 | 2021-12-21 | 国网冀北电力有限公司电力科学研究院 | Output decision method and device for power market, computer equipment and storage medium |
CN113822707B (en) * | 2021-09-10 | 2024-04-30 | 国网冀北电力有限公司电力科学研究院 | Output decision method and device for electric power market, computer equipment and storage medium |
CN113904352A (en) * | 2021-10-11 | 2022-01-07 | 国网福建省电力有限公司 | Power distribution optimization method and terminal for hybrid energy storage system |
CN114142473A (en) * | 2021-12-06 | 2022-03-04 | 浙江华云电力工程设计咨询有限公司 | Micro-grid configuration method and system under condition of unknown wind and solar capacity |
CN114142473B (en) * | 2021-12-06 | 2024-02-02 | 浙江华云电力工程设计咨询有限公司 | Micro-grid configuration method and system under condition of unknown wind-solar capacity |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106998072A (en) | A kind of mixed energy storage system capacity configuration optimizing method for optimizing operation towards power distribution network | |
Wong et al. | Review on the optimal placement, sizing and control of an energy storage system in the distribution network | |
Falahati et al. | A new smart charging method for EVs for frequency control of smart grid | |
Brekken et al. | Optimal energy storage sizing and control for wind power applications | |
CN109301853A (en) | A kind of micro-capacitance sensor Multiple Time Scales energy management method for stabilizing power swing | |
Sarvi et al. | An optimized fuzzy logic controller by water cycle algorithm for power management of stand-alone hybrid green power generation | |
Sujatha et al. | Enhancement of PQ in grid connected PV system using hybrid technique | |
Bogaraj et al. | A novel energy management scheme using ANFIS for independent microgrid | |
Dehbozorgi et al. | Line loss reduction and voltage profile improvement in radial distribution networks using battery energy storage system | |
CN109586325A (en) | A kind of new energy energy storage Optimal Configuration Method | |
CN104734166A (en) | Hybrid energy storage system and wind power generation power smooth control method | |
CN104537258A (en) | Cone optimization modeling method for allowing distributed stored energy to participate in running adjustment of active power distribution network | |
CN105186499A (en) | Multi-target probabilistically optimal power flow fuzzy modelling and solving method for power distribution network | |
Akshya et al. | Grid integration for electric vehicle and photovoltaic panel for a smart home | |
CN106058900A (en) | Hybrid energy storage capacity optimization configuration method in grid connected wind-solar generation | |
Justin et al. | Power quality improvement for vehicle‐to‐grid and grid‐to‐vehicle technology in a microgrid | |
Grisales-Noreña et al. | Optimal operation of battery storage systems in standalone and grid-connected DC microgrids using parallel metaheuristic optimization algorithms | |
Xu et al. | Dynamic reconfiguration of distribution network containing distributed generation | |
Ma et al. | Coordination of generation and transmission planning for power system with large wind farms | |
Bonì et al. | A simulation analysis for assessing the reliability of AC/DC hybrid microgrids–Part II: Port area and residential area | |
Alayi et al. | Optimization of Renewable energy consumption in charging electric vehicles using intelligent algorithms | |
Liu et al. | Multi-objective optimization dispatch of PV-MG considering demand response actions | |
Krim et al. | Joint optimisation of sizing and fuzzy logic power management of a hybrid storage system considering economic reliability indices | |
CN114759616A (en) | Micro-grid robust optimization scheduling method considering characteristics of power electronic devices | |
Gong et al. | Economic dispatching strategy of double lead-acid battery packs considering various factors |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170801 |
|
RJ01 | Rejection of invention patent application after publication |