EP1966580A2 - Detection d'une fuite de fluide dans un circuit par intercorrelation - Google Patents

Detection d'une fuite de fluide dans un circuit par intercorrelation

Info

Publication number
EP1966580A2
EP1966580A2 EP06847154A EP06847154A EP1966580A2 EP 1966580 A2 EP1966580 A2 EP 1966580A2 EP 06847154 A EP06847154 A EP 06847154A EP 06847154 A EP06847154 A EP 06847154A EP 1966580 A2 EP1966580 A2 EP 1966580A2
Authority
EP
European Patent Office
Prior art keywords
circuit
fluid
intercorrelation
measurement signals
leak
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
EP06847154A
Other languages
German (de)
English (en)
Inventor
Jaouad Zemmouri
Bruno Garnier
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Optical System and Research for Industry and Science OSYRIS SA
Original Assignee
Optical System and Research for Industry and Science OSYRIS SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Optical System and Research for Industry and Science OSYRIS SA filed Critical Optical System and Research for Industry and Science OSYRIS SA
Publication of EP1966580A2 publication Critical patent/EP1966580A2/fr
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00Investigating fluid-tightness of structures
    • G01M3/02Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/26Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors
    • G01M3/28Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds
    • G01M3/2807Investigating fluid-tightness of structures by using fluid or vacuum by measuring rate of loss or gain of fluid, e.g. by pressure-responsive devices, by flow detectors for pipes, cables or tubes; for pipe joints or seals; for valves ; for welds for pipes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M3/00Investigating fluid-tightness of structures
    • G01M3/02Investigating fluid-tightness of structures by using fluid or vacuum
    • G01M3/04Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point
    • G01M3/24Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using infrasonic, sonic, or ultrasonic vibrations
    • G01M3/243Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point using infrasonic, sonic, or ultrasonic vibrations for pipes

Definitions

  • the present invention relates to the automatic detection of fluid leaks (liquid or gas) in a circuit. It relates more particularly, but not exclusively, to the detection of a leakage of fluid that may have a very low flow rate in a circuit that may be subject to disturbances, in particular thermal disturbances, which have an influence on the flow rate of the fluid transported.
  • the invention advantageously allows a reliable detection and with a very good sensitivity of the appearance of a leak in a coolant circuit, such as a fluid circuit used to cool or to heat an enclosure.
  • a coolant circuit such as a fluid circuit used to cool or to heat an enclosure.
  • One of the preferred, but not exclusive, applications of the invention lies in the detection of water leaks in the cooling circuits of the enclosures of electric arc furnaces used in the iron and steel industry for the melting of scrap or the refining of steels.
  • Another application of the invention is the detection of leaks in pipelines for transporting and / or distributing fluids (pipelines, pipelines, etc.).
  • a method widely used to date to detect a leakage of fluid in a circuit consists in placing two flowmeters respectively at the output and the input of the circuit, and continuously calculating, by means of these two flowmeters, the difference between the inlet fluid flow rate and the output fluid flow rate. A leak is detected when this difference in flow rates exceeds a predefined threshold.
  • any disturbance of the circulating fluid resulting in a change in the volume of transported fluid is likely to falsify the detection (second disadvantage). It is for example and mainly a thermal disturbance (heating or cooling of the fluid), or a hydraulic disturbance (for example closure or opening of a pipe in a circuit forming a more or less complex network of several pipes in parallel). It is understood that in the case of an increase in the volume of fluid transported, for example under the effect of an increase in the temperature of the fluid between its input and its output of the circuit, the output flow rate becomes greater than the flow rate. input, which prevents the detection of any leak below this increase in flow.
  • the tank and more specifically the vault and the panels of arc furnaces which are used for the melting of scrap or the refining of steels.
  • the furnace tank is charged with the various ferrous materials to be melted, then the temperature of the scrap is raised to their melting point (typically between 1500 ° C. and 2000 ° C.) mainly by means of electric arcs generated inside the tank.
  • the molten steels are then poured out of the furnace for further processing.
  • the cooling of this type of furnace is obtained by means of at least one cooling circuit formed by a more or less complex network and more or less long tubing, inside which is circulated a liquid heat transfer fluid (generally water) which allows to effectively cool the tubings.
  • a liquid heat transfer fluid generally water
  • an arc furnace comprises in practice several independent cooling circuits, generally including at least one cooling circuit for the side panels of the tank and a cooling circuit for the vault of the furnace tank.
  • Temperature correction solutions are proposed by the manufacturers of flow meters, and are implemented in particular in the field of arc furnaces mentioned above.
  • a method for detecting nitrogen leakage in a circuit forming a barrier between the primary and secondary spaces for thermal insulation is provided.
  • a cryogenic tank said method being based on the differential nitrogen flow rate method at the outlet and the inlet of the circuit.
  • it is taught to provide temperature corrections and pressure corrections to flow measurements by measuring the temperature and pressure of the flowing fluid.
  • 0 188 911 a detection method which is similar to the aforementioned method of the flow differential, and which in this document is used to detect a leakage of fluid, and in particular a gas leak, in a pipe of the pipeline type. This method is essentially based on the calculation of the FD parameter below:
  • Q 1n is the amount of gas entering the pipeline and Q out being the amount of gas leaving the pipeline.
  • this publication it is also proposed, to reduce the errors on the measurement signal and improve the accuracy of the measurement, to increase the integration time and to remove the average value of the input signal over the period of time. integration by calculating the following FD * parameter:
  • This method of detection by simulation has the disadvantage of being very complicated, and unsuitable for complex fluid circuits comprising several possible fluid paths, such as for example the cooling circuits of iron-arc furnaces. is performed only at the end of the integration period, which disadvantageously induces a delay in the detection. But in a prejudicial way this delay is all the more important as the period of integration chosen is great for reasons of precision of the measurement.
  • German patent application DE 26 03 715 It has also been proposed in German patent application DE 26 03 715 to perform a leak detection, by measuring, for example, the pressure at the inlet of a circuit and the pressure at the outlet of the circuit, by withdrawing at each measurement. input and output the same reference value, which is calculated from the input and output signals, and calculating a cross-correlation of the signals.
  • this publication does not describe the means by which inter-correlation automatically detects a leak.
  • this method does not allow to obtain a sufficient sensitivity of detection, and at the Applicants know this relatively old method (published in 1977) would never have been used industrially.
  • the present invention has the general objective of proposing a new technical solution for automatically detecting leakage of a fluid (liquid and / or gas) flowing in a circuit.
  • Another more particular object of the invention is to propose a new technical solution which makes it possible to automatically detect leakage of a fluid (liquid and / or gas) flowing in a circuit, with improved detection sensitivity.
  • Another more particular objective of the invention is to propose a new technical solution which makes it possible to automatically detect a leak of a fluid (liquid and / or gas) flowing in a circuit subjected to disturbances (in particular thermal disturbances and / or hydraulic systems) affecting the volume of the flowing fluid.
  • Another more particular object of the invention is to propose a new technical solution which makes it possible to automatically detect a leakage of a fluid (liquid and / or gas) in flow in a circuit subjected to electromagnetic disturbances.
  • Another more particular object of the invention is to propose a new technical solution which makes it possible to automatically detect a leakage of a fluid (liquid and / or gas) flowing in a circuit subjected to acoustic disturbances.
  • Another more particular object of the invention is to propose a new technical solution which makes it possible to automatically detect a leakage of a fluid (liquid and / or gas) in flow in a complex circuit comprising several different possible paths for the flowing fluid. . Summary of the invention
  • the invention whose first object is a method of controlling a circuit in which a fluid circulates, said method comprising the following steps: a) the circulating fluid is detected by means of an input sensor and an output sensor respectively positioned at the input and at the output of the circuit, and respectively delivering two measurement signals E (t) and S (t) characteristics of the flow or fluid flow variations,
  • an offset correction is performed by removing the measurement signals [E (t) or E k ; S (t) or S k ] their mean value (Emoy;
  • the method comprises the following additional and optional features, taken alone or, where appropriate, in combination: in step (b) of processing the signals E (t) and S (t), the measurement signals E (t) and S (t) are sampled, and the cross-correlation function Icc is a function discrete calculated using the following formula:
  • E k are the values resulting from the sampling of the measurement signal E (t);
  • S k are the values resulting from the sampling of the measurement signal S (t); n is the total number of samples E k and S k in a predefined calculation window; k is an integer from 0 to (n-1); j is an integer taking the values between - (n-1) and (n-1);
  • Iccj is a vector of size 2n-1;
  • step (b) for processing the measurement signals an intercorrelation of the two measurement signals is calculated, minus the autocorrelation of the output measurement signal;
  • step (b) of processing the measurement signals in step (b) of processing the measurement signals, the amplitude variations of the central peak (P) of intercorrelation are detected; more particularly, the variations of the amplitude of the central cross-correlation peak (P) are compared with at least one predefined threshold, and the appearance of a leak is detected when said variation of the amplitude is negative and becomes, in absolute value, greater than this threshold;
  • said threshold is preferably auto-adaptive and is calculated as a function of the amplitudes of the previous central cross-correlation peaks; in another variant embodiment, in step (b) of treatment measurement signals, the integral of the inter-correlation is calculated, and the result of this integral is compared with a predefined threshold (s).
  • a predefined threshold s.
  • an intercorrelation of the two measurement signals (with or without prior offset correction) is calculated, and the amplitude variations of the central intercorrelation peak (P) are detected, or the integral of the cross-correlation is calculated, after possibly subtracting the autocorrelation of the output or input measurement signal.
  • the variations of the amplitude of the central intercorrelation peak (P) are compared with at least one predefined threshold, and detection is made the appearance of a leak when said variation of the amplitude is negative and becomes, in absolute value, greater than this threshold.
  • the result of this integral is compared with a threshold (s) predefined.
  • the invention also has for its second object the use of the method referred to above to detect the appearance of a leak in a cooling circuit or in a heating circuit, and more particularly, to detect the appearance of a leak in a cooling circuit of an electric arc furnace.
  • the third object of the invention is also the use of the method referred to above, for detecting the appearance of a leak in a pipe or a network of pipes for the transport and / or distribution of at least one fluid.
  • the invention also has for its fourth object an installation comprising at least one circuit, means for circulating a fluid in this circuit, an input sensor and an output sensor respectively positioned at the input and at the output. of the circuit, and delivering respectively two measurement signals E (t) and S (t) characteristic of the flow or fluid flow variations, and electronic means for processing the measurement signals E (t) and S (t).
  • Said electronic means (32) are designed to implement the measurement signal processing step (b) E (t) and (St) which is defined in the method referred to above.
  • the sensors are preferably non-invasive type flow meters.
  • the circuit is a cooling circuit or a heating circuit. More particularly, the installation is for example constituted by an electric arc furnace.
  • the circuit is a pipeline or pipe network for the transport and / or distribution of at least one fluid.
  • the fifth subject of the invention is a program recorded on a medium or in a memory and which, when executed by a programmable processing unit, performs the processing of two measurement signals E (t) and S (t). characteristics of the flow or fluid flow variations respectively at the inlet and the outlet of a circuit, said processing of the measurement signals E (t) and S (t) being carried out in accordance with step ( b) previously defined treatment.
  • FIG. 1 schematically represents an installation of the invention
  • FIG. 2 schematically represents an installation of the invention, which has been modified for experimental purposes to provoke a fluid leak
  • FIG. 3 schematically represents an installation of the invention, which has been modified for experimental purposes to cause thermal disturbances
  • FIG. 4 schematically represents an installation of the invention, which has been modified to experimental purposes to cause hydraulic disturbances
  • FIG. 5 is a functional diagram illustrating the main successive steps of a preferred variant embodiment of step (b) of processing the measurement signals according to the invention
  • FIGS. 6 to 26 are curves resulting from tests carried out by means of the installations of FIGS. 2 to 4, and are detailed below,
  • FIG. 27 is a block diagram illustrating the main successive steps of another variant embodiment of step (b) of processing the measurement signals according to the invention.
  • FIG. 1 diagrammatically shows a circuit 1, which is fed with a fluid (f) (symbolized by arrows), from a set 2 comprising a buffer tank, for the storage of this fluid, associated a pump for forced circulation of this fluid.
  • Circuit 1 comprises:
  • a main evacuation pipe 12 whose inlet is sealingly connected to the outlet of the network 11 common to all of the pipes 110, and the outlet of which is connected in a sealed manner at the inlet of the set 2 for forced circulation of the fluid (f).
  • the fluid (f) circulates in a closed circuit, being introduced into the main inlet pipe 10, circulates in the pipes 110 of the network 11, and is then conveyed at the outlet of the circuit 1 towards the assembly 2, via the evacuation pipe 12.
  • the circuit 1 comprises a complex network 11 of tubings 110 defining for the fluid (f) several possible paths between the input and the output of the circuit 1.
  • the invention is particularly interesting for automatically detecting a leak in this type of complex circuit 1, it is however not limited to this type of circuit, and can also be applied to a simpler single path circuit structure and constituted for example by a single tube.
  • said circuit is equipped with a control device 3 which is in accordance with the invention.
  • This control device 3 comprises:
  • a sensor 30 which is mounted at the inlet of the circuit 1, on the main intake duct 10, and which in operation delivers an electrical measurement signal E (t), - a sensor 31 which is mounted at the output of the circuit 1, on the discharge pipe 12, and which in operation delivers an electrical measurement signal S (t),
  • the electronic means 32 for processing the measurement signals E (t) and S (t) delivered by the sensors 30 and 31.
  • the sensors 30 and 31 are identical and may in general be constituted by any sensor, which allows to detect and characterize the flow or fluid flow variations.
  • a measurement signal E (t) or S (t) characterizes the quantity of fluid per unit of time passing through the sensor and / or characterizes the material waves generated by the fluid flow variations.
  • the measurement signals E (t) and S (t) delivered by the sensors 30 and 31 can be of analog type, or be of digital type when the sensors integrate an analog / digital converter.
  • the sensors 30 and 31 may be of the invasive type (for example measurement probe introduced inside the corresponding conduit 10 or 12) or of the non-invasive type (such as in the example of sensors illustrated in FIG. 1).
  • Non-invasive sensors are however preferred for the implementation of the invention, because the invasive sensors are detrimental to the source of a pressure drop on the fluid.
  • the sensors 30 and 31 are flow meters selected from the following list:
  • the first two types of flowmeters above are of the non-invasive type and are therefore preferred over the last two types of flowmeters which are of the invasive type.
  • the electromagnetic flowmeters are preferably used, because they have the best performances to date (precision, measurement range, robustness, etc.).
  • the electronic processing means 32 may be implemented in different forms, knowing that what is important for the invention lies in the method of processing the measurement signals E (t) and S (t) which is detailed below.
  • the electronic processing means 32 may be implemented:
  • a programmable processing unit such as a microcomputer executing a signal processing program
  • E (t) and S (t) according to the invention loaded in random access memory, or in the form of a specific electronic card whose electronic architecture comprises a microprocessor or microcontroller capable of executing an on-board signal processing program E (t) and S (t) according to the invention, or - in the form an electronic card comprising a specific electronic circuit ASIC type, specially designed to perform a signal processing program E (t) and S (t) according to the invention.
  • a specific electronic card whose electronic architecture comprises a microprocessor or microcontroller capable of executing an on-board signal processing program E (t) and S (t) according to the invention, or - in the form an electronic card comprising a specific electronic circuit ASIC type, specially designed to perform a signal processing program E (t) and S (t) according to the invention.
  • FIG. 5 shows the main steps of a preferred variant embodiment of the invention for processing the measurement signals E (t) and S (t) delivered by the sensors 30 and 31.
  • the signal processing is based on an intercorrelation of the measurement signals E (t) and S (t) [block 43] and on a detection [block 44] and a follow-up [block 45] of the central peak ( P) resulting from the cross-correlation of the measurement signals.
  • the intercorrelation function is a mathematical function that is known per se. In the continuous time domain, this intercorrelation function lcc (t), when applied to the measurement signals E (t) and S (t), is defined by the relation:
  • the cross-correlation ultimately consists of calculating the overlap integral between the signals E (t) and S (t) by temporally shifting one signal relative to the other.
  • the intercorrelation which is calculated in the algorithm of FIG. 5 [block 43] is a discretized intercorrelation function.
  • the discrete intercorrelation function I C q implemented is preferably defined by the following relation:
  • S k are the values resulting from the sampling of the measurement signal S (t) [vector of size n]; n is the total number of samples E k and S k in a predefined calculation window; k is an integer from 0 to (n-1); j is an integer taking the values between - (n-1) and (n-1) and thus takes the following successive values: - (n-1); - (n-2); ...; -2; -1; 0; 1; 2; ...; (n-1);
  • Iccj is a vector of size 2n-1.
  • the continuous measurement signals E (t) and S (t) are sampled with a predefined sampling frequency (fe).
  • sampling operations are usually carried out using analog / digital converters. These converters can be integrated in the sensors 30, 31 or can be integrated in the electronic processing means 32.
  • Blocks (41) and (42) Offset Correction When the system speed is set, the input flow rate E (t) and output S (t) signals oscillate around a mean value.
  • This average value (Emoy and Smoy) is calculated in parallel for each signal [block 41].
  • the offset correction [block 42] consists in subtracting from each sample E k the average value Emoy of the corresponding signal, and each sample Sk the average value Smoy of the corresponding signal. We obtains the corrected signals: (E k - Emoy) and (S k - Smoy)
  • Emoy and Smoy average values are calculated once and for all and stored in memory.
  • these average values Emoy and Smoy are recalculated as and when in a sliding time window of predefined duration T moy .
  • the value of this duration T m0 y must be greater than the duration of occurrence of a leak in the circuit, and must be sufficiently small to take into account drift and / or disturbances of the system.
  • This duration T moy is fixed case by case by the skilled person.
  • this duration is chosen so that in practice the corrected signals (E k - Emoy) and (S k - Smoy) oscillate substantially around zero.
  • This offset correction step is important because it makes it possible, in combination with the cross-correlation of the corrected signals (block 43), to improve the sensitivity of the leak detection.
  • E k and S k are those directly derived from the sampling (variant embodiment without offset correction) or are the corrected samples obtained after subtraction of the average value of the signal (embodiment variant with offset correction).
  • this offset correction could be performed on the analog measurement signals E (t) and S (t), (before sampling) by means of an analog subtractor.
  • the computation of the intercorrelation is carried out with 1000 successive samples E k and S k (k varying from 1 to n and n being 1000), with a sampling frequency fe of 1 kHz.
  • an intercorrelation vector Icq is thus calculated every second.
  • n samples Ek and S k are performed repetitively in successive sliding windows of n samples Ek and S k .
  • These successive windows (or series) of n samples E k or S k can be non-overlapping in time (we take n samples, then the following n samples without overlapping between the series of samples), or can on the contrary overlap in part.
  • the two calculation windows of n samples E k and S k that are used can be defined on the same time interval, without time offset between the windows (in this case the samples E k and S k successive of each window were all sampled at the same times), or on the contrary can be defined with a more or less significant temporal shift between the calculation windows.
  • the sampling frequency (fe) and the number n of samples characterize a parameter T, corresponding to the integration period for the discretized intercorrelation.
  • Block 45 Tracking and alarm
  • the variations of the central peak P are controlled by comparing the amplitude of the central peak (P) of the intercorrelation with at least one threshold (s).
  • This threshold (s) may, depending on the case, be a predefined threshold of constant value.
  • the self-adaptive threshold which is a function of the values E k and S k , and more precisely successive values calculated for the central peak of the intercorrelation; for example, the self-adaptive threshold corresponds to the average value over a predefined number N of previously calculated central peak values.
  • a leak is automatically detected when the variation in the amplitude of the central peak of the intercorrelation is negative and is in absolute value greater than a predefined threshold (s), and an alarm (audible, visual, sending) is triggered if necessary. automatic alarm message by any known telecommunication means, etc.)
  • the variation of amplitude of the central peak is negative and exceeds in absolute value one of the thresholds (si) to (s3), one triggers a characteristic alarm of this threshold.
  • the fluid (f) flowing in the circuit 1 is a liquid (in this case water); the sensors 30 and 31 used are electromagnetic flowmeters having an accuracy of the order of 0.5%.
  • the duration T moy for the offset correction was 30 seconds.
  • the sampling frequency (fe) of the measurement signals E (t) and S (t) is equal to 1 KHz.
  • the number of samples (n) for intercorrelation is equal to 1000.
  • the discretized intercorrelation function Icq is calculated on (n) samples Ek and Sk acquired simultaneously and not temporally shifted, and in sliding windows of n samples which are juxtaposed and do not overlap.
  • Leak Detection FIG. 2 shows a fluid circuit 1 modified for purely experimental purposes to cause fluid leakage.
  • This circuit 1 of FIG. 2 differs from that of FIG. 1 only in that on one of the tubes 110 of this circuit, a bypass 13 has been made, on which a leakage valve 14 is mounted. a receptacle 15. When the valve 14 is open, a very small portion of the fluid (f) is taken from the circuit 1 and is not redirected at the output of the circuit 1 to the assembly 2, but flows into the receptacle 15 , which makes it possible to cause a leak in the circuit 1 in the context of experiments of the invention.
  • these means 13, 14, 15 for causing a leak in the circuit 1 are provided solely for experimental purposes to test the invention, and are not found in the context of a final operational and non-experimental installation comprising the circuit 1 and its means 2 for supplying fluid.
  • FIG. 2 In a first step, the installation of FIG. 2 was operated without leakage, by supplying the circuit 1 with water with an average flow rate of the order of 295 m 3 / h, the leakage valve 14 being closed.
  • FIGS. 6 and 7 show, over a short period of operation time without leakage, the raw data ER and Sk coming from the sampling of the measurement signals E (t) and S (t) (signals at the output of FIG. block 42 after offset correction).
  • FIG. 8 shows the cross-correlation of the sampled measurement signals ER and S k of FIGS. 6 and 7.
  • a leak with a flow rate of the order of 1.2 m 3 / h is caused by opening the valve 14.
  • FIGS. 9 and 10 show, for a short period of operating time taken at the beginning of leakage, the raw data E k and S k coming from the sampling of the measurement signals E (t) and S (t). (output signals of block 42 after offset correction).
  • FIG. 11 shows the inter-correlation of the sampled measurement signals E k and S k of FIGS. 9 and 10.
  • FIGS. 12 and 13 show the measurement signals E (t) and S (t) delivered by the sensors 30 and 31 for the period of time between 470s and 570s.
  • FIG. 14 shows the time course of the amplitude of the central peaks of the intercorrelation (output of block 44) for the period of time between 0s and 730s.
  • the parts of the curve for discriminating the five successive leaks are indicated by the arrows F1, F2, F3, F4 and F5.
  • FIG. 14 makes it possible to show a sudden negative variation in the amplitude of the central peak of the intercorrelation at each occurrence of a leak. It also makes it possible to show that this amplitude variation, which makes it possible to characterize and detect the appearance of a leak, increases in absolute value with the flow rate of the leak.
  • FIG. 15 shows the amplitude variation (in absolute value) of the central cross-correlation peak characteristic of a leak as a function of the relative value of the leakage caused (in m 3 / h). This figure shows that the variation of the amplitude of the central cross-correlation peak is proportional to the relative value of the leak.
  • the intercorrelation can therefore not only be used to detect the appearance of a leak in the circuit 1, but also advantageously makes it possible to characterize the relative value of the leak that has been detected.
  • FIG. 3 shows a fluid circuit 1 modified for experimental purposes only to cause thermal disturbances in circuit 1.
  • This circuit 1 of FIG. 3 differs from that of FIG. 2 only in that added on one of the tubes 110 of the circuit 1:
  • a heating element 16 for example a heating resistor which makes it possible to locally heat the tubing 110 and thereby the portion of fluid circulating locally in this tubing 110;
  • FIG. 3 In a first phase, the installation of FIG. 3 is operated without causing a leak (valve 14 closed).
  • FIGS. 16 and 17 respectively show the measurement signals E (t) and S (t) during this test with fluid heating and without leakage, and in FIG. 18 the temperature of the fluid as a function of time.
  • FIG. 19 shows the evolution over time of the amplitude of the central cross-correlation peaks. This curve of FIG. 19 shows that, at the start of heating, the variation in the amplitude of the central intercorrelation peaks increases sharply, but in the opposite direction (positive variation) to the variation of the amplitude of the central peaks of FIG. intercorrelation in case of leakage.
  • a second phase the heating of the tubing is stopped 110 for a time sufficient to return to substantially identical fluid temperature conditions throughout circuit 1.
  • FIG. 20 shows the temperature of the fluid as a function of time (substantially constant) measured upstream of the heating element 16 by the sensor 17a (curve A) and the temperature of the fluid as a function of the time measured downstream of the heating element 16 by the sensor 17b (curve B) at the end of the second phase and during the third heating phase.
  • FIGS. 21 and 22 respectively show the measurement signals E (t) and S (t) during this leakage heating test, and in FIG. 23 the evolution over time of the amplitude of the central peaks of FIG. intercorrelation.
  • the curve 23 there is a sudden negative variation in the amplitude of the central peaks of the intercorrelation very shortly after the outbreak of the leak. This sudden variation in the amplitude of the central intercorrelation peaks thus makes it possible to always detect the appearance of the leak, despite the local thermal disturbance in the circuit 1, which was caused by the heating of the tubing 110 by means of the heating element 16. Hydraulic disturbances
  • FIG. 4 shows a fluid circuit 1 modified for purely experimental purposes only to cause hydraulic disturbances in the circuit 1.
  • This circuit 1 of FIG. 4 differs from that of FIG. 2 only in that has added on one of the tubings 110 of the circuit 1 a valve 18 for closing or opening said tubing 110.
  • This curve of FIG. 26 shows that very shortly after the closing of the valve 18 (hydraulic disturbance in the circuit 1), the variation in the amplitude of the central intercorrelation peaks increases sharply, but in the opposite direction (positive variation ) the variation of the amplitude of the central intercorrelation peaks in case of leakage. This behavior is comparable to what has been previously described for thermal disturbances, with particular reference to FIG. 19.
  • the detection of the central peak of correlation (peak of larger amplitude) (block 44 of FIG. 5) can be replaced by a calculation of the sum of all the samples I CC j resulting from the intercorrelation (which amounts to calculating the integral of the intercorrelation).
  • this sum is close to zero.
  • the result of this sum becomes strongly negative.
  • a negative threshold (s) is chosen in this case.
  • the detection (block 45) is carried out either directly on the result of the intercorrelation, but on the result of the reduced intercorrelation (block 47) of the result of the autocorrelation (block 46) of the corrected output signal (S k -S avg ).
  • the invention finds its application in the detection of the appearance of fluid leak (s) in any circuit within which a fluid circulates, said fluid may be a liquid or a gas, or a mixture gas / liquid.
  • the invention is of interest in all applications: where the flow rate of the leaks to be detected can be very small compared to the flow rate of the fluid at the inlet of the circuit, and / or
  • a particularly interesting example of application of the invention lies in the detection of leaks in a cooling circuit or in a heating circuit, inside which circulates a coolant.
  • the invention may advantageously be used in the field of iron and steel to automatically detect reliably and quickly the appearance of leaks in the cooling circuits of an electric arc furnace which form a complex network of tubing.
  • the circuit 1 of Figure 1 is for example the cooling circuit fitted to the vault of an electric arc furnace, or the cooling circuit panels of an electric arc furnace.
  • a particularly interesting example of application of the invention lies in the detection of leaks in a pipe or in a network of pipes for the transport and / or distribution of fluids (pipelines, pipelines, etc.).

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Examining Or Testing Airtightness (AREA)

Abstract

On contrôle automatiquement un circuit dans lequel circule un fluide, en mettant en oevre les étapes suivantes : (a) on détecte le fluide en circulation au moyen d'un capteur d'entrée et d'un capteur de sortie positionnés respectivement à l'entrée et à la sortie du circuit, et délivrant respectivement deux signaux de mesure E(t) et S(t) caractéristiques de l'écoulement ou des variations d'écoulement du fluide, (b) on traite les deux signaux de mesure E(t) , S(t) dans le but de détecter le cas échéant une fuite de fluide dans le circuit. A l'étape (b) de traitement des signaux E(t) et S(t), on effectue une correction offset en retirant au signaux de mesure [E(t) ou E<SUB>k</SUB> ; S(t) ou S<SUB>k</SUB>] leur valeur moyenne (Emoy ; Smoy), et on calcule une intercorrélation des deux signaux de mesure corrigés. Plus particulièrement, dans une variante de réalisation on détecte les pics centraux (P) d'intercorrélation, et on compare les variations de l'amplitude du pic central d'intercorrélation (P) avec au moins un seuil prédéfini (s). On détecte l'apparition d'une fuite lorsque ladite variation de l'amplitude est négative et devient, en valeur absolue, supérieure à ce seuil.

Description

DETECTION D'UNE FUITE DE FLUIDE DANS UN CIRCUIT PAR
INTERCORRELATION
Domaine technique La présente invention concerne la détection automatique de fuites de fluide (liquide ou gaz) dans un circuit. Elle concerne plus particulièrement, mais non exclusivement, la détection d'une fuite de fluide pouvant présenter un très faible débit dans un circuit pouvant être sujet à des perturbations, notamment des perturbations thermiques, qui ont une influence sur le débit de fluide transporté. Par exemple, l'invention permet avantageusement une détection fiable et avec une très bonne sensibilité de l'apparition d'une fuite dans un circuit de fluide caloporteur, tel qu'un circuit de fluide utilisé pour refroidir ou pour chauffer une enceinte. Une des applications préférentielles, mais non exclusive, de l'invention réside dans la détection de fuites d'eau dans les circuits de refroidissement des enceintes des fours électriques à arcs utilisés en sidérurgie pour la fusion des ferrailles ou l'affinage des aciers. Une autre application de l'invention est la détection de fuites dans des canalisations de transport et/ou de distribution de fluides (pipelines, gazoducs, etc..) Art antérieur
Une méthode largement utilisée à ce jour pour détecter une fuite de fluide dans un circuit, tel que par exemple un circuit de refroidissement ou de chauffage d'une enceinte, consiste à placer deux débitmètres respectivement à la sortie et à l'entrée du circuit, et à calculer en permanence, au moyen de ces deux débitmètres, la différence entre le débit de fluide en entrée et le débit de fluide en sortie. Une fuite est détectée lorsque cette différence de débits excède un seuil prédéfini.
Cette méthode désignée ci-après « méthode du différentiel des débits » présente au moins trois inconvénients majeurs. Elle est limitée par la précision de la mesure des débitmètres utilisés
(1er inconvénient). Moins la précision des débitmètres utilisés est bonne, et plus le débit de la fuite de fluide doit être important pour que celle-ci puisse être détectée. Par exemple, et à titre indicatif les débitmètres électromagnétiques actuellement commercialisés ont une précision de l'ordre de 0,3 % à 1%. Ainsi, avec ce type de débitmètres, la méthode du différentiel des débits peut être utilisée à ce jour au mieux pour détecter des débits de fuites qui sont au moins supérieurs selon le cas de 0,6% à 2% du débit de circulation du fluide dans le circuit.
Egalement, toute perturbation du fluide en circulation se traduisant par une modification du volume de fluide transporté est susceptible de venir fausser la détection (deuxième inconvénient). Il s'agit par exemple et principalement d'une perturbation thermique (chauffage ou refroidissement du fluide), ou encore d'une perturbation hydraulique (par exemple fermeture ou ouverture d'une canalisation dans un circuit formant un réseau plus ou moins complexe de plusieurs canalisations en parallèle). On comprend en effet que dans le cas d'une augmentation du volume de fluide transporté, par exemple sous l'effet d'une augmentation de la température du fluide entre son entrée et sa sortie du circuit, le débit de sortie devient supérieur au débit d'entrée, ce qui empêche de détecter toute fuite inférieure à cette augmentation de débit. A l'inverse, dans le cas d'une diminution du volume de fluide transporté, par exemple sous l'effet d'une forte baisse de la température du fluide entre son entrée et sa sortie du circuit, le débit de sortie devient inférieur au débit d'entrée, ce qui peut déclencher une fausse détection de fuite, alors qu'aucune fuite n'est présente.
La présence de perturbations électromagnétiques peut également perturber les signaux de mesure de débit, et de ce fait fausser la détection de fuite (3eme inconvénient).
D'une manière générale, les trois inconvénients précités rendent cette méthode du différentiel des débits peu fiable et le cas échéant inadaptée pour détecter des fuites de fluide présentant un très faible débit et/ou des fuites d'un fluide soumis à des perturbations, notamment thermiques, hydrauliques, affectant le volume de fluide transporté, ou électromagnétiques.
Or il existe de nombreuse applications industrielles dans lesquelles un fluide est en circulation dans un circuit, et il est nécessaire de détecter de manière fiable et de manière suffisamment précoce toute fuite, même très faible, et de surcroît dans un environnement dans lequel le fluide peut subir des perturbations non stationnaires, et notamment des perturbations thermiques plus ou moins importantes dans le temps.
Parmi ces nombreuses applications industrielles, on peut plus particulièrement citer toutes les applications industrielles dans lesquelles un fluide circule dans un circuit de refroidissement ou dans un circuit de chauffage d'une installation.
Plus particulièrement, dans le domaine de la sidérurgie, il est usuel de refroidir la cuve et plus spécifiquement la voûte et les panneaux des fours à arcs, qui sont utilisés pour la fusion des ferrailles ou l'affinage des aciers. Dans ce type de four, on charge la cuve du four avec les différents matériaux ferreux à fondre, puis la température des ferrailles est élevée jusqu'à leur point de fusion (typiquement entre 15000C et 2000°c) principalement au moyen d'arcs électriques générés à l'intérieur de la cuve. Les aciers en fusion sont ensuite coulés en dehors du four en vue de leur traitement ultérieur.
Le refroidissement de ce type de four est obtenu au moyen d'au moins un circuit de refroidissement formé par un réseau plus ou moins complexe et plus ou moins long de tubulures, à l'intérieur desquelles est mis en circulation forcée un liquide caloporteur (généralement de l'eau) qui permet de refroidir efficacement les tubulures. Typiquement, un four à arcs comporte en pratique plusieurs circuits de refroidissement indépendants, dont généralement au moins un circuit de refroidissement pour les panneaux latéraux de la cuve et un circuit de refroidissement pour la voûte de la cuve du four. Ce type de four sidérurgique est particulièrement vulnérable aux fuites de liquide dans son circuit de refroidissement, pour au moins deux raisons :
- si à cause d'une fuite dans le circuit de refroidissement, de l'eau ruisselle à l'intérieur du four, elle peut atteindre le fond du four qui est garni de briques réfractaires ; dans ce cas, les briques réfractaires se détériorent très rapidement, et cette détérioration aboutit à un percement du fond de la cuve du four ; dans ce cas, le métal en fusion s'échappe hors de la cuve avec potentiellement de graves conséquences humaines et financières ;
- si à cause d'une fuite dans le circuit de refroidissement, de l'eau tombe à l'intérieur de la cuve d'un four en fonctionnement, et vient au contact du métal en fusion, elle peut se décomposer en oxygène et en hydrogène dans certaines conditions à cause de la température très élevée. Cette décomposition engendre un risque important d'explosion accidentelle. On comprend en conséquence que dans ce type particulier d'application les enjeux pour trouver une détection automatique de fuite qui soit fiable et précoce sont extrêmement importants, tant sur le plan humain pour éviter des accidents pouvant être dans certains cas mortels, que sur le plan financier pour éviter une immobilisation longue et coûteuse du four en cas d'accident.
Aujourd'hui, les méthodes de détection de fuite qui sont mises en oeuvre industriellement sur ces fours sidérurgiques sont essentiellement des méthodes basées sur une mesure du différentiel des débits d'entrée et de sortie, et ne donnent donc pas entièrement satisfaction compte tenu des inconvénients précités inhérents à ce type de méthode de détection.
Pour essayer de corriger les effets des perturbations thermiques, sur le volume de fluide transporté, et tenter de pallier aux problèmes précités de fausse détection ou de détection manquée qui peuvent en découler, des perfectionnements de la méthode précitée du différentiel des débits ont été proposés à ce jour. Ces perfectionnements sont basés sur une correction des mesures des débitmètres à partir notamment de la température ou de la pression du flux de fluide.
Des solutions de correction de température sont notamment proposées par les fabricants de débitmètres, et sont mises en oeuvre notamment dans le domaine des fours à arcs précités. Par exemple également, dans un autre domaine d'application décrit dans la demande de brevet français FR 2 509 839, on propose une méthode de détection de fuite d'azote dans un circuit formant une barrière entre les espaces primaire et secondaire de calorifugeage d'un réservoir cryogénique, ladite méthode étant basée sur la méthode du différentiel des débits d'azote en sortie et entrée du circuit. Dans cette publication, il est enseigné d'apporter des corrections de température et des corrections de pression aux mesures de débit en mesurant la température et la pression du fluide en écoulement.
D'une part ce type de corrections n'est en définitive pas très fiable et impose en outre un calibrage des mesures de température et de pression qui peut être délicat à réaliser, et qui doit être effectué régulièrement. D'autre part et surtout, ces corrections de température ou de pression ne permettent pas de pallier au premier inconvénient ci-dessus lié à Ia faible précision de mesure des débitmètres. II est par ailleurs proposé dans la demande de brevet européen EP
0 188 911 une méthode de détection qui s'apparente à la méthode précitée du différentiel des débits, et qui dans ce document est utilisée pour détecter une fuite de fluide, et notamment une fuite de gaz, dans une canalisation du type pipeline. Cette méthode est basée essentiellement sur le calcul du paramètre FD ci-après :
FD
Q1n étant la quantité de gaz qui entre dans le pipeline et Qout étant la quantité de gaz qui sort du pipeline. Dans cette publication, il est également proposé, pour réduire les erreurs sur le signal de mesure et améliorer la précision de la mesure, d'augmenter le temps d'intégration et de retirer la valeur moyenne du signal d'entrée sur la période d'intégration en calculant le paramètre FD* suivant :
Le calcul du paramètre FD ou du paramètre FD* présente tous les inconvénients précités de la méthode du différentiel des débits.
Par ailleurs, dans cette publication, il est enseigné de détecter une fuite lorsque le paramètre calculé (FD ou FD*) dépasse une valeur seuil calculée par simulation. Cette simulation prend en compte de nombreux paramètres dont le diamètre de la canalisation dans laquelle circule le fluide, et les éventuels changements de section, la pression du flux de fluide, les frictions du fluide, la gravité,...
Cette méthode de détection par simulation présente l'inconvénient d'être très compliquée, et inadaptée à des circuits de fluide complexes comportant plusieurs chemins de fluides possibles, tels que par exemple les circuits de refroidissement des fours sidérurgiques à arcs En outre la détection de fuite n'est réalisée qu'à la fin de la période d'intégration, ce qui de manière désavantageuse induit un retard dans la détection. Or de manière préjudiciable ce retard est d'autant plus important que la période d'intégration choisie est grande pour des raisons de précision de la mesure.
On a également proposé dans la demande de brevet allemand DE 26 03 715 de réaliser une détection de fuite, en mesurant, par exemple, la pression à l'entrée d'un circuit et la pression en sortie du circuit, en retirant à chaque mesure d'entrée et de sortie une même valeur de référence, qui est calculée à partir des signaux d'entrée et de sortie, et en calculant une intercorrélation des signaux. D'une part, dans cette publication, on ne décrit pas le moyen permettant à partir de Pintercorrélation de détecter automatiquement une fuite. D'autre part, il semble que cette méthode ne permette pas d'obtenir une sensibilité de détection suffisante, et à la connaissance des demanderesses cette méthode relativement ancienne (publiée en 1977) n'aurait jamais été utilisée industriellement.
On a également déjà proposé à ce jour de localiser une fuite de fluide dans un circuit, en détectant le bruit de la fuite au moyen de trois capteurs acoustiques (qui ne permettent pas de détecter l'écoulement ou les variations d'écoulement du fluide dans le circuit), en calculant des intercorrélations des signaux de mesure acoustique, et en déterminant la position physique de la fuite à partir des positions relatives des pics d'intercorrélation de plus grande amplitude. Cette méthode de localisation d'une fuite est décrite par exemple dans le brevet US 5 544 074. Cette méthode présente au moins deux inconvénients majeures. Elle peut être mise en oeuvre uniquement dans un circuit de fluide simple comportant un unique chemin pour le fluide, et cette détection de fuite est inutilisable dans des circuits soumis à de fortes perturbations acoustiques. Objectifs de l'invention
La présente invention a pour objectif général de proposer une nouvelle solution technique permettant de détecter automatiquement une fuite d'un fluide (liquide et/ou gaz) en écoulement dans un circuit.
Un autre objectif plus particulier de l'invention est de proposer une nouvelle solution technique qui permet de détecter automatiquement une fuite d'un fluide (liquide et/ou gaz) en écoulement dans un circuit, avec une sensibilité de détection améliorée.
Un autre objectif plus particulier de l'invention est de proposer une nouvelle solution technique qui permet de détecter automatiquement une fuite d'un fluide (liquide et/ou gaz) en écoulement dans un circuit soumis à des perturbations (notamment perturbations thermiques et/ou hydrauliques) affectant le volume du fluide en écoulement.
Un autre objectif plus particulier de l'invention est de proposer une nouvelle solution technique qui permet de détecter automatiquement une fuite d'un fluide (liquide et/ou gaz) en écoulement dans un circuit soumis à des perturbations électromagnétiques. Un autre objectif plus particulier de l'invention est de proposer une nouvelle solution technique qui permet de détecter automatiquement une fuite d'un fluide (liquide et/ou gaz) en écoulement dans un circuit soumis à des perturbations acoustiques. Un autre objectif plus particulier de l'invention est de proposer une nouvelle solution technique qui permet de détecter automatiquement une fuite d'un fluide (liquide et/ou gaz) en écoulement dans un circuit complexe comportant plusieurs trajets différents possibles pour le fluide en écoulement. Résumé de l'invention
Tout ou partie de ces objectifs, dont au moins l'objectif général précité, est atteint par l'invention, qui a pour premier objet un procédé de contrôle d'un circuit dans lequel circule un fluide, ledit procédé comportant les étapes suivantes : (a) on détecte le fluide en circulation au moyen d'un capteur d'entrée et d'un capteur de sortie positionnés respectivement à l'entrée et à la sortie du circuit, et délivrant respectivement deux signaux de mesure E(t) et S(t) caractéristiques de l'écoulement ou des variations d'écoulement du fluide,
(b) on traite les deux signaux de mesure E(t) , S(t) dans le but de détecter le cas échéant une fuite de fluide dans le circuit.
Selon un premier aspect de l'invention, à l'étape (b) de traitement des signaux E(t) et S(t), on effectue une correction offset en retirant aux signaux de mesure [E(t) ou Ek ; S(t) ou Sk] leur valeur moyenne (Emoy ;
Smoy), et. on calcule une intercorrélation des deux signaux de mesure corrigés.
Plus particulièrement, le procédé comporte les caractéristiques additionnelles et facultatives ci-après, prises isolément ou le cas échéant en combinaison : - à l'étape (b) de traitement des signaux E(t) et S(t), on réalise un échantillonnage des signaux de mesure E(t) et S(t), et la fonction d'intercorrélation Icc est une fonction discrète calculée au moyen de la formule suivante :
où :
Ek sont les valeurs résultant de l'échantillonnage du signal de mesure E(t);
Sk sont les valeurs résultant de l'échantillonnage du signal de mesure S(t); n est le nombre d'échantillons total Ek et Sk dans une fenêtre de calcul prédéfinie ; k est un nombre entier compris entre 0 et (n-1) ; j est un nombre entier prenant les valeurs entre -(n-1) et (n- 1) ;
Iccj est un vecteur de taille 2n-1 ;
- dans une variante de réalisation, à l'étape (b) de traitement des signaux de mesure, on calcule une intercorrélation des deux signaux de mesure, diminuée de l'autocorrélation du signal de mesure de sortie ;
- dans une variante de réalisation, à l'étape (b) de traitement des signaux de mesure, on détecte les variations d'amplitude du pic central (P) d'intercorrélation ; plus particulièrement, on compare les variations de l'amplitude du pic central d'intercorrélation (P) avec au moins un seuil prédéfini, et on détecte l'apparition d'une fuite lorsque ladite variation de l'amplitude est négative et devient, en valeur absolue, supérieure à ce seuil ;
- ledit seuil est de préférence auto-adaptatif et est calculé en fonction des amplitudes des pics centraux d'intercorrélation précédents ; - dans une autre variante de réalisation, à l'étape (b) de traitement des signaux de mesure, on calcule l'intégrale de Pintercorrélation, et on compare le résultat de cette intégrale avec un seuil (s) prédéfini. Selon un deuxième aspect de l'invention, à l'étape (b) de traitement des signaux E(t) et S(t), on calcule une intercorrélation des deux signaux de mesure (avec ou sans correction d'offset préalable), et on détecte les variations d'amplitude du pic central (P) d'intercorrélation, ou on calcule l'intégrale de l'intercorrélation, après éventuellement soustraction de l'autocorrélation du signal de mesure de sortie ou d'entrée. Plus particulièrement, dans Ia variante avec détection des variations d'amplitude du pic central (P) d'intercorrélation, on compare les variations de l'amplitude du pic central d'intercorrélation (P) avec au moins un seuil prédéfini, et on détecte l'apparition d'une fuite lorsque ladite variation de l'amplitude est négative et devient, en valeur absolue, supérieure à ce seuil.
Plus particulièrement, dans la variante avec calcul de l'intégrale de l'intercorrélation, pour détecter l'apparition d'une une fuite, on compare le résultat de cette intégrale avec un seuil (s) prédéfini.
L'invention à également pour deuxième objet l'utilisation du procédé visé précédemment pour détecter l'apparition d'une fuite dans un circuit de refroidissement ou dans un circuit de chauffage, et plus particulièrement, pour détecter l'apparition d'une fuite dans un circuit de refroidissement d'un four électrique à arcs.
L'invention a également pour troisième objet l'utilisation du procédé visé précédemment, pour détecter l'apparition d'une fuite dans une canalisation ou un réseau de canalisations de transport et/ou de distribution d'au moins un fluide.
L'invention a également pour quatrième objet une installation comportant au moins un circuit, des moyens de mise en circulation d'un fluide dans ce circuit, un capteur d'entrée et un capteur de sortie positionnés respectivement à l'entrée et à la sortie du circuit, et délivrant respectivement deux signaux de mesure E(t) et S(t) caractéristiques de l'écoulement ou des variations d'écoulement du fluide, et des moyens électroniques de traitement des signaux de mesure E(t) et S(t). Lesdits moyens électroniques (32) sont conçus pour mettre en œuvre l'étape (b) de traitement de signaux de mesure E(t) et (St) qui est définie dans le procédé visé précédemment.
Dans une variante particulière de réalisation, les capteurs sont de préférence des débitmètres de type non invasif.
Dans une application particulière de l'invention, le circuit est un circuit de refroidissement ou un circuit de chauffage. Plus particulièrement, l'installation est par exemple constituée par un four électrique à arcs.
Dans une autre application, le circuit est une canalisation ou un réseau de canalisations pour le transport et/ou la distribution d'au moins un fluide. L'invention a pour cinquième objet un programme enregistré sur un support ou dans une mémoire, et qui, lorsqu'il est exécuté par une unité de traitement programmable, réalise le traitement de deux signaux de mesure E(t) et S(t) caractéristiques de l'écoulement ou des variations d'écoulement du fluide respectivement à l'entrée et à la sortie d'un circuit, ledit traitement des signaux de mesure E(t) et S(t) étant réalisé conformément à l'étape (b) de traitement définie précédemment.
Brève description des dessins
D'autres caractéristiques et avantages de l'invention apparaîtront plus clairement à la lecture de la description détaillée ci-après de plusieurs variantes préférées de réalisation de l'invention, laquelle description est donnée à titre d'exemple non limitatif et non exhaustif de l'invention, et en référence aux dessins annexés sur lesquels :
- la figure 1 représente de manière schématique une installation de l'invention, - la figure 2 représente de manière schématique une installation de l'invention, qui été modifiée à des fins expérimentales pour provoquer une fuite de fluide,
- la figure 3 représente de manière schématique une installation de l'invention, qui a été modifiée à des fins expérimentales pour provoquer des perturbations thermiques, - la figure 4 représente de manière schématique une installation de l'invention, qui a été modifiée à des fins expérimentales pour provoquer des perturbations hydrauliques,
- La figure 5 est un schéma fonctionnel illustrant les principales étapes successives d'une variante préférée de réalisation de l'étape (b) de traitement des signaux de mesure selon l'invention,
- les figures 6 à 26 sont des courbes résultant de tests réalisés au moyen des installations des figures 2 à 4, et sont détaillées ci-après,
- la figure 27 est un schéma fonctionnel illustrant les principales étapes successives d'une autre variante de réalisation de l'étape (b) de traitement des signaux de mesure selon l'invention.
Description détaillée
On a représenté de manière schématique sur la figure 1 un circuit 1 , qui est alimenté avec un fluide (f) (symbolisé par des flèches), à partir d'un ensemble 2 comprenant un réservoir tampon, pour le stockage de ce fluide, associé à une pompe de mise en circulation forcée de ce fluide. Le circuit 1 comporte :
- une conduite d'admission principale 10, dont l'entrée est raccordée de manière étanche à la sortie de l'ensemble 2 (réservoir/pompe),
- un réseau 11 complexe de tubulures 110 en parallèle, dont l'entrée commune à toutes les tubulures 110, est raccordée de manière étanche à la sortie de la conduite d'admission principale 10, lesdites tubulures 110 pouvant présenter des caractéristiques géométriques différentes (longueurs, diamètres, etc.),
- une conduite d'évacuation principale 12, dont l'entrée est raccordée de manière étanche à la sortie du réseau 11 commune à toutes les de tubulures 110, et dont la sortie est raccordée de manière étanche à l'entrée de l'ensemble 2 de mise en circulation forcée du fluide (f).
Lorsque la pompe de l'ensemble 2 fonctionne, le fluide (f) circule en circuit fermé, en étant introduit dans la conduite d'admission principale 10, circule dans les tubulures 110 du réseau 11, puis est acheminé, en sortie du circuit 1 en direction de l'ensemble 2, via la conduite d'évacuation 12.
Dans l'exemple particulier illustré, le circuit 1 comporte un réseau 11 complexe de tubulures 110 définissant pour le fluide (f) plusieurs chemins possibles entre l'entrée et la sortie du circuit 1. Bien que l'invention soit particulièrement intéressante pour détecter automatiquement une fuite dans ce type de circuit 1 complexe, elle n'est toutefois pas limitée à ce type de circuit, et peut également être appliquée à un circuit de structure plus simple à chemin unique et constitué par exemple par un simple tube.
Dans le but de détecter automatiquement une fuite de fluide (f) dans le circuit 1, ledit circuit est équipé d'un dispositif de contrôle 3 qui est conforme à l'invention.
Ce dispositif de contrôle 3 comporte :
- un capteur 30 qui est monté à l'entrée du circuit 1 , sur la conduite principale d'admission 10, et qui en fonctionnement délivre un signal électrique de mesure E(t), - un capteur 31 qui est monté en sortie du circuit 1, sur la conduite d'évacuation 12, et qui en fonctionnement délivre un signal électrique de mesure S(t),
- des moyens électroniques 32 pour le traitement des signaux de mesure E(t) et S(t) délivrés par les capteurs 30 et 31. Les capteurs 30 et 31 sont identiques et peuvent d'une manière générale être constitués par tout capteur, qui permet de détecter et de caractériser l'écoulement ou les variations d'écoulement du fluide. Selon l'application, un signal de mesure E(t) ou S(t) caractérise la quantité de fluide par unité de temps passant au droit du capteur et/ou caractérise les ondes de matière engendrées par les variations d'écoulement du fluide.
Les signaux de mesure E(t) et S(t) délivrés par les capteurs 30 et 31 peuvent être de type analogique, ou être de type numérique lorsque les capteurs intègrent un convertisseur analogique/numérique.
Les capteurs 30 et 31 peuvent être de type invasif (par exemple sonde de mesure introduite à l'intérieur de la conduite correspondante 10 ou 12) ou de type non invasif (tel que dans l'exemple de capteurs illustré sur la figure 1). Les capteurs non invasifs sont toutefois préférés pour la mise en œuvre de l'invention, car les capteurs invasifs sont de manière préjudiciable la source d'une perte de charge sur le fluide.
De préférence, et de manière non limitative et non exhaustive de l'invention, les capteurs 30 et 31 sont des débitmètres choisis parmi la liste suivante :
- débitmètres électromagnétiques
- débitmètres par ultrasons
- débitmètres massiques (utilisant le principe de Coriolis) - débitmètres vortex
Les deux premiers types ci-dessus de débitmètres sont de type non invasif et sont donc préférentiels par rapport aux deux derniers types de débitmètres qui sont de type invasifs.
Plus particulièrement, pour la mise en oeuvre de l'invention, parmi ces débitmètres, on utilise de préférence les débitmètres électromagnétiques, car ils possèdent à ce jour les meilleures performances (précision, gamme de mesure, robustesse...).
Les moyens électroniques de traitement 32 peuvent être implémentés sous différente formes, sachant que l'important pour l'invention réside dans la méthode de traitement des signaux de mesure E(t) et S(t) qui est détaillée ci-après. Par exemple et de manière non limitative de l'invention, les moyens électroniques de traitement 32 peuvent être implémentés :
- sous la forme d'une unité de traitement programmable telle qu'un microordinateur exécutant un programme de traitement des signaux
E(t) et S(t) conforme à l'invention et chargé en mémoire vive, ou - sous la forme d'une carte électronique spécifique dont l'architecture électronique comporte un microprocesseur ou microcontrôleur apte à exécuter un programme embarqué de traitement des signaux E(t) et S(t) conforme à l'invention, ou - sous la forme d'une carte électronique comportant un circuit électronique spécifique de type ASIC, spécialement conçu pour exécuter un programme de traitement des signaux E(t) et S(t) conforme à l'invention.
On a représenté sur la figure 5, les principales étapes d'une variante préférée de mise en oeuvre de l'invention pour le traitement des signaux de mesure E(t) et S(t) délivrés par les capteurs 30 et 31.
D'une manière générale, le traitement des signaux repose sur une intercorrélation des signaux de mesure E(t) et S(t) [bloc 43] et sur une détection [bloc 44] et un suivi [bloc 45] du pic central (P) résultant de l'intercorrélation des signaux de mesure.
La fonction d'intercorrélation est une fonction mathématique qui est connue en soi. Dans le domaine temporel continu, cette fonction d'intercorrélation lcc(t), lorsqu'elle est appliquée aux signaux de mesure E(t) et S(t), est définie par la relation :
(1 ) Icc (t) = E(t) ® S(t) = + τ)dτ
L'intercorrélation consiste en définitive à calculer l'intégrale de recouvrement entre les signaux E(t) et S(t) en décalant temporellement un signal par rapport à l'autre.
En pratique, l'intercorrélation qui est calculée dans l'algorithme de la figure 5 [bloc 43] est une fonction d'intercorrélation discrétisée.
La fonction d'intercorrélation discrète lCq mise en œuvre est de préférence définie par la relation suivante :
(2) IcCj = ∑EkSJ+k k=û où : Ek sont les valeurs résultant de l'échantillonnage du signal de mesure E(t) [vecteur de taille n] ;
Sk sont les valeurs résultant de l'échantillonnage du signal de mesure S(t) [vecteur de taille n ]; n est le nombre d'échantillons total Ek et Sk dans une fenêtre de calcul prédéfinie ; k est un nombre entier compris entre 0 et (n-1) ; j est un nombre entier prenant les valeurs entre -(n-1) et (n-1) et prend ainsi les valeurs successives suivantes : -(n-1) ; -(n-2) ; ... ; -2 ; -1 ; 0 ; 1 ; 2 ; ... ; (n-1) ;
Iccj est un vecteur de taille 2n-1.
L'algorithme de la figure 5 va à présent être détaillé. Blocs (40) : Conversion Analogique/numérique
Pour permettre le calcul d'une fonction d'intercorrélation discrétisée, on réalise un échantillonnage des signaux de mesure continus E(t) et S(t), avec une fréquence d'échantillonnage (fe) prédéfinie.
Ces opérations d'échantillonnage sont de manière usuelle réalisées au moyen de convertisseurs analogiques/numériques. Ces convertisseurs peuvent être intégrés aux capteurs 30,31 ou être intégrés dans les moyens électroniques de traitement 32.
En sortie de cette étape d'échantillonnage, on obtient, au rythme de Ia fréquence d'échantillonnage (fe), une succession d'échantillons discrets Ek et Sk.
Blocs (41) et (42) : Correction offset Lorsque le régime du système est établi, les signaux de mesure de débit d'entrée E(t) et de sortie S(t) oscillent autour d'une valeur moyenne.
Cette valeur moyenne (Emoy et Smoy) est calculée en parallèle pour chaque signal [bloc 41].
La correction offset [bloc 42] consiste à retrancher à chaque échantillon Ek la valeur moyenne Emoy du signal correspondant, et à chaque échantillon Sk la valeur moyenne Smoy du signal correspondant. On obtient en sortie les signaux corrigés : (Ek - Emoy) et (Sk - Smoy)
Si le système ne dérive pas et/ou n'est pas fortement perturbé, ces valeurs moyennes Emoy et Smoy sont calculées une fois pour toutes et stockées en mémoire. En revanche, si le système dérive et/ou est fortement perturbé, ces valeurs moyennes Emoy et Smoy sont recalculées au fur et à mesure dans une fenêtre temporelle glissante de durée Tmoy prédéfinie. La valeur de cette durée Tm0y doit être supérieure à la durée d'apparition d'une fuite dans le circuit, et doit être suffisamment faible pour prendre en compte la dérive et/ou les perturbations du système. Cette durée Tmoy est fixée au cas par cas par l'homme du métier. De préférence, cette durée est choisie de telle sorte qu'en pratique les signaux corrigés (Ek - Emoy) et (Sk - Smoy) oscillent sensiblement autour de zéro.
Cette étape de correction d'offset est importante car elle permet, en combinaison avec l'intercorrélation des signaux corrigés (bloc 43) d'améliorer la sensibilité de la détection de fuite.
Dans la suite de la description, on considérera que les échantillons
Ek et Sk sont ceux issus directement de l'échantillonnage (variante de réalisation sans correction d'offset) ou sont les échantillons corrigés obtenus après soustraction de la valeur moyenne du signal (variante de réalisation avec correction offset).
Dans une autre variante de réalisation, cette correction offset pourrait être réalisée sur les signaux de mesure analogiques E(t) et S(t), (avant échantillonnage) au moyen d'un soustracteur analogique. Blocs 43 et 44 : Intercorrélation/Pic central Pour calculer l'intercorrélation entre les deux signaux discrétisés Ek et Sk, on stocke en mémoire à chaque calcul, un nombre n prédéfini d'échantillons successifs, et on calcule la fonction d'intercorrélation ICCJ discrétisée au moyen de la formule (2) précédemment décrite, c'est-à-dire (2n-1) valeurs successives lCc-(n-i) à lcc(n-i). Par exemple, et de manière non limitative de l'invention, on effectue le calcul de l'intercorrélation avec 1000 échantillons successifs Ek et Sk (k variant de 1 à n et n valant 1000), avec une fréquence d'échantillonnage fe de 1 kHz. Pour chaque série de n échantillons, on calcule ainsi un vecteur d'intercorrélation Icq toutes les secondes. Dans cet exemple, on calcule ainsi à chaque seconde 1999 valeurs lcc-(n-i)à lcc(n-i): lcc-(n-i) = Eo x Sn-I lθC-(n-2) = Eo x Sn-2 + El * Sn-I
So + En-I * Si Pour chaque série de (2n-1) valeurs d'échantillons Icq, on détermine automatiquement et de manière connue en soi le pic central p de l'intercorrélation (pic de plus grande amplitude).
Les étapes précitées d'intercorrélation et de détection du pic central, qui correspondent respectivement aux blocs 43 et 44 de la figure 5, sont effectuées de manière répétitive dans des fenêtres successives glissantes de n échantillons Ek et Sk. Ces fenêtres (ou séries) successives de n échantillons Ek ou Sk peuvent être sans chevauchement dans Ie temps (on prend n échantillons, puis les n échantillons suivants sans recoupement entre les séries d'échantillons), ou peuvent au contraire se chevaucher en partie.
Egalement, pour chaque calcul d'intercorrélation Icq, les deux fenêtres de calcul de n échantillons Ek et Sk qui sont utilisées peuvent être définies sur le même intervalle de temps, sans décalage temporel entre les fenêtres (dans ce cas les échantillons Ek et Sk successifs de chaque fenêtre ont tous été échantillonnés aux même instants), ou au contraire peuvent être définies avec un décalage temporel plus ou moins important entre les fenêtres de calcul.
La fréquence d'échantillonnage (fe) et Ie nombre n d'échantillons caractérisent un paramètre T, correspondant à la période d'intégration pour l'intercorrélation discrétisée. Ce paramètre T est défini par la relation suivante : T = — fe Dans l'exemple précité, avec une fréquence d'échantillonnage de 1KHz, et 1000 échantillons dans chaque fenêtre successive, la période d'intégration T est de 1 s.
Pour chaque application donnée, il revient à l'homme du métier de fixer judicieusement les valeurs de la fréquence d'échantillonnage (fe) et du nombre n d'échantillons, de telle sorte que la période d'intégration T soit compatible avec les temps caractéristiques du système. Bloc 45 : Suivi et alarme
On contrôle les variations du pic central P en comparant l'amplitude du pic central (P) de l'intercorrélation avec au moins un seuil (s).
Ce seuil (s) peut selon le cas être un seuil prédéfini de valeur constante.
Il peut également être constitué par un seuil auto-adaptif qui est fonction des valeurs Ek et Sk, et plus précisément des valeurs successives calculées pour le pic central de l'intercorrélation ; par exemple, le seuil auto- adaptif correspond à la valeur moyenne sur un nombre N prédéfini de valeurs de pic central précédemment calculées.
On détecte automatiquement une fuite, lorsque la variation de l'amplitude du pic central de l'intercorrélation est négative et est en valeur absolue supérieure à un seuil (s) prédéfini, et on déclenche le cas échéant une alarme (sonore, visuelle, envoi automatique d'un message d'alarme par tout moyen de télécommunication connu, etc..)
Il est également possible de calculer plusieurs seuils de détection (s) auto-adaptatifs, qui se différencient par la valeur de la variable N (nombre de valeurs de pic central utilisées pour calculer le seuil (s). Par exemple, un seuil (si) à court terme est calculé sur les 10 derniers pics centraux ( N= 10) ; un seuil (s2) à moyen terme est calculé sur les 100 derniers pics centraux (N=100) ; un seuil (s3) à long terme est calculé sur les 500 derniers pics centraux (N = 500). Lorsque la variation d'amplitude du pic central est négative et dépasse en valeur absolue l'un des seuils (si) à (s3), on déclenche une alarme caractéristique de ce seuil.
Résultats expérimentaux
Dans tous les exemples de mise en œuvre de l'invention décrits dans les tests expérimentaux ci-après, le fluide (f) circulant dans le circuit 1 est un liquide (en l'occurrence de l'eau) ; les capteurs 30 et 31 utilisés sont des débitmètres électromagnétiques ayant une précision de l'ordre de 0,5 %. La durée Tmoy pour la correction d'offset valait 30 secondes. La fréquence d'échantillonnage (fe) des signaux de mesure E(t) et S(t) est égale à 1 KHz. Le nombre d'échantillons (n) pour l'intercorrélation est égal à 1000. La fonction d'intercorrélation discrétisée Icq est calculée sur (n) échantillons Ek et Sk acquis simultanément et non décalés temporellement, et dans des fenêtres glissantes de n échantillons qui sont juxtaposées et ne se chevauchent pas. Détection d'une fuite On a représenté sur la figure 2, un circuit 1 de fluide modifié à des fins uniquement expérimentales pour provoquer des fuites de fluide.
Ce circuit 1 de la figure 2 se différencie de celui de la figure 1 uniquement en ce que sur une des tubulures 110 de ce circuit, on a réalisé une dérivation 13, sur laquelle est montée une vanne de fuite 14. Cette dérivation 13 débouche dans un réceptacle 15. Lorsque la vanne 14 est ouverte, une très faible partie du fluide (f) est prélevée du circuit 1 et n'est pas redirigée en sortie du circuit 1 vers l'ensemble 2, mais s'écoule dans le réceptacle 15, ce qui permet de provoquer une fuite dans le circuit 1 dans le cadre d'expérimentations de l'invention. Bien entendu, ces moyens 13, 14, 15 pour provoquer une fuite dans le circuit 1 sont prévus uniquement à des fins expérimentales pour tester l'invention, et ne se retrouvent pas dans la cadre d'une installation opérationnelle et non expérimentale finale comportant le circuit 1 et ses moyens 2 d'alimentation en fluide.
Dans un premier temps, on a fait fonctionner l'installation de la figure 2 sans fuite, en alimentant le circuit 1 avec de l'eau avec un débit moyen de l'ordre de 295 m3/h, la vanne de fuite 14 étant fermée. On a représenté sur les figures 6 et 7, sur une courte période de temps de fonctionnement sans fuite, les données brutes ER et Sk issues de l'échantillonnage des signaux de mesure E(t) et S(t) (signaux en sortie du bloc 42 après correction offset). On a représenté sur la figure 8, l' intercorrélation des signaux de mesure échantillonnés ER et Sk des figures 6 et 7.
Dans un deuxième temps, on provoque une fuite ayant un débit de l'ordre de 1,2m3/h, en ouvrant la vanne 14.
On a représenté sur les figures 9 et 10, sur une courte période de temps de fonctionnement prise en début de fuite, les données brutes Ek et Sk issues de l'échantillonnage des signaux de mesure E(t) et S(t) (signaux en sortie du bloc 42 après correction offset). On a représenté sur Ia figure 11, Pintercorrélation des signaux de mesure échantillonnés Ek et Sk des figures 9 et 10. La comparaison d'une part des figures 6 et 7 (signaux de mesure Ek et SR sans fuite) avec d'autre part les figures 9 et 10 (signaux de mesure Ek et Sk avec fuite) montre qu'il est très difficile de discriminer dans ces signaux un changement permettant de caractériser l'apparition d'une fuite. En revanche, la comparaison des figures 8 et 11 résultant de Pintercorrélation respectivement sans fuite et avec fuite, montre une modification notable du résultat de la fonction d'intercorrélation, et notamment une variation négative du minimum, qui permet de caractériser et détecter l'apparition de la fuite. Sensibilité de la détection de fuite
Pour caractériser la sensibilité de la détection de fuite, on a provoqué, en ouvrant de manière appropriée la vanne de fuite 14 de l'installation de la figure 2, cinq fuites successives présentant un débit croissant.
Ces fuites successives sont caractérisées dans le tableau I ci-après. Tableau I
On a représenté sur les figures 12 et 13, les signaux de mesure E(t) et S(t) délivrés par les capteurs 30 et 31 pour la période de temps entre 470s et 570s. On a représenté sur la figure 14, l'évolution dans le temps de l'amplitude des pics centraux de l'intercorrélation (sortie du bloc 44) pour la période de temps entre 0s et 730s. Sur cette figure 14, les parties de la courbe permettant de discriminer les cinq fuites successives sont indiquées par les flèches F1 , F2, F3, F4 et F5.
L'analyse de la figure 14, permet de montrer une variation négative brusque de l'amplitude du pic central de l'intercorrélation à chaque apparition d'une fuite. Elle permet également de montrer que cette variation d'amplitude, qui permet de caractériser et de détecter l'apparition d'une fuite, augmente en valeur absolue avec le débit de Ia fuite.
On a représenté sur la figue 15 la variation d'amplitude (en valeur absolue) du pic central d'intercorrélation caractéristique d'une fuite en fonction de la valeur relative de la fuite provoquée (en m3/h). Cette figure 15 montre que la variation de l'amplitude du pic central d'intercorrélation est proportionnelle à la valeur relative de la fuite. L'intercorrélation peut donc non seulement être utilisée pour détecter l'apparition d'une fuite dans le circuit 1, mais permet également avantageusement de caractériser la valeur relative de la fuite qui a été détectée.
Enfin, cette expérience a permis de valider qu'il était possible de détecter des fuites de très faible amplitude par rapport au débit principal du fluide à l'entrée du circuit 1. Dans le cas de la première fuite provoquée qui présente la plus faible amplitude, il a été possible de détecter cette fuite avec une sensibilité de l'ordre du millième (1.9.10"3), alors que comparativement la précision de mesure des débitmètres utilisés était de l'ordre de 0,5 %. Perturbations thermiques
On représenté sur la figure 3, un circuit 1 de fluide modifié à des fins uniquement expérimentales pour provoquer des perturbations thermiques dans le circuit 1. Ce circuit 1 de la figure 3 se différencie de celui de la figure 2 uniquement en ce qu'on a ajouté sur l'une des tubulures 110 du circuit 1 :
- un élément chauffant 16 (par exemple une résistance chauffante) qui permet de chauffer localement la tubulure 110 et par là-même la partie de fluide circulant localement dans cette tubulure 110 ;
- deux capteurs de température 17a et 17b permettant de mesurer les variations de température du fluide circulant dans ladite tubulure 110 équipée de l'élément chauffant 16.
Dans une première phase, on fait fonctionner l'installation de la figure 3, sans provoquer de fuite (vanne 14 fermée).
On a représenté sur les figures 16 et 17 respectivement les signaux de mesure E(t) et S(t) pendant cet essai avec chauffage du fluide et sans fuite, et sur la figure 18, la température du fluide en fonction du temps
(sensiblement constante) mesurée en amont de l'élément chauffant 16 par le capteur 17a (courbe A) et la température du fluide en fonction du temps mesurée en aval de l'élément chauffant 16 par le capteur 17b (courbe B). On a représenté sur la figure 19, l'évolution dans le temps de l'amplitude des pics centraux d'intercorrélation. Cette courbe de la figure 19 montre, qu'en début de chauffage, la variation de l'amplitude des pics centraux d'intercorrélation augmente brusquement, mais dans le sens opposé (variation positive) à la variation de l'amplitude des pics centraux d'intercorrélation en cas de fuite.
Dans une deuxième phase, on arrête le chauffage de la tubulure 110 pendant une durée suffisante pour revenir à des conditions de température du fluide sensiblement identiques dans tout le circuit 1.
Dans une troisième phase, on chauffe de nouveau la tubulure 110 au moyen de l'élément chauffant 16, puis on génère une fuite en ouvrant la vanne 14.
Les caractéristiques de la fuite sont données le tableau II ci- dessous : Tableau II
On a représenté sur Ia figure 20, la température du fluide en fonction du temps (sensiblement constante) mesurée en amont de l'élément chauffant 16 par le capteur 17a (courbe A) et la température du fluide en fonction du temps mesurée en aval de l'élément chauffant 16 par le capteur 17b (courbe B), en fin de deuxième phase et au cours de la troisième phase de chauffage.
On a représenté sur les figures 21 et 22 respectivement les signaux de mesure E(t) et S(t) pendant cet essai de chauffage avec fuite, et sur la figure 23 l'évolution dans le temps de l'amplitude des pics centraux d'intercorrélation. En référence à la courbe 23, on constate une brusque variation négative de l'amplitude des pics centraux de l'intercorrélation très peu de temps après le déclenchement de la fuite. Cette variation brusque de l'amplitude des pics centraux d'intercorrélation permet ainsi de toujours détecter l'apparition de la fuite, malgré la perturbation thermique locale dans le circuit 1 , qui a été provoquée par le chauffage de la tubulure 110 au moyen de l'élément chauffant 16. Perturbations hydrauliques
On a représenté sur la figure 4, un circuit 1 de fluide modifié à des fins uniquement expérimentales pour provoquer des perturbations hydrauliques dans le circuit 1. Ce circuit 1 de la figure 4 se différencie de celui de la figure 2 uniquement en ce qu'on a ajouté sur l'une des tubulures 110 du circuit 1 une vanne 18 permettant de fermer ou d'ouvrir ladite tubulure 110.
Pour générer une perturbation hydraulique, on ferme la vanne 18. On a représenté sur les figures 24 et 25 respectivement les signaux de mesure E(t) et S(t) pendant cet essai avec fermeture de la vanne 18, et sur la figure 26 l'évolution dans le temps de l'amplitude des pics centraux de l'intercorrélation.
Cette courbe de la figure 26 montre que très peu de temps après la fermeture de la vanne 18 (perturbation hydraulique dans le circuit 1) la variation de l'amplitude des pics centraux d' intercorrélation augmente brusquement, mais dans le sens opposé (variation positive) à la variation de l'amplitude des pics centraux d'intercorrélation en cas de fuite. Ce comportement est comparable à ce qui a été précédemment décrit pour les perturbations thermiques, en référence notamment à la figure 19. Sachant que l'effet de la perturbation hydraulique sur les résultats de l'intercorrélation (variation positive de l'amplitude des pics centraux d'intercorrélation) est inversé par rapport à l'effet d'une fuite (variation négative de l'amplitude des pics centraux d'intercorrélation), on peut en conclure que le suivi de l'évolution dans le temps de l'amplitude des pics centraux résultant de l'intercorrélation des signaux de mesure E(t) et S(t) permet toujours de détecter l'apparition d'une fuite dans le circuit, même en présence de perturbations hydrauliques. Autres variantes de réalisations
L'invention n'est pas limitée à la variante préférée de réalisation qui vient d'être décrite en référence aux figures 1 à 26. En particulier et de manière non exhaustive, la détection du pic central de corrélation (pic de plus grande amplitude) ( bloc 44 de la figure 5) peut être remplacée par un calcul de la somme de tous les échantillons lCCj issus de l'intercorrélation (ce qui revient à calculer l'intégrale de l'intercorrélation). En l'absence de fuite (voir figure 8), cette somme est proche de zéro. En présence d'une fuite (voir figure 11), le résultat de cette somme devient fortement négatif. Pour la détection de fuite (bloc 45), on choisit dans ce cas un seuil (s) négatif. Lorsque le résultat négatif de cette somme des échantillons lCCj issus de l'intercorrélation devient inférieur à ce seuil (s) négatif prédéfini, dans ce cas on détecte automatiquement l'apparition d'une fuite. Dans une autre variante de réalisation, illustrée sur la figure 27, on réalise la détection (bloc 45) non plus directement sur Ie résultat de l'intercorrélation, mais sur le résultat de l'intercorrélation diminué (bloc 47) du résultat de l'autocorrélation (bloc 46) du signal de sortie corrigé (Sk-Smoy). Dans un autre variante, on peut également réaliser la détection (bloc 45) sur le résultat de l'intercorrélation diminué (bloc 47) du résultat de l'autocorrélation (bloc 46) du signal d'entrée corrigé (Ek-Em0y). Ceci permet avantageusement d'améliorer le rapport signal sur bruit et d'améliorer la sensibilité de détection de fuite. Applications de l'invention L'invention trouve son application à la détection de l'apparition de fuite(s) de fluide dans tout circuit à l'intérieur duquel circule un fluide, ledit fluide pouvant être un liquide ou un gaz, ou un mélange gaz/liquide.
De préférence, l'invention trouve son intérêt dans toutes les applications: - où le débit des fuites à détecter peut être très faible rapport au débit de fluide à l'entrée du circuit, et/ou
- le fluide subit des perturbations modifiant son volume (notamment perturbations thermiques, hydrauliques, ...). Un exemple particulièrement intéressant d'application de l'invention, réside dans la détection de fuites dans un circuit de refroidissement ou dans un circuit de chauffage, à l'intérieur duquel circule un fluide caloporteur.
Plus particulièrement, l'invention peut avantageusement être utilisée dans le domaine de la sidérurgie pour détecter automatiquement de manière fiable et rapide l'apparition de fuite(s) dans les circuits de refroidissement d'un four électrique à arcs qui forment un réseau complexe de tubulures. Dans ce cas, le circuit 1 de la figure 1 est par exemple le circuit de refroidissement équipant la voûte d'un four électrique à arcs, ou le circuit de refroidissement des panneaux d'un four électrique à arcs.
Un exemple particulièrement intéressant d'application de l'invention réside dans la détection de fuites dans une canalisation ou dans un réseau de canalisations pour le transport et/ou de distribution de fluides (pipelines, gazoducs, etc..)

Claims

REVENDICATIONS
1. Procédé de contrôle d'un circuit (1) dans lequel circule un fluide (f), ledit procédé comportant les étapes suivantes : (a) on détecte le fluide en circulation au moyen d'un capteur d'entrée (30) et d'un capteur de sortie (31) positionnés respectivement à l'entrée et à la sortie du circuit, et délivrant respectivement deux signaux de mesure E(t) et S(t) caractéristiques de l'écoulement ou des variations d'écoulement du fluide,
(b) on traite les deux signaux de mesure E(t) , S(t) dans le but de détecter le cas échéant une fuite de fluide dans le circuit, caractérisé en ce qu'à l'étape (b) de traitement des signaux E(t) et S(t), on effectue une correction offset en retirant aux signaux de mesure [E(t) ou ER ; S(t) ou Sk] leur valeur moyenne (Emoy ; Smoy), et on calcule une intercorrélation des deux signaux de mesure corrigés.
2. Procédé selon la revendication 1 , caractérisé en ce qu'à l'étape (b) de traitement des signaux E(t) et S(t), on réalise un échantillonnage des signaux de mesure E(t) et S(t), et la fonction d'intercorrélation
Icc est une fonction discrète calculée au moyen de la formule suivante :
IcCj = ∑EkSJ+k
où : Ek sont les valeurs résultant de l'échantillonnage du signal de mesure E(t);
Sk sont les valeurs résultant de l'échantillonnage du signal de mesure S(t); n est le nombre d'échantillons total Ek et Sk dans une fenêtre de calcul prédéfinie ; k est un nombre entier compris entre O et (n-1) ; j est un nombre entier prenant les valeurs entre ~(n-1) et (n-
1) ;
Iccj est un vecteur de taille 2n-1.
3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce qu'à l'étape (b) de traitement des signaux de mesure, on calcule une intercorrélation des deux signaux de mesure, diminuée de l'autocorrélation du signal de mesure de sortie ou d'entrée.
4. Procédé selon l'une des revendications 1 à 3, caractérisé en ce qu'à l'étape (b) de traitement des signaux de mesure, on détecte les variations d'amplitude du pic central (P) d'intercorrélation.
5. Procédé selon la revendication 4, caractérisé en ce qu'à l'étape (b), de traitement des signaux de mesure, on compare les variations de l'amplitude du pic central d'intercorrélation (P) avec au moins un seuil prédéfini (s), et on détecte l'apparition d'une fuite lorsque ladite variation de l'amplitude est négative et devient, en valeur absolue, supérieure à ce seuil.
6. Procédé selon la revendication 5, caractérisé en ce que ledit seuil (si , s2 ou s3) est auto-adaptatif et est calculé en fonction des amplitudes des pics centraux d'intercorrélation précédents.
7. Procédé selon l'une des revendications 1 à 3, caractérisé en ce qu'à l'étape (b) de traitement des signaux de mesure, on calcule l'intégrale de l'intercorrélation, et on compare le résultat de cette intégrale avec un seuil (s) prédéfini.
8. Utilisation du procédé visé à l'une des revendications 1 à 7, pour détecter l'apparition d'une fuite dans un circuit de refroidissement ou dans un circuit de chauffage.
9. Utilisation selon la revendication 8, pour détecter l'apparition d'une fuite dans un circuit de refroidissement d'un four électrique à arcs.
10. Utilisation du procédé visé à l'une des revendications 1 à 7, pour détecter l'apparition d'une fuite dans une canalisation ou un réseau de canalisations de transport et/ou de distribution d'au moins un fluide.
11. Installation comportant au moins un circuit (1), des moyens (2) de mise en circulation d'un fluide (f) dans ce circuit (1), un capteur d'entrée (30) et un capteur de sortie (31) positionnés respectivement à l'entrée et à la sortie du circuit (1) , et délivrant respectivement deux signaux de mesure E(t) et S(t) caractéristiques de l'écoulement ou des variations d'écoulement du fluide, et des moyens électroniques (32) de traitement des signaux de mesure E(t) et S(t), caractérisée en ce que lesdits moyens électroniques (32) sont conçus pour mettre en œuvre l'étape (b) de traitement de signaux de mesure E(t) et (St) qui est définie dans le procédé visé à l'une des revendications 1 à 7.
12. Installation selon la revendication 11 , caractérisée en ce que les capteurs (30,31) sont des débitmètres de type non invasif.
13. Installation selon la revendication 11 ou 12, caractérisée en ce que le circuit (1) est un circuit de refroidissement ou un circuit de chauffage.
14. Installation selon la revendication 13, caractérisée en ce qu'elle est constituée par un four électrique à arcs.
15. Installation selon la revendication 11 ou 12, caractérisée en ce que le circuit (1) est une canalisation ou un réseau de canalisations pour le transport et/ou la distribution d'au moins un fluide.
16. Programme enregistré sur un support ou dans une mémoire, et qui, lorsqu'il est exécuté par une unité de traitement programmable (32), réalise automatiquement le procédé de contrôle visé à l'une des revendications 1 à 7.
EP06847154A 2005-12-28 2006-12-28 Detection d'une fuite de fluide dans un circuit par intercorrelation Withdrawn EP1966580A2 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0513405A FR2895508B1 (fr) 2005-12-28 2005-12-28 Detection d'une fuite de fluide dans un circuit par intercorrelation
PCT/FR2006/002896 WO2007074244A2 (fr) 2005-12-28 2006-12-28 Detection d'une fuite de fluide dans un circuit par intercorrelation

Publications (1)

Publication Number Publication Date
EP1966580A2 true EP1966580A2 (fr) 2008-09-10

Family

ID=37451189

Family Applications (1)

Application Number Title Priority Date Filing Date
EP06847154A Withdrawn EP1966580A2 (fr) 2005-12-28 2006-12-28 Detection d'une fuite de fluide dans un circuit par intercorrelation

Country Status (3)

Country Link
EP (1) EP1966580A2 (fr)
FR (1) FR2895508B1 (fr)
WO (1) WO2007074244A2 (fr)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8850871B2 (en) 2010-09-30 2014-10-07 Siemens Aktiengesellschaft Pipeline leak location using ultrasonic flowmeters
CN102305552B (zh) * 2011-08-16 2013-11-13 东北大学 钢包底吹喷粉漏钢检测装置及漏钢检测方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2603715C3 (de) * 1976-01-31 1980-03-27 Rolf Prof. Dr.-Ing. 6100 Darmstadt Isermann Verfahren zur Leckerkennung und Leckortung bei Rohrleitungen
DE3607913C1 (en) * 1986-03-10 1987-08-13 Fraunhofer Ges Forschung Method for locating leaks in pipelines
DE69305190T2 (de) * 1992-01-16 1997-03-13 Toshiba Kawasaki Kk Verfahren und Vorrichtung zur Positionserkennung einer anormalen Stelle einer erdverlegten Rohrleitung
US5343737A (en) * 1992-09-22 1994-09-06 Joseph Baumoel Method and apparatus for leak detection and pipeline temperature modelling method and apparatus
FR2727205A1 (fr) * 1994-11-22 1996-05-24 Gaz De France Procede et dispositif de detection et de localisation de fuites dans une conduite parcourue par un fluide
US5565764A (en) * 1995-05-05 1996-10-15 Texas Instruments Incorporated Digital processing method for parameter estimation of synchronous, asynchronous, coherent or non-coherent signals
JP3857670B2 (ja) * 2003-07-04 2006-12-13 株式会社山武 不明水発生分布推定装置、方法およびプログラム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO2007074244A3 *

Also Published As

Publication number Publication date
FR2895508A1 (fr) 2007-06-29
WO2007074244A2 (fr) 2007-07-05
FR2895508B1 (fr) 2008-03-07
WO2007074244A3 (fr) 2007-08-23

Similar Documents

Publication Publication Date Title
EP2693176B1 (fr) Procédé de détection de défauts d&#39;un roulement par analyse vibratoire
EP1799869B1 (fr) Procede de contrôle du bullage d&#39;un bain de metal en fusion dans un recipient metallurgique et installation de mise en oeuvre
EP0184521B1 (fr) Procédé et dispositif de détection des fuites dans un échangeur de chaleur en fonctionnement
EP3325934B1 (fr) Procédé de contrôle de l&#39;étanchéité de produits scellés et installation de détection de fuites
FR2647549A1 (fr) Procede et dispositif pour mesurer des qualites d&#39;un fluide polyphasique
FR2805042A1 (fr) Procede et dispositif non intrusif pour caracteriser les perturbations d&#39;ecoulement d&#39;un fluide a l&#39;interieur d&#39;une canalisation
EP3658763B1 (fr) Procédé et dispositif de détection d&#39;allumage d&#39;une chambre à combustion de moteur-fusée, procédé de démarrage d&#39;un moteur-fusée, programme d&#39;ordinateur, support d&#39;enregistrement et moteur-fusée
EP2936145A1 (fr) Dispositif et procede de discrimination d&#39;un gaz dans un echantillon
WO2007074244A2 (fr) Detection d&#39;une fuite de fluide dans un circuit par intercorrelation
EP3023766A1 (fr) Procédé et dispositif de mesure de perméation par spectrométrie de masse
CN115510568A (zh) 管道故障检测方法及系统、电子设备和存储介质
WO2008081089A2 (fr) Procede et dispositif de detection et/ou de quantification de fuites d&#39;eau
FR2895509A1 (fr) Detection d&#39;une fuite de fluide dans un circuit par intercorrelation
WO2000043762A1 (fr) Procede et dispositif de detection de la formation d&#39;un depot de matiere sur une face d&#39;un capteur de flux thermique
EP0080406A1 (fr) Procédé d&#39;interrogation d&#39;un détecteur de teneur en gaz combustible et dispositif pour sa mise en oeuvre
FR2946704A1 (fr) Methode pour reduire la perte de charge d&#39;un liquide en ecoulement dans une conduite en tenant compte de la degradation d&#39;agents reducteurs de trainee.
WO2016042060A1 (fr) Systeme et procede de detection de particules
EP4065959B1 (fr) Procede de caracterisation de composes d&#39;intérêt dans une chambre de mesure présentant une variation d&#39;humidité relative
EP3901587B1 (fr) Procédé de mesure ultrasonique avec prise en compte de la quantité de bulles gazeuses
EP3593147A1 (fr) Procede de mesure d&#39;une vitesse d&#39;un fluide
EP3377862B1 (fr) Procédé de détection du niveau d&#39;un fluide dans un état liquide et/ou solide dans un réservoir et système associé
EP3942375B1 (fr) Méthode de détection d&#39;anomalies dans une installation de traitement des eaux mettant en oeuvre un appareil d&#39;injection d&#39;oxygène dans un bassin d&#39;épuration
WO2009150353A2 (fr) Procede et installation de determination d&#39;une valeur de debit d&#39;un fluide
JP2004333280A (ja) 濃度測定装置
FR3075970A1 (fr) Dispositif de verification de l’integrite d’un systeme de detection de court-circuit

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20080701

AK Designated contracting states

Kind code of ref document: A2

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LI LT LU LV MC NL PL PT RO SE SI SK TR

17Q First examination report despatched

Effective date: 20090226

DAX Request for extension of the european patent (deleted)
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20120703