This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

On Efficient Aggregation of Distributed Energy Resources

Zuguang Gao, Khaled Alshehri, and John R. Birge Z. Gao and J. R. Birge are with Booth School of Business, The University of Chicago, Chicago, IL 60637, USA. {zuguang.gao, john.birge}@chicagobooth.eduK. Alshehri is with the Systems Engineering Department, King Fahd University of Petroleum and Minerals (KFUPM), Dhahran, Saudi Arabia. [email protected]
Abstract

The rapid expansion of distributed energy resources (DERs) is one of the most significant changes to electricity systems around the world. Examples of DERs include solar panels, small natural gas-fueled generators, combined heat and power plants, etc. Due to the small supply capacities of these DERs, it is impractical for them to participate directly in the wholesale electricity market. We study in this paper an efficient aggregation model where a profit-maximizing aggregator procures electricity from DERs, and sells them in the wholesale market. The interaction between the aggregator and the DER owners is modeled as a Stackelberg game: the aggregator adopts two-part pricing by announcing a participation fee and a per-unit price of procurement for each DER owner, and the DER owner responds by choosing her payoff-maximizing energy supplies. We show that our proposed model preserves full market efficiency, i.e., the social welfare achieved by the aggregation model is the same as that when DERs participate directly in the wholesale market. We also note that two-part pricing is critical for market efficiency, and illustrate via an example that with one-part pricing, there will be an efficiency loss from DER aggregation, due to the profit-seeking behavior of the aggregator.

I INTRODUCTION

The increasing adoption of distributed energy resources (DERs) challenges fundamental assumptions in existing electricity markets design and operation, as end-consumers not only become proactive, but also small-scale producers who can have significant impacts on the power system [1]. According to the North American Electric Reliability Corporation (NERC), DER is “any resource on the distribution system that produces electricity and is not otherwise included in the formal NERC definition of the Bulk Electric System (BES)”[2]. The existence of DERs causes a fundamental shift, because electric power demand is largely assumed to be fixed from independent market operators’ (ISOs) perspective [3, 4, 5, 6, 7]. This fundamental shift calls for the quantification of the impacts of DERs on wholesale markets design. Such quantification is not easy; ISOs do not have oversight over the distribution power network, and hence, cannot include DER owners as market participants. Furthermore, even if ISOs can oversee the distribution power system, it would be a significant burden and impractical to include DER supply directly into the wholesale electricity market operations through ISOs, due to the communicational, computational, and operational complexity. 111The number of market participants would grow from hundreds to millions.

Refer to caption
Figure 1: Overall interactions in the proposed efficient aggregation model

Different models have been proposed to include DER supply into wholesale electricity markets. One possibility is to have a Distribution System Operator (DSO) acting as a market manager at the distribution level, and finding socially-optimal dispatch, similar to ISOs [8, 9, 10, 11, 12, 13]. Another possibility is to have fully-distributed electricity market designs, in which end-consumers can trade energy among themselves [14, 15]. The third model, which we adopt here, is to have an aggregating company whose role is to collect energy from DER owners, and participate in the wholesale market as a producer of electricity. This model has been adopted by California ISO and New York ISO [16, 17], and seems to be realistic for practical implementation, especially with FERC’s recent Order No. 2222 in September, 2020 [18]. The aggregator here buys DER supply from their owners, and bids directly to the ISO similar to generating companies. The relationship between the aggregator and prosumers in the same geographical footprint is naturally monopolistic, where aggregators become price-making in retail electricity markets as they can send price offers to prosumers in order to collect DER supply. Such price offers need to be high-enough so that DER owners are attracted to sell, but small-enough so that the aggregator can maximize her profits. This profit-seeking behavior can impact the overall electricity market efficiency, but at the same time, due to the impracticality of direct DER participation into wholesale markets, aggregators are necessary and important players. This gives rise to the following important question: In the presence of a profit-seeking and a monopolistic aggregator, is there an aggregation model that can attain a socially optimal (efficient) market outcome? The presence of such a mechanism can in fact be significant. First, in reality, aggregators are mostly profit-seeking, and are often monopolistic, which makes the markets prone to efficiency losses as demonstrated in [19, 20]. Second, it is infeasible for DER owners to participate in wholesale markets, so the presence of such intermediaries is inevitable. Third, if such a mechanism can be designed, it would address various debates surrounding whether or not DER aggregation need to be done by profit-making entities, or social intermediaries.

In this paper, we address the above question by proposing a DER aggregation model that yields efficient market outcomes, even when a DER aggregator 𝒜{\cal A} is monopolistic and profit-seeking. We analytically show that with our proposed model, the optimal DER capacities being integrated through 𝒜{\cal A} are equivalent to a benchmark (ideal, not realistic) case in which DER owners can directly participate in the wholesale market. Briefly, our aggregation model utilizes a two-price offers (𝐏,𝐩)\bf(P,p) from 𝒜{\cal A} to the prosumers, where 𝐏\bf P corresponds to a fixed DER owner participation costs (connection charges), and the other prices 𝐩\bf p are for marginal acquisition of DER capacities. We remark here that having one-price offer 𝐩\bf p would not yield efficient market outcomes due to the monopolistic behavior of 𝒜{\cal A}, as shown in [19, 20]. For an overall illustration, refer to Figure 1. Here, we adopt a game-theoretic analysis of the relationship between the prosumers and 𝒜{\cal A}, and show that the equilibrium prices (𝐏,𝐩)(\bf P^{*},p^{*}) and DER capacities 𝐱(𝐏,𝐩)\bf{x}^{*}(\bf P^{*},p^{*}) attain efficient market outcomes.

The paper is organized as follows: Section II discusses the benchmark ideal model. In Section III, we propose an efficient DER integration model, and state the main theoretical results of the paper. Via illustrative examples, we conduct simulations in Section IV to gain some practical insights. The paper concludes in Section V.

Notation

We let \mathbb{R} denote the set of real numbers, and +\mathbb{R}_{+} (resp. ++\mathbb{R}_{++}) denote the set of nonnegative (resp. positive) numbers. For zz\in\mathbb{R}, we let [z]+:=max{z,0}[z]^{+}:=\max\{z,0\}. In any optimization problem, a decision variable xx at optimality is denoted by xx^{*}. We use boldface to distinguish a vector 𝐱\bf{x} from a scalar xx. For an event 𝒳{\cal X}, we let 𝟙{𝒳}\mathds{1}\{{\cal X}\} denote the indicator function.

II DIRECT PROSUMER PARTICIPATION MODEL (BENCHMARK)

In this section, we introduce an ideal model where prosumers can participate (buy or sell) directly in the wholesale market, and no aggregator is present. This model, though not realistic, serves as a benchmark for evaluating the market efficiency of our following models. There are three parties in this model: prosumers who can buy and sell energy directly in the wholesale market, conventional generators who generate and sell electricity in the wholesale market, and an independent system operator (ISO) who clears the wholesale market. In the following, we describe the optimization problems solved by each of these three parties. Figure 2 illustrates the interactions between prosumers and the ISO.

Refer to caption
Figure 2: Benchmark model

II-A Prosumer’s Problem

Consider a power network with nn nodes (locations). At each location ii, there are NiN_{i} number of prosumers. We assume that each prosumer jj at location ii can purchase energy at the wholesale price λi\lambda^{i}. Furthermore, prosumer jj at location ii is endowed with a capacity Cji0{C}^{i}_{j}\geq 0 of power production222CjiC^{i}_{j} can be zero for consumers with no DER supply, so the model includes those regular consumers. from a collection of resources, such as solar panels, wall-mounted batteries, and plug-in electric vehicles. The power produced can be consumed locally by prosumer jj or sold (partially) back to the wholesale market, again at the wholesale price λi\lambda^{i}. The prosumer observes λi\lambda^{i} and decides the amount of energy to buy/sell. Let ujiu^{i}_{j} be prosumer jj’s utility of power consumption. We make the following assumption on prosumers’ utility of consumption.

Assumption 1.

Each prosumer’s utility of consumption ujiu^{i}_{j} is increasing, strictly concave, and differentiable. Furthermore, the domain of ujiu^{i}_{j} is [0,Z][0,Z] where Z>CjiZ>C^{i}_{j} is some (large) upper bound of the amount of energy a prosumer can consume. We assume that uji(z)z\frac{\partial u^{i}_{j}(z)}{\partial z}\to\infty as z0z\to 0, and uji(z)z0\frac{\partial u^{i}_{j}(z)}{\partial z}\to 0 as zZz\to Z.

We can then write prosumer jj’s optimization problem as

maxxjidjiπji(xji,dji):=λi(xjidji)+uji(Cji+djixji)\displaystyle\max_{x^{i}_{j}-d^{i}_{j}}\pi^{i}_{j}(x^{i}_{j},d^{i}_{j}):=\lambda^{i}(x^{i}_{j}-d^{i}_{j})+u^{i}_{j}(C^{i}_{j}+d^{i}_{j}-x^{i}_{j}) (1)
s.t. CjiZxjidjiCji,\displaystyle\quad\text{s.t. }\qquad C^{i}_{j}-Z\leq x^{i}_{j}-d^{i}_{j}\leq C^{i}_{j},

where djid^{i}_{j} is the amount of energy prosumer jj purchases and xjix^{i}_{j} is the amount of energy she sells. For selling xjix^{i}_{j} at the wholesale price, prosumer jj receives λixji\lambda^{i}x^{i}_{j}; for buying djid^{i}_{j} at the wholesale price, she is charged λidji\lambda^{i}d^{i}_{j}. The prosumer leaves with Cji+djixjiC^{i}_{j}+d^{i}_{j}-x^{i}_{j} to consume and her utility from consumption would be uji(Cji+djixji)u^{i}_{j}(C^{i}_{j}+d^{i}_{j}-x^{i}_{j}). While Assumption 1 further imposes strict concavity of prosumers’ utilities, our analysis throughout this paper remains largely applicable to generic concave utilities, but strict concavity allows us to derive unique analytical solutions and gain deep insights.

We note that since prosumers at location ii can buy and sell energy at the same price λi\lambda^{i}, prosumer jj essentially has xjidjix^{i}_{j}-d^{i}_{j} as the single decision to make in solving (1), while we keep both xjix^{i}_{j} and djid^{i}_{j} in the notation to clearly represent the amount of selling and buying, and to be consistent with the notation in the following sections.

Lemma 1.

Under Assumption 1, given any λi\lambda^{i}, there exists a unique optimal solution (xjidji)(x^{i}_{j}-d^{i}_{j})^{*} for the prosumer’s problem (1) which satisfies

uji(z)z|z=Cji(xjidji)=λi.\displaystyle\frac{\partial u^{i}_{j}(z)}{\partial z}\bigg{|}_{z=C^{i}_{j}-(x^{i}_{j}-d^{i}_{j})^{*}}=\lambda^{i}. (2)

The above Lemma directly follows from the properties of problem (1). Without loss of generality, and for ease of exposition, we restrict our attention to the case in which at most one of xjix^{i}_{j} and djid^{i}_{j} can be nonzero. We first solve (1) for (xjidji)(x^{i}_{j}-d^{i}_{j})^{*}. Then, we let xji=[(xjidji)]+{x^{i}_{j}}^{*}=\left[(x^{i}_{j}-d^{i}_{j})^{*}\right]^{+} and dji=[(xjidji)]+{d^{i}_{j}}^{*}=\left[-(x^{i}_{j}-d^{i}_{j})^{*}\right]^{+}. We also use the notation xji(λi){x^{i}_{j}}^{*}(\lambda^{i}) and dji(λi){d^{i}_{j}}^{*}(\lambda^{i}) to denote the optimal response of prosumer jj at location ii for a given wholesale market price λi\lambda^{i}.

II-B Generator’s Problem

Let GiG_{i} denote the number of dispatchable conventional generators at location ii. Generator jj at location ii chooses to supply yji[y¯ji,y¯ji]y^{i}_{j}\in\left[\underline{y}^{i}_{j},\overline{y}^{i}_{j}\right]. Let cji(yji)c^{i}_{j}(y^{i}_{j}) be denote the production cost. We have the following assumption.

Assumption 2.

Each generator’s cost function cjic^{i}_{j} is increasing, strictly convex, and differentiable in [y¯ji,y¯ji]\left[\underline{y}^{i}_{j},\overline{y}^{i}_{j}\right]. Furthermore, we let cji(yji)yji0\frac{\partial c^{i}_{j}(y^{i}_{j})}{\partial y^{i}_{j}}\to 0 as yjiy¯jiy^{i}_{j}\to\underline{y}^{i}_{j} and cji(yji)yji\frac{\partial c^{i}_{j}(y^{i}_{j})}{\partial y^{i}_{j}}\to\infty as yjiy¯jiy^{i}_{j}\to\overline{y}^{i}_{j}.

By selling yjiy^{i}_{j}, generator jj earns a compensation λiyi\lambda^{i}y^{i}. Given a wholesale price λi\lambda^{i}, generator jj maximizes its payoff by solving

maxyji[y¯ji,y¯ji]π^ji(yji):=λjiyjicji(yji).\displaystyle\max_{y^{i}_{j}\in\left[\underline{y}^{i}_{j},\overline{y}^{i}_{j}\right]}\hat{\pi}^{i}_{j}(y^{i}_{j}):=\lambda^{i}_{j}y^{i}_{j}-c^{i}_{j}(y^{i}_{j}). (3)
Lemma 2.

Under Assumption 2, given any λi\lambda^{i}, there exists a unique optimal solution yji{y^{i}_{j}}^{*} for the generator’s problem (3) which satisfies

cji(yji)yji|yji=yji=λi.\displaystyle\frac{\partial c^{i}_{j}(y^{i}_{j})}{\partial y^{i}_{j}}\bigg{|}_{y^{i}_{j}={y^{i}_{j}}^{*}}=\lambda^{i}. (4)

The above Lemma directly follows from the properties of problem (3). We also use the notation yji(λi){y^{i}_{j}}^{*}(\lambda^{i}) to denote the optimal response of generator jj at location ii at a given wholesale price λi\lambda^{i}.

II-C The Economic Dispatch Problem

Wholesale electricity markets in the United States and other countries are managed by independent system operators (ISOs). An ISO clears the market by matching supply and demand via social welfare maximization (in practice, this is often done by production cost minimization to meet fixed system demands [4, 5, 6]), while ensuring that the power flows satisfy the network and line capacity constraints. Specifically, let Xi=j[Ni]xjiX^{i}=\sum_{j\in[N_{i}]}x^{i}_{j} be the total power supply at node ii from prosumers; let Yi=j[Gi]yjiY^{i}=\sum_{j\in[G_{i}]}y^{i}_{j} be the total power supply at node ii from conventional generators; and let Di=j[Ni]djiD^{i}=\sum_{j\in[N_{i}]}d^{i}_{j} be the total demand at node ii. Furthermore, we let 𝐟{\bf f} be the vector of capacities of transmission lines in the power network, and let 𝐁{\bf B} be the line-to-node incidence matrix. The system operator chooses a vector 𝐡{\bf h}, where each element hih^{i} is the net injection to node ii. We have the following network constraints:

𝐡=𝐃𝐘𝐗,\displaystyle{\bf h=D-Y-X}, (5a)
𝟏T𝐡=0,𝐁𝐡𝐟.\displaystyle{\bf 1}^{T}{\bf h}=0,\qquad{\bf Bh}\leq{\bf f}. (5b)

We note that (5a) ensures the total supply matches the total demand at each node; (5b) ensures that the total net injection by the system operator is zero over the power network (here, 𝟏{\bf 1} is a vector of ones), and the total power transmission at each line does not exceed its capacity. In addition to the network constraints, the ISO needs to also consider all participant-specific constraints described earlier in problems (1) and (3):

𝐂𝐙𝐱𝐝𝐂,𝐲¯𝐲𝐲¯{\bf C-Z\leq x-d\leq C,\qquad\underline{y}\leq y\leq\overline{y}} (6)

The objective of the system operator is to maximize the social welfare, which includes the prosumer surplus (PS), generator surplus (GS), and merchandizing surplus (MS):

𝖯𝖲\displaystyle{\sf PS} :=i[n](j[Ni]uji(djixji+Cji)λi(hi+Yi)),\displaystyle:=\sum_{i\in[n]}\left(\sum_{j\in[N_{i}]}u^{i}_{j}(d^{i}_{j}-x^{i}_{j}+{C}^{i}_{j})-\lambda^{i}(h^{i}+Y^{i})\right), (7a)
𝖦𝖲\displaystyle{\sf GS} :=i[n]j[Gi](λiyjicji(yji)),𝖬𝖲:=i[n]λihi.\displaystyle:=\sum_{i\in[n]}\sum_{j\in[G_{i}]}\left(\lambda^{i}y^{i}_{j}-c^{i}_{j}(y^{i}_{j})\right),\ {\sf MS}:=\sum_{i\in[n]}\lambda^{i}h^{i}. (7b)

where we have imposed the relationship (5a) in deriving (7a). The social welfare that the system operator optimizes is the sum of PS, GS, and MS. After canceling terms, the social welfare can be written as

𝒲B:=i[n](j[Ni]uji(dji+Cjixji)j[Gi]cji(yji)).\displaystyle\mathcal{W}_{B}:=\sum_{i\in[n]}\left(\sum_{j\in[N_{i}]}u^{i}_{j}(d^{i}_{j}+{C}^{i}_{j}-x^{i}_{j})-\sum_{j\in[G_{i}]}c^{i}_{j}(y^{i}_{j})\right). (8)

The system operator’s economic dispatch problem is then:

max𝒲B(𝐡,𝐱𝐝,𝐲)\displaystyle\max\qquad\mathcal{W}_{B}({\bf h,x-d,y}) (9)
subject to(5)(6).\displaystyle\text{subject to}\qquad\eqref{eq:directcons}-\eqref{eq:directcons_p}.
Assumption 3.

The system operator’s economic dispatch problem (9) is feasible.

Proposition 1 (Competitive Equilibrium).

Under Assumptions 1-3, there exists a unique optimal solution (𝐡,(𝐱𝐝),𝐲)\left({\bf h^{*},(x-d)^{*},y^{*}}\right) to (9). Denote the optimal Lagrange multipliers of constraints (5a) by 𝛌{\boldsymbol{\lambda}}. Then, the following statements are true:

  • (𝐱𝐝)\left(\bf x-d\right)^{*} and 𝝀\boldsymbol{\lambda} satisfy (2).

  • 𝐲{\bf y}^{*} and 𝝀\boldsymbol{\lambda} satisfy (4).

The above proposition states that solving the system operator’s problem (9) leads to a competitive equilibrium. From the prospective of prosumer jj at node ii, this means that given the wholesale market price λi\lambda^{i} (the optimal Lagrange multiplier of  (5a) for the same node), the corresponding solution to her problem, which satisfies (2), is the same as the optimal decision made by the ISO via solving  (9). This is also true for all other prosumers and generators. Having all market participants being satisfied with the the competitive equilibrium, it serves as a good benchamrk for market efficiency.

III PROPOSED EFFICIENT AGGREGATION MODEL

The direct participation model introduced in the previous section is a benchmark: prosumers are allowed to participate directly and sell their production in the wholesale market. In reality, the supply capacities of prosumers are typically too small for consideration in the wholesale market. 333California ISO requires a minimum of 0.5𝖬𝖶0.5{\sf MW} for a DER aggregator to participate [21], and New York ISO requires 0.1𝖬𝖶0.1{\sf MW}. Also, computing the dispatch and settlement for a large number of prosumers raises an untenable computational burden on the system operator. The presence of DER aggregators brings benefits to the system, as they open the door for DER owners to participate and bring more flexibility to the grid. However, the profit-seeking nature of these aggregators can cause efficiency losses [19, 20]. To resolve this, we propose in this section an efficient aggregation model under two-part pricing. In this model, prosumers sell part of their DER supply productions to an aggregator 𝒜\mathcal{A}, based on the price offers made by 𝒜\mathcal{A}. The interactions between 𝒜\mathcal{A} and prosumers are modeled as a Stackelberg game [22]. The aggregator acts as a leader and announces a price pair (Pji,pji)(P^{i}_{j},p^{i}_{j}) for each prosumer jj at location ii. The prosumer follows by choosing the amount of energy to sell. If the prosumer decides to sell a nonzero fraction of her capacity to the aggregator, she pays to the aggregator a participation fee PjiP^{i}_{j}, and earns the price pjip^{i}_{j} for each unit of energy sold. The aggregator 𝒜\mathcal{A} then sells all these procured capacity to the wholesale market at the wholesale price λi\lambda^{i}. The goal of the aggregator 𝒜\mathcal{A} is to choose prices (Pji,pji)(P^{i}_{j},p^{i}_{j}) that maximize her profit, while anticipating how DER owners would respond. Note that here, DER owners have access only to one aggregator, so 𝒜{\cal A} is in fact monopolistic, which further signifies the importance of our mechanism as it yields socially-optimal outcomes.

We also note that differential pricing is allowed in this model, i.e., the price pair (Pji,pji)(P^{i}_{j},p^{i}_{j}) can be set differently for different prosumers. While it is reasonable to have varying prices depending on locations [23], there exist some debates on whether prices should be allowed to set differently for prosumers from the same location. The legal issue is not the main focus of this paper. Though we assumed differential pricing in the model, the equilibrium we obtain in the end has the same marginal price at each location, that is, pji=λi,j[Ni]p^{i}_{j}=\lambda^{i},\ \forall j\in[N_{i}]. The participation fee can be differentiated by, for example, mailing different coupons to different prosumers to encourage their participations, which is arguably more justifiable. In the remainder of this section, we show that this aggregation model achieves the same socially-optimal market outcomes as in the direct participation model.

III-A Prosumer’s Problem

Consider prosumer jj at location ii. Upon seeing the prices (Pji,pji)(P^{i}_{j},p^{i}_{j}) announced by 𝒜\mathcal{A}, prosumer jj decides if she would sell part of her capacity to 𝒜\mathcal{A}. If she chooses so, she would pay a fee PjiP^{i}_{j} to 𝒜\mathcal{A}, and receives pjip^{i}_{j} for each unit of energy sold. We may write prosumer jj’s payoff as

πji(xji,dji):=\displaystyle\pi_{j}^{i}(x^{i}_{j},d^{i}_{j}):=
{pjixjiPji+uji(dji+Cjixji)λidji, if xji>0,uji(dji+Cji)λidji, if xji=0,\displaystyle\qquad\begin{cases}p^{i}_{j}x^{i}_{j}-P^{i}_{j}+u^{i}_{j}\left(d^{i}_{j}+{C}^{i}_{j}-x^{i}_{j}\right)-\lambda^{i}d^{i}_{j},\text{ if }x^{i}_{j}>0,\\ u^{i}_{j}\left(d^{i}_{j}+{C}^{i}_{j}\right)-\lambda^{i}d^{i}_{j},\hfill\text{ if }x^{i}_{j}=0,\end{cases} (10)

Given (Pji,pji,λi)(P^{i}_{j},p^{i}_{j},\lambda^{i}), prosumer jj solves:

maxxji[0,Cji],dji[0,ZCji+xji]πji(xji,dji),\displaystyle\max_{x^{i}_{j}\in[0,C^{i}_{j}],d^{i}_{j}\in[0,Z-C^{i}_{j}+x^{i}_{j}]}\ \pi_{j}^{i}(x^{i}_{j},d^{i}_{j}), (11)

where djid^{i}_{j} is the amount of energy prosumer jj purchases at wholesale market price λi\lambda^{i}, 444In real markets, end-use customers are represented in the wholesale market via utility companies. Here, for ease of exposition, we omit the analysis of the relationship between utility companies and end-use customers and assume that they represent the true preferences of their customers in the wholesale market without modifications. This allows us to focus on the strategic aggregation component of the problem. Nevertheless, our insights remain largely applicable if this assumption is relaxed. and xjix^{i}_{j} is the amount of energy she sells to the aggregator. For buying djid^{i}_{j} at the wholesale price, prosumer jj is charged λidji\lambda^{i}d^{i}_{j}. If the prosumer does not sell (xji=0x^{i}_{j}=0), she has a total of dji+Cjid^{i}_{j}+C^{i}_{j} to consume, and her utility from consumption would be uji(dji+Cji)u^{i}_{j}(d^{i}_{j}+C^{i}_{j}). If the prosumer chooses to sell xji>0x^{i}_{j}>0 to 𝒜\mathcal{A}, she is charged a participation fee PjiP^{i}_{j}, and receives a compensation pjixjip^{i}_{j}x^{i}_{j}. The prosumer would have dji+Cjixjid^{i}_{j}+C^{i}_{j}-x^{i}_{j} to consume in this case and her utility from consumption would be uji(dji+Cjixji)u^{i}_{j}(d^{i}_{j}+C^{i}_{j}-x^{i}_{j}).

Let xji(Pji,pji,λi){x^{i}_{j}}^{*}(P^{i}_{j},p^{i}_{j},\lambda^{i}) and dji(Pji,pji,λi){d^{i}_{j}}^{*}(P^{i}_{j},p^{i}_{j},\lambda^{i}) denote the optimal response of prosumer jj given aggregator’s announced prices (Pji,pji)(P^{i}_{j},p^{i}_{j}) and the wholesale market price λi\lambda^{i}.555Sometimes, we drop arguments for simplicity. Note that if pji>λip^{i}_{j}>\lambda^{i}, the prosumer can arbitrage by buying at the price λi\lambda^{i} and selling at a higher price pjip^{i}_{j}. This will result in the prosumer earning infinite payoff and aggregator losing infinite profit, which would be avoided by the aggregator. Therefore, we may without loss of generality restrict our discussions to the case when pjiλi,i[n],j[Ni]p^{i}_{j}\leq\lambda^{i},\forall i\in[n],j\in[N_{i}]. In the case pji=λip^{i}_{j}=\lambda^{i}, similar to the direct participation model, we may enforce that xjix^{i}_{j} and djid^{i}_{j} cannot both be nonzero. We then have the following lemma on the optimal response of prosumers.

Lemma 3.

Consider an arbitrary prosumer jj at location ii. Let (z1,z2)(z_{1},z_{2}) be such that

uji(z)z|z=z1=λi,uji(z)z|z=z2=pji.\displaystyle\frac{\partial u^{i}_{j}(z)}{\partial z}\Big{|}_{z=z_{1}}=\lambda^{i},\qquad\frac{\partial u^{i}_{j}(z)}{\partial z}\Big{|}_{z=z_{2}}=p^{i}_{j}. (12)

Then, under Assumption 1, both (z1,z2)(z_{1},z_{2}) exist and are unique. Furthermore, prosumer jj’s optimal response can be described as follows.

  • If Cjiz2C^{i}_{j}\leq z_{2}, then, we have

    dji=[z1Cji]+,xji=0.{d^{i}_{j}}^{*}=\left[z_{1}-C^{i}_{j}\right]^{+},\qquad{x^{i}_{j}}^{*}=0.
  • If Cji>z2C^{i}_{j}>z_{2}, then, we have

    xji=(Cjiz2)𝟙{𝒳},dji=0,{x^{i}_{j}}^{*}=\left(C^{i}_{j}-z_{2}\right)\cdot\mathds{1}\left\{{\cal X}\right\},\qquad{d^{i}_{j}}^{*}=0,

    where

    𝒳:={Pjipji(Cjiz2)+uji(z2)uji(Cji)}.{\cal X}:=\{P^{i}_{j}\leq p^{i}_{j}\left(C^{i}_{j}-z_{2}\right)+u^{i}_{j}\left(z_{2}\right)-u^{i}_{j}\left(C^{i}_{j}\right)\}.

The above Lemma states that if the capacity exceeds a certain value z2z_{2} (note that z2z_{2} is the value at which the marginal utility of consumption is equal to the aggregator’s marginal price offer pjip^{i}_{j}) and the upfront fee PjiP^{i}_{j} is not too high, then prosumers have an incentive to sell DER supply to the aggregator. If the capacity is too small or the upfront fee is too high, then prosumers would not sell DER supply and prefer to consume it locally.

III-B Aggregator’s Problem

The DER aggregator 𝒜\mathcal{A} collects power from prosumers and sells it in the wholesale market. By offering the prices (Pji,pji)(P^{i}_{j},p^{i}_{j}) to each prosumer jj at location ii, 𝒜\mathcal{A} procures a total capacity of j[Ni]xji(Pji,pji)\sum_{j\in[N_{i}]}{x^{i}_{j}}^{*}(P^{i}_{j},p^{i}_{j}) from location ii. 𝒜\mathcal{A} then sells it at the wholesale market price λi\lambda^{i}. Given the wholesale market price λi\lambda^{i}, the aggregator’s profit from prosumer jj at location ii is

Πji(Pji,pji):=\displaystyle\Pi^{i}_{j}\left(P^{i}_{j},p^{i}_{j}\right):= (13)
Pji𝟙{xji(Pji,pji)>0}+(λipji)xji(Pji,pji).\displaystyle\qquad P^{i}_{j}\mathds{1}\left\{{x^{i}_{j}}^{*}(P^{i}_{j},p^{i}_{j})>0\right\}+(\lambda^{i}-p^{i}_{j}){x^{i}_{j}}^{*}(P^{i}_{j},p^{i}_{j}).

Anticipating the response functions xji(Pji,pji){x^{i}_{j}}^{*}(P^{i}_{j},p^{i}_{j})’s, she seeks to maximize her overall profit:

max𝐏𝟎,𝐩𝟎i[n]j[Ni]Πji(Pji,pji)=:Π^(𝐏,𝐩),\displaystyle\max_{\bf P\geq 0,p\geq 0}\underbrace{\sum_{i\in[n]}\sum_{j\in[N_{i}]}\Pi^{i}_{j}\left(P^{i}_{j},p^{i}_{j}\right)}_{=:\hat{\Pi}(\bf P,p)}, (14)

Aggregator 𝒜\mathcal{A}’s profit is composed of two parts; the total participation fees charged to prosumers who sell positive amount of energy and the profits earned by reselling the procured energy in the wholesale market. We note that 𝒜{\cal A}’s profit from prosumer jj depends only on the response and the prices of prosumer jj, and there is no coupling among prosumers’ problems. It should be clear that the Aggregator’s problem (14) can be decomposed to optimizing (Pji,pji)(P^{i}_{j},p^{i}_{j}) for each prosumer at each location:

maxPji0,pji0Πji(Pji,pji).\displaystyle\max_{P^{i}_{j}\geq 0,p^{i}_{j}\geq 0}\Pi^{i}_{j}\left(P^{i}_{j},p^{i}_{j}\right). (15)

Solving (15) for each i[n]i\in[n] and j[Ni]j\in[N_{i}] would lead to the vectors (𝐏,𝐩)\bf(P^{*},p^{*}), which constitute an optimal solution for (14). We have the following result on the optimal decisions of the aggregator.

Lemma 4.

Under Assumption 1, consider an arbitrary prosumer jj at location ii, with z1z_{1} as in (12), and wholesale market price  λi\lambda^{i}, When z1<Cjiz_{1}<C^{i}_{j}, the aggregator’s optimal pricing decision is

pji\displaystyle{p^{i}_{j}}^{*} =λi,\displaystyle=\lambda^{i}, (16a)
Pji\displaystyle{P^{i}_{j}}^{*} =λi(Cjiz1)+uji(z1)uji(Cji).\displaystyle=\lambda^{i}(C^{i}_{j}-z_{1})+u^{i}_{j}(z_{1})-u^{i}_{j}(C^{i}_{j}). (16b)

When z1Cjiz_{1}\geq C^{i}_{j}, prosumer jj will not sell DER supply to 𝒜{\cal A}, i.e., Πji(Pji,pji)=0\Pi^{i}_{j}\left(P^{i}_{j},p^{i}_{j}\right)=0, for any (Pji,pji)+2({P^{i}_{j}},{p^{i}_{j}})\in\mathbb{R}_{+}^{2}.

Lemma 4 provides an optimal solution (Pji,pji)({P^{i}_{j}}^{*},{p^{i}_{j}}^{*}) to (15), and the collection  (𝐏,𝐩)\bf(P^{*},p^{*}) form an optimal solution to (14). We also note that the optimal solution may not be unique in general: 𝒜{\cal A} may deviate from (16) by further decreasing pji{p^{i}_{j}}^{*} and increasing Pji{P^{i}_{j}}^{*} to earn the same profit, while keeping the response of the prosumer jj being unchanged. This lemma states that there is an optimal pricing scheme which sets the marginal price pjip^{i}_{j} to be the wholesale market price λi\lambda^{i}, and all prosumers at the same location will be offered the same pji=λip^{i}_{j}=\lambda^{i}. The participation fees PjiP^{i}_{j}, however, will be charged differently for prosumers with different utility functions. A keen reader would observe that in view of Lemmas 1 and 3, the optimal DER supply to the wholesale market is the same with and without the aggregator 𝒜{\cal A}, thus making this aggregation model economically efficient with socially optimal market outcomes.

III-C Aggregator-Prosumers Interaction as a Stackelberg Game

Let 𝒢(𝝀)\mathcal{G}(\boldsymbol{\lambda}) denote the Stackelberg game among the aggregator and the prosumers for a given vector of wholesale market prices 𝝀\boldsymbol{\lambda}. Aggregator 𝒜\mathcal{A} acts as a Stackelberg leader and sets the prices (𝐏,𝐩)\bf(P,p). The prosumer follows by responding with xji(Pji,pji,λi){x^{i}_{j}}^{*}(P^{i}_{j},p^{i}_{j},\lambda^{i}). We now define the equilibrium of the game.

Definition 1.

(𝐏,𝐩,𝐱(𝐏,𝐩))\left(\bf P^{*},p^{*},{x}^{*}(P^{*},p^{*})\right) constitutes a Stackelberg equilibrium of the game 𝒢(𝛌)\mathcal{G}(\boldsymbol{\lambda}) if

  • Prosumers: For any xji[0,Cji]x^{i}_{j}\in\left[0,C^{i}_{j}\right], we have

    πji(xji(Pji,pji),Pji,pji,λi)πji(xji,Pji,pji,λi)\pi_{j}^{i}\left({x^{i}_{j}}^{*}({P^{i}_{j}},{p^{i}_{j}}),{P^{i}_{j}},{p^{i}_{j}},\lambda^{i}\right)\geq\pi_{j}^{i}\left({x^{i}_{j}},{P^{i}_{j}}^{*},{p^{i}_{j}}^{*},\lambda^{i}\right)

    for all i[n]i\in[n] and j[Ni]j\in[N_{i}].

  • Aggregator: For all 𝐏𝟎,𝐩𝟎{\bf P\geq 0,p\geq 0}, we have

    Π^(𝐏,𝐩,𝐱(𝐏,𝐩),𝝀)Π^(𝐏,𝐩,𝐱(𝐏,𝐩),𝝀),\hat{\Pi}({\bf P^{*},p^{*},x^{*}(P^{*},p^{*})},\boldsymbol{\lambda})\geq\hat{\Pi}(\bf{P,p,x^{*}(P,p)},\boldsymbol{\lambda}),

We then have the following Stackelberg equilibrium for the game 𝒢(𝝀){\cal G}(\boldsymbol{\lambda}), which follows directly from the prosumer’s optimal response (Lemma 3) and the aggregator’s optimal pricing (Lemma 4).

Proposition 2.

Under Assumption 1, the game 𝒢(𝛌)\mathcal{G}(\boldsymbol{\lambda}) admits a Stackelberg equilibrium that satisfies

Pji=[λi(Cjizji)+uji(zji)uji(Cji)]+\displaystyle{P^{i}_{j}}^{*}=\left[\lambda^{i}(C^{i}_{j}-z^{i}_{j})+u^{i}_{j}\left(z^{i}_{j}\right)-u^{i}_{j}(C^{i}_{j})\right]^{+} (17a)
pji=λi\displaystyle{p^{i}_{j}}^{*}=\lambda^{i} (17b)
xji(Pji,pji)=[Cjizji]+\displaystyle{x^{i}_{j}}^{*}({P^{i}_{j}}^{*},{p^{i}_{j}}^{*})=\left[C^{i}_{j}-z^{i}_{j}\right]^{+} (17c)

for each prosumer jj at each location ii, where zjiz^{i}_{j} satisfies uji(z)z|z=zji=λi\frac{\partial u^{i}_{j}(z)}{\partial z}\Big{|}_{z=z^{i}_{j}}=\lambda^{i}, and dji=[zjiCji]+{d^{i}_{j}}^{*}=\left[z^{i}_{j}-C^{i}_{j}\right]^{+}.

We note that the game 𝒢(𝝀){\cal G}(\boldsymbol{\lambda}) may admit other Stackelberg equilibria, but Proposition 2 provides the most economically efficient one. In case of non-uniqueness, this economically efficient equilibrium corresponds to the case where prosumers slightly prefer participation, i.e., if selling xjix^{i}_{j} amount of energy earns the prosumer jj the same πji\pi^{i}_{j} as not selling, then she chooses to sell this xjix^{i}_{j}. 666Alternatively, one can impose a slight perturbation of Pji{P^{i}_{j}}^{*} to Pjiϵ{P^{i}_{j}}^{*}-\epsilon, where ϵ>0\epsilon>0, to ensure maximum DER supply. The above equilibrium is quite intuitive, 𝒜{\cal A} passes the location marginal price λi\lambda^{i} obtained from wholesale market outcomes as is to prosumers, so 𝒜{\cal A} has no marginal profits from 𝐩{\bf p}^{*}. Instead, 𝒜{\cal A} makes all of the profits from the upfront participation fees 𝐏{\bf P^{*}}.

III-D Generator’s Problem

For a given wholesale market price λi\lambda^{i}, the conventional generators solve the same problem as described in Section II-B, and the result of Lemma 2 still applies under the current aggregation model.

III-E The Economic Dispatch Problem

The system operator solves an optimization problem similar to that in the direct participation model as described in Section II-C. The network constraints (5) remain valid. Besides, the ISO also considers participant-specific constraints:

𝟎𝐱𝐂,𝟎𝐝𝐙𝐂+𝐱,𝐲¯𝐲𝐲¯{\bf 0\leq x\leq C,\quad 0\leq d\leq Z-C+x,\quad\underline{y}\leq y\leq\overline{y}} (18)

The objective of the system operator is to maximize the social welfare, which now includes the prosumer surplus (PS), aggregator surplus (AS), generator surplus (GS), and merchandizing surplus (MS):

𝖯𝖲\displaystyle{\sf PS} :=i[n]j[Ni](uji(djixji+Cji)λidji\displaystyle:=\sum_{i\in[n]}\sum_{j\in[N_{i}]}\Bigg{(}u^{i}_{j}(d^{i}_{j}-x^{i}_{j}+{C}^{i}_{j})-\lambda^{i}d^{i}_{j}
+pjixjiPji𝟙{xji>0})\displaystyle\qquad\qquad\qquad\qquad+p^{i}_{j}x^{i}_{j}-P^{i}_{j}\mathds{1}\left\{x^{i}_{j}>0\right\}\Bigg{)} (19a)
𝖠𝖲\displaystyle{\sf AS} :=i[n]j[Ni](Pji𝟙{xji>0}+λixjipjixji)\displaystyle:=\sum_{i\in[n]}\sum_{j\in[N_{i}]}\Bigg{(}P^{i}_{j}\mathds{1}\left\{x^{i}_{j}>0\right\}+\lambda^{i}x^{i}_{j}-p^{i}_{j}x^{i}_{j}\Bigg{)} (19b)
𝖦𝖲\displaystyle{\sf GS} :=i[n]j[Gi](λiyjicji(yji)),𝖬𝖲:=i[n]λihi.\displaystyle:=\sum_{i\in[n]}\sum_{j\in[G_{i}]}\left(\lambda^{i}y^{i}_{j}-c^{i}_{j}(y^{i}_{j})\right),{\sf MS}:=\sum_{i\in[n]}\lambda^{i}h^{i}. (19c)

The social welfare is the sum of the above four terms. By the supply-demand balance 𝐡=𝐃𝐘𝐗{\bf h=D-Y-X}, and after canceling terms, we write the social welfare as

𝒲A:=i[n](j[Ni]uji(dji+Cjixji)j[Gi]cji(yji)),\displaystyle{\cal W}_{A}:=\sum_{i\in[n]}\left(\sum_{j\in[N_{i}]}u^{i}_{j}(d^{i}_{j}+{C}^{i}_{j}-x^{i}_{j})-\sum_{j\in[G_{i}]}c^{i}_{j}(y^{i}_{j})\right), (20)

which is the same as 𝒲B{\cal W}_{B} in (8). The system operator’s economic dispatch problem is then:

max𝒲A(𝐡,𝐱,𝐝,𝐲)\displaystyle\max\qquad\mathcal{W}_{A}({\bf h,x,d,y}) (21)
subject to(5),(18).\displaystyle\text{subject to}\qquad\eqref{eq:directcons},\ \eqref{eq:aggcons_p}.
Assumption 4.

The system operator’s economic dispatch problem (21) is feasible.

As the system operator solves (21), the wholesale market prices 𝝀\boldsymbol{\lambda} are given by the optimal Lagrange multiplier of the constraint (5a). We then have the following proposition.

Proposition 3 (Competitive Equilibrium).

Under Assumptions 124, there exists an optimal solution (𝐡,𝐱,𝐝,𝐲)\left({\bf h^{*},x^{*},d^{*},y^{*}}\right) to (21). Let 𝛌{\boldsymbol{\lambda}} denote the corresponding optimal Lagrange multipliers of constraints (5a). Then, the following statements are true:

  • (𝐱,𝐝)(\bf x^{*},\bf{d}^{*}) are consistent with Lemma 3, given (𝐏,𝐩)\bf(P^{*},p^{*}) and 𝝀\boldsymbol{\lambda}.

  • (𝐏,𝐩){\bf(P^{*},p^{*})} are consistent with Lemma 4, given 𝐱{\bf x^{*}} and 𝝀\boldsymbol{\lambda}.

  • 𝐲{\bf y}^{*} satisfies (4), given 𝝀\boldsymbol{\lambda}.

We now present the following theorem, which states that our proposed aggregation model achieves the same market efficiency as the benchmark direct participation model.

Theorem 1 (Main Result).

Let 𝒲A{\cal W}_{A}^{*} be the optimal social welfare of (21), and let 𝒲B{\cal W}_{B}^{*} be the optimal social welfare of (9). Then, we have that

𝒲A=𝒲B.\displaystyle{\cal W}_{A}^{*}={\cal W}_{B}^{*}. (22)

Further, we have that the optimal 𝐱,𝐝,𝐲{\bf x^{*},d^{*},y^{*}} (from Proposition 3) solving (21) are the same as those solving (9).

In summary, under the proposed aggregation model, the aggregator procures energy from prosumers using two-part pricing, the aggregator would optimally pay the wholesale market price to the prosumers for each unit of energy procured, while the participation fee is differently charged to each prosumer as the additional consumer surplus when she sells those energy compared with not selling (which is dependent on her utility of consumption). Theorem 1, together with Proposition 3, implies that under the aggregator’s two-part differential pricing scheme, the prosumers’ optimal buying and selling behavior is exactly the same as those in the direct participation model. As a result, the social welfare achieved under the aggregation model matches that of the direct participation model, i.e., there is no loss of market efficiency from the aggregation.

Remark 1.

The significance of Theorem 1 follows from the fact that via our proposed aggregation model, DER aggregation through a profit-seeking aggregator 𝒜{\cal A} is equivalent to solving a socially-optimal economic dispatch model where 𝒜{\cal A} is absent. Hence, the potential efficiency loss due to presence of a monopolistic profit-seeking aggregator 𝒜{\cal A} is off-set by two-part pricing. This is not possible via one-part pricing, as demonstrated in [19, 20] and in Section IV, where efficiency loss arises from the profit-seeking behavior of 𝒜{\cal A}.

IV ILLUSTRATIVE EXAMPLES

Refer to caption
Figure 3: Left: Efficient aggregation vs. no DER integration. Adding more prosumers attains a higher social welfare. Middle: Comparison between the two extremes (efficient aggregation vs. no DER) and the one-part pricing model (inefficient). Right: Quantification of the efficiency loss for the one-part pricing model, via the metric PoAg.

In this section, we illustrate our results via stylized examples, with the goal of revealing some key insights. For ease of exposition, we consider a power system with 1 node, and hence, there are no network constraints. We remark that however, our theoretical results still hold for any power network. We also assume that there is a fixed demand D¯\bar{D}, and thus (5a) becomes

D+D¯YX=0.\displaystyle{D+\bar{D}-Y-X=0}. (23)

We also consider one conventional power generator with cost

c(y)=αy2+βy.c(y)=\alpha y^{2}+\beta y.

With one prosumer, consider the isoelastic utility function uu, with risk-aversion parameter η\eta [24]:

u(z)={z1η11ηη0,η1,ln(z)η=1.\displaystyle u(z)=\begin{cases}\frac{z^{1-\eta}-1}{1-\eta}&\eta\geq 0,\eta\neq 1,\\ \ln(z)&\eta=1.\end{cases} (24)

When η=0\eta=0, the prosumer is risk-neutral, for η>0\eta>0, the prosumer is risk-averse, and increasing η\eta implies more risk-aversion. For C>p1/ηC>p^{-1/\eta}, from Proposition 2, the Stackelberg equilirium prices of the game 𝒢(λ){\cal G}({\lambda}) are

p=λ,P=[λ(Cλ1/η)+u(λ1/η)u(C)]+,{p}^{*}=\lambda,\quad P^{*}=\left[\lambda\left(C-\lambda^{-1/\eta}\right)+u\left(\lambda^{-1/\eta}\right)-u(C)\right]^{+},

and the equilibrium response is

x(P,p)=Cλ1/η.{x}^{*}({P}^{*},{p}^{*})=C-\lambda^{-1/\eta}.

We note that by Lemma 3, the case when Cλ1/ηC\leq\lambda^{-1/\eta} corresponds to x=0x^{*}=0 and d=λ1/ηCd^{*}=\lambda^{-1/\eta}-C. The ISO’s problem for this single-prosumer case is then:

max𝒲A(x,d,y):=u(Cx+d)c(y)\displaystyle\max\qquad\mathcal{W}_{A}({x,d,y}):=u(C-x+d)-c(y) (25)
subject toD¯+dxy=0,0yy¯,\displaystyle\text{subject to}\qquad\bar{D}+d-x-y=0,\quad 0\leq y\leq\bar{y},
0xC,0dZC+x.\displaystyle\qquad\qquad\qquad 0\leq x\leq C,\quad 0\leq d\leq Z-C+x.

We compare the optimal value of (25), i.e., the optimal social welfare, with the following alternative models.

1) No DER participation model

In this model, the prosumer is restricted to not selling her energy, i.e., x=0x=0. The ISO’s problem is the same as (25) with the additional constraint that x=0x=0, i.e.,

max𝒲N(d,y):=u(C+d)c(y)\displaystyle\max\qquad\mathcal{W}_{N}({d,y}):=u(C+d)-c(y) (26)
subject toD¯+dy=0,0yy¯,\displaystyle\text{subject to}\qquad\bar{D}+d-y=0,\quad 0\leq y\leq\bar{y},
0dZC.\displaystyle\qquad\qquad\qquad 0\leq d\leq Z-C.

2) One-part pricing model (P=0P=0)

In this model, the prosumer may sell her energy to the aggregator 𝒜{\cal A} at a fixed marginal price pp set by 𝒜{\cal A}, with no participation fees. In this model, the prosumer solves

maxx,dπ(x,d)=u(C+dx)λd+px\displaystyle\max_{x,d}\pi(x,d)=u(C+d-x)-\lambda d+px (27)
s.t. 0xC,0dZC+x.\displaystyle\text{s.t. }\quad 0\leq x\leq C,\quad 0\leq d\leq Z-C+x.

With the isoelastic utility function (24), an optimal response of the prosumer is

d=[λ1/ηC]+andx=[Cp1/η]+.d^{*}=\left[\lambda^{-1/\eta}-C\right]^{+}\quad\text{and}\qquad x^{*}=\left[C-p^{-1/\eta}\right]^{+}.

Note that if Cλ1/ηC\leq\lambda^{-1/\eta}, then d0d^{*}\geq 0 and x=0x^{*}=0, which leads to zero profit of the aggregator for any 0pλ0\leq p\leq\lambda, leading to no DER dispatch at the wholesale market level. If C>λ1/ηC>\lambda^{-1/\eta}, then the aggregator would choose a pp such that C>p1/ηC>p^{-1/\eta}, and she solves

max0pλΠ=(λp)x(p)=(λp)(Cp1/η).\max_{0\leq p\leq\lambda}\Pi=(\lambda-p)x^{*}(p)=(\lambda-p)\left(C-p^{-1/\eta}\right).

In this case, the Stackelberg equilibrium of the aggregator-prosumer game 𝒢^(λ)\hat{{\cal G}}({\lambda}) (this is a different game from 𝒢(λ){\cal G}({\lambda}), but similarly defined) is given by (x(p),p)(x^{*}(p^{*}),p^{*}) that satisfy

(1η)p+ηCp1+1/η=λ,x(p)=Cp1/η.(1-\eta)p^{*}+\eta C{p^{*}}^{1+1/\eta}=\lambda,\quad x^{*}(p^{*})=C-p^{-1/\eta}.

If we let η=1\eta=1 (logarithmic utility), we then have

p(λ)=λC,x(p)=Cp1.p^{*}(\lambda)=\sqrt{\frac{\lambda}{C}},\quad x^{*}(p^{*})=C-p^{-1}.

Now, the next question is, how can we define the corresponding social welfare. Note that here, the Stackelberg equilibrium (x(p),p)(x^{*}(p^{*}),p^{*}) does not yield an efficient wholesale market outcome, so instead of using prosumer’s utility in ISO’s economic dispatch problem, we need to construct an induced function. First, we note that

x(p(λ))=CCλ,x^{*}(p^{*}(\lambda))=C-\sqrt{\frac{C}{\lambda}},

and hence, the corresponding inverse supply function for the prosumer is

pA(x)=C(Cx)2.p_{A}(x)=\frac{C}{(C-x)^{2}}.

Now, the system operator solves the economic dispatch problem given by

max𝒲O(x,y):=c(y)0xpA(x^)𝑑x^\displaystyle\max\qquad{\cal W}_{O}(x,y):=-c(y)-\int_{0}^{x}p_{A}(\hat{x})d\hat{x} (28)
subject toD¯xy=0,0yy¯,\displaystyle\text{subject to}\qquad\bar{D}-x-y=0,\quad 0\leq y\leq\bar{y},
0xC.\displaystyle\qquad\qquad\qquad 0\leq x\leq C.

Since there is no elastic demand dd, the above problem is equivalent to the following cost minimization problem

min𝒞O(x,y):=c(y)+xCx\displaystyle\min\qquad{\cal C}_{O}(x,y):=c(y)+\frac{x}{C-x} (29)
subject toD¯xy=0,0yy¯,\displaystyle\text{subject to}\qquad\bar{D}-x-y=0,\quad 0\leq y\leq\bar{y},
0xC.\displaystyle\qquad\qquad\qquad 0\leq x\leq C.

Anaolgously, we can define the cost 𝒞A:=𝒲A{\cal C}_{A}:=-{\cal W}_{A} for the efficient dispatch model. To capture the efficiency loss for this model, we adopt the Price of Aggregation (PoAg) metric introduced in [19, 20], which is defined as

PoAg:=𝒞O𝒞A.\text{\sf PoAg}:=\frac{{\cal C}^{*}_{O}}{{\cal C}^{*}_{A}}.

One can look at PoAg as similar to the Price-of-Anarchy in noncooperative games, but PoAg naturally captures the hierarchical nature of electricity markets. By definition, PoAg1\text{\sf PoAg}\geq 1, and when it is equal to 11, it means that the aggregation model yields a socially-optimal market outcome. Note that PoAg can be analogously defined for social welfare maximization, to compare any inefficient model to 𝒲A{\cal W}_{A}^{*}.

Figure 3 illustrates the market outcomes for the above models. In our simulations, we vary the capacity CC from 0 to D¯=100\bar{D}=100. One can think of this as a proxy for renewable energy integration, C=0C=0 implies %0 renewable integration, and C=D¯C=\bar{D} implies %100 of the total inflexible load might be met exclusively with DER capacity. We also assume the following parameters

η=1,α=0.01,β=1,Z=1000,y¯=1000.\eta=1,\ \alpha=0.01,\ \beta=1,\ Z=1000,\ \bar{y}=1000.

The above parameters are picked such that it is more expensive to use the conventional generator than DERs, and conventional generators can fully meet the total system’s demand. DER aggregation is expected to improve the social welfare, and this is demonstrated in Figure 3 (Left). Since our model is efficient, the improvements shown are the maximum possible ones. We then increase the number of prosumers to two, having the same capacity and η\eta, and observe that the efficient welfare further improves. This is natural: the more DER capacities available, the more cheaper resources are available to the ISO. Next, we fix the number of prosumers to one, and study other aspects. The one-part pricing model discussed above is also expected to improve the social welfare, compared to no DER participation, but remains inefficient, as Figure 3 (Middle) shows. We quantify the efficiency loss via the PoAg metric and plot it in Figure 3 (Right). Here, PoAg increases nonlinearly as CC increases.

V CONCLUSIONS AND FUTURE WORKS

To include DER capacities in the wholesale electricity market, we proposed an efficient DER aggregation model. The aggregator announces a two-part price offer to each prosumer (DER owner): a participation fee and a per-unit price paid to the prosumer for the energy procured. This interaction is modeled as a Stackelberg game, which admits an efficient equilibrium. Under this equilibrium, the per-unit price is equal to the wholesale market price, and prosumers choose to sell the same amount of energy as if they were participating directly in the wholesale market (ideal case). As a result, the social welfare achieved under this aggregation model is the same as the social welfare achieved when prosumers can sell directly, i.e., there is no efficiency loss even with a monopolistic aggregator. The two-part pricing mechanism is essential in achieving this efficiency result. As illustrated via an example, there may be efficiency loss if no participation fees can be charged. Under the current model, the aggregator charges different participation fees to the prosumers. There may be some debates on whether differential pricing on the participation fee should be allowed. If the aggregator has to choose a single participation fee for all prosumers at the same location, we would no longer have the efficiency result, and characterizing the efficiency loss in this aggregation process would be a promising direction to work on. Another direction for future work is to consider stochastic DER capacities under two-part pricing, and see how the stochastic nature of the DER supply [25] would affect the aggregator-prosumer game and the overall market efficiency.

References

  • [1] I. Ritzenhofen, J. R. Birge, and S. Spinler, “The structural impact of renewable portfolio standards and feed-in tariffs on electricity markets,” European Journal of Operational Research, vol. 255, no. 1, pp. 224–242, 2016.
  • [2] “Distributed Energy Resouces: Connection Modeling and Reliability Considerations,” North American Electric Reliability Corporation, Tech. Rep., 2017.
  • [3] D. Bertsimas, E. Litvinov, X. A. Sun, J. Zhao, and T. Zheng, “Adaptive robust optimization for the security constrained unit commitment problem,” IEEE Transactions on Power Systems, vol. 28, no. 1, pp. 52–63, 2013.
  • [4] D. Gan and E. Litvinov, “Energy and reserve market designs with explicit consideration to lost opportunity costs,” IEEE Transactions on Power Systems, vol. 18, no. 1, pp. 53–59, 2003.
  • [5] E. Litvinov, T. Zheng, G. Rosenwald, and P. Shamsollahi, “Marginal loss modeling in LMP calculation,” IEEE Transactions on Power Systems, vol. 19, no. 2, pp. 880–888, 2004.
  • [6] T. Zheng and E. Litvinov, “Ex post pricing in the co-optimized energy and reserve market,” IEEE Trans. Power Systems, vol. 21, no. 4, pp. 1528 – 1538, 2006.
  • [7] F. Zhao, T. Zheng, and E. Litvinov, “Constructing demand curves in forward capacity market,” IEEE Trans. Power Systems, vol. 33, no. 1, pp. 525–535, 2018.
  • [8] J. Lian, D. Wu, K. Kalsi, and H. Chen, “Theoretical framework for integrating distributed energy resources into distribution systems,” in 2017 IEEE Power Energy Society General Meeting, July 2017, pp. 1–5.
  • [9] E. Ntakou and M. Caramanis, “Price discovery in dynamic power markets with low-voltage distribution-network participants,” in 2014 IEEE PES T&D Conference and Exposition, April 2014, pp. 1–5.
  • [10] S. D. Manshadi and M. E. Khodayar, “A hierarchical electricity market structure for the smart grid paradigm,” IEEE Transactions on Smart Grid, vol. 7, no. 4, pp. 1866–1875, July 2016.
  • [11] P. M. Sotkiewicz and J. M. Vignolo, “Nodal pricing for distribution networks: efficient pricing for efficiency enhancing DG,” IEEE Transactions on Power Systems, vol. 21, no. 2, pp. 1013–1014, May 2006.
  • [12] J. Terra, R. Ferreira, C. Borges, and M. Carvalho, “Using distribution-level locational marginal pricing to value distributed generation: Impacts on revenues captured by generation agents,” in 2017 IEEE PES Innovative Smart Grid Technologies Conference - Latin America (ISGT Latin America), Sep. 2017, pp. 1–6.
  • [13] S. Huang, Q. Wu, S. S. Oren, R. Li, and Z. Liu, “Distribution locational marginal pricing through quadratic programming for congestion management in distribution networks,” IEEE Transactions on Power Systems, vol. 30, no. 4, pp. 2170–2178, July 2015.
  • [14] F. Moret and P. Pinson, “Energy collectives: a community and fairness based approach to future electricity markets,” IEEE Transactions on Power Systems, pp. 1–1, 2018.
  • [15] F. A. Rahimi and A. Ipakchi, “Transactive energy techniques: Closing the gap between wholesale and retail markets,” The Electricity Journal, vol. 25, no. 8, pp. 29 – 35, 2012.
  • [16] J. Gundlach and R. Webb, “Distributed energy resource participation in wholesale markets: Lessons from the California ISO,” Energy Law Journal, vol. 39, no. 1, pp. 47–77, 2018.
  • [17] M. Lavillotti, “DER market design: Aggregations,” New York ISO, Tech. Rep., 2018.
  • [18] FERC, “Participation of distributed energy resource aggregations in markets operated by regional transmission organizations and independent system operators,” 2020, Order No. 2222.
  • [19] K. Alshehri, M. Ndrio, S. Bose, and T. Başar, “The impact of aggregating distributed energy resources on electricity market efficiency,” Proc. 53rd Annual Conference on Information Systems and Sciences, Johns Hopkins University, Baltimore, Maryland, March 20-22, 2019.
  • [20] ——, “Quantifying market efficiency impacts of aggregated distributed energy resources,” IEEE Transactions on Power Systems, vol. 35, no. 5, 2020.
  • [21] “Distributed energy resource provider participation guide with checklist,” California ISO, Tech. Rep., 2016.
  • [22] T. Başar and G. J. Olsder, Dynamic Noncooperative Game Theory.   SIAM, 1999.
  • [23] J. R. Birge, A. Hortaçsu, and J. M. Pavlin, “Inverse optimization for the recovery of market structure from market outcomes: An application to the MISO electricity market,” Operations Research, vol. 65, no. 4, pp. 837–855, 2017.
  • [24] L. Ljungqvist and T. J. Sargent, Recursive Macroeconomic Theory.   MIT Press, 2000.
  • [25] N. Sunar and J. R. Birge, “Strategic commitment to a production schedule with uncertain supply and demand: Renewable energy in day-ahead electricity markets,” Management Science, vol. 65, no. 2, pp. 714–734, 2019.

APPENDIX

Proof of Proposition 1.

First note that (9) can be equivalently written as:

max𝒲B(𝐱𝐝,𝐲)\displaystyle\max\qquad\mathcal{W}_{B}({\bf x-d,y}) (30)
subject to 1(𝐗𝐃+𝐘)=𝟎,𝐁(𝐗𝐃+𝐘)𝐟,\displaystyle\text{subject to}\ \ {\bf 1}^{\top}\left({\bf X-D+Y}\right)=\bf 0,\ {\bf-B(X-D+Y)}\leq{\bf f},
𝐂𝐙𝐱𝐝𝐂,𝐲¯𝐲𝐲¯.\displaystyle\qquad\qquad\ \ {\bf C-Z\leq x-d\leq C},\ {\bf\underline{y}\leq y\leq\overline{y}.}

By Assumption 1 and Assumption 2, each ujiu^{i}_{j} is strictly concave, and each cjic^{i}_{j} is strictly convex. Thus, 𝒲B{\cal W}_{B} given by (8) is strictly concave. It can be seen from (30) that the feasible set is compact, and is nonempty by Assumption 3. Therefore, there exists a unique optimal solution ((𝐱𝐝),𝐲)\left({\bf(x-d)^{*},y^{*}}\right) to (30). By letting 𝐡=𝐘(𝐗𝐃){\bf h^{*}}={\bf-Y^{*}-(X-D)^{*}}, we have a unique optimal solution (𝐡,(𝐱𝐝),𝐲)\left({\bf h^{*},(x-d)^{*},y^{*}}\right) to (9).

We write the Lagrangian of (9) as

=\displaystyle{\cal L}= i[n](j[Ni]uji(Cji(xjidji))j[Gi]cji(yji)\displaystyle\sum_{i\in[n]}\Bigg{(}\sum_{j\in[N_{i}]}u^{i}_{j}\left({C}^{i}_{j}-(x^{i}_{j}-d^{i}_{j})\right)-\sum_{j\in[G_{i}]}c^{i}_{j}(y^{i}_{j})
+λi(hi+XiDi+Yi)+j[Ni](μ¯ji(Cji(xjidji))\displaystyle+\lambda^{i}\left(h^{i}+X^{i}-D^{i}+Y^{i}\right)+\sum_{j\in[N_{i}]}\Big{(}\overline{\mu}^{i}_{j}\left(C^{i}_{j}-(x^{i}_{j}-d^{i}_{j})\right)
+μ¯ji(xjidjiCji+Z))+j[Gi](ν¯ji(y¯jiyji)\displaystyle+\underline{\mu}^{i}_{j}\left(x^{i}_{j}-d^{i}_{j}-C^{i}_{j}+Z\right)\Big{)}+\sum_{j\in[G_{i}]}\Big{(}\overline{\nu}^{i}_{j}(\overline{y}^{i}_{j}-y^{i}_{j})
+ν¯ji(yjiy¯ji)))+extra terms,\displaystyle+\underline{\nu}^{i}_{j}(y^{i}_{j}-\underline{y}^{i}_{j})\Big{)}\Bigg{)}+\text{extra terms}, (31)

where λi,μ¯ji,μ¯ji,ν¯ji,ν¯ji\lambda^{i},\overline{\mu}^{i}_{j},\underline{\mu}^{i}_{j},\overline{\nu}^{i}_{j},\underline{\nu}^{i}_{j} are the Lagrange multipliers corresponding to constraints (5a) and (6), and the “extra terms” correspond to constraint (5b). Under Assumptions 1 and 2, it should be clear that the optimal solution (xjidji)Cji(x^{i}_{j}-d^{i}_{j})^{*}\neq C^{i}_{j}, (xjidji)CjiZ(x^{i}_{j}-d^{i}_{j})^{*}\neq C^{i}_{j}-Z, and yjiy¯ji{y^{i}_{j}}^{*}\neq\overline{y}^{i}_{j}, yjiy¯ji{y^{i}_{j}}^{*}\neq\underline{y}^{i}_{j}. Then, from the KKT optimality conditions, we have that μ¯ji=μ¯ji=ν¯ji=ν¯ji=0\overline{\mu}^{i}_{j}=\underline{\mu}^{i}_{j}=\overline{\nu}^{i}_{j}=\underline{\nu}^{i}_{j}=0. Further, we have that

(xjidji)\displaystyle\nabla_{(x^{i}_{j}-d^{i}_{j})}{\cal L} =uji(xjidji)|(xjidji)+λi=0,\displaystyle=\frac{\partial u^{i}_{j}}{\partial(x^{i}_{j}-d^{i}_{j})}\Big{|}_{(x^{i}_{j}-d^{i}_{j})^{*}}+\lambda^{i}=0, (32a)
yji\displaystyle\nabla_{y^{i}_{j}}{\cal L} =cjiyji|yji+λi=0,\displaystyle=-\frac{\partial c^{i}_{j}}{\partial y^{i}_{j}}\Big{|}_{{y^{i}_{j}}^{*}}+\lambda^{i}=0, (32b)

where (32a) is equivalent to (2) and (32b) is equivalent to (4). ∎

Proof of Lemma 3.

Consider prosumer jj at location ii. By Assumption 1, uji(z)z\frac{\partial u^{i}_{j}(z)}{\partial z} is continuous and ranges from \infty to 0. By Intermediate Value Theorem, for any given λi\lambda^{i} and pjip^{i}_{j}, there exist a (z1,z2)(z_{1},z_{2}) such that (12) holds. By strict concavity of ujiu^{i}_{j}, uji(z)z\frac{\partial u^{i}_{j}(z)}{\partial z} is strictly decreasing, so the (z1,z2)(z_{1},z_{2}) is unique.

We next find the optimal solution to (III-A). We look at the two cases when xji>0x^{i}_{j}>0 and xji=0x^{i}_{j}=0 separately. Recall that pjiλip^{i}_{j}\leq\lambda^{i}. Then for any xji>0x^{i}_{j}>0 and djixjid^{i}_{j}\geq x^{i}_{j}, we have that

πji(xji,dji)πji(0,djixji).\displaystyle\pi^{i}_{j}\left(x^{i}_{j},d^{i}_{j}\right)\leq\pi^{i}_{j}\left(0,d^{i}_{j}-x^{i}_{j}\right).

For any xji>0x^{i}_{j}>0 and 0<dji<xji0<d^{i}_{j}<x^{i}_{j}, we have that

πji(xji,dji)πji(xjidji,0).\displaystyle\pi^{i}_{j}\left(x^{i}_{j},d^{i}_{j}\right)\leq\pi^{i}_{j}\left(x^{i}_{j}-d^{i}_{j},0\right).

Therefore, we may without loss of generality restrict to the case when xjix^{i}_{j} and djid^{i}_{j} are not both strictly positive. We can rewrite (III-A) as

πji(xji,dji)={pjixjiPji+uji(Cjixji), if xji>0,uji(dji+Cji)λidji, if xji=0.\displaystyle\pi_{j}^{i}(x^{i}_{j},d^{i}_{j})=\begin{cases}p^{i}_{j}x^{i}_{j}-P^{i}_{j}+u^{i}_{j}\left({C}^{i}_{j}-x^{i}_{j}\right),\text{ if }x^{i}_{j}>0,\\ u^{i}_{j}\left(d^{i}_{j}+{C}^{i}_{j}\right)-\lambda^{i}d^{i}_{j},\hfill\text{ if }x^{i}_{j}=0.\end{cases} (33)
  • If Cjiz2C^{i}_{j}\leq z_{2}, then uji(Cjixji)+pjixjiuji(Cji)u^{i}_{j}\left(C^{i}_{j}-x^{i}_{j}\right)+p^{i}_{j}x^{i}_{j}\leq u^{i}_{j}\left(C^{i}_{j}\right) for any xji>0x^{i}_{j}>0. Thus we have the optimal xji=0{x^{i}_{j}}^{*}=0. If we further have that Cji<z1C^{i}_{j}<z_{1}, then the first order condition of uji(dji+Cji)λidjiu^{i}_{j}\left(d^{i}_{j}+C^{i}_{j}\right)-\lambda^{i}d^{i}_{j} leads to dji=z1Cji{d^{i}_{j}}^{*}=z_{1}-C^{i}_{j}. If we instead have Cjiz1C^{i}_{j}\geq z_{1}, then uji(dji+Cji)λidjiuji(Cji)u^{i}_{j}\left(d^{i}_{j}+C^{i}_{j}\right)-\lambda^{i}d^{i}_{j}\leq u^{i}_{j}\left(C^{i}_{j}\right) for any dji>0d^{i}_{j}>0.

  • If Cji>z2C^{i}_{j}>z_{2}, then uji(dji+Cji)λidji<uji(Cji)u^{i}_{j}\left(d^{i}_{j}+C^{i}_{j}\right)-\lambda^{i}d^{i}_{j}<u^{i}_{j}\left(C^{i}_{j}\right) for any dji>0d^{i}_{j}>0. We thus have dji=0{d^{i}_{j}}^{*}=0. When xji>0x^{i}_{j}>0, the first order condition of pjixjiPji+uji(Cjixji)p^{i}_{j}x^{i}_{j}-P^{i}_{j}+u^{i}_{j}\left({C}^{i}_{j}-x^{i}_{j}\right) leads to xji=Cjiz2{x^{i}_{j}}=C^{i}_{j}-z_{2}. Moreover, it is only optimal to have xji=Cjiz2>0{x^{i}_{j}}^{*}=C^{i}_{j}-z_{2}>0 if pjixjiPji+uji(Cjixji)uji(Cji)p^{i}_{j}x^{i}_{j}-P^{i}_{j}+u^{i}_{j}\left({C}^{i}_{j}-x^{i}_{j}\right)\geq u^{i}_{j}\left(C^{i}_{j}\right), or equivalently, Pjipji(Cjiz2)+uji(z2)uji(Cji)P^{i}_{j}\leq p^{i}_{j}\left(C^{i}_{j}-z_{2}\right)+u^{i}_{j}\left(z_{2}\right)-u^{i}_{j}\left(C^{i}_{j}\right). Therefore, we have that xji=(Cjiz2)𝟙{𝒳}{x^{i}_{j}}^{*}=\left(C^{i}_{j}-z_{2}\right)\cdot\mathds{1}\left\{{\cal X}\right\}, where 𝒳={Pjipji(Cjiz2)+uji(z2)uji(Cji)}{\cal X}=\{P^{i}_{j}\leq p^{i}_{j}\left(C^{i}_{j}-z_{2}\right)+u^{i}_{j}\left(z_{2}\right)-u^{i}_{j}\left(C^{i}_{j}\right)\}.

Proof of Lemma 4.

Consider prosumer jj at location ii. Let (z1,z2)(z_{1},z_{2}) be as in (12). Then z1z2z_{1}\leq z_{2} since pjiλip^{i}_{j}\leq\lambda^{i}. First note that when Cjiz1z2C^{i}_{j}\leq z_{1}\leq z_{2}, xji=0{x^{i}_{j}}^{*}=0 by Lemma 3, and Πji(Pji,pji)=0\Pi^{i}_{j}\left(P^{i}_{j},p^{i}_{j}\right)=0 for any (Pji,pji)+2(P^{i}_{j},p^{i}_{j})\in\mathbb{R}^{2}_{+}, i.e., the aggregator 𝒜{\cal A} earns zero profit from the prosumer regardless of her pricing decisions.

For those prosumers with Cji>z1C^{i}_{j}>z_{1}, we have that dji=0{d^{i}_{j}}^{*}=0 from Lemma 3. Since the aggregator earns zero profit if the prosumer does not sell her energy, the aggregator would choose a pjip^{i}_{j} such that Cji>z2C^{i}_{j}>z_{2} (because otherwise xji=0{x^{i}_{j}}^{*}=0) and thus xji=Cjiz2{x^{i}_{j}}^{*}=C^{i}_{j}-z_{2}. We can add in the aggregator’s problem the constraint that the prosumers would benefit from selling:

maxPji,pji0Pji+(λipji)xji\displaystyle\max_{P^{i}_{j},p^{i}_{j}\geq 0}P^{i}_{j}+\left(\lambda^{i}-p^{i}_{j}\right){x^{i}_{j}}^{*} (34)
s.t. pjixji+uji(Cjixji)Pjiuji(Cji),\displaystyle\text{s.t. }p^{i}_{j}{x^{i}_{j}}^{*}+u^{i}_{j}\left(C^{i}_{j}-{x^{i}_{j}}^{*}\right)-P^{i}_{j}\geq u^{i}_{j}\left(C^{i}_{j}\right),

where xji{x^{i}_{j}}^{*} is the optimal response of the prosumer, as a function of (Pji,pji)(P^{i}_{j},p^{i}_{j}). We observe from (34) that, the optimal PjiP^{i}_{j} should satisfy Pji=pjixji+uji(Cjixji)uji(Cji)P^{i}_{j}=p^{i}_{j}{x^{i}_{j}}^{*}+u^{i}_{j}\left(C^{i}_{j}-{x^{i}_{j}}^{*}\right)-u^{i}_{j}\left(C^{i}_{j}\right) in the maximization problem. Thus, (34) can be rewritten as

maxpji0pjixji+uji(Cjixji)uji(Cji)+(λipji)xji\displaystyle\max_{p^{i}_{j}\geq 0}p^{i}_{j}{x^{i}_{j}}^{*}+u^{i}_{j}\left(C^{i}_{j}-{x^{i}_{j}}^{*}\right)-u^{i}_{j}\left(C^{i}_{j}\right)+\left(\lambda^{i}-p^{i}_{j}\right){x^{i}_{j}}^{*}
=\displaystyle= maxz20uji(z2)uji(Cji)+λi(Cjiz2),\displaystyle\max_{z_{2}\geq 0}u^{i}_{j}(z_{2})-u^{i}_{j}\left(C^{i}_{j}\right)+\lambda^{i}\left(C^{i}_{j}-z_{2}\right), (35)

and the first order condition of (35) leads to

uji(z2)z2|z2=z2=λi.\displaystyle\frac{\partial u^{i}_{j}(z_{2})}{\partial z_{2}}\Big{|}_{z_{2}=z_{2}^{*}}=\lambda^{i}. (36)

By (12), we also have that uji(z)z|z=z2=pji\frac{\partial u^{i}_{j}(z)}{\partial z}\Big{|}_{z=z_{2}^{*}}={p^{i}_{j}}^{*}. We thus conclude that pji=λi{p^{i}_{j}}^{*}=\lambda^{i}. This further leads to z2=z1z_{2}^{*}=z_{1}, and thus

Pji\displaystyle{P^{i}_{j}}^{*} =pji(Cjiz2)+uji(z2)uji(Cji)\displaystyle={p^{i}_{j}}^{*}\left(C^{i}_{j}-z_{2}^{*}\right)+u^{i}_{j}\left(z_{2}^{*}\right)-u^{i}_{j}\left(C^{i}_{j}\right)
=λi(Cjiz1)+uji(Cjiz1)uji(Cji).\displaystyle=\lambda^{i}\left(C^{i}_{j}-z_{1}\right)+u^{i}_{j}\left(C^{i}_{j}-z_{1}\right)-u^{i}_{j}\left(C^{i}_{j}\right).

Proof of Proposition 3.

In this proof, we show that the optimal solution to (9) can be used to construct an optimal solution to (21).

First note that the objectives of both problems have the same expression, i.e., 𝒲A(𝐡,𝐱,𝐝,𝐲)=𝒲B(𝐡,𝐱𝐝,𝐲)=i[n](j[Ni]uji(dji+Cjixji)j[Gi]cji(yji)).{\cal W}_{A}({\bf h,x,d,y})={\cal W}_{B}({\bf h,x-d,y})=\sum_{i\in[n]}\left(\sum_{j\in[N_{i}]}u^{i}_{j}(d^{i}_{j}+{C}^{i}_{j}-x^{i}_{j})-\sum_{j\in[G_{i}]}c^{i}_{j}(y^{i}_{j})\right). While we took (𝐱𝐝)(\bf x-d) as a single vector of decision variables in solving (9) and obtained a unique optimal solution (𝐡,(𝐱𝐝),𝐲)\left({\bf h^{*},(x-d)^{*},y^{*}}\right) in Proposition 1, we can equivalently consider 𝐱{\bf x} and 𝐝\bf d as two separate vectors of decision variables. Consider the constraints (6) and constraints (18). For any 𝐱,𝐝,𝐲{\bf x,d,y} satisfying constraints (18), we have that 𝐱𝐝𝐂{\bf x-d\leq C} and 𝐱𝐝𝐂𝐙{\bf x-d\geq C-Z}, which implies constraints (6). With the other constraints being the same, the feasible region of (21) is a subset of the feasible region of (9), which implies that 𝒲B𝒲A{\cal W}_{B}^{*}\geq{\cal W}_{A}^{*}.

Let (𝐡,(𝐱𝐝),𝐲)\left({\bf h^{*},(x-d)^{*},y^{*}}\right) be the optimal solution to (9). Let 𝐱=[(𝐱𝐝)]+\bf x^{*}=\left[(x-d)^{*}\right]^{+} and 𝐝=[(𝐱𝐝)]+\bf d^{*}=\left[-(x-d)^{*}\right]^{+}, then 𝐱𝐝=(𝐱𝐝)\bf x^{*}-d^{*}=(x-d)^{*}. We show that (𝐡,𝐱,𝐝,𝐲)\left({\bf h^{*},x^{*},d^{*},y^{*}}\right) is feasible to (21). It suffices to show that 𝟎𝐱𝐂\bf 0\leq x^{*}\leq C and 𝟎𝐝𝐙𝐂+𝐱\bf 0\leq d^{*}\leq Z-C+x^{*}. By definition, 𝐱𝟎\bf x^{*}\geq 0 and 𝐝𝟎\bf d^{*}\geq 0. If (𝐱𝐝)<𝟎\bf(x-d)^{*}<0, then 𝐱=𝟎\bf x^{*}=0 and 𝐝=(𝐱𝐝)𝐙𝐂=𝐙𝐂+𝐱\bf d^{*}=-(x-d)^{*}\leq Z-C=Z-C+x^{*}. If (𝐱𝐝)𝟎\bf(x-d)^{*}\geq 0, then 𝐱=(𝐱𝐝)𝐂\bf x^{*}=(x-d)^{*}\leq C and 𝐝=𝟎\bf d^{*}=0. In either case, we have that 𝐱,𝐝\bf x^{*},d^{*} are feasible to (18), and thus (𝐡,𝐱,𝐝,𝐲)\left({\bf h^{*},x^{*},d^{*},y^{*}}\right) is feasible to (21). Therefore, this set of (𝐡,𝐱,𝐝,𝐲)\left({\bf h^{*},x^{*},d^{*},y^{*}}\right) is an optimal solution to (21), and thus 𝒲A=𝒲B{\cal W}_{A}^{*}={\cal W}_{B}^{*}. Further, the corresponding optimal Lagrange multipliers of constraints (5a) are the same in problem (9) and problem (21), i.e., the wholesale market prices 𝝀{\boldsymbol{\lambda}} are the same under both models.

From Proposition 1, we know that yji{y^{i}_{j}}^{*} and λi\lambda_{i} satisfy (4), i[n],j[Gi]\forall i\in[n],\forall j\in[G_{i}]. Moreover, (xjidji)\left(x^{i}_{j}-d^{i}_{j}\right)^{*} and λi\lambda^{i} satisfy (2), i[n],j[Ni]\forall i\in[n],\forall j\in[N_{i}]. Consider an arbitrary prosumer jj at an arbitrary location ii. With pji=λip^{i}_{j}=\lambda^{i}, (z1,z2)(z_{1},z_{2}) as defined in (12), and Pji=[λi(Cjiz1)+uji(z1)uji(Cji)]+{P^{i}_{j}}=\left[\lambda^{i}(C^{i}_{j}-z_{1})+u^{i}_{j}(z_{1})-u^{i}_{j}(C^{i}_{j})\right]^{+}, we have that z1=z2=Cji(xjidji)z_{1}=z_{2}=C^{i}_{j}-(x^{i}_{j}-d^{i}_{j})^{*}. If Cjiz2=Cji(xjidji)C^{i}_{j}\leq z_{2}=C^{i}_{j}-(x^{i}_{j}-d^{i}_{j})^{*}, then (xjidji)<0(x^{i}_{j}-d^{i}_{j})^{*}<0, which implies that dji=[(xjidji)]+=[z1Cji]+{d^{i}_{j}}^{*}=\left[-(x^{i}_{j}-d^{i}_{j})^{*}\right]^{+}=\left[z_{1}-C^{i}_{j}\right]^{+} and xji=[(xjidji)]+=0{x^{i}_{j}}^{*}=\left[(x^{i}_{j}-d^{i}_{j})^{*}\right]^{+}=0. If Cji>z2=Cji(xjidji)C^{i}_{j}>z_{2}=C^{i}_{j}-(x^{i}_{j}-d^{i}_{j})^{*}, then (xjidji)>0(x^{i}_{j}-d^{i}_{j})^{*}>0, which implies that dji=[(xjidji)]+=0{d^{i}_{j}}^{*}=\left[-(x^{i}_{j}-d^{i}_{j})^{*}\right]^{+}=0 and xji=[(xjidji)]+=Cjiz2{x^{i}_{j}}^{*}=\left[(x^{i}_{j}-d^{i}_{j})^{*}\right]^{+}=C^{i}_{j}-z_{2}. In either case, the optimal response xji{x^{i}_{j}}^{*} and dji{d^{i}_{j}}^{*} matches those described in Lemma 3. Further, by Lemma 4, the choice of pji=λip^{i}_{j}=\lambda^{i} and Pji=[λi(Cjiz1)+uji(z1)uji(Cji)]+{P^{i}_{j}}=\left[\lambda^{i}(C^{i}_{j}-z_{1})+u^{i}_{j}(z_{1})-u^{i}_{j}(C^{i}_{j})\right]^{+} are optimal for the aggregator to maximize her profit from the prosumer jj. We have thus verified all statements listed in Proposition 3.∎

Proof of Theorem 1.

Immediately follows from the proof of Proposition 3. ∎