Brevet FR2604081A1 — publié le 19880325

- 1 - La présente demande concerne un procédé et un dis- positif d'exploration par échographie ultrasonore d'un milieu en mouvement, permettant de mesurer la vitesse axiale dudit mouvement, en projection sur l'axe 'z" d'un faisceau d'excita- vions ultrasonores périodiques émis, à une période de récur- rence "T", par au moins un transducteur ultrasonore. L'objet de la présente demande trouve une applica- tion particulièrement avantageuse dans le domaine de l'explo- ration échographique d'organes en mouvement, tels que les pa- rois du coeur, et d'écoulements sanguins. Le problème technique à résoudre par tout procédé et tout dispositif d'exploration par échographie ultrasonore d'un milieu en mouvement, permettant de mesurer la Vwtes = axiale dudit mouvement, en projection sur l'axe "z" d'un.ais- ceau d'excitations ultrasonores périodiques émis, à une pério- de de récurrence "T", par au moins un transducteur ultrasonc- re, est d'obtenir une estimation aussi exacte que possible de la vitesse axiale du mouvement étudié afin de réaliser, à l'aide de dispositifs de visualisation, des images précises des organes et des écoulements sanguins soumis à une explora- tion échographique ultrasonore. Depuis plusieurs années, diverses solutions à ce problème technique ont été proposées, en particulier les sys- tèmes DOPPLER ultrasonores à onde pulsée, couramment utilisés pour mesurer en un point donné les vitesses d'écoulements san- guins, 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. - 2- 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 san- guins. La demande de brevet européen n O 092 841, par exem- ple, concerne un tel appareil utilisant la mesure du décalage de phase entre les échos successifs rétrodiffusés, en réponse à une excitation récurrente, par des cibles en mouvement. Cependant, les appareils de mise en oeuvre de cet- te méthode connue utilisant le décalage de phase sont limités par une relation d'incertitude liant la résolution axiale Az et la précision de la mesure de vitesse AV/V à la longueur d'onde A: AV A Az =- V 2 Cette relation, citée dans le chapitre II, para- graphe 2.3-a, de l'ouvrage "Doppler Ultrasound and Its Use in Clinical Measurement", P. ATKINSON and J.P. WOODCOCK, Academ.- Press, 1982, impose donc un compromis entre la résolution axiale et la précision de la mesure de vitesse, ce qui est in- compatible avec une mesure précise d'un profil de vitesse ou d'une image d'écoulements sanguins. En ce. sens, la demande de brevet français n 85 17 851, déposée le 03.12. 1985, décrit un autre type de traitement du signal échographique dans lequel la précision de la mesure de la vitesse n'est pas limitée par la résolution spatiale. Il s'agit d'un procédé d'analyse temporelle qui con- siste à interpréter la rétrodiffusion des signaux ultrasono- res, non plus en termes de décalage de fréquence ou de phase, mais en termes de décalage temporel des signaux échographiques à la suite de chaque envoi de signaux impulsionnels. Ce procé- dé repose sur le principe suivant: soit un milieu en mouve- ment à la vitesse axiale Vz(t), en projection sur l'axe "z d'un faisceau d'excitations ultrasonores périodiques émis, à une période de récurrence "T", par un transducteur ultrasono- re. Si à l'instant t=O, alors que ledit milieu est à une dis- 3 - tance "z" dudit transducteur, celui-ci émet une première exci- tation ultrasonore, l'écho el(t) rétrodiffusé par le milieu sera reçu à l'instant t1 = 2z/C, C étant la vitesse de propa- gation de l'onde ultrasonore. Puis, si à l'instant t=T, le OS transducteur émet une deuxième excitation ultrasonore, l'écho e2(t) renvoyé par le milieu en mouvement sera détecté par le transducteur à l'instant t = T + 2[z + Vz(t)T]/C, le milieu ayant parcouru la distance Vz(t)T pendant la période "T". Prenant pour origine du temps pour le deuxième écho l'instant t=T (origine de l'excitation correspondante), on obtient la relation e2(t) = el[t-T(t)] avec T(t) = 2Vz(t)T/C. Cette re- lation est très générale et peut s'écrire: ei+l(t) = ei[t-Ti(t)] (1) Ti(t) = 2Viz(t)T/C (2) Ti(t) étant alors le décalage temporel induit par le dépla- cement du milieu entre les tirs "i" et "i+1". On voit, par conséquent, sur la formule (2) que l'on peut mesurer la vitesse axiale Viz(t) à partir du déca- lage temporel Ti(t), celui-ci étant extrait de la formule (1) par un traitement approprié. La demande de brevet français n 85 17 851 propose comme méthode d'extraction du décalage temporel Ti(t) l'usa- ge des fonctions d'intercorrélation, le décalage temporel cherché étant celui pour lequel la fonction d'intercorrélation entre deux échos successifs ei(t) et ei+l(t) est maxima- le. Néanmoins, cette méthode, si elle a l'avantage de donner la vitesse axiale Viz(t) de façon précise, présente l'incon- vénient d'un dispositif de mise en oeuvre complexe exigeant, en plus d'un étage d'émission et d'un étage de réception et de traitement des signaux échographiques renvoyés vers le trans- ducteur, de nombreux circuits corrélateurs, autant de circuits moyenneurs, et un circuit d'interpolation à microprocesseur ou à logique câblée. 4 - Le problème technique général à résoudre par l'ob- jet de la présente demande est de proposer un procédé et un dispositif d'exploration par échographie ultrasonore d'un mi- lieu en mouvement, permettant de mesurer la vitesse axiale du- dit mouvement, en projection sur l'axe 'z" d'un faisceau d'ex- citations ultrasonores périodiques émis, à.une période de ré- currence "T", par au moins un transducteur ultrasonore, ledit dispositif comprenant, en outre, un étage d'émission, un étage de réception et de traitement des signaux échographiques ren- voyés vers le transducteur, procédé et dispositif grâce aux- quels on obtient, à la fois, une mesure précise de la vitesse axiale sans limitation due à la résolution axiale, et une mise en oeuvre aisée à partir de circuits électroniques simples. Dans un mode particulier de réalisation relatif aux écoule- ments sanguins, il s'avère que le signal rétrodiffusé par le sang est très faible par rapport aux échos fixes. Le disposi- tif proposé devra alors permettre d'extraire le signal recher- ché malcré la présence de ces échos fixes parasites. La solution au problème technique général à résou- dre consiste, selon le procédé objet de la présente demande, en ce qu'on réalise une analyse temporelle du signal échogra- phique reçu par ledit transducteur ultrasonore de façon à en extraire le décalage temporel Ti(t) induit par le mouvement du milieu entre deux échos successifs ei(t) et ei+1(t), à l'aide d'un développement limité de la formule: ei+1 = eitt-Ti(t)], reliant les deux dits échos, la vitesse axiale recherchée Viz(t) étant égale à Ti(t)C/2T, o 'C" représente la vi- tesse de propagation de l'onde ultrasonore. Corrélativement, le dispositif selon l'invention est remarquable en ce que le- dit étage de réception et de traitement des signaux échogra- phiques comprend notamment un circuit d'estimation de ladite vitesse axiale Viz(t) comportant, d'une part, un circuit - d'extraction du décalage temporel Ti(t) entre deux échos consécutifs ei(t) et ei+1(t), constitué pour résoudre l'équation en Ti(t): P (-l)n (n) n ei+i(t) =O n!E -ei (t).Ti (t) (3) n=O n! o ei(n)(t) est la dérivée d'ordre "n" par rapport au temps de l'écho ei(t), et, d'autre part, un circuit de multiplica- tion par C/2T de façon à retenir Viz(t) à partir de Ti(t) selon la formule (2). A titre d'exemple, on effectue ledit développement limité à l'ordre 1, le décalage temporel Ti(t) étant alors donné par: Ti(t) = [ei(t) - ei+1(t)] /ei(1 (t) (4) On remarque sur cette formule que l'on peut ainsi déterminer le signe de Ti(t) et donc le signe de Viz(t), ce qui permet d'obtenir le sens de la vitesse axiale. En ef- fet, un décalage te:nporel négatif représente un mouvement se rapprochant du transducteur, tandis qu'à l'inverse, un décala- ge positif est caractéristique d'un mouvement s'éloignant du transducteur. Par ailleurs, le procédé, selon l'invention, lié à l'utilisation de la formule (4) peut être mis en oeuvre de façon très simplifiée à l'aide d'un dispositif tel que ledit circuit d'extraction du décalage temporel Ti(t) est composé d'une ligne à retard, d'un soustracteur, d'un circuit de cal- cul de dérivées premières par rapport au temps et d'un divi- seur. Ce dispositif présente, en plus de sa simplicité, l'avantage de pouvoir être réalisé de façon entièrement analo- gique. Dans le mode particulier de réalisation relatif aux écoulements sanguins, le circuit d'estimation de la vites- se axiale est précédé par un circuit d'élimination des échos fixes et d'atténuation des échos liés à des mouvements lents, -- 6 -- comportant notamment un soustracteur de deux lignes échogra- phiques consécutives. Le signal résultant de cette différence est ensuite traité par le circuit d'estimation de la vitesse axiale. Comme on le verra plus loin en détail, la soustraction de deux lignes échographiques consécutives permet de s'affran- chir des échos rigoureusement fixes et de limiter les effets des échos correspondant à des mouvements lents. La description qui va suivre en regard des dessins annexés, donnés à titre d'exemples non limitatifs, fera bien comprendre en quoi consiste l'invention et comment elle peut être réalisée. La figure 1 montre un exemple de réalisation d'un dispositif de mise en oeuvre du procédé selon l'invention. La figure 2 représente un mode d'exécution parti- culier de l'étage d'émission du dispositif de la figure 1. La figure 3 est un mode de réalisation préféren- tiel d'un circuit d'estimation de la vitesse axiale. La figure 4 montre un schéma d'un circuit de sup- pression des échos fixes et d'atténuation des échos liés à des mouvements lents. La figure 5 est un diagramme donnant une loi de limitation des valeurs du décalage temporel Ti(t). La figure 6 représente un circuit discriminateur du dispositif de la figure 1. La figure 1 montre un schéma d'un dispositif de mise en oeuvre d'un procédé d'exploration par échographie ul- trasonore d'un milieu en mouvement, permettant de mesurer la vitesse axiale dudit mouvement, en projection sur l'axe "z" d'un faisceau d'excitations ultrasonores périodiques émis, à une période de récurrence "T", par un transducteur ultrasonore 10. Ce dispositif comprend, en outre, un étage d'émission 20, un étage 30 de réception et de traitement des signaux échogra- phiques renvoyés vers le transducteur 10, ainsi qu'un disposi- tif 40 de commande de balayage mécanique du transducteur. A la place de ce transducteur, pourrait bien entendu être uti- - 7' - lisé 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étail- lée sur la figure 2, l'étage d'émission 20 comprend un généra- teur 21 de signaux électriques d'excitation, envoyés vers le transducteur 10 qui les convertit en trains périodiques de si- gnaux 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 "Fe -de l'ordre de 5 kilohertz par exempledéterminée par un séquenceur comprenant successivement un oscillateur 22, ici de fréquence 32 méga-- hertz, et un diviseur de fréquence 23. Ce diviseur délivre les signaux d'horloge sur la connexion 102, ainsi que d'autres si- gnaux 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 comman- dent en particulier le dispositif 40, pour le balayage du transducteur. Un séparateur 24 des étages d'émission 2u 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'aveugle- ment 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 profon- deur 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. Comme on peut le voir à la figure 1, ledit étage 30 de réception et de traitement des signaux échographiques - 8 - comprend un circuit 330 d'estimation de la vitesse axiale per- mettant, à partir d'une analyse temporelle du signal, d'ex- traire dans un premier temps le décalage temporel Ti(t) in- duit par le mouvement du milieu entre deux échos successifs ei(t) et ei+ l(t) à l'aide d'un développement limité de la relation (1): ei+1(t) = ei[t - Ti(t)] reliant les deux échos, soit: p (-l)n (n) n ei+l(t) = E- ei (t).Ti (t) (3) n=O n! ei+l(t) pouvant être mesuré directement, et les dérivées ei(n t) pouvant être calculées-électroniquement, la relation (3) représente une équation en Yi(t) qu'il suffit de résou- dre pour en extraire le décalage temporel Ti(t) cherché. La vitesse axiale Viz(t) s'en déduit par application de la for- mule suivante: Viz(t) = Ti(t)C/2T (5) o C représente la vitesse de propagation de l'onde ultrasono- re. C'est pourquoi, conformément à la figure 3, le circuit 330 d'estimation de la vitesse axiale comporte, d'une part, un circuit 340 d'extraction du décalage temporel Ti(t) constitué pour résoudre l'équation (3), et d'autre part, un circuit 350 de multiplication par C/2T. En pratique, le décalage temporel Ti(t) est très faible vis-à-vis des temps caractéristiques de variation du signal échographique. Par exemple, l'application de la formule (2) avec: T = 200 ps, C = 1 500 m/s et Viz(t) = 5 cm/s (écoulements sanguins, mouvements des parois du coeur) conduit à: i(t) = 13,3 ns -9- Pour un signal échographique centré à 4 MHz, le décalage temporel est effectivement inférieur au dizième de la période du signal, soit 250 ns. En conséquence, un développe- ment limité au premier ordre de la relation (1) est justifié. L'équation (3) s'écrit alors avec n = 1: ei+1(t) = ei(t) - ei (1t)Ti(t) (6) Soit: ti(t) = [ei(t) - ei+l(t)]'/ei(1 (t) (4) Le circuit 340 d'extraction du décalage temporel, montré à la figure 3, représente un exemple, entièrement réa- lisé en analogique, de calcul de l'expression (6). A cet ef- fet, le circuit 340 est composé d'une ligne à retard 341, d'un soustracteur 342, d'un circuit 343 de calcul de la dérivée première par rapport au temps et d'un diviseur 344. Une ligne à retard possédant les qualités requises: retard important (200 ps), très stable (250 ps) avec une grande dynamique (79 dB) et une bande passante élevée (4 MHz) est décrite dan:- Le brevet français 78 014 62. Cependant, son implantation numérique est égale- ment concevable. Toutefois, le calcul de la dérivée nécessite quelques précautions. En effet, en échantillonnant le signal échographique au pas AT du signal 106 (60 ns dans l'exemple décrit), les points obtenus sont trop distants les uns des au- tres pour obtenir une estimée satisfaisante de la dérivée par la relation: ei(1)(kAT) = [ei((k+1)AT) - ei(kAT)]/AT La solution consiste à décaler, de façon analogi- que, l'écho ei(kAT) d'un temps 'e" petit devant AT (e = 5 ns par exemple), puis à soustraire ei(kT-e) de eikAT et, en- fin, à diviser la différence retenue par e, conformément à la formule: ei()(kAT) = [ei(kAT) - ei(kAT-c)]/e - 10 - En pratique, de façon à réduire le bruit d'estima- tion, il est avantageux de faire la moyenne du décalage tempo- rel Ti(t), calculé par la relation générale (3) et en-parti- culier par la relation (4), sur une fenêtre temporelle de lar- geur W afin d'évaluer un décalage temporel moyen Ti(t) défi- ni par: I t+W Ti(t): - $ti(u)du W t i Cette opération est effectuée par un circuit intégrateur-moyenneur 332 situé entre le circuit 340 d'extrac- tion du décalage temporel et le circuit 350 de multiplication par C/2T. On en déduit alors une vitesse axiale moyenne cor- respondante Viz(t) par: Viz(t) = Ti(t)C/2T Après avoir calculé Ti(t) et avant multiplica- tior. par C/2T. on peut également faire une moyenne de Ti(t) sur N tirs i consécutifs, à une cadence donnée par le signal 104, ce qui donne ici une moyenne sur N = 5 tirs. Cette moyen- ne est réalisée par le circuit 333 de la figure 3, en sortie de ce circuit on obtient donc: 1 N T(t) = k(t) N k=1 d'o la vitesse axiale: Vz(t) = T(t) C/2T Par ailleurs, afin de rester dans le domaine de validité du développement limité défini par la relation (6) qui suppose Ti(t) suffisamment petit, on a avantage à effec- tuer une limitation des valeurs de Ti(t), représentée par une fonction de distribution f(T) s'annulant lorsque T dépasse en valeur absolue une valeur maximale Tmax- La figure 5 mon- - il - tre un exemple d'une telle fonction de distribution: quand T est compris entre -Tmax et Tmax, f(T) = T et quand I T J>Tmax, f(T) = 0. Tmax peut être pris égal au dixième de la période P du signal échographique soit 25 ns pour P = 250 ns. Dans ce cas, le décalage temporel moyen Ti(t) est donné par: 1 t+W Ti(t) = tf[Ti(u)]du Dans ce but, le circuit 330 d'estimation de la vitesse axiale comporte un circuit 331 de limitation des valeurs du décalage temporel, situé entre le circuit 340 d'extraction du décalage temporel et le circuit intégrateur-moyenneur 332. Le signal de sortie du circuit 330 d'estimation de la vitesse axiale est alors validé ou non par un circuit dis- criminateur 360 montré à la figure 6, puis les valeurs ainsi confirmées sont envoyées vers le dispositif de vu.ualisation 312 par l'intermédiaire u dispositif de codage ii couleur 370. La présence du circuit discriminateur 360 est in- dispensable. En effet, si les lignes échographiques successi- ves, obtenues au rythme du déclenchement de l'excitation par le signal 102 à la fréquence F = 1/T, sont issues de cibles parfaitement fixes, le résultat de la différence entre ces deux lignes n'est que du bruit. Une ligne échographique s'é- crit, d'une manière générale: Yi(t) = f(t) + ei(t) o f(t) est le signal du aux cibles fixes et ei(t) l'écho produit par le milieu en mouvement. - 12 - La différence di(t) entre deux lianes consécuti- ves vaut donc: di(t) = vi+l(t) - Vi(t) = ei+1(t) - ei(t) (7) Si, Pendant un intervalle de temps donné, les li- qnes échoqraphiques vi(t) ne provient que de cibles fixes, on voit d'après (7) que di(t) = O au bruit près. Donc, le résultat fourni par le circuit 330 d'estimation de la vitesse axiale, aui traite ce bruit, n'est pas l'indication d'une vi- tesse nulle. et il devient nécessaire de valider ou non ce ré- sultat. Le circuit 360 comprend, à cet effet, en série, un multiDlieur 361 qui recoit le sianal de sortie di(t) du soustracteur 342 et élève ce siqnal de différence au carré, un sommateur 362 qui permet le calcul de l'énerqie locale sur une fenêtre de larqeur W' (éventuellement éqale à W) selon la for- mule: t+W' Liât) = t di2(u)du t Un circuit 364. 365 de calcul de moyenne aui. com- me dans le cas des circuits 344 et 345. est un accumulateur comprenant un additionneur 364 et une liane à retard 365 de retard T (ou un multiple de T). Permet d'effectuer la movenne de l'éneraie locale sur M tirs c'est-à-dire (M - 1) différen- ces selon l'expression: 1 i=M-1 E(t) = - E Ei(t) (.12) M - 1 i=1. avant d'envoyer la valeur ainsi obtenue vers un circuit de va- lidation comprenant un comparateur 461 recevant d'une part. sur une première entrée. la sortie de l'accumulateur 364. 365, ou directement celle du sommateur 362 dans le cas o le *- 13 - circuit de calcul de moyenne ne serait pas prévu, et d'autre part, sur une deuxième entrée 462, une tension de référence formant seuil, la sortie du comparateur étant ici au niveau logique O ou 1 selon que la tension reçue sur sa première en- trée est respectivement inférieure ou supérieure au seuil de référence "N(t) proportionnel au niveau de bruit N(t). Un mul- tiplieur 463, qui reçoit le signal de sortie VZ(t) du cir- cuit 330 sur une première entrée, transmet ce signal, noté dé- sormais V'z(t), sur une sortie ou transmet simplement des valeurs nulles selon que le signal de validation délivré sur une deuxième entrée par le comparateur 461 est 1 ou 0 respec- tivement. En effet, en dehors des zones d'écoulements vrais, l'énergie moyenne calculée en sortie du circuit 364, 365 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 appro- priée du seuil, c'est-à-dire que N(t) est également donné par: t+w' [ d2(u)du t i en dehors de toute excitation. Le niveau effectif du seuil est donc déterminé par le coefficient a, au choix de l'utilisa- teur. Au contraire, en présence de signaux rétrodiffusés par les cibles en mouvement, l'énergie moyenne du signal di(t) est supérieure à celle du bruit seul, ce qui autorise la vali- dation des signaux délivrés par le circuit 330 d'estimation de la vitesse axiale. Remarquons que la valeur de a peut servir aussi, à la visualisation, à fixer une frontière entre un affichage couleur et un affichage gris; si E(t) est supérieur au seuil "N(t), l'affichage sera en couleur. Par contre, il sera grissi E(t) est inférieur au seuil. La sortie du circuit discriminateur 360 est envoyée vers le dispositif de mémorisation, conversion de -14- balaYaqe et codaqe de couleur 370, par lequel transite éaale- ment, avant visualisation, le siqnal de sortie de l'amplifica- teur 311 de la voie de traitement 301. Un tel dispositif est décrit, par exemple, dans la demande de brevet européen EP-A- 0 100 094. La ficure 3 de ce document montre en effet, entre les bornes A, B, C, et ER, EG. EB. un exemmle des cir- cuits aui Peuvent être mis en place dans le cas de l'inven- tion, la borne A recevant le signal échoaraphique classique et les bornes B et C les paramètres caractéristiques du milieu en mouvement examiné. Ce dispositif 370 et le dispositif 312 per- mettent alors la visualisation en temps réel des écoulements ou des déplacements superposés à l'image de réflectivité échoqraphique classique. La figure 4 donne le schéma d'un circuit 320 de suppression des échos fixes et d'atténuation des échos liés à des mouvements lents, particulièrement nécessaire dans le cas d'étude d'écoulements sanquins. Le circuit numérique 320 de suppression des éck ', représenté sur la figure 4, comprend lui-même, dans cet exem- ple de réalisation, 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 pour- rait être de plusieurs Périodes T = 1/F, mais il est avanta- aeux que ce retard soit le plus faible possible et égal à T. Le soustracteur 322 effectue donc la différence di(t) entre deux liGnes échoqraphiaues successives vi(t) et Yi+1(t). C'est pourquoi, lorsque ce circuit 320 est pré- sent, sa sortie peut être directement connectée à l'entrée du circuit discriminateur 360 qui exige aussi la connaissance de di(t). D'autre part, la relation de récurrence: ei+l(t) = eift - i(t)], - 15 - étant éqalement vérifiée par di(t), le sianal de sortie du circuit 320 sert de siqnal d'entrée au circuit 330 d'estima- tion de la vitesse axiale. 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 étu- diés. Cette présence d'échos fixes est qênante du fait de leur amplitude bien plus élevée (de l'ordre de +40 dB dans le cas des écoulements sanquins) que celle des siqnaux utiles, les siqnaux 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 siqnal de commande d'échantillonnaqe à la fréquence de 16 MHz. - 16 - REVENDICATIONS: 1. Dispositif d'exploration par échographie ultraso- nore d'un milieu en mouvement, permettant de mesurer la vites- se axiale dudit mouvement, en projection sur l'axe "z' d'un faisceau d'excitations ultrasonores périodiques émis, à une période de récurrence "T", par au moins un transducteur ultrasonore (10), comprenant, en outre, un étage d'émission (20) et un étage (30) de réception et de traitement des signaux échographiques renvoyés vers le transducteur (10), caractérisé en ce que ledit étage (30) de réception et de traitement des signaux échographiques comprend notamment un circuit (330) d'estimation de ladite vitesse axiale Viz(t) comportant, d'une part, un circuit (340) d'extraction du décalage temporel Ti(t) entre deux échos successifs ei(t) et ei+l(t), constitué pour résoudre l'équation en Ti(t): p (-1)n (n) n ei+ l(t) = - ei (t).Ti (t) n=O n! o ei(n)(t) est la dérivée d'ordre "n" par rapport au temps de l'écho ei(t), et, d'autre part, un circuit (350) de multiplication par C/2T de façon à obtenir Viz(t) à partir de Ti(t), C étant la vitesse de propagation de l'onde ultra- sonore. 2. Dispositif, selon la revendication 1, caractérisé en ce que, l'ordre de dérivation "p" étant égal à 1, ledit circuit (340) d'extraction du décalage temporel li(t) est composé d'une ligne à retard (341), d'un soustracteur (342), d'un circuit (343) de calcul de dérivée première par rapport au temps et d'un diviseur (344). 3. Dispositif, selon la revendication 2, caractérisé en ce que, le signal échographique étant numérisé par échan- tillonnage au pas AT, le circuit de calcul de la dérivée pre- - 17 - mière par rapport au temps comporte une ligne à retard déca- lant l'écho ei(kAT) d'un temps e petit devant AT, un soustracteur réalisant la différence ei(kAT) - ei(kAT - e), et un diviseur par c de façon à obtenir la dérivée ei(1)(kAT). 4. Dispositif, selon l'une des revendications 1 à 3, caractérisé en ce que ledit circuit (330) d'estimation de la vitesse axiale Viz(t) comporte également un circuit intégra- teur-moyenneur (332) de Ti(t) sur'une fenêtre de largeur W, situé entre le circuit (340) d'extraction du décalage tempo- rel et le circuit (350) de multiplication par C/2T. 5. Dispositif, selon la revendication 4, caractérisé en ce que le circuit (330) d'estimation de la vitesse axiale Viz(t) comporte, en outre, un circuit (331) de limitation des valeurs du décalage temporel Ti(t), situé entre-le cir- cuit (340) d'extraction du décalage temporel et le circuit intégrateur-moyenneur (ian). 6. Dispositif, -2lon l'une des revendications 1 à 5, caractérisé en ce que le circuit d'estimation de la vitesse axiale est suivi d'un circuit discriminateur (360) comprenant un circuit (361) d'élévation au carré de la différence di(t) = ei(t) - ei+l(t) entre deux échos successifs, un circuit sommateur (362) sur une fenêtre de largeur W' donnant: t+W' E(t) = S di2(u)du t et un circuit de validation composé d'un comparateur de E(t) à un seuil aN(t) proportionnel au niveau de bruit N(t) et d'un multiplicateur (463) des valeurs de Viz(t) par 1 ou O suivant que E(t) est supérieur ou inférieur à aN(t). 7. Dispositif, selon l'une des revendications 1 à 6, caractérisé en ce que le circuit d'estimation (330) de la vitesse axiale est précédée par un circuit d'élimination (320) 260408_1 - 18 - des échos fixes et d'atténuation des échos liés à des mouvements lents, comportant notamment un soustracteur (322) de deux lignes échographiques consécutives.