FR 3040216 A1 20170224 FR 1557848 A 20150820 1 Procédé d'estimation d'un temps de rétention dans une colonne de chromatographie Description DOMAINE TECHNIQUE Le domaine technique de cette demande est l'analyse d'espèces chimiques ou biologiques par chromatographie, et plus particulièrement la modélisation d'un signal de chromatographie. ART ANTERIEUR L'utilisation de la chromatographie est une technique très répandue pour l'analyse d'espèces chimiques en milieu liquide ou gazeux. Cette technique d'analyse est basée sur des colonnes de chromatographie, dont le fonctionnement est bien connu : une particule parcourt un canal, entre une entrée et une sortie, en étant portée par un fluide, dit fluide porteur, également désigné par le terme phase mobile. La paroi du canal comporte un revêtement, appelé phase stationnaire, avec lequel la particule présente une affinité, de telle sorte que la particule est apte à être momentanément adsorbée, puis désorbée. En fonction de l'affinité avec la phase stationnaire, le parcours de la particule à travers le canal peut être plus ou moins long. Une colonne de chromatographie comporte également un détecteur, placé à la sortie du canal, pour détecter la présence de la particule lorsque cette dernière sort du canal. Généralement, le signal détecté par le détecteur prend la forme d'un histogramme représentant le nombre de détections en fonction du temps, cet histogramme étant désigné par le terme 20 chromatogramme. Lorsque différentes particules, de différentes espèces, sont injectées en même temps dans une colonne, le temps de parcours de chaque particule dépend de son affinité avec la phase stationnaire, dépendant de la nature chimique de la particule. Aussi, le chromatogramme présente différents pics, chaque pic représentant la distribution du temps de parcours de 25 particules de la même espèce dans la colonne. Le chromatogramme peut être modélisé, en fonction de chaque espèce de particule, et cela a été largement décrit dans la littérature. Cette modélisation correspond à l'établissement d'un modèle dit modèle direct, permettant d'estimer un chromatogramme connaissant la nature des particules introduites. Cela est par exemple décrit dans la demande de brevet W0201026228, 3040216 2 cette dernière décrivant la modélisation d'un chromatogramme basé sur la prise en compte d'un modèle spatio-temporel. Le parcours de la particule à travers la colonne dépend de paramètres de fonctionnement de la colonne, ces paramètres ayant une influence sur l'adsorption ou la désorption de la particule. Il 5 peut s'agir de la température de la colonne ou de la composition des phases mobile ou stationnaire. Afin d'obtenir une meilleure séparation entre particules d'espèces chimiques différentes, ces paramètres de fonctionnement peuvent être modulés. Par exemple, la température ou la composition de la phase mobile peuvent varier temporellement entre l'entrée et la sortie de la particule. Cela peut également concerner la phase stationnaire, dont 10 la composition peut varier spatialement, si bien que durant le parcours de la particule, les propriétés d'adsorption et/ou de désorption évoluent temporellement. Du fait de cette modulation, qu'elle soit temporelle ou spatiale, les propriétés d'adsorption et de désorption de la phase stationnaire vis-à-vis de chaque type de particule sont modifiées. Or, les modélisations réalisées selon l'art antérieur sont basées sur des comportements moyens 15 de particules, établis selon des modèles déterministes. Ces approches ne sont pas adaptées à une situation dans laquelle un paramètre de fonctionnement de la colonne évolue durant le parcours de la particule dans la colonne. Les inventeurs ont déterminé un procédé de modélisation permettant d'établir un modèle direct pouvant s'appliquer à des colonnes de chromatographie dont un, voire plusieurs, paramètres de fonctionnement évolue durant le parcours d'une particule dans la colonne. EXPOSE DE L'INVENTION Un objet de l'invention est un procédé d'estimation d'un temps de rétention d'une particule dans une colonne de chromatographie, ladite colonne de chromatographie comportant un canal dont la paroi est formée d'un revêtement, dit phase stationnaire, apte à adsorber et désorber ladite particule, ce dernier étant porté par un fluide porteur, dit phase mobile, parcourant la dite colonne entre une entrée et une sortie, définissant un temps mort, le procédé comportant les étapes suivantes : a) estimation d'un temps passé dans la phase mobile, après l'entrée de la particule dans la colonne et jusqu'à une première adsorption dans la phase stationnaire, à l'aide d'une loi de probabilité du temps passé par la particule dans ladite phase mobile, dite première loi de probabilité ; 3040216 3 b) suite à chaque adsorption de la particule dans la phase stationnaire, estimation d'un temps passé par la particule dans ladite phase stationnaire jusqu'à une désorption, à l'aide d'une loi de probabilité du temps passé par la particule dans la phase stationnaire, dite deuxième loi de probabilité ; 5 c) suite à chaque désorption de la phase stationnaire, estimation d'un temps passé de la particule dans la phase mobile jusqu'à une nouvelle adsorption dans la phase stationnaire, à l'aide de ladite première loi de probabilité ; d) répétition des étapes b) et c) jusqu'à la ce que la somme des estimations du temps passé dans la phase mobile soit supérieure audit temps mort ; 10 e) sommation des estimations des temps passés par la particule dans la phase stationnaire, afin d'estimer le temps de rétention de ladite particule dans la colonne. Selon un mode de réalisation, durant le parcours de la particule à travers la colonne, un paramètre régissant le fonctionnement de la colonne varie, de telle sorte que la première loi ou la deuxième loi évolue en fonction dudit paramètre. 15 Le paramètre peut alors être choisi parmi - une température de la colonne, - une composition de la phase mobile, en particulier en phase liquide, auquel cas le paramètre évolue en fonction du temps durant le parcours de la particule, selon une pluralité de paliers. 20 Le paramètre peut être une composition de la phase stationnaire, cette composition variant entre l'entrée et la sortie de la colonne. Selon un mode de réalisation, durant le parcours de la particule, un paramètre est porté successivement à différents paliers, à chaque palier correspondant une valeur du paramètre, à chaque palier de température étant attribué un rang, l'étape b) comprenant alors les sous- 25 étapes suivantes : i) estimation de l'instant auquel a lieu ladite adsorption, et détermination de la température de la colonne à cet instant ; ii) estimation d'un temps passé par la particule dans la phase stationnaire jusqu'à une désorption, à l'aide d'une loi de probabilité du temps passé dans la phase stationnaire, 30 cette loi dépendant d'une valeur dudit paramètre audit instant d'adsorption. 3040216 4 Ce mode de réalisation peut également comprendre les étapes de : iii) estimation d'un instant de désorption de ladite particule, à partir de l'estimation établie lors de la sous-étape ii) ; iv) détermination d'une valeur du paramètre de la colonne correspondant audit instant de 5 désorption ; v) lorsque ladite valeur audit instant de désorption est différente de la valeur considérée lors de la sous-étape ii), réitération des sous-étapes ii) à iv) en considérant que l'instant d'adsorption correspond au début du palier de rang 10 suivant ; lorsque ladite valeur audit instant de désorption correspond à la valeur considérée lors de la sous-étape ii), détermination d'un temps passé dans la phase stationnaire au cours de ladite adsorption. Le temps passé dans la phase stationnaire, au cours de l'adsorption i, peut alors être obtenu en 15 fonction d'au moins un des paramètres suivants : - la durée de chaque palier entre l'instant d'adsorption, déterminé lors de l'étape i) et l'avant-dernière itération des sous-étapes ii) à v), - l'estimation du temps passé dans la phase stationnaire, déterminé lors la dernière itération des sous-étapes ii) à v). 20 Ladite première loi de probabilité peut suivre une loi exponentielle, paramétrée par un taux d'adsorption, ce dernier représentant une probabilité par unité de temps que ladite particule, se trouvant dans la phase mobile, soit adsorbée dans la phase stationnaire. Ladite loi de probabilité du temps passé dans la phase stationnaire suit une loi exponentielle, paramétrée par un taux de désorption, ce dernier représentant une probabilité par unité de 25 temps que ladite particule, se trouvant dans la phase stationnaire, soit désorbée de ladite la phase stationnaire. Ledit temps mort peut notamment correspondre au temps moyen passé par une particule, transportée par le fluide porteur, dans ladite colonne sans adsorption dans la phase stationnaire. 3040216 5 Selon un mode de réalisation, les étapes a) à e) sont réalisées pour une pluralité de particules de même espèce, de façon à obtenir une distribution d'un temps de rétention de ladite espèce dans la colonne. Selon un mode de réalisation, les étapes a) à e) sont réalisées en considérant une pluralité de 5 particules de différentes espèces de façon à obtenir une distribution d'un temps de rétention pour chacune de ces espèces. Le procédé peut alors comprendre une étape f) de détermination d'un indicateur statistique d'au moins une distribution du temps de rétention d'une espèce, par exemple une moyenne, ou un moment. 10 Un autre objet de l'invention est un support d'enregistrement d'informations, comportant des instructions pour l'exécution d'un procédé tel que décrit dans cette demande, ces instructions étant aptes à être exécutées par un processeur. L'invention repose sur une description du comportement individuel d'une molécule durant son passage à travers la colonne, selon une technique dite de sondage. Cette approche 15 microscopique permet de tenir compte d'une évolution d'un paramètre régissant le fonctionnement de la colonne durant le parcours d'une molécule. En décrivant le comportement d'un ensemble de molécules de la même espèce, comportant un nombre de molécules considéré comme suffisamment représentatif, l'invention permet d'établir une distribution du temps de rétention correspondant à ladite espèce, à partir de laquelle on peut 20 déterminer des indicateurs statistiques, en particulier la moyenne, la variance ou autres moments. L'invention sera mieux comprise par la lecture des exemples de réalisation décrits ci-après, lesquels se basent sur les figures suivantes. FIGURES 25 La figure 1 représente une colonne de chromatographie. La figure 2 représente les étapes d'un procédé selon un premier mode de réalisation de l'invention. Les figures 3A et 3B présentent respectivement les étapes d'un procédé selon un deuxième et un troisième mode de réalisation de l'invention. 3040216 6 La figure 4A représente le résultat de simulations du parcours de 500 particules d'une même espèce. Elle représente l'évolution du temps cumulé passé, par chaque molécule, dans la phase mobile, entre chaque adsorption. La figure 4B représente la distribution du temps de rétention dans la colonne pour l'ensemble des parcours modélisés. 5 La figure 5 représente des estimations du temps de rétention ajusté moyen de différentes molécules, obtenues selon un procédé selon l'invention, confrontées à un chromatogramme expérimental d'un mélange comprenant lesdites molécules. EXPOSE DE MODES DE REALISATION PARTICULIERS La figure 1 représente une colonne de chromatographie 1, comportant un canal 10 et un 10 détecteur 20. La paroi du canal 10 comporte un revêtement 12, dit phase stationnaire S. Le canal comporte également une partie centrale 14 dans lequel un fluide porteur est apte à circuler, entre une entrée in et une sortie out. Un tel fluide porteur est généralement désigné par le terme de phase mobile M. Le fluide porteur peut être un gaz ou un liquide, dont l'affinité avec la phase stationnaire S est 15 négligeable, si bien que le fluide porteur n'interagit pas avec la phase stationnaire durant son parcours dans la colonne. Son temps de parcours dans la colonne, c'est-à-dire entre l'entrée et la sortie de la colonne, est désigné par le terme temps mort, et noté to. Ce temps mort to correspond au temps de parcours d'une particule n'interagissant pas avec la phase stationnaire. La fonction du fluide porteur est de transporter une particule x entre l'entrée et la sortie du 20 canal. Durant son parcours dans la colonne, la particule x interagit avec la phase stationnaire S, ce qui entraîne une pluralité d'adsorptions / désorptions. Entre chaque adsorption, la particule est portée par la phase mobile M, et évolue à la même vitesse que cette dernière. Sur la figure 1, le trajet de la particule x, dans la colonne 1, est représenté par une succession de flèches, chaque flèche montrant le parcours de la particule x entre deux adsorptions successives. Les 25 temps passés entre chaque adsorption sont également représentés. Ils sont désignés par (entre l'entrée et la première adsorption), tim., (entre la première et la deuxième adsorption), , t3,m , t4,m , ts,in , t6,m, tuv,,, D'une façon générale, le temps tim désigne le temps passé dans la phase mobile par une particule x suite à une eme adsorption, i étant un entier positif compris entre 0 et Ni , Ni représentant un nombre d'adsorptions subies par la 30 particule x dans la colonne. 3040216 7 A la sortie du canal est placé un détecteur 20, apte à émettre un signal représentatif du nombre de particules sortant de la colonne en fonction du temps. Ce signal correspond au chromatogramme évoqué dans la description de l'art antérieur. Le détecteur 20 est apte à être connecté à un processeur 30, ce dernier étant relié à une 5 mémoire 32 comportant des instructions, ces dernières pouvant être exécutées par le processeur 30, pour mettre en oeuvre le procédé représenté sur les figures 2, 3A ou 3B, et décrit ci-après. Ces instructions peuvent être sauvegardées sur un support d'enregistrement, lisible par un processeur, de type disque dur, CDROM ou autre type de mémoire. Les inventeurs ont choisi de modéliser le chromatogramme, c'est-à-dire de déterminer un 10 modèle direct de la colonne, en suivant une approche originale, basée sur un modèle stochastique, en simulant le parcours individuel d'une particule, de façon séquentielle et aléatoire, en se basant sur des lois de probabilités décrivant le comportement de la particule dans la colonne. La simulation d'un grand nombre de particules x d'une même espèce X permet alors d'obtenir une distribution statistique relative à leur temps de parcours dans la colonne, 15 également appelé temps de rétention tR . Une telle distribution constitue le modèle direct de la colonne vis-à-vis de ladite espèce de particule. Par particule, on entend une entité chimique ou biologique apte à être détectée par le détecteur. Il s'agit par exemple d'une molécule, d'une protéine, d'un peptide d'un agrégat moléculaire, d'un complexe moléculaire, d'une nanoparticule. 20 Par particule de même espèce, on entend une particule de même nature chimique, ou constituée par les mêmes éléments chimiques. Dans la suite de la description, l'indice X désigne une espèce de particule, tandis que l'indice x désigne une particule individuelle. La modélisation est basée sur l'utilisation : d'une loi de probabilité du temps passé par la particule dans la phase mobile, avant ou 25 après chaque adsorption dans la phase stationnaire, cette loi étant désignée par la suite par le terme première loi de probabilité. d'une loi de probabilité du temps passé par la particule dans la phase stationnaire, au cours de chaque adsorption, cette loi étant désignée par la suite par le terme deuxième loi de probabilité. 30 La première et la deuxième lois de probabilités dépendent de l'espèce X de chaque particule et de paramètres P régissant le fonctionnement de la colonne de chromatographie. 3040216 8 D'une façon générale, le parcours d'une particule dans la colonne est décrit par des tirages aléatoires successifs, respectivement selon la première et la deuxième loi de probabilité, afin de déterminer des temps passés : entre chaque adsorption, dans la phase mobile, 5 au cours de chaque adsorption, dans la phase stationnaire. Comme précédemment évoqué en lien avec l'art antérieur, le fonctionnement de la colonne 1 est régi par des paramètres de fonctionnement P, dits paramètres, ces derniers ayant une influence sur l'adsorption et la désorption de la particule. De tels paramètres de fonctionnement sont par exemple : 10 la température T de la colonne ; la composition de la phase mobile, en particulier dans le cas où il s'agit d'un liquide ; la composition de la phase stationnaire ; Ces paramètres peuvent évoluer lorsque la particule traverse la colonne. La température et/ou la composition de la phase mobile peuvent être modulées temporellement, tandis que la 15 composition de la phase stationnaire peut être modulée spatialement, entre l'entrée in et la sortie out de la colonne. Dans les modes de réalisation décrits ci-dessous, la première et la deuxième loi de probabilité suivent des lois exponentielles, respectivement paramétrées par un taux d'adsorption ka xp et un taux de désorption kd xp. Le taux d'adsorption ka xp correspond à la probabilité par unité 20 de temps qu'une particule x de l'espèce X, présente dans la phase mobile, soit adsorbée dans la phase stationnaire, en considérant des paramètres de fonctionnement P. Son inverse correspond à une durée moyenne passée par la particule dans la phase mobile, entre deux adsorptions successives. Le taux de désorption kd xp correspond à la probabilité par unité de temps qu'une particule x de l'espèce X, présente dans la phase stationnaire, soit désorbée, en 25 considérant des paramètres de fonctionnement P. Son inverse 1correspond à une durée ka x,p moyenne passée par la particule dans la phase stationnaire, suite à une adsorption. Les taux d'adsorption et de désorption sont connus et peuvent être mémorisés dans la mémoire 32. Ainsi, les première et deuxième lois de probabilité s'écrivent respectivement : - fr,,(t, ka x,p) = ka X,P X e- xt kaxp (1) 30 - fd (t, kd xr, r) x e txkdx,p = kd X,P (2) 1 ka X,P 3040216 9 On comprend que l'indice P désigne soit un, soit plusieurs paramètres de fonctionnement. La figure 2 représente un premier mode de réalisation d'un procédé selon l'invention, dans lequel les paramètres de fonctionnement de la colonne sont constants durant le parcours d'une particule. La particule x est, dans cet exemple, une molécule gazeuse. 5 Étape 100: initialisation. La molécule x est introduite dans la colonne. Étape 110 : détermination du temps de parcours avant la première adsorption. Cette étape comprend un tirage aléatoire conditionné par la première loi de probabilité, de façon à établir un temps tom passé dans la phase mobile, avant une première adsorption. Le procédé suit alors un processus itératif, l'indice i désignant le rang de l'itération. 10 Étape 120 : lors de chaque adsorption i, réalisation d'un tirage conditionné par la deuxième loi de probabilité, de façon à établir un temps tis, passé par la molécule dans la phase stationnaire au cours de ladite adsorption. Étape 130: lors de chaque désorption suivant une adsorption i, réalisation d'un tirage aléatoire conditionné par la première loi de probabilité, de façon à établir un temps tim, passé dans la 15 phase mobile suite à la en' adsorption. Étape 140 : réitération ou sortie de l'algorithme. Lors de l'étape 140, l'indice i est incrémenté et les étapes 120 et 130 sont réitérées, jusqu'à ce que le cumul du temps passé dans la phase mobile atteigne le temps mort to. Cette condition de réitération correspond à E1=0 t0 (3) 20 Lorsque cette condition est atteinte, l'étape 140 détermine le temps passé par la particule x dans la colonne, dit temps de rétention, et noté tR, selon l'expression tO = tR (4) où N, représente le nombre d'itérations réalisées, c'est-à-dire le nombre d'adsorptions dans la phase stationnaire. 25 On peut également définir un temps de rétention tR dit ajusté, correspondant au temps cumulé passé par la particule x dans la phase stationnaire, soit = tR (41 On va à présent décrire un deuxième mode de réalisation, en lien avec la figure 3A, selon lequel la colonne est modulée en température. Plus précisément, durant le parcours de la molécule x 3040216 10 à travers la colonne, la température T évolue selon une fonction temporelle prédéterminée. A chaque instant t, correspond une température de colonne T telle que T = f(t), f étant une fonction décrivant l'évolution temporelle de la température. Par exemple, la colonne est portée successivement à NK paliers de température Tk, k étant un 5 entier strictement positif représentant le rang de chaque palier de température. Une modulation en température correspond à une succession de paliers de température dans un ordre prédéterminé. La fonction d'évolution de température f peut alors être exprimée selon l'expression : f (t) = Tk = TL + 77 T vIV N K i ticK1" u (t C X i\-") ( ) 10 où : Ti est une température initiale de la modulation en température, cette température correspondant par exemple à l'entrée de la particule dans la colonne ; - tf est la durée de la modulation en température de la colonne ; Tf est la température finale de la colonne, lors du dernier palier ; 15 H est la fonction de Heaviside, avec H(t) = 0 si t O. Selon ce mode de réalisation, on fait l'hypothèse que le temps mort to est indépendant de l'évolution de la température. Aussi, selon cette première approche, on considère que la première loi de probabilité ne dépend pas de la température, le coefficient d'adsorption ka étant considéré comme indépendant de la température. La deuxième loi de probabilité dépend de la 20 température, du fait de l'évolution du ratio -ka en fonction de la température, exprimé par ka l'équation de Van't Hoff : ka 1 -AH-ELS - = -e RT R kd /3 où : AH représente l'enthalpie d'activation, dont la valeur dépend du matériau composant 25 la phase stationnaire et de la température ; AS représente l'entropie d'activation dont la valeur dépend du matériau composant la phase stationnaire et de la température ; R est la constante des gaz parfaits (6) 3040216 11 - fl est le rapport des phases, c'est-à-dire le volume de la phase mobile sur le volume de la phase stationnaire. AH AS Il vient : kd = flicaer+ (6'). Les différentes étapes de ce mode de réalisation sont les suivantes : 5 Étape 200: initialisation : la température Tk est initialisée à Tk.i. Étape 210 : détermination du temps de parcours dans la phase mobile avant une première adsorption. Cette étape comprend un tirage aléatoire conditionné par la première loi de probabilité, de façon à établir un temps to,m passé dans la phase mobile, avant une première adsorption. 10 Le procédé suit alors un processus itératif, l'indice i désignant le rang de l'itération. Étape 220 : lors de chaque adsorption i, détermination d'un temps t,,,, passé par la particule dans la phase stationnaire au cours de ladite adsorption. Cette étape prend en compte le fait que durant une adsorption i, la température de la phase stationnaire est susceptible d'évoluer. Elle englobe différentes sous-étapes, effectuées de façon itératives, l'indice j désignant le rang 15 de ces itérations, cet indice étant initialisé à 1. A chaque indice j est associé un palier de température Tk. L'indice j = 1 correspond au palier de température Tk durant lequel a lieu l'adsorption i, les indices suivants correspondant aux paliers de température TkA4 , Tk+2... suivants. Sous-Étape 221 : détermination de l'instant, noté tif 1 , auquel a lieu l'adsorption de rang i, et détermination de la température Tk correspondant à cet instant, par l'application de la fonction 20 d'évolution de température Tk = f (tif-1). L'instant ti 1 auquel a lieu d'adsorption i est obtenu par t. -1 = tt,n, + tu (7) où : tl,m correspond à la somme des temps ti,n passés dans la phase mobile jusqu'à l'adsorption i ; lors de la première adsorption (i = /), ce terme est égal à tom-, . 25 tu correspond à la somme des temps tu passés dans la phase stationnaire jusqu'à l'adsorption i ; lors de la première adsorption (i = /), ce terme est nul. Sous-Étape 222 : estimation d'une durée d'adsorption dans la phase stationnaire, au cours de l'adsorption i, lors du palier de température Tk associé à l'indice d'itération j, en effectuant un tirage aléatoire conditionné par la deuxième loi de probabilité, de façon à établir un temps ti,s,k 3040216 12 passé dans la phase stationnaire, lors de l'adsorption i, au cours du palier de température Tk. Cette estimation est réalisée, par tirage aléatoire selon la deuxième loi de probabilité : lorsque j =1, à l'instant auquel la en' adsorption a lieu, ce dernier étant déterminé lors de la sous-étape 221 ; 5 lorsque j >1, à l'instant tif correspondant au début du palier de température Tk associé à l'itération j. Sous-Etape 223 : A partir du tirage de estimation d'un instant de désorption tid, associé à l'itération j, selon l'expression : ti d = ti,s,k (8) 10 Cette estimation n'est validée que si, à l'instant de désorption tifd, la température de la colonne correspond à la température Tk prise en compte pour estimer le temps ti,s,k passé dans la phase stationnaire. Autrement dit, cette estimation est validée si il n'y a pas de changement de température entre les instants tifd et C'est l'objet des étapes 224 et 226 Sous-Etape 224 : estimation de la température de la colonne, à l'instant de désorption tifd 15 déterminé lors de l'étape 223, par application de la fonction d'évolution de la température f(tifd) Sous-Etape 225 : validation de l'estimation de l'instant de désorption estimé lors de l'étape 223, en comparant la température f(tifd) à l'instant de désorption estimé ainsi que la température f(tif) à laquelle cet instant a été estimé : Si f(tifd) # f(tif), c'est-à-dire si la température à l'instant de désorption 20 ti d = ti + ti,s,k est différente de la température Tk à laquelle la durée d'adsorption t,,,,k a été estimée, l'estimation de l'instant de désorption associéee à l'itération j n'est pas validée. En effet, cela signifie que la température de la colonne a évolué entre l'instant tif, auquel l'estimation de la durée d'adsorption ti,s,k est effectuée, et l'instant de désorption tifd. 25 Dans ce cas, une nouvelle itération j +1 est effectuée, en repartant de la sous-étape 222, en considérant que l'instant d'adsorption trElcorrespond au début du palier de température Tk+i suivant le palier de température Tk correspondant à l'itération j. Autrement dit, tif+1 = Eqq:k1 Atq (9) 3040216 13 où Atq représente la durée du palier de température de rang q, et k représente le rang de l'échelon de température correspondant à l'itération courante : Tk = f(t') Les sous-étapes 221 à 225 sont alors réitérées, en prenant en compte une température de colonne T = Tk+i f(ti-1+1). 5 Ainsi, du fait du décalage temporel effectué en réactualisant le temps d'adsorption selon l'expression (9), lors de chaque nouvelle itération j, on considère que la particule vient d'être adsorbée dans une phase stationnaire, portée à une température Tk = f(t) et dans laquelle la durée d'adsorption est conditionnée par la loi de probabilité fd (t, kd x,p=T) = kd X,T x e xj- La colonne de chromatographie est alors 10 considérée comme étant la succession de NK colonnes élémentaires indépendantes, chaque colonne élémentaire étant portée à une température Tk. Si f(tid) = f(t), c'est-à-dire si la température à l'instant de désorption estimé tic/ correspond à la température à laquelle la durée d'adsorption ti,s,k a été estimée, l'estimation de l'instant de désorption est validée. On passe alors à l'étape 230, en 15 considérant que le temps passé dans la phase stationnaire S au cours de l'adsorption i est tel que ti,s = tic/ - ti d1 (10) Étape 230 : lors de chaque désorption suivant une adsorption i, réalisation d'un tirage aléatoire conditionné par la première loi de probabilité, pour estimer un temps t,,', passé dans la phase mobile entre l'adsorption de rang i et l'adsorption de rang i +1. 20 Étape 240: itération suivante ou sortie de l'algorithme. Les étapes 220 et 230 sont réitérées, en incrémentant l'indice i, jusqu'à ce que le cumul du temps passé dans la phase mobile atteigne le temps mort to, de même que dans le mode de réalisation précédent. Cette condition correspond à to 1=0 25 Lorsque cette condition est remplie, l'étape 240 détermine le temps passé par la particule dans la colonne, dit temps de rétention, et noté tR, tel que tO = tR 1=1 3040216 14 On peut également définir un temps de rétention tR dit ajusté, correspondant au temps passé par la particule dans la phase stationnaire, soit : EL t1 = tR' . Ce procédé est mis en oeuvre pour une pluralité de particules x de même espèce X, de façon à déterminer une distribution du temps de rétention tR ou du temps de rétention ajusté t;? pour 5 ladite espèce X. On peut alors estimer des paramètres statistiques relatifs à cette distribution, tels qu'une moyenne, un paramètre représentatif d'une dispersion, par exemple la variance ou l'écart-type, ou d'autres paramètres. Le procédé permet alors une estimation d'un chromatogramme correspondant à ladite espèce X sans recours à une expression analytique. Le recours à un procédé stochastique permet d'effectuer, lors de chaque adsorption, un tirage 10 selon une loi de probabilité tenant compte de l'évolution de la température de la colonne. On comprend qu'un tel tirage est effectué de façon séquentielle, ce qui permet une prise en compte des conditions de température à l'instant où le tirage est réalisé. Le procédé est donc particulièrement adapté à une évolution de température de la phase stationnaire durant le parcours de la particule à travers la colonne. 15 De façon plus générale, le procédé est adapté à toute variation temporelle ou spatiale d'un paramètre de la colonne durant le parcours d'une particule. Cette variation entraîne une évolution modélisable de la première loi et/ou de la deuxième loi de probabilité précédemment décrites au cours de ce parcours. Il est alors possible de modéliser un temps de rétention par une estimation du temps passé par la particule dans la phase stationnaire lors de chaque 20 adsorption. On comprend qu'outre une variation temporelle de la température de la colonne, le procédé est adapté à modéliser une variation temporelle de la composition C de la phase mobile M (c'est-à-dire du fluide porteur), ou encore une variation spatiale de la composition de la phase stationnaire S. Selon un troisième mode de réalisation, la composition C de la phase mobile M évolue en 25 fonction du temps, selon différents paliers. Dans ce cas, une fonction d'évolution temporelle g décrit l'évolution de cette composition en fonction du temps. Ainsi, selon ce mode de réalisation, la composition de la phase mobile évolue selon des paliers Ck tels que Ck =g(t). Les principales étapes de ce mode de réalisation sont décrites ci-dessous, en lien avec la figure 3B. L'étape 300 est une étape d'initialisation au cours de laquelle la composition de la phase mobile 30 est initialisée de telle sorte que C = Ck.0. L'étape 310 correspond à l'étape 210 préalablement décrite. A l'issue de cette étape, le procédé suit un processus itératif, analogue à celui décrit en 3040216 15 lien avec la figure 3A, l'indice i désignant le rang de l'itération. Selon ce mode de réalisation, la première loi peut varier en fonction de la composition de la phase mobile. Étape 320 : lors de chaque adsorption i, détermination du temps t,,,, passé par la particule dans la phase stationnaire. Cette étape prend en compte le fait que durant une adsorption i, la 5 composition de la phase mobile est susceptible d'évoluer. Elle englobe différentes sous-étapes, effectuées de façon itératives, de façon similaire au précédent mode de réalisation. L'indice j désigne le rang de l'itération, cet indice étant initialisé à 1. A chaque indice j est associé une composition de la phase mobile Ck. L'indice j = 1 correspond à la composition de la phase mobile lors de l'adsorption i, les indices suivants correspondant aux paliers de composition 10 suivants Ck+1, Ck+2... Sous-Étape 321 : détermination de l'instant, noté q-1- , auquel a lieu l'adsorption de rang i, et détermination de la composition de la phase mobile à cet instant, par l'application de la fonction d'évolution de la phase mobile Ck = g (tif-1). Cet instant d'adsorption est obtenu de la même façon que durant la sous-étape 221 précédemment décrite. 15 Sous-Étape 322 : estimation d'une durée d'adsorption dans la phase stationnaire, au cours de l'adsorption i, lorsque la phase mobile a la composition Ck, en effectuant un tirage aléatoire conditionné par la deuxième loi de probabilité, de façon à établir un temps ti,s,k passé dans la phase stationnaire, lors de l'adsorption i alors que la phase mobile a pour composition Ck. Cette estimation est réalisée, par tirage aléatoire selon la deuxième loi de probabilité : 20 lorsque j =1, à l'instant tl auquel la Cm' adsorption a lieu, obtenu lors de la sous- étape 321 ; - lorsque j >1, à l'instant ri correspondant au début de du palier Ck associé à l'itération j. Sous-Étape 323 : A partir du tirage de ti,s,k, estimation d'un instant de désorption tid, associé à l'itération j, selon l'expression : 25 ti = ti + ti,s,k (8)d Cette estimation n'est validée que si, à l'instant de désorption tifd, la composition de la phase mobile correspond à la composition Ck prise en compte pour estimer le temps ti,s,k passé dans la phase stationnaire. Autrement dit, cette estimation est validée si il n'y a pas de changement de composition de la phase mobile entre qd et tif. C'est l'objet des sous-étapes 324 et326. 3040216 16 Sous-Etape 324 : estimation de la composition de la phase mobile de la colonne à l'instant de désorption tifd déterminé lors de la sous-étape 323 par application de la fonction d'évolution de la composition de la phase mobile g(ti- d). Sous-Etape 325 : validation de l'estimation de l'instant de désorption tifd estimé lors de l'étape 5 324, en comparant la composition de la phase mobile g(ti- I d) à l'instant de désorption estimé ti!ci et la composition de la phase mobile g(t!) en fonction de laquelle cet instant a été estimé : Si g(ti ! d) # g(ti!), c'est-à-dire si la composition de la phase mobile à l'instant de désorption estimé tifg est différente de la composition Ck à laquelle la durée d'adsorption t,,,,k a été estimée, l'estimation de l'instant de désorption tifd associée à 10 l'itération j n'est pas validée. En effet, cela signifie que la composition de la phase mobile a évolué entre l'instant d'adsorption ti!, auquel l'estimation de la durée d'adsorption ti,s,k est effectuée, et l'instant de désorption qd. Dans ce cas, une nouvelle itération j +1 est effectuée, en repartant de la sous-étape 322, en considérant que l'instant d'adsorption trElcorrespond au début de l'échelon de 15 composition Ck+i suivant l'échelon de composition Ck correspondant à l'itération j. Autrement dit, trn = Eqq:k1Atq , où Atq représente la durée de l'échelon de composition de rang q, et k représente le rang de l'échelon de composition correspondant à l'itération courante : Ck = g(ti!) Les étapes 321 à 325 sont alors réitérées, en prenant en compte une composition de la 20 phase mobile telle que Ck+i = g(trE1). Si g(t! d) = g(q), c'est-à-dire si la composition de la phase mobile à l'instant de désorption estimé ti!ci correspond à la composition de la phase mobile à laquelle la durée d'adsorption t,,,,k a été estimée, l'estimation de l'instant de désorption est validée. On passe alors à l'étape 330, en considérant que le temps passé dans la phase 25 1 stationnaire S au cours de l'adsorption i est tel que ti,s = tig - tig Étape 330 : lors de chaque désorption suivant une adsorption i, réalisation d'un tirage conditionné par la première loi de probabilité, de façon à établir un temps t,,', passé dans la phase mobile entre l'adsorption de rang i et l'adsorption de rang i +1. La première loi de 3040216 17 probabilité peut dépendre de la composition de la phase mobile, du fait de la dépendance du taux d'adsorption ka. Étape 340 : itération suivante ou sortie de l'algorithme. Les étapes 320 à 330 sont réitérées, en incrémentant l'indice i, jusqu'à ce que le cumul du temps passé dans la phase mobile atteigne 5 le temps mort to, de même que dans le mode de réalisation précédent, de façon analogue à l'étape 240. Lorsque cette condition est remplie, l'étape 340 détermine le temps passé par la particule dans la colonne, ou le temps de rétention ajusté, en fonction du cumul des temps passés dans la phase stationnaire. Selon un quatrième mode de réalisation, la composition Cs de la phase stationnaire S évolue 10 spatialement en fonction d'une distance par rapport à l'entrée de la colonne. Ainsi, on peut déterminer une fonction h d'évolution de la phase stationnaire, telle que Cs = h(D,), D, représentant une distance à laquelle se produit une adsorption i par rapport à l'entrée de la colonne. Lors de chaque adsorption, on utilise une deuxième loi de probabilité correspondant à la composition Cs de la phase stationnaire correspondant à la distance D. 15 Essais Des simulations ont été effectuées, en simulant le parcours de 500 molécules d'Anthracène, la modulation en température de la colonne étant réalisée selon dix paliers de température de même durée, entre une température initiale T, = 60 °C à t = 0 min et une température finale Tf = 240 °C à t = 4 min. Les molécules sont introduites dans la colonne à t = 0 min. La figure 4A 20 représente, sur le même graphique : la fonction d'évolution de température Tk = f(t), k étant ici compris entre 1 et 10. le temps cumulé passé, pour chaque réalisation r, dans la phase mobile, jusqu'à une adsorption i, l'indice i évoluant en fonction du temps t. Ce temps cumulé passé dans la phase mobile est noté tinr i , et est déterminé selon l'expression tinr i = EL0 tinr +(11) , 25 tinr 1 désignant le temps dans la phase mobile suite à une adsorption 1, tnir j=o désignant le temps passé dans la phase mobile avant la première adsorption. L'accolade placée sur le graphe représente la plage selon laquelle s'étendent les temps de rétention tr? de chaque réalisation, ces derniers correspondant au temps d'expérience t lorsque le temps désorbé cumulé est égal au temps mort qui, dans cette colonne, est de 60s. 3040216 18 Lors de cet essai, on a considéré une colonne de longueur 30 mètres, de diamètre 0.25 mm, et dont la phase stationnaire, d'épaisseur 0.25 um, a la composition suivante : 5% phényle arylène, 95% diméthylepolysiloxane. La figure 4B représente la distribution du temps de rétention dans la colonne établie à l'aide des 5 500 réalisations. Il s'agit de l'histogramme des temps de rétention tr? pour chaque réalisation r, avec 1 r 500. Cette distribution constitue un modèle direct de la colonne, selon cette évolution de température, pour la molécule d'Anthracène. Plus le nombre de réalisations augmente, plus la précision statistique du modèle direct s'améliore. Par ailleurs, des simulations ont été réalisées en mettant en oeuvre cinq composés de type 10 hydrocarbures aromatiques polycycliques (HAP) : l'acenaphthène (ACE), le fluoranthène (FTN), le benzo(a)pyrène (B(a)P), l'anthracène (ANT) et l'indéno(1,2,3-cd)pyrène (IND). Pour chacune de ces molécules, on a considéré une colonne de longueur 30m, de diamètre 250 um, la phase stationnaire étant constituée d'un mélange 5% phényle arylène, 95% diméthylepolysiloxane d'épaisseur 0,25 um. Pour chaque espèce, on a établi une distribution des temps de rétention 15 ajustés, à partir de laquelle on a estimé une moyenne du temps de rétention ajusté. Ces temps de rétention moyens sont reportés, sur la figure 5, pour chaque molécule, sous la forme de tirets verticaux. Des essais expérimentaux ont été réalisés, en mettant en oeuvre une solution de méthanol comportant les cinq composés précédemment cités, la concentration de chacun d'entre eux 20 étant de 100 ug/mL. La colonne de chromatographie a été couplée à un détecteur par ionisation de flamme, dit Flame Ionisation Detector. Dans les simulations, comme dans les essais expérimentaux, la colonne a été soumise à un gradient de température de 5°C par minute, entre une température initiale de 50°C et une température finale de 250 °C. La durée de chaque palier est d'environ 5s. Le chromatogramme, 25 corrigé du temps mort de la colonne, est représenté sur la figure 5. On remarque la bonne corrélation du temps de rétention ajusté de chaque espèce, entre les essais expérimentaux et les résultats des modélisations. Le procédé décrit ci-dessus s'applique à la modélisation de colonnes de chromatographies aussi bien en phase liquide qu'en phase gazeuse. Il peut également être appliqué à la 30 chromatographie multidimensionnelle, consistant à disposer en série différentes colonnes. Ce procédé peut être mis en oeuvre lors du dimensionnement d'une colonne de chromatographie 3040216 19 ou lors de l'optimisation de paramètres de fonctionnement d'une colonne. Il peut également être utilisé pour établir un modèle direct de la colonne utilisé lors du dépouillement de mesures. REVENDICATIONS1. Procédé d'estimation d'un temps de rétention (tR, tR,) d'une particule (x) dans une colonne de chromatographie (1), ladite colonne de chromatographie comportant un canal (10) dont une paroi est formée d'un revêtement (12), dit phase stationnaire (S), apte à adsorber et désorber ladite particule, cette dernière étant portée par un fluide porteur, dit phase mobile (M), parcourant la dite colonne entre une entrée (in) et une sortie (out), définissant un temps mort (to) le procédé comportant les étapes suivantes : a) estimation (110, 210) d'un temps passé dans la phase mobile (44 après l'entrée de la particule dans la colonne et jusqu'à une première adsorption dans la phase stationnaire, à l'aide d'une loi de probabilité du temps passé par la particule dans ladite phase mobile, dite première loi de probabilité (fm(t, xp)); b) suite à chaque adsorption (i) de la particule dans la phase stationnaire, estimation (120, 220) d'un temps passé (t,,,,) par la particule dans ladite phase stationnaire jusqu'à une désorption, à l'aide d'une loi de probabilité du temps passé par la particule dans la phase stationnaire, dite deuxième loi de probabilité (f.,(t, kd xp)); c) suite à chaque désorption de la phase stationnaire, estimation d'un temps passé de la particule dans la phase mobile (t,,') jusqu'à une nouvelle adsorption dans la phase stationnaire, à l'aide de ladite première loi de probabilité (fm(t, ka x,p)) ; d) répétition des étapes b) et c) jusqu'à la ce que la somme des estimations du temps passé dans la phase mobile soit supérieure audit temps mort (to) ; e) sommation des estimations des temps passés par la particule dans la phase stationnaire (t,,,), afin d'estimer un temps de rétention (tR) de ladite particule dans la colonne. 2. Procédé selon la revendication 1, tel que durant le parcours de la particule (x) entre l'entrée (in) et la sortie (out), un paramètre (P) régissant le fonctionnement de la colonne varie, de telle sorte que la première loi (fin(t, Ica xp))et/ou la deuxième loi (f.,(t, kd xp)) évolue en fonction dudit paramètre. 3. Procédé selon la revendication 2, dans lequel le paramètre (P) peut alors être choisi parmi : la température (T) de la colonne, la composition (C) de la phase mobile (M), ledit paramètre évoluant en fonction du temps selon une pluralité de paliers (Tk, Ck). 3040216 21 4. Procédé selon la revendication 2, dans lequel le paramètre est la composition de la phase stationnaire (S), cette dernière variant entre l'entrée et la sortie de la colonne. 5. Procédé selon l'une quelconque des revendications 2 ou 3, dans lequel durant le parcours de ladite particule, un paramètre (T, C) de la colonne est porté successivement à différents paliers, 5 à chaque palier correspondant une valeur du paramètre (Tk, Ck), à chaque palier étant attribué un rang (k), l'étape b) comprenant alors les sous-étapes suivantes : i) estimation (221, 321) de l'instant (t1-1) auquel a lieu ladite adsorption, et détermination dudit paramètre (Tk, Ck) de la colonne à cet instant ; ii) estimation (222, 322) d'un temps passé (ti,,,k) de la particule dans la phase stationnaire 10 jusqu'à une désorption, à l'aide d'une loi de probabilité du temps passé dans la phase stationnaire, cette loi dépendant d'une valeur du paramètre (Tk, Ck) audit instant d'adsorption. 6. Procédé selon la revendication 5, caractérisé en ce que, l'étape b) comporte également les sous-étapes suivantes : 15 i) estimation (223, 323) d'un instant de désorption (td) de ladite particule, à partir de l'estimation établie lors de la sous-étape ii) ; iv) détermination (224, 324) d'une valeur du paramètre correspondant audit instant de désorption ; v) 20 Lorsque ladite valeur du paramètre audit instant de désorption est différente de la valeur considérée lors de la sous-étape ii), réitération des sous-étapes ii) à iv) en considérant que l'instant d'adsorption (t +1) correspond au début du palier de rang (k+1) suivant lorsque la valeur du paramètre audit instant de désorption correspond à la valeur 25 considérée lors de la sous-étape ii), détermination d'un temps passé (ti,$) dans la phase stationnaire au cours de ladite adsorption. 7. Procédé selon la revendication 6, dans lequel, lors de la sous-étape v), ledit temps passé dans la phase stationnaire est obtenu en fonction : - de la durée de chaque palier (k) entre l'instant d'adsorption, déterminé lors de l'étape i) 30 et l'avant-dernière itération des sous-étapes ii) à v), 3040216 22 - et de l'estimation du temps passé (ti,,,k) dans la phase stationnaire, déterminé lors de ladite dernière itération des sous- étapes ii) à v). 8. Procédé selon l'une quelconque des revendications précédentes, dans lequel lors de l'étape a) ou c), ladite première loi de probabilité suit une loi exponentielle, paramétrée par un taux 5 d'adsorption (ka), ce dernier représentant une probabilité par unité de temps que ladite particule, se trouvant dans la phase mobile, soit adsorbée dans la phase stationnaire. 9. Procédé selon l'une quelconque des revendications précédentes, dans lequel lors de l'étape b), ladite loi de probabilité du temps passé dans la phase stationnaire suit une loi exponentielle, paramétrée par un taux de désorption (kd), ce dernier représentant une probabilité par unité de 10 temps que ladite particule, se trouvant dans la phase stationnaire, soit désorbée de ladite la phase stationnaire. 10. Procédé selon l'une quelconque des revendications précédentes, dans lequel ledit temps mort (to) correspond au temps moyen passé par une particule, transportée par le fluide porteur, dans ladite colonne sans adsorption dans la phase stationnaire. 15 11. Procédé selon l'une quelconque des revendications précédentes, dans lequel les étapes a) à e) sont réalisées pour une pluralité de particules (x) de même espèce (X), de façon à obtenir une distribution d'un temps de rétention (tR, tR') de ladite espèce dans la colonne. 12. Procédé selon l'une quelconque des revendications précédentes, dans lequel les étapes a) à e) sont réalisées en considérant une pluralité de particules de différentes espèces de façon à 20 obtenir une distribution d'un temps de rétention pour chacune de ces espèces. 13. Procédé selon l'une quelconque des revendications 11 ou 12, comportant une étape f) de détermination d'un indicateur statistique d'au moins une distribution du temps de rétention d'une espèce, par exemple une moyenne. 14. Support d'enregistrement d'informations, comportant des instructions pour l'exécution 25 d'un procédé selon l'une quelconque des revendications 1 à 13, ces instructions étant aptes à être exécutées par un processeur (30).