WO1995020141A1 - Transmitter with improved compensation - Google Patents

Transmitter with improved compensation Download PDF

Info

Publication number
WO1995020141A1
WO1995020141A1 PCT/US1995/000835 US9500835W WO9520141A1 WO 1995020141 A1 WO1995020141 A1 WO 1995020141A1 US 9500835 W US9500835 W US 9500835W WO 9520141 A1 WO9520141 A1 WO 9520141A1
Authority
WO
WIPO (PCT)
Prior art keywords
transmitter
compensation
value
span
membership
Prior art date
Application number
PCT/US1995/000835
Other languages
French (fr)
Inventor
Jogesh Warrior
Scott E. Brigham
Gary A. Lenz
Original Assignee
Rosemount Inc.
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Rosemount Inc. filed Critical Rosemount Inc.
Priority to BR9506548A priority Critical patent/BR9506548A/en
Priority to DE69504036T priority patent/DE69504036T2/en
Priority to EP95908616A priority patent/EP0741858B1/en
Priority to MX9602017A priority patent/MX9602017A/en
Priority to JP51969195A priority patent/JP3557213B2/en
Priority to RU96116924A priority patent/RU2138781C1/en
Publication of WO1995020141A1 publication Critical patent/WO1995020141A1/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D3/00Indicating or recording apparatus with provision for the special purposes referred to in the subgroups
    • G01D3/028Indicating or recording apparatus with provision for the special purposes referred to in the subgroups mitigating undesired influences, e.g. temperature, pressure
    • G01D3/036Indicating or recording apparatus with provision for the special purposes referred to in the subgroups mitigating undesired influences, e.g. temperature, pressure on measuring arrangements themselves
    • G01D3/0365Indicating or recording apparatus with provision for the special purposes referred to in the subgroups mitigating undesired influences, e.g. temperature, pressure on measuring arrangements themselves the undesired influence being measured using a separate sensor, which produces an influence related signal
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D18/00Testing or calibrating apparatus or arrangements provided for in groups G01D1/00 - G01D15/00
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D3/00Indicating or recording apparatus with provision for the special purposes referred to in the subgroups
    • G01D3/02Indicating or recording apparatus with provision for the special purposes referred to in the subgroups with provision for altering or correcting the law of variation
    • G01D3/022Indicating or recording apparatus with provision for the special purposes referred to in the subgroups with provision for altering or correcting the law of variation having an ideal characteristic, map or correction data stored in a digital memory
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S706/00Data processing: artificial intelligence
    • Y10S706/90Fuzzy logic
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10STECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S706/00Data processing: artificial intelligence
    • Y10S706/902Application using ai with detail of the ai system
    • Y10S706/903Control
    • Y10S706/906Process plant

Definitions

  • This invention relates to a technique for compensating a sensed variable, where the variable can be representative of position as in a process automation application, or representative of some other physical variable such as pressure, temperature, pH, optical intensity as in a process control industry application. More particularly, the invention applies to devices, such as transmitters, actuators and positioners, which compensate a sensed variable to provide an output representative of the variable.
  • Measurement transmitters sense process variables such as pressure, temperature, flow, pH, position, displacement, velocity and the like in a process control or process automation installation.
  • Transmitters have analog-to-digital (A/D) converters for digitizing sensor outputs representative of sensed process variable and a compensation circuit for compensating the repeatable errors in the digitized process variable outputs. Temperature is one of the main sources of the error.
  • the compensation circuit typically comprises a microprocessor which calculates the compensated process variable output with long polynomial functions selected to fit the error characteristics of the sensor over a span of pressures.
  • Constants in the long polynomial function are individually selected to each sensor. During manufacture, individual testing of each sensor generates a set of characterization constants related to the sensor errors which is later stored in a transmitter EEPROM. Using this compensation scheme, process variables can typically be corrected to an accuracy of .05% over the span of the primary process variable which the transmitter measures. For example, known pressure transmitters having a span of 0 to 150 inches of water provide corrected pressures within .05% accuracy. Limited electrical power and limited time to compute the output make it difficult to complete more complex computation needed to improve accuracy. Errors in the operating characteristic of the sensor can be a complex, sometimes non-linear function of many variables.
  • the primary variable (the variable which is compensated), contributes directly to the error, while secondary process variables (which affect the measurement of the primary process variables) contribute indirectly to the error.
  • secondary process variables which affect the measurement of the primary process variables
  • contributions of secondary variables become significant.
  • Current approaches solve this quandary with high order polynomials in multiple process variables, but the resulting equation is arithmetically ill-conditioned and sensitive to the manner in which the polynomial is computed, in that overflows may occur.
  • One transmitter compensation equation is an eleventh order polynomial with approximately 100 terms in three variables, which must be calculated each time the transmitter outputs a process variable. Generating characterization constants for these high order polynomials is costly and time consuming. Furthermore, this approach cannot optimally capture the real behavior of the non-linear process variables, which interact nonlinearly.
  • a measurement transmitter has a sensor for sensing a process variable (PV) such as pressure and digitizing means for digitizing an output representative " of the sensed PV.
  • the sensor senses the PV within a span of PV values.
  • a memory inside the transmitter stores at least two membership functions, each membership function having a non-zero value over a predetermined region of the PV span and a substantially zero value over the remainder of the span.
  • the memory also stores a set of compensation formulas, each formula corresponding to a membership function.
  • a selection circuit in the transmitter selects those membership functions which have a non-zero ordinate at the value of the digitized PV and a correction circuit provides at least one correction value, each correction value calculated from a compensation formula corresponding to a selected membership function.
  • a weighting circuit weights each correction value by the ordinate of the corresponding selected membership function, and combines the multiplicands to provide a compensated PV.
  • the compensated PV is coupled to a control circuit connecting the transmitter to a control system.
  • a second embodiment includes a sensor for sensing a primary PV such as differential pressure, and other sensors for sensing secondary PVs such as line pressure and temperature.
  • a set of converters digitize the sensed PVs.
  • Each of the variables is assigned at least one membership function, with at least one of the variables having assigned at least two single dimensional membership functions.
  • the membership functions having a substantially non-zero ordinate at the digitized PV values are selected, and compensation formulas corresponding to the selected membership functions are retrieved from a memory.
  • An FAND circuit forms all unique three element combinations of the ordinates and provides the "rule strength" or minimum of the each of the combinations.
  • a weighting circuit function perform in substantially the same way as described above to provide a compensated primary PV, which is formatted and coupled to a two wire circuit.
  • FIGS. 3A-C are plots of the three membership functions A-C respectively and FIG. 3D is a plot of the all three membership functions A-C, all shown as a function of uncompensated normalized pressure;
  • FIG. 4 is a flowchart of compensation circuit 58 in FIG. 2;
  • FIG. 5 is a block diagram of compensation circuit 58 with an alternative embodiment of membership function selection circuit 64;
  • FIG. 6 is a plot of a multidimensional membership functions
  • FIG. 7 is a plot of the error as a function of pressure for two differential pressure sensors A and B. TABLE 1 shows constants K ⁇ through K 10 for each of the three regions.
  • a pressure transmitter shown generally at 2 transmits an output representative of pressure to a digital control system (DCS) 4 via a two wire current loop shown generally at 6.
  • a fluid 8 in a tank 10 flows through pipe 12 into a series of other pipes 14, 16 and 18, all containing fluid 8.
  • Measurement transmitter 2 senses the pressure difference across an orifice plate 20 situated in the flow of fluid 8. The pressure difference is representative of the flow rate of fluid 8 in pipe 12.
  • a valve 22 located downstream from transmitter 2 controls the flow in pipe 12 as a function of commands received from DCS unit 4 over another two wire loop 24.
  • DCS unit 4 is typically located in a control room away from the process control field installation and in an explosion-proof and intrinsically safe area, whereas transmitter 2 and valve 22 are mounted directly onto- pipe 12 in the field.
  • transmitter 2 is shown with two terminals 50, 52 which are couplable to two terminals of DCS 4 over twisted wire pair 6.
  • DCS 4 is modeled as a resistance and a power supply in series and is shown generally at 4.
  • Transmitter 2 has a sensor section 54 including a capacitance based differential pressure sensor 54A, an absolute pressure sensor 54B and a temperature sensor 54C.
  • Transmitter 2 senses differential pressures between 0 and 250 inches of water.
  • the types of process variables which transmitter 2 measures may include ones representative of position, volumetric flow, mass flow, temperature, level, density, displacement, pH, turbidity, dissolved oxygen and ion concentration.
  • Analog output from sensors 54A-C is coupled to converter circuit 56, which includes voltage or capacitance based analog-to-digital (A/D) converters which can be of the type disclosed in U.S. Patents 4,878,012, 5,083,091, 5,119,033 and 5,155,455, assigned to the same assignee as the present invention.
  • A/D analog-to-digital
  • Each of converters 56A-C generates a serial bitstream of 10 to 16 bits representative of the corresponding digitized process variable (PV) onto a bus connected to compensation circuit 58.
  • Compensation circuit 58 uses fuzzy logic to provide an output representing a compensated PV and typically comprises a microprocessor such as a Motorola 68HC05 with integrated memory. Circuit 58 compensates the errors in the digitized signal representing differential pressure with the digitized signals representing absolute pressure, temperature and differential pressure. Compensation circuit 58 is based on the premise that compensation is most accurately modelled by segmenting the variables to be compensated into multiple regions which overlap each other, where each region has assigned to it a simplified compensation formula optimized for that region and a membership function which can be multidimensional. The "strength" of the formula in the region is variable throughout the region and is described by the ordinate of the membership function at the value of the variable to be compensated.
  • the ordinate of the membership function is typically a number between 0 and 100 percent, indicating the extent to which the value of the variable to be compensated can be modeled by the compensation formula assigned to the selected region.
  • Compensation is determined by first selecting the regions which include the value of the variable to be compensated, and selecting the membership functions and compensation formulas corresponding to each selected region. The next step is to provide a set of correction values, by calculating each of the compensation formulas at the value of the variable to be compensated, and determining the strength of each correction value from the corresponding membership function. Finally, a compensation value is provided by combining the correction values, as weighted by the strength of the membership function at the variable value to be compensated.
  • a membership function selection circuit 64 selects which membership function is non-zero at the digitized P,T,L value and outputs signals representative of the selected membership functions on bus 64B. Circuit 64 also outputs ordinates of the selected membership functions at the digitized P,T,L values (the "rule strengths") at bus 64A.
  • compensation circuit 58 includes at least two single dimensional membership functions for differential pressure, each overlapping the other. If more than one variable is used for compensation, there has to be at least two membership functions for one of the variables.
  • FIGS. 3A-C show differential pressure membership functions A, B and C, each of which have a non-zero value over a different predetermined range of uncompensated pressures within the span.
  • the variable to be compensated (differential pressure) is compensated by the all three variables (P, T and L), but only P is assigned membership functions. (In the most general case, each variable is assigned multiple membership functions.) Membership function A, shown as a solid line in FIG. 3A, has a non-zero value between 0 and 50% span and a zero value thereafter. Membership function B, shown as a dotted line in FIG. 3B, has a non-zero value between 0 and 100% span and a zero value elsewhere. Membership function C, shown as a solid line in FIG. 3C, has a non-zero value between 50% and 100% span and a zero elsewhere.
  • FIG. 3D shows membership functions A, B and C plotted as a function of normalized pressure span.
  • the non-zero segments of membership functions A, B and C define Regions 1, 2 and 3, respectively.
  • the form of the equations need not be the same for each of the regions.
  • the preferred form of the compensation formula for Regions 1-3 to meet the required, accuracy with the metal cell DP sensor is given by Equation 1, which has a second order term as its highest term and requires no more that ten characterization constants.
  • Compensation formula evaluation circuit 66 evaluates and provides a correction value for each of the compensation formulas corresponding to the selected membership functions.
  • the set of characterization constants for each of Regions 1-3 are stored in memory 68 and given below in TABLE 1.
  • Memory 68 is a non-volatile memory containing membership functions, compensation formulas and characterization constants for the compensation formulas.
  • Combining function circuit 70 receives the correction values and the rule strengths arid provides a compensated P process variable according to the equation given by:
  • Output circuit 62 receives and formats the compensated differential pressure PV and couples it to terminals 50, 52 for transmission over process control loop 6.
  • Output circuit 62 may be realized in several ways.
  • a first alternative is a digital-to-analog circuit where the compensated PV is converted to an analog current representative of the compensated PV and is thereafter coupled onto current loop 6.
  • a second alternative is a fully digital transmission, such as Fieldbus, of the compensated PV onto loop 6.
  • a third implementation superimposes a digital signal representative of the PV on an analog current also representative of the PV, such as in the HART® protocol.
  • the number and the functional form of the membership functions are determined by the compensation accuracy required (e.g. .05% accuracy) and the sensor's operating characteristics. For example, a sensor with a significant amount of error which must be compensated requires more membership functions than does a sensor which substantially meets the required amount of accuracy. Membership functions for the sensor which needs more compensation may each have a different functional form (e.g. exponential, gaussian, polynomial, constant, cubic spline, gaussian and logarithmic).
  • a pressure of approximately 30% of span corresponding to an applied pressure of 75.0 inches of water, indicated on FIG. 3D by a solid vertical line and included in the non-zero segments of membership function A and B.
  • Membership functions A and B, corresponding to Regions 1 and 2 are the "selected membership functions".
  • the values of the two membership functions at 30% of span are .359 and .641, respectively.
  • the compensation formulas for Region 1 and 2 are given by Equation 3 and 5: f p (P,T,L) —2.512+278.5154P-4.137T+2.4908 -3.4611P 2
  • the compensated pressure is provided by a combining function, given by
  • Equation 2 above and is 75.112 inches of water, simplified from: p - -359(75.188)+.641(75.070) (7) comp" .359+.641
  • T and L values substituted into the above equation correspond to room temperature and atmospheric line pressure.
  • the resulting correction value from the second order function is insensitive to the manner in which computation takes place (e.g. no overflows), requires less execution time, takes fewer characterization constants and provides more space in memory for additional software functionality in transmitter 2.
  • Another benefit of a fuzzy logic implementation of compensation circuit 58 is to capture the effect of non-linear interaction between variables, which is difficult to model in a prior art single polynomial compensation scheme.
  • the types of variables adapted for use in the disclosed compensation scheme are not limited to sensed PVs.
  • the variable may be a time dependent variable, such as the first or second derivative, or the integral, of the variable.
  • the corresponding membership function would be arranged to provide minimal compensation when the derivative is large (i.e. the magnitude of the compensation is insignificant compared to the magnitude of the pressure change, so it is adequate to approximately compensate the primary PV) .
  • Optimal value stem actuation by a positioner or actuator requires a sensed position and may include a velocity and an acceleration.
  • Another type of variable is a "history dependent" variable, where effects of hysteresis are taken into account. History dependent PVs include information about the previous measurements taken with the specific sensor in transmitter 2. For example, extreme overpressurization of a capacitive based pressure sensor modifies its capacitance as a function of pressure in subsequent measurements.
  • variable is a "position dependent" variable, where the value of the variable changes with position, such as in a diaphragm having one stiffness when bowed and another stiffness in the absence of applied pressure.
  • variable is a "device dependent” variable, where the membership functions and compensation formulas change with the materials used to manufacture transmitter 2. For example, a sensor sensing pressure within a low pressure- ' range has different compensation requirements than does a high range pressure sensor. Similarly, a pressure sensor with a diaphragm made of HASTELLOY® has different error characteristics, and hence requires different compensation, than does one made of MONEL®.
  • the present invention solves inaccuracies in a prior art compensation technique called piecewise linear fitting.
  • piecewise linear fitting the span of the variable of interest is segmented into two or more ranges, and a linear equation is selected for each range which optimally fits each of the ranges.
  • the present compensation scheme with the overlapping membership functions, provides a smooth transition between ranges of the variable of interest.
  • FIG. 4 a flowchart of the functions in compensation circuit 58 is disclosed.
  • the process variables P,T,L are sensed and digitized in blocks 200 and 202 respectively.
  • a counter for counting the number of regions is initialized in block 204.
  • a decision block 206 retrieves the ith membership function from a memory block 208 and determines whether the digitized P,T,L value is in the ith region described by the ith membership function. If the digitized point is included in the region, a computation block 210 retrieves appropriate compensation formulas and characterization constants from memory 208 to compute the ordinate value of a membership function £ mj _(P,T,L) and a correction value f c i(P,T,L) computed from the ith compensation formula, or otherwise increments the region counter i. Decision block 212 causes the loop to re-execute until all the regions which include the digitized P,T,L point are selected. Then block 214 computes the compensated differential pressure as indicated.
  • FIG. 5 details an alternative embodiment of membership function selection circuit 64.
  • fuzzy compensation circuit 58 receives digitized differential pressure (P), digitized absolute line pressure (L) and digitized temperature (T), and uses those three variables to provide a compensated differential pressure.
  • the three main functional blocks are a rule strength circuit 302, a compensation formula evaluation circuit 304 and a combining circuit 306.
  • all of the three variables (P,T,L) are assigned multiple membership functions.
  • differential pressure is assigned four membership functions defined as f p i f p 2 / fp 3 and fp4, # temperature is assigned three membership functions defined as f- i/ ⁇ t2' ano - ⁇ t3' an ⁇ - absolute pressure is assigned two membership functions defined as f ⁇ l and f j ⁇ *
  • Circuit 58 is preferably implemented in a CMOS microprocessor (with adequate on-chip memory), so as to conserve power in the transmitter, which receives power solely from the current loop.
  • Circuit 310 receives the digitized P value and selects those membership functions which have a non-zero ordinate at the digitized P value.
  • circuit 310 is the ordinate of each of the selected membership functions corresponding to the digitized P value, and is labelled at 310A. For example, if the digitized P value were included in the non-zero portion of three of the four P membership functions, then circuit 310 outputs three values, each value being an ordinate of the three selected membership functions corresponding to the digitized P value.
  • bus 310A includes the ordinates: [ f p 2(Po)' f p3 (Pn), f p 4(p ⁇ ) ]•
  • circuit 312 receives the digitized T value and selects temperature membership functions having a non-zero value at the digitized T value. If the digitized T value were included in the non-zero portion of two of the three T membership functions, then circuit 312 outputs two values on bus 312A, each value being an ordinate of a selected membership function.
  • bus 312A includes the ordinate ⁇ : [ f 2( t ⁇ )' ⁇ t S ⁇ o) 3 • Ir ⁇ similar fashion, circuit 314 receives the digitized L value and selects absolute pressure membership functions having a non-zero value at the digitized L value. If the digitized L value were included in both of the two L membership functions, then circuit 314 outputs two values on bus 314A, each value being an ordinate of a selected membership function.
  • bus 314A includes the ordinates: [ fn(-l ⁇ )/ £12( ⁇ -0) 3 • Fuzzy AND circuit 316 forms all unique three element combinations of the ordinates it receives from circuits 310-314 (where each combination includes one value from each of the three busses 310A, 312A and 314A) and outputs the fuzzy AND (the minimum) of each of the unique combinations on a bus 316A.
  • the set of unique membership function ordinate combinations is: f p2 ⁇ P ⁇ ft2 ( t 0 f ndo )
  • the effect of the fuzzy AND circuit 316 is to take single variable membership functions for P, T and L and create multivariable membership functions in P-T-L space. Although it cannot be rendered graphically, circuit 316 creates in P-T-L space a set of 24 three- variable membership functions from the four P, three T and two L single-dimensional membership functions. There are 24 compensation formulas corresponding to the 24 membership functions. In general, the number of multivariable membership functions created is equal to the product of the number of membership functions defined for each individual variable.
  • FIG. 6 gives an example of multivariable membership functions in two variables, P and T. Twelve overlapping pentahedrally shaped two-variable membership functions are defined in P-T space from four triangularly shaped P membership functions and three triangularly T membership functions. Each multivariable membership function corresponds to a compensation formula, and the ordinate of the multivariable membership function (the output of the fuzzy AND) is called a "rule strength" which describes the extent to which the compensated pressure can be modelled with the corresponding
  • Circuit 316 selects those compensation formulas corresponding to each "rule strength" output on bus 316B.
  • Bus 316B has as many signals in it as there are compensation formulas.
  • a "one" value corresponding to a specific compensation formula indicates that it is selected for use in compensation formula evaluation circuit 304.
  • each of the twelve rule strengths defines a point on the surface of twelve. separate pentahedrons, so that twelve compensation formulas (out of a total of 24) are selected.
  • Memory 308 stores the form and the characterization constants for each of the compensation formulas.
  • Compensation formula evaluation circuit 304 retrieves the constants for the selected compensation formulas indicated via bus 316B from memory 308, and calculates a correction value corresponding to each of the selected compensation formulas.
  • Combining circuit 306 receives the correction values and the rule strengths for each of the selected regions and weights the correction values by the appropriate rule strength.
  • the weighted average is given by Equation 4.
  • the characterization constants stored in memory 308 are the result of a weighted least squares fit between the actual operating characteristics of the sensor and the chosen form of the compensation formula for that compensation formula. (The weighted least squares fit is performed during manufacture, rather than operation of the unit. )
  • the weighted least squares fit is given by: b-p- B (8) where b is a nxl vector of calculated characterization coefficients, P is the nxn weighted covariance matrix of the input data matrix X and s is the nxl weighted covariance vector of X with y.
  • the data matrix X is of dimension mxn where each row is one of data vectors representing one of the m (P,T,L) characterization points.
  • FAND circuit 316 is obviated and membership function circuits 310-314 are replaced by three explicitly defined three dimensional membership functions having the form of a radial basis function given generally by:
  • X is a three dimensional vector whose components are the digitized P, T and L values
  • x ⁇ is a three dimensional vector defining the center of the function in P-T-L space
  • controls the width of the function.
  • FIG. 7 shows the sensor error on the respective y axes 400,402 plotted as a function of sensed differential pressure on x axes 404,406 for two pressure sensors A and B (labelled), each connected as shown for pressure sensor 54A in FIG. 2.
  • Sensor A senses a wide range of pressures between 0 and 1000 PSI, while sensor B senses pressure over a tenth of the other sensor's span; from 0 to 100 PSI.
  • the error for sensor A is greater at any given pressure than the error for sensor B at ' the same pressure.
  • a dual sensor transmitter as described here has an output representative of the converted output from sensor B at low pressures, but switches to an output representative of the converted output from sensor A over higher pressures.
  • the present compensation scheme provides a smooth transmitter output when the transmitter switches between the sensors A and B.
  • the output from sensor A is treated as one process variable and output from sensor B is treated as another process variable.
  • each process variable has assigned to it a membership function and a compensation formula, which indicate the extent to which the process variable can be modelled by the compensation formula.
  • a correction value is provided from computing each of the two compensation formulas, and a combining function weights the correction values and provides a compensated pressure.
  • This applicability of the present compensation scheme to dual sensors applies equally well to transmitters having multiple sensors sensing the same process variable, and to transmitters with redundant sensors where each sensor senses a range of PVs substantially the same as the other.
  • the present invention can be applied to devices outside of the process control and process automation industry, and for example could be used to compensate control surface position in an airplane.
  • the type of variables used in the compensation circuit can be other than PVs
  • the compensation formulas and membership functions can be of forms other than polynomials
  • the combining function can be a non-linear averaging function.

Abstract

A circuit (58) for compensating variables in a measurement transmitter (2). Within the transmitter, a sensor (54) senses a primary process variable such as differential pressure and a converter (56) digitizes the sensed process variable. The sensor (54) senses pressures within a span of pressures values. A memory (68) inside the transmitter stores at least two membership functions. The memory (68) also stores a set of compensation formulas, each formula corresponding to a membership function. A selection circuit (64) in the transmitter (2) selects those membership functions which have a non-zero value at the digitized PV, and a correction circuit (66) provides at least one correction value calculated from a compensation formula corresponding to a selected membership function. A weighting circuit (70) multiplies each correction value by its corresponding selected membership function, and combines the multiplicands to provide a compensated PV. The compensated PV is coupled to a control circuit (6) connecting the transmitter (2) to a control system (4).

Description

TRANSMITTER WITH IMPROVED COMPENSATION
BACKGROUND OF THE INVENTION This invention relates to a technique for compensating a sensed variable, where the variable can be representative of position as in a process automation application, or representative of some other physical variable such as pressure, temperature, pH, optical intensity as in a process control industry application. More particularly, the invention applies to devices, such as transmitters, actuators and positioners, which compensate a sensed variable to provide an output representative of the variable.
There is a need to improve the accuracy with which measurement transmitters and devices with actuated outputs, such as a positioner, compensate outputs representative of process variables. Measurement transmitters sense process variables such as pressure, temperature, flow, pH, position, displacement, velocity and the like in a process control or process automation installation. Transmitters have analog-to-digital (A/D) converters for digitizing sensor outputs representative of sensed process variable and a compensation circuit for compensating the repeatable errors in the digitized process variable outputs. Temperature is one of the main sources of the error. The compensation circuit typically comprises a microprocessor which calculates the compensated process variable output with long polynomial functions selected to fit the error characteristics of the sensor over a span of pressures. Constants in the long polynomial function are individually selected to each sensor. During manufacture, individual testing of each sensor generates a set of characterization constants related to the sensor errors which is later stored in a transmitter EEPROM. Using this compensation scheme, process variables can typically be corrected to an accuracy of .05% over the span of the primary process variable which the transmitter measures. For example, known pressure transmitters having a span of 0 to 150 inches of water provide corrected pressures within .05% accuracy. Limited electrical power and limited time to compute the output make it difficult to complete more complex computation needed to improve accuracy. Errors in the operating characteristic of the sensor can be a complex, sometimes non-linear function of many variables. The primary variable (the variable which is compensated), contributes directly to the error, while secondary process variables (which affect the measurement of the primary process variables) contribute indirectly to the error. As the need for accuracy increases, contributions of secondary variables become significant. Current approaches solve this quandary with high order polynomials in multiple process variables, but the resulting equation is arithmetically ill-conditioned and sensitive to the manner in which the polynomial is computed, in that overflows may occur. One transmitter compensation equation is an eleventh order polynomial with approximately 100 terms in three variables, which must be calculated each time the transmitter outputs a process variable. Generating characterization constants for these high order polynomials is costly and time consuming. Furthermore, this approach cannot optimally capture the real behavior of the non-linear process variables, which interact nonlinearly.
In addition to concerns of software and computational complexity, power consumption is critical for transmitters which receive all their operating power over the same wires used for communication. Furthermore, some "intrinsically safe" areas where transmitters are installed limit the transmitter's available power. The finite current budget not only limits the number and complexity of the calculations, but impacts the functionality able to be incorporated in the transmitter. For example, A/D converters could convert digitized sensor outputs more rapidly if more power were available, thereby increasing the transmitter update rate. An EEPROM large enough to accommodate all the characterization constants also consumes power which would otherwise provide additional functionality.
There is thus a need for an accurate method for compensating process variables which is computationally simple and requires small numbers of stored characterization constants, so as to consume a reduced amount of power and provide excess power for additional functionality and increased update rates in the transmitter. SUMMARY OF THE INVENTION
In an embodiment, a measurement transmitter has a sensor for sensing a process variable (PV) such as pressure and digitizing means for digitizing an output representative"of the sensed PV. The sensor senses the PV within a span of PV values. A memory inside the transmitter stores at least two membership functions, each membership function having a non-zero value over a predetermined region of the PV span and a substantially zero value over the remainder of the span. The memory also stores a set of compensation formulas, each formula corresponding to a membership function. A selection circuit in the transmitter selects those membership functions which have a non-zero ordinate at the value of the digitized PV and a correction circuit provides at least one correction value, each correction value calculated from a compensation formula corresponding to a selected membership function. A weighting circuit weights each correction value by the ordinate of the corresponding selected membership function, and combines the multiplicands to provide a compensated PV. The compensated PV is coupled to a control circuit connecting the transmitter to a control system.
A second embodiment includes a sensor for sensing a primary PV such as differential pressure, and other sensors for sensing secondary PVs such as line pressure and temperature. A set of converters digitize the sensed PVs. Each of the variables is assigned at least one membership function, with at least one of the variables having assigned at least two single dimensional membership functions. The membership functions having a substantially non-zero ordinate at the digitized PV values are selected, and compensation formulas corresponding to the selected membership functions are retrieved from a memory. An FAND circuit forms all unique three element combinations of the ordinates and provides the "rule strength" or minimum of the each of the combinations. A weighting circuit function perform in substantially the same way as described above to provide a compensated primary PV, which is formatted and coupled to a two wire circuit. BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a sketch of a field mounted transmitter shown in a process control installation; FIG. 2 is a block diagram of a transmitter made according to the present invention;
FIGS. 3A-C are plots of the three membership functions A-C respectively and FIG. 3D is a plot of the all three membership functions A-C, all shown as a function of uncompensated normalized pressure;
FIG. 4 is a flowchart of compensation circuit 58 in FIG. 2; FIG. 5 is a block diagram of compensation circuit 58 with an alternative embodiment of membership function selection circuit 64;
FIG. 6 is a plot of a multidimensional membership functions; FIG. 7 is a plot of the error as a function of pressure for two differential pressure sensors A and B. TABLE 1 shows constants K^ through K10 for each of the three regions.
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS In FIG. 1, a pressure transmitter shown generally at 2 transmits an output representative of pressure to a digital control system (DCS) 4 via a two wire current loop shown generally at 6. A fluid 8 in a tank 10 flows through pipe 12 into a series of other pipes 14, 16 and 18, all containing fluid 8. Measurement transmitter 2 senses the pressure difference across an orifice plate 20 situated in the flow of fluid 8. The pressure difference is representative of the flow rate of fluid 8 in pipe 12. A valve 22 located downstream from transmitter 2 controls the flow in pipe 12 as a function of commands received from DCS unit 4 over another two wire loop 24. DCS unit 4 is typically located in a control room away from the process control field installation and in an explosion-proof and intrinsically safe area, whereas transmitter 2 and valve 22 are mounted directly onto- pipe 12 in the field.
In FIG. 2, transmitter 2 is shown with two terminals 50, 52 which are couplable to two terminals of DCS 4 over twisted wire pair 6. DCS 4 is modeled as a resistance and a power supply in series and is shown generally at 4. Transmitter 2 has a sensor section 54 including a capacitance based differential pressure sensor 54A, an absolute pressure sensor 54B and a temperature sensor 54C. Transmitter 2 senses differential pressures between 0 and 250 inches of water. However, the types of process variables which transmitter 2 measures may include ones representative of position, volumetric flow, mass flow, temperature, level, density, displacement, pH, turbidity, dissolved oxygen and ion concentration. Analog output from sensors 54A-C is coupled to converter circuit 56, which includes voltage or capacitance based analog-to-digital (A/D) converters which can be of the type disclosed in U.S. Patents 4,878,012, 5,083,091, 5,119,033 and 5,155,455, assigned to the same assignee as the present invention. Each of converters 56A-C generates a serial bitstream of 10 to 16 bits representative of the corresponding digitized process variable (PV) onto a bus connected to compensation circuit 58.
Compensation circuit 58 uses fuzzy logic to provide an output representing a compensated PV and typically comprises a microprocessor such as a Motorola 68HC05 with integrated memory. Circuit 58 compensates the errors in the digitized signal representing differential pressure with the digitized signals representing absolute pressure, temperature and differential pressure. Compensation circuit 58 is based on the premise that compensation is most accurately modelled by segmenting the variables to be compensated into multiple regions which overlap each other, where each region has assigned to it a simplified compensation formula optimized for that region and a membership function which can be multidimensional. The "strength" of the formula in the region is variable throughout the region and is described by the ordinate of the membership function at the value of the variable to be compensated. The ordinate of the membership function is typically a number between 0 and 100 percent, indicating the extent to which the value of the variable to be compensated can be modeled by the compensation formula assigned to the selected region. Compensation is determined by first selecting the regions which include the value of the variable to be compensated, and selecting the membership functions and compensation formulas corresponding to each selected region. The next step is to provide a set of correction values, by calculating each of the compensation formulas at the value of the variable to be compensated, and determining the strength of each correction value from the corresponding membership function. Finally, a compensation value is provided by combining the correction values, as weighted by the strength of the membership function at the variable value to be compensated.
A membership function selection circuit 64 selects which membership function is non-zero at the digitized P,T,L value and outputs signals representative of the selected membership functions on bus 64B. Circuit 64 also outputs ordinates of the selected membership functions at the digitized P,T,L values (the "rule strengths") at bus 64A. As a general rule, compensation circuit 58 includes at least two single dimensional membership functions for differential pressure, each overlapping the other. If more than one variable is used for compensation, there has to be at least two membership functions for one of the variables. FIGS. 3A-C show differential pressure membership functions A, B and C, each of which have a non-zero value over a different predetermined range of uncompensated pressures within the span. The variable to be compensated (differential pressure) is compensated by the all three variables (P, T and L), but only P is assigned membership functions. (In the most general case, each variable is assigned multiple membership functions.) Membership function A, shown as a solid line in FIG. 3A, has a non-zero value between 0 and 50% span and a zero value thereafter. Membership function B, shown as a dotted line in FIG. 3B, has a non-zero value between 0 and 100% span and a zero value elsewhere. Membership function C, shown as a solid line in FIG. 3C, has a non-zero value between 50% and 100% span and a zero elsewhere. FIG. 3D shows membership functions A, B and C plotted as a function of normalized pressure span. The non-zero segments of membership functions A, B and C define Regions 1, 2 and 3, respectively. The form of the equations need not be the same for each of the regions. The preferred form of the compensation formula for Regions 1-3 to meet the required, accuracy with the metal cell DP sensor is given by Equation 1, which has a second order term as its highest term and requires no more that ten characterization constants.
P∞RRECTED ( P. T. L) -K1+K2P+K2 TtKtLtK5P2tK6 T2
(1)
+KJL2+KePL+K9 TP+K10LT
Compensation formula evaluation circuit 66 evaluates and provides a correction value for each of the compensation formulas corresponding to the selected membership functions. The set of characterization constants for each of Regions 1-3 are stored in memory 68 and given below in TABLE 1.
REGION 1 REGION 2 REGION 3
Kl -2.5152 -3.4206 -7.1604
*2 278.5154 283.4241 293.4994 κ3 -4.1357 -2.3884 -0.3094
K4 2.4908 2.5038 2.7488 κ5 -3.4611 -10.5786 -17.4490 κ6 -4.1901 -5.6594 -6.9354 κ7 -0.1319 -0.1589 -0.2082 κ8 11.9573 11.8335 11.4431
Kg -9.3189 -10.3664 -11.5712
κ10 1.1318 1.2281 1.3502
Memory 68 is a non-volatile memory containing membership functions, compensation formulas and characterization constants for the compensation formulas. Combining function circuit 70 receives the correction values and the rule strengths arid provides a compensated P process variable according to the equation given by:
Figure imgf000011_0001
where N is the number of selected regions, W_ is the rule strength for the ith region, f^(P,T,L) is the correction value from the compensation formula corresponding to the ith region and Pcomp represents the compensated differential pressure. Output circuit 62 receives and formats the compensated differential pressure PV and couples it to terminals 50, 52 for transmission over process control loop 6. Output circuit 62 may be realized in several ways. A first alternative is a digital-to-analog circuit where the compensated PV is converted to an analog current representative of the compensated PV and is thereafter coupled onto current loop 6. A second alternative is a fully digital transmission, such as Fieldbus, of the compensated PV onto loop 6. A third implementation superimposes a digital signal representative of the PV on an analog current also representative of the PV, such as in the HART® protocol.
The number and the functional form of the membership functions are determined by the compensation accuracy required (e.g. .05% accuracy) and the sensor's operating characteristics. For example, a sensor with a significant amount of error which must be compensated requires more membership functions than does a sensor which substantially meets the required amount of accuracy. Membership functions for the sensor which needs more compensation may each have a different functional form (e.g. exponential, gaussian, polynomial, constant, cubic spline, gaussian and logarithmic). Consider a pressure of approximately 30% of span, corresponding to an applied pressure of 75.0 inches of water, indicated on FIG. 3D by a solid vertical line and included in the non-zero segments of membership function A and B. Membership functions A and B, corresponding to Regions 1 and 2 are the "selected membership functions". The values of the two membership functions at 30% of span are .359 and .641, respectively. The compensation formulas for Region 1 and 2 are given by Equation 3 and 5: fp (P,T,L) —2.512+278.5154P-4.137T+2.4908 -3.4611P2
(3)
-4.1901-T2-0.1319 **2 + 11.9573P -9.3189-TP+l .1318 T f„ (P, T, L) --3.4206+283.4241P-2.3884T+2.5038 -10.5786P2
(5) -5.694Γ**2-0.1589L**2+11.8335P -10.3664ΓP+1.2281 Γ
Correction values from Equation 3 and 5 are 75.188 and
75.070 inches of water, respectively. The compensated pressure is provided by a combining function, given by
Equation 2 above, and is 75.112 inches of water, simplified from: p - -359(75.188)+.641(75.070) (7) comp" .359+.641
The T and L values substituted into the above equation correspond to room temperature and atmospheric line pressure.
Rather than executing a single eleventh order polynomial as in the prior art, only two second order polynomials are computed. The resulting correction value from the second order function is insensitive to the manner in which computation takes place (e.g. no overflows), requires less execution time, takes fewer characterization constants and provides more space in memory for additional software functionality in transmitter 2. Another benefit of a fuzzy logic implementation of compensation circuit 58 is to capture the effect of non-linear interaction between variables, which is difficult to model in a prior art single polynomial compensation scheme. The types of variables adapted for use in the disclosed compensation scheme are not limited to sensed PVs. The variable may be a time dependent variable, such as the first or second derivative, or the integral, of the variable. In this case, the corresponding membership function would be arranged to provide minimal compensation when the derivative is large (i.e. the magnitude of the compensation is insignificant compared to the magnitude of the pressure change, so it is adequate to approximately compensate the primary PV) . Optimal value stem actuation by a positioner or actuator, such as in a pick and place machine, requires a sensed position and may include a velocity and an acceleration. Another type of variable is a "history dependent" variable, where effects of hysteresis are taken into account. History dependent PVs include information about the previous measurements taken with the specific sensor in transmitter 2. For example, extreme overpressurization of a capacitive based pressure sensor modifies its capacitance as a function of pressure in subsequent measurements. Different compensation formulas apply depending on the severity and frequency of the overpressurization. Another type of variable is a "position dependent" variable, where the value of the variable changes with position, such as in a diaphragm having one stiffness when bowed and another stiffness in the absence of applied pressure. Another type of variable is a "device dependent" variable, where the membership functions and compensation formulas change with the materials used to manufacture transmitter 2. For example, a sensor sensing pressure within a low pressure-' range has different compensation requirements than does a high range pressure sensor. Similarly, a pressure sensor with a diaphragm made of HASTELLOY® has different error characteristics, and hence requires different compensation, than does one made of MONEL®.
The present invention solves inaccuracies in a prior art compensation technique called piecewise linear fitting. In piecewise linear fitting, the span of the variable of interest is segmented into two or more ranges, and a linear equation is selected for each range which optimally fits each of the ranges. Unfortunately, there are typically small. discontinuities, or mismatches, at the boundaries between the separately compensated ranges. The present compensation scheme, with the overlapping membership functions, provides a smooth transition between ranges of the variable of interest. In FIG. 4, a flowchart of the functions in compensation circuit 58 is disclosed. The process variables P,T,L are sensed and digitized in blocks 200 and 202 respectively. A counter for counting the number of regions is initialized in block 204. A decision block 206 retrieves the ith membership function from a memory block 208 and determines whether the digitized P,T,L value is in the ith region described by the ith membership function. If the digitized point is included in the region, a computation block 210 retrieves appropriate compensation formulas and characterization constants from memory 208 to compute the ordinate value of a membership function £mj_(P,T,L) and a correction value fci(P,T,L) computed from the ith compensation formula, or otherwise increments the region counter i. Decision block 212 causes the loop to re-execute until all the regions which include the digitized P,T,L point are selected. Then block 214 computes the compensated differential pressure as indicated.
FIG. 5 details an alternative embodiment of membership function selection circuit 64. Exactly as in FIG. 2, fuzzy compensation circuit 58 receives digitized differential pressure (P), digitized absolute line pressure (L) and digitized temperature (T), and uses those three variables to provide a compensated differential pressure. The three main functional blocks are a rule strength circuit 302, a compensation formula evaluation circuit 304 and a combining circuit 306. However, in this alternative embodiment, all of the three variables (P,T,L) are assigned multiple membership functions. In particular, differential pressure is assigned four membership functions defined as fpi fp2/ fp3 and fp4,# temperature is assigned three membership functions defined as f- i/ ^t2' ano- ^t3' anα- absolute pressure is assigned two membership functions defined as f^l and fj^* Circuit 58 is preferably implemented in a CMOS microprocessor (with adequate on-chip memory), so as to conserve power in the transmitter, which receives power solely from the current loop. Circuit 310 receives the digitized P value and selects those membership functions which have a non-zero ordinate at the digitized P value. Because the non-zero portions of the membership functions may overlap, there is usually more that one selected membership function for each digitized PV. When the membership functions overlap each other by 50%, 2N equations are computed where N is the number of variables which are divided into more than one membership function. The output of circuit 310 is the ordinate of each of the selected membership functions corresponding to the digitized P value, and is labelled at 310A. For example, if the digitized P value were included in the non-zero portion of three of the four P membership functions, then circuit 310 outputs three values, each value being an ordinate of the three selected membership functions corresponding to the digitized P value. Specifically for P = pø, bus 310A includes the ordinates: [ fp2(Po)' fp3(Pn), fp4(pø) ]• At about the same time so as to be effectively simultaneous, circuit 312 receives the digitized T value and selects temperature membership functions having a non-zero value at the digitized T value. If the digitized T value were included in the non-zero portion of two of the three T membership functions, then circuit 312 outputs two values on bus 312A, each value being an ordinate of a selected membership function. Specifically for T = tø, bus 312A includes the ordinateε: [ f 2(t θ)' ^tS^o) 3 • Irϊ similar fashion, circuit 314 receives the digitized L value and selects absolute pressure membership functions having a non-zero value at the digitized L value. If the digitized L value were included in both of the two L membership functions, then circuit 314 outputs two values on bus 314A, each value being an ordinate of a selected membership function. Specifically for L = lø, bus 314A includes the ordinates: [ fn(-lø)/ £12(^-0) 3 • Fuzzy AND circuit 316 forms all unique three element combinations of the ordinates it receives from circuits 310-314 (where each combination includes one value from each of the three busses 310A, 312A and 314A) and outputs the fuzzy AND (the minimum) of each of the unique combinations on a bus 316A. For the set of P, T and L values from the example above, the set of unique membership function ordinate combinations is: fp2<Pθ ft2(t0 fndo)
Figure imgf000017_0001
The effect of the fuzzy AND circuit 316 is to take single variable membership functions for P, T and L and create multivariable membership functions in P-T-L space. Although it cannot be rendered graphically, circuit 316 creates in P-T-L space a set of 24 three- variable membership functions from the four P, three T and two L single-dimensional membership functions. There are 24 compensation formulas corresponding to the 24 membership functions. In general, the number of multivariable membership functions created is equal to the product of the number of membership functions defined for each individual variable. FIG. 6 gives an example of multivariable membership functions in two variables, P and T. Twelve overlapping pentahedrally shaped two-variable membership functions are defined in P-T space from four triangularly shaped P membership functions and three triangularly T membership functions. Each multivariable membership function corresponds to a compensation formula, and the ordinate of the multivariable membership function (the output of the fuzzy AND) is called a "rule strength" which describes the extent to which the compensated pressure can be modelled with the corresponding compensation formula.
Circuit 316 selects those compensation formulas corresponding to each "rule strength" output on bus 316B. Bus 316B has as many signals in it as there are compensation formulas. A "one" value corresponding to a specific compensation formula indicates that it is selected for use in compensation formula evaluation circuit 304. In our specific example, each of the twelve rule strengths defines a point on the surface of twelve. separate pentahedrons, so that twelve compensation formulas (out of a total of 24) are selected. Memory 308 stores the form and the characterization constants for each of the compensation formulas. Compensation formula evaluation circuit 304 retrieves the constants for the selected compensation formulas indicated via bus 316B from memory 308, and calculates a correction value corresponding to each of the selected compensation formulas. Combining circuit 306 receives the correction values and the rule strengths for each of the selected regions and weights the correction values by the appropriate rule strength. The weighted average is given by Equation 4. The characterization constants stored in memory 308 are the result of a weighted least squares fit between the actual operating characteristics of the sensor and the chosen form of the compensation formula for that compensation formula. (The weighted least squares fit is performed during manufacture, rather than operation of the unit. ) The weighted least squares fit is given by: b-p- B (8) where b is a nxl vector of calculated characterization coefficients, P is the nxn weighted covariance matrix of the input data matrix X and s is the nxl weighted covariance vector of X with y. The data matrix X is of dimension mxn where each row is one of data vectors representing one of the m (P,T,L) characterization points.
In an alternate embodiment of compensation circuit 58 shown in FIG. 5, FAND circuit 316 is obviated and membership function circuits 310-314 are replaced by three explicitly defined three dimensional membership functions having the form of a radial basis function given generally by:
Figure imgf000020_0001
In the radial basis function, X is a three dimensional vector whose components are the digitized P, T and L values, x^ is a three dimensional vector defining the center of the function in P-T-L space, and σ controls the width of the function. A set of multidimensional membership functions, such as with radial basis functions, effectively replaces the function of FAND circuit 316, since the FAND circuit provides a set of multidimensional membership functions from sets of single dimensional membership functions.
The present invention is particularly suitable when used in a transmitter with dual differential pressure sensors. FIG. 7 shows the sensor error on the respective y axes 400,402 plotted as a function of sensed differential pressure on x axes 404,406 for two pressure sensors A and B (labelled), each connected as shown for pressure sensor 54A in FIG. 2. Sensor A senses a wide range of pressures between 0 and 1000 PSI, while sensor B senses pressure over a tenth of the other sensor's span; from 0 to 100 PSI. The error for sensor A is greater at any given pressure than the error for sensor B at 'the same pressure. A dual sensor transmitter as described here has an output representative of the converted output from sensor B at low pressures, but switches to an output representative of the converted output from sensor A over higher pressures. The present compensation scheme provides a smooth transmitter output when the transmitter switches between the sensors A and B. In the same fashion as disclosed in FIG3A-D, the output from sensor A is treated as one process variable and output from sensor B is treated as another process variable. As disclosed, each process variable has assigned to it a membership function and a compensation formula, which indicate the extent to which the process variable can be modelled by the compensation formula. A correction value is provided from computing each of the two compensation formulas, and a combining function weights the correction values and provides a compensated pressure. This is a preferred compensation scheme for dual sensor transmitters in that output from both sensors is used throughout a switchover range of pressures, (i.e. no data is discarded for pressures measured within the switchover range) with the relative weighting of the output from each sensor defined by each sensor's membership function. This applicability of the present compensation scheme to dual sensors applies equally well to transmitters having multiple sensors sensing the same process variable, and to transmitters with redundant sensors where each sensor senses a range of PVs substantially the same as the other. Although the present invention has been described with reference to preferred embodiments, workers skilled in the art will recognize that changes may be made in form and detail without departing from the spirit and scope of the invention. The present invention can be applied to devices outside of the process control and process automation industry, and for example could be used to compensate control surface position in an airplane. The type of variables used in the compensation circuit can be other than PVs, the compensation formulas and membership functions can be of forms other than polynomials, and the combining function can be a non-linear averaging function.

Claims

WHAT IS CLAIMED IS:
1. A measurement device related to controlling a process, comprising: sensing means for sensing a PV within a span of PV values and for providing a signal representative of the sensed PV; storing means for storing at least two membership functions, each membership function having a non-zero value over a predetermined region of the PV span and a substantially zero value over a remainder of the span, each membership function having a compensation formula corresponding thereto; compensation means for calculating correction values according to the compensation formula for those membership functions which are non-zero at the value of the sensed PV, the compensation means adjusting each correction value by the corresponding membership function value to provide a weighted correction value and combining the weighted correction values to provide a compensated PV; and an output circuit for coupling the compensated PV to a loop circuit.
2. A measurement transmitter, comprising: conversion means for sensing a PV within a span of PV values and for providing a digitized output representative of the sensed PV; a memory for storing at least two membership functions, each membership function having a non-zero value over a predetermined region of the PV span and a substantially zero value over the remainder of the span, and for storing a set of compensation formulas, each formula corresponding to a membership function; selection means for selecting those membership functions which have a non-zero value at the digitized PV; correction means for providing at least one correction value, each correction value calculated from a compensation formula corresponding to a selected membership function; weighting means for weighting each correction value by its corresponding selected membership function, and for combining the weighted correction values to provide a compensated PV; and an output circuit for coupling the compensated PV to a control circuit.
3. The transmitter of claim 2 where the non-zero regions of two of the membership functions overlap each other.
4. The transmitter of claim 2 where at least one membership function is triangularly shaped.
5. The transmitter of claim 2 where at least one membership function is a gaussian function.
6. The transmitter of claim 2 where the weighting means combine the correction values according to a weighted average.
7. The transmitter of claim 2 where at least one of the compensation formulas is a polynomial function.
8. The transmitter of claim 2 where the PV is one of a set of PVs representative of differential pressure, position, volumetric flow, mass flow, temperature,- level, density, displacement, pH, turbidity, dissolved oxygen and ion concentration.
9. The transmitter of claim 8 where there are three membership functions for the differential pressure, the membership functions each having a center point corresponding to a maximum value of the membership function, where the center points are evenly spaced along the span of PV values.
10. The transmitter of claim 8 where there are three membership functions each overlapping at least one other membership function by 50%.
11. A method for calculating compensated process variables, comprising: sensing a PV representative of a process, the PV taking on values within a predetermined span of PV values; converting the sensed PV to a digitized PV; storing at least two membership functions, each membership function having a non¬ zero value over a predetermined region of the PV span and a substantially zero value over the remainder of the span; storing a set of compensation formulas, each formula corresponding to a membership function; selecting those membership functions which have a non-zero value at the digitized PV; providing at least one correction value, each correction value calculated from a compensation formula corresponding to a selected membership function; weighting each correction value by its corresponding selected membership function, and combining the weighted correction values to provide a compensated PV; and coupling the compensated PV to a control circuit.
12. A method for calculating a compensated variable related to an automated process, comprising: sensing a PV representative of the process, the PV taking on values within a predetermined span of PV values; converting the sensed PV to a digitized PV; storing at least two membership functions, each membership function having a non¬ zero value over a predetermined region of the PV span and a substantially zero value over the remainder of the span; storing a set of compensation formulas, each formula corresponding to a membership function; selecting those membership functions which have a non-zero value at the digitized PV; providing at least one correction value, each correction value calculated from a compensation formula corresponding to a selected membership function; weighting each correction value by its corresponding selected membership function, and combining the weighted correction values to provide a compensated PV; and coupling the compensated PV to a control circuit.
13. A transmitter in a process control system comprising: a sensor for sensing a physical variable in the system, the physical variable taking on values within a span of values, the sensor having an output representative of the physical variable within the span; digitizing means for digitizing the sensor output; selection and storing means for selecting a region from a set of regions, each region related to a different portion of the span, a region becoming selected when the sensor output is included within the region, each region having a compensation formula and a membership function corresponding thereto; computation means for retrieving the membership functions and the compensation formulas corresponding to the selected regions, and for providing the ordinates of the membership functions at the value of the sensor output and for providing the value of the compensation formulas at the value of the sensor output; and compensation means for weighting the value of the compensation formulas by the ordinate of the corresponding membership functions so as to provide a compensated transmitter output.
14. The transmitter of claim 13 where each variable is assigned at least one membership function and the weighting means further comprise: forming means for forming a combination of the ordinates of each of the selected membership functions, each combination comprising one ordinate related to each variable; an FAND circuit for selecting the minimum value in each of the combinations; where the weighting means weight each correction value by the ordinate of the corresponding membership functions to provide a compensated transmitter output.
15. A measurement transmitter, comprising: sensing means comprising at least two sensors, one sensor adapted to sense PVs within a first range of PVs and the other sensor adapted to sense PVs within a second range of PV values, each sensor providing an output signal representative of the sensed PV, the ranges including a maximum and a minimum PV value, the minimum and the maximum defining a transmitter span; storing means for storing at least two membership functions, each membership function having a non-zero value over a predetermined region of the transmitter span and a substantially zero value over a remainder of the transmitter span, each membership function having a compensation formula corresponding thereto; compensation means for calculating correction values according to the compensation formula for those membership functions which are non-zero at the value of the sensed PV, the compensation means adjusting each correction value by the corresponding membership function value to provide a weighted correction value and combining the weighted correction values to provide a compensated PV; and an output circuit for coupling the compensated PV to a loop circuit.
16. The measurement transmitter of claim 15 where the first range is substantially the same as the second range.
17. The measurement transmitter of claim 15 where the first range is larger than the second range.
PCT/US1995/000835 1994-01-25 1995-01-17 Transmitter with improved compensation WO1995020141A1 (en)

Priority Applications (6)

Application Number Priority Date Filing Date Title
BR9506548A BR9506548A (en) 1994-01-25 1995-01-17 Measuring device related to the control of a measurement transmitting process and the process of calculating compensated process variables
DE69504036T DE69504036T2 (en) 1994-01-25 1995-01-17 CONVERTER WITH IMPROVED COMPENSATION
EP95908616A EP0741858B1 (en) 1994-01-25 1995-01-17 Transmitter with improved compensation
MX9602017A MX9602017A (en) 1994-01-25 1995-01-17 Transmitter with improved compensation.
JP51969195A JP3557213B2 (en) 1994-01-25 1995-01-17 Transmitter with improved compensation function
RU96116924A RU2138781C1 (en) 1994-01-25 1995-01-17 Transducer with improved compensation

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US08/186,288 1994-01-25
US08/186,288 US5642301A (en) 1994-01-25 1994-01-25 Transmitter with improved compensation

Publications (1)

Publication Number Publication Date
WO1995020141A1 true WO1995020141A1 (en) 1995-07-27

Family

ID=22684353

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US1995/000835 WO1995020141A1 (en) 1994-01-25 1995-01-17 Transmitter with improved compensation

Country Status (11)

Country Link
US (2) US5642301A (en)
EP (1) EP0741858B1 (en)
JP (1) JP3557213B2 (en)
CN (1) CN1139483A (en)
BR (1) BR9506548A (en)
CA (1) CA2178809A1 (en)
DE (1) DE69504036T2 (en)
MX (1) MX9602017A (en)
RU (1) RU2138781C1 (en)
SG (1) SG44457A1 (en)
WO (1) WO1995020141A1 (en)

Families Citing this family (69)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5642301A (en) * 1994-01-25 1997-06-24 Rosemount Inc. Transmitter with improved compensation
US7254518B2 (en) * 1996-03-28 2007-08-07 Rosemount Inc. Pressure transmitter with diagnostics
JP3790619B2 (en) * 1996-11-29 2006-06-28 松下電器産業株式会社 Processor capable of suitably performing rounding processing including positive value processing and saturation calculation processing
US6104875A (en) * 1997-12-18 2000-08-15 Honeywell Inc. Method for field programming an industrial process transmitter
US6850806B2 (en) * 1999-04-16 2005-02-01 Siemens Energy & Automation, Inc. Method and apparatus for determining calibration options in a motion control system
US6594602B1 (en) * 1999-04-23 2003-07-15 Halliburton Energy Services, Inc. Methods of calibrating pressure and temperature transducers and associated apparatus
US6295875B1 (en) 1999-05-14 2001-10-02 Rosemount Inc. Process pressure measurement devices with improved error compensation
US6508131B2 (en) 1999-05-14 2003-01-21 Rosemount Inc. Process sensor module having a single ungrounded input/output conductor
US6264774B1 (en) 1999-08-23 2001-07-24 Fargo Electronics, Inc. Card laminator with regulatory control
US6865499B2 (en) * 2001-04-26 2005-03-08 Siemens Energy & Automation, Inc. Method and apparatus for tuning compensation parameters in a motion control system associated with a mechanical member
US6859747B2 (en) * 2001-04-26 2005-02-22 Siemens Energy & Automation, Inc. Method and apparatus for self-calibrating a motion control system
US6516672B2 (en) 2001-05-21 2003-02-11 Rosemount Inc. Sigma-delta analog to digital converter for capacitive pressure sensor and process transmitter
ITTO20010705A1 (en) * 2001-07-18 2003-01-18 St Microelectronics Srl SELF-CALIBRATING OVER-SAMPLING ELECTROMECHANICAL MODULATOR AND RELATED SELF-CALIBRATION METHOD.
US7506328B2 (en) * 2002-02-11 2009-03-17 Xerox Corporation Method and system for optimizing performance of an apparatus
US6839546B2 (en) 2002-04-22 2005-01-04 Rosemount Inc. Process transmitter with wireless communication link
US6910381B2 (en) * 2002-05-31 2005-06-28 Mykrolis Corporation System and method of operation of an embedded system for a digital capacitance diaphragm gauge
US7053787B2 (en) * 2002-07-02 2006-05-30 Halliburton Energy Services, Inc. Slickline signal filtering apparatus and methods
US7440735B2 (en) * 2002-10-23 2008-10-21 Rosemount Inc. Virtual wireless transmitter
US6834258B2 (en) * 2002-12-31 2004-12-21 Rosemount, Inc. Field transmitter with diagnostic self-test mode
US6935156B2 (en) * 2003-09-30 2005-08-30 Rosemount Inc. Characterization of process pressure sensor
RU2003133241A (en) * 2003-11-17 2005-04-20 Александр Владимирович Лопатин (RU) DEVICE FOR DETERMINING VALUE OF VARIABLE TIME OF A MEASURED VALUE
WO2005064280A2 (en) * 2003-12-29 2005-07-14 Madison Technologies Limited Measurement results data processing apparatus, system, method, and program
CN1954138B (en) * 2004-03-02 2011-02-16 罗斯蒙德公司 Process device with improved power generation
US8538560B2 (en) * 2004-04-29 2013-09-17 Rosemount Inc. Wireless power and communication unit for process field devices
US7164355B2 (en) * 2004-05-19 2007-01-16 Rosemount Inc. Process transmitter with a plurality of operating modes
US8145180B2 (en) 2004-05-21 2012-03-27 Rosemount Inc. Power generation for process devices
US7262693B2 (en) * 2004-06-28 2007-08-28 Rosemount Inc. Process field device with radio frequency communication
US8160535B2 (en) * 2004-06-28 2012-04-17 Rosemount Inc. RF adapter for field device
WO2006026749A2 (en) * 2004-08-31 2006-03-09 Watlow Electric Manufacturing Company Operations system distributed diagnostic system
US7680460B2 (en) * 2005-01-03 2010-03-16 Rosemount Inc. Wireless process field device diagnostics
US9184364B2 (en) * 2005-03-02 2015-11-10 Rosemount Inc. Pipeline thermoelectric generator assembly
US7334484B2 (en) * 2005-05-27 2008-02-26 Rosemount Inc. Line pressure measurement using differential pressure sensor
US8452255B2 (en) 2005-06-27 2013-05-28 Rosemount Inc. Field device with dynamically adjustable power consumption radio frequency communication
US7379792B2 (en) * 2005-09-29 2008-05-27 Rosemount Inc. Pressure transmitter with acoustic pressure sensor
US7415886B2 (en) * 2005-12-20 2008-08-26 Rosemount Inc. Pressure sensor with deflectable diaphragm
US7308830B2 (en) * 2006-01-26 2007-12-18 Rosemount Inc. Pressure sensor fault detection
US7913566B2 (en) * 2006-05-23 2011-03-29 Rosemount Inc. Industrial process device utilizing magnetic induction
US8188359B2 (en) * 2006-09-28 2012-05-29 Rosemount Inc. Thermoelectric generator assembly for field process devices
US7808379B2 (en) * 2007-03-05 2010-10-05 Rosemount Inc. Mode selectable field transmitter
WO2008144169A2 (en) * 2007-05-14 2008-11-27 Dow Global Technologies, Inc. Faced fiber insulation batt and method of making same
US8898036B2 (en) 2007-08-06 2014-11-25 Rosemount Inc. Process variable transmitter with acceleration sensor
US7484416B1 (en) 2007-10-15 2009-02-03 Rosemount Inc. Process control transmitter with vibration sensor
WO2009154748A2 (en) 2008-06-17 2009-12-23 Rosemount Inc. Rf adapter for field device with low voltage intrinsic safety clamping
US8250924B2 (en) 2008-04-22 2012-08-28 Rosemount Inc. Industrial process device utilizing piezoelectric transducer
RU2453931C1 (en) * 2008-05-27 2012-06-20 Роузмаунт, Инк. Improved temperature compensation of multivariate pressure gage
JP5255698B2 (en) 2008-06-17 2013-08-07 ローズマウント インコーポレイテッド Wireless adapter for field devices with variable voltage drop
US8694060B2 (en) 2008-06-17 2014-04-08 Rosemount Inc. Form factor and electromagnetic interference protection for process device wireless adapters
US8929948B2 (en) 2008-06-17 2015-01-06 Rosemount Inc. Wireless communication adapter for field devices
CN102084626B (en) * 2008-06-17 2013-09-18 罗斯蒙德公司 RF adapter for field device with loop current bypass
US7977924B2 (en) 2008-11-03 2011-07-12 Rosemount Inc. Industrial process power scavenging device and method of deriving process device power from an industrial process
US7870791B2 (en) * 2008-12-03 2011-01-18 Rosemount Inc. Method and apparatus for pressure measurement using quartz crystal
US7954383B2 (en) * 2008-12-03 2011-06-07 Rosemount Inc. Method and apparatus for pressure measurement using fill tube
US8327713B2 (en) 2008-12-03 2012-12-11 Rosemount Inc. Method and apparatus for pressure measurement using magnetic property
US8626087B2 (en) 2009-06-16 2014-01-07 Rosemount Inc. Wire harness for field devices used in a hazardous locations
US9674976B2 (en) * 2009-06-16 2017-06-06 Rosemount Inc. Wireless process communication adapter with improved encapsulation
US8429978B2 (en) 2010-03-30 2013-04-30 Rosemount Inc. Resonant frequency based pressure sensor
US8234927B2 (en) 2010-06-08 2012-08-07 Rosemount Inc. Differential pressure sensor with line pressure measurement
JP2012002741A (en) 2010-06-18 2012-01-05 Yamatake Corp Physical quantity sensor
US8132464B2 (en) 2010-07-12 2012-03-13 Rosemount Inc. Differential pressure transmitter with complimentary dual absolute pressure sensors
CN101936746B (en) * 2010-08-02 2012-09-05 北京华控技术有限责任公司 Non-electric quantity determination system and method
US9057634B2 (en) * 2010-08-11 2015-06-16 Rosemount Inc. Noise detection and avoidance
US10761524B2 (en) 2010-08-12 2020-09-01 Rosemount Inc. Wireless adapter with process diagnostics
US9310794B2 (en) 2011-10-27 2016-04-12 Rosemount Inc. Power supply for industrial process field device
US8752433B2 (en) 2012-06-19 2014-06-17 Rosemount Inc. Differential pressure transmitter with pressure sensor
US9048901B2 (en) 2013-03-15 2015-06-02 Rosemount Inc. Wireless interface within transmitter
EP3066442A1 (en) * 2013-11-04 2016-09-14 Nextnav, LLC Determining calibrated measurements of pressure for different sensors
US10635064B2 (en) * 2014-06-30 2020-04-28 Rosemount Inc. Hysteretic process variable sensor compensation
US9616876B2 (en) 2015-06-29 2017-04-11 Westinghouse Air Brake Technologies Corporation Verification system and method for an end-of-train unit of a train and an improved end-of-train unit for a train
DE102018122014A1 (en) * 2018-09-10 2020-03-12 Endress + Hauser Flowtec Ag Measuring system and measuring arrangement thus formed

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0363958A2 (en) * 1988-10-14 1990-04-18 Hitachi, Ltd. Method and apparatus for controlling the fuel injection for internal combustion engines
JPH02133215A (en) * 1988-11-11 1990-05-22 Mitsubishi Electric Corp Shock absorber control device

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4161880A (en) * 1978-01-05 1979-07-24 Electromedics, Inc. Linearized digital thermometer
US4437164A (en) * 1981-03-05 1984-03-13 Bristol Babcock Inc. Ridge circuit compensation for environmental effects
JPH0774961B2 (en) * 1988-04-07 1995-08-09 株式会社日立製作所 Auto tuning PID controller
US4926364A (en) * 1988-07-25 1990-05-15 Westinghouse Electric Corp. Method and apparatus for determining weighted average of process variable
JPH02132502A (en) * 1988-07-28 1990-05-22 Omron Tateisi Electron Co Working method and adjusting device for fuzzy control device
US5083288A (en) * 1989-02-24 1992-01-21 Arpad Somlyody Apparatus for configuring a remote process sensor and signal transmitter
JP2647217B2 (en) * 1989-12-28 1997-08-27 出光興産株式会社 Composite control method
US5161110A (en) * 1990-02-27 1992-11-03 Atlantic Richfield Company Hierarchical process control system and method
EP0452824A3 (en) * 1990-04-17 1992-12-02 Apt Instruments Corp. Approximate reasoning apparatus
JP2816758B2 (en) * 1990-09-07 1998-10-27 株式会社日立製作所 Apparatus and method for measuring flow rate using fuzzy inference
JPH04119814A (en) * 1990-09-10 1992-04-21 Nissei Plastics Ind Co Temperature control method of injection molding machine
US5251124A (en) * 1991-04-30 1993-10-05 Omron Corporation Fuzzy controller apparatus and method for steady state control
CA2049618A1 (en) * 1991-07-18 1993-01-19 Christopher J. O'brien Integrated transmitter and controller
JPH07105194B2 (en) * 1991-09-24 1995-11-13 岩崎電気株式会社 Fuzzy control device for glass tube sealing process for lamps
JP2654736B2 (en) * 1992-05-20 1997-09-17 株式会社荏原製作所 Dry sludge melting furnace equipment
US5642301A (en) * 1994-01-25 1997-06-24 Rosemount Inc. Transmitter with improved compensation

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0363958A2 (en) * 1988-10-14 1990-04-18 Hitachi, Ltd. Method and apparatus for controlling the fuel injection for internal combustion engines
JPH02133215A (en) * 1988-11-11 1990-05-22 Mitsubishi Electric Corp Shock absorber control device

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
PATENT ABSTRACTS OF JAPAN vol. 14, no. 367 (M - 1008) 9 August 1990 (1990-08-09) *

Also Published As

Publication number Publication date
US5960375A (en) 1999-09-28
CN1139483A (en) 1997-01-01
BR9506548A (en) 1997-08-19
EP0741858A1 (en) 1996-11-13
US5642301A (en) 1997-06-24
RU2138781C1 (en) 1999-09-27
DE69504036T2 (en) 1999-04-15
DE69504036D1 (en) 1998-09-17
MX9602017A (en) 1997-05-31
SG44457A1 (en) 1997-12-19
JP3557213B2 (en) 2004-08-25
JPH09508210A (en) 1997-08-19
CA2178809A1 (en) 1995-07-27
EP0741858B1 (en) 1998-08-12

Similar Documents

Publication Publication Date Title
US5642301A (en) Transmitter with improved compensation
EP0608919B1 (en) Terminal element
US5828567A (en) Diagnostics for resistance based transmitter
EP0757808B1 (en) Smart valve positioner
US6519546B1 (en) Auto correcting temperature transmitter with resistance based sensor
KR100718209B1 (en) Mass flow controller, method for regulating the flow of mass through a mass flow controller and method for calibrating a mass flow controller
CN115950557B (en) Intelligent temperature compensation method based on pressure transmitter
AU608421B2 (en) Process control system
US4569012A (en) Method and system for controlling fluid transportation in a pipe network
EP2483644B1 (en) Process device with sampling skew error correction
CN110346012A (en) A kind of high-lift lock water-level gauge automatic calibrating method
EP1214633B1 (en) Process transmitter with orthogonal-polynomial fitting
EP0384709B1 (en) Learning Machine
EP1256862B1 (en) Iterative method of aircraft sideslip compensation for multi-function probe air data systems
EP0916079A1 (en) Method and apparatus for increasing update rates in measurement instruments
US8401824B2 (en) Method for determining an output value of a sensor of automation technology
EP1259860B1 (en) Method for finding optimal set-points for machines and processes
US7584063B2 (en) Multivariable transmitter and computation processing method of the same
Luceño et al. Computing optimal adjustment schemes for the general tool-wear problem
Cagan FOR ROBUST OPTIMALITY
Svinolupov et al. Pressure control devices with the microprocessor error correction
CN110263465A (en) One kind being used for the probabilistic improvement gray scale algorithm of engineering structure
WO1998057133A1 (en) A method for compensating hysteresis effects in a force measuring device

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 95191314.X

Country of ref document: CN

AK Designated states

Kind code of ref document: A1

Designated state(s): BR CA CN JP MX RU

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH DE DK ES FR GB GR IE IT LU MC NL PT SE

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
WWE Wipo information: entry into national phase

Ref document number: PA/a/1996/002017

Country of ref document: MX

WWE Wipo information: entry into national phase

Ref document number: 2178809

Country of ref document: CA

WWE Wipo information: entry into national phase

Ref document number: 1995908616

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 1995908616

Country of ref document: EP

WWG Wipo information: grant in national office

Ref document number: 1995908616

Country of ref document: EP