Next Article in Journal
Continuous Control Set Predictive Control with Affine Registration Technique for Permanent Magnet Synchronous Motor Drive
Previous Article in Journal
Municipal Solid Waste as a Renewable Energy Source: Advances in Thermochemical Conversion Technologies and Environmental Impacts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Energy System Dispatch Considering Carbon Trading Mechanisms and Refined Demand Response for Electricity, Heat, and Gas

School of Electrical and Electronic Engineering, Changchun University of Technology, Changchun 130012, China
*
Author to whom correspondence should be addressed.
Energies 2024, 17(18), 4705; https://doi.org/10.3390/en17184705
Submission received: 26 August 2024 / Revised: 18 September 2024 / Accepted: 20 September 2024 / Published: 21 September 2024
(This article belongs to the Section C: Energy Economics and Policy)

Abstract

:
To realize a carbon-efficient and economically optimized dispatch of the integrated energy system (IES), this paper introduces a highly efficient dispatch strategy that integrates demand response within a tiered carbon trading mechanism. Firstly, an efficient dispatch model making use of CHP and P2G technologies is developed to strengthen the flexibility of the IES. Secondly, an improved demand response model based on the price elasticity matrix and the capacity for the substitution of energy supply modes is constructed, taking into account three different kinds of loads: heat, gas, and electricity. Subsequently, the implementation of a reward and penalty-based tiered carbon trading mechanism regulates the system’s carbon trading costs and emissions. Ultimately, the goal of the objective function is to minimize the overall costs, encompassing energy purchase, operation and maintenance, carbon trading, and compensation. The original problem is reformulated into a mixed-integer linear programming problem, which is solved using CPLEX. The simulation results from four example scenarios demonstrate that, compared with the conventional carbon trading approach, the aggregate system costs are reduced by 2.44% and carbon emissions are reduced by 3.93% when incorporating the tiered carbon trading mechanism. Subsequent to the adoption of demand response, there is a 2.47% decrease in the total system cost. The proposed scheduling strategy is validated as valuable to ensure the low-carbon and economically efficient functioning of the integrated energy system.

1. Introduction

Amidst rapid economic development, the reliance on fossil fuels has led to significant issues, including the diminishment of non-renewable energy sources and the aggravation of climate change worldwide. Consequently, governments have prioritized the advancement of renewable energy sources as a crucial strategy to lead the transition in energy and develop an environmentally conscious society in addressing the fossil energy issues and environmental damage [1,2]. In recent years, the advancement and implementation of IESs have been pivotal in enhancing the consumption of renewable energy and improving energy utilization [3]. Optimized scheduling is an essential prerequisite for energy production, consumption, and balancing supply and demand in an IES. It is also a core element to ensure coordination, complementarity, and the economical operation of the diverse energy subsystems within the system [4,5].
Demand response (DR) can facilitate the synergistic correlation between supply and demand in an integrated energy system, mitigate the disparities between peak and off-peak loads, and enhance the economic efficiency of IES operations [6,7]. Common load demand response can be categorized into substitutable loads, transferable loads, and curtailable loads based on customer-side demand characteristics. Furthermore, based on the response type, it can be grouped into price-based and incentive-based demand responses [8,9]. In an IES, electrical, thermal, gas, and cooling loads exhibit distinct characteristics. Even within the same load category, response behaviors may differ significantly when exposed to varying incentives. Thus, the integration of DR can offer greater flexibility in IES scheduling [10].
In [11], the authors proposed the incorporation of incentive-based DR in microgrid energy management to maximize benefits for participating users. In [12,13], demand response according to time-of-day pricing was incorporated into IES scheduling to facilitate the rational coordination of multiple energy forms to meet consumption needs. In [14], the incorporation of time-of-day tariff-based demand response in an IES with CHP units can alleviate pressure on the power supply and enhance system stability. However, [12,14] only analyzed the demand response of electric loads, without considering other forms of energy loads. In [15,16], a demand response mechanism considering both electric and thermal loads is introduced to optimize integrated energy system dispatch, satisfying users’ energy requirements and enhancing renewable energy utilization. However, the demand response modeling method needs further optimization. In [17], demand response modeling was optimized, but only electric and thermal loads were considered. In [18], a flexible renewable resource model was developed and utilized in the efficient allocation of an IES. However, there are issues with the insufficient refinement of flexible resource classification in demand response and the oversimplified substitution relationships between different energy loads. In [19], it is pointed out that introducing a carbon trading mechanism can enable the low-carbon and economic operation of an IES. In [20], CHP and P2G were incorporated into the IES structure within a carbon trading mechanism, enhancing system’s operational flexibility. While CHP technology markedly enhances the overall efficiency of energy utilization through the simultaneous generation of electricity and heat, P2G technology offers a highly effective solution for conversion and energy storage, proving particularly advantageous for the integration of renewable energy sources. In comparison to other emerging technologies, such as hydrogen, P2G not only offers flexible energy storage capabilities but also assumes a critical role in multi-energy complementary systems. In [21], it is demonstrated that a dispatch scheme combining carbon trading mechanisms and demand response can simultaneously reduce greenhouse gas emissions and overall operational expenses. In [22,23], the traditional carbon trading mechanism was introduced in IES scheduling; however, this mechanism is overly simplistic for modeling carbon trading and costs. Therefore, [24] introduced a tiered carbon trading system in the scheduling plan to stratify carbon trading costs and further reduce emissions while ensuring economic operation. In [25], it was highlighted that the reward and penalty ladder carbon trading mechanism can more effectively stimulate user participation to adopt emission reduction measures. In reference [26], a dispatch model for IESs addressing electricity, heat, and gas under a tiered carbon trading approach was proposed. However, the modeling of carbon trading costs and equipment emissions in this study is insufficiently comprehensive. In [27], the modeling of carbon transaction costs and equipment emissions was optimized and the tiered carbon trading approach with incentives and penalties was introduced. However, neither [26] nor [27] considered demand response, resulting in an inability to regulate the user side. In [28], the DR of electricity, heat, and gas loads within the incentive–penalty ladder carbon trading mechanism was considered, but only load reduction and leveling were addressed, without considering the mutual substitution between loads.
To address the issues highlighted above, this paper suggests an integrated energy system dispatch model that incorporates a carbon trading mechanism and electricity–heat–gas DR, making the following key contributions:
  • The integration of CHP and P2G technologies into the IES scheduling model significantly enhances the operational flexibility of the system. CHP (combined heat and power) and P2G (power-to-gas) technologies improve the efficiency of converting and distributing electricity, heat, and gas, thereby optimizing the overall operational efficiency of the system.
  • The incorporation of a carbon trading mechanism into IES scheduling thoroughly accounts for the impact of each kind of equipment on carbon emissions. A reward-and-penalty-based tiered carbon trading model is established to effectively mitigate carbon emissions and trading costs. This mechanism incentivizes users to implement additional emission reduction measures, thereby achieving low-carbon operation through tiered rewards and penalties.
  • A refined DR model based on the price elasticity matrix is developed, further considering the interplay between electricity, heat, and gas load substitutions. By accounting for each load’s reduction and leveling capabilities, the model maximizes the regulatory potential of flexible resources. This model can dynamically adjust load allocation in response to price signals, thereby enhancing the system’s responsiveness and economic performance.
The subsequent sections of this paper are arranged as follows: Section 2 introduces the structure of the integrated energy system and the demand response model. Section 3 introduces a reward-and-penalty-based tiered carbon trading model. Section 4 presents the scheduling model and solution method for the IES. Section 5 provides a case study with four different scenarios to illustrate the effectiveness of the suggested scheduling strategy. Section 6 concludes the paper.

2. IES Structure

The proposed IES optimization framework is fully aligned with both international and national climate objectives, including the Paris Agreement and China’s dual goals of “carbon peaking” and “carbon neutrality”. Through the integration of carbon trading mechanisms and demand response strategies, the framework effectively contributes to carbon emission reductions across multiple energy sectors, including electricity, heat, and gas. The incorporation of technologies such as CHP and P2G enables the more efficient utilization of energy resources, facilitating the transition toward cleaner, more sustainable energy systems. Additionally, the framework’s capacity to optimize energy dispatch while factoring in carbon costs provides a practical and effective solution for managing carbon emissions, consistent with the national target of peaking emissions by 2030 and attaining carbon neutrality by 2060. This framework not only assists policymakers in designing low-carbon energy systems but also equips industrial sectors with strategic tools to achieve their emission reduction targets, underscoring its practical importance in advancing climate action. The structural framework of the proposed IES is illustrated in Figure 1. The constructed IES structure includes the superior grid, superior gas grid, wind turbines (WTs), and photovoltaic panels (PVs) on the energy supply side, and the combined heat and power (CHP), electric boiler (EB), gas boiler (GB), and power-to-gas (P2G) on the energy conversion side. Specifically, the CHP comprises a gas turbine (GT), an organic Rankine cycle (ORC), and a waste heat boiler (WHB).

2.1. Demand Response Model

The DR model proposed in this study simultaneously addresses the DR of three types of loads: electricity, heat, and gas, thereby optimizing the utilization of flexible resources. These loads are categorized into curtailable loads, shiftable loads, and replaceable loads. Based on load response properties, DR is classified into price-based DR and substitution-based DR. Price-driven demand response loads are further categorized into curtailable loads (CLs) and shiftable loads (SLs).

2.1.1. Price-Based Demand Response

CLs determine whether to scale down their load based on the price at certain points prior to and following the DR. The price elasticity of the demand matrix for curtailable loads is defined as E C L ( i , j ) :
E C L ( i , j ) = e 11 C L 0 0 0 e 22 C L 0 0 0 e i j C L
e i j C L = Δ P i C L P i C L .0 Δ ϕ j ϕ j 0 1
e i j C L is the element within the elasticity matrix E C L ( i , j ) , representing the elasticity coefficient available of the CL at moment i in relation to the price at moment j. P i C L .0 denotes the initial CL and Δ P i C L denotes the variation in CL following DR. ϕ j 0 is the original price, and Δ ϕ j is the variation in price after DR. ρ j 0 represents the initial tariff, while ρ j represents the TOU tariff. Therefore, the change in CL after DR at moment i, denoted as Δ P C L , i , can be described as follows:
Δ P C L , i = P C L , i 0 [ j = 1 24 E C L ( i , j ) ρ j ρ j 0 ρ j 0 ]
Shiftable loads (SLs) refer to loads where users respond to prices according to their own needs, enabling flexible adjustments in working hours. The time-of-use (TOU) pricing is used as a signal to guide users in shifting their load from high-demand periods to low-demand periods. The price elasticity matrix for shiftable loads is denoted as E S L ( i , j ) :
E S L ( i , j ) = e 11 S L e 12 S L e 1 j S L e 21 S L e 22 S L e 2 j S L e i 1 S L e i 2 S L e i j S L
e i j S L = Δ P i S L P i S L .0 Δ ϕ j ϕ j 0 1
e i j S L is the element of the elasticity matrix E S L ( i , j ) , representing the elasticity coefficient available of the SL at moment i in relation to the price at moment j. P i S L , 0 represents the initial SL, and Δ P i S L denotes the variation in SL following DR. Therefore, the change in SL after DR, denoted as Δ P S L , i , can be described as follows:
Δ P S L , i = P S L , i 0 [ j = 1 24 E S L ( i , j ) ρ j ρ j 0 ρ j 0 ]

2.1.2. Replacement Demand Response

The alternative demand response strategy seeks to optimize overall energy efficiency by coordinating the substitution of electric, heating, and gas loads. The core concept involves leveraging the price elasticity of each load to flexibly adjust their timing and usage, thereby reducing total operating costs and enhancing system stability. The following description specifically addresses the mutual substitution of electric, heating, and gas loads.
(1)
Mutual Replacement of Electric and Heating Loads
Mutual replacement of electric and heating loads occurs when users, guided by the prices of electricity and heat, adjust between electric and heat loads. For example, users can increase the electric load and replace part of the heat load with electric equipment during periods of low electricity prices. Conversely, when heat prices are low, users can cut down on electrical load and boost the use of heating load. The replaceable load (RL) can be described as follows:
Δ P i R L , e = ξ e , h Δ P i R L , h
ξ e , h = η e λ e η h λ h
where Δ P i R L , e and Δ P i R L , h represent the amounts of electric RL and heating RL, respectively. Here, ξ e , h represents the electric heat substitution factor, η e and η h denote the substitution efficiencies for electric and heating energy, and λ e and λ h indicate the thermal values per unit of electric and heating energy, respectively.
(2)
Mutual Replacement of Electric and Gas loads
Mutual replacement of electric and gas loads occurs when users, guided by electricity and gas prices, adjust between electric and gas loads. For example, when electricity prices are low, users can increase electric loads and use electrical equipment to replace some of the gas loads. Conversely, when gas prices are low, users can reduce the electric load and increase the use of gas loads. The RL can be described as follows:
Δ P i R L , e = ξ e , g Δ P i R L , g
ξ e , g = η e λ e η g λ g
where Δ P i R L , e and Δ P i R L , g represent the amounts of electric RL and gas RL, respectively. ξ e , g represents the electric gas substitution factor. η e and η g denote the electric and gas substitution efficiencies, while λ e and λ g are the thermal values of electric and gas energy, respectively.

3. Carbon Trading Mechanism Model

The carbon emissions trading mechanism is a policy instrument designed to control and reduce greenhouse gas emissions through market-based tools. The core of the mechanism involves the establishment of legal carbon emission credits and their trade in a specialized carbon market. In practice, the administration initially distributes a specific amount of emission allowances to emission sources, such as energy suppliers, free of charge. If the actual emissions of carbon from the energy suppliers are lower than the issued credits, they can sell the excess trading credits at the prevailing carbon market price and obtain financial gains. This mechanism regulates the cost of carbon emissions through market supply and demand, incentivizes enterprises to reduce emissions, and ultimately drives down overall carbon emissions.

3.1. Initial Carbon Emission Assignment

The allocation of initial carbon allocation in an IES comprises four components: the electricity acquired from the superior grid, GT, GB, and gas load. Consequently, the carbon credits of the IES can be expressed as follows:
E I = E g r i d + E G T + E G B + E g l o a d
E g r i d = ε e i = 1 T P i g r i d
E G T = ε h i = 1 T ( ε e , h P i G T , e + P i G T , h )
E G B = ε h i = 1 T P i G B
E g l o a d = ε g l o a d i = 1 T P i g l o a d
Here, E I represents the initial carbon emission allowance; E g r i d , E G B , E G T , and E g l o a d denote the carbon emission allowances for power purchased from the superior grid, GB, GT, and gas loads, respectively. P i g r i d is the power purchased from the higher level, and ε e is the unremunerated carbon emission allowance per unit of electricity. P i G T , e and P i G T , h are the electric power and heat power of the GT, respectively; ε h is the unremunerated carbon emission entitlement for each unit of heat; ε e , h is the conversion factor from electricity generation to heat generation; P i G B is the GB heat power; P i g l o a d is the power of the gas load; and ε g l o a d is the unremunerated carbon emission entitlement of the gas load.

3.2. Actual Carbon Emissions

The consumption of electricity and gas serves as a pivotal determinant of the carbon emissions across the entire energy framework. Moreover, in power-to-gas (P2G) technology, carbon dioxide is utilized as a fundamental raw material during converting electricity into natural gas. The carbon dioxide used in the P2G process should be encompassed within the carbon trading system. Additionally, carbon emissions resulting from gas loads are also considered. Consequently, the effective carbon emissions of the IES can be articulated as described below:
E A = E g r i d , a + E g l o a d , a + ( G G T + G G B G P 2 G ) ε g
G G T = i = 1 T P i e , G T η G T , e λ g a s + i = 1 T P i h , G T η G T , h λ g a s
G G B = i = 1 T P i G B η G B λ g a s
G P 2 G = i = 1 T P i P 2 G η P 2 G λ g a s
E g r i d , a = ε p i = 1 T P i g r i d
E g l o a d , a = ε q i = 1 T P i g l o o d
where E A denotes the actual carbon emissions; ε g , ε p , and ε q represent the carbon emission coefficients for gas consumption, purchased electricity, and gas load, respectively; G G T indicates the gas consumption of the GT; η G T , e and η G T , h denote the efficiencies of the GT in electricity and heat production; λ g a s refers to the calorific value of natural gas; G G B signifies the gas consumption of the GB; η G B denotes the heat production efficiency of the GB; G P 2 G represents the amount of absorbed gas in the P2G process; P i P 2 G is the electricity consumed by the P2G process; η P 2 G refers to the electricity-to-gas efficiency of the P2G process; E g r i d , a denotes the actual carbon emissions associated with purchased electricity; and E g l o a d , a represents the actual carbon emissions from the gas load.

3.3. Laddered Carbon Trading Mechanism

The study establishes a tiered carbon trading system with rewards and penalties grounded on a unified ladder carbon trading framework, predicated on the correlation between E A and E I . To mitigate carbon dioxide emissions, the difference between E A and E I is bifurcated into several intervals. When negative, it indicates that E A is below E I , allowing the IES to sell the surplus quota for a profit; when positive, it signifies that E A is above E I , necessitating that the IES purchase additional quotas to meet its demand. The greater the excess, the higher the carbon trading cost. In summary, the model of the tiered carbon trading system with rewards and penalties is articulated as follows:
C C O 2 = γ ( 2 + 3 ω ) υ + γ ( 1 + 3 ω ) ( E A E I + 2 υ ) , E A E I 2 υ γ ( 1 + ω ) υ + γ ( 1 + 2 ω ) ( E A E I + υ ) , 2 υ < E A E I υ γ ( 1 + ω ) E A E I , υ < E A E I 0 γ E A E I , 0 < E A E I υ γ ( 1 + ρ ) ( E A E I υ ) + γ υ , υ < E A E I 2 υ γ ( 1 + 2 ρ ) ( E A E I 2 υ ) + ( 2 + ρ ) γ υ , 2 υ < E A E I 3 υ γ ( 1 + 3 ρ ) ( E A E I 3 υ ) + ( 3 + 3 ρ ) γ υ , 3 υ < E A E I 4 υ γ ( 1 + 4 ρ ) ( E A E I 4 υ ) + ( 4 + 6 ρ ) γ υ , 4 υ < E A E I
C C O 2 denotes the carbon trading cost; γ signifies the base rate for carbon trading; ρ represents the rate of price advancement; υ indicates the length of the carbon emission interval; and ω is the incentive factor.

4. IES Scheduling Model

4.1. Objective Function

In this paper, the purpose of lowering the total operating cost is achieved through the implementation of DR and a stepped carbon trading mechanism. The total operating cost comprises energy purchase costs, operation and maintenance costs, carbon trading costs, and compensation costs.
min f = f b u y + f o p + f C O 2 + f D R
(1)
Cost of purchasing energy
f b u y = i = 1 T p e , i P e , b u y ( i ) + i = 1 T p g , i P g , b u y ( i )
where P e , b u y ( i ) and P g , b u y ( i ) represent the power of electricity and gas purchased; p e , i and p g , i denote the electricity and gas price, respectively.
(2)
Operation and maintenance costs
f o p = i = 1 T t = 1 9 ω t P t , i
where t takes values from 1, 2, …, 8, and 9 to represent the GT, WHB, ORC, EB, GB, ES, GS, HS, and the P2G process, respectively. Here, ω t is the O&M cost coefficient of equipment t, and P t , i is the output power of equipment t.
(3)
Carbon trade costs
Carbon trade costs are detailed in Equation (22).
(4)
Compensation costs
f D R = i = 1 T ζ R Δ P i R L , e + Δ P i R L , h + Δ P i R L , g + ζ S Δ P SL , i e + Δ P SL , i g + Δ P SL , i h
where ζ R and ζ S represent the compensation cost coefficients for shiftable and substitutable loads, respectively.

4.2. Constraint

(1)
Electric power balance constraints
P i W T + P i P V + P i e , b u y + P i e , C H P + P i e , c h P i e , P 2 G P i e , E B P i e , d i s = P i e , l o a d + Δ P C L , i e + Δ P S L , i e + Δ P i R L , e
P i W T and P i P V represent the power generated by the WT and PV, respectively. P i e , C H P is the power produced by the CHP unit; P i e , c h and P i e , d i s are the ES’s charge and discharge electric power, respectively; P i e , P 2 G is the power consumed by the P2G process; P i e , E B is the power consumed by the EB; P i e , l o a d is the amount of electric loads before Demand response (DR); and Δ P CL , i e , Δ P S L , i e , and Δ P i R L , e are the electric load amounts of the CL, SL, and RL after DR, respectively.
(2)
Heat power balance constraints
P i h , C H P + P i h , E B + P i h , G B + P i h , c h P i h , d i s = P i h , l o a d + Δ P C L , i h + Δ P S L , i h + Δ P i R L , h
P i h , C H P represents the heat power of the CHP unit; P i h , E B and P i h , G B are the heat production powers of the EB and GB, respectively; P i h , c h and P i h , d i s are the HS’s charge and discharge heat power, respectively; P i h , l o a d is the amount of heat load before DR; and Δ P C L , i h , Δ P S L , i h , and Δ P i R L , h are the heat load amounts of the CL, SL, and RL after DR, respectively.
(3)
Gas power balance constraints
P i g , b u y + P i g , c h + P i g , P 2 G P i g , C H P P i g , G B P i g , d i s = P i g , l o a d + Δ P C L , i g + Δ P S L , i g + Δ P i R L , g
P i g , b u y represents the purchased power; P i g , c h and P i g , d i s are the charging and discharging powers of the GS; P i g , C H P is the gas consumption power of the CHP unit; P i g , P 2 G is the gas production power of the P2G system; P i g , G B is the gas consumption power of the GB; P i g , l o a d is the gas loading capacity before DR; and Δ P C L , i g , Δ P S L , i g , and Δ P i R L , g are the gas load amounts of the CL, SL, and RL after DR, respectively.
(4)
CHP constraints
P i e , C H P = P i e , G T + P i e , O R C
P i h , C H P = P i h , G T β i τ W H B
P i e , O R C = P i h , G T α i δ O R C
P i e , G T = P i g 2 e , C H P τ G T e
P i h , G T = P i g 2 h , C H P τ G T h
P i g , C H P = P i g 2 e , C H P + P i g 2 h , C H P
α i + β i = 1
0 α i , β i 1
The constraint on the electric power production of the CHP unit, denoted as P i e , C H P , encompasses two components: the power generated by the GT and the power produced by the ORC. The constraint on the heat production of the CHP unit, represented as P i h , C H P , includes the heat generated by the WHB, which is accounted for as part of the CHP’s total heat production. The electric power production of the ORC at time i is indicated by P i e , O R C . The heat power production of the GT at time i is denoted as P i h , G T . The proportion of the residual heat from the GT allocated to the ORC is represented by α i . The electric power generation efficiency of the ORC is denoted by δ O R C . The heat conversion efficiency of the WHB is indicated by τ W H B . The gas consumption power for gas-to-electricity and gas-to-heat conversions are denoted by P i g 2 e , C H P and P i g 2 e , C H P , respectively. The efficiencies of gas-to-electricity and gas-to-heat conversions are represented by τ G T e and τ G T h , respectively. The calorific value of natural gas, denoted as λ g a s , is 9.88 kWh/m3.
(5)
P2G constraints
0 P i P 2 G P P 2 G , max
P i P 2 G P i 1 P 2 G r P 2 G
The parameter P P 2 G , max represents the upper output limit of the P2G process; the parameter r P 2 G denotes the maximum ramp rate of the P2G process.
(6)
EB constraints
0 P i E B P E B , max
P i E B P i 1 E B r E B
The parameter P E B , max represents the upper output limit of the EB; the parameter r E B denotes the maximum ramp rate of the EB.
(7)
GB constraints
0 P i G B P G B , max
P i G B P i 1 G B r G B
The parameter P G B , max represents the upper output limit of the GB; the parameter r G B denotes the maximum ramp rate of the GB.
(8)
GT constraints
0 P i G T P G T , max
P i G T P i 1 G T r G T
The parameter P G T , max represents the upper output limit of the GT; The parameter r G T denotes the maximum ramp rate of the GT.
(9)
ORC constraints
0 P i O R C P O R C , max
P i O R C P i 1 O R C r O R C
The parameter P O R C , max represents the upper output limit of the ORC; the parameter r O R C denotes the maximum ramp rate of the ORC.
(10)
Energy storage constraints
Three different kinds of energy storage devices are incorporated into the IES to further increase its operating flexibility: electrical, thermal, and gas. These storage systems are modeled uniformly due to the inherent similarities in their respective energy storage mechanisms.
0 P i n , d i s P i n , d i s , max τ i d i s
0 P i n , c h P i n , c h , max τ i c h
τ i c h + τ i d i s 1
S i + 1 n = S i n + η n , c h P i n , c h η n , d i s P i n , d i s
S n , min S i n S n , max
S 0 n = S T n
Here, the letter n stands for electrical, thermal, and gas energy storage devices, respectively; the maximum power for charging and discharging is represented by P i n , d i s , max and P i n , c h , max ; the binary variables τ i d i s and τ i c h represent the energy storage device’s charging and discharging states; the energy stored in the energy storage device is denoted by the symbol S i n ; η n , c h and η n , d i s denote the charging and discharging efficiencies; S n , min and S n , max indicate the energy storage device’s top and lower capacity limitations.
The comprehensive structure of the constraints is depicted in Figure 2.

4.3. Solution Method

The IES scheduling model, which is stated as a mixed-integer linear programming (MILP) problem, incorporates DR and a tiered carbon trading system. It is discussed in this study. Therefore, MATLAB and CPLEX are utilized to solve the model. First, the established IES model is supplemented with data on wind power and photovoltaic inputs, as well as heat, gas, and electricity demands. Subsequently, the parameters of each unit are incorporated. Finally, the model aims to minimize the overall cost, which derives the best scheduling strategy for each unit given the restrictions of power balance and unit-specific constraints. The program flow chart is shown in Figure 3.

5. Case Analysis

5.1. Parameter Setting

An industrial park in the north that operates for 24 h a day, seven days a week, with an hourly unit of operation is the focus of this study. Table 1 presents the parameters of the energy conversion equipment within the IES, whereas Table 2 presents the specs of the energy storage equipment. Parameters pertaining to the carbon trading mechanism are also included in Table 3. Figure 4 illustrates the electrical load, thermal load, gas load, and the power generated by wind and photovoltaic sources. The TOU pricing for electricity procured from the main grid and the gas obtained from the primary gas grid is depicted in Figure 5a and Figure 5b, respectively. Prior to incorporating DR, the IES sells electricity at 0.78 CNY/kWh, gas at 0.36 CNY/kWh, and heat at 0.67 CNY/kWh. Fixed loads constitute 50% of the total load, with the SL, RL, and CL comprising 35%, 10%, and 5%, respectively. Four different scenarios have been created for comparison analysis in order to look at how the IES developed in this study will be affected by the inclusion of DR and the carbon trading system.
Scenario 1 utilizes the traditional carbon trading mechanism without DR consideration.
Scenario 2 utilizes the traditional carbon trading mechanism with DR consideration.
Scenario 3 employs the laddered carbon trading mechanism without DR consideration.
Scenario 4 employs the laddered carbon trading mechanism with DR consideration.

5.2. Scenario Comparison

Figure 6 illustrates the initial and TOU prices of energy purchased by users from the IES before and after considering DR. As depicted in Figure 6a, the peak hours for gas prices are from 6:00 to 12:00 and 18:00 to 21:00, while the valley hours are from 13:00 to 17:00 and 22:00 to 24:00. As shown in Figure 6b, the peak hours for heat prices occur between 4:00 and 6:00, 12:00 and 14:00, and 19:00 and 24:00, with valley hours occur from 1:00 to 3:00. As illustrated in Figure 6c, the peak hours for electricity prices are from 8:00 to 11:00 and 20:00 to 22:00, with valley hours span from 1:00 to 7:00.
The gas load, heat load, and electricity load, respectively, are shown in Figure 7a–c), both before and after taking DR into account, as well as the variations in each kind of load. Some of the gas load is moved to low gas price hours during high gas price hours (6:00–12:00 and 18:00–21:00). During periods of high gas prices and low electricity prices (6:00–7:00, 12:00, and 19:00), some gas loads are converted to electricity loads. Conversely, during periods of low gas prices and high electricity prices (22:00), some electricity loads are converted to gas loads. During high heat price periods (4:00–6:00, 12:00–14:00, and 19:00–24:00), a portion of the heat load is shifted to low heat price periods. During periods of high heat prices and low electricity prices (4:00–5:00 and 13:00–14:00), some heat loads are converted to electricity loads. During periods of low heat prices and high electricity prices (8:00–11:00), a portion of the electric load is converted to heat load. During high tariff hours (8:00–11:00 and 20:00–22:00), a part of the electricity demand is moved to hours with reduced tariffs, resulting in a curtailment of electric load.
The overall balances of thermal, electrical, and gas power in Scenario 4 are shown in Figure 8a, Figure 8b and Figure 8c, respectively. As shown in Figure 8a, The thermal energy requirement is mostly satisfied by the CHP and GB during the hours of 1:00–3:00 and 7:00–18:00, with the HS storing any spare thermal energy. During the periods of 4:00–6:00 and 19:00–24:00, the thermal energy demand is at its highest. Therefore, to satisfy the thermal energy requirement, EB heat generation and HS heat release are also required in addition to the CHP and GB.
As depicted in Figure 8b, during the periods of 1:00–3:00, 7:00–8:00, and 12:00–17:00, the power demand is lower, and it can be met by the CHP, WT, and PV alone, with excess power stored by the P2G system and the ES. During the peak periods of 4:00–6:00, 9:00–11:00, and 19:00–24:00, in addition to the CHP, PV, and WT, the power consumption of the EB and peak demand can be met by purchasing power and discharging the ES. At 22:00, the P2G system needs to store energy for the conversion from electricity to gas. As shown in Figure 8c, the gas demand, along with the GB and CHP consumption, is primarily met by purchased gas, while the P2G process and the GS help reduce the amount of purchased gas.
The power balance of thermal, electrical, and gaseous energy in Scenario 3 are shown in Figure 9a, Figure 9b and Figure 9c, respectively. Each load is left as a raw load without load peak changing since DR is not taken into account. As a consequence, every device’s output is impacted, raising the overall cost.
The simulation results for the four scenarios are presented in Table 4, with the total cost for Scenario 1 reaching CNY 50,198.48. This is attributed to the lack of consideration for DR and the ladder carbon trading mechanism, resulting in the inflexible scheduling of the IES. In comparison to Scenario 1, Scenario 2 incorporates DR, which alters users’ energy consumption habits and shifts part of the load to periods with lower energy prices. Additionally, electric, thermal, and gas loads can be substituted with lower-cost alternatives, enhancing the flexibility of IES scheduling, reducing operational pressure on equipment, and lowering O&M costs by 8% and total costs by 2.2%. When comparing Scenario 3 to Scenario 1, the stepped carbon trading mechanism in Scenario 3 reduces carbon emissions by 3.94%, decreases carbon trading costs by 18%, and lowers total costs by 2.44%. When DR and the stepped carbon trading system are coupled, Scenario 4’s total costs and carbon emissions are reduced by 2.24% and 3.42%, respectively, compared to Scenario 3’s.
In summary, this paper’s suggested stepped carbon trading method successfully lowers emissions of carbon and carbon trading expenses. Additionally, the introduction of DR rationalizes the scheduling of each device within the IES, thereby significantly reducing its operation and maintenance costs. Thus, the scheduling technique described in this study effectively meets low-carbon and economic goals.

6. Discussion

The discussion is divided into two key sections: addressing the uncertainty of user behavior and examining the limitations of the current model.
A significant challenge for the practical implementation of demand response models lies in the inherent unpredictability of user behavior. Users’ energy consumption patterns are influenced by a multitude of factors, such as price volatility, climate variability, and individual preferences. These behavioral fluctuations can result in suboptimal demand response outcomes, thereby diminishing dispatch efficiency and compromising the carbon reduction goals of the integrated energy system. To address this uncertainty, future research could explore the integration of probabilistic models or stochastic optimization techniques to simulate varying user demand scenarios, thereby enhancing the model’s robustness. Moreover, leveraging machine learning algorithms and advanced load forecasting models could offer deeper insights into the complex patterns of user behavior and significantly enhance forecasting accuracy. Finally, the development of more flexible and targeted demand response incentives may help mitigate the scheduling challenges posed by unpredictable user behavior to a certain extent.
While the proposed IES optimization framework offers promising results, several limitations should be acknowledged. First, the use of MILP may introduce significant computational complexity, especially when scaling up to larger systems or incorporating more detailed constraints. This could potentially limit the real-time adaptability of the model in dynamic or fast-changing environments. Furthermore, the real-world application of the proposed model requires comprehensive and accurate data on energy demand, supply, and carbon pricing, which may not always be readily available or reliable. Data uncertainty or unavailability could affect the robustness and accuracy of the optimization results. Addressing these limitations in future research, such as through the use of advanced algorithms or machine learning techniques for load forecasting, could enhance both the computational efficiency and the practical applicability of the system.

7. Conclusions

To enhance the economic efficiency of IES operation and optimize the overall carbon emissions of the system, this paper proposes a scheduling strategy that integrates DR and a stepped carbon trading mechanism. The precise findings drawn from the simulation results mentioned above are as follows:
  • The mutual substitution between the loads of electricity, heat, and gas as well as the shiftable and curtailable loads of each are covered by the suggested DR approach. This method accomplishes peak shaving and valley filling goals while lowering the operating strain on every piece of IES equipment.
  • In comparison to conventional carbon trading schemes, the ladder mechanism decreases carbon emissions more efficiently and makes low-carbon system operation easier.
  • By integrating DR into the ladder trading of the carbon mechanism, the IES’s operating and maintenance expenses as well as energy purchase prices may be decreased, all while lowering the system’s carbon emissions. The IES performs more economically and in terms of carbon emissions thanks to this modification.
Future research will explore the integration of carbon capture technologies to further optimize the IES. Additionally, multi-objective optimization algorithms will be employed to balance and optimize multiple objectives. Furthermore, the application of machine learning techniques will be explored to improve the load forecasting accuracy in the DR model.

Author Contributions

Conceptualization, L.G., S.Y., J.G. and N.C.; methodology, L.G., S.Y. and N.C.; software, S.Y.; validation, S.Y. and J.G.; data curation, S.Y. and J.G.; writing—original draft preparation, L.G., S.Y. and N.C.; writing—review and editing, L.G., S.Y. and N.C.; visualization, S.Y.; supervision, L.G. and N.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Science and Technology Research Program of Jilin Provincial Department of Education, grant number JJKH20240854KJ and the APC was funded by Changchun University of Technology.

Data Availability Statement

Data are contained within paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. You, C.; Kim, J. Optimal Design and Global Sensitivity Analysis of a 100% Renewable Energy Sources Based Smart Energy Network for Electrified and Hydrogen Cities. Energy Convers. Manag. 2020, 223, 113252. [Google Scholar] [CrossRef]
  2. Wang, S.; Meng, Z.; Yuan, S. IEC 61970 Standard Based Common Information Model Extension of Electricity-Gas-Heat Integrated Energy System. Int. J. Electr. Power Energy Syst. 2020, 118, 10584. [Google Scholar] [CrossRef]
  3. Li, J.; Huang, Y.; Zhu, M. Gradient Descent Iterative Method for Energy Flow of Integrated Energy System Considering Multiple Modes of Compressors. Energy Convers. Manag. 2020, 207, 112534. [Google Scholar] [CrossRef]
  4. Wang, Z.; Hu, J.; Liu, B. Stochastic Optimal Dispatching Strategy of Electricity-Hydrogen-Gas-Heat Integrated Energy System Based on Improved Spectral Clustering Method. Int. J. Electr. Power Energy Syst. 2021, 126, 106495. [Google Scholar] [CrossRef]
  5. Alabi, T.M.; Aghimien, E.I.; Agbajor, F.D.; Yang, Z.; Lu, L.; Adeoye, A.R.; Gopaluni, B. A Review on the Integrated Optimization Techniques and Machine Learning Approaches for Modeling, Prediction, and Decision Making on Integrated Energy Systems. Renew. Energy 2022, 194, 822–849. [Google Scholar] [CrossRef]
  6. Li, P.; Wang, Z.; Wang, J.; Yang, W.; Guo, T.; Yin, Y. Two-Stage Optimal Operation of Integrated Energy System Considering Multiple Uncertainties and Integrated Demand Response. Energy 2021, 225, 120256. [Google Scholar] [CrossRef]
  7. Bahrami, S.; Sheikhi, A. From Demand Response in Smart Grid Toward Integrated Demand Response in Smart Energy Hub. IEEE Trans. Smart Grid 2016, 7, 650–658. [Google Scholar] [CrossRef]
  8. Toolabi Moghadam, A.; Soheyli, F.; Sanei, S.; Akbari, E.; Khorramdel, H.; Ghadamyari, M. Bi-level Optimization of the Integrated Energy Systems in the Deregulated Energy Markets Considering the Prediction of Uncertain Parameters and Price-based Demand Response Program. Energy Sci. Eng. 2022, 10, 2772–2793. [Google Scholar] [CrossRef]
  9. Zheng, S.; Sun, Y.; Qi, B.; Li, B. Incentive-Based Integrated Demand Response Considering S&C Effect in Demand Side with Incomplete Information. IEEE Trans. Smart Grid 2022, 13, 4465–4482. [Google Scholar]
  10. Wang, J.; Zhong, H.; Ma, Z.; Xia, Q.; Kang, C. Review and Prospect of Integrated Demand Response in the Multi-Energy System. Appl. Energy 2017, 202, 772–782. [Google Scholar] [CrossRef]
  11. Harsh, P.; Das, D. Optimal Coordination Strategy of Demand Response and Electric Vehicle Aggregators for the Energy Management of Reconfigured Grid-Connected Microgrid. Renew. Sustain. Energy Rev. 2022, 160, 112251. [Google Scholar] [CrossRef]
  12. Luo, Y.; Zhang, X.; Yang, D.; Sun, Q.; Zhang, H. Optimal Operation and Cost–Benefit Allocation for Multi-participant Cooperation of Integrated Energy System. IET Gener. Transm. Distrib. 2019, 13, 5239–5247. [Google Scholar] [CrossRef]
  13. Kiptoo, M.K.; Lotfy, M.E.; Adewuyi, O.B.; Conteh, A.; Howlader, A.M.; Senjyu, T. Integrated Approach for Optimal Techno-Economic Planning for High Renewable Energy-Based Isolated Microgrid Considering Cost of Energy Storage and Demand Response Strategies. Energy Convers. Manag. 2020, 215, 112917. [Google Scholar] [CrossRef]
  14. Majidi, M.; Mohammadi-Ivatloo, B.; Anvari-Moghaddam, A. Optimal Robust Operation of Combined Heat and Power Systems with Demand Response Programs. Appl. Therm. Eng. 2019, 149, 1359–1369. [Google Scholar] [CrossRef]
  15. Shao, C.; Ding, Y.; Wang, J.; Song, Y. Modeling and Integration of Flexible Demand in Heat and Electricity Integrated Energy System. IEEE Trans. Sustain. Energy 2018, 9, 361–370. [Google Scholar] [CrossRef]
  16. Yang, Y.; Yan, G.; Mu, G. Bi-Level Decentralized Control of Electric Heating Loads Considering Wind Power Accommodation in Real-Time Electricity Market. Int. J. Electr. Power Energy Syst. 2022, 135, 107536. [Google Scholar] [CrossRef]
  17. Luo, Y.; Hao, H.; Yang, D.; Zhou, B. Multi-Objective Optimization of Integrated Energy Systems Considering Ladder-Type Carbon Emission Trading and Refined Load Demand Response. J. Mod. Power Syst. Clean Energy 2024, 12, 828–839. [Google Scholar] [CrossRef]
  18. Wang, W.; Huang, S.; Zhang, G.; Liu, J.; Chen, Z. Optimal Operation of an Integrated Electricity-Heat Energy System Considering Flexible Resources Dispatch for Renewable Integration. J. Mod. Power Syst. Clean Energy 2021, 9, 699–710. [Google Scholar] [CrossRef]
  19. Sun, P.; Hao, X.; Wang, J.; Shen, D.; Tian, L. Low-carbon Economic Operation for Integrated Energy System Considering Carbon Trading Mechanism. Energy Sci. Eng. 2021, 9, 2064–2078. [Google Scholar] [CrossRef]
  20. Li, Y.; Bu, F.; Gao, J.; Li, G. Optimal Dispatch of Low-Carbon Integrated Energy System Considering Nuclear Heating and Carbon Trading. J. Clean. Prod. 2022, 378, 134540. [Google Scholar] [CrossRef]
  21. Zhang, N.; Hu, Z.; Dai, D.; Dang, S.; Yao, M.; Zhou, Y. Unit Commitment Model in Smart Grid Environment Considering Carbon Emissions Trading. IEEE Trans. Smart Grid 2016, 7, 420–427. [Google Scholar] [CrossRef]
  22. Luo, Y.; Zhang, X.; Yang, D.; Sun, Q. Emission Trading Based Optimal Scheduling Strategy of Energy Hub with Energy Storage and Integrated Electric Vehicles. J. Mod. Power Syst. Clean Energy 2020, 8, 267–275. [Google Scholar] [CrossRef]
  23. Luo, Y.; Yang, D.; Yin, Z.; Zhou, B.; Sun, Q. Optimal Configuration of Hybrid-energy Microgrid Considering the Correlation and Randomness of the Wind Power and Photovoltaic Power. IET Renew. Power Gener. 2020, 14, 616–627. [Google Scholar] [CrossRef]
  24. Wang, L.; Dong, H.; Lin, J.; Zeng, M. Multi-Objective Optimal Scheduling Model with IGDT Method of Integrated Energy System Considering Ladder-Type Carbon Trading Mechanism. Int. J. Electr. Power Energy Syst. 2022, 143, 108386. [Google Scholar] [CrossRef]
  25. Wang, R.; Wen, X.; Wang, X.; Fu, Y.; Zhang, Y. Low Carbon Optimal Operation of Integrated Energy System Based on Carbon Capture Technology, LCA Carbon Emissions and Ladder-Type Carbon Trading. Appl. Energy 2022, 311, 118664. [Google Scholar] [CrossRef]
  26. Zhang, X.; Liu, X.; Zhong, J.; Gao, W. Electricity–Gas-integrated Energy Planning Based on Reward and Penalty Ladder-type Carbon Trading Cost. IET Gener. Transm. Distrib. 2019, 13, 5263–5270. [Google Scholar]
  27. Zhang, Y.; Han, Y.; Liu, D.; Dong, X. Low-Carbon Economic Dispatch of Electricity-Heat-Gas Integrated Energy Systems Based on Deep Reinforcement Learning. J. Mod. Power Syst. Clean Energy 2023, 11, 1827–1841. [Google Scholar] [CrossRef]
  28. Gao, L.; Fei, F.; Jia, Y.; Wen, P.; Zhao, X.; Shao, H.; Feng, T.; Huo, L. Optimal Dispatching of Integrated Agricultural Energy System Considering Ladder-Type Carbon Trading Mechanism and Demand Response. Int. J. Electr. Power Energy Syst. 2024, 156, 109693. [Google Scholar] [CrossRef]
Figure 1. Structure of the IES.
Figure 1. Structure of the IES.
Energies 17 04705 g001
Figure 2. Comprehensive structure of the constraints.
Figure 2. Comprehensive structure of the constraints.
Energies 17 04705 g002
Figure 3. IES program flow chart.
Figure 3. IES program flow chart.
Energies 17 04705 g003
Figure 4. The power generated by different load in the IES.
Figure 4. The power generated by different load in the IES.
Energies 17 04705 g004
Figure 5. Prices for energy purchased from the upper grid and gas grid: (a) electricity prices, (b) gas prices.
Figure 5. Prices for energy purchased from the upper grid and gas grid: (a) electricity prices, (b) gas prices.
Energies 17 04705 g005
Figure 6. Initial and time-of-use tariffs for three loads: (a) gas prices, (b) heat prices, and (c) electricity prices.
Figure 6. Initial and time-of-use tariffs for three loads: (a) gas prices, (b) heat prices, and (c) electricity prices.
Energies 17 04705 g006
Figure 7. Change in load before and after DR in Scenario 4: (a) change in gas load, (b) change in heat load, and (c) change in electric load.
Figure 7. Change in load before and after DR in Scenario 4: (a) change in gas load, (b) change in heat load, and (c) change in electric load.
Energies 17 04705 g007
Figure 8. Balance of each power under Scenario 4: (a) gas power balance, (b) heat power balance, and (c) electric power balance.
Figure 8. Balance of each power under Scenario 4: (a) gas power balance, (b) heat power balance, and (c) electric power balance.
Energies 17 04705 g008
Figure 9. Balance of each power under Scenario 3: (a) gas power balance, (b) heat power balance, and (c) electric power balance.
Figure 9. Balance of each power under Scenario 3: (a) gas power balance, (b) heat power balance, and (c) electric power balance.
Energies 17 04705 g009
Table 1. Parameters of energy conversion equipment.
Table 1. Parameters of energy conversion equipment.
Equipment NamePower Upper and Lower Limits (kW)Energy Conversion EfficiencyOperational Costs (CNY/kW)Maximum Ramp Rate (kW/min)
GT[0, 4000]0.3/0.450.0415
GB[0, 1500]0.80.02515
EB[0, 1000]0.80.02510
P2G[0, 600]0.60.025
ORC[0, 600]0.80.0155
WHB[0, 500]0.80.0155
Table 2. Parameters of energy storage equipment.
Table 2. Parameters of energy storage equipment.
Equipment NameInitial Capacity (kWh)Maximum Capacity (kWh)Maximum Charge/Discharge Power (kW)Charge/Discharge EfficiencyOperational Costs (CNY/kW)
ES804003000.95/0.900.018
GS504003000.95/0.900.016
HS504003000.95/0.900.016
Table 3. Carbon trading-related parameters.
Table 3. Carbon trading-related parameters.
ParameterValueParameterValue
ε e 424 g/kwh ε p 968 g/kwh
ε h 432 g/kwh ε g 504 g/kwh
ε g l o a d 210 g/kwh ε q 320 g/kwh
υ 2000 kg γ 0.25 CNY/kg
ρ 0.25 ω 0.2
Table 4. Costs and carbon emissions for the four scenarios.
Table 4. Costs and carbon emissions for the four scenarios.
ScenarioTotal Cost (CNY)Energy Purchase Cost (CNY)Operation and Maintenance Cost (CNY)Carbon Trading Cost (CNY)Compensation Cost (CNY)Carbon Emissions (kg)
150,198.4842,943.323106.014149.15041,728.50
249,091.9041,846.352855.853908.94480.7640,616.06
348,973.4742,582.872989.673400.93040,085.71
447,875.1141,422.732833.513258.62360.2538,715.46
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, L.; Yang, S.; Chen, N.; Gao, J. Integrated Energy System Dispatch Considering Carbon Trading Mechanisms and Refined Demand Response for Electricity, Heat, and Gas. Energies 2024, 17, 4705. https://doi.org/10.3390/en17184705

AMA Style

Gao L, Yang S, Chen N, Gao J. Integrated Energy System Dispatch Considering Carbon Trading Mechanisms and Refined Demand Response for Electricity, Heat, and Gas. Energies. 2024; 17(18):4705. https://doi.org/10.3390/en17184705

Chicago/Turabian Style

Gao, Lihui, Shuanghao Yang, Nan Chen, and Junheng Gao. 2024. "Integrated Energy System Dispatch Considering Carbon Trading Mechanisms and Refined Demand Response for Electricity, Heat, and Gas" Energies 17, no. 18: 4705. https://doi.org/10.3390/en17184705

APA Style

Gao, L., Yang, S., Chen, N., & Gao, J. (2024). Integrated Energy System Dispatch Considering Carbon Trading Mechanisms and Refined Demand Response for Electricity, Heat, and Gas. Energies, 17(18), 4705. https://doi.org/10.3390/en17184705

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop