Simulation of Coupled Power and Gas Systems with Hydrogen-Enriched Natural Gas

: Due to the increasing share of renewable energy sources in the electrical network, the focus on decarbonization has extended into other energy sectors. The gas sector is of special interest because it can offer seasonal storage capacity and additional ﬂexibility to the electricity sector. In this paper, we present a new simulation method designed for hydrogen-enriched natural gas network simulation. It can handle different gas compositions and is thus able to accurately analyze the impact of hydrogen injections into natural gas pipelines. After describing the newly deﬁned simulation method, we demonstrate how the simulation tool can be used to analyze a hydrogen-enriched gas pipeline network. An exemplary co-simulation of coupled power and gas networks shows that hydrogen injections are severely constrained by the gas pipeline network, highlighting the importance and necessity of considering different gas compositions in the simulation.


Introduction
The mitigation of climate change requires the rapid decarbonization of all sectors [1].In recent decades, great progress has been made in the electricity sector by increasing the share of renewable energies [2].The decarbonization process, however, is not exclusively limited to the electrical energy sector.The focus is therefore now shifting to other sectors, for example, industry and transport, in order to reach the Paris goals on climate change.In the past, different energy sectors mainly operated separately or independently.However, with the increasing use of inter-sectoral flexibility options, different energy sectors are set to be coupled more closely and should therefore be analyzed as a whole [3,4].An integrated assessment of the electricity and gas sectors is therefore particularly important.In this paper, we establish a simulation tool framework that is capable of solving hydrogenenriched natural gas flow simulation problems.

The Role of Natural Gas (NG) in Future Energy Systems
Natural gas is an important energy carrier used to generate heat or electricity.In 2019, 10.2% of the total electricity consumption in Germany was covered by natural gas-based generation [5].In 2020, natural gas storage in Germany amounted to around 23.9 billion cubic meters, which corresponds to roughly 250 TWh of energy [6].Due to the increasing share of renewable power generation, long-term storage with power-to-gas (PtG) will most likely become necessary and cost-efficient [7].The gas sector therefore has great potential to provide enormous storage capacity for the power network.Moreover, PtG technologies can provide flexibility to the electrical network [8].The conversion of electrical energy into chemical energy is an essential part of every PtG technology and is achieved using water electrolysis [9].In addition to the impacts on the power sector, hydrogen generated from electrolysis can be a key factor in the decarbonization of heavy industries such as steel or chemistry [10].

Blending Hydrogen into Existing NG Infrastructure
Hydrogen generated using electrolysis can either be directly stored or be transported using pipeline network systems.It can also be further converted into methane using methanation reactions.In order to predominantly utilize the existing gas transmission infrastructure, the feasibility of blending hydrogen into the natural gas network is being intensively investigated [11].Due to different infrastructure states and regulations, the maximum permitted levels of hydrogen concentration in the natural gas transmission system vary from country to country [12].German regulations currently allow for up to 10% hydrogen concentration in the natural gas network.In the future, up to 20% is planned, and the German Technical and Scientific Association for Gas and Water (DVGW) estimates that up to 50% hydrogen concentration is feasible [13].In a recent study, the impact of higher and fluctuating hydrogen concentrations of up to 50% on a variety of industrial combustion systems was investigated [14].Therefore, in order to simulate future gas networks with uncertain hydrogen concentrations, it is important to take gas mixture properties into consideration.

Simulation of Hydrogen-Enriched NG Network
Natural gas is typically a mixture of gases, in which the major component is methane.Depending on sources and locations, the compositions of natural gas can vary significantly.In the study [15], it is stated that only around 23% of industrial plants have a real-time natural gas mixture quality measurement, while more than 50% of them only have weekly or monthly data available.In addition, more than 70% of the customers suffer from poor natural gas supply quality about once a month [15].Natural gas composition must therefore be considered as a factor that has a significant impact on network operation safety.However, in many publications such as [16,17], the properties of natural gas are not listed, which makes it difficult to reproduce their results.When considering the injection of hydrogen into the natural gas network, the impact of gas composition is even greater.
Several research studies were carried out with respect to the steady-state simulation methods for hydrogen-enriched natural gas networks.For example, Abdolahi et al. used several equations of state (EOS) to model the natural gas mixture, whereas the heating values were not considered [18].In the work by Giulio et al., an analysis was performed to evaluate the gas heating values with respect to the hydrogen injection into the natural gas pipelines, while the impact of the gas mixture properties on the gas flow calculation was not the main focus [19].Pellegrino et al. established a simulation framework for the hydrogen-enriched natural gas network using the virial EOS (cf.[20]), but the height difference between the pipeline inlet and outlet was not considered [21].
In addition, there are a number of tools available that can be used to analyze the hydrogen-enriched natural gas network.PSS ® SINCAL (cf.[22]), SAInt (cf.[23]), and MYNTS (cf.[24]), for example, are well established tools that are capable to run complex simulations for gas pipeline networks.Nevertheless, they are closed-source software and therefore hard to extend with functionalities for the specific requirements of analysis.TransiEnt is a Modelica library that is capable of steady-state and dynamic simulations of coupled-power, gas, and heat networks [25].However, it is currently only supported in Dymola, which is a commercial software development environment for Modelica.There is a relatively new Python package-so-called pandapipes-that can handle different gas mixture properties for gas network simulation.The nodal difference of gas mixture compositions, however, has not yet been considered [26].

Modeling of Gas Pipelines
A gas network typically consists of pipelines, compressor stations, valves, and regulator stations.For simulation purposes, a number of additional virtual components are necessary, for example, short pipes and fictitious resistances.However, the working principle of these components is similar to that of normal pipelines.As the pipeline is the dominant component in a gas network, the modeling of pipelines is explained in detail in this section.
Steady-state gas flows can be calculated using pipeline equations.There are several pipeline equations that can be used to calculate volumetric gas flow rates, although they are mostly derived from the isothermal Euler equation (Equation (1), cf.[27,28]).
After calculating the integral and substituting the constant parameters with real values, the steady-state volumetric flow rate in a pipeline can be calculated as in Equation (2) (cf.[16,29]).Unless stated otherwise, the gas volumetric flow rates will be converted to the ones under standard reference conditions of 15 °C (288.15K) and 1 bar (101,325 Pa), which correspond to T st and P st in the pipeline equation [30].The unit of the standardized volumetric flow rate is sm 3 /s.
where C = π R 16M air is a constant, which is around 13.29, D is the pipe diameter, η is the pipe efficiency, f is the friction factor, G is the gas specific gravity, L is the pipe length, P st is the reference pressure, P i is the inlet pressure, P j is the outlet pressure, Q is the volumetric flow rate, T a is the average temperature, T st is the reference temperature, and Z is the compressibility factor.In Equation ( 2), E represents the potential where H i and H j are the inlet and outlet height, respectively.The average pressure P a and average temperature T a can be calculated using the following formulas (cf.[16]): The friction factor is a very important variable for calculating gas flow rates.There are various ways to calculate the pipeline friction factor that involve different levels of computational complexity and accuracy.However, the friction factor does not change significantly in the fully turbulent zone [31], which is also the case in this paper.Because the actual friction factor needs to be calculated or calibrated using pipeline efficiency, this paper uses the simplest method for pipeline friction, which is only related to the pipeline diameter D (cf.[16,29]).
η is the efficiency of the pipe to convert the theoretical friction factor into an actual one, taking other sources of friction into account, for example, valves and tees.The aging of the pipeline and corrosion or rust inside also contribute additional friction to the gas transmission system.Therefore, the real gas flow rate in a pipe is generally lower than the one calculated by flow equations where η = 1.To account for such extra flow reductions, the efficiency factor η is usually chosen between 0.75 and 0.95, while experience shows that for an old pipeline it can be reduced to lower than 0.7 [32].However, Mohitpour et al. suggested η values between 0.92 and 0.97, which are much higher than the ones mentioned above (cf.[33]).In principle, η should be chosen according to actual gas network pipelines.In the remaining part of this paper, all pipeline friction factors are set to 0.85, which is a good assumption for a common pipeline operating status.
By reviewing Equation ( 2), it can be seen that most constants are preset pipeline parameters.It therefore can be simplified further using algebraic transformations, as shown below in Equation ( 7): (7) where C pipe reflects constants and pipeline parameters Typically, as shown in Equation ( 2), pipeline network simulations are based on volumetric flow rate conservation, which means the sum of all gas volumetric flows injected into or flowing out of one node equals zero.However, considering variant gas compositions in pipelines makes this assumption invalid.To deal with this issue, mass flow conservation is adopted in this work.

Thermal Formulation
Temperature is another important state variable in a gas network simulation.To calculate the temperature profile alongside a gas pipeline, an extra equation based on the first law of thermodynamics is needed (Equation (9), cf.[27,28]): where Q m is the gas mass flow rate, h is its specific enthalpy, v is the velocity of gas flow, U l is the heat transfer coefficient of the pipeline, g is the gravitational acceleration, and sin θ = When considering a steady-state simulation, the change of velocity can be ignored ( dv dx = 0).By further assuming the gas enthalpy as a function of pressure and temperature, the change of enthalpy dh can be rewritten as Equation ( 10) [27]: , Equation ( 9) can be rewritten as Equation (11).
To simplify the calculation, here the potential energy term is ignored.By substituting dp dx with the help of Equations ( 1) and ( 2), Equation ( 11) can be written in the form of Equation ( 12) (cf.[16,21]): After calculating the integral and simplification, the gas pipeline outlet temperature can then be calculated using the following equation: where

Calculation of Gas Mixture Properties
As we can see in Equation ( 7), there are two properties that are directly related to the gas composition: the gas specific gravity G and the compressibility factor Z. For the purpose of simulating the gas network, these values are typically obtained using empirical models.In the case of the compressibility factor, for example, Papay's equation [17,21] (Equation ( 14)) or AGA [34] (Equation ( 15)) are used, where p r = p/p c and T r = T/T c .These approaches are suitable when calculating the compressibility factor of typical natural gas.However, the critical point conditions (p c and T c ) used in the formula are constants, assuming that the gas mixture properties do not change during network operation.However, this assumption is no longer valid when considering hydrogen blended into the natural gas network.Therefore, an automated calculation with variant gas composition is not possible using these empirical models.
Now consider the gas specific gravity and heating value of the gas mixture.According to the DVGW technical regulation [35], natural gas is defined as H-gas and L-gas with respect to its composition.The properties of these two types of gases are shown in Table 1.It can be seen that heating values and specific gravity vary considerably, which makes an accurate calculation more difficult.Therefore, a comprehensive gas mixture property calculation is used in this paper.
Table 1.Properties of natural gas in Germany [30].

Gas Properties Units L-Gas H-Gas
Heating value MJ/sm 3 To calculate the gas mixture properties, the "thermo" python package is used [36].The gas mixture properties are calculated using the Peng-Robinson EOS (PREOS) for a mixture of any number of compounds.The mathematical formulations can be found in the package repository or in the referenced literature [37,38].
The gas mixture heating values are calculated using another package named "Cantera" [39], based on the combustion data GRI-Mech 3.0 [40].Limited to the available gas species data in this source, hydrocarbon species with over 3 carbon atoms (e.g., butane, pentane) are considered as methane in the heating value calculation.
It is first used to balance all the chemical reaction equations of complete combustion based on the composition of the gas mixture.The corresponding heating value (HV) can be subsequently calculated using Equation ( 16) (cf.[41]), where R is the set of all reactants (all gas species in the gas mixture and O 2 ) and P represents all products.H stands for the enthalpy, and x is the mole fraction of a single gas species in R or P. In this paper, the higher heating value (HHV) of the gas mixture is used, which assumes all water content in the end product is in a liquid state.
The natural gas composition used in this work is listed in Table A1 [42].

Solution Flow
In this tool, a gas network is modeled as a set of nodes and pipelines.The network nodes are classified into three types: reference nodes, supply nodes, and demand nodes.For the reference nodes, pressures are known and gas flows need to be calculated to balance network demand.For supply nodes and demand nodes, the gas flows are known and pressures need to be calculated.The passive sign convention system is used in this tool, with negative flows representing supply and positive ones representing consumption.
The simulation method presented in this paper combines conventional pipeline equations with comprehensive modeling of gas mixture and utilizes the Newton-Raphson (NR) method to solve steady-state gas flow calculation problems iteratively.As can be seen in Figure 1, the first step of the NR method is the initialization of network variables, which in this case are the initial estimates of unknown nodal pressures.In an electrical power flow simulation, a "flat-start" initialization is typically applied.In contrast, the pressures at both ends of a pipeline must not be the same; otherwise, this will result in zero flow in the pipeline and thus a poorly conditioned Jacobian matrix.From experience, pressures at pipeline outlets are initialized as 0.98 of those at inlets.To avoid the duplicated assignment of initial estimates, the initialization algorithm assigns pipeline outlet pressures starting from reference nodes until each node has a pressure value.Based on the initial estimates, the physical properties of the gas mixture inside pipelines can be calculated.For each iteration step, the Jacobian matrix is updated, which is later used to calculate the state variables.Afterward, the new gas mixture composition and its properties are calculated.Since the approach presented here is aimed at solving the static flow problem, flows entering a node with different compositions are assumed to be completely mixed at the node.Therefore, the gas mixture composition in a pipeline is always considered to be the same as the one at its inlet.As long as the error after one iteration step is bigger than the set tolerance, the program iterates over to set pressures to the nodes and simultaneously updates the gas mixture properties to be used for the pipelines.If the error remains within the tolerance, the simulation has converged and results are saved.

Study Case 1: Impacts of Different Calculation Methods on a Single Pipeline with Variable Gas Mixture Composition
Since the pipeline is the dominant component in a gas network, in the first study case the simplest case is considered, which is the analysis of different gas compositions in one single pipeline.The main purpose of this study case is to stress the importance of taking the gas composition into account when performing a pipeline flow calculation with respect to hydrogen blended into natural gas.To achieve this, PREOS is used in this study to calculate the gas mixture properties, which is proved to be one of the most accurate methods to model the natural gas mixture [43].Since the assumption made for Equation ( 6) is only valid for systems with higher pressures, only high-pressure gas pipelines are considered in this section.In Germany, the pressure level of high-pressure gas networks is defined as an operational range between 1 and 100 bar [44].Therefore, in this paper, the analyzed pressure range is also set between 1 and 100 bar.As described in previous sections, multiple pipeline flow equations are available.A comparison of the most popular methods is shown in Figure 2a to give an overview of these different methods.It can be noted that at lower pressures (from 1 to 10 bar), the calculation results of different methods do not differ significantly.Within a higher pressure range (from 10 to 100 bar), the Panhandle B method tends to overestimate the flow rates, while the Weymouth method underestimates them.The flow rates calculated with the method used in this tool lie in the middle range of all analyzed methods, both in the lower pressure range and in the higher pressure range.
Typically, simulations of gas pipeline networks are based on constant gas properties, which is no longer accurate when the properties change significantly, for example, when the operation conditions deviate strongly from the nominal value or when complex systems with different gas mixtures are simulated.As shown in Figure 2b, a comparison between the Papay's equation and the PREOS is made.The calculation results from these two alternatives appear to be quite similar.However, at a higher pressure, the calculated flow rates with Papay's equation, which is based on constant critical conditions, are around 2% lower than those using PREOS which is a more accurate method for calculating gas properties.The pipeline transmission is generally restricted by the minimum and maximum permitted volumetric flow rates and the maximum permitted pressure.Therefore, it is important to investigate how the flow rate and pressure change when different gas compositions are applied to the pipeline.Now consider a binary gas mixture of CH 4 and H 2 .The simplified pipeline equation (Equation ( 7)) shows that the gas flow rate is dependent on the compressibility factor Z and the specific gravity G of the gas mixture.As shown in Figure 3a, the gas compressibility factor increases in a non-linear fashion with increasing hydrogen concentration.If we consider a 20% hydrogen concentration in the gas mixture, the corresponding compressibility factor is around 2% greater than the one calculated with the linear mole fraction model.When using Papay's equation, it is difficult to calculate the T c and P c of the gas mixture; thus, a simple linear mole fraction method is used.As shown in Figure 3a, at low hydrogen concentrations, this method works well, but when the hydrogen concentration is increased further, the error also becomes bigger.Obuba et al. showed that Papay's method tends to underestimate the Z-factor [45], an observation that also corresponds to the results shown in this figure .Turning now to the gas specific gravity G.As shown in Figure 3b, the gas-specific gravity is linearly related to the hydrogen concentration, and the calculation results obtained from using PREOS and the linear mole fraction model are nearly identical.After analyzing the gas properties of gas mixtures with different hydrogen concentrations, the pipeline flow calculation results are shown here.In the practical gas network operation, the energy requirement of the demand side has to be met.Therefore, the actual gas flow rate may vary due to variations in the heating values of the gas mixture.In this study case, a constant energy demand at the pipeline outlet is assumed in order to analyze the pressure drops throughout the pipeline with respect to different hydrogen concentrations.
As shown in Figure 4, when considering a certain amount of energy demand, the pressure drop in the pipeline increases with increasing hydrogen concentration.At around 90% hydrogen concentration, however, there is a turning point in the pressure drop curve.Looking back to Equation ( 7) and Figure 3a, the reason for this turning point is that the increase of the compressibility factor Z slows down when the hydrogen concentration increases.Besides the nodal pressure and the flow rate inside pipelines, the temperature of the transmitted gas mixture is also an important variable for gas network simulation.In Figure 5a, simulations are run to illustrate the temperature profile inside a pipeline with respect to the different inlet temperatures.It can be noted that the gas mixture temperature tends to be the same as its ambient temperature after a long-distance transmission.In Figure 5b, several simulations are run assuming different hydrogen concentration rates inside the pipeline.In the case of natural gas, the gas mixture temperature is assumed to be the same as the ambient temperature if the pipeline is longer than 25 km.However, with the increase of hydrogen concentration inside the pipe, the descent rate of the temperature becomes lower and the temperature difference between the pipe inlet and outlet has to be taken into consideration in certain cases.
Gas transmission in the pipeline network is based on pressure control.The results presented in this section show that in order to keep the same energy flow in one pipeline, gas flow rates and pressure drops also change when the hydrogen concentration changes in the pipeline.
Since the volumetric flow rate is the most important variable to be analyzed, only the volumetric flow will be analyzed and shown in the rest of the paper.

Study Case 2: A Simple Gas Network with Hydrogen Injection
In the previous study case, simulations are performed on a simple pipeline model to illustrate the importance and necessity of taking gas mixture properties into account as part of a hydrogen-enriched natural gas network simulation.Considering the possibility of blending hydrogen into the natural gas network, if the hydrogen injection occurs at the demand side, the hydrogen concentration could be gradually increased and therefore cause risky conditions on a local level.In this section, a number of simulations are performed on a hydrogen-enriched natural gas network, to analyze the impacts caused by different hydrogen concentrations in the network or nodal hydrogen injections into the network.
To reproduce the network behavior with the utmost accuracy, a compatible network size for both gas and power grids should be considered.In this work, the gas grid is synthetically generated based on the CIGRE high-voltage transmission benchmark grid [46], ignoring some buses and connections.
The synthetically generated simple gas network model consists of 8 nodes and 8 branches as shown in Figure 6.Node 1 and node 2 are set as reference nodes, where the nodal pressure values are known.At these nodes, the pressures remain constant, while the flow rates are variable to be able to balance the network demands.The other nodes are defined as demand nodes, where the flow rates are predefined.The demand can be also converted into energy demand using the gas mixture heating values.As shown in Figure 7, when nodal energy demands stay the same, the overall pressure at all demand nodes will be lower when the hydrogen concentration is higher.However, from 80% to 100% hydrogen concentration, the nodal pressures increase slightly, which corresponds to the results shown in Figure 4.
Water electrolysis is one of the most promising electrical flexibility options so far as it can generate hydrogen by electrolyzing water.To analyze the impacts of electrolyzer operation on the natural gas network, it is assumed that a known amount of hydrogen is fed into the natural gas network at the network nodes.In this study case, the distribution of injected hydrogen in the network is analyzed to see if the maximum permitted hydrogen concentration in the pipeline is exceeded.Now assume there are three electrolyzers, which are accordingly located at nodes 3, 4, and 5.The energy flow rates of hydrogen blending into the natural gas grid are assumed as 35 MW, which corresponds to around 2.9 sm 3 /s.In this case, different gas injections are considered to be fully mixed at conjunction nodes.The corresponding simulation results are shown in Figure 8.As shown in the figure, a relatively high hydrogen concentration of around 29.6% occurs in pipeline 3.In pipelines 2, 4, and 6, the hydrogen concentration also reaches around 13%, which is already critical with respect to current technical regulations.It should be noted, however, that although the assumed hydrogen injection rate is relatively high given the current status of real-world implementation, the total energy delivered with hydrogen covers only around 1.88% of the total energy demand in this system.Therefore, when the total hydrogen share in the gas transmission system is increased, it is important to monitor the hydrogen concentration in pipelines in order to ensure a secure operation.

Study Case 3: A Simple Example of a Simulation of Coupled Power and Gas Networks
A number of studies analyzed the technical and economical performances of water electrolyzers considering different PV and wind penetration levels [47,48].However, these studies focused mainly on the electrolyzer and its power consumption rather than on the grid perspective.In this part, a simple analysis is performed by coupling the power grid (Figure 9) and gas grid (Figure 6).To this end, a quasi-dynamic simulation is performed based on a time-series simulation.The energy demands in both networks are calculated using standard demand profiles, which are taken from [49,50] and assume an annual demand of 400 MWh and 600 MWh, based on the scenario defined in the original power system model [46].To include the generation of renewable energy, weather data from 2015 are used [51].The wind and PV generation profiles are calculated using the abovementioned weather data and formulas taken from [52,53].The total installed generation capacity of wind and PV generation units is set at 500 MW with a wind and PV ratio of 50-50%.Furthermore, renewable energy generation is considered to be equally distributed at each node.In Figure 10, the PV and wind generation are plotted as box plots [54] to show their seasonal behavior.The orange lines in the middle represent the median value of the power generation of each month.The lower and upper bounds of the boxes represent the lower and upper quartiles of monthly generation.The whiskers extend from the box and cover the major range of the data.The outliers refer to the data past the end of the whiskers, which are the peaks and valleys of power generation.It can be seen that although PV has a peak generation in the summer and the maximum total generation of renewable energy usually occurs in warmer seasons, the relatively high continuous generation of renewable energy occurs during winter.Therefore, two scenarios are assumed based on the above-mentioned renewable generation, which correspond to June and December, meaning that the PV-rich and wind-rich seasons can be considered.
To investigate the impact of coupling both networks, the combination of three water electrolyzers and one gas-fired power plant is analyzed (see Figures 6 and 9).The water electrolyzers are considered to be installed at nodes 3, 4, and 5, each with a maximum power consumption rate of 20 MW and an efficiency of 80%.A simple scheduling method is applied to ensure that the electrolyzers are only operated at times when the generation of renewable energies exceeds its nodal power demand.The gas-fired power plant is assumed at node 8, which consumes gas and produces electrical power with an efficiency of 60%.Here, instead of changing settings on the power grid side, a constant volumetric flow rate of gas consumption is assumed so that the impact of the time-variant hydrogen concentration in the natural gas grid can be illustrated.Two simulations are run separately for June and December 2015 with a sampling time of one hour.The results are shown in Figures 11 and 12.By comparing the results from Figure 11a,b, it is clear that the hydrogen concentrations in June are much higher than those in December.This result can be explained by the fact that the gas demand is usually rather low in summer while generation spikes of renewable energies are rather common.Because natural gas is mostly used for space heating purposes, the gas demand increases during colder days.Although more hydrogen is produced in the winter, the energy supplied with hydrogen does not play a major role.
In Figure 12, the power generation of the gas-fired power plant at node 8 is analyzed.It should be noted that the power generation in June is not stable, which can vary up to 20%.This is caused by the high concentration of hydrogen that occurs in pipeline 7, which directly supplies node 8.This further indicates that directly blending hydrogen into the natural gas network is theoretically impractical during summer.Although the power generation in December is relatively more stable, a difference of over 5% is still possible.Therefore, to ensure a secure operation of the power grid coupled with the hydrogen-enriched natural gas network, the heating values of the gas mixture should be taken into consideration.

Conclusions
The gas sector can provide seasonal storage capacity and additional flexibility to the power sector in an integrated power and gas system.To accurately analyze the impacts caused by coupling power and gas networks, appropriate modeling and simulation approaches are needed.In this paper, a simulation method and a corresponding tool are presented that are capable to simulate hydrogen-enriched natural gas networks.To achieve this, the Peng-Robinson equation of states is applied to calculate the gas mixture properties.
In the previous section, several study cases were developed and simulated to demonstrate the importance and necessity of modeling variant gas compositions.First, different modeling approaches were analyzed by applying them to a single pipeline.By analyzing the results, it can be concluded that hydrogen injection can pose challenges to the gas pipeline network.The gas mixture composition has a significant impact on the static gas flow calculation, which can represent the long-term gas network operation.With a potentially high penetration of hydrogen in the natural gas network, the gas volumetric flow rate and nodal pressure might vary due to the different hydrogen concentrations in gas pipelines.In addition, when the hydrogen concentration in the natural gas network system increases, the heating value of the transported gas mixture decreases.Therefore, in order to maintain a sufficient energy supply, the volumetric flow rates have to be increased, which in turn leads to increasing pressure drops throughout the pipelines Therefore, local pipeline congestion may occur in the future, which should be analyzed in advance during network planning.
The method is then applied to a simple network, where nodal hydrogen injections take place.The simulation results further confirm the findings from the previous study case.Moreover, the presented method and tool can also be used to calculate hydrogen concentration in pipelines, which is an important state variable in the analysis of hydrogen-enriched natural gas networks.Furthermore, the simulation results also provide information about the heating value of the gas mixture, which is very helpful for real-world applications, where very little real-time measurement data are currently available.
Finally, a coupled power and gas network is simulated and analyzed.This study case shows that the power and gas demand varies greatly on a seasonal level.The feasibility of blending hydrogen directly into the natural gas grid is therefore investigated using the developed simulation tool.Due to the relatively lower gas demand in summer, the permitted amount of synthetic hydrogen in the gas grid is rather limited.Therefore, a more reasonable option is to store the synthetically generated hydrogen separately or together with natural gas in warmer seasons.In winter, the gas demand is higher due to the greater heating demand.Blending hydrogen into the gas grid is therefore more feasible.

Figure 1 .
Figure 1.Flow chart of the solving process.

Figure 2 .
Figure 2. (a) Comparison between different pipeline flow rate calculation methods.(b) Comparison of flow rate calculation with constant gas properties and those calculated using PREOS.

Figure 3 .
Figure 3. (a) Gas mixture compressibility in the relation to the injection of hydrogen.(b) Gas mixture specific gravity in relation to the injection of hydrogen.

Figure 4 .
Figure 4. Pressure drop throughout a pipeline in relation to different hydrogen concentrations.

Figure 5 .
Figure 5. (a) Temperature profile throughout a gas pipe with different pipeline inlet temperatures.(b) Temperature profile throughout a gas pipe with different hydrogen concentrations.

Figure 7 .
Figure 7. Pressure at demand nodes with respect to different hydrogen concentrations in the network.

Figure 8 .
Figure 8. Hydrogen concentration in pipelines with respect to the injection of hydrogen at three nodes.

Figure 10 .
Figure 10.(a) Relative PV power generation in relation to installed generation capacity.(b) Relative wind power generation in relation to installed generation capacity.

Figure 11 .
Figure 11.(a) Hydrogen concentration in pipelines in June.(b) Hydrogen concentration in pipelines in December.

Figure 12 .
Figure 12.Gas-fired power plant generation with constant volumetric flow consumption.