CN105260618B - A kind of energy-storage system typical curve method for digging based on cloud model - Google Patents
A kind of energy-storage system typical curve method for digging based on cloud model Download PDFInfo
- Publication number
- CN105260618B CN105260618B CN201510741815.XA CN201510741815A CN105260618B CN 105260618 B CN105260618 B CN 105260618B CN 201510741815 A CN201510741815 A CN 201510741815A CN 105260618 B CN105260618 B CN 105260618B
- Authority
- CN
- China
- Prior art keywords
- mrow
- energy
- storage system
- charge
- msub
- 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.)
- Expired - Fee Related
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention belongs to power system energy-storage system technical field, more particularly to a kind of energy-storage system typical curve method for digging based on cloud model, selectes the determination application scenarios of energy-storage system and suitable control strategy, statistics energy-storage system longitudinal direction sequential power.The frequency distribution of energy-storage system longitudinal direction sequential performance number is decomposed into several varigrained cloud models, the expectation weighted sum to normal cloud model group, provides energy-storage system typical case's operation curve.The present invention can greatly reduce amount of calculation, be easy to that energy-storage system charge and discharge power is integrated with understanding cognition;Backward cloud generator in cloud model can be converted to precise information the cloud model represented with numerical characteristic, realize extraction and conversion of the numerical value to concept, that is, realize that energy-storage system charge and discharge power data excavates energy-storage system typical case's operation curve.
Description
Technical field
The invention belongs to power system energy-storage system technical field, more particularly to a kind of energy-storage system allusion quotation based on cloud model
Type curve method for digging.
Background technology
Fluctuation, intermittence and the unpredictability of renewable energy power generation can give the operation of existing power system to bring one
Fixed impact.Energy-storage system has quick response and dynamic regulation ability, can effectively improve the friendly of New-energy power system.
The research for energy-storage system charge and discharge power all concentrates on optimal control and capacity configuration etc. at present, for allusion quotation
The excacation concern of type operation curve is less, there is not clear enough to energy-storage system charge and discharge power overall cognitive, and calculates
Measure the problems such as huge.Excavate and determine that energy-storage system typical case's operation curve under application scenarios is that energy-storage system operation characteristic is carried
Take and summarize, it is possible to achieve for the overall cognitive of energy-storage system charge and discharge power situation, and then aid in solving asking for capacity configuration
Topic.The typical operation curve extracted is utilized to save amount of calculation with amount of compressed data.
The content of the invention
The problem of existing for the research of current energy-storage system charge and discharge power, the present invention proposes a kind of storage based on cloud model
Energy system typical curve method for digging, including:
Step 1, the application scenarios and control strategy for determining energy-storage system, the charge and discharge power to energy-storage system carry out data and adopted
Collection, note jth day k moment energy-storage system charge and discharges performance number is Pjk, j=1,2 ..., n, n is to count total number of days, k=1,2 ..., s, s
For daily sampled point number, k moment charge and discharge power maximum is P in n dayskmax, then mark of the energy-storage system charge and discharge power at the k moment
The one value is:Pjk_pu=Pjk/Pkmax, Pjk_pu∈[0,1];
Step 2, the perunit charge and discharge performance number at different moments for counting n days are simultaneously stored in matrix Ms*n=[p1,p2,…,pk,…,
ps]T, pkIt is the perunit charge and discharge vector power at k moment in energy-storage system n days, pk=[P1k_pu,P2k_pu,…,Pnk_pu];
Step 3, statistics perunit charge and discharge vector power pkFrequency distribution, obtain pkCurve of frequency distribution f (x)k, and root
According to curve of frequency distribution f (x)kPeak FkjCorresponding abscissa obtains the expectation E of cloud modelxki, i.e. the position of centre of gravity of cloud;i
=1,2 ..., m, initial value i=1;
Step 4, the expectation E with cloud modelxkiCentered on, first trough position of the right and left appearance is found out, is divided
X is not designated as itleftWith xright, it is determined that the water dust scope d of first cloud modeli=min (Exki-xleft,xright-Exki);Utilize
Reverse cloud algorithm without degree of certainty information calculates entropy ETkiWith super entropy Heki, then calculate cloud model data distribution function
fcloud-ki;
Step 5, from curve of frequency distribution f (x)kIn subtract cloud model distribution function fcloud-ki, obtain new frequency distribution
Curve f ' (x)k, and find f ' (x)kPeak Fki+1;Calculate peak FkiWith peak Fki+1Difference Distanceki=Fki-
Fki+1;
If step 6, DistancekiMore than or equal to threshold epsilon, then peak F is found outki+1Corresponding abscissa Exki+1, then i=i+
1 and 4~step 5 of repeat step, until DistancekmLess than threshold epsilon, and obtain corresponding abscissa Exki;If
DistancekiLess than threshold epsilon, then k=k+1, that is, enter the normal cloud model group for finding subsequent time;
Step 7, the representative charge and discharge performance number vector at k moment is obtained according to step 4~step 6:(Exk1,Exk2,…,
Exki,…,Exkm);Energy-storage system is as follows in the computational methods of the typical charge and discharge performance number at k moment:
Step 8, k+1 is calculated successively, k+2 ..., energy-storage system typical case's charge and discharge performance number at s moment, obtain energy-storage system allusion quotation
Type operation curve P=[P1, P2..., Ps]。
Reverse cloud algorithm in the step 4 is:
Step 401, by scope in (Exk1-d1, Exk1+d1) in power data be designated as xk1i, i=1,2 ..., T;According to sample
This xk1iSample is calculated it is expected:
Step 402, calculate sample variance:
Step 403, calculate entropy:
Step 404, calculate super entropy:
The data distribution function computational methods of cloud model in the step 4 are:
Beneficial effects of the present invention are as follows:Applied to peak load shifting, the plan of tracking wind-powered electricity generation is contributed, and stabilizes wind-powered electricity generation or photovoltaic
The application scenarios such as the power swing of generating, amount of calculation can be greatly reduced, energy-storage system charge and discharge power can be integrated with understanding cognition;
Backward cloud generator in cloud model can be converted to precise information the cloud model represented with numerical characteristic C (Ex, En He),
Extraction and conversion of the numerical value to concept are realized, that is, it is bent to realize that energy-storage system charge and discharge power data excavates energy-storage system typical case's operation
Line.
Brief description of the drawings
Fig. 1 is the energy-storage system typical curve mining algorithm flow chart based on cloud model;
Energy-storage system control strategy figure when Fig. 2 is the power swing for stabilizing photovoltaic plant;
Fig. 3 is the energy-storage system typical case's operation curve figure excavated in embodiment;
Fig. 4 is energy-storage system power distribution and nonparametric probability figure in embodiment.
Embodiment
Below in conjunction with the accompanying drawings, embodiment is elaborated.
The present invention proposes a kind of energy-storage system typical curve method for digging based on cloud model, as shown in figure 1, including:
Step 1, set power swing of the application scenarios of energy-storage system to stabilize photovoltaic plant.Setting time span is 1
Year.Selection control strategy is first-order low-pass ripple.As shown in Fig. 2 the light stores up association system mainly by photovoltaic plant, energy-storage system
Formed with control system.PsFor photovoltaic plant power output, PbIt is energy-storage system power output, PwIt is that light storage association system is always defeated
Go out power.Their existing relations are shown below (after discretization):
Pout(k)=α Pout(k-1)+(1-α)·Ps(k)
Ps(k)+Pb(k)=Pout(k)
Wherein, k walks for time discretization, and the corresponding moment is tk=t0+ k Δs t, t0For initial time, Δ t is when sampling
Between be spaced.τ is the time constant of low pass filter, and α is the constant relevant with τ, referred to as filter constant.
According to above-mentioned control strategy and equation, the energy-storage system charge and discharge power of 1 year is calculated.
Charge and discharge power to energy-storage system carries out data acquisition, and note jth day k moment energy-storage system charge and discharges performance number is Pjk, j
=1,2 ..., n, n is count total number of days, and k=1,2 ..., s, s is daily sampled point number, and k moment charge and discharges power is maximum in n days
It is worth for Pkmax, then perunit value of the energy-storage system charge and discharge power at the k moment be:Pjk_pu=Pjk/Pkmax, Pjk_pu∈[0,1];
Step 2, the perunit charge and discharge performance number at different moments for counting n days are simultaneously stored in matrix Ms*n=[p1,p2,…,pk,…,
ps]T, pkIt is the perunit charge and discharge vector power at k moment in energy-storage system n days, pk=[P1k_pu,P2k_pu,…,Pnk_pu];
Step 3, statistics perunit charge and discharge vector power pkFrequency distribution, obtain pkCurve of frequency distribution f (x)k, and root
According to curve of frequency distribution f (x)kPeak FkjCorresponding abscissa obtains the expectation E of cloud modelxki, i.e. the position of centre of gravity of cloud;i
=1,2 ..., m, initial value i=1;
Step 4, the expectation E with cloud modelxkiCentered on, first trough position of the right and left appearance is found out, is divided
X is not designated as itleftWith xright, it is determined that the water dust scope d of first cloud modeli=min (Exki-xleft,xright-Exki);Utilize
Reverse cloud algorithm without degree of certainty information calculates entropy ETkiWith super entropy Heki, then calculate cloud model data distribution function
fcloud-ki;
Step 5, from curve of frequency distribution f (x)kIn subtract cloud model distribution function fcloud-ki, obtain new frequency distribution
Curve f ' (x)k, and find f ' (x)kPeak Fki+1;Calculate peak FkiWith peak Fki+1Difference Distanceki=Fki-
Fki+1;
If step 6, DistancekiMore than or equal to threshold epsilon, then peak F is found outki+1Corresponding abscissa Exki+1, then i=i+
1 and 4~step 5 of repeat step, until DistancekmLess than threshold epsilon, and obtain corresponding abscissa Exki;If
DistancekiLess than threshold epsilon, then k=k+1, that is, enter the normal cloud model group for finding subsequent time;
Step 7, the representative charge and discharge performance number vector at k moment is obtained according to step 4~step 6:(Exk1,Exk2,…,
Exki,…,Exkm);Energy-storage system is as follows in the computational methods of the typical charge and discharge performance number at k moment:
Step 8, k+1 is calculated successively, k+2 ..., energy-storage system typical case's charge and discharge performance number at s moment, obtain energy-storage system allusion quotation
Type operation curve P=[P1, P2..., Ps]。
Reverse cloud algorithm in the step 4 is:
Step 401, by scope in (Exk1-d1, Exk1+d1) in power data be designated as xk1i, i=1,2 ..., T;According to sample
This xk1iSample is calculated it is expected:
Step 402, calculate sample variance:
Step 403, calculate entropy:
Step 404, calculate super entropy:
The data distribution function computational methods of cloud model in the step 4 are:
Energy-storage system typical case's operation curve is excavated according to the algorithm proposed, as shown in Figure 3.Its charge and discharge power is united
Meter with nonparametric probability be fitted, and with 6 kinds of typical weathers (be sunny, fine (thin cloud) by typical weather type set, it is fine
Turn cloudy, cloudy turn to fine, cloudy, 6 kinds of overcast and rainy snow etc., and each typical weather type corresponds to a kind of photovoltaic power situation, each
Kind photovoltaic power power situation corresponds to a kind of energy-storage system charge and discharge power features.) under energy-storage system charge-discharge electric power be distributed it is non-
The matched curve of parameter Density Estimator contrasts, as shown in Figure 4.From fig. 4, it can be seen that the norm of nonparametric kernel density of typical operation curve
Evaluation fitting curve is bonded preferably with energy storage charge and discharge power matched curve under 6 kinds of typical weathers.Typical case's fortune shown in this explanation Fig. 3
Row curve has the feature of energy storage charge and discharge power curve under 6 kinds of typical weathers concurrently, can characterize energy-storage system and stabilize photovoltaic plant
Power feature under power swing application scenarios.
This embodiment is only the preferable embodiment of the present invention, but protection scope of the present invention is not limited thereto,
Any one skilled in the art the invention discloses technical scope in, the change or replacement that can readily occur in,
It should all be included within the scope of the present invention.Therefore, protection scope of the present invention should be with scope of the claims
It is defined.
Claims (3)
- A kind of 1. energy-storage system typical curve method for digging based on cloud model, it is characterised in that including:Step 1, the application scenarios and control strategy for determining energy-storage system, the charge and discharge power to energy-storage system carry out data acquisition, Remember that jth day k moment energy-storage system charge and discharges performance number is Pjk, j=1,2 ..., n, to count total number of days, k=1,2 ..., s, s is n Daily sampled point number, k moment charge and discharge power maximum is P in n dayskmax, then perunit of the energy-storage system charge and discharge power at the k moment It is worth and is:Pjk_pu=Pjk/Pkmax, Pjk_pu∈[0,1];Step 2, the perunit charge and discharge performance number at different moments for counting n days are simultaneously stored in matrix Ms*n=[p1,p2,…,pk,…,ps]T, pkIt is the perunit charge and discharge vector power at k moment in energy-storage system n days, pk=[P1k_pu,P2k_pu,…,Pnk_pu];Step 3, statistics perunit charge and discharge vector power pkFrequency distribution, obtain pkCurve of frequency distribution f (x)k, and according to frequency Rate distribution curve f (x)kPeak FkiCorresponding abscissa obtains the expectation E of cloud modelxki, i.e. the position of centre of gravity of cloud;I=1, 2 ..., m, initial value i=1;Step 4, the expectation E with cloud modelxkiCentered on, first trough position of the right and left appearance is found out, is remembered respectively For xleftWith xright, it is determined that the water dust scope d of first cloud modeli=min (Exki-xleft,xright-Exki);Using need not The reverse cloud algorithm of degree of certainty information calculates entropy ETkiWith super entropy Heki, then calculate cloud model data distribution function fcloud-ki;Step 5, from curve of frequency distribution f (x)kIn subtract cloud model distribution function fcloud-ki, obtain new curve of frequency distribution f’(x)k, and find f ' (x)kPeak Fk(i+1);Calculate peak FkiWith peak Fk(i+1)Difference Distanceki=Fki- Fk(i+1);If step 6, DistancekiMore than or equal to threshold epsilon, then peak F is found outk(i+1)Corresponding abscissa Exk(i+1), then i=i+1 And 4~step 5 of repeat step, until DistancekiLess than threshold epsilon, and obtain corresponding abscissa Exki;If DistancekiLess than threshold epsilon, then k=k+1, that is, enter the normal cloud model group for finding subsequent time;Step 7, the representative charge and discharge performance number vector at k moment is obtained according to step 4~step 6:(Exk1,Exk2,…,Exki,…, Exkm);Energy-storage system is as follows in the computational methods of the typical charge and discharge performance number at k moment:<mrow> <msub> <mi>P</mi> <mi>k</mi> </msub> <mo>=</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <mfrac> <msub> <mi>F</mi> <mrow> <mi>k</mi> <mi>i</mi> </mrow> </msub> <mrow> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>F</mi> <mrow> <mi>k</mi> <mi>i</mi> </mrow> </msub> </mrow> </mfrac> <mo>&CenterDot;</mo> <msub> <mi>E</mi> <mrow> <mi>x</mi> <mi>k</mi> <mi>i</mi> </mrow> </msub> </mrow>Step 8, k+1 is calculated successively, k+2 ..., energy-storage system typical case's charge and discharge performance number at s moment, obtain energy-storage system typical case's fortune Row curve P=[P1, P2..., Ps]。
- 2. method according to claim 1, it is characterised in that the reverse cloud algorithm in the step 4 is:Step 401, by scope in (Exki-d1, Exki+d1) in power data be designated as xk1i, i=1,2 ..., T;According to sample xk1i Sample is calculated it is expected:Step 402, calculate sample variance:Step 403, calculate entropy:Step 404, calculate super entropy:
- 3. method according to claim 2, it is characterised in that the data distribution function of the cloud model in the step 4 calculates Method is:<mrow> <msub> <mi>f</mi> <mrow> <mi>c</mi> <mi>l</mi> <mi>o</mi> <mi>u</mi> <mi>d</mi> <mo>-</mo> <mi>k</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mi>exp</mi> <mo>&lsqb;</mo> <mo>-</mo> <mfrac> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>k</mi> <mi>l</mi> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>E</mi> <mi>x</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <msubsup> <mi>E</mi> <mrow> <mi>T</mi> <mi>k</mi> <mi>i</mi> </mrow> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>&rsqb;</mo> <mo>.</mo> </mrow> 2
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510741815.XA CN105260618B (en) | 2015-11-04 | 2015-11-04 | A kind of energy-storage system typical curve method for digging based on cloud model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510741815.XA CN105260618B (en) | 2015-11-04 | 2015-11-04 | A kind of energy-storage system typical curve method for digging based on cloud model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105260618A CN105260618A (en) | 2016-01-20 |
CN105260618B true CN105260618B (en) | 2017-11-28 |
Family
ID=55100307
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510741815.XA Expired - Fee Related CN105260618B (en) | 2015-11-04 | 2015-11-04 | A kind of energy-storage system typical curve method for digging based on cloud model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105260618B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106779340B (en) * | 2016-12-01 | 2023-08-15 | 中国电力科学研究院 | Extraction method and evaluation system of typical working condition curve of energy storage system |
CN107194495B (en) * | 2017-04-21 | 2020-05-12 | 北京信息科技大学 | Photovoltaic power longitudinal prediction method based on historical data mining |
CN113991711B (en) * | 2021-11-16 | 2023-06-16 | 国网甘肃省电力公司电力科学研究院 | Capacity configuration method for energy storage system of photovoltaic power station |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102684222A (en) * | 2012-05-14 | 2012-09-19 | 华北电力大学 | Method for smoothly controlling wind power generation power based on energy storage technology |
CN103427429A (en) * | 2013-08-21 | 2013-12-04 | 华北电力大学 | Real-time power distribution system and method of user-side distributed energy storage device |
CN103488865A (en) * | 2013-07-04 | 2014-01-01 | 清华大学 | Energy storing device capacity computing method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013144577A2 (en) * | 2012-03-30 | 2013-10-03 | Sony Corporation | Energy storage |
-
2015
- 2015-11-04 CN CN201510741815.XA patent/CN105260618B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102684222A (en) * | 2012-05-14 | 2012-09-19 | 华北电力大学 | Method for smoothly controlling wind power generation power based on energy storage technology |
CN103488865A (en) * | 2013-07-04 | 2014-01-01 | 清华大学 | Energy storing device capacity computing method |
CN103427429A (en) * | 2013-08-21 | 2013-12-04 | 华北电力大学 | Real-time power distribution system and method of user-side distributed energy storage device |
Non-Patent Citations (1)
Title |
---|
基于负荷-光伏等效负荷曲线动态分段的配电线路联络开关优化配置;陈禹等;《电力自动化设备》;20150331;第35卷(第3期);第47-53页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105260618A (en) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Teske et al. | High renewable energy penetration scenarios and their implications for urban energy and transport systems | |
Liu et al. | Deriving optimal refill rules for multi-purpose reservoir operation | |
CN103942728B (en) | Cascade hydropower station group daily power generation plan making method | |
Rajkumar et al. | Techno-economical optimization of hybrid pv/wind/battery system using Neuro-Fuzzy | |
CN104701880B (en) | Method for calculating maximum photovoltaic capacity accepted by power grid based on peak regulation constraint | |
Capizzi et al. | Recurrent neural network-based control strategy for battery energy storage in generation systems with intermittent renewable energy sources | |
CN105260797A (en) | Microgrid energy storage power station program control method based on economical evaluation | |
CN105260618B (en) | A kind of energy-storage system typical curve method for digging based on cloud model | |
CN105244920B (en) | Consider the energy-storage system multi objective control method and its system of cell health state | |
CN103218757A (en) | Method for determining credibility of photovoltaic power capacities | |
CN106548253A (en) | Method based on the wind power prediction of nonparametric probability | |
CN104319807A (en) | Method for obtaining multiple-wind farm capacity credit based on Copula function | |
CN105071423A (en) | Energy storage system capacity allocation method taking wind limiting characteristic and economy into consideration | |
CN103986193B (en) | A kind of method that maximum wind grid connection capacity obtains | |
CN105958520B (en) | The operation control strategy of energy-storage system of accumulator in a kind of power distribution network | |
CN105515042A (en) | Photovoltaic access capacity optimization algorithm of wind-photovoltaic combination farm project | |
CN106684913A (en) | Energy storage power station tracking generation plan control system and method based on multiple agents | |
CN106096807A (en) | A kind of complementary microgrid economical operation evaluation methodology considering small power station | |
CN104066165A (en) | Wireless communication power allocation method based on energy collection mode | |
CN106779340A (en) | A kind of extracting method and its evaluation system of energy-storage system typical condition curve | |
CN104598687B (en) | The optimization construction method of the water quality monitoring photovoltaic storage battery electric power system of Miniature Buoy power source | |
CN106655266B (en) | It is a kind of access new energy area power grid can flexible modulation power configuration method | |
CN110323779B (en) | Method and system for dynamically aggregating power of distributed power generation and energy storage device | |
CN111626645A (en) | Method and system for measuring and calculating future-year energy storage configuration capacity | |
Khatibi et al. | An analysis for increasing the penetration of renewable energies by optimal sizing of pumped-storage power plants |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171128 Termination date: 20181104 |
|
CF01 | Termination of patent right due to non-payment of annual fee |