CN107134810B - Independent micro-energy-grid energy storage system optimal configuration solving method - Google Patents

Independent micro-energy-grid energy storage system optimal configuration solving method Download PDF

Info

Publication number
CN107134810B
CN107134810B CN201710431745.7A CN201710431745A CN107134810B CN 107134810 B CN107134810 B CN 107134810B CN 201710431745 A CN201710431745 A CN 201710431745A CN 107134810 B CN107134810 B CN 107134810B
Authority
CN
China
Prior art keywords
storage system
energy storage
power
unit
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.)
Active
Application number
CN201710431745.7A
Other languages
Chinese (zh)
Other versions
CN107134810A (en
Inventor
崔明勇
王楚通
王玉翠
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yanshan University
Original Assignee
Yanshan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yanshan University filed Critical Yanshan University
Priority to CN201710431745.7A priority Critical patent/CN107134810B/en
Publication of CN107134810A publication Critical patent/CN107134810A/en
Application granted granted Critical
Publication of CN107134810B publication Critical patent/CN107134810B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/008Circuit arrangements for ac mains or ac distribution networks involving trading of energy or energy transmission rights
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • H02J3/386
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/002Flicker reduction, e.g. compensation of flicker introduced by non-linear load
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/003Load forecast, e.g. methods or systems for forecasting future load demand
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Battery Electrode And Active Subsutance (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses an optimal configuration solving method for an independent micro energy network energy storage system, which comprises the following steps: establishing a self-sufficient probability demand of the micro energy network; establishing an energy storage system optimal configuration model of the micro energy network; aiming at an optimal configuration model of the energy storage system, estimating the service life of the storage battery by adopting a rain flow counting method to obtain an investment coefficient; and (3) adopting a proper operation strategy in the simulation operation process by predicting annual load data and wind power data, and selecting a traversal algorithm or a particle swarm algorithm according to an optimization range to seek the optimal energy storage system configuration. The invention provides self-sufficient probability requirements and pursues system stability while meeting the economy; the system stability is maintained while the economic optimization is achieved, and a balance point is found among the investment cost, the operation cost and the pollution control cost.

Description

Independent micro-energy-grid energy storage system optimal configuration solving method
Technical Field
The invention relates to the technical field of planning of independent micro energy network energy storage systems, in particular to an optimal configuration solving method for two aspects of independent micro energy network electricity energy storage and heat energy storage, and particularly relates to an optimal configuration solving method for an independent micro energy network energy storage system.
Background
With the gradual depletion of traditional fossil energy, environmental problems and the increasing severity of global warming problems, low-carbon new energy represented by wind and light is vigorously developed, and the improvement of the permeability of renewable energy of the existing power grid becomes one of important ways for solving the problems. Thus, the energy internet concept proposed by jirimm-rifugin has attracted widespread attention. The micro-energy network can be applied to micro-energy networks in areas such as factories, large buildings, urban and rural concentrated living areas, isolated islands and the like, and is one of the trends of future energy system development as an important composition form of an energy internet. The concept of the micro energy network is developed on the concept of a micro power grid, generally comprises 4 energy forms of cold, heat, electricity and gas, and all energy supply equipment in a region are uniformly integrated and scheduled by using the internet of things technology and the information technology so as to optimize energy supply for the cold, heat and electricity loads in the region and improve the energy utilization efficiency.
However, the renewable energy represented by wind and light has strong intermittent and random fluctuation, and phenomena such as wind and light abandonment are often caused. Especially during the heat supply period, the peak regulation capacity of the whole micro-energy grid electric energy can be reduced by the 'fixing the power by heat' operation mode of the cogeneration unit, and even a large amount of 'abandoned wind' is caused. To take up renewable energy and enhance the flexibility of micro energy networks, we introduce energy storage systems in micro energy networks that can be used for multi-energy storage. The invention relates to an independent micro-energy-grid energy storage system which comprises two forms of electricity storage and heat storage.
In summary, in the existing micro energy grid framework, it is necessary to invent an optimal configuration method for the energy storage system, so as to achieve the purposes of achieving optimal economy and reducing loss of the whole system under a certain stability.
Disclosure of Invention
The invention aims to provide an optimal configuration solving method for an independent micro energy network energy storage system, so that the independent micro energy network energy storage system can achieve economic optimization in the independent micro energy network, and meanwhile, the stability of the system is maintained.
Energy storage systems are divided into electrical energy storage systems (batteries) and thermal energy storage systems, and the configuration of the energy storage system includes a power configuration and a capacity configuration.
The storage battery energy storage system comprises a storage battery, a converter and other equipment, so that investment cost is settled in two forms of power and capacity respectively. The heat storage system comprises a heat storage tank, a heat conduction material and the like, so the investment cost is also settled in two forms of power and capacity.
In order to solve the technical problems, the following technical scheme is proposed to realize:
in the independent micro-energy-grid energy storage system, the configuration of the energy storage system influences the investment cost of the energy storage system, the operation cost of the whole micro-grid and the pollution control cost; the lower energy storage system can not achieve the expected economical efficiency and stability of the system, the operation cost can not be effectively reduced, and the discharged CO can not be effectively reduced2And harmful gasesThe volume content is higher; the higher energy storage system is configured, so that the investment cost is higher, and the overall maintenance cost is also relatively higher; therefore, in the method for solving the optimal configuration of the independent micro energy network energy storage system, the configuration selection of the optimal energy storage system can reach the balance among the investment cost, the operation cost and the pollution control cost, and the sum of the investment cost, the operation cost and the pollution control cost is found at a balance point, namely the configuration of the energy storage system with the minimum total cost.
The method comprises the following steps:
step 1, establishing self-sufficiency probability requirements of a micro energy network;
step 2, establishing an energy storage system optimal configuration model of the micro energy network;
step 3, aiming at the optimal configuration model of the energy storage system, estimating the service life of the storage battery by adopting a rain flow counting method to obtain an investment coefficient;
and 4, adopting a proper operation strategy in the simulation operation process by predicting the annual load data and the wind power data, and selecting a traversal algorithm or a particle swarm algorithm according to the optimization range to seek the optimal energy storage system configuration.
Further, in step 1, the self-sufficient probability requirement is modeled as:
Figure GDA0002624803370000031
Figure GDA0002624803370000032
wherein, PSSeAnd PSShThe self-sufficient probabilities of the electric load and the heat load in the micro energy network are respectively, and the delta w, the delta d and the delta h respectively meet the wind output prediction error of normal distribution.
Further, in step 2, the specific process of establishing the energy storage system optimization configuration model of the micro energy network includes the following steps:
the comprehensive total cost of investment, operation and pollution control environment and fuel cost needs to be considered in the model, and the model is established by taking the minimum total cost as an objective function:
Min(IC+OC+PC) (3)
Figure GDA0002624803370000033
Figure GDA0002624803370000041
Figure GDA0002624803370000042
wherein, IC is the investment cost of the energy storage system, OC is the operation cost of the micro energy network, PC is the pollution control cost of the micro energy network, alpha is the unit power investment coefficient of the storage battery, beta is the unit capacity investment coefficient of the storage battery, chi is the unit power investment coefficient of the heat storage system, delta is the unit capacity investment coefficient of the heat storage system,
Figure GDA0002624803370000043
is the maximum power of the storage battery,
Figure GDA0002624803370000044
is the maximum capacity of the storage battery,
Figure GDA0002624803370000045
is the maximum power of the heat storage system,
Figure GDA0002624803370000046
is the maximum capacity of the heat storage system, NT is the total number of days, NH is the total number of hours, NG is the total number of conventional thermal power units, NL is the total number of cogeneration units, PithIs the power, F, generated by a conventional unit during a certain period of timeeIs a function of the power and cost, IithIs the state index of whether the distributed power supply works, the work is 1, the non-work is 0, PlthIs produced by a certain cogeneration unit in a certain period of timePower of generation, FhIs a function of the power and cost, LlthThe index is the state index of whether the distributed power supply works, the work is 1, and the non-work is 0; SUth,SDthStarting or stopping costs, alpha, of the generator set, respectivelyKRepresenting the cost coefficient of treatment, beta, of different pollutantsKExpressing the emission coefficients of different pollutants, and NK expresses the total amount of the pollutants;
wherein, the unit capacity investment coefficient of the storage battery is as follows:
Figure GDA0002624803370000047
in the formula, CEIs the total electrical energy storage unit capacity investment cost, CmIs the sum of the repair, maintenance and equipment disposal costs per unit volume of equipment; the expenses are evenly distributed in the period life of the energy storage system, so that a unit investment coefficient in a planning period is obtained;
the coal consumption cost can be expressed as a quadratic function form of the generated power; the power cost functions of the conventional thermal power generating unit and the cogeneration unit are respectively as follows:
Figure GDA0002624803370000051
Figure DEST_PATH_FDA0002624803360000033
the constraint conditions of the energy storage system optimization configuration model of the micro energy network are as follows:
electric power balance constraint, heat supply balance constraint, wind power output constraint, unit constraint, power storage system constraint and heat storage system constraint; the unit constraints comprise unit output upper and lower limit constraints, extraction type unit heat output upper and lower limit constraints, extraction type unit net power generation output upper and lower limit constraints, total power climbing constraints of the unit and heat climbing constraints of the extraction type unit; the specific content of the constraint condition is as follows:
(1) electric power balance constraint:
Figure GDA0002624803370000053
where NR is the amount of new energy, PrthIs the power generated by the new energy source, PESSThe charging or discharging power of the energy storage system is provided, the charging is a negative value, and the discharging is a positive value; pload,thThe power required by the load for the period of time;
(2) and (3) heat supply balance constraint:
Figure GDA0002624803370000054
in the formula, hlthThe thermal power of the thermoelectric power unit l in the time interval; h ishsThe storage and heat release power of the heat storage tank in the time period; heat release is positive, heat storage is negative, hload,thThe thermal load of the system for the period of time; NL is the number of all thermoelectric units;
(3) wind power output restraint:
Figure GDA0002624803370000055
in the formula (I), the compound is shown in the specification,
Figure GDA0002624803370000061
is the rated power, v, of the wind turbineCI,vRAnd vCORespectively representing the cut-in wind speed, the rated wind speed and the cut-out wind speed of the fan, vhtIs the wind speed for a certain period of time;
(4) unit restraint:
1) and (3) restraining the upper and lower limits of the unit output:
pi,min≤pi,t≤pi,max (13)
in the formula, pi,min、pi,maxRespectively the minimum output and the maximum output of the unit under the pure condensation working condition;
2) the upper limit and the lower limit of the thermal output of the thermoelectric unit are restricted:
0≤hi,t≤hi,max (14)
in the formula, hi,maxThe maximum limit value of the thermal output of the unit i is determined, and the value is mainly determined by the capacity of the heat exchanger;
3) and (3) total power climbing constraint of the unit:
Pith-Pit(h-1)≤URi(1-yith)+Pi minyith (15)
Pith(h-1)-Pith≤DRi(1-zith)+Pi minzith (16)
in the formula, URiIs a ramp-up limit, yithIs whether the unit is started or not, Pi minIs the minimum power generation amount, DR, of the generator setiIs a ramp down limit, zithIs the amount of whether the unit is in a shutdown state or not;
4) thermal ramp restraint of the thermoelectric unit:
hi,t-hi,t-1≤Δhu,i (17)
hi,t-1-hi,t≤Δhd,i (18)
in the formula,. DELTA.hu,i、Δhd,iRespectively the maximum variation of the thermal power of the extraction type unit in unit time;
(5) and (3) electric storage system constraint:
the charging process is
Figure GDA0002624803370000062
The discharge process is
Figure GDA0002624803370000071
SOC (t) is the residual electric quantity of the energy storage system at the end of the tth period; SOC (t-1) is the residual electric quantity of the energy storage system at the end of the t-1 th time period; delta is the self-discharge rate of the energy storage system; pc、PdRespectively charging and discharging power of the energy storage system; etac、ηdRespectively the charging and discharging efficiency of the electricity storage system;
Figure GDA0002624803370000072
rated capacity of the power storage system;
(6) and (3) heat storage system constraint:
Figure GDA0002624803370000073
in the formula, HHS(t) thermal energy storage capacity for time period t; mu is heat storage and heat dissipation loss rate; qHS_ch(t)、QHS_dis(t) and ηhch、ηhdisRespectively the heat absorption and discharge power and efficiency in the time period t.
Further, in step 3, the life of the storage battery is estimated by using a rain flow counting method, which specifically comprises the following steps:
the method is mainly used in engineering, and is particularly widely applied to fatigue life calculation. The strain-time history data recording is rotated by 90 degrees, the time coordinate axis is vertical downwards, the data recording is like a series of roofs, and rainwater flows downwards along the roofs, so the method is called a rainwater flow counting method.
The service life of the energy storage battery is influenced by factors including the discharge depth, the rate performance, the charge-discharge cut-off voltage and the ambient temperature of the battery; the influence of the multiplying power performance of the battery on the service life of the battery is not considered, and the rated value is taken as the maximum power of the energy storage battery; the influence of the charge-discharge cut-off voltage of the battery on the service life of the battery is not considered, and the capacity value range of the energy storage battery is set; the ambient temperature is taken as room temperature, regardless of the effect of temperature on battery life; through simplification, the life of the energy storage battery is estimated by using the following formula:
Figure GDA0002624803370000081
T=ceil(1/365·dloss) (23)
wherein d islossThe damage rate of the service life of the electric energy storage system is one day, theta is a period coefficient, the full period is 1, and the half period is 0.5; cyciThe maximum cycle number corresponding to the ith cycle period; t is the life cycle, ceil is the ceiling function; the service life of the energy storage system can be estimated by constructing a charging and discharging curve of the energy storage system in one day.
Further, in step 4, the energy storage system configuration that is optimized according to the optimization range is selected through a traversal algorithm or a particle swarm optimization, and the specific process is as follows:
determining an adopted method according to the configuration range of the energy storage system, and if the range is smaller and the configuration types are less, adopting a traversal algorithm; if the range is large, the particle swarm algorithm is adopted to seek the optimal energy storage system configuration to respectively obtain the capacity and the power of the energy storage system and the capacity and the power of the heat storage system; when the optimizing range is smaller, selecting a configuration calculation cost from the configuration range of the energy storage system meeting the self-sufficient probability condition, judging whether the configuration used in the given range is completed or not, completing the calculation process if the traversal of all the configurations in the given range is completed, and continuing the calculation of the next configuration if the traversal is not completed; the method is comprehensive and is suitable for the optimization problem with a small range;
when the optimizing range is large, a particle swarm optimization algorithm is adopted, wherein the algorithm is an iterative-based optimization tool; the algorithm starts from a group of random solutions and searches an optimal solution through continuous iteration; in each iteration, updating the population by tracking individual extrema and global extrema; the algorithm has the advantages of simplicity in implementation, high speed in convergence, high precision and the like, and is widely applied to engineering practice.
The particle swarm algorithm comprises the following basic steps:
1) randomly initializing the position and speed of each particle in the population;
2) evaluating the fitness of each particle, storing the current position and the adaptive value of each particle in pbest of each particle, and storing the position and the adaptive value of an individual with the optimal adaptive value in all pbest in gbest;
3) by the formula (17) hi,t-hi,t-1≤Δhu,iUpdate the velocity and displacement of the particles:
Figure GDA0002624803370000091
4) for each particle, comparing its fitness value with its experienced best position, and if so, taking it as the current best position;
5) comparing the values of all pbest and gbest at present, and updating the gbest;
6) if the preset stopping condition of the operation precision or the iteration times is met, stopping the search and outputting the result, otherwise, returning to the step 3 and continuing the search.
Compared with the prior art, the invention has the following beneficial effects:
1. the invention considers the factors comprehensively and provides a method for configuring the energy local area network multi-energy storage system under the background of the energy Internet; applying a background update taking into account the lifetime of the electrical energy storage system;
2. the invention provides self-sufficient probability requirements and pursues system stability while meeting the economy; the system stability is maintained while the economic optimization is achieved, and a balance point is found among the investment cost, the operation cost and the pollution control cost;
3. in the aspect of algorithm selection, one algorithm is specifically selected according to the optimizing range, and the calculation result is accurate;
4. the invention better accords with the characteristics of the independent micro energy network, has strong practicability, achieves the economic optimization in the independent micro energy network, and simultaneously maintains the stability of the system; the purposes of realizing optimal economy and reducing loss of the whole system under certain stability are achieved.
Drawings
FIG. 1 is a basic schematic of the method of the present invention;
FIG. 2 is a basic flow diagram of the present invention;
fig. 3 is an explanatory view of the rain flow counting method.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
The invention discloses an optimal configuration solving method for an independent micro-energy-grid energy storage system, which is characterized in that a basic schematic diagram is shown in figure 1, and the configuration of the energy storage system influences three aspects of investment cost of the energy storage system, operation cost of the whole micro-grid and pollution control cost; the lower energy storage system can not achieve the expected economical efficiency and stability of the system, the operation cost can not be effectively reduced, and the discharged CO can not be effectively reduced2And some harmful gases are higher in content; the higher energy storage system is configured, so that the investment cost is higher, and the overall maintenance cost is also relatively higher; therefore, the configuration selection of the optimal energy storage system can reach balance among investment cost, operation cost and pollution control cost, and the sum of the investment cost, the operation cost and the pollution control cost is found at a balance point, namely the configuration of the energy storage system with the minimum total cost is reached;
the basic flow chart of the method of the invention is shown in fig. 2, and the steps of the method are as follows:
step 1, establishing self-sufficient probability requirements of a micro energy network, and establishing a self-sufficient probability requirement model; the self-sufficient probability demand model is as follows:
Figure GDA0002624803370000101
Figure GDA0002624803370000102
wherein, PSSeAnd PSShThe self-sufficient probabilities of the electric load and the heat load in the micro energy network are respectively, and the delta w, the delta d and the delta h respectively meet the wind output prediction error of normal distribution.
Step 2, establishing an energy storage system optimal configuration model of the micro energy network;
the method for establishing the energy storage system optimization configuration model of the micro energy network comprises the following specific processes:
the comprehensive total cost of investment, operation and pollution control environment and fuel cost needs to be considered in the model, the minimum total cost is taken as an objective function, and the model is established as follows:
Min(IC+OC+PC) (3)
Figure GDA0002624803370000111
Figure GDA0002624803370000112
Figure GDA0002624803370000113
in the formula, IC is the investment cost of the energy storage system, OC is the operation cost of the micro-grid, PC is the pollution control cost of the micro-grid, alpha is the unit power investment coefficient of the storage battery, beta is the unit capacity investment coefficient of the storage battery, chi is the unit power investment coefficient of the heat storage system, delta is the unit capacity investment coefficient of the heat storage system,
Figure GDA0002624803370000114
is the maximum power of the storage battery,
Figure GDA0002624803370000115
is the maximum capacity of the storage battery,
Figure GDA0002624803370000116
is the maximum power of the heat storage system,
Figure GDA0002624803370000117
is the maximum capacity of the heat storage system, NT is the total number of days, NH is the total number of hours, NG is the total number of conventional thermal power units, NL is the total heat and powerNumber of co-production units, PithIs the power, F, generated by a conventional unit during a certain period of timeeIs a function of the power and cost, IithIs the state index of whether the distributed power supply works, the work is 1, the non-work is 0, PlthIs the power, F, generated by a cogeneration unit during a certain period of timehIs a function of the power and cost, LlthThe index is the state index of whether the distributed power supply works, the work is 1, and the non-work is 0; SUth, SDthStarting or stopping costs, alpha, of the generator set, respectivelyKRepresenting the cost coefficient of treatment, beta, of different pollutantsKExpressing the emission coefficients of different pollutants, and NK expresses the total amount of the pollutants;
wherein the investment coefficient per unit capacity for the storage battery is:
Figure GDA0002624803370000121
in the formula, CEIs the total electrical energy storage unit capacity investment cost, CmIs the sum of the repair, maintenance and equipment disposal costs per unit volume of equipment; the expenses are evenly distributed in the period life of the energy storage system, so that a unit investment coefficient in a planning period is obtained;
the coal consumption cost can be expressed as a quadratic function form of the generated power; the power cost functions of the conventional thermal power generating unit and the cogeneration unit are respectively as follows:
Figure GDA0002624803370000122
Figure 678112DEST_PATH_FDA0002624803360000033
the constraint conditions of the energy storage system optimization configuration model of the micro energy network are as follows:
electric power balance constraint, heat supply balance constraint, wind power output constraint, unit constraint, power storage system constraint and heat storage system constraint; the unit constraints comprise unit output upper and lower limit constraints, extraction type unit heat output upper and lower limit constraints, extraction type unit net power generation output upper and lower limit constraints, total power climbing constraints of the unit and heat climbing constraints of the extraction type unit; the specific content of the constraint condition is as follows:
(1) electric power balance constraint:
Figure GDA0002624803370000124
where NR is the amount of new energy, PrthIs the power generated by the new energy source, PESSThe charging or discharging power of the energy storage system is provided, the charging is a negative value, and the discharging is a positive value; pload,thThe power required by the load for the period of time;
(2) and (3) heat supply balance constraint:
Figure GDA0002624803370000125
in the formula, hlthThe thermal power of the thermoelectric power unit l in the time interval; h ishsThe storage and heat release power of the heat storage tank in the time period; heat release is positive, heat storage is negative, hload,thThe thermal load of the system for the period of time; NL is the number of all thermoelectric units;
(3) wind power output restraint:
Figure GDA0002624803370000131
in the formula (I), the compound is shown in the specification,
Figure GDA0002624803370000132
is the rated power, v, of the wind turbineCI,vRAnd vCORespectively representing the cut-in wind speed, the rated wind speed and the cut-out wind speed of the fan, vhtIs the wind speed for a certain period of time;
(4) unit restraint:
1) and (3) restraining the upper and lower limits of the unit output:
pi,min≤pi,t≤pi,max (13)
in the formula, pi,min、pi,maxRespectively the minimum output and the maximum output of the unit under the pure condensation working condition;
2) the upper limit and the lower limit of the thermal output of the thermoelectric unit are restricted:
0≤hi,t≤hi,max (14)
in the formula, hi,maxThe maximum limit value of the thermal output of the unit i is determined, and the value is mainly determined by the capacity of the heat exchanger;
3) and (3) total power climbing constraint of the unit:
Pith-Pit(h-1)≤URi(1-yith)+Pi minyith (15)
Pith(h-1)-Pith≤DRi(1-zith)+Pi minzith (16)
in the formula, URiIs a ramp-up limit, yithIs whether the unit is started or not, Pi minIs the minimum power generation amount, DR, of the generator setiIs a ramp down limit, zithIs the amount of whether the unit is in a shutdown state or not;
4) thermal ramp restraint of the thermoelectric unit:
hi,t-hi,t-1≤Δhu,i (17)
hi,t-1-hi,t≤Δhd,i (18)
in the formula,. DELTA.hu,i、Δhd,iRespectively the maximum variation of the thermal power of the extraction type unit in unit time;
(5) and (3) electric storage system constraint:
the charging process is
Figure GDA0002624803370000141
The discharge process is
Figure GDA0002624803370000142
SOC (t) is the residual electric quantity of the energy storage system at the end of the tth period; SOC (t-1) is the residual electric quantity of the energy storage system at the end of the t-1 th time period; delta is the self-discharge rate of the energy storage system; pc、PdRespectively charging and discharging power of the energy storage system; etac、ηdRespectively the charging and discharging efficiency of the electricity storage system;
Figure GDA0002624803370000143
rated capacity of the power storage system;
(6) and (3) heat storage system constraint:
Figure GDA0002624803370000144
in the formula, HHS(t) thermal energy storage capacity for time period t; mu is heat storage and heat dissipation loss rate; qHS_ch(t)、QHS_dis(t) and ηhch、ηhdisRespectively the heat absorption and discharge power and efficiency in the time period t.
Step 3, aiming at the optimal configuration model of the energy storage system, estimating the service life of the storage battery by adopting a rain flow counting method to obtain an investment coefficient;
the method for estimating the service life of the storage battery by adopting the rain flow counting method comprises the following specific steps:
the method is mainly used in engineering, and is particularly widely applied to fatigue life calculation. The strain-time history data recording is rotated by 90 degrees, the time coordinate axis is vertical downwards, the data recording is like a series of roofs, and rainwater flows downwards along the roofs, so the method is called a rainwater flow counting method.
The service life of the energy storage battery is influenced by factors including the discharge depth, the rate performance, the charge-discharge cut-off voltage and the ambient temperature of the battery; the influence of the multiplying power performance of the battery on the service life of the battery is not considered, and the rated value is taken as the maximum power of the energy storage battery; the influence of the charge-discharge cut-off voltage of the battery on the service life of the battery is not considered, and the capacity value range of the energy storage battery is set; the ambient temperature is taken as room temperature, regardless of the effect of temperature on battery life;
1) the SOC-time curve is rotated clockwise, the rain flow starts at the beginning of the recording and consequently at the inner edge of each peak;
2) rain drops vertically down at the peak (i.e., eave), one to the opposite with a maximum more positive than the initial maximum, or one to the opposite with a minimum more negative than the initial minimum;
3) when the rain flow meets the rain flowing from the roof above, the flow is stopped, and a cycle is formed;
4) drawing each cycle according to the starting point and the end point of the raindrop flow, taking out all the cycles one by one, and recording the peak-valley value of the cycles;
5) the horizontal length of each rain stream is taken as the depth of discharge for that cycle.
As shown in FIG. 3, the gray solid line (A-B-C-D-E-F-G) represents the battery SOC variation curve. By using a rain flow counting method, a cycle counting period 1(B-C-B ', the depth of discharge is 0.075), a cycle 2(E-F-E', the depth of discharge is 0.05), a cycle counting half period 3(A-B-B '-D, the depth of discharge is 0.28) and a half period 4(D-E-E' -G, the depth of discharge is 0.2) can be obtained.
Through simplification, the life of the energy storage battery is estimated by using the following formula:
Figure GDA0002624803370000151
T=ceil(1/365·dloss) (23)
wherein d islossThe damage rate of the service life of the electric energy storage system is one day, theta is a period coefficient, the full period is 1, and the half period is 0.5; cyciIs the ith cycleMaximum cycle number corresponding to the cycle of the loop; t is the life cycle, ceil is the ceiling function; the service life of the energy storage system can be estimated by constructing a charging and discharging curve of the energy storage system in one day.
And 4, adopting a proper operation strategy in the simulation operation process by predicting the annual load data and the wind power data, and selecting a traversal algorithm or a particle swarm algorithm according to the optimization range to seek the optimal energy storage system configuration.
The method comprises the following specific steps of selecting a traversal algorithm or a particle swarm algorithm according to an optimization range to seek the optimized energy storage system configuration, wherein the specific steps are as follows: determining an adopted method according to the configuration range of the energy storage system, and if the range is smaller and the configuration types are less, adopting a traversal algorithm; if the range is large, the particle swarm algorithm is adopted to seek the optimal energy storage system configuration to respectively obtain the capacity and the power of the energy storage system and the capacity and the power of the heat storage system; when the optimizing range is smaller, selecting one configuration from the configuration ranges of the energy storage system meeting the self-sufficient probability condition, calculating the cost and judging whether the configuration used in the given range is completed or not, completing the calculation process if the traversal of all the configurations in the given range is completed, and continuing the calculation of the next configuration if the traversal is not completed; the method is comprehensive and is suitable for the optimization problem with a small range;
when the optimizing range is large, a particle swarm optimization algorithm is adopted, wherein the algorithm is an iterative-based optimization tool; the algorithm starts from a group of random solutions and searches an optimal solution through continuous iteration; in each iteration, updating the population by tracking individual extrema and global extrema; the algorithm has the advantages of simplicity in implementation, high speed in convergence, high precision and the like, and is widely applied to engineering practice.
The particle swarm algorithm comprises the following basic steps:
1) randomly initializing the position and speed of each particle in the population;
2) evaluating the fitness of each particle, storing the current position and the adaptive value of each particle in pbest of each particle, and storing the position and the adaptive value of an individual with the optimal adaptive value in all pbest in gbest;
3) by the formula (17) hi,t-hi,t-1≤Δhu,iUpdate the velocity and displacement of the particles:
Figure GDA0002624803370000171
4) for each particle, comparing its fitness value with its experienced best position, and if so, taking it as the current best position;
5) comparing the values of all pbest and gbest at present, and updating the gbest;
6) if the preset stopping condition of the operation precision or the iteration times is met, stopping the search and outputting the result, otherwise, returning to the step 3 and continuing the search.
The above-mentioned embodiments are merely illustrative of the preferred embodiments of the present invention, and do not limit the scope of the present invention, and various modifications and improvements of the technical solution of the present invention by those skilled in the art should fall within the protection scope defined by the claims of the present invention without departing from the spirit of the present invention.

Claims (3)

1. An independent micro energy network energy storage system optimal configuration solving method is characterized in that:
the configuration of the energy storage system influences the investment cost of the energy storage system, the operation cost of the whole microgrid and the pollution control cost; the lower energy storage system can not achieve the expected economical efficiency and stability of the system, the operation cost can not be effectively reduced, and the discharged CO can not be effectively reduced2And the content of harmful gases is higher; the higher energy storage system is configured, so that the investment cost is higher, and the overall maintenance cost is also relatively higher; therefore, the configuration of the optimal energy storage system is selected to be balanced among the investment cost, the operation cost and the pollution control cost, and the sum of the investment cost, the operation cost and the pollution control cost is found at a balance point, namely the configuration of the energy storage system with the minimum total cost is obtained;
the method comprises the following steps:
step 1, establishing self-sufficiency probability requirements of a micro energy network;
step 2, establishing an energy storage system optimal configuration model of the micro energy network;
step 3, aiming at the optimal configuration model of the energy storage system, estimating the service life of the storage battery by adopting a rain flow counting method to obtain an investment coefficient;
step 4, by predicting load data and wind power data of one year, adopting a proper operation strategy in the simulation operation process, and selecting a traversal algorithm or a particle swarm algorithm according to an optimization range to seek the optimal energy storage system configuration;
in step 1, the self-sufficient probability requirement is modeled as:
Figure FDA0002624803360000011
Figure FDA0002624803360000012
wherein, PSSeAnd PSShThe self-sufficient probabilities of the electric load and the heat load in the micro energy network are respectively, and the delta w, the delta d and the delta h respectively meet the wind output prediction error of normal distribution;
in step 2, the specific process of establishing the energy storage system optimization configuration model of the micro energy network includes the following steps:
the comprehensive total cost of investment, operation and pollution control environment and fuel cost needs to be considered in the model, and the model is established by taking the minimum total cost as an objective function:
Min(IC+OC+PC) (3)
Figure FDA0002624803360000021
Figure FDA0002624803360000022
Figure FDA0002624803360000023
wherein, IC is the investment cost of the energy storage system, OC is the operation cost of the micro energy network, PC is the pollution control cost of the micro energy network, alpha is the unit power investment coefficient of the storage battery, beta is the unit capacity investment coefficient of the storage battery, chi is the unit power investment coefficient of the heat storage system, delta is the unit capacity investment coefficient of the heat storage system,
Figure FDA0002624803360000024
is the maximum power of the storage battery,
Figure FDA0002624803360000025
is the maximum capacity of the storage battery,
Figure FDA0002624803360000026
is the maximum power of the heat storage system,
Figure FDA0002624803360000027
is the maximum capacity of the heat storage system, NT is the total number of days, NH is the total number of hours, NG is the total number of conventional thermal power units, NL is the total number of cogeneration units, PithIs the power, F, generated by a conventional unit during a certain period of timeeIs a function of the power and cost, IithIs the state index of whether the distributed power supply works, the work is 1, the non-work is 0, PlthIs the power, F, generated by a cogeneration unit during a certain period of timehIs a function of the power and cost, LlthThe index is the state index of whether the distributed power supply works, the work is 1, and the non-work is 0; SUth,SDthStarting or stopping costs, alpha, of the generator set, respectivelyKRepresenting the cost coefficient of treatment, beta, of different pollutantsKExpressing the emission coefficients of different pollutants, and NK expresses the total amount of the pollutants;
wherein, the unit capacity investment coefficient of the storage battery is as follows:
Figure FDA0002624803360000031
in the formula, CEIs the total electrical energy storage unit capacity investment cost, CmIs the sum of the repair, maintenance and equipment disposal costs per unit volume of equipment; the expenses are evenly distributed in the period life of the energy storage system, so that a unit investment coefficient in a planning period is obtained;
the coal consumption cost is expressed as a quadratic function form of the generated power; the power cost functions of the conventional thermal power generating unit and the cogeneration unit are respectively as follows:
Figure FDA0002624803360000032
Figure FDA0002624803360000033
the constraint conditions of the energy storage system optimization configuration model of the micro energy network are as follows:
electric power balance constraint, heat supply balance constraint, wind power output constraint, unit constraint, power storage system constraint and heat storage system constraint; the unit constraints comprise unit output upper and lower limit constraints, extraction type unit heat output upper and lower limit constraints, extraction type unit net power generation output upper and lower limit constraints, total power climbing constraints of the unit and heat climbing constraints of the extraction type unit; the specific content of the constraint condition is as follows:
(1) electric power balance constraint:
Figure FDA0002624803360000034
where NR is the amount of new energy, PrthIs the power generated by the new energy source, PESSThe charging or discharging power of the energy storage system is provided, the charging is a negative value, and the discharging is a positive value; pload,thThe power required by the load for the period of time;
(2) and (3) heat supply balance constraint:
Figure FDA0002624803360000035
in the formula, hlthThe thermal power of the thermoelectric power unit l in the time interval; h ishsThe storage and heat release power of the heat storage tank in the time period; heat release is positive, heat storage is negative, hload,thThe thermal load of the system for the period of time; NL is the number of all thermoelectric units;
(3) wind power output restraint:
Figure FDA0002624803360000041
in the formula (I), the compound is shown in the specification,
Figure FDA0002624803360000042
is the rated power, v, of the wind turbineCI,vRAnd vCORespectively representing the cut-in wind speed, the rated wind speed and the cut-out wind speed of the fan, vhtIs the wind speed for a certain period of time;
(4) unit restraint:
1) and (3) restraining the upper and lower limits of the unit output:
pi,min≤pi,t≤pi,max (13)
in the formula, pi,min、pi,maxRespectively the minimum output and the maximum output of the unit under the pure condensation working condition;
2) the upper limit and the lower limit of the thermal output of the thermoelectric unit are restricted:
0≤hi,t≤hi,max (14)
in the formula, hi,maxThe maximum limit value of the heat output of the unit is determined, and the value is mainly determined by the capacity of the heat exchanger;
3) and (3) total power climbing constraint of the unit:
Pith-Pit(h-1)≤URi(1-yith)+Pi minyith (15)
Pit(h-1)-Pith≤DRi(1-zith)+Pi minzith (16)
in the formula, URiIs a ramp-up limit, yithIs whether the unit is started or not, Pi minIs the minimum power generation amount, DR, of the generator setiIs a ramp down limit, zithIs the amount of whether the unit is in a shutdown state or not;
4) thermal ramp restraint of the thermoelectric unit:
hi,t-hi,t-1≤Δhu,i (17)
hi,t-1-hi,t≤Δhd,i (18)
in the formula,. DELTA.hu,i、Δhd,iRespectively the maximum variation of the thermal power of the extraction type unit in unit time;
(5) and (3) electric storage system constraint:
the charging process is
Figure FDA0002624803360000051
The discharge process is
Figure FDA0002624803360000052
SOC (t) is the residual electric quantity of the energy storage system at the end of the tth period; SOC (t-1) is the residual electric quantity of the energy storage system at the end of the t-1 th time period; delta is the self-discharge rate of the energy storage system; pc、PdRespectively charging and discharging power of the energy storage system; etac、ηdRespectively the charging and discharging efficiency of the electricity storage system;
Figure FDA0002624803360000053
rated capacity of the power storage system;
(6) and (3) heat storage system constraint:
Figure FDA0002624803360000054
in the formula, HHS(t) thermal energy storage capacity for time period t; mu is heat storage and heat dissipation loss rate; qHS_ch(t)、QHS_dis(t) and ηhch、ηhdisRespectively the heat absorption and discharge power and efficiency in the time period t.
2. The method for solving the optimal configuration of the independent micro energy network energy storage system according to claim 1, wherein the method comprises the following steps: in step 3, estimating the service life of the storage battery by using a rain flow counting method, which comprises the following specific steps:
the service life of the energy storage battery is influenced by factors including the discharge depth, the rate performance, the charge-discharge cut-off voltage and the ambient temperature of the battery; the influence of the multiplying power performance of the battery on the service life of the battery is not considered, and the rated value is taken as the maximum power of the energy storage battery; the influence of the charge-discharge cut-off voltage of the battery on the service life of the battery is not considered, and the capacity value range of the energy storage battery is set; the ambient temperature is taken as room temperature, regardless of the effect of temperature on battery life; estimating the life of the energy storage battery by using the following formula:
Figure FDA0002624803360000061
T=ceil(1/365·dloss) (23)
wherein d islossThe damage rate of the service life of the electric energy storage system is one day, theta is a period coefficient, the full period is 1, and the half period is 0.5; cyciThe maximum cycle number corresponding to the ith cycle period; t is the life cycle, ceil is the ceiling function; the service life of the energy storage system can be estimated by constructing a charging and discharging curve of the energy storage system in one day.
3. The method for solving the optimal configuration of the independent micro energy network energy storage system according to claim 1, wherein the method comprises the following steps: in step 4, the energy storage system configuration which is optimized according to the optimization range selection traversal algorithm or the particle swarm optimization algorithm is searched, and the specific process is as follows:
determining an adopted method according to the configuration range of the energy storage system, and if the range is smaller and the configuration types are less, adopting a traversal algorithm; if the range is large, the particle swarm algorithm is adopted to seek the optimal energy storage system configuration to respectively obtain the capacity and the power of the energy storage system and the capacity and the power of the heat storage system; when the optimizing range is smaller, selecting a configuration calculation cost from the configuration range of the energy storage system meeting the self-sufficient probability condition, judging whether the configuration used in the given range is completed or not, completing the calculation process if the traversal of all the configurations in the given range is completed, and continuing the calculation of the next configuration if the traversal is not completed; the method is comprehensive and is suitable for the optimization problem with a small range;
when the optimizing range is large, a particle swarm optimization algorithm is adopted, wherein the algorithm is an iterative-based optimization tool; the algorithm starts from a group of random solutions and searches an optimal solution through continuous iteration; in each iteration, updating the population by tracking individual extrema and global extrema; the particle swarm algorithm comprises the following basic steps:
1) randomly initializing the position and speed of each particle in the population;
2) evaluating the fitness of each particle, storing the current position and the adaptive value of each particle in pbest of each particle, and storing the position and the adaptive value of an individual with the optimal adaptive value in all pbest in gbest;
3) by the formula hl,t(h-1)-hl,th≤Δhd,lUpdate the velocity and displacement of the particles:
Figure FDA0002624803360000071
4) for each particle, comparing its fitness value with its experienced best position, and if so, taking it as the current best position;
5) comparing the values of all pbest and gbest at present, and updating the gbest;
6) if the preset stopping condition of the operation precision or the iteration times is met, stopping the search and outputting the result, otherwise, returning to the step 3 and continuing the search.
CN201710431745.7A 2017-06-09 2017-06-09 Independent micro-energy-grid energy storage system optimal configuration solving method Active CN107134810B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710431745.7A CN107134810B (en) 2017-06-09 2017-06-09 Independent micro-energy-grid energy storage system optimal configuration solving method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710431745.7A CN107134810B (en) 2017-06-09 2017-06-09 Independent micro-energy-grid energy storage system optimal configuration solving method

Publications (2)

Publication Number Publication Date
CN107134810A CN107134810A (en) 2017-09-05
CN107134810B true CN107134810B (en) 2021-01-08

Family

ID=59733939

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710431745.7A Active CN107134810B (en) 2017-06-09 2017-06-09 Independent micro-energy-grid energy storage system optimal configuration solving method

Country Status (1)

Country Link
CN (1) CN107134810B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107633330A (en) * 2017-09-21 2018-01-26 深圳低碳城综合能源有限公司 System and method are coordinated in the multipotency interconnection that optimization is searched based on tree
CN108270244B (en) * 2018-01-26 2022-08-09 国网辽宁省电力有限公司 Method for scheduling new energy power system in multiple control domain operation modes
CN108649553A (en) * 2018-03-23 2018-10-12 中国电力科学研究院有限公司 A kind of the configuration selection method and device of micro- energy net energy supplying system
CN110323786A (en) * 2018-03-28 2019-10-11 华北电力大学 Dispatching method and device based on micro-capacitance sensor
CN109919478B (en) * 2019-02-28 2023-04-07 天津大学 Comprehensive energy microgrid planning method considering comprehensive energy supply reliability
CN110429642B (en) * 2019-07-19 2022-07-12 昆明理工大学 Micro-grid distributed scheduling method considering renewable energy
CN110365058A (en) * 2019-08-15 2019-10-22 中国人民解放军陆军勤务学院 Isolated micro-capacitance sensor running optimizatin method based on simulated annealing particle swarm algorithm
CN111489020A (en) * 2020-03-31 2020-08-04 杭州鸿晟电力设计咨询有限公司 Independent type comprehensive energy grid electricity-gas energy storage system optimal configuration solving method
CN111523204B (en) * 2020-03-31 2023-09-22 杭州鸿晟电力设计咨询有限公司 Optimal configuration solving method for grid-connected comprehensive energy grid electricity-gas energy storage system
CN112488370B (en) * 2020-11-20 2023-09-22 西安热工研究院有限公司 Energy storage optimization method adopting demand side response planning
CN113762643B (en) * 2021-09-23 2023-07-14 国网湖南省电力有限公司 Energy storage capacity optimization configuration method of regional comprehensive energy system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105574620A (en) * 2016-01-06 2016-05-11 重庆邮电大学 Micro-grid dynamic optimization scheduling method combined with double master control cooperation and MPSO algorithm
CN105787605A (en) * 2016-03-24 2016-07-20 上海电力学院 Micro-grid economic and optimal operation and scheduling method based on improved quantum genetic algorithm
CN106022503A (en) * 2016-03-17 2016-10-12 北京睿新科技有限公司 Micro-grid capacity programming method meeting coupling type electric cold and heat demand

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105574620A (en) * 2016-01-06 2016-05-11 重庆邮电大学 Micro-grid dynamic optimization scheduling method combined with double master control cooperation and MPSO algorithm
CN106022503A (en) * 2016-03-17 2016-10-12 北京睿新科技有限公司 Micro-grid capacity programming method meeting coupling type electric cold and heat demand
CN105787605A (en) * 2016-03-24 2016-07-20 上海电力学院 Micro-grid economic and optimal operation and scheduling method based on improved quantum genetic algorithm

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
考虑储热装置与碳捕集设备的风电消纳低碳经济调度;卢志刚等;《电工技术学报》;20160930;第31卷(第17期);全文 *

Also Published As

Publication number Publication date
CN107134810A (en) 2017-09-05

Similar Documents

Publication Publication Date Title
CN107134810B (en) Independent micro-energy-grid energy storage system optimal configuration solving method
CN109858759B (en) Industrial park comprehensive energy balance scheduling method
CN109767078B (en) Multi-type power supply maintenance arrangement method based on mixed integer programming
CN111144620B (en) Electric hydrogen comprehensive energy system considering seasonal hydrogen storage and robust planning method thereof
CN109919399B (en) Day-ahead economic dispatching method and system for comprehensive energy system
CN109103929B (en) Power distribution network economic optimization scheduling method based on improved dynamic kriging model
CN111737884B (en) Multi-target random planning method for micro-energy network containing multiple clean energy sources
CN110084465B (en) Wind power generation system cost/power supply reliability evaluation method based on energy storage
CN111445107A (en) Multi-objective optimization configuration method for cold-heat-power combined supply type micro-grid
CN111030094A (en) Method for evaluating reliability of energy supply of multi-energy complementary park
CN116402210A (en) Multi-objective optimization method, system, equipment and medium for comprehensive energy system
CN117081143A (en) Method for promoting coordination and optimization operation of park comprehensive energy system for distributed photovoltaic on-site digestion
CN109861292B (en) Method for improving clean energy consumption based on multi-energy storage system
CN113435659B (en) Scene analysis-based two-stage optimized operation method and system for comprehensive energy system
Monterrat et al. Integration Cost of Variable Renewable Resources to Power Systems–A Techno-economic Assessment in European Countries
Dong et al. Sizing of a stand-alone photovoltaic/wind energy system with hydrogen and battery storage based on improved ant colony algorithm
CN116205458A (en) Method for establishing comprehensive energy system optimization scheduling model of hydroelectric hydrogen production by considering uncertainty
CN113962528A (en) Multi-scene optimized mixed renewable energy system commissioning decision method and system
Ma et al. Multi-Point Layout Planning of Multi-Energy Power Supplies Based on Time-series Production Simulation
CN113536581A (en) Energy storage system multi-state reliability modeling method considering operation strategy
Karemore et al. Multi-objective Design procedure for hybrid (wind–photovoltaic) system by GA
Yu et al. Cooperative operation of chemical-free energy storage system with solar photovoltaic for resilient power distribution in buildings—A case study
Venayagamoorthy et al. Grid operations with high penetration of photovoltaic systems
Thaherkhani et al. Modelling Optimal PV System Sizing for Zero Energy Buildings
Tang et al. Research on Capacity Configuration Optimization of Multi-Energy Complementary System Using Deep Reinforce Learning

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