CN110084465A - Wind generator system cost/Reliability Estimation Method based on energy storage - Google Patents
Wind generator system cost/Reliability Estimation Method based on energy storage Download PDFInfo
- Publication number
- CN110084465A CN110084465A CN201910181619.XA CN201910181619A CN110084465A CN 110084465 A CN110084465 A CN 110084465A CN 201910181619 A CN201910181619 A CN 201910181619A CN 110084465 A CN110084465 A CN 110084465A
- Authority
- CN
- China
- Prior art keywords
- indicates
- power
- indicate
- load
- cost
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004146 energy storage Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims abstract description 19
- 230000005611 electricity Effects 0.000 claims abstract description 17
- 238000005457 optimization Methods 0.000 claims abstract description 9
- 238000013210 evaluation model Methods 0.000 claims abstract description 6
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 4
- 238000007906 compression Methods 0.000 claims description 40
- 230000006835 compression Effects 0.000 claims description 38
- 238000009826 distribution Methods 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000010248 power generation Methods 0.000 claims description 6
- 238000012423 maintenance Methods 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 description 9
- 238000011160 research Methods 0.000 description 6
- 239000007789 gas Substances 0.000 description 5
- 238000007600 charging Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000004140 cleaning Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000002803 fossil fuel Substances 0.000 description 2
- 230000035800 maturation Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000012552 review Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000556 factor analysis Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000005338 heat storage Methods 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Biodiversity & Conservation Biology (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Wind Motors (AREA)
Abstract
The invention discloses a kind of wind generator system cost/Reliability Estimation Method based on energy storage, comprising: step 1: the operation of AA-CAES energy-accumulating power station and wind generator system generated output model are established.Step 2: cost/Reliability Evaluation Model is established.Step 3: objective function is established according to step 1 and step 2, optimization aim is to obtain fan capacity when the sum of system year cost and load power-off damages are minimum value.The present invention is using power-off damages as the assessment object of system reliability, AA-CAES energy-accumulating power station charge-discharge electric power is optimized to balance the randomness of wind-power electricity generation using dynamic programming algorithm, is influenced caused by the cost and reliability of system by changing power system capacity analysis different capabilities.
Description
Technical field
The present invention relates to Operation of Electric Systems fields, are related to a kind of wind generator system based on adiabatic compression air energy storage
Cost/Reliability Estimation Method.
Technical background
In recent years, energy and environment problem is increasingly severe, and countries in the world have started to focus and energy-saving and emission-reduction and renewable energy
The development and utilization in source.Factors shadows such as randomness, system flexible modulation power supply however the consumption of renewable energy power is contributed by it
It rings[1].The appearance of extensive energy storage technology be realize smooth renewable energy power output randomness, improve flexible peaking power source and
One of the major measure of the comprehensive consumption of renewable energy multipotency manifold formula [2-3].Wherein, water-storage, battery and compression are empty
Gas energy storage (Compressed air energy storage, CAES) is three kinds of energy storage skills suitable for storing large-scale electric energy
Art [4].Water-storage technology maturation, efficiency are higher, cost is relatively low, still, since it is limited by landform and water resources condition,
Often far from load, electrical energy transportation loss is big for addressing.Batteries to store energy high conversion efficiency, response are rapid, but by cost, service life
And the problems such as environment, influences, large-scale application is restricted.Traditional power station CAES is because capacity is big, the service life is long, starting is fast, response
Flexibly the advantages that, has much application prospect, CAES technology be considered as at present most the extensive energy storage technology of development potentiality it
One, have received widespread attention [5].
Twentieth century the seventies and eighties, international fuel price are substantially increased, under the low background of gas peak regulation unit economy,
America and Europe takes the lead in having carried out the technical research of CAES, successively builds up German Huntorf (1978) and U.S. McIntosh (1991
Year) two comercial operation power stations, it is used to peak-load regulating and power backup (black starting-up), cycle efficieny is much higher than gas peak regulation machine
Group, so far operational excellence [6-7].However, the tradition power station CAES is primarily present following two drawbacks: 1) fossil fuel is depended on,
Application in the deficient area of natural gas resource and current low-carbon energy system is limited.2) running efficiency of system is lower, it is difficult to aobvious
Write lifting system performance driving economy.With the maturation application of heat-storage technology, scholars propose advanced adiabatic compression air energy storage
(Advanced adiabatic compressed air energy storage, AA-CAES) this concept solves above-mentioned two
Item drawback.AA-CAES is a kind of by recycling and reusing air compression thermal energy, abandons conventional CAES technology fossil fuel and mends right ring
Section improves the cleaning energy storage technology [8] of system circulation efficiency.Currently, the power station AA-CAES demonstration project ADELE is German emerging
It builds, power station rated power 90MW, rated capacity 360MWh, target circulation efficiency 70% [9].In China, in May, 2017,
National Energy Board's approval project verification seat of honour 50MW grades of AA-CAES national level energy storage demonstrative project --- Community of Jin Tan County salt cave compressed air
Energy-storing and power-generating system [10];In July in the same year, reply is set up the project, and China is first to be shown by " internet+" wisdom energy of core of AA-CAES
Typical term mesh --- high ditch cable base intelligent micro-grid of doing nothing;In March, 2018, National Energy Board is " working energy refers within 2018
Lead opinion " in point out, actively push forward Community of Jin Tan County compressed-air energy storage project, research promotes 100MW compressed-air energy storage power station.
Cost/the Reliability Analysis Research currently distributed rationally to consideration AA-CAES energy-accumulating power station is blank.In this background
Under, the present invention has developed to be planned based on AA-CAES station capacity in cost/reliability assessment system.First, it is contemplated that wind-force
The randomness and the power station AA-CAES operation characteristic of power generation, the operation of the building power station AA-CAES and wind generator system generated output mould
Type is emulated using operating condition of the Monte Carlo Method to the power station AA-CAES, wind-driven generator.Secondly, user is made
For market element, the damages of interruptible load load is calculated, and is advised with powering off the minimum target of damages using dynamic
The method of drawing optimizes the power station AA-CAES charge-discharge electric power, establishes cost/Reliability Evaluation Model.Finally by simulation analysis AA-CAES
The influence that station capacity generates system economy and power supply reliability.
Bibliography
[1] Shu Yinbiao, Zhang Zhigang, Guo Jianbo are waited in the consumption critical factor analysis of new energy and solution research [J]
State's electrical engineering journal, 2017,37 (1): 1-8.
[2] Hu Juan, Yang Shuili, Hou Chaoyong, the As-Is analysis for waiting scale energy storage technology demoncal ration to apply and enlightenment
[J] electric power network technique, 2015,39 (4): 879-885.
[3] Madaeni S H, Sioshansi R, Denholm P.How Thermal Energy Storage
Enhances the Economic Viability of Concentrating Solar Power[J].Proceedings
Of the IEEE, 2012,100 (2): 335-347.
[4] Xu Yujie, Chen Haisheng, Liu Jia wait the compressed-air energy storage and generating integrated system characteristic point of wind light mutual complementing
Analyse [J] Proceedings of the CSEE, 2012,32 (20): 88-95.
[5] Liu Chang, Xu Yujie, Hu Shan wait compressed-air energy storage power station technology economic analysis [J] energy storage science and skill
Art, 2015,4 (2): 158-168.
[6]International Energy Agency.Prospects for Large-Scale Energy
Storage in Decarbonised Power Grids [R] .Paris:International Energy Agency,
2009.
[7] D í az-Gonz á lez F, Sumper A, Gomis-Bellmunt O, et al.A review of energy
storage technologies for wind power applications[J].Renewable&Sustainable
Energy Reviews, 2012,16 (4): 2154-2171.
[8] Jubeh N M, Najjar Y S H.Green solution for power generation by
Adoption of adiabatic CAES system [J] .Applied Thermal Engineering, 2012,44 (44):
85-89.
[9] Luo X, Wang J, Krupke C, et al.Modelling study, efficiency analysis
and optimisation of large-scale Adiabatic Compressed Air Energy Storage
Systems with low-temperature thermal storage [J] .Applied Energy, 2016,162:589-
600.
[10] Mei Shengwei, public cyclopentadienyl fine jade, the state of Qin is good, wait advanced adiabatic compression air energy storage technology of the based on the gas storage of salt cave and
Application prospect [J] electric power network technique, 2017,41 (10): 3392-3399.
Summary of the invention
Problem to be solved by this invention be consider running of wind generating set characteristic, AA-CAES charge-discharge electric power and power-off at
Under the premise of this, a kind of wind generator system cost/Reliability Evaluation based on advanced adiabatic compression air energy storage is proposed
Method.
The technical solution adopted by the present invention to solve the technical problems is: the wind-force based on advanced adiabatic compression air energy storage
Electricity generation system cost/Reliability Estimation Method.Include the following steps:
Step 1: the advanced adiabatic compression air power station energy storage AA-CAES operation and wind generator system generated output are established
Model.
Formula (1) indicates the constraint of the power station AA-CAES compression horsepower.Wherein, Pc,tIndicate the compression horsepower of period t;ηcIndicate pressure
Compression process efficiency;Indicate that period t flows into the flow of compressor;γ indicates air specific heat ratio;RgIndicate ideal gas constant;nc
Indicate compressor series;Tc,l,inWithRespectively indicate the air themperature into l grades of compressors and final compressor;βc,lWithRespectively indicate l grades of compressors specified compression ratio and final compressor period t compression ratio, with air storage chamber air pressure
It is related.
Pc,minνc,t≤Pc,t≤Pc,maxνc,t (2)
Formula (2) indicates the constraint of compression horsepower bound.Wherein, Pc,minAnd Pc,maxRespectively indicate the bound of compression horsepower;
νc,tFor binary variable, for indicating whether the power station AA-CAES is in compression condition, when the power station AA-CAES is in compression condition
When, νc,t=1, on the contrary νc,t=0.
Formula (3) indicates the constraint of expanding machine generated output.Wherein, Pg,tIndicate the generated output of period t;ηgExpression generated electricity
Journey efficiency;Indicate that period t flows into the flow of expanding machine;Tg,j,in,tIndicate that the t period enters the air themperature of j-th stage expanding machine,
It is related to air storage chamber temperature;ngIndicate expansion series;βg,jIndicate j-th stage expanding machine nominal expansion ratio.
Pg,minνg,t≤Pg,t≤Pg,maxνg,t (4)
Formula (4) indicates the constraint of generated output bound.Wherein, Pg,minAnd Pg,maxRespectively indicate the bound of generated output;
νg,tFor binary variable, for indicating whether the power station AA-CAES is in generating operation mode, when the power station AA-CAES is in generating operation mode
When, νg,t=1, on the contrary νg,t=0.
νc,tνg,t=0 (5)
Formula (5) indicates the constraint of the power station AA-CAES operating condition.The constraint is for guaranteeing that the power station AA-CAES does not work at the same time
In compression condition and generating operation mode.
Formula (6) indicates that wind-driven generator considers the uncertainty that its generated energy changes with wind speed, its is hourly average
Wind speed is described as obeying Weibull distribution.Wherein, ν indicates mean hourly wind speed, unit m/s;C indicates annual mean wind speed;k
The distribution and shape of wind speed are described.
Formula (7) indicates the correlation for ignoring generated output between the mechanical loss and blower of blower, doubly-fed wind turbine
Generated output.Wherein, PWTIndicate the wind energy power of blower capture;ρ indicates atmospheric density, unit kg/m3;S indicates that wind wheel is swept
Wind area, unit m2;Formula (8) indicates power coefficient.
The calculation method of parameter in formula (9) expression (8).The calculation method of parameter in formula (10) expression (9).Wherein,
θpIndicate propeller pitch angle;λ indicates tip speed ratio;ω indicates blade end angular speed;R indicates blade radius.When wind speed is cut less than blower
When entering wind speed, it is 0 that blower, which issues power,;When wind speed is between incision wind speed and rated wind speed, θp=0, by formula (6)-(9)
Calculate power of fan;When wind speed is higher than rated wind speed, blower variable blade control system work makes blower issue the specified function of power
Rate;When wind speed is greater than blower cut-out wind speed, blower is out of service.
Step 2: cost/Reliability Evaluation Model is established.
The fixed cost of formula (11) expression wind generator system and AA-CAES energy-accumulating power station.Wherein, CAIndicate initial outlay
Convert the fixed cost summation of equal years value: A1、A2Respectively indicate being fixed into for wind generator system and AA-CAES energy-accumulating power station
This;α indicates Annual Percentage Rate;N indicates loan year.
CT=CA+CM (12)
Formula (12) indicates year comprehensive cost.Wherein, CMIndicate the maintenance of wind generator system and AA-CAES energy-accumulating power station at
This,.CMNumerical value is the 2% of wind generator system year fixed cost.
Formula (13) indicates that average load powers off cost.Wherein, fAVE(h, t) indicates h type load and load power-off time t's
Average load powers off cost function;fc(h, t) indicates that load of the h type load at load power-off time t powers off cost;Table
Show the peak load of h type load.
Formula (14) indicates that synthetic load powers off cost coefficient.Wherein, fcom(t) indicate that the synthesis under load power-off time t is negative
Lotus powers off cost coefficient, reflects the relationship of electric system unit quantity of electricity load power-off cost and load section power-off time t, Ltype
Indicate load type set, σhIndicate the electricity consumption ratio of load h, LhIndicate the rate of load condensate of load h.
Formula (15) indicates the power-off of i-th load at damages.Wherein, CiIndicate that i-th load powers off damages,
QEENS,hIndicate the electricity shortage amount of h type load.
Formula (16) indicates that annual load powers off damages.Wherein, N indicates always to power off number, Q in 1 yearEENSIt indicates
Total electricity shortage amount.
Step 3: establishing objective function according to step 1 and step 2, emulates wind turbine power generation function using Monte carlo algorithm
Rate obtains fan capacity when the sum of system year cost and load power-off damages are minimum value by optimization method.
Preferably, optimization method is to be filled by dynamic programming to advanced adiabatic compression air energy storage systems in step 3
Discharge power optimizes.
This patent, which is different from existing research work, following features:
1) using power-off damages as the assessment object of system reliability
2) AA-CAES energy-accumulating power station charge-discharge electric power is optimized to balance wind-power electricity generation using dynamic programming algorithm
Randomness.
3) system analyzes different capabilities to the cost of system and can by changing power system capacity under the premise of 1) 2) meeting
By being influenced caused by property
Detailed description of the invention
Fig. 1 typical case's two stages of compression double expansion AA-CAES Power station structure figure
Fig. 2 overall plan flow chart
Fig. 3 AA-CAES charge-discharge electric power optimizes dynamic programming path
Specific embodiment
A kind of wind generator system cost/Reliability Estimation Method based on energy storage of the present invention, this method is specifically such as
Under:
Step 1: the operation of the power station AA-CAES and wind generator system generated output model are established.
The structure chart of the energy storage model is as shown in Figure 1;Formula (1) indicates the constraint of the power station AA-CAES compression horsepower.Wherein, Pc,t
Indicate the compression horsepower of period t;ηcIndicate compression process efficiency;Indicate that period t flows into the flow of compressor;γ indicates air
Specific heat ratio;RgIndicate ideal gas constant;ncIndicate compressor series;Tc,l,inWithIt respectively indicates and enters l grades of compressors
With the air themperature of final compressor;βc,lWithThe specified compression ratio and final compressor for respectively indicating l grades of compressors exist
The compression ratio of period t, it is related with air storage chamber air pressure.
Pc,minνc,t≤Pc,t≤Pc,maxνc,t (2)
Formula (2) indicates the constraint of compression horsepower bound.Wherein, Pc,minAnd Pc,maxRespectively indicate the bound of compression horsepower;
νc,tFor binary variable, for indicating whether the power station AA-CAES is in compression condition, when the power station AA-CAES is in compression condition
When, νc,t=1, on the contrary νc,t=0.
Formula (3) indicates the constraint of expanding machine generated output.Wherein, Pg,tIndicate the generated output of period t;ηgExpression generated electricity
Journey efficiency;Indicate that period t flows into the flow of expanding machine;Tg,j,in,tIndicate that the t period enters the air themperature of j-th stage expanding machine,
It is related to air storage chamber temperature;ngIndicate expansion series;βg,jIndicate j-th stage expanding machine nominal expansion ratio.
Pg,minνg,t≤Pg,t≤Pg,maxνg,t (4)
Formula (4) indicates the constraint of generated output bound.Wherein, Pg,minAnd Pg,maxRespectively indicate the bound of generated output;
νg,tFor binary variable, for indicating whether the power station AA-CAES is in generating operation mode, when the power station AA-CAES is in generating operation mode
When, νg,t=1, on the contrary νg,t=0.
νc,tνg,t=0 (5)
Formula (5) indicates the constraint of the power station AA-CAES operating condition.The constraint is for guaranteeing that the power station AA-CAES does not work at the same time
In compression condition and generating operation mode.
Formula (6) is indicated by taking doubly-fed wind turbine as an example, considers the uncertainty that its generated energy changes with wind speed, it is every
The mean wind speed of hour is described as obeying Weibull distribution.Wherein, ν (m/s) indicates mean hourly wind speed;C indicates annual
Wind speed;K describes the distribution and shape of wind speed.
Formula (7) indicates the correlation for ignoring generated output between the mechanical loss and blower of blower, doubly-fed wind turbine
Generated output.Wherein, PWTIndicate the wind energy power of blower capture;ρ indicates that atmospheric density, unit are (kg/m3);S indicates wind wheel
Wind sweeping area (m2);Formula (8) indicates power coefficient.
The calculation method of parameter in formula (9) expression (8).The calculation method of parameter in formula (10) expression (9).Wherein,
θpIndicate propeller pitch angle;λ indicates tip speed ratio;ω indicates blade end angular speed;R indicates blade radius.When wind speed is cut less than blower
When entering wind speed, it is 0 that blower, which issues power,;When wind speed is between incision wind speed and rated wind speed, θp=0, by formula (6)-(9)
Calculate power of fan;When wind speed is higher than rated wind speed, blower variable blade control system work makes blower issue the specified function of power
Rate;When wind speed is greater than blower cut-out wind speed, blower is out of service.
Step 2: cost/Reliability Evaluation model is established.
Formula (11) indicate wind generator system and AA-CAES energy-accumulating power station fixed cost include design, mounting cost and
Cost of land etc..Wherein, CAIndicate initial outlay conversion to the fixed cost summation for waiting years value: A1、A2Respectively indicate wind-power electricity generation
The fixed cost of system and AA-CAES energy-accumulating power station;α indicates Annual Percentage Rate;N indicates loan year.
CT=CA+CM (12)
Formula (12) indicates year comprehensive cost.Wherein, CMIndicate the maintenance of wind generator system and AA-CAES energy-accumulating power station at
This, including equipment replacement, cleaning, lubrication etc..CMNumerical value be about wind generator system year fixed cost 2%.
Formula (13) indicates that average load powers off cost.Wherein, fAVE(h, t) indicates h type load and load power-off time t's
Average load powers off cost function;fc(h, t) indicates that load of the h type load at load power-off time t powers off cost;Table
Show the peak load of h type load.
Formula (14) indicates that synthetic load powers off cost coefficient.Wherein, fcom(t) indicate that the synthesis under load power-off time t is negative
Lotus powers off cost coefficient, reflects the relationship of electric system unit quantity of electricity load power-off cost and load section power-off time t, Ltype
Indicate load type set, σhIndicate the electricity consumption ratio of load h, LhIndicate the rate of load condensate of load h.
Formula (15) indicates the power-off of i-th load at damages.Wherein, CiIndicate that i-th load is powered off into damages,
QEENS,hIndicate the electricity shortage amount of h type load.
Formula (16) indicates that annual load powers off damages.Wherein, N indicates always to power off number, Q in 1 yearEENSIt indicates
Total electricity shortage amount.
When establishing cost/Reliability Evaluation model, need to carry out AA-CAES energy-accumulating power station charge-discharge electric power excellent
Change.In the present invention, the rated capacity of AA-CAES energy-accumulating power station and efficiency for charge-discharge use E respectivelyCAES(MWh) and ηCAESTable
Show.Here with represented by formula (17) in 1 year the minimum target of equivalent load quadratic sum to AA-CAES energy-accumulating power station per hour
Charge-discharge electric power optimize:
Wherein, PLIt (t) is the load power of t moment;PCAESIt (t) is the charging (compression) of AA-CAES energy-accumulating power station or electric discharge
(expansion) power, charging are positive, and electric discharge is negative.Battery charge and discharge process need to meet following constraint conditions:
The minimax capacity-constrained of formula (18-19) expression battery each period.Wherein, EmaxAnd EminIt is respectively maximum and
Minimum capacity;ECAESFor AA-CAES stored energy capacitance;Period t and t-1 period energy-accumulating power station residue is held respectively by E (t) and E (t-1)
Amount;Δ t indicates time interval (being set as 1h herein).Here battery charging and discharging strategy is optimized using dynamic programming, per small
When AA-CAES energy-accumulating power station charge and discharge process regard a state in Dynamic Programming as, original state be E (0)=
0.05ECAES, after first hour, original state becomes E (0)+P from E (0)CAES(1), the recurrence equation after t hours are as follows:
Wherein PCAESIt (t) is the charge-discharge electric power of t hours AA-CAES energy-accumulating power stations;υnFor t hours decision variables;
U (t) is by PCAES(t) and the permission decision set that determines such as capacity-constrained;υn(PCAES(t), u (t)) it is equivalent load square
With.After sequence finds the minimum value of formula (17), then the charge-discharge electric power of each period AA-CAES energy-accumulating power station of search conversely, such as
Shown in Fig. 3;
Step 3: as shown in Fig. 2, establishing objective function according to step 1 and step 2, optimization aim is to obtain system year
The minimum value of the sum of cost and load power-off damages and corresponding fan capacity.
The present invention assumes that each hour various state parameters are constant using 1 year 8736 hour as time scale, when
Certain variable is value of the variable in moment t in the value of moment t.Assuming that this paper load type only includes industrial load and breaks
Electric damages is assumed to be 3$/kWh.
Every kilowatt of construction cost of blower of the present invention is 8000$/kW.The annual interest fixed rate of interest is 0.1, and repaying is limited to 20 years in year.If wind
Blower single-machine capacity in is 5MW, and incision wind speed is 3m/s, rated wind speed 12m/s, cut-out wind speed 25m/s, and air is close
Degree is ρ=1.25kg/m3, blade radius is 37.5m, and the control system of three blades, wind-driven generator can obtain optimal wind energy
Usage factor CP, then tip speed ratio λ=8.5.The construction cost in the power station AA-CAES is 300$/kW, and the service life is 40 years.Annual fortune
Row maintenance cost is the 2% of construction cost.
To compare capacity scale to the influence difference of system performance driving economy and power supply reliability and in wind resource difference
Region, capacity scale to the influence difference of system performance driving economy and power supply reliability, the present invention be arranged following 4 scenes into
Row emulation and comparative analysis.
1 scene research of table
System AA-CAES capacity is 50MW in scene 1, and separate unit fan capacity is 5MW, and the parameter of Weibull distribution is c=
7, k=2;System AA-CAES capacity is 100MW, the Weibull distribution parameter and field of separate unit fan capacity and wind speed in scene 2
Scape 1 is identical;System AA-CAES capacity is 100MW in scene 3, and separate unit fan capacity is 5MW, and the parameter of Weibull distribution is c
=5, k=2;System AA-CAES capacity is 200MW in scene 4, the Weibull distribution parameter of separate unit fan capacity and wind speed with
Scene 1 is identical.
Show that power system capacity scale can impact cost and power supply reliability: when power system capacity scale by emulation
Bigger, power supply reliability is higher, and economy is better;In the present invention, mean hourly wind speed obeys Weibull distribution within 1 year.Work as shape
When shape parameter changes, the result of optimization is also different, i.e. the mean wind speed of system location is higher, and annual total cost is smaller,
System power supply reliability is better, and the wind generator system capacity that need to be installed is smaller.Under difference AA-CAES capacity as shown in Table 2
Total cost and optimization fan capacity can be intuitive find out conclusion.
Total cost and optimization fan capacity under 2 difference AA-CAES capacity of table
Claims (2)
1. wind generator system cost/Reliability Estimation Method based on energy storage, which is characterized in that the specific step of this method
Suddenly it is:
Step 1: the advanced adiabatic compression air power station energy storage AA-CAES operation and wind generator system generated output model are established;
Formula (1) indicates the constraint of the power station AA-CAES compression horsepower;Wherein, Pc,tIndicate the compression horsepower of period t;ηcIndicate compressed
Journey efficiency;Indicate that period t flows into the flow of compressor;γ indicates air specific heat ratio;RgIndicate ideal gas constant;ncIt indicates
Compressor series;Tc,l,inAnd Tc,nc,inRespectively indicate the air themperature into l grades of compressors and final compressor;βc,lWith
βc,nc,tRespectively indicate l grades of compressors specified compression ratio and final compressor period t compression ratio, with air storage chamber gas
It is pressed with pass;
Pc,minνc,t≤Pc,t≤Pc,maxνc,t (2)
Formula (2) indicates the constraint of compression horsepower bound;Wherein, Pc,minAnd Pc,maxRespectively indicate the bound of compression horsepower;νc,tFor
Binary variable, for indicating whether the power station AA-CAES is in compression condition, when the power station AA-CAES is in compression condition,
νc,t=1, on the contrary νc,t=0;
Formula (3) indicates the constraint of expanding machine generated output;Wherein, Pg,tIndicate the generated output of period t;ηgIndicate power generation process effect
Rate;Indicate that period t flows into the flow of expanding machine;Tg,j,in,tIndicate that the t period enters the air themperature of j-th stage expanding machine, with
Air storage chamber temperature is related;ngIndicate expansion series;βg,jIndicate j-th stage expanding machine nominal expansion ratio;
Pg,minνg,t≤Pg,t≤Pg,maxνg,t (4)
Formula (4) indicates the constraint of generated output bound;Wherein, Pg,minAnd Pg,maxRespectively indicate the bound of generated output;νg,tFor
Binary variable, for indicating whether the power station AA-CAES is in generating operation mode, when the power station AA-CAES is in generating operation mode,
νg,t=1, on the contrary νg,t=0;
νc,tνg,t=0 (5)
Formula (5) indicates the constraint of the power station AA-CAES operating condition;The constraint is being pressed for guaranteeing that the power station AA-CAES does not work at the same time
Contracting operating condition and generating operation mode;
Formula (6) indicates that wind-driven generator considers the uncertainty that its generated energy changes with wind speed, by its mean wind speed hourly
It is described as obeying Weibull distribution;Wherein, ν indicates mean hourly wind speed, unit m/s;C indicates annual mean wind speed;K description
The distribution and shape of wind speed;
Formula (7) indicates the correlation for ignoring generated output between the mechanical loss and blower of blower, the power generation of doubly-fed wind turbine
Power;Wherein, PWTIndicate the wind energy power of blower capture;ρ indicates atmospheric density, unit kg/m3;S indicates wind wheel swing flap face
Product, unit m2;Formula (8) indicates power coefficient;
The calculation method of parameter in formula (9) expression (8);The calculation method of parameter in formula (10) expression (9);Wherein, θpIt indicates
Propeller pitch angle;λ indicates tip speed ratio;ω indicates blade end angular speed;R indicates blade radius;Wind speed is cut when wind speed is less than blower
When, it is 0 that blower, which issues power,;When wind speed is between incision wind speed and rated wind speed, θp=0, wind is calculated by formula (6)-(9)
Machine power;When wind speed is higher than rated wind speed, blower variable blade control system work makes blower issue power rated power;When
When wind speed is greater than blower cut-out wind speed, blower is out of service;
Step 2: cost/Reliability Evaluation Model is established;
The fixed cost of formula (11) expression wind generator system and AA-CAES energy-accumulating power station;Wherein, CAIndicate initial outlay conversion
To the fixed cost summation of equal years value: A1、A2Respectively indicate the fixed cost of wind generator system and AA-CAES energy-accumulating power station;α
Indicate Annual Percentage Rate;N indicates loan year;
CT=CA+CM (12)
Formula (12) indicates year comprehensive cost;Wherein, CMIndicate the maintenance cost of wind generator system and AA-CAES energy-accumulating power station,;
CMNumerical value is the 2% of wind generator system year fixed cost;
Formula (13) indicates that average load powers off cost;Wherein, fAVE(h, t) expression h type load is averaged with load power-off time t's
Load powers off cost function;fc(h, t) indicates that load of the h type load at load power-off time t powers off cost;Indicate h class
The peak load of load;
Formula (14) indicates that synthetic load powers off cost coefficient;Wherein, fcom(t) indicate that the synthetic load under load power-off time t is disconnected
Electric cost coefficient reflects the relationship of electric system unit quantity of electricity load power-off cost and load section power-off time t, LtypeIt indicates
Load type set, σhIndicate the electricity consumption ratio of load h, LhIndicate the rate of load condensate of load h;
Formula (15) indicates the power-off of i-th load at damages;Wherein, CiIndicate that i-th load powers off damages, QEENS,hTable
Show the electricity shortage amount of h type load;
Formula (16) indicates that annual load powers off damages;Wherein, N indicates always to power off number, Q in 1 yearEENSIndicate total power supply
It is in shortage;
Step 3: establishing objective function according to step 1 and step 2, emulates wind turbine power generation power using Monte carlo algorithm, leads to
It crosses optimization method and obtains fan capacity when the sum of system year cost and load power-off damages are minimum value.
2. wind generator system cost/Reliability Estimation Method according to claim 1 based on energy storage, feature
Be: optimization method in step 3 carries out advanced adiabatic compression air energy storage systems charge-discharge electric power by dynamic programming excellent
Change.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910181619.XA CN110084465B (en) | 2019-03-11 | 2019-03-11 | Wind power generation system cost/power supply reliability evaluation method based on energy storage |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910181619.XA CN110084465B (en) | 2019-03-11 | 2019-03-11 | Wind power generation system cost/power supply reliability evaluation method based on energy storage |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110084465A true CN110084465A (en) | 2019-08-02 |
CN110084465B CN110084465B (en) | 2022-03-18 |
Family
ID=67412378
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910181619.XA Active CN110084465B (en) | 2019-03-11 | 2019-03-11 | Wind power generation system cost/power supply reliability evaluation method based on energy storage |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110084465B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110580394A (en) * | 2019-09-11 | 2019-12-17 | 广西大学 | Parallel multilayer Monte Carlo double-fed asynchronous fan parameter optimization method |
CN110635496A (en) * | 2019-09-30 | 2019-12-31 | 杭州电子科技大学 | Day-ahead robust scheduling method considering load fluctuation of power distribution system |
CN111666664A (en) * | 2020-05-22 | 2020-09-15 | 湖北省电力勘测设计院有限公司 | Method for evaluating confidence capacity of newly added energy supply element of electricity-gas integrated energy system |
CN112598539A (en) * | 2020-12-28 | 2021-04-02 | 江苏徐工信息技术股份有限公司 | Wind power curve optimization calculation and abnormal value detection method for wind generating set |
CN117252046A (en) * | 2023-11-20 | 2023-12-19 | 中国电力工程顾问集团有限公司 | Dynamic simulation method and device for compressed air energy storage system |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103488864A (en) * | 2013-07-04 | 2014-01-01 | 清华大学 | Evaluation method for risk of wind power, photovoltaic power generation and energy storage combined power generation system |
CN104851053A (en) * | 2015-05-14 | 2015-08-19 | 上海电力学院 | Wind-photovoltaic-energy-storage-contained method for power supply reliability evaluation method of distribution network |
CN105406470A (en) * | 2015-12-21 | 2016-03-16 | 国家电网公司 | Reliability evaluation method for active power distribution network based on switch boundary subarea division |
WO2017163561A1 (en) * | 2016-03-25 | 2017-09-28 | 株式会社日立製作所 | Operation support device and wind power system |
CN107482638A (en) * | 2017-07-21 | 2017-12-15 | 杭州电子科技大学 | Supply of cooling, heating and electrical powers type micro-capacitance sensor multiobjective Dynamic Optimization dispatching method |
CN108090666A (en) * | 2017-12-13 | 2018-05-29 | 华中科技大学 | A kind of coordinated dispatching method of power grid electric energy and spare capacity containing AA-CAES |
-
2019
- 2019-03-11 CN CN201910181619.XA patent/CN110084465B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103488864A (en) * | 2013-07-04 | 2014-01-01 | 清华大学 | Evaluation method for risk of wind power, photovoltaic power generation and energy storage combined power generation system |
CN104851053A (en) * | 2015-05-14 | 2015-08-19 | 上海电力学院 | Wind-photovoltaic-energy-storage-contained method for power supply reliability evaluation method of distribution network |
CN105406470A (en) * | 2015-12-21 | 2016-03-16 | 国家电网公司 | Reliability evaluation method for active power distribution network based on switch boundary subarea division |
WO2017163561A1 (en) * | 2016-03-25 | 2017-09-28 | 株式会社日立製作所 | Operation support device and wind power system |
CN107482638A (en) * | 2017-07-21 | 2017-12-15 | 杭州电子科技大学 | Supply of cooling, heating and electrical powers type micro-capacitance sensor multiobjective Dynamic Optimization dispatching method |
CN108090666A (en) * | 2017-12-13 | 2018-05-29 | 华中科技大学 | A kind of coordinated dispatching method of power grid electric energy and spare capacity containing AA-CAES |
Non-Patent Citations (2)
Title |
---|
RAJESH KARKI等: "A simplified wind power generation model for reliability evaluation", 《IEEE TRANSACTIONS ON ENERGY CONVERSION》 * |
刘畅 等: "压缩空气储能电站技术经济性分析", 《储能科学与技术》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110580394A (en) * | 2019-09-11 | 2019-12-17 | 广西大学 | Parallel multilayer Monte Carlo double-fed asynchronous fan parameter optimization method |
CN110580394B (en) * | 2019-09-11 | 2022-08-12 | 广西大学 | Parallel multilayer Monte Carlo double-fed asynchronous fan parameter optimization method |
CN110635496A (en) * | 2019-09-30 | 2019-12-31 | 杭州电子科技大学 | Day-ahead robust scheduling method considering load fluctuation of power distribution system |
CN110635496B (en) * | 2019-09-30 | 2021-04-13 | 杭州电子科技大学 | Day-ahead robust scheduling method considering load fluctuation of power distribution system |
CN111666664A (en) * | 2020-05-22 | 2020-09-15 | 湖北省电力勘测设计院有限公司 | Method for evaluating confidence capacity of newly added energy supply element of electricity-gas integrated energy system |
CN111666664B (en) * | 2020-05-22 | 2023-10-27 | 湖北省电力勘测设计院有限公司 | Assessment method for confidence capacity of newly-added energy element of electric-gas comprehensive energy system |
CN112598539A (en) * | 2020-12-28 | 2021-04-02 | 江苏徐工信息技术股份有限公司 | Wind power curve optimization calculation and abnormal value detection method for wind generating set |
CN112598539B (en) * | 2020-12-28 | 2024-01-30 | 徐工汉云技术股份有限公司 | Wind power curve optimization calculation and outlier detection method for wind generating set |
CN117252046A (en) * | 2023-11-20 | 2023-12-19 | 中国电力工程顾问集团有限公司 | Dynamic simulation method and device for compressed air energy storage system |
CN117252046B (en) * | 2023-11-20 | 2024-02-27 | 中国电力工程顾问集团有限公司 | Dynamic simulation method and device for compressed air energy storage system |
Also Published As
Publication number | Publication date |
---|---|
CN110084465B (en) | 2022-03-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110084465A (en) | Wind generator system cost/Reliability Estimation Method based on energy storage | |
Adefarati et al. | Energizing renewable energy systems and distribution generation | |
Deng et al. | System modeling and optimization of microgrid using genetic algorithm | |
CN111737884B (en) | Multi-target random planning method for micro-energy network containing multiple clean energy sources | |
Chen et al. | Sizing of energy storage system for microgrids | |
Jintao et al. | Optimized operation of multi-energy system in the industrial park based on integrated demand response strategy | |
CN112307603A (en) | Hybrid energy storage capacity optimal configuration method and system considering large-scale wind power access | |
CN115882483A (en) | Method for realizing optimal energy storage capacity configuration of system by using capacity elasticity | |
CN112561120B (en) | Microgrid-based optimized operation method for day-ahead market clearing system | |
Yang et al. | Coordinated Optimal Configuration Method of Hybrid Energy Storage Systems in Energy Internet System | |
Qing et al. | Trading model of virtual power plant considering shapley contribution based revenue distribution | |
CN111668882A (en) | Method and device for optimizing output of micro power supply in intelligent energy ring network | |
Si et al. | A multi-objective optimization strategy for combined heat and power systems of the energy internet | |
Chen et al. | A hybrid energy storage power system dispatch strategy for demand response | |
Mei et al. | Hierarchical optimal scheduling of regional integrated energy power system based on multi-objective particle swarm optimization algorithm | |
Kurinjimalar et al. | Analysis of Vehicle to Grid (V2G) Integration in Renewable based Micro grid | |
CN115374999B (en) | Water-electricity hydrogen production optimal configuration method suitable for start-stop characteristics of hydrogen production equipment | |
Wang et al. | Optimization Configuration for Energy Storage of Renewable Energy Power Stations Participating in Primary Frequency Regulation Service Based on Markov Decision Process | |
Tao et al. | Aggregating Energy Storage in Virtual Power Plant and Its Application in Unit Commitment | |
Tian et al. | Research on planning method of park-level integrated energy system considering multi-objective collaboration | |
Ahmed et al. | Integration of Hybrid Wind/Battery System into Diesel Power Plant | |
GAO et al. | MICROGRID MULTI-OBJECTIVE DYNAMIC OPTIMISATION SCHEDULING MODEL OF LOW POLLUTANT EMISSION. | |
Du et al. | Integrated energy micro-grid capacity optimization allocation method considering electricity-thermal-gas coupling | |
Wang et al. | An Aggregation Model of Multiple Energy Storage Units Based on Effect Approximation Method | |
Liu | Optimal Dispatching of Microgrid Considering Demand Response and Low-Carbon Hydrogen Production |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |