Brevet FR2481798A1 — publié le 19811106

La présente invention concerne un procédé permettant de s'assurer de l'état de l'écoulement d'un métal liquide et ce dans une large mesure à partir de la direction de dé- placement ou analogue d'un maximum d'une fonction de corré- lation mutuelle de signaux de fluctuation apparaissant aux sorties de capteurs électromagnétiques disposés dans le métal liquide en écoulement. Sans toutefois y être limitée, la présente invention concerne notamment un procédé spécia- lement approprié, à titre d'exemple, pour la mesure de grandeurs caractéristiques (mesure de la vitesse d'écoule- ment) du coeur d'un réacteur nucléaire refroidi par un métal liquide, pour la détection d'anomalies du coeur et pour le diagnostic de la cause de l'anomalie (détection de vides, auto-étalonnage/diagnostic du détecteur) et analogues. Jusqu'à présent on connaissait par exemple un cap- teur d'écoulement ou de circulation de courants de Foucault, en tant qu'instrument permettant la mesure élec- tromagnétique de la vitesse d'écoulement d'un fluide élec- triquement conducteur. Un tel capteur d'écoulement utilise, comme cela est représenté sur la figure annexée à la présente demande, un détecteur 7 constitué par une bobine d'excitation 6 et par deux bobines de détection 4, 5 disposées de part et d'autre de la bobine d'excitation 6 (c'est-à-dire formant un ensem- ble à trois bobines), et le détecteur 7 est placé suivant la direction de l'écoulement d'un fluide conducteur. On connaît également ce qu'on appelle un système à cinq bobines, dans lequel des bobines d'excitation sont en outre disposées à l'extérieur des bobines de détection. Dans tous les cas, un champ magnétique à courant alternatif produit par la bobine d'excitation 6 est tel que le flux magnétique 0 est réparti de façon symétrique par rapport à l'axe lorsque le fluide F est statique (courbe 1). Mais, lorsque le fluide F se dé- place suivant la direction indiquée par une flèche, le flux magnétique 0 est déformé sur le côté aval (courbe 2). Par conséquent une tension S1 induite dans la bobine de détection 4 sur le côté amont diminue lorsque la vitesse d'écoulement augmente, mais augmente lorsque la vitesse d'écoulement diminue, tandis qu'une tension S2 induite dans la bobine de détection 5 sur le côté aval présente un comportement inver- se par rapport à la tension Si. La différence entre S2 et S1 fournit un signal dépendant de la vitesse d'écoulement, c'est-à-dire un signal eF de vitesse d'écoulement. C'est le principe de la mesure dans le capteur d'écoulement classique. Cependant il existe une "fluctuation" dans les si- gnaux de sortie S1 et S2 des bobines de détection 4 et 5 et cette fluctuation est pulsatoire en rapport avec la moyenne dans le temps du signal. En déterminant la fonction de corrélation mutuelle R1,2 (T) des fluctuations de S1 et S on a trouvé que la vitesse d'écoulement pouvait être obtenue en divisant la distance réelle L entre les deux bobi- nes de détection par un temps de retard T0 au niveau du maximum de R1,2 (T), en se basant sur le concept selon le- quel le temps de retard T est inversement proportionnel à la vitesse d'écoulement. C'est sur cette base que la deman- deresse a proposé un nouveau capteur d'écoulement (demande de brevet japonais déposée sous le numéro provisoire 128363/1978). Conformément à ce système, la vitesse d'écoule- ment peut être obtenue indépendamment de l'amplitude du signal lui-même et par conséquent peut être appliquée pour réaliser l'auto-étalonnage du capteur d'écoulement. Cepen- dant, étant donné que l'on ne prend en compte que la vites- se d'écoulement, on n'a encore découvert aucun procédé per- mettant de s'assurer dans une large mesure de l'état de l'écoulement, incluant l'état des vides. De façon classique, on connaît un procédé de détec- tion des vides, permettant de détecter ces derniers en utilisant la valeur efficace de la fluctuation du signal de vitesse d'écoulement. Cependant ce procédé a pour incon- vénients le fait que la perturbation, qui n'a aucune rela- tion avec les vides, ne peut pas être éliminée et que par conséquent le niveau du bruit de fond est susceptible d'aug- menter, ce qui réduit la sensibilité de détection des vides. Au contraire, par exemple dans le domaine des réacteurs nucléaires refroidis par métal liquide, on a assis- té à une demande pressante concernant le développement d'une telle méthode, rendant possible de s'assurer de façon précise de l'état d'écoulement d'un métal liquide, au moyen d'une association de la détection de la vitesse d'écoulement à la détection des vides, en vue de contrôler la vitesse d'écoulement dans un sous-ensemble du combustible, de détec- ter une anomalie du coeur ou bien de diagnostiquer la cause d'une telle anomalie, ladite méthode permettant en outre un auto-diagnostic du détecteur. C'est pourquoi un objet de la présente invention est de fournir un procédé de contrôle de l'état d'écoulement d'un métal liquide, rendant possible la détection simultanée de la vitesse d'écoulement et des vides dans le métal liqui- de en association réciproque. Un autre objet de la présente invention est de four- nir un procédé rendant possible une détection correcte de l'état d'écoulement d'un métal liquide par détection de la vitesse d'écoulement du métal liquide et mesure simultanée de la présence ou de l'absence de vides en utilisant la méthode de corrélation, le pourcentage des vides, la varia- tion de la vitesse d'écoulement, etc. Un autre objet de la présente invention est de four- nir un procédé permettant de contrôler l'état d'écoulement d'un métal liquide et permettant de s'assurer dans une large mesure de l'état d'écoulement et d'effectuer un auto-diagnos- tic des défauts d'un détecteur. On va maintenant supposer que la conductivité dimi- nue de façon importante en raison de la présence de vides dans un fluide électriquement conducteur, et que par exemple le courant de Foucault circulant à travers le fluide diminue de telle sorte que le flux magnétique 0 dans les bobines de détection 4, 5 augmente, de la manière représentée par la courbe 3 sur la figure 1. Par conséquent les deux tensions S1, S2 induites dans les bobines de détection 4, 5, augmen- tent. Il faut remarquer de façon spécifique que, alors que l'incrément de S1 et le décrément de S2 résultant de la vitesse d'écoulement sont réciproquement inverses (ou possè- dent des polarités opposées),Vl'incrément de S1 et le décré- ment de S2 résultant de la présence des vides sont récipro- quement identiques (ou possèdent la même polarité). Ainsi on s'est rendu compte que différentes informations concer- nant l'état d'écoulement du fluide pouvaient être obtenues à partir du phénomène mentionné ci-dessus, au moyen de la détermination de la fonction de corrélation mutuelle R 1,2 (T) des fluctuations de S1 et S2. C'est cette idée qui est à la base de la présente invention. La figure 2 représente un diagramme de la fonction de corrélation mutuelle R1 2 (T) des fluctuations des signaux de sortie S1 et S2 des bobines de détection et montre letemps cd retard entre S1 et S2. A ce sujet, en ce qui con- cerne la vitesse d'écoulement lorsqu'aucun vide n'existe, la perturbation locale de la vitesse, présente dans l'écou- lement turbulent, se déplace suivant la direction aval en même temps que le fluide et des variations de polarités opposées agissent sur Sl et S2 comme cela a été décrit r précédemment. Par conséquent la fonction R1 2 (T) présente une corrélation négative et se présente sous la forme d'une courbe telle que la courbe 11 de la figure 2. Le temps de retard T0 au point A, auquel la valeur absolue de R 1,2 (T) est maximum, est inversement proportionnel à la vitesse d'écoulement, de sorte que l'on peut obtenir cette dernière en divisant la longueur effective L entre les bobines de détection par T. C'est le principe du capteur mentionné ci- dessus de circulation de courants de Foucault, qui avait déjà été proposé par la demanderesse. Conformément à cette invention, on peut obtenir en supplément les informations suivantes: (1) Dans le cas o il existe des vides Lorsque les vides se déplacent en même temps que l'écoulement, il apparaît des fluctuations dans les signaux de sortie Si, S2 des bobines de détection. Comme cela a été décrit précédemment, étant donné que l'on a trouvé que les vides provoquent une variation de même polarité à la fois pour S1 et S2, la fonction de corrélation mutuelle R1 2 (T) prend une valeur positive, comme représenté par une courbe 12 sur la figure 2. S'il n'existe aucun vide, la fonction de corrélation mutuelle prend une valeur négative comme re- présenté par la courbe il sur la figure 2. Par conséquent on peut déterminer la présence ou l'absence des vides en détec- tant la polarité du maximum de R1,2 (T). Conformément à cette méthode, les informations non corrélées sont réduites au minimum et peuvent être distinguées d'un signal de fluctuation dû à une fluctuation de l'écoulement, par le fait que la polarité de R1,2 (X) est opposée. Par conséquent cette méthode peut améliorer la sensibilité de détection des vides par rapport à la méthode classique mentionnée précédemment qui utilise la valeur efficace de la fluctua- tion des signaux de vitesse d'écoulement. Lorsque le pourcentage des vides augmente, la courbe 12 s'accroit comme représenté sous la forme de la courbe 13 dans la figure 2. En effet le point A se déplace suivant la direction positive (direction de C) en raison du mélange des vides et sa valeur dépend du pourcentage des vides. Par conséquent on peut détecter le pourcentage des vides au moyen de l'amplitude du maximum de polarité positive. Lorsque le pourcentage des vides augmente, la valeur maximum (PIC = IR1,2 (T-)I) de R1, 2 (T) devient suffisamment supé- rieure à une valeur (PA = jR1,2 (T0) j) obtenue en l'absence de vides. Par conséquent on peut obtenir la vitesse de dé- placement des vides en divisant la distance réelle L par un temps de retard T' qui fournit le point maximum C de pola- rité positive de la courbe 12. On peut déterminer un taux de glissement par comparaison de cette vitesse des vides à la vitesse d'écoulement. (2) Vérification de l'état d'écoulement Lorsque la vitesse d'écoulement augmente simplement, la fluctuation due à l'écoulement turbulent augmente égale- ment, tandis que le tempos de retard devient plus faible de sorte que le point A se déplace vers le point B'. Au contrai- re, lorsque la vitesse d'écoulement diminue, le point A se déplace vers le point B". Lorsque seule la turbulence de l'écoulement augmente en raison de modifications des condi- tions de surface de la paroi d'un canal d'écoulement, sans aucune modification de la vitesse moyenne d'écoulement, le point A se déplace vers le point B. Lorsque les vides inter- fèrent, le point A se déplace en direction du point C, comme cela a été décrit précédemment. Par conséquent on peut déterminer l'état de l'écoulement en examinant dans quelle direction et avec quelle amplitude le maximum A de la fonction de corrélation mutuelle se déplace. (3) Auto-diagnostic du défaut du détecteur Une anomalie d'un circuit électrique ou des bobines peut être déterminée par comparaison (par exemple par divi- sion) du signal de vitesse d'écoulement eF obtenu à partir du fonctionnement normal du détecteur à (par) une vitesse d'écoulement e obtenue à partir du temps de retard T qui fournit le maximum A de la fonction de corrélation mutuelle de S1 et S2. En d'autres termes, bien que la modification du courant d'excitation ou de l'impédance des bobines de détection provoque des variations dans le signal de sortie eF du capteur d'écoulement, elle ne provoque aucune modifi- cation de la position du maximum (T) de R1,2 (T). C'est pour cette raison qu'il est possible d'établir une distinc- tion entre le défaut ou l'anomalie du détecteur et la varia- tion de la vitesse d'écoulement elle-même, en déterminant si le rapport entre eF et e0 est égal ou non à une valeur prédéterminée. Comme cela a été décrit ci-dessus, à partir des modifications de la pDsition (temps de retard T) et de l'amplitude (intensité de corrélation PA) du maximum A de la fonction de corrélation mutuelle R 1,2 (T) des fluctuations des signaux de sortie S1 et S2 des bobines de détection, il est possible d'obtenir toutes les informations non seulement sur la vitesse d'écoulement, mais également de l'absence ou sur la présence des vides dans le fluide, sur l'état d'écoule- ment et de l'auto-diagnostic du détecteur. A titre d'exemple on a décrit ci-dessous et illus- tré schématiquement aux dessins annexés un mode d'exécution du procédé selon l'invention. La figure 1, dont il a été fait mention, est un diagramme illustrant le principe d'un capteur d'écoulement de courants de Foucault de l'art antérieur. La figure 2 est un diagramme illustrant une fonc- tion de corrélation mutuelle, considérée à la base du procé- dé selon la présente invention. La figure 3 représente un schéma-bloc d'un système permettant de mettre en oeuvre le procédé conforme à 248 17f98 l'invention. La figure 4 est un schéma-bloc d'un système utilisé lorsque le procédé selon la présente invention est appliqué à un système de contrôle prévu pour un réacteur nucléaire refroidi par un métal liquide. La figure 3 montre une forme de réalisation spécifi- que de la constitution de l'appareil requis pour la mise en oeuvre du procédé conforme à la présente invention. Le signal de sortie d'un détecteur (D) 21, qui contient une bobine d'excitation et des bobines de détection disposées comme représenté sur la figure 1, est amplifié globalement par un préamplificateur (PA) 22. La différence entre les signaux S1 et S2 des bobines de détection, après amplifica- tion, est obtenue au moyen d'un amplificateur principal (MA) 23 et délivrée en tant que signal de la vitesse d'écoulement eF. Cette structure est la même que la structure de l'art antérieur. La partie désignée par la référence 24 concerne l'objet proprement dit de la présente invention et comporte des amplificateurs (FPA) 25 de prétraitement du signal de fluctuation, un corrélateur (CORR) 26, un circuit de calcul (CALC) et un circuit de diagnostic (DIAG) 28. Les signaux de sortie S1 et S2 provenant du préamplificateur 22 sont introduits dans les amplificateurs 25 de prétraitement du signal de fluctuation de manière à sélectionner le signal de fluctuation d'une bande de fréquences nécessaire et à amplifier le faible signal de fluctuation pour l'amener à une amplitude suffisamment élevée. Ensuite la fonction de corrélation mutuelle R1,2 (X) des deux signaux est obtenue au moyen du corrélateur 26 et le temps de retard To du maximum et la variation de l'amplitude du maximum sont cal- culés à partir de la fonction R1,2 (T) par le circuit de calcul 27. La vitesse d'écoulement e est également calcu- lée à partir de T O' C'est dans le circuit de diagnostic 28 que se trouvent effectués les diagnostics suivants (i) et (ii), moyennant l'utilisation de la valeur e de la vitesse d'écoulement et de l'amplitude du maximum fournies par le circuit de calcul 27, et du signal e F de vitesse d'écoule- ment obtenu sur la base du procédé de fonctionnement clas- sique, et un signal de sortie (e x) est alors produit pour des appareils externes (par exemple un calculateur et diffé- rents dispositifs d'affichage). (i) La valeur k = eF/eO est déterminée et, moyennant l'utilisation de cette valeur en tant que coefficient d'éta- lonnage, le capteur d'écoulement de courants de Foucault est auto-étalonné. Si la valeur k devient nettement supérieu- re à la valeur initiale et si la variation est supérieure à la dérive ou modification attendue des caractéristiques du détecteur en fonction du temps, on passe à l'évaluation indiquée à la rubrique suivante (ii). Si la modification ne provient pas des vides, le détecteur est diagnostiqué comme étant anormal. (ii) La présence ou l'absence des vides (ou le taux de vides) est mesurée à partir du mouvement du maximum (point A sur la figure 2) de la fonction de corrélation. Si cela est nécessaire, on mesure la vitesse de déplacement des vides. De même on distingue des modifications dans l'état d'écoulement. La figure 4 montre une autre forme de réalisation, dans laquelle le procédé selon la présente invention est appliqué à un système de contrôle du coeur d'un réacteur refroidi par un métal liquide. Le coeur 32 du caisson 31 d'un réacteur possède une densité d'énergie élevée et la tempéra- ture du sodium liquide 33 servant de fluide de refroidisse- ment est élevée de sorte que les ensembles de combustible doivent être contrôlés individuellement à des fins de sécu- rité. Des instruments 34, tels que des thermomètres et des capteurs d'écoulement, sont généralement prévus à cet effet au-dessus du coeur. Si en tant qu'appareillage monté au- dessus du coeur, on utilise des capteurs d'écoulement de courants de Foucault, les signaux provenant des détecteurs atteignent le préamplificateur (PA) 36 par l'intermédiaire des conducteurs 35 et les signaux de vitesse d'écoulement eF sonit délivrés par des amplificateurs principaux (MA) 37. D'une manière générale ces signaux de sortie eF sont explo- rés de façon séquentielle par un multiplexeur (MPX) 38 et ils sont délivrés à un calculateur (COM) 39 par l'intermé- diaire d'un conducteur 44. La structure décrite précédemment est la même que la structure de l'art antérieur. Lorsqu'on met en oeuvre la présente invention, des signaux de sortie du préamplificateur (signaux S1, S2 des bobines de détection) sont envoyés directement au multi- plexeur 38 et sont transmis en même temps que le signal de vitesse d'écoulement e F à un système (X) 40 de contrôle du coeur par l'intermédiaire de conducteurs 41, 42, de sorte que le fonctionnement du multiplexeur 38 est commandé dans son ensemble par le calculateur 39. Le système 40 de con- trôle du coeur correspond à un système désigné par la réfé- rence 24 sur la figure 3. Sur la figure 4, en général les détecteurs 34 sont prévus en un nombre de quelques douzai- nes à quelques centaines, nombre qui correspond au nombre des sous-ensembles de combustible. Par conséquent, les nom- bres des préamplificateurs 36 et des amplificateurs princi- paux 37 sont identiques à ces nombres. Cependant, étant donné que le multiplexeur 38 peut explorer de façon séquen- tielle un grand nombre de signaux, la sortie du multiplexeur peut être au moins un canal et au maximum plusieurs canaux, pour chaque conducteur 41, 42 ou 44. Par conséquent le sys- tème 40 de contrôle du coeur selon la présente invention peut être prévu au moins en un exemplaire et au maximum en plusieurs exemplaires afin d'accroître la vitesse de traite- ment et d'obtenir une redondance. Comme cela a été décrit précédemment, le système de contrôle du coeur selon la présente invention fournit une quantité importante d'informations avantageuses pour la sécurité du réacteur ainsi que pour le fonctionnement et le contrôle de ce dernier. Cependant, en général les procédu- res de traitement des données prennent un certain temps de sorte que l'on utilise de préférence les procédures suivan- tes pour des applications pratiques: (a) Pendant le fonctionnement normal Les signaux de vitesse d'écoulement eF provenant des quelques douzaines ou centaines d'amplificateurs princi- paux 37 sont explorés séquentiellement par le multiplexeur 38 et sont envoyés au calculateur 39, mais ne sont pas en- voyés au système 40 de contrôle du coeur. Ceci est identique au système classique. (b) Inspection périodique et étalonnage du détecteur 248 1798 Ces opérations doivent être effectuées au bout d'un intervalle de temps prédéterminé (allant par exemple d'une semaine à plusieurs mois). Etant donné que ces procédures ne concernant pas la sécurité du réacteur, on pourrait utiliser une durée plus importante de mesure (quelques minutes pour chaque sous-ensemble) afin d'améliorer suffisamment la pré- cision de la mesure. Les signaux provenant des détecteurs situés à l'extérieur de chaque sous-ensemble du coeur sont explorés séquentiellement par le multiplexeur 38 sur la base d'un signal de commande 43 provenant du calculateur 39, de sorte que les signaux de sortie S11 S2 des préamplifica- teurs 36 et le signal de sortie e F des amplificateurs princi- paux 37 sont délivrés au système 40 de contrôle du coeur, et l'étalonnage et le diagnostic de panne du détecteur cor- respondant sont effectués pour chaque sous-ensemble. Les résultats sont envoyés au calculateur 39 en vue d'un calcul des caractéristiques du coeur ou bien en vue d'un affichage pour un opérateur. (c) Contrôle des anomalies du coeur Comme cela a été décrit à la rubrique (a) "pendant le fonctionnement normal", le signal de sortie e F du capteur d'écoulement est exploré séquentiellement pour chaque sous- ensemble, ce qui permet un affichage de la vitesse d'écoule- ment du fluide de refroidissement traversant chaque sous- ensemble. Séparément, les signaux S1 et S2 des bobines de dé- tection de chaque capteur d'écoulement et le signal de la vitesse d'écoulement eF sont envoyés au système 40 de con- trôle du coeur séquentiellement pour chaque sous-ensemble. A partir du rapport eF/eO et de la direction de déplacement ainsi que de l'amplitude de déplacement du maximum de la fonction R1,2 (T), on peut réaliser par conséquent une détermination telle que la présence ou l'absence d'une panne du détecteur, l'anomalie du débit (ou de la vitesse d'écou- lement) du fluide de refroidissement pour les sous-ensembles et une détection précoce de la défaillance de barreaux de combustible par détection des vides, et les signaux de sortie correspondants sont envoyés au calculateur. Etant donné que le système de contrôle du coeur selon l'invention possède un temps de traitement relativement long, il est également possible d'introduire les signaux de sortie de détection S1, S2 d'un sous-ensemble correspondant,au système de contrôle du coeur et de les faire analyser par ce sys- tème, par l'intermédiaire du multiplexeur 38 sur la base des signaux de commande 43, 45 provenant du calculateur 39, uniquement lorsqu'une anomalie est observée dans le signal de sortie eF du capteur d'écoulement, au lieu de soumettre chaque fois des données au contrôle. Etant donné qu'un fonctionnement de sécurité et un contrôle du coeur peuvent être effectués inclusivement par le système de contrôle du coeur selon la présente invention, comme décrit ci-dessus, on peut obtenir une détection pré- coce d'une anomalie et la fiabilité d'une telle détection peut être remarquablement accrue. Comme cela a été décrit précédemment, la présente invention est conçue de manière à détecter la position, la polarité et l'amplitude du maximum de la-fonction de corré- lation mutuelle des fluctuations de chaque signal de sortie d'une bobine de détection. Par conséquent il est possible de mesurer la vitesse d'écoulement du métal liquide, de façon associée, au moyen de la détection des vides dans le métal liquide. En outre la mesure et la détection peuvent être effectuées simultanément. Par conséquent l'état de l'écoulement (variation de la vitesse d'écoulement, varia- tion d'une fluctuation de 1' écoulement.--uélange de vides, etc.) peut être analysé plus correctement et peut être mesuré in- dépendamment de perturbations (bruits, etc.) autres que des vides. En conséquence la sensibilité de détection des vides peut être nettement améliorée par rapport au procédé classique qui fournit simplement la valeur efficace du signal de fluctuation. Lorsque le taux de vides est impor- tant (plus de quelques pour cent), on peut également mesu- rer la vitesse des vides. En outre, on peut réaliser une dis- crimination de défauts du détecteur en comparant la vitesse d'écoulement obtenue en utilisant le procédé selon la pré- sente invention et le signal de sortie du capteur d'écoule- ment obtenu au moyen du procédé classique. Sur la base de ces caractéristiques remarquables, on peut utiliser effica- cement la présente invention pour le contrôle de la vitesse d'écoulement du sous-ensemble, pour la détection de l'anoma- lie et pour la discrimination des causes de l'anomalie dans un réacteur refroidi par un métal liquide. Bien que la présente invention ait été décrite en rapport avec des formes de réalisation préférées, il appa- raîtra de façon évidente au spécialiste que de nombreuses modifications peuvent y être apportées sans sortir du cadre de la présente invention. REVENDICATIONS 1. Procédé pour contrôler l'état de l'écoulement d'un métal liquide, du type consistant à disposer suivant la direction d'écoulement du métal liquide un détecteur (7 21) constitué par une bobine d'excitation (6) parcourue par un courant alternatif, et au moins deux bobines de'détec- tion (4, 5) disposées des deux côtés de la bobine d'excita- tion (6), et à détecter un maximum (A) de la fonction de corrélation mutuelle (R1,2 (T)) de signaux de fluctuation apparaissant dans lesdites bobines de détection, caractérisé en ce que l'on détecte la direction de déplacement du maximum (A) de la fonction de corrélation mutuelle de manière à déterminer une modification de l'état de l'écoulement du métal liquide. 2. Procédé selon la revendication 1, caractérisé en ce que la variation de l'état de l'écoulement est une varia- tion de la vitesse de l'écoulement et/ou du pourcentage des vides dans le métal liquide. 3. Procédé pour contrôler l'état de l'écoulement d'un métal liquide, du type consistant à disposer suivant la direction d'écoulement du métal liquide un détecteur (7 21) constitué par une bobine d'excitation (6) parcourue par un courant alternatif, et au moins deux bobines de détec- tion (4, 5) disposées des deux côtés de la bobine d'excita- tion (6), et à détecter un maximum (A) de la fonction de corrélation mutuelle (R1,2 (T)) de signaux de fluctuation apparaissant dans lesdites bobines de détection, caractérisé en ce que l'on détecte la polarité du maximum (A) de la fonction de corrélation mutuelle de manière à déterminer ainsi la présence ou l'absence de vides. 4. Procédé selon la revendication 3, caractérisé en ce que l'on détecte l'amplitude du maximum (C) de polarité positive de manière à déterminer ainsi le pourcentage des vides. 5. Procédé selon la revendication 3, caractérisé en ce que l'on détecte un temps de retard (To) de manière à ob- tenir le maximum (C) de polarité positive et que l'on déter- mine la vitesse d'écoulement des vides dans le métal liquide à partir du temps de retard (T) ainsi détecté et de la dis- tance (L) présente entre lesdites bobines de détection (4,5).