Two-Stage Robust Economic Dispatch of Regional Integrated Energy System Considering Source-Load Uncertainty Based on Carbon Neutral Vision

: A regional integrated energy system is an important carrier of the energy Internet. It is a major challenge for the operation of a regional integrated energy system to deal with the uncertainty of distributed energy and multiple loads by using the coupling characteristics of equipment in a regional integrated energy system. In this paper, a two-stage robust economic dispatch model of a regional integrated energy system is proposed considering the source-load uncertainty. Firstly, the basic architecture of the regional integrated energy system is introduced. Based on the extreme scenario of uncertain power supply and load, the uncertainty set was established, the two-stage robust optimization model of regional integrated energy system was constructed and the column-and-constraint generation algorithm was used to solve the model. The effectiveness of the two-stage robust optimization model in improving the economy and robustness of the system was analyzed.


Introduction
China proposes to achieve the development goal of carbon neutrality by 2060 to contribute to the sustainable economic development of the Asia-Pacific region [1]. To achieve this, renewable energy penetration must be increased. However, with the increase in the penetration rate of renewable energy, large-scale renewable energy is difficult to absorb, and the phenomenon of "abandoning wind and light" still exists. With the increase in energy consumption, the failure of the power grid often endangers other energy systems. The problem of independent planning and operation of renewable energy and various energy systems is becoming increasingly prominent. Faced with this challenge, the concept of the energy Internet has emerged and become the focus of academic and industrial circles [2]. The outstanding feature of the energy Internet is the combination of new energy technology and information technology, and the complementary coupling of cold, heat and electricity energy. According to the scale, integrated energy systems can be divided into trans-regional level, regional level and user level [3]. A regional integrated energy system (RIES) is the manifestation of regional distribution and system characteristics of an integrated energy system. The RIES is involved in energy production, conversion, transportation, storage and use; the depth of fusion power, intermittent renewable energy sources such as wind, solar, natural gas, biomass energy and environmental energy, such as thermal energy through various types of complementary energy between coupling and optimization, and renewable energy in order to meet the user's cold, heat, electricity and other energy demand. At the same time, efficient and clean utilization of energy can be achieved [4][5][6].
How to utilize the multi-energy coupling of RIESs and the flexibility of multi-type energy storage to cope with the uncertainty of electricity, gas and heat loads and renewable energy is a frequently studied topic in the research of integrated energy systems, as is research on regional integrated energy [7][8][9][10]. Ref. [7] discussed the transformation of regional integrated energy with the control of carbon emission reduction, and took Dalian City, China, as an example for analysis. Ref. [9] studied the market transaction and management of regional integrated energy considering the application of high proportion of renewable energy.
Ref. [11] considers that the use of battery energy storage systems (BESSs) as backup power affects the operating costs of regional integrated energy systems (RIESs) in different situations and describes a regional integrated energy system that includes wind turbines, photovoltaic generators, gas turbines and battery energy storage systems. Refs. [8,12] proposes a multi-objective stochastic programming method based on the generation method of multi-dimensional associated scenario sets for RIES expansion plans. The scenario generation method considers the characteristics, timing, autocorrelation and cross-correlation of renewable energy and multi-energy loads. Considering the importance of the energy pipeline, the pipeline risk index of energy network expansion planning is defined. Meanwhile, a multi-objective stochastic energy network-based opportunity-constrained programming model is proposed to reduce investment costs and pipeline risks. Ref. [13] extended the concept of demand response (DR) to the RIES and proposed an optimized operation model of the RIES considering the DR mechanism on energy prices. Based on the DR modeling of the RIES, the operation optimization model with environmental benefits, economic benefits and energy supply reliability as objective functions was established in detail for the first time. In the configuration of the rated capacity and power of various energy storage devices in the RIES [14], off-grid and grid-connected operation modes are proposed and a configuration optimization model is established. The RIES is divided into four parts: power supply, energy conversion, energy storage and load. Based on the concept of the energy hub, the four parts are modeled separately.
However, it is difficult to obtain the probability distribution of random parameters in practical application, which limits the application of the above reference algorithm. The global optimality of solution efficiency and solution cannot be satisfied simultaneously because of the existence of non-convex constraints in chance-constrained programming. The computational efficiency and accuracy of the probabilistic scenario method are easily affected by the number of generated scenarios and the method of scene generation and reduction.
Robust optimization does not require presupposition of the distribution followed by uncertain parameters. After the fluctuation range of the uncertain parameters is determined in advance, when the value of the uncertain parameters falls in the given uncertain set, the constraints of the robust optimization model must be satisfied, and the deterministic feasible solution [15][16][17] can be obtained. Robust optimization has developed into an important and very common tool to deal with uncertain optimization problems and has been widely used in power system optimization scheduling [18][19][20][21][22][23][24]. In view of the existence of various uncertainties, Ref. [23] proposed an uncertainty quantification method for energy management. In the rescheduling stage, decision variables representing the charge and discharge state of energy storage systems were considered, while the existing model places them in the pre-scheduling stage to improve the operational flexibility of energy storage systems in energy management.
Two-stage robust optimization (TRO) can ensure the normal operation of the system in the most extreme case where only the fluctuation range information of uncertain parameters is obtained. At the same time, TRO can adjust the conservatism of the system through dynamic iteration strategy. This addresses the deficiency of the traditional robust optimization method.
The contributions of this study can be summarized as follows.
(1) The uncertainties of wind power, photovoltaic output and multiple types of loads were considered, and a two-stage robust and optimized economic scheduling model was established with the goal of minimizing operation cost and regulation cost with dayahead scheduling. First, the model was built taking into account the electrical, gas and thermal network characteristics of the RIES as well as energy conversion equipment, energy storage equipment and various composite characteristics. Considering the extreme cases of wind and photovoltaic power and load forecasting errors in dayahead scheduling, the TRO model was constructed. (2) Combined with the strong duality theory, the column-and-constraint generation algorithm was used to transform the model into a mixed integer linear programming problem. (3) A regional integrated energy system in summer was taken as an example to verify the effectiveness of the model. The difference between the two-stage robust optimization model with the deterministic scheduling method and the worst case with different uncertain budget values was compared and analyzed. The results show that the twostage robust economic optimization model can significantly improve the operating economy of the system.
The remainder of this paper is organized as follows. In Section 2, the structural framework of the regional integrated energy system is introduced. Section 3 mainly lists the output model of the main equipment of the regional integrated energy system. In Section 4, the construction of a two-stage robust optimal scheduling model is detailed, the objective function is proposed to minimize the system adjustment cost and the constraints of the system are listed. Section 5 proposes column-and-constraint generation to solve the model. In Section 6, a case study is outlined to prove the validity and economy of the proposed model. Concluding remarks are given in Section 7.

Regional Integrated Energy System Structure
In structure, a RIES can be regarded as an energy unit with multiple inputs, outputs, conversions and storage devices. Based on common forms of regional integrated energy systems, a RIES containing energy coupling units, distributed power sources and energy storage devices was constructed. The structure of the RIES is shown in Figure 1. The energy coupling unit includes the gas turbine cogeneration unit, gas boiler, electric boiler, electric refrigeration equipment and absorption refrigeration equipment. Distributed power includes distributed wind power and distributed photovoltaic power generation units. Energy storage equipment includes electric energy storage and thermal energy storage. The RIES purchases online energy from the upstream external power grid and heat network and natural gas through the regional electricity, heat and gas vendors. The energy is reasonably distributed through conversion and storage devices to meet various needs of users.  Figure 1. Schematic diagram of regional integrated energy system.

GAS Turbine Cogeneration Unit
GAS Turbine Cogeneration Unit is an important part of regional Integrated Energy system.
where P GT t and H GT t are the electrical power and thermal power output by the gas turbine at time period t, respectively; G GT t is the amount of natural gas consumed by gas turbine in time period t; η GT,p is the power generation efficiency of gas turbine; η GT,h is the waste heat recovery efficiency of gas turbine; α GT,p t is the scheduling factor of gas turbine in time period t; P GT,min and P GT,max are the minimum and maximum output power of gas turbine, respectively; H GT,min and H GT,max are the minimum and maximum thermal output power of gas turbine, respectively.

Gas Boiler
When the heat energy recovered by waste heat recovery device cannot meet the demand of the heat load, gas-fired boiler can be used to supply energy to part of the heat load. The thermal power output constraint of the gas-fired boiler can be expressed as follows.
where H GB t is the thermal power output by gas boiler in time period t; G GB t is the amount of natural gas consumed by gas boiler in time period t and η GB is the heating efficiency of gas boiler. H GB,min and H GB,max are the minimum and maximum thermal power of the gas-fired boiler, respectively.

Electric Boiler
The power constraints of electric heating boilers are shown below.
where H EB t is the thermal power output by the electric boiler in time period t; P EB t is the amount of natural gas consumed by the electric boiler in time period t and η EB is the power generation efficiency of the electric boiler; H GB,min and H GB,max are the minimum and maximum thermal power of electric boiler, respectively.

Absorption Refrigeration
Absorption chillers use thermal energy to drive the refrigeration cycle and consume thermal energy for refrigeration. The power constraints of an absorption chiller are described as follows.
where C AC t is the output power of the absorption chiller in time period t, H AC t is the thermal power absorbed by the absorption chiller in time period t, η AC is the efficiency of heat exchange device, C AC,min and C AC,max are the maximum and minimum output power of the absorption chiller, respectively.

Electric Refrigeration
The physical model of an electric refrigerator is described as follows: where C EC t is the output power of the electric refrigerator in time period t, P EC t is the thermal power absorbed by the electric refrigerator in time period t, η EC is the conversion efficiency of the electric refrigerator, C EC,min and C EC,max are the maximum and minimum output power of the electric refrigerator, respectively.

Heat Storage Equipment
The heat storage device can store heat and release heat when the system heat is insufficient. It is assumed that there is no energy dissipation over time in the heat storage device, and the energy storage model of the heat storage device can be expressed as: where P ES,ch t and P ES,dis t are respectively the charge and discharge power of electric energy storage equipment in time period t; P ES,ch,min and P ES,ch,max are the minimum and maximum charging power of electric energy storage equipment, respectively; P ES,dis,min and P ES,dis,max are the minimum and maximum discharge power of electric energy storage equipment, respectively; u ES,ch t and u ES,dis t are charge and discharge state variables of the electric energy storage device at time period t, respectively; S E t is the charged state of the electric energy storage device at time period t; S E,min and S E,max are the minimum and maximum state of charge of electric energy storage equipment, respectively. η ES,ch and η ES,dis respectively refer to the charging and discharging efficiency of an electric energy storage device; Cap ES is the rated capacity of an electric energy storage device.

Modeling of Uncertainties
Considering the actual output of distributed power sources such as distributed wind power and distributed photovoltaic power, and that all kinds of user loads such as cold, heat, gas and electricity fluctuate around the predicted output, in this study, the uncertainty sets of distributed power supply and multi-class load were established to describe the related uncertainties, and a two-stage robust optimal scheduling model of integrated energy was established.
This study analyzed the historical data of wind power and photovoltaic output. According to the historical data, a range was selected to cover all wind power or photovoltaic output scenarios in a given period of time. Thus, the uncertain performance of wind and photovoltaic power can be expressed as follows with the uncertainty set constrained by a certain radius: where U x refers to the uncertain set of wind and photovoltaic power and user load such as cold, heat, gas and electricity; x ∈ {W, PV, HL, CL, GL, PL} is energy and load type. P x t represents uncertain variables in the system; P x, f t represents the predicted value of distributed power and load. P x t and P x t represent the upper and lower limits of fluctuation deviation of distributed power supply and load, respectively; z x,1+ t and z x,1− t are auxiliary 0-1 variables describing the distributed power supply and the actual load value at the upper or lower boundary. Γ x is the uncertain budget value, which is used to adjust the conservatism of the optimal solution; the larger its value, the more conservative the final scheduling scheme.

Regional Integrated Energy System Two-Stage Robust Optimal Scheduling Model Considering Bilateral Uncertainties
Based on the RIES model, a two-stage robust optimization RIES model was established. The first phase aims to determine the optimal scheduling scheme for wind, photovoltaic and multi-energy loads in the day-ahead forecast scenario. Based on the optimization results of the first stage, the output of each unit is adjusted again in the second stage to cope with the "most extreme" scenario with the largest prediction error with the day-ahead scheduling scheme, and the results are fed back to the first stage. Through the iterations of the first and second stages, all extreme scenarios can satisfy the given constraints, and the corresponding day-ahead scheduling strategy can be obtained.

The Objective Function
A two-stage robust optimization RIES model is presented in this paper, and the objective function is shown in Equations (29)-(32). The main objective of the first stage is to minimize the total energy consumption cost and operation and maintenance cost of the RIES with the given day-ahead forecast, and the max-min objective of the second stage ensures the lowest system adjustment cost in the uncertain scenario of intra-day scheduling.
where ∆• represents the deviation of each scheduling value in intra-day scheduling relative to intra-day scheduling; b • t represents the penalty cost of unit deviation of each conversion equipment and purchased energy.

The Constraints
Constraints are divided into first-stage constraints and second-stage constraints. The first-stage constraints include energy network and coupling device constraints in the dayahead scheduling process. The second stage constraints include the adjustment constraints of electricity, gas, heat networks and coupling equipment with the uncertain parameters of the day.

First Stage Constraints
The constraints in the first stage include Equations (1)-(24), discarding wind discard light constraint (33)-(34) and power balance constraints in the prediction scenario, as shown in (35)

Second Stage Constraints
Because the constraint of the first stage does not consider the influence of the uncertain scenario in the real-time operation stage, in the second stage, distributed power supply and multiple types of loads are introduced into the constraints, and rescheduling constraints in uncertain scenarios are considered in the day-ahead scheduling of the RIES. The secondstage constraints are shown below.

Column-and-Constraint Generation Algorithm
The two-stage robust optimization model established in this study is a two-stage, three-layer large-scale min-max-min problem with robust constraints, which is difficult to solve simultaneously. Therefore, the column-and-constraint generation (CCG) algorithm was used to transform the model into a two-level optimization model containing main problem (MP) and subproblem (SP), and iterative solutions were performed.
The bilayer robust optimization model constructed in this study is expressed in compact form as follows: Obj = min x c x x + max u min y c y y (63) In the first stage, the outer layer is minimized, and the optimization variable is x; in the second stage, the inner layer is minimum-maximized, and the optimization variables are u and y; A 1 , A 2 , A 3 , G 1 , G 2 , D 2 , D 3 , J 2 are coefficient matrices; c x and c y are the coefficient column vectors of the first and second stages, respectively; b 1 , b 2 , b 3 , h 1 , h 2 are constant sequence vectors.
The CCG algorithm divides the original problem into main problem and subproblem. When solving the master problem, the constraints and variables related to the subproblem are always introduced to solve MP and SP iteratively. When solving the proposed problem, it is necessary to initialize the uncertain variable u as the initial worst scenario. The value of the uncertain variable in the worst case is obtained by solving the subproblem, and the main problem is solved according to the worst case obtained by the subproblem.
MP is a relaxation of the original problem, providing a lower bound on the optimal solution. SP obtains the value of the uncertain variable in the worst case and gives the upper bound of the optimal solution. The main problem can be expressed as follows: where y k is the scheduling mechanism after the k-th iteration; u k 1 and u k 2 represent the worst scenario obtained by the k-th iteration.
By solving the subproblem, the uncertain variables in the worst scenario are passed to MP. The dual problem of the subproblem can be expressed as follows: The subproblem can be transformed into a mixed integer linear programming problem by introducing auxiliary variables and the large M method. The CCG algorithm is used to solve the problem through the above transformation. The solution process is shown in Figure 2.

Case
Taking a typical summer day of an integrated energy system in Northern China as an example, the RIES included fans, photovoltaics, cogeneration units, electric boilers, electric refrigeration, gas boilers, absorption refrigeration and energy storage devices. The scheduling period is 24 h. Based on historical data, the fluctuation deviation of wind and photovoltaic power output is 15%, and the fluctuation deviation of electric load, cooling load, heat load and natural gas load is 10%. Figures 3 and 4 shows the deviation ranges of photovoltaic and wind power and various loads.

Case
Taking a typical summer day of an integrated energy system in Northern China as an example, the RIES included fans, photovoltaics, cogeneration units, electric boilers, electric refrigeration, gas boilers, absorption refrigeration and energy storage devices. The scheduling period is 24 h. Based on historical data, the fluctuation deviation of wind and photovoltaic power output is 15%, and the fluctuation deviation of electric load, cooling load, heat load and natural gas load is 10%. Figures 3 and 4 show the deviation ranges of photovoltaic and wind power and various loads.

Case
Taking a typical summer day of an integrated energy system in Northern China as an example, the RIES included fans, photovoltaics, cogeneration units, electric boilers, electric refrigeration, gas boilers, absorption refrigeration and energy storage devices. The scheduling period is 24 h. Based on historical data, the fluctuation deviation of wind and photovoltaic power output is 15%, and the fluctuation deviation of electric load, cooling load, heat load and natural gas load is 10%. Figures 3 and 4 shows the deviation ranges of photovoltaic and wind power and various loads.

Optimization Results of Regional Integrated Energy System
The scheduling results of equipment in the RIES are shown in Figures 3-8. During 1-6 and 23-24 h, the unit electricity price is lower, and the RIES mainly purchases power from the grid. At this time, the gas turbine runs at a lower power, and the electric boiler and electric refrigerator provide a small amount of heat energy and cold energy, respectively. At 7-12 h the electricity price rises, while the simultaneous electric load and cooling load gradually rise. At this time, the gas turbine power increases, while due to the increase in electricity price, the RIES purchases a lower proportion of electricity from the grid and flexibly relies on the discharge of electric energy storage equipment to balance the power load demand of the system. As a cogeneration unit, the heat energy provided by gas turbine rises with the increase in electric power output. At this time, the cooling load of the system is mainly generated by the absorption chiller absorbing heat energy, and a small amount of cooling load is provided by electric refrigeration equipment. From 12 to 17 h, the price of electricity decreases when it reaches a high peak value, the cooling load and heat load rise further and the electric load reaches the peak value at 16 h. Therefore, the proportion of purchased electricity and purchased electricity both rise, and the electric energy storage is charged to deal with the peak load in the late 18-22 h period. During 18-22 h, the price of electricity rises, while the heat and cold loads peak. At this time, the output of the gas turbine increases, and the thermal output of the gas turbine increases. The heat load at 19 and 20 h is completely provided by the gas turbine. The cooling load is mainly provided by the absorption chiller, and part of the cooling load is provided by the electric chiller. From 23 to 24 h, the electricity price falls to the lowest value, and the cold and heat loads decrease. Most of the electric loads are provided by the external purchasing power of the grid, and the proportion of electric refrigeration and heating increases.

Optimization Results of Regional Integrated Energy System
The scheduling results of equipment in the RIES are shown in Figures 3-6. During 1-6 and 23-24 h, the unit electricity price is lower, and the RIES mainly purchases power from the grid. At this time, the gas turbine runs at a lower power, and the electric boiler and electric refrigerator provide a small amount of heat energy and cold energy, respectively. At 7-12 h the electricity price rises, while the simultaneous electric load and cooling load gradually rise. At this time, the gas turbine power increases, while due to the increase in electricity price, the RIES purchases a lower proportion of electricity from the grid and flexibly relies on the discharge of electric energy storage equipment to balance the power load demand of the system. As a cogeneration unit, the heat energy provided by gas turbine rises with the increase in electric power output. At this time, the cooling load of the system is mainly generated by the absorption chiller absorbing heat energy, and a small amount of cooling load is provided by electric refrigeration equipment. From 12 to 17 h, the price of electricity decreases when it reaches a high peak value, the cooling load and heat load rise further and the electric load reaches the peak value at 16 h. Therefore, the proportion of purchased electricity and purchased electricity both rise, and the electric energy storage is charged to deal with the peak load in the late 18-22 h period. During 18-22 h, the price of electricity rises, while the heat and cold loads peak. At this time, the output of the gas turbine increases, and the thermal output of the gas turbine increases. The heat load at 19 and 20 h is completely provided by the gas turbine. The cooling load is mainly provided by the absorption chiller, and part of the cooling load is provided by the electric chiller. From 23 to 24 h, the electricity price falls to the lowest value, and the cold and heat loads decrease. Most of the electric loads are provided by the external purchasing power of the grid, and the proportion of electric refrigeration and heating increases.    In order to analyze the economy of the proposed method, the results of the two-stage robust optimal scheduling method were compared with those of the traditional deterministic scheduling method. Given that intra-day output deviation of wind and photovoltaic power is negative 15% and load deviation is positive 10%, in this scenario, the deterministic scheduling was compared with the two-stage scheduling conclusion proposed in this    In order to analyze the economy of the proposed method, the results of the two-stage robust optimal scheduling method were compared with those of the traditional deterministic scheduling method. Given that intra-day output deviation of wind and photovoltaic power is negative 15% and load deviation is positive 10%, in this scenario, the deterministic scheduling was compared with the two-stage scheduling conclusion proposed in this    In order to analyze the economy of the proposed method, the results of the two-stage robust optimal scheduling method were compared with those of the traditional deterministic scheduling method. Given that intra-day output deviation of wind and photovoltaic power is negative 15% and load deviation is positive 10%, in this scenario, the deterministic scheduling was compared with the two-stage scheduling conclusion proposed in this In order to analyze the economy of the proposed method, the results of the twostage robust optimal scheduling method were compared with those of the traditional deterministic scheduling method. Given that intra-day output deviation of wind and photovoltaic power is negative 15% and load deviation is positive 10%, in this scenario, the deterministic scheduling was compared with the two-stage scheduling conclusion proposed in this study, and day-ahead scheduling cost and intra-day scheduling energy purchase deviation cost were calculated as shown in Table 1. The purpose of the two-stage robust optimization constructed in this study is to minimize the total operating cost of the system in extreme scenarios. In this scenario, the day-ahead scheduling plan needs to satisfy the system constraints when the system uncertainty fluctuates the most so that the two-stage robust optimization day-ahead scheduling strategy is conservative. It can be seen from Table 1 that the day-ahead scheduling cost of two-stage robust optimization is higher than that of deterministic optimization method. In the diurnal scheduling stage, the difference between diurnal scheduling and real-time energy supply needs to be compensated by adjusting unit output. For the two-stage robust optimization, the conservative day-ahead scheduling makes it more flexible and the intraday adjustment cost is relatively low. Therefore, the total operating cost of the operation strategy system proposed by the two-stage robust optimization is lower than that of the deterministic optimization. Therefore, compared with the deterministic optimization, the scheduling strategy based on the two-stage robust optimization has stronger robustness and has stronger ability to resist the fluctuating risks of wind and photovoltaic power and multi-energy loads.

Scenario Analysis
In order to verify the capability of two-stage robust optimization in dealing with wind, light and load uncertainties, four scenarios were designed for comparative analysis.
Scenario 1 is the basic case without considering the uncertainties of wind, light and load, and the robustness coefficient is Γ G = Γ L = 0; Scenario 2 only considers the uncertainty of the power generation side, and the robustness coefficient is set as Γ G = 16, Γ L = 0; In scenario 3, only load-side uncertainty is considered, and the robustness coefficient is set as Γ G = 0, Γ L = 10; In scenario 4, both generator-side and load-side uncertainties are considered, and the robustness coefficient is set as Γ G = 16, Γ L = 10.
Day-ahead scheduling results, intra-day scheduling results and deviation coefficients in the four scenarios are shown in Table 2. When both source and load uncertainties are ignored, the day-ahead optimal operating cost is higher than that when load uncertainties are not considered. This because of the ignored load uncertainty and uncertain energy sources such as wind and light that depend on the forecast for scheduling, and there may be a positive deviation, increasing the system scheduling capacity. The system without scheduling at the present stage in response to increased possible load cost increases; therefore, scenario 2 has a scheduling stage cost comparable to scenario 1. Similarly, the pre-dispatch cost of scenario 4 decreases compared with scenario 3. The cost of day-ahead scheduling in scenario 3 is higher than that in scenario 1 because scenario 3 considers the load uncertainty in day-ahead scheduling, and the system needs to formulate strategies in day-ahead scheduling to deal with the deviation caused by load fluctuation. Therefore, the cost of day-ahead dispatching rises.

Conclusions
In this paper, a regional integrated energy system optimization method based on two-stage robust optimization is proposed, which takes into account the uncertainty of source-load on the east side of the RIES; establishes the combined optimization problem of energy conversion equipment, electricity storage and heat storage devices and gives the day-ahead scheduling and energy purchase scheme of the RIES. The validity of the method was verified by case analysis.
On the energy supply side, compared with previous studies that considered single energy, this study comprehensively considered the coupling characteristics of various energy sources, improved the absorption capacity of new energy and utilized the coupling characteristics of energy sources to suppress the fluctuation of the new energy grid-connected mode. On the load side, the uncertainty of the load side is added into the constraint, and the uncertainty of the load side is considered as well as the uncertainty of the energy supply side. In terms of model solving, it is difficult to obtain the probability distribution of random parameters in practical application, which limits the application of chance constrained programming, the probabilistic scenario method and the point estimation method. The global optimality of the solution efficiency and solution cannot be satisfied simultaneously because of the existence of non-convex constraints in chance-constrained programming. The computational efficiency and accuracy of the probabilistic scenario method are easily affected by the number of generated scenarios and the method of scene generation and reduction. In order to solve these problems, this paper proposes a column constraint algorithm to solve the model. The validity of the method was verified by case analysis.
The results show that: (1) The cost of day-purchase deviation obtained by the two-stage robust optimization method is lower than that obtained by the deterministic optimization method, while the traditional scheduling method does not consider the uncertain factors in advance, resulting in a large day-purchase deviation. The total cost of two-stage robust optimal scheduling is better than that of traditional deterministic scheduling methods. (2) Through the multi-energy coupling device and multi-energy coordination joint optimization, the proposed day-ahead scheduling strategy can effectively cope with the source-charge uncertainty in the day-ahead scheduling stage and enhance the robustness of the system. (3) By adjusting the uncertain budget value, RIES operators can flexibly adjust the scheduling conservatism of the RIES, which provides a reference for RIES operators to weigh the robustness and economy of the system.