FR 2948438 A1 20110128 FR 0955230 A 20090727 La présente invention concerne un procédé d'estimation du volume d'un réservoir de gaz sous pression. L'invention peut notamment être utilisée dans une station de remplissage de bouteilles ou de réservoirs de gaz sous pression. L'invention peut s'appliquer avantageusement dans les stations de remplissage en gaz carburant de véhicules (hydrogène notamment). Les stations de remplissage de réservoirs de gaz de véhicules automobiles fonctionnant à l'hydrogène sont généralement adaptées à un ou des types déterminés de réservoirs à remplir. II est en effet important de connaître le volume du réservoir à remplir pour réaliser un remplissage optimisé et sûr d'un point de vue de la sécurité. De plus, généralement la mesure de la pression et de la température à l'intérieur du réservoir n'est pas possible ou n'est pas appropriée. Le document WO2005080201A1 décrit une station de remplissage qui estime la densité initiale du gaz dans un réservoir au moyen d'une mesure préliminaire de pression et d'une injection de gaz dont la quantité est déterminée. La densité est estimée via un modèle d'équations obtenues par régression. La station décrite dans ce document réalise également une estimation du volume du réservoir en deux temps : deux jets de gaz sont introduits et, à chaque fois, les pressions sont mesurées et comparées à des rampes de pression de référence pour estimer le volume du réservoir. Par ailleurs, cette station connue contrôle le remplissage uniquement par la masse de gaz transférée. L'estimation du volume réalisée selon ce procédé présente cependant des inconvénients. Outre le fait qu'il est nécessaire de réaliser deux injections de gaz et deux estimations successives, selon ce procédé le volume est estimé par comparaison de la densité estimée avec des tables exprimant le volume en fonction de courbes de densité de référence. De plus, la densité de gaz estimée selon ce document ne permet pas de connaître avec une précision suffisante la quantité de gaz dans le réservoir. Ainsi, la détermination du niveau initial dans le réservoir est faite de façon indirecte. Cela nécessite de connaître à l'avance le type de réservoir (volume, matériau) afin de pouvoir choisir cette référence. Par ailleurs, la quantité initiale de gaz (en nombre de molécules) n'est pas calculée. Le document EP1693611A2 décrit quant à lui une station de remplissage qui détermine la densité de gaz dans le réservoir par une mesure directe dans le réservoir. Cette méthode présente les mêmes inconvénients que ceux discutés précédemment (mesures dans le réservoir inappropriées). Pour réaliser le remplissage selon ce document, le volume du réservoir doit être connu. Ainsi, ce procédé est inutilisable pour un réservoir de volume inconnu. Un but de la présente invention est de pallier tout ou partie des inconvénients de l'art antérieur relevés ci-dessus. A cette fin, le procédé selon l'invention, par ailleurs conforme à la définition générique qu'en donne le préambule ci-dessus, est essentiellement caractérisé en ce que qu'il comprend : - une étape de connexion sélective du réservoir avec au moins une source de gaz sous pression, - une première étape de détermination de la pression (Ptk(tl )) dans le réservoir, - une première étape de détermination de la quantité (ntk(ti)) de gaz contenu dans le réservoir à un premier instant (t=tl), puis - une étape de transfert d'un jet de gaz de contrôle depuis une source vers l'intérieur du réservoir, le jet de contrôle étant réalisé durant une durée discrète (t2-t1), - une seconde étape de détermination de la pression (Ptk(t2)) dans le réservoir après transfert du jet de gaz de contrôle, - une étape de détermination de la quantité (Q=ntk(t2)-ntk(tl )) de gaz transféré de la source vers l'intérieur du réservoir lors du jet de gaz de contrôle, - une étape de calcul du volume (Vtk) du réservoir en appliquant la loi de conservation de l'enthalpie (H) du gaz (l'enthalpie du gaz dans le réservoir après le transfert du jet de gaz de contrôle (ntk(t2).(Htk(t2)) est égale à la somme d'une part de l'enthalpie du gaz dans le réservoir avant le transfert du jet de gaz de contrôle (ntk(tl).Htk(tl )) et, d'autre part, de l'enthalpie du gaz transféré de la source dans le réservoir lors du transfert du jet de gaz de contrôle ((ntk(t2) û ntk(tl)).Hsi(tl )), en considérant que le transfert du jet de gaz de contrôle est adiabatique (Ttk(tl ))=Ttk(t2) et Tsi(tl)=Tsi(t2)), en exprimant l'enthalpie du gaz dans le réservoir comme une fonction uniquement de la température Ttk(ti) du gaz et de la pression Ptk(ti) du gaz (Htk=fonction de (Ttk(ti) ; Ptk(ti)), en utilisant la loi des gaz parfait ou une équation d'état de Van de Waals exprimant le volume du réservoir Vtk(ti) (en mètre cube) uniquement en fonction des variables connues ou estimées précédemment ((Ptk(ti) en pascal, température Ttk(ti) en kelvin et quantité ntk(ti) de gaz en mole), Par ailleurs, des modes de réalisation de l'invention peuvent comporter l'une ou plusieurs des caractéristiques suivantes : - la loi de conservation de l'enthalpie (H) du gaz utilise une modélisation connue de l'enthalpie donnée par l'équation suivante : H=m.T(ti).P(ti) + q avec m=A.T(ti) + B et q = C.T(ti) +D dans lequel T(ti) = température (en Pascal) du gaz au temps t=ti ; P(ti) = pression du gaz (en Pascal) au temps t=ti ; A (en J.Pa 1.K-'), B (en J.Pa'), C (en J.K-1) et D (en J) des coefficients déterminés exprimant l'enthalpie H du gaz en fonction de sa pression et de sa température, - pour le gaz hydrogène, les coefficients A, B, C et D sont donnés par : A est compris entre 10-5 et 10 et de préférence est égal à 4,2131250.10-4, B est compris entre 10-5 et 10 et de préférence est égal à 0,60347812, C est compris entre 10-5 et 100 et de préférence est égal à 14,700936, D est compris entre -1000 et 1000 et de préférence est égal à -507,57537. - la au moins une source de gaz sous pression comprend au moins un conteneur de gaz sous pression de volume connu destiné à transférer du gaz dans le réservoir par équilibrage de pression, on utilise la loi des gaz parfaits pour remplacer la variable température en fonction de la pression du volume et de la quantité dans la loi de conservation de l'enthalpie H, le volume (Vtk) du réservoir étant donné par la formule : Vtk = ù ntk(tl). [B.Ptk(tl,) + D] ù Q. (A.Psi(tl) + C). Psi(tl).Vs i + B.Psi(tl) + D nsi(tl).R (A.Ptk(tl) + C). Ptk(tl) _ (A.Ptk(t2) + C). Ptk(t2) ntk(t2) . [B.Ptk(t2) + D ] (A.Ptk(tl) + C ). Ptk((ttl) _ (A.Ptk(t2) + C ). Ptk((tt2) avec ntk(tl) =la quantité de gaz (en mole) contenu dans le réservoir (1) avant transfert du jet de gaz de contrôle, ntk(t2) =la quantité de gaz (en mole) contenu dans le réservoir (1) après transfert du jet de gaz de contrôle, Ptk(tl)= la pression dans le réservoir (en pascal) déterminée avant transfert du jet de gaz de contrôle, Ptk(t2)= la pression dans le réservoir (en pascal) déterminée après transfert du jet de gaz de contrôle, Q=la quantité (en mole) déterminée de gaz transféré de la source vers l'intérieur du réservoir (1) lors du jet de gaz de contrôle, Psi(tl )= la pression (en pascal) dans la source (Si) déterminée avant transfert du jet de gaz de contrôle, Vsi= le volume connu (en mètre cube) de la source (Si), R = la constante des gaz parfaits en J.K-'.mol-'. - la au moins une source de gaz sous pression comprend au moins un conteneur de gaz sous pression de volume connu destiné à transférer du gaz dans le réservoir par équilibrage de pression, et en ce que, on utilise l'équation d'état de Van der Waals pour exprimer la variable température en fonction de la pression du volume et de la quantité dans la loi de conservation de l'enthalpie H, + - le procédé comporte une étape d'estimation de la température (Ttk(ti)) du gaz dans le réservoir à l'instant (t=ti) en utilisant la volume (Vtk) estimé du réservoir et la loi des gaz parfaits :Ptk(ti).Vtk = ntk(ti).R.Ttk - la première, respectivement la seconde, étape de détermination de la pression dans le réservoir comprend une mesure de pression au niveau d'une conduite reliant sélectivement la source au réservoir, la mesure de pression étant réalisée à l'entrée du réservoir et considérée être la pression à l'intérieur du réservoir, - l'étape de transfert du jet de gaz de contrôle est réalisée pendant une durée (t2-t1) comprise entre 0,1 et 5 seconde et de préférence une durée (t2-t1) comprise entre 0,1 et 3 seconde et encore plus préférentiellement pendant une durée (t2-t1) comprise entre 0,1 et 1 seconde, - l'étape de transfert du jet de gaz de contrôle depuis la source vers l'intérieur du réservoir réalise une montée de pression dans le réservoir compris entre 5.104 Pa et 106 Pa et de préférence entre 2.105 Pa et 8,5.105 Pa, - la au moins une source de gaz sous pression comprend au moins un conteneur de gaz sous pression de volume connu destiné à transférer du gaz dans le réservoir par équilibrage de pression, l'étape de détermination de la quantité (Q=ntk(t2)-ntk(tl )) de gaz transféré de la source vers l'intérieur du réservoir lors du transfert du jet de gaz de contrôle comprenant : une détermination de la pression dans le conteneur (Si) avant (Psi(tl )) et après (Psi(t2)) le transfert du jet de gaz de contrôle, une détermination de la température (Tsi(tl � dans le conteneur (Si) avant le transfert du jet de gaz de contrôle, un calcul de la quantité Q de gaz (en mol) délivré par le conteneur (S1, S2, S3, ...Si) lors du jet de gaz de contrôle en appliquant la loi des gaz parfaits au conteneur (Q=nsi(t2)- nsi(t1)=(Psi(tl ).Vsi)/(R.Tsi(tl )) - (Psi(tl ).Vsi)/(R.Tsi(tl )) et en considérant que le transfert du jet de contrôle est adiabatique pour le conteneur (Tsi(t2))=Tsi(tl)), la quantité Q de gaz (en mole) délivré par le conteneur (S1, S2, 53,...Si) lors du jet de gaz de contrôle étant calculée selon la formule : Q=((Psi(t1)-(Psi(t2)).Vsi)/(R.Tsi(tl )) avec (Psi(tl )) et (Psi(t2)) les pressions (en Pa) dans le conteneur (Si) déterminées respectivement avant (t=t1) et après (t=t2) le transfert du jet de gaz de contrôle, R= la constante des gaz parfaits (en J.K-'.mol-'), Tsi(tl) = la température (en kelvin) dans le conteneur déterminée avant le transfert du jet de contrôle, - la au moins une source de gaz sous pression comprend au moins un conteneur de gaz sous pression de volume connu destiné à transférer du gaz dans le réservoir par équilibrage de pression, l'étape de détermination de la quantité (Qin=ntk(tl)-ntk(tO)) de gaz transféré de la source vers l'intérieur du réservoir lors du jet initial comprenant : une détermination de la pression dans le conteneur (Si) avant (Psi(tO)) et après (Psi(tl )) le transfert du jet initial, une détermination de la température (Tsi(tO)) dans le conteneur (Si) avant le transfert du jet initial, un calcul de la quantité Qin de gaz (en mol) délivré par le conteneur (Si) lors du jet initial en appliquant la loi des gaz parfaits au conteneur : (Qin=nsi(tO)-nsi(tl)=(Psi(tO).Vsi)/(R.Tsi(to)) - (Psi(tl ).Vsi)/(R.Tsi(tl )) et en considérant que le transfert du jet initial est isotherme dans le conteneur (température du gaz dans le conteneur inchangée avant et après le transfert du jet initial Tsi(tO))=Tsi(tl)), la quantité Qin de gaz (en mole) délivré par le conteneur (S1, S2, 53,...Si) lors du jet initial étant calculée selon la formule : Qin=((Psi(t0)-(Psi(tl )).Vsi)/(R.Tsi(to)) avec (Psi(tO)) et (Psi(tl )) les pressions (en Pa) dans le conteneur (S1, S2, S3,...Si) déterminée respectivement avant (t=tO) et après (t=tl) le transfert du jet initial R= la constante des gaz parfaits (en J.K-'.mol-'), Tsi(to) = la température (en kelvin) dans le conteneur (Si) déterminée avant le transfert du jet initial, - la quantité de gaz (ntk(ti)) contenu dans le réservoir est calculée à tout moment (t=ti) en ajoutant à la quantité initiale (ntk(tO)) de gaz contenu dans le réservoir déterminée, la ou les quantités Qi déterminées de gaz transférées par la ou le sources : ntk(ti) = ntk(t0) +Q(i), - l'étape de détermination de la quantité (Q=ntk(t2)-ntk(tl)) de gaz transférée de la source vers l'intérieur du réservoir lors du transfert du jet de contrôle comprend une mesure via un débit mètre au niveau de la connexion entre le réservoir et la source, - l'étape de détermination de la pression initiale (Ptk(tO)) dans le réservoir comprend une mesure de pression (2) au niveau d'une conduite reliant sélectivement une source (S1, S2, S3,...Si) au réservoir (1), la mesure (2) de pression étant réalisée à l'entrée du réservoir (1) et considérée être la pression à l'intérieur du réservoir (1). - l'étape de détermination de la pression (Ptk(tl )) dans le réservoir après (t=tl) le transfert du jet initial comprend une mesure de pression au niveau d'une conduite reliant sélectivement une source au réservoir, la mesure de pression étant réalisée à l'entrée du réservoir et considérée être la pression à l'intérieur du réservoir, - la détermination de la température (Tsi(tO)) dans le conteneur avant le transfert du jet initial comprend une mesure de température ambiante à proximité du conteneur, la température mesurée étant considérée être la température (Tsi(t)) à l'intérieur du conteneur, - la détermination de la pression dans le conteneur avant (Psi(tO)) et après (Psi(tl )) le transfert du jet initial est réalisée par une mesure de pression au niveau d'une conduite reliant sélectivement le conteneur au réservoir, la mesure de pression étant réalisée à la sortie du conteneur et considérée être la pression à l'intérieur du conteneur, - l'étape de détermination de la quantité (Qin=ntk(tl)-ntk(tO)) de gaz transférée de la source vers l'intérieur du réservoir lors du jet initial comprend une mesure via un débit mètre au niveau de la connexion entre le réservoir et la source, - l'étape de transfert du jet initial de gaz depuis une source vers l'intérieur du réservoir réalise une montée de pression dans le réservoir compris entre 5.105 Pa et 106 Pa et de préférence entre 2.105 Pa et 8,5.105.Pa, la connexion du réservoir avec la source comprend au moins une vanne de commande du transfert de gaz et en ce que l'étape de détermination de la quantité (Qin=ntk(tl)-ntk(tO)) de gaz transféré de la source vers l'intérieur du réservoir lors du jet initial comprend une estimation de débit transitant à travers la vanne en fonction de son ouverture, - la connexion du réservoir avec la source comprend au moins une vanne de commande du transfert de gaz et en ce que l'étape de détermination de la quantité (Q=ntk(t2)-ntk(tl)) de gaz transféré de la source vers l'intérieur du réservoir lors du jet de contrôle comprend une estimation de débit transitant à travers la vanne en fonction de son ouverture L'invention peut concerner également tout dispositif ou procédé alternatif comprenant toute combinaison des caractéristiques ci-dessus ou ci-dessous. En particulier, l'invention peut concerner une station de remplissage comprenant un organe électronique de contrôle de calcul et de commande utilisant le procédé selon l'une quelconque des caractéristiques ci-dessus ou ci-après. D'autres particularités et avantages apparaîtront à la lecture de la description ci-après, faite en référence aux figures dans lesquelles : - la figure 1 représente une vue schématique et partielle illustrant la structure et le fonctionnement d'une station de remplissage selon un exemple de réalisation de l'invention, - la figure 2 représente de façon partielle et schématique une succession d'étapes susceptibles d'être mises en oeuvre selon l'invention, - la figure 3 représente de façon partielle et schématique une succession d'étapes susceptibles d'être mises en oeuvre selon l'invention. Comme représenté à la figure 1, le procédé peut être mis en oeuvre par une station de remplissage comprenant, par exemple, plusieurs conteneurs source de gaz sous pression S1, S2...,Si. Classiquement, les conteneurs S1 à Si peuvent être utilisés pour remplir un réservoir 1 par équilibrage de pression (� en cascade � c'est-à-dire en utilisant plusieurs conteneurs successivement ayant des pressions de stockage croissantes). En variante ou en plus, la station peut utiliser un compresseur de gaz ou tout organe susceptible de fournir du gaz à haute pression (entre 200 et 850 bar par exemple). Les conteneurs source S1 à Si sont reliés en parallèle à une conduite 5 commune de remplissage sélectivement connectable à un ou des réservoirs à remplir. Des vannes 7 respectives peuvent être prévues pour contrôler l'alimentation en gaz des conteneurs S1 à Si à la ligne commune 5. De même, une ou plusieurs vannes 6 peut être prévue pour contrôler le flux de gaz dans la conduite 5. Une logique électronique 8 de commande peut être prévue pour réaliser tout ou partie du processus. Pour calculer ou estimer le volume Vtk du réservoir 1, le procédé utilise la quantité ntk(ti) de gaz présent dans le réservoir 1. Cette quantité de gaz ntk(ti) de gaz (par exemple en mole) déjà présent dans le réservoir 1 peut être connue par avance ou peut être calculée préalablement à la méthode d'estimation du volume du réservoir Vtk. Pour calculer le volume Vtk du réservoir le procédé nécessite de transférer jet de gaz dit �jet de gaz de contrôle � ou �jet de contrôle � d'une source Si vers le réservoir 1. Dans un premier temps la station détermine 2 de la pression Ptk(tl) dans le réservoir 1 (cf. étape 201, figure 3). Par exemple, cette détermination de la pression Ptk(tl) dans le réservoir 1 comprend de préférence une mesure de pression (capteur 2) au niveau de la conduite 5 reliée au réservoir 1. C'est-à-dire que la mesure 2 de pression est réalisée à l'entrée/sortie du réservoir 1 et est considérée être la pression à l'intérieur du réservoir 1. Comme explicité ci-dessus, la quantité ntk(tl) de gaz contenu dans le réservoir 1 avant le transfert du jet de contrôle est supposée connue de la station (étape 202, figure 3). La station transfert ensuite le jet de gaz de contrôle depuis une source Si vers l'intérieur du réservoir 1 (étape 203, figure 3). Le jet de contrôle est réalisé durant une durée discrète (t2-t1) comprise par exemple entre 0,1 et 5 seconde ou de préférence entre 0,1 et 3 seconde ou encore plus préférentiellement entre 0,1 et une seconde. Par exemple, le transfert du jet de gaz de contrôle réalise une montée de pression dans le réservoir compris entre 5.104 Pa et 106 Pa et de préférence entre 2.105 Pa et 8,5.105 Pa. La station réalise ensuite une nouvelle détermination 2 de la pression Ptk(t2) dans le réservoir 1 après transfert du jet de gaz de contrôle (étape 204, figure 3). Cette nouvelle détermination 2 de la pression Ptk(t2) dans le réservoir 1 peut bien entendu être réalisée de la même façon que lors de l'étape 201 (avant transfert du jet de contrôle). La station va ensuite déterminer la quantité Q=ntk(t2)-ntk(tl) de gaz transféré de la source Si vers l'intérieur du réservoir 1 lors du jet de gaz de contrôle (étape 205), figure 3). Cette détermination peut être réalisée par différentes méthodes explicitées ci-dessous. Par exemple, la station peut déterminer la quantité Q=ntk(t2)-ntk(tl) de gaz transféré (en mole) de la source Si vers l'intérieur du réservoir 1 lors du jet de contrôle via un débitmètre au niveau de la connexion entre le réservoir 1 et la source Si (par exemple au niveau de la conduite 5). En variante, la détermination de la quantité Q de gaz transféré lors du jet de contrôle est estimée en fonction d'un modèle de vanne 6 (l'ouverture de la vanne 6 permet d'estimer le débit transféré). Selon une autre possibilité, la détermination de la quantité Q transférée lors du jet de contrôle peut être réalisée par calcul en utilisant une équation d'état connue pour le gaz (loi des gaz parfait ou équation d'état de Van der Waals). L'équation d'état permet d'exprimer la quantité de gaz (en mole) en fonction de sa pression (en Pa) et de son volume (en m3) et de sa température (en kelvin). Selon une particularité avantageuse, on considère que le transfert du jet initial est isotherme dans le conteneur (Tsi(tO))=Tsi(tl )). De cette façon, la quantité Q transférée lors du jet de contrôle peut être exprimée uniquement en fonction de variables connues ou estimées précédemment : - Psi(tl )) et Psi(t2)= les pressions (en Pa) dans le conteneur Si déterminées respectivement avant (t=tl) et après (t=t2) le transfert du jet de contrôle, - Tsi(tl) = la température (en kelvin) dans le conteneur déterminée pendant le transfert du jet initial, - Vi le volume (en m3) connu du conteneur Si, - R= la constante des gaz parfaits (en J.K-'.mol-'). La détermination de la température Tsi(tl) dans le conteneur Si pendant le transfert du jet de contrôle peut être réalisée, par exemple, en mesurant 3 la température ambiante à proximité du conteneur Si. Cette température ambiante est considérée être la température Tsi(tl) à l'intérieur du conteneur Si. Bien entendu, cette mesure de température ambiante peut être réalisée pendant ou après le transfert du jet de contrôle (la température peu avant, pendant et peu après le transfert du jet initial étant considérée constante). Comme précédemment pour le réservoir 1, la détermination de la pression dans le conteneur Si avant Psi(tl) et après Psi(t2) le transfert du jet de contrôle peuvent être réalisées chacune par une mesure de pression 4 au niveau de l'entrée/sortie du conteneur Si et sont considérées être les pressions à l'intérieur du conteneur. Ainsi, en utilisant par exemple l'équation d'état des gaz parfaits au niveau du conteneur réalisant le jet initial on obtient : Q= nsi(tl)-nsi(t2)= ((Psi(t1)-(Psi(t2)).Vsi)/(R.Tsi(tl )) avec : R= la constante des gaz parfaits (en J.K-'.mol-'), nsi(tl) et nsi(t2) =les quantités de gaz (en mole) dans le conteneur respectivement avant et après transfert du jet initial Psi(tl) et Psi(t2) = la pression (en Pa) dans le conteneur Si avant et après le transfert du jet initial, Tsi(tl) = la température (en kelvin) dans le conteneur Si avant (=pendant) le transfert du jet initial. Ainsi, cette quantité Q peut être calculée facilement à partir de variables facilement mesurables ou disponibles. Bien entendu, cette quantité de gaz transféré Q transférée lors du jet de contrôle peut être calculée en utilisant l'équation d'état de Van der Waals pour le gaz. La quantité Q de gaz transféré lors du jet de contrôle va permettre de calculer ensuite le volume Vtk du réservoir 1 (étape 206, figure 3). Pour estimer par calcul ce volume Vtk, le procédé utilise la loi de conservation de l'enthalpie H du gaz. C'est-à-dire que l'enthalpie du gaz dans le réservoir 1 après le transfert du jet de gaz de contrôle ntk(t2).Htk(t2) est égale à la somme d'une part de l'enthalpie du gaz dans le réservoir 1 avant le transfert du jet de gaz de contrôle ntk(tl).Htk(tl) et, d'autre part, de l'enthalpie du gaz transféré de la source Si dans le réservoir lors du transfert du jet de gaz de contrôle (ntk(t2) ù ntk(tl)).Hsi(tl ). De plus, selon le procédé, le transfert du jet de gaz de contrôle est considéré adiabatique (Ttk(tl ))=Ttk(t2) et Tsi(tl)=Tsi(t2)). En exprimant classiquement l'enthalpie du gaz dans le réservoir Htk comme une fonction uniquement de la température Ttk(ti) du gaz et de la pression Ptk(ti) du gaz (Htk=fonction de (Ttk(ti) ; Ptk(ti)) et en utilisant la loi des gaz parfait ou une équation d'état de Van de Waals, le procédé permet d'exprimer le volume Vtk du réservoir uniquement en fonction des variables connues ou estimées précédemment. Par exemple, une modélisation connue de l'enthalpie d'un gaz est donnée par l'expression suivante : H=m.T(ti).P(ti) + q avec m=AT(ti) + B et q = CT(ti) +D ((Ptk(ti)=pression en pascal, Ttk(ti)= température en kelvin et ntk(ti) = quantité en mole). dans lequel T(ti) = température (en Pascal) du gaz au temps t=ti ; P(ti) = pression du gaz (en Pascal) au temps t=ti ; A (en J.Pa 1.K-'), B (en J.Pa'), C (en J.K-1) et D (en J) étant des coefficients calculés ou bien connus par des modèles ou tables exprimant l'enthalpie H du gaz en fonction de sa pression et de sa température. Par exemple, pour le gaz hydrogène (données � NIST �), les coefficients A, B, C et D peuvent être donnés par A est compris entre 10-5 et 10 et de préférence est égal à 4,2131250.10-4, B est compris entre 10-5 et 10 et de préférence est égal à 0,60347812, C est compris entre 10-5 et 100 et de préférence est égal à 14,700936, D est compris entre -1000 et 1000 et de préférence est égal à -507,57537. En utilisant cette expression de l'enthalpie avant (t=tl) et après (t=t2) le jet de contrôle on a pour le réservoir (Htk(ti)= enthalpie du réservoir à l'instant ti) Htk(tl)=(A.Ptk((tl) + C).Ttk(tl) + B.Ptk(tl) + D Htk(t2) = (A.Ptk(t2) + C).Ttk(t2) + B.Ptk(t2) + D et pour le conteneur Si (Hsi(ti)= enthalpie du conteneur à l'instant ti) : HSi(tl)=(A.Psi((tl) + C).Tsi(tl) + b.Psi(tl) + D HSi(t2) = (A.Psi(t2) + C).Tsi(t2) + B.Psi(t2) + D Comme vu ci-dessus, la conservation de l'enthalpie peut s'écrire : ntk(t2).Htk(t2) = ntk(tl).Htk(tl) + (ntk(t2) ûntk(tl)).Hsi(tl ) En remplaçant l'enthalpie H par son expression ci-dessus (H= mT.P +q) et en remplaçant la température en fonction de la pression du volume et de la quantité (équation des gaz parfaits T=P.V/R.T), on obtient Vtk = ù ntk(tl). [B.Ptk(tl,) + D] ù Q. (A.Psi(tl) + C). Psi(tl).Vs i + B.Psi(tl) + D nsi(tl).R (A.Ptk(tl) + C). Ptk(tl) _ (A.Ptk(t2) + C). Ptk(t2) + ntk(t2) . [B.Ptk(t2) + D ] (A.Ptk(tl) + C ). Ptk((ttl) _ (A.Ptk(t2) + C ). Ptk((tt2) avec ntk(tl) =la quantité de gaz (en mole) contenu dans le réservoir (1) avant transfert du jet de gaz de contrôle, ntk(t2) =la quantité de gaz (en mole) contenu dans le réservoir (1) après transfert du jet de gaz de contrôle, Ptk(tl)= la pression dans le réservoir (en pascal) déterminée avant transfert du jet de gaz de contrôle, Ptk(t2)= la pression dans le réservoir (en pascal) déterminée après transfert du jet de gaz de contrôle, Q=la quantité (en mole) déterminée de gaz transféré de la source vers l'intérieur du réservoir (1) lors du jet de gaz de contrôle, Psi(tl )= la pression (en pascal) dans la source (Si) déterminée avant transfert du jet de gaz de contrôle, Vsi= le volume connu (en mètre cube) de la source (Si), R = la constante des gaz parfaits en J.K-'.mol-'. De cette façon, le procédé permet ainsi d'estimer avec une bonne précision le volume Vtk du réservoir à remplir. Connaissant le volume Vtk du réservoir, la station de remplissage peut aisément calculer la température Ttk(ti) du gaz dans le réservoir à tout instant (avec l'équation d'état choisie et via la pression (Ptk(ti) et la quantité de gaz dans le réservoir (ntkti) précédemment calculées ou mesurées). En variante, dans la loi de conservation de l'enthalpie, on peut utiliser l'équation d'état de Van der Waals pour remplacer la température en fonction de la pression du volume et de la quantité. Dans ce cas, le volume Vtk du réservoir 1 est donné par une formule faisant intervenir les variables connues ou déterminées précédemment. Sa résolution pourra être réalisée par exemple par calcule numérique. On comprend que le procédé selon l'invention permet d'adapter la station de remplissage à tout type de réservoir ou bouteille à remplir, même des réservoirs inconnus (volume). 10 Le procédé permet ainsi d'optimiser les remplissages (en taux de remplissage et en vitesse de remplissage), en permettant également de contrôler la température estimée dans le réservoir. Le procédé ne nécessite aucune communication (mesures de température, de volume, de pression) à l'intérieur d'un réservoir ou d'un véhicule portant ce réservoir. Un seul jet de gaz permet d'estimer avec précision le volume. Dans le procédé ci-dessus on a indiqué que la quantité de gaz contenu dans le réservoir (avant transfert du jet de contrôle) devait être connue. On va décrire ci-dessous une méthode avantageuse (et facultative) d'estimation de la quantité de gaz ntk(ti) contenu dans un réservoir 1 de gaz sous pression de volume Vtk inconnu. Estimation préalable de la quantité de qaz contenu dans le réservoir : Cette méthode peut notamment être utilisée avant la mise en oeuvre du procédé de calcul du volume Vtk décrit ci-dessus (c'est-à-dire avant le transfert du jet de contrôle). Pour estimer la quantité de gaz contenu dans le réservoir, le procédé peut comprendre une étape de détermination 2 de la pression initiale Ptk(t0) dans le réservoir 1 (étape 101, figure 2). Cette détermination de la pression initiale Ptk(t0) dans le réservoir 1 peut comprendre une mesure de pression (capteur 2) au niveau de la conduite 5 reliée au réservoir 1. C'est-à-dire que la mesure 2 de pression est réalisée à l'entrée/sortie du réservoir 1 et est considérée être la pression à l'intérieur du réservoir 1. Le réservoir 1 est connecté sélectivement à au moins une source Si de gaz sous pression et la station réalise un transfert d'un jet initial de gaz depuis la source Si vers l'intérieur du réservoir 1 (étape 102, figure 2). A l'instant t=t0, le jet initial est réalisé durant une durée discrète (t1-t0) comprise entre 0,1 et 5 seconde et de préférence entre 0,1 et 3 seconde et encore plus préférentiellement entre 0,1 et 1 seconde. Par exemple, l'étape 102 de transfert du jet initial de gaz depuis la source Si vers l'intérieur du réservoir 1 réalise une montée de pression dans le réservoir 1 limitée et comprise entre 5.104 Pa et 106 Pa et de préférence entre 2.105 Pa et 8,5.105 Pa. La station détermine ensuite la pression Ptk(tl) dans le réservoir 1 après transfert du jet initial (étape 13, figure 2). Comme précédemment, l'étape de détermination de la pression Ptk(tl) dans le réservoir 1 après (t=tl) le transfert du jet initial peut comprendre une mesure de pression 2 au niveau de l'entrée/sortie du réservoir 1. La station détermine d'abord la quantité Qin=ntk(tl)-ntk(t0) de gaz transféré de la source Si vers l'intérieur du réservoir 1 lors du jet initial (étape 104, figure 2). Cette étape de détermination de la quantité Qin de gaz transférée de la source Si vers l'intérieur du réservoir 1 peut également être réalisée via un débitmètre au niveau de la connexion entre le réservoir 1 et la source Si (par exemple au niveau de la conduite 5). En variante, la détermination de la quantité Qin de gaz transféré lors du jet initial est estimée en fonction d'un modèle de vanne 6 (l'ouverture de la vanne 6 permet d'estimer le débit transféré). Selon une autre possibilité, la détermination de la quantité Qin transférée lors du jet initial peut être réalisée par calcul en utilisant une équation d'état connue pour le gaz (loi des gaz parfait ou équation d'état de Van der Waals). L'équation d'état permet d'exprimer la quantité de gaz (en mole) en fonction de sa pression (en Pa) et de son volume (en m3) et de sa température (en kelvin). Selon une particularité avantageuse, on considère que le transfert du jet initial est isotherme dans le conteneur (Tsi(tO))=Tsi(tl)). De cette façon, la quantité Qin transférée lors du jet initial peut être exprimée uniquement en fonction de variables connues ou estimées précédemment : - Psi(tO)) et Psi(tl )= les pressions (en Pa) dans le conteneur Si déterminées respectivement avant (t=tO) et après (t=tl) le transfert du jet initial, - Tsi(t0) = la température (en kelvin) dans le conteneur déterminée avant le transfert du jet initial, - Vi le volume (en m3) connu du conteneur Si, - R= la constante des gaz parfaits (en J.K-'.mol-'). La détermination de la température Tsi(t0) dans le conteneur Si avant le transfert du jet initial peut être réalisée, par exemple, en mesurant la température ambiante à proximité du conteneur Si. Cette température ambiante est considérée être la température Tsi(t0) à l'intérieur du conteneur Si. Bien entendu, cette mesure de température ambiante peut être réalisée pendant ou après le transfert du jet initial (la température peu avant, pendant et peu après le transfert du jet initial étant considérée constante). Comme précédemment pour le réservoir 1, la détermination de la pression dans le conteneur Si avant Psi(t0) et après Psi(tl) le transfert du jet initial peuvent être réalisées chacune par une mesure de pression 4 au niveau de l'entrée/sortie du conteneur Si et sont considérées être les pressions à l'intérieur du conteneur. Ainsi, en utilisant par exemple l'équation d'état des gaz parfaits au niveau du conteneur réalisant le jet initial on obtient : Qin= nsi(t0)-nsi(tl )= ((Psi(t0)-(Psi(tl )).Vsi)/(R.Tsi(tO)) (équation I) avec : R= la constante des gaz parfaits (en J.K-1.mor1), nsi(t0) et nsi(tl) =les quantités de gaz (en mole) dans le conteneur respectivement avant et après transfert du jet initial Psi(t0) et Psi(tl) = la pression (en Pa) dans le conteneur Si avant et après le transfert du jet initial, Tsi(t0) = la température (en kelvin) dans le conteneur Si avant (=pendant) le transfert du jet initial. Bien entendu, cette quantité de gaz transféré Qin transférée lors du jet initial peut être calculée en utilisant l'équation d'état de Van der Waals pour le gaz. Lorsque cette quantité de gaz transféré Qin est disponible, la station peut ensuite calculer la quantité initiale ntk(t0) de gaz contenu dans le réservoir 1 avant le transfert du jet initial (étape 105, figure 2). En effet, cette quantité initiale ntk(t0) peut être calculée en appliquant une équation d'état connue pour le gaz (équation des gaz parfaits ou équation d'état de Van der Waals) exprimant la quantité de gaz en mole en fonction de sa pression (en Pa) de son volume (en m3) et de sa température (en kelvin). Le procédé considère de plus que le transfert du jet initial est isotherme dans le réservoir Ttk(tO))=Ttk(tl) et utilise la quantité Qin de gaz transférée déterminée à l'étape précédente. Plus précisément, la quantité initiale ntk(t0) de gaz (en mole) peut être calculée uniquement en fonction de variables connues ou déjà estimées : Ptk(t0)= la pression initiale (en Pascal) dans le réservoir Qin = la quantité déterminée de gaz (en mole) transféré lors du jet initial, Ptk(tl)= la pression de gaz (en Pascal) dans le réservoir déterminée après transfert du jet initial de gaz. Par exemple, en utilisant l'équation des gaz parfaits pour le gaz dans le réservoir aux instants tO et t1, on peut écrire Ptk(tO).Vtk =ntk(tO).R.Ttk(tO) et Ptk(tl ).Vtk =ntk(tl ).R.Ttk(tl) avec Vtk = le volume (en m3) inconnu du réservoir Ttk(ti) = la température (en kelvin) inconnue dans le réservoir à l'instant ti ntk(ti) = la quantité (en mole) de gaz dans le réservoir à l'instant ti. Du fait du transfert isotherme lors du jet initial (Ttk(tO))=Ttk(tl )), on a dès lors : ntk(tl )- ntk(tO)=((Ptk(t1)-Ptk(tO)).Vtk)/(R.Ttk(tO)) (équation II) En revenant à l'équation des gaz parfaits dans le réservoir et en remplaçant Vtk/Ttk(tO) en fonction de Qin et (Psi(tl)-Psi(tO)), on obtient : ntk(tO)=(Ptk(tO).Vtk)/(R.Ttk(tO)) = (Qin.Ptk(tO))/(Ptk(tl)-Ptk(tO)) (équation III) En utilisant l'équation (I) et l'équation (III) on peut donc calculer la quantité de gaz dans le réservoir1 à tout instant t=ti en ajoutant à la quantité initiale calculée ntk(t0) la ou les quantités ajoutées par le ou les conteneurs Si. En variante, le calcul de la quantité initiale ntk(t0) de gaz contenu dans le réservoir 1 avant le transfert du jet initial peut être réalisé avec l'équation d'état de Van der Waals : _ ntk (ti )'R'T k (ti ) _ avdw Pk (ti )' (ti l Vtk ù bvdw'ntk ) Vtk ntk (ti ) avec aVdW = un coefficient connu, compris pour le cas de hydrogène entre ... et... et de préférence égal à 1,1900982.10-2 Pa.m2.mol-'. 2 10 15 bVdW = un coefficient connu, compris pour le cas de l'hydrogène entre ... et... et de préférence égal à 1,7936793.10-5 Pa.m3.mo1-'. Dans ce cas également, la quantité initiale ntk(t0) de gaz (en mole) peut être calculée uniquement en fonction de variables connues ou déjà estimées Avec l'équation d'état de Van der Waals, on peut également exprimer la pression Ptk(ti) du gaz au temps ti dans le réservoir 1 en fonction de la quantité ntk(ti) de gaz, du volume Vtk et de la température Ttk(ti) selon la formule : pk (ti) = ntk (ti ).R.l tk (ti )'Vtk2 ù avdw'ntk (ti )2 (Vtk ù bvdw'ntk (ti )) Vtk 2 . (Vtk ù bvdW .ntk (ti )) On peut alors exprimer la quantité de gaz par l'expression polynomiale suivante : `ti )3 .bvdw .( k `ti ).Vtk 2 + avdw) ù tk `ti ) 2.avdw .Vtk + ntk `ti ).Vtk 2.R.Ttk ù Ptk `ti ).Vtk3 = 0 En utilisant les désignations suivantes : 2 a = bvdw•\Pk (ti)•Vtk + avdw ; b = ùavdw. d = ùPk(ti).Vk3; .x = ntk(ti) On obtient l'expression explicite suivante : 1 9 {: b c. 2 .d 27a 4.9a bc.--.2b t'l 20 Comme précédemment, en utilisant la quantité Qin transférée lors du jet de gaz initial, la quantité initiale ntk(t0) de gaz contenu dans le réservoir 1 avant le transfert du jet initial peut être exprimée en fonction uniquement de variables connues ou déjà estimées (Ptk(tO), Qin, Ptk(tl)). Connaissant la quantité de gaz initiale (ntk(tO), la station peut ensuite calculer la quantité contenue dans le réservoir 1 à tout instant ti ntk(ti). En effet, il suffit d'ajoutant à la quantité initiale 25 ntk(t0) de gaz estimée, la ou les quantités de gaz transférées par la ou le sources (Si, S2, S3, ...,Si). On comprend que le procédé selon l'invention permet d'adapter la station de remplissage à tout type de réservoir ou bouteille à remplir, même des réservoirs inconnus (volume et contenu initial inconnus). Le procédé permet ainsi d'optimiser les remplissages (en taux de remplissage et en vitesse de remplissage), en permettant également de contrôler la température estimée dans le réservoir. Le procédé ne nécessite aucune communication (mesures de température, de volume, de pression) à l'intérieur d'un réservoir ou d'un véhicule portant ce réservoir. Un seul jet de gaz permet d'estimer la quantité initiale contenue dans le réservoir et un seul jet de gaz permet d'estimer le volume. REVENDICATIONS1. Procédé d'estimation du volume (Vtk) d'un réservoir (1) de gaz sous pression comprenant : - une étape de connexion sélective du réservoir (1) avec au moins une source (S1, S2, S3, ...Si) de gaz sous pression, - une première étape de détermination (2) de la pression (Ptk(tl )) dans le réservoir (1), - une première étape de détermination de la quantité (ntk(ti)) de gaz contenu dans le réservoir (1) à un premier instant (t=t1), puis - une étape de transfert d'un jet de gaz de contrôle depuis une source (S1, S2, S3, ...,Si) vers l'intérieur du réservoir (1), le jet de contrôle étant réalisé durant une durée discrète (t2-t1), - une seconde étape de détermination (2) de la pression (Ptk(t2)) dans le réservoir (1) après transfert du jet de gaz de contrôle, - une étape de détermination de la quantité (Q=ntk(t2)-ntk(tl)) de gaz transféré de la source (S1, S2, S3, ...,Si) vers l'intérieur du réservoir (1) lors du jet de gaz de contrôle, - une étape de calcul du volume (Vtk) du réservoir (1) en appliquant la loi de conservation de l'enthalpie (H) du gaz (l'enthalpie du gaz dans le réservoir (1) après le transfert du jet de gaz de contrôle (ntk(t2).(Htk(t2)) est égale à la somme d'une part de l'enthalpie du gaz dans le réservoir (1) avant le transfert du jet de gaz de contrôle (ntk(tl).Htk(tl )) et, d'autre part, de l'enthalpie du gaz transféré de la source (Si) dans le réservoir lors du transfert du jet de gaz de contrôle ((ntk(t2) ûntk(tl)).Hsi(tl )), en considérant que le transfert du jet de gaz de contrôle est adiabatique (Ttk(tl ))=Ttk(t2) et Tsi(tl)=Tsi(t2)), en exprimant l'enthalpie du gaz dans le réservoir comme une fonction uniquement de la température Ttk(ti) du gaz et de la pression Ptk(ti) du gaz (Htk=fonction de (Ttk(ti) ; Ptk(ti)), en utilisant la loi des gaz parfait ou une équation d'état de Van de Waals exprimant le volume du réservoir Vtk(ti) (en mètre cube) uniquement en fonction des variables connues ou estimées précédemment ((Ptk(ti) en pascal, température Ttk(ti) en kelvin et quantité ntk(ti) de gaz en mole). 2. Procédé selon la revendication 1, caractérisé en ce que la loi de conservation de l'enthalpie (H) du gaz utilise une modélisation connue de l'enthalpie donnée par l'équation suivante : H=m.T(ti).P(ti) + q avec m=A.T(ti) + B et q = C.T(ti) +D dans lequel T(ti) = température (en Pascal) du gaz au temps t=ti ; P(ti) = pression du gaz (en Pascal) au temps t=ti ; A (en J.Pa 1.K-'), B (en J.Pa'), C (en J.K-1) et D (en J) des coefficients déterminés exprimant l'enthalpie H du gaz en fonction de sa pression et de sa température. 3. Procédé selon la revendication 2, caractérisé en ce que, pour le gaz hydrogène, les coefficients A, B, C et D sont donnés par A est compris entre 10-5 et 10 et de préférence est égal à 4,2131250.10-4, B est compris entre 10-5 et 10 et de préférence est égal à 0,60347812, C est compris entre 10-5 et 100 et de préférence est égal à 14,700936, D est compris entre -1000 et 1000 et de préférence est égal à -507,57537. 4. Procédé selon l'une quelconque des revendications 2 ou 3, caractérisé en ce que la au moins une source de gaz sous pression comprend au moins un conteneur (Si, S2, S3, ...Si) de gaz sous pression de volume connu (Vs1, Vs2, Vs3, ...Vsi) destiné à transférer du gaz dans le réservoir (1) par équilibrage de pression, et en ce que, on utilise la loi des gaz parfaits pour remplacer la variable température en fonction de la pression du volume et de la quantité dans la loi de conservation de l'enthalpie H, le volume (Vtk) du réservoir (1) étant donné par la formule : Vtk = ù ntk(tl). [B.Ptk(tl 1) + D ]ù Q. (A.Psi(tl) + C). Psi(tl).Vs + B.Psi(tl) + D nsi(tl).R (A.Ptk(tl) + C). PtkR(tl) _ (A.Ptk(t2) + C). PtkR(t2) ntk(t2) . [B.Ptk(t2) + D ] (A.Ptk(tl) + C) . Ptk(ttl) ù (A.Ptk(t2) + C ). Ptk(tt2) avec ntk(tl) =la quantité de gaz (en mole) contenu dans le réservoir (1) avant transfert du jet de 20 gaz de contrôle, ntk(t2) =la quantité de gaz (en mole) contenu dans le réservoir (1) après transfert du jet de gaz de contrôle, Ptk(t1)= la pression dans le réservoir (en pascal) déterminée avant transfert du jet de gaz de contrôle, 25 Ptk(t2)= la pression dans le réservoir (en pascal) déterminée après transfert du jet de gaz de contrôle, Q=la quantité (en mole) déterminée de gaz transféré de la source vers l'intérieur du réservoir (1) lors du jet de gaz de contrôle, Psi(t1)= la pression (en pascal) dans la source (Si) déterminée avant transfert du jet de gaz 30 de contrôle, Vsi= le volume connu (en mètre cube) de la source (Si), + R = la constante des gaz parfaits en J.K-'.mol-'. 5. Procédé selon l'une quelconque des revendications 2 ou 3, caractérisé en ce que la au moins une source de gaz sous pression comprend au moins un conteneur (Si, S2, S3, ...Si) de gaz sous pression de volume connu (Vsi, Vs2, Vs3, ...Vsi) destiné à transférer du gaz dans le réservoir (1) par équilibrage de pression, et en ce que, on utilise l'équation d'état de Van der Waals pour exprimer la variable température en fonction de la pression du volume et de la quantité dans la loi de conservation de l'enthalpie H. 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce qu'il comporte une étape d'estimation de la température (Ttk(ti)) du gaz dans le réservoir (1) à l'instant (t=ti) en utilisant la volume (Vtk) estimé du réservoir (1) et la loi des gaz parfaits : Ptk(ti).Vtk = ntk(ti).R.Ttk. 7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que la première, respectivement la seconde, étape de détermination de la pression dans le réservoir (1) comprend une mesure de pression (2) au niveau d'une conduite reliant sélectivement la source (S1, S2, S3, ...Si) au réservoir (1), la mesure (2) de pression étant réalisée à l'entrée du réservoir (1) et considérée être la pression à l'intérieur du réservoir (1). 8. Procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce que l'étape de transfert du jet de gaz de contrôle est réalisée pendant une durée (t2-t1) comprise entre 0,1 et 5 seconde et de préférence une durée (t2-t1) comprise entre 0,1 et 3 seconde et encore plus préférentiellement pendant une durée (t2-t1) comprise entre 0,1 et 1 seconde. 9. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que l'étape de transfert du jet de gaz de contrôle depuis la source vers l'intérieur du réservoir (1) réalise une montée de pression dans le réservoir (1) compris entre 5.104 Pa et 106 Pa et de préférence entre 2.105 Pa et 8,5.105 Pa. 10. Procédé selon l'une quelconque des revendications 1 à 9, caractérisé en ce que la au moins une source de gaz sous pression comprend au moins un conteneur (S1, S2, S3, ...Si) de gaz sous pression de volume connu (Vsi, Vs2, Vs3, Vsi) destiné à transférer du gaz dans le réservoir par équilibrage de pression, l'étape de détermination de la quantité (Q=ntk(t2)-ntk(t1)) de gaz transféré de la source (S1, S2, S3, ...Si) vers l'intérieur du réservoir (1) lors du transfert du jet de gaz de contrôle comprenant : - une détermination de la pression dans le conteneur (Si) avant (Psi(t1)) et après (Psi(t2)) le transfert du jet de gaz de contrôle, - une détermination de la température (Tsi(t1)) dans le conteneur (Si) avant le transfert du jet de gaz de contrôle, - un calcul de la quantité Q de gaz (en mol) délivré par le conteneur (S1, S2, S3, ...Si) lors du jet de gaz de contrôle en appliquant la loi des gaz parfaits au conteneur : (Q=nsi(t2)-nsi(t1)=(Psi(t1).Vsi)/(R.Tsi(t1)) - (Psi(t1).Vsi)/(R.Tsi(t1)) et en considérant que le transfert du jet de contrôle est adiabatique pour le conteneur (Tsi(t2))=Tsi(t1)), la quantité Q de gaz (en mole) délivré par le conteneur (Si, S2, 53,...Si) lors du jet de gaz de contrôle étant calculée selon la formule : Q=((Psi(t1)-(Psi(t2)).Vsi)/(R.Tsi(t1)) avec (Psi(t1)) et (Psi(t2)) les pressions (en Pa) dans le conteneur (Si) déterminées respectivement avant (t=t1) et après (t=t2) le transfert du jet de gaz de contrôle, R= la constante des gaz parfaits (en J.K-'.mol-'), Tsi(t1) = la température (en kelvin) dans le conteneur déterminée avant le transfert du jet de contrôle. 11. Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce que la première étape de détermination de la quantité (ntk(t1)) de gaz contenu dans le réservoir (1) à un premier instant (t=t1) comprend : - une connexion sélective du réservoir (1) avec au moins une source (S1, S2, 53,...Si) de gaz sous pression, - une détermination de la pression initiale (Ptk(tO)) dans le réservoir, - un transfert d'un jet de gaz initial depuis une source (S1, S2, 53,...Si) vers l'intérieur du réservoir (1), le jet initial étant réalisé durant une durée discrète (t1-t0), - une détermination de la pression (Ptk(t1)) dans le réservoir après transfert du jet initial, - une détermination de la quantité (Qin=ntk(t1)-ntk(tO)) de gaz transféré de la source (Si, S2, 53,...Si) vers l'intérieur du réservoir (1) lors du jet initial, - une étape de calcul de la quantité initiale (ntk(tO)) de gaz contenu dans le réservoir (1) avant le transfert du jet initial, la quantité initiale (ntk(tO)) de gaz contenu dans le réservoir (1) étant calculée en appliquant la loi des gaz parfaits au réservoir (ntk(t0)=(Ptk(tO).Vtk)/(R.Ttk(to)) et en considérant que le transfert du jet initial est isotherme pour le réservoir (Ttk(t0))=Ttk(t1)) et en utilisant la quantité (Qin) de gaz transféré déterminée à l'étape précédente, la quantité initiale (ntk(tO)) de gaz étant calculée selon une formule : ntk(tO)=(Ptk(tO).Qin)/(Ptk(t1)-Ptk(tO)), dans laquelle ntk(t0) = la quantité (en mole) de gaz dans le réservoir avant transfert du jet initial de gaz, Ptk(t0)= la pression initiale (en Pascal) dans le réservoir déterminée avant transfert du jet initial de gaz ; Ptk(t1)= la pression de gaz (en Pascal) dans le réservoir déterminée après transfert du jet initial de gaz et Qin= la quantité déterminée de gaz (en mole) transféré de la source (S1, S2, 53,...Si) vers l'intérieur du réservoir (1) lors du jet initial. 12. Procédé selon la revendication 10, caractérisé en ce que la au moins une source de gaz sous pression comprend au moins un conteneur (S1, S2, 53,...Si) de gaz sous pression de volume connu (Vs1, Vs2, Vs3, ...Vsi) destiné à transférer du gaz dans le réservoir par équilibrage de pression, l'étape de détermination de la quantité (Qin=ntk(t1)-ntk(tO)) de gaz transféré de la source vers l'intérieur du réservoir (1) lors du jet initial comprenant : - une détermination de la pression dans le conteneur (Si) avant (Psi(tO)) et après (Psi(tl )) le transfert du jet initial, - une détermination de la température (Tsi(tO)) dans le conteneur (Si) avant le transfert du jet initial, - un calcul de la quantité Qin de gaz (en mol) délivré par le conteneur (Si) lors du jet initial en appliquant la loi des gaz parfaits au conteneur : (Qin=nsi(t0)-nsi(t1)=(Psi(tO).Vsi)/(R.Tsi(to)) - (Psi(tl ).Vsi)/(R.Tsi(tl )) et en considérant que le transfert du jet initial est isotherme dans le conteneur (température du gaz dans le conteneur inchangée avant et après le transfert du jet initial Tsi(tO))=Tsi(tl )), la quantité Qin de gaz (en mole) délivré par le conteneur (S1, S2, 53,...Si) lors du jet initial étant calculée selon la formule : Qin=((Psi(t0)-(Psi(tl )).Vsi)/(R.Tsi(to)) avec (Psi(tO)) et (Psi(tl)) les pressions (en Pa) dans le conteneur (S1, S2, 53,...Si) déterminée respectivement avant (t=tO) et après (t=tl) le transfert du jet initial R= la constante des gaz parfaits (en J.K-'.mol-'), Tsi(to) = la température (en kelvin) dans le conteneur (Si) déterminée avant le transfert du jet initial. 13. Procédé selon l'une quelconque des revendications 1 à 12, caractérisé en ce que la quantité de gaz (ntk(ti)) contenu dans le réservoir (1) est calculée à tout moment (t=ti) en ajoutant à la quantité initiale (ntk(tO)) de gaz contenu dans le réservoir (1) déterminée, la ou les quantités (Q(i)) déterminées de gaz transférées par la ou le sources : ntk(ti) = ntk(t0) +Q(i). 14. Procédé selon l'une quelconque des revendications 1 à 13, caractérisé en ce que l'étape de détermination de la quantité (Q=ntk(t2)-ntk(tl )) de gaz transférée de la source (S1, S2, 53,...Si) vers l'intérieur du réservoir (1) lors du transfert du jet de contrôle comprend une mesure via un débit mètre au niveau de la connexion entre le réservoir (1) et la source (S1, S2, 53,...Si).