Robust Optimization Model for Energy Purchase and Sale of Electric–Gas Interconnection System in Multi-Energy Market

With the increasing coupling of the power system and the natural gas system, the electric–gas interconnection system has become a typical form of comprehensive energy utilization. Through the energy conversion function of the coupling unit, the system can flexibly participate in the bidding for purchasing and selling energy in a power market and a natural gas market on the premise of meeting the internal demand of multiple loads. To solve the internal coordination and optimization problem and the external flexible bidding problem in the multi-energy market, this paper proposes a robust optimization model of energy purchase and sale for the electric–gas interconnection system in a multi-energy market. Firstly, the basic structure of the electric–gas interconnection system is introduced, and the steady-state model of energy flow in the system is built based on the energy hub model. Secondly, considering the uncertainty of energy prices and the output power of renewable energy units in the system, a bidding model for energy purchase and sale of the electric–gas interconnection system in multi-energy market based on the idea of robust optimization is constructed in the framework of the Nordic energy market. Finally, empirical analysis based on the actual data is carried out, and the results prove the validity and superiority of the model. In this paper, aiming at the uncertainty of energy price, a large number of scenes are generated by Latin hypercube sampling (LHS), and then a k-means algorithm is used to reduce the scenes, so as to simulate typical scenes. Aiming at the uncertainty of the output power of the renewable energy unit in the system, a cardinal uncertainty set is used to control deviation between the actual output power and predicted output power, so that the overall robustness of the model can be controlled. The proposed model can make decision-making independent of the accurate probability distribution of uncertainty factors, and is suitable for complex multi energy systems. Meanwhile, the model possesses excellent robustness, which can effectively reduce the risk of bidding loss in the process of energy purchase and sale.


Introduction
An integrated energy system deepens the coupling of the power system and other energy systems by utilizing energy conversion technology, thereby integrating multiple heterogeneous energy sources to achieve multi-energy complementarity and energy cascade utilization [1][2][3]. The participation of multiple energy sources makes an integrated energy system extremely flexible in participating in market transactions [4,5]. However, due to the differences in output power characteristics, market structure, price mechanism and price volatility of heterogeneous energy sources, the decision-makersʹ bidding decisions should not only consider the internal influences of output power characteristics and the coupling situation of multi-type heterogeneous energy sources, but also consider the external influence of the multi-dimensional market structure and the multi-level price system, which will bring enormous challenges for bidding decision-makers [6]. Therefore, it is urgent to carry out research on the bidding strategy of the integrated energy system in a multi-energy market, which will provide a scientific and reasonable decision-making basis for both the internal multi-energy collaborative operation management and the external multi energy market bidding.
The electric-gas interconnection system, mostly based on industrial parks, is an integrated energy system with broad development prospects [7,8]. The system is connected with a public power network and natural gas pipeline network externally, and can realize combined gas-heating-coolingpower supply through an energy conversion unit internally, whose internal energy conversion and utilization are relatively independent [9,10]. Considering that the electric-gas interconnection system involves the coupling and conversion of gas-heating-cooling-power energy sources, it is necessary to build the steady-state model of energy flow about multi-type heterogeneous energy in decisionmaking optimization [11][12][13][14]. An energy hub model is widely used in energy flow analysis of the multi-energy systems, especially in the research about planning and operation optimization of multienergy systems, which is defined as an input-output port model for describing the exchange and coupling relationship of energy, load and network under the multi-energy systems. So far, however, the general energy hub models seldom consider the influence of system uncertainties, which greatly reduce the practical application effect of the model.
In the interaction process with the external multi-energy market, the role of the electric-gas interconnection system is not only the energy producer, but also the energy consumer [15]. When its internal energy supply cannot meet demand, it will purchase energy from the external natural gas market or the power market, and when its internal energy supply exceeds the demand, it can sell the surplus to the natural gas market or the power market. The particularity of its positioning makes its bidding behavior different from that of the traditional generation side or the user side market [11][12][13][14][15]. Meanwhile, due to the fact that the electric-gas interconnection system can flexibly participate in the purchase and sale of energy in the multi-energy system, it is necessary to make a coordinated bidding decision in the multi-energy market. So far, there have been many research works on the bidding strategy optimization for both generation side or user side in the power system, and the bidding decision-making subjects mainly include single energy generator, joint generator of renewable energy and traditional energy, joint generator of renewable energy and energy conversion unit, virtual power plant considering demand side resources, electricity seller of multi-type load aggregation, large consumer, micro-grid subject of paralleling in the network or off-network flexibly, etc.. All the bidding decision-making subjects except the micro-grid only participate in market transactions as a single buyer or a single seller. Due to the limitations of the decision-making subjects and market background, there are few studies on the optimization of energy bidding decisionmaking which consider the duality of the multi energy system's identity in the market and the differences of the multi energy market.
When research into the bidding decision-making optimization method is carried out, according to the treatment of uncertain variables, optimization methods can be generally divided into two categories: stochastic optimization and robust optimization [16]. Stochastic optimization, depending on the probability distribution of uncertain variables, can transform uncertain optimization problems into deterministic optimization problems by describing the continuous or discrete probability distribution characteristics of uncertain variables. However, there are many problems in this method, such as the fact that it easily falls into dimension disaster, a high accuracy requirement for probability estimation, poor anti-noise ability and so on. Robust optimization is an optimization method that considers the influence of the worst-case on system decision-making. Compared with stochastic optimization, robust optimization focuses on the boundary condition of uncertain variables, so it is not necessary to obtain the accurate probability distribution of uncertain variables, and its conservatism can be controlled by adjusting the scale of uncertainty set [17]. A robust optimization method can effectively solve the problems such as high sensitivity of the model to the parameters and large influence of parameter estimation error on the model, which is widely used in the power system unit commitment [18] and economic dispatch [19][20][21], but seldom used in the multi-energy system.
Based on the above background, this paper proposes a robust optimization model of energy purchase and sale for the electric-gas interconnection system in the multi-energy market. Firstly, the basic structure of the electric-gas interconnection system is introduced, and the steady-state model of energy flow in the system is built based on the energy hub model. Secondly, in the framework of Nordic energy market, a bidding model for energy purchase and sale of the electric-gas interconnection system in the multi-energy market based on the idea of robust optimization is constructed, considering the uncertainty of energy price and the output power of renewable energy units in the system. Finally, empirical analysis based on the actual data is carried out, and the results prove the validity and superiority of the model. In this paper, aiming at the uncertainty of energy prices, a large number of scenes are generated by Latin hypercube sampling (LHS), and then k-means algorithm is used to reduce the scenes, so as to simulate typical scenes. Aiming at the uncertainty of the output power of the renewable energy unit in the system, cardinal uncertainty set is used to control the deviation between the actual output power and the predicted output power, so that the overall robustness of the model can be controlled. The proposed model can make decision-making independent of the accurate probability distribution of uncertainty factors, and is suitable for complex multi-energy systems. Meanwhile, the model possesses excellent robustness, which can effectively reduce the risk of bidding loss in the process of energy purchase and sale.

Typical Structure of the Electric-Gas Interconnection System
The electric-gas interconnection system consists of the natural gas system, the heating system, the cooling system and the power system distributed in the same area. The system components can be divided into external energy supply unit, internal energy supply unit, energy coupling unit, energy storage unit, energy supply network and load unit. The typical structure of the electric-gas interconnection system is shown in Figure 1, and the return water network for cooling and heating is not shown in this figure.  Power load in the power system is supplied by the distributed or cluster renewable energy system (RES), the combined heat and power (CHP) system, the energy storage (ES) system and the upper grid. Heat load in the heating system is supplied by the CHP system, the heat pump (HP) system and the heat storage system [22]. Cooling load in the cooling system is transformed and supplied by the absorption refrigeration (AR) system, the refrigerating machine (RM) system and the cooling storage system. Natural gas load in the gas supply system is supplied by the power-to-gas (P2G) system, external gas sources and gas storage tank. The natural gas pipeline network, as the low-pressure gas distribution network (7.5 × 103pa and below), is usually equipped with a pressure regulating device at the air sources, without considering the compressor and pressure regulating equipment temporarily. In conclusion, due to the introduction of a series of coupling elements, the system enhances the strong multi energy coupling relationship, which requires joint analysis.

Energy Hub Model of Energy Flow in the Electric-Gas interconnection system
The establishment of steady-state model of energy flow of multi-type heterogeneous energy in the electric-gas interconnection system is the basic premise for decision-makers to make decision optimization. The energy hub model is widely used in the steady-state modeling of energy flow [14,16]. In this model, for an energy conversion device with only single input and single output, the relationship between input and output is as follows: where, L  and P  are the input and output of the steady-state system;    ， is the coupling coefficient between the input and output. The integrated energy system consists of multiple energy conversion devices and various energy forms. The coupling relationship between input and output in the system can be described by the energy hub model. The formula is as follows: where, vectors P and L are the input and output variables of the system respectively. Each coupling coefficient in the coupling matrix  corresponds to the coupling relationship between a specific input and a specific output.
By combining the energy concentrator conversion matrix and the flow direction of different energy sources in the electric-gas interconnection system shown in Figure 1, the corresponding energy input and output relations of natural gas, heat, cooling and the power systems in the electricgas interconnection system can be obtained. Due to the bidding period =1 t h , the power and energy (electric quantity and natural gas quantity) of each bidding period are the same in numerical value. Without specification, the bidding power (Unit: MW) is adopted below to describe the model. The energy conversion model of the natural gas system is as follows: where, G,t L (Unit: MW) represents the output power of the natural gas system in the t period, and subscript t represents the time scale variable, the same below.    The energy conversion model of the heating system in the electric-gas interconnection system is as follows: represents the power allocated by the power system to the HP system.
The energy conversion model of the cooling system is as follows: where, C,t L (Unit: MW) represents the output power of the cooling system, The energy conversion model of the power supply system is as follows: where, G,t L represents the output power of the power system.
represents the predicted value of the actual power purchased or sold, which is positive when sold and negative when purchased.
represents the output power of the RE system.
ch arge E,t P represents the energy storage power of the ES system, and discharge E,t P represents the discharge power of the ES system.

Dynamic Energy Charging and Discharging Model of the Energy Storage (ES) System
Energy storage equipment is an important equipment of energy hub, which can transfer energy in a time dimension. Excess energy or cheap energy can be stored by energy storage equipment and released when the energy demand is large or the energy price is high, which can make efficient use of energy and reduce energy cost.
In a broad sense, energy storage includes not only electricity storage (battery), but also heat storage (the hot water storage system), gas storage (gas storage tank), cooling storage (ice storage, water storage, etc.). Their processes of energy charging and discharging are similar, which need to consider the constraints of multiple periods, mainly including the constraints of charging and discharging state, energy storage and so on. The general dynamic model of generalized energy storage can be expressed as follows: ch arg e disch arg e ch arg e ch arg e ch arg e ch arg e ch arg e min max t t t disch arg e disch arg e disch arg e disch arg e disch arg e min max t t t u P P u P   (9) disch arg e storage storage storage ch arg e ch arg e 1 disch arg e = 1-+ where, ch arg e t u and disch arg e t u represent the charging and discharging states of energy storage equipment respectively, which are 0-1 variables. 1 represents the working state, and 0 represents the non-working state. Equation (7)  Equation (12) shows that the energy storage of the energy storage equipment is equal at the end time and the initial time.

Market Structure and Bidding Process
Assuming that the electric-gas interconnection system is the price acceptor in both the power market and the natural gas market, the decision maker of the system only needs to optimize his bidding volume according to the predicted market transaction price when bidding. In the day-ahead power market and the day-ahead natural gas market, the bidding is conducted in hours, that is, there are 24 bidding periods in each trading day, and the decision maker of the system can choose to declare as the buyer or the seller in each period.
Considering the relative maturity of the Nordic power market, the research work of this paper is based on the simplified market mechanism of Nord pool. The bidding process of the Nordic power market is shown in Figure 2 The bidding process of the Nordic power market is as follows: the day-ahead market on the Dth day will be closed at 12:00 pm on the -1 D th ( ) day, and all participants need to submit the bidding information (combination of price and electric quantity) of T ( T = 24) trading periods on the Dth day to the dispatching authority before that.
The time interval from the close of the day-ahead market to the start of the real-time dispatching on the Dth day is long, and the predicted power supply and demand situation will generally deviate from the actual situation. The system dispatching operation organization can balance the supply and demand deviation through the intra-day market, and carry out unbalanced electric quantity settlement after the end of the dayʹs operation.
The generation of renewable energy in the electric-gas interconnection system is uncontrollable, and the bidding subject of the electric-gas interconnection system may use the price difference between the power market and the natural gas market for arbitrage, so the actual output power of the electric-gas interconnection system in the next day may deviate from the declared value of the previous output power. The independent dispatching organization needs to mobilize other adjustable units to balance the output power deviation, resulting in the re-dispatching cost of the system. In the Nordic power market, although the pricing mechanism of unbalanced electric quantity in different countries is slightly different, its formulation principle is to ensure that the income of unbalanced electric quantity is lower than the trading income in the day-ahead market or in the intraday market. Therefore, this paper charges the unbalanced electric quantity with the penalty price.
The Nordic natural gas trading market mainly includes within day contracts market, day contracts market and month contracts market. The natural gas buyer or seller submits the purchased or sold quantities and prices of natural gas to the natural gas trading hub, and the trading system automatically matches the order according to the bid of the buyer and the seller. At the end of the day, the trading system calculates the average trading price of the day by weighting according to the quantities and prices of natural gas on that day. This paper discusses the bidding strategy of the electric-gas interconnection system in the power market and the natural gas market. Since the transaction price of natural gas cannot be determined before the transaction actually occurs, and the fluctuation of the natural gas clearing price is small, the average transaction price of natural gas is taken as the price of natural gas purchase or sale approximately.

Treatment of Uncertain Factors
When the electric-gas interconnection system participates in the bidding of the multi-energy market, the uncertain variables involved include the clearing price of the power market, the average clearing price of the natural gas market, the penalty price of unbalanced electric quantity, and the actual output power of the system in the next day. Two kinds of uncertain variables, price and output power, are dealt with by the scenario analysis method and cardinal uncertainty set method respectively.
The clearing price of the power market, the clearing price of the natural gas market, and the penalty price of unbalanced electricity quantity will directly affect the income of the electric-gas interconnection system, while the clearing price is affected by the bidding behavior of the market subjects, load fluctuation and other factors, which has considerable volatility and uncertainty. When using the scenario analysis method to describe the uncertainty of these factors, it can deal with uncertainty factors and carry out risk analysis by calculating the maximum (or minimum) expectation of the model objective function in various scenarios.
First of all, according to the historical data of clearing price in the power market and the natural gas market, a specific probability distribution can be obtained by data fitting. Then, a large number of scenes are obtained by sampling. The larger the number of scenes are, the more accurate the results are, but the larger the calculation amount is. In order to balance the calculation amount and accuracy, firstly, according to the probability distribution of market clearing price, a large number of scenarios are generated by Latin hypercube sampling (LHS). Then, the k-means algorithm is used to reduce and merge the scene, and the probability of the scene after merging is calculated. The flow of the kmeans algorithm is shown in the  For the uncertainty of the next day's actual output power of the electric-gas interconnection system, select the appropriate 'uncertainty set' to be expressed. At present, the uncertain sets are often described in the form of box type, polyhedron type, ellipsoid type, cardinal type, etc. The characteristic of the cardinal uncertainty set used in this paper is that it can control the relative value of the uncertainty parameter offset and measure the parameter fluctuation accurately. The cardinal uncertainty set of the next day's actual output power of the electric-gas interconnection system is expressed as follows: where,  represents the uncertainty set of the next dayʹs actual output power of the electric-gas interconnection system.

Robust Model for Energy Purchase and Sale of the Electric-Gas Interconnection System
The model assumes that: (1) Both the power market and the natural gas market involved in the transaction are bilateral markets, and energy can be purchased or sold flexibly. (2) The scale of the electric-gas interconnection system is not enough to influence the power market and the natural gas market. (3) The positive deviation and negative deviation of the actual intra-day purchase or sale of electric quantity of the electric-gas interconnection system use the same penalty price to settle the unbalanced electric quantity. (4) In the natural gas market, the price is relatively stable and there is no deviation penalty mechanism.
Robust optimization is an optimization method that makes a decision based on the uncertainty of influencing factors. Its goal is to get the optimal decision under the worst case of uncertain factors. Therefore, when using robust optimization for decision-making, it is necessary to make clear what the worst impact of uncertain factors on profits is. Although the worst situation may not occur in reality, it is still necessary to consider and make decisions based on the worst case, so as to avoid the serious consequences brought by the worst-case as much as possible. Therefore, this paper establishes the purchase and sale decision-making model of the electric-gas interconnection system as follows where, E,t

Model Transformation
In the process of solving the above robust optimization model, considering that some components in Formula (14) are independent of variables, Formula (14) can be simplified as follows: The solution of robust optimization problem belongs to NP-hard problem mathematically, so this paper decomposes the original problem into main problem and sub problem. The main problem takes bid E,t P and C as decision variables, and its goal is to maximize the profit of the electric-gas interconnection system in all the worst cases of the sub problem. The model is as follows: Where, K in the main problem represents the index set of the iteration. E, ,

Parameter Setting
In order to simplify the model, the following assumptions are made: the daily clearing prices of the power market and the natural gas market in each period meet the normal distribution in a certain range; the mean value is taken from the actual market data of European energy exchange (EEX) and Nord pool, respectively; and the standard deviation is 3 EUR/MWh. Under the above data settings, 1000 daily clearing price scenarios of the natural gas market and the power market are generated respectively through LHS. After reduction, typical daily clearing price scenarios are obtained. The penalty price is set as 50% of the daily clearing price of the day-ahead market. The price data of typical daily clearing price scenarios are shown in Figure 4, the specific price data can be found in supplementary materials The structure of the electric-gas interconnection system in the example analysis of this paper is shown in Figure 1, the specific load data can be seen in supplementary materials. The load data of power, heat, cooling and natural gas in the calculation example is shown in Figure 5. The conversion efficiency and the upper and lower limits of each coupling unit output power are shown in Table 1.

Analysis of Optimized Operation Results
In order to verify the feasibility of the model and algorithm, this paper uses Matlab R2017a to compile the model mentioned above, and carries out trial calculation, the specific MATLAB code can be seen in supplementary materials When  = 6, the declared strategy and the optimal output power of each unit of the electric-gas interconnection system in the natural gas market and the power market in each period can be obtained, as shown in Figures 6-8.  Figure 6 shows the declared strategy of the electric-gas interconnection system in each period of the natural gas market and the power market, the specific decision result data can be found in supplementary materials, with positive value representing sale and negative value representing purchase. As can be seen from the figure, the electric-gas interconnection system sells natural gas in the natural gas market in the periods 1-5, purchases natural gas from the natural gas market in the periods 6-24, and sells power in all the periods. In general, due to the relatively low price of natural gas, the system decision-makers in this case tend to participate in the multi energy market transaction by purchasing natural gas and selling surplus electric energy, that is, internal demand is mostly transformed into supply through purchased natural gas. The uncontrollable component is the RE system, and the output power of electric energy in the system comes from the renewable energy units. After the output power of some renewable energy units is converted into other energy forms or supplies power to the coupling units, the remaining part is the declared value of the output power of electric energy in the system for bidding. Therefore, the optimal output power of the renewable energy units is actually the output power of the renewable energy unit ensuring the optimal purchase and sale strategy of the system, and it is not the actual output of the next day. If the decision-makers declare according to the declared strategy of energy purchase and sale corresponding to the optimal output power of the renewable energy units, the deviation of system output power will be minimized and the income will be maximized. In Figure 7, the upper and lower curves respectively correspond to the upper and lower limits of the RES output power in the system, and this value can be predicted according to the historical actual data. The shaded part is the predicted confidence interval of the RES output power within 24 hours. The red curve in the middle represents the optimized RES output power strategy. It can be seen that the output power decision of RES at all times is within the confidence interval of the RES output power within 24 hours, and the optimal value of the RES output power in the periods 8-9 is the upper limit value of allowable output power. According to  there are two factors that influence the optimal output to reach the upper limit: one is that the system power load is at the peak in the periods 8-9; the other is that the power price is at the peak in the periods 8-9. On the basis of meeting the internal energy demand, the system decision-makers sell power to the power market with the maximum capacity due to the profit seeking mentality. In the above scenario, the output power of energy coupling unit and energy storage unit in each period of the electric-gas interconnection system is as shown in Figure 8. It can be seen from the figure that P2G equipment is in full load state all day, HP and RM are in full load state except for some periods, CHP and AR fail to reach 50% of the rated power rate in most of the periods and stay in auxiliary state, and ES is charged in the periods 3-5 and discharged in the periods 9-10 and 20.

Sensitivity Analysis
RES is the internal energy supply unit of the electric-gas interconnection system. This section discusses the influence that the change of RES output power restriction has caused on the declared strategy of the energy purchase and sale and the output power of the internal coupling unit. Under each scenario, let the upper limit of RES output power change by 20%, 10%, ₋10% and ₋20% in each period, and compare the original optimization results with the four scenarios set in this section. The expected profits of the electric-gas interconnection system participating in the market bidding under the four scenarios are shown in Table 2. The expected profit is the expected value of the profit under the probability distribution of each price sample, while the actual profit is the average value of the obtained profit assuming that a certain sample in the price sample set will actually occur. It can be seen from the table that the expected profit generated by the system participating in the market bidding is positively related to the change of the upper limit of RES output power, and when the upper limit of RES output power changes in a negative direction, the change range of the expected profit is larger. When the upper limit of RES output power changes positively, only the optimal output power of RES in periods 8 and 9 changes (as shown in Figure 9), and when the upper limit of RES output power changes negatively, the optimal output power of RES in the periods 7-9, 11 and 19 changes simultaneously (as shown in Figure 9). Compared with the positive change, the negative change of the upper limit of RES output power is more significant, which is one of the reasons why the expected profit is more sensitive to the negative change of RES output power.  Figures 10-13 show the bidding strategy of the electric-gas interconnection system in four scenarios and the variance ratio of the declared strategy in each period compared with the control scenario. It can be seen from the figure that under the parameter setting in this paper, when the upper limit of the internal power supply unit output changes, the internal power supply unit output changes correspondingly, but the power declared strategy does not change. Considering that the power transaction price is more advantageous than the natural gas transaction price, the system decision-makers are more inclined to reduce the natural gas declaration volume in periods 8-9 on the premise of meeting the internal load demand, so as to ensure the optimal declared strategy of energy purchase and sale. In scenario 1 and scenario 2, the power consumption of internal energy supply units is promoted by reducing the purchase of natural gas in general, and the power generation consumption of the internal energy supply unit is reduced more in the 8th period when the natural gas purchase amount is relatively large, which directly leads to the increase of expected profits. In scenario 3 and scenario 4, the power generation reduction of the internal energy supply unit is generally compensated by increasing the purchase of natural gas, and the purchase of natural gas increases substantially in the periods 7-9 and 19, which also directly leads to the reduction of expected profits.     Figures 14-17 show the output variation ratio of the internal coupling unit in each period of the electric-gas interconnection system compared with the control scenario. It can be seen from the figure that under the parameter settings in this paper, when the upper limit of the output of the internal power supply unit in the electric-gas interconnection system changes, the output of the internal power supply unit changes correspondingly, the output of AR, P2G, RM does not change, and the output of HP and CHP changes correspondingly. In scenario 1 and scenario 2, the output of HP increases in periods 8 and 9, and the change proportion is positively related to the upper limit change of power supply unitʹs output. The reason for the change is that the RES output increases in periods 8 and 9, the power declared strategy remains unchanged, and the natural gas purchase volume is reduced. When the natural gas flow and electric energy flow change, HP and CHP will change their output correspondingly to maintain the stability of heat energy flow. In scenario 3 and scenario 4, due to the increase of gas purchase volume in some time periods, CHP and HP increase or decrease their output respectively in corresponding time periods to maintain the supply-demand balance of the system. In particular, the P2G output decreases in the eighth period, and the reason is that in the controlled situation CHP output has been close to the upper limit of output in the eighth period and the increased amplitude is compressed, so P2G needs to be called to maintain the balance of supply and demand in the system.

Robustness Test
In order to verify the robustness of the declared strategy calculated by the robust optimization model constructed in this paper, this section compares the expected profit of the system with the profit conducted by determining the actual output of the next day of the system, and analyzes the economic risks of the electric-gas interconnection system under different scenarios. Without considering the output uncertainty of the electric-gas interconnection system, the predicted value of the actual system output is simply set as the average value of the historical data. When  = 0, is the deterministic parameter. When considering the uncertainty of the system output, the robustness factor  > 0,  is set to 0, 6, 12 respectively, and the actual output of the next day is assumed to be As shown in Table 3, when  increases, the expected profit decreases, the profit deviation decreases gradually and the actual profit increases. When  = 0, i.e. without considering the uncertainty, if the actual output of the electric-gas interconnection system deviates greatly from the power generation plan, the system may lose more profits in the day-ahead bidding, and the profit deviation will disturb the normal operation of the electric-gas interconnection system. When  is 6 or 12, although the power deviation cannot be avoided and the profit deviation still exists, the profit deviation is obviously low and the actual profit is also significantly improved, which shows that considering the uncertainty in the optimization process can reduce the risk brought by the uncertainty factors to the profit. However, when considering the uncertainty of the system output, with the increase of the robustness factor  , the expected profit obtained by robust optimization will be reduced, so the decision-maker needs to make a reasonable decision according to his own risk tolerance.

Conclusions
Due to energy diversity and high coupling of the system, the bidding behavior of the electricgas interconnection system in the market is more flexible. Considering the complementary coupling characteristics of the power system and the natural gas system, this paper uses the idea of robust optimization to establish a robust optimization model for the purchase and sale of energy in the electric-gas interconnection system in the environment of the multi-energy market based on the energy hub model. In dealing with the uncertain parameters, the cardinal uncertainty set method is used to control the deviation between the actual output and the predicted output, so as to make the overall robustness of the model controllable. The results show that: (1) Under the condition that the price of power is superior to that of natural gas, the system decisionmakers tend to sell power and purchase natural gas to maintain the balance of supply and demand in the system and ensure the maximum bidding revenue. (2) The optimal bidding strategy for purchasing and selling energy, the output of the coupling unit and the expected profit of the electric-gas interconnection system are more sensitive when the upper limit of output of the internal power supply unit changes in the negative direction than when it changes in the positive direction, so the corresponding margin should be kept in the system planning. (3) When the robust parameter  increases, the expected profit decreases and the profit deviation decreases accordingly. The decision-maker should make a reasonable decision by adjusting the robust parameter according to his own risk tolerance. This paper focuses on the bidding strategy of the external energy purchase and sale and internal energy management of the electric-gas interconnection system, as well as the changing relationship between them, without considering the economic scheduling of the internal coupling unit. Therefore, in the next study, a refined bidding model of energy purchase and sale based on the economic scheduling of the coupling unit can be established accordingly.