CN115526401A - Novel power supply optimal planning method for electric power system based on digital twinning - Google Patents
Novel power supply optimal planning method for electric power system based on digital twinning Download PDFInfo
- Publication number
- CN115526401A CN115526401A CN202211201159.0A CN202211201159A CN115526401A CN 115526401 A CN115526401 A CN 115526401A CN 202211201159 A CN202211201159 A CN 202211201159A CN 115526401 A CN115526401 A CN 115526401A
- Authority
- CN
- China
- Prior art keywords
- power
- planning
- wind
- power supply
- constraint
- 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
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000004422 calculation algorithm Methods 0.000 claims description 29
- 239000002245 particle Substances 0.000 claims description 22
- 230000002068 genetic effect Effects 0.000 claims description 10
- 238000005457 optimization Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 238000010276 construction Methods 0.000 claims description 8
- 210000000349 chromosome Anatomy 0.000 claims description 6
- 238000004146 energy storage Methods 0.000 claims description 6
- 230000007246 mechanism Effects 0.000 claims description 6
- 210000004027 cell Anatomy 0.000 claims description 5
- 230000007613 environmental effect Effects 0.000 claims description 5
- 238000010248 power generation Methods 0.000 claims description 5
- 230000009194 climbing Effects 0.000 claims description 4
- 238000005265 energy consumption Methods 0.000 claims description 4
- 239000000446 fuel Substances 0.000 claims description 4
- KFLQGJQSLUYUBF-WOJBJXKFSA-N Phyllanthin Chemical compound C([C@H](COC)[C@@H](COC)CC=1C=C(OC)C(OC)=CC=1)C1=CC=C(OC)C(OC)=C1 KFLQGJQSLUYUBF-WOJBJXKFSA-N 0.000 claims description 3
- KFLQGJQSLUYUBF-PMACEKPBSA-N Phyllanthin Natural products C([C@@H](COC)[C@H](COC)CC=1C=C(OC)C(OC)=CC=1)C1=CC=C(OC)C(OC)=C1 KFLQGJQSLUYUBF-PMACEKPBSA-N 0.000 claims description 3
- 238000004140 cleaning Methods 0.000 claims description 3
- 238000011161 development Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims description 3
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 2
- 229910052799 carbon Inorganic materials 0.000 claims description 2
- 238000003912 environmental pollution Methods 0.000 claims description 2
- 238000012423 maintenance Methods 0.000 claims description 2
- 108090000623 proteins and genes Proteins 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 239000007789 gas Substances 0.000 claims 2
- 229930002875 chlorophyll Natural products 0.000 claims 1
- 235000019804 chlorophyll Nutrition 0.000 claims 1
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 claims 1
- 238000007599 discharging Methods 0.000 claims 1
- 230000003287 optical effect Effects 0.000 claims 1
- 239000002912 waste gas Substances 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 6
- 238000004088 simulation Methods 0.000 abstract description 6
- 238000010521 absorption reaction Methods 0.000 abstract description 3
- 238000004891 communication Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 230000003993 interaction Effects 0.000 abstract description 2
- 238000013528 artificial neural network Methods 0.000 description 9
- 230000008859 change Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000002411 adverse Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241001130943 Phyllanthus <Aves> Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- 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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Software Systems (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Artificial Intelligence (AREA)
- Development Economics (AREA)
- Molecular Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Mathematical Physics (AREA)
- Game Theory and Decision Science (AREA)
- Computing Systems (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Computational Linguistics (AREA)
- Biomedical Technology (AREA)
- Evolutionary Biology (AREA)
- Educational Administration (AREA)
- Bioinformatics & Computational Biology (AREA)
- Water Supply & Treatment (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Primary Health Care (AREA)
Abstract
The invention provides a novel power supply optimal planning method for an electric power system based on digital twins, which solves the load fluctuation of a power grid and improves the absorption capacity of the power grid by proportioning the power supply capacity. The strategy solves the problem of uncertainty of new energy output by establishing a five-dimensional new power system model and modeling wind power and photovoltaic units by using a mechanism-data driven hybrid model, establishes a structure planning model and a system planning model of a virtual power system, simulates the system planning modeling by digital twins, feeds back a calculation result to a planning main body to guide planning decision, realizes automatic updating iteration of model structure and parameters by real-time communication and interaction of sensors and physical entity information, verifies the effect of a planning scheme by combining indexes such as system planning cost, wind curtailment quantity and the like, provides detailed information of operating points of the power system, and searches for the optimal occupation ratio of each type of power supply by using multiple times of simulation after updating parameters.
Description
Technical Field
The invention relates to a novel power supply optimal planning method for an electric power system, so as to stabilize load fluctuation of a power grid and improve the consumption capacity of the power grid.
Background
The core problem of the novel power system is the multilayer mixed integer nonlinear programming problem, and the problems of construction cost, resource assessment, environmental cost and the like of a plurality of systems need to be considered, wherein the most important problem is to solve the problem of multi-system coupling. The novel power system is different from the traditional power system in that the ratio of the new energy to the power grid is higher than fifty percent. Because of the influence of factors such as wind speed, wind direction and solar radiation intensity, huge uncertainty and volatility are brought to wind power and photovoltaic power, huge problems are brought in the planning of new energy access to a power grid, and the problems are mainly reflected in poor electric energy quality, operation stability of a power system, reduction of safety, difficulty in power transmission and the like, so that accurate prediction of output of new energy is guaranteed to be the basis of power supply planning of a novel power system. The novel power system planning comprises two parts, wherein one part is used for planning the occupation ratio of various power supply types of a power system, and the other part is used for considering the type of the commissioning equipment from space and time after the occupation ratio is determined.
The digital twin is a virtual model constructed from multiple physical fields, multiple space-time, multiple angles, multiple subjects fused and the like in a digital space, the virtual model monitors the dynamic evolution process of the physical model in real time through a sensor, data analysis, processing, prediction, diagnosis and the like are carried out by combining big data, and the operation and decision of a physical entity are controlled in the whole life cycle of a product through the digital model.
A novel power system digital model is established in a digital twin, and the digital twin has main functions in the planning of a novel power system: modeling and simulating system planning, and feeding back a result to a main body planning to guide planning decision, wherein the model interactively stores a large amount of data in real time along with communication-physical entities, the physical information comprises regional perennial power supply planning information, regional power supply installed capacity and quantity, regional user electricity utilization information, wind power photovoltaic power information, regional generator set information and the like, a structure planning model and a system planning model are established in a virtual model, continuous self-learning is carried out in big data, and virtual-real updating iteration is carried out to obtain a wind-light output predicted value; load prediction and power output balance, system stability and wind-abandoning light-abandoning quantity of a project area are calculated in a structure planning model, an optimal power occupation scheme of various types of power supplies of a novel power system is obtained, the consumption capacity of the power system is improved, and safer and more stable operation of a power grid system is guaranteed.
Disclosure of Invention
Firstly, constructing a novel power system framework under a digital twin, wherein the novel power system framework consists of five-dimensional models: the method comprises the steps of establishing a structure planning model and a system planning model of the virtual power system, establishing a digital twin model for system planning modeling simulation, feeding back a calculation result to a planning main body to guide planning decision, realizing automatic updating and iteration of a model structure and parameters through real-time communication interaction of a sensor and physical entity information, calculating wind-light power prediction by using an algorithm and combining a mass data optimization virtual model, determining the proportion of each type of power supply, evaluating the operation stability of a novel power system planning scheme by using the digital twin model, checking the effect of the planning scheme by combining indexes such as system planning cost, light curtailment quantity and the like, providing detailed information of the operation point of the power system, and searching the optimal proportion of each type of power supply by using multiple times of simulation after parameter updating.
The technical problem of the invention is mainly realized by the following technical scheme:
a power optimal planning method frame of a novel power system power optimal planning method based on digital twinning is characterized in that: comprises the following steps of (a) carrying out,
step 1: step 1, establishing a digital twin five-dimensional space model, as shown in fig. 1:
U DT ＝(PE,VM,DD,CN,SS) (1)
in the formula: PE is a physical entity, VM is a virtual model, DD is a digital twin database, CN is data connection, and SS is a service; a physical entity system may consist of a device or product, a physical system, an activity process, or even an organization. They act according to the laws of physics and cope with uncertain circumstances. Firstly, obtaining novel power system structure planning and system planning information, wherein the information comprises installed capacity and quantity of regional power supplies, power system frequency, voltage amplitude, phase angle, power grid architecture, load type and the like. The virtual model is real-time mapping of a physical entity and can be formulated by a geometric model, a behavior model or a rule model, the digital twin database comprises static attribute data and dynamic data attributes, part of data is generated by virtual simulation, a simulation result is reflected, and a system planning main body is guided by the result to carry out further planning; the service module comprises application services such as load side demand and power system monitoring, power transmission line diagnosis, prediction and health management.
Step 2: modeling a virtual power system by using a mechanism-data hybrid model, modeling mechanisms of a wind turbine generator and a photovoltaic generator, and correcting model errors by data driving in a parallel mode. The wind turbine generator mainly comprises blades, a rotating wheel, a controller, a generator and the like, and according to a leaf element momentum theory, the wind turbine generator modeling method is used for modeling the torque of the wind wheel:
in the formula:
wind wheel torques of all wind turbine generators in the region y; r and r0 are respectively the radius of the wind wheel and the radius of the hub; c is the chord length of the phyllanthus; a is the inclination angle of the wind direction relative to the wind wheel; ρ is the gas density; v tot Is the relative velocity of the leaf element; l 1 ,l d The gas lift coefficient and the resistance coefficient of the phyllanthin are taken as the parameters;
the mechanism modeling of the photovoltaic unit is as follows:
in the formula: i is the theoretical output current of the photovoltaic cell; i is sc ,U oc ,I m ,U m ,P max ,K I ,K V Short-circuit current, open-circuit voltage, maximum power point electric power, maximum power point voltage, maximum power and electric temperature system tested by photovoltaic cell manufacturers under IEC conditionsNumber and voltage temperature coefficients, etc.;
in practical engineering, due to the influence of adverse conditions such as climate and unit service life, the deviation between the output result and the actual value is large because an operation curve of a photovoltaic unit model and a wind turbine unit model cannot be predicted, so that an empirical model can be established for correcting the deviation in data driving in a parallel mode, the data driving is modeled by a neural network, a hybrid model parallel mode is shown in fig. 2, and the operation parameter of physical equipment is Y 0 The primary output parameter of the mechanism model is Y 1 The data-driven model output is d y ，Y 2 And outputting the mixed model. In the whole model operation, let ERROR coefficient ERROR = Y 0 -Y 1 -d y As minimum as possible, i.e. the function satisfies Y as much as possible 0 ＝d y +Y 1 ＝Y 2 And when the neural network training is finished, the output of the wind turbine generator and the photovoltaic generator after error correction is finally obtained.
The main performance parameters of the wind turbine generator are as follows: wind wheel rotational speed, generator rotational speed and generator power, photovoltaic unit main parameter is: outputting voltage and current as network output layer; the wind turbine generator system has the following influence parameters: wind speed, wind direction, blade cross-section, horizontal axis wind speed, photovoltaic unit influence parameter: illumination intensity, humidity, photovoltaic cell life, module connection loss, the above-mentioned parameter is the network input layer.
The updating process of the electric power system model based on the digital twins is mainly divided into two processes, wherein in one process, a physical system object directly transmits digital twins information to the model through a sensor technology, the model carries out super-parameter iterative computation under an intelligent neural network algorithm, so that the regular change of electric quantities such as voltage, frequency, transmission line loss and the like of a novel electric power system full life cycle is realized through a real-time evolution model, a digital twins database is compared with measured data through historical operating data, and the electric power system planning parameters in a certain area are predicted by using the intelligent neural network algorithm, so that the process reflects the predictability of the digital twins; the other process is that the influence of a power system model under extreme conditions is deduced manually, the boundary conditions which are not beneficial to the stable operation of the power system under the rare conditions in reality are assumed manually, the key information of the power system under various adverse conditions is changed through an intelligent neural algorithm, the information obtained under the assumed conditions is compared in the intelligent artificial algorithm, and an optimal planning scheme is selected, and the process reflects the hypothesis of digital twinning;
and 3, step 3: the novel structure planning modeling of the power system establishes a target function of a structure planning model of the virtual power system by taking economy, abandoned wind and abandoned light quantity and system fluctuation rate as references, and the total operation cost is as follows:
in the formula: t is the planning interval year;
the system construction, system operation and maintenance, unit fuel and environmental pollution cost in the t year are respectively expressed as follows:
in the formula: c is a novel power system power supply type;
O t,c,i representing new increase of a certain class in planning yearInstalled capacity and total installed capacity of a unit i of a type power supply c;
represents the construction cost of a certain type of newly added power supply capacity;
represents the operating cost of a certain type of power supply and the i variable operating power cost;
representing the cost of the power consumption fuel of the power supply c unit i;
representing the cost of electricity consumption carbon emission and the environmental cleaning cost of a certain power supply c unit i; f t,c,i ＝B t,c,i O t,c,i Representing the annual output of a certain type of power supply m unit i; b is t,c,i Representing the total number of annual operating hours of a certain type of power m unit i.
The wind and light abandoning amount is:
wherein, the first and the second end of the pipe are connected with each other,
is the abandoned wind power of the wind power plant after the installed capacity after T interval years,
the light power is the light power of the installed capacity of the photoelectric plant after T interval years;
the fluctuation rate of the wind, light, fire and storage combined power generation system is as follows:
wherein the content of the first and second substances,
and the new installed capacity generating capacity of the wind power plant, the photovoltaic power plant and the stored energy after T intervals.
The structure planning of novel electric power system establishes the constraint condition, its characterized in that, resource development constraint:
K Gt +Y m,i ≥K pt +Y m,i+1 (1.1)
in the formula: k Gt 、K pt Representing the total power consumption and the maximum consumption of the unit in a certain place; y is m,i 、Y m,i+1 Representing the construction time of a certain type of power supply c unit i and i + 1;
the environmental constraints are:
in the formula: f e,t,c,i 、F tot Representing an upper limit on the amount of exhaust gas discharged and dischargeable.
And (3) new energy ratio constraint:
in the formula:
representing the energy consumption of the new energy unit;
representing the new energy ratio; s. the tot Represents the total amount of energy consumption;
the power balance constraint is:
wherein, P load,t Is the total load value of the regional power grid,
for newly adding the output data of the installed capacity of the thermal power in the T interval years, and the method comprises the following steps:
in the formula:
adding a predicted value and a predicted error of the capacity of the wind power system;
the capacity of the newly added photovoltaic power generation system is predicted value and prediction error;
the unit output constraint conditions are as follows:
wherein d is 0 、d 1 、d 2 Respectively is the state coefficient P of peak regulation of the thermal power generating unit in a certain area i (a)、P i (b) Respectively is the lower limit and the upper limit of the output of the peak regulating generator set in different states of the thermal power generating unit, P i (max)、P i (min) respectively representing the upper limit and the lower limit of the output of the regional normal peak regulating thermal power generating unit;
the minimum start-stop time constraint of the thermal power generating unit is as follows:
wherein, the first and the second end of the pipe are connected with each other,
is the current working operation time of t-1 moments of a first group of thermal power generating units in a certain area,
for the current work stop time of the ith group of thermal power generating units at the t-1 th moment in a certain region,
is the starting operation instantaneous time of the ith thermal power generating unit at the tth moment of a certain region, namely the minimum starting time,
is the minimum time from running to stopping of the first group of thermal power generating units in a certain area, u i,t The state variation of the ith thermal power generating unit in the operation period at the tth moment in a certain area is shown, wherein 1 represents the operation of the thermal power generating unit, and 0 represents the shutdown of the thermal power generating unit;
and (3) climbing restraint of the thermal power generating unit:
in the formula:
respectively providing maximum downward and upward climbing power of a newly added thermal power unit in a certain area during capacity working;
the upper limit and the lower limit of the charge and discharge power of the energy storage battery are as follows:
in the formula:
respectively representing the minimum charge-discharge and maximum power of each energy storage battery in a certain area;
novel reserve restraint of electric power system:
in the formula: r 1,t 、R 0,t Respectively representing the positive spare capacity and the negative spare capacity required by the generator sets of other power systems after wind and light power converters are ignored in a certain area; r res,1,t 、R res,0,t Dividing the capacity into a positive spare capacity and a negative spare capacity required during the safe and stable operation of a novel power system representing a certain area;
respectively, represent the confidence level of the unit equipment backup requirement that the new power system must meet.
And 4, step 4: and carrying out simulation analysis by combining the steps, further optimizing indexes such as economy, wind curtailment and light curtailment of the structure planning by utilizing an optimization algorithm based on digital twin modeling driving, carrying out weighted sum on each index by considering multiple indexes such as the economy, the wind curtailment and the system fluctuation rate level of the structure planning model, and optimizing decision variables such as the model number and the number of the equipment and the optimized operation mode of the system equipment under different working conditions. On the basis of inheriting a nonlinear model after system coupling, firstly, coding decision variables by using an intelligent optimization algorithm, searching according to algorithm rules in a solving area, further calculating to obtain decision variables corresponding to positions solved for multiple times, wherein the decision variables represent calculation basic conditions provided by power planning based on digital twins, and finally, carrying out next calculation by using an optimization algorithm according to indexes such as calculation construction cost, absorption capacity and the like of the system planning model after checking is completed through checking structural planning constraints. As shown in fig. 3.
During solving, the virtual model is solved by adopting an improved particle swarm-genetic algorithm, and the algorithm solving flow is as follows:
(1) Firstly, initializing each parameter and particle population;
(2) Each particle performs a large range of motion in the search space;
(3) Judging the constraint condition, if the particle does not meet the constraint condition, returning to the step (2), recalculating the speed and the searching direction of the particle, and updating the particle;
(4) Comparing and updating the individual optimal solution and the global optimal solution of all the particles;
(5) Judging whether the convergence condition of the particle swarm algorithm is met, if the global optimal solution cannot be updated after continuous 25 iterations, determining convergence, indicating that the optimal solution is obtained, transmitting the optimal solution to a genetic algorithm for next-step solving, wherein the genetic algorithm takes all particles as chromosomes, the dimensionalities of the particles as genes of the chromosomes, inputting the historical optimal solutions of all the particles, and moving the initial chromosome set which is taken as the genetic algorithm to the genetic algorithm; if not, returning to the step (2);
(6) Simulating the planning scheme of the virtual-novel power system built by the cloudsps platform to operate so as to select and operate the planning scheme;
(7) Starting cross and variation operation on the virtual-novel power system planning scheme;
(8) Judging whether the constraint condition is met, if so, returning to the step (7);
(9) Judging whether the convergence condition of the genetic algorithm is met, when the global optimal solution is not carried out after continuous 25 iterations, considering the convergence, and if the convergence condition is not met, returning to the step (6); and if so, outputting a final result.
Therefore, the invention has the following advantages: (1) The method comprises the following steps that a visualization model built based on a Cloudsps platform can simulate the operation and energy flow of novel power system equipment such as a photovoltaic generator set, a wind generating set and a thermal generator set according to real-time collected actual scene data; (2) According to the real-time data of the twin database, the actual condition of the result load is optimized, and the planning scheme of the novel power system can be reflected in real time. (3) By adopting a data driving method and combining real-time data of a twin database, the output of the wind generating set and the photovoltaic set can be fully reflected according to the change characteristics of the ring factors, the expression of the model is more accurate, the energy equipment can more flexibly deal with the load change, the investment cost is further optimized, and the fluctuation rate of a novel power system is reduced.
Drawings
Fig. 1 is a digital twinning five-dimensional space diagram.
Fig. 2 is a hybrid model parallel mode.
FIG. 3 is an optimization algorithm based on a digital twin modeling drive.
Fig. 4 is a load prediction diagram.
FIG. 5 is a wind-solar prediction graph.
Fig. 6 is a new energy occupancy map.
Fig. 7 is a graph showing the change in the wind curtailment rate.
Detailed Description
The technical solution of the present invention is further described below by way of examples with reference to the accompanying drawings.
Firstly, the existing installed capacity of new energy, the future installed capacity, power type parameters, wind and light prediction parameters and environment cleaning cost are assumed, wind and light power prediction is contrastively analyzed in a digital twin database, and a conversion model of numerical weather forecast and wind power and photovoltaic output prediction is constructed by adopting a GLSTM neural network. And comparing data in a digital twin database, calculating the values to obtain a prediction deviation if searching to obtain the actual wind-solar output value and the wind-solar output prediction under similar conditions, and compensating the weight by the prediction deviation. And further optimizing the neural network algorithm for the digital twin modeling drive by combining equipment constraint and the like, calculating the optimal combination of the power supply in the novel power system, obtaining the optimal planning scheme of the novel power system, and improving the absorption capacity, the power grid safety and the stability of the novel power system.
To verify the advantages of the present invention, the following simulation analysis was performed:
taking a certain area in China as an example for calculation, assuming that the planning year is 2035 years, the maximum load of the area in 2022 years is 1678 thousands KW, the average load is 1022 thousands KW, and by inputting the characteristic parameters of the highest temperature, the lowest temperature and the weather, the load value of the area within the planning year range is predicted by using a neural network, as shown in FIG. 4, assuming that the spare rate of power supply equipment of the area is set to 8%.
Relevant data such as wind direction, temperature, wind speed and the like are selected according to local numerical weather forecast and input into a digital twin database, and output prediction is carried out on certain time period of wind power and photovoltaic by combining a neural network prediction scheme, as shown in fig. 5. After prediction, planning information is transmitted to a digital twin database through a communication network, data processing, data normalization and data analysis are further transmitted to a digital twin model, decision variables of the model are optimized in a neural network algorithm, power supply occupation ratios of various types every year in a planning year are obtained through digital twin optimization, as shown in fig. 6, and changes of new energy occupation ratios and abandoned wind and abandoned light quantities are obtained through analysis, as shown in fig. 7.
The above embodiments are only intended to illustrate the technical solution of the present invention and not to limit the same, and those skilled in the art to which the present invention pertains can make various modifications or compensations to the described embodiments without departing from the scope of the present invention as defined by the appended claims.
Claims (7)
1. A novel power supply optimal planning method of a power system based on digital twinning is characterized by comprising the following steps: comprises the following steps of (a) carrying out,
establishing a digital twin five-dimensional space model;
modeling a virtual power system by using a mechanism-data hybrid model, modeling mechanisms of a wind turbine generator and a photovoltaic generator, and correcting model errors by data driving in a parallel mode;
establishing a target function and a constraint condition of a structure planning model of the virtual power system according to the economy, the wind curtailment and light curtailment amount and the system fluctuation rate as references;
and solving an objective function by using an optimization algorithm based on digital twin modeling drive, and performing weighted sum on each index, wherein optimization decision variables comprise equipment models, quantity and system equipment optimization operation modes under different working conditions.
2. The optimal planning method for the power supply of the novel power system based on the digital twin as claimed in claim 1, is characterized in that: establishing digital twin five-dimensional space model
U DT ＝(PE,VM,DD,CN,SS) (1)
In the formula: PE is a physical entity, VM is a virtual model, DD is a digital twin database, CN is a data connection, and SS is a service.
3. The optimal planning method for the power supply of the novel power system based on the digital twin as claimed in claim 1, is characterized in that: modeling a virtual power system by using a mechanism-data hybrid model, modeling mechanisms of a wind turbine generator and a photovoltaic generator, correcting model errors by data driving in a parallel mode, and modeling wind wheel torque according to a leaf element momentum theory:
in the formula:
wind wheel torques of all wind turbine generators in the region y; r and r0 are respectively the radius of the wind wheel and the radius of the hub; c is the chord length of the chlorophyll; a is the inclination angle of the wind direction relative to the wind wheel; ρ is the gas density; v tot Is the relative velocity of the phyllanthin; l 1 ,l d The gas lift coefficient and the resistance coefficient of the phyllanthin are taken as the parameters;
the mechanism modeling of the photovoltaic unit is as follows:
in the formula: i is the theoretical output current of the photovoltaic cell; i is sc ,U oc ,I m ,U m ,P max ,K I ,K V Short circuit current, open circuit voltage, maximum for photovoltaic cell manufacturers testing under IEC conditionsPower point power, maximum power point voltage, maximum power, electrical temperature coefficient, and voltage temperature coefficient.
4. The optimal planning method for the power supply of the novel power system based on the digital twin as claimed in claim 1 is characterized in that: establishing an objective function of a structure planning model of the virtual power system according to the reference of economy, the light curtailment amount of the abandoned wind and the system fluctuation rate, wherein the total operation cost is as follows:
in the formula: t is the planning interval year;
the system construction, system operation and maintenance, unit fuel and environmental pollution cost in the t year are respectively expressed as follows:
in the formula: c is a novel power system power supply type;
O t,c,i the installed capacity and the total installed capacity of a certain type of power supply c unit i are newly increased in the planning year;
represents the construction cost of a certain type of newly added power supply capacity;
represents the operating cost of a certain type of power supply and the i variable operating power cost;
representing the cost of the power consumption fuel of the power supply c unit i;
representing the cost of electricity consumption carbon emission and the cost of environmental cleaning of a certain power supply c unit i; f t,c,i ＝B t,c,i O t,c,i Representing the annual output of a certain type of power supply m unit i; b is t,c,i Represents the total annual operating hours of a certain type of power supply m unit i;
the wind and light abandoning amount is:
wherein, P i w,aban (T) is the power of the abandoned wind after the installed capacity of the wind power plant after T years of interval, P i pv,aban (T) is the abandoned optical power of the installed capacity of the photovoltaic plant after T interval years;
the fluctuation rate of the wind, light, fire and storage combined power generation system is as follows:
wherein, P i w (t)、P i pv (t)、P i C (T) New wind, photovoltaic and energy storage plants after T years intervalsAnd (4) installed capacity power generation.
5. The optimal planning method for the power supply of the novel power system based on the digital twin as claimed in claim 1, is characterized in that: and (3) establishing constraint conditions for the structural planning of the novel power system, wherein the constraint conditions comprise resource development constraint, environment constraint, new energy ratio constraint, power balance constraint, unit output constraint, thermal power unit minimum start-stop time constraint, thermal power unit climbing constraint, energy storage battery charge-discharge power upper and lower limit constraint and novel power system standby constraint.
6. The optimal planning method for the power supply of the novel power system based on the digital twin as claimed in claim 1, is characterized in that: the resource development constraints are based on the following formulas:
K Gt +Y m,i ≥K pt +Y m,i+1 (0.1)
in the formula: k Gt 、K pt Representing the total power consumption and the maximum consumption of the unit in a certain place; y is m,i 、Y m,i+1 Representing the construction time of a certain type of power supply c unit i and i + 1;
the environmental constraints are based on the following formula:
in the formula: f e,t,c,i 、F tot Represents the discharge amount and the upper limit of the discharge amount of the waste gas;
the new energy ratio constraint is based on the following formula:
in the formula:
representing the energy consumption of the new energy unit;
representing the new energy ratio; s tot Represents the total amount of energy consumption;
the power balance constraint is based on the following equation:
wherein, P load,t For the total load value, P, of the regional power grid i g (T) is the output data of newly-increased thermal power installed capacity in T interval years, and:
in the formula: p i w,yc (t)、ΔP i w,wc (t) a predicted value and a predicted error of the capacity of the newly-increased wind power system are obtained; p i pv,yc (t)、ΔP i pv,wc (t) a predicted value and a predicted error of the capacity of the newly added photovoltaic power generation system are obtained;
the unit output constraint is based on the following formula:
wherein d is 0 、d 1 、d 2 Respectively is the state coefficient P of peak regulation of thermal power generating unit in a certain region i (a)、P i (b) Respectively as the lower limit and the upper limit of the output power of the peak regulating generator set in different states of the thermal power generating unit P i (max)、P i (min) respectively representing the upper limit and the lower limit of the output of the normal peak-shaving thermal power generating unit in the region;
the minimum start-stop time constraint of the thermal power generating unit is based on the following formula:
wherein the content of the first and second substances,
is the current working operation time at t-1 moments of a first group of thermal power generating units in a certain area,
the current work stop time of the ith thermal power generating unit at the t-1 th moment in a certain region,
is the starting operation instantaneous time of the ith thermal power generating unit at the tth moment of a certain region, namely the minimum starting time,
is the minimum time from running to stopping of the first group of thermal power generating units in a certain area, u i,t The state variation of the ith thermal power generating unit in the operation period at the tth moment in a certain area is shown, wherein 1 represents the operation of the thermal power generating unit, and 0 represents the shutdown of the thermal power generating unit;
the ramp constraint of the thermal power generating unit is based on the following formula:
in the formula:
respectively providing maximum downward and upward climbing power of a newly added thermal power unit in a certain area during capacity working;
the upper limit and the lower limit of the charging and discharging power of the energy storage battery are based on the following formula:
in the formula:
respectively representing the minimum charge-discharge power and the maximum power of each energy storage battery in a certain area;
the novel power system backup constraints are based on the following formula:
in the formula: r 1,t 、R 0,t Respectively representing the positive spare capacity and the negative spare capacity required by the generator sets of other power systems after wind and light power converters are ignored in a certain area; r res,1,t 、R res,0,t Dividing the capacity into a positive spare capacity and a negative spare capacity required during the safe and stable operation of a novel power system representing a certain area;
respectively, represent the confidence level of the unit equipment backup requirement that the new power system must meet.
7. The power optimal planning method of the novel power system based on the digital twin power supply optimal planning method according to claim 1, wherein in the solving, the virtual model is solved by adopting an improved particle swarm-genetic algorithm, and the algorithm solving process is as follows:
(1) Firstly, initializing each parameter and particle population;
(2) Each particle performs a large range of motion in the search space;
(3) Judging the constraint condition, if the particle does not meet the constraint condition, returning to the step (2), recalculating the speed and the searching direction of the particle, and updating the particle;
(4) Comparing and updating the individual optimal solution and the global optimal solution of all the particles;
(5) Judging whether the convergence condition of the particle swarm algorithm is met, if the global optimal solution cannot be updated after continuous 25 iterations, determining convergence, indicating that the optimal solution is obtained, transmitting the optimal solution to a genetic algorithm for next-step solving, wherein the genetic algorithm takes all particles as chromosomes, the dimensionalities of the particles as genes of the chromosomes, inputting the historical optimal solutions of all the particles, and moving the initial chromosome set which is taken as the genetic algorithm to the genetic algorithm; if not, returning to the step (2);
(6) Simulating and operating a planning scheme of a virtual-novel power system built by the cloudsps platform so as to select and operate the planning scheme;
(7) Starting cross and variation operation on the virtual-novel power system planning scheme;
(8) Judging whether the constraint condition is met, if so, returning to the step (7);
(9) Judging whether the convergence condition of the genetic algorithm is met, when the global optimal solution is not carried out after continuous 25 iterations, considering the convergence, and if the convergence condition is not met, returning to the step (6); and if so, outputting a final result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211201159.0A CN115526401A (en) | 2022-09-29 | 2022-09-29 | Novel power supply optimal planning method for electric power system based on digital twinning |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211201159.0A CN115526401A (en) | 2022-09-29 | 2022-09-29 | Novel power supply optimal planning method for electric power system based on digital twinning |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115526401A true CN115526401A (en) | 2022-12-27 |
Family
ID=84698792
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211201159.0A Pending CN115526401A (en) | 2022-09-29 | 2022-09-29 | Novel power supply optimal planning method for electric power system based on digital twinning |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115526401A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116109335A (en) * | 2023-04-10 | 2023-05-12 | 国网浙江省电力有限公司 | New energy dynamic electricity price data management system based on digital twin |
| CN116579902A (en) * | 2023-04-07 | 2023-08-11 | 南京电力设计研究院有限公司 | Digital park electric carbon data mapping method, system, equipment and storage medium |
| CN117318110A (en) * | 2023-11-28 | 2023-12-29 | 惠州市长晟科技有限公司 | Distributed micro electric energy storage system, method, computer equipment and storage medium |
| CN117592386A (en) * | 2024-01-19 | 2024-02-23 | 新立讯科技股份有限公司 | New energy BA standardization implementation quality assessment method and system |
| CN118071387A (en) * | 2024-04-22 | 2024-05-24 | 厦门市盛迅信息技术股份有限公司 | Electric power facility operation cost prediction method and system based on digital twin |
-
2022
- 2022-09-29 CN CN202211201159.0A patent/CN115526401A/en active Pending
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116579902A (en) * | 2023-04-07 | 2023-08-11 | 南京电力设计研究院有限公司 | Digital park electric carbon data mapping method, system, equipment and storage medium |
| CN116579902B (en) * | 2023-04-07 | 2023-12-12 | 南京电力设计研究院有限公司 | Digital park electric carbon data mapping method, system, equipment and storage medium |
| CN116109335A (en) * | 2023-04-10 | 2023-05-12 | 国网浙江省电力有限公司 | New energy dynamic electricity price data management system based on digital twin |
| CN116109335B (en) * | 2023-04-10 | 2023-09-08 | 国网浙江省电力有限公司 | New energy dynamic electricity price data management system based on digital twin |
| CN117318110A (en) * | 2023-11-28 | 2023-12-29 | 惠州市长晟科技有限公司 | Distributed micro electric energy storage system, method, computer equipment and storage medium |
| CN117318110B (en) * | 2023-11-28 | 2024-03-08 | 惠州市长晟科技有限公司 | Distributed micro electric energy storage system, method, computer equipment and storage medium |
| CN117592386A (en) * | 2024-01-19 | 2024-02-23 | 新立讯科技股份有限公司 | New energy BA standardization implementation quality assessment method and system |
| CN117592386B (en) * | 2024-01-19 | 2024-04-05 | 新立讯科技股份有限公司 | New energy BA standardization implementation quality assessment method and system |
| CN118071387A (en) * | 2024-04-22 | 2024-05-24 | 厦门市盛迅信息技术股份有限公司 | Electric power facility operation cost prediction method and system based on digital twin |
| CN118071387B (en) * | 2024-04-22 | 2024-06-21 | 厦门市盛迅信息技术股份有限公司 | Electric power facility operation cost prediction method and system based on digital twin |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115526401A (en) | Novel power supply optimal planning method for electric power system based on digital twinning | |
| CN106327006A (en) | Comprehensive benefit analysis-based micro-power-grid optimal configuration method | |
| CN103138256B (en) | A kind of new energy electric power reduction panorama analytic system and method | |
| CN111210079B (en) | Operation optimization method and system for distributed energy virtual power plant | |
| CN111009895B (en) | Microgrid optimal scheduling method, system and equipment | |
| CN110518580A (en) | Active power distribution network operation optimization method considering micro-grid active optimization | |
| CN112001598A (en) | Energy storage configuration evaluation and operation optimization method for different users based on energy storage type selection | |
| CN114243791A (en) | Multi-objective optimization configuration method, system and storage medium for wind-solar-hydrogen storage system | |
| CN113258561A (en) | Multi-attribute decision-based multi-distributed power supply micro-grid multi-objective optimization scheduling method | |
| CN112836849A (en) | Virtual power plant scheduling method considering wind power uncertainty | |
| CN114301081B (en) | Micro-grid optimization method considering storage battery energy storage life loss and demand response | |
| CN112072643A (en) | Light-storage system online scheduling method based on depth certainty gradient strategy | |
| CN112036735B (en) | Energy storage capacity planning method and system for energy storage system of photovoltaic power station | |
| CN113300400A (en) | Distributed micro-grid scheduling method | |
| CN117277327A (en) | Grid-connected micro-grid optimal energy management method based on intelligent agent | |
| CN115693787B (en) | Method for analyzing new energy acceptance of optical storage and distribution power grid in consideration of source load randomness | |
| CN111342456A (en) | Method and system for modeling energy system of transformer area | |
| CN116961008A (en) | Micro-grid capacity double-layer optimization method considering power spring and load demand response | |
| CN116995740A (en) | Distributed wind power and energy storage optimal configuration method and system for power distribution network | |
| CN114268124B (en) | Distributed power supply credible capacity assessment method based on equal power supply reliability | |
| CN115459323A (en) | Improved genetic algorithm-based microgrid optimization scheduling method | |
| CN112580256B (en) | Distributed power supply location and volume fixing method considering fault rate influence on electric automobile | |
| CN111293719B (en) | AC/DC hybrid micro-grid optimized operation method based on multi-factor evolution algorithm | |
| Gong et al. | Economic dispatching strategy of double lead-acid battery packs considering various factors | |
| CN113591224A (en) | Urban power grid cascading failure risk assessment method and 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 |