Abstract
Load serving entities (LSEs) maximise their profit by increasing the difference between revenue from electricity sale to consumers and the cost of procuring that electricity. Considering the offered sale prices, consumers minimise energy bills by scheduling their energy consumption. Varying spot market prices and consumer behaviour impact LSE’s electricity procurement and sale price decisions. The two conflicting objectives of LSE profit maximization and consumer cost minimization can be modelled effectively by hierarchical bilevel programming. Additionally, LSE has to consider spot market price uncertainty and renewable availability during different hours of the day. This paper considers these issues for LSE’s riskbased profit maximization decision making model under RE availability by proposing a bilevel framework to determine optimal dynamic sale prices and energy procurement decisions. The upper level considers riskbased profit maximization for LSE and the lower level addresses the consumer’s objective of cost minimization. This work considers Conditional Value at Risk (CVaR) to model spot market price risk for pragmatic characterisation of LSE’s riskaverse behaviour. A case study on the PJM market shows the effectiveness of the proposed approach.
Introduction
Load serving entity is a retailer or distribution utility that procures energy from wholesale electricity markets to sell it to its consumers. LSE makes optimal decisions for energy procurement and sale prices to maximise its profit and minimise associated risk. The risk in wholesale market trading arises from volatile and uncertain spot market prices [1]. Varying but flexible consumer demand and renewable energy (RE) penetration makes this a challenging problem for LSE. LSE offers dynamic prices to manage flexible demand [2]. In response to the offered dynamic prices, consumers alter/ optimise energy consumption to minimise their energy bill. This optimisation by consumers influences LSE’s energy procurement decisions and eventually, its profit and risk management strategies. While making such riskbased decisions, variability and availability of RE generation creates challenges for its efficient utilisation.
LSE maximises its profit by minimising procurement cost and maximising revenue from electricity sale [1]. Generally, LSE performs various types of forward contracts to eliminate the possibility of high procurement cost arising from unknown future variations in spot market prices [3]. This restricts LSE’s ability to buy at a lower cost when spot market prices are low. As spot prices are uncertain, it introduces variability in procurement cost, which needs to be minimised while making procurement decisions. Financial risk posed by uncertain spot prices necessitates understanding LSE’s behaviour, considering its risk preferences [3]. For such analysis, Conditional Value at Risk (CVaR) is the most advanced and welladopted risk measure [3,4,5,6]. LSE’s revenue can be maximised by offering optimal fixed or timevarying/ dynamic sale prices. Time of use (ToU) and realtime pricing (RTP) are the most popular variants of timevarying pricing schemes [7,8,9].
For LSE’s profit maximisation, the consumer’s price responsive behaviour is often modelled by the price quota curve (PQC) or market share function (MSF) [3, 4]. Consumers’ responsiveness to sale price change is modelled by price elasticity of demand [4, 10]. However, these approaches cannot capture consumers’ response to sale prices when consumers receive sale price information a day or few hours prior to energy consumption. Dynamic prices such as dynamic ToU and RTP signals are communicated well ahead of realtime. Consumers make their consumption decision rationally rather than based on price elasticity, to minimise energy bill [11]. Consequently, LSE’s decisions would change with modified consumer demand. A hierarchical decisionmaking framework can effectively model consumer behaviour in response to offered sale prices under a bilevel programmingbased approach [12].
In bilevel entirety, two conflicting objectives, LSE’s profit maximisation and consumers’ energy bill minimisation, become a joint optimisation problem [13, 14]. Consumers’ optimisation problem is expressed in the lower level of the bilevel programming model. Suppressing/shifting of consumer’s flexible demand modifies its welfare associated with comfort/utility cost. Considering flexible demand and sale prices, consumers tradeoff between energy bill minimisation and welfare/comfort maximisation [11, 14]. Consumers being sensitive to (sale) price change, seek to minimise energy bill irrespective of discomfort [15]. Such LSE’s decisionmaking problems are considered for short or mediumterm planning horizon.
LSE’s decisionmaking problem is expressed in the upper level of the bilevel model. For shortterm planning, LSE manages energy procurement from spot and regulation markets [11, 14]. For mediumterm planning, LSE manages forward contracts portfolio along with procurement from spot market [3, 16]. Additionally, distributed generation (DG) including RE as procurement options would influence LSE’s sale price and procurement decisions, and eventually its profit and risk management strategy [15, 17, 18]. LSE’s sale price and procurement decisionmaking problem under hierarchical framework are investigated targeting optimal utilisation of thermal DG [15]. However, it neglects demand shifting modelling, which is important to exploit economic benefits for both LSE and consumers [18]. Decision of residential consumer is modelled at the appliances level by incorporating typical operational cycle of these appliances [19]. Such a model requires accurate and detailed historical information of each consumer’s energy consumption. With a focus on electric fleet, RTP prices are determined by clustering demand profile [20]. Multiple consumer classes such as plugin electric vehicles and residential consumers are considered together in LSE’s decisionmaking [21]. Consideration of the classes together is practically unacceptable as it leads to the same sale prices for both classes, irrespective of their different demand profile. Though dynamic sale prices in a bilevel framework are determined, spot price uncertainty and risk modelling are neglected [22]. LSE can target an assured level of uncertainty to maximise its profit, considering the fractile model [23]. However, this consideration cannot minimise the impact of profit risk.
Spot market and bilateral contracts are imperative sources for energy procurement in LSE’s decisionmaking. RE sources can be considered as selfgeneration to supply a part of electricity. RE characteristics of availability and intermittency may potentially impact LSE’s decisionmaking. RE being zerocost energy sources, static ToU prices are determined using PQC [8]. Notably, static ToU prices remain stable for a longer duration and are communicated at lower frequencies [19]. Focus of this work is to determine dynamic sale prices for rational consumers. In a bilevel framework, RE’s investment cost is considered to highlight its impact on fixed sale prices [24]. For such analysis, spot price risk and consumer behaviour modelling, which significantly impacts sale price dynamics, are neglected. As a whole, the impact of RE on dynamic sale price optimisation in a hierarchical decisionmaking framework is yet to be investigated. Moreover, for such investigation, analysis on procurement decisions along with price risk modelling is essential to comprehend the riskneutral and riskaverse behaviour of LSE.
In this context, the proposed work presents a hierarchical decisionmaking model to determine optimal dynamic sale prices and energy procurement decisions for an LSE. This considers consumers’ flexible demand under RE availability and spot market price uncertainty. RE cost is considered to optimise dynamic sale prices. The influence of rational consumers’ price responsive behaviour on LSE’s profit maximisation is considered by bilevel optimisation. Spot market, bilateral contracts, and selfgeneration are considered for procuring electricity. Thermal generation unit and RE are considered as selfgeneration. Risk of spot market prices is modelled by CVaR approach to determine riskbased strategies for LSE. Formulated bilevel optimisation problem is transformed into an equivalent singlelevel problem using the classical KarushKuhnTucker (KKT) approach. Obtained nonlinear formulation is linearised by auxiliary variables and strong duality theory. Results from the proposed model help LSE to decide its procurement and selling price strategies, and consumers to determine their optimal consumption pattern.
In this view, key contributions of the paper are

Dynamic sale prices determination considering RE availability and flexible demand in hierarchical decisionmaking framework

Optimal dynamic sale price determination on an hourly basis to manage consumers’ flexible demand

Modelling of hierarchical decisionmaking framework under spot market price risk using CVaR to comprehend riskneutral and riskaverse decisions of LSE.
LSE’s DecisionMaking Framework
This work considers energy trading by a retailer or distribution utility as LSE’s operation. LSE determines optimal energy procurement and dynamic sale prices to be offered to consumers in the mediumterm, to maximise its profit. Considered LSE is responsible for supplying only energy and does not operate distribution networks. Hence, LSE’s economic operation is only targeted in this work. LSE procures energy from the wholesale electricity market through bilateral contracts and spot market. A part of its energy requirement is procured from available thermal and RE selfgeneration. Uncertainty of spot market price introduces volatility in procurement cost, eventually making the profit risky. LSE aims to minimise profit risk according to its risk aversion level. Hence, a riskaverse decisionmaking model in addition to riskneutral model for LSE is considered. Further, consumers can adjust their consumption based on offered sale prices. This priceresponsive behaviour of consumers influences LSE’s decisionmaking problem. Thus, consumers as decisionmakers optimise their consumption pattern to minimise energy procurement cost (energy bill) based on offered sale price dynamics. The decisionmaking problem consists of two decisionmakers (LSE and consumers), each trying to optimise their respective objectives, subject to respective constraints. LSE’s decisionmaking problem is formulated as a bilevel programming problem, as illustrated in Fig. 1.
The upperlevel of bilevel problem describes LSE’s riskbased profit maximisation. To analyse LSE’s riskaverse behaviour, tradeoff between risk and profit for LSE has been modelled using CVaR. The lower level of bilevel is consumers’ cost minimisation problem. Both LSE and consumers optimise their objectives over a jointly dependent set.
The bilevel optimisation problem is solved by transforming it into an equivalent singlelevel problem. For this, the lower level problem is replaced with its KKT optimality conditions. The KKT condition appears as Lagrangian and complementarity slackness constraints in the equivalent singlelevel optimisation problem. The resulting problem is nonlinear, which is linearised by introducing auxiliary variables and strong duality theory.
It is considered that the aggregated outcome from RE would be provided to LSE. RE generation uncertainty consideration is avoided due to aggregation [25]. However, the impact of RE variability and spot market prices is high; hence both are considered in dynamic sale price determination. This work also assumes that consumers respond rationally to the dynamic prices offered by LSE.
BiLevel Mathematical Formulation
The proposed bilevel problem discussed in “LSE’s DecisionMaking Framework” section is mathematically formulated in this section.
Upper Level: LSE’s Problem
LSE’s tradeoff problem, i.e., profit maximisation and risk minimisation, subject to procurement and sale price constraints, is mathematically formulated in the following subsections.
Procurement Cost
Total cost for energy procurement in hour t is equal to the sum of energy cost from bilateral contracts \( {C}_t^{bc}, \) thermal generating units \( {C}_t^{TG} \), renewable generation \( {C}_t^{PV} \), and spot market \( {C}_t^S \).

(a)
Bilateral contract
LSE signs bilateral contracts with electricity suppliers. It procures energy at a mutually agreed price \( {\lambda}_{i,t}^{bc} \) for a specified time t from i^{th} bilateral contract. Energy procurement cost incurred due to bilateral contracts in hour t is given by (2). Constraint (3) imposes minimum and maximum energy procurement limits from individual bilateral contracts if the contract is exercised as indicated by a binary variable δ_{i, t}.

(b)
Selfgeneration
LSE utilises available distributed energy resources (DER) as a procurement option. Conventional thermal generation and renewable generation from solar photovoltaic (PV) systems are considered as available DER facilities to LSE. At time t, total energy procured from selfgeneration \( {P}_t^{SG} \) would be equal to energy procured from thermal and PV generation (4).

(i)
Thermal generation
The operational cost of thermal generation \( {C}_t^{TG} \) considers fuel cost, startup and shut down costs. Cost of thermal selfgeneration is
The fuel cost of thermal generating units is modelled by a quadratic function and is given by (6). Piecewise linearisation of cost function is shown in Appendix 1. This linearisation is required to formulate a linear bilevel problem.
\( {C}_t^{TG} \) is subject to the following constraints
Constraints (7) and (8) decide start up and shut down cost of generating units. Constraints (9) and (10) decide ramp up and ramp down limits. Minimum up and down time of generating units is given by constraints (11) and (12). Minimum & maximum generation limits are decided by constraint (13). (14) is a nonnegativity constraint. (15) is a variable declaration constraint.

(ii)
Renewable generation
RE generation from PV system is considered. Considering the nondispatchable characteristics of RE, this generation is considered a must run. It provides \( {PV}_t^{act} \)generation. The energy procurement cost from PV systems is considered using (16). This cost recovers the investment cost of PV systems. LSE incurs some cost when it procures RE through contracts. RE cost is considered in this work considering these factors though fuel cost for RE generation is negligible. The value of \( {\lambda}_t^{PV} \) can be mathematically calculated using bottomup or any other methodology by considering relevant cost parameters of PV systems. Calculation of \( {\lambda}_t^{PV} \) is beyond the scope of the work, and so its value is assumed. Available RE generation could be curtailed during hours of low demand using (17). This constraint is required as available RE generation and demand cannot be matched by the limited flexible demand of consumers.

(iii)
Spot market
LSE considers spot market as one of the procurement options. As future prices in the spot market are uncertain, its forecasted/ expected values are considered. Expected energy purchase cost from the spot market for each hour t is given by (18). This work considers procurement from the spot market only.
Constraint (19) represents energy balance equation to ensure that energy procured from various sources meets the consumer demand in each hour t. Here \( {L}_t^{act} \) is consumers’ optimised demand, obtained from the lowerlevel problem.
Revenue
LSE resells acquired energy to its consumers at a per unit sale price \( {\lambda}_t^{sale} \). Revenue for each time period is determined by (20), subject to constraints (21)–(23).
Constraints (21) and (22) restrict the sale price within a certain percentage of nominal sale price at hour t. The constraint given by (23) is used to impose an upper bound on the sale prices. This constraint ensures sufficient number of low sale prices hours. It also helps to alleviate LSE’s market power [14]. In the absence of (23), LSE’s optimal choice will touch the maximum limit of sale price, i.e.,\( {\lambda}_t^{sale}=\left(1+{z}_t^{max}\right){\lambda}_t^{nsale} \), which seems to be unfair to consumers. It is assumed that LSE and consumers have mutually agreed upon considered minimum, maximum and average sale prices.
Risk Modelling
Spot market prices are highly uncertain and can not be predicted with certainty beforehand during decisionmaking. Uncertainty of spot prices introduces risk for LSE. This risk in LSE’s energy procurement is quantified and measured using CVaR. Risk measure CVaR is defined as the conditional expectation of losses that exceeds a threshold value: value at risk (VaR). VaR provides the expected maximum loss over a target horizon within a confidence interval. CVaR is computed by taking a weighted average of the losses exceeding VaR values. In this work, historical simulation method of CVaR computation is applied [26]. CVaR computation algorithm is shown in Appendix 2. Risk incorporated by spot market procurement is accounted by (24). \( CVaR\left({\lambda}_t^S\right) \) represents calculated CVaR for per unit procurement for given confidence from the historical simulation method.
LSE’s Profit and Risk Tradeoff
LSE’s objective is to maximise expected profit and minimise risk imposed by spot prices. LSE determines a tradeoff between expected profit and risk considering its risk preferences. This tradeoff is formulated based on the CVaR criterion in the upper level of the bilevel problem. LSE’s multiobjective function is formulated by combining profit and risk objectives using a risk weighting parameter β. The value of β describes risk preferences, i.e., riskaverse and riskneutral behaviour of LSE. Riskneutral behaviour describes LSE’s concerns for only profit maximisation and is obtained by assigning β = 0. The riskaverse behaviour of LSE indicates its emphasis on risk and is accomplished by assigning higher values of β. The LSE’s multiobjective function is given by
The first term in (25), E(Prof_{t}) represents LSE’s expected profit and is obtained by subtraction of expected procurement cost (1) from revenue generated (20), as given by (26). The second term in (25) represents risk, which is to be minimised.
Lower Level: Consumers’ Problem
Consumers are the second decisionmaker of the proposed bilevel problem. Hence, at the lower level, consumers determine their optimal energy consumption based on dynamic sale prices offered by LSE. Consumers received dynamic price signals from LSE. LSE determines these prices in the upper level of the proposed bilevel problem. A group of consumers with similar characteristics can be represented by different classes. For these consumer classes, mathematical formulation remains the same. However, the available flexible demand and nominal demand profiles of these consumer classes will be different. This is reflected in the mathematical model through appropriate factors. Moreover, consumer is aware of its nominal/ initial demand profile and available flexible demand. Consumer objective is mathematically formulated as
subject to constraints
In (27), \( {L}_t^{act} \) is revised energy demand, which is a decision variable of considered consumer class. Inequality constraints (28) and (29) describe the minimum and maximum limits on consumer demand in terms of their initial demand [14]. These constraints represent a relation between revised demand, initial demand, and available flexible demand of consumers. The minimum and maximum limit on energy (electricity) consumption is determined by multiplying a factor x^{l} and x^{u} with initial demand, respectively. These factors represent consumer’s flexible demand in percentage terms. As mentioned before, consumers have information about their initial demand profile and flexible demand. Consumers aim to redistribute their demand to minimise energy bill. Hence, this work considers that consumers do not reduce or suppress their total demand. However, they can shift it to other hours considering their flexible demand while maintaining the total demand at its original value. In other words, the sum of L_{t }over 24 hours equals the sum of \( {L}_t^{act} \). This demand shift is mathematically formulated as (30).
Equivalent SingleLevel Programming Problem
The bilevel problem formulated for LSE in (1)–(30) can be solved by transforming it to one level. This transformation is done by adding equivalent KKT optimality conditions [27] of lowerlevel problems (27)–(30) as constraints to the upperlevel problem. This includes both necessary and sufficient conditions for optimality [27]. Resultant problem contains complementary and slackness conditions, making it nonlinear. The nonlinear complementary conditions are replaced by their equivalent linear expressions to make the problem linear.
KKT Optimality Condition for LowerLevel Problem
The KKT optimality conditions are obtained by determining Lagrangian functions for the lower level problem (27)–(30), given by (31) [17, 27]. σ, γ_{t} and μ_{t} in (31) represent the Lagrange multiplier associated with the lowerlevel problem constraints (28)–(30), respectively.
KKT necessary optimality conditions of the lower level problem are obtained by partial derivatives of the Lagrangian function (31). KKT necessary optimality conditions consist of stationarity (32), dual feasibility (33)–(34) and complementary slackness conditions (35)–(36), respectively. These conditions are in addition to the primal feasibility constraints (28)–(30). Therefore, the lowerlevel problem is represented by (28)–(30) and (35)–(36). To this end, the bilevel optimization problem is transformed into a nonlinear complementary model.
Equivalent Linear Expression
The nonlinear complementary slackness conditions (35)–(36) obtained from KKT conditions are linearised by introducing an auxiliary binary variable (w_{t }and υ_{t}) and a sufficiently positive large constant (M_{1}and M_{2}) [27]. This results in two constraints per complementarity constraints as
Constraints (37)–(39) replace (35) and constraints (40)–(42) replace (36).
Product of consumer demand and sale prices, i.e., \( {L}_t^{act}{\lambda}_t^{sale} \) in (27) introduces nonlinearity in the problem. Its equivalent linear expression is determined by utilising duality theory. The dual of lower level problem (27)–(30) is obtained and shown in (43).
where σ, γ_{t} and μ_{t} are the dual variables associated with the constraints (28)–(30). Optimal solution is obtained by equating primal (27) and dual (43) as
It is worth pointing out here that (44) does not connect lower and upper levels. The nonlinear term \( {L}_t^{act}{\lambda}_t^{sale} \) in (26) is replaced by its equivalent linear term represented in (44). Eq. (26) can be rewritten as
Objective function of the equivalent single level problem is obtained by substituting (45) in (26) as
The equivalent singlelevel formulation of the bilevel problem consists of is a Mixed Integer Linear Programming problem (MILP). Objective function (46) is to be maximised subject to procurement (1)–(19), sale price (21)–(23), CVaR (24), consumers constraints (28)–(30), and KKT (32)–(34), (37)–(42).
Case Study
Data
The proposed formulation for an LSE is illustrated via a case study. Historical spot price data from the PJM electricity market are considered for this study [28]. RE generation profile for a solar based PV system in the Pennsylvania region [29] is considered (shown in Fig. 2). A fixed price of 38$/MWh is considered for energy procurement from RE. It is assumed that one thermal selfgenerating unit is available with LSE. Specification for considered unit is shown in Table 1. It is worth noting that only one thermal generating unit is there to highlight the impact of RE availability. One bilateral contract at 35$/MWh is considered with maximum and minimum procurement limits of 300 MW and 30 MW, respectively. Dynamic sale prices and energy procurement are determined for each hour of a day for LSE.
This case study considers one type of consumer class to illustrate performance of the proposed model. The model can be applied to other consumer classes with their respective demand profiles and flexibility. Forecasted demand profile of considered consumer class is depicted in Fig. 3. It is considered that consumers can vary their demand up to 15% of their original demand at each hour. Nominal sale price for any hour is set 5% higher than corresponding spot market price. An average sale price is taken equal to 38 $/MWh. To keep sale prices within a reasonable limit, LSE considers minimum value of sale prices equal to nominal sale prices and the maximum sale price 20% higher than the nominal sale price at each hour.
Statistical Calculations
Expected values of spot market prices for each hour are calculated as a mean of historical data (from 1/12/2016 to 11/03/2017) corresponding to that hour, over the entire planning period. CVaR value for 95% confidence level is computed as per unit value of historical spot market price data. The historical simulation method, as described in Appendix 2, is utilised to obtain per unit CVaR values. Fig. 4 depicts calculated expected spot prices and CVaR values over 24 hours.
Simulations
The objective function (46), subject to constraints (1)–(24), (28)–(30), (32)–(34), and (37)–(42) is simulated for different values of risk weighting parameter β to maximise LSE’s profit. The resulting MILP problem has been solved using CPLEX solver commercially available in GAMS® software on Intel®, Core^{(}™^{)}, i5 CPU, 3.20 GHz and 4GB of RAM system [30]. CPLEX is a highperformance mathematical programming solver for mixed integer linear programming (MILP) [30]. It uses branch and cut algorithm for solving MILP problems. The solution of the proposed problem generates 508 variables and 96 discrete variables. Average execution and computation time to solve the optimisation problem are 0.089 s and 2.88 s, respectively.
Results
Considering LSE’s risk preferences, i.e., riskneutral (β = 0) and riskaverse (β > 0), results evaluated from simulations are illustrated in Figs. 5, 6, 7, 8 and 9. Highest riskaverse behaviour of LSE is obtained for β = 1.5. Results include dynamic sale prices, optimised consumer demand and energy procurement portfolio.
Dynamics of optimised sale prices and flexible consumer demand is shown in Figs. 5 and 6, respectively, for riskneutral and riskaverse LSE. Curves for sale price and demand labelled with β = 0 and β = 1.5 represent riskneutral and riskaverse behaviour of LSE. The dynamics of consumers’ optimised demand indicates that consumers modify their energy consumption during each hour to minimise their overall energy cost (Figs. 5 and 6). In other words, consumer response follows the sale price signal offered by the LSE, which provides them with a minimum cost (energy bill). It shifts most of the flexible demand to hours where sale prices are low. This demand shift also incetivizes LSE to procure most of its demand during low market price hours.
Figure 5 shows that the dynamics of offered sale prices for riskneutral and riskaverse case differs. This difference in sale price dynamics is governed by wholesale prices and spot price risk. As riskneutral LSE (β = 0) ignore spot price risk, it offers higher value of sale prices at hours when procurement cost tends to be high. However, a riskaverse LSE considers a given level of spot price risk and tries to offer high sale prices at hours when spot price risk is high. In both the cases, LSE does this to shift its consumers’ flexible demand from these hours to other beneficial hours. Consumers respond to these prices, and this results in different consumption patterns for riskneutral and riskaverse cases (Fig. 6). At certain hours, sale prices are low, though consumers reduce their demand (hours 20, 21 for riskneutral case). This happens because consumers can benefit by shifting demand to hours when sale prices are comparatively low (hour 15–17) (Figs. 5 and 6).
Energy procurement from the spot market, bilateral contracts, and thermal and RE selfgeneration obtained for riskneutral and riskaverse LSE are shown in Figs. 7 and 8, respectively. These energy allocations are determined for corresponding optimised demand of consumers (Fig. 6). Figures 7 and 8 show that riskneutral LSE procures most of the energy from spot markets, whereas riskaverse LSE procures most of the energy from sources without uncertainty. It can be observed from obtained energy allocation that available RE generation is completely utilised in both cases. This utilisation significantly impacts energy procurement from costeffective sources during its availability. To decrease overall procurement cost, riskneutral LSE tries to increase consumer demand during RE availability hours where spot prices are low (Fig. 7). For instance, at hours 13–17, consumer demand increases and is procured from the spot market due to their low prices. LSE offers lower sale prices to consumers to increase consumer demand during these hours (Fig. 5). Riskaverse LSE also tries to increase consumer demand during RE availability hours, to minimise risk imposed by uncertain spot prices. This increases procurement from low uncertainty sources, which is low during RE availability (Fig. 8). LSE offers lower sale prices during RE availability and increases with high risk (Fig. 5).
Figure 9 illustrates the efficient frontier depicting riskaverse behaviour of LSE for different values of β. Frontier implies that for a high value of β, LSE’s expected profit and associated risk decreases. Depending upon the willingness to accept spot price risk, LSE may select decisions corresponding to that. In case the value of βis too high, (following the frontier pattern, it is implied that) it would reduce risk and expected profit to zero, which indicates that LSE’s participation in retail activities is not recommended. This means riskaverse LSE can reduce risk up to a certain level.
Conclusion
This paper models hierarchical interaction between LSE and consumers using bilevel programming. Dynamic sale prices and decisions on energy procurement are determined subject to spot market price uncertainty and RE availability. This hierarchical interaction is modelled using bilevel programming. LSE’s objective of profit maximisation and risk minimisation is considered at the upper level. Consumers’ objective of optimising cost by altering its consumption pattern is considered at the lower level. Obtained sale prices and procurement decisions are analysed for riskneutral and riskaverse LSE.
It is observed from the case study that riskneutral LSE maximises its profit to provide price incentives (low sale prices) to consumers to shift their demand to low spot market prices periods. However, riskaverse LSE provides price incentives when risk is low. Sale price dynamics is also impacted by RE utilisation. Both riskneutral and riskaverse LSE offer sale prices to increase consumer demand during RE availability duration to utilise RE generation. However, as risk averseness increases, LSE shifts its procurement from risky spot market to riskfree contracts. Thus, LSE effectively utilises RE by actively using consumers’ flexible demand to achieve its objective.
This work is formulated for an LSE, and its market power is alleviated by restricting the sale price. The proposed model could potentially be extended by incorporating multiple LSEs.
Abbreviations
 g, G :

Index and set for generating units
 i, I :

Index and set for bilateral contracts
 t, T :

Index and set for time
 k :

Index for blocks for piecewise linear cost function
 a _{g,} b _{g,} c _{ g } :

Cost coefficients of g^{th} generating unit
 \( {C}_{g,k}^{ini},{C}_{g,k}^{fin} \) :

Initial and final costs of k^{th} block of linearised cost of g^{th }generating unit [$]
 \( {C}_g^{su},{C}_g^{sd} \) :

Constant startup and shut down cost of g^{th }generating unit [$]
 L _{ t } :

Consumers’ forecasted demand in hour t [MWh]
 n :

Number of blocks in piecewise linear fuel cost function
 \( {P}_{g,k}^{ini},{P}_{g,k}^{fin} \) :

Initial and final power of k^{th} block of g^{th }generating unit [MW]
 \( {P}_g^{min},{P}_g^{max} \) :

Min. and max. Capacity of g^{th }generating unit [MW]
 PV _{ t } :

PV generation in hour t [MW]
 \( {R}_g^u,{R}_g^d \) :

Ramp up and down limit of g^{th}generating unit [MW/h]
 s _{g, k} :

Slope of cost in k^{th} block of linearised cost for g^{th }generating unit [$/MW]
 \( {T}_g^u,{T}_g^d \) :

Minimum up and down time of g^{th }generating unit [h]
 \( {z}_t^{min},{\mathrm{z}}_t^{max} \) :

Min. and max. limits on LSE’s sale price in hour t[%]
 β :

Risk weighting factor
 \( {\lambda}_{i,t}^{bc} \) :

Price of i^{th} bilateral contract in hour t [$/MWh]
 \( {\lambda}_t^{nsale} \) :

Nominal sale price in hour t [$/MWh]
 \( {\lambda}_t^{PV} \) :

Price of PV generation in hour t [$/MWh]
 \( {\lambda}_t^S \) :

Expected/Forecasted spot price in hour t [$/MWh]
 ΔP _{g, k} :

Length of k^{th} block of linearised cost of g^{th }generating unit [MW]
 \( {C}_t^{bc} \) :

Cost of electricity from bilateral contracts in hour t
 \( {C}_t^S \) :

Cost of electricity procurement from spot market in hour t
 \( {C}_t^{TG} \) :

Cost of electricity procurement from thermal selfgeneration in hour t
 \( {C}_t^{PV} \) :

Cost of electricity procurement from PV generation in hour t
 \( {c}_{g,t}^{su},{c}_{g,t}^{sd} \) :

Startup and shut down costs of g^{th }generating unit in hour t [$]
 \( {L}_t^{act} \) :

Consumers’ optimised demand in hour t [MWh]
 \( {P}_{i,t}^{bc} \) :

Energy procured from i^{th }bilateral contract in hour t [MWh]
 P _{g, t} :

Energy generated from g^{th }generation unit in hour t [MWh]
 \( {P}_t^S \) :

Energy procured from spot market in hour t [MWh]
 \( {PV}_t^{sch} \) :

Actual energy scheduled by LSE from PV generation in hour t [MWh]
 \( {PV}_t^{cur} \) :

PV generation curtailment in hour t [MWh]
 \( {X}_{g,t}^{ON},{X}_{g,t}^{OFF} \) :

Variable representing ON/OFF time of g^{th }generating unit in hour t
 \( {\lambda}_t^{sale} \) :

Sale price offered to consumers in hour t [$/MWh]
 u _{g, t} :

ON/OFF (0/1) status of g^{th }generation unit in hour t
 δ _{i, t} :

Binary variable, equal to 1 if i^{th} bilateral contract in hour t is exercised, otherwise zero
References
 1.
Yang J, Zhao J, Luo F et al (2018) Decisionmaking for electricity retailer: a brief survey. IEEE Trans Smart Grid 9(5):4140–4153
 2.
O’Connell N, Pinson P, Madsen H et al (2014) Benefits and challenges of electrical demand response: a critical review. Renew Sust Energ Rev 39:686–699
 3.
Carrion M, Conejo AJ, Arroyo JM (2007) Forward contracting and selling price determination for a retailer. IEEE Trans Power Syst 22(4):2105–2114
 4.
Hatami A, Seifi H, SheikhElEslami MK (2011) A stochasticbased decisionmaking framework for an electricity retailer: timeofuse pricing and electricity portfolio optimisation. IEEE Trans Power Syst 26(4):1808–1816
 5.
Jamshidi M, Kebriaei H, SheikhElEslami MK (2018) An intervalbased stochastic dominance approach for decision making in forward contracts of electricity market. Energy. 158:383–395
 6.
Zhao H, Song S, Zhang Y, Gupta JND, Devlin AG (2019) Optimal decisions of a supply chain with a riskaverse retailer and portfolio contracts. IEEE Access 7:123877–123892
 7.
Joseph S, Erakkath Abdu J (2018) Realtime retail price determination in smart grid from realtime load profile. Int Trans Electr Energ Syst:e2509. https://doi.org/10.1002/etep.2509
 8.
Nojavan S, Zare K, MohammadiIvatloo B (2017) Optimal stochastic energy management of retailer based on selling price determination under smart grid environment in the presence of demand response program. Appl Energy 187:449–464
 9.
Nojavan S, Zare K (2018) Optimal energy pricing for consumers by electricity retailer. Int J Electr Power Energy Syst 102:401–412
 10.
Yusta JM, RamirezRosado IJ, DominguezNavarro JA et al (2005) Optimal electricity price calculation model for retailers in a deregulated market. Int J Elect Power Energy Syst 27(5–6):437–447
 11.
Sekizaki S, Nishizaki I, Hayashida T (2016) Electricity retail market model with flexible price settings and elastic pricebased demand responses by consumers in distribution network. Int J Electr Power Energy Syst 81:371–386
 12.
Sinha A, Malo P, Deb K (2018) A review on bilevel optimisation: from classical to evolutionary approaches and applications. IEEE Trans Evol Comp 22(2):276–295
 13.
Carrión M, Arroyo JM, Conejo AJ (2009) A bilevel stochastic programming approach for retailer futures market trading. IEEE Trans Power Syst 24(3):1446–1456
 14.
Zugno M, Morales JM, Pinson P, Madsen H (2013) A bilevel model for electricity retailer participation in a demand response market environment. Energy Econ 36:182–197
 15.
Khojasteh M, Jadid S (2015) Decisionmaking framework for supplying electricity from distributed generationowning retailers to pricesensitive customers. Util Policy 37:1–12
 16.
Nojavan S, MohammadiIvatloo B, Zare K (2015) Optimal bidding strategy of electricity retailers using robust optimisation approach considering timeofuse rate demand response programs under market price uncertainties. IET Genr Tran Dist 9(4):328–338
 17.
Khazaei H, Vahidi B, Hosseinian SH, Rastegar H (2015) Twolevel decisionmaking model for a distribution company in dayahead market. IET Genr Tran Dist 9(12):1308–1315
 18.
Wei W, Liu F, Mei S (2015) Energy pricing and dispatch for smart grid retailers under demand response and market price uncertainty. IEEE Trans Smart Grid 6(3):1364–1374
 19.
Alves MJ, Antunes CH (2018) A semivectorial bilevel programming approach to optimise electricity dynamic timeofuse retail pricing. Comput Oper Res 92:130–144
 20.
Badri A, Hoseinpour LK (2016) A shortterm optimal decision making framework of an electricity retailer considering optimised EVs charging model. Int Trans Electr Energ Syst 26(8):1705–1724
 21.
RashidizadehKermani H, VahedipourDahraie M, AnvariMoghaddam A, Guerrero JM (2019) A stochastic bilevel decisionmaking framework for a loadserving entity in dayahead and balancing markets. Int Trans Electr Energ Syst:e12109. https://doi.org/10.1002/20507038.12109
 22.
Alekseeva E, Brotcorne L, Lepaul S et al (2018) A bilevel approach to optimise electricity prices. Yug J Opl Res 29(1):9–30
 23.
Sekizaki S, Nishizaki I (2019) Decision making of electricity retailer with multiple channels of purchase based on fractile criterion with rational responses of consumers. Int J Electr Power Energy Syst 105:877–893
 24.
Golmohamadi H, Keypour R (2018) A bilevel robust optimisation model to determine retail electricity price in presence of a significant number of invisible solar sites. Sust Energy, Grids Net 13:93–111
 25.
Bird L, Milligan M, Lew D (2013) Integrating variable renewable energy: Challenges and solutions. National Renew. Energy Lab., NREL/TP6A20–60451.
 26.
Dahlgren R, Liu C, Lawarree J (2003) Risk assessment in energy trading. IEEE Trans Power Syst 18(2):503–511
 27.
FortunyAmat J, McCarl B (1981) A representation and economic interpretation of a twolevel programming problem. J Opl Res Soc 32(9):783–792
 28.
PJM Price data (2020). [Online]. Available: http://www.pjm.com/marketsandoperations/energy.aspx. Accessed 1 Nov 2020
 29.
Solar Data (2020). [Online]. Available; https://www.renewables.ninja/. Accessed 1 Nov 2020
 30.
The GAMS info (2020). [Online]. Available: http://www.gams.com/. Accessed 1 Nov 2020
 31.
Soroudi A (2017) Introduction to programming in GAMS. Springer. https://doi.org/10.1007/9783319623504
Author information
Affiliations
Corresponding author
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
ESM 1
(DOCX 102 kb)
Rights and permissions
About this article
Cite this article
Chawda, S., Mathuria, P. & Bhakar, . BiLevel Approach for Load Serving Entity’s Sale Price Determination under Spot Price Uncertainty and Renewable Availability. Technol Econ Smart Grids Sustain Energy 6, 18 (2021). https://doi.org/10.1007/s4086602100112z
Received:
Accepted:
Published:
Keywords
 Dynamic sale price
 Hierarchical decision making
 Load serving entity
 Price uncertainty
 Renewable energy