Method for calibrating a differential pressure fluid flow measuring system

(57) Abrégé/Abstract:

In a flow measurement system (11) comprising a primary flow sensor (10) and a secondary electronic conversion apparatus for converting high and low fluid pressure signals from the primary flow sensor to an electronic signal output representative of fluid flow, a process for calibrating the system as an integral unit, comprising the steps of applying a plurality of known reference fluid flow rates to the primary flow sensor over the system's operative range of flow conditions, sensing the electrical output of the secondary electronic conversion apparatus for each of the said plurality of reference fluid flow rates, establishing a correction factor for each of said values of electrical output which will tend to linearize the relationship between fluid flow rate and the electrical output of the measuring system, and storing the correction factors in non-volatile memory in a microprocessor which is part of the secondary electronic conversion apparatus.

Canada

<img file="CA2247407C_D0001.tif" />

http://opic.gc.ca Ottawa-Hull K1A 0C9 · http://cipo.gc.ca

OPIC-CIPO 191

CA 02247407 1998-08-27

<img file="CA2247407C_D0002.tif" />

PP”T WORLD INTELLECTUAL PROPERTY ORGANIZATION ·*· International Bureau

INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT)

<td> (51) International Patent Classification 6 _; G01F 25/00</td><td> A3</td><td colspan="2"> (11) International Publication Number: WO ,97/46851 (43) International Publication Date: 11 December 1997 (11.12.97)</td>

<td colspan="3"> (21) International Application Number: PÇT/US97/08193 (22) International Filing Date: 15.May 1997 (15.05.97) (30) Priority Data: 08/650,140 17 May 1996 (17.05.96) US (71) Applicant: DIETERICH TECHNOLOGY HOLDING CORP. [US/US]; 5601 North 71st Street, P.O. Box 9000, Boulder, CO 80301 (US). (72) Inventors: SHANAHAN, Eugene, ’Miles; 8417 Saw Tooth Lane, Niwot, CO 80503 (US). GARNETT, John, Everett; 858 Gapter Road, Boulder, CO 80303 (US). GRAY, Robert, James; 941 Eldorado Lane, Louisville, CO 80027 (US). DRAGNEA, Gabriel; 1248 Atwood Court, Longmont, CO 80501 (US). (74) Agent: HANES, Richard, W.; Hanes &amp; Schute, P.C., 7222 Commerce Center Drive #243, Colorado Springs, CO 80919 (US).</td><td> (81) Designated States: AU, BR, CA, CN, JP, KP, MX, RU, SG. Published With international search report. (88) Date of publication of the international search report: 12 February 1998 (12.02.98)</td>

<td colspan="4"> (54) Title: METHOD FOR CALIBRATING A FLOW MEASURING SYSTEM (57) Abstract In a flow measurement system (11) comprising a primary flow sensor (10) and a secondary electronic conversion apparatus for converting high and low fluid pressure signals from the primary flow sensor to an electronic signal output representative of fluid flow, a process for calibrating the system as an integral unit, comprising the steps of applying a plurality of known reference fluid flow rates to the primary flow sensor over the system's operative range of flow conditions, sensing the electrical output of the secondary electronic conversion apparatus for each of the said plurality of reference fluid flow rates, establishing a correction factor for each of said values of electrical output which will tend to linearize the relationship between fluid flow rate and the electrical output of the measuring system, and storing the correction factors in non-volatile memory in a microprocessor which is part of the secondary electronic conversion apparatus.</td>

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193

METHOD FOR CALIBRATING A DIFFERENTIAL PRESSURE FLUID FLOW MEASURING SYSTEM

The present invention relates to differential pressure based fluid flow measurement, including rate, volume and mass, and more specifically to a method for improving the accuracy of measurement by calibrating the system as a whole, over its operating range, with known fluid flow references.

BACKGROUND OF THE INVENTION

Differential pressure based fluid flow meters, or measurement systems, for fluid filling a confined channel, such as a pipe, typically include two basic elements. A system comprises a fluid immersed sensor, as the primary element of the system, and a secondary element comprising mechanical and electronic means to convert the sensed fluid dynamics into a format usable for obtaining needed information, such as rate of flow, volume of flow and mass of flowing fluid.

Primary sensors which detect differential fluid pressure in the flow include such specific types as averaging pitot tubes, orifice plates, flow nozzles and venturi meters.

The secondary portion of the system may include a differential pressure (DP) transducer, a static pressure transducer, a fluid temperature transducer and a flow computer. The raw electrical output of these transducers is typically processed into a signal which is transmitted for subsequent electronic

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193 processing in additional apparatus, such as a flow computer. The combination of the transducer and the subsequent electronic means, such as A/D converters, amplifiers, and a central processing unit is often referred to in the trade as a DP transmitter and that term will also be used in this specification. Often the primary sensor and the secondary, or DP transmitter, are physically widely separated and are commonly treated as separate instrumentation entities, meaning that the two elements are calibrated separately and present their own individual accuracy standards.

Preliminary to an understanding of the current state of the art of fluid flow measurement, and the accuracy attainable with available systems, vis-à-vis the present invention, it is useful to review the alternative methods of rating the accuracy of measurement devices.

There are two methods of expressing accuracy. The first is usually expressed as “percent of reading” or “percent of value.” This method identifies the error at a specific reading. In measurement devices the use of this method implies that the device produces an error statement based on its readings over the entire operating range. The second method is expressed as the “percent of full scale” (percent of fs) accuracy standard. This method relates the error in a device when it is measuring a quantity that represents 100% of the output. The equivalent “percent of value” error figure for a device calibrated with a “percent of full scale” accuracy statement is calculated using the following equation:

CA 02247407 1998-08-27

WO 97/46851 % Reading Error = °/o FullScale Error VoofScale

X 100

PCT/US97/08193 equation # 1

The differences between these two accuracy ranking methods are not obvious until they are plotted. Figure 1 shows the results of each of the two methods for expressing measurement error. Curve “A” shows a 1% of reading error while curve “B” shows a 1% of full scale error statement. Both are plotted over a 10:1 turndown. In the example shown, a “±1% of full scale” error is equivalent to a “±5% of value” error at a 20% of full range level (5:1 turndown). Flow turndown is the ratio of the highest flow rate expected to be measured by the system to the lowest flow rate expected to be measured, at io some stated accuracy. This quantity is typically expressed on one line with a colon, such as 10:1 for a turndown ratio of 10. Most flow measurement devices have a minimum flow and a maximum flow that can be measured within stated accuracy limits. A typical flow rate falls between these values.

In differential pressure measurement devices, flow rate Q is derived in accordance with the formula Q = C /DP, where C is a constant and DP is the difference between the sensed high and low fluid pressures. Different types of differential pressure sensors have their own benefits and drawbacks, but all of them suffer from an inherent inability to provide fluid pressure outputs which are absolutely related to a fluid flow rate in accordance with the stated .mathematical formula over a wide operating range. The induced errors and ultimate performance characteristics of these differential pressure flow sensors are defined by and vary with the installation parameters of each device, such as

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193 the shape and dimension of the sensor, location and number of its plenum apertures and the sensor’s position inside the fluid carrying pipe. Figure 3 is a diagrammatic curve illustrating in solid line the mathematically ideal relationship between fluid flow Q and differential pressure DP. The dotted line curve illustrates the same relationship but with the inherent errors of a primary sensing element taken into consideration. An averaging pitot tube type of primary sensor or other DP sensor can be made to exhibit a close similarity to the stated mathematical ideal over a wide flow range, but even those types of sensors cannot offer much better than ±1% of reading over a 10:1 turndown.

io In addition to the error introduced by the primary, the components of the secondary of a flow measurement system are each responsible for introducing error into the system. Although error is contributed by all of the secondary components, a large portion of it is induced by the non-linearity of the transducers. Linearity, in this context, is the ability of the components to produce, or closely approximate, a linear relationship between the actual physical input and the DP transmitter output. Other components of the secondary also contribute to non-linearity, resulting in secondary error. The solid line curve of Figure 4 illustrates ideal linearity of the DP transmitter, according to the formula: m.a. = B x DP, where m.a. is the secondary output in milliamperes, B is a constant and DP is the differential fluid pressure. The dotted line curve of Figure 4 diagrammatically illustrates the non-linear performance of the secondary of the flow meter (without reference to the static

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193 pressure or temperature of the fluid) according to the formula: m.a. = B x DP + E_S(DP) where E_s is the % error in differential pressure due to the secondary. While there may be several different representations of linearity, the term is assumed to represent the independent linearity of the secondary of the flow meter. A straight line is used in Figure 5 to minimize the maximum deviation of the actual characteristic. The graph of this Figure has been normalized to show linearity as a percent of reading error deviation from the ideal zero error curve, which would be a straight horizontal line. Current flow measurement practice recognizes the non-linearity of the secondary and attempts to compensate for it by “characterizing” the DP transmitter. This is most often done by employing dead weights, or some other form of calibrated DP source, to simulate various pressures on the transducer diaphragm and then modifying the electrical output to correct for the observed discrepancies. In this type of process however, neither actual fluid flow nor the aberrant characteristics of the primary sensor are given consideration in the linearizing attempt.

Accordingly, when independently calibrated primary and secondary elements are combined, the resulting measurement system accuracy must be determined by combining the error of the primary and secondary elements over the selected range of operation in accordance with the “square root of the sum of the squares” rule:

% System Error - + Es equation #2 where: Ep = % Error in flow due to the primary; and

Eg = % Error in flow due to the secondary, i.e.

S'

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193

Eg = Ej)p2 +Ep2 + Εψ2 where :Εθρ = Error due to the differential pressure transducer,

Ep = Error due to the static pressure transducer, and Εψ = Error due to the temperature transducer.

The primary element’s error contributes directly to the error in flow, whereas the secondary device contributes to error in the differential pressure (DP), which is the square of the error in flow. The percent error in DP will depend upon the percent scale, as shown in Figure 1. For a DP transmitter or meter, the contribution to the secondary’s flow calculation error is:

7oEfs % Eqp — ± 1 + -77.7. -1 X100 equation #3 where: °/aEfs = the % of full scale error (accuracy) of the DP transmitter;

°/oDP = the % of scale at which the DP transmitter is operating.

Ep and Εψ are similarly calculated.

To determine the error in flow at the minimum scale for the desired flow turndown, equation 3 must be substituted into equation 2. To better illustrate the point, reference is made to Figure 2 where the reading flow system accuracy is plotted as a function of turndown for different types of flow meters. An example of a flow meter with a % of reading accuracy, where there is no secondary element contributing to the system error, is shown in solid lines. The error of an orifice plate or other type of DP sensor working with a DP transmitter is calculated using equations 2 and 3, above, and is plotted in the graph of Figure 2 with dotted lines. In this exemplary curve the orifice plate

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193 sensor itself contributes a ±1% of reading error (Ep) and the secondary element has a ±0.1% full scale (fs) error (Es).

After deriving the dotted line curve of Figure 2 , the problem of the prior art custom of separately calibrating the primary and secondary flow meter elements is apparent. The errors of the separate elements are additive within the flow system, creating a combination where the total error can never be less than the least accurate component of the system.

It is therefore the primary object of the invention to provide a synergistic method of obtaining higher measurement system accuracy by calibrating the io assembled combination of primary and secondary system elements with actual reference fluid flow.

A further object of the invention is to provide a method of measurement system calibration that results in unexpectedly better accuracy than that obtained by the square root of the sum of the squares of the individually calibrated system elements.

Another object of the invention is to provide a method for measurement system calibration which results in a wider flow turndown with higher accuracy.

Another object of the invention is provide a process for calibrating a flow measurement system where the flexibility of the secondary element’s linearization potential is advantageously employed to correct for errors in the ?

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193 flow primary, as well as to correct for the nonlinearity in the secondary’s own components.

Other and still further objects, features and advantages of the invention will become apparent upon a reading of the detailed description of a preferred form of the invention.

The most pertinent known prior art is listed as follows;

Dieterich Standard Annubar® Flow Measurement Systems brochure, where, on pages 2 and 3 (unnumbered), a flow measurement system of the type described in this specification is shown, and on page 9, pressure transmitters io and their accuracy are discussed. This measurement system shown does not contemplate the calibration and linearization method of the present invention.

Dieterich Standard Diamond II Annubar® Flow sensors catalog where, on page 3 a diagram of a flow measuring system is shown which is similar to the general system referred to in this specification, but without the linearization method of the present invention, and on page 7 thereof accuracy of the sensor primary is discussed.

U.S. Patent No. 5,111,827 to Bôrje Rantala, disclosing, in a respiratory sampling device, a microprocessor to adjust a flow rate in accordance with the composition of a gas, as determined by a gas analyzer.

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193

U.S. Patent No 4,836,019 to Floyd W. Hagen, disclosing an air data sensor having a computer programmed to take various compensation tables from calibration curves for various sensors. The computer is programmed to provide an appropriate correction factor to the pressure outputs.

U.S. Patent No. 2,869,367 to D.W. Moore, describing a system utilizing pressure responsive diaphragms to change electrical resistances in order to linearize the response characteristic of a system.

THE DRAWINGS

Figure 1 is a graph showing the relationship between “percent of reading” and “percent of full scale” systems of defining measurement accuracy.

Figure 2 is a graph showing the percent of the reading flow system accuracy as a function of turndown for two different types of fluid flow meters.

Figure 3 is a diagrammatic curve illustrating in solid line the mathematically ideal relationship between fluid flow Q and differential pressure DP. The dotted line illustrates the error of primary flow meter measurements.

Figure 4 is a graph illustrative of a typical characteristic curve of the output of a differential pressure transmitter (secondary) as related to differential pressure. The solid line represents a linear relationship between differential pressure and electrical system output. The dashed line is exemplary of an

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193 actual curve, taking into consideration the error Eg introduced by non-linearity in the secondary.

Figure 5 is a normalized curve showing system output linearity as a percent of reading error deviation from an ideal linear curve which would be represented by the straight horizontal line with zero error.

Figure 6 is a block diagram of a typical DP fluid flow rate measurement system adapted to use the calibration method of the present invention.

Figure 7 is a diagrammatic block flow diagram of the initial calibration of the DP transmitter used in a system shown in Figure 6.

io Figure 8 is an exemplary table showing percentage values for differential pressure taken in ten percent increments over the operating range of a DP transmitter where URL is the Upper Range Limit of the transmitter measured in inches of water column. The table shows corresponding values of the initial transmitter characterization, done as shown in Figure 7 (Calib. DP) and the

Corrected DP (DP_corrected)’ shown in Figure 6.

Figure 9 is a block flow diagram of the calibration method of the present invention, with exemplary tables of values taken from indicated points in the process.

/0

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193

Figure 10 is a graph showing the relationship between Dp_measurecj₃ as a percentage of the upper range limit (URL) and the normalized DP, or DP', as a percentage of URL

SUMMARY OF THE INVENTION

The method of the present invention calibrates and characterizes a differential pressure flow measurement system as a unitary assembly, using an actual calibrated flow on the primary sensor. By calibrating the system as a whole, with actual reference fluid flow, instead of relying on the calibration of the primary and secondary elements individually, an unexpected and synergistic io effect is produced whereby the system accuracy and turndown is better than that obtained by the square root of the sum of the squares of the individually calibrated system elements.

DETAILED DESCRIPTION

Before addressing the details of the calibration and characterization method of the present invention, a brief description will be made of a flow measurement system in which the calibration method would be used in order to establish the background of the invention.

Figure 6 is a block diagram of a typical differential pressure dependent flow measurement system of the type to which the calibration method of the present invention is particularly suited. In this measurement system a pitot tube type of differential pressure flow sensor 4 functions as the primary element, //

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193 however traditional differential pressure devices such as an orifice plate, venturi tube or flow nozzle can also act as the differential pressure sensing device. The pitot tube shown is of the type described in U.S. Patent No.

4,154,100 to James B. Harbaugh et al. and U.S. Patent No. 4,559,836 to Darrel

F. Coleman et al., to which patents reference is made for a more complete explanation of the DP flow sensor 4 of Figure 6.

The upstream facing side of the pitot tube 4 senses the average sensor impact pressure of flowing fluid to establish a high pressure value. The downstream facing side of the pitot tube senses low pressure. The high and low io fluid pressures are conducted from the plenums 6 and 8 of the flow sensor 4 to a pressure transducer 10, the first stage of the differential pressure transmitter 11. The pressure transducer transforms the respective high and low fluid pressures issuing out of the primary sensor into an electrical signal whose character is a function of the differential pressure (DP), that is the difference between the sensed high and low fluid pressures. A typical transducer would be one equipped with a sensing element comprising a silicon diaphragm, into the surfaces of which are diffused piezo resistors which comprise an electrical bridge whose output is analogous to the differential pressure applied to the transducer. Other types of DP transmitters may be used, such as capacitance or variable reluctance. The electrical analog signal from the pressure transducer 10 is transformed in the A/D converter 12 to a digital signal for input into a central processing unit (CPU) 16. The CPU 16 performs the square root

CA 02247407 1998-08-27

WO 97/46851 PCT/US97i’08l93 function as well as other calculations, including the system error corrections required for producing an accurate final output signal of the flow measurement system.

After calibration, in accordance with the method of the present invention, and during operation in actual flow measuring conditions, the electronic signal, representing the fluid differential pressure, addresses stored calibration information in an EEPROM memory 18, which information may be in the form of a look-up table or a polynomial equation 19, to provide collective corrections to the errors in that signal which occur as a result of non-linearity in io the transmitter 11 and aberrations in the operation of the primary sensor 4. Obtaining the information to be stored in the EEPROM to implement these corrections is the result of the calibration and characterization process of the present invention and will subsequently be described in detail. The corrected DP signal is further processed by the CPU in the step labeled 27 to determine the square root of the corrected DP. That value is used to determine the flow Q which is then read out in an appropriate instrument 30 or applied to some further process or instrumentation.

Having generally described the function of the flow meter in actual operation, the method of arriving at the calibration coefficients which are stored in the EEPROM will now be described. As mentioned earlier, both the primary sensor and the secondary element contribute error and each adversely effects the accuracy of the system. By providing correction for both of these sources /3

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193 of error at the same time the final flow system signal will have greater accuracy than it would otherwise have by calibrating the primary and secondary elements individually.

With respect to an understanding of the error contributed by the primary sensor, refer again to Figure 3. As for the secondary, the DP transmitter’s nonlinear output is initially linearized against a calibrated DP source, such as a dead weight tester. See Figure 7 for a diagrammatic illustration of an initial linearization process of the DP transmitter by itself. At this point it should be noted that linearization is a form of calibration that is used when the desired io output is a linear function of the input of a device. There are situations, however when the desired output is exponential or a polynomial function or some other function of the input. Therefore, to generically define the calibration method of the present invention, which covers ail such functions, the terms “characterizing” or “characterization” may be used to broadly include is linearization, as well as the others mentioned.

The initial linearization process characterizes the transmitter’s output from the known differential forces supplied by a calibrated DP source. As a result of this initial characterization, the transmitter output will be referred to as

DPcalibrated. As an example of the result of such initial calibration, reference is made to Figure 8 where DP_cai}brated is shown for increments of differential pressure values from zero to 100% of the upper range limit (URL) of the transmitter, as supplied by the calibrated DP source. As part of the initial /7

CA 02247407 1998-08-27

WO 97/46851 PCT/US97/08193 calibration process, these DP_cai values are stored in the non- volatile memory of the transmitter’s CPU. See Figure 6 for a showing of the central processing unit.

However, in accordance with the objects of the present invention, it is 5 not enough to linearize the DP transmitter alone. The objects of the invention are accomplished by characterizing the system DP output against the transmitter’s DP_caijb_rated output, plus the flow coefficient, K, for the primary sensor, that is: DP^al K -^5, DP corrected.

In other words, the initial characterization of the DP transmitter is 10 further enhanced in the context of the total measurement system by taking into consideration in the calibration, the K value variation of the primary sensor.

This is accomplished by storing one or more correction coefficients in the

EEPROM memory 18 and using those coefficients to characterize the measurement system as an integral whole.

Reference to Figure 9 will facilitate an understanding of how the correction coefficients 21 which are stored in the EEPROM are derived. The numbers in the tables which are referenced to different portions of the diagram are exemplary only, for the purpose of better understanding the process. First, the K values for a chosen different number of flows are determined. The number of fluid flow values chosen for the calibration process is a judgment call. For purposes of this description, eleven will be chosen, covering the operating range of the system in ten percent increments, starting with zero.

...... /-U

CA 02247407 ,008-00-27 pQ-f/gg 9//0819)

IPEA/US 17 DEC '97

Although reference fluid flows may be obtained through a number of different techniques, a preferred form involves the use of a highly accurate laboratory weigh tank, which determines mass flow directly. A number of different chosen reference flows are directed into a pipe having a primary sensor, such as the averaging pitot tube 4 of Figure 6. The sensed high and low pressures are directed to the already initially calibrated DP transmitter, from which a differential pressure signal DP_measurecj (DPm) is obtained. This signal, along with an input signal Q, representing a calibrated reference fluid flow from the weigh tank calculates the actual value K of the primary sensor, K_act, according to the formula: K_act = ,---- · From the K_act value for each of the chosen V DPm flow rates an average K (K_aVg) is determined by known statistical analysis methods. In the next step, a nominal DP or DP' is determined according to the ( Kact Ί² formula J x DPm = DP'. Following determination of the nominal DP for each of the chosen flow rates, a statistical regression analysis is used to arrive at one or more correction coefficients for the system which are then stored in the EEPROM. Preferably, the regression uses the least squares method to fit a line through a set of points, each point representing a respective nominal DP value for one of the measured DP values. The dependent variable in the analysis is the measured DP (DPm), while the independent variable is the nominal DP, DP¹. In accordance with well known

AMENDED SHEET

CA 02247407 1998-08-27

WO 97/46851 PCT/US97M8193 techniques the regression analysis can be expanded for increased accuracy by using additional independent variables, DP’2, DP'3, DP'4 and DP'5. The regression analysis yields the X Variables, shown in the table of Figure 9, with their exemplary coefficient values. These are the coefficients which need to be 5 applied to the initial calibration values stored into the transmitter’s EEPROM to determine the new calibration values:

DP_Cor ⁼ Intercept ⁺ Xvarl(DPcai)⁺ Xvar2(DPcal)2 + •••Xvar5(F⁾Fcal)^ where: X_Vam *^st^^e regression coefficient for the nth independent variable;

intercept is the point where the straight line crosses the Y axis.

Depending on the desired accuracy, 2, 3 or 6 independent variables could be used, instead of the five shown in this example.

When the memory is addressed in actual system operation the final corrected DP value (DP_cor) is used to determine the correct flow output.

Calibration is carried out over the full flow range of the measurement system, establishing an accuracy of at least ±0.5% over the full calibrated flow range. The calculated correction values (corrected DP’s) replace the corrections in the EEPROM which were established during the initial characterization of the transmitter, as shown in Figure 7. Characterizing the flow meter with this method of calibration permits the electrical output 27 of the meter to track the true reference flow, thereby compensating for all errors /7

CA 02247407 1998-08-27

WO 97/46851

PCT/ÜS97/08193 due to non-linearity, and the flow coefficient of the primary system, as shown in Figure 10.

CA 02247407 2005-12-01