Is the Large-Scale Development of Wind-PV with Hydro-pumped Storage Economically Feasible in Greece?

: The achievement of the long-term national energy targets in Greece for large-scale integration of wind and solar energy may be facilitated by the development of hydro-pumped storage projects. In light of the above, technical aspects related with the operation of the Greek power system and its ability to absorb renewable energy are analyzed in connection with the role of hydro-pumped storage and relative economic aspects. The aim of this work is to assess the potential contribution of hydro-pumped storage projects and estimate the capacity magnitude order to support large-scale wind and photovoltaic (PV) integration in Greece. For this purpose, scenarios for the Greek power system with focus on Wind and PV development, in conjunction with hydro-pumped storage capacity, are developed, and results for current situation and reference years 2030 and 2050 are presented. For the simulation, among others, high resolution mesoscale wind data for a typical year in the whole Greek territory are used for the steady state simulation of the Greek power system, in order to better estimate the power that could be generated from installed wind turbines, taking into consideration technical characteristics of a typical commercial wind turbine. Results indicate the need of gradual development of hydro-pumped storage in parallel with the large-scale integration of wind and PV capacity into the Greek power system. In addition, the feasibility of the examined scenarios is supported from the low cost of wind and PV generation. In the case of Greece, thanks to the complex morphology and hydraulic conditions of the country, hydro-pumped storage composes an efﬁcient and low-cost storage solution.


Background
The consequences of climate change due to energy production from traditional energy sources, the growing energy demand and the EU decarbonization targets lead to a need for replacement of conventional resources with renewable ones.
However, in the power sector, large-scale integration of intermittent output renewable sources is a great challenge, especially in non-interconnected or saturated grids. For this reason, the growth of large-scale wind and solar integration is a prerequisite for the simultaneous development of energy storage infrastructure.
Storage infrastructure may be a decisive factor for the size of wind and solar energy integration in the power system, as energy storage can act as a balancing component of the system. The energy storage technology that has met the biggest development, is applicable at large scale and has a considerable rate of efficiency is hydro-pumped storage.

Recent Situation and Trends
In 2019, 60 GW of wind power capacity was added globally, with the cumulative capacity reaching 650 GW. Costs have fallen rapidly for both onshore and offshore wind power, increasing dramatically the amount of bid prices in auctions around the world, especially in the past year. This development can be explained considering the constant Appl. Sci. 2021, 11, 2368 3 of 21 technology innovation, the reduced financing costs and the ever-growing competition in the industry, with electric utilities and large oil and gas companies continuing to invest further in that direction. Additionally, incentives such as robust regulatory reforms, wind energy's cost-competitiveness and its potential environmental benefits show that there is no expected deceleration in the wind power sector progress. Finally, the world's first commercial floating wind park has been commissioned in Scotland, while the sizes of turbines continued to increase, with some manufacturers producing turbines of up to 10 MW.
Meanwhile, development in solar PV sector is even faster, with almost 100 GW of new capacity added and increased the global capacity to 580 GW at the end of 2019. This growth can be explained mainly due to the raised awareness of the technology's potential to alleviate pollution, reduce carbon dioxide emissions and provide energy access to developing countries. Subsequently, intense competition, increasing efficiencies and reduction in energy costs have led to record-low auction prices.
Hydro-pumped storage advancements over the past few years continue to upgrade an already proven and reliable technology that represented more than 95% of all energy storage solutions globally, in terms of cumulative capacity. These include improved efficiencies with modern reversible pump-turbines, adjustable-speed pumped turbines, advanced equipment controls such as static frequency converters and generator insulation systems, as well as innovative underground construction methods and design capabilities, leading to faster response times of load follow for intermittent renewables more efficiently and cost effectively. Globally, there are approximately 270 pumped storage plants either operating or under construction, representing a combined generating capacity of over 127,000 MW. Of these total installations, 36 units consist of adjustable speed machines, 17 of which are currently in operation, totaling 3569 MW and 19 of which are under construction, totaling 4558 MW [23].

Current Power System Situation in Greece
The annual electricity demand in the Interconnected System of Greece, with reference year 2019, was 52.17 TWh. The total nominal capacity of conventional plants was 8806 MW (December 2019), comprising 3904 MW of lignite power plants and 4902 MW of natural gas [24]. Conventional power plants possess a significant share of electricity production which reaches approximately 51% on the annual electricity production (2019). Wind installations represent a capacity of 3283 MW, while the capacity of photovoltaics including roof-top PVs reaches 2639 MW [25]. Both technologies contributed in 2019 with 14% of generated power [24]. Moreover, the nominal capacity of large and small hydroelectric power plants was 3411 MW, biomass units accounted for 87 MW, while Combined Heat and Power (CHP) 105.47 MW [25]. During 2019, the interconnections balance, i.e., net imports, contributed to the 19% of electricity consumption in Greece. Overall, lignite units are gradually being phased out, while the share of electricity produced by natural gas and lignite power plants combined has been reduced compared to past years. Renewable sources have a smaller share in the energy mix, increased compared to the past, and will further increase in the future [25].

Prospects for the Future
The National Energy and Climate Planning (NECP) was issued in 2019 with the view to present the long-term energy targets set for the country's energy sector, in accordance with the recent European Directives. Among the targets set is the reduction of greenhouse gases up to 56% (with reference levels of year 2005), the gradual decommissioning of lignite power plants by 2028 and the increase of the share of Renewable Energy Sources (RES) in gross electricity consumption up to 60% [26].
The installed capacity of hydropower plants is not expected to increase significantly, since hydropotential has been already exploited to a great extent and environmental Appl. Sci. 2021, 11, 2368 4 of 21 constraints prevent its further development in Greece. However, there is a potential for development of reverse hydro-pumped storage, since reverse pumped storage units could be constructed within the facilities of existing hydropower plants without raising major environmental issues. The installation of storage facilities is essential, in view of the excess energy due to the higher penetration of intermittent renewable energy sources and may provide more flexibility to the system. The interconnection of mainland's power system with the small autonomous ones on the islands is essential, in order to construct a unified system, exploit further the abundant Renewable Energy Sources (RES) potential, withdraw local oil stations and secure power supply on the islands. In parallel, it is expected that the power supply in Greece will be based more on renewable energy sources, while the role of lignite will be reduced and eventually eliminated by 2028. Overall, the main target set towards 2050 is the gradual decarbonization of the power generation sector and its transition to a more sustainable future through the higher integration of renewable sources [26].

Wind and PV Development
Renewable energy sources shall substantially contribute towards the achievement of decarbonization targets set for the energy sector. One of the main objectives of energy planning is the increase of renewable energy sources' share in the energy mix, in order to substitute part of the conventional units' production. In the renewable sources capacity development scenarios, reference is mainly made to large penetration of wind and photovoltaic systems. In the corresponding studies, wind installations' capacity is projected to reach, according to Ministry of Environment and Energy in 2030, approximately 7 GW, while by 2050 11-18 GW (off-shore wind installation included) [27]. According to a study conducted by European Commission, by 2030, wind installations' capacity may reach 6 GW and by 2050 7.8 GW [28]. In addition, World Wide Fund for Nature (WWF) estimates that the wind installations will account for 5 GW, in 2030, and 6.7 GW, in 2050 [29]. Moreover, the increase of PV capacity is projected, according to the Energy Roadmap by the Ministry of Energy (the National Energy and Climate Plan and the Long-term strategy towards 2050), for 2030 at 7.7 GW, while for 2050 it is estimated at 8-12 GW, depending on the scenario considered [26,27]. According to a study conducted by the European Commission for Greece, in 2030, it is estimated that 5.6 GW of PV will have been constructed, while by 2050 9 GW [28]. Furthermore, according to WWF's projections, PV capacity will reach 4.8 GW by 2030 and 7.1 GW by 2050 [29]. Table 1  There are significant differences between the figures predicted for wind and PV capacity by the competent bodies. In order to achieve a comparative assessment of the various development scenarios, the ratio of the forecasted wind and photovoltaic capacity by the annual mean load of the country is presented in Figure 1. In 2030 the wind and PV capacity lies between 150-205% of the annual mean load demand, while in 2050 between 200-255%. The ratio between Wind and PV in terms of installed capacity is almost 1 by 1 in most of the forecasts. More PV than Wind capacity is expected by 2050 according to the Energy Roadmap of the Ministry of Energy. 200-255%. The ratio between Wind and PV in terms of installed capacity is almost 1 by 1 in most of the forecasts. More PV than Wind capacity is expected by 2050 according to the Energy Roadmap of the Ministry of Energy.

A simplified simulation of the Greek Power System
In this connection, a simplified simulation of the Greek power system is proposed to calculate the RES electricity that can be absorbed directly, the energy surplus and the ability of Hydro-pumped storage to exploit it.
Methods for solving units commitment and optimization algorithms for economic dispatch in power systems with high wind penetration have been extensively analysed and proposed in previous works [30][31][32][33]. In most of these cases, wind generation is obtained by the use of forecasting tools and unit scheduling is estimated in order to meet demand with a high probability over the scheduling horizon. It is out of the aims of the current approach to optimize the load dispatch and units commitment. The proposed methodology is based on the analysis of the steady-state operation of the Greek power system and takes into account the specific characteristics of demand, the technical features of conventional and hydropower plants and the technical constraints for the smooth and safe operation of the system. Hourly time series are used instead of a probabilistic approach [34], due to the fact that the sequence of events is essential for the simulation of hydro-pumped storage [35].
The large-scale Wind and PV integration is challenging for the power system operator due to their variable power output, the difficulty to predict wind power output accurately and the limited capacity of Greek power system with neighboring countries. However, flexible conventional and hydro units are able to cover any sudden or expected deficit, which may occur to the system. Additionally, the landscape of Greece constitutes an ideal topography for the construction of hydro-pumped storage units which could contribute significantly to the stability of the system and increase renewable peak demand supply.
For the simulation of the Greek power system a steady state analysis is performed based on hourly timeseries data. Alternative scenarios for reference years 2019, 2030 and 2050 are comparably presented in terms of energy flows and Levelized Cost of Energy (LCOE).
For the simplified simulation all the corresponding data of the power units (load factors, technical constraints, maintenance periods, stochasticity) are taken into consideration. The main target of this simulation is the accomplishment of energy balance, giving priority to renewable energy sources. All the relative assumptions are conservative in order to ensure the safe operation of the power system, with respect to all the technical constraints. The contribution of interconnections is not taken into consideration due to their

A Simplified Simulation of the Greek Power System
In this connection, a simplified simulation of the Greek power system is proposed to calculate the RES electricity that can be absorbed directly, the energy surplus and the ability of Hydro-pumped storage to exploit it.
Methods for solving units commitment and optimization algorithms for economic dispatch in power systems with high wind penetration have been extensively analysed and proposed in previous works [30][31][32][33]. In most of these cases, wind generation is obtained by the use of forecasting tools and unit scheduling is estimated in order to meet demand with a high probability over the scheduling horizon. It is out of the aims of the current approach to optimize the load dispatch and units commitment. The proposed methodology is based on the analysis of the steady-state operation of the Greek power system and takes into account the specific characteristics of demand, the technical features of conventional and hydropower plants and the technical constraints for the smooth and safe operation of the system. Hourly time series are used instead of a probabilistic approach [34], due to the fact that the sequence of events is essential for the simulation of hydro-pumped storage [35].
The large-scale Wind and PV integration is challenging for the power system operator due to their variable power output, the difficulty to predict wind power output accurately and the limited capacity of Greek power system with neighboring countries. However, flexible conventional and hydro units are able to cover any sudden or expected deficit, which may occur to the system. Additionally, the landscape of Greece constitutes an ideal topography for the construction of hydro-pumped storage units which could contribute significantly to the stability of the system and increase renewable peak demand supply.
For the simulation of the Greek power system a steady state analysis is performed based on hourly timeseries data. Alternative scenarios for reference years 2019, 2030 and 2050 are comparably presented in terms of energy flows and Levelized Cost of Energy (LCOE).
For the simplified simulation all the corresponding data of the power units (load factors, technical constraints, maintenance periods, stochasticity) are taken into consideration. The main target of this simulation is the accomplishment of energy balance, giving priority to renewable energy sources. All the relative assumptions are conservative in order to ensure the safe operation of the power system, with respect to all the technical constraints. The contribution of interconnections is not taken into consideration due to their limited capacity, and the Greek power system is considered as a remote power system. This assumption gives results on the safe side, while in reality, the management of the system could be easier with interconnections.
The innovation of the analysis is based on the use of simultaneous mesoscale typical wind year data for the representation of the wind resource in the whole Greek territory. Ninety points are used to provide high geographical resolution and cover current and future possible sites for wind farms. At each point, hourly wind power output is calculated for one year on the basis of hourly wind speed and wind power installed capacity. Then the aggregated hourly wind power in the whole Greek territory is used as an input for the simulation of the Greek power system. Typical wind turbines power curves are used to calculate wind power output, given the wind installed capacity and the time series of wind speed at each point under consideration. Then the aggregated wind power output P W prod for every hour is resulted.
The main principles of the methodology, units' commitment and load dispatch are presented in this section in four steps. In the first step PV power is absorbed in priority; in the second one, conventional units commitment is defined; in the third one, wind power absorption is calculated, and in the last one, Hydro and Hydro-pumped storage output is resulted.

Step 1. PV Power Absorbed in Priority
In the step 1 PV power is absorbed in priority. Wind and PV will gradually represent the mainstream of electricity in the Greek power system; then the management rules and operational principles of the power system should be reformed in comparison with the recent practices [35]. One of the prime assumptions of the current approach is that PV energy is absorbed in priority. PV production is considered as a predictable source of energy, and it is allowed to be absorbed in priority without any constraints, unless it surpasses the load demand. In a large area under consideration, like Greece, geographical distribution of PV has a positive effect on the smoothing of PV power output fluctuations and on the predictability of PV aggregated power output. For every hour i of the year (i = 1 to 8760), given the initial actual demand P D and subtracting the PV power output P PV , the residual load P D-PV is calculated (Equation (1)) to be used as the basis for the commitment of conventional units:

Step 2. Conventional Units Commitment
Lignite power plants are used today as a base load units. They are not flexible to switch on and switch off, and their operation should be scheduled in medium term. The number of lignite power plants to be committed is estimated with respect to the expected demand of next days, taking into consideration their non-flexible features and their high technical minimums of the order of 50%. The criterion is that the technical minimums of lignite units to be committed should not exceed the minimum expected residual load of the following seven days. Initially, there are set to operate at their technical minimums or more in case of power deficit. Given this simulation approach, the operation of lignite power units in the Greek power system will be gradually reduced as far as the PV capacity increases. It should be noted that lignite power plants are considered only for reference year 2019, since by 2028 according to national targets the lignite units shall have been decommissioned.
Natural gas units are more flexible than the lignite units, and as a result, they are dispatched in order to support the balance between supply and demand in the short term. The number of natural gas units operating is defined in order to be able to cover the mean residual load of the next twenty-four hours. Initially, they are set to operate at their technical minimums (30% of the nominal output) or more in case of power deficit. The flexibility of natural gas units allows the efficient and rapid response to the fluctuations of power supply.
The combined heat and power (CHP) plants and the biomass units are dispatched hereafter. The total nominal power of these units is relatively small, and their contribution is limited. The produced energy of both is defined using an average capacity factor of 70%.

Step 3. Wind Power Absorbed and Curtailed
After the definition of conventional units commitment, the absorbed P Wabs and curtailed P Wcurt wind energy are calculated. The wind power allowed to be absorbed by the system is limited by technical restrictions, in order to ensure the stability of the system. A percentage of instantaneous penetration δ of wind power to the grid is applied in the calculations [35,36]. This is a dynamic limit which ensures the safe operation of the system in case that all the wind power is lost, and it is related with the ability of other units (i.e., conventional, hydro and storage) to supply load demand [36]. On the other hand, wind power absorption should respect the technical minimums of the committed conventional power plants [35]. Technical minimums of lignite, natural gas, CHP and biomass are taken into consideration. These two technical restrictions are related with power system and define the ability of the power system to absorb wind power. Finally, the actual wind power absorbed P W abs cannot exceed the available wind power P W prod , as it has been calculated using the mesoscale data. Then, the actual absorbed wind energy is resulted (Equation (2)): Meanwhile, the wind energy curtailment is calculated (Equation (3)):

Step 4. Hydropower Output and Hydro-Pumped Storage
Hydropower plants supply the peaks of the load curve. The power produced annually from hydropower plants is highly dependent on the hydrological year. In a good year in terms of hydrological conditions, the annual production could reach 5 TWh. A special methodological approach is used for hydropower [35]. Before the integration of PV in the Greek power system, hydropower was used for peak demand supply. At that time peak demand occurred especially in summer period. Today, hydropower has a more complex role. The integration of PV has shifted the need for peak supply from summer to winter. Hydropower will continue to provide peak power supply whenever peak load demand occurs. Additionally, hydropower will balance the variability of rest RES power generation. Renewable power surplus and wind curtailment may occur in low demand periods or in windy ones. During low demand hours hydropower plants are switched off. During peak demand periods, if there is wind power surplus, hydropower plants may reduce their operation saving water for peak demand periods of low wind. So, wind power plants could save water in the hydro plants' reservoirs and hydro generation will not constrain wind power absorption. Good and bad hydrologic years occur in Greece, with an average annual hydro energy output in Greece of 5 TWh. This energy amount is split into the hours of peak demand using a simplified iterative approach, which finally defines the threshold of power demand for hydropower generation. The operation of hydroelectric power plants is defined by the "peak-shaving" method; it is defined by a minimum power demand threshold (P T_H ), when this load demand limit is surpassed hydroelectric power is produced. This limit is determined after an iterative procedure (goal seek function in MS Excel is used) for the whole year, so that the integral of hydro annual production will be equal with 5 TWh (typical hydraulic year). In this step, the remaining load demand after subtracting PV, lignite, natural gas, CHP, biomass and wind power absorbed is used (P D-PV-CONV-W ). Then, hydropower output is defined for every hour of the year (Equations (4) and (5)) with respect to the nominal power of hydro plants in Greece. This approach is very close to the reality, since existing hydro-plants are operated today in the way of seasonal storage and peak supply. Additionally, the PV's summer peak supply match well with the dry period of Greece from May to September.
By this approach, the power surplus (wind curtailment and PV surplus) P surplus is minimized in most of the conservative scenarios. Obviously, for large-scale wind and PV integration, due to the technical constraints, power surplus occurs and could be transformed through hydro-pumped storage to useful peak demand supply. The aforementioned excess energy is stored in hydro-pumped storage units, by pumping water to the upper reservoir when there is a surplus of energy, while it is recovered through the hydro-turbines operation when other renewable sources are not available.
Due to the distribution of the curtailed power, it is not economically feasible to exploit 100% of this energy. Such a scenario would require enormous installed capacity of pumps and volume of the upper reservoirs which would be used only few hours per year. A sufficient degree of annual exploitation of excess system's power can be considered to be 70% [6]. On this basis, the required installed capacity of pumps P P_nom is defined in order to achieve the target of 70% exploitation of the annual energy surplus. The hydro turbines capacity is considered to be equal with the pumps' nominal output. Therefore, the nominal power of the pumps is calculated (Equations (6) and (7)): where P pump is the power used for pumping, P P_nom is the nominal power of pumps and P surplus is the available power surplus (wind curtailment and PV surplus). The volume of the upper reservoir is defined to ensure a 24-h operation of the turbines at their nominal capacity.

Wind Data-Mesoscale Modeling
Simultaneous information on wind statistics over every potential area for wind farm development is required for this analysis. Even if a large number of wind measurements are available, it is practically difficult to represent simultaneous data series and cover every potential area of interest. Installation of a mast network for this purpose could lead to rather prohibitive technical and economic restrictions. Additionally, existing wind monitoring networks are relatively large and can provide large spatial coverage but not necessarily high resolution [37]. On the other hand, use of wind potential maps is not a solution since they only provide an estimation of the spatial distribution of the mean wind speed without any information on its temporal variation. Application of a Numerical Weather Prediction (NWP) model can effectively provide the information required.
In this connection, high resolution analytical wind data timeseries for typical wind year are used. These data have been produced by the systematic application of a numerical weather prediction model. Analytical presentation and description of the approach was given in [38]. In Figure 2, the relative high resolution wind atlas of Greece is presented in terms of power density and parameter c of Weibull distribution, which are the most common ways to present Aeolian maps. The wind atlas of Greece was based on a typical wind year and 12 months of weather model simulations for grid boxes 2 × 2 km 2 in size. The numerical weather prediction model used is "MM5" which is run operationally at the National Observatory of Athens since 2002 [39] and has been verified [40,41] for its forecast skill over the area of interest.
production. However, as the spatial dispersion is increased and more wind farms are introduced, the time periods with cumulative zero or rated production are reduced [34]. The aggregated hourly wind power output is calculated on the basis of the mesoscale wind data, and the selected points in the whole Greek territory for the current wind energy development (2019) and the scenarios of installed capacity in each point under consideration for 2030 and 2050 (Figure 4). In Figure 2, high resolution wind atlas are presented for the typical wind year [38]. The indexes of Power density (W/m 2 ) and Parameter c of Weibull distribution (m/s) are presented. Both are widely used in Aeolian maps for representation of the wind potential [42,43]    The selection of 90 representative points in the Greek territory is based on the location of wind farms as it is expressed by investors' interest and it is depicted in the Regulatory Authority for Energy (RAE) geographical information system in the Greek territory ( Figure 3). For future scenarios, spatial dispersion of wind farms and existing plans for interconnections of Greek islands with the power system of mainland Greece are considered. Data for the selected 90 points are presented in Table 2. The duration curve of wind power output is dependent on the spatial distribution of wind farms. Historical data of single wind turbine power production typically show extended time periods with zero or rated production. However, as the spatial dispersion is increased and more wind farms are introduced, the time periods with cumulative zero or rated production are reduced [34]. The aggregated hourly wind power output is calculated on the basis of the mesoscale wind data, and the selected points in the whole Greek territory for the current wind energy development (2019) and the scenarios of installed capacity in each point under consideration for 2030 and 2050 ( Figure 4). In Figure 2, high resolution wind atlas are presented for the typical wind year [38]. The indexes of Power density (W/m 2 ) and Parameter c of Weibull distribution (m/s) are presented. Both are widely used in Aeolian maps for representation of the wind potential [42,43] Figure 3 shows the overview of wind farms applications in the Greek territory and selection of 90 representative points in the whole Greek territory.
In Table 2 the details of the 90 representative points are presented (name, location, k and c of the Weibull distribution and wind power density) [38]. Appl. Sci. 2021, 11, x FOR PEER REVIEW 10 of 21  In Table 2 the details of the 90 representative points are presented (name, location, k and c of the Weibull distribution and wind power density) [38].   In Table 2 the details of the 90 representative points are presented (name, location, k and c of the Weibull distribution and wind power density) [38].

Load Demand Data
Actual time series of load demand data for the interconnected power system have been used (Power Public Corporation S.A. data, https://www.dei.gr/en, accessed on 15 January 2021). Corresponding adjustments to the demand time series were realized in order to formulate the corresponding timeseries for the years 2030 and 2050. The base year for the load demand time-series is 2006. This is the last year before the start-up of PV development in Greece. The forecasts for annual electricity demand and peak power demand are based on relevant studies that have been conducted for the power system of Greece. The comparative study of the researches carried out for Greece's future demand concludes that, in 2030, a modest estimation is considered to be 57.2 TWh with a peak demand of 10.5 GW, while for 2050, 74 TWh with a peak demand of 13 GW. The studies that were taken into consideration were conducted by the Ministry of Environment and Energy [27], the European Commission [28] and WWF [29]. A comparative representation of the electricity demand forecasts is presented in Figure 5.

Load Demand Data
Actual time series of load demand data for the interconnected power system have been used (Power Public Corporation S.A. data, https://www.dei.gr/en). Corresponding adjustments to the demand time series were realized in order to formulate the corresponding timeseries for the years 2030 and 2050. The base year for the load demand time-series is 2006. This is the last year before the start-up of PV development in Greece. The forecasts for annual electricity demand and peak power demand are based on relevant studies that have been conducted for the power system of Greece. The comparative study of the researches carried out for Greece's future demand concludes that, in 2030, a modest estimation is considered to be 57.2 TWh with a peak demand of 10.5 GW, while for 2050, 74 TWh with a peak demand of 13 GW. The studies that were taken into consideration were conducted by the Ministry of Environment and Energy [27], the European Commission [28] and WWF [29]. A comparative representation of the electricity demand forecasts is presented in Figure 5.

PV Data
In terms of PV production, actual time series provided by CRES (Centre of Renewable Energy Sources) were used with the appropriate adjustments. The adjustment is based on PVGIS estimated annual output [44]. Table 3 present the information used from PVGIS and representative duration curves of PV power output for four cases. Table 3. Annual PV production [44].

Scenarios
By 2030, the electricity demand in the country is expected to reach 57.15 TWh, with a peak of 10 GW and by 2050, 74 TWh with a peak of 11.3 GW [27].
Lignite power plants are expected to be decommissioned by 2028, in view of the achievement of decarbonization targets. The nominal output of natural gas units for reference years 2030 and 2050 is expected to reach 6.97 GW [27] and 7.1 GW [29].

PV Data
In terms of PV production, actual time series provided by CRES (Centre of Renewable Energy Sources) were used with the appropriate adjustments. The adjustment is based on PVGIS estimated annual output [44]. Table 3 present the information used from PVGIS and representative duration curves of PV power output for four cases. Table 3. Annual PV production [44].

Scenarios
By 2030, the electricity demand in the country is expected to reach 57.15 TWh, with a peak of 10 GW and by 2050, 74 TWh with a peak of 11.3 GW [27].
Lignite power plants are expected to be decommissioned by 2028, in view of the achievement of decarbonization targets. The nominal output of natural gas units for reference years 2030 and 2050 is expected to reach 6.97 GW [27] and 7.1 GW [29].
In this connection, the reference scenario by 2030 for wind and PV capacity refers to a cumulative capacity of up to 16 GW (8 GW wind, 8 GW PV). By 2050, it is assumed that the renewable energy sources' capacity, i.e., wind and PV, will account for 18 GW aggregated, with wind installations of 9 GW and photovoltaic of 9 GW. The normalized wind and PV capacities (by the average annual demand) are 1.83 by 2030 and 2.48 by 2050, close to the relative normalized capacities in the corresponding studies discussed and presented in Figure 1.
The percentage of instantaneous wind penetration (δ) is considered to be 50% for 2030 and 60% for 2050. The increase of the instantaneous penetration (δ) in 2050 is based on the fact that the management of the grid renewable energy will have matured, and weather load forecast models will be widely used operationally. However, both figures are considered as very conservative approaches, which will keep the results on the safe side.
According to IPTO adequacy study for 2020-2030, as of December 2019, 700 MW of hydroelectric plants have been licensed, including 590 MW of hydro-pumped storage facilities [45].
The nominal output of power plants using as feedstock biomass is considered to reach 300 MW by 2030 and 600 MW by 2050 [27].

Energy Mix by 2030 and 2050
By 2030 and 2050, considering the developments and assumptions mentioned above, the final energy mix could be formed as depicted in Figure 6.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 14 of 21 In this connection, the reference scenario by 2030 for wind and PV capacity refers to a cumulative capacity of up to 16 GW (8 GW wind, 8 GW PV). By 2050, it is assumed that the renewable energy sources' capacity, i.e., wind and PV, will account for 18 GW aggregated, with wind installations of 9 GW and photovoltaic of 9 GW. The normalized wind and PV capacities (by the average annual demand) are 1.83 by 2030 and 2.48 by 2050, close to the relative normalized capacities in the corresponding studies discussed and presented in Figure 1.
The percentage of instantaneous wind penetration (δ) is considered to be 50% for 2030 and 60% for 2050. The increase of the instantaneous penetration (δ) in 2050 is based on the fact that the management of the grid renewable energy will have matured, and weather load forecast models will be widely used operationally. However, both figures are considered as very conservative approaches, which will keep the results on the safe side.
According to IPTO adequacy study for 2020-2030, as of December 2019, 700 MW of hydroelectric plants have been licensed, including 590 MW of hydro-pumped storage facilities [45].
The nominal output of power plants using as feedstock biomass is considered to reach 300 MW by 2030 and 600 MW by 2050 [27].

Energy Mix by 2030 and 2050
By 2030 and 2050, considering the developments and assumptions mentioned above, the final energy mix could be formed as depicted in Figure 6. The estimated nominal output of hydro-pumped storage units reaches 1500 MW for the reference year 2030. The capacity factor of pumped-storage units, expressed by the rated power of the turbines, reaches 19%. This value of capacity factor is below the lower acceptable limit which renders the investment economically viable (a benchmark of the capacity factor could be 25% [6]). In this study, this value is considered acceptable, since, according to current data for Greece, large hydro units' capacity factor varies to similar levels (15-17%); therefore, this value could also be considered for reverse hydro-pumped storage projects. The capacity factor of hydropower units may be lower than benchmark, because they are used as a safety net, for cases of emergency, so they are over-dimensioned in order to be able to support the system when a deficit occurs.
The energy contribution of wind and PV installations reaches 49%, and the energy contribution of conventional units is lowered to 31%. The electricity sector is characterized by a higher renewable share, and electricity production is partially decarbonized, compared to the current situation. Hydro-pumped storage units contribute with 6% to the annual electricity demand, by exploiting 70% of curtailed energy.
In 2050, the installed capacity of wind and photovoltaic installations is considered to be 24 GW in total, which results in a reduction of the conventional power plants energy contribution to 25%. Renewable energy sources possess the highest share of production, accounting for approximately 60%. The energy surplus increases due to the higher integration of renewables, and the hydro-pumped storage capacity is estimated at 2700 MW. The share of hydro-pumped storage projects reaches 7% of the total energy produced. The capacity factor of hydro-pumped storage units is estimated at 23%. This higher capacity factor could be attributed to the higher integration of RES, which results in the increase of curtailed power. Consequently, the higher integration of RES contributes to the enhancement of the economically viable operation of hydro-pumped storage units.

Economic Assessment
To answer the question whether large-scale integration of Wind-PV with the parallel development of hydro-pumped storage is feasible, the levelized cost of energy is used. The levelized cost of energy is calculated as follows: where, CAPEX is the capital expenditure, CRF the annuity factor, E the annual energy produced, OPEX the operation and maintenance (O&M) expenditures, DR the discount rate and N the lifetime of the investment. The levelized cost of energy is a valuable cost indicator because of its relative simplicity which allows the comparison of different technologies. At the same time, it is a useful tool for benchmarking the cost of different technology units, considering the lifetime of the units.
The levelized cost of energy has been calculated for the current power supply system and for the examined scenarios for 2030 and 2050, based on the simulation results and the data collected [26][27][28][29][46][47][48][49][50][51][52][53][54][55][56][57][58][59][60][61][62][63][64][65]. Studies conducted and demonstrating projections for costs data, such as CAPEX and OPEX, have been considered for the calculation of the LCOE. The final values used constitute an average value of the costs identified in other publications. A discount rate DR equal to 7% is considered in all cases to secure that estimations of LCOE are modest and relative risks are considered. After a long period of economic uncertainty in Greece due to the economic crisis, in the last 3 years, Greek government bonds had constantly relatively low rates, less than 3%. Tables 4-6 show the details about LCOE components, data, assumptions and results for the years 2019, 2030 and 2050, respectively. All the prices presented in the three tables are not harmonized. In some cases, different sources are available and mean values are used. The extracted value of LCOE evaluates not only the relative CAPEX and OPEX for each type of units but also their utilization, as it is introduced by the energy output in the denominator. Then, the results of LCOE are useful for the specific cases studies in Greece, under the assumptions of the current methodological approach and conclusions cannot be universalized. Additionally, the various parameters of cost in case of hydro and hydro-pumped storage are site dependent, and then, results are not useful for other case studies.  In Figure 7, the size of the bubbles corresponds to the amount of energy generated by each type of technology, whereas on the horizontal axis the LCOE of each unit type is presented, and on the vertical axis, the share of the cost of the different types of technologies in the total cost of the system is depicted.
In the current power supply system, the participation of lignite and natural gas units is crucial in order to meet electricity demand. The LCOE of the system is estimated to 107.59 €/MWh, and the largest share of the cost is due to the conventional units. Although CHP and biomass units are characterized by high levelized costs due to their limited nominal capacity, they do not participate significantly to the final composition of the cost of the system.
The weighted average LCOE of the power system in Greece by 2030, according to the calculations, is decreased to 79 €/MWh. This reduction could be attributed to the decommissioning of lignite units, which present a higher LCOE compared to RES (wind and PV installations) in 2019. Natural gas units gradually replace lignite ones by 2030. Carbon dioxide allowances are expected to rise to a higher level by 2030; therefore, the operational expenditures of conventional units are increased, resulting in a higher LCOE for natural gas units. The expected decrease of the CAPEX required for RES installations contribute to the significant decrease of the weighted average LCOE of the system. By 2050, a further reduction of the system's cost is expected; the levelized cost is approximately 75 €/MWh, reduced by 30% compared to current levels. Increasing the participation of renewable energy sources in the energy production with the simultaneous significant reduction of their CAPEX leads to a decrease in the levelized cost of the system. Hydro-pumped storage is among the technologies presenting the lowest LCOE, following wind and PV. The CO 2 prices are significantly higher compared to 2030 resulting in the increase of the levelized cost of energy of natural gas power plants. In the current power supply system, the participation of lignite and natural gas units is crucial in order to meet electricity demand. The LCOE of the system is estimated to 107.59 €/MWh, and the largest share of the cost is due to the conventional units. Although CHP and biomass units are characterized by high levelized costs due to their limited nominal capacity, they do not participate significantly to the final composition of the cost of the system.
The weighted average LCOE of the power system in Greece by 2030, according to the calculations, is decreased to 79 €/MWh. This reduction could be attributed to the decommissioning of lignite units, which present a higher LCOE compared to RES (wind and PV installations) in 2019. Natural gas units gradually replace lignite ones by 2030. Carbon dioxide allowances are expected to rise to a higher level by 2030; therefore, the operational expenditures of conventional units are increased, resulting in a higher LCOE for natural gas units. The expected decrease of the CAPEX required for RES installations contribute to the significant decrease of the weighted average LCOE of the system.
By 2050, a further reduction of the system's cost is expected; the levelized cost is approximately 75 €/MWh, reduced by 30% compared to current levels. Increasing the participation of renewable energy sources in the energy production with the simultaneous significant reduction of their CAPEX leads to a decrease in the levelized cost of the system. Hydro-pumped storage is among the technologies presenting the lowest LCOE, following

Discussion and Conclusions
High renewable penetration could lead to a significant reduction of the system's cost, while hydro-pumped storage systems may contribute to the peak demand supply.
The size of the reverse hydro units should be proportional to the development of Wind and PV installed capacity. The over-dimensioning of reverse hydroelectric projects may lead to installed hydroelectric pumps and turbines which will be used only for a few hours per year without ensuring their economic feasibility, while under-dimensioning will lead to insufficient exploitation of the potential. Furthermore, the higher integration of renewable sources (wind and PV) across the country (spatial dispersion) results in significant energy surplus and curtailment during the year, therefore energy that could be stored in hydro-pumped storage systems is increased. In this case, hydro-pumped storage units can contribute significantly to the energy balance, as analyzed in the scenarios.
The economic assessment of the future power system of Greece, based on the scenarios examined in this connection, demonstrates that the increase of RES contribution may lead to lower costs for the system. Greater penetration of photovoltaic and wind installations results in a reduction of the levelized cost of energy of the system, since capital expenditures required for wind and PV systems are expected to decrease significantly over the years. The cost of thermal units in reference years 2030 and 2050 is expected to be higher due to the increase in the cost of carbon dioxide emission allowances. The reduction of the energy production share of the latter in combination with the higher integration of RES prevents the increase of the weighted average LCOE of the power system. Author Contributions: Conceptualization, G.C. and A.Z.; methodology, G.C. and A.D.; software and formal analysis, A.D.; investigation, A.D.; data curation, V.K. and K.L.; writing-original draft preparation, A.D., T.C.; writing-review and editing, A.D. and G.C.; supervision, G.C.; All authors have read and agreed to the published version of the manuscript.