Effects of hydrogen blending on hydraulic and thermal characteristics of natural gas pipeline and pipe network

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 Hydrogen-Blended 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.


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 low-carbon-intensity fuels (such as natural gas, hydrogen, etc.) are of great significance to the energy conservation, emission reduction and energy structure optimization [1][2][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 long-distance 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][5][6][7]. It is also noted that the transportation of hydrogen through natural gas pipelines is a low-cost option and is beneficial for emission reduction [8][9][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][12][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 long-distance 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 Hadj-Taïeb [15] and Elaoud et al. [16] assessed the pressure variation of a small-scale 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 long-distance 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 gas-hydrogen 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 Hydrogen-Blended 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 long-distance natural gas pipeline with the centrifugal compressor and a ringed multi-loop 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 long-distance pipeline and the performance of a centrifugal compressor; Section 5 discusses the impact of hydrogen blending on the hydraulic and thermal characteristics of a ring-shaped multi-loop pipe network; Concluding remarks of this work are summarized in Section 6.

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], where A is the cross-section area of the pipeline, m 2 ; p is the pressure, Pa; q is the gas density, kg/m 3 ; w is the gas velocity, m/s; k is the coefficient of hydraulic friction; g is the gravitational acceleration, g = 9.81 m/s 2 ; h 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 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 SRK-EOS [24], PR-EOS [25], BWRS-EOS [26], and GERG2008-EOS [27], etc. In this study, the widely used BWRS-EOS in natural gas engineering is applied, where p is the gas pressure, kPa; T is the gas temperature, K; q 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, a and c represent the model parameters in BWRS-EOS. 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 long-distance 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, 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, where n denotes the rotational speed, rpm; , c e2 and c e3 are the parameters to be fitted in the pressure ratio curve; , c g2 and c g3 are the parameters to be determined in the polytropic efficiency curve.

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, 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.
To discretize the problem in space and time, the spatial step is set as Dx and the time step is set as Dt. 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.
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 second-order central difference scheme. While for the thermal governing equation, the convection term is discretized by the first-order upwind scheme. Readers of interest can refer to [28] for more detail. The discrete hydraulic governing equation can be arranged as follows, 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], where 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 Graphical-Processing-Unit (GPU)-accelerated algorithm [30].
4 Effects of hydrogen blending on longdistance natural gas pipeline

Pipeline parameters
Taking a long-distance 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.
To quantitatively analyze the influence of hydrogen blending on hydraulic and thermal states of the long-distance 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, where U 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 hydraulic-thermal 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

Results and discussion
Firstly, the influence of hydrogen blending on gas transmission capacity of a long-distance 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 long-distance pipeline is almost the same under different soil temperatures at pipeline buried depth.
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. 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 = KpD/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 Table 4. Volume flow rate and its deviation of a long-distance pipeline at different hydrogen blending ratios when the pressures at the start and the end of the pipeline are fixed. 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.
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 long-distance 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. 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.

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.
The centrifugal compressor in Figure 2 takes the PCL-804 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 ratio-volume flow rate curve and the shaft power-volume 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.
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.
5 Effects of hydrogen blending on natural gas pipe network

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 ring-shaped multi-loop 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.

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 long-distance 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.
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.

Conclusion
In this paper, the influences of hydrogen blending on hydraulic and thermal characteristics of a long-distance natural gas pipeline and a ring-shaped multi-loop pipe network are numerically investigated. For the long-distance 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 ring-shaped multi-loop 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: 1. For the long-distance 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. 2. 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 long-distance pipeline is more significant. 3. 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. 4. 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.