Regular Article
Effects of hydrogen blending on hydraulic and thermal characteristics of natural gas pipeline and pipe network
^{1}
School of Mechanical Engineering, Beijing Key Laboratory of Pipeline Critical Technology and Equipment for Deepwater Oil and Gas Development, Beijing Institute of Petrochemical Technology, Beijing 102617, PR China
^{2}
Beijing Key Laboratory of Process Fluid Filtration and Separation, College of Mechanical and Transportation Engineering, China University of Petroleum (Beijing), Beijing 102249, PR China
^{*} Corresponding author: lijingfa@bipt.edu.cn
Received:
20
July
2021
Accepted:
20
September
2021
Blending a fraction of hydrogen into the natural gas pipeline or urban pipe network is an efficient approach for hydrogen delivery. In this paper, the mathematical model of HydrogenBlended Natural Gas (HBNG) transportation is established, and the influences of hydrogen blending on hydraulic and thermal characteristics of natural gas pipeline and pipe network are numerically investigated. The impact of hydrogen blending ratio on the performance of centrifugal compressor and the operating point for joint operation of pipeline and compressor is discussed. Results illustrate that compared with natural gas without hydrogen, the hydrogen blending can reduce the pipeline friction resistance and increase the volume flow rate. However, due to the lower volumetric calorific value of HBNG, the energy flow rate actually decreases under the same transportation condition. Meanwhile, the temperature drop along the pipeline slows down due to the blended hydrogen. The performance degradation of centrifugal compressor occurs with the increasing hydrogen blending ratio, and the operating point for joint operation of pipeline and centrifugal compressor moves to the direction of higher volume flow rate and lower pressure. This study is expected to shed a light on the hydrogen delivery by natural gas pipelines and pipe networks.
© H. Zhang et al., published by IFP Energies nouvelles, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Air pollution and global warming are two major environmental threats to human beings nowadays, which are tightly related to the overuse of fossil fuels. The development and utilization of lowcarbonintensity fuels (such as natural gas, hydrogen, etc.) are of great significance to the energy conservation, emission reduction and energy structure optimization [1–3]. As a kind of secondary energy, hydrogen has been recognized as a clean energy for its zerocarbon merit as only water is produced in its combustion. However, currently the utilization of hydrogen is heavily limited by the insufficient development of hydrogen transportation and distribution facilities. Pipeline transport is an economic and efficient way for longdistance and largescale natural gas transmission, and the relevant technologies have been maturely developed in past years. Literature studies have preliminarily shown that blending a certain amount of hydrogen into natural gas and using the existing natural gas pipelines or urban pipe networks for transportation is an efficient method for hydrogen delivery [4–7]. It is also noted that the transportation of hydrogen through natural gas pipelines is a lowcost option and is beneficial for emission reduction [8–10]. Therefore, it is meaningful to study the transportation of hydrogen through natural gas pipelines and pipe networks.
Due to the substantial difference in physical and thermal properties of hydrogen and natural gas, the influences of hydrogen blending on hydraulic and thermal characteristics of natural gas pipeline and pipe network should be fully considered in the transportation technology. At present, researchers at home and abroad have devoted effort to exploring the effects of hydrogen blending on natural gas pipeline transmission [11–13]. For example, Tabkhi et al. [14] applied an isothermal model of natural gas pipeline transportation to investigate the impact of hydrogen blending on hydraulic conditions of a longdistance pipeline. It was found adding hydrogen into natural gas significantly reduces the pipeline transmission capacity, and the maximum allowable hydrogen blending ratio is about 6% (mass fraction) without changing any facilities in the pipeline system. Elaoud and HadjTaïeb [15] and Elaoud et al. [16] assessed the pressure variation of a smallscale ringed natural gas pipe network associated with hydrogen blending at steady and transient states by using isothermal models. Results showed the pressure oscillation of pipe network turns out to be more and more obvious with the increasing hydrogen blending ratio at transient states. Guandalini et al. [17] analyzed the effect of hydrogen blending on hydraulic conditions of a natural gas pipeline, and found the impact of low hydrogen blending ratio (5% volume fraction) on pressure drop (0.1%) is negligible, but the mixed hydrogen was able to increase the volume flow rate of the pipeline.
Later, Wang et al. [18] investigated the influence of hydrogen blending on both hydraulic and thermal characteristics of a longdistance pipeline as well as on the performance of a centrifugal compressor. The effect of temperature change on pipeline transportation capacity was taken into account. It was found the pipeline transportation capacity decreases with the increase of temperature, and the volume flow rate of the pipeline increases with the increasing hydrogen blending ratio, indicating the hydrogen blending reduces the pipeline friction resistance. During the joint operation of the pipeline and centrifugal compressor, both the pressure and flow rate corresponding to the operating point were reduced with the increase of hydrogen blending ratio. Hafsi et al. [19, 20] studied the steady and transient flow of natural gashydrogen mixture in a simple triangular pipe network by using the isothermal model ignoring the temperature variation. Results demonstrated that the pressure drop at steady state is enhanced by hydrogen blending, and the flow velocity in the pipe increases significantly as well. In addition, the overpressure at transient condition was proportional to the mass fraction of mixed hydrogen. Wu [21] conducted an investigation on the variation of node flow rate and pressure of Chongqing natural gas pipe networks when 10 hydrogen blending sources were considered. It was revealed that the pressure at each node shows a decline trend with the increasing hydrogen blending ratio, and the nodes with larger flow rate suffer a higher pressure drop. Besides, the pressure drop at the node close to the gas source was greater than that far away from the gas source. The research of Uilhoorn [22] demonstrated when the volume flow rate is fixed, the blended hydrogen can alleviate the pressure drop and temperature drop of the gas pipeline, and the compressor power is significantly reduced as well; however, when the energy flow rate is fixed, the pressure drop increases and the temperature drop decreases, as well as the compressor power raises significantly.
To sum up, the previous researches on the hydraulic condition of HydrogenBlended Natural Gas (HBNG) pipeline mainly focused on a single pipeline or a single pipe network, lacking of comprehensive consideration covering both the pipeline and pipe network. The study on the influence of hydrogen blending on the performance of centrifugal compressor is still insufficient. In addition, almost all of these researches are based on the isothermal model of pipeline transportation, the effect of temperature change on the hydraulic characteristics of HBNG transportation is rarely considered. At the same time, the impact of hydrogen blending on the thermal condition of natural gas pipeline is still needed to be further clarified. To fill the gap, in this study the influences of hydrogen blending on hydraulic and thermal characteristics of a longdistance natural gas pipeline with the centrifugal compressor and a ringed multiloop pipe network are systematically investigated. The effects of hydrogen blending ratio on the performance of the centrifugal compressor and the joint operating point of pipeline and compressor are discussed. This study is expected to provide valuable reference for the development of pipeline transportation technology of HBNG.
The structure of this work is organized as follows: Section 2 introduces the mathematical model describing the hydraulic and thermal conditions of natural gas pipelines and pipe networks; Section 3 presents the numerical method for hydraulic and thermal simulations of gas pipeline transportation based on the linearization technique; Section 4 studies the influences of hydrogen blending on the hydraulic and thermal characteristics of a longdistance pipeline and the performance of a centrifugal compressor; Section 5 discusses the impact of hydrogen blending on the hydraulic and thermal characteristics of a ringshaped multiloop pipe network; Concluding remarks of this work are summarized in Section 6.
2 Mathematical model for HBNG pipeline transportation
The flow and heat transfer of HBNG in pipeline obey “Laws of mass conservation, momentum conservation, and energy conservation”. Correspondingly, the hydraulic and thermal governing equations of HBNG pipeline transportation can be expressed as follows [23],
(3)where A is the crosssection area of the pipeline, m^{2}; p is the pressure, Pa; ρ is the gas density, kg/m^{3}; w is the gas velocity, m/s; λ is the coefficient of hydraulic friction; g is the gravitational acceleration, g = 9.81 m/s^{2}; θ is the inclination angle of the pipeline, rad; c_{v} is the specific heat capacity at constant volume, J/(kg·°C); K is the total heat transfer coefficient of the pipeline, W/(m^{2}·°C); T_{g} is the gas temperature, °C; T_{g} is the temperature of ambient media around the pipeline, °C; D is the outer diameter of the pipeline, m.
The state parameters (PVT relationship) of HBNG can be calculated by the gas Equation Of State (EOS), such as the SRKEOS [24], PREOS [25], BWRSEOS [26], and GERG2008EOS [27], etc. In this study, the widely used BWRSEOS in natural gas engineering is applied,
(4)where p is the gas pressure, kPa; T is the gas temperature, K; ρ is the gas density, kmol/m^{3}; R is the gas constant, R = 8.3145 kJ/(kmol·K). A_{0}, B_{0}, C_{0}, D_{0}, E_{0}, a, b, c, d, α and γ represent the model parameters in BWRSEOS. For a pure gas component, these 11 parameters can be calculated according to the critical pressure, critical density and acentric factor. For HBNG mixture, the 11 parameters can be calculated based on that of each pure gas component by the mixing rules.
Equations (1)–(4) are the basic mathematical models for the HBNG in pipelines and pipe networks. Compared with the longdistance pipeline, the pipe network is composed of multiple pipelines and nodes. At each junction node, the flow rate, pressure, and outlet temperature should satisfy following relationships,
(7)where N represents the number of junction nodes; m is the mass flow rate, kg/s; c_{p} is the specific heat capacity at constant pressure, J/(kg·°C); the subscripts “in” and “out” denote the inlet and outlet of junction nodes, respectively.
In a HBNG pipeline system, the centrifugal compressor is an important equipment providing the energy head for transportation, and the fitting of compressor performance curve is needed in the calculation. In this work, the following relations are applied to fit the pressure ratio and polytropic efficiency of a centrifugal compressor,
(9)where n denotes the rotational speed, rpm; , , , , , , , , , and are the parameters to be fitted in the pressure ratio curve; , , , , , , , , , , and are the parameters to be determined in the polytropic efficiency curve.
3 Numerical methods
The hydraulic and thermal governing equations (1)–(3) are typical nonlinear hyperbolic equations. For the convenience of calculation, first the linearization technique [28] is used to reformulate equations (1)–(3) into the following linear form,
(10)where the expression of U, B(U) and F(U) is listed in Table 1, in which the first row denotes the hydraulic equation and the second row stands for the thermal equation.
The coefficients of linear equation system.
To discretize the problem in space and time, the spatial step is set as Δx and the time step is set as Δt. In this study, the uniform mesh presented in Figure 1 is employed to discretize the pipeline into N segments, and the nodes are labeled as 1, 2, …, N + 1. The nodes 1 and N + 1 denote the boundary of the pipeline.
Fig. 1 Discretization in space and time. 
In this work, the finite difference method is applied to discretize equation (10) on the mesh shown in Figure 1. Specifically, the unsteady term of hydraulic governing equation is discretized by the forward difference scheme, and the convection term is discretized by the secondorder central difference scheme. While for the thermal governing equation, the convection term is discretized by the firstorder upwind scheme. Readers of interest can refer to [28] for more detail. The discrete hydraulic governing equation can be arranged as follows,
(11)where , , , and I is a 2 × 2 identity matrix, while the matrix G is given by,
Similarly, the discrete thermal governing equation can be written as below [29],
It can be observed that the coefficient matrix of linear equations (11) and (13) is in a Tridiagonal form and is a sparse matrix. Therefore, the Tridiagonal Matrix Algorithm (TDMA) is adopted to efficiently solve the linear systems (11) and (13). Besides, the acceleration techniques can be applied to speed up the solution of linear systems, such as the GraphicalProcessingUnit (GPU)accelerated algorithm [30].
4 Effects of hydrogen blending on longdistance natural gas pipeline
4.1 Pipeline parameters
Taking a longdistance pipeline as an example, the influence of hydrogen blending on hydraulic and thermal characteristics of the pipeline transportation system is investigated, such as the flow rate, transportation capacity, outlet pressure, pressure drop, and temperature drop, etc. Besides, the impact of hydrogen blending on the performance of centrifugal compressor is also explored. The detailed parameters and sketch map of the pipeline system are shown in Table 2 and Figure 2, and the composition of natural gas studied in this work is presented in Table 3.
Fig. 2 Sketch map of the studied longdistance pipeline. 
Parameters of the longdistance natural gas pipeline.
Composition of natural gas.
To quantitatively analyze the influence of hydrogen blending on hydraulic and thermal states of the longdistance pipeline, a proper comparison standard should be established beforehand. Due to the large difference in calorific value between hydrogen and natural gas, mixing hydrogen into natural gas changes the volumetric calorific value of natural gas, thus the volume flow rate cannot truly reflect the gas transmission capacity of the pipeline. Therefore, it is more reasonable to use energy flow rate than volume flow rate to characterize the gas transmission capacity, which is defined as,
(14)where Φ is the energy flow rate, MJ/h; Q is the volume flow rate, m^{3}/h; H is the higher calorific value of natural gas mixture, MJ/m^{3}.
In hydraulicthermal simulations, the initial and boundary conditions should be supplemented to close the governing equations. In this study, the following two boundary settings are adopted: (1) control the pressure at the start and the end of the pipeline; (2) control the energy flow rate and the pressure at the end of the pipeline. It is worth pointing out the literature study indicates that when the upper limit of hydrogen blending ratio is 20%–30%, the HBNG can be used normally without modifying domestic gas cooking appliances. Thus the concerned hydrogen blending ratio (vol.%) in this work is set as 0%, 5%, 10%, 15%, 20% and 30%, respectively. In addition, owning to that the soil temperature at pipeline buried depth is different in different seasons, the temperature of 8.4 °C in winter, 25.2 °C in summer, and the annual average temperature of 15 °C, are considered in simulations.
4.2 Results and discussion
Firstly, the influence of hydrogen blending on gas transmission capacity of a longdistance pipeline is studied. For the convenience of analysis, the effect of hydrogen blending on compressor performance is discussed in Section 4.3, only the pipeline segment 2–3 in Figure 2 is focused on here. Assuming that the outlet pressure and outlet temperature of the compressor are 9.6 MPa and 50 °C, and the pressure at the pipeline end is 7.4 MPa.
Table 4 shows the volume flow rate (Q) of the longdistance pipeline at different hydrogen blending ratios and soil temperatures at pipeline buried depth. It is found if pressures at the start and the end of pipeline 2–3 are fixed, the gas volume flow rate gradually increases with the increasing hydrogen blending ratio whether in winter or in summer. When the hydrogen blending ratio reaches 30%, the volume flow rate of the HBNG pipeline is increased by 11.13% on average compared with the pipeline without hydrogen, indicating the hydrogen blending has certain influences on gas transmission capacity. The reason can be attributed to properties difference between hydrogen and natural gas, which leads to the reduction of friction resistance along the pipeline with mixed hydrogen, and then the gas volume flow rate increases when the pressures at the start and the end of the pipeline remain unchanged. Table 4 also illustrates that the effect of hydrogen blending on gas volume flow rate of the longdistance pipeline is almost the same under different soil temperatures at pipeline buried depth.
Volume flow rate and its deviation of a longdistance pipeline at different hydrogen blending ratios when the pressures at the start and the end of the pipeline are fixed.
For further discussion, Figure 3a displays the profiles of gas volume flow rate against the hydrogen blending ratio under different soil temperatures at buried depth. It is seen that the three profiles corresponding to summer temperature, winter temperature and annual average temperature are approximately parallel to each other. That is, the influence of hydrogen blending on gas volume flow rate is basically the same in different seasons. It is also found from Figure 3a that at the same hydrogen blending ratio, the gas volume flow rate in summer is lower than that in winter. It indicates that with the increase of soil temperature at pipeline buried depth, the gas volume flow rate gradually decreases. It is worth noting that the gas transmission capacity of pipeline segment 2–3 is not improved although the gas volume flow rate increases with the rise of hydrogen blending ratio when the pressures at the start and the end of pipeline 2–3 are fixed. Figure 3b shows the corresponding energy flow rate of the pipeline under the conditions listed in Table 4. Compared with Figures 3a and 3b, it is found although the gas volume flow rate increases with the increasing hydrogen blending ratio, the energy flow rate of the pipeline declines almost linearly. The primary reason lies in the fact that under the same condition, the calorific value of hydrogen is about 1/3 of that of natural gas. Thus the mixed hydrogen decreases the calorific value of natural gas, and the increase of volume flow rate is not enough to offset the effect of hydrogen blending on the decrease of calorific value of natural gas.
Fig. 3 Variation of volume flow rate and energy flow rate against hydrogen blending ratio under fixed pressures at the start and the end of the pipeline. (a) Volume flow rate; (b) Energy flow rate. 
Figure 4 displays the temperature drop along the pipeline under different soil temperatures at buried depth and different hydrogen blending ratios when the pressures at the start and the end of the pipeline 2–3 are fixed. It is seen that the variation trend of temperature drop is consistent under different soil temperatures at buried depth. That is, the gas temperature gradually decreases with the increase of pipeline mileage when the gas source temperature remains unchanged. The lower the soil temperature at pipeline buried depth, the greater the total temperature drop along the pipeline, where the temperature drop is the largest in winter and reaches the minimum value in summer. Under the same soil temperature at pipeline buried depth and with the increase of hydrogen blending ratio, the variation trend of temperature drop becomes slower and the total temperature drop of the pipeline decreases gradually. For instance, when the hydrogen blending ratio is 30%, the temperature at the pipeline end is about 2 °C higher than that of the pipeline without hydrogen. The reason for this phenomenon can be attributed to the decrease of mass flow rate as hydrogen blending ratio increases, which leads to a larger coefficient a = KπD/mc_{p} in Sukhov formula for calculating the temperature distribution along the pipeline. In addition, the Joule–Thomson coefficient of natural gas decreases with the increasing hydrogen blending ratio. Compared with the natural gas without hydrogen, the mixed hydrogen has a similar effect with the thermal insulation on the pipeline. The above factors eventually result in that the temperature drop of HBNG is smaller than that of natural gas without hydrogen.
Fig. 4 Temperature drop at different soil temperatures and hydrogen blending ratios when pressures at the start and the end of the pipeline are fixed. (a) Soil temperature is 8.4 °C; (b) Soil temperature is 15 °C; (c) Soil temperature is 25.2 °C. 
When the energy flow rate and the end pressure of pipeline 2–3 are fixed, the influence of hydrogen blending on outlet pressure, pressure drop and temperature drop along the pipeline is systematically analyzed below. For the sake of simplicity, only the soil temperature of 15 °C at pipeline buried depth and the corresponding energy flow rate when transporting natural gas without hydrogen are taken for comparison. The temperature at compressor outlet is assumed to be 50 °C.
Table 5 and Figure 5a demonstrate the outlet pressure of the longdistance pipeline at different hydrogen blending ratios when the energy flow rate and terminal pressure are fixed. It can be seen when the hydrogen blending ratio rises from 0% to 30%, the outlet pressure of HBNG pipeline is 5.43% higher than that of traditional natural gas pipeline without hydrogen. The main reason is that for the longdistance pipeline, the adding of hydrogen into the natural gas breaks the original stable operation condition of the pipeline, and the outlet pressure of the compressor will increase to keep the energy flow rate and terminal pressure of the pipeline unchanged. Therefore, it is essential to increase the outlet pressure of the compressor when the energy flow rate and terminal pressure of the pipeline are fixed.
Fig. 5 Variation of outlet pressure and volume flow rate at different hydrogen blending ratios for fixed energy flow rate and end pressure of the pipeline. (a) Outlet pressure; (b) Volume flow rate. 
Outlet pressure and its deviation of a longdistance pipeline at different hydrogen blending ratios when the energy flow rate and the end pressure of the pipeline are fixed.
Figure 5a illustrates that the outlet pressure curves corresponding to different soil temperatures at pipeline buried depth have similar variation tendency, indicating the effect of soil temperature on outlet pressure is the same at same hydrogen blending ratio. In addition, since the volumetric calorific value of hydrogen is lower than that of natural gas, when the energy flow rate of pipeline remains unchanged, the volume flow rate of pipeline increases along with the rising hydrogen blending ratio, as shown in Figure 5b. The increased volume flow rate leads to a higher pressure drop along the pipeline. Therefore, when the energy flow rate and terminal pressure of the pipeline are fixed, the pressure drop increases with the raise of hydrogen blending ratio, and thus the outlet pressure of the compressor must be increased, as shown in Figure 6. Similarly, when the energy flow rate and terminal pressure of the pipeline are fixed, the temperature drop declines with the increase of hydrogen blending ratio, as displayed in Figure 7. The reasons for the effect of hydrogen blending ratio and soil temperature at buried depth on the pressure drop and temperature drop along the pipeline are basically the same, which is not analyzed here for brevity.
Fig. 6 Pressure drop at different soil temperatures and hydrogen blending ratios for fixed energy flow rate and end pressure of the pipeline. (a) Soil temperature is 8.4 °C; (b) Soil temperature is 15 °C; (c) Soil temperature is 25.2 °C. 
Fig. 7 Temperature drop at different soil temperatures and hydrogen blending ratios for fixed energy flow rate and end pressure of the pipeline. (a) Soil temperature is 8.4 °C; (b) Soil temperature is 15 °C; (c) Soil temperature is 25.2 °C. 
4.3 The influence of hydrogen blending on centrifugal compressor
Centrifugal compressor is the key equipment in the pipeline transportation system as it provides the required pressure head to overcome the gas flow resistance. When hydrogen is mixed into natural gas, the properties of HBNG are different from that of natural gas without hydrogen, which exerts substantial influences on the performance of centrifugal compressor. Thus it is of great significance to reveal the effect of hydrogen blending on compressor performance. In Table 6, the gas constant and adiabatic index of HBNG at different hydrogen blending ratios are presented. It can be obviously seen that the adiabatic index is weakly affected by the hydrogen blending, with an increase of 2.56% when the hydrogen blending ratio is 30%. Therefore, the effect of hydrogen blending on the adiabatic index is negligible when the hydrogen blending ratio is small. For the selection of similarity criterion and translation methods of compressor performance curves, literature study [31] found that under the same condition, the translation result of the first kind of similarity method is smaller than the experimental data, while that of the second type of similarity method is larger than the experimental results. However, the relative deviation of the former is much smaller. Hence in this work, the first kind of similarity criterion and translation method is applied to translate the performance curves of centrifugal compressor at different hydrogen blending ratios for convenient comparison.
Gas constant and adiabatic index of natural gas at different hydrogen blending ratios.
The centrifugal compressor in Figure 2 takes the PCL804 compressor as an example. Figures 8 and 9 show the translation curves of shaft power (W) and pressure ratio (p_{r}) at different hydrogen blending ratios and the rotational speed of 4880 rpm and 5490 rpm, respectively. It is seen at two different rotational speeds, the pressure ratiovolume flow rate curve and the shaft powervolume flow rate curve move down with the increasing hydrogen blending ratio. The reason for this phenomenon is that with the increase of hydrogen blending ratio, the density of HBNG decreases, resulting in the increase of gas constant R_{g} and the decrease of pressure ratio. Then the outlet pressure, total energy head and shaft power of the centrifugal compressor decrease as well.
Fig. 8 Performance curve of centrifugal compressor at different hydrogen blending ratios and the rotational speed of 4880 rpm. (a) Shaft power; (b) Pressure ratio. 
Fig. 9 Performance curve of centrifugal compressor at different hydrogen blending ratios and the rotational speed of 5490 rpm. (a) Shaft power; (b) Pressure ratio. 
Taking pipeline 1–3 in Figure 2 as an example, the influences of hydrogen blending on the joint operation of centrifugal compressor and pipeline are explored. Assuming that only the soil temperature of 15 °C at pipeline buried depth is considered, the volume flow rate is 3.84 × 10^{7} Nm^{3}/d when the pipeline transmission capacity is fixed, and the gas source temperature and pressure are assumed to be 34 °C and 6.74 MPa. Figure 10 displays the pressure drop and temperature drop when the pipeline is operated jointly with the centrifugal compressor. It is found both the outlet pressure and the outlet temperature gradually decrease with the increase of hydrogen blending ratio. Figure 11 shows the joint operation map of the pipeline and centrifugal compressor. It indicates with the increase of hydrogen blending ratio, the characteristic curves of both centrifugal compressor and pipeline move downward, and the operating point moves to the direction of higher volume flow rate and lower pressure. The primary reason for the performance degradation is that the HBNG has a smaller density compared with natural gas without hydrogen, which results in the rise of volume flow rate in the compressor inlet. Then the operating point moves towards the higher volume flow rate and lower pressure at the same time, which ultimately reduces the outlet pressure of centrifugal compressor.
Fig. 10 Pressure and temperature drop for joint operation of pipeline and compressor at 15 °C. (a) pressure drop; (b) temperature drop. 
Fig. 11 Performance curve for joint operation of pipeline and compressor. 
5 Effects of hydrogen blending on natural gas pipe network
5.1 Parameters of pipe network
In this part, the impacts of hydrogen blending on hydraulic and thermal characteristics of natural gas pipe network are studied by taking the ringshaped multiloop pipe network as an example. The topology of the pipe network is shown in Figure 12, with the pipe parameters and node flow information shown in Tables 7 and 8, where the negative value represents the outflow. The node 9 is the gas supply source, and the pressure of the gas supply source is 360 kPa. The composition of natural gas in the pipe network is the same as that in Table 3.
Fig. 12 Topology structure of pipe network. 
Parameters of pipe network.
Volume flow rate of each node in the pipe network.
5.2 Results and discussion
To investigate the influence of hydrogen blending on hydraulic and thermal characteristics of natural gas pipe network, the gas supply pressure and energy flow rate of the pipe network are kept constant. The variations of pressure and volume flow rate at each node of the pipe network under different hydrogen blending ratios are presented in Table 9 and Figure 13. As can be seen from Table 9 and Figure 13a, with the increase of hydrogen blending ratio, the volume flow rate at each node of the pipe network increases. When the hydrogen blending ratio reaches 30%, the volume flow rate at each node increases by about 25.48% compared with that without hydrogen. However, it should be noted although the volume flow rate at each node rises with the increasing hydrogen blending ratio when the gas supply pressure and energy flow rate of the pipe network are fixed, the pressure at each node decreases, as shown in Figure 13b. It is also found the influence of hydrogen blending on the pressure at different nodes varies a lot. That is, under the same hydrogen blending ratio, the pressure drop at the node far away from the gas source (node 9) is significantly larger than that at the node close to the gas source. Figure 13c demonstrates that the temperature drop of the pipe network gradually decreases with the increasing hydrogen blending ratio, showing the same variation trend as that of longdistance pipeline. However, owning to the short distance and low pressure of the pipe network, the temperature variation of pipe network is not obvious, which can be neglected in engineering practice. Thus the natural gas pipe network can be treated as an isothermal model in numerical simulations.
Fig. 13 Variations of node pressure, volume flow rate and gas temperature against hydrogen blending ratio at constant gas supply pressure and energy flow rate. (a) Node pressure; (b) Node volume flow rate; (c) Node temperature. 
Node pressure and volume flow rate under different hydrogen blending ratios at constant gas supply pressure and energy flow rate.
To further explore the impact of hydrogen blending, it is only assumed that the energy flow rate of the pipe network remains constant. Then the variations of pressure, flow rate and temperature at each node under different hydrogen blending ratios are obtained, as shown in Tables 10, 11 and Figure 14. It is seen with the increase of hydrogen blending ratio, most of the pressure at the distribution nodes of the pipe network shows a downward trend, while the pressure of the gas supply source and the volume flow rate of all nodes increase. It can be inferred with the increase of hydrogen blending ratio, the larger volume flow rate at each node and higher pressure of gas supply source are required to keep the energy flow rate of the pipe network constant. A similar temperature drop tendency can be observed from Figure 14c, that is, the temperature drop of the pipe network gradually decreases with the increasing hydrogen blending ratio. However, as illustrated above, the temperature change of the pipe network is negligible in engineering practice.
Fig. 14 Variations of node pressure, volume flow rate and gas temperature against hydrogen blending ratios at constant energy flow rate. (a) Node pressure; (b) Node volume flow rate; (c) Node temperature. 
Node pressure (kPa) under different hydrogen blending ratios at constant energy flow rate.
Node volume flow rate (10^{4} Nm^{3}/d) under different hydrogen blending ratios at constant energy flow rate.
6 Conclusion
In this paper, the influences of hydrogen blending on hydraulic and thermal characteristics of a longdistance natural gas pipeline and a ringshaped multiloop pipe network are numerically investigated. For the longdistance pipeline, the volume flow rate, energy flow rate, outlet pressure, pressure drop and temperature drop along the pipeline, and performance characteristics of the centrifugal compressor are analyzed at different hydrogen blending ratios. In study of a ringshaped multiloop pipe network, the variation trend of node pressure, node volume flow rate, node energy flow rate, and node temperature against the hydrogen blending ratio is comprehensively explored. The main findings of this study are summarized as follows:
For the longdistance natural gas pipeline with fixed pressures at the start and the end of the pipeline, the volume flow rate increases with the raise of hydrogen blending ratio, but the energy flow rate shows opposite trend. To keep the energy flow rate constant, the pressure at the compressor outlet should be increased. When the energy flow rate and pressure at the end of the pipeline are fixed, the outlet pressure gradually raises with the increasing hydrogen blending ratio, and the pressure drop along the pipeline increases as well.
The temperature drop along the pipeline is negatively related with the soil temperature at pipeline buried depth, that is, the lower the soil temperature, the greater the total temperature drop along the pipeline. Under the same soil temperature, the temperature drop along the pipeline and pipe network slows down with the increase of hydrogen blending ratio, thus the temperature drop becomes smaller for HBNG. Additionally, compared with the temperature drop of gas pipe network which is negligible, the temperature drop along a longdistance pipeline is more significant.
With the increase of hydrogen blending ratio, both the pressure ratio and the shaft power of the centrifugal compressor decrease, resulting in the operating point of the joint operation of compressor and pipeline moves to the direction of higher volume flow rate and lower pressure. Thus it is obvious to know with the increase of hydrogen blending ratio, the rotational speed of centrifugal compressor should be increased to satisfy the transportation requirement at a constant energy flow rate.
For the pipe network with fixed gas supply pressure and energy flow rate, the volume flow rate at the nodes increases with the raise of hydrogen blending ratio, while the pressure decreases in this process. At the same hydrogen blending ratio, the pressure drop at the nodes close to the gas supply source is smaller than that of the nodes far away from the gas supply source. When the energy flow rate is fixed, the pressure at the distribution nodes drops and the pressure of gas supply source raises with the increase of hydrogen blending ratio. Meanwhile, the volume flow rate of both distribution nodes and gas supply source increases.
Conflict of interest
The authors declare no competing financial interest.
Acknowledgments
This study is supported by the National Natural Science Foundation of China (Nos. 51904031, 51936001), the Natural Science Foundation of Beijing Municipality (No. 3204038), and the Award Cultivation Foundation from Beijing Institute of Petrochemical Technology (No. BIPTACF002).
References
 Solomon B.D., Krishna K. (2011) The coming sustainable energy transition: History, strategies, and outlook, Energy Policy 39, 11, 7422–7431. [Google Scholar]
 Ye B., Zhang K., Jiang J.J., Miao L.X., Li J. (2017) Towards a 90% renewable energy future: A case study of an island in the South China Sea, Energ. Convers. Manage. 142, 28–41. [Google Scholar]
 Liu B., Liu S.X., Guo S.S., Zhang S.X. (2019) Economic study of a largescale renewable hydrogen application utilizing surplus renewable energy and natural gas pipeline transportation in China, Int. J. Hydrogen Energ. 45, 3, 1385–1398. [Google Scholar]
 Öney F., Veziroğlu T.N., Dülger Z. (1994) Evaluation of pipeline transportation of hydrogen and natural gas mixtures, Int. J. Hydrogen Energ. 19, 10, 813–822. [Google Scholar]
 Asadnia M., Mehrpooya M. (2017) A novel hydrogen liquefaction process configuration with combined mixed refrigerant systems, Int. J. Hydrogen Energ. 42, 23, 15564–15585. [Google Scholar]
 Adam P., Engelshove S., Heunemann F., Thiemann T., Bussche C.V.D. (2020) Hydrogen infrastructurethe pillar of energy transition: The practical conversion of longdistance gas networks to hydrogen operation, Siemens Energy, Gascade Gastransport GmbH, Nowega GmbH. [Google Scholar]
 Haeseldonckx D., D’haeseleer W. (2007) The use of the naturalgas pipeline infrastructure for hydrogen transport in a changing market structure, Int. J. Hydrogen Energ. 32, 10–11, 1381–1386. [Google Scholar]
 Yang C., Ogden J. (2007) Determining the lowestcost hydrogen delivery mode, Int. J. Hydrogen Energ. 32, 2, 268–286. [Google Scholar]
 Messaoudani Z.L., Rigas F., Hamid M.D.B., Hassan C.R.C. (2016) Hazards, safety and knowledge gaps on hydrogen transmission via natural gas grid: A critical review, Int. J. Hydrogen Energ. 41, 39, 17511–17525. [Google Scholar]
 International Energy Agency (2019) The future of hydrogen: seizing today’s opportunities. [Google Scholar]
 Tiekstra G.C., Gasunie N.V.N., Folkert P.K., Gasunie N.V.N. (2008) The NATURALHY project: first step in assessing the potential of the existing natural gas network for hydrogen delivery, IGRC, Paris. [Google Scholar]
 Kippers M.J., De Laat J.C., Hermkens R.J.M., Overdiep J.J., Van Der Molen A., Van Erp W.C., Van Der Meer A. (2001) Pilot project on hydrogen injection in natural gas on island of Ameland in the Netherlands, Proc. Int. Gas Research Conf. 2, 1163–1177. [Google Scholar]
 Klopffer M.H., Berne P., Espuche É. (2015) Development of innovating materials for distributing mixtures of hydrogen and natural gas. Study of the barrier properties and durability of polymer pipes, Oil Gas Sci. Technol.Rev. IFP Energies nouvelles 70, 2, 305–315. [Google Scholar]
 Tabkhi F., AzzaroPantel C., Pibouleau L., Domenech S. (2008) A mathematical framework for modelling and evaluating natural gas pipeline networks under hydrogen injection, Int. J. Hydrogen Energ. 33, 11, 6222–6231. [Google Scholar]
 Elaoud S., HadjTaïeb E. (2008) Transient flow in pipelines of highpressure hydrogennatural gas mixtures, Int. J. Hydrogen Energ. 33, 18, 4824–4832. [Google Scholar]
 Elaoud S., Hafsi Z., HadjTaïeb L. (2017) Numerical modelling of hydrogennatural gas mixtures flows in looped networks, J. Petrol. Sci. Eng. 159, 532–541. [Google Scholar]
 Guandalini G., Colbertaldo P., Campanari S. (2017) Dynamic modeling of natural gas quality within transport pipelines in presence of hydrogen injections, Appl. Energ. 185, 2, 1712–1723. [Google Scholar]
 Wang W., Wang Q.Y., Deng H.Q., Cheng G.X., Li Y. (2020) Feasibility analysis on the transportation of hydrogennatural gas mixtures in natural gas pipelines, Nat. Gas Ind. 403, 130–136 (in Chinese). [Google Scholar]
 Hafsi Z., Elaoud S., Mishra M. (2019) A computational modelling of natural gas flow in looped network: Effect of upstream hydrogen injection on the structural integrity of gas pipelines, J. Nat. Gas Sci. Eng. 64, 107–117. [Google Scholar]
 Hafsi Z., Elaoud S., Akrout M., HadjTaïeb E. (2017) Numerical approach for steady state analysis of hydrogennatural gas mixtures flows in looped network, Arab. J. Sci. Eng. 42, 5, 1941–1950. [Google Scholar]
 Wu C. (2018) Feasibility study on blending hydrogen into natural gas distribution network, Chongqing University, Chongqiang, China (in Chinese). [Google Scholar]
 Uilhoorn F. (2009) Dynamic behaviour of nonisothermal compressible natural gases mixed with hydrogen in pipelines, Int. J. Hydrogen Energ. 34, 16, 6722–6729. [Google Scholar]
 Li Y.X., Yao G.Z. (2009) Design and management of natural gas pipeline, 2nd edn., Academic Press, Dongying, China (in Chinese). [Google Scholar]
 Soave G. (1972) Equilibrium constants from a modified RedlichKwong equation of state, Chem. Eng. Sci. 27, 6, 1197–1203. [Google Scholar]
 Peng D.Y., Robinson D.B. (1976) A new twoconstant equation of state, Ind. Eng. Chem. Fundam. 15, 1, 59–64. [Google Scholar]
 Soave G.S. (1999) An effective modification of the Benedict–Webb–Rubin equation of state, Fluid Phase Equilibr. 164, 2, 157–172. [Google Scholar]
 Varzandeh F., Stenby E.H., Yan W. (2017) Comparison of GERG2008 and simpler EoS models in calculation of phase equilibrium and physical properties of natural gas related systems, Fluid Phase Equilibr. 434, 21–43. [Google Scholar]
 Wang P., Yu B., Deng Y., Zhao Y. (2015) Comparison study on the accuracy and efficiency of the four forms of hydraulic equation of a natural gas pipeline based on linearized solution, J. Nat. Gas Sci. Eng. 22, 235–244. [Google Scholar]
 Wang P., Ao S.M., Yu B., Han D.X., Xiang Y. (2019) An efficiently decoupled implicit method for complex natural gas pipeline network simulation, Energies 12, 8, 1516. [Google Scholar]
 Xiang Y., Wang P., Yu B., Sun D.L. (2020) GPUaccelerated hydraulic simulations of largescale natural gas pipeline networks based on a twolevel parallel process, Oil Gas Sci. Technol.Rev. IFP Energies nouvelles 75, 86. [Google Scholar]
 Wang L.Y., Wang P., Cao Z.Z., Yu B., Li W. (2017) Similarity conversion of centrifugal natural gas compressors based on predictorcorrector, Procedia Comput. Sci. 108, 1973–1981. [Google Scholar]
All Tables
Volume flow rate and its deviation of a longdistance pipeline at different hydrogen blending ratios when the pressures at the start and the end of the pipeline are fixed.
Outlet pressure and its deviation of a longdistance pipeline at different hydrogen blending ratios when the energy flow rate and the end pressure of the pipeline are fixed.
Gas constant and adiabatic index of natural gas at different hydrogen blending ratios.
Node pressure and volume flow rate under different hydrogen blending ratios at constant gas supply pressure and energy flow rate.
Node pressure (kPa) under different hydrogen blending ratios at constant energy flow rate.
Node volume flow rate (10^{4} Nm^{3}/d) under different hydrogen blending ratios at constant energy flow rate.
All Figures
Fig. 1 Discretization in space and time. 

In the text 
Fig. 2 Sketch map of the studied longdistance pipeline. 

In the text 
Fig. 3 Variation of volume flow rate and energy flow rate against hydrogen blending ratio under fixed pressures at the start and the end of the pipeline. (a) Volume flow rate; (b) Energy flow rate. 

In the text 
Fig. 4 Temperature drop at different soil temperatures and hydrogen blending ratios when pressures at the start and the end of the pipeline are fixed. (a) Soil temperature is 8.4 °C; (b) Soil temperature is 15 °C; (c) Soil temperature is 25.2 °C. 

In the text 
Fig. 5 Variation of outlet pressure and volume flow rate at different hydrogen blending ratios for fixed energy flow rate and end pressure of the pipeline. (a) Outlet pressure; (b) Volume flow rate. 

In the text 
Fig. 6 Pressure drop at different soil temperatures and hydrogen blending ratios for fixed energy flow rate and end pressure of the pipeline. (a) Soil temperature is 8.4 °C; (b) Soil temperature is 15 °C; (c) Soil temperature is 25.2 °C. 

In the text 
Fig. 7 Temperature drop at different soil temperatures and hydrogen blending ratios for fixed energy flow rate and end pressure of the pipeline. (a) Soil temperature is 8.4 °C; (b) Soil temperature is 15 °C; (c) Soil temperature is 25.2 °C. 

In the text 
Fig. 8 Performance curve of centrifugal compressor at different hydrogen blending ratios and the rotational speed of 4880 rpm. (a) Shaft power; (b) Pressure ratio. 

In the text 
Fig. 9 Performance curve of centrifugal compressor at different hydrogen blending ratios and the rotational speed of 5490 rpm. (a) Shaft power; (b) Pressure ratio. 

In the text 
Fig. 10 Pressure and temperature drop for joint operation of pipeline and compressor at 15 °C. (a) pressure drop; (b) temperature drop. 

In the text 
Fig. 11 Performance curve for joint operation of pipeline and compressor. 

In the text 
Fig. 12 Topology structure of pipe network. 

In the text 
Fig. 13 Variations of node pressure, volume flow rate and gas temperature against hydrogen blending ratio at constant gas supply pressure and energy flow rate. (a) Node pressure; (b) Node volume flow rate; (c) Node temperature. 

In the text 
Fig. 14 Variations of node pressure, volume flow rate and gas temperature against hydrogen blending ratios at constant energy flow rate. (a) Node pressure; (b) Node volume flow rate; (c) Node temperature. 

In the text 