CN113343478B - Independent microgrid capacity optimal configuration method considering uncertainty and demand response - Google Patents
Independent microgrid capacity optimal configuration method considering uncertainty and demand response Download PDFInfo
- Publication number
- CN113343478B CN113343478B CN202110704011.8A CN202110704011A CN113343478B CN 113343478 B CN113343478 B CN 113343478B CN 202110704011 A CN202110704011 A CN 202110704011A CN 113343478 B CN113343478 B CN 113343478B
- Authority
- CN
- China
- Prior art keywords
- cost
- load
- time
- power
- formula
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000004044 response Effects 0.000 title claims abstract description 20
- 238000005457 optimization Methods 0.000 claims abstract description 48
- 230000005611 electricity Effects 0.000 claims abstract description 33
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 230000002068 genetic effect Effects 0.000 claims abstract description 5
- 102100031186 Chromogranin-A Human genes 0.000 claims description 45
- 101000993094 Homo sapiens Chromogranin-A Proteins 0.000 claims description 45
- 238000004146 energy storage Methods 0.000 claims description 29
- 230000005855 radiation Effects 0.000 claims description 22
- 238000010248 power generation Methods 0.000 claims description 16
- 238000005338 heat storage Methods 0.000 claims description 15
- 239000000446 fuel Substances 0.000 claims description 9
- 238000009434 installation Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 7
- 230000001932 seasonal effect Effects 0.000 claims description 7
- 238000007599 discharging Methods 0.000 claims description 6
- 238000013178 mathematical model Methods 0.000 claims description 5
- 238000012423 maintenance Methods 0.000 claims description 4
- 238000012546 transfer Methods 0.000 claims description 4
- 230000014759 maintenance of location Effects 0.000 claims description 3
- 239000002283 diesel fuel Substances 0.000 claims description 2
- 239000012530 fluid Substances 0.000 claims description 2
- 230000008569 process Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 229910052731 fluorine Inorganic materials 0.000 claims 1
- 125000001153 fluoro group Chemical group F* 0.000 claims 1
- 238000010586 diagram Methods 0.000 abstract description 9
- 230000008901 benefit Effects 0.000 abstract description 3
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 7
- 229910052744 lithium Inorganic materials 0.000 description 7
- 230000000694 effects Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- HBBGRARXTFLTSG-UHFFFAOYSA-N Lithium ion Chemical compound [Li+] HBBGRARXTFLTSG-UHFFFAOYSA-N 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910001416 lithium ion Inorganic materials 0.000 description 2
- 102000010029 Homer Scaffolding Proteins Human genes 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
Images
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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- 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
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/30—Wind power
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Geometry (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (AREA)
- Evolutionary Computation (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention relates to an independent microgrid capacity optimal configuration method considering uncertainty and demand response, which is characterized by comprising the following steps of: firstly, giving a basic structure diagram of an independent micro-grid containing heat pump electricity storage, and explaining the working principle of heat pump electricity storage; secondly, establishing an uncertain scene set with adjustable conservative degree based on a robust optimization idea for representing the influence of source load uncertainty on the operation cost of the independent microgrid; then, an independent micro-grid double-layer robust optimization configuration model is constructed, an outer layer optimization model is used for configuring and optimizing the capacity of the wind turbine generator, the photovoltaic array, the diesel generator and the heat pump for power storage by taking the annual value of the total cost and the like as a target, an inner layer optimization model is used for considering the demand response on the basis of the established uncertain scene set and is used for optimizing the operation by taking the annual operation cost of the system as the lowest target; and finally, solving the model by using a genetic algorithm of embedded mixed integer linear programming. The method has the advantages of comprehensive consideration of elements, clear model hierarchy and high calculation efficiency.
Description
Technical Field
The invention relates to the field of planning and design of independent microgrids, in particular to an independent microgrid capacity optimal configuration method considering uncertainty and demand response.
Background
The independent microgrid is a small-sized power system which comprises a distributed power supply, a load and an energy storage and is not connected with a main network, and has complete power generation and distribution functions, so that the independent microgrid gradually becomes a key technology for solving the power supply problem in remote areas. However, wind power and photovoltaic output are obviously intermittent and fluctuating under the influence of weather, and load cannot be accurately predicted under the influence of various uncertain factors, so that the influence of source load uncertainty (wind power output uncertainty and load uncertainty) is considered during the optimal configuration of the microgrid, and an obtained configuration result can better adapt to the actual operation scene in the future. Secondly, when the optimal configuration is carried out, Demand Response (DR) is considered in the microgrid optimal configuration model, the timeshiftable potential of timeshiftable load is exerted, the energy storage capacity required by building an independent microgrid can be reduced, and therefore the investment cost of the microgrid can be further reduced. Moreover, due to high investment and construction costs of the microgrid, particularly, an energy storage device which is one of important components of the independent microgrid is very expensive, a new energy storage technology, namely heat pump electricity storage, is worthy of paying attention to provide a new idea for solving the problem of high energy storage cost, compared with pumped storage and compressed air energy storage, the heat pump electricity storage has higher energy density, so that the unit capacity cost is reduced, and the heat pump electricity storage is basically not limited by geographical conditions.
In summary, the invention provides an independent microgrid capacity optimization configuration method considering uncertainty and demand response. Firstly, an independent micro-grid basic structure containing heat pump electricity storage is established, the operation principle of the heat pump electricity storage is emphatically explained, and a simplified mathematical model of the heat pump electricity storage is provided. And secondly, constructing a microgrid robust operation scene considering source load uncertainty, and laying a foundation for establishing an optimized configuration model in the next step. And then, establishing a microgrid double-layer optimization configuration model, optimizing the quantity or capacity of each device by taking the lowest annual value of the total cost as a target through an outer layer model, considering source load uncertainty and excitation type demand response through an inner layer model, and optimizing each operation variable by taking the lowest annual operation cost of the system as a target. Finally, the reasonability and effectiveness of the method provided by the invention are illustrated by examples.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and establish a scientific and reasonable independent micro-grid capacity optimal configuration method which is strong in applicability and good in effect and can optimize annual values of investment cost and the like on the premise of meeting the reliability requirement, the non-renewable energy power generation ratio requirement and the renewable energy permeability requirement.
The technical scheme adopted for realizing the purpose of the invention is that the independent microgrid capacity optimal configuration method considering uncertainty and demand response is characterized by comprising the following steps of:
1) independent micro-grid basic structure for providing heat pump-containing electricity storage
The independent microgrid comprises a Wind Turbine (WT), a photovoltaic array (PV), a diesel generator (DEG), a load and heat pump electricity storage (PHES) system, wherein the wind turbine and the photovoltaic array are connected to an alternating current bus through different types of converters; the load considers two types of loads, namely rigid load and time-shifting load;
the PHES system comprises a heat storage device, an electric heat conversion device and Working Fluid (WF), wherein the heat storage device comprises a heat storage tank and a cold storage tank, and heat storage media are respectively filled in the heat storage tank and the cold storage tank; the electric-heat conversion device comprises a group of coaxial compressors and expanders, a group of heat exchangers, a group of motors and generators;
the simplified mathematical model of the PHES system is established by referring to a mathematical model of battery energy storage (1):
in the formula: t represents a time period, and the value of t is 1,2, … and 24; s SOC,t The state of charge of the PHES system in the t period; epsilon is the energy self-loss rate; eta in And η out Efficiency of the PHES system in charge and discharge states, respectively;
and
charging and discharging power for the PHES system, respectively; delta t is the time interval of adjacent time periods, and 1h is taken;
and
minimum and maximum values of state of charge, respectively;
and
minimum and maximum powers to charge or discharge, respectively, the PHES system;
and
is a binary variable, represents the working state of the PHES,
a value of 1 indicates that PHES is in a charging state,
when the discharge voltage is 1, the PHES is in a discharge state, and the PHES cannot be simultaneously 1;
2) establishing independent microgrid operation scene considering source load uncertainty
Processing wind and light, namely uncertainty and load of a source in a source load, namely uncertainty of the load in the source load by adopting a robust optimization method, wherein the key of robust optimization is the construction of an uncertainty set; the ranges of wind speed, solar radiation intensity and load can be represented by predicted values plus prediction error; in summary, for the independent microgrid, the indeterminate set is represented by formula (2);
in the formula: the value of i is 1,2 and 3, which respectively represent three uncertain variables of wind speed, solar radiation intensity and load; t represents a time interval, and the values of t are 1,2, … and 24; u. of i,t The value of the ith uncertain variable in the t period after the uncertainty is considered;
the predicted value of the ith uncertain variable in the t period is obtained; x is the number of i,t The variable is a binary variable, when the variable is 1, the variable indicates that the ith uncertain variable obtains an upper limit value or a lower limit value in a t period, and when the variable is 0, the variable indicates that a predicted value is obtained; Δ u i,t Setting the maximum prediction error of the ith uncertain variable in the t time period according to the historical prediction deviation value; gamma-shaped i The conservative coefficient is introduced, represents the number of the ith uncertain variable which is taken as an upper limit value or a lower limit value in one day, and is set as an integer from 0 to 24;
3) establishing independent micro-grid double-layer optimization configuration model
The optimization configuration process not only relates to the planning problem of a long-time scale, namely the determination of a wind turbine generator, a photovoltaic array, a diesel generator and energy storage capacity, but also relates to the operation problem of a short-time scale, namely the determination of abandoned wind power, abandoned light power, output power of the diesel generator and energy storage at different moments, the determination of a time-shiftable load transfer period and the determination of transfer power, and an independent microgrid double-layer optimization configuration model is constructed, wherein the outer layer is responsible for solving a long-time scale variable, and the inner layer is responsible for solving a short-time scale variable;
establishing an outer layer optimization model:
a. establishing an objective function of an outer layer optimization model, wherein the outer layer is an annual value C of the total cost of the operators for constructing the independent micro-grid under the condition of meeting the construction and operation requirements of the independent micro-grid total The lowest is the target; the decision variable is the installation quantity or capacity of each device, the single machine capacity of the wind turbine generator and the photovoltaic module is given, and the variable needing to be optimized is the installation quantity; for the heat storage devices and the electrothermal conversion devices of the PHES system, the number of the installed devices is 1, and the variable to be optimized is the installed capacity, the objective functions are equations (3) to (6):
C total ＝min(C inv +C ope ) (3)
in the formula: c inv The equal annual cost of equipment investment, including the annual average investment cost and annual maintenance cost of each equipment; c ope The system annual running cost is calculated by an inner layer optimization model; k is the type of equipment in the independent micro-grid, and 1 to 4 of k respectively represent a wind turbine generator, a photovoltaic array and a heat storage/cooler and an electric heat conversion device of the PHES system; c' f,k Initial investment cost of original unit of kth equipment; considering that the service life of some equipment is lower than the system operation life and needs to be replaced after the operation time reaches the service life, for the convenience of calculation, the unit initial investment cost of each equipment is converted into the equivalent unit initial investment cost C which does not need to be replaced any more by the formula (5) f,k I.e. C f,k Initial investment cost of an equivalent unit of the kth equipment; c m,k Annual maintenance cost per unit quantity or capacity for class k devices; a is k The installation quantity or capacity of the kth equipment; k 1 The equipment is a set of equipment which does not need to be replaced; k 2 The equipment set which needs to be replaced within the operating life of the system is provided; l k Indicating that class k devices need to be in class l k Replacing every year; g k Replacing a set of years for a class k device; r is an equal-year value operator which is used for converting the current value into the equal-year value without considering the inflation of the currency; r is the discount rate; l is the system operation life;
b. establishing constraint conditions of outer layer optimization model
Due to the limited area of construction within the area for installation of various types of equipment, the configurable number or capacity of distributed power sources and heat pump electricity storage systems should satisfy the constraint of equation (7):
in the formula: a is k The installation quantity or capacity of the kth equipment;
for maximum installation of devices of the kth classThe quantity or capacity;
step of establishing inner layer optimizing model
a. Establishing an objective function of an inner layer optimization model, wherein the inner layer considers Demand Response (DR), and the lowest annual operation cost of a system under an operation scene considering source load uncertainty is taken as a target; the objective function is formula (8):
in the formula: c DO,d The oil consumption cost of the diesel engine set on the d day; c DR,d A compensation cost is given to the user participating in incentive type DR for day d; c Q,d The electricity cost is abandoned for the day d; c T,d Punishment cost for the electricity shortage on day d;
aiming at the characteristic that the wind speed, the solar radiation intensity and the load have seasonality, the seasonal typical days of the wind speed, the solar radiation intensity and the load are used for replacing wind speed, the solar radiation intensity and the load data of each day in the season, and then an uncertain set of uncertain variables of each seasonal section is constructed based on the formula (2); thereby, formula (8) is rewritten as:
in the formula: s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; n is a radical of an alkyl radical s Days represented by the s typical day; c DO,s 、C DR,s 、C Q,s And C T,s Respectively the fuel consumption cost of the diesel engine set on the s-th typical day, the compensation cost given to the users participating in DR, the electricity abandoning cost and the electricity shortage punishment cost;
cost 1: calculation of DEG Fuel consumption cost, the fuel consumption cost C of the diesel engine set on the s-th typical day DO,s Calculated from the formulae (10) to (11);
in the formula: s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; t represents a time period, and the value of t is 1,2, … and 24; delta t is the time interval of adjacent time periods, and 1h is taken; rho is the price of diesel oil, and 6.13 yuan/L is taken; f s,t DEG oil consumption for the t-th typical day; p D,s,t DEG output power for a period t;
and
respectively a lowest allowable output force and a highest allowable output force of the DEG; a and b are respectively the slope and intercept coefficient of the DEG power-fuel curve, and are respectively 0.084 and 0.246;
cost 2: calculation of DR cost without considering load P of t-th typical day when DR L,s,t Mainly by the rigid load P during t L0,s,t And time-shiftable loads P L1,s,t Two parts are shown as formula (12);
P L,s,t ＝P L0,s,t +P L1,s,t (12)
after DR is implemented, the operator may load P on timeshifts L1,s,t The water is utilized; the considered incentive DR is realized by an agreement between an operator and a user, and the operator can adjust the working time interval of the time-shiftable load in the load through the controller when needed and compensate the user according to the adjustment amount;
considering that the working period of adjusting the load can influence the comfort of the user, the user is allowed to shift the load forwards and backwards at most for a period of time; therefore, the time-shiftable load P 'at a certain time period after DR is implemented' L1,s,t Represented by formulas (13) to (14):
in the formula: s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; t represents a time period, and the value of t is 1,2, … and 24; p is L1,s,t Time-shiftable loads for the t-th typical day;
and
is a binary variable when
A time of 1 indicates that there is a time-shiftable load to shift to the current slot,
a time shift load of 1 indicates that the time shift load of the time interval is transferred to other time intervals, and obviously, the two time shift loads cannot be simultaneously 1;
and
loads respectively transferred to the t-1 time interval and the t +1 time interval;
and
load transferred from the time period t to the time period t-1 and the time period t +1 respectively;
the DR cost is mainly a compensation cost required for the user to adjust the time-shiftable load, and is proportional to the adjustment amount of the time-shiftable load, and is represented by equation (15):
in the formula: s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; t represents a time period, and the value of t is 1,2, … and 24; delta t is the time interval of adjacent time periods, and 1h is taken; c DR,s DR cost for the s typical day; alpha is the cost for transferring unit electric quantity to compensate for the user, and the invention takes 0.32 yuan/kWh;
and
load transferred from the time period t to the time period t-1 and the time period t +1 respectively;
cost 3: calculating the electricity abandoning cost, which consists of two parts, namely wind abandoning cost and light abandoning cost, and is expressed as formula (16);
in the formula: s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; t represents a time period, and the value of t is 1,2, … and 24; delta t is the time interval of adjacent time periods, and 1h is taken; beta is the cost of unit electricity abandonment, and 0.6 yuan/kWh is taken; n is WT,s,t And n PV,s,t The number of fans and the number of photovoltaic groups which are disconnected with the independent micro-grid at the t-th typical day are respectively; p WT0,s,t And P PV0,s,t The power generation power of a single fan and a single group of photovoltaic are respectively in the t period of the s typical day;
cost 4: calculating the punishment cost of insufficient electric quantity, wherein the operator is punished when the load power supply is insufficient, the punishment cost is in direct proportion to the electric quantity shortage and is expressed by the following formulas (17) to (20):
P′ L,s,t ＝P L0,s,t +P′ L1,s,t (19)
in the formula, s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; t represents a time period, and the value of t is 1,2, … and 24; delta t is the time interval of adjacent time periods, and 1h is taken; gamma is the unit power shortage punishment cost, and 1.3 yuan/kWh is taken; p is T,s,t Power shortage for the t-th typical day; p' L,s,t The load value after the load participates in demand response for the t-th typical day;
and
charging power and discharging power of PHES at the t period of the s typical day are respectively; p WT,s,t And P PV,s,t Wind power and photovoltaic power at the t-th typical day respectively; p D,s,t DEG output power for a period t; p L0,s,t Is the rigid load of the t period of the s typical day; p' L1,s,t Time-shiftable loads for the t-th typical day after DR are implemented; n is WT And n PV The number of the configured fans and the number of the photovoltaic sets are respectively given by the outer layer; n is WT,s,t And n PV,s,t The number of fans and the number of photovoltaic groups which are disconnected with the independent micro-grid at the t-th typical day are respectively; p WT0,s,t And P PV0,s,t The power generation power of a single fan and a single group of photovoltaic are respectively in the t period of the s typical day;
b. determining constraints of an inner optimization model
Constraint 1: a power balance constraint, equation (21);
in the formula, s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; t represents a time period, and the value of t is 1,2, … and 24; delta t is the time interval of adjacent time periods, and 1h is taken; p' L,s,t The load value after the load participates in demand response for the t-th typical day;
and
charging power and discharging power of the PHES at the tth typical day time period, respectively; p WT,s,t And P PV,s,t Wind power and photovoltaic power at the t-th typical day respectively; p is D,s,t The output power for the t-th typical day period DEG; p T,s,t Power shortage for the t-th typical day;
constraint 2: the DEG power generation capacity is restricted, and for the independent microgrid, the DEG power generation capacity should not exceed 20%, as shown in a formula (22);
in the formula, s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; t represents a time interval, and the values of t are 1,2, … and 24; delta t is the time interval of adjacent time periods, and 1h is taken; n is s Days represented by the s typical day; p D,s,t DEG output power for a period t; p' L,s,t The load value after the load participates in demand response for the t-th typical day;
constraint 3: the power abandonment rate lambda cannot be higher than a certain limit value lambda in order to ensure the utilization level of renewable energy resources 0 Represented by formula (23);
in the formula, lambda is the electricity abandoning rate; lambda [ alpha ] 0 The upper limit of the power abandonment rate is; s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; t represents a time interval, and the values of t are 1,2, … and 24; delta t is the time interval of adjacent time periods, and 1h is taken; n is s Days represented by the s typical day; n is a radical of an alkyl radical WT And n PV The number of the configured fans and the number of the photovoltaic sets are respectively; n is WT,s,t And n PV,s,t The number of fans and the number of photovoltaic groups which are disconnected with the independent micro-grid at the t-th typical day are respectively; p WT0,s,t And P PV0,s,t The power generation power of a single fan and a single group of photovoltaic are respectively in the t period of the s typical day;
constraint 4: reliability constraint, system power shortage rate Zeta can not be higher than a certain limit value Zeta 0 ；
In the formula, zeta is the power shortage rate; zeta 0 Is the upper limit of the power shortage; s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; t represents a time period, and the value of t is 1,2, … and 24; delta t is the time interval of adjacent time periods, and 1h is taken; n is s Days represented by the s typical day; p T,s,t Power shortage for the t-th typical day; p' L,s,t The load value after the load participates in demand response for the t-th typical day;
solving of independent micro-grid double-layer optimization configuration model
Solving the outer layer optimization model by using a genetic algorithm with an elite retention strategy; after the configuration result is transmitted to the inner layer by the outer layer, the inner layer optimization model is a mixed integer nonlinear programming model and is difficult to directly solve due to the fact that the nonlinear expression of the oil consumption cost of the diesel generator in the formula (11) causes the inner layer optimization model to be a mixed integer nonlinear programming model, therefore, the formula (11) needs to be equivalently converted into a linear expression (25) by adopting a Big-M method, and then the high-efficiency solution can be realized by programming in a Matlab environment by adopting a Yalmip tool box and a commercial solver Cplex;
in the formula: s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; t represents a time period, and the value of t is 1,2, … and 24; sigma 1,s,t And σ 2,s,t Are binary variables, and the sum of the two is 1, sigma 1,s,t A time of 1 indicates that the diesel unit is in a shutdown state for the s typical day t period, σ 1,s,t When the time is 0, the diesel engine set is in a power generation state in the s typical day t period, sigma 2,s,t A time of 1 indicates that the diesel engine set is in a power generation state in the s typical day t period, sigma 2,s,t When the value is 0, the diesel engine set is in a shutdown state in the t period of the s typical day; m 1 、M 2 、M 3 And M 4 For sufficiently large numbers, M in actual calculation 1 And M 4 The rated oil consumption per hour of the diesel generator can be taken, M 2 And M 3 The rated power of the diesel generator can be taken; f s,t DEG oil consumption for the s typical day t period; p is D,s,t The output power for the t-th typical day period DEG;
and
respectively a lowest allowable output force and a highest allowable output force of the DEG; a and b are the slope and intercept coefficient, respectively, of the DEG power-fuel curve.
The invention relates to an independent microgrid capacity optimal configuration method considering uncertainty and demand response, which comprises the steps of firstly providing a basic structure diagram of an independent microgrid containing heat pump electricity storage, and explaining the working principle of heat pump electricity storage; secondly, establishing an uncertain scene set with adjustable conservative degree based on a robust optimization idea for representing the influence of source load uncertainty on the operation cost of the independent microgrid; then, an independent micro-grid double-layer robust optimization configuration model is constructed, wherein the outer layer model configures the capacity of the wind turbine generator, the photovoltaic array, the diesel generator and the heat pump for power storage by taking the annual value of the total cost as a target, the inner layer model considers the demand response on the basis of the established uncertain scene set, and the operation optimization is carried out by taking the annual operation cost of the system as the target; finally, the model is solved by a genetic algorithm with embedded mixed integer linear programming. The method has the advantages of being scientific and reasonable, comprehensive in considered elements, clear in model level, high in calculation efficiency, strong in applicability and good in effect.
Drawings
Fig. 1 is a basic structure diagram of a microgrid when a PHES is in a charging state;
fig. 2 is a basic structure diagram of the microgrid when the PHES is in a discharge state;
FIG. 3 is a graph of annual wind speed over time within an independent microgrid;
FIG. 4 is a graph of annual solar radiation intensity over time within an individual piconet;
FIG. 5 is a graph of annual load in an isolated microgrid over time;
FIG. 6 is a graph of wind speed within an individual micro-net at different typical times of day;
FIG. 7 is a graph of solar radiation intensity within individual micro-nets at different typical times of day;
FIG. 8 is a graph of load in individual micro-nets at different typical times of day;
fig. 9 is a flowchart of a solution of an independent microgrid double-layer optimization configuration model;
FIG. 10 is a diagram of the operation of the spring and autumn typical day under scheme 2;
FIG. 11 is a diagram illustrating the operation of the spring and autumn typical day in case 3;
FIG. 12 is a diagram of the operation of the typical spring and autumn day under scenario 4;
FIG. 13 is a diagram of the operation of a typical winter day for scenario 4;
fig. 14 is a diagram of the operation condition of a typical winter day in the case of the scheme 5.
Detailed Description
The invention relates to an independent microgrid capacity optimal configuration method considering uncertainty and demand response, which comprises the following steps:
1) independent micro-grid basic structure for providing heat pump-containing electricity storage
The independent micro-grid mainly comprises a wind turbine generator, a photovoltaic array, a diesel generator, a load and a PHES system, wherein the wind turbine generator and the photovoltaic array are connected to an alternating current bus through converters of different types; the load takes into account both rigid and time-shiftable types of load. In order to conveniently describe the operation principle of the PHES in different working states, independent microgrid basic structures of the PHES in a charging state and a discharging state are constructed, which correspond to fig. 1 (charging state) and fig. 2 (discharging state), respectively.
2) Establishing independent microgrid operation scene considering source load uncertainty
The meteorological data and the load data of a certain area in a certain year are selected for example analysis, wherein the hour-scale data of the wind speed and the load are actually measured data, the solar radiation intensity data are generated by combining average daily total solar radiation of the area in a month with Homer software, the annual wind speed change curve with time is shown in figure 3, the annual solar radiation intensity change curve with time is shown in figure 4, the annual load change curve with time is shown in figure 5, the annual average wind speed is 7.80m/s, the daily average solar radiation is 4.13kWh/m2, the average load is 1.02MW, and the maximum load is 1.77 MW.
Because 365 days of the year are simultaneously optimized, the optimized variables are more, the calculated amount is larger, and therefore, aiming at the characteristic that the wind speed, the solar radiation intensity and the load are seasonal, the seasonal typical days of the wind speed, the solar radiation intensity and the load are used for replacing the wind speed, the solar radiation intensity and the load data of each day in the season, and then the uncertain set of the uncertain variables of each seasonal typical day is constructed based on the formula (2). When an uncertain set of typical day uncertain variables of each season is constructed, the conservative coefficient of wind speed, solar radiation intensity and load is calculated by taking 6 as an example; the maximum prediction errors of the wind speed, the solar radiation intensity and the load are calculated by taking 10% as an example, and the wind speed typical day scene of different season segments with uncertainty is obtained and is respectively shown in fig. 6, the solar radiation intensity typical day scene of different season segments with uncertainty is obtained and is shown in fig. 7, and the load typical day scene of different season segments with uncertainty is obtained and is shown in fig. 8.
3) Establishing independent micro-grid double-layer optimization configuration model
An independent microgrid double-layer optimization configuration model is constructed by using formulas (3) to (24), and all parameters in the model are set as: economic parameters of the distributed power supply are shown in table 1; the economic and technical parameters of lithium battery energy storage and PHES are shown in Table 2; the system operation year is 20 years; the current sticking rate is 5%; the upper limit of the electricity abandonment rate is 10 percent; the upper limit of the power shortage rate is 0.3 percent; each time interval can be time-shifted by 10% of the load of the time interval; the initial charge state of the energy storage is 0.2, and the maximum charge-discharge power is 20% of the rated capacity.
TABLE 1 economic and technical parameters of distributed power supplies
4) Solving independent micro-grid double-layer optimization configuration model
The solving flow chart is shown in fig. 9. Solving the outer layer optimization model by using a genetic algorithm with an elite retention strategy; after the configuration result is transmitted to the inner layer by the outer layer, the inner layer optimization model is a mixed integer nonlinear programming model and is difficult to directly solve due to the fact that the nonlinear expression of the oil consumption cost of the diesel generator in the formula (11) causes the inner layer optimization model to be a mixed integer nonlinear programming model, therefore, the formula (11) needs to be equivalently converted into a linear form by adopting a Big-M method, as shown in the formula (25), and then the high-efficiency solution can be realized by programming in a Matlab environment by adopting a Yalmip tool box and a commercial solver Cplex.
In order to fully analyze and consider the influence of source charge uncertainty and DR on a configuration result and use the heat pump electricity storage system as an economic advantage of a system energy storage scheme compared with a battery energy storage scheme, the invention sets the following 5 schemes for comparative analysis.
Scheme 1: selecting a lithium ion battery with the highest global application ratio in the current battery energy storage as a system energy storage scheme, and considering source charge uncertainty and DR, wherein the parameters of the lithium ion battery are shown in a table 2;
scheme 2: selecting PHES as a system energy storage scheme, and not considering source load uncertainty and DR;
scheme 3: under the configuration result of the scheme 2, calculating each economic cost of the microgrid when source load uncertainty is considered and the minimum system annual operating cost is taken as a target;
scheme 4: selecting PHES as a system energy storage scheme, and calculating an optimal configuration result when source load uncertainty is considered;
scheme 5: PHES is selected as a system energy storage scheme, and DR is further considered on the basis of considering source load uncertainty, namely the model of the invention.
The optimal configuration results under each scheme are shown in table 3, and the configuration results are analyzed in different ways.
Firstly, analyzing configuration results under different energy storage schemes
The configuration results under the lithium battery energy storage scheme and the PHES energy storage scheme respectively correspond to the scheme 1 and the scheme 2 in the table 3, and comparison shows that: although the energy storage capacity required to be configured when the PHES is adopted is higher than 417kWh when the lithium battery is adopted, the total cost of the PHES energy storage scheme is still reduced by 192.3 ten thousand yuan compared with the total cost of the lithium battery energy storage scheme. This is because: although the cycle efficiency of the PHES is 20% lower than that of the lithium battery, the PHES uses gravel with lower cost as a heat storage medium, so that the unit capacity cost of the PHES is far lower than that of the lithium battery, and the service life of the PHES is longer than that of the lithium battery, so that the replacement cost of equipment is reduced.
TABLE 3 comparison of the results of the optimized configurations for different scenarios
Analysis of influence of uncertainty of source load on configuration result
Comparing scheme 2 and scheme 3, it can be seen that: the same configuration results, the annual operating cost of the system of the scheme 3 is increased by 73.8 ten thousand yuan, and the DEG power generation capacity accounts for and rises by 3.56%. This is because: when uncertainty is not considered, the obtained configuration scheme cannot well cope with the adverse effect (the wind-light output is reduced, the load is increased) brought by source charge uncertainty to the system operation, the use of the diesel generator must be increased in order to meet the requirements of the system power shortage rate and the power abandon rate, and as a result, the system operation cost is greatly increased. Taking the typical spring and autumn days as an example, the operation conditions of the two schemes on the typical days are respectively shown in fig. 10 and fig. 11, and it can be found that: on this typical day, when uncertainty is not taken into account, it is sufficient to satisfy the operating conditions without starting the diesel generator (corresponding to fig. 10), whereas after the uncertainty is added, it is necessary to start the diesel generator only by the fact that the energy storage has failed to satisfy the requirements, thus increasing the cost.
In case of the scheme 4, source load uncertainty is taken into consideration during planning, and compared with the scheme 3, the annual value of investment cost and the like of the scheme 4 is increased by 8.9 ten thousand yuan, the annual running cost of the system is reduced by 27.8 ten thousand yuan, the total cost is reduced by 18.9 ten thousand yuan, and the proportion of the DEG generated energy, the power abandonment rate and the power shortage rate are all reduced. The operation conditions in spring and autumn at the typical day under the scheme 4 are shown in fig. 12, and the diesel generator can be not started and the operation requirements can be met.
In conclusion, although the configuration result considering the source load uncertainty increases the investment cost, the configuration result can better cope with the uncertainty operation scene which may appear in the future, and the economical efficiency is better than the scheme without considering the uncertainty in the long term.
Analysis of influence of DR on configuration results
It is thus demonstrated that if the system operator considers DR after the microgrid has been built, a more rational deployment scenario can be obtained by considering the impact of DR during planning.
While the present invention has been described in detail and with reference to specific embodiments thereof, it will be apparent to one skilled in the art that various changes and modifications can be made therein without departing from the spirit and scope thereof as defined in the appended claims.
Claims (1)
1. An independent microgrid capacity optimization configuration method considering uncertainty and demand response is characterized by comprising the following steps:
1) independent micro-grid basic structure with heat pump for electricity storage
The independent microgrid comprises a Wind Turbine (WT), a photovoltaic array (PV), a diesel generator (DEG), a load and heat pump electricity storage (PHES) system, wherein the wind turbine and the photovoltaic array are connected to an alternating current bus through different types of converters; the load considers two types of loads, namely rigid load and time-shifting load;
the PHES system comprises a heat storage device, an electric heat conversion device and Working Fluid (WF), wherein the heat storage device comprises a heat storage tank and a cold storage tank, and heat storage media are respectively filled in the heat storage tank and the cold storage tank; the electric-heat conversion device comprises a group of coaxial compressors and expanders, a group of heat exchangers, a group of motors and generators;
the simplified mathematical model of the PHES system is established by referring to a mathematical model of battery energy storage (1):
in the formula: t represents a time period, and the value of t is 1,2, … and 24; s SOC,t The state of charge of the PHES system in the t period; epsilon is the energy self-loss rate; eta in And η out Efficiency of the PHES system in charge and discharge states, respectively;
and
charging and discharging power for the PHES system, respectively; delta t is the time interval of adjacent time periods, and 1h is taken;
and
minimum and maximum values of state of charge, respectively;
and
minimum and maximum powers to charge or discharge, respectively, the PHES system;
and
is a binary variable, represents the working state of the PHES,
a value of 1 indicates that PHES is in a charging state,
when the voltage is 1, the PHES is in a discharge state, and the PHES cannot be simultaneously 1;
2) establishing independent microgrid operation scene considering source load uncertainty
Processing wind and light, namely uncertainty and load of a source in a source load, namely uncertainty of the load in the source load by adopting a robust optimization method, wherein the key of robust optimization is the construction of an uncertainty set; the range of wind speed, solar radiation intensity and load can be represented by a predicted value plus a prediction error; in summary, for the independent piconets, the uncertainty set is represented by equation (2);
in the formula: the value of i is 1,2 and 3, which respectively represent three uncertain variables of wind speed, solar radiation intensity and load; u. of i,t The value of the ith uncertain variable in the t period after the uncertainty is considered;
the predicted value of the ith uncertain variable in the t period is obtained; x is a radical of a fluorine atom i,t The variable is a binary variable, when the variable is 1, the variable indicates that the ith uncertain variable obtains an upper limit value or a lower limit value in a t period, and when the variable is 0, the variable indicates that a predicted value is obtained; Δ u i,t Setting the maximum prediction error of the ith uncertain variable in the t time period according to the historical prediction deviation value; gamma-shaped i The conservative coefficient is introduced, represents the number of the ith uncertain variable which is taken as an upper limit value or a lower limit value in one day, and is set as an integer from 0 to 24;
3) establishing independent micro-grid double-layer optimization configuration model
The optimization configuration process not only relates to the planning problem of a long time scale, namely the determination of a wind turbine generator, a photovoltaic array, a diesel generator and energy storage capacity, but also relates to the operation problem of a short time scale, namely the determination of abandoned wind power, abandoned light power, output power of the diesel generator and energy storage at different moments, the determination of a time-shiftable load transfer period and the determination of transfer power, and an independent micro-grid double-layer optimization configuration model is constructed, wherein the outer layer is responsible for solving a long time scale variable, and the inner layer is responsible for solving a short time scale variable;
establishing an outer layer optimization model:
a. establishing an objective function of an outer layer optimization model, wherein the outer layer is an annual value C of the total cost of the operators for constructing the independent micro-grid under the condition of meeting the construction and operation requirements of the independent micro-grid total The lowest is the target; the decision variable is the installation quantity or capacity of each device, the single machine capacity of the wind turbine generator and the photovoltaic module is given, and the variable needing to be optimized is the installation quantity; for the heat storage devices and the electrothermal conversion devices of the PHES system, the number of the installed devices is 1, and the variable to be optimized is the installed capacity, the objective functions are equations (3) to (6):
C total ＝min(C inv +C ope ) (3)
in the formula: c inv The equal annual cost of equipment investment, including the annual average investment cost and annual maintenance cost of each equipment; c ope The annual running cost of the system is calculated by an inner layer optimization model; k is the type of equipment in the independent micro-grid, and 1 to 4 of k respectively represent a wind turbine generator, a photovoltaic array and a heat storage/cooler and an electric heat conversion device of the PHES system; c' f,k Is class kInitial investment cost of original unit of equipment; considering that the service life of some equipment is lower than the system operation life and needs to be replaced after the operation time reaches the service life, for the convenience of calculation, the unit initial investment cost of each equipment is converted into the equivalent unit initial investment cost C which does not need to be replaced any more by the formula (5) f,k I.e. C f,k Initial investment cost of an equivalent unit of the kth equipment; c m,k Annual maintenance cost per unit quantity or capacity for class k devices; a is k The installation quantity or capacity of the kth equipment; k 1 The equipment is a set of equipment which does not need to be replaced; k 2 The equipment set which needs to be replaced within the operating life of the system is provided; l k Indicating that class k devices need to be in class l k Replacing every year; g k Replacing a set of years for a class k device; r is an equal-year value operator which is used for converting the current value into the equal-year value without considering the inflation of the currency; r is the discount rate; l is the system operation life;
b. establishing constraint conditions of outer layer optimization model
Due to the limited area of construction for installing various types of equipment in a region, the configurable number or capacity of distributed power sources and heat pump electricity storage systems should satisfy the constraint of equation (7):
in the formula: a is k The installed number or capacity of the kth class of equipment;
the maximum installed number or capacity of the kth class of equipment;
step of establishing inner layer optimization model
a. Establishing an objective function of an inner-layer optimization model, wherein the inner layer considers Demand Response (DR) and aims to minimize the annual operation cost of a system under an operation scene considering source load uncertainty; the objective function is formula (8):
in the formula: c DO,d The oil consumption cost of the diesel engine set on the d day; c DR,d A compensation cost is given to the user participating in incentive type DR for day d; c Q,d The electricity cost is abandoned for the day d; c T,d Punishment cost for the electricity shortage on day d;
aiming at the characteristic that the wind speed, the solar radiation intensity and the load have seasonality, the seasonal typical days of the wind speed, the solar radiation intensity and the load are used for replacing wind speed, the solar radiation intensity and the load data of each day in the season, and then an uncertain set of uncertain variables of each seasonal section is constructed based on the formula (2); thus, equation (8) is rewritten as:
in the formula: s represents typical days of different season segments, and s takes 1,2 and 3 to represent typical days of spring and autumn, summer and winter respectively; n is a radical of an alkyl radical s Days represented by the s typical day; c DO,s 、C DR,s 、C Q,s And C T,s Respectively the fuel consumption cost of the diesel engine set on the s-th typical day, the compensation cost given to the users participating in DR, the electricity abandoning cost and the electricity shortage punishment cost;
cost 1: calculation of DEG Fuel consumption cost, the fuel consumption cost C of the diesel engine set on the s-th typical day DO,s Calculated from the formulae (10) to (11);
in the formula: rho is the price of diesel oil, and 6.13 yuan/L is taken; f s,t DEG oil consumption for the t-th typical day; p D,s,t The output power for the t-th typical day period DEG;
and
respectively a lowest allowable output force and a highest allowable output force of the DEG; a and b are respectively the slope and intercept coefficient of the DEG power-fuel curve, and are respectively 0.084 and 0.246;
cost 2: calculation of DR cost without considering load P of t-th typical day when DR L,s,t Mainly by the rigid load P during t L0,s,t And time-shiftable loads P L1,s,t Two parts are shown as formula (12);
P L,s,t ＝P L0,s,t +P L1,s,t (12)
after DR is implemented, the operator may load P on timeshifts L1,s,t The water is utilized; the considered incentive DR is realized by an agreement between an operator and a user, and the operator can adjust the working time interval of the time-shiftable load in the load through the controller when needed and compensate the user according to the adjustment amount;
considering that the working period of adjusting the load can influence the comfort of the user, the user is allowed to shift the load forwards and backwards at most for a period of time; therefore, the time-shiftable load P 'at a certain time period after DR is implemented' L1,s,t Represented by formulas (13) to (14):
in the formula: p L1,s,t Time-shiftable loads for the t-th typical day;
and
is a binary variable when
A time of 1 indicates that there is a time-shiftable load to shift to the current slot,
a time shift load of 1 indicates that the time shift load of the time interval is transferred to other time intervals, and obviously, the two time shift loads cannot be simultaneously 1;
and
loads respectively transferred to the t-1 time interval and the t +1 time interval;
and
load transferred from the time period t to the time period t-1 and the time period t +1 respectively;
the DR cost is mainly a compensation cost required for the user to adjust the time-shiftable load, and is proportional to the adjustment amount of the time-shiftable load, and is represented by equation (15):
in the formula: c DR,s DR cost for the s typical day; alpha is the cost for transferring unit electric quantity compensation to the user, and 0.32 yuan/kWh is taken;
cost 3: calculating the electricity abandoning cost, which consists of two parts, namely wind abandoning cost and light abandoning cost, and is expressed as formula (16);
in the formula: beta is the cost of unit electricity abandonment, and 0.6 yuan/kWh is taken; n is WT,s,t And n PV,s,t The number of fans and the number of photovoltaic groups which are disconnected with the independent micro-grid at the t-th typical day are respectively; p WT0,s,t And P PV0,s,t The power generation power of a single fan and a single group of photovoltaic are respectively in the t period of the s typical day;
cost 4: calculating the penalty cost for insufficient power, wherein the operator is penalized when the load power supply is insufficient, the penalty cost is in direct proportion to the power shortage, and the penalty cost is expressed in the following formulas (17) to (20):
P′ L,s,t ＝P L0,s,t +P′ L1,s,t (19)
in the formula, gamma is the penalty cost of unit power shortage, and 1.3 yuan/kWh is taken; p T,s,t Power shortage for the t-th typical day; p' L,s,t The load value after the load participates in demand response for the t-th typical day;
and
respectively, the charging work of PHES during the t period of the s typical dayRate and discharge power; p WT,s,t And P PV,s,t Wind power and photovoltaic power at the t-th typical day respectively; p L0,s,t Is the rigid load of the t period of the s typical day; p' L1,s,t Time-shiftable loads for the t-th typical day after DR are implemented; n is WT And n PV The number of the configured fans and the number of the photovoltaic sets are respectively given by the outer layer; n is WT,s,t And n PV,s,t The number of fans and the number of photovoltaic groups which are disconnected with the independent micro-grid at the t-th typical day are respectively; p is WT0,s,t And P PV0,s,t The power generation power of a single fan and a single group of photovoltaic are respectively in the t period of the s typical day;
b. determining constraints of an inner optimization model
Constraint 1: a power balance constraint, equation (21);
constraint 2: the DEG power generation capacity is restricted, and for the independent microgrid, the DEG power generation capacity should not exceed 20%, as shown in a formula (22);
constraint 3: the power abandonment rate lambda cannot be higher than a certain limit value lambda in order to ensure the utilization level of renewable energy resources 0 Represented by formula (23);
in the formula, lambda is the electricity abandoning rate; lambda [ alpha ] 0 The upper limit of the power abandonment rate is;
constraint 4: reliability constraint, system power shortage rate Zeta can not be higher than a certain limit value Zeta 0 ；
In the formula, zeta is the power shortage rate; zeta 0 Is the upper limit of the power shortage;
solving of independent micro-grid double-layer optimization configuration model
Solving the outer layer optimization model by using a genetic algorithm with an elite retention strategy; after the configuration result is transmitted to the inner layer by the outer layer, the inner layer optimization model is a mixed integer nonlinear programming model and is difficult to solve directly due to the fact that the nonlinear expression of the oil consumption cost of the diesel generator in the formula (11) causes that the inner layer optimization model is a mixed integer nonlinear programming model, so that the formula (11) needs to be equivalently converted into a linear expression (25) by adopting a Big-M method, and then a Yalmip tool box and a commercial solver Cplex are adopted to realize efficient solution in a programming mode under a Matlab environment;
in the formula: sigma 1,s,t And σ 2,s,t Are binary variables, and the sum of the two is 1, sigma 1,s,t A time of 1 indicates that the diesel unit is in a shutdown state at t time period of the s typical day, sigma 1,s,t When the time is 0, the diesel engine set is in a power generation state in the s typical day t period, sigma 2,s,t When the time is 1, the diesel engine set is in a power generation state at t time interval of the s typical day, sigma 2,s,t When the value is 0, the diesel engine set is in a shutdown state in the t period of the s typical day; m 1 、M 2 、M 3 And M 4 For sufficiently large numbers, M in actual calculation 1 And M 4 The rated oil consumption per hour of the diesel generator can be taken, M 2 And M 3 The rated power of the diesel generator can be taken.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110704011.8A CN113343478B (en) | 2021-06-24 | 2021-06-24 | Independent microgrid capacity optimal configuration method considering uncertainty and demand response |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110704011.8A CN113343478B (en) | 2021-06-24 | 2021-06-24 | Independent microgrid capacity optimal configuration method considering uncertainty and demand response |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113343478A CN113343478A (en) | 2021-09-03 |
| CN113343478B true CN113343478B (en) | 2022-09-20 |
Family
ID=77478249
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110704011.8A Active CN113343478B (en) | 2021-06-24 | 2021-06-24 | Independent microgrid capacity optimal configuration method considering uncertainty and demand response |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113343478B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114243694B (en) * | 2021-12-15 | 2023-09-15 | 东北电力大学 | Grid-connected micro-grid optimal configuration method considering ladder carbon transaction and demand response |
| CN114256884A (en) * | 2021-12-21 | 2022-03-29 | 国网山东省电力公司烟台供电公司 | Photovoltaic uncertainty-considered diesel engine energy storage micro-grid vehicle optimal configuration method and system |
| CN115374999A (en) * | 2022-07-06 | 2022-11-22 | 国网浙江省电力有限公司丽水供电公司 | Hydropower hydrogen production optimal configuration method suitable for starting and stopping characteristics of hydrogen production equipment |
| CN115566701B (en) * | 2022-08-26 | 2024-05-10 | 中国长江三峡集团有限公司 | Method and device for determining constant-speed pumping and accumulating capacity and variable-speed pumping and accumulating capacity |
| CN115660187B (en) * | 2022-11-02 | 2024-04-30 | 国家电网有限公司 | Low-carbon town ground source heat pump capacity optimization configuration method |
| CN116722595A (en) * | 2023-04-03 | 2023-09-08 | 中国长江电力股份有限公司 | Demand response and source network collaborative optimization configuration model considering multiple uncertainty factors and solving method thereof |
| CN116914860B (en) * | 2023-09-12 | 2024-01-05 | 中国电建集团西北勘测设计研究院有限公司 | Power supply capacity configuration optimization method and system for multi-energy complementary power generation system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103151798A (en) * | 2013-03-27 | 2013-06-12 | 浙江省电力公司电力科学研究院 | Optimizing method of independent microgrid system |
| CN106487036A (en) * | 2016-11-04 | 2017-03-08 | 浙江工业大学 | A kind of independent photovoltaic dragging system capacity collocation method based on multi-objective optimization algorithm |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107017617B (en) * | 2017-05-26 | 2019-07-26 | 太原理工大学 | The adaptive droop control method of mixed energy storage system in improved direct-current grid |
| CN108009693B (en) * | 2018-01-03 | 2021-09-07 | 上海电力学院 | Grid-connected micro-grid double-layer optimization method based on two-stage demand response |
| CN109713673B (en) * | 2019-02-25 | 2022-12-09 | 上海电力学院 | Method for configuring and optimizing operation of grid-connected micro-grid system in electricity selling environment |
| CN109842158A (en) * | 2019-03-28 | 2019-06-04 | 广东工业大学 | A kind of micro-capacitance sensor Optimal Configuration Method |
| CN110206598B (en) * | 2019-06-04 | 2022-04-01 | 中国科学院工程热物理研究所 | Heat pump energy storage power generation system based on indirect cold storage and heat storage |
| CN111340299B (en) * | 2020-02-29 | 2024-03-22 | 上海电力大学 | Multi-objective optimal scheduling method for micro-grid |
-
2021
- 2021-06-24 CN CN202110704011.8A patent/CN113343478B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103151798A (en) * | 2013-03-27 | 2013-06-12 | 浙江省电力公司电力科学研究院 | Optimizing method of independent microgrid system |
| CN106487036A (en) * | 2016-11-04 | 2017-03-08 | 浙江工业大学 | A kind of independent photovoltaic dragging system capacity collocation method based on multi-objective optimization algorithm |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113343478A (en) | 2021-09-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113343478B (en) | Independent microgrid capacity optimal configuration method considering uncertainty and demand response | |
| Zhong et al. | Towards a 100% renewable energy electricity generation system in Sweden | |
| CN109767105B (en) | Wind, light, water, fire and storage combined system-based multi-energy complementary coordinated power generation scheduling method | |
| Muselli et al. | Design of hybrid-photovoltaic power generator, with optimization of energy management | |
| Nejabatkhah et al. | Optimal design and operation of a remote hybrid microgrid | |
| Blake et al. | Optimization of distributed energy resources in an industrial microgrid | |
| CN111541244A (en) | Power grid side energy storage device capacity calculation method considering power consumption cost of energy storage device | |
| CN110932261A (en) | Multi-energy system combined installation planning method based on global benefit maximization | |
| CN113541195B (en) | Method for consuming high-proportion renewable energy in future power system | |
| CN104616071A (en) | Wind-solar storage complementary generation system configuration optimization method | |
| CN116914860B (en) | Power supply capacity configuration optimization method and system for multi-energy complementary power generation system | |
| Brandeis et al. | Analysis of electrical energy storage technologies for future electric grids | |
| CN110768303A (en) | Optimization configuration method for equipment capacity of island-type energy system | |
| Rehman et al. | Wind farm-battery energy storage assessment in grid-connected microgrids | |
| Hammad et al. | Numerical Analysis Of PV Water Pumping System In Different Geographical Locations In Egypt. | |
| CN110571868B (en) | Optimal configuration method of micro-grid | |
| Sameti | Electrical energy efficient building through distributed generation | |
| CN111786419A (en) | Regional power grid energy configuration system | |
| Yang et al. | Optimal economic dispatch of microgrid with battery energy storage considering uncertainty of renewable generation | |
| CN117791744B (en) | Power supply capacity configuration method of multi-energy complementary power generation system based on hydroelectric energy storage factory | |
| Maharjan et al. | Distributed Energy Systems Design Optimization Based on Life Cycle Environmental and Economic Impacts | |
| Truong et al. | Maximizing solar home battery systems' contribution to the energy transition of the power system | |
| Chen et al. | Mathematical programming model for energy system design | |
| Oksenych et al. | Methodology for Designing the Capacity of Solar Power Plants for an Offline Home Network | |
| Li et al. | Day-head Peak-shaving Model for Coordinated Wind-photovoltaic-pumped-storage-hydropower Generation Systems |
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 |