CN114744658A - Battery energy storage system size division method and system based on micro-grid - Google Patents
Battery energy storage system size division method and system based on micro-grid Download PDFInfo
- Publication number
- CN114744658A CN114744658A CN202210409307.1A CN202210409307A CN114744658A CN 114744658 A CN114744658 A CN 114744658A CN 202210409307 A CN202210409307 A CN 202210409307A CN 114744658 A CN114744658 A CN 114744658A
- Authority
- CN
- China
- Prior art keywords
- energy storage
- battery energy
- storage system
- wolf
- micro
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004146 energy storage Methods 0.000 title claims abstract description 213
- 238000000034 method Methods 0.000 title claims abstract description 25
- 241000282461 Canis lupus Species 0.000 claims abstract description 183
- 239000011159 matrix material Substances 0.000 claims abstract description 33
- 238000005457 optimization Methods 0.000 claims abstract description 17
- 241000282421 Canidae Species 0.000 claims description 32
- 230000000694 effects Effects 0.000 claims description 24
- 239000013598 vector Substances 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 230000007123 defense Effects 0.000 claims description 3
- 238000004513 sizing Methods 0.000 claims description 2
- 238000010248 power generation Methods 0.000 abstract description 8
- 230000002159 abnormal effect Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 241000282465 Canis Species 0.000 description 1
- 244000062662 apex predator Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000000474 nursing effect Effects 0.000 description 1
- 230000036578 sleeping time Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/28—Arrangements for balancing of the load in a network by storage of energy
- H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/10—Flexible AC transmission systems [FACTS]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
Abstract
The invention discloses a battery energy storage system size division method and system based on a microgrid, and relates to the field of the microgrid, wherein an encoding matrix is obtained according to parameter data of each battery energy storage system in the microgrid and the microgrid, and input parameters of a gray wolf optimization algorithm and the position of a wolf group are set according to the encoding matrix; the wolf pack is divided through a fitness objective function, the position of the last alpha wolf is output as the target size of the battery energy storage system after the iteration of the grey wolf optimization algorithm is entered, the grey wolf optimization algorithm is embedded into the battery energy storage system size dividing method, the battery energy storage system with the target size is obtained, and the goal of reasonably controlling the cost of the battery energy storage system is achieved when the micro-grid fails and the frequency is unstable due to imbalance between power generation and load requirements, enough frequency support is provided for the micro-grid to keep the safety of the power system, so that the power system can stably run under normal and abnormal conditions.
Description
Technical Field
The invention relates to the field of micro-grids, in particular to a battery energy storage system size division method and system based on a micro-grid.
Background
As a single controllable network, microgrid cluster loads and micro-sources (SGs and NS-RES) can provide reliable power supply to remote communities (including islands) and industrial complexes. Extensive substitution of Conventional Power Plants (CPPs) by NS-RES has a significant impact on the frequency response of the power system, which can significantly alter the frequency behavior of the system. In other words, power systems with limited inertia have a significant impact on the frequency behavior of the system, for example, loss of a large generator set may result in unacceptable frequency levels, possibly leading to grid collapse, and therefore stable control of the system frequency is particularly important when the penetration level of NS-RES is high.
In addition, in the interconnection mode of the microgrid, it controls and maintains the voltage and frequency of the microgrid within the allowable ranges of frequency and voltage deviations, mainly through the utility grid. However, in the case of independent micro-grids, the NS-RES penetration of the grid is high, and there is an imbalance between the power generation and load demands, so that the frequency deviation of the independent micro-grids is fast.
Currently, with the development of energy materials, Battery Energy Storage Systems (BESS) are widely used in power grids, which provide frequency support by injecting transient power into a microgrid, however, an oversize of a battery energy storage system increases costs, and an undersize results in insufficient battery capacity and thus poor performance. Therefore, the invention provides a battery energy storage system with a proper size, so that when the frequency is unstable due to the failure of the micro-grid and the imbalance between the power generation and the load demand, enough frequency support is provided for the micro-grid to stabilize the frequency response of the micro-grid, the safety of the power system is kept, the power system can stably run under normal and abnormal conditions, and the cost of the battery energy storage system is controlled.
Disclosure of Invention
In order to control the cost of a battery energy storage system while providing enough frequency support for a power grid to stabilize the frequency response of the microgrid when the frequency is unstable due to the fact that the microgrid has a fault and the requirements for power generation and load are unbalanced, the invention provides a battery energy storage system size division method based on the microgrid, which obtains the target size of the battery energy storage system through a wolf optimization algorithm so as to stabilize the frequency response of the microgrid when the microgrid has a fault; the size division method comprises the following steps:
s1: acquiring a coding matrix according to parameter data of each battery energy storage system in the microgrid and the microgrid, wherein elements in the coding matrixThe optional size of a battery energy storage system in the micro-grid is represented, X in elements represents the battery energy storage system, n in elements represents the number of the battery energy storage systems, and p represents the position of the battery energy storage system X;
s2: setting the number of wolves in a wolve group, namely the number of the battery energy storage systems X, the maximum iteration times and the value range of the size parameter variable of the battery energy storage system, namely the range formed by the minimum value of the position of the wolve group and the maximum value of the position of the wolve group through the coding matrix, and initializing the position of the wolve group, namely the position of the battery energy storage system X through the value range;
s3: judging whether each battery energy storage system meets a preset constraint condition or not; if yes, the next step is carried out, if not, the step S2 is returned to reinitialize the wolf pack position;
s4: acquiring a fitness objective function value corresponding to the position of each wolf through the fitness objective function, namely acquiring a frequency deviation effect objective function value of the micro-grid after the micro-grid is subjected to frequency modulation through a battery energy storage system;
s5: dividing a wolf group into alpha, beta, delta and omega according to a fitness objective function value, wherein the alpha wolf is a wolf with the optimal frequency modulation effect of the battery energy storage system, the beta wolf is a wolf with the suboptimal frequency modulation effect of the battery energy storage system, the delta wolf is a wolf with the third best frequency modulation effect of the battery energy storage system, and the omega wolf is the rest wolf;
s6: initializing iteration times and starting counting the iteration times;
s7: updating the position of the wolf group;
s8: checking whether the updated wolf meets the preset constraint condition, if so, entering the next step, and if not, returning to the step S7;
s9: and adding 1 to the iteration times, judging whether the iteration times are larger than or equal to the maximum iteration times, if not, returning to the step S7, and if so, outputting the position of the alpha wolf, namely the target size of the battery energy storage system.
Further, the formula expression of the coding matrix in step S1 is:
in the formula (II)The optional size of the battery energy storage system in the microgrid is represented, X represents the battery energy storage system, n represents the number of transverse elements in the coding matrix, namely the number of the battery energy storage systems, and p represents the position of the battery energy storage system X.
Further, the formula for initializing the wolf pack position in step S2 is:
Position=Pmin+rand()*(Pmax-Pmin);
in the formula, PminRepresents the minimum value of the position of the wolf pack, PmaxRepresents the maximum value of the Position of the wolf pack, rand () represents a random function, and Position represents the random Position of the wolf pack.
Further, the preset constraint condition includes:
constraint equation 1:
0≤Pd(t)≤Pd;
in the formula, Pd is the standard output power of the battery energy storage system, and Pd (t) is the output power of the battery energy storage system at the time t;
constraint equation 2:
S1≤SOC(t)≤S2;
in the formula, S1 and S2 respectively represent a standard upper limit and a standard lower limit of the energy storage charge amount of the battery energy storage system, and soc (t) represents the energy storage state of charge of the battery energy storage system at time t;
constraint condition formula 3, that is, the relationship between the energy storage state of charge and the output power of the battery energy storage system needs to satisfy the equation:
in the formula, SOC (0) is the energy storage state of charge at the moment of fault occurrence, etafThe discharge efficiency of the battery energy storage system; edThe energy storage capacity of the battery energy storage system is obtained;
constraint equation 4:
where Δ f is the frequency deviation of the battery energy storage system and Δ fmaxThe maximum frequency deviation of the battery energy storage system;is the slip of the battery energy storage system,the maximum slip of the battery energy storage system;
constraint equation 4:
wherein H is an inertia time constant, KLAdjusting the comprehensive frequency of the load, wherein the delta P is the maximum possible power shortage fault disturbance quantity under the defense standard in the micro-grid, and the delta PdThe rated discharge power of the battery energy storage system.
Further, the formula expression of the fitness objective function in step S4 is as follows:
wherein N represents the frequency of sampling the battery energy storage system within a preset time period, fiRepresenting the frequency f of the battery energy storage system obtained by the ith sampling0And F is a fitness objective function value of the battery energy storage system.
Further, the step S7 of updating the position of the wolf pack includes the steps of:
s71: the distances between each wolf and alpha, beta and delta wolfs are obtained, and the obtaining formula is as follows:
in the formula,respectively representing the distance between the position of each wolf and the positions of alpha, beta and delta wolfs;respectively representing the current positions of alpha, beta and delta wolfs; are all random vectors, saidWhereinIs a random value between (0, 1);represents the position of the current wolf;
s72: updating the positions of the wolf groups according to the distances between each wolf and alpha, beta and delta wolfs respectively, wherein the formula is as follows:
in the formula,wherein,is a constant vector that decays from 2 to 0 with the number of iterations,is a random value between (0,1),is a random coefficient;the updated position of each wolf according to alpha, beta and delta wolfs,is the updated position of the wolf.
The invention also provides a battery energy storage system size division system based on the microgrid, which obtains the target size of the battery energy storage system through a wolf optimization algorithm so as to stabilize the frequency response of the microgrid when a fault occurs; the system comprises:
the coding matrix acquisition module is used for acquiring the energy storage system of each battery in the micro-grid according to the micro-gridObtaining a coding matrix from the parametric data, the elements of the coding matrixThe optional size of a battery energy storage system in the micro-grid is represented, X in elements represents the battery energy storage system, n in elements represents the number of the battery energy storage systems, and p represents the position of the battery energy storage system X;
the setting module is used for setting the number of wolves in the wolves, namely the number and the maximum iteration times of the battery energy storage system X, and the value range of the size parameter variable of the battery energy storage system, namely the range formed by the minimum value of the positions of the wolves and the maximum value of the positions of the wolves through the coding matrix, and initializing the positions of the wolves, namely the positions of the battery energy storage system X through the value range;
the judgment module is used for reinitializing the wolf pack positions through the setting module when each battery energy storage system does not meet the preset constraint condition;
the fitness objective function value acquisition module is used for acquiring a fitness objective function value corresponding to the position of each wolf through the fitness objective function when each battery energy storage system meets a preset constraint condition, namely a frequency deviation effect objective function value of the micro-grid after the micro-grid passes through the battery energy storage system for frequency modulation;
the dividing module is used for dividing the wolf group into alpha, beta, delta and omega according to the fitness objective function value, wherein the alpha wolf is the wolf with the optimal frequency modulation effect of the battery energy storage system, the beta wolf is the wolf with the suboptimal frequency modulation effect of the battery energy storage system, the delta wolf is the wolf with the third best frequency modulation effect of the battery energy storage system, and the omega wolf is the rest wolf;
the iteration module is used for initializing the iteration times and starting counting the iteration times;
the position updating module is used for updating the position of the wolf pack;
the checking module is used for updating the position of the wolf pack again through the position updating module when the updated wolf does not meet the preset constraint condition;
and the output module is used for adding 1 to the iteration times when the updated wolf meets the preset constraint condition, outputting the position of the alpha wolf when the current iteration times is greater than or equal to the maximum iteration times, namely the target size of the battery energy storage system, and updating the position of the wolf cluster through the position updating module when the current iteration times is less than the maximum iteration times.
Compared with the prior art, the invention at least has the following beneficial effects:
(1) the method comprises the steps of obtaining an encoding matrix according to parameter data of a micro-grid and each battery energy storage system in the micro-grid, setting input parameters (the number, the maximum iteration times and the value range of a wolf) of a wolf optimization algorithm according to the encoding matrix, and initializing the position of a wolf group; the method comprises the steps of dividing a wolf group through a fitness objective function, outputting the position of the last alpha wolf as the target size of a battery energy storage system after entering the iteration of a wolf optimization algorithm, embedding the wolf optimization algorithm into a battery energy storage system size dividing method to obtain a battery energy storage system with a proper target size, and achieving the purpose of reasonably controlling the cost of the battery energy storage system while providing enough frequency support for a micro-grid to stabilize the frequency response of the micro-grid and keeping the safety of the power system when the micro-grid fails and the frequency is unstable due to imbalance between power generation and load requirements;
(2) in the size division method, whether each battery energy storage system meets the preset constraint condition or not is judged, and after the position of the wolf pack is updated, whether the updated wolf meets the preset constraint condition or not is checked, so that the accuracy of the target size is further improved;
(3) according to the invention, the battery energy storage system with the optimal size (namely the target size) is obtained to provide frequency support for the micro-grid, so that the conventional power generation burden participating in inertial response and primary frequency control is reduced while high-frequency fluctuation is reduced;
(4) the battery energy storage system with the optimal size obtained by the invention stabilizes the frequency response of the micro-grid and improves the robustness of the power system.
Drawings
Fig. 1 is a flow chart of a method for partitioning the size of a battery energy storage system based on a microgrid;
FIG. 2 is a system block diagram of a microgrid-based battery energy storage system sizing system;
FIG. 3 is a gray wolf level diagram in the gray wolf optimization algorithm.
Detailed Description
The following are specific embodiments of the present invention and are further described with reference to the accompanying drawings, but the present invention is not limited to these embodiments.
Example one
Battery Energy Storage Systems (BESS) are widely used in electrical grids to provide frequency support by injecting transient power into the microgrid, however, overdimensioning of battery energy storage systems increases cost, while undersizing results in insufficient battery capacity and thus poor performance. In order to solve the technical problem, as shown in fig. 1, the invention provides a battery energy storage system size division method based on a microgrid, which obtains a target size of the battery energy storage system through a wolf optimization algorithm so as to stabilize the frequency response of the microgrid when a fault occurs; the size division method comprises the following steps:
s1: acquiring a coding matrix according to the parameter data of the microgrid and each battery energy storage system in the microgrid, wherein elements in the coding matrixThe optional size of a battery energy storage system in the micro-grid is represented, X in elements represents the battery energy storage system, n in elements represents the number of the battery energy storage systems, and p represents the position of the battery energy storage system X;
the formula expression of the coding matrix in step S1 is:
in the formula (II)Indicating selectable size of battery energy storage system in microgridX denotes the battery energy storage system, n denotes the number of transverse elements in the coding matrix, i.e. the number of battery energy storage systems, and p denotes the position of the battery energy storage system X.
S2: setting the number of wolves in a wolve group, namely the number of the battery energy storage systems X, the maximum iteration times and the value range of the size parameter variable of the battery energy storage system, namely the range formed by the minimum value of the position of the wolve group and the maximum value of the position of the wolve group through the coding matrix, and initializing the position of the wolve group, namely the position of the battery energy storage system X through the value range;
the value range of the size parameter variable of the battery energy storage system is determined by the optional size of the battery energy storage system in the microgrid.
The formula for initializing the wolf pack position in step S2 is:
Position=Pmin+rand()*(Pmax-Pmin);
in the formula, PminRepresents the minimum value of the position of the wolf pack, PmaxRepresents the maximum value of the wolf pack Position, rand () represents a random function, and Position represents the wolf pack random Position.
S3: judging whether each battery energy storage system meets a preset constraint condition or not; if yes, the next step is carried out, if not, the step S2 is returned to reinitialize the wolf pack position;
s4: acquiring a fitness objective function value corresponding to the position of each wolf through the fitness objective function, namely acquiring a frequency deviation effect objective function value of the micro-grid after the micro-grid is subjected to frequency modulation through a battery energy storage system;
the formula expression of the fitness objective function in step S4 is as follows:
wherein N represents the number of times the frequency of the battery energy storage system is sampled within a preset time period, fiRepresenting the frequency f of the battery energy storage system obtained by the ith sampling0And F is a fitness objective function value of the battery energy storage system.
The battery energy storage system is subjected to N times of frequency sampling within a preset time period, and then the fitness objective function value corresponding to the battery energy storage system is obtained through the fitness objective function designed by the invention, so that wolves of different levels are accurately divided, and the target size obtaining precision is improved.
For the gray wolf optimization algorithm (GWO), it should be noted that:
the sirius is a canine, considered a top predator, who is at the top of the biosphere food chain. The wolf is most like the group, and there are 5-12 wolfs on average in each group. Of particular interest are those having a very strict hierarchy of social levels, as shown in fig. 3. The first level of the pyramid is the leader in the population, called α. In the wolf pack, α is an individual with management ability, and is mainly responsible for matters related to various decisions in hunting, sleeping time and place, food distribution and other groups. The second level of the pyramid is the brain team of α, called β, which is primarily responsible for assisting α in making decisions. When alpha of the whole wolf group is vacant, beta will take over the position of alpha. The dominance of beta in the wolf group is next to alpha, and the wolf group gives the command of alpha to other members and feeds back the execution conditions of the other members to alpha to play the role of a bridge. The third layer of the pyramid is delta, which follows decision commands of alpha and beta and is mainly responsible for investigation, sentry, nursing and other matters. Poorly adapted α and β will also be reduced to δ. The bottom layer of the pyramid is omega, which is mainly responsible for the balance of the relationships within the population. Wolves have several tricky processes such as tracking, chasing and embracing, all of which must be performed accurately to create conditions for a perfect hunt.
GWO has the advantage of having a high level of robustness and the ability to handle high dimensional, non-linear and non-convex problems, and is an ideal tool for achieving the target dimensions of the present invention.
S5: dividing a wolf group into alpha, beta, delta and omega according to the fitness objective function value, wherein the alpha wolf is the wolf with the optimal frequency modulation effect of the battery energy storage system, the beta wolf is the wolf with the suboptimal frequency modulation effect of the battery energy storage system, the delta wolf is the wolf with the third best frequency modulation effect of the battery energy storage system, and the omega wolf is the rest wolf;
s6: initializing iteration times and starting counting the iteration times;
the initial value of the number of iterations in this example is 1.
S7: updating the position of the wolf group;
the updating of the position of the wolf pack in the step S7 includes the steps of:
s71: obtaining the distance between each wolf and alpha, beta and delta wolfs respectively, wherein the obtaining formula is as follows:
in the formula,respectively representing the distance between the position of each wolf and the positions of alpha, beta and delta wolfs;respectively representing the current positions of alpha, beta and delta wolfs; are all random vectors, saidWhereinIs a random value between (0, 1);represents the position of the current wolf;
s72: updating the positions of the wolf clusters according to the distances between each wolf and alpha, beta and delta wolf respectively, wherein the formula is as follows:
in the formula,wherein,is a constant vector that decays from 2 to 0 with the number of iterations,is a random value between (0,1),is a random coefficient;for each wolf, the updated position is based on alpha, beta, delta wolf,to updateThe position of the rear wolf.
S8: checking whether the updated wolf meets the preset constraint condition, if so, entering the next step, and if not, returning to the step S7;
in the size division method, whether each battery energy storage system meets the preset constraint condition or not is judged, and after the position of the wolf pack is updated, whether the updated wolf meets the preset constraint condition or not is checked, so that the accuracy of the target size is further improved.
S9: and adding 1 to the iteration times, judging whether the iteration times are larger than or equal to the maximum iteration times, if not, returning to the step S7, and if so, outputting the position of the alpha wolf, namely the target size of the battery energy storage system.
The preset constraint conditions comprise:
constraint equation 1:
0≤Pd(t)≤Pd;
in the formula, t represents time t, Pd is the standard output power of the battery energy storage system, and Pd (t) is the output power of the battery energy storage system at time t;
constraint equation 2:
S1≤SOC(t)≤S2;
in the formula, S1 and S2 respectively represent a standard upper limit and a standard lower limit of the energy storage charge amount of the battery energy storage system, and soc (t) represents the energy storage state of charge of the battery energy storage system at time t;
constraint condition formula 3, that is, the relationship between the energy storage state of charge and the output power of the battery energy storage system needs to satisfy the equation:
in the formula, SOC (0) is the energy storage state of charge at the moment of fault occurrence, etafThe discharge efficiency of the battery energy storage system; edIs electricityThe energy storage capacity of the battery energy storage system is large;
constraint equation 4:
where Δ f is the frequency deviation of the battery energy storage system and Δ fmaxThe maximum frequency deviation of the battery energy storage system;is the slip of the battery energy storage system,the maximum slip of the battery energy storage system;
constraint equation 4:
wherein H is an inertia time constant, KLThe comprehensive frequency regulation coefficient of the load is adopted, the delta P is the maximum possible power shortage fault disturbance quantity under the defense standard in the micro-grid, and the delta PdThe rated discharge power of the battery energy storage system.
The method comprises the steps of obtaining an encoding matrix according to parameter data of a micro-grid and each battery energy storage system in the micro-grid, setting input parameters (the number, the maximum iteration times and the value range of a wolf) of a wolf optimization algorithm according to the encoding matrix, and initializing the position of a wolf group; the wolf pack is divided through a fitness objective function, the position of the last alpha wolf is output as the target size of the battery energy storage system after iteration of a wolf optimization algorithm is entered, the wolf optimization algorithm is embedded into the battery energy storage system size dividing method, the battery energy storage system with the proper target size is obtained, and the purposes that when the micro-grid fails, and frequency instability is caused by imbalance between power generation and load requirements, enough frequency support is provided for the micro-grid to stabilize the frequency response of the micro-grid, the safety of the power system is kept, the power system can stably run under normal and abnormal conditions, and meanwhile, the cost of the battery energy storage system is reasonably controlled are achieved.
Example two
The invention also provides a battery energy storage system size division system based on the microgrid, which obtains the target size of the battery energy storage system through a wolf optimization algorithm so as to stabilize the frequency response of the microgrid when a fault occurs; the system comprises:
the coding matrix acquisition module is used for acquiring a coding matrix according to the parameter data of the microgrid and each battery energy storage system in the microgrid, and elements in the coding matrixThe optional size of a battery energy storage system in the micro-grid is represented, X in elements represents the battery energy storage system, n in elements represents the number of the battery energy storage systems, and p represents the position of the battery energy storage system X;
the setting module is used for setting the number of wolves in the wolves, namely the number and the maximum iteration times of the battery energy storage system X, and the value range of the size parameter variable of the battery energy storage system, namely the range formed by the minimum value of the positions of the wolves and the maximum value of the positions of the wolves through the coding matrix, and initializing the positions of the wolves, namely the positions of the battery energy storage system X through the value range;
the judgment module is used for reinitializing the wolf pack positions through the setting module when each battery energy storage system does not meet the preset constraint condition;
the fitness objective function value acquisition module is used for acquiring a fitness objective function value corresponding to the position of each wolf through the fitness objective function when each battery energy storage system meets a preset constraint condition, namely a frequency deviation effect objective function value of the micro-grid after the micro-grid passes through the battery energy storage system for frequency modulation;
the dividing module is used for dividing the wolf group into alpha, beta, delta and omega according to the fitness objective function value, wherein the alpha wolf is the wolf with the optimal frequency modulation effect of the battery energy storage system, the beta wolf is the wolf with the suboptimal frequency modulation effect of the battery energy storage system, the delta wolf is the wolf with the third best frequency modulation effect of the battery energy storage system, and the omega wolf is the rest wolf;
the iteration module is used for initializing the iteration times and starting counting the iteration times;
the position updating module is used for updating the position of the wolf pack;
the checking module is used for updating the position of the wolf pack again through the position updating module when the updated wolf does not meet the preset constraint condition;
and the output module is used for adding 1 to the iteration times when the updated wolf meets the preset constraint condition, outputting the position of the alpha wolf when the current iteration times is greater than or equal to the maximum iteration times, namely the target size of the battery energy storage system, and updating the position of the wolf cluster through the position updating module when the current iteration times is less than the maximum iteration times.
The system provides frequency support for the microgrid by acquiring a battery energy storage system with the optimal size (namely the target size), so that the conventional power generation burden of participating in inertial response and primary frequency control is reduced while high-frequency fluctuation is reduced.
It should be noted that all the directional indicators (such as up, down, left, right, front, and rear … …) in the embodiment of the present invention are only used to explain the relative position relationship between the components, the movement situation, etc. in a specific posture (as shown in the drawing), and if the specific posture is changed, the directional indicator is changed accordingly.
Moreover, descriptions of the present invention as relating to "first," "second," "a," etc. are for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicit ly indicating a number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of the present invention, "a plurality" means at least two, e.g., two, three, etc., unless specifically limited otherwise.
In the present invention, unless otherwise expressly stated or limited, the terms "connected," "secured," and the like are to be construed broadly, and for example, "secured" may be a fixed connection, a removable connection, or an integral part; can be mechanically or electrically connected; they may be directly connected or indirectly connected through intervening media, or they may be connected internally or in any other suitable relationship, unless expressly stated otherwise. The specific meanings of the above terms in the present invention can be understood by those skilled in the art according to specific situations.
In addition, the technical solutions in the embodiments of the present invention may be combined with each other, but it must be based on the realization of those skilled in the art, and when the technical solutions are contradictory or cannot be realized, such a combination of technical solutions should not be considered to exist, and is not within the protection scope of the present invention.
Claims (7)
1. A battery energy storage system size division method based on a micro-grid is characterized in that a target size of the battery energy storage system is obtained through a wolf optimization algorithm so as to stabilize the frequency response of the micro-grid when a fault occurs; the size division method comprises the following steps:
s1: acquiring a coding matrix according to parameter data of each battery energy storage system in the microgrid and the microgrid, wherein elements in the coding matrixThe optional size of a battery energy storage system in the micro-grid is represented, X in elements represents the battery energy storage system, n in elements represents the number of the battery energy storage systems, and p represents the position of the battery energy storage system X;
s2: setting the number of wolves in a wolve group, namely the number of the battery energy storage systems X, the maximum iteration times and the value range of the size parameter variable of the battery energy storage system, namely the range formed by the minimum value of the position of the wolve group and the maximum value of the position of the wolve group through the coding matrix, and initializing the position of the wolve group, namely the position of the battery energy storage system X through the value range;
s3: judging whether each battery energy storage system meets a preset constraint condition or not; if yes, the next step is carried out, if not, the step S2 is returned to reinitialize the wolf pack position;
s4: acquiring a fitness objective function value corresponding to the position of each wolf through the fitness objective function, namely a frequency deviation effect objective function value of the micro-grid after frequency modulation through the battery energy storage system;
s5: dividing a wolf group into alpha, beta, delta and omega according to the fitness objective function value, wherein the alpha wolf is the wolf with the optimal frequency modulation effect of the battery energy storage system, the beta wolf is the wolf with the suboptimal frequency modulation effect of the battery energy storage system, the delta wolf is the wolf with the third best frequency modulation effect of the battery energy storage system, and the omega wolf is the rest wolf;
s6: initializing iteration times and starting counting the iteration times;
s7: updating the position of the wolf group;
s8: checking whether the updated wolf meets the preset constraint condition, if so, entering the next step, and if not, returning to the step S7;
s9: and adding 1 to the iteration times, judging whether the iteration times are larger than or equal to the maximum iteration times, if not, returning to the step S7, and if so, outputting the position of the alpha wolf, namely the target size of the battery energy storage system.
2. The microgrid-based battery energy storage system size division method as claimed in claim 1, wherein the formula expression of the coding matrix in the step S1 is as follows:
in the formula (II)The optional size of the battery energy storage system in the micro-grid is represented, X represents the battery energy storage system, n represents the number of transverse elements in the coding matrix, namely the number of the battery energy storage systems, and p represents the position of the battery energy storage system X.
3. The microgrid-based battery energy storage system size division method as claimed in claim 2, wherein the formula for initializing the wolf pack position in step S2 is as follows:
Position=Pmin+rand()*(Pmax-Pmin);
in the formula, PminRepresents the minimum value of the wolf group position, PmaxRepresents the maximum value of the Position of the wolf pack, rand () represents a random function, and Position represents the random Position of the wolf pack.
4. The microgrid-based battery energy storage system size division method of claim 1, wherein the preset constraint conditions include:
constraint equation 1:
0≤Pd(t)≤Pd;
in the formula, Pd is the standard output power of the battery energy storage system, and Pd (t) is the output power of the battery energy storage system at the time t;
constraint equation 2:
S1≤SOC(t)≤S2;
in the formula, S1 and S2 respectively represent a standard upper limit and a standard lower limit of the energy storage charge amount of the battery energy storage system, and soc (t) represents the energy storage state of charge of the battery energy storage system at time t;
constraint condition formula 3, that is, the relationship between the energy storage state of charge and the output power of the battery energy storage system needs to satisfy the equation:
in the formula, SOC (0) is the energy storage state of charge at the moment of fault occurrence, etafThe discharge efficiency of the battery energy storage system; edThe energy storage capacity of the battery energy storage system is set;
constraint equation 4:
where Δ f is the frequency offset of the battery energy storage systemDifference,. DELTA.fmaxThe maximum frequency deviation of the battery energy storage system;is the slip of the battery energy storage system,the maximum slip of the battery energy storage system;
constraint equation 4:
wherein H is an inertia time constant, KLThe comprehensive frequency regulation coefficient of the load is adopted, the delta P is the maximum possible power shortage fault disturbance quantity under the defense standard in the micro-grid, and the delta PdThe rated discharge power of the battery energy storage system.
5. The microgrid-based battery energy storage system size division method of claim 3, wherein the formula expression of the fitness objective function in the step S4 is as follows:
wherein N represents the number of times the frequency of the battery energy storage system is sampled within a preset time period, fiRepresenting the frequency f of the battery energy storage system obtained by the ith sampling0And F is a fitness objective function value of the battery energy storage system.
6. The microgrid-based battery energy storage system sizing method of claim 5, wherein the step S7 of updating the location of the wolf pack comprises the steps of:
s71: the distances between each wolf and alpha, beta and delta wolfs are obtained, and the obtaining formula is as follows:
in the formula,respectively representing the distance between the position of each wolf and the positions of alpha, beta and delta wolfs;respectively representing the current positions of alpha, beta and delta wolfs; are all random vectors, saidWhereinIs a random value between (0, 1);indicating the position of the current wolf;
s72: updating the positions of the wolf clusters according to the distances between each wolf and alpha, beta and delta wolf respectively, wherein the formula is as follows:
7. A battery energy storage system size division system based on a micro-grid is characterized in that a target size of the battery energy storage system is obtained through a wolf optimization algorithm so as to stabilize the frequency response of the micro-grid when a fault occurs; the system comprises:
the coding matrix acquisition module is used for acquiring a coding matrix according to the parameter data of the micro-grid and each battery energy storage system in the micro-grid, and elements in the coding matrixThe optional size of a battery energy storage system in the micro-grid is represented, X in elements represents the battery energy storage system, n in elements represents the number of the battery energy storage systems, and p represents the position of the battery energy storage system X;
the setting module is used for setting the number of wolves in the wolves, namely the number and the maximum iteration times of the battery energy storage system X, and the value range of the size parameter variable of the battery energy storage system, namely the range formed by the minimum value of the positions of the wolves and the maximum value of the positions of the wolves through the coding matrix, and initializing the positions of the wolves, namely the positions of the battery energy storage system X through the value range;
the judgment module is used for reinitializing the wolf pack positions through the setting module when each battery energy storage system does not meet the preset constraint condition;
the fitness objective function value acquisition module is used for acquiring a fitness objective function value corresponding to the position of each wolf through the fitness objective function when each battery energy storage system meets a preset constraint condition, namely a frequency deviation effect objective function value of the micro-grid after the micro-grid passes through the battery energy storage system for frequency modulation;
the dividing module is used for dividing the wolf group into alpha, beta, delta and omega according to the fitness objective function value, wherein the alpha wolf is the wolf with the optimal frequency modulation effect of the battery energy storage system, the beta wolf is the wolf with the suboptimal frequency modulation effect of the battery energy storage system, the delta wolf is the wolf with the third best frequency modulation effect of the battery energy storage system, and the omega wolf is the rest wolf;
the iteration module is used for initializing the iteration times and starting counting the iteration times;
the position updating module is used for updating the position of the wolf pack;
the checking module is used for updating the position of the wolf pack again through the position updating module when the updated wolf does not meet the preset constraint condition;
and the output module is used for adding 1 to the iteration times when the updated wolf meets the preset constraint condition, outputting the position of the alpha wolf when the current iteration times is greater than or equal to the maximum iteration times, namely the target size of the battery energy storage system, and updating the position of the wolf cluster through the position updating module when the current iteration times is less than the maximum iteration times.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210409307.1A CN114744658B (en) | 2022-04-19 | 2022-04-19 | Battery energy storage system size division method and system based on micro-grid |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210409307.1A CN114744658B (en) | 2022-04-19 | 2022-04-19 | Battery energy storage system size division method and system based on micro-grid |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114744658A true CN114744658A (en) | 2022-07-12 |
CN114744658B CN114744658B (en) | 2024-05-14 |
Family
ID=82281339
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210409307.1A Active CN114744658B (en) | 2022-04-19 | 2022-04-19 | Battery energy storage system size division method and system based on micro-grid |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114744658B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170102675A1 (en) * | 2015-10-08 | 2017-04-13 | Johnson Controls Technology Company | Photovoltaic energy system with preemptive ramp rate control |
US20170104337A1 (en) * | 2015-10-08 | 2017-04-13 | Johnson Controls Technology Company | Photovoltaic energy system with value function optimization |
WO2017062908A1 (en) * | 2015-10-08 | 2017-04-13 | Johnson Controls Technology Company | Photovoltaic energy system with solar intensity prediction |
CN106896716A (en) * | 2017-04-17 | 2017-06-27 | 华北电力大学(保定) | Micro-capacitance sensor alternating current-direct current section transverter pid parameter optimization method based on grey wolf algorithm |
CN110401209A (en) * | 2019-05-14 | 2019-11-01 | 东华大学 | The energy management method of peak load shifting based on more random composite optimization grey wolf algorithms |
CN113471989A (en) * | 2021-07-02 | 2021-10-01 | 东北电力大学 | Intelligent micro-grid secondary frequency control method based on wolf optimization method |
CN114172176A (en) * | 2021-11-19 | 2022-03-11 | 国网内蒙古东部电力有限公司电力科学研究院 | Flywheel energy storage array control method and system applied to wind generating set |
-
2022
- 2022-04-19 CN CN202210409307.1A patent/CN114744658B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170102675A1 (en) * | 2015-10-08 | 2017-04-13 | Johnson Controls Technology Company | Photovoltaic energy system with preemptive ramp rate control |
US20170104337A1 (en) * | 2015-10-08 | 2017-04-13 | Johnson Controls Technology Company | Photovoltaic energy system with value function optimization |
WO2017062908A1 (en) * | 2015-10-08 | 2017-04-13 | Johnson Controls Technology Company | Photovoltaic energy system with solar intensity prediction |
CN106896716A (en) * | 2017-04-17 | 2017-06-27 | 华北电力大学(保定) | Micro-capacitance sensor alternating current-direct current section transverter pid parameter optimization method based on grey wolf algorithm |
CN110401209A (en) * | 2019-05-14 | 2019-11-01 | 东华大学 | The energy management method of peak load shifting based on more random composite optimization grey wolf algorithms |
CN113471989A (en) * | 2021-07-02 | 2021-10-01 | 东北电力大学 | Intelligent micro-grid secondary frequency control method based on wolf optimization method |
CN114172176A (en) * | 2021-11-19 | 2022-03-11 | 国网内蒙古东部电力有限公司电力科学研究院 | Flywheel energy storage array control method and system applied to wind generating set |
Non-Patent Citations (2)
Title |
---|
杨灿;王勇;: "基于灰狼算法的分布式混合能源优化配置方法", 微型电脑应用, no. 12, 19 December 2018 (2018-12-19) * |
牛成玉;傅仲文;: "基于灰狼优化算法的微电网并网控制器参数优化设计", 电器与能效管理技术, no. 20, 30 October 2017 (2017-10-30) * |
Also Published As
Publication number | Publication date |
---|---|
CN114744658B (en) | 2024-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101919134B (en) | Event-based control system for wind turbine generators and control method thereof | |
EP2921698A1 (en) | System and method for automatic generation control in wind farms | |
KR101778309B1 (en) | Lead acid storage battery and lead acid storage battery system for a natural energy using system | |
CN112103980B (en) | Energy management method of hybrid energy storage system combining AGC frequency modulation of thermal power generating unit | |
CN104779607B (en) | One of direct-current micro-grid distributed and coordinated control method and system | |
US9843190B2 (en) | Power system management device and method | |
JP6899807B2 (en) | Aggregator device, program and supply / demand adjustment method to control the charge / discharge power of each storage battery | |
WO2017006371A1 (en) | Renewable power system and sizing method for controllable plant associated with renewable power system | |
CN111682563A (en) | Micro-grid intelligent frequency control method based on electric energy storage system | |
CN113054688A (en) | Renewable energy cluster output distributed coordination control method considering uncertainty | |
CN116231765B (en) | Virtual power plant output control method | |
JP2015027158A (en) | Power storage device and charge/discharge method for power storage device | |
CN111049170B (en) | Peak-shaving energy storage system considering temperature loss, operation optimization method and device | |
CN114744658A (en) | Battery energy storage system size division method and system based on micro-grid | |
CN112186764B (en) | Access optimization method and device for power distribution network equipment and electronic equipment | |
CN114006410A (en) | Large-scale offshore wind power access point optimization method based on opportunity constraint planning | |
US20190319464A1 (en) | Method for operating ess at optimal efficiency | |
CN114744657A (en) | Battery energy storage system size optimization method based on micro-grid | |
CN114330938B (en) | Distributed energy storage planning method and system for power distribution network | |
CN107230994B (en) | A kind of method and device of control energy-storage system power output | |
CN108336769B (en) | Method and system for inhibiting active power output error of distributed power supply under communication interference | |
WO2024057581A1 (en) | Charge-state calculation device | |
US12000901B2 (en) | Battery state estimation based on multiple rates of hysteresis transit | |
CN112751343B (en) | Power distribution network double-layer optimization method based on distributed cooperative control | |
WO2023033685A1 (en) | Support of an electrical power grid |
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 |