Prediction of Physical Properties and Phase Characteristics of Ethane and Ethane Mixture in the Ethane Pipeline

: In order to realize the safe transportation of liqueﬁed ethane pipeline in the Oilﬁeld of China, it is necessary to fully study the process of pipeline replacement, operation and shutdown. The accurate calculation of physical property parameters and critical parameters is the basis of studying the gas-liquid two-phase ﬂow and heat and mass transfer process of liqueﬁed ethane in the pipeline. In this paper, different equations of states (EOSs) were used to predict the physical properties (such as density, dew point and dynamic viscosity) of ethane or ethane mixture, and the predicted results were compared with the corresponding experimental data from the literature. The prediction performance of different EOSs were evaluated by using two evaluation indicators, including average absolute deviation (AAD) and average relative deviation (ARD). The results showed that the PR-Peneloux EOS has the best performance for predicting the density of CH 4 -C 2 H 6 -N 2 mixture with an ARD value of 4.46%; for predicting the dew point, the BWRS EOS exhibits the superior performance with an ARD value of 0.58%; and for predicting dynamic viscosity, the SuperTRAPP formula has the smallest calculation error, with an ARD value of 1.33%. Considering the comparison results of the calculation accuracy of density, dew point and dynamic viscosity of ethane or ethane mixture by using different EOSs, PR-Peneloux EOS was recommended to calculate the phase characteristics in the process of ethane pipeline replacement operation. The phase characteristics of ethane for pipeline transport in the oilﬁeld of China were obtained. The critical temperature is 32.79°C and the critical pressure is 4.97 MPa.


Introduction
Ethane (C 2 H 6 ) is used as a raw material to produce ethylene, which is the most important basic raw material for the petrochemical industry [1][2][3][4]. It can be used to produce almost all chemical products with a critical pressure and critical temperature of 4.88 MPa and 32.25°C, respectively, which can be transported by a liquid phase or gas phase in pipelines. Compared with the gas phase ethane transportation method, the liquid phase ethane transportation method has the advantages of large transportation capacity and high economy [5][6][7]. Liquid phase ethane pipelines have been built in some countries [8][9][10]. In China, Petrochina Company proposes to transport ethane by means of liquefaction, and plans to build a long-distance ethane pipeline. Since the liquefied ethane is a saturated liquid in the pipeline, it is easy to vaporize due to external temperature disturbance or internal pressure fluctuation, forming a two-phase flow in the pipeline, which is obviously different from the conventional crude oil and natural gas pipelines [11,12]. In addition, this can cause some transportation problems such as pipeline vibration, ultra-pressure and equipment damage. At present, the production process and methods for natural gas, refined oil, liquefied petroleum gas and other media cannot solve the technical problems of the liquid phase ethane transportation pipeline.
In order to realize the safe transportation of liquefied ethane pipelines, it is necessary to fully study the process of pipeline replacement, operation and shutdown. Pressure and temperature are the basis for the safe transportation of liquefied ethane pipelines. When the transportation pressure of liquefied ethane is lower than the critical pressure of commercial ethane or the temperature is higher than its critical temperature, vaporization of liquefied ethane in the pipeline will occur, which will endanger the safe transportation of the pipeline. The accurate calculation of ethane physical property parameters is the key factor to accurately predict the flow parameters of liquefied ethane pipelines, such as temperature, pressure and flow rate.
Commonly, the equation of state (EOS) is used to obtain the physical properties of pure gas or a mixture of gases. For example, PR EOS and GERG-2008 EOS were employed to calculate the bubble point of pure ethane, and the PR EOS showed a superior performance with an average deviation of 0.41% [13]. Later, GERG-2008 EOS was used to calculate the density of the mixture of propane and hydrogen with different mixed ratios, and the results also showed that the GERG-2008 EOS can effectively predict the density of the mixture gas [14]. In addition, Vitali et al. [15] employed GERG-2008 EOS to predict the vapor-liquid equilibrium and density of the CO 2 -rich mixture. A quantitative estimation was conducted, and the results showed that the GERG-2008 EOS was more accurate in the description of both vapor-liquid equilibrium and density when compared with cubic EOS. Except for the prediction of the bubble point, density and vapor-liquid equilibrium of the above gases, Seo et al. [16] used extended Redlich-Kwong-Peng-Robinson (eRK-PR) EOS to predict the thermodynamic properties of JP-10, which is an important fuel in the petroleum and aviation industries. The results showed that the eRK-PR EOS predicts the thermodynamic properties of JP-10 more accurately than the other EOSs. The above results show that EOS can be used to calculate the key physical parameters with good accuracy. However, the prediction performance of EOSs for the physical properties and phase characteristics of ethane mixture are still unclear. The accurate calculation of physical property parameters depends on the applicability of the EOS to the calculation of ethane parameters. In addition, the EOS is also the basis of the study of ethane phase characteristics, which determines the accuracy of the calculation of critical pressure and critical temperature of ethane, and further affects the safe replacement of the liquefied ethane pipeline into production, operation and shutdown.
This paper conducted relevant research to address the issues of inaccurate prediction accuracy of physical property parameters in CH 4 -C 2 H 6 -N 2 mixture and unclear phase state change law during the replacement process of liquid-phase ethane pipeline. For the first time, seven EOSs (PR, PR-Peneloux and SRK) were evaluated based on experimental data for predicting density and dew point of CH 4 -C 2 H 6 -N 2 mixture. In addition, LBC, Pedersen and Super TRAPP equations were also evaluated for predicting viscosity. Furthermore, utilizing the preferred equation of state, we were able to predict the physical property parameters of the mixture under varying concentration ratios, pressures and temperatures during the operation of the liquid ethane pipeline. Finally, the key physical property parameters of liquid ethane in relevant oil fields and the phase transition rules during pipeline operation are presented. The research findings provide a theoretical foundation for safe displacement operations of liquid-phase ethane pipelines and subsequent technical standards formulation for easily vaporized substance pipelines.

Cubic-Type PVT EOS
For the cubic-type PVT EOS, if the incompressibility of the liquid is considered, as V approaches some parameter b, P goes to infinity, and then the general cubic form of P is: The parameters such as b, θ, η, δ and ε in Equation (1) can be constants including zero, or they can vary with temperature or composition. The values of common cubic-type PVT equation of state parameters are listed in Tables 1 and 2.
where V is molar volume, m 3 /mol; R is the gas constant, 8.314 J/(mol·k); P is pressure, Pa; T is pressure, K; Z is the compression factor; a and b are characteristic parameters; y i is the molar fraction of component i; K ij is the bivariate interaction coefficient; T c is the critical temperature of natural gas, K; and P a is the volume translation coefficient.

Non-Cubic PVT EOS
(1) BWRS EOS Based on a large number of experimental data, Starling et al. [21] modified the BWR equation and proposed a BWRS equation with higher accuracy for predicting the thermodynamic parameters of various hydrocarbon mixtures, which is expressed as follows where P is pressure, Pa; T is pressure, K; ρ is density kg/m 3 ; and R is the gas constant, 8.314 J/(mol·k). A 0 , B 0 , C 0 , D 0 , E 0 , a, b, c, d, α, γ are state parameters.
(2) LKP EOS On the basis of the multi-parameter BWRS EOS, LKP EOS was proposed [22], which can simultaneously calculate the fugacity coefficients of vapor-liquid two-phase where z is compression factor; T is eccentricity factor; p r is pressure ratio; V r is molar volume ratio; T r is temperature ratio; the superscript "0" is a simple fluid parameter; and the superscript "R" is the reference fluid parameter.
(3) GERG-2008 EOS The GERG-2008 EOS is also used as a standard to calculate the thermodynamic properties of natural gas. GERG-2008 EOS is generally expressed by Helmholtz free energy, which are given by [23] α ρ, T, X = α 0 ρ, T, X + α r ρ, T, X where X is the molar fraction of each component in the mixture, %; δ is density variable; and τ is temperature variable. The above EOS adopts the Van der Waals mixing rule when analyzing and calculating the physical properties of the gas mixture. Considering that non-polar molecules (such as, water) are involved in the process of replacing N 2 with ethane in the liquid phase, the whole system is in a non-polar system. Therefore, the Van der Waals mixing rule can be used for correlation analysis and the calculation of the physical properties of the mixed gas. For hydrocarbon mixture systems, Van der Waals mixing rules are generally used to calculate the gravitational parameters and volume terms, which are given by [24]: The RK equation can obtain a certain degree of accuracy when calculating the heat capacity of pure components and mixtures, but the accuracy of the RK equation is often not ideal when considering the calculation of multi-component gas-liquid equilibrium. In addition, the American Gas Association (AGA) proposed the AgA8-92DC equation to calculate the compressibility coefficient of natural gas (namely, AGA8 equation), but this equation is mainly applicable to the calculation condition of methane content being not less than 70% and the temperature being higher than −10°C. Therefore, RK equation and AGA8 equation are not considered in the calculation of this paper. The CPA equation of state is also often used to complete the physical property calculation, but is no different from the SRK equation of state in the calculation of non-polar molecules such as CH 4 , C 2 H 6 and N 2 , so it is not considered separately.
In summary, it is preliminarily concluded that SRK, PR, BWRS, LKP and GREG2008 EOS can be used to calculate the phase characteristics of C 2 H 6 and C 2 H 6 -N 2 mixed component gases. Next, the applicability and accuracy of each equation of state for mixed gases will be further explored, thus laying a foundation for the next phase characteristics research.

EOS Evaluation for Predicting Physical Properties
To evaluate the calculation accuracy of each EOS on the physical properties of liquid ethane, it is necessary to calculate the data of different physical properties, such as density and viscosity, dew point, etc. and compared with the corresponding experimental data, so as to select the most accurate equation of state and lay the foundation for the subsequent research on the phase characteristic of ethane. Based on the above considerations, in order to quantitatively analyze the prediction accuracy of each equation of state, some commonly used model evaluation parameters, such as average absolute deviation (AAD) and average relative deviation (ARD), are defined as follows where N is the total number of the data samples; y cal i and y exp i represent the calculated value and experimental data, respectively. The smaller the calculated values of AAD and ARD, the better the prediction performance of the EOSs.

Evaluation of EOSs for Predicting Density
In this study, the empirical formula, PR EOS, PR-Peneloux EOS, SRK EOS, SRK-Peneloux EOS, BWRS EOS, LKP EOS, GREG-2008 EOS are used to predict the density of the CH 4 -C 2 H 6 -N 2 mixture. The calculation accuracy of the density of the CH 4 -C 2 H 6 -N 2 mixture is validated and analyzed. Experimental data of different mixed concentrations of CH 4 -C 2 H 6 -N 2 mixture under different pressure and temperature conditions are derived from Funke's work, and it is employed to validate the calculation results of the above empirical formulas. The density data of CH 4 -C 2 H 6 -N 2 mixture under different temperatures, pressures and molar concentrations in the experiments by Funke et al. [25] is used in this work, which is listed in Table 3. It can be clearly seen from Table 1 that the concentration of C 2 H 6 plays a leading role in the density change of the CH 4 -C 2 H 6 -N 2 mixture. The higher the concentration of C 2 H 6 , the higher the density of the CH 4 -C 2 H 6 -N 2 mixture. In addition, it can be found that when the pressure and temperature are constant, the density of the CH 4 -C 2 H 6 -N 2 mixture increases by 1.48% for every 1% increase in the C 2 H 6 concentration. The density increase of the CH 4 -C 2 H 6 -N 2 mixture is mainly because the density of C 2 H 6 is higher than that of CH 4 and N 2 (under standard conditions, the density of C 2 H 6 is 1.356 kg/m 3 , the density of CH 4 is 0.717 kg/m 3 and the density of N 2 is 1.251 kg/m 3 ). Therefore, when the pressure and temperature are constant or the range of change is low, the density of CH 4 -C 2 H 6 -N 2 mixture is shown to increase with the increase of the concentration of C 2 H 6 .
The density of CH 4 -C 2 H 6 -N 2 mixture of twenty samples in Table 3 are calculated using PR EOS, PR-Peneloux EOS, SRK EOS, SRK-Peneloux EOS, BWRS EOS, LKP EOS and GREG-2008 EOS at different experimental temperatures, pressures and molar concentrations. The empirical formulas were used to calculate the mixture density, and the results were compared with experimental data from Funke et al. [25], as shown in Figure 1. It can be found that the calculated value of density of CH 4 -C 2 H 6 -N 2 mixture using different empirical formula possess different accuracies. Furthermore, it also can be seen that in most working conditions (temperature in the range of −179~−181 • C, pressure in the range of 0.1~1.1 MPa), the prediction of density of the CH 4 -C 2 H 6 -N 2 mixture by the above empirical formulas have a certain degree of accuracy. On the contrary, the PR-Peneloux EOS has the best performance for predicting the density of CH4-C2H6-N2 mixture. This is because compared with the PR EOS, the volume translation coefficient c is introduced into the PR-Peneloux EOS, which can further modify the characteristic parameters in the equation and improve the accuracy of density calculation of the PR EOS in the low temperature liquid state. Figure 2b shows the ARD of the density of CH4-C2H6-N2 mixture using the empirical formula, PR EOS, PR-Peneloux EOS, SRK EOS, SRK-Peneloux EOS, BWRS EOS, LKP EOS and GREG-2008 EOS compared with the experimental data. It can be found that the value of ARD for all EOSs is less than 20%. This illustrates that the PR EOS, PR-Peneloux EOS, SRK EOS, SRK-Peneloux EOS, BWRS EOS, LKP EOS and GREG-2008 EOS are capable of accurately predicting the density of CH4-C2H6-N2 mixture in the conditions of −179~−181 °C and 0.1~1.1 MPa. PR-Peneloux EOS can give the best prediction of the density of CH4-C2H6-N2 mixture with an ARD of 4.46%, which is smaller than 8.89% of SRK EOS, 10.07% of BWRS EOS, 11.08% of SRK-Peneloux EOS, 12.66% of GREG-2008 EOS, 13.06% of LKP EOS and 17.5% of PR EOS. In order to further evaluate the calculation accuracy of the above empirical formulas, the AAD and ARD are calculated and analyzed, as shown in Figure 1. It can be clearly seen from Figure 2a  Obviously, the PR EOS has the worst performance for predicting the density of CH 4 -C 2 H 6 -N 2 mixture; this is because the gas to be calculated is a mixture of ethane, methane and N 2 , and the coefficient a, reflecting the attraction between different molecules, and the coefficient b, reflecting the repulsive force in the PR EOS, are out of alignment when calculating the density parameters of the mixed components, resulting in the density values calculated by the PR EOS all being higher than the corresponding experimental values. Based on the above analysis, we can draw the conclusion that the above EOSs can efficiency predict the density of CH4-C2H6-N2 mixture and the prediction accuracy is in the order of PR-Peneloux EOS > SRK EOS > BWRS EOS > SRK-Peneloux EOS > GREG-2008 EOS > LKP EOS > PR EOS.

Evaluation of EOSs for Predicting Dew Point
The phase state envelope diagram of C2H6-N2 mixed components is composed of a dew point curve and bubble point curve, which is an important basis for analyzing and judging the phase state of C2H6, N2 and C2H6-N2 mixed components. The experimental data of dew point of ethane-N2 mixed components by Bier et al. [26] and Syed et al. [27] are collected, which can be divided into 5 groups with 64 sets of data of different temperatures and pressure working conditions, as shown in Table 4. In the same way, these experimental data are used to validate the calculation results of the EOSs. On the contrary, the PR-Peneloux EOS has the best performance for predicting the density of CH 4 -C 2 H 6 -N 2 mixture. This is because compared with the PR EOS, the volume translation coefficient c is introduced into the PR-Peneloux EOS, which can further modify the characteristic parameters in the equation and improve the accuracy of density calculation of the PR EOS in the low temperature liquid state. Figure 2b shows the ARD of the density of CH 4 -C 2 H 6 -N 2 mixture using the empirical formula, PR EOS, PR-Peneloux EOS, SRK EOS, SRK-Peneloux EOS, BWRS EOS, LKP EOS and GREG-2008 EOS compared with the experimental data. It can be found that the value of ARD for all EOSs is less than 20%. This illustrates that the PR EOS, PR-Peneloux EOS, SRK EOS, SRK-Peneloux EOS, BWRS EOS, LKP EOS and GREG-2008 EOS are capable of accurately predicting the density of CH 4 -C 2 H 6 -N 2 mixture in the conditions of −179~−181 • C and 0.1~1.1 MPa. PR-Peneloux EOS can give the best prediction of the density of CH 4 -C 2 H 6 -N 2 mixture with an ARD of 4.46%, which is smaller than 8.89% of SRK EOS, 10.07% of BWRS EOS, 11.08% of SRK-Peneloux EOS, 12.66% of GREG-2008 EOS, 13.06% of LKP EOS and 17.5% of PR EOS.
Based on the above analysis, we can draw the conclusion that the above EOSs can efficiency predict the density of CH 4 -C 2 H 6 -N 2 mixture and the prediction accuracy is in the order of PR-Peneloux EOS > SRK EOS > BWRS EOS > SRK-Peneloux EOS > GREG-2008 EOS > LKP EOS > PR EOS.

Evaluation of EOSs for Predicting Dew Point
The phase state envelope diagram of C 2 H 6 -N 2 mixed components is composed of a dew point curve and bubble point curve, which is an important basis for analyzing and judging the phase state of C 2 H 6 , N 2 and C 2 H 6 -N 2 mixed components. The experimental data of dew point of ethane-N 2 mixed components by Bier et al. [26] and Syed et al. [27] are collected, which can be divided into 5 groups with 64 sets of data of different temperatures and pressure working conditions, as shown in Table 4. In the same way, these experimental data are used to validate the calculation results of the EOSs.  Figure 3 shows the comparison of dew point of C 2 H 6 -N 2 mixture between experimental data and calculation results by different EOSs under different working conditions (pressure: 0.52~13.01 MPa, temperature 192.52~301.7 K). Obviously, all the predicted values of all EOSs almost coincide with the corresponding experimental data when the concentration of ethane is higher than 84%, as shown in Figure 3a,b. However, with the continuous decrease of ethane concentration (ranging from 49.82% to 84.99%), some deviations between the predicted values of the EOSs and the experimental data began to appear, and the smaller the concentration of ethane, the greater the deviation, as shown in Figure 3c,d. However, when the concentration of ethane is lower 49.82%, the deviation between the predicted value and the experimental value of the equation of state begins to decrease. In addition, it also can be found from Figure 3 that BWRS EOS can accurately predict the dew point value of C 2 H 6 -N 2 mixture for all ethane and N 2 mixed ratios.   Figure 3 shows the comparison of dew point of C2H6-N2 mixture between experimental data and calculation results by different EOSs under different working conditions (pressure: 0.52~13.01 MPa, temperature 192.52~301.7 K). Obviously, all the predicted values of all EOSs almost coincide with the corresponding experimental data when the concentration of ethane is higher than 84%, as shown in Figure 3a,b. However, with the continuous decrease of ethane concentration (ranging from 49.82% to 84.99%), some deviations between the predicted values of the EOSs and the experimental data began to appear, and the smaller the concentration of ethane, the greater the deviation, as shown in Figure 3c,d. However, when the concentration of ethane is lower 49.82%, the deviation between the predicted value and the experimental value of the equation of state begins to decrease. In addition, it also can be found from Figure 3 that BWRS EOS can accurately predict the dew point value of C2H6-N2 mixture for all ethane and N2 mixed ratios.  To better evaluate the performance of EOSs for predicting the dew point of C2H6-N2 mixture, the AAD and ARD of different EOSs are calculated, as shown in Figure 4. From Figure 4, it is seen that the AAD and ARD values of SRK EOS are the smallest among all EOSs under the working conditions of 1# (C =95.02, C = 4.98) and 2# (C =84.99, C = 15.01). This illustrates that the SRK EOS has the best performance for predicting the dew point of C2H6-N2 mixture with a higher concentration of C2H6 (>84.99%). However, when the concentration of C2H6 is less than 84.99%, the AAD and ARD values of BWRS EOS are the smallest among all EOSs. That is to say, the BWRS EOS has the best calculation accuracy for dew point of C2H6-N2 mixture under the above conditions. To better evaluate the performance of EOSs for predicting the dew point of C 2 H 6 -N 2 mixture, the AAD and ARD of different EOSs are calculated, as shown in Figure 4. From Figure 4, it is seen that the AAD and ARD values of SRK EOS are the smallest among all EOSs under the working conditions of 1# (C C 2 H 6 =95.02, C N 2 = 4.98) and 2# (C C 2 H 6 = 84.99, C N 2 = 15.01). This illustrates that the SRK EOS has the best performance for predicting the dew point of C 2 H 6 -N 2 mixture with a higher concentration of C 2 H 6 (>84.99%). However, when the concentration of C 2 H 6 is less than 84.99%, the AAD and ARD values of BWRS EOS are the smallest among all EOSs. That is to say, the BWRS EOS has the best calculation accuracy for dew point of C 2 H 6 -N 2 mixture under the above conditions.  To better evaluate the performance of EOSs for predicting the dew poi mixture, the AAD and ARD of different EOSs are calculated, as shown in F Figure 4, it is seen that the AAD and ARD values of SRK EOS are the small EOSs under the working conditions of 1# (C =95.02, C = 4.98) and 2# C = 15.01). This illustrates that the SRK EOS has the best performance the dew point of C2H6-N2 mixture with a higher concentration of C2H6 (>8 ever, when the concentration of C2H6 is less than 84.99%, the AAD and A BWRS EOS are the smallest among all EOSs. That is to say, the BWRS EO calculation accuracy for dew point of C2H6-N2 mixture under the above con  In order to select the optimal EOS applicable to all working conditions AAD and ARD under above five working conditions are averaged again, as ure 5. Obviously, the average AAD and ARD values of BWRS EOS are the values of 1.40 and 0.58%, respectively. Namely, the BWRS EOS has the bes for predicting the dew point of C2H6-N2 mixture regardless of the mixed rat N2. On the contrary, the LKP EOS exhibits the worst performance for predi point of C2H6-N2 mixture, and the corresponding values of AAD and ARD 1.50%, respectively. The reason for this is that in the process of calculating of C2H6-N2 mixture, the mixing rule used in the LKP EOS to calculate the m tual critical parameters has a deviation, which lead to a large deviation in th  In order to select the optimal EOS applicable to all working conditions, the values of AAD and ARD under above five working conditions are averaged again, as shown in Figure 5. Obviously, the average AAD and ARD values of BWRS EOS are the smallest, with values of 1.40 and 0.58%, respectively. Namely, the BWRS EOS has the best performance for predicting the dew point of C 2 H 6 -N 2 mixture regardless of the mixed ratio of C 2 H 6 and N 2 . On the contrary, the LKP EOS exhibits the worst performance for predicting the dew point of C 2 H 6 -N 2 mixture, and the corresponding values of AAD and ARD are 3.75 and 1.50%, respectively. The reason for this is that in the process of calculating the dew point of C 2 H 6 -N 2 mixture, the mixing rule used in the LKP EOS to calculate the multi-comp virtual critical parameters has a deviation, which lead to a large deviation in the LKP EOS. In order to select the optimal EOS applicable to all working conditions, the values of AAD and ARD under above five working conditions are averaged again, as shown in Figure 5. Obviously, the average AAD and ARD values of BWRS EOS are the smallest, with values of 1.40 and 0.58%, respectively. Namely, the BWRS EOS has the best performance for predicting the dew point of C2H6-N2 mixture regardless of the mixed ratio of C2H6 and N2. On the contrary, the LKP EOS exhibits the worst performance for predicting the dew point of C2H6-N2 mixture, and the corresponding values of AAD and ARD are 3.75 and 1.50%, respectively. The reason for this is that in the process of calculating the dew point of C2H6-N2 mixture, the mixing rule used in the LKP EOS to calculate the multi-comp virtual critical parameters has a deviation, which lead to a large deviation in the LKP EOS.

Evaluation of EOSs for Predicting Dynamic Viscosity
The dynamic viscosity calculation is very different from the density calculation when calculating the physical properties of the components. The density value can be directly obtained by relying on the EOSs. However, due to the complexity of the gaseous components in the oilfield and the dynamic viscosity also being a parameter related to density, pressure and temperature, if the EOS is directly used for dynamic viscosity calculation, there is a certain error in the calculation of density which will lead to the deviation of

Evaluation of EOSs for Predicting Dynamic Viscosity
The dynamic viscosity calculation is very different from the density calculation when calculating the physical properties of the components. The density value can be directly obtained by relying on the EOSs. However, due to the complexity of the gaseous components in the oilfield and the dynamic viscosity also being a parameter related to density, pressure and temperature, if the EOS is directly used for dynamic viscosity calculation, there is a certain error in the calculation of density which will lead to the deviation of dynamic viscosity calculation results. Therefore, for predicting dynamic viscosity, a special empirical formula or semi-empirical formula is often used in filed engineering to solve the dynamic viscosity value, such as mixing rule, LBC (Lohrenz-Bray-Clark) [28], Pedersen [29] and SuperTRAPP equations based on TRAPP equations, etc. These empirical equations have built-in reference gases such as methane, propane, etc. In summary, the EOS and the viscosity equation are combined for calculating the parameters used in the viscosity equation. The experimental data of dynamic viscosities from Carr et al.'s [30] work are employed to validate the prediction performance of the EOSs, which is listed in Table 5.  Figure 6 presents the comparison of dynamic viscosity between the predicted results and experimental data at a different molar concentration ratio of C 2 H 6 and N 2 . Clearly, all empirical formulas predict dynamic viscosity with good accuracy under the conditions of atmospheric pressure and temperature of 300 K. The AAD and ARD results for the prediction of dynamic viscosity using different empirical formulas is shown in Figure 7.

Phase Characteristic of C2H6 or N2 Component
According to the definition of the U.S. Energy Information Administration (EIA), commercial ethane is ethane with a content of 95% or more used in the market. The ethane It can be found from Figure 7 that among the above empirical formulas, the Super-TRAPP formula has the smallest calculation error, with AAD and ARD values of 0.16 and 1.33%, respectively. The reason for this is that the SuperTRAPP formula is an extended prediction model of propane as a reference fluid. The viscosity values of propane and ethane are similar under the same working condition. Therefore, in the concentration system dominated by ethane, the calculated viscosity is closer to the experimental value. On the contrary, the Pedersen formula has the largest calculation error, with AAD and ARD values of 0.57 and 4.36%, respectively. The SuperTRAPP formula is recommended for the calculation of dynamic viscosity. The calculated viscosity of the Pedersen equation is higher than the experimental value when the ethane concentration is less than 10%, but with the increase of ethane concentration, the calculated viscosity is lower than the experimental value. This may be because the Pedersen formula is a relative EOS, which is based on the viscosity change rule of the reference substance liquid methane to predict the viscosity of other substances.

Phase Characteristic of C 2 H 6 or N 2 Component
According to the definition of the U.S. Energy Information Administration (EIA), commercial ethane is ethane with a content of 95% or more used in the market. The ethane product produced by the Oilfield of China is composed of 97% ethane, 1.5% methane and 1.5% propane, which belongs to commercial ethane. The phase state of commercial ethane in the Oilfield of China and pure N 2 are calculated using PR-Peneloux-EOS, as shown in Figures 8 and 9. temperature and a higher critical pressure, so it is difficult to form gas-liquid two-phase flow under normal pipeline transportation conditions. In the process of replacing N2 with liquid-phase ethane, the content of N2 drops from 100% to 0%, while the content of ethane rises from 0% to 100%. In this replacement process, the proportion of each component of the mixed gas is complicated, so it is still necessary to further study the phase characteristics of the mixed components of ethane and N2 based on the study of the phase states of ethane and N2.

Phase Characteristic of C2H6-N2 Mixed Component
In the process of replacing nitrogen with liquid-phase ethane, due to the large difference in the physical properties of the two media inside the pipeline, high purity liquidphase ethane will vaporize after contacting nitrogen to form gas-phase ethane, methane and other components. These gases will be mixed with nitrogen to form a new mixed gas, and with the continuous progress of the ethane replacement process, the proportion of nitrogen components in the mixture gas will continue to decrease. In addition, components such as ethane and methane in high-purity ethane will continue to increase until the proportion of nitrogen in the mixture gas inside the pipeline decreases to 0% and the proportion of high-purity ethane becomes 100%. Then, it can be considered that the nitrogen replacement is completed. At this time, the pipeline is full of high-purity ethane, as shown in Figure 10. It can be seen from Figure 8 that the critical temperature and pressure of commercial ethane in the Oilfield of China are 32.79 • C and 4.97 MPa, respectively. Meanwhile, the critical temperature and pressure of N 2 are −145.39°C and 3.30 MPa, respectively, as can be found in Figure 9. In pipeline replacement production state, liquefied ethane is more likely to undergo a gas-liquid phase change than natural gas. The gas-phase ethane formed by vaporization may also be liquefied again, resulting in the occurrence of a gas-liquid two-phase flow, which greatly reduces the safety and pipeline transportation efficiency in the process of transportation. Figure 9 reveals that pure N 2 has a lower critical temperature and a higher critical pressure, so it is difficult to form gas-liquid two-phase flow under normal pipeline transportation conditions. In the process of replacing N 2 with liquid-phase ethane, the content of N 2 drops from 100% to 0%, while the content of ethane rises from 0% to 100%. In this replacement process, the proportion of each component of the mixed gas is complicated, so it is still necessary to further study the phase characteristics of the mixed components of ethane and N 2 based on the study of the phase states of ethane and N 2 .

Phase Characteristic of C 2 H 6 -N 2 Mixed Component
In the process of replacing nitrogen with liquid-phase ethane, due to the large difference in the physical properties of the two media inside the pipeline, high purity liquid-phase ethane will vaporize after contacting nitrogen to form gas-phase ethane, methane and other components. These gases will be mixed with nitrogen to form a new mixed gas, and with the continuous progress of the ethane replacement process, the proportion of nitrogen components in the mixture gas will continue to decrease. In addition, components such as ethane and methane in high-purity ethane will continue to increase until the proportion of nitrogen in the mixture gas inside the pipeline decreases to 0% and the proportion of high-purity ethane becomes 100%. Then, it can be considered that the nitrogen replacement is completed. At this time, the pipeline is full of high-purity ethane, as shown in Figure 10.

Phase Characteristic of C2H6-N2 Mixed Component
In the process of replacing nitrogen with liquid-phase ethane, due to the large difference in the physical properties of the two media inside the pipeline, high purity liquidphase ethane will vaporize after contacting nitrogen to form gas-phase ethane, methane and other components. These gases will be mixed with nitrogen to form a new mixed gas, and with the continuous progress of the ethane replacement process, the proportion of nitrogen components in the mixture gas will continue to decrease. In addition, components such as ethane and methane in high-purity ethane will continue to increase until the proportion of nitrogen in the mixture gas inside the pipeline decreases to 0% and the proportion of high-purity ethane becomes 100%. Then, it can be considered that the nitrogen replacement is completed. At this time, the pipeline is full of high-purity ethane, as shown in Figure 10. During the whole replacement process, with the continuous progress of replacement construction, the proportion of gas components inside the pipeline is constantly changing. Therefore, it is necessary to carry out phase-state characteristics research on the high purity ethane and nitrogen gas with different mixing ratios and to lay the foundation for the subsequent research on the mixing state and flow process of high-purity ethane and nitrogen gas in the replacement process. During the whole replacement process, with the continuous progress of replacement construction, the proportion of gas components inside the pipeline is constantly changing. Therefore, it is necessary to carry out phase-state characteristics research on the high purity ethane and nitrogen gas with different mixing ratios and to lay the foundation for the subsequent research on the mixing state and flow process of high-purity ethane and nitrogen gas in the replacement process. Figure 11 shows the phase diagram of C 2 H 6 -N 2 with different mixing ratios calculated using PR-Peneloux EOS. As can be seen, the physical properties of ethane are changed after being mixed with nitrogen. Its virtual critical pressure keeps increasing, while the virtual critical temperature keeps decreasing, which results in the enlarging of the gas-liquid two-phase region. The normal transportation pressure in the liquid-phase ethane pipeline is 5~8 MPa, while in the ethane replacement operation, the pipeline pressure ranges from 0 MPa to 15 MPa and the temperature is between −50~40°C. Therefore, based on this consideration, the phase diagram of C 2 H 6 -N 2 in actual operating conditions is given, as shown in Figure 11b. According to Figure 1b, the critical parameters of C 2 H 6 -N 2 mixed gas with different mixing ratios within the actual working conditions can be clarified and the phase state of the pipeline transportation medium can be determined according to the pressure and temperature at different stages of displacement production. In addition, although nitrogen can be dissolved in ethane, the dissolution condition is very harsh. When the pressure is normal pressure and the temperature drops to 90 K, the molar proportion of nitrogen after dissolution is only 7%, and when the temperature rises to 110 K, the molar proportion is only 3%. Therefore, the component of nitrogen dissolved in liquidphase ethane can be ignored when calculating the physical parameters of the pipeline replacement medium and flow process under the pressure and temperature conditions of the displacement production process.
harsh. When the pressure is normal pressure and the temperature drops to 90 K, the molar proportion of nitrogen after dissolution is only 7%, and when the temperature rises to 110 K, the molar proportion is only 3%. Therefore, the component of nitrogen dissolved in liquid-phase ethane can be ignored when calculating the physical parameters of the pipeline replacement medium and flow process under the pressure and temperature conditions of the displacement production process.

Conclusions
In this paper, different EOSs were used to predict the physical properties (such as density, dew point and dynamic viscosity) of ethane or ethane mixture, and the prediction performance was evaluated by two evaluation indicators, including average absolute deviation (AAD) and average relative deviation (ARD). Then, the phase characteristics of ethane or ethane mixture in the process of ethane pipeline replacement operation were obtained based on the optimal EOS. The main conclusions are as follows: (2) All EOSs can accurately predict the dew point of the C2H6-N2 mixture when the concentration of ethane is greater than 84%. The BWRS EOS has the best performance with an average ARD value of 0.58%, regardless of the mixed ratio of C2H6 and N2.

Conclusions
In this paper, different EOSs were used to predict the physical properties (such as density, dew point and dynamic viscosity) of ethane or ethane mixture, and the prediction performance was evaluated by two evaluation indicators, including average absolute deviation (AAD) and average relative deviation (ARD). Then, the phase characteristics of ethane or ethane mixture in the process of ethane pipeline replacement operation were obtained based on the optimal EOS. The main conclusions are as follows:  Data Availability Statement: Data will be made available on request.

Conflicts of Interest:
The authors declare no conflict of interest.