<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ICA</journal-id><journal-title-group><journal-title>Intelligent Control and Automation</journal-title></journal-title-group><issn pub-type="epub">2153-0653</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ica.2013.41006</article-id><article-id pub-id-type="publisher-id">ICA-27704</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  Fault Detection and Isolation in Industrial Systems Based on Spectral Analysis Diagnosis
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>hmed</surname><given-names>Hafaifa</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Mouloud</surname><given-names>Guemana</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Attia</surname><given-names>Daoudi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Faculty of Science and Technology, University of Medea, Medea, Algeria</addr-line></aff><aff id="aff1"><addr-line>Electrical Engineering Field, Faculty of Science and Technology, University of Djelfa DZ, Djelfa, Algeria</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>hafaifa@hotmail.com(HH)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>05</day><month>02</month><year>2013</year></pub-date><volume>04</volume><issue>01</issue><fpage>36</fpage><lpage>41</lpage><history><date date-type="received"><day>November</day>	<month>9,</month>	<year>2012</year></date><date date-type="rev-recd"><day>December</day>	<month>21,</month>	<year>2012</year>	</date><date date-type="accepted"><day>December</day>	<month>30,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   The diagnoses in industrial systems represent an important economic objective in process industrial automation area. To guarantee the safety and the continuity in production exploitation and to record the useful events with the feedback experience for the curative maintenance. We propose in this work to examine and illustrate the application ability of the spectral analysis approach, in the area of fault detection and isolation industrial systems. In this work, we use a combined analysis diagram of time-frequency, in order to make this approach exploitable in the proposed supervision strategy with decision making module. The obtained results, show clearly how to guarantee a reliable and sure exploitation in industrial system, thus allowing better performances at the time of its exploitation on the supervision strategy. 
 
</p></abstract><kwd-group><kwd>Diagnosis; Spectral Analyzes; Faults Detection and Isolation; Condition Monitoring</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The development of diagnosis system appears additional expensive solution in the industrial investment, but it deadens during the production phase. During the execution of the strategy of diagnosis, the measurable signals can provide significant information, on the defects appearance to facilitate the determination of the supervision parameters, representing the defects signatures, their degree and persistence. From these parameters, the decisional system can conceive powerful diagnosis approach [1-3]. To realise this prospect, we proposes in this work to examine and illustrate the application ability of the spectral analysis approach, in the area of fault detection and isolation industrial systems. Thus it is essential to maintain the exploitation system apart from this instability zone. The signals temporal representation does not give a good sensitivity of fault detection on defective components, while the frequencies representation given by Fourier transform does not allow the temporal localization of these components [4-6]. In this work, we use a combined analysis diagram of time-frequency, in order to make this approach exploitable in the proposed supervision strategy with decision making module.</p><p>The obtained results, show clearly how to guarantee a reliable and sure exploitation of the studied industrial system, thus allowing better performances at the time of its exploitation on the supervision strategy. The application results of the investigated approach, validated by experimental test, allowing that the use of the combined time-frequency method give excellent results on the tested models with better performances. The proposed approach, allows the detection of the defects at the beginning of their appearance. The studied approach ensures an effective monitoring, allows detecting the defects to avoid any degradation and also to ensure an operation reliable and increase the working time of this industrial system.</p></sec><sec id="s2"><title>2. Diagnosis Based on Spectral Analysis Approach</title><p>The detection and the identification of the defects in the industrial systems it is an important subject in modern automation engineering companies [7-9]. Indeed, in large industrial applications it is often posed the question, how guarantee the people security and preserving their environment?</p><p>The diversity of the developed approaches for the diagnosis of the industrial systems, seem to be the result of various applications aimed by specifications which result from it. Thus, the nature of information available on the system or the defects type to be detected led to the implementation of this specific strategy [3,6,7]. In this work, we present the application of the spectral analysis method to the diagnosis of the complex industrial systems, based on Fourier transform; we consider the signal <img src="6-7900237\a59c2ce7-c618-4da6-a94e-fc214f079383.jpg" /> absolutely integrable on<img src="6-7900237\ebde9e80-13c2-4b3b-a6f7-4a337fa7c242.jpg" />, the transform of Fourier of this signal is to give by:</p><disp-formula id="scirp.27704-formula127642"><label>(1)</label><graphic position="anchor" xlink:href="6-7900237\a805ad96-5a13-491a-877f-3da4404cf089.jpg"  xlink:type="simple"/></disp-formula><p>This transform is invertible by:</p><disp-formula id="scirp.27704-formula127643"><label>(2)</label><graphic position="anchor" xlink:href="6-7900237\808b0dfd-b578-4555-b86b-34cb08158bfd.jpg"  xlink:type="simple"/></disp-formula><p>While applying this transform, we pass from the temporal field to the frequencies presentation, the energy is preserved in both representations. It is noted that the signal <img src="6-7900237\24f28969-e4dc-4f68-bbaa-806252cb2b7b.jpg" /> tends towards the zero when t tends towards the infinite, for practical reasons, this signal is sampling in the following form [<xref ref-type="bibr" rid="scirp.27704-ref8">8</xref>]:</p><disp-formula id="scirp.27704-formula127644"><label>(3)</label><graphic position="anchor" xlink:href="6-7900237\c489cbe4-f91d-4f80-abf4-1f35fccb69b4.jpg"  xlink:type="simple"/></disp-formula><p>With <img src="6-7900237\62d272f0-a0ea-421a-9f7b-709c3d119d27.jpg" /> is an approximation of <img src="6-7900237\b9f5ca55-8fc9-43e8-afb5-f40a41fc8b47.jpg" /> at the points<img src="6-7900237\cac8f75f-3872-4c8c-adec-194f08307709.jpg" />. The power density spectral is defined as being the square module of the Fourier transform of<img src="6-7900237\4a31ba8c-e94b-4fac-ac3f-16ff16171ac6.jpg" />.</p><p>In our investigation, we add a Gaussian white signal disturbed, as it’s shown on the <xref ref-type="fig" rid="fig1">Figure 1</xref>, this white vi-</p><p>bration is a realized with a random process, were the power density spectral is the same for all frequencies. This white vibration is a normal law of average and variance given by the following equation [<xref ref-type="bibr" rid="scirp.27704-ref10">10</xref>]:</p><disp-formula id="scirp.27704-formula127645"><label>(4)</label><graphic position="anchor" xlink:href="6-7900237\d0e64e8d-adc3-4a29-800d-289c56e849a7.jpg"  xlink:type="simple"/></disp-formula><p>In literary, there exist many methods to estimate the faults [4,5,11]. In our investigation, we chose the power density spectral, the estimation being carried out directly starting from the Fourier transform of the signal without concerned the models parametric in consideration. Under the hypothesis of stationary and ergodic signals, the function of autocorrelation of the power density spectral is defined by:</p><disp-formula id="scirp.27704-formula127646"><label>(5)</label><graphic position="anchor" xlink:href="6-7900237\d70f6fae-ab51-433f-8b53-b6667290c911.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.27704-formula127647"><label>(6)</label><graphic position="anchor" xlink:href="6-7900237\5302c07e-0de4-44d8-bb34-aa597e560414.jpg"  xlink:type="simple"/></disp-formula><p>One supposes we have N samples <img src="6-7900237\76bbf8b3-20a1-4300-968e-6cd531171b44.jpg" /> of the signal to be analyzed and we thus seek to consider the power density spectral starting from the measured data. Using a nonparametric estimate of the spectrum, each method being related to the equalities in the Equaion (6). The representation of the power density spectral using the FFT is shown on <xref ref-type="fig" rid="fig2">Figure 2</xref>, this resolution method is nominal one Hz equal to the inverse number of multiplies sampling by the sampling frequency; this resolution decreases the important noise in our system.</p><sec id="s2_1"><title>2.1. Estimation of the Spectrum by the Periodogram Method</title><p>Now we will use the spectrum estimation by the method of periodogram [<xref ref-type="bibr" rid="scirp.27704-ref12">12</xref>], in this method we used the signal directly. The spectrum is estimated by the following relation:</p><disp-formula id="scirp.27704-formula127648"><label>(7)</label><graphic position="anchor" xlink:href="6-7900237\62ae611e-594e-4d96-b684-46e0a06ae7d3.jpg"  xlink:type="simple"/></disp-formula><p>The periodogram is the convolution of the spectrum by a window in cardinal sine, there resolution is about<img src="6-7900237\40624521-a051-4672-bd5c-e4adbf9b8ef5.jpg" />the periodogram is a skewed estimator, given by [<xref ref-type="bibr" rid="scirp.27704-ref13">13</xref>]:</p><disp-formula id="scirp.27704-formula127649"><label>(8)</label><graphic position="anchor" xlink:href="6-7900237\96531959-84bb-415f-9aa8-cafcbaadaa1c.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.27704-formula127650"><label>(9)</label><graphic position="anchor" xlink:href="6-7900237\06843abd-8d6e-4aaa-b567-6564242e9744.jpg"  xlink:type="simple"/></disp-formula><p>This periodogram representation is the triangular window of Fourier transform, when<img src="6-7900237\908e095e-56fb-4a32-8e6f-17ef1a789f3f.jpg" />, the presentation of the power spectrum using the periodogram method, is shown on the <xref ref-type="fig" rid="fig3">Figure 3</xref>, the skew becomes null, that show us that the variance is practically independent of N and proportional to the spectrum, given by [<xref ref-type="bibr" rid="scirp.27704-ref6">6</xref>]:</p><disp-formula id="scirp.27704-formula127651"><label>(10)</label><graphic position="anchor" xlink:href="6-7900237\be52778e-2ff5-43b5-a201-ad9b3f9cbf2b.jpg"  xlink:type="simple"/></disp-formula><p>The periodogram is not consistent an estimator for the DSP, to decrease the variance of this estimator, we can used the periodogram modify. It is possible to balance the estimator of the autocorrelation function by a window function <img src="6-7900237\e4ca4942-b817-4baf-ab44-ce445b8f55ad.jpg" /> on N measures. The spectral estimator must be calculated as follows:</p><disp-formula id="scirp.27704-formula127652"><label>(11)</label><graphic position="anchor" xlink:href="6-7900237\28961485-c065-43e3-9766-6ff021f555e1.jpg"  xlink:type="simple"/></disp-formula><p>Using a normalized function U to eliminate the skew introduced by the window function<img src="6-7900237\4e86f649-b82f-4571-9ec6-a55f21ab2af9.jpg" />, given by:</p><disp-formula id="scirp.27704-formula127653"><label>(12)</label><graphic position="anchor" xlink:href="6-7900237\0f7cbffb-45ae-447f-9ca2-44de90035c4b.jpg"  xlink:type="simple"/></disp-formula><p>This new estimator is as skewed by comparison with the precedent; the applied window function to the samples is obtained before calculation of FFT. His advantage is reduced its variance, because the effect of smoothing of the power density spectral is realized by the convolution with the window function. The representation of the power density spectral, by using of the periodogram modify method is shown on <xref ref-type="fig" rid="fig4">Figure 4</xref>. The choice of the window function becomes important, when the power density spectral of the studied signal includes superimposed lines of broad band noise, like the vibration in our application.</p></sec><sec id="s2_2"><title>2.2. Power Density Spectral Estimated by PRONY</title><p>The proposed approach of detection and localisation of defects using PRONY [<xref ref-type="bibr" rid="scirp.27704-ref12">12</xref>] method is an exponential sum of the probability defects model to approach at the follows<img src="6-7900237\5695a71f-8096-4700-bd1f-3ad244a29f1e.jpg" />, given by:</p><disp-formula id="scirp.27704-formula127654"><label>(13)</label><graphic position="anchor" xlink:href="6-7900237\5d4c43c4-17a9-4b09-ad08-b02312cb85e5.jpg"  xlink:type="simple"/></disp-formula><p>In this model, the exponential <img src="6-7900237\2fb3d74d-2d27-49da-91c1-fd8ad6b2c339.jpg" /> noted the poles of model are carrying the information of the nature of defects, known in signal treatment and in system modelling by:</p><p><img src="6-7900237\f7da4ada-f24d-4806-aae8-add0dca23d38.jpg" />and <img src="6-7900237\d37ebe11-dcba-46ed-96a0-63663da74fc6.jpg" />&#160;&#160;&#160; (14)</p><p>With <img src="6-7900237\c96ca338-9baf-4417-a4b4-48cc6f46c73e.jpg" /> are respectively the amplitude and the phase of the damping vector and <img src="6-7900237\45658791-3ad8-478d-b627-4aef126e4720.jpg" /> is the frequency and the step of sampling. To determine the frequencies and the damping coefficients, in our application we have used the following equations:</p><disp-formula id="scirp.27704-formula127655"><label>(15)</label><graphic position="anchor" xlink:href="6-7900237\c41a9500-a5ec-49e2-bee8-d7e3c5ef5996.jpg"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.27704-formula127656"><label>(16)</label><graphic position="anchor" xlink:href="6-7900237\bc1b8bfe-6834-48b6-9fcb-a5b0846d141c.jpg"  xlink:type="simple"/></disp-formula><p>For the determination of the amplitudes and the phases we based on the Equations (15) and (16):</p><disp-formula id="scirp.27704-formula127657"><label>(17)</label><graphic position="anchor" xlink:href="6-7900237\9e2847a6-c81c-4bfc-99d6-a604e932e49c.jpg"  xlink:type="simple"/></disp-formula><p>The approximation function becomes then:</p><disp-formula id="scirp.27704-formula127658"><label>(18)</label><graphic position="anchor" xlink:href="6-7900237\910c7946-6bee-4589-a817-09833bf3b199.jpg"  xlink:type="simple"/></disp-formula><p>The power density spectral estimated by PRONY is then given by:</p><disp-formula id="scirp.27704-formula127659"><label>(19)</label><graphic position="anchor" xlink:href="6-7900237\a150e6f3-c335-41e3-9f9a-de253eae446a.jpg"  xlink:type="simple"/></disp-formula><p>The spectral analysis consists in raising the fault signal measures on the examined system and to proceed has a systematic analysis, to seek the presence of fault situation of the whole of the defects suitable for affect the installation considered. That gives access to the diagnosis; to identify with precision the nature of the anomaly and to specify revolves it with a good localization of these defects. The representation of the power density spectral with the use of Prony method is shown on the Figures 5(a)-(c).</p></sec></sec><sec id="s3"><title>3. Conclusion</title><p>In the presented results in this paper, we have showed that the use of the analysis spectral method, the power of the signals models using Prony method, these responses shows that the technique of Prony makes it possible to distinguish the frequencies <img src="6-7900237\dc09858c-92e5-4bb7-93bd-a6aabe8ebe5a.jpg" /> and <img src="6-7900237\62f6b926-8602-49ec-97d0-6eaa3b638e95.jpg" /> in spite of their close values, and in the presence of satisfactory noise signal<img src="6-7900237\7f525826-e568-47dd-aeb5-552acee8d76e.jpg" />. For our investigation the starting and stopping control of the studied sytsem is supervised by the distributed control system, the fault measurements in the used system is done on-line, these measurements are taken periodically by using the various measurements sensors.</p></sec><sec id="s4"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.27704-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">C. M. Xia, J. Howell and N. F. Thornhill, “Detecting and Isolating Multiple Plant-Wide Oscillations via Spectral Independent Component Analysis,” Automatica, Vol. 41, No. 12, 2005, pp. 2067-2075.  
doi:10.1016/j.automatica.2005.02.011</mixed-citation></ref><ref id="scirp.27704-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">O. Coutier-Delgosha, R. Fortes-Patella and J. L. Reboud, “Evaluation of the Turbulence Model Influence on the Numerical Simulations of Unsteady Cavitation,” Journal of Fluids Engineering, Vol. 125, No. 1, 2003, pp. 38-45.  
doi:10.1115/1.1524584</mixed-citation></ref><ref id="scirp.27704-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">A. Hafaifa, A. Daoudi and K. Laroussi, “Application of Fuzzy Diagnosis in Fault Detection and Isolation to the Compression System Protection,” Control and Intelligent Systems ACTA Press, Vol. 39, No. 3, 2011, pp. 151-158.</mixed-citation></ref><ref id="scirp.27704-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">A. Hafaifa, F. Laaouad and K. Laroussi, “A Numerical Structural Approach to Surge Detection and Isolation in Compression Systems Using Fuzzy Logic Controller,” International Journal of Control, Automation, and Systems, Vol. 9, No. 1, 2011, pp. 69-79.  
doi:10.1007/s12555-011-0109-3</mixed-citation></ref><ref id="scirp.27704-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Y. Tsujimoto, Y. Yoshida, Y. Maekawa, S. Watanabe, and T. Hashimoto, “Observations of Oscillating Cavitation of an Inducer,” Journal of Fluids Engineering, Vol. 119, No. 1, 1997, pp. 775-781. doi:10.1115/1.2819497</mixed-citation></ref><ref id="scirp.27704-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">J. L. Weiss, “Threshold Computations for Detection of Failures in SISO Systems with Transfer Function Errors,” Proceedings of the American Control Conference, Atlanta, 15-17 June 1988, pp. 2213-2218.</mixed-citation></ref><ref id="scirp.27704-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">A. S. Willsky, “A Survey of Design Methods for Failure Detection in Dynamic Systems,” Automatica Journal, Vol. 12, 1976, pp. 601-611.  
doi:10.1016/0005-1098(76)90041-8</mixed-citation></ref><ref id="scirp.27704-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">A. Hafaifa, F. Laaouad and K. Laroussi, “Fuzzy Logic Approach Applied to the Surge Detection and Isolation in Centrifugal Compressor,” Automatic Control and Computer Sciences, Vol. 44, No. 1, 2010, pp. 53-59.  
doi:10.3103/S0146411610010074</mixed-citation></ref><ref id="scirp.27704-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">M. Guemana, S. Aissani and A. Hafaifa, “Use a New Calibration Method for Gas Pipelines: An Advanced Method Improves Calibrating Orifice Flowmeters While Reducing Maintenance Costs,” Hydrocarbon Processing Journal, Vol. 90, No. 8, 2011, pp. 63-68.</mixed-citation></ref><ref id="scirp.27704-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">D. Van Schrick, “Estimator Schemes for Instrument Fault Detection and Isolation,” IEEE Proceeding of the International Conference on System, Man and Cybernetic, Le-Touquet, 17-20 October 1993, pp. 406-411.  
doi:10.1109/ICSMC.1993.385045</mixed-citation></ref><ref id="scirp.27704-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">B. K. Walker and E. Gai, “Fault Detection Threshold Determination Techniques Using Markov Theory,” International Journal Guidance, Control and Dynamics, Vol. 2, No. 4, 1979, pp. 313-319. doi:10.2514/3.55881</mixed-citation></ref><ref id="scirp.27704-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">P. Weber, S. Gentil, P. Ripoll and L. Foulloy, “Multiple Fault Detection and Isolation,” Proceeding of the 14th IFAC World Congress, Vol. P, No. 7e-092, 1999, pp. 223-228.</mixed-citation></ref><ref id="scirp.27704-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">H. L. Jiang, M. A. A. Shoukat Choudhury and L. Sirish Shah, “Detection and Diagnosis of Plant-Wide Oscillations from Industrial Data Using the Spectral Envelope Method,” Journal of Process Control, Vol. 17, No. 2, 2007, pp. 143-155. doi:10.1016/j.jprocont.2006.09.006</mixed-citation></ref></ref-list></back></article>