FR 3050030 A1 20171013 FR 1653142 A 20160408 1 DOMAINE TECHNIQUE DE L'INVENTION La présente invention porte de manière générale sur le magnétisme appliqué aux géosciences. Elle concerne la quantification de l'aimantation d'un échantillon sédimentaire et plus particulièrement la quantification de la variation des aimantations rémanentes le long d'une portion de carotte sédimentaire. L'invention trouve par exemple pour application particulièrement avantageuse les domaines de la paléo-climatologie, des paléo-environnements et de la stratigraphie sédimentaire au sens large, incluant la détection de niveaux riches en particules issues de polluants métalliques ÉTAT DE LA TECHNIQUE Le magnétisme environnemental est de plus en plus utilisé dans des études scientifiques très variées par exemple dans les domaines de la paléo-climatologie, des paléo-environnements et de la stratigraphie sédimentaire au sens large. Il est souvent très utile d'estimer les aimantations d'échantillons qui se présentent sous la forme d'une carotte sédimentaire. Une carotte sédimentaire se présente souvent sous forme d'une carotte s'étendant selon une direction principale, désignée par la suite axe long, et qui présente une section constante le long de cet axe. La carotte sédimentaire comprend au moins une portion sédimentaire obtenue par carottage. Bien souvent, elle comprend également une gouttière, par exemple en plastique, de section en U. Les côtés du U peuvent présenter une longueur identique ou non. Lorsque c'est le cas, chaque côté du U mesure par exemple 2 cm (10-2 mètre). La longueur de la carotte sédimentaire, c'est-à-dire sa dimension selon son axe long mesure environ 1,5 m. La gouttière est remplie par le sédiment. Elle est présente également un couvercle. La gouttière et le couvercle forment alors une enveloppe de protection, de préférence fermée, pour encapsuler le sédiment. Une telle carotte sédimentaire est communément qualifiée pas le vocable anglais « U-channel », signifiant canal en U. La variation de l'aimantation rémanente le long d'une carotte sédimentaire dépend des conditions de dépôts sédimentaire et donc des conditions climatiques, tectoniques et de la composition des grains se déposant. 3050030 2 Il est souvent très utile d'estimer les aimantations rémanentes d'une manière continue le long de carottes sédimentaires pour caractériser les porteurs de l'aimantation dans les couches sédimentaires successives qu'ils contiennent. L'aimantation rémanente est par exemple générée en passant une carotte 5 sédimentaire dans un anneau de Halbach. Habituellement, les techniques mises en oeuvre pour estimer les aimantations rémanentes utilisent des magnétomètres supraconducteurs. Ces techniques présentent plusieurs inconvénients. Les magnétomètres existants combinent une acquisition d'aimantations 10 rémanentes parallèles à l'axe long de la carotte sédimentaire et un système de mesures avec des sondes placées loin de la carotte sédimentaire (sondes supraconductrices). L'acquisition de l'aimantation rémanente parallèlement à l'axe long de la portion de carotte sédimentaire donc perpendiculairement aux laminations 15 sédimentaires implique un lissage du signal. De plus, la distance importante des sondes ne permet pas de caractériser l'aimantation rémanente avec une résolution spatiale satisfaisante. Un objet de la présente invention est de proposer une solution pour fournir une acquisition et une mesure d'aimantations rémanentes dont la résolution spatiale 20 est améliorée. Un autre objet de la présente invention est d'atteindre ce but tout en réduisant le coût ou sans augmenter de manière significative le coût des solutions existantes. Les autres objets, caractéristiques et avantages de la présente invention apparaîtront à l'examen de la description suivante et des dessins d'accompagnement. 25 Il est entendu que d'autres avantages peuvent être incorporés. RÉSUMÉ DE L'INVENTION Selon un mode de réalisation, l'invention porte sur un banc de mesure du champ magnétique d'une carotte sédimentaire s'étendant selon une direction 30 principale, désignée axe long, le banc de mesure comprenant au moins: -un dispositif de déplacement en translation de la carotte sédimentaire au sein du banc de mesure, -une chambre de mesure comprenant une sonde à vanne de flux agencée pour être disposée en regard et de préférence au contact d'une portion de carotte sédimentaire de manière mesurer le champ magnétique de ladite portion de carotte sédimentaire située en regard de la sonde à vanne de flux. 3050030 3 Le dispositif de déplacement est configuré de manière à faire translater la carotte sédimentaire par rapport à la sonde à vanne de flux de manière à ce que la sonde à vanne de flux mesure successivement le champ magnétique d'une succession de portions de la carotte sédimentaire. 5 Un autre aspect de la présente invention concerne un procédé d'estimation des aimantations rémanentes de portions d'une carotte sédimentaire s'étendant selon une direction principale, le procédé comprenant au moins les étapes suivantes: i. Définition d'une succession de N volumes élémentaires (correspondant aux 10 laminations sédimentaires) et de leurs dimensions, chaque volume élémentaire correspondant à une portion de carotte sédimentaire, les volumes élémentaires étant adjacents et successivement disposés selon ladite direction principale de la carotte sédimentaire; ii. Translation de la carotte sédimentaire relativement à une sonde à vanne de 15 flux de manière à faire passer chaque volume élémentaire en regard de la sonde à vanne de flux ; iii. Mesure par la sonde à vanne de flux du champ magnétique produit par chaque volume élémentaire ; iv. En utilisant les mesures de champ magnétique de la sonde à vanne de flux : 20 détermination d'une aimantation rémanente M_x propre à chaque volume élémentaire, cette étape comprenant; - une première identification « directe » d'une pluralité de valeurs possibles d'aimantation rémanente M_x, cette première identification comprenant: 25 o Le calcul de valeurs de champ magnétique B_x à partir d'une pluralité de valeurs d'aimantations rémanentes M_x choisies arbitrairement ; o l'identification, parmi ladite pluralité de valeurs d'aimantations rémanentes M_x choisies arbitrairement, d'une ou plusieurs 30 valeurs d'aimantation rémanente M_x pour lesquelles le champ magnétique B_x calculé est compris entre 0.7 et 1.3 fois, de préférence entre 0.8 et 1.2 fois, et de préférence entre 0.9 et 1.1 fois, le champ magnétique B_x mesuré par la sonde à vanne de flux ; 3050030 4 - une deuxième identification « inverse » d'une ou plusieurs valeurs d'aimantation rémanente M_x, cette deuxième identification comprenant: o le calcul d'une valeur M_x d'aimantation rémanente en fonction 5 de la mesure du champ magnétique mesuré par la sonde à vanne de flux. Ainsi, la valeur M_x d'aimantation rémanente est obtenue par inversion des mesures du champ magnétique mesuré par la sonde à vanne de flux. De préférence, les étapes de translation et de mesure sont effectuées en 10 utilisant le banc de mesure selon l'invention. De manière particulièrement avantageuse, l'utilisation d'une sonde à vanne de flux en combinaison avec les autres caractéristiques de la présente invention permet d'améliorer la résolution spatiale des carottes sédimentaires. Dans le cadre du développement de la présente invention il a été remarqué 15 que les bancs de mesure existants ne permettent pas d'effectuer des mesures fiables et précises sur des carottes sédimentaire possédant une forte aimantation, qu'elle soit ponctuelle ou affectant toute la carotte. En effet, ces aimantations intenses saturent les magnétomètres supraconducteurs. En effet, lors d'études ayant mené à la présente invention, les aimantations 20 rémanentes, n'ont pu être acquises que de manière partielle, les sections de la carotte sédimentaire les plus concentrées en minéraux magnétiques saturant le magnétomètre supraconducteur. Ainsi, le banc de mesure automatisé selon l'invention comprend une sonde à vanne de flux (« fluxgate »). Ce type de sonde est habituellement utilisé pour la mesure 25 de champs magnétiques ambiants mais n'est pas utilisé pour la mesure d'aimantations rémanentes. Lors d'essais pratiqués dans le cadre du développement de la présente invention, la sonde vanne de flux a montré une sensibilité excellente de l'ordre de 1 nT. Cette valeur de champ est générée par des « U-channels » présentant des moments 30 magnétiques de l'ordre de 10-8 Am2, des valeurs de détection largement suffisantes pour mesurer des aimantations rémanentes isothermes. Cette sensibilité permet donc de mesurer des aimantations rémanentes isothermes avec une grande précision pour la plupart des types de sédiments. Un résultat inattendu de l'étude préliminaire correspond au fait que la résolution spatiale de la sonde vanne de flux est meilleure 35 d'un facteur deux par rapport aux magnétomètres supraconducteurs. 3050030 5 Par ailleurs cette solution induit un coût significativement moins élevé que celui des solutions existantes, en particulier des solutions mettant en oeuvre des magnétomètres supraconducteurs. Par ailleurs les capteurs supraconducteurs sont fragiles et peu transportables. 5 Le banc de mesure selon l'invention est quant à lui robuste et bien plus facilement transportable, notamment du fait de l'absence de tout équipement cryogénique. Par ailleurs, on peut séparer la chambre de mesure comprenant la sonde à vanne de flux des profilés sur lesquels la carotte sédimentaires se déplace pour passer en regard de la sonde à vanne de flux. Cette démontrabilité facilite grandement 10 la transportabilité du banc de mesure. Par ailleurs la sonde à vanne de flux de sature pas ou sature beaucoup moins lorsque les aimantations sont fortes. Dans la présente invention la sonde à vanne de flux couplée aux étapes du 15 procédé permet d'estimer avec une résolution spatiale améliorée les aimantations rémanentes des différents volumes élémentaires formant chacun une portion de la carotte sédimentaire. Selon un autre aspect, la présente invention porte sur un produit programme 20 d'ordinateur ou sur un média non-transitoire lisible par un ordinateur, comprenant des instructions, qui lorsqu'elles sont effectuées par au moins un processeur, exécute les étapes du procédé selon l'invention mentionné ci-dessus, en particulier l'étape iii mentionnée ci-dessus. 25 De manière avantageuse, le procédé selon l'invention peut présenter au moins l'une des caractéristiques et étapes facultatives suivantes prises seules ou en combinaison. ladite étape de définition est effectuée de manière à ce que chaque volume élémentaire de la carotte sédimentaire présente une même épaisseur, 30 l'épaisseur étant mesurée selon ladite direction principale. alternativement, ladite étape de définition est effectuée de manière à ce que chaque volume élémentaire de carotte sédimentaire présente une épaisseur indépendante de celle des autres volumes élémentaires, 'épaisseur étant mesurée selon ladite direction principale. 35 De préférence ladite étape de définition est effectuée de manière à ce que chaque volume élémentaire de carotte sédimentaire présente au moins un 3050030 6 critère d'homogénéité également désigné critère sédimentologique significativement homogène, ledit critère étant pris parmi : une couleur, une texture. Selon un mode de réalisation, le procédé comprend une étape préalable de 5 fourniture d'une carotte sédimentaire présentant des aimantations rémanentes naturelles, puis une étape de transformation des aimantations rémanentes en aimantations rémanentes artificielles, ladite étape de transformation comprenant un passage de la carotte sédimentaire dans un anneau de Halbach. 10 Un autre aspect de la présente invention concerne un appareil d'estimation des aimantations rémanentes de portions d'un échantillon sédimentaire formant une carotte sédimentaire de section rectangulaire (« u-channel ») et s'étendant selon une direction principale, l'appareil comprenant un banc de mesure selon l'invention et un module de 15 calcul agencé pour mettre en oeuvre les étapes de procédé mentionnée ci-dessus et en particulier l'étape iv. BRÈVE DESCRIPTION DES FIGURES Les buts, objets, ainsi que les caractéristiques et avantages de l'invention 20 ressortiront mieux de la description détaillée d'un mode de réalisation de cette dernière qui est illustré par les dessins d'accompagnement suivants dans lesquels : La FIGURE 1 est une vue schématique en perspective d'un exemple de banc de mesure selon l'invention. La FIGURE 2 est une vue agrandie de la zone A de la figure 1. 25 La FIGURE 3 est une vue du dessus du banc de mesure illustré en figure 1. La FIGURE 4 est une vue de côté du banc de mesure illustré en figure 1. La FIGURE 5 est une vue selon la coupe B-B du banc de mesure illustré en figure 1 La FIGURE 6 est une vue selon la coupe C-C du banc de mesure illustré en 30 figure 1 La FIGURE 7 illustre des étapes d'un exemple de méthode d'estimation des aimantations rémanentes selon l'invention. La FIGURE 8 est une représentation schématique d'un échantillon de sédiments formant une carotte sédimentaire. 3050030 7 La FIGURE 9 est une représentation schématique d'un échantillon de sédiments formant une carotte sédimentaire après passage de la carotte sédimentaire dans un anneau de Halbach. Les FIGURES 10 à 13 illustrent des étapes d'une méthode d'inversion utilisée 5 dans le cadre de la présente invention. Les dessins sont donnés à titre d'exemples et ne sont pas limitatifs de l'invention. Ils constituent des représentations schématiques de principe destinées à faciliter la compréhension de l'invention et ne sont pas nécessairement à l'échelle des applications pratiques. En particulier les épaisseurs relatives des différentes pièces ne 10 sont pas forcément représentatives de la réalité. DESCRIPTION DÉTAILLÉE DE L'INVENTION Avant d'entamer une revue détaillée des modes de réalisation de l'invention, 15 sont énoncées ci-après des caractéristiques optionnelles qui peuvent éventuellement être utilisées en association ou alternativement : la chambre de mesure du banc de mesure comprend un bouclier magnétique formant un cylindre ouvert à chacune de ses extrémités à l'intérieur duquel sont situés la sonde à vanne de flux et ladite portion de carotte sédimentaire, de 20 préférence située en contact ou à proximité immédiate de la sonde à vanne de flux. le bouclier magnétique comprend en deux couches de mu-métal. Chaque couche forme un manchon. Les deux manchons sont concentriques. le dispositif de déplacement comprend une courroie de transport, entrainée en 25 déplacement par un moteur, et dans lequel la courroie de transport traverse la chambre de mesure. le moteur est un moteur pas à pas. Les aimantations rémanentes sont des aimantations rémanentes artificielles, de préférence isothermes. 30 La carotte sédimentaire présente une section carrée de 2 cm de côté et une longueur égale ou inférieure à 1,5 m. la sonde à vanne de flux est agencée pour être disposée au contact de la carotte sédimentaire. la sonde à vanne de flux est agencée pour être disposée perpendiculairement à 35 l'axe long de la carotte sédimentaire de manière à mesurer le champ magnétique produit localement par la portion de carotte sédimentaire 3050030 8 l'échantillon. La sonde à vanne de flux est agencée pour mesurer l'aimantation perpendiculairement à l'axe long de la carotte. La sonde à vanne de flux est configurée pour mesurer, perpendiculairement à ladite direction principale, le champ magnétique généré par ladite portion de 5 carotte sédimentaire. la sonde à vanne de flux est agencée de manière à mesurer perpendiculairement à ladite direction principale le champ magnétique généré par ladite portion de carotte sédimentaire située en regard de la sonde à vanne de flux. 10 la prise de mesure par la sonde à vanne de flux est synchronisée avec le déplacement en translation de la carotte sédimentaire. Ainsi la sonde à vanne de flux ne prend une mesure que lorsque qu'un déplacement de la carotte sédimentaire est achevé pour amener un volume élémentaire au droit de la sonde à vanne de flux. 15 le banc de mesure permet une prise de mesures ponctuelles. Alternativement, il permet une prise de mesure en mode scanner. Un même banc de mesure permet de prendre ces deux types de mesure. la sonde à vanne de flux est située verticalement en dessous de la carotte sédimentaire lorsque la carotte sédimentaire est située en partie au moins dans 20 la chambre de mesure. le dispositif de déplacement est configure de manière à faire translater la carotte sédimentaire par rapport à la sonde à vanne de flux de manière à ce que la sonde à vanne de flux mesure le champ magnétique .produit par les couches sédimentaires successives 25 avant la mesure de champ magnétique, la portion de carotte sédimentaire est aimantée perpendiculairement à son axe long à l'aide d'un anneau de Halbach. Dans la présente description les termes échantillon, échantillon sédimentaire ou carotte sédimentaire sont interchangeables et désignent, comme indiqué dans la 30 section relative à l'art antérieur, une portion de sédiment s'étendant selon une direction principale et qui peut être partiellement ou entièrement encapsulée dans une enveloppe. Un exemple de banc de mesure selon l'invention va maintenant être décrit en 35 détail en référence aux figures 1 à 6, puis un procédé pour estimer les aimantations rémanentes sera ensuite présenté en référence aux figures 7 à 13. 3050030 9 Le banc de mesure 1 comprend globalement une chambre amagnétique 30 à l'intérieur de laquelle est amenée une carotte sédimentaire 100 pour en mesurer les aimantations rémanentes. La chambre amagnétique 30 est portée par un cadre 10 qui comporte également un dispositif de déplacement de la carotte sédimentaire 100 afin 5 de déplacer cette dernière depuis l'extérieur jusqu'à l'intérieur de la chambre amagnétique 30. Ces différents éléments du banc de mesure 1 vont maintenant être détaillés. Le cadre 10 comprend des éléments verticaux 13 prenant appui sur le plan de 10 travail. Ces éléments sont dans cet exemple constitués de profilés, typiquement en aluminium. Ces profilés verticaux 13 supportent une poutre supérieure 11, s'étendant de préférence horizontalement et destinée à supporter la carotte sédimentaire 100 et sa courroie de transport 21. De préférence, au moins une autre poutre 12, de préférence horizontale et inférieure, est solidaire des profilés verticaux pour assurer la 15 robustesse et la stabilité du banc de mesure 1. Le dispositif de déplacement de la carotte sédimentaire 100 comprend la courroie de transport 21 qui s'étend sur et le long de la poutre 11 ainsi qu'un chariot 25 fixé à la courroie de transport 21 est sur lequel la carotte sédimentaire 100 est destinée reposer par gravité. De préférence le chariot 25 est inélastique et peu déformable. Le 20 chariot 25 s'étend sur toute la longueur la carotte sédimentaire 100 et dépasse de cette dernière comme illustré en figure 1 et 3. Cela permet de fixer les extrémités du chariot 25 à la courroie de transport 21. Selon un autre mode de réalisation la carotte sédimentaire 100 repose directement sur la courroie 21. 25 Cette courroie de transport 21 est fermée ou ouverte lorsqu'elle est fixée aux extrémités du chariot 25. Elle est entraînée en déplacement par un moteur 22 de préférence pas à pas, solidaire du cadre 10. De préférence, des poulies 23, 24 sont disposées sur le cadre 10 pour guider la courroie de transport 21 en faisant un renvoi d'angle depuis la poutre supérieure 11 30 vers la poutre inférieure 12. Ainsi, la courroie de transport 21, lorsqu'elle arrive à l'extrémité de la poutre supérieure 11, repart entre les poutres supérieure 11 et inférieure 12, en sens inverse, jusqu'au moteur 22. La courroie de transport 21 permet ainsi une translation de la carotte sédimentaire 100. Avantageusement, le banc de mesure comprend des guides 26 configurés 35 pour guider la translation de la carotte sédimentaire 100 en évitant son déplacement transversalement à la direction prévue de translation. Deux guides 26 sont ainsi 3050030 10 disposés de part et d'autre de la courroie 21 et du chariot 25 comme illustré en figure 1 et 3. De préférence, au milieu de la poutre supérieure 11, le cadre 10 supporte la chambre amagnétique 30. La chambre amagnétique 30 comprend une entrée et une 5 sortie. Elle est traversée de part en part par la carotte sédimentaire 100 lorsque cette dernière translate. La chambre amagnétique 30 comprend un bouclier magnétique 31 configuré pour empêcher aux perturbations magnétiques externes d'interférer avec l'intérieur de la chambre amagnétique 30. Ce bouclier 31 comprend typiquement une couche interne 10 31 et une couche externe 32 formant chacune un manchon cylindrique. De préférence, chacune des couches 31, 32 est en mumétal. Ces manchons cylindriques apparaissent en coupe sur les figures 5 et 6. Au niveau de la chambre amagnétique 30, la poutre supérieure 11 est discontinue. La chambre amagnétique 30 forme un logement 14 configuré pour 15 contenir un capteur de champ magnétique également désigné magnétomètre. Ce magnétomètre est une sonde à vanne de flux 40, habituellement désignée en anglais « fluxgate ». Sur les figures 5 et 6, pour des raisons de clarté, seul un boitier de la sonde à vanne de flux est illustré. La sonde à vanne de flux 40 est située verticalement en dessous de la carotte 20 sédimentaire 100 lorsque la carotte sédimentaire est située en partie au moins dans la chambre de mesure 30. Ce type de sonde est utilisé pour mesurer les champs magnétiques ambiants, par exemple le champ magnétique ambiant d'une pièce. En revanche, comme indiqué ci-dessus, dans les solutions connues ce type de sonde n'est jamais utilisé pour 25 mesurer le champ magnétique d'une carotte sédimentaire « U-channel » et en déduire les aimantations rémanentes. Comme cela apparaîtra dans la description qui suit, l'utilisation de cette sonde à vanne de flux confère de nombreux avantages. Avantageusement, la prise de mesure par la sonde à vanne de flux est synchronisée avec le déplacement en translation de la carotte sédimentaire. Ainsi en 30 fin de déplacement de la carotte sédimentaire pour amener un volume élémentaire au droit de la sonde à vanne de flux, cette dernière effectue une mesure. Le banc de mesure permet sélectivement une prise de mesure ponctuelle ou en mode scanner. 3050030 11 En référence à la figure 7, des étapes principales non forcément indispensables d'un exemple de procédé selon l'invention vont maintenant être détaillées. Une première étape (étape 710) consiste à se procurer un échantillon sous 5 forme d'une carotte sédimentaire par exemple de section continue. Cette carotte sédimentaire est typiquement obtenue par carottage dans une couche sédimentaire, puis encapsulation dans une enveloppe protectrice. Une telle carotte sédimentaire est illustrée en figure 8. Sur cette figure 8 on a également représenté un repère (x, y, z) dans lequel : 10 - x, axe positif vers le bas et perpendiculaire (transversal) à la direction principale selon laquelle s'étend la carotte sédimentaire 100; - y, axe positif vers l'utilisateur de l'invention, perpendiculaire (transversal) à la direction principale selon laquelle s'étend la carotte sédimentaire 100 ; - z, direction principale selon laquelle s'étend la carotte sédimentaire 100, 15 positif dans le sens contraire au déplacement de la carotte sédimentaire 100 lors de la mesure. Le sens de déplacement 102 de la carotte sédimentaire 100 est illustré en figure 8. 20 A l'étape 730 illustrée en figure 7, la carotte sédimentaire 100 est ensuite traité de manière à lui faire perdre son aimantation naturelle et à la transformer en une aimantation rémanente. Pour cela, on passe la carotte sédimentaire dans un anneau de Halbach, dont les propriétés bien connues de l'homme du métier du magnétisme environnemental, sont par exemple décrites dans la publication suivante Rochette, P., 25 Vadeboin, F., Clochard, L., 2001. Rock magnetic applications of Halbach cylinders, Physics of the Earth and Planetary Interiors 126, 109-117. Le passage de l'échantillon dans l'anneau de Halbach créé l'aimantation rémanente, c'est donc un passage préliminaire à la mesure du champ généré par cette 30 aimantation artificielle. Lors du passage de la carotte sédimentaire dans l'anneau de Halbach, ce dernier entraine une saturation de l'aimantation des sédiments dans une direction prédéterminée (ici -x). Les autres aimantations sont effacées par la ré-aimantation dans l'anneau de Halbach, sauf s'il existe un signal de haute coercivité. Bien souvent l'aimantation rémanente naturelle est effacée lors de traitements arrivant 35 bien avant le passage dans l'anneau de Halbach. 3050030 12 La figure 9 illustre la distribution et les intensités des aimantations dans la carotte sédimentaire après passage dans l'anneau de Halbach. Chaque volume élémentaire 101 est considéré comme ayant une aimantation uniforme. Chaque vecteur d'aimantation est au centre d'un volume élémentaire 101. La 5 hauteur de chaque flèche illustre l'intensité du champ magnétique mesuré pour le volume élémentaire qui lui est associé. Ensuite, la carotte sédimentaire 100, lorsqu'elle est présente dans la chambre amagnétique 30, est disposée au contact de la vanne de flux. La distance entre la sonde à vanne de flux et le sédiment correspond ainsi à 10 l'épaisseur de l'enveloppe qui encapsule le sédiment. Par exemple, cette enveloppe en plastique présente une épaisseur de 2 mm (10-3 mètres). Cette distance est notée « d » en figure 6. Si la carotte sédimentaire 100, n'est pas encapsulée dans une enveloppe alors c'est le sédiment qui est directement au contact de la sonde à vanne de flux. On 15 a alors d = O. L'aimantation rémanente est acquise perpendiculairement à la face de mesure grâce à l'anneau de Halbach. Ainsi, la face principale 41 est disposée au droit d'une portion 101 de la carotte sédimentaire 100. La section de la carotte sédimentaire 100 étant 20 rectangulaire, chaque portion forme un volume élémentaire rectangulaire. Les volumes élémentaires ou portions sont référencées en figure 9. Dans l'exemple illustré les volumes élémentaires forment chacun un prisme. Pour une position relative de la carotte sédimentaire 100 par rapport à la sonde à vanne de flux 40, cette dernière mesure un champ magnétique : celui 25 correspondant principalement à l'effet de l'aimantation du volume élémentaire 101 situé au droit de la face principale 41 de la sonde 40. Mais l'aimantation des volumes élémentaires adjacents va aussi influencer la valeur de la mesure. Typiquement, les dimensions de la face principale 41 sont les suivantes : - Largeur de 2 cm 30 - Longueur de 1,5 m La largeur et la longueur sont prises selon des dimensions respectivement perpendiculaire et parallèle à la direction principale selon laquelle s'étend la carotte sédimentaire 100 (l'axe long de la carotte sédimentaire). Le bouclier magnétique 31 présente un diamètre interne suffisamment grand 35 pour loger à la fois la sonde à vanne de flux 40 et la portion 101 de la carotte sédimentaire 100 située au droit de la face principale 41 de la sonde 40. 3050030 13 La carotte sédimentaire 100 est ensuite disposé sur le banc de mesure et orienté de manière à ce que la sonde à vanne de flux 40 mesure le champ magnétique de la portion 101 de la carotte sédimentaire 100 disposée au droit de sa face principale 41. La carotte sédimentaire 100 traverse la chambre amagnétique 31. La sonde 40 5 mesure alors le champ magnétique selon un pas de mesures prédéfini, chaque mesure étant donc influencée par la portion 101 de la carotte sédimentaire 100 (au droit du capteur), mais aussi par les portions adjacentes (étape 740). Les mesures de champ magnétique effectuées par la sonde à vanne de flux 40 sont ensuite transmises au module de calcul. Ce dernier peut être intégré dans le 10 banc de mesure, par exemple dans un module de commande comprenant une carte de commande 50 du banc de mesure 1 ou alternativement être déporté de ce dernier en étant en communication filaire ou sans fil avec ce dernier. Les mesures du champ magnétique effectuées par la sonde à vanne de flux 40 constituent ainsi un premier type de données entrantes pour que le module de 15 calcul détermine une estimation des aimantations rémanentes. Un deuxième type de données entrantes est constitué par des données issues d'une analyse sédimentologique de la carotte sédimentaire formant « U-channel ». Ce deuxième type de données porte sur l'épaisseur de volumes unitaires successifs 20 formant la carotte sédimentaire. L'épaisseur des volumes unitaires est prise selon la direction principale selon laquelle la carotte sédimentaire s'étend, c'est-à-dire selon l'axe z sur le repère (x, y, z) de la figure 8. Dans la suite de la description les termes volumes élémentaires et prismes seront utilisés de manière interchangeable. 25 Dans le procédé qui sera détaillé par la suite on considère que chaque volume élémentaire présente une aimantation homogène unique, chaque volume élémentaire étant sensé correspondre à une séquence sédimentologique unique Ces étapes de définition de l'épaisseur des volumes élémentaires sont illustrées en 750 et 760 sur la figure 7. 30 La définition de ces volumes unitaires ou volumes élémentaires est par exemple basée sur une analyse de paramètres géophysiques, géochimiques, minéralogiques et sédimentologiques si disponibles. Ainsi, le sédimentologue, en analysant visuellement l'échantillon, détermine une succession de portions 101 adjacentes, également désignées primes 101, qui ensemble forment la carotte 35 sédimentaire 100. L'épaisseur de chaque volume élémentaire 101 peut être définie indépendamment des épaisseurs des autres volumes élémentaires 101. 3050030 14 Alternativement, tous les volumes élémentaires101 présentent la même épaisseur. Dans la mesure du possible, chaque section présente une certaine homogénéité pour le ou les paramètres utilisés pour délimiter l'épaisseur de chaque volume élémentaire. Ainsi, chaque volume élémentaire 101 présente une certaine homogénéité de couleur 5 ou de texture, et présente une couleur ou une texture différente de celle de la section adjacente. Cette définition des volumes élémentaires basée sur une analyse de la carotte permet d'optimiser le pas de mesure. Elle permet ainsi d'optimiser le temps de mesure. 10 Le module de calcul fournit ensuite une estimation de l'aimantation rémanente de chaque volume élémentaire 101 de la carotte sédimentaire 100. Cette estimation est basée sur une méthode qui va maintenant être décrite en détail en référence aux figures 10 à 13. 15 La sonde à vanne de flux 40 uniaxe mesure le champ magnétique selon le sens (-x), soit positivement vers le haut. On utilisera le terme B_x pour désigner le champ magnétique mesuré par la sonde à vanne de flux selon ce sens. Comme indiqué ci-dessus en référence à l'étape 740 de la figure 7 et à la 20 figure 9, lors du passage de la carotte sédimentaire 100 dans l'anneau de Halbach ce dernier entraine une saturation de l'aimantation des sédiments dans la direction (x), mais positivement vers le haut, soit M_, Les autres aimantations sont effacées par la ré-aimantation dans l'anneau de Halbach, sauf s'il existe un signal de haute coercivité. 25 La figure 10 illustre le rendu des mesures de champ magnétique (B_x) le long de la carotte sédimentaire 100. Comme indiqué ci-dessus, chaque séquence sédimentaire sera assimilée à un volume élémentaire i rectangulaire (une section de la carotte sédimentaire), dans lequel on considère une aimantation homogène unique (M_x)i. Dans cet exemple les 30 volumes élémentaires forment chacun un prisme rectangulaire. Par défaut, dans l'approche ci-dessous on considère N volumes élémentaires de même épaisseur. Néanmoins et comme indiqué précédemment, l'invention s'étend aux modes de réalisation permettant de contraindre ces épaisseurs par une analyse visuelle préalable comme cela a été indiqué en référence aux étapes 750, 760 de la figure 7 et comme 35 cela sera détaillé par la suite. L'épaisseur standard est égale à l'intervalle entre deux 3050030 15 mesures, avec le centre de chaque volume élémentaire placé sous chaque mesure, donc N correspond au nombre de mesures. La méthode selon l'invention prévoit ensuite deux étapes principales : une 5 étape de modélisation directe et une étape d'inversion pour retrouver une série de moments dipolaires associés à chaque volume élémentaire rectangulaire. Etape 1 - modélisation directe La première étape consiste à modéliser directement les variations 10 d'aimantation de la carotte sédimentaire dans chaque volume élémentaire en utilisant une méthode basique 'essai et erreur' (en anglais trial and error). Parmi cifférentes possibilités de séries de valeurs (0,exy ouxy, (ux)3, (m_es ) ), l'algorithme va trouver laquelle prédit le mieux les observations B. L'équation reliant une mesure lax et un volume élémentaire rectangulaire aimanté (ici selon (-x)) est la suivante : 15 Equation (1) 1..), s.incttin , I? avec i, j, k, des indices égaux à 1 ou 2 et délimitant les arêtes du volume élémentaire dans les directions ( -y, -z, x), respectivement. Le coefficient s correspond 20 au produit (s,s,sk) avec s1=-1 et 52=+1, alors que (a, b, c) sont les coordonnées des coins du volume élémentaire par rapport au point de mesure, toujours selon ( -y, -z, x), et Ruk=(a2+b2+c2)0.5. Cette équation est tirée de la publication suivante : Plouff, D., 1976. Gravity and magnetic fields of polygonal prisms and application to magnetic terrain corrections, 25 Geophysics, 41(4), 727-741. Le repère (-y,-z,x) provient du fait que, dans cette publication, les équations sont définies dans un repère (N, E, vertical) alors que l'invention prévoit de préférence un repère spécifique à l'utilisation d'un échantillon U-channel, ce repère spécifique 30 étant défini ci-dessus en référence à la figure 8. Chaque mesure lax peut être influencée par plusieurs volumes élémentaires adjacents. Les différentes possibilités d'aimantation pour tous les volumes 3050030 16 élémentaires correspondent à un total de nmx" cas, avec nmx, le nombre de valeurs testées pour chaque volumes élémentaires (par exemple : 4 valeurs si on teste IVI_x = 0.1 1, 10 OU 100 A/m). Le problème devient très couteux en temps de calcul si on augmente le 5 nombre de cas à tester (par exemple un pas de 0.1) afin d'obtenir une meilleure résolution de recherche. Néanmoins, en restant raisonnable, on peut rapidement effectuer un premier balayage des modèles possibles et restreindre ensuite les intervalles de M_x possibles pour chaque volume élémentaire. 10 La figure 11 illustre les modèles a priori et finaux de cette étape dite 'directe'. Sur cette figure les volumes élémentaires ont la même épaisseur. On obtient alors une série de N valeurs de M_x permettant d'expliquer assez bien toutes les mesures B. Cette étape de modélisation directe (i.e. on part des paramètres M_x pour prédire les données 13_x) est avantageuse pour ensuite effectuer 15 une inversion plus rapide et robuste. Sur cette figure 11 comme sur les figures 12 et 13, la courbe en pointillé et les flèches dont le corps est en pointillé sont des valeurs estimées ou prédites et la courbe et les flèches en trait pleins sont des valeurs mesurées par la sonde à vanne de flux. La figure du haut illustre un exemple de modélisation à exclure. En effet, les 20 données prédites et les données mesurées sont trop différentes. La figure du bas illustre un exemple de modélisation possible et améliorable. En effet, les données prédites et les données mesurées sont plus proches que sur la figure du haut. 25 Etape 2 - Inversion pour retrouver une série de moments dipolaires La seconde étape consiste en effet à déterminer précisément toutes les valeurs M_x de la suite de volumes élémentaires qui expliquera le mieux les données la X, en utilisant une méthode inverse. Il s'agit de partir des données lax (entrées) pour retrouver les paramètres M_x (sorties). L'équation (1) n'étant pas inversable 30 simplement, le procédé selon l'invention prévoit d'utiliser l'équation d'un dipôle : Equation (2) 2 -1 '111 3050030 17 avec Cm, un facteur de proportionnalité égal à 10-7 SI, m_x, composante selon -x du moment dipolaire (en A.m2) du dipôle, r-x, composante selon -x du vecteur r allant du dipôle à la mesure (r = (r_x*r_ x+ry*ry+rz*rz)°5). 5 Le moment dipolaire correspond au produit de l'aimantation (en A/m) par le volume considéré (en m3, ici le volume d'un volume élémentaire rectangulaire). L'analogie avec un dipôle correspond en fait à une concentration de l'aimantation au centre de chaque volume élémentaire. Ceci est un inconvénient en considérant la réalité des séquences sédimentaires, mais le pas de mesure, donc l'espacement entre 10 deux dipôles, est suffisamment petit pour annihiler l'incohérence entre un volume élémentaire et un dipôle. De plus, le dipôle présente l'avantage d'une équation aisément dérivable, même si la relation entre et m_x reste non-linéaire. Rappelons encore que chaque mesure le long de l'échantillon U- channel correspond à la somme de plusieurs B_x, chacun lié à un dipôle. 15 L'ensemble du problème mathématique est résolu par l'utilisation de la méthode d'inversion non- linéaire généralisée fonctionnant par itérations et développée par Tarantola et Valette dans les publications suivantes : Tarantola, A. et Valette, B., 1982. Generalized Nonlinear Inverse Problems 20 Solved Using the Least Squares Criterion, Rev. Geophys. Sp. Phys., 20,219-232 Quesnel, Y., 2006. Interprétation des données magnétiques martiennes : contraintes sur l'évolution primitive de Mars. Thèse de doctorat, Université de Nantes. Cet algorithme linéarise à chaque itération l'équation (2) en utilisant son 25 développement de Taylor au r degré (dérivées partielles). Lors d'une itération, des valeurs possibles de paramètres (ici la série de m_x) sont proposées pour prédire des valeurs de données (une série de B_x prédits) qui sont sensées de mieux en mieux expliquer les véritables observations (la série de B_x mesurés). Le critère de convergence est calculé par les moindres carrés (calcul du x2). Cette approche 30 fonctionne si le nombre de données initiales est au moins égal ou supérieur au nombre de paramètres à rechercher et des contraintes sont apportées sur certains paramètres. Cela est le cas avec le procédé selon l'invention car le nombre de mesures est égal au nombre de dipôles m_x, dont nous connaissons la position parfaitement. L'équation principale du problème est : 35 3050030 18 Equation (3) + G p; 5 avec pk.i et po, ensembles des valeurs de paramètres rnx à l'itération k+1 et au départ de l'inversion, respectivement ; Gk, matrice des dérivées partielles de l'équation (1) pour l'itération k, Cdodo et Cpopo, matrices de covariance associées aux données B_x et aux paramètres m_x ; do et g(pk), ensembles des données observées et prédites par les paramètres de l'itération k précédente, respectivement. 10 Le critère des moindres carrés, x2, se calcule à chaque pas via : Equation (4) X2= 15 avec N, le nombre de données. Les itérations s'arrêtent lorsque ce critère est faible et stable depuis quelques itérations, et que les paramètres ne varient plus de façon significative. Selon un mode de réalisation l'invention prévoit différents moyens de contrôle 20 de la convergence et du temps de calcul de l'inversion. Parmi ces moyens, le procédé prévoit par exemple l'ajustement des écarts-types initiaux associés aux données lax et aux paramètres m_x, le calcul de la distribution finale des résidus entre données et prédictions, etc. Mais c'est surtout le résultat de l'étape préalable de modélisation directe qui permet de faire converger rapidement l'algorithme vers une solution 25 robuste, en proposant des valeurs de rnx au départ de l'inversion qui sont déjà quasi- idéales, en explorant seulement 'localement' les valeurs possibles. Selon un autre mode de réalisation le procédé selon l'invention, l'algorithme travaille de façon autonome en optimisant au maximum tous ces leviers de contrôle 30 afin de garantir un résultat robuste du premier coup. 3050030 19 Le produit final 'mathématique' de l'inversion correspond donc à une série de dipôles dont on connaît les positions (préalablement déterminées par la position des données) et les moments (en A.m2) selon (-x) déterminés par l'inversion. Si besoin, l'algorithme convertit alors les moments dipolaires en aimantation (en A/m) en divisant 5 chaque moment par le volume d'un volume élémentaire. La figure 12 illustre cette deuxième étape. Selon un mode de réalisation alternatif, le procédé prend en compte des 10 épaisseurs de volumes élémentaires définies par l'utilisateur. Comme indiqué aux étapes 750, 760 de la figure 7, on effectue, préalablement à l'estimation, une description préalable de la carotte formant la carotte sédimentaire. Cette description permet de délimiter parfaitement toutes les séquences de dépôts sédimentaires dans la carotte sédimentaire. Le nombre de volumes élémentaires (et donc de dipôles) sera 15 différent du nombre de mesures. On exclue de préférence les cas où la résolution d'analyse permet d'obtenir plus de volumes élémentaires que de mesures possibles. La plupart du temps ce travail visuel permet surtout de synthétiser la carotte sédimentaire en seulement quelques séquences, ce qui implique un nombre N de données largement supérieur au nombre Np de volumes élémentaires. L'étape 1 se 20 déroule comme décrite plus haut, excepté le fait que le nombre de cas à tester est réduit car nmx" devient nme. La figure 13 illustre ce mode de réalisation alternatif. La figure du bas illustre un exemple de modélisation à exclure. En effet, les données prédites et les données mesurées sont trop différentes. 25 La figure du bas illustre un exemple de modélisation possible et améliorable. En effet, les données prédites et les données mesurées sont plus proches que sur la figure du bas. L'étape 2 d'inversion est toute aussi similaire que la description présentée plus 30 haut, avec Np paramètres pour N mesures. Si les volumes élémentaires issus de l'analyse sont peu nombreux, alors le fait de focaliser toute l'aimantation en 1 seul point (dipôle) n'est pas forcément idéal, mais cela permet aussi d'augmenter la rapidité de convergence de l'inversion. 35 Au vu de la description qui précède il ressort clairement que l'invention permet de fournir une estimation fiable, reproductible et simple de l'aimantation rémanente de 3050030 20 carottes sédimentaires, en particulier de carottes sédimentaires présentant une forte aimantation. Un autre objet de la présente invention est de proposer telle solution pour tout en réduisant le coût 5 L'invention n'est pas limitée aux modes de réalisations précédemment décrits et s'étend à tous les modes de réalisation couverts par les revendications. En particulier elle s'étend à tous les types de carottes sédimentaires qui s'étendent principalement selon une direction privilégiée. La section de la carotte le 10 long de cet axe peut être constante ou non. 3050030 21 REFERENCES 1 Banc de mesure 10 Cadre 5 11 Poutre horizontale supérieure 12 Poutre horizontale inférieure 13 Profilés verticaux 21 Courroie de transport 10 22 Moteur 23 Poulie de renvoi supérieur 24 Poulie de renvoi inférieur 25 Chariot 26 Guide 15 30 Chambre amagnétique 31 Bouclier magnétique 32 Couche interne 33 Couche externe Sonde à vanne de flux 20 41 Face principale 50 Carte de commande 100 Carotte sédimentaire 101 Portion d'échantillon ou volume élémentaire REVENDICATIONS1 Banc de mesure (1) du champ magnétique d'une carotte sédimentaire (100) s'étendant selon une direction principale, caractérisé en ce que le banc de mesure (1) comprend au moins: - un dispositif de déplacement en translation de la carotte sédimentaire (100) au sein du banc de mesure (1), - une chambre de mesure (30) comprenant une sonde à vanne de flux (40) agencée pour être disposée en regard et de préférence au contact d'une portion (101) de carotte sédimentaire (100) de manière mesurer, perpendiculairement à ladite direction principale, le champ magnétique généré par ladite portion (101) de carotte sédimentaire (100) située en regard de la sonde à vanne de flux (40); et ce que le dispositif de déplacement est configure de manière à faire translater la carotte sédimentaire (100) par rapport à la sonde à vanne de flux (40) de manière à ce que la sonde à vanne de flux (40) mesure successivement le champ magnétique d'une succession de portions (101) de la carotte sédimentaire (100). 2. Banc de mesure (1) selon la revendication précédente dans lequel la sonde à vanne de flux (40) est agencée de manière à mesurer perpendiculairement à ladite direction principale le champ magnétique généré par ladite portion (101) de carotte sédimentaire (100) située en regard de la sonde à vanne de flux (40). 3. Banc de mesure (1) selon l'une quelconque des revendications précédentes dans lequel la sonde à vanne de flux (40) est située verticalement en dessous de la carotte sédimentaire (100) lorsque la carotte sédimentaire (100) est située en partie au moins dans la chambre de mesure (30). 4. Banc de mesure (1) selon l'une quelconque des revendications précédentes dans lequel la prise de mesure par la sonde à vanne de flux (40) est synchronisée avec le déplacement en translation de la carotte sédimentaire (100) de manière à déclencher une mesure par la sonde à vanne de flux (40) lorsque qu'un déplacement de la carotte sédimentaire est achevé pour amener un volume élémentaire (101) au droit de la sonde à vanne de flux (40). 3050030 23 5. Banc de mesure (1) selon l'une quelconque des revendications précédentes dans lequel la chambre de mesure (30) comprend un bouclier (31) magnétique formant un manchon à l'intérieur duquel sont situées la sonde à vanne de flux (40) et ladite portion (101) de carotte sédimentaire (100) située en regard de la 5 sonde à vanne de flux (40). 6. Banc de mesure (1) selon la revendication précédente dans lequel le bouclier (31) magnétique comprend au moins un manchon et de préférence deux manchons concentriques en mu-métal. 7 Banc de mesure (1) selon l'une quelconque des revendications précédentes dans lequel le dispositif de déplacement comprend une courroie de transport (21), entrainée en déplacement par un moteur (22), et dans lequel la courroie de transport (21) traverse la chambre de mesure (30). 8. Banc de mesure (1) selon la revendication précédente dans lequel le moteur (22) est un moteur pas à pas. 9. Appareil d'estimation des aimantations rémanentes de portions (101) d'un 20 échantillon formant une carotte sédimentaire (100) de section rectangulaire (« u-channel ») et s'étendant selon une direction principale, l'appareil comprenant un banc de mesure (1) selon l'une quelconque des revendications précédentes et un module de calcul agencé pour mettre en oeuvre au moins les étapes suivantes : 25 Définition d'une succession de N volumes élémentaires (101) et de leurs dimensions, chaque volume élémentaire (101) correspondant à une portion de carotte sédimentaire (100), les volumes élémentaires (101) étant adjacents et successivement disposés selon ladite direction principale de la carotte sédimentaire (100); 30 Mesure par la sonde à vanne de flux du champ magnétique produit par chaque volume élémentaire (101) ; En utilisant les mesures de champ magnétique de la sonde à vanne de flux : détermination d'une aimantation rémanente M_x propre à chaque volume élémentaire (101), cette étape comprenant; 35 - une première identification « directe » d'une pluralité de valeurs possibles d'aimantation rémanente M_x, cette première identification comprenant: 10 15 3050030 24 Le calcul de valeurs de champ magnétique B_, à partir d'une pluralité de valeurs d'aimantations rémanentes M_, choisies arbitrairement ; l'identification, parmi ladite pluralité de valeurs d'aimantations rémanentes M_, choisies arbitrairement, d'une ou plusieurs valeurs 5 d'aimantation rémanente M_, pour lesquelles le champ magnétique B_, calculé est compris entre 0.7 et 1.3 fois, de préférence entre 0.8 et 1.2 fois, et de préférence entre 0.9 et 1.1 fois, le champ magnétique B_, mesuré par la sonde à vanne de flux ; - une deuxième identification « inverse » d'une ou plusieurs valeurs 10 d'aimantation rémanente M_x, cette deuxième identification comprenant: - le calcul d'une valeur M_x d'aimantation rémanente en fonction du champ magnétique mesuré par la sonde à vanne de flux. 10. Procédé d'estimation des aimantations rémanentes de portions d'une carotte 15 sédimentaire (100) s'étendant selon une direction principale, le procédé comprenant au moins les étapes suivantes: i. Définition d'une succession de N volumes élémentaires (101) et de leurs dimensions, chaque volume élémentaire (101) correspondant à une portion de carotte sédimentaire (100), les volumes élémentaires (101) étant adjacents et 20 successivement disposés selon ladite direction principale de la carotte sédimentaire (100); ii. Translation de la carotte sédimentaire (100) relativement à une sonde à vanne de flux de manière à faire passer chaque volume élémentaire (101) en regard de la sonde à vanne de flux ; 25 iii. Mesure par la sonde à vanne de flux du champ magnétique produit par chaque volume élémentaire (101) ; iv. En utilisant les mesures de champ magnétique de la sonde à vanne de flux : détermination d'une aimantation rémanente M_x propre à chaque volume élémentaire (101), cette étape comprenant; 30 - une première identification « directe » d'une pluralité de valeurs possibles d'aimantation rémanente M-x, cette première identification comprenant: Le calcul de valeurs de champ magnétique B_x à partir d'une pluralité de valeurs d'aimantations rémanentes M_x choisies arbitrairement ; l'identification, parmi ladite pluralité de valeurs d'aimantations 35 rémanentes M_x choisies arbitrairement, d'une ou plusieurs valeurs d'aimantation rémanente M_x pour lesquelles le champ magnétique B_x 3050030 25 calculé est compris entre 0.7 et 1.3 fois, de préférence entre 0.8 et 1.2 fois, et de préférence entre 0.9 et 1.1 fois, le champ magnétique B_, mesuré par la sonde à vanne de flux ; - une deuxième identification « inverse » d'une ou plusieurs valeurs 5 d'aimantation rémanente M_x, cette deuxième identification comprenant: - le calcul d'une valeur M_x d'aimantation rémanente en fonction du champ magnétique mesuré par la sonde à vanne de flux. 11. Procédé selon la revendication précédente dans lequel ladite étape de 10 définition est effectuée de manière à ce que chaque volume élémentaire (101) de carotte sédimentaire (100) présente une même épaisseur, l'épaisseur étant mesurée selon ladite direction principale. 12. Procédé selon la revendication 10 dans lequel ladite étape de définition est 15 effectuée de manière à ce que chaque volume élémentaire (101) de carotte sédimentaire (100) présente une épaisseur indépendante de celle des autres volumes élémentaires (101), l'épaisseur étant mesurée selon ladite direction principale. 20 13. Procédé selon la revendication précédente dans lequel l'épaisseur est déterminée par analyse visuelle de la carotte selon un critère d'homogénéité. 14. Procédé selon la revendication précédente dans lequel l'épaisseur le critère d'homogénéité est pris parmi : la couleur et la texture. 25 15. Procédé selon l'une quelconque des revendications 10 à 14 comprenant une étape préalable de fourniture d'une carotte sédimentaire (100) présentant des aimantations rémanentes naturelles, puis une étape de transformation des aimantations rémanentes en aimantations rémanentes artificielles, ladite étape 30 de transformation comprenant un passage de la carotte sédimentaire (100) dans un anneau de Halbach. 16. Produit programme d'ordinateur comprenant des instructions, qui lorsqu'elles sont effectuées par au moins un processeur, exécute au moins l'étape iv. du 35 procédé selon l'une quelconque des revendications 10 à 15.