GB2539254A

GB2539254A - Method and apparatus for monitoring pipeline

Info

Publication number: GB2539254A
Application number: GB1510267.6A
Authority: GB
Inventors: Hanisch Christoph
Original assignee: PIMON GmbH
Current assignee: PIMON GmbH
Priority date: 2015-06-12
Filing date: 2015-06-12
Publication date: 2016-12-14
Also published as: GB201510267D0; EP3332272A1; WO2016198683A1

Abstract

An optical fibre sensor system is used to determine the type and location of disturbances along a pipeline, such as digging, drilling, tapping, sabotage etc. The system includes an optical fibre 20 buried along the pipeline 35, a radiation source 30 (e.g. infra-red source), a detector 40 and a signal processor 50. The processor analyses values corresponding to the condition of the optical fibre to create a pattern corresponding to a disturbance, and compares the pattern to recorded patterns of known disturbances.

Description

Title: Method and Apparatus for monitoring pipeline

BACKGROUND OF THE INVENTION

Field of the Invention

[0001] The invention relates to an optical fibre sensor system and a method for determining a location and type of a disturbance along a pipeline.

Introduction to the Invention

[0002] Oil, water, gas and other product pipelines form a critical network in every part of the world and are an easy target for intruders. The pipelines are also susceptible to earthquakes, tsunamis and to other geohazard incidents. Monitoring the pipelines in order to keep the pipelines safe from damage is a major challenge. The long distances, often through remote and hostile territory, make the costs of most conventional monitoring systems prohibitive. If an oil or gas pipeline is damaged it can have devastating impacts on human life and health due to explosions, fire and contamination; environment due to poisoning of flora and fauna, as well as the associated financial losses and damage to both image and reputation.

[0003] The pipelines are susceptible to various types of third-party interference, such as deliberate acts (illegal tapping, sabotage) or unintended disruption (construction work, farming). This third-party interference can cause huge financial and environmental damage and loss of reputation to the pipeline operators. A reliable, real-time pipeline monitoring system is therefore required in order to protect nature, human health and economic interests.

[0004] According to a survey in 2009, about 36% of all worldwide pipeline leaks were caused by third party interference (TPI), such as illegal tapping, sabotage or construction work. Examples of the third party interference includes: 04 October 2001 Fairbanks Alaska (USA), 990t oil due to sabotage; 2000 Tschernigow, Ukraine, loss of 500,000 Liter diesel due to illegal tapping; 10 June 1999 Bellingham, Whatcorn Creek, Washington (USA), 880t petrol due to construction work, financial damage USD 45 million.

[0005] Illegal tapping is a major problem in emerging countries like India, China and South America. In 2011 for example, Petroleos Mexicanos (PEMEX) counted 1,324 cases of illegal tapping in Mexico. Every day Pemex estimates 40,000 liters of oil and gas, which sums up to an annual damage of more than 1 billion US$ of lose. Such damage could be avoided or substantially reduced, if third party interference was detected before or close to the occurrence of the interference happens.

[0006] Many times the pipeline monitoring is done by walking, driving and flyingalong the pipeline. The annual costs for walking and driving the line vary from 100 to 350 per kilometre. The additional costs for flying-the-line amount to 4.50 per kilometre. Usually two inspections are performed per month, which adds up to 108 per kilometre ( 4.50/km x 2 flights/month x 12 months). The annual costs are between 208 and 458 per kilome-tre.

Description of the Related Art

[0007] A number of systems for the sensing of an acoustic disturbance are known. The acoustic disturbance can be representative of damage to the pipelines through third party interference or geohazards, as described in the introduction. These systems involve the use of an optical fibre laid alongside the pipeline, which acts as a sensor and detects changes in the pattern of back-scattered radiation in order to sense an acoustic disturbance. For exam- ple, Pimon GmbH, Munich, Germany, sells an apparatus PMS2500-vibrO that utilizes dis-tributed fibre optical sensing technology to detect the acoustic disturbance. The PMS2500-vibrO system combines an optical time domain reflectometer (OTDR) with an analysis and pattern recognition software and offers a customized interface with geographic information system (GIS) mapping.

[0008] International Patent Applications No WO 2011/05813, WO 2011/ 015812 and WO 2011/059501 (QinetiQ Ltd) all teach various aspects of using a distributed fibre optic sensing system for establishing events of interest. Similarly UK Patent Application No GB 2 491 658 also teaches a method and system for locating an acoustic disturbance. These patent applications all have in common that the systems analyse the back-scattered radiation from the optical fibre to establish the event of interest. Such systems are useful in determining an event of interest from the acoustic disturbance, but the systems are known to produce "false positives" in which events are identified that are of no interest and fail to identify some events of interest, in particular when such events have not been seen before.

[0009] One solution to the issue of incorrect identification of events would be to use an artificial neural networks (ANN) to train the system to recognise the event. The ANNs are computational models and are inspired by animal central nervous systems, in particular the brain, that are capable of machine learning and pattern recognition. The ANNs are usually presented as a system of nodes or "neurons" connected by "synapses" that can compute values from inputs, by feeding information from the inputs through the ANN. The synapses are the mechanism by which one of the neurons passes a signal to another one of the neurons.

[0010] One example of the ANN is for the recognition of handwriting. A set of input neurons may be activated by pixels in a camera of an input image representing a letter or a digit. The activations of these input neurons are then passed on, weighted and transformed by some function determined by a designer of the ANN to other neurons, etc. until finally an output neuron is activated that determines which character (letter or digit) was imaged. ANNs have been used to solve a wide variety of tasks that are hard to solve using ordinary rule-based programming, including computer vision and speech recognition.

[0011] There is no single formal definition of an ANN. Commonly a class of statistical models will be termed "neural" if the class consists of sets of adaptive weights (numerical parameters that are tuned by a learning algorithm) and are capable of approximating nonlinear functions of the inputs of the statistical models. The adaptive weights can be thought of as the strength of the connections (synapses) between the neurons.

[0012] The ANNs have to be trained in order to produce understandable results. There are three major learning paradigms: supervised learning, unsupervised learning and reinforcement learning.

[0013] In a supervised learning, the learning paradigms all have in common that a set of pre-analyzed data, for example a waveform, is analyzed by the ANN and the weights of the connections (synapses) between the neurons in the ANN are adapted such that the output of the ANN is correlated with a known event. There is a cost involved in this training. An improvement in the efficiency of the results of the ANN can be obtained by using a greater number of data items representing the known event in a training set. The greater number of data items requires, however, an increase in computational power and time for the analysis in order to get the correct results. There is therefore a trade-off that needs to be established between the time taken to train the ANN and the accuracy of the results.

[0014] Recent developments in the ANNs involve so-called 'deep learning'. Deep-learning is a set of algorithms that attempt to use layered models of inputs. Jeffrey Heaton, University of Toronto, has discussed deep learning in a review article entitled 'Learning Multiple Layers of Representation' published in Trends in Cognitive Sciences, vol. 11, No. 10, pag-es 428 to 434, 2007. This publication describes multi-layer neural networks that contain top-down connections and training of the multilayer neural networks one layer at a time to generate sensory data, rather than merely classifying the data.

[0015] Neuron activity in the prior art ANNs is computed for a series of discrete time steps and not by using a continuous parameter. The activity level of the neuron is usually defined by a so-called "activity value", which is set to be either 0 or 1, and which describes an action potential' at a time step t. The connections between the neurons, i.e. the synapses, are weighted with a weighting coefficient, which is usually chosen have a value in the in- terval [-LO, + 1.0]. Negative values of the weighting coefficient represent "inhibitory syn-apses" and positive values of the weighting coefficient indicate "excitatory values". The computation of the activity value in the ANNs uses a simple linear summation model in which weighted ones of some or all of the active inputs received on the synapses at a neuron are compared with a (fixed) threshold value of the neuron. If the summation results in a value that is greater than the threshold value, the following neuron is activated.

[0016] One example of a learning system is described in international patent application No. WO 1998 027 511 (Geiger), which teaches a method of detecting image characteris-tics, irrespective of size or position. The method involves using several signal-generating devices, whose outputs represent image information in the form of characteristics evaluated using non-linear combination functions.

[0017] International patent application No. WO 2003 017252 relates to a method for rec- ognizing a phonetic sound sequence or character sequence. The phonetic sound sequence or character sequence is initially fed to the neural network and a sequence of characteristics is formed from the phonetic sequence or the character sequence by taking into consideration stored phonetic and/or lexical information, which is based on a character string se-quence. The device recognizes the phonetic and the character sequences by using a large knowledge store having been previously programmed.

[0018] An article by Hans Geiger and Thomas Waschulzak entitled Theorie and Anwendung strukturierte konnektioni sti sche Systeme', published in Informatik-Fachrei chte, Springer-Verlag, 1990, pages 143 -152 also describes an implementation of a neural net-work. The neurons in the ANN of this article have activity values between zero and 255. The activity values of each one of the neurons changes with time such that, even if the inputs to the neuron remain unchanged. The output activity value of the neuron would change over time. This article teaches the concept that the activity value of any one of the nodes is dependent at least partly on the results of earlier activities. The article also in-cludes brief details of the ways in which system may be developed.

[0019] Multiple tests have shown that the systems as described in the prior art do not deliver reliable results. A system and a method for determining a location and type of a dis- turbance along a pipeline delivering reliable results is needed therefore, which is compara-tively easy to install, run, and maintain.

SUMMARY OF THE INVENTION

0020] The principal of the method and apparatus of determining locations and types of disturbances along a pipeline, as described in this disclosure, is based upon analysing back-scattered radiation in an optical fibre using a signal processor.

0021] More specifically, a fibre sensor system for determining a location and type of a disturbance along a pipeline is provided. The fibre sensor system comprises an optical fibre, a radiation source for launching radiation into the optical fibre, a detector for detecting back-scattered radiation back-scattered from the optical fibre, and a signal processor connected to the detector for generating a plurality of values from the back-scattered radiation, corresponding to the condition of the optical fibre at a certain time at a certain location of the optical fibre. The signal processor comprises a recorder to record the values corresponding to the condition of the optical fibre at a specific repetition rate, an analyser to analyse the change of those values per repetition over time and create a value over time pattern corresponding to the specific disturbance.

Q022] The radiation source as mentioned above is in one aspect designed to produce a series of optical pulses.

0023] In another preferred embodiment, the fibre sensor system comprises a memory for storing a plurality of patterns representative of known disturbances, and a means to compare the value over time patterns corresponding to the specific disturbance as created with the recorded patterns representative of known disturbances.

[0024] In another aspect of the disclosure, the radiation source produces infrared light in a series of optical pulses.

[0025] In a further aspect of the disclosure, the specific repetition rate at which the values corresponding to the condition of the optical fibre are recorded, is 250 Hz.

[0026] Also disclosed is a method for determining a location and type of a disturbance along a pipeline, the method comprising the steps of launching radiation into an optical fibre buried along a pipeline, detecting back-scattered radiation back-scattered from the optical fibre, generating a plu-rality of values from the back-scattered radiation, corresponding to the condition of the optical fibre at a certain time at a certain location of the optical fibre; recording the values corresponding to the condition of the optical fibre at a specific repeti-tion rate, analysing the development of those values per repetition over time; and creating a value over time pattern corresponding to the specific disturbance.

[0027] In another aspect of the disclosure, the method for determining a location and type of a disturbance further comprises the step of comparing the value over time pattern with recorded patters of known disturbance types.

[0028] In a further aspect of the disclosure, the values corresponding to the condition of the optical fibre are recorded at a preferred specific repetition rate of 250 Hz. Currently, the preferred spatial sampling rate is 1 meter (this corresponds to 100 MHz ADC sampling freuctenzy). Another temporal and spacial sampling rates may be used, as the sampling frequencies may be varied within a certain range.

[0029] The optical fibre sensor system taught in this disclosure is able to provide: alert in case of emergency with an accurate GIS location, identify threats in real time in order to reduce the risk from leakages caused by digging, drilling, tapping, sabotage, earthquake etc., monitor around 50 km of pipeline with one measuring unit without any additional power supply and with scalability to cover deployments of any length by adding additional units. The system reduces false or irrelevant alarms and is capable of differenti-ating multiple events down to 5-15 meter resolution; as well as track potential intruders, moving in different velocities (vehicles or on foot) along the pipeline. The system can detect direction and speed of the possible intruders.

DESCRIPTION OF THE FIGURES

[0030] Fig. 1 shows an example of an optical fibre sensor used in the system of the disclosure [0031] Fig. 2 to 6 illustrate graphs showing various outputs generated by processing of a back-scattered signal in a processor.

[0032] Fig. 7 is a flow chart, describing the steps of the method of this disclosure.

[0033] Fig. 8 is a graph showing peaks in the frequency power spectra typical for diesel engines.

DETAILED DESCRIPTION OF THE INVENTION

[0034] The invention is described on the basis of the drawings. It will be understood that the embodiments and aspects of the invention described herein are only examples and do not limit the protective scope of the claims in any way. The invention is defined by the claims and their equivalents. It will be understood that features of one aspect or embodi-ment of the invention can be combined with a feature of a different aspect or aspects and/or embodiments of the invention.

[0035] The optical fibre sensor system 10 of this disclosure is based on distributed fibre optical sensing and integrates an optical time domain reflectometer (OTDR) with detection software, as well as a customized interface with geographic information system (GIS) mapping. The fibre optic cables are often already buried near to or attached to the pipeline for telecommunication purposes. The system turns a standard single mode telecommunication fibre into a listening device for events of interest.

[0036] The optical time-domain reflectometer (OTDR) is an optoelectronic instrument used to characterize optical fibres. The reflectometer injects series of optical pulses into the optical fibre that is being tested and detects at the same fibre end the returning ("back-scattered") light that has been scattered or reflected back from points along the fibre. The OTDR produces a so-called reflectogram from this measurement and performs an efficient, precise and wide analysis of the fibre characteristics.

[0037] A laser diode launches radiation in the form of a plurality of optical pulses into one end of the optical fibre and a photodiode measures the returning light. The photo diode is connected to the OTDR. The back-scattered signal detected by the photodiode provides relevant information about the events of interest along the optical fibre. The back-scattered signal is influenced by both the attenuation and reflections, which the injected laser pulse experiences on its way through the optical fibre. If a certain area within the optical fibre has a higher degree of attenuation or reflection, for instance caused by a bending or a fibre connector, this higher degree of attenuation or reflection is detected by differences in the back-scattered signals. The velocity of the optical pulse is known and the exact location of the event of interest can therefore be determined on the basis of the time difference between the injection of the optical pulse and the return of the back-scattered signal.

[0038] Suppose now that a disturbance occurs near the optical fibre. The pressure on the optical fibre will cause at least a small degree of bending of the optical fibre and will result in a change of the returned back-scattered signal. This change can be detected in the OTDR. Suppose that the disturbance is a moving object, then the slight bending of the op- tical fibre will change location, enabling the speed and direction of movement of the dis-turbance to be established.

[0039] The optical fibres make a good sensor, as the optical fibres are able to measure vibrations over long distances. The reflections of the signals in the optical fibres change on vibrations (e.g. caused by vehicles, footsteps, digging, drilling), on temperature alterations (e.g. caused by escaping pressurized gas), or when optical fibres get strained, bent, kinked, or cut off

Example

[0040] Fig. 1 shows an example of the optical fibre sensor system 10 using the teachings of this disclosure. An optical fibre 20 used in this sensor system 10 and can be a specially laid optical fibre 20, which has been placed in a region of interest. The optical fibre 20 can also be a standard telecommunications optical fibre, which generally carries data. The optical fibre 20 can be laid near or directly adjacent to a pipeline 25, as outlined in the introduction

to the description.

[0041] The optical fibre 20 is connected to a radiation source 30, which in one non-limiting aspect of the invention is a semiconductor laser producing radiation 35 comprising a plurality of optical pulses at 1.55pm and launches the radiation 35 into the optical fibre at repetition rates up to 4 kHz in the present embodiment.

[0042] The returning back-scattered signal 37 from the optical fibre 20 can be observed at the same end of the optical fibre 20 in form of a reflectogram by a detector 40 in which the time delay of signal incidence is shown as distance on the optical fibre 20. A highly coherent laser as the radiation source 30 is used in order to increase the sensitivity. Due to this high level of coherence, an interference pattern can be observed in the back-scattered sig- nal 37 (fingerprint). Vibrations cause bending in the optical fibre 20 and result in temporary changes in this interference pattern. The interference pattern can be used to detect a disturbance 45, for example a third-party interference caused by manual or machine digging. Every kind of incident causes certain multiple vibrations and changes the back-scattered signals 37 leading to differences in the interference pattern. The back-scattered signals 37 are received and the signal patterns are digitalized. The signal in the detector 40 passes through to a signal processor 50. The signal processor 50 is able to classify the disturbance 45 and send an alarm to the system operator, if required.

[0043] The signal processor 50 is designed to separate all "regular" sounds from possible disturbances 45 in the back-scattered signal as well as any inoffensive or other irrelevant incidents.

[0044] The processor 50 further includes an input device (not shown) that is used to input information items relating to the reflectogram. The information items may include a name or a label generally attached to the reflectogram and/or to one or more features in the reflectogram. The input device is connected to a signal processor (not shown) having a memory. The signal processor compares the characteristics relating to a particular feature in a stored reflectogram with the inputted information and can associate the particular stored reflectogram (or a portion thereof) with the inputted information. This association is memorized so that if an unknown pattern or feature in a reflectogram is detected by the detector 40, the processor 50 can determine that this unknown pattern is in fact a known pattern and output the associated item of information to the user at an output.

[0045] The pattern recognition system 10 can be trained to recognize a large number of patterns in the reflectogram using an unsupervised leaning process. This unsupervised learning process is carried out using a set of prestored patterns representing features of interest from existing reflectograms and running the prestored patterns through the processor 50.

[0046] The system and method of the current disclosure can be used to determine and clas-sift' unknown disturbances 45, as shown in Fig. 2 to 6.

[0047] In this example of the system and method the radiation 35 is launched into the optical fibre 20. The detector 40 receives the back-scattered signal 37 and passes the back-scattered signal 37 to the processor 50.

[0048] A label can be associated with the one or more features in the back-scattered signal.

[0049] The procedure is repeated for a different feature. This different feature creates a different structure within the plurality of layers 180. The learning procedure can then pro-teed using different ones of the features.

[0050] Processing of the back-scattered signal 37 in the processor 50 is described with respect to Figs. 2 to 6.

[0051] The measuring data to be evaluated comprises a value matrix spanned by sampled fibre range and the temporal sample steps within a certain observation time. It must be distinguished between detection of location and detection of time of the fibre (20). For each laser pulse injected into the fibre (20), the back-scattered signal is detected. Via the time elapsed between injecting the pulse and detection of a certain part of the back scattered signal, the location is determined by means of the run time of the light (OTDR-Method), from which the scattered signal emanates. An analogue-digital-converter scans the back scattered signal at 100 NIEL, so that for each meter of fibre a back scattered value can be allocated. The back scattered values for the entire fibre (e.g. 20000 values for 20 km) per injected laser pulse is called a fingerprint. A fingerprint is a line in the matrix to be evalu-ated. The temporal development of the fibre is obtained by repeated scanning of the fibre, i.e. by repeatedly injecting laser pulses with a repetition rate of e.g. 250 Hz. The temporal development of a fibre location is derived from evaluating subsequent fingerprints at this location. The values at such a location form the columns of the matrix and are called timeslice.

[0052] In Fig. 2 the evaluation is carried out for the example of seven pickaxe beatings 201 to 207 carried out at a sample installation of the device according to the disclosure. Due to the specific physical layout of the sample installation, the same action is detected in three positions, namely at 175 m, 230 m and 530 m from the measuring distance [0053] The classification of disturbances on the optical fibre 20, which is due to activities around the optical fibre 20, is carried out by mainly investigating the temporal develop-ment or by recognising patterns within a single one of the time slices, wherein intermediate results of adjoining timeslices are also considered (stribe condition). For the detection of manual excavation activities the following steps are carried out: [0054] Step I: Use of a Laplace operator. This operator generates positive or negative peaks at those locations, at which the input signal shows a substantial change. Small changing values are those values close to zero and steady input values are those values that are precisely zero. This Laplace operator improves the signal-noise relation of the measured values. The one dimensional Laplace operator is a vector of coefficients which are each to be multiplied with one out of an equal number of subsequent input values. Similar to an scalar product the sum of all products ensues the output value at the position the center coefficient has been applied to. The number of coefficients is odd. The coefficients comply three requirements: The middle one has a positive maximum value. The coefficients are symmetric with regard to the central one and behave like a Gaussian or a Lorentz distribution. The sum of all coefficients is zero.

[0055] In Fig. 3 the impact of the Laplace operator is shown. The Graph in Fig. 3c shows the Laplace data, which arise from the signal of a pickaxe. The amount is formed by the result of the Laplace operator, so that one works with positive peaks only, see the Graph in Fig. 3d.

[0056] Step 2: The Laplace output often contains solitary peaks which are not caused by a digging stroke. Stroke events instead show a whole cluster of peaks that enables one to locate the time of the hit by locating the highest density of peaks in such a cluster. In order to neglect single peaks and to identify the most likely hit time the amount of Laplace out-put within a time slice is averaged as described below. By doing so single peaks are dwarfed and one can assign a hit time to each peak cluster. The latter is a local maxima of the averaged timeslice. In order to achieve the above-quoted averaging of the absolute Laplace data, a sliding average is formed for every single point within the timeslice, e.g. a value is summed using the five direct preceding values and five direct succeeding values.

After a division by the number of all those values, here eleven values, the result of the slid-ing average substitutes the previous value. By doing so, the graphs are smoothed, and artifacts, such as isolated peaks, can be supressed, as can be seen in Fig. 3b.

[0057] Step 3: Definition of a threshold as a function of the noise background, by multiply-ing the noise average with an eligible factor, and sorting out the over threshold values.

[0058] Step 4: Use of a "step condition": Values which are over the defined threshold values can also be generated by a person walking around over the optical fibre 20 or by a car near the optical fibre. To be considered to be a blow of a pickaxe or a spade, over threshold values must fulfil two additional prerequisites, the step condition and the stripe condition.

[0059] The step condition: Starting point or symmetric point is the temporal position of values which are over the threshold values found in the 3rd step. In a defined time distance from this temporal position two so-called value areas that are equal in size are looked at which are located in time before and after the main peak at the temporal position. For ex-ample, suppose that the peak is located at a temporal peak at time t=200, then the two value areas chosen are a right value area stretching from t=220 to t=320 and a left value area stretching from the left from t=80 to t=180. For these value areas an average is formed in form of the values in the raw signal over the value area (measuring values before step 1) in each case. The step condition requires that the amount of the difference of both averages may not surpass a certain value. If it does, the peak is at or close to a clear step in the basic level of the raw signal Steps can be seen in the signal of a pedestrian. They are, however, extremely rare for a pickaxe signal, which is shown the Graph in Fig. 3a.

[0060] Step 5: Use of a "stripe condition": the remaining values are now compared with other, neighbouring values, in the value matrix, meaning that a value will be compared to other values in this time slice, but also with values from the neighbouring ones of the time slices. For pickaxe signals, it is typical that the blow of the pickaxe in the ground generates almost always several values in a relatively narrow stripe (see Fig. 2). This stripe is given by an eligible delta t and an eligible delta x. A value fulfils the condition if there are sufficient other values next to it in the stripe area as defined by delta t and delta x, so that a defined number of values is reached in a "cluster", e.g. 5. Values 201 to 207 shown in Fig. 2 have gone through all those steps. These values can be supposed to have been created in the raw signal due to an excavation or a pickaxe activity.

[0061] Signals of pedestrian can be classified into those signals created by walking across the sensor cable and those signals originating from walking along the sensor cable. If someone walks across the sensor cable, at right angles to the sensor cable, patterns of chronologically separated steps can be recognised in the affected timeslice (see Fig. 4 401). Even if the a pedestrian walks on the same spot, the measured signal looks signal 401 in Fig. 4. A clear deviation from the raw signal, which usually shows a length of about 0.5 seconds for single steps is, typical. This deviation in the raw signal can be formulated as an additional condition in the algorithm described above for the detection of a blow by a pick-axe. If the step condition described above was not fulfilled, and the event was rejected as a blow, a pedestrian step event can be ascertained by examining this additional condition of the deviation from the raw signal. This is valid only for single pedestrian steps. If the pedestrian runs along the sensor cable, the steps are not so clearly separated from each other.

Then the movement of the pedestrian can be detected by use of the following algorithm.

[0062] To show movements of the pedestrians or the vehicles over the optical fibre 20, a spline algorithm is used. In contrast to the narrow peaks in the raw data relating to the excavation activities, steps in the raw data or displacements are typical for other movements in the signal. The spline algorithm is a mathematical function, which interpolates the de-velopment of the measuring data. Narrow peaks in the signal as well as noise are ignored. The algorithm comprises the following steps (see Fig. 4 and 5): [0063] Step 1: Filtering out or suppression of narrow peaks. Should a value differ clearly in the raw signal from the average of the values in close proximity, the value is substituted with the average. This is illustrated in the graph 402 and 502 in Fig. 4 for two single steps of a pedestrian and Fig. 5 for a pedestrian walking along the optical fibre, in which a sub-stitution graph is constructed from the raw data in the upper graph 401 and 501.

[0064] Step 2: Smoothing of the data by use of a gliding average as in graph 403 and 503).

[0065] Step 3: Calculation of a spline, which returns only the substantial changes of the signal, graphs 404 and 504.

[0066] Step 4: The amount of the derivative (slope) of the spline (graph 404 in Fig. 4 and 504 in Fig. 5) is always high when there are clear steps in the raw signal [0067] Step 5: Display of the values, which are larger than a certain threshold value 410 and 510, as a signal of the pedestrian, graphs 405 and 505.

[0068] Furthermore, an alternative algorithm was developed, which is suitable to detect pedestrian movement. This algorithm has a stronger response to digging activities than the Spline-algorithm and may only be used effectively, if by means of the above algorithm (step and stripe condition) a digging activity can be excluded. Similar to the step condition described within each timeslice for each point in time the passed raw signal is compared with the following raw signal. These two areas are distributed symmetrically around the position and are separated by a gap. E.G.: Position = 200; Values between 80 and 180 form an area 1; values between 220 and 320 form an area 2. The difference value of the averages of the two areas is a function with particularly large amplitudes for pedestrian activities. To supress noise values below a threshold are set to zero and only values above the threshold are displayed. A maximum of this function appears for steps in the basic signal, in particular when the step condition above excludes a hit event.

[0069] In Fig. 6 the signal of the pedestrian is shown that walks along the sensor cable in one direction and then turns back.

[0070] Any stationary vehicles with diesel engines (tractors, JCBs, trucks) can be detected very simply by means of a Fourier transformation. To do so, the Fourier transformation is calculated in a calculation window (segment of a time slice, e.g., 512 values) and the power spectrum derived therefrom. The basic oscillation and harmonics of the engines are visi- ble as peaks around the power spectrum. For every monitored fibre position the Fourier- transform-algorithm is applied. This means the previous 512 values of each timeslice are taken as the input vector of Fast Fourier Transform. This gives an output vector of 512 complex values. Deriving the square of the absolute output values gives the frequency-power spectra of the harmonic components in the input between zero and the Nyquist-frequency. This is half the frequency that has been used when recording the input values, that means for a 250 Hz temporal sampling rate the Nyquist-frequency is 125 Hz. The high peak at zero Hz is caused by the constant input signal offset and is to deduct or to neglect. The harmonic vibration of a diesel engines generates clear peaks in the frequency power spectra as displayed in the plot shown in Fig. 8. The fundamental frequency and its higher- order harmonics are to observe. For setting an alarm state a suitable threshold can be de-fined that the peak in the power spectra has to surpass in order to switch to an alarm state.

[0071] Moving vehicles can be distinguished from the pedestrians as a function of their weight and their speed by the size of a segment of the optical fibre in which the moving vehicles cause a signal. While the pedestrian crossing the signal generated in the sensor cable is seen in a response area of 3 to 5 metres along the sensor cable, the signal generated in the sensor cable relating to a crossing car is between 10 to 12 metres of the response area. The signal relating to a truck has a response area which is bigger than 20 metres.

[0072] For the moving vehicles, about half of the affected time slices show a typical wave-form, like half an oscillation of a sinus wave, of the measuring signal. After subjecting the timeslice to a gliding average (averaging over 50 to 100 values) the wave form is preserved while the background noise and isolated high peaks disappear. Such a wave form is not included in the pedestrian's signal and, hence, points with all probability to a moving car.

The FWHM of the wave form points to the time needed to cross the cable.

[0073] Should the system 10 be unable to identify the feature, then the system 10 can give an appropriate warning and human intervention can be initiated in order to review the disturbance and classify the unknown feature or to resolve in the other conflicts.

[0074] Fig 7. is a flow chart, describing the steps of the method of this disclosure. The method for determining a location and type of a disturbance along a pipeline, comprises the following steps: [0075] The first step 701 comprises the injection of radiation into an optical fibre burried along a pipeline, followed by the detection of the back-scattered radiation in the next step 702, which was back-scattered from the optical fibre and which is generating a plurality of values, corresponding to the condition of the optical fibre at a certain time at a certain location of the optical fibre; the third step 703 comprises the recording of the generated values corresponding to the condition of the optical fibre at a specific repetition rate, analysing the development of those values per repetition over time; and creating a value over time pat-tern corresponding to the specific disturbance; the location and the type of a disturbance is determined in the fourth step 704 by comparing the value over time pattern with recorded patterns of known disturbance types.

Claims

Claims
2.
3.
4
5.
6.

An optical fibre sensor system (10) for determining a location and type of a disturbance (85) along a pipeline comprising: an optical fibre (20) to be buried along the pipeline; a radiation source (30) for launching radiation (35) into the optical fibre (20); a detector (40) for detecting back-scattered radiation (37) back-scattered from the optical fibre (20); a signal processor (50) connected to the detector (40) for generating a plurality of values (55) from the back-scattered radiation (37), corresponding to the condition of the optical fibre at a certain time at a certain location of the optical fibre; wherein the signal processor (50) comprises means a recorder to record the values corresponding to the condition of the optical fibre (20) at a specific spatial and temporal sampling rate, an analyser to analyse the development of the values per repetition over time and create a value-over-time pattern corresponding to the specific disturbance, and a comparator to compare this value over time pattern with recorded patterns of known disturbance types.

The optical fibre sensor system (10) of claim 1, wherein the radiation source (30) launches a series of optical pulses (35p).

The optical fibre sensor system (10) of claim 2, wherein the radiation source (30) launches infrared pulses The optical fibre sensor system (10) of any of the above claims, further comprising a memory (190) for storing a plurality of patterns (196) representative of disturbances (45.

The optical fibre sensor system (10) of any of the above claims, wherein the specifis repetition rate at which the values corresponding to the condition of the optical fibre are recorded is 250 Hz.

A method for determining a location and type of a disturbance (85) along a pipe-line, the method comprising: launching (205) radiation (35) into an optical fibre (20) buried along a pipeline; detecting (210) back-scattered radiation (37) back-scattered from the optical fibre (20); generating a plurality of values (55) from the back-scattered radiation (37), corre-sponding to the condition of the optical fibre (20) at a certain time at a certain location of the optical fibre (20); recording the values (55) corresponding to the condition of the optical fibre (20) at a specific repetition rate; analysing the development of those values (55) per repetition over time; and creating a value over time pattern corresponding to the specific disturbance.
7 The method for determining a location and type of a disturbance (85) of claim 6, further comprising the step of comparing the value over time pattern with recorded patterns of known disturbance types.
8. The method for determining a location and type of a disturbance (85) of claims 6 or 7, wherein the values corresponding to the condition of the optical fibre are recorded at a specific repetition rate of 250 Hz.