Brevet FR2740561A1 — publié le 19970430

La présente invention concerne une méthode pour évaluer la variation d'intensité d'un rayonnement polychromatique ayant un spectre de fréquence connu, après traversée d'une certaine épaisseur d'un corps constitué d'une ou plusieurs phases de composition connue, dans le but de corriger des erreurs liées à des artefacts. La méthode selon l'invention permet de tenir compte par exemple de différents facteurs: composition chimique du corps, la distribution en énergie du faisceau polychromatique, éléments constituants ledit corps, variation du coefficient d'absorption du rayonnement etc. La méthode selon l'invention convient particulièrement pour modéliser et/ou corriger le durcissement d'un faisceau provenant d'un rayonnement polychromatique traversant un corps absorbant. Dans toute la description, on désigne sous l'expression "durcissement d'un faisceau", l'augmentation de l'énergie moyenne d'un faisceau pendant la traversée d'un corps absorbant". La présente invention s'applique particulièrement bien pour simuler le durcissement d'un faisceau de rayons X et corriger les erreurs introduites par l'apparition d'un artefact associé à ce durcissement. La méthode selon l'invention, trouve des applications notamment pour mesurer l'épaisseur "effective" de matière traversée par un faisceau de rayon X polychromé dans un corps absorbant de composition chimique connue. Elle est appliquée, par exemple, pour déterminer la répartition et le pourcentage de phases dans un composé dont on connaît la nature des phases le constituant. La méthode selon l'invention trouve des applications notamment dans le cadre de la caractérisation de roches (roches réservoirs) dans le domaine pétrolier. Lorsqu'un faisceau d'un rayonnement polychromatique traverse un corps absorbant, il subit un changement dû notamment à son interaction avec le ou les éléments constituants ce corps. Ainsi le spectre en énergie du faisceau après la traversée de ce corps diffère de celui qu'il possède à l'entrée. Ce changement est fonction notamment des énergies contenues dans le spectre et peut varier selon la position de la zone traversée par le faisceau et son épaisseur. Ainsi, les plus faibles valeurs d'énergie étant associées à un fort coefficient d'absorption produisent dans la majorité des cas les artefacts de durcissement les plus importants. Dans des applications à caractère quantitatif, il s'avère intéressant de posséder une méthode permettant de s'affranchir des incertitudes introduites par la nature polychromatique du rayonnement. I1 existe un certain nombre de méthodes décrites dans l'art antérieur pour minimiser, voire éliminer les incertitudes introduites par la polychromaticité du rayonnement. Une des méthodes consiste à utiliser un filtre, par exemple une plaque de matériau absorbant une partie au moins du rayonnement et plus particulièrement les plus faibles énergies. Cette méthode présente néanmoins l'inconvénient d'absorber une partie du faisceau et de conduire à une perte d'intensité. Selon une autre méthode, on réalise un étalonnage en utilisant un échantillon modèle ou de référence. Néanmoins un tel échantillon doit être homogène en composition et présenter une atténuation proche de celle de l'objet à analyser et il n'est pas toujours possible de respecter ces deux derniers critères de sélection. Ces méthodes ne permettent donc pas de s'affranchir de manière simple des phénomènes introduits par la polychromaticité du faisceau. La méthode selon l'invention a pour objet d'évaluer la variation d'intensité d'un rayonnement polychromatique ayant un spectre de fréquence connu, après traversée d'une certaine épaisseur d'un corps ou objet constitué d'une ou plusieurs phases de composition connue, dans le but de corriger des erreurs liées à des artefacts. Elle est caractérisée en ce que l'on évalue cette variation d'intensité par une modélisation de l'influence d'un certain nombre de paramètres définis liés à la nature et à des propriétés physico-chimiques de chacune des phases, influant sur le coefficient d'absorption du corps dans une ou plusieurs bandes spectrales du rayonnement polychromatique. Le modèle permettant l'évaluation de la variation d'intensité, est défini par exemple par la fonction de modélisation I(E) suivante: 1 H intensité du rayonnement x Ej H énergie du canal i du spectre X Lk H longueur traversée dans la phase k Pk H densité de la phase k j.k H fraction de 1' élément j dans la phase coefficient d' atténuation massique pour chaque élément j de la phase k Pour réaliser le modèle, on modélise par exemple les influences respectives de paramètres caractéristiques d'une pluralité de couches successivement traversées par le rayonnement. Suivant un mode de réalisation, on évalue la variation d'intensité sélectivement pour une ou plusieurs bandes de fréquence du rayonnement. I1 est possible de sélectionner l'énergie du rayonnement incident, de manière par exemple à éliminer les trop faibles énergies, qui introduisent des incertitudes sur les mesures. Suivant un mode de mise en oeuvre, on compare la variation d'intensité mesurée après traversée d'une certaine zone d'un corps à analyser, avec des valeurs de variations déduites du modèle pour déterminer au moins un paramètre caractéristique du corps analysé. Par cette comparaison, on peut déterminer l'épaisseur effective de la zone traversée par un rayonnement connaissant la composition du corps analysé, ou bien, si l'épaisseur du corps analysé est connue, déterminer la concentration de la zone traversée par le rayonnement en un ou plusieurs de ses constituants. Lorsque l'on cherche à déterminer la texture d'un corps ou objet, on effectue plusieurs mesures d'intensité du rayonnement ayant traversé au moins une zone du corps analysé en faisant varier l'angle d'incidence du rayonnement par rapport à cette zone, on associe à chacune desdites valeurs mesurées une valeur d'intensité déterminée à partir d'un rayonnement polychromatique puis une valeur d'intensité déterminée pour un rayonnement monochromatique, et à partir desdites valeurs on déduit la composition de ladite zone traversée par le rayonnement. On utilise avantageusement un rayonnement polychromatique, suivant plusieurs directions et on détermine la distribution de constituants du corps dans au moins une zone traversée par ledit rayonnement. La méthode peut être utilisée par exemple pour simuler le durcissement d'un rayonnement polychromatique ayant traversé un corps, ou bien encore pour caractériser un échantillon, tel qu'un échantillon de roche ou un composite, dans le domaine pétrolier. Avantageusement, la méthode selon l'invention trouve des applications dans tout le domaine des analyses basées sur l'exploitation de radiographies X. Elle permet notamment d'augmenter le pouvoir de discrimination des phases actuellement obtenu par les techniques classiques. La méthode selon l'invention trouve ainsi avantageusement son application pour étudier les milieux poreux, par exemple dans des domaines dont la liste non exhaustive est donnée ci-après à titre indicatif: - en minéralogie pour identifier des minéraux au sein d'une roche et étudier leur distribution éventuelle, - en mécanique des fluides en milieux poreux, par exemple pour définir la géométrie des grains ou de pores et leur connexité, ces paramètres ayant une influence non négligeable sur le comportement du fluide, - en mécanique des roches, par exemple pour vérifier l'état de micro-fissuration des échantillons de roche. Dans des domaines autres que celui de l'industrie pétrolière, la méthode selon l'invention est avantageusement utilisée pour contrôler par tomographie la texture des catalyseurs, de matériaux composites ou roches réservoirs par exemple. D'autres caractéristiques et avantages de la méthode selon l'invention, apparaîtront à la lecture de la description donnée ci-après, à titre illustratif et nullement limitatif, en se référant aux dessins annexés où: les figures 1A et 1B représentent des tables de correspondance entre des valeurs de transmittance observées pour un faisceau polychromatique (TRP) et des transmittances monochromatiques (TRM) pour deux valeurs d'intensité; les figures.2, 3 montrent schématiquement un exemple de réalisation pour un montage permettant de valider la méthode selon l'invention, la figure.4 montre schématiquement les spectres Sf respectifs du rayonnement avant traversée d'un échantillon (I1) et après sa traversée (I2); la figure.5 montre schématiquement les spectres Sf respectifs du rayonnement polychromatique de rayons X, après traversée de l'échantillon, le premier (I2) étant celui mesuré et le deuxième (I3) étant déduit du modèle; la figure.6 montre schématiquement les spectres Sf respectifs du rayonnement polychromatique de rayons X, le premier (Ix) étant celui de la source, le deuxième (Ixf) étant celui filtré à travers une plaque de graphite; la figure.7 montre schématiquement les spectres Sf respectifs du rayonnement polychromatique de rayons X, le premier (Ie) étant celui obtenu expérimentalement après passage à travers un filtre et l'échantillon analysé, le deuxième (IM) étant le spectre correspondant obtenu par le modèle; et la figure 8 représente schématiquement un exemple d'application de la méthode pour réaliser la tomographie d'un échantillon Le spectre en énergie d'un faisceau issu d'un rayonnement polychromatique traversant un corps absorbant est modifié du fait notamment de son interaction avec les éléments chimiques constituant le corps. Cette modification dépend aussi des valeurs d'énergie du faisceau, car le coefficient d'absorption ou coefficient d'atténuation du corps varie en fonction de l'énergie. Le corps est composé de plusieurs espèces chimiques et aussi de plusieurs phases ayant des natures différentes. Pour un faisceau monochromatique traversant un corps constitué d'éléments chimiques différents et de plusieurs phases, il est possible de déterminer l'intensité du faisceau après la traversée du corps en faisant appel à la loi de Beer-Lambert qui s'exprime par exemple de la manière suivante: I = lo exp(-t) où I est l'intensité du faisceau monochromatique ayant traversé le corps, Io l'intensité du faisceau avant traversée du corps et A est le coefficient d'absorption lié à l'énergie du faisceau monochromatique, déterminé en tenant compte des sections efficaces d'absorption et des densités atomiques des éléments. Néanmoins, cette formule ne peut être appliquée de façon valable pour des rayonnements polychromatiques du fait de la distribution en énergie du spectre de tels faisceaux. La mise en oeuvre de la méthode selon l'invention est donnée ci-après à titre illustratif et nullement limitatif en relation avec un rayonnement polychromatique, par exemple un faisceau de rayons X traversant un corps constitué d'éléments chimiques différents et éventuellement ayant des phases présentant des natures différentes. On sait de plus que les artefacts introduits ou incertitudes introduites dans la mesure sont fonction de l'épaisseur de l'échantillon. On procède par exemple aux étapes suivantes: En supposant connue la composition chimique et minéralogique du corps traversé par le faisceau de rayons X, on mémorise les paramètres indicatifs de cette composition dans des mémoires, celles par exemple d'un micro-contrôleur, sous forme de tables de données, connaissant les éléments chimiques constitutifs des phases, la densité de ces phases et le coefficient massique d'atténuation pour chaque élément chimique j d'une phase k, on a accès au coefficient d'atténuation du faisceau dans chacune de ces phases, à partir d'une hypothèse sur la répartition des phases k ou d'une connaissance préalable de leur répartition, il est possible d'évaluer la distance (l'épaisseur) effectivement traversée par le faisceau, c'est-à-dire la distance tenant compte de l'angle d'incidence du faisceau par rapport à un axe du corps Lk, par ailleurs, on a mémorisé le spectre d'énergie du faisceau indicent, avant traversée du corps, par exemple sous forme d'un tableau de valeurs intensité/énergie. A partir de ces données, on détermine la valeur d'intensité I du faisceau ayant traversé au moins une phase k sur une longueur Lk ceci pour une valeur d'énergie Ei donnée, par exemple à l'aide de la fonctionsuivante: l H intensité du rayonnement x E H énergie du canal i du spectre X Lk e longueur traversée dans la phase k Pk H densité de la phase k fi k coefficient d' atténuation massique pour chaque élément j de la phase k Dans l'exemple donné ci-dessus de mise en oeuvre de la méthode selon l'invention, on modélise le corps en le subdivisant en fonction des phases de nature différente qu'il comporte. I1 en résulte une série de valeurs de transmittance polychromatique pour chaque couche qui est fonction notamment de l'épaisseur effective traversée par le faisceau. Parallèlement, la transmittance pour plusieurs faisceaux monochromatiques d'énergies choisie est calculée pour chaque épaisseur, ou couche C1. Le pas de subdivision est choisi notamment par rapport à la précision que l'on souhaite obtenir et/ou en fonction des dimensions caractéristiques des phases. L'intensité initiale du faisceau prise en compte dans l'expression peut être une intensité obtenue pour la couche précédente à partir de la fonction de simulation 1. Par exemple sur la figure 2B, on a schématisé un corps absorbant 1 modélisé par un empilement de couches ayant des épaisseurs Cl, par exemple correspondant chacune à une phase de nature différente. La méthode selon l'invention permet par exemple de simuler le durcissement du faisceau d'un rayonnement polychromatique. A l'aide de la fonction de simulation (1), on détermine la transmittance cumulée par intégration à chaque couche. On peut ainsi dresser une table de correspondance entre la transmittance polychromatique et la transmittance monochromatique. Les figures 1A et 1B montrent graphiquement deux exemples de table de conversion ainsi obtenues respectivement pour des faisceaux monochromatiques d'énergie de 14 keV et de 10 keV pour un échantillon de quartz et une source X sur une cible en or. L'utilisation de la méthode pour de la simulation a été validée en utilisant un microscope électronique à balayage de type JEOL JSM-35CF équipé d'une diode Si(Li), associé à un système d'analyse Tracor TN 5500. Le principe de la méthode selon l'invention consiste à créer un rayonnement X au moyen d'un faisceau d'électrons frappant une cible métallique, que l'on applique à un corps ou objet à analyser et à utiliser un détecteur tel qu'un spectromètre pour mesurer le spectre de la source X (utilisé comme donnée d'entrée par la modélisation), puis dans un deuxième temps le spectre des rayons X après traversée du corps absorbant que l'on cherche à comparer à un spectre obtenu en partant du spectre de la source X auquel on a appliqué la fonction de simulation (1) pour toutes les valeurs d'intensité. Sur le schéma de principe de la figure 2, la source de rayonnement X polychromatique qui est par exemple un rayonnement de type synchrotron, ou encore produite par bombardement d'un faisceau électronique 1 sur une cible 2, génère des rayons X en direction de l'échantillon ou matériau absorbant 3. Le rayonnement issu de l'échantillon, est reçu par un détecteur 4 tel qu'un spectromètre positionné de façon à recueillir le spectre du rayonnement X ayant traversé le corps absorbant. Le corps 3 situé entre la source et le détecteur 4, présente une surface d'entrée Se située du côté du rayonnement X et une surface de sortie Ss située du côté du détecteur 4 (Fig.3). Le faisceau de rayons X est envoyé sous un angle d'incidence a par rapport à un plan du corps, celui de sa face d'entrée par exemple. Sur la figure 3, le rayonnement X traverse dans une direction perpendiculaire à la surface d'entrée Se du corps, ressort modifié par la surface de sortie Ss du corps avant d'être détecté par le spectromètre 4. Par ailleurs, on acquiert le spectre d'intensité du rayonnement X avant sa traversée de l'échantillon et on le mémorise dans les mémoires du micro-contrôleur 11. Les figures 4 et 5 montrent les résultats de simulation obtenus lorsque le corps absorbant est un échantillon de polyméthacrylate de méthyle (PPMA) de 1 mm d'épaisseur. La tension d'accélération des électrons est, dans cet exemple, de 25 kV et la cible choisie est une feuille de cuivre de 1 Wm d'épaisseur. Le cuivre fournit une raie principale à 8 keV, énergie qui présente une atténuation d'environ 47%. La figure 4 montre néanmoins que l'atténuation des énergies inférieures à 5 keV est quasi totale, et que l'atténuation du spectre est d'environ 36 %. Sur la figure 4, sont représentés les spectres mesurés de la source de rayons X avant par la courbe I1 et après la traversée, courbe I2 d'un échantillon de PPMA ayant une épaisseur de 1 mm. On applique la fonction de modélisation (1) à chacune des valeurs d'énergie de la courbe I1 et on détermine les valeurs d'intensités corrigées représentées par la courbe I3 sur la figure 5. On remarque que les courbes I2 et I3 représentées sur la figure 5, sont sensiblement identiques et ainsi que les valeurs déterminées à l'aide de la fonction de modélisation (1) sont sensiblement identiques aux valeurs d'intensité mesurées par le spectromètre 4, après la traversée du corps, ce qui valide la méthode de simulation. Selon un autre mode de mise en oeuvre de la méthode, il est possible d'utiliser un filtre permettant d'atténuer les faibles énergies du spectre, ou filtre passe-haut, ces dernières étant responsables des artefacts les plus importants, donc du durcissement le plus important. On procède selon une manière sensiblement identique à la méthode décrite ci-dessus, le corps modélisé étant constitué du corps absorbant et du filtre. Avantageusement, une des mises en oeuvre de la méthode permet de déterminer l'épaisseur d'un échantillon dont on connaît la composition chimique et minéralogique. On procède par exemple de la manière suivante: on établit le modèle de référence en déterminant des séries de valeurs d'intensité du rayonnement monochromatique à l'aide de la fonction de modélisation (1) pour différentes valeurs d'énergie, et pour différentes couches Ci, et l'on mémorise ces données sont mémorisées, par exemple dans le micro-contrôleur, on effectue ensuite une série de mesures d'intensité du faisceau polychromatique de rayons X après la traversée de l'échantillon, et on mémorise ces données dans le micro contrôleur, ceci pour plusieurs zones de l'échantillon, on compare ensuite chacune des valeurs d'intensité mesurées aux valeurs déterminées pour le modèle et on lui associe, par exemple la valeur d'intensité la plus proche fournie par lui, et sachant qu'à une valeur d'intensité est associée une épaisseur ou couche CP, on en déduit la valeur de l'épaisseur effective de la zone traversée par le faisceau. Par épaisseur effective, on entend l'épaisseur réelle de solides traversés en tenant compte de l'angle d'incidence du rayonnement X. En connaissant l'angle d'incidence du faisceau, il est possible de remonter à la valeur de l'épaisseur traversée. La mise en oeuvre de la méthode selon l'invention permet aussi de déterminer la concentration d'un élément et/ou d'une phase contenue dans un corps dont on connaît la nature chimique et/ou minéralogique. Par la méthode, on peut déterminer par exemple, la concentration de platine dans un catalyseur ayant comme matrice l'alumine. On fait l'hypothèse que l'épaisseur est connue, et on calcul les valeurs de l'intensité transmise correspondant à des valeurs de concentration du platine et/ou de l'alumine correspondant à une valeur d'intensité transmise. Ces valeurs d'intensité déterminées sont ensuite mémorisées par exemple dans un tableau sous forme d'un couple (intensité du faisceau de rayons X, composition d'une phase choisie). On effectue ensuite une mesure de l'intensité du rayonnement X ayant traversé le corps à analyser, à savoir la matrice et l'élément, et on compare cette valeur mesurée aux valeurs déterminées et mémorisées dans le tableau de données. D'une façon identique à la méthode décrite ci-dessus, on associe à la valeur d'intensité mesurée la valeur la plus proche d'intensité déterminée et on en déduit la valeur de concentration pour un élément donné. Cette manière de procéder peut aussi être étendue sans sortir du cadre de l'invention pour déterminer le pourcentage d'une phase dans une matrice connue, par exemple lorsque l'on cherche à déterminer la quantité de feldspath dans une roche réservoir. La méthode de détermination de la variation l'intensité d'un faisceau de rayons X ayant traversé un échantillon peut être avantageusement utilisée dans des applications de tomographie, et en particulier pour réaliser des microtomographies d'un échantillon. La plupart des logiciels connus de reconstruction tomographique sont adaptés à des faisceaux monochromatiques. Il s'avère donc nécessaire de s'affranchir de la polychromaticité du faisceau de rayonnement, par exemple en effectuant une correction des données mesurées à l'aide de la méthode précédemment décrite. La figure 8 schématise un microtomographe qui présente notamment comme avantage d'obtenir des précisions au niveau de la résolution dans la gamme du micromètre. Le microtomographe schématisé à la figure 8 est constitué de deux parties, une colonne 5 et un tiroir 6, par exemple. La colonne est par exemple une colonne d'un microscope électronique classique et ne sera donc pas décrite en détail. Le tiroir 6 comporte les cibles 2 bombardées par le faisceau électronique, qui sont la source du rayonnement X, présentant des caractéristiques sensiblement identiques à celles données à la figure 4 (spectre I1). A l'intérieur de ce tiroir 6, se trouvent un porte-échantillon 7 solidaire d'un axe de rotation 8, sur lequel on a positionné l'échantillon 9 L'axe de rotation 8 permet le basculement de l'échantillon selon trois directions repérées sur la figure par trois axes orthogonaux. La rotation de l'échantillon peut se faire manuellement par des multiples de 3 ou plus, par exemple ce qui correspond à un maximum de 120 rotations par tour complet. Une autre manière de procéder consiste à automatiser la rotation à l'aide d'un dispositif approprié, cette variante permettant d'augmenter le nombre de pas et de passer à 1200 pas par tour complet. Le dispositif de détection 10 du faisceau ayant traversé l'échantillon comprend une caméra CCD pouvant être associée à tout dispositif permettant de s'affranchir d'effets parasites et d'augmenter la qualité et la précision des mesures. Le dispositif de détection 10 est associé à un système de traitement 11 programmé pour mettre en oeuvre les méthodes classiques de tomographie et par ailleurs la méthode d'évaluation de variation d'intensité selon l'invention L'énergie du faisceau d'électrons est choisie, par exemple de façon à maximiser la production des photons de la raie choisie, la limite supérieure étant donnée par la capacité de la cuve haute tension de la colonne électronique dans cet exemple de réalisation. Le courant du faisceau est réglé pour obtenir un faisceau intense tout en conservant une sonde électronique de faibles dimensions. Le choix du pas de rotation de l'échantillon est, quant à lui, fixé pour obtenir une résolution spatiale de l'ordre de celle des projections, le nombre de projections M obéissant à la loi suivante M est supérieur ou égal à IJ2 *N où N est le nombre de pixels d'une ligne des projections. La méthode consiste à acquérir ensuite un ensemble de valeurs d'intensité associées à des angles de rotation de l'échantillon, pour un faisceau de rayons X. On obtient ainsi un tableau de données regroupant des valeurs d'intensité mesurées Mi et des valeurs d'angle Bi. Simultanément, on a déterminé les valeurs d'intensité IPi correspondant au spectre incident et à différentes épaisseurs d'échantillons. On a de plus déterminé les valeurs d'intensité IMi obtenues pour plusieurs faisceaux monochromatiques, dans une gamme de valeurs d'énergie comprenant au moins les valeurs d'énergie du spectre du rayonnement envoyé sur l'échantillon et à différentes épaisseurs d'échantillons, de façon à constituer le modèle de référence. Ces deux ensembles de valeurs sont regroupés, par exemple en un seul tableau de données et mémorisés dans le microcontrôleur. Suivant la méthode, on compare chacune des valeurs d'intensité mesurées Mi aux valeurs d'intensité déterminées pour le faisceau polychromatique IPi, et après comparaison, on remplace la valeur mesurée Mi par la valeur déterminée la plus proche. On compare ensuite la valeur IPi aux valeurs obtenues pour des faisceaux monochromatiques et on conserve la valeur d'intensité que l'on aurait obtenue à l'aide d'un faisceau monochromatique, la plus proche par exemple. C'est cette valeur qui va être utilisée dans le logiciel de tomographie pour déterminer et caractériser la zone d'échantillon traversée. Les méthodes éventuelles de traitement permettant d'améliorer les résolutions des images habituellement utilisées en tomographies peuvent être avantageusement utilisées sans sortir du cadre de l'invention. Cette manière de procéder permet avantageusement d'obtenir des précisions spatiales de l'ordre, ou inférieure, à la dizaine de microns en comparaison aux résolutions offertes en mode radiographique qui est inférieure à 3 Rm. On a décrit l'application de la méthode à l'analyse tomographique d'objets ou d'échantillons. On ne sortirait pas du cadre de l'invention toutefois en appliquant la méthode selon l'invention à la corrections d'artefacts dans des appareils de tomographies du corps humain par exemple, tels que des scanners X. REVENDICATIONS 1) Méthode pour évaluer la variation d'intensité d'un rayonnement polychromatique ayant un spectre de fréquence connu, après traversée d'une certaine épaisseur d'un corps constitué d'une ou plusieurs phases de composition connue, dans le but de corriger des erreurs liées à des artefacts, caractérisée en ce que l'on évalue cette variation d'intensité par une modélisation de l'influence d'un certain nombre de paramètres définis liés à la nature et à des propriétés physico-chimiques de chacune des phases, influant sur le coefficient d'absorption du corps dans une ou plusieurs bandes spectrales du rayonnement polychromatique. 2) Méthode selon la revendication 1, caractérisée en ce que le modèle permettant l'évaluation de la variation d'intensité, est défini par la relation suivante: l H intensité du rayonnement x Ej énergie du canal i du spectre X Lk H longueur traversée dans la phase k Pk H densité de la phase k fi k coefficient d' atténuation massique pour chaque élément j de la phase k 3) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on modélise les influences respectives de paramètres caractéristiques d'une pluralité de couches successivement traversées par le rayonnement. 4) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on évalue la variation d'intensité sélectivement pour une ou plusieurs bandes de fréquence du rayonnement. 5) Méthode selon l'une des revendications précédentes, caractérisée en ce que l'on compare la variation d'intensité mesurée après traversée d'une certaine zone d'un corps à analyser, avec des valeurs de variations déduites du modèle pour déterminer au moins un paramètre caractéristique du corps analysé. 6) Méthode selon la revendication précédente, caractérisée en ce que, la composition du corps analysé étant connue, on détermine par comparaison l'épaisseur effective de la zone traversée par le rayonnement. 7) Méthode selon la revendication 5, caractérisée en ce que, l'épaisseur du corps analysé étant connue, on détermine par comparaison la concentration de la zone traversée par le rayonnement en un ou plusieurs de ses constituants. 8 - Méthode selon la revendication 5, caractérisée en ce que l'on effectue plusieurs mesures d'intensité du rayonnement ayant traversé au moins la dite zone du corps analysé en faisant varier l'angle d'incidence du rayonnement par rapport à cette zone, on associe à chacune desdites valeurs mesurées une valeur d'intensité déterminée à partir d'un rayonnement polychromatique puis une valeur d'intensité déterminée pour un rayonnement monochromatique, et à partir desdites valeurs on déduit la composition de ladite zone traversée par le rayonnement. 9 - Méthode selon la revendication 8, caractérisée en ce que l'on irradie le corps étudié au moyen d'un faisceau de rayonnement X polychromatique, suivant plusieurs directions et on détermine la distribution de constituants du corps dans au moins une zone traversée par ledit rayonnement. 10 - Utilisation de la méthode selon l'une des revendications précédentes, pour simuler le durcissement d'un rayonnement polychromatique ayant traversé un corps. 1 1 - Application de la méthode selon l'une des revendications précédentes pour caractériser un échantillon, tel qu'un échantillon de roche ou un composite, dans le domaine pétrolier.