CN115276055A - Energy storage configuration method and system based on power grid frequency spatial distribution characteristics - Google Patents
Energy storage configuration method and system based on power grid frequency spatial distribution characteristics Download PDFInfo
- Publication number
- CN115276055A CN115276055A CN202211054711.8A CN202211054711A CN115276055A CN 115276055 A CN115276055 A CN 115276055A CN 202211054711 A CN202211054711 A CN 202211054711A CN 115276055 A CN115276055 A CN 115276055A
- Authority
- CN
- China
- Prior art keywords
- energy storage
- power
- fault
- unbalanced
- capacity
- 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
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- 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
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Power Engineering (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses an energy storage configuration method and system based on the frequency space distribution characteristics of a power grid, and belongs to the technical field of electrical engineering. The method comprises the following steps: the method comprises the steps of calculating the total power capacity of energy storage configuration and the system fault probability in the population iteration process based on a non-inferior ranking genetic algorithm to obtain the optimal total power capacity of energy storage configuration and the optimal system fault probability, wherein the total power capacity of energy storage configuration is used as an economic index for measuring the power system, the fault probability is used as a frequency stability index for measuring the power system, and an energy storage capacity configuration scheme which is most consistent with the expectation of an investor is obtained according to the economic performance of the energy storage configuration and the system frequency stability weight coefficient of the investor. The method can give consideration to both economy and power grid frequency stability, quickly realize the site selection and constant volume of the inertial support energy storage system, ensure the transient stability of the power grid frequency, and solve the problem that the existing energy storage configuration method is difficult to give consideration to both frequency space distribution characteristics and solving difficulty.
Description
Technical Field
The invention belongs to the technical field of electrical engineering, and particularly relates to an energy storage configuration method and system based on power grid frequency spatial distribution characteristics.
Background
In recent years, as new energy generating sets are continuously connected to a power system, the starting proportion of the traditional synchronous generating sets is continuously reduced, the inertia of the system is reduced, and the frequency stability faces a great challenge. And the energy storage can participate in the inertial support link of the frequency of the power system through a proper inertial support control mode by virtue of the bidirectional and rapid adjusting characteristic of the energy storage, so that the frequency transient stability of the power system is improved. However, the price of energy storage is relatively high, and an energy storage configuration method for improving the inertia of the system is the key for ensuring the stability and the economy of the system.
The existing energy storage configuration method mainly takes the frequency of an inertia central point as the frequency of the whole network, neglects the influence of the power grid, and weights and equates all generators of the system into an equivalent unit according to the inertia size, so that the dynamic frequency change of the whole system can be well reflected, the calculation amount is small, the configuration is easy, but the model cannot reflect the dynamic frequency difference of different areas of the system, and when the energy storage inertia support is configured, the model can only analyze the influence of the energy storage capacity on the frequency support, cannot analyze the effect of the energy storage access point on the frequency support, and has certain limitation. In order to analyze the dynamic influence of the units in different areas, an all-state analysis method is mainly adopted at present, numerical iteration solution is carried out by establishing detailed mathematical models and network equations of the units, but the method has large calculation amount, the simulation mode adopted for optimal configuration has long time consumption, and the method is not suitable for being used as a frequency model for optimal configuration. In addition, the configuration methods in the prior art do not consider the economic efficiency and stability index of the power system at the same time.
Disclosure of Invention
Aiming at the defects and the improvement requirements of the prior art, the invention provides an energy storage configuration method and an energy storage configuration system based on the frequency-space distribution characteristics of a power grid, and aims to simultaneously consider the frequency-space distribution characteristics of energy storage and reduce the calculation amount of energy storage configuration.
In order to achieve the above object, according to an aspect of the present invention, there is provided an energy storage configuration method based on grid frequency spatial distribution characteristics, including:
s1, setting upper and lower limits of energy storage power capacity of a new energy source side and a network side and a node range set No accessed by a network side centralized energy storage power station e ；
S2, collecting No by upper and lower limits of energy storage power capacity and node range e For input, a genetic algorithm is used to randomly generate N pop Configuring an initial population by each energy storage;
s3, calculating the total energy storage power capacity configured for each individual in the initial population; and calculating an N-1 fault accident set A based on a direct current power flow frequency response model N-1 The frequency change rate of the lower non-fault synchronous generator is limited according to the frequency change rate to obtain the fault switching rate of the system;
s4, according to the total energy storage power capacity and the fault tripping rate, performing non-dominated sorting and congestion degree sorting on the energy storage configuration population to obtain a Pareto grade and a congestion degree;
s5, selecting an energy storage configuration parent according to the Pareto grade and the crowding degree, generating energy storage configuration offspring through cross variation, calculating the total energy storage power capacity and the fault probability of individuals in the energy storage configuration offspring, and combining two indexes of the offspring and the parent to obtain the current total energy storage power capacity and the system fault probability; if the iteration times do not reach the preset value, taking the current total energy storage power capacity and the system failure rate as input, and returning to the step S4, otherwise, taking a set consisting of individuals with Pareto grades of 1 in the current total energy storage power capacity and the system failure rate as a Pareto optimal solution set to output;
s6, normalizing the Pareto optimal solution set, and presetting a weight coefficient lambda 1 And λ 2 And obtaining an optimal energy storage configuration scheme.
Further, the total energy storage power capacity y configured for each individual in the population 1i Comprises the following steps:
wherein i is more than or equal to 1 and less than or equal to N pop When j is more than or equal to 1 and less than or equal to N e When x is ij The energy storage power capacity of the ith energy storage configuration result in the population configured on the j new energy sources is represented, and S is satisfied eminj ≤x ij ≤S emaxj (ii) a When j = N e At +1, x ij The energy storage power capacity of the ith energy storage configuration result in the population is centrally configured on the network side and meets the requirement of S eminj ≤x ij ≤S emaxj ；S emaxj Represents the upper limit of the energy storage power capacity of the new energy source side and the network side, S eminj Representing the new energy source side and the net side energy storage power capacity lower limit.
Further, calculating the probability of individual fault switching in the population, including:
under the constraint of upper and lower limits of energy storage power capacity, calculating a direct current power flow frequency response model N-1 fault accident set A N-1 Fault generating unbalanced power vector Δ P in (1) d Generating an unbalanced power vector Δ P based on the fault d The relation satisfied between the electromagnetic power and the unbalanced electromagnetic power is calculated, and the unbalanced electromagnetic power delta P born by the stored energy in the power system is calculated E Unbalanced electromagnetic power Δ P borne by non-faulty synchronous generator nodes G ；
Judging the unbalanced electromagnetic power delta P born by the stored energy of all individuals in the population E Whether the configured stored energy power is exceededAnd if the unbalanced power borne by the mth stored energy exceeds the configured capacity limit, merging the mth stored energy into the fault-generated disturbance and load node set d as a disturbance node, and recalculating the unbalanced power delta P borne by all the stored energy except the mth stored energy E ' unbalanced power DeltaP ' borne by non-faulty synchronous generator ' G Again determine Δ P' E Whether the configured capacity limit is exceeded or not is judged until all the unbalanced power born by the stored energy meets the configured capacity limit, and the final unbalanced power born by the non-fault synchronous generator is obtained; wherein m is more than or equal to 1 and less than or equal to N e +1，N e The energy storage number of the inertial supports on the new energy side;
and calculating the frequency change rate of the non-fault synchronous generator according to the final unbalanced power born by the non-fault synchronous generator, and limiting according to the frequency change rate to obtain the fault switching rate.
Further, unbalanced electromagnetic power borne by energy storage and non-fault synchronous generator nodes in the power system and a fault-generated unbalanced power vector delta P d Satisfies the relationship:
wherein, Δ P GE ＝[ΔP E ，ΔP G ] T ，
In order to only consider a power grid susceptance matrix of power grid inductance, subscript GE represents an energy storage and non-fault synchronous generator node set, and subscript d represents disturbance generated by faults and a load node set; b is GE-GE 、B GE-d 、B d-GE And B d-d The block matrixes are all power grid susceptance matrixes;
is B d-d The inverse matrix of (c).
Further, the frequency change rate of the non-fault synchronous generator and the unbalanced power borne by the non-fault synchronous generator meet the following conditions:
wherein, Δ P G Unbalanced electromagnetic power vector, Δ P, assumed for a non-faulty synchronous generator node G ＝[ΔP G1 ，ΔP G2 ，…，ΔP GNg ] T Ng is the number of non-fault synchronous generators in the power system; f. of 1 …f Ng For each non-faulted synchronous generator rotor frequency; j is a unit of 1 …J Ng Is the time constant of the inertia of the rotor,
representing the derivative operator.
Further, the probability of fail-over y 2i Comprises the following steps:
wherein N is AN-1 Indicates the number of N-1 faults, S a Capacity of the a-th non-faulty synchronous generator, e ab Indicating whether the frequency change rate of the a-th non-fault generator exceeds the limit in the b-th fault scene, and if the frequency change rate of the generator exceeds the limit, e ab =1, otherwise, e ab ＝0。
Further, in step S6, the two normalized indexes of the Pareto optimal solution set are weighted and then the minimum value is obtained, so as to obtain the weight coefficient λ of the two indexes 1 And λ 2 。
According to another aspect of the present invention, there is provided an energy storage configuration system based on grid frequency spatial distribution characteristics, including:
a parameter setting module for setting the upper and lower limits of the energy storage power capacity of the new energy source side and the network side and the node range set No accessed by the network side centralized energy storage power station e ；
An initial population generation module for storing energy and powerUpper and lower limit of quantity and node range set No e For input, an NSGA-II algorithm is adopted to randomly generate N pop Configuring an initial population by each energy storage;
the energy storage total power capacity and fault probability switching calculation module comprises an energy storage total power capacity calculation unit and a fault probability switching calculation unit, wherein the energy storage total power capacity calculation unit is used for calculating the energy storage total power capacity configured for each individual in the initial population; the fault probability switching calculation unit is used for calculating an N-1 fault accident set A based on a direct current power flow frequency response model N-1 The frequency change rate of the lower non-fault synchronous generator is limited according to the frequency change rate to obtain the fault switching rate of the system;
the sorting module is used for carrying out non-dominated sorting and congestion degree sorting on the energy storage configuration population according to two indexes of the total energy storage power capacity and the fault tripping rate to obtain a Pareto grade and a congestion degree;
the Pareto optimal solution set generation module is used for selecting an energy storage configuration parent according to the Pareto grade and the crowding degree, generating energy storage configuration offspring through cross variation, respectively calculating the energy storage total power capacity and the fault probability of an individual in the energy storage configuration offspring through an energy storage total power capacity calculation unit and a fault probability calculation unit, and combining two indexes of the offspring and the parent to obtain the current energy storage total power capacity and the current system fault probability; if the iteration times do not reach the preset value, the current total energy storage power capacity and the system fault probability are used as input, the step S4 is returned, otherwise, a set formed by individuals with Pareto grades of 1 in the current total energy storage power capacity and the system fault probability is used as a Pareto optimal solution set to be output;
an optimal energy storage configuration scheme generation module, configured to normalize the Pareto optimal solution set, and preset a weight coefficient λ 1 And λ 2 And obtaining an optimal energy storage configuration scheme.
Further, the failure tripping rate calculation unit includes:
the unbalanced electromagnetic power calculation node is used for calculating a direct current power flow frequency response model N-1 fault accident set A under the constraint of upper and lower limits of energy storage power capacity under the condition of not counting N-1 Fault generating unbalanced power vector Δ P in (1) d Generating an unbalanced power vector Δ P based on the fault d The relation satisfied between the electromagnetic power and the unbalanced electromagnetic power is calculated, and the unbalanced electromagnetic power delta P born by the stored energy in the power system is calculated E Unbalanced electromagnetic power Δ P borne by non-faulty synchronous generator nodes G ；
The unbalanced electromagnetic power updating node is used for judging the unbalanced electromagnetic power delta P borne by all individual stored energy in the population E If the unbalanced power born by the mth stored energy exceeds the configured capacity limit, the mth stored energy is taken as a disturbance node and is merged into a disturbance and load node set d generated by a fault, and the unbalanced power delta P born by all the stored energy except the mth stored energy is recalculated according to the unbalanced electromagnetic power calculation node E ' unbalanced power DeltaP ' borne by non-faulty synchronous generator ' G Again determine Δ P' E Whether the configured capacity limit is exceeded or not is judged until all the unbalanced power born by the stored energy meets the configured capacity limit, and the final unbalanced power born by the non-fault synchronous generator is obtained; wherein m is more than or equal to 1 and less than or equal to N e +1，N e The energy storage number of the inertial supports on the new energy side;
and the fault switching probability calculation node is used for calculating the frequency change rate of the non-fault synchronous generator according to the final unbalanced power born by the non-fault synchronous generator, and obtaining the fault switching probability according to the limitation of the frequency change rate.
According to a third aspect of the present invention, there is provided a computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the method according to any one of the first aspects.
Generally, by the above technical solution conceived by the present invention, the following beneficial effects can be obtained:
(1) The method is based on a genetic algorithm, the total power capacity of energy storage configuration and the system fault probability of switching are calculated in the population iteration process, the optimal total power capacity of energy storage configuration and the optimal system fault probability of switching are obtained, wherein the total power capacity of energy storage configuration is used as an economic index for measuring the power system, the fault switching rate is used as a frequency stability index for measuring the power system, and an energy storage capacity configuration scheme which is most consistent with the expectation of an investor is obtained according to the economic performance of the energy storage configuration and the frequency stability weight coefficient of the system of the investor.
(2) Further, when the system fault probability switching calculation is carried out, the frequency of the power grid has a spatial distribution characteristic when the system is disturbed, a frequency response analysis model based on direct current flow is adopted, unbalanced electromagnetic power borne by energy storage and non-fault synchronous generator nodes (namely, all inertial nodes) in the power system is calculated, whether unbalanced electromagnetic power borne by all individual energy storage in a population exceeds the configured energy storage power capacity (namely, whether the unbalanced electromagnetic power is saturated or not) is judged, the nodes with the saturated energy storage power capacity are merged into the non-inertial nodes, the unbalanced power of all the inertial nodes is recalculated until all the unbalanced power borne by the energy storage meets the configured capacity limit, the final unbalanced power of all the inertial nodes can be obtained, the frequency change rate of each unit under the non-fault condition is calculated based on the unbalanced power borne by the non-fault synchronous generator, the frequency change rate reflects the frequency dynamic difference of different areas of the system, namely, the spatial distribution characteristic of the power grid frequency is considered, and meanwhile, compared with the mode of a large number of full-state time domain analysis models in the prior art, the frequency response calculation process with the energy storage power system is simplified, the frequency change rate calculation process is calculated, and the time is saved effectively.
In summary, the energy storage configuration method and system based on the power grid frequency spatial distribution characteristics of the invention includes that an outer layer is an energy storage multi-objective optimization configuration calculation model based on a genetic algorithm, and the optimal position and capacity of the configured energy storage are searched and configured aiming at the set economic and stability targets; the inner layer is a frequency response model which considers the frequency space distribution characteristic and is used for improving the inertia of the power system, and the frequency response model is used for analyzing the frequency stability of the power system with the inertia support energy storage device under the fault condition; the method can give consideration to both economy and power grid frequency stability, quickly realize the site selection and constant volume of the inertial support energy storage system, ensure the transient stability of the power grid frequency, and solve the problem that the existing energy storage configuration method is difficult to give consideration to both frequency space distribution characteristics and solving difficulty.
Drawings
FIG. 1 is a flow chart of energy storage optimization configuration considering spatial distribution characteristics of grid frequency according to the present invention;
FIG. 2 is a schematic diagram of a four-machine two-zone power system according to an embodiment of the present invention;
fig. 3 is a Pareto optimal solution set diagram of an energy storage optimization configuration result provided by the embodiment of the present invention;
fig. 4 is a comparison diagram of simulation results of different schemes provided by the embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and do not limit the invention. In addition, the technical features involved in the embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
In the present invention, the terms "first", "second", and the like in the description and the drawings are used for distinguishing similar objects, and are not necessarily used for describing a particular order or sequence.
As shown in fig. 1, the present invention provides an energy storage configuration method considering the spatial distribution characteristics of the grid frequency, which includes the following steps:
s1, according to the number N of new energy power stations of the power system e Setting the energy storage number N of the new energy side inertial supports e And upper and lower limits S of energy storage power capacity at each location eminj 、S emaxj And setting a node range set No which can be accessed by the network side centralized energy storage power station according to the traditional generator set and the new energy access node e (ii) a The method comprises the following steps of firstly, obtaining a new energy source side and a network side, wherein the new energy source side and the network side comprise upper and lower limits of energy storage power capacity;
s2, using the upper and lower limits S of the energy storage power capacity eminj 、S emaxj And an access nodePoint Range set No e For input, a genetic algorithm is adopted to randomly generate an energy storage configuration initial population X = [) 1 ,X 2 ,…,X Npop ] T (ii) a Wherein the initial population X comprises N pop Individual configuration results, N pop For a predetermined number of populations, X i Configuring an ith configuration result individual in the initial population X for energy storage, wherein the expression is as follows:
X i ＝[x i1 x i2 …x i(Ne+2) ]
wherein, when j is more than or equal to 1 and less than or equal to N e When x is ij The energy storage power capacity of the ith energy storage configuration result in the initial population X configured on the side of the j new energy station is shown, and S is satisfied eminj ≤x ij ≤S emaxj ，S emaxj Represents the upper limit of the energy storage power capacity of the new energy source side and the network side, S eminj Representing the lower limit of the energy storage power capacity of the new energy source side and the network side; when j = N e At +1, x ij The energy storage power capacity of the ith energy storage configuration result in the initial population centrally configured on the network side is shown, and S is satisfied eminj ≤x ij ≤S emaxj (ii) a When j = N e At +2, x ij The position of an energy storage node, x, of the ith energy storage configuration result in the initial population, which is centrally configured on the network side ij ∈No e ，No e And a node range set which can be accessed by the network side centralized energy storage power station.
In this embodiment, a non-inferiority ranking genetic algorithm (NSGA-II) is adopted, and in other embodiments, a genetic algorithm of other multi-objective optimization scenarios may also be adopted to generate the initial population for energy storage configuration.
S3, calculating each individual X in the randomly generated initial population X i Configured total power capacity y of energy storage 1i (ii) a N-1 fault accident set A calculated based on direct current power flow frequency response model N-1 The frequency change rate of the lower non-fault synchronous generator set is limited according to the frequency change rate to obtain the system fault switching rate y 2i Wherein i is more than or equal to 1 and less than or equal to N pop ；
Specifically, each individual X i Configured total energy storage power capacity y 1i Comprises the following steps:
wherein, the N-1 failure accident set A N-1 The method comprises the following steps: power generator tripping faults, load disconnection faults and transmission line disconnection faults.
The process of calculating the frequency change rate of the synchronous generator under the N-1 fault accident set based on the direct current power flow frequency response model comprises the following steps:
step S31, under the condition of not counting the limit of energy storage power capacity, calculating an N-1 fault accident set A N-1 Fault generating unbalanced power vector Δ P in (1) d Generating an unbalanced power vector Δ P based on the fault d Calculating unbalanced electromagnetic power delta P born by energy storage and non-fault synchronous generator nodes in electric power system GE Wherein the nodes of the energy storage and non-fault synchronous generator are inertia nodes, delta P GE Comprises the following steps:
in the formula (I), the compound is shown in the specification,
a power grid susceptance matrix which only considers power grid inductance; subscript GE represents the set of energy storage and non-fault synchronous generator nodes; subscript d represents the disturbance and load node set generated by the fault; b GE-GE 、B GE-d 、B d-GE And B d-d The block matrixes are all power grid susceptance matrixes;
is B d-d The inverse matrix of (d); delta P d An unbalanced power vector is generated for the fault.
ΔP GE ＝[ΔP E ，ΔP G ] T ，ΔP E Unbalanced electromagnetic power vector, Δ P, to be absorbed for energy storage G Unbalanced electromagnetic power assumed by non-failing synchronous generator nodesVector quantity; wherein, Δ P E ＝[ΔP E1 ,ΔP E2 ,…,ΔP ENe ,ΔP E(Ne+1) ] T Ne +1 represents N e +1 stored energy; delta P G ＝[ΔP G1 ,ΔP G2 ,…,ΔP GNg ] T And Ng represents the number of non-faulty synchronous generators.
Step S32, judging the unbalanced power delta P born by the stored energy of each individual in the population E Whether the configured energy storage power capacity x is exceeded ij (j≤N e + 1), if m is m (1. Ltoreq. M. Ltoreq.N e + 1) unbalanced power Δ P borne by the stored energy Em Exceeding the corresponding capacity limit x im Then, the mth stored energy is taken as a disturbance node and merged into a set d, wherein the set d is a disturbance and load node set generated by a fault, namely a set of non-inertial nodes; recalculating the unbalanced power delta P 'borne by the remaining energy storage (i.e. all other energy storage except the mth energy storage) and the non-faulty synchronous generator' GE Judging the unbalanced power delta P 'born by the stored energy of each individual again' E Whether the configured energy storage power capacity is exceeded until N e The unbalanced power born by +1 stored energy all meets the corresponding capacity limit to obtain N e +1 unbalanced power stored and unbalanced power born by all non-fault synchronous generators;
wherein, delta P' GE Can be expressed as:
wherein sgn (—) is a sign function; the matrix variable denoted by the band' is the matrix variable after the node where the mth energy storage is located is merged into the node set d.
Step S33, calculating the frequency change rate of the non-fault synchronous generator according to the unbalanced power born by all the non-fault synchronous generators, wherein the expression is as follows:
ng is the number of non-fault synchronous generators in the power system; f. of 1 …f Ng For each non-faulted synchronous generator rotor frequency; j. the design is a square 1 …J Ng Is the time constant of the inertia of the rotor,
representing the derivative operator.
Step S33, according to the frequency change rate limitation, obtaining the fault switching probability y 2i ：
Wherein N is AN-1 Representing the number of N-1 faults; s a Capacity of the a-th non-fault synchronous generator; e.g. of a cylinder ab Expressed as whether the frequency change rate of the a-th non-fault generator set exceeds the limit in the b-th fault scene, and when the frequency change rate of the generator set exceeds the limit, e ab And =1, otherwise 0.
S4, according to the configured total energy storage power capacity and the fault probability, performing non-dominated sorting and congestion degree sorting on the energy storage configuration population to obtain a Pareto grade and congestion degree;
s5, selecting an energy storage configuration parent by adopting a binary competitive bidding competition according to the sorted Pareto grades and the crowdedness, simulating binary intersection and polynomial variation according to a preset NSGA-II intersection probability and a variation probability, and generating energy storage configuration offspring; calculating the total energy storage capacity and the system fault switching probability of each individual configuration in the energy storage configuration filial generation, and combining the total energy storage capacity and the system fault switching probability of the energy storage configuration parent generation with the total energy storage capacity and the system fault switching probability of the filial generation to obtain the current total energy storage capacity and the system fault switching probability;
s6, iterating and executing the S4 and the S5, wherein in the iteration process of each time, the non-dominated sorting and the crowding degree sorting are carried out again, the energy storage configuration result elite is reserved according to the pareto grade and the crowding degree, a new energy storage configuration parent is screened out, and the total energy storage power capacity and the total system fault tripping rate of each individual configuration in each new energy storage configuration parent are combined with the total energy storage power capacity and the total system fault tripping rate of each individual configuration in a child generated by cross variation of the new energy storage configuration parent, so that the total energy storage power capacity and the total system fault tripping rate in the current iteration process are obtained;
circularly iterating until the preset iteration number G max Selecting a set consisting of individuals with a Pareto grade of 1 in the total energy storage power capacity and the system failure probability as a Pareto optimal solution set for output, and recording to finally obtain the optimal total energy storage power capacity y in each individual 1i [1]And the system fail-over probability y 2i [1]；
S7, carrying out optimal energy storage total power capacity y in Pareto optimal solution set 1i [1]And the system fault probability y 2i [1]Normalizing and presetting corresponding weight coefficient lambda 1 And λ 2 Obtaining an energy storage inertial support power grid optimal configuration scheme;
specifically, the optimal energy storage total power capacity index is normalized as follows:
the optimal failover probability index is normalized to:
wherein, y 1i [1]，y 2i [1]And respectively representing the total energy storage power capacity and the fault probability of the ith energy storage configuration result in the energy storage configuration Pareto optimal solution set.
In this embodiment, two normalization indexes in the Pareto optimal solution set are weighted and then the minimum value is obtained to obtain the weight coefficient λ 1 And λ 2 And further obtaining an optimal energy storage configuration scheme, wherein the expression is as follows
Wherein λ is 1 、λ 2 And respectively obtaining a weight coefficient of the set energy storage total power capacity index and a weight coefficient of the fault switching probability index to obtain an energy storage inertial support power grid optimal configuration scheme.
The invention also provides an energy storage configuration system based on the frequency space distribution characteristic of the power grid, which comprises the following components:
a parameter setting module for setting the upper and lower limits of the energy storage power capacity of the new energy source side and the network side and the node range set No accessed by the network side centralized energy storage power station e ；
An initial population generation module for collecting No with the upper and lower limits of energy storage power capacity and node range e For input, N is randomly generated by adopting an NSGA-II algorithm pop Configuring an initial population by each energy storage;
the energy storage total power capacity and fault probability switching calculation module comprises an energy storage total power capacity calculation unit and a fault probability switching calculation unit, wherein the energy storage total power capacity calculation unit is used for calculating the energy storage total power capacity configured by each individual in the initial population; the fault switching probability calculation unit is used for calculating an N-1 fault accident set A based on a direct current power flow frequency response model N-1 The frequency change rate of the lower non-fault synchronous generator is limited according to the frequency change rate to obtain the fault switching rate of the system;
the sorting module is used for sorting the energy storage configuration populations in a non-dominated sorting mode and a congestion degree sorting mode according to two indexes of the total energy storage power capacity and the fault tripping rate to obtain a Pareto grade and the congestion degree;
the Pareto optimal solution set generation module is used for selecting an energy storage configuration parent according to the Pareto grade and the crowding degree, generating energy storage configuration offspring through cross variation, respectively calculating the energy storage total power capacity and the fault probability of individuals in the energy storage configuration offspring through an energy storage total power capacity calculation unit and a fault probability calculation unit, and combining two indexes of the offspring and the parent to obtain the current energy storage total power capacity and the system fault probability; if the iteration times do not reach the preset value, the current total energy storage power capacity and the system fault probability are used as input, the step S4 is returned, otherwise, a set formed by individuals with Pareto grades of 1 in the current total energy storage power capacity and the system fault probability is used as a Pareto optimal solution set to be output;
an optimal energy storage configuration scheme generation module, configured to normalize the Pareto optimal solution set, and preset a weight coefficient λ 1 And λ 2 And obtaining an optimal energy storage configuration scheme.
Energy storage total power capacity y configured for each individual in population 1i Comprises the following steps:
wherein i is more than or equal to 1 and less than or equal to N pop When j is more than or equal to 1 and less than or equal to N e When x is ij The energy storage power capacity of the ith energy storage configuration result in the population configured on the j new energy sides is represented, and S is satisfied eminj ≤x ij ≤S emaxj (ii) a When j = N e +1 time, x ij The energy storage power capacity of the ith energy storage configuration result in the population is centrally configured on the network side and meets the requirement of S eminj ≤x ij ≤S emaxj ；S emaxj Represents the upper limit of the energy storage power capacity of the new energy source side and the network side, S eminj Representing the new energy source side and the net side energy storage power capacity lower limit.
The failure tripping rate calculation unit comprises:
the unbalanced electromagnetic power calculation node is used for calculating a direct current power flow frequency response model N-1 fault accident set A under the constraint of upper and lower limits of energy storage power capacity under the condition of not counting N-1 Fault generating unbalanced power vector Δ P in (1) d Generating an unbalanced power vector Δ P according to a fault d The relation satisfied between the electromagnetic power and the unbalanced electromagnetic power is calculated, and the unbalanced electromagnetic power delta P born by the stored energy in the power system is calculated E Unbalanced electromagnetic power Δ P borne by non-faulty synchronous generator nodes G ；
An unbalanced electromagnetic power updating node for judging the unbalanced electromagnetic power delta P born by all individual stored energy in the population E Whether or not to exceed configuredAnd energy storage power capacity, if the unbalanced power borne by the mth energy storage exceeds the configured capacity limit, merging the mth energy storage into the fault-generated disturbance and load node set d as a disturbance node, and recalculating the unbalanced power delta P 'borne by all the energy storage except the mth energy storage according to the unbalanced electromagnetic power calculation node' E And unbalanced power delta P 'borne by non-faulty synchronous generator' G Again determine Δ P' E Whether the configured capacity limit is exceeded or not is judged until all the unbalanced power born by the stored energy meets the configured capacity limit, and the final unbalanced power born by the non-fault synchronous generator is obtained; wherein m is more than or equal to 1 and less than or equal to N e +1，N e The energy storage number of the inertial supports on the new energy side;
and the fault switching probability calculation node is used for calculating the frequency change rate of the non-fault synchronous generator according to the final unbalanced power born by the non-fault synchronous generator, and obtaining the fault switching probability according to the limitation of the frequency change rate.
Example 1:
in this embodiment, a typical two-zone four-machine system is taken as an example for explanation, as shown in fig. 2, in the typical two-zone four-machine system, 1 to 11 in the figure are different bus nodes in a power grid, a synchronous generator 2 is replaced by a wind farm, inertia constants of three equivalent units are respectively 50s, 4.2s and 3.58s, and capacity of each unit is 1000MVA. And selecting to configure an energy storage system on the wind power plant side and the power grid side respectively to support the inertia of the power grid. The power capacities of the two energy storage systems and the access point of the energy storage system on the network side are variables to be optimized. The initial population size is set to 100, the number of iterations G max The setting is 1000, the cross probability is 0.9, the mutation probability is 1/3, the binary cross parameter is 10, and the mutation parameter is 10.
According to the provided optimal configuration method of the energy inertial support, a Pareto optimal solution set shown in fig. 3 can be obtained. The optimal solution set is an energy storage optimization configuration result which can be selected by an investor according to the attention degree of the investor on two indexes of the frequency stability of the power grid and the energy storage investment cost. The results of the partial configuration of the system are shown in table 1. Scheme 1 sets up rights for investors in accordance with the method of the inventionCoefficient of gravity λ 1 ＝0.5，λ 2 The energy storage configuration result is 0.5, and the capacity ratio of the generator tripping machine can be greatly reduced with smaller energy storage investment cost; scheme 2 for investors to set a weighting factor lambda according to the method of the invention 1 ＝0，λ 2 =1 energy storage configuration result, the tripping capacity ratio of the system can be minimized under the existing constraint conditions. Fig. 4 is a comparison of the survival rates of the units without energy storage configuration and with two energy storage configuration schemes. As can be seen from fig. 4, both the scheme 1 and the scheme 2 can ensure that the frequency change rate of the system equivalent unit 1 is not out of limit under the condition of N-1 fault, and because the inertia constant of the synchronous generator 3 is large, when only one energy storage power station is configured on the network side, the synchronous generator 3 is preferentially ensured not to be switched off, and the survival rate of the synchronous generator 3 is higher than that of the synchronous generator 4. Because the fault scenes are few, the survival capacity of the unit can be greatly reduced by one fault. And because only one energy storage system is configured on the network side, the configuration scheme can not ensure that the frequency change rate of all the units is within 0.5Hz/s in the transient process of the N-1 fault. If the survival rate of each unit in the system is further ensured, a proper amount of stored energy can be configured on the side of the weak-inertia synchronous generator, so that the inertia of the area is improved, and the frequency of the unit tripping under the fault condition is reduced. Table 2, table 3, and table 4 show the frequency change rate of each unit in the N-1 fault scenario after no energy storage and energy storage configuration according to the scheme 1 and the scheme 2, respectively. After energy storage is configured according to the scheme 1, the frequency change rate of the generator 1 in the fault scene 10 is reduced from 0.71Hz/s to 0.5Hz/s, so that the unit can tolerate the frequency change rate range. After energy storage is configured according to the scheme 2, the frequency change rate of the generator set 4 in the fault scene 1 is reduced from-0.57 Hz/s to-0.5 Hz/s in the scheme 2. The method can fully play the role of the energy storage inertia support and the reasonability of the energy storage inertia support configuration.
TABLE 1
TABLE 2
TABLE 3
TABLE 4
In general, based on a non-inferior ranking genetic algorithm, the energy storage optimal configuration Pareto optimal solution set which gives consideration to economy and frequency transient can be solved according to the set economic index of the total energy storage configuration power capacity and the stability index of the fault switching probability; and simultaneously, per-unit combining the two indexes of the Pareto optimal solution set based on a normalization method, and obtaining an energy storage power capacity configuration scheme which best meets the expectation of an investor according to the energy storage configuration economy of the investor and a system frequency stability weight coefficient.
According to the method, when the system is disturbed, the frequency of the power grid has a spatial distribution characteristic, the frequency response analysis model based on the direct current load flow is adopted, the frequency change rate of different units can be rapidly calculated, compared with the mode of analyzing the model through a large number of full-state time domains in the prior art, the frequency response calculation process of the power system with the energy storage is simplified, the calculation solving time of an algorithm can be effectively saved, and the frequency response analysis model considers the problems of the frequency spatial distribution characteristic and the solving difficulty.
Furthermore, those skilled in the art will appreciate that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, systems and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It will be understood by those skilled in the art that the foregoing is only an exemplary embodiment of the present invention, and is not intended to limit the invention to the particular forms disclosed, since various modifications, substitutions and improvements within the spirit and scope of the invention are possible and within the scope of the appended claims.
Claims (10)
1. An energy storage configuration method based on grid frequency spatial distribution characteristics is characterized by comprising the following steps:
s1, setting upper and lower limits of energy storage power capacity of a new energy source side and a network side and setting node range set No accessed by a network side centralized energy storage power station e ；
S2, collecting No by using upper and lower limits of energy storage power capacity and node range e Randomly generating N for input by genetic algorithm pop Configuring an initial population by each energy storage;
s3, calculating the total energy storage power capacity configured by each individual in the initial population; and calculating an N-1 fault accident set A based on a direct current power flow frequency response model N-1 The frequency change rate of the lower non-fault synchronous generator is limited according to the frequency change rate to obtain the fault switching rate of the system;
s4, according to the total energy storage power capacity and the fault tripping rate, performing non-dominated sorting and congestion degree sorting on the energy storage configuration population to obtain a Pareto grade and a congestion degree;
s5, selecting an energy storage configuration parent according to the Pareto grade and the crowding degree, generating energy storage configuration offspring through cross variation, calculating the total energy storage power capacity and the fault switching probability of individuals in the energy storage configuration offspring, and combining two indexes of the offspring and the parent to obtain the current total energy storage power capacity and the system fault switching probability; if the iteration times do not reach the preset value, the current total energy storage power capacity and the system fault probability are used as input, the step S4 is returned, otherwise, a set formed by individuals with Pareto grades of 1 in the current total energy storage power capacity and the system fault probability is used as a Pareto optimal solution set to be output;
s6, normalizing the Pareto optimal solution set and presetting a weight coefficient lambda 1 And λ 2 And obtaining an optimal energy storage configuration scheme.
2. The method of claim 1, wherein each individual in the population is matchedTotal power capacity y of the energy storage 1i Comprises the following steps:
wherein i is more than or equal to 1 and less than or equal to N pop When j is more than or equal to 1 and less than or equal to N e When x ij The energy storage power capacity of the ith energy storage configuration result in the population configured on the j new energy sources is represented, and S is satisfied eminj ≤x ij ≤S emaxj (ii) a When j = N e At +1, x ij The energy storage power capacity of the ith energy storage configuration result in the population is centrally configured on the network side and meets the requirement of S eminj ≤x ij ≤S emaxj ；S emaxj Representing the upper limit of the energy storage power capacity of the new energy source side and the network side, S eminj And representing the lower limit of the energy storage power capacity of the new energy source side and the network side.
3. The method of claim 2, wherein calculating the probability of an individual failing in the population comprises:
under the constraint of upper and lower limits of energy storage power capacity, calculating a direct current power flow frequency response model N-1 fault accident set A N-1 Fault generating unbalanced power vector Δ P in (1) d Generating an unbalanced power vector Δ P based on the fault d The relation satisfied between the electromagnetic power and the unbalanced electromagnetic power is calculated, and the unbalanced electromagnetic power delta P born by the stored energy in the power system is calculated E Unbalanced electromagnetic power Δ P borne by non-faulty synchronous generator nodes G ；
Judging unbalanced electromagnetic power delta P born by energy storage of all individuals in population E If the unbalanced power born by the mth stored energy exceeds the configured capacity limit, the mth stored energy is taken as a disturbance node and is merged into a disturbance and load node set d generated by a fault, and the unbalanced power delta P born by all the stored energy except the mth stored energy is recalculated E ' unbalanced power Δ P ' borne by non-faulted synchronous generator ' G Judging again the Δ P E ' whether the configured capacity limit is exceeded, or notAll the unbalanced power born by the stored energy meets the configured capacity limit, and the final unbalanced power born by the non-fault synchronous generator is obtained; wherein m is more than or equal to 1 and less than or equal to N e +1，N e The energy storage number of the inertial supports on the new energy side;
and calculating the frequency change rate of the non-fault synchronous generator according to the final unbalanced power born by the non-fault synchronous generator, and limiting according to the frequency change rate to obtain the fault switching probability.
4. A method according to claim 3, characterized in that the unbalanced electromagnetic power assumed by the energy storage and non-faulty synchronous generator nodes in the power system and the fault-producing unbalanced power vector Δ Ρ d Satisfies the following relation:
wherein, Δ P GE ＝[ΔP E ，ΔP G ] T ，
In order to only consider a power grid susceptance matrix of power grid inductance, subscript GE represents an energy storage and non-fault synchronous generator node set, and subscript d represents disturbance generated by faults and a load node set; b is GE-GE 、B GE-d 、B d-GE And B d-d All are block matrixes of the power grid susceptance matrix;
is B d-d The inverse matrix of (c).
5. The method of claim 4, wherein the ratio of frequency change of the non-faulty synchronous generator to the unbalanced power assumed by the non-faulty synchronous generator satisfies:
wherein, Δ P G Unbalanced electromagnetic power vector, Δ P, assumed for a non-faulty synchronous generator node G ＝[ΔP G1 ，ΔP G2 ，…，ΔP GNg ] T Ng is the number of non-fault synchronous generators in the power system; f. of 1 …f Ng For each non-faulted synchronous generator rotor frequency; j. the design is a square 1 …J Ng Is the time constant of the inertia of the rotor,
representing the derivative operator.
6. The method of claim 5, wherein the probability of fail-over y is 2i Comprises the following steps:
wherein N is AN-1 Indicates the number of N-1 faults, S a Capacity of the a-th non-failing synchronous generator, e ab Indicating whether the frequency change rate of the a-th non-fault generator exceeds the limit in the b-th fault scene, and if the frequency change rate of the generator exceeds the limit, e ab =1, otherwise, e ab ＝0。
7. The method according to any one of claims 1 to 6, wherein in step S6, the two normalized indexes of the Pareto optimal solution set are weighted and then the weighted two indexes are subjected to minimum value calculation to obtain the weight coefficient λ of the two indexes 1 And λ 2 。
8. An energy storage configuration system based on grid frequency spatial distribution characteristics, comprising:
a parameter setting module for setting the upper and lower limits of the energy storage power capacity of the new energy source side and the network side and the node range set No accessed by the network side centralized energy storage power station e ；
An initial population generation module for collecting No with the upper and lower limits of energy storage power capacity and node range e For input, an NSGA-II algorithm is adopted to randomly generate N pop Configuring an initial population by each energy storage;
the energy storage total power capacity and fault probability switching calculation module comprises an energy storage total power capacity calculation unit and a fault probability switching calculation unit, wherein the energy storage total power capacity calculation unit is used for calculating the energy storage total power capacity configured by each individual in the initial population; the fault switching probability calculation unit is used for calculating an N-1 fault accident set A based on a direct current power flow frequency response model N-1 The frequency change rate of the lower non-fault synchronous generator is limited according to the frequency change rate to obtain the fault switching rate of the system;
the sorting module is used for sorting the energy storage configuration populations in a non-dominated sorting mode and a congestion degree sorting mode according to two indexes of the total energy storage power capacity and the fault tripping rate to obtain a Pareto grade and the congestion degree;
the Pareto optimal solution set generation module is used for selecting an energy storage configuration parent according to the Pareto grade and the crowding degree, generating energy storage configuration offspring through cross variation, respectively calculating the energy storage total power capacity and the fault probability of an individual in the energy storage configuration offspring through an energy storage total power capacity calculation unit and a fault probability calculation unit, and combining two indexes of the offspring and the parent to obtain the current energy storage total power capacity and the current system fault probability; if the iteration times do not reach the preset value, the current total energy storage power capacity and the system fault probability are used as input, the step S4 is returned, otherwise, a set formed by individuals with Pareto grades of 1 in the current total energy storage power capacity and the system fault probability is used as a Pareto optimal solution set to be output;
an optimal energy storage configuration scheme generation module, configured to normalize the Pareto optimal solution set, and preset a weight coefficient λ 1 And λ 2 And obtaining an optimal energy storage configuration scheme.
9. The system of claim 8, wherein the fault tripping rate calculation unit comprises:
the unbalanced electromagnetic power calculation node is used for calculating a direct current power flow frequency response model N-1 fault accident set A under the constraint of upper and lower limits of energy storage power capacity under the disregard of N-1 Fault generating unbalanced power vector Δ P in (1) d Generating an unbalanced power vector Δ P based on the fault d The relation satisfied between the electromagnetic power and the unbalanced electromagnetic power is calculated, and the unbalanced electromagnetic power delta P born by the stored energy in the power system is calculated E Unbalanced electromagnetic power delta P borne by non-faulted synchronous generator nodes G ；
An unbalanced electromagnetic power updating node for judging the unbalanced electromagnetic power delta P born by all individual stored energy in the population E If the unbalanced power borne by the mth stored energy exceeds the configured capacity limit, merging the mth stored energy into the fault-generated disturbance and load node set d as a disturbance node, and recalculating the unbalanced power delta P 'borne by all the stored energy except the mth stored energy according to the unbalanced electromagnetic power calculation node' E And unbalanced power delta P 'borne by non-faulted synchronous generator' G Again determine Δ P' E Whether the configured capacity limit is exceeded or not is judged until all the unbalanced power born by the stored energy meets the configured capacity limit, and the final unbalanced power born by the non-fault synchronous generator is obtained; wherein m is more than or equal to 1 and less than or equal to N e +1，N e The energy storage number of the inertial supports on the new energy side;
and the fault switching probability calculation node is used for calculating the frequency change rate of the non-fault synchronous generator according to the final unbalanced power born by the non-fault synchronous generator, and obtaining the fault switching probability according to the limitation of the frequency change rate.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211054711.8A CN115276055B (en) | 2022-08-31 | 2022-08-31 | Energy storage configuration method and system based on power grid frequency spatial distribution characteristics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211054711.8A CN115276055B (en) | 2022-08-31 | 2022-08-31 | Energy storage configuration method and system based on power grid frequency spatial distribution characteristics |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115276055A true CN115276055A (en) | 2022-11-01 |
| CN115276055B CN115276055B (en) | 2023-04-18 |
Family
ID=83755031
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211054711.8A Active CN115276055B (en) | 2022-08-31 | 2022-08-31 | Energy storage configuration method and system based on power grid frequency spatial distribution characteristics |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115276055B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160099567A1 (en) * | 2014-10-02 | 2016-04-07 | Mitsubishi Electric Research Laboratories, Inc. | Dynamic and Adaptive Configurable Power Distribution System |
| CN109347127A (en) * | 2018-10-27 | 2019-02-15 | 华中科技大学 | A kind of energy storage Optimal Configuration Method for coping with block supply line fault |
| CN110417012A (en) * | 2019-08-01 | 2019-11-05 | 国电南瑞科技股份有限公司 | Power grid energy accumulation capacity configuration and system under a kind of constraint of frequency security |
| CN111146820A (en) * | 2019-12-31 | 2020-05-12 | 国网浙江省电力有限公司嘉兴供电公司 | Hybrid energy storage system optimal configuration method considering wind power uncertainty |
| CN112906200A (en) * | 2021-01-25 | 2021-06-04 | 国网山西省电力公司电力科学研究院 | Power system energy storage configuration method and device, computer equipment and storage medium |
-
2022
- 2022-08-31 CN CN202211054711.8A patent/CN115276055B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160099567A1 (en) * | 2014-10-02 | 2016-04-07 | Mitsubishi Electric Research Laboratories, Inc. | Dynamic and Adaptive Configurable Power Distribution System |
| CN109347127A (en) * | 2018-10-27 | 2019-02-15 | 华中科技大学 | A kind of energy storage Optimal Configuration Method for coping with block supply line fault |
| CN110417012A (en) * | 2019-08-01 | 2019-11-05 | 国电南瑞科技股份有限公司 | Power grid energy accumulation capacity configuration and system under a kind of constraint of frequency security |
| CN111146820A (en) * | 2019-12-31 | 2020-05-12 | 国网浙江省电力有限公司嘉兴供电公司 | Hybrid energy storage system optimal configuration method considering wind power uncertainty |
| CN112906200A (en) * | 2021-01-25 | 2021-06-04 | 国网山西省电力公司电力科学研究院 | Power system energy storage configuration method and device, computer equipment and storage medium |
Non-Patent Citations (3)
| Title |
|---|
| XIA CHEN: "A Model Predictive Control Method for Hybrid", 《CSEE JOURNAL OF POWER AND ENERGY SYSTEMS》 * |
| 朱俊澎: "考虑输电网潮流约束的时序生产模拟及新能源消纳能力评估", 《电网技术》 * |
| 谢石骁等: "基于机会约束规划的混合储能优化配置方法", 《电网技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115276055B (en) | 2023-04-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Duman et al. | Optimal power flow of power systems with controllable wind‐photovoltaic energy systems via differential evolutionary particle swarm optimization | |
| Muthuswamy et al. | Environmental and economic power dispatch of thermal generators using modified NSGA‐II algorithm | |
| CN108233359B (en) | Partition optimization method and system for extra-high voltage alternating current-direct current feed-in receiving end power grid | |
| CN104835018B (en) | The reliability estimation method and device of medium voltage distribution network rack | |
| CN107994582B (en) | Method and system for reconstructing power distribution network containing distributed power supply | |
| CN110689186B (en) | Power distribution network power supply reliability assessment method and system based on wind power output randomness | |
| CN104866919A (en) | Multi-target planning method for power grid of wind farms based on improved NSGA-II | |
| Li et al. | Optimization of loop-network reconfiguration strategies to eliminate transmission line overloads in power system restoration process with wind power integration | |
| Ali et al. | Optimal planning of RDGs in electrical distribution networks using hybrid SAPSO algorithm | |
| CN115207935B (en) | Reactive power coordination optimization method for improving transient voltage stability of voltage weak area | |
| Rupolo et al. | A new parallel and decomposition approach to solve the medium-and low-voltage planning of large-scale power distribution systems | |
| CN109888817B (en) | Method for carrying out position deployment and capacity planning on photovoltaic power station and data center | |
| Rezaeian‐Marjani et al. | Probabilistic assessment of D‐STATCOM operation considering correlated uncertain variables | |
| Wei et al. | Transmission network planning with N-1 security criterion based on improved multi-objective genetic algorithm | |
| Wang et al. | Pareto optimization of power system reconstruction using NSGA-II algorithm | |
| CN105069517A (en) | Power distribution network multi-objective fault recovery method based on hybrid algorithm | |
| CN115276055B (en) | Energy storage configuration method and system based on power grid frequency spatial distribution characteristics | |
| CN110751328A (en) | High-proportion renewable energy power grid adaptive planning method based on joint weighted entropy | |
| Hosseini et al. | Optimal placement of dgs in distribution system including different load models for loss reduction using genetic algorithm | |
| CN113629769B (en) | Line weight-based power grid partition searching method and system | |
| CN115764859A (en) | Fault recovery method, system, equipment and medium for multi-energy coupling power distribution network | |
| Shuaishuai et al. | Optimal testing resource allocation for modular software systems based-on multi-objective evolutionary algorithms with effective local search strategy | |
| Ali et al. | Distributed generation sizing and placement using computational intelligence | |
| Sattianadan et al. | Power loss minimization by the placement of DG in distribution system using PSO | |
| CN114204613A (en) | Reactive compensation method and system for offshore wind farm access power system |
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 |