Brevet FR3022668A1 — publié le 20151225

Système de tomographie dynamique dit « 4D » Le domaine général de l'invention est celui de la tomographie. Les domaines plus particuliers de l'invention concernent l'étude de l'évolution temporelle de matériaux soumis in situ à des sollicitations contrôlées. A titre d'exemples, on peut citer le cas des écoulements diphasiques immiscibles dans un milieu poreux, à différentes saturations, ou encore la déformation mécanique de matériaux soumis à des efforts croissants. La tomographie est une technique qui consiste à reconstruire le volume d'un objet à partir d'une série de radiographies bidimensionnelles prises sous différents angles. Généralement, les rayonnements utilisés sont des rayons X choisis pour leur aptitude à traverser la matière. Cette technique est basée sur la différence d'absorption du rayonnement par les éléments de volume de la pièce radiographiée. Elle est ainsi capable de révéler la présence d'inhomogénéités comme des porosités ou des grains de différentes matières. Le principal avantage de cette technique est de pouvoir observer à l'intérieur de l'échantillon de manière non destructive. Le synoptique d'un tomographe est représenté en figure 1. Il comprend principalement : - Un générateur GX de rayons X dont l'énergie et l'intensité sont adaptées à la taille et à l'absorption de l'échantillon E à analyser. Le faisceau de rayons X peut traverser une épaisseur comprise entre quelques millimètres et quelques centimètres, selon la nature de la matière traversée et l'énergie ou la fréquence des rayons X; - Un ensemble de platines motorisées PM en translation et en rotation portant l'échantillon E. Les courses sont généralement de quelques dizaines de centimètres. Cet ensemble comporte en particulier une platine en rotation permettant de faire tourner l'échantillon E d'un angle 0 sur 360 degrés avec une grande précision ; - Un détecteur matriciel D dont la définition en ligne et en colonne est généralement de plusieurs milliers de pixels. A titre d'exemple, le détecteur peut être constitué d'un scintillateur au césium assurant la conversion des rayons X en photons visibles, couplé à un ensemble de 3022668 2 photodétection dans cette bande de longueurs d'onde visibles. Cet ensemble peut être à base de détecteurs « CMOS », acronyme signifiant « Complementary Metal Oxide Semiconductor » ou « CCD », acronyme signifiant « Charge Coupled Deviced ». On note r - (r I, r2) un repère pris 5 dans le plan du détecteur matriciel et I(r , 0) l'intensité au point r de la radiographie lorsque l'échantillon E a subi un angle de rotation 0; - Des moyens de pilotage et de commande du générateur et des différentes platines non représentés sur la figure 1 ; - Des moyens d'acquisition et de traitement des différentes 10 données et de reconstruction tridimensionnelle de l'échantillon à partir des différentes radiographies. Ces moyens sont représentés par un calculateur C sur la figure 1. On utilise le cologarithme de l'intensité et non pas l'intensité elle-15 même pour linéariser l'absorption qui, par nature, est une fonction exponentielle décroissante de l'épaisseur traversée de l'échantillon E. On définit le sinogramme s(r , 0) par : s(r , 0) -log I(r, 0) où I(r , 0) représente l'intensité ayant traversé I 0(r) l'échantillon et reçue par le détecteur et Io(r) l'intensité reçue en l'absence 20 d'échantillon. L'image tridimensionnelle de l'échantillon analysé est représentée par la densité radiographique j(x) au point x. La « reconstruction » est l'opération mathématique linéaire ou transformée de Radon inverse, notée R, qui permet de calculer f(x) à partir de l'ensemble des radiographies : 25 j(x) = R[s(r, , 0)] Les points x sont également appelés « voxels », contraction de « VOlumetric piXELs », le terme « pixel » étant, lui-même, la contraction de « PICTure ELement ». Selon l'application, leur taille varie entre une fraction de micromètre et quelques centaines de micromètres. Compte-tenu de la 30 grande quantité de données à traiter pour effectuer cette reconstruction tomographique, il est aujourd'hui souvent fait appel à une implémentation de l'algorithme de reconstruction sur processeurs graphiques « GPU », acronyme signifiant « Graphics Processing Unit ». 3022668 3 Actuellement, lorsque l'on souhaite suivre l'évolution temporelle d'un milieu soumis par exemple à des sollicitations mécaniques extérieures, il convient d'acquérir une série d'images tomographiques, et de les analyser individuellement. 5 Les écoulements diphasiques immiscibles en milieux poreux constituent une classe importante de problèmes où l'imagerie tridimensionnelle résolue en temps est utilisée pour caractériser ces milieux, par exemple pour l'industrie pétrolière. Ainsi, pour caractériser à différents niveaux de saturation les perméabilités partielles et la pression capillaire, il 10 faut imager, pour différents volumes de fluide injectés, l'ensemble du milieu poreux et segmenter les deux phases fluides. La démarche est illustrée par la figure 2. Dans un premier temps, on réalise une tomographie de l'échantillon E que l'on souhaite analyser dans son état de référence. En appelant si(r,0) le sinogramme de l'échantillon imagé, on calcule une 15 première reconstruction tridimensionnelle/1(x), fi(x) = R[s 0)] Dans un second temps, on applique la sollicitation souhaitée sur l'échantillon. La sollicitation S est symbolisée par une flèche noire en chevrons sur la figure 2. On réalise une seconde tomographie de l'échantillon 20 E. En appelant s2(r,0) le nouveau sinogramme de l'échantillon, on calcule une seconde reconstruction tridimensionnelle f2(x) qui vaut : .f2(x) = R[s2(r, 0)] La connaissance de ces deux fonctions permet de caractériser la transformation T due à la sollicitation S sur l'échantillon E. 25 A titre d'exemple, les figures 3, 4 et 5 représentent schématiquement la même section d'un milieu poreux P à différents stades de l'invasion de l'espace poral de cette section par un fluide. Ainsi, la figure 3 montre une section du milieu poreux P dans son état initial à l'instant to. L'espace poreux est supposé saturé par un fluide a schématisé par un 30 niveau de blanc et le squelette solide du milieu poreux est en gris sombre. Lors de l'invasion du milieu poreux par un second fluide b représenté en noir sur les figures 4 et 5, celui-ci occupe progressivement la section comme on le voit sur ces deux figures représentant le milieu poreux aux instants ti et '2. Chacune de ces figures nécessite une grande quantité d'informations pour 3022668 4 être caractérisée. La figure 6, qui montre la différence entre les deux états t1 et to, nécessite une bien moindre quantité d'informations. Tous les phénomènes de croissance ou de propagation relèvent de cette même classe. Ainsi, la propagation d'une fissure dans un solide où 5 l'on image directement la variation de contraste liée à l'ouverture de la fissure peut être décrite par ce même formalisme. La solidification donnant lieu à une séparation de phases constitue un autre exemple de ces phénomènes. Pour analyser, à partir de tomographies, les champs de déplacements et de déformations d'un échantillon lorsqu'il est soumis à une 10 sollicitation mécanique S, on utilise généralement une technique dite de corrélation d'images volumiques, plus connue sous l'acronyme anglais « DVC », signifiant « Digital Volume Correlation ». Son principe est exposé en figure 7. Sous l'effet de la sollicitation, les points .v à l'intérieur de 15 l'échantillon se sont déplacés aux points X(x). L'hypothèse la plus courante de la DVC est de supposer que les coefficients d'absorption locaux n'ont subi qu'un mouvement sans autre altération, soit : f2(X) =fi(x) que l'on peut aussi écrire : 20 .f2M= .fe - U(X)) équation (1) où U(X) est la forme eulérienne du champ de déplacement. Le champ de déplacement U(X) peut être décomposé sur une base de champs de déplacements élémentaires g : U(X) = aiyfi(X) 25 Les champs g sont, par exemple, des fonctions de forme d'éléments finis. Les amplitudes a, sont déterminées au moyen de la minimisation de la fonctionnelle suivante : {a}= argmin tb} (f2 (X - (X - b1yi1(X)))2 dX équation (2) On obtient ainsi le champ des déplacements U(X), à partir duquel 30 il est possible de procéder à une exploitation quantitative des déformations, des éventuelles discontinuités, ou tout autre traitement interprétatif renseignant sur les propriétés mécaniques de l'échantillon. Lors d'un essai mécanique, chaque état de sollicitation demande l'acquisition de 3022668 5 radiographies et la reconstruction d'une image tridimensionnelle correspondante. Un des inconvénients de cette méthode est qu'elle nécessite l'enregistrement et l'exploitation d'une très grande quantité de données. Par 5 exemple, si l'on souhaite avoir une définition de l'image qui soit de l'ordre de 1000 pixels dans chaque direction de l'espace et en faisant l'hypothèse courante que l'échantillonnage angulaire est uniforme et de pas tel que le plus grand déplacement d'un point matériel de l'échantillon entre deux radiographies consécutives est inférieur à un pixel sur le détecteur, il est 10 nécessaire d'enregistrer de l'ordre de 3100 radiographies de l'échantillon, chaque image ayant une définition de 1000 x 1000 pixels, soit au total l'enregistrement de plus de trois milliards de pixels. Dans la plupart des applications, l'enregistrement demande donc une certaine durée, souvent comprise entre une et quelques heures. Pendant ce laps de temps, il est, 15 bien entendu, nécessaire que l'échantillon n'évolue pas. Or, pour certaines applications, il est impossible de garantir la stabilité de l'échantillon si la durée d'enregistrement est trop élevée. La puissance de la source de rayons X, ou la sensibilité du détecteur sont le plus généralement les facteurs limitants pour augmenter la cadence d'acquisition. 20 Par conséquent, la durée d'enregistrement limite le champ des applications de la DVC par tomographie. En particulier, il est impossible d'enregistrer des phénomènes dynamiques. De façon générale, pour suivre l'évolution de l'image d'un 25 échantillon de matériau soumis à une transformation T connue au cours du temps t, il est naturel de considérer, pour chaque pas de temps, l'acquisition d'une nouvelle image, soit f(x; t), dont le contenu d'information est donné par celui d'une image tridimensionnelle multipliée par le nombre d'acquisitions temporelles N1. L'image à l'instant t1 peut être décrite comme l'effet d'une 30 transformation T de l'image à l'instant to, soit encore : f(x; - Tex; to»pJ Equation (3) Dans cette écriture, la dénomination collective p rassemble l'ensemble des paramètres qui définissent quantitativement et de manière exhaustive la transformation T. Selon cette hypothèse, la connaissance de .f(x; et des paramètres p suffit à définir entièrement f(x; t1). Or, il existe de 3022668 6 très nombreux exemples où l'information nécessaire à caractériser les paramètres p est très sensiblement moindre que celle définissant f(x; to). Le principe du tomographe selon l'invention est de mesurer p plutôt que f(v; ro, et ce en exploitant un nombre de projections limité à deux pour chaque 5 instant t. Il est à noter que le fait de conduire l'analyse sur un grand nombre d'instants peut rendre l'analyse plus robuste et plus performante. Dans le cas où la lecture de l'image est « directe », comme pour l'exemple des écoulements diphasiques en milieux poreux, il est avantageux de ne s'intéresser qu'à la différence entre deux états successifs, par exemple 10 les instants t1 et to ou les instants t2 et t1. Ainsi, la figure 6 représente la variation de l'état du milieu poreux entre les instants t1 et t0. En effet, dans ce cas, la différence ne concerne qu'un nombre limité de voxels, ici d'une part restreint à l'espace poral, et d'autre part confiné dans une région de faible extension. Cette propriété est qualifiée de « parcimonie » et est d'autant plus 15 facilement satisfaite que la fréquence d'acquisition est élevée, limitant ainsi les différences entre deux images. Le fait de disposer de l'image initiale permet de s'abstraire de la complexité du squelette solide qui reste invariant. Les seules variations attendues sont des transitions entre le niveau de gris du fluide envahi a vers celui du fluide envahisseur b. Ainsi, la différence 20 f(x; -f(x; to) ne peut prendre que les deux valeurs 0 ou la différence des niveaux de gris des deux fluides pour une invasion irréversible, et ce exclusivement si à l'instant to, le voxel x était occupé par le fluide a. Il s'agit d'un problème de reconstruction binaire parcimonieuse. Dans ce cas, la transformation T se résume à l'addition du champ f(x; ti) -f(x; to) binaire et 25 parcimonieux et ce sont ces deux caractères qui réduisent la complexité du problème. De plus, le front d'invasion est limité à l'interface entre les deux fluides, ce qui réduit encore un peu plus la difficulté du problème. L'acquisition d'une série temporelle de paires de projections combinée à une rotation continue du milieu poreux permet de contraindre le problème de 30 manière supplémentaire. En effet, une fois une zone envahie, elle ne peut plus évoluer. Dans ce cas, il peut être intéressant de représenter l'ensemble des transformations T, par un champ supporté par l'espace poral et dont la valeur est le temps auquel le pore a été envahi par le second fluide. Un seul temps suffit si cette invasion est irréversible. De plus, la continuité du fluide 3022668 7 envahisseur impose que ce temps d'invasion ne présente pas de minima locaux. Comme nous l'avons vu précédemment, une autre classe de 5 transformation décrit l'évolution temporelle d'un objet ou d'une pièce mécanique par son seul champ de déplacement. Ce champ peut être un simple mouvement de corps rigide ou plus complexe et attaché à un ensemble de pièces mobiles, donnant lieu à un ou plusieurs mécanismes et qui peut être caractérisé par quelques degrés de liberté cinématiques 10 décrivant le mouvement d'objets rigides les uns par rapport aux autres. Plus généralement, les objets déformables peuvent être décrits au premier ordre par un champ de déplacement tridimensionnel dont la complexité est génériquement bien moindre que celle de la microstructure de l'objet lui-même. Dans ce cas, la transformation T est décrite par l'équation 15 (1), et les paramètres p sont ici les composantes du champ de déplacement. Comme on l'a vu, le champ de déplacement U(X) peut être décomposé sur une base de champs de déplacements élémentaires g : U(X) = et les amplitudes {a} constituent les paramètres p de la transformation. Les champs g sont par exemple des fonctions de forme 20 d'éléments finis. Là encore, la base de champs de déplacements élémentaires g comporte beaucoup moins d'inconnues que la base de données d'enregistrement. Ainsi, en reprenant l'exemple précédent, chaque élément fini peut avoir des dimensions 10 fois supérieures à celles d'un voxel. En prenant en compte que le champ de déplacement est calculé selon 25 les trois directions de l'espace, il comporte donc environ 330 fois moins de données que l'enregistrement primaire des voxels, soit 3 millions de données pour U(X) contre un milliard pour f(X). Le système de tomographie selon l'invention exploite cette 30 propriété. Comme on l'a vu, pour analyser les déplacements d'un échantillon sous l'effet d'une sollicitation, le nombre de données nécessaires est bien moins important que celui nécessaire à la tomographie initiale de l'échantillon. En pratique, dans un grand nombre d'applications, deux images de l'échantillon déformé suffisent. L'algorithme de détermination du champ 35 de déplacement nécessite quelques aménagements par rapport à la 3022668 8 technique DVC décrite ci-dessus. Il sera dénommé P-DVC pour « Projectionbased Digital Volume Correlation » dans la suite. Ainsi, le tomographe selon l'invention comporte les moyens nécessaires permettant de réaliser deux radiographies simultanées, 5 permettant de suivre l'évolution de la déformation de l'échantillon en temps réel. On appelle ce type de tomographie : « tomographie dynamique ou tomographie 4D ». Plus précisément, l'invention a pour objet un tomographe comportant au moins une première source d'émission, un premier détecteur 10 matriciel et des moyens de calcul agencés pour réaliser une tomographie initiale d'un objet à partir de radiographies issues du premier détecteur matriciel prises sous différents angles, caractérisé en ce que : le tomographe comporte une seconde source d'émission, un second détecteur matriciel et ; 15 lorsque l'objet est soumis à une sollicitation connue à un instant donné, les moyens de calcul déterminent les évolutions de l'objet soumis à ladite sollicitation à partir de la seule connaissance d'une première radiographie de l'objet sollicité issue du premier détecteur matriciel, d'une seconde radiographie de l'objet sollicité issue du second détecteur matriciel 20 et de la tomographie initiale, la première radiographie et la seconde radiographie étant prises au même instant donné. Avantageusement, le tomographe comporte des moyens agencés de façon à enregistrer une succession temporelle de couples de premières et de secondes radiographies de façon à reconstituer l'évolution de l'objet dans 25 le temps. Avantageusement, l'objet étant un milieu poreux, la sollicitation est l'écoulement d'un second fluide dans ledit milieu poreux comportant un premier fluide. Avantageusement, la sollicitation est une sollicitation mécanique, 30 l'objet se déformant sous cette sollicitation. Avantageusement, les moyens de calcul mettent en oeuvre un algorithme de type « P-DVC » ou « Projection-based Digital Volume Correlation », le champ de déplacement de l'objet sollicité étant déterminé à partir de la minimisation de fonctions dépendant de la première et seconde 3022668 9 radiographies de l'objet sollicité, et de la tomographie initiale effectuée dans l'état de référence. Avantageusement, la première et la seconde source d'émission sont des sources de rayonnement X. 5 L'invention sera mieux comprise et d'autres avantages apparaîtront à la lecture de la description qui va suivre donnée à titre non limitatif et grâce aux figures annexées parmi lesquelles : La figure 1 déjà commentée représente le synoptique d'un 10 tomographe selon l'art antérieur ; La figure 2 déjà commentée représente le principe général d'une transformation entre deux états analysée par l'utilisation d'images 3D de ces différents états ; Les figures 3, 4 et 5 représentent schématiquement trois états 15 successifs d'une section de milieu poreux lors d'un écoulement diphasique immiscible ; La figure 6 représente la différence de structure du milieu poreux entre deux instants, ceux des figures 3 et 4 de l'écoulement du fluide ; La figure 7 déjà commentée représente le principe général d'une 20 transformation entre deux états analysée par l'utilisation d'images 3D de ces différents états dans le cas particulier d'une déformation mécanique analysée par DVC ; La figure 8 représente le synoptique d'un tomographe selon l'invention ; 25 La figure 9 représente le principe de la DVC modifiée selon l'invention. A titre d'exemple non limitatif, un tomographe selon l'invention est représenté en figure 8. Il comprend principalement : 30 - Deux générateurs notés GX1 et GX2 de rayons X dont la puissance est adaptée à la taille et à l'absorption de l'échantillon E à analyser. Ces générateurs ont sensiblement les mêmes caractéristiques techniques ; - Deux détecteurs matriciels D1 et D2, chaque détecteur étant 35 associé à un des deux générateurs de rayons X. Les deux détecteurs sont 3022668 10 disposés dans deux plans différents. Généralement, les deux plans sont perpendiculaires entre eux. Ces deux détecteurs ont sensiblement les mêmes caractéristiques techniques. Leur définition en ligne et en colonne est généralement de plusieurs milliers de pixels. A titre d'exemple, la technologie 5 employée est à base de détecteurs « CMOS », acronyme signifiant « Complementary Metal Oxide Semiconductor » avec scintillateur en césium. - Un ensemble de platines motorisées PM en translation et en rotation portant l'échantillon E. Cet ensemble comporte en particulier une platine en rotation permettant de faire tourner l'échantillon E d'un angle 0 sur 10 360 degrés avec une grande précision. La platine est disposée de façon à ne pas faire écran entre les générateurs et les détecteurs. L'échantillon est situé à des distances sensiblement identiques des deux détecteurs et des deux générateurs ; - Des moyens de pilotage et de commande des deux 15 générateurs et des différentes platines non représentés sur la figure 8; - Des moyens d'acquisition et de traitement des différentes données et de reconstruction tridimensionnelle de l'échantillon à partir des différentes radiographies. Ces moyens sont représentés par un calculateur C sur la figure 8. 20 Comme on l'a dit, le tomographe selon l'invention a un grand nombre d'applications : - Il peut être utilisé pour caractériser les écoulements diphasiques dans les milieux poreux. La technique utilisée est 25 une variante de la reconstruction adaptée aux champs binaires et parcimonieux ; - Il peut être utilisé pour réaliser des mesures de déplacement d'un échantillon d'un matériau soumis à une sollicitation connue. La technique utilisée est une variante de la technique 30 dite « Projection-based Digital Volume Correlation » précédemment exposée. Dans le premier cas, un exemple d'algorithme de reconstruction binaire est présenté dans l'article intitulé « Efficient Binary Tomographic 3022 66 8 11 Reconstruction » de S. Roux, H. Leclerc et F. Hild, J. Math Imaging Vis. (2014) 49:335-351. Dans le second cas, la technique s'inspire des principes décrits dans l'article intitulé « Projection savings in CT-based Digital Volume 5 Correlation » de H. Leclerc, S. Roux et F. Hild publié dans « Experimental Mechanics »(D01 10.1007/s11340-014-9871-5). La différence fondamentale entre le dispositif décrit dans cet article et le tomographe selon l'invention est que le dispositif décrit comporte un et un seul générateur de rayon X et un et un seul détecteur. Ainsi, pour réaliser deux radiographies sous deux angles 10 différents, il est nécessaire de faire tourner l'échantillon avec les contraintes que cela entraîne, la principale étant l'impossibilité de suivre l'évolution de l'échantillon en temps réel sans perturbations. Comme précédemment, la mesure comporte deux étapes. Son 15 principe est exposé en figure 9. Dans une première étape, on réalise classiquement une première tomographie de l'échantillon E que l'on souhaite analyser. En appelant s1(r,0) les images tomographiques de l'échantillon, on calcule une première fonction de reconstruction tridimensionnelle fi qui vaut, comme précédemment : 20 fi(x) = R[si(r,0)] Cette tomographie est réalisée par les deux ensembles fonctionnant simultanément. Les deux séries d'images obtenues sont mises en correspondance, permettant d'associer les deux référentiels des détecteurs. 25 Dans une seconde étape, on applique la sollicitation S souhaitée sur l'échantillon E. La sollicitation S est symbolisée par une flèche noire en chevrons sur la figure 9. On réalise alors uniquement deux radiographies simultanées de l'échantillon E sollicité, la première prise par le premier ensemble comportant le premier générateur GX1 et le premier détecteur D1, 30 la seconde prise par le second ensemble comportant le second générateur GX2 et le second détecteur D2. Par construction, ces deux radiographies sont prises sous deux angles 0 différents sans qu'il soit nécessaire de faire tourner l'échantillon. On note s2(r, G) ces deux images. On sait, par ce qui a été dit précédemment, que les données 35 issues de deux radiographies sont potentiellement suffisantes pour calculer 3022668 12 les paramètres p de la transformation T. Ainsi, dans le cas particulier de la déformation d'un matériau on calcule les amplitudes {a} du champ de déplacement : U(X)- a ,v ,(X) 5 sur la base de champs g(x) selon l'algorithme P-DVC ou équivalent. Classiquement, les champs g sont des fonctions de forme d'éléments finis. Les amplitudes a, sont déterminés au moyen de la minimisation de la fonction suivante, voisine de l'équation (2), en utilisant les mêmes notations que précédemment : 10 argmin lb} (s2 (r, 0) - P0 [f,(X - b,yJ,(X))1)2 dr où Po correspond à la projection de la reconstruction tridimensionnelle ji corrigée du déplacement sur les plans des détecteurs. Dans le cadre de l'invention, les calculs sont faits à partir des informations issues uniquement des deux images s2(r 15 Selon une réalisation possible du calcul des champs de déplacement selon l'algorithme P-DVC, l'estimation se fait par itérations successives n. A chaque itération, on construit une image de référence artificiellement déformée:11(X). On a la relation : :11(X) = - (X)) 20 Entre deux itérations, on a la relation simple : n+1 7 I a, = a, + aa, Les termes de correction da,n' sont obtenus à partir de la linéarisation des différences dans les images projetées. On définit : g/.1 (r' 0)= P0 v:fin (x) 25 Le calcul des termes de correction da:1+1 est obtenu par la résolution du système linéaire suivant : M nda = dB," avec M,n, =Iff g , (Ag (r , 0)dr 30 et dB, if(s 2(r , 0) - Po[j"(X)]).g; (r , 0)dr 3022668 13 Le processus d'itérations est arrêté dès que les termes de correction sont inférieurs à un seuil prédéterminé. En pratique, la convergence est obtenue après une dizaine d'itérations. 5 Ainsi, on diminue de façon très significative le nombre de radiographies nécessaires à la détermination du champ de déplacement. De plus, en enregistrant une succession temporelle de couples de premières et de secondes radiographies, il est possible de reconstituer l'évolution du champ de déplacement de l'objet dans le temps d'où la terminologie de 10 « tomographie dynamique » ou « 4D ». L'avantage principal du tomographe selon l'invention est de permettre un suivi dynamique de l'évolution d'un échantillon soumis à une sollicitation donnée. On peut, à titre de premier exemple, l'utiliser pour suivre 15 l'évolution dynamique du front entre deux fluides immiscibles dans un écoulement dynamique sous l'effet d'une pression appliquée. A titre de second exemple, on peut également analyser la dépendance temporelle de la cinématique d'un échantillon soumis à une sollicitation mécanique. Un tomographe selon l'invention à deux ensembles 20 sources/détecteurs a nécessairement un coût supérieur à un tomographe comportant un seul ensemble. Cependant, l'augmentation de coût est assez limitée dans la mesure où de nombreux éléments du système tomographique peuvent être mutualisés comme les écrans de protection contre les radiations, les platines motorisées, les dispositifs antivibratoires ou encore 25 les systèmes de climatisation. REVENDICATIONS1. Tomographe comportant au moins une première source d'émission (GX1), un premier détecteur matriciel (D1) et des moyens de calcul (C) agencés pour réaliser une tomographie initiale d'un objet (E) à partir de radiographies issues du premier détecteur matriciel prises sous différents angles, caractérisé en ce que : le tomographe comporte une seconde source d'émission (GX2), un second détecteur matriciel (D2) et, lorsque l'objet est soumis à une sollicitation dont le cadre descriptif est connu à un instant donné, les moyens de calcul déterminent les évolutions de l'objet soumis à ladite sollicitation à partir de la seule connaissance d'une première radiographie de l'objet sollicité issue du premier détecteur matriciel, d'une seconde radiographie de l'objet sollicité issue du second détecteur matriciel et de la tomographie initiale, la première radiographie et la seconde radiographie étant prises au même instant donné. 2. Tomographe selon la revendication 1, caractérisé en ce que, le tomographe comporte des moyens agencés de façon à enregistrer une succession temporelle de couples de premières et de secondes radiographies de façon à reconstituer l'évolution de l'objet dans le temps. 3. Tomographe selon l'une des revendications précédentes, caractérisé en ce que l'objet étant un milieu poreux, la sollicitation est un écoulement diphasique d'un second fluide dans ledit milieu poreux comportant un premier fluide. 4. Tomographe selon l'une des revendications 1 ou 2, caractérisé en ce que la sollicitation est une sollicitation mécanique, l'objet se déformant sous la sollicitation. 30 5. Tomographe selon la revendication 4, caractérisé en ce que les moyens de calcul mettent en oeuvre un algorithme de type « P-DVC » ou « Projection-based Digital Volume Correlation », le champ de déplacement 25 3022668 15 de l'objet sollicité étant déterminé à partir de la minimisation de fonctions dépendant de la première et seconde radiographies de l'objet sollicité, et de la tomographie initiale effectuée dans l'état de référence. 5 6. Tomographe selon l'une des revendications précédentes, caractérisé en ce que la première et la seconde source d'émission sont des sources de rayonnement X.