1. Technical Field
This invention relates to the estimation of a signal with a physical time delay, and more particularly to a system of reducing or compensating for the time delay caused by a physical impediment in a signal using an infinite impulse response (“IIR”) filter.
2. Background Information
Signals from sensors for a variety of physical phenomena (such as pressure, temperature, flow, acceleration, heat flux, and optical intensity) may be delayed by a physical impediment. For example, in order to measure the temperature of a process fluid flowing through a conduit, a temperature sensor may be positioned in the fluid flow. However, it is often necessary to physically separate the temperature sensor from the fluid flow, e.g., due to compatibility issues. For example, the process fluid to be measured may be chemically incompatible with metallic temperature sensors, e.g., resulting in chemical attack or contamination of the solution and/or electrodes. In addition, the fluid may damage the temperature sensor, and build up of process fluid on the sensor may decrease the sensor's sensitivity. The fluid may also be part of a sanitary process, in which foreign objects such as sensors should not contact the process fluid. These issues may thus tend to preclude the placement of conventional temperature detectors in direct contact with the process fluid.
A conventional approach is to place the temperature sensor within a protective casing. With such a casing, the temperature sensor may be placed within the process fluid flow, while being protected from the process fluid by the casing. This approach relies on thermal conduction through the casing wall to the temperature sensor, to obtain temperature data. A drawback of this traditional solution is that the casing acts as a temperature insulator, thus impeding the sensor's ability to detect temperature change.
In typical examples, casings for temperature detectors to be placed in conduits containing corrosive fluids are fabricated from polymers such as PFA (perfluoroalkoxy polymer resin), PTFE (polytetrafluoroethylene), polyvinyl chloride (PVC), or various combinations thereof, such as perfluoroalkoxy-polytetrafluoroethylene co-polymer. The relatively poor thermal conductivity of these materials tends to adversely affect the accuracy and response time provided by such external temperature detection approaches. Some techniques for compensating for the inaccuracy and delay of such temperature sensors involve the use of additional temperature sensors, including sensors positioned on the outside of the conduit for the process fluid flow. Differences between the signals captured from these multiple sensors may be used to help estimate or otherwise compensate for the time delay. These techniques, however, may be impractical for many applications, such as those involving relatively complicated, expensive casings, such as those which may contain other devices in addition to sensors. It may thus be cost prohibitive to use multiples of these relatively expensive, complicated casings on the conduit.
Referring to the chart of
Turning to
As can be seen in
A need therefore exists for a system that compensates for, or otherwise mitigates the effect of time-related impediments to accurate quality measurements, without the need for multiple sensors.
In an aspect of the invention, a system for predicting, in real time, a physical quality with an impediment to accurate measurement, includes a sensor configured to detect a physical quality (Qdetect), wherein measurement of the physical quality is subject to an impediment. The system includes an infinite impulse response filter (IIR) configured to filter Qdetect and to output a first filtered quality measurement (Qfiltered1) in real time. A processor is configured to calculate the estimated quality Qestimate using Qdetect and Qfiltered1 .
In a variation of the foregoing aspect, the physical quality is temperature of a process fluid, as detected by a resistive temperature detector (RTD), with the physical impediment being a thermally insulative protective casing.
In another aspect of the invention, a method for transforming raw physical quality data into an estimated measured quality, includes obtaining raw data representing a physical quality (Qdetect), in which an impediment exists to the detection of the physical quality. The method also includes filtering Qdetect through an IIR, so that the IIR outputs Qfiltered1; and calculating, the a processor in real time, the estimated measured quality (Qestimate) using Qfiltered1 and Qdetect.
In variations of each of the foregoing aspects, multiple IIR filters may be used to enhance output accuracy.
The above and other features and advantages of this invention will be more readily apparent from a reading of the following detailed description of various aspects of the invention taken in conjunction with the accompanying drawings, in which:
In the following detailed description, reference is made to the accompanying drawings that form a part hereof, and in which is shown by way of illustration, specific embodiments in which the invention may be practiced. These embodiments are described in sufficient detail to enable those skilled in the art to practice the invention, and it is to be understood that other embodiments may be utilized. It is also to be understood that structural, procedural and system changes may be made without departing from the spirit and scope of the present invention. The following detailed description is, therefore, not to be taken in a limiting sense, and the scope of the present invention is defined by the appended claims and their equivalents. For clarity of exposition, like features shown in the accompanying drawings shall be indicated with like reference numerals and similar features as shown in alternate embodiments in the drawings shall be indicated with similar reference numerals.
Briefly, the inventor discovered that an otherwise conventional infinite impulse response (IIR) filter (e.g. “Kalman” or other low pass filters) may be configured to mimic, in real time, a physical system in which a physical impediment (e.g., a time-related impediment) is responsible for unwanted time delays in a sensor's detection of a physical property. As a specific example, the inventor discovered that such an IIR filter may be configured to mimic the effects that the aforementioned RTD (or thermister or thermocouple, etc.) casing would generate in delaying the process temperature from reaching the sensor. In particular, the inventor hypothesized that programming an IIR with a time constant correlated to that of the physical impediment would mimic this physical impediment. The inventor then discovered that this specially configured IIR filter may be used to effectively predict the temperature that would ultimately reach the RTD, to thus substantially eliminate the time delay to provide enhanced, i.e., nearly instantaneous, real time temperature detection.
As alluded to above, using a filter in an attempt to reduce time delay is counterintuitive, since filters themselves tend to slow down and attenuate the actual transmission of a signal passing therethrough (see, e.g.,
These approaches thus capture multiple, sequential outputs from a single temperature sensor, and effectively predict where the temperature of this single temperature sensor will ultimately settle, for improved temperature response times.
Optionally, a second infinite impulse filter may be used in conjunction with the first filter, with the difference between the two filters being used to make further corrections. As will be discussed in greater detail below, the second filter may be a substantial duplicate of the first filter, but which in particular embodiments uses a time constant that may be greater than that used in the first filter. (The time constants used may vary depending on the particular material e.g., casing, through which the temperature must propagate.)
Functionally, the filters are similar to conventional RC circuits in which some part of a new measurement and some part of the previous reading are mixed together linearly (e.g., some fraction of the new measurement and the complementary fraction of the last reading). The response is thus similar to the exponential response of an RC electrical circuit, to effectively provide an exponential, or ‘infinite’ response, from a single sensor. An aspect of the invention is thus the realization that the response coming through the casing to the temperature sensor from a change in the measured temperature may be modeled substantially accurately using this single or double-filtering approach. With the correct time constants, these embodiments enable one to effectively see how much of the remaining (exponential) heat transfer hasn't yet arrived, and then optionally filter again with a longer time constant so the difference is effectively corrected.
As used herein, the term “real time” refers to operations effected nearly simultaneously with actual events, such as to provide results which are delayed nominally only by the execution speed of the processor(s) used.
Turning to
Qestimate=2×Qdetected−Qfiltered1 Eq. 1
Although a single IIR filter may be used as shown, it should be recognized that additional IIR filters may be used, to provide enhanced results in many applications, as will be discussed below.
In particular examples, quality Q may be temperature and constant τ may be a time constant which correlates to a time delay associated with a barrier that is a relatively poor thermal conductor, as discussed below with respect to
Referring now to the embodiment shown as system 102 of
With reference now to
Processor 26 is thus configured to receive both Qfiltered1 20, from IIR1 18, and Qfiltered2 24, from IIR2 22. Processor 26 is also configured to process these inputs to produce estimated temperature, Qestimate 28. In particular embodiments, for example, the processor may be configured to calculate Qestimate 28 by the following Equation 2, i.e., by subtracting the value of Qfiltered2 24 from twice the value of Qfiltered1 20, or by Equation 3, both of which have been found to yield similar results in many applications.
Qestimate=(2×Qfiltered1)−Qfiltered2 Eq. 2
Qestimate=Qfiltered2+(Qfiltered2−Qfiltered1) Eq. 3
Optionally, additional IIRs may be added, such as shown in phantom in
Qestimate=Qfiltered1+(Qfiltered1−Qfiltered2)+(Qfiltered1+Qfilteredadd) Eq. 4
As discussed above, any of the embodiments disclosed herein may be configured in which the physical quantity to be measured is temperature. Such a system is shown at 202 in
In this embodiment, the second IIR filter, IIR2 22, is configured to receive and re-filter Tfiltered1 20, to output Tfiltered2 24. Processor 26 is configured to receive both Tfiltered1 20, from IIR1 18, and Tfiltered2 24, from IIR2 22, and to generate estimated temperature, Testimate 28. Thus, in this example, a raw signal from a single temperature sensor is processed in series through both the first and second IIR filters, creating an “in series” output, and is also processed separately through the first IIR filter, creating a “first filter” output. The “in series” output and the “first filter” output may be combined to determine the estimated temperature.
As mentioned above, in particular embodiments, the processor 26 may be configured to combine the filtered outputs to produce Testimate 28 using the aforementioned Equation 1, 2, 3 and/or 4, in which Q=temperature T.
As discussed above, the time constants are matched to the casing material, e.g., they are set to one-half that of the casing in particular embodiments. The time constant for a particular casing may be discovered empirically, e.g., by sending process fluid of a known temperature through a conduit, and measuring the time required for a temperature sensor in the casing to reach the known fluid temperature.
The various embodiments of the subject invention may provide improved speed and accuracy, relative to conventional approaches by modeling the expected delay to effectively predict where the measured quality (e.g., temperature) will settle. Moreover, any of these embodiments may be further configured to use Equation 5, discussed hereinbelow, to further reduce the time delay for enhanced, nominally real time output.
Embodiments of the claimed invention also involve methods for predicting a physical quality in applications involving a time-related impediment, such as predicting the temperature of a process fluid using a temperature sensor disposed within a casing. For convenience, these methods will be shown and described with respect to temperature measured by a sensor disposed within an insulative casing or with any other barrier disposed between the sensor and the substance being measured. It should be recognized, however, that any one or more of the embodiments shown and described herein may be used to analyze substantially any quality for which speed of detection may be hindered by some physical impediment. Some non-exclusive examples of such qualities may include temperature, pressure, flow, density, concentration, pH, ORP, index of refraction, turbidity, weight, mass, luminosity, position, etc. Moreover, it should be recognized that these qualities may be measured with substantially any type of sensor, including electronic, mechanical, electro-mechanical, chemical, and/or electro-chemical sensors, etc.
Turning now to
Turning now to
Method 500 of
Turning to
In this example, the above-described Testimate 28 more closely tracks the solid process line 30. Although Testimate 28 slightly lags the process 30, Testimate 28 comes close to the amplitude of the actual process temperature. The chart indicates that Testimate 28 reached 100° C. at about 50 seconds, and fell to about 15° C. at about 95 seconds. Since the process fluid 30 was rising in temperature at 95 seconds, Testimate 28 did not fall in temperature to the low of 0° Celsius.
This computer model also allowed the inventor to compare Testimate 28 to the temperature (Tdetected) 16 detected by the temperature sensor, (represented in this
Turning now to
Qestimatefinal=(SlopeQdetected×K)+(Qdetected) Eq. 5
As can be seen from the embodiment of
In this particular example, the value used for the constant K was 40. It should be recognized, however, that constant K is related to the thermal time constant τ discussed above. As such, the value of constant K is expected to change based on the particular application, and may be determined by empirical testing as discussed hereinabove.
It may be seen that in this particular example, slope estimate 34, rather than forming a clear line, appears to be somewhat scattered around the process temperature 30. This may reflect noise which has become part of the slope estimate. Therefore, processor 26 may be further programmed with a conventional smoothing algorithm, e.g., which averages or uses standard deviations of the data to smooth out the estimate 34 for enhanced clarity.
Various embodiments of the present invention have been shown and described herein with reference to use of an electronic sensor 14, such as an RTD. It should be recognized, however, that substantially any type of sensor, including mechanical (e.g., pneumatic), electro-mechanical, and/or electro-chemical sensors / control systems, may be used without departing from the scope of the present invention.
In the preceding specification, the invention has been described with reference to specific exemplary embodiments for the purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to the precise form disclosed. Many modifications and variations are possible in light of this disclosure. It is intended that the scope of the invention be limited not by this detailed description, but rather by the claims appended hereto.
Number | Name | Date | Kind |
---|---|---|---|
3417617 | Rall | Dec 1968 | A |
3487213 | Horan et al. | Dec 1969 | A |
3493949 | Servos et al. | Feb 1970 | A |
3759083 | Erickson et al. | Sep 1973 | A |
3878724 | Allen | Apr 1975 | A |
3892281 | Brown | Jul 1975 | A |
3969943 | Ohno et al. | Jul 1976 | A |
3972237 | Turner | Aug 1976 | A |
4092863 | Turner | Jun 1978 | A |
4140999 | Conway | Feb 1979 | A |
4178798 | Wessel | Dec 1979 | A |
4221125 | Oliver et al. | Sep 1980 | A |
4232682 | Veth | Nov 1980 | A |
4380155 | Paddock et al. | Apr 1983 | A |
4470710 | Crane et al. | Sep 1984 | A |
4537067 | Sharp et al. | Aug 1985 | A |
4556329 | Benson et al. | Dec 1985 | A |
4562554 | Stixrud et al. | Dec 1985 | A |
4605314 | Stenmark | Aug 1986 | A |
4815323 | Ellinger et al. | Mar 1989 | A |
20040023559 | Wu | Feb 2004 | A1 |
20040024559 | Down et al. | Feb 2004 | A1 |
20040192196 | Kim et al. | Sep 2004 | A1 |
20050145273 | Atwood et al. | Jul 2005 | A1 |
20060156830 | Drahm et al. | Jul 2006 | A1 |
20060209921 | Brown et al. | Sep 2006 | A1 |
20070008200 | Hsieh et al. | Jan 2007 | A1 |
20070027639 | Henry et al. | Feb 2007 | A1 |
20070069931 | Hsieh et al. | Mar 2007 | A1 |
20070200193 | Horn | Aug 2007 | A1 |
20070252633 | Frankel et al. | Nov 2007 | A1 |
20080025369 | Patel et al. | Jan 2008 | A1 |
20090277271 | Seppa et al. | Nov 2009 | A1 |
20090277281 | Henry et al. | Nov 2009 | A1 |
Number | Date | Country | |
---|---|---|---|
20110224940 A1 | Sep 2011 | US |