2.1. Simulation Approach
Three promising locations for grid-supporting electrolysis in the distribution grid of ‘Energienetze Steiermark’ (located in Styria, Austria) are considered.
Figure 2 illustrates the network area of ‘Energienetze Steiermark’, in line with the plans of the H2 Roadmap [
22] and the European Hydrogen Backbone (EHB) [
23]. Within this network area, three locations were selected. For confidentiality reasons from side of the DSO, the exact location of the three P2G systems cannot not be disclosed, but for electrolysis to make sense as a substitute to conventional grid-enforcement measures, potential locations need to fulfill a number of requirements. First and most obvious, grid-bottlenecks caused by future infeed and limitation of current grid capacity are expected and need to be solved. Second, and this is probably related to the first condition, these locations will be close to areas with high expected deployment of intermittent renewable generation. Third, proximity to the future hydrogen network is considered to be essential as a sink for the produced hydrogen. Put differently, electrolysis will be positioned at the intersection of i) high-voltage electricity grid, ii) envisaged RES build and iii) the envisaged hydrogen network. In these three selected locations, the surrounding grid was modelled and the P2G systems were positioned.
The overall simulation is carried out for a timeframe of 30 years (2025-2055) to show the increasing need of grid reinforcement due to increasing implementation of RES generation (mostly PV, but also wind), while also considering increasing electricity demand (private EV and private heat pump rollout scenarios where considered). To represent a discrete expansion of P2G capacity to counteract grid congestion by increasing feed-in, a five-year expansion cycle was implemented. In the beginning of each five-year simulation cycle, the maximum grid congestion was taken as reference for expansion capacity of the upcoming five years.
The modelling was undertaken against historical data of 2023 with respect to RES generation and load profiles, electricity prices and grid charges. Strong expansion of RES was assumed in accordance with [
1,
24], as further outlined in
Section 2.2. Despite some efforts to predict and estimate future prices or costs of hydrogen (see for example [
25,
26,
27]), we consider the hydrogen price as the most uncertain variable in this modelling exercise. Therefore, the hydrogen price is exogenously given and varied in a range between 1 and 10 EUR/kg.
1 In all locations, it is assumed that the P2G is grid-connected (front of meter) and is not working co-located being directly connected to the plant (behind the meter). This implies that grid charges are applicable to their full extent.
2.2. Grid Simulation
As a first step, a grid simulation was undertaken to determine the grid-supportive operation. The grid-supportive operation was defined as the minimal required operation of the electrolyzer, to compensate peaks of renewables production which would have led to violation of grid constraints and as a consequence of infeed curtailment. It is a rule-based strategy and is determined on the basis of grid simulation and the respective technical limits of network components.
This is said to be the minimal required operation of the P2G to consume renewable generation and thus maintain valid operational state of the grid. Hence, the P2G is designed in order to compensate the maximum grid-violation that can be observed within the modelled timeframe. On the other hand, the P2G load is constrained from the upper side the grid is overloaded from the demand side. This typically happens when renewable infeed is low or zero and consumption is high. Considering both limitations from the upper and lower side, the P2G is left with an operational window, reflecting the minimum load required to absorb production peaks and constraints from the upper side reflecting demand peaks.
Figure 3 depicts such an operational window in a simplified way by meeting primary substation capacity limit as an example of a bottleneck of a distribution grid. In this example, grid components would be overloaded from excess feed-in between 9:00 and 17:00. The maximum overload is 10 MW at 13:00. Therefore, the P2G is designed at a capacity of 10 MW. However, when running at full capacity, it would lead to a congestion between hour 2 and 7 as well as 19 and 23, hence, the operational window is reduced during these timeframes by the required amount. Typical bottlenecks of distribution grids are line loading limits (predominantly in network level 3), primary substation capacity (predominantly for network level 4 – see
Figure 3), or grid voltage limits (predominantly network level 5).
The grid simulation was built in Digsilent PowerFactory
2 and the whole calculation was automatized by the PowerFactory Python interface, so future scenarios and controls where setup using Python, and PowerFactory was fully controlled via Python. For yearly simulation, yearly consumption and generation profiles were configured in the grid, and the PowerFactory Quasi-Dynamic-Simulation function was used to calculate one year within 15 min timesteps with high performance. Yearly consumption and generation profiles were provided by the DSO, containing real measurement data from primary substation transformers and feeders in the historical timeframe from 2023-01-01 to 2023-12-31 in 15-minutes granularity.
Figure 4 shows the grid topology of all three locations. All three locations are characterized by different grid layouts and therefore, the P2G is connected at different network levels spanning from network level 3 to network level 5.
Figure 5 illustrates the assumptions on RES expansion and
Figure 6 depicts the assumptions for growth of demand due to electrification of heat and mobility in all three locations. Location A and Location C are characterised by low-voltage (LV) and medium-voltage (MV) photovoltaic rollout, while Location B experiences also some increase in wind capacities. The assumptions are based on multiple sources, e.g., [
28,
29] and were fine-tuned based on several discussions with the DSO, being aware of the local conditions. In all locations, significant expansion of RES as well as additional demand from heatpumps and electric vehicles (EVs) is expected.
Location A
The P2G is supposed to relieve the substation on network level 4. The current demand/injection (historical data 2023) spans between +9 MW (demand) and -17 MW (injection). The transformer station capacity (considering n-1 security) lies between +50 MW and -100 MW. Strong growth of solar PV on low voltage (LV) and medium voltage (LV) is assumed.
Location B
The P2G is supposed to relieve a grid branch on 110 kV level (network level 3). The current demand/injection spans between +50 MW and -15 MW. The line capacities are between +100 MVA and -200 MVA for “lower” line and +200 and -400 MVA for “upper” line. Strong growth of wind and solar PV is anticipated.
Location C
One feeder of a 20 kV medium voltage (MV) grid is modelled that is characterized by two main branches. Current peak load in this feeder is around 2.3 MW and maximal RES generation at minimal load causes a reverse power flow -0.63 MW. To increase hosting capacity to enable the expected amount of PV integration, one P2G was located in each of the branches on network level 5. While one of the branches is around 19 km long and the other around 22 km, the P2G were placed on the most optimal location in the feeder to compensate future PV infeed most effectively. According to current common practice in network planning, network level 5 consumers are not allowed to consume more than 8 to max. 10 MW, which constitutes a non-technical limit for the size of the P2G.
2.3. Market Based Operation
The optimal operation based on market signals is determined by solving an optimization problem, implemented in the IESopt framework [
30]. IESopt is an open-source optimization framework built on top of the JuMP modeling language [
31], written in Julia [
32] with a Python interface.
For the market based operation, the P2G system is modelled as a price taker against the electricity market. The optimal market based operation is a result of the marginal costs of hydrogen production and the hydrogen price. While the hydrogen price is assumed to be exogenous, the marginal costs of hydrogen production is determined by the electricity price, the grid charges and the conversion ratio of the P2G. For the electricity price, quarter-hourly historical spot market prices for Austria in the year 2023 are considered. The grid charges are taken from historical data for the network area of Styria and consist of an energy component in EUR/MWh and a peak-load component in EUR/MW for the monthly quarter-hourly consumption peak, as described in
Table 1. The optimization maximizes yearly profits by choosing electricity consumption of the P2G, being constrained by the size of the P2G and facing electricity prices and grid charges as costs. The output value of the P2G is determined by the exogenously assumed hydrogen price.
The modeling assumes perfect foresight on the electricity prices over one entire month, because peak-load grid charges are settled monthly in Austria. Hence, market based operation is determined month-ahead on a quarter-hourly resolution. Since the market based operation is implemented as linear program (LP), the conversion ratio of the P2G as produced hydrogen mass per electric input is considered at a constant value of ~18.3 kg/MWh. The optimization always chooses full load as long as the marginal costs of production fall below the hydrogen price, considering fixed costs arising from peak load grid charges. Therefore, no minimum load arising from the technical specifications of PEMEL is considered in the optimization model.
The monthly operation schedule determined by the optimization model is then used as in input in the technical simulation model. The technical simulation combines the operational schedules of the grid simulation and the market based operation. The technical simulation makes sure that the operational schedule of the P2G lies within the operational window, as illustrated in
Section 3. In addition, the technical simulation model considers the non-linearities and varying conversion factors as a function of load, which provides a preciser inputs to the cost-benefit analysis when the P2G is working in partial load. This is particularly the case when the P2G is working in grid-supportive mode only or the market based operation is limited from the upper side due to grid constraints.
2.4. Technical Simulation
The technical simulation corrects the results of the market optimisation and aims to adequately reflect the technical characteristics of the P2G. To properly absorb renewable production peaks, the selected electrolysis technology must be appropriately flexible. Therefore, due to its higher flexibility in dealing with power ramping, PEM technology was chosen for the electrolyzer [
33,
34]. To represent the technical lifetime of the PEMEL, it is assumed that after end of life, the electrolyzer must be removed or replaced. Regarding lifetime, literature diverges quite extremely by giving values of 20,000 to 90,000 hours of operation or 10 to 25 years [
14,
35,
36,
37]. Therefore, a lump-sum lifetime of 15 years was assumed, even though it might seem a little overestimated for a possible full load operation (e.g., a total lifetime of 90,000 hours leads to a capacity factor of around 0.68 for a lifetime of 15 years). Assuming a lifetime of 15 years when considering a timeframe between 2025 and 2055 implies reinvestment around 2040. Reinvestment is then evaluated at the relevant CAPEX value of the respective year.
The system simulated in this study is designed to provide a realistic representation of a P2G system. Starting at the coupling point to the electrical distribution grid, alternating current is converted into direct current using an ACDC converter. For location B, an additional transformer is needed at the substation due to technical specifications.
3 Direct current is fed into the PEMEL. The output of the electrolyzer would normally be buffered in a small-scale hydrogen storage tank (with a storage time of about two hours based on the electrolyzer output) to ensure a constant or quasi-constant back pressure to the hydrogen output stream of the electrolyzer. Due to the limited losses of the storage tank and its small size, which makes it only usable for technical purposes, it is neglected in the simulated system. The compressor required to increase the pressure of the hydrogen gas before feeding it into the pipeline is assigned to the gas grid. The compressor therefore has no influence on the P2G in technical or economic terms and is also neglected.
The technical modelling is carried out using TESCA, an internal AIT framework for deterministic time series simulations [
38]. The simulation framework provides models for various components of energy systems in the electricity, hydrogen and heat sectors. For this study, the electricity and hydrogen grids were implemented as unlimited source and unlimited sink, respectively. Therefore, the P2G system was only restricted by its own technical boundaries, like nominal power of the components. The transformer and the ACDC converter were implemented with load-dependent efficiency curves.
In contrast to the models of transformer and ACDC converter, integration of the core component of the P2G, the PEMEL, into the entire system is based on in-depth modeling of the electrochemical process from [
39]. The efficiency of electrolysis
is calculated as a combination of voltage efficiency
and faraday efficiency
. Voltage efficiency on the one hand accounts for electric losses in the PEM cell. It is computed as outlined in (
1), where thermoneutral voltage is
and total cell voltage
as defined in (
2).
The reversible voltage depends on both the operational temperature and pressure, and activation overpotentials of anode and cathode accounting for kinetics at the electrodes, and ohmic overvoltage accounting for ohmic losses in the membrane. As activation overpotentials increase with the natural logarithm of the applied current and ohmic losses linearly to the applied current, voltage losses are increasing for increasing applied electrical power. Faraday efficiency on the other hand accounts for back-diffusion of hydrogen through the membrane by the pressure difference, reducing the actual hydrogen throughput of the electrolyzer. Hydrogen back-diffusion is not depending on the theoretical throughput, but rather on the pressure difference between cathode and anode side (which are assumed constant as 30 bar at the cathode and 1 bar at the anode side in this study). With constant hydrogen back-diffusion the faraday efficiency increases to 100 % asymptotically in dependency of the load. For auxiliary technology such as water pumps, water treatment and cooling system, a constant auxiliary power of 5 % of the electrolysis nominal power is assumed, reducing the stack efficiency even further to an overall system efficiency of the electrolysis system . Due to the constant nature ot the balance of plant (BoP) power the influence on the overall system efficiency is decreasing with increasing load of the electrolyzer.
2.5. Cost-Benefit Analysis
A cost-benefit analysis (CBA) is undertaken over the entire simulation time frame (e.g., between 2025 and 2055, 30 years). To evaluate measures of the P2G against conventional grid enforcement, a baseline scenario for the entire time frame is defined. In the baseline scenario, the surplus renewables generation is curtailed so that the network can still be operated in its technical boundaries. Then, both scenarios (classic grid enforcement-scenario and P2G-scenario) can be compared to the baseline scenario.
Figure 7 shows the main characteristics of the baseline scenario: curtailed energy in GWh and lost revenues in thousand euros from sales on the electricity market. The curtailed generation is a direct consequence of the renewables expansion in combination with the congested grid components and the lost revenues on the electricity market consist of the product of the curtailed generation and the historical electricity prices 2023.
At this point, it is worth mentioning that the cost-benefit analysis takes a global perspective, without a specific view on the ownership structure of the components. For example, lost revenues on the electricity market concern the owner of the renewable plants, not the DSO or the owner of the electrolyzer. However, in a global cost-benefit analysis, revenues and expenses can be added without specific attribution to actors. The composition of costs and benefits from a global perspective are illustrated in
Figure 8. First, it can be concluded that, independently of the hydrogen price, the recovered revenues from renewables generation (that would be curtailed in the baseline scenario) contribute to the overall benefits. Assuming a hydrogen price of zero, no market based operation could be observed, but an (unprofitable) grid-supportive operation would maintain network operation within its limits and ensure that all renewables production could be injected into the grid. The computation of the lost revenues is shown in (
3), where
represents the curtailed electricity in year
y at time
t and
the corresponding electricity price in the respective time interval
t. This calculation is repeated and then summed up for each year
y between 2025 and 2055, whereby a discount factor
r (assumed at 4% in real terms) is applied. Index
t represents all quarterly hours of the year and therefore runs from 1 to 35,040.
However, it is obvious that global costs and benefits depend on the assumed hydrogen price. With an increasing hydrogen price (equivalent to willingness to pay for the output of the P2G), revenues from hydrogen sales increase. This is a result of two sources: first, hydrogen output is evaluated at a higher price and, second, the output of the P2G increases since market based operation finds more profitable hours. At some point, we can observe a ‘break-even hydrogen price’ where benefits start to outweigh the costs. Please note that the curves depicted in
Figure 8 are drawn for illustrative purposes only. In fact, both curves are non-linear and proportions are very different.
For each location, costs and revenues are accumulated and discounted over the entire time horizon from 2025 to 2050. Equation (
4) represents revenues from hydrogen, where
denotes the hydrogen output and
h the hydrogen price. Total costs of hydrogen production
are illustrated in (
5) on the basis of annual costs for hydrogen production
in (
6). Here,
denotes the electrical input energy and
the electricity price.
Of course, there are some interdependencies in these equations. Electrical input
depends on the hydrogen price
h and is determined by optimization (based on a linear relationship between
and
) to maximize profits from operation of the electrolyzer. In the technical simulation,
is then determined as a (non-linear) function of
, considering a number of efficiency-curves. The investment costs in a given future year
y (
) are determined by the electrolyzer capacity multiplied by the specific CAPEX value expressed in EUR/kW illustrated in
Table 2. The costs for maintenance in the respective year (
) are compound by the value of 2 % of CAPEX per year.
In the case modeled here, it is assumed that the P2G is connected to the electricity grid (front of meter), feeds into the hydrogen network and does therefore
not consume exclusively renewable electricity. Grid costs (GC) for peak load consumption are illustrated in (
7), where
denotes the peak consumption (maximum load value) within month
m and
denotes the peak-load grid tariff. In Austria, they are formed by the average across all monthly peak consumption values. Grid costs for energy consumption are illustrated in (
8) which are composed of the total energy consumption multiplied by the variable grid tariff
. Specific values for the peak and variable grid tariff are illustrated in
Table 1. The data is based on grid costs for Styria [
40] and refers to the second amendment in 2023. In this version, reductions were granted from the state to the consumer in light of the energy crises in the sharply rising costs for electricity during 2022. For simplification, we have not assumed grid tariffs for injection.
Electricity prices
p consist of historical data for the Day-ahead (15M) EXAA Spot auction from the year 2023 [
41] and were kept constant over the entire lifetime. It is clear that this is a strong assumption, however, electricity prices, in particular the upper and the lower end of the price range is notoriously hard to project with long term fundamental models. The average price is around 100 EUR/MWh with prices ranging between 75 EUR/MWh and 130 EUR/MWh in around 50 % of time.
Table 2 illustrates the cost for investment and operation of the PEMEL. Assumptions were taken from the ‘Technology Catalogue for Renewable Fuels’ from the Danish Energy Agency [
36]. CAPEX and OPEX assumptions are chosen for a system size around 10 MW. The value is actually applied to higher capacity values between 20-70 MW, but the slightly higher CAPEX due to missing scale-effects should reflect a more conservative approach in light of the fact that CAPEX of electrolysis has been underestimated in the past, see [
42]. Actual values for investment in any future year are interpolated from the anchor year values in 2025, 2030 and 2040. For years beyond 2040, the value of 2040 was kept stable.
In a last step, the results of the CBA are compared to the cost of conventional grid enforcement measures. For this purpose, the DSO ‘Energienetze Steiermark’ has estimated the cost of undertaking conventional measures, the values are illustrated in
Table 3. These numbers represent the cost of grid enforcement from the point of view in 2023 to ensure that all anticipated renewable generation can be absorbed by the grid. The numbers can be understood as net present values in real 2023 money terms and can be compared to the lost revenues computed by (
3) and illustrated in
Figure 8. Note that a comparison of grid enforcement costs and the economic value of the additional feed-in (
) is not necessarily in favor of grid enforcement.
All historical economic data (electricity grid enforcement, electricity prices, and grid charges) refer to the year 2023. Therefore, the economic valuation is defined in real money terms for the year 2023. The weighted average cost of capital (WACC) is assumed at 4 % in real terms. The hydrogen price is considered the most uncertain variable in this modeling exercise. Therefore, the hydrogen price is varied in a range between 1 and 10 EUR/kg and all results are illustrated (where relevant) in relation to the hydrogen price.