Next Article in Journal
Use of Available Daylight to Improve Short-Term Load Forecasting Accuracy
Next Article in Special Issue
A Comprehensive Analysis of the Voltage Unbalance Factor in PV and EV Rich Non-Synthetic Low Voltage Distribution Networks
Previous Article in Journal
Bioenergy Conversion Potential of Decaying Hardwoods
Previous Article in Special Issue
Reliability Assessment of a Fault-Tolerant PV Multistring Inverter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Integrated Profit Maximization Model for Retailers under Varied Penetration Levels of Photovoltaic Systems

by
Ioannis P. Panapakidis
1,*,
Nikolaos Koltsaklis
2 and
Georgios C. Christoforidis
2
1
Department of Electrical and Computer Engineering, University of Thessaly, 38221 Volos, Greece
2
Department of Electrical and Computer Engineering, University of Western Macedonia, 50100 Kozani, Greece
*
Author to whom correspondence should be addressed.
Energies 2021, 14(1), 92; https://doi.org/10.3390/en14010092
Submission received: 17 November 2020 / Revised: 21 December 2020 / Accepted: 23 December 2020 / Published: 26 December 2020
(This article belongs to the Special Issue Smart Photovoltaic Energy Systems for a Sustainable Future)

Abstract

:
In contemporary energy markets, the Retailer acts as the intermediate between the generation and demand sectors. The scope of the Retailer is to maximize its profits by selecting the appropriate procurement mechanism and selling price to the consumers. The wholesale market operation influences the profits since the mix of generation plants determines the system marginal price (SMP). In the related literature, the SMP is treated as a stochastic variable, and the wholesale market conditions are not taken into account. The present paper presents a novel methodology that aims at connecting the wholesale and retail market operations from a Retailer’s perspective. A wholesale market clearing problem is formulated and solved. The scope is to examine how different photovoltaics (PV) penetration levels in the generation side influences the profits of the Retailer and the selling prices to the consumers. The resulting SMPs are used as inputs in a retailer profit maximization problem. This approach allows the Retailer to minimize economic risks and maximize profits. The results indicate that different PV implementation levels on the generation side highly influences the profits and the selling prices.

1. Introduction

The gradual liberalization and development of a single and competitive energy market has been one of the main pillars of the European Union (EU) energy policy, in the context of the wider changes in the energy sector in past few years [1]. The creation of a liberalized electricity market was considered a key priority by the EU institutions, as it was considered that is an important step towards the completion of the internal energy market, more efficient production, transmission and distribution of electricity, enhancing the security of electricity supply and strengthening the competitiveness of the European economy in conjunction with the protection of the environment [2]. For the last 20 years, the EU has sought a methodical opening of the energy market by distinguishing between competitive, e.g., electricity supply to consumers and non-competitive e.g., grid operation activities, third-party access to their infrastructure, the liberalization of electricity supply activity, the gradual lifting of restrictions on consumer choice of retailer, and the establishment of independent market regulators [3,4,5].
Formerly, the organization of energy companies under the status of legal monopoly, with vertically integrated enterprises, identified the multilateralism of generation, transmission, and distribution activities with the one-sidedness of only one provider of all of these services. With the liberalization of energy markets, an accelerating polymerization manifested into the competing activities of generation and distribution. Commonly, the transition from the monopolizing energy industry in the liberalized energy market is done without the activities of the industry to change numerically, but establishing now the possibility of changes in the number of energy actors in the generation and distribution sectors. The result of these changes was an increase of competitive activities, and the transformation of the state into a regulatory supervisor, who sets the criteria and grants the license to deal with these activities for everyone [6,7].
The Retailer is a new entity in the competitive electricity market. Acting as the intermediate agent between generation and supply, the Retailer has only rights for supplying the demand and not purchase energy through interconnections and bid it in the wholesale market [8]. Within a competitive retail market environment, the scope of the Retailer is to offer competitive energy services to its customers to maximize its profits. The Retailer faces many sources of uncertainties such as the selection of electricity procurement mechanism, market prices, consumers’ behavior, and others. The maximization of profits can be accomplished with the selection of the optimal procurement mechanism, and the definition of the optimal selling price offered to the consumers. The maximization of the profit is an applied mathematical optimization task. Subject to competition, the selling price should lead simultaneously to increased profits and consumers satisfaction [9,10].
There is a variety of studies in the literature that investigate the problem of profit maximization. A Monte Carlo based simulation to extract future loads in order to analyze a set of retailer strategies for balancing settlement risks is used in [11]. The scope is to determine the different levels of demand responses to the profit. The problem of optimal contract design referring to prices and quantities, both at supply and end-user levels is studied in [12]. The proposed methodology allows the Retailer to make robust contractual decisions referred to prices and power quantities. Pool market prices are simulated via an ARIMA model in [13]. Apart from the pool market, the procurement mechanisms include forward contracts (FCs). The model is applied to three types of consumers. A sensitivity analysis of the factors that influence profit takes place in [14]. The parameters are price strategy, i.e., variations between periods, upper price limits, consumer elasticity on prices, and others. An investigation of how the retailer’s decisions related to the risk of pool prices influence the profit is carried out in [15]. The conclusions of the study denote that as the Retailer relies more on the pool market, the profits are increased. This is due to the fact that FC prices are higher than the average pool market price. Pool market prices are derived by a Generalized Autoregressive Conditional Heteroskedasticity (GARCH) model in [16]. The Retailer’s approach towards risk is modeled using the conditional value-at-risk (CVaR) measure. The authors examine how the profit changes taking into account different combinations of procurement mechanisms. In [17], the authors focus on the Chinese electricity market. The Retailer model considers fixed prices and dynamic prices. The procurement mechanisms include bilateral mid-long-term contracts and pool market. The authors conclude that dynamic prices lead to higher profits. However, the difference in profits using fixed and dynamic prices becomes lower when the Retailer increases the amount of purchased electricity from the contracts. In [18], the CVaR measure is employed to balance the amounts between the pool and forward markets. The authors consider the linear price/demand function. For each hour, market prices are calculated as the average of historical prices for that hour, over the entire planning period. In [19], the pricing mechanism is time-of-use (TOU) rates. Market prices are modeled as stochastic variables. The paper introduces the time-series responses of the consumers to the time-varying TOU prices in medium-term planning. In [20], both consumers’ demand and system marginal prices (SMPs) are extracted via a SARIMA model. To reduce the number of scenarios, the K-means clustering algorithm is used. The pool market price series are divided per time period within the day, namely peak, shoulder, and valley periods. The authors conclude that purchasing electricity from various mechanisms results in higher profits than utilizing sole purchase mechanisms. In [21], the focus is placed in the comparison of Value-at-Risk (Var) and CVaR measures considering different FC structures. The CVaR wins the competition and leads to lower profit losses. In [22], the examined electricity market is PJM in the USA. The SARIMA model is utilized to develop 500 scenarios for demand and 500 for price. A discussion is provided on how profit is modified between the scenarios. In [23], market prices and consumers’ demand are assumed known variables and not scenarios-based modeling takes place. No FCs are regarded. The paper considers different Retailers that provide different fixed prices. In [24], the Markowitz theory is regarded to optimize Retailer’s portfolio. The scope is to determine the optimal price per consumer categories, namely, residential, industrial, street lighting, small and large commercial. The results show that the retailer, by being more realistic in choosing consumers to its portfolio, can offer more competitive tariffs to key consumers and keep the portfolio optimal and stable in relation to the risk–return ratio. In [25], the aim is to define the optimal selling price. The Retailer owns a hydrogen storage system containing an electrolyzer, hydrogen storage tanks, and fuel cell. Demand and prices are modeled via a series of scenarios. The power unit, together with the storage, are used as additional mechanisms to cover the demand. In [26], prices and loads are stochastic variables and are extracted through a SARIMA model. The initial 500 scenarios are reduced to 100. Apart from CVaR, the authors examine the application of minimax regret and chance-constrained measures.
It should be mentioned that most of the previous studies consider a fixed selling price. A more accurate pricing mechanism that is able to transfer the costs of generated electricity to the end consumers is real-time pricing (RTP) [27]. RTP is a concept that belongs to the family of price-based demand response (DR) methods. Contrary to an incentive-based DR, the consumers participation and response in a price-based DR is not mandatory [28]. The RTP scheme is examined considering various price/demand functions in [29]. The authors examine the influence of the selection of different functions in the profits. Load profiling information is used in [30,31]. The general approach of these studies is: A clustering is applied to a group of consumers, and clusters of consumers with similar characteristics are formed. The profit maximization problem per cluster is solved, and the selling price per consumer clusters is defined.
Based on the previous literature overview, the main conclusions that can be drawn are:
(i)
All the studies concerned with the Retailer profit maximization problem do not involve wholesale market clearing problems. The wholesale operation is not examined, and thus potential changes in the electricity generation mix are not validated in terms of retail market conditions. More specifically, there is little evidence in the literature on how increased capacity of renewably energy influences the cost of electricity and the RTP offered to the consumers that reflect this cost.
(ii)
The pool market prices are simulated as stochastic variables through a scenarios based approach. By incorporating stochastic programming simulations in mixed-integer linear problems, the overall complexity of the optimization problem increases.
(iii)
The topic of dynamic pricing has not received considerable attention. In the majority of the studies the selling price is fixed and not time-variant.
In the present study, the focus is to replace the stochastic programming approach by providing a market-clearing solution on the wholesale side. An integrated profit maximization model is proposed. The term integration refers to the connection between wholesale and retail markets. The proposed model includes a wholesale side and a retail side models. A market-clearing problem is solved to extract the SMP for various photovoltaics (PV) implementation scenarios. More specifically, the market-clearing problem is solved hourly for a full year. The scope of the scenarios is to assess how the implementation of different renewable generation capacity of an energy system influences the Retailer actions and the selling price limits. Based on the DR model that simulates the consumers’ behavior on the price, the final selling price is derived and sent to the consumers. It is assumed that the consumers react rationally and modify their load accordingly. Then, the profitability, i.e., the ratio of revenue to income is calculated. The proposed model is flexible; different wholesale market conditions can be simulated in detail; different procurement mechanisms can be exploited together with various types of dynamic pricing.

2. Methodology

2.1. General Framework

The schematic representation of the proposed model is shown in Figure 1. Apart from the generation mix, interconnections interchanges, and others, the outputs of the cost-optimal market-clearing problem include the hourly SMPs for a full year. The number of SMP time series equals the number of PV implementation scenarios. Two procurement mechanisms are regarded, namely the future market and pool market. The outputs of the profit maximization problem include the RTPs transferred to the consumers.

2.2. Mathematical Formulation of the Wholesale Market Optimization Model

This part deals with the cost-optimal market-clearing of a power system, considering the power units’ operational constraints and interconnections with neighboring power systems. The examined period is a daily one, with an hourly time step t . The power units taken into consideration include thermal power ones t h , nuclear n u , hydroelectric h , and renewable energy ones r . Thermal power units, along with nuclear power ones, are commonly referred to as type- l units. Additional sources of supply are the electricity imports m from neighboring power systems, while electricity exports x can also be carried to the same systems. Both supply and consumption entities offer/bid specific quantities at certain prices. Each offer/bid is also divided into a number of blocks with a non-decreasing (increasing) order for offers (bids). In particular:
  • price ( C l , b l , t u n i t )—quantity ( C B l , b l , t ) pair for the energy supply of type- l power units,
  • price ( C m , b l , t i m p )—quantity ( C B m , b l , t ) pair for the energy supply of electricity imports m ,
  • price ( C x , b l , t e x p )—quantity ( C B x , b l , t ) pair for the energy consumption of electricity exports x , and
  • price ( C d , b l , t d m )—quantity ( C B d , b l , t ) pair for the energy consumption of electricity demand bids d .
A specific cost ( C h , t h r ) is assigned with the cleared electricity supply from hydroelectric units, which is dependent on the assumed hydrological conditions of the specific date d t ( E d t h y d ), and the hourly availability of each hydroelectric unit h ( P h , t m a x ). Concerning the renewable energy units r , a certain cost ( C r , t h r ) is associated with the cleared electricity supply from renewable energy units, based on the assumed meteorological conditions for each type of resource r and time period t ( A V r , t ), and the installed capacity of each renewable energy source r ( P r m a x ).
Each energy supply unit l is characterized by a specific technical maximum ( P l m a x ) and minimum ( P l m i n ) for their operation. Specific limits exist in their ability to increase (up) or decrease (down) their outputs, commonly referred to as ramp-up ( R U l ) and ramp-down limits ( R D l ). Moreover, type- l units are identified based on their minimum uptimes ( T l U T ) and downtimes ( T l D T ).
When they provide secondary reserve, they should operate under Automatic Generation Control (AGC), in which they have other technical maximums ( P l m a x , A G C ) and minimums ( P l m i n , A G C ). Type- l units can also provide reserves to the systems that can be primary-up and -down, secondary-up and –down, and tertiary upward spinning and non-spinning. This reserve supply has a specific cost based on the type of reserve, the type of supplier l , and the time period t , namely, C l , t 1 u p , C l , t 1 d n , C l , t 2 u p , C l , t 2 d n , and C l , t 3 for primary-up, primary–down, secondary-up, secondary-down, and tertiary reserve. Further, there are upper bounds for the provision of each type of reserve per unit l , namely, R l 1 u p for primary-up, R l 1 d n for primary-down, R l 3 s for tertiary spinning, and R l 3 n s for tertiary non-spinning reserve. The maximum level for the provision of the secondary-up and –down reserve is related to the ramp-up ( R U l A G C ) and down limits ( R U l A G C ) of each unit l when operating under AGC. Two other cost components are associated with the start-up ( C l s t u ) and the shut-down ( C l s h d ) decision making of each type- l unit. Furthermore, the cases of unmet energy ( C t u e ) or inability to cover part of reserve requirements per type have been assigned with specific cost penalties ( C t u r ).
The system’s requirements for each reserve type include:
  • primary-up reserve requirements ( D R t 1 u p ),
  • primary-down reserve requirements ( D R t 1 d n ),
  • secondary-up reserve requirements ( D R t 2 u p ),
  • secondary-down reserve requirements ( D R t 2 d n ), and
  • tertiary reserve requirements ( D R t 3 ).

2.2.1. Objective Function of the Wholesale Market Optimization Problem

The mathematical model’s objective function refers to minimizing a given power system’s daily energy and reserves cost as provided by the following relationship (1).
M i n   C t o t a l = l b l t C l , b l , t u n i t · b l , b l , t T h e r m a l N u c l e a r + h r t C h r , t h r · p h r , t H y d r o R E S + m b l t C m , b l , t i m p · b m , b l , t I m p o r t s + l t C l s t u · y l , t S t a r t u p   c o s t + l t C l s h d · z l , t S h u t d o w n   c o s t + t n d t · C t u e E n e r g y   u n m e t   c o s t x b l t C x , b l , t e x p · b x , b l , t E x p o r t s   r e v e n u e s d b l t C d , b l , t d m · b d , b l , t D e m a n d   r e v e n u e s + l t C l , t 1 u p · r l , t 1 u p P r i m a r y u p   c o s t + l t C l , t 1 d n · r l , t 1 d n P r i m a r y d o w n   c o s t + l t C l , t 2 u p · r l , t 2 u p S e c o n d a r y u p   c o s t + l t C l , t 2 d n · r l , t 2 d n S e c o n d a r y d o w n   c o s t + l t C l , t 3 · r l , t 3 T e r t i a r y   c o s t + t C t u r · ( n v t 1 u p + n v t 1 d n + n v t 2 u p + n v t 2 d n + n v t 3 ) R e s e r v e s   u n m e t   c o s t

2.2.2. Energy Demand Balance

This section provides the mathematical set-up of the energy demand balance.
t h p t h , t + n u p n u , t h p h , t + r p r , t + m p m , t + n d t = d p d , t + x p x , t                 t
Equation (2) states that the total energy supply from all types of suppliers, including thermal t h , nuclear n u , hydroelectric h , renewable energy units r and imports m plus the amount of unmet demand must meet the cleared demand bids d and energy exports x in each time period t .

2.2.3. Operational Constraints

This set of constraints refers to operational limits. The available capacity of each energy supply entity l , m , and consumption entity x and d is divided into a number of blocks b l , each of which has a unique price-quantity pair.
b l , b l , t C B l , b l , t     l , b l , t    
b m , b l , t C B m , b l , t     m , b l , t
b x , b l , t C B x , b l , t       x , b l , t
b d , b l , t C B d , b l , t       d , b l , t      
p l , t = b l b l , b l , t         l , t
p m , t = b l b m , b l , t     m , t
p x , t = b l b x , b l , t       x , t    
p d , t = b l b d , b l , t     d , t
p l , t + r l , t 2 u p P l m a x · ( x l , t x l , t A G C ) + P l m a x , A G C · x l , t A G C   l , t
p l , t r l , t 2 d n P l m i n · ( x l , t x l , t A G C ) + P l m i n , A G C · x l , t A G C       l , t
p l , t + r l , t 1 u p + r l , t 2 u p + r l , t 3 s P l m a x · x l , t       l , t
p l , t r l , t 1 d n r l , t 2 d n P l m i n · x l , t                 l , t
The cleared energy supply of each supply entity l ( m ) in each block b l and time period t has the upper limit of the corresponding offered quantity, according to constraints (3) and (4). Constraints (5) and (6) defines that the cleared energy exports to each interconnected power system x (demand bids d ) in each block b l and time period t must not exceed the corresponding available quantity. Equations (7) and (8) set that the total cleared energy supply from each entity l ( m ) in each time period t amounts to the sum of all blocks of the respective cleared energy supply. The total cleared electricity exports to each interconnected power system x (demand bids d) in each interval t amount to the sum of all blocks of the corresponding cleared energy exports (demand bids), according to Equations (9) and (10). Constraints (11)–(14) model the operational limits, upper ((11) and (13)) and lower ((12) and (14)) respectively, of each unit l in each interval t , also considering their reserve provision capabilities (primary-up, primary-down, secondary-up and –down, and tertiary spinning).

2.2.4. Reserve Provision Limits

This section provides the mathematical set-up of the reserve provision limits.
r l , t 1 u p R l 1 u p · x l , t               l , t
r l , t 1 d n R l 1 d n · x l , t             l , t
r l , t 2 u p 15 · R U l A G C · x l , t A G C                 l , t  
r l , t 2 d n 15 · R D l A G C · x l , t A G C                 l , t  
r l , t 3 s R l 3 s · x l , t                                                       l , t
r l , t 3 n s R l 3 n s · x l , t 3 n s                                           l , t  
r l , t 3 n s P l m i n · x l , t 3 n s                                           l , t  
r l , t 3 n s R l 3 n s · ( 1 x l , t )                       l , t
p l , t p l , t 1 R U l · 60                         l , t
p l , t 1 p l , t R D l · 60                         l , t
This set of constraints (15)–(22) models the reserve provision limits of each reserve supplier l in each time period t , by imposing upper limits for each reserve type, subject to the entity’s dispatch ( x l , t ) or not ( 1 x l , t ), and operation under AGC or not ( x l , t A G C ). More specifically, constraint (15) puts an upper bound on the primary-up reserve supply, constraint (16) on the primary-down reserve supply, constraint (17) on the secondary-up reserve supply, constraint (18) on the secondary-down reserve supply, constraint (19) on the tertiary-spinning reserve supply, and constraints (20)–(22) on the tertiary non-spinning reserve supply. Constraints (23) and (24) define the ramp-up ( R U l ) and -down ( R D l ) limits, for each entity l in each interval t , namely the rate at which a supplying entity can increase (up) or decrease (down) its output.

2.2.5. System’s Reserve Requirements

This section provides the mathematical set-up of the system’s reserve requirements.
h t r h t , t 1 u p + n v t 1 u p D R t 1 u p               t
h t r h t , t 1 d n + n v t 1 d n D R t 1 d n           t
h t r h t , t 2 u p + n v t 2 u p D R t 2 u p             t  
h t r h t , t 2 d n + n v t 2 d n D R t 2 d n           t  
h t ( r h t , t 3 s + r h t , t 3 n s ) + n v t 3 = r h t , t 3 D R t 3           t  
The following group of constraints (25)–(29) defines the system-wide requirements for each reserve type, primary-up (25), primary-down (26), secondary-up (27), secondary-down (28), and tertiary spinning and non-spinning (29), respectively, which can be provided by hydrothermal units h t , namely hydroelectric and thermal power units plus the amount of unmet demand per reserve type.

2.2.6. Renewable Energy Generation

This section provides the mathematical set-up of renewable energy generation.
p r , t + n r r , t A V r , t · P r m a x           r , t  
Equation (30) imposes that the energy supply of each renewable energy technology r in each interval t ( p r , t ) plus the amount of the corresponding curtailed energy ( n r r , t ) has an upper bound related to the availability of each renewable energy technology r in each interval t ( A V r , t ) multiplied with its installed capacity in each interval t ( C C r , t ).

2.2.7. Hydroelectric Power Generation

This section provides the mathematical set-up of hydroelectric power generation.
h t d t p h , t E d t h y d               d t
p h , t P h , t m a x                                             h , t
Constraint (31) imposes that the daily energy supply from all hydropower units h in each date d t of the year ( h t d t p h , t ) must not exceed a specific maximum bound ( E d t h y d ). In addition, Constraint (32) states that the daily maximum must not exceed the total available installed capacity.

2.2.8. Time-Related Constraints

This section provides the mathematical set-up of the time-related constraints.
t = t T l U T + 1 t y l , t x l , t             l , t
t = t T l D T + 1 t z l , t 1 x l , t               l , t
y l , t z l , t = x l , t x l , t 1                 l , t  
x l , t A G C x l , t               l , t
Constraint (33) models the minimum uptime of each supply unit l , based on which it must be online in each interval t ( x l , t = 1 ) if and only if it has started-up ( y l , t = 1 ) during the previous ( T l U T 1 ) intervals. Analogously, Constraint (34) defines the minimum downtime of each supply unit l , according to which it must remain offline in each interval t ( x l , t = 0 ) if it has been down ( z l , t = 1 ) during the previous ( T l D T 1 ) intervals. Moreover, Equation (35) formulates the logical relationship between start-up ( y l , t ) and shut-down ( z l , t ) decisions of each supply entity l in each time period t . Furthermore, a supply entity l can operate under AGC if and only if it operates in the dispatch phase, according to Constraint (36).

2.2.9. Net Electricity Trading

This section provides the mathematical set-up of the net electricity trading.
p m , t p x , t   b l b m , b l , t               ( m , x )   M X , t
p x , t p m , t b l b x , b l , t                 ( m , x )   M X , t
Constraint (37) states that net electricity imports ( p m , t p x , t ) of each interconnection (imports) m in each time period t , also considering their reserve supply capability, must not exceed the corresponding interconnection capacity ( b l b m , b l , t ). Further, Constraint (38) sets the same condition for net electricity exports ( p x , t p m , t ).
The overall problem is formulated as a mixed-integer linear programming (MILP) problem, involving the cost minimization objective function (1), and subject to Constraints and Equations (2)–(38).

2.3. Mathematical Formulation of the Retail Market Optimization Model

The decision making actions of the Retailer are classified into medium-term and short-term actions. The medium-term actions determine the electricity purchased from the future market via forwards contracts (FCs). After this decision, the short-term actions define the amount purchased from the pool market and the selling price. Below the formulation of the retail market optimization model is analyzed.

2.3.1. Future Market

In order to form the FCs structure, it is necessary to categorize the 24 h period into various intra-day periods, according to the following:
  • Period#1={01:00,02:00,03:00,04:00,05:00,06:00,07:00,08:00,09:00,10:00,11:00,12:00,13:00,14:00}
  • Period#2={15:00,16:00,17:00,18:00}
  • Period#3={19:00, 20:00,21:00, 22:00}
  • Period#4={23:00,24:00}
Period#3 refers to peak hours while the remaining to off-peak hours. Based on these periods five FCs are considered, which differ intra-day wise. Each FC is composed of five blocks that correspond to each power amount that differ in the price. The block limits are 20, 40, 60, 80, and 100 kW. Let   f = 1 , 2 , , F be the number of FCs and j = 1 , 2 , , N j the number of blocks. The cost C t F that is related to the FCs in period t   is calculated by
C t F = f F j = 1 N j λ f j t F P f j t F , t
0 P f j t F P f j t , max F , f , j , t
P f t F = j = 1 N j P f t j F , f , t
where λ f j t F and P f j t F are the price and amount purchased from f -th contract and N j -th block, respectively, at period t = 1 , 2 , .. , T . Parameter P f j t , max F denotes the block size. Equation (40) sets a limit of the amount to be purchased of each block. The term refers to the maximum power offered by the specific FC, i.e., it refers to the block size. Equation (41) refers to the total power of each contract P f j t F which is given as the sum of the powers of each block.

2.3.2. Pool Market

After the solution of the market-clearing problem, the SMPs are available. Let λ t P be the price in the pool market at period   t = 1 , 2 , .. , T . The cost of energy purchased from the pool C t P is expressed as
C t P = λ t P E t P , t
where E t P is the energy purchased for the pool at period t .

2.3.3. Demand Response

It is considered that the Retailer serves one group of consumers of the same type. The responsiveness of the consumers to the selling price is modeled through a linear price/demand function. This means that when consumers are offered a high selling price, they linearly decrease their demand.
Initially, the consumers are charged with a nominal selling price r o at period t = 1 , 2 , .. , T . After the solution of the retail optimization problem, a new price is offered that maximizes Retailer’s profit. Let r ( t ) be a random selling price at period t = 1 , 2 , .. , T after the solution of the problem. Due to the nominal price r o the consumers’ nominal demand is indicated as d o . The demand response of the consumers to the selling price r ( t ) is denoted as d ( t ) . The linear price/demand function is expressed as
d ( r ( t ) ) = d o { 1 + β [ r ( t ) r o ] r o }
The parameter that controls the level of the demand variation is the elasticity β .

2.3.4. Demand Balance

To reach out into a feasible solution to the problem, an energy balance constraint should be imposed. For each time instant t the demand of the consumers E t D should be covered by two procurement mechanisms:
E t D = E t P + f F P f t F , t
To correspond to formulation of the problem, Equation (43) is rewritten as
E t D = E o t D { 1 + β ( λ t D λ o t D ) λ o t D } , t
where E o t D is the nominal demand, λ o t D is the nominal selling price and λ t D is the selling price at period t = 1 , 2 , .. , T .

2.3.5. Selling Price Limit

Theoretically, the higher the offered price the higher the profits will be. In order for the Retailer to become competitive and follow the rules of the market set by the regulatory authorities, a price limitation is needed. If the consumers elasticity would be extremely low, the selling price would correspond to unrealistic high values without price limitation. The price limitation is given by the following equation:
t = 1 Τ E t D λ t D λ max D τ = 1 Τ E t D
where λ max D is the upper limit of the selling price.

2.3.6. Expected Profit and Network Access Costs

The expected profit is obtained by the difference between the revenues and the costs. The revenues R t refer to the product between the selling price and the consumers’ load:
R t = λ t D E t D , t
Apart from the costs related with the electricity procurement, the retailer is charged with transmission and distribution system access tariffs C N :
C N = t ( λ e t E t D + λ p t P t )
where λ e t and λ p t are the energy and power elements of the network access tariff, respectively and P t Is the contracted power.
The Retailer’s profit R P is expressed as
R P = t ( R t C t P ) t C t F C N = t [ ( λ t D E t D λ t P Ε t P λ e t E t D λ p t P t ) f F j = 1 N j λ f j t F P f j t F ]
The profit function is composed of elements that refer to incomes and expenses.

2.3.7. Objective Function of the Retail Market Optimization Problem

The profit maximization problem is formulated as
Maximize P f j t F , λ t D , E t P t = 1 N T [ ( λ t D E t D λ t P Ε t P λ e t E t D λ p t P t ) f F j = 1 N j λ f j t F P f j t F ]
Subject to: (40), (44) and (46).

3. Results

Both the cost-optimal market-clearing and profit maximization problems are formulated as mixed-integer nonlinear problems and are solved by commercial software [32]. The proposed methodological framework for the wholesale market clearing problem has been assessed on an illustrative case study, the key characteristics of which are provided in Table 1 and Table 2. Table 1 summarizes the studied power system’s economic data, and Table 2, the key technical data. The examined power system includes four LIGnite-fired units (LIG), six Natural Gas Combined Cycle (NGCC) units, two Natural Gas-fired Gas Turbines (NGGT), and hydroelectric units. It additionally includes wind turbines and solar panels, each of which with a capacity of 1 GW. The available interconnection capacity for both imports and exports is assumed to be 2 GW.
The Hellenic energy system serves as the test system where both models are applied [33]. A general test case with a Retailer covering the demand of 84 residential consumers is regarded. The load series refer to the aggregated load of the consumers, and the network access tariffs correspond to the ones of the Hellenic system. The nominal electricity selling tariff r o is obtained by the sum of the SMP, transmission, and distribution network access tariffs increased by 10%. Following this approach, the Retailer charges the consumers an amount of 10% higher than the nominal, i.e., initial costs related to the pool market. This amount is considered in order to cover some of the operational expenses of the Retailer. The elasticity parameter is β = 1.50. Hence, it is considered that all consumers are very elastic in modifying the demand for small price changes. For simplicity reasons, it is considered that the elasticity is the same for all hours of the day and that all consumers are characterized by the same elasticity. It is assumed that the Retailer is a limited portfolio, i.e., its decisions in the retailer market do not influence retail and wholesale markets operation and conditions.
To fully examine the influence of various implementation rates of PV capacity in the wholesale market, eight scenarios are regarded, denoted as S s , s = 1 , 2 , , 8 . Table 3 shows the solutions of the market-clearing problem for a full year. It is shown the amount of annual PV installed capacity. Scenario S1 refers to 1000 MW while S8 to 8564 MW. Due to these different PV capacity mixes, there are visible variations in the capacities of the other technologies and the imports/exports balance. As the PVs shares in the electricity mix increases, thermal units are displaced. More specifically, there is a decrease in the generation of lignite and NG fueled units. Moreover, for 3000 MW generation and higher RES curtailments appear. RES technologies, apart from PVs, include onshore wind turbines and small hydroelectric plants. The last row of the Table presents the annual average SMP. It can be noticed that from S1 to S8, the mean SMP is decreased by 32.60%, i.e., from 49.17 €/MWh to 33.14 €/MWh. While expensive technologies are displaced, the cost of electricity decreases. It is assumed that the PV generation is compensated by a feed-in-tariff mechanism, i.e., all the generated amount is fed into the grid covering a specific part of the demand. The variations of the daily average SMPs time series are depicted in Figure 2. The figure presents the mean annual SMP per scenario.
To thoroughly investigate the influence of the various PV generation shares on the Retailer’s profits and the levels of electricity rates to the end-consumers, four cases are taken into account. The retail market optimization problem is solved for four different days in the year, one per different season, namely for 2 January 2019, 25 April 2019, 14 August 2019, and 9 October 2019. The first test day, i.e., 2 January 2019, is the next day to New Year Day. This means that it is close to an official bank holiday in Greece. Many professional activities are in halt during this day, and thus, it corresponds to low consumption, especially in the commercial and industrial sectors. This is also the case for 25 April 2019 and 14 August 2019. They refer to previous days of two official bank holidays in spring and summer, respectively. Day 9 October 2019 is a working day. Therefore, the selection of days allows the user to examine the profit maximization problem in working days and days with special conditions. Every scenario is examined for each test day, leading to a total number of 32 test cases. Figure 3 presents the SMP series for the various test cases as extracted by the execution of the cost-optimal market-coupling solution. It can be noticed that in periods with large PV generated capacity the SMP is zero. For instance, in 9 October 2019 and S8, the SMP is zero in 10:00–14:00 h. Because of the increased daylight in the noon hours, the PV generation is high leading to low generation costs. Thus, the nominal price r o is low in the specific period and includes only the network access costs increased by 10%.
Table 4 present the structure of the FCs. It should be noted that currently, there is no active FC in the Hellenic energy market. The number of blocks N j , block sizes P f j t , max F and prices λ f j t F are values set by the authors tailored to the conditions of the profit optimization problem.
The scope of the FC is to decrease the economic risks associated with the fluctuations of the SMP in the day-ahead market. It is a bilateral agreement between a potential producer and the Retailer to purchase future amounts of electricity at a fixed price. According to Table 4, each FC applies on four different intra-day periods. Prices among those periods differ. The higher prices are noticed at the period 19:00–22:00 h corresponding to the evening peak. As block size decreases, prices become higher. For instance, a producer that establishes a FC with the Retailer motivates the latter to buy high amounts of power. Additional amounts of the same block size is slightly more expensive. Each FC refers to different block sizes and the number of blocks is 5. The electricity purchase from one FC to another refers to higher prices.
The income, profit, and profitability for 2 January 2019, 25 April 2019, 14 August 2019, and 9 October 2019 are presented in Table 5, Table 6, Table 7 and Table 8, respectively. The highest profits are met in January. During winter, the PV production is minimal; thus, generation costs are increased. This fact leads the Retailer to offer high prices to the consumers. As the PV generation increases between the Scenarios, the profits reduce. Contrary to the decreased generation costs, and thus lower prices, August refers to increased profits. This is due to the fact that in the summer months, the annual demand peaks of Hellenic system are observed. There are no great changes in the profitability between the scenarios of the same test day. The profitability provides information regarding the relationship between the obtained revenues and real profits. While the upper threshold of selling price to the consumers is limited from market regulation to keep the retail market competition robust, the Retailer, in order to increase the profitability ratio, needs to establish more favorable FCs with the producer.
Figure 4, Figure 5, Figure 6 and Figure 7 illustrate the nominal price and RTPs per scenario for 2 January 2019, 25 April 2019, 14 August 2019, and 9 October 2019, respectively.
According to the results, both the initial price, i.e., the initial RTPs and final RTPs obtained by the solution of the profit maximization problem, follow the trends of the hourly generation cost as expressed by the SMP. Considering a linear relationship between the offered price and the demand response, an increment of price increase results in a demand decrease. The final demand is lower than the nominal one in all hours. Figure 8, Figure 9, Figure 10 and Figure 11 show the initial and final demand demands, i.e., demand responses for 2 January 2019, 25 April 2019, 14 August 2019, and 9 October 2019, respectively.
The RTP scheme leads to total energy demand reduction. In a theoretical case, where the Retailer is obliged by the system operator to reduce the demand of its clients, a RTP scheme can be applied. Additionally, it is assumed that the initial demand is known. For instance, this information can be obtained by a short-term load forecasting routine.
As mentioned before, it is assumed that the same elasticity value characterizes all consumers. According to the figures, the selling prices series follow the nominal ones. The observed deviations from the nominal price depend on the procurement mechanism and the selling price limitation. While the Retailer does not influence the wholesale market price determination and the upper limit of the selling price limitation in the retail market, the only parameter for open negotiation is the FCs. The selling price limitation is usually set by the regulation authority. The deviations between the nominal and final demand are analogous to the deviations of prices.
The linear price/demand function is expressed as a straight line with a negative slope. The elasticity expresses the slope. In order to examine the influence of the elasticity variation in the profit, different values are considered for a sensitivity analysis. The elasticity values under consideration are β = { 0.50 , 0.75 , 1.00 , 1.50 } . For the analysis, test day 2 January 2019 and S1 are used. This test case refers to the highest profit and profitability compared to the others. Table 9 presents the respective results. Values near zero correspond to high inelastic behavior. It can be noticed that the inelastic demand pattern corresponds to profit and profitability increments. As consumers become more elastic, they take advantage of the low electricity prices and adjust their demand accordingly.
Figure 12 shows the RTPs per elasticity value. For β = 0.50 and β = 0.75 selling prices are very high. This fact may lead to the consumers’ dissatisfaction with the pricing policy. Figure 13 present the initial and final demands per elasticity value. For extreme selling prices, large reductions are expected.
Table 10 presents an example of the contribution of the two procurement mechanisms for the 24 h period. Test day 2 January 2019 is regarded and scenarios S1 and S8. The scope is to compare the deviations in the demand coverage between the two extreme scenarios. During early morning hours the cost of electricity is low and thus, all demand is covered by the pool market. The two scenarios differ in the period 11:00–16:00 h. Since this period usually corresponds to high solar irradiation and hence, PV generation, the SMP is lower. In the case of scenario S8 all electricity is purchased from the pool market.

4. Conclusions

This paper presented a methodology that is built in two optimization models. The methodology aims at integrating the two competitive markets; the scope is to study the influence of the wholesale market conditions in the retail market, both for the Retailer and consumers. In the literature, the SMP is obtained through a set of scenarios. In this paper, the SMP is obtained by the solution of a market-clearing problem. Therefore it is based on actual data in the generation side and not extracted via a simulation based approach. Moreover, the CVaR method is not necessary for modeling stochastic processes in the proposed methodology.
The main conclusions drawn by the study can be summarized in the following: the different PV implementation shares on the generation side highly influences the profits and the levels of the selling prices. An increased PV capacity results in lower generation costs due to the fact that expensive marginal units are displaced. This eventually leads to lower profits since the RTPs are directly connected with the SMPs. Lower selling prices can be witnessed not only between the seasons but also between the periods of the day. When daily PV generation is high, the SMP can be zero. To overcome this effect, it is proposed to form a basic fixed selling price is related to network access costs. This approach prevents the Retailer from having negative profits. Moreover, the higher the deviation is between SMPS and RTP, the higher profits can be accomplished. This finding is based only on a predefined value of the elasticity parameter. When the demand becomes inelastic, the prices are increased together with profits. This is due to the fact that the Retailer needs to increase the prices in cases of low demand so that to obtain profits.
The contribution of the paper can be summarized in the following:
  • The methodology of the paper unifies the wholesale and retail market operations through two optimization models. The first one refers to a market-learning price problem where the objective is to minimize electricity generation cost. The decision variables of the problem are the schedules of the power plants and interconnection exchanges. The retail market model aims at the maximization of the Retailer’s profit. Here the decision variables are the procurement mechanism and the real-time selling prices to the consumers.
  • The unified consideration of the two markets minimizes the Retailer’s economic risks and allows more robust and reliable strategic decisions of the Retailer in the competitive retail market. Wholesale market prices are not modeled as stochastic variables but are a product of the solution of the wholesale market operation.
  • Both models are flexible in terms of further expansions and additions. In particular, the retail market model can take into consideration different price/demand function, elasticity values, and electricity tariff structures.
  • The methodology of the paper quantifies and evaluates how different PV penetration levels affect the selling prices to the consumers. RTPs are implemented in the model in order to connect wholesale market conditions to the retail market. The actual electricity costs are transferred to the consumers. By modeling the consumer’s response to the selling prices, load management techniques can be manifested, leading to benefits on the retail side, e.g., lower electricity costs for consumers, and on the wholesale side, e.g., congestion management in the transmission and distribution levels and others.
Future expansions of the methodology will regard different types of price/demand functions and different schemes of load management. Through the RTPs developed in the paper, total demand reduction is feasible. The set-up of the retail market problem can be modified to include other schemes like peak clipping or valley filling. Finally, other RES technologies implementation scenarios will be investigated.

Author Contributions

I.P.P. performed the research, developed the retail market optimization model and wrote the paper. N.K. provided guidance for the research and developed the wholesale market optimization model. G.C.C. provided guidance for the research and revised the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Sets
l L Set of thermal and nuclear power units
h H Set of hydroelectric units
h r H R Set of hydroelectric and renewable energy units
h t H T Set of hydrothermal units
b l B L Set of blocks of energy supply and consumption functions
d D Set of demand entities
m M Set of interconnection for energy imports
n u N U Set of nuclear units
r R Set of renewable energy sources
t T Set of time periods
t h T H Set of thermal units
x X Set of interconnection for energy exports
Parameters
C l , b l , t u n i t Cost of each block b l of each type- l unit’s energy supply function in period t [€/MW]
A V r , t Availability factor of unit r in period t [pu]
C h r , t h r Cost of each unit’s h r energy supply offer in period t [€/MW]
C B d , b l , t Available quantity of block b l of each demand entity’s d energy consumption function in period t [MW]
C B l , b l , t Available quantity of block b l of each type- l unit’s energy supply function in period t [MW]
C B m , b l , t Available quantity of block b l of each interconnection’s m energy supply function in period t [MW]
C B x , b l , t Available quantity of block b l of each interconnection’s x energy consumption function in period t [MW]
C d , b l , t d m Cost of each block b l of each demand entity’s d energy consumption function in period t [€/MW]
C l , t 1 d n Primary-down reserve supply cost of type-l unit in period t [€/MW]
C l , t 1 u p Primary-up reserve supply cost of type-l unit in period t [€/MW]
C l , t 2 d n Secondary -down reserve supply cost of type-l unit in period t [€/MW]
C l , t 2 u p Secondary-up reserve supply cost of type-l unit in period t [€/MW]
C l , t 3 Tertiary reserve supply cost of type-l unit in period t [€/MW]
C l s h d Shut-down cost of type- l unit [€/shut-down]
C l s t u Start-up cost of type- l unit [€/start-up]
C m , b l , t i m p Cost of each block b l of each interconnection’s m energy supply function in period t [€/MW]
C t u e Unmet energy cost in period t [€/MW]
C t u r Unmet reserves cost in period t [€/MW]
C x , b l , t e x p Cost of each block b l of each interconnection’s x energy consumption function in period t [€/MW]
D R t 1 d n Primary-down reserve requirements in period t [MW]
D R t 1 u p Primary-up reserve requirements in period t [MW]
D R t 2 d n Secondary-down reserve requirements in period t [MW]
D R t 2 u p Secondary-up reserve requirements in period t [MW]
D R t 3 Tertiary reserve requirements in period t [MW]
E d t h y d Maximum daily amount of hydroelectric generation in date d t [MWh]
P h , t m a x Available maximum capacity of unit h in period t [MW]
P l m a x , A G C Technical maximum of type- l unit under automatic generation control [MW]
P l m a x Technical maximum of type- l unit [MW]
P l m i n , A G C Technical minimum of type- l unit under automatic generation control [MW]
P l m i n Technical minimum of type- l unit [MW]
P r m a x Installed capacity of unit r [MW]
R D l Ramp-down limit of type- l unit in period t [MW/min]
R D l A G C Ramp-down limit of type- l unit under automatic generation control in period t [MW/min]
R U l Ramp-up limit of type- l unit in period t [MW/min]
R U l A G C Ramp-up limit of type- l unit under automatic generation control in period t [MW/min]
R l 1 d n Maximum primary-down reserve supply of type- l unit [MW]
R l 1 u p Maximum primary-up reserve supply of type- l unit [MW]
R l 3 n s Maximum tertiary non-spinning reserve supply of type- l unit [MW]
R l 3 s Maximum tertiary spinning reserve supply of type- l unit [MW]
T l D T Minimum downtime of type- l unit [h]
T l U T Minimum uptime of type- l unit [h]
Continuous variables
b l , b l , t Cleared quantity of block b l of each type- l unit’s energy supply function in period t (MW)
b d , b l , t Cleared quantity of block b l of each demand entity’s d energy consumption function in period t (MW)
b m , b l , t Cleared quantity of block b l of each interconnection’s m energy supply function in period t (MW)
b x , b l , t Cleared quantity of block b l of each interconnection’s x energy consumption function in period t (MW)
n d t Cleared quantity of unmet energy demand in period t [MW]
n r r , t Cleared quantity of the curtailed energy of unit r in time period t [MW]
n v t 1 d n Cleared quantity of unmet primary-down reserve requirements in period t [MW]
n v t 1 u p Cleared quantity of unmet primary-up reserve requirements in period t [MW]
n v t 2 d n Cleared quantity of unmet secondary-down reserve requirements in period t [MW]
n v t 2 u p Cleared quantity of unmet secondary-up reserve requirements in period t [MW]
n v t 3 Cleared quantity of unmet tertiary reserve requirements in period t [MW]
p h , t Cleared power output of unit h in time period t [MW]
p h r , t Cleared power output of unit h r in time period t [MW]
p d , t Cleared power consumption of demand entity d in time period t [MW]
p m , t Cleared power output of interconnection m in time period t [MW]
p n u , t Cleared power output of unit n u in time period t [MW]
p r , t Cleared power output of unit r in time period t [MW]
p t h , t Cleared power output of unit t h in time period t [MW]
p x , t Cleared power consumption of interconnection x in time period t [MW]
r l , t 1 d n Cleared primary-down reserve provision of type- l unit in period t [MW]
r l , t 1 u p Cleared primary-up reserve provision of type- l unit in period t [MW]
r l , t 2 d n Cleared secondary-down reserve provision of type- l unit in period t [MW]
r l , t 2 u p Cleared secondary-up reserve provision of type- l unit in period t [MW]
r l , t 3 Cleared tertiary reserve provision of type- l unit in period t [MW]
r l , t 3 n s Cleared tertiary non-spinning reserve provision of type- l unit in period t [MW]
r l , t 3 s Cleared tertiary spinning reserve provision of type- l unit in period t [MW]
Binary variables
x l , t 1, if type- l unit operates in dispatch phase in period t
0, otherwise
x l , t 3 n s 1, if type- l unit provides non-spinning tertiary reserve in period t
0, otherwise
x l , t A G C 1, if type- l unit operates under automatic generation control in period t
0, otherwise
y l , t 1, if type- l unit starts-up in period t
0, otherwise
z l , t 1, if type- l unit shuts-down in period t
0, otherwise

References

  1. Siksnelyte, I.; Zavadskas, E.M. Achievements of the European Union countries in seeking a sustainable electricity sector. Energies 2019, 12, 2254. [Google Scholar] [CrossRef] [Green Version]
  2. Hyland, M. Restructuring European electricity markets—A panel data analysis. Util. Pol. 2016, 30, 33–42. [Google Scholar] [CrossRef] [Green Version]
  3. Newbery, D.; Strbac, G.; Viehoff, I. The benefits of integrating European electricity markets. Energy Pol. 2016, 94, 253–263. [Google Scholar] [CrossRef] [Green Version]
  4. Ringler, P.; Keles, D.; Fichtner, W. How to benefit from a common European electricity market design. Energy Pol. 2017, 101, 629–643. [Google Scholar] [CrossRef]
  5. Sencar, M.; Pozeb, V.; Krope, T. Development of EU (European Union) energy market agenda and security of supply. Energy 2014, 77, 117–124. [Google Scholar] [CrossRef]
  6. Correljé, A.F.; De Vries, L.J. Competitive Electricity Markets Design, Implementation, Performance. Elsevier Global Energy Policy and Economics Series, 1st ed.; Elsevier Science: Amsterdam, The Netherlands, 2008. [Google Scholar]
  7. Sioshansi, F.; Pfaffenberger, W. Electricity Market Reform: An International Perspective, 1st ed.; Elsevier Science: Amsterdam, The Netherlands, 2006. [Google Scholar]
  8. Bae, M.; Kim, H.; Kim, E.; Chung, A.Y.; Kim, H.; Roh, J.H. Toward electricity retail competition: Survey and case study on technical infrastructure for advanced electricity market system. Appl. Energy 2014, 133, 252–273. [Google Scholar] [CrossRef]
  9. Boroumand, R.H.; Goutte, S.; Porcher, S.; Porcher, T. Hedging strategies in energy markets: The case of electricity retailers. Energy Econ. 2015, 51, 503–509. [Google Scholar] [CrossRef]
  10. Guo, Y.; Shao, P.; Wang, J.; Dou, X.; Zhao, W. Purchase strategies for power retailers considering load deviation and CVaR. Energy Proc. 2019, 158, 6658–6663. [Google Scholar] [CrossRef]
  11. Gabriel, S.A.; Genc, M.F.; Balakrishnan, S. A simulation approach to balancing annual risk and reward in retail electrical power markets. IEEE Trans. Power Syst. 2002, 17, 1050–1057. [Google Scholar] [CrossRef]
  12. Gabriel, S.A.; Conejo, A.J.; Plazas, M.A.; Balakrishnan, S. Optimal price and quantity determination for retail electric power contracts. IEEE Trans. Power Syst. 2006, 21, 180–187. [Google Scholar] [CrossRef]
  13. Carrión, M.; Conejo, A.J.; Arroyo, J.M. Forward contracting and selling price determination for a retailer. IEEE Trans. Power Syst. 2007, 22, 2105–2114. [Google Scholar] [CrossRef]
  14. Yusta, J.M.; Rosado, I.J.R.; Navarro, J.A.D.; Vidal, J.M.P. Optimal electricity price calculation model for retailers in a deregulated market. Int. J. Electr. Power Energy Syst. 2005, 27, 437–447. [Google Scholar] [CrossRef]
  15. Carrión, M.; Arroyo, J.M.; Conejo, A.J. A bilevel stochastic programming approach for retailer futures market trading. IEEE Trans. Power Syst. 2009, 24, 1446–1456. [Google Scholar] [CrossRef]
  16. Hatami, R.A.; Seifi, H.; Sheikh-El-Eslami, M.K. Optimal selling price and energy procurement strategies for a retailer in an electricity market. Electr. Power Syst. Res. 2009, 79, 246–254. [Google Scholar] [CrossRef]
  17. Sun, B.; Wang, F.; Xie, J.; Sun, X. Electricity Retailer trading portfolio optimization considering risk assessment in Chinese electricity market. Electr. Power Syst. Res. 2021, 190, 1–12. [Google Scholar] [CrossRef]
  18. Chawda, S.; Mathuria, P.; Bhakar, R. Risk-based retailer profit maximization: Time of use price setting for elastic demand. Int. Trans. Electr. Energ. Syst. 2019, 29, 1–22. [Google Scholar] [CrossRef]
  19. Sekizaki, S.; Nishizaki, I.; Hayashida, T. 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. 2019, 105, 877–893. [Google Scholar] [CrossRef]
  20. Charwand, M.; Gitizadeh, M.; Siano, P. A new active portfolio risk management for an electricity retailer based on a drawdown risk preference. Energy 2017, 118, 387–398. [Google Scholar] [CrossRef]
  21. Boroumanda, R.H.; Goutteb, S.; Guesmid, K.; Porchera, T. Potential benefits of optimal intra-day electricity hedging for the environment: The perspective of electricity retailers. Energy Pol. 2019, 132, 1120–1129. [Google Scholar] [CrossRef] [Green Version]
  22. do Prado, J.C.; Qiao, W. A stochastic decision-making model for an electricity retailer with intermittent renewable energy and short-term demand response. IEEE Trans. Smart Grid 2019, 10, 2581–2592. [Google Scholar] [CrossRef]
  23. Fateh, H.; Safari, A.; Bahramara, S. A bi-level optimization approach for optimal operation of distribution networks with retailers and micro-grids. J. Oper. Autom. Power Eng. 2020, 8, 15–21. [Google Scholar]
  24. Algarvio, H.; Lopes, F.; Sousa, J.; Lagarto, J. Multi-agent electricity markets: Retailer portfolio optimization using Markowitz theory. Electr. Power Syst. Res. 2017, 148, 282–294. [Google Scholar] [CrossRef]
  25. Nojavan, S.; Zare, K. Optimal energy pricing for consumers by electricity retailer. Int. J. Electr. Power Energy Syst. 2018, 102, 401–412. [Google Scholar] [CrossRef]
  26. Charwand, M.; Gitizadeh, M. Risk-based procurement strategy for electricity retailers: Different scenario-based methods. IEEE Trans. Eng. Manag. 2020, 67, 1–11. [Google Scholar] [CrossRef]
  27. Qadrdan, M.; Cheng, M.; Wu, J.; Jenkins, N. Benefits of demand-side response in combined gas and electricity networks. Appl. Energy 2017, 19215, 360–369. [Google Scholar] [CrossRef] [Green Version]
  28. Alipour, M.; Zare, K.; Seyedi, H.; Jalali, M. Real-time price-based demand response model for combined heat and power systems. Energy 2019, 1681, 1119–1127. [Google Scholar] [CrossRef]
  29. Yusta, J.M.; Khodr, H.M.; Urdaneta, A.J. Optimal pricing of default customers in electrical distribution systems: Effect behavior performance of demand response models. Electr. Power Syst. Res. 2007, 77, 548–558. [Google Scholar] [CrossRef]
  30. Mahmoudi-Kohan, N.; Parsa Moghaddam, M.; Sheikh-El-Eslami, M.K. An annual framework for clustering-based pricing for an electricity retailer. Electr. Power Syst. Res. 2010, 80, 1042–1048. [Google Scholar] [CrossRef]
  31. Mahmoudi-Kohan, N.; Parsa Moghaddam, M.; Sheikh-El-Eslami, M.K.; Shayesteh, E. A three-stage strategy for optimal price offering by a retailer based on clustering techniques. Int. J. Electr. Power Energy Syst. 2010, 32, 1135–1142. [Google Scholar] [CrossRef]
  32. Brooke, A.; Kendrick, D.; Meeraus, A.; Raman, R. GAMS: A User’s Guide; GAMS Development Corporation: Washington, DC, USA, 1998. [Google Scholar]
  33. Energy Exchange Group, S.A. Available online: http://www.enexgroup.gr/nc/en/home/ (accessed on 25 December 2020).
Figure 1. Visual description of the proposed model. (source: Authors).
Figure 1. Visual description of the proposed model. (source: Authors).
Energies 14 00092 g001
Figure 2. Daily average system marginal price (SMP) series per scenario. (source: Authors).
Figure 2. Daily average system marginal price (SMP) series per scenario. (source: Authors).
Energies 14 00092 g002
Figure 3. Daily average SMP per test day: (a) 2 January 2019, (b) 25 April 2019, (c) 14 August 2019, and (d) 9 October 2019. (source: Authors).
Figure 3. Daily average SMP per test day: (a) 2 January 2019, (b) 25 April 2019, (c) 14 August 2019, and (d) 9 October 2019. (source: Authors).
Energies 14 00092 g003
Figure 4. Nominal price and RTP corresponding to day 2 January 2019 for: (a) S1, (b) S2, (c) S3, (d) S4, (e) S5, (f) S6, (g) S7, and (h) S8. (source: Authors).
Figure 4. Nominal price and RTP corresponding to day 2 January 2019 for: (a) S1, (b) S2, (c) S3, (d) S4, (e) S5, (f) S6, (g) S7, and (h) S8. (source: Authors).
Energies 14 00092 g004
Figure 5. Nominal price and real-time pricing (RTP) corresponding to day 25 April 2019 for: (a) S1, (b) S2, (c) S3, (d) S4, (e) S5, (f) S6, (g) S7, and (h) S8. (source: Authors).
Figure 5. Nominal price and real-time pricing (RTP) corresponding to day 25 April 2019 for: (a) S1, (b) S2, (c) S3, (d) S4, (e) S5, (f) S6, (g) S7, and (h) S8. (source: Authors).
Energies 14 00092 g005
Figure 6. Nominal price and RTP corresponding to day 16 August 2019 for: (a) S1, (b) S2, (c) S3, (d) S4, (e) S5, (f) S6, (g) S7, and (h) S8. (source: Authors).
Figure 6. Nominal price and RTP corresponding to day 16 August 2019 for: (a) S1, (b) S2, (c) S3, (d) S4, (e) S5, (f) S6, (g) S7, and (h) S8. (source: Authors).
Energies 14 00092 g006
Figure 7. Nominal price and RTP corresponding to day 9 October 2019 for: (a) S1, (b) S2, (c) S3, (d) S4, (e) S5, (f) S6, (g) S7, and (h) S8. (source: Authors).
Figure 7. Nominal price and RTP corresponding to day 9 October 2019 for: (a) S1, (b) S2, (c) S3, (d) S4, (e) S5, (f) S6, (g) S7, and (h) S8. (source: Authors).
Energies 14 00092 g007
Figure 8. Nominal demand and demand per scenario for 2 January 2019. (source: Authors).
Figure 8. Nominal demand and demand per scenario for 2 January 2019. (source: Authors).
Energies 14 00092 g008
Figure 9. Nominal demand and demand per scenario for 25 April 2019. (source: Authors).
Figure 9. Nominal demand and demand per scenario for 25 April 2019. (source: Authors).
Energies 14 00092 g009
Figure 10. Nominal demand and demand per scenario for 14 August 2019. (source: Authors).
Figure 10. Nominal demand and demand per scenario for 14 August 2019. (source: Authors).
Energies 14 00092 g010
Figure 11. Nominal demand and demand per scenario for 9 October 2019. (source: Authors).
Figure 11. Nominal demand and demand per scenario for 9 October 2019. (source: Authors).
Energies 14 00092 g011
Figure 12. Nominal price and RTPs per elasticity values. (source: Authors).
Figure 12. Nominal price and RTPs per elasticity values. (source: Authors).
Energies 14 00092 g012
Figure 13. Nominal demand and final demand per elasticity values. (source: Authors).
Figure 13. Nominal demand and final demand per elasticity values. (source: Authors).
Energies 14 00092 g013
Table 1. Main economic data of the studied power system.
Table 1. Main economic data of the studied power system.
UnitStart-Up Cost
(€/Start-Up)
Shut-Down Cost
(€/Shut-Down)
CO2 Emission Factor
(tnCO2/MWh)
CO2 Price
(€/tnCO2)
Fuel Price
(€/t for Lignite and €/MWhth for Natural Gas)
Fuel Heating Value (GJ/t for Lignite)Efficiency
(p.u.)
O&M Cost
(€/MWh)
Minimum Average Variable Cost
(€/MWh)
LIG-150,00020,0001.52519.564.880.382.6578.50
LIG-250,00020,0001.4257.504.290.331.5055.46
LIG-350,00010,0001.32521.007.330.382.1462.10
LIG-450,00010,0001.22515.774.770.331.8567.43
NGCC-120,00010,0000.372515.350.000.581.2037.04
NGCC-220,00010,0000.372515.350.000.571.2037.58
NGCC-320,00010,0000.372515.350.000.581.2037.08
NGCC-420,00010,0000.372515.350.000.601.2036.03
NGCC-520,00010,0000.372515.350.000.581.2036.99
NGCC-620,00010,0000.372515.350.000.511.2040.76
NGGT-1250010000.532530.000.000.401.2089.45
NGGT-2250010000.532530.000.000.401.2089.45
HYDRO000250.000.000.000.003.00
Table 2. Main technical data of the studied power system.
Table 2. Main technical data of the studied power system.
UnitTechnical Maximum (MW)Technical Minimum (MW)Technical Maximum under AGC (MW)Technical Minimum under AGC (MW)Primary-Up Reserve Capability (MW)Ramp-Up Limit under AGC (MW/min)Ramp-Down Limit under AGC (MW/min)Tertiary Spinning Reserve Capability (MW)Tertiary Non-Spinning Reserve Capability (MW)Ramp-Up Limit (MW/min)Ramp-Down Limit (MW/min)Minimum Uptime (h)Minimum Downtime (h)
LIG-134218800280045044161
LIG-228915100280045044161
LIG-325619500280045044161
LIG-427315000280045044161
NGCC-13782203602403612123403401212123
NGCC-24221953902203612123803801212123
NGCC-33902203702403610.510.535135110.510.5123
NGCC-44331953902203614143903901414123
NGCC-54171823902003619.219.23753752424123
NGCC-655094550943612124954951212123
NGGT-115050120638771201208811
NGGT-216055128698771251258811
HYDRO1000010002501501501000100015015000
Table 3. Solutions of the market-clearing problem.
Table 3. Solutions of the market-clearing problem.
ScenarioS1S2S3S4S5S6S7S8
PV (MW)10002000300040005000600070008564
Lignite (GWh)1042.02975.23764.25521.05391.56302.43341.851751.242
Natural gas (GWh)13,944.4114,163.3413,923.1413,583.9013,425.5413,060.1613,537.8112,970.15
Hydro (GWh)1026.921026.921026.921026.921026.921026.921026.921026.92
RES (GWh)3748.205183.426614.807981.518907.639427.459294.219029.44
Imports (GWh)1523.43983.0232811.01525.06314.15294.04180.33141.94
Exports (GWh)−5652.59−6699.54−7507.73−8006.06−8433.40−8478.60−8748.74−9287.30
Demand (GWh)15,632.4115,632.4115,632.4115,632.4115,632.4115,632.4115,632.4115,632.41
RES curtailment (GWh)003.8472.35581.451496.853065.315576.03
Unmet energy (GWh)00000000
Marginal Price (€/MWh)49.1748.4046.8444.0240.1737.6435.4433.14
Table 4. Structure of the forward contracts (FCs).
Table 4. Structure of the forward contracts (FCs).
FCsBlocksBlock Size (kW) 01:00–14:00 h
(€/MWh)
15:00–18:00 h
(€/MWh)
19:00–22:00 h
(€/MWh)
23:00–00:00 h
(€/MWh)
1110039.1337.6251.6242.23
1210043.0541.3856.7846.45
1310047.3545.5262.4651.09
1410052.0950.0768.7056.20
1510057.3055.0875.5761.82
218039.7938.2552.4842.94
228043.7742.0857.7347.23
238048.1546.2963.5151.95
248052.9650.9169.8657.15
258058.2656.0176.8462.86
316040.4638.9053.3743.66
326044.5142.7858.7048.02
336048.9647.0664.5752.83
346053.8551.7771.0358.11
356059.2456.9578.1363.92
414041.1439.5554.2644.39
424045.2543.5059.6948.83
434049.7847.8565.6653.71
444054.7652.6472.2359.09
454060.2357.9079.4564.99
512041.8340.2155.1845.14
522046.0244.2360.6949.65
532050.6248.6666.7654.62
542055.6853.5273.4460.08
552061.2558.8880.7866.09
Table 5. Income, profit, and profitability for the various scenarios for 2 January 2019.
Table 5. Income, profit, and profitability for the various scenarios for 2 January 2019.
ScenarioIncome (€)Profit (€)Profitability
#11193.892641.9190.538
#21156.077609.6690.527
#3984.374477.0380.485
#4964.551468.4240.486
#5888.334440.6580.496
#6886.789435.8210.491
#7798.472387.7620.486
#8785.271380.0890.484
Table 6. Income, profit, and profitability for the various scenarios for 25 April 2019.
Table 6. Income, profit, and profitability for the various scenarios for 25 April 2019.
ScenarioIncome (€)Profit (€)Profitability
#1614.113205.1500.334
#2612.801204.3990.334
#3608.985202.6050.333
#4570.626189.9220.333
#5497.872164.5630.333
#6477.856158.0670.331
#7446.772148.9190.333
#8403.441132.4270.328
Table 7. Income, profit, and profitability for the various scenarios for 14 August 2019.
Table 7. Income, profit, and profitability for the various scenarios for 14 August 2019.
ScenarioIncome (€)Profit (€)Profitability
#1878.211380.9580.434
#2876.408379.5720.433
#3842.977355.6100.422
#4755.901301.6040.399
#5648.294244.4340.377
#6595.642228.0650.383
#7535.919193.4530.361
#8502.484182.8250.364
Table 8. Income, profit, and profitability for the various scenarios for 9 October 2019.
Table 8. Income, profit, and profitability for the various scenarios for 9 October 2019.
ScenarioIncome (€)Profit (€)Profitability
#1823.845346.5280.421
#2772.802312.8140.405
#3728.675284.3310.390
#4680.723270.5030.397
#5618.500241.2960.390
#6570.156219.8970.386
#7546.718210.7900.386
#8530.865325.5150.387
Table 9. Income, profit, and profitability for the various elasticity values.
Table 9. Income, profit, and profitability for the various elasticity values.
ElasticityIncome (€)Profit (€)Profitability
−1.501193.892641.9190.538
−1.001192.364711.8430.597
−0.751244.569800.3770.643
−0.501404.487996.4120.709
Table 10. Contribution of the two procurement mechanisms for S1 and S8.
Table 10. Contribution of the two procurement mechanisms for S1 and S8.
Hour (h)FCs (%)Pool Market (%)Total Energy (%)Hour (h)FCs (%)Pool Market (%)Total Energy (%)
10 (S1) 0 (S8)100 (S1) 100 (S8)10013100 (S1) 0 (S8)0 (S1) 100 (S8)100
20 (S1) 0 (S8)100 (S1) 100 (S8)10014100 (S1) 0 (S8)0 (S1) 100 (S8)100
30 (S1) 0 (S8)100 (S1) 100 (S8)10015100 (S1) 0 (S8)0 (S1) 100 (S8)100
40 (S1) 0 (S8)100 (S1) 100 (S8)10016100 (S1) 0 (S8)0 (S1) 100 (S8)100
50 (S1) 0 (S8)100 (S1) 100 (S8)100170 (S1) 0 (S8)100 (S1) 100 (S8)100
60 (S1) 0 (S8)100 (S1) 100 (S8)100180 (S1) 0 (S8)100 (S1) 100 (S8)100
7100 (S1) 100 (S8)0 (S1) 0(S8)100190 (S1) 0 (S8)100 (S1) 100 (S8)100
8100 (S1) 100 (S8)0 (S1) 0(S8)100200 (S1) 0 (S8)100 (S1) 100 (S8)100
9100 (S1) 100 (S8)0 (S1) 0(S8)100210 (S1) 0 (S8)100 (S1) 100 (S8)100
100 (S1) 0 (S8)100 (S1) 100 (S8)100220 (S1) 0 (S8)100 (S1) 100 (S8)100
11100 (S1) 0 (S8)0 (S1) 100 (S8)100230 (S1) 0 (S8)100 (S1) 100 (S8)100
12100 (S1) 0 (S8)0 (S1) 100 (S8)100240 (S1) 0 (S8)100 (S1) 100 (S8)100
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Panapakidis, I.P.; Koltsaklis, N.; Christoforidis, G.C. A Novel Integrated Profit Maximization Model for Retailers under Varied Penetration Levels of Photovoltaic Systems. Energies 2021, 14, 92. https://doi.org/10.3390/en14010092

AMA Style

Panapakidis IP, Koltsaklis N, Christoforidis GC. A Novel Integrated Profit Maximization Model for Retailers under Varied Penetration Levels of Photovoltaic Systems. Energies. 2021; 14(1):92. https://doi.org/10.3390/en14010092

Chicago/Turabian Style

Panapakidis, Ioannis P., Nikolaos Koltsaklis, and Georgios C. Christoforidis. 2021. "A Novel Integrated Profit Maximization Model for Retailers under Varied Penetration Levels of Photovoltaic Systems" Energies 14, no. 1: 92. https://doi.org/10.3390/en14010092

APA Style

Panapakidis, I. P., Koltsaklis, N., & Christoforidis, G. C. (2021). A Novel Integrated Profit Maximization Model for Retailers under Varied Penetration Levels of Photovoltaic Systems. Energies, 14(1), 92. https://doi.org/10.3390/en14010092

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop