A flow loop is a pipe configuration in which a single channel of fluid flow is split into two separate flow channels and later combined into a single channel again. Flow loops exhibit self-balancing characteristics when exposed to a constant external force in a specific orientation. The external physical force may be applied to a flow loop such that a differentiation or separation of content of the fluid flow occurs near the splitting section of the flow loop. The separation may cause the fluid flow through each of the separate channels to have differing characteristics from the original single channel. The difference in the content of the separate channels may be used for measurement and other forms of analysis in order to identify the nature of the fluid flowing through the original single channel. This technique is useful where the original fluid has characteristics that render these types of measurements or other forms of analysis impossible.
The external physical force or energy applied to separate the fluid flow may be extracted as the separate channels of fluid flow combine. When the external energy application and extraction are balanced with a flow loop configuration that induces a desired content separation after the split, the flow loop exhibits a unique ability to adapt to changes in the input flow conditions. If the inlet pressure increases for a short period of time, the merging outlet flow will exhibit a proportional pressure increase with a slight time delay while the flow loop maintains the content of each flow path. For example, application of a magnetic field to a flow loop will separate a fluid containing negatively ionized particles from a remainder of the input fluid.
Also, application of a gravity field to a properly structured flow loop will separate fluids having two different densities. Specifically, the denser fluid will flow in the direction of the gravity field through one separate flow path, while the less dense fluid will flow in the opposite direction through the other separate flow path. The two fluids will be combined at the outlet of the flow loop to regain the original state of the input fluid flow. Where the two fluids have significantly different densities, the buoyancy force will also act on the less dense fluid to guide it in the direction opposite to the gravity field. The more dense fluid will lose potential energy as it flows through a lower separate flow path, but the lost potential energy will be converted into either kinetic energy in the form of an increased flow velocity or another form of potential energy in the form of an increased fluid pressure. The opposite energy conversion occurs for the less dense fluid in a higher separate flow path. An overall energy conservation will be maintained for the flow loop after accounting for energy loss due to friction and subsequent heat generation. The pressure is also balanced in the flow loop under the gravity field.
In one embodiment, a method of measuring a flowing volume and density of two liquid phases in a mixture with a gas phase is disclosed. The method comprises providing a flow loop having an inlet vertical column and an outlet vertical column interconnected by a top horizontal section and a bottom horizontal section and feeding a three phase mixture into the inlet vertical column, and wherein a density of the first liquid phase is lower than a density of the second liquid phase. The method includes separating the gas phase from the first and second liquid phases in the inlet vertical column, such that the gas phase flows to an upper section of the inlet vertical column and through the top horizontal section and a liquid mixture of the first and second liquid phases flows to a lower section of the inlet vertical column and through the bottom horizontal section, and measuring a differential pressure of the liquid mixture in the lower section of the inlet vertical column. The method further comprises measuring a flowing volume of the liquid mixture, and calculating a density value for the liquid mixture using the differential pressure. In one embodiment, the diameter of the inlet vertical column is larger than a diameter of the top horizontal section, a diameter of the bottom horizontal section, and a diameter of the outlet vertical column. The differential pressure may be measured using a differential pressure sensor with remote diaphragm seals. The conversion of the differential pressure to the density value for the liquid mixture may be accomplished by dividing the differential pressure by the height between the measuring points of the differential sensor and the gravitational acceleration constant [ρmix=ΔP/(g×h)]. 9The method may further comprise determining a volume percentage of the liquid mixture for the first liquid phase and determining a volume percentage of the liquid mixture for the second liquid phase, and wherein the volume percentage of the liquid mixture for the first liquid phase may be determined by calculating percentage based on the density of the liquid mixture density and the reference density of the first liquid phase and the reference density of the second liquid phase.
In one embodiment, the step of determining the volume percentage of the liquid mixture of the first liquid phase, the volume percentage of the liquid mixture for said first liquid phase is determined by: V1/V=(ρmix−ρ2)/(ρ1−ρ2), wherein: V is the first and second volume; V1 is the first volume; ρmix is the density value for liquid mixture; ρ1 is the first liquid density; and ρ2 is the second liquid density. Also, determining the volume percentage of the liquid mixture for the second liquid phase may be determined by calculating: V2/V=1−V1/V.
In another embodiment, the method may further comprise determining a mass percentage of said first liquid phase in said liquid mixture (m1) and determining a mass percentage of the second liquid phase in the liquid mixture (m1), and wherein the mass percentage of the liquid mixture for the first liquid phase may be determined by calculating: m1/mmix=(V1/V)*(ρmix/ρ1), wherein: mmix is the mass of liquid 1 and liquid 2 mixture in the vertical column volume of V; V1 is the first volume; V is the first and second volume; ρmix is the density for liquid mixture; ρ1 is the first liquid density; and, m1 is the mass percentage of said first liquid phase. The mass percentage of the liquid mixture for the second liquid phase may be determined by calculating: m2/mmix=(V2/V)*(ρmix/ρ2), wherein: mmix is the mass of liquid 1 and liquid 2; m2 is the mass percentage of the second liquid phase; V is the first and second volume; V2 is the second volume; ρmix is the density value for liquid mixture; and, ρ2 is the second liquid density.
As more fully set-out below, the three phase mixture may be under turbulent flow conditions, and wherein the method includes feeding the three phase mixture into a plurality of horizontal pipe sections and beginning the separation of the gas phase from the first and second liquid phases in the horizontal pipe sections before feeding the three phase mixture into the inlet vertical column of the flow loop. The horizontal pipe sections may be configured in a series upstream of the flow loop. Also, in another embodiment, the horizontal pipe sections may be positioned parallel to one another between the split section and the convergence section, and the method includes feeding the three phase mixture through the split section and into the plurality of horizontal pipe sections, and beginning the separation of the gas phase from the first and second liquid phases in the horizontal pipe sections before feeding the three phase mixture through the convergence section and into the inlet vertical section of the flow loop.
In yet another embodiment, a method of measuring a density of two liquid phases in a mixture with a gas phase is disclosed which includes the steps of providing a flow loop, wherein a diameter of an inlet vertical column is larger than a diameter of a top horizontal section, a diameter of a bottom horizontal section, and a diameter of an outlet vertical column. The method further includes feeding a three phase mixture into a series of horizontal pipe sections and beginning to separate a gas phase from a first and second liquid phases in the series of horizontal pipe sections. The three phase mixture is feed into the inlet vertical column of flow loop and the method further comprises continuing to separate the gas phase from the first and second liquid phases in the inlet vertical column, such that the gas phase flows to an upper section of the inlet vertical column and through the top horizontal section and a liquid mixture of the first and second liquid phases flows to a lower section of the inlet vertical column and through the bottom horizontal section. A differential pressure of the liquid mixture in the lower section of the inlet vertical column is measured, and a density value for the liquid mixture using the differential pressure is calculated. Under either turbulent flow and/or laminar conditions, the method further comprises adjusting the differential pressure of the liquid mixture for a friction pressure loss, wherein the friction pressure loss is a pressure drop caused by friction forces in the liquid mixture. Also, the method may include measuring a flow rate of the liquid mixture and calculating a flow velocity of the liquid mixture, estimating the friction pressure loss using the flow velocity, and calculating a gravity differential pressure by subtracting the friction pressure loss from the differential pressure, and then calculating of the density value for the liquid mixture uses the gravity differential pressure. With this embodiment, the method may include measuring a second differential pressure of the liquid mixture in a lower section of the outlet vertical column; and calculating a gravity differential pressure by subtracting the second differential pressure from the differential pressure, then dividing the difference in half and then calculating the density value for the liquid mixture using the gravity differential pressure.
With reference to
Differential pressure measurement system 24 may be positioned on lower section 26 of inlet vertical column 12. Differential pressure measurement system 24 may be a remote seal pressure transmitter. Alternatively, differential pressure measurement system 24 may be achieved by having accurate static (as opposed to differential) pressure sensor with long term stability. According to the teachings of this disclosure, the true differential pressure of the two separate heights in the vertical column 12 is desired. Gas flow meter 28 and valve 30 may be positioned on top horizontal section 16. Gas flow meter 28 may be an eTube™ based gas flow measurement system as described in U.S. Pat. No. 7,653,489 and U.S. Pat. No. 7,623,975, which are both incorporated herein by reference. In another embodiment, gas flow meter 28 may be another differential pressure meter such as a Venturi meter, or could be a linear meter such as a gas turbine meter. Valve 30 may be a ball valve, v-notch ball valve, needle valve or other standard control valve. Liquid flow meter 32 may be positioned on bottom horizontal section 18. Liquid flow meter 32 may be an eTube™ flow meter. In another embodiment, liquid flow meter 32 may be another differential pressure meter such as a Venturi meter, or could be a linear meter such as a turbine meter.
A three phase mixture may be fed through inlet 20 into inlet vertical column 12. The three phase mixture may include a gas phase, a first liquid phase, and a second liquid phase. The first liquid phase may have a lower density than the second liquid phase. Under laminar flow conditions, the gas phase will separate from a liquid mixture of the first and second liquid phases in inlet vertical column 12 due to the force of gravity and the buoyancy force. The gas phase may rise into upper section 34 of inlet vertical column 12, while the liquid mixture may flow into lower section 26. The greater diameter of inlet vertical column 12 may encourage laminar flow for a given flow rate and, in turn, more efficient separation of the gas phase from the liquid mixture. A differential pressure may be measured for the liquid mixture at first vertical position 36 and second vertical position 38 in inlet vertical column 12 using differential pressure measurement system 24. The measured differential pressure value is directly proportional to the density of the liquid mixture. The measured differential pressure value may be used to approximate the density of the liquid mixture, using the following formula:
ρmix=ΔP/(g×h)
wherein
ρmix=density of the liquid mixture,
ΔP=differential pressure,
g=gravitational acceleration,
h=vertical height between measurement points,
more particularly, where ρ is the density of the liquid mixture, ΔP is the differential pressure of the liquid mixture between first vertical position 36 and second vertical position 38, g is the acceleration due to gravity, h is the height difference between first vertical position 36 and second vertical position 38.
In one embodiment, an objective for the liquid side is to measure the volume flow of the two liquids. The density of each liquid will be predetermined by sampling the liquid stream and determining the densities through laboratory analysis. The calculated density will be used (1) in the measurement of the liquid flow rate, and (2) along with known densities of the individual liquids to determine the density percentage of the two liquids.
A physical sample of the liquid mixture may be drawn at the upstream of inlet 20. The laboratory analysis of liquid sample presents the density of the first liquid and the second liquid (ρL1 and ρL2). The two individual density values and the measured mixture density (ρmix) may be used to determine density percentages of the first liquid phase and the second liquid phase volume in the liquid mixture. In one preferred embodiment, the system measures the volume percentages. Mass values can be determined by multiplying by the liquid densities, but are not required. When the fluids with two (2) separate densities are well mixed, they will separate after a period of time in the static condition as follows.
As illustrated in
ρmix h=ρ1 h1+ρ2 h2.
The volume can be set to be;
V=A h
V1=A h1
V2=A h2 where A is the cross sectional area of the cylinder.
In addition;
V=V1+V2. The equation A.
Now further modifying the density relationship using V=A h.
ρmix V=ρ1+V1+ρ2 V. The equation B.
Notice the following.
V is the defined volume in the column for the density measurement.
ρ1 (first liquid density) and ρ2 (second liquid density) are known values from the laboratory analysis.
ρmix is the calculated and known parameter from ρmix=ΔP/(g×h)
Thus, between the equation A and B, there are only 2 unknowns, namely, V1 and V2. They are:
V1=V(ρmix−ρ2)/(ρ1−ρ2)
V2=V−V1=V(1−(ρmix−ρ2)/(ρ1−ρ2)).
Under dynamic condition where the mixed fluids are flowing, ΔP term will include friction factors.
Assume that ΔP=Pbottom−Ptop
Where Ptop=pressure at the top of the column,
Ptop=Pdymanic_top+1/2ρmix v2+Ztop
Pbottom=Pdymanic_bottom+1/2ρmix v2+Zbottom.
Pdymanic_top; the incremental change in static pressure at the measurement point at the top due to dynamic flow (36).
Pdymanic_bottom; the incremental change in static pressure at the measurement point at the bottom due to dynamic flow (38).
ρmix: the average density of the flowing media,
v: the average flow velocity,
1/2ρmix v2: kinetic energy of the fluid,
Ztop: potential energy of the flowing fluid at the top (36),
Zbottom: potential energy of the flowing fluid at the bottom (38).
Therefore, 1/2ρmix v2 is cancelled and results with:
ΔP=Pdymanic_bottom−Pdymanic_top+Zbottom−Ztop.
When the flow rate is small, then the dynamic pressure terms remain small to negligible. When it is fully static, ΔP=Zbottom−Ztop=ρmix g h
Thus under flowing condition, the more generic of the measurement ΔP is:
ΔP=ρmix g h−(Pdymanic_bottom−Pdymanic_top)
ΔP=ρmix g h−∈
Where ∈=(Pdymanic_bottom−Pdymanic_top), and denotes the incremental changes caused by frictions etc. ∈ term consists of the combination of fluid friction and the flow loop physical configuration. It is expected to remain a minor term for laminar flow and at the turbulent flow, it will have to be empirically assessed to ensure the influence of fluid property and mechanical geometry.
Referring again to
In one embodiment, the mixture density is used as a parameter in the flow rate measurement performed by meter 32, when meter 32 is of the differential pressure type such as the eTube. Volume accumulation can be performed for appropriate periods of time, such as Hourly and Daily periods. Using the Volume percentages, accumulated volume can be determined for Liquid 1 and for Liquid 2. As noted earlier, a feature of one preferred embodiment is the measuring of oil and water volumes.
In one embodiment, inlet 20 and outlet 22 may be positioned higher than a midpoint on inlet vertical column 12 and outlet vertical column 14, respectively. This arrangement may allow flow loop 10 to accommodate a larger volume of the liquid mixture than the gas phase. The gas phase may flow through top horizontal section 16, while the liquid mixture may flow through bottom horizontal section 18. In other embodiments, inlet 20 and outlet 22 may be positioned at other heights on inlet vertical column 12 and outlet vertical column 14 to accommodate a different expected volume ratio of the gas phase to the liquid mixture.
Gas flow meter 28 may measure the flow rate of the gas phase flowing through top horizontal section 16. Liquid flow meter 32 may measure the flow rate of the liquid mixture flowing through bottom horizontal section 18. Valve 30 may be used to adjust the amount of the gas phase allowed to flow through top horizontal section 16. This adjustment may be used to ensure that none of the liquid mixture will flow through top horizontal section 16. Valve 30 may even be used to shut off flow through top horizontal section 16 in such a situation. Also, this adjustment may be necessary to ensure that the level of the liquid phase in inlet vertical column 12 does not drop below a minimum liquid level required for accurate operation of differential pressure measurement system 24.
The gas phase may again combine with the liquid mixture in outlet vertical column 14, and the three phase mixture may flow out of outlet vertical column 14 through outlet 22. The energy change in the gas phase and the liquid mixture may balance in the outlet vertical column 14 before the three phase mixture flows out through outlet 22.
Referring to
As seen in
A three phase mixture may be fed into split section 42 in which the three phase mixture is divided into separate flow paths, namely each of the plurality of pre-separation horizontal pipe sections 40. In the plurality of pre-separation horizontal pipe sections 40, a gas phase of the three phase mixture may begin to separate from two liquid phases due to gravity forces, buoyancy forces, and the larger diameter of pre-separation horizontal pipe sections 40 than the surrounding flow lines. In one embodiment, each of pre-separation horizontal pipe sections 40 may have a diameter of approximately six inches. The divided flow of the three phase mixture may be combined in convergence section 44, and the three phase mixture may then flow through inlet 20 into inlet vertical column 12 where the separation of the gas phase from the two liquid phases may continue. As described above, the gas phase may rise to upper section 34 of inlet vertical column 12, while a liquid mixture of the first and second liquid phases may flow to lower section 26.
A differential pressure may be measured for the liquid mixture at first vertical position 36 and second vertical position 38 in lower section 26 of inlet vertical column 12 using differential pressure measurement system 24. The measured differential pressure value may be adjusted for a friction pressure loss caused by frictional forces associated with the turbulent or laminar flow conditions of the liquid mixture in inlet vertical column 12. This adjustment for friction pressure loss may be accomplished by measuring a flow velocity of the liquid mixture in lower section 26 of inlet vertical column 12, and estimating the friction pressure loss using the measured flow velocity. The flow velocity may be calculated from the flow rate determined from meter 32 and the area of the inside of pipe 18,
Flow Velocity=Flow Rate/Pipe Area
The friction pressure loss may be estimated using the friction factor vs. Reynolds Number for pipe flow according to Moody (for instance, see L. F. Moody, Trans. ASME 66, 671 (1944)) The graph covers both the laminar and turbulent flow regime. In one preferred embodiment, only meter 32 is provided with the system. The differential pressure associated with gravity, also referred to as the gravity differential pressure, may be calculated by subtracting the friction pressure loss estimation from the measured differential pressure. The friction pressure loss adjustment may be accomplished by using an eTube™ flow meter as described in U.S. Pat. No. 7,653,489 and U.S. Pat. No. 7,623,975, which are both incorporated herein by reference.
The gravity differential pressure may be used to calculate the density of the liquid mixture using the formula: ρmix=(ΔP+∈)/(g×h). A volume percentage or mass percentage of each of first and second liquid phases in the liquid mixture may be calculated, and a volume may be calculated for each of first and second liquid phases, as described above.
Referring to
The friction pressure loss adjustment required for turbulent flow conditions may be accomplished by measuring the first differential pressure of the liquid mixture in inlet vertical column 12 and measuring the second differential pressure of the liquid mixture in outlet vertical column 14. The first differential pressure of the liquid mixture will be a positive value, which will include a positive first gravity differential pressure (potential energy) and a negative first friction differential pressure. In other words, the first differential pressure of the liquid mixture will include a first gravity pressure increase and a first friction pressure loss. The second differential pressure of the liquid mixture will be a negative value, which will include a negative second gravity differential pressure and a negative second friction differential pressure. In other words, the second differential pressure of the liquid mixture will include a second gravity pressure loss and a second friction pressure loss. The first gravity differential pressure may be calculated (between first vertical position 36 and second vertical position 38 in inlet vertical column 12) by subtracting the second differential pressure value from the first differential pressure value, then dividing the difference in half. This calculation may be expressed by the following formula:
ΔPg1=ρmix g h=1/2(ΔP1−ΔP2)
where ΔPg1 is the first gravity differential pressure (pressure at position 38 and at position 36), ΔP1 is the first differential pressure, and ΔP2 is the second differential pressure. This formula is derived from the following formulas for the first differential pressure and the second differential pressure:
(+)ΔP1=(+)ΔPg1+(−)ΔPf1=ρmix g h−∈
(−)ΔP2=(−)ΔPg2+(−)ΔPf2=ρmix g(−h)−∈
where ΔPf1 is the first friction differential pressure (the term is negative as the pressure drops in the direction of the flow), ΔPg2 is the second gravity differential pressure (the measurement is always negative for pressure at position 52 and at position 50), and ΔPf2 is the second friction differential pressure. Subtracting the second differential pressure from the first differential pressure yields the following formula:
ΔP1−ΔP2=2 ρmix g h
The first friction differential pressure, ΔPf1(∈), and the second friction differential pressure, ΔPf2(∈) will cancel one another in the above formula, leaving the following formula, where ΔPg1 is approximately equal to ΔPg2 (the same as ρmix g h):
ΔP1−ΔP2=(+)ΔPg1+(+)ΔPg2=2(ΔPg1)
The calculated first gravity differential pressure may be used to calculate the density of the liquid mixture using the formula ρmix=ΔP/(g×h) described above. A volume percentage or mass percentage of each of first and second liquid phases in the liquid mixture may be calculated using the pre-determined density for each of first and second liquid phases, as described above.
While preferred embodiments of the present invention have been described, it is to be understood that the embodiments are illustrative only and that the scope of the invention is to be defined solely by the appended claims when accorded a full range of equivalents, many variations and modifications naturally occurring to those skilled in the art from a review hereof.