CA2101979A1 - Gas density meter and method - Google Patents
Gas density meter and methodInfo
- Publication number
- CA2101979A1 CA2101979A1 CA002101979A CA2101979A CA2101979A1 CA 2101979 A1 CA2101979 A1 CA 2101979A1 CA 002101979 A CA002101979 A CA 002101979A CA 2101979 A CA2101979 A CA 2101979A CA 2101979 A1 CA2101979 A1 CA 2101979A1
- Authority
- CA
- Canada
- Prior art keywords
- pressure
- chamber
- gas
- sample gas
- flowing
- 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.)
- Abandoned
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F15/00—Details of, or accessories for, apparatus of groups G01F1/00 - G01F13/00 insofar as such details or appliances are not adapted to particular types of such apparatus
- G01F15/02—Compensating or correcting for variations in pressure, density or temperature
- G01F15/04—Compensating or correcting for variations in pressure, density or temperature of gases to be measured
- G01F15/043—Compensating or correcting for variations in pressure, density or temperature of gases to be measured using electrical means
- G01F15/046—Compensating or correcting for variations in pressure, density or temperature of gases to be measured using electrical means involving digital counting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N7/00—Analysing materials by measuring the pressure or volume of a gas or vapour
Abstract
Abstract of the Disclosure A method and apparatus for determining a volume correction ratio for a gas (12) by measuring the temperature and pressure of sample gas (20) tapped from a pipeline (14) as it undergoes controlled expansion. Sample (20) gas flows through first chamber (22) and mechanical pressure regulator (24), where the pressure of the sample gas (20) is reduced. The pressure decay rate of sample gas (20) is measured in the second chamber (26) and compared to a reference gas (31) pressure decay rate. The volume correction ratio is determined from the rate of change of pressure in the first chamber (22), the absolute pressure of the pipeline gas (12), and the absolute temperature of the pipeline gas (12). The volume correction ratio is used to determine the flowing condition density Pf from a measured base condition density Pb.
Description
GAS DENSITY METER AND METHOD
Back~round of tlle Invention In one aspect, the present invention relates to a method and instrumentation that measures, in real time, the tensity of a gas flowing in a pipeline. In this aspect, the invention can also measure a base condition density of the same gas where the base condition densitv corresponds to a ` 5 density (i.e. mass/volume) of the gas determined as if the gas were at some defined base temperature and pressure condition. In a broader aspect, the Pf present invention determines a ratio p of the measured flowing condition density compared to the base condition density. This ratio can be referred to as a volume correction ratio and can be used to translate a measured pipeline gas volumetric flowrate to a corresponding base condition volumetric flow rate.
Heretofore, accurate measurement of gas flowrates has often been difficult because accurate determination of gas density is important when measuring gas flowrates and gas supercompressibility effects gas density in ways that are not detectable by volumetric flowmeters. For instance, with differerltial pressure flowmeters (e.g., an orifice plate meter), the density term appears directly in the volume flow equation (i.e. ~). Also, with other flowmeters (such as turbine flowmeters), as well as with differential pressure flowmeters, the translation of volumetric flowrate to base conditions (i.e. volume correction ratio) is directly related to the density of the flowing gas.
Yet, at elevated gas pressure, computing gas density from the ideal gas - law is inadequate and can be substantially wrong because the ideal gas law -'~ does not properly account for supercompressibility. If gas properties are well 25 known, the difficulty can be partially alleviated because gas . supercompressibility can be directly computed and accommodated in the - ~ calculations. However, if the gas is a variable mixture of unlike molecules, gas properties necessary to compute supercompressibility are unknown and difficult to measure.
.....
., ..... .
.
~, , . .
~'~':- , . .
:
, 210~979 Gas densities, volume correction ratios, and volumetric flowrates at base conditions are normally calculated using flow computers from contemporaneous measurements of several gas parameters. Generally, for determining volumetric flowrates at base conditions, the pipeline gas volumetric flowrate at pipeline conditions is measured, the gas temperature - and pressure at pipeline conditions are measured, and the composition of the gas is measured. Gas composition is normally measured by a chromatograph. From this measured data, gas supercompressibility at both pipeline and base conditions is calculated and from that the density of the gas at pipeline and base conditions is calculated. When operating pressure is elevated (and supercompressibility effects cannot be ignored) gas supercompressibility is usually estimated from either virial equations of state or from correlations such as NX-19. Using virial coefHcients is severely limited because the virial coefficients are functions of temperature, `` 15 pressure, and composition, are largely unknown, and have significant real time uncertainties. Also, using correlations such as NX-19 is often inaccurate because the correlations can be inaccurate for many compositions.
An alternate method for determining volumetric flowrates at base conditions involves the Gerg Equations. The Gerg Equations estimate - supercompressibility and density from knowledge of the heating value, the density of the gas at base conditions and the percentage of carbon dioxide and nitrogen in the exhaust of burned gas. Such measurements can be - made using PMl's GB 3000 product (Precision Measurements, Inc., Tulsa, . 25 Oklahoma). The Gerg Equations allow more rapid computation of ;~ supercompressibility than the above described composition methods and are, therefore, preferred in applications at normal natural gas pipeline ~` pressures. However, as with composition methods, extenshe use of a flow ; computer is required to solve the Gerg Equations.
` 30 Moreover, each of the measurements recited above in describing both the composition method or the Gerg Equations method introduce the potential for measurement error. The aggregation of such errors can substantially influence accuracy. Because of this, it is common practice to frequently calibrate and mainhin each individual measurement device.
,~., ,, ":,.
., :
, . , :' , .. .
.
-, 210197~
In co-pending patent applications, Serial No. 07/793,753, filed on November 18,1991, and Serial No. 07/787,188, filed on November 4,1991, inventions are disclosed that among other things can determine the volume correction ratio of a pipeline gas from energy type measurements.
These inventions involve the measurement of energy flowrate and energy content of a sample of gas tapped from the pipeline. A base condition volumetric flowrate of the sample gas can then be determined by the ratio of the energy flowrate of the sample gas to the energy content of the sample gas. These inventions also measure the ratio of the mass flowrate of the pipeline gas in the pipeline compared to the mass flowrate of the sample gas hpped from the pipeline. The volume correction ratio (or the base condition volumetric flowrate of the pipeline gas) can be calculated easily from the mass flowrate ratio and the base condition volumetric flowrate of the sample gas.
While it is desirable to measure energy flowrate and energy content in many applications, in applications where energy flowrate or energy content are not required (but where base condition volurnetric flowrates are desirable), the energy measurement equipment provides excessive expenæ.
Further, the inventions in these two co-pending patent applications are directed to monitoring combustible pipeline gas and are thus inappropriate for monitoring flows of non-combustible gas.
,. ~
~.
The object of the present invention is to provide an improved ` 25 method and apparatus for measuring a volume correction ratio and a flowing condition density of a pipeline gas.
Another object of the present invention is to make these measurements without making energy type measuremenh.
~-~ The preænt invention is a method and apparatus that can determine 30 a volume correction ratio for a gas (i.e. a ratio of a flowing condition density compared to a baæ condition density p--), and in another aspect can also determine a flowing condition density of a gas. It can make these ~! determinations accurately in real time, without making several highly .-:, .. .
-`` 21~1979 variable measurements, without making any energy type measurements, and without making excessive numeric calculations or correlations in a flow computer.
The present invention measures the volume correction ratio by measuring the temperature and pressure of sample gas tapped from a pipeline as it undergoes controlled expansion. Before the expansion, the sample gas is trapped to assure that it has the same composition as the pipeline gas. The rate of expansion of the trapped sample gas can be controlled in several ways, but it is preferred that the expansion be controlled by establishing mass flow at sonic conditions through a sonic nozzle. It is important that the expansion of the sample gas occur while the trapped sample gas is at substantially the same temperature as the pipeline gas being monitored. In a preferred embodiment, this is achieved by immersing a portion of the meter within the pipeline.
In particular, sample gas is hpped from the pipeline and flows to a first chamber having a fixed volume. A ffrst valve controls the flow of sample gas into the first chamber. When the first valve is opened, sample gas flows into the first chamber preferably until the pressure in the first charnber equals the pressure in the pipeline.
The first valve is then closed and stops further pipeline gas from flowing into the first chamber. In order to measure the volume correction ~, ratio, the sample gas temperature must be maintained at subshntially the same temperature as the pipeline gas temperature. Sample gas flows from ~ ~ the first chamber and then through a flow restrictor located downstream of ;~ 25 the first chamber. The flow restrictor is preferably a sonic nozzle. The sample gas pressure at the inlet of the flow restrictor is established by a ~- pressure regulator so that gas flows from the first chamber at a substantially constant flowrate. If a sonic nozzle is used, the established pressure is sufficient so that the mass flow through the sonic nozzle is critical (i.e., thesonic nozzle inlet pressure is sufficient so that the mass flowrate through `~- the sonic nozzle is proportional to the density and sonic velity of the gas ` ~- and nozzle cross-sectional area).
~, Flowing sample gas from the first chamber when the first valve is closed causes the pressure in the first chamber to reduce over time. The rate , .
of change of the pressure in the first chamber is measured as the pressure in ~',.
, .
"
~' `
the first chamber decays. The temperature of pipeline gas being monitored is also measured.
The volume correction ratio is determined from the rate of change of pressure in the first chamber, the absolute pressure of the pipeline gas, and 5 the absolute temperature of the pipeline gas. The determination can be done without resorting to estimates and approximations for supercompressibility.
The present invention eliminates the need to estimate or consider the effects of supercompressibility when determining gas density because 10 the pressure rate of change measurements in the first chamber are made and used in density calculations directly.
In another aspect of the invention, the volume correction ratio (i.e.
Pf the ratio p of the flowing condition density pf to the base condition density Pb) is used to determine the aowing condition density pf from a measured 15 base condition density Pb. The preferred method and apparatus for measuring the base condition density Pb are disclosed generally in U.S.
Patent No. 4,677,~41 to Kennedy, et al. It is preferred that the sonic nozzle ofthe present invention for measuring the volume correction ratio also be used as a pore where the square of the mass flow rate through the pore is 20 inversely proportional to the density of gas flowing through the pore when the present invention measures the base condition density Pb.
To determine the base condition density Pb, it is preferred that the .-first valve be maintained open to allow sample gas to aow into the first ~`chamber and therefrom; through the pressure regulator; and then into a 2S second chamber also having a fixed volume. Also preferably, the flow into the second chamber is controlled by a second valve. The flow of sample gas from the second chamber is restricted by the pore. The mass flowrate through the pore is measured and from that mass flowrate measurement and a similar measurement of a reference gas which is made periodically, 30 the base condition density Pb of the sample gas can be determined.
~^In another preferred embodiment of the invention, volumetric aow - ~signals are received from a volumetric flowmeter in the pipeline and the volumetric and mass flowrate in the pipeline at base conditions or pipeline conditions can be determined.
. .
,.
21~1~79 The aspect of the invention that measures the base condition density Pb ~i.e., U.S. Patent No. 4,677,841 in conjunction with a reference gas) has proven to be accurate. Also, using a reference gas avoids the need to force the sample gas temperature and pressure to base conditions at the pore, such 5as 14.71 psia and 60F, in order to make a base condition density Pb measurement. This is because the density of the gas flowing through the pore is measured at ambient conditions and then translated to the base condition density through a prior knowledge of the reference gas base condition.
10The foregoing objects and advantages of the present invention will appear from the following description. In the description, references are made to the accompanying drawings which form a part hereof and in which a preferred embodiment of the present invention is shown by way of illustration. Such embodiment does not necessarily represent the full scope 15of the invention however.
Brief Description of the Drawin~s Fig. 1 is a schematic drawing of the fundamental apparatus of the ; preferred embodiment of the present invention;
Fig. 2 is a plot of the pressure in a first fixed-volume chamber as a 20function of time when the present invention as shown in Fig. 1 is operating;
Fig. 3 is a plot of the normalized rate of change of pressure in the first chamber as a function of absolute pressure in the first chamber when the present invention as shown in Fig. 1 is operating;
~^j 25Fig. 4 is a schematic drawing of the preferred way of using a`' differential pressure volumetric flowmeter with the present invention as shown in Fig.1;
Fig. 5 is a schematic drawing of the preferred way of using a linear volumetric flowmeter with the preænt invention;
~-, 30Fig. 6 is a schematic drawing showing the preferred method of installing the present invention;
Fig. 7 is a detailed diagram of a preferred embodiment of a first valve of the present invention;
,~, .
:' 2101~79 Fig. 8 is a detailed cross-sectional view of a preferred embodiment of a pore/sonic nozzle of the present invention; and Fig. 9 shows plots of experimental data for ethane gas and methane gas that relates to a preferred embodiment of the present invention.
Detailed Description of the Preferred Embodiments The present invention includes many aspects, including the aspect of measuring a volume correction ratio [Pf) without measuring a base condition density Pb- The preferred aspect of the present invention, however, involves the measurement of Pb and is, therefore, discussed first In this description of the preferred embodiments of the invention, the apparatus of the preferred embodiment and its operation is first described. Then an analysis of the invention in scientific and mathematical terms is presented. Finally, operation of the present invention with a volumetric flowmeter is described.
Apparahls of the Preferred Embodiments In Fig. 1, a gas density meter of the present invention is designated generally as 10. The solid lines connecting elements of the gas meter 10 represent flow of sample gas and the dashed lines represent electrical signals.
The gas density meter 10 has two modes of operation: a base condition mode and a flowing condition mode. In the base condition mode, the meter 10 measures a base condition density Pb of a pipeline gas 12 flowing through a pipeline 14 (i.e., the density of the gas at some designated " base temperature, pressure, and composition). When the meter 10 is operating in the base condition mode, valve 16 is open and valve 18 cycles open and closed. In the flowing condition mode, the meter 10 measures a ~ flowing condition density pf of the pipeline gas 12 flowing through the pipeline 14 (i.e., the density of the gas 12 at the temperature and pressure that it is flowing through the pipeline 14). When the meter 10 operates in ~- 30 flowing condition mode, valve 16 is closed and valve 18 is open. Acomplete measurement cycle constitutes operation in both the base condition and the flowing condition modes. Note that in applications " 2~0~79 where only a volume correction ratio ~Pf)is desired and it is not necessary to determine either the base or the flowing condition densities, it is not necessar,v to operate in the base condition mode.
Referring in particular to Fig. 1, the sample gas 20 is tapped from the 5 pipeline 14 and flowed into a first fixed-volume chamber 22 when valve 16 is open. The volume of the first chamber 22 is small, about ~10 cubic centimeters. The sample 20 gas fills the first chamber 22 until the pressure in the first chamber 22 is equal to the pipeline gas 12 pressure in the pipeline14. The sample gas 20 then flows to a mechanical pressure regulator 24 10 where the pressure of the sample gas 20 is reduced; and if valve 18 is open, flows into a second fixed-volume chamber 26. The sample gas 20 then flows from the second chamber 26 through a pore 28.
- When operating in the base condition mode, the flow of sample gas 20 through the pore 28 is controlled using valve 18. The pore 28 is operated 15 in a special operating mode so that the square of the sample gas 20 mass `~ flowrate through the pore 28 is inversely proportional to the density of the ` ~s sample gas 20 at the time it is flowing through the pore 28. This special operating mode was discovered by Kennedy and is described in U.S. Patent No. 4,677,841.
There are many ways to determine the sample gas 20 mass flowrate through the pore 28, but the preferred way is to measure a sample gas 20 ;. pressure decay rate in the second chamber 26 and compare this to a reference gas 31 pressure decay rate, as now explained.
Opening valve 18 allows sample gas 20 into the second chamber 26.
i 25 Closing valve 18 allows the pressure in the second chamber 26 to decay due to the flow through the pore 28. The pressure in the second chamber 26 is measured by pressure sensor 34. The rate of change of pressure in the . secont chamber 26 can be determined by comparing signals from the pressure ænsor 34. The sample gas 20 pressure decay rate is then compared ` ~ 30 t~ a pressure decay rate of a standard or reference gas 31.
The reference gas 31 has a known composition and a known density -: Pref at base conditions (i.e. 60F, 14.71 psia). About once a day, valve 18 is closed and valve 29 is opened for a referencing process. The referencing .. .
, ,;.
; ' .
21~1~7~
.
process usually requires a few cyc1es of operation (about 30 to 200 seconds).
In the referencing process, the pressure decay rate of the reference gas 31 is measured in the same manner as the sample gas 20 pressure decay rate.
Also, the temperature of the reference gas 31 in the second chamber 26 is 5 measured by temperature sensor 35. The values of the reference gas 31 temperature and pressure decay rate are stored in the microprocessor 36.
The base condition density Pb of the sample gas 20 can be determined relative to the known base condition density Pref for the reference gas 31 by multiplying by the ratio of sample gas 20 pressure decay rate compared to 10 the reference gas 31 pressure decay rate. Also, changes in temperature from the time of the referencing process can be accounted for by measuring the - temperature of the sample gas 20 in the second chamber 26, comparing the sample gas 20 temperature to the reference gas 31 temperature, and making a linear adjustment.
The above-described referencing process could be replaced by initial absolute calibration. However, frequent referencing eliminates longer term calibration drift.
~ After the base condition density Pb has been determined during a ^' measurement cycle, the valve 16 before the first chamber 22 is closed so that ^~ 20 the flowing condition density pf can be determined. At the point in time ;- when the valve 16 is closed, pressure in the first chamber 22 is initially at pipeline pressure. The pressure in the first chamber 22 is measured using a - pressure sensor 30. The flow of sample gas 20 from the first chamber 22 is ~ controlled by a non-venting mechanical pressure regulator 24 which : ~ 25 reduces the sample gas pressure to a lower level, about 20-30 psig, after the sample gas 20 flows from the first chamber 22. Following the mechanical - ~ regulator 24, sample gas 20 flows into the second chamber 26 since the valve . 18 should be open.
- ~ With the valve 18 open, the pressure applied to the pore 28 (i.e., the ` 30 pressure in the second chamber 26) raises to the level of the regulator 24 output and is maintained at roughly that level by the regulator 24. The sample gas 20 flows through the pore 28 which now operates as a sonic nozzle 28 because of the higher applied pressure. The mechanical regulator 24 is set such that the pressure applied to the pore 28 (i.e sonic nozzle 28) is 210~97~
sufficient to ensure critical or sonic flow during the flowing condition mode.
Volumetric flow through the sonic nozzle 28 is therefore constant throughout the flowing condition mode of each measurement cycle (i.e. the cycle of measured pressure decay in the first chamber 22). This is because the composition x of the sample gas 20 is fixed since closing valve 16 traps sample gas 20 in the first chamber 22. The composition x of the sample gas 20 that flows through the sonic nozzle 28 (i.e. pore 28) can and will vary from measurement cycle to cycle, but in any single measurement cycle, the composition x is fixed (because valve 16 closes), and the sonic nozzle 28 flowrate is therefore constant during that cycle.
As explained later under the Scientific and Mathematical Analysis, critical flow through a sonic nozzle 28 is the preferred manner of ` determining the mass flowrate of sample gas 20 at flowing conditions and at - 15 base conditions that exits the first chamber 22. Other types of flow restrictors, besides a sonic nozzle 28, such as orifices, capillaries, and venturis could be used, however. The term flow restrictor as used herein refers to a fluid device where the flow through the device corresponds to - the applied gas pressure. It should also be noted that the analysis associated ` 20 with Eqs. (11) through (15) is specifically for a sonic nozzle 28.
` As molecules of sample gas 20 flow from the first chamber 22, the pressure in the first chamber 22 reduces and the pressure sensor 30 measures the changing pressure.
The sample gas 20 must remain at pipeline gas temperature Tf when ' 25 the invention operates in flowing condition mode. Referring to Fig. 6, this ~-~ may be done by inserting the unit into the pipeline 14.
In particular, the meter 10 shown in Fig. 6 is built into a meter body 44. The entire unit (i.e. the meter 10 built into the body 44) can then be - installed easily within a pipeline, such as 14.
In Fig. 6, sample gas 20 enters the meter 10 through a sample gas inlet 46 when valve 16 is open. The sample gas inlet 46 is immersed within the pipeline gas 12 flowing through the meter body 44. Since the pipeline gas 12 flowing through the pipeline 14 may be at a high pressure, a diaphragm 47 is used to actuate valve 16. In particular, the diaphragm 47 is located in casing 48 and splits the volume inside the casing 48 into an upper poffion 50 and a ,-:
2~197~
lower portion 52. When the solenoid valve 57 is open, gas 54 at pipeline pressure flows from within the pipeline 14 (or the meter body 44) through line 56, through regulator 59, into the upper portion 50 of the casing 48. The regulator 59 reduces the gas pressure to about 50 psi. The pressure in the 5 upper portion 50 of the casing 48 pushes on the diaphragm 47 and drives down arm 58 to open valve 16. The opening and closing of the solenoid valve 58 and, thus valve 16 through diaphragm 47 and arm 58, is controlled by electrical signals. Cavity 50 has an atmospheric bleed.
Alternatively, the source of pressure in cavity 50 could be air pressure 10 supplied by an external instrument air supply.
Valve 16 is shown in detail in Fig. 7. When valve 16 is closed, a spring 59 pushes a conical shaped plug 61 into an annular channel 63 to ;~ prevent flow through the channel 63. When the solenoid valve 57 is opened and diaphragm 47 drives down arm 58, arm 58 displaces plug 61 - 15 downward against the force of the spring 59 and opens valve 16. With the plug 61 displaced downward, sample gas 20 is able to flow into the valve 16 through the annular channel 63 and exit the valve 16 through valve outlet 65 which is in fluid communication with the annular cylinder 63.
, When valve 16 is open, sample gas 20 flows into the first chamber 22.
In Fig. 6, the first chamber 22 consists mostly of a hollow coil 23 that spiralsaround the inside perimeter of the meter body 44. The sample gas 20 flows from the first chamber 22 through the mechanical pressure regulator 24.
The regulator 24 operates in a normal gas regulator fashion to reduce the sample gas 20 pressure (i.e. a regulator valve (not shown) is positioned by an ; ~; 25 arm 68 that is displaced by a combination of a spring 60 and a diaphragm 62 within a casing 64.
After flowing through the pressure regulator 24, the sample gas 20 ~- flows to the second chamber 26 if valve 18 is open. In most circumstances, the pressure in line 70 following the pressure regulator 24 is sufficiently low (i.e. 20-30 psig) so that valve 18 can be a solenoid valve rather than a diaphragm driven valve like valve 16.
After the sample gas 20 flows through the second chamber 26, it flows through the pore 28. The pore 28 is shown in detail in Fig. 8. It is preferred that pore 28 be a 0.002" diameter hole drilled through a piece of sapphire 72.
Such a pore 28 can be used to accurately and reliably measure base condition 21~197.~
density Pb using the invention described in U. S. Patent No. 4,677,841 and a reference gas apparatus as discussed above. The pore 28 operates to measure the density of sample gas flowing through the pore 28 (as inversely proportional to the square of the mass flow rate of sample gas 20 through 5 the pore 28) when the valve 18 is closed and the meter 10 is operating in the ` base condition mode. When valve 18 is open and the meter 10 is operating in the flowing condition mode, the pore 28 acts as a sonic nozzle 28 because the ratio of the sample gas pressure applied to the pore 28 (i.e. sonic nozzle 28), which is the regulated pressure in line 70 and the second volume 26, `` 10 compared to the pressure downstream of the pore 28 (i.e. sonic nozzle 28), which is ambient pressure, is sufficient so that there is constant sample gas mass flow through the pore 28 (i.e. sonic nozzle 28).
After the sample gas 20 flows through the pore 28, it exhausts through opening 74. The exhaust gas can be catalytically burned or returned 15 to the pipeline 14 if desired.
The volume of the first chamber 22 is the volume along the sample gas flow path from the point 67 in the valve 16 where plug 61 meets channel 63 to the point (not shown) in the pressure regulator 24 where the regulator valve reduces the sample gas 20 pressure. Since intake valve 16, 20 the hollow coil 23, and the regulator 24 are immersed within the pipeline ~', gas stream 12, the temperature of the sample gas 20 in the first chamber 22 is `~ maintained at the flow temperature Tf of the pipeline gas 12. Pressure sensor 30 measures the sample gas pressure in the first chamber 22 and temperature ser 32 measures the temperature in the first chamber 22.
25 Pressure sensor 34 measures sample gas pressure in the second chamber 34.
Signals from the temperature sensor 32 and pressure sensors 30 and 34 are sent to the microprocessor 36 to calculate the volume correction ratio(Pf~
and the flowing condition density pf.
~'' .
2~0197~
-1~
Scientific and Mathematical Analysis The following analysis is recited to further describe the present invention and to describe the preferred nature of the role of the microprocessor 36. Assuming that the molecular weight (or composition x) 5 and volume of the gas is the same at flowing conditions as at base conditions, the real gas law allows the volume correction ratio Pf to be written as:
.
~, p, P, Zb Tb Pb b f I (1) where the subscript "f" refers to flowing conditions, subscript "b" refers to - 10 base conditions, "p" is gas density"~Pn is gas pressure, "Z" issupercompressibility, and "T" is gas temperature. In the present invention, the composition x of the sample gas 20 is the same for conditions at base and for flowing conditions because valve 16 closes and, thus, traps a cerhin population of sample gas 20 molecules in the first chamber 22. Also, the -' 15 base condition volume is the same as the flowing condition volume ~, because the volume of the first chamber 22 remains constant.
Therefore, the volume correction ratio (Pf) can be determined by measuring three ratios pbf, Tb, and zb. The base condition temperature Tb and pressure Pb are assigned (e.g. 60F, 14.7 psia) and need not be measured.
20 The absolute flowing condition temperature Tf and pressure Pf are measured by temperature sensor 32 and pressure sensor 30, respectively.
-~ Sensors 30 and 32 measure conditions within the first chamber 22, but the temperature of the sample gas 20 in the first chamber 22 is mainhined at Tf and the pressure of the sample gas 20 in the first chamber 22 is preferably Pf 25 at the point in time when valve 16 closes. Note that it would be possible to measure Tf and Pf directly from the pipeline 14, although this is not necessary.
21~1979 While both the T and P measurements seem to form ratios, the transducer constants, or slope calibrations do not cancel because the values of both Pb and Tb are fixed and not measuret. Accuracy of pressure sensor 30 and temperature sensor 32 therefore must be maintained.
In order to measure the supercompressibility ratio zb, the time rate of change of pressure in the first chamber 22 and the pressure of the sample gas at the sonic nozzle 28 inlet are measured. The following analysis explains .~ Zb in detail the preferred way of determining the supercompressibility ratio zf.
The time rate of change of pressure in a volume with gas flowing 10 therefrom at a mass flowrate (D is:
P ~RT~ z2 ~
MV) I1-c~,x~ (2) ;
where c(T,x) is a third virial pressure coefficient which is a function of i, temperature T and composition x, and Z is a supercompressibility factor ~;~ which can be represented by the pressure virial æries: Z(T,x) = 1 + b(T,x)P +
15 c(T,x)P2; where b(T,x) is a second virial pressure coefficient and also a function of temperature T and composition x. In the gas density meter 10 as ~;~ shown in Figs. 1 and 6, the volume V appearing in Eq. (2) is the volume of the first chamber 22. Since the first chamber 22 is immersed in the pipeline gas stream 12, the sample gas temperature T in the first volume is the 20 pipeline flow condition temperature Tf and can be considered to be constant for one measurement cycle which is the period between successive closings of valve 16 intended to be about 30 to 60 seconds. Also, the composition x of sample gas 20 in the first chamber 22 is constant over a measurement cycle because of the above described trapping by closing valve 16.
The invention measures the rate of change of pressure as described in - Eq. (2) at several pressures as the pressure in the first chamber 22 decays over a measurement cycle. Referring to Fig. 3, these several pressures include pressures approaching the flowing condition pressure Pf and the base condition pressure Pb.
:, .
~- 21~1~79 Since the quantity (MV/RT) is constant as pressure in the first chamber 22 decays over a single measurement cycle, (i.e. the sample gas 20 ; in the first chamber has constant molecular weight M, temperature T, and is ~ 5 within the same volume V), Eq. (2) can be solved at two different pressures x P1 and P2-Z2,(TpX~) ~2Z22( Pl [ 1 -c(Tyx~ p~2 ] P2[ 1 -c(T,x~ P2] (3) where P represents dt. Note both sides of Eq. (3) are expressed in terms of the flowing temperature Tf. This is because the design condition whereby the sample gas 20 in the first chamber 22 is held at ~e flowing temperature Tf when the meter 10 is operating in flowing condition mode. This ~ apparent error in measurement conditions is accounted for by a '. .;'i Z22(T"x~) te nperature transfer term Z2(Tvx~
Solving Eq. (3) for Zl/z2 yields:
Z~ D2~P~ -C~
Z2tr~ IJ~p;Jll-ccr~ )p2J
:, Multiplying each side of Eq. (4) by the temperature transfer term ¦Z
results in ;`' Z12~ 'Z22tr~ 1)2~1-1-C~ )P121 IPI~
Z2trvX~) IZ2crvx~)J ~ -CCr~X~)P2~ IP2~ (5) .... .
~ Substituting Eq. (5) into Eq. (1) and designating Pl as Pf and P2 as Pb provides Pf ~ the volume correction ratio Pb:
21~79 -1~
p IP ~ ~T ~ ¦Pb~ ~f~ ¦-zbcTvx~ Ctl~rxf)pbll ~`` p ~PbJ ~Tf) ~p J ~(I)bJ ¦zb~Y f)~ 'P f ~ (6) Since the base condition pressure Pb is usually close to 1 Bar and "c(Tf,xf)" is usually small (s3xlO~ Bar2), Eq. (6) can be fur~er modified to:
Pf ~Pf~ ~Tb~ P;~ ¦'C3f~ ¦Zb(TVXf)l' r 2 : ~ 5 Pb ~PbJ ~Tf) PfJ ~tDb~ ¦Zb~pXf)J ¦2-C(~pXf)Pf¦
's ~e temperature transfer term [~1' can be stated in terms of : . the virial series as:
~, Zb (Tb~Xf) 1 +b CI b~Xf)Pb+C CrVXf)P2b Zb~Xf) 1 +b(Tf,X~Pb+cClyxf)pb (8) Since the product cI,x) Pb is very small, Eq. (8) may be written as:
(9) Substituting Eq. (9) into Eq. (7) results in:
Pf (Pf) (Tb) (P") ~ f~ { [b k 2 ]
, ' ".
21~1~7~
Equation (10) now becomes the basis for measuring the volume correction ratio p--b. The last term in Eq. (10) involving the third virial coefficient "c" is nearly unity. The average value of "c" is about 3x10~ Bar-2.
Using this value for "c", the flowing condition pressure Pf can reach a 5 maximurn of 55 Bar and the error associated with approximate1y the last term as unity will be less than 0.5%. Alternatively, an average value of "cn can be adopted as a compensating constant, and the maximum pressure Pf limit can be raised to about 85 Bar for 0.5% accuracy. In addition, as will be shown later, it may be possible to calculate the value "cn to further raise the 10 maximum operating pressure Pf (See Eqs. 22-26).
The term 1 + lb (Tb,x) - b (Tf,x)l Pb in Eq. (10) is also nearly unity for gas flowing at near ambient conditions and can be ignored. For example, the value of ~'bn (the second virial pressure coefficient) is generally about 2x10-3because natural gaæs are composed mostly of methane. The value of ~'bn 15 varies with temperature, but the variation is small about 40 PPM/C for natural gases. Natural gas pipelines usually operate at 40F to 100F (i.e., 4Cto 38C), and Pb is typically about 1 Bar. Thus, for natural gas pipelines, the maximum temperature difference (Tb-Tf) should be about i10C so the maximum expected variation in the term 1 + lb(Tb x) - (Tf,x)] Pb iS about 20 iO.02 t~ 0.04%.
The term ( f)l/2 in Eq. (10) must be determined by either measuring the mass flow rates, ~f for flowing conditions and c~b for base conditions; or holding the rates ~f or ~Db constant. If the mass flowrate through the sonic nozzle 28 (i.e., pore 28) is held constant, the ratio of mass flowrates ( f ) / iS
25 unity. This can happen, however, only if pressure regulator 24 is perfect.
But, most pressure regulators cannot maintain perfect output pressure when inlet pressure varies. In the real world, where mass flowrate through the sonic nozzle 28 is not constant due to real regulator conditions, the mass flowrate ratio ~ f~ / iS measured. This measurement or correction for ~J
30 changes in mass flowrate during a measurement cycle can be done by measuring the sonic nozzle 28 inlet pressure using a pressure transducer 34 210.~97~
-1~
to measure the pressure in the second chamber 26. The following analysis explains this method of measuring the mass flowrate ratio ( f) / .
The sonic velocity of a fluid v can be accurately represented as:
2 ~dP~
! ld p~. (11) ~, where P is absolute pressure, p is density, and the process operates at constant entropy s. For an ideal gas:
~' ( P ) t t r (lv .
:
10 where ~ is the weight density and k is the specific heat ratio. Combining Eqs.
(11) and (12) for an ideal gas gives:
V=~/kgcRT (13) where R is the universal gas constant, T is the absolute temperature and gc `~a' ~ 15 is the gravity constant. The ideal gas equation (i.e., Eq. (13)) is independent ~, of pressure.
~ If supercompressibility exists (i.e., the gas is not ideal), the relation in "~ Eq. (13) is not accurate. But Eq. (13) canbe modified for non-ideal gases to an equa'don of the form:
V=~/k"gcZRT (14) :,, where ks is an isentropic exponent. Note that Z is a function of pressure.
Gasunie's handbook (cite) uses a formula similar to Eq. (14) (i.e.
.~
.
:;~
,:
,`,. , ' `
. . ~ . ' ' ."" ' ' ~ ' ~
, ~'~ ^ ' ' ` ' ' .
' 21~1~7~
.~1000 k Z R T/M ), and the Gasunie handbook indicates that it fits data well.
In Eq. (14), the isentropic constant k9 is a number which forces the P-v relation to accurately approximate the local real gas relation for a small isentropic expansion or compression. The constant k9 is not a thermodynamic property and shoutd not be associated with the specific heat ratio for real gases, although comparing Eq. (13) to Eq. (14) shows that ks = k for ideal gases (i.e. Z=1). At most, ks seems to be inversely proportional molecutar weight M (note the Gasunie represenhtion .~1000kZRT/M )-Experimental testing has confirmed Eq. (14). The testing was carried out on a pore of the sarne type as pore 28 (i.e. sonic nozzle 28) but at higher idet pressures than is preferred in the present invention. Sonic flow through the pore (i.e. sonic nozzle) was, therefore, present. The testing was done for three gases: nitrogen, methane, and ethane. Ethane exhibited the largest Z effects. The actual sonic velocity at the pore inlet for both ethane and methane have been calculated and the results confirm Eq. (14). Refer in particular to Fig. 9 in which normalized sonic velocity through the pore was plotted versus pressure, along with Z and Z1/2 for both methane and ethane. From Fig. 9, it is easy to see that normalized velocity data parallels Zl/2, not Z; and this confirms Eq. (14).
Using Eq. (14), the mass flowrate ON of sample gas 20 exiting the nozzle 28 (i.e., pore 28), under critical pressure conditions, can be stated as:
~I)N = Vl AN Pl = AN Pl Ml (Z RgT ) (15) :
- where the subscript "N" indicates a nozzle mechanical condition, "In - 25 represents a dynamic idet condition, and AN is the cross-sectional area of the sonic nozzle 28. As used herein, the term critical flow means a condition where the ratio of pressure upstream from the sonic nozzle 28 compared to the pressure downstream of the sonic nozzle 28 is sufficient so that the relations in Eqs. (14) and (15) for non-ideal gases are substantially accurate.
2iO1973 -2~
The mass flowrate ratio (~--) for two states, where pressure in the Arst ` chamber 26 is P1 or P2, can be represented as:
~ ~1 Z /2 Pll ; (D2 ZIt P12 (16) 5 Referring to Eqs. (15) and (16), the temperature ratio (~) is nearly constant over a single measurement cycle (due to the heat sinking of the assembly) Tu~ /2 ~` and the ratio Tu can therefore be ignored.
` Since the sonic nozzle inlet pressure is relatively low (~ 30 psig), the supercompressibility Z in Eq. (16) can be represented as a truncated virial 10 pressure series:
~..
Z=1~P and ~=l~P
2 (17) .
Due to the small nature of "bn, the supercompressibility ratio in Eq. (16) can be represented by:
= 1 + 2 ( P~2 -PIl ) ll (18) , Equation (18) shows that the low pressure operating conditions at the ~: ~ inlet of the sonic nozzle 28 (i.e., pore 28) are accompanied by ~`'r~ supercompressibility effects too small to be of concern when determining ~ 20 the mass flowrate ratio (a~l) In Eq. (18), the value of the second virial - -~ pressure coefficient ~b~v is about 0.002 Bar-l and the expected variation of the ~, sonic nozzle inlet pressure is about 1/3 Bar. Thus, the second term in Eq.
., ~
.~
, .
, '~' :
.
. ~
2101~7~
(19) can be expected to vary about 0.03% and can be neglected. This means that the mass flowrate ratio (--) in Eq. (16) can be determined by measuring, only the absolute pressure at the sonic nozzle 28 inlet. The mass flow ratio ) in Eq. 16) can be represented as:
~1 = PLI
~D2 P L 2 (l9) Equation (10) can now be rewAtten as:
) ( ) ~ /2 ( )1/2 Pb Pb Tf ~P ~ PLb 12'C(T~Xf)P,¦ (20) 10 where (p--) is substituted for (~nd the term ~1 +~b~VXf)~bCrvxf)]pb~ is ignored. To solve Eq. (20), the pressure decay rate ratio (p~) must be determined. The following discussion explains the determination of the ;` P /2 pressure decay rate ratio Pf Figures 2 and 3 are the results of a computer simulation. Figure 2 represents a typical plot of pressure in the first chamber 22 as it decays over a measurement cycle. The pressure transducer 30 that senses the pressure in the first chamber 22 should have sufficient bandwidth to accurately measure the pressure in the first chamber 22 throughout the entire measurement ~` cycle. The miaoprocessor 36 records several pressure measurements from . 20 transducer 30 as the pressure decays and the time that each pressure measurement occurs. The time rate of change of the pressure in the first chamber 22 (i.e. P') can then be computed for several pressures during the 2101~79 measurement cycle. Figure 3 shows the variation in P' versus pressure in the first chamber 22. The curvature of P' is due to supercompressibilit,v of the gas.
The value of P' at flowing conditions Pf and at base conditions Pb cannot be measured directly without additional apparatus. To directly measure P' at Pf would require an initial or starting pressure in the first chamber 22 that is higher than the pipeline gas 12 pressure because the P' measurement should be centeret about the pipeline gas 12 pressure Pf. For the same reason, directly measuring P'b would probably require exhausting the first chamber 22 into a partial vacuum. Direct measurement of P'b and P'f does not therefore seem practical.
However, intermediate P' values can be measured accurately.
Numeric methods can then be used to obhin accurate values of P'f and P'b.
For example, considering the two left most dah points in Fig. 3, (PO~P~o) and P1,P'l), the time rate of change of pressure at the base condition is:
p ¦(Pb Pl) Po (Pj ) .¦ (21) A similar expression for P'f also exists. One condition on the accuracy of Eq.
(21) is that the function p be relatively linear in the region local to Pb (or Pf if determining P'f). Figure 3 shows that the rate of change of pressure P' with pressure or time is a slowly varying function, so Eq. (21) is accurate for interpolation.
Measurement of pressure decay rates P' required to solve Eq. (21) for P'b and P'f can be accomplished in several fashions. The preferred method ; 25 is to divide the time between two pressure measurements ~t into the difference of the measured pressures ~P to obtain values of P' and then use Eq. (21) to interpolate the values of P'b and P'f. This method has the - advantage of stability and has good resolution because time measurement resolution is very good.
Another method is to record the absolute pressures in a time series (i.e. record the pressure in the first chamber 22 vs. time as in Fig. 2).
.
~ ' i 2iO1~79 Numeric differentiation methods can then be applied to transform the pressure time series into normalized P' versus P data as in Fig. 3. The numeric differentiation process can be carried out using Newton's formulas for data spaced equally in time, or using LaGrange's formula for data spaced 5 unequally in time. If numeric differentiation is used, it i9 preferred that a large number of points be employed because numeric differentiation can be a noise amplif,ving process.
Analogue differentiation using operational amplifier tecluliques could also be used to measure P', but the relatively slow nature of the 10 measurement cycle would be a severe complication.
~Zb Another method for measuring the supercompressibility ratio ~z in Eq. (1) is to compute the virial coefficients of the gas from pressure data or from the derivative of the pressure data; and, then compute the ratio (z--) using the computed virial coefficients (i.e. Z = 1 + bP + cP2). The value of a 15 single pressure decay rate "n" in a series of pressure decay rate measurements can be described as:
~ ~MV~ Zn Idn Pn ~ RT )n (l-CP2n) (22) where the dependence of Z and c on temperature and composition can be 20 ignored because of the constant composition and relatively constant temperature environment in the first chamber 22. The measured value at another point ~n+1n is:
~MV~ Zn~l ~n~l Pn~l ~ RT Jn~l (l~cPn~l) (23) 25 Since Z can be represented by a virial series as Z = 1 + bP +~p2 Z2/(l-cP2) = 1 + 2b P + (3c+b2)P2. Therefore, in light of Eq. (19), dividing Eq.(22) by (23) results in:
2iO~7~
(~J(P~n~l)-152b(Pn-pn~1)+(3c+b2)(pn-pn~l) (24) because Tn = Tn+l due to the relatively constant temperature in the first chamber 22, Mn = Mn+1 due to trapping (i.e. valve 22 is closed), and Vn =
Vn+1 due to the fixed volume of chamber 22.
With another set of pressure measurements, the following similar ~ relation can be formed:
r~`` (p~ ~ (p~ n~ l = 2b (Pn~ l ~ Pn~2) + (3c + b2) (p2n ~ l - p2 2) (25) , , . ., ~- 10 And, with a third set of pressure measurements, the following relation:
:'`'l nJ 1=2b(P 2-P 3)+(3c+b2)(P2nl2-p2n~3) (26) Therefore, by maWng pressure measurements in the first chamber 22 with f~ pressure sensor 30 (i.e. Pn~ Pn+l, etc.) and in the second chamber 26 with pressure sensor 34 (i.e. Pl, n; Pl, n+l; etc.) at a multiplicity of points throughout the measurement cycle, values for "b" and "c" can be computed.
. Note that it might be statistically advanhgeous to solve for "b~ and "(3c +
b2)~ and then compute ~c". Values Zb and Zf can be computed from a virial expansion using these values of "b't and ~c", and known Pb or measured Pf, ~ ~ respectively. (Note that this method of determining the virial coefficients - 20 may also be used to reduce errors associated with the last term ~2~1',x ,~p2]
in Eq. (10) at high operating pressures Pf.) eration With a Volumetric Flowmeter Referring to Figs. 4 and 5, the invention can be used to determine a base condition volumetric flowrate Qb Of the pipeline gas 12 flowing '`
:'.
2iO1973 ` -25-through the pipeline 14 from a flowing condition volumetric flowrate Qf measured by a volumetric flowmeter (e.g. 38 in Pig. 4 and 40 in Pig. 5):
Q (Pf)Q
P b (27) -., 5In Figs. 4 and 5, flow of sample gas 20 and reference gas is shown by ` solid lines and electrical connections are shown by dashed lines. Referring in particular to Fig. 4, flow of pipeline gas 12 through an orifice plate 38 produces a pressure drop ~Pf and a DP cell 42 measures the pressure drop ~Pf. A signal from the DP cell 42 representing ~Pf is sent to the 10 microprocessor 36 to determine the mass flowrate ~f of the pipeline gas 12 through the pipeline 14. The mass flowrate ~f through the orifice plate 38 (i.e. through the pipeline 14) is described by the orifice equation:
~3f=KCDfYf ~/j;;~ (28) 15 where the values of pf is calculated within the microprocessor 36 from previously gathered data; Yf is the upstream expansion factor that accounts for changes in pressure as the pipeline gas 12 flows through the orifice 38 (e.g. AGA-3 standard); and CDf and K are known constants for the orifice plate 38. The microprocessor 36 can then determine the flowing condition 20 flowrate Qf for a differential pressure meter by:
=~l)f=K CDfYf~ (29) and the base condition flowrate Qb by:
K CDfYf ~
P b P b (30) ~' ,`:
,~' ., -,~ , .
, -2~
Figure 5 shows operation of the invention with a linear volumetric flowmeter, such as a turbine meter, a vortex meter, or any other meter with an output signal linearly proportional to the volumetric flow of the gas 12 through the pipeline. If a turbine meter 40 is used as shown in Fig. 5, the 5 flowing condition volumetric flow rate Qf is given by the calibration equation:
Q = f.
Kt (31) where ft is the frequency of rotation and Kt is a scaling constant. The frequency of rotation ft is detected by a frequency detector 41 located on the 10 outer surface of pipeline 14 near impeller 43 of the turbine meter 40. The signal from the frequency detector 41 is then relayed to the microprocessor, along with the dah used to determine the volume correction ratio Pf. The microprocessor 36 then calculates the base volumetric flowrate by:
Q _P f f b-- K (32) Many modifications and variations of the preferred embodiment that are within the spirit and scope of the invention will be apparent to those with ordinary skill in the art.
.
Back~round of tlle Invention In one aspect, the present invention relates to a method and instrumentation that measures, in real time, the tensity of a gas flowing in a pipeline. In this aspect, the invention can also measure a base condition density of the same gas where the base condition densitv corresponds to a ` 5 density (i.e. mass/volume) of the gas determined as if the gas were at some defined base temperature and pressure condition. In a broader aspect, the Pf present invention determines a ratio p of the measured flowing condition density compared to the base condition density. This ratio can be referred to as a volume correction ratio and can be used to translate a measured pipeline gas volumetric flowrate to a corresponding base condition volumetric flow rate.
Heretofore, accurate measurement of gas flowrates has often been difficult because accurate determination of gas density is important when measuring gas flowrates and gas supercompressibility effects gas density in ways that are not detectable by volumetric flowmeters. For instance, with differerltial pressure flowmeters (e.g., an orifice plate meter), the density term appears directly in the volume flow equation (i.e. ~). Also, with other flowmeters (such as turbine flowmeters), as well as with differential pressure flowmeters, the translation of volumetric flowrate to base conditions (i.e. volume correction ratio) is directly related to the density of the flowing gas.
Yet, at elevated gas pressure, computing gas density from the ideal gas - law is inadequate and can be substantially wrong because the ideal gas law -'~ does not properly account for supercompressibility. If gas properties are well 25 known, the difficulty can be partially alleviated because gas . supercompressibility can be directly computed and accommodated in the - ~ calculations. However, if the gas is a variable mixture of unlike molecules, gas properties necessary to compute supercompressibility are unknown and difficult to measure.
.....
., ..... .
.
~, , . .
~'~':- , . .
:
, 210~979 Gas densities, volume correction ratios, and volumetric flowrates at base conditions are normally calculated using flow computers from contemporaneous measurements of several gas parameters. Generally, for determining volumetric flowrates at base conditions, the pipeline gas volumetric flowrate at pipeline conditions is measured, the gas temperature - and pressure at pipeline conditions are measured, and the composition of the gas is measured. Gas composition is normally measured by a chromatograph. From this measured data, gas supercompressibility at both pipeline and base conditions is calculated and from that the density of the gas at pipeline and base conditions is calculated. When operating pressure is elevated (and supercompressibility effects cannot be ignored) gas supercompressibility is usually estimated from either virial equations of state or from correlations such as NX-19. Using virial coefHcients is severely limited because the virial coefficients are functions of temperature, `` 15 pressure, and composition, are largely unknown, and have significant real time uncertainties. Also, using correlations such as NX-19 is often inaccurate because the correlations can be inaccurate for many compositions.
An alternate method for determining volumetric flowrates at base conditions involves the Gerg Equations. The Gerg Equations estimate - supercompressibility and density from knowledge of the heating value, the density of the gas at base conditions and the percentage of carbon dioxide and nitrogen in the exhaust of burned gas. Such measurements can be - made using PMl's GB 3000 product (Precision Measurements, Inc., Tulsa, . 25 Oklahoma). The Gerg Equations allow more rapid computation of ;~ supercompressibility than the above described composition methods and are, therefore, preferred in applications at normal natural gas pipeline ~` pressures. However, as with composition methods, extenshe use of a flow ; computer is required to solve the Gerg Equations.
` 30 Moreover, each of the measurements recited above in describing both the composition method or the Gerg Equations method introduce the potential for measurement error. The aggregation of such errors can substantially influence accuracy. Because of this, it is common practice to frequently calibrate and mainhin each individual measurement device.
,~., ,, ":,.
., :
, . , :' , .. .
.
-, 210197~
In co-pending patent applications, Serial No. 07/793,753, filed on November 18,1991, and Serial No. 07/787,188, filed on November 4,1991, inventions are disclosed that among other things can determine the volume correction ratio of a pipeline gas from energy type measurements.
These inventions involve the measurement of energy flowrate and energy content of a sample of gas tapped from the pipeline. A base condition volumetric flowrate of the sample gas can then be determined by the ratio of the energy flowrate of the sample gas to the energy content of the sample gas. These inventions also measure the ratio of the mass flowrate of the pipeline gas in the pipeline compared to the mass flowrate of the sample gas hpped from the pipeline. The volume correction ratio (or the base condition volumetric flowrate of the pipeline gas) can be calculated easily from the mass flowrate ratio and the base condition volumetric flowrate of the sample gas.
While it is desirable to measure energy flowrate and energy content in many applications, in applications where energy flowrate or energy content are not required (but where base condition volurnetric flowrates are desirable), the energy measurement equipment provides excessive expenæ.
Further, the inventions in these two co-pending patent applications are directed to monitoring combustible pipeline gas and are thus inappropriate for monitoring flows of non-combustible gas.
,. ~
~.
The object of the present invention is to provide an improved ` 25 method and apparatus for measuring a volume correction ratio and a flowing condition density of a pipeline gas.
Another object of the present invention is to make these measurements without making energy type measuremenh.
~-~ The preænt invention is a method and apparatus that can determine 30 a volume correction ratio for a gas (i.e. a ratio of a flowing condition density compared to a baæ condition density p--), and in another aspect can also determine a flowing condition density of a gas. It can make these ~! determinations accurately in real time, without making several highly .-:, .. .
-`` 21~1979 variable measurements, without making any energy type measurements, and without making excessive numeric calculations or correlations in a flow computer.
The present invention measures the volume correction ratio by measuring the temperature and pressure of sample gas tapped from a pipeline as it undergoes controlled expansion. Before the expansion, the sample gas is trapped to assure that it has the same composition as the pipeline gas. The rate of expansion of the trapped sample gas can be controlled in several ways, but it is preferred that the expansion be controlled by establishing mass flow at sonic conditions through a sonic nozzle. It is important that the expansion of the sample gas occur while the trapped sample gas is at substantially the same temperature as the pipeline gas being monitored. In a preferred embodiment, this is achieved by immersing a portion of the meter within the pipeline.
In particular, sample gas is hpped from the pipeline and flows to a first chamber having a fixed volume. A ffrst valve controls the flow of sample gas into the first chamber. When the first valve is opened, sample gas flows into the first chamber preferably until the pressure in the first charnber equals the pressure in the pipeline.
The first valve is then closed and stops further pipeline gas from flowing into the first chamber. In order to measure the volume correction ~, ratio, the sample gas temperature must be maintained at subshntially the same temperature as the pipeline gas temperature. Sample gas flows from ~ ~ the first chamber and then through a flow restrictor located downstream of ;~ 25 the first chamber. The flow restrictor is preferably a sonic nozzle. The sample gas pressure at the inlet of the flow restrictor is established by a ~- pressure regulator so that gas flows from the first chamber at a substantially constant flowrate. If a sonic nozzle is used, the established pressure is sufficient so that the mass flow through the sonic nozzle is critical (i.e., thesonic nozzle inlet pressure is sufficient so that the mass flowrate through `~- the sonic nozzle is proportional to the density and sonic velity of the gas ` ~- and nozzle cross-sectional area).
~, Flowing sample gas from the first chamber when the first valve is closed causes the pressure in the first chamber to reduce over time. The rate , .
of change of the pressure in the first chamber is measured as the pressure in ~',.
, .
"
~' `
the first chamber decays. The temperature of pipeline gas being monitored is also measured.
The volume correction ratio is determined from the rate of change of pressure in the first chamber, the absolute pressure of the pipeline gas, and 5 the absolute temperature of the pipeline gas. The determination can be done without resorting to estimates and approximations for supercompressibility.
The present invention eliminates the need to estimate or consider the effects of supercompressibility when determining gas density because 10 the pressure rate of change measurements in the first chamber are made and used in density calculations directly.
In another aspect of the invention, the volume correction ratio (i.e.
Pf the ratio p of the flowing condition density pf to the base condition density Pb) is used to determine the aowing condition density pf from a measured 15 base condition density Pb. The preferred method and apparatus for measuring the base condition density Pb are disclosed generally in U.S.
Patent No. 4,677,~41 to Kennedy, et al. It is preferred that the sonic nozzle ofthe present invention for measuring the volume correction ratio also be used as a pore where the square of the mass flow rate through the pore is 20 inversely proportional to the density of gas flowing through the pore when the present invention measures the base condition density Pb.
To determine the base condition density Pb, it is preferred that the .-first valve be maintained open to allow sample gas to aow into the first ~`chamber and therefrom; through the pressure regulator; and then into a 2S second chamber also having a fixed volume. Also preferably, the flow into the second chamber is controlled by a second valve. The flow of sample gas from the second chamber is restricted by the pore. The mass flowrate through the pore is measured and from that mass flowrate measurement and a similar measurement of a reference gas which is made periodically, 30 the base condition density Pb of the sample gas can be determined.
~^In another preferred embodiment of the invention, volumetric aow - ~signals are received from a volumetric flowmeter in the pipeline and the volumetric and mass flowrate in the pipeline at base conditions or pipeline conditions can be determined.
. .
,.
21~1~79 The aspect of the invention that measures the base condition density Pb ~i.e., U.S. Patent No. 4,677,841 in conjunction with a reference gas) has proven to be accurate. Also, using a reference gas avoids the need to force the sample gas temperature and pressure to base conditions at the pore, such 5as 14.71 psia and 60F, in order to make a base condition density Pb measurement. This is because the density of the gas flowing through the pore is measured at ambient conditions and then translated to the base condition density through a prior knowledge of the reference gas base condition.
10The foregoing objects and advantages of the present invention will appear from the following description. In the description, references are made to the accompanying drawings which form a part hereof and in which a preferred embodiment of the present invention is shown by way of illustration. Such embodiment does not necessarily represent the full scope 15of the invention however.
Brief Description of the Drawin~s Fig. 1 is a schematic drawing of the fundamental apparatus of the ; preferred embodiment of the present invention;
Fig. 2 is a plot of the pressure in a first fixed-volume chamber as a 20function of time when the present invention as shown in Fig. 1 is operating;
Fig. 3 is a plot of the normalized rate of change of pressure in the first chamber as a function of absolute pressure in the first chamber when the present invention as shown in Fig. 1 is operating;
~^j 25Fig. 4 is a schematic drawing of the preferred way of using a`' differential pressure volumetric flowmeter with the present invention as shown in Fig.1;
Fig. 5 is a schematic drawing of the preferred way of using a linear volumetric flowmeter with the preænt invention;
~-, 30Fig. 6 is a schematic drawing showing the preferred method of installing the present invention;
Fig. 7 is a detailed diagram of a preferred embodiment of a first valve of the present invention;
,~, .
:' 2101~79 Fig. 8 is a detailed cross-sectional view of a preferred embodiment of a pore/sonic nozzle of the present invention; and Fig. 9 shows plots of experimental data for ethane gas and methane gas that relates to a preferred embodiment of the present invention.
Detailed Description of the Preferred Embodiments The present invention includes many aspects, including the aspect of measuring a volume correction ratio [Pf) without measuring a base condition density Pb- The preferred aspect of the present invention, however, involves the measurement of Pb and is, therefore, discussed first In this description of the preferred embodiments of the invention, the apparatus of the preferred embodiment and its operation is first described. Then an analysis of the invention in scientific and mathematical terms is presented. Finally, operation of the present invention with a volumetric flowmeter is described.
Apparahls of the Preferred Embodiments In Fig. 1, a gas density meter of the present invention is designated generally as 10. The solid lines connecting elements of the gas meter 10 represent flow of sample gas and the dashed lines represent electrical signals.
The gas density meter 10 has two modes of operation: a base condition mode and a flowing condition mode. In the base condition mode, the meter 10 measures a base condition density Pb of a pipeline gas 12 flowing through a pipeline 14 (i.e., the density of the gas at some designated " base temperature, pressure, and composition). When the meter 10 is operating in the base condition mode, valve 16 is open and valve 18 cycles open and closed. In the flowing condition mode, the meter 10 measures a ~ flowing condition density pf of the pipeline gas 12 flowing through the pipeline 14 (i.e., the density of the gas 12 at the temperature and pressure that it is flowing through the pipeline 14). When the meter 10 operates in ~- 30 flowing condition mode, valve 16 is closed and valve 18 is open. Acomplete measurement cycle constitutes operation in both the base condition and the flowing condition modes. Note that in applications " 2~0~79 where only a volume correction ratio ~Pf)is desired and it is not necessary to determine either the base or the flowing condition densities, it is not necessar,v to operate in the base condition mode.
Referring in particular to Fig. 1, the sample gas 20 is tapped from the 5 pipeline 14 and flowed into a first fixed-volume chamber 22 when valve 16 is open. The volume of the first chamber 22 is small, about ~10 cubic centimeters. The sample 20 gas fills the first chamber 22 until the pressure in the first chamber 22 is equal to the pipeline gas 12 pressure in the pipeline14. The sample gas 20 then flows to a mechanical pressure regulator 24 10 where the pressure of the sample gas 20 is reduced; and if valve 18 is open, flows into a second fixed-volume chamber 26. The sample gas 20 then flows from the second chamber 26 through a pore 28.
- When operating in the base condition mode, the flow of sample gas 20 through the pore 28 is controlled using valve 18. The pore 28 is operated 15 in a special operating mode so that the square of the sample gas 20 mass `~ flowrate through the pore 28 is inversely proportional to the density of the ` ~s sample gas 20 at the time it is flowing through the pore 28. This special operating mode was discovered by Kennedy and is described in U.S. Patent No. 4,677,841.
There are many ways to determine the sample gas 20 mass flowrate through the pore 28, but the preferred way is to measure a sample gas 20 ;. pressure decay rate in the second chamber 26 and compare this to a reference gas 31 pressure decay rate, as now explained.
Opening valve 18 allows sample gas 20 into the second chamber 26.
i 25 Closing valve 18 allows the pressure in the second chamber 26 to decay due to the flow through the pore 28. The pressure in the second chamber 26 is measured by pressure sensor 34. The rate of change of pressure in the . secont chamber 26 can be determined by comparing signals from the pressure ænsor 34. The sample gas 20 pressure decay rate is then compared ` ~ 30 t~ a pressure decay rate of a standard or reference gas 31.
The reference gas 31 has a known composition and a known density -: Pref at base conditions (i.e. 60F, 14.71 psia). About once a day, valve 18 is closed and valve 29 is opened for a referencing process. The referencing .. .
, ,;.
; ' .
21~1~7~
.
process usually requires a few cyc1es of operation (about 30 to 200 seconds).
In the referencing process, the pressure decay rate of the reference gas 31 is measured in the same manner as the sample gas 20 pressure decay rate.
Also, the temperature of the reference gas 31 in the second chamber 26 is 5 measured by temperature sensor 35. The values of the reference gas 31 temperature and pressure decay rate are stored in the microprocessor 36.
The base condition density Pb of the sample gas 20 can be determined relative to the known base condition density Pref for the reference gas 31 by multiplying by the ratio of sample gas 20 pressure decay rate compared to 10 the reference gas 31 pressure decay rate. Also, changes in temperature from the time of the referencing process can be accounted for by measuring the - temperature of the sample gas 20 in the second chamber 26, comparing the sample gas 20 temperature to the reference gas 31 temperature, and making a linear adjustment.
The above-described referencing process could be replaced by initial absolute calibration. However, frequent referencing eliminates longer term calibration drift.
~ After the base condition density Pb has been determined during a ^' measurement cycle, the valve 16 before the first chamber 22 is closed so that ^~ 20 the flowing condition density pf can be determined. At the point in time ;- when the valve 16 is closed, pressure in the first chamber 22 is initially at pipeline pressure. The pressure in the first chamber 22 is measured using a - pressure sensor 30. The flow of sample gas 20 from the first chamber 22 is ~ controlled by a non-venting mechanical pressure regulator 24 which : ~ 25 reduces the sample gas pressure to a lower level, about 20-30 psig, after the sample gas 20 flows from the first chamber 22. Following the mechanical - ~ regulator 24, sample gas 20 flows into the second chamber 26 since the valve . 18 should be open.
- ~ With the valve 18 open, the pressure applied to the pore 28 (i.e., the ` 30 pressure in the second chamber 26) raises to the level of the regulator 24 output and is maintained at roughly that level by the regulator 24. The sample gas 20 flows through the pore 28 which now operates as a sonic nozzle 28 because of the higher applied pressure. The mechanical regulator 24 is set such that the pressure applied to the pore 28 (i.e sonic nozzle 28) is 210~97~
sufficient to ensure critical or sonic flow during the flowing condition mode.
Volumetric flow through the sonic nozzle 28 is therefore constant throughout the flowing condition mode of each measurement cycle (i.e. the cycle of measured pressure decay in the first chamber 22). This is because the composition x of the sample gas 20 is fixed since closing valve 16 traps sample gas 20 in the first chamber 22. The composition x of the sample gas 20 that flows through the sonic nozzle 28 (i.e. pore 28) can and will vary from measurement cycle to cycle, but in any single measurement cycle, the composition x is fixed (because valve 16 closes), and the sonic nozzle 28 flowrate is therefore constant during that cycle.
As explained later under the Scientific and Mathematical Analysis, critical flow through a sonic nozzle 28 is the preferred manner of ` determining the mass flowrate of sample gas 20 at flowing conditions and at - 15 base conditions that exits the first chamber 22. Other types of flow restrictors, besides a sonic nozzle 28, such as orifices, capillaries, and venturis could be used, however. The term flow restrictor as used herein refers to a fluid device where the flow through the device corresponds to - the applied gas pressure. It should also be noted that the analysis associated ` 20 with Eqs. (11) through (15) is specifically for a sonic nozzle 28.
` As molecules of sample gas 20 flow from the first chamber 22, the pressure in the first chamber 22 reduces and the pressure sensor 30 measures the changing pressure.
The sample gas 20 must remain at pipeline gas temperature Tf when ' 25 the invention operates in flowing condition mode. Referring to Fig. 6, this ~-~ may be done by inserting the unit into the pipeline 14.
In particular, the meter 10 shown in Fig. 6 is built into a meter body 44. The entire unit (i.e. the meter 10 built into the body 44) can then be - installed easily within a pipeline, such as 14.
In Fig. 6, sample gas 20 enters the meter 10 through a sample gas inlet 46 when valve 16 is open. The sample gas inlet 46 is immersed within the pipeline gas 12 flowing through the meter body 44. Since the pipeline gas 12 flowing through the pipeline 14 may be at a high pressure, a diaphragm 47 is used to actuate valve 16. In particular, the diaphragm 47 is located in casing 48 and splits the volume inside the casing 48 into an upper poffion 50 and a ,-:
2~197~
lower portion 52. When the solenoid valve 57 is open, gas 54 at pipeline pressure flows from within the pipeline 14 (or the meter body 44) through line 56, through regulator 59, into the upper portion 50 of the casing 48. The regulator 59 reduces the gas pressure to about 50 psi. The pressure in the 5 upper portion 50 of the casing 48 pushes on the diaphragm 47 and drives down arm 58 to open valve 16. The opening and closing of the solenoid valve 58 and, thus valve 16 through diaphragm 47 and arm 58, is controlled by electrical signals. Cavity 50 has an atmospheric bleed.
Alternatively, the source of pressure in cavity 50 could be air pressure 10 supplied by an external instrument air supply.
Valve 16 is shown in detail in Fig. 7. When valve 16 is closed, a spring 59 pushes a conical shaped plug 61 into an annular channel 63 to ;~ prevent flow through the channel 63. When the solenoid valve 57 is opened and diaphragm 47 drives down arm 58, arm 58 displaces plug 61 - 15 downward against the force of the spring 59 and opens valve 16. With the plug 61 displaced downward, sample gas 20 is able to flow into the valve 16 through the annular channel 63 and exit the valve 16 through valve outlet 65 which is in fluid communication with the annular cylinder 63.
, When valve 16 is open, sample gas 20 flows into the first chamber 22.
In Fig. 6, the first chamber 22 consists mostly of a hollow coil 23 that spiralsaround the inside perimeter of the meter body 44. The sample gas 20 flows from the first chamber 22 through the mechanical pressure regulator 24.
The regulator 24 operates in a normal gas regulator fashion to reduce the sample gas 20 pressure (i.e. a regulator valve (not shown) is positioned by an ; ~; 25 arm 68 that is displaced by a combination of a spring 60 and a diaphragm 62 within a casing 64.
After flowing through the pressure regulator 24, the sample gas 20 ~- flows to the second chamber 26 if valve 18 is open. In most circumstances, the pressure in line 70 following the pressure regulator 24 is sufficiently low (i.e. 20-30 psig) so that valve 18 can be a solenoid valve rather than a diaphragm driven valve like valve 16.
After the sample gas 20 flows through the second chamber 26, it flows through the pore 28. The pore 28 is shown in detail in Fig. 8. It is preferred that pore 28 be a 0.002" diameter hole drilled through a piece of sapphire 72.
Such a pore 28 can be used to accurately and reliably measure base condition 21~197.~
density Pb using the invention described in U. S. Patent No. 4,677,841 and a reference gas apparatus as discussed above. The pore 28 operates to measure the density of sample gas flowing through the pore 28 (as inversely proportional to the square of the mass flow rate of sample gas 20 through 5 the pore 28) when the valve 18 is closed and the meter 10 is operating in the ` base condition mode. When valve 18 is open and the meter 10 is operating in the flowing condition mode, the pore 28 acts as a sonic nozzle 28 because the ratio of the sample gas pressure applied to the pore 28 (i.e. sonic nozzle 28), which is the regulated pressure in line 70 and the second volume 26, `` 10 compared to the pressure downstream of the pore 28 (i.e. sonic nozzle 28), which is ambient pressure, is sufficient so that there is constant sample gas mass flow through the pore 28 (i.e. sonic nozzle 28).
After the sample gas 20 flows through the pore 28, it exhausts through opening 74. The exhaust gas can be catalytically burned or returned 15 to the pipeline 14 if desired.
The volume of the first chamber 22 is the volume along the sample gas flow path from the point 67 in the valve 16 where plug 61 meets channel 63 to the point (not shown) in the pressure regulator 24 where the regulator valve reduces the sample gas 20 pressure. Since intake valve 16, 20 the hollow coil 23, and the regulator 24 are immersed within the pipeline ~', gas stream 12, the temperature of the sample gas 20 in the first chamber 22 is `~ maintained at the flow temperature Tf of the pipeline gas 12. Pressure sensor 30 measures the sample gas pressure in the first chamber 22 and temperature ser 32 measures the temperature in the first chamber 22.
25 Pressure sensor 34 measures sample gas pressure in the second chamber 34.
Signals from the temperature sensor 32 and pressure sensors 30 and 34 are sent to the microprocessor 36 to calculate the volume correction ratio(Pf~
and the flowing condition density pf.
~'' .
2~0197~
-1~
Scientific and Mathematical Analysis The following analysis is recited to further describe the present invention and to describe the preferred nature of the role of the microprocessor 36. Assuming that the molecular weight (or composition x) 5 and volume of the gas is the same at flowing conditions as at base conditions, the real gas law allows the volume correction ratio Pf to be written as:
.
~, p, P, Zb Tb Pb b f I (1) where the subscript "f" refers to flowing conditions, subscript "b" refers to - 10 base conditions, "p" is gas density"~Pn is gas pressure, "Z" issupercompressibility, and "T" is gas temperature. In the present invention, the composition x of the sample gas 20 is the same for conditions at base and for flowing conditions because valve 16 closes and, thus, traps a cerhin population of sample gas 20 molecules in the first chamber 22. Also, the -' 15 base condition volume is the same as the flowing condition volume ~, because the volume of the first chamber 22 remains constant.
Therefore, the volume correction ratio (Pf) can be determined by measuring three ratios pbf, Tb, and zb. The base condition temperature Tb and pressure Pb are assigned (e.g. 60F, 14.7 psia) and need not be measured.
20 The absolute flowing condition temperature Tf and pressure Pf are measured by temperature sensor 32 and pressure sensor 30, respectively.
-~ Sensors 30 and 32 measure conditions within the first chamber 22, but the temperature of the sample gas 20 in the first chamber 22 is mainhined at Tf and the pressure of the sample gas 20 in the first chamber 22 is preferably Pf 25 at the point in time when valve 16 closes. Note that it would be possible to measure Tf and Pf directly from the pipeline 14, although this is not necessary.
21~1979 While both the T and P measurements seem to form ratios, the transducer constants, or slope calibrations do not cancel because the values of both Pb and Tb are fixed and not measuret. Accuracy of pressure sensor 30 and temperature sensor 32 therefore must be maintained.
In order to measure the supercompressibility ratio zb, the time rate of change of pressure in the first chamber 22 and the pressure of the sample gas at the sonic nozzle 28 inlet are measured. The following analysis explains .~ Zb in detail the preferred way of determining the supercompressibility ratio zf.
The time rate of change of pressure in a volume with gas flowing 10 therefrom at a mass flowrate (D is:
P ~RT~ z2 ~
MV) I1-c~,x~ (2) ;
where c(T,x) is a third virial pressure coefficient which is a function of i, temperature T and composition x, and Z is a supercompressibility factor ~;~ which can be represented by the pressure virial æries: Z(T,x) = 1 + b(T,x)P +
15 c(T,x)P2; where b(T,x) is a second virial pressure coefficient and also a function of temperature T and composition x. In the gas density meter 10 as ~;~ shown in Figs. 1 and 6, the volume V appearing in Eq. (2) is the volume of the first chamber 22. Since the first chamber 22 is immersed in the pipeline gas stream 12, the sample gas temperature T in the first volume is the 20 pipeline flow condition temperature Tf and can be considered to be constant for one measurement cycle which is the period between successive closings of valve 16 intended to be about 30 to 60 seconds. Also, the composition x of sample gas 20 in the first chamber 22 is constant over a measurement cycle because of the above described trapping by closing valve 16.
The invention measures the rate of change of pressure as described in - Eq. (2) at several pressures as the pressure in the first chamber 22 decays over a measurement cycle. Referring to Fig. 3, these several pressures include pressures approaching the flowing condition pressure Pf and the base condition pressure Pb.
:, .
~- 21~1~79 Since the quantity (MV/RT) is constant as pressure in the first chamber 22 decays over a single measurement cycle, (i.e. the sample gas 20 ; in the first chamber has constant molecular weight M, temperature T, and is ~ 5 within the same volume V), Eq. (2) can be solved at two different pressures x P1 and P2-Z2,(TpX~) ~2Z22( Pl [ 1 -c(Tyx~ p~2 ] P2[ 1 -c(T,x~ P2] (3) where P represents dt. Note both sides of Eq. (3) are expressed in terms of the flowing temperature Tf. This is because the design condition whereby the sample gas 20 in the first chamber 22 is held at ~e flowing temperature Tf when the meter 10 is operating in flowing condition mode. This ~ apparent error in measurement conditions is accounted for by a '. .;'i Z22(T"x~) te nperature transfer term Z2(Tvx~
Solving Eq. (3) for Zl/z2 yields:
Z~ D2~P~ -C~
Z2tr~ IJ~p;Jll-ccr~ )p2J
:, Multiplying each side of Eq. (4) by the temperature transfer term ¦Z
results in ;`' Z12~ 'Z22tr~ 1)2~1-1-C~ )P121 IPI~
Z2trvX~) IZ2crvx~)J ~ -CCr~X~)P2~ IP2~ (5) .... .
~ Substituting Eq. (5) into Eq. (1) and designating Pl as Pf and P2 as Pb provides Pf ~ the volume correction ratio Pb:
21~79 -1~
p IP ~ ~T ~ ¦Pb~ ~f~ ¦-zbcTvx~ Ctl~rxf)pbll ~`` p ~PbJ ~Tf) ~p J ~(I)bJ ¦zb~Y f)~ 'P f ~ (6) Since the base condition pressure Pb is usually close to 1 Bar and "c(Tf,xf)" is usually small (s3xlO~ Bar2), Eq. (6) can be fur~er modified to:
Pf ~Pf~ ~Tb~ P;~ ¦'C3f~ ¦Zb(TVXf)l' r 2 : ~ 5 Pb ~PbJ ~Tf) PfJ ~tDb~ ¦Zb~pXf)J ¦2-C(~pXf)Pf¦
's ~e temperature transfer term [~1' can be stated in terms of : . the virial series as:
~, Zb (Tb~Xf) 1 +b CI b~Xf)Pb+C CrVXf)P2b Zb~Xf) 1 +b(Tf,X~Pb+cClyxf)pb (8) Since the product cI,x) Pb is very small, Eq. (8) may be written as:
(9) Substituting Eq. (9) into Eq. (7) results in:
Pf (Pf) (Tb) (P") ~ f~ { [b k 2 ]
, ' ".
21~1~7~
Equation (10) now becomes the basis for measuring the volume correction ratio p--b. The last term in Eq. (10) involving the third virial coefficient "c" is nearly unity. The average value of "c" is about 3x10~ Bar-2.
Using this value for "c", the flowing condition pressure Pf can reach a 5 maximurn of 55 Bar and the error associated with approximate1y the last term as unity will be less than 0.5%. Alternatively, an average value of "cn can be adopted as a compensating constant, and the maximum pressure Pf limit can be raised to about 85 Bar for 0.5% accuracy. In addition, as will be shown later, it may be possible to calculate the value "cn to further raise the 10 maximum operating pressure Pf (See Eqs. 22-26).
The term 1 + lb (Tb,x) - b (Tf,x)l Pb in Eq. (10) is also nearly unity for gas flowing at near ambient conditions and can be ignored. For example, the value of ~'bn (the second virial pressure coefficient) is generally about 2x10-3because natural gaæs are composed mostly of methane. The value of ~'bn 15 varies with temperature, but the variation is small about 40 PPM/C for natural gases. Natural gas pipelines usually operate at 40F to 100F (i.e., 4Cto 38C), and Pb is typically about 1 Bar. Thus, for natural gas pipelines, the maximum temperature difference (Tb-Tf) should be about i10C so the maximum expected variation in the term 1 + lb(Tb x) - (Tf,x)] Pb iS about 20 iO.02 t~ 0.04%.
The term ( f)l/2 in Eq. (10) must be determined by either measuring the mass flow rates, ~f for flowing conditions and c~b for base conditions; or holding the rates ~f or ~Db constant. If the mass flowrate through the sonic nozzle 28 (i.e., pore 28) is held constant, the ratio of mass flowrates ( f ) / iS
25 unity. This can happen, however, only if pressure regulator 24 is perfect.
But, most pressure regulators cannot maintain perfect output pressure when inlet pressure varies. In the real world, where mass flowrate through the sonic nozzle 28 is not constant due to real regulator conditions, the mass flowrate ratio ~ f~ / iS measured. This measurement or correction for ~J
30 changes in mass flowrate during a measurement cycle can be done by measuring the sonic nozzle 28 inlet pressure using a pressure transducer 34 210.~97~
-1~
to measure the pressure in the second chamber 26. The following analysis explains this method of measuring the mass flowrate ratio ( f) / .
The sonic velocity of a fluid v can be accurately represented as:
2 ~dP~
! ld p~. (11) ~, where P is absolute pressure, p is density, and the process operates at constant entropy s. For an ideal gas:
~' ( P ) t t r (lv .
:
10 where ~ is the weight density and k is the specific heat ratio. Combining Eqs.
(11) and (12) for an ideal gas gives:
V=~/kgcRT (13) where R is the universal gas constant, T is the absolute temperature and gc `~a' ~ 15 is the gravity constant. The ideal gas equation (i.e., Eq. (13)) is independent ~, of pressure.
~ If supercompressibility exists (i.e., the gas is not ideal), the relation in "~ Eq. (13) is not accurate. But Eq. (13) canbe modified for non-ideal gases to an equa'don of the form:
V=~/k"gcZRT (14) :,, where ks is an isentropic exponent. Note that Z is a function of pressure.
Gasunie's handbook (cite) uses a formula similar to Eq. (14) (i.e.
.~
.
:;~
,:
,`,. , ' `
. . ~ . ' ' ."" ' ' ~ ' ~
, ~'~ ^ ' ' ` ' ' .
' 21~1~7~
.~1000 k Z R T/M ), and the Gasunie handbook indicates that it fits data well.
In Eq. (14), the isentropic constant k9 is a number which forces the P-v relation to accurately approximate the local real gas relation for a small isentropic expansion or compression. The constant k9 is not a thermodynamic property and shoutd not be associated with the specific heat ratio for real gases, although comparing Eq. (13) to Eq. (14) shows that ks = k for ideal gases (i.e. Z=1). At most, ks seems to be inversely proportional molecutar weight M (note the Gasunie represenhtion .~1000kZRT/M )-Experimental testing has confirmed Eq. (14). The testing was carried out on a pore of the sarne type as pore 28 (i.e. sonic nozzle 28) but at higher idet pressures than is preferred in the present invention. Sonic flow through the pore (i.e. sonic nozzle) was, therefore, present. The testing was done for three gases: nitrogen, methane, and ethane. Ethane exhibited the largest Z effects. The actual sonic velocity at the pore inlet for both ethane and methane have been calculated and the results confirm Eq. (14). Refer in particular to Fig. 9 in which normalized sonic velocity through the pore was plotted versus pressure, along with Z and Z1/2 for both methane and ethane. From Fig. 9, it is easy to see that normalized velocity data parallels Zl/2, not Z; and this confirms Eq. (14).
Using Eq. (14), the mass flowrate ON of sample gas 20 exiting the nozzle 28 (i.e., pore 28), under critical pressure conditions, can be stated as:
~I)N = Vl AN Pl = AN Pl Ml (Z RgT ) (15) :
- where the subscript "N" indicates a nozzle mechanical condition, "In - 25 represents a dynamic idet condition, and AN is the cross-sectional area of the sonic nozzle 28. As used herein, the term critical flow means a condition where the ratio of pressure upstream from the sonic nozzle 28 compared to the pressure downstream of the sonic nozzle 28 is sufficient so that the relations in Eqs. (14) and (15) for non-ideal gases are substantially accurate.
2iO1973 -2~
The mass flowrate ratio (~--) for two states, where pressure in the Arst ` chamber 26 is P1 or P2, can be represented as:
~ ~1 Z /2 Pll ; (D2 ZIt P12 (16) 5 Referring to Eqs. (15) and (16), the temperature ratio (~) is nearly constant over a single measurement cycle (due to the heat sinking of the assembly) Tu~ /2 ~` and the ratio Tu can therefore be ignored.
` Since the sonic nozzle inlet pressure is relatively low (~ 30 psig), the supercompressibility Z in Eq. (16) can be represented as a truncated virial 10 pressure series:
~..
Z=1~P and ~=l~P
2 (17) .
Due to the small nature of "bn, the supercompressibility ratio in Eq. (16) can be represented by:
= 1 + 2 ( P~2 -PIl ) ll (18) , Equation (18) shows that the low pressure operating conditions at the ~: ~ inlet of the sonic nozzle 28 (i.e., pore 28) are accompanied by ~`'r~ supercompressibility effects too small to be of concern when determining ~ 20 the mass flowrate ratio (a~l) In Eq. (18), the value of the second virial - -~ pressure coefficient ~b~v is about 0.002 Bar-l and the expected variation of the ~, sonic nozzle inlet pressure is about 1/3 Bar. Thus, the second term in Eq.
., ~
.~
, .
, '~' :
.
. ~
2101~7~
(19) can be expected to vary about 0.03% and can be neglected. This means that the mass flowrate ratio (--) in Eq. (16) can be determined by measuring, only the absolute pressure at the sonic nozzle 28 inlet. The mass flow ratio ) in Eq. 16) can be represented as:
~1 = PLI
~D2 P L 2 (l9) Equation (10) can now be rewAtten as:
) ( ) ~ /2 ( )1/2 Pb Pb Tf ~P ~ PLb 12'C(T~Xf)P,¦ (20) 10 where (p--) is substituted for (~nd the term ~1 +~b~VXf)~bCrvxf)]pb~ is ignored. To solve Eq. (20), the pressure decay rate ratio (p~) must be determined. The following discussion explains the determination of the ;` P /2 pressure decay rate ratio Pf Figures 2 and 3 are the results of a computer simulation. Figure 2 represents a typical plot of pressure in the first chamber 22 as it decays over a measurement cycle. The pressure transducer 30 that senses the pressure in the first chamber 22 should have sufficient bandwidth to accurately measure the pressure in the first chamber 22 throughout the entire measurement ~` cycle. The miaoprocessor 36 records several pressure measurements from . 20 transducer 30 as the pressure decays and the time that each pressure measurement occurs. The time rate of change of the pressure in the first chamber 22 (i.e. P') can then be computed for several pressures during the 2101~79 measurement cycle. Figure 3 shows the variation in P' versus pressure in the first chamber 22. The curvature of P' is due to supercompressibilit,v of the gas.
The value of P' at flowing conditions Pf and at base conditions Pb cannot be measured directly without additional apparatus. To directly measure P' at Pf would require an initial or starting pressure in the first chamber 22 that is higher than the pipeline gas 12 pressure because the P' measurement should be centeret about the pipeline gas 12 pressure Pf. For the same reason, directly measuring P'b would probably require exhausting the first chamber 22 into a partial vacuum. Direct measurement of P'b and P'f does not therefore seem practical.
However, intermediate P' values can be measured accurately.
Numeric methods can then be used to obhin accurate values of P'f and P'b.
For example, considering the two left most dah points in Fig. 3, (PO~P~o) and P1,P'l), the time rate of change of pressure at the base condition is:
p ¦(Pb Pl) Po (Pj ) .¦ (21) A similar expression for P'f also exists. One condition on the accuracy of Eq.
(21) is that the function p be relatively linear in the region local to Pb (or Pf if determining P'f). Figure 3 shows that the rate of change of pressure P' with pressure or time is a slowly varying function, so Eq. (21) is accurate for interpolation.
Measurement of pressure decay rates P' required to solve Eq. (21) for P'b and P'f can be accomplished in several fashions. The preferred method ; 25 is to divide the time between two pressure measurements ~t into the difference of the measured pressures ~P to obtain values of P' and then use Eq. (21) to interpolate the values of P'b and P'f. This method has the - advantage of stability and has good resolution because time measurement resolution is very good.
Another method is to record the absolute pressures in a time series (i.e. record the pressure in the first chamber 22 vs. time as in Fig. 2).
.
~ ' i 2iO1~79 Numeric differentiation methods can then be applied to transform the pressure time series into normalized P' versus P data as in Fig. 3. The numeric differentiation process can be carried out using Newton's formulas for data spaced equally in time, or using LaGrange's formula for data spaced 5 unequally in time. If numeric differentiation is used, it i9 preferred that a large number of points be employed because numeric differentiation can be a noise amplif,ving process.
Analogue differentiation using operational amplifier tecluliques could also be used to measure P', but the relatively slow nature of the 10 measurement cycle would be a severe complication.
~Zb Another method for measuring the supercompressibility ratio ~z in Eq. (1) is to compute the virial coefficients of the gas from pressure data or from the derivative of the pressure data; and, then compute the ratio (z--) using the computed virial coefficients (i.e. Z = 1 + bP + cP2). The value of a 15 single pressure decay rate "n" in a series of pressure decay rate measurements can be described as:
~ ~MV~ Zn Idn Pn ~ RT )n (l-CP2n) (22) where the dependence of Z and c on temperature and composition can be 20 ignored because of the constant composition and relatively constant temperature environment in the first chamber 22. The measured value at another point ~n+1n is:
~MV~ Zn~l ~n~l Pn~l ~ RT Jn~l (l~cPn~l) (23) 25 Since Z can be represented by a virial series as Z = 1 + bP +~p2 Z2/(l-cP2) = 1 + 2b P + (3c+b2)P2. Therefore, in light of Eq. (19), dividing Eq.(22) by (23) results in:
2iO~7~
(~J(P~n~l)-152b(Pn-pn~1)+(3c+b2)(pn-pn~l) (24) because Tn = Tn+l due to the relatively constant temperature in the first chamber 22, Mn = Mn+1 due to trapping (i.e. valve 22 is closed), and Vn =
Vn+1 due to the fixed volume of chamber 22.
With another set of pressure measurements, the following similar ~ relation can be formed:
r~`` (p~ ~ (p~ n~ l = 2b (Pn~ l ~ Pn~2) + (3c + b2) (p2n ~ l - p2 2) (25) , , . ., ~- 10 And, with a third set of pressure measurements, the following relation:
:'`'l nJ 1=2b(P 2-P 3)+(3c+b2)(P2nl2-p2n~3) (26) Therefore, by maWng pressure measurements in the first chamber 22 with f~ pressure sensor 30 (i.e. Pn~ Pn+l, etc.) and in the second chamber 26 with pressure sensor 34 (i.e. Pl, n; Pl, n+l; etc.) at a multiplicity of points throughout the measurement cycle, values for "b" and "c" can be computed.
. Note that it might be statistically advanhgeous to solve for "b~ and "(3c +
b2)~ and then compute ~c". Values Zb and Zf can be computed from a virial expansion using these values of "b't and ~c", and known Pb or measured Pf, ~ ~ respectively. (Note that this method of determining the virial coefficients - 20 may also be used to reduce errors associated with the last term ~2~1',x ,~p2]
in Eq. (10) at high operating pressures Pf.) eration With a Volumetric Flowmeter Referring to Figs. 4 and 5, the invention can be used to determine a base condition volumetric flowrate Qb Of the pipeline gas 12 flowing '`
:'.
2iO1973 ` -25-through the pipeline 14 from a flowing condition volumetric flowrate Qf measured by a volumetric flowmeter (e.g. 38 in Pig. 4 and 40 in Pig. 5):
Q (Pf)Q
P b (27) -., 5In Figs. 4 and 5, flow of sample gas 20 and reference gas is shown by ` solid lines and electrical connections are shown by dashed lines. Referring in particular to Fig. 4, flow of pipeline gas 12 through an orifice plate 38 produces a pressure drop ~Pf and a DP cell 42 measures the pressure drop ~Pf. A signal from the DP cell 42 representing ~Pf is sent to the 10 microprocessor 36 to determine the mass flowrate ~f of the pipeline gas 12 through the pipeline 14. The mass flowrate ~f through the orifice plate 38 (i.e. through the pipeline 14) is described by the orifice equation:
~3f=KCDfYf ~/j;;~ (28) 15 where the values of pf is calculated within the microprocessor 36 from previously gathered data; Yf is the upstream expansion factor that accounts for changes in pressure as the pipeline gas 12 flows through the orifice 38 (e.g. AGA-3 standard); and CDf and K are known constants for the orifice plate 38. The microprocessor 36 can then determine the flowing condition 20 flowrate Qf for a differential pressure meter by:
=~l)f=K CDfYf~ (29) and the base condition flowrate Qb by:
K CDfYf ~
P b P b (30) ~' ,`:
,~' ., -,~ , .
, -2~
Figure 5 shows operation of the invention with a linear volumetric flowmeter, such as a turbine meter, a vortex meter, or any other meter with an output signal linearly proportional to the volumetric flow of the gas 12 through the pipeline. If a turbine meter 40 is used as shown in Fig. 5, the 5 flowing condition volumetric flow rate Qf is given by the calibration equation:
Q = f.
Kt (31) where ft is the frequency of rotation and Kt is a scaling constant. The frequency of rotation ft is detected by a frequency detector 41 located on the 10 outer surface of pipeline 14 near impeller 43 of the turbine meter 40. The signal from the frequency detector 41 is then relayed to the microprocessor, along with the dah used to determine the volume correction ratio Pf. The microprocessor 36 then calculates the base volumetric flowrate by:
Q _P f f b-- K (32) Many modifications and variations of the preferred embodiment that are within the spirit and scope of the invention will be apparent to those with ordinary skill in the art.
.
Claims (10)
1. Apparatus for correction of volumetric gas flows, the apparatus being characterized by:
means (46) for tapping sample gas (20) from a pipeline (14) through which pipeline gas (12) flows at a flowing temperature and pressure in contrast to a base temperature and pressure;
a first chamber (22) of fixed volume for receiving sample gas (20) tapped from the pipeline (14), the sample gas (20) being maintained at substantially the same temperature as the pipeline gas (12) flowing temperature when contained in the first chamber;
a first pressure sensor (30) for measuring the pressure of the sample gas (20) in the first chamber (22);
a first temperature sensor (32) for measuring the pipeline gas flowing temperature;
a first valve (16) for controlling the flow of sample gas (20) to the first chamber (22);
a flow restrictor (28) located downstream of the first chamber (22) through which sample gas flowing from the first chamber flows;
a pressure regulator (24) located downstream of the first chamber (22) and upstream of the flow restrictor (28) for reducing the pressure of sample gas flowing from the first chamber (22) before it flows to the flow restrictor (28) so that flow from the first chamber (22) through the flow restrictor (28) is maintained at a substantially constant rate; and a controller (36) for receiving signals from the first pressure sensor (30) and the temperature sensor (32) and for determining a ratio of a density of the pipeline gas at the flowing temperature and pressure compared to a density of the pipeline gas at the base temperature and pressure.
means (46) for tapping sample gas (20) from a pipeline (14) through which pipeline gas (12) flows at a flowing temperature and pressure in contrast to a base temperature and pressure;
a first chamber (22) of fixed volume for receiving sample gas (20) tapped from the pipeline (14), the sample gas (20) being maintained at substantially the same temperature as the pipeline gas (12) flowing temperature when contained in the first chamber;
a first pressure sensor (30) for measuring the pressure of the sample gas (20) in the first chamber (22);
a first temperature sensor (32) for measuring the pipeline gas flowing temperature;
a first valve (16) for controlling the flow of sample gas (20) to the first chamber (22);
a flow restrictor (28) located downstream of the first chamber (22) through which sample gas flowing from the first chamber flows;
a pressure regulator (24) located downstream of the first chamber (22) and upstream of the flow restrictor (28) for reducing the pressure of sample gas flowing from the first chamber (22) before it flows to the flow restrictor (28) so that flow from the first chamber (22) through the flow restrictor (28) is maintained at a substantially constant rate; and a controller (36) for receiving signals from the first pressure sensor (30) and the temperature sensor (32) and for determining a ratio of a density of the pipeline gas at the flowing temperature and pressure compared to a density of the pipeline gas at the base temperature and pressure.
2. The apparatus of claim 1, further characterized in that the flow restrictor (28) is a sonic nozzle having an inlet and further characterized in that the pressure regulator (24) reduces the pressure of the sample gas flowing from the first chamber to a pressure that is sufficient to maintain critical flow through the sonic nozzle.
3. The apparatus of claim 2 further characterized by:
a second pressure sensor (34) for measuring the pressure of the sample gas at the sonic nozzle inlet; and further characterized in that the controller (36) receives a signal from the second pressure sensor (34).
a second pressure sensor (34) for measuring the pressure of the sample gas at the sonic nozzle inlet; and further characterized in that the controller (36) receives a signal from the second pressure sensor (34).
4. The apparatus of claims 1, 2, or 3, further characterized in that the pressure regulator (24) reduces the pressure of the sample gas flowing from the first chamber (22) to a pressure that is sufficient to maintain critical flow through the flow restrictor (28).
5. The apparatus of claims 1, 2, 3 or 4, further characterized by:
a volumetric flowmeter (38, 40) for measuring a flowing condition volumetric flow rate of the pipeline gas Qf flowing through the pipeline (14) at a flowing temperature and pressure; and the controller (36) being further characterized by deriving the base condition volumetric flow rate of the pipeline gas Qb by multiplying the flowing condition volumetric flow rate of the pipeline gas Qf by the ratio .
a volumetric flowmeter (38, 40) for measuring a flowing condition volumetric flow rate of the pipeline gas Qf flowing through the pipeline (14) at a flowing temperature and pressure; and the controller (36) being further characterized by deriving the base condition volumetric flow rate of the pipeline gas Qb by multiplying the flowing condition volumetric flow rate of the pipeline gas Qf by the ratio .
6. A method for determining a ratio of a density Pf of a pipeline gas (12) in a pipeline (14) at a flowing temperature and pressure compared to a density Pb of the pipeline gas at a base temperature and pressure, the method being characterized by:
flowing sample gas (20) tapped from the pipeline to a first chamber (22) of fixed volume, the sample gas (20) being maintained at substantially the same temperature as the pipeline gas flowing temperature when contained in the first chamber;
measuring the pressure of the sample gas (20) in the first chamber (22);
measuring the temperature of the sample gas (20) in the first chamber (22) which is substantially the same as the pipeline gas flowing temperature;
stopping the flow of sample gas (20) to the first chamber (22) when the pressure in the first chamber (22) reaches the pipeline gas flow pressure;
flowing sample gas (20) from the first chamber (22) at a substantially constant mass flowrate;
determining the time rate of change of the pressure in the first chamber (22) and deriving the ratio from the time rate of change of pressure in the first chamber (22) and the flowing temperature and pressure.
flowing sample gas (20) tapped from the pipeline to a first chamber (22) of fixed volume, the sample gas (20) being maintained at substantially the same temperature as the pipeline gas flowing temperature when contained in the first chamber;
measuring the pressure of the sample gas (20) in the first chamber (22);
measuring the temperature of the sample gas (20) in the first chamber (22) which is substantially the same as the pipeline gas flowing temperature;
stopping the flow of sample gas (20) to the first chamber (22) when the pressure in the first chamber (22) reaches the pipeline gas flow pressure;
flowing sample gas (20) from the first chamber (22) at a substantially constant mass flowrate;
determining the time rate of change of the pressure in the first chamber (22) and deriving the ratio from the time rate of change of pressure in the first chamber (22) and the flowing temperature and pressure.
7. The method of claim 6, further characterized in that sample gas (20) is flowed from the first chamber (22) at a substantially constant mass flowrate by flowing sample gas (20) from the first chamber (22) to a sonic nozzle (28) at a sample gas pressure sufficient to maintain critical flow through the sonic nozzle.
8. The method of claims 6 or 7, further characterized by the steps of:
measuring the pressure of the sample gas at the inlet to the sonic nozzle (28); and correcting for variations in the sample gas mass flowrate from the first chamber (22).
measuring the pressure of the sample gas at the inlet to the sonic nozzle (28); and correcting for variations in the sample gas mass flowrate from the first chamber (22).
9. The method of claims 6, 7 or 8, further characterized by the steps of:
measuring a base condition density Pb of the pipeline gas; and deriving the flowing condition density Pf of the pipeline gas from the flowing temperature and pressure, the time rate of change of pressure in the first chamber (22), and the base condition density Pb.
measuring a base condition density Pb of the pipeline gas; and deriving the flowing condition density Pf of the pipeline gas from the flowing temperature and pressure, the time rate of change of pressure in the first chamber (22), and the base condition density Pb.
10. The method of claim 9, further characterized in that the base condition density Pb of the pipeline gas is measured by:
flowing sample gas to the first chamber (22) and allowing the sample gas to flow through into a second chamber (26);
stopping the flow of sample gas into the second chamber (26);
flowing sample gas from the second chamber (26) through a pore (28), the pore (28) being sized so that the square of the sample gas mass flow rate through the pore (28) is inversely proportional to the density of the sample gas flowing through the pore (28);
measuring the time rate of change of the pressure in the second chamber (26) as sample gas flows through the pore (28);
deriving a pore density of the sample gas as proportional to the inverse of the square of the time rate of change of pressure in the second chamber (26); and comparing the pore density to the density of a reference gas measured in a like manner.
flowing sample gas to the first chamber (22) and allowing the sample gas to flow through into a second chamber (26);
stopping the flow of sample gas into the second chamber (26);
flowing sample gas from the second chamber (26) through a pore (28), the pore (28) being sized so that the square of the sample gas mass flow rate through the pore (28) is inversely proportional to the density of the sample gas flowing through the pore (28);
measuring the time rate of change of the pressure in the second chamber (26) as sample gas flows through the pore (28);
deriving a pore density of the sample gas as proportional to the inverse of the square of the time rate of change of pressure in the second chamber (26); and comparing the pore density to the density of a reference gas measured in a like manner.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US07/956,143 US5307668A (en) | 1992-10-05 | 1992-10-05 | Gas density meter and method |
| US07/956,143 | 1992-10-05 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CA2101979A1 true CA2101979A1 (en) | 1994-04-06 |
Family
ID=25497811
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002101979A Abandoned CA2101979A1 (en) | 1992-10-05 | 1993-08-05 | Gas density meter and method |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5307668A (en) |
| EP (1) | EP0591639A3 (en) |
| CA (1) | CA2101979A1 (en) |
| MX (1) | MX9305480A (en) |
Families Citing this family (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5540077A (en) * | 1994-06-10 | 1996-07-30 | Scott Specialty Gases, Inc. | Method and gas mixture for calibrating an analyzer |
| US5551282A (en) * | 1995-03-20 | 1996-09-03 | Badger Meter, Inc. | Method and apparatus for measuring volume correction using molar quantities |
| US5684245A (en) * | 1995-11-17 | 1997-11-04 | Mks Instruments, Inc. | Apparatus for mass flow measurement of a gas |
| US5902927A (en) | 1995-12-01 | 1999-05-11 | Perception Incorporated | Fluid metering apparatus and method |
| US6076542A (en) * | 1995-12-01 | 2000-06-20 | Perception Incorporated | Fluid metering method |
| US5734093A (en) * | 1995-12-19 | 1998-03-31 | The Dow Chemical Company | Method and apparatus for determining physical properties of a gas for use in rheometry |
| US5656784A (en) * | 1997-01-03 | 1997-08-12 | Eagle Research Corp. | Fluid density variation compensation for fluid flow volume measurement |
| FR2765966B1 (en) * | 1997-07-09 | 1999-08-20 | Gaz De France | METHOD AND DEVICE FOR DETERMINING THE DENSITY OF A GAS |
| DE59710962D1 (en) * | 1997-08-15 | 2003-12-11 | Alstom Switzerland Ltd | Pipe system for the controlled distribution of a flowing medium and method for operating such a pipe system |
| WO1999036768A1 (en) * | 1998-01-16 | 1999-07-22 | Bg Intellectual Property Limited | Method and apparatus for measuring the relative density of a gas |
| US6058761A (en) * | 1998-01-30 | 2000-05-09 | Badger Meter, Inc. | Measurement of relative density of combustible gases |
| US6119710A (en) | 1999-05-26 | 2000-09-19 | Cyber Instrument Technologies Llc | Method for wide range gas flow system with real time flow measurement and correction |
| FI107840B (en) * | 1999-12-09 | 2001-10-15 | Vaisala Oyj | Measurement method and system for humidity or gas sensor |
| GB0227109D0 (en) * | 2002-11-20 | 2002-12-24 | Air Prod & Chem | Volume flow controller |
| US8949233B2 (en) * | 2008-04-28 | 2015-02-03 | Alexandria Investment Research and Technology, Inc. | Adaptive knowledge platform |
| PT2458357E (en) * | 2010-11-29 | 2014-06-11 | Air Prod & Chem | Method of, and apparatus for, measuring the pressure of a gas |
| ES2749877T3 (en) | 2010-11-29 | 2020-03-24 | Air Prod & Chem | Method and apparatus for measuring the molecular weight of a gas |
| PL2458348T3 (en) | 2010-11-29 | 2014-01-31 | Air Prod & Chem | Method of, and apparatus for, measuring the mass flow rate of a gas |
| ES2659146T3 (en) | 2012-05-24 | 2018-03-14 | Air Products And Chemicals, Inc. | Method and apparatus for providing a gas mixture |
| ES2845173T3 (en) | 2012-05-24 | 2021-07-26 | Air Prod & Chem | Procedure and apparatus for regulating the mass flow rate of a gas |
| PL2667277T3 (en) | 2012-05-24 | 2018-05-30 | Air Products And Chemicals, Inc. | Method of, and apparatus for, providing a gas mixture |
| ES2556783T3 (en) | 2012-05-24 | 2016-01-20 | Air Products And Chemicals, Inc. | Method and apparatus for measuring the physical properties of biphasic fluids |
| EP2667176B1 (en) | 2012-05-24 | 2015-02-25 | Air Products And Chemicals, Inc. | Apparatus for measuring the true contents of a cylinder of gas under pressure |
| ES2905838T3 (en) | 2012-05-24 | 2022-04-12 | Air Prod & Chem | Method and apparatus for measuring the mass flow rate of a gas |
| WO2014151475A1 (en) | 2013-03-15 | 2014-09-25 | Watkins Bobby G Ii | Flow control and gas metering process |
| EP2806271B1 (en) * | 2013-05-24 | 2017-04-26 | Mems Ag | Method and measuring device for determining physical gas properties |
| CN105717159A (en) * | 2016-04-18 | 2016-06-29 | 广东海洋大学 | Method and device for measuring specific heat ratio of gas |
| CN108572014B (en) * | 2017-12-31 | 2023-03-17 | 深圳市前海海洋仪表科技有限公司 | Ultrasonic water meter and method for acquiring water flow temperature by means of ultrasonic water meter |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3701280A (en) * | 1970-03-18 | 1972-10-31 | Daniel Ind Inc | Method and apparatus for determining the supercompressibility factor of natural gas |
| US4285245A (en) * | 1979-12-06 | 1981-08-25 | Precision Machine Products, Inc. | Method and apparatus for measuring and controlling volumetric flow rate of gases in a line |
| US4379402A (en) * | 1981-01-22 | 1983-04-12 | Beckman Instruments, Inc. | Gas analysis instrument having flow rate compensation |
| US4527418A (en) * | 1984-02-09 | 1985-07-09 | Honeywell Inc. | Method of measuring specific gravity and apparatus utilizing the same |
| US4677841A (en) * | 1984-04-05 | 1987-07-07 | Precision Measurement, Inc. | Method and apparatus for measuring the relative density of gases |
| GB8719105D0 (en) * | 1987-08-12 | 1987-09-16 | Schlumberger Electronics Uk | Fluid metering |
| US5092159A (en) * | 1990-10-15 | 1992-03-03 | Texaco Inc. | Method and apparatus utilizing a single nozzle for effecting measurment of steam characteristics |
| US5226728A (en) * | 1991-11-04 | 1993-07-13 | Badger Meter, Inc. | Method and apparatus for measuring mass flow and energy content using a differential pressure meter |
| US5201581A (en) * | 1991-11-18 | 1993-04-13 | Badger Meter, Inc. | Method and apparatus for measuring mass flow and energy content using a linear flow meter |
-
1992
- 1992-10-05 US US07/956,143 patent/US5307668A/en not_active Expired - Fee Related
-
1993
- 1993-07-31 EP EP93112322A patent/EP0591639A3/en not_active Ceased
- 1993-08-05 CA CA002101979A patent/CA2101979A1/en not_active Abandoned
- 1993-09-07 MX MX9305480A patent/MX9305480A/en not_active IP Right Cessation
Also Published As
| Publication number | Publication date |
|---|---|
| MX9305480A (en) | 1994-04-29 |
| US5307668A (en) | 1994-05-03 |
| EP0591639A2 (en) | 1994-04-13 |
| EP0591639A3 (en) | 1995-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2101979A1 (en) | Gas density meter and method | |
| US5201581A (en) | Method and apparatus for measuring mass flow and energy content using a linear flow meter | |
| US5323657A (en) | Volumetric flow corrector and method | |
| US5401162A (en) | Microbridge-based combustion control | |
| US5226728A (en) | Method and apparatus for measuring mass flow and energy content using a differential pressure meter | |
| CN100419385C (en) | Differential pressure type flowmeter and differential pressure type flowmeter controller | |
| CN1688948B (en) | Methods and apparatus for pressure compensation in a mass flow controller | |
| US4767259A (en) | Cooling air flow controlling apparatus for gas turbine | |
| EP0554095A2 (en) | Determination of fuel characteristics | |
| JP2008524585A (en) | System and method for measuring flow rate | |
| US5016482A (en) | Method of taking a proportional sample of flowing gas in a line | |
| US5115687A (en) | Method and apparatus for taking a proportional sample of flowing gas in a line | |
| RU2206075C2 (en) | Means for measuring power consumption | |
| US5551282A (en) | Method and apparatus for measuring volume correction using molar quantities | |
| US5357809A (en) | Volumetric flow corrector having a densitometer | |
| US3653399A (en) | Gas flow controlling system | |
| Aumanand et al. | A novel method of using a control valve for measurement and control of flow | |
| EP0608736B1 (en) | Volumetric flow corrector having a densitometer | |
| EP0069173B1 (en) | Method and apparatus for measuring and controlling volumetric flow rate of gases in a line | |
| US3365932A (en) | Densitometer | |
| JPS61223517A (en) | Fluidic flowmeter | |
| Miralles | Preliminary considerations in the use of industrial sonic nozzles | |
| GB2069726A (en) | Fluid flow control apparatus and method | |
| Levine et al. | A primary standard high vacuum calibration station for industrial applications | |
| Steven et al. | Meeting the Challenges of CO2 Measurement with a new kind of Orifice Meter |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FZDE | Dead |