Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
In Tube Integrated Electronic Nose System on a Flexible Polymer Substrate
Next Article in Special Issue
Sensing Performance of Precisely Ordered TiO2 Nanowire Gas Sensors Fabricated by Electron-Beam Lithography
Previous Article in Journal
Optimization of Passive Low Power Wireless Electromagnetic Energy Harvesters
Previous Article in Special Issue
The Effect of the Thickness of the Sensitive Layer on the Performance of the Accumulating NOx Sensor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Overcoming the Slow Recovery of MOX Gas Sensors through a System Modeling Approach

by
Javier G. Monroy
1,*,
Javier González-Jiménez
1 and
Jose Luis Blanco
2
1
Department of System Engineering and Automation, University of Málaga, Campus de Teatinos, 29071 Málaga, Spain
2
Department of Engineering, University of Almería, 04120 La Cañada, Almería, Spain
*
Author to whom correspondence should be addressed.
Sensors 2012, 12(10), 13664-13680; https://doi.org/10.3390/s121013664
Submission received: 2 August 2012 / Revised: 2 October 2012 / Accepted: 5 October 2012 / Published: 11 October 2012
(This article belongs to the Special Issue Gas Sensors - 2013)

Abstract

: Metal Oxide Semiconductor (MOX) gas transducers are one of the preferable technologies to build electronic noses because of their high sensitivity and low price. In this paper we present an approach to overcome to a certain extent one of their major disadvantages: their slow recovery time (tens of seconds), which limits their suitability to applications where the sensor is exposed to rapid changes of the gas concentration. Our proposal consists of exploiting a double first-order model of the MOX-based sensor from which a steady-state output is anticipated in real time given measurements of the transient state signal. This approach assumes that the nature of the volatile is known and requires a precalibration of the system time constants for each substance, an issue that is also described in the paper. The applicability of the proposed approach is validated with several experiments in real, uncontrolled scenarios with a mobile robot bearing an e-nose.

1. Introduction

The deployment of olfactory sensors is becoming an increasing practice in many industrial and environmental applications due to advances in the gas sensing technology. The exploitation of olfactory sensors can be classified into two main groups according to the level of control over the measurement conditions: Closed Sampling Systems (CSS), where the gas sensors are usually hosted in test chambers with controlled airflow, volatile exposure times, temperature and humidity, etc., and Open Sampling Systems (OSS), with no control over the sensing conditions. Our interest is in the latter, which are more flexible and practical for field applications. Examples of such uses are environmental exploration [1], gas distribution modeling [2], buried land mine detection [3] or pollution monitoring [4]. Some of these applications are usually accomplished with the help of a mobile robot carrying the sensors on board, which makes the sensing task even more challenging.

Within the different technologies and materials available for gas sensor fabrication [5], MOX (Metal Oxide Semiconductor) transducers are one of the most popular and widely employed in mobile robotics olfaction, due to their high sensitivity and low prices. However they present some shortcomings including poor selectivity, response drift (age factor), influence by environmental factors such as humidity and temperature [6] and major limitations in their response speed [7]. These limitations come from the sensing mechanism underlying MOX technology, that is, the exchange of oxygen molecules between the volatile and the metal film [8,9].

Among these drawbacks, the long duration of the acquisition cycles (up to tens of seconds) is of special concern for OSS, since inaccurate readings are inevitable when measuring rapid changes of gas concentration, as illustrated in Figure 1. Observe how this limitation is particularly noticeable in the decay phase, when the output recovers to the baseline level (the steady output value given by a gas sensor when exposed to clean air). As a consequence of this slow dynamic response and because of the intermittent and chaotic nature of turbulent airflow in OSS [10,11], steady state values are rarely reached, and therefore gas sensing based on MOX technology must deal with the transient information of the signals [12]. This problem becomes crucial when the sensors are carried on a vehicle (typically a mobile robot) to provide measurements along the way. The adopted solution in such cases is to reduce the vehicle velocity to a few cm/s, such in [13]. This proposal, however, is not acceptable in many applications since the sampling of space must be as quickly as possible to cope with the rapid dynamics intrinsic to gas propagation.

To overcome this shortcoming of MOX-based electronic noses, we propose to estimate the steady state sensor output from the noisy and distorted transient signal, which corresponds to the search for an inverse dynamical model. In this work we follow the method to obtain, first, a forward model of the system and then invert this model. For that, we rely on a linear sensor model that comprises two first-order systems corresponding to the rise and recovery phases, respectively, with variable time constants. The applicability of the proposed approach to Open Sampling Systems is validated through three different experiments in real scenarios, performed with a mobile robot bearing an MOX-based e-nose.

2. Related Works

In general, a model seeks to represent a system (empirical objects, phenomena, and physical processes) in a logical, objective and simplified way, allowing to predict the output of the system provided the input. Thus, a model of an MOX sensor must predict the sensor resistance (transient and steady state) when exposed to a certain gas concentration profile. Our interest in having such a model is to use it in a reverse way: given a sequence of measurements from the transient response of the MOX sensor, we seek the exciting gas distribution through the estimation of the steady state sensor resistance.

The modeling of the dynamics of MOX sensors has been addressed in the literature for a variety of purposes. Gardner et al. [14] proposed a non-linear diffusion-reaction model to obtain the theoretical transient and steady state responses based on the reactions taking place at semiconductor level. This model is not applicable to our approach since we pretend to obtain the gas distribution from the sensor readings rather than from the electrical and physical properties of the sensors. Later, aiming at increasing the response speed of gas sensors, T. Yamanaka et al. [15] reported a two-phases (corresponding to the rise and decay phases) second order linear model to describe the transient response of a semiconductor gas sensor from a visualized gas distribution image. Despite its success, this model requires the use of a CCD camera as a gas detector, which is neither our case.

More recently, E. Llobet [16] reviewed the principal methods for dynamic analysis of the gas sensor response. Interestingly, the main use of these methods is to perform gas classification based on the transient response, from which a feature vector is extracted. For example, in [17] Box-Jenkins linear filters were applied to model an array of MOX sensors in the presence of four alcohols and water vapor with the aim of reducing the effect of the sensor drift in a classification process, and in [18] a multi-exponential transient spectroscopy (METS) method is proposed to improve the selectivity of chemical sensors in the analysis of gas mixtures.

Special mention deserves some works that rely, as it is our case, on modeling the sensor response to predict steady state values from the initial part of the transient. In [19], a so-called ARMA and multi-exponential models are proposed for reducing the time necessary to calibrate a sensor array. Nevertheless, since the focus is on the calibration of MOX sensors, the dynamic models are only applied to the rise transient signals recorded in Closed Sampling Systems over long time periods (over 800 s), while we aim to predict the gas distribution profile in real time and in OSS. In [20], A. Pardo et al. propose and compare different nonlinear inverse dynamic models of gas sensing systems for quantitative measurements. However, the considered dynamic conditions differ from those of OSS. First, a measurement chamber is used to obtain the gas sensor readings, which implicitly modifies the dynamic properties of the measured signals, and second, the acquisition frequency is too low (one sample per minute) to reflect the fast and highly dynamic changes of the gas concentration in OSS.

Based on the multi-exponential model proposed elsewhere [19,21], and taking into account the differences between rise and decay phases of MOX sensors, in this work we exploit a simplified version of it, where only one exponential is considered to model each phase (see Figure 2). According to this model, mathematically expressed in Equation (1), three phases can be considered in the output of a typical MOX gas sensor when exposed to an ideal step in concentration: baseline, gas measurement (rise) and recovery.

R ( t ) = { R 0 , t < t s R 0 + ( R max - R 0 ) ( 1 - e - ( t - t s ) τ r ) , t s < t < t e R 0 + ( R max - R 0 ) e - ( t - t s - Δ t ) τ d , t e < t
where τr and τd are the time constants for the rise and recovery phases respectively, ts and te represent the starting and ending times of the step excitation, R0 and R 0 are the sensor response level before and after the stimulus, Rmax is the saturation level, and R max is the maximum response level during the gas measurement phase. Notice that R max is usually lower than Rmax for short input pulses, as depicted in Figure 2.

This model was used by [22] for characterizing the response of an MOX-based e-nose carried by a robot. In this work we also exploit this model but making use of its inverted form, that is, to predict the gas distribution that the sensor is exposed to from its readings.

3. The Proposed MOX Model

As shown in the block diagram of Figure 3, three different sub-processes can be distinguished in an MOX-based gas sensing process:

  • A non-linear static system representing the measurement electronic circuit (Ω to V).

  • A transformation (Ω to Ω) that captures the non-linear rate at which the sensor resistance varies over time (although we have separated the transduction and dynamic phases to explicitly denote both functionalities, both stages take place within the MOX transducer and they are most likely coupled).

  • A signal transduction mechanism (ppm to Ω) which results from the chemical interaction between the sensor sensitive surface and the molecules of reducing gases.

Next, we model each of these stages and invert them to come up with a complete inverted MOX sensor model. It is important to remark that the proposed inverted model does not aim to recover the gas concentration (X[ppm]), which corresponds to a sensor calibration problem not addressed here, but only the gas distribution, providing relative results proportional to the gas concentration (R1).

3.1. Measurement Circuit

This stage stands for the electronic circuit in charge of measuring the changes in the sensor resistance, which typically consists of a simple voltage divider:

Y = V R L ( t ) = V CC × R L R 2 ( t ) + R L
where VCC is the circuit voltage and RL is the load resistance.

Since the magnitude measured is the output voltage Y (i), we can easily recover the MOX resistance changes R2(i) by solving Equation (3).

R 2 ( i ) = R L × ( V CC - Y ( i ) ) Y ( i )

3.2. Transient Behavior Stage

As reported by previous authors [7,22], the transient response of an MOX sensor can be expressed by two first-order systems, as depicted in Equation (1). It clearly resembles a low-pass filter response with the particularity that the filter cutoff frequency ( f c = 1 2 π τ ) is different for each phase the sensor is working at (rise or recovery). Equation (4) presents the two phases transfer function of this stage in the Laplace domain, where a common static delay (td) has additionally been considered in both phases to compensate for the delay introduced by the pneumatic circuit used to draw in the gas and flow it through the sensors:

R 2 ( s ) R 1 ( s ) = { A τ r s + 1 e - t d s , for rise phases A τ d s + 1 e - t d s , for recovery phases
where s is the Laplace variable, A is the filter gain, τrandτd are the filter time constants in the rise or recovery phase respectively, and td (in seconds) is the system delay. For nomenclature clarification τrecovery has been denoted as τd, where the subindex d stands for decay.

Since our interest remains in estimating the volatile distribution that the sensor is being exposed to, given the sensor readings, we work out R1 as a function of the sensor resistance measurements R2 by applying the inverse Laplace transform to Equation (4). As the transfer function has different expressions according to the phase (rise or recovery) that the sensor is working at, the parameters of the resulting differential equation will have to switch accordingly. Approximating the derivative of the measured sensor resistance R 2 by a backward first order finite difference, such differential equation can be written as:

R 1 ( i - N ) R 2 ( i ) + τ R 2 ( i ) - R 2 ( i - 1 ) Δ t
where R1(i) is the unknown steady resistance value for the given gas concentration at the time step i, R2 is the measured sensor resistance, τ is the time constant for either the rise (τr) or the recovery phase (τd), is the number of samples for the system delay, and Δt is the time between samples. Notice also that the scale factor A from Equation (4) has been dropped.

As observed in Equation (5), at each time step this dynamic model requires to know the value of parameters τr and τd, and the phase the sensor is working at. The latter may be determined from the slope of the measured MOX resistance R2, considering that the sensor is working under rise phase for positive values of the derivative, and under decay phase for negative values of it.

Notice that instead of working with the sensor resistance (which is inversely proportional to the volatile concentration), it may be more convenient to deal with the sensor conductance (GMOX = 1/RMOX), which is proportional to the gas concentration.

For estimating the values of τr and τd, a common practice in system modelling is that of identifying the parameters of the system upon its step response [23]. The problem with such procedure is that the time constants, in practice, depend not only on the type of volatile but also on its concentration [20]. In this paper this dependency is made explicit by adjusting a polynomial regression model over a sequence of concentration pulses with different amplitudes for each target gas. In Section 5.1, an example of such relation is depicted for the target gas ethanol.

3.3. Transduction Stage

The transduction stage is commonly defined by the sensitivity characteristics and the temperature and humidity dependencies of the transducer. For the case of Closed Sampling Systems, those characteristics are usually provided by the sensor manufacturer, relating the volatile concentration [ppm] to the sensor resistance ratio RS/R0 (sensor resistance in gas over sensor resistance in air) for different target gases and test conditions. Nevertheless, those sensitivity characteristics are obtained by measuring steady state values of the MOX sensor resistance after very long and constant exposure times, which are not applicable to Open Sampling Systems and thus not considered in this paper.

4. Signal Conditioning and Preprocessing

The proposed model-based approach given by Equation (3) and Equation (5) relies on the sensor readings to obtain an estimation of the gas. As can be appreciated, a first order derivative needs to be computed to obtain such estimation, which notably degrades the signal-to-noise ratio and consequently the accuracy in the estimation (see Figure 4). Additionally, MOX sensors are susceptible to long and short term drift [5], gradually changing the sensor resistance even if exposed to exactly the same gas concentration under identical environmental conditions.

It then becomes necessary to carry out a signal conditioning to prepare the sensor readings to the posterior estimation process. Initially, for the purpose of drift compensation and dynamic range enhancement, the raw sensor readings R2(i) are divided by the sensor baseline resistance at t = 0, that is, R2(0). This transformation is known as relative baseline manipulation [5]. Later, in order to mitigate the noise effects on the model, a low pass filter followed by a sub-sampling process are applied to the signal, as depicted in Figure 4. The cutoff frequency of the filter and the down-sampling rate have been determined experimentally according to the sampling frequency.

5. Experimental Results

This section presents three different experiments designed with increasing complexity to test how the proposed model can anticipate the steady state values of the sensor resistance from transient measurements in Open Sampling Systems. We start by testing our approach in a scenario where airflow and volatile distribution were well controlled. Then, two experiments of gas distribution mapping in unmodified environments are described.

5.1. Train of Gas Pulses in a Controlled Scenario

This experiment is designed to validate, in the simplest possible way, the ability of the proposed approach, given by Equation (3) and Equation (5), to estimate the volatile distribution that the MOX gas sensors are being exposed to.

Since knowing the ground truth distribution of a volatile in an OSS is a tough task (we can even say utopian), a specific setup was designed to keep its concentration as constant as possible and to confine it in a predefined region, avoiding the dispersion of gas particles to undesired locations. The experiment setup is depicted in Figure 5. The gas source, composed by a small cup filled with acetone, was placed inside a cardboard box with a small upper aperture. A fan located at the bottom of the box (beneath the gas source) was used to generate a constant upstream airflow pushing the volatile through the chamber aperture, while a second fan, placed about 20 cm over the box, sucked up the air from the box. This setup allows us to keep the volatile confined in an approximated gas column between both fans.

A Patrolbot robotic base [24] was then placed in front of the gas source and commanded to rotate with constant angular speed ω. The e-nose, placed at the right side of the mobile base, then described a circular trajectory being exposed to the volatile only when passing through the gas column. The resulting excitation signal was then a train of identical duration pulses. The position of the source at each lap (and consequently the duration of each excitation pulse) was precisely obtained from the radial scan provided by a laser range scanner (SICK LMS500) carried by the Patrolbot.

Strictly speaking, because the sensor needs some time to entirely enter into the odor column, the excitation signal will not be a perfect pulse, but the deviation is small and can be neglected.

The pre-calibration of the time constants with the gas concentration, as depicted in Section 3.2, has been achieved by a polynomial regression over more than 50 short pulses with different amplitudes of the same volatile, in this case ethanol. The resulting dependency is depicted in Figure 6. Please, note that since the gas concentration is not available, we are considering the sensor conductance (1/R2) instead.

Figure 7 presents the results of this experiment for a robot angular speed of 30°/s, that is, the e-nose is exposed to the same pulse every 12 s. The width of the pulse (in seconds) is computed from the robot rotational speed and the angular references detected from the laser scan. The sequence of rise and recovery phases in the sensor output (solid red line) can be observed. It is also clear how the long tails of the recovery phases are interrupted by the next rising phase. (This phase overlap produces an interesting accumulation effect in the sensor output, increasing the sensor reading value on each exposition to the gas source even though the volatile concentration was constant along the experiment duration. The study of this effect is out of the scope of this work.)

As can be noticed, the estimation of the gas distribution (dashed blue line) as provided by applying the proposed model is a more accurate estimation of the train of pulses than the raw sensor readings. The improvement is significant for the recovering phases, not only because the overlap between phases is avoided (the recovery phase is interrupted by the next rise phase before reaching the baseline), but because our estimation provides values of the gas distribution more consistent with the reality.

5.2. 1D Gas Distribution Mapping

The motivation of this experiment is to demonstrate the utility of the proposed approach in the generation of gas distribution maps, a challenging problem in robotic olfaction [25,26].

As described in the introduction, the slow recovery of MOX sensors becomes a serious drawback when the sensors are carried on a vehicle (typically a mobile robot). In such cases, the adopted solution is to reduce the vehicle velocity to a few cm/s, increasing considerably the execution time of the task. We demonstrate here that our estimation of the gas distribution to a certain extent removes the need to reduce the robot velocity for maintaining the map accuracy.

For such a goal, two different gas distribution maps of an unmodified long indoor corridor are generated by driving the robot at two different velocities. Figure 8 shows a comparison of the gas maps generated using the Kernel-based method proposed in [27] for a reduced robot speed of 0.1 m/s, while Figure 9 represents the results of a similar experiment using a robot speed of 0.4 m/s.

In both cases, the robot travels the corridor twice (round trip), passing over a cap filled with acetone placed in the middle of it. The map evolves as the robot moves and collects new sensory data, so we decided to analyze the map as built, at two different points of the robot path: at the end of the corridor (top row sub-figures), and at the initial position after the task has been completed (bottom row sub-figures).

As can be seen from the results, the “tail” effect produced by the slow recovery of MOX sensors leads to a gas distribution map with high concentrations along the robot path once the source has been hit. As expected, this effect is more harmful when increasing the robot speed, one of the reasons why the speed of olfactory robots is usually kept small. Please, note how this effect is substantially palliated when applying the proposed gas estimation (see Figure 9(e) and Figure 9(f)), even when the robot speed is high.

As a conclusion, our approach not only presents a noticeable improvement in the correct localization of the gas source but also provides confident values of the gas distribution along the robot path. This is important since in most cases it is desirable to not only know the location of the gas source but also know how the volatile emanating from it has spread in the surroundings.

5.3. 2D Gas Distribution Mapping

In this last experiment, we pursue to test our model-based estimation of the gas concentration in a more complex environment. The idea is to accomplish a complete olfaction mission in a realistic scenario where no alterations of the environment are done. Typical examples are the localization of leaks or the declaration of areas of high concentration levels of harmful gases.

The testing scenario consisted of two adjacent rooms communicated through a small corridor (see Figure 10). The robot inspected both rooms following the predefined path marked as solid red line, at a speed of 0.3m/s. The gas source, composed by a cardboard plate impregnated in ethanol, was placed on the floor along the robot path. To be able to compare the maps generated from the raw sensor readings with those provided after applying the proposed MOX sensor model, doors and windows were kept closed to avoid uncontrolled airflows.

Figures 11 and 12 depict the mean and predictive variance gas distribution maps generated by the KernelDM+V algorithm [2], when fed with the MOX sensor readings and the model-based gas estimation, respectively.

Focusing on the mean maps, it can be seen that in both cases the maximum concentration falls near the real source location (marked as a white circle). However, when the map is built from the raw sensor readings, a high concentration area appears along the corridor connecting both rooms. The shape reassembles that of a gas plume generated by a dominant airflow, which we know is not possible since doors and windows were kept closed during the experiment. This “fake” plume is attributed to the actual robot path and the slow recovery of MOX sensors. When building the map from the gas concentration estimated by the proposed model (Figure 12), this “fake” plume does not appear, correctly representing the gas distribution in the inspected area.

Focusing now on the predictive variance maps, it has been previously reported that the variance of a set of gas concentration measurements has been suggested as a feature that can identify the location of a source of gas [27,28]. Following this principle, it can be seen that the gas source cannot be precisely located in the first case (Figure 11) since there are high variance values not only nearby the real location of the source but also along the corridor. Again, the application of the proposed MOX model leads to a predictive variance map where all high variance values fall near the real source location, allowing a correct estimation of it.

6. Conclusions and Future Works

In this paper, a modeling approach has been reported for the improvement of open sampling system olfaction applications based on MOX gas sensors. The proposed approach compensates the slow dynamic behavior of the MOX sensor, by forecasting the steady state values of the sensor resistance from a sequence of transient measurements. The exploited model is based on two first order systems (rise and recovery) with time constants that depend on the sensor reading amplitude. Our approach has been validated in different scenarios, demonstrating its utility in applications like gas source localization or gas distribution mapping. Additionally, we have proved that our approach enables a considerable increase in the speed at which a mobile base carrying the MOX-based e-nose can inspect the environment, which directly implies an important reduction in the execution times of the olfaction task.

Future work includes modeling additional parameters as temperature and humidity in order to improve the gas distribution estimation in real complex scenarios. Besides, the study of the transduction stage (calibration) for open sampling systems will be considered to complete the MOX sensor model.

Acknowledgments

This work was partly supported by the Regional Government of Andalucía under research contract P08-TEP-4016.

References

  1. Trincavelli, M.; Coradeschi, S.; Loutfi, A. Online Classification of Gases for Environmental Exploration. Proceedings of the IEEE /RSJ International Conference on Intelligent Robots and Systems (IROS), St. Louis, MO, USA, 11– 15 October 2009; pp. 3311–3316.
  2. Lilienthal, A.J.; Reggente, M.; Trincavelli, M.; Blanco, J.L.; Gonzalez, J. A Statistical Approach to Gas Distribution Modelling with Mobile Robots The Kernel DM+V Algorithm. Proceedings of the IEEE /RSJ International Conference on Intelligent Robots and Systems (IROS), St. Louis, MO, USA, 11– 15 October 2009; pp. 570–576.
  3. Rachkov, M.; Marques, L.; De Almeida, A. Multisensor demining robot. Autono. Robots 2005, 18, 275–291. [Google Scholar]
  4. Trincavelli, M.; Reggente, M.; Coradeschi, S.; Ishida, H.; Loutfi, A.; Lilienthal, A.J. Towards environmental monitoring with mobile robots. Proceedings of the IEEE /RSJ International Conference on Intelligent Robots and Systems (IROS), Nice, France, 22– 26 September 2008; pp. 2210–2215.
  5. Pearce, T., Schiffman, S., Nagle, H., Gardner, J., Eds.; Handbook of Machine Olfaction, Electronic Nose Technology; Wiley-VCH: Berlin, Germany, 2003.
  6. Moseley, P., Tofield, B., Eds.; Solid State Gas Sensors; Adam Hilger: Bristol, UK, 1987.
  7. Gonzalez-Jimenez, J.; Monroy, J.G.; Blanco, J.L. The multi-camber electronic nose—An improved olfaction sensor for mobile robotics. Sensors 2011, 11, 6145–6164. [Google Scholar]
  8. Albert, K.; Lewis, N.; Schauer, C.; Sotzing, G.; Stitzel, S.; Vaid, T.; Walt, D. Cross reactive chemical sensor arrays. Chem. Rev. 2000, 100, 2595–2626. [Google Scholar]
  9. Sberveglieri, G. Classical and novel techniques for the preparation of SnO2 thin-film gas sensors. Sens. Actuators B: Chem. 1992, 6, 239–247. [Google Scholar]
  10. Shraiman, B.I.; Siggia, E.D. Scalar Turbulence. Nature 2000, 405, 639–646. [Google Scholar]
  11. Lilienthal, A.J.; Wandel, M.R.; Weimar, U.; Zell, A. Experiences Using Gas Sensors on an Autonomous Mobile Robot. Proceedings of 4th European Workshop on Advanced Mobile Robots, EUROBOT 2001, Lund, Sweden, 19– 21 September 2001; pp. 1–8.
  12. Trincavelli, M.; Coradeschi, S.; Loutfi, A. Odour classification system for continuous monitoring applications. Sens. Actuators B: Chem. 2009, 58, 265–273. [Google Scholar]
  13. Ishida, H.; Suetsugu, K.; Nakamoto, T.; Moriizumi, T. Study of autonomous mobile sensing system for localization of odor source using gas sensors and anemometric sensors. Sens. Actuators A: Phys. 1994, 45, 153–157. [Google Scholar]
  14. Gardner, J.W. A non-linear diffusion-reaction model of electrical conduction in semiconductor gas sensors. Sens. Actuators B: Chem. 1990, 1, 166–170. [Google Scholar]
  15. Yamanaka, T.; Ishida, H.; Nakamoto, T.; Moriizumi, T. Analysis of gas sensor transient response by visualizing instantaneous gas concentration using smoke. Sens. Actuators A: Phys. 1998, 69, 77–81. [Google Scholar]
  16. Lobet, E. Dynamic Pattern Recognition Methods and System Identification. In Handbook of Machine Olfaction; Pearce, T., Schiffman, S., Nagle, H., Gardner, J., Eds.; Wiley-VCH: Weinheim, Germany, 2000; Chapter 12; pp. 293–321. [Google Scholar]
  17. Holmberg, M.; Winquist, F.; Lundstrom, I.; Davide, F.; DiNatale, C.; DAmico, A. Drift counteraction for an electronic nose. Sens. Actuators B: Chem. 1996, 36, 528–535. [Google Scholar]
  18. Samitier, J.; López-Villegas, J.; Marco, S.; Cámara, L.; Pardo, A.; Ruiz, O.; Morante, J. A new method to analyze signal transients in chemical sensors. Sens. Actuators B: Chem. 1994, 18, 308–312. [Google Scholar]
  19. Natale, C.D.; Marco, S.; Davide, F.; D'Amico, A. Sensor-array calibration time reduction by dynamic modelling. Sens. Actuators B: Chem. 1995, 25, 578–583. [Google Scholar]
  20. Pardo, A.; Marco, S.; Samitier, J. Nonlinear inverse dynamic models of gas sensing systems based on chemical sensor arrays for quantitative measurements. IEEE Trans. Instrum. Meas. 1998, 47, 644–651. [Google Scholar]
  21. Gutierrez-Osuna, R.; Nagle, H.T.; Schiffman, S.S. Transient response analysis of an electronic nose using multi-exponential models. Sens. Actuators B: Chem. 1999, 61, 170–182. [Google Scholar]
  22. Lilienthal, A.J. Gas Distribution Mapping and Gas Source Localisation with a Mobile Robot. Ph.D. Thesis, Wilhelm-Schickard Institute, University of Tübingen, Tübingen, Germany, 2004. [Google Scholar]
  23. Ogata, K. Modern Control Engineering; Prentice-Hall Electrical Engineering Series. Instrumentation and Controls Series; Prentice Hall: Englewood Cliffs, NJ, USA, 2010. [Google Scholar]
  24. MobileRobots Inc. Available online: http://www.mobilerobots.com (accessed on 10 October 2012).
  25. Loutfi, A.; Coradeschi, S.; Lilienthal, A.J.; Gonzalez, J. Gas distribution mapping of multiple odour sources using a mobile robot. Robotica 2009, 27, 311–319. [Google Scholar]
  26. Lilienthal, A.J.; Streichert, F.; Zell, A. Model-based Shape Analysis of Gas Concentration Gridmaps for Improved Gas Source Localisation. Proceedings of the IEEE International Conference on Robotics and Automation (ICRA), Barcelona, Spain, 18– 22 April 2005; pp. 3575–3580.
  27. Lilienthal, A.J.; Duckett, T. Building Gas Concentration Gridmaps with a Mobile Robot. Robot. Auton. Syst. 2004, 48, 3–16. [Google Scholar]
  28. Hernandez Bennets, V.; Lilienthal, A.J.; Neumann, P.; Trincavelli, M. Mobile robots for localizing gas emission sources on landfill sites: Is bio-inspiration the way to go? Front. Neuroeng. 2012, 4, 20. [Google Scholar]
Figure 1. Rise and recovery phases of MOX sensor response to a step gas concentration. Figure (a) shows a 2D plot of the sensor response over time. The shaded blue region denotes the sensor exposure to the analyte. Figure (b) depicts a 3D gas distribution map generated from the reading of an MOX sensor carried by a mobile robot along a corridor. Observe how the recovery phase after the gas exposure is several time longer than the rise one.
Figure 1. Rise and recovery phases of MOX sensor response to a step gas concentration. Figure (a) shows a 2D plot of the sensor response over time. The shaded blue region denotes the sensor exposure to the analyte. Figure (b) depicts a 3D gas distribution map generated from the reading of an MOX sensor carried by a mobile robot along a corridor. Observe how the recovery phase after the gas exposure is several time longer than the rise one.
Sensors 12 13664f1 1024
Figure 2. Ideal response of an MOX sensor (solid red line) when excited with a step gas concentratration (dashed blue line). The curve shows the three phases of a measurement: (I) baseline, (II) GAS measurement, and (III) recovery phase.
Figure 2. Ideal response of an MOX sensor (solid red line) when excited with a step gas concentratration (dashed blue line). The curve shows the three phases of a measurement: (I) baseline, (II) GAS measurement, and (III) recovery phase.
Sensors 12 13664f2 1024
Figure 3. Block diagram of a smelling process with an MOX gas sensor. The sensor is excited by a volatile [ppm] producing a variation in the sensor resistance that is measured as an electrical signal [V] by means of a measurement circuit.
Figure 3. Block diagram of a smelling process with an MOX gas sensor. The sensor is excited by a volatile [ppm] producing a variation in the sensor resistance that is measured as an electrical signal [V] by means of a measurement circuit.
Sensors 12 13664f3 1024
Figure 4. Noisy sensor conductance readings and its derivative (solid red line), and the corresponding filtered versions (dashed blue line).
Figure 4. Noisy sensor conductance readings and its derivative (solid red line), and the corresponding filtered versions (dashed blue line).
Sensors 12 13664f4 1024
Figure 5. Picture of the Patrolbot mobile robot, the e-nose and the scheme of the gas source setup.
Figure 5. Picture of the Patrolbot mobile robot, the e-nose and the scheme of the gas source setup.
Sensors 12 13664f5 1024
Figure 6. Pre-calibration of the time constants τr and τd with the gas concentration (i.e., sensor conductance). Red circles and green squares represent the values of τr and τd respectively, while the lines represent their liner regression.
Figure 6. Pre-calibration of the time constants τr and τd with the gas concentration (i.e., sensor conductance). Red circles and green squares represent the values of τr and τd respectively, while the lines represent their liner regression.
Sensors 12 13664f6 1024
Figure 7. Results of the experiment in the controlled scenario for a robot angular speed ω = 30°/s. The train of acetone pulses (ground truth) is plotted as green shaded bars, the normalized raw sensor readings (Y) (See Figure 3) are plotted as a solid red line and the gas distribution estimated by the proposed approach is plotted as a dashed blue line.
Figure 7. Results of the experiment in the controlled scenario for a robot angular speed ω = 30°/s. The train of acetone pulses (ground truth) is plotted as green shaded bars, the normalized raw sensor readings (Y) (See Figure 3) are plotted as a solid red line and the gas distribution estimated by the proposed approach is plotted as a dashed blue line.
Sensors 12 13664f7 1024
Figure 8. Snapshots of the maps generated in the 1D gas distribution mapping experiment at the middle and end of the robot path (left/right column, respectively), for a robot speed of 0.1m/s. Top row sub-figures illustrate a comparison between the normalized sensor readings and the gas distribution estimated here. Middle and bottom row sub-figures show a 3D reconstruction of the maps generated with the raw MOX readings (middle row), or the estimation provided by the proposed model (bottom row).
Figure 8. Snapshots of the maps generated in the 1D gas distribution mapping experiment at the middle and end of the robot path (left/right column, respectively), for a robot speed of 0.1m/s. Top row sub-figures illustrate a comparison between the normalized sensor readings and the gas distribution estimated here. Middle and bottom row sub-figures show a 3D reconstruction of the maps generated with the raw MOX readings (middle row), or the estimation provided by the proposed model (bottom row).
Sensors 12 13664f8 1024
Figure 9. Snapshots of the maps generated in the 1D gas distribution mapping experiment at the middle and end of the robot path (left/right column, respectively), for a robot speed of 0.4m/s. Top row sub-figures illustrate a comparison between the normalized sensor readings and the gas distribution estimated here. The middle and bottom row sub-figures show a 3D reconstruction of the maps generated with the raw MOX readings (middle row), or the estimation provided by the proposed model (bottom row).
Figure 9. Snapshots of the maps generated in the 1D gas distribution mapping experiment at the middle and end of the robot path (left/right column, respectively), for a robot speed of 0.4m/s. Top row sub-figures illustrate a comparison between the normalized sensor readings and the gas distribution estimated here. The middle and bottom row sub-figures show a 3D reconstruction of the maps generated with the raw MOX readings (middle row), or the estimation provided by the proposed model (bottom row).
Sensors 12 13664f9 1024
Figure 10. Map of the experimental area used in the 2D mapping experiment. Blue points represent obstacles detected by the onboard SICK laser scan, the robot path is marked as red solid line, and the gas source position is pointed with a green circle.
Figure 10. Map of the experimental area used in the 2D mapping experiment. Blue points represent obstacles detected by the onboard SICK laser scan, the robot path is marked as red solid line, and the gas source position is pointed with a green circle.
Sensors 12 13664f10 1024
Figure 11. Mean and predictive variance gas distribution maps of the inspected area generated by the KernelDM+V algorithm when fed with the sensor readings (after baseline manipulation and delay correction). The gas source location has been marked as a white circle.
Figure 11. Mean and predictive variance gas distribution maps of the inspected area generated by the KernelDM+V algorithm when fed with the sensor readings (after baseline manipulation and delay correction). The gas source location has been marked as a white circle.
Sensors 12 13664f11 1024
Figure 12. Mean and predictive variance s distribution maps of the inspected area generated by the KernelDM+V algorithm after applying the proposed model-base estimation of the gas concentration. The gas source location has been marked as a white circle.
Figure 12. Mean and predictive variance s distribution maps of the inspected area generated by the KernelDM+V algorithm after applying the proposed model-base estimation of the gas concentration. The gas source location has been marked as a white circle.
Sensors 12 13664f12 1024

Share and Cite

MDPI and ACS Style

Monroy, J.G.; González-Jiménez, J.; Blanco, J.L. Overcoming the Slow Recovery of MOX Gas Sensors through a System Modeling Approach. Sensors 2012, 12, 13664-13680. https://doi.org/10.3390/s121013664

AMA Style

Monroy JG, González-Jiménez J, Blanco JL. Overcoming the Slow Recovery of MOX Gas Sensors through a System Modeling Approach. Sensors. 2012; 12(10):13664-13680. https://doi.org/10.3390/s121013664

Chicago/Turabian Style

Monroy, Javier G., Javier González-Jiménez, and Jose Luis Blanco. 2012. "Overcoming the Slow Recovery of MOX Gas Sensors through a System Modeling Approach" Sensors 12, no. 10: 13664-13680. https://doi.org/10.3390/s121013664

Article Metrics

Back to TopTop