Next Article in Journal
Cold-Pressed Grape Seed Oil Encapsulation Using a Submerged-Nozzle Dispersion Encapsulation Process
Next Article in Special Issue
Future Prospects of MeOH and EtOH Blending in Gasoline: A Comparative Study on Fossil, Biomass, and Renewable Energy Sources Considering Economic and Environmental Factors
Previous Article in Journal
Evaluation and Simulation Study of Rock Drillability Grades in the Tazhong-Tabei Block
Previous Article in Special Issue
The Evolving Technological Framework and Emerging Trends in Electrical Intelligence within Nuclear Power Facilities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Economic Dispatch of Integrated Energy Systems Considering Wind–Photovoltaic Uncertainty and Efficient Utilization of Electrolyzer Thermal Energy

1
Engineering Research Center of Renewable Energy Power Generation and Grid-Connected Control, Ministry of Education, Xinjiang University, Urumqi 830017, China
2
Electric Power Research Institute of State Grid Xinjiang Electric Power Co., Ltd., Urumqi 830011, China
*
Authors to whom correspondence should be addressed.
Processes 2024, 12(8), 1627; https://doi.org/10.3390/pr12081627
Submission received: 25 June 2024 / Revised: 29 July 2024 / Accepted: 30 July 2024 / Published: 2 August 2024

Abstract

:
Currently, high levels of output stochasticity in renewable energy and inefficient electrolyzer operation plague IESs when combined with hydrogen energy. To address the aforementioned issues, an IGDT-based economic scheduling strategy for integrated energy systems is put forth. Firstly, this strategy establishes an IES consisting of coupled electricity, heat, hydrogen, and gas taking the hydrogen production electrolyzer’s thermal energy utilization into account. Second, to minimize the system’s overall operating costs, a deterministic scheduling model of the IES is built by taking into account the stepped carbon trading mechanism and the integrated demand response. Lastly, an optimal dispatch model is built using the information gap decision theory under the two strategies of risk aversion and risk seeking, taking into account the uncertainty of renewable energy generation. CPLEX is the solver used to solve the proposed model. After taking into account the effective use of thermal energy from the electrolyzer and loads demand response, the results show that the system carbon emission is reduced by 2597.68 kg and the operating cost is lowered by 44.65%. The IES scheduling model based on IGDT can effectively manage costs while maintaining system risk control, all while accommodating decision-makers’ varying risk preferences. This study can provide a useful reference for the research related to the scheduling of the IES low-carbon economy.

1. Introduction

Promoting the use of renewable energy sources and developing various uses for hydrogen energy is crucial from a practical standpoint when it comes to the “dual-carbon” aim and the green energy transition [1]. Several heterogeneous energy sources and conversion equipment can be combined by integrated energy systems (IESs) to satisfy the various demands of loads and offer a solution to encourage the use of energy efficiently and renewable energy sources [2,3].
Because of its high energy density, huge storage capacity, and high utilization rate, hydrogen energy can be used as a long-term energy storage system [4,5]. Simultaneously, IES carbon reduction goals can be achieved by the use of hydrogen energy, which is clean and environmentally benign [6]. Reference [7] proposes a wind–hydrogen energy system model, consisting of a wind turbine power generation system and a hydrogen storage system, and optimizes the configured capacity of the hydrogen storage system with the objective function of minimizing the system investment cost, operation, and maintenance cost. Reference [8] proposes an integrated energy system integrating hydrogen and electricity, considering both the operating costs of the system and the environmental costs associated with the buyback for constructing an optimal scheduling model for the coupled electricity–hydrogen–carbon system. Reference [9], an integrated energy system with wind turbines (WT), solar panels, and hydrogen storage systems, was built. A multi-objective day-ahead optimal scheduling model for the system was established, considering both the system’s operating costs and environmental costs comprehensively. NSGA-Ⅱ was then used to solve the problem. Reference [10] reduces the carbon emissions of the system by diversifying the utilization of hydrogen energy within the framework of green certificate trading mechanisms, focusing on energy applications. Reference [11] built an integrated electricity, heat, hydrogen, and gas energy system to lower carbon emissions and increase the use of renewable energy sources. This study did just that by utilizing hydrogen energy to enhance the system’s use of renewable energy sources and to realize the synergistic interaction between various energy sources. Utilizing the extended duration and high capacity that hydrogen energy storage can provide, the aforementioned study elaborates on integrating hydrogen energy storage into IESs to achieve hydrogen–electricity hybridization. However, during operation, the hydrogen fuel cell (HFC) and electrolyzer (EL) will unavoidably produce heat, and none of the literature [7,8,9,10,11] addressed the use of this portion of the heat, leading to energy waste in the system.
Reference [12] examined various methods of utilizing hydrogen energy, along with the physical attributes of each type of hydrogen storage mechanism. For the integrated thermoelectric hydrogen system, a two-layer optimization model is suggested, with the lower layer optimizing the system operation scheme and the upper layer determining the equipment capacity of the hydrogen energy system. Reference [13] describes the construction of an optimal scheduling model for a wind–photovoltaic–electricity–hydrogen storage integrated energy system. This model takes into account waste heat recovery from the HFC to promote the application of clean energy, such as wind–photovoltaic–hydrogen energy, within the park’s integrated energy system. Reference [14] suggests designing an integrated energy system in conjunction with hydrogen energy to achieve a synergistic supply of heat, cold, and electricity. It also establishes the ideal configuration for the integrated energy system by thoroughly assessing its energy, economics, environment, and grid dependence. The waste heat recovery from fuel cells was examined in reference [15], which also created an IES for combined biomass–photovoltaic–hydrogen energy and assessed the system’s performance from the energy–exergy–economy–environment standpoint.
As of right now, fuzzy optimization [16], stochastic optimization [17], and robust optimization [18] are the three primary types of approaches used to address the internal uncertainties of IESs. To solve numerous uncertainty problems in wind power integrated multi-energy systems, reference [19,20] developed an IES scheduling model based on fuzzy planning. However, fuzzy optimization is inherently subjective due to the necessity of selecting the membership function to describe uncertainty and its potential implications. A stochastic optimization approach was presented in reference [21] to address the unpredictability of renewable energy units in a microgrid and enhance system performance. A stochastic optimization-based optimal scheduling model for a combined wind–hydro system was created in reference [22] to address the deterministic nature of wind power output. However, stochastic optimization heavily relies on probabilistic models of uncertain quantities because it is grounded in probabilistic analysis. To address the uncertainties surrounding renewable energy sources, reference [23] created a coordination network and energy planning model for an IES utilizing robust optimization. A two-stage resilient optimal configuration model is proposed in reference [24] that examines the effects of various system uncertainties on both the configuration cost and the system’s running costs. References [25,26] employ a two-stage distributionally robust optimization approach to manage uncertainties in renewable energy and load demand within the system. Based on the range of variations in a given uncertainty, robust optimization selects the worst perturbation scenarios, resulting in overly conservative outcomes and suboptimal system economics. In contrast to the previously mentioned techniques, Information Gap Decision Theory (IGDT) exhibits the ability to establish the intended objective, optimize the range of fluctuation for uncertain parameters, and demonstrate strong adaptability when the probability distribution or fuzzy affiliation function of uncertain parameters is not known [27,28].
This paper aims to address the aforementioned issues by enhancing the energy utilization efficiency of hydrogen energy equipment within the system and mitigating the impact of source-load uncertainty on IES scheduling. Firstly, the operation framework of the IES is established by considering demand response, the stepped carbon trading mechanism, and the thermoelectric characteristics of EL. Second, the IES low-carbon economic dispatch model was established with the goal of minimizing the system costs associated with energy purchase, stepped carbon trading, and penalties for wind and light abandonment. Lastly, the IGDT method was used to discuss the uncertainty of WT and PV, and the risk-seeking and risk aversion models were created based on the decision-makers’ respective levels of risk preference. The original model must be converted into a mixed-integer linear model and solved in MATLAB using the CPLEX solver because the generated model is nonlinear. The effectiveness of the proposed model is verified by setting up different scenarios to analyze the scheduling results.

2. Operational Framework for IES

2.1. Structure of IES

The IES’s operating framework is depicted in Figure 1, with an emphasis on the EL’s efficient use of thermal energy. The system integrates energy supply equipment with multiple energy sources to meet its internal energy requirements. The three primary components of the IES model developed in this study are the energy supply unit, energy conversion and storage unit, and energy consumption unit, as seen in Figure 1. The energy supply unit of the system is made up of wind turbine (WT), photovoltaic (PV), the electricity grid, and the gas network. The energy conversion unit contains an EL, HFC, methane reactor (MR), combined heat and power (CHP), and gas boiler (GB). The energy storage unit has battery energy storage (BES), a hydrogen storage tank (HST), thermal storage tank (TST), and gas storage tank (GST).
The majority of the energy used to produce hydrogen in the EL and electricity in HFC is lost as heat, making the current electricity–hydrogen–electricity conversion process inefficient. As a result, recovering the working waste heat from the EL and HFC using a waste heat recovery device is crucial to increasing their energy conversion efficiency and making the most of hydrogen energy’s low-carbon and clean qualities.
The IES’s demand for thermal energy will be met by CHP, GB, EL, and HFC through the realization of hydrogen energy equipment cogeneration. This would lower the GB and CHP’s need for natural gas during the peak hours of thermal energy demand. The heat storage tank efficiently increases the stability and flexibility of the system’s energy supply and can further lower the system’s carbon emissions if the heat generated by the system is more than the heat load requirement.

2.2. Mathematical Models

2.2.1. Electrolyzer

The EL is a crucial device for converting electrical energy into hydrogen by employing the chemical principle of water electrolysis. Throughout this process, a portion of the energy dissipates as heat. Hence, the waste heat produced during electrolysis can be captured with a waste heat recovery device to enhance the energy efficiency of the EL [29]. At the moment, alkaline anion exchange membranes, high-temperature solid oxide with alkaline water electrolysis, and proton exchange membrane water electrolysis are the three primary technical pathways for producing hydrogen from electrolyzed water. Due to its advantages of mature technology, low cost, and long service life, the alkaline EL has become the most widely used equipment for producing hydrogen. In this paper, an alkaline EL is used for water electrolysis to produce hydrogen. The alkaline EL mathematical model is represented by the following equation [30].
P EL ,   H 2 ( t ) = η EL P e , EL ( t ) P EL ,   h ( t ) = ( 1 η EL ) P e , EL ( t ) P e , EL min P e , EL ( t ) P e , EL max Δ P e , EL min P e , EL ( t + 1 ) P e , EL ( t ) Δ P e , EL max
where P e , EL ( t ) is the electrical energy input to the EL, η EL is the energy conversion efficiency of the EL, in this paper, the value is 0.87; P EL ,   H 2 ( t ) is the hydrogen production power of the EL; P EL ,   h ( t ) is the thermal energy output from the EL; P e , EL min and P e , EL max are the lower and upper limits of the electrical energy input to the EL, respectively; Δ P e , EL max , Δ P e , EL min are the upper and lower limits of the climbing power of the EL, respectively.

2.2.2. Hydrogen Fuel Cell

As a hydrogen–electricity conversion device within the system, HFCs can be utilized for combined heat and power generation, as they produce significant heat energy alongside electrical output. Given the adjustable electric and thermal energy conversion efficiencies of HFC, this study aims to enhance energy supply flexibility by constructing a mathematical model that incorporates an adjustable electric-to-thermal ratio [31]. The mathematical model is shown below [30].
P HFC ,   e ( t ) = η HFC e P H 2 ,   HFC ( t ) P HFC ,   h ( t ) = η HFC h P H 2 ,   HFC ( t ) P H 2 ,   HFC min P H 2 ,   HFC ( t ) P H 2 ,   HFC max Δ P H 2 ,   HFC min P H 2 ,   HFC ( t + 1 ) + P H 2 ,   HFC ( t ) Δ P H 2 ,   HFC max κ HFC min P HFC ,   h ( t ) / P HFC ,   e ( t ) κ HFC max
where P HFC ,   e ( t ) , P HFC ,   h ( t ) are the electrical energy and thermal energy output from the HFC, respectively; η HFC e , η HFC h are the power generation efficiency and thermal efficiency of the HFC, respectively, the values are 0.55, 0.4, respectively; P H 2 ,   HFC ( t ) is the hydrogen consumption power of HFC; P H 2 ,   HFC max ( t ) , P H 2 ,   HFC min ( t ) are the upper and lower limits of the hydrogen input to the HFC, respectively; Δ P H 2 ,   HFC max ( t ) , Δ P H 2 ,   HFC min ( t ) are the upper and lower limits of the climbing power of the HFC, respectively; κ HFC min and κ HFC max are the lower and upper limits of the thermoelectric ratio of the HFC.

2.2.3. Methane Reactor

The MR utilizes hydrogen from the EL and CO2 from the IES for methanation, which is a mathematical model represented as follows.
P MR ,   g ( t ) = η MR P H 2 ,   MR ( t ) P H 2 ,   MR min P H 2 , MR ( t ) P H 2 ,   MR max Δ P H 2 ,   MR min P H 2 ,   MR ( t + 1 ) P H 2 , MR ( t ) Δ P H 2 ,   MR max
where P MR ,   g ( t ) is the natural gas power output from the MR; η MR is the energy conversion efficiency of the MR, the value is 0.6; P H 2 ,   MR ( t ) is the hydrogen consumption power of the MR; P H 2 ,   MR max , P H 2 ,   MR min are the upper and lower limits of the hydrogen energy input to the MR, respectively; and Δ P H 2 ,   MR max , Δ P H 2 ,   MR min are the upper and lower limits of the creep power of the MR, respectively.

2.2.4. Adjustable Thermoelectric Ratio CHP

Natural gas is burned in CHP for cogeneration. Because conventional CHP may only operate in “heat by electricity” or “heat by electricity” modes, their fixed cogeneration ratio restricts their operational flexibility. Adjustable thermoelectric ratio CHP can adjust the thermoelectric ratio according to the real-time situation of thermal and electric loads by controlling the pumping ratio of the turbine and the angle of the inlet guide vane, the mathematical model of which is shown below.
P CHP ,   e ( t ) = η CHP e P g ,   CHP ( t ) P CHP ,   h ( t ) = η CHP h P g ,   CHP ( t ) P g ,   CHP min P g , CHP ( t ) P g ,   CHP max Δ P g ,   CHP min P g ,   CHP ( t + 1 ) P g ,   CHP ( t ) Δ P g ,   CHP max κ CHP min P CHP ,   h ( t ) / P CHP ,   e ( t ) κ CHP max
where P CHP ,   e ( t ) , P CHP ,   h ( t ) are the electrical and thermal energy output of the CHP; η CHP e , η CHP h are the electrical and thermal energy conversion efficiency of the CHP; P g ,   CHP ( t ) is the natural gas power input to the CHP; P g ,   CHP max , P g ,   CHP min are the upper and lower limits of the natural gas power input to the CHP; Δ P g ,   CHP max , Δ P g ,   CHP min are the upper and lower limits of the ramp power of the CHP; κ CHP max , κ CHP min are the upper and lower limits of the CHP’s thermoelectricity ratio, respectively.

2.2.5. Gas Boiler

The GB produces a large amount of heat energy by burning natural gas, which is mathematically modeled as shown in the following equation.
P G B , h ( t ) = η G B h P g , G B ( t ) P g , G B min P g , G B ( t ) P g , G B max Δ P g , G B min P g , G B ( t + 1 ) P g , G B ( t ) Δ P g , G B max
where P G B , h ( t ) stands for the heat energy output of the GB; η G B h for the heat efficiency of the GB, the value is 0.92; P g , G B ( t ) for the GB consumption of natural gas power; P g , G B min , P g , G B max for the lower and upper limits of the heat power output of the GB; Δ P g , G B min , Δ P g , G B max for the lower and upper limits of the climbing power of the GB.

2.2.6. Energy Storage System

Energy storage devices in the IES include BES, HST, GST, and TST, which are modeled using a generic approach due to the similarity of their energy charging and discharging processes, as shown in the following equation.
  W s ( t ) = W s ( t 1 ) ( 1 δ s ) + ( η s , cha P s , cha ( t ) P s , dis ( t ) η s , dis ( t ) ) Δ t
where s is the type of energy, W s ( t ) is the energy stored in the energy storage device; δ s is the energy self-loss rate; η s , cha , η s , dis are the charging and discharging efficiency of the energy storage device, respectively; Δ t is the time interval.

2.3. Stepped Carbon Trading Model

A carbon trading mechanism is a system of trading that permits the market-based purchase and sale of carbon emission rights [32]. To reduce greenhouse gas emissions, the mechanism can support the commercialization of carbon credits, which will optimize the energy mix of the IES. Each source of carbon emissions is first given a specific amount of carbon credits by the regulator. Producers generate and release energy based on the credits that have been assigned to them. Producers can sell the leftover credits on the carbon market if their actual carbon emissions are less than the credits they have been allotted; if their actual emissions are more than the credits they have been allotted, they must buy more carbon credits. The carbon emission right quota model, the real carbon emission model, and the stepped carbon emission trading model are the three key parts of the stepped carbon trading model [33].

2.3.1. Carbon Emission Allowance Model

The IES developed in this study has three different kinds of carbon emission sources: GB, CHP, and purchased power from the upper grid. This study’s unpaid quota carbon quota approach is predicated on the idea that all power purchases from the upper grid originate from coal-fired units that produce energy.
E IES = E e , bay + E CHP + E GB E e , bay = χ e t = 1 T P e , bay ( t ) E CHP = χ g t = 1 T ( P CHP , e ( t ) + P CHP , h ( t ) ) E GB = χ g t = 1 T P GB , h ( t )
where E IES , E e , bay , E CHP , E GB are the carbon credits of the IES, purchased power from the upper grid, CHP, and GB, respectively; χ e , χ g are the carbon credits per unit of electricity consumed by coal-fired units and per unit of gas consumed by natural gas-fired units, respectively; P e , bay ( t ) is the amount of power purchased from the upper grid; and T is the scheduling cycle.

2.3.2. Real Carbon Emissions Modeling

The real carbon emissions are predicted as indicated by the following equation, as part of the CO2 created by the IES during the hydrogen-to-natural-gas conversion process is absorbed by the MR.
E IES , a = E e , buy , a + E total , a E MR , a E e , buy , a = t = 1 T ( a 1 + b 1 P e , buy ( t ) + c 1 P e , buy 2 ( t ) ) E total , a = t = 1 T ( a 2 + b 2 P total ( t ) + c 2 P total 2 ( t ) ) P total ( t ) = P CHP ,   e ( t ) + P CHP ,   h ( t ) + P GB ,   h ( t ) E MR ,   a = t = 1 T ϖ P MR ,   g ( t )
where E IES , a , E e , buy , a are the actual carbon emission from the IES and power purchased from the higher grid, respectively; E total , a is the total actual carbon emission from CHP, GB, and MR; E MR , a is the amount of CO2 absorbed by the MR during its operation; P total ( t ) is the equivalent output power of CHP, GB, and MR; a 1 , b 1 , c 1 and a 2 , b 2 , c 2 are the parameters for calculating the carbon emissions of coal-fired units and natural gas-consuming energy equipment, respectively, and are the parameters for the CO2 absorbed during the conversion of natural gas by the MR.

2.3.3. Stepped Carbon Emissions Trading Model

The actual carbon allowances for the IES to participate in the carbon trading market can be derived from the IES’s carbon credit allowances and actual carbon emissions, as shown in the following equation.
E IES ,   t = E IES ,   a E IES
where E IES ,   t is the amount of carbon credits traded by the IES.
This study uses a ladder pricing technique to separate various purchasing intervals and further limit the carbon emissions of the IES. The price of the response interval increases with the amount of carbon permits purchased. The stepped carbon transaction costs are shown in the following equation.
C CO 2 = λ E I E S , t E I E S , t l λ ( 1 + α ) ( E I E S , t l ) + λ l l E I E S , t 2 l λ ( 1 + 2 α ) ( E I E S , t 2 l ) + λ ( 2 + α ) l 2 l E I E S , t 3 l λ ( 1 + 3 α ) ( E I E S , t 3 l ) + λ ( 3 + 3 α ) l 3 l E I E S , t 4 l λ ( 1 + 4 α ) ( E I E S , t 4 l ) + λ ( 4 + 4 α ) l E I E S , t 4 l  
where C CO 2 is the stepped carbon transaction cost; λ is the carbon transaction base price; l is the length of the carbon emission interval; and α is the price growth rate.

2.4. Demand Response Model

The loads in the IES consist of heat, electricity, and gas. The flexibility of gas loads is not taken into account in this article. Electricity and heat loads include base loads and transferable loads. The base load is an uncontrollable load, which responds totally to the user’s demand, and the system cannot change the mode and time of its energy usage. To ensure that the total loads of the entire cycle after the transfer stay the same as it did before, the transferable loads can be dynamically altered for each period of energy consumption.

2.4.1. Demand Response for Electric Loads

Demand response refers to regulating the system load in the power system through pricing mechanisms and incentives to optimize the load curve and ensure stable system operation. In order to achieve the goal of reducing or shifting electric loads, this study employs price-based demand response for electric loads, charging fixed tariffs differently according to the period of energy consumption, and assisting users in changing their power usage habits [34]. The relationship between energy prices and demand is depicted by a price elasticity matrix. The energy usage cycle is split into three time periods: peak, flat, and valley. In this paper, demand response is modeled using the price elasticity matrix method, where the price elasticity matrix is expressed as follows [35].
E = ε 11 ε 12 ε 1 t ε 21 ε 22 ε 2 t ε t 1 ε t 2 ε t t
ε i i = Δ L i L i ( Δ p i p i ) 1
ε i j = Δ L i L i ( Δ p j p j ) 1
where ε i i and ε i j are the coefficients of auto-elasticity and cross-elasticity, respectively; L i represents the quantity of electricity demanded in time i ; Δ L i denotes the relative increment of quantity demanded; p i and Δ p i stand for the tariffs and their increments in time i ; p j and Δ p j denote the tariffs and their increments in time j .
Customer demand for electric loads can be categorized into fixed electric loads and transferable electric loads. Combined with the price elasticity matrix shown in Equation (14), the demand response model for electric loads is presented below.
P e _ load ( t ) = P e _ load fhl ( t ) + P e _ load shift ( t ) = L i + L i E Δ p i p i 0 P e _ load shift ( t ) P e _ load ,   max shift ( t ) t = 1 T P e _ load shift ( t ) Δ t = W e _ load shift
where P e _ load fhl ( t ) represents the fixed electric load; P e _ load shift ( t ) denotes the transferable electric load; P e _ load ,   max shift ( t ) is the upper limit of the transferable electric load; W e _ load shift signifies the total amount of transferable electric load; and Δ t is the time step, set at 1 h.

2.4.2. Demand Response for Heat Loads

Hot water loads are the primary ones for which demand response is taken into account since heat loads involve temporal delay and fuzzy perception characteristics. The primary heat load regulator that keeps the water’s temperature within an acceptable range is temperature.
P h _ load ,   min ( t ) P h _ load ( t ) P h _ load ,   max ( t ) P h _ load ,   min ( t ) = γ ρ w V c ( t ) T h , min T h , i n Δ t P h _ load ,   max ( t ) = γ ρ w V c ( t ) T h , max T h , i n Δ t
where P h _ load ( t ) is the heat loads power; P h _ load ,   max ( t ) , P h _ load ,   min ( t ) are the upper and lower limits of the heat loads power, respectively; γ is the specific heat capacity of water; ρ w is the density of water; V c ( t ) is the volume of cold water injected; T h , max , T h , min are the upper and lower limits of water temperature, respectively; T h , i n is the initial water temperature.
The heat loads can be categorized into fixed heat loads and transferable heat loads with the expression shown in the following equation.
P h _ load ( t ) = P h _ load fhl ( t ) + P h _ load shift ( t ) 0 P h _ load shift ( t ) P h _ load ,   max shift ( t )
t = 1 T P h _ load shift ( t ) Δ t = W h _ load shift
where P h _ load fhl ( t ) is the fixed heat load; P h _ load shift ( t ) is the transferable heat load; P h _ load ,   max shift ( t ) is the upper limit of the transferable heat loads; W h _ load shift is the total amount of transferable heat loads.

3. IES Optimal Economic Scheduling Model

3.1. Objective Function

The primary objective of this study is to assess the overall operational cost of the IES, which includes considerations of the system’s energy procurement costs, stepped carbon trading expenses, and costs associated with wind and photovoltaic curtailment, as shown in the following equation.
C = min ( C buy + C CO 2 + C cut )
where C buy is the cost of purchasing energy; C CO 2 is the cost of carbon trading; C cut is the wind and photovoltaic curtailment cost.
The combined cost of purchased gas and electricity constitutes the system’s purchased energy costs, which are determined using the following equation.
C buy = C e ,   buy + C g ,   buy
C e ,   buy = t = 1 T α t P e , buy ( t ) C g ,   buy = t = 1 T β t P g , buy ( t )
where α t is the electricity price in period t ; P e , buy ( t ) is the quantity of electricity purchased; β t is the gas price in period t ; and P g , buy ( t ) is the quantity of gas purchased.
Equation (10) displays the system carbon trading cost calculation, while the following shows the wind and photovoltaic curtailment cost.
C cut = δ WT t = 1 T P WT , cut ( t ) + δ PV t = 1 T P PV , cut ( t )
where δ WT , δ PV are the unit cost of wind curtailment and unit cost of photovoltaic curtailment; P WT , cut ( t ) , P PV , cut ( t ) are the wind curtailment cost and photovoltaic curtailment cost, respectively.

3.2. Restrictive Condition

3.2.1. WT and PV Unit Output Constraints

0 P WT ( t ) P WT max 0 P PV ( t ) P PV max
where P WT max , P PV max are the upper limits of WT and PV output, respectively.

3.2.2. Energy Storage Unit Operational Constraints

W s ,   min W s ( t ) W s ,   max 0 x s ,   cha ( t ) + x s ,   dis ( t ) 1 x s ,   cha ( t ) P s ,   cha min P t s , cha ( t ) x s ,   cha ( t ) P s ,   cha max x s ,   dis ( t ) P s ,   dis min P t s ,   dis ( t ) x s ,   dis ( t ) P s ,   dis max x s ,   cha ( t ) P s ,   cha down P s ,   cha ( t ) P s ,   cha ( t 1 ) x s ,   cha ( t ) P s ,   cha up x s ,   dia ( t ) P s ,   dia down P s ,   dia ( t ) P s ,   dia ( t 1 ) x s ,   dai ( t ) P s ,   dia up W 1 s = W T s
where W s ,   max , W s ,   min are the upper and lower energy storage limits of the energy storage device; x s ,   cha , x s ,   dis are the charging flag bit and discharging flag bit of the energy storage device; P s ,   cha max , P s ,   cha min are the upper and lower energy storage power limits; P s ,   dis max , P s ,   dis min are the upper and lower energy discharge power limits; P s ,   cha up   , P s ,   cha down are the upper and lower energy storage climbing power limits; P s ,   dis up   , P s ,   dis down are the upper and lower energy discharge climbing power limits; W 1 s , W T s are the initial value and the end value of the energy storage device in a scheduling cycle, respectively.

3.2.3. Purchased Energy Power Constraints

0 P e ,   buy t P e ,   buy max 0 P g ,   buy t P g ,   buy max

3.2.4. Electric Power Balance Constraints

P e , buy t + P WT t + P PV t + P CHP , e t + P HFC , e t + P e , dis t = P e _ load fhl ( t ) + P e _ load shift ( t ) + P e , EL t + P e , cha t
where P e , cha t , P e , dis t are the electric storage charging power and discharging power, respectively.

3.2.5. Thermal Power Balance Constraints

P HFC , h t + P CHP , h t + P GB , h t + P EL , h t + P h , dis t = P h _ load fhl ( t ) + P h _ load shift ( t ) + P h , cha t
where P h , cha t , P h , dis t are the TST heat storage power and exothermic power, respectively.

3.2.6. Gas Power Balance Constraints

P g , buy t + P MR , g t + P g , dis t = P g , Load t + P g , cha t + P g , CHP t +   P g , GB t
where P g , cha t , P g , dis t are the GST filling power and deflating power, respectively; P g , Load t is for the gas loads.

3.2.7. Hydrogen Equilibrium Constraints

P EL ,   H 2 ( t ) + P H 2 , dis ( t ) = P H 2 , MR ( t ) + P H 2 , HFC ( t ) + P H 2 , cha ( t )
where P H 2 , cha t , P H 2 , dis t are the hydrogen charging power and hydrogen discharging power of the HST, respectively.

4. IGDT-Based IES Optimized Scheduling Model

This study uses IGDT to predict the uncertainty in the amplitude of wind power fluctuation, to address the stochastic and intermittent character of wind and solar power output. The aim is to study the effect of uncertain parameters on the IES and to quantify the error between the predicted and actual values of uncertain parameters, provided that the preset objectives are met. The equation below illustrates a typical IGDT model [36,37].
min C ( U , x ) s . t .   H ( U , x ) = 0 G ( U , x ) 0
where U is the uncertainty parameter; x is the decision variable; C ( U , x ) is the objective function; H ( U , x ) and G ( U , x ) are the equality constraints and inequality constraints, respectively.
This study characterizes the uncertainty of WT and PV output using the envelope constraint. The equation below displays the interval uncertainty model for solar and wind power.
U ( α WT , P WT ,   pr ( t ) ) = P WT ( t ) : ( 1 α WT ) P WT ,   pr ( t ) P WT ( t ) ( 1 + α WT ) P WT ,   pr ( t ) U ( α PV , P PV ,   pr ( t ) ) = P PV ( t ) : ( 1 α PV ) P PV ,   pr ( t ) P PV ( t ) ( 1 + α PV ) P PV ,   pr ( t )
where P WT ,   pr ( t ) , P PV ,   pr ( t ) are the predicted power of WT and PV; P WT ( t ) , P PV ( t ) are the actual power of WT and PV; α WT , α PV are the power fluctuation of WT and PV, respectively; the formula for both are shown below.
α WT = ( P WT ,   pr P WT ) P WT ,   pr α PV = ( P PV ,   pr P PV ) P PV ,   pr
In this paper, considering uncertainties from both WT and PV, IGDT leads to an increase in the number of objectives in the model when dealing with multiple uncertain variables. This consequently escalates the complexity of solving the model. Therefore, this paper employs a linear weighting approach to aggregate the uncertainties of WT and PV, as shown in Equation (32).
φ = κ WT α WT + κ PV α PV
where φ is the combined uncertainty; κ WT , κ PV are the weighting factors for WT uncertainty and PV uncertainty, respectively. The impact of wind uncertainty and PV uncertainty on the system is proportional to the total power they generate; specifically, the larger the total power generated, the greater the weighting factor. The formula for calculating the weighting factor is shown in the following equation.
κ WT = t = 1 T P WT ( t ) t = 1 T ( P WT ( t ) + P WT ( t ) ) κ PV = t = 1 T P PV ( t ) t = 1 T ( P PV ( t ) + P PV ( t ) )

4.1. IGDT Scheduling Model Based on Risk Aversion

The goal of the risk aversion approach is to enhance the degree of uncertainty while gaining the ability to minimize risk at the expense of increased operating costs. The IDGT scheduling model based on risk aversion is shown in the following equation.
max φ max C ( 1 + β RM ) C 0 P WT ture ( t ) U ( α WT , P WT ( t ) ) P PV ture ( t ) U ( α PV , P PV ( t ) ) H ( U , x ) = 0 G ( U , x ) 0
where C 0 is the optimal economic operating cost of the IES with uncertainty equal to 0; β RM is the cost deviation parameter under the risk aversion strategy.
The two-layer planning model in Equation (34) shows the maximum uncertainty when the costs do not exceed expectations, and the lower layer represents the operating costs that must not exceed the expected costs when the actual output of renewable energy varies within the uncertainty interval. Operating costs rise with increasing uncertainty. Therefore, the two-tier planning model can be changed to a single-tier planning model, as shown in the following equation.
max φ max C ( 1 + β RM ) C 0 P WT ture ( t ) = U ( α WT , P WT ( t ) ) P PV ture ( t ) = U ( α PV , P PV ( t ) ) H ( U , x ) = 0 G ( U , x ) 0

4.2. IGDT Scheduling Model Based on Risk Seeking

The goal of the risk-seeking strategy is to enable the decision-maker to use uncertainty as a tool for cost cutting. Considering the need to satisfy scheduling decision-makers while seeking opportunities for cost reduction despite increased risks, the IGDT optimization model is obtained, as shown in the following equation.
max φ max C ( 1 β OM ) C 0 P WT ture ( t ) U ( α WT , P WT ( t ) ) P PV ture ( t ) U ( α PV , P PV ( t ) ) H ( U , x ) = 0 G ( U , x ) 0
where β OM is the cost deviation parameter under the risk-seeking strategy.
By applying the model transformation procedure in Section 4.1, the model presented in Equation (36) can be changed into a single-layer optimization model, as indicated by the following equation:
max φ max C ( 1 β OM ) C 0 P WT ture ( t ) = U ( α WT , P WT ( t ) ) P PV ture ( t ) = U ( α PV , P PV ( t ) ) H ( U , x ) = 0 G ( U , x ) 0

4.3. Model Linearization Process

The economic dispatch model of the IES, developed in this study to consider uncertainties in wind and solar power generation and efficient heat utilization of the EL, is formulated as a mixed-integer nonlinear model. It needs to be transformed into a mixed-integer linear model. For Equation (8), linearization is performed using piecewise linear approximation. The linearization process is outlined as follows.
Step 1: The original function is divided into Q intervals by taking the Q + 1 segmentation point r 1 , r 2 , , r Q + 1 according to the desired accuracy.
Step 2: Add the Q + 1 continuous auxiliary variable w 1 , w 2 , , w Q + 1 with the Q binary-type auxiliary variable z 1 , z 2 , , z Q + 1 and satisfy Equation (38).
w 1 + w 2 + + w Q + 1 = 1 z 1 + z 2 + + z Q = 1 w 1 0 , w 2 0 , , w Q + 1 0 w 1 z 1 , w 2 z 1 + z 2 , , w Q + 1 z Q
Step 3: Replace the nonlinear function with the linear expression shown in Equation (39).
P e , buy = q = 1 Q + 1 w q r q E e , buy , a = q = 1 Q + 1 w q E e , buy , a ( r q )
Since Equation (10) is already a piecewise function, step 1 above can be omitted, and steps 2 and 3 can be combined for linear transformation.
The model proposed in this paper has been transformed into a mixed-integer linear model through segmented linearization. Therefore, the proposed model will be implemented in the MATLAB R2022b environment, and the CPLEX solver will be utilized to solve the problem. Compared to other solvers, CPLEX exhibits superior accuracy and robustness in solving linear programming problems, enabling rapid identification of optimal solutions for the optimization model.

5. Calculus Analysis

5.1. Parameter Setup

An arithmetic example is set up for verification in order to confirm the validity of the economic dispatch model proposed in this study considering the IES. Figure 2 and Figure 3 depict the load data and scenery forecasts within the IES. The optimized scheduling is executed in a 24 h cycle. Table 1 displays the operational characteristics of the system’s equipment; Table 2 displays the cost of time-shared power; Table 3 displays the real carbon emission model; and Table 4 displays the price of natural gas, which is 0.35 CNY/kw·h. The cost per unit of wind and photovoltaic curtailment is 0.2 CNY/kW·h, and the carbon emission allowance per unit of electricity consumed by a coal-fired unit is χ e = 0.798 kg/kW·h. Similarly, the carbon emission allowance per unit of natural gas consumed by a natural gas-fired unit is χ g = 0.385 kg/kW·h [30,38].
The following six scenarios are set up for comparison in order to assess the superiority of the IES scheduling model put forward in this research.
Scenario 1: The IES does not take into account the working waste heat utilization of the EL, the demand response, or the uncertainty of the scenery output;
Scenario 2: The IES considers the utilization of operating waste heat from the EL without considering the demand response and without considering the uncertainty of the scenery output;
Scenario 3: The IES does not take into account the working waste heat utilization of the EL and considers the demand response without taking into account the scenario output uncertainty;
Scenario 4: The IES considers the utilization of operating waste heat from the electrolysis tanks, taking into account the demand response and disregarding the scenery output uncertainty;
Scenario 5: The IES considers the use of operating waste heat from the EL, considers the demand response, takes into account the WT and PV output uncertainty, and uses a risk aversion model for optimization;
Scenario 6: The IES considers the use of operating waste heat from the EL, considers demand response, takes into account the WT and PV output uncertainty, and uses a risk-seeking type for optimization.

5.2. Analysis of Deterministic Model Simulation Results

Scenarios 1 to 4 are solved under the deterministic model, and the results obtained are shown in Table 4.

5.2.1. Analysis of the Impact of EL Heat Use on the IES

The effective use of heat energy in the EL is taken into account in Scenarios 2 and 4, as indicated in Table 4. In comparison to Scenario 1, Scenario 2 lowers carbon emissions by 1413.76 kg, lowers the cost of carbon trading by 37.64%, lowers the cost of buying gas by 27.96%, raises the cost of buying electricity by 12.9%, and lowers the overall cost by 26.59%. Scenario 4 represents a reduction of 43.52% in carbon emissions, 47.58% in carbon trading expenses, 28.45% in gas purchase expenses, and 31.51% in total expenditures as compared to Scenario 3.
The efficient utilization of heat energy from the EL effectively reduces system carbon emission costs, as well as minimizes expenses associated with purchasing gas and penalties for WT and PV curtailment. This leads to a significant reduction in the total operating costs of the system.
Figure 4 and Figure 5 depict the scheduling outcomes for Scenarios 1 and 2, respectively. Figure 4b indicates that the primary heat production equipment of the IES are the HFC and the CHP. When the TST releases heat energy and the GB is used as a backup heat source, the heat demand of the IES is comparatively high between the hours of 1 and 10. The IES’s thermal energy requirement is mostly satisfied by CHP, HFC, and EL, as shown in Figure 5b. During the 2–7 h when there is a high demand for thermal energy, the GB only produces a tiny amount of energy.
The primary cause of the reduction in Scenario 2 carbon emissions and system expenses is the consideration of the EL’s heat utilization, which lowers the output of the system’s gas-consuming equipment while simultaneously increasing the EL’s efficiency of use and electrical energy demand, both of which enhance the system’s capacity to use wind energy while raising the system’s power purchases.

5.2.2. Analysis of Demand Response on the IES

Figure 6 and Figure 7 depict the variations in thermal and electric loads following IES demand response. After considering price-based demand response for electric loads, users can adjust their energy consumption strategies based on electricity prices during different time periods. This allows them to select suitable time periods to meet their energy requirements, thereby enhancing both their own energy efficiency and the overall power utilization rate of the system. The electric loads undergo a temporal shift due to price-based demand response; however, the overall loads before and after the shift remain unchanged. Figure 3 illustrates how a price-based demand response controls the electric loads profile. The lowest power costs occur during the 1–7 h and 24 h when the IES’s electric energy demand increases by 13.03% and 10.57%, respectively. Customers adjust their energy usage strategies to shift loads to periods with lower electricity prices during the 8–11 and 17–22 periods, as these are times when electricity prices are higher. As can be observed from Figure 8a, sufficient PV and WT generation occurs during the 12–16 h to increase the system’s wind energy consumption rate, leading to an increase in the electric load demand of the IES during this period following optimization.
This study sets the transferable heat loads to account for 10% of the total heat loads to ensure the comfort of users’ energy use, since users have greater reliability requirements for heat loads. The total loads in the loads curve’s peak hours drop by 6.32%, the total loads in the trough hours rise by 7.50%, and the peak-to-valley difference drops from 388.8 kW to 293.98 kW following the demand response of the thermal loads, as shown in Figure 7. This demonstrates that a tendency of peak shaving follows the demand response in the heat loads curve.
Table 4 illustrates that, when the demand response of the electric and heat loads is taken into consideration, Scenario 4 lowers its carbon emissions by 1183.92 kg when compared to Scenario 2, lowers the cost of bought gas by 3.90%, and lowers the total cost by 24.60%. When Figure 5b and Figure 8b are compared, it is evident that the IES, which takes into account the high desired utilization of EL heat, reduces the system’s natural gas demand because heat loads shift in the 2–6 h, which removes the need for GB outlets and increases the number of HFC outlets in the 12–16 h. It has been demonstrated that demand response can maximize the benefits of CHP equipment and improve the system’s capacity for carbon reduction.

5.3. IGDT-Based IES Scheduling Analysis

The deviation coefficients have an impact on the IES optimal scheduling findings that are based on IGDT. Table 5 presents the operating costs and uncertainty solution outcomes for two techniques: one based on a risk-seeking strategy and the other on a risk aversion strategy for decision optimization, each with varying deviation coefficients.
With a deviation coefficient of 0.02 as an example, when the IES’s dispatch is biased toward conservatism—that is, when risk aversion is selected—the corresponding uncertainty is 0.0184 and the operating cost is 4438.78 CNY. This means that when the real-time scenery output fluctuates within the positive 0.0184 range and below, it can guarantee that the IES’s operating cost is not greater than 4438.78 CNY at this deviation coefficient.
The corresponding uncertainty is 0.0186, and the operating cost is 4264.71 CNY when IES scheduling is biased toward opportunistic and a risk-seeking strategy is selected. In other words, when real-time scenery output fluctuates within the range of negative 0.0186 and above, it can be guaranteed that the IES operating cost is not higher than 4264.71 CNY at this deviation factor.
From a system operation perspective, using the operating cost of 4351.75 CNY at a deviation factor of 0 as a benchmark, it is evident that under a risk-seeking strategy, the minimum acceptable uncertainty of the IES increases with rising deviation coefficient, while the intervals of wind turbine and PV unit outputs continue to narrow. This risk-seeking behavior escalates system operational risk, concurrently reducing operational costs to a minimum of 3916.57 CNY. The compressed intervals reflect the balance between increased risk and lowered costs. Under the risk-seeking strategy, the decision-maker can release up to an upper limit of 435.18 CNY in the cost space. If the operating costs continue to decrease beyond this threshold, the system may fail due to changes in the operating mechanism stemming from excessive risk. The upper limit of the IES’s operating costs under the risk aversion strategy is 4786.92 CNY, with decision-makers reserving a maximum of 435.71 CNY for effective risk management.
Figure 9 and Figure 10 illustrate the trends in operational cost and uncertainty with deviation factor. An increase in the system’s scenic output uncertainty has a favorable effect on reducing the operating cost of the system under the risk-seeking strategy. Here, uncertainty is inversely related to operating cost and positively correlated with the deviation coefficient. Simultaneously, as uncertainty grows, so does the system’s risk.
The deviation factor rises together with the level of uncertainty and system operation expenses under the risk aversion strategy. The reason for this is that when utilizing a risk aversion scheduling technique, the decision-maker fights the impact of uncertainty on the system by raising operating costs since they feel that uncertainty may adversely move the objective under that direction of deviation.

6. Discussion

This paper focuses on enhancing the energy efficiency of an EL within an IES by recovering and utilizing its waste heat to meet heat load demands. It leverages the flexible characteristics of electric and thermal loads, allowing them to participate in demand response as flexible loads. This reduces pressure on energy supply during peak consumption hours, effectively improving the system’s economic and low-carbon characteristics. The study considers uncertainties associated with WT and PV, using IGDT theory to offer decision-makers with varying risk preferences a theoretical basis for selecting optimal scheduling strategies.
This paper serves as a valuable reference for research in the field of low-carbon economic dispatch in an IES, considering the impact of uncertainty. However, there are still aspects that require further research and expansion.
  • This paper presents a basic simulation of thermal energy recovery from an EL and an HFC. The IES optimal scheduling model developed in this study will be enhanced in the future by more accurately modeling the thermal energy utilization of the EL and HFC, and by developing a comprehensive approach to thermal energy recovery.
  • The paper does not consider the ladder gas price; however, in our future work, we plan to include it in the study and conduct a comprehensive analysis of the system’s carbon emissions and economic costs.
  • In the future, we will continue to investigate uncertainties on both the source and load sides in the IES, thoroughly evaluate the advantages and disadvantages of current uncertainty approaches, and more comprehensively integrate theoretical findings with real-world implementations.
  • In this paper, the optimal scheduling of the system under several typical day situations is not examined, as only one set of wind and load data is chosen. Different typical day scenarios exhibit varying wind and solar production, along with cooling, heating, and electric demands. These variations lead to different IES scheduling strategies. To gain a deeper understanding of IES scheduling algorithms across various typical day scenarios, we intend to expand the dataset in future work.
  • In this paper, day-ahead scheduling has been conducted with a time scale of 1 h. Future research will explore intraday scheduling, incorporating real-time phases to optimize the IES across multiple time scales.

7. Conclusions

The integrated demand response of loads and the effective use of thermal energy from electrolysis tanks are taken into account in this study’s proposal for an IES low-carbon economic dispatch strategy. An IES dispatch model is constructed based on the IGDT under the two strategies of risk aversion and risk seeking, which can fulfill the risk preferences of different decision-makers, in order to balance the equilibrium of IES risk management and operating expenses. The simulation leads to the following conclusions:
  • The efficient use of waste heat generated by the EL in the IES can improve energy efficiency, relieve pressure on the heat supply during peak heat load hours, and lessen the IES’s reliance on GB and CHP, all of which lower the system’s carbon emissions. When using the EL thermal energy-efficient utilization model, the system’s carbon output is decreased by 1413.76 kg and its operating costs are lowered by 26.59% when compared to the standard EL utilization model;
  • When heat, electricity, gas, and hydrogen energy sources are coupled internally in the IES, the electric loads and heat loads act as flexible loads that can engage in demand response. This allows the electric loads to be transferred in a way that is reasonable based on the price of electricity, effectively reducing the heat loads’ peak and valley differences and increasing the system’s energy use flexibility. After considering the demand response of the loads, the carbon emission of the system is reduced by 1138.92 kg and the operating cost is reduced by 24.60%;
  • IGDT is employed to simulate the uncertainty of WT and PV fluctuations and to quantify the error between predicted and actual values of WT and PV outputs. The IGDT-based IES optimal scheduling model can derive two scheduling strategies, risk seeking and risk aversion, that satisfy the risk preferences of decision-makers. This model achieves a balance between IES risk control and operational costs while effectively addressing the stochastic nature of WT and PV outputs.

Author Contributions

Methodology, L.X. and A.B.; formal analysis, Y.K.; investigation, Y.Z.; data curation, J.L.; writing—original draft, J.L.; writing—review and editing, Z.Y. and L.X.; supervision, Z.Y. and W.L.; project administration, W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article, and further inquiries can be directed to the corresponding authors.

Conflicts of Interest

Authors Ji Li, Lei Xu, Yuying Zhang, Yang Kou, Weile Liang and Alihan Bieerke were employed by the Electric Power Research Institute of State Grid Xinjiang Electric Power Co., Ltd. The remaining author declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. Electric Power Research Institute of State Grid Xinjiang Electric Power Co., Ltd. had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Nomenclature

AbbreviationsFull name
IESintegrated energy systems
IGDTinformation gap decision theory
HFChydrogen fuel cell
ELelectrolyzer
WTwind turbine
PVphotovoltaic
MRmethane reactor
CHPcombined heat and power
GBgas boiler
BESbattery energy storage
HSThydrogen storage tank
TSTthermal storage tank
GSTgas storage tank
Technical Parameters
ParametersValueParametersValue
η EL 0.87 η G B h 0.95
η HFC e 0.55 η s , cha 0.8
η HFC h 0.4 η s , dis 0.8
η MR 0.6 χ e 0.798 kg/kW·h
η CHP e 0.4 χ g 0.385 kg/kW·h
η CHP h 0.52 δ WT , δ PV 0.2 CNY/kW·h
Variable Descriptions
VariablesDescription
P EL ,   H 2 ( t ) the hydrogen production power of the EL (kW)
P e ,   EL ( t ) the electrical energy input to the EL (kW)
P EL ,   h ( t ) the thermal energy output from the EL (kW)
P HFC ,   e ( t ) the electrical energy from the HFC (kW)
P H 2 ,   HFC ( t ) the hydrogen input to the HFC (kW)
P HFC ,   h ( t ) the thermal energy output from the HFC (kW)
P MR ,   g ( t ) the natural gas power output from the MR (kW)
P H 2 ,   MR ( t ) the hydrogen energy input to the MR (kW)
P CHP ,   e ( t )   P CHP ,   h ( t ) the electrical and thermal energy output of the CHP (kW)
P g ,   CHP ( t ) the natural gas power input to the CHP (kW)
P G B , h ( t ) the heat energy output of the GB (kW)
P g , G B ( t ) the GB consumption of natural gas power (kW)
W s ( t ) the energy stored in the energy storage device (kW·h)
E IES the carbon credits of the IES (kg)
E e , bay the carbon credits of purchased power (kg)
E CHP the carbon credits of CHP (kg)
E GB the carbon credits of GB (kg)
E IES , a the actual carbon emission from the IES (kg)
E e , buy , a the actual carbon emission from power purchased (kg)
E MR , a the amount of CO2 absorbed by the MR (kg)
E total , a the total actual carbon emission from CHP, GB, and MR (kg)
P total ( t ) the equivalent output power of CHP, GB, and MR (kW)
P e _ load ( t ) , P h _ load ( t ) the electrical and heat loads (kW)
P e _ load shift ( t ) , P h _ load shift ( t ) the transferable electrical load and heat load (kW)
P e _ load fhl ( t ) , P h _ load fhl ( t ) the fixed electrical and heat loads (kW)
W e _ load shift , W h _ load shift the total amount of transferable electrical load and heat loads (kW)
C the overall operational cost (CNY)
C buy the cost of purchasing energy (CNY)
C CO 2 the cost of stepped carbon trading (CNY)
C cut WT and PV curtailment cost (CNY)

References

  1. Liu, D.; Luo, Z.; Qin, J.; Wang, H.; Wang, G.; Li, Z.; Zhao, W.; Shen, X. Low-carbon dispatch of multi-district integrated energy systems considering carbon emission trading and green certificate trading. Renew. Energy 2023, 218, 119312. [Google Scholar] [CrossRef]
  2. Ma, K.; Zhang, R.; Yang, J.; Song, D. Collaborative optimization scheduling of integrated energy system considering user dissatisfaction. Energy 2023, 274, 127311. [Google Scholar] [CrossRef]
  3. Wu, J.; Gu, J.; Liu, S.; Jin, Z. Strategies for Improving Resilience of Regional Integrated Energy Systems in the Prevention-Resistance Phase of Integration. Prot. Control Mod. Power Syst. 2023, 8, 29. [Google Scholar] [CrossRef]
  4. He, J.; Wu, Y.; Yong, Z.; Tan, Q.; Liu, F. Bi-level optimization of a near-zero-emission integrated energy system considering electricity-hydrogen-gas nexus, a two-stage framework aiming at economic and environmental benefits. Energy Convers. Manag. 2022, 274, 116434. [Google Scholar] [CrossRef]
  5. Li, Z.; Xia, Y.; Bo, Y.; Wei, W. Optimal planning for electricity-hydrogen integrated energy system considering multiple timescale operations and representative time-period selection. Appl. Energy 2024, 362, 122965. [Google Scholar] [CrossRef]
  6. Chang, P.; Li, C.; Zhu, Q.; Zhu, T.; Shi, J. Optimal scheduling of electricity and hydrogen integrated energy system considering multiple uncertainties. iScience 2024, 27, 109654. [Google Scholar] [CrossRef] [PubMed]
  7. Dong, Y.; Han, Z.; Li, C.; Ma, S.; Ma, Z. Research on the optimal planning method of hydrogen-storage units in wind–hydrogen energy system considering hydrogen energy source. Energy Rep. 2023, 9, 1258–1264. [Google Scholar] [CrossRef]
  8. Bu, F.; Wang, S.; Bai, H.; Wang, Y.; Yu, L.; Liu, H. An integrated demand response dispatch strategy for low-carbon energy supply park considering electricity–hydrogen–carbon coordination. Energy Rep. 2023, 9, 1092–1101. [Google Scholar] [CrossRef]
  9. Fang, R. Multi-objective optimized operation of integrated energy system with hydrogen storage. Int. J. Hydrogen Energy 2019, 44, 29409–29417. [Google Scholar]
  10. Rong, H.; Kuang, H. A cooperative model of photovoltaic and electricity-to-hydrogen including green certificate trading under the conditional value at risk. Glob. Energy Interconnect. 2023, 6, 403–417. [Google Scholar] [CrossRef]
  11. 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]
  12. Zeng, G.; Liu, M.; Lei, Z.; Huang, X. Bi-level robust planning of hydrogen energy system for integrated electricity–heat–hydrogen energy system considering multimode utilization of hydrogen. Energy 2024, 303, 132029. [Google Scholar] [CrossRef]
  13. Zhan, Z.; Xu, J.; Zhang, T.; Kong, X.; Li, H.; Zhang, B. Optimal scheduling of integrated wind-photovoltaic-hydrogen energy system considering hydrogen application and waste heat recovery. Energy Rep. 2024, 11, 3684–3694. [Google Scholar] [CrossRef]
  14. Wang, J.; Xue, K.; Guo, Y.; Ma, J.; Zhuo, X.; Liu, M.; Yan, J. Multi-objective capacity programming and operation optimization of an integrated energy system considering hydrogen energy storage for collective energy communities. Energy Convers. Manag. 2022, 268, 116057. [Google Scholar] [CrossRef]
  15. Wang, Y.; Shi, L.; Song, M.; Jia, M.; Li, B. Evaluating the energy-exergy-economy-environment performance of the biomass-photovoltaic-hydrogen integrated energy system based on hybrid multi-criterion decision-making model. Renew. Energy 2024, 224, 120220. [Google Scholar] [CrossRef]
  16. Yang, X.; Wang, X.; Deng, Y.; Mei, L.; Deng, F.; Zhang, Z. Integrated energy system scheduling model based on non-complete interval multi-objective fuzzy optimization. Renew. Energy 2023, 218, 119289. [Google Scholar] [CrossRef]
  17. Alvarez, G.E. Stochastic optimization considering the uncertainties in the electricity demand, natural gas infrastructures, photovoltaic units, and wind generation. Comput. Chem. Eng. 2022, 160, 107712. [Google Scholar] [CrossRef]
  18. Li, F.; Wang, D.; Guo, H.; Zhang, J. Distributionally Robust Optimization for integrated energy system accounting for refinement utilization of hydrogen and ladder-type carbon trading mechanism. Appl. Energy 2024, 367, 123391. [Google Scholar] [CrossRef]
  19. Mohammadi, M.; Noorollahi, Y.; Mohammadi-Ivatloo, B. Fuzzy-based scheduling of wind integrated multi-energy systems under multiple uncertainties. Sustain. Energy Technol. Assess. 2020, 37, 100602. [Google Scholar] [CrossRef]
  20. Gholami, K.; Karimi, S.; Rastgou, A. Fuzzy risk-based framework for scheduling of energy storage systems in photovoltaic-rich networks. J. Energy Storage 2022, 52, 104902. [Google Scholar] [CrossRef]
  21. Bornapour, M.; Hooshmand, R.A.; Khodabakhshian, A.; Parastegari, M. Optimal stochastic coordinated scheduling of proton exchange membrane fuel cell-combined heat and power, wind and photovoltaic units in micro grids considering hydrogen storage. Appl. Energy 2017, 202, 308–322. [Google Scholar] [CrossRef]
  22. Lu, N.; Su, C.; Guo, C.; Wang, P.; Wu, H.; Sui, Q. Stochastic optimal scheduling of wind power and pumped-storage hydropower complementary systems with multiple uncertainties. J. Energy Storage 2024, 78, 110060. [Google Scholar] [CrossRef]
  23. Chen, L.; Yang, D.; Cai, J.; Yan, Y. Robust optimization based coordinated network and source planning of integrated energy systems. Int. J. Electr. Power Energy Syst. 2024, 157, 109864. [Google Scholar] [CrossRef]
  24. Yang, M.; Liu, Y. A two-stage robust configuration optimization framework for integrated energy system considering multiple uncertainties. Sustain. Cities Soc. 2024, 101, 105120. [Google Scholar] [CrossRef]
  25. Fan, W.; Ju, L.; Tan, Z.; Li, X.; Zhang, A.; Li, X.; Wang, Y. Two-stage distributionally robust optimization model of integrated energy system group considering energy sharing and carbon transfer. Appl. Energy 2023, 331, 120426. [Google Scholar] [CrossRef]
  26. Shi, Z.; Yang, Y.; Xu, Q.; Wu, C.; Hua, K. A low-carbon economic dispatch for integrated energy systems with CCUS considering multi-time-scale allocation of carbon allowance. Appl. Energy 2023, 351, 121841. [Google Scholar] [CrossRef]
  27. Peng, B.; Kong, X.; Tian, C.; Zhang, F.; Ma, X. An IGDT-based a low-carbon dispatch strategy of urban integrated energy system considering intermittent features of renewable energy. Energy Rep. 2023, 10, 4390–4401. [Google Scholar] [CrossRef]
  28. Boroumandfar, G.; Khajehzadeh, A.; Eslami, M.; Syah, R.B. Information gap decision theory with risk aversion strategy for robust planning of hybrid photovoltaic/wind/battery storage system in distribution networks considering uncertainty. Energy 2023, 278, 127778. [Google Scholar] [CrossRef]
  29. Qiu, Y.; Li, Q.; Ai, Y.; Chen, W.; Benouzid, M.; Liu, S. Two-stage distributionally robust optimization-based coordinated scheduling of integrated energy system with electricity-hydrogen hybrid energy storage. Prot. Control Mod. Power Syst. 2023, 8, 1–14. [Google Scholar] [CrossRef]
  30. Wu, Q.; Li, C. Modeling and operation optimization of hydrogen-based integrated energy system with refined power-to-gas and carbon-capture-storage technologies under carbon trading. Energy 2023, 270, 126832. [Google Scholar] [CrossRef]
  31. Gao, J.; Meng, Q.; Liu, J.; Wang, Z. Thermoelectric optimization of integrated energy system considering wind-photovoltaic uncertainty, two-stage power-to-gas and ladder-type carbon trading. Renew. Energy 2024, 221, 119806. [Google Scholar] [CrossRef]
  32. Qin, M.; Yang, Y.; Zhao, X.; Xu, Q.; Yuan, L. Low-carbon economic multi-objective dispatch of integrated energy system considering the price fluctuation of natural gas and carbon emission accounting. Prot. Control Mod. Power Syst. 2023, 8, 61. [Google Scholar] [CrossRef]
  33. Zhang, J.; Liu, Z. Low carbon economic scheduling model for a park integrated energy system considering integrated demand response, ladder-type carbon trading and fine utilization of hydrogen. Energy 2024, 290, 130311. [Google Scholar] [CrossRef]
  34. Yang, H.; Zhang, X.; Ma, Y.; Zhang, D. Critical peak rebate strategy and application to demand response. Prot. Control Mod. Power Syst. 2021, 6, 1–14. [Google Scholar] [CrossRef]
  35. Li, W.; He, M.; Cai, T. Low-carbon economic dispatch of hydrogen-containing integrated energy system considering stepped demand response. Energy Rep. 2024, 11, 4222–4232. [Google Scholar] [CrossRef]
  36. Feng, B.; Fu, Y.; Huang, Q.; Ma, C.; Sun, Q.; Ronald, W. Multi-objective optimization of an integrated energy system with high proportion of renewable energy under multiple uncertainties. Energy Rep. 2023, 9, 695–701. [Google Scholar] [CrossRef]
  37. 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]
  38. Zhu, X.; Xue, J.; Hu, M.; Liu, Z.; Gao, X.; Huang, W. Low-carbon economy dispatching of integrated energy system with P2G-HGT coupling wind power absorption based on stepped Carbon emission trading. Energy Rep. 2023, 10, 1753–1764. [Google Scholar] [CrossRef]
Figure 1. IES Architecture.
Figure 1. IES Architecture.
Processes 12 01627 g001
Figure 2. Electric, heat, and gas load curves.
Figure 2. Electric, heat, and gas load curves.
Processes 12 01627 g002
Figure 3. WT and PV forecast curves.
Figure 3. WT and PV forecast curves.
Processes 12 01627 g003
Figure 4. Scenario 1 Scheduling Results. (a) Scenario 1 Electricity Scheduling Diagram. (b) Scenario 1 Heat Scheduling Diagram.
Figure 4. Scenario 1 Scheduling Results. (a) Scenario 1 Electricity Scheduling Diagram. (b) Scenario 1 Heat Scheduling Diagram.
Processes 12 01627 g004
Figure 5. Scenario 2 Scheduling Results. (a) Scenario 2 Electricity Scheduling Diagram. (b) Scenario 2 Heat Scheduling Diagram.
Figure 5. Scenario 2 Scheduling Results. (a) Scenario 2 Electricity Scheduling Diagram. (b) Scenario 2 Heat Scheduling Diagram.
Processes 12 01627 g005
Figure 6. Scenario 4 electric loads curve before and after demand response.
Figure 6. Scenario 4 electric loads curve before and after demand response.
Processes 12 01627 g006
Figure 7. Scenario 4 heat loads curve before and after demand response.
Figure 7. Scenario 4 heat loads curve before and after demand response.
Processes 12 01627 g007
Figure 8. Scenario 4 Scheduling Results. (a) Scenario 4 Electricity Scheduling Diagram. (b) Scenario 4 Heat Scheduling Diagram.
Figure 8. Scenario 4 Scheduling Results. (a) Scenario 4 Electricity Scheduling Diagram. (b) Scenario 4 Heat Scheduling Diagram.
Processes 12 01627 g008
Figure 9. Running costs and uncertainty under risk-seeking strategy.
Figure 9. Running costs and uncertainty under risk-seeking strategy.
Processes 12 01627 g009
Figure 10. Running costs and uncertainty under risk aversion strategy.
Figure 10. Running costs and uncertainty under risk aversion strategy.
Processes 12 01627 g010
Table 1. Main equipment parameters.
Table 1. Main equipment parameters.
DeviceCapacity (kW)Climbing Constraints (%)
EL50020
HFC25020
CHP60020
MR25020
GB30020
BES45020
TST45020
HST20020
GST20020
Table 2. Time-sharing tariff.
Table 2. Time-sharing tariff.
Times (h)Price of Electricity (CNY)
01:00–07:00, 23:00–24:000.38
08:00–11:00, 15:00–18:000.68
12:00–14:00, 19:00–22:001.20
Table 3. Modeling of real carbon emissions.
Table 3. Modeling of real carbon emissions.
Power Consumption TypeGas Consuming Type
a 1 b 1 c 1 a 2 b 2 c 2
36−0.380.00343−0.0040.001
Table 4. Deterministic model optimization results.
Table 4. Deterministic model optimization results.
Scenario 1Scenario 2Scenario 3Scenario 4
E IES (kg)4233.702819.942896.871636.02
C CO 2 (CNY)1212.64756.22780.27409.00
C e ,   buy (CNY)715.68808.2600
C g ,   buy (CNY)5695.114102.785510.463942.74
C cut (CNY)238.68104.6563.300
C (CNY)7862.115771.926354.034351.75
Table 5. Calculation of the coefficient of deviation.
Table 5. Calculation of the coefficient of deviation.
Risk-Seeking Strategy Risk Aversion Strategy
Operational Cost (CNY)UncertaintyDeviation FactorUncertaintyOperational Cost (CNY)
4351.750004351.75
4264.710.01860.020.01844438.78
4177.680.03760.040.03644525.82
4090.640.05830.060.05454612.85
4003.610.08510.080.07124699.89
3916.570.13270.100.08774786.92
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

Li, J.; Xu, L.; Zhang, Y.; Kou, Y.; Liang, W.; Bieerke, A.; Yuan, Z. Economic Dispatch of Integrated Energy Systems Considering Wind–Photovoltaic Uncertainty and Efficient Utilization of Electrolyzer Thermal Energy. Processes 2024, 12, 1627. https://doi.org/10.3390/pr12081627

AMA Style

Li J, Xu L, Zhang Y, Kou Y, Liang W, Bieerke A, Yuan Z. Economic Dispatch of Integrated Energy Systems Considering Wind–Photovoltaic Uncertainty and Efficient Utilization of Electrolyzer Thermal Energy. Processes. 2024; 12(8):1627. https://doi.org/10.3390/pr12081627

Chicago/Turabian Style

Li, Ji, Lei Xu, Yuying Zhang, Yang Kou, Weile Liang, Alihan Bieerke, and Zhi Yuan. 2024. "Economic Dispatch of Integrated Energy Systems Considering Wind–Photovoltaic Uncertainty and Efficient Utilization of Electrolyzer Thermal Energy" Processes 12, no. 8: 1627. https://doi.org/10.3390/pr12081627

APA Style

Li, J., Xu, L., Zhang, Y., Kou, Y., Liang, W., Bieerke, A., & Yuan, Z. (2024). Economic Dispatch of Integrated Energy Systems Considering Wind–Photovoltaic Uncertainty and Efficient Utilization of Electrolyzer Thermal Energy. Processes, 12(8), 1627. https://doi.org/10.3390/pr12081627

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