In processing plants, gases and liquids are routed through various conduits and tanks. In order to control the operation of the plant, the current status of the liquids and gases in the various portions of the plant as well as the state of controlled devices such as valves and burners must be monitored. The variables monitored to determine these states are referred to generally as process variables and can include pressures, differential pressures, temperatures, tank levels, valve positions, and the like. These process variables are determined using one or more sensors which generate an electrical signal indicative of the process variable. Because processing plants are dynamic, the electrical signals produced for the process variables are also dynamic and tend to vary significantly over time.
In the past, filtering had been used to remove or dampen the variations in the process signals. More recently, statistical process monitoring (SPM) has been introduced in which the dynamics of the process variables are used as one measure of the state of the process variable. Under SPM, a median or mean of the process variable signal can be determined to represent the process variable itself while a standard deviation or coefficient of variation can be calculated to provide a measure of the dynamic nature of the process variable.
Some process variables have a normal amount of variation or a normal standard deviation. When the standard deviation increases above or significantly below the normal amount of standard deviation, it can be indicative of abnormal conditions within a portion of the processing plant. For example, if the standard deviation of a differential pressure sensor drops, it can be indicative of a plugged impulse line.
The statistical values, collectively referred to as SPM data, are generally determined by electronics within field devices such as sensors, transmitters and valve controls that are located within the processing environment. These field devices can be configured to generate alerts when the SPM data crosses a threshold indicating abnormal operation. These alerts are sent to a host system, which can make the alert visible to an operator. Alternatively, the field device can transmit the SPM data to the host system via a digital communication protocol. The host system can be configured to generate alarms when the SPM data crosses various thresholds. In addition, multi-variant statistical analysis can be performed on the SPM data at the host system.
A process field device includes a sensor providing a signal indicative of a process variable and a temperature sensor providing a signal indicative of a temperature of the sensor. Converter circuitry converts the sensor signal into sensor values and the temperature sensor signal into temperature values. A processor implements a digital filter that filters the sensor values to produce filtered sensor values and executes statistical process monitoring instructions that compute a statistical value from the filtered sensor values, wherein the digital filter is temperature-dependent such that the filter changes as the temperature values change.
A sensing system includes a filter construction module that constructs a high pass filter for filtering sensor values indicative of a process variable, the filter construction module setting values for parameters of the filter based on a temperature value indicative of a temperature of the sensor that produced the sensor values.
In further embodiments, a control system includes a temperature-dependent filter that filters a series of sensor values in response to a temperature value to produce filtered values. The control system further includes a statistical process monitoring unit that determines a statistical value from the filtered values and that issues an alert when the statistical value crosses a threshold.
In a method, coefficients for temperature-dependent filter parameters are trained. A sensor signal and a temperature signal are received and a filter is constructed using the temperature-dependent filter parameters and the temperature signal. The sensor signal is applied to the filter to produce a filtered signal and a statistical value is generated from the filtered signal. An alarm is activated when the statistical value crosses a threshold.
The portion 50 of the process plant illustrated in
In any event, one or more user interfaces or computers 72 and 74 (which may be any types of personal computers, workstations, etc.) accessible by plant personnel such as configuration engineers, process control operators, maintenance personnel, plant managers, supervisors, etc. are coupled to the process controllers 60 via a communication line or bus 76 which may be implemented using any desired hardwired or wireless communication structure, and using any desired or suitable communication protocol such as, for example, an Ethernet protocol. In addition, a database 78 may be connected to the communication bus 76 to operate as a data historian that collects and stores configuration information as well as online process variable data, parameter data, status data, and other data associated with the process controllers 60 and field devices 64 and 66 within the process plant 10. Thus, the database 78 may operate as a configuration database to store the current configuration, including process configuration modules, as well as control configuration information for the process control system 54 as downloaded to and stored within the process controllers 60 and the field devices 64 and 66. Likewise, the database 78 may store historical abnormal situation prevention data, including statistical data collected and/or generated by the field devices 64 and 66 within the process plant or statistical data determined from process variables collected by the field devices 64 and 66.
While the process controllers 60, I/O devices 68 and 70, and field devices 64 and 66 are typically located within and distributed throughout the sometimes harsh plant environment, the workstations 72 and 74, and the database 78 are usually located in control rooms, maintenance rooms or other less harsh environments easily accessible by operators, maintenance personnel, etc.
Generally speaking, the process controllers 60 store and execute one or more controller applications that implement control strategies using a number of different, independently executed, control modules or blocks. The control modules may each be made up of what are commonly referred to as function blocks, wherein each function block is a part or a subroutine of an overall control routine and operates in conjunction with other function blocks (via communications called links) to implement process control loops within the process plant. As is well known, function blocks, which may be objects in an object-oriented programming protocol, typically perform one of an input function, such as that associated with a transmitter, a sensor or other process parameter measurement device, a control function, such as that associated with a control routine that performs PID, fuzzy logic, etc. control, or an output function, which controls the operation of some device, such as a valve, to perform some physical function within the process plant. Of course, hybrid and other types of complex function blocks exist, such as model predictive controllers (MPCs), optimizers, etc. It is to be understood that while the Fieldbus protocol and the DeltaV™ system protocol use control modules and function blocks designed and implemented in an object-oriented programming protocol, the control modules may be designed using any desired control programming scheme including, for example, sequential function blocks, ladder logic, etc., and are not limited to being designed using function blocks or any other particular programming technique.
As illustrated in
Additionally, as shown in
Generally speaking, the blocks 80 and 82 or sub-elements of these blocks, collect data, such a process variable data, and perform statistical processing or analysis on the data. For example, the block 80 includes a set of four statistical process monitoring (SPM) blocks or units SPM1-SPM4 which may collect process variable or other data and perform one or more statistical calculations on the collected data to determine, for example, a mean, a median, a standard deviation, a root-mean-square (RMS), a rate of change, a range, a minimum, a maximum, etc. of the collected data and/or to detect events such as drift, bias, noise, spikes, etc., in the collected data. Neither the specific statistical data generated, nor the method in which it is generated is critical. Thus, different types of statistical data can be generated in addition to, or instead of, the specific types described above. Additionally, a variety of techniques, including known techniques, can be used to generate such data. The term statistical process monitoring (SPM) block is used herein to describe functionality that performs statistical process monitoring on at least one process variable or other process parameter, and may be performed by any desired software, firmware or hardware within the device or even outside of a device for which data is collected.
As another example, the block 82 of
The parameters of the SPM blocks within the field devices may be made available to an external client such as workstation 74 through the bus or communication network 76 and the controller 60. Additionally or in the alternative, the parameters and other information gathered by or generated by the SPM blocks within blocks 80 and 82 may be made available to the workstation 74 through, for example, a Communication Server 89. This connection may be a wireless connection, a hardwired connection, an intermittent connection (such as one that uses one or more handheld devices) or any other desired communication connection using any desired or appropriate communication protocol. Of course, any of the communication connections described herein may use a communication server to integrate data received from different types of devices in a common or consistent format.
In current SPM based facilities, alerts are set by setting thresholds for the various SPM parameters such as median/mean and standard deviation. It has been known that field devices set to operate at different sensing ranges should use different thresholds for their alarm limits. In general, sensors that are operating in ranges used to detect smaller changes in values of a process variable should have smaller alarm limits than sensors operating at larger ranges.
The present inventors have further discovered that to take full advantage of SPM the effects of transmitter temperature, range, and manufacturing variability and the effects of impulse line temperature should be minimized. This would provide more repeatable measurements and better SPM alerts under different operating conditions.
Standard deviation values are also affected by the frequency of the process variations within the conduit or tank. At higher frequencies, the sensors are incapable of responding quickly enough to sense the changes in the process variable and as a result, the standard deviation of the process signal at higher frequencies tends to drop. In other words, the standard deviation has non-constant frequency response. In addition, the frequency response of the standard deviation changes with the temperature of the sensor. In general, the sensors respond more slowly when they are colder resulting in a lower standard deviation at colder temperatures across all frequencies.
The frequency response of the standard deviation produced by the sensors is undesirable because it has a significant value at 0 Hz. Because of this, if the frequency response of the standard deviation is used directly, step changes in the process would trigger alerts since the standard deviation would respond strongly to such step changes. To avoid this, the art has applied a high pass filter to the sensor signal before computing the standard deviation so that the frequency response of the standard deviation is 0 at 0 Hz.
Embodiments described herein provide a dynamic high pass filter designed to normalize the standard deviation frequency response to remove variations due to temperature, sensor range, and manufacturing variability. By normalizing the standard deviation frequency response, embodiments described herein make it possible to set common standard deviation alert thresholds for sensors operating at different temperatures, at different ranges, and for sensors manufactured in different ways.
According to an embodiment of the invention,
In this embodiment, sensor 246 senses a process variable such as static pressure, differential pressure, temperature, and valve position, for example. Although only a single sensor is shown for simplicity, device 64 may have multiple different sensors. Analog output from sensor 246 is transmitted to conditioning circuitry 252, which amplifies and conditions (e.g. filters) the signals. Converter circuitry 254 converts the analog signals generated by sensor 246 into digital signals usable by microprocessor 258. As shown in
Microprocessor 264 analyzes the sensor signals to determine a process state. In particular, memory 266 (which may be non-volatile random access memory (NVRM)) includes lookup tables in which are stored algorithm coefficients which are used to determine particular values of process state, such as mass flow rate, based on the magnitude of the sensed process variables. Additionally, hardware parameters and process fluid parameters, such as the type and the bore diameter of a pipe carrying the process fluid and the viscosity and density of the process fluid, are uploaded into memory 266 through control loop 218. In other embodiments, data relating to hardware parameters and fluid parameters are directly entered into transmitter 264 through a user interface (not indicated in
Using hardware parameters, process fluid parameters, sensed process variables, and algorithm coefficients, microprocessor 264 performs a run-time calculation. The runtime calculation evaluates a process condition equation, to determine a process condition or state of the fluid, such as the mass flow rate. A process condition signal representing the calculated process condition and the sensor signals are transmitted to workstation 74 over control loop 218 using communication circuitry 268. Communication circuitry 268 includes voltage regulator 268A, modulator circuitry 268B, loop current controller 268C and a protocol receiver, such as a 4-20 mA HART® receiver or transceiver 268D, to enable transmitter circuitry module 234 to communicate with communication interface 280 of workstation 74.
NVRAM 266 also contains instructions for statistical process monitoring (SPM). These instructions are executed in parallel with the instructions used to generate the process state. The SPM instructions compute statistical values from the sensor signals such as the mean/median, standard direction, etc. of one or more of the sensor signals. In particular, and as detailed below, the SPM instructions cause microprocessor 264 to perform steps in which a high pass filter is dynamically constructed to provide a desired frequency response to the statistical values created by the SPM. In accordance with some embodiments, the constructed high pass filter compensates for temperature changes at the sensor, changes in the sensor's frequency response due to the range setting of the sensor, and manufacturing variabilities in the sensor.
In accordance with some embodiments, NVRAM 266 also holds alert thresholds and microprocessor 264 is able to determine when a statistical value crosses a threshold, when this occurs, microprocessor conveys the alert along with a Time Stamp to workstation 74 using communication circuitry 268.
In the alternative, or in addition, the statistical values produced by microprocessor 264 are sent to workstation 74 through communication circuitry 268. Workstation 74 can store the statistical values in a memory 276 and can display historical statistical values on an output monitor 274. Microprocessor 270 of workstation 74 can compare the statistical values to thresholds stored in memory 276 and can issue an alert when a statistical value crosses a threshold. Note that a threshold may be crossed by going above or below a given threshold.
Workstation 74 typically comprises a personal computer, such as a portable computer, installed near an industrial process to monitor and regulate process variables and process conditions of the industrial process. Workstation 74 includes inputs 272, such a keyboard, mouse, or other user interface, that enable operators to input process control inputs into memory 276. Workstation 74 also includes outputs 274, such as a monitor, which enable data to be extracted from workstation 274. Peripherals 278, such as a printer or other commonly used devices, can also be connected to workstation 74 to input and extract data from microprocessor 270. Workstation 74 includes communications interface 280, which transmits data to and receives data from microprocessor 264. Workstation 74 is configured to provide complete evaluation and analysis of the capabilities of module 234 and the sensor signals generated by sensor 246.
Frequency response measurement unit 910 then instructs transducer 914 to change the frequency of the fluid variation and the sensor's resulting measurements are converted into a standard deviation once again. This new standard deviation value is then stored by frequency response measurement unit 910. In this way, frequency response measurement unit 910 is able to determine a standard deviation value for a plurality of frequencies and thus can determine the frequency response of the standard deviation for the sensor at a given temperature. When the complete frequency response has been determined at a selected temperature, frequency response measurement unit 910 instructs sensor temperature control 912 to change the temperature of the sensor to a new temperature value. At the new temperature value, frequency response measurement unit 910 once again instructs transducer 914 to vary the fluid at various different frequencies to obtain a frequency response of the standard deviation at the new temperature. By repeating these steps for each temperature, frequency response measurement unit 910 produces sensor frequency responses 916, which includes a frequency response for the sensor at each selected temperature.
At step 1002, a target standard deviation frequency response 918 is retrieved or computed by microprocessor 264. In accordance with one embodiment, the target frequency response is the standard deviation frequency response of a filtered sensor signal produced by a model sensor operating at a selected temperature and at the same range as the sensor in the field device and filtered by a high pass difference filter. An example of such a target frequency response is shown in
In other embodiments, the target frequency response is the standard deviation frequency response of a filtered sensor signal produced by a model sensor operating at a selected temperature and in a selected range regardless of the range of the present sensor and filtered by a high pass difference filter. Because the model sensor is operating at a selected range, the range of the model sensor will be different from the range of the sensor in the field device at times. By selecting a target frequency response produced by a model sensor operating at a selected temperature and at a selected range without regards to the range of the sensor in the field device, the frequency response of the standard deviation can be normalized so that variations due to temperature and operating range of the sensor in the device are removed.
In still further embodiments, the target frequency response can be the standard deviation frequency response formed by a combination of a model sensor with a constant frequency response and a model high pass filter such as a high pass difference filter.
At step 1004, a filter fitting unit 922 constructs filter parameters 924 for a high pass filter constructed from the target frequency response and the sensor frequency response at each temperature. Thus, filter fitting unit 922 constructs a separate high pass filter for each temperature.
Filter fitting 922 constructs the filter by first generating a frequency response for the filter as:
Where D(ω) is the target standard deviation frequency response, A(ω,T) is the sensor frequency response at temperature T and H(ω,T) is the frequency response for the constructed filter at temperature T.
Once the frequency response for the filter has been determined using EQ. 1 above, parameters for the high pass filter are determine to achieve the frequency response for the filter. In accordance with one embodiment, the high pass filter is constructed as a finite impulse response filter which takes the form of:
Where N is the order of the filter and b0, b1, b2, . . . bN are the coefficients of the filter.
Within the class of FIR filters there is a set of filters known as Generalized Linear Phase (GLP) filters. The linear phase property is often desirable in digital filters because it reduces the amount of distortion in, for example, audio signals. There are four types of GLP filters. The GLP Type IV filter is generally well suited for diagnostics because it has exactly 0 response at DC (0 Hz), but can be given any arbitrary response at the Nyquist frequency which is one half the sampling frequency. The GLP Type IV filter requires an odd order N which results in an even number of filter coefficients. The GLP Type IV filter also has a property of anti-symmetry in the coefficients meaning that bi=−bN-i. The high pass differencing filter
is actually a very simple GLP Type IV filter with order N=1. Using a heuristic approach, the present inventors have discovered that an order 5 GLP Type IV filter provides good curve fitting to the standard deviation frequency response computed using EQ. 1 at each of the different temperatures.
An order 5 GLP Type IV filter is expressed as:
yk=boxk+b1xk-1+b2xk-2−b2xk-3−b1xk-4−b0xk-5 EQ. 3
Where yk is the output of the filter at xk, xk-1, xk-2, xk-3, xk-4, and xk-5 are the sensor measurement values at xk, xk-1, xk-2, xk-3, xk-4, and xk-5, respectively, and b0, b1, and b2 are the parameters of the high pass filter. Note that although there are six multiply operations in EQ. 3, there are only three unique coefficients b0, b1, and b2. As a result, filter fitting unit 922 only has to determine the values of b0, b1, and b2 when constructing a filter to achieve the standard deviation frequency response computed by EQ. 1.
Any known numerical digital filter design algorithm, such as Window Design, Frequency Sampling, Weighted Least Squares, Equiripple, Remez Exchange, and Parks-McClellan, for example, may be used to identify the values of that best fit the filter computed in Eq. 1.
The values of b0, b1, and b2 determined at step 1004 are stored as filter parameters 924 at step 1006 for each temperature.
The vertical axis in
At step 1008, a parameter curve fitting 926 fits a curve to the parameter values shown in
bi=ai,0+ai,1t+ai,2t2+ai,3t3,for i=0,1,2 EQ. 5
Where t represents the normalized temperature, ai,0, ai,1, ai,2, and ai,3 are the coefficients for the filter parameter curves and bi is the filter parameter.
At step 1010, a coefficient matrix is constructed from the temperature curves to produce coefficient matrix 928. An example of such a coefficient matrix is:
Note that the coefficient matrix is not necessary but aids in the understanding and the computation of the filter parameters.
After the coefficient matrix has been constructed, the field device is prepared to dynamically construct temperature-dependent high pass filters to reach the target frequency response selected at step 1002 for any temperature of the sensor. At step 1012, temperature value 932 representing the temperature of the sensor measured by PRT 256 is provided to filter construction unit 930. At step 1014, filter construction unit 930 uses EQ. 4 to determine a normalized temperature. At step 1016, filter construction unit 930 uses the coefficient matrix 928 and the normalized temperature to calculate the values of the filter parameters, b0, b1, and b2. In particular, the coefficient matrix of EQ. 6 multiplies a temperature matrix of:
Where T is the temperature matrix and t is the normalized temperature.
At step 1018, the filter parameters are stored as high pass filter 902 and are applied to measurements 900, which may be stored in NVRAM 266 or may be processed as they are received by processor 264 from processor 258. The measurements 900 comprise a series of sensor values that together represent a sensor signal. The output of high pass filter 902 is filtered signal 904, which is applied to standard deviation calculation 906 at step 1020 to form a standard deviation 908 that has the target frequency response. At step 1022, the standard deviation is compared to the alarm thresholds and if the standard deviation crosses the alarm threshold, an alarm is activated.
Although reference to standard deviation is made above, other statistical metrics may be used such as maximum and minimum, for example, in place of the standard deviation.
Although the temperature of the sensor is used above, in other embodiments, the temperature of different parts of the system may be used in place of the temperature of the sensor. For example, in some embodiments, the temperature of an impulse line connected to a sensor may be used to construct the high pass filter. In accordance with one particular embodiment, a pressure sensor is positioned near a valve to sense cavitation at the valve. In such embodiments, the pressure sensor is connected to the valve through an impulse line. By measuring the temperature of the impulse line, the high pass filter can be constructed to compensate for temperature changes in the impulse line that affect the SPM data generated from the pressure sensor signal. This is an example of compensating on a system level versus on a transmitter or field device level.
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.
Number | Name | Date | Kind |
---|---|---|---|
4607530 | Chow | Aug 1986 | A |
5735604 | Ewals | Apr 1998 | A |
6446444 | Chase | Sep 2002 | B1 |
7180211 | Sinha | Feb 2007 | B2 |
7949495 | Wiklund et al. | May 2011 | B2 |
20050149264 | Tarvin | Jul 2005 | A1 |
20060047480 | Lenz et al. | Mar 2006 | A1 |
20070010900 | Kavaklioglu | Jan 2007 | A1 |
20070186657 | Sato | Aug 2007 | A1 |
20080082294 | Pihlaja | Apr 2008 | A1 |
20080208527 | Kavaklioglu | Aug 2008 | A1 |
20090055457 | Miller | Feb 2009 | A1 |
20090299543 | Cox | Dec 2009 | A1 |
20120303301 | Park | Nov 2012 | A1 |
Number | Date | Country |
---|---|---|
101078913 | Nov 2007 | CN |
1 285 841 AL | Feb 2003 | EP |
2009-505276 | Feb 2009 | JP |
WO 2006107933 | Oct 2006 | WO |
WO 2010039225 | Apr 2010 | WO |
Entry |
---|
First Chinese Office Action for Chinese Patent Application No. 201420191899.5, dated Jul. 4, 2014, 2 pages. |
Notification of Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, PCT/US2014/049017, dated Nov. 10, 2014. |
Communication from European Patent Application No. EP 14755223.6, dated May 2, 2016. |
Office Action from Canadian Patent Application No. 2,923,136, dated May 30, 2016. |
Office Action from Chinese Patent Application No. EP 201410157723.2 dated Oct. 8, 2016. |
Number | Date | Country | |
---|---|---|---|
20150088279 A1 | Mar 2015 | US |