Domaine technique

La présente invention concerne le domaine de la science du sol et des géosciences de l’environnement. Plus particulièrement la caractérisation du carbone contenu dans les formations superficielles, et notamment dans les sols.

Afin de répondre aux enjeux écologiques, ou de respecter certaines législations/directives environnementales, les acteurs de la science du sol et des géosciences de l’environnement (laboratoires de recherche, bureaux d’étude, agences environnementales, exploitants agricoles) sont de plus en plus souvent amenés à mettre en place des protocoles de suivi des impacts des activités humaines sur les stocks de carbone des sols et des éco-agro- systèmes. Ces suivis et ces études d’impact impliquent de pouvoir étudier de larges séries d’échantillons dans des délais relativement courts au regard des méthodologies classiquement employées. De plus, ces méthodes s’accompagnent souvent de contraintes environnementales et sécuritaires qui augmentent les délais et les coûts analytiques, et nécessitent souvent de faire appel à des prestataires spécialisés (par exemple des laboratoires d’analyse).

Les formes organiques du carbone stockées dans les formations superficielles, et notamment dans les sols, représentent un enjeu majeur pour l’agriculture et le climat. Elles jouent en effet un rôle primordial sur la qualité structurelle et la valeur fertilisante des sols, mais interviennent surtout dans le cycle du carbone, en représentant le plus important réservoir de carbone organique à la surface de la Terre. Assurer la sécurité alimentaire des populations et favoriser la séquestration du carbone organique dans les sols nécessite de trouver un compromis entre maintenir des stocks de matière organique suffisamment labiles pour libérer les nutriments nécessaires à la croissance végétale tout en favorisant le stockage de formes de carbone suffisamment résistantes pour atténuer, à long terme, les émissions anthropogéniques de gaz à effet de serre.

Les formes inorganiques de carbone sont de minéralogie variable. Essentiellement représentées par des carbonates et des oxalates de calcium (Ca) dans les sols, elles peuvent parfois intégrer d’autres cations (Fe dans la sidérite, Mg dans la dolomite, etc.), notamment en domaine tropical. Indépendamment de leur rôle dans les cycles biogéochimiques (notamment comme puits de carbone), les formes minérales du carbone posent des problèmes techniques pour l’analyse des formes organiques. En effet, aucune méthode de routine ne permet de caractériser les différentes formes de carbone à l’aide d’une seule et unique mesure, et les approches classiquement employées reposent sur des extractions et des séparations spécifiques (fumigation acide, méthode Walkley and Black, calcimètrie Bernard) suivies d’analyses ciblées à l’aide des équipements scientifiques spécialisés (analyse élémentaire CHN, spectrométrie IR, GC-MS, ¹³ C-RMN).

On connaît aussi des méthodes d'analyse thermique de la matière organique des sols reposant sur des mesures de quantités de composés hydrocarbonés (HC), de monoxyde de carbone (CO) et/ou de dioxyde de carbone (CO2) libérées au cours du temps par un échantillon soumis à une séquence de chauffe en atmosphère inerte puis/ou à une séquence de chauffe en atmosphère oxydante. Ces méthodes ont été développées initialement dans le domaine de l'industrie pétrolière, à des fins de caractérisation de la fraction organique des roches sédimentaires. On connaît ainsi le dispositif ROCK-EVAL® (IFP Energies nouvelles, France), développé par la demanderesse et décrit notamment dans les brevets FR 2227797 (US 3953171) et FR 2472754 (US 4352673), qui comprend un four de pyrolyse distinct d'un four d'oxydation, un détecteur du type à ionisation de flamme (FID) pour détecter les composés hydrocarbonés (HC) et un détecteur du type infrarouge (IR) pour détecter le monoxyde de carbone (CO) et/ou le dioxyde de carbone (CO2). On connaît en outre des procédés développés pour des applications particulières du domaine de l'industrie pétrolière, ayant chacun leur propre séquence de températures de chauffe en pyrolyse et/ou de chauffe en atmosphère oxydante. Notamment, on connaît la méthode "ROCK-EVAL® BULK ROCK", dédiée plus particulièrement aux échantillons conventionnels de roches mères (Behar et al., 2001). La séquence de chauffe en atmosphère inerte de cette méthode est caractérisée par une température initiale T1 du four de pyrolyse généralement comprise entre 300°C et 350°C, température qui est maintenue pendant une durée prédéterminée de quelques minutes. C’est durant cette phase que sont libérés les hydrocarbures dits « libres » (correspondant en réalité à des hydrocarbures de poids moléculaire léger à lourd) initialement contenus dans l’échantillon de roche. Leur quantité est estimée via la mesure de la surface d’un premier pic, noté Si, de la courbe (aussi appelée thermogramme) représentant la quantité de composés hydrocarbonés libérés pendant la séquence de chauffe en atmosphère inerte. Puis, la température de pyrolyse est augmentée progressivement jusqu’à une température T2, généralement de 650°C. Durant cette phase, on assiste à la volatilisation des composés hydrocarbonés très lourds, ainsi qu’au craquage de la matière organique non volatile (kérogène). La quantité de composés hydrocarbonés libérés durant cette phase de craquage thermique est estimée via la mesure de la surface d’un second pic, noté S2. En parallèle, les quantités de CO et de CO2 sont mesurées et représentées également sous la forme de courbes. Ces courbes présentent deux pics classiquement notés S3CO (respectivement S3C02), qui est considéré comme correspondant au CO (respectivement CO2) généré par le craquage de la matière organique de l’échantillon pendant la chauffe sous atmosphère inerte, et S3’CO (respectivement S3'C02) qui est considéré comme correspondant au CO (respectivement CO2) généré par la décomposition thermique des formes carbonatées (notamment la calcite) pendant la chauffe sous atmosphère inerte. Puis le résidu de l'échantillon issu de la chauffe en atmosphère inerte est soumis à une chauffe en atmosphère oxydante : à partir d’une température comprise entre environ 300°C et 400°C, et valant de préférence 300°C, on élève la température du résidu de l’échantillon considéré selon un gradient de température compris entre 20 et 40°C/minute, jusqu’à une température de fin d’oxydation comprise entre 750 et 950°C, et valant de préférence 850°C. Durant cette séquence de chauffe en atmosphère oxydante, les quantités de CO et de CO2 libérées par le résidu de l'échantillon sont mesurées et représentées sous la forme de courbes, conduisant à un pic classiquement noté S4CO (respectivement S4C02) qui est considéré comme correspondant à la quantité de CO (respectivement CO2) généré par la combustion de la matière organique pendant le cycle d’oxydation. A partir de ces mesures, cette méthode définit un certain nombre de paramètres standards, en particulier le paramètre noté TOC (pour "Total Organic Carbon" en anglais) qui correspond à la teneur en carbone de l'échantillon, déterminée à partir de la quantité totale de HC libérée par l'échantillon et des quantités de CO et de CO2 libérées en-dessous de températures-seuils pendant la phase de pyrolyse et la phase d'oxydation ; et le paramètre noté MinC (pour "Minerai Carbon" en anglais) qui correspond à la teneur en carbone minéral de l'échantillon, déterminée à partir des quantités de CO et de CO2 libérées par l'échantillon au-dessus de températures-seuils pendant la phase de pyrolyse et la phase d'oxydation.

Technique antérieure

Les documents suivants seront cités au cours de la description :

Behar F., Beaumont V., De B., Penteado H.L. (2001) Rock-Eval 6 Technology: Performances and Developments, Oil & Gas Science and Technology 56, 111-134.

Disnar, J. R., Guillet, B., Keravis, D., Di-Giovanni, C., Sebag, D., 2003. Soil organic matter (SOM) characterization by Rock-Eval pyrolysis: scope and limitations. Organic Geochemistry 34, 327-343.

Malou, O. P., Sebag, D., Moulin P., Chevallier, T., Badiane-Ndour, N.Y., Thiam, A., Chapuis- Lardy, L. 2020. The Rock-Eval® signature of soil organic carbon in Arenosols of the Senegalese groundnut Basin. How do agricultural practices matter? » Agriculture, Ecosystems & Environment 301 : 107030. https://doi.Org/10.1016/j.agee.2020.107030. Pillot, D., Deville, E., Prinzhofer, A., 2014. Identification and Quantification of Carbonate Species Using Rock-Eval Pyrolysis. Oil & Gas Science and Technology - Revue d’IFP Energies nouvelles 69, 341-349.

Sebag, D., Disnar, J. R., Guillet, B., Di Giovanni, C., Verrecchia, E.P., Durand, A., 2006. Monitoring organic matter dynamics in soil profiles by “Rock-Eval pyrolysis”: bulk characterization and quantification of dégradation. European Journal of Soil Science 57, 344-355.

Sebag, D., Verrecchia, E.P., Cécillon, L., Adatte, T., Albrecht, R., Aubert, M., Bureau, F., Cailleau, G., Copard, Y., Decaëns, T., Disnar, J. -R., Hetényi, M., Nyilas, T., Trombino, L., 2016. Dynamics of soil organic matter based on new Rock-Eval indices. Geoderma 284, 185-203.

Depuis les années 2000, ce type d'analyses a également été utilisé et adapté pour étudier la fraction organique des formations superficielles et notamment des sols.

Par exemple on connaît la méthode décrite dans le document (Disnar et al., 2003), dans lequel est proposée une adaptation de la séquence de chauffe en pyrolyse consistant en une température initiale de la séquence de chauffe en atmosphère inerte valant 200°C, au lieu de 300°C dans la méthode "ROCK-EVAL® BULK ROCK" décrite précédemment. Cette température initiale plus basse permet de d’étendre le pic S2 défini précédemment aux températures de craquage des composants organiques les plus thermolabiles. Ces constituants thermolabiles sont en effet beaucoup plus abondants dans les formations superficielles que dans les roches sédimentaires. Pour estimer leur contribution et ainsi évaluer la stabilité thermique de la matière organique, cette méthode définit le paramètre R400 qui mesure la proportion relative du pic S2 correspondant aux températures inférieures à 400°C. Par ailleurs, ce même document met en évidence un écart négatif entre le paramètre TOC tel que défini dans la méthode ROCK-EVAL® BULK ROCK et les teneurs en carbone organique mesurées avec des méthodes normalisées (par exemple des analyses élémentaires) et recommande une correction statistique (c'est-à-dire l'application d’un coefficient correcteur établi sur un panel représentatif d’échantillons de sols) pour corriger le paramètre TOC, défini pour le domaine pétrolier, afin qu'il soit vraiment représentatif de la teneur en carbone organique, notamment pour les échantillons riches en matières organiques peu décomposées par rapport à leurs précurseurs biogéniques (litières, composts, tourbes, etc.). On connaît également le document (Sebag et al., 2006) qui utilise cette séquence de chauffe adaptée, et qui propose en complément une déconvolution du pic S2 pour évaluer le degré de décomposition des constituants organiques. En effet, les thermogrammes des formations superficielles, en particulier des échantillons organiques, présentent une distribution plurimodale dont les modes principaux sont toujours situés dans des intervalles de température particuliers (300-320°C ; 360-380°C ; 420-440°C ; 470-490°C et 540-560°C). La méthode décrite dans ce document consiste à décomposer le pic S2 en cinq distributions gaussiennes élémentaires centrées sur ces modes. Cette déconvolution mathématique revient à estimer la contribution relative au pic S2 que chaque distribution élémentaire considérée comme relative à une classe de constituants définis par leur seule température de craquage. Ces contributions élémentaires sont ensuite utilisées pour calculer de nouveaux paramètres qui mesurent la stabilité thermique globale de la matière organique (R- index) et le degré de décomposition de la fraction thermolabile (l-index). Cependant, la méthode de déconvolution repose sur une approche itérative pour ajuster la position et la surface de chaque distribution élémentaire en fonction de paramètres statistiques fixés arbitrairement, ce qui réduit la reproductibilité de la décomposition.

On connaît également le document (Sebag et al., 2016) qui décrit une approche alternative en montrant qu’une intégration des thermogrammes représentatifs de la quantité de HC contenus dans l'échantillon par bandes de température (200-340°C ; 340-400°C ; 400- 460°C ; 460-520°C et 520-650°C) aboutit à des résultats comparables à la déconvolution mathématique, et sont, qui plus est, parfaitement reproductibles. Cette approche a été utilisée pour de nombreuses applications, mais souffre de quelques limites intrinsèques. Tout d’abord, elle repose sur la définition de températures-seuils définies empiriquement pour que la surface de chaque bande de température approche au mieux la surface de la distribution élémentaire correspondante obtenue par déconvolution. En outre, ces surfaces et les paramètres qui en découlent ne fournissent qu’une information qualitative dans la mesure où les unes comme les autres correspondent à des proportions relatives ou des rapports de proportions relatives. Enfin, cette approche se limite à l’examen des quantités de HC émises au cours de la phase de pyrolyse, alors que la majorité du carbone est représentée par le CO et le CO2 émis au cours des phases de pyrolyse et d’oxydation.

En outre, malgré ces tentatives d'adaptation ou d’amélioration, les méthodes connues ne fournissent pas des valeurs des paramètres standards TOC et MINC qui soient représentatives respectivement de la teneur d'un échantillon en carbone organique et en carbone minéral dans le cas de formations superficielles. En effet, le paramètre TOC (respectivement MINC) de la méthode "ROCK-EVAL® BULK ROCK" mesure la teneur en carbone à partir des quantités de CO et de CO2 émises en-dessous (respectivement au- dessus) de températures-seuils considérées comme des limites thermiques entre les formes organiques et minérales de carbone.

Or cette définition est erronée que ce soit pour des formations carbonatées ou non carbonatées. En effet, comme montré dans le document (Malou et al. 2020) cette définition implique notamment le calcul systématique d’un paramètre MINC, même pour les échantillons sans forme minérale de carbone. Mais plus généralement, comme constaté par la demanderesse, ceci implique que, même dans les échantillons contenant des formes minérales de carbone, une part du MINC provient du craquage thermique des constituants organiques au-dessus des températures-seuils. Enfin, comme montré dans le document (Pillot et al. 2014), certaines formes minérales de carbone (comme la sidérite ou les oxalates) sont susceptibles de se décomposer à des températures inférieures aux températures-seuils utilisées pour les calculs des paramètres TOC et de MINC standards.

La présente invention vise à pallier les inconvénients liés à l’utilisation d'un procédé initialement développé pour le domaine pétrolier pour analyser la matière organique des formations superficielles, et notamment des sols. En particulier, la présente invention vise à obtenir des thermogrammes présentant des pics mieux différenciés et communs à chaque type d’effluents (HC, CO et CO2). Cette invention permet ainsi de mieux discriminer les différentes classes de composés en fonction de leur résistance thermique. Elle conduit également à statuer sur la nature organique et/ou minérale des formes de carbone et ainsi à définir des nouveaux paramètres plus représentatifs de la teneur en carbone organique et en carbone minéral d'échantillons analysés par des méthodes normalisées.

Plus généralement, la présente invention permet une caractérisation rapide et fiable des formes de carbone présentes dans un échantillon pour déterminer si elles contiennent des formes de nature minérale ou non, pour préciser le degré de décomposition et/ou la stabilité thermique de la part de matière organique la plus réactive (liée à l’atome d’hydrogène), et pour quantifier plus précisément les teneurs en carbone (organique versus minéral, thermolabile versus thermostable).

Résumé de l’invention

La présente invention concerne un procédé pour caractériser et quantifier le carbone présent dans une formation superficielle, à partir d’un échantillon représentatif de ladite formation superficielle, dans lequel on applique au moins les étapes suivantes audit échantillon : A. on chauffe ledit échantillon selon une première séquence de chauffe sous atmosphère inerte, et on mesure en continu une quantité de composés hydrocarbonés, une quantité de CO et une quantité de CO2 libérés pendant ladite première séquence de chauffe, ladite première séquence de chauffe sous atmosphère inerte comprenant au moins une succession d'au moins six paliers isothermes d'une durée prédéterminée, ladite succession desdits au moins six paliers isothermes comprenant un premier palier isotherme à une première température comprise entre 80 et 200°C (TO), un deuxième palier isotherme à une deuxième température comprise entre 340 et 380°C (T1), un troisième palier isotherme à une troisième température comprise entre 400 et 440°C (T2), un quatrième palier isotherme à une quatrième température comprise entre 450 et 490°C (T3), un cinquième palier isotherme à une cinquième température comprise entre 500 et 540°C (T4), un sixième palier isotherme à une sixième température comprise entre 580 et 650°C (T5), lesdits paliers isothermes étant reliés entre eux par un gradient thermique ;

B. on chauffe un résidu dudit échantillon issu de ladite première séquence de chauffe selon une seconde séquence de chauffe sous atmosphère oxydante, et on mesure une quantité de CO et une quantité de CO2 libérés pendant ladite seconde séquence de chauffe, ladite deuxième séquence de chauffe sous atmosphère oxydante débutant à une température minimale comprise entre 150 et 300°C, se terminant à une température maximale comprise entre 850 et 1200°C, et suivant un gradient thermique

C. on détermine une valeur d’un paramètre SCmin représentatif d’une proportion en carbone minéral par rapport au carbone total dudit échantillon à partir d’un ratio entre ladite quantité de CO2 libéré par ledit résidu dudit échantillon au-delà d’une température intermédiaire de ladite deuxième séquence de chauffe comprise entre 620 et 680°C, et ladite quantité de CO2 libéré par ledit résidu pendant ladite deuxième séquence de chauffe ;

D. on quantifie une teneur en carbone organique et/ou une teneur en carbone minéral dudit échantillon de la manière suivante : i. si ladite valeur dudit paramètre SCmin est inférieure ou égale à une valeur seuil prédéfinie dudit paramètre SCmin, ladite teneur en carbone minéral est nulle et ladite teneur en carbone organique est fonction de la totalité desdites quantités de HC, de CO et de C02 libérés pendant ladite première séquence de chauffe, et desdites quantités de CO et de C02 libérés pendant ladite deuxième séquence de chauffe jusqu'à ladite température intermédiaire de ladite deuxième séquence de chauffe. ii. si ladite valeur dudit paramètre SCmin est supérieure à ladite valeur seuil prédéfinie dudit paramètre SCmin :

- on estime une quantité de CO2 libéré par craquage thermique de la matière organique dudit échantillon pendant ladite première séquence de chauffe au-delà de ladite troisième température (T2) à partir de ladite quantité de CO2 libéré pendant ladite première séquence de chauffe jusqu'à ladite troisième température (T2) de ladite première séquence de chauffe, multipliée par un facteur k prédéfini ; et/ou

- on détermine ladite teneur en carbone organique et/ou ladite teneur en carbone minéral en fonction desdites quantités de HC, de CO et de CO2 libérés pendant ladite première séquence de chauffe et desdites quantités de CO et de CO2 libérés pendant ladite deuxième séquence de chauffe jusqu'à ladite température intermédiaire de ladite deuxième séquence de chauffe, et en prenant en compte ladite quantité estimée CO2 libéré par craquage thermique de la matière organique dudit échantillon pendant ladite première séquence de chauffe au-delà de ladite troisième température (T2).

Selon une mise en oeuvre de l’invention, si ladite valeur dudit paramètre SCmin est supérieure à ladite valeur seuil prédéfinie dudit paramètre SCmin :

- ladite teneur en carbone organique peut être fonction de la totalité de ladite quantité de HC libérés pendant ladite première séquence de chauffe, de ladite quantité de CO libéré pendant ladite première séquence de chauffe jusqu'à ladite cinquième température (T4) à laquelle on ajoute la moitié de ladite quantité de CO libéré pendant ladite première séquence de chauffe au-delà de ladite cinquième température (T4), de ladite quantité de CO2 libéré pendant ladite première séquence de chauffe jusqu'à ladite troisième température (T2) à laquelle on ajoute ladite quantité estimée de CO2 libéré par craquage thermique de la matière organique dudit échantillon pendant ladite première séquence de chauffe au-delà de ladite troisième température (T2), et desdites quantités de CO et de C02 libérés pendant ladite deuxième séquence de chauffe jusqu'à ladite température intermédiaire de ladite deuxième séquence de chauffe, et/ou

- ladite teneur en carbone minéral peut être fonction de ladite quantité de CO2 libéré pendant ladite première séquence de chauffe au-delà de ladite troisième température (T2) de ladite première séquence de chauffe à laquelle on retire ladite quantité estimée de CO2 libéré par craquage thermique de la matière organique dudit échantillon pendant ladite première séquence de chauffe au-delà de ladite troisième température (T2), de ladite moitié de ladite quantité de CO libéré pendant ladite première séquence de chauffe au-delà de ladite cinquième température (T4), et de ladite quantité de CO2 libéré pendant ladite deuxième séquence de chauffe au-delà de ladite température intermédiaire de ladite deuxième séquence de chauffe.

Selon une mise en oeuvre de l’invention, ledit facteur k prédéfini peut être compris entre 1.3 et 1 .4, et vaut de préférence 1 .3724.

Selon une mise en oeuvre de l’invention, au moyen d’une pluralité d’échantillons provenant de formations superficielles non carbonatées, on peut appliquer au moins ladite première séquence de chauffe à chacun desdits échantillons de ladite pluralité d’échantillons, on peut mesurer une quantité de CO2 libéré par chacun desdits échantillons pendant ladite première séquence de chauffe, et on peut déterminer ledit facteur k au moyen d’une régression linéaire appliquée entre lesdites quantités de CO2 libérés par ladite pluralité d’échantillons jusqu’à ladite troisième température (T2) de ladite première séquence de chauffe et lesdites quantités de CO2 libérés par ladite pluralité d’échantillons au-delà de ladite troisième température (T2) de ladite première séquence de chauffe.

Selon une mise en oeuvre de l’invention, ladite valeur seuil prédéfinie dudit paramètre SCmin peut être comprise entre 0.03 et 0.05 et vaut de préférence 0.04.

Selon une mise en oeuvre de l’invention, ladite première température (TO) de la dite première séquence de chauffe peut être comprise entre 80°C et 150°C, et vaut de préférence 150°C.

Selon une mise en oeuvre de l’invention, ladite température minimale de ladite deuxième séquence de chauffe peut être égale à ladite première température (TO) de la dite première séquence de chauffe.

Selon une mise en oeuvre de l’invention, ladite durée prédéterminée d’un desdits paliers isothermes de ladite première séquence de chauffe peut être comprise entre 3 et 5 minutes.

Selon une mise en oeuvre de l’invention, l’un desdits gradients thermiques de ladite première séquence de chauffe peut être compris entre 1 et 50 ^o C.min ^_1 , de préférence entre 20 et 25 ⁰ C.min ^_1 .

Selon une mise en oeuvre de l’invention, ledit gradient thermique de ladite deuxième séquence de chauffe peut être compris entre 20 et 40 ^o C.min ^_1 , de préférence entre 20 et 25 ⁰ C.min ^_1 .

Selon une mise en oeuvre de l’invention, si ladite valeur dudit paramètre SCmin est inférieure ou égale à ladite valeur seuil prédéfinie dudit paramètre SCmin : ladite teneur en carbone minéral peut être définie par un paramètre MinCfs tel que MinCfs=0 ; et/ou

- ladite teneur en carbone organique peut être définie par un paramètre TOCfs s’exprimant selon une formule du type :

TOCfs = (S2 _t x 0.083) + (S3C02 _t x— ) + (S3CO _t x — )+(S4C02 _a x — ) + (S4CO _a x—) où S2 _t , S3C02 _t , S3CO _t représentent lesdites quantités respectivement de HC, de CO ₂ et de CO libérés pendant ladite première séquence de chauffe, S4C02 _a et S4CO _a représentent respectivement lesdites quantités de CO ₂ et de CO libérés pendant ladite deuxième séquence de chauffe et jusqu’à ladite température intermédiaire de ladite deuxième séquence de chauffe.

Selon une mise en oeuvre de l’invention, si ladite valeur dudit paramètre SCmin est supérieure à ladite valeur seuil prédéfinie dudit paramètre SCmin :

- ladite teneur en carbone minéral peut être définie par un paramètre MinCfs s’exprimant selon une formule du type : où S3C02 _b représente ladite quantité de CO ₂ libéré pendant ladite première séquence de chauffe au-delà de ladite troisième température (T2), S3C02 _c représente ladite quantité estimée de CO ₂ libéré par craquage thermique de la matière organique dudit échantillon pendant ladite première séquence de chauffe au-delà de ladite troisième température (T2), S3CO _b représente ladite quantité de CO libéré pendant ladite première séquence de chauffe au-delà de ladite cinquième température (T4), et S4C02 _b représente ladite quantité de CO ₂ libéré pendant ladite deuxième séquence de chauffe au-delà de ladite température intermédiaire de ladite deuxième séquence de chauffe ; et/ou ladite teneur en carbone organique peut être définie par un paramètre TOCfs s’exprimant selon une formule du type : où S2 _t représente ladite quantité de HC libérés pendant ladite première séquence de chauffe, S3C02 _a représente ladite quantité de CO ₂ libéré pendant ladite première séquence de chauffe jusqu'à ladite troisième température (T2), S3C02 _c représente ladite quantité estimée de CO ₂ libéré par craquage thermique de la matière organique dudit échantillon pendant ladite première séquence de chauffe au-delà de ladite troisième température (T2), S3CO _a et S3CO _b représentent lesdites quantités de CO libérés pendant ladite première séquence de chauffe respectivement jusqu'à et au-delà de ladite cinquième température (T4), S C02 _a et S CO _a représentent respectivement lesdites quantités de CO2 et de CO libérés pendant ladite deuxième séquence de chauffe et jusqu’à ladite température intermédiaire de ladite deuxième séquence de chauffe.

D'autres caractéristiques et avantages du procédé selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après.

Liste des figures

La figure 1 illustre de manière schématique la séquence de chauffe sous atmosphère inerte selon l’invention.

Les figures 2A et 2B illustrent des thermogrammes obtenus au moyen d’une mise en oeuvre de la séquence de chauffe sous atmosphère inerte selon l’invention, obtenus respectivement pour un échantillon de sol non carbonaté et pour un échantillon de sol carbonaté.

Les figures 2C et 2D illustrent des thermogrammes obtenus au moyen d'une séquence de chauffe selon l'art antérieur, obtenus respectivement pour les échantillons de sol non carbonaté et de sol carbonaté des figures 2A et 2B.

Les figures 3A et 3B illustrent des thermogrammes obtenus au moyen d’une mise en oeuvre de la séquence de chauffe sous atmosphère oxydante selon l’invention, obtenus respectivement pour l'échantillon de sol non carbonaté des figures 2A et 2C et pour l'échantillon de sol carbonaté des figures 2B et 2D.

Les figures 4A et 4B présentent des histogrammes comparant les valeurs respectivement des teneurs en carbone organique Corg et des teneurs en carbone minéral Cmin déterminées par le procédé selon l'invention et par l'art antérieur pour 9 échantillons de sols.

Les figures 5A et 5B présentent des histogrammes comparant les valeurs respectivement des teneurs en carbone organique Corg et des teneurs en carbone minéral Cmin déterminées par le procédé selon l'invention et par l'art antérieur pour 7 échantillons de tourbes.

Description des modes de réalisation

L’invention concerne un procédé pour caractériser et quantifier le carbone (ou autrement dit les formes de carbone) présentes dans une formation superficielle. Par "formation superficielle" ("superficial deposit" en anglais), on entend une formation continentale ou littorale, meuble ou secondairement consolidée, provenant de la désagrégation mécanique et/ou chimique de roches préexistantes, et formée à l’interface lithosphère/biosphère/atmosphère. On distingue (i) les "formations superficielles allochtones" (comme les colluvions, alluvions, lœss, etc.) qui ont subi ou subissent encore des déplacements proches ou lointains, et ne reposent plus sur leur matériel-parent, et les "formations superficielles autochtones" (comme les arènes, altérités, argiles à silex, etc.) qui ont évolué sur place à partir d’un matériel parent qui constitue encore leur substrat.

Par sol ("soil" en anglais), on entend l’ensemble des couches externes des formations superficielles, dont les propriétés sont directement contrôlées par les actions mutuelles de l’eau, de l’air et des organismes vivants et morts, voire des activités humaines pour les périodes les plus récentes.

Ces termes sont notamment définis dans le document de référence (Dictionnaire encyclopédique de Science du Sol, Mathieu & Lozet, Lavoisier, 2011).

Le procédé selon l'invention requiert de disposer d'au moins un échantillon représentatif de la formation superficielle : cet échantillon peut avoir été prélevé manuellement dans une fosse ou par carottage à l’aide d’une tarière. Avantageusement, l’échantillon tel que prélevé est tamisé au moyen d'un tamis dont les orifices ont un diamètre de 2 mm, séché à une température inférieure à 40°C, puis broyé jusqu'à obtenir des fragments ayant des dimensions inférieures à 200 pm.

Le procédé selon l’invention peut être avantageusement mais non limitativement mis en œuvre au moyen du dispositif ROCK-EVAL® (IFP Energies nouvelles, France), tel que décrit dans les brevets FR 2227797 (US 3953171) et FR 2472754 (US 4352673). En effet, le dispositif ROCK-EVAL® comprend au moins :

- un four de pyrolyse en atmosphère non oxydante,

- des moyens de transfert des résidus de pyrolyse dans un four d'oxydation,

- un four d’oxydation en atmosphère oxydante,

- des moyens de mesure de la quantité de composés hydrocarbonés (HC) libérés au cours de la pyrolyse,

- des moyens de mesure du monoxyde de carbone (CO) et du dioxyde de carbone (C0 ₂ ).

Le procédé peut également être mis en œuvre au moyen d’un unique four de pyrolyse, pouvant fonctionner à la fois en atmosphère non oxydante et en atmosphère oxydante, coopérant avec un dispositif de mesures de la quantité de composés hydrocarbonés libérés au cours de la pyrolyse, et un dispositif de mesure du monoxyde de carbone et du dioxyde de carbone.

Le procédé selon l’invention comprend au moins les étapes suivantes :

1- Séquence de chauffe sous atmosphère inerte (pyrolyse)

2- Séquence de chauffe sous atmosphère oxydante (oxydation)

3- Caractérisation et quantification du carbone présent dans l'échantillon

Les étapes du procédé selon l’invention sont détaillées ci-après de manière non limitative pour un échantillon de sol. Les étapes du procédé selon l’invention peuvent en effet tout aussi bien être appliquées à un échantillon provenant d'une autre couche d'une formation superficielle.

1. Séquence de chauffe sous atmosphère inerte (pyrolyse)

Au cours de cette étape, un échantillon de sol est chauffé sous atmosphère inerte (comme par exemple sous un flux d’azote, d’hélium) selon une séquence de températures prédéfinies, variables dans le temps.

Selon l’invention, la séquence de chauffe sous atmosphère inerte comprend au moins une succession de paliers isothermes, chaque palier isotherme ayant une durée prédéterminée, et deux paliers isothermes consécutifs dans cette succession de paliers isothermes étant reliés entre eux par une rampe se présentant sous la forme d'un gradient thermique prédéterminé. Autrement dit, la séquence de chauffe sous atmosphère inerte comprend une succession de paliers isothermes reliés entre eux par une rampe, la durée de chaque palier isotherme pouvant être différente de la durée des autres paliers isothermes, et la pente d'une rampe pouvant être différente de la pente des autres rampes thermiques.

Selon l'invention, la succession de paliers isothermes comprend au moins six paliers isothermes reliés entre eux par cinq rampes. Selon l'invention, la succession de paliers isothermes comprend un premier palier isotherme à une température (notée TO par la suite) comprise entre 80 et 200°C, un deuxième palier isotherme à une température (notée T1 par la suite) comprise entre 340 et 380°C, un troisième palier isotherme à une température (notée T2 par la suite) comprise entre 400 et 440°C, un quatrième palier isotherme à une température (notée T3 par la suite) comprise entre 450 et 490°C, un cinquième palier isotherme à une température (notée T4 par la suite) comprise entre 500 et 540°C, un sixième palier isotherme à une température (notée T5 par la suite) comprise entre 580 et 650°C.

La figure 1 illustre de manière schématique une telle séquence de températures, en représentant l'évolution de la température T en fonction du temps t, la séquence de températures présentant six paliers isothermes aux températures TO, T1 , T2, T3, T4 et T5 telles que définies ci-dessus.

Selon l'invention, la durée prédéterminée des paliers isothermes est non nulle (par exemple supérieure à une demi-minute), et peut être préférentiellement comprise entre 3 et 5 minutes. De telles durées permettent de considérer que le craquage des composés ayant une température de craquage proche de la température du palier isotherme est terminé. De manière générale, la durée d’un palier isotherme peut être différente de la durée des autres paliers isothermes dans la séquence de chauffe en atmosphère inerte.

Avantageusement, le gradient thermique entre deux paliers isothermes peut être compris entre 1 et 50°C.min ¹ , de préférence entre 20° et 25°C.min ¹ . De telles valeurs constituent des compromis permettant de débuter le craquage thermique des composés associés au palier isotherme ayant la température la plus élevée parmi ces deux paliers isothermes, tout en limitant la durée de mise en oeuvre du procédé. De manière générale, la valeur d’un gradient thermique entre deux paliers isothermes peut être différente des autres gradients thermiques dans la séquence de chauffe en atmosphère inerte.

Selon l’invention, on mesure également, en continu, une quantité de composés hydrocarbonés libérés durant la chauffe sous atmosphère inerte, ainsi que la quantité de CO2 et de CO contenus dans l’effluent résultant de ladite chauffe. Autrement dit, au cours de cette séquence, on mesure en continu la quantité de HC, de CO et de CO2 libérés par l'échantillon par craquage thermique de la matière organique et par la décomposition thermique des minéraux carbonatés. Ainsi, on obtient à l’issue de cette étape appliquée à un échantillon donné, une première courbe représentative de la quantité de composés hydrocarbonés libérés au cours du temps pendant au moins une partie de la phase de pyrolyse, ainsi que deux autres courbes représentatives de la quantité de CO et CO2 libérés au cours du temps, pendant la phase de pyrolyse. La mesure de la quantité de composés hydrocarbonés peut être réalisée au moyen d’un détecteur du type à ionisation de flamme (FID). La mesure de la quantité de CO et CO2 libérés peut être réalisée au moyen d’un détecteur du type infrarouge (IR). A noter que de tels capteurs mesurent un flux de HC, de CO et/ou de CO2 _, et donnent des valeurs mesurées en millivolts (mV). De manière classique, on peut déterminer une quantité de HC, de CO et/ou de CC>2en déterminant une aire sous la courbe mesurée (éventuellement entre des températures prédéfinies) par ces capteurs, et en divisant cette aire par la masse en mg de l’échantillon. En variante, d’autres moyens de mesure de la quantité de HC, de CO et/ou de CO2 peuvent être utilisés.

De manière générale, cette séquence particulière de chauffe sous atmosphère inerte permet de séparer thermiquement des classes de composés caractérisés par leur température de craquage propre. Ainsi, les composés craquant entre TO et T1 correspondent à des « composés très labiles thermiquement » qui sont particulièrement abondants dans les tissus biologiques frais ; les composés craquant entre T1 et T2 correspondent aux « composés labiles thermiquement » , qui sont majoritaires dans les échantillons organiques comme les litières ou les tourbes ; les composés craquant entre T2 et T3 correspondent aux « composés résistant thermiquement » et majoritaires dans les échantillons organo-minéraux (sols) ou minéraux (alluvions, colluvions) ; les composés craquant entre T3 et T4 correspondent aux « composés réfractaires thermiquement »; et les composés craquant entre T4 et T5 correspondent aux « composés très réfractaires thermiquement » et sont présents en plus grandes proportions dans des résidus de décomposition ou des fractions exogènes, comme les matières organiques pyrogéniques ou pétrogéniques.

Selon une mise en oeuvre de l'invention, le premier palier isotherme peut être de préférence à une température comprise entre 80°C et 150°C, de manière à permettre de récupérer les contributions des composés organiques les plus labiles présents dans un échantillon de sol. De préférence, le premier palier isotherme est à une valeur de 150°C, ce qui est une température suffisante pour récupérer les contributions des composés organiques les plus labiles présents dans la plupart des formations superficielles, et notamment les sols.

Selon une mise en oeuvre de l'invention, la température maximale T5 de la séquence de chauffe sous atmosphère inerte peut être de préférence 600°C. Une telle température permet d'éviter d'obtenir des courbes de CO et de CO2 présentant des pics incomplets en fin de pyrolyse, notamment lorsque la température maximale atteint 650°C, et en particulier sur des échantillons végétaux (litières, tourbes, composts).

Selon une mise en oeuvre de l'invention, la succession de paliers isothermes selon l'invention telle que décrite ci-dessus peut être précédée par une phase de montée en température du four de pyrolyse, qui peut être sous la forme d’un gradient thermique, par exemple compris entre 1 et 50 ^o C.min ^_1 , de préférence entre 20° et 25°C.min ^_1 , ou de tout autre forme de courbe de montée en température du four de pyrolyse. Cette phase préliminaire de montée en température du four de pyrolyse permet de mettre le four de pyrolyse à la température du premier palier isotherme de la séquence de chauffe en atmosphère inerte selon l’invention. Cette phase préliminaire peut contribuer à débuter le craquage thermique des composés dont la température de craquage est inférieure à la température du premier palier isotherme, notamment dans le cas de tissus biologiques frais.

Selon une mise en oeuvre de l'invention, la succession de paliers isothermes selon l'invention telle que décrite ci-dessus peut être suivie par une phase d’abaissement de la température du four de pyrolyse, qui peut être sous la forme d'un gradient thermique, par exemple compris entre -1 et -50 ^o C.min ^_1 , de préférence entre -20° et -25°C.min ^_1 , ou de tout autre forme de courbe de baisse en température du four de pyrolyse. Cette phase terminale d’abaissement de la température du four de pyrolyse permet, si besoin, de terminer le craquage thermique des composés associés au dernier palier isotherme de la séquence de chauffe sous atmosphère inerte selon l'invention.

Ainsi, de manière générale, cette séquence particulière de chauffe sous atmosphère inerte permet de séparer thermiquement des classes de composés définis par leur température de craquage spécifique. Ce fractionnement thermique permet notamment d'obtenir des thermogrammes présentant des pics dissociés, et ainsi d'être en mesure de quantifier et d'analyser indépendamment les flux de HC, CO et C02 libérés chaque classe de composés. Les figures 2A et 2B illustrent des thermogrammes obtenus au moyen de la séquence de chauffe sous atmosphère inerte selon l'invention, respectivement pour un échantillon de sol non carbonaté et pour un échantillon de sol carbonaté. Les figures 2C et 2D illustrent des thermogrammes obtenus au moyen d'une séquence de chauffe selon l'art antérieur, ne présentant pas de paliers isothermes, obtenus respectivement pour les échantillons de sol non carbonaté et de sol carbonaté des figures 2A et 2B. Plus précisément, sur ces figures, la courbe T représente l'évolution en fonction du temps t de la température du four de pyrolyse pendant cette étape, ou autrement dit la courbe T représente la séquence de chauffe en atmosphère inerte à laquelle sont soumis les échantillons. On peut observer sur cette figure que la séquence de chauffe en atmosphère inerte mise en oeuvre pour ces échantillons comprend une phase préliminaire de mise en température du four de pyrolyse jusqu’à une température du premier palier isotherme valant 200°C, ainsi qu’une phase terminale d’abaissement de la température du four de pyrolyse à la fin du dernier palier isotherme qui est à une température de 650°C. Par ailleurs, la courbe HC représente l'évolution de l’intensité Q (en mV) du signal d’un capteur FID mesurant la quantité de composés hydrocarbonés libérés pendant la séquence de chauffe en atmosphère inerte et la courbe C02 représente l'évolution de l’intensité Q (en mV) du signal d’un capteur IR mesurant la quantité de CO2 libéré pendant la séquence de chauffe en atmosphère inerte. On peut observer sur les figures 2A et 2B que les courbes HC et C02 présente des pics bien dissociés, contrairement aux courbes des figures 2C et 2D, ce qui démontre que la séquence de chauffe multi-paliers sous atmosphère inerte permet bien un fractionnement thermique des composés contenus dans un échantillon en classes de composés ayant leur température de craquage propre.

A noter qu'une analyse de sensibilité, menée sur un panel diversifié de formations superficielles comprenant des sols et des sédiments récents, montre que la séquence de chauffe sous atmosphère inerte selon l'invention ne modifie pas significativement les quantités de HC, de CO et de CO2 libérées au cours de la pyrolyse. Autrement dit, la quantité totale de HC, de CO et de CO2 est conservée par la séquence de chauffe sous atmosphère inerte selon l'invention. Par contre, la séquence de chauffe multi-paliers selon l'invention permet de mieux dissocier les pics relatifs à des classes de composés distinctes comme discuté ci-dessus.

Par la suite, on utilisera les paramètres définis dans le tableau 1 ci-après. Plus précisément, les paramètres S 2 _t (respectivement S3C02, et S3CO, _, ), avec i variant de 0 à 6 correspondent aux quantités de HC (respectivement CO2 et CO) libérés pendant les gammes de températures DT, définies de la manière suivante : DTo = T0- ; TO < DT1 £ T1 ; T1 < DT2 £ T2 ; T2 < DT ₃ £ T3; T3 < DT £ T4; T4 < DT ₅ £ T5 ; DT ₆ = T5+ ; et où TO- correspond au premier palier isotherme à la température TO de la séquence de chauffe en atmosphère inerte, éventuellement précédé d’une phase préliminaire de montée en température du four de pyrolyse ; et où T5+ correspond au dernier palier isotherme à la température T5 de la séquence de chauffe en atmosphère inerte, éventuellement suivi d’une phase terminale d’abaissement de la température du four de pyrolyse.

[Table 1] 2) Séquence de chauffe sous atmosphère oxydante (oxydation)

Au cours de cette deuxième étape, le résidu solide de l'échantillon obtenu à l’issue de la séquence de de chauffe sous atmosphère inerte telle que décrite à l’étape 1 ci-dessus est soumis à une oxydation selon un programme de températures prédéfinies, variables dans le temps.

Le programme de températures de la séquence de chauffe sous atmosphère oxydante selon l'invention est le suivant : à partir d’une température (T'min) comprise entre 150°C et 300°C, de préférence comprise entre 150°C et 200°C, et valant de manière préférentielle 150°C ou de manière très préférée étant égale à la température du premier palier isotherme de la séquence de chauffe en atmosphère inerte de manière à pouvoir effectuer des comparaisons, on élève la température du résidu de l’échantillon issu de l’étape 1) selon un gradient de température compris entre 20 et 40 ^o C.min ^_1 , de préférence compris entre 20° et 25°C.min ^_1 , jusqu’à une température de fin d’oxydation (T'max) comprise entre 850 et 1200°C, et valant de préférence 900°C qui est une température suffisante pour épuiser le stock de carbone minéral pour la plupart des échantillons de formations superficielles.

Selon l’invention, on mesure, en continu, une quantité représentative de CO et CO2 libérés pendant cette seconde séquence de chauffe. Selon une mise en oeuvre de l’invention, cette mesure peut être réalisée au moyen d’un détecteur du type infrarouge (IR). A noter qu’un tel capteur mesure un flux de CO et/ou de CO2, et fournit des valeurs mesurées en millivolts (mV). De manière classique, on détermine une quantité de CO et/ou de CÜ2 en déterminant une aire sous la courbe mesurée (éventuellement entre des températures prédéfinies) par ce capteur, et on divise cette aire par la masse en mg de l’échantillon. En variante, d’autres moyens de mesure de la quantité de CO et/ou de CO2 peuvent être utilisés.

De manière générale, la gamme de températures préférentielle pour la température initiale de la séquence de chauffe en atmosphère oxydante, plus basse que les températures initiales connues de l'art antérieur (en général 300°C), permet d'éviter des épisodes de combustion instantanée du résidu de l'échantillon au début du cycle d’oxydation.

Les figures 3A et 3B illustrent des thermogrammes obtenus au moyen de la séquence de chauffe sous atmosphère oxydante selon l'invention, respectivement pour l'échantillon de sol non carbonaté des figures 2A et 2C et pour l'échantillon de sol carbonaté des figures 2B et 2D. Plus précisément, sur ces figures, la courbe T représente l'évolution en fonction du temps t de la température du four d'oxydation pendant cette étape, ou autrement dit la courbe T représente la séquence de chauffe en atmosphère oxydante à laquelle sont soumis les résidus des échantillons. Par ailleurs, la courbe C02 (respectivement CO) représente l'évolution de l’intensité QC02 en mV (respectivement QCO) du signal d’un capteur IR mesurant la quantité de CO2 (respectivement CO) libéré au cours du temps par le résidu issu de l’étape 1) et soumis à la séquence de chauffe sous atmosphère oxydante T. On peut observer sur ces figures les flux de CO2 au-dessus d’une température d’environ 650°C qui sont spécifiques aux sols carbonatés et témoignent de la présence de formes minérales de carbone.

Par la suite, on utilisera les paramètres définis dans le tableau 2 ci-après. Plus précisément, les paramètres S4C02, et S4CO, _, , avec i variant de 0 et 4 correspondent aux quantités de CÛ2et CO libérés pendant les gammes de températures DT', définies de la manière suivante : DT'O £ T'min ; T'min < DT'1 £ T'intl ; T'intl < DT'2 £ T'int2 ; T'int2 < DT'3 £ T'int3 et T’int3 <DT’4 £ T’max, où T'intl , T'int2, et T'int'3 sont des températures intermédiaires comprises entre T’min et T’max et sont telles que T’intl est comprise entre 420 et 480°C et vaut de préférence 450°C, T’int2 est comprise entre 520 et 580°C et vaut de préférence 550°C, et T’int3 est comprise entre 620 et 680°C et vaut de préférence 650°C. Ces températures intermédiaires correspondent à des minimas relatifs observés par la demanderesse sur des courbes de CO et de CO2 mesurées pour une pluralité d’échantillons de formations superficielles, et notamment de sols, de nature et d’origine variées. Considérant que, selon la littérature, la température de combustion est une approximation de la stabilité thermique, les flux de CO et C02 mesurés dans les gammes de température DT’1 , DT’2 et DT’3 peuvent être rapportés à des classes de constituants organiques de stabilité croissante. Considérant que la limite de stabilité thermique de la calcite est proche de T3, les flux de C02 mesurés dans la gamme de température DT’4 peuvent être rapportés à des formes minérales de carbone.

[Table 2]

3) Caractérisation et quantification du carbone présent dans l'échantillon

A cours de cette étape, il s'agit de caractériser et de quantifier le carbone présent dans l'échantillon de sol, et notamment de déterminer une teneur en carbone organique et/ou une teneur en carbone minéral. Cette étape comprend au moins les deux sous-étapes détaillées ci-dessous.

3.1) Détermination d'un paramètre caractérisant la proportion en carbone minéral dans l’échantillon

Selon l'invention, cette première sous-étape vise à déterminer un paramètre caractérisant la proportion de carbone minéral par rapport au carbone total présent dans l’échantillon considéré, à partir de la quantité de CO2 mesurée à l'étape 2 décrite ci-dessus.

Plus précisément, selon l’invention, on détermine un paramètre noté SCmin, représentatif d’une proportion en carbone minéral par rapport au carbone total dans ledit échantillon à partir d’un ratio entre la quantité de CO2 libéré par le résidu de l’échantillon au-delà de la température intermédiaire T’int3 de la séquence de chauffe sous atmosphère oxydante comprise entre 620 et 680°C (définie à l’étape 2), et la totalité de la quantité de CO2 libéré par ce résidu pendant la séquence de chauffe sous atmosphère oxydante. Ainsi, le paramètre SCmin est défini par la proportion de CO2 libéré par décomposition thermique du carbone minéral contenu l’échantillon au cours de la séquence de chauffe en atmosphère oxydante par rapport à la quantité totale de CO2 émise au cours de cette séquence de chauffe en atmosphère.

Selon une mise en oeuvre de l'invention, on détermine le paramètre SCmin selon une formule du type : où les S C02i , avec i de 1 et 4, sont tels que définis à l'étape précédente. De manière générale, le terme S4C02 correspond très majoritairement au CO2 émis par décomposition des espèce minérales carbonatées et S C02 _t correspond au flux total de CÜ2au cours de la phase d’oxydation.

Selon l'invention, on définit une valeur seuil du paramètre SCmin en deçà de laquelle on peut considérer que l'échantillon de sol est dépourvu de forme minérale (c’est-à-dire que l’échantillon est non carbonaté), ou autrement dit, en deçà de laquelle on peut considérer que toute forme de carbone contenue dans cet échantillon est de nature organique.

Selon une mise en oeuvre de l'invention, la valeur seuil du paramètre Scmin, notée SCmin_seuil, selon laquelle l'échantillon de sol considéré est dépourvu de forme minérale peut être comprise entre 0.03 et 0.05, et vaut de préférence 0.04. Ces valeurs correspondent aux erreurs liées à la mise en oeuvre du procédé selon l'invention, dues à l'appareil de mesure lui-même et à la détermination des quantités de CO2 libérées dans les gammes de températures DT', définies dans le tableau 2.

Selon une mise en oeuvre de l'invention, on peut en outre caractériser le carbone d'un échantillon en définissant quatre classes d’abondance des formes minérales de carbone d'un échantillon, les classes d'abondance pouvant être définies en fonction de la valeur du paramètre SCmin selon l'invention de la manière suivante :

Si SCmin < 0.04 , l'échantillon ne comporte pas de forme minérale de carbone ;

Si 0.04 < SCmin < 0.2 : l'échantillon comporte des formes minérales de carbone sous forme de traces ;

Si 0.2 < SCmin < 0.6 : l'échantillon comporte des formes minérales de carbone ;

Si SCmin > 0.6 : l'échantillon comporte des formes minérales de carbone en abondance.

3.2) Détermination d’une teneur en carbone organique et/ou une teneur en carbone minéral

Au cours de cette étape, on détermine une teneur en carbone organique et/ou une teneur en carbone minéral à partir des quantités de HC, de CO et de CO2 mesurées lors des étapes 1) et 2) décrites ci-dessus, et en fonction de la valeur du paramètre SCmin décrit à la sous- étape 3.1) précédente.

Selon une mise en oeuvre de l'invention, on définit des nouveaux paramètres, notés TOCfs et MinCfs par la suite, spécifiques aux séquences de chauffe sous atmosphère inerte et sous atmosphère oxydante selon l'invention et dédiés aux échantillons provenant d'une formation superficielle. Les paramètres TOCfs et MinCfs selon l'invention représentent respectivement la teneur en carbone organique et la teneur en carbone minéral d'un échantillon d'une formation superficielle, contrairement aux paramètres classiques TOC et MinC définis dans l'art antérieur, qui ne sont pas adaptés à des échantillons de formations superficielles. De plus, selon l'invention, les formules pour déterminer les paramètres TOCfs et MinCfs ainsi définies sont fonction de la valeur du paramètre SCmin par rapport à la valeur seuil SCmin_seuil définie à la sous-étape précédente. Selon l’invention, on définit deux cas possibles en fonction de la valeur du paramètre SCmin : a) Premier cas : Scmin < SCmin_seuil

Selon l’invention, si la valeur du paramètre SCmin est inférieure à la valeur seuil SCmin_seuil définie à la sous-étape précédente (c’est-à-dire si l’échantillon considéré est non carbonaté), la teneur en carbone minéral est nulle et la teneur en carbone organique est fonction de la totalité des quantités de HC, de CO et de CO2 libérés pendant la séquence de chauffe sous atmosphère inerte, et des quantités de CO et de CO2 libérés pendant la séquence de chauffe sous atmosphère oxydante jusqu'à la température intermédiaire T’int3 définie précédemment.

Ainsi, contrairement à l’enseignement de l’art antérieur qui calculait un paramètre MinC non nul même dans le cas d’une absence de carbone minéral dans l’échantillon, la présente invention permet de déterminer une teneur en carbone minéral représentative de la teneur en carbone minéral réelle d’un échantillon d’une formation superficielle sans forme minérale de carbone. De plus, la teneur en carbone organique selon l’invention est déterminée en prenant en compte, contrairement à l’art antérieur, la totalité des quantités de CO et de CO2 libérés pendant les séquences de chauffe en atmosphère inerte et en atmosphère oxydante. En effet, la détermination de la teneur en carbone organique selon l’art antérieur ne considérait que les quantités de CO et CO2 libérés en deçà de la température T2 de la séquence de chauffe sous atmosphère inerte définie à l’étape 1), la moitié de la quantité de CO libéré au-delà de la température T2 de la séquence de chauffe sous atmosphère inerte définie à l’étape 1 , et les quantités en CO et CO2 libérés en deçà de la température intermédiaire T’int3 de la séquence de chauffe sous atmosphère oxydante définie à l’étape 2) de la séquence de chauffe sous atmosphère oxydante.

Selon une mise en oeuvre de l’invention, si la valeur du paramètre SCmin est inférieure à la valeur seuil SCmin_seuil définie à la sous-étape précédente, on peut en outre définir des paramètres MinCfs et TOCfs, représentant respectivement la teneur en carbone minéral et la teneur en carbone organique d’un échantillon d’une formation superficielle, selon les formules suivantes :

- MinCfs = 0 ;

- TOCîs = PCîs + RCîs , où PCfs est un paramètre représentatif de la teneur en carbone organique pyrolysé d'une formation superficielle et RCfs est un paramètre représentatif de la teneur en carbone organique résiduelle d'une formation superficielle, ces teneurs étant définies par des formules spécifiques aux séquences de chauffe sous atmosphère inerte et sous atmosphère oxydante selon l'invention et pouvant s'écrire de la manière suivante : b) Deuxième cas : Scmin > SCmin_seuil

Selon l’invention, si la valeur du paramètre SCmin est supérieure à la valeur seuil SCmin_seuil définie à la sous-étape précédente (c’est-à-dire si l’échantillon considéré est carbonaté), on estime préalablement une quantité de CO2 libéré par craquage thermique de la matière organique contenue dans l’échantillon pendant la séquence de chauffe sous atmosphère inerte au-delà de la température T2 définie à l’étape 1) ci-dessus.

Selon l’invention, cette quantité de CO2 libéré par craquage thermique de la matière organique pendant la phase de pyrolyse au-delà de la température T2, notée S3C02 _c par la suite, est estimée à partir de la quantité de CO2 libéré pendant la séquence de chauffe sous atmosphère inerte jusqu'à la température T2, notée S3C02 _a par la suite, multipliée par un facteur k prédéfini, ce qui peut s’écrire selon une formule du type :

S3C02 _c = k x S3C02 _a

Selon une première mise en oeuvre de l’invention, la valeur du facteur k prédéfini est compris entre 1.3 et 1 .4, et vaut de préférence 1 .3724.

Selon une deuxième mise en oeuvre de l’invention, la valeur du facteur k prédéfini peut être déterminée de manière expérimentale, au moyen d’une pluralité (au moins 10, de préférence au moins 50) d’échantillons provenant de formations superficielles non carbonatées, et pour lesquels on applique au moins la séquence de chauffe sous atmosphère inerte selon l’invention à chacun des échantillons de la pluralité d’échantillons, on mesure en continu une quantité de CO2 libéré par chacun des échantillons pendant cette séquence de chauffe sous atmosphère inerte, et on détermine le facteur k au moyen d’une régression linéaire appliquée entre les quantités de CO2 libérés par la pluralité d’échantillons jusqu’à la température T2 de la séquence de chauffe sous atmosphère inerte définie à l’étape 1) et les quantités de CO2 libérés par la pluralité d’échantillons au-delà de cette température T2. En effet, les échantillons de formations superficielles sans forme minérale de carbone présentent une propriété remarquable : les quantités de CO2 émises au cours des séquences de chauffe en atmosphère inerte au-dessus et en-dessous de la température limite T2 sont linéairement corrélées. Autrement dit, la quantité totale de CO2 émise à une température supérieure à la température limite T2 vaut k fois la quantité totale de CO2 émise à une température inférieure à la température limite T3, ce qui peut s’exprimer selon une formule du type : OÙ k est établi pour des échantillons présentant un SCmin £ SCmin_seuil. Selon cette mise en oeuvre de l’invention, la valeur de ce facteur k, prédéfinie à partir d’une pluralité d’échantillons sans forme minérale, peut être généralisée au cas d’échantillons avec des formes minérales de carbone (c’est-à-dire si SCmin est supérieure à la valeur seuil SCmin_seuil ) . En effet, considérant que la source du carbone organique est en première approximation comparable pour les sols non carbonatés et les sols carbonatés, il semble raisonnable de considérer que la répartition des composés labiles par rapport aux composés thermiquement stables au-dessus des températures seuils sont comparables pour les deux types d’échantillons.

Dans le cas de d’un échantillon présentant des formes minérales de carbone, la demanderesse a montré que les quantités totales de CO2 libérées au-dessus de la température T2 de la phase de pyrolyse (c'est-à-dire le paramètre S3CO _b ) cumulent les flux issus du craquage thermique de la matière organique (c'est-à-dire le paramètre S3CO _c ) et ceux issus de la décomposition thermique de ces formes minérales. Par ailleurs, plus les teneurs en carbone minéral sont élevées, plus la contribution des formes organiques thermiquement stables est négligeable.

Selon l’invention, si la valeur du paramètre SCmin est supérieure à la valeur seuil SCmin_seuil définie à la sous-étape précédente, on détermine une teneur en carbone organique et/ou en carbone minéral au moyen des quantités de HC, de CO, et de C02 mesurées pendant la séquence de chauffe en atmosphère inerte et des quantités de CO et de C02 mesurées pendant la séquence de chauffe en atmosphère oxydante, et en prenant en compte la quantité S3C02 _c de CO2 libéré par craquage thermique de la matière organique de l’échantillon pendant la séquence de chauffe sous atmosphère inerte au-delà de la température T2, préalablement estimée tel que décrit ci-dessus.

Selon une mise en oeuvre de l’invention, si la valeur du paramètre SCmin est supérieure à la valeur seuil SCmin_seuil définie à la sous-étape précédente, la teneur en carbone organique peut être déterminée à partir de la totalité de la quantité de HC libérés pendant la séquence de chauffe sous atmosphère inerte, de la quantité de CO libéré pendant la séquence de chauffe sous atmosphère inerte jusqu'à la température T4 à laquelle on ajoute la moitié de la quantité de CO libéré pendant la séquence de chauffe sous atmosphère inerte au-delà cette température T4, de la quantité de CO2 libéré pendant la séquence de chauffe sous atmosphère inerte jusqu'à la température T2 à laquelle on ajoute la quantité estimée de CO2 libéré par craquage thermique de la matière organique pendant la séquence de chauffe sous atmosphère inerte au-delà de la température T2, et des quantités de CO et de CO2 libérés pendant la séquence de chauffe sous atmosphère oxydante jusqu'à la température intermédiaire T’int3 prédéfinie à l’étape 2).

Selon une mise en oeuvre de l’invention, si la valeur du paramètre SCmin est supérieure à la valeur seuil SCmin_seuil définie à la sous-étape précédente, on peut définir un paramètre TOCfs représentant la teneur en carbone organique d’un échantillon d’une formation superficielle selon une formule du type :

TOCîs = PCfs + RCîs , où PCfs est un paramètre représentatif de la teneur en carbone organique pyrolysé d'un échantillon d'une formation superficielle et RCfs est un paramètre représentatif de la teneur en carbone organique résiduelle d'une formation superficielle, ces teneurs étant définies par des formules spécifiques aux séquences de chauffe sous atmosphère inerte et sous atmosphère oxydante selon l'invention, et pouvant s'écrire de la manière suivante :

Selon une mise en oeuvre de l’invention, si la valeur du paramètre SCmin est supérieure à la valeur seuil SCmin_seuil définie à la sous-étape précédente, la teneur en carbone minéral peut être déterminée à partir de de la quantité de CO2 libéré pendant la séquence de chauffe sous atmosphère inerte au-delà de la température T2 à laquelle on retire la quantité estimée de CO2 libéré par craquage thermique de la matière organique pendant la séquence de chauffe sous atmosphère inerte au-delà de la température T2, de la moitié de la quantité de CO libéré pendant la séquence de chauffe sous atmosphère inerte au-delà de la température T4, et de la quantité de CO2 libéré pendant la séquence de chauffe sous atmosphère oxydante au-delà de la température intermédiaire T’int3. Selon une mise en œuvre de l’invention, si la valeur du paramètre SCmin est supérieure à la valeur seuil SCmin_seuil définie à la sous-étape précédente, on peut définir le paramètre MinCfs représentant la teneur en carbone minéral d’un échantillon d’une formation superficielle selon une formule du type :

Ainsi, les définitions des teneurs en carbone minéral et en carbone organique selon l’invention tiennent compte des formes de carbone présentes de la manière suivante : (i) toutes les formes de carbone sont organiques dans les échantillons tels que SCmin £ SCmin_seuil et (ii) les autres classes d'abondance de carbone minéral comprennent une part de carbone organique thermiquement stable à laquelle s’ajoute une part de carbone minéral, variable en nature comme en proportion.

Ainsi, à l'issue de ces deux sous-étapes, on obtient une caractérisation et une quantification fiable du carbone présent dans un échantillon d'une formation superficielle, via la détermination d’un paramètre SCmin représentant une proportion en carbone minéral dans l’échantillon, et la détermination de la teneur en carbone organique et/ou en carbone minéral de l’échantillon, par exemple au moyen de paramètres TOCfs et MinCfs définis spécifiquement pour des formation superficielles carbonatées et non carbonatées.

3.3) Détermination du statut thermique des formes organiques de carbone

Au cours de cette sous-étape, qui est optionnelle, il s'agit de déterminer le statut thermique des formes organiques de carbone présentes dans l'échantillon. Par statut thermique des formes organiques de carbone présentes dans l'échantillon, on entend un indicateur de la répartition des différentes classes de composés présentes dans l’échantillon et définies chacune par leur température de craquage.

La matière organique d'un sol est un mélange hétérogène complexe comprenant des constituants de nature et d’origine diverses : résidus de la décomposition graduelle des constituants biogéniques les plus labiles, particules libres et constituants impliqués dans les complexes organo-minéraux, constituants pyrogéniques ou pétrogéniques. Le procédé de fractionnement thermique selon l'invention ne permet pas de séparer ces constituants spécifiques. En revanche, il permet de mesurer les contributions des classes de composés définies chacune par leur température de craquage. La stabilité thermique étant considérée comme une variable liée à la stabilité biogéochimique (i.e. la résistance à la décomposition par les microorganismes), les contributions des classes de composés ainsi définies peuvent être utilisées pour décrire l’hétérogénéité de la matière organique des sols.

Selon une première réalisation de cette variante de l'invention comprenant une sous-étape de détermination du statut thermique des formes de carbone de l'échantillon, on peut définir un indice de décomposition ID et un indice de stabilité IS à partir des contributions des classes de composés définies aux étapes 1 et 2 précédentes selon des formules du type :

- Indice de décomposition ID = log [S2 ₁ + 52 ₂ ] / 52 ₃

- Indice de stabilité IS = [S2 ₃ + S2 ₄ + S2 ₅ + 52 ₆ ] / 100

Ainsi, ces deux indices sont directement liés à la fraction la plus réactive du carbone organique (c'est-à-dire la part de carbone pyrolysé sous forme de composés hydrocarbonés).

L’indice de stabilité IS mesure les contributions relatives des classes de composés thermiquement stables (52 ₃ , S2 ₄ , S2 ₅ et S2 ₆ ) qui sont particulièrement abondantes dans les formations superficielles et les couches profondes des sols, par opposition aux classes de composés plus labiles (S2 _X et S2 ₂ ) et plus abondantes dans les tissus végétaux peu décomposés présents dans les couches organiques et les couches superficielles des sols. L’indice de décomposition ID mesure le rapport entre ces classes de composés les plus labiles (52 _x et 52 ₂ ) et la classe de composés intermédiaire (52 ₃ ) qui est particulièrement abondante dans les couches organiques et organo-minérales des sols. De manière générale, l'indice de décomposition ID mesure le degré de transformation de la matière organique au fur et à mesure que les composés des classes de composés les plus labiles sont décomposés et que les composés des classes de composés les plus stables s’accumulent.

Selon une deuxième réalisation de cette variante de l'invention comprenant une sous-étape de détermination du statut thermique des formes de carbone de l'échantillon, on peut déterminer le statut thermique des formes de carbone d'un échantillon à partir des contributions relatives des différentes classes de composés d'un échantillon pour calculer une teneur en carbone organique thermiquement labile (c'est-à-dire résultant du craquage thermique et de la combustion en-dessous de T3) et une teneur en carbone organique thermiquement stable (c'est-à-dire résultant du craquage thermique et de la combustion en- dessous de T3).

Selon une mise en oeuvre de l'invention, on peut déterminer au moins un des paramètres caractérisant le statut thermique de l'échantillon suivants :

- une teneur en carbone organique pyrolysé thermiquement labile, notée COPL, définie selon une formule du type : avec KCOPL compris entre 1 et 2, et valant préférentiellement 1.3359.

- une teneur en carbone organique pyrolysé thermiquement stable, notée COPS, définie selon une formule du type : avec KCOPS compris entre 0.1 et 1 , et valant préférentiellement 0.6274.

- une teneur en carbone organique résiduel thermiquement labile, notée CORL, définie selon une formule du type :

- une teneur en carbone organique résiduel thermiquement stable, notée CORS, définie selon une formule du type :

Ces teneurs partielles sont particulièrement utiles pour comprendre la dynamique de la matière organique dans les sols car elles permettent de comparer les stocks totaux de carbone organique avec les contributions des différentes formes de carbone, et d’effectuer un suivi de l’évolution de ces stocks au cours du temps. Elles permettent notamment d’explorer les relations entre les stocks en carbone organique thermiquement labile et en en carbone organique thermiquement stable qui sont directement corrélés dans les sols, ce qui indique que le transfert d’un stock à l’autre s’opère effectivement dans les sols.

Exemples

Les caractéristiques et avantages du procédé selon l’invention apparaîtront plus clairement à la lecture de l’exemple d'application ci-après.

Le procédé selon l'invention a été appliqué à une série de 9 échantillons de sols distincts, notés par la suite ECHSi avec i variant de 1 à 9, les échantillons ECHS2, ECHS3, ECHS8, ECHS9 provenant de sols carbonatés, et les échantillons ECHS1 , ECHS4, ECHS5, ECHS6, ECHS7 provenant de sols non carbonatés. Le procédé selon l'invention a été en outre appliqué à une série de 7 échantillons de tourbes distinctes, notés par la suite ECHTi avec i=1 à 7 par la suite, provenant d’une même séquence tourbeuse et datés de 340, 2850, 4540, 5220, 6050, 8850 et 9640 ans.

La figure 4A (respectivement la figure 4B) présente un histogramme comparant les valeurs des teneurs en carbone organique Corg (respectivement en carbone minéral Cmin) déterminées par le procédé selon l'invention et par l'art antérieur pour chacun des 9 échantillons de sols ECHS1 à ECHS9. Plus précisément, la figure 4A (respectivement la figure 4B) présente les valeurs (barres de droite en gris clair) du paramètre TOCfs (respectivement MinCf) selon l'invention et les valeurs (barres de gauche en gris foncé) du paramètre TOC (respectivement MinC) selon l'art antérieur pour chacun de ces 9 échantillons. De manière générale, les valeurs de TOCfs sont globalement plus élevées (environ de 8%) que les valeurs de TOC selon l’art antérieur, ce qui correspond à la proportion de carbone organique considéré à tort comme minéral dans l’art antérieur. On peut aussi observer que les valeurs de MinCfs sont nettement plus basses que les valeurs MinC pour les échantillons de sols non carbonatés pour lesquels la valeur de MinC selon l’art antérieur correspond à des flux de CO et de CO2 émis au cours de la pyrolyse et de l’oxydation de formes organiques de carbone. On peut également observer que les valeurs de MinCfs sont plus proches, voire identiques à celle de MinC pour les échantillons de sols carbonatés pour lesquels la proportion de carbone organique comptabilisé à tort dans le MinC selon l’art antérieur est moins importante que dans les sols non carbonatés, voire négligeable si les formes de carbone organique sont entièrement pyrolysées et oxydées en dessous des températures-seuils utilisées pour calculer MINC selon l’art antérieur.

La figure 5A (respectivement la figure 5B) présente un histogramme comparant les valeurs des teneurs en carbone organique Corg (respectivement en carbone minéral Cmin) déterminées par le procédé selon l'invention et par l'art antérieur pour chacun des 7 échantillons de tourbes ECHT1 à ECHT7. Plus précisément, la figure 5A (respectivement la figure 5B) présente les valeurs (barres de droite en gris clair) du paramètre TOCfs (respectivement MinCf) selon l'invention et les valeurs (barres de gauche en gris foncé) du paramètre TOC (respectivement MinC) selon l'art antérieur pour chacun de ces 7 échantillons. On peut observer sur cette figure que les valeurs de TOC (respectivement MinC) et de TOCfs (respectivement MinCfs) sont nettement plus élevées (respectivement basses) puisque, pour ces échantillons organiques, la quasi-totalité du carbone pris en compte pour le calcul du MINC selon l’art antérieur provient en fait de la pyrolyse et de l’oxydation de composés organiques. Ainsi, la présente invention permet la caractérisation et la quantification fiable des formes de carbone présentes dans un échantillon d’une formation superficielle, que la formation superficielle soit carbonatée ou non. En particulier, la présente invention définit des séquences de chauffe en atmosphère inerte et en atmosphère oxydante appropriées aux formations superficielles, ainsi que des teneurs en carbone organique et/ou en carbone minéral prenant en compte une quantité estimée de CO2 libéré par craquage thermique de la matière organique de l’échantillon pendant la séquence de chauffe en atmosphère inerte au- delà de la température T2.