Abstract:
Systems and methods are disclosed for measuring densities and flow rates of gas-liquid fluid mixtures. In the systems and methods, the fluid mixture is caused to exhibit swirling flow as it flows through a conduit that includes a constriction, a first pressure difference is measured between two vertically-spaced measurement positions in the conduit, a second pressure difference is measured between two horizontally-spaced measurement positions in the conduit, the first horizontally-spaced measurement position being at the constriction region and the second horizontally-spaced measurement position being upstream or downstream of the constriction region, and one or more of the pressure differences is used to determine a density or a flow rate of the gas-liquid fluid mixture.

Description:
FIELD OF THE INVENTION  
       [0001]    The present invention relates to a method of measuring the density of fluid flow, and more particularly to a method of measuring the flow rate of a gas-liquid fluid mixture. 
       BACKGROUND OF THE INVENTION 
       [0002]    The determination of gas and liquid flow rates and mixture densities in gas-liquid fluid mixtures is important in the oil and gas industry. 
         [0003]    An example of an apparatus for measuring such flow rates is Schlumberger&#39;s Vx™ system (see e.g. I. Atkinson, M. Berard, B.-V. Hanssen, G. Ségéral, 17 th  International North Sea Flow Measurement Workshop, Oslo, Norway 25-28 Oct. 1999 “New Generation Multiphase Flowmeters from Schlumberger and Framo Engineering AS”) which comprises a vertically mounted Venturi flow meter, a dual energy gamma-ray hold up measuring device and associated processors. This system allows the simultaneous calculation of gas, water and oil volumetric flow rates in multi phase flows. 
         [0004]    Although providing proven performance, the Vx™ system and other conventional multiphase flow meters are relatively expensive, which tends to preclude their application in “brown” field sites (i.e. oil and gas wells where capacity has fallen below about 1000 barrels/day (0.0018 m 3 /sec)) and other low hydrocarbon producers. However, such sites probably account for around 2-3 million oil and gas wells worldwide. 
       SUMMARY OF THE INVENTION  
       [0005]    The present invention is at least partly based on the realisation that it is possible to measure a vertical pressure difference across a substantially horizontally conduit conveying a fluid flow, and to use that pressure difference to determine the fluid density. The fluid density can then be used to calculate a flow rate, such as the liquid or the gas flow rate in the flow of a gas-liquid fluid mixture, or for other purposes. 
         [0006]    Thus, in general terms, the present invention provides a method and a corresponding apparatus for determining the density of a fluid, such as a gas-liquid fluid mixture. 
         [0007]    A first aspect of the present invention provides an apparatus for measuring flow properties of a multiphase mixture flowing in a conduit, including the density of the multiphase mixture, where the multiphase mixture comprises a gas phase and a liquid phase, the apparatus comprising: 
         [0008]    to induce the multiphase mixture to exhibit swirling flow through the conduit: 
         [0009]    a Venturi disposed downstream of the swirl element; and 
         [0010]    a differential pressure meter for measuring a pressure difference across the phase separated multiphase mixture between two vertically-spaced measurement positions in a throat of the Venturi. 
         [0011]    Measuring the pressure difference in a direction perpendicular to the flow direction has an advantage in that the frictional pressure loss caused by the measurement will generally be zero. Frictional pressure losses can be difficult to predict, in particular for multiphase flows, and can introduce significant errors in density measurements. 
         [0012]    Density is an important parameter in process control. In multiphase flow it can be used (along with the known or measured single phase densities) to determine the fractional hold up. This is important if the flow proceeds e.g. to a device that can only handle a certain gas fraction, such as a pump. 
         [0013]    Unlike the Vx™ system, which requires gamma-ray-determined hold-ups to calculate flow rates, the method for measuring a flow rate uses the first pressure difference (which is effectively a measurement of the mixture density or gravitational pressure head across the conduit) in the determination of flow rates. As a pressure difference measurement can be obtained using relatively inexpensive, conventional and robust technology, a multiphase flow meter that applies the method may be suitable for use in brown field sites. A further advantage relative to the Vx™ system is that health and safety issues pertaining to the use of gamma-ray sources can be avoided. 
         [0014]    In general, in order to calculate a fluid density from the first pressure difference, and also to calculate a flow rate using that density and the second pressure difference, the fluid mixture should exhibit a predetermined type of flow at the measurement positions. For example, if it is known that the mixture is e.g. stratified, churning, or homogenised, appropriate adjustments can be made to relations used in the calculation of the fluid density and flow rate. Preferably, however, the mixture is conditioned to exhibit swirling flow, which separates the liquid from the gas at the measurement positions. For example, the conduit may have a swirl element, such as a helical insert or vane assembly, for inducing the mixture to exhibit swirling flow at these positions. The swirl element may include one or more spiral-shaped members extending along the conduit in the direction of fluid flow. Preferably, the spiral shaped members are positioned at the wall of the conduit and, when viewed along the axis the conduit, leave a central core of the conduit unimpeded (i.e. they do not extend radially inwards as far as the central axis of the conduit). Alternatively, the swirl element may be formed by a tangential flow inlet to the conduit. 
         [0015]    An advantage of swirling flow is that it is relatively easy to induce and sustain (unlike stratified or homogenised flow which may be unstable over typical measurement distances). Further, modelling the characteristics of swirling flow through a Venturi is relatively straightforward, compared to e.g. modelling stratified or churning flow. Also, swirling flow is symmetrical about the flow axis, resulting in the second pressure difference being independent of angular orientation. 
         [0016]    Inducing the mixture to exhibit swirling flow separates the liquid and gas phases of the mixture. The swirling flow causes the liquid of the mixture to be displaced to the wall of the conduit, e.g. to form an annulus of liquid adjacent the wall of the conduit, leaving a gas core at the centre of the conduit. 
         [0017]    Swirling flow in the constriction region will have increased centrifugal acceleration relative to swirling flow outside the constriction region (e.g. at the entrance to the constriction region). This increased centrifugal acceleration can enhance the displacement of liquid to the wall of the conduit. Thus preferably, the swirling flow in the constriction region is induced by provoking swirling flow in the conduit upstream of the constriction region. 
         [0018]    Within the constriction region, the separated liquid layer tends to flow more slowly than the gas, which increases the liquid hold up, making it easier to investigate the properties of the liquid. Further, enhanced centrifugal separation at the constriction region can reduce the amount of entrained gas in the liquid, improving estimates of gas or liquid hold-up. 
         [0019]    Typically, the constriction region is provided by a Venturi. The constriction region may be the throat of the Venturi. The second horizontally-spaced measurement position can be at the inlet to the Venturi. 
         [0020]    The method may further include the step of measuring the absolute pressure of the flow, e.g. at one of the measurement positions, but preferably on a transverse cross-section of the conduit containing the vertically-spaced measurement positions. The absolute pressure may also then be used in the determination of the flow rate. The method may further include the step of measuring the temperature of the flow e.g. at one of the measurement positions, but preferably on the transverse cross-section of the conduit containing the vertically-spaced measurement positions. The temperature may also be used in the determination of the flow rate. 
         [0021]    The method may further including the step of measuring a pressure difference between two further vertically-spaced measurement positions in the conduit (e.g. at facing upper and lower wall portions of the conduit) spaced upstream or downstream of the vertically-spaced measurement positions where the first pressure difference is measured; 
         [0022]    wherein the first, second and third pressure differences are used to determine said flow rate. Typically, and conveniently, the further vertically-spaced measurement positions can be on a conduit transverse cross-section which also contains one of the horizontally-spaced measurement positions. In particular, when the vertically-spaced measurement positions for the first pressure difference are on a conduit transverse cross-section which contains one of the horizontally-spaced measurement positions, the further vertically-spaced measurement positions for the third pressure difference can be on a conduit transverse cross-section which contains the other of the horizontally-spaced measurement positions. In this case, that other horizontally-spaced measurement position can coincide with one of the further vertically-spaced measurement positions. 
         [0023]    The additional pressure difference, which like the first pressure difference is effectively the measurement of the mixture density or gravitational pressure head across the conduit, can be used to compensate for slip in the flow of the mixture when determining the flow rate. 
         [0024]    To compensate for irregularities in the flow and to reduce the effect of noise in the measurements, the measured pressure difference(s) may be time-averaged pressure difference(s). 
         [0025]    The method may further include the step of measuring the absolute pressure of the flow on a transverse cross-section of the conduit containing the further vertically-spaced measurement positions and preferably at one of the further vertically-spaced measurement positions, wherein this absolute pressure is also used in the determination of the flow rate. The method may further include the step of measuring the temperature of the flow on the transverse cross-section of the conduit containing the further vertically-spaced measurement positions and preferably at one of the further vertically-spaced measurement positions, wherein this temperature is also used in the determination of the flow rate. 
         [0026]    The liquid of the mixture may comprise oil and/or water. The gas may comprise natural gas. Thus the gas-liquid fluid mixture may be a mixture of natural gas, condensate and optionally water. 
         [0027]    A second aspect of the present invention provides a method for determining flow properties of a multiphase mixture flowing through a conduit, the flow properties including a density of the multiphase mixture, where the multiphase mixture comprising a gas phase and a liquid phase, the method comprising: 
         [0028]    generating a swirling flow of the multiphase mixture in the conduit; 
         [0029]    passing the swirling flow of the multiphase mixture through a Venturi; and 
         [0030]    measuring a differential pressure between two vertically-spaced measurement positions across a throat of the Venuri. 
         [0031]    Other aspects of the present invention provide an apparatus for providing measurements useable in determining a flow rate of a gas-liquid fluid mixture, the apparatus including: 
         [0032]    a conduit through which the gas-liquid fluid mixture can flow, the conduit extending substantially horizontally and having a constriction region having a reduced conduit cross-section; 
         [0033]    a first pressure meter for measuring a first pressure difference between two vertically-spaced measurement positions in the conduit; and 
         [0034]    a second pressure meter for measuring a second pressure difference between two horizontally-spaced measurement positions in the conduit, the first horizontally-spaced measurement position being at the constriction region and the second horizontally-spaced measurement position being upstream or downstream of the constriction region; 
         [0035]    wherein the first and the second pressure differences are usable to determine a flow rate (such as the gas flow rate and/or the liquid flow rate) of the gas-liquid fluid mixture. Thus the first pressure meter of the method for providing measurements useable in determining a flow rate corresponds to the pressure meter of the more general method for providing measurements useable in determining the density of a fluid. The vertically-spaced measurement positions may be on a conduit transverse cross-section which also contains one of the horizontally-spaced measurement positions. In this case, that horizontally-spaced measurement position can coincide with one of the vertically-spaced measurement positions. 
         [0036]    For example, the apparatus may further include a conditioning element which conditions the fluid mixture to exhibit a predetermined type of flow at the measurement positions. The conditioning element may be a swirl element for inducing the mixture to exhibit swirling flow at the measurement positions. The swirl element may include one or more spiral-shaped members extending along the conduit in the direction of fluid flow. 
         [0037]    The second horizontally-spaced measurement position may be upstream of the constriction region. 
         [0038]    The constriction region may be provided by a Venturi. The constriction region may be the throat of the Venturi. The second horizontally-spaced measurement position may be at the inlet to the Venturi. 
         [0039]    The apparatus may further include a device for measuring the absolute pressure of the flow, e.g. at one of the measurement positions, but preferably on a transverse cross-section of the conduit containing the vertically-spaced measurement positions (conveniently this device can be incorporated in the first pressure meter itself). The apparatus may further include a device for measuring the temperature of the flow e.g. at one of the measurement positions, but preferably on the transverse cross-section of the conduit containing the vertically-spaced measurement positions. 
         [0040]    The apparatus may further include a third pressure meter for measuring a third pressure difference between two further vertically-spaced measurement positions in the conduit (e.g. at facing upper and lower wall portions of the conduit) spaced upstream or downstream of the vertically-spaced measurement positions where the first pressure difference is measured. The further vertically-spaced measurement positions can be on a conduit transverse cross-section which also contains one of the horizontally-spaced measurement positions. In particular, when the vertically-spaced measurement positions for the first pressure difference are on a conduit transverse cross-section which contains one of the horizontally-spaced measurement positions, the further vertically-spaced measurement positions for the third pressure difference can be on a conduit transverse cross-section which contains the other of the horizontally-spaced measurement positions. In this case, that other horizontally-spaced measurement position can coincide with one of the further vertically-spaced measurement positions. 
         [0041]    The pressure meter(s) may measure time-averaged pressure differences. 
         [0042]    The apparatus may further include a device for measuring the absolute pressure of the flow on a transverse cross-section of the conduit containing the further vertically-spaced measurement positions and preferably at one of the further vertically-spaced measurement positions (conveniently this device can be incorporated in the third pressure meter itself). The apparatus may further include a device for measuring the temperature of the flow on the transverse cross-section of the conduit containing the further vertically-spaced measurement positions and preferably at one of the further vertically-spaced measurement positions. 
         [0043]    A third aspect of the present invention provides a meter (e.g. a densitometer or a hold up meter) including the apparatus for measuring flow properties of a multiphase mixture and a processor arranged to determine a density of the fluid using the pressure difference measured by the pressure meter of the apparatus. 
         [0044]    The processor may calculate respective time-averaged pressure difference(s) from sequences of pressure differences measured by the pressure meter(s), the time-averaged pressure difference(s) then being used to determine the density or the flow rate. 
         [0045]    A further aspect of the present invention provides an oil well pipeline or a gas well pipeline including an apparatus according to the second aspect or a meter according to the third aspect. 
         [0046]    A further aspect of the present invention provides an apparatus according to the second aspect when conveying a gas-liquid fluid mixture, or a meter according to the third aspect when conveying a gas-liquid fluid mixture. 
         [0047]    Another aspect of the present invention provides a section of conduit having a conduit wall and having upper and lower fluid passages which each traverse the wall; 
         [0048]    wherein each fluid passage extends between an inner opening on the internal face of the wall and an outer opening on the external face of the wall, and the passages are arranged such that, when the section of conduit is viewed along its axial direction with the inner opening of the upper passage directly above the inner opening of the lower passage, the passages extend substantially horizontally between their respective inner and outer openings. 
         [0049]    Thus the section of conduit may be part of the conduit of the first or second aspect, the inner openings of the fluid passages providing the vertically-spaced measurement positions at which the first pressure difference is measured, or the further vertically-spaced positions at which the third pressure difference is measured. 
         [0050]    In use, the section of conduit will be arranged with the inner opening of the upper passage vertically above the inner opening of the lower passage. Because the passages will then extend substantially horizontally across the wall of the conduit section between their respective inner and outer openings, unnecessary fluid columns above the inner opening of the upper passage and below the inner opening of the lower passage can be avoided. This in turn can improve the sensitivity of a pressure meter connected to the outer openings of the fluid passages to changes in the gravitational pressure head between the inner openings. 
         [0051]    Typically, the section of conduit has a circular transverse cross-section. Preferably, the inner openings of the fluid passages are on the same transverse cross-section of the section of conduit. Preferably, the inner openings are at facing portions of the wall (i.e. on a section of conduit which has a circular transverse cross-section, the inner openings are preferably at diametrically opposing positions). Typically, the passages extend horizontally in a direction which is perpendicular to the axial direction of the section of conduit. 
         [0052]    The conduit section may further have a pressure meter connected across the outer openings of the fluid passages for measuring a pressure difference between the inner openings. 
         [0053]    Although in use the passages extend substantially horizontally across the wall of the conduit, each passage may extend initially from its inner opening in a direction which is perpendicular to that part of the internal face of the wall where the inner opening is located. Typically, this results in a short section of passage which is non-horizontal, or even vertical. When the section of conduit conveys a swirling flow of gas-liquid mixture, the liquid at the wall will have a substantial circumferential velocity component. However, by having the passage extend initially from the inner opening in a direction which is perpendicular to the internal face of the wall, this velocity component is prevented from forcing liquid from the swirling flow directly into the passages. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS  
         [0054]    Embodiments of the present invention will be described by way of example with reference to the accompanying drawings, in which: 
           [0055]      FIG. 1  shows schematically a longitudinal section through a first embodiment of an apparatus for providing measurements useable in determining a flow rate of a gas-liquid fluid mixture; 
           [0056]      FIG. 2  shows schematically the gas and liquid distribution for the section of  FIG. 1 ; 
           [0057]      FIG. 3  shows schematically a longitudinal section through a second embodiment of an apparatus for providing measurements useable in determining a flow rate of a gas-liquid fluid mixture; 
           [0058]      FIG. 4  is a plot of gas volume fraction against gas hold up; 
           [0059]      FIG. 5  is a schematic transverse cross-section of the conduit of  FIG. 1  or  3  at the measurement positions for ΔP 1 ; and 
           [0060]      FIG. 6  shows a schematic transverse cross-section through the apparatus of either  FIG. 1  or  3 . 
       
    
    
     DETAILED DESCRIPTION 
       [0061]      FIG. 1  shows schematically a longitudinal section through a first embodiment of an apparatus for providing measurements useable in determining a flow rate of a gas-liquid fluid mixture. 
         [0062]    The apparatus comprises a substantially horizontal conduit  1  of circular cross-section. The conduit has a Venturi  2 , the throat  3  of the Venturi providing a constriction region in the conduit. 
         [0063]    A gas-liquid fluid mixture flows through the conduit in the direction indicated by arrow  4 . A swirl element (not shown) induces the mixture to exhibit swirling flow as indicated by arrow  5 . An effect of this swirling flow is that liquid from the mixture is displaced to the wall of the conduit to form a liquid annulus around a gas core, shown schematically in  FIG. 2 . In the throat of the Venturi, centrifugal acceleration enhances the displacement of liquid to the wall of the conduit. 
         [0064]    Returning to  FIG. 1 , on a transverse cross-section at the Venturi throat, a first pressure meter  7  measures the pressure difference, ΔP 1 , between facing measurement positions on upper and lower wall portions of the throat. A second pressure meter  6  measures the pressure difference, ΔP 2 , between a measurement position in a throat of the Venturi at the transverse cross-section for ΔP 1 , and a measurement position in the conduit at the inlet to the Venturi. The diameter of the conduit at the Venturi inlet cross-section containing the upstream measurement position for ΔP 2  is D, and the diameter of the Venturi throat is d. 
         [0065]    Assuming there is no slip between the liquid and gas phases, the following analysis allows the gas and the liquid flow rates to be calculated. 
         [0000]    Notation: q=volumetric flow rate
       α=hold up   GVF=gas volume fraction   ρ=density   P=pressure   ΔP=pressure difference   v=velocity   Subscripts   T=total   L=liquid   G=gas   M=mixture   H=homogeneous   Superscripts   t=throat   i=inlet       
 
         [0081]    The total volumetric flow rate in the conduit is given by the Venturi equation: 
         [0000]    
       
         
           
             
               q 
               T 
             
             = 
             
               k 
                
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       P 
                       2 
                     
                   
                   
                     ρ 
                     M 
                   
                 
               
             
           
         
       
     
         [0000]    where k is a constant. 
         [0082]    Experience with the Vx™ system indicates that using the mixture density at the throat of the Venturi, ρ M   t , in the above equation provides more accurate results. This density can be obtained from the pressure difference across the throat, ΔP 1 =f(ρ M   t ,d). 
         [0083]    For example, it can be shown (see the Annex) that for swirling flow through a conduit of circular cross-section that: 
         [0000]      ΔP 1 ρ M   t dg 
         [0000]    where g is acceleration due to gravity. 
         [0084]    Thus the Venturi equation can be rewritten as: 
         [0000]    
       
         
           
             
               q 
               T 
               t 
             
             = 
             
               k 
                
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     
                       P 
                       2 
                     
                   
                   
                     
                       ρ 
                       M 
                       t 
                     
                      
                     
                       f 
                        
                       
                         ( 
                         
                           
                             Δ 
                              
                             
                                 
                             
                              
                             
                               P 
                               1 
                             
                           
                           , 
                           d 
                         
                         ) 
                       
                     
                   
                 
               
             
           
         
       
     
         [0085]    The liquid and gas densities, ρ L  and ρ G , are generally known at line conditions or can be estimated. For example, at typical line conditions of a mixture of natural gas and condensate, ρ L  is effectively a constant, and ρ G  varies in a known way (i.e. according to the ideal gas law) with pressure and temperature. If the pressure and/or temperature in the conduit are not known, they can be measured by suitable devices. For example, pressure meter  7  may incorporate a device for measuring the absolute pressure in the conduit at the throat, P t . A further device (not shown) may be provided for measuring the temperature at that position. 
         [0086]    Under the no slip assumption, the liquid and gas densities allow the gas hold up, α G , to be determined from the expression: 
         [0000]      ρ M =α G ρ G +(1−α G )ρ L    
         [0000]    which can be rearranged as: 
         [0000]    
       
         
           
             
               α 
               G 
             
             = 
             
               
                 
                   ρ 
                   L 
                 
                 - 
                 
                   ρ 
                   M 
                 
               
               
                 
                   ρ 
                   L 
                 
                 - 
                 
                   ρ 
                   G 
                 
               
             
           
         
       
     
         [0000]    whereby: 
         [0000]      q G   t =q T   t α G   t    
         [0000]    
       
      
       q 
       L 
       t 
       =q 
       T 
       t 
       −q 
       G 
       t  
      
     
         [0087]    Thus the volumetric flow rates of the gas and liquid in the throat can be calculated from knowledge of the liquid and gas densities in the throat, and from the measurements ΔP 1  and ΔP 2 . Evidently, these volumetric flow rates can be converted into mass flow rates, or into volumetric flow rates at other positions in the conduit (e.g. applying the ideal gas law, but ignoring the effect of temperature which will be small, q G   i  can be calculated from the expression 
         [0000]        q   G   i   =q   G   t ( P   t /(P t   +ΔP   2 ))). 
         [0088]    Similar analyses can be performed for types of flow other than swirling flow, but in such analyses the constant k in the Venturi equation and the expression relating ΔP 1  to the mixture density may be different. 
         [0089]      FIG. 3  shows schematically a longitudinal section through a second embodiment of an apparatus for providing measurements useable in determining a flow rate of a gas-liquid fluid mixture. The apparatus is similar to the apparatus of the first embodiment, and the same numbers are used in  FIGS. 1 and 3  for common features of the two embodiments. However, the second embodiment differs from the first embodiment in that on the inlet transverse cross-section containing the upstream measurement position for ΔP 2 , a third pressure meter  8  measures the pressure difference, ΔP 3 , between facing positions on upper and lower wall portions of the conduit. 
         [0090]    The following analysis allows the gas and the liquid flow rates to be calculated, with the additional pressure difference, ΔP 3 , enabling slip between the gas and liquid phases to be accounted for. As in the previous analysis, it assumes that ρ L  is a known constant, and ρ G  varies in a known way with pressure and temperature. 
         [0091]    The pressure differences ΔP 1  and ΔP 3  are related to the densities of the fluid mixture at respectively the inlet to the Venturi and the throat, i.e. ΔP 3 =f(ρ M   i ,D)and ΔP 1 =f(ρ M   t ,d). 
         [0092]    Thus, as in the previous analysis, the Venturi equation can be rewritten as: 
         [0000]    
       
         
           
             
               
                 
                   
                     q 
                     T 
                     t 
                   
                   = 
                   
                     k 
                      
                     
                       
                         
                           Δ 
                            
                           
                               
                           
                            
                           
                             P 
                             2 
                           
                         
                         
                           
                             ρ 
                             M 
                             t 
                           
                            
                           
                             f 
                              
                             
                               ( 
                               
                                 
                                   Δ 
                                    
                                   
                                       
                                   
                                    
                                   
                                     P 
                                     1 
                                   
                                 
                                 , 
                                 d 
                               
                               ) 
                             
                           
                         
                       
                     
                   
                 
               
               
                 
                   Equation 
                    
                   
                       
                   
                    
                   
                     ( 
                     1 
                     ) 
                   
                 
               
             
           
         
       
     
         [0093]    Further, the respective gas hold ups at the inlet to the Venturi and the throat can be determined from: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       α 
                       G 
                       i 
                     
                     = 
                     
                       
                         ( 
                         
                           
                             ρ 
                             L 
                           
                           - 
                           
                             ρ 
                             M 
                             i 
                           
                         
                         ) 
                       
                       
                         ( 
                         
                           
                             ρ 
                             L 
                           
                           - 
                           
                             ρ 
                             G 
                             i 
                           
                         
                         ) 
                       
                     
                   
                    
                   
                     
 
                   
                    
                   
                     
                       α 
                       G 
                       t 
                     
                     = 
                     
                       
                         ( 
                         
                           
                             ρ 
                             L 
                           
                           - 
                           
                             ρ 
                             M 
                             t 
                           
                         
                         ) 
                       
                       
                         ( 
                         
                           
                             ρ 
                             L 
                           
                           - 
                           
                             ρ 
                             G 
                             t 
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   Equation 
                    
                   
                       
                   
                    
                   
                     ( 
                     2 
                     ) 
                   
                 
               
             
           
         
       
     
         [0000]    whereby the change in gas hold up from the inlet to the throat of the Venturi can be calculated: 
         [0000]      Δα G =α G   t −α G   i    Equation (3) 
         [0094]    Next, we seek an expression for the change in GVF from the inlet to the throat of the Venturi. Firstly, applying the ideal gas law: 
         [0000]    
       
         
           
             
               q 
               G 
               t 
             
             = 
             
               
                 q 
                 G 
                 i 
               
                
               
                 
                   P 
                   i 
                 
                 
                   P 
                   t 
                 
               
                
               
                 
                   T 
                   t 
                 
                 
                   T 
                   i 
                 
               
             
           
         
       
     
         [0095]    Because typical variations in absolute temperature are small, the T t /T i  term will be close to unity and can be disregarded. Conveniently, pressure meters  7 ,  8  incorporate devices to measure P i  and P t  as well as ΔP 1  and ΔP 3 . Thus the change in gas flow rate from the inlet to the throat of the Venturi can be expressed as: 
         [0000]    
       
         
           
             
               Δ 
                
               
                   
               
                
               
                 q 
                 G 
               
             
             = 
             
               
                 
                   q 
                   G 
                   t 
                 
                 - 
                 
                   q 
                   G 
                   i 
                 
               
               = 
               
                 
                   q 
                   G 
                   t 
                 
                  
                 
                   ( 
                   
                     1 
                     - 
                     
                       
                         P 
                         t 
                       
                       
                         P 
                         i 
                       
                     
                   
                   ) 
                 
               
             
           
         
       
     
         [0096]    Now, GVF=q G /(q G +q L ), whereby, on the reasonable assumption that q L  is invariant: 
         [0000]    
       
         
           
             
               
                 Δ 
                  
                 
                     
                 
                  
                 G 
                  
                 
                     
                 
                  
                 V 
                  
                 
                     
                 
                  
                 F 
               
               
                 Δ 
                  
                 
                     
                 
                  
                 
                   q 
                   G 
                 
               
             
             = 
             
               
                 
                   q 
                   L 
                 
                 
                   
                     ( 
                     
                       
                         q 
                         G 
                       
                       + 
                       
                         q 
                         L 
                       
                     
                     ) 
                   
                   2 
                 
               
               = 
               
                 
                   q 
                   L 
                 
                 
                   q 
                   T 
                   2 
                 
               
             
           
         
       
     
         [0097]    Thus: 
         [0000]    
       
         
           
             
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     G 
                      
                     
                         
                     
                      
                     V 
                      
                     
                         
                     
                      
                     F 
                   
                   = 
                   
                     
                       
                         q 
                         L 
                       
                       
                         q 
                         T 
                         
                           t 
                           2 
                         
                       
                     
                     · 
                     
                       
                         q 
                         G 
                         t 
                       
                        
                       
                         ( 
                         
                           1 
                           - 
                           
                             
                               P 
                               t 
                             
                             
                               P 
                               i 
                             
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   Equation 
                    
                   
                       
                   
                    
                   
                     ( 
                     4 
                     ) 
                   
                 
               
             
           
         
       
     
         [0098]    Slip can be expressed in the relation between GVF and α G . That is: 
         [0000]    
       
         
           
             
               G 
                
               
                   
               
                
               V 
                
               
                   
               
                
               F 
             
             = 
             
               
                 
                   q 
                   G 
                 
                 
                   q 
                   T 
                 
               
               = 
               
                 
                   α 
                   G 
                 
                  
                 
                   
                     V 
                     G 
                   
                   
                     V 
                     H 
                   
                 
               
             
           
         
       
     
         [0099]    In general, where there is slip, GVF≧α G , although when there is single phase flow GVF=α G =1 and GVF=α G =0. Thus, as shown in  FIG. 4 , the slip law can be approximated by two straight lines A and B representing a deviation from a situation where there is no slip (i.e. GVF=α G ). 
         [0100]    For line A: 
         [0000]    
       
         
           
             
               
                 
                   Δ 
                    
                   
                       
                   
                    
                   G 
                    
                   
                       
                   
                    
                   V 
                    
                   
                       
                   
                    
                   F 
                 
                 
                   Δ 
                    
                   
                       
                   
                    
                   
                     α 
                     G 
                   
                 
               
               ≥ 
               1 
             
             , 
             
               
                 and 
                  
                 
                     
                 
                  
                 G 
                  
                 
                     
                 
                  
                 V 
                  
                 
                     
                 
                  
                 F 
               
               = 
               
                 
                   0 
                    
                   
                       
                   
                    
                   when 
                    
                   
                       
                   
                    
                   
                     α 
                     G 
                   
                 
                 = 
                 
                   
                     0 
                      
                     
                       
 
                     
                     ⇒ 
                     
                       G 
                        
                       
                           
                       
                        
                       V 
                        
                       
                           
                       
                        
                       F 
                     
                   
                   = 
                   
                     
                       
                         Δ 
                          
                         
                             
                         
                          
                         G 
                          
                         
                             
                         
                          
                         V 
                          
                         
                             
                         
                          
                         F 
                       
                       
                         Δα 
                         G 
                       
                     
                      
                     
                       α 
                       G 
                     
                   
                 
               
             
           
         
       
     
         [0000]    while for line B: 
         [0000]    
       
         
           
             
               
                 
                   Δ 
                    
                   
                       
                   
                    
                   G 
                    
                   
                       
                   
                    
                   V 
                    
                   
                       
                   
                    
                   F 
                 
                 
                   Δ 
                    
                   
                       
                   
                    
                   
                     α 
                     G 
                   
                 
               
               &lt; 
               1 
             
             , 
             
               
                 and 
                  
                 
                     
                 
                  
                 G 
                  
                 
                     
                 
                  
                 V 
                  
                 
                     
                 
                  
                 F 
               
               = 
               
                 
                   1 
                    
                   
                       
                   
                    
                   when 
                    
                   
                       
                   
                    
                   
                     α 
                     G 
                   
                 
                 = 
                 
                   
                     1 
                      
                     
                       
 
                     
                     ⇒ 
                     
                       G 
                        
                       
                           
                       
                        
                       V 
                        
                       
                           
                       
                        
                       F 
                     
                   
                   = 
                   
                     
                       
                         
                           Δ 
                            
                           
                               
                           
                            
                           G 
                            
                           
                               
                           
                            
                           V 
                            
                           
                               
                           
                            
                           F 
                         
                         
                           Δα 
                           G 
                         
                       
                        
                       
                         ( 
                         
                           
                             α 
                             G 
                           
                           - 
                           1 
                         
                         ) 
                       
                     
                     + 
                     1 
                   
                 
               
             
           
         
       
     
         [0101]    Applying this slip law, it is then possible to iteratively determine the GVF and hence to calculate the gas and liquid flow rates. The steps of the iteration are as follows:
       1) Determine q T   t  from Equation (1)   2) Determine α G   t  from Equation (2)   3) Determine Δα G  from Equation (3)   4) Assuming no slip, set GVF t =α G   t      5) q G   t =GVF t ·q T   t q L =q T   t −q G   t      6) Determine ΔGVF from Equation (4)   7) Calculate       
 
         [0000]    
       
         
           
             
               Δ 
                
               
                   
               
                
               G 
                
               
                   
               
                
               V 
                
               
                   
               
                
               F 
             
             
               Δ 
                
               
                   
               
                
               
                 α 
                 G 
               
             
           
         
       
       
         
           
             8) If 
           
         
       
     
         [0000]    
       
         
           
             
               
                 Δ 
                  
                 
                     
                 
                  
                 G 
                  
                 
                     
                 
                  
                 V 
                  
                 
                     
                 
                  
                 F 
               
               
                 Δα 
                 G 
               
             
             ≥ 
             1 
           
         
       
     
         [0000]    then reset 
         [0000]    
       
         
           
             
               G 
                
               
                   
               
                
               V 
                
               
                   
               
                
               
                 F 
                 t 
               
             
             = 
             
               
                 
                   Δ 
                    
                   
                       
                   
                    
                   G 
                    
                   
                       
                   
                    
                   V 
                    
                   
                       
                   
                    
                   F 
                 
                 
                   Δα 
                   G 
                 
               
                
               
                 α 
                 G 
                 t 
               
             
           
         
       
     
         [0000]    (slip law line A), whereas if 
         [0000]    
       
         
           
             
               
                 Δ 
                  
                 
                     
                 
                  
                 G 
                  
                 
                     
                 
                  
                 V 
                  
                 
                     
                 
                  
                 F 
               
               
                 Δα 
                 G 
               
             
             &lt; 
             1 
           
         
       
     
         [0000]    then reset 
         [0000]    
       
         
           
             
               G 
                
               
                   
               
                
               V 
                
               
                   
               
                
               
                 F 
                 t 
               
             
             = 
             
               
                 
                   
                     Δ 
                      
                     
                         
                     
                      
                     G 
                      
                     
                         
                     
                      
                     V 
                      
                     
                         
                     
                      
                     F 
                   
                   
                     Δα 
                     G 
                   
                 
                  
                 
                   ( 
                   
                     
                       α 
                       G 
                       t 
                     
                     - 
                     1 
                   
                   ) 
                 
               
               + 
               1 
             
           
         
       
     
         [0000]    (slip law line B)
       9) Iterate around 5) to 8) until GVF t  has converged   10)q G   t =GVF t ·q T   t q L =q T   t −q G   t          
 
         [0112]    Again, these volumetric flow rates can be converted into mass flow rates, or into volumetric flow rates at other positions in the conduit. 
         [0113]    The analyses described above in relation to both the first and second embodiments can be performed by a suitably arranged processor (not shown in  FIGS. 1 and 3 ) which receives the pressure difference measurements from the pressure meters of the respective apparatus (and optionally receives absolute pressure and temperature measurements). The flow rate calculations resulting from the analyses can be then be transmitted, stored and/or displayed. Because the swirling flow may not be steady, and to reduce the effect of noise in the pressure difference measurements, it can be advantageous for the analyses to be performed on time-averaged pressure difference measurements. 
         [0114]    Although the pressure differences ΔP 1  and ΔP 3  will generally be small, it is possible to measure the differences with enough accuracy. For example, putting approximate values of d=0.1 m and g=10 m s −2  into the expression ΔP 1 =ρ M   t dg, and assuming a maximum value for ρ M   t =1000 kg m −3 , gives an approximate upper value for ΔP 1 =1000 Pa. Honeywell&#39;s ST 3000™ Series 100 differential pressure meter model STD120 is an example of a conventional pressure meter that could be used to measure ΔP 1  and ΔP 3 . This meter should have a sufficient accuracy of about ±1 Pa over the span 0 to 1000 Pa. 
         [0115]    Indeed, because the ΔP 1  term in the Venturi equation (Equation (1) above), appears in a square root, the impact of errors in the measurement of ΔP 1  on the calculated flow rates is reduced. 
         [0116]    Precautions can also be taken to improve the accuracy of the ΔP 1  and ΔP 3  measurements. For example, as shown in  FIG. 5 , which is a schematic transverse cross-section of the conduit of  FIG. 1  or  3  at the measurement positions for ΔP 1 , fluid-filled passages  9 ,  10  connecting to pressure meter  7  extend substantially horizontally through the wall  11  of the conduit. By avoiding unnecessary fluid columns above and below the locations where the passages open into the conduit at facing upper and lower wall portions, the sensitivity of the pressure meter to changes in the gravitational pressure head between the passages openings is improved. Nonetheless, the passages do have short vertical sections  12 ,  13  adjacent the openings into the conduit. These short sections help to prevent liquid from the swirling flow  5  being forced into the passages by the circumferential velocity component of the flow. Preferably, pressure meter  7  is isolated from the conduit fluid by isolation bellows  14 ,  15  in passages  9 ,  10 . Suitable bellows are available from e.g. Honeywell in their ST 3000™ Series 100. 
         [0117]    While the invention has been described in conjunction with the exemplary embodiments described above, many equivalent modifications and variations will be apparent to those skilled in the art when given this disclosure. Accordingly, the exemplary embodiments of the invention set forth above are considered to be illustrative and not limiting. Various changes to the described embodiments may be made without departing from the scope of the invention.