Brevet FR3156204A1 — publié le 20250606

L’invention concerne un procédé pour caractériser le carbone d’un échantillon d’une formation superficielle. Le procédé comprend une chauffe sous atmosphère oxydante de l’échantillon considéré entre une température initiale comprise entre 100 et 300 °C et une température finale comprise 650 et 1000°C, et au cours de laquelle on mesure en continu des quantités de CO et de CO 2 libérés. Puis à partir de ces quantités et d’une température limite de la séquence de chauffe comprise entre 540 et 600°C valant préférentiellement 580°C, on caractérise le carbone organique présent dans ledit échantillon en déterminant un indicateur représentatif de la stabilité thermique du carbone organique et/ou un paramètre représentatif de la labilité thermique du carbone organique. Figure 3 à publier. Procédé pour la caractérisation du carbone organique dans un échantillon d’une formation superficielle La présente invention concerne le domaine de la science du sol et des géosciences de l’environnement, plus précisément la présente invention concerne le domaine de la caractérisation du carbone contenu dans les formations superficielles, et notamment dans les sols. La présente invention trouve une application particulière dans le cas de sols à faible et très faible teneur en carbone organique. 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. Les sols à faibles teneurs en carbone organique méritent une attention particulière, car leur réhabilitation et leur restauration permettraient le stockage d’importantes quantités de carbone organique et ainsi de contribuer aux efforts d'atténuation des émissions anthropiques de gaz à effet de serre, mais aussi à l'amélioration de la qualité et de la santé des sols. Toutefois, quantifier et caractériser les faibles teneurs en carbone organique (noté Corg par la suite) des sols en présence de carbone minéral (noté Cmin par la suite) est problématique et représente un véritable défi technique. La quantification du carbone organique pour les sols dépourvus de carbone minéral est assez simple avec les outils d’analyse standards, comme un analyseur élémentaire. En revanche, elle est plus complexe dans le cas de sols en présence de Cmin. En effet, dans ce cas de figure, la détermination du Corg total nécessite des prétraitements des échantillons (calcimétrie ou décarbonatation). Pour les méthodes analytiques nécessitant de passer par une étape de calcimétrie (mesure du volume de CO2 dégagé par action de l'acide chlorhydrique (HCl) sur les minéraux carbonatés d'un échantillon, par exemple selon la méthode normalisée NF ISO 10693), la quantification des formes de carbone consiste d’une part à doser le carbone total (noté Ctotal par la suite) par combustion sèche et d’autre à déterminer le Cmin total par calcimétrie. Le Corg est alors obtenu la différence entre le Ctotal et le Cmin total. Cette méthode présente cependant des limites, car elle n’est pas adaptée aux échantillons de sols riches en Cmin et à faibles teneurs en Corg. Ces limitent sont liées aux incertitudes inhérentes aux valeurs élevées de Ctotal et de Cmin total qui se cumulent et conduisent à de mauvaises mesures de Corg, voire à des teneurs négatives. Pour les méthodes analytiques nécessitant de passer par une étape de décarbonatation (analyse consistant à mettre en contact un échantillon de sol avec un volume de HCl dilué, afin d’éliminer les différentes formes de carbonates présentes, en libérant du dioxyde de carbone), la quantification des formes de carbone consiste d’une part à doser le carbone total par combustion sèche et d’autre part à déterminer le Corg total après décarbonatation de l’échantillon. Le Cmin total est alors obtenu par la différence entre le Ctotal et le Corg total. Ce type de prétraitement conduit souvent à altérer la fraction organique la plus labile et, en conséquence, à sous-estimer les teneurs en Corg et sur-estimer les teneurs en Cmin. Les documents suivants seront cités au cours de la description : Behar, F., Beaumont, V., Penteado, H.D.B., 2001. Rock-Eval 6 technology: performances and developments. Oil Gas Sci. Technol. 56, 111–134. https://doi.org/10.2516/ogst:2001013. Disnar, J.-R., Guillet, B., Kéravis, D., Di-Giovanni, C., Sebag, D., 2003. Soil organic matter (SOM) characterization by Rock-Eval pyrolysis: scope and limitations. Org. Geochem. 34, 327–343. https://doi.org/10.1016/S0146-6380(02)00239-5. Plante, Alain F.; Fernández, José M.; Leifeld, J. (2009) Application of thermal analysis techniques in soil science. In : Geoderma, vol. 153, n° 1-2, p. 1–10. DOI: 10.1016/j.geoderma.2009.08.016. Saenger, A., Cécillon, L., Sebag, D., Brun, J.-J., 2013. Soil organic carbon quantity, chemistry and thermal stability in a mountainous landscape: A Rock–Eval pyrolysis survey. Org. Geochem. 54, 101–114. https://doi.org/10.1016/j.orggeochem.2012.10.008 Sebag, D., Disnar, J.R., Guillet, B., Giovanni, C.D., Verrecchia, E.P., Durand, A., 2006. Monitoring organic matter dynamics in soil profiles by ‘Rock-Eval pyrolysis’: bulk characterization and quantification of degradation. Eur. J. Soil Sci. 57, 344–355. https://doi.org/10.1111/j.1365-2389.2005.00745. Sebag, D., Verrecchia, E.P., Cécillon, L., Adatte, T., Albrecht, R., Aubert, M., Bureau, F., Cailleau, G., Copard, Y., Decaens, 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. https://doi.org/10.1016/j.geoderma.2016.08.025. Concernant plus particulièrement la caractérisation du carbone organique, on connait les documents (Disnar et al., 2003 ; Sebag et al., 2006 ; Saenger et al., 2013 ; Plante et al., 2009) qui décrivent des indicateurs de stabilité thermique déterminés à partir d’une analyse thermique de la matière organique des sols. Plus précisément, ces documents décrivent la mise en œuvre l’analyse thermique "ROCK-EVAL® BULK ROCK", développée initialement pour des échantillons de roches mères (Behar et al., 2001), et qui comprend des mesures de quantités de composés hydrocarbonés (HC), de monoxyde de carbone (CO) et/ou de dioxyde de carbone (CO 2 ) libérés au cours du temps par un échantillon soumis à une séquence de chauffe en atmosphère inerte puis à une séquence de chauffe en atmosphère oxydante. Selon cette méthode, 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. Puis, la température de pyrolyse est augmentée progressivement jusqu’à une température T2, généralement de 650°C. 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 CO 2 sont mesurées et représentées également sous la forme de courbes. La courbe de mesure du CO (respectivement CO 2 ) présentent deux pics, un premier pic classiquement noté S3CO (respectivement S3CO2) et qui est considéré comme correspondant au CO (respectivement CO 2 ) généré par le craquage de la matière organique de l’échantillon pendant la chauffe sous atmosphère inerte, et un deuxième pic classiquement noté S3’CO (respectivement S3'CO2) et qui est considéré comme correspondant au CO (respectivement CO 2 ) 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 CO 2 libérés 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 S4CO2) qui est considéré comme correspondant à la quantité de CO (respectivement CO 2 ) généré par la combustion de la matière organique pendant le cycle d’oxydation. Toutefois, le procédé décrit dans les documents (Disnar et al., 2003 ; Sebag et al., 2006) exploite uniquement le thermogramme S2 correspondant à aux HC libérés pendant la phase de pyrolyse. Plus précisément, le procédé décrit dans ces documents repose sur une déconvolution pour subdiviser le pic S2 (aussi noté HCpyr) en aires reflétant différents degrés de stabilité thermique de la matière organique des sols. Ces documents aboutissent à la distinction entre les « biopolymères » qui sont les composés organiques pyrolysés à moins de 400 °C et dominants dans les échantillons biologiques frais, et les « géopolymères » qui correspondent à diverses classes de composés pyrolysés au-dessus de 400 °C et dominants dans les échantillons non biologiques, c’est-à-dire les sols et les roches sédimentaires. On connait en outre le document (Sebag et al., 2016) qui décrit deux indicateurs également établis à partir d’une analyse thermique telle que décrite ci-dessus, comprenant une séquence de chauffe en deux phases : une phase de pyrolyse et une phase d’oxydation. Ces deux indicateurs, notés I et R, sont également calculés à partir du thermogramme S2 qui renseigne sur la quantité totale de composés hydrocarbonés libérés pendant la phase de pyrolyse. Plus précisément, les auteurs délimitent sur le thermogramme S2 cinq surfaces, notées A1, A2, A3, A4 et A5, correspondant chacune à des gammes de températures. Les surfaces A1 et A2, dont les gammes de températures sont comprises respectivement entre 200-340°C, et 340-400°C, permettent de considérer les classes thermiques de carbone organique hautement labiles (A1) ou labiles (A2). La surface A3, entre 400 et 460°C, correspond à une classe thermique de carbone organique plus résistant. Et pour les surfaces A4 et A5, dont les gammes de températures sont comprises respectivement entre 460-520°C et 520-650°C, les formes de carbone sont considérées respectivement comme réfractaires (A4) et hautement réfractaires (A5). A partir de ces surfaces, ce document définit les indicateurs I et R suivants : (1) (2) Ainsi, l’indicateur I utilise les surfaces A1 et A2 pour mettre en évidence le degré de transformation de la fraction organique immature, et est donc représentatif de la labilité du carbone organique d’un échantillon. L’indicateur R repose sur les surfaces A3, A4 et A5, et représente la proportion de composés carbonés thermiquement stables, et est donc représentatif de la stabilité thermique du carbone organique d’un échantillon. Par ailleurs, ce document définit aussi le diagramme I=f(R), qui permet d’appréhender la dynamique de la fraction la plus réactive (riche en liaison H-C) de la matière organique du sol, représentant moins de 20% de la matière organique totale. Comme cela sera démontré ci-après dans l’exemple d’application comparatif, ces indicateurs ne sont pas toujours représentatifs de toutes les formes de carbone organique. De plus, sur les sols à faibles teneurs en Corg, les thermogrammes S2 sont très bruités et peu reproductibles, ce qui impacte directement la fiabilité des indicateurs I et R calculés pour des échantillons présentant de très faibles teneurs de Corg. La présente invention vise à pallier les inconvénients de l’art antérieur. En particulier, la présente invention définit un nouveau protocole d’analyse thermique, comprenant une unique séquence de chauffe, en l’espèce sous atmosphère oxydante. Le procédé selon l’invention comprend en outre la définition de nouveaux indicateurs, adaptés au nouveau protocole, pour permettre la caractérisation du carbone organique, dans le cas d’échantillons à forte et à faible teneur en carbone organique. Ainsi, la présente invention permet une caractérisation plus rapide, plus simple et sur la totalité du carbone organique, et ce pour des échantillons caractérisés par une large gamme de teneurs en carbone organique. L'invention concerne un procédé pour caractériser le carbone organique présent dans une formation superficielle, à partir d’un échantillon représentatif de ladite formation superficielle, ledit échantillon n’ayant pas été chauffé en atmosphère inerte préalablement. Selon l’invention, le procédé comprend au moins les étapes suivantes : A) on chauffe ledit échantillon en atmosphère oxydante selon une séquence de températures prédéfinie dont une température initiale (T0’) est comprise entre 100 et 300 °C, et une température finale (TF’) est comprise entre 650 et 1000°C, ladite séquence de températures comprenant au moins un gradient thermique compris entre 1°C/min et 50°C/min, et on mesure en continu au moins une quantité de CO et une quantité de CO 2 libérés pendant ladite séquence de chauffe en atmosphère oxydante ; B) on caractérise le carbone organique présent dans ledit échantillon en déterminant, à partir desdites mesures en continu desdites quantités de CO et de CO 2 libérés pendant ladite séquence de chauffe en atmosphère oxydante, au moins un indicateur représentatif de la stabilité thermique du carbone organique défini par une formule du type : et/ou un paramètre représentatif de la labilité thermique du carbone organique défini par une formule du type : où , , , , et représentent la teneur en carbone libérée sous forme de CO 2 et de CO par ledit échantillon pendant ladite séquence de températures dans une gamme de températures comprises respectivement entre ladite température initiale (T0’) et 340°C, entre 340 et 400°C, entre 400 et 460°C, entre 460 et 520°C, et entre 520 et une température limite, ladite température limite étant comprise entre 540°C et 600°C et valant préférentiellement 580°C. Selon une mise en œuvre de l’invention, on peut déterminer ladite teneur en carbone libérée sous forme de CO 2 et de CO, i variant entre 1 et 5, selon une formule du type : où QCO2 et QCO représentent l’évolution dans le temps de ladite quantité respectivement de CO 2 et de CO libérés pendant ladite séquence de températures entre une température minimale et une température maximale de ladite gamme de températures de ladite teneur et m est la masse dudit échantillon. Selon une mise en œuvre de l’invention, ladite température initiale (T0’) peut valoir 200°C. Selon une mise en œuvre de l’invention, ladite température finale (TF’) peut valoir 850°C. Selon une mise en œuvre de l’invention, on peut déterminer en outre au moins un paramètre représentatif d’une teneur en carbone organique défini par une formule du type : où et représentent lesdites quantités respectivement de CO 2 et de CO libérés pendant ladite séquence de températures jusqu’à ladite température limite, et où S4'CO représente ladite quantité de CO libéré pendant ladite séquence de températures au-delà ladite température limite. Selon une mise en œuvre ded l’invention, si ladite température finale est supérieure ou égale à 850°C, on peut déterminer en outre au moins un paramètre représentatif d’une teneur en carbone minéral défini par une formule du type : où S5 représente ladite quantité de CO 2 libéré pendant ladite séquence de températures au-delà ladite température limite. Selon une mise en œuvre de l’invention, ladite séquence de températures peut comprendre un palier isotherme d'une durée prédéterminée à ladite température initiale (T’0) de ladite séquence de températures, ladite durée prédéterminée dudit palier isotherme à ladite température initiale (T’0) de ladite séquence de températures étant comprise entre 1 à 5 minutes, et valant préférentiellement 3 minutes. Selon une mise en œuvre de l’invention, ledit au moins un gradient thermique de ladite séquence de températures peut être compris entre 15°C/min et 40°C/min, et valoir très préférentiellement 25°C/min. Selon une mise en œuvre de l’invention, on peut déterminer ladite température limite en recherchant un minimum dans une courbe représentant une évolution de la quantité de carbone QC libéré sous forme de CO et de CO 2 pendant ladite séquence de températures selon une formule du type : où QCO2 et QCO représentent ladite évolution dans le temps de ladite quantité respectivement de CO 2 et de CO libérés pendant ladite séquence de températures. Selon une mise en œuvre de l’invention, on peut déterminer ledit minimum de ladite courbe représentant ladite évolution de ladite quantité de carbone QC libéré sous forme de CO et de CO 2 pendant ladite séquence de températures de manière graphique ou numérique. Liste des figures [Fig 1] La FIG. 1 illustre de manière schématique un exemple de mise en œuvre de la séquence de températures sous atmosphère oxydante du procédé selon l’invention. [Fig 2] La FIG. 2 illustre un exemple de thermogramme obtenu au moyen d’une mise en œuvre de la séquence de températures sous atmosphère oxydante selon l’invention. [Fig 3] [Fig 4] La FIG. 3 (respectivement la FIG. 4 ) illustre un exemple de courbe représentant l’évolution d’une quantité de carbone libérée sous forme de CO 2 lors de la chauffe sous atmosphère oxydante selon l’invention pour un échantillon à forte teneur en carbone organique (respectivement pour un échantillon à faible teneur en carbone organique). [Fig 5A] [Fig 5B] La FIG. 5 (respectivement la FIG. 5 ) présente sous la forme de carrés noirs la répartition des valeurs d’un indicateur de labilité en fonction des valeurs d’un indicateur de stabilité thermique déterminés par le procédé selon l’invention pour des échantillons à faible teneur en carbone organique (respectivement à forte teneur en carbone organique). [Fig 6] La FIG. 6 présente, pour les échantillons à faible teneur en carbone organique de la FIG. 5 , la répartition des valeurs de l’indicateur de labilité en fonction des valeurs de l’indicateur de stabilité thermique déterminés par le procédé selon l’invention sous la forme de carrés noirs, ainsi que la répartition des valeurs d’un indicateur de labilité en fonction des valeurs d’un indicateur de stabilité thermique déterminés par un procédé selon l’art antérieur sous la forme de cercles grisés. [Fig 7] La FIG. 7 présente, pour les échantillons à forte teneur en carbone organique de la FIG. 5 , la répartition des valeurs de l’indicateur de labilité en fonction des valeurs de l’indicateur de stabilité thermique déterminés par le procédé selon l’invention sous la forme de carrés noirs, ainsi que la répartition des valeurs d’un indicateur de labilité en fonction des valeurs d’un indicateur de stabilité thermique déterminés par un procédé selon l’art antérieur sous la forme de cercles grisés. L’invention concerne un procédé pour caractériser le carbone organique présent 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érites, argiles à silex, etc.) qui ont évolué in situ à 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. 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 µm. Ainsi, par la suite, on appelle « échantillon représentatif de la formation superficielle » un échantillon prélevé dans la formation superficielle étudiée, et ayant subi optionnellement une préparation comprenant un tamisage, un séchage (non sous atmosphère inerte) et un broyage. 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 d’oxydation en atmosphère oxydante, des moyens de mesure du monoxyde de carbone (CO) et du dioxyde de carbone (CO 2 ). Le procédé peut être alternativement mis en œuvre au moyen de tout four permettant une chauffe en atmosphère oxydante, coopérant avec un ou des dispositifs de mesures du monoxyde de carbone et/ou du dioxyde de carbone. Ainsi, contrairement aux procédés selon l’art antérieur, le procédé selon l’invention ne nécessite pas de séquence de chauffe sous atmosphère inerte (pyrolyse). En effet, la Demanderesse a observé, par des analyses réalisées sur de nombreux échantillons à faibles teneurs en carbone organique, que les courbes représentatives des composés hydrocarbonés, de CO et de CO 2 libérés durant la phase de pyrolyse présentent une faible reproductibilité, introduisant un bruit aléatoire et générant des quantifications et des caractérisations du carbone organique avec une possible erreur de mesure. La Demanderesse a ensuite observé qu’à partir de la forme des thermogrammes des sols, il est possible d’établir des indicateurs pour caractériser la matière organique des sols à partir des signaux d’oxydation, notamment à partir d’un signal combiné correspondant à la somme du CO et du CO 2 émis durant la phase unique d’oxydation. L’échantillon à analyser n’ayant pas été chauffé en atmosphère inerte préalablement (autrement dit, le procédé est appliqué directement à un échantillon représentatif de la formation superficielle tel que décrit ci-dessus), le procédé selon l’invention comprend au moins les étapes suivantes : Séquence de chauffe sous atmosphère oxydante Caractérisation du carbone organique présent dans l'échantillon Les étapes du procédé selon l’invention sont décrites 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. Séquence de chauffe sous atmosphère oxydante Au cours de cette étape, l’échantillon est soumis à une oxydation selon une séquence de températures prédéfinie, dont la température initiale (notée T0’ par la suite) est comprise entre 100 et 300 °C et vaut préférentiellement 200°C, et la température finale (notée TF’ par la suite) est comprise 650 et 1000°C, et vaut de préférence 850°C. Selon l’invention, la séquence de températures de cette chauffe sous atmosphère oxydante comprend au moins un gradient thermique compris entre 1 et 50°C/min, de préférence entre 15° et 40°C/min, et valant de préférence 25°C/min. Ainsi, au cours de cette étape, une chauffe sous atmosphère oxydante est appliquée à l’échantillon lui-même, et non pas au résidu de l’échantillon issu d’une chauffe en atmosphère inerte comme décrit dans l’art antérieur. Autrement dit, le procédé selon l’invention est caractérisé, par rapport à l’art antérieur, par une seule et unique séquence de chauffe (en l’espèce, sous atmosphère oxydante). Ou autrement dit encore, selon l’invention, l’échantillon n’est pas soumis à une chauffe en atmosphère inerte, préalablement à la chauffe en atmosphère oxydante. De manière générale, la gamme de températures pour la température initiale T0’ de la séquence de températures sous atmosphère oxydante permet d'éviter des épisodes de combustion instantanée du résidu de l'échantillon au début de l’oxydation. De manière générale, la gamme de températures pour la température finale TF’ de la séquence de températures sous atmosphère oxydante permet d’assurer une décomposition totale de toutes les formes de carbone organique, même celles qui sont très réfractaires. Si la température finale TF’ est choisie supérieure à 850°C, la séquence de températures sous atmosphère oxydante permet également d’assurer une décomposition totale de toutes les formes de carbone inorganique, notamment les minéraux carbonatés. De plus, comme cela sera discuté ci-après, cela permet d’obtenir un thermogramme présentant un creux, et ainsi de pouvoir déterminer la température limite du procédé selon l’invention. Selon une mise en œuvre de l’invention, la séquence de températures sous atmosphère oxydante peut en outre comprendre un palier isotherme à la température initiale T0’ d’une durée prédéterminée non nulle (par exemple supérieure à une demi-minute), préférentiellement comprise entre 1 et 3 minutes, et valant très préférentiellement 1 minutes. Un tel palier isotherme permet de libérer les composés labiles et/ou volatiles absorbés à la surface des particules constituants la matrice minérale du sol. La FIG. 1 illustre de manière schématique un exemple de mise en œuvre de la séquence de températures sous atmosphère oxydante du procédé selon l’invention, dans laquelle la température varie entre une température initiale T0’ maintenue pendant une durée prédéterminée non nulle et une température finale TF’, les deux températures étant reliées par un gradient thermique. Selon l’invention, on mesure, en continu, des quantités de CO et de CO 2 libérés pendant cette séquence de chauffe. Selon une mise en œuvre 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 CO 2 , et fournit des valeurs mesurées en millivolts (mV). De manière classique, on détermine une quantité de CO et une quantité de CO 2 en déterminant une aire sous les courbes respectivement de CO et de CO 2 mesurées (éventuellement entre des températures intermédiaires comme cela sera défini ci-après) par ce capteur, et on divise l’aire par la masse en mg de l’échantillon. En variante, d’autres moyens de mesure de la quantité de CO et/ou de CO 2 peuvent être utilisés. La FIG. 2 illustre un exemple de thermogramme obtenu au moyen de la séquence de chauffe sous atmosphère oxydante selon l’invention. Plus précisément, la courbe T’ représente l'évolution en fonction du temps t de la température T du four d'oxydation pendant cette étape, la courbe CO2 (respectivement la courbe CO) représente l'évolution de l’intensité QCO 2 en mV (respectivement QCO) du signal d’un capteur IR, à partir de laquelle on détermine la quantité de CO 2 (respectivement CO) libéré au cours du temps pendant la séquence de chauffe sous atmosphère oxydante T’. De manière classique, on appelle S4CO2 (respectivement S4CO) la portion de la courbe CO2 (respectivement CO) à gauche d’une température limite d’oxydation TLO, ici sensiblement égale à 600 °C, et qui est considérée classiquement comme correspondant au CO 2 (respectivement au CO) ayant été généré par le craquage de la matière organique de l’échantillon pendant la chauffe sous atmosphère oxydante. De manière classique, on appelle S5 (respectivement S4’CO) la portion de la courbe CO2 (respectivement CO) à droite de la température limite d’oxydation TLO. Lorsque les valeurs de la portion S4’CO sont nulles ou proches de zéro, la portion S5 de la courbe CO2 est considérée comme correspondant au CO 2 ayant été généré par la décomposition des minéraux carbonatés de l’échantillon pendant la chauffe sous atmosphère oxydante. Lorsque les valeurs de la portion S4’CO sont non nulles, la portion S5 (respectivement la portion S4’CO) de la courbe CO2 (respectivement CO) est considérée comme correspondant au CO 2 (respectivement CO) ayant été généré par combustion des composés organiques thermo-résistants présents dans l’échantillon pendant la chauffe sous atmosphère oxydante. Avantageusement, la séquence de températures de la chauffe sous atmosphère oxydante peut en outre comprendre un ou plusieurs paliers isothermes intermédiaires, à une température comprise entre les températures initiale et finale de la séquence de températures de la chauffe sous atmosphère oxydante. Selon une mise en œuvre de l’invention, la séquence de températures de la chauffe sous atmosphère oxydante peut comprendre un premier palier isotherme intermédiaire à une température comprise entre 200 et 360 °C, et valant préférentiellement 340°C, et/ou un deuxième palier isotherme intermédiaire à une température comprise entre 380 et 420 °C, et valant préférentiellement 400°C, et/ou un troisième palier intermédiaire à une température comprise entre 440 et 480 °C, et valant préférentiellement 460°C et/ou un quatrième palier intermédiaire à une température comprise entre 520 et 600 °C, et valant préférentiellement 580°C. Ces paliers isothermes permettent de mieux séparer, dans une courbe représentative de l’évolution en fonction de la température de la quantité de CO et de CO 2 libérés pendant la chauffe sous atmosphère oxydante, les différents classes thermiques de carbone organique, plus précisément respectivement, une classe correspondant au carbone organique thermiquement très labile, une classe correspondant au carbone organique thermiquement labile, une classe correspondant au carbone organique thermiquement résistant et une classe correspondant au carbone organique thermiquement réfractaire. Cela peut notamment permettre de contribuer à améliorer le résultat de l’étape 2) du procédé selon l’invention décrite ci-dessous car cela permet réaliser la combustion complète d’une classe avant de commencer la combustion de la suivante. Le ou les paliers isothermes intermédiaires peuvent être d’une durée prédéterminée non nulle (par exemple supérieure à une demi-minute), préférentiellement comprise entre 1 et 5 minutes, et valant très préférentiellement 3 minutes. De telles durées sont suffisantes pour libérer l’ensemble du carbone associé à une classe thermique de carbone donnée notamment en fin du cycle d’oxydation. Selon une mise en œuvre de l’invention, la séquence de températures de la chauffe sous atmosphère oxydante peut comprendre un nombre de gradients thermiques NG défini par NG = NII+1 où NII est le nombre de paliers isothermes intermédiaires de la séquence de températures. Ainsi, le ou les paliers isothermes intermédiaires sont reliés entre eux par des gradients thermiques, et le palier isotherme intermédiaire à la température la plus basse (respectivement la plus élevée) est relié également par un gradient thermique à la température initiale (respectivement la température finale) de la séquence de températures. Selon ce mode de réalisation, le ou les gradients thermiques peuvent être compris entre 1°C et 50°C/min, de préférence entre 15 et 40°C/min, et valoir très préférentiellement 25°C/min. 2) Caractérisation du carbone organique présent dans l'échantillon Au cours de cette étape, il s'agit de caractériser le carbone organique présent dans l'échantillon de sol, à partir des quantités de CO et de CO 2 mesurées en continu lors de l’étape 1) décrite ci-dessus. Plus précisément selon l’invention, il s’agit de déterminer au moins un indicateur représentatif de la stabilité thermique du carbone organique de l’échantillon considéré, noté par la suite, et/ou un indicateur représentatif de la labilité du carbone organique de l’échantillon considéré, noté par la suite, selon les formules suivantes : (3) (4) où , , , , et correspondent à la teneur en carbone libérée sous forme de CO 2 et de CO pendant la séquence de chauffe dans une gamme de températures comprises respectivement entre la température initiale de la séquence de chauffe (c’est-à-dire T0’) et 340°C, entre 340 et 400°C, entre 400 et 460°C, entre 460 et 520°C, et entre 520°C et une température limite prédéfinie. Selon l’invention, la température limite est comprise entre 540°C et 600°C et vaut préférentiellement 580°C. La température limite selon l’invention correspond à la température à laquelle se termine la libération du carbone organique (pour des températures inférieures à la température limite) et à laquelle débute (pour des températures supérieures à la température limite) la libération du carbone minéral contenu dans un échantillon. Ainsi, la gamme de températures utilisée pour déterminer le terme A5 décrit ci-dessus a, par rapport à l’art antérieur, été abaissée de 650°C à 600°C. En effet, la Demanderesse a pu observer, après des analyses réalisées sur de nombreux échantillons à teneurs en carbone variées, que la température limite au-delà de laquelle on libère du carbone minéral et non plus du carbone organique est de 600 °C lorsque l’échantillon subit une unique chauffe en atmosphère oxydante, contrairement au protocole selon l’art antérieur où la limite entre le CO 2 d’origine organique et le CO 2 d’origine minéral en phase d’oxydation après une phase de pyrolyse est établie à 650 °C (Behar et al. 2001). La FIG. 3 (respectivement la FIG. 4 ) illustre un exemple de courbe représentant l’évolution de la quantité de carbone QC-CO2 (en mgC g -1 ) libérée sous forme de CO 2 lors de la chauffe sous atmosphère oxydante selon l’invention pour un échantillon à forte teneur en carbone organique (respectivement pour un échantillon à faible teneur en carbone organique), ainsi qu’une alternance de zones grisées et blanches, notées A1 à A5, délimitant les portions de la courbe de CO 2 intervenant dans le calcul des termes A1 à A5 définis ci-dessus par leurs températures minimales et maximales. On peut observer aussi bien sur la FIG. 3 que sur la FIG. 4 un creux dans les courbes à une température d’environ 600°C. Cela illustre le fait que, quelle que soit la teneur en carbone de l’échantillon, la température limite au-delà de laquelle l’échantillon libère du carbone minéral et non plus du carbone organique est de 600 °C lorsque l’échantillon subit une unique chauffe en atmosphère oxydante. Selon une mise en œuvre de l’invention, on peut déterminer la température limite selon l’invention en recherchant un minimum dans une courbe représentant l’évolution de la quantité de carbone QC libéré sous forme de CO et de CO2 pendant la chauffe sous atmosphère oxydante s’exprimant selon une formule du type : (5) où QCO et QCO2 représentent l’évolution dans le temps (et donc en fonction de la température) de la quantité respectivement de CO et de CO 2 libérés (mesurées en mV) pendant la séquence de températures, et les coefficients et correspondent à la masse molaire du carbone divisée par la masse molaire respectivement du CO et du CO 2. Cette détermination peut être réalisée de manière graphique ou, de préférence, de manière numérique (par une fonction de recherche de minimum). . Il est bien clair que le minimum recherché est compris entre 540 et 600°C. . Avantageusement pour cette mise en œuvre, la température finale de la séquence de températures selon l’invention est supérieure ou égale à 850°C, de manière à ce qu’un creux puisse être bien distingué dans la courbe, et ainsi faciliter la détermination du minimum de la courbe, que ce soit de manière graphique ou numérique. Selon une mise en œuvre de l’invention, on peut déterminer la teneur en carbone libérée sous forme de CO 2 et de CO associée à un terme Ai selon l’invention, i variant entre 1 et 5, de la manière suivante : (6) où QCO et QCO2 représentent l’évolution dans le temps (et donc en fonction de la température) de la quantité respectivement de CO et de CO 2 libérés (mesurées en mV) pendant la séquence de températures entre les températures minimale et maximale associées au terme Ai selon l’invention, et m est la masse de l’échantillon (en g). Autrement dit, on détermine l’aire sous la courbe de mesure en continu du CO et sous la courbe de mesure du CO 2 entre les températures minimale et maximale associées à chacun des termes Ai selon l’invention définis aux équations (3) et (4), pondéré par la masse de l’échantillon, et des coefficients et correspondant à la masse molaire du carbone divisée par la masse molaire respectivement du CO et du CO 2. Par ailleurs, selon l’invention, on détermine au moins un indicateur représentatif de la qualité du carbone organique en tenant compte de la quantité de CO et de CO 2 mesurés en continu. En effet, même si a majeure partie (plus de 90%) du carbone organique total est émise sous forme de CO 2 , la quantité de CO émise n’est pas négligeable. La présente invention permet ainsi de prendre en compte la totalité du carbone organique. Cela signifie aussi que, par rapport à l’art antérieur, les indicateurs selon l’invention prennent en compte la totalité de la fraction organique du sol, et pas seulement les seuls composés hydrocarbonés libérés par pyrolyse de la matière organique du sol. De manière classique, les termes A1 et A2, associés à des températures comprises respectivement entre 200 et 340°C, et entre 340 et 400°C permettent de considérer des classes (ou « pools » en anglais) de composés carbonés hautement labiles (A1) ou labiles (A2). Le terme A3, associé à des températures comprises entre 400 et 460°C correspond à des classes de composés carbonés plus résistants. Et pour les termes A4 et A5, associés à des températures comprises respectivement entre 460 et 520°C, et entre 520 et une température limite comprise entre 540 et 600°C, les formes de carbone sont considérées comme réfractaires (A4) et hautement réfractaires (A5). Selon une mise en œuvre préférée de l’invention, on peut en outre déterminer au moins une teneur en carbone organique à partir des quantités de CO et de CO 2 mesurées en continu lors de l’étape 1). Selon cette mise en œuvre préférée de l’invention, on peut déterminer un paramètre représentatif de la teneur en carbone organique libéré par l’échantillon selon une formule du type : où S4CO2 et S4CO représentent les quantités respectivement de CO 2 et de CO libérés pendant la séquence de chauffe en atmosphère oxydante jusqu’à la température limite selon l’invention comprise entre 540 et 600°C, et valant préférentiellement 580 °C, et S4’CO représente la quantité de CO libéré pendant la séquence de chauffe en atmosphère oxydante au-delà de la température limite selon l’invention. Selon une autre mise en œuvre préférée de l’invention, si la température finale TF’ de la séquence de températures sous atmosphère oxydante est supérieure ou égale à 850°C, on peut en outre déterminer une teneur en carbone minéral de l'échantillon à partir des quantités de CO et de CO 2 mesurées en continu lors de l’étape 1). Selon cette mise en œuvre préférée de l’invention, on peut déterminer un paramètre représentatif de la teneur en carbone minéral libéré par l’échantillon selon une formule du type : où S5 représente la quantité de CO 2 libéré pendant la séquence de températures en atmosphère oxydante au-delà de la température limite selon l’invention. Les caractéristiques et avantages du procédé selon l’invention apparaitront plus clairement à la lecture de l’exemple d’application ci-après. Une cinquantaine d’échantillons de sols à faibles et fortes teneurs en Corg ont été analysés. Plus précisément, pour les sols à faible teneur en Corg, une vingtaine d’échantillons de sols a été sélectionnée à partir de leurs teneurs en Corg classées en classes (teneur en Corg inférieure ≤ 0,01%, inférieure à 0,02%, comprise entre 0,03 et 0,04%, et comprise entre 0,05 et 0,09%). Pour les sols à fortes teneurs en Corg, une trentaine d’échantillons de sols a été sélectionnée, dont les teneurs en Corg varient entre 0,8 et 14,5 %, avec une teneur moyenne de 3 %. Préalablement à leur analyse, les échantillons ont été broyés finement (dimension inférieure 200 µm) préalablement séchés à 40 °C jusqu’à stabilisation de leur poids et tamisés avec un tamis de 2 mm. Chaque échantillon a été analysé avec : - le procédé décrit dans le document (Sebag et al., 2016), pour déterminer les indicateurs R et I représentatifs respectivement de la stabilité thermique du carbone organique et de la labilité du carbone organique, et dont les formules sont données aux équations (1) et (2) ci-dessus. Pour rappel, le procédé décrit dans (Sebag et al., 2016) repose sur un protocole de chauffe de l’échantillon comprenant, en amont d’une phase d’oxydation, une phase de pyrolyse dont la température finale est de 650°C. Par ailleurs, la température maximale associée à la détermination du terme A5 dans l’équation (2) est définie à 650°C dans (Sebag et al., 2016). - le procédé selon l’invention, pour déterminer les indicateurs S et L représentatifs respectivement de la stabilité thermique du carbone organique et de la labilité du carbone organique, et dont les formules sont données aux équations (3) et (4) ci-dessus. La FIG. 5 présente (respectivement la FIG. 5 ) présente sous la forme de carrés noirs la répartition des valeurs de l’indicateur L en fonction des valeurs de l’indicateur S déterminés par le procédé selon l’invention pour les échantillons à faible teneur en carbone organique (respectivement à forte teneur en carbone organique). On peut observer que l’indicateur L (Corg labile du sol) varie entre 0,34 et 0,77 avec une valeur moyenne de 0,53, alors que l’indicateur S (Corg stable du sol) varie entre 0,29 et 0,67 avec une valeur moyenne de 0,45 dans les sols à faibles teneurs en Corg ( FIG. 5 ). Par ailleurs, on peut observer que l’indicateur L varie entre -0,07 et 0,69 avec une valeur moyenne de 0,36, alors que l’’indicateur S varie entre 0,30 et 0,87 avec une valeur moyenne de 0,48 dans les sols à fortes teneurs en Corg ( FIG. 5 ). La FIG. 6 présente, pour les échantillons à faible teneur carbone organique, sous la forme de carrés noirs, la répartition des valeurs de l’indicateur L en fonction des valeurs de l’indicateur S déterminés par le procédé selon l’invention, ainsi que, sous la forme de cercles grisés, la répartition des valeurs de l’indicateur I en fonction des valeurs de l’indicateur R déterminés par le procédé selon l’art antérieur. On peut observer que, pour un même échantillon, l’indicateur de labilité L selon l’invention est plus élevé que l’indicateur de labilité I selon l’art antérieur. Sur les indicateurs de stabilité du Corg, le contraire est observé, c’est-à-dire que l’indicateur R selon l’art antérieur est plus élevé que l’indicateur S selon l’invention. Ainsi, la FIG. 6 met en évidence des signatures de stabilité thermique de la matière organique des sols différentes entre la signature mesurée avec les indicateurs I et R selon l’art antérieur et celle mesurée avec les indicateurs L et S selon l’invention. Par ailleurs, les indicateurs de labilité selon l’art antérieur et selon l’invention (respectivement I et L) sont très peu liés positivement, avec un coefficient de détermination R² = 0,3, de même que les indicateurs de stabilité selon l’art antérieur et selon l’invention (respectivement R et S), avec un coefficient de détermination R² = 0,1. En d’autres termes, l’indicateur de labilité selon l’art antérieur I contribue à 30% à l’indicateur selon l’invention L, et l’indicateur selon l’art antérieur R ne contribue qu’à 10% à l’indicateur selon l’invention S. On a par ailleurs pu déterminer que l’indicateur I selon l’art antérieur est anti-corrélé avec l’indicateur S selon l’invention (r = -0,29) et l’indicateur R selon l’art antérieur est anti-corrélé avec l’indicateur L selon l’invention (r = -0,53). Ceci peut s’expliquer par la construction différente des indicateurs de labilité et de stabilité. La FIG. 7 présente, pour les échantillons à forte teneur carbone organique, sous la forme de carrés noirs, la répartition des valeurs de l’indicateur L en fonction des valeurs de l’indicateur S déterminés par le procédé selon l’invention, ainsi que, sous la forme de cercles grisés, la répartition des valeurs de l’indicateur I en fonction des valeurs de l’indicateur R déterminés par le procédé selon l’art antérieur. On peut observer que, pour un même échantillon, l’indicateur de labilité L selon l’invention est plus élevé que l’indicateur de labilité I selon l’art antérieur. Sur les indicateurs de stabilité du Corg, le contraire est observé, c’est-à-dire que l’indicateur R selon l’art antérieur est plus élevé que l’indicateur S selon l’invention. Ainsi, la FIG. 7 met en évidence des signatures de stabilité thermique de la matière organique des sols différentes entre la signature mesurée avec les indicateurs I et R selon l’art antérieur et celle mesurée avec les indicateurs L et S selon l’invention. Par ailleurs, les indicateurs de labilité selon l’art antérieur I et selon l’invention L sont très peu liés négativement, avec un coefficient de détermination R² = 0,2, de même que les indicateurs de stabilité selon l’art antérieur R et selon l’invention S, avec un coefficient de détermination R² = 0,1. Cela veut dire que plus l’indicateur de labilité selon l’invention L augmente, plus l’indicateur de labilité selon l’art antérieur I diminue, et plus l’indicateur de stabilité thermique selon l’invention S augmente, plus l’indicateur de stabilité thermique selon l’art antérieur R diminue. Par ailleurs, on a pu observer que les indicateurs I et R et les indicateurs L et S sont divergents et très peu liés. En effet, les coefficients de corrélations de Pearson sont de -0,41 entre indicateur L et indicateur I, et de -0,26 entre indicateur S et indicateur R. La divergence entre ces quatre indicateurs s’illustre aussi par des corrélations positives entre l’indicateur L et l’indicateur R (r = 0,45) et entre indicateur S et l’indicateur I (r = 0,32). Ceci peut s’expliquer par le fait que les indicateurs I et R sont obtenus à partir d’un principe analytique différent (phase de pyrolyse versus phase d’oxydation) sur une faible fraction du carbone organique alors que les indicateurs L et S sont calculés à partir d’un signal différent (pics S4CO et S4CO2 versus pic S2) obtenu par une seule phase d’oxydation et qui représente la totalité du carbone organique. Ainsi, la présente invention permet la caractérisation du carbone organique présent dans un échantillon d’une formation superficielle de manière simple, rapide et fiable. Le procédé selon l’invention est notamment plus simple, plus rapide et donc plus économique que les procédés selon l’art antérieur car il ne nécessite pas de séquence de chauffe sous pyrolyse (économie d’un four de pyrolyse ou d’un four d’oxydation permettant une pyrolyse, de détecteurs de HC, d’énergie pour la chauffe en pyrolyse notamment, durée de 40 minutes environ au lieu de 1h30 selon l’art antérieur). Le procédé selon l’invention permet par ailleurs une caractérisation plus fiable du carbone organique que selon l’art antérieur. En effet, la pertinence de ces nouveaux descripteurs réside dans leur capacité à caractériser la stabilité thermique de toute la partie organique de la matière organique du sol, et pas uniquement la partie organique liée aux quantités de composés hydrocarbonés de la matière organique du sol comme les descripteurs selon l’art antérieur. Procédé pour caractériser le carbone organique présent dans une formation superficielle, à partir d’un échantillon représentatif de ladite formation superficielle, caractérisé en ce que, ledit échantillon n’ayant pas été chauffé en atmosphère inerte préalablement, il comprend au moins les étapes suivantes : A) on chauffe ledit échantillon en atmosphère oxydante selon une séquence de températures prédéfinie dont une température initiale (T0’) est comprise entre 100 et 300 °C, et une température finale (TF’) est comprise entre 650 et 1000°C, ladite séquence de températures comprenant au moins un gradient thermique compris entre 1°C/min et 50°C/min, et on mesure en continu au moins une quantité de CO et une quantité de CO 2 libérés pendant ladite séquence de chauffe en atmosphère oxydante ; B) on caractérise le carbone organique présent dans ledit échantillon en déterminant, à partir desdites mesures en continu desdites quantités de CO et de CO 2 libérés pendant ladite séquence de chauffe en atmosphère oxydante, au moins un indicateur représentatif de la stabilité thermique du carbone organique défini par une formule du type : et/ou un paramètre représentatif de la labilité thermique du carbone organique défini par une formule du type : où , , , , et représentent la teneur en carbone libérée sous forme de CO 2 et de CO par ledit échantillon pendant ladite séquence de températures dans une gamme de températures comprises respectivement entre ladite température initiale (T0’) et 340°C, entre 340 et 400°C, entre 400 et 460°C, entre 460 et 520°C, et entre 520 et une température limite, ladite température limite étant comprise entre 540°C et 600°C et valant préférentiellement 580°C. Procédé selon la revendication 1, dans lequel on détermine ladite teneur en carbone libérée sous forme de CO 2 et de CO, i variant entre 1 et 5, selon une formule du type : où QCO2 et QCO représentent l’évolution dans le temps de ladite quantité respectivement de CO 2 et de CO libérés pendant ladite séquence de températures entre une température minimale et une température maximale de ladite gamme de températures de ladite teneur et m est la masse dudit échantillon. Procédé selon l’une des revendications précédentes, dans lequel ladite température initiale (T0’) vaut 200°C. Procédé selon l’une des revendications précédentes, dans lequel ladite température finale (TF’) vaut p 850°C. Procédé selon l’une des revendications précédentes, dans lequel, on détermine en outre au moins un paramètre représentatif d’une teneur en carbone organique défini par une formule du type : où et représentent lesdites quantités respectivement de CO 2 et de CO libérés pendant ladite séquence de températures jusqu’à ladite température limite, et où S4'CO représente ladite quantité de CO libéré pendant ladite séquence de températures au-delà ladite température limite. Procédé selon l’une des revendications précédentes, dans lequel, si ladite température finale est supérieure ou égale à 850°C, on détermine en outre au moins un paramètre représentatif d’une teneur en carbone minéral défini par une formule du type : où S5 représente ladite quantité de CO 2 libéré pendant ladite séquence de températures au-delà ladite température limite. Procédé selon l’une des revendications précédentes, dans lequel ladite séquence de températures comprend un palier isotherme d'une durée prédéterminée à ladite température initiale (T’0) de ladite séquence de températures, ladite durée prédéterminée dudit palier isotherme à ladite température initiale (T’0) de ladite séquence de températures étant comprise entre 1 à 5 minutes, et valant préférentiellement 3 minutes. Procédé selon l’une des revendications précédentes, dans lequel ledit au moins un gradient thermique de ladite séquence de températures est compris entre 15°C/min et 40°C/min, et vaut très préférentiellement 25°C/min. Procédé selon l’une des revendications précédentes, dans lequel on détermine ladite température limite en recherchant un minimum dans une courbe représentant une évolution de la quantité de carbone QC libéré sous forme de CO et de CO 2 pendant ladite séquence de températures selon une formule du type : où QCO2 et QCO représentent l’évolution dans le temps de ladite quantité respectivement de CO 2 et de CO libérés pendant ladite séquence de températures. Procédé selon la revendication 9, dans lequel on détermine ledit minimum de ladite courbe représentant ladite évolution de ladite quantité de carbone QC libéré sous forme de CO et de CO 2 pendant ladite séquence de températures de manière graphique ou numérique.