Simplified Interval Observer Scheme: A New Approach for Fault Diagnosis in Instruments

1. Introduction

Processes supervision systems for operators have evolved as new techniques of detection and isolation of faults have appeared. Research in this field has also grown as the complexity of industrial processes has grown and this has motivated the development of different focuses for FDI system design.

The diagnosis of faults can be done using observers. One great advantage of the diagnosis schemes based on observers is that in comparison with other methods they are very large schemes. The high level of complexity in current industrial processes has led to a situation where the amount of information generated by these processes can overcome the capacity of analysis of human operators, which hinders decision making [1,2]. The most recent supervision systems have the capacity to carry out diagnosis and maintenance functions, in order to guarantee the correct functionality of highly complex processes [3–8].

As an example of a process that is hard to supervise, in this work the production of biogas in an anaerobic reactor is used as a case study in which faults are diagnosed and isolated using a scheme based on observers. In many publications about non linear observers for the design of FDI systems, the residuals are based in the error of the estimation obtained by the observer [9]. In biological processes, due to their non linear nature, in the majority of cases they are not completely observable, therefore it is more appropriate to consider some relationships among parameters, instead of attempting to estimate them individually [10,11]. The work presented in [11] explores a methodology to determine the global state and parameters of biological reactors. The method proposed in the article allows one to formalize the design of asymptotic observers, which are capable of evaluating certain variables of state, which are not measured for the anaerobic digestion process, despite certain doubts about the kinetics of the process. We should mention that asymptotic observers need the information of all the input variables of the process, which it is practically impossible to achieve under normal conditions in the operation of an anaerobic reactor. Overcoming this problem implies knowing all the input variables of an anaerobic reactor and for this purpose an observer by intervals was developed in [12]. The main characteristic of intervals observers is that they are capable of providing estimations of guaranteed intervals of non measured state variables instead of an exact estimation, if a superior and inferior limit is provided for each one of the input processes.

2. Fault Diagnosis Schemes

Diagnosis schemes based on observers can be classified according the type of fault detected: sensor faults (Instrument Fault Detection or IFD), actuator faults (Actuator Fault Detection or AFD), and component faults (Component Fault Detection or CFD). Diagnosis schemes can also be classified according the number of observers that are used. There are schemes with one observer: a Direct scheme is a scheme of just one observer of complete order. The Simplified Observer Scheme (SOS), is a scheme of one observer of reduced order. For sensor faults (IFD), the only observer in this scheme uses all the inputs and one output, which only provides simple redundancy and only allows the localization of faults in one sensor. In AFD, the only observer uses all the outputs and just one input. When several observers constitute a bank of observers of reduced order we have a Dedicated Observer Scheme (DOS). For faults in sensors (IFD), each observer uses all the inputs and just one output. The number of observers equals the number of outputs (sensors). For actuator faults (AFD) each observer uses one input and all the outputs. It should be mentioned that the DOS scheme allows the localization of multiple faults, either in sensors (IFD) or in actuators (AFD). The Generalized Observer Scheme (GOS) is formed by a bank of observers of reduced order. For faults in sensors (IFD), each observer uses all the inputs and m-1 outputs, where m is the number of outputs. For actuator faults (AFD), each observer uses all the outputs and n-1 inputs, n being the number of inputs.

3. Design and Implementation of a New Diagnosis Scheme: SIOS-IFD

The SIOS-IFD is a scheme with just one interval observator, of reduced order, for faults in sensors. The main advantage of the SIOS-IFD scheme over all the previously presented schemes, is the fact that no input measurements are is required; it is only necessary to have knowledge of the interval of values that the named inputs can reach. SIOS-IFD only allows the localization of faults in one sensor, because it requires the in line measurement of just one output.

Figure 1 shows a block diagram of the SIOS-IFD. In this figure it can be observed that the SIOS-IFD does not use the inputs u_i (I = 1,2,3,4,…,n), but rather just uses one output (y₃) to estimate the other two outputs (ŷ₁ e ŷ₂) and in that way be able to generate the responses: r₁ = y₁ – ŷ₁ y r₂ = y₂ – ŷ₁.

4. Mathematical Model of the Anaerobic Reactor Used as Case Study

Next the simplified ADM1 mathematic model of the UASB reactor of the Instituto Tecnológico de Orizaba, Veracruz, México (Figure 2), which is the case study of the project of this article, is presented.

Simplified ADM1 Model:

5. Interval Observer Designed for the SIOS-IFD Scheme

In this section the design of the interval observer designed for the IOS – IFD scheme is presented. The designed interval observer is capable of stimating value x₁ y Q_CH4 from the in-line measurement of s₁. It should be emphasized that the designed interval observer does not require any measurement of the reactor input variables. The first necessary condition for the design of an interval observer is that a hypothetic observer of known inputs must exist, called base observer. In order to satisfy this first condition, the asymptotic observer presented next was designed.

The model described by the set of differential non linear Equation (1) can be rewritten in the following way:

The asymptotic observer is designed under the assumption that all inputs are known, and m measures states on-line. Thereby the the space of states can be divided in such way that Equation (2) can be rewritten as:

Equation (3) represents the asymptonic observer that was designed for the anaerobic reactor described for the Equation (1):

In the design of this asymptonic observer, it was assumed that the rate of dilution D, and the inputs to the digestor are considered as knowntherefore A(t) y b(t) are known ∀ t ≥ 0. The asymptonic observer 1 was designed assuming s₁ as the only state measured on-line, meaning, v₂(t) = [s₁] (dim v₂(t) = 1), while x₁ y s₁ were assumed as estimated states, meaning v₁(t) = [x₁ Q_CH4]^T (dim v₁(t) = 2). For the calculation of the matrix N, N₁ was arbitrarily chosen as an identity matrix, to calculate N₂ in the following way: N₂= –N₁C₁C₂^§, (where C₂^§ is the pseudo reverse generalization of C₂), thus obtaining the following result:

Substituting N₁ and N₂ in Equation (3), the following is obtained:

Finally, W_e⁺ was calculated using the minimum value of the dilution rate D(t) = 0.01, obtaining:

W_e⁻ was calculated using the maximun value of the dilution value D(t) = 1, obtaining:

The eigenvalues of W_e⁺ y W_e⁻ are, respectively:

Now, in order to make the observer asymptotically stable, the following conditions must be fulfilled:

• W_e⁻_i,j ≥ 0 ∀i ≠ j.

• W_e⁺ y W_e⁻ must be stable.

Since both conditions are fulfilled, the observer is asymptotically stable. Next the design of an interval observer based on the designed asymptotic observer is presented. Once the asymptotic observer that will work as base observer have been designed, we continued with the design of the interval observer, which is based in the supposition that the values of the input vector b(t) are unknown, but their upper limit b⁺(t) and lower limit b⁻(t) are known, so b⁻(t) ≤ b(t) ≤ b⁺(t). Thus, the asymptotic observer presented in Equation (3) works as a structure base to build the next interval observer:

For the upper limit:

For the lower limit:

The convergence of the interval observer is based in the principal of cooperation defined by [13] (see Lema 1). Let ẽ⁺(t) = v̂₁⁺ – v₁ y ẽ⁻(t) = v₁ – v̂₁⁻ be the errors of estimation associated with Equations (4) and (5), respectively. The follwoing expression represents the dynamics of those errors:

It uses ẽ*(t) to refer to any of the errors ẽ⁺(t) o ẽ⁻(t) so that their dynamics have the same mathematical structure. In the previous equation Ω* = N₁⁻¹ N(b⁺(t) – b(t)) for the case of the upper limit, while Ω* = N₁⁻¹ N(b(t) – b⁻(t)) for the case of lower limit

Lema 1 [11]

In ω̇ = f (ω, t). The system it is cooperative if $\frac{\partial f_h(\omega ,t)}{\partial {\omega }_d}\ge 0 \forall h \ne d$.

This implies that if ω(0) ≥ 0 then ω(t) ≥ 0 ∀t ≥ 0. Therefore, if the initial conditions of the estimated variables are unknown but have their limits organized as: v̂₁⁻(0) ≤ v₁(0) ≤ v̂₁⁺ (0) then ẽ*(0) ≥ 0. Besides, if the system (6) is cooperative, the matrix W_e is stable (Hurwitz) and if Ω* is positive (or zero), it can be guaranteed that ẽ*(t) ≥ 0 ∀t ≥ 0, and in consequence: v̂₁⁻(t) ≤ v₁(t) ≤ v̂₁⁺(t) ∀t ≥ 0 [12].

6. Experimental Results of the Interval Observer

To experimentally verify the operation of the developed interval observer, the anaerobic digester was fed with wastewater from a brewery, which had maximum values of 3 gCOD/L and minimum values of 2 g COD/L, these values are based to the digester being accustomed to consume 3 gCOD/L. If the concentration of the water is greater than this, it can be diluted to achieve the desired value, but if it is lesser, it will be difficult to remove excess water to achieve the desired value, being in the worst case a concentration of 2 gCOD/L. The dilution rate is bounded by a maximum value of 0.74 d⁻¹ and minimum value of 0.26 d⁻¹ (these parameters were obtained experimentally by applying to the plant a positive bounded control as explained by Zavala [14], which sets the value of D by reference to the values of 3 g COD/L and 2 g COD/L for S1, these being the maximum and minimum value that is generated from D). The concentration of organic matter in the effluent (s1) was measured on-line to estimate x1 and Q_CH4. Figure 3 shows the convergence of the range observer for the concentration of anaerobic mass x1; Figure 4 shows the curves, and on-line measurement for the methane biogas flow Q_CH4.

7. Experimental Results of the SIOS-IFD

The developed SIOS-IFD scheme is able to detect unique, sudden and permanent faults in sensors s₁, x₁ and Q_CH4, from response r₁ = y₁ – ŷ₁ and response r₂ = y₂ – ŷ₂, being y₁ = x₁ and y₂ = Q_CH4.

Figure 5 shows the responses to a single, sudden and permanent fault, +5% in the sensor s₁, on day 45 of experimentation. This figure shows that both responses have faults on day 45, because both ◯₁ and Q̂_CH4, have been estimated from s₁ so the two estimates are wrong from the moment in which the fault of s₁ occurs.

Figure 6 shows that each of the responses generated by the diagnostic scheme developed are filtered through a moving window averaging filter. Subsequently, each filteredresponse is evaluated by a threshold detector. The output of each threshold detector is connected to a display, which can only display one of the following values: −1, 0 or 1. The value “−1” indicates that the response reached a value averaged less than the threshold, so, the response had a negative average change. Similarly, the value “1” indicates that the response had an average positive change. Finally, the value “0” indicates that either the responses does not have any media change; or rather variations in their average values were below preset thresholds.

The array of symptoms presented in Figure 6 [11] is the result of the evaluation of the responses generated by the fault diagnosis scheme developed for the example of a single, sudden and permanent fault of +5% in sensors₁. These responses were shown in Figure 5.

Table 2 presents the Structured Diagnostic Matrix for unique faults in the sensors s₁, x₁ y Q_CH4, obtained from the developed response assessment scheme.

Shaded In Table 2 is symptom vector [11], which is the result of the evaluation of the response generated by the fault diagnosis scheme developed for the example of a single, sudden and permanent fault of +5% in the sensor s₁. In this same Table it can be seen that there is a different symptom vector for each single, sudden and permanent fault in sensors s₁, x₁ y Q_CH4.

8. Conclusions

In this article experimental and simulation results of a novel system of diagnosis of faults in sensors are presented. It has been named as SIOS-IFD (Simplified Observer Interval Scheme—Instrument Fault Detection). The main advantage of the SIOS-IFD scheme in comparison with all the other schemes presented above is the fact that the measurement does not require any input, only knowledge of the range of values that such targets can achieve is required. The SIOS-IFD scheme only allows the location of a single sensor fault, for which it requires the in-line measurement in a single output.