CN110137955A - A kind of decision-making technique counted and the robust Unit Combination of CVaR is dispatched - Google Patents
A kind of decision-making technique counted and the robust Unit Combination of CVaR is dispatched Download PDFInfo
- Publication number
- CN110137955A CN110137955A CN201910423718.4A CN201910423718A CN110137955A CN 110137955 A CN110137955 A CN 110137955A CN 201910423718 A CN201910423718 A CN 201910423718A CN 110137955 A CN110137955 A CN 110137955A
- Authority
- CN
- China
- Prior art keywords
- constraint
- risk
- unit
- wind
- stage
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000005457 optimization Methods 0.000 claims abstract description 41
- 238000010276 construction Methods 0.000 claims abstract description 5
- 238000000691 measurement method Methods 0.000 claims abstract description 5
- 238000004146 energy storage Methods 0.000 claims description 61
- 238000007599 discharging Methods 0.000 claims description 23
- 230000033228 biological regulation Effects 0.000 claims description 13
- 230000005540 biological transmission Effects 0.000 claims description 10
- 230000014509 gene expression Effects 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 230000009194 climbing Effects 0.000 claims description 8
- 230000005611 electricity Effects 0.000 abstract description 3
- 230000007246 mechanism Effects 0.000 abstract description 3
- 238000005520 cutting process Methods 0.000 abstract description 2
- 230000006978 adaptation Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 21
- 238000004364 calculation method Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000009977 dual effect Effects 0.000 description 6
- 230000009021 linear effect Effects 0.000 description 6
- 238000009826 distribution Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000010248 power generation Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000011217 control strategy Methods 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 239000002356 single layer Substances 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000003462 Bender reaction Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000029087 digestion Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000011541 reaction mixture Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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/386—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of decision-making techniques of the robust Unit Combination scheduling of meter and CVaR, comprising steps of from section, time, the uncertain set of Spatial Dimension construction polyhedron；System system risk that may be present is measured using CVaR measurement method；Based on the uncertain set of polyhedron and robust optimization, building meter and the robust Unit Combination Optimized model that wind-powered electricity generation section can be received；Solving model, and the admissible wind-powered electricity generation range of system that exports optimization, risk of loss, unit open/stop state, abandonment and cutting load amount decision.Method of the present invention compared to given uncertain set boundaries, the boundary of constructed uncertain set is optimization gained, by adjusting uncertainty parameter, the control to Robust Optimization Model conservative can be achieved, and model is based on current electric power system dispatching framework, the adaptation mechanism between Unit Combination and robust optimization is deferred to, can be directly used for the Optimization of Unit Commitment By Improved for solving current electric system.
Description
Technical Field
The invention relates to the technical field of analysis and scheduling of a power system, in particular to a decision method for robust unit combination scheduling considering a conditional value at Risk (CVaR).
Background
The wind power generation is influenced by practical requirements in the aspects of environmental protection, energy and the like, and the wind power generation shows a rapid development trend in China. The renewable energy and environmental pollution-free property of wind power generation are the inherent motivations for the development of the wind power generation, and in addition, the wind power station has the advantages of short construction time, flexible investment and the like, and the construction of the energy is promoted.
However, at present, wind power with high proportion and uncertain attribute is merged into a power grid, so that a plurality of problems are brought to the operation of a power system, the traditional unit combination concept and method are difficult to effectively deal with, the decision is conservative or not aggressive, and the wind abandon and load shedding are caused or more cost is paid. In order to effectively deal with the uncertainty problem of large-scale wind power, consume more wind power and achieve the best energy-saving and emission-reducing effect, the following key problems must be solved:
①, the power system resolves the boundary decision problem of the wind power uncertainty interval;
② measurement of the magnitude of the risk caused by wind power uncertainty exceeding the safety margin that the system can dispatch.
In view of the above, it is desirable to provide a decision method that can be directly used for solving the unit combination scheduling problem of the current power system based on the current power system architecture.
Disclosure of Invention
In order to solve the technical problem, the technical scheme adopted by the invention is to provide a decision method for combining and scheduling a CVaR-based robust unit, which comprises the following steps:
constructing a polyhedral uncertain set from interval, time and space dimensions by balancing unit combination cost and risk loss cost;
based on an interval robust optimization thought and a risk decision theory, a CVaR measurement method is adopted to measure system risks possibly existing in the system;
constructing a robust generator set combination optimization model considering an acceptable wind power interval based on a polyhedron uncertain set and a robust optimization method;
and solving the model, and outputting the wind power range, risk loss, start/stop state of the unit, wind curtailment and load shedding amount decision which can be accepted by the optimized system.
In the above method, the measuring the system risk possibly existing in the system by using the CVaR measurement method includes:
the risk of wind curtailment due to underestimation and the risk of load loss due to overestimation.
In the method, the construction of the robust generator set combination optimization model considering the receivable wind power interval based on the polyhedron uncertain set and the robust optimization method comprises the following steps:
establishing a stage 1 model for the unit combination problem of the day before the unit start/stop state and the operation risk based on the decision variables;
according to the established stage 1 model, establishing a stage 2 model for the system economic operation problem corresponding to the worst scenario in the wind power uncertain set of the unit output, the load shedding amount and the wind abandoning amount at each moment based on a decision variable;
and determining a two-stage robust optimization model considering the acceptable wind power interval.
In the method, the solution model adopts a modified C & CG algorithm.
In the above method, the set of polyhedral uncertainties is:
in the formula, Wwt、Respectively obtaining an actual value, a predicted value and a lower boundary and an upper boundary of an uncertain interval of the wind power plant w at a time t;andrespectively representing auxiliary {0,1} variables of the fluctuation condition of the wind power output; gamma-shapedTAnd ΓSUncertainty parameters of the uncertainty set in space and time are respectively.
In the above method, the CVaR value expression of the wind curtailment risk due to underestimation is as follows:
the CVaR value expression for the risk of loss of load due to overestimation is:
in the formula,and the wind power prediction error is obtained.
In the above method, the stage 1 model specifically includes:
an objective function:
in the formula,is a quadratic cost function of the unit and is linearized in three sections in a model, ag、bg、cgRespectively the coefficients of the quadratic cost function,the output power of the unit g in the 1 st stage in the time period t;
is a typical exponential start-up cost function;
respectively representing the cost coefficients of charging and discharging of the energy storage system e in the time period t;
g and W are the number of the units and the wind power plant in the system respectively; k is a penalty coefficient;
the constraint conditions include:
the method comprises the following steps of unit output power upper and lower limit constraint, unit climbing rate constraint, unit minimum on-off time constraint, active power balance constraint, node active power balance and transmission capacity constraint, risk constraint, energy storage system charging/discharging state constraint, energy storage system charging/discharging power constraint, energy storage system capacity constraint and energy storage system charging/discharging regulation strategy constraint.
In the above method, the phase 2 model specifically includes:
an objective function:
in the formula, cwt、cdtRespectively the cost of wind abandoning and load shedding;
ΔDdtis the d thThe load shedding amount of the load in the time t;
ΔWwtthe wind curtailment quantity of the w wind power plant in the time period t is obtained;
the constraint conditions include: the method comprises the following steps of unit output power upper and lower limit constraint, unit climbing rate constraint, active power balance constraint, load shedding amount constraint, air curtailment amount constraint and node active power balance constraint.
In the method, the two-stage robust optimization model considering the acceptable wind power interval specifically comprises the following steps:
in the formula, F1Is the objective function of stage 1; f2Is the objective function of the 2 nd stage;
C1represents the constraints satisfied by the stage 1 variables; c2 represents the constraint satisfied by the stage 2 variable;
ugta decision variable 0/1 in the 1 st stage represents the running state of the unit g in the time period t, 1 represents running, and 0 represents shutdown;andrespectively representing auxiliary {0,1} variables of the fluctuation condition of the wind power output;
Pgt、ΔDdtand Δ WwtIs a continuous variable of stage 2, wherein PgtThe output power of the unit g in the time period t is obtained; delta DdtThe load shedding amount of the d load in the time period t; Δ WwtThe wind curtailment quantity of the w wind power plant in the time period t is obtained;
Wwt、respectively the actual value of the wind farm w output in the time period t, the lower boundary and the upper boundary of the uncertainty interval.
Compared with a method for giving an uncertain set boundary, the boundary of the established uncertain set is obtained through optimization, conservatism control of a robust optimization model can be realized by adjusting uncertainty parameters, and the model is based on a current power system scheduling framework, follows a self-adaptive mechanism between unit combination and robust optimization, and can be directly used for solving the unit combination problem of the current power system.
Drawings
FIG. 1 is a flow chart provided by the present invention;
FIG. 2 is a probability density function normal distribution diagram of a wind power plant prediction error in the invention;
FIG. 3 is a flow chart for solving a C & CG algorithm for solving a two-stage robust optimization unit combination model in the invention;
FIG. 4 is a wiring diagram of a 6-node system in case analysis according to the present invention;
FIG. 5 is a graph of a wind power prediction value, a wind power prediction error band, and a range of extinction of an uncertain set in the case of the present invention;
FIG. 6 is a wind power acceptance range curve diagram under different node systems in the case of the invention;
(a) the wind power receiving range curve diagram under the 6-node system, and the wind power receiving range curve diagram under the 6-node system;
fig. 7 is a graph of the scheduling total cost and risk cost for different node systems and different risk thresholds in the case of the present invention;
(a) the operation cost and risk cost curve graphs of the 6-node system and the systems under different risk thresholds, (b) the operation cost and risk cost curve graphs of the 118-node system and the systems under different risk thresholds;
fig. 8 is a diagram illustrating a scheduling result of whether to consider the energy storage control strategy in the present invention;
(a) the energy storage scheduling result diagram in the 1 st stage, and (b) and (c) are energy storage scheduling result diagrams with or without constraint formulas (34) and (35), respectively.
Detailed Description
The invention provides a two-stage robust unit combination decision model considering CVaR (conditional value at risk value) based on the idea that an energy storage system deals with wind power uncertainty to reduce system risk, and the invention is described in detail by combining a specific implementation mode and an attached drawing of the specification.
As shown in fig. 1, the present invention provides a decision method for combining and scheduling a CVaR-based robust unit, which includes the following steps:
s1, building a dimensional polyhedral uncertainty set from an interval, time and space by balancing unit combination cost and risk loss cost, and realizing flexibility adjustment of the uncertainty set, thereby avoiding conservatism of a robust scheduling result.
The polyhedron uncertainty set is specifically as follows:
in the formula, Wwt、Respectively obtaining an actual value, a predicted value and a lower boundary and an upper boundary of an uncertain interval of the wind power plant w at a time t;andrespectively representing auxiliary {0,1} variables of the fluctuation condition of the wind power output; gamma-shapedTAnd ΓSUncertainty parameters of the uncertainty set in space and time are respectively.
The formula (1) represents the actual wind power output, which is an uncertain set represented by a polyhedral structure
The formula (2) and the formula (3) are used for controlling the scene covered by the uncertain set, and reasonably selecting the parameter gammaTAnd ΓSTo reduce the conservatism of robust scheduling.
The equation (4) represents pole constraint of a polyhedron, and the same wind power can only reach an uncertain upper limit value or lower limit value at the same time but cannot reach the upper limit value or lower limit value at the same time.
S2, based on the interval robust optimization thought and the risk decision theory, measuring the system risk possibly existing in the system by adopting a CVaR measuring method; the system risk comprises a wind abandoning risk and a load shedding risk, and specifically comprises the following steps:
in this embodiment, the wind power prediction error is:
assuming that the mean obedience of wind power prediction errors is 0 and the variance is sigmaw 2The probability density function of the normal distribution of (2) is shown in fig. 2.
In the embodiment, the maximum wind power uncertainty region capable of being accepted by the system is represented by adopting the CVaR calculation formulaThe average potential loss caused by the situation of (1) is, as shown in the shaded part in fig. 2, that is, under the condition of fully exploiting the maximum regulation capacity of the system, the average loss caused by the wind power fluctuation exceeding the regulation and control range of the system is defined as CVaR accepted by the wind power of the system.
If the actual output power of the wind power exceeds the upper boundary of the regulation and control range of the system, the safety of the system operation can be ensured by taking a wind abandoning measure. Therefore, the CVaR value expression of the wind curtailment risk due to underestimation is:
in the formula,the expected wind curtailment risk average corresponding to the shaded portion on the right in fig. 2.
Similarly, if the actual output power of the wind power is lower than the lower boundary of the regulation and control range of the system, a load shedding measure can be taken to ensure the safety of the operation of the system. At this time, the CVaR value expression of the risk of loss of load due to overestimation is:
in the formula (8), the reaction mixture is,corresponding to the expected load loss risk average for the left shaded portion in fig. 2.
However, both the formula (7) and the formula (8) are non-linear integral expressions, which are difficult to solve directly by a commercial solver, and thus need to be linearized; therefore, a linearization means is adopted to carry out linearization processing on the CVaR value expression of the wind power abandoned wind risk and the load loss risk. Expression after linear transformation:
in the formula,representing a CVaR value of wind power wind abandon risk caused by underestimation of a wind power plant w in a time period t;CVaR values representing the risk of load shedding due to overestimation of the wind farm w during the time period t, respectively.
In summary, the equations (9) to (12) establish a functional correspondence relationship between the wind power uncertainty interval receivable by the system and the CVaR value faced by the system.
S3, constructing a robust generator set combination optimization model considering an acceptable wind power interval based on a polyhedron uncertain set and a robust optimization method; the method specifically comprises the following steps:
the embodiment is based on a two-stage robust optimization problem established by a polyhedron uncertain set:
s31, establishing a stage 1 model for the day-ahead unit combination problem of the unit start/stop state and the operation risk based on the decision variables;
s32, according to the established stage 1 model, establishing a stage 2 model for the system economic operation problem corresponding to the worst situation in the wind power uncertain set of the unit output, the load shedding amount and the wind abandoning amount at each moment based on the decision variables.
And S33, determining a two-stage robust optimization model considering the receivable wind power interval according to the step S31 and the step S32. Wherein,
(1) stage 1 model
In the stage 1, the starting and stopping states of the unit, the wind power output range accepted by the system, the energy storage related decision quantity and the like are decided.
The objective function of the stage 1 is to minimize the start-stop cost, the operation cost, the energy storage charge-discharge cost and the operation risk cost of the system of the unit; namely, it is
An objective function:
in the formula,is a quadratic cost function of the unit and is linearized in three sections in a model, ag、bg、cgRespectively, the secondary costCoefficients of the function;the output power of the unit g in the 1 st stage in the time period t;
the method is a typical exponential type starting cost function, and can be linearized by adopting a three-section piecewise linearization method for convenient analysis;
respectively representing the cost coefficients of charging and discharging of the energy storage system e in the time period t;
g and W are the number of the units and the wind power plant in the system respectively; and K is a penalty coefficient.
The decision maker in the embodiment can balance the reliability and the running economy of the system by adjusting the penalty coefficient K, so that the range of the wind power which can be consumed by the system is determined;
the constraint conditions include:
① upper and lower limit constraint of unit output power
In the formula,the maximum output power and the minimum output power allowed by the unit g are respectively.
② unit ramp rate constraint
In the formula,the upward and downward climbing rates of the unit g are respectively.
③ unit minimum on-off time constraints
In the formula,respectively the time that the unit g is started and stopped at the initial moment;the minimum start-up and shut-down time of the unit g.
④ active power balance constraints
In the formula, DdtRepresenting the active power predicted value of the load d in the time period t; d is the total load.
⑤ node active power balance and transmission capacity constraint
In the formula, BijIs the line admittance between node i and node j; thetaitIs the phase angle of node i during time t; f. ofij,tThe active transmission power of the power transmission line between the node i and the node j is obtained;the limit value of the active transmission power of the power transmission line; b and L are the number of nodes and transmission lines in the system respectively.
⑥ Risk constraints
In the formula, RiskdaA day-ahead risk target threshold value;representing the installed capacity of the wind farm w.
In this embodiment, the charging and discharging behavior of the energy storage system, that is, the energy storage may be used as a power supply to provide power to the system, or may absorb power from the system as a load. The energy storage system has operation flexibility, so that the wind power uncertainty can be effectively dealt with by matching with a conventional thermal power generating unit so as to reduce the operation risk of the system. Furthermore, from the perspective of spatial location, the decentralized access of renewable energy sources inevitably leads to the spatially decentralized configuration of energy storage that deals with its uncertainty, which is also an energy source that effectively deals with network constraints. From the time process, the stored energy is also a bridge for energy transfer, and the function of the stored energy is essentially standby transmission. Therefore, the scheduling strategy of the energy storage not only affects the consumption capacity of the system to the wind power, but also affects the safety of the system.
The energy storage scheduling model is as follows:
⑦ energy storage system charge/discharge state constraints
In the formula,respectively, the status signs that the energy storage system is in discharging or charging. The constraint may be applied to ensure that the energy storage system is only in a charged or discharged state at a time.
⑧ energy storage system charge/discharge power constraints
In the formula,respectively, the discharging power and the charging power of the energy storage system in the phase 1 model and the allowed maximum discharging power and maximum charging power thereof.
⑨ energy storage system capacity constraints
In the formula, EetStoring energy value of energy storage system for time period t ηch、ηdRespectively charging and discharging efficiencies of the energy storage system;maximum and minimum values allowed by energy storage of the energy storage system; ee,int、Ee,TThe electric quantity values of the energy storage system in the initial period and the final period are respectively. In order to ensure that the energy storage system can function normally in the next scheduling period, the electric quantity value of the energy storage system at the end of each period is required to be equal to that of the initial period.
⑩ energy storage system charge/discharge regulation strategy constraints
Considering that the wind power output is higher than the predicted value at the peak load period or lower than the predicted value at the valley load period, the influence on the peak regulation of the power grid is not obvious. The embodiment limits the charging and discharging modes of the energy storage system under the two conditions to avoid the condition of discharging in a low-ebb period or charging in a high-peak period, thereby reducing the loss caused by the disordered charging and discharging of the energy storage system.
In the formula,the discharge power and the charge power of the energy storage system in the second-stage model are respectively.
(1) Stage 2 model
And the stage 2 is to minimize wind curtailment and load shedding cost corresponding to the worst case described by the uncertain set when the unit start-stop state obtained by the decision in the stage 1 and the boundary of the wind power receivable interval are given. A two-stage robust optimization model considering an acceptable wind power interval is constructed based on adaptive robust optimization and comprises the following steps:
an objective function:
in the formula, cwt、cdtThe cost of the wind curtailment and the load shedding are respectively, and the values of the embodiment are 10$/MW and 1000 $/MWh.
Constraint conditions are as follows:
① upper and lower limit constraint of unit output power
② unit ramp rate constraint
③ active power balance constraints
④ load shedding amount constraint
⑤ air curtailment restriction
⑥ node active power balance constraint
The 2 nd-stage model is formed by equations (36) to (43), energy storage system regulation and control strategy constraints (30) to (35), uncertain concentration constraints (1) to (5) and transmission capacity constraints (21) to (24).
Based on the thought, a two-stage robust optimization model considering the acceptable wind power interval is given:
in the formula, F1For the purpose of stage 1A scaling function; f2Is the objective function of the 2 nd stage;
C1represents the constraints satisfied by the stage 1 variables; c2 represents the constraint satisfied by the stage 2 variable;
ugta decision variable 0/1 in the 1 st stage represents the running state of the unit g in the time period t, 1 represents running, and 0 represents shutdown;
Pgt、ΔDdtand Δ WwtIs a continuous variable of stage 2, wherein PgtThe output power of the unit g in the time period t is obtained; delta DdtThe load shedding amount of the d load in the time period t; Δ WwtThe wind curtailment quantity of the w wind power plant in the time period t is obtained.
And S4, solving the model, and outputting the wind power range, risk loss, start/stop state of the unit, wind abandon and load shedding amount decisions which can be accepted by the optimized system.
The constructed two-stage problem forms a CVaR-considered two-stage robust optimization unit combination model, but the model has a multilayer structure and cannot be directly solved. For convenience of expression, a model expressed in a matrix form is given, and the model is solved by using a modified C & CG (column and Constraint Generation) algorithm.
Stage 1: main Problem (MP)
Decision variables for stage 1 include ugt、Andbecause no wind abandon and load shedding phenomena occur in the 1 st stage, u can be ensuredgt、Feasibility of the solution andexistence of a solution.
Wherein x represents a binary state vector of the generator set;
vector quantityAnd vector y represents a continuous vector of generator output and a phase angle vector of each node respectively;
u represents a binary state vector of the charge-discharge state of the energy storage system;
and z represents the charge and discharge power and capacity vector of the energy storage system;
w represents a wind farm output boundary vector;
q represents an operational risk vector;
a, B, C, D, E, A, B, C, D, E and F are all corresponding constant coefficient matrixes; for example, the matrix a can be derived from the constraint equations (14) to (20). Compared with the traditional robust unit combination, the proposed model introduces more decision variables and constraints, such as w, u,and z, which will add significantly to the computational complexity. The main problem is the MILP problem, and thus can be solved by the existing commercial solver CPLEX.
Stage 2: sub-problem (SP)
The decision variables of stage 2 include Pgt、ΔDdtAnd Δ WwtAnd uncertain set of wind power
In the formula, s represents a curtailment wind and a tangential load vector, and v describes a binary vector of the wind power uncertain set;
f, G, H, G, H, I, J, K, L, M and N are all corresponding constant coefficient matrixes;
representing the Hadamard product.
In this embodiment, since the sub-problem is a max-min structure, the model of this structure cannot be directly solved, and the max-min structure model needs to be converted into a single-layer MILP problem for convenient solution.
Firstly, a mode of solving the dual problem of the internal minimization problem is converted into a max problem, according to Wangdan and the like, in the report of China Motor engineering in 2014, "household temperature control load demand response and energy efficiency power plant modeling considering user comfort constraint" are provided, the mentioned strong dual theory is used for solving the dual problem of the min problem of the inner layer, then the dual problem is combined with the problem of the outer layer, and finally the dual problem is converted into a single-layer max problem. The matrix form of the transformed subproblem is expressed as follows:
λT≥0 (49)
Nv≤h (50)
in the formula, λ is a dual variable of the constraint formulae (37) to (43). It was observed that the objective function (47) contained a bilinear term λTv, an external approximation method provided by ' theoretical research of unit combination in power market ' provided by doctor thesis in 2006 or a characteristic and research framework of automatic demand response in intelligent power utilization ' provided by prosperous et al in an automation journal of a power system in 2013 can be adopted, and a large M linear method is used for solving.
In some cases, the outer approximation may not be able to find a globally optimal solution. Therefore, the bilinear term lambda is solved by adopting a large M linearization methodTv：
By the above process, the sub-problem (46) translates to the standard single layer MILP problem, as shown in equation (55):
in the formula,is an auxiliary vector. q is a continuous vector;
equations (52) through (54) are auxiliary constraints generated by the large M method, so that the bilinear terms can be equivalently converted into the following equations:
after the above linear transformation, the main problem (45) and the sub-problem (55) constitute a standard two-stage MILP problem, and this example uses Ross t.mewton et al, 2011 to propose "Green power volume research: price elasticity and policy analysis, and an improved C & CG algorithm is used for solving the two-stage unit combination model. Compared with Benders decomposition algorithm, the algorithm has the advantages that the iteration times are greatly reduced, and the convergence rate is greatly improved. The overall flow of the algorithm solution is shown in fig. 3, and the steps are as follows:
a1, initialization
Setting the lower bound LB as 0, the upper bound UB as + ∞, and the convergence error epsilon as more than or equal to 0; the iteration number k is equal to 0, and the solution space is equal to O; go to step A2;
a2, solving the following main problem formula (45);
ψ≥fTs (57)
get the optimal solutionUpdatingGo to step A3;
a3, solving the subproblems
Based on given unit combination state xk+1、uk+1、wk+1Solving a subproblem (55) to obtain the worst fluctuation scene of wind power outputAnd the air abandoning quantity and the load cutting quantity sk+1；
Updating the Upper bound
If UB-LB is less than or equal to epsilon, stopping iteration and outputting an optimal solution;
otherwise, setting l to l +1, and solving the sub-problem to the worst wind power output scene corresponding to the optimal solutionPassing to the main question, go to step A4;
a4, adding variables and constraints
Increasing the variable yk+1And sk+1And adding constraint equations (59) and (60) and returning the constraint equations to the main problem, updating l +1 and O ∪ { k +1}, solving the outer layer problem of the updated sub-problem and returning to the step A2.
ψ≥fTsk+1(59)
The beneficial effects of the embodiment are that:
the embodiment provides a two-stage robust unit combination optimization model considering CVaR, and has the following beneficial effects:
1) compared with the method for giving the uncertain set boundary, the boundary of the uncertain set constructed by the embodiment is obtained by optimization, and the conservative control of the robust optimization model can be realized by adjusting the uncertainty parameter; meanwhile, a proper risk threshold value can be selected, or compromise is carried out between the system operation risk and the system operation cost, so that an optimal scheduling solution can be obtained.
2) The configuration of the energy storage increases the flexibility resources of the thermal power generating unit, improves the capability of the system for coping with wind power uncertainty, and relieves the limitation of the climbing rate of the system and the limitation of the network restriction.
3) The model is based on the current power system scheduling architecture, follows a self-adaptive mechanism between unit combination and robust optimization, can be directly used for solving the unit combination problem of the current power system, and verifies the correctness of the model. Meanwhile, aiming at the characteristics of the model of the embodiment, the C & CG algorithm is adopted for solving, and the advantages of the method in the aspects of calculation efficiency and calculation speed are verified.
The present embodiment will be described below by way of specific examples.
The present case takes a 6-node system and a modified IEEE 118-node system as examples, and analyzes the effectiveness of the model proposed in the above embodiments. As shown in fig. 4, the system is a 6-node system, and the system includes 1 wind farm with a capacity of 250MW, 1 energy storage system, 3 thermal power generating units, and 7 lines. Technical data of capacity, climbing rate, lines and the like of a thermal Power generating unit refer to Qiaoyan Bian et al and provide 'distribution robust solution to the redundant scheduling with partial wind Power information' on IEEE Transactions on Power Systems (International Electrical and electronics Engineers institute of Electrical and electronics Engineers) in 2015 (a robust solution method for distribution of partial wind Power information is considered in backup scheduling).
The modified IEEE118 node system comprises 53 generators, 91 load nodes, 186 lines, 14 capacitors and 9 tap changers; 3 wind power plants with the capacity of 250MWh are respectively connected with nodes 59, 66 and 94. The detailed IEEE118 node system data may refer to motor.ec.it.edu/data/scuc _ 118. The parameters of the energy storage system are shown in the table 1, and the charging/discharging price of the energy storage system is 0.4/0.6($/kWh) respectively. The test calculation adopts Visual Studio 2016C + + software to call a CPLEX12.8 solver to solve, and the computer is configured with a Win10 system, Intel Core i7-8700k series, a main frequency of 3.0GHz and a memory of 16G. The simulation time scale was 1 day divided into 24 periods.
TABLE 1 energy storage System parameters
The wind power predicted value and the wind power predicted error band adopted by the calculation example are shown in FIG. 5. in the calculation example, the selected confidence levels are β respectivelyT95% and βS95%, i.e. in a 6-node systemAnd
in a 118-node systemWind power prediction error delta wwtThe mean value of (2) is 0 and the variance is obtained as follows.
(1) Scheduling result considering combination of CVaR robust unit
Taking a 6-node test system as an example, the scheduling results of the scheduling model taking into account the risk robust unit combination model proposed by the model of the present embodiment and the document Negash a.i. et al in ieee power & Energy Society General Meeting (the institute of electrical and electronics engineers, council of electrical and Energy), on "Optimizing and responding to the price and quantity", in 2014, are compared and analyzed, see table 2.
TABLE 2 comparison of scheduling results
As can be seen from table 2, comparing the scheduling results of the literature, the most economical unit G1 in the model of the present embodiment always operates in the research period, the most expensive unit G2 exits from operating in 11: 00-12: 00 time intervals, the unit G3 is in an operating state only in 9:00, 11: 00-12: 00 time intervals of the load peak, the corresponding total system cost is 94236.95$, and the savings of 94236.95$ -94632.321 $ 395.37; the peak clipping and valley filling functions of the energy storage system are explained, on one hand, the thermal power generating unit can provide enough flexible resources to deal with the uncertainty of wind power, and therefore the wind power consumption capability of the system is improved. On the other hand, the pressure of standby configuration of the thermal power generating unit at the peak moment can be reduced, the starting/stopping times of the unit are reduced, the unit is used as an extra flexibility source, and the UC scheduling cost is reduced.
(2) Comparison with other unit combination models
Now, the model of the present embodiment is discussed to be respectively compared and analyzed with deterministic unit combinations (10% uncertainty interval), given robust unit combinations of symmetric uncertainty sets, and robust unit combinations considering asymmetric uncertainty sets considering risks, and scheduling results obtained by these 4 models are shown in table 3 below. In 2001, "Using utility information formation to a customer demand management information referring modules" (model for calibrating customer demand management behavior Using utility information) was proposed by fahrioglum et al in IEEE Transactions on Power Systems (the institute of electrical and electronics engineers electrical Power Systems).
TABLE 3 scheduling results of different unit combination models
As can be seen from table 3, the operation cost and risk cost of the model of the present embodiment are lower than those of the other three models, which indicates that the model of the present embodiment has better operation flexibility capability and capability of reducing the operation risk of the system. In addition, even if the configuration of the energy storage system increases the flexibility adjustment capability of the system, some wind abandoning or load shedding amount still exists, which is the problem that the system cannot cope with the intermittence of high-proportion wind power at a certain moment due to the limitation of the climbing rate of the thermal power generating unit.
By comparing the iteration times and the running time, taking the 118 node as an example, the total calculation time of the model of the embodiment is 8.677 seconds, and the iteration times is 3 times, which shows that the method of the embodiment has high calculation efficiency due to the fact that the linear property of the model is maintained, and can meet the calculation efficiency requirement of solving the actual system.
(3) Flexibility uncertainty set and range of absorption
The uncertain set of wind power of the model of the embodiment is flexible, adjustable and asymmetric, and mainly comprises two reasons:
(1) in order to limit the conservatism of the robust method, the robustness of the system is limited from a time dimension formula (2) and a space dimension formula (3) respectively.
(2) The asymmetry of the wind power absorption range, namely the asymmetry of the wind power uncertain set, is caused by setting different wind abandon penalty costs and load shedding penalty costs. As can be seen from table 3, even though the load shedding amount is much smaller than the wind curtailment amount, the load shedding cost is larger than the wind curtailment cost because we set a higher load shedding penalty cost, thereby making the load shedding risk lower.
The range of consumption obtained by optimizing the model in this embodiment is compared with the scheduling result of the risk robust unit combination model, as shown in fig. 6.
As can be seen from fig. 6, compared with the literature, the upper and lower boundaries of the digestion wind power decided by the model of the present embodiment are smaller in most of the time period. For example, in a 6-node system, time periods 1: 00-3: 00, 5:00, 7:00, 9:00, 12:00, 14:00, 16:00, 18:00, 21: 00-22: 00;
in the 118-node system, the time period is 2:00, 8:00, 11: 00-16: 00. The uncertainty of the system with more schedulable flexible resources for handling the wind power output is proved by configuring the energy storage system, so that the conservatism of the traditional robust method is reduced on one hand, and the running risk of the system is also reduced while the safe running of the system is ensured on the other hand.
(4) Uncertain set conservation analysis
Adjusting an uncertainty parameter ΓTAnd ΓSThe robustness of the model can be controlled, so that the conservatism of the scheduling result is reduced. To analyze the effect of uncertainty parameter variation on the conservation of the results obtained in the model of this example, Table 4 lists the results for ΓTAnd ΓSDifferent combinations of the calculated results.
TABLE 4 gammaTAnd ΓSScheduling results under different combinations
As can be seen from Table 4, at fixed gammaSInvariant, with ΓTThe total operating cost is gradually increased, but the risk cost is gradually reduced, and the magnitude of the operating risk is always smaller than or equal to a given risk threshold value. It can also be seen that when fixing ΓTThe wider the wind farm (Γ)SThe larger the total operation cost is), the smaller the total operation cost is, the strategy of wind power plant geographical distribution is increased, the network constraint limitation is relaxed, and the robustness of the system for wind power uncertainty can be improved; but the risk cost is increasing gradually because: in order to improve the robustness of the system to wind power uncertainty, an over-conservative strategy is adopted, and a suboptimal solution is generated. Obviously, ΓTAnd ΓSThe value of (A) has a decisive role in the robustness and conservation of the robust optimization result. Gamma-shapedTAnd ΓSIf the value is too large, the robustness of the optimization result is good, but the optimization result is extremely conservative and has poor economy; and gamma isTAnd ΓSThe values are made too small, so that too few wind farms have the allowable output reaching the prediction boundary, which will not sufficiently reflect the uncertainty of the wind power.
Therefore, by reasonably setting the uncertainty parameter, the adjustable robust optimization method can eliminate the situation of extremely low probability of the uncertainty parameter, so that the model is more practical and more suitable for engineering practice. Table 4 also lists the uncertainty set Γ at different wind farmsTAnd ΓSThe following calculation time again demonstrates the effectiveness of the method.
(5) Risk threshold RiskdaInfluence of (2)
In an embodiment, the Risk threshold RiskdaIs also an important parameter, which affects not only the scheduling strategy but also the feasibility of the scheduling solution. In practice, the selection may be made based on historical data, risk preferences of the decision maker, electricity contract, and the like. FIG. 7 shows Risk at different Risk thresholdsdaThe running cost of the system and the trend of the risk cost curve.
As shown in FIG. 7, taking the 118-node system as an example, with the Risk threshold RiskdaGradually, the operating cost of the system gradually decreases, while the risk cost gradually increases. Risk threshold value Risk of systemdaWhen the flow rate is about 1000, the running cost and the risk cost of the system tend to be smooth. Working systemRisk threshold value RiskdaThe risk cost of the system is slightly reduced when the cost is more than 1400 deg. Conversely, when the Risk threshold of the system is lowered to RiskdaWhen the value is less than or equal to 10, the system has no feasible solution. This means that the minimum feasible risk threshold for the system is 10. At the same time, it can be observed that the running cost and Risk cost of the system and the Risk threshold RiskdaThe relationship between is not strictly linear. It is also possible to derive upper and lower limits for the risk threshold for a feasible solution to exist for the system.
(6) Effect of energy storage charging/discharging control strategy
Taking a 6-node system as an example, fig. 8 shows the roles of constraints (34) and (35) of the energy storage charging/discharging regulation strategy.
Fig. 8(a) shows the energy storage scheduling result in phase 1. Therefore, during the peak period of power utilization, the stored energy is used for responding to the increased load demand through discharging, and during the valley period of power utilization, the reverse peak regulation characteristic of wind power is responded through charging the stored energy. The load curve is smoothed by the peak clipping and valley filling effects of the stored energy. Fig. 8(b) and 8(c) are comparisons of the energy storage scheduling results of the presence or absence of the constraint equations (34) and (35). It can be seen that the energy storage charge/discharge regulation strategy constrains the charge/discharge state of the energy storage system during the net load peak and valley periods to be limited: the unit combination plan before the day of the 1 st stage requires that the energy storage system can only be charged in the valley period and can only be discharged in the peak period. In the real-time scheduling link of the 2 nd stage, if the actual output of the wind power is higher than the predicted output, redundant wind power can be charged for the energy storage system, however, if the actual output of the wind power is lower than the predicted output, the charging power can be reduced, and discharging cannot be performed. The constraint not only ensures that the scheduling strategy of the energy storage system is consistent with the peak regulation of the power grid, but also avoids the loss caused by disordered charging and discharging of the energy storage.
The present invention is not limited to the above-mentioned preferred embodiments, and any structural changes made under the teaching of the present invention shall fall within the protection scope of the present invention, which has the same or similar technical solutions as the present invention.
Claims (9)
1. A decision method for considering CVaR robust unit combination scheduling is characterized by comprising the following steps:
constructing a polyhedral uncertain set from interval, time and space dimensions by balancing unit combination cost and risk loss cost;
based on an interval robust optimization thought and a risk decision theory, a CVaR measurement method is adopted to measure system risks possibly existing in the system;
constructing a robust generator set combination optimization model considering an acceptable wind power interval based on a polyhedron uncertain set and a robust optimization method;
and solving the model, and outputting the wind power range, risk loss, start/stop state of the unit, wind curtailment and load shedding amount decision which can be accepted by the optimized system.
2. The decision method as claimed in claim 1, wherein the measuring the system risk that the system may have using the CVaR measurement method comprises:
the risk of wind curtailment due to underestimation and the risk of load loss due to overestimation.
3. The decision method according to claim 1, wherein the construction of the robust wind turbine combination optimization model taking into account the acceptable wind power interval based on the polyhedral uncertain set and the robust optimization method comprises the following steps:
establishing a stage 1 model for the unit combination problem of the day before the unit start/stop state and the operation risk based on the decision variables;
according to the established stage 1 model, establishing a stage 2 model for the system economic operation problem corresponding to the worst scenario in the wind power uncertain set of the unit output, the load shedding amount and the wind abandoning amount at each moment based on a decision variable;
and determining a two-stage robust optimization model considering the acceptable wind power interval.
4. The decision method of claim 1, wherein the solution model employs a modified C & CG algorithm.
5. The decision method of claim 2, wherein the set of polyhedral uncertainties is:
in the formula, Wwt、Respectively obtaining an actual value, a predicted value and a lower boundary and an upper boundary of an uncertain interval of the wind power plant w at a time t;andrespectively representing auxiliary {0,1} variables of the fluctuation condition of the wind power output; gamma-shapedTAnd ΓSUncertainty parameters of the uncertainty set in space and time are respectively.
6. Decision method according to claim 5 characterized in that the CVaR value expression of the wind curtailment risk due to underestimation is:
the CVaR value expression for the risk of loss of load due to overestimation is:
in the formula,and the wind power prediction error is obtained.
7. The decision method according to claim 3, wherein the stage 1 model is specifically:
an objective function:
in the formula,is a quadratic cost function of the unit and is linearized in three sections in a model, ag、bg、cgRespectively the coefficients of the quadratic cost function,the output power of the unit g in the 1 st stage in the time period t;
is a typical exponential start-up cost function;
respectively representing the cost coefficients of charging and discharging of the energy storage system e in the time period t;
g and W are the number of the units and the wind power plant in the system respectively; k is a penalty coefficient;
the constraint conditions include:
the method comprises the following steps of unit output power upper and lower limit constraint, unit climbing rate constraint, unit minimum on-off time constraint, active power balance constraint, node active power balance and transmission capacity constraint, risk constraint, energy storage system charging/discharging state constraint, energy storage system charging/discharging power constraint, energy storage system capacity constraint and energy storage system charging/discharging regulation strategy constraint.
8. The decision method according to claim 7, wherein the stage 2 model is specifically:
an objective function:
in the formula, cwt、cdtRespectively the cost of wind abandoning and load shedding;
ΔDdtthe load shedding amount of the d load in the time period t;
ΔWwtthe wind curtailment quantity of the w wind power plant in the time period t is obtained;
the constraint conditions include: the method comprises the following steps of unit output power upper and lower limit constraint, unit climbing rate constraint, active power balance constraint, load shedding amount constraint, air curtailment amount constraint and node active power balance constraint.
9. The decision method according to claim 3, wherein the two-stage robust optimization model taking into account the acceptable wind power interval is specifically:
in the formula, F1Is the objective function of stage 1; f2Is the objective function of the 2 nd stage;
C1represents the constraints satisfied by the stage 1 variables; c2 represents the constraint satisfied by the stage 2 variable;
ugta decision variable 0/1 in the 1 st stage represents the running state of the unit g in the time period t, 1 represents running, and 0 represents shutdown;andrespectively representing auxiliary {0,1} variables of the fluctuation condition of the wind power output;
Pgt、ΔDdtand Δ WwtIs a continuous variable of stage 2, wherein PgtThe output power of the unit g in the time period t is obtained; delta DdtThe load shedding amount of the d load in the time period t; Δ WwtThe wind curtailment quantity of the w wind power plant in the time period t is obtained;
Wwt、respectively the actual value of the wind farm w output in the time period t, the lower boundary and the upper boundary of the uncertainty interval.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910423718.4A CN110137955B (en) | 2019-05-21 | 2019-05-21 | Decision method for robust unit combination scheduling considering CVaR |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910423718.4A CN110137955B (en) | 2019-05-21 | 2019-05-21 | Decision method for robust unit combination scheduling considering CVaR |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110137955A true CN110137955A (en) | 2019-08-16 |
| CN110137955B CN110137955B (en) | 2021-04-02 |
Family
ID=67572142
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910423718.4A Active CN110137955B (en) | 2019-05-21 | 2019-05-21 | Decision method for robust unit combination scheduling considering CVaR |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110137955B (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110492534A (en) * | 2019-08-23 | 2019-11-22 | 国网新疆电力有限公司经济技术研究院 | Meter and the random optimization dispatching method of electric system containing wind-powered electricity generation of flexibility |
| CN110889634A (en) * | 2019-11-29 | 2020-03-17 | 四川大学 | Water resource global optimization configuration method based on two-layer multi-follower robust optimization |
| CN110909959A (en) * | 2019-12-16 | 2020-03-24 | 国网青海省电力公司电力科学研究院 | Wind power operation risk-considering multi-energy complementary power system robust optimization method |
| CN111126664A (en) * | 2019-11-25 | 2020-05-08 | 广西电网有限责任公司 | Active power distribution network alternating current power flow management method based on robust optimization |
| CN112054554A (en) * | 2020-08-18 | 2020-12-08 | 国网山东省电力公司临沂供电公司 | Non-parameter statistics-based adaptive distribution robust unit combination method and system |
| CN113381400A (en) * | 2021-05-28 | 2021-09-10 | 国网青海省电力公司 | Method and device for evaluating capacity of storing and storing new energy |
| CN113393054A (en) * | 2021-07-05 | 2021-09-14 | 华北电力大学 | Optimal scheduling method and optimal scheduling system of wind storage combined system |
| CN114123316A (en) * | 2021-11-10 | 2022-03-01 | 云南电网有限责任公司 | SCUC model optimization method and device considering wind power uncertainty |
| CN114372655A (en) * | 2021-11-10 | 2022-04-19 | 北京交通大学 | Emergency rescue material dynamic allocation method based on distributed robust optimization technology |
| CN114611754A (en) * | 2022-02-09 | 2022-06-10 | 上海奉贤燃机发电有限公司 | Distributed power supply and energy storage virtual power plant robustness optimization method considering risks |
| CN117039924A (en) * | 2023-07-12 | 2023-11-10 | 三峡大学 | Micro-grid group two-stage low-frequency load shedding method considering risk avoidance |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106127389A (en) * | 2016-06-24 | 2016-11-16 | 河海大学 | A kind of virtual plant combined heat and power scheduling Robust Optimization Model |
| CN107239863A (en) * | 2017-04-12 | 2017-10-10 | 广东电网有限责任公司电力调度控制中心 | The robust Unit Combination method of power system security constraints |
-
2019
- 2019-05-21 CN CN201910423718.4A patent/CN110137955B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106127389A (en) * | 2016-06-24 | 2016-11-16 | 河海大学 | A kind of virtual plant combined heat and power scheduling Robust Optimization Model |
| CN107239863A (en) * | 2017-04-12 | 2017-10-10 | 广东电网有限责任公司电力调度控制中心 | The robust Unit Combination method of power system security constraints |
Non-Patent Citations (2)
| Title |
|---|
| YUMIN ZHANG: "Risk–Based Admissibility Analysis of Wind Power Integration into Power System With Energy Storage System", 《ELECTRIC POWER COMPONENTS AND SYSTEMS》 * |
| 张玉敏: "基于狄利克雷模型的分布鲁棒机组组合", 《中国电机工程学报》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110492534A (en) * | 2019-08-23 | 2019-11-22 | 国网新疆电力有限公司经济技术研究院 | Meter and the random optimization dispatching method of electric system containing wind-powered electricity generation of flexibility |
| CN111126664A (en) * | 2019-11-25 | 2020-05-08 | 广西电网有限责任公司 | Active power distribution network alternating current power flow management method based on robust optimization |
| CN110889634B (en) * | 2019-11-29 | 2022-06-14 | 四川大学 | Water resource global optimization configuration method based on two-layer multi-follower robust optimization |
| CN110889634A (en) * | 2019-11-29 | 2020-03-17 | 四川大学 | Water resource global optimization configuration method based on two-layer multi-follower robust optimization |
| CN110909959A (en) * | 2019-12-16 | 2020-03-24 | 国网青海省电力公司电力科学研究院 | Wind power operation risk-considering multi-energy complementary power system robust optimization method |
| CN112054554A (en) * | 2020-08-18 | 2020-12-08 | 国网山东省电力公司临沂供电公司 | Non-parameter statistics-based adaptive distribution robust unit combination method and system |
| CN113381400A (en) * | 2021-05-28 | 2021-09-10 | 国网青海省电力公司 | Method and device for evaluating capacity of storing and storing new energy |
| CN113381400B (en) * | 2021-05-28 | 2022-09-27 | 国网青海省电力公司 | Method and device for evaluating capacity of storing and storing new energy |
| CN113393054A (en) * | 2021-07-05 | 2021-09-14 | 华北电力大学 | Optimal scheduling method and optimal scheduling system of wind storage combined system |
| CN113393054B (en) * | 2021-07-05 | 2023-11-24 | 华北电力大学 | Optimal scheduling method and optimal scheduling system for wind-storage combined system |
| CN114372655A (en) * | 2021-11-10 | 2022-04-19 | 北京交通大学 | Emergency rescue material dynamic allocation method based on distributed robust optimization technology |
| CN114123316B (en) * | 2021-11-10 | 2023-06-13 | 云南电网有限责任公司 | SCUC model optimization method and equipment considering wind power uncertainty |
| CN114123316A (en) * | 2021-11-10 | 2022-03-01 | 云南电网有限责任公司 | SCUC model optimization method and device considering wind power uncertainty |
| CN114611754A (en) * | 2022-02-09 | 2022-06-10 | 上海奉贤燃机发电有限公司 | Distributed power supply and energy storage virtual power plant robustness optimization method considering risks |
| CN117039924A (en) * | 2023-07-12 | 2023-11-10 | 三峡大学 | Micro-grid group two-stage low-frequency load shedding method considering risk avoidance |
| CN117039924B (en) * | 2023-07-12 | 2024-05-28 | 三峡大学 | Micro-grid group two-stage low-frequency load shedding method considering risk avoidance |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110137955B (en) | 2021-04-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110137955B (en) | Decision method for robust unit combination scheduling considering CVaR | |
| Sedhom et al. | IoT-based optimal demand side management and control scheme for smart microgrid | |
| US9489701B2 (en) | Adaptive energy management system | |
| Ding et al. | Optimal dispatching strategy for user-side integrated energy system considering multiservice of energy storage | |
| CN110690702B (en) | Active power distribution network optimal scheduling and operation method considering comprehensive bearing capacity | |
| Mohamed et al. | Microgrid online management and balancing using multiobjective optimization | |
| CN108564305A (en) | Energy Internet advertising distribution integration requirement side response method and device | |
| CN108197766A (en) | A kind of active distribution network Optimal Operation Model for including micro-capacitance sensor group | |
| US11693376B1 (en) | Smart green power node | |
| CN108964048B (en) | Microgrid capacity configuration and operation scheduling method and device | |
| Lu et al. | Two-stage robust scheduling and real-time load control of community microgrid with multiple uncertainties | |
| Verrilli et al. | Stochastic model predictive control for optimal energy management of district heating power plants | |
| Quijano et al. | Probabilistic rolling-optimization control for coordinating the operation of electric springs in microgrids with renewable distributed generation | |
| Zhang et al. | An innovative real-time price based distributed optimal energy management of multi-microgrids in a smart distribution system | |
| CN116961008A (en) | Micro-grid capacity double-layer optimization method considering power spring and load demand response | |
| CN115935619A (en) | Demand response-based day-ahead low-carbon scheduling method and device for active power distribution network | |
| Zhou et al. | Power System Flexibility: Modeling, Optimization and Mechanism Design | |
| Liang et al. | Adaptive robust energy management strategy for campus-based commercial buildings considering comprehensive comfort levels | |
| CN114498625A (en) | Income prediction method and system of wind-solar-storage integrated power supply | |
| CN114188942A (en) | Power grid dispatching method comprising large-scale new energy base | |
| CN114204599A (en) | Distributed energy storage stationing operation control method for inhibiting voltage fluctuation of power distribution network | |
| Zhang et al. | An Optimal Dispatching Model of Active Distribution Network Based on the Collaboration of Source-Load-Storage | |
| CN111724259B (en) | Energy and rotary reserve market clearing method considering multiple uncertainties | |
| CN114219148B (en) | Industrial micro-grid demand response optimization method and device | |
| Mohamed et al. | Environmental/Economic Power Dispatch of Micro Grid Using Multiobjective Optimization |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |