FR 2596270 A1 19871002 FR 8604225 A 19860325 APPAREIL D'EXPLORATION PAR ECHOGRAPHIE ULTRASONORE DE MILIEUX EN MOUVEMENT ET NOTAMMENT D'ECOULEMENTS SANGUINS OU D'ORGANES TELS QUE LE COEUR La présente invention concerne un appareil d'exploration par échographie ultrasonore de milieux en mouvement et notamment d'écoulements sanguins ou d'organes tels que le coeur, en vue de la détermination des paramètres de mouvement de tels milieux, ledit appareil comprenant au moins un transducteur ultrasonore associé à un étage d'émission périodique d'un signal impulsionnel de fréquence de récurrence déterminée F = 1/T et à un étage de réception des signaux échographiques renvoyés vers le transducteur et de traitement des signaux reçus. Depuis plusieurs années, des systèmes Doppler ultrasonores h ondes pulsées sont couramment utilisés pour mesurer en un point donné les vitesses d'écoulements sanguins, ou du moins la projection de telles vitesses sur l'axe du faisceau émis par les transducteurs ultrasonores. Plus récemment sont apparus des appareils permettant d'obtenir en temps réel la distribution des vitesses d'écoulements sur le trajet parcouru par l'onde ultrasonore et, plus encore, sur le plan de coupe obtenu en balayant le transducteur. La plupart de ces systèmes exploitent le décalage de fréquence ou le décalage de phase du signal rétrodiffusé par les cibles mobiles pour en déduire la vitesse axiale des écoulements sanguins. La demande de brevet européen n 0092841, par exemple, a trait à un tel appareil. Cependant, alors que ladite demande utilise la mesure du décalage de phase entre les échos successifs rétrodiffusés, en réponse à une excitation récurrente, par les cibles en mouvement, le but de l'invention est de proposer un autre type de traitement de signal dans le domaine temporel, avec lequel on interprète la rétrodiffusion des signaux ultrasonores non plus en ces termes de décalage de fréquence ou de phase mais en termes de décalage temporel progressif des signaux échographiques à la suite de chaque envoi de signaux impulsionnels. A cet effet, l'appareil selon l'invention est caractérisé en ce qu'il comprend une voie de traitement numérique composée successivement d'un circuit de suppression des échos fixes, d'un circuit d'estimation de paramètres d'écoulement, d'un circuit discriminateur, d'un dispositif de mémorisation, de conversion de balayage et de codage de couleur, et d'un dispositif de visualisation. Dans cette structure le principe d'estimation en temps réel des paramètres d'écoulement repose sur la formulation temporelle des signaux ultrasonores rétrodiffusés par les cibles en mouvement. Les signaux échographiques successivement obtenus en réponse aux excitations récurrentes du transducteur fixe peuvent se mettre sous la forme : ei+î(t) = ei(t - T(t)) (2) où i est le rang de la ligne échographique et où T(t) traduit le décalage temporel induit par le déplacement des cibles entre deux tirs. Le décalage est lié à la vitesse axiale locale VZ(t) (projection sur l'axe Z de propagation de l'onde ultrasonore) par la relation T(t) = 2 . Vz(t) . T (3) c où c est la vitesse de propagation des ultrasons, et T la période de récurrence des excitations.La localisation de la vitesse d'écoulement est traduite par la relation classique en échographie liant le temps t à la profondeur Z 1 Z = 2 . ct (4) La relation de récurrence (2) est applicable également aux signaux obtenus en ne prenant en compte que le signe des signaux échographiques correspondant aux lignes échographiques, ce qui donne si+l(t) = signe de (ei+l (t)) signe de (ei+î (t- T(t))) = Si (t- T(t)) (5) L'information contenue dans ces signaux réduits au signe si(t) est en effet suffisante pour estimer localement des vitesses moyennes, car les fronts montants ou descendants d'un tel signal si+l(t) se révèlent être des répliques, retardées de T (t), du signal précédent si(t) On dispose donc d'une estimation T de T(t) pour chaque front montant ou descendant du signal Si en mesurant le décalage temporel entre ce front, montant ou descendant, et le front correspondant du signal suivant si+l. Une information est ainsi disponible toutes les demi-périodes du signal échographique. Une multiplication 1 bit s'avère d'ailleurs être particulièrement adaptée à la mesure de T lorsque les signaux sont différents, le résultat de la multiplication 1 bit est 0, alors qu'il est égal à 1 si ces signaux sont identiques. De cette manière, la fonction fi(to) définie ci-dessous fi(to) = si+î(t).dt (6) où to et to + w encadrent une fenêtre temporelle de longueur w, vérifie la relation où V est la vitesse estimée et vO la fréquence centrale. En fait, ce calcul ne permet d'obtenir que l'amplitude de la vitesse et non son sens qui représente pourtant une information très importante pour le clinicien. Pour lever cette ambiguité quant au sens de la vitesse, il faut considérer que fi(to) est la valeur prise en O par la fonction de corrélation locale, définie par C(t0,u) = #i+1 (t+u).dt (8) Cette fonction a les propriétés suivantes. D'une part, si la vitesse est positive, T(t) est supérieur à O et l'expression C(t0,u) est une fonction croissante au voisinage de u = O, tandis que cette même expression C(to,u) est une fonction décroissante au voisinage de u = O si la vitesse est au contraire négative. D'autre part, au voisinage du pic de corrélation, C(to, u) est triangulaire. Ainsi, le calcul de la pente A de C(to,u) au voisinage de u = O permet de déterminer le sens de la vitesse ainsi que la valeur du pic de corrélation Cmax. Cette valeur est en effet utilisable pour estimer la variance a de la vitesse selon la relation # # (1 - Cmax /2 (9) Cmax est en fait calculé de la manière suivante ou encore Une unité de traitement comprenant un corrélateur 1 bit à trois voies de calcul et une unité de calcul fournira donc une estimation des paramètres V et aV. Dans l'appareil selon l'invention, qui peut comprendre une voie de traitement analogique traditionnelle pour la visualisation de plans de coupe explorés selon le principe de l'échographie classique, le circuit de suppression des échos fixes comprend, de préférence, successivement un convertisseur analogique-numérique, commandé par des signaux de com mande du rythme d'échantillonnage à une fréquence f = 1/ a t, et un soustracteur recevant sur sa première entrée de signe déterminé directement la sortie du convertisseur et sur sa deuxième entrée de signe contraire cette même sortie mais re tardée d'un nombre entier de périodes a t d'échantillonnage dans un circuit à retard, et délivrant un signal dit de différence.Le signal de différence ainsi délivré est utilisé à la fois dans le circuit suivant d'estimation de paramètres d'écoulement et dans le circuit discriminateur. Dans le mode de réalisation proposé, le circuit d'estimation de paramètres d'écoulement comprend d'une part un circuit d'intercorrélation qui, à partir desdits signaux de différence relatifs à deux lignes échographiques successives, délivre un nombre déterminé de valeurs de fonction de corrélation, et d'autre part un circuit d'interpolation qui, à partir desdites valeurs, délivre des paramètres représentatifs de la vitesse des milieux explorés le long de l'axe de propagation des ondes ultrasonores. Le circuit d'intercorrélation comprend trois voies en parallèle composées chacune d'un corrélateur 1 bit commandé par lesdits signaux de commande à la fréquence d'échantillonnage F, lesdits corrélateurs recevant sur une première entrée directement ledit signal de différence délivré par le circuit de suppression des échos fixes et sur une deuxième entrée ce même signal mais retardé dans des lignes à retard dont les retards prennent trois valeurs distinctes T - At, T, T + qat, et délivrant trois valeurs de fonction de corrélation, et des circuits de calcul de moyenne étant avantageusement prévus dans chacune des trois voies en parallèle du circuit d'intercorrélation, en sortie des corrélateurs. Dans la réalisation proposée, le circuit discriminateur comprend successivement un circuit d'élévation au carré du signal de différence délivré par le circuit de suppression des échos fixes, un sommateur de calcul de l'énergie locale du signal de différence, et un circuit de validation, ou non, des signaux de sortie du circuit d'estimation de paramètres d'écoulement selon la valeur calculée de ladite énergie locale. Il est en outre avantageux de prévoir dans ce circuit discriminateur un autre circuit de validation, ou non, des paramètres d'écoulement selon la valeur de l'un au moins de ces paramètres. Les particularités et avantages de l'invention apparaîtront maintenant de façon plus précise dans la description qui suit et dans les dessins annexés, donnés à titre d'exemples et dans lesquels - la figure 1 montre un exemple de réalisation préférentiel de l'appareil selon l'invention - les figures 2, 3, 4 représentent des exemples particuliers de réalisation respectivement de l'étage d'émission, du circuit de suppression des échos fixes, et du circuit d'estimation de paramètres d'écoulement - la figure 5 met en évidence le principe de fonctionnement de l'unité de calcul. L'appareil représenté sur la figure 1 comprend, de façon classique, un transducteur ultrasonore 10 associé d'une part à un étage d'émission 20 et d'autre part à un étage 30 de réception et de traitement, ainsi qu'à un dispositif 40 de commande de balayage mécanique du transducteur. A la place de ce transducteur pourrait bien entendu être utilisé un réseau de transducteurs, associé alors à un dispositif de commande de balayage électronique. Dans la réalisation décrite de façon plus détaillée sur la figure 2, l'étage d'émission 20 comprend un générateur 21 de signaux électriques d'excitation, envoyés vers le transducteur 10 qui les convertit en trains périodiques de signaux impulsionnels ultrasonores. Cette émission est commandée par des signaux d'horloge disponibles sur une connexion 102 et délivrés à une fréquence de récurrence F déterminée -de l'ordre de 5 kilohertz par exemple- par un séquenceur comprenant successivement un oscillateur 22, ici de fréquence 32 mégahertz, et un diviseur de fréquence 23. Ce divi-seur délivre les signaux d'horloge sur la connexion 102, ainsi que d'autres signaux de commande respectivement sur des connexions 104 et 106, à 1 kilohertz et 16 mégahertz dans l'exemple ici décrit.Les signaux de commande présents sur la connexion 104 commandent en particulier le dispositif 40, pour le balayage du transducteur. Un séparateur 24 des étages d'émission 20 et de réception et de traitement 30 est inséré entre le transducteur 10, le générateur 21, et ledit étage 30, et évite l'aveuglement des circuits de réception par les signaux d'émission. L'étage de réception et de traitement 30 comprend, en sortie du séparateur 24, un amplificateur haute fréquence 300 incluant la compensation de gain en fonction de la profondeur et suivi de deux voies de traitement 301 et 302 prévues en parallèle. La voie 301, de type traditionnel, comprend ici, en série, un détecteur d'enveloppe 310, un amplificateur de compression logarithmique 311, un dispositif de mémorisation et de conversion de balayage 370 incluant aussi une fonction de codage de couleur, et un dispositif 312 de visualisation. Cette voie 301 permet l'obtention, sous forme d'images en échelle de gris, de plans de coupe des milieux explorés selon le principe de l'échographie classique. Conformément à l'invention, la voie 302 comprend en série un circuit 320 de suppression des échos fixes, un circuit 630 d'estimation de paramètres d'écoulement, un circuit discriminateur 360, le dispositif de mémorisation, de conversion de balayage et de codage de couleur 370, et le dispositif de visualisation 312. Le circuit numérique 320 de suppression des échos fixes comprend lui-même, dans l'exemple de réalisation représenté sur la figure 3, un convertisseur analogique-numérique 321 dont la sortie est reliée d'une part directement à l'entrée négative d'un soustracteur 322 et d'autre part à l'entrée positive de ce soustracteur par l'intermédiaire d'un circuit à retard 323. Le retard apporté par ce dernier circuit 323 pourrait être de plusieurs périodes T = 1/F, mais il est avantageux que ce retard soit le plus faible possible et égal à T. Ce circuit 320 est prévu pour l'élimination de tous les échos fixes, en particulier de ceux dont l'apparition est provoquée par la réflexion des ondes ultrasonores contre les parois des vaisseaux où se produisent les écoulements étudiés. Cette présence d'échos fixes est gênante du fait de leur amplitude bien plus élevée (de l'ordre de +40 dB dans le cas des écoulements sanguins) que celle des signaux utiles, les signaux rétrodiffusés par les cibles mobiles. Le circuit 320 est également commandé, par l'intermédiaire de la connexion 106, par le diviseur de fréquence 23 du séquenceur qui lui fournit le signal de commande d'échantillonnage à la fréquence de 16 mégahertz. Le circuit 630 d'estimation de paramètres d'écoulement, représenté sur la figure 4, comprend ici un circuit d'intercorrélation et une unité de calcul 650. A partir des signaux de différence entre deux lignes échographiques successives d'échantillons di(t), di+l(t), etc... (où i représente le rang de ce signal) successivement fournis par le circuit 320 de suppression des échos fixes, le circuit d'intercorrélation délivre trois valeur de fonction de corrélation. A partir de ces trois valeurs, l'unité de calcul délivre des paramètres caractérisant les différents écoulements rencontrés le long de l'axe de propagation de l'onde ultrasonore. Ces paramètres sont ici la composante axiale de la vitesse locale moyenne V7 et la variance locale a 2 de celleci, le mot local étant ici employé au sens de la localisation en profondeur le long dudit axe Z de propagation. Plus précisément, le circuit d'intercorrélation 630 comprend trois corréla#teurs 642 qui, sur une première entrée, reçoivent chacun directement la sortie di+i(t) du circuit 320 de suppression d'échos fixes et, sur une deuxième entrée, reçoivent cette même sortie du circuit 320, mais retardée par des lignes à retard 641 et correspondant donc au signal précédent di(t). De plus, chacune de ces lignes 641 présente, pour permettre le calcul des trois valeurs de la fonction de corrélation, un retard distinct prenant trois valeurs T - At, T, T + At, où at représente la période d'échantillonnage imposée par la connexion 106. Ce calcul, effectué en parallèle, des trois valeurs de la fonction de corrélation utilise K échantillons successifs des deux signaux d'entrée des corrélateurs. Les groupes de K échantillons définissent des fenêtres temporelles successives de longueur K At se décalant progressivement au rythme de la fréquence imposée par la connexion 106.La fonction de corrélation est définie par une expression du type k=K fi(J,P) = E di((k+J) #t). di+i( k=I dans laquelle - J détermine le début de la fenêtre temporelle de longueur K at - P représente le décalage temporel introduit entre di et di+1 et pour lequel on calcule la valeur de la fonction de corrélation (P variant de -1 à +1) par pas de 1) - i représente le rang de la différence entre deux lignes échographiques successives ei et ex+1. Les corrélateurs 642, commandés par l'intermédiaire de la connexion 106 de sortie du diviseur de fréquence 23 du séquenceur, sont, on l'a vu, des corrélateurs 1 bit.L'unité de calcul 650 peut être une unité de calcul programmée, incluant un microprocesseur, ou de préférence une unité de calcul câblée.Cette unité de calcul dispose de trois valeurs de la fonction de corrélation C1 = C(t0, - At) (12) C2 = C(to, O) (13) C3 = C(to, Lt) (14) Le signe de la pente de la fonction de corrélation, auquel correspond le sens de la vitesse, est déterminé par :: s = signe (C3 - C1) (15) En utilisant la relation (7), la vitesse est donnée par l'expression La valeur A de la pente en O de la fonction de corrélation est obtenue selon le tableau ci-après C2 b 0,5 C2 0,5 C3-Cl u O A = (C2-C1)/ At A = (C3-C2)/ #t C3-C1 O A = (C3-C2)/ At A = (C2-C1)/ Qt Ce tableau fait l'inventaire, en association avec les figures Sa à 5d correspondantes, des quatre types de configuration possible décrivant les valeurs Cl, C2, C3. On en déduit la valeur CMAX du pic de corrélation ainsi que la variance c Entre les trois corrélateurs et les trois entrées correspondantes de l'unité de calcul 650 peuvent être prévus, et le sont ici de préférence, trois circuits de calcul de moyenne qui sont en fait des accumulateurs comprenant chacun un additionneur 344 et une ligne à retard 345 de retard T (ou un multiple de T). Ces circuits de calcul de moyenne permettent d'accumuler les valeurs de fonction de corrélation sur N lignes échographiques successives et d'en effectuer la moyenne avant envoi vers l'unité de calcul. Les additionneurs 344 et les lignes à retard 345 sont reliés au séquenceur (22, 23) par l'intermédiaire de la connexion 104, pour leur remise à zéro à intervalles N x T réguliers.La fréquence de ces remises à zéro imposées par la connexion 104 est celle des signaux présents sur la connexion 102, divisée par N. Les signaux de sortie du circuit 630 d'estimation des paramètres d'écoulement sont alors validés, ou non, par le circuit discriminateur 360, puis les valeurs ainsi confirmées sont envoyées vers le dispositif de visualisation 312 par l'intermédiaire du dispositif de codage de couleur 370. La présence du circuit discriminateur 360 est indispensable. En effet, en dehors des zones d'écoulement, le signal de sortie du circuit 320 de suppression des échos fixes est essentiellement du bruit. Le résultat alors fourni par le circuit 630 d'estimation des paramètres de l'écoulement, qui traite ce bruit, n'est pas l'indication d'une vitesse nulle, et il devient nécessaire de valider, ou non, ce résultat.Le circuit 360 comprend à cet effet en série un multiplieur qui reçoit le signal de sortie di du circuit 320 de suppression des échos fixes et élève ce signal de différence au carré, un sommateur qui permet le calcul de l'énergie locale de ce signal de différence selon la formule k=K Ei(J) = E d2 ((k+J) #t) (11) k=1 un circuit de calcul de moyenne qui, comme dans le cas des circuits 344 et 345, est un accumulateur comprenant un additionneur et une ligne à retard de retard T (ou un multiple de T) et permet d'effectuer la moyenne de l'énergie locale sur N tirs c'est-à-dire (N-l) différences selon l'expression i = N-1 E(J) = E Ei(J) (12) i=1 avant d'envoyer la valeur ainsi obtenue vers un circuit de validation. Ce circuit de validation comprend un comparateur recevant d'une part sur une première entrée la sortie de l'accumulateur, ou directement celle du sommateur dans le cas où le circuit de calcul de moyenne ne serait pas prévu, et d'autre part sur une deuxième entrée une tension de référence formant seuil, la sortie du comparateur étant ici au niveau logique 0 ou 1 selon que la tension reçue sur sa première entrée est respectivement inférieure ou supérieure au seuil de référence. Deux multiplieurs qui reçoivent respectivement les signaux de sortie Vz et o2 du circuit 630 sur leur première entrée#transmettent ces deux signaux, notés désormais V'z et a'2, sur leur sortie respective ou transmettent simplement des valeurs nulles selon que le signal de validation délivré sur leur deuxième entrée par le comparateur est 1 ou 0 respectivement.En effet, en dehors des zones d'écoulement vrais, l'énergie moyenne calculée en sortie de l'accumulateur est celle du bruit seul, et peut être mesurée seule, a priori, en l'absence d'excitation, en vue de déterminer la valeur appropriée du seuil. Au. contraire, en présence de signaux rétrodiffusés par les cibles en mouvement, l'énergie moyenne du signal di est supérieure à celle du bruit seul, ce qui autorise la validation des signaux délivrés par le circuit 630 d'estimation des paramètres de l'écoulement. Dans une variante de réalisation de ce circuit discriminateur 360, le signal de variance lui-même peut être utilisé comme signal de discrimination. En effet, comme d'un tir à l'autre le bruit est décorrélé, si la variance calculée est grande, c'est parce' qu'on se trouve en présence de bruit seul. Au contraire, en présence de signaux rétrodiffusés, la variance est nettement plus faible et, la encore, la présence d'un autre comparateur, recevant sur une première entrée le signal de sortie c@2 et sur une deuxième entrée une tension de référence constituant un seuil, et de deux autres multiplieurs permet de valider ou non, comme précédemment, les signaux appelés alors VYZ et c"2 en sortie de ces multiplieurs.Bien entendu, au lieu d'être situé en aval des premiers multiplieurs du circuit de validation, ce nouveau circuit de validation pourrait être placé en amont du premier circuit de validation, les deux signaux de vitesse et de variance éventuellement validés par lui-même étant alors envoyés sur les entrées correspondantes des nouveaux multiplieurs 463 et 464. A la suite des N excitations commandées par le signal 102 et ayant permis la mesure des paramètres d'écoulement ou de déplacement des structures situées sur l'axe du transducteur, un déplacement électronique ou mécanique dudit transducteur est commandé par le signal de commande de balayage 104. Le procédé d'estimation est alors réinitialisé dans une autre direction d'observation. Les deux sorties du circuit discriminateur 360 sont envoyées vers le dispositif de mémorisation, conversion de balayage et codage de couleur 370, par lequel transite également, avant visualisation, le signal de sortie de l'amplificateur 311 de la voie de traitement 301. Un tel dispositif est décrit par exemple dans la demande de brevet européen EP-A--0100094. La figure 3 de ce document montre en effet, entre les bornes A, B, C et ER, EG, Eg, un exemple des circuits qui peuvent être mis en place dans le cas de la présente invention, la borne A recevant le signal échographique classique et les bornes B et C les paramètres caractéristiques du milieu en mouvement examiné. Ce dispositif 370 et le dispositif 312 permettent alors la visualisation en temps réel des écoulements ou des déplacements superposés à l'image de réflectivité échographique classique. REVENDICATIONS 1. Appareil d'exploration de milieux en mouvement par échographie ultrasonore, et notamment d'écoulements sanguins, en vue de la détermination des paramètres de mouvement de tels milieux, ledit appareil comprenant au moins un transducteur ultrasonore associé à un étage d'émission périodique d'un signal impulsionnel de fréquence de récurrence déterminée F = 1/T et à un étage de réception des signaux échographiques renvoyés vers le transducteur et de traitement de ces signaux re çus, et étant caractérisé en ce qu'il comprend une voie de traitement numérique composée successivement d'un circuit de suppression des échos fixes, d'un circuit d'estimation de paramètres d'écoulement, d'un circuit discriminateur, d'un dispositif de mémorisation, de conversion de balayage et de codage de couleur, et d'un dispositif de visualisation. 2. Appareil selon la revendication 1, caractérisé en ce qu'il lui est associé une voie de traitement analogique traditionnelle, pour la visualisation de plans de coupe explorées selon le principe de l'échographie classique. 3. Appareil selon l'une des revendications 1 et 2, caractérisé en ce que le circuit de suppression des échos fixes comprend successivement un convertisseur analogique-numérique,commandé par des signaux de commande du rythme d'échantillonnage à une fréquence f = 1/ Qt, et un soustracteur recevant sur sa première entrée de signe déterminé directement la sortie du convertisseur et sur sa deuxième entrée de signe contraire cette même sortie mais retardée d'un nombre entier de périodes Qt d'échantillonnage dans un circuit à retard, et délivrant un signal dit de différence. 4. Appareil selon la revendication 3, caractérisé en ce que le circuit d'estimation de paramètres d'écoulement comprend d'une part un circuit d'intercorrélation qui, à partir desdits signaux de différence relatifs à deux lignes échographiques successives, délivre trois valeurs de fonction de corrélation, et d'autre part une unité de calcul qui, à partir desdites valeurs, délivre des paramètres représentatifs de la vitesse des milieux explorés le long de l'axe de propagation des ondes ultrasonores. 5. Appareil selon la revendication 4, caractérisé en ce que le circuit d'intercorrélation comprend trois voies en parallèle composées chacune d'un corrélateur 1 bit commandé par lesdits signaux de commande à la fréquence d'échantillonnage F, lesdits corrélateurs recevant sur une première entrée directement ledit signal de différence délivré par le circuit de suppression des échos fixes et sur une deuxième entrée ce même signal mais retardé dans des lignes à retard dont les retards prennent trois valeurs distinctes T - Qt, T, T + Qt, et délivrant trois valeurs de fonction de corrélation. 6. Appareil selon la revendication 5, caractérisé en ce qu'un circuit de calcul de moyenne est inséré dans chacune des trois voies en parallèle du circuit d'intercorrélation, entre la sortie des corrélateurs et les entrées respectives de l'unité de calcul. 7. Appareil selon la revendication 6, caractérisé en ce que les trois circuits de calcul de moyenne comprennent chacun un accumulateur constitué successivement d'un additionneur, recevant sur une première entrée la sortie du corrélateur correspondant, et d'une ligne à retard de retard égal à T ou un multiple de T dont la sortie est renvoyée sur la deuxième entrée de l'additionneur, les sorties des additionneurs constituant les sorties des trois voies en parallèle. 8. Appareil selon l'une des revendications 6 et 7, caractérisé en ce que le circuit discriminateur comprend successivement un circuit d'élévation au carré du signal de différence délivré par le circuit de suppression des échos fixes, un sommateur de calcul de l'énergie locale du signal de différence, et un circuit de validation, ou non, des signaux de sortie du circuit d'estimation de paramètres d'écoulement selon la valeur calculée de ladite énergie locale. 9. Appareil selon la revendication 8, caractérisé en ce qu'il comprend un deuxième circuit de validation, ou non, desdits signaux de sortie du circuit d'estimation de paramètres d'écoulement selon la valeur de l'un au moins desdits paramètres d'écoulement. 10. Appareil selon l'une des revendications 8 et 9, caractérisé en ce que chaque circuit de validation comprend un comparateur à seuil et deux multiplieurs placés en série sur chaque voie de transfert des paramètres estimés pour transmettre ces deux paramètres ou simplement des valeurs nulles selon que le signal de validation délivré par le comparateur est 1 ou O respectivement. 11. Appareil selon l'une des revendications 8 à 10, caractérisé en ce que le circuit discriminateur comprend un accumulateur inséré en série entre la sortie du sommateur et l'entrée du premier circuit de validation pour calculer l'énergie locale moyenne du signal de différence.