Low-carbon optimal operation of the integrated energy system considering integrated demand response and oxygen-rich combustion capture technology

In view of the current operating constraints and environmental pollution problems of traditional units, in this article, oxygen-rich combustion capture technology is introduced to transform gas-ﬁ red units, demand response technology is used on the load side, the energy conversion equipment such as power-to-gas equipment is combined to form an integrated energy system, and then, a low-carbon optimization approach of the integrated energy system is proposed. First, the system architecture is constructed, and a model with an oxygen-rich combustion unit and integrated demand response is established. Second, a power-to-gas equipment model considering reaction waste heat utilization and oxygen recovery is established. Finally, a stepped carbon trading mechanism is introduced to establish a low-carbon economic scheduling model for the integrated energy system with the goal of minimizing the operating cost of the integrated energy system. The simulation results show that the total cost and carbon emissions of the integrated energy system are reduced by 6.44% and 44.24%, respectively, under this model. At the same time, the operation adjustment capability and the oxygen production ef ﬁ ciency of the internal units of the system are improved.


Introduction
The issues of environmental pollution and energy shortage have accelerated the global development trend of low-carbon emission and clean energy.In order to achieve the grand development goal of "carbon peaking and carbon neutrality," the complementary and coordinated optimization of an integrated energy system (IES) has become a key technical direction for addressing the contradiction between energy supply and demand and improving energy conversion capabilities.Research on its optimized operation has become a hot topic (Han et al., 2021).
At present, the main measures taken to achieve the low-carbon IES can be classified into two categories: policy, such as establishing carbon trading mechanisms, and technology, such as carbon capture technology and demand-side response technology.The low-carbon economic operation of the entire system can be achieved by coordinating the above two main measures.
Currently, the main research method for achieving low-carbon scheduling is to introduce carbon trading mechanisms into the IES.Qin et al. (2018) introduced carbon trading mechanisms into the electricity-gas-heat integrated energy system, which can be conducive to the low-carbon economic operation of the system, promoting the development of low-carbon units.Zhou et al. (2018) introduced carbon trading mechanisms into the IES, and the results showed that carbon trading mechanisms can effectively promote the consumption of new energy, improving the low-carbon and economic efficiency of the IES.Wang et al. (2019) used the lifecycle analysis method to analyze the carbon emissions generated by different energy chain migration and transformation processes in the IES.The results showed that reasonable planning of the carbon market can promote the low-carbon development of the IES.
In the above research, the carbon trading mechanism can only function on gas-fired units within the IES.As thermal power serves as the external supply base load of the IES, the carbon trading mechanism cannot be effectively optimized.Therefore, the coordination of technology and policy should be considered.Lu et al. (2018) established a minimum coordination relationship model between the heat storage and release rate and the electric boiler power under the condition of extreme wind power consumption.The results showed the effectiveness of the heat storage and coordinated heating relationship between extreme wind power consumption and the electric boiler.Gao et al. (2019) introduced power-to-gas (P2G) carbon capture power plants into the IES.The results showed that power-to-gas-carbon capture power plants can effectively improve the absorption capacity of wind power and reduce system carbon emissions and operating costs.He et al. (2018) analyzed the flexible operation modes of carbon capture systems and P2G equipment in the IES.The results showed that the proposed model and flexible operation mode can effectively reduce carbon emissions and operating costs of integrated energy systems.
It is worth mentioning that currently, limited research exists on the introduction of carbon capture technology in the IES.Among them, some scholars have adopted post-combustion capture technology for the low-carbon transformation of gas turbines, but there are shortcomings such as a huge covering area for the equipment, weak capture capacity, and high capture cost, which are suitable for thermal power units (Liu et al., 2010).The precombustion capture technology requires significant modifications to the unit, resulting in low applicability.The oxygen-rich combustion capture technology combines the advantages of pre-and postcombustion capture technologies, with little impact on the power generation process of the units and strong carbon capture ability.However, it has high requirements for the cleanliness of the fuel itself and is suitable for gas-fired units, with good development prospects (Dai et al., 2023).
Current research at home and abroad is mainly focused on the economics and operating characteristics of oxygen-rich combustion technology.Gao et al. (2014) analyzed the impact of changes in unit operating parameters and different flue gas re-circulation methods on air separation oxygen generation equipment and gas capture equipment and then established a calculation model for the operating energy consumption of the related equipment.Cui et al. (2021) established a dynamic cycle model and a dynamic change model of flue gas composition for oxygen-rich combustion flue gas and analyzed the energy consumption of power plants at different oxygen concentrations in oxygen-rich combustion scenarios.Oboirien et al. (2014) analyzed the energy consumption and economic performance of oxygen-rich combustion equipment under different parameters.Koiwanit et al. (2014) established a full process model of the oxygen-rich combustion system and proposed a utilization plan for residual heat based on thermodynamic principles to improve the operational economy of the system.Zaharia et al. (2019), Yang et al. (2022), andGuo et al. (2023) conducted an economic analysis of the investment cost, operating cost, and CO 2 emission reduction cost of oxygen-rich combustion power plants.
With the continuous enrichment and expansion of resources on the demand side, considering only oxygen-rich combustion technology and carbon trading mechanisms cannot meet the needs of the low-carbon operation of the IES.At present, utilizing demand-side resources for carbon reduction has become a research hotspot, but most of it focuses on the demand-side resource interaction from the perspective of "electricity."Zhou et al. (2018) and Ji et al. (2022) introduced the "electricity price" demand response (DR) model into the IES, and the results showed that the proposed model can effectively reduce carbon emissions and operating costs of the IES.Cui et al. (2022) proposed a new carbon reduction mechanism for power systems, which guides users to actively respond and reduce carbon emissions from a "carbon perspective" (Zhou et al., 2018;Malehmirchegini et al., 2022).Overall, whether from the perspective of users or the system, demand response mechanisms can bring carbon reduction benefits.
In summary, this paper proposes an optimal operation scheduling model of the IES that considers integrated demand response and oxygen-rich combustion capture technology.The main contributions of this paper are as follows: 1.An oxygen-rich combustion technology is introduced to transform the gas-fired unit, and a model of the oxygen-rich combustion capture unit is established for solving the operational constraints and environmental pollution problems of traditional units within the IES. 2. An integrated demand response model of price and subsidy incentives is constructed based on the price elasticity matrix principle and response subsidy policy to provide full play to the adjustment ability of demand-side response resources.3. A stepped carbon trading mechanism is proposed, the impact of the base price of the carbon trading mechanism on system costs and carbon emissions is discussed, and the low-carbon operation model of the IES considering the stepped carbon trading mechanism is constructed.
2 Architecture of an integrated electricity-heat-gas-cooling energy system The structure of the IES is shown in Figure 1.Here, the gas-fired unit adopts the oxygen-rich combustion capture technology for lowcarbon transformation; electric power is supplied by the combined heat and power unit, wind power, and purchased electricity; the heat load is supplied by the combined heat and power unit and the residual heat from the P2G converter; the gas load is supplied by the external gas network and the P2G converter; the cold load is supplied by the air-conditioning and refrigeration, the residual heat from lithium bromide, and the liquefied natural gas tanks in the gasification process to supply the loads; and the oxygen required by the oxygen-rich combustion capture unit is supplied by the P2G converter, air oxygen generator, and oxygen storage tank.Oxygen for the oxygen-rich combustion capture unit is supplied by electric gas conversion equipment, air oxygen generating equipment, and oxygen storage tanks (Kang et al., 2022).

2.1
The modeling and combined heatingsupply principle of the oxygen-rich combustion unit

Modeling of oxygen-rich combustion capture units
Oxygen-rich combustion capture technology increases the concentration of CO 2 in flue gas after combustion using highpurity oxygen or pure oxygen, which can be helpful for purification or storage after capture.The energy flow diagram of the oxygen-rich combustion unit is shown in Figure 2.
The main energy destinations of the oxygen-rich combustion capture unit are the system electrical load, carbon capture equipment, air separation oxygen generation equipment, and system heat load (Zhu et al., 2022).
Oxygen-rich combustion capture technology has a higher degree of fuel cleanliness, and back-pressure gas-fired units are selected as the object of transformation.The oxygen-rich combustion capture units in the following paragraph refer to the gas-fired units transformed by the oxygen-rich combustion capture technology.The relationship between electricity generation P GU t and natural gas consumption of a gas-fired unit V GU t is shown as follows: where η GU is the gas-electricity conversion efficiency of the gasfired unit; L CH4 is the calorific value of methane, which is usually taken as 36 MJ/m 3 ; r is the conversion rate of the calorific value and power, which is usually taken as 3.6%; P GU max and P GU min are the upper and lower limits of the gas-fired unit, respectively; and ΔR GU d and ΔR GU u are the upper and lower limits of climbing constraints in the gas-fired unit, respectively.
The heating power of the oxygen-rich combustion capture unit is the heating power of the gas-fired unit, and the electricitygenerating power and heat-producing power of the back-pressure Integrated energy system structure.

Frontiers in Energy Research
frontiersin.org unit present a linear mathematical relationship in an approximate manner, with the specific expression as follows: where E OCC t is the heat-producing power of the oxygen-rich combustion capture unit at time t and A GU is the heat-to-power ratio of the gas-fired unit.
The introduction of the oxygen-rich combustion capture system will change the structure of the gas-fired unit to cause part of the energy loss, known as the baseline energy consumption P OCC B , which does not change with the change in the operating status of the unit and can be considered a fixed value.In addition, the unit needs to consume part of the electrical energy to maintain the operation of the oxygen-rich combustion capture system, and its operational energy consumption mainly includes the air separation oxygen generation equipment and carbon capture equipment; then, the energy consumption of the oxygen-rich combustion capture system P OCC t is shown as follows: where t 1, 2, 3..., T is the operational phase, the value of which is taken as 24; P CCE t and P ASO t are the power consumption of the carbon capture equipment and the air separation oxygen generation equipment at time t, respectively; m CCE t is the CO 2 capture mass of the carbon capture equipment at time t; λ CCE and λ ASO are the unit operating energy consumption of the carbon capture equipment and the air separation oxygen generation equipment, respectively; β CCE t and e G are the carbon capture level of the carbon capture equipment and the carbon emission intensity of the gas turbine at time t, respectively; and O ASO t is the amount of oxygen produced by the air separation oxygen generation equipment at time t.
The oxygen consumption O OCC t in the oxygen-rich combustion capture unit t is shown as follows: where α OCC is the oxygen consumption per unit of the output power when the oxygen-rich combustion capture unit is in operation.Zhang et al. (2022) proved the process flow and equipment operation steps of traditional air technology to be basically the same as those of oxygen-rich combustion capture technology.Technically, the two operation modes can be converted to each other, and the switching time is usually less than 1/3 h.Thus, the influence of different modes of conversion is ignored.The restart time of the air separation oxygen generation equipment is longer with larger start-up energy consumption and minimum operating power.Therefore, the carbon capture equipment will stop running in the air operation mode, the unit exits the oxygenrich combustion operation state, and the oxygen produced by the air separation oxygen generation equipment is stored in the oxygen tank.Then, the oxygen-rich combustion capture system collectively consumes energy as follows: where a ASO min and P ASO max are the minimum operating coefficient and the upper limit of the electrical power consumption of the air separation oxygen generation equipment, respectively.Energy flow diagram for oxygen-rich combustion capture units.

Frontiers in Energy Research
frontiersin.org The net output power of the oxygen-rich combustion capture unit P OCC t,N represents the overall external power generation of the unit.The net output power at time period t is shown as follows: The CO 2 generated from the operation of the unit is mostly absorbed by the carbon capture equipment, and a small amount of it is emitted into the atmosphere.The net carbon emissions from the oxygen-rich combustion capture unit at time t are shown as follows: Zhu et al. ( 2022) defined the correlation between the net external power generation output and the net carbon emissions of a carbon capture unit as the "electro-carbon characteristic" of the unit, which is expressed as follows for an oxygen-rich combustion capture unit: where β CCE max is the limiting capture level of the carbon capture equipment.
When the oxygen-rich combustion capture unit is in operation, the CO 2 captured by the carbon capture equipment can be used either as a reaction feedstock for the power-to-gas equipment or for carbon sequestration, which can be expressed as follows: where m CCE t,P2G is the mass supplied by the carbon capture equipment to the power-to-gas equipment at time t and m CCE t,CS is the carbon mass sequestered at time t.
Based on the above analysis, the upper and lower limits of the net output power of the oxygen-rich combustion capture unit are introduced as follows: In addition, the following power constraints need to be met when operating the air separation oxygen generation equipment: Oxygen storage tanks can achieve the use of oxygen across different time periods and store the oxygen in the form of liquid.Regarding the storage or discharge period, the form of oxygen is usually gaseous.In addition, the model for oxygen storage tanks in the gaseous form is presented as follows: In Figure 3, assuming that the amount of oxygen produced by the air separation oxygen generation equipment at point N exactly meets the demand for oxygen-rich operation, when approaching from point N to point H, the carbon emission increases, the air separation oxygen generation equipment that maintains the minimum power cannot meet the demand for oxygen-rich unit operation, and the difference in oxygen can be supplied by the oxygen storage tank.Therefore, ABCDNGE is the oxygen-rich operation area of the unit when using the oxygen storage tank.At this time, the phenomena of discarding oxygen occur.Pentahedron NHGEF presents the extra oxygen-rich operation area in the unit equipped with the oxygen storage tank, which indicates that when the net output of the oxygenrich combustion capture unit is larger, the energy consumption of the air separation oxygen equipment is reduced, and the oxygen storage tank can be used to supplement for the difference in oxygen.When the net output of the unit is smaller, the energy consumption of the air oxygen equipment can be improved to add oxygen into the oxygen storage tank.The analysis shows that the oxygen-rich combustion capture unit equipped with the oxygen storage tank can expand the scope of operation.
Therefore, the operation and maintenance (O&M) cost of the oxygen-rich combustion capture unit is shown as follows: Electric carbon characteristics of oxygen-rich combustion capture units.
where s ASO yw 、s CCE yw 、s GU yw 、s OT yw are the O&M cost factors for the air separation oxygen generation equipment, carbon capture equipment, gas-fired units, and oxygen storage tanks, respectively.
Detailed modeling of the internal equipment in the IES can be found in Supplementary Material.Air liquid energy storage architecture is shown in Figure 4.
3 Low-carbon strategies based on lowcarbon demand response DR resources mainly include price-based demand response (PDR) and incentive-based demand response (IDR).PDR guides users to engage in reasonable electricity consumption behavior by changing electricity prices, thereby adjusting electricity consumption plans.IDR includes interruptible loads, demand-side bidding, emergency demand-side response, and direct load control.The power department usually first signs a contract with the load agency (aggregator) to classify and integrate various users who can improve DR resources, and finally, it is uniformly regulated by the dispatch center (Tian et al., 2023).
DR resources can optimize the output plans of units with different carbon emission intensities by changing the load curve (Zhang et al., 2021), thereby achieving the effect of reducing system carbon emissions.The low-carbon principle of price-based demand response is shown in Figure 5.
As shown in Figure 5, the demand response can shift some load from the peak period to the valley period.During the peak period, when the net output of low-carbon units has reached the upper limit, the output of high-cost and high-carbon emission thermal power units supplies the difference in the load.After the price-based demand response, this part of the load is supplied by the output of low-cost, low-carbon emission carbon capture units during the valley period, or additional consumption of low-marginal cost, carbon-neutral wind power generation, which effectively reduces the carbon emissions of the system.

Price-based demand response models
The electricity load demand response model in this article adopts the price demand response elasticity matrix method.The electricity load change rate and electricity price change rate can be characterized by the elasticity index of electricity consumption and electricity price, and its expression is as follows: where m is the elasticity index of electricity consumption and electricity price; c p and Δc p are the differences in the electricity price in the peak-valley period and the fixed electricity price; and L and ΔL are the amount of electricity consumption and load response before the load response, respectively.
Based on the ratio of the time-of-use electricity price to the fixed electricity price, the following elasticity matrix is established: where N ii is the self-elasticity coefficient; N e is the elasticity matrix of electricity consumption and electricity price; L i and ΔL i are the amount of electricity and its variation for the user in time period i, respectively; N ij is the cross-elasticity coefficient; and c pi and Δc pi and c pj and Δc pj are the electricity price and its variation in time periods i and j, respectively.
Based on the electricity elasticity matrix M e , the customer load response is obtained as follows: L e,0,f 0 0 where L e,0,f , L e,0,p , and L e,0,g denote the electricity consumption in the peak, valley, and flat periods before the price-based incentive, respectively; ΔL e (t) is the amount of electricity use transferred after the price-based incentive; and c e,p0,f , c e,p0,p , and c e,p0,g and Δc e,pf , Δc e,pp , and Δc e,pg denote the difference between the fixed electricity price and time-of-use electricity price, respectively.
Since gas energy has the same commodity property as electricity, the gas load can be adjusted based on the published time-of-day gas price, so the IDR model of the gas load can also be modeled by the method of the price-type elasticity matrix.Analogous to the electricity price-based IDR model, the relationship between the gas load and the time-of-day gas price also exists, as shown in Eq. ( 29), which is expressed as follows: where ΔL g (t) is the amount of the gas load shifted after the price-based incentive for the gas; L g,0,f , L g,0,p , and L g,0,g denote the amount of gas used in the peak, valley, and flat periods before the price-based incentive, respectively; c g,p0,f , c g,p0,p , and c g,p0,g and Δc g,pf , Δc g,pp , and Δc g,pg represent the difference between fixed and time-of-day gas prices, respectively; and N g denotes the elasticity matrix of gas price and gas consumption.

Incentive-based demand response
Incentive-based demand response refers to the direct reduction or interruption of a portion of the load by users based on compensation and incentive mechanisms established by energy operators during the peak period of load or in emergency situations.The incentive demand response constructed in this article includes electricity load, gas load, heat load, and cooling load.

Incentive-based electricity load and gas load
When the system is in peak operation, the electricity/gas consumer has a prior agreement with the energy operator to curtail part of the load, and its mathematical model is shown as follows: where ε e and ε g are reduction allowance factors for electricity and gas loads, respectively; L e,cut (t) and L g,cut (t) are the amount of electricity and gas loads curtailed by the user, respectively; and E IDR,e (t) and E IDR,g (t) are the costs of incentive subsidies for users to curtail their electricity and gas loads, respectively.
For the reduction in electricity and gas loads, in order to prevent too much impact on the life of users, the amount of adjustment of electricity and gas loads at moment t should be within a certain range, and the total load change within 1 day should also meet the limit as follows: where L e,0 (t) and L g,0 (t) are the initial electricity load and initial gas load, respectively; φ e,cut and φ g,cut are the upper limits of the reduction rates for electricity and gas loads, respectively; and φ e total and φ g total are the upper limits of the total reduction rates for electricity and gas loads, respectively.

Modeling of heat/cooling load
Considering the time delay and perceptual ambiguity of heat and cooling loads, adjusting the temperature within the comfort range will not have a significant impact on users.The quantitative relationship between the heat/cooling load demand and indoor and outdoor temperatures is shown as follows: where ε is the indoor heat loss per unit of the floor area under the condition of a difference in temperature; S is the indoor area of a building; C is the specific heat capacity per unit of the floor area; T h,in (t) and T h,out (t) are the indoor and outdoor temperatures of the heating system, respectively; and T c,in (t) and T c,out (t) are the indoor and outdoor temperatures of the cooling system, respectively.
The mathematical model of the demand response for the heat load and cooling load can be expressed as follows: where ΔT h,in (t) and ΔT c,in (t) are the outdoor temperature changes for heating and cooling systems, respectively; ΔL h (t) and ΔL c (t) are the changes in the heat and cooling loads of the user at time t, respectively; ΔL h, max and ΔL c, max are the upper limits of the amount of variation in the heat and cooling loads of the users, respectively; and T max and T min are the upper and lower temperature limits where the user is in the comfort range, respectively.
Due to the time delay and ambiguity of user requirements for indoor temperature comfort, the impact of small changes in the heat/cooling load on user comfort is relatively small.However, when the heat/cooling load changes significantly, the impact on user comfort is significant, so the impact of the heat/cooling load change on user comfort is not linearly correlated.Therefore, this article adopts a stepped compensation method to encourage users to change the heat/cooling load.The compensation cost for the user heat/cooling load at time t is shown as follows: where E h,IDR (t) and E c,IDR (t) are the cost of incentive-based IDR subsidies for heat and cooling loads, respectively and ε h,IDR and ε c,IDR are the subsidy factors for changes in heat and cooling loads, respectively.
Taking the heat load incentive compensation as an example, the greater the deviation degree of heat load, the greater the impact on user comfort.This article adopts a stepped subsidy coefficient, which varies in a stepped manner based on the impact of the actual load on user comfort.The greater the deviation degree, the greater the subsidy coefficient, which can be expressed as follows: where λ h,IDR is the penalty factor of the subsidy coefficient for the heat load changes and α 1 and α 2 are the dividing boundaries for the amount of the heat load change, respectively.
Similarly, the subsidy coefficients for changes in cooling loads can be expressed as follows: where λ c,IDR is the penalty factor of the subsidy coefficient for the cooling load changes and β 1 and β 2 are the dividing boundaries for the amount of the cooling load change, respectively.

Stepped carbon trading models
At present, the domestic carbon emission trading market has started trial operation, guiding various industries to achieve lowcarbon emissions through policies.The main methods of allocating carbon emission reduction quotas include free allocation, paid allocation, and mixed allocation.Currently, free allocation is mainly used in China.The stepped carbon trading mechanism mainly consists of three parts: initial carbon emission quotas, actual carbon emissions, and carbon trading costs (Guan et al., 2018;Chrispim, 2021;Liang et al., 2021).

Modeling of carbon emission quotas
The main source of carbon emissions in the system consists of two components, namely, the oxygen-rich combustion unit and the equivalent emissions from power purchases.The model of carbon emission quotas is shown as follows: where G IES t,pe , G buy t,pe , and G OCC t,pe are carbon emission quotas for the IES, purchased electricity, and oxygen-rich combustion units, respectively; g p and g h are the carbon emission intensity per unit of electrical and thermal power, respectively, which are usually taken as 0.728 kg CO 2 /(KWh) and 0.102 kg CO 2 /MJ; P buy t is the amount of purchased electricity of the system at time t; and g p−h is the electric heat conversion factor, which is usually taken as 6 MJ CO 2 /(KWh).

Modeling of actual carbon emissions
In this system, it is assumed that all purchased electricity comes from coal-fired units and that the P2G converter in the system will reduce some of the carbon emissions.Therefore, the model of the actual carbon emissions is shown as follows: where G IES t,pf , G buy t,pf , and G buy t,pf are the actual carbon emissions of the IES and purchased electricity, respectively; a 1 、b 1 、c 1 are the carbon emission factors for coal power; and a 2 、b 2 、c 2 are the carbon emission coefficients for natural gas energy supply.

Modeling of the stepped carbon transaction cost
In order to ensure a reasonable control system of carbon emissions, this paper adopts a stepped carbon transaction cost model.The stepped pricing mechanism divides multiple purchase intervals, and the more carbon emission rights to be purchased, the higher the purchase price of the corresponding interval.The stepped carbon transaction cost model can be expressed as follows: where F IES t,CO2 is the systematic carbon transaction cost at time t; λ is the carbon trading base price; l is the length of the carbon emission interval; and α is the growth rate of carbon trading prices.
5 Optimization model of the integrated energy system

Objective function
The model proposed in this article comprehensively considers the cost of purchasing natural gas from the system, the cost of purchasing electricity from the power grid, the cost of operation and maintenance, the cost of carbon storage, the cost of carbon trading, the cost of wind abandonment penalty, and the cost of load loss penalty, with the objective function of minimizing the total operating cost for optimization.The objective function is shown as follows: where F IES t,buy is the cost of purchasing energy for the system at time t; F IES t,wh is the cost of operating and maintaining the system at time t; F CO2 t,CS is the cost of carbon sequestration at time t; and F FD t,cf and F IES t,fh are the system wind abandonment penalty cost and the electric load loss penalty cost at time t, respectively.

Cost of the purchased electricity
The system purchased energy cost includes two components: the purchased electricity cost and purchased gas cost.In addition, the price of the purchased electricity is decided by using the time-of-day electricity price mechanism: is the price at which the system purchases electricity from the grid at time t and s buy t,gas and V buy t are the price and volume of gas purchased by the system at time t, respectively.are the O&M costs for the P2G equipment, wind power, liquefied natural gas (LNG) tanks, lithium bromide absorption refrigerant, and electric refrigeration and air conditioning, respectively.

Cost of carbon sequestration
Carbon sequestration costs are required to sequester the remaining CO 2 when the P2G equipment is unable to consume the captured CO 2 , and the cost of carbon sequestration can be expressed as follows: where s CO2 CS is the cost coefficient of carbon sequestration.

Cost of wind abandonment penalty and load loss penalty
In order to make the system realize the full consumption of new energy and reduce the occurrence of the wind abandonment phenomenon, the wind abandonment penalty cost is used to improve the wind power consumption rate of the system.The cost of the wind abandonment penalty is shown as follows: where s WP cf is the unit cost of wind abandonment penalties; P yc t and P WP t denote the predicted wind power output and wind power consumption power at time t, respectively; s WP fh is the unit cost of load loss penalties; and P yc,Eload t and P Load t denote the predicted power of the load and the actual electrical load at time t, respectively.

Cost of the demand response
5.2 Constraints 5.2.1 System balance constraints are the discharge power of the energy storage system, the electrical energy released by the LNG gasification station in the gasification of natural gas, the power of electrical energy absorbed by the energy storage system, and the power of the system electrical load loss, respectively.
• Gas power balance constraints where V LNG t and V load t are the natural gas power released from the LNG gasification station and the system natural gas load, respectively.
• Heat power balance constraint where E Load t is the heat load of the system and E XHL t is the heat power consumed by the residual heat of lithium bromide refrigeration.
• Cold power balance constraint where , and L Load t are the cold power released from the LNG gasification station, the cold power released from the energy storage unit, air-conditioning refrigeration, lithium bromide refrigeration, and the system cold load, respectively.
• Oxygen balance constraint where O loss t is the amount of oxygen discarded at time t.
• Wind power operation constraint

LAES constraint
The liquid air energy storage (LAES) and discharge constraints and storage capacity constraints are shown as follows: Comprehensive energy system load and wind power forecasts.
This article constructs a complete low-carbon economic optimization model for the IES, taking into account the system energy purchase cost, operation and maintenance cost, carbon trading cost, carbon sequestration cost, demand response cost, wind abandonment penalty, and load loss penalty.The optimal output plan for each unit and equipment with the minimum sum of multiple costs during the scheduling cycle can be obtained, resulting in the lowest total operating cost of the system.Power balance diagram of the integrated energy system.

Parameters of the algorithm
The structure of the algorithm is shown in Figure 1.The carbon emission intensity of the gas turbine is 1.964 kg/m3 , the climbing rate is 10,000 kW/h, and the operation and maintenance cost is 0.02 yuan/KW; the maximum storage/discharge rate and efficiency of the oxygen storage tank are 2,500 m 3 /h and 95%, respectively, the maximum oxygen storage capacity is 10,000 m 3 /h, the dissipation coefficient is 0.02%, and the operation and maintenance coefficient is 0.02 yuan/m 3 ; the base price of carbon emission is taken as 200 CNY/ton, the interval length is taken as 2t, and the growth rate is taken as 0.25; the operation and maintenance cost of the P2G equipment is taken as 0.02 yuan/ kWh; the penalty for wind abandonment of wind power is 0.6 CNY/KWh, and the operation and maintenance cost is 0.2 yuan/kWh; and the cost of the load loss penalty is 10 CNY/ kWh. Figure 6 shows the power load prediction curve of the IES, which contains the power of each load and the predicted power of wind power.Table 1 shows the time-of-use electricity price and time-of-use natural gas price; Table 2 shows the parameters of the IES equipment; Table 3 shows the basic parameters of oxygen-rich combustion capture technology.In the calculation analysis, t means the time period of t in the scheduling day.For example, when t = 1, it means the first scheduling time period of the day.

Analysis of basic operational results
In order to verify that the proposed strategy can satisfy the system supply and demand balance, the operational results obtained from the proposed model are analyzed.Carbon dioxide emissions from integrated energy systems.Power-to-gas (P2G) device operation results.

Scenario
Demand response Oxygen-rich combustion capture technology P2G residual heat and pure oxygen recovery Figure 7A shows that the IES only purchases electricity from the grid at the time periods of T = 11, 13, 14, and 18, and the electric load at other times mainly comes from the oxygen-rich combustion capture unit and wind power, which indicates that the IES has high self-sufficiency in electric energy.Figure 7B shows that the residual heat from the reaction of the oxygen-rich combustion unit and the P2G equipment can meet the system heat load.Figure 7C shows that the gas load mainly consists of the upper gas network, liquefied natural gas, and natural gas produced by machination, which can meet the gas turbine and system gas load.Similarly, Figure 7D shows that the output of each unit in the system can meet the system cooling load.Load comparison and after integrated energy system demand response.
Frontiers in Energy Research frontiersin.org14 Ji et al. 10.3389/fenrg.2024.1296709 Figure 8 shows that most of the CO 2 generated by the IES is sequestered after the carbon quota, while the other part is used as raw material for the reaction of the P2G equipment, which reduces the cost of carbon trading.In addition, the natural gas generated from it can be used as a raw material for gas turbine units, achieving internal carbon resource circulation and improving the clean efficiency of system operation.
Figure 9 shows that the higher electric power consumption of the P2G equipment at night reduces the wind loss caused by the reverse peaking characteristic of wind power, and the total amount of methane generated by the equipment accounts for approximately 6.05% of the total consumption of the system, which, to a certain extent, reduces the dependence of the system on the external gas source, improves the internal energy conversion capability of the system, and further improves the carbon trading revenue of the system due to the carbon reduction characteristic of the system itself.In addition, the amount of the residual heat heating volume of the P2G equipment accounts for approximately 1.57% of the total heat energy demand, indicating that the P2G equipment also shows certain heating potential.Thus, the P2G equipment has manifold Air separation oxygen generation unit operation.

FIGURE 12
Capacity of the oxygen storage tank and oxygen mixing ratio of the oxygen storage tank.
Frontiers in Energy Research frontiersin.orgpotentials, which can effectively improve the operational flexibility and economic revenues of the system.

Effectiveness analysis of the operational strategies
In order to verify the effectiveness of the proposed strategy, four different operational scenarios are set up for comparison and verification, and details of the scenarios are shown in Table 4.
Table 5 shows that the carbon emission and total cost of scenario 2 are higher than those of scenario 1, indicating that the demand response can change the energy use habits of users on the load side, promote the consumption of wind power, and improve the economy and low-carbon of the IES.The carbon transaction cost and carbon emissions of scenarios 3 and 4 are much better than those of scenarios 1 and 2, indicating that the introduction of oxygen-rich combustion capture technology and demand response technology not only reduces the carbon emissions of the system to obtain higher economic benefits but also improves the flexibility of the system and reduces the amount of the system wind abandonment.
Figure 10 shows that the electricity, gas, heat, and cooling loads change in scenarios 2 and 3 before and after introducing the demand response and oxygen-rich combustion capture technology.Under the guidance of time-of-use pricing, the electric and gas loads are transferred from the peak period of energy consumption to the valley period, achieving the effect of load "peak shaving and valley filling."On the other hand, the cooling and heat loads are reasonably adjusted within the comfort level according to the indoor temperature, water temperature changes, and external heat disturbance factors, and certain adjustment incentive subsidies are provided.For electricity and gas loads, after considering the incentive-based IDR, different reductions occur during the later time periods of electricity and gas consumption, respectively.After the introduction of oxygen-rich combustion capture technology and demand response technology, the electricity load is guided by the time-of-use electricity price to transfer to the peak period of the wind power output.The oxygen-rich combustion capture technology regulates the output of gas units in the IES, and the proportion of clean energy increases significantly, indicating that the power supply capacity within the system has increased.In addition, the carbon emissions have been reduced significantly so as to realize the low-carbon and economic operation of the IES.
Combined with Table 6 and Figure 11, it can be seen that in the time period from 5 to 10 o'clock, the power consumption of the air separation oxygen generation equipment in scenario 4 is significantly lower than that in scenario 3, and the oxygen supply of the P2G equipment in scenario 4 accounts for approximately 4.17% of the total oxygen generation of the system, which effectively reduces the oxygen consumption of the gas turbine.Moreover, the gas turbine output of scenario 4 accounts for approximately 75.55% of the total output of the turbine, which is approximately 2.1% higher than that of scenario 3, which shows that the oxygen recovery of the power conversion equipment can effectively reduce the oxygen consumption, improve the net output level of the gas turbine, and reduce the operating cost of the integrated energy system.
In summary, this strategy can improve the output level of the system units, improve the consumption of wind power, tap the operating potential of the P2G equipment and the oxygen equipment, and take into account the low-carbon and economic operation of the system, which verifies the effectiveness of this strategy in the low-carbon and economic operation of the IES.

Capacity configuration analysis of the air separation oxygen generation equipment and oxygen storage tank
Figure 12 shows that the change trend of the total cost of the system first decreases and then increases with the capacity increase of the oxygen storage tank and air separation oxygen generation equipment, and there is a minimum point (30%, 10,000, 454370.14).

FIGURE 13
Analysis of the impact of the carbon trading base price on operating costs and carbon emissions.
Frontiers in Energy Research frontiersin.org The reason is that with the increase in the capacity of the oxygen storage tank and air separation oxygen generation equipment, the maintenance cost of the equipment gradually increases, but the carbon capture volume is also in the upward trend.The oxygen storage tank capacity has little effect on the carbon emission, while air separation capacity has a greater effect on it, and the overall trend of the carbon emission is to decrease with the increase in air separation capacity.In addition, the trend levels are off after 10,000 m 3 because the capacity of the installed 1-MW thermal power unit requires up to 10,000 m 3 /h oxygen.

Analysis of the carbon trading price
With the proposal of the "dual-carbon" goal, low-carbon emission has become a development trend, and carbon trading prices have changed accordingly.Different carbon trading prices are analyzed in this article, and the analysis of the system economy and low-carbon emission under four options is conducted, which proves the rationality of introducing demand response and oxygen-rich combustion units.
Figure 13 shows that the costs of option 1 and option 2 have been continuously increasing, and the capture of oxygen-rich combustion first increases and then decreases after combustion.The reason is that with the increase in carbon trading, oxygen-rich combustion capture units after combustion can do carbon capture and obtain profits.The reason for the high cost of oxygen-rich combustion units at 10-50 CNY/t is that the oxygen generation unit and P2G equipment produce oxygen at low-carbon trading prices but do not consume carbon dioxide.The oxygen-rich combustion capture units start to do the carbon capture when the carbon trading price is over 50 CNY/t, resulting in a decrease in the cost.When the carbon trading price is 120 CNY/t, the system cost of oxygen-rich combustion power plants is lower than that of conventional power plants, while the system cost capture of power plants after combustion are not lower than that of conventional power plants, proving the economic advantages of oxygen-rich combustion power plants.The P2G equipment in option 4 consumes a portion of CO 2 , resulting in a lower carbon trading price than that in option 3.
It can be seen that in most cases, the carbon emissions of the oxygen-rich combustion power plant system are lower than that of the post-combustion capture power plant.The reason for the high carbon trading price of 10-50 CNY/t is that they have not undergone carbon capture, and the oxygen-rich combustion power plant also needs to supply air separation and oxygen production facilities for operation, consuming additional electricity and generating additional CO 2 .When the carbon trading price is 30 CNY/t, the carbon emissions of the oxygenrich combustion power plant system decrease first compared to the post-combustion capture system, and there is more room for the decrease, proving the advantages of oxygen-rich combustion power plants in carbon capture.

Conclusion
This article introduces the oxygen-rich combustion capture unit and demand response into the IES, considering the residual heat recovery and oxygen utilization of the P2G equipment.In addition, an IES operation optimization method considering carbon trading efficiency is proposed.The following conclusions can be drawn through the example analysis: 1.The introduction of oxygen-rich combustion capture technology and demand response can meet the multi-energy requirements within the system, achieve the bi-directional conversion of multiple energies, improve the operational flexibility of the system, and have good carbon reduction effects and economic efficiency.2. The capacity of the air separation oxygen generation equipment and oxygen storage tank has the most cost advantage.However, at this time, the carbon emissions are not optimal and need to be configured according to the actual situation.3.Under the action of the integrated demand response, the IES has achieved the effect of load "peak shaving and valley filling," with a significant increase in the proportion of clean energy and an increase in the output of low-carbon units within the system, improving the low-carbon and economic efficiency of the IES. 4. By comparing the costs and carbon emissions under different carbon trading prices, it has been proven that the scheduling scheme of the oxygen-rich combustion power plant system can achieve the optimal effect in terms of economy and low-carbon performance.

FIGURE 4
FIGURE 4Air liquid energy storage architecture.
FIGURE 5Demand response low-carbon principles.

TABLE 1
Time-of-use electricity price and time-of-use gas price.

TABLE 2
Basic parameters of the integrated energy system.

TABLE 3
Basic parameters of oxygen-rich combustion capture technology.

TABLE 4
Settings of different scenarios.

TABLE 5
Operation strategy in different scenarios.

TABLE 6
Oxygen blending ratio and oxygen consumption of air separation oxygen generation in scenario 4.