Next Article in Journal
Quantifying Circadian Aspects of Mobility-Related Behavior in Older Adults by Body-Worn Sensors—An “Active Period Analysis”
Next Article in Special Issue
Hepatocellular Carcinoma Automatic Diagnosis within CEUS and B-Mode Ultrasound Images Using Advanced Machine Learning Methods
Previous Article in Journal
A Comparative Study on Supervised Machine Learning Algorithms for Copper Recovery Quality Prediction in a Leaching Process
Previous Article in Special Issue
An Estimation Method of Continuous Non-Invasive Arterial Blood Pressure Waveform Using Photoplethysmography: A U-Net Architecture-Based Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Tsinghua Shenzhen International Graduate School, Tsinghua University, Shenzhen 518055, China
2
Shenzhen Institute of Information Technology, Shenzhen 518172, China
3
Research Institute of Tsinghua, Pearl River Delta, Guangzhou 510700, China
4
Wenliu Oil Production Plant of Zhongyuan Oilfield, SINOPEC, Puyang 457001, China
5
Shenzhen LeEngSTAR Technology Co., Ltd., Shenzhen 518055, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(6), 2120; https://doi.org/10.3390/s21062120
Submission received: 29 January 2021 / Revised: 9 March 2021 / Accepted: 12 March 2021 / Published: 17 March 2021
(This article belongs to the Special Issue Machine Learning for Sensing and Healthcare 2020–2021)

Abstract

:
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.

1. 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 d p 1 measurement and a throat section for rear differential pressure d p 2 measurement, as well as a pressure sensor at the inlet and a temperature sensor at the outlet for measuring p and T respectively. The gas and liquid flowrate can be calculated from the above four variables and a direct fitting method has been proposed by Fincke himself [3,4]. Thereafter, multiple variables nonlinear regression [5,6], artificial neural network (ANN), support vector machine (SVM) [7,8] and fluctuating-property-based (FPB) [9,10] methods have been employed to increase its accuracy. However, despite its simple structure, its theoretical model and working principle have not been fully established and explained yet. For example, both the homogeneous flow model [7,8] and separated-flow model [9,10] have been established for the throat section of ETV, but no matter the gas-liquid flow is ideally mixed or separated, there will be no difference between the convergent and throat sections to be exploited for gas/liquid flowrate measurement. Methods such as ANN, SVM, and FPM tend to require large numbers of reference data as the training set, which are usually difficult and expensive to acquire in the oil fields, and their ability to extrapolate beyond their test range still need to be proved. Therefore, practical algorithms of ETV for wet-gas metering are largely lacking in the open literature.
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 (∂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.

2. Theoretical Model

2.1. Convergent Section

When pure gas flows through an ETV tube, and the measured front differential pressure is d p 1 , then its the volume flowrate can be directly calculated by:
Q g o = C d ε 1 β 4 A 2 Δ p 1 ρ g
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 t p , and the “actual” gas flowrate Q g should be calculated by Q g = Q t p / ϕ g , where ϕ g is referred to as the “gas over-reading”. ϕ g is always larger than one because the front differential pressure d p 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:
ϕ g 1 = a 1 X + b 1
where X = Q l Q g ρ l ρ g , whereas Chisholm [23], de Leeuw [24], Reader-Harris [25] proposed a general model like:
ϕ g 1 = 1 + C 1 X + X 2
where:
C 1 = S 1 ρ g ρ l + 1 S 1 ρ l ρ g = ( ρ g ρ l ) n + ( ρ l ρ g ) n
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:
A d p t p = d [ ρ g Q g 2 A α + ρ l Q l 2 A ( 1 α ) ]
By assuming the cross-sectional void fraction α as a constant and integrating Equation (5) from the Venturi inlet to the throat, we have:
d p 1 = 1 2 1 β 4 A 2 [ ρ g Q g 2 α + ρ l Q l 2 ( 1 α ) ]
By using similar method for gas and liquid single-phase, we have:
d p 1 = 1 2 1 β 4 A 2 ρ g Q g 2 α 2 = 1 2 1 β 4 A 2 ρ l Q l 2 ( 1 α ) 2
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 F r 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:
2 K 1 2 d p 1 = F r g 2 + ( 1 S 1 ρ l ρ g + S 1 ρ g ρ l ) F r g F r l + F r l 2
where F r g is the gas Froude number, F r g = Q g A D g D ρ g ρ l ρ g ; F r l is the liquid Froude number, F r l = Q l A D g D ρ l ρ l ρ g , K 1 = C d ε 1 β 4 β 2 g D ( ρ l ρ g ) and A D = π D 2 / 4 . It is worth mentioning that X = F r l / F r g .
Therefore, the gradient of the front differential pressure d p 1 can be written as:
( d p 1 ) = F r g 2 K 1 2 [ 2 + ( 1 S 1 ρ l ρ g + S 1 ρ g ρ l ) X ( 1 S 1 ρ l ρ g + S 1 ρ g ρ l ) + 2 X ]
From Equation (9) it can be noted that the gradient of d p 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 F r g and the parameter K 1 .

2.2. 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:
Δ p g o = τ w g P g o L A = λ L d ρ g 2 ( Q g A ) 2
where τ w g is the shear stress of gas on wall, τ w g = λ ρ g u g o 2 8 ; P g o is the perimeter of gas at pipe wall, P g o = π 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:
Q t p 2 = 1 λ L / d A 2 Δ p 2 ρ g
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]:
A l d p t p d z = τ w l P l S F
for liquid, and:
A g d p t p d z = τ w g P g + S F
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 χ = 4 P i / A i are assumed the same for the gas and liquid, then divide Equation (12) by Equation (13), we have:
Z 2 = A l d p 2 d z + S F A g d p 2 d z S F = ρ l ρ g ( u l u g ) 2
After substituting Equation (14) into Equation (13) and cancelling S F , we have:
d p t p d z = λ L d ρ g 2 ( Q g A g ) 2 Z 2 + A g A l A g A l + 1
From Equations (15) and (10), it can be derived that:
ϕ g 2 = 1 + ( Z + 1 Z ) X + X 2
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 d p 2 can also be written as similar form as Equation (9):
( d p 2 ) = F r g 2 K 2 2 [ 2 + ( 1 S 2 ρ l ρ g + S 2 ρ g ρ l ) X ( 1 S 2 ρ l ρ g + S 2 ρ g ρ l ) + 2 X ]
where K 2 = 1 λ L / d β 2 g D ( ρ l ρ g ) .

2.3. 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 ( d p 1 ) and ( d p 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 d p 1 and rear differential pressure d p 2 are shown in Figure 2a. From Figure 2a, it can be noted that there is always an intersection angle between ( d p 1 ) and ( d p 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 F r 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 .

3. 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 m3/h to 650 m3/h, and the liquid flow test range is from 0 to 8 m3/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 m3/h and 0~4.1 m3/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, where the dash-dotted line represents the Lockhart–Martinelli (L-M) number of 1.0 whereas the dashed line represents the L-M number of 0.4. Therefore, it is notable that the L-M range of the test points is 0~0.4, and the flow patterns are mainly annular and slug flows.
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:
MAE = 1 n i = 1 n | y ^ i y i |
MAPE = 100 % n i = 1 n | y ^ i y i y i |
where y ^ 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 m a x ) whereas MAPE is mainly used for the errors of flowrates (e.g., Q g and Q l ).

3.1. Direct Fitting Method

3.1.1. Linear Model

Simultaneous solution of the linear over-reading Equation (2) of the convergent and throat sections leads to the linear models as follows:
[ Q g t 1 Q g t 2 ] = [ a 1 ρ l ρ g b 1 a 2 ρ l ρ g b 2 ] [ Q l Q g ]
Matrix inversion can be applied to Equation (20) directly and the result is:
[ Q l Q g ] = [ a 1 ρ g ρ l b 1 ρ g ρ l a 2 b 2 ] [ Q g t 1 Q g t 2 ]
where a 1 = b 2 a 1 b 2 a 2 b 1 , a 2 = a 2 a 1 b 2 a 2 b 1 , b 1 = b 1 a 1 b 2 a 2 b 1 , b 2 = a 1 a 1 b 2 a 2 b 1 .
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 g t 1 and Q g t 2 , 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 F r g ); (3) fix the “bad values” or “blind spots” of some complicated algorithms.

3.1.2. 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]:
ϕ g 1 = C 3 1 + C 1 X + X 2
ϕ g 2 = C 4 1 + C 2 X + X 2
By solving the above two equations simultaneously, we can get:
Q g = a 1 Q g t p 1 2 + a 2 Q g t p 2 2 ± ( a 1 Q g t p 1 2 + a 2 Q g t p 2 2 ) 2 ( a 3 Q g t p 1 2 + a 4 Q g t p 2 2 ) 2
Q l = a 1 Q l t p 1 2 + a 2 Q l t p 2 2 ± ( a 1 Q l t p 1 2 + a 2 Q l t p 2 2 ) 2 ( a 3 Q l t p 1 2 + a 4 Q l t p 2 2 ) 2
where a 1 = C 2 2 C 3 2 ( C 1 C 2 ) , a 2 = C 1 2 C 4 2 ( C 1 C 2 ) , a 3 = 1 C 3 2 ( C 1 C 2 ) , a 4 = 1 C 4 2 ( C 1 C 2 ) .
Similar to the linear model, variable substitutions can also be implemented to prevent the errors introduced during the solving process.

3.2. Iteration Method

3.2.1. 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 over-reading 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 over-reading model. Currently, the most common over-reading models for wet gas measurement includes:
  • 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:
C Hom = ρ g ρ l + ρ l ρ g
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.
2.
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:
n = { 0.606 ( 1 e 0.746 F r g ) if   F r g 1.5 0.41 if   0.5 F r g < 1.5
It is notable that the de Leeuw model needs to know F r g first and calculates n and C 1 latter, so ϕ ^ g should be determined through an iteration method.
3.
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:
n = max [ 0.583 0.18 β 2 0.578 exp ( 0.8 F r g H ) , 0.392 0.18 β 2 ]
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:
C d = 1 0.0463 exp ( 0.05 F r g ) min ( 1 , X 0.016 )
Because C d 1 is assumed for other over-reading models, so the ϕ g RH calculated by Equations (3) and (4) should be divided by C d to facilitate its comparison with other models. The data fitting results of the homogeneous model, de Leeuw model and Reader-Harris model are shown in Figure 4, where the MAE-H, MAE-D and MAE-R in the titles correspond to the MAE (defined in Equation (18)) of Homogeneous, de Leeuw and Reader-Harris models respectively.
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.

3.2.2. 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 ϖ d p , where d ϖ represents the pressure loss across a classical Venturi tube, d p 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 p 2 and front differential pressure d p 1 respectively, as shown in Figure 1.
The pressure loss ratio should be normalized first:
Y Y max = Δ p 2 Δ p 1 Δ p 2 Δ p 1 | dry ( Δ p 2 Δ p 1 Δ p 2 Δ p 1 | dry ) | X = 0.3
and then data fitted in the following form:
Y Y max = 1 exp [ a X c exp ( b F r gas H ) ]
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]:
Y Y max = 1 exp [ 35 X 0.75 exp ( 0.28 F r gas H ) ]
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:
Y Y max = 1 exp [ 5.5883 X 0.439 exp ( 0.2586 F r gas ) ]
where Δ p 2 Δ p 1 | dry and Δ p 2 Δ p 1 | X = 0.4 are the minimum and maximum of the test data respectively, because the test data covered the two extreme cases of dry gas and X = 0.4 .
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.

3.2.3. 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 F r g to be known, this algorithm assumes X ^ = 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:
  • First, the initial values of the L-M number X ^ is assumed.
  • Then the gas over-reading ϕ ^ g is estimated according to the homogenous model (Equation (26)).
  • Then the gas/liquid flowrates ( Q ^ g and Q ^ l ) are calculated from the definitions of gas over-reading ϕ ^ g and L-M number X ^ respectively.
  • Then F r ^ g and F r ^ l are calculated by its definition.
  • Then the gas over-reading ϕ ^ g is calculated by the de Leeuw model (Equations (3), (4) and (27)). After ϕ ^ g is calculated, ϕ ^ g and X ^ are substituted into the Step 3 to continue the following process until ϕ ^ g converges.
  • Then the L-M number X ^ is calculated by the pressure loss ratio correlation of the Reader-Harris model (Equations (30) and (33)). After X ^ is calculated, it is then substituted into the Step 2 to continue the following process until X ^ converges.
Finally, the last Q ^ g and Q ^ l are recorded and used as the estimates of the gas and liquid flowrates, respectively.

3.3. Comparison of the Direct Fitting and Iteration Method

The gas and liquid flowrate estimates of the linear, general and iteration models are shown in Figure 7 and Figure 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.
It can be noted from Figure 7 and Figure 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 F r 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.

4. Sensitivity Analysis and Error Ddistributions

4.1. 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 p 1 (solid line) and each contour line of d p 2 (dashed line) has only one intersection point, the x and y coordinates of which corresponds to the gas flowrate Q g (or F r g ) and liquid flowrate Q l (or F r 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 d p 1 and d p 2 combination, we can easily find the corresponding F r l and F r 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 d p 1 , whilst the red solid line (tangent L 2 ) represents the tangent of the contour of the rear DP d p 2 . Therefore, Point 0 corresponds to the working condition and its coordinates correspond to the gas/liquid flowrates. As d p 1 and d p 2 are measured by the same type of pressure transmitter, then the relative error of d p 1 and d p 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 | .
From Figure 11, it can be clearly noted that if the measured differential pressure has an absolute error of δ ( d p i ) = σ · d p 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 δ ( d p i ) = ( d p i ) · d s = d p i x d x + d p i y d y , then the expression of δ i ( i = 1   or   2 ) in Figure 11 can be written as: δ i = δ ( d p i ) / d p i y . Similarly, because the tangents and the gradients of the differential pressures ( d p i ) are mutually perpendicular, then the slope of the tangents can be expressed as: k i = d p i x / d p 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:
{ δ g = δ 1 + δ 2 k 1 k 2 δ l = k 2 δ 1 + k 1 δ 2 k 1 k 2
After substituting the expressions of k i and δ i into Equation (34) and introducing the definitions of relative errors σ g = δ g / F r g , σ l = δ l / F r l , we have:
{ σ g = σ F r g d p 1 d p 2 y + d p 2 d p 1 y d p 1 x d p 2 y d p 2 x d p 1 y σ l = σ F r l d p 1 d p 2 x + d p 2 d p 1 x d p 1 x d p 2 y d p 2 x d p 1 y
From Equation (35), it can be further derived that:
σ l σ g = 1 X d p 1 d p 2 x + d p 2 d p 1 x d p 1 d p 2 y + d p 2 d p 1 y ( ρ l ρ g + ρ g ρ l ) + 2 / X ( ρ l ρ g + ρ g ρ l ) + 2 X
From Equation (36), it can be noted that σ l / σ g is only determined by the L-M number X , the d p 1 and d p 2 , and their local gradients. The rightmost part of Equation (36) is based on the assumption that ( d p 2 ) ( d p 1 ) = ( d p 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 .

4.2. 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 p . 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 p 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).

4.3. 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 p . 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 d p 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 p 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 p 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 p error on the ETV is much less than the classical one, which suggests the ETV may be more accurate than the classical one under the same conditions. Meanwhile, the average velocity is much higher in the throat section of an ETV, which means the ETV can be designed to be more compact without undermining the signal intensity.

5. Conclusions

In this paper, theoretical models are established for the convergent and throat sections of an ETV. Several direct fitting and iteration flowrate algorithms are proposed and compared with the existing ones. The gradients of front and rear differential pressures are derived and the relationship between the differential pressure gradient and the flowrate relative error is studied analytically. Finally, the relative error distributions of an ETV are obtained and compared with the ones of a classical Venturi tube. The following important conclusions can be obtained through this research:
  • The iteration algorithm is more accurate than the direct fitting algorithm because it considers the influence of more parameters and therefore, its contour maps of differential pressures are closer to the reality than those of the other algorithms.
  • The gas flowrate error of the ETV iteration algorithm increases with the liquid content X whilst the liquid flowrate error of the ETV iteration algorithm decreases with the liquid content X .
  • The relative errors of liquid flowrates tend to be 2 to 3 times larger than those of the gas flowrates, which can be explained by the theoretical model and is in good agreement with the experimental results.
  • The ETV tube tends to be more accurate than the classical one. Additionally, it can be designed more compactly under the same signal intensity due to its significantly higher velocity in the throat section.

Author Contributions

Conceptualization, H.X. and P.Y.; methodology, H.X.; software, H.X.; validation, H.X. and P.Y.; formal analysis, H.X.; writing—original draft preparation, H.X.; writing—review and editing, M.Z., G.W., H.Z. and E.W.; supervision, M.Z. and X.Z.; project administration, Y.L.; funding acquisition, Y.L. and X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received funding from National Natural Science Foundation of China (Grant No. 61571252) and China National Key Research Scheme (Grant No. 2016YFC0303706).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors gratefully acknowledge the Chengdu Verification Branch for Oil & Gas Large Flowrate Measurement station (CVB) for providing the experimental facility and Shenzhen LeEngSTAR Technology Co., Ltd. for providing the wet gas sensors and measurement units.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. ISO 5167-4. Measurement of Fluid Flow by Means of Pressure Differential Devices Inserted in Circular Cross-Section Conduits Running Full—Part 4: Venturi Tubes; International Organization for Standardization: Geneva, Switzerland, 2003; pp. 1–24. [Google Scholar]
  2. ISO/TR 11583:2012(E). Measurement of Wet Gas Flow by Means of Pressure Differential Devices Inserted in Circular Cross-Section Conduits; International Organization for Standardization: Geneva, Switzerland, 2012. [Google Scholar]
  3. Fincke, J.R. Method and System for Measuring Multiphase Flow Using Multiple Pressure Differentials; Bechtel BWXT Idaho: Scoville, ID, USA, 2001; pp. 1–13. [Google Scholar]
  4. Fincke, J.R.; Ronnenkamp, C.; Kruse, D.; Krogue, J.; Householder, D. Performance characteristics of an extended throat flow nozzle for the measurement of high void fraction multi-phase flows. In Proceedings of the 4th International Symposium Fluid Flow Measurement, Denver, CO, USA, 28–30 June 1999. [Google Scholar]
  5. Xu, Y.; Yuan, C.; Long, Z.; Zhang, Q.; Li, Z.; Zhang, T. Research the Wet Gas Flow Measurement Based on Dual-Throttle Device. Flow Meas. Instrum. 2013, 34, 68–75. [Google Scholar] [CrossRef]
  6. Zhang, Q.; Xu, Y.; Zhang, T. Wet gas metering based on dual differential pressure of long throat Venturi tube. Trans. Tianjin Univ. 2012, 45, 147–153. [Google Scholar]
  7. Xu, L.; Li, H.; Xu, J. Wet gas metering by using a long-throat Venturi meter. In Seventh International Symposium on Instrumentation and Control Technology: Sensors and Instruments, Computer Simulation, and Artificial Intelligence; International Society for Optics and Photonics: Beijing, China, 2008; pp. 1–9. [Google Scholar]
  8. Xu, L.; Zhou, W.; Li, X.; Tang, S. Wet gas metering using a revised Venturi meter and soft-computing approximation techniques. IEEE Trans. Instrum. Meas. 2010, 60, 947–956. [Google Scholar] [CrossRef]
  9. Xu, L.; Zhou, W.; Li, X.; Wang, M. Wet-gas flow modeling for the straight section of throat-extended Venturi meter. IEEE Trans. Instrum. Meas. 2011, 60, 2080–2087. [Google Scholar]
  10. Zhou, W.; Xu, L.; Li, X. Wet gas flow modeling for the straight section of throat-extended Venturi meter. In Proceedings of the 2010 IEEE Instrumentation & Measurement Technology Conference Proceedings, Austin, TX, USA, 3–6 May 2010. [Google Scholar]
  11. Almutairi, Z.; Al-Alweet, F.M.; Alghamdi, Y.A.; Almisned, O.A.; Alothman, O.Y. Investigating the Characteristics of Two-Phase Flow Using Electrical Capacitance Tomography (ECT) for Three Pipe Orientations. Processes 2020, 8, 51. [Google Scholar] [CrossRef] [Green Version]
  12. Xu, Z.; Wu, F.; Yang, X.; Li, Y. Measurement of Gas-Oil Two-Phase Flow Patterns by Using CNN Algorithm Based on Dual ECT Sensors with Venturi Tube. Sensors 2020, 20, 1200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Roshani, M.; Phan, G.T.T.; Ali, P.J.M.; Roshani, G.H.; Hanus, R.; Duong, T.; Corniani, E.; Nazemi, E.; Kalmoun, E.M. Evaluation of flow pattern recognition and void fraction measurement in two phase flow independent of oil pipeline’s scale layer thickness. Alex. Eng. J. 2021, 60, 1955–1966. [Google Scholar] [CrossRef]
  14. Fang, L.; Zeng, Q.; Wang, F.; Faraj, Y.; Zhao, Y.; Lang, Y.; Wei, Z. Identification of two-phase flow regime using ultrasonic phased array. Flow Meas. Instrum. 2020, 72, 101726. [Google Scholar] [CrossRef]
  15. Zhang, Y.; Azman, A.N.; Xu, K.W.; Kang, C.; Kim, H.B. Two-phase flow regime identification based on the liquid-phase velocity information and machine learning. Exp. Fluids 2020, 61, 1–16. [Google Scholar] [CrossRef]
  16. Johansen, E.S.; Hall AR, W.; Ünalmis Ö, H.; Rodriguez, D.J.; Vera, A.; Ramakrishnan, V. A prototype wet-gas and multiphase flowmeter. In Proceedings of the 25th International North Sea Flow Measurement Workshop, Oslo, Norway, 16–18 October 2007. [Google Scholar]
  17. Johansen, E.S.; Unalmis, O.H.; Lievois, J. WET-GAS FLOWMETER. U.S. Patent 7654155B2, 2 February 2007. Weatherford/Lamb, Inc.: Houston, TX, USA. [Google Scholar]
  18. Arsalan, M.; Karimi, M.A.; Shamim, A. Multiphase Flow Meter Combining Extended Throat Venturi with Microwave Resonators; WIPO: New York, NY, USA, 2020; pp. 1–26. [Google Scholar]
  19. ISO 12213-2:2006. Natural Gas—Calculation of Compression Factor—Part 2: Calculation Using Molar-Composition Analysis; International Organization for Standardization: Geneva, Switzerland, 2006; pp. 1–32. [Google Scholar]
  20. Murdock, J. Two-phase flow measurement with orifices. J. Basic Eng. 1962, 84, 419–432. [Google Scholar] [CrossRef]
  21. Bizon, E. Two-Phase Flow Measurement with Sharp-Edged Orifices and Venturis; Atomic Energy of Canada Ltd.: Chalk River, ON, Canada, 1965. [Google Scholar]
  22. Lin, Z. Two-phase flow measurements with sharp-edged orifices. Int. J. Multiph. Flow 1982, 8, 683–693. [Google Scholar] [CrossRef]
  23. Chisholm, D. A theoretical basis for the Lockhart-Martinelli correlation for two-phase flow. Int. J. Heat Mass Transf. 1967, 10, 1767–1778. [Google Scholar] [CrossRef]
  24. De Leeuw, R. Liquid correction of Venturi meter readings in wet gas flow. In Proceedings of the North Sea Flow Measurement Workshop 1997, Kristiansand, Norway, 27–30 October 1997. [Google Scholar]
  25. Reader-Harris, M.; Graham, E. An improved model for Venturi-tube over-reading in wet gas. In Proceedings of the 27th International North Sea Flow Measurement Workshop, Tonsberg, Norway, 20–23 October 2009. [Google Scholar]
  26. Yu, P.; Xu, Y.; Zhang, T.; Zhu, Z.; Ba, X.; Li, J.; Qin, Z. A study on the modeling of static pressure distribution of wet gas in Venturi. AIChE J. 2015, 61, 699–708. [Google Scholar] [CrossRef]
  27. Chisholm, D.; Sutherland, L. Prediction of pressure gradients in pipeline systems during two-phase flow. Proc. Inst. Mech. Eng. 1969, 184, 24–32. [Google Scholar] [CrossRef]
  28. Bjørnera, M.; Fosbølb, P.; Lisberga, M.; Lisberga, H. Generalized Wet Gas Venturi Meter Correlations: Assessment and Improvement. In Proceedings of the 35th North Sea Flow Measurement Workshop, Tønsberg, Norway, 24–26 October 2017; pp. 1–25. [Google Scholar]
  29. Chang, H.; Duan, J. Natural gas measurement technology system and its prospect in China. Nat. Gas Ind. B 2020, 7, 370–379. [Google Scholar] [CrossRef]
  30. Duan, J. The Technical Status and Prospect of Natural Gas Flow Traceability System in China. In Proceedings of the FLOMEKO 2019, Lisbon, Portugal, 26–28 June 2019. [Google Scholar]
  31. Collins, A.; Clark, S. Evolution of Wet Gas Venturi Metering and Wet Gas Correction Algorithms. Meas. Control 2013, 46, 15–20. [Google Scholar] [CrossRef] [Green Version]
  32. Collins, A.; Tudge, M.; Wade, C. Evaluating and improving wet gas corrections for horizontal Venturi meters. In Proceedings of the 33rd International North Sea Flow Measurement Workshop, Tønsberg, Norway, 20–23 October 2015; pp. 1–21. [Google Scholar]
  33. Reader-Harris, M.; Hodges, D.; Gibson, J. Venturi-tube performance in wet gas using different test fluids. In Proceedings of the 24th International North Sea Flow Measurement Workshop, St. Andrews, UK, 24–27 October 2006. [Google Scholar]
Figure 1. Schematic diagram of an extended-throat Venturi (ETV) tube.
Figure 1. Schematic diagram of an extended-throat Venturi (ETV) tube.
Sensors 21 02120 g001
Figure 2. Schematic diagram of the working principle: (a) the contours of d p 1 and d p 2 in the ideal case; (b) the variations of index n of several common over-reading models in the real case. (The filled area corresponds to the Reader-Harris model under β = 0.4 ~ 0.75 , where the black area represents the oil whilst the blue one represents the water. n = 0.5 corresponds to the homogenous model, n = 0.25 corresponds to the Chisholm orifice plate model, whereas n = 0 corresponds to the stratified model or dense phase model ( ρ g = ρ l )) [28].
Figure 2. Schematic diagram of the working principle: (a) the contours of d p 1 and d p 2 in the ideal case; (b) the variations of index n of several common over-reading models in the real case. (The filled area corresponds to the Reader-Harris model under β = 0.4 ~ 0.75 , where the black area represents the oil whilst the blue one represents the water. n = 0.5 corresponds to the homogenous model, n = 0.25 corresponds to the Chisholm orifice plate model, whereas n = 0 corresponds to the stratified model or dense phase model ( ρ g = ρ l )) [28].
Sensors 21 02120 g002
Figure 3. Schematic diagram of the multiphase flow experimental facility and test data: (a) experimental facility; (b) test data.
Figure 3. Schematic diagram of the multiphase flow experimental facility and test data: (a) experimental facility; (b) test data.
Sensors 21 02120 g003
Figure 4. The prediction results of several common over-reading equations for the gas over-reading of an ETV.
Figure 4. The prediction results of several common over-reading equations for the gas over-reading of an ETV.
Sensors 21 02120 g004
Figure 5. Schematic diagram of a classical Venturi.
Figure 5. Schematic diagram of a classical Venturi.
Sensors 21 02120 g005
Figure 6. The estimates of the original and improved Reader-Harris correlation for the normalized differential pressure ratio of an ETV.
Figure 6. The estimates of the original and improved Reader-Harris correlation for the normalized differential pressure ratio of an ETV.
Sensors 21 02120 g006
Figure 7. Gas flowrate estimates of the linear, general and iteration models: (a) absolute error; (b) relative error. (The central black line denotes the ideal case with zero error, the upper and lower red lines denote the 10% relative error range, and the upper and lower black lines denote the 20% relative error range. The vertical distance between the point and the central black line denote the error.
Figure 7. Gas flowrate estimates of the linear, general and iteration models: (a) absolute error; (b) relative error. (The central black line denotes the ideal case with zero error, the upper and lower red lines denote the 10% relative error range, and the upper and lower black lines denote the 20% relative error range. The vertical distance between the point and the central black line denote the error.
Sensors 21 02120 g007
Figure 8. Liquid flowrate estimates of the linear, general and iteration models: (a) absolute error; (b) relative error. (The central black line denotes the ideal case with zero error, the upper and lower red lines denote the 10% relative error range, and the upper and lower black lines denote the 20% relative error range. The vertical distance between the point and the central black line denote the error.
Figure 8. Liquid flowrate estimates of the linear, general and iteration models: (a) absolute error; (b) relative error. (The central black line denotes the ideal case with zero error, the upper and lower red lines denote the 10% relative error range, and the upper and lower black lines denote the 20% relative error range. The vertical distance between the point and the central black line denote the error.
Sensors 21 02120 g008
Figure 9. The comparisons of differential pressure (DP) contours of the linear model and the iteration model: (a) linear model; (b) iteration model.
Figure 9. The comparisons of differential pressure (DP) contours of the linear model and the iteration model: (a) linear model; (b) iteration model.
Sensors 21 02120 g009
Figure 10. The relative errors of gas and liquid flowrates as a function of the Lockhart-Martinelli (L-M) number X : (a) gas flowrate; (b) liquid flowrate.
Figure 10. The relative errors of gas and liquid flowrates as a function of the Lockhart-Martinelli (L-M) number X : (a) gas flowrate; (b) liquid flowrate.
Sensors 21 02120 g010
Figure 11. The method of estimating the relative errors of flowrates with the differential pressure contour maps.
Figure 11. The method of estimating the relative errors of flowrates with the differential pressure contour maps.
Sensors 21 02120 g011
Figure 12. The relative error distributions of an ETV based on the d p 1 and d p 2 contour maps: (a) gas flowrate; (b) liquid flowrate.
Figure 12. The relative error distributions of an ETV based on the d p 1 and d p 2 contour maps: (a) gas flowrate; (b) liquid flowrate.
Sensors 21 02120 g012
Figure 13. The comparison of the classical and the ETV: (a) the gas relative error distributions of the classical Venturi tube from the ISO/TR 11583; (b) the liquid relative error distributions of the classical Venturi tube from the ISO/TR 11583; (c) the gas relative error distributions of the ETV from the experiments; (d) the liquid relative error distributions of the ETV from the experiments.
Figure 13. The comparison of the classical and the ETV: (a) the gas relative error distributions of the classical Venturi tube from the ISO/TR 11583; (b) the liquid relative error distributions of the classical Venturi tube from the ISO/TR 11583; (c) the gas relative error distributions of the ETV from the experiments; (d) the liquid relative error distributions of the ETV from the experiments.
Sensors 21 02120 g013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xue, H.; Yu, P.; Zhang, M.; Zhang, H.; Wang, E.; Wu, G.; Li, Y.; Zheng, X. A Wet Gas Metering System Based on the Extended-Throat Venturi Tube. Sensors 2021, 21, 2120. https://doi.org/10.3390/s21062120

AMA Style

Xue H, Yu P, Zhang M, Zhang H, Wang E, Wu G, Li Y, Zheng X. A Wet Gas Metering System Based on the Extended-Throat Venturi Tube. Sensors. 2021; 21(6):2120. https://doi.org/10.3390/s21062120

Chicago/Turabian Style

Xue, Haobai, Peining Yu, Maomao Zhang, Haifeng Zhang, Encheng Wang, Guozhu Wu, Yi Li, and Xiangyuan Zheng. 2021. "A Wet Gas Metering System Based on the Extended-Throat Venturi Tube" Sensors 21, no. 6: 2120. https://doi.org/10.3390/s21062120

APA Style

Xue, H., Yu, P., Zhang, M., Zhang, H., Wang, E., Wu, G., Li, Y., & Zheng, X. (2021). A Wet Gas Metering System Based on the Extended-Throat Venturi Tube. Sensors, 21(6), 2120. https://doi.org/10.3390/s21062120

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop