CC BY
56
International Journal of Open Information Technologies ISSN: 2307-8162 vol. 3, no. 9, 2015
Towards Events Recognition in a Distributed Fiber-Optic Sensor System: Kolmogorov-Zurbenko Filtering
Aleksey Fedorov, Maxim Anufriev, Andrey Zhirnov, Evgeniy Nesterov, Dmitry Namiot, Alexey Pnev, and Valery Karasik
Abstract—The paper is about de-noising procedures aimed on events recognition in signals from a distributed fiber-optic vibration sensor system based on the phase-sensitive optical time-domain reflectometry. We report experimental results on recognition of several classes of events in a seismic background. A de-noising procedure uses the framework of the time-series analysis and Kolmogorov-Zurbenko filtering. We demonstrate that this approach allows revealing signatures of several classes of events.
Keywords—de-noising, optical reflectometry, Kolmogorov-Zurbenko filter, time series analysis.
I. INTRODUCTION
Real-time monitoring systems with the use of distributed fiber-optic vibration sensor systems based on the phase-sensitive optical time-domain reflectometry technique [1-6] have a fascinating prospective for applications. Examples include a control for access on protected areas (e.g., securing national borders), oil and gas pipelines, communications lines, and structural health monitoring [79].
The core of such system is phase-sensitive optical time-domain reflectometry technique, which has sufficiently high sensitivity and spatial resolution [10-18]. A main feature of this type of reflectometry is a sufficiently large coherence length of the employed optical pulse. Signals reflected from centers of the Rayleigh backscattering exhibit the coherent summation of their complex wave amplitudes.
On the one hand, monitoring systems based on the phase-sensitive optical time-domain reflectometry are sensitive
Manuscript received July 30, 2015.
Aleksey Fedorov is a student at Bauman Moscow State Technical University (e-mail: akfedorov@student.bmstu.ru).
Maxim Anufriev is an engineer at Bauman Moscow State Technical University (makc.anufriev@gmail.com).
Andrey Zhirnov is a PhD student in Scientific-Educational Center “Photonics and Infrared Technology” at Bauman Moscow State Technical University (aaamizhirnov@mail.ru).
Evgeniy Nesterov is a scientist in Scientific-Educational Center “Photonics and Infrared Technology” at Bauman Moscow State Technical University (evgeny.t.nesterov@gmail.com).
Dmitry Namiot is a senior scientist at Lomonosov Moscow State University (e-mail: dnamiot@gmail.com).
Alexey Pnev is a senior scientist and head of the laboratory in Scientific-Educational Center “Photonics and Infrared Technology” at Bauman Moscow State Technical University (apniov@gmail.com).
Valery Karasik is a head of the Scientific-Educational Center “Photonics and Infrared Technology” at Bauman Moscow State Technical University (karassik@bmstu.ru).
enough to register sufficiently small fluctuations [6]. On the other hand, an algorithm, which allows revealing the nature of fluctuations, should supplement such systems. Indeed, the crucial challenge here is to reveal: are fluctuations caused by natural changes of background or by any kind of activities? In other words, a non-trivial problem of de-noising comes to the fore.
To make a decision about the nature of fluctuations, one can continuously analyze signals from the system in time or frequency domains. Due to a sufficiently complex structure of signals from the monitoring system this problem is rather challenging. Because of very promising application of such systems, this problem has been extensively studied during last decade (see [18-20] and reference therein). However, at this moment there is no universal (device-independent and external-condition-independent) solution for recognition of events for vibration sensor systems based on the phase-sensitive optical time-domain reflectometry technique.
Mathematically speaking, signals from distributed fiberoptic vibration sensor systems are time series [21]. Then for their analysis various types of time series analysis can be applied. In this paper, we present experimental results on an application of the Kolmogorov-Zurbenko filtering [22-26] for signal de-noising in a fiber-optic distributed vibration sensor system based on the phase-sensitive optical time-domain reflectometry. In the considered case, the main goal of the system is a control for access on the protected area. Consequently, we are confronted with the problem of events recognition (e.g., human passage, human group passage, or car travel) in a seismic background.
The paper is organized as follows. In Section II, we describe our setup for collecting experimental data and parameters of the fiber-optic distributed vibration sensor system. In Section III, we describe basic properties of signals from the system. In Section VI, we describe the basic de-noising procedure based on the Kolmogorov-Zurbenko filter. We also show results on the application of the Kolmogorov-Zurbenko filter based de-noising procedure to measured signals. In Section V, we give our conclusion.
II. Setup for collecting experimental data
The fiber-optic distributed vibration sensor system for collecting of data is located on the Bauman Moscow State Technical University polygon in Moscow Region. The setup is presented on Fig. 1: 1 is the primary light source (laser), 2 is the acoustic-optic modulator, 3 is the circulator, 4 is the
16
International Journal of Open Information Technologies ISSN: 2307-8162 vol. 3, no. 9, 2015
fiber-optic sensor, 5 is the preamplifier, 6 is the optical filter, 7 is the detector, 8 is the ADC converter, 9 is the programmable logic device, 10 is the computer.
Probe signal has the wavelength 1550 nm, probe pulse of duration 200 ns, and the signal from the semiconductor laser of power 300-500 mW is launched into the standard optical fiber. Probe signal has ultra-narrow linewidth, which is less than 1 MHz. In our experiments, length of the optical fiber cable l is approximately 50 km.
Fig.1: Fiber-optic distributed sensor system based on the phase-sensitive optical time-domain reflectometry: setup for collecting of experimental data. The sensing element of the system is the standard single mode fiber (fiber-optic cable).
The idea of its work can be presented as follows. In case of a vibration impact, the intensity of backscattering light changes according to the level of the impact. The circulator is used for launching of the probe signal into the optical fiber and the backscattering signal to a detector. Signals in the system are the sum of all scattered signals during the time of pulse with taking into account their phases.
III. Signals from the system: recognition of events
A typical result of the measurement using the setup is presented in Fig. 2a. This figure presents an overlap between signals measured by our setup (Section II) on the region of the cable with length 0.5 km during one second.
As it was mentioned above, the crucial challenge is to recognize any kind of deliberate activity in these signals. It is clear that this problem can be essentially divided into two related sub-problems. The first part is to register the event in a (seismic) background. It is seen from Fig. 2b that an event can be registered by the system via measurement of the difference between signals in neighboring moments of time. From such a procedure, one can find, e.g., a position of the event in the cable. However, due to natural fluctuation of the background such simple procedure leads to a sufficiently large number of false positives.
E(l) 6(1)
0 l 500 0 l 500
Fig. 2: An overlap of typical signals measured by the system on the region of the cable with length 0.5 km during one second (left) and the maximum deviation (right): a signature of an event.
The second problem is to classify registered event. In this work, we are guided by the following basic classification of events: (i) a single event, which is localized in space and in time; (ii) a single event, which is delocalized in space and localized in time; (iii) a single event, which is localized in
space and delocalized in time; (iv) a single event, which is delocalized in space and delocalized in time.
For our setup, an important part is using of preliminary tests. These tests consist of collecting of a large number experimental data corresponding to a background and typical types of activates. Preliminary tests have been organized on the polygon at similar experimental conditions (importantly, the same climate conditions). From preliminary tests, we use a sufficiently large number of experimentally measured signals of type (i)-(iv) to obtain minimum, maximum, and average characteristics values (characteristic time scales of events and characteristic length scales) as well as parameters of the background.
Fig. 3. Measured signals from the system: (a) Background signals experimentally measured by the system; (b) sequence of space localized short-term events with additive noise form background; (c) sequence of space delocalized short-term events with additive noise form background; (d) space localized long-term event with additive noise form background; (e) space delocalized long-term event with additive noise form background; results of the de-noising procedure (f)-(j) for the signals (a)-(e).
The background noise is the waterfall form, i.e., intensity as a function of time (t) and position in the fiber-optic cable (l), is presented in Fig. 3a. A sequence of space-time localized events is presented in Fig. 3b. Fig. 3c shows a sequence of delocalized in space and localized in time events. In Fig. 3d, one can see a sequence of localized in space and delocalized in time events. Fig. 3e shows a sequence of space-time delocalized events.
IV. Basic algorithm
Signals from a distributed fiber-optic vibration sensor systems are time series. Therefore, for their analysis various types of time series analysis can be employed. one of the possible solutions is to use statistical apparatus of the time series analysis.
17
International Journal of Open Information Technologies ISSN: 2307-8162 vol. 3, no. 9, 2015
For instants, an event can be detected in the background in the difference in the value of signals between time moments t=j and t=j+1. Toward this end, we suggest implementing the spatial low-pass filter of the Kolmogorov-Zurbenko type [22-26]. Similar approaches to a de-noising procedure based on moving average calculation have been recently discussed [27,28].
Fig. 4. 3D visualization of the signal measured from the system after the de-noising procedure based on the basic algorithm (1)-(4) and classified by the system as sequence of space and time localized events, which can correspond to human passage along the cable.
In general, the Kolmogorov-Zurbenko filter is a series of iterations of the moving average filter. The first iteration of Kolmogorov-Zurbenko filter is an application of the moving average filter over a given process.
(ж-i):
KZ[JS'('-')],U., - E A"
sj)x —. m
(1)
Here, X(l,t) is the signal (real-valued time series), m is the time window, and k is the filter order [22]. The iteration process of a simple operation of moving average is very computationally efficient. The Kolmogorov-Zurbenko filter has been employed in investigations of climate fluctuations and seminal studies of the turbulence problem [22-26].
Using the de-noising procedure based on the application of the Kolmogorov-Zurbenko filter (1) to signals, the sensor system can detect an event in the background. Indeed, if the moving average deviation of the signal X(l,t)
Sf4 = |x(/,/)-KZ[X{V)]| (2)
exceeds the critical value, then a potential event is detected. The key question for any implementation of such an approach is about determination of the critical value (2).
We suggest using an adaptive critical value calculation in a rather straightforward way. Let us consider the difference between two signals with a time window j in the following form:
Sj (A/) = \X(t + j, /)-X(t, l)\. (3)
By integrating this difference over the cable length
jSj (l)dl = S*
(4)
one can obtain the critical value, which characterizes the integrated difference between two signals from the system in neighboring countdowns of time. On practice, we obtain critical values on the base of preliminary tests of the system (see Section III). This critical value is continuously measured by the monitoring system during preliminary tests.
Finally, from the de-noising procedure with taking into account the critical value, we obtain the de-noised signal as follows:
We optimize the de-noising procedure with respect to parameters of the Kolmogorov-Zurbenko filter: the size of the time window and order of the filter. This allows obtaining the maximal level of recognized events.
The results of our application of the de-noising procedure (1)-(4) for experimentally measured signals of type (i)-(iv) are presented in Fig. 3. In spite of the simplicity of used de-noising procedure, events can be clearly detected. The results are presented in Fig. 3. Furthermore, patterns of different classes of events can be reliably revealed (See Fig. 4).
V. RESULTS AND CONCLUSION
Controllability of a workflow is the most important requirement for its efficient implementation. Then a problem of a design of real-time monitoring systems for nondestructive testing has a paramount importance for technological proces ses in a broad class of applications in industry. Mainly, the goal of a monitoring system is to measure and analyze specific characteristics of a process, which allows to make a decision about the correctness of its functioning.
In this work, we have investigated the basic algorithmic solution for recognition of events in seismic backgrounds. Performance of the suggested algorithmic solution has been demonstrated on experimentally measured signals.
REFERENCES
[1] Barnoski M.K., Jensen S.M. Fiber waveguides: a novel technique for investigating attenuation characteristics // Applied Optics - 1976 -Vol. 15 - pp. 2112-2115.
[2] Costa B., Sordo B., Menaglia U., Piccari L., Grasso G. Attenuation measurements performed by backscattering technique // Electronics Letters - 1980 - Vol.16 - pp. 352-353.
[3] Barnoski M.K., Rourke M.D., Jensen S.M., Melville R.T. Optical time domain reflectometer // Applied Optics - 1977- Vol. 16 - pp. 2375-2380.
[4] Healey P., Malyon D.J. OTDR in single-mode fibre at 1.5-um using heterodyne detection // Electronics Letters - 1982 - Vol. 18 - pp. 862-863.
[5] Healey P. Statistic of Rayleigh backscatter from a single-mode optical fibre // Electronics Letters - 1985 - Vol. 21- pp. 226-228.
[6] Healey P. Instrumentation principles for optical time domain reflectometry // Journal of Physics E: Scientific Instruments - 1986 -Vol. 19 - pp. 334-341.
[7] Raj B., Jayakumar T., Thavasimuthu M., Practical non-destructive testing. Woodhead Publishing, 2002.
[8] Cau F., Fanni A., Montisci A., Testoni P., Usai M. A signalprocessing tool for non-destructive testing of inaccessible pipes //
18
International Journal of Open Information Technologies ISSN: 2307-8162 vol. 3, no. 9, 2015
Engineering Applications of Artificial Intelligence - 2006 - Vol. 19 -pp. 753-760.
[9] Bao X., Chen L. Recent progress in distributed fiber optic sensors, Sensors - 2012 - Vol. 12 - p. 8601-8639.
[10] Tai K., Hasegawa A., Tomita A. Observation of modulational instability in optical fibers // Physical Review Letters - 1986 - Vol. 56 - pp. 135-138.
[11] Takada K., Himeno A., Yukimatsu K. Phasenoise and shotnoise limited operations of low coherence optical time domain reflectometry // Applied Physics Letters - 1991 - Vol. 59 - P. 24832485.
[12] Juskaitis R., Mamedov A.M., Potapov V.T., Shatalin S.V. Distributed interferometric fiber sensor system // Optics Letters - 1992 - Vol. 17
- pp. 1623-1625.
[13] Kulchin Yu.N., Vitrik O.B., Dyshlyuk A.V., Shalagin A.M., Babin S.A., Vlasov A.A. Application of optical time-domain reflectometry for the interrogation of fiber Bragg sensors // Laser Physics - 2007 -Vol. 17 - pp. 1335-1339.
[14] Koyamada Y., Imahama M., Kubota K., Hogari K. Fiber-optic distributed strain and temperature sensing with very high measure and resolution over long range using coherent OTDR // IEEE Journal of Lightwave Technology - 2009 - vol. 27 - pp. 1142-1146.
[15] Martins H.F., Martin-Lopez S., Corredera P., Salgado P., Frazao O., Gonzalez-Herraez M. Modulation instability-induced fading in phase-sensitive optical time-domain reflectometry // Optics Letters - 2013 -Vol. 38 - pp. 872-874.
[16] Alekseev A.E., Vdovenko V.S., Gorshkov B.G., Potapov V.T., Simikin D.E. A phase-sensitive optical time-domain reflectometer with dual-pulse phase modulated probe signal // Laser Physics - 2014
- Vol. 24 - P. 115106.
[17] Duan N., Peng F., Rao Y.-J., Du J., Lin Y. Field test for real-time position and speed monitoring of trains using phase-sensitive optical time domain reflectometry (Ф-OTDR) // SPIE Proceedings - 2014 -Vol. 9157 - 91577A.
[18] Shi Y., Feng H., An Y., Feng X., Zeng Z. Research on wavelet analysis for pipeline pre-warning system based on phase-sensitive optical time domain reflectometry // IEEE/ASME International Conference on Advanced Intelligent Mechatronics (AIM) Besancon, France.
[19] Li Q., Zhang C., Li L., Zhong X. Localization mechanisms and location methods of the disturbance sensor based on phase-sensitive OTDR // Optik - 2014 - Vol. 125 - pp. 2099-2103.
[20] Lyons W. B., Lewis E. Neural networks and pattern recognition techniques applied to optical fibre sensors // Transactions of the Institute of Measurement and Control - 2000 - Vol. 22 - pp. 385-404.
[21] Kuklin S.N., Namiot D.E. Ob odnom podhode k obrabotke vremennyh rjadov // International Journal of Open Information Technologies - 2015 - Vol. 3 - pp.18-22.
[22] Zurbenko I. The spectral analysis of time series, North-Holland Series in Statistics and Probability.
[23] Zurbenko I., Porter P., Rao S., Ku J., Gui R., Eskridge R. Detecting discontinuities in time series of upper air data: Development and demonstration of an adaptive filter technique // Journal of Climate -1996 - Vol. 9 - pp. 3548-3560.
[24] Yang W., Zurbenko I. Kolmogorov-Zurbenko filters // Wiley Interdisciplinary Reviews: Computational Statistics - 2010 - Vol. 2 -pp. 340-351 (2010).
[25] Zurbenko I., Luo M. Restoration of time-spatial scales in global temperature data // American Journal of Climate Change - 2012 -Vol. 1 - pp. 154-163.
[26] Zurbenko I.G., Potrzeba A.L. Periods of excess energy in extreme weather events // Journal of Climatology 2013 - Vol. 2013 - P. 410898.
[27] Yang J. Toward physical activity diary: motion recognition using simple acceleration features with mobile phones // Proceeding of the 1st Int. Workshop on IMCE at ACM Multimedia 2009, ACM Press, 2009, pp. 1-10.
[28] Khan A.M., Siddiqi M.H., Lee S.-W. Exploratory data analysis of acceleration signals to select light-weight and accurate features for real-time activity recognition on smartphones // Sensors - 2013 - Vol. 13 - pp. 13099-13122.
19
