Brevet FR2922991A1 — publié le 20090501

La présente invention concerne un procédé d'estimation de paramètres caractéristiques d'un réservoir cryogénique et notamment de paramètres géométriques du réservoir. L'invention permet notamment l'amélioration de la mesure de niveau dans les réservoirs cryogéniques afin d'améliorer l'efficacité de la chaîne logistique d'approvisionnement en liquide de ces réservoirs. Les réservoirs concernés comprennent un réservoir interne de stockage du fluide (ou enveloppe interne) disposé dans un réservoir externe (ou enveloppe externe). Ces deux enveloppes sont séparées par une couche d'isolant. Les réservoirs stockent des liquides cryogéniques tels que l'oxygène, l'argon, l'azote avec des capacités de 100 litres à 100000 litres par exemple. Les pressions de stockage peuvent être comprises entre 3 bar et 35 bar. Les paramètres géométriques d'un réservoir (interne) sont nécessaires notamment pour notamment estimer le niveau de liquide et la quantité livrable dans le réservoir. Parmi les paramètres utiles, on peut citer notamment (pour le cas d'un réservoir cylindrique d'extrémités elliptiques): le rayon R , la hauteur totale du réservoir (htot ), la hauteur F de la partie elliptique (fond), la hauteur maximale de liquide Hmax . Pour de nombreux réservoirs cryogéniques ces paramètres sont inconnus ou ne sont identifiables qu'au prix d'efforts importants. 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 qu'il comprend une étape de calcul d'un premier paramètre géométrique (R) du réservoir en fonction : - d'une masse de liquide livrée (mdelivered) déterminée au cours d'un remplissage, - de la différence entre les différentiels de pression (DPmes) mesurés avant et après remplissage (DPmes after ùDPmes before ), chaque différentiel de pression mesurant le différentiel de pression (entre les parties haute et basse du réservoir, -des masses volumiques du gaz et du liquide (pg , pl) dans le réservoir. Par ailleurs, des modes de réalisation de l'invention peuvent comporter l'une ou plusieurs des caractéristiques suivantes : - le réservoir comporte une portion cylindrique et au moins une extrémité ayant une portion elliptique de hauteur (F) déterminée et en ce que le premier paramètre géométrique et le rayon (R) du cylindre, - le rayon (R) est calculé selon une équation du type : R = mdelivered g [1ù (MAvo)] 5 (DPmes _ after ù DPmes _before dans laquelle g est l'accélération de la pesanteur terrestre et MAVO est un coefficient de masse volumique fonction de la masse volumique du liquide et du 10 gaz dans le réservoir et éventuellement dans des conduits de mesure de pression lorsque la mesure du différentiel de pression est réalisée par au moins un capteur de pression déporté relié aux parties haute et basse du réservoir via des tuyaux de mesure respectifs, - pour l'estimation du coefficient de masse volumique (MAVO) le procédé utilise l'une au moins des hypothèses suivantes : - la masse volumique du liquide dans le réservoir pl est considérée comme étant égale à la moyenne des masses volumiques de liquide avant pl before et après pl after remplissage, - la masse volumique du liquide avant remplissage pl before est égale à la masse volumique à l'équilibre à la pression du réservoir, - la masse volumique du liquide après remplissage pl after est égale à la moyenne entre d'une part la masse volumique du liquide à l'équilibre à la pression du réservoir avant remplissage pondérée par la fraction du volume occupé par ce liquide et, d'autre part, la masse volumique du liquide dans le camion considérée à l'équilibre à la pression du camion pondérée par la fraction du volume disponible dans le réservoir avant remplissage, - le volume du liquide dans le réservoir avant remplissage est considéré égal à une fraction connue (par exemple 30%) du volume maximal de liquide dans le réservoir, -la masse volumique du gaz dans le réservoir pg est calculée à la pression du réservoir Ptank et pour une température augmentée par rapport à la température d'équilibre à la pression du réservoir (par exemple augmentée de 20 K), - le différentiel de pression DPmes mesuré est corrigé en tenant compte d'une valeur de différence de pression additionnelle (DPpipe) créée par le gaz présent dans les tuyaux de mesure en cas de mesure déportée, -le volume du liquide V1 dans un réservoir cylindrique de rayon R et ayant une extrémité elliptique de hauteur F est donné par la relation Si hi F = Vi=itR2[hi-3] R2 Sihi maximale de liquide Hmax, la hauteur hg before de gaz avant remplissage en étant déduite comme le complément : hl _before = FS H. hg _before = htot-FS H. - le procédé comporte une étape de calcul d'un quatrième paramètre géométrique constitué par la hauteur de liquide maximale Hmax, celle-ci étant est déduite la hauteur totale du réservoir calculée à partir de l'équation suivante : F Vmax C 2F 1 Hmax =ù+ htot - -J 3 Vtot 3 - le procédé comporte une étape de calcul d'un cinquième paramètre géométrique constitué par la perte thermique du réservoir, ladite perte thermique exprimée en pourcentage d'oxygène (%02) d'oxygène perdu par jour étant approximée selon une relation du type : % 02 perte par jour = pi Vtopz dans lequel Vtot est le volume total du réservoir pl est un coefficient de l'ordre de 0,6 et de préférence égal à 0,65273 et p2 un coefficient de l'ordre de 0,3 et de préférence égal à 0,37149, - le procédé comporte une étape de collecte d'une pluralité de valeurs de masse de liquide livrée (mdelivered) déterminées au cours respectivement de plusieurs remplissages du réservoir, une étape de calcul, pour chaque valeur de masse de liquide livrée (mdelivered) du rapport r de remplissage r= mdelivered = constante DPmes_after ùD mes_before et en ce que ne sont utilisées dans le procédé d'estimation que les valeurs de masse de liquide livrée (mdelivered) pour lesquelles la valeur absolue du rapport25 r de remplissage ne présente pas une différence supérieure à un seuil fixé vis-à-vis d'une valeur constante de référence déterminée. - la valeur constante de référence déterminée pour rapport r de remplissage est constituée de la moyenne de ce rapport de remplissage r calculée pour plusieurs remplissages. 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 illustrant un premier exemple de réservoir cryogénique de mise en oeuvre de l'invention (conduits externes aux parois du réservoir), - la figure 2 représente une vue schématique illustrant un second exemple de réservoir cryogénique de mise en oeuvre de l'invention (conduits internes aux parois du réservoir). La méthode qui va être décrite ci-après peut être mise en oeuvre par un ordinateur d'un système de contrôle de réservoir (local ou déporté). Cette méthode comprend une mesure de pression et de différence de pression DPmes et peut comprendre une transmission à distance de données. Les pressions sont mesurées via des conduits 11, 12 qui peuvent être dans l'espace inter parois du réservoir (figure 2) ou à l'extérieur 11 (figure 1). Le réservoir 1 peut comprendre un dispositif de pressurisation tel qu'un réchauffeur 3 de vaporisation apte à prélever du liquide pour le vaporiser et le réinjecter dans le réservoir. Ce réchauffeur 3 régule classiquement la pression au sein du réservoir 1. Par souci de simplification, le réservoir intérieur qui stocke le fluide sera désigné ci-après uniquement par le terme � réservoir �. Le liquide apporté par camion de livraison lors des remplissages peut également être considéré à l'état d'équilibre (gamme de température de 10 K autour de l'équilibre, par exemple 77,2 à 87,9 K pour de l'azote). La pression du liquide dans le camion de livraison est choisie, selon la pression du réservoir, entre 1 et 2 bar. Le liquide est introduit dans le réservoir par pompage. Entre deux remplissages le réservoir 1 est soumis aux phénomènes suivants : - la quantité de liquide décroit (consommation de l'utilisateur) et la baisse de pression correspondante est corrigée par le réchauffeur 3, - de la chaleur entre dans le réservoir à travers les parois du réservoir et le réchauffeur (conduction, radiation). Après un certain temps d'équilibre, du liquide se vaporise dans le réservoir ce qui contribue à une perte de liquide. De plus, la masse volumique du liquide diminue lorsque le liquide se réchauffe et par conséquent le niveau de liquide est plus élevé que s'il avait conservé sa température de livraison. Selon une particularité avantageuse, on considère des températures propres au gaz et au liquide dans le réservoir mais sans que ces températures ne soient fonction de la localisation dans le réservoir. C'est-à-dire que dans la suite les températures de gaz Tg et de liquide Tl sont des températures moyennes. Le niveau de liquide estimé est basé sur le différentiel de pression mesuré DPmes entre les extrémités basse et haute du réservoir. Selon la méthode actuelle, la hauteur de liquide calculée h11 est calculée selon la formule (équation 1) : 105 DPmes Pli g Avec pl, une valeur de masse volumique de liquide de calibration constante (mais modifiable par un opérateur) ; g étant l'accélération de la pesanteur terrestre. Du fait que le réservoir n'est pas un cylindre géométriquement parfait (extrémités elliptiques, cf. figures 1 et 2), le volume VI de liquide utilise deux équations suivant que le liquide a un niveau situé au-dessous ou au-dessus de la partie F elliptique (équations 2) : Si hllest au-dessus de la zone F elliptique F alors VI =n R2 hti --3 h11 = sinon= =2 ùTc FR2ù7c (Fûhm 2 R2 ù R 2 (Fùh11)2 3F' R étant le rayon du réservoir (au niveau de sa portion cylindrique). La masse de liquide contenue dans le réservoir ml est déduite en utilisant la masse volumique du liquide pl, (équation 3) : m1 = Pu V1 La masse de liquide mi peut être exprimée en fonction du différentiel de pression mesuré DPmes . De préférence, selon une caractéristique avantageuse possible de l'invention, le niveau de liquide hi calculé est corrigé en tenant compte d'une valeur de différence de pression additionnelle DPpipe créée par le gaz présent dans les tuyaux de mesure 11, 12, aussi bien lorsque les tuyaux 11 sont situés dans le réservoir (figure 2) que hors du réservoir (figure 1). C'est-à-dire que les capteurs 4 de pression sont déportés et � lisent � des pressions influencées par le fluide dans les conduits 11, 12 les reliant aux parties hautes et basses du réservoir. Le différentiel de pression DPmes mesuré de façon déportée entre les parties hautes et basses du réservoir étant relié au différentiel de pression dit � réel � DPreal entre les parties hautes et basses du réservoir selon la formule : DPmes = DPreal ù DPpipe Cas où les conduits sont à l'extérieur de la paroi du réservoir (fiqure 1) : DPwaii est le différentiel de pression entre les deux extrémités du conduit vertical traversant l'inter-paroi (en haut ou en bas). DPtot length est la différence de pression due à la pression de gaz dans la partie de conduite 11 reliant le point le plus haut à l'organe 4 de mesure (capteur) déporté. DPamb est la différence de pression due à la pression de gaz dans la partie de conduite 11 reliant le point le plus bas à l'organe 4 de mesure (capteur) déporté. Le différentiel de pression DPwaii entre les deux extrémités du conduit vertical traversant l'inter-paroi (en haut ou en bas) peut être considéré sensiblement identique aux parties haute et basse (uniquement le fait de gaz dans le conduit). En considérant la forme du conduit inférieur 12 dans l'inter-paroi : le conduit s'étend près de l'enveloppe extérieure pour � capter � les calories extérieures au réservoir et vaporiser complètement le fluide dans le conduit 12 de mesure. Entre les extrémités supérieure et inférieure de cette portion, la pression est sensiblement la même (différentiel de 0,5 bar au maximum). Cas où les conduits sont dans l'inter-paroi (fiqure 2) : DPside gaz est la différence de pression dans la partie du tuyau relié à la partie supérieure du réservoir et se trouvant en face du gaz dans le réservoir (contient du gaz), DPside hq est la différence de pression dans la partie du tuyau supérieur se trouvant en face du liquide dans le réservoir (contient du liquide). La masse totale mt0 de fluide au sein du réservoir (liquide et gaz) peut être exprimée en fonction de données parmi notamment : - du différentiel de pression DPmes mesuré entre les parties haute et basse du réservoir, -de la masse volumique du liquide dans le réservoir pl, - de la masse volumique du gaz dans le réservoir pg , - de la masse volumique du gaz dans le tuyau 11 en face de la partie gaz dans le réservoir mesurant la pression dans la partie haute du réservoir pside_gas , - de la masse volumique du gaz dans le tuyau 11 en face de la partie liquide dans le réservoir mesurant la pression dans la partie haute du réservoir pside_hqu~d , - de l'accélération de la pesanteur terrestre g, - du rayon R du réservoir, - de la hauteur de la partie elliptique F , - de la hauteur totale du réservoir ht0 - La masse totale mt0 de fluide au sein du réservoir peut être exprimée sous la forme d'une équation du type ci-dessous qui sera justifiée plus en détail ci-après : 1 Os DP R2 mes P1û Pg û F (Pl +P)+ Pside gas P1û Pside liquid Pg m tot - -7c g tot g P1 ù Pg + Pside gas ù Pside liquid / 3 P1 ù Pg + Pside gas ù Pside liquid (équation 101) Cette équation peut être appliquée aussi bien avant qu'après un remplissage du réservoir. On peut faire l'hypothèse que les masses volumiques du gaz et du liquide sont constantes avant remplissage et après remplissage et égales à leurs valeurs moyennes. Ainsi, en appliquant cette formule 101 aux états avant mtot(avant) et après remplissage mtot(après) on peut exprimer la masse de liquide livrée mdelivered = mtot(après) - mtot(avant) et éliminer ainsi les inconnues géométriques F et htot.. De cette façon, la masse livrée (connue en principe lors de la livraison) peut être exprimée en fonction uniquement des différentiels de pression DPmes mesurés avant et après remplissage DPmes before et DPmes after et du rayon R (inconnu). Ainsi, le rayon R peut être exprimé en fonction uniquement de la masse livrée, des masses volumiques de gaz et de liquide et des différentiels de pression sous la forme suivante (équation 102): P sida liquid ù P sida gas Pl û Pg 10s 7[ DPmes after ù DPmes before Il faut cependant faire remarquer que, en pratique, les bulletins de livraison des opérateurs ne sont pas toujours fiables en ce qui concerne la masse de liquide réellement livrée mdelivered. En effet, des erreurs peuvent être dues à une transcription incorrecte par l'opérateur et/ou à des pertes de liquide lors de la manipulation. Ainsi, une imprécision de la masse livrée mdelivered peut fausser l'estimation du rayon R du réservoir. Pour pallier ce problème, une méthode 25 particulière décrite ci-dessous peut être utilisée pour ne retenir que les valeurs de masse livrée mdelivered fiables. Si l'on considère que les masses volumiques de gaz pg et du liquide pl sont constantes avant et après remplissage et égales à leurs valeurs moyennes, le ( delivered g 1 ù R rapport r entre d'une part la masse livrée et, d'autre part, la différence entre les différentiels de pression DPmes mesurés avant et après remplissage (DPmes after. ùDPmes before) est constant. C'est-à-dire que (équation 103) : 1) calculer ce rapport pour chaque remplissage 2) calculer la différence relative entre ce rapport r et la valeur moyenne 10 rmean de ce rapport 3) sélectionner les bulletins de livraison (masses livrées) qui, en valeur absolue, ne s'écartent pas trop de la moyenne rmean (par exemple pas plus de 10% de divergence) De préférence seules ces données sont utilisées pour le calcul du rayon 15 (équation 102). La masse volumique du liquide dans le réservoir pl est considérée comme étant égale à la moyenne des masses volumiques de liquide avant pl before et après Pi after remplissage (équation 104): p i before + p i after pi = - 2 - 20 On fait l'hypothèse également que la masse volumique du liquide avant remplissage pi before est égale à la masse volumique à l'équilibre à la pression du réservoir (équation 105) : pl_before = Pl_eq_ tank On fait l'hypothèse que la masse volumique du liquide après remplissage 25 pl after est égale à la moyenne entre d'une part la masse volumique du liquide à l'équilibre à la pression du réservoir avant remplissage pondérée par la fraction du volume occupé par ce liquide et, d'autre part, la masse volumique du liquide dans r_ mdelivered = constante DPmes _ after ù DPmes _ before Ainsi, pour sélectionner les masses livrées fiables, pour chaque bulletin de livraison la méthode peut : le camion considérée à l'équilibre à la pression du camion pondérée par la fraction du volume disponible dans le réservoir avant remplissage. Ceci conduit à (équation 106) : Pl_after = 0.7 * Pl_eq_truck + 0.3 * Pl_eq_tank La masse volumique du gaz dans le réservoir pg est calculée à la pression du réservoir Ptank et pour une température supérieure de 20 K à la température d'équilibre. En effet, le gaz dans le réservoir se trouve réchauffé après remplissage par rapport à sa température d'équilibre (environ 40 K au dessus de l'équilibre) juste avant le remplissage suivant. La valeur de 20 K est une moyenne qui peut être choisie avantageusement. On arrive alors à l'expression suivante (équation 107 qui donne la masse volumique du gaz comme étant calculée en fonction de plusieurs paramètre) : p g = p g (Tg = Teq tank + 20K, Ptank ) La masse volumique pslde_gas du gaz dans la conduite 11 reliant la partie haute du réservoir et située dans le réservoir (dans l'inter-parois), c'est-à-dire le gaz situé dans le conduit en face du gaz dans le réservoir peut être calculée à la pression du réservoir et à une température Tgg donnée par la formule suivante (équation 108) : T =T + d -Pipe `Tamb T) si la conduite est à l'extérieur gg g w _ length g Tgg =Tamb si la conduite est à l'intérieur Avec : - d _ pipe= la distance (espacement) entre le conduit supérieur 11 et la paroi du réservoir intérieur, - w_length = l'épaisseur de l'isolation du réservoir interne et - Tamb = la température ambiante autour du réservoir. Pour le cas d'un conduit 11 situé à l'intérieur (dans l'espace inter-parois) on peut considérer un profil de température linéaire dans l'épaisseur de l'isolation. Dans le cas d'un conduit 11 extérieur, la température du gaz dans le conduit Il est considérée égale à celle de la température ambiante. La masse volumique Pside_/iquid du gaz dans le conduit 11 et faisant face à la phase liquide dans le réservoir est calculée à la pression courante du réservoir et à une température Tg/ selon la relation ci-dessous (équation 109) : Tgt = TI + d _pipe (Tamb ù Tl) si le conduit est dans l'inter - parois w _ length Tgl = Tamb si le conduit est à l' extérieur des parois De la même façon, pour une conduite 11 dans le réservoir et faisant face à la phase liquide contenue dans le réservoir, on considère un profil de température linéaire entre le liquide situé dans le réservoir et la température ambiante à l'extérieur autour du réservoir. Le volume du liquide VI dans un réservoir cylindrique de rayon R et ayant une extrémité elliptique de hauteur F est donné par la relation (équation 110) : Si hi F VI =7c R2 hiù3 2 Sihi 2F 3 i 13 Lorsque hl est supérieure ou égale à F (la plupart des cas), en utilisant les équations 110 et 113, la formule donnant la masse peut se simplifier pour donner (équation 114) : 2 10 5 DPmes F mtot = 'Tc R ù 3 (Pl + P g )+(Pside_liquid ù Pside_gas ) hi + Pside_gas htot g La hauteur de liquide dans le réservoir est donc reliée à la mesure de pression différentielle DPmes selon la relation suivante (équation 115) : 105 DPmes l û (Pg û Pride_gas)htot g Pl ù Pg +P side_gas ù Pside_liquid Ainsi, les équations 14 et 15 conduisent à l'équation 101 présentée ci-avant. 10 La rayon R du réservoir peut ainsi être calculé et estimé pour chaque valeur plausible de masse livrée lors d'un remplissage. Le rayon moyen peut ainsi être calculé sur la base de ces multiples calculs. Il s'agit du premier paramètre déterminé à partir de la seule mesure du différentiel de pression DP de la masse de liquide livrée et de quelques approximations concernant les masses 15 volumiques. Calcul de la hauteur F de la partie elliptique d'extrémité. Cette hauteur de fond F se déduit directement du rayon R estimé, le rapport entre ces deux paramètres géométriques est considéré sensiblement constant pour tous les fabricants de réservoirs. Ce second paramètre géométrique peut être 20 déduit (cf. les équations 116 ci-dessous pour deux exemples de constructeurs). R =1.9 pour des réservoirs de la société "Cryolor" F ù = 2 pour des réservoirs de la société "Chart" Dans le cas où le fabricant est inconnu, l'approximation 1,95 peut être utilisée par exemple. Estimation de la hauteur maximale de liquide et de la hauteur totale H,( et hl = R F 25 ht0t ). Pour la plupart des réservoirs cryogéniques le ratio entre le volume de liquide maximal Vmax et le volume total du liquide Vtot est constant et dépend uniquement du niveau de pression dans le réservoir. Ce ratio est de 0,95 pour les réservoirs avec des pressions basses et moyennes (comprises entrel et 15 bar) et 5 de 0,90 pour des pressions hautes (supérieures à 15 bar), cf. équation 117) : Vmax = 0.95 pour des réservoirs à basse et moyenne pression Vmax =7c R 2 [Hmax 3 On obtient (équation 119) F Vmax [hù__] 2F 10 Hmax =ù+ 3 Vtot 3 Ainsi, une fois que la hauteur totale maximale de liquide Hmax peut être déduite en utilisant cette dernière équation (en connaissant Vmax) pour ensuite estimer F . Vtot La hauteur totale du réservoir est déterminée et estimée à partir du 15 différentiel de pression mesuré juste après un remplissage DPmes after en faisant l'hypothèse que le remplissage est total. Dans ce cas, le différentiel de pression mesuré peut être exprimé sous la forme (équation 120) : DPmes _ after = P l _ after Hmax g + P g _ after (htot ù Hmax g (Pgg_before hg_before + Pgl_before hl_before) g 20 Dans cette équation 120, les masses volumiques de gaz dans le conduit 11 de mesure Pside gas_before et pside liqu id _before sont calculées à la pression du réservoir après remplissage mais avec des valeurs de températures de gaz dans les conduits 11 avant remplissage (équation 121) : V tot Sachant que (équation 118) : R2rh 2F1 3 Vtot =7c tot V tot Vmax = 0.90 pour des réservoirs à haute pression htot est déterminée, la hauteur P Bide gas _before = P g (Tgg _ before, Ptank _ after ) Pside liquid_before = P g (Tgl_before, Ptank _after) Ceci est le résultat de l'inertie thermique de l'isolation du réservoir située dans l'espace inter-parois. Dans les faits, le temps caractéristique de conduction de chaleur pour cette épaisseur d'isolation (0,045m de perlite par exemple) peut être calculé selon l'équation suivante (équation 122) : 2 e _ 0.0452 = 0.66 heure a 8.6 x 10-7 Avec e = l'épaisseur de l'isolation et a = la diffusivité thermique de l'isolation. Le temps requis pour une stabilité thermique du gaz dans le conduit 11 situé à l'intérieur est au moins deux fois la durée t (de l'ordre de 1, 33 heure environ) ce qui est supérieur au temps moyen d'un remplissage (0,4 heure environ). La hauteur de liquide dans le réservoir avant remplissage est estimée à 30% de la hauteur maximale de liquide et la hauteur de gaz peut en être déduite par la suite (équations 123) : hl _before =0.3 H. hg _before = htot_ 0.3 H. Détection des remplissaqes complets Pour détecter si un remplissage est complet parmi une liste de données couvrant plusieurs remplissages, on peut utiliser la procédure suivante en faisant l'hypothèse qu'au moins un remplissage est complet dans la liste de remplissages dont on possède les données. 1) on détermine le différentiel de pression DP maximum juste après un remplissage et pour tous les remplissages 2) un remplissage est considéré total ou complet si la différence relative entre le différentiel de pression juste après remplissage et le différentiel de pression maximum est inférieur à un seuil (par exemple 5%). De préférence seuls les remplissages considérés complets sont utilisés pour déterminer la hauteur totale du réservoir. Il faut souligner que, plus le volume du réservoir est grand, plus la probabilité d'avoir un remplissage incomplet est grande. Ceci s'explique par le fait que les volumes de livraisons par camion sont limités par le stockage du camion. Par conséquent, plus le volume à remplir est grand, plus il sera nécessaire de disposer d'un grand nombre de données de remplissages. Hauteur totale du réservoir et estimation du niveau de liquide maximal. (htot et Hmax )• Pour chaque remplissage considéré complet, la hauteur totale est déterminée en utilisant l'équation 120. Ensuite la valeur moyenne est calculée et donne le troisième paramètre géométrique du réservoir. La hauteur de liquide maximale Hmax est déduite quant à elle à partir de la hauteur totale du réservoir à partir de l'équation 119. Ceci donne un quatrième paramètre géométrique. Pour les réservoirs européens ayant un volume supérieur ou égale à 50m3, le remplissage est rarement complet. Dans ce cas, la détermination de la hauteur totale du réservoir peut être basée sur une estimation globale utilisant l'équation suivante 124 basée sur l'équation précédente 11. htot = Vtot + 2F n R2 3 Ainsi la hauteur maximale de liquide Hmax peut être déduite de l'équation 119. Estimation des pertes thermiques. Les pertes thermiques d'un réservoir cryogénique sont exprimées en général en % d'oxygène perdu par jour. Selon l'invention il apparaît suffisant de réaliser une estimation globale de cette perte (ce paramètre est moins sensible ou important que la hauteur totale). Cette perte est fonction du volume total du réservoir, elle décroit lorsque le 25 volume augmente. Par exemple pour un réservoir de la société Cryolor, une approximation correcte est (équation 125). Oz perte par jour =0.65273 Vt-137149 En connaissant le volume du réservoir en m3, cette dernière équation estime les pertes thermiques journalières ce qui constitue le cinquième paramètre du 30 réservoir. Cette équation peut être utilisée pour d'autres types de réservoirs (d'autres fabricants). Le procédé décrit ci-dessus permet d'estimer avec une bonne précision les paramètres géométriques et de perte thermique à partir de simples mesures de différentiel de pression, de pression et de la masse livrée lors des livraisons. REVENDICATIONS 1. Procédé d'estimation de paramètres caractéristiques d'un réservoir (1) cryogénique et notamment de paramètres géométriques comprenant une étape de calcul d'un premier paramètre géométrique (R) du réservoir en fonction : - d'une masse de liquide livrée (mdelivered) déterminée au cours d'un remplissage, - de la différence entre les différentiels de pression (DPmes) mesurés avant et après remplissage (DPmes afrer ùDPmes before ), chaque différentiel de pression mesurant le différentiel de pression (entre les parties haute et basse du réservoir, - des masses volumiques du gaz et du liquide (pg, pl) dans le réservoir. 2. Procédé selon la revendication 1, caractérisé en ce que le réservoir comporte une portion cylindrique et au moins une extrémité ayant une portion elliptique de hauteur (F) déterminée et en ce que le premier paramètre géométrique et le rayon (R) du cylindre. 3. Procédé selon la revendication 2, caractérisé en ce que le rayon (R ) est calculé selon une équation du type : R = mdelivered g [1ù (MAvo)] 5 10 7L (DPmes _ after ù DPmes _ before dans laquelle g est l'accélération de la pesanteur terrestre et MAVO est un coefficient de masse volumique fonction de la masse volumique du liquide et du gaz dans le réservoir et éventuellement dans des conduits (11) de mesure de pression lorsque la mesure du différentiel de pression est réalisée par au moins un capteur de pression déporté relié aux parties haute et basse du réservoir via des tuyaux (11) de mesure respectifs. 4. Procédé selon la revendication 3, caractérisé en ce que pour l'estimation du coefficient de masse volumique (MAVO) le procédé utilise l'une au moins des hypothèses suivantes :- la masse volumique du liquide dans le réservoir pl est considérée comme étant égale à la moyenne des masses volumiques de liquide avant pi before et après pi after remplissage, - la masse volumique du liquide avant remplissage pl_before est égale à la masse volumique à l'équilibre à la pression du réservoir , - la masse volumique du liquide après remplissage pl after est égale à la moyenne entre d'une part la masse volumique du liquide à l'équilibre à la pression du réservoir avant remplissage pondérée par la fraction du volume occupé par ce liquide et, d'autre part, la masse volumique du liquide dans le camion considérée à l'équilibre à la pression du camion pondérée par la fraction du volume disponible dans le réservoir avant remplissage, - le volume du liquide dans le réservoir avant remplissage est considéré égal à une fraction connue (par exemple 30%) du volume maximal de liquide dans le réservoir, - la masse volumique du gaz dans le réservoir pg est calculée à la pression du réservoir Ptank et pour une température augmentée par rapport à la température d'équilibre à la pression du réservoir (par exemple augmentée de 20 K), 20 - le différentiel de pression DPmes mesuré est corrigé en tenant compte d'une valeur de différence de pression additionnelle (DPplpe ) créée par le gaz présent dans les tuyaux (11) de mesure en cas de mesure déportée. 5. Procédé selon l'une quelconque des revendications 2 à 4, caractérisé 25 en ce que le volume du liquide VI dans un réservoir cylindrique de rayon R et ayant une extrémité elliptique de hauteur F est donné par la relation Si hi F Vi=7LR2~hiù3 2 Sihi 6. Procédé selon l'une quelconque des revendication 2 à 5, caractérisé en ce qu'il comporte une étape de calcul d'un second paramètre géométrique constitué de la hauteur (F) de la portion elliptique du réservoir en fonction de la valeur de rayon (R) calculée, la valeur de la hauteur (F) de la portion elliptique étant donnée par l'équation suivante : R =K F K étant une constante connue représentative d'un type de réservoir ou une constante choisie arbitrairement telle que 1,95 environ. 7. Procédé selon l'une quelconque des revendications 2 à 6, caractérisé en ce qu'il comporte une étape de calcul d'un troisième paramètre géométrique constitué de la hauteur totale du réservoir (htot) à partir du différentiel de pression mesuré juste après un remplissage DPmes after en faisant l'hypothèse que le remplissage est total, ce différentiel de pression mesuré étant exprimé sous la forme suivante : DPmes _ after = A0 H. + Al htot + A2 hg _ before + A3 hl _before Dans laquelle Ao,A1,A2, A3 sont des coefficients dépendant des masses volumiques du gaz et du liquide avant et après remplissage, Hmax est la hauteur maximale de liquide dans le réservoir, hg before étant la hauteur de gaz dans le réservoir avant remplissage, hl before étant la hauteur de liquide dans le réservoir avant remplissage, et en utilisant l'hypothèse suivante : la hauteur de liquide dans le réservoir avant remplissage hl before est estimée à un seuil FS déterminé connu exprimé en pourcentage de la hauteur maximale de liquide Hmax, la hauteur hg before de gaz avant remplissage en étant déduite comme le complément : hl before = FS Hmax hg _before = htot-FS H. 8. Procédé selon la revendication 7, caractérisé en ce que qu'il comporte une étape de calcul d'un quatrième paramètre géométrique constitué par la hauteur de liquide maximale Hmax, celle-ci étant est déduite la hauteur totale du réservoir calculée à partir de l'équation suivante :F Vmax C 2F 1 Hmax =_+ htot ù 3 Vtot 3 9. Procédé selon l'une quelconque des revendication 1 à 8, caractérisé en ce que qu'il comporte une étape de calcul d'un cinquième paramètre géométrique constitué par la perte thermique du réservoir, ladite perte thermique exprimée en pourcentage d'oxygène (%02) d'oxygène perdu par jour étant approximée selon une relation du type : % 02 perte par jour = pi Vtopz dans lequel Vtot est le volume total du réservoir pi est un coefficient de l'ordre de 0,6 et de préférence égal à 0,65273 et p2 un coefficient de l'ordre de 0,3 et de préférence égal à 0,37149 10. Procédé selon l'une quelconque des revendications 1 à 9, caractérisé en ce qu'il comporte une étape de collecte d'une pluralité de valeurs de masse de liquide livrée (mdelivered) déterminées au cours respectivement de plusieurs remplissages du réservoir, une étape de calcul, pour chaque valeur de masse de liquide livrée (mdelivered) du rapport r de remplissage r= mdelivered = constante DPmes_after ù DPmes_before et en ce que ne sont utilisées dans le procédé d'estimation que les valeurs de masse de liquide livrée (mdelivered) pour lesquelles la valeur absolue du rapport r de remplissage ne présente pas une différence supérieure à un seuil fixé vis-à-vis d'une valeur constante de référence déterminée. 11. Procédé selon la revendication 10, caractérisé en ce que vise la valeur constante de référence déterminée pour rapport r de remplissage est constituée de la moyenne de ce rapport de remplissage r calculée pour plusieurs remplissages. 21