Optimal Scheduling of an Regional Integrated Energy System with Energy Storage Systems for Service Regulation

Ancillary services are critical to maintaining the safe and stable operation of power systems that contain a high penetration level of renewable energy resources. As a high-quality regulation resource, the regional integrated energy system (RIES) with energy storage system (ESS) can effectively adjust the non-negligible frequency offset caused by the renewable energy integration into the power system, and help solve the problem of power system frequency stability. In this paper, the optimization model aiming at regional integrated energy system as a participant in the regulation market based on pay-for-performance is established. Meanwhile YALMIP + CPLEX is used to simulate and analyze the total operating cost under different dispatch modes. This paper uses the actual operation model of the PJM regulation market to guide the optimal allocation of regulation resource in the regional integrated energy system, and provides a balance between the power trading revenue and regulation market revenue in order to achieve the maximum profit.


Introduction
With a large number of Renewable Energy Sources (RES) integrated into the power grid and the construction of the Energy Internet, it can be expected that by 2050, renewable energy generation as a proportion of total energy consumption may reach 80%. The World Wide Fund for Nature (WWF) predicts that the world will be powered by 100% renewable energy by the middle of this century. Such a change in power and energy systems will bring significant benefits in terms of economic efficiency, energy supply security, and environmental protection [1]. Renewable energy usually has the characteristics of intermittence and volatility. Non-negligible frequency offset may appear under high penetration level of renewable energy, which causes power system frequency stability issues. Therefore, additional regulation capability is required to provide regulation services. The Regional Integrated Energy System (RIES) with electrical energy storage system (ESS), a concept originated from the Energy Internet, has been employed as an effective technology to solve this problem [2,3].
The Energy Internet can be divided into two parts, energy primary network and regional integrated energy system [4][5][6][7][8]. The initial purpose of using ESS in the RIES is to generate profits by utilizing energy price arbitrage opportunities, known as load shifting or peak shaving [9][10][11]. It solves the renewable energy consumption problem and makes arbitrage by charging energy when prices are 2. Regional Integrated Energy System 2.1. Structure of the Regional Integrated Energy System ETH Zurich proposed the concept of Energy Hub and Energy Interconnect [17]. The Energy Hub model can represent the generation and transformation between various energies. Energy Hub is an abstract model of multiple practical devices and their relationships. Energy Hub has been widely used in many different fields, such as optimal design and optimal scheduling [18,19]. An Energy Hub is considered as a unit where multiple energy resources can be converted, conditioned, and stored. It represents an interface between different energy infrastructures and/or loads. Energy hubs consume power from their input ports which are connected to generation units such as electricity and natural gas infrastructures, and provide certain required energy services such as electricity, heating, cooling, compressed air, etc. through the output ports. The RIES is a multi-energy flow supply system that connects natural gas networks, distributed power generation, CCHP systems, energy storage devices, etc. through the use of information and communication technologies, so Energy Hub can be used to describe energy coupling relationship in RIES. In Figure 1, the structure diagram of Energy Hub of the RIES includes electrical/gas energy inputs and power/thermal energy outputs [20]. The figure shows the connection between the various devices in the RIES and the energy flow. In RIES, not only various types of distributed generation equipment (photo-voltaic Cell, PV; wind turbine, WT; gas turbine, GT) exist, but also it contains various forms of energy input and output (gas boiler, GB). In addition, the system also contains energy storage equipment (energy storage, ES; thermal storage, TS) and energy conversion equipment (waste heat boiler, WHB; electric refrigerator, ER; electric boiler, EB). In order to accurately describe the composition of the RIES and the energy flow between devices, the expression of Wang et al. [16] is used for reference. Classifying the buses according to the different types of energy transfer media, buses included in this paper are electrical buses, air buses, smoke buses, hot water buses and steam buses respectively. The bus structure of a gas driving CCHP is shown in Figure 2. In order to accurately describe the composition of the RIES and the energy flow between devices, the expression of Wang et al. [16] is used for reference. Classifying the buses according to the different types of energy transfer media, buses included in this paper are electrical buses, air buses, smoke buses, hot water buses and steam buses respectively. The bus structure of a gas driving CCHP is shown in Figure 2.

The Role of the RIES Energy Storage Equipment in the Regulation Market
The frequency of the system can be adjusted by controlling the charge and discharge process ESS of RIES [21]. Different from conventional thermal power plants, ESS can be economically operated in the non-full load mode and provide the adjusting ability twice of its capacity. This feature is ideal for providing the power system regulation service. The comparison of regulation ability between EES and the traditional unit is shown in Figure 3. Compared with traditional generators,

The Role of the RIES Energy Storage Equipment in the Regulation Market
The frequency of the system can be adjusted by controlling the charge and discharge process ESS of RIES [21]. Different from conventional thermal power plants, ESS can be economically operated in the non-full load mode and provide the adjusting ability twice of its capacity. This feature is ideal for providing the power system regulation service. The comparison of regulation ability between EES and the traditional unit is shown in Figure 3. Compared with traditional generators, there are some significant advantages of employing energy storage equipment to provide regulation services. First is the fast response and high ramp rate. It also has a high adjustment accuracy. It can be flexibly switched between the output state and the fully discharging state (or the charging state) with a high accuracy of actions. Hence some malfunctions can be avoided such as adjusting to the opposite direction when the thermal power unit respond to the automatic generation control (AGC) signal [22].

The Role of the RIES Energy Storage Equipment in the Regulation Market
The frequency of the system can be adjusted by controlling the charge and discharge process ESS of RIES [21]. Different from conventional thermal power plants, ESS can be economically operated in the non-full load mode and provide the adjusting ability twice of its capacity. This feature is ideal for providing the power system regulation service. The comparison of regulation ability between EES and the traditional unit is shown in Figure 3. Compared with traditional generators, there are some significant advantages of employing energy storage equipment to provide regulation services. First is the fast response and high ramp rate. It also has a high adjustment accuracy. It can be flexibly switched between the output state and the fully discharging state (or the charging state) with a high accuracy of actions. Hence some malfunctions can be avoided such as adjusting to the opposite direction when the thermal power unit respond to the automatic generation control (AGC) signal [22].

The Regulation Service Based on Pay-for-Performance
Regulation service is a reliability product in the ancillary services market that gives marketbased compensation to resources. Regulation service can adjust output or consumption in response to an automatic generation control (AGC) signal. Its main goal is to help to maintain interconnection

The Regulation Service Based on Pay-for-Performance
Regulation service is a reliability product in the ancillary services market that gives market-based compensation to resources. Regulation service can adjust output or consumption in response to an automatic generation control (AGC) signal. Its main goal is to help to maintain interconnection frequency, as well as to help track moment-to-moment fluctuations in customer loads, and to correct the unintended fluctuations in generation and manage differences between actual and scheduled power flow between control areas. The regulation market based on pay-for-performance provides a market for regulation resources. All kinds of regulation resources participants into this market for the purchase and sale of the regulation service to make revenue based on their regulation performance.
Besides, comparing with the old way of compensation (the regulation service based on pay-for-capacity), the regulation service based on pay-for-performance means that the regulation service is compensated depending on the regulation capacity and performance both. So, the Benefits Factor η b f , the Regulation Mileage Rate η mileage , and the Historical Performance Score η h_ps are introduced to adjust the quotation. When there are multiple regulation resources, the market uses the regulation Benefits Factor and the Historical Performance Score to adjust the quoted price of the regulation resource, so that the regulation resource with the best performance is more likely to be selected. The framework of regulation service process based on pay-for-performance is illustrated in Figure 4. The regulation market clearing engine optimizes resource energy schedules and the regulation market clearing prices for operation hour while respecting appropriate data from the regulation service participant and the PJM market. compensated depending on the regulation capacity and performance both. So, the Benefits Factor , the Regulation Mileage Rate , and the Historical Performance Score ℎ_ are introduced to adjust the quotation. When there are multiple regulation resources, the market uses the regulation Benefits Factor and the Historical Performance Score to adjust the quoted price of the regulation resource, so that the regulation resource with the best performance is more likely to be selected. The framework of regulation service process based on pay-for-performance is illustrated in Figure 4. The regulation market clearing engine optimizes resource energy schedules and the regulation market clearing prices for operation hour while respecting appropriate data from the regulation service participant and the PJM market.

Regulation Signal Type
In the regulation market based on pay-for-performance, there are two different types of ACG control signal, RegA and RegD. RegA is for the traditional regulation resources with the slow regulation rate. It uses traditional power generation to carry out regulation service. Its physical characteristics limit the ramp rate. RegD is for the regulation resources with faster regulation rate, which uses the real-time dynamic signal to carry out regulation service and no ramp rate limitation exists. Hydro power can fulfill the characteristics of both signal types at the same time.
The compensation for regulation resources participate in the ancillary service market is composed of two parts. The first part depends on the regulation capacity and the second part depends on the regulation performance which is determined by regulation mileage, quality and clearing price. Mileage is the summation of movement that a resource is following requested by the regulation control signal. It is the absolute sum of movement of the regulation signal in a given time period, which can be expressed as: (1) Figure 4. The regulation service process base pay-for-performance framework.

Regulation Signal Type
In the regulation market based on pay-for-performance, there are two different types of ACG control signal, RegA and RegD. RegA is for the traditional regulation resources with the slow regulation rate. It uses traditional power generation to carry out regulation service. Its physical characteristics limit the ramp rate. RegD is for the regulation resources with faster regulation rate, which uses the real-time dynamic signal to carry out regulation service and no ramp rate limitation exists. Hydro power can fulfill the characteristics of both signal types at the same time.
The compensation for regulation resources participate in the ancillary service market is composed of two parts. The first part depends on the regulation capacity and the second part depends on the regulation performance which is determined by regulation mileage, quality and clearing price. Mileage is the summation of movement that a resource is following requested by the regulation control signal. It is the absolute sum of movement of the regulation signal in a given time period, which can be expressed as: where M RegA and M RegD represent regulation mileage of different regulation resources. P RegA i and P RegA i−1 are the power differences of the regulation resource in the adjacent time interval. Resources following the dynamic signal (RegD) can change more dramatically than those on traditional signal (RegA). The mileage ratio β mileage is defined as (2). It is the measurement of the relative work (movement) of RegD resources relative to RegA (unit-less). The mileage shows ratio the regulation capability of different regulation resources:

Benefits Factor
The benefits factor η b f shows the mutual substitution rate between two regulation resources (RegD & RegA), so that they can be quantified. The η b f transforms a fast moving resource's MWs into traditional MWs or Effective MWs. The offer price of RegD resources can be adjusted based on η b f . The employment of η b f can be used to sort the regulation resources. Based on the ranking results, the resource with good performance and low cost is preferred. In default situation, the Benefits Factor for all RegA resources is 1. η b f can be calculated as follows [23]: 3.2.1. Calculation of the Regulation Capacity Adjusted Based on the Historical Performance Score where P reg_per f , MW/h is the regulation capacity adjusted based on the historical performance score. P reg_cap , MW/h is the regulation capacity provided by the regulation service resources. η h_ps is the historical performance score for different types of regulation service published by the regulation market based on pay for performance. It can be written as [24]: The performance score η h ps is calculated from three components: the accuracy score, S a , the delay score, and the precision score [15]. Accuracy denotes the correlation or degree of relationship between control signal and system response. Delay score, S d , denotes the time delay between control signal and point of highest correlation. Precision score, S p , denotes the instantaneous error between the control signal and the regulation. At this time, all three components above are weighted equally. Here, a, b, and c are the weight of each factors; in this paper, a = b = c = 1/3.

Calculation of Initial Adjusted Total
Offer where p o f f er_total is the initial adjusted total offer of regulation service. p cap_o f f er , p per f _o f f er are the offer given by the regulation resources based on its regulation capacity and performance respectively. c loc,total is the lost opportunity cost, which is a calculated value based on forgone revenue when a resource is dispatched uneconomically in order to provide regulation service. Its calculation method will be discussed later: where p o f f er represents the initial adjusted total offer of regulation resource that participates in the regulation service. Assuming η b f = 1 for this step of the benefits factor calculation.

Ranking Resources
Put resources in ascending rank order based on Initial Adjusted Total Offer. RegD resources with p o f f er_total equal to zero.

Calculation of Benefits Factor
The benefits factor is a piecewise function. When the benefits factor of the RegD regulation resource is higher, the offer becomes cheaper and the more efficient resources that can participate in the regulation service: Here, α RegD is the proportion of RegD regulation resource in all regulation resources. P reg_req is the market regulation demand.

Calculation of Effective Regulation
Capacity Here , P e f f ective is the effective regulation capacity. In accordance with the ranking, regulation resources are called until the cumulative effective regulation service capacity meet the market regulation demand.

The Lost Opportunity Cost
In general, Lost Opportunity Cost is the difference in net compensation from the energy market between what a resource receives when providing regulation service and what it would have received for providing energy only. The purpose is to quantify the foregone revenue or increased costs incurred by a resource when it is dispatched uneconomically to provide regulation service.
c loc,total consists of three parts, which can be expressed as [25]: c loc,total = c loc,be f ore + c loc + c loc,a f ter (9) c loc,be f ore is the pre-regulation cost, which refers to the increased cost that occurs due to change the output of the regulation resources from the optimal economic operating point to meet the system demand, before the regulation resource participate in regulation. c loc is the cost during the regulation. Because of the output of the regulation resource will change to non-economic operation model during it participation of the regulation. c loc,a f ter is the post-regulation cost, which is the increased that occurs to change the output of the regulation resources from the regulation operating point to optimal economic operating point: Here p f orecasted is the day-ahead location marginal price (LMP). c reg is the generating cost before the regulation resource participates into regulation. When RIES is in optimal economic operating state, the exchange power between RIES and the power grid is P eco . P reg is the exchange power between RIES and the power grid when the regulation resource receives the AGC signal. t is the time required for the regulation resources to adjust from the economic operation model to the regulation model.
Assuming the ramp rate of the regulation resources is ∆P, t can be expressed as: Since the ramp rate of electrical ESS tends to be infinite, t = 0, c loc,be f ore = c loc,a f ter = 0, so the pre-regulation cost and post-regulation cost can be ignored.

Objective Function
The problem of maximizing revenue from a RIES is formulated as a MILP optimization problem [17]. The RIES model presented above is combined with a cost function in order to maximize the revenue from the energy arbitrage and the regulation market. The objective function c total when the RIES participates into arbitrage and regulation is given by: c total = min c e + c gas + c es + c ts − r reg (13) Energies 2018, 11, 195 8 of 19 where c total is the system overall operating cost. c e , c gas is the electricity and natural gas purchase price, respectively. c es is the electrical ESS operating cost. c ts is the thermal ESS operating cost. r reg is the revenue from the regulation market: where T is the number of time periods for the dispatch cycle. ∆t is the time period length of each dispatch period, h; P t grid is the power exchange between RIES and the power grid at time t = 1,2, . . . ,H. p t grid is the locational marginal price of the electricity market, $/MW·h. When the value of p t grid is positive. It indicates that RIES buys electricity from power grid. When it is negative, then RIES should sells electricity to obtain profit. Thus we have: where p t grid/b , p t grid/s represents electricity purchasing price and selling price of RIES system, respectively. Their value is related to the day-ahead price: where the symbol p gas is the unit heat value price of the natural gas, $/MBTu. P n gt and Q m gb represent generated power of nth cogeneration gas turbines and heating power of mth gas boiler, MW.
α MBtu−MWh is the conversion coefficient that convert megawatt hours to British thermal units, the value is 3.4121, from 1 MWh = 3.4121 MBTu. η gb is the heat efficiency of the gas boiler. η gt_e is the electrical efficiency of the gas turbines. It is assumed that the storage device does not reduce its efficiency during operation in the whole operation cycle, the cost of a single complete charge and discharge cycle c cycle is calculated by the following equation. The procurement cost is c purchase . M cycles is the number of reuse times: The operation cost of one dispatch cycle is: where c es−cap , kW·h is the capacity. P t es,d/es,c is charge and discharge power of the ESS. The major part of the heat storage operation cost is the electricity charge for heat water pump, which can be expressed as: Here, the term c ts represents the operating cost of thermal ESS, c t grid is the electricity price of this period in $/MW·h. Q t hs,d/hs,c is the charge and discharge thermal energy of the thermal ESS. In the regulation market based on pay for performance, revenue is settled according to the actual contribution of each participant. In the bidding session, each service resource provider reports the callable capacity, response speed, callable time, bid price and other information. The regulation market selects preferred ancillary service provider according to market rules. Based on the settlement rules Energies 2018, 11, 195 9 of 19 and data of the regulation market based on pay for performance, the revenue objective function of RIES participating in the regulation service is: r t per f = P t reg β t mailage c t rmpcp (22) In arbitrary period t, the revenue of RIES participating into regulation mainly from two aspects, regulation capacity benefits r t cap which provide regulation spare capacity and the regulation performance benefits r t per f , which is shown as in Equations (20)- (22). Here C t reg is the regulation spare capacity provided by RIES in MW. c t rmccp is the regulation market capacity clearing price (RMCCP) of the regulation market. P t reg is the power exchange of RIES participating in the regulation service. When the electrical ESS in RIES is used as regulation service, there are two operating modes, which are charging and discharging. Therefore, we should take the absolute value of P t reg . Symbol β t mailage is the mileage rate of the regulation resource. c t rmpcp is the regulation market performance clearing price (RMCCP) of regulation market.

Constraints
The objective function Equation (13) is subjected to a set of system constraints. All the constraints are summarized as follows:

1.
Electrical bus balance constraint P grid + P gt + P pv + P wt + P es,d = L e + P er + P eb + P es,c + P reg 2.
Smoke bus balance constraint P gt

Air bus balance constraint
In Equations (23)- (27), P grid is the power exchange of RIES and power grid. P gt , P pv , P wt , P es,d , L e are the generate power of GT, output of PV, output of WT, discharge power of EES, electrical load respectively. The symbols P er , P eb , P es,c , P reg , η gt_t are the consumed power of ER, consumed power of EB, charge power of EES, capacity of regulation, thermal power production efficiency of GT, respectively. Q whb , η whb , Q gb , Q ts,d , Q hx , Q ts,c represent the output of WHB, thermal power production efficiency of WHB, thermal production power of GB, thermal power of thermal ESS, heat exchanger thermal output, heat discharge power of thermal ESS respectively. η hx , L th are the heat exchanger efficiency, thermal load. η er represents the refrigeration efficiency of ER, which is the ratio of refrigerating capacity and input electricity power. Q er is the refrigeration power, L c is the cooling load.

Constraints of the Units Output Upper and Lower Limit
represent the electric power, lower and upper electric power limit, thermal power, lower and upper thermal power limit of the device i respectively.

Electricity EES constraints
The S OC (T + 1) and S OC (T) are the state of charge (SOC) at instant T + 1 and T for the electrical ESS, respectively. S OC_min and S OC_max are the lower and the upper limits of the electrical ESS SOC. σ es , η es,c , η es,d denote the self-discharge rate, charging efficiency, discharging efficiency, respectively. C es−cap , P es,c , P es,d are the capacity of the electrical ESS, the power effectively delivered to the grid by the electrical ESS, the power effectively drawn from the grid by the electrical ESS, respectively. I es,c and I es,d are the charge-discharge state variables of the electrical ESS (0-1 variables).

Charge-discharge power constraint of electrical ESS
I es,c P es,c min + I es,d P es,d min ≤ P es,c/d ≤ I es,c P es,c max + I es,d P es,d max (30) Here the terms P es,c min , P es,d min , P es,c max , P es,d max represent the minimum and maximum charge-discharge power of the electricity ESS, respectively.  (31) where E ts (T + 1) and E ts (T) are the storage state at instant T + 1 and T for the thermal ESS, respectively. The terms E ts min , E ts max are the capacity lower and the upper limits of the thermal ESS, respectively. σ ts , η ts,c , η ts,d represent the thermal dissipation rate, charging and discharging efficiency respectively., respectively. Q ts,c and Q ts,d are the thermal charging and discharging power of thermal ESS. I ts,c , I ts,d are the variables of the thermal ESS charge-discharge state (0-1 variables).

Charge-discharge power constraint of thermal ESS
I ts,c Q ts,c min + I ts,d Q ts,d min ≤ Q ts,c/d ≤ I hs,c Q ts,c max + I ts,d Q ts,d max (32) Here, Q ts,c min , Q ts,d min are the minimum thermal charging and discharging power respectively. Q ts,c max , Q ts,d max represent the maximum thermal charge and discharge power of the thermal ESS respectively.

Constraints of the Gas Turbine Output Upper and Lower Limit
Here, P min gt and P max gt are the upper limit and lower limit of gas turbine output, respectively.

Constraint of the Gas Turbine Ramp Rate
Here, the term ∆P gt represents the output difference of the gas turbine between time t and t + 1.

Reserve Constraint
Considering the possibility that the RIES is likely to be disconnected from the primary power gird at any time, to ensure uninterrupted operation of important loads, the maximum output of the RIES in period k should be equal to or higher than the needs of the important load. The requirement of important load can be described as: Simultaneously, to ensure that the important load can run continuously for a period of time t s , the maximum energy provided by RIES in the period ts should be greater than the energy demand of the important load in this period, which can be written as:

Solution
The standard solution model of a mixed integer linear programming problem can be written as: Here, the variables to be optimized include output, conversion equipment input, ESS input and output, electricity purchase and sell. The constraint equation is an energy balance equation of the bus and the ESS. The inequality constraints are the operation constraints of each unit. For the model above, this paper transforms it into a mixed integer linear programming problem (MILP), which is solved by MATLAB by using YALMIP + CPLEX.

Basic Data
A framework is outlined in this paper for calculating the maximum revenue from a RIES that participates in a day-ahead market, i.e., energy arbitrage, and a regulation market based on pay-for-performance.
An actual operating RIES that is located in Guizhou Province, China, was selected as a simulation case study in this paper. The system includes PV, wind power, gas boiler, electrical/thermal ESS and three gas turbines. Among the three gas turbines, two of them (1#, 2#) are located at the high-speed rail station. The minimum down/up time of these two turbines is 2 h/1 h. One of the turbines is located at a hospital (3#), and its minimum down/up time is 3 h/2 h. The length of dispatching time period is 24 h. the length of dispatching time is 1 h. The capacity of each device is shown in Table 1.
In this simulation, we take a Siemens SGT-100 industrial gas turbine as an example. The power generation of this gas turbine is 5.05 MW. Its electrical efficiency is 30.2% and heat rate is 11,292 Btu/kWh, which means that the fuel cost coefficient is 11.292 MBTu/MWh [26]. As for gas purchase price, we use the average price of gas purchase price in the first two quarters of 2017 which is released by the US Energy Information Administration. The purchase price of natural gas for power generation is 3.738 $/MBTu [27].
The capacity, initial energy, the ESS upper and lower limits of the ESS are presented in Table 2. The procurement cost of electrical ESS is 3.03 × 10 6 $ and the number of charge-discharge cycle is 6000 times.
The relevant parameters of the main device are shown in Table 3 [28][29][30][31][32]. Due to the self-discharge rate of lithium battery storage is only about 3% per month [31], 0% can be taken as the self-discharge rate. Renewable energy output curve, RIES electrical/thermal load forecast curve are shown in Figures 5 and 6. Electricity and regulation market price is from the US PJM power market. The price curve of 17 August 2017 is shown in Figure 7. The dispatching time length is 1 h, and the specific values are shown in Table A1 of Appendix A [33]. The historical performance score of different types of regulation resources published by PJM is shown in Table 4.
The sensitive load in RIES is 5 MW, the important load should operate uninterruptedly for at least one hour after the system disconnects from the primary power grid.

Examples Analysis
and can be obtained after taking related parameters of different types of regulation resources into Equations (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12). Then we can rank the results of the and according to the calculations. In this paper, the regulation resource is electrical ESS and gas turbine. The specific parameters are shown in Table 5. In order to reflect the and of different regulation resources, hydropower and the Demand Side Resource (DSR) datas are added for comparison into Table 5.

Examples Analysis
and P e f f ective can be obtained after taking related parameters of different types of regulation resources into Equations (1)- (12). Then we can rank the results of the η b f and P e f f ective according to the calculations. In this paper, the regulation resource is electrical ESS and gas turbine. The specific parameters are shown in Table 5. In order to reflect the η b f and P e f f ective of different regulation resources, hydropower and the Demand Side Resource (DSR) datas are added for comparison into Table 5. In order to show the cost difference between electrical ESS of RIES and other types of regulation resources that have participated into regulation service, three operating modes are selected as examples in this paper. Operating Cost of RIES under there three modes are shown in Table 6. The output of each unit from the RIES under Mode 1 is shown in Figure 8. RIES utilizes the electrical ESS to decouple "electric power defined by its quantity of thermal produced" coupling relationship. The power output limit of the gas turbine is removed.
In order to show the cost difference between electrical ESS of RIES and other types of regulation resources that have participated into regulation service, three operating modes are selected as examples in this paper.

•
Mode 1: RIES is in economic operation mode and does not participate in regulation service. • Mode 2: Gas turbines and electrical ESS in RIES participate in regulation services. • Mode 3: Only the electrical ESS in RIES participates in regulation services.
Operating Cost of RIES under there three modes are shown in Table 6. The output of each unit from the RIES under Mode 1 is shown in Figure 8. RIES utilizes the electrical ESS to decouple "electric power defined by its quantity of thermal produced" coupling relationship. The power output limit of the gas turbine is removed. The objective function with lowest electricity generation cost is connected to the objective function with lowest heat generation cost. Obtain the optimal dispatching strategy of each unit under the premise of satisfying the operating constraints, while with the minimum cost of electricity/ heat generation. In Mode 1, RIES achieves the electric-thermal joint dispatch, but does not participate into the regulation market.
In Mode 2, when electrical ESS of RIES and gas turbines are used as regulation resources, their ramp rate of them are infinite, which means that t = 0, , = , = 0. And the noneconomic dispatching cost in regulation service = 0. However, due to the gas turbine acts as a regulation unit with the down/up time constraints, it needs to ensure a certain period of operating The objective function with lowest electricity generation cost is connected to the objective function with lowest heat generation cost. Obtain the optimal dispatching strategy of each unit under the premise of satisfying the operating constraints, while with the minimum cost of electricity/heat generation. In Mode 1, RIES achieves the electric-thermal joint dispatch, but does not participate into the regulation market.
In Mode 2, when electrical ESS of RIES and gas turbines are used as regulation resources, their ramp rate of them are infinite, which means that t = 0, c loc,be f ore = c loc,a f ter = 0. And the non-economic dispatching cost in regulation service c loc = 0. However, due to the gas turbine acts as a regulation unit with the down/up time constraints, it needs to ensure a certain period of operating time. Besides, the fuel cost is too high, so the overall operating costs increased.
In Mode 3, only the electrical ESS in RIES participates in regulation services. System operating cost can be significantly reduced. By rational use of the compensation mechanism of electrical ESS to participate in regulation service, and deeply engaged in regulation on the basis of promoting the renewable energy consumption, the fact that electrical ESS participates into regulation services can bring additional revenue can be fully reflected. Simulation result of Mode 3 is shown as Figure 9. The objective function with lowest electricity generation cost is connected to the objective function with lowest heat generation cost. Obtain the optimal dispatching strategy of each unit under the premise of satisfying the operating constraints, while with the minimum cost of electricity/ heat generation. In Mode 1, RIES achieves the electric-thermal joint dispatch, but does not participate into the regulation market.
In Mode 2, when electrical ESS of RIES and gas turbines are used as regulation resources, their ramp rate of them are infinite, which means that t = 0, , = , = 0. And the noneconomic dispatching cost in regulation service = 0. However, due to the gas turbine acts as a regulation unit with the down/up time constraints, it needs to ensure a certain period of operating time. Besides, the fuel cost is too high, so the overall operating costs increased.
In Mode 3, only the electrical ESS in RIES participates in regulation services. System operating cost can be significantly reduced. By rational use of the compensation mechanism of electrical ESS to participate in regulation service, and deeply engaged in regulation on the basis of promoting the renewable energy consumption, the fact that electrical ESS participates into regulation services can bring additional revenue can be fully reflected. Simulation result of Mode 3 is shown as Figure 9.  Combined with the price curve to analyze, the operating characteristics of the electric storage does not participate in regulation service can be discussed as follows: Figure 8a shows the electrical load balance condition and the ES-SOC curve in the economic model. At 1: 00, when the day-ahead LMP is high, RIES releases energy to the power grid to arbitrage and reduce the operating cost. At 1: 00-10: 00, the LMP and load are in the trough period. The renewable energy output gradually increased, the electric load is supplied by the power grid and the gas turbine, and the electric ESS is charged. At 11: 00-17: 00, the load gradually increased and the day-ahead LMP is relatively stable. RIES gains more power from the grid. At 17:00-18:00, the day-ahead LMP reaches the peak. Electrical ESS starts to release energy to reduce system operating cost. At 20:00-21:00, the procedure is the same as 17:00-18:00. Figure 8a shows the thermal load balance condition and the TS-SOC curve in the economic model. Due to the presence of the thermal ESS and electric boiler, thermal load is not completely provided by GTs. The flexibility of system operation has been improved. Figure 9 shows the electrical load balance condition and ES-SOC curve in the regulation model. The operating characteristics of electrical ESS which participates in the regulation service are summarized as follows: the regulation capacity clearing price and the regulation performance clearing price is the highest at 2:00, 11:00, 14:00, 17:00, 22:00. The electrical ESS should be charged before these windows. With the rising of the regulation capacity clearing price and the regulation performance clearing price, the electrical ESS begins to release energy and makes RIES profit from the regulation market.

Conclusions
To increase the proportion of renewable energy in the energy system is not only the future direction but also an important strategy for the development of the global energy industry. However, the great amount of the renewable energy integration into the power system will cause non-negligible frequency offsets. In this paper, an optimization model for RIES to participate in the regulation market based on pay-for-performance is proposed and three operating modes are selected as the example. To achieve the maximum profit by making a balance between the power trading revenue and regulation market revenue, this paper uses the actual operation model and the data of the PJM regulation market to guide the optimal allocation of regulation resource in the regional integrated energy system. The simulation results in MATLAB are presented s to verify the performance of the optimal scheduling. The optimization model of RIES with electrical ESS based on pay-for-performance as proposed in this paper shows promising prospects for future development.