A Wet Gas Metering System Based on the Extended-Throat Venturi Tube

Although the use of a classical Venturi tube for wet gas metering has been extensively studied in the literature, the use of an extended-throat Venturi (ETV) tube has rarely been reported since its first proposal by J. R. Fincke in 1999. The structure of an ETV is very simple, but due to the complexity of multiphase flow, its theoretical model has not been fully established yet. Therefore, in this paper theoretical models have been developed for the convergent and throat sections of an ETV, and the gradients of front and rear differential pressures are derived analytically. Several flowrate algorithms have been proposed and compared with the existing ones. Among them, the iteration algorithm is found to be the best. A reasonable explanation is provided for its performance. The relationship between the differential pressure gradient and the flowrate relative error is also studied, such that the relative error distributions varying with ETV measured flowrates can be derived. The gas flowrate error of ETV increases with the liquid content whilst the liquid flowrate error of ETV decreases with the liquid content, and the relative errors of liquid flowrate are generally 2 to 3 times larger than that of the gas flowrate. Finally, the ETV tends to be more accurate than the classical Venturi tube. The ETV can be designed more compact under the same signal intensity due to its significantly higher velocity in the throat section.


Introduction
Wet gas typically refers to a type of two-phase flow with the gas volume fraction (GVF) larger than 95% or the Lockhard Martinelli (L-M) number less than 0.3. Compared with the oil-gas-water multiphase flow, the flow pattern and governing rule of wet gas are relatively stable and simple, and therefore, there are many companies dedicated to wet gas metering and the most commonly employed device is the classical Venturi tube. The structure and algorithm of a classical Venturi tube for single-phase metering are well known and included in the ISO 5167-4 [1], and its algorithms for wet gas metering have been extensively studied in the literature and finally included in the ISO/TR 11583 [2]. In contrast, the study of an extended-throat Venturi (ETV) tube for wet gas metering is rarely reported in the literature since its first proposal by J. R. Fincke [3,4]. As shown in Figure 1, the basic structure of an ETV tube is very simple. It includes a convergent section for front differential pressure dp 1 measurement and a throat section for rear differential pressure dp 2 measurement, as well as a pressure sensor at the inlet and a temperature flowrate. These two examples essentially exploited the difference between a differential pressure (DP) flowmeter and a velocity (volumetric) flowmeter to yield the gas and liquid flowrates, which is different from the conventional ETV that uses the difference between two DP flowmeters for flowrate measurement. Therefore, these types of ETV are beyond the scope of this paper and only conventional ETV will be discussed in the following parts.
In this paper, the theoretical models of the convergent and throat section of an ETV tube are established, and the working principle of the wet gas metering system is explained. Thereafter, several flowrate algorithms based on direct fitting method and iteration method are proposed and compared with the existing ones. It is found that the iteration method enjoys higher prediction accuracy than the linear fitting method, and a reasonable explanation is provided. Finally, the gradients of the front differential pressure (∂dp1/∂Qg, ∂dp1/∂Ql) and rear differential pressure (∂dp2/∂Qg, ∂dp2/∂Ql) are derived, and their relationships with the relative errors of gas and liquid flowrate are studied. As a result, the error distribution maps of an ETV are obtained from the test results, and these maps are compared with the ones of a classical Venturi tube found in the literature. The results show that the ETV tends to have higher measurement accuracy than the classical Venturi tube. In addition, the ETV can be designed to be more compact due to its significantly higher flow velocity in the throat section.  In recent years, other measuring techniques, such as electrical capacitance [11,12], photon (gamma and X-ray) attenuation [13] and ultrasound [14,15], have been increasingly used in ETV for determining the void fraction and flow regimes of multiphase flow. For example, a wet-gas flowmeter based on an ETV tube and a sonar flowmeter has been proposed by Weatherford [16,17]. It uses the sonar flowmeter to determine the convection velocity and thus the volumetric flowrate. Another example is a multiphase flowmeter that combines an ETV with two microwave resonators [18]. The time-based responses of these two sensors are correlated to extract the average flow velocity and the volumetric flowrate. These two examples essentially exploited the difference between a differential pressure (DP) flowmeter and a velocity (volumetric) flowmeter to yield the gas and liquid flowrates, which is different from the conventional ETV that uses the difference between two DP flowmeters for flowrate measurement. Therefore, these types of ETV are beyond the scope of this paper and only conventional ETV will be discussed in the following parts.
In this paper, the theoretical models of the convergent and throat section of an ETV tube are established, and the working principle of the wet gas metering system is explained. Thereafter, several flowrate algorithms based on direct fitting method and iteration method are proposed and compared with the existing ones. It is found that the iteration method enjoys higher prediction accuracy than the linear fitting method, and a reasonable explanation is provided. Finally, the gradients of the front differential pressure (∂dp 1 /∂Q g , ∂dp 1 /∂Q l ) and rear differential pressure (∂dp 2 /∂Q g , ∂dp 2 /∂Q l ) are derived, and their relationships with the relative errors of gas and liquid flowrate are studied. As a result, the error distribution maps of an ETV are obtained from the test results, and these maps are compared with the ones of a classical Venturi tube found in the literature. The results show that the ETV tends to have higher measurement accuracy than the classical Venturi tube. In addition, the ETV can be designed to be more compact due to its significantly higher flow velocity in the throat section.

Convergent Section
When pure gas flows through an ETV tube, and the measured front differential pressure is dp 1 , then its the volume flowrate can be directly calculated by: where β = d/D, A = πd 2 /4, C d is the discharge coefficient, ε is the expansion factor, and the values of C d and ε can be determined from ISO 5167-4 [1]. ρ g denotes the gas density, which can be calculated from the measured pressure p and temperature T according to ISO 12213-2 [19]. However, when wet gas flows through an ETV, the flowrate obtained from Equation (1) is referred to as the "indicated" gas flowrate Q tp , and the "actual" gas flowrate Q g should be calculated by Q g = Q tp /φ g , where φ g is referred to as the "gas over-reading". φ g is always larger than one because the front differential pressure dp 1 of wet gas is contributed by the gas and liquid together. φ g is usually calculated from the Lockhard-Martinelli number X, and the equation employed is referred to as the over-reading equation. For example, Murdock [20], Bizon [21], Lin et al. [22] proposed a linear model like: where X = Q l Q g ρ l ρ g , whereas Chisholm [23], de Leeuw [24], Reader-Harris [25] proposed a general model like: where: where S 1 = u g /u l is referred to as the slip ratio. Equations (3) and (4) are referred to as the general model because its basic form can be theoretically derived. For example, if the pressure drop of convergent section is assumed to be dominated by the accelerate pressure drop [26], then the differential form of the pressure drop can be written as: By assuming the cross-sectional void fraction α as a constant and integrating Equation (5) from the Venturi inlet to the throat, we have: By using similar method for gas and liquid single-phase, we have: From Equations (6) and (1), it can be derived that the gas over-reading is in the form of Equations (3) and (4). Meanwhile, it can be derived from Equation (7) that S 1 = ρ l ρ g for this scenario, and this model is referred to as the stratified model as no friction force is assumed between the gas and liquid. On the other hand, if the force between gas and liquid is very strong so that no velocity slip exists, which means S 1 = 1, then this model is referred to as the homogeneous model. The stratified and homogeneous model can be viewed as two extreme scenarios of S 1 , therefore, S 1 usually varies between 1 and ρ l ρ g in actual working conditions, and its exact value is also dependent on the gas Froude number Fr g , the types of fluids and the structure parameter of a Venturi tube (e.g., diameter ratio β).
If the Equation (1) and the definitions of X and φ g are substituted into Equation (3) and simplified, then the following expression can be obtained: where Fr g is the gas Froude number, Fr l is the liquid Froude number, gD(ρ l −ρ g ) and A D = πD 2 /4. It is worth mentioning that X = Fr l /Fr g . Therefore, the gradient of the front differential pressure dp 1 can be written as: From Equation (9) it can be noted that the gradient of dp 1 is a function of the slip ratio S 1 , the density ratio ρ l ρ g and the Lockhard-Martinelli number X. Its direction is not affected by the gas Froude number Fr g and the parameter K 1 .

Throat Section
When pure gas flows through a throat tube with a diameter of d and length of L, its pressure drop can be calculated by the Fanning's law of friction: where τ wg is the shear stress of gas on wall, τ wg = λρ g u 2 go 8 ; P go is the perimeter of gas at pipe wall, P go = πd for pure gas; A is the cross-sectional area of the throat, A = πd 2 4 . λ is the pipe friction factor, if the pipe wall is rough and flow is turbulent, then λ is a constant which does not change with the Reynold number Re.
After some modifications of Equation (10), the indicated gas flowrate of the throat section can be obtained: When wet gas flows through the same throat tube, the pressure drop of liquid and gas can be derived from the momentum conservation equation [23]: for liquid, and: − A g dp tp dz = τ wg P g + S F (13) for gas, where S F is the shear force per unit length at the interface between the phases. If the pipe friction factor λ and the hydraulic diameter χ = 4P i /A i are assumed the same for the gas and liquid, then divide Equation (12) by Equation (13), we have: After substituting Equation (14) into Equation (13) and cancelling S F , we have: From Equations (15) and (10), it can be derived that: where Z = 1 S 2 ρ l ρ g , as shown in Equation (14). The expression of Equation (16) for the throat tube is identical with Equation (3) for the convergent one, but it is generally believed that the flow within it is more or less homogenous, so that S 2 ≈ 1 [27].
Similar to the convergent section, the throat section also has its own linear model as Equation (2) and general model as Equation (3), and the gradient of rear differential pressure dp 2 can also be written as similar form as Equation (9): where gD(ρ l −ρ g ) .

Working Principle
By comparing Equations (9) and (17), it can be noted that for the ETV tube, the L-M number X and density ratio ρ l ρ g of the convergent section and the throat section are always the same, therefore, the only possible reason for the intersection angle between ∇(dp 1 ) and ∇(dp 2 ) is the difference in slip ratio S. Some scholars believe that the flow in the throat section is very close to the homogeneous flow (e.g., S 2 = 1 or n 2 = 0.5) [27], whilst the flow in the convergent section may have a higher slip ratio S 1 due to the existence of the acceleration pressure drop (e.g., S 1 > 1, or n 1 < 0.5), and the ETV uses this difference to calculate the gas and liquid flowrates.
In the ideal case, if there is only acceleration pressure drop in the convergent section whilst there is only frictional pressure drop in the throat section, then the contours of the non-dimensional front differential pressure dp 1 and rear differential pressure dp 2 are shown in Figure 2a. From Figure 2a, it can be noted that there is always an intersection angle between ∇(dp 1 ) and ∇(dp 2 ) and the intersection angle varies with X, as shown in Equations (9) and (17). In the real case, the index n of Equation (3) is obtained from data fitting [28]. The variations of index n of several common over-reading models with gas Froude number Fr g are shown in Figure 2b [28], from which it can be noted that the index n is also affected by the fluid type and diameter ratio β. Therefore, it will be difficult to continue using the analytical method for the derivation of index n. is also affected by the fluid type and diameter ratio . Therefore, it will be difficult to continue using the analytical method for the derivation of index .

Flowrate Algorithms
The experiments of ETV were carried out at the wet gas test facility of Chengdu Verification Branch of National Oil and Gas Large Flowrate Measurement Station of China (CVB) which is designed for gas-liquid two phase flow studies consisting of water and natural gas. The schematic diagram of the CVB wet gas test facility is shown in Figure 3a, from which it can be noted that this facility is mainly composed of the gas and liquid reference meters, the liquid injection and regulation system and the gas-liquid separation system. The gas used for the wet gas experiments comes from a high pressure pipeline, and is metered by a high accuracy ultrasonic flowmeter before mixing. Then the natural gas flows into the test section where the liquid is injected and mixed. The flowrate of the liquid is regulated by the liquid injection and regulation system. A 0.5 inch and a 1 inch Coriolis mass flow meters work as liquid reference meters on the liquid injection line. The natural gas drives the liquid through the test section and then flow into the gas-liquid separation system where the gas and liquid are separated. The separated liquid is returned to the water tank and recirculated by the liquid pump, whereas the separated natural gas is discharged into a low-pressure pipeline. The test pressure range of the CVB is from 1.5 to 4.0 MPa, the gas flow test range is from 8 m 3 /h to 650 m 3 /h, and the liquid flow test range is from 0 to 8 m 3 /h. The CVB is a Class A metering station with a measurement uncertainty ( = 2) as low as 0.05%~0.07%, and more relevant information can be found in [29,30].
In order to simulate the flow conditions of the Yakela gas fields in northwestern China, the whole system is pressurized to 2.0 MPa, and the gas and liquid flowrate ranges are 32~156 m 3 /h and 0~4.1 m 3 /h respectively. There are 31 test points in total, where 5 of them represent the pure gas working conditions, whilst the rest 26 points represent the wet gas ones and their distribution are shown in Figure 3b

Flowrate Algorithms
The experiments of ETV were carried out at the wet gas test facility of Chengdu Verification Branch of National Oil and Gas Large Flowrate Measurement Station of China (CVB) which is designed for gas-liquid two phase flow studies consisting of water and natural gas. The schematic diagram of the CVB wet gas test facility is shown in Figure 3a, from which it can be noted that this facility is mainly composed of the gas and liquid reference meters, the liquid injection and regulation system and the gas-liquid separation system. The gas used for the wet gas experiments comes from a high pressure pipeline, and is metered by a high accuracy ultrasonic flowmeter before mixing. Then the natural gas flows into the test section where the liquid is injected and mixed. The flowrate of the liquid is regulated by the liquid injection and regulation system. A 0.5 inch and a 1 inch Coriolis mass flow meters work as liquid reference meters on the liquid injection line. The natural gas drives the liquid through the test section and then flow into the gas-liquid separation system where the gas and liquid are separated. The separated liquid is returned to the water tank and recirculated by the liquid pump, whereas the separated natural gas is discharged into a low-pressure pipeline. The test pressure range of the CVB is from 1.5 to 4.0 MPa, the gas flow test range is from 8 m 3 /h to 650 m 3 /h, and the liquid flow test range is from 0 to 8 m 3 /h. The CVB is a Class A metering station with a measurement uncertainty (k = 2) as low as 0.05%~0.07%, and more relevant information can be found in [29,30].
In order to simulate the flow conditions of the Yakela gas fields in northwestern China, the whole system is pressurized to 2.0 MPa, and the gas and liquid flowrate ranges are 32~156 m 3 /h and 0~4.1 m 3 /h respectively. There are 31 test points in total, where 5 of them represent the pure gas working conditions, whilst the rest 26 points represent the wet gas ones and their distribution are shown in Figure 3b  All of the measured ETV have the same diameter ratio of = ⁄ = 0.5, whereas their inlet diameters can be 50 or 80 mm to satisfy the requirements of different flowrate ranges. As the effect of on the flowrate algorithm is negligible, only the test results of = 50 mm are shown below to avoid redundancy. The detection frequency of the ETV tube is = 10 Hz.  For each measurement, there is an associated error. Therefore, it is sometimes more convenient to use mean absolute error (MAE) or mean absolute percentage error (MAPE) to compare the average values of errors of different data sets. The definitions of MAE and MAPE are as follows: where is the measured value, is the reference value and is the number of measurements of a dataset. The ranges of MAE of MAPE are both from 0 to infinity and larger MAE and MAPE mean lower accuracy. However, MAE is usually used for the errors of ratios (e.g., and ⁄ ) whereas MAPE is mainly used for the errors of flowrates (e.g., and ).

Linear Model
Simultaneous solution of the linear over-reading Equation (2) of the convergent and throat sections leads to the linear models as follows: Matrix inversion can be applied to Equation (20) directly and the result is: All of the measured ETV have the same diameter ratio of β = d/D = 0.5, whereas their inlet diameters D can be 50 or 80 mm to satisfy the requirements of different flowrate ranges. As the effect of D on the flowrate algorithm is negligible, only the test results of D = 50 mm are shown below to avoid redundancy. The detection frequency of the ETV tube is f = 10 Hz.
For each measurement, there is an associated error. Therefore, it is sometimes more convenient to use mean absolute error (MAE) or mean absolute percentage error (MAPE) to compare the average values of errors of different data sets. The definitions of MAE and MAPE are as follows: whereŷ i is the measured value, y i is the reference value and n is the number of measurements of a dataset. The ranges of MAE of MAPE are both from 0 to infinity and larger MAE and MAPE mean lower accuracy. However, MAE is usually used for the errors of ratios (e.g., φ g and Y/Y max ) whereas MAPE is mainly used for the errors of flowrates (e.g., Q g and Q l ).

Linear Model
Simultaneous solution of the linear over-reading Equation (2) of the convergent and throat sections leads to the linear models as follows: Matrix inversion can be applied to Equation (20) directly and the result is: where However, it is worth mentioning that the matrix of Equation (20) is relatively ill conditioned, so large errors may be generated during the matrix inversion. Therefore, variable substitution method may be used instead and the gas and liquid flowrates are directly fitted with the indicated gas flowrate Q gt1 and Q gt2 , so as to prevent this matrix inversion process.
The linear model is easy to solve but has a limited application range, so it is usually used together with the classification method. However, as the linear model enjoys the merits of simplicity and stability, it can be used for: (1) providing initial values for iteration; (2) providing reference for classification (such as estimating the value of Fr g ); (3) fix the "bad values" or "blind spots" of some complicated algorithms.

General Model
Reader-Harris and Graham improves the conventional Chisholm model [25]. They believe in the gas-liquid two-phase flow, the discharge coefficient of the Equation (1) changes, so they added a correction factor before the Chisholm model to achieve better data fitting results [25]: By solving the above two equations simultaneously, we can get: where Similar to the linear model, variable substitutions can also be implemented to prevent the errors introduced during the solving process.

Over-Reading Equations
At present, there are a large number of over-reading equations related to wet gas in the literature. Solartron ISA summarizes this [31,32], and most of the common overreading equations can be written as the form of Equations (3) and (4). This equation is also referred to as the general model as its basic form can be derived theoretically, but the index n is usually determined experimentally and it has become the key part for a good overreading model. Currently, the most common over-reading models for wet gas measurement includes:

1.
Homogeneous Model The homogeneous model assumes the gas and liquid are uniformly mixed without any slip, so the wet gas is treated as an pseudo-single phase, and the constant C 1 in Equations (3) and (4) becomes: So n = 0.5 and it is notable that the index n of homogeneous model is a constant. Therefore, if the L-M number is known, this model can obtainφ g directly and then estimate the gas/liquid flowrate. In this paper, the homogeneous model will be used for providing the initial values for other more complicated models. de Leeuw Model de Leeuw tested a Venturi tube with an inner diameter of four inches and β = 0.4 under different pressures with nitrogen and diesel oil as the flow medium. It is found that the index n varies between 0.41 and 0.606, and its specific form is: It is notable that the de Leeuw model needs to know Fr g first and calculates n and C 1 latter, soφ g should be determined through an iteration method.

Reader-Harris Model
Reader-Harris conducted a series of experiments for a Venturi tube with 0.4 ≤ β ≤ 0.75, and obtained the index n as follows: where H depends on the liquid types, which is 1.0 for hydrocarbon liquid, 1.35 for water at ambient and 0.79 for liquid water in a wet-steam flow.
In addition, the discharge coefficient is also changed as follows: Because C d ≈ 1 is assumed for other over-reading models, so the φ RH g calculated by Equations (3) and (4)   From Figure 4, it can be noted that the data fitting results of Reader-Harris model (ISO/TR 11583 model) are slightly better than the de Leeuw model, whilst the data fitting results of the de Leeuw model are slightly better than the homogenous model. However, considering the complexity of the model, the number of parameters involved and the risk of overfitting, de Leeuw model is employed in this paper to calculate the gas over-reading φ g and its initial value is provided by the homogeneous model.

Determination of X
The liquid content of wet gas is usually described by the L-M number X = Q l Q g ρ l ρ g . According to the research by de Leeuw [24] and Reader-Harris [25,33], X is usually obtained from the data fitting with the pressure loss ratio d dp , where d represents the pressure loss across a classical Venturi tube, dp represents the acceleration pressure drop through the convergent section of Venturi, as shown in Figure 5. For the ETV, the above two parameters correspond to the rear differential pressure dp 2 and front differential pressure dp 1 respectively, as shown in Figure 1.
ing of an ETV.
From Figure 4, it can be noted that the data fitting results of Reader-Harris model (ISO/TR 11583 model) are slightly better than the de Leeuw model, whilst the data fitting results of the de Leeuw model are slightly better than the homogenous model. However, considering the complexity of the model, the number of parameters involved and the risk of overfitting, de Leeuw model is employed in this paper to calculate the gas over-reading and its initial value is provided by the homogeneous model.

Determination of
The liquid content of wet gas is usually described by the L-M number = .
According to the research by de Leeuw [24] and Reader-Harris [25,33], is usually obtained from the data fitting with the pressure loss ratio , where d represents the pressure loss across a classical Venturi tube, d represents the acceleration pressure drop through the convergent section of Venturi, as shown in Figure 5. For the ETV, the above two parameters correspond to the rear differential pressure d and front differential pressure d respectively, as shown in Figure 1. dp D dϖ d Figure 5. Schematic diagram of a classical Venturi.
The pressure loss ratio should be normalized first: The pressure loss ratio should be normalized first: and then data fitted in the following form: For the classical Venturi tube, data fitting are carried out by Reader-Harris and the final form are embodied by the ISO/TR 11583 [2]: As there is obvious structure difference between the classical Venturi tube and the ETV, so if Equation (32) is applied to the ETV directly then serious deviations will occur, as shown in Figure 6. So, it is necessary to refit Equations (30) and (31) for the ETV and the results are as follows: because the test data covered the two extreme cases of dry gas and = 0.4.
The fitting results of the improved Reader-Harris correlation for the normalize ferential pressure ratio of an ETV are also shown in Figure 6, where the MAE (defin Equation (18)) in the title refers to the MAE of the improved Reader-Harris correlati Figure 6. The estimates of the original and improved Reader-Harris correlation for the norma differential pressure ratio of an ETV.

Calculation Procedures
The iteration method combines the over-reading equation of the de Leeuw m (Equations (3), (4) and (27)) with the pressure loss ratio correlation of the Reader-H model (Equations (31) and (33)), and calculates the gas and liquid flowrate through tion. As both the over-reading equation and the pressure loss ratio correlation requi Y/Y max Figure 6. The estimates of the original and improved Reader-Harris correlation for the normalized differential pressure ratio of an ETV.
The fitting results of the improved Reader-Harris correlation for the normalized differential pressure ratio of an ETV are also shown in Figure 6, where the MAE (defined in Equation (18)) in the title refers to the MAE of the improved Reader-Harris correlation.

Calculation Procedures
The iteration method combines the over-reading equation of the de Leeuw model (Equations (3), (4) and (27)) with the pressure loss ratio correlation of the Reader-Harris model (Equations (31) and (33)), and calculates the gas and liquid flowrate through iteration. As both the over-reading equation and the pressure loss ratio correlation require the gas Froude number Fr g to be known, this algorithm assumesX = 0 and uses the homogeneous model (Equations (3), (4) and (26)) for the initial values. The specific procedures of the iteration method are as follows:

1.
First, the initial values of the L-M numberX is assumed.

2.
Then the gas over-readingφ g is estimated according to the homogenous model (Equation (26)).

3.
Then the gas/liquid flowrates (Q g andQ l ) are calculated from the definitions of gas over-readingφ g and L-M numberX respectively. 4.
ThenFr g andFr l are calculated by its definition.

5.
Then the gas over-readingφ g is calculated by the de Leeuw model (Equations (3), (4) and (27)). Afterφ g is calculated,φ g andX are substituted into the Step 3 to continue the following process untilφ g converges. 6.
Then the L-M numberX is calculated by the pressure loss ratio correlation of the Reader-Harris model (Equations (30) and (33)). AfterX is calculated, it is then substituted into the Step 2 to continue the following process untilX converges.
Finally, the lastQ g andQ l are recorded and used as the estimates of the gas and liquid flowrates, respectively.

Comparison of the Direct Fitting and Iteration Method
The gas and liquid flowrate estimates of the linear, general and iteration models are shown in Figures 7 and 8. The black circles denote the linear model results, the blue squares denote the general model results, and the red asterisks denote the iteration model results. The central black line denotes the ideal case when the estimated value is equal to the reference value so that the error is always zero. The upper and lower red lines specify a 10% relative error range so that points within this range have relative errors less than 10%. Similarly, the upper and lower black lines specify a 20% relative error range. The error is represented by the the vertical distance between the test point and the central black line. The MAPE L , MAPE G and MAPE I in the titles correspond to the MAPE (defined in Equation (19)) of linear, general and iteration models respectively.

Then the L-M number
is calculated by the pressure loss ratio correlation of the Reader-Harris model (Equations (30) and (33)). After is calculated, it is then substituted into the Step 2 to continue the following process until converges.
Finally, the last and are recorded and used as the estimates of the gas and liquid flowrates, respectively.

Comparison of the Direct Fitting and Iteration Method
The gas and liquid flowrate estimates of the linear, general and iteration models are shown in Figures 7 and 8. The black circles denote the linear model results, the blue squares denote the general model results, and the red asterisks denote the iteration model results. The central black line denotes the ideal case when the estimated value is equal to the reference value so that the error is always zero. The upper and lower red lines specify a 10% relative error range so that points within this range have relative errors less than 10%. Similarly, the upper and lower black lines specify a 20% relative error range. The error is represented by the the vertical distance between the test point and the central black line. The MAPE , MAPE and MAPE in the titles correspond to the MAPE (defined in Equation (19)) of linear, general and iteration models respectively. It can be noted from Figures 7 and 8 that the iteration model has the highest accuracy of the three methods, whilst for the direct fitting models, the general model is more accurate than the linear one. This is because as the models become more complicated, the impacts of more parameters are considered, and the contour maps of the differential pressures become closer to the real ones, as shown in Figure 9.
The relative errors of gas and liquid flowrates as a function of the L-M number X are shown in Figure 10a,b respectively, from which it can be noted that the relative errors of gas flowrate increase slowly with X, whilst the relative errors of liquid flowrate decrease rapidly with X. These results are generally in accord with the relative errors distributions of ETV, which will be shown in Section 4.2. It is worth mentioning that the influence of gas Froude number Fr g is not shown in Figure 10, and the trends revealed by the test results are subject to the randomness caused by the small numbers of samples. It can be noted from Figures 7 and 8 that the iteration model has the highest accuracy of the three methods, whilst for the direct fitting models, the general model is more accurate than the linear one. This is because as the models become more complicated, the impacts of more parameters are considered, and the contour maps of the differential pressures become closer to the real ones, as shown in Figure 9.  Figure 10a,b respectively, from which it can be noted that the relative errors of gas flowrate increase slowly with , whilst the relative errors of liquid flowrate decrease rapidly with . These results are generally in accord with the relative errors distributions of ETV, which will be shown in Section 4.2. It is worth mentioning that the influence of gas Froude number is not shown in Figure 10, and the trends revealed by the test results are subject to the randomness caused by the small numbers of samples.  It can be noted from Figures 7 and 8 that the iteration model has the highest accuracy of the three methods, whilst for the direct fitting models, the general model is more accurate than the linear one. This is because as the models become more complicated, the impacts of more parameters are considered, and the contour maps of the differential pressures become closer to the real ones, as shown in Figure 9.  Figure 10a,b respectively, from which it can be noted that the relative errors of gas flowrate increase slowly with , whilst the relative errors of liquid flowrate decrease rapidly with . These results are generally in accord with the relative errors distributions of ETV, which will be shown in Section 4.2. It is worth mentioning that the influence of gas Froude number is not shown in Figure 10, and the trends revealed by the test results are subject to the randomness caused by the small numbers of samples.

Sensitivity Analysis
Although the iteration method can provide accurate results, it is relatively complicated and not very intuitive. Therefore, contour maps like Figure 9b can be used to facilitate our understanding. From Figure 9b, it can be noted that each contour line of d (solid line) and each contour line of d (dashed line) has only one intersection point, the and coordinates of which corresponds to the gas flowrate (or ) and liquid flowrate (or ) respectively, as long as the pressure (or ) and geometric parameters does not change. Therefore, by looking for the intersection points of each d and d combination, we can easily find the corresponding and from Figure 9b, and then calculated the and . In addition to calculating and , the contour maps can also be used to estimate the relative error of the and , and the working principle is shown in Figure 11. The blue solid line (tangent ) in Figure 11 represents the tangent of the contour of the front DP d , whilst the red solid line (tangent ) represents the tangent of the contour of the rear DP d . Therefore, Point 0 corresponds to the working condition and its coordinates correspond to the gas/liquid flowrates. As d and d are measured by the same type of pressure transmitter, then the relative error of d and d are assumed the same, which is = 1 % in this paper. The slopes of and are assumed to be and respectively, and it can be noted From Figure 11 that | | > | |.

Sensitivity Analysis
Although the iteration method can provide accurate results, it is relatively complicated and not very intuitive. Therefore, contour maps like Figure 9b can be used to facilitate our understanding. From Figure 9b, it can be noted that each contour line of dp 1 (solid line) and each contour line of dp 2 (dashed line) has only one intersection point, the x and y coordinates of which corresponds to the gas flowrate Q g (or Fr g ) and liquid flowrate Q l (or Fr l ) respectively, as long as the pressure p (or ρ g ) and geometric parameters does not change. Therefore, by looking for the intersection points of each dp 1 and dp 2 combination, we can easily find the corresponding Fr l and Fr g from Figure 9b, and then calculated the Q l and Q g .
In addition to calculating Q l and Q g , the contour maps can also be used to estimate the relative error of the Q l and Q g , and the working principle is shown in Figure 11. The blue solid line (tangent L 1 ) in Figure 11 represents the tangent of the contour of the front DP dp 1 , whilst the red solid line (tangent L 2 ) represents the tangent of the contour of the rear DP dp 2 . Therefore, Point 0 corresponds to the working condition and its coordinates correspond to the gas/liquid flowrates. As dp 1 and dp 2 are measured by the same type of pressure transmitter, then the relative error of dp 1 and dp 2 are assumed the same, which is σ = 1% in this paper. The slopes of L 1 and L 2 are assumed to be k 1 and k 2 respectively, and it can be noted From Figure 11 that |k 1 | > |k 2 |. 021, 21, 2120 15 of 19 Figure 11. The method of estimating the relative errors of flowrates with the differential pressure contour maps.
From Figure 11, it can be clearly noted that if the measured differential pressure has an absolute error of (d ) = • d ( = 1 or 2), then it will cause an absolute error of and for the gas and liquid flowrate respectively. As the total derivative of differential pressure is (d ) = ∇(d ) • d = d + d , then the expression of ( = 1 or 2) in Figure 11 can be written as: . Similarly, because the tangents and the gradients of the differential pressures ∇(d ) are mutually perpendicular, then the slope of the tangents can be expressed as: = .
From Figure 11, it can also be noted that and are mainly determined by Points 3 and 4, the coordinates of which are the solutions of two sets of binary functions. After some simplifications, the following expressions can be derived: After substituting the expressions of and into Equation (34) and introducing the definitions of relative errors = ⁄ , = ⁄ , we have: From Equation (35), it can be further derived that: Variable Y (Liquid Froude Number Fr l ) Figure 11. The method of estimating the relative errors of flowrates with the differential pressure contour maps.
From Figure 11, it can be clearly noted that if the measured differential pressure has an absolute error of δ(dp i ) = σ·dp i (i = 1 or 2), then it will cause an absolute error of δ g and δ l for the gas and liquid flowrate respectively. As the total derivative of differential pressure is δ(dp i ) = ∇(dp i )·ds = ∂dp i ∂x dx + ∂dp i ∂y dy, then the expression of δ i (i = 1 or 2) in Figure 11 can be written as: δ i = δ(dp i )/ ∂dp i ∂y . Similarly, because the tangents and the gradients of the differential pressures ∇(dp i ) are mutually perpendicular, then the slope of the tangents can be expressed as:k i = ∂dp i ∂x / ∂dp i ∂y . From Figure 11, it can also be noted that δ g and δ l are mainly determined by Points 3 and 4, the coordinates of which are the solutions of two sets of binary functions. After some simplifications, the following expressions can be derived: After substituting the expressions of k i and δ i into Equation (34) and introducing the definitions of relative errors σ g = δ g /Fr g , σ l = δ l /Fr l , we have: ∂x + dp 2 ∂dp 1 ∂x dp 1 ∂dp 2 ∂y + dp 2 From Equation (36), it can be noted that σ l /σ g is only determined by the L-M number X, the dp 1 and dp 2 , and their local gradients. The rightmost part of Equation (36) is based on the assumption that ∇(dp 2 ) ≈ ∇(dp 1 ) = ∇(dp Hom ), from which it can be noted that the ratio between the relative error of liquid and gas flowrate is a constant, which depends on the density ratio ρ l /ρ g and the L-M number X.

Relative Error Distributions
The relative error distributions of an ETV tube are shown in Figure 12, where the color scale bar represents the local relative error of gas/liquid flowrate caused by 1% relative error of dp. From Equation (36), it can be noted that ⁄ is only determined by the L-M number , the d and d , and their local gradients. The rightmost part of Equation (36) is based on the assumption that ∇(d ) ≈ ∇(d ) = ∇(d ), from which it can be noted that the ratio between the relative error of liquid and gas flowrate is a constant, which depends on the density ratio ⁄ and the L-M number .

Relative Error Distributions
The relative error distributions of an ETV tube are shown in Figure 12, where the color scale bar represents the local relative error of gas/liquid flowrate caused by 1% relative error of d . The black dash-dotted line represents the L-M number of one, the black dashed line represents the L-M number of 0.4, whereas the red dash-dotted line represents the L-M number of zero. The results presented in the titles denote the mean value of relative errors of gas/liquid flowrates within the L-M range 0 to 0.4.
From Figure 12, it can be noted that 1% relative error of d will cause roughly 2.8% of relative error for the gas flowrate and 6.4% of relative error for the liquid flowrate. In addition, the relative errors of gas and liquid flowrate generally increase with the L-M number and decreases with the total flowrate , which suggests this device and algorithm are not very suitable for metering multiphase flow with high or low . Under the same conditions, the liquid flowrate error is roughly 2~3 times the gas flowrate error, which is consistent with the predictions of Equation (36).

Comparison of the Classical and the ETV Tubes
By using a similar method, the contour maps of a classical Venturi tube can be generated from the empirical correlations provided in ISO/TR 11583 [2]. With the contour maps of the classical and ETV, it is possible to make a comparison between these two devices and the results are shown in Figure 13, where the color scale bar still represents the local relative error of gas/liquid flowrate caused by 1% relative error of d . The black dashed line represents the L-M number of 0.3, and the results presented in the titles refer to the mean value of relative errors of gas/liquid flowrates within the L-M range 0 to 0.3. From Figure 13 it is notable that the ranges of are both set as 0 to 0.3, and the relative error of all d are set as 1% to make the comparison as fair as possible.
The relative error distributions of a classical Venturi tube are shown in Figure 13a, from which it can be noted that 1% relative error of d will cause roughly 7.1% relative error for the gas flowrate and 20.2% relative error for the liquid flowrate. As a reference, the relative error distributions of an ETV are shown in Figure 13b, from which it can be noted that 1% relative error of d will cause roughly 2.4% relative error for the gas flowrate and 6.2% relative error for the liquid flowrate. Therefore, it can be concluded that the impact of the d error on the ETV is much less than the classical one, which suggests Figure 12. The relative error distributions of an ETV based on the dp 1 and dp 2 contour maps: (a) gas flowrate; (b) liquid flowrate.
From Figure 12, it can be noted that 1% relative error of dp will cause roughly 2.8% of relative error for the gas flowrate and 6.4% of relative error for the liquid flowrate. In addition, the relative errors of gas and liquid flowrate generally increase with the L-M number X and decreases with the total flowrate Q tot , which suggests this device and algorithm are not very suitable for metering multiphase flow with high X or low Q tot . Under the same conditions, the liquid flowrate error is roughly 2~3 times the gas flowrate error, which is consistent with the predictions of Equation (36).

Comparison of the Classical and the ETV Tubes
By using a similar method, the contour maps of a classical Venturi tube can be generated from the empirical correlations provided in ISO/TR 11583 [2]. With the contour maps of the classical and ETV, it is possible to make a comparison between these two devices and the results are shown in Figure 13, where the color scale bar still represents the local relative error of gas/liquid flowrate caused by 1% relative error of dp. The black dashed line represents the L-M number of 0.3, and the results presented in the titles refer to the mean value of relative errors of gas/liquid flowrates within the L-M range 0 to 0.3. From Figure 13 it is notable that the ranges of X are both set as 0 to 0.3, and the relative error of all dp are set as 1% to make the comparison as fair as possible.