Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Market Power and Withholding Behavior of Energy Storage Units

Yiqian Wu1, Bolun Xu2 and James Anderson1 1Y. Wu and J. Anderson are with the Department of Electrical Engineering, Columbia University, New York, NY 10027, USA {yiqian.wu2,james.anderson}@columbia.edu2B. Xu is with the Department of Earth and Environmental Engineering, Columbia University, New York, NY 10027, USA bx2177@columbia.edu
Abstract

Electricity markets are experiencing a rapid increase in energy storage unit participation. Unlike conventional generation resources, quantifying the competitive operation and identifying if a storage unit is exercising market power is challenging, particularly in the context of multi-interval bidding strategies. We present a framework to differentiate strategic capacity withholding behaviors attributed to market power from inherent competitive bidding in storage unit strategies. Our framework evaluates the profitability of strategic storage unit participation, analyzing bidding behaviors as both price takers and price makers using a self-scheduling model, and investigates how they leverage market inefficiencies. Specifically, we propose a price sensitivity model derived from the linear supply function equilibrium model to examine the price-anticipating bidding strategy, effectively capturing the influence of market power. We introduce a sufficient ex-post analysis for market operators to identify potential exploitative behaviors by monitoring instances of withholding within the bidding profiles, ensuring market resilience and competitiveness. We discuss and verify applicability of the proposed framework to realistic settings. Our analysis substantiates commonly observed economic bidding behaviors of storage units. Furthermore, it demonstrates that significant price volatility offers considerable profit opportunities not only for participants possessing market power but also for typical strategic profit seekers.

I INTRODUCTION

Electricity markets are seeing a surging amount of battery energy storage unit participants. Jointly facilitated by the reduced cost of battery cells and removed market participation barriers [1], battery energy storage is becoming increasingly competitive. In California, the capacity of grid-scale battery energy storage increased by ten times in three years, from less than 500 MW in 2020 to 5,000 MW in mid-2023, and is projected to reach 10,000 MW in 2025 [2]. On the other hand, participation in wholesale electricity markets to arbitrage price differences is becoming the main grid service for storage units [3], surpassing frequency regulation which is a specialized service that can only accommodate very limited storage capacity [4].

Energy storage units participate in wholesale electricity markets either by self-scheduling [5] or submitting competitive economic bids [6]. In self-scheduling, storage units design their charging and discharging schedule ahead of wholesale market clearance. In economic bidding, storage units submit separate bids to charge and discharge. In both cases, storage units often perform private optimization according to their operation characteristics, degradation cost, and opportunity cost. Unlike thermal generators that directly design bid values based on their fuel cost and heat rate curve, using optimization to design self-schedule bids is critical for storage units to systematically account for factors such as future price volatility and state-of-charge constraints [7, 8, 9].

Accompanying optimized market participation is the capacity withholding of storage units. Capacity withholding occurs when a resource is purposefully limiting its supply despite the current price being higher than its real marginal production cost. Capacity withholding is often a critical sign that a participant is exercising market power by limiting the supply of a given resource in order to drive up the price and obtain higher profit. Hence this conduct is strictly monitored and regulated [10]. However, the same regulation does not necessarily apply to energy storage units. Due to their limited energy capacity, these units must strategically target its timing to charge and discharge accounting for prices over a sequence of time, rather than focusing on the price at any single specific interval. Harvey and Hogan [11] noted that “If a unit is energy limited, its offer price will exceed the unit’s incremental cost.”, and “This economic withholding can conceptually be distinguished from the exercise of market power.”.

Nevertheless, as a regular participant in electricity markets, storage units can use the same withholding strategy to exercise market power as conventional generators. While storage units can conduct capacity withholding for justifiable motivations to seek charge and discharge arbitrage opportunities, it can further increase its withholding to exercise market power. To this end, it is extremely difficult to identify whether storage unit is exercising market power. Yet, there have not been systematic studies on the intricacies of multi-interval bidding patterns, nor have they sufficiently considered the constraints of energy-limited generation resources, including energy storage systems.

We propose a systematic framework to understand how to differentiate market power from storage bids that are legitimate withholding. The contributions of this work are two-fold.

  1. 1.

    We evaluate the profitability of strategic storage unit participants in the electricity market. More precisely, we analyze their bidding behavior through a self-scheduling profit-maximization model following two sets of market roles and bidding strategies that exploit market inefficiencies: conducting competitive arbitrage as price takers or exercising market power as price makers. Uniquely, we design a price sensitivity model derived from the linear supply function equilibrium model to examine the price-anticipating bidding strategy, offering a clear view of how the exercise of market power impacts prices.

  2. 2.

    From the perspective of the market operator, we provide insights into identifying market exploitation behavior accounting for the instances of withholding by analyzing the bidding profiles and the corresponding price series. This result proves to be a sufficient ex-post tool for market efficiency monitoring. By applying our framework, we resolve a conjecture made by Harvey and Hogan [11] pertaining to maximized energy utilization of energy limited units during the peak (or valley for energy storage units) price periods. Additionally, we reveal the considerable profit margin exposed to both price takers and price makers attributed to high price volatility. These findings are useful for maintaining market resilience and providing robust competitiveness.

The main results are validated using historical price data from NYISO.

II MODEL FORMULATION

II-A Competition and Market Power

Refer to caption
(a)
Refer to caption
(b)
Figure 1: Bidding behavior of participants as a price taker (q(λ)superscript𝑞𝜆q^{\prime}(\lambda)italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_λ ), solid line) and a price maker (q′′(λ)superscript𝑞′′𝜆q^{\prime\prime}(\lambda)italic_q start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ( italic_λ ), dashed line) and the corresponding impacts on the market outcome: (a) bid supply curves, (b) market clearing results. The axis Q𝑄Qitalic_Q measures power output. The supply function q~(λ)~𝑞𝜆\tilde{q}(\lambda)over~ start_ARG italic_q end_ARG ( italic_λ ) represents the remainder of the aggregated supply within the system. Demand is considered inelastic at D𝐷Ditalic_D, ΔqtΔsubscript𝑞𝑡\Delta q_{t}roman_Δ italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT indicates the equivalent capacity withholding. q(λ)superscript𝑞𝜆q^{\prime}(\lambda)italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_λ ) is the optimal solution to problem (3), and q′′(λ)superscript𝑞′′𝜆q^{\prime\prime}(\lambda)italic_q start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ( italic_λ ) is that to problem (4).

Competitive arbitrage and the exercise of market power are both mechanisms that leverage market inefficiencies. Arbitrage involves market participants leveraging price difference and the ability to leverage this at no risk. In contrast, the exercise of market power involves market participants deliberately reducing their production below competitive levels in order to increase market clearing prices, thereby enhancing their profitability [11]. In this work, we aim to analyze the consequences of these bidding strategies by modeling the optimal strategy making process for multi-interval bidding from the perspective of a single storage unit. The underlying theory of market power is that the market clearing price will be influenced by the offers or bids from certain participants. We explicitly model such impact on the price at time t𝑡titalic_t using a non-negative price sensitivity parameter αtsubscript𝛼𝑡\alpha_{t}italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT. The price at time t𝑡titalic_t is modeled as

λt=λ¯tαtqtsubscript𝜆𝑡subscript¯𝜆𝑡subscript𝛼𝑡subscript𝑞𝑡\displaystyle\lambda_{t}=\bar{\lambda}_{t}-\alpha_{t}q_{t}italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT (1)

where λtsubscript𝜆𝑡\lambda_{t}italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT is the influenced price, λ¯tsubscript¯𝜆𝑡\bar{\lambda}_{t}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT is the nominal clearing price and qtsubscript𝑞𝑡q_{t}italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT is the dispatch decision. It is important to clarify that λ¯tsubscript¯𝜆𝑡\bar{\lambda}_{t}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT does not represent the competitive price; rather, it signifies the price in the absence of the participation of a specific market participant.111A similar price sensitivity model is introduced in [12]. However, the model in [12] quantifies price impact via the amount of withheld production, λt=λ~t+αtΔqtsubscript𝜆𝑡subscript~𝜆𝑡subscript𝛼𝑡Δsubscript𝑞𝑡\lambda_{t}=\tilde{\lambda}_{t}+\alpha_{t}\Delta q_{t}italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = over~ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT roman_Δ italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT. This approach does not capture the immediate price drop following the entry of a price maker. Additionally, it calculates withheld production based on the competitive level, which complicates the analysis.

Equation (1) captures the fact that, as a market participant possessing market power, their participation in the market (qt0subscript𝑞𝑡0q_{t}\neq 0italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ≠ 0) naturally brings down the market price (λt<λ¯tsubscript𝜆𝑡subscript¯𝜆𝑡\lambda_{t}<\bar{\lambda}_{t}italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT). This effect is particularly notable when the participant serves as a large-scale pivotal supplier, whose power capacity is essential for easing potential power shortages in the system. Hence, the nominal price λ¯tsubscript¯𝜆𝑡\bar{\lambda}_{t}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT functions as an indicator of the system’s load level or supply capacity. Additionally, offering more power to the market lowers the price. Therefore, practices such as physical capacity withholding (qt<P¯subscript𝑞𝑡¯𝑃q_{t}<\bar{P}italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG, where P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG denotes the power capacity) can raise market prices, generating a profit margin for the withholding entity and possibly others but also diminishing social welfare by affecting the overall market equilibrium.

To understand the bidding behavior of participants anticipating their impact on prices and the consequent market outcomes, we consider a bid-based market that clears by relying on affine supply functions [13]:

qt(λt)=atλt+btsubscript𝑞𝑡subscript𝜆𝑡subscript𝑎𝑡subscript𝜆𝑡subscript𝑏𝑡\displaystyle q_{t}(\lambda_{t})=a_{t}\lambda_{t}+b_{t}italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) = italic_a start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT (2)

where the parameter pair (at,btsubscript𝑎𝑡subscript𝑏𝑡a_{t},b_{t}italic_a start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT) characterizes the offer chosen by the participant for time interval t𝑡titalic_t.

Generally, a rational market participant designs their bid supply function qt(λt)subscript𝑞𝑡subscript𝜆𝑡q_{t}(\lambda_{t})italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) by solving a profit-maximization problem using price forecast λ^tsubscript^𝜆𝑡\hat{\lambda}_{t}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT over a future time period t𝒯𝑡𝒯t\in~{}\mathcal{T}italic_t ∈ caligraphic_T:

qt(λt)=argmaxat,btπ(qt(λt);λ^t)s.t.(2)formulae-sequencesubscript𝑞𝑡subscript𝜆𝑡subscript𝑎𝑡subscript𝑏𝑡argmax𝜋subscript𝑞𝑡subscript𝜆𝑡subscript^𝜆𝑡stitalic-(2italic-)\displaystyle q_{t}(\lambda_{t})=\underset{a_{t},b_{t}}{\mathrm{argmax}}\ \pi(% q_{t}(\lambda_{t});\hat{\lambda}_{t})\quad\mathrm{s.t.}\ \eqref{eq:bidding_curve}italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) = start_UNDERACCENT italic_a start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_UNDERACCENT start_ARG roman_argmax end_ARG italic_π ( italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) ; over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) roman_s . roman_t . italic_( italic_) (3)

where π=t𝒯λ^tqt(λt)C(qt(λt))𝜋subscript𝑡𝒯subscript^𝜆𝑡subscript𝑞𝑡subscript𝜆𝑡𝐶subscript𝑞𝑡subscript𝜆𝑡\pi=\sum_{t\in\mathcal{T}}\hat{\lambda}_{t}q_{t}(\lambda_{t})-C\left(q_{t}(% \lambda_{t})\right)italic_π = ∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) - italic_C ( italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) ) represents the total profit and C()𝐶C(\cdot)italic_C ( ⋅ ) is a given strictly convex cost function.

Ideally, the supply function qt(λt)subscript𝑞𝑡subscript𝜆𝑡q_{t}(\lambda_{t})italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) should reflect the actual operation costs of production in a competitive market. However, bids that incorporate price sensitivity, as determined by solving

qt(λt)=argmaxat,btπ(qt(λt),λ¯tαqt(λt))s.t.(2)formulae-sequencesubscript𝑞𝑡subscript𝜆𝑡subscript𝑎𝑡subscript𝑏𝑡argmax𝜋subscript𝑞𝑡subscript𝜆𝑡subscript¯𝜆𝑡𝛼subscript𝑞𝑡subscript𝜆𝑡stitalic-(2italic-)\displaystyle q_{t}(\lambda_{t})=\underset{a_{t},b_{t}}{\mathrm{argmax}}\ \pi(% q_{t}(\lambda_{t}),\bar{\lambda}_{t}-\alpha q_{t}(\lambda_{t}))\quad\mathrm{s.% t.}\ \eqref{eq:bidding_curve}italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) = start_UNDERACCENT italic_a start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_UNDERACCENT start_ARG roman_argmax end_ARG italic_π ( italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) , over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_α italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) ) roman_s . roman_t . italic_( italic_) (4)

may obscure the true cost information. Assuming a quadratic cost function C()𝐶C(\cdot)italic_C ( ⋅ ), we illustrate the variation in bidding behavior and the consequent market outcome in Fig. 1 (for simplicity, the subscript t𝑡titalic_t in the bidding functions is neglected in the figures). q(λ)superscript𝑞𝜆q^{\prime}(\lambda)italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_λ ) is the optimal solution to problem (3), and q′′(λ)superscript𝑞′′𝜆q^{\prime\prime}(\lambda)italic_q start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ( italic_λ ) is that to problem (4). As shown in Fig. 1(1(a)), the bidding curve shifts from q(λ)superscript𝑞𝜆q^{\prime}(\lambda)italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_λ ) to q′′(λ)superscript𝑞′′𝜆q^{\prime\prime}(\lambda)italic_q start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT ( italic_λ ) due to the consideration of price sensitivity. This shift results in the market clearing price increasing from λtsubscriptsuperscript𝜆𝑡\lambda^{\prime}_{t}italic_λ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT to λt′′subscriptsuperscript𝜆′′𝑡\lambda^{\prime\prime}_{t}italic_λ start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT, while maintaining the same level of power output qtsubscript𝑞𝑡q_{t}italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT, indicating a clear reduction in economic efficiency. This alteration in the bidding curve can be seen as cost-related withholding, i.e., economic withholding. The market clearing outcome related to this behavior is shown in Fig. 1(1(b)). Here, we use q~(λ)~𝑞𝜆\tilde{q}(\lambda)over~ start_ARG italic_q end_ARG ( italic_λ ) to denote the remainder of the aggregated supply function within the system. It can be seen that, as two primary approaches of conducting capacity withholding or further exercising market power [11], economic withholding through the bid supply function is equivalent to physical withholding of the amount ΔqtΔsubscript𝑞𝑡\Delta q_{t}roman_Δ italic_q start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT from the perspective of market outcomes. For simplicity, we will focus on the effects of withholding practices and refer to this as capacity withholding throughout the remainder of the paper, without loss of generality.

Additionally, it can be inferred from Fig. 1(1(b)) that the determination of the price sensitivity parameter αtsubscript𝛼𝑡\alpha_{t}italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT is influenced by the slopes of the system supply function (q~(λ)~𝑞𝜆\tilde{q}(\lambda)over~ start_ARG italic_q end_ARG ( italic_λ ) in the figure), the demand function (D𝐷Ditalic_D in the figure), and the unit supply function. Although delving into the specifics of this process extends beyond the scope of this work, it highlights the necessity for market participants with market power to obtain the knowledge of the system demand level and the bids from other participants to accurately estimate their influence on market outcomes and making corresponding bidding decisions. Interested readers may refer to [12] for more details.

We define market participants based on their bidding strategies and their impacts on the market clearing prices as follows:

Definition 1 (Price Taker)

A market participant is a price taker if it accepts the existing prices as given and lacks the market share to influence market prices on its own.

Definition 2 (Price Maker)

A market participant is a price maker if they anticipate the influence of their bids on market prices. A price maker has sufficient knowledge of the system status, such as demand levels and the bids of other participants.

A more thorough analysis of the bidding strategies associated with these two categories of market participants will be discussed in later sections.

Remark 1

Our price sensitivity model can be generalized to scenarios where the bidding entity operates with other forms of cost functions, or participates in the market with constant cost bids or quantity bids.

II-B Strategic Bidding of Energy Storage Units

We adopt a convex self-scheduling model to characterize the bidding strategies of energy storage units, utilizing a series of price point forecasts λ^tsubscript^𝜆𝑡\hat{\lambda}_{t}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT for all t𝒯𝑡𝒯t\in\mathcal{T}italic_t ∈ caligraphic_T, where 𝒯={1,2,,T}𝒯12𝑇\mathcal{T}=\{1,2,\ldots,T\}caligraphic_T = { 1 , 2 , … , italic_T } [14]. The profit from future intervals can be regarded as the opportunity value at the current time interval. Most existing market designs require storage unit operators to make dispatch decisions based on future price estimates. This strategy reflects the allocation of limited output into a profile that maximizes total profit over the scheduling period. The solution to the profit-maximization problem forms the control policy or decision of the storage units. For ease of exposition, we assume that the storage unit accurately estimates their impact on market prices with αtsubscript𝛼𝑡\alpha_{t}italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT:

maximizept,bt,etsubscript𝑝𝑡subscript𝑏𝑡subscript𝑒𝑡maximize\displaystyle\underset{p_{t},b_{t},e_{t}}{\text{maximize}}start_UNDERACCENT italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_UNDERACCENT start_ARG maximize end_ARG t𝒯λ^t(ptbt)subscript𝑡𝒯subscript^𝜆𝑡subscript𝑝𝑡subscript𝑏𝑡\displaystyle\quad\sum_{t\in\mathcal{T}}\hat{\lambda}_{t}(p_{t}-b_{t})∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) (5a)
s.t.formulae-sequencest\displaystyle\mathrm{s.t.}roman_s . roman_t . 0pt,btP¯,t𝒯formulae-sequence0subscript𝑝𝑡formulae-sequencesubscript𝑏𝑡¯𝑃for-all𝑡𝒯\displaystyle\quad 0\leq p_{t},b_{t}\leq\bar{P},\,\quad\forall t\in\mathcal{T}0 ≤ italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ≤ over¯ start_ARG italic_P end_ARG , ∀ italic_t ∈ caligraphic_T (5b)
pt=0 if λ^t<0,t𝒯formulae-sequencesubscript𝑝𝑡0 if subscript^𝜆𝑡0for-all𝑡𝒯\displaystyle\quad p_{t}=0\textrm{ if }\hat{\lambda}_{t}<0,\,\quad\forall t\in% \mathcal{T}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = 0 if over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < 0 , ∀ italic_t ∈ caligraphic_T (5c)
etet1=ptη+btη,t𝒯formulae-sequencesubscript𝑒𝑡subscript𝑒𝑡1subscript𝑝𝑡𝜂subscript𝑏𝑡𝜂for-all𝑡𝒯\displaystyle\quad e_{t}-e_{t-1}=-\frac{p_{t}}{\eta}+b_{t}\eta,\,\quad\forall t% \in\mathcal{T}italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_e start_POSTSUBSCRIPT italic_t - 1 end_POSTSUBSCRIPT = - divide start_ARG italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG start_ARG italic_η end_ARG + italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_η , ∀ italic_t ∈ caligraphic_T (5d)
0etE,t𝒯formulae-sequence0subscript𝑒𝑡𝐸for-all𝑡𝒯\displaystyle\quad 0\leq e_{t}\leq E,\,\quad\forall t\in\mathcal{T}0 ≤ italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ≤ italic_E , ∀ italic_t ∈ caligraphic_T (5e)

where ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT and btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT denote the amount of energy discharged and charged respectively over time interval t𝑡titalic_t. We refer to intervals when 0<pt<P¯0subscript𝑝𝑡¯𝑃0<p_{t}<\bar{P}0 < italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG or 0<bt<P¯0subscript𝑏𝑡¯𝑃0<b_{t}<\bar{P}0 < italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG for any t𝒯𝑡𝒯t\in\mathcal{T}italic_t ∈ caligraphic_T as withholding intervals. Objective (5a) represents the total profit π𝜋\piitalic_π over the future period t𝒯𝑡𝒯t\in\mathcal{T}italic_t ∈ caligraphic_T. For price makers, it follows from (1) that the price forecast takes the form based on the price sensitivity: λ^t=λ¯tαt(ptbt)subscript^𝜆𝑡subscript¯𝜆𝑡subscript𝛼𝑡subscript𝑝𝑡subscript𝑏𝑡\hat{\lambda}_{t}=\bar{\lambda}_{t}-\alpha_{t}(p_{t}-b_{t})over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ). Note that incorporating both forms of the price for profit calculation in (5a) preserves the convexity of model (5). Constraint (5b) captures the charging/discharging power lower and upper bounds. Constraint (5c) ensures that the storage unit does not discharge during periods of negative pricing, a sufficient condition to preclude simultaneous charging and discharging [14]. The inter-temporal relationship of the state of charge (SoC) etsubscript𝑒𝑡e_{t}italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT is defined in (5d) with charging and discharging efficiency parameter η(0,1]𝜂01\eta\in(0,1]italic_η ∈ ( 0 , 1 ]. Inequality (5e) models energy storage capacity. We define storage units as operating in idle scenarios if they remain on standby according to model (5) throughout the scheduling period 𝒯𝒯\mathcal{T}caligraphic_T, i.e., pt=0,bt=0formulae-sequencesubscript𝑝𝑡0subscript𝑏𝑡0p_{t}=0,~{}b_{t}=0italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = 0 , italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = 0, for all t𝒯𝑡𝒯t\in\mathcal{T}italic_t ∈ caligraphic_T.

Remark 2

Although the framework for energy storage unit bidding in this paper is based on self-scheduling and neglects the operating cost C()𝐶C(\cdot)italic_C ( ⋅ ) for profit calculation, it can be easily adapted to allow economic bidding [5], accounting for the costs associated with storage discharge, such as those related to unit degradation, operation and maintenance, in its profit calculation. The results discussed in this paper would still apply in this extended context.

Note that our strategic bidding model is deterministic. We have intentionally ignored stochastic effects in order to simplify our analysis. Introducing uncertainty would necessitate more sophisticated bidding strategies which would unnecessarily complicate the analysis and obfuscate the underlying mechanisms at play. Consequently, we plan to explore these more realistic settings in future work.

III MAIN RESULTS

We now present our framework for energy storage bidding, demonstrating the distinction between a storage unit exercising market power and competitive capacity withholding. In the sequel, we examine the bidding strategies of both price takers and price makers.

Theorem 1

Assume the energy storage unit behaves rationally and designs its bid by solving the profit-maximization problem (5) using price forecast λ^tsubscript^𝜆𝑡\hat{\lambda}_{t}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT over a horizon of length NT𝑁𝑇N*Titalic_N ∗ italic_T for N𝑁Nitalic_N bidding scheduling periods, with each period of length T𝑇Titalic_T. Given a series of observed storage power output profiles {pt,bt}subscript𝑝𝑡subscript𝑏𝑡\{p_{t},~{}b_{t}\}{ italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT } and market clearing prices {λt}subscript𝜆𝑡\{{\lambda}_{t}\}{ italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT } for all t𝒯~𝑡~𝒯t\in\tilde{\mathcal{T}}italic_t ∈ over~ start_ARG caligraphic_T end_ARG, where 𝒯~={1,2,,NT}~𝒯12𝑁𝑇\tilde{\mathcal{T}}=\{1,2,\ldots,NT\}over~ start_ARG caligraphic_T end_ARG = { 1 , 2 , … , italic_N italic_T }, the storage unit is not evidently exercising market power, if the following conditions are satisfied:

  1. 1.

    The number of withholding intervals is no larger than the number of non-idle scheduling periods Nsuperscript𝑁N^{\prime}italic_N start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT, i.e., t𝒯~𝟙{0<pt<P¯}+t𝒯~𝟙{0<bt<P¯}NNsubscript𝑡~𝒯subscript10subscript𝑝𝑡¯𝑃subscript𝑡~𝒯subscript10subscript𝑏𝑡¯𝑃superscript𝑁𝑁\sum_{t\in\tilde{\mathcal{T}}}\mathds{1}_{\{0<p_{t}<\bar{P}\}}+\sum_{t\in% \tilde{\mathcal{T}}}\mathds{1}_{\{0<b_{t}<\bar{P}\}}\leq N^{\prime}\leq N∑ start_POSTSUBSCRIPT italic_t ∈ over~ start_ARG caligraphic_T end_ARG end_POSTSUBSCRIPT blackboard_1 start_POSTSUBSCRIPT { 0 < italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG } end_POSTSUBSCRIPT + ∑ start_POSTSUBSCRIPT italic_t ∈ over~ start_ARG caligraphic_T end_ARG end_POSTSUBSCRIPT blackboard_1 start_POSTSUBSCRIPT { 0 < italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG } end_POSTSUBSCRIPT ≤ italic_N start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ≤ italic_N,

  2. 2.

    The price-decision relationship is consistent with Proposition 4.

Proposition 4 reveals the relationship between the control decisions of a price taker and the prevailing market prices. Details of Proposition 4 will be provided in subsequent subsections.

Theorem 1 suggests that the market operator can conduct ex-post analysis according to the observed storage power output profiles and market clearing price series to examine the exercise of market power with the aim of monitoring market efficiency. This method aligns with the “after-the-fact standard” approach for analysis as recommended by Harvey and Hogan in [11]. Note that from the market operator’s perspective, that actual competitive bids and the resulting clearing prices are not accessible. The available information for conducting market efficiency analysis comprises solely the effective storage power output profiles and the corresponding prices. The presence of partial operation intervals during the observation period indicates instances of storage unit conducting capacity withholding. Merely counting these instances below a predefined threshold and examining the prices associated with these instances offers a practical and efficient method to differentiate the causes of capacity withholding – whether or not it results from market power exercise or competitive behavior.

The proof of Theorem 1 will compromise an examination of bidding strategies derived from simplified individual profit maximization problems discussed in subsequent subsections. Strictly speaking, we should directly deal with the bidding model (5); however, we will initially introduce a set of assumptions to facilitate preliminary discussions and propositions. These assumptions will later be revisited and removed, allowing us to generalize the conditions and validate the applicability of Theorem 1.

In the simplified profit maximization model, we focus on scenarios with short clearing periods, e.g., T=2,3𝑇23T=2,3italic_T = 2 , 3. Given the short duration of storage charging and discharging, we assume the storage capacity constraint (5e) is strictly satisfied throughout the scheduling period 𝒯𝒯\mathcal{T}caligraphic_T and neglect it in the model formulation:

Assumption 1 (Unsaturated SoC)

During the clearing period, the energy storage SoC remains strictly feasible, i.e., 0<et<E0subscript𝑒𝑡𝐸0<e_{t}<E0 < italic_e start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < italic_E for all t𝒯𝑡𝒯t\in\mathcal{T}italic_t ∈ caligraphic_T.

In light of Assumption 1, the SoC evolution constraint (5d) reduces to:

t𝒯ptηbtη=0.subscript𝑡𝒯subscript𝑝𝑡𝜂subscript𝑏𝑡𝜂0\displaystyle\sum_{t\in\mathcal{T}}\frac{p_{t}}{\eta}-b_{t}\eta=0.∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT divide start_ARG italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG start_ARG italic_η end_ARG - italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_η = 0 . (6)
Assumption 2 (Positive profit operation)

The storage unit remains idle, i.e., pt=0,bt=0formulae-sequencesubscript𝑝𝑡0subscript𝑏𝑡0p_{t}=0,~{}b_{t}=0italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = 0 , italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = 0, for all t𝒯𝑡𝒯t\in\mathcal{T}italic_t ∈ caligraphic_T, if the induced profit is zero, i.e., π=0𝜋0\pi=0italic_π = 0.

Assumption 3 (Non-negative pricing)

Prices remain non-negative, λt0subscript𝜆𝑡0\lambda_{t}\geq 0italic_λ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ≥ 0, for all t𝒯𝑡𝒯t\in\mathcal{T}italic_t ∈ caligraphic_T.

Based on Assumption 3, we omit the constraint (5c) in the bidding model, however the non-simultaneous charging and discharging condition does still hold.

III-A Withholding as Price Taker

First, we aim to analyze the withholding bidding behavior assuming that the energy storage unit participates in the market as a price taker. From Assumptions 1-3, the corresponding bidding model is:

maximizept,btsubscript𝑝𝑡subscript𝑏𝑡maximize\displaystyle\underset{p_{t},b_{t}}{\text{maximize}}start_UNDERACCENT italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_UNDERACCENT start_ARG maximize end_ARG t𝒯λ^t(ptbt)subscript𝑡𝒯subscript^𝜆𝑡subscript𝑝𝑡subscript𝑏𝑡\displaystyle\quad\sum_{t\in\mathcal{T}}\hat{\lambda}_{t}(p_{t}-b_{t})∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) (7a)
s.t.formulae-sequencest\displaystyle\mathrm{s.t.}roman_s . roman_t . t𝒯ptηbtη=0:θ:subscript𝑡𝒯subscript𝑝𝑡𝜂subscript𝑏𝑡𝜂0𝜃\displaystyle\quad\sum_{t\in\mathcal{T}}\frac{p_{t}}{\eta}-b_{t}\eta=0:\,\theta∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT divide start_ARG italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG start_ARG italic_η end_ARG - italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_η = 0 : italic_θ (7b)
0ptP¯:δt,δt+,t𝒯:0subscript𝑝𝑡¯𝑃superscriptsubscript𝛿𝑡superscriptsubscript𝛿𝑡for-all𝑡𝒯\displaystyle\quad 0\leq p_{t}\leq\bar{P}:\,\delta_{t}^{-},\delta_{t}^{+},\,% \forall t\in\mathcal{T}0 ≤ italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ≤ over¯ start_ARG italic_P end_ARG : italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , ∀ italic_t ∈ caligraphic_T (7c)
0btP¯:βt,βt+,t𝒯:0subscript𝑏𝑡¯𝑃superscriptsubscript𝛽𝑡superscriptsubscript𝛽𝑡for-all𝑡𝒯\displaystyle\quad 0\leq b_{t}\leq\bar{P}:\,\beta_{t}^{-},\beta_{t}^{+},\,% \forall t\in\mathcal{T}0 ≤ italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ≤ over¯ start_ARG italic_P end_ARG : italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , ∀ italic_t ∈ caligraphic_T (7d)

where the corresponding dual variable is defined after each constraint. Note that although we formulate the bidding model of price takers using λ^tsubscript^𝜆𝑡\hat{\lambda}_{t}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT in the objective (7a), in the presence of price makers, λ^tsubscript^𝜆𝑡\hat{\lambda}_{t}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT will be replaced by the influenced price based on  (1) to determine the ultimate dispatch decision and profit.

The Lagrangian function of (7) is:

L=𝐿absent\displaystyle L=italic_L = t𝒯(λ^t(ptbt)+θ(ptηbtη)\displaystyle\sum_{t\in\mathcal{T}}\Big{(}\hat{\lambda}_{t}(p_{t}-b_{t})+% \theta(\frac{p_{t}}{\eta}-b_{t}\eta)∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT ( over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) + italic_θ ( divide start_ARG italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG start_ARG italic_η end_ARG - italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_η ) (8)
+δtptδt+(ptP¯)+βtbtβt+(btP¯)),\displaystyle+\delta_{t}^{-}p_{t}-\delta_{t}^{+}(p_{t}-\bar{P})+\beta_{t}^{-}b% _{t}-\beta_{t}^{+}(b_{t}-\bar{P})\Big{)},+ italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ( italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - over¯ start_ARG italic_P end_ARG ) + italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ( italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - over¯ start_ARG italic_P end_ARG ) ) ,

and the corresponding KKT conditions are (t𝒯for-all𝑡𝒯\forall t\in\mathcal{T}∀ italic_t ∈ caligraphic_T):

Lpt=λ^t+θη+δtδt+=0𝐿subscript𝑝𝑡subscript^𝜆𝑡𝜃𝜂superscriptsubscript𝛿𝑡superscriptsubscript𝛿𝑡0\displaystyle\frac{\partial L}{\partial p_{t}}=\hat{\lambda}_{t}+\frac{\theta}% {\eta}+\delta_{t}^{-}-\delta_{t}^{+}=0divide start_ARG ∂ italic_L end_ARG start_ARG ∂ italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG = over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG + italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT - italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = 0 (9a)
Lbt=λ^tθη+βtβt+=0𝐿subscript𝑏𝑡subscript^𝜆𝑡𝜃𝜂superscriptsubscript𝛽𝑡superscriptsubscript𝛽𝑡0\displaystyle\frac{\partial L}{\partial b_{t}}=-\hat{\lambda}_{t}-\theta\eta+% \beta_{t}^{-}-\beta_{t}^{+}=0divide start_ARG ∂ italic_L end_ARG start_ARG ∂ italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG = - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_θ italic_η + italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT - italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = 0 (9b)
t𝒯ptηbtη=0subscript𝑡𝒯subscript𝑝𝑡𝜂subscript𝑏𝑡𝜂0\displaystyle\sum_{t\in\mathcal{T}}\frac{p_{t}}{\eta}-b_{t}\eta=0∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT divide start_ARG italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG start_ARG italic_η end_ARG - italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_η = 0 (9c)
0pt,btP¯formulae-sequence0subscript𝑝𝑡subscript𝑏𝑡¯𝑃\displaystyle 0\leq p_{t},~{}b_{t}\leq\bar{P}0 ≤ italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ≤ over¯ start_ARG italic_P end_ARG (9d)
δt,δt+,βt,βt+0superscriptsubscript𝛿𝑡superscriptsubscript𝛿𝑡superscriptsubscript𝛽𝑡superscriptsubscript𝛽𝑡0\displaystyle\delta_{t}^{-},~{}\delta_{t}^{+},~{}\beta_{t}^{-},~{}\beta_{t}^{+% }~{}\geq 0italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ≥ 0 (9e)
δtpt=0,δt+(ptP¯)=0formulae-sequencesuperscriptsubscript𝛿𝑡subscript𝑝𝑡0superscriptsubscript𝛿𝑡subscript𝑝𝑡¯𝑃0\displaystyle\delta_{t}^{-}p_{t}=0,\,\delta_{t}^{+}(p_{t}-\bar{P})=0italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = 0 , italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ( italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - over¯ start_ARG italic_P end_ARG ) = 0 (9f)
βtbt=0,βt+(btP¯)=0formulae-sequencesuperscriptsubscript𝛽𝑡subscript𝑏𝑡0superscriptsubscript𝛽𝑡subscript𝑏𝑡¯𝑃0\displaystyle\beta_{t}^{-}b_{t}=0,\,\beta_{t}^{+}(b_{t}-\bar{P})=0italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = 0 , italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT ( italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - over¯ start_ARG italic_P end_ARG ) = 0 (9g)

where Eqs. (9a) and (9b) correspond to stationarity conditions, constraints (9c) and (9d) to primal feasibility, constraint (9e) to dual feasibility, Eqs. (9f) and (9g) to complementary slackness. Given that the optimization model (7) is a linear program, there is zero duality gap, yielding KKT conditions both necessary and sufficient for optimality.

Proposition 1

For a strategic price taker making bidding decisions based on model (7), the bidding decisions {pt,bt}subscriptsuperscript𝑝𝑡subscriptsuperscript𝑏𝑡\{p^{*}_{t},~{}b^{*}_{t}\}{ italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT } throughout the period 𝒯𝒯\mathcal{T}caligraphic_T will include at least one interval t𝑡titalic_t at capacity, i.e., pt=P¯subscriptsuperscript𝑝𝑡¯𝑃p^{*}_{t}=\bar{P}italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG or bt=P¯subscriptsuperscript𝑏𝑡¯𝑃b^{*}_{t}=\bar{P}italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG, except for the idle scenarios.

Proof:

Suppose the unit withholds a certain amount of capacity throughout the period 𝒯𝒯\mathcal{T}caligraphic_T: 0<pt,bt<P¯,t𝒯𝒯formulae-sequence0subscriptsuperscript𝑝𝑡formulae-sequencesubscriptsuperscript𝑏𝑡¯𝑃for-all𝑡superscript𝒯𝒯0<p^{*}_{t},~{}b^{*}_{t}<\bar{P},\,\forall t\in\mathcal{T}^{\prime}\subseteq% \mathcal{T}0 < italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG , ∀ italic_t ∈ caligraphic_T start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ⊆ caligraphic_T, where 𝒯superscript𝒯\mathcal{T}^{\prime}caligraphic_T start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT is the set of non-idle intervals. According to the complementary constraints (9f) and (9g), we have δt+,βt+=0,t𝒯formulae-sequencesuperscriptsuperscriptsubscript𝛿𝑡superscriptsuperscriptsubscript𝛽𝑡0for-all𝑡superscript𝒯{\delta_{t}^{+}}^{*},~{}{\beta_{t}^{+}}^{*}=0,\,\forall t\in\mathcal{T}^{\prime}italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = 0 , ∀ italic_t ∈ caligraphic_T start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT. The non-idle condition suggests that there exists ti,tj𝒯subscript𝑡𝑖subscript𝑡𝑗superscript𝒯t_{i},t_{j}\in\mathcal{T}^{\prime}italic_t start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT , italic_t start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ∈ caligraphic_T start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT, that pti>0superscriptsubscript𝑝subscript𝑡𝑖0p_{t_{i}}^{*}>0italic_p start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT > 0 and btj>0superscriptsubscript𝑏subscript𝑡𝑗0b_{t_{j}}^{*}>0italic_b start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT > 0, thus we have δti|pti>0=0evaluated-atsuperscriptsuperscriptsubscript𝛿subscript𝑡𝑖superscriptsubscript𝑝subscript𝑡𝑖00{\delta_{t_{i}}^{-}}^{*}|_{p_{t_{i}}^{*}>0}=0italic_δ start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT | start_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT > 0 end_POSTSUBSCRIPT = 0 and βtj|btj>0=0evaluated-atsuperscriptsuperscriptsubscript𝛽subscript𝑡𝑗superscriptsubscript𝑏subscript𝑡𝑗00{\beta_{t_{j}}^{-}}^{*}|_{b_{t_{j}}^{*}>0}=0italic_β start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT | start_POSTSUBSCRIPT italic_b start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT > 0 end_POSTSUBSCRIPT = 0 based on (9f) and (9g), then the corresponding stationary conditions (9a), (9b) for these intervals are

{λ^ti+θη=0,ti{ti|pti>0,ti𝒯},λ^tjθη=0,tj{tj|btj>0,tj𝒯}.\displaystyle\left\{\begin{aligned} &\hat{\lambda}_{t_{i}}+\frac{\theta}{\eta}% =0,\forall t_{i}\in\{t_{i}~{}|~{}p_{t_{i}}^{*}>0,t_{i}\in\mathcal{T}^{\prime}% \},\\ &-\hat{\lambda}_{t_{j}}-\theta\eta=0,\forall t_{j}\in\{t_{j}~{}|~{}b_{t_{j}}^{% *}>0,t_{j}\in\mathcal{T}^{\prime}\}.\end{aligned}\right.{ start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG = 0 , ∀ italic_t start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∈ { italic_t start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT | italic_p start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT > 0 , italic_t start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ∈ caligraphic_T start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT } , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUBSCRIPT - italic_θ italic_η = 0 , ∀ italic_t start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ∈ { italic_t start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT | italic_b start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT > 0 , italic_t start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ∈ caligraphic_T start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT } . end_CELL end_ROW (10)

Clearly, solutions to (10) need to satisfy λ^tiη2=λ^tjsubscript^𝜆subscript𝑡𝑖superscript𝜂2subscript^𝜆subscript𝑡𝑗\hat{\lambda}_{t_{i}}\eta^{2}=\hat{\lambda}_{t_{j}}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUBSCRIPT. Given the profit is calculated as (7a), where ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT and btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT satisfy (7b), prices satisfying λ^tiη2=λ^tjsubscript^𝜆subscript𝑡𝑖superscript𝜂2subscript^𝜆subscript𝑡𝑗\hat{\lambda}_{t_{i}}\eta^{2}=\hat{\lambda}_{t_{j}}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT end_POSTSUBSCRIPT result in a profit π=0𝜋0\pi=0italic_π = 0 . Based on Assumption 2, the considered intervals keep idle, which violates the non-idle condition. ∎

Proposition 2

For a strategic price taker making bidding decisions based on model (7), assuming the storage unit is not idle throughout the scheduling period 𝒯𝒯\mathcal{T}caligraphic_T, then

  1. 1.

    there must exist λ^min<λ^maxsubscript^𝜆minsubscript^𝜆max\hat{\lambda}_{\mathrm{min}}<~{}\hat{\lambda}_{\mathrm{max}}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT < over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT, where λ^min=mint𝒯{λ^t}subscript^𝜆minsubscriptmin𝑡𝒯subscript^𝜆𝑡\hat{\lambda}_{\mathrm{min}}=\mathrm{min}_{t\in\mathcal{T}}\{\hat{\lambda}_{t}\}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT = roman_min start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT { over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT } and λ^max=maxt𝒯{λ^t}subscript^𝜆maxsubscriptmax𝑡𝒯subscript^𝜆𝑡\hat{\lambda}_{\mathrm{max}}=\mathrm{max}_{t\in\mathcal{T}}\{\hat{\lambda}_{t}\}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT = roman_max start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT { over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT },

  2. 2.

    the storage unit operates at capacity either at the highest price interval tλ^maxsubscriptsubscript𝑡^𝜆max{t_{\hat{\lambda}}}_{\mathrm{max}}italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT or the lowest price interval tλ^minsubscriptsubscript𝑡^𝜆min{t_{\hat{\lambda}}}_{\mathrm{min}}italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT: ptλ^max=P¯subscriptsubscriptsuperscript𝑝subscript𝑡^𝜆max¯𝑃{p^{*}_{t_{\hat{\lambda}}}}_{\mathrm{max}}=\bar{P}italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG or btλ^min=P¯subscriptsubscriptsuperscript𝑏subscript𝑡^𝜆min¯𝑃{b^{*}_{t_{\hat{\lambda}}}}_{\mathrm{min}}=\bar{P}italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG.

Proof:

We first prove 1) that the prices throughout the scheduling period 𝒯𝒯\mathcal{T}caligraphic_T are not identical. This can be done by assuming that all prices are identical: λ^=λ^1=λ^2==λ^tT^𝜆subscript^𝜆1subscript^𝜆2subscript^𝜆subscript𝑡𝑇\hat{\lambda}=\hat{\lambda}_{1}=\hat{\lambda}_{2}=\ldots=\hat{\lambda}_{t_{T}}over^ start_ARG italic_λ end_ARG = over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = … = over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT end_POSTSUBSCRIPT. Thus the profit is calculated as π=t𝒯λ^t(ptbt)=t𝒯λ^(ptbt)𝜋subscript𝑡𝒯subscript^𝜆𝑡subscriptsuperscript𝑝𝑡subscriptsuperscript𝑏𝑡subscript𝑡𝒯^𝜆subscriptsuperscript𝑝𝑡subscriptsuperscript𝑏𝑡\pi=\sum_{t\in\mathcal{T}}\hat{\lambda}_{t}(p^{*}_{t}-b^{*}_{t})=\sum_{t\in% \mathcal{T}}\hat{\lambda}(p^{*}_{t}-b^{*}_{t})italic_π = ∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) = ∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG ( italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ). Given the storage unit charge and discharge reset constraint (7b), we have t𝒯pt=t𝒯btη2subscript𝑡𝒯subscriptsuperscript𝑝𝑡subscript𝑡𝒯subscriptsuperscript𝑏𝑡superscript𝜂2\sum_{t\in\mathcal{T}}p^{*}_{t}=\sum_{t\in\mathcal{T}}b^{*}_{t}\eta^{2}∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = ∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, which yields a negative profit π𝜋\piitalic_π, violating the non-idle condition. Thus, there must exist distinct minimum and maximum prices λ^minsubscript^𝜆min\hat{\lambda}_{\mathrm{min}}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT and λ^maxsubscript^𝜆max\hat{\lambda}_{\mathrm{max}}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT.

Next, we prove the storage unit operates at capacity either at the highest or the lowest price interval as 2). The counter conditions include i) both of these two intervals are at partial capacity, ii) one of these two intervals is at partial capacity and another is idle, or iii) both of these two intervals are idle. For case i), given ptλ^max,btλ^min(0,P¯)subscriptsubscriptsuperscript𝑝subscript𝑡^𝜆maxsubscriptsubscriptsuperscript𝑏subscript𝑡^𝜆min0¯𝑃{p^{*}_{t_{\hat{\lambda}}}}_{\mathrm{max}},~{}{b^{*}_{t_{\hat{\lambda}}}}_{% \mathrm{min}}~{}\in(0,\bar{P})italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT , italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT ∈ ( 0 , over¯ start_ARG italic_P end_ARG ), we have δtλ^max,δ+tλ^max,βtλ^min,β+tλ^min=0subscriptsubscriptsuperscriptsuperscript𝛿subscript𝑡^𝜆maxsubscriptsubscriptsuperscriptsuperscript𝛿subscript𝑡^𝜆maxsubscriptsubscriptsuperscriptsuperscript𝛽subscript𝑡^𝜆minsubscriptsubscriptsuperscriptsuperscript𝛽subscript𝑡^𝜆min0{{\delta^{-}}^{*}_{t_{\hat{\lambda}}}}_{\mathrm{max}},~{}{{\delta^{+}}^{*}_{t_% {\hat{\lambda}}}}_{\mathrm{max}},~{}{{\beta^{-}}^{*}_{t_{\hat{\lambda}}}}_{% \mathrm{min}},~{}{{\beta^{+}}^{*}_{t_{\hat{\lambda}}}}_{\mathrm{min}}=0italic_δ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT , italic_δ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT , italic_β start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT , italic_β start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT = 0 based on (9f) and (9g). Following Proposition 1, we assume the units operates at capacity at interval t0subscript𝑡0t_{0}italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. Thus we have for these intervals:

ifpt0=P¯,then{λ^max+θη=0,λ^minθη=0,λ^t0+θηδt0+=0,\displaystyle\mathrm{if}~{}{p^{*}_{t_{0}}}=\bar{P},~{}\mathrm{then}\left\{% \begin{aligned} &\hat{\lambda}_{\mathrm{max}}+\frac{\theta}{\eta}=0,\\ &-\hat{\lambda}_{\mathrm{min}}-\theta\eta=0,\\ &\hat{\lambda}_{t_{0}}+\frac{\theta}{\eta}-\delta_{t_{0}}^{+}=0,\end{aligned}\right.roman_if italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG , roman_then { start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT - italic_θ italic_η = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG - italic_δ start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = 0 , end_CELL end_ROW (11)
oror\displaystyle\mathrm{or}roman_or
ifbt0=P¯,then{λ^max+θη=0,λ^minθη=0,λ^t0θηβt0+=0.\displaystyle\mathrm{if}~{}{b^{*}_{t_{0}}}=\bar{P},~{}\mathrm{then}\left\{% \begin{aligned} &\hat{\lambda}_{\mathrm{max}}+\frac{\theta}{\eta}=0,\\ &-\hat{\lambda}_{\mathrm{min}}-\theta\eta=0,\\ &-\hat{\lambda}_{t_{0}}-\theta\eta-\beta_{t_{0}}^{+}=0.\end{aligned}\right.roman_if italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG , roman_then { start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT - italic_θ italic_η = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - italic_θ italic_η - italic_β start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = 0 . end_CELL end_ROW (12)

Solutions to (11) need to satisfy λ^max=λ^t0δt0+subscript^𝜆maxsubscript^𝜆subscript𝑡0superscriptsubscript𝛿subscript𝑡0\hat{\lambda}_{\mathrm{max}}=\hat{\lambda}_{t_{0}}-\delta_{t_{0}}^{+}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT = over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - italic_δ start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT. Given non-negative δt0+superscriptsubscript𝛿subscript𝑡0\delta_{t_{0}}^{+}italic_δ start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT, such solutions don’t exist; similar for solutions to (12). Therefore, case i) doesn’t exist. Similar proof can be provided for case ii) and case iii). Therefore, the storage unit must be at capacity either at tλ^maxsubscriptsubscript𝑡^𝜆max{t_{\hat{\lambda}}}_{\mathrm{max}}italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT or tλ^minsubscriptsubscript𝑡^𝜆min{t_{\hat{\lambda}}}_{\mathrm{min}}italic_t start_POSTSUBSCRIPT over^ start_ARG italic_λ end_ARG end_POSTSUBSCRIPT start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT. ∎

Remark 3

Propositions 1 and 2 suggest that during the scheduling period, the storage unit will operate at its capacity for at least one interval, either at the highest or the lowest price point. This insight is readily applicable to scenarios where the SoC of the storage unit is constrained throughout the period. In situations where the storage unit’s maximum output is limited by its capacity, it can still be asserted, without loss of generality, that the storage unit will maximize its energy utilization during the period of highest or lowest prices. This principle aligns with the conjecture made in Section III.A of [11], which stated that “for energy limited unit, efficient pricing would fully utilize the energy of the unit in the highest price hours over the period of the limitation”. Our findings validate this conjecture and also extend it by exploring its implications for the bidding strategies of energy storage units.

Proposition 3

For a strategic price taker making bidding decisions based on model (7), given strictly heterogeneous prices λ^1λ^2λ^Tsubscript^𝜆1subscript^𝜆2subscript^𝜆𝑇\hat{\lambda}_{1}\neq\hat{\lambda}_{2}\ldots\neq\hat{\lambda}_{T}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≠ over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT … ≠ over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT, the bidding decisions {pt,bt}subscriptsuperscript𝑝𝑡subscriptsuperscript𝑏𝑡\{p^{*}_{t},~{}b^{*}_{t}\}{ italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT } throughout the period 𝒯𝒯\mathcal{T}caligraphic_T should include one and only one interval t𝑡titalic_t at partial capacity, 0<pt<P¯0subscriptsuperscript𝑝𝑡¯𝑃0<p^{*}_{t}<\bar{P}0 < italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG or 0<bt<P¯0subscriptsuperscript𝑏𝑡¯𝑃0<b^{*}_{t}<\bar{P}0 < italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG, except for the idle scenarios.

Proof:

To prove that there exists one and only one partial interval throughout the bidding period, we first prove there is at least one interval with partial capacity and then prove at most one interval with partial capacity.

First, we prove at least one interval operates at partial capacity. For non-idle intervals 𝒯𝒯superscript𝒯𝒯\mathcal{T}^{\prime}\subseteq\mathcal{T}caligraphic_T start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ⊆ caligraphic_T, suppose the unit operates at full capacity pt=P¯subscriptsuperscript𝑝𝑡¯𝑃p^{*}_{t}=\bar{P}italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG or bt=P¯,t𝒯formulae-sequencesubscriptsuperscript𝑏𝑡¯𝑃for-all𝑡superscript𝒯b^{*}_{t}=\bar{P},\,\forall t\in\mathcal{T}^{\prime}italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG , ∀ italic_t ∈ caligraphic_T start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT. Obviously, it violates the charging and discharging balance t𝒯ptηbtη=0subscript𝑡superscript𝒯superscriptsubscript𝑝𝑡𝜂superscriptsubscript𝑏𝑡𝜂0\sum_{t\in\mathcal{T}^{\prime}}\frac{p_{t}^{*}}{\eta}-b_{t}^{*}\eta=0∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT divide start_ARG italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT end_ARG start_ARG italic_η end_ARG - italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT italic_η = 0.

Second, we prove at most one interval operates at partial capacity. Suppose there are more than one interval that operate at partial capacity. For partial intervals, we have δt+,βt+=0superscriptsuperscriptsubscript𝛿𝑡superscriptsuperscriptsubscript𝛽𝑡0{\delta_{t}^{+}}^{*},~{}{\beta_{t}^{+}}^{*}=0italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT , italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = 0 and δt=0superscriptsuperscriptsubscript𝛿𝑡0{\delta_{t}^{-}}^{*}=0italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = 0 or βt=0,t{t|0<pt<P¯|| 0<bt<P¯,t𝒯}formulae-sequencesuperscriptsuperscriptsubscript𝛽𝑡0for-all𝑡conditional-set𝑡formulae-sequence0brasubscriptsuperscript𝑝𝑡bra¯𝑃 0subscriptsuperscript𝑏𝑡¯𝑃𝑡𝒯{\beta_{t}^{-}}^{*}=0,\,\forall t\in\{t|0<p^{*}_{t}<\bar{P}\ ||\ 0<b^{*}_{t}<% \bar{P},t\in\mathcal{T}\}italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = 0 , ∀ italic_t ∈ { italic_t | 0 < italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG | | 0 < italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG , italic_t ∈ caligraphic_T } based on (9f) and (9g), thus

{λ^t+θη+δt=0,λ^tθη=0,or{λ^t+θη=0,λ^tθη+βt=0,\left\{\begin{aligned} &\hat{\lambda}_{t}+\frac{\theta}{\eta}+\delta_{t}^{-}=0% ,\\ &-\hat{\lambda}_{t}-\theta\eta=0,\end{aligned}\right.\quad\text{or}\quad\left% \{\begin{aligned} &\hat{\lambda}_{t}+\frac{\theta}{\eta}=0,\\ &-\hat{\lambda}_{t}-\theta\eta+\beta_{t}^{-}=0,\end{aligned}\right.{ start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG + italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_θ italic_η = 0 , end_CELL end_ROW or { start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_θ italic_η + italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = 0 , end_CELL end_ROW

leading to

{λ^t1+θη=0,λ^t2θη=0,or{λ^t1+θη=0,λ^t2+θη=0,or{λ^t1θη=0,λ^t2θη=0.\left\{\begin{aligned} &\hat{\lambda}_{t_{1}}+\frac{\theta}{\eta}=0,\\ &-\hat{\lambda}_{t_{2}}-\theta\eta=0,\end{aligned}\right.\ \text{or}\ \left\{% \begin{aligned} &\hat{\lambda}_{t_{1}}+\frac{\theta}{\eta}=0,\\ &\hat{\lambda}_{t_{2}}+\frac{\theta}{\eta}=0,\\ \end{aligned}\right.\ \text{or}\ \left\{\begin{aligned} &-\hat{\lambda}_{t_{1}% }-\theta\eta=0,\\ &-\hat{\lambda}_{t_{2}}-\theta\eta=0.\end{aligned}\right.{ start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - italic_θ italic_η = 0 , end_CELL end_ROW or { start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG = 0 , end_CELL end_ROW or { start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - italic_θ italic_η = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - italic_θ italic_η = 0 . end_CELL end_ROW (13)

All pairs in (13) would have led to idle intervals as discussed in the proof for Proposition 1, or require equal prices, which violate either the non-idle or the strictly heterogeneous price assumption. ∎

In competitive markets, i.e., market participants bid as price takers according to model (7), we have the following price-decision relationship:

Proposition 4

Given a series of prices λ^tsubscript^𝜆𝑡\hat{\lambda}_{t}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT throughout the period 𝒯𝒯\mathcal{T}caligraphic_T, a strategic price taker makes bidding decisions {pt,bt}subscriptsuperscript𝑝𝑡subscriptsuperscript𝑏𝑡\{p^{*}_{t},~{}b^{*}_{t}\}{ italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT } based on model (7). Denote the set of discharge withholding intervals {u𝒯|𝟙{0<pu<P¯}=1}conditional-set𝑢𝒯subscript10subscript𝑝𝑢¯𝑃1\{u\in\mathcal{T}|\mathds{1}_{\{0<p_{u}<\bar{P}\}}=1\}{ italic_u ∈ caligraphic_T | blackboard_1 start_POSTSUBSCRIPT { 0 < italic_p start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG } end_POSTSUBSCRIPT = 1 } and charge withholding intervals {v𝒯|𝟙{0<bv<P¯}=1}conditional-set𝑣𝒯subscript10subscript𝑏𝑣¯𝑃1\{v\in\mathcal{T}|\mathds{1}_{\{0<b_{v}<\bar{P}\}}=1\}{ italic_v ∈ caligraphic_T | blackboard_1 start_POSTSUBSCRIPT { 0 < italic_b start_POSTSUBSCRIPT italic_v end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG } end_POSTSUBSCRIPT = 1 }, then the bidding decisions satisfy:

  1. 1.

    if the unit discharges at capacity during interval x𝑥xitalic_x, i.e., px=P¯subscriptsuperscript𝑝𝑥¯𝑃p^{*}_{x}=\bar{P}italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG, then λ^x>λ^usubscript^𝜆𝑥subscript^𝜆𝑢\hat{\lambda}_{x}>\hat{\lambda}_{u}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT > over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT and λ^x>λ^vη2subscript^𝜆𝑥subscript^𝜆𝑣superscript𝜂2\hat{\lambda}_{x}>\frac{\hat{\lambda}_{v}}{\eta^{2}}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT > divide start_ARG over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_v end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ,

  2. 2.

    if the unit charges at capacity during interval y𝑦yitalic_y, i.e., by=P¯subscriptsuperscript𝑏𝑦¯𝑃b^{*}_{y}=\bar{P}italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG, then λ^u>λ^yη2subscript^𝜆𝑢subscript^𝜆𝑦superscript𝜂2\hat{\lambda}_{u}>\frac{\hat{\lambda}_{y}}{\eta^{2}}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT > divide start_ARG over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG and λ^v>λ^ysubscript^𝜆𝑣subscript^𝜆𝑦\hat{\lambda}_{v}>\hat{\lambda}_{y}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_v end_POSTSUBSCRIPT > over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT,

  3. 3.

    if the unit keeps idle during interval z𝑧zitalic_z, i.e., pz=bz=0subscriptsuperscript𝑝𝑧subscriptsuperscript𝑏𝑧0p^{*}_{z}=b^{*}_{z}=0italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0, then λ^zη2>λ^u>λ^zsubscript^𝜆𝑧superscript𝜂2subscript^𝜆𝑢subscript^𝜆𝑧\frac{\hat{\lambda}_{z}}{\eta^{2}}>\hat{\lambda}_{u}>\hat{\lambda}_{z}divide start_ARG over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG > over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT > over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT and λ^z>λ^v>λ^zη2subscript^𝜆𝑧subscript^𝜆𝑣subscript^𝜆𝑧superscript𝜂2\hat{\lambda}_{z}>\hat{\lambda}_{v}>\hat{\lambda}_{z}\eta^{2}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT > over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_v end_POSTSUBSCRIPT > over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT.

Proof:

Following the KKT conditions in (9) and similar to the previous discussion, we have:

for a discharge withholding interval u𝑢uitalic_u
{λ^u+θη=0,λ^uθη+βu=0,\left\{\begin{aligned} &\hat{\lambda}_{u}+\frac{\theta}{\eta}=0,\\ &-\hat{\lambda}_{u}-\theta\eta+\beta_{u}^{-}=0,\end{aligned}\right.{ start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT - italic_θ italic_η + italic_β start_POSTSUBSCRIPT italic_u end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = 0 , end_CELL end_ROW (14)
for a charge withholding interval v𝑣vitalic_v
{λ^v+θη+δv=0,λ^vθη=0,\left\{\begin{aligned} &\hat{\lambda}_{v}+\frac{\theta}{\eta}+\delta_{v}^{-}=0% ,\\ &-\hat{\lambda}_{v}-\theta\eta=0,\end{aligned}\right.{ start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_v end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG + italic_δ start_POSTSUBSCRIPT italic_v end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_v end_POSTSUBSCRIPT - italic_θ italic_η = 0 , end_CELL end_ROW
for a discharge at-capacity interval x𝑥xitalic_x
{λ^x+θηδx+=0,λ^xθη+βx=0,\left\{\begin{aligned} &\hat{\lambda}_{x}+\frac{\theta}{\eta}-\delta^{+}_{x}=0% ,\\ &-\hat{\lambda}_{x}-\theta\eta+\beta_{x}^{-}=0,\end{aligned}\right.{ start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG - italic_δ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT - italic_θ italic_η + italic_β start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = 0 , end_CELL end_ROW
for a charge at-capacity interval y𝑦yitalic_y
{λ^y+θη+δy=0,λ^yθηβy+=0,\left\{\begin{aligned} &\hat{\lambda}_{y}+\frac{\theta}{\eta}+\delta^{-}_{y}=0% ,\\ &-\hat{\lambda}_{y}-\theta\eta-\beta_{y}^{+}=0,\end{aligned}\right.{ start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG + italic_δ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT - italic_θ italic_η - italic_β start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = 0 , end_CELL end_ROW
and for an idle interval z𝑧zitalic_z
{λ^z+θη+δz=0,λ^zθη+βz=0.\left\{\begin{aligned} &\hat{\lambda}_{z}+\frac{\theta}{\eta}+\delta^{-}_{z}=0% ,\\ &-\hat{\lambda}_{z}-\theta\eta+\beta_{z}^{-}=0.\end{aligned}\right.{ start_ROW start_CELL end_CELL start_CELL over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG + italic_δ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = 0 , end_CELL end_ROW start_ROW start_CELL end_CELL start_CELL - over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT - italic_θ italic_η + italic_β start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = 0 . end_CELL end_ROW

Given non-negative dual variables δ(),δ()+,β(),β()+superscriptsubscript𝛿superscriptsubscript𝛿superscriptsubscript𝛽superscriptsubscript𝛽\delta_{(\cdot)}^{-},\delta_{(\cdot)}^{+},\beta_{(\cdot)}^{-},\beta_{(\cdot)}^% {+}italic_δ start_POSTSUBSCRIPT ( ⋅ ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_δ start_POSTSUBSCRIPT ( ⋅ ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT , italic_β start_POSTSUBSCRIPT ( ⋅ ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT , italic_β start_POSTSUBSCRIPT ( ⋅ ) end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT, where ()(\cdot)( ⋅ ) represents intervals {u,v,x,y,z}𝑢𝑣𝑥𝑦𝑧\{u,v,x,y,z\}{ italic_u , italic_v , italic_x , italic_y , italic_z }, the proof is straightforward. ∎

TABLE I: Storage Unit Control Policy as Price Taker in Two-Interval Bidding
Scenario Interval 1111 Interval 2222
p1superscriptsubscript𝑝1p_{1}^{*}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT b1superscriptsubscript𝑏1b_{1}^{*}italic_b start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT p2superscriptsubscript𝑝2p_{2}^{*}italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT b2superscriptsubscript𝑏2b_{2}^{*}italic_b start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT
λ^1>λ^2η2subscript^𝜆1subscript^𝜆2superscript𝜂2\hat{\lambda}_{1}>\frac{\hat{\lambda}_{2}}{\eta^{2}}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > divide start_ARG over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG P¯η2¯𝑃superscript𝜂2\bar{P}\eta^{2}over¯ start_ARG italic_P end_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT 0 0 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
λ^2η2λ^1λ^2η2subscript^𝜆2superscript𝜂2subscript^𝜆1subscript^𝜆2superscript𝜂2\hat{\lambda}_{2}\eta^{2}\leq\hat{\lambda}_{1}\leq\frac{\hat{\lambda}_{2}}{% \eta^{2}}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≤ over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≤ divide start_ARG over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG 0 0 0 0
λ^1<λ^2η2subscript^𝜆1subscript^𝜆2superscript𝜂2\hat{\lambda}_{1}<\hat{\lambda}_{2}\eta^{2}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT 0 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯η2¯𝑃superscript𝜂2\bar{P}\eta^{2}over¯ start_ARG italic_P end_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT 0
TABLE II: Storage Unit Control Policy as Price Maker in Two-Interval Bidding
Scenario Interval 1111 Interval 2222
p1superscriptsubscript𝑝1p_{1}^{*}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT b1superscriptsubscript𝑏1b_{1}^{*}italic_b start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT p2superscriptsubscript𝑝2p_{2}^{*}italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT b2superscriptsubscript𝑏2b_{2}^{*}italic_b start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT
λ¯1>λ¯2η2subscript¯𝜆1subscript¯𝜆2superscript𝜂2\bar{\lambda}_{1}>\frac{\bar{\lambda}_{2}}{\eta^{2}}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG λ¯12α1P¯η2λ¯2+2α2P¯η2subscript¯𝜆12subscript𝛼1¯𝑃superscript𝜂2subscript¯𝜆22subscript𝛼2¯𝑃superscript𝜂2\bar{\lambda}_{1}-2\alpha_{1}\bar{P}\eta^{2}\geq\frac{\bar{\lambda}_{2}+2% \alpha_{2}\bar{P}}{\eta^{2}}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≥ divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 2 italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG P¯η2¯𝑃superscript𝜂2\bar{P}\eta^{2}over¯ start_ARG italic_P end_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT 0 0 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
λ¯12α1P¯η2<λ¯2+2α2P¯η2subscript¯𝜆12subscript𝛼1¯𝑃superscript𝜂2subscript¯𝜆22subscript𝛼2¯𝑃superscript𝜂2\bar{\lambda}_{1}-2\alpha_{1}\bar{P}\eta^{2}<\frac{\bar{\lambda}_{2}+2\alpha_{% 2}\bar{P}}{\eta^{2}}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT < divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 2 italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG λ¯1λ¯2η22(α1+α2η4)subscript¯𝜆1subscript¯𝜆2superscript𝜂22subscript𝛼1subscript𝛼2superscript𝜂4\frac{\bar{\lambda}_{1}-\frac{\bar{\lambda}_{2}}{\eta^{2}}}{2(\alpha_{1}+\frac% {\alpha_{2}}{\eta^{4}})}divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG 2 ( italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + divide start_ARG italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) end_ARG 0 0 λ¯1λ¯2η22(α1+α2η4)η2subscript¯𝜆1subscript¯𝜆2superscript𝜂22subscript𝛼1subscript𝛼2superscript𝜂4superscript𝜂2\frac{\bar{\lambda}_{1}-\frac{\bar{\lambda}_{2}}{\eta^{2}}}{2(\alpha_{1}+\frac% {\alpha_{2}}{\eta^{4}})\eta^{2}}divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG start_ARG 2 ( italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + divide start_ARG italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT end_ARG ) italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG
λ¯2η2λ¯1λ¯2η2subscript¯𝜆2superscript𝜂2subscript¯𝜆1subscript¯𝜆2superscript𝜂2\bar{\lambda}_{2}\eta^{2}\leq\bar{\lambda}_{1}\leq\frac{\bar{\lambda}_{2}}{% \eta^{2}}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≤ over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≤ divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG 0 0 0 0
λ¯1<λ¯2η2subscript¯𝜆1subscript¯𝜆2superscript𝜂2\bar{\lambda}_{1}<\bar{\lambda}_{2}\eta^{2}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT λ¯1+2α1P¯η2>λ¯22α2P¯η2subscript¯𝜆12subscript𝛼1¯𝑃superscript𝜂2subscript¯𝜆22subscript𝛼2¯𝑃superscript𝜂2\frac{\bar{\lambda}_{1}+2\alpha_{1}\bar{P}}{\eta^{2}}>\bar{\lambda}_{2}-2% \alpha_{2}\bar{P}\eta^{2}divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG > over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT 0 λ¯1+λ¯2η22(α1+α2η4)subscript¯𝜆1subscript¯𝜆2superscript𝜂22subscript𝛼1subscript𝛼2superscript𝜂4\frac{-\bar{\lambda}_{1}+\bar{\lambda}_{2}\eta^{2}}{2(\alpha_{1}+\alpha_{2}% \eta^{4})}divide start_ARG - over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ) end_ARG (λ¯1+λ¯2η2)η22(α1+α2η4)subscript¯𝜆1subscript¯𝜆2superscript𝜂2superscript𝜂22subscript𝛼1subscript𝛼2superscript𝜂4\frac{(-\bar{\lambda}_{1}+\bar{\lambda}_{2}\eta^{2})\eta^{2}}{2(\alpha_{1}+% \alpha_{2}\eta^{4})}divide start_ARG ( - over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 ( italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT ) end_ARG 0
λ¯1+2α1P¯η2λ¯22α2P¯η2subscript¯𝜆12subscript𝛼1¯𝑃superscript𝜂2subscript¯𝜆22subscript𝛼2¯𝑃superscript𝜂2\frac{\bar{\lambda}_{1}+2\alpha_{1}\bar{P}}{\eta^{2}}\leq\bar{\lambda}_{2}-2% \alpha_{2}\bar{P}\eta^{2}divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 2 italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ≤ over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT 0 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯η2¯𝑃superscript𝜂2\bar{P}\eta^{2}over¯ start_ARG italic_P end_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT 0

We illustrate the bidding strategies of price takers based on model (7) through the two-interval bidding scenario, with the control policy detailed in Table I. The decision-making process is primarily influenced by the price difference between the two intervals. Considering the energy loss during the charging/discharging cycle, the storage unit operation is justified only if the losses are offset by the profits derived from arbitraging the price difference. This arbitrage strategy indicates that the storage unit will fully utilize its capacity for profit as long as it remains economically viable, without engaging in unjustified withholding, i.e., no partial operation below P¯η2¯𝑃superscript𝜂2\bar{P}\eta^{2}over¯ start_ARG italic_P end_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT.

III-B Withholding as Price Maker

We now explore the scenario in which the storage unit operates as a price maker, anticipating the effects of its bidding behavior on the market clearing price, using the following bidding model:

maximizept,btsubscript𝑝𝑡subscript𝑏𝑡maximize\displaystyle\underset{p_{t},b_{t}}{\text{maximize}}start_UNDERACCENT italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_UNDERACCENT start_ARG maximize end_ARG t𝒯(λ¯tαt(ptbt))(ptbt)subscript𝑡𝒯subscript¯𝜆𝑡subscript𝛼𝑡subscript𝑝𝑡subscript𝑏𝑡subscript𝑝𝑡subscript𝑏𝑡\displaystyle\quad\sum_{t\in\mathcal{T}}(\bar{\lambda}_{t}-\alpha_{t}(p_{t}-b_% {t}))(p_{t}-b_{t})∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT ( over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ( italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) ) ( italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ) (15)
s.t.formulae-sequencest\displaystyle\mathrm{s.t.}roman_s . roman_t . (7b)(7d).italic-(7bitalic-)italic-(7ditalic-)\displaystyle\quad\eqref{eq:power_bal}-\eqref{eq:b_bar}.italic_( italic_) - italic_( italic_) .

Incorporating the impacts of market power on market clearing prices into the price maker’s bidding model introduces quadratic terms into the objective function (15), in contrast to the linear terms found in the price-taker’s model as (7a). However, both problems are still convex.

Thus, the corresponding KKT conditions are (t𝒯for-all𝑡𝒯\forall t\in\mathcal{T}∀ italic_t ∈ caligraphic_T):

Lpt=λ¯t2αtpt+2αtbt+θη+δtδt+=0𝐿subscript𝑝𝑡subscript¯𝜆𝑡2subscript𝛼𝑡subscript𝑝𝑡2subscript𝛼𝑡subscript𝑏𝑡𝜃𝜂superscriptsubscript𝛿𝑡superscriptsubscript𝛿𝑡0\displaystyle\frac{\partial L}{\partial p_{t}}=\bar{\lambda}_{t}-2\alpha_{t}p_% {t}+2\alpha_{t}b_{t}+\frac{\theta}{\eta}+\delta_{t}^{-}-\delta_{t}^{+}=0divide start_ARG ∂ italic_L end_ARG start_ARG ∂ italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG = over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + 2 italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG + italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT - italic_δ start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = 0 (16a)
Lbt=λ¯t+2αtpt2αtbtθη+βtβt+=0𝐿subscript𝑏𝑡subscript¯𝜆𝑡2subscript𝛼𝑡subscript𝑝𝑡2subscript𝛼𝑡subscript𝑏𝑡𝜃𝜂superscriptsubscript𝛽𝑡superscriptsubscript𝛽𝑡0\displaystyle\frac{\partial L}{\partial b_{t}}=-\bar{\lambda}_{t}+2\alpha_{t}p% _{t}-2\alpha_{t}b_{t}-\theta\eta+\beta_{t}^{-}-\beta_{t}^{+}=0divide start_ARG ∂ italic_L end_ARG start_ARG ∂ italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT end_ARG = - over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT + 2 italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT - italic_θ italic_η + italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT - italic_β start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = 0 (16b)
(9c)(9e)italic-(9citalic-)italic-(9eitalic-)\displaystyle\eqref{eq:kkt_bal}-\eqref{eq:kkt_dual}italic_( italic_) - italic_( italic_)

where Eqs.(16a) and (16b) correspond to stationarity conditions, and the rest are the same as that of the optimality conditions in (9).

Proposition 5

For a strategic price maker making bidding decisions based on model (15), throughout the period 𝒯𝒯\mathcal{T}caligraphic_T

  1. 1.

    the bidding decisions {pt,bt,t𝒯}subscriptsuperscript𝑝𝑡subscriptsuperscript𝑏𝑡for-all𝑡𝒯\{p^{*}_{t},~{}b^{*}_{t},\,\forall t\in\mathcal{T}\}{ italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , ∀ italic_t ∈ caligraphic_T } should include at least one interval t𝑡titalic_t at partial capacity, 0<pt<P¯0subscriptsuperscript𝑝𝑡¯𝑃0<p^{*}_{t}<\bar{P}0 < italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG or 0<bt<P¯0subscriptsuperscript𝑏𝑡¯𝑃0<b^{*}_{t}<\bar{P}0 < italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG, except for the idle scenarios,

  2. 2.

    there might exist multiple full or partial intervals.

Proof:

1) follows similar arguments to Proposition 3; 2) is apparent following the proof for Proposition 3 based on the KKT conditions (16). ∎

The two-interval bidding scenario for price makers based on model (15) is summarized in Table II. Compared to the strategies of price takers as in Table I, the control decisions for each non-idle scenario further splits into two sub-scenarios, either λ¯1>λ¯2η2subscript¯𝜆1subscript¯𝜆2superscript𝜂2\bar{\lambda}_{1}>\frac{\bar{\lambda}_{2}}{\eta^{2}}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG or λ¯1<λ¯2η2subscript¯𝜆1subscript¯𝜆2superscript𝜂2\bar{\lambda}_{1}<\bar{\lambda}_{2}\eta^{2}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Such decisions are established accounting for the economic losses associated with both charging/discharging energy loss and their potential impact on market prices, including both immediate and future contrast effects. Take the scenario where λ¯1>λ¯2η2subscript¯𝜆1subscript¯𝜆2superscript𝜂2\bar{\lambda}_{1}>\frac{\bar{\lambda}_{2}}{\eta^{2}}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG for example. In addition to comparing prices based on charging and discharging efficiencies, i.e., λ¯1>λ¯2η2subscript¯𝜆1subscript¯𝜆2superscript𝜂2\bar{\lambda}_{1}>\frac{\bar{\lambda}_{2}}{\eta^{2}}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG, price makers also consider the marginal revenue from operational intervals. If the projected profit be substantial, i.e., with high marginal revenue as λ¯12α1P¯η2λ¯2+2α2P¯η2subscript¯𝜆12subscript𝛼1¯𝑃superscript𝜂2subscript¯𝜆22subscript𝛼2¯𝑃superscript𝜂2\bar{\lambda}_{1}-2\alpha_{1}\bar{P}\eta^{2}\geq\frac{\bar{\lambda}_{2}+2% \alpha_{2}\bar{P}}{\eta^{2}}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≥ divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 2 italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , the strategy aligns with that of price takers as operating at full capacity, suggesting that a larger price difference leads to greater profitability, sufficient even for arbitrage by price takers. Moreover, price makers, facing potentially lower marginal revenues, i.e., λ¯12α1P¯η2<λ¯2+2α2P¯η2subscript¯𝜆12subscript𝛼1¯𝑃superscript𝜂2subscript¯𝜆22subscript𝛼2¯𝑃superscript𝜂2\bar{\lambda}_{1}-2\alpha_{1}\bar{P}\eta^{2}<\frac{\bar{\lambda}_{2}+2\alpha_{% 2}\bar{P}}{\eta^{2}}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT < divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 2 italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG, still possess the strategic flexibility to optimize control decisions for maximizing profit, exploiting their influence on market outcome.

III-C Three-Interval Bidding Scenario

We further examine the three-interval bidding scenario to help illustrate the bidding pattern of price makers. Compared to the two-interval bidding discussed earlier, the three-interval scenario demonstrates how energy storage units strategically allocate their limited resources into a portfolio that maximizes total profit based on a series of prices over the scheduling period. The three-interval bidding establishes a foundation for analyzing longer-term bidding scenarios.

Define the number of full capacity intervals throughout the scheduling period 𝒯𝒯\mathcal{T}caligraphic_T: o=t𝒯𝟙{pt=P¯}+t𝒯𝟙{bt=P¯}𝑜subscript𝑡𝒯subscript1subscript𝑝𝑡¯𝑃subscript𝑡𝒯subscript1subscript𝑏𝑡¯𝑃o=\sum_{t\in{\mathcal{T}}}\mathds{1}_{\{p_{t}=\bar{P}\}}+\sum_{t\in{\mathcal{T% }}}\mathds{1}_{\{b_{t}=\bar{P}\}}italic_o = ∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT blackboard_1 start_POSTSUBSCRIPT { italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG } end_POSTSUBSCRIPT + ∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT blackboard_1 start_POSTSUBSCRIPT { italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = over¯ start_ARG italic_P end_ARG } end_POSTSUBSCRIPT, and correspondingly for partial intervals: s=t𝒯𝟙{0<pt<P¯}+t𝒯𝟙{0<bt<P¯}𝑠subscript𝑡𝒯subscript10subscript𝑝𝑡¯𝑃subscript𝑡𝒯subscript10subscript𝑏𝑡¯𝑃s=\sum_{t\in{\mathcal{T}}}\mathds{1}_{\{0<p_{t}<\bar{P}\}}+\sum_{t\in{\mathcal% {T}}}\mathds{1}_{\{0<b_{t}<\bar{P}\}}italic_s = ∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT blackboard_1 start_POSTSUBSCRIPT { 0 < italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG } end_POSTSUBSCRIPT + ∑ start_POSTSUBSCRIPT italic_t ∈ caligraphic_T end_POSTSUBSCRIPT blackboard_1 start_POSTSUBSCRIPT { 0 < italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG } end_POSTSUBSCRIPT.

Corollary 1

For a strategic price maker making bidding decisions based on model (15), consider the bidding case consisting three intervals, then if there exist two intervals j1,j2subscript𝑗1subscript𝑗2j_{1},~{}j_{2}italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT operating with partial capacity, s=2𝑠2s=2italic_s = 2, and one full t0subscript𝑡0t_{0}italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, o=1𝑜1o=1italic_o = 1 and λ¯j2>λ¯j1,λ¯j2>λ¯t0formulae-sequencesubscript¯𝜆subscript𝑗2subscript¯𝜆subscript𝑗1subscript¯𝜆subscript𝑗2subscript¯𝜆subscript𝑡0\bar{\lambda}_{j_{2}}>\bar{\lambda}_{j_{1}},\bar{\lambda}_{j_{2}}>\bar{\lambda% }_{t_{0}}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT > over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT , over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT > over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT, a sufficient condition for this scenario is λ¯j1+2αj1P¯η2>λ¯j22αj2P¯subscript¯𝜆subscript𝑗12subscript𝛼subscript𝑗1¯𝑃superscript𝜂2subscript¯𝜆subscript𝑗22subscript𝛼subscript𝑗2¯𝑃\frac{\bar{\lambda}_{j_{1}}+2\alpha_{j_{1}}\bar{P}}{\eta^{2}}>\bar{\lambda}_{j% _{2}}-2\alpha_{j_{2}}\bar{P}divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + 2 italic_α start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG > over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG.

Proof:

Note that charging and discharging intervals appear in pairs, then based on KKT conditions (9f), (9g), (16a) and (16b) we have

λ¯t02αt0P¯θηβt0+=0(full),subscript¯𝜆subscript𝑡02subscript𝛼subscript𝑡0¯𝑃𝜃𝜂superscriptsubscript𝛽subscript𝑡00full\displaystyle-\bar{\lambda}_{t_{0}}-2\alpha_{t_{0}}\bar{P}-\theta\eta-\beta_{t% _{0}}^{+}=0\ (\mathrm{full}),- over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG - italic_θ italic_η - italic_β start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = 0 ( roman_full ) , (17a)
λ¯j12αj1bj1θη=0(partial),subscript¯𝜆subscript𝑗12subscript𝛼subscript𝑗1subscript𝑏subscript𝑗1𝜃𝜂0partial\displaystyle-\bar{\lambda}_{j_{1}}-2\alpha_{j_{1}}b_{j_{1}}-\theta\eta=0\ (% \mathrm{partial}),- over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_b start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - italic_θ italic_η = 0 ( roman_partial ) , (17b)
λ¯j22αj2pj2+θη=0(partial).subscript¯𝜆subscript𝑗22subscript𝛼subscript𝑗2subscript𝑝subscript𝑗2𝜃𝜂0partial\displaystyle\bar{\lambda}_{j_{2}}-2\alpha_{j_{2}}p_{j_{2}}+\frac{\theta}{\eta% }=0\ (\mathrm{partial}).over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + divide start_ARG italic_θ end_ARG start_ARG italic_η end_ARG = 0 ( roman_partial ) . (17c)

Combining (17b) and (17c), we have

λ¯j1+2αj1bj1η2=λ¯j22αj2pj2.subscript¯𝜆subscript𝑗12subscript𝛼subscript𝑗1subscriptsuperscript𝑏subscript𝑗1superscript𝜂2subscript¯𝜆subscript𝑗22subscript𝛼subscript𝑗2subscriptsuperscript𝑝subscript𝑗2\displaystyle\frac{\bar{\lambda}_{j_{1}}+2\alpha_{j_{1}}b^{*}_{j_{1}}}{\eta^{2% }}=\bar{\lambda}_{j_{2}}-2\alpha_{j_{2}}p^{*}_{j_{2}}.divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + 2 italic_α start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT . (18)

Note that intervals j1,j2subscript𝑗1subscript𝑗2j_{1},j_{2}italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT , italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT are partial, bj1,pj2<P¯subscriptsuperscript𝑏subscript𝑗1subscriptsuperscript𝑝subscript𝑗2¯𝑃b^{*}_{j_{1}},p^{*}_{j_{2}}<\bar{P}italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT , italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT < over¯ start_ARG italic_P end_ARG. Equation (18) becomes

λ¯j1+2αj1P¯η2>λ¯j22αj2P¯.subscript¯𝜆subscript𝑗12subscript𝛼subscript𝑗1¯𝑃superscript𝜂2subscript¯𝜆subscript𝑗22subscript𝛼subscript𝑗2¯𝑃\displaystyle\frac{\bar{\lambda}_{j_{1}}+2\alpha_{j_{1}}\bar{P}}{\eta^{2}}>% \bar{\lambda}_{j_{2}}-2\alpha_{j_{2}}\bar{P}.divide start_ARG over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT + 2 italic_α start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG > over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT - 2 italic_α start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT over¯ start_ARG italic_P end_ARG . (19)

For three-interval bidding, the inequality constraint (19) provides a sufficient condition for the scenario where s=2,o=1formulae-sequence𝑠2𝑜1s=2,~{}o=1italic_s = 2 , italic_o = 1 when λ¯j2>λ¯j1,λ¯j2>λ¯t0formulae-sequencesubscript¯𝜆subscript𝑗2subscript¯𝜆subscript𝑗1subscript¯𝜆subscript𝑗2subscript¯𝜆subscript𝑡0\bar{\lambda}_{j_{2}}>\bar{\lambda}_{j_{1}},~{}\bar{\lambda}_{j_{2}}>\bar{% \lambda}_{t_{0}}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT > over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT , over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT > over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT. It also implies when the ultimate price difference is larger between intervals j1subscript𝑗1j_{1}italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and j2subscript𝑗2j_{2}italic_j start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, there can be no less full intervals than partial. Such analysis can be easily extended to other scenarios.

Remark 4

Higher price difference, resulting in higher profitability, complicate the distinction between withholding behavior due to strategic arbitrage versus market power exercise. Both strategies exploit market inefficiencies, however, high profit margin can mask the source of earnings, whether derived from arbitrage or price manipulation.

III-D Understanding Multi-Interval Bidding

TABLE III: Scenarios Price Maker Exercising Market Power by Conducting One-Interval Capacity Withholding
Scenario Taker Maker
asuperscript𝑎a^{\dagger}italic_a start_POSTSUPERSCRIPT † end_POSTSUPERSCRIPT b𝑏bitalic_b c𝑐citalic_c d𝑑ditalic_d e𝑒eitalic_e a𝑎aitalic_a b𝑏bitalic_b c𝑐citalic_c d𝑑ditalic_d e𝑒eitalic_e
b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG downgrade ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG p^superscript^𝑝\hat{p}^{\ddagger}over^ start_ARG italic_p end_ARG start_POSTSUPERSCRIPT ‡ end_POSTSUPERSCRIPT
btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT b^superscript^𝑏\hat{b}^{\ddagger}over^ start_ARG italic_b end_ARG start_POSTSUPERSCRIPT ‡ end_POSTSUPERSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG 00 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG downgrade ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG 00
btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG upgrade ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG upgrade ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG
btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG switch ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG 00
btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT 00 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG
ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG 00 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG 00
btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT 00 b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG 00 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG 00
btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG 00
p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG switch ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG 00
btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT 00 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG 00 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG 00 p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG
btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT 00 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG 00 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
ptsubscript𝑝𝑡p_{t}italic_p start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT 00 p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
btsubscript𝑏𝑡b_{t}italic_b start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT 00 P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG P¯¯𝑃\bar{P}over¯ start_ARG italic_P end_ARG
{a,b,c,d,e}𝑎𝑏𝑐𝑑𝑒\{a,b,c,d,e\}{ italic_a , italic_b , italic_c , italic_d , italic_e } represent different random intervals.
p^^𝑝\hat{p}over^ start_ARG italic_p end_ARG and b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG give the generic notation for the withheld power output, i.e., p^,b^(0,P¯)^𝑝^𝑏0¯𝑃\hat{p},\,\hat{b}\in(0,\bar{P})over^ start_ARG italic_p end_ARG , over^ start_ARG italic_b end_ARG ∈ ( 0 , over¯ start_ARG italic_P end_ARG ).

The ex-post analysis is established in two phases. First, the market operator counts and compares the number of full or partial intervals, as illustrated in Fig.2 summarizing the discussion from the previous subsections. For price takers, as suggested in Propositions 1 and 3, they typically participate in the market with the profiles where the number of full intervals are no more the that of partial ones, including the idle scenarios. On the other hand, for price makers, Proposition 5 and Corollary 1 indicate that most scenarios feature the number of partial intervals no less than that of full ones, suggesting the exercise of market power through capacity withholding. Exceptions occur where the number of full intervals exceed that of partial ones, notably in situations with significant price differentials, as detailed in Remark 4. Secondly, examining price-decision relationship across the intervals reveals the bidding strategy adopted by the market participant, indicating whether it involves the exercise of market power. With this analysis in place, we revisit Theorem 1 and present its proof as follows:

Refer to caption
Figure 2: Bidding pattern regarding capacity withholding as a price taker or a price maker given strictly heterogeneous prices.
Proof:

We first discuss the scenarios for price takers and price makers within the framework of Assumptions 13, then explore more practical cases following the relaxation of these assumptions.

For price takers, from Proposition 3 we’ve covered that given strictly heterogeneous prices, there exists only one partial interval, which establishes a sufficient condition for Theorem 1. Relaxing this constraint to allow for instances of identical prices potentially leads to the presence of additional partial intervals within the scheduling period. Therefore, Theorem 1 is applicable to price takers given these assumptions.

For price makers, Proposition 5 states that the number of partial intervals is always no less than one, i.e., s1𝑠1s\geq 1italic_s ≥ 1. When s=1𝑠1s=1italic_s = 1, assuming small price sensitivity parameter αtsubscript𝛼𝑡\alpha_{t}italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT, we can conclude that the bidding profiles from price makers are identical to those observed among price takers, thus the price maker is not evidently exercising market power. However, scenarios exist where, despite the profile showing only one withholding interval, the price maker is exercising market power as summarized in Table III. We take the scenario b^^𝑏\hat{b}over^ start_ARG italic_b end_ARG downgrade from Table III for example to explain how Theorem 1 can identify the one-interval-withholding market power exercise:

For the price taker, in order to derive the bidding decisions pb=P¯,bc=b^formulae-sequencesuperscriptsubscript𝑝𝑏¯𝑃superscriptsubscript𝑏𝑐^𝑏p_{b}^{*}=\bar{P},\,~{}b_{c}^{*}=\hat{b}italic_p start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = over¯ start_ARG italic_P end_ARG , italic_b start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = over^ start_ARG italic_b end_ARG, the corresponding clearing prices need to satisfy λ^b>λ^cη2subscript^𝜆𝑏subscript^𝜆𝑐superscript𝜂2\hat{\lambda}_{b}>\frac{\hat{\lambda}_{c}}{\eta^{2}}over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT > divide start_ARG over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG as illustrated in Proposition 4. Similarly, for the price maker, following the KKT conditions (16), the bidding decisions pb=p^,bc=0formulae-sequencesuperscriptsubscript𝑝𝑏^𝑝superscriptsubscript𝑏𝑐0p_{b}^{*}=\hat{p},\,~{}b_{c}^{*}=0italic_p start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = over^ start_ARG italic_p end_ARG , italic_b start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = 0 are derived under the condition:

λ^cαc<λ^b+αb(P¯2p^)<λ^cαcη2,subscript^𝜆𝑐subscript𝛼𝑐subscript^𝜆𝑏subscript𝛼𝑏¯𝑃2^𝑝subscript^𝜆𝑐subscript𝛼𝑐superscript𝜂2\displaystyle\hat{\lambda}_{c}-\alpha_{c}<\hat{\lambda}_{b}+\alpha_{b}(\bar{P}% -2\hat{p})<\frac{\hat{\lambda}_{c}-\alpha_{c}}{\eta^{2}},over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT - italic_α start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT < over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT + italic_α start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( over¯ start_ARG italic_P end_ARG - 2 over^ start_ARG italic_p end_ARG ) < divide start_ARG over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT - italic_α start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (20)

and the resulting clearing prices are:

λ^b=λ^b+αb(P¯p^),subscriptsuperscript^𝜆𝑏subscript^𝜆𝑏subscript𝛼𝑏¯𝑃^𝑝\displaystyle\hat{\lambda}^{\prime}_{b}=\hat{\lambda}_{b}+\alpha_{b}(\bar{P}-% \hat{p}),over^ start_ARG italic_λ end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT + italic_α start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( over¯ start_ARG italic_P end_ARG - over^ start_ARG italic_p end_ARG ) ,
λ^c=λ^cαcb^.subscriptsuperscript^𝜆𝑐subscript^𝜆𝑐subscript𝛼𝑐^𝑏\displaystyle\hat{\lambda}^{\prime}_{c}=\hat{\lambda}_{c}-\alpha_{c}\hat{b}.over^ start_ARG italic_λ end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT = over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT - italic_α start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT over^ start_ARG italic_b end_ARG .

Given αb(P¯p^)>0subscript𝛼𝑏¯𝑃^𝑝0\alpha_{b}(\bar{P}-\hat{p})>0italic_α start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT ( over¯ start_ARG italic_P end_ARG - over^ start_ARG italic_p end_ARG ) > 0 and αcb^>0subscript𝛼𝑐^𝑏0\alpha_{c}\hat{b}>0italic_α start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT over^ start_ARG italic_b end_ARG > 0, we have λ^b>λ^cη2subscriptsuperscript^𝜆𝑏superscriptsubscript^𝜆𝑐superscript𝜂2\hat{\lambda}^{\prime}_{b}>\frac{\hat{\lambda}_{c}^{\prime}}{\eta^{2}}over^ start_ARG italic_λ end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT > divide start_ARG over^ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG. Thus, the price-decision relationship between (λ^b,λ^c)subscriptsuperscript^𝜆𝑏subscriptsuperscript^𝜆𝑐(\hat{\lambda}^{\prime}_{b},\hat{\lambda}^{\prime}_{c})( over^ start_ARG italic_λ end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT , over^ start_ARG italic_λ end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT ) and (pb=p^,bc=0)formulae-sequencesuperscriptsubscript𝑝𝑏^𝑝superscriptsubscript𝑏𝑐0(p_{b}^{*}=\hat{p},b_{c}^{*}=0)( italic_p start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = over^ start_ARG italic_p end_ARG , italic_b start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT = 0 ) violates Proposition 4. The proof for the remaining scenarios follow similar arguments.

Therefore, Theorem 1 holds for price makers under the given Assumptions 13. Specially, if the unit is held idle, no market power is exercised.

Next, we generalize the conditions following the relaxation of the initial assumptions. Dropping Assumption 1 generates scenarios where energy storage unit output is limited by its energy capacity, potentially leading to the occurrence of additional partial blocks. This modification does not compromise the validity of Theorem 1. Regarding Assumption 2, incorporating scenarios of active operation during zero-profit periods allows for flexible dispatch decisions pt,bt[0,P¯]subscriptsuperscript𝑝𝑡subscriptsuperscript𝑏𝑡0¯𝑃p^{*}_{t},b^{*}_{t}\in[0,\bar{P}]italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT , italic_b start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT ∈ [ 0 , over¯ start_ARG italic_P end_ARG ], which also does not alter the conclusions drawn in Theorem 1. Assumption 3 is initially established to simplify the analysis of KKT optimality conditions. Negative prices enforce storage unit no discharging according to (5c), therefore pt=0subscriptsuperscript𝑝𝑡0p^{*}_{t}=0italic_p start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT = 0 for certain intervals. The removal of such assumption does not impact the conclusions reached in prior subsections. We defer further discussion of this topic to readers with an interest in exploring it in greater detail. ∎

IV NUMERICAL EXPERIMENTS

We validate the proposed framework for energy storage unit bidding by simulating a 24-interval dispatch in the day-ahead market. Our simulation features a 2.5MW/10MWh storage unit that starts and ends at 50% SoC, with charging and discharging efficiencies set at 0.9. Note that the energy capacity constraint (5e) is taken into account in the simulation. We compare the bidding profiles and resulting profits between price takers and price makers, using price data from a winter day in New York City in 2018, as reported by NYISO [15]. This price series serves as the benchmark for competitive market clearing prices. Nominal clearing prices λ¯tsubscript¯𝜆𝑡\bar{\lambda}_{t}over¯ start_ARG italic_λ end_ARG start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT are calculated by assuming that non-withholding power supply from the price taker is contributed to the market, thereby leading to competitive prices. The price sensitivity parameter αtsubscript𝛼𝑡\alpha_{t}italic_α start_POSTSUBSCRIPT italic_t end_POSTSUBSCRIPT is modeled to be linearly proportional to the competitive price level, reflecting the load level and the slope of the remaining supply function, with an average value of $1.00/MWh/MWh for low market power price makers and $2.00/MWh/MWh for high market power cases. Note that this parameter is amplified for the purposes of validating our framework in the simulation and would be expected to be significantly lower in practical scenarios [12].

The market clearing results are presented in Fig. 3. Figure 3(3(a)) illustrates the capacity withholding behavior of price makers with both low and high market power in comparison to price takers. Figure 3(3(b)) displays the resulting shifts in market clearing prices from the competitive benchmark, highlighting the exercise of market power.

Refer to caption
(a) Storage power output
Refer to caption
(b) Market clearing price
Figure 3: Storage unit control policy and the resulting market clearing price considering different market participants as price maker with low market power (makerlsubscriptmaker𝑙\mathrm{maker}_{l}roman_maker start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT), high market power (makerhsubscriptmaker\mathrm{maker}_{h}roman_maker start_POSTSUBSCRIPT italic_h end_POSTSUBSCRIPT), and price taker: (a) storage power output (positive values mean that unit is discharging), (b) market clearing price.

Figure 3(3(a)) illustrates that, as a price maker, capacity withholding occurs mainly during periods of peak or valley prices to arbitrage price differences. For example, during peak hours, price makers may withhold storage discharging to raise prices, e.g., 5-9pm, as reflected by the corresponding market clearing prices in Fig. 3(3(b)). Conversely, during off-peak hours, they might withhold charging to lower prices. Furthermore, a price maker with a higher degree of market power can exert a more significant influence on the market. This is demonstrated by the more aggressive withholding behavior and the subsequent alterations in clearing prices during both peak and off-peak periods. Notably, even when applying the price maker’s bidding model (15), there are instances where the price maker bids similarly to the price taker, leading to identical power outputs and market clearing prices, e.g., 8am-5pm. This indicates that they are not actively exercising market power during these times. For the price taker, there is only one interval of withholding, i.e., 8-9pm, consistent with Theorem 1. This theorem provides a criterion for identifying evident instances of market power exercise. Specially, in some cases, such as from 6-7pm, the power production from the price maker might exceed that of the price taker, which could seem counterintuitive at first glance but could be interpreted as strategic practices aimed at maximizing profit. The profits generated under these scenarios are summarized in Table IV.

TABLE IV: Profit of Market Participants under Different Level of Market Power
Scenario Price Taker ($currency-dollar\$$) Price Maker ($currency-dollar\$$)
No market power 37.95
Low market power 47.50 42.02
High market power 66.66 49.11

Table IV indicates that the price maker exercises market power to gain additional profit. When possessing a higher level of market power, the price maker’s earnings are 76% greater than in scenarios where they act as a price taker ($49.11 vs $37.95), and 25% greater with a lower level of market power ($42.02 vs $37.95). Interestingly, in markets where a price maker is present, price takers may achieve higher profits than the price maker ($42.02 vs $47.50 and $49.11 vs $66.66), aligning with findings discussed in [16]. This observation highlights the vulnerability of market inefficiencies to strategic exploitation.

V CONCLUSIONS

We examine the multi-interval strategic bidding withholding of energy storage units adopting a self-scheduling model, considering both price takers and price makers. For price makers, we introduce a bidding strategy that anticipates market prices through a price sensitivity analysis. The proposed framework serves as an ex-post market monitoring tool, allowing market operators to distinguish market power exercise from competitive withholding by observing the withholding instances within the bidding profiles and the corresponding clearing prices. Our findings support the economic bidding behaviors commonly seen among storage units and reveal that significant price fluctuations present substantial profit opportunities for both market power holders and strategic profit-seekers. In future work, we aim to enhance our model by incorporating uncertainty in price forecasts into the framework and develop corresponding criteria for market power assessment.

VI Acknowledgements

JA acknowledges funding from NSF grants ECCS 2144634 and ECCS 2047213 as well as the Columbia Data Science Institute. Bolun Xu acknowledges funding from NSF grants ECCS 2239046.

References

  • [1] A. Sakti, A. Botterud, and F. O’Sullivan, “Review of wholesale markets and regulations for advanced energy storage services in the united states: Current status and path forward,” Energy policy, vol. 120, pp. 569–579, 2018.
  • [2] California ISO , “Special report on energy storage,” 2023. [Online]. Available: https://www.caiso.com/Documents/2022-Special-Report-on-Battery-Storage-Jul-7-2023.pdf
  • [3] U.S. Energy Information Administration, “Battery systems on the u.s. power grid are increasingly used to respond to price,” 2022. [Online]. Available: https://www.eia.gov/todayinenergy/detail.php?id=53199
  • [4] W. Ma and B. Xu, “A data-driven nonlinear recharge controller for energy storage in frequency regulation,” in 2021 IEEE Power & Energy Society General Meeting (PESGM), 2021, pp. 1–5.
  • [5] H. Mohsenian-Rad, “Coordinated price-maker operation of large energy storage units in nodal energy markets,” IEEE Transactions on Power Systems, vol. 31, no. 1, pp. 786–797, 2015.
  • [6] Y. Baker, N. Zheng, and B. Xu, “Transferable energy storage bidder,” IEEE Transactions on Power Systems, 2023.
  • [7] C. Chen, L. Tong, and Y. Guo, “Pricing energy storage in real-time market,” in 2021 IEEE Power & Energy Society General Meeting (PESGM), 2021, pp. 1–5.
  • [8] D. Krishnamurthy, C. Uckun, Z. Zhou, P. R. Thimmapuram, and A. Botterud, “Energy storage arbitrage under day-ahead and real-time price uncertainty,” IEEE Transactions on Power Systems, vol. 33, no. 1, pp. 84–93, 2017.
  • [9] R. K. Bansal, P. You, D. F. Gayme, and E. Mallada, “A market mechanism for truthful bidding with energy storage,” Electric Power Systems Research, vol. 211, p. 108284, 2022.
  • [10] R. Baldick, U. Helman, B. F. Hobbs, and R. P. O’Neill, “Design of efficient generation markets,” Proceedings of the IEEE, vol. 93, no. 11, pp. 1998–2012, 2005.
  • [11] S. M. Harvey and W. W. Hogan, “Market power and withholding,” Harvard University, Cambridge, MA, 2001.
  • [12] W. Reinisch and T. Tezuka, “Market power and trading strategies on the electricity market: a market design view,” IEEE Transactions on Power Systems, vol. 21, no. 3, pp. 1180–1190, 2006.
  • [13] R. Baldick, R. Grant, and E. Kahn, “Theory and application of linear supply function equilibrium in electricity markets,” Journal of Regulatory Economics, vol. 25, pp. 143–167, 2004.
  • [14] B. Xu, M. Korpås, and A. Botterud, “Operational valuation of energy storage under multi-stage price uncertainties,” in 2020 59th IEEE Conference on Decision and Control (CDC), 2020, pp. 55–60.
  • [15] New York ISO, “Energy market & operational data,” 2024. [Online]. Available: https://www.nyiso.com/energy-market-operational-data
  • [16] Y. Wu, J. Kim, and J. Anderson, “Mitigation-aware bidding strategies in electricity markets,” in 2023 IEEE Power & Energy Society General Meeting (PESGM), 2023, pp. 1–5.