So that the manner in which the above recited features of the present invention can be understood in detail, a more particular description of the invention, briefly summarized above, may be had by reference to embodiments, some of which are illustrated in the appended drawings. It is to be noted, however, that the appended drawings illustrate only typical embodiments of this invention and are therefore not to be considered limiting of its scope, for the invention may admit to other equally effective embodiments.
Embodiments of the invention relate to multiphase flowmeters capable of determining phase fractions within a multiphase fluid mixture under a broad range of flow conditions such as wet-gas flow. Combining this phase fraction information with a total combined flow rate of the fluid mixture determined based on a sensed velocity of the fluid mixture through a given area enables resolving flow rates for the phase fractions. The phase fractions and hence phase flow rates determined can include oil, gas, and water phases individually or gas and liquid (e.g., oil/water) phases.
Various differential pressure devices utilize a flow nozzle (e.g., the Venturi-based meter 102), an orifice plate, or V-cone to produce changes in velocity and pressure of the fluid flow 110 according to conservation of energy and mass as the fluid flow 110 passes through such devices. Any of these differential pressure devices can therefore provide a differential pressure measurement suitable for applying in equations relating to the conservation of energy and mass in order to determine a property, such as density, of the fluid flow 110. Calculations shown hereinafter refer to the Venturi-based meter 102 as an example way to determine the density while similar modified equations can derive the same results with other ones of the differential pressure devices.
The Venturi-based meter 102 includes first and second ports 116, 118 exposed to pressures of the fluid flow 110 that traverses a constriction formed by a converging inner diameter portion 120 of the conduit 108. As a further example of a different type of flow nozzle, the inner diameter can diverge in some embodiments instead of converge to create a measurable pressure difference. The Venturi-based meter 102 defines a differential pressure sensing meter between the first port 116 disposed upstream of the converging inner diameter portion 120 and the second port 118 located in a throat section downstream of the converging inner diameter portion 120.
The Sonar-based meter 104 can include first, second and third pressure sensing elements 122, 123, 124 distributed along a length of the conduit 108. Spacing between the sensing elements 122-124 enables sensing short duration local pressure variations traveling with the fluid flow (referred to as “flow velocity sensing”) and can also enable sensing acoustic signals traveling at the speed of sound through the fluid flow 110 within the conduit 108 (referred to as “acoustic sensing”). For some embodiments, coils of optical fiber wrapped around the conduit 108 define each of the sensing elements 122-124. Other pressure measuring devices such as piezoelectric or polyvinylidene fluoride (PVDF) based detectors can provide pressure time-varying signals with the Sonar-based meter 104. The acoustic signals and/or the local pressure variations can originate from naturally occurring phenomenon as the fluid flow 110 travels through the conduit 108.
Regardless of the type of the sensing elements 122-124 utilized, interpretation of these signals from the Sonar-based meter 104 enables determination of at least the mixture flow velocity (νm) of the fluid flow 110. For some embodiments, interpreting the signals from the Sonar-based meter 104 permits determination of the speed of sound (SOS, αmix) of the fluid flow 110. U.S. Pat. Nos. 6,354,147 and 6,782,150, which are herein incorporated by reference, describe examples of appropriate calculations for determining the amix and the velocity with similar apparatus that are suitable examples of the Sonar-based meter 104 with the sensing elements 122-124.
The WLR meter 106 can operate based on principles of spectroscopy by relying on differences in absorption between oil and water of near infrared light. In some embodiments, an intrusive probe of the WLR meter 106 within the fluid flow 110 provides a sample region 126 in which input light passes through a portion of the fluid flow 110 and is detected thereafter. Absorption of the input light by the fluid flow 110 attenuates the input light and depends in a wavelength conditioned manner on the contents of the fluid flow 110. As a suitable example of the WLR meter 106, U.S. Patent Publication No. 2006/0186340 and U.S. patent application Ser. No. 11/625,427 (WEAT/0641.P1), which are herein incorporated by reference, describe an infrared optical fiber system capable of determining, for example, the percentage of water and the percentage of oil.
where ρg is known gas density, d is a diameter of the conduit 108 along a throat section at which the Sonar-based meter 104 is disposed, and μg is known viscosity of gas. Using the throat Reynolds number (Red) and the Sonar velocity (νS) that is measured enables calculation of the mixture flow velocity (νm) according to:
νm=νS(1+C0+C1Red−C
where C0, C1, C2, are sonar calibration coefficients. The mixture flow velocity (νm) is typically 1-5% lower than the measured Sonar velocity (νS).
A differential pressure measurement occurs at differential pressure step 204 using the Venturi-based meter 102. Density step 206 calculates a mixture density (ρm) of the fluid flow 110 utilizing the differential pressure measurement from step 204 plugged into an algorithm such as follows based on conservation of energy and mass equations given a total flow rate (Qt) determined by the flow velocity (νm) from step 202. As a prerequisite to determining the total flow rate (Qt), calculation of the extended throat area (At) occurs next using:
Multiplying the results from Equations 2 and 3 yields the total flow rate (Qt) defined by:
Qt=νmAt (Equation 4).
Thereafter, the algorithm proceeds in setting up a final calculation for the density step 206 as represented by forthcoming Equation 12. First, calculation of a Venturi diameter ratio (β) includes:
where D is a diameter of the conduit 108 at an inlet of the Venturi-based meter 102. Second, calculation of a velocity of approach factor (E) entails:
E=(1−β4)−0.5 (Equation 6).
Third, calculation of a pressure ratio (τ) at the inlet and the throat includes:
where P1 is pressure at the first port 116, P2 is pressure at the second port 118, PS is static pressure, and ΔP is the differential pressure. Fourth, a dry-gas Reynolds number (ReD) based on inlet diameter differs from the throat Reynolds number (Re) by factor β according to:
ReD=Redβ (Equation 8).
Fifth, calculation of a dry-gas discharge coefficient (valid for ISA 1932 Nozzle as an example) involves:
Sixth, calculation of a fluid expansibility coefficient (ε) refers to the isentropic expansion of gas at a change in pressure as defined by:
where κ is a known gas isentropic coefficient. Seventh, calculation of a dry gas flow coefficient (Kg) includes:
Kg=CdEε (Equation 11).
Eighth, calculation of the mixture density (ρm) from conservation of energy equation provides:
At liquid fraction determination step 208, using the mixture density (ρm) as determined in Equation 12 by inputting the velocity information (i.e., step 202) from the Sonar-based meter 104 with the response (i.e., step 204) from the Venturi-based meter 102 enables calculating the liquid content (liquid holdup, HL) as follows:
where the gas density (ρg) and a liquid density (ρt) are known or approximated values. An output step 210 includes calculating from results in Equations 4 and 13 a liquid flow rate (Ql) and a gas volume flow rate (GVF, Qg) using the following equations:
Ql=HLQt (Equation 14)
Q
g=(1×HL)Qt (Equation 15).
In addition to the mixture flow velocity (νm), the Sonar-based meter 104 measures at SOS step 305 a SOS of the fluid flow 110. At liquid fraction determination step 308, using the mixture density (ρm) and SOS enables solving for liquid content (liquid holdup, HL) and water-in-liquid ratio (WLR) using the following two equations having these two values as the only unknowns:
ρm=(1−WLR)(HL)ρo+(WLR)(HL)ρw+(1−HL)ρg (Equation 17),
where ρo, ρw and ρg are known densities of oil, water and gas, respectively, αo, αw and αg are a known speed of sound of oil, water and gas, respectively, d is a diameter of the conduit, E is the Young's modulus for the material of the conduit 108, and t is a wall thickness of the conduit 108. Upon determining the liquid content (HL) and the water-in-liquid ratio (WLR), an output step 310 includes calculating an oil flow rate (Qo), a water flow rate (Qw), and a gas volume flow rate (GVF, Qg) using Equation 15 above and the following equations:
Q
o=(1−WLR)(HL)Qt (Equation 18)
Q
w
=WLR(HL)Qt (Equation 19).
At water content step 405, the WLR based meter 106 measures a water-in-liquid ratio (WLR). Next, liquid fraction determination step 408 calculates liquid content (liquid holdup, HL) using the mixture density (ρm) and the water-in-liquid ratio (WLR) according to:
For some embodiments, the liquid content (HL) may be measured directly with the WLR based meter 106 such as described in the previously incorporated U.S. Patent Publication No. 2006/0186340 and U.S. patent application Ser. No. 11/625,427 (WEAT/0641.P1). An output step 410 utilizes the liquid content (HL) and the water-in-liquid ratio (WLR) in Equations 15, 18 and 19 to calculate an oil flow rate (Qo), a water flow rate (Qw), and a gas-volume-fraction flow rate (Qg).
A
i=αoixo+αwixw+αgixg+S (Equation 21)
where Ai is total absorbance at wavelength i and includes chemical (absorption) and physical (scattering) effects, aoi, awi, and agi are absorption coefficients for oil, water and gas respectively at wavelength i, xo, xw, and xg are pathlengths (and hence phase fractions) of oil, water, and gas respectively, and S is a scatter contribution (wavelength independent) to overall absorbance. Making four separate absorbance measurements for four different wavelengths enables solving for four unknowns (xo, xw, xg, and S) in Equation 21.
In a total flow rate measurement step 504, one of either or both of the differential pressure based meter 102 and the Sonar-based meter 104 provide fluid pressure related data to enable determination of a total flow rate of the mixture. For example, Equations 1-4 illustrate use of the Sonar-based meter 104 to determine the total flow rate. The differential pressure based meter 102 relies on solving Equation 12 for the total flow rate with the mixture density determined by the respective phase fractions measured in the phase fraction determination step 502 and known corresponding individual densities. An output step 506 utilizes the individual phase fractions and the total flow rate to calculate an oil flow rate (Qo), a water flow rate (Qw), and a gas-volume-fraction flow rate (Qg) based on the teachings herein.
The system 100 improves accuracy and certainty in the flow rate measurements as described heretofore due to unique aspects of each of the meters 102, 104, 106, which combine to provide sensitivity to certain measurement parameters and insensitivity to other parameters. For example, the WLR based meter 106 provides its results that are independent of GVF and flow rate. The Sonar-based meter 104 measures total volumetric flow rate independent of the WLR and GVF in wet-gas flows. Further, the Venturi-based meter 102 measures momentum dependent on volumetric flow rate and GVF but independent on WLR in wet-gas flows.
Providing multiple ways to utilize at the signal interface circuitry 112 the data received from the meters 102, 104, 106 establishes redundancy that can be beneficial for meter diagnostic purposes or in the event of failure of either the WLR based meter 106 or the Venturi-based meter 102 since options exist to make flow rate measurements using all the meters 102, 104, 106, the Venturi-based meter 102 and the Sonar-based meter 104, or the WLR based meter 106 and the Sonar-based meter 104. For example, a check of the meters 102, 104, 106 in a subsea application of the system 100 can occur utilizing this redundancy. Operation can proceed with functional ones of the meters 102, 104, 106 based on the outcome of such check.
In addition, certain types of flows may lend themselves to a particular configuration of the system 100 relative to other operational modes of the system 100. The Lockhart-Martinelli number describes the wetness of a wet-gas as provided in the following definition and typical values for different wet-gas types:
where νSL and νSg are superficial liquid and gas velocities, and ρL and ρg are liquid and gas densities, respectively. As the value of Lockhart-Martinelli number decreases, the flow approaches to a single-phase gas flow. Larger values refer to flows with increased liquid content. Analysis of the wet-gas types and the capabilities of the meters 102, 104, 106 reveals that when the flow is close to a single-phase gas flow (i.e., Wet-gas Type 1), the speed of sound measurements are expected to be excellent. In this case, one of the dual-configuration solutions (i.e., the Venturi-based meter 102 and the Sonar-based meter 104 or the WLR based meter 106 and the Sonar-based meter 104) can resolve the phase flow rates. When the liquid content is high (i.e., Wet-gas Type 3), the accuracy in the speed of sound measurements tends to decrease, meaning that a triple-configuration (all the meters 102, 104, 106) can be preferred to resolve all the phase flow rates or two phase flow rates can be determined with the combination of the Venturi-based meter 102 and the Sonar-based meter 104 (see,
Ordering and placement of the meters 102, 104, 106 along the conduit 108 can change from the arrangement shown in
While the Sonar-based meter 104 represents one type of flow velocity meter, other devices may be employed to measure flow velocity for utilization in some embodiments of the invention based on the foregoing description. For example, at least two optical sensors spatially distributed along a length of the conduit and designed to detect light interactions with the fluid mixture such that detected time-varying signals can be processed via cross-correlation or an array processing algorithm may provide desired flow velocity information. U.S. patent application Ser. No. 11/421,700, which is herein incorporated by reference, describes such an exemplary multiphase flowmeter.
While the foregoing is directed to embodiments of the present invention, other and further embodiments of the invention may be devised without departing from the basic scope thereof, and the scope thereof is determined by the claims that follow.
This application claims benefit of U.S. provisional patent application Ser. No. 60/826,180, filed Sep. 19, 2006, which is herein incorporated by reference.
Number | Date | Country | |
---|---|---|---|
60826180 | Sep 2006 | US |