CA2436965C - Method and apparatus for determining component flow rates for a multiphase flow - Google Patents
Method and apparatus for determining component flow rates for a multiphase flow Download PDFInfo
- Publication number
- CA2436965C CA2436965C CA002436965A CA2436965A CA2436965C CA 2436965 C CA2436965 C CA 2436965C CA 002436965 A CA002436965 A CA 002436965A CA 2436965 A CA2436965 A CA 2436965A CA 2436965 C CA2436965 C CA 2436965C
- Authority
- CA
- Canada
- Prior art keywords
- values
- measured
- fluid
- flow rates
- conduit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/74—Devices for measuring flow of a fluid or flow of a fluent solid material in suspension in another fluid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/66—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
Abstract
A method for determining component flow rates of a multiphase fluid in a conduit (the fluid consisting of at least three known components) includes the steps of: measuring at each of two different positions along the conduit at least four mixture quantities, typically the sound speed, the flow velocity of the multiphase fluid, the pressure and the temperature; providing a speed of sound in each of the components at the measured pressures and temperatures;
providing a trial value for each of either the component flow rates or the phase fractions; using a predetermined model to calculate values for the measured mixture quantities based on the trial values for each of either the component flow rates or the phase fractions; using a predetermined error function to determine an error value; and using a predetermined optimizing algorithm to determine whether the calculated values are acceptable, and, if they are not, to provide a new trial value for each of either the component flow rates or the phase fractions. In some applications, the error function is the sum of the squares of the difference between the measured and calculated values at each point.
providing a trial value for each of either the component flow rates or the phase fractions; using a predetermined model to calculate values for the measured mixture quantities based on the trial values for each of either the component flow rates or the phase fractions; using a predetermined error function to determine an error value; and using a predetermined optimizing algorithm to determine whether the calculated values are acceptable, and, if they are not, to provide a new trial value for each of either the component flow rates or the phase fractions. In some applications, the error function is the sum of the squares of the difference between the measured and calculated values at each point.
Description
"Method and Apparatus for Determining Component Flow Rates for a Multiphase Flow"
The present invention concerns techniques for the measurement of the flow rate of each component of a multiphase fluid. More particularly, the present invention concerns techniques based on measurements at different locations along a conduit for determining the flow rates of each component of a multiphase fluid within the conduit.
In extracting oil and gas from a formation, it is advantageous to monitor the flow -rates of the different components of the production fluid, usually gas, oil and water.
Such a three-component mixture presents a problem that is often more difficult to solve than that of determining the flow rates for binary mixtures. For a broad class of well mixed binary mixtures flowing in a conduit of known cross sectional area, for which the density and the sound speed of each component are known, measuring the sound speed of the mixture and the convection velocity of the mixture provides sufficient information to determine the flow rates of each of the components. The sound speed of the mixture is directly related to the phase fractions of each component in the mixture.
More specifically, to determine the four unknowns, the phase fraction of each component and the flow rate of each component, four equations are used. One equation expresses that the sum of the component (phase) fractions is unity. Another is that the speed of sound in the mixture depends on the speed of sound and density of the individual components and the phase fraction of each constituent. Finally, there are two equations relating the flow velocities of the individual components to the overall mixture flow rate.
The approach of using a single sound speed measurement and a mixture flow rate measurement used in determining the component velocities in the case of a binary mixture, however, cannot be extended to solve the problem of finding the component velocities of a fluid with more than two components. Additional information is required beyond what is provided by a measurement of the speed of sound in the mixture and a measurement of the flow rate of the mixture.
The present invention concerns techniques for the measurement of the flow rate of each component of a multiphase fluid. More particularly, the present invention concerns techniques based on measurements at different locations along a conduit for determining the flow rates of each component of a multiphase fluid within the conduit.
In extracting oil and gas from a formation, it is advantageous to monitor the flow -rates of the different components of the production fluid, usually gas, oil and water.
Such a three-component mixture presents a problem that is often more difficult to solve than that of determining the flow rates for binary mixtures. For a broad class of well mixed binary mixtures flowing in a conduit of known cross sectional area, for which the density and the sound speed of each component are known, measuring the sound speed of the mixture and the convection velocity of the mixture provides sufficient information to determine the flow rates of each of the components. The sound speed of the mixture is directly related to the phase fractions of each component in the mixture.
More specifically, to determine the four unknowns, the phase fraction of each component and the flow rate of each component, four equations are used. One equation expresses that the sum of the component (phase) fractions is unity. Another is that the speed of sound in the mixture depends on the speed of sound and density of the individual components and the phase fraction of each constituent. Finally, there are two equations relating the flow velocities of the individual components to the overall mixture flow rate.
The approach of using a single sound speed measurement and a mixture flow rate measurement used in determining the component velocities in the case of a binary mixture, however, cannot be extended to solve the problem of finding the component velocities of a fluid with more than two components. Additional information is required beyond what is provided by a measurement of the speed of sound in the mixture and a measurement of the flow rate of the mixture.
2 To provide the required additional information, the prior art teaches multiphase flow meters that typically rely on several so called orthogonal measurement systems (said to be a orthogonal because each measurement system provides information that is at least partially independent of the information provided by the other measurement systems). Multiphase flow meters accord'vng to such an orthogonal system approach include meters based on multiple-energy nuclear sources, ultraviolet measurements, capacitance measurements, venturi effect measurements, and infrared measurements.
An alternative approach to determining the component flow velocities of a multiphase fluid in a conduit is to determine the additional required information from multiple point measurements, i.e. from measurements of the same information at different places along the conduit. For example, measurements of the speed of sound in the fluid, the flow rate of the fluid, and the pressure and temperature of the fluid would be made at two or more locations along the conduit. In addition, a multiphase flow model is used to relate the values of the parameters that are measured at one location to those measured at another location. These relationships can provide the additional constraints required to solve for additional parameters. Typically, the equations are nonlinear in several variables. Many methods can be employed to "solve" for the flow parameters of interest.
One class of methods defines an error function based on the -calculated values of the parameters compared to the measured values of the parameters. The flow-related parameters sought (i.e. the component velocities or component phase fractions) are adjusted -iteratively until the error is minimized. In this context, the flow parameters that result in minimizing the value of error function are assumed to be correct.
One company, previously Loke of Norway (now owned by FMC/KOS), is considered by many to have pioneered the general approach based on multiple point measurements for production allocation measurements in oil and gas production facilities. One implementation by Loke of the general approach is their software called Idun, which uses conventional pressure and temperature measurements along with a choke position measurement and irnowledge of the fluid property characteristics to M&C Folio No !'31858W0 3 estimate the component flow rates (or phase fractions). The overall accuracy and robustaesss of the Idun approach is directly influenced by the type and quality of sensors available. *. Aspects of the Idun approach are disclosed in "A Cost Effective Technique for Production WeIl Testini', G. Nerby et al., Offshore Technology Gonference, Houston, Texas, US,1-4 May 1995.
Another approach to determining component flow rates in a multiphase fluid based on making rriessurements at multiple locations along a conduit canying the multiphase fluid is that based on a gradio-venturi system, which includes a ventari meter and employs a remote pressuie sensor located several hundred feet above the venturi. . The pressure difference 'betweett the pressure at the venturi and that at the remote transducer can be related to the flow rate and composition through a multiphase flow model, and can be used in conjunction with the pressure difference due to the flow thorough the vertturi to estimate the component flow rates. Such an approach hag several drawbacks. It requires a venturi, which is intrusive, and it has an accuracy limited in two-phase flow to t10 l of the total flow. Moreover, the accuracy degrades substantiaAy in the presence of any sigaiftcant entrained gas.
What is needed is a system of ineasluing the component velocities aad phase fractions of a fluid that includes at least three components, and that is nonintrusive, sufficiently acaurate, and that does not provide spurious solutions because of insufficient information.
Accordingly, the present invention provides an apparatus and corresponding method for determining component flow rates of a multiphase fluid in a conduit, the fluid consisting of at least three lmown components, the method including the steps of measuring at each of two different positions along the conduit at least four mixture quantities; providing a speed of sound in each of the components at the measured pressures and temperatures; providing a trial value for each of either the component flow rates or the phase fractions; using a predetermined model to calculate values for the measured mutnre quantities based on the trial values for each of either the component flow rates or the phase fractions; using a predetertnined error fitnction to determine an,error value; and using a predetermined optimizing algorithm to detennirite Empf.zei t:27f0I1 AMENDED SHEET
whether the calculated values are acceptable, and, if they are not, to provide a new trial value for each of either the component flow rates or the phase fractions.
In another aspect, the invention provides a method for determining component flow rates of a multiphase fluid in a conduit, the fluid consisting of at least three known components, the fluid having a pressure and temperature at each location of the fluid, the method comprising the steps of:
a) measuring at each of two different positions along the conduit at least four mixture quantities;
b) providing speed of sound values for each of the components at the pressures and temperatures at which the four different mixture quantities are measured;
c) providing a trial value for each of either the component flow rates or phase fractions of the fluid;
d) using a predetermined model to calculate values for the measured mixture quantities based on the speed of sound values and the trial values for each of either the component flow rates or the phase fractions;
e) using a predetennined error function to determine an error value based on differences between the measured mixture quantities and the calculated values for the measured mixture quantities;
f) determining whether the calculated values are acceptable based on the error value; and g) if the calculated values are not acceptable, using a predetermined optimizing algorithm to provide a new trial value for each of either the component flow rates or the phase fractions and repeating steps d through f.
In another aspect, the invention provides an apparatus for determining component flow rates of a multiphase fluid in a conduit, the fluid consisting of at least three known components, the fluid having a pressure and temperature at each location of the fluid, the apparatus comprising:
sensors for measuring at each of two different positions along tile conduit at least four mixture quantities;
a modeler for using a predetermined model to calculate values for the measured mixture quantities based on speed of sound values for each of the components at the pressures and temperatures at which the four different mixture quantities are measured and trial values for 4a each of either the component flow rates or the phase fractions, the trial values and speed of sound values provided to the modeler;
an error function evaluator for using a predetermined error function to determine an error value based an differences between the measured mixture quantities and the calculated values for the measured mixture quantities; and an optimizer for using a predetennined optimizing algorithm to determine whether the calculated values are acceptable based on the error value and, if the calculated values are not acceptable, provide a new trial value to the modeler for each of either the component flow rates or the phase fractions.
In a fiufiher aspect of the invention, the error function is the sum of the squares of the difference between the measured and calculated values at each point.
In another, further aspect of the invention, the four mixture quantities are the sound speed, the flow velocity of the multiphase fluid, the pressure and the temperature.
The above and other objects, features and advantages of the invention will become apparent from a consideration of the subsequent detailed description of a preferred embodiment of the invention given, by way of example, with reference to the accompanying drawings, in which:
Figure 1 is a schematic/block diagram of a system for measuring component flow velocities for a multiphase fluid in a conduit, according to the present invention;
Figure 2 is a set of plots showing the results of a simulation of use of the present invention in case of no free gas being present at an upstream position on a conduit;
Figure 3 is a set of plots showing the results of a simulation of use of the present invention in case of 20% (a phase fraction) free gas being present at an upstream position on a conduit; and Figure 4 is a set of plots showing the results of a simulation of use of the present invention in case of 40% free gas being present at an upstream position on a conduit.
4b The present invention for determining component flow rates (or, equivalently, phase fractions) for a multiphase fluid in a conduit is based on making measurements at different points along the conduit, as indicated in Figure 1. The term flow rate is used here in reference to a point along a conduit to mean mass per unit time flowing past the point along the conduit. The flow rate at a point is related to the average flow velocity at the point by the cross sectional area of the conduit at the point. The tenn phase fraction of a component of a multiphase fluid in a conduit (at a point along the conduit) is used here to indicate the fraction of the total mass of a sample of the fluid (at the point along the conduit) that is due to the presence of the component. Knowing the flow rate of each component is equivalent to knowing the phase fractions and the mixture flow rate. The flow rate for a component, at a point along a conduit, is simply the product of the component flow velocity, the phase fraction for the component, and the total mass per unit length at the point along the conduit.
In the particular embodiment of a system according to the invention indicated in Figure 1, four mixture quantities, namely the sound speed, the (mixture) flow velocity, the pressure and the temperature, are measured at two different positions along the conduit 11, at position A and at position B, using a measurement module 12 at A and a measurement module 13 at B, respectively. In the preferred embodiment, it is assumed that the flow can be modeled as time stationary so that measurements at the two positions at the same time are considered to be related. In other embodiments, however, the assumption that the flow is time stationary can be relaxed. In such embodiments, measurements at position B (downstream of position A) would be assumed to be related to measurements at position A only after a time interval depending on the flow velocity measured at position A, and possibly also at position B.
Still referring to Figure 1, according to the invention values of the measured quantities are provided to a multiphase flow modeler 15, which also receives _as input the speed of sound in each of the individual components. Such component speed of sound values are determined based on pressure and temperature measurements at positions A and B and on an assumed PVT model. The multiphase flow modeler 15 then relates the measurements at position B to the measurements at position A, according to the methodology described below, which involves using successive trial values of the component flow rates (or the phase fractions) and an error function that compares the measured mixture quantities with calculated mixture quantities based on the trial values, and iterating until the value of the error function is acceptably small. In the preferred embodiment of the invention, the system includes an error function evaluator 16 that computes, for each new trial value of the component quantities (component flow rates or phase fractions), each new corresponding value of the error function. An optimizer 17 then detennines whether the calculated values are close enough to the measured values based on a predetermined tolerance, and if not, determines new trial values for the component quantities. The component rates (or the phase fractions) that minimize the error are ultimately provided by the optimizer 17 as the output of the system.
For an N-component mixture in a conduit, the component densities p; , phase fractions O; , and component sound speeds ci are related to the mixture sound speed cmix and mixture density pm,;, by the equations, N
/Omix - ~ OfOi ~ (1) i=1 z 2 =1 and (2) , /PmixCmix 1=1 ~iCf a N
oi -1 . (3) i=1 For fluids in a flexible conduit, the propagation velocity of sound, or any other compressional wave, is influenced by the structural properties of the conduit.
For a fluid in a conduit that is a vacuum-backed tube of radius R, thickness t and Young's modulus E, the measured mixture sound speed Cmeas is related to the infinite domain sound speed cm;x and to the structural properties of the tube via the relation, 2 = 2 + . (4) AmixCmeas /~mixCmix Et However, as mentioned above, for a mixture of three or more components (called here a multi-component mixture, as opposed to a binary component mixture), the speed of sound no longer uniquely identifies the phase fractions of the different components.
Additional information is required.
The above relations combine to provide a relation between measured sound speed c,,te,,s , input fluid and structural properties represented by 6= Et , and unknown phase fractions 9, (where the index refers to one or another of the components), a relation that can be expressed in the form, 1 N Y'ia C2 /Omix,a ~j C2 + O"a (5) max,a i /~ia aa where the index a refers to either position A or position B, and the index i refers to the different components of the fluid.
The invention will now be described in an embodiment that is particularly advantageous for determining the component flow rates (or, equivalently, the phase fractions) for a production fluid being pumped from a formation containing oil, gas, and water. It should however be understood that the invention has more general applicability; as will be clear, nothing about the invention restricts it to fluids of oil, gas, and water.
In the particular case of a three-phase (three-component) fluid consisting of gas, oil, and water, equation (5) becomes, 2 = Pmix,a 2+ Y'oa 2+ ~wa2 + 6a (6) Cmix,a ,PgaCga /PoaCoa waCwa with Oga +Ooa +Owa =1, and where the index a refers to either position A or to position B along a conduit bearing the production fluid.
The invention imposes, as one condition, mass flow continuity of the individual phases, and uses a model for the evolution of gas from the oil. As the production fluid is pumped from position A up to position B, closer to the surface, gas escapes out of the oil/gas solution. The present invention uses a model that calculates the amount of gas at position B based on, among other factors, the amount of gas and oil at position A.
Conceptually, the model tracks the multiphase flow as a four-component mixture, a mixture consisting of water, oil, gas present at position A, and gas evolved from the oil between position A and position B.
The invention will now be described in the particular case of well-mixed flow of a fluid consisting of oil, gas, and water. The extension to the general case, i.e. to the case of flow that is not well-mixed, will be apparent to one skilled in the art.
In the case of well-mixed flow, there is no slip among the individual phases, i.e.
the mixture is homogeneous. In such a mixture, the condition of mass flow continuity from position A to position B becomes, /Omix,A'4AUA /Omix,BABuB . (7) where Aa is the cross section of the conduit at the position indicated by the subscript a, and uQ is the corresponding flow rate. In terms of the individual phases, assuming that the cross sections at position A and position B are equal, the condition of mass flow continuity expressed by equation (7) becomes, UA (~OAPOA + 0gA/0gA + Y'wA/OivA 1- uB \0oBAB + Y'gB/OgB + Y'wB/OwB J = (8) Conservation of the individual phases results in three equations, one equation for each of the three components. For the oil and water components, the invention in the preferred embodiment assumes that the density of the water and oil components is the same at position A and position B, yielding, UAY'oA/OoA - uBY'oB/OoB ~ (9) and uAOwA/0wA = uB5"wB/OwB = (10) These simplifying assumptions can be relaxed without altering the basis of the approach of the present invention., The third equation derives from the conservation of the gas phase. The present invention models two 'mechanisms to account for the change in volumetric flow rate of the gas. First, it is assumed that the gas present at position A under goes an expansion characterized by its compressibility factor. Secondly, it is assumed that amount of gas that evolves from the fluid as it moves from position A to position B is proportional to the amount of oil present at position A. Expressing these relations in terms of volumetric flow rates at standard temperature and pressure (indicated by the subscript R)_ yields the following relations, UAY'gAAA 1 PA ~,R +NR\PB -PA/1'oA ZI AAA =uBOgBAB Z PB T (11) ZA PR A B R B
where f3 is the so-called gas evolution constant, which is defined as the amount of gas (in standard cubic feet) evolved from a barrel of oil as the pressure is lowered one pound per square inch (psi), and Z is the gas compressibility factor, defmed as, Z_ VPTR (12) VRPRT
The gas evolution constant 6 can be determined from a standard differential vaporization test in which the amount of dissolved gas (in standard cubic feet) per barrel of residual oil is reported as a function of pressure (see e.g. PVT and Phase Behavior of Petroleum Reservoir Fluids, by Ali Danesh, Elsevier, 1998). The gas compressibility factor Z is also reported in typical PVT field analyses of black oils.
The first term in equation (11) represents the amount of gas present at position A, and the second term represents the amount of gas evolved from the oil as the fluid moves from position A to position B.
The invention therefore constructs a system of seven equations. The sound speed mixing laws, given by equation (6), provide two equations. The phase conservation conditions, given by equations (9-11), provide three equations.
These seven independent equations are in terms of six unknowns (the component flow rates or the phase fractions at both positions A and B) and so provide more than enough information to uniquely determine the six unknowns. Since the system of equations is over determined, instead of attempting to solve the system of equations (which may include inconsistent conditions), the invention seeks a solution that minimizes an error function. In the preferred embodiment, the invention uses as an error function a simple sum of the squares of the differences between the calculated values of a quantity measured at the two positions A and B and the measured values, i.e. using F to represent a quantity that is measured at position A and position B, where the index i ranges over N such quantities, the invention uses as the error function, E-EI Fcalc,i - Fmeas,i I (13) i=1 The application of this general approach can take many forms, i.e. many different algorithms or procedures may be used to minimize the error function E; how, specifically, the error function is minimized is not the subject of the present invention.
The accuracy of the general approach used in the present invention, that of minimizing an error function as explained above, depends on many factors, including:
the fidelity of the parametric model; the quality, quantity and types of measurements;
the accuracy of the calibration of the model input parameters; the accuracy of the well parameters, including for example the radius R, thickness t and Young's modulus E
of the conduit; and the accuracy of the fluid parameters, including for example the component speeds of sound.
The inventors conducted several Monte Carlo simulations to evaluate the multiphase flow measurement approach developed above. The well parameters were selected to be representative of a producing black oil. The parameters used in the simulation are listed below.
Table 1. Parameters of components.
Fluid Density sound speed (kg/m~3) (ft/see) oil 700 3500 water 1000 5000 gas 1.2*P/Patõ 1200 Table 2. Parameters of the mixture and conduit.
Component-independent parameters Pressure at A 3000 psi Pressure at B 2800 psi Temperature at A and B 100 C
(3, gas evolution constant 0.316 scf/stb/psi Inner diameter of conduit 3.0 inches Wall thickness at A 0.50 inches Wall thickness at B 0.25 inches Figures 2-4 present the results of the Monte Carlo simulations, in which the input parameters and measured parameters are varied in a random, but bounded, manner. The three figures illustrate the ability of the invention, even using the simplifying assumptions noted above, to measure the multiphase flow parameters as a function of oil and water cut for three levels of free gas present at position A (namely 0% for Figure 2, 20% for Figure 3, and 40% for Figure 4) in the presence of the uncertainty for each of the parameters indicated in Table 3 below. The pressure was not directly varied in the simulations; the effect of any uncertainty in pressure is accounted for in the variations of the other model parameters.
Table 3. Uncertainty in model parameters.
Parameter Uncertainty (3, gas evolution constant 2%
Oil sound speed 2%
Gas sound speed 1%
Water sound speed 1%
Oil density 1%
Gas density 2%
Water density 0.5%
Measured sound speed at A and B 0.25%
Measured flow velocity at A and B 0.25%
For each of the simulations indicated in Figures 2-4, a single wall sub (measurement station) is assumed, with a thickness at position A of 0.5", and a thickness at position B
of 0.25", as indicated in Table 2.
The simulations indicate that the invention provides reasonable estimates of oil, gas, and water phase fractions, within 10% of reference values for a reasonable level of uncertainty in both measured and input parameters for free gas volume fractions ranging from 0% at position A to 40% at position A.
The solution appears to be robust to uncertainty; gradually increasing the uncertainty in the input and measured parameters appears to translate into only a gradual degrading in flow measurement accuracy.
It is to be understood that the above-described arrangements are only illustrative of the application of the principles of the present invention. Numerous other modifications and alternative arrangements may be devised by those skilled in the art without departing from the spirit and scope of the present invention, and the appended claims are intended to cover such modifications and arrangements.
An alternative approach to determining the component flow velocities of a multiphase fluid in a conduit is to determine the additional required information from multiple point measurements, i.e. from measurements of the same information at different places along the conduit. For example, measurements of the speed of sound in the fluid, the flow rate of the fluid, and the pressure and temperature of the fluid would be made at two or more locations along the conduit. In addition, a multiphase flow model is used to relate the values of the parameters that are measured at one location to those measured at another location. These relationships can provide the additional constraints required to solve for additional parameters. Typically, the equations are nonlinear in several variables. Many methods can be employed to "solve" for the flow parameters of interest.
One class of methods defines an error function based on the -calculated values of the parameters compared to the measured values of the parameters. The flow-related parameters sought (i.e. the component velocities or component phase fractions) are adjusted -iteratively until the error is minimized. In this context, the flow parameters that result in minimizing the value of error function are assumed to be correct.
One company, previously Loke of Norway (now owned by FMC/KOS), is considered by many to have pioneered the general approach based on multiple point measurements for production allocation measurements in oil and gas production facilities. One implementation by Loke of the general approach is their software called Idun, which uses conventional pressure and temperature measurements along with a choke position measurement and irnowledge of the fluid property characteristics to M&C Folio No !'31858W0 3 estimate the component flow rates (or phase fractions). The overall accuracy and robustaesss of the Idun approach is directly influenced by the type and quality of sensors available. *. Aspects of the Idun approach are disclosed in "A Cost Effective Technique for Production WeIl Testini', G. Nerby et al., Offshore Technology Gonference, Houston, Texas, US,1-4 May 1995.
Another approach to determining component flow rates in a multiphase fluid based on making rriessurements at multiple locations along a conduit canying the multiphase fluid is that based on a gradio-venturi system, which includes a ventari meter and employs a remote pressuie sensor located several hundred feet above the venturi. . The pressure difference 'betweett the pressure at the venturi and that at the remote transducer can be related to the flow rate and composition through a multiphase flow model, and can be used in conjunction with the pressure difference due to the flow thorough the vertturi to estimate the component flow rates. Such an approach hag several drawbacks. It requires a venturi, which is intrusive, and it has an accuracy limited in two-phase flow to t10 l of the total flow. Moreover, the accuracy degrades substantiaAy in the presence of any sigaiftcant entrained gas.
What is needed is a system of ineasluing the component velocities aad phase fractions of a fluid that includes at least three components, and that is nonintrusive, sufficiently acaurate, and that does not provide spurious solutions because of insufficient information.
Accordingly, the present invention provides an apparatus and corresponding method for determining component flow rates of a multiphase fluid in a conduit, the fluid consisting of at least three lmown components, the method including the steps of measuring at each of two different positions along the conduit at least four mixture quantities; providing a speed of sound in each of the components at the measured pressures and temperatures; providing a trial value for each of either the component flow rates or the phase fractions; using a predetermined model to calculate values for the measured mutnre quantities based on the trial values for each of either the component flow rates or the phase fractions; using a predetertnined error fitnction to determine an,error value; and using a predetermined optimizing algorithm to detennirite Empf.zei t:27f0I1 AMENDED SHEET
whether the calculated values are acceptable, and, if they are not, to provide a new trial value for each of either the component flow rates or the phase fractions.
In another aspect, the invention provides a method for determining component flow rates of a multiphase fluid in a conduit, the fluid consisting of at least three known components, the fluid having a pressure and temperature at each location of the fluid, the method comprising the steps of:
a) measuring at each of two different positions along the conduit at least four mixture quantities;
b) providing speed of sound values for each of the components at the pressures and temperatures at which the four different mixture quantities are measured;
c) providing a trial value for each of either the component flow rates or phase fractions of the fluid;
d) using a predetermined model to calculate values for the measured mixture quantities based on the speed of sound values and the trial values for each of either the component flow rates or the phase fractions;
e) using a predetennined error function to determine an error value based on differences between the measured mixture quantities and the calculated values for the measured mixture quantities;
f) determining whether the calculated values are acceptable based on the error value; and g) if the calculated values are not acceptable, using a predetermined optimizing algorithm to provide a new trial value for each of either the component flow rates or the phase fractions and repeating steps d through f.
In another aspect, the invention provides an apparatus for determining component flow rates of a multiphase fluid in a conduit, the fluid consisting of at least three known components, the fluid having a pressure and temperature at each location of the fluid, the apparatus comprising:
sensors for measuring at each of two different positions along tile conduit at least four mixture quantities;
a modeler for using a predetermined model to calculate values for the measured mixture quantities based on speed of sound values for each of the components at the pressures and temperatures at which the four different mixture quantities are measured and trial values for 4a each of either the component flow rates or the phase fractions, the trial values and speed of sound values provided to the modeler;
an error function evaluator for using a predetermined error function to determine an error value based an differences between the measured mixture quantities and the calculated values for the measured mixture quantities; and an optimizer for using a predetennined optimizing algorithm to determine whether the calculated values are acceptable based on the error value and, if the calculated values are not acceptable, provide a new trial value to the modeler for each of either the component flow rates or the phase fractions.
In a fiufiher aspect of the invention, the error function is the sum of the squares of the difference between the measured and calculated values at each point.
In another, further aspect of the invention, the four mixture quantities are the sound speed, the flow velocity of the multiphase fluid, the pressure and the temperature.
The above and other objects, features and advantages of the invention will become apparent from a consideration of the subsequent detailed description of a preferred embodiment of the invention given, by way of example, with reference to the accompanying drawings, in which:
Figure 1 is a schematic/block diagram of a system for measuring component flow velocities for a multiphase fluid in a conduit, according to the present invention;
Figure 2 is a set of plots showing the results of a simulation of use of the present invention in case of no free gas being present at an upstream position on a conduit;
Figure 3 is a set of plots showing the results of a simulation of use of the present invention in case of 20% (a phase fraction) free gas being present at an upstream position on a conduit; and Figure 4 is a set of plots showing the results of a simulation of use of the present invention in case of 40% free gas being present at an upstream position on a conduit.
4b The present invention for determining component flow rates (or, equivalently, phase fractions) for a multiphase fluid in a conduit is based on making measurements at different points along the conduit, as indicated in Figure 1. The term flow rate is used here in reference to a point along a conduit to mean mass per unit time flowing past the point along the conduit. The flow rate at a point is related to the average flow velocity at the point by the cross sectional area of the conduit at the point. The tenn phase fraction of a component of a multiphase fluid in a conduit (at a point along the conduit) is used here to indicate the fraction of the total mass of a sample of the fluid (at the point along the conduit) that is due to the presence of the component. Knowing the flow rate of each component is equivalent to knowing the phase fractions and the mixture flow rate. The flow rate for a component, at a point along a conduit, is simply the product of the component flow velocity, the phase fraction for the component, and the total mass per unit length at the point along the conduit.
In the particular embodiment of a system according to the invention indicated in Figure 1, four mixture quantities, namely the sound speed, the (mixture) flow velocity, the pressure and the temperature, are measured at two different positions along the conduit 11, at position A and at position B, using a measurement module 12 at A and a measurement module 13 at B, respectively. In the preferred embodiment, it is assumed that the flow can be modeled as time stationary so that measurements at the two positions at the same time are considered to be related. In other embodiments, however, the assumption that the flow is time stationary can be relaxed. In such embodiments, measurements at position B (downstream of position A) would be assumed to be related to measurements at position A only after a time interval depending on the flow velocity measured at position A, and possibly also at position B.
Still referring to Figure 1, according to the invention values of the measured quantities are provided to a multiphase flow modeler 15, which also receives _as input the speed of sound in each of the individual components. Such component speed of sound values are determined based on pressure and temperature measurements at positions A and B and on an assumed PVT model. The multiphase flow modeler 15 then relates the measurements at position B to the measurements at position A, according to the methodology described below, which involves using successive trial values of the component flow rates (or the phase fractions) and an error function that compares the measured mixture quantities with calculated mixture quantities based on the trial values, and iterating until the value of the error function is acceptably small. In the preferred embodiment of the invention, the system includes an error function evaluator 16 that computes, for each new trial value of the component quantities (component flow rates or phase fractions), each new corresponding value of the error function. An optimizer 17 then detennines whether the calculated values are close enough to the measured values based on a predetermined tolerance, and if not, determines new trial values for the component quantities. The component rates (or the phase fractions) that minimize the error are ultimately provided by the optimizer 17 as the output of the system.
For an N-component mixture in a conduit, the component densities p; , phase fractions O; , and component sound speeds ci are related to the mixture sound speed cmix and mixture density pm,;, by the equations, N
/Omix - ~ OfOi ~ (1) i=1 z 2 =1 and (2) , /PmixCmix 1=1 ~iCf a N
oi -1 . (3) i=1 For fluids in a flexible conduit, the propagation velocity of sound, or any other compressional wave, is influenced by the structural properties of the conduit.
For a fluid in a conduit that is a vacuum-backed tube of radius R, thickness t and Young's modulus E, the measured mixture sound speed Cmeas is related to the infinite domain sound speed cm;x and to the structural properties of the tube via the relation, 2 = 2 + . (4) AmixCmeas /~mixCmix Et However, as mentioned above, for a mixture of three or more components (called here a multi-component mixture, as opposed to a binary component mixture), the speed of sound no longer uniquely identifies the phase fractions of the different components.
Additional information is required.
The above relations combine to provide a relation between measured sound speed c,,te,,s , input fluid and structural properties represented by 6= Et , and unknown phase fractions 9, (where the index refers to one or another of the components), a relation that can be expressed in the form, 1 N Y'ia C2 /Omix,a ~j C2 + O"a (5) max,a i /~ia aa where the index a refers to either position A or position B, and the index i refers to the different components of the fluid.
The invention will now be described in an embodiment that is particularly advantageous for determining the component flow rates (or, equivalently, the phase fractions) for a production fluid being pumped from a formation containing oil, gas, and water. It should however be understood that the invention has more general applicability; as will be clear, nothing about the invention restricts it to fluids of oil, gas, and water.
In the particular case of a three-phase (three-component) fluid consisting of gas, oil, and water, equation (5) becomes, 2 = Pmix,a 2+ Y'oa 2+ ~wa2 + 6a (6) Cmix,a ,PgaCga /PoaCoa waCwa with Oga +Ooa +Owa =1, and where the index a refers to either position A or to position B along a conduit bearing the production fluid.
The invention imposes, as one condition, mass flow continuity of the individual phases, and uses a model for the evolution of gas from the oil. As the production fluid is pumped from position A up to position B, closer to the surface, gas escapes out of the oil/gas solution. The present invention uses a model that calculates the amount of gas at position B based on, among other factors, the amount of gas and oil at position A.
Conceptually, the model tracks the multiphase flow as a four-component mixture, a mixture consisting of water, oil, gas present at position A, and gas evolved from the oil between position A and position B.
The invention will now be described in the particular case of well-mixed flow of a fluid consisting of oil, gas, and water. The extension to the general case, i.e. to the case of flow that is not well-mixed, will be apparent to one skilled in the art.
In the case of well-mixed flow, there is no slip among the individual phases, i.e.
the mixture is homogeneous. In such a mixture, the condition of mass flow continuity from position A to position B becomes, /Omix,A'4AUA /Omix,BABuB . (7) where Aa is the cross section of the conduit at the position indicated by the subscript a, and uQ is the corresponding flow rate. In terms of the individual phases, assuming that the cross sections at position A and position B are equal, the condition of mass flow continuity expressed by equation (7) becomes, UA (~OAPOA + 0gA/0gA + Y'wA/OivA 1- uB \0oBAB + Y'gB/OgB + Y'wB/OwB J = (8) Conservation of the individual phases results in three equations, one equation for each of the three components. For the oil and water components, the invention in the preferred embodiment assumes that the density of the water and oil components is the same at position A and position B, yielding, UAY'oA/OoA - uBY'oB/OoB ~ (9) and uAOwA/0wA = uB5"wB/OwB = (10) These simplifying assumptions can be relaxed without altering the basis of the approach of the present invention., The third equation derives from the conservation of the gas phase. The present invention models two 'mechanisms to account for the change in volumetric flow rate of the gas. First, it is assumed that the gas present at position A under goes an expansion characterized by its compressibility factor. Secondly, it is assumed that amount of gas that evolves from the fluid as it moves from position A to position B is proportional to the amount of oil present at position A. Expressing these relations in terms of volumetric flow rates at standard temperature and pressure (indicated by the subscript R)_ yields the following relations, UAY'gAAA 1 PA ~,R +NR\PB -PA/1'oA ZI AAA =uBOgBAB Z PB T (11) ZA PR A B R B
where f3 is the so-called gas evolution constant, which is defined as the amount of gas (in standard cubic feet) evolved from a barrel of oil as the pressure is lowered one pound per square inch (psi), and Z is the gas compressibility factor, defmed as, Z_ VPTR (12) VRPRT
The gas evolution constant 6 can be determined from a standard differential vaporization test in which the amount of dissolved gas (in standard cubic feet) per barrel of residual oil is reported as a function of pressure (see e.g. PVT and Phase Behavior of Petroleum Reservoir Fluids, by Ali Danesh, Elsevier, 1998). The gas compressibility factor Z is also reported in typical PVT field analyses of black oils.
The first term in equation (11) represents the amount of gas present at position A, and the second term represents the amount of gas evolved from the oil as the fluid moves from position A to position B.
The invention therefore constructs a system of seven equations. The sound speed mixing laws, given by equation (6), provide two equations. The phase conservation conditions, given by equations (9-11), provide three equations.
These seven independent equations are in terms of six unknowns (the component flow rates or the phase fractions at both positions A and B) and so provide more than enough information to uniquely determine the six unknowns. Since the system of equations is over determined, instead of attempting to solve the system of equations (which may include inconsistent conditions), the invention seeks a solution that minimizes an error function. In the preferred embodiment, the invention uses as an error function a simple sum of the squares of the differences between the calculated values of a quantity measured at the two positions A and B and the measured values, i.e. using F to represent a quantity that is measured at position A and position B, where the index i ranges over N such quantities, the invention uses as the error function, E-EI Fcalc,i - Fmeas,i I (13) i=1 The application of this general approach can take many forms, i.e. many different algorithms or procedures may be used to minimize the error function E; how, specifically, the error function is minimized is not the subject of the present invention.
The accuracy of the general approach used in the present invention, that of minimizing an error function as explained above, depends on many factors, including:
the fidelity of the parametric model; the quality, quantity and types of measurements;
the accuracy of the calibration of the model input parameters; the accuracy of the well parameters, including for example the radius R, thickness t and Young's modulus E
of the conduit; and the accuracy of the fluid parameters, including for example the component speeds of sound.
The inventors conducted several Monte Carlo simulations to evaluate the multiphase flow measurement approach developed above. The well parameters were selected to be representative of a producing black oil. The parameters used in the simulation are listed below.
Table 1. Parameters of components.
Fluid Density sound speed (kg/m~3) (ft/see) oil 700 3500 water 1000 5000 gas 1.2*P/Patõ 1200 Table 2. Parameters of the mixture and conduit.
Component-independent parameters Pressure at A 3000 psi Pressure at B 2800 psi Temperature at A and B 100 C
(3, gas evolution constant 0.316 scf/stb/psi Inner diameter of conduit 3.0 inches Wall thickness at A 0.50 inches Wall thickness at B 0.25 inches Figures 2-4 present the results of the Monte Carlo simulations, in which the input parameters and measured parameters are varied in a random, but bounded, manner. The three figures illustrate the ability of the invention, even using the simplifying assumptions noted above, to measure the multiphase flow parameters as a function of oil and water cut for three levels of free gas present at position A (namely 0% for Figure 2, 20% for Figure 3, and 40% for Figure 4) in the presence of the uncertainty for each of the parameters indicated in Table 3 below. The pressure was not directly varied in the simulations; the effect of any uncertainty in pressure is accounted for in the variations of the other model parameters.
Table 3. Uncertainty in model parameters.
Parameter Uncertainty (3, gas evolution constant 2%
Oil sound speed 2%
Gas sound speed 1%
Water sound speed 1%
Oil density 1%
Gas density 2%
Water density 0.5%
Measured sound speed at A and B 0.25%
Measured flow velocity at A and B 0.25%
For each of the simulations indicated in Figures 2-4, a single wall sub (measurement station) is assumed, with a thickness at position A of 0.5", and a thickness at position B
of 0.25", as indicated in Table 2.
The simulations indicate that the invention provides reasonable estimates of oil, gas, and water phase fractions, within 10% of reference values for a reasonable level of uncertainty in both measured and input parameters for free gas volume fractions ranging from 0% at position A to 40% at position A.
The solution appears to be robust to uncertainty; gradually increasing the uncertainty in the input and measured parameters appears to translate into only a gradual degrading in flow measurement accuracy.
It is to be understood that the above-described arrangements are only illustrative of the application of the principles of the present invention. Numerous other modifications and alternative arrangements may be devised by those skilled in the art without departing from the spirit and scope of the present invention, and the appended claims are intended to cover such modifications and arrangements.
Claims (6)
1. A method for determining component flow rates of a multiphase fluid in a conduit, the fluid consisting of at least three known components, the fluid having a pressure and temperature at each location of the fluid, the method comprising the steps of:
a) measuring at each of two different positions along the conduit at least four mixture quantities;
b) providing speed of sound values for each of the components at the pressures and temperatures at which the four different mixture quantities are measured;
c) providing a trial value for each of either the component flow rates or phase fractions of the fluid;
d) using a predetermined model to calculate values for the measured mixture quantities based on the speed of sound values and the trial values for each of either the component flow rates or the phase fractions;
e) using a predetermined error function to determine an error value based on differences between the measured mixture quantities and the calculated values for the measured mixture quantities;
f) determining whether the calculated values are acceptable based on the error value;
and g) if the calculated values are not acceptable, using a predetermined optimizing algorithm to provide a new trial value for each of either the component flow rates or the phase fractions and repeating steps d through f.
a) measuring at each of two different positions along the conduit at least four mixture quantities;
b) providing speed of sound values for each of the components at the pressures and temperatures at which the four different mixture quantities are measured;
c) providing a trial value for each of either the component flow rates or phase fractions of the fluid;
d) using a predetermined model to calculate values for the measured mixture quantities based on the speed of sound values and the trial values for each of either the component flow rates or the phase fractions;
e) using a predetermined error function to determine an error value based on differences between the measured mixture quantities and the calculated values for the measured mixture quantities;
f) determining whether the calculated values are acceptable based on the error value;
and g) if the calculated values are not acceptable, using a predetermined optimizing algorithm to provide a new trial value for each of either the component flow rates or the phase fractions and repeating steps d through f.
2. A method as in claim 1, wherein the error function is the sum of the squares of the difference between the measured and calculated values at each point.
3. A method as in claim 1 or 2, wherein the four mixture quantities are the sound speed, the flow velocity of the multiphase fluid, the pressure and the temperature.
4. An apparatus for determining component flow rates of a multiphase fluid in a conduit, the fluid consisting of at least three known components, the fluid having a pressure and temperature at each location of the fluid, the apparatus comprising:
sensors for measuring at each of two different positions along tile conduit at least four mixture quantities;
a modeler for using a predetermined model to calculate values for the measured mixture quantities based on speed of sound values for each of the components at the pressures and temperatures at which the four different mixture quantities are measured and trial values for each of either the component flow rates or the phase fractions, the trial values and speed of sound values provided to the modeler;
an error function evaluator for using a predetermined error function to determine an error value based an differences between the measured mixture quantities and the calculated values for the measured mixture quantities; and an optimizer for using a predetermined optimizing algorithm to determine whether the calculated values are acceptable based on the error value and, if the calculated values are not acceptable, provide a new trial value to the modeler for each of either the component flow rates or the phase fractions.
sensors for measuring at each of two different positions along tile conduit at least four mixture quantities;
a modeler for using a predetermined model to calculate values for the measured mixture quantities based on speed of sound values for each of the components at the pressures and temperatures at which the four different mixture quantities are measured and trial values for each of either the component flow rates or the phase fractions, the trial values and speed of sound values provided to the modeler;
an error function evaluator for using a predetermined error function to determine an error value based an differences between the measured mixture quantities and the calculated values for the measured mixture quantities; and an optimizer for using a predetermined optimizing algorithm to determine whether the calculated values are acceptable based on the error value and, if the calculated values are not acceptable, provide a new trial value to the modeler for each of either the component flow rates or the phase fractions.
5. An apparatus as in claim 4, wherein the error function is the sum of the squares of the difference between the measured and calculated values at each point.
6. An apparatus as in claim 4 or 5, wherein the four mixture quantities are the sound speed, the flow velocity of the multiphase fluid, the pressure and the temperature.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US25119100P | 2000-12-04 | 2000-12-04 | |
US60/251,191 | 2000-12-04 | ||
PCT/GB2001/005340 WO2002046705A2 (en) | 2000-12-04 | 2001-12-03 | Method and apparatus for determining component flow rates for a multiphase flow |
Publications (2)
Publication Number | Publication Date |
---|---|
CA2436965A1 CA2436965A1 (en) | 2002-06-13 |
CA2436965C true CA2436965C (en) | 2008-11-25 |
Family
ID=22950871
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA002436965A Expired - Fee Related CA2436965C (en) | 2000-12-04 | 2001-12-03 | Method and apparatus for determining component flow rates for a multiphase flow |
Country Status (6)
Country | Link |
---|---|
US (1) | US6898541B2 (en) |
EP (1) | EP1340055B1 (en) |
AU (1) | AU2002220872A1 (en) |
CA (1) | CA2436965C (en) |
DE (1) | DE60141158D1 (en) |
WO (1) | WO2002046705A2 (en) |
Families Citing this family (76)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7359803B2 (en) * | 2002-01-23 | 2008-04-15 | Cidra Corporation | Apparatus and method for measuring parameters of a mixture having solid particles suspended in a fluid flowing in a pipe |
US7328624B2 (en) * | 2002-01-23 | 2008-02-12 | Cidra Corporation | Probe for measuring parameters of a flowing fluid and/or multiphase mixture |
US7275421B2 (en) * | 2002-01-23 | 2007-10-02 | Cidra Corporation | Apparatus and method for measuring parameters of a mixture having solid particles suspended in a fluid flowing in a pipe |
US7032432B2 (en) * | 2002-01-23 | 2006-04-25 | Cidra Corporation | Apparatus and method for measuring parameters of a mixture having liquid droplets suspended in a vapor flowing in a pipe |
US7212928B2 (en) * | 2002-09-06 | 2007-05-01 | Invensys Systems, Inc. | Multi-measurement vortex flow meter |
CA2513105C (en) * | 2002-11-12 | 2013-01-08 | Cidra Corporation | An apparatus having an array of piezoelectric film sensors for measuring parameters of a process flow within a pipe |
US7165464B2 (en) * | 2002-11-15 | 2007-01-23 | Cidra Corporation | Apparatus and method for providing a flow measurement compensated for entrained gas |
EP1565709B1 (en) * | 2002-11-15 | 2014-03-19 | CiDra Corporation | An apparatus and method for providing a flow measurement compensated for entrained gas |
WO2004048906A2 (en) * | 2002-11-22 | 2004-06-10 | Cidra Corporation | Method for calibrating a flow meter having an array of sensors |
US7096719B2 (en) | 2003-01-13 | 2006-08-29 | Cidra Corporation | Apparatus for measuring parameters of a flowing multiphase mixture |
WO2004063675A2 (en) * | 2003-01-13 | 2004-07-29 | Cidra Corporation | Apparatus and method using an array of ultrasonic sensors for determining the velocity of a fluid within a pipe |
US7343818B2 (en) * | 2003-01-21 | 2008-03-18 | Cidra Corporation | Apparatus and method of measuring gas volume fraction of a fluid flowing within a pipe |
EP1590637B1 (en) | 2003-01-21 | 2008-11-05 | Expro Meters, Inc. | An apparatus and method of measuring gas volume fraction of a fluid flowing within a pipe |
US7058549B2 (en) | 2003-01-21 | 2006-06-06 | C1Dra Corporation | Apparatus and method for measuring unsteady pressures within a large diameter pipe |
US6945095B2 (en) * | 2003-01-21 | 2005-09-20 | Weatherford/Lamb, Inc. | Non-intrusive multiphase flow meter |
US20060048583A1 (en) * | 2004-08-16 | 2006-03-09 | Gysling Daniel L | Total gas meter using speed of sound and velocity measurements |
CA2514696C (en) | 2003-01-21 | 2012-12-11 | Cidra Corporation | Measurement of entrained and dissolved gases in process flow lines |
CA2515551C (en) * | 2003-01-27 | 2013-09-17 | Cidra Corporation | An apparatus and method for providing a flow measurement compensated for entrained gas |
EP1599705B1 (en) | 2003-03-04 | 2019-01-02 | CiDra Corporation | An apparatus having a multi-band sensor assembly for measuring a parameter of a fluid flow flowing within a pipe |
EP1631797A2 (en) | 2003-06-05 | 2006-03-08 | CiDra Corporation | Apparatus for measuring velocity and flow rate of a fluid having a non-negligible axial mach number using an array of sensors |
US7121152B2 (en) * | 2003-06-06 | 2006-10-17 | Cidra Corporation | Portable flow measurement apparatus having an array of sensors |
US7197938B2 (en) * | 2003-06-24 | 2007-04-03 | Cidra Corporation | Contact-based transducers for characterizing unsteady pressures in pipes |
US7245385B2 (en) * | 2003-06-24 | 2007-07-17 | Cidra Corporation | Characterizing unsteady pressures in pipes using optical measurement devices |
US20050050956A1 (en) * | 2003-06-24 | 2005-03-10 | Gysling Daniel L. | Contact-based transducers for characterizing unsteady pressures in pipes |
CA2530596C (en) | 2003-06-24 | 2013-05-21 | Cidra Corporation | System and method for operating a flow process |
US7150202B2 (en) * | 2003-07-08 | 2006-12-19 | Cidra Corporation | Method and apparatus for measuring characteristics of core-annular flow |
WO2005010468A2 (en) * | 2003-07-15 | 2005-02-03 | Cidra Corporation | A configurable multi-function flow measurement apparatus having an array of sensors |
US7152460B2 (en) | 2003-07-15 | 2006-12-26 | Cidra Corporation | Apparatus and method for compensating a coriolis meter |
WO2005010469A2 (en) * | 2003-07-15 | 2005-02-03 | Cidra Corporation | A dual function flow measurement apparatus having an array of sensors |
US7134320B2 (en) * | 2003-07-15 | 2006-11-14 | Cidra Corporation | Apparatus and method for providing a density measurement augmented for entrained gas |
WO2005012843A2 (en) | 2003-08-01 | 2005-02-10 | Cidra Corporation | Method and apparatus for measuring parameters of a fluid flowing within a pipe using a configurable array of sensors |
US7882750B2 (en) * | 2003-08-01 | 2011-02-08 | Cidra Corporate Services, Inc. | Method and apparatus for measuring parameters of a fluid flowing within a pipe using a configurable array of sensors |
US7322251B2 (en) * | 2003-08-01 | 2008-01-29 | Cidra Corporation | Method and apparatus for measuring a parameter of a high temperature fluid flowing within a pipe using an array of piezoelectric based flow sensors |
WO2005015135A2 (en) * | 2003-08-08 | 2005-02-17 | Cidra Corporation | Piezocable based sensor for measuring unsteady pressures inside a pipe |
US7054764B2 (en) * | 2003-09-29 | 2006-05-30 | Air Products And Chemicals, Inc. | Flow monitoring using flow control device |
US7110893B2 (en) * | 2003-10-09 | 2006-09-19 | Cidra Corporation | Method and apparatus for measuring a parameter of a fluid flowing within a pipe using an array of sensors |
US7237440B2 (en) * | 2003-10-10 | 2007-07-03 | Cidra Corporation | Flow measurement apparatus having strain-based sensors and ultrasonic sensors |
US7171315B2 (en) * | 2003-11-25 | 2007-01-30 | Cidra Corporation | Method and apparatus for measuring a parameter of a fluid flowing within a pipe using sub-array processing |
WO2005059479A1 (en) * | 2003-12-11 | 2005-06-30 | Cidra Corporation | Method and apparatus for determining a quality metric of a measurement of a fluid parameter |
US7330797B2 (en) * | 2004-03-10 | 2008-02-12 | Cidra Corporation | Apparatus and method for measuring settlement of solids in a multiphase flow |
WO2005088262A2 (en) * | 2004-03-10 | 2005-09-22 | Cidra Corporation | Method and apparatus for measuring parameters of a stratified flow |
US7367239B2 (en) * | 2004-03-23 | 2008-05-06 | Cidra Corporation | Piezocable based sensor for measuring unsteady pressures inside a pipe |
US7426852B1 (en) | 2004-04-26 | 2008-09-23 | Expro Meters, Inc. | Submersible meter for measuring a parameter of gas hold-up of a fluid |
CA2568349C (en) * | 2004-05-17 | 2013-07-16 | Cidra Corporation | Apparatus and method for measuring compositional parameters of a mixture |
WO2006112878A2 (en) | 2004-09-16 | 2006-10-26 | Cidra Corporation | Apparatus and method for providing a fluid cut measurement of a multi-liquid mixture compensated for entrained gas |
US7389687B2 (en) * | 2004-11-05 | 2008-06-24 | Cidra Corporation | System for measuring a parameter of an aerated multi-phase mixture flowing in a pipe |
FR2880121B1 (en) * | 2004-12-23 | 2007-02-23 | Inst Francais Du Petrole | METHOD FOR DETERMINING THE COMPOSITION OF A POLYPHASIC FLUID |
US7561203B2 (en) * | 2005-01-10 | 2009-07-14 | Nokia Corporation | User input device |
US7725270B2 (en) * | 2005-03-10 | 2010-05-25 | Expro Meters, Inc. | Industrial flow meter having an accessible digital interface |
US7962293B2 (en) * | 2005-03-10 | 2011-06-14 | Expro Meters, Inc. | Apparatus and method for providing a stratification metric of a multiphase fluid flowing within a pipe |
WO2010120258A2 (en) * | 2005-03-17 | 2010-10-21 | Cidra Corporation | An apparatus and method of processing data to improve the performance of a flow monitoring system |
US7657392B2 (en) * | 2005-05-16 | 2010-02-02 | Cidra Corporate Services, Inc. | Method and apparatus for detecting and characterizing particles in a multiphase fluid |
US7526966B2 (en) * | 2005-05-27 | 2009-05-05 | Expro Meters, Inc. | Apparatus and method for measuring a parameter of a multiphase flow |
WO2006130499A2 (en) * | 2005-05-27 | 2006-12-07 | Cidra Corporation | An apparatus and method for fiscal measuring of an aerated fluid |
US7249525B1 (en) | 2005-06-22 | 2007-07-31 | Cidra Corporation | Apparatus for measuring parameters of a fluid in a lined pipe |
EP1899686B1 (en) * | 2005-07-07 | 2011-09-28 | CiDra Corporation | Wet gas metering using a differential pressure based flow meter with a sonar based flow meter |
US7603916B2 (en) | 2005-07-07 | 2009-10-20 | Expro Meters, Inc. | Wet gas metering using a differential pressure and a sonar based flow meter |
WO2007009097A1 (en) * | 2005-07-13 | 2007-01-18 | Cidra Corporation | Method and apparatus for measuring parameters of a fluid flow using an array of sensors |
US20070055464A1 (en) * | 2005-08-17 | 2007-03-08 | Gysling Daniel L | System and method for providing a compositional measurement of a mixture having entrained gas |
US7454981B2 (en) * | 2006-05-16 | 2008-11-25 | Expro Meters. Inc. | Apparatus and method for determining a parameter in a wet gas flow |
US7624650B2 (en) | 2006-07-27 | 2009-12-01 | Expro Meters, Inc. | Apparatus and method for attenuating acoustic waves propagating within a pipe wall |
US7624651B2 (en) | 2006-10-30 | 2009-12-01 | Expro Meters, Inc. | Apparatus and method for attenuating acoustic waves in pipe walls for clamp-on ultrasonic flow meter |
US7673526B2 (en) | 2006-11-01 | 2010-03-09 | Expro Meters, Inc. | Apparatus and method of lensing an ultrasonic beam for an ultrasonic flow meter |
EP2092278A2 (en) | 2006-11-09 | 2009-08-26 | Expro Meters, Inc. | Apparatus and method for measuring a fluid flow parameter within an internal passage of an elongated body |
CN101779105B (en) * | 2007-06-12 | 2012-03-21 | 松下电器产业株式会社 | Flowmeter and fluid supply system |
US7954362B2 (en) * | 2007-07-27 | 2011-06-07 | Multiphase Flow International Llc | Ultrasound multiphase fraction meter and method for determining phase fractions in a multiphase fluid |
DE102008043327A1 (en) * | 2008-10-30 | 2010-05-06 | Endress + Hauser Flowtec Ag | Method and thermal flow meter for determining and / or monitoring at least one, at least dependent on the chemical composition of a measuring medium size |
CA2763461C (en) * | 2009-05-26 | 2017-08-29 | Expro Meters, Inc. | Method and apparatus for monitoring multiphase fluid flow |
WO2011019829A1 (en) | 2009-08-11 | 2011-02-17 | Expro Meters, Inc. | Method and apparatus for monitoring multiphase fluid flow |
DE102010001646A1 (en) * | 2010-02-05 | 2011-08-11 | Endress + Hauser Flowtec Ag | Flow meter for use in measuring system to measure e.g. temperature of carbon dioxide in laboratory, has data processing unit computing amount flow of mass flow and generating signal based on determined amount flow of component of medium |
US8701461B2 (en) | 2011-02-22 | 2014-04-22 | Southern Methodist University | Calibration tube for multiphase flowmeters |
US9383476B2 (en) * | 2012-07-09 | 2016-07-05 | Weatherford Technology Holdings, Llc | In-well full-bore multiphase flowmeter for horizontal wellbores |
CA2846656C (en) * | 2012-11-19 | 2022-06-14 | Manus P. Henry | Net oil and gas well test system |
DE102012112516A1 (en) * | 2012-12-18 | 2014-06-18 | Endress + Hauser Flowtec Ag | Method for verifying the reliability of measured data measured by an ultrasonic flow measurement according to the transit time difference method and ultrasonic flowmeter |
EP2749334B1 (en) | 2012-12-28 | 2018-10-24 | Service Pétroliers Schlumberger | Method and device for determining the liquid volume fraction of entrained liquid |
US9778091B2 (en) * | 2014-09-29 | 2017-10-03 | Schlumberger Technology Corporation | Systems and methods for analyzing fluid from a separator |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4445389A (en) * | 1981-09-10 | 1984-05-01 | The United States Of America As Represented By The Secretary Of Commerce | Long wavelength acoustic flowmeter |
EP0634961A1 (en) * | 1992-04-06 | 1995-01-25 | MOUNTFORD, Norman Duncan Gerard | Ultrasonic treatment of liquids in particular metal melts |
US5343737A (en) * | 1992-09-22 | 1994-09-06 | Joseph Baumoel | Method and apparatus for leak detection and pipeline temperature modelling method and apparatus |
US5415024A (en) * | 1992-12-16 | 1995-05-16 | Marathon Oil Company | Composition analyzer for determining composition of multiphase multicomponent fluid mixture |
RU2250438C9 (en) * | 1998-06-26 | 2005-08-27 | Сидрэ Копэрейшн | Method and device for measuring parameters of fluids in pipes |
US6354147B1 (en) * | 1998-06-26 | 2002-03-12 | Cidra Corporation | Fluid parameter measurement in pipes using acoustic pressures |
US6216091B1 (en) | 1998-09-25 | 2001-04-10 | Panametrics, Inc. | Ultrasonic measurement system with molecular weight determination |
-
2001
- 2001-11-28 US US09/996,626 patent/US6898541B2/en not_active Expired - Lifetime
- 2001-12-03 WO PCT/GB2001/005340 patent/WO2002046705A2/en active Search and Examination
- 2001-12-03 DE DE60141158T patent/DE60141158D1/en not_active Expired - Lifetime
- 2001-12-03 EP EP01999789A patent/EP1340055B1/en not_active Expired - Lifetime
- 2001-12-03 AU AU2002220872A patent/AU2002220872A1/en not_active Abandoned
- 2001-12-03 CA CA002436965A patent/CA2436965C/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
WO2002046705A2 (en) | 2002-06-13 |
DE60141158D1 (en) | 2010-03-11 |
WO2002046705A3 (en) | 2002-10-17 |
CA2436965A1 (en) | 2002-06-13 |
EP1340055A2 (en) | 2003-09-03 |
AU2002220872A1 (en) | 2002-06-18 |
US20020123852A1 (en) | 2002-09-05 |
EP1340055B1 (en) | 2010-01-20 |
US6898541B2 (en) | 2005-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2436965C (en) | Method and apparatus for determining component flow rates for a multiphase flow | |
AU2010282333B2 (en) | Method of measuring multi-phase fluid flow downhole | |
AU2009325212B2 (en) | A method and apparatus for wet gas flow measurements and measurement of gas properties | |
US20090308601A1 (en) | Evaluating multiphase fluid flow in a wellbore using temperature and pressure measurements | |
EP0805959A1 (en) | Method for determination of flow rate in a fluid | |
CA2455250A1 (en) | Non-intrusive multiphase flow meter | |
Guillot et al. | Using simultaneous downhole flow-rate and pressure measurements to improve analysis of well tests | |
Peffer et al. | An improved method for calculating bottomhole pressures in flowing gas wells with liquid present | |
CA2823688C (en) | Method for in-situ calibrating a differential pressure plus sonar flow meter system using dry gas conditions | |
Ünalmis | Sound speed in downhole flow measurement | |
Gudmundsson et al. | Gas-liquid metering using pressure-pulse technology | |
Mursaliyev | Implementation of Virtual Flow Metering Concept in Kashagan Field | |
US7240537B2 (en) | Method for the determination of the wall friction profile along pipes by pressure transients measurements | |
Collins et al. | Evolution of wet gas Venturi metering and wet gas correction algorithms | |
Henry et al. | New applications for Coriolis meter-based multiphase flow metering in the oil and gas industries | |
Richmond et al. | Comparison of implicit and explicit methods for interpreting displacement data | |
Ünalmis | Production Rate Measurement Optimization using Test Separator and In-Well Sound Speed | |
Haugen et al. | Instrumented oilfields: downhole production allocation by passive acoustic listening | |
Al-Khamis et al. | Evaluation of PhaseWatcher Multiphase Flow Meter (MPFM) Performance in Sour Environments | |
Salimov et al. | Allocation Of Zonal Production In Smart Wells: Offshore Abu-Dhabi Case Study | |
Mactaggart | A sour gas leak detection system implementation | |
Selim et al. | Tracking Subsea Gas Wells Performance Without Periodic Production Testing on Test Separator | |
Speed | SPE-201313-MS | |
Gregerson et al. | New developments in wet gas metering to meet offshore challenges | |
Al-Amri et al. | Advanced Utilization of Downhole Sensors for Water-cut and Flow Rate Allocation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request | ||
MKLA | Lapsed |
Effective date: 20171204 |