CROSS-REFERENCE TO RELATED APPLICATIONS

[0001] This application claims the benefit of U.S. Provisional Application No. 63/382,753 entitled “Carbon Dioxide Multiphase Flow Measurement Based on Dielectric Permittivity,” filed November 8, 2022. This application also claims the benefit of U.S. Non-Provisional Application No. 18/063,158 entitled “Carbon Dioxide Multiphase Flow Measurement Based on Dielectric Permittivity,” filed December 8, 2022.

FIELD

[0002] This patent application relates to apparatus and methods for measuring mass flow rate of a CO2 rich stream. Specifically, apparatus and methods for measuring mass flow rate of a multi-phase CO2 rich stream using dielectric permittivity are described.

BACKGROUND

[0003] Carbon capture, utilization, storage (CCUS), or Carbon capture, sequestration (CCS) applications at surface or subsea rely on accurate measurements of CO2 mass flow rate to manage reporting, inventory, distribution, and other aspects of such operations. Accurate measurements are desired for process control, leak detection, and verification of CO2 quantities for regulatory purposes. In one instance, mass flow rate error of no more than 2.5% may be required.

[0004] Currently, ultrasonic meters, turbine meters and Coriolis meters have been used and/or proposed for mass flow measurement of CO2 rich streams. The ultrasonic meters and turbine meters are used to measure volumetric flow in a dense liquid or supercritical phase, and pressure-temperature-volume (PVT) or equation-of- state (EoS) models or gamma densitometers are used to determine or measure flow density in order to calculate mass flow rate. Such methods are complicated by the complexity in the changing phase behavior of CO2 rich streams with different impurities and in the difficulty in controlling fluid flow phase which can cause multiphase conditions to arise near e.g. CO2 critical point pressure (P _c =73.77 bar) and temperature (T _c =30.98 degC).

[0005] Coriolis flow meters have been tested, and have been determined, in one case, to have 5% error for two-phase flow of a CO2 rich stream. Use of machine learning techniques to classify flow patterns, along with individual phase models have been shown to improve the result obtained from Coriolis flow meters to an error rate of 1.5-2%, but Coriolis devices for large pipes are relatively bulky and expensive.

[0006] Improved methods and apparatus for accurate measurement of mass flow rate of a CO2 rich stream are needed.

SUMMARY

[0007] Embodiments described herein provide a method, comprising measuring a dielectric permittivity of a CO2 rich stream; determining a density of the CO2 rich stream from the measured dielectric permittivity; determining a viscosity of the CO2 rich stream from the measured dielectric permittivity; measuring a pressure drop of the CO2 rich stream flowing through a flow restriction; and determining mass flow rate of the CO2 rich stream using the measured pressure drop, the determined density, and the determined viscosity.

BRIEF DESCRIPTION OF THE DRAWINGS

[0008] Figs. 1A-1 D are schematic views of flow measurement devices that can be used to practice the methods described herein to determine mass flow of a CO2 rich stream.

DETAILED DESCRIPTION

[0009] Mass flow rate measurement with relative error less than about 2.5% is needed for many CO2 flow CCUS or CCS applications at topside or subsea. Such error in flow rates can be realized by using a flow restriction (such as an orifice plate, a flow nozzle, or a Venturi device) to detect a pressure drop across the flow restriction of a flowing CO2-rich stream, and co-locating with the flow restriction a bulk fluid dielectric permittivity sensor. Flow rates from pressure drop measured across a flow restriction are available using known relations. For example, using a Venturi flow restriction, mass flow rate is given by the following known equation: where Q _m is the mass flow rate, C _d is the discharge coefficient, which is a function of Reynolds Number Re, A _T is the Venturi throat cross-sectional area, and is the ratio of throat diameter to inlet diameter. Density of the fluid is p and pressure drop across the Venturi flow restriction is AP. Reynolds Number, at a Venturi flow restriction, is related to mass flow rate by the following equation: where D is the Venturi throat inner diameter, p is the dynamic viscosity, V is the flow velocity at the Venturi throat section, ri is the kinematic viscosity. Discharge coefficients are available as tables or equations incorporating Reynolds Number. Thus, if flow pressure drop and fluid bulk density are known, an iterative process can be followed to calculate the Reynolds Number and discharge coefficient and to converge upon the mass flow rate of a stream.

[0010] Fluid bulk density of a CO2 rich stream can be ascertained from measuring bulk (relative) dielectric permittivity or dielectric constant e _r of the CO2 rich stream. Research shows that, for CO2, bulk fluid dielectric permittivity of CO2 can be fitted to a quadratic function of bulk density to high precision, largely independent of flow pressure and temperature, as follows:

E _r = 1 + Ap + Bp ² (3) where A = 5.099x1 O’ ⁴ m ³ /kg and B = 1 .189x1 O’ ⁷ m ⁶ /kg ² . Equation 3 is shown to be valid at pressures of 1 bar to 300 bar and temperatures of 0 °C to 80 °C, covering the pressure and temperature operating range of the CCUS applications. The critical point CO2 density is p _c = 464 kg/m ³ and the corresponding dielectric permittivity is s _r ,c =1 .262. Research also shows that, unlike the dynamic viscosity (p), the kinematic viscosity of CC has a similarly good correlation with the bulk fluid dielectric permittivity, largely independent of flow pressure and temperature, as follows:

+ ij _c , for Sr Sr.c , gas conditions (4a) , for Sr > sr.c , liquid, dense, supercritical conditions (4b) where Y] _Q is an initial gas-phase CO2 kinematic viscosity at a low density available from literature, L is a regression coefficient of experimental data available from literature, is the kinematic viscosity of CO2 at critical conditions, also available from literature. Dynamic and kinematic viscosity are related as (see, also, Eq.2):

The presence of impurities (such as N2) in CO2 rich streams has little effect on the permittivity-versus-density correlation (Eq.3) and on the permittivity-versus-kinematic viscosity correlation (Eq.4). Thus, if dielectric permittivity of a CO2 rich stream can be ascertained, density and viscosity can be calculated and mass flow rate can be determined from pressure drop measured across a standard flow restriction, such as the Venturi flow restriction exemplified above.

[0011] Generally, microwave transmission and reflection properties of fluids can be used to determine bulk dielectric permittivity of a CO2 rich stream. A microwave transmitter-receiver pair can be used to detect a cut-off frequency of microwaves propagating within a measurement pipe section (a circular waveguide) with a flowing CO2 rich stream. The transmitter-receiver pair are installed such that the transmitter transmits microwaves into the flowing fluid, and the receiver can be disposed at a location opposite from the transmitter, such that the transmitter-receiver pair is oriented along a diameter of the flow pathway of the measurement pipe cross section. The transmitter and receiver can also be disposed in locations that are not directly opposite, one from the other. The transmitter may transmit microwaves at a suitable range of frequencies in the (circular waveguide) measurement pipe section, and attenuation of the signals is measured by the receiver at the corresponding range of frequencies. The peak frequency at which signal attenuation is minimized, sometimes referred to as a “cutoff frequency” (of a dominant propagation mode such as the TE11 mode), can be related to bulk dielectric permittivity using one of the following equations: f _c = K _mode - U /c=y, (6a) where f _c is a dominant mode cutoff frequency, K _mode is the dominant mode coefficient (K _mode = 0.586 for the TEn mode, K _mode = 0.97 for the TE21 mode,), f _{c 0} is the cutoff frequency measured at a standard condition where permittivity is equal to unity, such as in an empty pipe, c ₀ is the vacuum speed of light, D is the flow pathway diameter of the measurement pipe section, and s _r is the bulk fluid dielectric permittivity of the fluid. Equations 6a and 6b can be used together, or only one of equations 6a and 6b can be used. Equation 6b can be used if cutoff frequency at the standard condition is known. Otherwise, equation 6a can be used to determine the bulk fluid dielectric permittivity. As an example, with a sufficiently high-frequency scan resolution, the measurement of cutoff frequency (and hence the permittivity and density determination) can be obtained with accuracy in the range of 489 MHz to 692 MHz for the TEn mode (with sr from 2 to 1 ) for a 10-inch diameter pipe.

[0012] Two microwave transmitter-receiver pairs can be used for an accurate determination of fluid dielectric permittivity from a combined drift-immune transmission attenuation measurement. In a CO2 rich flow measurement pipe section, two microwave transmitter-receiver pairs can be installed at two different pipe cross sections, or at the same pipe cross section, with the said cross section(s) being at the upstream, or at the restriction, or at the downstream of a flow restriction, such as at the inlet section, or the throat section, or the outlet section of a Venturi device. In one embodiment of two microwave transmitter-receiver pairs being installed at two different pipe cross sections, the two transmitter-receiver pairs can be installed along two intersecting diameters of the flow pathway or along two non-intersecting diameters of the flow pathway. The four transmission measurement data of the two pairs can be used to obtain one transmission measurement to improve accuracy and stability by compensating for instrument gain drift. A first transmitter-receiver pair T1 and R1 can be installed at a first diameter of the flow pathway and a second transmitter-receiver pair T2 and R2 can be installed at a second diameter of the flow pathway that does not intersect with the first diameter. At a suitably high transmitting frequency (that is below the cutoff frequency), the two transmitter-receiver pairs enable four measurements of attenuation at two substantially different transmitterreceiver (far and near) spacings, a T1-R1 measurement (directly across the flow pathway with a near-spacing), a T1-R2 measurement f ₁₂ (across the flow pathway and axially displaced with a far-spacing), a T2-R2 measurement f ₂₂ (directly across the flow pathway with substantially the same near-spacing), and a T2-R1 measurement f ₂₁ (directly across the flow pathway and axially displaced with substantially the same far-spacing). A compensated differential measurement can be determined from the four measurements, as follows:

Dielectric permittivity is quadratically correlated to the compensated differential measurement, as follows: where f _cdm = - Adm (empty pipe), a and b can be determined by correlation of the modeling results from 3D electromagnetic simulations or from experiments. Note that a and b are dependent on the transmitting frequency, the transmitter-receiver antenna types (magnetic dipole or electric dipole), and the transmitter-receiver antenna spacings, and the pipe diameter.

[0013] Dielectric permittivity can also be determined using low-frequency capacitance sensors, as is known in the art. Such sensors can be used, optionally with electrical capacitance tomography techniques known in the art, to determine dielectric permittivity instead of, or in addition to, measurements using microwave sensors.

[0014] Research shows that the above relations are durable where some impurities are present in a carbon-captured CO2 stream. The relations above have been shown to hold closely for pure CO2 and pure N2 over the pressure and temperature ranges normally encountered in a CCUS process. A post-combustion and pre-com bustion captured CO2 rich stream typically has at least 95 vol% CO2, where N2 and O2 are present in a small quantity up to 1.3 vol% and where water is present in a quantity up to 600 ppmv. Where a continuous flow of a CO2 rich stream contains, or is expected to contain, free water, the free water (having a high permittivity s _r = 40 to 80, depending on salinity and temperature) can be detected using a microwave reflection sensor. The microwave reflection sensor is installed at an internal wall of the flow pathway to resolve a large dielectric permittivity of any liquid water along the wall. A probe of the microwave reflection sensor may be optionally configured to contact any liquid water flowing along the wall to provide near-wall permittivity and conductivity readings that can be used to detect the presence of free water that may cause flow-assurance issues (such as the risks of formation of CO2 ice-like hydrates or the risks of carbonic-acid pipe corrosion), and to resolve water salinity (calculated from near-wall water-rich permittivity and conductivity readings). Dielectric permittivity of a CO2 stream with dispersed water can be ascertained using a suitable dielectric mixing model, such as a simplified Ramu-Rao model, as follows: where a is the volumetric water holdup (e,g, up to 600 ppmv) that may be determined from the s _r measurement from the microwave reflection sensor, or from the microwave transmission cutoff-frequency or from the compensated differential attenuation, and s _C02 is the dielectric permittivity of the non-water portion of the CO2 stream that can then be used to calculate mass flow rate using the relations above. A microwave reflection sensor can be used with any combination of the other microwave sensors described above. The microwave reflection sensor may be located substantially in the same plane as one or more microwave transmitterreceiver pairs, or in a different plane.

[0015] Figs. 1A-1 D are schematic views of flow measurement devices that can be used to practice the methods described herein to determine mass flow of a CO2 rich stream. Fig. 1A shows a flow measurement device 100 that uses a single microwave transmitter-receiver pair 102. The microwave transmitter-receiver pair 102 is shown located at an outlet side of a flow restriction 104 of the flow measurement device. An optional microwave reflection sensor 106 is included in the flow measurement device 100. Here, the microwave reflection sensor 106 is located at an inlet side of the flow restriction 104. The flow restriction 104 may be a flow nozzle flow restriction or another flow restriction such as a Venturi-style or an orifice plate whose pressure drop properties are closely related to flow rate.

[0016] Fig. 1 B shows a flow measurement device 120 that uses two pairs of microwave transmitter-receivers 122A and 122B. In this case, the pairs 122A and 122B are located along intersecting diameters of the flow measurement device 120. The view of Fig. 1 B is slightly vertically angled to show the positional relationship of the two pairs 122A and 122B.

[0017] Fig. 1 C shows a flow measurement device 130 that uses the two pairs of microwave transmitter-receivers 122A and 122B located along non-intersecting diameters of the flow measurement device 130. The flow measurement devices 120 and 130 use the same flow restriction as the flow measurement device 100, and both also have an optional microwave reflection sensor 106, as in the flow measurement device 100. As noted above, each of the flow measurement devices 100, 120, and 130 can use a Venturi-style flow restriction as illustrated in Fig. 1 D where the two pairs of microwave transmitter-receivers 122A and 122B are located along nonintersecting diameters of the flow measurement device 140 at the restriction (throat section) of a Venturi-style device. Each of the flow measurement devices 100, 120, 130 and 140 has an optional temperature sensor 160 located at an outlet side of the flow restriction in each case. Each device 100, 120, 130 and 140 also has a differential pressure instrument 170 to measure the pressure and pressure drop across the flow restriction. A digital processing system can be configured to receive signals from the various sensors of the flow measurement devices 100, 120, 130 and 140 and to calculate mass flow rate of the CO2 rich stream flowing therein using the relations described above.

[0018] The preceding description has been presented with reference to present embodiments. Persons skilled in the art and technology to which this disclosure pertains will appreciate that alterations and changes in the described structures and methods of operation can be practiced without meaningfully departing from the principle, and scope of this present disclosure. Accordingly, the foregoing description should not be read as pertaining only to the precise structures described and shown in the accompanying drawings, but rather should be read as consistent with and as support for the following claims, which are to have their fullest and fairest scope.