# Patent

## Title
Multiphase flowmeter using a combination of pressure differentials and ultrasound doppler readings

## Abstract

A multiphase flowmeter, and a method of analyzing and measuring multiphase flow is described. The body of the multiphase flowmeter consists of the following components. It has a straight flow pipe in which two orifice plates a lodged near the extremities. At each of these two locations, a pressure differential transmitter measures the pressure drop that occurs as the fluid passes through the orifice plates. Along the length between the orifice plates, an ultrasound Doppler sensor measures the volume flow rate of the multiphase fluid, a pressure transmitter measures the in-pipe pressure and a temperature sensor measures the temperature. From the signals obtained from the above sensors, an overall analysis of the multiphase fluid flow is performed providing a complete set of flow measurement data for the mixture and for each of the three phases consisting of oil, water and gas.

## Background
CROSS REFERENCE TO RELATED APPLICATION
This application is a National Stage of International Application No. PCT/CN2007/003068, filed Oct. 29, 2007. This application claims the benefit and priority of Chinese Application No. 200710044625.8, filed Aug. 7, 2007. The entire disclosure of each of the above applications is incorporated herein by reference.
FIELD OF INVENTION
The present invention relates to an apparatus and method of analysis for the measurement of a multiphase fluid flow. The invention is particularly, but not exclusively, suitable for measuring the flow of multiphase fluids that are produced in oil and gas wells, on shore and off shore installations, pipe lines and in refineries.
BACKGROUND OF THE INVENTION
The oil and gas industry is increasingly calling for the development of a new technology that is compact, light weight and most importantly affordable, which can be installed “in-line” and produce a complete set of accurate measurements on each component of multiphase flow.
Various devices are used in the industry, however none is considered capable of meeting satisfactorily all the requirements that are posed by wide ranging set of field conditions in the industry. No single device handles the many flow and fluid conditions that are encountered due to both geological conditions and industrial methods of extraction, and/or meets the precision requirements over the full extent of the large range of flow rates, water cut, and gas fraction that occur in the field. As a result, the devices that exist tend to find niche applications only, or they are combined with other equipments and devices in order to fully service the flow measurements of multiphase fluids. In addition, several devices make use of radioactive sources, which hold significant disadvantages. The use of these sources imposes careful and important containment requirements to mitigate the possibility of contamination, and as a result the industry is reluctant to fully accept methods making use of radioisotopes.
Commercialized MPFMs (multiphase flowmeters) can roughly be divided into two categories:
The first category is based on a pre-requisite step of separation of the liquid and gas phases. Once the liquid and gas phases are separated, the flow measurements are conducted upon the liquid phase and the gas phase separately. The principles of such a method are simple and well known at large and good precision in measurement is generally achieved. This category of MPFM has been widely accepted in the oil industry worldwide. One example is U.S. Pat. No. 6,338,276 B1. The separation of liquid and gas phases are usually achieved using gravity or centrifugal forces. The equipment, called separators, is generally large, difficult to install and relatively costly. The degree of precision/imprecision of the measurements is directly affected by the efficiency of the separation process and compounded by the inherent precision/imprecision of the individual and separate measurements of the liquid and gas phases that are performed after the separation. Room for improvement of the achievable precision inherent to this method is therefore limited. In addition measurements are neither performed in-line or on a real-time.
The second category of MPFMs does not require any fluid pre-separation. It measures directly the various parameters of the multiphase flow. It generally uses an orifice or venturi flowmeter for flow rate measurement. It also uses Gamma radioisotopes, microwaves or capacitance/impedance watercut meters to independently measure watercut and GOR. Generally this method can perform inline real-time measurements of multiphase flow with an acceptable level of precision/imprecision, and it has achieved a degree of acceptance by the industry. An example of this category of device is provided by the Framo Phase Watch VX and U.S. Pat. No. 6,935,189 B2. However, because the method used by such devices is sensitive to interferences in flow regimes and in patterns of multiphase flow, the stability of the measurements represent an issue that is not fully resolved by the technology. In addition, in many areas of the world the water fraction contained in the fluid is disturbed with varying degrees of salinity, affecting the density of the water and the precision of measurements that use impedance/capacitance techniques. The structure of the MPFMs in this category is generally complex and costly. The calibration, installation and maintenance all poses significant operational challenges to the industry.
The main objective of the present invention is to meet the many requirements of the industry with a device that is completely self-sufficient, flow regime independent, compact, light-weight, easy to install and accurate in measuring multiphase flow.
SUMMARY OF THE INVENTION
By performing a complete analysis of the multiphase fluid on the basis of the pressure and ultrasound sensors within the device, this invention solves the problems associated with radio-active sources and with the need to make use of combined technologies.
The technology presented here is therefore different with other existing and published technologies and methods. The following aspects, either individually or combined, constitute this differentiation:
The present invention does not require any pre-conditioning of the multiphase flow, such as pre-separation of the gas and fluid and/or mixing.
The present invention does not make use of upstream or downstream devices to measure water fraction or perform separately other two-phase in well measurements such as gas-oil ratio. In other words it is self sufficient in the sense that it analyses and measures multiphase flow completely from direct readings performed by the device's sensors.
The present invention does not make use of any radioactive sources.
The present invention is simple, accurate with measurement stability and repeatability.
First the pressure differentials are obtained from the transmitters located at the two orifice plates, a volume flow rate of the mixture is obtained from the ultrasound Doppler sensor. From these data, the microcomputer calculates the gas-liquid ratio, the watercut and the flow rates of the gas, water and oil phases. The in-pipe gas flow rate is then converted to a gas flow rate under standard atmospheric pressure.
The invention provides the measurements of volume and mass flow rates and respective cumulative for each individual phase of oil, water and gas, and for the total mixture, as well as the watercut and the gas-liquid ratio.



BRIEF DESCRIPTION OF THE DRAWINGS
These and other aspects of the invention will become apparent from the following description, when taken in combination with accompanying drawings, in which:
FIG. 1 is a diagrammatic view of the multiphase mass flowmeter;
FIG. 2 illustrates the design principal of the multiphase mass flowmeter, and
FIG. 3 is the working flowchart of the multiphase mass flowmeter.



SYMBOLS USED IN THE DIAGRAMS


    
    
        10 is the uniform pipe
        1 is the upstream orifice plate
        2 is the downstream orifice plate
        3 is the differential pressure transmitter at the upstream orifice plate
        4 is the differential pressure transmitter ate the downstream orifice plate
        5 is the ultrasound Doppler sensor
        6 is the pressure transmitter
        7 is the temperature transmitter
        8 is a microcomputer
        9 is the display module
        D is the inner diameter of the uniform pipe
        d1 is the diameter of orifice plate 1
        d2 is the diameter of orifice plate 2
        X is the direction of flow
    
    


DETAILED DESCRIPTION OF THE INVENTION
The following parameters are defined for use in the equations that follow in this description:
cd1=flow rate coefficients of orifice plate 1
cd2=flow rate coefficients of orifice plate 2
D=diameter of uniform pipe
d1=inner diameter of orifice plate 1
d2=inner diameter of orifice plate 2
dp1=pressure differential at orifice plate 1
dp2=pressure differential at orifice plate 2
n1=gas liquid ratio (GLR)
QS=flow rate obtained with the ultrasound Doppler
QL1=flow rate of liquid
Qm1=mass flow rate at orifice plate 1
Qm2=mass flow rate at orifice plate 2
Q1=volume flow rate at orifice plate 1
Q2=volume flow rate at orifice plate 2
Qg1=gas phase volume flow rate
Qw1=water phase volume flow rate
Qo1=oil phase volume flow rate
ρ1=density of total mixture at orifice plate 1
ρ2=density of total mixture at orifice plate 2
WR1=water cut
Reference is made to FIG. 1 of the drawings which depicts a uniform flow pipe 10 of diameter D. Along the direction of the flow, as indicated by the arrow X, two orifice plates 1 and 2 are placed at defined locations. The upstream orifice plate 1 and the downstream orifice plate 2 have inner diameters d1 and d2 respectively. At a defined location between the orifice plates 1 and 2, an ultrasound Doppler 5 is installed. At arbitrary locations between the two orifice plates, a pressure transmitter 6 and a temperature transmitter 7 are installed.
When multiphase flow passes through the uniform pipe 10, along the direction X, a pressure drop takes place at orifice plate 1 due to the localized constraining of the flow area. This pressure drop dp1 is measured by the differential pressure transmitter 3. Similarly, at orifice plate 2, the pressure drop dp2 takes place and is measured by the differential pressure transmitter 4. In between these two points, the ultrasound Doppler flowmeter 5 measures the flow rate Qs, the pressure transmitter 6 measures the pressure P and the temperature transmitter 7 measures the temperature T. This invention and its method for analyzing multiphase flow and for measuring the individual flow rates for the phases (oil, water and gas) makes use of the above parameters, which are measured directly.
The mass flow rates of the total mixture when passing through orifice plates 1 and 2 are obtained independently as follows:




  
    
      
        
          
            Q
            
              m
              ⁢
              
                
              
              ⁢
              1
            
          
          =
          
            
              
                
                  c
                  1
                
                ⁢
                
                  c
                  
                    d
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
                ⁢
                
                  d
                  1
                  2
                
              
              
                
                  1
                  -
                  
                    
                      (
                      
                        
                          d
                          1
                        
                        D
                      
                      )
                    
                    4
                  
                
              
            
            ⁢
            
              
                
                  ρ
                  1
                
                ⁢
                
                  dp
                  1
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            Q
            
              m
              ⁢
              
                
              
              ⁢
              2
            
          
          =
          
            
              
                
                  c
                  2
                
                ⁢
                
                  c
                  
                    d
                    ⁢
                    
                      
                    
                    ⁢
                    2
                  
                
                ⁢
                
                  d
                  2
                  2
                
              
              
                
                  1
                  -
                  
                    
                      (
                      
                        
                          d
                          2
                        
                        D
                      
                      )
                    
                    4
                  
                
              
            
            ⁢
            
              
                
                  ρ
                  2
                
                ⁢
                
                  dp
                  2
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        where
        ⁢
        
          

        
        ⁢
        
          
            c
            1
          
          =
          
            
              c
              2
            
            =
            
              
                π
                4
              
              ⁢
              
                2
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          and
          ⁢
          
            
          
          ⁢
          let
        
      
    
    
      
        (
        1
        )
      
    
  
  
    
      
        
          
            
              k
              1
            
            =
            
              
                d
                1
                4
              
              
                1
                -
                
                  
                    (
                    
                      
                        d
                        1
                      
                      D
                    
                    )
                  
                  4
                
              
            
          
          ,
          
            
          
          ⁢
          
            
              k
              2
            
            =
            
              
                d
                2
                4
              
              
                1
                -
                
                  
                    (
                    
                      
                        d
                        2
                      
                      D
                    
                    )
                  
                  4
                
              
            
          
        
        ⁢
        
          

        
      
    
    
      
        (
        2
        )
      
    
  




From the mass flow rate conservation, then:

Qm1=Qm2 

Applying equation (1) to the above, one gets




  
    
      
        
          G
          =
          
            
              
                
                  k
                  2
                
                ⁢
                
                  dp
                  2
                
              
              
                
                  k
                  1
                
                ⁢
                
                  dp
                  1
                
              
            
            =
            
              
                
                  ρ
                  1
                
                
                  ρ
                  2
                
              
              =
              
                
                  Q
                  2
                
                
                  Q
                  1
                
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          
            where
            ⁢
            
              
            
            ⁢
            
              Q
              1
            
          
          =
          
            
              
                
                  Q
                  
                    m
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
                
                  ρ
                  1
                
              
              ⁢
              
                
              
              ⁢
              and
              ⁢
              
                
              
              ⁢
              
                Q
                2
              
            
            =
            
              
                Q
                
                  m
                  ⁢
                  
                    
                  
                  ⁢
                  2
                
              
              
                ρ
                2
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




Here using subscripts o, w and g to represent oil, water and gas, then:

Q1=Qo1+Qw1+Qg1 

Q2=Qo2+Qw2+Qg2  (4)

From orifice plate 1 to 2, the volume flow rates of oil, water and gas in the multiphase flow will satisfy the following linear relationships,

Qo2=aQo2, Qw2=bQw1, Qg2=cQg1 

Where a, b and c are the coefficients of linearity for each phase relating the volume flow rate in orifice plate 1 to the volume flow rate in orifice plate 2. From Equation (3),




  
    
      
        G
        =
        
          
            
              
                k
                2
              
              ⁢
              
                dp
                2
              
            
            
              
                k
                1
              
              ⁢
              
                dp
                1
              
            
          
          =
          
            
              
                aQ
                
                  o
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              +
              
                bQ
                
                  w
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              +
              
                cQ
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
            
            
              
                Q
                
                  o
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              +
              
                Q
                
                  w
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              +
              
                Q
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




After the multiphase flow passes orifice plate 1, the watercut WR1 and the gas liquid ratio (GLR) defined as n1 are:




  
    
      
        
          
            WR
            1
          
          =
          
            
              Q
              
                w
                ⁢
                
                  
                
                ⁢
                1
              
            
            
              
                Q
                
                  o
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              +
              
                Q
                
                  w
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
            
          
        
        ,
        
          
        
        ⁢
        
          
            n
            1
          
          =
          
            
              Q
              
                g
                ⁢
                
                  
                
                ⁢
                1
              
            
            
              
                Q
                
                  o
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              +
              
                Q
                
                  w
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  




From Equation (5):




  
    
      
        G
        =
        
          
            
              
                k
                2
              
              ⁢
              
                dp
                2
              
            
            
              
                k
                1
              
              ⁢
              
                dp
                1
              
            
          
          =
          
            
              
                a
                ⁡
                
                  (
                  
                    1
                    -
                    
                      WR
                      1
                    
                  
                  )
                
              
              +
              
                bWR
                1
              
              +
              
                cn
                1
              
            
            
              1
              +
              
                n
                1
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




Obviously, the coefficients a, b and c can be expressed by the values of G at pure oil, pure water and pure gas phase states:

a=G|pure oil b=G|pure water c=G|pure gas 

On the other hand, the volume flow rate Qs obtained by the ultrasound Doppler sensor and the liquid phase flow rate QL1 have the following relationship:




  
    
      
        
          
            Q
            s
          
          =
          
            
              1
              s
            
            ⁢
            
              Q
              
                L
                ⁢
                
                  
                
                ⁢
                1
              
            
          
        
        ⁢
        
          

        
        ⁢
        
          Therefore
          ,
        
      
    
    
      
        (
        8
        )
      
    
  
  
    
      
        
          
            Q
            
              o
              ⁢
              
                
              
              ⁢
              1
            
          
          +
          
            W
            
              w
              ⁢
              
                
              
              ⁢
              1
            
          
          +
          
            Q
            
              g
              ⁢
              
                
              
              ⁢
              1
            
          
        
        =
        
          
            
              Q
              
                L
                ⁢
                
                  
                
                ⁢
                1
              
            
            +
            
              Q
              
                g
                ⁢
                
                  
                
                ⁢
                1
              
            
          
          =
          
            
              sQ
              s
            
            +
            
              Q
              
                g
                ⁢
                
                  
                
                ⁢
                1
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  
  
    
      
        
          
            
              Q
              1
            
            -
            
              Q
              
                L
                ⁢
                
                  
                
                ⁢
                1
              
            
          
          
            Q
            
              L
              ⁢
              
                
              
              ⁢
              1
            
          
        
        =
        
          
            
              Q
              
                g
                ⁢
                
                  
                
                ⁢
                1
              
            
            
              Q
              
                L
                ⁢
                
                  
                
                ⁢
                1
              
            
          
          =
          
            
              
                Q
                
                  g
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
              
              
                sQ
                s
              
            
            =
            
              n
              1
            
          
        
      
    
    
      
        (
        10
        )
      
    
  




Combining Equations (7) and (9), the GLR n1 and watercut WR1 can be expressed using the ultrasound Doppler flow rate Qs the coefficients a, b and c, and the various G values under pure phase states.
Especially under the condition of pure liquid flow (only oil and water phases combined), Qg1=0, or when GOR n1=0, from Equation (7), the watercut WR1 can be directly expressed using G as follows:




  
    
      
        G
        =
        
          
            
              
                k
                2
              
              ⁢
              
                dp
                2
              
            
            
              
                k
                1
              
              ⁢
              
                dp
                1
              
            
          
          =
          
            
              
                (
                
                  b
                  -
                  a
                
                )
              
              ⁢
              
                WR
                1
              
            
            +
            a
          
        
      
    
    
      
        (
        11
        )
      
    
  




Finally, after the liquid phase volume flow rate QL1 is determined by Equation (8), the volume flow rates of gas, water and oil phases can be expressed as:

Qg1=n1QL1 

Qw1=WR1·QL1 

Qn1=(1−WR1)QL1  (12)

FIG. 3 gives the flowchart of the multiphase mass flowmeter. The pressure differential readings dp1 and dp2 from the pressure differential transmitters 3 and 4, as well as the volume flow rate Qs from the ultrasound Doppler flowmeter 5, are inputted into the intelligent microcomputer 8, where the software calculates GOR a1 and watercut WR1 using Equations (10) and (7). Finally, the software will calculate the individual phase flow rates of gas, water and oil using Equation (12). In FIG. 3, the pressure P and temperature T from transmitters 6 and 7 enable the conversion of the above gas phase flow rate into a standard atmospheric value. As per the design principle illustrated in FIG. 3, the parameters computed by the invention are displayed on the display module 9, including:
Flow pressure
Flow temperature
Pressure differential across the device
Actual oil flow rate
Standard oil flow rate
Actual water flow rate
Actual water cut
Standard water cut
Actual gas flow rate
Standard gas flow rate
Mix density
Mix velocity
Actual gas volume fraction
Standard gas volume fraction
Standard gas oil ratio
Accumulated oil volume in actual condition
Accumulated oil volume in standard condition
Accumulated water volume in actual condition
Accumulated water volume in standard condition
Accumulated gas volume in actual condition
Accumulated gas volume in standard condition
It is evident that components of the body of the flowmeter could be disposed in a variety of configurations without departing from the scope of the invention. Although all flowmeters are shown as a compact structure with means of creating and measuring pressure differentials in a single conduit, it will be appreciated by persons skilled in the art that the components can be disposed in different orders or widely spaced with respect to each other.
The following are the principal advantages of the present invention. Radioactive sources are not used. The device is self sufficient in the sense that it analyses and measures multiphase flow without requiring other devices to be installed upstream or downstream in order to measure water cut or perform other two-phase in well measurements. It does not need other upstream devices to perform fluid separation or fluid conditioning. The system is compact and easily installed on any single conduit at any point immediately downstream of the well head.

## Claims

1. A multiphase flowmeter for measuring the flow of a multiphase fluid including oil, water and gas flowing through an in-well conduit, said flowmeter comprising:
two orifices plates having inner diameters respectively where disposed in said conduit at locations near the extremities of a uniform flow pipe where pressure differentials are measured with two pressure differential transducers;
an ultrasound Doppler sensor between said two orifices plates;
a pressure transmitter; and
a temperature transmitter.

2. A flowmeter as in claim 1 where the ultrasound Doppler sensor, pressure transmitter and temperature transmitter are disposed over a plurality of tubular members or in different order or locations.

3. A flowmeter as in claim 1 where the intelligent microcomputer calculates the individual phase flow rates of gas, water and oil with a pressure differential at one orifice plate (dp1) and a pressure differential at the other orifice plate (dp2) from the pressure differential transmitters, the volume flow rate (Qs) from the ultrasound Doppler sensor as well as the pressure (P) and temperature (T) from transmitters.

4. A flowmeter as in claim 3 where a display unit is connected to the intelligent microcomputer to receive and display the values of the individual phase flow rates of gas, water and oil that have been calculated.

5. A method of measuring the flow rate of components of the multiphase fluid including oil, water and gas flowing in an in-well conduit, said method comprising the steps of:
measuring the pressure differentials at the two orifice plates having inner diameters respectively;
measuring the mixture volume flow rate at with the ultrasound Doppler sensor between said two orifices plates; and
measuring in-pipe pressure and temperature at any point between the two orifice plates, and processing the two pressure differentials at the two orifice plates, the flow rate at the ultrasound Doppler sensor, the in-pipe pressure and temperature measurements to provide a volume and mass flow rate of the total mixture and of each component of said multiphase fluid comprising oil, water and gas, a watercut measurement and a gas liquid ratio measurement.

6. A method as in claim 5 where the intelligent microcomputer calculates the individual phase flow rates of gas, water and oil with a pressure differential at one orifice plate (dp1) and a pressure differential at the other orifice plate (dp2) from the pressure differential transmitters, the volume flow rate (Qs) from the ultrasound Doppler sensor as well as the pressure (P) and temperature (T) from transmitters.

7. A method as in claim 6 where a display unit is connected to the intelligent microcomputer to receive and display the values of the individual phase flow rates of gas, water and oil that have been calculated.