Low-Carbon Economic Optimization of Integrated Energy System Considering Reﬁned Utilization of Hydrogen Energy and Generalized Energy Storage

: In order to improve the level of new energy consumption in the system and utilize the clean and efﬁcient characteristics of hydrogen energy, an integrated energy system (IES) scheduling model considering reﬁned utilization of hydrogen energy and generalized energy storage is proposed. Firstly, the two-stage hydrogen energy utilization model of power-to-gas (P2G) is ﬁnely modeled, and the waste heat of the P2G methanation reaction is innovatively coupled with the Kalina cycle to improve the thermoelectric decoupling capability of the combined heat and power (CHP) unit. Secondly, integrated demand response, electric vehicles, and hydrogen-containing multi-source energy storage equipment are used as generalized energy storage resources to cut peaks and ﬁll valleys. Then, on the basis of considering the ladder-type carbon trading mechanism, the IES conventional operation model is constructed with the minimum operating cost of the system as the objective function. Furthermore, considering the source-load uncertainty of IES operation, a multi-energy complementary optimal scheduling model of hydrogen-containing IES based on conditional value-at-risk was established. Through simulation analysis, it can be seen that the proposed model takes into account both economic and environmental beneﬁts and improves the system’s ability to “peak cutting and valley ﬁlling” and measure risk levels.


Introduction
With the further promotion of the "dual carbon" goal, a higher proportion of clean energy is connected to the power grid [1,2]. How to build an efficient, low-carbon, and clean energy system has become the key to solving environmental pollution and the energy crisis [3]. The IES integrates data processing [4], energy transmission, and other technologies that can achieve multi-energy collaborative optimization scheduling, facilitate energy cascade utilization, promote new energy accommodations, and reduce carbon emissions [5].
As a clean energy source with no carbon emissions and high combustion heat, hydrogen energy can be combined with different energy sources in IES to achieve low-carbon and efficient scheduling. Through P2G technology, the surplus electric energy can be converted into methane [6], and the reference [7] refines the P2G operation model on this basis by adding electric to hydrogen and hydrogen to electricity links to realize the interaction between electrical energy and hydrogen energy. In [8], the waste heat from hydrogen to methane is recovered and used as heating power, while the hydrogen-containing natural gas pipeline network and hydrogen fuel cell vehicles are added to further improve the utilization of hydrogen energy. In [9,10], IES recovers methanation waste heat by using carbon capture equipment to improve the energy efficiency of the system and introduces a green certificate-carbon trading mechanism to increase hydrogen production and wind

IES Structure Considering GES and Hydrogen Energy Utilization
In this paper, IES includes three links: the energy supply side, the energy convers side, and the load demand side. At the same time, the regulation role of actual ene storage and virtual energy storage is considered to realize the multi-energy coupling timization operation of IES. The supply side includes wind turbines (WT), photovol (PV), gas grids, and the superior power grid. On the energy conversion side, the elect ity-hydrogen coupling link includes methane reactors (MR), electrolyzers (EL), hydrog fuel cells (HFC), and P2G waste heat recovery devices. In addition, the GT and Kal cycles can be used to generate electricity. Waste heat boilers (WHB), gas boilers (GB), a power-to-heat (P2H) are used for heat production. Electricity storage (ES), gas stor (GS), thermal storage (TS), and hydrogen storage (HS) devices can be used as actual ergy storage devices. On the load demand side, EV charging and discharging models a IDR can be used as virtual energy storage resources to flexibly schedule multiple ene sources in IES. The system structure diagram is shown in Figure 1.  EL produces hydrogen by electrolysis of water, which is the key equipment in electricity-hydrogen coupling link. EL is beneficial for IES to utilize the pollution-free a good combustion performance characteristics of hydrogen energy. Meanwhile, the c rent research shows that the energy conversion efficiency of EL electrolytic hydrogen p duction can reach 80%, which is much higher than the efficiency of electrolytic natural production. Therefore, this paper improves the cleanliness and efficiency of IES ene use through the electricity-hydrogen coupling link.

Electricity-Hydrogen Coupling Link
The electricity-hydrogen coupling link models the two-stage hydrogen energy u zation model of P2G, including electrolytic hydrogen production and hydrogen ene utilization. The hydrogen energy utilization link includes the hydrogen to methane c version process of MR and the thermal and electrical conversion process of HFC, wh can realize the cascade utilization of hydrogen energy. In addition, in order to impr the efficiency of energy utilization, this paper introduces a P2G waste heat recovery dev to recycle the waste heat of the hydrogen methanation reaction.
1. EL equipment model EL uses water as the raw material and converts electric energy to hydrogen ene through the electrolytic hydrogen production process. The resulting hydrogen ene EL produces hydrogen by electrolysis of water, which is the key equipment in the electricity-hydrogen coupling link. EL is beneficial for IES to utilize the pollution-free and good combustion performance characteristics of hydrogen energy. Meanwhile, the current research shows that the energy conversion efficiency of EL electrolytic hydrogen production can reach 80%, which is much higher than the efficiency of electrolytic natural gas production. Therefore, this paper improves the cleanliness and efficiency of IES energy use through the electricity-hydrogen coupling link.

Energy Conversion Equipment Model 2.2.1. Electricity-Hydrogen Coupling Link
The electricity-hydrogen coupling link models the two-stage hydrogen energy utilization model of P2G, including electrolytic hydrogen production and hydrogen energy utilization. The hydrogen energy utilization link includes the hydrogen to methane conversion process of MR and the thermal and electrical conversion process of HFC, which can realize the cascade utilization of hydrogen energy. In addition, in order to improve the efficiency of energy utilization, this paper introduces a P2G waste heat recovery device to recycle the waste heat of the hydrogen methanation reaction. EL equipment model EL uses water as the raw material and converts electric energy to hydrogen energy through the electrolytic hydrogen production process. The resulting hydrogen energy flows into MR, HFC, and HS, respectively. The equipment model for EL is as follows: where P EL,t and P H,EL,t denote the input electric power and output hydrogen energy of EL at time t; η EL is the hydrogen production efficiency of EL; P max EL and P min EL are the maximum and minimum power consumption of the EL; ∆P max EL and ∆P min EL are the upper and lower limits of the ramp rate of EL.

2.
MR equipment model MR converts the hydrogen energy input from EL into CH 4 through the methanation reaction and supplies it to the gas load and gas unit. The MR equipment model is as follows: where P H,MR,t and P g,MR,t denote the input hydrogen power and output gas power of MR at time t; η MR is the CH 4 where P HFC,t , P e,HFC,t , and P h,HFC,t are, respectively, the hydrogen energy input and the electricity and heat energy output of HFC at time t; η HFC is the energy conversion efficiency of HFC; P max HFC and P min HFC are the maximum and minimum hydrogen consumption power of HFC; ∆P max HFC and ∆P min HFC are the upper and lower limits of the ramp rate of HFC. k max HFC and k min HFC are the upper and lower limits of the adjustable heat-to-electric ratio of HFC. 4.
P2G waste heat recovery device model The chemical reaction process of P2G in two-stage operation is shown in Equations (4) and (5).
where ∆H is the entropy of a chemical reaction. When its value is negative, it means that the chemical reaction is exothermic. Formula (5) shows that P2G releases heat during the methanation reaction. In order to avoid heat waste and improve energy efficiency, the waste heat from the P2G reaction is Energies 2023, 16, 5700 5 of 23 recovered by the waste heat recovery device. The heat released by the P2G methanation reaction in unit time H CH 4 is calculated by the following formula: where ∆H and P EL are methanation reaction heat release and EL power consumption; v H 2 is the rate of electrolytic hydrogen production; ρ H 2 and M H 2 are the density and molar mass of hydrogen; η CH 4 is the proportion of reaction heat recovered by the heating network. According to Equation (4), P2G can recover 0.1188 kW·h of reaction heat for every 1 kW·h of electric energy consumed. Based on this, the thermal power of P2G waste heat recovery is obtained as follows: where η P2G and P h,P2G,t are P2G waste heat recovery efficiency and heat power recovered at time t.

Thermoelectric Decoupling Model of CHP Unit
Traditional CHP units are prone to high forced power output when meeting thermal loads, resulting in issues such as determining power by heat and energy waste. In this paper, the thermoelectric decoupling model of CHP units is introduced by coupling the P2G waste heat recovery device with the Kalina cycle.

1.
Gas turbine where P GT,t , P max GT , and P min GT are the input gas power of GT and its upper and lower limits; P e,GT,t and P h,GT,t are the electric power and thermal power output by GT; P WHB,t and P h,KLN,t denote the thermal power flowing into WHB and Karina cycles, respectively; η gte and η gth are the gas-to-electricity and gas-to-thermal efficiency of GT, respectively. ∆P max GT and ∆P min GT are the upper and lower limits of the ramp rate of GT.

2.
Waste heat boiler model where P WHB,t , P max WHB , and P min WHB denote the input power of WHB and its upper and lower limits; η whb and P h,WHB,t denote the waste heat recovery efficiency and output thermal power of WHB; P h,CHP,t denotes the thermal power output of the CHP.

3.
Coupling model of Kalina cycle and P2G waste heat recovery device In this paper, the waste heat of GT power generation and P2G methanation waste heat are recovered, and waste heat power generation is carried out through the Kalina cycle to realize secondary utilization of energy, further improve the energy supply flexibility of CHP units, and achieve the optimal ratio of heating and power supply.    P e,KLN,t = η kln (P h,KLN,t + P h,P2G,t ) P min KLN ≤ P e,KLN,t ≤ P max KLN P e,CHP,t = P e,GT,t + P e,KLN,t (10)  where P e,KLN,t , P max KLN , and P min KLN denote the output power of Kalina cycle and its upper and lower limits; η kln denote the thermoelectric conversion efficiency of the Kalina cycle; and P e,CHP,t is the output power of CHP.

P2H and GB Model
P2H converts electric energy into heat energy during periods of low electricity prices, which can effectively utilize the curtailment of wind and photovoltaic power and adjust the power load curve. GB outputs thermal power by burning natural gas and compensates for thermal energy when waste heat boiler heat is insufficient. The power model of P2H and GB is shown in Equation (11).
where P i,t and P h,i,t denote the input power and output thermal power of device i at time t, respectively, i ∈ {P2H, GB}; η i is the energy conversion efficiency of device i; P max The multi-source energy storage model in this paper includes ES, TS, GS, and HS with similar charging and discharging models, which can be expressed by the following equations.
where i is the type of actual energy storage device, i ∈ {ES, TS, GS, HS}; S i,t and µ i are the capacity and energy self-loss coefficient of the energy storage device at time t; P i,ch,t and P i,dis,t are the charging and discharging power of the energy storage device at time t; S i,0 and S i,24 are the energy of the energy storage device at the beginning and end of a day, respectively.

IDR Model
Different flexible loads have different responses to electricity price changes and incentive subsidies. In this paper, the electric load and gas load are affected by price factors. Therefore, the IDR of electric load and gas load includes price-based, incentive-based, and alternative DR. Heat load DR only considers transferable, flexible heat loads.

1.
Price-based DR The price elasticity matrix is used to realize the time shift of flexible electric load and gas load, which promotes users to realize "peak cutting and valley filling". The price elasticity coefficient m i,t,k is calculated as follows: where i is the type of energy, e denotes electric energy, and g denotes gas energy; m i,t,k denotes the energy price elasticity coefficient between time t and time k of load i; P i,0 and ∆P i,t are the initial load and the load response after DR; q i,k,0 and ∆q i,k are the initial energy price at time k and the change of energy price after DR. Furthermore, the mathematical model of price-type electrical load and gas load DR is obtained.
where P i,t,0 and ∆P i,TS,t denote the energy i load before the price-based DR and the load change after the DR; ∆P max i,TS,t is the upper limit of the change in the DR of energy i at time t; M i denotes the price elasticity matrix of energy i, which is composed of m i,t,k ; T is the scheduling period, which is set to 24 h.

2.
Incentive DR (1) Electric load and gas load Incentive DR Users can sign a contract with the power system to reduce a certain amount of electric and gas load within a specified time range and obtain compensation, the formula is as follows: where ∆P i,cut,t and P i,cut,t are the total reduction and load of energy i after DR at time t; ϕ c is the upper limit of the total change rate of the load.
(2) Heat load DR Transferable flexible heat loads are considered in the heat load DR, such as the airconditioning heat load, which can adjust its operation period according to the peak and valley characteristics of the load. The mathematical model of the flexible heat load is as follows: P h,load,t = P h,load,t0 − P h,out,t + P h,in,t P h,out,t P h,in,t = 0 (20) where P h,load,t and P h,load,t0 denote the heat load before and after the DR at time t; P h,out,t and P h,in,t denote the heat load transferred out and transferred in at time t, respectively. IES provides incentive compensation to users according to the peak shaving effect of heat load DR, and the formula is as follows: where C DR,h,t is the total compensation cost of heat load DR; ζ 1 is the compensation cost per unit of peak regulation power.

Alternative demand response
The alternative demand response does not affect the user's energy use experience, and the user can choose the appropriate energy use mode according to the electricity price and gas price data published by the operator. The formula is as follows: ∆P re,g,t = −ε e,g ∆P re,h,t (22) ε e,g = v e ϕ e v g ϕ g ∆P min re,i,t ≤ ∆P re,i,t ≤ ∆P max re,i,t , i ∈ {e, g} where ∆P re,g,t and ∆P re,h,t denote the corresponding increments of gas load and electric load at time t; v e and v g are the unit calorific value of electric energy and gas energy, respectively; ϕ e and ϕ g are the energy efficiency of electric energy and gas energy; ε e,g is the conversion coefficient of electric energy and gas energy; ∆P min re,i,t and ∆P max re,i,t are the lower limit and upper limit of the alternative quantity of load i. Therefore, the electric load and gas load can be expressed by Equation (25).

Charging and Discharging Models of EV
EV grid-connected models include disordered charging and ordered charge-discharge models. Among them, EV disordered charging models are unscheduled resources, and the probability density function of the EV return time and the all-day driving distances obey the normal distribution. Therefore, the load calculation method for EV disordered charging can be referred to [24].
EVs based on V2G technology can be used as virtual energy storage devices and are regarded as flexible power loads to regulate load fluctuations [25]. The IES dispatching center predicts the travel demand of EV users based on historical statistical data and formulates reasonable charging and discharging strategies that are conducive to peak shaving and valley filling of the power grid. This paper mainly studies the electric energy interaction models with IES charging stations when EVs are connected to the grid based on V2G technology, as shown in Equation (26).
where δ ch,k,t and δ dis,k,t are 0-1 variables, when δ ch,k,t and δ dis,k,t take 1, it means EV is in charging and discharging states, respectively; otherwise, it takes 0; σ k,t is the state of EV. In addition, σ k,t = 1 and σ k,t = 0 denote the grid-connected and off-grid state of EV k at time t; P EV ch,k,t and P EV dis,k,t denote the charging and discharging power of EV k at time t; P EV,max ch,k and P EV,max dis,k are the maximum charging and discharging power of the EV k; S EV k,t is the battery capacity of the EV k at time t; η EV ch and η EV dis are the charging and discharging efficiency of EV. S EV,min k and S EV,max k are the minimum and maximum battery capacities of EV k. The subsidy incentive C EV mot for EVs discharge is shown in Equation (27).
where p e is the subsidized unit price of IES for EVs discharge power.

Scenario Analysis Considering Source-Load Uncertainty
The scenario analysis method includes scenario generation and scenario reduction. In the scenario generation, the actual values of renewable energy output and load power are superimposed by their predicted values and error values, where the prediction error obeys the normal distribution, that is, ε ∼ N 0, σ 2 . A large number of scenarios satisfying the prediction error probability distribution can be generated by the Latin Hypercube sampling method with high sampling accuracy. Then, the scenario reduction is carried out by the Energies 2023, 16, 5700 9 of 23 backward reduction method based on Kantorovich distance [26], and finally, k typical scenarios of wind-photovoltaic load are obtained.

Stepped-Type Carbon Trading Mechanism
Compared with the traditional carbon trading mechanism, the stepped-type carbon trading mechanism adopted in this paper sets multiple carbon emission intervals with increasing prices. The larger the carbon emission rights purchased, the higher the purchase price of the corresponding interval. The method of carbon emission quota allocation in IES is shown as follows: where P buy,t is the amount of electricity purchased at time t; β e and β h denote the carbon emission quota of the unit power and heat. β e,h denotes the electric-thermal conversion coefficient.
The cost of ladder-type carbon trading is shown in Formula (29) [27].
where F CO 2 is the ladder-type carbon transaction cost; E a is the actual carbon emissions; λ and ρ are the base price and price growth rate of carbon trading; d is the length of the carbon emission interval.

IES Basic Scheduling Model
This paper comprehensively considers the low-carbon and economy of IES operation, so the objective function is composed of IES operation cost and carbon transaction cost. The operation cost of IES includes energy purchase cost F BUY , IDR cost F IDR , and equipment operation and maintenance cost F OM . The objective function is shown as follows: where F BUY , F IDR , and F OM are shown in Formulas (31)-(33); F CO 2 is shown in Equation (29).
where χ e,t and χ h,t are the unit prices of electricity and gas at time t; P h,buy,t is the amount of natural gas purchased at time t.
where c i,TS , c i,re , and c i,cut are the unit power compensation prices of users participating in IDR.
where U denotes the set of all devices in IES; x is the device type; P x,t and ψ x are the unit operation and maintenance costs of device x and the output of device x at time t.

IES Economic Scheduling Model Considering CVaR
CVaR can be used to measure the risks in the course of IES operations and assess the loss range of IES under specific risks. The scenario probability obtained by the scenario generation method is ϕ k . Therefore, the expected operating cost of IES under scenario k is as follows: where F BUY,k , F IDR,k , F OM,k , and F CO 2 ,k are the energy purchase cost, IDR cost, equipment operation and maintenance cost, and carbon transaction cost in scenario k, respectively. The CVaR mathematical model of IES economic scheduling is as follows: where α is the VaR value; β is the confidence level; K is the number of scenarios; Through the risk preference coefficient δ, the expected costs of IES are combined with the average loss costs, δ ∈ (0, 1). The overall objective function is as follows: where F EX is the expected operating costs of IES.

1.
Electrical power balance constraints 0 ≤ P buy,t ≤ P max buy (37) P buy,t + P WT,t + P PV,t = P e,Load,t + P e,ch,t − P e,dis,t + P P2H,t + P EL,t − P e,HFC,t − P e,GT,t − P e,KLN,t + n * (P EV ch,k,t − P EV dis,k,t ) where P max buy is the upper limit of purchased power; P WT,t and P PV,t are the wind power and photovoltaic power at time t; n is the number of scheduled EVs in IES.

2.
Gas power balance constraints 0 ≤ P g,buy,t ≤ P max g,buy (39) where P max g,buy is the upper limit of purchased gas power.

4.
Hydrogen power balance constraint

5.
WT and PV constraints where P max WT and P max PV are the upper limit values of WT and PV output.
where i ∈ {ES, TS, GS, HS}; U i,ch,t and U i,dis,t are 0-1 variables, when U i,ch,t and U i,dis,t take 1, it means the EV is in charging and discharging states, respectively; otherwise, it takes 0.

Solution Process
The two-stage optimization model constructed in this paper considers the refined utilization of hydrogen energy in the first stage. The second stage considers the role of IDR and EV virtual energy storage on the basis of the first stage. The model solution uses the YALMIP toolbox of the Matlab platform for compilation, and the solver part uses Gurobi. The solution flow chart is shown in Figure 2.

Solution Process
The two-stage optimization model constructed in this paper considers the refined utilization of hydrogen energy in the first stage. The second stage considers the role of IDR and EV virtual energy storage on the basis of the first stage. The model solution uses the YALMIP toolbox of the Matlab platform for compilation, and the solver part uses Gurobi. The solution flow chart is shown in Figure 2.  Figure 2. Flow chart of model solution.
As shown in Figure 2, in the first stage, the P2G two-stage operation is first modeled. Secondly, the Kalina cycle is added to realize the thermoelectric decoupling of the CHP unit. Finally, the P2G waste heat recovery device is coupled with the Kalina cycle to realize the cascading utilization of hydrogen energy. In the second stage, the charging and discharging models of EVs and IDRs are regarded as virtual energy storage resources to realize peak cutting and valley filling. On this basis, typical scenarios are set up to analyze the conventional operation model. Finally, considering the influence of source-load uncertainty, the IES operation model considering CVaR is analyzed.

Parameter Settings
In this paper, a park-level integrated energy system in northern China during the winter heating period is selected as the research object. The IES topology adopted is shown in Figure 1. Among them, equipment parameters and actual energy storage parameters are shown in Tables A1 and A2, respectively. The PV, WT, and multi-source load As shown in Figure 2, in the first stage, the P2G two-stage operation is first modeled. Secondly, the Kalina cycle is added to realize the thermoelectric decoupling of the CHP unit. Finally, the P2G waste heat recovery device is coupled with the Kalina cycle to realize the cascading utilization of hydrogen energy. In the second stage, the charging and discharging models of EVs and IDRs are regarded as virtual energy storage resources to realize peak cutting and valley filling. On this basis, typical scenarios are set up to analyze the conventional operation model. Finally, considering the influence of source-load uncertainty, the IES operation model considering CVaR is analyzed.

Parameter Settings
In this paper, a park-level integrated energy system in northern China during the winter heating period is selected as the research object. The IES topology adopted is shown in Figure 1. Among them, equipment parameters and actual energy storage parameters are shown in Tables A1 and A2, respectively. The PV, WT, and multi-source load forecasting curves are shown in Figure 3, and the time-of-use electricity price and gas price curves are shown in Figure 4. The park contains 60 scheduled EVs, and the operating parameters of EVs are shown in Table A3. In the ladder-type carbon trading mechanism, the base price of carbon trading is 0.250 yuan/kg, the interval length is 2000 kg, the price growth rate of each ladder is 25%, and the actual carbon emissions calculation parameters can be referred to [28]. forecasting curves are shown in Figure 3, and the time-of-use electricity price and gas price curves are shown in Figure 4. The park contains 60 scheduled EVs, and the operating parameters of EVs are shown in Table A3. In the ladder-type carbon trading mechanism, the base price of carbon trading is 0.250 yuan/kg, the interval length is 2000kg, the price growth rate of each ladder is 25%, and the actual carbon emissions calculation parameters can be referred to [28].  As shown in Figures 3 and 4, the wind power output is large at night, and the electric load demand is small, which makes it easy to cause wind curtailment. Meanwhile, the heat load output is large, and the energy supply pressure is large. During the daytime, the demand for electricity and gas is large, and the price of electricity and gas is also high. Therefore, the cost of purchasing electricity and gas is high. Meanwhile, the disordered EV charging load will further increase the power supply pressure.

Analysis of the Combined Operation of Kalina Cycle and Electricity-Hydrogen Coupling Link
To verify the validity of the model proposed in this paper, the influence of uncertain factors on IES is first ignored, and the conventional basic operation scenarios are analyzed by using renewable energy and load forecasting power. On the basis of considering the ladder-type carbon trading mechanism and the disordered charging load of EVs, the following four scenarios are set to analyze the joint operation of the Kalina cycle and the electricity-hydrogen coupling link.   forecasting curves are shown in Figure 3, and the time-of-use electricity price and gas price curves are shown in Figure 4. The park contains 60 scheduled EVs, and the operating parameters of EVs are shown in Table A3. In the ladder-type carbon trading mechanism, the base price of carbon trading is 0.250 yuan/kg, the interval length is 2000kg, the price growth rate of each ladder is 25%, and the actual carbon emissions calculation parameters can be referred to [28].  As shown in Figures 3 and 4, the wind power output is large at night, and the electric load demand is small, which makes it easy to cause wind curtailment. Meanwhile, the heat load output is large, and the energy supply pressure is large. During the daytime, the demand for electricity and gas is large, and the price of electricity and gas is also high. Therefore, the cost of purchasing electricity and gas is high. Meanwhile, the disordered EV charging load will further increase the power supply pressure.

Analysis of the Combined Operation of Kalina Cycle and Electricity-Hydrogen Coupling Link
To verify the validity of the model proposed in this paper, the influence of uncertain factors on IES is first ignored, and the conventional basic operation scenarios are analyzed by using renewable energy and load forecasting power. On the basis of considering the ladder-type carbon trading mechanism and the disordered charging load of EVs, the following four scenarios are set to analyze the joint operation of the Kalina cycle and the electricity-hydrogen coupling link.  As shown in Figures 3 and 4, the wind power output is large at night, and the electric load demand is small, which makes it easy to cause wind curtailment. Meanwhile, the heat load output is large, and the energy supply pressure is large. During the daytime, the demand for electricity and gas is large, and the price of electricity and gas is also high. Therefore, the cost of purchasing electricity and gas is high. Meanwhile, the disordered EV charging load will further increase the power supply pressure.

Analysis of the Combined Operation of Kalina Cycle and Electricity-Hydrogen Coupling Link
To verify the validity of the model proposed in this paper, the influence of uncertain factors on IES is first ignored, and the conventional basic operation scenarios are analyzed by using renewable energy and load forecasting power. On the basis of considering the ladder-type carbon trading mechanism and the disordered charging load of EVs, the following four scenarios are set to analyze the joint operation of the Kalina cycle and the electricity-hydrogen coupling link.
Scenario 1: Only the electricity-hydrogen coupling link is considered. Scenario 2: The Kalina cycle is added on the basis of Scenario 1. Scenario 3: A P2G waste heat recovery  Table 1. Compared with Scenario 1, Scenario 2 considers the Kalina cycle and dynamically adjusts the electrothermal ratio of CHP units through the Kalina cycle to achieve thermoelectric decoupling of CHP units. The GT output thermal power distribution diagram in Scenario 2 is shown in Figure 5. Scenario 1: Only the electricity-hydrogen coupling link is considered. Scenario 2: The Kalina cycle is added on the basis of Scenario 1. Scenario 3: A P2G waste heat recovery device is added and coupled with the Kalina cycle. The running results of scenarios 1 to 3 are shown in Table 1. Compared with Scenario 1, Scenario 2 considers the Kalina cycle and dynamically adjusts the electrothermal ratio of CHP units through the Kalina cycle to achieve thermoelectric decoupling of CHP units. The GT output thermal power distribution diagram in Scenario 2 is shown in Figure 5. As can be seen from Figure 5, during the low electricity price period, the heat load demand is large, and all the heat generated by GT is recovered by WHB. Therefore, the electric power of the CHP unit is output by GT. In order to reduce the power purchase cost of the system during the peak period of electricity prices, the proportion of GT output thermal power flowing into the Kalina cycle continues to increase, so that it can reach maximum power generation. When wind curtailment occurs, the Kalina cycle can extract part of the steam for heat production, increasing the new energy consumption space by reducing the forced electrical output of the units. According to Table 1, compared with Scenario 1, the total operating cost of Scenario 2 is reduced by 812.3 yuan, the carbon emissions are reduced by 508.27kg, and the consumption rate of new energy is increased by 5.18%.

Comparative analysis of Scenario 2 and Scenario 3
Based on Scenario 2, Scenario 3 considers the coupling effect of the Kalina cycle and P2G waste heat recovery device, and the electric power of the CHP unit in Scenario 3 is shown in Figure 6. As can be seen from Figure 5, during the low electricity price period, the heat load demand is large, and all the heat generated by GT is recovered by WHB. Therefore, the electric power of the CHP unit is output by GT. In order to reduce the power purchase cost of the system during the peak period of electricity prices, the proportion of GT output thermal power flowing into the Kalina cycle continues to increase, so that it can reach maximum power generation. When wind curtailment occurs, the Kalina cycle can extract part of the steam for heat production, increasing the new energy consumption space by reducing the forced electrical output of the units. According to Table 1, compared with Scenario 1, the total operating cost of Scenario 2 is reduced by 812.3 yuan, the carbon emissions are reduced by 508.27 kg, and the consumption rate of new energy is increased by 5.18%.

2.
Comparative analysis of Scenario 2 and Scenario 3 Based on Scenario 2, Scenario 3 considers the coupling effect of the Kalina cycle and P2G waste heat recovery device, and the electric power of the CHP unit in Scenario 3 is shown in Figure 6. As can be seen from Figure 6, the law of GT thermal power flowing into the Kalina cycle is the same as that in Scenario 2. However, in Scenario 3, the Kalina cycle can also use the waste heat output power of the methanation reaction in P2G during the night low period from 23:00 to 5:00 to increase the output power of the CHP unit. In Scenario 3, the Kalina cycle increases the output power of the CHP unit while ensuring that the output power of WHB and GT is not reduced. Compared to Scenario 2, Scenario 3 reduced the electricity purchase cost by 160.69 yuan. Therefore, in Scenario 3, the scheduling enthusiasm of P2G is improved, and the power of EL electrolytic hydrogen production is improved, which is conducive to the consumption of wind power at night. It can be seen from Table 1 that, compared with Scenario 2, the total cost of Scenario 3 is reduced by 184.01 yuan and the carbon emissions are reduced by 180.75 kg. In addition, Scenario 3 effectively reduces operating costs by 996.31 yuan compared to Scenario 1. Compared with the method of directly utilizing P2G waste heat for heating in the literature [8], the strategy proposed in this paper effectively reduces the operating cost by 115.01 yuan and the carbon emission by 134.41kg, which proves that the method proposed in this paper has obvious advantages.

Ladder-Type Carbon Trading and GES Operation Analysis
In order to analyze the impact of ladder-type carbon trading mechanisms and GES on the operation of IES, Scenarios 4 to 6 are set. Scenario 4: Based on Scenario 3, only the traditional carbon trading mechanism is considered. Scenario 5: IDR is added on the basis of Scenario 3. Scenario 6: Based on Scenario 5, the EV charging and discharging model is considered. The running results of scenarios 4 to 6 are shown in Table 2. Compared with Scenario 4, Scenario 3 considers the ladder-type carbon trading mechanism, and the cost of electricity purchase has increased, while the cost of gas As can be seen from Figure 6, the law of GT thermal power flowing into the Kalina cycle is the same as that in Scenario 2. However, in Scenario 3, the Kalina cycle can also use the waste heat output power of the methanation reaction in P2G during the night low period from 23:00 to 5:00 to increase the output power of the CHP unit. In Scenario 3, the Kalina cycle increases the output power of the CHP unit while ensuring that the output power of WHB and GT is not reduced. Compared to Scenario 2, Scenario 3 reduced the electricity purchase cost by 160.69 yuan. Therefore, in Scenario 3, the scheduling enthusiasm of P2G is improved, and the power of EL electrolytic hydrogen production is improved, which is conducive to the consumption of wind power at night. It can be seen from Table 1 that, compared with Scenario 2, the total cost of Scenario 3 is reduced by 184.01 yuan and the carbon emissions are reduced by 180.75 kg. In addition, Scenario 3 effectively reduces operating costs by 996.31 yuan compared to Scenario 1. Compared with the method of directly utilizing P2G waste heat for heating in the literature [8], the strategy proposed in this paper effectively reduces the operating cost by 115.01 yuan and the carbon emission by 134.41 kg, which proves that the method proposed in this paper has obvious advantages.

Ladder-Type Carbon Trading and GES Operation Analysis
In order to analyze the impact of ladder-type carbon trading mechanisms and GES on the operation of IES, Scenarios 4 to 6 are set. Scenario 4: Based on Scenario 3, only the traditional carbon trading mechanism is considered. Scenario 5: IDR is added on the basis of Scenario 3. Scenario 6: Based on Scenario 5, the EV charging and discharging model is considered. The running results of scenarios 4 to 6 are shown in Table 2. Compared with Scenario 4, Scenario 3 considers the ladder-type carbon trading mechanism, and the cost of electricity purchase has increased, while the cost of gas purchase has decreased significantly. The reason is that in Scenario 4, the gas unit is in a state of high energy consumption, so the amount of gas purchased is large, while in Scenario 3, the carbon price of carbon emissions beyond the rated range presents a step growth. Therefore, IES increases the proportion of power purchased to reduce the output of gas units with large carbon emissions. As can be seen from Table 2, the equivalent carbon price of Scenario 4 is lower, so the carbon transaction cost of Scenario 3 is only reduced by 91.64 yuan, but the carbon emission of Scenario 3 is reduced by 8499.98 kg, and the total cost is reduced by 2907.92 yuan.

2.
Comparative analysis of Scenario 3, Scenario 5 and Scenario 6 Based on Scenario 3, Scenario 5 proposes IDR strategies to improve the sensitivity of users to electricity and gas prices released by operators, encourage users to adjust energy use mode and energy use period, relieve energy supply pressure during peak hours, and realize peak cutting and valley filling. As can be seen from Tables 1 and 2, compared with Scenario 3, Scenario 5 significantly reduces the cost of electricity and gas purchases by 3261.56 yuan and carbon emissions by 286.11 kg. Although Scenario 5 increases IDR cost by 960.17 yuan, it is far less than the economic and environmental benefits brought by IDR. Scenario 6 considers EV charging and discharging models on the basis of Scenario 5. As a virtual energy storage device, EV can not only reduce the electricity load during peak hours but also implement discharge operations, further reduce the cost of electricity purchase, and play the roles of peak cutting and valley filling. In Table 2, the IDR cost includes the EV subsidy incentive of 382.97 yuan. Compared with Scenario 3, Scenario 6 effectively reduces operating costs by 3407.47 yuan and carbon emissions by 620.13 kg. It can be seen from Table 2 that even if the IDR cost in Scenario 6 is higher, Scenario 6 still has the best operating economy and the lowest carbon emissions, which proves the effectiveness of the strategy adopted in this paper. Figure 7 shows the power balance diagram of each energy source in Scenario 6. Combined with the optimization results of multi-source energy storage in Figure 8, the following analysis can be performed.

Analysis of IES Scheduling Results
As shown in Figures 7 and 8, from 23:00 to 5:00, the valley period of electricity prices, IES mainly meets the electricity demand through power purchase and wind power generation. Both ES and EVs are charged, while EL produces hydrogen to reduce costs. At the same time, the heat load demand at night is large, the CHP units and the P2H output are large, and the TS stores heat during this period. The gas load demand is met by gas purchases, GS, and MR, and the hydrogen energy required by MR is provided by EL and HS. ES, EV, and HFC output electric energy, and GT output gradually increases during the peak periods of electricity prices from 06:00 to 12:00 and 19:00 to 22:00. As a result, the amount of electricity purchased is gradually decreasing. The heat load demand is reduced, and the heat is mainly provided by CHP and GB. Meanwhile, TS and HFC produce thermal energy. The gas purchase of IES increases, but the gas price is high during this time, so GS mainly releases energy. In 13:00 to 18:00, which is the normal period for electricity prices, the proportion of photovoltaic power generation in IES gradually increases, and ES stores excess electric energy. At this time, the heat load demand is reduced, and the gas purchase amount is also reduced. GS stores gas from 13:00 to 14:00, which is the normal period for gas prices. As can be seen from Figure 8, the state of charge (SOC) of ES and HS tends to be the same because, in periods of low electricity prices, the wind power output is large, so both ES and HS are stored energy. During the peak period of energy supply, ES discharges. Meanwhile, the HFC outputs electric and thermal energy, and the hydrogen consumption is provided by HS. In summary, the IES scheduling model proposed in this paper reduces energy purchase costs and makes multi-source energy storage equipment play a role in peak regulation. the strategy adopted in this paper. Figure 7 shows the power balance diagram of each energy source in Scenario 6. Combined with the optimization results of multi-source energy storage in Figure 8, the following analysis can be performed.   As shown in Figure A1, the optimization curves of electricity, heat, and gas loads before and after IDR are shown. Through IDR, the electricity and gas load in the peak period of electricity and gas prices are transferred to the low period and normal period to reduce the operating costs of IES. At the same time, the heat load with higher power at night is partially transferred to daytime by IDR. After IDR optimization, the peak-to-valley differences of electricity, gas, and heat loads are reduced by 6.51%, 21.28%, and 7.13%, respectively, so IDR has the effect of suppressing load peaks and "cutting peaks and filling valleys".

The impact of Electric Vehicles on GES
According to Table 2, EV mainly changes the electric power balance of IES. Therefore, by comparing and analyzing the power dispatching situations of Scenario 5 and Scenario 6, the further effect of EV on GES can be seen. Figure 9 shows the EV schedule in different ways, as shown in Figure A2, which is the electric balance diagram of Scenario 5. same time, the heat load demand at night is large, the CHP units and the P2H output are large, and the TS stores heat during this period. The gas load demand is met by gas purchases, GS, and MR, and the hydrogen energy required by MR is provided by EL and HS. ES, EV, and HFC output electric energy, and GT output gradually increases during the peak periods of electricity prices from 06:00 to 12:00 and 19:00 to 22:00. As a result, the amount of electricity purchased is gradually decreasing. The heat load demand is reduced, and the heat is mainly provided by CHP and GB. Meanwhile, TS and HFC produce thermal energy. The gas purchase of IES increases, but the gas price is high during this time, so GS mainly releases energy. In 13:00 to 18:00, which is the normal period for electricity prices, the proportion of photovoltaic power generation in IES gradually increases, and ES stores excess electric energy. At this time, the heat load demand is reduced, and the gas purchase amount is also reduced. GS stores gas from 13:00 to 14:00, which is the normal period for gas prices. As can be seen from Figure 8, the state of charge (SOC) of ES and HS tends to be the same because, in periods of low electricity prices, the wind power output is large, so both ES and HS are stored energy. During the peak period of energy supply, ES discharges. Meanwhile, the HFC outputs electric and thermal energy, and the hydrogen consumption is provided by HS. In summary, the IES scheduling model proposed in this paper reduces energy purchase costs and makes multi-source energy storage equipment play a role in peak regulation. As shown in Figure A1, the optimization curves of electricity, heat, and gas loads before and after IDR are shown. Through IDR, the electricity and gas load in the peak period of electricity and gas prices are transferred to the low period and normal period to reduce the operating costs of IES. At the same time, the heat load with higher power at night is partially transferred to daytime by IDR. After IDR optimization, the peak-to-valley differences of electricity, gas, and heat loads are reduced by 6.51%, 21.28%, and 7.13%, respectively, so IDR has the effect of suppressing load peaks and "cutting peaks and filling valleys".

The impact of Electric Vehicles on GES
According to Table 2, EV mainly changes the electric power balance of IES. Therefore, by comparing and analyzing the power dispatching situations of Scenario 5 and Scenario 6, the further effect of EV on GES can be seen. Figure 9 shows the EV schedule in different ways, as shown in Figure A2, which is the electric balance diagram of Scenario 5. As is shown in Figures 9 and A2, compared with the electric power balance diagram of Scenario 6 in Figure 7a, the peak period of disordered EV charging is from 15:00 to 22:00, during which the demand for electric load is large and the electricity price is high. As is shown in Figures 9 and A2, compared with the electric power balance diagram of Scenario 6 in Figure 7a, the peak period of disordered EV charging is from 15:00 to 22:00, during which the demand for electric load is large and the electricity price is high. The disordered EV charging load increases the power load. IES needs additional power purchases and increases the output of CHP units to meet the power load demand. It can be seen from Table 2 that the power purchase cost of Scenario 5 is 1148.60 yuan higher than that of Scenario 6. EVs, as generalized energy storage devices, operate charging and discharging during valley and peak periods of electricity prices, respectively, reducing the output of coal-fired units and GTs and reducing carbon emissions by 334.02 kg. By adjusting the number of EVs, the impact of EVs on GES is further analyzed, as shown in Figure 10. The disordered EV charging load increases the power load. IES needs additional power purchases and increases the output of CHP units to meet the power load demand. It can be seen from Table 2 that the power purchase cost of Scenario 5 is 1148.60 yuan higher than that of Scenario 6. EVs, as generalized energy storage devices, operate charging and discharging during valley and peak periods of electricity prices, respectively, reducing the output of coal-fired units and GTs and reducing carbon emissions by 334.02kg. By adjusting the number of EVs, the impact of EVs on GES is further analyzed, as shown in Figure 10.  As shown in Figure 10a, the IES operating cost increases as the number of EVs increases. Meanwhile, the difference between the operating cost of EV ordered charging and disordered charging increases gradually, and the advantage of EV ordered charging and discharging is more obvious. When the number of EVs rises from 60 to 70, the IES operating cost increases substantially in both EV operating modes. According to Figure 10b, the carbon emissions of different EV quantities also have the same trend. Therefore, by charging and discharging EVs in an orderly manner and in the range of 60 units, IES can ensure efficient operation of each unit and realize optimum economy.

Uncertainty Analysis Considering CVaR
Combined with the new energy output and load power forecast data for the region, 2000 random scenarios are generated by LHS, and 100 typical scenarios are obtained after scenario reduction. To consider the effect of source-load uncertainty, β is taken as 0.90, and the relationship curve between IES expected costs and CVaR values under different risk preference coefficients is obtained by combining typical scenarios, as shown in Figure 11.  As shown in Figure 10a, the IES operating cost increases as the number of EVs increases. Meanwhile, the difference between the operating cost of EV ordered charging and disordered charging increases gradually, and the advantage of EV ordered charging and discharging is more obvious. When the number of EVs rises from 60 to 70, the IES operating cost increases substantially in both EV operating modes. According to Figure 10b, the carbon emissions of different EV quantities also have the same trend. Therefore, by charging and discharging EVs in an orderly manner and in the range of 60 units, IES can ensure efficient operation of each unit and realize optimum economy.

Uncertainty Analysis Considering CVaR
Combined with the new energy output and load power forecast data for the region, 2000 random scenarios are generated by LHS, and 100 typical scenarios are obtained after scenario reduction. To consider the effect of source-load uncertainty, β is taken as 0.90, and the relationship curve between IES expected costs and CVaR values under different risk preference coefficients is obtained by combining typical scenarios, as shown in Figure 11.
As shown in Figure 11, the corresponding CVaR values of IES operating costs decrease as the risk preference coefficients δ increase, and the IES expected operating cost and the CVaR value can be obtained from Equations (34)-(36) on the basis of the generated scenarios. When δ is small (δ takes 0.1, 0.2), it represents that IES operators are more aggressive investors, preferring to take greater operational risks in order to obtain lower expected costs. When δ is large (δ takes 0.9, 0.8), it represents that the investment preference of IES operators is conservative, raising the expected costs to obtain smaller risk losses. The inflection point of the relationship curve occurs when δ = 0.4. If δ is reduced at this point, the expected costs will be stabilized. Therefore, choosing the appropriate risk preference coefficients is beneficial for IES to reduce the operating risks and the operating costs of the system.  As shown in Figure 11, the corresponding CVaR values of IES operating costs d crease as the risk preference coefficients δ increase, and the IES expected operating co and the CVaR value can be obtained from Equations (34)-(36) on the basis of the generat scenarios. When δ is small ( δ takes 0.1, 0.2), it represents that IES operators are mo aggressive investors, preferring to take greater operational risks in order to obtain low expected costs. When δ is large ( δ takes 0.9, 0.8), it represents that the investment pr erence of IES operators is conservative, raising the expected costs to obtain smaller ri losses. The inflection point of the relationship curve occurs when 0.4 δ = . If δ is r duced at this point, the expected costs will be stabilized. Therefore, choosing the appr priate risk preference coefficients is beneficial for IES to reduce the operating risks and t operating costs of the system.

Conclusions
This paper improves the low-carbon and economic performance of multi-energy co pled IES. Kalina cycles and GES are introduced on the supply side and demand side, r spectively, to participate in optimal dispatching, and the basic scheduling model is esta lished. Meanwhile, considering the influence of source-load uncertainty, a multi-ener complementary optimal scheduling model of hydrogen-containing IES based on CV was established. After analysis, the following conclusions are obtained: (1) In this paper, the hydrogen energy utilization model is refined, and the waste heat the P2G methanation reaction is recovered and coupled with the Kalina cycle to furth improve the thermoelectric decoupling capacity of CHP units. While increasing t consumption rate of new energy to 93.88%, the operation cost was reduced by 996. yuan, effectively improving the energy utilization rate and operation economy of IE (2) In the GES, IDR and EV can cooperate with the actual energy storage equipment adjust the load peak-to-valley difference of the system. After IDR optimization, t peak-to-valley differences of electricity, gas, and heat loads are reduced by 6.51 21.28%, and 7.13%, respectively, realizing the "peak cutting and valley filling" of t

Conclusions
This paper improves the low-carbon and economic performance of multi-energy coupled IES. Kalina cycles and GES are introduced on the supply side and demand side, respectively, to participate in optimal dispatching, and the basic scheduling model is established. Meanwhile, considering the influence of source-load uncertainty, a multienergy complementary optimal scheduling model of hydrogen-containing IES based on CVaR was established. After analysis, the following conclusions are obtained: (1) In this paper, the hydrogen energy utilization model is refined, and the waste heat of the P2G methanation reaction is recovered and coupled with the Kalina cycle to further improve the thermoelectric decoupling capacity of CHP units. While increasing the consumption rate of new energy to 93.88%, the operation cost was reduced by 996.31 yuan, effectively improving the energy utilization rate and operation economy of IES. (2) In the GES, IDR and EV can cooperate with the actual energy storage equipment to adjust the load peak-to-valley difference of the system. After IDR optimization, the peak-to-valley differences of electricity, gas, and heat loads are reduced by 6.51%, 21.28%, and 7.13%, respectively, realizing the "peak cutting and valley filling" of the load. Although the IES needs to pay 1343.23 yuan for the user's virtual energy storage scheduling, the total cost is reduced by 3407.47 yuan. Therefore, virtual energy storage improves the users' sensitivity to energy prices, reduces IES operating costs, and achieves an all-win situation for both the users and the IES operators. (3) CVaR is used to measure the risk losses from uncertainty in IES. IES operators can choose appropriate risk appetite coefficients to reduce the risk losses and system operation costs of IES. The IES scheduling model established in this paper provides a strategy for IES to achieve refined utilization of hydrogen energy and low-carbon economic operation. However, due to the influence of users' travel and consumption habits, both EV travel and user response behavior have uncertainty and randomness. In the future, further research will be conducted to address the uncertainties surrounding integrated demand response and EV travel.

Conflicts of Interest:
The authors declare no conflict of interest.