FR 2925690 A1 20090626 FR 0709045 A 20071221 SETVAL54.FRD Contrôle non destructif, en particulier pour des tubes en cours de fabrication ou à 5 l'état fini. L'invention concerne le contrôle non destructif de matériaux, notamment pour des tubes en fabrication. 10 On connaît déjà différentes propositions, sur lesquelles on reviendra, tendant à utiliser les réseaux neuronaux dans le cadre du contrôle non destructif de matériaux. Mais cet existant n'est pas de nature à fonctionner en milieu industriel, sur des équipements déjà en service, en temps réel, et tout en permettant une classification au vol des imperfections selon leur nature, de sorte que l'on puisse remédier rapidement à un 15 problème survenu en phase de production. La demande de brevet français n° 0605923 non publiée traite du contrôle non destructif. Un objet de l'invention est d'améliorer la situation vers un système qui puisse: 20 - être utilisé en milieu industriel et implanté facilement sur des équipements déjà existants dans ce milieu. - être utilisé en temps réel, c'est-à-dire donner un diagnostic rapide, notamment à une vitesse suffisante pour ne pas ralentir la vitesse de production globale, et -permettre une classification des imperfections selon leur nature, à partir d'une faible 25 quantité d'informations, de manière à connaître leur gravité et à permettre la détermination de la cause technique à l'origine de l'imperfection et ainsi remédier rapidement au problème en phase de production. Selon un premier aspect, il est proposé un dispositif qui forme outil d'aide à 30 l'exploitation, pour le contrôle non destructif de tubes (ou autres produits sidérurgiques) en cours et en sortie de fabrication. Un tel outil est destiné à tirer des informations sur d'éventuelles imperfections du produit. Des capteurs ultrasonores émetteurs sont excités sélectivement selon une loi de temps choisie. Des signaux de retour sont captés par des capteurs ultrasonores récepteurs formant un agencement de géométrie choisie, monté en couplage ultrasonore avec le tube par l'intermédiaire d'un milieu liquide. Enfin, il y a généralement un mouvement relatif de rotation/translation entre le produit et l'agencement de transducteurs. L'outil d'aide à l'exploitation proposé comprend : - un convertisseur, capable d'isoler sélectivement une représentation numérique de possibles échos dans des fenêtres temporelles désignées, en fonction du mouvement relatif de rotation/translation, et d'en tirer une image d'imperfections possibles dans le produit, ladite représentation comprenant l'amplitude et le temps de vol d'au moins un écho, et de générer un graphe 3D parallélépipédique, - un bloc de transformation capable de générer une image 3D d'imperfections possibles dans le tube à partir du graphe 3D et d'une base de données, - un filtre, capable de déterminer, dans les images, des zones d'imperfection présumée, ainsi que des propriétés de chaque imperfection présumée, - un étage de sortie configuré pour générer un signal de conformité ou de non conformité d'un produit. L'invention peut aussi se placer au niveau d'un dispositif de contrôle non destructif de tubes (ou autres produits sidérurgiques) en cours ou en sortie de fabrication, qui comprend : - un agencement de transducteurs ultrasonores de géométrie choisie, monté en couplage ultrasonore avec le tube par l'intermédiaire d'un milieu couplant, avec mouvement relatif de rotation/translation entre le tube et l'agencement de transducteurs, - des circuits pour exciter sélectivement ces éléments transducteurs selon une loi de temps choisie, et pour recueillir les signaux de retour qu'ils captent, et - un outil d'aide à l'exploitation tel que ci-dessus. Un autre aspect de l'invention s'exprime sous forme de procédé de contrôle non destructif de tubes (ou autres produits sidérurgiques) en cours ou en sortie de fabrication, comprenant les étapes suivantes : a. prévoir un agencement de transducteurs ultrasonores de géométrie choisie, monté en couplage ultrasonore avec le tube par l'intermédiaire d'un milieu couplant, avec mouvement relatif de rotation/translation entre le tube et l'agencement de transducteurs, b. exciter sélectivement ces éléments transducteurs selon une loi de temps choisie, c. recueillir les signaux de retour qu'ils captent, afin d'analyser sélectivement ces signaux de retour, pour en tirer des informations sur d'éventuelles imperfections du tube, lesdites informations comprenant l'amplitude et le temps de vol d'au moins un écho et générer un graphe 3D parallélépipédique. d. isoler sélectivement une représentation numérique de possibles échos dans des fenêtres temporelles désignées, en fonction du mouvement relatif de rotation/translation, et en tirer une image 3D d'imperfections possibles dans le tube à partir du graphe 3D parallélépipédique et d'une base de données, e. générer un signal de conformité ou de non conformité d'un produit. L'étape e peut comprendre: el. filtrer les images selon des critères de filtrage choisis, afin d'y déterminer des zones d'imperfection présumée Zcur, ainsi que des propriétés de chaque imperfection présumée, e2. former des entrées numériques de travail, à partir d'un extrait des images correspondant à une zone d'imperfection présumée Zcur, de propriétés de l'imperfection présumée dans la même zone, issues du filtre; et de données de contexte, e3. appliquer les entrées ainsi formées à au moins un agencement du genre circuit neuronal, e4. traiter numériquement la sortie de l'agencement du genre circuit neuronal selon des critères de décision choisis, pour en tirer une décision et/ou une alarme, et e5. écarter et marquer des tubes décidés non conformes par l'étape e4. D'autres aspects, caractéristiques et avantages de l'invention apparaîtront à l'examen de la description détaillée ci-après de quelques modes de réalisation non limitatifs, ainsi que des dessins annexés sur lesquels : - la figure 1 est une vue schématique en perspective d'un tube, ayant des imperfections ou défauts dits étalons; - la figure 2 est une vue latérale schématique illustrant un exemple d'installation du type "contrôle à tête tournante" sur un tube en sortie de fabrication ; - les figures 3A à 3C sont des détails de différents types de mesure d'épaisseur et de contrôle d'imperfections longitudinales et transverses ; - la figure 4 est le schéma de principe de l'électronique associée à un capteur ultrasonore en contrôle non destructif dans une installation classique ; -les figures 5A et 5B sont une vue en bout et une vue latérale d'un type particulier de 10 cellule de contrôle non destructif, communément appelé "à tête tournante" et représenté schématiquement ; - la figure 6 montre la complexité des trajets ultrasonores rencontrés dans un tube, sur un exemple simple ; - les figures 6A et 6B sont des diagrammes temporels schématiques de signaux 15 ultrasonores, pour un capteur sous incidence oblique, et pour un capteur sous incidence normale (perpendiculaire), respectivement ; - la figure 7 est un graphe montrant une représentation classique de la sélectivité d'une installation de contrôle ; - la figure 8 est un schéma de principe de l'électronique associée à un capteur 20 ultrasonore en contrôle non destructif dans un exemple d'installation susceptible de mettre en oeuvre l'invention ; - la figure 8A est un schéma fonctionnel plus détaillé d'une partie de la figure 8; - la figure 8B est un autre schéma fonctionnel plus détaillé d'une partie de la figure 8; - la figure 9 est une copie d'écran schématisée illustrant des images ultrasonores 25 numérisées d'imperfections potentielles dans un tube; - la figure 9A est une copie d'écran selon une autre orientation ; - les figures 10A à 10D sont des représentations schématiques de différents types d'imperfections selon la classification API (American Petroleum Institute) et qui constituent les données de sortie du réseau neuronal tendant à déterminer le type de 30 l'imperfection ; - la figure 11 est le schéma fonctionnel plus détaillé d'une autre partie de la figure 8; - la figure 11A est une vue de détail du bloc de transformation de la figure 11; - la figure 12 est un diagramme séquentiel illustrant le traitement d'imperfections potentielles successives dans une image ; - la figure 13 est le schéma fonctionnel d'un système de filtres; - la figure 14 est le schéma fonctionnel d'un montage à réseau neuronal tendant à 5 déterminer le type d'une imperfection dans un tube ; - la figure 15 est le schéma fonctionnel d'un montage à réseau neuronal tendant à déterminer le degré de gravité d'une imperfection dans un tube ; - la figure 16 est le schéma fonctionnel du modèle de neurone ; - la figure 17 est un exemple de fonction de transfert d'un neurone élémentaire ; et 10 - la figure 18 est le schéma général d'une installation pour la détection des défauts sur différents types de capteurs. Les dessins contiennent des éléments de caractère certain. Ils pourront donc non seulement servir à mieux faire comprendre la présente invention, mais aussi contribuer à sa définition, le cas échéant. 15 Dans la suite du présent texte, un capteur ultrasonore pourra être désigné indifféremment par les termes capteur, ou palpeur ou transducteur, bien connus de l'homme du métier. 20 Circuits neuronaux L'utilisation de réseaux neuronaux dans le cadre du contrôle non destructif de matériaux a fait l'objet de nombreuses publications, la plupart du temps assez théoriques, que l'on considérera maintenant. 25 L'article "Localization and Shape Classification of Defects using the Finite Element Method and the Neural Networks" de ZAOUI, MARCHAND et RAZEK (NDT.NETAOUT 1999, vol. IV n° abrégé 8) formule des propositions dans ce domaine. Cependant, ces propositions sont faites dans le cadre de manipulations en laboratoire, et 30 l'application décrite ne permet pas une mise en oeuvre en ligne, en milieu industriel. En outre, seule la détection par courants de Foucault est traitée, ce qui est souvent insuffisant. L'article "Automatic Detection of Defect in Industrial ultrasound images using a neural Network" de Lawson et Parker (Proc. of Int. Symposium on Lasers, Optics, and Vision for Productivity in Manufacturing I (Vision Systems: Applications), June 1996, Proc. of SPIE vol. 2786, pp. 37-47, 1996), décrit l'application du traitement d'images et des réseaux neuronaux à l'interprétation de ce que l'on appelle "scan TOFD". La méthode dite TOFD (Time of Flight Diffraction) consiste à repérer les positions du capteur ultrasonore où l'on peut observer une diffraction du faisceau sur les bords de l'imperfection, ce qui permet par la suite de dimensionner l'imperfection. Cette méthode est difficilement adaptable aux équipements de contrôle non destructif déjà existants, en particulier en milieu industriel. L'article "Shape Classification of Flaw Indications in 3-Dimensional ultrasonic Images" de Dunlop et McNab (IEE Proceedings - Science, Measurement and Technology -- July 1995 -- Volume 142, Issue 4, p. 307-312) concerne le diagnostic en termes de corrosion de pipeline. Le système permet un contrôle non destructif en profondeur et permet une étude dans les trois dimensions et en temps réel. Cependant, le système est très lent. Cela rend son utilisation en milieu industriel relativement difficile. L'article "Application of neuro-fuzzy techniques in oil pipelines ultrasonic nondestructive testing" de Ravanbod (NDT&E International 38 (2005) p 643-653) suggère que les algorithmes de détection de l'imperfection peuvent être améliorés par l'utilisation d'éléments de logique floue, mêlés au réseau de neurones. Cependant, les techniques étudiées concernent là aussi l'inspection d'imperfections de pipeline et un diagnostic sur des imperfections de corrosion. DE 42 01 502 C2 décrit une méthode pour créer un signal destiné à un réseau de neurones, mais n'apporte que peu ou pas d'enseignements sur l'interprétation des résultats, en termes de diagnostic. En outre, une nouvelle fois, seule la détection par courants de Foucault est traitée. La publication de brevet japonais 11-002626 concerne la détection d'imperfections longitudinales uniquement, et seulement par courants de Foucault. La publication de brevet No. 08-110323 se contente d'une étude en fréquence des 5 signaux obtenus par ultrasons. La publication de brevet No. 2003-279550 décrit un programme pour faire la différence entre une zone qualifiée de saine et une zone mauvaise d'un produit en utilisant un réseau de neurones. Ce programme ne va pas plus loin, et ne permet pas la classification 10 ni la localisation d'imperfections. En conséquence, l'application de ce programme peut fréquemment conduire au rejet de pièces qui seraient pourtant considérées comme bonnes si les résultats étaient interprétés par un opérateur humain. Contrôle non destructif de tubes La description détaillée ci-après est donnée essentiellement dans le cadre du contrôle non destructif de tubes en sortie de fabrication, à titre non limitatif. Comme indiqué sur la figure 1, les imperfections dans un tube T peuvent être 20 distinguées d'après leur position. Ainsi, les imperfections de surface, interne ou externe, comprennent les imperfections longitudinales LD, et les imperfections circonférentielles (ou transverses ou travers ou transversales) CD et les imperfections obliques ou inclinées ID ; par différents agencements de capteurs, on cherche à les détecter dès qu'elles s'étendent sur une longueur et une profondeur définies selon les normes ou les 25 spécifications ou cahiers de charges des clients (à titre d'exemple une valeur de longueur d'imperfection citée dans les normes est de 1/2 pouce, soit 12,7 mm environ avec une profondeur d'environ 5 % de l'épaisseur du produit contrôlé). On s'intéresse aussi aux imperfections � dans la paroi �, c'est-à-dire dans la masse MD (non visibles sur la figure 1), qui correspondent souvent à des inclusions ou des dédoublures, que l'on 30 cherche à détecter en même temps que l'on fait une mesure d'épaisseur. Les faisceaux ultrasonores sont représentés divergents sur la figure 1 pour faire comprendre la détection d'imperfections. En pratique, ils seront plutôt convergents, comme on le verra. 15 De manière classique, en contrôle non destructif par ultrasons, on utilise l'un des trois types d'installations suivants : les installations dites "à tête tournante", les installations dites "à tube tournant", et les installations à capteur encerclant multiéléments, toutes bien connues de l'homme du métier. Dans le cas d'utilisation de capteurs opérant en balayage électronique, la rotation relative tube/capteurs est virtuelle. Telle qu'utilisée ici, l'expression "mouvement relatif de rotation/translation entre le tube et l'agencement de transducteurs" couvre le cas où la rotation relative est virtuelle. Sur la figure 2, la machine de contrôle non destructif à tête tournante comprend un dispositif à ultrasons proprement dit, monté sur une enceinte à eau ou "boîte à eau" 100, que traverse le tube T à la vitesse v = 0,5 mètre par seconde, par exemple. Les capteurs ou palpeurs à ultrasons émettent des ondes longitudinales dans l'eau. Un capteur donné travaille par exemple à 1 ou quelques MHz. Il est excité répétitivement par impulsions, de forme d'onde choisie, à une cadence (ou fréquence) de récurrence Fr qui est de l'ordre de quelques kHz ou dizaines de kHz, par exemple de 10 kHz. Par ailleurs, un transducteur à ultrasons possède : - un rayonnement de champ proche, pratiquement parallèle, dans une zone dite de 20 Fresnel, siège de nombreuses interférences, dont la longueur dans l'axe du faisceau est N=0,25 D2/1 où D est le diamètre de la pastille active du transducteur, et X sa longueur d'onde de travail, et - un rayonnement de champ lointain, dans une zone dite de Fraunhofer, selon un 25 faisceau divergent d'angle 2 a, avec sin a = 1,22 X / D Les figures 3A, 3B, 3C représentent des capteurs rendus convergents au moyen d'une lentille (à ultrasons) concave, tels qu'ils sont couramment utilisés dans les applications 30 de contrôle des tubes. On utilise de préférence la zone de Fraunhofer, moins perturbée. Ainsi, pour des capteurs tels que P11 et P12, le faisceau d'ultrasons, qui est en général focalisé, s'étend au voisinage d'un plan perpendiculaire à l'axe du tube T. La détection se fait donc sensiblement en section droite. Leurs rôles sont les suivants: - ou bien leur faisceau est également perpendiculaire à l'axe du tube T dans la section 5 droite, et ils servent à la mesure d'épaisseur (par exemple P1, figure 3A); on parle alors de "palpage droit"; - ou bien leur faisceau est incident sur l'axe du tube T, en section droite, et ils servent à détecter les imperfections longitudinales (par exemple P11, figure 3B). Dans ce cas, l'angle d'incidence dans la section droite est, de préférence, choisi pour n'engendrer 10 dans le tube que des ondes ultrasonores transversales ou de cisaillement, compte tenu des caractéristiques de l'interface eau/métal du tube (en principe eau/acier). On prévoit généralement deux capteurs P11 et P12, d'incidences opposées par rapport à l'axe du tube (figure 2). 15 La machine comprend aussi des capteurs tels que P21 et P22 dont, par contre, le faisceau d'ultrasons, qui est également focalisé en règle générale, s'étend au voisinage d'un plan passant par l'axe du tube, mais incident par rapport au plan perpendiculaire à l'axe du tube T (voir capteur P21, Figure 3C). Dans ce cas, l'angle d'incidence par rapport au plan perpendiculaire à l'axe du tube est de préférence choisi pour 20 n'engendrer dans le tube que des ondes ultrasonores transversales ou de cisaillement, compte tenu des caractéristiques de l'interface eau/métal du tube (en principe eau/acier). Ces capteurs servent à détecter les imperfections transverses. On prévoit généralement deux capteurs P21 et P22, d'incidences opposées par rapport au plan perpendiculaire à l'axe du tube (figure 2). 25 Le contrôle des imperfections se fait généralement en focalisant le faisceau. Le point de focalisation se mesure par rapport au "bond", qui correspond au premier trajet aller-retour des ultrasons dans l'épaisseur du tube. Ainsi, le capteur de la figure 3A est focalisé au demi-bond, tandis que les capteurs des figures 3B et 3C sont focalisés à 30 trois-quarts de bond. Par ailleurs, le contrôle des imperfections externes se fait généralement au bond, et celui des imperfections internes au demi-bond. On note Ta le temps de présence requis pour que le palpeur puisse recevoir correctement le retour du faisceau d'ultrasons, représentatif d'une éventuelle imperfection. Ce temps Ta dépend de la somme des deux temps suivants: - d'une part le temps de propagation aller-retour d'ondes ultrasonores longitudinales, sur 5 la hauteur de la "colonne d'eau" présente entre le palpeur et le tube, sur le trajet des ultrasons, -d'autre part le temps de propagation d'ondes ultrasonores transversales, tel qu'il est requis à l'intérieur du tube pour faire le contrôle non destructif lui-même. Ce temps dépend principalement d'un choix du nombre souhaité de réflexions des ondes 10 transversales à l'intérieur de la paroi du tube. Classiquement, les palpeurs sont entraînés en rotation autour de l'axe du tube, par des moyens non représentés, à une vitesse T de l'ordre de plusieurs milliers de tours par minute, (6000 tpm par exemple). Dans le cas également connu de l'homme du métier 15 où c'est le tube qui est entraîné en rotation alors que les palpeurs ne sont pas entraînés en rotation (installation dite "à tube tournant") la vitesse de rotation du tube est de l'ordre de quelques dizaines à quelques milliers de tours par minute. On peut appeler "cellule" chaque ensemble capteur - milieu de transmission (eau) - tube. 20 Pour une cellule, on doit considérer en outre l'ouverture de faisceau Od des palpeurs à ultrasons en détection. Une ouverture peut être définie avec deux composantes (figure 1), l'une Od 1 dans la section droite du tube, l'autre Od2 dans le plan passant par l'axe du tube et le palpeur. 25 Le réglage de l'installation (en fonction de la vitesse de rotation, de la vitesse de défilement, des dimensions Od1 et Od2 et du nombre de palpeurs) doit garantir un balayage par les faisceaux ultrasonores de l'ensemble des surfaces et du volume du tube à contrôler. 30 Il est à noter que certaines normes ou cahiers des charges ou spécifications des clients imposent un recouvrement des zones balayées. Le temps d'analyse Ta est donc défini par un compromis entre: - la cadence (ou fréquence) de récurrence Fr, - en section droite du tube, la vitesse de rotation w, compte tenu de l'ouverture en détection Odl des palpeurs à ultrasons (autrement dit, compte tenu de la rotation des capteurs, la composante Odl de l'ouverture de faisceau doit permettre un temps de présence de l'imperfection en face des capteurs qui soit au moins égal à Ta ), - le long du tube, la vitesse de défilement v de celui-ci, compte tenu de l'ouverture en détection Od2 d'un palpeur à ultrasons, et du nombre NFi de palpeurs consacrés à la même fonction Fi (qui constituent donc un groupe de palpeurs), sur la périphérie du tube (autrement dit, compte tenu de l'avance du tube, la composante Od2 de l'ouverture de faisceau doit permettre un temps de présence de l'imperfection en face du capteur (ou du groupe de capteurs) qui soit au moins égal à Ta). - le nombre de palpeurs dédiés au même rôle (c'est à dire à la même fonction) , et - les temps de propagation des ondes tels définis précédemment. Classiquement, la machine comprend typiquement en tout deux capteurs tels que P11, P12 pour le contrôle des imperfections de type LD, et éventuellement ID, deux capteurs tels que P21, P22 pour le contrôle des imperfections de type CD, auxquels s'ajoute en principe un capteur tel que P1, pour la mesure de l'épaisseur du produit et le contrôle des imperfections de type MD. Chaque capteur peut être en fait un groupe de capteurs travaillant ensemble, comme on le verra. La machine possède, de façon intégrée ou à part, une électronique d'excitation et de détection associée à chacun des capteurs. Elle comprend (figure 4) un émetteur 70 d'impulsions, par exemple à 250 Volts pour l'excitation du palpeur PO monté sur la boîte à eau 100. Comme partie intégrante du système de contrôle non destructif, le palpeur à ultrasons P0, ici émetteur/récepteur, reçoit les échos consécutifs à cette excitation. Les lignes 700 et 710 transmettent respectivement l'impulsion d'excitation et le signal aux bornes du palpeur à un amplificateur 73. La sortie de l'amplificateur 73 sert à la visualisation pour l'opérateur et/ou au pilotage d'un automate de tri, capable d'écarter (en aval) les tubes non conformes. La visualisation est par exemple réalisée par un oscilloscope 750, qui reçoit comme signal la sortie de l'amplificateur 73, et comme base de temps 752 un signal d'un étage de synchronisation 753 en provenance de l'émetteur 70. Un étage de seuil 754 évite l'aveuglement de l'oscilloscope au moment de l'impulsion d'émission. Une autre sortie de l'amplificateur 73 va vers un étage de traitement de signal 760. Ce traitement comprend généralement un redressement, un lissage et un filtrage. Il est suivi d'un étage de détection ou sélecteur 762, capable d'isoler les échos significatifs, de manière connue. En détection de l'imperfection, c'est la présence d'un écho, avec son amplitude ou sa durée (donc son énergie), qui sont significatifs, dans certains créneaux temporels, essentiellement le demi-bond et le bond. Pour la détection en épaisseur, on vérifie que l'équivalent-distance de l'écart temporel entre les échos de fond respectifs correspond bien à l'épaisseur voulue du tube. Les anomalies détectées selon ces critères peuvent servir à émettre une alarme en 764, et/ou à piloter un automate 766 de tri qui évacue les tubes non conformes, en les marquant selon la ou les anomalies détectées. Matériellement dans le cas d'une installation à tête tournante (figures 5A et 5B), la cellule comprend encore, sur un support mécanique 80, la boîte à eau 100, qui loge un ensemble capteur P0, avec une connexion 701, qui réunit les lignes 700 et 710 de la figure 4. On prévoit par exemple trois roulements 81 à 83 pour centrer le tube T. Selon la technique connue (machine vendue par exemple par la société allemande GE NUTRONIK anciennement NUKEM), l'ensemble capteur PO comprend des capteurs qui tournent à quelques milliers de tours/minute autour du tube. On peut également utiliser une pluralité de capteurs répartis en anneau autour du tube. L'anneau comprend par exemple 6 secteurs de 128 capteurs ultrasonores, répartis autour de la périphérie. Les secteurs de capteurs sont alternativement légèrement décalés dans le sens de l'axe du tube. Ceci permet d'avoir un recouvrement entre deux secteurs de capteurs consécutifs longitudinalement, et réduit également les problèmes d'interférences. L'interférence intervient lorsqu'un capteur donné reçoit des échos dus à un tir effectué sur un autre capteur. A cela s'ajoute un banc (non représenté) de guidage du tube en amont et en aval du poste de contrôle non destructif, pour bien positionner le tube qui défile en continu, par rapport aux capteurs à ultrasons. Le contrôle non destructif doit se faire sur toute la périphérie du tube. Mais il est également essentiel que ce contrôle suive la vitesse linéaire v du tube en sortie de fabrication. On en vient donc à un compromis entre la vitesse linéaire v du tube, la cadence (ou fréquence) de récurrence Fr, le temps d'analyse Ta, l'ouverture de travail Od du palpeur à ultrasons en détection, et la vitesse de rotation co, le nombre de capteurs assurant la même fonction et la vitesse de propagation des ondes ultrasonores. Il est également souhaitable que la même installation puisse travailler sur toute une gamme de diamètres de tubes (et aussi d'épaisseurs de tubes), couvrant la gamme de production. Il est alors fréquent de prévoir plusieurs valeurs de la vitesse de rotation co et de la fréquence de récurrence Fr, valeurs que l'on sélectionne en fonction du diamètre du tube à traiter. On notera enfin que tout changement de fabrication implique un nouveau réglage des angles d'attaque des ultrasons de chaque capteur sur la périphérie du tube. Cette opération délicate, effectuée manuellement, prend couramment de l'ordre d'une demi-heure, temps pendant lequel la production de tubes est arrêtée. Telles sont les conditions dans lesquelles s'effectue actuellement le contrôle non destructif par ultrasons des tubes, ou d'autres produits profilés et/ou à parois minces, en sortie de fabrication. Dans le domaine du contrôle non destructif par ultrasons, on utilise souvent la terminologie suivante : - � balayage � (ou � scan �) désigne une suite de positions relatives tube/capteur, - � incrément � désigne le pas de balayage (inversement proportionnel à la fréquence de récurrence ou fréquence de tirs ultrasonores), - �Ascan � désigne le graphe de la tension électrique mesurée aux bornes d'un capteur ultrasonore, avec en abscisse le temps de vol et en ordonnée une représentation de la tension électrique, appelée également amplitude ultrasonore, - �Bscan � désigne une image relative à une valeur donnée de l'incrément, avec, en abscisse, le balayage correspondant au tir ultrasonore, éventuellement exprimé en degré d'angle du capteur par rapport à la pièce à inspecter, en ordonnée le temps de vol, et en chaque point l'amplitude ultrasonore convertie en dégradé de gris ou de couleur, - � Echodynamique � désigne un graphe avec en abscisse l'indication du tir ultrasonore et en ordonnée l'amplitude maximale relevée dans un sélecteur temporel du Ascan pour le tir correspondant, - � Cscan � désigne une image avec en abscisse et en ordonnée la position équivalente dans un espace plan du point de tir de l'onde ultrasonore et représentant, convertie en dégradé de gris, l'amplitude ultrasonore maximale pour ce tir relevée dans le sélecteur temporel considéré du Ascan ("amplitude de l'image"). Dans le cas d'un tube, un point de l'abscisse du Cscan correspond à une position sur la longueur du tube et un point de l'ordonnée correspond à une position sur la circonférence du tube. Dans le cas d'un produit plat, un point de l'abscisse du Cscan correspond à une position sur la longueur du produit plat et un point de l'ordonnée correspond à une position sur la largeur du produit plat. Par ailleurs, la demanderesse utilise dans la suite de la description les termes suivants: - � Bscan 3D parallélépipédique� qui désigne une représentation 3D comportant en outre la position du capteur sur l'axe du tube, la représentation étant considérée comme brute, la forme du tube n'apparaissant pas, - � Bscan 3D réduit� qui désigne un Bscan 3D parallélépipédique limité à une zone avec indication ultrasonore de défaut probable à l'issue des filtrages, - � Bscan 3D tube� qui possède les mêmes dimensions que le Bscan 3D parallélépipédique, les données étant représentées dans le tube inspecté, l'amplitude pouvant éventuellement constituer une dimension supplémentaire. La figure 6 est une vue en coupe longitudinale schématique d'un système formé d'un capteur, de sa colonne d'eau et du tube, et avec illustration de différents trajets ultrasonores formant échos. Elle permet de bien comprendre la complexité de ces trajets, et la difficulté de l'analyse. La figure 6A est un diagramme amplitude/temps schématique du signal ultrasonore au niveau d'un capteur qui travaille sous incidence oblique. A partir de l'instant Texcit d'excitation du capteur, on trouve un écho d'interface eau-tube à l'instant Tinterf (que l'on peut également noter TphiExter0). Sont marqués ensuite (trait pointillé vertical) l'instant TphiInter où le faisceau d'ultrasons atteint la peau interne du tube, où il se réfléchit et réfracte, ainsi que l'instant TphiExterl où le faisceau d'ultrasons atteint la peau externe du tube. Du fait de l'incidence oblique, il n'y a pas d'écho réfléchi significatif qui revient au capteur en TphiInter en l'absence d'imperfection à cet endroit. Ceci vaut aussi en TphiExterl. La figure 6B est un diagramme amplitude/temps schématique du signal ultrasonore au niveau d'un capteur qui travaille sous incidence normale. La chronologie générale des signaux est la même que pour la figure 6A (à un facteur près, lié à l'incidence). Par contre, sous incidence normale, il y a des échos significatifs en TphiInter et en TphiExterl, même en l'absence d'imperfection aux endroits concernés du tube. Actuellement, les systèmes de contrôle non destructif utilisés en production de tubes fonctionnent en faisant le rapport K entre : -l'amplitude As d'un signal provenant du tube à inspecter, et - l'amplitude A0 du signal provenant d'un défaut étalon de référence, pour le type de contrôle considéré. Ce � défaut étalon de référence � est en général défini sur un tube étalon muni d'un défaut artificiel (par exemple une entaille en U ou en V) de caractéristiques dimensionnelles choisies, par exemple conformément à une norme de contrôle non destructif, et/ou au cahier des charges d'un client. L'hypothèse implicite est que cette amplitude de signal est proportionnelle à la criticité de l'imperfection, c'est-à-dire à sa profondeur (DD). Le graphe de la figure 7 (bien connu de l'homme du métier, voir Nondestructive Testing Handbook - chapitre statistics du volume 7 publié par l'ASNT - American Society for Nondestructive Testing) représente la répartition réelle K = f(DD). Il montre qu'en réalité, la corrélation est très mauvaise (de l'ordre de 0.3 à 0.4 pour le contrôle par ultrasons). Plus précisément, sur le graphe de la figure 7, si l'on fixe l'amplitude de référence A0 (K=1) sur la valeur XL (profondeur de l'imperfection maximum acceptable) au centre de la distribution (elle- même axée sur l'oblique TDis), on voit que l'on peut encore trouver des imperfections à K = 0,5 de profondeur DD supérieure à XL. Il s'ensuit que, par prudence, on est amené à fixer A0 pour une valeur nettement plus basse que XL. Par conséquent, on écarte en production des tubes qui, pourtant, seraient en fait satisfaisants. C'est d'autant plus néfaste, économiquement, que les techniques de fabrication de tubes demeurent assez lourdes, tant en complexité qu'en besoins en énergie. La Demanderesse s'est donc attachée à améliorer la situation. La figure 8 montre un dispositif amélioré par rapport à celui de la figure 4. La sortie de l'amplificateur 73 est appliquée à un étage 761, qui numérise l'amplitude du signal issu de l'amplificateur 73, et travaille sur ce signal numérisé. Ce traitement sera décrit ci-après en référence à la figure 11. On peut conserver ensuite des étages 764 et 766 fonctionnellement semblables à ceux de la figure 4. Le signal brut du capteur, tel que visible sur l'oscilloscope 750, est dénommé A-Scan par les hommes du métier. Il comprend des échos selon le schéma défini par la figure 6. Il est souhaitable de passer à une imagerie des imperfections du tube, à l'aide des signaux d'ultrasons. On décrira maintenant l'obtention d'une image. En pratique, une image est obtenue en considérant plusieurs explorations successives du tube par un capteur Px, sous des angles successifs qui couvrent sensiblement une section droite du tube. Il est possible de le faire avec des tirs successifs par un seul capteur, en utilisant la rotation relative tube/capteur. On se place ici, à titre d'exemple non limitatif, dans le cas d'une installation du genre dit à tête tournante. Sur la figure 8A, on considère un capteur Px, qui peut être de l'un des types P1, Pli, P12, P21 et P22 précité. Dans l'exemple illustré, ce capteur Px comprend en fait n capteurs élémentaires Px-1, ..., Px-i, ..., Px-n, qui sont alignés selon l'axe longitudinal du tube, et qui font l'objet d'un tir ultrasonore au même moment. Sur la figure 8A, ce qui est entre les capteurs élémentaires et le graphe 3D de sortie 769 peut être considéré comme un convertisseur. Le signal Ascan du premier capteur élémentaire Px-1 est appliqué à un amplificateur 73-1, suivi de deux voies parallèles : celle du sélecteur 763-1A et celle du sélecteur 763- 1B. Chaque sélecteur 763-1A comprend deux sorties respectivement d'amplitude maximale et de temps de vol. La sortie d'amplitude maximale est reliée à un numériseur de ligne 765-lAa. La sortie de temps de vol est reliée à un numériseur de ligne 765-lAt. La sortie de numériseur de ligne 765-lAa d'amplitude maximale est reliée à un stockage tampon de données 768-Aa recueillant les données provenant des numériseurs de ligne d'amplitude maximale 765-iAa d'indice i allant de 1 à n. La sortie de numériseur de ligne 765-1At de temps de vol est reliée à un stockage tampon de données 768-At recueillant les données provenant des numériseurs de ligne de temps de vol 765-iAt d'indice i allant de 1 à n. La sortie de numériseur de ligne 765-1Ba d'amplitude maximale est reliée à un stockage tampon de données 768-Ba recueillant les données provenant des numériseurs de ligne d'amplitude maximale 765-iBa d'indice i allant de 1 à n. La sortie de numériseur de ligne 765-1Bt de temps de vol est reliée à un stockage tampon de données 768-Bt recueillant les données provenant des numériseurs de ligne de temps de vol 765-iBt d'indice i allant de 1 à n. Sur la base des informations obtenues lors du passage du tube étalon, l'opérateur peut saisir dans les stockages tampons 768-Aa et 768-At l'information T_lA correspondant à une indication de position et de largeur temporelle, qui lui désigne, en fonction de la géométrie connue du tube, les instants où il va trouver un � écho de peau interne �, relatif à l'intérieur du tube par exemple le premier écho Intl de la figure 6. La figure 6A montre plus nettement la fenêtre temporelle � Int � correspondante, autour de Tphilnter. De même, sur la base des informations obtenues lors du passage du tube étalon, l'opérateur peut saisir dans les stockages tampons 768-Ba et 768-Bt l'information Ti B correspondant à une indication de position et de largeur temporelle, qui lui désigne, en fonction de la géométrie connue du tube, les instants où il va trouver un � écho de peau externe �, relatif à l'extérieur du tube, par exemple le premier écho Exil de la figure 6. La figure 6A montre plus nettement la fenêtre temporelle � Ext � correspondante, autour de TphiExter. Le schéma se répète pour les autres capteurs Px-2, ... Px-i, ... Px-n. Ainsi, chaque sélecteur temporel 763 définit des fenêtres temporelles compte tenu de l'instant d'émission des ultrasons, et des intervalles de temps prédéterminables où l'on peut attendre des échos concernant ce sélecteur. L'illustration des figures 6 montre comment l'on peut définir les intervalles de temps intéressants, compte tenu de l'angle d'attaque du faisceau d'ultrasons sur le tube, ainsi que du diamètre (interne ou externe) et de l'épaisseur du tube. Un intervalle de temps donné correspond à un écho donné en un point donné du tube, pour une position relative donnée entre le tube et le capteur. Pour simplifier, on admet ici que les instants de tirs sont synchronisés sur la rotation relative tube/capteurs, de sorte qu'un capteur élémentaire travaille toujours sur la même génératrice longitudinale du tube. La sortie de son sélecteur fournit donc une suite espacée d'échantillons analogiques de signal, qui correspondent chacun à l'amplitude d'un écho attendu sur une paroi du tube. Ces échantillons du capteur Px-1 (par exemple) sont numérisés en 765. Le synchronisme avec l'émission peut être assuré par une liaison (non représentée) avec l'émetteur 70, ou avec son déclencheur, le circuit de synchronisation 753, ou sa base de temps 752 (figure 8). L'affichage 750 peut être maintenu, si on le souhaite. Le système peut fonctionner sur tube tournant à vitesse sensiblement constante. Dans ce cas, la vitesse angulaire et l'avance du tube peuvent être mesurées à l'aide d'un codeur angulaire précis, par exemple le modèle RS0550168 fourni par la société Hengstler, et d'un vélocimètre laser, par exemple le modèle LSV 065 fourni par la société Polytec. Le tube peut aussi ne pas tourner, tandis que c'est le système de capteurs qui tourne. Dans ce cas, le vélocimètre laser suffit pour mesurer l'avance du tube, tandis que la vitesse de rotation des capteurs peut être connue au moyen d'un codeur angulaire. Pour un tir donné, l'ensemble des capteurs Px-1 à Px-n fournit une ligne d'une image, 10 qui correspond à une section droite du tube. Dans l'autre dimension de l'image, un capteur élémentaire donné fournit une ligne qui correspond à une génératrice du tube. Les numériseurs 765-lAa, 765-2Aa, ...,765-iAa, ..., 765-nAa et 765-lAt, 765-2At, ...,765-iAt, ..., 765-nAt permettent de remplir une image � interne �, relative à la peau 15 interne du tube. Les numériseurs 765-1Ba, 765-2Ba, ...,765-iBa, ..., 765-nBa et 765-lBt, 765-2Bt, .,765-iBt, ..., 765-nBt permettent de remplir une image � externe �, relative à la peau externe du tube, avec Tvol max le temps de vol de l'écho d'amplitude maximale...DTD: 20 Le graphe 3D parallélépipédique mémorisé en 769 vaut pour le capteur ou groupe de capteurs Px considéré. Chaque point de cette image correspond, transposé en nuances de gris, a une valeur d'amplitude de l'écho dû à la réflexion du signal ultrasons sur une imperfection éventuelle de la zone du tube considérée. Cette valeur peut également représenter le rapport entre l'amplitude maximale du signal ultrasonore capté sur le tube 25 en cours de test et l'amplitude maximale du signal ultrasonore obtenu avec un � défaut étalon de référence � artificiel, tel que défini plus haut. Le graphe 3D parallélépipédique est une représentation du Bscan 3D préparatoire numérisé en 769 - préparatoire en ce sens qu'il sert de base à la génération du Bscan 3D tube. La forme du graphe 3D est généralement distincte de la forme du produit examiné, notamment pour les tubes. 30 Les informations du graphe 3D parallélépipédique peuvent comprendre l'ensemble des couples (temps de vol, amplitude) de la courbe AScan sur une durée de numérisation déterminée. Les graphes 3D parallélépipédiques numérisés en 769 comprennent les graphes 3D parallélépipédiques 891 construits à partir des données provenant d'un groupe de capteurs P11 et les graphes parallélépipédiques 892 construits à partir des données provenant d'un groupe de capteurs P12 et respectivement P21 et P22 tel que représenté sur la figure 11. Cette image correspond maintenant à une zone du tube, obtenue par la réunion des zones sensiblement annulaires du tube qui correspondent à chacune des lignes numérisées. En fait, il s'agit de zones annulaires ou hélicoïdales si le faisceau d'ultrasons est appliqué sensiblement perpendiculairement à l'axe du tube. On sait qu'il en est différemment selon le mouvement relatif tube/capteur. Les zones sont alors plutôt elliptiques et, en fait, gauches ou "tordues" dans l'espace. Dans la présente description, l'expression "zones annulaires" couvre ces différentes possibilités. Il est à noter que pour obtenir une reconstitution complète du graphe 3D, l'information supplémentaire de positionnement du capteur par rapport au tube est nécessaire. Elle est disponible en entrée séparée 740. Cette information vient d'un codeur ou d'un ensemble de lasers permettant de mesurer le positionnement spatial. Comme le tube peut être assimilé à un cylindre sans épaisseur, l'information de position peut être réduite à deux dimensions. On comprend que la mise en oeuvre de l'invention sur un banc existant de contrôle par ultrasons ("UT bench") implique : - l'accessibilité aux données brutes de contrôle par ultrasons ("UT Raw Data"), qui se fait par exemple à l'aide d'une carte d'acquisition, comme le modèle NI 6024 série E ou NI 6251 série M de la société National Instrument, ou par accès direct aux données numériques d'une électronique de contrôle du banc, - la disponibilité d'une information en ligne sur la vitesse de rotation (du tube ou capteur) ou la position angulaire relative du tube par rapport au capteur, et - la disponibilité d'une information en ligne sur la vitesse d'avance du tube ou la position linéaire relative du capteur projetée sur l'axe. Le schéma de la figure 8A peut être appliqué: - en parallèle à un capteur de type P 1 l et à un capteur de type P12, qui observent la même zone du tube selon deux directions différentes. Chaque capteur va permettre d'obtenir une image interne et une image externe. Ensuite, l'une des images pourra être choisie en fonction d'une commande notée � Int/Ext �. -en parallèle à un capteur de type P21 et à un capteur de type P22, qui, là aussi, vont chacun permettre d'obtenir une image interne et une image externe. Le schéma de la figure 8A peut aussi être appliqué à un capteur de type Pl, auquel cas on prévoit trois voies parallèles derrière chaque amplificateur (au moins virtuellement). L'une des voies fonctionne sur un créneau temporel répétitif positionné comme indiqué en � Volum. � sur la figure 6B. Cette voie permet un contrôle d'imperfections en volume, c'est-à-dire dans l'épaisseur du tube. Les deux autres voies peuvent fonctionner respectivement sur les créneaux temporels répétitifs positionnés comme indiqué en � WphiExterO � et en � WphiInterl � sur la figure 6B. Ces deux autres voies permettent une mesure de l'épaisseur du tube. La distinction entre les 3 voies est purement fonctionnelle (virtuelle). En effet, lesdites deux autres voies peuvent être physiquement la même, dans laquelle on discrimine les instants ou créneaux � WphiExterO � et � WphiInterl �. On peut aussi utiliser une seule voie physique, dans laquelle on discrimine les instants ou créneaux � WphiExterO �, � Volum. � et � Wphilnterl �. Il est représentatif de décrire plus en détail le cas d'un capteur de type P11 avec un capteur de type P12. C'est ce que l'on fera maintenant. Il est rappelé que ces deux groupes de capteurs P11 et P12 servent à la détection d'imperfections longitudinales sur les tubes. Le contrôle ultrasonore est réalisé avec des tirs d'ultrasons (US) dans deux directions privilégiées ("dock wise" - "counter dock wise"): - Un capteur ou groupe de capteurs P11 fournit une image ultrasonore du tube dans une direction de travail ("dock wise"). - Un deuxième capteur ou groupe de capteurs P12 fournit une image ultrasonore du même tube dans une autre direction de travail ("counter dock wise"). Ainsi, les imperfections longitudinales sont avantageusement détectées avec 2 capteurs ou groupes de capteurs dont les axes de faisceau sont inclinés symétriquement par rapport à un plan perpendiculaire à l'axe du tube. L'inclinaison est par exemple d'environ +/- 17°. Ceci fournit un exemple d'application du système à deux capteurs, ou deux groupes de capteurs, tel que mentionné plus haut. Dans le mode de réalisation de la figure 8B, chaque fenêtre de numérisation 782 issue d'un amplificateur 781 peut être caractérisée par un début, une durée et une fréquence de numérisation qui définissent un nombre n de points du signal AScan pris en considération. Chaque fenêtre de numérisation 782 fournit alors un nombre n de couples d'information (Amplitude, Temps de vol), pour chaque tir ultrasonore. Le Buffer/Multiplexeur 788 remet l'ensemble des données ainsi recueillies dans le graphe 3D parallélépipédique 769 en tenant compte des positions respectives des capteurs au moment où le signal a été reçu, ceci à la fois grâce à la connaissance de la configuration géométrique des capteurs les uns par rapport aux autres, et grâce à l'information de positionnement tube/capteur au moment du tir ultrasonore 740. Il est maintenant fait référence à la figure 9. Pour le premier sens de contrôle (onglet � sens 1 � sélectionné), les images 903 et 904 sont des vues en coupe (respectivement transversale et longitudinale) du Bscan 3D Tube, 3D avec la géométrie du tube, tel que décrit plus loin, provenant des capteurs P11. Le positionnement de ces coupes est fixé grâce aux paramètres � coupe transversale à (mm) � et � coupe longitudinale à (degrés) �. Les images 905 (interne) et 906 (externe) sont des CScans, tels que définis précédemment, l'image 905 (respectivement 906) se concentrant sur une zone temporelle du Ascan où les imperfections en peau interne (respectivement externe) sont censées être détectées. Les informations nécessaires à la reconstruction des images 905 et 906 proviennent du BScan 3D parallélépipédique 891 de la figure 11. L'image 901 est une représentation 3D en transparence du Bscan 3D Tube d'une portion du produit à contrôler, portion dans laquelle sont identifiées des zones potentiellement intéressantes, telles que décrites plus loin. Les mêmes images 903 bis, 904 bis, 905bis, 906 bis et 902 sont reconstituées pour le deuxième sens de contrôle (onglet � sens 2 � activé), voir figure 9A. Nous rappelons ici que la description précédente concerne la détection de défauts à orientation longitudinale. La même démarche s'applique pour la recherche de défauts transversaux (avec les groupes de capteurs P21 et P22). Il est maintenant fait référence à la figure 11. Les blocs d'images 901 et 902 sont obtenues à partir des graphes 3D parallélépipédiques 891 et 892 au moyen du bloc de transformation 930 tel que détaillé en figure 11A. Le bloc convertisseur 891 de la figure 11 correspond au montage de la figure 8A, appliqué au capteur P11. De même, le bloc convertisseur 892 correspond lui aussi au montage de la figure 8A, mais appliqué au capteur P12. Les blocs convertisseurs 891 et 892 utilisent les données de contexte tube/capteurs du bloc 740. Ces données sont relatives aux caractéristiques du tube en cours d'examen et des capteurs en cours d'utilisation. Le bloc de transformation 930 est disposé en aval des graphes 3 D parallélépipédiques 891 et 892 et peut présenter la structure illustrée sur la figure 11A. Le bloc de transformation 930 effectue un calcul temporel du parcours de la propagation des ondes dans le tube en prenant en compte la conversion de modes au moment de l'impact d'une onde ultrasonore sur un défaut. A l'impact, une onde transversale peut se transformer en onde longitudinale et vice versa. Le bloc de transformation 930 peut estimer la propagation de l'énergie du faisceau acoustique à partir de calculs de coefficients de transmission et de réflexion. Une analyse du spectre fréquentiel de l'Ascan peut être réalisée. Le bloc de transformation 930 peut comprendre une base de données 939 d'essais réels ou simulés permettant une comparaison avec les graphes 3D reçus. Le bloc de transformation 930 peur reconstituer l'image Bscan 3D avec la géométrie du tube. Comme illustré sur la figure 11A, le bloc de transformation 930 comprend deux blocs 931 et 932 d'élimination de zones de Bscans 3D non utiles à partir d'un graphe 3D, le bloc 931 traitant les données des images 3D 891 et le bloc 932 traitant les données des images 3D 892, deux blocs 933 et 934 de filtrage par application d'une fenêtre temporelle simulée, en aval respectivement des blocs 931 et 932, un bloc de simulation théorique 935, un bloc de calcul de tolérance 937 alimentant un bloc d'algorithme inverse 936, le bloc 936 fournissant les images 901 et 902 définies précédemment. L'élimination par les blocs 931 et 932 permet de réduire la quantité d'information traitée en conservant les zones potentiellement intéressantes à représenter de façon tridimensionnelle. Le filtrage peut être effectué en long à partir d'un Cscan. La longueur choisie peut être supérieure à la longueur d'une zone d'amplitude supérieure à un seuil. On peut ensuite traiter les Bscans 3D parallélépipédiques incluant une zone à imperfection potentielle. Le filtrage par les blocs 933 et 934 peut être effectué en bornant la fenêtre temporelle par les échos d'interface et de fond. Ces blocs de filtrage peuvent également borner la zone angulaire du tube potentiellement intéressante et si nécessaire décaler ces zones afin de cerner et reconstituer complètement la zone potentiellement intéressante. Les images issues des blocs 933 et 934 sont des Bscan 3D réduits. Le bloc de simulation théorique 935 peut comprendre une base de données de simulations, par exemple d'Ascans ou de Bscans 3D en fonction des types et position des défauts. La base de données peut comprendre des résultats simulés et/ou des résultats d'essais sur des défauts naturels et/ou artificiels. Le bloc d'algorithme inverse 936 peut comparer des Ascans ou Bscans 3D théoriques provenant du bloc de simulation théorique 935 et des Ascans ou Bscans 3D obtenus lors de l'inspection afin de déterminer l'Ascan ou le Bscan 3D théorique le plus proche et, par conséquent, le(s) défaut(s) le(s) plus probable(s). A titre d'exemple, le bloc d'algorithme inverse 936 compare un Ascan expérimental filtré correspondant à une position en longueur et à une position angulaire avec les Ascans théoriques sur cette même position en longueur et en développée. A titre d'autre exemple, le bloc d'algorithme inverse 936 compare un Bscan 3D issu d'un Bscan 3D réduit correspondant à une position en longueur avec les Bscans 3D théoriques sur cette même position en longueur. Les deux comparaisons peuvent être effectuées. Le meilleur ensemble de représentations théoriques des échos est alors l'ensemble présentant le plus petit écart par rapport aux données expérimentales. Après le bloc de transformation 930 sont illustrés les filtres 921 et 922, voir figure 11, qui permettent notamment de réaliser des extraits des images, et de leurs données de préparation, en tant que données d'entrée réunies par le bloc combinateur 960 pour le traitement neuronal ou expert 970. Dans le mode de réalisation décrit, le filtre 921 possède : une sortie de signal Zcur désignant une zone de travail dans l'image. Cette sortie est utilisée par une fonction d'extraction 951 qui réalise en conséquence un extrait de l'image (Cscan) pour la zone Zcur, et un accès à la préparation d'image 891 pour y obtenir des informations mémorisées (dites Ascan) relatives à la même zone Zcur. L'ensemble de ces données est transmis par la fonction d'extraction 951 au combinateur 960, comme entrées du traitement neuronal ou expert 970, une sortie fournissant des informations obtenues par filtrage, certaines au moins 25 relatives à la zone Zcur, qu'il transmet comme entrée du traitement neuronal ou expert, en option (trait tireté) des sorties de données complémentaires filtrées vers une mémoire 990. 30 Il en est de même pour le filtre 922, avec la fonction d'extraction 952, pour la même zone courante Zcur. Le système neuronal 970 alimente une logique de décision et d'alarme 992, qui pilote un automatisme de tri et marquage 994. Il peut être prévu une interface d'interprétation 996 par un opérateur, laquelle peut présenter tout ou partie des données contenues dans la mémoire 990, en rapport avec la portion de tube en cours d'examen. Les données contenues dans la mémoire 990 proviennent des filtres 921 et 922. En plus de sa prédiction (origine, type et gravité de l'indication) le système neuronal 970 fournit une évaluation de la confiance que l'on peut apporter à cette prédiction. Cette information est accessible aux opérateurs qui disposent également de données complémentaires plus qualitatives telles que l'historique de la commande en cours ou les problèmes qui ont eu lieu au cours de l'élaboration du produit. L'opérateur, ou un spécialiste peut alors intervenir pour pondérer les prédictions en conséquence. Ici, la figure 11 traite des informations provenant au minimum de deux groupes de capteurs assurant la même fonction ou destinés au même type de contrôle (les 2 groupes P 11 et P12 ou les 2 groupes P21 et P22). Le même schéma peut servir à traiter les informations provenant d'un plus grand nombre de groupes de capteurs destinés à des contrôles de type différent. Le nombre d'images traitées en même temps est augmenté d'autant. La fonction primaire des filtres 921 et 922 est de déterminer des zones d'imperfections dans les images 901 et 902. De façon générale, le filtrage est agencé pour repérer les zones à analyser, et y distinguer les imperfections des autres indications. Le filtrage travaille sur deux portions homologues de deux images. Les deux filtres peuvent travailler conjointement. Par balayage de l'image numérique, on localise d'abord les endroits de l'image où existent des imperfections potentielles. A cet effet, il est possible d'appliquer un seuil fixe établi par étalonnage.30 On peut utiliser un seuil qui s'adapte au niveau de bruit en cours dans l'image. La méthode est fondée sur la théorie de la détection d'un signal dans un bruit blanc, qui peut se fonder sur deux hypothèses : Hypothèse HO : mesure = bruit blanc de moyenne m_b et d'écart-type std_b 5 Hypothèse Hl : mesure = signal + bruit blanc On procède à des tests statistiques, qui permettent de déterminer si l'on est dans le cadre de l'hypothèse HO, ou de l'hypothèse Hl. Ces calculs statistiques sont effectués en temps réel sur n points glissants de l'image correspondant à des tirs consécutifs, le nombre n pouvant être déterminé par l'apprentissage. 10 Selon cette méthode (cas dit "additif gaussien"), on peut par exemple utiliser le critère de Neyman-Pearson pour déterminer un seuil de détection selon une probabilité de fausse-alarme (pfa) donnée. Cela s'exprime par la formule [21] annexée. On utilise la fonction cumulative gaussienne, nommée en général Q (ou encore la fonction d'erreur 15 erf), qu'il faut inverser pour obtenir le seuil, selon la formule [22] annexée. En pratique, on constate fréquemment la présence de bruit de fond pouvant avoir plusieurs origines (par exemple: présence d'eau à l'intérieur du tube, ronflement électrique, phénomènes acoustiques dus à la structure de la matière du produit 20 contrôlé). L'usage d'un seuil variable évite les fausses alarmes qui se produisent si l'on applique un seuil fixe. Parmi les autres fausses indications susceptibles d'apparaître, les parasites se 25 manifestent par des pics très brefs dans le signal d'ultrasons. Ces parasites peuvent être écartés par des algorithmes simples que l'on peut appeler algorithmes de comptage cumulatif ou encore intégrateurs (exemple: "n coups avant alarme" ou "double seuil"). La demanderesse a encore considéré la "spire", qui est le trajet suivi par le capteur le 30 long de la surface cylindrique à laquelle le tube est assimilé. Un filtrage peut être effectué le long de chaque spire pour encore réduire le taux de fausses alarmes. On utilise à cet effet par exemple un filtre de Butterworth, et/ou une transformée de Fourier discrète, telle qu'une transformée de Fourier rapide. Cette technique est appliquée à chaque ligne numérique. Le même genre d'algorithme peut être appliqué dans le sens de la longueur du tube. 5 Ainsi, des imperfections potentielles sont localisées. Lorsqu'une imperfection est repérée, sa position correspond à la position analysée dans les images de la figure 9 (par exemple), avec une image 3D, une coupe transversale et une coupe axiale. Les indications de position radiale/épaisseur (ou plus simplement de situation interne, 10 externe ou en masse de l'imperfection) peuvent être représentées comme des attributs des points de l'image. On aura ainsi: deux images 2D représentant les imperfections éventuelles en peau externe du tube, deux images 2D représentant les imperfections éventuelles en peau interne du 15 tube, et une image 2D représentant les imperfections éventuelles dans l'épaisseur du tube. On considère maintenant les imperfections � confirmées �, après élimination des 20 parasites et des fausses alarmes, notamment. Pour la suite, la Demanderesse a choisi actuellement de travailler sur une zone d'image de taille fixe. Il faut donc cadrer cette zone sur les données d'existence d'imperfection que l'on vient d'obtenir. Autrement dit, il y a lieu de positionner les points repérés comme étant supérieurs au seuil pour déterminer la zone complète autour d'une imperfection. C'est un besoin par exemple, si l'on souhaite déterminer l'obliquité d'une imperfection. 30 L'algorithme s'articule autour de différentes étapes : - détection de contour (gradient de Roberts par exemple), - dilatation (rassemblement des contours proches), 25 - érosion, puis fermeture, ce qui permet de déterminer un masque autour des imperfections, - une dernière étape d'entourage permet de localiser complètement l'imperfection. Pour chaque imperfection, on obtient ainsi les coordonnées de la zone d'image correspondante, qui seront utiles pour l'analyse par réseau de neurones qui intervient ensuite. La figure 12 illustre ce traitement des zones d'image sous la forme d'un diagramme de 10 flux. En début d'images (801), on a de zéro à p zones d'image à traiter, comme présentant une imperfection confirmée. L'opération 803 suppose qu'il existe au moins une première zone, laquelle sert de zone courante à traiter Zcur en 805. Pour cette zone 15 Zcur: l'opération 807 extrait sélectivement les données des images 901 et 902 qui correspondent à cette zone (définie par ses coordonnées dans l'image). l'opération 809 extrait sélectivement des données qui sont intervenues dans la préparation des images 901 et 902, et qui correspondent à la zone Zcur. Des 20 exemples de ces données seront donnés ci-après. l'opération 811 réalise le traitement neuronal ou expert proprement dit, sur lequel on reviendra. Les résultats obtenus pour la zone Zcur sont mémorisés sélectivement en 813, en correspondance d'une désignation de la zone Zcur. 25 Le test 820 recherche s'il existe une autre zone à traiter dans l'image, auquel cas on recommence en 805 avec cette autre zone comme indiqué en 821 ; sinon le traitement de la ou des images en cours est terminé (822). Dans le cas du traitement du capteur P1, il n'y a qu'une seule image, ce qui change le 30 nombre de paramètres d'entrée. A part cela, le traitement peut être généralement le même. Après la détermination de chaque zone d'intérêt Zcur, le filtrage peut comporter d'autres fonctions. Pour ces autres fonctions, la figure 13 illustre de façon schématique l'interaction entre le filtrage et la suite des opérations illustrées sur la figure 11. La figure 13 est semblable à la figure 11, mais seulement pour l'image 901. Elle fait apparaître : les éléments de contexte tube- capteurs du bloc 740, l'extracteur 951 qui trouve les données pour la zone Zcur, dans l'image 901 et sa préparation 891, un bloc interne/externe 7410, indiquant si l'imperfection dans la zone Zcur considérée est située en peau interne ou en peau externe. Ce que le filtrage ajoute aux données de base, est défini en plus de détails à savoir, pour chaque zone Zcur (bloc 805), comme l'indique le contenu du cadre en trait tireté : 15 - une recherche de l'angle d'obliquité en 941, une indication de longueur d'imperfection 942, Il peut s'y ajouter encore, notamment : - une indication d'alignement en C-Scan, en 945, et - en 946, une indication sur l'existence d'autres imperfections dans la même section 20 droite du tube. Dans le mode de réalisation décrit, les données telles que 945 et 946 vont vers la mémoire 990. Les autres données vont vers les réseaux neuronaux ou systèmes experts 970. Ceux-ci sont ici séparés en deux fonctions, comme on le verra maintenant. 25 Une imperfection dans un tube peut être définie par sa position, son type, et sa gravité souvent assimilée à sa profondeur. Dans le mode de réalisation décrit, le type et le degré de profondeur d'une imperfection de tube sont déterminés séparément à l'aide de deux processus neuronaux de même structure générale, que l'on détaillera maintenant sur un 30 exemple. Le cas du type de l'imperfection est traité selon la figure 14, tandis que le cas de la gravité est traité selon la figure 15. Les types peuvent être définis par exemple comme illustré sur les figures 10A à 10D. Ces figures illustrent quatre types, constituant un choix simplificateur par rapport à la liste des imperfections fournies par l'API et pouvant être produits par les processus d'élaboration du tube. Les intitulés en français et en anglais sont ceux utilisés par les hommes du métier pour désigner la nature des imperfections. On observe que les imperfections des types 1 et 3 sont droites, celles des figures 2 et 4 arquées (à "chord"). Une correspondance entre les imperfections réelles et les quatre types ci-dessus peut être définie comme suit : Nom en français Nom en anglais Affectation Entaille Notch TYPE 1 Tapure Crack TYPE 1 Paille/repliure perpendiculaire ou Seam (perpendicular) TYPE 1 droite (laminage) Paille/repliure (laminage) Seam (arcuate), "overlap" TYPE 2 Gravelure Sliver TYPE 3 Origine billette Rolled-in-slug TYPE 4 Rayure Gouge TYPE 4 Inclusion Inclusion TYPE 4 Manque de matière ("défourni") Bore-slug TYPE 4 Chevauchement/recouvrement/repliure Lap TYPE 4 Ici, les figures 14 et 15 utilisent toutes deux des circuits neuronaux à trois neurones 15 intermédiaires (ou "neurones cachés"), notés NC 121 à NC 123 pour la figure 14, et NC 141 à NC 143 pour la figure 15. Les figures 14 et 15 ont en commun un certain nombre d'entrées. Pour tenter de faciliter la compréhension, les entrées sont illustrées par des types de traits différents. Les traits doubles indiquent que les entrées sont multiples, c'est-à-dire répétées pour chaque point de la zone Zcur. Tout d'abord, en 7410, il provient, conformément à l'état considéré des sélecteurs 763 concernés, une information indiquant s'il s'agit de traiter d'une imperfection située en peau interne ou en peau externe de la paroi du tube. Cette information peut également être obtenue sur le BScan 3D. La seconde catégorie de grandeurs d'entrée communes comprend les grandeurs de 10 contexte, qui viennent du bloc 740 (figure 13) : - en 7401, WT/OD, qui est le rapport de l'épaisseur de paroi au diamètre du tube, - en 7402, Freq, qui est la fréquence de travail des sondes à ultrasons, - en 7403, ProbDiam, qui est le diamètre utile des sondes à ultrasons. 15 La troisième catégorie de grandeurs communes comprend des quantités issues du filtrage qui peuvent être considérées comme communes aux deux capteurs 921 et 922 (ou plus). On fait par exemple la moyenne des résultats sur les deux capteurs, ou bien on prend le résultat le plus représentatif (maximum/minimum, suivant le cas) Ces quantités sont les grandeurs en 9201, l'obliquité du défaut, et en 9202 sa longueur. Ces 20 deux grandeurs sont aisément repérables sur les deux images de la figure 9, qui ont une symétrie miroir. Il est maintenant fait référence à la figure 14 seulement. La catégorie suivante de grandeurs comprend des grandeurs de mesures distinctes pour chacun des deux capteurs 25 (ou groupe de capteurs), et pour chacune des zones Zcur, ce qui est reflété sur le dessin par l'usage d'un trait double. Pour un premier capteur, on a : - en 9511, K], qui est le rapport entre l'amplitude maximale du signal ultrasonore 30 rencontré dans la zone Zcur et sur l'image 901, par rapport à l'amplitude maximale du � défaut étalon de référence � précité. En fait, dans l'exemple, l'amplitude en chaque pixel de l'image 901 est définie par ce rapport ; K1 est alors simplement le maximum d'amplitude rencontré dans la zone Zcur de l'image 901 ; on note Pmaxl le point de la zone Zcur où ce maximum est rencontré. - en 9512, QBE1 qui est une grandeur du C-Scan dite QuantBumpsEchodyn, représentant le nombre de maxima locaux rencontrés dans la zone Zcur de l'image 901 à proximité du point Pmaxl d'amplitude maximale. Ce nombre QBE1 est limité aux maximas locaux rencontrés au voisinage de Pmaxl, de part et d'autre, mais sans que l'amplitude du signal ne soit redescendue en dessous d'un niveau correspondant au bruit de fond. QBE1 va généralement prendre soit la valeur 1, soit la valeur 2. Ces deux grandeurs proviennent de l'image 901, via l'extracteur 951, ce que reflète la notation 951(901) sur le dessin. Il s'y ajoute : - en 9518, RT1 qui est une grandeur représentant le temps de montée de l'écho dans le signal ultrasonore natif dit A-Scan, (il s'agit de l'écart entre le moment où le signal est maximal et le dernier moment antérieur où le signal est au niveau du bruit de fond, exprimé couramment en micro-secondes). Cette grandeur RT1 a été antérieurement mesurée en sortie de l'amplificateur 73 concerné (figure 8A) ; elle a été stockée, par exemple en 891, en correspondance du point du tube qu'elle concerne. C'est ainsi qu'elle peut être récupérée sélectivement par l'extracteur 951. La grandeur RT1 peut maintenant être directement mesurée par l'opérateur sur l'image 903 de la figure 9, ou encore sur le BScan 3D parallélépipédique. Pour le second capteur, on a : - en 9521, K2, qui est défini comme K1, mais pour l'image 902 au lieu de l'image 901. Dans l'exemple, K2 est simplement le maximum d'amplitude rencontré dans la zone Zcur de l'image 902 ; on note Pmax2 le point de la zone Zcur où ce maximum est rencontré. - en 9522, QBE2 est défini comme QBE1, mais dans l'image 902 au lieu de l'image 901, et au voisinage de Pmax2. Là aussi, QBE2 va généralement prendre soit la valeur 1, soit la valeur 2. Ces deux grandeurs proviennent de l'image 902, via l'extracteur 952. Il s'y ajoute : - en 9528, RT2 qui est une grandeur représentant le temps de montée de l'écho dans le signal natif dit A-Scan. Comme précédemment, cette grandeur RT2 a été antérieurement mesurée en sortie de l'amplificateur 73 concerné (figure 8A) ; elle a été stockée, par exemple en 892, en correspondance du point du tube qu'elle concerne. C'est ainsi qu'elle peut être récupérée sélectivement par l'extracteur 952. La grandeur RT2 peut maintenant être directement mesurée par l'opérateur sur l'image 903A de la figure 9, ou encore sur le BScan 3D parallélépipédique. La dernière entrée 958 du réseau neuronal est une valeur constante, notée ConstantA, 10 qui représente une constante déterminée lors du calage du modèle et résultant de l'apprentissage. La sortie 998 de la figure 14 est une grandeur indicative du type de l'imperfection et son inclinaison moyenne (définie en fonction du type). 15 Le cas du degré de profondeur (ou gravité) de l'imperfection est traité selon la figure 15. Les entrées sont les mêmes que pour la figure 14, sauf : pour le premier capteur, le bloc 9512 est remplacé par un bloc 9513, qui traite une grandeur EW_l, ou EchodynWidth, qui est la largeur à mi-hauteur (50%) de la 20 forme d'onde échodynamique, pour ce premier capteur. Cette grandeur EW_l est tirée du Cscan. de même, pour le second capteur, le bloc 9522 est remplacé par un bloc 9523, qui traite la grandeur EW 2, ou EchodynWidth, qui est la largeur à mi-hauteur (50%) de la forme d'onde échodynamique, pour ce second capteur. 25 en 959, la constante, notée maintenant ConstantB, est différente. la sortie 999 est une indication de gravité d'imperfection, notée DD. On note que, dans les deux cas (figures 14 et 15), un circuit neuronal 970 donné traite un extrait d'image 951 pour l'un des groupes de capteurs ultrasonores, ainsi qu'un extrait 30 d'image 952 correspondant à la même zone, mais provenant d'un autre groupe de capteurs. La Demanderesse a observé qu'il était possible d'obtenir des résultats très satisfaisants, sous réserve d'un ajustement convenable des paramètres d'un système expert, par exemple des circuits neuronaux, et éventuellement de leur nombre, pour optimiser la prédiction. De plus, la demanderesse a constaté que par combinaison des informations recueillies par les différents réseaux de neurones il était possible d'affiner encore la prédiction. Globalement, les paramètres d'entrée du réseau de neurones ou du système expert sont 10 alors des caractéristiques des deux images 3D (rapport de l'amplitude max par rapport à l'amplitude de l'étalon, largeur d'écho, orientation de l'écho représentative de l'obliquité de l'imperfection...) et du contrôle (capteur, dimensions du tube...). Les paramètres de sorties sont les caractéristiques de l'imperfection (profondeur, inclinaison/ type). La décision et/ou alarme 992 peut se prendre automatiquement à 15 l'aide de critères de décision choisis, à base de seuils, assortis d'une marge de sécurité selon les besoins. Pour définir ces seuils, on peut s'aider des résultats de l'apprentissage. Il est maintenant fait référence à la Figure 16, qui est un modèle du circuit neuronal élémentaire des figures 14 ou 15, pour deux capteurs. 20 Ce modèle comprend un niveau ou couche d'entrée IL ("Input Layer"), qui regroupe tous les paramètres d'entrée (souvent appelés "neurones d'entrée"). Pour ne pas surcharger la figure, sont seuls représentés trois neurones El à E3, plus une constante, qui peut être considérée elle aussi comme un neurone E0. Cette constante est le plus 25 souvent appelée � biais �. En pratique, les neurones d'entrée sont plus nombreux, conformément à la figure 14 ou à la figure 15, selon le cas. Ensuite est prévu au moins un niveau ou couche HL ("Hidden Layer" ou "couche cachée"), qui comprend k neurones (dont seulement 2 sont représentés pour ne pas 30 surcharger le dessin).5 Enfin vient le neurone de sortie S1, qui fournit la décision, sous la forme d'une valeur représentative de l'importance d'une imperfection du tube, par exemple une imperfection longitudinale. Cette sortie correspond au bloc 998 sur la figure 14 et 999 sur la figure 15. On notera que le "neurone" - constante EO intervient pour pondérer non seulement la ou les couches cachées HL, mais aussi le neurone de sortie (couche OL ou "Output Layer"). 10 Le comportement général d'un circuit neuronal tel qu'utilisé ici est donné par la formule [Il] de l'Annexe 1, où wij est le poids affecté au signal Xi présent à l'entrée du neurone J• Dans le circuit prévu ici, un neurone élémentaire se comporte selon la formule [12], 15 comme cela est schématisé sur la figure 17. La sortie S1 de la figure 16 fournit une valeur estimée qui correspond à la formule [13] de l'annexe 1. Par apprentissage, la Demanderesse a ajusté les neurones cachés et leurs poids de sorte 20 que la fonction f soit une fonction non-linéaire, continue, dérivable et bornée. L'exemple actuellement préféré est la fonction arc-tangente. On sait qu'un réseau neuronal détermine ses coefficients w;i communément appelés synapses par apprentissage. Cet apprentissage doit faire intervenir typiquement 3 à 10 25 fois plus d'exemples qu'il y a de poids à calculer, tout en couvrant correctement la plage des conditions de travail désirées. Partant d'exemples Ep (p = 1 à m), on détermine pour chaque exemple l'écart Dp entre la valeur Sp donnée par le circuit neuronal et la valeur réelle Rp mesurée ou définie 30 expérimentalement. C'est ce que rappelle la formule [14].5 La qualité de fonctionnement du circuit neuronal est définie par une grandeur globale d'écart Cg, dite "coût". Elle peut s'exprimer par exemple selon la formule [15], comme une grandeur globale d'écart quadratique pondéré. L'apprentissage pose différents problèmes dans un cas comme celui du contrôle des imperfections dans les tubes, notamment du fait qu'il s'agit de techniques lourdes, comme déjà indiqué. La demanderesse a tout d'abord conduit un premier apprentissage sur simulation. On peut utiliser à cet effet le logiciel CIVA développé et commercialisé par le Commissariat à l'Energie Atomique, France. Ce premier apprentissage a permis de repérer les paramètres influents, et de construire une première version du réseau de neurones sur la base d'imperfections virtuelles. La fonction de coût a été optimisée. La demanderesse a ensuite conduit un second apprentissage combinant les résultats obtenus sur simulation et des imperfections artificielles, c'est-à-dire créées intentionnellement sur des tubes réels. Ce second apprentissage permet de construire une seconde version du réseau de neurones, dont la fonction de coût a également été optimisée. La demanderesse a ensuite combiné les résultats obtenus sur des imperfections artificielles, et sur un ensemble d'imperfections présentes sur des tubes réels, ces imperfections étant connues avec précision par des mesures faites a posteriori hors la chaîne de fabrication. Cette troisième phase a permis de valider la dernière version du réseau de neurones. Cette version s'est avérée opérationnelle pour la surveillance en fabrication. Toutefois, lors de son implantation sur une installation nouvelle ou modifiée, il convient actuellement de lui faire subir un "calage", à l'aide d'une dizaine d'échantillons artificiels couvrant l'ensemble de la gamme des imperfections à traiter. Il s'ensuit naturellement une optimisation. Les figures 11, 12, 14 et 15 ont été décrites dans le cadre des capteurs P 11 et P12. Le même principe peut s'appliquer au groupe de capteurs P1. Dans ce cas, il n'y aura pas d'image 2, et le réseau construit a moins de paramètres d'entrée, comme déjà indiqué. Les circuits décrits pour deux capteurs peuvent être utilisés pour un seul, mais sans paramètres d'entrée pour la partie � Image 2 �. Le même principe peut s'appliquer aussi aux deux groupes de capteurs P21 et P22, chargés de détecter les imperfections transversales, en tenant compte du fait que les capteurs sont pour cette détection inclinés (par exemple de +/-17°) dans un plan passant par l'axe du tube. On comprendra que, dans chaque cas, il intervient un traitement numérique du type défini par la figure 11, éléments 992 à 996 exceptés. Ce traitement est globalement désigné par 763, conformément à la figure 8, où il est suivi des blocs 764 et 766. 15 On obtient ainsi un ensemble tel que représenté par la figure 18, avec : pour le capteur P1, un traitement 763-1, suivi d'une phase de décision et alarme 764-1 ; pour les capteurs PI l et P12, un traitement 763-10, suivi d'une phase de décision et alarme 764-10 ; 20 pour les capteurs P21 et P22, un traitement 763-20, suivi d'une phase de décision et alarme 764-20 ; les trois phases 764-1, 764-10 et 764-20 étant interprétées conjointement par l'automate de tri et d'alarme 767. Une variante de la figure 18, non représentée, consiste à ne prévoir qu'une seule phase 25 � Décision & alarme �, utilisant directement les sorties des trois traitements 763-1, 763-10 et 763-20. Le contrôle non destructif proprement dit se fait "à la volée", c'est-à-dire au fur et à mesure que le tube défile dans l'installation de contrôle. La décision issue du traitement 30 des informations décrit plus haut peut être prise elle aussi soit au fur et à mesure que le tube défile dans l'installation de contrôle (avec décision-alarme et marquage "à la volée") ; une variante consiste à prendre cette décision après que toute la longueur du 10 tube ait été inspectée, ou même encore plus tard (après le contrôle de l'ensemble d'un lot de tubes par exemple), chaque tube étant repéré/identifié (N° d'ordre par exemple). Dans ce cas, il est nécessaire que les informations obtenues soient enregistrées (mémorisées). Les enregistrements peuvent faire l'objet d'une analyse postérieure par un opérateur habilité à prendre une décision après avoir analysé les résultats enregistrés et traités par le(s) réseau(x) de neurones. Bien entendu, compte-tenu des propriétés des circuits neuronaux, il est possible de regrouper au moins partiellement l'ensemble des circuits neuronaux (contenus dans les traitements 763-1, 763-10 et 763-20) en un seul circuit neuronal, ayant toutes les entrées voulues. Le mode de réalisation décrit utilise directement des réseaux neuronaux à titre d'exemple de systèmes experts. L'invention n'est pas limitée à ce genre de réalisation. Ici, l'expression "agencement du genre circuit neuronal" peut couvrir d'autres techniques d'analyse de données non linéaires, avec ou sans circuits neuronaux. De façon générale, le convertisseur peut comprendre une entrée d'amplitude maximale dans un sélecteur et une entrée de temps de vol correspondant. Lesdites entrées peuvent fournir les données suffisantes pour la décision de conformité ou de non conformité d'un produit. Le bloc de transformation peut comprendre un élément d'élimination de données inutiles, un élément de filtrage de zones repérées, un simulateur et un élément d'interprétation. La réduction de la quantité d'information permet une vitesse de traitement plus élevée. Le simulateur peut comprendre un élément de simulation théorique, un calculateur de tolérance et un algorithme inverse. L'étage de sortie peut comprendre:30 - un combinateur, agencé pour préparer des entrées numériques de circuit neuronal, à partir d'un extrait des images correspondant à une zone d'imperfection présumée, de propriétés de l'imperfection présumée dans la même zone, issues du filtre, et de données de contexte, - au moins un circuit neuronal, qui reçoit des entrées issues du combinateur, - un étage numérique de décision et alarme, opérant sur la base de la sortie du circuit neuronal, et - un automate de tri et de marquage, agencé pour écarter et marquer des tubes décidés non conformes par l'étage numérique de décision et alarme. Le système ici proposé a été décrit dans le cas du contrôle non destructif lors de la fabrication de tubes sans soudure, cas auquel l'invention s'applique particulièrement bien. Les mêmes techniques peuvent s'appliquer notamment à des produits sidérurgiques longs non nécessairement tubulaires. Dans le cas de tubes soudés ou autres produits soudés (comme par exemple des tôles ou des plaques), le système s'avère capable de déterminer en plus les limites du cordon de soudure, et par conséquent de localiser les éventuelles imperfections dans le cordon de soudure, qui peuvent être à surveiller. De leur côté, les imperfections situées en dehors des limites du cordon de soudure, qui peuvent correspondre à des inclusions déjà présentes dans le feuillard (ou produit) de base, sont à considérer différemment. Annexe . Section 1 = F(E wz.iXi) N S1 = F (E Ei wi + wa ) =1 k S = ESiwi+wô i=i DP=SP ù Rp (14) P_M D2 Cg = ~P=1 P 2M (12) (13) Section 2 a p f a 1 e ' at ° dx = Q mb feuil 2~r stdb (seuil ù stdb seuil = stdb Q'1(pfa) + mb (22) (15) (21) Revendications 1. Dispositif formant outil d'aide à l'exploitation, pour le contrôle non destructif, en cours ou en sortie de fabrication, de produits sidérurgiques, tels que des tubes ou autres produits longs, cet outil étant destiné à tirer des informations sur d'éventuelles imperfections du produit, à partir de signaux de retour que captent (73), consécutivement à l'excitation sélective (70) de capteurs ultrasonores émetteurs selon une loi de temps choisie, des capteurs ultrasonores récepteurs formant un agencement de géométrie choisie, monté en couplage ultrasonore avec le produit par l'intermédiaire d'un milieu liquide, avec mouvement relatif de rotation/translation entre le tube et l'agencement de transducteurs, ledit outil d'exploitation étant caractérisé en ce qu'il comprend : - un convertisseur (891 ; 892) capable d'isoler sélectivement une représentation numérique de possibles échos dans des fenêtres temporelles désignées, en fonction du mouvement relatif de rotation/translation, ladite représentation comprenant l'amplitude et le temps de vol d'au moins un écho, et de générer un graphe 3D parallélépipédique, - un bloc de transformation (930) capable de générer une image 3D (901 ; 902) d'imperfections possibles dans le tube à partir du graphe 3D et d'une base de données, - un filtre (921 ; 922), capable de déterminer, dans les images (901 ; 902), des zones d'imperfection présumée (Zcur), ainsi que des propriétés de chaque imperfection présumée, et - un étage de sortie configuré pour générer un signal de conformité ou de non conformité d'un produit. 2. Dispositif selon la revendication 1, dans lequel le convertisseur (891 ; 892) 25 comprend une entrée d'amplitude maximale dans un sélecteur et une entrée de temps de vol correspondant. 3. Dispositif selon la revendication 1 ou 2, dans lequel le bloc de transformation comprend un élément d'élimination de données inutiles, un élément de filtrage de zones repérées, un simulateur et un élément d'interprétation. 30 4. Dispositif selon la revendication 3, dans lequel le simulateur comprend un élément de simulation théorique, un calculateur de tolérance et un algorithme inverse. 42 5. Dispositif selon l'une des revendications précédentes, dans lequel l'étage de sortie comprend: - un combinateur (960), agencé pour préparer des entrées numériques de travail, à partir d'un extrait (951 ; 952) des images correspondant à une zone d'imperfection présumée (Zcur), et de propriétés de l'imperfection présumée dans la même zone, issues du filtre (921 ; 922), - au moins un agencement du genre circuit neuronal (970), qui reçoit des entrées de travail issues du combinateur (960), - un étage numérique de décision et alarme (992) opérant sur la base de la sortie de 10 l'agencement du genre circuit neuronal (970), et - un automate de tri et de marquage (994) agencé pour écarter et marquer des produits décidés non conformes par l'étage numérique de décision et alarme (992). 6. Dispositif selon l'une des revendications précédentes, dans lequel ledit outil d'exploitation comprend deux convertisseurs (891,892) respectivement consacrés à 15 deux agencements de transducteurs (P 1 1,P12;P21,P22) ultrasonores de géométrie choisie (P11,P12;P21,P22), montés en couplage ultrasonore sensiblement selon une symétrie miroir de la direction de leurs faisceaux ultrasonores respectifs, et le combinateur (960) est agencé pour opérer sélectivement sur les échos de peau interne ou sur les échos de peau externe ou sur les échos intervenant dans la masse du tube, mais 20 en même temps sur les données relatives à l'un et l'autre des deux agencements transducteurs. 7. Dispositif selon l'une des revendications précédentes, dans lequel le convertisseur (891 ; 892) est agencé pour isoler sélectivement une représentation numérique de possibles maximums d'échos dans des fenêtres temporelles désignées 25 correspondant à des échos de peau interne, à des échos de peau externe, ainsi que des échos provenant de la masse du tube, respectivement, et le combinateur (960) est agencé pour opérer sélectivement sur les échos de peau interne ou sur les échos de peau externe ou sur les échos intervenus dans la masse. 8. Dispositif selon l'une des revendications précédentes, dans lequel le 30 combinateur (960) reçoit au moins une entrée (9511 ;9521) relative à un extremum d'amplitude de l'image dans la zone d'imperfection présumée. 9. Dispositif selon l'une des revendications précédentes, dans lequel le filtre (921 ; 922) est agencé pour produire, en tant que propriétés de chaque imperfection présumée, son obliquité et sa longueur, tandis que le combinateur (960) reçoit des entrées correspondantes d'obliquité d'imperfection (931) et de longueur d'imperfection (932). 10. Dispositif selon l'une des revendications précédentes, dans lequel le filtre (921 ; 922), le combinateur (960), le circuit neuronal (970) et l'étage numérique de décision et alarme (992) sont agencés pour opérer itérativement sur une suite de zones d'imperfection présumée (Zcur), déterminées par ledit filtre (921 ; 922). 11. Dispositif selon la revendication 10, dans lequel le filtre (921 ; 922), le combinateur (960), le circuit neuronal (970) et l'étage numérique de décision et alarme (992) sont agencés pour opérer alternativement sur la peau interne et la peau externe du tube. 12. Dispositif selon l'une des revendications précédentes, dans lequel ledit agencement du genre circuit neuronal comprend : - un premier circuit neuronal (NC 121- NC 123) propre à évaluer la nature d'une imperfection parmi une pluralité de classes prédéfinies, et - un second circuit neuronal (NC141- NC143) propre à évaluer la gravité d'une imperfection. 13. Dispositif selon la revendication 12, dans lequel les deux circuits neuronaux ont 20 des entrées qui diffèrent par : - une entrée (9512 ; 9522) de nombre de maxima voisins pour le premier circuit neuronal, et -une entrée (9513 ; 9523) de largeur d'écho pour le second circuit neuronal. 14. Dispositif selon l'une des revendications 12 et 13, dans lequel les sorties des 25 deux circuits neuronaux sont combinées pour affiner la prédiction. 15. Dispositif selon l'une des revendications précédentes, dans lequel l'émission et la réception des signaux ultrasonores sont réalisées à chaque fois par un même transducteur, pour une partie au moins de l'agencement de capteurs. 16. Dispositif de contrôle non destructif de tubes en cours ou en sortie de 30 fabrication, comprenant:- un agencement de transducteurs ultrasonores de géométrie choisie, monté en couplage ultrasonore avec le tube par l'intermédiaire d'un milieu couplant, avec mouvement relatif de rotation/translation entre le tube et l'agencement de transducteurs, - des circuits pour exciter sélectivement (70) ces éléments transducteurs selon une loi de temps choisie, et pour recueillir (73) les signaux de retour qu'ils captent, et - un outil d'aide à l'exploitation selon l'une des revendications précédentes. 17. Procédé de contrôle non destructif de produits sidérurgiques, tels que des tubes ou autres produits longs, en cours ou en sortie de fabrication, comprenant les étapes suivantes : a. prévoir un agencement de transducteurs ultrasonores de géométrie choisie, monté en couplage ultrasonore avec le tube par l'intermédiaire d'un milieu couplant, avec mouvement relatif de rotation/translation entre le tube et l'agencement de transducteurs, b. exciter sélectivement (70) ces éléments transducteurs selon une loi de temps choisie, c. recueillir (73) les signaux de retour qu'ils captent, afin d'analyser sélectivement ces signaux de retour (760-766), pour en tirer des informations sur d'éventuelles imperfections du tube, lesdites informations comprenant l'amplitude et le temps de vol d'au moins un écho, et générer un graphe 3D parallélépipédique, d. isoler sélectivement une représentation numérique de possibles échos dans des fenêtres temporelles désignées, en fonction du mouvement relatif de rotation/translation (891 ; 892), et en tirer une image 3D (901 ; 902) d'imperfections possibles dans le tube à partir du graphe 3D et d'une base de données, e. générer un signal de conformité ou de non conformité d'un produit. 18. Procédé selon la revendication 17, dans lequel l'étape e., comprend: el. filtrer (921 ; 922) les images (901 ; 902) selon des critères de filtrage choisis, afin d'y déterminer des zones d'imperfection présumée (Zcur), ainsi que des e2. propriétés de chaque imperfection présumée, 952) former (960) des entrées numériques de travail, à partir d'un extrait (951 ; des images correspondant à une zone d'imperfection présumée (Zcur), de ,propriétés de l'imperfection présumée dans la même zone, issues du filtre (921 ; 922), et de données de contexte (740), e3. appliquer les entrées ainsi formées (960) à au moins un agencement du genre circuit neuronal (970), e4. traiter numériquement la sortie de l'agencement du genre circuit neuronal (970) selon des critères de décision choisis, pour en tirer une décision et/ou une alarme (992), et e5. écarter et marquer (994) des tubes décidés non conformes par l'étape e4.