Next Article in Journal
An Elastic Charging Service Fee-Based Load Guiding Strategy for Fast Charging Stations
Next Article in Special Issue
Implementation of a Real-Time Microgrid Simulation Platform Based on Centralized and Distributed Management
Previous Article in Journal
Mathematical Modelling of Mooring Systems for Wave Energy Converters—A Review
Previous Article in Special Issue
Aggregators’ Optimal Bidding Strategy in Sequential Day-Ahead and Intraday Electricity Spot Markets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Energy Trading and Pricing in Microgrids with Uncertain Energy Supply: A Three-Stage Hierarchical Game Approach

1
School of Electrical Engineering, Yanshan University, Qinhuangdao 066004, China
2
Department of Energy Technology, Aalborg University, Aalborg 9100, Denmark
*
Author to whom correspondence should be addressed.
Current address: Hebei Avenue 438, Qinhuangdao 066004, China
Energies 2017, 10(5), 670; https://doi.org/10.3390/en10050670
Submission received: 5 March 2017 / Revised: 22 April 2017 / Accepted: 6 May 2017 / Published: 11 May 2017
(This article belongs to the Special Issue Distributed Energy Resources Management)

Abstract

:
This paper studies an energy trading and pricing problem for microgrids with uncertain energy supply. The energy provider with the renewable energy (RE) generation (wind power) determines the energy purchase from the electricity markets and the pricing strategy for consumers to maximize its profit, and then the consumers determine their energy demands to maximize their payoffs. The hierarchical game is established between the energy provider and the consumers. The energy provider is the leader and the consumers are the followers in the hierarchical game. We consider two types of consumers according to their response to the price, i.e., the price-taking consumers and the price-anticipating consumers. We derive the equilibrium point of the hierarchical game through the backward induction method. Comparing the two types of consumers, we study the influence of the types of consumers on the equilibrium point. In particular, the uncertainty of the energy supply from the energy provider is considered. Simulation results show that the energy provider can obtain more profit using the proposed decision-making scheme.

1. Introduction

In the microgrid, energy trading is an important segment [1,2]. The energy provider determines the energy purchase to meet the consumer demands. Meanwhile, in order to increase its profit, the energy provider faces a problem of pricing decision. With the development of renewable energy (RE), it is reasonable for the energy provider to use the renewable energy as supply [3,4]. Due to the introduction of the renewable energy, the energy provider has to predict the generating capacity of the renewable energy system, and then decides how much energy it needs to purchase from the electricity markets. The energy provider’s prediction can have a deviation from the actual generation, which leads to the uncertainty (UC) of the energy supply [5,6,7].
There are some works in literature related to the interactions among the consumers. A non-cooperative game (NCG) was formulated among the consumers in [8,9,10,11,12,13,14]. In [9], the price-taking consumers and the price-anticipating consumers were considered. The price-taking consumers assume that their energy consumption cannot affect the electricity price, whereas the price-anticipating consumers believe that their energy consumption can change the electricity price. Recently, the Stackelberg game (SG) is formulated between the energy provider and the consumers in [15,16,17,18,19,20]. In addition, the authors in [20,21] took into account a two-stage Stackelberg game between the power station and the consumers. In order to reduce the cost of the energy purchased from the electricity markets, the renewable energy was taken into account in [22]. Most of these works mainly focus on the price-taking consumers, and they seldom take into account the renewable energy generation, thereby the uncertainty of the energy supply is not involved. The differences of the proposed work with the above literature are shown in Table 1.
In this paper, we consider the uncertainty of the energy supply caused by the wind power generation. Furthermore, we both consider the price-taking consumers and the price-anticipating consumers. We model the interactions between the energy provider and the consumers as a three-stage hierarchical game. The energy provider, which is the hierarchical game’s leader, determines the price and the energy purchase to maximize its profit. Finally, we obtain the equilibrium of the hierarchical game through the backward induction method [23,24].
The rest of the paper is organized as follows. Section 2 introduces the problem formulation. Section 3 shows the backward induction method of the three-stage hierarchical game for the energy provider and the price-taking consumers. In Section 4, the backward induction method of the three-stage game for the energy provider and the price-anticipating consumers is described. Section 5 gives the simulation and comparison results. Finally, the conclusions are summarized in Section 6.

2. Problem Formulation

We consider the energy trading and pricing problem in the microgrid consisting of one energy provider and a set N = { 1 , . . . , N } of consumers. The energy provider and the consumers are integrated into a microgrid with renewable energy generation. The energy provider purchases energy from the electricity markets when the renewable energy supply is not enough. In that case, the energy supply includes the energy generated from the renewable energy source and the energy purchased from the electricity markets. In the microgrid, the system structure of the energy trading is given in Figure 1. Because of the uncertainty of the energy generated from renewable energy sources, the energy purchased from the electricity markets is uncertain. According to the interaction between the energy provider and the consumers, we establish a hierarchical game as below.
  • Leader: the energy provider determines the energy purchase and the pricing strategy to maximize its profit.
  • Followers: the consumers determine the energy demands to maximize their payoffs.
According to the types of the consumers, we consider two scenarios in this paper. In scenario A, there is no competition among the price-taking consumers, i.e., the consumers’ energy consumption cannot affect the price announced by the energy provider. In scenario B, the interactions among the price-anticipating consumers are formulated into a non-cooperative game, i.e., the consumers’ energy consumption can change the price announced by the energy provider [9].

3. Wind Power Generation Model

There exists a large body of literature on wind power forecasting, and the day-ahead wind forecast based on numerical weather prediction (NWP) models can enable relatively accurate wind forecasts [25,26]. Because the operating time moves closer to the near term, the computation complexity often renders NWP models intractable at a high spatial resolution [26]. An adaptive wavelet neural network was proposed for mapping the NWP’s wind speed and wind direction forecasts to wind power forecasts in [27]. The authors in [28] proposed a novel statistical wind power forecast framework, which leverages the spatio-temporal correlation in wind speed and direction data among geographically dispersed wind farms. In [29], the author developed a feed-forward neural network approach for wind power generation forecasting to improve the wind forecasting accuracy. However, the wind power forecast is relatively complex, and the forecast errors cannot be avoided. Generally, the wind speed can be approximated as the Gamma distribution [30], inverse Gaussian [31], log-normal [32], and Weibull [33,34,35,36]. Alternatively, copula theory has recently been applied to wind speed and wind power as a way of modeling nonlinear dependence structures [37]. According to the wind speed, we establish the wind power generation model adopting the mixed copula function in this paper. Firstly, we introduce the copula theory.
The copula theory was proposed by Sklar in 1959 [38]. Supposing that F ( x 1 , x 2 , , x N ) is a joint distribution function whose marginal distributions are F 1 ( x 1 ) , F 2 ( x 2 ) , , F N ( x N ) , then there is a copula function C satisfies [39]:
F ( x 1 , x 2 , , x N ) = C ( F 1 ( x 1 ) , F 2 ( x 2 ) , , F N ( x N ) ) .
Defining that F 1 1 ( u 1 ) , F 2 1 ( u 2 ) , , F N 1 ( u N ) are the pseudo inverse functions of F 1 ( x 1 ) , F 2 ( x 2 ) , , F N ( x N ) , respectively. Therefore, the copula function C can be obtained as the following:
C ( u 1 , u 2 , , u N ) = F ( F 1 1 ( u 1 ) , F 2 1 ( u 2 ) , , F N 1 ( u N ) ) ,
where the marginal distributions of C ( u 1 , u 2 , , u N ) follow uniform distribution in [ 0 , 1 ] [40].
When N = 2 , the copula function C is a binary function. H ( x 1 , x 2 ) is a joint distribution function whose marginal distributions are F ( x 1 ) and W ( x 2 ) , and F 1 ( u ) and W 1 ( v ) are the pseudo inverse function of F ( x 1 ) and W ( x 2 ) , respectively. Therefore, the copula function C is expressed as the following:
C ( u , v ) = H ( F 1 ( u ) , W 1 ( v ) ) .
The mixed copula function was further proposed in [41,42]:
C M ( u , v ) = λ 1 C 1 ( u , v , γ 1 ) + λ 2 C 2 ( u , v , γ 2 ) + λ 3 C 3 ( u , v , γ 3 ) ,
where C M ( u , v ) is the mixed copula function that is composed of C 1 ( u , v , γ 1 ) , C 2 ( u , v , γ 2 ) , and C 3 ( u , v , γ 3 ) . λ 1 , λ 2 , and λ 3 are weight coefficients of C 1 ( u , v , γ 1 ) , C 2 ( u , v , γ 2 ) , and C 3 ( u , v , γ 3 ) , respectively, and satisfy λ 1 + λ 2 + λ 3 = 1 . γ 1 , γ 2 , and γ 3 are correlation coefficients and can measure the correlation degree of variables.
The results of [41] showed that the relevant structures of the mixed copula function are more flexible than a single copula function. The wind power generation model is established by the mixed copula function (see details in [38]).

4. Scenario A: The Three-Stage Game for Price-Taking Consumers

In this section, we prove the existence and uniqueness of the hierarchical equilibrium by using the backward induction method. First, we analyze the consumers’ demands given the energy provider’s pricing strategy and energy purchase. Then, we study the energy provider’s pricing strategy given the consumers’ energy demands and the energy purchase. Finally, we analyze the energy purchased from the electricity markets in the case of uncertain renewable generation, and then obtain the maximum profit of the energy provider.

4.1. Consumer’s Energy Demands in Stage III

In Stage I, the energy provider needs to determine the energy purchased from the electricity markets. In Stage II, the energy provider announces the price to the consumers. In Stage III, the consumers determine their energy demands given the unit price p announced by the energy provider in Stage II. The payoff of consumer i is defined as the difference between the satisfaction level and the payment for energy purchase, i.e.,
u i ( p , c i ) = h ( c i c i ) 2 p c i ,
where h is the consumers’ cost coefficient, c i is the actual energy demand of the consumer i, and c i is the energy demand of the consumer i to maintain normal operation of appliances. The consumers determine their demands to maximize their payoffs. p is a fixed price announced by the energy provider.
The first derivative of u i ( p , c i ) with respect to c i is:
u i ( p , c i ) c i = 2 h ( c i c i ) p .
Letting u i ( p , c i ) / c i = 0 , we obtain:
2 h ( c i c i ) p = 0 , i N .
From Equation (7), we can obtain the energy demands of the consumer i:
c i ( p ) = c i p 2 h .
Adding Equation (8) from 1 to N, we obtain the total energy consumption of all consumers:
i N c i * ( p ) = i N c i p N 2 h .
We assume that Q = i N c i for convenience. When all consumers’ total demands are i N c i * ( p ) , the payoffs of the consumers are at a maximum. Next, we consider how the energy provider makes the purchase strategy and pricing strategy in Stages I and II based on the total energy demands, respectively. In particular, we show that the energy provider can determine a price in Stage II such that the total energy demands (as a function of price) cannot exceed the total energy supply.

4.2. Optimal Pricing Strategy in Stage II

In Stage II, the energy provider determines the pricing strategy according to the consumers’ energy demands, given the energy purchase in Stage I. The profit of the energy provider is denoted by:
W ( p s , β ) = min ( p i N c i * ( p ) , p ( β + β 0 ) p s ) μ β 0 p s p w β p s ,
which is the difference between the revenue and the total cost. We assume that p s is the summation of the energy generated from the renewable energy source and the energy purchased from the electricity markets. β 0 p s indicates the energy purchase, β is the uncertainty factor of the renewable energy source, β p s is the wind power generation, and μ is the energy provider’s cost coefficient. In this paper, the energy provider’s cost comes from the energy purchase and the wind power generation. p w is the cost coefficient and p w β p s is the wind power generating cost. Equation (10) denotes the fact that the revenue of the energy provider is determined by the consumers’ demands subject to its available supply. In Stage II, the objective is to find the optimal price p that maximizes the energy provider’s profit, i.e.,
W I I ( p s , β ) = max p 0 W ( p s , β ) ,
where W I I ( p s , β ) denotes the maximum profit of the energy provider in Stage II. Since the energy supply p s from the energy provider is given in this stage, the total cost μ β 0 p s is already fixed. Therefore, the maximum revenue can be achieved by optimizing the price:
max p 0 min ( p i N c i * ( p ) , p ( β + β 0 ) p s ) .
Let us define the consumers’ energy demands D ( p ) = p i N c i * and the energy supply S ( p ) = p β p s . Then, we have:
D ( p ) = N p 2 2 h + Q p ,
S ( p ) = p ( β + β 0 ) p s .
From the above equations, we observe that D ( p ) is a quadratic function, and S ( p ) is a linear function. Thus, we can obtain the maximum point of D ( p ) at p d = Q h / N . The relationships between S ( p ) and D ( p ) are described in Figure 2.
  • S 1 ( p ) (excessive supply): S 1 ( p ) doesn’t intersect with D ( p ) , p * = p d ;
  • S 2 ( p ) (excessive supply): S 2 ( p ) has one intersection with D ( p ) , where D ( p ) has a non-negative slope, p * = p d ;
  • S 3 ( p ) (conservative supply): S 3 ( p ) has one intersection with D ( p ) , where D ( p ) has a negative slope, p * = p h , where p h is the intersection point of D ( p ) and S ( p ) and p * is the optimal price announced by the energy provider.
Letting D ( p ) = S ( p ) , we obtain a quadratic function with respect to p and make it equal to zero:
N p 2 2 h + ( Q ( β + β 0 ) p s ) p = 0 .
Solving the above Equation (15), we obtain the intersection point of D ( p ) and S ( p ) :
p h = 2 h ( Q ( β + β 0 ) p s ) N .
In the excessive supply regime, the maximum profit of the energy provider is at p = p d :
W I I E S = h Q 2 2 N μ β 0 p s p w β p s .
In the conservative supply regime, the maximum profit of the energy provider is at p = p h :
W I I C S = 2 h ( Q ( β + β 0 ) p s ) N ( β + β 0 ) p s μ β 0 p s p w β p s .
The optimal pricing decision and the corresponding optimal profit at Stage II are given in Table 2

4.3. Energy Supply Strategy in Stage I

In Stage I, the energy provider determines the energy purchase to maximize its profit by taking into account the uncertainty factor of the energy supply β [15]. The profit of the energy provider in the Stage I is given as follows:
W I = max p s 0 W I I ( p s , β ) ,
where W I I ( p s , β ) is the energy provider’s profit functions with respect to p s and the uncertain factor β obtained in Stage II.
We assume that the wind power generation P = β p s , and the minimum power and maximum power of the wind power generation are P min and P max , respectively. The probability density function of the wind power f W P ( P ) can be obtained in [38]. From Figure 2, we can obtain that the maximum consumers’ demands i N c i * ( p ) is Q / 2 when the price p is Q h / N . Thus, we consider the following two intervals:
(1)
Interval I: p s [ 0 , Q 2 ] . In this interval, the energy provider’s profit function is:
W I I 1 ( p s ) = E P [ P min , P max ] [ W I I C S ( P ) ] = P min P max W I I C S ( P ) f W P ( P ) d P .
(2)
Interval II: p s [ Q 2 , ] . The energy provider’s profit function is:
W I I 2 ( p s ) = E P [ P min , Q 2 ] [ W I I C S ( P ) ] + E P [ Q 2 , P max ] [ W I I E S ( P ) ] = P min Q 2 W I I C S ( P ) f W P ( P ) d P + Q 2 P max W I I E S ( P ) f W P ( P ) d P .
By comparing Interval I with Interval II, we can obtain the maximum profit of the energy provider and the optimal energy purchase in scenario A.

5. Scenario B: The Three-Stage Game for Price-Anticipating Consumers

Since the price is set by the energy provider based on the total energy consumption, the consumers are interactive with each other. Thus, we formulate a non-cooperative game among the consumers. The non-cooperative game has a unique Nash equilibrium if p ( c ) is a linear rotational symmetric function, and p ( c ) is formulated as follows [8]:
p ( c ) = ω i N c i + p 0 ,
where ω is a positive parameter to implement the elastic pricing, c i is the actual energy demands of the consumer i, and p 0 is a basic price.

5.1. Consumer’s Energy Demands in Stage III

In Stage I and Stage II, the energy provider determines the energy purchased from the electricity markets and the pricing strategy for the consumers, respectively. In Stage III, similar to Equation (5), we formulate the payoff of price-anticipating consumer i given the unit price p ( c ) announced by the energy provider as follows:
u i ( p ( c ) , c i ) = ( h ( c i c i ) 2 + b i ) p ( c ) c i ,
where h and c i were defined in Equation (5), and b i is a base value of the satisfaction level of consumer i and is different for each consumer, which reflects the flexibility of the consumers. The first derivative of u i ( c ) with respect to c i is:
u i ( c ) c i = 2 h ( c i c i ) ω c i ω i N c i p 0 .
Letting u i ( c ) / c i = 0 , we obtain:
2 h ( c i c i ) ω c i ω i N c i p 0 = 0 , i N .
Adding Equation (25) from 1 to N, we have:
2 h i N c i + 2 h i N c i ω i N c i ω N i N c i i N p 0 = 0 ,
from which we obtain the total energy consumption of all consumers:
i N c i * ( ω ) = 2 h i N c i N p 0 2 h + ω ( N + 1 ) .
To simplify the calculation process, we make:
2 h i N c i i N p 0 = G ,
and then we have:
i N c i * ( ω ) = G 2 h + ω ( N + 1 ) .

5.2. Optimal Pricing Strategy in Stage II

In Stage II, the energy provider determines the pricing strategy according to the consumers’ energy demands, given the energy purchase in Stage I. The profit of the energy provider is:
W ( p s , β ) = min ( p ( c ) i N c i * ( ω ) , p ( c ) ( β + β 0 ) p s ) μ β 0 p s p w β p s ,
and the maximum profit of the energy provider is:
W I I ( p s , β ) = max ω 0 W ( p s , β ) ,
where W I I ( p s , β ) denotes the maximum profit of the energy provider in Stage II. We can maximize the revenue of the energy provider by optimizing the price:
max ω 0 min ( p ( c ) i N c i * ( ω ) , p ( c ) ( β + β 0 ) p s ) .
Let us define the consumers’ total energy demands D ( ω ) = p ( c ) i N c i * and the energy supply S ( ω ) = p ( c ) ( β + β 0 ) p s . Then,
D ( ω ) = ω G 2 [ 2 h + ω ( N + 1 ) ] 2 + p 0 G 2 h + ω ( N + 1 ) ,
and the intersection point of D ( ω ) , and the y-axis is p 0 G / 2 h .
The first derivative of D ( ω ) with respect to ω is:
D ( ω ) ω = G 2 ( 2 h + ω ( N + 1 ) ) 2 2 ω ( N + 1 ) G 2 ( 2 h + ω ( N + 1 ) ) 3 p 0 ( N + 1 ) G ( 2 h + ω ( N + 1 ) ) 2 = 2 h G 2 ω ( N + 1 ) G 2 2 h p 0 ( N + 1 ) G p 0 ω ( N + 1 ) 2 G ( 2 h + ω ( N + 1 ) ) 3 .
When
ω < 2 h ( G p 0 ( N + 1 ) ) G ( N + 1 ) + p 0 ( N + 1 ) 2 ,
D ( ω ) ω > 0 , so D ( ω ) is an increasing function. When
ω > 2 h ( G p 0 ( N + 1 ) ) G ( N + 1 ) + p 0 ( N + 1 ) 2 ,
D ( ω ) ω < 0 , so D ( ω ) is an decreasing function.
Letting D ( ω ) / ω = 0 , we obtain:
ω 0 = 2 h ( G p 0 ( N + 1 ) ) G ( N + 1 ) + p 0 ( N + 1 ) 2
and
S ( ω ) = ω G ( β + β 0 ) p s 2 h + ω ( N + 1 ) + p 0 ( β + β 0 ) p s ,
and the intersection point of S ( ω ) , and the y-axis is p 0 ( β + β 0 ) p s .
The first derivative of S ( ω ) with respect to ω is:
S ( ω ) ω = G ( β + β 0 ) p s 2 h + ω ( N + 1 ) ω G ( β + β 0 ) p s ( N + 1 ) ( 2 h + ω ( N + 1 ) ) 2 = 2 h G ( β + β 0 ) p s ( 2 h + ω ( N + 1 ) ) 2 .
Since i N c i * > 0 and S ( ω ) / ω > 0 , S ( ω ) is an increasing function. The relationships between S ( ω ) and D ( ω ) are described in Figure 3.
(a) When G p 0 ( N + 1 ) < 0 , we can obtain the following conclusions from Figure 3a:
  • S 1 ( ω ) (excessive supply): p 0 ( β + β 0 ) p s p 0 G / 2 h , ω * = 0 ,
  • S 2 ( ω ) (conservative supply): p 0 ( β + β 0 ) p s < p 0 G / 2 h , ω * = ω p ,
where ω p is the intersection point of D ( ω ) and S ( ω ) , and ω * is the optimal parameter of the elastic price.
Because p ( c ) is a linear rotational symmetric function, the case with ω * = 0 is neglected.
(b) When G p 0 ( N + 1 ) 0 and p 0 ( β + β 0 ) p s p 0 G / 2 h , we have the conclusions by analyzing Figure 3b:
  • S 3 ( ω ) (excessive supply): S 3 ( ω ) has one intersection with D ( ω ) , where D ( ω ) has a non-negative slope, ω * = ω 0 ,
  • S 4 ( ω ) (conservative supply): S 4 ( ω ) has three intersections with D ( ω ) , ω * = ω p ,
  • S 5 ( ω ) (conservative supply): S 5 ( ω ) has one intersection with D ( ω ) , where D ( ω ) has a negative slope, ω * = ω p .
(c) When G p 0 ( N + 1 ) 0 and p 0 ( β + β 0 ) p s < p 0 G / 2 h , we can get the conclusions from Figure 3c:
  • S 6 ( ω ) (excessive supply): S 6 ( ω ) doesn’t intersect with D ( ω ) , ω * = ω 0 ,
  • S 7 ( ω ) (excessive supply): S 7 ( ω ) has one or two intersections with D ( ω ) , where both intersections are located in the increasing interval of D ( ω ) , ω * = ω 0 ,
  • S 8 ( ω ) (conservative supply): S 8 ( ω ) has two intersections with D ( ω ) , where both intersections are located in the both sides of ω 0 , respectively, ω * = ω p .
Letting D ( ω ) = S ( ω ) , we obtain a quadratic function with respect to ω and make it equal to zero:
( G ( N + 1 ) + p 0 ( N + 1 ) 2 ) ( β + β 0 ) p s ω 2 + [ ( 4 h p 0 ( N + 1 ) + 2 h G ) ( β + β 0 ) p s ( G p 0 ( N + 1 ) + G 2 ) ] ω + 4 h 2 p 0 ( β + β 0 ) p s 2 h p 0 G = 0 .
For convenience, we define:
A = G ( N + 1 ) + p 0 ( N + 1 ) 2 , B = 4 h p 0 ( N + 1 ) + 2 h G , C = G p 0 ( N + 1 ) + G 2 , D = 4 h 2 p 0 , E = 2 h p 0 G , Δ = ( B ( β + β 0 ) p s C ) 2 4 A ( β + β 0 ) p s ( D ( β + β 0 ) p s E ) .
Solving the above Equation (40), we obtain the intersection point of D ( ω ) and S ( ω ) :
ω p = ( B ( β + β 0 ) p s C ) + Δ 2 A ( β + β 0 ) p s .
In the excessive supply regime, the maximum profit of the energy provider is at ω = ω 0 :
W I I E S = [ G + p 0 ( N + 1 ) ] 2 8 h ( N + 1 ) μ β 0 p s p w β p s .
In the conservative supply regime, the maximum profit of the energy provider is at ω = ω p :
W I I C S = ( C B ( β + β 0 ) p s + Δ ) G ( β + β 0 ) p s 4 A h ( β + β 0 ) p s + ( C B ( β + β 0 ) p s + Δ ) ( N + 1 ) + p 0 ( β + β 0 ) p s μ β 0 p s p w β p s .
The optimal pricing decision and the corresponding optimal profit in Stage II are given in Table 3.

5.3. Energy Supply Strategy in Stage I

In Stage I, the energy provider also determines the energy purchase to maximize its profit by taking into account the uncertainty of the energy supply. The profit of the energy provider in the Stage I is given by:
W I = max p s 0 W I I ( p s , β ) ,
where W I I ( p s , β ) is the energy provider’s profit function with respect to p s and the uncertain factor β obtained in Stage II.
We assume that the wind power generation P = β p s , and the minimum power and maximum power of the wind power generation are P min and P max , respectively. The probability density function of the wind power f W P ( P ) can be obtained in [38]. From Figure 3, we can obtain that the maximum consumers’ demands i N c i * ( ω ) is A / 4 h ( N + 1 ) when ω = ω 0 . For convenience, we assume that L = A / 4 h ( N + 1 ) and consider the following two intervals:
(1)
Interval I: p s [ 0 , A 4 h ( N + 1 ) ] . In this interval, the energy provider’s profit function is:
W I I 1 ( p s ) = E P [ P min , P max ] [ W I I C S ( P ) ] = P min P max W I I C S ( P ) f W P ( P ) d P .
(2)
Interval II: p s [ A 4 h ( N + 1 ) , ] . The energy provider’s profit function is:
W I I 2 ( p s ) = E P [ P min , A 4 h ( N + 1 ) ] [ W I I C S ( P ) ] + E P [ A 4 h ( N + 1 ) , P max ] [ W I I E S ( P ) ] = P min A 4 h ( N + 1 ) W I I C S ( P ) f W P ( P ) d P + A 4 h ( N + 1 ) P max W I I E S ( P ) f W P ( P ) d P .
Similar to scenario A, we can obtain the maximum profit of the energy provider and the optimal amount of energy purchased from the electricity markets.

6. Simulation Results

This section presents simulation studies of the proposed scheme using MATLAB 7.11.0 (MathWorks, Natick, MA, USA). In the simulations, we assume that the wind power generation follows a uniform distribution in [ P min , P max ] , and we select that h = 0 . 04 , p 0 = 0 . 1 , N = 100 , and μ = 10 . For the parameter c i , we select a set of stochastic values. Then, we can obtain the profit of the energy provider under different β 0 for two scenarios as shown in Figure 4 and Figure 5, respectively.
From Figure 4 and Figure 5, we observe that the optimal energy supply decreases with the increase of the β 0 and the maximum profit is changed from Interval II to Interval I. In general, the wind power generation cost is less than the cost of purchasing energy. From the profit function of the energy provider, when β 0 increases, only by decreasing p s can the profit of the energy provider be maximized. Thus, it is verified that the proposed method is effective, and the simulation values of Figure 4 and Figure 5 are shown in Table 4.
Taking β 0 = 0 . 1 as an example, the comparisons between the two scenarios are shown in Figure 6. From Figure 6 and Table 4, we observe that the energy provider can obtain more profit in scenario A.
To explain the effect of the uncertainty, taking β 0 = 0 . 6 under scenario A as an example, we show the profit of the energy provider under the certain and uncertain energy supply in Figure 7. It is observed that the energy provider can achieve the higher profit under the certain energy supply. In reality, the uncertainty of the energy supply is necessary because the energy generated from the renewable energy sources is uncertain.

7. Conclusions

In this paper, we establish a model for energy trading and pricing in the microgrid. We formulate a hierarchical game between the energy provider with the renewable energy generation and the consumers, e.g., the price-taking consumers and the price-anticipating consumers. In the hierarchical game, the energy provider acts as the leader and the consumers act as the followers. The equilibrium point of the hierarchical game is obtained through the backward induction method. Furthermore, we also consider the uncertainty of the energy supply in the problem formulation. The simulation results show that the optimal energy supply can be obtained based on the reasonable pricing strategy and purchase strategy. Comparing the price-taking consumers with the price-anticipating consumers, we can obtain that the energy provider obtains more profit from the price-taking consumers. From the simulation results, we also can obtain that the energy provider’s profit reduces because of the uncertainty of the energy supply.
However, we do not consider that the consumers can sell the energy to the energy provider when the consumers have photovoltaic (PV) panels and a storage system. In that case, the energy demands of the consumers will be uncertain, and the payoff of the consumer includes two additional parts: one part is the PV generation cost, and the other part is the uncertainty of the energy demands. In order to compute the payoff of the consumer, we need to know the distribution that the PV generation follows. Then, we can get the average payoff of the consumer by expectation, and the optimal energy demands are obtained by the derivation method. Simultaneously, the profit of the energy provider needs to introduce two additional parts that denote buying the energy from the consumers and selling the energy to the electricity markets, which will be considered in the future.

Acknowledgments

This research was supported in part by the National Natural Science Foundation of China under Grants 61573303, 61503324, and 61573300, in part by the Natural Science Foundation of Hebei Province under Grants F2016203438, E2017203284, and E2016203092, in part by a project funded by the China Postdoctoral Science Foundation under Grants 2015M570233 and 2016M601282, in part by a project funded by the Hebei Education Department under Grant BJ2016052, and in part by the Technology Foundation for Selected Overseas Chinese Scholar under Grant C2015003052

Author Contributions

Kai Ma and Shubing Hu wrote the paper; Jie Yang analyzed the data; and Chunxia Dou and Josep M. Guerrero contributed the analysis tools.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, H.; Li, Y.; Han, Z.; Vucetic, B. A stackelberg game-based energy trading scheme for power beacon-assisted wireless-powered communication. In Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Brisbane, Australia, 19–24 April 2015; pp. 3177–3181. [Google Scholar]
  2. Misra, S.; Bera, S.; Ojha, T.; Mouftah, H.T.; Anpalagan, A. ENTRUST: Energy trading under uncertainty in smart grid systems. Comput. Netw. 2016, 110, 232–242. [Google Scholar] [CrossRef]
  3. Belgana, A.; Rimal, B.P.; Maier, M. Open energy market strategies in microgrids: A Stackelberg game approach based on a hybrid multiobjective evolutionary algorithm. IEEE Trans. Smart Grid 2015, 6, 1243–1252. [Google Scholar] [CrossRef]
  4. Jia, L.; Tong, L. Dynamic pricing and distributed energy management for demand response. IEEE Trans. Smart Grid 2016, 7, 1128–1136. [Google Scholar] [CrossRef]
  5. Duan, L.; Huang, J.; Shou, B. Investment and pricing with spectrum uncertainty: A cognitive operator’s perspective. IEEE Trans. Mob. Comput. 2011, 10, 1590–1604. [Google Scholar] [CrossRef]
  6. Hu, M.C.; Lu, S.Y.; Chen, Y.H. Stochastic–multiobjective market equilibrium analysis of a demand response program in energy market under uncertainty. Appl. Energy 2016, 182, 500–506. [Google Scholar] [CrossRef]
  7. Nie, S.; Huang, C.Z.; Huang, G.H.; Li, Y.P.; Chen, J.P.; Fan, Y.R.; Cheng, G.H. Planning renewable energy in electric power system for sustainable development under uncertainty—A case study of Beijing. Appl. Energy 2016, 162, 772–786. [Google Scholar] [CrossRef]
  8. Ma, K.; Hu, G.; Spanos, C.J. Distributed energy consumption control via real-time pricing feedback in smart grid. IEEE Trans. Control Syst. Technol. 2014, 22, 1907–1914. [Google Scholar]
  9. Ma, K.; Hu, G.; Spanos, C.J. A cooperative demand response scheme using punishment mechanism and application to industrial refrigerated warehouses. IEEE Trans. Ind. Inform. 2015, 11, 1520–1531. [Google Scholar] [CrossRef]
  10. Mohsenian-Rad, A.H.; Wong, V.W.S.; Jatskevich, J.; Schober, R.; Leon-Garcia, A. Autonomous demand-side management based on game-theoretic energy consumption scheduling for the future smart grid. IEEE Trans. Smart Grid 2010, 1, 320–331. [Google Scholar] [CrossRef]
  11. Baharlouei, Z.; Hashemi, M.; Narimani, H.; Mohsenian-Rad, H. Achieving optimality and fairness in autonomous demand response: Benchmarks and billing mechanisms. IEEE Trans. Smart Grid 2013, 4, 968–975. [Google Scholar] [CrossRef]
  12. Yu, M.; Hong, S.H. Supply–demand balancing for power management in smart grid: A Stackelberg game approach. Appl. Energy 2016, 164, 702–710. [Google Scholar] [CrossRef]
  13. Gao, B.; Ma, T.; Tang, Y. Power transmission scheduling for generators in a deregulated environment based on a game-theoretic approach. Energies 2015, 8, 13879–13893. [Google Scholar] [CrossRef]
  14. Liu, N.; Wang, C.; Lin, X.; Lei, J. Multi-party energy management for clusters of roof leased PV prosumers: A game theoretical approach. Energies 2016, 9, 536. [Google Scholar] [CrossRef]
  15. Bu, S.; Yu, F.R. A game-theoretical scheme in the smart grid with demand-side management: Towards a smart cyber-physical power infrastructure. IEEE Trans. Emerg. Top. Comput. 2013, 1, 22–32. [Google Scholar] [CrossRef]
  16. Maharjan, S.; Zhu, Q.; Zhang, Y.; Gjessing, S. Dependable demand response management in the smart grid: A stackelberg game approach. IEEE Trans. Smart Grid 2013, 4, 120–132. [Google Scholar] [CrossRef]
  17. Soliman, H.M.; Leon-Garcia, A. Game-theoretic demand-side management with storage devices for the future smart grid. IEEE Trans. Smart Grid 2014, 5, 1475–1485. [Google Scholar] [CrossRef]
  18. Maharjan, S.; Zhu, Q.; Zhang, Y.; Gjessing, S.; Basar, T. Demand response management in the smart grid in a large population regime. IEEE Trans. Smart Grid 2016, 7, 189–199. [Google Scholar] [CrossRef]
  19. Lee, J.; Guo, J.; Choi, J.K.; Zukerman, M. Distributed energy trading in microgrids: a game-theoretic model and its equilibrium analysis. IEEE Trans. Ind. Electron. 2015, 62, 3524–3533. [Google Scholar] [CrossRef]
  20. Yoon, S.G.; Choi, Y.J.; Park, J.K.; Bahk, S. Demand response design based on a Stackelberg game in smart grid. In Proceedings of the International Conference on ICT Convergence, Jeju, Korea, 14–16 October 2013; pp. 177–178. [Google Scholar]
  21. Fadlullah, Z.M.; Quan, D.M.; Kato, N.; Stojmenovic, I. GTES: An optimized game-theoretic demand-side management scheme for smart grid. IEEE Syst. J. 2014, 8, 588–597. [Google Scholar] [CrossRef]
  22. Yang, B.; Li, J.; Han, Q.; He, T.; Chen, C.; Guan, X. Distributed control for charging multiple electric vehicles with overload limitation. IEEE Trans. Parallel Distrib. Syst. 2016, 27, 3441–3454. [Google Scholar] [CrossRef]
  23. Abegaz, B.W.; Mahajan, S.M. Optimal dispatch control of energy storage systems using forward-backward induction. In Proceedings of the 2015 International Conference on Clean Electrical Power (ICCEP), Taormina, Italy, 16–18 June 2015; pp. 731–736. [Google Scholar]
  24. Cho, J.; Kleit, A.N. Energy storage systems in energy and ancillary markets: A backwards induction approach. Appl. Energy 2015, 147, 176–183. [Google Scholar] [CrossRef]
  25. Mahoney, W.P.; Parks, K.; Wiener, G.; Liu, Y.; Myers, W.L.; Sun, J.; Monache, L.D.; Hopson, T.; Johnson, D.; Haupt, S.E. A wind power forecasting system to optimize grid integration. IEEE Trans. Sustain. Energy 2012, 3, 670–682. [Google Scholar] [CrossRef]
  26. Constantinescu, E.M.; Zavala, V.M.; Rocklin, M.; Lee, S.; Anitescu, M. A computational framework for uncertainty quantification and stochastic optimization in unit commitment with wind power generation. IEEE Trans. Power Syst. 2011, 26, 431–441. [Google Scholar] [CrossRef]
  27. Kanna, B.; Singh, S.N. Long term wind power forecast using adaptive wavelet neural network. In Proceedings of the 2016 IEEE Uttar Pradesh Section International Conference on Electrical, Computer and Electronics Engineering (UPCON), Varanasi, India, 9–11 December 2016; pp. 671–676. [Google Scholar]
  28. Xie, L.; Gu, Y.; Zhu, X.; Genton, M.G. Short-term spatio-temporal wind power forecast in robust look-ahead power system dispatch. IEEE Trans. Smart Grid 2014, 5, 511–520. [Google Scholar] [CrossRef]
  29. Finamore, A.R.; Galdi, V.; Calderaro, V.; Piccolo, A.; Conio, G.; Grasso, S. Artificial neural network application in wind forecasting: An one-hour-ahead wind speed prediction. In Proceedings of the 5th IET International Conference on Renewable Power Generation (RPG), London, UK, 21–23 September 2016; pp. 1–6. [Google Scholar]
  30. Sherlock, R.H. Analyzing winds for frequency and duration on atmospheric pollution. Am. Meteorol. Soc. 1951, 4, 42–49. [Google Scholar]
  31. Bardsley, W.E. Note on the use of the inverse Gaussian distribution for wind energy applications. J. Appl. Meteorol. 1980, 19, 1126–1130. [Google Scholar] [CrossRef]
  32. Luna, R.E.; Church, H.W. Estimation of long-term concentrations using a `universal’ wind speed distribution. J. Appl. Meteorol. 1974, 13, 910–916. [Google Scholar] [CrossRef]
  33. Hennessey, J.P.J. A comparison of the Weibull and Rayleigh distributions for estimating wind power potential. Wind Eng. 1978, 2, 156–164. [Google Scholar]
  34. Justus, C.G.; Hargraves, W.R.; Yalcin, A. Nationwide assessment of potential output from wind-powered generators. J. Appl. Meteorol. 1976, 15, 673–678. [Google Scholar] [CrossRef]
  35. Stewart, D.A.; Essenwanger, O.M. Frequency distribution of wind speed near the surface. J. Appl. Meteorol. 1978, 17, 1633–1642. [Google Scholar] [CrossRef]
  36. Takle, E.S.; Brown, J.M. Note on the use of Weibull statistics to characterize wind-speed data. J. Appl. Meteorol. 1978, 17, 556–559. [Google Scholar] [CrossRef]
  37. Liu, S.; Li, G.; Xie, H.; Wang, X. Correlation characteristic analysis for wind speed in different geographical hierarchies. Energies 2017, 10, 237. [Google Scholar] [CrossRef]
  38. Pan, X.; Wang, L.; Xu, Y.; Zhang, L.; Liu, W.; Wu, R. A wind farm power modeling method based on mixed Copula. Dianli Xitong Zidonghua/Autom. Electr. Power Syst. 2014, 38, 17–22. [Google Scholar]
  39. Du, M.; Yi, J.; Mazidi, P.; Cheng, L.; Guo, J. A parameter selection method for wind turbine health management through SCADA data. Energies 2017, 10, 253. [Google Scholar] [CrossRef]
  40. Ma, C. Robust exponential stability of reaction-diffusion generalized Cohen-Grossberg neural networks with distributed delays. J. Xinjiang Norm. Univ. 2007, 26, 18–24. [Google Scholar]
  41. Li, X.M.; Shi, D.J. Research on dependence structure between shanghai and shenzhen stock markets. Appl. Stat. Manag. 2006, 25, 729–736. [Google Scholar]
  42. Hu, L. Dependence patterns across financial markets: A mixed copula approach. Appl. Financ. Econ. 2006, 16, 717–729. [Google Scholar] [CrossRef]
Figure 1. The system structure of energy trading.
Figure 1. The system structure of energy trading.
Energies 10 00670 g001
Figure 2. The relationships between S ( p ) and D ( p ) .
Figure 2. The relationships between S ( p ) and D ( p ) .
Energies 10 00670 g002
Figure 3. The relationships between S ( ω ) and D ( ω ) under the different conditions.
Figure 3. The relationships between S ( ω ) and D ( ω ) under the different conditions.
Energies 10 00670 g003
Figure 4. Scenario A: the profit of the energy provider under different β 0 .
Figure 4. Scenario A: the profit of the energy provider under different β 0 .
Energies 10 00670 g004
Figure 5. Scenario B: the profit of the energy provider under different β 0 .
Figure 5. Scenario B: the profit of the energy provider under different β 0 .
Energies 10 00670 g005
Figure 6. Comparisons between Scenario A with Scenario B.
Figure 6. Comparisons between Scenario A with Scenario B.
Energies 10 00670 g006
Figure 7. The effect of the uncertainty of the energy supply.
Figure 7. The effect of the uncertainty of the energy supply.
Energies 10 00670 g007
Table 1. Differences of the proposed work with the literature.
Table 1. Differences of the proposed work with the literature.
IndexesREUCNCGSG
[3,4]×××
[5,6,7]×××
[8,9,10,11,12,13,14]×××
[15,16,17,18,19,20,21]×××
This work
Table 2. Optimal pricing decision and profit in Stage II in scenario A.
Table 2. Optimal pricing decision and profit in Stage II in scenario A.
Total Energy ObtainedOptimal PriceOptimal Profit
in Stages I and II p * ( p s , β ) W II ( p s , β )
Excessive Supply Regime: p s Q 2 p E S = p d W I I E S ( p s , β ) in Equation (17)
Conservative Supply Regime: p s < Q 2 p C S = p h W I I C S ( p s , β ) in Equation (18)
Table 3. Optimal pricing decision and profit in Stage II in scenario B.
Table 3. Optimal pricing decision and profit in Stage II in scenario B.
Total Energy ObtainedOptimal ParameterOptimal Profit
in Stages I and II p * ( p s , β ) W II ( p s , β )
Excessive Supply Regime: p s A 4 h ( N + 1 ) ω E S = ω 0 W I I E S ( p s , β ) in Equation (43)
Conservative Supply Regime: p s < A 4 h ( N + 1 ) ω C S = ω p W I I C S ( p s , β ) in Equation (44)
Table 4. Simulation values of the two scenarios.
Table 4. Simulation values of the two scenarios.
Scenario AScenario B
β 0 p s Profit p s Profit
0.134933.79399.631.52
0.324634.61288.633.27
0.618635.220935.13

Share and Cite

MDPI and ACS Style

Ma, K.; Hu, S.; Yang, J.; Dou, C.; Guerrero, J.M. Energy Trading and Pricing in Microgrids with Uncertain Energy Supply: A Three-Stage Hierarchical Game Approach. Energies 2017, 10, 670. https://doi.org/10.3390/en10050670

AMA Style

Ma K, Hu S, Yang J, Dou C, Guerrero JM. Energy Trading and Pricing in Microgrids with Uncertain Energy Supply: A Three-Stage Hierarchical Game Approach. Energies. 2017; 10(5):670. https://doi.org/10.3390/en10050670

Chicago/Turabian Style

Ma, Kai, Shubing Hu, Jie Yang, Chunxia Dou, and Josep M. Guerrero. 2017. "Energy Trading and Pricing in Microgrids with Uncertain Energy Supply: A Three-Stage Hierarchical Game Approach" Energies 10, no. 5: 670. https://doi.org/10.3390/en10050670

APA Style

Ma, K., Hu, S., Yang, J., Dou, C., & Guerrero, J. M. (2017). Energy Trading and Pricing in Microgrids with Uncertain Energy Supply: A Three-Stage Hierarchical Game Approach. Energies, 10(5), 670. https://doi.org/10.3390/en10050670

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