CN116205337A - Annual operation evaluation method and device for source network charge storage project planning - Google Patents
Annual operation evaluation method and device for source network charge storage project planning Download PDFInfo
- Publication number
- CN116205337A CN116205337A CN202211714372.1A CN202211714372A CN116205337A CN 116205337 A CN116205337 A CN 116205337A CN 202211714372 A CN202211714372 A CN 202211714372A CN 116205337 A CN116205337 A CN 116205337A
- Authority
- CN
- China
- Prior art keywords
- annual
- power
- source network
- new energy
- sum
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims description 20
- 239000002245 particle Substances 0.000 claims abstract description 92
- 238000004146 energy storage Methods 0.000 claims abstract description 53
- 230000005611 electricity Effects 0.000 claims abstract description 50
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 37
- 230000006870 function Effects 0.000 claims description 37
- 238000004364 calculation method Methods 0.000 claims description 27
- 238000010248 power generation Methods 0.000 claims description 20
- 239000002699 waste material Substances 0.000 claims description 14
- 238000007599 discharging Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 28
- 238000012423 maintenance Methods 0.000 abstract description 5
- 238000004088 simulation Methods 0.000 abstract description 5
- 230000002349 favourable effect Effects 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 235000013305 food Nutrition 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000001580 bacterial effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003592 biomimetic effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002431 foraging effect Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002922 simulated annealing 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
-
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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/381—Dispersed generators
-
- 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]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- General Health & Medical Sciences (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Educational Administration (AREA)
- Power Engineering (AREA)
- Game Theory and Decision Science (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Primary Health Care (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Water Supply & Treatment (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a method and a device for evaluating annual operation of source network charge storage project planning, which consider the charge and discharge plan of an energy storage device in each period in the operation process of the source network charge storage project according to the characteristics of combined operation of wind power, photovoltaic, energy storage, a power grid and loads. The method is characterized by taking the total annual electricity purchasing cost, the energy storage device loss cost and the minimum equipment operation and maintenance cost as targets, taking into consideration the operation constraints of wind power, photovoltaic, a power grid and energy storage equipment, the charging rate constraint of energy storage and the like, and providing a source network charging and storage annual 8760-hour operation simulation method based on a particle swarm algorithm. And finally, the method can output the economic index of the source network in any period and throughout the year, the technical indexes such as the new energy electricity discarding rate, the new energy electricity ratio and the like, and provides a favorable support for decision.
Description
Technical Field
The invention belongs to the technical field of power supply system simulation, and particularly relates to a method and a device for evaluating annual operation of source network load storage project planning.
Background
The source network charge storage system is a novel network structure consisting of a distributed power supply, energy storage equipment, electric equipment, a main power grid and a control device. The source network charge storage system is an important means for solving a plurality of problems of the power system, and needs to exert huge economic benefits brought by the source network charge storage system as much as possible. Thus, the economical operation of the source network charge storage system is a key to attract users and enable it to be generalized in the power system. However, as the types of devices in the source network charge storage system increase and the scale is continuously enlarged, the problem of optimizing the economic operation of the micro-grid becomes more and more complex.
The source network load storage integrated project comprises a power supply such as wind power, photovoltaic and the like, a power grid, a load and an energy storage system; the running targets comprise low cost, high reliability, low waste wind and light rate and the like of the whole system life cycle, and the constraint conditions comprise total electric quantity balance, real-time electric power balance, energy storage device power constraint, energy storage device electric quantity constraint, new energy electric quantity duty ratio constraint, adjustable characteristic constraint of load and the like. The source network charge storage economic operation problem is a complex nonlinear mixed integer random programming problem, and has the characteristics of high variable dimension, large calculation scale, continuous variable discrete variable inclusion, multiple constraint conditions, nonlinear objective function and the like. Based on the difficulty of optimizing the problem, the traditional solution is easy to fall into dimension disasters in a limited scheduling period, so that the problem can be solved theoretically, which takes years or cannot be solved at all.
Intelligent optimization algorithms are increasingly accepted as compared to traditional planning algorithms. Intelligent computation is also called soft computation, which is a solution algorithm designed according to the principle of the intelligent computation and imitating the law of the intelligent computation after people are inspired by the development law of the nature or the biology. In recent years, many biomimetic intelligent optimization algorithms, which are quite different from the principles of classical mathematical programming and which are expected to solve complex optimization problems by simulating natural ecosystems, are proposed one after the other, like simulated annealing algorithms, artificial immune algorithms, bacterial foraging algorithms and particle swarm algorithms, genetic algorithms, etc. These algorithms are each advantageous and disadvantageous and are widely used by people with unique flow characteristics and search characteristics.
Among them, the particle swarm algorithm (Particle Swarm Optimization, PSO) is a subject that has emerged in the 90 th century, and is known for its simple concept, convenient implementation, and fast convergence. The basic idea of the particle swarm algorithm is to simulate the predation behavior of a bird swarm randomly searching for food, and the bird swarm adjusts its search path through own experience and communication between the swarms, so as to find the place with the largest food. Where the position/path of each bird is then the combination of arguments, the food density at each arrival location, i.e. the function value. Each search adjusts its own search direction and speed according to its own experience (the optimal place searched by its own history) and population communication (the optimal place searched by its population history), which is called tracking extremum, so as to find the optimal solution.
Referring to fig. 1, the basic steps of the particle swarm algorithm are:
step1, determining the motion state of a particle is described by two parameters of position and speed, and thus the two parameters are initialized;
step2, the result (function value) of each search is particle fitness, and then the individual historical optimal position of each particle and the historical optimal position of the group are recorded;
step3, the historical optimal position of the individual and the historical optimal position of the group are equivalent to generating two forces, and the motion state of the particles is influenced by the combination of the inertia of the particles, so that the position and the speed of the particles are updated.
The particle swarm algorithm has the following characteristics:
(1) compared with the traditional algorithm, the method has the advantages that the calculation speed is very high, and the global searching capability is also very high;
(2) PSO is not very sensitive to the size of the population, so the initial population is set to 500-1000, and the influence of speed is not great;
(3) the particle swarm algorithm is suitable for the problem of continuous function extremum, and has stronger global searching capability for the problem of nonlinearity and multimodal.
The method is very suitable for solving the problems of high variable dimension, large calculation scale, continuous variable discrete variable inclusion, various constraint conditions and nonlinear objective function of the source network charge storage optimization operation based on the characteristics of high calculation speed of a particle swarm algorithm and strong global searching capability for nonlinear and multimodal problems.
Most of the current methods for evaluating the planning schemes adopt 4 typical day methods, which are respectively calculated and verified in spring, summer, autumn and winter in one day, and the typical day methods have higher speed, but because of fewer included scenes, the difference between the calculated result and the real situation is larger, so that investment waste is easily caused, or the running effect cannot be expected.
Disclosure of Invention
The invention provides a method and a device for evaluating annual operation of source network charge storage project planning, so that the planned values of running cost, new energy electric quantity ratio, new energy electricity rejection rate and the like in all periods of the year corresponding to the planned wind power, photovoltaic and energy storage capacity can be counted, and the source network charge storage project planning scheme can be evaluated according to the running cost, the new energy electric quantity ratio and the new energy electricity rejection rate, and the evaluation result is more accurate.
In order to achieve the purpose, the annual operation evaluation method for planning the source network load storage project comprises the following steps:
step1, constructing annual hourly power generation capacity of wind power and photovoltaic according to capacity configuration of wind power and photovoltaic and annual hourly wind power and photovoltaic power output determined by a project planning scheme; according to the load information, constructing the annual hour-by-hour load power consumption;
step2, constructing constraint conditions according to a real-time balance principle of electric power;
step3, constructing an objective function;
and step 7, evaluating the planning scheme of the source network charge storage project according to the technical indexes of the new energy.
Further, in step2, the constraint conditions include electric power and power balance constraint, operation constraint of wind power equipment, operation constraint of photovoltaic equipment, operation constraint of grid equipment and charging rate constraint of energy storage.
Further, in step3, the objective function is the lowest total annual cost.
Further, in step 4, the number of groups is 200, and the number of iterations is 10 to 15.
Further, step 5 includes the steps of:
step 5.1, initializing all particles, wherein each particle is an array formed by the hours of the current month, and each number represents the charge and discharge plan of the energy storage device at the current moment;
step 5.2, calculating an objective function for all particles to obtain objective function values of different particles;
step 5.3, taking the objective function value with the smallest objective function value as the population optimal function value of the round, and taking the particle corresponding to the optimal function value as the optimal particle of the round;
step 5.4, generating N particles of the next round by adopting the following formula:
v f (t+1)=v f (t)+c 1 ×rand×(pbest f -SOC f )+c 2 ×rand×(gbest-SOC f )
SOC f (t+1)=SOC f +v f (t+1)
wherein v is f (t+1) represents the speed of the f-th particle at the time of updating, v f (t) represents the speed of the f-th particle at this time, c 1 Represents a self-learning factor, rand is a random number between 0 and 1, pbest f For the f-th particle self-historic optimal position, SOC f An array of the present time for the f-th particle, c 2 Representing a population learning factor, wherein gbest is the historical optimal particle position of the population;
step 5.5, judging the position of each particle, and when the position or the speed of the particle exceeds a limit value, pulling the particle back to the boundary;
step 5.6, repeating the steps 5.2 to 5.5 until the maximum iteration times are reached;
and 5.7, outputting the final population optimal particle position as a result of a particle swarm algorithm, namely an optimal charging and discharging plan of the energy storage device.
Further, in step 6, the electricity purchasing cost of the power grid is as follows: the sum of the product of the main grid power and its electricity price per hour throughout the year.
Further, in step 6, the new energy power consumption P is calculated by hour-by-hour calculation or sum calculation abd,t ,
The calculation formula of the hour-by-hour algorithm is as follows:
P abd,t =P WT,t +P PV,t +P BESS,t -P Load,t
wherein P is WT,t For the t hour wind power, P PV,t For the t hour photovoltaic power, P BESS,t Power of energy storage device at t hour, P Load,t The power used for loading at the t hour;
the sum calculation method is as follows:
P abd,sum =P WT,sum +P PV,sum +P Grid,sum -P Load,sum
wherein P is abd,sum Represents the new energy waste amount of the whole year, P WT,sum Represents total annual wind power generation amount, P PV,sum Representing total photovoltaic power generation amount, P Grid,sum Representing total power purchase amount from power grid all year around, P Load,sum Representing the total annual load power consumption.
Further, in step 6, the new energy power rejection rate γ abd The calculation method of (1) is as follows:
wherein P is abd,sum Represents the new energy waste amount of the whole year, P WT,sum Represents total annual wind power generation amount, P PV,sum Representing the total amount of photovoltaic power generation throughout the year.
Further, in step 6, the new energy electric quantity is the ratio gamma green The calculation method comprises the following steps:
in the above, P abd,sum Represents the new energy waste amount of the whole year, P WT, Represents total annual wind power generation amount, P PV, Representing total photovoltaic power generation amount, P Load, Representing the total annual load power consumption.
An annual operation evaluation device for planning a source network charge storage project comprises:
the input module is used for acquiring basic data, wherein the basic data comprises capacity configuration of wind power and photovoltaic, and wind power, photovoltaic output and load information of each year and every hour;
the processing module is used for storing an objective function and constraint conditions, and solving the objective function and constraint conditions according to a particle swarm algorithm to obtain an optimal charge-discharge plan of the energy storage device;
and the output module is used for calculating new energy technical indexes according to the basic data and the optimal charging and discharging plan of the energy storage device, wherein the new energy technical indexes comprise the power grid electricity purchasing cost, the new energy electricity discarding quantity, the new energy electricity discarding rate and the new energy electricity occupying ratio.
Compared with the prior art, the invention has at least the following beneficial technical effects:
according to the characteristics of wind power, photovoltaic, energy storage, power grid and load combined operation, according to the capacity configuration of wind power and photovoltaic determined by a source network charge storage project planning scheme, the charge and discharge plan of an energy storage device in each period of the whole year is considered in the source network charge storage project simulation operation process, new energy technical indexes of the source network charge storage project in the whole year under the charge and discharge plan are calculated, the new energy technical indexes comprise a scene of 365 days of the whole year, the accuracy and the reliability of the new energy technical indexes are obviously higher than those of a method of calculating 4 typical days of data only in 4 days, and the source network charge storage project planning scheme can be effectively evaluated, so that investment is guided, and waste or unexpected operation effects are avoided.
The method is characterized in that the aim of lowest annual operation cost is achieved, real-time balance of electric power, the upper and lower power limits of all equipment, the upper and lower charge rate limits of energy storage and the equal charge rate of energy storage at the initial and end moments are taken as constraint conditions, the full-feasible-area optimization solution based on a particle swarm algorithm is completed, and finally the economic indexes such as annual power grid electricity purchasing cost and the technical indexes such as new energy waste, new energy waste rate and new energy electricity ratio are obtained. Has the following advantages:
first: the ultra-large scale (8760 hours) data can be calculated, and the calculation method section by section is adopted, so that the calculation difficulty is greatly reduced, the calculation speed is improved, and the 365-day source network load storage operation evaluation based on the particle swarm algorithm is realized for the first time under the condition that the result accuracy is not affected.
Second,: in the calculation process and result, the charge and discharge cost of the energy storage device and the power grid electricity purchasing cost are considered, and in the calculation of the power grid electricity purchasing cost, the peak, valley and flat electricity price are fully considered. The cost of calculation is accurate, and the method provides favorable support for decision.
Third,: the electric quantity, the electric quantity and the green electric duty ratio of any time and any time period can be calculated. And a quantitative basis is provided for whether the planning scheme of the project can meet the rigid requirements of the policy.
Drawings
FIG. 1 is a flow chart of a particle swarm algorithm;
FIG. 2 is an hour-by-hour output curve of a month of January photovoltaic;
FIG. 3 is an hour-by-hour output curve of the January wind power;
FIG. 4 is a strategy block diagram of an annual hourly operation evaluation method for source network load storage project planning;
FIG. 5 is an example of a velocity generation process for a particle swarm algorithm;
FIG. 6 is a graph of power grid electricity purchases from year to year;
FIG. 7 is a new energy power ratio year by year;
FIG. 8 shows new energy power rejection rate month by month throughout the year;
fig. 9 is a schematic diagram of an annual 8760 hour operation evaluation device for a source network load storage project.
Detailed Description
In order to make the purpose and technical scheme of the invention clearer and easier to understand. The present invention will now be described in further detail with reference to the drawings and examples, which are given for the purpose of illustration only and are not intended to limit the invention thereto.
In the description of the present invention, it should be understood that the terms "center", "longitudinal", "lateral", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", etc. indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, are merely for convenience in describing the present invention and simplifying the description, and do not indicate or imply that the devices or elements referred to must have a specific orientation, be configured and operated in a specific orientation, and thus should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include one or more such feature. In the description of the present invention, unless otherwise indicated, the meaning of "a plurality" is two or more. In the description of the present invention, it should be noted that, unless explicitly specified and limited otherwise, the terms "mounted," "connected," and "connected" are to be construed broadly, and may be either fixedly connected, detachably connected, or integrally connected, for example; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the above terms in the present invention will be understood in specific cases by those of ordinary skill in the art.
The invention provides a rapid and efficient annual hourly (8760 hours) source network load storage operation simulation method which can cope with various complex constraint conditions, so that the planned annual operation cost, new energy electric quantity ratio, new energy electricity rejection rate and other numerical values corresponding to wind power, photovoltaic and energy storage capacity can be counted. Providing quantized data support for decision makers.
Example 1
The annual operation evaluation method for planning the source network charge storage project comprises the following steps:
and step1, constructing an annual hourly generating capacity curve of the wind power and the photovoltaic according to the capacity configuration of the wind power and the photovoltaic and the annual hourly wind power and photovoltaic generating capacity determined by the project planning scheme. And constructing a load power consumption curve hour by hour all year according to the load information. Fig. 2 is an hour-by-hour power generation curve of the photovoltaic in 1 month, and fig. 3 is an hour-by-hour power generation curve of the wind power in 1 month.
Step2, constructing constraint conditions
According to the real-time balance principle of the power, constructing power balance constraint:
P WT,t +P PV,t +P Grid,t +P BESS,t =P Load,t
in the above, P WT,t For the t hour wind power, P PV,t For the t hour photovoltaic power, P Grid,t Main power grid power, P for the t hour BESS,t The power of the energy storage device at the t hour, wherein the power of the energy storage device is either positive or negative, representing power supply to the load, or represents charging. P (P) Load,t And (5) using electric power at t hours as a load.
According to the characteristics of wind power, photovoltaic, power grid and energy storage, constructing operation constraint of each device:
P WT,min ≤P WT,t ≤P WT,max
P PV,min ≤P PV,t ≤P PV,max
P Grid,min ≤P Grid,t ≤P Grid,max
P BESS,min ≤P BESS,t ≤P BESS,max
in the above, P WT,max Representing the upper power output limit of wind power, P WT,min Representing wind powerLower power output limit, P PV,max Representing the upper power output limit of photovoltaic, P PV,min Representing the lower power output limit, P, of the photovoltaic Grid,max Representing the upper power output limit, P, of a large power grid Grid,min The lower limit of the power output of the wind power is represented, and the source network charge storage does not allow power to be transmitted to the power grid according to the policy requirement, so P Grid,min ≥0,P BESS,max Represents the upper power output limit of stored energy, P BESS,min Representing the lower power output limit of the stored energy.
According to the energy storage operation characteristics, constructing the charge rate (SOC) constraint of the energy storage:
SOC BESS,min ≤SOC BESS,t ≤SOC BESS,max
in the above, SOC BESS,t Representing the charge rate (percent of remaining charge) of the energy storage device at t hours, SOC BESS,max Representing the upper limit of the charge rate of the energy storage device, SOC BESS,min Representing the lower limit of the charge rate of the energy storage device.
SOC BESS,start =SOC BESS,End
In the above, SOC BESS,start Representing the charge rate and SOC of the energy storage device at the beginning of the scheduling period BESS,End Representing the charge rate of the energy storage device at the end of the scheduling period.
Step3, constructing an objective function
The operation and maintenance cost, loss cost, electricity purchasing cost from a main power grid and the like of the functional equipment are fully considered, and the economic cost is taken as an objective function, as follows:
in the above, minF 1 F is the objective function value Op,t Representing the operation and maintenance cost at the t hour, F Loss,t Represents the electric energy loss cost of the energy storage device at the t hour, F Grid,t Representing the cost of purchasing electricity from the main grid at hour t.
F Op, = WT ×C WT + PV ×C PV + BESS ×C BESS
In the above, gamma WT Representing the operation and maintenance coefficient (unit is yuan/kW) of wind power, C WT Representing the installed capacity of wind power, gamma PV Representing the operation and maintenance coefficient (unit is yuan/kW) of the photovoltaic, C PV Representing the installed capacity of the photovoltaic, S BESS Operation and maintenance coefficient (unit is yuan/kW) representing energy storage, C BESS Representing the installed capacity of the stored energy.
F Loss, = ch, ×E ch, ×Cost t + dis,t ×E dis,t ×Cost t
In the above, eta ch, Represents the efficiency of energy storage and charging at t hours, E ch, Represents the energy of energy storage and charge at t hour, eta dis,t Represents the efficiency of energy storage and discharge at t hours, E dis,t Represents the energy of the energy storage discharge at the t hour, cost t Representing the electricity costs at hour t.
F Grid, =rice Buy,t ×E Grid,
In the above formula, price Buy,t Representing the price of electricity purchased from the power grid at the t hour E Grid, Representing the energy purchased from the grid at hour t.
The population number N of the particle swarm, the iteration number Iter of the particle swarm and the self-learning factor c of the particle swarm are required to be set 1 Population learning factor c of particle swarm 2 Upper speed limit v of particle swarm max Lower speed limit v of particle swarm min . The greater the population number and the number of iterations, the greater the accuracy of the final solution, but the longer the computational time spent. The smaller the population number and the number of iterations, the lower the accuracy of the final solution, but the shorter the calculation time taken. The population number proposal is selected as 200, and the iteration number proposal is selected as 10-15.
In order to realize operation acceleration, the annual time is equally divided into 12 sections within 8760 hours, each section is calculated by adopting the following method, and finally the statistics are summed up.
1) Initializing all particles, wherein each particle is an array formed by the hours of the current month, each number represents the charge and discharge plan of the energy storage device at the current moment, and all particles are uniformly and randomly distributed in a solution space in the initialization process, and the solution space is a multidimensional space meeting constraint conditions.
2) Calculating an objective function for all particles to obtain a function value;
3) And selecting the optimal function value as the population optimal function value of the round, and selecting the particle corresponding to the optimal function value as the optimal particle of the round.
4) The following formula is used to generate the next round of N particles:
v f (t+1)=v f (t)+c 1 ×rand×(pbest f -SOC f )+c 2 ×rand×(gbest-SOC f )
SOC f (t+1)=SOC f +v f (t+1)
in the above, v f (t+1) represents the speed of the f-th particle at the time of updating, v f (t) represents the speed of the f-th particle at this time, c 1 Represents a self-learning factor, rand is a random number between 0 and 1, pbest f For the f-th particle self-historic optimal position, SOC f An array of the present time for the f-th particle, c 2 Representing a population learning factor, gbest is the historical optimal particle position of the population.
By the above formula, the velocity of each particle can be calculated. SOC (State of Charge) f (t+1) represents the updated array of the f-th particle. For a fixed value point in the initialization process, the value is maintained unchanged. The concept is shown in fig. 4. For each iteration, the optimal value pbest of each particle and the optimal value gbest of the population are updated.
5) Determining the position of each particle
When the position or velocity of the particle exceeds a limit value, it is necessary to pull it back to the boundary.
(1) When the particle velocity v f (t) > particle groupUpper speed limit v of (2) max When let v f (t)=v max ;
When the particle velocity v f (t) < upper speed limit v of particle swarm min When let v f (t)=v min ;
(2) When SOC is f >SOC max At the time, let SOC f =SOC max ;
When SOC is f <SOCmin, let SOC f =SOCmin。
Wherein SOC is max The SOCmin is the minimum value of the charge rate of the stored energy; the programming method adopted is as follows:
6) Repeating the iteration until reaching the maximum iteration number
Repeating 2), 3), 4), 5) until the maximum number of iterations Iter is reached.
7) And outputting a final population optimal function value gbest which is an array and represents the optimal charge and discharge plan of the energy storage device as a result of a particle swarm algorithm.
Results P calculated according to particle swarm algorithm BESS,t The formula for real-time power balance can be combined:
P WT,t +P PV,t +P Grid,t +P BESS,t =P Load,t
at this time P is known WT,t 、P PV,t 、P Load,t 、P BESS,t The power P of the main power grid at the t hour can be obtained Grid,t 。
And Price Buy,t Is a price that takes into account different periods of peak-to-valley, for example:
price will be Price Buy,t And P Grid,t Multiplying can result in the cost of purchasing electricity from the main grid. The total cost of the annual green electricity can be calculated by flattening the electricity cost. The total cost of annual power supply is the total cost of annual power supply, namely the cost of annual power purchase from a main power grid plus the cost of green power.
According to the charge and discharge quantity of the energy storage device in each period and the charge and discharge efficiency of the energy storage device, the charge and discharge loss of the energy storage device can be calculated, the charge and discharge loss of the energy storage device in each period is multiplied by the electricity consumption cost, and then the energy storage device and the electricity consumption cost are summed together, so that the annual energy storage device loss cost can be obtained.
The operation and maintenance cost of the equipment is generally a fixed proportion multiplied by the investment of the equipment, and wind power, photovoltaic and energy storage equipment respectively have different operation and maintenance cost coefficients.
The annual total cost can be obtained by summing the annual electricity purchasing cost, the energy storage device loss cost and the equipment operation and maintenance cost.
Because permit to sell the electricity to the electric wire netting, consequently when new energy generated energy is greater than load power consumption, and energy storage device also can not absorb, then can cause new energy to abandon the electricity, calculate new energy and abandon the electric quantity all year round and have 2 modes:
the first is to calculate the following formula on an hour-by-hour basis:
P abd,t =P WT,t +P PV,t +P BESS,t -P Load,t
in the above, P abd,t Represents tThe new energy waste amount of the hour. For all P's greater than 0 abd,t And adding to obtain the total electric quantity abandoned all the year round.
The second method is a sum-total calculation method,
P abd,sum = WT, + PV, + Grid, -Load,
in the above, P abd,sum Representing the new energy waste amount of the whole year. P (P) WT, Represents total annual wind power generation amount, P PV, Representing total photovoltaic power generation amount, P Grid, Representing total power purchase amount from power grid all year around, P Load, Representing the total annual load power consumption.
Calculating new energy electricity rejection rate:
in the above, gamma abd The power rejection rate is used as new energy.
The annual green duty cycle of the project was calculated:
in the above, gamma green Is annualGreen electrical duty cycle.
The capacity configuration of wind power, photovoltaic and energy storage in the calculation example is as follows:
the photovoltaic power generation amount of 1-12 months is as follows:
the wind power generation amount of 1-12 months is as follows:
referring to fig. 6, the main grid electricity purchases for 1-12 months are as follows:
month of month | Generating capacity (Yi degree) | Month of month | Generating capacity (Yi degree) |
1 | 5.88 | 7 | 5.05 |
2 | 6.07 | 8 | 4.98 |
3 | 5.91 | 9 | 4.92 |
4 | 3.95 | 10 | 5.09 |
5 | 4.02 | 11 | 6.02 |
6 | 4.97 | 12 | 7.10 |
Referring to fig. 7, the new energy power ratio is 1-12 months
Referring to fig. 8, the new energy electricity rejection rate is 1-12 months
And 8, evaluating the planning scheme of the source network charge storage project according to the technical index of the new energy.
Example 2
Referring to fig. 9, an annual 8760 hour operation evaluation device for a source network load storage item is characterized by comprising:
the input module is used for acquiring basic data, wherein the basic data comprises capacity configuration of wind power and photovoltaic, and wind power, photovoltaic output and load information of each year and every hour;
the calculation module is used for storing an objective function and constraint conditions, and solving the objective function and constraint conditions according to a particle swarm algorithm to obtain an optimal charge-discharge plan of the energy storage device;
and the output module is used for calculating new energy technical indexes according to the basic data and the optimal charging and discharging plan of the energy storage device, wherein the new energy technical indexes comprise the power grid electricity purchasing cost, the new energy electricity discarding quantity, the new energy electricity discarding rate and the new energy electricity occupying ratio.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a non-volatile computer readable storage medium executable by a processor. Based on this understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
According to the characteristics of wind power, photovoltaic, energy storage, power grid and load combined operation, the charge and discharge plan of the energy storage device in each period is considered in the operation process of a source network charge storage project. The annual electricity purchasing cost, the energy storage device loss cost and the equipment operation and maintenance cost are taken as targets, and the annual 8760-hour operation simulation method of the source network charge storage based on the particle swarm algorithm is provided by taking the operation constraints of wind power, photovoltaic, a power grid and energy storage equipment, the charge rate (SOC) constraint of energy storage and the like into consideration. And finally, the method can output the technical indexes such as any period of the charge storage of the source network, the annual economic index, the new energy electricity discarding rate, the new energy electricity ratio and the like.
The above is only for illustrating the technical idea of the present invention, and the protection scope of the present invention is not limited by this, and any modification made on the basis of the technical scheme according to the technical idea of the present invention falls within the protection scope of the claims of the present invention.
Claims (10)
1. The annual operation evaluation method for planning the source network charge storage project is characterized by comprising the following steps of:
step1, constructing annual hourly power generation capacity of wind power and photovoltaic according to capacity configuration of wind power and photovoltaic and annual hourly wind power and photovoltaic output determined by a source network charge storage project planning scheme; according to the load information, constructing the annual hour-by-hour load power consumption;
step2, constructing constraint conditions according to a real-time balance principle of electric power;
step3, constructing an objective function;
step 4, setting the group number and the iteration number of the particle swarm algorithm;
step 5, solving an equation set consisting of constraint conditions and objective functions by using a particle swarm algorithm to obtain an optimal charge-discharge plan of the energy storage device;
step 6, calculating new energy technical indexes including power grid electricity purchasing cost, new energy electricity discarding amount, new energy electricity discarding rate and new energy electricity occupying ratio according to the annual hourly generated energy of wind power and photovoltaic, annual hourly load power consumption and the optimal charging and discharging plan of the energy storage device;
and step 7, evaluating the planning scheme of the source network charge storage project according to the technical indexes of the new energy.
2. The annual operation evaluation method for planning a source network charge storage project according to claim 1, wherein in the step2, the constraint conditions include an electric power balance constraint, an operation constraint of wind power equipment, an operation constraint of photovoltaic equipment, an operation constraint of grid equipment and a charging rate constraint of stored energy.
3. The annual operation evaluation method for planning a source network load storage project according to claim 2, wherein in the step3, the objective function is the minimum annual total cost.
4. The annual operation evaluation method for planning source network load storage projects according to claim 1, wherein in the step 4, the number of groups is 200, and the number of iterations is 10-15.
5. The annual hourly operation evaluation method of source network charge storage project planning according to claim 1, wherein said step 5 comprises the steps of:
step 5.1, initializing all particles, wherein each particle is an array formed by the hours of the current month, and each number represents the charge and discharge plan of the energy storage device at the current moment;
step 5.2, calculating an objective function for all particles to obtain objective function values of different particles;
step 5.3, taking the objective function value with the smallest objective function value as the population optimal function value of the round, and taking the particle corresponding to the optimal function value as the optimal particle of the round;
step 5.4, generating N particles of the next round by adopting the following formula:
v f (t+1)=V f (t)+c 1 ×rand×(pbest f -SOC f )+c 2 ×rand×(gbest-SOC f )SOC f (t+1)=SOC f +V f (t+1)
wherein v is f (t+1) represents the speed of the f-th particle at the time of updating, v f (t) represents the speed of the f-th particle at this time, c 1 Represents a self-learning factor, rand is a random number between 0 and 1, pbest f For the f-th particle self-historic optimal position, SOC f For the f-th particle, c2 represents a population learning factor, and gbest is the historical optimal particle position of the population;
step 5.5, judging the position of each particle, and when the position or the speed of the particle exceeds a limit value, pulling the particle back to the boundary;
step 5.6, repeating the steps 5.2 to 5.5 until the maximum iteration times are reached;
and 5.7, outputting the final population optimal particle position as a result of a particle swarm algorithm, namely an optimal charging and discharging plan of the energy storage device.
6. The annual operation evaluation method for planning a source network charge storage project according to claim 1, wherein in the step 6, the power grid purchase cost is as follows: the sum of the product of the main grid power per hour and the electricity price per hour throughout the year.
7. The annual operation evaluation method for planning of source network load storage project according to claim 1, wherein in said step 6, new energy waste amount P is calculated by hour calculation or sum calculation abd,t ,
The calculation formula of the hour-by-hour algorithm is as follows:
P abd,t =P WT,t +P PV,t +P BESS,t -P Load,t
wherein P is WT,t For the t hour wind power, P PV,t Is the firstt hours photovoltaic power, P BESS,t Power of energy storage device at t hour, P Load,t The power used for loading at the t hour;
the sum calculation method is as follows:
P abd,sum =P WT,sum +P PV,sum+ P Grid,sum -P Load,sum
wherein P is abd , sum Represents the new energy waste amount of the whole year, P WT,sum Represents total annual wind power generation amount, P PV,sum Representing total photovoltaic power generation amount, P Grid,sum Representing total power purchase amount from power grid throughout the year, PL oad,sum Representing the total annual load power consumption.
8. The annual operation evaluation method for planning a source network charge storage project according to claim 1, wherein in the step 6, the new energy power rejection rate γ is abd The calculation method of (1) is as follows:
wherein P is abd,sum Represents the new energy waste amount of the whole year, P WT,sum Represents total annual wind power generation amount, P PV,sum Representing the total amount of photovoltaic power generation throughout the year.
9. The annual operation evaluation method for planning a source network charge storage project according to claim 1, wherein in the step 6, the new energy electric quantity is a ratio gamma green The calculation method comprises the following steps:
wherein P is abd,sum Represents the new energy waste amount of the whole year, P WT,sum Represents total annual wind power generation amount, P PV,sum Representing total photovoltaic power generation amount, P Load,sum Representing the total annual load power consumption.
10. An annual operation evaluation device for planning a source network load storage project is characterized by comprising:
the input module is used for acquiring basic data, wherein the basic data comprises capacity configuration of wind power and photovoltaic, and wind power, photovoltaic output and load information of each year and every hour;
the processing module is used for storing an objective function and constraint conditions, and solving the objective function and constraint conditions according to a particle swarm algorithm to obtain an optimal charge-discharge plan of the energy storage device;
and the output module is used for calculating new energy technical indexes according to the basic data and the optimal charging and discharging plan of the energy storage device, wherein the new energy technical indexes comprise the power grid electricity purchasing cost, the new energy electricity discarding quantity, the new energy electricity discarding rate and the new energy electricity occupying ratio.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211714372.1A CN116205337A (en) | 2022-12-29 | 2022-12-29 | Annual operation evaluation method and device for source network charge storage project planning |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211714372.1A CN116205337A (en) | 2022-12-29 | 2022-12-29 | Annual operation evaluation method and device for source network charge storage project planning |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116205337A true CN116205337A (en) | 2023-06-02 |
Family
ID=86512120
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211714372.1A Pending CN116205337A (en) | 2022-12-29 | 2022-12-29 | Annual operation evaluation method and device for source network charge storage project planning |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116205337A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117852927A (en) * | 2024-03-06 | 2024-04-09 | 中能智新科技产业发展有限公司 | Source network charge storage integrated planning production simulation method, system and storage medium |
-
2022
- 2022-12-29 CN CN202211714372.1A patent/CN116205337A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117852927A (en) * | 2024-03-06 | 2024-04-09 | 中能智新科技产业发展有限公司 | Source network charge storage integrated planning production simulation method, system and storage medium |
CN117852927B (en) * | 2024-03-06 | 2024-06-04 | 中能智新科技产业发展有限公司 | Source network charge storage integrated planning production simulation method, system and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112103946B (en) | Particle swarm algorithm-based microgrid energy storage optimization configuration method | |
CN108960510B (en) | Virtual power plant optimization trading strategy device based on two-stage random planning | |
WO2023274425A1 (en) | Multi-energy capacity optimization configuration method for wind-solar-water-fire storage system | |
CN112131733B (en) | Distributed power supply planning method considering influence of charging load of electric automobile | |
CN112084652B (en) | Capacity configuration calculation method and system for wind-solar-energy-storage complementary power generation system | |
CN112508221A (en) | Day-ahead scheduling decision method considering source-load uncertainty under limited energy storage | |
CN111401604B (en) | Power system load power prediction method and energy storage power station power distribution method | |
CN111404193B (en) | Data-driven-based microgrid random robust optimization scheduling method | |
CN109492748A (en) | A kind of Mid-long term load forecasting method for establishing model of the electric system based on convolutional neural networks | |
CN113507110A (en) | Distributed energy storage cluster optimization control method for improving new energy consumption in power distribution network | |
CN116205337A (en) | Annual operation evaluation method and device for source network charge storage project planning | |
CN113011101B (en) | Control method and system for energy storage to participate in frequency modulation auxiliary service optimization | |
CN109617048B (en) | Power grid planning typical scene selection method based on multi-target linear programming | |
CN116388252B (en) | Wind farm energy storage capacity optimal configuration method, system, computer equipment and medium | |
CN113675890A (en) | TD 3-based new energy microgrid optimization method | |
CN113794199A (en) | Maximum profit optimization method of wind power energy storage system considering electric power market fluctuation | |
CN116826752A (en) | Multi-objective low-carbon loss reduction optimization scheduling strategy method for energy consumption of transformer area | |
CN113469412A (en) | Real-time operation strategy optimization method and system for comprehensive energy system | |
CN114971048A (en) | Zero-carbon park comprehensive energy planning simulation method and system and computer equipment | |
CN112670982B (en) | Active power scheduling control method and system for micro-grid based on reward mechanism | |
CN116054270A (en) | Power grid dispatching optimization method and device and nonvolatile storage medium | |
CN109447233A (en) | Electric car charge and discharge dispatching method and system | |
CN111092448B (en) | Dynamic optimization method for optimal battery energy storage capacity configuration based on user side load | |
CN113408886B (en) | Energy storage capacity configuration method and device | |
CN115618258B (en) | Method and system for extracting key operation modes of power system planning |
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 |