FR2895508A1 - Circuit e.g. furnace roof/panel cooling circuit, controlling method for e.g. electric arc furnace, involves processing signals of fluid flow/flow variation by calculating intercorrelation function of signals to detect fluid leakage - Google Patents
Circuit e.g. furnace roof/panel cooling circuit, controlling method for e.g. electric arc furnace, involves processing signals of fluid flow/flow variation by calculating intercorrelation function of signals to detect fluid leakage Download PDFInfo
- Publication number
- FR2895508A1 FR2895508A1 FR0513405A FR0513405A FR2895508A1 FR 2895508 A1 FR2895508 A1 FR 2895508A1 FR 0513405 A FR0513405 A FR 0513405A FR 0513405 A FR0513405 A FR 0513405A FR 2895508 A1 FR2895508 A1 FR 2895508A1
- Authority
- FR
- France
- Prior art keywords
- circuit
- fluid
- signals
- leak
- measurement signals
- 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.)
- Granted
Links
- 239000012530 fluid Substances 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000012545 processing Methods 0.000 title claims description 30
- 238000001816 cooling Methods 0.000 title claims description 24
- 238000010891 electric arc Methods 0.000 title claims description 8
- 238000005259 measurement Methods 0.000 claims abstract description 65
- 230000006870 function Effects 0.000 claims abstract description 19
- 238000009434 installation Methods 0.000 claims abstract description 17
- 238000010438 heat treatment Methods 0.000 claims description 22
- 238000005070 sampling Methods 0.000 claims description 19
- 238000012937 correction Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000001514 detection method Methods 0.000 description 29
- 239000007789 gas Substances 0.000 description 11
- 239000007788 liquid Substances 0.000 description 11
- 229910000831 Steel Inorganic materials 0.000 description 8
- 239000010959 steel Substances 0.000 description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 230000010354 integration Effects 0.000 description 7
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 6
- 230000035945 sensitivity Effects 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- 229910052742 iron Inorganic materials 0.000 description 3
- 238000002844 melting Methods 0.000 description 3
- 230000008018 melting Effects 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 239000011449 brick Substances 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000002826 coolant Substances 0.000 description 2
- 238000005314 correlation function Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- 229910052751 metal Inorganic materials 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 238000007670 refining Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000001960 triggered effect Effects 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- CWYNVVGOOAEACU-UHFFFAOYSA-N Fe2+ Chemical compound [Fe+2] CWYNVVGOOAEACU-UHFFFAOYSA-N 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000009529 body temperature measurement Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000013529 heat transfer fluid Substances 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007620 mathematical function Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M3/00—Investigating fluid-tightness of structures
- G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
- G01M3/26—Investigating 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/28—Investigating 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/2807—Investigating 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M3/00—Investigating fluid-tightness of structures
- G01M3/02—Investigating fluid-tightness of structures by using fluid or vacuum
- G01M3/04—Investigating fluid-tightness of structures by using fluid or vacuum by detecting the presence of fluid at the leakage point
- G01M3/24—Investigating 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/243—Investigating 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
Description
DETECTION D'UNE FUITE DE FLUIDE DANS UN CIRCUIT PAR INTERCORRELATIONDETECTION OF A FLUID LEAK IN A CIRCUIT BY 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. Technical Field 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. For example, 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. 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 transport and / or fluid distribution pipes (pipelines, pipelines, etc.). Prior art A method widely used to date for detecting a fluid leak in a circuit, such as for example a circuit for cooling or heating an enclosure, consists in placing two flowmeters respectively at the output and at the input of the circuit, and in 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. This method, hereinafter referred to as the flow differential method, has at least three major disadvantages.
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 25 perturber les signaux de mesure de débit, et de ce fait fausser la détection de fuite (3ème 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 30 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. It is limited by the precision of the measurement of the flowmeters used (1st disadvantage). The lower the accuracy of the flow meters used, the greater the flow rate of the fluid leak must be so that it can be detected. For example, and as an indication, the currently marketed electromagnetic flowmeters have an accuracy of the order of 0.3% to 1%. Thus, with this type of flow meter, the differential flow rate method can best be used today to detect leak rates that are at least greater, depending on the case, from 0.6% to 2% of the flow rate of the fluid. in the circuit. Also, 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. Conversely, in the case of a reduction in the volume of fluid transported, for example under the effect of a sharp drop in the temperature of the fluid between its inlet and its outlet of the circuit, the output flow rate becomes lower than the input flow, which can trigger a false leak detection, while no leak is present. The presence of electromagnetic disturbances can also disturb the flow measurement signals, and thus distort the leak detection (3rd disadvantage). In general, the three aforementioned drawbacks make this flow differential method unreliable and, if necessary, unsuitable for detecting fluid leaks having a very low flow rate and / or leakage of a fluid subjected to disturbances. in particular thermal, hydraulic, affecting the volume of fluid transported, or electromagnetic. However, there are numerous industrial applications in which a fluid is circulating in a circuit, and it is necessary to reliably and sufficiently early detect any leak, even very low, and moreover in an environment in which the fluid can to undergo non-stationary disturbances, and in particular more or less significant thermal disturbances over time. Among these many industrial applications, we can particularly mention all industrial applications in which a fluid flows in a cooling circuit or in a heating circuit of an installation. More particularly, in the field of iron and steel, it is usual to cool 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.
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 1500 C 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. In this type of furnace, the furnace tank is loaded 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. Typically, 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.
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. This type of steel furnace is particularly vulnerable to liquid leaks in its cooling circuit, for at least two reasons: if because of a leak in the cooling circuit, water flows inside the furnace, it can reach the bottom of the oven which is lined with refractory bricks; in this case, the refractory bricks deteriorate very rapidly, and this deterioration results in a piercing of the bottom of the furnace tank; in this case, the molten metal escapes out of the tank with potentially serious human and financial consequences; - if, because of a leak in the cooling circuit, water falls inside the tank of a furnace in operation, and comes into contact with the molten metal, it can decompose into oxygen and hydrogen under certain conditions because of the very high temperature. This decomposition generates a significant risk of accidental explosion.
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. It is therefore understandable that in this particular type of application the stakes for finding an automatic leak detection that is reliable and early are extremely important, both on the human level to avoid accidents that can be fatal in certain cases, and on the plane. to avoid long and costly stopping of the furnace in the event of an accident. Today, the leak detection methods that are implemented industrially on these steel furnaces are essentially methods based on a measurement of the differential of the flow rates of entry and exit, and therefore do not give entirely satisfactory considering the disadvantages mentioned above inherent to this type of detection method. To try to correct the effects of thermal disturbances, the volume of fluid transported, and try to overcome the aforementioned problems of false detection or missed detection that may result, improvements of the aforementioned method of differential flow rates have been proposed to this day. These improvements are based on a correction of the flowmeter measurements from, in particular, the temperature or the pressure of the fluid flow. Temperature correction solutions are proposed by the manufacturers of flow meters, and are implemented in particular in the field of arc furnaces mentioned above.
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é à la faible précision de mesure des débitmètres. For example also, in another field of application described in the French patent application FR 2,509,839, there is provided a method for detecting nitrogen leakage in a circuit forming a barrier between the primary and secondary spaces for thermal insulation. a cryogenic tank, said method being based on the differential nitrogen flow rate method at the outlet and the inlet of the circuit. In this publication, it is taught to provide temperature corrections and pressure corrections to flow measurements by measuring the temperature and pressure of the flowing fluid. On the one hand this type of corrections is ultimately not very reliable and also imposes a calibration of temperature and pressure measurements that can be difficult to achieve, and must be performed regularly. On the other hand and above all, these corrections of temperature or pressure do not make it possible to overcome the first disadvantage mentioned above related to the low measurement accuracy of the flow meters.
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 = E f Qdt = Q. dt - jeu, dt Q,n étant la quantité de gaz qui entre dans le pipeline et eu, étant la 30 quantité de gaz qui sort du pipeline. It is furthermore proposed in the European patent application EP 0 188 911 a detection method which is similar to the aforementioned method of differential flow, and which in this document is used to detect a fluid leak, and in particular a leak of gas, in a pipeline type pipeline. This method is based essentially on the calculation of the FD parameter below: FD = E f Qdt = Q. dt - set, dt Q, n being the quantity of gas entering the pipeline and eu, being the quantity of gas coming out of the 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 : FD* _ f Q,n dt - fQou, dt JQ,ndt 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. 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 par rapport à la méthode précitée du différentiel des débits. In 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: FD * _ f Q, n dt - fQou, dt JQ, ndt The calculation of the parameter FD or the parameter FD * has all the aforementioned drawbacks of the differential flow rate method. Moreover, in this publication, it is taught to detect a leak when the calculated parameter (FD or FD *) exceeds a threshold value calculated by simulation. This simulation takes into account many parameters including the diameter of the pipe in which the fluid circulates, and any changes of section, the pressure of the fluid flow, fluid friction, gravity, ... 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. 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. OBJECTS OF THE INVENTION The general objective of the present invention is to propose a new technical solution that makes it possible to automatically detect a 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 a leakage of a fluid (liquid and / or gas) flowing in a circuit, with an improved detection sensitivity compared to the aforementioned method of differential flow.
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. 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.
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. 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. .
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. Conformément à l'invention, à l'étape (b) de traitement des 30 signaux E(t) et S(t), on calcule une fonction d'intercorrélation Icc des deux signaux de mesure E(t) et S(t). SUMMARY OF THE INVENTION All or some of these objectives, including at least the general objective mentioned above, are achieved by the invention, which has for its first object 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, (b) the two measurement signals E (t), S (t) are processed for the purpose of detecting, if necessary, a leakage of fluid in the circuit. According to the invention, in step (b) of processing the signals E (t) and S (t), an intercorrelation function Icc of the two measurement signals E (t) and S (t) is calculated. .
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 : More particularly, the method comprises the following additional and optional features, taken alone or, where appropriate, in combination:
- à 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 : n-1 Icci = EkS/+k k=0 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: n-1 Icci = EkS / + kk = 0
où :or :
Ek sont les valeurs résultant de l'échantillonnage du signal de mesure E(t); Ek are the values resulting from the sampling of the measurement signal E (t);
Sk sont les valeurs résultant de l'échantillonnage du signal de mesure S(t); Sk are the values resulting from the sampling of the measurement signal S (t);
n est le nombre d'échantillons total Ek et Sk dans une fenêtre de calcul prédéfinie ; n is the total number of samples Ek and Sk in a predefined calculation window;
k est un nombre entier compris entre 0 et (n-1) ; k is an integer from 0 to (n-1);
j est un nombre entier prenant les valeurs entre -(n-1) et (n-1) ; j is an integer taking the values between - (n-1) and (n-1);
Iccj est un vecteur de taille 2n-1 ; Iccj is a vector of size 2n-1;
20 à l'étape (b) de traitement des signaux de mesure, on réalise une correction offset avant calcul de I'intercorrélation, en retirant au signaux de mesure [E(t) ou Ek ; S(t) ou Sk] leur valeur moyenne (Emoy ; Smoy) ; In the measuring signal processing step (b), an offset correction is performed before calculating the inter-correlation, by removing from the measurement signals [E (t) or Ek; S (t) or Sk] their mean value (Emoy, Smoy);
à l'étape (b) de traitement des signaux de mesure, on détecte les 25 pics centraux (P) d'intercorrélation ; in step (b) of processing the measurement signals, the central peaks (P) of intercorrelation are detected;
- à 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, et on détecte l'apparition d'une fuite lorsque ladite variation de l'amplitude est négative et devient, en 30 valeur absolue, supérieure à ce seuil ; 15 - au moins un seuil est auto-adaptatif et est calculé en fonction des amplitudes des pics centraux d'intercorrélation précédents. 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 oeuvre 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. Dns 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. in step (b), for processing the measurement signals, 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; At least one threshold is auto-adaptive and is calculated according to the amplitudes of the previous central cross-correlation peaks. 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 respectively delivering 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. In a particular variant embodiment, the sensors are preferably non-invasive type flow meters. In a particular application of the invention, the circuit is a cooling circuit or a heating circuit. More particularly, the installation is for example constituted by an electric arc furnace. In another application, the circuit is a pipe or pipe network for transporting and / or dispensing at least one fluid.
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 15 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, 20 - 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 25 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. Description détaillée On a représenté de manière schématique sur la figure 1 un circuit 30 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. 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. BRIEF DESCRIPTION OF THE DRAWINGS Other characteristics and advantages of the invention will appear more clearly on reading the detailed description below of several preferred embodiments of the invention, which description is given by way of non-limiting example and Non-exhaustive of the invention, and with reference to the accompanying drawings in which: Figure 1 schematically shows an installation of the invention, Figure 2 schematically shows an installation of the invention, which has been modified to Fig. 3 shows schematically an installation of the invention, which has been modified for experimental purposes to cause thermal disturbances, Fig. 4 schematically shows an installation of the invention. invention, which has been modified for experimental purposes to cause hydraulic disturbances FIG. 5 is a block diagram illustrating the main successive steps of a preferred embodiment of step (b) of processing the measurement signals according to the invention. FIGS. 6 to 26 are curves resulting from FIG. tests performed using the facilities of Figures 2 to 4, and are detailed below. DETAILED DESCRIPTION FIG. 1 diagrammatically shows a circuit 30 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. fluid, associated with a forced circulation pump of this fluid. The circuit 1 comprises: a main inlet duct 10, the inlet of which is connected in leaktight manner to the outlet of the assembly 2 (reservoir / pump), a complex network 11 of pipes 110 in parallel, of which common inlet to all the pipes 110, is sealingly connected to the outlet of the main inlet pipe 10, said pipes 110 may have different geometrical characteristics (lengths, diameters, etc ...), main discharge 12, the inlet of which is connected in leaktight manner to the outlet of the network 11 common to all of the pipes 110, and whose outlet is sealingly connected to the inlet of the circulation assembly 2 forced fluid (f). When the pump of the assembly 2 operates, 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 discharge pipe 12. In the particular example illustrated, the circuit 1 comprises a complex network 11 of pipes 110 defining for the fluid (f) several possible paths between the inlet and the circuit output 1. Although 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 circuit of more structure single path single and constituted for example by a single 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. In order to automatically detect a fluid leak (f) in the circuit 1, 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 input of the circuit 1, on the main intake pipe 10, and which in operation delivers an electrical measurement signal E (t), a sensor 31 which is mounted at the outlet of the circuit 1, on the discharge pipe 12, and which in operation delivers an electrical measurement signal S (t), - electronic means 32 for the processing of the measurement signals E (t) and S (t ) delivered by the sensors 30 and 31. The sensors 30 and 31 are identical and can in general be constituted by any sensor, which makes it possible to detect and characterize the flow or the fluid flow variations. Depending on the application, 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 may be of analog type, or be of digital type when the sensors integrate an analog / digital converter.
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 oeuvre 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 30 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 30 d'intercorrélation Icc(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) = JE(r)S(t +z)di 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. 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. Preferably, and in a nonlimiting and non-exhaustive manner of the invention, the sensors 30 and 31 are flow meters selected from the following list: electromagnetic flowmeters ultrasonic flowmeters mass flowmeters (using the Coriolis principle) vortex flowmeters The first two types The above flowmeters are of the non-invasive type and are therefore preferred over the last two types of flow meters which are of the invasive type. More particularly, for the implementation of the invention, among these flowmeters, 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. For example and in a nonlimiting manner of the invention, the electronic processing means 32 may be implemented: in the form of a programmable processing unit such as a microcomputer executing a signal processing program E (t) and S (t) according to the invention and 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 onboard program for processing signals E (t) and S (t) according to the invention, or - in the form of an electronic card comprising a specific electronic circuit of the ASIC type, specially designed to execute a signal processing program E (t) and S (t) in accordance with the invention. 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. in general, 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 tracking [block 45] of the central peak (P) resulting from the intercorrelation of the measurement signals. The intercorrelation function is a mathematical function that is known per se. In the continuous time domain, this cross-correlation function Icc (t), when applied to the measurement signals E (t) and S (t), is defined by the relation: (1) Icc (t) = E (t) S (t) = JE (r) S (t + z) di The intercorrelation finally consists of calculating the overlap integral between the signals E (t) and S (t) by temporally shifting a signal compared to each other.
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 Icci mise en oeuvre est de préférence définie par la relation suivante : n-I (2) Icci = EEkSi+k k=0 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, 25 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 la fréquence d'échantillonnage (fe), une succession d'échantillons discrets Ek et Sk. In practice, the intercorrelation which is calculated in the algorithm of FIG. 5 [block 43] is a discretized intercorrelation function. The discrete intercorrelation function Icci implemented is preferably defined by the following relation: nI (2) Icci = EEkSi + kk = 0 where: Ek are the values resulting from the sampling of the measurement signal E (t) [ vector size n]; Sk are the values resulting from the sampling of the measurement signal S (t) [vector of size n]; n is the total number of samples Ek and Sk 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 algorithm of Figure 5 will now be detailed. Blocks (40): Analog / Digital Conversion In order to calculate a discretized intercorrelation function, continuous measurement signals E (t) and S (t) are sampled at a sampling frequency (f). ) predefined. These sampling operations are usually carried out using analog / digital converters. These converters may be integrated in the sensors 30, 31 or may be integrated in the electronic processing means 32. At the output of this sampling step, a succession of samples is obtained at the sampling rate (f e). discreet Ek and 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]. 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].
Si le système ne dérive pas, 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, ces valeurs moyennes Emoy et Smoy sont recalculées au fur et à mesure. 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) Cette étape de correction d'offset est préférentielle car elle permet d'améliorer la sensibilité de la détection de fuite, mais est toutefois facultative. Dans la suite de la description, on considérera indifféremment 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 : lntercorré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 Icc-(n-1) à Icc(n-1). 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 Icc-(n-1) à Icc (n-1): Icc-(n-1) = Eo x Sn-1 Icc-(n-2) = E0 x Sn-2 + E1 x Sn-1 lCC(n-2) = En-2 X SO + En-1 X S1 ICC(n-1) = En-1 X SO Pour chaque série de (2n-1) valeurs d'échantillons IcCj, 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 le 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. If the system does not drift, these Emoy and Smoy average values are calculated once and for all and stored in memory. On the other hand, if the system drifts, these average values Emoy and Smoy are recalculated as and when. The offset correction [block 42] consists in subtracting from each sample Ek the average value Emoy of the corresponding signal, and each sample Sk the average value Smoy of the corresponding signal. The corrected signals are obtained as follows: (Ek ù Emoy) and (Sk ù Smoy) This offset correction step is preferred because it makes it possible to improve the sensitivity of the leak detection, but it is optional. In the remainder of the description, it will be considered indifferently that the samples Ek and Sk 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 (variant of realization with offset correction). In another variant embodiment, this offset correction could be performed on the analog measurement signals E (t) and S (t), (before sampling) by means of an analog subtractor. Blocks 43 and 44: Intercorrelation / Central Peak To calculate the cross-correlation between the two discretized signals Ek and Sk, a predefined number n of successive samples is stored in memory at each calculation, and the discretized ICCJ intercorrelation function is calculated. by means of the above-described formula (2), i.e. (2n-1) successive values Icc- (n-1) to Icc (n-1). For example, and without limitation to the invention, the computation of the intercorrelation is carried out with 1000 successive samples Ek and Sk (k varying from 1 to n and n being 1000), with a sampling frequency fe of 1 kHz. For each series of n samples, an intercorrelation vector ICq is thus calculated every second. In this example, every second value 1999 is calculated Icc- (n-1) to Icc (n-1): Icc- (n-1) = Eo x Sn-1 Icc- (n-2) = E0 x Sn-2 + E1 x Sn-1 lCC (n-2) = En-2 X SO + En-1 X S1 ICC (n-1) = En-1 X SO For each series of (2n-1) d-values IcCj samples, the central peak p of the intercorrelation (peak of greater amplitude) is determined automatically and in a manner known per se. The aforementioned steps of intercorrelation and detection of the central peak, which respectively correspond to the blocks 43 and 44 of FIG. 5, are performed repetitively in successive sliding windows of n samples Ek and Sk. These successive windows (or series) n samples Ek or Sk can be non-overlapping in time (we take n samples, then n subsequent samples without overlapping between sets of samples), or can instead overlap in part. Also, for each ICq intercorrelation computation, the two calculation windows of n samples Ek and Sk that are used can be defined on the same time interval, without temporal offset between the windows (in this case the successive samples Ek and Sk of each window were all sampled at the same time), or on the contrary can be defined with a more or less significant temporal shift between the calculation windows.
La fréquence d'échantillonnage (fe) et le 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 = n 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 1s. The sampling frequency (fe) and the number n of samples characterize a parameter T, corresponding to the integration period for the discretized intercorrelation. This parameter T is defined by the following relation: T = n fe In the above example, with a sampling frequency of 1KHz, and 1000 samples in each successive window, the integration period T is 1s.
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. For each given application, it is up to the skilled person to judiciously set the values of the sampling frequency (fe) and the number n of samples, so that the integration period T is compatible with the times. system features.
Bloc 45 : Suivi et alarmeBlock 45: Tracking and alarm
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). 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).
Ce seuil (s) peut selon le cas être un seuil prédéfini de valeur constante. This threshold (s) may, depending on the case, be a predefined threshold of constant value.
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. It can also be constituted by a self-adaptive threshold which is a function of the values Ek and Sk, 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.
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...) 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 ...)
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 5 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 oeuvre de l'invention décrits dans les tests expérimentaux ci-après, le fluide (f) circulant dans le circuit 1 10 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 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 loci est calculée sur 15 (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 20 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 25 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 30 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. It is also possible to calculate several auto-adaptive detection thresholds, which are differentiated by the value of the variable N (number of central peak values used to calculate the threshold (s).) For example, a threshold (if ) in the short term is calculated on the last 10 central peaks (N = 10), a middle-term threshold (s2) is calculated on the last 100 central peaks (N = 100), a long-term threshold (s3) is calculated on the last 500 central peaks (N = 500) When the amplitude variation of the central peak is negative and 5 exceeds in absolute value one of the thresholds (si) to (s3), an alarm is triggered characteristic of this threshold Experimental Results In all the examples of implementation of the invention described in the experimental tests below, the fluid (f) flowing in the circuit 1 10 is a liquid (in this case water); sensors 30 and 31 used are electromagnetic flowmeters having a precision of the order of e 0.5% The sampling frequency (f e) of the measurement signals E (t) and S (t) is equal to 1 KHz. The number of samples (n) for the intercorrelation is equal to 1000. The discretized loci intercorrelation function is calculated on 15 (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 modified fluid circuit 1 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 allows to cause a leak in the circuit 1 in the context of experiments of the invention. Of course, 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. In a first step, the installation of FIG. 2 was operated without leakage, by supplying circuit 1 with water with an average flow rate of the order of 295 m 3 / h, the leakage valve 14 being closed. .
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 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 la figure 8, l'intercorrélation des signaux de mesure échantillonnés Ek 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 la figure 11, l'intercorré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 Sk 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 l'intercorré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. FIGS. 6 and 7 show, for a short period of operation time without leakage, the raw data Ek and Sk coming from the sampling of the measurement signals E (t) and S (t) (signals at the output of the block 42 after offset correction). FIG. 8 shows the cross-correlation of the sampled measurement signals Ek and Sk of FIGS. 6 and 7. In a second step, a leak with a flow rate of the order of 1.2 m 3 / h is generated, 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 Ek and Sk coming from the sampling of the measurement signals E (t) and S ( t) (signals output from block 42 after offset correction). FIG. 11 shows the cross-correlation of the sampled measurement signals Ek and Sk of FIGS. 9 and 10. The comparison on the one hand of FIGS. 6 and 7 (measurement signals Ek and Sk without leakage) with another Figures 9 and 10 (measurement signals Ek and Sk with leakage) show that it is very difficult to discriminate in these signals a change to characterize the appearance of a leak. On the other hand, the comparison of FIGS. 8 and 11 resulting from the leak-free and leak-free intercorrelation shows a notable modification of the result of the intercorrelation function, and in particular a negative variation of the minimum, which makes it possible to characterize and detect the appearance of the leak. Sensitivity of the Leak Detection In order to characterize the sensitivity of the leak detection, it has been caused, by appropriately opening the leakage valve 14 of the installation of FIG. 2, five successive leaks having an increasing flow rate.
Ces fuites successives sont caractérisées dans le tableau I ci-après. Tableau I Fuite Début Fin Débit du Quantité Débit de Rapport fuite fuite 3 d'eau fuite débit de (s) (s) circuit (m /h) pendant la m /h fuite/débit fuite ( 3 ) du circuit F1 247 297 295 8 L 0,58 1,9.10"3 F2 372 412 295 10 L 0,72 3,0.10"3 F3 497 527 295 10 L 1,2 4,1.10-3 F4 597 620 295 10 L 1,6 5,3.10"3 _ 697 712 295 10 L 2,4 8,1.10"3 F5 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 Os 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 la 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. These successive leaks are characterized in Table I below. Table I Leakage Start End Quantity Flow Rate of Leakage Water leak leakage rate of flow (s) circuit (m / h) during m / h leak / leak flow (3) of F1 circuit 247 297 295 8 L 0.58 1.9.10 "3 F2 372 412 295 10 L 0.72 3.0.10" 3 F3 497 527 295 10 L 1.2 4.1.10-3 F4 597 620 295 10 L 1.6 5.3.10 In FIGS. 12 and 13, the measurement signals E (t) and S (t) delivered by the sensors 30 and 31 for the period are shown in FIGS. 12 and 13. of time between 470s and 570s. FIG. 14 shows the evolution over time of the amplitude of the central peaks of the intercorrelation (output of block 44) for the time period between Os and 730s. In this figure 14, the parts of the curve for discriminating the five successive leaks are indicated by the arrows F1, F2, F3, F4 and F5. The analysis of 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 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 m3 / 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.
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). Finally, this experiment made it possible to validate that it was possible to detect leaks of very small amplitude with respect to the main flow of the fluid at the inlet of the circuit 1. In the case of the first induced leak which has the lowest amplitude it was possible to detect this leak with a sensitivity of the order of a thousandth (1.9.10 "3), whereas comparatively the measuring accuracy of the flowmeters used was of the order of 0.5%. FIG. 3 shows a fluid circuit 1 modified for experimental purposes only to cause thermal disturbances in the circuit 1. This circuit 1 of FIG. 3 differs from that of FIG. 2 only in that it has been added on one of the pipes 110 of the circuit 1: a heating element 16 (for example a heating resistor) which makes it possible to locally heat the pipe 110 and thereby the portion of fluid circulating locally in this pipe. 110, two temperature sensors 17a and 17b making it possible to measure the temperature variations of the fluid flowing in said tubing 110 equipped with the heating element 16. In a first phase, the installation of FIG. cause leakage (valve 14 closed).
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. 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 (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 time measured downstream of the heating element 16 by the sensor 17b (curve B). 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.
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 Début Fin fuite Débit du Quantité Débit de Rapport fuite (s) circuit (m3/h) d'eau fuite débit de (s) pendant la (m3/h) fuite/débit fuite du circuit 510 525 0,7 m3/h 0,44 L 0,11 1,7.10"' On a représenté sur la 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. In a second phase, the heating of the tubing 110 is stopped for a time sufficient to return to substantially identical fluid temperature conditions throughout the circuit 1. In a third phase, the tubing 110 is heated again by means of the heating element 16, then a leak is generated by opening the valve 14. The characteristics of the leak are given in Table II below: Table II Start End leak Quantity Flow Rate Report leak (s) circuit (m3 / h ) water leak flow rate (s) during the (m3 / h) leakage / leak rate of the circuit 510 525 0.7 m3 / h 0.44 L 0.11 1.7.10 "Figure 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 fluid temperature as a function of 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 In FIGS. 21 and 22, respectively, the measurement signals E (t) and S (t) are shown during this leakage heating test, and in FIG. 23 the amplitude evolution over time. central peaks of intercorrelation. Referring to 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.
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. Applications de l'invention L'invention trouve son application à la détection de l'apparition de 30 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. Hydraulic Disturbances FIG. 4 shows a fluid circuit 1 modified for 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 a valve 18 has been added to one of the tubings 110 of the circuit 1, making it possible to close or open the said tubing 110. To generate a hydraulic disturbance, the valve 18 is closed. FIGS. the measurement signals E (t) and S (t) during this test with closing of the valve 18, and in FIG. 26 the evolution over time of the amplitude of the central peaks of the intercorrelation. 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. Knowing that the effect of the hydraulic disturbance on the results of intercorrelation (positive variation of the amplitude of the central peaks cross-correlation) is inversed with respect to the effect of a leak (negative variation in the amplitude of the central intercorrelation peaks), it can be concluded that the follow-up of the evolution over time of the amplitude of the Central peaks resulting from the cross-correlation of the measurement signals E (t) and S (t) still make it possible to detect the occurrence of a leak in the circuit, even in the presence of hydraulic disturbances. Applications of the invention The invention finds application in the detection of the appearance of leakage (s) of fluid in any circuit within which a fluid circulates, said fluid being able to be a liquid or a gas, or a gas / liquid mixture.
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. Preferably, 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 the fluid undergoes disturbances modifying its volume (in particular thermal disturbances, hydraulic, ...). 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. More particularly, 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.
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...) In this case, 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.)
Claims (15)
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0513405A FR2895508B1 (en) | 2005-12-28 | 2005-12-28 | DETECTION OF A FLUID LEAK IN A CIRCUIT BY INTERCORRELATION |
| EP06847154A EP1966580A2 (en) | 2005-12-28 | 2006-12-28 | Detection of a fluid leakage in a circuit by intercorrelation |
| PCT/FR2006/002896 WO2007074244A2 (en) | 2005-12-28 | 2006-12-28 | Detection of a fluid leakage in a circuit by intercorrelation |
| FR0611471A FR2895509B1 (en) | 2005-12-28 | 2006-12-28 | DETECTION OF A FLUID LEAK IN A CIRCUIT BY INTERCORRELATION |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0513405A FR2895508B1 (en) | 2005-12-28 | 2005-12-28 | DETECTION OF A FLUID LEAK IN A CIRCUIT BY INTERCORRELATION |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2895508A1 true FR2895508A1 (en) | 2007-06-29 |
| FR2895508B1 FR2895508B1 (en) | 2008-03-07 |
Family
ID=37451189
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR0513405A Expired - Fee Related FR2895508B1 (en) | 2005-12-28 | 2005-12-28 | DETECTION OF A FLUID LEAK IN A CIRCUIT BY INTERCORRELATION |
Country Status (3)
| Country | Link |
|---|---|
| EP (1) | EP1966580A2 (en) |
| FR (1) | FR2895508B1 (en) |
| WO (1) | WO2007074244A2 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102305552A (en) * | 2011-08-16 | 2012-01-04 | 东北大学 | Steel ladle bottom blowing dusting bleed-out detection device and bleed-out detection method |
| WO2012044435A1 (en) * | 2010-09-30 | 2012-04-05 | Siemens Aktiengesellschaft | Pipeline leak location system and method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE2603715A1 (en) * | 1976-01-31 | 1977-08-11 | Rolf Prof Dr Ing Isermann | Pipeline leakage detection and location - by measuring pressure gradients on section ends and pressure differential using stated formula |
| DE3607913C1 (en) * | 1986-03-10 | 1987-08-13 | Fraunhofer Ges Forschung | Method for locating leaks in pipelines |
| US5343737A (en) * | 1992-09-22 | 1994-09-06 | Joseph Baumoel | Method and apparatus for leak detection and pipeline temperature modelling method and apparatus |
| US5544074A (en) * | 1992-01-16 | 1996-08-06 | Kabushiki Kaisha Toshiba | Method and apparatus for detecting the position of an abnormal site of a buried pipe |
| JP2005023763A (en) * | 2003-07-04 | 2005-01-27 | Yamatake Corp | Unknown water generation distribution inferring device, its method, and program |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2727205A1 (en) * | 1994-11-22 | 1996-05-24 | Gaz De France | METHOD AND DEVICE FOR DETECTING AND LOCATING LEAKS IN A FLUID-CONDUCTED DRIVE |
| 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 |
-
2005
- 2005-12-28 FR FR0513405A patent/FR2895508B1/en not_active Expired - Fee Related
-
2006
- 2006-12-28 EP EP06847154A patent/EP1966580A2/en not_active Withdrawn
- 2006-12-28 WO PCT/FR2006/002896 patent/WO2007074244A2/en active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE2603715A1 (en) * | 1976-01-31 | 1977-08-11 | Rolf Prof Dr Ing Isermann | Pipeline leakage detection and location - by measuring pressure gradients on section ends and pressure differential using stated formula |
| DE3607913C1 (en) * | 1986-03-10 | 1987-08-13 | Fraunhofer Ges Forschung | Method for locating leaks in pipelines |
| US5544074A (en) * | 1992-01-16 | 1996-08-06 | Kabushiki Kaisha Toshiba | Method and apparatus for detecting the position of an abnormal site of a buried pipe |
| US5343737A (en) * | 1992-09-22 | 1994-09-06 | Joseph Baumoel | Method and apparatus for leak detection and pipeline temperature modelling method and apparatus |
| JP2005023763A (en) * | 2003-07-04 | 2005-01-27 | Yamatake Corp | Unknown water generation distribution inferring device, its method, and program |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012044435A1 (en) * | 2010-09-30 | 2012-04-05 | Siemens Aktiengesellschaft | Pipeline leak location system and method |
| US8850871B2 (en) | 2010-09-30 | 2014-10-07 | Siemens Aktiengesellschaft | Pipeline leak location using ultrasonic flowmeters |
| US9791345B2 (en) | 2010-09-30 | 2017-10-17 | Siemens Aktiengesellschaft | Pipeline leak location using ultrasonic flowmeters |
| CN102305552A (en) * | 2011-08-16 | 2012-01-04 | 东北大学 | Steel ladle bottom blowing dusting bleed-out detection device and bleed-out detection method |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2007074244A3 (en) | 2007-08-23 |
| WO2007074244A2 (en) | 2007-07-05 |
| FR2895508B1 (en) | 2008-03-07 |
| EP1966580A2 (en) | 2008-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR2574545A1 (en) | METHOD FOR DETECTING LEAKS IN AN OPERATING HEAT EXCHANGER | |
| FR2647549A1 (en) | METHOD AND DEVICE FOR MEASURING QUALITIES OF A POLYPHASIC FLUID | |
| EP3658763B1 (en) | Method and device for detecting ignition in a combustion chamber of a rocket motor, method for starting a rocket motor, computer program, recording medium and rocket motor | |
| FR2461249A1 (en) | METHOD FOR MEASURING LEAKAGE FLOWS | |
| EP3325934A1 (en) | Method for checking the tightness of sealed products and leak-detection facility | |
| EP0449753B1 (en) | Method and apparatus for inspecting metallic pipes by means of eddy currents | |
| EP1799869B1 (en) | Method for controlling a molten metal bath bubbling in a metallurgical vessel and a device for carrying out said method | |
| FR3028950A1 (en) | METHOD AND DEVICE FOR PERMEATION MEASUREMENT BY MASS SPECTROMETRY | |
| FR2803032A1 (en) | METHOD AND DEVICE FOR MEASURING CIRCULATING FLUID FLOW IN A PIPE | |
| EP1966580A2 (en) | Detection of a fluid leakage in a circuit by intercorrelation | |
| EP2097729A2 (en) | Method and device for detecting and/or quantifying water leaks | |
| FR2895509A1 (en) | Fluid circuit e.g. electric arc furnace panel cooling circuit, control method for e.g. iron and steel industry, involves carrying offset correction by removing mean value at measurement signals, and calculating inter-correlation of signals | |
| WO2000043762A1 (en) | Method and device for detecting the formation of a material deposit on a heat flow sensor | |
| EP0080406A1 (en) | Method and apparatus for interrogating a combustible gas concentration sensor | |
| EP3032553A1 (en) | Diagnostic device for immersed electrical transformer and electrical transformer comprising such a device | |
| WO2018162249A1 (en) | Method for measuring a speed of a fluid | |
| CH655781A5 (en) | METHOD AND DEVICE FOR MONITORING AND CONTROLLING AN EVAPORATOR. | |
| EP3593147B1 (en) | Method for measuring a speed of a fluid | |
| FR2989459A1 (en) | FUEL CONSUMPTION MEASUREMENT SYSTEM OF AN ENGINE AND METHOD OF CONTROLLING SUCH A SYSTEM | |
| EP3942375B1 (en) | Method for detecting anomalies in a water treatment facility using an apparatus for injecting oxygen into a water treatment retention basin | |
| FR2549576A1 (en) | METHOD AND DEVICE FOR LEAKAGE CONTROL ON A FLUID DISTRIBUTION NETWORK | |
| EP2489995A1 (en) | Detector of the presence of a liquid. | |
| EP3728935B1 (en) | Method for detecting the presence of gas hydrate in a pipe intended for transporting hydrocarbon fluid | |
| CA2351633A1 (en) | Installation for controlling the sealed condition of water-gas heat exchangers for industrial furnaces | |
| FR3075970A1 (en) | DEVICE FOR VERIFYING THE INTEGRITY OF A SHORT-CIRCUIT DETECTION SYSTEM |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| TP | Transmission of property | ||
| ST | Notification of lapse |
Effective date: 20120831 |