Next Article in Journal
Magnetohydrodynamics Flow Past a Moving Vertical Thin Needle in a Nanofluid with Stability Analysis
Next Article in Special Issue
An Integration Mechanism between Demand and Supply Side Management of Electricity Markets
Previous Article in Journal
Derating Guidelines for Lithium-Ion Batteries
Previous Article in Special Issue
The Impact of Substituting Production Technologies on the Economic Demand Response Potential in Industrial Processes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Role of Demand Response Aggregators and the Effect of GenCos Strategic Bidding on the Flexibility of Demand

1
Department of Electrical Engineering and Electronic Engineering, Chittagong University of Engineering and Technology, Chittagong 4349, Bangladesh
2
School of Electrical Engineering and Computer Science, Queensland University of Technology, Brisbane 4000, Australia
*
Author to whom correspondence should be addressed.
Energies 2018, 11(12), 3296; https://doi.org/10.3390/en11123296
Submission received: 30 October 2018 / Revised: 15 November 2018 / Accepted: 20 November 2018 / Published: 26 November 2018
(This article belongs to the Special Issue Demand Response in Electricity Markets)

Abstract

:
This paper presents an interactive trading decision between an electricity market operator, generation companies (GenCos), and the aggregators having demand response (DR) capable loads. Decisions are made hierarchically. At the upper-level, an electricity market operator (EMO) aims to minimise generation supply cost considering a DR transaction cost, which is essentially the cost of load curtailment. A DR exchange operator aims to minimise this transaction cost upon receiving the DR offer from the multiple aggregators at the lower level. The solution at this level determines the optimal DR amount and the load curtailment price. The DR considers the end-user’s willingness to reduce demand. Lagrangian duality theory is used to solve the bi-level optimisation. The usefulness of the proposed market model is demonstrated on interconnection of the Pennsylvania-New Jersey-Maryland (PJM) 5-Bus benchmark power system model under several plausible cases. It is found that the peak electricity price and grid-wise operation expenses under this DR trading scheme are reduced.

1. Introduction

Demand-side management (DSM) with energy security and affordability objective has opened up new challenges and opportunities for the electricity market participants [1]. The term, DSM, was first introduced by the Electric Power Research Institute (EPRI) in the 1980s [2]. It refers to the change in the usual electricity consumption to optimized consumption by coordinating the demand response (DR) capable loads, renewable resources, on-site generation, and energy storage [3].
Specifically, a variety of demand response (DR) programs as part of the DSM initiatives have been used as potential resources to balance supply and demand, reduce operation cost, and to enhance power system efficiency [4,5]. The DR is defined as a change in electricity usages by the end-user customer from a usual consumption, in response to a changing electricity price [6,7]. It can be promoted via monetary incentives to lower energy consumption when the market prices are high [8]. The market price becomes higher because of power system network congestion and a power supply deficit [9].
Development of DR is triggered by policies, market mechanisms, and implementation frameworks. The Federal Electricity Regulatory Commission (FERC) in the USA assessed that DSM could reduce 9.2% of the peak demand, which is accounted to be 72,000 MW [10]. The ever-increasing DR potential is attributed to new and expanded programs. It urges market reformation toward end-user compensation and economic reward from DR service providers.
The DSM undertaken by the active participation of end-user customers is a key aspect of a smart power system [11,12,13,14,15,16,17,18,19,20]. In optimized consumption, peak load demand clips, and the valleys load demand fill [5]. From the grid operator perspective, the DSM measures substantially minimize the number of expensive peaking generating plants’ activation [21]. Furthermore, electricity networks become more sustainable, efficient, and secure.
The third-party aggregator as a profit-making liaising agent could be a key driver to make the DSM smarter via technology adoption and enhanced end-users’ interaction [22]. Aggregators will have to possess the necessary enabling technologies to manage varieties of the DR programs for the end-user customers while enticing them with appropriate compensation. The aggregators, in turn, receive rewards from the market operator for the cumulative DR services and help in shaping the load demand [23].
This paper proposes a DR integrated bi-level optimization framework. The bi-level optimization is a special kind of mathematical program, where an optimization problem contains another optimization problem as a constraint [24,25,26]. The outer and inner optimization tasks in the bi-level set are commonly referred to as the upper-level and lower-level optimization problems, respectively [27]. The bi-level model in this paper consists of an electricity market operator’s operation cost minimization problem at the upper-level, and a DR transaction cost (DRTC) minimization problem at the lower-level. These problems use a hierarchical structure for the decision making of two independent decision makers. The transaction cost accounts for the costs associated with the disutility of the DR capable loads. Disutility, in this paper, refers to the equivalent monetary value of the DR that could not be achieved [22,28].
The DR integration by dynamic pricing design is investigated in [29], and an end-user’s disutility is introduced in the aggregator’s DR supply offer curve in [30]. Although a demand response exchange (DRX) is introduced, its interaction with an energy market operator (EMO) is ignored [29,31]. Authors in [32] proposed varieties of DR programs to improve day-ahead market efficiency by using a multi-attribute decision-making approach. The DR purchase and supply offers bids are managed by a DRX operator (DRXO) to create a competitive DR market environment. The DRXO clears the DR bidding and aims at minimizing the DR transaction cost. The importance of aggregator and market modelling involving multiple aggregators as an independent mediator between the users and EMO is a challenging task as recognized by researchers in [33,34,35]. Automated DR for an industry in a smart grid paradigm is an exciting research area. Considering the practical scheduling constraints of the industry in the DR program can be challenging, the widespread implementation of these programs is restricted as reported [36]. Overall, the effect of DR bidding on the wholesale electricity prices and operating expense at grid-level has been largely neglected in previous works.
The objectives of this research are to: (1) Propose a DRX trading framework and integrate it into a security constraint electricity market clearing model; (2) manage supply deficits by using demand-side resources in a more engaging way, and from an economic point of view, using monetary value as a key operational parameter; (3) quantify how the aggregator and end-users may benefit financially when the GenCos exercise their conflicting economic interest to uplift market price; (4) evaluate the DR compensation price and quantity settlement for the participants in DR trading and the end-user customers; and (5) investigate a variety of bidding behaviors, such as strategic, regulated, and competitive, and its impact on the aggregator’s profit.
The main contribution of this paper is the introducing of a DR bidding mechanism at the lower-level with competing aggregators and its impact on the bidding strategy of generation companies (GenCos). We provide a concrete basis of why the approach of utilizing a DRX and having a bi-level framework rather than integrating with the existing ancillary service structure of the wholesale market is useful. Additionally, if a DRX mechanism is integrated, then how do we set its involvement margin? Overall, a comprehensive oversight of supply and demand-side management from a network perspective and from an economic point of view is investigated in this research.
In the DRX market, the aggregators compete by submitting their bids in the form of supply functions [28]. Aggregators participate in the DRX market while purchasing DR resources from the end-user customers. At the upper-level, this paper leverage the locational marginal prices (LMPs) based method where the LMP is found at each power system buses [37]. The changes in the LMP depends on the system conditions, such as transmission constraints, generation capacity, supply offer price, loading level, and so forth [38]. The aggregator gets paid at the LMP rate and the payment received by the aggregator is indeed a portion of operation cost saving due to the reduced market clearing price.
The end-user’s inconvenience due to participation in the DR is captured by a disutility cost function, which increases with the offered DR [39], and receives a disutility compensation, which is the part of the LMP payment received by the aggregator from the EMO. The end-users reduce energy consumption cost, anticipating the compensation from aggregators, while the EMO reduces the system overall operation cost. This compensation payment is a cost to the aggregators and the DR transaction is advantageous as long as the compensation given to the end-users does not outweigh the operation cost. Therefore, the EMO needs to trade-off between the DR transaction and operation cost in a competitive manner.
In the proposed market, aggregators offer DR supply on behalf of the end-user customers. Upon procuring the required DR, the EMO modifies the load demand at different power system buses and the GenCos bid supply offer to cater the modified load demand. Then, the EMO recalculates the new supply share among GenCos, and the market is settled for minimal operation cost by solving economic dispatch. The economic dispatch problem is subject to the power balance at every bus, thermal limits of transmission lines, power limits for the GenCos, and the aggregated DR limit at DR capable buses.
The rest of this paper is organized as follows: The problem formulation for the proposed market model is described in Section 2. The problem is formulated as bi-level programming with the EMO’s objective at the upper-level and DRX’s objective at the lower-level. Sections 3 presents a computational setup considering benchmark network and data. Several cases are formed to investigate the model performance. The simulation results in Section 4 provide the analysis of decisions of the market player in both the levels. Finally, Section 5 provides the concluding remarks of this research.

2. Problem Formulation

The interactions among several market participants at different levels can be summarised using Figure 1. The upper-level deals a security constraint economic dispatch (SCED) while the lower-level is about the DRX and the aggregators. The division of the market mechanism in DRX can be either horizontal among aggregators or vertical among the end-users and the aggregator. In the proposed DRX mechanism, the users get paid by the aggregator for their DR capable loads. To solve the DRX integrated market clearing model, a bi-level optimization technique has been used [26,40].
The DR supply offer is made available in a separate DRX market at the lower-level. The user’s inconvenience cost is minimized while achieving the targeted DR at the upper-level. The interaction between the aggregators and the users in the DRX market is modelled as a mechanism design problem. The aggregators make a demand reduction bid, and the DRX operator (DRXO) determines the DR share among the aggregators to minimize the DR transaction cost. At the upper-level, the interaction between GenCos and EMO is modelled as a SCED problem, where the generation supply offer bid is settled by the EMO, aiming to minimize operation cost based on the reduced demand and transaction cost at the lower-level.
The bi-level interaction is formulated as a mathematical programming problem with complementarity constraints [41,42,43]. The aggregator’s payoff comprises of the difference of DR selling revenue at the locational marginal price (LMP) and compensation provided to the customers. The LMP and compensation price is obtained by solving the lower and upper-level problem, respectively. A Lagrangian duality with associated Karush Kuhn Tucker (KKT) optimality conditions is used to solve the bi-level optimization. A gradient descent algorithm is also identified to solve bi-level electricity market model as reported in [44]. The effectiveness of the proposed market model is demonstrated on a PJM 5-Bus benchmark power system considering several case studies. The proposed DR integrated market model minimizes the operation cost, and maximizes the aggregator’s payoff within a certain incremental DR level.
Let us consider a set, T := {0, …, Tk − 1}, is a time period, and assume that the number of time interval Tk satisfies Tk = 24/k (k = 1 h duration interval). Therefore, the set, T , has 24 components. The proposed market model can be framed as a bi-level programming problem given by Equations (1)–(16). Equations (1)–(10) and (11)–(16), respectively, represent the upper-level and lower-level problems.
Minimize Ψ n k ( P g , P w ) n N c n ( P g n k ) P g n k + m N 𝓇 λ m k d d m k  
Subject to:
n N P g n k = i N 𝒷 ( 1 χ w ) D i k + ( i , j ) N 𝓁 B b ( ϑ i k ϑ j k ) λ i k S  
B f ( ϑ i k ϑ j k ) F i j   ( i ,   j ) N 𝓁 ,   k
F i j min F i j F i j max ν l k F   ( i ,   j ) N 𝒷
P g n k min P g n k P g n k max ,     { λ g , k S ,   min , λ g , k S ,   max }   n N ,   k
R n d n P g n k P g n k 1 R n u p ,   k ,   n N ,   k
ϑ i min ϑ ϑ i max ,   i
ϑ i = 1 = 0   when   i   is   the   slack   bus
0 < D i k , χ w R ,   k ,   i N 𝒷
where   0 < λ m k d , 0 < d m k   arg   Minimize {  
DRTC = m | ( m , j ) N 𝒶 s m , u ( d m j k ) d m j k  
Subject to:
m | ( m , j ) N 𝒶 A a d m j k = j | ( m , j ) N 𝓇 A r d j m k ,   λ m k d   m ,   j
j N 𝓇 A r d j m k = χ w D i k ,   λ j , k d ,   m ,   k ,   N 𝒷
r m j min d m j k r m j max ,     { λ m j , k d , min , λ m j , k d , max } ,   k ,   ( m j ) A  
| d m j k d m j k 1 |   Δ d m j k ,   k ,   ( m j ) A
0 d m k d m max ,   k ,   m N 𝓊 .   }  

2.1. SCED’s Problem at Upper-Level

A power system network represents an optimal power flow (OPF) model with a number of Nb buses, and Nl transmission lines are considered [25]. In the model, 𝒩𝒷 is the set of busses and 𝒩𝓁 the sets of the transmission lines, respectively. Further, define 𝒩:= {1, 2,…, Ng} is the set for the generation companies. The demand at the ith bus is denoted by Di,∀i𝒩𝒷. The operation cost in (1) consists of two terms. The first term, cn ( P g n k ), refers to the offers cost of the generation P g n k units. The second term, λ m k d , refers to the demand reduction cost for the d m k is for the DR amount provided by aggregator m in time k. The equality constraint (2) represents the grid wise supply-demand balance. The first term on the right-hand side of this constraint is a DR regulating factor. The Bb is for the bus admittance matrix. The term, (ϑikϑjk), is for the voltage phase angles. Representing xij for the reactance of the transmission line from the bus i and j, the diagonal element of the Bb is the sum of all 1/xij. The off-diagonal element of the Bb is equal to the negative of 1/xij, between the bus i and j. The power flow in the lines connecting the bus i and j are represented by Equation (3). The details refer to [25,45,46]. The transmission line and generation supply limits are provided in Equations (4)–(6), respectively. The associated dual variables contracting to the constraints are presented by the double-headed arrow. The ramping rate constraints are given in Equation (6). The constraint Equation (7) imposes voltage phase angle limits. The nonzero decision variables in Equation (10) are found by solving the lower-level DRX optimization model.
The few key assumptions in this model are that the offer pricing scheme in the market clearing model does not consider any source of non-convexities, such as start-up costs of the generating plants or their up/downtime constraints, for simplicity reasons. Also, a direct current approximation of the power system is used where the voltage magnitude in buses of the network is assumed to be unity, the voltage phase differences small, and the transmission lines are lossless, and no security criterion, like outage of transmission lines or generators, is considered.

2.2. DRX’s Problem at Lower-Level

The lower-level problem fixes the optimal DR amount to be traded and the DR transaction cost. The DRX market can be constructed as a directed graph, : = (𝒩, 𝒜), indicated by the node set, 𝒩, composed with an arc set, 𝒜 ⊂ (m, j), joining the trading node m and j as presented in [25]. Two different nodes in L, namely 𝒩𝒶: = {1, 2, …, Na} and 𝒩𝓇: = {1, 2, …, Nr}, corresponding for the DR aggregator and DR purchaser (as see Figure 2). The DR buyers are simplified as those price-sensitive load buses where the unit load reduction save more operation cost.
The dealings between buyers and aggregators in are denoted using the node-arc incidence matrix, A, organized of the vertical concatenation of matrix Aa and Ar. For three aggregators (Na = 3) and two buyers, (Nr = 2), the matrices, Aa and Ar, in Equation (14), take the arrangement as shown in Figure 2.
The quantity of rows in Aa and Ar is equal to the number of aggregators and buyers, respectively. It is assumed the aggregator can sell the DR to several buyers, and each of the DR buyers can buy the DR from several aggregators. The total number of possible deals between the DR buyer and aggregators denote the number of columns of the A matrix. Any DR quantity, dmj (dmj > 0), associated with arc, (m, j) ∈ A, represent the DR transaction from the aggregator node, m𝒩𝒶, to the DR buyer node, j𝒩𝓇. The objective function (13) specifies the DR resource cost function (sm,,u) of the end-user customer, u, under the aggregator, m. The bidding price rises with the growing level DR as represented by:
s m , u = ω m ( α m d m j 2 + β m ( 1 Γ u ) d m j )  
The term, ωm, is the aggregator’s bidding factor, which is used to characterize aggregator price-taking behaviour [26]. The lower the ωm value; the lesser is the price escalation and vice versa. The ωm value can be determined by tracking the LMPs and end-user’s consumption flexibility [47]. The quadratic (αm) and linear (βm) coefficient terms are assumed to be unity. The Γu, ∀i is called end-user type, which takes values between 0 and 1. The non-negativity of type value scales user wiliness to adjust compensation price [48,49,50]. At a higher value of Γu, users are likely to waive compensation for the equal quantity of DR. This intern reduces the DR transaction cost.
For a set of optimal decision variable, dmj (also called primal variable), of length, Na × Nr, the objective function (11) returns a scalar DR transaction cost. The constraint (10) specifies the sum of DR from all aggregators is equal to the total DR provided to DR buses. The coupling constraint (12) refers to the summation of the DR from all the DR buses is as large as the DR limits imposed by the operator.
The term, χw (χw ∈ R), introduced in the upper-level is a regulation parameter for the DR level settings, which limits load demand, Di, at a power system bus, i, i𝒩𝒷. Based on the exogenous variables, such as the adjustment required for the renewables, the total DR from the reference consumption is to be increased or decreased.
The aggregator’s DR collecting behaviour, which can be regulated by the DRXO, may increase the number of end-users by providing a higher compensation rate. The constraint (13) set by the DRXO relates the minimum and maximum DR limits over any transaction, (m, j) ∈ 𝒜. The constraint (14) limits the pick-up/drop-off rate. The constraint (18) specifies that the DR from an aggregator does not exceed its maximum capacity limits, d m m a x .
Assume a vector of dispatch decisions at kth step Ψnk: = [Pg1k, Pg2k, …, PgNgk]’ to represent the generation profile. Define further dispatch decisions considering the DR are ΨnkDR: = [Pg1kDR, Pg2kDR, …, PgNgkDR]’ accordingly. Where the Ψnk to ΨnkDR, respectively, denote generation without and with DR, respectively. The aggregator’s profit (0 < Rm, ∀m) is the difference of a percentage of total cost saving due to the DR and the compensation given to the end-user customers.
R m = γ m ( n , k ) { λ n k S ( Ψ n k ) Ψ n k λ n k S ( Ψ n k DR ) Ψ n k DR } ( m , k ) s m ( d m k k )  
where the λSnknk) and λSnknkDR) are the LMP without and with the DR, respectively. The term in the bracket in Equation (18) corresponds to the total LMP reduction. The term, αm, is for a percentage of monetary reward; aggregator obtains from the operator. The term, γm, equals to unity reflect the aggregator gets paid at a rate of the LMP. The γm > 1 means the payment is made at a rate that is higher than the LMP. This indirectly subsidizes the end-user to increase the incremental curtailment. The γm < 1 refers a partial amount of the LMP is given to limit the DR trading.
Note that the DR bid influences the LMP and consequently reduces the electricity operation cost. For a given value of DR, the actual reduction in the operation cost and the LMP may be different depending on the strategic supply bidding game of the GenCos. The aggregator’s payoff would be cost-effective when the end-user’s DR payment cost is lesser and the LMP is greater.
The GenCo has its own supply cost ($) curve, usually in quadratic form. Indeed, GenCo submits an hourly supply offer at the marginal price ($/MWh). The marginal price is of the form, anPgn + bn, where Pgn is the generation output of unit n, and an, bn are the cost coefficients. It is customary to submit generation blocks, q > 0, ∀q ∈{1, 2, 3, …, Q}; the seller wants to sell at a price, constitute stepwise price-generation block pair expressed by:
c n ( P g n q ) = {   a n 1 ( P g n P g n 1 )   + b n 1 , a n q 1 ( P g n P g n q 1 )   + a n k ( P g n P g n q ) + b n q , P g n P g n 1 b n q 1 ,   P g n 1 < P g n P g n q 1 P g n q < P g n Q .  
The number of blocks and its size rely on the generation capacities. In this study, the quantity index, q, is later on replaced by k to calculate a varying generation profile at kth time step. Table 1 shows the cost coefficients and the generation limit values [26].
In a competitive market environment; the GenCos are a price taker. Hence, during the auction, the supply offer price is made beyond the true marginal cost. The variation of selling offer around the marginal cost is modelled by incorporating a bidding parameter (αn) as shown in Figure 3. The GenCo manipulate αn to increase its revenue when the system load level changes from one critical load level to another.
The solution process of the proposed bi-level model consists of three steps, which are as follows: (Step 1) Implement a Lagrangian duality theorem [51] to convert the lower-level problem by replacing its KKT optimality conditions; (Step 2) consider the decision vector of the lower-level problem as an input parameter for the upper-level problem; and (Step 3) transform the overall problem as a mathematical problem with equilibrium constraints with the technique provided in [52].

3. System Model

The performance of the proposed market model is investigated on a PJM 5-Bus test power system (Figure 4) consisting of three low-cost (G1 and G2 at Bus#1, G5 at Bus#5); two high-cost generation units (G3 at Bus#3, G4 at Bus#4); and the aggregated loads (D2, D3, and D4, at Bus#2, #3, and #4 respectively). The PJM is a regional electric power trading pool market that dispatches the generation and coordinates’ day-ahead capacity and real-time balancing market. It determines the LMP and the generation supply share considering bid-based pricing in a competitive manner. The LMP represents the value of the electricity at the specific location. In the day-ahead capacity market, hourly LMP is determined based on forecasted loads for each hour of the following day. The over or under-estimated loads are adjusted in the real-time balancing market. The test system data, such as generation limits [53] and transmission line parameters [47], are listed in Table 1. The offer prices for all the possible capacities are given in Table 2. The colored numbers in 7th coloumn are the marginal cost (strategic pursuit to achieve higher profit) to the G1 and G4 to deliver an additional unit of electricity at a specific location of the system. The line data, generation capacity, and load demand are at a base of 100 MVA. The upper bound of each aggregator and the DR bidding price for the users under each aggregator are provided in Table 3. The EMO collects the supply offer bid and the load demand data. Further, it settles the generation supply share and the LMP. The LMP also was known as the marginal cost to the EMO to deliver an additional unit of electricity at a specific location of a power system. The LMPs ($/MWh) are Lagrangian multipliers determined at the upper-level. The solution of the model assigns the required system load demand among the generation units in such a way that minimized the operation cost.

4. Numerical Result

In this section, three cases are formulated to validate the efficacy of the proposed market model.
Case#1 discusses the congestion-free and congested conditions and the identification of the critical load buses.
Case#2 discusses the competitive bidding (offering bids at a true marginal price) by GenCos and its impact on LMP and the aggregator’s payoff.
Case#3 discusses the strategic bidding (includes a rebidding at higher price and generation capacity withholding when the system load demands exceed crucial load levels) and its impact at the lower-level market clearing.
To enable the DR market mechanism, the critical load buses need to be identified. The critical load buses reflect those buses where a unit load reduction comparatively saves higher operating cost than the rest of the buses. At this end, a cost sensitivity is investigated for different loading levels. The effect of transmission capacity constraints on locational marginal price (LMP) is analyzed.

4.1. Case#1 (Sensitivity Analysis)

The transmission lines in a power system are planned to operate at safe margins with respect to their capacity limit; congestion only happens at some critical lines and are usually well identified. In the congestion-free condition, the transmission lines are not capacity constrained, i.e., they can transfer any amount of power between the buses, and the given demand is met with no locational variation in LMP. For the given test system, the LMP without transmission line constraints is $30.17/MWh, to serve a 3.00 p.u. of the load at Bus#2, Bus#3, and Bus#4 with a total operating cost of k$12.22. The generator, G1, G2, G3, and G5, provide 1.10, 1.00, 0.90, and 6.00 p.u., respectively. The expensive unit, G4, does not get dispatched in this instance. Non-zero Lagrange multipliers of value $16.58/MWh and $10.83/MWh are found at the Bus#1 and the Bus#5 as generators G1, G2, and G5 hit their upper limits, respectively. The results are summarized in Table 4.
In the congested condition, when line 5–4 is capacity constrained and the maximum power transfer capacity is limited to 2.4 p.u., the market is cleared with high LMP as the cheaper generation cannot be delivered due to the transmission constraint. The LMPs at those five buses are determined to be $28.19/MWh, $32.86/MWh, $30.60/MWh, $31.23/MWh, and $10.60/MWh, as shown in Table 4. The G3 increases from 0.90 to 4.06 p.u., while the G5 reduces from 6.00 to 2.84 p.u. The non-zero Lagrangian multipliers due to hitting the upper limits of G5 are now zero as they are within their limits. The upper limit multipliers for the G1 and G2 still exist. The transmission line constraint allowed the G3 to increase its output and does not hit its upper limit anymore.
The load bus Bus#4 is found to be most sensitive, the Bus#3 is medium, and the Bus#2 is less sensitive. The aggregators are assumed to supply the DR to the most sensitive, Bus#4 and Bus#3. The aggregators bid into the DRX for the incremental DR level to mitigate the LMP spikes. Each of the aggregators has three type/group of users having a different disutility price. The user’s type value defined in aggregator’s DR offer cost function in Equation (18). The DR offer prices are different between two user’s groups. The DR bidding prices are given in Table 3. The offer price continues to rise from the order of low price to high price. However, the bidding parameters are assumed to be constant for a specific DR level.
The DR bidding outcomes determine the transaction cost and compensation price, which can be considered as payback price for the end-users. The load curtailment amount in the DRX market is regulated by the LMP determined in the upper-level market. The aggregators get paid at the LMP rate of those buses where they provide the DR. The aggregator net payoff comprises of the difference of the LMP rate they get paid and the compensation cost provided to the end-users. The DRX market is considered during the peak demand hours, from 6 am to 9 am in the morning and from 4 pm to 8 pm in the evening.

4.2. Case#2 (Competitive Bidding)

Operation Cost and LMPs: The market model is simulated twice, without and with a DRX mechanism, to quantify the benefits of reduced LMP and operating cost. The operation cost without a DRX is determined to be k$747.49 and an average (over a day) LMP at Bus#4 is found to be $56.49/MWh. The reduction of the LMP ($49.14/MWh LMP at Bus#4) and the operation cost (k$724.16) is evident even with a modest amount of incremental DR (1.95%) for few hours. The effect of different levels of DR on the operating cost with, without DRX transaction cost, is summarized in Table 5 and Figure 5. The operation cost without considering DRX gradually decreases due to serving a reduced load demand. It is interesting to note that with DRX, the operation cost reduces with the DR up to 15.28% and increases for any further increase in DR levels. This is because, at higher DR levels, the DR transaction cost outweighs savings from the reduction in the market clearing cost (the DR transaction cost and the DR amount among the aggregators and users are the decision variables determined by solving the market at the lower-level). The DR transaction cost is reasonably small for up to 12.94% DR and increases sharply beyond it. The DR transaction cost progressively increases with the increasing level of DR. A larger value of DR yields a higher selling revenue for the DRX participants. The DR transaction cost at the lower-level equivalently provides monetary benefit for the end-user customers. The end-user proportionately shares this based on the disutility price they offered. However, the increased DR transaction cost exceeds the overall operation cost provide an argument for limiting the DR payments to periods when the LMP are likely to exceed a specific threshold. The emission without DRX was 19,077 ton a day. With 1.95% and 5.92% DR, the emission reduced to 18,806 and 18,167 ton, respectively. A 1.42% and 6.88% emission, respectively, reduces in each day.
Table 6 presents average LMP at different buses in the test system. It is observed that the LMP at Bus#4 is highest and at Bus#5 is lowest (the lower LMP reflects the existence of low-cost generation units to meet demand). The LMP without DR was found to be maximum at Bus#4, and minimum at Bus#5. The LMP at Bus#3 and #4 without DR were $55.25/MWh, and $56.49/MWh, respectively. At 1.95%DR, it reduces to $52.30/MWh and $53.46/MWh, respectively. With the increased DR, the LMP reduces, however, till a certain level and does not further reduce. The LMP was found minimal when the DR level 10.60%. The LMPs, usually, are the basis for payment to generators, aggregators, and payments by the retailers.
The GenCos and the aggregators are paid at the LMP rate at their respective buses. The radar charts in Figure 6 compare the hourly LMP variation for a day. The radial axis starting from the centre shows the LMP value, while the equiangular distance along the peripheries shows an hour of a day. The hourly LMP is depicted by the marker on the radial axis. The individual line colour represents the LMP with different levels of DR.
The non-peak hour LMP at Bus#2, #3, and #4 does not significantly change. Further, the LMP at the Bus#5 only changes between $10.76 and $24.01. In Bus#5, the unit G5 get paid as its bid. These values are almost the same irrespective of the DR level. At Bus#2, #3, and #4, the LMPs of the peak demand period are responsive to the incremental DR level. The peak hour LMP at Bus#4, without DR, is found to be $82.67, which reduces to $67.45 at 10.60% DR. At Bus#3, those values are $80.85 and $66.10. At Bus#2, the LMPs are $75.86 and $62.40. Increasing the DR level decreases the peak hour durations. With a DR level of up to 5.92%, the duration of the LMP spike reduces to 4 h in a day and is further reduced to only 3 h for DR level of 8.26% and can be totally avoided with the DR level of 10.60%.
Aggregators’ Payoff: The lower-level model allows the aggregators to participate in the DR exchange in the same way that supply-side Genco’s bid into day-ahead electricity markets. The proposed DRX integrated market clearing model promotes the true market value of DR in daily scheduling intervals. Figure 7 illustrates the aggregators’ payoff with different incremental DR levels and Table 7 presents the DR amount supplied by the aggregators. The aggregator’s payoff not only depends on the amount of DR traded, but also on the LMP. The aggregator payoff is a difference of the DR selling revenue at LMP, and the compensation provided to end-users and is determined using Equation (18). Assuming, the reward scale factor, γm, as unity, with a 1.95% DR, the aggregator, A1, A2, and A3, earn k$15.15, k$13.74, and k$14.41, offering the DR amounts of 2.69, 2.66, and 2.66, respectively, in the DR trading period in a day. Irrespective of the equal DR provided by both the J2 and J3, the payoff difference happens due to different compensation price among the user group. The payoff rises and becomes maximum at the 8.26% DR level. After a 12.94% DR, a payoff reduction is observed, because of the decrease in the LMP and increase in the rate of compensation price. At the 15.28% DR level, A1 becomes marginally profitable and both A2, A3 outweigh compensation paid to the end-users than the LMP at which the A2 and A3 get paid. For any DR levels greater than 17.72%, all the aggregators lose. The polynomial fit (second order) of the payoff is also shown in the figure. As observed, with the increasing DR level, the payoff increases, reaches a maximum value, then decreases before finally becoming zero, at which point, the EMO also loses. Until such a critical DR level, the DR payoff resulting from new revenue generation driven by DRX market is considerable.
DR and Compensation Share among the End-Users: In the proposed model, aggregators compensate the users for the DR that they provide. The transaction cost and compensation price are determined in DRX markets. The magnitude of possible DR transaction cost depends on how much the users can afford and the marginal disutility price. Figure 8 shows a disutility price for different user groups. Three types/groups of users under each aggregator are considered. The disutility price for user group, U1, under the aggregator, A1, is the stepper. The DR price is smallest for the U1, while it is highest for the U3. The user, U1, under the A2 has minimum disutility ($21.71), while a maximum for the group, U3 ($26.43). The offer price for the U1 under aggregator, A3, is $19.2. The price for the rest of the two group, U2 and U3, is around $24.60. The aggregator, A1, has a maximum DR supply capacity of 2 p.u. and A2 and A3 whereas having the capacity of 1.75 p.u. The end-users get paid by the aggregator at a compensation price, which is a dual multiplier associated with Equation (12) of the lower-level problem. Due to receiving compensation (payback), the users reduce the energy consumption cost. However, from the EMO perspective, the DR benefit can be obtained if the operation cost reduction at the upper-level does not outweigh the compensation given to the end-users at the lower-level. This possibly occurs when the DR demand rises significantly higher if the user seeks higher compensation. Multiplying the DR price with the allocated DR is regarded as overall DR compensation benefit for the users.
Figure 9 presents the optimal DR amount provided by the user groups for each incremental DR level. The DR transaction cost, the amount supplied, and the compensation price is found by solving the lower-level problem. Until a 10.60% DR, the user group, U2, under the aggregators do not provide any DR amount. From 10.60 to 15.28% DR; both the U1 and U2 deliver with a varying amount. At a 17.62% DR and onward, all the user group under the aggregator, A3, contribute to supply.
The DR compensation is given in Table 8. At the 8.26% DR level, the compensation benefit provided to the user groups are determined to be k$9.06, k$10.76, and k$9.52, respectively. The DR transaction cost at 5.92%DR is $29349. Next, at 8.26% DR, transaction cost increases to $48762. A k$19.41 increase due to 2.34% DR level. The cost increases sharply as the DR requirement rises. At 17.72% DR, the transaction cost increases almost three times than that of the cost at 8.26% DR. The higher transaction cost due to rising DR compensation offer price is neither profitable for the aggregator nor to the EMO. As seen from the DRX integrated operation curve (Figure 5), the DRX integrated market cost gradually decreases with the incremental DR level. The cost reaches a minimum of 8.26% DR. Further, it started to increase, and about 15.28% DR, the DRX integrated operation cost become equal to the cost without DRX. From the economic point of view, the DR integration level should be less than 15.28% DR.
The payoff margin of the aggregators varies with the DR level (Figure 7) and the trend line shows that the maximum payoff is achieved at the 8.26% DR level. The payoff reduces as the DR level rises. At 15.28% DR, the payoff is marginal for the A1, while for the A2 and A3, DR transaction is profitable. Now, if the LMP scaled parameter, γm, is increased, then the payoff margin also increases to transact higher amount of DR, provided end-users kept their compensation price fixed.
Generation DispatchShare: The solution of the upper-level economic market clearing problem finds the optimal generation dispatch. The G1, G2, and G5 bid a lower price than others. In the market clearing model, the generation units’ output is arranged according to the transmission security and economy-based merit order. Figure 10, Figure 11 and Figure 12 compare the optimal generation dispatch without and with DRX (with two different DR level), respectively. Without DRX, the generation unit, G1 and G2, are at their maximum capacity, while the G3 and G4 change within its generation limits across the scheduling hour (as in Figure 10). The amount supplied by G5 is restricted by line capacity constraint, thereby, most of the time, its output changes around 2.85 p.u. The dotted line in red colour indicates the total load demand served by all the generators.
With the proposed DRX model and with a 1.95% DR participation, the total demand without and with DR is shown by the red and green colour in Figure 11. The DR amount, which is the difference between those two curves, becomes higher during the peak period. It is observed that the G4 is not dispatched at hours 7 am to 9 am as peak loads are curtailed. The G4 dispatches a few hours around evening peak. Further, in Figure 12, with 5.92% DR, the G4 only dispatches around evening peak periods. Such dispatch results are economical, since the most expensive unit got restriction. The reduced dispatch of G4 is compensated mainly by the second least expensive unit, G3, to meet the load.

4.3. Case#3 (Strategic Bidding)

In strategic bidding, the GenCo bids at a higher price by price distortion or generation capacity withholding deliberately. This is not a routine practice, rather adopt this strategy when the system load demand during the peak hours exceeds significant load levels. It would argue that such opportunistic behaviour cannot always be perceived as profitable [55]. Usually, the GenCos are prone to exercise market power and try to increase their revenue, eventually yielding to a higher overall operating cost [56]. The strategic bidding was mostly attributed to the expensive units. The impact on dispatch share, LMPs, operation cost, and aggregator’ payoff is investigated considering the following scenarios:
Scenario#1: The unit, G3, bids $83.85 instead of $66.10 in the hour 7 pm, 8 pm, and 9 pm. The unit, G4, and others bid as completive.
Scenario#2: The unit, G4, bids $103.92 instead of $82.67 in those hours indicated in Scenario#1. The unit, G3, and others bid as completive.
Scenario#3: Both the G3 and G4 bid simultaneously with the new offer price, while others bid as completive.
Therefore, Scenario#1, Scenario#2, and Scenario#3, respectively, consider strategic bidding of the unit G3 only, the unit G4 only, and both the G3 and G4 simultaneously. The rest of the conditions, like loading level, transmission constraints, and the capacity of the rest of the units, are the same. The impact on dispatch share, LMPs, and operation cost are investigated.
Table 9 shows the impact of GenCos strategic bidding on aggregator’s payoff. The payoff set of each aggregator at different DR levels is shown in a second and third row. The sum of the payoff for all aggregators is presented in the fourth and fifth row. At 5.92% DR, the payoff maximum was accounted to be k$67.84 in Scenario#1, while minimum at competitive bidding determines to be k$62.99. A relative payoff variation was determined to be $k2.2. In the same DR level at Scenario#3, the relative increase payoff rises to k$7.69. It is observed that at 8.26% DR, the payoff maximum accounted to be k$79.01 in Scenario#3, while minimum at competitive bidding determines to be k$62.69. A relative payoff variation is $k20.03. As seen, the corresponding payoff for the Scenario#1 and Scenario#2, comparison to the competitive case is almost the same accounted to be k$17.33. In conclusion, GenCos strategic bidding enhances aggregator’s payoff, the payoff is comparatively higher if the number of GenCo practice strategic bidding becomes higher.
The supply offer prices for strategic bidding are highlighted in Table 2. Table 10 compares the operation cost and the LMPs for competitive and strategic bidding. The operation costs at 5.92% and 8.26% DR are shown in this table. The second rows show the cost without DR, while the third and fourth are for the 5.92% and 8.26% DR. In Scenario#1, when the G3 is strategic, the operation cost rises to k$764.99. In Scenario#2 and #3, the cost is k$753.61 and k$772.08, respectively. The relative decrease in operating cost compared to the competitive case found to be maximum at Scenario#3 is 3.69%.
The strategic bidding must clear the market with a higher price. The LMP average over a day in Bus#4 without the DR for the Scenario#1, Scenario#2, and Scenario#3 was found to be $54.98, $56.49, and $56.49, respectively. The DR reduces the price anyway. In Scenario#3, with a 5.92% DR, the average LMP at Bus#4 was found to be $52.57, which was $50.92 in competitive Case#2. At this DR level, the strategic bidding increases the LMP $1.67/MWh.
In summary, the Scenario#3 is the worst-case bidding, which escalates the operating cost to a maximum, keeping the dispatch share the same as that of the competitive case. However, the aggregator’s payoff is higher when a higher number of generators adopt strategic bidding.

5. Conclusions

This paper presents a market framework to effectively integrate DR and reduce wholesale electricity prices. An exemplary DRX market for three aggregators, each having three types of customers, has been modeled. In the model, the aggregators can provide predictable DR services at short notice, which is manageable in real-time. Several plausible bidding behaviors of the market participants have been considered.
Simulation results compare the market clearing price, GenCos’ revenue, aggregator’s payoff, and the amount and cost of the transacted DR. Both competitive and strategic bidding scenarios with different levels of strategies were investigated. The proposed DRX scheme ensures that the least cost operation is achieved either by increasing DR when the GenCos bid strategically or decreasing the amount of DR when the GenCos bid competitively. In competitive bidding, the operation cost is the lowest, which means the DR amount can be reduced, resulting in reduced customer inconvenience. As a price maker, this strategy pushes the market clearing price higher.
The aggregators are price takers when the GenCos exercise strategic bidding. For the same amount of DR, the aggregators get a higher payoff, which would not have been possible under competitive bidding. The DR was found to be profitable as long as its transaction is cost-effective and economical. However, if the GenCos do not return from the price escalation strategy, they must reduce their dispatch share and will likely face a “missing money” problem where the dispatched energy from the GenCos is less and even reduces further.
In conclusion, the DR can be used to exclude GenCos price-lifting behaviour to some extent while maintaining the minimal energy consumption cost for the end-users. Moreover, the reduction in the overall operating cost, carbon emission, and congestion cost were achieved. However, beyond a critical DR level, the DR compensation price becomes higher, outweighing the benefits.
This paper focuses mainly on DR trading at the wholesale level and a detailed distribution market framework could be explored to allow aggregators to trade DR. Collecting DR from the residential customer under different DR programs in the distribution electricity market should be investigated in the future.

Author Contributions

N.M. established the model, implemented the simulation, and wrote this article; Y.M. guided the research, revised the paper, and refined the language.

Acknowledgments

This paper is one of the chapters of Nur’s PhD thesis. Authors pleased to acknowledge the scholarship support Nur received from Queensland University of Technology during his PhD tenure. Authors remember Mr. Paul McArdle, Managing Director of Global Roam Pty Ltd, for sharing his ideas on the role of demand response aggregators in competitive electricity markets. Heartiest gratitude to Professor Gerard Ledwich, who helped us to get updates on innovations and opportunities in demand-side management (DSM) by organizing a series of seminars.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

𝒩𝒶: m, Naset, index, number of aggregators
𝒩𝓇, j, Nrset, index and number of DR buyer
𝒩𝓁, Ngnumber and set of generation companies (GenCos)
𝒩, i, Nbset, index and number of power system buses
𝒩𝒷, l, Nlset, index and number of transmission lines
Bbmatrix of dimension Nb ×Nb for admittance
ϑikvoltage angle at bus i in time k
αm, βmDR offer cost coefficient for aggregator
ωm ∈ (0 ~ 1)aggregator bidding parameter/end-user type
an, bncost-coefficient for the GenCos n𝒩𝓁
Δdmjramp-down limit of the DR provided by m𝒩𝒶
d m k DR amount provided by aggregator m
d m k m a x ,   d m k m i n limit of DR provided by aggregator m
P g n k m a x , P g n k m i n limit of generation amount provided by GenCo
F i j k power flow from the bus i to j
F i j m a x , F i j m a x line maximum capacity limit between the bus i and j
R n d n ramp-down limit of generator
R n u p ramp-up limit of generator n
Dikload demand at power system bus i
dmjDR sells by aggregator
λ j k D Lagrangian multipliers of upper-level problem
λ m k d Lagrangian multipliers of lower-level problem
χw ∈ Rthe DR tunng parameter to deal wind-generated power variability
cn ( P g n k )cost of the GenCos, n
sm,u( · )aggregated DR selling offer cost of m
Ψ n k D ( Ψ n k ) LMP without DR at bus n, in time k [$/MWh]
λ n k D ( Ψ n k D R ) LMP with DR at bus n, in time k [$/MWh]

References

  1. Fang, X.; Hu, Q.; Li, F.; Wang, B.; Li, Y. Coupon-Based Demand Response Considering Wind Power Uncertainty: A Strategic Bidding Model for Load Serving Entities. IEEE Trans. Power Syst. 2015, 31, 1025–1037. [Google Scholar] [CrossRef]
  2. Faruqui, A.; George, S. Quantifying customer response to dynamic pricing. Electr. J. 2005, 18, 53–63. [Google Scholar] [CrossRef]
  3. Mohammad, N.; Mishra, Y. Transactive Market Clearing Model with Coordinated Integration of Large-Scale Solar PV Farms and Demand Response Capable Loads. In Proceedings of the 2017 Australasian Universities Power Engineering Conference, AUPEC, Melbourne, Australia, 21–23 November 2017; pp. 1–6. [Google Scholar]
  4. Vardakas, J.S.; Zorba, N.; Verikoukis, C.V. A Survey on Demand Response Programs in Smart Grids: Pricing Methods and Optimization Algorithms. IEEE Commun. Surv. Tutor. 2015, 17, 152–178. [Google Scholar] [CrossRef]
  5. Deng, R.; Yang, Z.; Chow, M. A Survey on Demand Response in Smart Grids. IEEE Trans. Ind. Inform. 2015, 11, 1–8. [Google Scholar] [CrossRef]
  6. Albadi, M.H.; El-Saadany, E.F. A summary of demand response in electricity markets. Electr. Power Syst. Res. 2008, 78, 1989–1996. [Google Scholar] [CrossRef]
  7. Dave, S.; Sooriyabandara, M.; Zhang, L. Application of a game-theoretic energy management algorithm in a hybrid predictive-adaptive scenario. In Proceedings of the 2nd IEEE PES International Conference and Exhibition on Innovative Smart Grid Technologies (ISGT Europe), Manchester, UK, 5–7 December 2011. [Google Scholar] [CrossRef]
  8. Märkle-Huß, J.; Feuerriegel, S.; Neumann, D. Large-scale demand response and its implications for spot prices, load and policies: Insights from the German-Austrian electricity market. Appl. Energy 2018, 210, 1290–1298. [Google Scholar] [CrossRef]
  9. Kiani, A.; Annaswamy, A. Wholesale energy market in a smart grid: Dynamic modeling and stability. In Proceedings of the 2011 50th IEEE Conference on Decision and Control and European Control Conference, Orlando, FL, USA, 12–15 December 2011; pp. 2202–2207. [Google Scholar]
  10. Ben, F.; David Burns, J.G.D.K.; Samin, M.P.L.; Peirovi, S. Demand Response and Advanced Metering. Available online: https://www.ferc.gov/legal/staff-reports/2017/DR-AM-Report2017.pdf (accessed on 5 August 2018).
  11. Asimakopoulou, G.E.; Vlachos, A.G.; Hatziargyriou, N.D. Hierarchical Decision Making for Aggregated Energy Management of Distributed Resources. Power Syst. IEEE Trans. 2015, 30, 3255–3264. [Google Scholar] [CrossRef]
  12. Khan, A.A.; Razzaq, S.; Khan, A.; Khursheed, F. Owais HEMSs and enabled demand response in electricity market: An overview. Renew. Sustain. Energy Rev. 2015, 42, 773–785. [Google Scholar] [CrossRef]
  13. Setlhaolo, D.; Xia, X. Optimal scheduling of household appliances with a battery storage system and coordination. Energy Build. 2015, 94, 61–70. [Google Scholar] [CrossRef]
  14. Zhao, Q.; Shen, Y.; Li, M. Control and Bidding Strategy for Virtual Power Plants with Renewable Generation and Inelastic Demand in Electricity Markets. IEEE Trans. Sustain. Energy 2016, 7, 562–575. [Google Scholar] [CrossRef]
  15. Paterakis, N.G.; Erdinc, O.; Bakirtzis, A.; Catalao, J.P. Optimal Household Appliances Scheduling under Day-Ahead Pricing and Load-Shaping Demand Response Strategies. IEEE Trans. Ind. Inform. 2015, 11, 1509–1519. [Google Scholar] [CrossRef]
  16. Shariatzadeh, F.; Mandal, P.; Srivastava, A.K. Demand response for sustainable energy systems: A review, application and implementation strategy. Renew. Sustain. Energy Rev. 2015, 45, 343–350. [Google Scholar] [CrossRef]
  17. Nunna, K.H.S.V.S.; Doolla, S. Responsive end-user-based demand side management in multimicrogrid environment. IEEE Trans. Ind. Inform. 2014, 10, 1262–1272. [Google Scholar] [CrossRef]
  18. Cappers, P.; Goldman, C.; Kathan, D. Demand response in U.S. electricity markets: Empirical evidence. Energy 2010, 35, 1526–1535. [Google Scholar] [CrossRef] [Green Version]
  19. Vivekananthan, C.; Mishra, Y.; Ledwich, G.; Li, F. Demand response for residential appliances via customer reward scheme. IEEE Trans. Smart Grid 2014, 5, 809–820. [Google Scholar] [CrossRef]
  20. Wang, Q.; Zhang, C.; Ding, Y.; Xydis, G.; Wang, J.; Østergaard, J. Review of real-time electricity markets for integrating Distributed Energy Resources and Demand Response. Appl. Energy 2015, 138, 695–706. [Google Scholar] [CrossRef]
  21. Parvania, M.; Fotuhi-Firuzabad, M. Integrating load reduction into wholesale energy market with application to wind power integration. IEEE Syst. J. 2012, 6, 35–45. [Google Scholar] [CrossRef]
  22. Gkatzikis, L.; Koutsopoulos, I.; Salonidis, T. The role of aggregators in smart grid demand response markets. IEEE J. Sel. Areas Commun. 2013, 31, 1247–1257. [Google Scholar] [CrossRef]
  23. Ali, M.; Alahäivälä, A.; Malik, F.; Humayun, M.; Safdarian, A.; Lehtonen, M. A market-oriented hierarchical framework for residential demand response. Int. J. Electr. Power Energy Syst. 2015, 69, 257–263. [Google Scholar] [CrossRef]
  24. Valinejad, J.; Barforoshi, T.; Marzband, M.; Pouresmaeil, E.; Godina, R.; Catalão, J.P.S. Investment Incentives in Competitive Electricity Markets. Appl. Sci. 2018, 8, 1978. [Google Scholar] [CrossRef]
  25. Mohammad, N.; Mishra, Y. Coordination of wind generation and demand response to minimise operation cost in day-ahead electricity markets using bi-level optimisation framework. IET Gener. Transm. Distrib. 2018, 12, 3793–3802. [Google Scholar] [CrossRef]
  26. Mohammad, N.; Mishra, Y. Competition Driven Bi-Level Supply Offer Strategies in Day Ahead Electricity Market. In Proceedings of the Australasian Universities Power Engineering Conference, Brisbane, Australia, 25–28 September 2016; pp. 1–6. [Google Scholar]
  27. Mohammad, N.; Mishra, Y. Book Chapter, Demand Side Management and Demand Response for Smart Grid. In Handbook of Smart Grid Communication Systems; Kabalci, K., Ersan, Y., Eds.; Springier Nature: Singapore, 2019; pp. 197–231. ISBN 978-981-13-1768-2. [Google Scholar]
  28. Xu, Y.; Li, N.; Low, S.H. Demand Response with Capacity Constrained Supply Function Bidding. IEEE Trans. Power Syst. 2015, 31, 1–12. [Google Scholar] [CrossRef]
  29. Nguyen, D.T.; Nguyen, H.T.; Le, L.B. Dynamic Pricing Design for Demand Response Integration in Power Distribution Networks. IEEE Trans. Power Syst. 2016, 31, 3457–3472. [Google Scholar] [CrossRef]
  30. Saez-Gallego, J.; Kohansal, M.; Sadeghi-Mobarakeh, A.; Morales, J.M. Optimal Price-Energy Demand Bids for Aggregate Price-Responsive Loads. IEEE Trans. Smart Grid 2017, 9, 5005–5013. [Google Scholar] [CrossRef]
  31. Nguyen, T.; Negnevitsky, M.; Groot, M. De Pool-based Demand Response Exchange: Concept and modeling. IEEE Trans. Power Syst. 2011, 26, 1677–1685. [Google Scholar] [CrossRef]
  32. Shafie-Khah, M.; Shoreh, M.H.; Siano, P.; Fitiwi, D.Z.; Godina, R.; Osorio, G.J.; Lujano-Rojas, J.; Catalao, J.P.S. Optimal Demand Response Programs for improving the efficiency of day-ahead electricity markets using a multi attribute decision making approach. In Proceedings of the 2016 IEEE International Energy Conference, ENERGYCON, Leuven, Belgium, 4–8 April 2016. [Google Scholar]
  33. Ayón, X.; Gruber, J.K.; Hayes, B.P.; Usaola, J.; Prodanović, M. An optimal day-ahead load scheduling approach based on the flexibility of aggregate demands. Appl. Energy 2017, 198, 1–11. [Google Scholar] [CrossRef]
  34. Nan, S.; Zhou, M.; Li, G. Optimal residential community demand response scheduling in smart grid. Appl. Energy 2018, 210, 1280–1289. [Google Scholar] [CrossRef]
  35. Ampimah, B.C.; Sun, M.; Han, D.; Wang, X. Optimizing sheddable and shiftable residential electricity consumption by incentivized peak and off-peak credit function approach. Appl. Energy 2018, 210, 1299–1309. [Google Scholar] [CrossRef]
  36. Shoreh, M.H.; Siano, P.; Shafie-khah, M.; Loia, V.; Catalão, J.P.S. A survey of industrial applications of Demand Response. Electr. Power Syst. Res. 2016, 141, 31–49. [Google Scholar] [CrossRef]
  37. Zimmerman, R.D.; Murillo-Sanchez, C.E.; Thomas, R.J. Matpower: Steady-State Operations, Planning, and Analysis Tools for Power Systems Research and Education. IEEE Trans. Power Syst. 2011, 26, 12–19. [Google Scholar] [CrossRef]
  38. Sood, Y.R.; Padhy, N.P.; Gupta, H.O. Deregulated model and locational marginal pricing. Electr. Power Syst. Res. 2007, 77, 574–582. [Google Scholar] [CrossRef]
  39. 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]
  40. Wang, B.; Fang, X.; Zhao, X.; Chen, H. Bi-level optimization for available transfer capability evaluation in deregulated electricity market. Energies 2015, 8, 13344–13360. [Google Scholar] [CrossRef]
  41. Fernandez-Blanco, R.; Arroyo, J.M.; Alguacil, N.; Guan, X. Incorporating Price-Responsive Demand in Energy Scheduling Based on Consumer Payment Minimization. IEEE Trans. Smart Grid 2016, 7, 817–826. [Google Scholar] [CrossRef]
  42. Kardakos, E.G.; Simoglou, C.K.; Bakirtzis, A.G. Optimal Offering Strategy of a Virtual Power Plant: A Stochastic Bi-Level Approach. IEEE Trans. Smart Grid 2016, 7, 794–806. [Google Scholar] [CrossRef]
  43. Zugno, M.; Morales, J.M.; Pinson, P.; Madsen, H. A bilevel model for electricity retailers’ participation in a demand response market environment. Energy Econ. 2013, 36, 182–197. [Google Scholar] [CrossRef]
  44. Zhao, H.; Wang, Y.; Guo, S.; Zhao, M.; Zhang, C. Application of a Gradient Descent Continuous Actor-Critic Algorithm for Double-Side Day-Ahead Electricity Market Modeling. Energies 2016, 9, 725. [Google Scholar] [CrossRef]
  45. Bertsekas, D. Network Optimization: Continuous and Discrete Models; Athena Scientific: Belmont, MA, USA, 1998; ISBN 1886529027. [Google Scholar]
  46. Pasqualetti, F.; Bicchi, A.; Bullo, F. A graph-theoretical characterization of power network vulnerabilities. In Proceedings of the 2011 American Control Conference, San Francisco, CA, USA, 29 June–1 July 2011; pp. 3918–3923. [Google Scholar]
  47. Papavasiliou, A.; Oren, S.S. Large-Scale integration of deferrable demand and renewable energy sources. IEEE Trans. Power Syst. 2014, 29, 489–499. [Google Scholar] [CrossRef]
  48. Rahimiyan, M.; Baringo, L. Strategic Bidding for a Virtual Power Plant in the Day-Ahead and Real-Time Markets: A Price-Taker Robust Optimization Approach. IEEE Trans. Power Syst. 2015, 31, 2676–2687. [Google Scholar] [CrossRef]
  49. Ghasemi, A.; Mortazavi, S.S.; Mashhour, E. Hourly demand response and battery energy storage for imbalance reduction of smart distribution company embedded with electric vehicles and wind farms. Renew. Energy 2016, 85, 124–136. [Google Scholar] [CrossRef]
  50. Wu, H.; Shahidehpour, M.; Alabdulwahab, A. Demand Response Exchange in the Stochastic Day-Ahead Scheduling With Variable Renewable Generation. IEEE Trans. Sustain. Energy 2015, 6, 516–525. [Google Scholar] [CrossRef]
  51. Conejo, A.J.; Castillo, E.; Minguez, R.; Garcia-Bertrand, R. Decomposition Techniques in Mathematical Programming: Engineering and Science Applications; Springer: Dordrecht, Netherlands, 2008; ISBN 9783540276852. [Google Scholar]
  52. Fortuny-amat, A.J.; Mccarl, B.; Fortuny-amat, J.; Mccarl, B. A Representation and Economic Interpretation of a Two-Level Programming Problem. J. Oper. Res. Soc. 1981, 32, 783–792. [Google Scholar] [CrossRef]
  53. Hao, H.; Lin, Y.; Kowli, A.S.; Barooah, P.; Meyn, S. Ancillary Service to the grid through control of fans in commercial Building HVAC systems. IEEE Trans. Smart Grid 2014, 5, 2066–2074. [Google Scholar] [CrossRef]
  54. Abdollahi, A.; Parsa Moghaddam, M.; Rashidinejad, M.; Sheikh-El-Eslami, M.K. Investigation of economic and environmental-driven demand response measures incorporating UC. IEEE Trans. Smart Grid 2012, 3, 12–25. [Google Scholar] [CrossRef]
  55. Guan, X.; Ho, Y.C.; Pepyne, D.L. Gaming and price spikes in electric power markets. IEEE Trans. Power Syst. 2001, 16, 402–408. [Google Scholar] [CrossRef]
  56. Mohammad, N. Competitive Demand Response Trading in Electricity Markets: Aggregator and End–user Perspectives. PhD Thesis, Queensland University of Technology, 2018. Accession No. 119702. pp. 1–208. [Google Scholar] [CrossRef]
Figure 1. The schematic of the DR trading process among market participants.
Figure 1. The schematic of the DR trading process among market participants.
Energies 11 03296 g001
Figure 2. The DR transaction mechanism, (a) a trading network presenting two DR procuring points (Nr = 2) and three aggregators (Na = 3); (b) a segmentation of the A matrix indicating the DR aggregator node, Na, and the DR procuring node, Nr.
Figure 2. The DR transaction mechanism, (a) a trading network presenting two DR procuring points (Nr = 2) and three aggregators (Na = 3); (b) a segmentation of the A matrix indicating the DR aggregator node, Na, and the DR procuring node, Nr.
Energies 11 03296 g002
Figure 3. Marginal cost function and bid curve for the GenCo.
Figure 3. Marginal cost function and bid curve for the GenCo.
Energies 11 03296 g003
Figure 4. The 5-Bus PJM power network with aggregated loads are at Bus#2 and #3 and Bus#4. The line supply data are provided in Table 1.
Figure 4. The 5-Bus PJM power network with aggregated loads are at Bus#2 and #3 and Bus#4. The line supply data are provided in Table 1.
Energies 11 03296 g004
Figure 5. The operation cost with and without DR transaction.
Figure 5. The operation cost with and without DR transaction.
Energies 11 03296 g005
Figure 6. The effect of incremental DR levels on the LMP at different buses.
Figure 6. The effect of incremental DR levels on the LMP at different buses.
Energies 11 03296 g006aEnergies 11 03296 g006b
Figure 7. Aggregator’s payoff across the incremental DR.
Figure 7. Aggregator’s payoff across the incremental DR.
Energies 11 03296 g007
Figure 8. The baseline DR supply offer bidding prices, which increase if the DR rises.
Figure 8. The baseline DR supply offer bidding prices, which increase if the DR rises.
Energies 11 03296 g008
Figure 9. The DR provided by the users under each aggregator.
Figure 9. The DR provided by the users under each aggregator.
Energies 11 03296 g009
Figure 10. Generation dispatch mix among the units without DRX assuming no elasticity to the demand.
Figure 10. Generation dispatch mix among the units without DRX assuming no elasticity to the demand.
Energies 11 03296 g010
Figure 11. Generation dispatch mix in the DRX integrated market clearing with 1.95% DR participation.
Figure 11. Generation dispatch mix in the DRX integrated market clearing with 1.95% DR participation.
Energies 11 03296 g011
Figure 12. Generation dispatch mix in the DRX integrated market clearing with 5.92% DR participation.
Figure 12. Generation dispatch mix in the DRX integrated market clearing with 5.92% DR participation.
Energies 11 03296 g012
Table 1. Generation capacity limits, emission rate, lines, and load data.
Table 1. Generation capacity limits, emission rate, lines, and load data.
Gen (Fuel type)Capacity Limits[Pnmin Pnmax ] (p.u.)Carbon Emissions (kgCO2/kWh) [54]
G1 (Coal), Bus#1[0.25, 1.10]0.940
G2 (Diesel), Bus#1[0.25, 1.00]0.778
G3 (Coal), Bus#3[1.50, 5.20]0.940
G4 (Gas), Bus#4[0.50, 3.00]0.581
G5 (Coal), Bus#5[2.25, 6.00]0.940
Transmission Lines and Load Data
From Bus iTo Bus jReactance (xij)Capacity Limit [Fijmin Fijmax]Load Demand Di (p.u.)
122.81[−8.75, 8.75]Bus#2 = 3.00
231.08[−8.75, 8.75]Bus#3 = 3.00
432.97[−8.75, 8.75]Bus#4 = 3.00
542.97[−8.75, 2.40]
510.64[−8.75, 8.75]
143.04[−8.75, 8.75]
Table 2. The cost coefficient, bidding quantities, and prices for generation units.
Table 2. The cost coefficient, bidding quantities, and prices for generation units.
GenBlockGeneration Supply Offer Price
[an, bn]Size[1st Block][2nd Block][3rd Block][4th Block][5th Block]
($/MWh, $/MWh2)(MWh)
G1 [0.115, 11.50]50.0011.61, 0.517.36, 1.123.1128.8634.61
G2 [0.115, 11.50]50.0011.61, 0.517.36, 1.023.1128.8634.61
G3 [0.355, 12.50]100.012.85, 1.030.60, 2.048.35, 3.066.10, 4.083.85, 5.0
G4 [0.425, 18.50]60.0018.92, 0.640.17, 1.261.42, 1.882.67, 2.4103.92, 3.0
G5 [0.265, 10.50]120.010.76, 1.224.01, 2.437.26,3.650.51, 4.863.76, 6.0
Table 3. The DR capacity, cost coefficient, and bidding prices.
Table 3. The DR capacity, cost coefficient, and bidding prices.
AggregatorCapacityDR Cost CoefficientsThe DR Offer Price Segments for the User Group at ωm = 1
dmaxαmβmU1U2U3
A11.90048512.1516.9723.0627.42
A21.700.09210.3521.7124.7226.43
A31.700.06811.0519.2024.6524.66
The bidding parameter, ωm in Equation (18), is varied from 1:5.50 with a different increment to rescaled price coefficient, αm and βm, for the different DR levels.
Table 4. The generator output and the LMPs to serve load {D2, D3, D4} = {3.00, 3.00, 3.00} p.u.
Table 4. The generator output and the LMPs to serve load {D2, D3, D4} = {3.00, 3.00, 3.00} p.u.
Scenarios{G1, G2, G3, G4, G5}LMPBus # ii} ($/MWh)Total Cost (k$)
Without Line Constraints{1.10, 1.00, 0.90, 0.00, 6.00}{30.17, 30.17, 30.17, 30.17, 30.17}12.22
Line 5–4 is Constrained{1.10, 1.00, 4.06, 0.00, 2.84}{28.19, 28.86, 30.60, 31.23, 10.60}18.48
ΔD2 = ±0.01 p.u.{1.10, 1.00, 4.05 ± 0.01, 0.00, 2.84}±$28.86
ΔD3 = ±0.01 p.u.{1.10, 1.00, 4.05 ± 0.01, 0.00, 2.84}±$30.60
ΔD4 = ±0.01 p.u.{1.10, 1.00, 4.05 ± 0.01, 0.00, 2.84}±$31.23
Table 5. Operation cost, DRX cost, emission for different DR amount.
Table 5. Operation cost, DRX cost, emission for different DR amount.
DR (p.u.)Base Demand (p.u.)DR (%)Operation Cost with DRX (k$)Operation Cost w/o DRX (k$)DR Transaction Cost (k$)Net Emissions (tonCO2)
0202.260747.49747.49019077
4.02198.181.95724.16716.417.7518806
12.21194.085.92686.31656.9629.3418167
17.03189.258.26672.61623.8448.7617764
21.86184.4210.60679.22591.7874.3217335
26.69179.5912.94687.10560.78107.2816881
31.52174.7615.28734.44529.78178.4816427
36.35169.9417.72802.55498.79266.9215974
39.96166.3219.37909.64475.36382.3415634
43.57162.7121.211040.62451.93516.0815294
Table 6. Average LMPs ($/MWh) at different buses for incremental DR.
Table 6. Average LMPs ($/MWh) at different buses for incremental DR.
The DR AmountAverage LMP ($/MWh)
% DRBus#1Bus#2Bus#3Bus#4Bus#5
050.5151.8255.2556.4916.28
1.9547.9249.1452.3053.4616.28
5.9246.3047.4550.4651.5516.28
8.2645.7646.8949.8550.9216.28
10.6044.1445.2148.0049.0116.28
12.9444.1445.2148.0049.0116.28
15.2844.1445.2148.0049.0116.28
17.7244.1445.2148.0049.0116.28
19.3744.1445.2148.0049.0116.28
21.2144.1445.2148.0049.0116.28
Table 7. The DR supply share among the aggregators in the DR trading hours in a day.
Table 7. The DR supply share among the aggregators in the DR trading hours in a day.
The DR amountDR (p.u.) Provided by the AggregatorsThe Payoff (k$) Achieved by the Aggregators
% DRJ1J2J3J1J2J3
1.952.692.662.6615.1513.7514.41
5.924.273.973.9723.0519.3520.59
8.266.705.175.1726.0917.3319.27
10.608.996.436.4422.9315.8916.19
12.9411.567.577.5726.3715.9816.09
15.2813.988.778.772.92−1.65−1.47
Table 8. DR compensation provided to end-users with incremental DR.
Table 8. DR compensation provided to end-users with incremental DR.
The DR AmountAggregator, A1Aggregator, A2Aggregator, A3
% DRU1 (k$)U2(k$)U3(k$)U1(k$)U2(k$)U3(k$)U1(k$)U2(k$)U3 (k$)
1.952.31002.88002.55300
5.929.060010.76009.51900
8.2617.050016.820014.8800
10.6030.144.28025.371.23023.221.220
12.9436.8414.75030.826.59030.736.570
15.2852.6632.97042.3717.22042.2517.170
17.6271.3355.01053.9332.29053.7836.460.943
19.3795.6276.20069.3352.65069.1560.3410.98
21.21121.2896.33086.8876.941033184.5482.3730.01
Table 9. A comparison of aggregator’s payoff for GenCos’ different strategies.
Table 9. A comparison of aggregator’s payoff for GenCos’ different strategies.
All GenCos Competitive (Case#2)The G3 Strategic (Scenario#1)The G4 Strategic (Scenario#2)The G3, G4 Both Strategic (Scenario#3)
Aggregator’s Payoff ($) without DR0000
Aggregator’s Payoff (k$) with 5.92% DR{23.05, 19.35, 20.59}{23.42,19.70, 20.95,{21.69, 18.06, 19.31}{24.72, 20.94, 22.18}
Aggregator’s Payoff (k$) with 8.26% DR{26.09, 17.32, 19.27}{23.42, 19.70, 20.94}{30.16, 20.72, 22.67}{32.22, 22.43, 24.37}
Total Payoff with 5.92% DR62.9964.0759.0767.84
Total Payoff with 8.26% DR62.6973.8673.5679.01
Relative payoff variation with 5.92% DR (64.07 – 62.69) = 2.20(59.07 − 62.99) = −6.22(67.84 − 62.99) = 7.69
Relative payoff variation with 8.26% DR (73.86 − 62.99) = 17.33(73.56 − 62.69) = 17.33(79.01 − 62.69) = 26.03
Table 10. A comparison of operation costs for a different degree of the strategy adopted by GenCos.
Table 10. A comparison of operation costs for a different degree of the strategy adopted by GenCos.
All GenCos Competitive (Case#2)The G3 Strategic (Scenario#1)The G4 Strategic (Scenario#2)The G3, G4 Strategic (Scenario#3)
Operation cost (k$) without DR747.49772.07753.61764.99
Operation cost (k$) with 5.92% DR686.31706.04687.97711.65
Operation cost (k$) with 8.26% DR672.60690.65681.31697.08
Relative increase (%) in cost at 5.92% DR 2.870.243.69
Relative increase (%) in cost at 8.26% DR 2.681.293.64
Average LMP ($/MWh) in Bus#4 without DR54.7254.9856.4956.49
Average LMP ($/MWh) in Bus#4 at 5.92% DR51.5551.9453.3253.32
Average LMP ($/MWh) in Bus#4 at 8.26% DR50.9251.9452.5752.57

Share and Cite

MDPI and ACS Style

Mohammad, N.; Mishra, Y. The Role of Demand Response Aggregators and the Effect of GenCos Strategic Bidding on the Flexibility of Demand. Energies 2018, 11, 3296. https://doi.org/10.3390/en11123296

AMA Style

Mohammad N, Mishra Y. The Role of Demand Response Aggregators and the Effect of GenCos Strategic Bidding on the Flexibility of Demand. Energies. 2018; 11(12):3296. https://doi.org/10.3390/en11123296

Chicago/Turabian Style

Mohammad, Nur, and Yateendra Mishra. 2018. "The Role of Demand Response Aggregators and the Effect of GenCos Strategic Bidding on the Flexibility of Demand" Energies 11, no. 12: 3296. https://doi.org/10.3390/en11123296

APA Style

Mohammad, N., & Mishra, Y. (2018). The Role of Demand Response Aggregators and the Effect of GenCos Strategic Bidding on the Flexibility of Demand. Energies, 11(12), 3296. https://doi.org/10.3390/en11123296

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