Integrated Demand Response for Micro-Energy Grid Accounting for Dispatchable Loads

: Micro-energy networks are the smallest element of integrated energy systems, and tapping into the integrated demand response potential of micro-energy networks is conducive to improving energy use efficiency and promoting the development of new energy sources on a large scale. This paper proposes a day-ahead integrated demand response strategy for micro-energy grid that takes into account the dispatchable loads. Considering the gradient use of thermal energy, a typical micro-energy grid structure including electricity, gas, medium-grade heat, low-grade heat, and cold energy is constructed, a comprehensive energy equipment model is established, and the refined scheduling models of the dispatchable loads are given. On this basis, with the operating economy of the micro-energy grid as the optimization objective, the integrated demand response strategies of tariff-type and incentive-type are proposed. Through case study analysis, it is verified that the proposed strategy can optimize the energy consumption structure of the micro-energy grid under the guidance of time-of-use tariffs, reducing the operating costs. The proposed strategy fully exploits the demand response potential of the micro-energy grid through the dispatchable loads and the multi-energy complementarity of electricity, heat, and cold, realizes the comprehensive coordination and optimization of source-network-load-storage, provides a larger peak-regulating capacity, and exhibits practical applicability in engineering.


Introduction
Energy is crucial for economic and social development.Promoting efficient energy use and improving new energy consumption capacity is key to achieving clean and lowcarbon electric energy transformation [1].Integrated energy systems (IES) are widely considered an important way to improve social energy efficiency and promote large-scale new energy development [2][3][4]."Integrated energy system" refers to a regional, integrated energy production, supply, and marketing system based on advanced physical information technology in a certain region; it contains two or more types of energy, including electricity, gas, heat, and cold, and multiple energy sources coupled with each other, and achieves high efficiency, clean, low-carbon use of energy through the coordinated optimization, multienergy complementation, and gradient utilization of the energy production, transmission, conversion, storage, and use of the energy [5,6].Under the conditions of an integrated energy system with multi-energy coupling, traditional electric power demand response (DR) has been expanded to integrated demand response (IDR) [7].
Demand response (DR) is a smart grid technology that targets the demand-side active response support capability in the context of insufficient grid peaking capacity.It helps the grid to cut peaks and fill valleys, ensuring stable operation through optimal allocation of load-side resources [8].Currently, China does not have a fully developed electricity market.The method of integrated demand response, which involves guiding load participation in the scheduling plan, can be divided into two types: tariff-based and incentive-based [9].The tariff type helps users optimize their energy use by setting appropriate power price signals, such as time of use (TOU), real-time pricing (RTP), and critical peak pricing (CPP).The incentive type is achieved through the signing of a demand response contract between the power company and the user, as well as through measures such as direct financial compensation and penalties for breach of contract.Through the demand response contract, the power company incentivizes users to support the power grid by providing direct economic compensation and default penalties.This approach aims to guide users towards supporting the power grid.
Traditional demand response only considers the reduction or transfer of electric energy on a limited time scale, which can result in economic losses and discomfort for users in their daily lives.As a result, it has significant limitations and cannot fully explore or mobilize the potential of load-side response [7].Integrated demand response (IDR) is the expansion of demand response in the context of integrated energy development.It makes full use of the multi-energy coupling and complementary characteristics of the integrated electricity-gas-heat-cooling energy system.IDR considers not only the reduction or transfer of electricity-gas-heat-cooling loads in the time scale but also the transformation of loads in the scale of the type of energy consumed, further improving the demand-side flexibility [10].The authors of [11] address the uncertainty of IDR and propose a twolayer planning model that considers the optimal allocation of equipment and operation strategy to enhance the effect of source-load interaction.The authors of [12] establish an IDR integrated energy efficiency model and propose a multi-objective integrated energy system scheduling method that considers integrated energy efficiency.The authors of [13] propose a distributed IDR strategy that considers the thermal delay effect for an industrial park with CHP units and heat accumulators.IDR can enhance demand-side flexibility by switching the type of energy use while ensuring customer comfort and economy [14,15].
The micro-energy grid is the smallest component of the integrated energy system and is a crucial element [16].Currently, there are more studies on the participation of micro-energy grids in integrated demand response [17][18][19].The authors of [20] present a price-based integrated demand response (P-IDR) method for micro-energy systems, enhancing optimization potential through energy balance equation modelling, multi-energy substitution effects consideration, and application of price elasticity theory and discrete choice theory; they indicate improved renewable energy accommodation, reduced peak/offpeak energy load differences, and cost savings through energy substitution characteristics consideration.The authors of [21] present an IDR model to enhance operational flexibility while considering customer satisfaction and investigate the coordinated scheduling and optimal operation strategy of a coupled heat-power-gas microgrid.Case studies demonstrate the effectiveness of the proposed strategies across various energy dispatching scenarios.However, these studies are still dominated by electric power demand response, with less consideration given to the integrated energy system's multi-energy coupling and complementation.
According to the way in which customer loads participate in integrated demand response, they can be divided into three groups: curtailable loads, convertible loads, and transferable loads.These are collectively referred to as dispatchable loads [14].Curtailable load refers to loads that can reduce or increase part or all of their energy consumption during peak hours to participate in demand response; of these, interruptible load (IL) is widely adopted as a key peaking resource both domestically and internationally [22,23].Transferable loads are loads that can be shifted in time to regulate the behaviour of energy use, where the total energy consumption remains approximately constant, including electric vehicles, energy storage devices, and certain processes [24].The authors of [25] present a novel robust decentralized charging strategy for large-scale EV fleets, which relies on uncertainty set-based robust optimization and employs an extended Jacobi proximal alternating direction method of multipliers algorithm; they demonstrate cost reduction and solution robustness in the simulations, benefiting EV users and the power grid.Convertible loads optimize energy use by altering the energy supply structure of the loads, without requiring load reduction or regulation in real time [26].However, studies often consider simple and homogeneous energy devices, storage types, and dispatchable loads, which fail to fully exploit the integrated demand response potential of micro-energy grids.
In this regard, this paper proposes a comprehensive demand response strategy for micro-energy networks, taking into account dispatchable loads; in the context of the increasing coupling of electricity-gas-heat-cold multiple energy sources and the rapid growth of new energy penetration, we take the micro-energy network as a typical research object for comprehensive demand response within the integrated energy system.Thus, we present an in-depth study of the comprehensive demand response strategy for the key issues.Section 2 takes into account the graded and gradient utilization of thermal energy, constructing a typical micro-energy grid structure, including electricity, gas, medium-grade heat, low-grade heat, and cold energy.Here, we establish a more comprehensive energy equipment model, and give a refined scheduling model for curtailable loads, transferrable loads, and convertible loads; additionally, an energy balance model for a micro-energy grid that takes into account a variety of storage and dispatchable loads on the basis of the model is established.Section 3 takes the operation economy of a micro-energy grid as the optimization goal, and proposes a micro-energy grid participation price model that takes into account the dispatchable loads.In Section 3, the optimal operation method of microenergy grid participating in price-type integrated demand response and incentive-type integrated demand response, taking into account dispatchable loads, is proposed.Here, we have the optimization objective of a micro-energy grid operation economy.In Section 4, an example analysis is carried out.In Section 5, a summary of the research content and results of this paper is presented.

Energy Supply Structure of Micro-Energy Grid
A micro-energy grid is the smallest cell of the integrated energy system, and is an important part of the regional integrated energy system.A micro-energy grid can be involved in community, factory, building, and other user energy systems.It is a micro integrated energy interconnection system centred on the electric power system, that is closely coupled with the electric power system, the natural gas system, the heating system, and the cooling system.It horizontally realizes the multi-energy complementarity of electricity, gas, cooling, heat, and vertically realizes the high degree of coordination of the links of "source, grid, load, and storage".The micro-energy interconnection system is a highly coordinated system of "source, grid, load and storage" in the vertical direction [27].The typical structure of a micro-energy grid is shown in Figure 1 below, which includes a variety of energy loads, abundant energy equipment, and the coupling and conversion of multiple energy sources such as cold, heat, electricity, and gas.
Figure 1 shows that the micro-energy grid is connected to the grid, gas network, and heat network of the regional integrated energy system.The micro-energy grid has the ability to purchase electricity, natural gas, and medium-grade heat (170 • C~550 • C) steam from the regional integrated energy system.The micro-energy grid comprises interconnected electricity, medium-grade heat, low-grade heat (room temperature ~170 degrees Celsius), and cold multiple energy supply grids, along with their related energy loads.The main sources of electrical energy are the external grid, the gas turbines, and the photovoltaics, including battery storage.The primary sources of medium-grade heat are external heat networks and waste heat boilers.The primary sources of low-grade heat include waste heat boilers, heat exchangers, electric boilers, and heat storage.Cooling energy is mainly derived from absorption chillers, electric refrigeration air conditioners, and ice storage equipment.Figure 1 shows that the micro-energy grid is connected to the grid, gas network, and heat network of the regional integrated energy system.The micro-energy grid has the ability to purchase electricity, natural gas, and medium-grade heat (170 °C~550 °C) steam from the regional integrated energy system.The micro-energy grid comprises interconnected electricity, medium-grade heat, low-grade heat (room temperature ~170 degrees Celsius), and cold multiple energy supply grids, along with their related energy loads.The main sources of electrical energy are the external grid, the gas turbines, and the photovoltaics, including battery storage.The primary sources of medium-grade heat are external heat networks and waste heat boilers.The primary sources of low-grade heat include waste heat boilers, heat exchangers, electric boilers, and heat storage.Cooling energy is mainly derived from absorption chillers, electric refrigeration air conditioners, and ice storage equipment.

Gas Turbine
Gas turbines (GTs) produce electricity by consuming natural gas, and their energy conversion is modelled as follows:

GT
GT.e gas GT ( )= ( ) where PGT(t) and GGT(t) are the output electric power and input natural gas flow rate of the gas turbine at time period t, respectively; λgas is the calorific value of natural gas; and ηGT.e is the power generation efficiency of the gas turbine.

Waste Heat Boiler
Waste heat boilers (WHBs) can recover the high-temperature hot steam produced by gas turbines while generating electricity and converting it into medium-grade and lowgrade hot steam to supply the related energy conversion equipment and heat load.Its energy conversion model is as follows: where HM_WHB(t) and HL_WHB(t) are the output medium-grade and low-grade heat power of the waste heat boiler at time period t, respectively; ηM_WHB and ηL_WHB are the mediumgrade heat and low-grade heat recovery efficiency of the waste heat boiler, respectively.

Gas Turbine
Gas turbines (GTs) produce electricity by consuming natural gas, and their energy conversion is modelled as follows: where P GT (t) and G GT (t) are the output electric power and input natural gas flow rate of the gas turbine at time period t, respectively; λ gas is the calorific value of natural gas; and η GT.e is the power generation efficiency of the gas turbine.

Waste Heat Boiler
Waste heat boilers (WHBs) can recover the high-temperature hot steam produced by gas turbines while generating electricity and converting it into medium-grade and low-grade hot steam to supply the related energy conversion equipment and heat load.Its energy conversion model is as follows: where H M_WHB (t) and H L_WHB (t) are the output medium-grade and low-grade heat power of the waste heat boiler at time period t, respectively; η M_WHB and η L_WHB are the mediumgrade heat and low-grade heat recovery efficiency of the waste heat boiler, respectively.

Heat Exchanger
Heat exchangers (HEs) can transfer heat from fluid with high temperature to fluid with low temperature, and their energy conversion model is: where H M_HE (t) and H L_HE (t) are the input medium-grade heat power and output lowgrade heat power of the heat exchanger, respectively; η HE is the low-grade heat recovery efficiency of the heat exchanger.
where H L_EB (t) is the output low-grade heat power of electric boiler at time period t; P EB (t) is the consumed electric power; η EB is the heat production efficiency of electric boiler.

Electric Refrigeration Air Conditioner
Electrically refrigeration air conditioners (ACs) are commonly used room-temperatureregulating equipment with the following energy conversion model: where H C_AC (t) is the cold power output of electric refrigeration air conditioner in time period t; P AC (t) is the consumed electric power; η AC is the refrigeration efficiency of electric refrigeration air conditioner.

Absorption Chiller
Absorption chillers (ACs) can absorb the waste heat generated in the process of energy production and conversion into cold energy output, and their energy conversion model is as follows: where H C_AB (t) is the cold power output of the absorption chiller at time period t; H L_AB (t) is the low-grade heat power consumed; η AB is the conversion efficiency of the absorption chiller.

Battery Energy Storage System
Battery energy storage systems (BESSs) are common electrical energy storage devices; they commonly use lead-acid batteries, lithium batteries, etc., which have charging, discharging, and standby states.The energy relationship between the amount of power stored in the battery and its charging and discharging power is as follows: where S BS (t) is the amount of power stored in the battery storage at time t; P BS,c (t) and P BS,d (t) are the charging and discharging power of the battery at time t, respectively; σ BS is the self-depletion coefficient of the battery storage; η BS,c and η BS,d are the charging and discharging efficiencies of the battery storage, respectively; and ∆t is the length of the unit time period.

Ice Storage Air Conditioning
Ice storage air conditioners (ISACs) consist of an electric refrigeration unit and an ice storage tank, which can use low valley electricity to make water into ice and store it in the ice storage tank, and then release the cold energy by melting ice in the peak power period.ISACs have a variety of operating modes, and this paper considers that they can make ice while supplying cooling, and can release cold simultaneously by both chiller and ice melting.The energy conversion model of this ice storage air conditioner is as follows: P IS (t) = P FRG (t) + P TK (t) (11) Equation ( 8) indicates that the cold power H C_IS (t) provided by the ice storage air conditioner at moment t is the sum of the cold power H C_FRG (t) generated by the electric refrigeration unit and the cold power H C_TK (t) generated by melting ice in the ice storage tank.Equation (9) indicates that the cold power generated H C_FRG (t) needs to consume the electric power of P FRG (t); η FRG is the refrigeration efficiency of refrigeration unit in ice storage air conditioner.Equation (10) is the relationship between the ice storage volume in the ice storage tank and the ice-making power and ice-melting power; S TK (t) is the ice storage volume in the ice storage tank at time t; P TK (t) is the electric power consumed by the ice making in the ice storage tank; η TK,c is the energy efficiency ratio of the ice making in the ice storage tank; η TK,d is the ice melting efficiency in the ice storage tank; and σ TK is the self-depletion coefficient of the ice storage tank.In Equation (11), P IS (t) is the total electric power consumed by the ice storage air conditioner at the moment t.

Phase Change Heat Storage
Phase change heat storage (PCHS) is a device that stores and releases heat through phase change materials.In this paper, phase change heat storage is mainly used to store low-grade heat, and the relationship between the stored heat, the heat storage power, and the heat release power is as follows: where S HS (t) is the amount of heat stored in the heat storage in time period t; H L_HS,c (t) and H L_HS,d (t) are the heat storage and exothermic power of the heat storage in time period t, respectively; η HS,c and η HS,d are the heat storage and exothermic efficiencies of the heat storage, respectively; and σ HS is the self-loss coefficient of the heat storage.

Dispatchable Load Model
According to the way in which customer loads participate in the integrated demand response, they can be divided into three categories: curtailable loads, transferable loads, and convertible loads.These are collectively referred to as dispatchable loads [3].

Curtailable Loads
Curtailable loads are loads that can reduce or increase part or all of their energy consumption during peak hours to participate in demand response.In this paper, curtailable loads include electric loads, medium-grade heat loads, low-grade heat loads, and cooling loads, and their unified models are as follows: where i indicates the i-th curtailable load; n cut is the number of curtailable loads; L cut (t) is the actual power of the curtailable load in time period t; L max cut,i (t) indicates the maximum power of the curtailable load i in time period t; u cut,i (t) is a 0/1 state parameter, which indicates that the curtailable load i has been curtailed when it is 1, and indicates that it is not curtailed when it is 0; the curtailment ratio µ cut,i (t) of the curtailable load i takes values in the range [0, 1] and indicates complete curtailment of load i during time period t when it is 1.
The operating constraints of the curtailable load include the total curtailment duration constraint (Equation ( 14)), and the single curtailment duration limit (Equation ( 15)): where N max cut,i is the maximum number of cuttable time slots for cuttable load i in a dispatch cycle T; τ is the start time of cut; T min cut,i and T max cut,i are the single minimum and maximum duration of cut for cuttable load i, respectively.

Transferable Load
Transferable loads are loads that can be shifted in time to regulate the behaviour of energy use, where the total energy consumption remains approximately constant, including electric vehicles, energy storage devices, and certain processes.Generally transferable loads (except for levelling loads) do not have continuity and timing requirements, and only need to satisfy the load power limitation, transfer time constraints, and constant energy consumption constraints, which are modelled as follows: where L trf,i,o (t) and L trf,i (t) are the power magnitude of the transferable load i before and after the transfer at time t, respectively; u trf,i is the 0/1 state parameter of the transferable load, which indicates that the load is being transferred when the value is 1, and that it is not transferred when the value is 0; L min trf,i and L max trf,i are the upper and lower limits of the operating power of the transferable load, respectively; [t min trf,i , t max trf,i ] are the upper and lower limits of the operating power of the transferable load, respectively.
One type of transferable load is the shifter load, which includes household appliances such as rice cookers and washing machines.The shifter load is a special case of the transferable load, which is constrained by the production process and can only be shifted as a whole on a time scale, and is modelled as follows: where L shift,i (t) is the power size of the shifter load i at time t; u shift,i is the 0/1 state parameter of the shifter load, when its value is 1 it means that in the load is shifted, and when its value is 0 it means that it has not been shifted; and v shift,i is the shifting time of the shifter load i.Assuming that the allowable running time period of the shifter load i is [t min shift,i ,t max shift,i ], the start running time of the shifter load is τ i , and the duration running time is ω i , the shift time v shift,i needs to satisfy

Convertible Load
Convertible load is the unique dispatchable load in the integrated energy system due to the complementary nature of multi-energy coupling.Convertible loads optimize energy use by altering the energy supply structure of the loads, without requiring load reduction or regulation in real time; thus, they have no impact on user comfort.Convertible loads are loads with multiple energy supply modes in the integrated energy system that can be independently switched by the user to choose the energy supply mode.This term does not refer to a specific device but rather to a certain kind of energy demand.In the micro-energy grid system shown in Figure 1, without considering energy storage, the user's cold demand can be satisfied by two types of energy conversion equipment, namely electric refrigeration and absorption refrigeration, and its energy conversion model is as follows: where L c (t) is the size of the cold load at time t; P change (t) and H L,change (t) are the size of the consumed electric power and low-grade thermal power at time t, respectively; η e2c and η l2c represent the efficiency of the electric cooling conversion and the efficiency of the heat and cooling conversion, respectively; α e2c and α l2c represent the ratio of the energy consumption of the electric cooling conversion and the heat and cooling conversion, respectively; the energy conversion strategy within the micro-energy grid can be changed to obtain the energy conversion strategy, without considering the energy storage.The energy conversion strategy within the energy network is used to obtain different energy supply structures.

Objective Function
The main appeal of micro-energy grids to participate in integrated demand response is to reduce operating costs and obtain greater economic benefits.The objective function of micro-energy grid participation in tariff-based integrated demand response is as follows: where C e is the daily cost of electricity purchased in the micro-energy grid; C g is the daily cost of gas purchased in the micro-energy grid; C h is the daily cost of heat purchased in the micro-energy grid; C GT is the cost of gas turbine power generation in the micro-energy grid; C cut is the cost of curtailable load response in the micro-energy grid; C trf is the cost of shifting load response in the micro-energy grid; C shift is the cost of levelling load response in the micro-energy grid, measured in terms of the time of levelling and the amount of levelling of the load volume to measure the levelisable load response cost; C om is the equivalent daily operation and maintenance cost of energy equipment in the micro-energy grid; C se is the start-stop cost of energy equipment in the micro-energy grid; C depre is the equivalent daily charging and discharging depreciation cost of battery storage; c e (t) is the time-of-day tariff; P buy (t) is the micro-energy grid power purchased in time period t; c g is the price of natural gas; G buy (t) is the t the gas flow rate purchased by the micro-energy grid in time period t; c h is the price of thermal steam; H buy (t) is the thermal power purchased by the micro-energy grid in time period t; n GT is the number of gas turbines in the microenergy grid; a GT,i , b GT,I , and c GT,i are the generation cost factors of common generating units; c cuta,I and c cutb,i are the cost factors of curtailable load i when it participates in demand response; c trf,i is the cost factor of the cost coefficient per unit shifted power when shifted load i participates in demand response; c shift1,i and c shift2,i are the cost coefficients when levelling load i participates in demand response; c om,i is the equivalent daily operation and maintenance cost per unit of output power of equipment i; P equip,i (t) is the output power of equipment i in time period t; n equip is the number of equipment devices; c se,i is the single start-stop cost of device i; u i (t) is the operation 0/1 state parameter of device i in time period t.When its value is 1, it means that the device is running, and when its value is 0, it means that the device stops running.P e BS,d , D e BS and S e BS are the rated discharge power, the rated depth of discharge, and the rated capacity, respectively, of the battery storage energy; p BS,d (t) and q BS,d (t) are the actual discharge power and actual discharge capacity of the battery storage energy in time period t; D BS (t) is the actual discharge depth of the battery storage energy in time period t; C BS_buy is the purchase cost of the battery storage energy; δ 1 and δ 2 are the coefficients between the relationship between the number of battery cyclic charging and discharging times and the discharge depth.
When the micro-energy grid responds to the higher-level peak regulation demand, that is, participates in the incentivized comprehensive demand response, the net cost of the optimized operation of the micro-energy grid needs to be subtracted from the peak shaving subsidy income.This is obtained by participating in the incentivized demand response on the basis of Equation (22), and the objective function model is as follows: where the calculation formula of the peak shaving compensation I IDR obtained by the micro-energy grid is as follows: where c IDR is the compensation price for peak shaving, t IDR,s and t IDR,e are the start and end times of peak shaving, and P buy,o (t) and P buy (t) are the power purchase power before and after the participation of the micro-energy grid in the incentive comprehensive demand response, respectively.

Constraints
( Medium-grade thermal power balance constraint: Low-grade thermal power balance constraint: Cold power balance constraint: where L e (t), L M (t), L L (t), and L C (t) are the predicted values of the non-adjustable electric load, medium-grade heat load, low-grade heat load, and cooling load in the micro-energy grid in the t period, respectively.P buy (t) is the input power of the micro-energy gateway interface, P GT (t) is the output power of the gas turbine, and P PV (t) is the value of the photovoltaic output in the micro-energy grid.P change1 (t), P change2 (t), H M,change (t), and H L,change (t) are the input power of the energy conversion equipment corresponding to the convertible loads, respectively; H M_buy (t) is the thermal power purchased from outside the microenergy grid in time period t; η M2L , η e2C , and η M2C are the energy conversion efficiency of electric heating, heat exchange, electric refrigeration, and absorption refrigeration of the convertible equipment corresponding to the convertible load in the micro-energy grid, respectively [0, 1].
(2) Equipment operation constraints.The operational constraints of gas turbines include the upper and lower limits of output power (Equation (39)) and the maximum climbing rate constraints (Equation (40)).

P min
GT ≤ P GT (t) ≤ P max GT (39) where P min GT and P max GT are the upper and lower limits of the output power of the gas turbine, respectively, and r max GT is the maximum ramp rate of the gas turbine.The operational constraints of battery energy storage include SOC constraints (Equation (41)), power constraints (Equation (42)), operating state constraints (Equation (43)), and ramp-up constraints (Equation ( 44 where SOC(t) is the ratio of the battery storage power to the rated capacity at time t; SOC max and SOC min are the upper and lower limits of the battery running SOC; P max BS,c and P min BS,c are the upper and lower limits of the battery charging power; P max BS,d and P min BS,d are the upper and lower limits of battery charging power, respectively; u c (t) and u d (t) are the state parameters of charge and discharge 0/1, respectively.When u c = 1, it means that the battery is charged; when u d = 1, the battery is discharged; r max BS is the maximum ramp-up rate of battery energy storage.
The operational constraints of ice storage air conditioning include ice storage capacity constraints (Equation (45)), power constraints (Equation ( 46 where S max TK and S min TK are the upper and lower limits of the ice storage capacity of the ice storage tank, respectively; P max IS is the maximum power consumption of the ice storage air conditioner; H max C_TK is the maximum cooling power generated by ice melting in the ice tank; u e1 (t), u e2 (t), and u s (t) are 0/1 state parameters of ice storage air conditioning refrigeration mechanism cooling, ice storage tank ice making, and ice melting in ice storage tank; when it is equal to 1, it is the running state, and when it is equal to 0, it is the stop state; r max TK is the maximum climb rate of the ice tank.

S min
HS ≤ S HS (t) ≤ S max HS (49) where S min HS and S max HS are the lower and upper limits of heat storage capacity of the accumulator, respectively; H max L_HS,c and H min L_HS,c are the upper and lower limits of heat storage power, respectively; H min L_HS,d and H max L_HS,d are the upper and lower limits of the exothermic power, respectively; u hc (t) and u hd (t) are the state parameters of 0/1 for heat storage and exothermic heat storage, respectively; when u hc = 1, this represents the heat storage of the accumulator; when u hd = 1, the accumulator is exothermic.Equation (52) indicates that the accumulator cannot store and operate exothermically at the same time; r max HS is the maximum creep rate of the accumulator.
In ( 55), the confidence level c characterizes the satisfaction of interval constraints, that is, the probability of interval constraint (54).
(5) Constraints on the transmission power of the gateway contact line.
0 ≤ P buy (t) ≤ P max buy (56) where P max buy is the upper limit of the transmission power of the micro-energy gateway port.( 6) Incentive-based comprehensive demand response peak shaving power constraints.
where p IDR (t) is the peak shaving power requirement issued by the superior.

Model Solving
This paper presents an MINLP problem, where the model includes control variables that consist of both continuous variables-P PV (t) and H min L_HS,d ; L cut (t) has an upper limit of L max cut,i (t); L trf,i (t) has an upper and lower limit of L min trf,i and L max trf,i .The discrete variables u c (t) and u d (t) are the state parameters of charge and discharge 0/1, respectively, which represent battery charging when u c = 1, and battery discharging when u d = 1; u e1 (t), u e2 (t), and u s (t) represent the 0/1 state parameters for the ice storage air conditioner chiller refrigeration, ice production in the ice storage tank, and ice melting in the ice storage tank, respectively.A value of 1 indicates a running state, while a value of 0 indicates a stopping state.u hc (t) and u hd (t) are parameters that represent heat storage and exothermic 0/1 states of the phase change heat storage, respectively.When u hc (t) = 1, it represents heat storage in the accumulator; when u hd (t) = 1, the accumulator is exothermic; u cut,i (t) is a 0/1 state parameter, which when it is 1, this indicates that the curtailable loadiis curtailed; when it is 0, this indicates that it is not curtailed.u trf,i (t) is a 0/1 state parameter.When it is 1, this indicates that the curtailable load is curtailed; when it is 0, this indicates that it is not curtailed.
MATLAB was used to build the model for this paper.The CPLEX solver was called in MATLAB to solve the problem and obtain the distribution of equipment power in each time period of the scheduling study.This was performed using the PV prediction value, load power, equipment parameters, and other relevant data.

Case Overview
An industrial park micro-energy network is used as an example to validate the integrated demand response strategy for micro-energy networks, taking into account the dispatchable loads proposed in this paper.The optimization calculations were performed using MATLAB R2019b, CPLEX 12.8 optimizer running on a 64 bit windows 10 operating system, with an Intel Core i7-6500U CPU (manufactured by Intel, headquartered in Santa Clara, CA, USA) model running at 2.50 GHz and 4 GB of RAM.
The energy supply structure of the micro-energy grid is shown in Figure 1, and it considers the dispatchable load.Table 1 displays the energy equipment featured in the micro-energy network along with their respective parameters.Table 2 presents the types of dispatchable loads included in the micro-energy network, along with their performance and price parameters.Figure 2 illustrates the photovoltaic output, as well as the electricity, medium-grade heat, low-grade heat, and cooling loads of the micro-energy grid.

Price-Based Comprehensive Demand Response Results
The CPLEX12.8 optimization software was used to solve the results of the time-ofuse electricity price optimization operation of the micro-energy grid s participation in the price-based comprehensive demand response.Figure 3 shows the optimal operation results of the power supply system of the micro-energy grid, where time 1 corresponds to 0:00.It can be seen from the figure that photovoltaic power generation is 100% consumed, and the gas turbine is in the trough electricity price period from 22:00 to 8:00 the next day.From 11:00 to 13:00, the power generation is 691.75 kW, and the full power generation is 1000 kW during peak and peak hours.This is because the power generation cost of gas turbines is higher than the time-of-use electricity price during valley hours.Additionally, battery energy storage is charged at 3:00-8:00 and 11:00-13:00 during the trough hours, and at 22:00, during peak electricity prices, and discharged at peak 8:00-11:00, with 13:00-17:00 , through the electricity price peak-to-valley arbitrage to further reduce the operating cost of the micro-energy grid in the peak.At the peak hours of electricity prices, the overall power purchase of the micro-energy grid decreases.Figure 4 displays the optimal operational results of the micro-energy grid heating system.The medium-grade and low-grade thermal power output by the waste heat boiler primarily depend on the electrical power output by the gas turbine.The heat accumulator Additionally, Table 3 shows the adjusted time-of-use electricity price for Zhejiang Province (China) starting from 15 October 2021.Furthermore, the cost of medium-grade hot steam is CNY 0.465 per kilowatt-hour, while the cost of natural gas is CNY 2.5 per cubic meter, which is equivalent to CNY 0.349 per kilowatt-hour.The subsidy for grid demand response peak shaving is a power subsidy priced at CNY 4 per kilowatt-hour; the valley filling compensation is a capacity subsidy priced at CNY 5 per kilowatt-day.

Price-Based Comprehensive Demand Response Results
The CPLEX12.8 optimization software was used to solve the results of the time-of-use electricity price optimization operation of the micro-energy grid's participation in the price-based comprehensive demand response.Figure 3 shows the optimal operation results of the power supply system of the micro-energy grid, where time 1 corresponds to 0:00.It can be seen from the figure that photovoltaic power generation is 100% consumed, and the gas turbine is in the trough electricity price period from 22:00 to 8:00 the next day.From 11:00 to 13:00, the power generation is 691.75 kW, and the full power generation is 1000 kW during peak and peak hours.This is because the power generation cost of gas turbines is higher than the time-of-use electricity price during valley hours.Additionally, battery energy storage is charged at 3:00-8:00 and 11:00-13:00 during the trough hours, and at 22:00, during peak electricity prices, and discharged at peak 8:00-11:00, with 13:00-17:00, through the electricity price peak-to-valley arbitrage to further reduce the operating cost of the micro-energy grid in the peak.At the peak hours of electricity prices, the overall power purchase of the micro-energy grid decreases.
battery energy storage is charged at 3:00-8:00 and 11:00-13:00 during the trough hours, and at 22:00, during peak electricity prices, and discharged at peak 8:00-11:00, with 13:00-17:00 , through the electricity price peak-to-valley arbitrage to further reduce the operating cost of the micro-energy grid in the peak.At the peak hours of electricity prices, the overall power purchase of the micro-energy grid decreases.Figure 4 displays the optimal operational results of the micro-energy grid heating system.The medium-grade and low-grade thermal power output by the waste heat boiler primarily depend on the electrical power output by the gas turbine.The heat accumulator Figure 4 displays the optimal operational results of the micro-energy grid heating system.The medium-grade and low-grade thermal power output by the waste heat boiler primarily depend on the electrical power output by the gas turbine.The heat accumulator stores heat during the trough period of electricity price, releasing it from 8:00 to 11:00 during peak hours and from 13:00 to 14:00.The electric boiler operates during peak hours of electricity price from 8:00 to 11:00 and from 13:00 to 21:00.Heat exchangers do not produce low-grade heat.Instead, they only produce heat during peak and off-peak hours of electricity prices.This is because electric heating is more expensive during peak hours than the cost of using medium-grade heat exchange.Conversely, during off-peak hours, electric heating is cheaper than the cost of using medium-grade heat exchange.By using this complementary approach, the operating cost of the micro-energy grid can be further reduced.stores heat during the trough period of electricity price, releasing it from 8:00 to 11:00 during peak hours and from 13:00 to 14:00.The electric boiler operates during peak hours of electricity price from 8:00 to 11:00 and from 13:00 to 21:00.Heat exchangers do not produce low-grade heat.Instead, they only produce heat during peak and off-peak hours of electricity prices.This is because electric heating is more expensive during peak hours than the cost of using medium-grade heat exchange.Conversely, during off-peak hours, electric heating is cheaper than the cost of using medium-grade heat exchange.By using this complementary approach, the operating cost of the micro-energy grid can be further reduced.

Power(kW)
Electricity price (yuan/kW)  When comparing the results of the micro-energy grid with and without ice storage air conditioner participation, the daily operating cost was CNY 133,439.05 and CNY 133,452.61,respectively, after the price-based integrated demand response.Ice storage air conditioners generally have a lifespan of over 15 years.It has been determined that their lifespan is specifically 15 years, resulting in a time-sharing tariff arbitrage of CNY 74,241.Additionally, participating in peak demand response can reduce the peak by 130 kW, 30 times a year, resulting in a peak subsidy of CNY 234,000.The total estimated gain is CNY 308,241.Equipping the micro-energy grid with an ice storage air-conditioner costing CNY 200,000 can recover the investment cost and generate revenue.

Incentivized Comprehensive Demand Response Results
Considering the two types of peak shaving demand issued by the power grid: peak shaving of 1200 kW during the 15:00-16:00 period and peak shaving of 1500 kW during the 15:00-16:00 period, the comparison of the operation of the micro-energy grid before and after participating in the incentive demand response is shown in Figure 5 below.Figure 5 shows three bars for each item, representing the relevant equipment within the micro-energy network before participating in incentivized demand response, peak shaving of 1450 kW, and peak shaving of 1750 kW.It is evident that the electric refrigeration air conditioner generates significantly less cooling power after participating in the incentive-based integrated demand response.Additionally, peak shaving of 1750 kW results in a greater reduction in output power compared to peak shaving of 1450 kW, which effectively reduces power consumption.The decrease in cooling power of the electric refrigeration air conditioner is compensated for by the ice storage air conditioner and the absorption chiller.These systems melt and cool the ice at their maximum output power.The absorption chiller starts from its original shutdown cooling at full output power using low-grade heat.As a result, the purchased heat demand of the micro-energy network increases, the heat converted by the heat exchanger increases, and the heat accumulators release heat at maximum power.In regard to dispatchable loads, the cost of peak shaving for the electrically transferable load is CNY 1.38, with the lowest response cost.During the 15:00-16:00 period, 450 kW of electric load is transferred to the 0:00-1:00 tariff trough period, and the thermally levelling load does not contribute to peak shaving and does not respond to it.The electrically curtailable load is reduced by 726 kW during peak shaving of 1450 kW and by 900 kW during peak shaving of 1750 kW.The cold curtailable load has a higher response cost than the electric curtailable load and exceeds the peaking subsidy, so it is not curtailed at a peak shaving of 1450 kW.However, it is curtailed by 480 kW at a peak shaving of 1750 kW.
Based on the analysis above, it can be concluded that the micro-energy network described in this paper effectively utilizes its demand response potential through various energy storage methods, dispatchable loads, and multi-energy complementarity of electricity, heat, and cold.Additionally, it can reduce energy costs by participating in incentive-based demand response programs.

Conclusions
In this paper, we propose a comprehensive demand response optimization strategy for micro-energy grids considering the dispatchable load.Considering the cascade utili- Figure 5 shows three bars for each item, representing the relevant equipment within the micro-energy network before participating in incentivized demand response, peak shaving of 1450 kW, and peak shaving of 1750 kW.It is evident that the electric refrigeration air conditioner generates significantly less cooling power after participating in the incentivebased integrated demand response.Additionally, peak shaving of 1750 kW results in a greater reduction in output power compared to peak shaving of 1450 kW, which effectively reduces power consumption.The decrease in cooling power of the electric refrigeration air conditioner is compensated for by the ice storage air conditioner and the absorption chiller.These systems melt and cool the ice at their maximum output power.The absorption chiller starts from its original shutdown cooling at full output power using low-grade heat.As a result, the purchased heat demand of the micro-energy network increases, the heat converted by the heat exchanger increases, and the heat accumulators release heat at maximum power.In regard to dispatchable loads, the cost of peak shaving for the electrically transferable load is CNY 1.38, with the lowest response cost.During the 15:00-16:00 period, 450 kW of electric load is transferred to the 0:00-1:00 tariff trough period, and the thermally levelling load does not contribute to peak shaving and does not respond to it.The electrically curtailable load is reduced by 726 kW during peak shaving of 1450 kW and by 900 kW during peak shaving of 1750 kW.The cold curtailable load has a higher response cost than the electric curtailable load and exceeds the peaking subsidy, so it is not curtailed at a peak shaving of 1450 kW.However, it is curtailed by 480 kW at a peak shaving of 1750 kW.
Based on the analysis above, it can be concluded that the micro-energy network described in this paper effectively utilizes its demand response potential through various energy storage methods, dispatchable loads, and multi-energy complementarity of electricity, heat, and cold.Additionally, it can reduce energy costs by participating in incentive-based demand response programs.

Conclusions
In this paper, we propose a comprehensive demand response optimization strategy for micro-energy grids considering the dispatchable load.Considering the cascade utilization of thermal energy, a typical micro-energy grid structure including electricity, gas, medium-grade heat, low-grade heat, and cold energy is constructed.A comprehensive energy equipment model is established, including a variety of energy storage models: battery energy storage, phase change heat storage, ice-cold storage device, and a refined scheduling model that can reduce load and transfer load, and a convert load is given.On this basis, taking the operation economy of the micro-energy grid as the optimization goal, considering the constraints of new energy consumption rate and energy supply reliability, a comprehensive demand response strategy of the micro-energy grid participating in the electricity price type and incentive type is proposed.
Through the analysis of the examples, the following conclusions can be obtained: (1) The price-based comprehensive demand response strategy of micro-energy grid proposed in this paper can optimize the energy consumption structure of a microenergy grid under the guidance of time-of-use electricity price, reduce operating costs, and ensure the consumption of new energy.For the power grid, the peak-to-valley difference can be reduced to achieve a win-win situation.(2) The incentive-based comprehensive demand response strategy of the micro-energy grid proposed in this paper enables the micro-energy grid to minimize its own costs while responding to the peak regulation demand of the upper level and obtain the peak shaving subsidy to further reduce its own operating cost.The proposed strategy fully excavates the demand response potential of the micro-energy grid through dispatchable load; the complementarity of electricity, heat, and cold realizes the comprehensive coordination and optimization of source-grid-load-storage, provides greater peak regulation capacity, and obtains greater benefits when participating in the incentive-based comprehensive demand response.Thus, it provides theoretical reference and support for the majority of multi-energy users to participate in the comprehensive demand response.(3) The method proposed in this paper can estimate whether the investment cost can be recovered by the micro-energy grid during the life cycle of the phase change heat accumulator and ice-cold storage device; additionally, it can provide a reference for the energy storage planning of the micro-energy grid.For the case of the micro-energy grid in this example, compared with the case of only electrochemical energy storage, the economic benefit and lower net cost of the micro-energy grid participating in the comprehensive demand response with the configuration of ice storage device are greater.
This paper deeply explores the integrated demand response potential of micro-energy networks and proposes a method with good engineering practicality.However, the paper acknowledges that the method still has deficiencies.The paper presents an MINLP problem model that is challenging to solve.It can be applied to solve small-scale micro-energy network systems.However, as the system scale expands and equipment increases, the algorithm's complexity, computational resources, and required time also increase significantly.Therefore, when calculating demand response solutions for larger energy networks, we propose that one tries to consider the micro-energy network as a whole and establish a centralized-distributed mechanism.This will allow for joint optimization of the regional integrated energy system through participation with multiple micro-energy networks in incentive-based demand response.
Author Contributions: X.Z.conceived and designed the entire review and wrote the paper.H.W., M.Z., M.D. and S.D. reviewed and edited the manuscript.All authors read and approved the manuscript.All authors have read and agreed to the published version of the manuscript.

Figure 2 .
Figure 2. Micro-energy grid PV output and load forecast results before the day.

Figure 3 .
Figure 3.The optimization operation results of power supply system.

Figure 2 .
Figure 2. Micro-energy grid PV output and load forecast results before the day.

Figure 3 .
Figure 3.The optimization operation results of power supply system.

Figure 3 .
Figure 3.The optimization operation results of power supply system.

Figure 4 .
Figure 4.The optimization operation results of heating system.

Figure 4 .
Figure 4.The optimization operation results of heating system.

Energies 2024 , 20 Figure 5 .
Figure 5.Comparison chart of the operation of the micro-energy grid before and after participating in the incentivized comprehensive demand response.

Figure 5 .
Figure 5.Comparison chart of the operation of the micro-energy grid before and after participating in the incentivized comprehensive demand response.
, P buy (t), P GT (t), P BS,c (t), P BS,d (t), SOC(t), H M_buy (t).H M_WHB (t), H L_WHB (t), H M_HE (t), H L_HE (t), H L_EB (t), P AC (t), H C_AB (t), P IS (t), S TK (t), H C_TK (t), S HS (t), H L_HS,c (t), H L_HS,d (t), L cut (t), and L trf,i (t)-and discrete variables: u c (t), u d (t), u e1 (t), u e2 (t), u s (t), u hc (t), u hd (t), u cut,i (t), and u trf,i (t).The continuous variables P PV (t) have an upper limit of P PV * (t); and P buy (t) have an upper limit of P max buy ; P GT (t) has an upper and lower limit of P max GT and P min GT ; P BS,c (t) has an upper and lower limit of P max P BS,d (t) has an upper and lower limit of P max BS,d and P min BS,d ; SOC(t) has an upper and lower limit of SOC max and SOC min ; H M_buy (t) has an upper limit of H max M_buy ; H M_WHB (t) has an upper limit of H max M_WHB ; H L_WHB (t) has an upper and lower limit of H max L_WHB ; and H L_HE (t) has an upper limit of H max L_HE ; H L_EB (t) has an upper limit of H max L_EB ; P AC (t) has an upper limit of P max AC ; H C_AB (t) has an upper limit of H max L_AB ; P IS (t) has an upper limit of P max IS ; S TK (t) has an upper and lower limit of S max TK and S min TK ; H C_TK (t) has an upper limit of H max C_TK ; S HS (t) has an upper and lower limit of S min HS and S max HS ; H L_HS,c (t) has an upper and lower limit of H max L_HS,c and H min L_HS,c ; H L_HS,d (t) has an upper and lower limit of H max L_HS,d

Table 1 .
The micro-energy grid contains the device and its parameters.

Table 2 .
Types of dispatchable loads of micro-energy grid and their performance and price parameters.