FR 2978233 A1 20130125 FR 1156653 A 20110722 La présente invention concerne un procédé de remplissage d'un réservoir avec du gaz sous pression. L'invention concerne plus particulièrement un procédé de remplissage d'un réservoir avec du gaz sous pression, notamment de l'hydrogène gazeux, pour atteindre un taux de remplissage cible prédéterminé, le remplissage étant interrompu lorsqu'une grandeur physique mesurée ou estimée du gaz dans le réservoir correspond au taux de remplissage cible ou lorsque la température dans le réservoir atteint un seuil maximum déterminée, le procédé comprenant : - une étape de mesure de la pression initiale (P(t0)) dans le réservoir avant le remplissage, - une étape de détermination de la quantité (m(t0)) initiale de gaz dans le réservoir avant remplissage, - une étape mesure de la pression (P(ti)) courante dans le réservoir lors du remplissage, - une étape de détermination de la quantité (Q(ti)) courante de gaz transférée dans le réservoir lors du remplissage, - une étape de calcul de la quantité (m(ti)) courante de gaz dans le réservoir lors du remplissage, une étape de détermination de la température (T(ti)) courante du gaz dans le réservoir lors du remplissage. L'invention s'applique de façon préférentielle aux remplissages rapides (c'est-à-dire de l'ordre de trois à quinze minutes par exemple) de réservoirs de gaz contenant de l'hydrogène à des pressions élevées (par exemple comprise entre 300 et 850 bar). Le remplissage de réservoirs d'hydrogène gazeux à haute pression de véhicules automobiles est rendu difficile du fait de nombreuses contraintes techniques. En effet, le remplissage doit être optimisé aussi bien du point de vue de la quantité transférée que de la durée de remplissage, sans pour autant générer un échauffement du réservoir incompatible avec sa structure. De plus, pour assurer un remplissage optimal, il est préférable de connaître les caractéristiques géométriques et structurelles du réservoir (qui peuvent varier d'un véhicule à l'autre). Ainsi, dans le cas d'un remplissage « sans communication », les données notamment relatives à la géométrie du réservoir des véhicules et à la quantité de gaz restante dans le réservoir ne sont pas transmises à la station de remplissage. Il en résulte que la station doit être programmée pour un type prédéterminé de réservoir ou doit calculer et estimer des données manquantes (cf. par exemple les documents FR2948438A1 et FR2948437A1). Dans le cas d'un remplissage « avec communication », le véhicule transmet tout ou partie de ces informations permettant d'optimiser le remplissage (cf. par exemple le protocole décrit dans le document SAE J 2799) Un but de l'invention est de proposer un procédé de remplissage amélioré répondant aux contraintes connues et pouvant s'appliquer aussi bien aux remplissages avec ou sans communication. De nombreux documents tels que EP1205704A1, US5628349, US5752552 et EP1336795 décrivent des procédés de remplissage utilisant la densité comme paramètre de contrôle du remplissage. Le document US6786245 décrit un procédé de remplissage dans lequel la température et la densité du gaz sont calculées à partir de la température ainsi qu'à partir de la pression et de la composition du gaz. La densité est calculée à partir du facteur de compressibilité à partir des équations du second ordre utilisant des coefficients viriels appliquées à l'équation d'état du gaz (cette méthode est bien connue notamment de l'article « The equation of state of neon between 27 and 70 K » de R.M. Gibbons Gas Council, London Research Station, Michael Road, London SW6, UK (1969)). Le document US6786245 ne donne cependant pas de méthode satisfaisante de contrôle de la température dans le réservoir lors du remplissage. Il en résulte que la densité estimée de façon relativement complexe est non satisfaisante (notamment car cette estimation est basée sur une température prise égale à la température de l'hydrogène sortant au niveau du pistolet et n'utilisant par le facteur de compressibilité de gaz réels). 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. Les inventeurs ont mis au point une nouvelle méthode d'estimation simple et fiable de la température du gaz dans le réservoir lors du remplissage (cette température est généralement non mesurable pratiquement dans le réservoir). Les inventeurs proposent également d'utiliser cette valeur calculée de la température comme donnée de premier niveau qui est utilisée ensuite pour calculer ou estimer des données de second niveau tel que la densité par exemple. 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, pour l'étape de détermination de la température (T(ti)) courante du gaz dans le réservoir, ladite température (T(ti)) est exprimée et calculée en fonction uniquement des variables que sont la pression (P(ti)) courante dans le réservoir et de la quantité (m(ti)) courante de gaz dans le réservoir; l'expression de la température (T(ti)) courante en fonction de la pression (P(ti)) courante et de la quantité (m(ti)) courante de gaz dans le réservoir étant obtenue à partir de l'équation d'état des gaz réels dans le réservoir P(ti).V.105=Z.n.R.T(ti) dans laquelle P(ti) est la pression du gaz dans le réservoir à l'instant ti en bar , V le volume du réservoir en m3, R la constante des gaz parfait égale à 8,314 en J/(mol.K), T(ti) la température du gaz dans le réservoir à l'instant ti en Kelvin (K) et Z le facteur de compressibilité sans unité, ce facteur Z de compressibilité étant exprimé comme une fonction de la température T(ti) et de la pression P(ti) du gaz dans le réservoir selon une formule du premier degré : Z(ti)=(e.T(ti)+f).P(ti) + g, dans laquelle e, f et g sont des coefficients prédéterminés empiriquement avec e en bar'.K-', f en bar-', g sans unité. Par ailleurs, des modes de réalisation de l'invention peuvent comporter l'une ou plusieurs des caractéristiques suivantes : - le facteur Z de compressibilité est exprimé selon la formule Z(ti) =(e.T(ti)+f).P(ti) + g, avec g = 0,99651, e=-1,75724.10-6 et f=1,17735.10-3 - l'étape de mesure de la pression initiale (P(t0)) dans le réservoir avant le remplissage consiste à mesurer la pression (P) à l'entrée du réservoir, au niveau d'un conduit de remplissage du réservoir et, lorsque le réservoir possède au niveau de son orifice d'entrée/sortie un clapet anti-retour, à retrancher à cette valeur mesurée (P) le niveau de pression (Pcd) nécessaire pour ouvrir le clapet anti-retour (P(t0)=P-Pcd), - le procédé comporte une étape détermination de la température (T(t0)) initiale du gaz dans le réservoir avant remplissage, ladite température (T(t0)) initiale du gaz dans le réservoir étant choisie parmi : la température ambiante (Tamb) mesurée autour du réservoir, une mesure de température à l'intérieur du réservoir, une estimation de la température du gaz dans le réservoir à partir de la température ambiante et à partir de l'historique des valeurs de température (T(ti)) dans le réservoir, - l'étape de détermination de la quantité (Q(ti)) courante du gaz transférée dans le réservoir lors du remplissage utilise l'un au moins parmi : un débitmètre disposé en amont de l'entrée du réservoir pour mesurer la quantité de gaz transférer au réservoir lors du remplissage, un logique de calcul qui détermine cette quantité (Q(ti)) courante du gaz transférée dans le réservoir à partir de mesures de pression et de température en amont de l'entrée du réservoir, - l'étape de détermination de la quantité (m(t0)) initiale de gaz dans le réservoir avant remplissage utilise la formule de calcul de la masse m(t0) de gaz initiale dans le réservoir en kg: R . T(t0) . Ke . T(t0) + f ). P(t0) + g] P(t0) .105 . V.M m(t0) = dans laquelle les coefficients g= 0,99651 sans unité, e=-1,75724.10-6 bar'.K-' et f=1,17735.10-3 bar-'; M est la masse molaire du gaz en kg/mol, V le volume du réservoir en m3, la quantité (m(ti)) courante de gaz dans le réservoir lors du remplissage étant obtenue en ajoutant à cette quantité initiale m(t0) la quantité (Q(ti)) courante du gaz transférée dans le réservoir lors du remplissage : m(ti) = m(t0)+Q(ti), - la température T(ti) courante du gaz dans le réservoir (en K) est donnée par la formule : (f .P(ti) + g v + 4. e. P(ti)2.105. V. M R. m(ti) 2. e. P(ti) dans laquelle P(ti) est la pression courante du gaz dans le réservoir à l'instant ti en bar; g, e et f des coefficients avec e en bar'.K-', f en bar-', g sans unité donnés par g= 0,99651 sans unité, e= -1,75724.10-6 bar'.K-' et f= 1,17735.10-3 bar-', M la masse molaire du gaz en kg/mol, V le volume du réservoir 15 en m3 - le volume (V) du réservoir en m3 est connu ou calculé, - le procédé est mis en oeuvre par une station de remplissage de réservoirs de gaz hydrogène comprenant au moins une source d'hydrogène à haute pression, au moins une conduite de transfert de reliant sélectivement la source à 20 un réservoir, et une logique électronique de contrôle et de commande pilotant le transfert de gaz entre la source et le réservoir, la station de remplissage fournissant à la logique électronique de contrôle et de commande au moins l'un paramètre d'entrée parmi : la pression (P) mesurée dans la conduite en amont du réservoir, cette pression (P) mesurée dans la conduite à l'entrée du réservoir étant 25 assimilée à la pression P(ti) dans le réservoir, le débit massique (Q(ti)) de gaz courant dans la conduite de transfert, la température (T) courante du gaz dans la conduite de transfert, la durée (t) de remplissage, la pression nominale maximale (Pmax) du réservoir, le volume (V) du réservoir, la température ambiante (Tamb), - le procédé est mis en oeuvre par une station de remplissage de réservoirs 30 de gaz hydrogène comprenant au moins une source d'hydrogène à haute pression, au moins une conduite de transfert de reliant sélectivement la source à un réservoir, et une logique électronique de contrôle et de commande pilotant le transfert de gaz entre la source et le réservoir, la logique électronique étant programmée pour calculer au moins l'une des données de sortie suivantes : la 35 densité initiale (p(t0)) du gaz dans le réservoir avant le remplissage, la densité T(ti) = - (£ P(ti) + g )+ courante (p(ti)) du gaz dans le réservoir lors du remplissage, une densité cible (pf) déterminée dans le réservoir correspondant à un critère d'arrêt du remplissage, la température (T(i)) courante du gaz dans le réservoir, - le remplissage est contrôlé avec au moins l'un des paramètres de contrôle suivant: la pression (P(ti)) courante dans le réservoir lors du remplissage, la température (T(ti)) courante du gaz dans le réservoir, la densité courante (p(ti)) du gaz dans le réservoir lors du remplissage, - le remplissage est contrôlé avec la densité courante (p(ti)) du gaz dans le réservoir lors du remplissage, ladite densité courante (p(ti)) du gaz dans le réservoir étant calculée à partir de la température (T(i)) courante du gaz dans le réservoir calculée, à partir de la quantité m(ti)) courante de gaz dans le réservoir calculée et à partir du volume (V) du réservoir (3) selon la formule p(ti))=m(ti)/V ; avec p(ti) en kg/m3, m(ti) en kg, et V en m3, le remplissage étant contrôlé selon une courbe ou une droite de variation de densité (p(ti)) courante en fonction du temps prédéterminée, le remplissage étant interrompu lorsque la densité courante atteint une valeur cible déterminée - la vitesse de remplissage est contrôlé au moyen d'au moins un régulateur de pression de type proportionnel intégral, - le remplissage est contrôlé au moyen d'un premier régulateur de pression de type proportionnel intégral qui reçoit comme paramètre d'entrée la densité courante (p(ti)) du gaz dans le réservoir et régule le remplissage en fonction de cette densité courante (p(ti)), - le remplissage est contrôlée au moyen d'un second régulateur de pression de type proportionnel intégral qui reçoit comme paramètre d'entrée, la température (T(ti)) ou la pression (P(ti)) courante du gaz dans le réservoir et régule de le remplissage en fonction de cette température (T(ti)) courante, respectivement pression (P(ti) courante, - le procédé utilise un sélecteur électronique, le sélecteur réalisant une comparaison de la température courante (T(ti)) du gaz dans le réservoir avec une température maximum prédéfinie pour ce point de remplissage et, lorsque la température courante (T(ti)) du gaz dans le réservoir est supérieure à cette température maximum prédéfinie pour ce point de remplissage, le sélecteur active le second régulateur de pression de type proportionnel intégral (PI) et désactive le premier régulateur proportionnel intégral et, lorsque la température courante T(ti) du gaz dans le réservoir est inférieure à cette température maximum prédéfinie pour ce point de remplissage, le sélecteur active le premier régulateur de pression de type proportionnel intégral (PI) et désactive le second régulateur proportionnel intégral, - le procédé utilise un limiteur de débit de remplissage fourni au réservoir en dessous d'un seuil déterminé, par exemple de 60g/s. L'invention peut concerner également tout dispositif ou procédé alternatif comprenant toute combinaison des caractéristiques ci-dessus ou ci-dessous. 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 d'un diagramme logique illustrant une succession d'étapes possibles dans la mise en oeuvre d'un exemple de réalisation du procédé selon l'invention, - la figure 2 illustre de façon schématique et partielle un détail d'un exemple possible de structure de station de remplissage selon l'invention. L'invention peut concerner notamment le remplissage de réservoirs d'hydrogène de type adaptatif. Le terme adaptatif signifiant que le procédé de remplissage s'adapte à des paramètres du réservoir qui ne sont pas forcément totalement connus. Le procédé se déroule de préférence selon deux phases : une première phase d'estimation par calcul de plusieurs paramètres primaires, puis une seconde phase de remplissage contrôlée avec la température du gaz dans le réservoir et la quantité du gaz dans le réservoir (ou la densité qui donne la même information que la quantité lorsque le volume du réservoir est connu). Avant de débuter le remplissage le procédé mis en oeuvre par une station de remplissage détermine de préférence un ou des objets parmi : - le type d'équation d'état utilisé pour le gaz considéré, - la pression initiale, la température initiale dans le réservoir, - le type de réservoir (par exemple type III ou type IV, dimension, volume...), - la pression nominale d'utilisation du réservoir (« nominal working pressure » en anglais), - la quantité de gaz contenue initialement dans le réservoir, - la quantité de gaz à transférer dans le réservoir. Ces données sont estimées par la station ou communiquées à la station de préférence avant de calculer en temps réel la quantité (masse par exemple) et la température du gaz dans le réservoir lors du remplissage. L'utilisation d'une équation d'état simple et fiable basée sur des données empiriques est particulièrement avantageuse. Selon l'invention, l'équation d'état du gaz utilise le facteur de compressibilité Z (Z=P.M/(p.R.T) avec p la masse volumique, P la pression, T la température, M la masse molaire et R la constante des gaz parfaits). Le facteur de compressibilité Z, grandeur sans unité, est exprimé à partir de données empiriques ajustées à l'aide d'une fonction F qui dépend de la pression P(ti) courante (c'est-à-dire en temps réel) et de la température T(ti) courante : Z(ti) = F( P(ti), T(ti) ) Pour l'hydrogène, à partir des données l'Institut National des Standards et de Technologies « NIST », il a été mis au point une expression du facteur de compressibilité Z(ti) au temps ti selon la formule : Z = a(T(ti)).P(ti) + g Avec Z coefficient de compressibilité sans dimension, T(ti) exprimée en Kelvin (K), et g une constante sans dimension. De plus le coefficient a (en bar) est exprimé comme une fonction linéaire (du ter degré) de la température T(ti), c'est-à-dire : a=e.T(ti) + f Au final, Z=(e.T(ti)+ f).P(ti) + g (formule 1) dans lequel P(ti) est la pression (en bar) du gaz à l'instant ti, T(ti) la température du gaz dans le réservoir à l'instant ti en Kelvin (K) et avec e en bar'.K-', f en bar-' et g sans unité sont des coefficients prédéterminés empiriquement. De préférence g = 0,99651 sans unité, e=-1,75724.10-6 bar'.K-' et f=1.17735.10-3 bar-'. Les inventeurs ont constaté que la valeur du facteur Z ainsi calculée 20 présente une erreur maximale de 0,820/0 par rapport aux valeurs données par l'institut NIST pour l'hydrogène. L'équation d'état des gaz réels peut être donnée par : P.V=n.Z.R.T (P étant la pression en Pa, V le volume en m3, n le nombre de moles en mol, R la constante des gaz parfaits en J/(mol.K), T la température en 25 K et Z le facteur de compressibilité du gaz sans unité. n=m/M (m étant la masse en kg et M la masse molaire en kg/mol). Donc P.V.M=Z.m.R.T, en remplaçant la pression P en bar (au lieu des Pascal) : P(en bar).105 = P(en Pa) on obtient : 30 P(en bar).105. V.M= m.T.R.Z= m.T.R.((e.T+f).P(en bar)+g) On a donc pour chaque instant ti, P(ti).105.M.V=m(ti).T(ti).R.((e.T(ti)+f).P(ti)+g) (équation B) Avec P(ti) en bar, V(ti) en m3, T(ti) en K et R la constante des gaz parfaits égale à 8,314 en J/(mol.K). 35 Comme schématisé à la figure 2, le réservoir 3 à remplir est classiquement alimenté en gaz sous pression à partir d'au moins une source 1 et via une conduite 2 de remplissage munie d'au moins une vanne de contrôle 5. Une logique 4 électronique contrôle le remplissage à partir de mesure de pression P, température T, mesure de débit Q(ti) dans la conduite 2 de transfert. La logique 4 électronique reçoit également de préférence les paramètres suivants : le volume V du réservoir 3, la pression nominale de service ou pression maximale Pmax pour le réservoir, le temps t et la température ambiante . La pression initiale dans le réservoir T(t0) (avant remplissage) peut être approximée à la valeur de la pression P mesurée dans la conduite 2 de remplissage à l'entrée du réservoir 3. Dans le cas où le réservoir comprend un orifice pourvu d'un clapet anti-retour (« NRV »), on retranche le niveau de pression Pcd nécessaire pour ouvrir le clapet anti-retour à cette valeur mesurée P(t0)=PPcd (généralement de 0,5 à 2bar). La pression de travail maximum admissible du réservoir « MAWP » est généralement 1,25 fois la pression nominale de travail « NWP ». Cette pression nominale de travail peut être par exemple 350bar ou jusqu'à 750bar. Les connecteurs des conduites 2 de transfert peuvent être configurés en fonction d'une pression nominale de travail déterminée (chaque connecteur est 20 par exemple conformé pour une valeur de pression déterminée). La nature du réservoir 3 influence fortement les échanges thermiques lors du remplissage. De préférence, les caractéristiques thermodynamiques du réservoir 3 sont connues de la station de remplissage ou transmises automatiquement ou manuellement à la station (réservoir de type III ou de type IV 25 par exemple). La température initiale du gaz T(t0) dans le réservoir (avant remplissage), si elle ne peut être mesurée peut être approximée à la valeur de la température ambiante Tamb autour du réservoir, au niveau de la station de remplissage. La quantité Q(ti)-Q(t0) de gaz transférée dans le réservoir 3 au cours du 30 remplissage peut être mesurée via un débitmètre installé sur la conduite 2 de transfert. Un débitmètre est en effet préférable en terme de précision par rapport à un calcul de débit via un algorithme couplé à une vanne. Un débitmètre donne une bonne mesure à l'état stable. Puisque la quantité initiale de gaz Q(t0) et le volume d'hydrogène sont évalués pendant des 35 états transitoires, la performance du débitmètre pendant les phases transitoires doit de préférence être vérifiée pour correctement déterminer le volume et la quantité initiale d'hydrogène Q(t0). Un facteur de correction peut être prévu le cas échant pour corriger les données relevées lors de phases transitoires. La quantité (masse) de gaz initiale contenue dans le réservoir 3 (à l'instant t0 avant remplissage) Q(t0), si elle n'est pas mesurable directement, peut être obtenue à partir de l'équation B ci-dessus. C'est-à-dire que m(t0) _ (P(t0).105.V.M)/(R.T(t0).(e.T(t0)+f).P(t0)+g)) avec M (masse molaire) égale à 2*10-3 kg/mol pour l'hydrogène et la pression P(t0) en bar, e en bar'.K-', f en bar-', g sans unité et les autres paramètres étant en unités S.I . Le volume V du réservoir est supposé connu de la station par exemple par communication, à défaut une estimation du volume peut être calculée (cf. les remarques et les exemples de méthodes connues évoquées ci-dessus). A partir des formules ci-dessus il devient facile de calculer la donnée primaire constituée par la température T(ti) (et la quantité ou masse m(ti)) courante du gaz dans le réservoir 3. En effet, connaissant la quantité de gaz initiale m(t0) dans le réservoir, le volume V du réservoir 3 il est possible de calculer la masse de gaz présente dans le réservoir en temps réel (m(ti)=m(ti-1)+ la quantité ajoutée entre ti-1 et ti) et prédire ainsi la température T(ti). L'équation d'état des gaz réels peut notamment être utilisée pour déterminer la température T(ti) en temps réel. En partant de l'équation B, on a : P(ti).105.M.V = m(ti).R.T(ti).((e.T(ti)+f).P(ti)+g) _> (P(ti).105.M.V)/m(ti) = R.T(ti).((e.T(ti)+f).P(ti)+g) _> e.P(ti).T2(ti) + (f.P(ti)+g).T(ti) - (P(ti).105.M.V)/(R.m(ti)) = 0 La résolution de cette équation du second degré en T(ti) donne l'expression de la température courante T(ti) en fonction des paramètres connus. En effet, l'équation a x2 + b x + c = 0 donne : -b±.N/b2-4ac x= 2a dans notre cas on obtient : ( () ) ( () )2 4. e. P(ti)2.105. V. M - £Pti+ g + f.Pti+ g + 1 T(ti) _ R. m(ti) 2. e. P(ti) On obtient ainsi une expression simple et performante de la température T(ti) du gaz dans le réservoir à tout instant ti. La logique électronique 4 peut ainsi optimiser le remplissage tout en maintenant : - la température T(ti) du gaz dans le réservoir inférieure à un seuil de sécurité, par exemple 85°C, - la pression P(ti) dans le réservoir en dessous de la pression maximum de travail admissible (MAWP). Les paramètres d'entrée peuvent comprendre la pression P=P(ti) mesurée dans la conduite 2 de transfert en amont du réservoir 3, le débit massique Q(ti) de gaz dans la conduite 2 de transfert, la température T courante du gaz dans la conduite 2 de transfert, la durée t de remplissage, la pression nominale maximale Pmax du réservoir (3), les caractéristiques TK thermiques du réservoir 3, le volume V du réservoir 3, la température ambiante Tamb. Les données de sortie peuvent comprendre : la température T(i) courante du gaz dans le réservoir la densité initiale p(t0) du gaz dans le réservoir 3 avant le remplissage, la densité courante p(ti) du gaz dans le réservoir 3 lors du remplissage, une densité cible pf déterminée dans le réservoir 3 correspondant à un critère d'arrêt du remplissage et un contrôle de non franchissement d'une pression maximale de service (MAWP). Le critère de sécurité est par exemple contrôlé par un signal à un bit de la logique électronique. Le bit est par exemple : - égal à « 1 » quand la pression est en dessous du seuil maximum admissible et que la densité courante est inférieure à la densité cible pf et - égal à « 0 » dans le cas contraire. Le remplissage peut donc être contrôlé par la densité du gaz dans le réservoir ou un autre paramètre équivalent (la quantité de gaz par exemple). En déterminant une densité initiale p(t0), une densité cible finale pf et une durée de remplissage, il est possible de définir une rampe de remplissage paramétrée prédéterminée. Le contrôle simultané du remplissage via la densité instantanée p(ti) et via la température T(ti) instantanée peut être réalisé en utilisant deux organes de régulation, par exemple un premier régulateur de pression de type proportionnel intégral (PI) qui reçoit comme paramètre d'entrée la densité courante p(ti) du gaz dans le réservoir 3 et régule le remplissage en fonction de cette densité courante p(ti). En effet, les réservoirs d'hydrogène gazeux se comportent comme des capacités et un régulateur pour le remplissage est de préférence un régulateur de type PI (proportionnel intégral). Un second régulateur de pression par exemple également de type proportionnel intégral (PI) peut être prévu. Le second régulateur reçoit comme paramètre d'entrée, la température T(ti) ou la pression P(ti) courante du gaz dans le réservoir 3 et régule de le remplissage en fonction de cette température T(ti) courante, respectivement pression P(ti) courante. Ces deux régulateurs commandent l'ouverture et la fermeture d'une vanne 5 située sur la conduite 2 de transfert. Le premier régulateur reçoit en entrée - les densités de l'estimateur d'une rampe de densité prédéterminée et - la densité courante p(ti). En sortie, ce premier régulateur délivre un signal de commande fixant le pourcentage d'ouverture à appliquer à la vanne. Le second régulateur reçoit en entrée la température maximum admissible et la température courante T(ti). En sortie, le second régulateur délivre un signal de commande fixant le pourcentage d'ouverture à appliquer à la vanne. Les signaux de sortie des deux régulateurs peuvent être différents. Le premier régulateur peut définir un grand pourcentage d'ouverture tandis que le second régulateur peut donner par exemple un signal d'ouverture en pourcentage relativement plus petit. Un comparateur de température peut être prévu pour délivrer en sortie un signal de commande d'un sélecteur électronique. Le comparateur compare en entrée la température T(ti) du gaz dans le réservoir 3 en temps réel avec la pression maximale admissible. Si la température courante T(ti) excède la température maximale admissible le sélecteur commute en position de régulation de température. Si la température courante T(ti) est inférieure à la température maximale admissible le sélecteur commute en position de régulation de densité. Un limiteur peut être prévu pour limiter le débit de gaz maximum fourni au réservoir (par exemple 60g/s d'hydrogène conformément au standard SAEJ2601). La figure 1 illustre un exemple d'étapes pouvant être mises en oeuvre par une station de remplissage selon l'invention. Dans deux premières étapes, la température initiale T(t0) et la pression initiale P(t0) peuvent être mesurées/déterminées. Au cours d'une étape ultérieure, la station peut déterminer si une communication statique est possible, c'est-à-dire si des paramètres sont susceptibles d'être communiqués à la station par l'entité (véhicule cadre de bouteilles) comprenant le réservoir. Si c'est non (N), à l'étape 103 des paramètres du réservoir TK peuvent le cas échéant être mesurés ou communiqués et à l'étape 104, le volume V et la densité initiale p(t0) sont calculés ou estimés. Si c'est oui (Y), des paramètres TK et notamment le volume V du réservoir sont communiqués à la station (étape 105). Après cette communication la densité initiale p(t0) dans ce réservoir est calculée ou estimée (comme à l'étape 104). Qu'il y ait ou non communication, l'étape suivante 107 peut être une étape de calcul de la densité cible Pf finale pour le réservoir. La durée F du remplissage peut être choisie ou calculée à l'étape 108 pour permettre de définir une rampe RA d'augmentation de densité. Lors du remplissage, la densité p(ti) courante peut être calculée en temps réel (étape 109). Comme décrit ci-dessus, cette densité p(ti) courante peut être calculée à partir de l'expression du facteur de compressibilité Z comme décrit ci-dessus (étape 110). Ce facteur Z de compressibilité est ensuite utilisé avec la pression courante P(ti), la quantité courante m(ti) de gaz dans le réservoir pour calculer la température courante (étapes 111 et 112). La densité p(ti) courante est ensuite calculée à partir de cette valeur de température T(ti) (via le volume et la quantité courante m(ti)). La station peut réguler le remplissage à l'aide de la densité p(ti) courante en utilisant les conditions d'interruption appropriée (la température doit rester inférieure à la température Tmax et le taux de remplissage SOC doit être inférieur à 1000/0 (étape 113)). Lorsque les conditions d'interruption sont remplies, le remplissage est interrompu (étape 114). De préférence, le remplissage est du type à communication statique. Ce type de communication permet en effet de réduire les incertitudes dans l'estimation de paramètres (volume...). La communication de valeurs réelles à la station de remplissage permet d'améliorer l'optimisation du remplissage. De préférence, les données suivantes sont communiquées à la station : - la pression nominale de service du réservoir (NWP), - le volume V du réservoir, le nombre de réservoirs, le ratio de la surface du réservoir sur son volume, les propriétés thermodynamiques du réservoir, son type (III ou IV)... Ces données peuvent être transférées sans fil, via un code-barres, un dispositif d'identification émetteur du type à radiofréquence (RFID)... Pour le remplissage d'un réservoir de bus ayant une capacité de stockage de 20 à 40kg d'hydrogène, la fin du remplissage est réalisée de préférence avec un compresseur. En effet, dans ce cas le différentiel de pression entre la conduite de transfert et le réservoir est constant et mesurable facilement. Une mesure de température du gaz au niveau de la conduite de transfert peut permettre de calculer l'enthalpie du gaz entrant dans le réservoir. REVENDICATIONS1. Procédé de remplissage d'un réservoir avec du gaz sous pression, notamment de l'hydrogène gazeux, pour atteindre un taux de remplissage cible (Pr) prédéterminé, le remplissage étant interrompu lorsqu'une grandeur physique (p(ti)) mesurée ou estimée du gaz dans le réservoir correspond au taux de remplissage cible (pf) ou lorsque la température (T(ti)) dans le réservoir atteint un seuil maximum (Tmax) déterminée, le procédé comprenant : - une étape de mesure de la pression initiale (P(t0)) dans le réservoir avant le remplissage, - une étape de détermination de la quantité (m(t0)) initiale de gaz dans le réservoir avant remplissage, - une étape mesure de la pression (P(ti)) courante dans le réservoir lors du remplissage, - une étape de détermination de la quantité (Q(ti)) courante de gaz transférée dans le réservoir lors du remplissage, - une étape de calcul de la quantité m(ti)) courante de gaz dans le réservoir lors du remplissage, - une étape de détermination de la température (T(ti)) courante du gaz dans le réservoir lors du remplissage, le procédé étant caractérisé en ce que, pour l'étape de détermination de la température (T(ti)) courante du gaz dans le réservoir, ladite température (T(ti)) est exprimée et calculée en fonction uniquement des variables que sont la pression (P(ti)) courante dans le réservoir et de la quantité (m(ti)) courante de gaz dans le réservoir; l'expression de la température (T(ti)) courante en fonction de la pression (P(ti)) courante et de la quantité (m(ti)) courante de gaz dans le réservoir étant obtenue à partir de l'équation d'état des gaz réels dans le réservoir P(ti).V.105=Z.n.R.T(ti) dans laquelle P(ti) est la pression du gaz dans le réservoir à l'instant ti en bar , V le volume du réservoir en m3, R la constante des gaz parfait égale à 8,314 en J/(mol.K), T(ti) la température du gaz dans le réservoir à l'instant ti en Kelvin (K) et Z le facteur de compressibilité sans unité, ce facteur Z de compressibilité étant exprimé comme une fonction de la température T(ti) et de la pression P(ti) du gaz dans le réservoir selon une formule du premier degré : Z(ti)=(e.T(ti)+f).P(ti) + g, dans laquelle e, f et g sont des coefficients prédéterminés empiriquement avec e en bar'.K-', f en bar-', g sans unité. 2. Procédé selon la revendication 1, caractérisé en ce que le facteur Z de compressibilité est exprimé selon la formule Z(ti) =(e.T(ti)+f).P(ti) + g, avec g = 0,99651, e=-1,75724.10-6 et f=1,17735.10-3. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que l'étape de mesure de la pression initiale (P(t0)) dans le réservoir (3) avant le remplissage consiste à mesurer la pression (P) à l'entrée du réservoir, au niveau d'un conduit (2) de remplissage du réservoir (3) et, lorsque le réservoir (3) possède au niveau de son orifice d'entrée/sortie un clapet antiretour, à retrancher à cette valeur mesurée (P) le niveau de pression (Pcd) nécessaire pour ouvrir le clapet anti-retour (P(t0)=P-Pcd). 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce qu'il comporte une étape détermination de la température (T(t0)) initiale du gaz dans le réservoir (3) avant remplissage, ladite température (T(t0)) initiale du gaz dans le réservoir (3) étant choisie parmi : la température ambiante (Tamb) mesurée autour du réservoir (3), une mesure de température à l'intérieur du réservoir, une estimation de la température du gaz dans le réservoir à partir de la température ambiante et à partir de l'historique des valeurs de température (T(ti)) dans le réservoir. 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que l'étape de détermination de la quantité (Q(ti)) courante du gaz transférée dans le réservoir lors du remplissage utilise l'un au moins parmi : un débitmètre disposé en amont de l'entrée du réservoir pour mesurer la quantité de gaz transférer au réservoir lors du remplissage, un logique (4) de calcul qui détermine cette quantité (Q(ti)) courante du gaz transférée dans le réservoir à partir de mesures de pression et de température en amont de l'entrée du réservoir (3). 6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce que l'étape de détermination de la quantité (m(t0)) initiale de gaz dans le réservoir avant remplissage utilise la formule de calcul de la masse m(t0) de gaz initiale dans le réservoir en kg: R . T(t0) . Ke . T(t0) + f ). P(t0) + g] dans laquelle les coefficients g= 0,99651 sans unité, e=-1,75724.10-6 bar'.K-' et f=1,17735.10-3 bar-'; M est la masse molaire du gaz en kg/mol, V le volume du réservoir en m3, la quantité (m(ti)) courante de gaz dans le réservoir lors du remplissage étant obtenue en ajoutant à cette P(t0) .105 . V .M m(t0) =quantité initiale m(t0) la quantité (Q(ti)) courante du gaz transférée dans le réservoir lors du remplissage : m(ti) = m(t0)+Q(ti). 7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que la température T(ti) courante du gaz dans le réservoir 5 (en K) est donnée par la formule : ( () ) ( () )a 4. e. P(ti)2.105. V. M - £Pti+ g + f.Pti+ g + T(ti) = 1 R. m(ti) 2. e. P(ti) dans laquelle P(ti) est la pression courante du gaz dans le réservoir à l'instant ti en bar; g, e et f des coefficients avec e en bar'.K-', f en bar-', g sans unité donnés par g= 0,99651 sans unité, e= -1,75724.10-6 bar'.K-' et f= 10 1,17735.10-3 bar-', M la masse molaire du gaz en kg/mol, V le volume du réservoir en m3. 8. Procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce que le volume (V) du réservoir en m3 est connu ou calculé. 9. Procédé selon l'une quelconque des revendications 1 à 8, 15 caractérisé en ce qu'il est mis en oeuvre par une station de remplissage de réservoirs de gaz hydrogène comprenant au moins une source (1) d'hydrogène à haute pression, au moins une conduite (2) de transfert de reliant sélectivement la source (1) à un réservoir (3), et une logique (4) électronique de contrôle et de commande pilotant le transfert de gaz entre la 20 source (1) et le réservoir (3), caractérisé en ce que la station de remplissage fournit à logique (4) électronique de contrôle et de commande au moins l'un paramètre d'entrée parmi : la pression (P) mesurée dans la conduite (2) en amont du réservoir (3), cette pression (P) mesurée dans la conduite à l'entrée du réservoir (3) étant assimilée à la pression P(ti) dans le réservoir (3), le 25 débit massique (Q(ti)) de gaz courant dans la conduite (2) de transfert, la température (T) courante du gaz dans la conduite (2) de transfert, la durée (t) de remplissage, la pression nominale maximale (Pmax) du réservoir (3), le volume (V) du réservoir (3), la température ambiante (Tamb). 10. Procédé selon l'une quelconque des revendications 1 à 9, 30 caractérisé en ce qu'il est mis en oeuvre par une station de remplissage de réservoirs de gaz hydrogène comprenant au moins une source (1) d'hydrogène à haute pression, au moins une conduite (2) de transfert de reliant sélectivement la source (1) à un réservoir (3), et une logique (4) électronique de contrôle et de commande pilotant le transfert de gaz entre la 35 source (1) et le réservoir (3), caractérisé en ce que la logique (4) électronique est programmée pour calculer au moins l'une des données de sortie suivantes : la densité initiale (p(t0)) du gaz dans le réservoir (3) avant le remplissage, la densité courante (p(ti)) du gaz dans le réservoir (3) lors du remplissage, une densité cible (pf) déterminée dans le réservoir (3) correspondant à un critère d'arrêt du remplissage, la température (T(i)) courante du gaz dans le réservoir (3). 11. Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce que le remplissage est contrôlé avec au moins l'un des paramètres de contrôle suivant: la pression (P(ti)) courante dans le réservoir (3) lors du remplissage, la température (T(ti)) courante du gaz dans le réservoir, la densité courante (p(ti)) du gaz dans le réservoir (3) lors du remplissage. 12. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce que le remplissage est contrôlé avec la densité courante (p(ti)) du gaz dans le réservoir (3) lors du remplissage, ladite densité courante (p(ti)) du gaz dans le réservoir (3) étant calculée à partir de la température (T(i)) courante du gaz dans le réservoir calculée, à partir de la quantité m(ti)) courante de gaz dans le réservoir calculée et à partir du volume (V) du réservoir (3) selon la formule p(ti))=m(ti)/V ; avec p(ti) en kg/m3, m(ti) en kg, et V en m3, le remplissage étant contrôlé selon une courbe ou une droite de variation de densité (p(ti)) courante en fonction du temps prédéterminée, le remplissage étant interrompu lorsque la densité courante atteint une valeur cible déterminée (pf). 13. Procédé selon l'une quelconque des revendications 1 à 12, caractérisé en ce que la vitesse de remplissage est contrôlé au moyen d'au moins un régulateur de pression de type proportionnel intégral (PI).