FR2906612A1 - Procede et dispositif de tomographie par impedance electrique. - Google Patents

Abstract

Procédé de tomographie par impédance électrique comprenant :- une étape de mesure électrique au cours de laquelle on impose des conditions électriques prédéterminées à la surface (3) d'un milieu à imager, tout générant une perturbation mécanique en des points prédéfinis du milieu pour modifier localement l'impédance du milieu, et on mesure un paramètre électrique en plusieurs points (7) à la surface du milieu,- et une étape de calcul au cours de laquelle on détermine l'impédance électrique en plusieurs points dans le volume intérieur (2) du milieu, en prenant en compte les mesures effectuées pendant la perturbation, en fonction d'une loi de modification de l'impédance électrique par cette perturbation.

Description

1 Procédé et dispositif de tomographie par impédance électrique. La

présente invention est relative aux procédés et dispositifs de tomographie par impédance électrique. Plus particulièrement, l'invention concerne un procédé de tomographie par impédance électrique pour imager un milieu présentant un certain volume intérieur délimité par une surface extérieure, ce procédé comprenant : - au moins une étape de mesure électrique au cours de laquelle on impose des conditions électriques prédéterminées à la surface du milieu et on mesure au moins un paramètre électrique en plusieurs points à la surface du milieu, - et au moins une étape de calcul au cours de laquelle on détermine au moins un paramètre relié à l'impédance électrique en plusieurs points dans le volume intérieur du milieu. Le document US-A-4 291 708 décrit un exemple d'un 20 tel procédé. Les procédés connus de tomographie par impédance électrique souffrent d'une résolution médiocre (de l'ordre de quelques centimètres pour les applications médicales ou de l'ordre de quelques dizaines de mètres dans le domaine 25 de la géophysique) et dépendant de la profondeur par rapport à la surface du milieu à imager. La présente invention a notamment pour but de pallier ces inconvénients. A cet effet, selon l'invention, un procédé du genre 30 en question est caractérisé en ce qu'au cours d'au moins certaines étapes de mesure électrique, on génère une perturbation mécanique en des points prédéfinis du milieu en modifiant localement l'impédance du milieu, et en ce qu'au cours de l'étape de calcul, on détermine 35 ledit paramètre relié à l'impédance électrique en prenant 2906612 2 en compte les mesures effectuées pendant ladite perturbation, en fonction d'une loi de modification de l'impédance électrique par ladite perturbation. Grâce à ces dispositions, on peut accroître 5 considérablement la précision et la rapidité de mise en oeuvre du procédé selon l'invention, grâce au fait que la perturbation susmentionnée effectue un "marquage" local du milieu. En particulier, on peut obtenir des précisions millimétriques, y compris jusqu'à des profondeurs assez 10 grandes par rapport à la surface extérieure du milieu, et ce pour un coût de mise en oeuvre très modique. Dans divers modes de réalisation du procédé selon l'invention, on peut éventuellement avoir recours en outre à l'une et/ou à l'autre des dispositions suivantes : 15 - les conditions électriques imposées comprennent au moins un courant imposé en au moins un point à la surface du milieu, et ledit paramètre électrique mesuré est un potentiel électrique (il est bien entendu possible d'imposer un potentiel et de mesurer des courants) ; - ledit paramètre relié à l'impédance électrique est la conductivité ; - la perturbation mécanique est une onde focalisée en au moins un point du milieu ; - l'onde est une onde acoustique ; - l'onde acoustique est une onde ultrasonore (générée par exemple par un ensemble de transducteurs piézo-électriques) ; - l'onde acoustique correspond à un signal modulé en amplitude à une fréquence de modulation adaptée pour générer une force de radiation ultrasonore entraînant un déplacement local du milieu ; - l'onde est une onde élastique (mélange d'ondes de compression et d'ondes de cisaillement, généré par exemple par un ensemble de vibreurs mécaniques disposés à la surface du sol) ; 2906612 3 - l'onde correspond à un signal codé ; - au cours de l'étape de calcul, on résout l'équation : (8) y(z) = Ek(z) AVuk(z)•Vuk(z) 5 pour tout point z du milieu à imager, où . k est un indice désignant un ensemble d'au moins un courant électrique j,k appliqué à la surface du milieu, i étant un indice désignant chaque courant électrique de cet 10 ensemble, Ek(z) est une énergie correspondant à la perturbation engendrée par l'onde pendant l'application de l'ensemble de courants électriques j,k, uk(z) est le potentiel électrique au point z du 15 milieu, ^ et A est une matrice représentative de la forme d'une tache focale produite par l'onde autour du point sur lequel elle est focalisée ; - la tache focale est sphérique et A est la matrice 20 identité ; - au cours de l'étape de calcul, on détermine ladite énergie en fonction de ladite loi de modification de l'impédance électrique par ladite perturbation ; - ladite énergie est calculée par la formule : 25 (6) Ek(z)=EDk(z)Jk, où . . z désigne un point situé dans le milieu, . Dk(z) est une valeur représentative de la perturbation électrique mesurée en un point d'indice i à la surface du 30 milieu et générée par l'onde pendant l'application de l'ensemble de courants électriques j,k aux points i ; - Dk(z) correspond à la formule suivante : 2906612 4 (3) Dk(z)=yAVuk(z)•VG(z,i), où y est la conductivité et G(z,i) la fonction de Neumann du milieu ; - les valeurs D,k(z) sont calculées à partir des 5 mesures effectuées, en utilisant des ondes correspondant à différents signaux S1(t), 1 étant un indice compris entre 1 et L ; - L est égal à 2 et on utilise deux signaux S1(t) = S1.s(t) et S2(t) = S2.s(t) d'amplitudes respectives 10 S1 et S2, les valeurs Dk(z) étant calculées, lorsque la forme de la tache focale est un disque ou une sphère, par la formule : Dk(z) - (ut" (z)-uk)(u,2 k (z)-uk)(s, -s2) (4) d(s2uk.l(z)-S2uk -,Sluk'2(z)+S1uk)vl. où 15 ^ d est soit égal à 2 pour une imagerie à 2 dimensions, soit égal à 3 pour une imagerie à 3 dimensions, • IVI est soit la surface de la tache focale pour une imagerie à 2 dimensions, soit le volume de la tache focale pour une imagerie à 3 dimensions ; 20 -l'onde est une onde ultrasonore, L est égal à 2 et on utilise deux signaux S1(t) = Sl.s(t) et S2(t) = S2.s(t) d'amplitudes respectives S1 et S2, s(t) étant un signal modulé en amplitude à une fréquence de modulation adaptée pour générer une force de radiation 25 ultrasonore entraînant un déplacement local du milieu, les valeurs Dk(z) étant calculées, lorsque la forme de la tache focale est un disque ou une sphère, par la formule : ( k,l k)(u2 kk)( 2 2) Dk (z)= ui (z) ù ui 1 (z) ù "1 `S1 _s 2 (4) J d(S22uk.1 (z) - s22uk -s 2u".2(z) + s12uk )v où . 30 d est soit égal à 2 pour une imagerie à 2 dimensions, 2906612 5 soit égal à 3 pour une imagerie à 3 dimensions, . IVI est soit la surface de la tache focale pour une imagerie à 2 dimensions, soit le volume de la tache focale pour une imagerie à 3 dimensions ; 5 au cours de l'étape de calcul, on part d'une conductivité supposée y et on répète les sous-étapes suivantes . a) on résout numériquement l'équation suivante : fdiv(yvuk )= 0 en tout point z du milieu (9) yvuk • n = jk sur la surface extérieure 10 y étant la valeur précédemment estimée de la conductivité (initialement, y est donc la valeur supposée susmentionnée), b) on calcule une erreur estimée ek sur la conductivité, c) on résout l'équation suivante : div(yvvk )= ûdiv(ekVuk ) en tout point z du milieu, (11) yvuk • n = 0 sur la surface extérieure. où vk est la solution de l'équation (9) et uk la solution de l'équation (11), d) on met à jour la conductivité par la formule : 15 (12) Yk =ûY i 1 2A,vuk.vvk A2Vuk 2 +ek, 20 où vk est la solution de l'équation (11) et uk la solution de l'équation (9), 2906612 6 et on utilise Yk comme nouvelle valeur de conductivité y avec un autre ensemble de courants jk, générant des lignes de courant non parallèles à celles générées par l'ensemble de courants j,k, dans au moins une zone du 5 milieu, les sous-étapes a) à d) étant réitérées jusqu'à satisfaction d'un critère d'arrêt ; - au cours de la sous-étape b), on calcule une erreur estimée ek sur la conductivité par la formule : 10 (10) ek =Ek/AVuk •Quk ùy ; -au cours de la sous-étape b), on calcule une erreur estimée ek sur la conductivité par la formule : (10 ) ek = Ek / ADuk ' Quk - 15 4'k /AVuk • Duk ù y + 1) - au cours de la sous-étape (d) pour chaque point z du milieu, on cherche quel indice k de conditions électriques correspond à la plus 2 20 grande énergie A2Vuk ce qui donne une fonction k(z), on met à jour la conductivité comme suit : Yk(,) - -Y 2 A2V uk(-) • Vvk(-) 2 A2Vuk(_) et on utilise yk(z) comme valeur de conductivité y; - le milieu à imager est un tissu biologique ; 25 - le milieu à imager est un organe humain (par exemple : sein, foie, cerveau, ou autres) ; - le milieu à imager est le sous-sol terrestre.

2906612 7 Par ailleurs, l'invention a également pour objet un dispositif adapté pour la mise en oeuvre d'un procédé tel que défini ci-dessus. D'autres caractéristiques et avantages de 5 l'invention apparaîtront au cours de la description suivante d'un de ses modes de réalisation, donné à titre d'exemple non limitatif, en regard des dessins joints. Sur les dessins : - la figure 1 est une vue schématique d'un 10 dispositif de tomographie d'impédance électrique selon une forme de réalisation de l'invention, - la figure 2 est un graphe représentant le signal correspondant aux ondes ultrasonores appliquées par le dispositif de la figure 1, 15 - et la figure 3 est un graphe similaire à la figure 1, pour une variante de l'invention. Sur les différentes figures, les mêmes références désignent des éléments identiques ou similaires. La figure 1 montre un dispositif 1 de tomographie 20 d'impédance électrique, utilisé pour imager un milieu (par exemple un milieu biologique, par exemple une partie d'un corps humain telle qu'un sein, un foie, un cerveau ou tout autre organe) qui présente un certain volume intérieur 2 délimité par une surface extérieure 3.

25 Le dispositif 1 comporte une unité centrale 4 (UC) telle qu'un ordinateur ou similaire, qui peut être relié à divers périphériques tels qu'un écran 5 et d'autres interfaces d'entrée et sortie (non représentées). L'unité centrale 4 est reliée à une interface de 30 mesure électrique 6 (INT. 1) telles que celles utilisées classiquement en tomographie d'impédance électrique, reliée à une pluralité d'électrodes 7 disposées sur la surface 3 du milieu à imager. Les électrodes 7 sont au nombre de I (entier naturel non nul) et sont désignées chacune par un 35 indice i compris entre 1 et I.

2906612 8 L'interface de mesure électrique 6 est adaptée pour imposer des conditions électriques prédéterminées sur certaines des électrodes 7 et pour mesurer au moins un paramètre électrique au niveau de tout ou partie des 5 électrodes 7. Par exemple, l'interface de mesure électrique 6 peut être adaptée pour imposer au moins un courant prédéterminé j au niveau d'une des électrodes 7 et pour mesurer des tensions ui au niveau de toutes les électrodes 7. Plus particulièrement, l'unité centrale 4 est adaptée 10 pour commander l'interface de mesure électrique 6 pour qu'elle impose successivement plusieurs courants jl, ... Jk, jK (k étant un indice désignant chaque mesure et K un nombre entier naturel non nul désignant le nombre total de mesures avec des courants différents chaque courant 15 d'indice k peut différer des autres courants par l'électrode 7 à laquelle il est appliqué et/ou par le signal auquel il correspond) et pour mesurer les tensions ui pour chaque courant jk. Alternativement, l'interface de mesure électrique 6 20 pourrait être adaptée pour imposer une ou plusieurs tensions et mesurer des courants au niveau des électrodes 7. Les courants et tensions en question peuvent par exemple être alternatifs, par exemple avec une fréquence de 25 l'ordre du kHz. Par ailleurs, l'unité centrale 4 est reliée à une interface de traitement de signal 8 (INT. 2) qui commande un réseau 9 de transducteurs piézo-électriques ultrasonores (par exemple une barrette rectiligne de transducteurs) 30 appliqué sur la surface 3 du milieu à imager. On notera que le réseau 9 de transducteurs peut également être un réseau à 2 dimensions et/ou être monté sur un support mobile permettant de faire varier la position et/ou l'orientation du réseau.

35 L'interface de traitement de signal 8 est adaptée 2906612 9 pour faire générer, par le réseau 9 de transducteurs, des ondes ultrasonores focalisées successivement en des points prédéterminés z (au nombre de Z) situés dans le milieu à imager, au moins certains de ces points z étant situés dans 5 le volume intérieur 2 (les autres étant éventuellement sur la surface 3). Les ondes ultrasonores en question peuvent par exemple de fréquence comprise entre 0,5 et 15 MHz, notamment de l'ordre du MHz. Lorsqu'on souhaite imager le volume intérieur 2 du 10 milieu, l'unité centrale 4 fait appliquer successivement, par l'interface de mesure électrique 6, des courants jk sur tout ou partie des électrodes 7. Pour chaque courant l'interface de mesure électrique 6 mesure la tension (uk)1<,<1 de chaque électrode 7 15 d'indice i, en l'absence d'onde acoustique ultrasonore dans le milieu à imager. Puis l'unité centrale 4 fait émettre, par l'interface 8 de traitement de signal et le réseau de transducteurs 9, des ondes ultrasonores focalisées 20 successivement sur les points prédéterminés z susmentionnés, de façon à générer des perturbations mécaniques du milieu localisées en chaque point z, se traduisant par des perturbations localisées de l'impédance électrique (et notamment de la conductivité) du milieu.

25 L'onde ultrasonore est focalisée sur chaque point z pendant quelques centaines de périodes de l'onde ultrasonore. L'onde ultrasonore en question peut être une onde non modulée à basse fréquence, qui induit une variation 30 locale et infinitésimale de volume dans la zone focale du faisceau ultrasonore. La fréquence à laquelle se produit cette vibration est la fréquence d'excitation des ultrasons. On produit alors des perturbations locales de l'impédance électrique ayant la même fréquence que le 35 signal ultrasonore. Pour augmenter le rapport signal sur 2906612 10 bruit nécessaire pour détecter l'influence du faisceau ultrasonore sur les signaux électriques, on peut par exemple émettre un signal ultrasonore codé S(t), tel que par exemple celui représenté sur la figure 2.

5 Comme signal codé, on peut utiliser par exemple une fonction chirp S(t) = sin (211 (f0 + tf. t) t) où f0 est une fréquence et Af une largeur de bande en fréquence. En variante, on peut également utiliser à titre de codage une réalisation prédéterminée d'un bruit blanc.

10 Plus généralement, on peut faire émettre successivement des ondes ultrasonores correspondant à des signaux S1 (t) , ... S1 (t) , ... SL (t) , 1 étant un exposant compris entre 1 et L, L étant un entier naturel non nul (L = 1 si on n'utilise qu'un seul signal S(t)). Dans ce 15 cas, les ondes correspondant à chaque signal SI peuvent par exemple être focalisées successivement sur les différents points z, avant d'émettre et de focaliser les ondes correspondant au signal S1+1. On notera que les signaux s1(t) peuvent éventuellement différer entre eux uniquement 20 par leur amplitude. Pendant que l'on applique le courant jk et le signal SI et que l'onde ultrasonore correspondant au signal S1 est focalisée sur un point donné z, on mesure les tensions (uk" L aux I électrodes 7 d'indices i. Si le 25 signal S1(t) est codé comme indiqué ci-dessus, la sensibilité de la mesure peut être améliorée par k,l,z déconvolution du signal électrique u I<<1 par le code appliqué au signal S1(t). Au total, l'unité centrale mémorise les I.K mesures 30 de tension k u lq<1,I<k<K faites sans focalisation d'ondes acoustiques et les I . K. L. Z mesures de tension uk,!,_l JI<,<l,l<k<K,I<l<l. faites avec focalisation d'ondes acoustiques. A partir de ces mesures, on effectue un calcul de reconstruction consistant à retrouver la conductivité y(z) 2906612 11 en tout point z du milieu à imager. Cette reconstruction correspond au problème suivant . Jdiv(y(z)Vuk (z))= 0 en tout point z du milieu (1) y(z)Vuk • n = jk sur la surface extérieure 5 où uk est la tension électrique (potentiel) sur la surface extérieure 3, uk(z) est la tension électrique au point z dans le volume intérieur 2 et est la normale à la surface extérieure 3 du milieu. Pour résoudre ce problème mathématique, on peut 10 éventuellement utiliser une méthode standard consistant à tester différentes conductivités en minimisant l'écart aux données mesurées, par exemple par la méthode des moindres carrés. Plus avantageusement, on pourra utiliser la méthode 15 décrite ci-dessous, qui s'est avérée particulièrement avantageuse, précise et robuste. Cette méthode part du fait que, comme enseigné par les travaux de H. Ammari et H. Kang ( Reconstruction of Small Inhomogenei_ties from Boundary Measurements, Lecture 20 notes in mathemetics, volume 1846, Springer Verlag, Berlin, 2004), la perturbation électrique due à un changement de conductivité locale en un point z du milieu à imager est donnée au premier ordre par la formule : (2) uk!_ _uk =(yp y)m]uk(z)•VG(z,i), 25 où . - uk est le potentiel au point i et e est le potentiel dû à la perturbation ultrasonore qui correspond au signal 1, focalisée au point z, - Yp est la conductivité perturbée localement par 30 les ultrasons, - M est un facteur géométrique, le tenseur de polarisation, qui dépend de Yp et de la forme de la zone focale de l'onde ultrasonore (par exemple, la zone dans 2906612 12 laquelle l'amplitude de l'onde ultrasonore est supérieure à la moitié de l'amplitude maximale), - la fonction G est la fonction de Neumann du milieu de conductivité y et donc, inconnue.

5 On peut extraire de ces perturbations électriques, une matrice D représentative de la perturbation électrique au point z en présence du courant k: (3) Dk(z)=rAVuk(z)•VG(z,i), où A est une matrice définie positive connue dépendant 10 uniquement de la forme de la zone focale. Cette matrice A est reliée au tenseur de polarisation M par la formule A = 1/d(y/yp+(d-1))M, où d est la dimension de l'espace (d = 2 pour une imagerie à 2 dimensions et d = 3 pour une imagerie à 3 dimensions).

15 La forme de la zone focale de l'onde ultrasonore étant connue, cette matrice A est connue. Par exemple, pour une zone focale circulaire ou sphérique, A est égale à la matrice identité Id. La matrice D peut être calculée à partir des 20 mesures effectuées, en utilisant seulement les différences d'amplitudes des ondes ultrasonores correspondant à différents signaux S1(t). Par exemple, à partir de deux signaux S1(t) _ S1. s (t) et S2(t) = S2. s (t) d'amplitudes respectives S1 et S2, 25 la matrice en question peut être calculée, lorsque la forme de la tache focale est un disque ou une sphère, par la formule : k (uk>~(z)-uk/(uk2(z)-uk)\SI ù`s2) (4 ) Di (z~ = d(S2u" (z) - S2uk - Sluk.2 (z~ + Siuk )v . où . 30 - d est la dimension de l'espace (2 pour une imagerie à 2 dimensions et 3 pour une imagerie à 3 dimensions, - 1VI est la surface (pour une imagerie à 2 2906612 13 dimensions) ou le volume (pour une imagerie à 3 dimensions) de la zone focale ultrasonore, et S1 et S2 les amplitudes des ondes ultrasonores. Cette matrice D permet de calculer l'énergie 5 électrique Ek(z) équivalente aux perturbations acoustiques aux points z, pour chaque courant k. Cette énergie électrique est définie par la formule : 10 (5) Ek(z)=y(z)AVuk(z)•Duk(z), et calculée en pratique par la formule : Ek(z) = JDk (z, y)J k (y)dy (6) = E Dk (z)J,k On peut aussi utiliser d'autres formules de quadrature, par exemple lorsque les courants ne sont pas 15 mesurés au niveau des électrodes i mais en des points distincts à la surface extérieure 3 du milieu à imager. En reportant la formule (6) dans la formule (1), le problème mathématique à résoudre peut s'écrire comme suit : (7) ( k div k k Vuk = 0 en tout point z du volume intérieur 2 du milieu à imager Avu •vu yvuk n = jk sur la surface extérieure 3 du milieu à imager On notera que, le cas échéant, la conductivité y sur la surface extérieure 3 du milieu à imager peut être 25 déterminée au moyen de mesures indépendantes, sans recourir à cette méthode. On résout cette équation non linéaire qui ne comporte aucun coefficient inconnu ce qui donne uk(z) en tout point z du milieu à imager. On obtient ainsi 20 10 15 2906612 14 (8) y(z) = Ek (z) AVuk(z) •Duk (z) pour tout point z du milieu à imager. Pour résoudre cette équation non linéaire, on peut utiliser l'algorithme suivant: 5 1/ on part d'une conductivité supposée y par exemple y=l en tout point du milieu. 2/ on répète les étapes suivantes : a) on résout numériquement, avec un solveur linéaire standard pour un courant jk le problème suivant Jdiv(yVuk )= 0 en tout point z du milieu (9) yDuk .n= jk sur la surface extérieure b) on calcule l'erreur (10) ek = Ek / AVuk • Vuk -y , Ekétant l'énergie calculée par la formule (6) ; c) on résout le problème div(yvvk )= -div(ekVuk ) en tout point z du milieu, (11) - ; yVuk • n = 0 sur la surface extérieure. d) on met à jour la conductivité par la formule (12) Yk =-Y i A2Vuk .Ovk 2 2 +ek , A2Vuk où vk est un potentiel électrique solution de l'équation (11) et Uk un potentiel électrique solution de l'équation 20 (9), et on utilise yk donnée par la formule (12) comme valeur de conductivité y à l'itération suivante, pour un ensemble de 2906612 15 courants jk, approprié (générant des lignes de courant non parallèles à celles générées par jk, au moins dans certaines zones du milieu à imager). On réitère les étapes a) à d) jusqu'à satisfaction 5 d'un critère d'arrêt, par exemple : - lorsqu'une norme de l'erreur e devient petite, - ou lorsqu'une norme de Vvk devient petite. En pratique, une dizaine d'itérations des étapes a) à d) suffisent pour converger.

10 En variante, on peut utiliser des ondes ultrasonores correspondant à des signaux S1(t) modulés à relativement basse fréquence, par exemple avec une fréquence de modulation de quelques centaines de Hz, comme représenté sur la figure 3. Dans ce cas, le faisceau 15 ultrasonore induit une force localisée dans la zone focale qui pousse le milieu. Cette force, connue sous le nom de force de radiation ultrasonore, engendre un déplacement local du milieu dont les variations temporelles sont liées non pas à la fréquence ultrasonore, mais à l'enveloppe 20 basse fréquence du signal ultrasonore. Lors de l'utilisation de cette technique n 2, la force de radiation ultrasonore pourra elle aussi être codée dans le temps, en modulant dans le temps l'amplitude du signal ultrasonore de poussée.

25 Dans ce cas, on peut utiliser le même procédé que celui décrit ci-dessus, de préférence en remplaçant la formule (4) par la formule (4') suivante , lorsque la forme de la tache focale est un disque ou une sphère : (u"(z) )(u(z) )( 2 2) k _ - Zli - Zti iJ 1 - `J 2 (4') D` (z) d S 2uk'' Z S 2uk -S 2uk.2 z + S 2uk V 30 Par ailleurs, pour rendre l'algorithme plus stable et robuste, on peut remplacer la formule de l'erreur (10) par 2906612 16 (10') ek =(Ek / AVuk •Vuk ) ,Ek l AVuk • ouk -y +1) Dans ce cas, la convergence est plus lente, mais de très forts contrastes dans les matériaux sont détectables 5 sans instabilité. Par ailleurs, on peut aussi utiliser les données de plusieurs courants simultanément dans le procédé ci-dessus. A cet effet, au cours de l'étape (d) pour chaque point z du domaine on cherche quel 2 10 courant k correspond à la plus grande énergie A2Vuk ce qui donne une fonction k(z), pour chaque courant k(z) on met à jour la conductivité comme suit : Yk(:) = ûY A2Duk( ) •Ovk(z) 2 2 +ek(') , y étant la valeur A2Vuk(') 15 précédente de la conductivité, et on utilise Yk(Z) comme valeur de conductivité y à l'itération suivante. On notera que le procédé selon l'invention peut également être utilisé en géophysique. Dans ce cas, le 20 milieu à imager est le sol terrestre et les ondes ultrasonores sont remplacées par des ondes élastiques, notamment des ondes de compression ou de cisaillement à basse fréquence (par exemple, de 5 Hz à 5000 Hz). Les transducteurs piézo-électriques susmentionnés sont alors 25 remplacés par un ensemble de vibreurs ou vérins mécaniques disposés à la surface du sol.

Claims (25)

REVENDICATIONS

1. Procédé de tomographie par impédance électrique pour imager un milieu présentant un certain volume intérieur (2) délimité par une surface extérieure (3), ce procédé comprenant : - au moins une étape de mesure électrique au cours de laquelle on impose des conditions électriques prédéterminées à la surface (3) du milieu et on mesure au moins un paramètre électrique en plusieurs points (7) à la surface (3) du milieu, -et au moins une étape de calcul au cours de laquelle on détermine au moins un paramètre relié à l'impédance électrique en plusieurs points dans le volume intérieur (2) du milieu, caractérisé en ce qu'au cours d'au moins certaines étapes de mesure électrique, on génère une perturbation mécanique en des points prédéfinis du milieu en modifiant localement l'impédance du milieu, et en ce qu'au cours de l'étape de calcul, on détermine ledit paramètre relié à l'impédance électrique en prenant en compte les mesures effectuées pendant ladite perturbation, en fonction d'une loi de modification de l'impédance électrique par ladite perturbation.

2. Procédé selon la revendication 1, dans lequel les conditions électriques imposées comprennent au moins un courant imposé en au moins un point (7) à la surface du milieu, et ledit paramètre électrique mesuré est un potentiel électrique.

3. Procédé selon la revendication 1 ou la revendication 2, dans lequel ledit paramètre relié à l'impédance électrique est la conductivité.

4. Procédé selon l'une quelconque des revendications précédentes, dans lequel la perturbation 35 mécanique est une onde focalisée en au moins un point du 2906612 milieu.

5. Procédé selon la revendication 4, dans lequel l'onde est une onde acoustique.

6. Procédé selon la revendication 5, dans lequel 5 l'onde acoustique est une onde ultrasonore.

7. Procédé selon la revendication 6, dans lequel l'onde correspond à un signal modulé en amplitude à une fréquence de modulation adaptée pour générer une force de radiation ultrasonore entraînant un déplacement local du 10 milieu.

8. Procédé selon la revendication 4, dans lequel l'onde est une onde élastique.

9. Procédé selon l'une quelconque des revendications 4 à 8, dans lequel l'onde correspond à un 15 signal codé.

10. Procédé selon l'une quelconque des revendications 4 à 9, dans lequel, au cours de l'étape de calcul, on résout l'équation : (8) Y(z) = Ek (z) AVuk(z)•Vuk(z) 20 pour tout point z du milieu à imager, où . k est un indice désignant un ensemble d'au moins un courant électrique jk à la surface du milieu, i étant un indice désignant chaque courant de cet ensemble, 25 ^ Ek(z) est une énergie correspondant à la perturbation engendrée par l'onde pendant l'application de l'ensemble de courants électriques jk à la surface du milieu, uk(z) est le potentiel électrique au point z du milieu, 30 • et A est une matrice représentative de la forme d'une tache focale produite par l'onde autour du point sur lequel elle est focalisée.

11. Procédé selon la revendication 10, dans lequel la tache focale est sphérique et A est la matrice identité. 18 2906612 19

12. Procédé selon la revendication 10 ou la revendication 11, dans lequel au cours de l'étape de calcul, on détermine ladite énergie en fonction de ladite loi de modification de l'impédance électrique par ladite 5 perturbation.

13. Procédé selon l'une quelconque des revendications 10 à 12, dans lequel ladite énergie est calculée par la formule : (\ / \ (6) Ek(z)=ED,k(z)jk 10 où . . z désigne un point situé dans le milieu, . D,k(z) est une valeur représentative de la perturbation électrique mesurée en un point d'indice i à la surface (3) du milieu et générée par l'onde pendant l'application de 15 l'ensemble de courants électriques j,k aux points i.

14. Procédé selon la revendication 13, dans lequel Dk(z) correspond à la formule suivante : (3) Dk(z)=yAVuk(z)•VG(z,i), où y est la conductivité et G(z,i) la fonction de Neumann 20 du milieu.

15. Procédé selon la revendication 13 ou la revendication 14, dans lequel les valeurs Dk(z) sont calculées à partir des mesures effectuées, en utilisant des ondes correspondant à différents signaux S1(t), 1 étant un 25 indice compris entre 1 et L.

16. Procédé selon la revendication 15, dans lequel L est égal à 2 et on utilise deux signaux S1(t) = S1.s(t) et S2 (t) = S2. s (t) d'amplitudes respectives S1 et S2, les valeurs Dk(z) étant calculées, lorsque la forme de la tache 30 focale est un disque ou une sphère, par la formule (ulk,l(z)_u1k)(ulk,2(z)_u1k)(si - s2 ) (4) D k(Z) - d(S2uk"(z)ùS2uk ùSluk.2(z)+Sllfk)vl. 2906612 où . d est soit égal à 2 pour une imagerie à 2 dimensions, soit égal à 3 pour une imagerie à 3 dimensions, IVI est soit la surface de la tache focale pour une 5 imagerie à 2 dimensions, soit le volume de la tache focale pour une imagerie à 3 dimensions.

17. Procédé selon la revendication 15, dans lequel l'onde est une onde ultrasonore, L est égal à 2 et on utilise deux signaux S1 (t) = S1. s (t) et S2(t) = S2. s (t) 10 d'amplitudes respectives S1 et S2, s(t) étant un signal modulé en amplitude à une fréquence de modulation adaptée pour générer une force de radiation ultrasonore entraînant un déplacement local du milieu, les valeurs Dk(z) étant calculées, lorsque la forme de la tache focale est un 15 disque ou une sphère, par la formule : ' Drk(Z) = (u.I(z)_u)(u,k2(Z)_U)(s12 _S22) (4 ) d(S22uk'I (z) - S22uk - SI2uk'2 (z) + SI2uk )V . où . d est soit égal à 2 pour une imagerie à 2 dimensions, 20 soit égal à 3 pour une imagerie à 3 dimensions, IVI est soit la surface de la tache focale pour une imagerie à 2 dimensions, soit le volume de la tache focale pour une imagerie à 3 dimensions.

18. Procédé selon l'une quelconque des 25 revendications 10 à 17, dans lequel, au cours de l'étape de calcul, on part d'une conductivité supposée y et on répète les sous-étapes suivantes : a) on résout numériquement l'équation suivante : Jdiv(yyDuk ) = 0 en tout point z du milieu (9) -yDuk n = jk sur la surface extérieure 30 y étant une valeur précédemment estimée de la conductivité, 20 2906612 21 b) on calcule une erreur estimée ek sur la conductivité, c) on résout l'équation suivante : (11) Jdiv(yVvk )_ -div(ekouk) en tout point z du milieu, yVuk • n = 0 sur la surface extérieure. d) on met à jour la conductivité par la formule : 5 (12) Yk =-Y 2 A2Vuk •Vvk 2 A2Vuk +ek, i où vk est la solution de l'équation (11) et uk est la solution de l'équation (9), et on utilise Yk comme nouvelle valeur de conductivité y, avec un autre ensemble de courants jk, générant des lignes 10 de courant non parallèles à celles générées par l'ensemble de courants jk dans au moins une zone du milieu, les sous- étapes a) à d) étant réitérées jusqu'à satisfaction d'un critère d'arrêt.

19. Procédé selon la revendication 18, dans lequel, 15 au cours de la sous-étape b), on calcule une erreur estimée k e sur la conductivité par la formule : (10) ek = Ek / AVuk • Vuk -y .

20. Procédé selon la revendication 18, dans lequel, au cours de la sous-étape b), on calcule une 20 erreur estimée ek sur la conductivité par la formule : (10' ) ek =(Ek/AQuk Quk ù /Ek /AVuk •Vuk -y +1).

21. Procédé selon l'une quelconque des revendications 18 à 20, dans lequel, au cours de la sous- 2906612 22 étape (d) . pour chaque point z du milieu, on cherche quel indice k de conditions électriques correspond à la plus grande énergie A 2V uk ce qui donne une fonction k(z), 7 5 on met à jour la conductivité comme suit ( 2 A'Vuk(') •Vvk(') 2 A'Vuk(:) i et on utilise Yk(z) comme valeur de conductivité y.

22. Procédé selon l'une quelconque des revendications précédentes, dans lequelle milieu à imager 10 est un tissu biologique.

23. Procédé selon la revendication 22, dans lequel le milieu à imager est un organe humain.

24. Procédé selon l'une quelconque des revendications 1 à 21, dans lequel le milieu à imager est 15 le sous-sol terrestre.

25. Dispositif adapté pour la mise en oeuvre d'un procédé selon l'une quelconque des revendications précédentes. Yk(:) _ Y