Next Article in Journal
Research and Performance Evaluation on Selective Absorption of H2S from Gas Mixtures by Using Secondary Alkanolamines
Next Article in Special Issue
Optimization Strategy of Hybrid Configuration for Volatility Energy Storage System in ADN
Previous Article in Journal
Green Extraction Process of Food Grade C-phycocyanin: Biological Effects and Metabolic Study in Mice
Previous Article in Special Issue
Layout Method of Met Mast Based on Macro Zoning and Micro Quantitative Siting in a Wind Farm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Community Integrated Energy System Multi-Energy Transaction Decision Considering User Interaction

School of Electrical Engineering, Xinjiang University, Urumqi 830017, China
*
Author to whom correspondence should be addressed.
Processes 2022, 10(9), 1794; https://doi.org/10.3390/pr10091794
Submission received: 16 July 2022 / Revised: 27 August 2022 / Accepted: 28 August 2022 / Published: 6 September 2022
(This article belongs to the Special Issue Modeling, Analysis and Control Processes of New Energy Power Systems)

Abstract

:
With the gradual liberalization of China’s energy market, the distributed characteristics of each entity in the community integrated energy system are more and more obvious, and the traditional centralized optimization is difficult to reveal the interaction between the entities. This paper aims to improve the profit of the community operator and the users’ value-added benefit of energy use, and proposes a multi-energy transaction decision of a community integrated energy system considering user interaction. First, a refined model of user interaction, including energy conversion, is established, and then the optimization model of multi-energy transaction decision between the community operator and the users is constructed based on the master–slave game. The upper layer aims to maximize the profit of the community operator according to the energy use strategies’ feedback from the users, decides the retail energy prices of the community operator to the users, and optimization variables include equipment output and energy purchased from the power grid and natural gas grid. The lower layer aims to maximize the value-added benefit of energy use for users. The users optimize their energy use strategies based on the retail energy prices published by the community operator. The model is solved by the differential evolution algorithm combined with the CPLEX solver. Finally, different scenarios are analyzed in a numerical example, and the results show that the strategy proposed in this paper to set community prices increases the community operator’s profit and profit margin by 5.9% and 7.5%, respectively, compared to using market energy prices directly. At the same time, the value-added benefit to users also increases by 15.2%. The community operator and users can achieve a win–win situation.

1. Introduction

With the disadvantages of low economic benefits and high energy consumption of traditional energy systems becoming increasingly prominent, integrated energy systems (IES) that can realize flexible energy conversion and efficient utilization have become the focus of energy research and development [1]. The community integrated energy system (CIES) near the user side contains a variety of energy coupling equipment, which couples and complements electricity, heat, natural gas, and other energy sources, enabling local consumption of renewable energy and providing users with comprehensive energy services. It is an important direction for the future development of the intelligent community [2]. Therefore, it has become a hot research issue to study how to improve the economics of the community integrated energy system, and to formulate transaction strategies between the community operator and the users to guide the users to rationally use energy to achieve a win–win situation [3,4].
At present, domestic and foreign scholars have focused on improving the economics of the integrated energy system, mainly on the refined modeling of equipment on the power supply side and demand side management. In terms of refined equipment modeling, this technology proposes a general dynamic energy efficiency model of an integrated energy system, which lays the foundation for the optimization of integrated energy system operation and the formulation of trading strategies [5]. Chen et al. established an optimization model considering coupled dynamic energy efficiency, and the results showed that considering the dynamic energy efficiency of equipment can improve the energy utilization rate, which is more consistent with the actual optimization results [6]. In demand side management, Wang Yongli et al. considered the electric and heating demand response, and the established source–load interactive model can reduce the operating cost of the service provider and improve wind power consumption [7]. Guo Zihao et al. considered the multi-energy flow coupling characteristics and the users’ flexible load and considered the source–load interaction to optimize the operator’s benefit [8]. Liu et al. considered the controllable degree of flexible load in the scheduling process and the constraint of user satisfaction to give full play to the demand response ability of the users. The results show that the proposed scheduling scheme can reduce the cost of the community operator [9]. It can be seen that coordinated optimization on both sides of supply and demand can improve the system economy. However, the references [7,8,9] are mostly centralized optimization and do not consider the energy trading and pricing problems of the operator in the market environment. There is a lack of research on the impact of the operator’s retail energy prices on users’ energy-use strategies.
With the development of the electricity market, the community operator can be regarded as a distribution-side entity or retail entity with self-production and self-selling capabilities [10]. It can guide users to participate in interaction by formulating reasonable retail energy prices, adjusting energy use strategies, and achieving demand side management. Aiming at the transaction problem between the operator and the users, Li Yuan et al. used the master–slave game method to establish an operator energy pricing model, including electric vehicles and P2G, which can improve the system economy [11]. References [12,13] analyzed the interaction mechanism between the community operator and the users in the electricity market based on the master–slave game model, with the operator as the leader and the users as the follower. In reference [14], the master–slave game model of the transaction between the community operator and consumers was established, and transaction strategies considering the demand response ability of consumers were proposed. Fu et al. constructed a user model containing four types of loads: electricity, heat, cold, and gas. Combined with the operator revenue optimization model, a master–slave game pricing mechanism between operators and users is proposed [15]. Fleischhacker et al. proposed an energy value allocation and stabilization algorithm based on a cooperative game. By investing in distributed energy, community operators can share value among their members [16]. Based on the master–slave game, Anoh et al. constructed an energy trading strategy between operators and consumers in the microgrid to optimize the interests of producers and consumers [17]. Wei et al. proposed a multi-leader and multi-follower Stackelberg game approach to solve the multi-energy trading problem. Multiple energy operators act as leaders to determine real-time energy prices, while multiple consumers act as followers to optimize their energy usage strategies [18]. However, the above models do not consider the role of convertible load in the process of user interaction. Li Peng et al. included convertible loads in consideration of integrated demand response, which improved user interaction but did not consider the impact of energy prices on convertible load [19]. Under the incentive of multiple retail energy prices, users will preferentially use energy with lower prices to meet the same energy demand. The actual amount of interaction will be influenced by consumer psychology [20], so considering the convertible load can further tap into the potential of user interaction.
Based on the above research, this paper establishes a refined model of user interaction considering energy conversion and constructs an optimization model of multi-energy transaction decisions between the community operator and the users based on the master–slave game. Taking the community operator as the leader, the optimization goal is the maximum daily profit, and the optimization variables are retail energy prices, equipment output, etc. The users are the followers, the optimization goal is the maximum value-added benefit of energy use, and the optimization variables are the users’ energy use strategies of electricity, heat, and natural gas. Finally, the validity of the proposed model is verified by an example.
The rest of this paper is organized as follows. Section 2 summarizes the models of community integrated energy system and of user interaction. Section 3 establishes the optimization model for the multi-energy transaction decision between the community operator and the users. Section 4 sets up different scenarios to analyze the trading strategy proposed in this paper. Finally, the conclusions are drawn in Section 5.

2. Models of Community Integrated Energy Systems and User Interactions

2.1. CIES Model Based on Energy Hub

The concept of an energy hub (EH) was first proposed by Geidl et al. [21], it simplifies the energy flow relationship. An EH with multiple input and output ports is modeled by a coupling matrix that can easily describe the transformation and coupling relationship between energy input and output [22]. Therefore, to analyze the energy coupling and input–output energy flow relationship in the system, the energy hub model is used to describe the CIES model abstractly, as shown in Figure 1. In this paper, the electricity, heat, and natural gas demanded by the users in winter are supplied by the community operator in CIES, who has certain renewable energy units according to natural and geographical conditions. In the actual operation process, the community operator purchases electric energy and natural gas energy from the energy market and uses the energy conversion equipment to convert the energy into the energy required by the users according to the multi-energy complementary characteristics. The renewable energy equipment of CIES includes wind turbines (WT) and photovoltaic (PV); energy conversion equipment includes combined heat and power (CHP) units, gas boilers (GB), and electric heat pumps (EHP); energy storage devices include electricity storage and heat storage.
According to the energy flow, the community energy supply model can be represented by the following matrix:
[ P out , t e P out , t h P out , t g ] = [ 1 r 1 0 r 2 η CHP , e r 1 η EHP 1 r 2 η CHP , h + r 3 η GB 0 0 1 r 2 r 3 ] [ P in , t e P in , t h P in , t g ] [ P c / d , e , t ES , n P c / d , h , t ES , n 0 ]
In the above formula: P out , t e , P out , t h , and P out , t g are the electricity, heat, and natural gas power supplied by the community to the users, respectively; P in , t e is the sum of purchased power P net , t e , wind power P w , t e , and photovoltaic power P v , t e ; P in , t h is the heat power purchased by the community. This paper considers that the heat energy is only supplied by the community, so it is taken as 0; P in , t g is the power of the community to purchase natural gas; P c / d , e , t ES , n and P c / d , h , t ES , n are the actual charging and discharging power of electric and thermal energy storage respectively; r 1 , r 2 , and r 3 are dispatch factors, which represent the proportion of heterogeneous energy flow into energy conversion equipment; η EHP and η G B are the efficiency of electric heat pump and gas boiler; η CHP , e and η CHP , h are the electrical and thermal efficiencies of the CHP unit.

2.2. The Models of User Interaction That Account for Energy Conversion

User interaction in energy communities makes sense and can improve the economy of the community operator [8,11]. In this paper, the community operator guides users to adjust energy use strategies by formulating reasonable retail energy prices. In general, user interaction strategies are often limited to responses in a single form of energy, such as load reduction and transfer, which have a greater impact on users’ actual energy use. With the development of user terminal equipment, the user side can realize the conversion of energy forms. When the community operator publishes retail energy prices, the multi-energy complementary users consider the energy conversion efficiency of the terminal equipment to obtain the difference in equivalent energy prices. Users can use the corresponding terminal equipment to achieve secondary energy conversion and choose the appropriate way to meet their own load demand. For example, users can choose to use electric heating or natural gas heating to achieve the same hot water demand according to the equivalent electricity price and the equivalent natural gas price, etc. At the same time point, the users’ actual energy use demands do not change, and the users’ actual energy use has little impact, which can improve the flexibility and economy of the users’ energy use. Therefore, it is of great significance to construct user-interactive models that consider the conversion of energy use.
According to the above analysis, multi-energy complementary users can choose to convert, reduce, transfer, and other ways to achieve interaction.

2.2.1. Convertible Load Model

This paper considers the users’ electricity–gas convertible load and improves the convertible model of reference [23]. That is, the influence of consumer psychology is considered when optimizing the convertible load model. Based on the principle of consumer psychology, the difference between the equivalent electricity price and the equivalent natural gas price affects the response of the user’s convertible load. The users’ interactive response range is divided into saturation zone, linear zone, and dead zone [24,25]. When the difference between the equivalent electricity price and the equivalent natural gas price is lower than the dead zone threshold, users are unwilling to respond interactively. When it exceeds the threshold, users start to respond interactively. In the linear zone, the users’ convertible load increases with the difference between the equivalent electricity price and the equivalent natural gas price, and it shows a linear upward trend. When the compensation limit is exceeded, the users’ electricity–gas convertible amount tends to be saturated. In this paper, the energy use rate λ con , t g , e of the users’ gas load to electric load and the energy use rate λ con , t e , g of the users’ electric load to gas load are used to characterize the influence of the difference between the equivalent electricity price and the equivalent natural gas price π con , t on the users’ mode of energy use.
λ con , t g , e = { 1 π con , t π con , max g , e π con , t + π con , min g , e π con , min g , e π con , max g , e π con , max g , e < π con , t < π con , min g , e 0 π con , min g , e π con , t 0
λ con , t e , g = { 0 0 π con , t π con , min e , g   π con , t π con , min e , g π con , max e , g π con , min e , g π con , min e , g < π con , t < π con , max e , g 1   π con , max e , g π con , t
In the above formula: π con , min g , e , π con , min e , g , π con , max g , e , and π con , max e , g are the dead zone threshold and saturation zone limit of the difference between the equivalent electricity price and the equivalent natural gas price when the users respond to convertible load.
According to the calorific value equivalence theorem and the energy conservation theorem, the constraints that the convertible load needs to satisfy are as follows:
{ L con , t e = L con , t e , n μ con , t e , g Δ L con , t e + μ con , t g , e Δ L con , t g , e L con , t g = L con , t g , n μ con , t g , e Δ L con , t g + μ con , t e , g Δ L con , t e , g Δ L con , t e = λ con , t e , g L con , t e , n Δ L con , t g = λ con , t g , e L con , t g , n Δ L con , t g , e = Δ L con , t g / I , Δ L con , t e , g = I Δ L con , t e μ con , t e , g + μ con , t g , e 1
In the above formula: L con , t e , n and L con , t g , n are the power before the convertible electrical load and convertible natural gas load response; L con , t e and L con , t g are the power after the response of the convertible electrical load and convertible natural gas load; μ con , t e , g and μ con , t g , e are 0–1 auxiliary variables for interactive response; Δ L con , t e and Δ L con , t g are convertible electrical load and convertible natural gas load response quantities; Δ L con , t g , e and Δ L con , t e , g are the increased power of the convertible electrical load and convertible natural gas load after the response; I is the electricity–gas conversion coefficient, which is taken as 1.25 in this paper.

2.2.2. Reducible Electrical Load Model

The reducible electrical load is the load that users can partially reduce [7]. The model is as follows:
{ L adj , t e = L adj , t e , n μ adj , t e Δ L adj , t e 0 Δ L adj , t e L adj , t e , n
In the above formula: L adj , t e , n is the load power that can be reduced before the users respond; μ adj , t e is a 0–1 variable of whether to reduce; Δ L adj , t e is the load power actually reduced by the users.
Considering that power load reduction has a great impact on user satisfaction, the maximum duration of power load reduction is constrained in this paper:
τ = t t + t adj , max e ( 1 μ adj , τ e ) 1 t = 1 , 2 , , T t adj , max e
In the above formula: t adj , max e is the maximum duration of electrical load reduction.

2.2.3. Transferable Electrical Load Model

After the community publishes the electricity price, users will transfer part of the electrical load from higher to lower hours to reduce the cost of energy use, such as washing machines, electric vehicles, and other loads. This part of the load is called a transferable power load [7]. The transferable electrical load model is as follows:
{ L tran , t e = L tran , t e , n + L tran , t e , in L tran , t e , out t = 1 T L tran , t e , in = t = 1 T L tran , t e , out 0 L tran , t e , in L tran , max e , in 0 L tran , t e , out L tran , max e , out
In the above formula: L tran , t e , n and L tran , t e are the power before and after the response of the transferable electrical load at time t; L tran , t e , in and L tran , t e , out are the actual transfer-in and transfer-out power of the transferable electrical load at time t; L tran , max e , in and L tran , max e , out are the maximum load that can be transferred in and out at time t.

2.2.4. Heat Load Model

Some users in the community have high requirements for thermal comfort, the indoor temperature cannot be adjusted, and they are willing to bear the additional cost of thermal comfort. This part of the heat load is a fixed heat load. Another part of the user is willing to adjust the thermal comfort range, which is an adjustable heat load.
The adjustable heat load adopts the first-order building thermodynamic model [26]:
L adj , t h = N 1 1 R ( T in , t + 1 e Δ t / τ 1 T in , t 1 e Δ t / τ 1 T out , t )
In the above formula: N 1 is the number of users with adjustable heating temperature; τ 1 = R C air , C air is the heat capacity of indoor air, which can be taken as 1.2 kWh/°C, R is the equivalent thermal resistance of the house, which can be taken as 6.8 °C/kW; T in , t is the indoor temperature of the heating that can be adjusted at time t; T out , t is the outdoor temperature.
The interactive response potential of heat load is mainly related to the human body’s heat-using psychology for temperature perception, which has a certain elasticity. To better describe the user’s thermal response potential, this paper introduces predicted mean vote (PMV), and the relationship between room temperature and PMV index value λ PMV , t is as follows [27]:
T in , t = T com , s M ( 2.43 λ PMV , t ) ( λ clo + 0.1 ) 3.76
In the above formula: T com , s is the average temperature of human skin in a comfortable state, which can be taken as 33.5 °C; λ clo and M are the thermal resistance of the clothes and the metabolic rate of the human body, take 0.11 (m2·°C)/W and 80 W/m2 respectively.
Considering the recommendations of the ISO-7730 standard and the daily routine of the users, this paper limits the time sharing of the PMV index value, which is expressed as:
{ | λ PMV , t | 1 , t [ 1 , 7 ] [ 21 , 24 ] | λ PMV , t | 0.5 , t [ 8 , 20 ]
The fixed heat load model is:
L fir , t h = N 2 1 R ( T set e Δ t / τ 1 T set 1 e Δ t / τ 1 T out , t )
In the above formula: N 2 is the number of users with non-adjustable heating temperature; T set is the most comfortable indoor temperature of the users who cannot be adjusted.
To prevent the indoor temperature of adjustable users from being lower than the most comfortable temperature, the following constraints are imposed on the indoor temperature:
t = 1 T T in , t T = T set
To sum up, the refined model of interaction considering user-side energy conversion is described as a matrix as follows:
[ L t e L t h L t g ] = [ P out , t e P out , t h P out , t g ] = [ L fir , t e L fir , t h L fir , t g ] + [ L adj , t e L adj , t h 0 ] + [ L tran , t e 0 0 ] + [ L con , t e 0 L con , t g ]
In the above formula: L t e , L t h and L t g are the electricity load, heat load, and natural gas load of the users; L fir , t e and L fir , t g are fixed electrical load and natural gas load; L fir , t h is the fixed heat load.

3. Optimization Model of Multi-Energy Transaction Decision between the Community Operator and the Users

3.1. Model Architecture of the Multi-Energy Transaction between the Community Operator and the Users Based on the Master–Slave Game

The master–slave game is an effective method to solve the problem of how to make decisions when there is an interest relation or conflict. The master–slave game is a dynamic non-cooperative game, and the unequal status of participants is the most fundamental difference between the master–slave game and the classical game. In the master–slave game, each subject has a different status and decision-making sequence. The leader has a leadership advantage and can occupy the first or advantageous position in the game, and the follower must follow the leader to make decisions. Not only does the retail price of energy set by the community operator affect consumer demand, but demand also affects price. Both parties have independent interest demands, and both make decisions with the goal of maximizing their own interests and influencing each other, and there is a master–slave game relationship.
The interaction relationship between the community operator and the users based on the master–slave game is shown in Figure 2. In the energy market environment, the community operator guides users to interactively use energy by setting reasonable retail energy prices. The energy purchase strategies of the energy market and equipment output are optimized under the maximum profit, and the users adjust the interactive energy use strategies according to the energy sales price of the community operator. In the master–slave game, the community operator, as the leader, guides users to interactively use energy by adjusting the retail prices of electricity and heat energy. As followers, users who receive retail prices of electricity and heat energy released by the community operator will change their energy usage habits to a certain extent and adjust their interactive energy usage strategies. At the same time, the users’ changes in their own energy use needs will affect the retail energy price formulation strategies of the community operator. Additionally, the cycle goes back and forth when neither the community operator nor the users can improve their own interests by changing their own strategies; the equilibrium solution of the master–slave game model is reached. Finally, the community operator determines the final retail energy prices, the energy purchased in the energy market, and the output of each device. The users determine the interactive energy use strategies according to the retail energy prices released by the community operator.

3.2. Pricing Model of Community Operator’s Retail Energy Prices

3.2.1. Objective Function

The community operator guides users to interact, optimize energy transactions, and maximize daily profit by formulating reasonable retail energy prices in the upper-level model. The objective function is:
C 1 = max ( C sell + C add C ob C mend C wv )
In the above formula: C 1 is the daily profit of the community operator; C sell is the energy retail income; C add is the additional income for thermal comfort; C ob is the cost of purchasing energy; C mend is the cost of equipment operation and maintenance; C wv is the cost of renewable power abandonment.
(1)
Energy retail income
C sell = t = 1 T i E F i , t L t i
In the above formula: F i , t is retail energy prices of the community; E is the load type of the users, E = { e , h , g } .
(2)
Additional income for thermal comfort
C add = t = 1 T F add L fir , t h
In the above formula: F add is the additional unit service price at which the heating load is not adjustable for the users to maintain the most comfortable temperature.
(3)
Cost of purchasing energy
C ob = t = 1 T [ F net , e , t P net , t e + F net , g , t P in , t g ]
In the above formula: F net , e , t and F net , g , t are the unit cost of purchasing electricity and natural gas at time t.
(4)
Cost of equipment operation and maintenance
C mend = t = 1 T j = 1 J F m , j P j , t
In the above formula: F m , j is the unit operation and maintenance costs of equipment j; P j , t is the output of equipment j at time t.
(5)
Cost of renewable power abandonment
C wv = t = 1 T [ Δ P w , t F w + Δ P v , t F v ]
In the above formula: Δ P w , t and Δ P v , t are the amount of abandoned wind and light at time t; F w and F v are the unit cost of abandoned wind and light.

3.2.2. Multiple Energy Price Constraints

To prevent the community operator from maliciously raising prices and ensuring the benefit of users, multiple energy prices need to be constrained:
F i , t , min F i , t F i , t , max
t = 1 T F i , t T F i , av
In the above formula: F i , t , max and F i , t , min are the maximum and minimum prices for electricity and heat prices to protect the benefits of the users and the operator; F i , av is the average price of electricity and heat in the energy market.

3.2.3. Energy Conversion Equipment Constraints

The traditional energy conversion equipment model considers the output efficiency of the equipment to be a fixed value. To ensure the accuracy of the equipment output, this paper adopts the dynamic energy-efficiency model of energy conversion equipment. The output efficiency of energy conversion equipment is mainly related to the load rate [28,29]. It is expressed by means of polynomial fitting, as shown in Formula (22).
η x , n = k = 0 n α x , k ( P x P x , N ) k
In the above formula: η x , n is the dynamic energy efficiency of equipment x using n-order polynomial fitting; α x , k is the fitting coefficient; P x and P x , N are the actual output power and rated output power of the equipment.
The input–output relationship of energy conversion equipment considering dynamic energy efficiency is as follows:
(1)
Combined heat and power unit
The output power efficiency of the CHP unit can be fitted by a fourth-order fitting [28]:
P CHP , t e = η CHP , 4 P CHP , t g
In the above formula: P CHP , t g is the natural gas power entering the CHP unit; P CHP , t e is the output electric power of the CHP unit.
In this paper, the strategy of determining the heat by electricity is adopted, and the thermoelectric ratio ψ CHP can be described by the second-order fitting of the electrical load rate of the CHP unit:
ψ CHP = k = 0 2 ( α ψ , k N CHP , e k )
In the above formula: N CHP , e is the electrical load rate of the CHP unit.
Therefore, the output thermal power P CHP , t h of the CHP unit is:
P CHP , t h = ψ CHP P CHP , t e
(2)
Gas boiler
The output thermal efficiency of the gas boiler can be fitted by the first-order:
P GB , t h = η GB , 1 P GB , t g
In the above formula: P GB , t g is the natural gas power entering the gas boiler; P GB , t h is the output thermal power of the gas boiler.
(3)
Electric heat pump
The output thermal efficiency of the electric heat pump is related to the load rate and temperature. In this paper, only the influence of the load rate is considered, and second-order fitting can be used [29]:
P EHP , t h = η EHP , 2 P EHP , t e
In the above formula: P EHP , t e is the electric power entering the electric heat pump; P EHP , t h is the output heat power of the electric heat pump.

3.2.4. Device Operation Constraints

(1)
Energy Conversion Equipment Constraints
{ P i min P i , t P i max P i down P i , t P i , t 1 P i up
In the above formula: P i max and P i min are the upper and lower limits of the output of the coupling device at time t; P i up and P i down are the upper and lower limits of the climbing power of the coupling device at time t.
(2)
Energy Storage Device Constraints
{ S i , t ES = ( 1 σ i ES ) S i , t 1 ES + P c / d , i , t ES , n Δ t P c / d , i , t ES , n = η c , i ES P c , i , t ES 1 η d , i ES P d , i , t ES C i , min ES S i , t ES C i , max ES 0 P c , i , t ES ε c , i , t P c , i , t max 0 P d , i , t ES ε d , i , t P d , i , t max S i , 0 ES = S i , T ES
In the above formula: S i , t ES is the energy storage value of the energy storage device at time t; σ i ES is the self-loss rate of the energy storage device; η c , i ES and η d , i ES are the charging and discharging efficiency of the energy storage device; P c , i , t ES and P d , i , t ES are the charging and discharging power of the energy storage device; C i , max ES and C i , min ES are the upper and lower limits of the capacity of the energy storage device; ε c , i , t and ε d , i , t are 0–1 auxiliary variables; P c , i , t max and P d , i , t max are the upper limit of energy storage charging and discharging power respectively.

3.2.5. Renewable Energy Output Constraints

In the below formula: P w , t e , pre and P v , t e , pre are the predicted power values of wind power and photovoltaic.
{ 0 P w , t e P w , t e , pre 0 P v , t e P v , t e , pre
In addition, the upper-layer model also needs to satisfy the power balance constraint of Equation (13). The above nonlinear model can be approximated by piecewise linearization.

3.3. Energy Use Strategies Model Considering User Interaction

According to the retail energy prices released by the community operator, users consider increasing or decreasing load, converting load, and transferring load. For details, see the refined model of user interaction considering energy conversion in Section 2. The users’ goal is to optimize their own interactive energy-use strategies and maximize the value-added benefit of energy use. This paper defines the value-added benefit of a users’ energy use, which is the users’ total energy use utility C ben minus the total cost. The total cost includes the cost of purchasing energy C ub and the additional cost of thermal comfort C app .
C 2 = max [ C ben ( C ub + C app ) ]
C ben = t = 1 T i E [ f 1 i L t i f 2 i 2 ( L t i ) 2 ]
C ub = C sell = t = 1 T i E F i , t L t i
C app = C add = t = 1 T F add L fir , t h
In the above formula: f 1 i and f 2 i are the constant coefficients of the users’ preference for type i energy, reflecting the users’ preference for energy demand.

3.4. Model-Solving Process

According to the reference [30], it is easy to prove that there is a unique equilibrium solution for the optimization model of multi-energy transaction decisions between the community operator and the users based on the master–slave game in this paper.
This paper uses MATLAB software programming, and the method of combining the differential evolution algorithm and CPLEX solver is used to solve the proposed model. The solution process is as follows:
(1)
Initialize the parameters of the community operator and the users, k = 0 , set the maximum number of iterations k max = 100 , use the differential evolution algorithm to randomly generate retail energy prices of 10 groups of the community operator, and transmit them to the energy use strategies model considering user interaction.
(2)
k = k + 1 .
(3)
The users receive retail energy prices published by the community operator. Use the CPLEX solver to solve the energy use strategies model and the optimal value-added benefit C 2 k , and return the energy use strategies to the model of the community operator.
(4)
The community operator optimizes the output of equipment and the amount of electricity and gas purchased in the market according to the energy use strategies of the users and calculates the optimal profit C 1 k of the community operator.
(5)
Use the variation and crossover of the differential evolution algorithm to generate a group of new retail energy prices and repeat the processes in (3) and (4). Additionally, calculate the optimal value-added benefit C 2 k of the users and the optimal profit C 1 k of the community operator under the new retail energy prices.
(6)
Perform selection operation: compare the optimal solutions of the community operator before and after mutation and crossover; if C 1 k C 1 k , then C 1 k + 1 = C 1 k , C 2 k + 1 = C 2 k ; if C 1 k < C 1 k , then C 1 k + 1 = C 1 k , C 2 k + 1 = C 2 k .
(6)
If k k max , end the program; otherwise, return to flow (2).

4. Case Analysis

4.1. Parameter Settings

This paper selects a community integrated energy system in the northern winter of China as the research object. The forecasting curves of the renewable energy output, initial load, and outdoor temperature curves are shown in Figure 3. Other required parameter data are shown in the table in the Appendix A. The equipment parameters of the community operator are shown in Table A1, which is used to deal with the community operator’s device model. The data parameters for the users are shown in Table A2, which is the parameter data required by the user model. The equipment efficiency fitted from Table A2 is shown in Figure A1. The time-of-use electricity price in the energy market is shown in Table A3 [31], it is the price that the community operator trades with the energy market. To protect the benefit of the users, the upper limit of the electricity price set by the community operator is not higher than the time-of-use electricity price in the energy market; the lower limit is not less than 0.2 RMB/(kWh). The heat price of the energy market is 0.35 RMB/(kWh) [32], and the thermal price range set by the community ranges from 0.2 to 0.5 RMB/(kWh). The natural gas price within the community is the same as the market gas price, which is 0.34 RMB/(kWh) [33].

4.2. Scene Settings

To illustrate that the multi-energy trading strategies proposed in this paper can improve the profit of the community operator and the value-added benefit of the users, and can improve renewable energy utilization, three scenarios are set for comparative analysis:
  • Scenario 1: Regardless of user interaction [6], the energy price sold by the community operator to the users is the market price.
  • Scenario 2: Considering user interaction [15] and ignoring the influence of the users’ electricity–gas convertible load, the community operator and the users compete to determine the electricity price and heat price of the community.
  • Scenario 3: Considering user interaction, using the refined model of user interaction considering energy conversion and considering the impact of electricity–gas convertible load, the community operator and the users play games to formulate the internal electricity price and heat price in the community.

4.3. Simulation Analysis

4.3.1. Analysis of the Multi-Energy Transaction Results of the Community Operator

Under the model proposed in this paper in Scenario 3, the electricity and heat prices traded between the community operator and the users are shown in Figure 4.
The electricity price trend set by the community operator basically follows the electricity price in the energy market, and the peak electricity price is in the peak energy use period, which is in line with the actual situation. The heat price set by the community operator fluctuates within the limit of the heat price. Judging from the electricity and heat prices traded by the community operator and the users in Figure 4, the average price of electricity set by the community operator in Scenario 3 is 0.5 RMB/(kWh), the average price of electricity in the energy market is 0.56 RMB/(kWh). The average price of heat within the community is 0.33 RMB/(kWh), and the average price of heat in the energy market is 0.35 RMB/(kWh). Based on the above analysis, the average electricity price and average heat price set by the community operator are 10.7% and 5.7% lower, respectively, than in the market. In Figure 4, in Scenario 3, the electricity and heat prices traded between the community operator and the users can protect the benefit of the users and promote energy transactions between the users and the community operator.
Figure 5 shows the results of electricity and gas transactions between the community operator and the energy market in different scenarios. Taking the peak period of electricity price as an example, the following analysis is made: from Figure 5, it can be seen that in Scenarios 2 and 3, compared with Scenario 1, the community operator purchases less electricity from the power grid during the peak period of electricity price, which effectively reduces the power supply pressure on the large power grid, indicating that user interaction can indirectly participate in the power market and reduce the peak power consumption of the large power grid. Compared with Scenario 2, Scenario 3 considers the convertible load, and part of the electricity load demand of the users is supplied by natural gas during the peak electricity price period. The energy market has the least amount of electricity purchased and the largest amount of gas during the corresponding period.

4.3.2. Analysis of Energy Use Strategies for User Interaction

Figure 6 shows the results of electricity and natural gas transactions between the community operator and the users in different scenarios. Scenario 1 does not consider user interaction, and its electricity and natural gas loads are all original loads. As shown in Figure 6, in both Scenarios 2 and 3, considering user interaction, the electrical load peak smoothed to varying degrees. Scenario 1 has the highest electrical load peak, followed by Scenario 2. Scenario 3 adopts the strategies of this paper, with the smallest electrical load peak and the largest gas load fluctuation.
Considering the limitations of space, this paper will focus on the analysis of the interactive energy use strategies of users’ electricity–gas convertible load. In Scenario 3, the amount of the convertible load interaction, including consumer psychology, is shown in Figure 7, which shows the relationship between the users’ actual convertible load response and the equivalent energy price difference. This is consistent with the optimized results in Figure 6. For users, when the electricity price is lower than the equivalent natural gas price (the natural gas price multiplied by the electricity–gas conversion coefficient). For example, to meet the same demand for hot water, the electricity cost of the users is lower than the natural gas cost. When the natural gas price is high, the cost of electricity is higher than the cost of natural gas. Therefore, when the electricity price is lower than the equivalent natural gas price, users replace part of the convertible natural gas load with electricity, and when the electricity price is higher than the equivalent natural gas price, users replace part of the convertible electricity load with natural gas. Compared with Scenario 2, the convertible load is considered in Scenario 3, and the electricity load curve during the peak period of the electricity price during the periods of 11:00~14:00 and 18:00~21:00 is lower than that of Scenario 2. In the low electricity price period from 1:00 to 7:00, the power load curve of Scenario 3 is higher than that of Scenario 2. After considering the convertible load, the electrical load curve of the users in Scenario 3 is smoother, and the outline of the electric load curve is optimized.
When the users’ heat load does not participate in the interaction, that is, the original heat load (Scenario 1), the indoor temperature is kept at the most comfortable temperature. When considering the users’ adjustable heat load, the indoor temperature is kept in a suitable range (Scenario 3 is used as an example). The heat load and indoor temperature of Scenarios 1 and 3 are shown in Figure 8.
According to Figure 8, the indoor temperature of users in Scenarios 1 and 3 cannot be adjusted at room temperature and is always maintained at about 22.6 degrees Celsius. The indoor temperature of users in Scenario 3 can be adjusted to maintain the indoor temperature within 20.5 to 24.5 degrees Celsius, ensuring the users’ thermal comfort. From time period 1 to time period 9, the heat load provided by Scenario 3 is lower than the original heat load, and coupled with heat loss, the temperature decreases. At the same time, to keep the temperature of the first and last sections consistent, the heat load adjustment in the last few periods is relatively large.

4.3.3. Cost-Benefit Analysis of the Community Operator and Users

To show the superiority of the model and trading strategy proposed in this paper, the profit and profit rate of the community operator, the cost of the users, and the value-added benefit are compared in different scenarios. The comparison results are shown in Table 1.
Overall, it can be seen from Table 1 that, compared with Scenarios 1 and 2, Scenario 3 is the establishment of the multi-energy prices within the community under the model proposed in this paper, the profit and profit margin of the operator, the value-added benefit of the users’ energy use, and the renewable energy utilization have increased significantly.
The specific analysis is as follows: In terms of on-site consumption of renewable energy, the on-site renewable energy utilization in Scenario 2 and Scenario 3 has increased by more than 5% compared with Scenario 1. Scenario 3 has the highest renewable energy utilization, followed by Scenario 2 and Scenario 1. Therefore, user interaction can improve on-site renewable energy utilization. In terms of the community operator’s profit and the energy value-added benefit of the users, compared with Scenario 1, Scenario 2 reduces the community operator’s profit by 130 RMB. The cost of the users’ energy use has been reduced by 1584 RMB, and the value-added benefit of the users’ energy use has been greatly improved. This is because Scenario 1 is sold to the users at the market energy prices, the retail energy prices are relatively high, and the users are not given preferential energy prices. Scenario 2 considers user interaction and gaming, the retail price of energy in the community is determined, and the energy prices are constrained within an appropriate range, which can greatly improve user satisfaction with energy use. Compared with Scenario 2, the community operator’s profit in Scenario 3 has increased by 391 RMB. After converting the load, the operator reduces the penalty cost of renewable energy curtailment and the power supply cost during peak power consumption. For the users, the total cost is reduced by 125 RMB after considering the convertible load. At the same time, the value-added benefit of the users’ energy use in Scenario 3 compared with Scenario 2 has been improved, and the satisfaction with energy use has been further improved. Compared with Scenarios 1 and 2, Scenario 3 has the highest profit for the community operator, the lowest total energy use cost for the users, and the highest value-added benefit. Although the community operator in Scenario 3 has the lowest turnover, it has the lowest cost. From the profit side, the profit in Scenario 3 is 5.9% and 9% higher than that in Scenario 1 and Scenario 2, respectively. In terms of profit margin, Scenario 3 has the highest profit margin, which is 7.5% and 4.7% higher than Scenario 1 and Scenario 2, respectively. Meanwhile, the value-added benefit of the users in Scenario 3 is 15.2% and 2.5% higher than that in Scenario 1 and Scenario 2, respectively. Scenario 3 can maximize the benefits of the community operator and the users. The community operator can actively guide user interaction during pricing by optimizing retail energy prices, improving the users’ value-added benefit of energy use, and at the same time increasing the community operator’s own profit and making the users more satisfied with the services provided by the community operator to consolidate and expand the user base and provide other value-added services. To sum up, Scenario 3 can consider the benefits of the community operator and the users and achieve a win–win situation.
To further clarify the work of this paper, Scenario 4 is added: a centralized optimization method is adopted to optimize the maximum profit of a community operator with a single objective [9], considering user interaction, but not considering the convertible load. The comparison results are shown in Table 2.
As can be seen from Table 2, Scenario 4 adopts the centralized optimization with the single goal of maximizing the profit of the community operator, without considering the interests of the users. From the optimization results, the profit of the community operator in Scenario 4 is higher than that in Scenario 3 because the community operator does not offer preferential energy prices to the users and fails to balance the interests of users in Scenario 4. Meanwhile, the cost of the community operator in Scenario 3 is lower than that in Scenario 4, and the profit margin of the community operator in Scenario 3 is 1.8% higher than that in Scenario 4. In Scenario 3, the community operator can optimize the energy prices and balance the value-added benefit of the users through appropriate profit sharing. Although the profit is reduced, it can improve the users’ satisfaction with the use of energy and improve the profitability of the community operator. In terms of user interaction, Scenario 3 considers the convertible load, which can enrich the way users use energy. Under the optimization method of the master–slave game in Scenario 3, the community energy prices can be optimized, and the win–win situation between the community operator and the users can be realized.
Through the above comparative analysis, using the multi-energy transaction decision optimization model of the community operator and the users considering user interaction proposed in this paper, it is possible to formulate reasonable retail energy prices for the community, guide the users to interact, and improve the community operator’s profit and value-added benefit of the users’ energy use to achieve a win–win situation.

5. Conclusions

The master–slave game model constructed in this paper describes the energy transaction between the community operator and the users and proposes the optimization model of the multi-energy transaction decision between the community operator and the users. The upper model considers the maximum profit of the community operator, while the lower model aims for the maximum value-added benefit of the users. The model is effectively analyzed by an example, and the relevant conclusions are as follows:
(1)
The retail energy prices of the community determined by the decision in this paper are reasonable and acceptable to the users. The average price of electricity and the average price of heat set by the community operator are 10.7% and 5.7% lower, respectively, than the market, which protects the interests of the users. The model is extensible. By modifying the corresponding model, more user groups can be promoted. For different countries and regions, the model of the upper community operator can be modified according to the actual situation, including the type of equipment and the type of renewable energy, which can be adjusted according to the needs. The lower user model can determine the types of user interaction load (including reducible, transferable, transferable, etc.) according to the living habits of residents in different countries and regions. At the same time, the lower model has variable data related to the users and can be applied to residential communities with different energy preferences.)
(2)
With the continuous improvement of the user side equipment, convertible load becomes possible. Users can choose appropriate energy modes to meet their energy needs according to different energy prices. The refined user interaction model that considers energy conversion constructed in this paper can reduce user costs.
(3)
The optimization model of the multi-energy transaction decision between the community operator and the users proposed in this paper considers the energy conversion on the user side, which can not only improve the profit of the community operator, but also increase the value-added benefit of energy use and realize a win–win situation for the community operator and the users. Using the strategy proposed in this paper to set the community prices increases the community operator’s profit and profit margin by 5.9% and 7.5%, respectively, compared to using market energy prices directly. At the same time, the value-added benefit to users also increases by 15.2%. In addition, user interaction can indirectly reduce the peak value of the grid, which is beneficial to grid security.
The model established in this paper mainly formulates the retail prices of community energy from the dimensions of the community and the users, ignoring the energy connection between the community and the community. At the same time, the impact of load and renewable energy uncertainty is ignored. The next step will continue to study the energy transaction strategies between the community and the community, and the impact of uncertainty. In terms of the community and the community, a single community may have an energy surplus or a shortage at some time. Multi-communities can trade surplus or shortage energy according to a certain mode to realize the efficient use of resources. Uncertainties also affect the decisions of the community operator. We will continue to study these two aspects in future research.

Author Contributions

Conceptualization and methodology, Y.L. and X.W.; simulation and analysis, Y.L.; investigation, X.W.; data curation, Y.L.; writing—original draft preparation, Y.L.; writing—review and editing, Y.L. and X.W.; supervision, X.W.; literature research, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Natural Science Foundation of Xinjiang Uygur Autonomous Region under Grant 2020D01C031.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The equipment parameters.
Table A1. The equipment parameters.
EquipmentParameter TypeParameter Value
CHPRated Capacity300 kW
Minimum output power100 kW
Electrical efficiency fitting coefficient α CHP , 0 = 0.09 ,   α CHP , 1 = 0.44 ,   α CHP , 2 = 0.14 ,   α CHP , 3 = 0.11 ,   α CHP , 4 = 0.06
Thermoelectric ratio fitting coefficient α ψ , 0 = 3.82 ,   α ψ , 1 = 5.84 ,   α ψ , 2 = 3.6
Operation and maintenance cost0.04 RMB/(kWh)
EHPRated Capacity200 kW
Thermal efficiency fitting coefficient α EHP , 0 = 2.61 ,   α EHP , 1 = 0.36 ,   α EHP , 2 = 0.026
Operation and maintenance cost0.06 RMB/kWh
GBRated Capacity1000 kW
Thermal efficiency fitting coefficient α GB , 0 = 0.81 ,   α GB , 1 = 0.13
Operation and maintenance cost0.02 RMB/kWh
Electricity storageRated Capacity500 kW
Charge/Discharge efficiency0.98
Attrition rate0.02
Operation and maintenance cost0.01 RMB/kWh
Heat storageRated Capacity500 kW
Charge/Discharge efficiency0.95
Attrition rate0.02
Operation and maintenance cost0.01 RMB/kWh
Table A2. Data parameters for users.
Table A2. Data parameters for users.
ParameterMeaningValueParameterMeaningValue
π con , min e , g Dead threshold for electrical energy conversion0 π con , max e , g Saturation value of electrical energy conversion0.2
π con , min g , e Dead threshold for natural gas conversion0 π con , max g , e Saturation value of natural gas conversion0.15
f 1 e First power coefficient of electrical energy preference1.5 f 2 e Quadratic coefficient of electrical energy preference0.0009
f 1 h First power coefficient of thermal energy preference1.1 f 2 h Quadratic coefficient of thermal energy preference0.0011
f 1 g First power coefficient of natural gas preference1.2 f 2 g Quadratic coefficient of natural gas preference0.001
N 1 Number of the users with adjustable heating temperature500 N 2 Number of the users with non-adjustable heating temperature300
L tran , max e , out Maximum load that can be transferred out80 kW L tran , max e , in Maximum load that can be transferred in80 kW
t adj , max e Maximum duration of electrical load reduction4 h T set The most comfortable indoor temperature22.6 °C
Table A3. Time of use electricity prices.
Table A3. Time of use electricity prices.
PeriodMarket Electricity Price (RMB/kWh)
Valley period: 01:00—07:000.35
Normal period: 08:00—10:00; 15:00—17:00; 22:00—24:000.5
Peak period: 11:00—14:00; 18:00—21:000.8
Figure A1. Fitting curve of equipment efficiency.
Figure A1. Fitting curve of equipment efficiency.
Processes 10 01794 g0a1

References

  1. Bai, L.; Li, F.; Cui, H.; Jiang, T.; Sun, H.; Zhu, J. Interval optimization based operating strategy for gas-electricity integrated energy systems considering demand response and wind uncertainty. Appl. Energy 2016, 167, 270–279. [Google Scholar] [CrossRef]
  2. Wang, L.; Hou, C.; Ye, B.; Wang, X.; Yin, C.; Cong, H. Optimal Operation Analysis of Integrated Community Energy System Considering the Uncertainty of Demand Response. IEEE Trans. Power Syst. 2021, 36, 3681–3691. [Google Scholar] [CrossRef]
  3. Gokcek, T.; Sengor, I.; Erdinc, O. A Novel Multi-Hierarchical Bidding Strategies for Peer-to-Peer Energy Trading among Communities. IEEE Access 2022, 10, 23798–23807. [Google Scholar] [CrossRef]
  4. Mediwaththe, C.P.; Stephens, E.R.; Smith, D.B.; Mahanti, A. Competitive Energy Trading Framework for Demand-Side Management in Neighborhood Area Networks. IEEE Trans. Smart Grid. 2018, 9, 4313–4322. [Google Scholar] [CrossRef]
  5. Li, J.; Tian, L.; Cheng, L.; Guo, J. Optimization Planning of Micro-energy System Considering the Characteristics of Variable Operating Conditions (1) Basic Model and Analysis. Autom. Electr. Power Syst. 2018, 42, 18–26+49. [Google Scholar]
  6. Chen, W.; Mu, Y.; Jia, H.; Wei, W.; He, W. Optimization Scheduling Method for Regional Integrated Energy System Considering Equipment Variable Operating Condition Characteristics. Power Syst. Technol. 2021, 45, 951–958. [Google Scholar]
  7. Wang, Y.; Ma, Y.; Song, F.; Ma, Y.; Qi, C.; Huang, F.; Xing, J.; Zhang, F. Economic and efficient multi-objective operation optimization of integrated energy system considering electro-thermal demand response. Energy 2020, 205, 118022. [Google Scholar] [CrossRef]
  8. Guo, Z.; Zhang, R.; Wang, L.; Zeng, S.; Li, Y. Optimal Operation of regional integrated energy system considering demand response. Appl. Therm. Eng. 2021, 191, 116860. [Google Scholar] [CrossRef]
  9. Liu, R.; Yang, T.; Sun, G.; Lin, S.; Muyeen, S.M.; Ma, T. Optimal Dispatch of Community Integrated Energy System Considering Comprehensive User Satisfaction Method. IEEJ Trans. Electr. Electron. Eng. 2022. [Google Scholar] [CrossRef]
  10. Miao, B.; Lin, J.; Li, H.; Liu, C.; Li, B.; Zhu, X.; Yang, J. Day-Ahead Energy Trading Strategies of Regional Integrated Energy System Considering Energy Cascade Utilization. IEEE Access 2020, 8, 138021–138035. [Google Scholar] [CrossRef]
  11. Li, Y.; Feng, C.; Wen, F.; Wang, K.; Huang, Y. Energy Pricing and Management of Park Energy Internet Including Electric Vehicles and Electricity-to-Gas Conversion. Autom. Electr. Power Syst. 2018, 42, 192–196. [Google Scholar]
  12. Wei, W.; Liu, F.; Mei, S. Energy Pricing and Dispatch for Smart Grid Retailers Under Demand Response and Market Price Uncertainty. IEEE Trans. Smart Grid. 2015, 6, 1364–1374. [Google Scholar] [CrossRef]
  13. Paudel, A.; Chaudhari, K.; Long, C.; Gooi, H.B. Peer-to-Peer Energy Trading in a Prosumer-Based Community Microgrid: A Game-Theoretic Model. IEEE Trans. Ind. Electron. 2019, 66, 6087–6097. [Google Scholar] [CrossRef]
  14. Ma, L.; Liu, N.; Zhang, J.; Wang, L. Real-Time Rolling Horizon Energy Management for the Energy-Hub-Coordinated Prosumer Community from a Cooperative Perspective. IEEE Trans. Power Syst. 2019, 34, 1227–1242. [Google Scholar] [CrossRef]
  15. Fu, Y.; Li, H. Energy Management and Pricing Strategy of Multi-energy Supply Service in Parks Considering User Utility Evaluation. Modern Electric Power. 2022, 39, 317–328. [Google Scholar]
  16. Fleischhacker, A.; Corinaldesi, C.; Lettner, G.; Auer, H.; Botterud, A. Stabilizing Energy Communities Through Energy Pricing or PV Expansion. IEEE Trans. Smart Grid. 2022, 13, 728–737. [Google Scholar] [CrossRef]
  17. Anoh, K.; Maharjan, S.; Ikpehai, A.; Zhang, Y.; Adebisi, B. Energy Peer-to-Peer Trading in Virtual Microgrids in Smart Grids: A Game-Theoretic Approach. IEEE Trans. Smart Grid. 2020, 11, 1264–1275. [Google Scholar] [CrossRef]
  18. Wei, F.; Jing, Z.X.; Wu, P.Z.; Wu, Q.H. A Stackelberg game approach for multiple energies trading in integrated energy systems. Appl. Energy 2017, 200, 315–329. [Google Scholar] [CrossRef]
  19. Li, P.; Wang, Z.; Wang, N.; Yang, W.; Li, M.; Zhou, X.; Yin, Y.; Wang, J.; Guo, T. Stochastic robust optimal operation of community integrated energy system based on integrated demand response. Int. J. Electr. Power Energy Syst. 2021, 128, 106735. [Google Scholar] [CrossRef]
  20. Ruan, W.; Wang, B.; Li, Y.; Yang, S. Research on User Response Behavior Under Peak and Valley Time-of-use Electricity Prices. Power Syst. Technol. 2012, 36, 86–93. [Google Scholar]
  21. Geidl, M.; Koeppel, G.; Favre-Perrod, P.; Klockl, B.; Andersson, G.; Frohlich, K. Energy hubs for the future. IEEE Power Energy Mag. 2007, 5, 24–30. [Google Scholar] [CrossRef]
  22. Wang, Y.; Zhang, N.; Kang, C.; Kirschen, D.S.; Yang, J.; Xia, Q. Standardized Matrix Modeling of Multiple Energy Systems. IEEE Trans. Smart Grid. 2019, 10, 257–270. [Google Scholar] [CrossRef]
  23. Cui, Y.; Guo, F.; Fu, X.; Zhao, Y.; Han, C. Coordinated and optimized source-charge dispatch of an integrated energy system that promotes wind power consumption through energy-supply conversion. Power Syst. Technol. 2022, 46, 1437–1447. [Google Scholar]
  24. Li, Y.; Zhang, Z.; Yang, L.; Tang, D.; Dong, C.; Ji, C. A Decision-making Model for Incremental Power Distribution and Sales Companies Considering Collaborative Optimization of Source, Load and Storage. Autom. Electr. Power Syst. 2021, 45, 125–132. [Google Scholar]
  25. Giannelos, S.; Konstantelos, I.; Strbac, G. Option Value of Demand-Side Response Schemes Under Decision-Dependent Uncertainty. IEEE Trans. Power Syst. 2018, 33, 5103–5113. [Google Scholar] [CrossRef]
  26. Wu, C.; Gu, W.; Xu, Y.; Jiang, P.; Lu, S.; Zhao, B. Bi-level optimization model for integrated energy system considering the thermal comfort of heat customers. Appl. Energy 2018, 232, 607–616. [Google Scholar] [CrossRef]
  27. Sun, Y.; Zhang, C.; Li, Z.; Zhang, X.; Li, F. Flexible Load Multi-power-level Regulation Strategies Considering multi-regional User Differentiated PMV. Proc. CSEE 2021, 41, 7574–7586. [Google Scholar]
  28. Mu, Y.; Chen, W.; Yu, X.; Jia, H.; Hou, K.; Wang, C.; Meng, X. A double-layer planning method for integrated community energy systems with varying energy conversion efficiencies. Appl. Energy 2020, 279, 115700. [Google Scholar] [CrossRef]
  29. Zhao, H.; Miao, S.; Li, C.; Zhang, D.; Tu, Q. Research on the Optimal Operation Strategies of the Integrated Energy System in the Park Considering the Coupled Response Characteristics of Cooling, Heating and Electricity Demand. Proc. CSEE 2022, 42, 573–589. [Google Scholar]
  30. Maharjan, S.; Zhu, Q.; Zhang, Y.; Gjessing, S.; Basar, T. Dependable Demand Response Management in the Smart Grid: A Stackelberg Game Approach. IEEE Trans. Smart Grid. 2013, 4, 120–132. [Google Scholar] [CrossRef]
  31. Liu, R.; Li, Y.; Yang, X.; Li, Y.; Sun, G.; Shi, S. Two-stage Optimal Scheduling of Community Integrated Energy System Considering Demand Response. Acta Energ. Sol. Sin. 2021, 42, 46–54. [Google Scholar]
  32. Li, L. Development and Problems of Residential Centralized Cooling (Heating) in Hefei. Anhui Archit. 2021, 28, 79–80. [Google Scholar]
  33. Zhao, Y.; Gao, H.; Wang, Z.; Zhang, R.; Zhong, L.; Li, Y.; Liu, J. Optimization Decision of Multi-energy Trading Game for Business Park Operators Considering User Power Substitution. Power Syst. Technol. 2021, 45, 1320–1331. [Google Scholar]
Figure 1. Community Integrated Energy System Model.
Figure 1. Community Integrated Energy System Model.
Processes 10 01794 g001
Figure 2. Multi-energy transaction decision-making model architecture between the community operator and the users.
Figure 2. Multi-energy transaction decision-making model architecture between the community operator and the users.
Processes 10 01794 g002
Figure 3. Forecasting curves of renewable energy output, initial load, and outdoor temperature.
Figure 3. Forecasting curves of renewable energy output, initial load, and outdoor temperature.
Processes 10 01794 g003
Figure 4. The price of electricity and heat that the operator trades with users.
Figure 4. The price of electricity and heat that the operator trades with users.
Processes 10 01794 g004
Figure 5. Electricity and natural gas trading results between the community operator and the energy market in different scenarios.
Figure 5. Electricity and natural gas trading results between the community operator and the energy market in different scenarios.
Processes 10 01794 g005
Figure 6. Electricity and natural gas loads of users in different scenarios.
Figure 6. Electricity and natural gas loads of users in different scenarios.
Processes 10 01794 g006
Figure 7. The relation between the convertible load and the equivalent energy price difference.
Figure 7. The relation between the convertible load and the equivalent energy price difference.
Processes 10 01794 g007
Figure 8. Heat load and indoor temperature in Scenarios 1 and 3.
Figure 8. Heat load and indoor temperature in Scenarios 1 and 3.
Processes 10 01794 g008
Table 1. Comparative Study 1.
Table 1. Comparative Study 1.
Compare ItemsScenario 1Scenario 2Scenario 3
Community operatorCost (RMB)13,94012,48611,976
Profit (RMB)445943294720
Profit margin (%)31.934.739.4
Renewable energy utilization (%)889396
Community usersTotal cost (RMB)18,39916,81516,690
Value-added benefit15,12116,99317,412
Table 2. Comparative Study 2.
Table 2. Comparative Study 2.
Compare ItemsScenario 3Scenario 4
Community operatorCost (RMB)11,97613,153
Profit (RMB)47204957
Profit margin (%)39.437.6
Community usersTotal cost (RMB)16,690-
Value-added benefit17,412-
The optimization methodThe master–slave gameThe centralized optimization
Whether user interaction is considered
Whether the convertible load is considered×
Whether the retail energy prices have been optimized×
Whether a win–win situation has been achieved×
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Y.; Wang, X. Community Integrated Energy System Multi-Energy Transaction Decision Considering User Interaction. Processes 2022, 10, 1794. https://doi.org/10.3390/pr10091794

AMA Style

Li Y, Wang X. Community Integrated Energy System Multi-Energy Transaction Decision Considering User Interaction. Processes. 2022; 10(9):1794. https://doi.org/10.3390/pr10091794

Chicago/Turabian Style

Li, Yuantian, and Xiaojing Wang. 2022. "Community Integrated Energy System Multi-Energy Transaction Decision Considering User Interaction" Processes 10, no. 9: 1794. https://doi.org/10.3390/pr10091794

APA Style

Li, Y., & Wang, X. (2022). Community Integrated Energy System Multi-Energy Transaction Decision Considering User Interaction. Processes, 10(9), 1794. https://doi.org/10.3390/pr10091794

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

Article Metrics

Back to TopTop