CN113793052A - Robust optimization scheduling method for regional comprehensive energy system - Google Patents

Robust optimization scheduling method for regional comprehensive energy system Download PDF

Info

Publication number
CN113793052A
CN113793052A CN202111115252.5A CN202111115252A CN113793052A CN 113793052 A CN113793052 A CN 113793052A CN 202111115252 A CN202111115252 A CN 202111115252A CN 113793052 A CN113793052 A CN 113793052A
Authority
CN
China
Prior art keywords
energy system
regional
robust
model
scheduling
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
Application number
CN202111115252.5A
Other languages
Chinese (zh)
Inventor
黄婧杰
刘科明
刘海军
袁亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changsha University of Science and Technology
Original Assignee
Changsha University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changsha University of Science and Technology filed Critical Changsha University of Science and Technology
Priority to CN202111115252.5A priority Critical patent/CN113793052A/en
Publication of CN113793052A publication Critical patent/CN113793052A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Tourism & Hospitality (AREA)
  • Data Mining & Analysis (AREA)
  • General Business, Economics & Management (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • Computational Mathematics (AREA)
  • Game Theory and Decision Science (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Software Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Algebra (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Databases & Information Systems (AREA)
  • Probability & Statistics with Applications (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Biology (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Educational Administration (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a robust optimization scheduling method for a regional integrated energy system, belonging to the technical field related to regional integrated energy system scheduling and comprising the following specific steps of: constructing a partial load performance function of each energy conversion equipment efficiency; constructing a regional comprehensive energy system model; constructing constraint conditions of a regional comprehensive energy system; constructing an optimized scheduling model; constructing a robust optimization scheduling model, and solving the robust optimization scheduling model to obtain an optimal scheduling scheme of the regional comprehensive energy system; and performing optimized dispatching on the regional comprehensive energy system according to the optimal dispatching scheme. The invention adopts the box-type uncertain set to effectively depict the fitting error and the range of the partial load performance function of the equipment, is applied to the optimization model, changes the scheduling scheme of the system, improves the operation robustness of the system and improves the reliability of heat supply.

Description

Robust optimization scheduling method for regional comprehensive energy system
Technical Field
The invention relates to the technical field of regional integrated energy system scheduling, in particular to a robust optimization scheduling method for a regional integrated energy system.
Background
A Regional Integrated Energy System (RIES) is one of important physical carriers of an energy internet, and is an important way to realize multi-energy complementation and high-utility energy. In actual operation, the efficiency of energy conversion plants such as gas turbines, Gas Boilers (GB) and the like has a significant Part Load Performance (PLP), i.e. the efficiency of the plant changes with the load factor when it is not operating under nominal conditions. The efficiency of most energy conversion devices decreases with decreasing load rate, for example, the power generation efficiency of a gas turbine at low electrical load rate is only 80% or less of full load operation. After the equipment efficiency PLP is considered, the prediction error of the operation cost of the multi-energy system is reduced by 13% compared with the PLP without the equipment efficiency. The PLP of the equipment efficiency can more accurately and effectively depict the running characteristics of the equipment, so that the RIES scheduling scheme is more accurate and reasonable. The impact of the accurate characterization of the device efficiency PLP on the scheduling scheme is not negligible. The function describing the plant efficiency PLP is usually obtained by fitting experimental data or manufacturer supplied test data, which results in an objective presence of fitting errors. In the prior art, a linear function and a power function are mostly adopted to fit the relationship between the heat-electricity ratio and the electric load rate of a Combined Heat and Power (CHP) cogeneration unit, or a linear and quadratic function is adopted to fit a coefficient of performance (COP) of a water chilling unit, which relates to the fitting error of a PLP function, and load loss can be caused by neglecting the error. Therefore, the research on a robust optimization scheduling method of the regional integrated energy system based on the fitting error of the partial load performance function of the equipment is a problem to be solved urgently by the technical personnel in the field.
Disclosure of Invention
In view of the above, the invention provides a robust optimization scheduling method for a regional integrated energy system, which overcomes the defects in the prior art.
In order to achieve the above purpose, the invention provides the following technical scheme:
a robust optimization scheduling method for a regional integrated energy system comprises the following specific steps:
constructing a partial load performance function of each energy conversion equipment efficiency according to the regional comprehensive energy system framework;
constructing a regional comprehensive energy system model according to the regional comprehensive energy system framework and the partial load performance function of each energy conversion equipment efficiency;
constructing a constraint condition of the regional comprehensive energy system for the regional comprehensive energy system model;
constructing an optimized scheduling model according to a partial load performance function of each energy conversion equipment efficiency, a regional comprehensive energy system model and a constraint condition of the regional comprehensive energy system;
introducing fitting errors of partial load performance functions of the efficiency of each energy conversion device into an optimized scheduling model to construct a robust optimized scheduling model;
solving an optimal scheduling scheme of the regional comprehensive energy system by the robust optimization scheduling model;
and performing optimized dispatching on the regional comprehensive energy system according to the optimal dispatching scheme.
Preferably, the expression of the part load performance function of the efficiency of the energy conversion device is:
Figure BDA0003275023880000021
in the formula (I), the compound is shown in the specification,
Figure BDA0003275023880000022
to the efficiency of the plant; and N is the load rate of the equipment, namely the ratio of the current output to the rated output of the equipment.
Preferably, the PLP of the CHP unit power generation efficiency is:
Figure BDA0003275023880000023
in the formula, kCHP,nFitting coefficients of PLP which are the generating efficiency of the CHP unit;
thermoelectric ratio of CHP
Figure BDA0003275023880000024
The relationship with the electrical load rate is:
Figure BDA0003275023880000025
in the formula, kλ,nFitting coefficient of PLP as thermoelectric ratio of CHP;
the fitting coefficients for PLP for thermal efficiency of GB are:
Figure BDA0003275023880000031
in the formula, kGB,nAre fitting coefficients.
The method has the advantages that the partial load performance function of the efficiency of the energy conversion equipment is considered, and the situation that the heat load requirement cannot be met is avoided.
Preferably, based on the regional integrated energy system framework, the regional integrated energy system optimization model is established with the aim of lowest daily operation cost of the regional integrated energy system.
Preferably, the constraints include power balance constraints, energy storage device related constraints, and energy storage device related constraints.
Preferably, the robust optimization scheduling model is constructed by the following steps:
describing a fitting error and an error range of a partial load performance function of the equipment efficiency by adopting an uncertain set;
introducing the uncertain set into an optimized scheduling model to form a compact form of a robust optimized scheduling model;
and converting the compact form of the robust optimized scheduling model into a mixed integer semi-positive definite programming model through Lagrange duality and piecewise linearization, and defining the mixed integer semi-positive definite programming model as the robust optimized scheduling model.
Preferably, the uncertainty set is a box uncertainty set.
Preferably, the method further comprises verifying the heat load loss rate of the robust optimization scheduling model by adopting a Monte Carlo method under the constraint condition.
According to the technical scheme, compared with the prior art, the robust optimization scheduling method for the regional comprehensive energy system is provided, after the fitting error of the efficiency PLP function of the energy conversion equipment is considered, robust scheduling is adopted, and the risk of load loss is reduced; the robust scheduling model based on the fitting error of the PLP function of the energy conversion equipment efficiency can be used for not only counting the influence of single equipment, but also expanding and calculating the uncertainty of the fitting error of the PLP function of the energy conversion equipment efficiency, and further considering the coupling and correlation among the uncertain errors, so that the RIES optimal scheduling scheme is more practical; the box-type uncertain set is adopted to effectively depict the fitting error and the range of the partial load performance function of the equipment, and is applied to an optimization model, so that the system scheduling scheme is changed, the system operation robustness is improved, and the heat supply reliability is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts.
FIG. 1 is a schematic diagram of a framework for a regional energy complex of the present invention;
FIG. 2 is a graph of the CHP thermoelectric ratio of the present invention fitted to its electrical loading rate;
FIG. 3 is a schematic illustration of predicted electrical and thermal loading at a given day in the area of the present invention;
FIG. 4 is a schematic view of the time of day electricity prices at a certain day in the area of the present invention;
FIG. 5 is a diagram illustrating the effect of the uncertainty parameter r on scheduling cost according to the present invention;
FIG. 6 is a schematic diagram of thermal power scheduling in the baseline and robust schemes of the present invention;
FIG. 7 is a schematic diagram of the scheduling of electric power in the baseline and robust schemes of the present invention;
FIG. 8 is a schematic flow chart of the method of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The embodiment of the invention discloses a robust optimization scheduling method for a regional comprehensive energy system, which comprises the following steps as shown in figure 8:
constructing a partial load performance function of each energy conversion equipment efficiency according to the regional comprehensive energy system framework;
constructing a regional comprehensive energy system model according to the regional comprehensive energy system framework and the partial load performance function of each energy conversion equipment efficiency;
constructing a constraint condition of the regional comprehensive energy system for the regional comprehensive energy system model;
constructing an optimized scheduling model according to a partial load performance function of each energy conversion equipment efficiency, a regional comprehensive energy system model and a constraint condition of the regional comprehensive energy system;
introducing fitting errors of partial load performance functions of the efficiency of each energy conversion device into an optimized scheduling model to construct a robust optimized scheduling model;
solving an optimal scheduling scheme of the regional comprehensive energy system by the robust optimization scheduling model;
and performing optimized dispatching on the regional comprehensive energy system according to the optimal dispatching scheme.
Example 1:
1. establishment of partial load performance function of efficiency of each energy conversion device
1.1 Regional Integrated Energy Systems (RIES) framework:
the regional power grid, the heat supply network and the like are coupled together, so that multi-energy complementation is realized, and the energy utilization rate is improved. The rines of the present embodiment mainly includes energy conversion devices such as CHP and GB, and energy storage devices such as a storage battery and a hot water storage tank, and a specific frame is as shown in fig. 1.
1.2 PLP function of energy conversion device efficiency
Considering that each energy conversion device normally operates under a partial load condition during actual operation, PLP of each energy conversion device is calculated by the following formula:
Figure BDA0003275023880000051
in the formula:
Figure BDA0003275023880000052
to the efficiency of the plant; and N is the load rate of the equipment, namely the ratio of the current output to the rated output of the equipment.
For a CHP unit, the PLP of the generating efficiency can be obtained by fitting a fourth-order polynomial, and the method specifically comprises the following steps:
Figure BDA0003275023880000061
in the formula, kCHP,nIs a fitting coefficient;
Figure BDA0003275023880000062
is the electrical loading rate of the CHP.
Thermoelectric ratio of CHP
Figure BDA0003275023880000063
The relation with the electrical load rate can be obtained by fitting a quadratic polynomial, and specifically comprises the following steps:
Figure BDA0003275023880000064
in the formula, kλ,nAre fitting coefficients.
The relationship between the thermal efficiency and the thermal load rate of GB can be obtained by fitting a quadratic polynomial, which specifically comprises the following steps: .
Figure BDA0003275023880000065
In the formula, kGB,nIs a fitting coefficient;
Figure BDA0003275023880000066
is the heat load rate of GB.
1.3 fitting error of PLP function
The plant efficiency PLP function is typically fit from a number of data points derived from either the manufacturer of the plant or from actual measurements from physical simulations. However, regardless of the source of the data, fitting the data points to the PLP function has fitting errors, such as: taking the performance of CHP thermoelectric ratio under different electrical load rates as an example, the rated power P under different load rates provided by Cummins company is quotedNThe fitting result of the data of the relationship between the thermoelectric ratio and the load rate of a plurality of CHP units with the power of 1000-2000kW is shown in FIG. 2 through quadratic function fitting, and as can be seen from FIG. 2, the PLP fitting function of the thermoelectric ratio of the CHP units has obvious fitting error, and actual data points are distributed on two sides of each fitting curve. When the electrical load rate is 90%, PNThe CHP train hot spot ratio, 1100kW, has the greatest deviation (1.56%) of the fitted value from the actual value.
For each energy conversion device PLP function, the fitting error cannot be determined based on the source of the data set used for fitting, the capacity size of the data set, the validity of the data, and the like, and these factors all affect the fitting error. Thus, the fitting error belongs to an uncertain variable that makes it difficult to obtain an accurate probability distribution.
2 RIES optimization model based on equipment PLP function
2.1 objective function
In this embodiment, the lowest daily running cost of the rees is targeted. RIE to establish energy inputs containing only two of electricity and natural gasS, the operation cost comprises the electricity purchasing cost C1Gas purchase cost C2
min=(C1+C2) (5)
Wherein
C1=CELE,tPELE,t (6)
Figure BDA0003275023880000071
In the formula, CELE,tAnd CNG,tRespectively purchasing electricity and gas prices at the time of t time period; pELE,tPurchasing electric power from the power grid for a period t; pCHP,tCHP electrical power for a period t; qGB,tGB thermal power in a t period;
Figure BDA0003275023880000072
and
Figure BDA0003275023880000073
the power generation efficiency of the CHP and the heat supply efficiency of the GB in a t period are respectively; l isHVThe low calorific value of natural gas is usually 9.7kWh/m3(ii) a And T belongs to T, and T is a scheduling period.
2.2 constraint Condition
1) Power balance constraint
PELE,t+PCHP,t+Pd,t=Pc,t+PL,t (8);
Figure BDA0003275023880000074
In the formula, Pc,t、Pd,tThe charging and discharging powers of the electricity storage device are respectively t time period; qc,t、Qd,tRespectively charging and discharging power of the heat storage device in a time period t; pL,t、QL,tRespectively, the electric load and the thermal load in the period t.
2) Energy storage device related constraints
The energy storage equipment for electricity storage, heat storage and the like has a universal mathematical model and is specifically expressed as follows
Figure BDA0003275023880000075
Wherein a first constraint ensures that the energy storage device is within a capacity range, EtRepresenting the capacity of the energy storage equipment in the t period; emin、EmaxThe minimum and maximum capacity of the energy storage device are respectively. The second and the third constraints limit the charge and discharge energy power range, wherein utThe variable is 0-1, the energy storage is only chargeable when the value is 1, and the energy storage is only dischargeable when the value is 0. PES,c,t、PES,d,tEnergy storage charging and discharging power is respectively stored for t time period; pES,c,min、PES,c,maxRespectively storing minimum energy storage power and maximum energy storage power; pES,d,min、PES,d,maxThe minimum energy storage power and the maximum energy discharge power are respectively. The fourth constraint is the relation between the energy storage capacity and the charging and discharging energy power, ec、ed、eselfRespectively charging energy efficiency, discharging energy efficiency and self-energy consumption rate for the energy storage equipment; final constraints guarantee the energy storage scheduling periodicity, Et-1The capacity of the energy storage equipment is in a t-1 time period; eTEnergy storage device capacity for a period of time T; e0Is the initial value of the energy storage capacity.
3) Device force range constraints
0≤PELE,t≤PELE,max (11)
QGB,min≤QGB,t≤QGB,max (12)
PCHP,min≤PCHP,t≤PCHP,max (13)
In the formula, PELE,maxFor the maximum power of purchasing electricity from the power grid, the RIES is specified not to allow electricity to be sold to the power grid; qGB,min、QGB,maxGB minimum and maximum thermal power respectively; pCHP,min、PCHP,maxCHP minimum and maximum electrical power, respectively.
Integrating the above objectives and constraints, a representation of the rees optimized scheduling model based on device efficiency PLP can be obtained as follows:
Figure BDA0003275023880000081
3 RIES robust optimization model based on equipment PLP function fitting error
3.1 Box aggregation robust optimization model
In the RIES optimization scheduling model based on the equipment efficiency PLP, fitting errors exist in the CHP unit heat-electricity ratio PLP function and variable, and the heat load requirement can not be met if the errors are not considered. And for the measurement of the error, a stochastic programming method can be adopted to solve the model formula (14). However, the accurate probability distribution of the fitting error is difficult to obtain, so that the uncertain quantity xi and the variation range of the fitting error are described by adopting an uncertain set U, and then the solution is carried out by adopting a robust optimization method. In this embodiment, a box-type uncertain set is selected to depict fitting errors and their ranges, and has the following structure:
Figure BDA0003275023880000082
in the formula (I), the compound is shown in the specification,
Figure BDA0003275023880000091
xi is the upper and lower limits of the uncertain quantity xi respectively, and a decision maker can adjust the upper and lower limits of the uncertain quantity according to actual conditions.
In another embodiment, fitting errors can be described in various forms such as an ellipsoid, a polyhedron and the like.
The compact form of the RIES robust optimization model based on the fitting error of the device PLP function corresponding to equation (14) is:
Figure BDA0003275023880000092
wherein x is a decision variable vector, and the specific expression is
x=[PELE,t,PCHP,t,Pd,t,Pc,t,QGB,t,Qd,t,Qc,t,Et,ut]T (16);
c is a coefficient column vector of the objective function formula (5); a (xi) is a coefficient matrix containing uncertain quantity xi; G. h is a coefficient matrix of x in the corresponding constraint; a. g and h are constant column vectors. Among the constraint conditions of the formula (15), the first constraint corresponds to the formula (9); the second constraint is an inequality constraint comprising the first three constraints of equation (10) and equations (11) - (13); the third is an equality constraint, which includes two constraints after the formulas (2) - (4), (8) and (10).
Because the first constraint comprises the uncertainty xi, the solution of the model in the formula (15) belongs to an NP difficult problem, and the model can be converted into a semi-definite dual model which is easy to solve through Lagrange dual.
3.2 Dual conversion and piecewise linearization processing
When the CHP electrical load factor is fixed, the thermoelectric ratio is not a definite value by the uncertain quantity xi, the thermoelectric ratio in the thermal power balance type (9) needs to be converted into the sum of the original PLP function expression and the uncertain quantity xi, and the method specifically comprises the following steps:
Figure BDA0003275023880000093
due to the addition of the uncertain parameters, equation (9) is relaxed as:
Figure BDA0003275023880000094
since the constraint of the above equation is satisfied for all possible xi ∈ U, the above equation can be equivalent to:
Figure BDA0003275023880000101
containing an indeterminate amount for processing formula (19)
Figure BDA0003275023880000102
Term, construct the corresponding lagrangian function as:
Figure BDA0003275023880000103
in the formula, gammat、δtThe dual variables are respectively corresponding to the upper and lower limits of the uncertain variable xi, and have no practical significance gamma, and delta is respectively a dual variable vector (the vector length is the scheduling period T) corresponding to the upper and lower limits of the uncertain variable xi.
The above equation is derived, making the derivative equal to zero.
Figure BDA0003275023880000104
The visible derivative function is a constant independent of xi, so when xi is zero vector, the original Lagrange function takes the minimum value, i.e.
Figure BDA0003275023880000105
According to the strong dual theorem, the method can be obtained,
Figure BDA0003275023880000106
thus, equation (19) is equivalent to:
Figure BDA0003275023880000107
the PLP function of the energy conversion device gives the model an intractable non-linear term. Therefore, in the case of the formula (19)
Figure BDA0003275023880000108
And in formula (7)
Figure BDA0003275023880000109
And processing the nonlinear terms by adopting a piecewise linearization method. By linearization
Figure BDA00032750238800001010
For example, it is divided into d segments,
Figure BDA00032750238800001011
the following constraints are attached:
Figure BDA0003275023880000111
in the formula, mkIs the slope of each segment; a iskIs the longitudinal axis intercept of each segment; omegakIs a continuous variable introduced; bkIs the kth segmentation point; bk-1Is the k-1 segment point; v. ofkIs a variable from 0 to 1, v k1 represents PCHP,tAt the k interval [ b ]k-1,bk]In (1).
And finally converting the robust optimization model of the formula (15) into the following optimization model through Lagrange duality and piecewise linearization:
Figure BDA0003275023880000112
performing segmentation processing on other nonlinear terms according to a method of an equation (25) to an equation (26) to obtain a robust optimization scheduling model; the model belongs to a mixed integer semi-definite programming model, a gurobi solver can be called in Matlab software to directly solve, an optimal scheduling scheme is obtained after solving, and robust optimization scheduling of the regional comprehensive energy system is completed.
According to the algorithm, the formula (27) is optimized, and an optimization model is obtained as follows:
Figure BDA0003275023880000121
example 2
Taking a typical RIES as an example, the area takes electricity and natural gas as energy input, and the electricity and the heat load are simultaneously supplied through an energy conversion device and a transmission network. In the regionThe predicted electric load, thermal load and time-of-use electricity price on a certain day are respectively shown in fig. 3 and fig. 4. Natural gas value 2 yuan/m3,PELE,maxThe detailed parameters of each device in 1MW, rees and the fitting parameters of the device efficiency PLP function are shown in table 1 and table 2, respectively. For the piecewise linear representation of the fitting function, considering the balance between linear approximation precision and calculation efficiency, the number d of the segments is generally between 30 and 70, so that d is taken to be 50.
Setting an uncertain parameter
Figure BDA0003275023880000131
The RIES scheduling problem when r is 0 is a deterministic optimization problem, i.e. the device PLP function fitting error is not considered. And setting the RIES optimization scheduling scheme as a reference scheme when r is 0, wherein the larger the r value is, the larger the uncertainty range of the fitting error of the PLP function of the equipment is shown to be, and the higher the tolerance of the RIES optimization scheduling to the error disturbance is.
TABLE 1 energy conversion device parameters
Figure BDA0003275023880000132
TABLE 2 energy storage device parameters
Figure BDA0003275023880000133
After the simulation calculation, the influence of different values of the uncertain parameter r on the scheduling cost is analyzed, and the simulation result is shown in fig. 5.
As can be seen from fig. 5, in the reference scenario, the operation cost of the ires is the lowest (13050.4 yuan), and as r increases, the cost also increases, i.e., the operation cost is positively correlated with the uncertain parameter r. The running cost of the robust scheduling scheme when r is 0.02 is increased by 0.63% compared with the reference scheme cost. This is because, after considering the error of the PLP function fitting, the RIES needs to increase the cost of energy input to reduce the sensitivity to random fitting errors and improve the robustness of the scheduling scheme.
The comparison results of the thermal power and electric power scheduling results of the reference scheme (r is 0) and the robust scheme (r is 0.02) are shown in fig. 6 and 7, respectively.
As can be seen from fig. 6, the robust scheme has a large difference from the reference scheme, and the output thermal power of the robust scheme GB in each period is slightly higher than the corresponding value of the reference scheme, because the robust scheme makes an optimal decision in the worst scenario where the negative thermal bias is the most severe (i.e., the CHP actual thermal output is the lowest) in most periods, and in this scenario, the GB thermal output needs to be increased to prevent the ches actual thermal output from being insufficient to cause the RIES thermal load loss.
Fig. 7 compares the electric power scheduling situations of the two schemes, and it can be seen that the overall difference of the electric power scheduling of the two schemes is not great, but the electric power output by the CHP under the 6 th-7 th period robust scheme is 16kW less than that of the reference scheme, and the corresponding power purchased from the power grid is 16kW more than that of the reference scheme. This is mainly because the electricity load is low valley in the period of 6-7, CHP is in the low load state, and it can be known from the formula (2) that CHP has low power generation efficiency at low electricity load rate and high power generation cost, while fig. 6 shows that GB is in the high heat load rate state at this period and high heat generation efficiency. The total energy cost by purchasing electricity at valley price from the grid, supplying electricity, respectively thermal load using GB high efficiency heating is lower during 6-7 hours compared to cogeneration with CHP, so that the power purchased from the grid during this hour is increased under a robust scheme.
Although the baseline solution runs at a slightly lower cost than the robust solution, there is a greater risk of heat load loss due to uncertainty in the fitting error not being considered, and the present embodiment employs the heat load loss probability to quantify this risk. Note that the robust scheme with r being 0.02 is only a scheduling scheme in a scene with the worst fitting error value, and actually, the value of the uncertainty ξ in each time period is a random value in the box set range. Thus, to compare the risk of heat load loss under the baseline scheme and the robust scheme (r ═ 0.02), in ξ [ -0.02,0.02]Under the constraint of (2), Monte Carlo simulation is used for randomly generating M (1000 scenes), and the actual CHP thermal power under different xi values of each scene is calculated and substituted into a scheduling scheme. For a reference or robust scheduling scheme, if the thermal energy input of a certain scene t period cannot meet the thermal load requirement, the thermal load requirement is metThere is a loss of thermal load under the scenario. The number of scenes with heat load loss is recorded as V, and the heat load loss probability R is obtained by calculating 100% V/Mh. The two scenario heat load loss probabilities are shown in table 3.
TABLE 3 probability of heat load loss for two schemes
Figure BDA0003275023880000151
As can be seen from Table 3, R for the robust schemehMuch smaller than the baseline scheme, i.e. the risk of heat load loss for the robust scheme is greatly reduced compared to the baseline scheme. This is because the robust scheme makes a decision to satisfy the original constraint in the scenario where the deviation of the thermoelectric ratio is the most severe, that is, it is ensured that the obtained robust scheduling scheme can satisfy the thermal load demand for any different fitting errors.
The invention discloses and provides a robust optimization scheduling method of a regional comprehensive energy system, which adopts robust scheduling after considering the fitting error of an efficiency PLP function of energy conversion equipment, thereby reducing the risk of load loss; the robust scheduling model based on the fitting error of the PLP function of the energy conversion equipment efficiency can be used for not only counting the influence of single equipment, but also expanding and calculating the uncertainty of the fitting error of the PLP function of the energy conversion equipment efficiency, and further considering the coupling and correlation among the uncertain errors, so that the RIES optimal scheduling scheme is more practical; and the system robustness and reliability can be further improved by adopting the function fitting error and applying the function fitting error to the optimization model.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. The device disclosed by the embodiment corresponds to the method disclosed by the embodiment, so that the description is simple, and the relevant points can be referred to the method part for description.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (8)

1. A robust optimization scheduling method for a regional integrated energy system is characterized by comprising the following specific steps:
constructing a partial load performance function of each energy conversion equipment efficiency according to the regional comprehensive energy system framework;
constructing a regional comprehensive energy system model according to the regional comprehensive energy system framework and the partial load performance function of each energy conversion equipment efficiency;
constructing a constraint condition of the regional comprehensive energy system for the regional comprehensive energy system model;
constructing an optimized scheduling model according to a partial load performance function of each energy conversion equipment efficiency, a regional comprehensive energy system model and a constraint condition of the regional comprehensive energy system;
introducing fitting errors of partial load performance functions of the efficiency of each energy conversion device into an optimized scheduling model to construct a robust optimized scheduling model;
solving an optimal scheduling scheme of the regional comprehensive energy system by the robust optimization scheduling model;
and performing optimized dispatching on the regional comprehensive energy system according to the optimal dispatching scheme.
2. The robust optimal scheduling method for the regional integrated energy system according to claim 1, wherein the expression of the partial load performance function of the energy conversion equipment efficiency is as follows:
Figure FDA0003275023870000011
in the formula,
Figure FDA0003275023870000012
To the efficiency of the plant; and N is the load rate of the equipment, namely the ratio of the current output to the rated output of the equipment.
3. The robust optimized dispatching method for regional integrated energy systems according to claim 2,
the PLP of the generating efficiency of the CHP unit is as follows:
Figure FDA0003275023870000013
in the formula, kCHP,nFitting coefficients of PLP which are the generating efficiency of the CHP unit;
thermoelectric ratio of CHP
Figure FDA0003275023870000014
The relationship with the electrical load rate is:
Figure FDA0003275023870000015
in the formula, kλ,nFitting coefficient of PLP as thermoelectric ratio of CHP;
the fitting coefficients for PLP for thermal efficiency of GB are:
Figure FDA0003275023870000021
in the formula, kGB,nAre fitting coefficients.
4. The robust optimal scheduling method for regional integrated energy system according to claim 1, wherein the regional integrated energy system optimization model is established based on the regional integrated energy system framework with the goal of lowest daily operation cost of the regional integrated energy system.
5. The robust optimized dispatching method for regional integrated energy systems according to claim 1, wherein the constraints comprise power balance constraints, energy storage device-related constraints, and energy storage device-related constraints.
6. The robust optimal scheduling method for the regional integrated energy system according to claim 1, wherein the robust optimal scheduling model is constructed by the following steps:
describing a fitting error and an error range of a partial load performance function of the equipment efficiency by adopting an uncertain set;
introducing the uncertain set into an optimized scheduling model to form a compact form of a robust optimized scheduling model;
and converting the compact form of the robust optimized scheduling model into a mixed integer semi-positive definite programming model through Lagrange duality and piecewise linearization, and defining the mixed integer semi-positive definite programming model as the robust optimized scheduling model.
7. The robust optimized dispatching method for regional integrated energy systems according to claim 6, wherein the uncertainty set is a box uncertainty set.
8. The robust optimized dispatching method for regional integrated energy systems according to claim 1, further comprising verifying a heat load loss rate of the robust optimized dispatching model by using a Monte Carlo method under a constraint condition.
CN202111115252.5A 2021-09-23 2021-09-23 Robust optimization scheduling method for regional comprehensive energy system Pending CN113793052A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111115252.5A CN113793052A (en) 2021-09-23 2021-09-23 Robust optimization scheduling method for regional comprehensive energy system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111115252.5A CN113793052A (en) 2021-09-23 2021-09-23 Robust optimization scheduling method for regional comprehensive energy system

Publications (1)

Publication Number Publication Date
CN113793052A true CN113793052A (en) 2021-12-14

Family

ID=79184295

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111115252.5A Pending CN113793052A (en) 2021-09-23 2021-09-23 Robust optimization scheduling method for regional comprehensive energy system

Country Status (1)

Country Link
CN (1) CN113793052A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115392052A (en) * 2022-09-22 2022-11-25 华北电力大学(保定) Comprehensive energy system interval optimization method considering variable working condition operation characteristics of equipment

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106097154A (en) * 2016-07-12 2016-11-09 广东电网有限责任公司电力科学研究院 The dispatching method of a kind of industrial park energy centre and system
CN107704966A (en) * 2017-10-17 2018-02-16 华南理工大学 A kind of Energy Load forecasting system and method based on weather big data
CN110110897A (en) * 2019-04-11 2019-08-09 华北电力大学 A kind of integrated energy system optimization method considering different storage energy operation strategies
CN111210054A (en) * 2019-12-22 2020-05-29 上海电力大学 Micro-energy network optimization scheduling method considering direct load control uncertainty
CN111697625A (en) * 2020-06-29 2020-09-22 中国人民解放军国防科技大学 Island micro-grid opportunity constraint energy scheduling method based on distributed robust optimization
CN112036040A (en) * 2020-09-02 2020-12-04 国网上海市电力公司 Random scheduling method for regional comprehensive energy system
CN112465263A (en) * 2020-12-11 2021-03-09 国网宁夏电力有限公司经济技术研究院 Comprehensive energy operation optimization method suitable for multiple scenes
CN112836882A (en) * 2021-02-08 2021-05-25 杭州意能电力技术有限公司 Regional comprehensive energy system operation optimization method considering equipment load rate change

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106097154A (en) * 2016-07-12 2016-11-09 广东电网有限责任公司电力科学研究院 The dispatching method of a kind of industrial park energy centre and system
CN107704966A (en) * 2017-10-17 2018-02-16 华南理工大学 A kind of Energy Load forecasting system and method based on weather big data
CN110110897A (en) * 2019-04-11 2019-08-09 华北电力大学 A kind of integrated energy system optimization method considering different storage energy operation strategies
CN111210054A (en) * 2019-12-22 2020-05-29 上海电力大学 Micro-energy network optimization scheduling method considering direct load control uncertainty
CN111697625A (en) * 2020-06-29 2020-09-22 中国人民解放军国防科技大学 Island micro-grid opportunity constraint energy scheduling method based on distributed robust optimization
CN112036040A (en) * 2020-09-02 2020-12-04 国网上海市电力公司 Random scheduling method for regional comprehensive energy system
CN112465263A (en) * 2020-12-11 2021-03-09 国网宁夏电力有限公司经济技术研究院 Comprehensive energy operation optimization method suitable for multiple scenes
CN112836882A (en) * 2021-02-08 2021-05-25 杭州意能电力技术有限公司 Regional comprehensive energy system operation optimization method considering equipment load rate change

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘科明等: "考虑设备部分负荷性能函数拟合误差的区域综合能源系统鲁棒优化调度", 《电力系统及其自动化学报》, vol. 34, no. 5, pages 11 - 17 *
周任军等: "碳捕集中考虑碳价概率分布拟合抽样误差的超分位数–鲁棒优化方法", 《中国电机工程学报》, vol. 40, no. 4, 20 February 2020 (2020-02-20), pages 1130 - 1136 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115392052A (en) * 2022-09-22 2022-11-25 华北电力大学(保定) Comprehensive energy system interval optimization method considering variable working condition operation characteristics of equipment

Similar Documents

Publication Publication Date Title
Yuan et al. Real-time pricing for smart grid with multi-energy microgrids and uncertain loads: a bilevel programming method
CN111815025A (en) Flexible optimization scheduling method for comprehensive energy system considering uncertainty of wind, light and load
CN110689189A (en) Combined cooling heating and power supply and demand balance optimization scheduling method considering energy supply side and demand side
Xia et al. Two‐stage robust optimisation of user‐side cloud energy storage configuration considering load fluctuation and energy storage loss
Abdolahi et al. Chance-constrained CAES and DRP scheduling to maximize wind power harvesting in congested transmission systems considering operational flexibility
Lei et al. Optimal allocation of a VRB energy storage system for wind power applications considering the dynamic efficiency and life of VRB in active distribution networks
Zeinal‐Kheiri et al. Real‐time energy management in a microgrid with renewable generation, energy storages, flexible loads and combined heat and power units using Lyapunov optimisation
CN111030094B (en) Method for evaluating reliability of energy supply of multi-energy complementary park
CN113162025A (en) Demand response-containing distributed low-carbon economic dispatching method for electrical interconnection network
Jithendranath et al. Stochastic planning of islanded microgrids with uncertain multi‐energy demands and renewable generations
Dimitroulis et al. Residential energy management system utilizing fuzzy based decision-making
CN111860966A (en) Energy storage-containing comprehensive energy system scheduling method based on fuzzy correlation opportunity planning
Li et al. Energy supply reliability assessment of the integrated energy system considering complementary and optimal operation during failure
Li et al. Research on optimal self‐scheduling horizon for the wind power and large‐scale CAES combined system
Thang Optimal sizing of distributed energy resources and battery energy storage system in planning of islanded micro-grids based on life cycle cost
Azizi et al. Daily operation of multi‐energy systems based on stochastic optimization considering prediction of renewable energy generation
CN113793052A (en) Robust optimization scheduling method for regional comprehensive energy system
Wang et al. Multi‐objectives combined electric heating dispatch model of wind power accommodation with heat storage device
Zhang et al. A two-stage robust operation approach for combined cooling, heat and power systems
Zhai et al. Low carbon economic dispatch of regional integrated energy system considering load uncertainty
CN114169941A (en) Comprehensive energy system low-carbon scheduling method considering carbon value uncertainty
Li et al. Multiobjective Optimization Model considering Demand Response and Uncertainty of Generation Side of Microgrid
Die et al. Optimal approach in the EMS of microgrids based on online forecasts
Salazar et al. Optimum energy management of islanded nanogrids through nonlinear stochastic dynamic programming
Qi et al. Coordinated optimization of integrated energy system considering demand response and energy storage

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