*Article* **Economic Dispatch of Multi-Energy System Considering Load Replaceability**

**Tao Zheng 1,2,3,\*, Zemei Dai 1,2,3, Jiahao Yao 4, Yufeng Yang 1,2,3 and Jing Cao 1,2,3**


Received: 28 July 2019; Accepted: 23 August 2019; Published: 28 August 2019

**Abstract:** By integrating gas, electricity, and cooling and heat networks, multi-energy system (MES) breaks the bondage of isolated planning and operation of independent energy systems. Appropriate scheduling of MES is critical to the operational economy, and it is essential to design scheduling strategies to achieve maximum economic benefits. In addition to the emergence of energy conversion systems, the other main novelty of MES is the multivariate of load, which offers a great optimization potential by changing load replaceability (flexibly adjusting the composition of loads). In this paper, by designing load replaceability index (LRI) of composite load in MES, its interaction mechanism with scheduling optimum is systematically analyzed. Through case studies, it is proven that the optimum can be improved by elevating load replaceability.

**Keywords:** multi-energy system; economic dispatch; load replaceability; multi-energy conversion

#### **1. Introduction**

Energy flexibility is crucial to the construction of more environmental-friendly and efficient power generation (consumption) patterns. Unlike conventional single-form energy-based energy system, which is usually detached from other types of energy systems from both the planning and operational perspectives, multi-energy system (MES) comprises various forms of subsystems including electricity, heat, cooling and gas networks, contributing to the interactions and corresponding energy flexibility through energy converters. Therefore, MES technologies gradually become irreplaceably key components in building the modern energy systems.

Compared with single-from energy-based system, the most significant difference of MES is the energy conversion function achieved by energy conversion equipments, which serve as in-between interfaces and assist in the formation of coupling relations. Energy converters are foundations of transformation, fusion and decomposition of energy flows among different subsystems. Because of the possibility of convertibility of different energy forms, MES outperforms systems of pure energy form by complementing the advantages of each subsystem.

An important concept in MES is the energy hub (EH), which gives an matrix model of production, conversion, consumption and storage of different energy carriers. Each of the matrix elements represent the abstract connection and conversion coefficients of internal components, which can be embodied by various types of multi-generators, among which combined heat and power (CHP) or combined heat cooling and power (CHCP) plants are the most widely investigated cogeneration or trigeneration plants.

Various research studies have recently begun to show intensive concern for the design and operation of MES, either from the perspective of microscopic (CHP and CHCP plants) or macroscopic (EHs). The main idea of optimal operation of CHP (CHCP) plants is similar to that of power systems, both starts from the optimization model building and ends with model solving by specific optimization algorithms like mathematical or evolutionary programming methods. Like economic electrical power dispatch, the first priority of multi-generation is economy. Some researchers leverage the strong computing capabilities of meta-heuristics to solve economic dispatch of CHP plants [1,2]. Apart from deterministic optimization, uncertainties stemming from volatile renewable energy are gradually introduced by stochastic programming. For example, a novel chance constrained programming model is used to formulate economic dispatch with CHP plants and wind power simultaneously [3]. Li groundbreakingly investigates how to use the temperature dynamics of heating networks to enhance the utilization of wind power, and a CHP dispatch model considering operational economy and wind power integration is built and solved by iterative methods [4]. As for economic dispatch of CHCP, a nonlinear economic optimal operation model is developed and solved by lingo solver [5]. Backward dynamic programming technique is used to obtain the optimal operational set-points of combined heat, cooling and power systems [6].

Besides designing operation schedules solely in pursuit of the economic potential elevation, optimization of (CHP) CCHP plants considering multiple operational output performance begins to capture researchers' interests. Under most circumstances, economy, energy utilization efficiency and environmental preservation are the three main optimization goals. By employing Weighting method and fuzzy optimum selection theory, the integrated performances of operation of CCHP plants considering these objectives are evaluated [7]. Dynamic object method-particle swarm optimization method is used to achieve the goal of carbon emission reduction, energy efficiency elevation and system cost minimization [8]. Additionally, some new evaluation indices evaluating the average useful output and total heat transfer area are studied, correspondingly, a multi-objective optimization-based schedule is designed [9]. Modified Bacterial Foraging Optimization method is adopted to minimize the operational cost and pollutant emission [10].

Optimal system-level operation of MES is usually implemented under EH framework by designing the synergies among generation, consumption and storage components. An optimal operation strategy of MES containing CHP plants, boiler, battery and water tank is studied [11], and particle swarm optimization (PSO) algorithm is used to obtain optimal scheduling of each device. An optimal expansion planning of EH containing CHP plants and natural gas furnaces is designed to minimize investment cost [12]. Under some circumstances MES have some uncertainties involving demand fluctuation, cost dynamics and converter efficiency change. In this situation deterministic framework no longer adapts to uncertain environments, hence robust optimization approach is used to increase the robustness towards uncertain bounded parameters [13]. In [14], uncertainties of demand, market price and renewable energy are considered by stochastic programming, and optimal operation set-points are obtained to minimize operation costs. Instead of treating EH as an individual system, studies of multi-EH systems take into consideration the mutual influence and restriction. By analyzing competitive and cooperative relation between multi-EHs, game theoretic optimal scheduling is established through quantum particle swarm optimization method [15]. Moreover, resilience of MES is researched in last several years, and coordination between different subsystems in post-disaster repair is investigated to reduce load shedding and repair duration [16].

Most of the existing MES research focuses on optimal planning or operation by considering multi-energy generation or conversion on the generation side. Recently, some research studies begin to consider the characteristics of multi-energy loads, which play essential roles in MES planning & operation as well as other applications like load prediction and demand response programs. Some researchers investigate the influence of multi-energy loads upon MES planning, and the capacity of gas turbine under different thermoelectric ratios is studied [17]. In respect to load prediction, load prediction for integrated cooling and heat system is implemented by considering the intrinsic coupling among multi-energy loads using multivariate phase space reconstruction technique [18]. By considering the operating attributes of loads, a quick cooling and heating load prediction model

is constructed [19]. Deep-learning technique is adopted to predict various types of load in regional energy system integration [20]. Some researchers also study how to adjust the consumption pattern to better participate in the demand response program [21,22].

Conventionally, the load composition from each energy network is assumed to be fixed without replaceability, which is characterized by the feasibility of substitution of one form of load with another one. Load considering replaceability can be regarded as a composite load, it characterizes the multi-energy consumption on the load side and is essential to flexible and efficient energy consumption. Composite load can adjust the proportion of heterogeneous energy consumption based on the external excitation signal such as market price, thus achieving simultaneous operation cost reduction and the load and generation balance.

In this paper, an optimal scheduling approach for MES consisting of gas, electricity, cooling and heat networks considering load replaceability is designed to maximize the economic potential of operation. For the convenience of analysis of the influence of load replaceability on the optimal scheduling, load replaceability index (LRI), which measures the overall load replaceability of MES, is analyzed in respect to the optimal scheduling results. The intrinsic relations between load replacement capability and optimum under various operating conditions are investigated to improve the scheduling performance by observing specific rules of load replaceability adjustment.

The remaining of the paper is organized as follows: Section 2 addresses load replaceability and the definition of LRI in MES; with the aid of load replaceability, Section 3 presents the economic scheduling model for an integrated electricity, gas, heat and cooling networks; Section 4 studies the relation between optimal scheduling results and load replacement capability; concluding remarks are given in Section 5.

#### **2. Load Replaceability in Multi-Energy System**

In this section, load replaceability and its relation with optimal scheduling of MES are systematically analyzed. The definition and calculation method of LRI are presented.

#### *2.1. Load Replaceability and Composite Load*

Load replaceability is used to characterize the selection margin of multi-energy consumption. If one specific load can choose two or more than two types of energy flow to meet its own load demand, and then this load is called replaceable load. From the perspective of consumers, the greater the number of energy flow types is and the greater the capacity of multi-generators is, the more the selection margin of multi-energy consumption is.

Generally speaking, loads in MES can be categorized as cooling, heat, power and gas load, corresponding to their respective type of energy network. Usually, due to the natural barrier between heterogenous energy systems in the isolated operating mode, the load can only be supplied by homogenous energy flow. Nevertheless, with the coupling and integration of heterogenous systems assisted by energy conversion systems, it should be emphasized that the load can be satisfied by heterogenous energy flow as well.

In this section, heat load is taken as an example to illustrate the replaceability phenomenon. In independent heat network, heat load can only be supported by the heat flow through steam pipelines. While heat load in MES can be supplied by various forms of energy. For example, on the one hand heat energy can still be transmitted through pipelines just like in independent heat network; on the other hand electric furnaces or gas heaters can transform electric or gas into heat energy. In this case, in addition to the thermal attribute, heat load has electric attribute (electric heater) and gaseous attribute (gas furnace). Electricity and gas are just intermediate transition states. Seen from the exterior, the final state is still heat load since heat flow is the only energy form in the consumption phase.

Though the final load state belongs to heat load in either segregated independent system or integrated dependent system, the introduction of intermediate transition states makes the load property change from monotype load in independent system into composite load in dependent system. The main

specialty of composite load is flexibility, especially the market response flexibility. The operator can adjust the proportion of each component in the composite load, based on the market price of each form of energy flow at different operating interval, thus reducing operation cost. The emergence of composite load enhances the optimization space in the consumption phase and further increases economic benefits. In order to vividly illustrate how load replaceability (composite load) influences the operation economy, an exemplary integrated gas power and cooling system in Figure 1 is briefly studied.

**Figure 1.** Schematic diagram of two-area classic system.

Suppose the current market prices of gas electricity and heat are 6 \$/kWh, 1 \$/kWh and 15 \$/kWh, respectively. And the amount of heat load demand (100 kWh) is much more than that of electricity or gas load demand, which is set by zero for simplicity. The maximal transmission capacity of each network is 50 kW. The following two scenarios are considered for operation cost calculation.


In Scenario 1, the operation cost is

$$(6 \times 50 \times 0.8 \times 0.4 + 15 \times (100 - 50 \times 0.8 \times 0.4) = 1356. \tag{1}$$

In Scenario 2, the operation cost is

$$1 \times 50 \times 0.5 + 6 \times 50 \times 0.5 + 15 \times (100 - 50) = 925. \tag{2}$$

As can be seen, due to the load replaceability, consumers can satisfy the demand by consuming much cheaper transformed energy. In this case it is the transformed heat flow from electric heaters and gas furnaces, thus reducing the overall operation cost. Nevertheless, the multi-generation plant (the CHP plant) can only supplement heat flow to a very limited degree. Heat flow is just a byproduct during the intermediate generation section and the multi-generation plant is not ultimately consumption-oriented, which cannot sufficiently use the price difference between different energy forms to reduce operation cost. As is shown by the results in (1) and (2), operation cost (925\$) considering load replaceability is obviously smaller than that is only with multi-generation plants.

#### *2.2. Load Replaceability Index*

In Section 2.1, basic information of load replaceability and composite load is addressed. From the deduction of a simple 2-scenario integrated energy system, it can be learned that load replaceability enhances the optimization space from the consumption side. In this section, instead of focusing on substitutes of one specific type of load, we give some quantitative indices to evaluate the broad system-level load replaceability. The LRI can reflect to what extent one type of load (the final load state) can be substituted by the other (the intermediate load state).

LRI can be defined from two perspectives: (1) the potential load replaceability and (2) the actual load replaceability. The former evaluates the potential of load substitution of the system before it is ever put into operation. The latter evaluates the actual load substitution ability during the operation.

#### 2.2.1. Load Replaceability Index Reflecting Potential Load Substitution Ability

As for an integrated gas, power, heat and cooling system, there exist four possible types of load in the system. Denote the maximal convertible load demand of the *i*th composite load by *Lmax <sup>i</sup>* , and denote the load component for the *i*th load from each network by *Xload <sup>i</sup>* , *<sup>P</sup>load <sup>i</sup>* and *<sup>G</sup>load <sup>i</sup>* . *<sup>X</sup>load <sup>i</sup>* can either be *Cload <sup>i</sup>* or *<sup>H</sup>load <sup>i</sup>* based on the final energy form. LRI of the *i*th load can be calculated by

$$\alpha\_{i} = \frac{\left(x\_{i}X\_{i}^{load} + u\_{i1}p\_{i}P\_{i}^{load} + u\_{i2}g\_{i}G\_{i}^{load} - L\_{i}^{\max}\right)}{L\_{i}^{\max}},\tag{3}$$

where *xi*, *pi*, and *gi* represents the energy efficiency ratio of respective intermediate load. Binary variables *ui*1, *ui*<sup>2</sup> ∈ {0, 1} represent *i*th load is (*ui*1, *ui*<sup>2</sup> = 1) or is not (*ui*1, *ui*<sup>2</sup> = 0) supplied by specific intermediate loads (*Pload <sup>i</sup>* or *<sup>G</sup>load <sup>i</sup>* ). Suppose there exist *m* loads in the integrated energy system, the overall LRI can be expressed as

$$\alpha = \frac{\sum\_{i=1}^{m} \alpha\_i L\_i^{\text{max}}}{\sum\_{i=1}^{m} L\_i^{\text{max}}}.\tag{4}$$

LRI in (4) actually uses the proportion of remaining alternative load reserve in the maximal load demand to quantify the substitution potential.

#### 2.2.2. Load Replaceability Index Reflecting Actual Load Substitution Ability

As can be seen from (4), all the parameters are actually the upper bounds of corresponding replaceable loads. That is to say, (4) reflects the remaining load conversion reserve at critical states (the maximum values). Nevertheless, replaceable loads in the actual scheduling would not certainly be equal to the maximum; values of replaceable loads should depend on the scheduling, which will be presented in detail in the next section. Therefore, a modified online scheduling-based LRI of the *i*th load is expressed by

$$\beta\_i = \frac{\mu\_{i1} p\_i P\_i^{load} + \mu\_{i2} g\_i G\_i^{load}}{L\_i},\tag{5}$$

where *Pload <sup>i</sup>* and *<sup>G</sup>load <sup>i</sup>* represent the actual intermediate load amount in the scheduling; *Li* represents the actual composite load demand. Then the system-level LRI can be expressed by

$$\beta = \frac{\sum\_{i=1}^{m} \beta\_i L\_i}{\sum\_{i=1}^{m} L\_i} \,. \tag{6}$$

Notice that the replaceable load term *xiXload <sup>i</sup>* disappears in (5); it is because that *xiXload <sup>i</sup>* + *ui*<sup>1</sup> *piPload <sup>i</sup>* + *ui*2*giGload <sup>i</sup>* = *Li* in the scheduling. If (3) is still used, LRI is always zero. In this case only the proportion of energy converted by intermediate loads in the composite load is used to quantify the degree of flexibility.

#### **3. Optimal Economic Scheduling of Multi-Energy System Considering Load Replaceability**

It can be learned from Section 2 that load replaceability offers redundant energy flow path and the optimization space of multi-energy complementary coordination. In response to the market price of different energy flow, the consumption side can participate in economic dispatch by adaptively adjusting the proportion of the intermediate load of various energy form in the composite load. Unlike conventional demand response program in power systems, because of the load replaceability capability through conversion systems, this multi-energy demand response in MES does not cause the change in the ultimate total energy consumption. The consumers only choose the intermediate energy source from which they buy to obtain the ultimate energy rather than change the consumption pattern like they do in electric utility demand response programs. That is to say, operation economy along with customer satisfaction on energy supply is guaranteed.

In this section, by considering load replaceability the economic dispatch model is established for MES in Figure 2 consisting of gas, power, and cooling and heat networks.

**Figure 2.** Schematic diagram of a multi-energy system (MES) consisting of gas, power, and cooling and heat networks. CHP = combined heat and power; ESS = energy storage system.

#### *3.1. Objective Function of Optimal Economic Scheduling of MES*

From the perspective of economy, operation cost mainly comes from three aspects including the cost of purchasing gas from external gas networks, the cost of purchasing electricity from external electric grids and depreciation expense of the multi-energy storage systems.

The cost of gas purchase *OG* is written by the accumulative product of gas price and purchase amount (which is procured from the external gas source shown in Figure 2) during the whole dispatching intervals

$$O\_G = \sum\_{t=1}^{T} r\_{\text{gas}}\left(t\right) \frac{G\_{\text{ext}}\left(t\right)}{R\_{\text{gas}}} t,\tag{7}$$

where *t* represents the dispatch interval; *T* represents the total number of dispatch intervals; *rgas* (*t*) represents gas price at time *t*; *Gext* (*t*) represents gas purchase amount at time *t*; *Rgas* = 9.77 *kw* · *h m*3 represents the calorific value.

The supplementary electricity is procured from the external main grid as shown in Figure 2. Similarly, we can compute the cost of electricity purchase *OP* as

$$O\_P = \sum\_{t=1}^{T} r\_{el\varepsilon}\left(t\right) P\_{\varepsilon\text{ext}}\left(t\right) t\_\prime \tag{8}$$

where *t* and *T* have the same meaning as in (7); *rele* (*t*) represents electricity price at time *t*; and *Pext* (*t*) represents electricity purchase amount at time *t*.

The depreciation cost steams from ageing of storage system during the repeated energy charge/discharge processes. In this paper, it is assumed that gas is either directly supplied to gas load or converted to other energy flow, and there is no gas storage system. For simplicity, it is expressed as the linear combination of storage amount of different energy flow:

$$O\_S = \sum\_{t=1}^{T} \left[ r\_{PS} P\_{PS}^S(t) + r\_{HS} H\_{HS}^S(t) + r\_{CS} \mathcal{C}\_{CS}^S(t) \right],\tag{9}$$

where *rPS*, *rHS*, and *rCS* represent the depreciation rate per unit quantity of power, and heat and cooling energy; *P<sup>S</sup> PS* (*t*), *<sup>H</sup><sup>S</sup> HS* (*t*) and *<sup>C</sup><sup>S</sup> CS* (*t*) represent the electricity, and heat and cooling energy stored at time *t*.

#### *3.2. Constraints of Optimal Economic Scheduling of MES*

#### 3.2.1. Energy Supply and Demand Balance

Similar to power balance equation in electric grids, the energy supply and demand balance should be guaranteed at the first place. Nevertheless, due to the coupling enabled by conversion systems, the energy generation side might contain transformed energy flow (the output of conversion systems), the energy consumption side might contain energy flow to be transformed (the input of conversion systems), which is exactly the intermediate load previously addressed in Section 2. Therefore, under the circumstances of MES, the load for each network can be categorized by convertible load and inconvertible load, which are equivalent to intermediate and ultimate load in Section 2.

The energy supply and demand balance equation for the cooling system can be written by

$$\mathbf{C}\_{AR}\left(t\right) + \mathbf{C}\_{ER}\left(t\right) + \mathbf{C}\_{\mathbf{CS}}^{R}\left(t\right) - \mathbf{C}\_{\mathbf{CS}}^{S}\left(t\right) = \mathbf{C}\_{load}^{D}\left(t\right) + \mathbf{C}\_{load}^{\text{PC}}\left(t\right) + \mathbf{C}\_{load}^{\text{GC}}\left(t\right),\tag{10}$$

where *t* represents the dispatch interval; *CAR* (*t*) represents cooling energy generated by the absorption refrigerator; *CER* (*t*) represents cooling energy generated by the electric refrigerator; *C<sup>R</sup> CS* (*t*) represents cooling energy released by cooling energy storage system; *C<sup>S</sup> CS* (*t*) represents cooling energy stored by cooling energy storage system; *C<sup>D</sup> load* (*t*) represents inconvertible cooling load; *<sup>C</sup>PC load* (*t*) represents convertible electricity-based cooling load; and *CGC load* (*t*) represents convertible gas-based cooling load.

Similarly, the energy supply and demand balance equation for the heat system can be written by

$$H\_{HE}\left(t\right) + H\_{GB}\left(t\right) + H\_{EH}\left(t\right) + H\_{HS}^R\left(t\right) - H\_{HS}^S\left(t\right) = H\_{load}^D\left(t\right) + H\_{load}^{PH}\left(t\right) + H\_{load}^{GH}\left(t\right), \tag{11}$$

where *HHE* (*t*) represents heat energy generated by heat exchanger; *HGB* (*t*) represents heat energy generated by gas boiler; *HEH* (*t*) represents heat energy generated by electric heater; *H<sup>R</sup> HS* (*t*) represents heat energy released by heat energy storage system; *H<sup>S</sup> HS* (*t*) represents heat energy stored by heat energy storage system; *H<sup>D</sup> load* (*t*) represents inconvertible heat load; *<sup>H</sup>PH load* (*t*) represents convertible power-based heat load; and *HGH load* (*t*) represents convertible gas-based heat load.

The energy supply and demand balance equation for the gas system can be written by

$$\left(\mathbf{G}\_{\rm crit}\left(\mathbf{t}\right) - \mathbf{G}\_{\rm CHF}\left(\mathbf{t}\right) - \mathbf{G}\_{\rm GT}\left(\mathbf{t}\right) - \mathbf{G}\_{\rm GB}\left(\mathbf{t}\right) = \mathbf{G}\_{\rm load}^{D}\left(\mathbf{t}\right) + \mathbf{G}\_{\rm load}^{\rm GC}\left(\mathbf{t}\right) + \mathbf{G}\_{\rm load}^{\rm GH}\left(\mathbf{t}\right),\tag{12}$$

where *GCHP* (*t*) represents gas energy consumed by the CHP plant; *GGT* (*t*) represents gas energy consumed by the gas turbine; *GGB* (*t*) represents gas energy consumed by the gas boiler; *G<sup>D</sup> load* (*t*) represents inconvertible gas load; *GGC load* (*t*) represents convertible cooling-oriented gas load; and *GGH load* (*t*) represents convertible heat-oriented gas load.

The energy supply and demand balance equation for the power system can be written by

$$P\_{\rm eff} \left( t \right) + P\_{\rm CHF} \left( t \right) + P\_{\rm GT} \left( t \right) + P\_{\rm PS}^{R} \left( t \right) - P\_{\rm PS}^{S} \left( t \right) - P\_{\rm EH} \left( t \right) - P\_{\rm ER} \left( t \right) = P\_{\rm load}^{D} \left( t \right) + P\_{\rm load}^{\rm PC} \left( t \right) + P\_{\rm load}^{PH} \left( t \right), \tag{13}$$

where *PCHP* (*t*) represents electricity energy generated by the CHP plant; *PGT* (*t*) represents electricity energy generated by the gas turbine; *P<sup>R</sup> PS* (*t*) represents electricity energy released by the electricity energy storage system; *P<sup>S</sup> PS* (*t*) represents electricity energy stored by the electricity energy storage system; *PEH* (*t*) represents electricity energy consumed by the electric heater; *PER* (*t*) represents electricity energy consumed by the electric refrigerator; *P<sup>D</sup> load* (*t*) represents inconvertible power load; *PPC load* (*t*) represents convertible cooling-oriented power load; and *<sup>P</sup>PH load* (*t*) represents convertible heat-oriented power load.

As can be seen from (10)–(13), there exist inconvertible load *C<sup>D</sup> load* (*t*), *<sup>H</sup><sup>D</sup> load* (*t*), *<sup>G</sup><sup>D</sup> load* (*t*), and *P<sup>D</sup> load* (*t*); meanwhile, the convertible load during the intermediate transition state are *<sup>C</sup>PC load* (*t*), *<sup>C</sup>GC load* (*t*), *HPH load* (*t*), *<sup>H</sup>GH load* (*t*), *<sup>G</sup>GC load* (*t*), *<sup>G</sup>GH load* (*t*), *<sup>P</sup>PC load* (*t*), and *<sup>P</sup>PH load* (*t*). Moreover, the first four *<sup>C</sup>PC load* (*t*), *<sup>C</sup>GC load* (*t*), *HPH load* (*t*), and *<sup>H</sup>GH load* (*t*) are actually the final state of the convertible load, while the last four *<sup>G</sup>GC load* (*t*), *GGH load* (*t*), *<sup>P</sup>PC load* (*t*), and *<sup>P</sup>PH load* (*t*) are actually the initial state of the convertible load.

The total amount of energy of composite load in the final energy form (cooling or heat) can thus be expressed by

$$\begin{cases} p\_{load}^{\text{PC}} P\_{load}^{\text{PC}}(t) + C\_{load}^{\text{PC}}(t) = L\_{load}^{\text{PC}}(t) \\ g\_{l@0f}^{\text{GC}} G\_{l@0d}^{\text{C}}(t) + C\_{l@0d}^{\text{GC}}(t) = L\_{load}^{\text{GC}}(t) \\ p\_{load}^{\text{ID}} P\_{load}^{\text{PR}}(t) + H\_{load}^{\text{HR}}(t) = L\_{load}^{\text{PR}}(t) \\ g\_{load}^{\text{CH}} G\_{load}^{\text{CH}}(t) + H\_{load}^{\text{CH}}(t) = L\_{load}^{\text{CH}}(t) \end{cases} \tag{14}$$

where *pPC load* represents the energy efficiency ratio of convertible cooling-oriented power load; *<sup>L</sup>PC load* (*t*) represents the composite power and cooling load; *gGC load* represents the energy efficiency ratio of convertible cooling-oriented gas load; *LGC load* (*t*) represents the composite power and gas load; *<sup>p</sup>PH load* represents the energy efficiency ratio of convertible heat-oriented power load; *LPH load* (*t*) represents the composite power and heat load; *gGH load* represents the energy efficiency ratio of convertible heat-oriented gas load; and *LGH load* (*t*) represents the composite gas and heat load.

#### 3.2.2. Energy Balance of Conversion Systems

As can be seen from Figure 2, the conversion systems comprise the CHP plant, gas turbine, gas boiler, electric heater and electric furnace. In this section, the input and output energy balance of respective conversion system is presented.

The CHP plant is the use of natural gas to generate electricity and useful heat simultaneously. The relation between input gas power and output electricity or heat power of the back pressure CHP plant is

$$P\_{\rm CHP}\left(t\right) = p\_{\rm CHP} G\_{\rm CHP}\left(t\right) \tag{15}$$

$$Q\_{CHP}\left(t\right) = q\_{QP} P\_{CHP}\left(t\right),\tag{16}$$

where *pCHP* represents the power generation efficiency of the CHP plant; *qQP* represents the ratio of electricity-to-heat; *GCHP* represents the input gas power of the CHP plant; and *QCHP* represents the surplus heat, which is transformed to cooling and heat energy through the absorption refrigerator and heat exchanger, respectively.

$$\mathbb{C}\_{AR}\left(t\right) = \eta\_{AR}\left(1 - \gamma\_{HC}\right)Q\_{CHP}\left(t\right),\tag{17}$$

$$H\_{HE}\left(t\right) = \eta\_{HE}\gamma\_{H\mathcal{C}}Q\_{CHP}\left(t\right). \tag{18}$$

where *ηAR* and *ηHE* represent the energy efficiency ratio of the absorption refrigerator and heat exchanger; *γHC* represents the proportion ratio.

The gas turbine is used to transform gas energy into electricity energy

$$P\_{GT}\left(t\right) = \eta\_{GT} G\_{GT}\left(t\right),\tag{19}$$

where *ηGT* represents the energy efficiency ratio of the gas turbine.

The gas boiler is used to transform gas energy into heat energy

$$H\_{GB}\left(t\right) = \eta\_{GB} G\_{GB}\left(t\right) \,. \tag{20}$$

where *ηGB* represents the energy efficiency ratio of the gas boiler.

The electric refrigerator is used to transform electricity energy into cooling energy

$$\mathcal{C}\_{ER}\left(t\right) = \eta\_{ER} P\_{ER}\left(t\right),\tag{21}$$

where *ηER* represents the energy efficiency ratio of the the electric refrigerator.

The electric heater is used to transform electricity energy into heat energy

$$H\_{EH}\left(t\right) = \eta\_{EH} P\_{EH}\left(t\right),\tag{22}$$

where *ηEH* represents the energy efficiency ratio of the electric heater.

3.2.3. Energy Storage Systems

Operating constraints of the energy storage system mainly comprise the storage state and some output constraints no matter it is the electricity, heat or cooling storage system. Take the electricity energy storage system as an example, the typical operating model can be written as

$$\begin{cases} \begin{array}{l} S\_{PS}\left(t\right) = \left(1 - \eta\_{PS}\right) S\_{PS}\left(t-1\right) + \eta\_{PS}^{S} P\_{PS}^{S}\left(t\right) - P\_{PS}^{R}\left(t\right) / \eta\_{PS}^{R} \\\ S\_{PS}^{\min} \le S\_{PS}\left(t\right) \le S\_{PS}^{\max} \\\ S\_{PS}\left(t\_{1}\right) = S\_{PS}\left(t\_{md}\right) \\\ \mathcal{U}\_{PS}^{S}\left(t\right) P\_{PS}^{\min} \le P\_{PS}^{S}\left(t\right) \le \mathcal{U}\_{PS}^{S}\left(t\right) P\_{PS}^{\max} \\\ \mathcal{U}\_{PS}^{R}\left(t\right) P\_{PS}^{\min} \le P\_{PS}^{R}\left(t\right) \le \mathcal{U}\_{PS}^{R}\left(t\right) P\_{PS}^{\max} \\\ \mathcal{U}\_{PS}^{R}\left(t\right) + \mathcal{U}\_{PS}^{S}\left(t\right) \le 1 \\\ \mathcal{U}\_{PS}^{R}\left(t\right), \mathcal{U}\_{PS}^{S}\left(t\right) \in \left(0, 1\right) \end{array} \tag{23}$$

where *SPS* (*t*) represents the remaining power in the electricity energy storage system; *ηPS* represents the attrition rate during the charge/discharge process; *η<sup>S</sup> PS* represents the charge efficiency; *<sup>η</sup><sup>R</sup> PS* represents the discharge efficiency; *U<sup>S</sup> PS* (*t*) represents the on/off charge state; *<sup>U</sup><sup>R</sup> PS* (*t*) represents the on/off discharge state; *S*min *PS* and *<sup>S</sup>*max *PS* represents the capacity limit; *SPS* (*t*1) represents the initial power in the beginning of one cycle of operation; and *SPS* (*tend*) represents the final power in the end of one cycle of operation. When considering the periodicity of the operation of ESS (energy storage system), the final power of the previous cycle should be equal to the initial power of the next cycle; hence, *SPS* (*tend*) = *SPS* (*t*1) should be satisfied.

#### 3.2.4. Other Constraints

In addition to constraints characterizing energy supply and demand balance, input and output energy relation of conversion systems, other mentionable operating constraints are those that describe the output limit of specific components in MES.

Firstly, the convertible load in each energy network has its lower/upper limit:

$$\begin{cases} \begin{aligned} &\frac{\text{P}^{\text{PC}}}{\text{load}} \le \text{P}^{\text{PC}}\_{\text{load}} \left(t\right) \le \overline{\text{P}}^{\text{PC}}\_{\text{load}}\\ &\underline{\text{C}^{\text{PC}}} \le \text{C}^{\text{PC}}\_{\text{load}} \left(t\right) \le \overline{\text{C}}^{\text{PC}}\_{\text{load}}\\ &\underline{\text{P}}^{\text{PH}} \le \text{P}^{\text{PH}}\_{\text{load}} \left(t\right) \le \overline{\text{P}}^{\text{PH}}\_{\text{load}}\\ &\underline{H}^{\text{PH}}\_{\text{load}} \le H^{\text{PH}}\_{\text{load}} \left(t\right) \le \overline{H}^{\text{PH}}\_{\text{load}}\\ &\underline{\text{C}}^{\text{CCC}}\_{\text{load}} \le G^{\text{CC}}\_{\text{load}} \le \overline{G}^{\text{CC}}\_{\text{load}}\\ &\underline{\text{C}}^{\text{CCC}}\_{\text{load}} \le G^{\text{CH}}\_{\text{load}} \left(t\right) \le \overline{\text{C}}^{\text{CC}}\_{\text{load}}\\ &\underline{\text{C}}^{\text{CH}}\_{\text{load}} \le G^{\text{CH}}\_{\text{load}} \le \overline{\text{C}}^{\text{PH}}\_{\text{load}}\\ &\underline{H}^{\text{CH}}\_{\text{load}} \le H^{\text{CH}}\_{\text{load}} \left(t\right) \le \overline{H}^{\text{CH}}\_{\text{load}} \end{aligned} \tag{24}$$

Secondly, the output of conversion system has its lower/upper limit:

$$\begin{cases} \underline{\rm D\_{\rm CHP}} \le P\_{\rm CHP} \left( t \right) \le \overline{P}\_{\rm CHP} \\\ \underline{\rm D\_{\rm GT}} \le P\_{\rm GT} \left( t \right) \le \overline{P}\_{\rm GT} \\\ \underline{\rm H\_{\rm GB}} \le H\_{\rm GB} \left( t \right) \le \overline{H}\_{\rm GB} \\\ \underline{\rm C\_{\rm ER}} \le C\_{\rm ER} \left( t \right) \le \overline{C}\_{\rm ER} \\\ \underline{H\_{\rm EH}} \le H\_{\rm EH} \left( t \right) \le \overline{H}\_{\rm EH} \end{cases} \tag{25}$$

Thirdly, the purchase amount of electricity *Pext* (*t*) and gas *Gext* (*t*) cannot be unlimited; it is assumed that they cannot surpass the allowable maximum.

$$\begin{cases} \ 0 \le P\_{\text{ext}}\left(t\right) \le P\_{\text{ext}}^{\text{max}}\\\ 0 \le G\_{\text{ext}}\left(t\right) \le G\_{\text{ext}}^{\text{max}} \end{cases} \tag{26}$$

Based on (26) it can be learned that the minimal purchase amount from external electricity or gas network is nonnegative. Particularly with respect to gas purchase, since gas is the only crude source which is used as the primitive input energy of all conversion systems and multi-generation plants, it cannot be less than zero. Otherwise, it means MES can autonomously inject gas toward external network out of nothing, which is against the law of conservation of energy. On the other hand, the electricity purchased takes the role of auxiliary, which can act as supplementary energy provider in due time. Because of the existence of gas load, electricity purchased cannot single-handedly solve the energy supply problem. Nevertheless, the energy supplier for heat and cooling load could be electricity purchased.

The model (7)–(26) belongs to mixed integer linear programming (MILP) and is solved by commercial optimization software package CPLEX under general algebraic modelling system (GAMS).

#### **4. Case Study**

In this section, MES in Figure 2 is used for numerical analyses. Specifically, the intrinsic relation between load replaceability (as is quantified by (4) and (6)) and optimal scheduling results is studied, and thus giving a priori knowledge of improving economic scheduling performance by reformulating the composite load.

#### *4.1. Simulation Settings*

Parameters of main system components are given in Table 1.


**Table 1.** Parameters and corresponding values of MES.

The load status described by loading rate at every scheduling interval is shown in Figure 3.

**Figure 3.** Hourly load status of gas, power, and heat and cooling load of MES.

The rated power of cooling, heat, and power and gas load is 11,000 kW, 8000 kW, 15,000 kW, and 10,000 kW, respectively. The day-ahead price information from the electricity and gas market is given in Table 2.



#### *4.2. Simulation Results and Discussions*

After presenting the basic simulation settings, optimal scheduling and its relation with load replaceability is studied by calculating the scheduling model in (7)–(9) in Section 3 and LRI in (3)–(6) in Section 2.2. As previously addressed in Section 2.2, there exist two ways of describing the load replacement capability of MES; therefore, by special setting of composite load in (14), we change the values of LRI both in (4) and (6), and analyze the resulting influence upon optimal scheduling results through a sensitivity analysis tests. Specifically, the following 2 scenarios are considered.

#### 4.2.1. Scenario 1

It is supposed that the composite load is of electricity-cooling and electricity-heat types. As for composite electricity-cooling load, the rated capacity of alternative electricity load is fixed to 2000 kW, while the capacity of alternative cooling load is increased from 600 kW to 6000 kW at the step size of 600 kW. Similarly, the rated capacity of alternative electricity load is fixed to 3750 kW in composite electricity-heat load, while the the capacity of alternative heat load is increased from 300 kW to 3000 kW at the step size of 300 kW. Consequently, 10 LRIs can be calculated based on (4) and (6), and the results are shown in Figure 4.

**Figure 4.** load replaceability index (LRI) with and without optimal scheduling in Scenario 1.

As can be seen from the blue curve in Figure 4, with the increase of alternative cooling and heat load, the remaining alternative load reserve is elevated from 900 kW to 9000 kW, thus enhancing LRI from 0.1 to 1. Meanwhile, the red curve shows that the proportion of energy converted by intermediate loads in the actual scheduling drops as the rated capacity of alternative heat (cooling) load increases, which means the operator chooses to obtain heat or cooling energy through redundant conversion paths. Take the electricity-cooling load as an example, the redundant conversion path is shown in dashed black line in Figure 5.

**Figure 5.** Redundant conversion path of electricity-cooling load in Scenario 1.

Due to the fact that energy efficiency ratio of electric refrigerator is much higher than that of alternative electricity load in the composite electricity and cooling load. The operator chooses redundant path such that more cooling energy can be harvested with equal amount of electricity energy. Similar conclusions can be made in respect to composite electricity and heat load. Therefore, the proportion of energy converted by intermediate electricity load in the composite load is decreased with the increase of transmission capacity of redundant paths, which is equivalent to the drop of LRI (from 0.86 to 0.35). Furthermore, the relation between LRI and optimal scheduling results in Scenario 1 is shown in Figures 6 and 7.

**Figure 6.** Relation between offline potential-based LRI and optimal scheduling results in Scenario 1. (**a**) The relation between LRI and optimal energy cost. (**b**) The relation between LRI and optimal cost of electricity purchase. (**c**) The relation between LRI and optimal cost of gas purchase.

**Figure 7.** Relation between online scheduling-based LRI and optimal scheduling results in Scenario 1. (**a**) The relation between LRI and optimal energy cost. (**b**) The relation between LRI and optimal cost of electricity purchase. (**c**) The relation between LRI and optimal cost of gas purchase.

As can be seen from Figure 6, the cost including cost of gas and electricity purchase decreases with the increase of offline potential-based LRI. Contrarily, there exists a positive correlation between energy cost and online scheduling-based LRI. It verifies that the optimization space can be enhanced by elevating the capacity of alternative load, thus enabling the decrease of minimal scheduling costs.

#### 4.2.2. Scenario 2

It is supposed that the composite load is of gas-cooling and gas-heat types. As for composite gas-cooling load, the rated capacity of alternative gas load is fixed to 5000 kW, while the capacity of alternative cooling load is increased from 600 kW to 6000 kW at the step size of 600 kW. Similarly, the rated capacity of alternative gas load is fixed to 6000 kW in composite gas-heat load, while the the capacity of alternative heat load is increased from 300 kW to 3000 kW at the step size of 300 kW. Consequently, 10 LRIs can be calculated based on (4) and (6), and the results are shown in Figure 8.

**Figure 8.** LRI with and without optimal scheduling in Scenario 2.

The relation between LRI and optimal scheduling results in Scenario 2 is shown in Figures 9 and 10. Similar to Scenario 1, it can be learned that the optimal energy cost can be decreased by elevating the offline potential based LRI, which quantifies the maximum alternative load capacity of the system. In addition, it can be found that the optimal cost does not decrease when LRI reaches the critical state. It means that there is no unrestricted expansion of optimization space by single-handedly elevating alternative load capacity, which could improve the expectancy of optimum under the joint force of other factors such as system configuration and market price signals.

**Figure 9.** Relation between offline potential-based LRI and optimal scheduling results in Scenario 2. (**a**) The relation between LRI and optimal energy cost. (**b**) The relation between LRI and optimal cost of electricity purchase. (**c**) The relation between LRI and optimal cost of gas purchase.

**Figure 10.** Relation between online scheduling-based LRI and optimal scheduling results in Scenario 2. (**a**) The relation between LRI and optimal energy cost. (**b**) The relation between LRI and optimal cost of electricity purchase. (**c**) The relation between LRI and optimal cost of gas purchase.

#### **5. Conclusions**

In this paper, optimal economic scheduling of multi-energy systems considering load replaceability is presented. Two types of LRI reflecting the offline potential and online scheduling performance are designed, and the relation between LRI and optimal scheduling results is analyzed. Through numerical analysis, it can be found that optimal scheduling performance can be improved by increasing the alternative load capacity (offline potential-based LRI). And this improvement is not unrestricted and could only maintain by changing load replaceability in a limited range, after which the optimum reaches the saturation phase.

**Author Contributions:** T.Z. and Z.D. conceived and designed the experiments; J.Y. performed the experiment; T.Z. wrote the paper; Y.Y. and J.C. reviewed the paper.

**Funding:** This work was supported by the National Key Research and Development Program of China (2018YFB0905000) and the project of 'Study and application of multi-energy micro-grid source-network-load-storage coordinated control system' of NARI Group Corporation (State Grid Electric Power Research Institute).

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

#### **References**


c 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
