EP2205899A2 - Procede de determination en temps reel du niveau de remplissage d'un reservoir cryogenique - Google Patents

Procede de determination en temps reel du niveau de remplissage d'un reservoir cryogenique

Info

Publication number
EP2205899A2
EP2205899A2 EP08840893A EP08840893A EP2205899A2 EP 2205899 A2 EP2205899 A2 EP 2205899A2 EP 08840893 A EP08840893 A EP 08840893A EP 08840893 A EP08840893 A EP 08840893A EP 2205899 A2 EP2205899 A2 EP 2205899A2
Authority
EP
European Patent Office
Prior art keywords
liquid
tank
gas
pressure
temperature
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
EP08840893A
Other languages
German (de)
English (en)
Inventor
Fouad Ammouri
Florence Boutemy
Jonathan Macron
Alain Donzel
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Air Liquide SA
LAir Liquide SA pour lEtude et lExploitation des Procedes Georges Claude
Original Assignee
Air Liquide SA
LAir Liquide SA pour lEtude et lExploitation des Procedes Georges Claude
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Air Liquide SA, LAir Liquide SA pour lEtude et lExploitation des Procedes Georges Claude filed Critical Air Liquide SA
Publication of EP2205899A2 publication Critical patent/EP2205899A2/fr
Withdrawn legal-status Critical Current

Links

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C13/00—Details of vessels or of the filling or discharging of vessels
    • F17C13/02—Special adaptations of indicating, measuring, or monitoring equipment
    • F17C13/021—Special adaptations of indicating, measuring, or monitoring equipment having the height as the parameter
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C13/00—Details of vessels or of the filling or discharging of vessels
    • F17C13/02—Special adaptations of indicating, measuring, or monitoring equipment
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C13/00—Details of vessels or of the filling or discharging of vessels
    • F17C13/02—Special adaptations of indicating, measuring, or monitoring equipment
    • F17C13/023—Special adaptations of indicating, measuring, or monitoring equipment having the mass as the parameter
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C13/00—Details of vessels or of the filling or discharging of vessels
    • F17C13/02—Special adaptations of indicating, measuring, or monitoring equipment
    • F17C13/025—Special adaptations of indicating, measuring, or monitoring equipment having the pressure as the parameter
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C13/00—Details of vessels or of the filling or discharging of vessels
    • F17C13/02—Special adaptations of indicating, measuring, or monitoring equipment
    • F17C13/028—Special adaptations of indicating, measuring, or monitoring equipment having the volume as the parameter
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C5/00—Methods or apparatus for filling containers with liquefied, solidified, or compressed gases under pressures
    • F17C5/02—Methods or apparatus for filling containers with liquefied, solidified, or compressed gases under pressures for filling with liquefied gases
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
    • G01F23/14—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measurement of pressure
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
    • G01F23/80—Arrangements for signal processing
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2221/00—Handled fluid, in particular type of fluid
    • F17C2221/01—Pure fluids
    • F17C2221/011—Oxygen
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2221/00—Handled fluid, in particular type of fluid
    • F17C2221/01—Pure fluids
    • F17C2221/014—Nitrogen
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2221/00—Handled fluid, in particular type of fluid
    • F17C2221/01—Pure fluids
    • F17C2221/016—Noble gases (Ar, Kr, Xe)
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2221/00—Handled fluid, in particular type of fluid
    • F17C2221/01—Pure fluids
    • F17C2221/016—Noble gases (Ar, Kr, Xe)
    • F17C2221/017—Helium
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2223/00—Handled fluid before transfer, i.e. state of fluid when stored in the vessel or before transfer from the vessel
    • F17C2223/01—Handled fluid before transfer, i.e. state of fluid when stored in the vessel or before transfer from the vessel characterised by the phase
    • F17C2223/0146—Two-phase
    • F17C2223/0153—Liquefied gas, e.g. LPG, GPL
    • F17C2223/0161—Liquefied gas, e.g. LPG, GPL cryogenic, e.g. LNG, GNL, PLNG
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2223/00—Handled fluid before transfer, i.e. state of fluid when stored in the vessel or before transfer from the vessel
    • F17C2223/03—Handled fluid before transfer, i.e. state of fluid when stored in the vessel or before transfer from the vessel characterised by the pressure level
    • F17C2223/035—High pressure (>10 bar)
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2250/00—Accessories; Control means; Indicating, measuring or monitoring of parameters
    • F17C2250/03—Control means
    • F17C2250/032—Control means using computers
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2250/00—Accessories; Control means; Indicating, measuring or monitoring of parameters
    • F17C2250/04—Indicating or measuring of parameters as input values
    • F17C2250/0404—Parameters indicated or measured
    • F17C2250/0408—Level of content in the vessel
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2250/00—Accessories; Control means; Indicating, measuring or monitoring of parameters
    • F17C2250/04—Indicating or measuring of parameters as input values
    • F17C2250/0404—Parameters indicated or measured
    • F17C2250/0421—Mass or weight of the content of the vessel
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2250/00—Accessories; Control means; Indicating, measuring or monitoring of parameters
    • F17C2250/04—Indicating or measuring of parameters as input values
    • F17C2250/0404—Parameters indicated or measured
    • F17C2250/0426—Volume
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2250/00—Accessories; Control means; Indicating, measuring or monitoring of parameters
    • F17C2250/04—Indicating or measuring of parameters as input values
    • F17C2250/0404—Parameters indicated or measured
    • F17C2250/043—Pressure
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2250/00—Accessories; Control means; Indicating, measuring or monitoring of parameters
    • F17C2250/04—Indicating or measuring of parameters as input values
    • F17C2250/0404—Parameters indicated or measured
    • F17C2250/043—Pressure
    • F17C2250/0434—Pressure difference
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2260/00—Purposes of gas storage and gas handling
    • F17C2260/02—Improving properties related to fluid or fluid transfer
    • F17C2260/024—Improving metering
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F17—STORING OR DISTRIBUTING GASES OR LIQUIDS
    • F17C—VESSELS FOR CONTAINING OR STORING COMPRESSED, LIQUEFIED OR SOLIDIFIED GASES; FIXED-CAPACITY GAS-HOLDERS; FILLING VESSELS WITH, OR DISCHARGING FROM VESSELS, COMPRESSED, LIQUEFIED, OR SOLIDIFIED GASES
    • F17C2270/00—Applications
    • F17C2270/01—Applications for fluid transport or storage

Definitions

  • the present invention relates to a method for determining in real time the filling level of a cryogenic tank.
  • the invention also relates to a method for determining in real time the quantity of fluid available at each moment.
  • This process requires knowing the parameters (geometry) of the tank. These parameters can be known to the system or estimated according to an estimation method (radius, height %) independent.
  • the invention relates more particularly to a method of determining in real time the filling level of a cryogenic tank intended to contain a two-phase liquid-gas mixture in which is calculated at each time step for the liquid and possibly for the gas at least one of: the level, the volume or the mass contained in the tank, the method measuring at each time step the pressure differential between the high and low parts of the tank and at least one of the pressures of said differential.
  • the invention thus relates to the improvement of level measurement in cryogenic tanks in order to improve the efficiency of the liquid supply logistics chain of these tanks.
  • the reservoirs concerned comprise an internal reservoir for storing the fluid (or inner casing) disposed in an external reservoir (or outer casing) with a vacuum between these reservoirs.
  • the tanks store cryogenic liquids such as oxygen, argon, nitrogen with capacities of 100 liters to 100000 liters for example. Storage pressures can be between 3 bar and 35 bar.
  • deliveries of cryogenic liquid by truck are based on two pieces of information: the time of delivery and the quantity available.
  • the delivery time is based on the liquid crossing of a fixed threshold (delivery threshold usually 30% of the tank capacity). This threshold incorporates measurement uncertainty to prevent the user from drying out. By making reliable the liquid level measurement in the tank it is possible to lower this trigger threshold and thus reduce the delivery frequency and therefore the associated costs.
  • the knowledge of the quantity of fluid available in each tank makes it possible to determine the routes of the delivery trucks. The more the expected quantity of delivery will be close to that actually delivered, the more the tour plan will be respected thus making it possible to fully exploit the potential benefits generated by the many improvements made on the logistics part.
  • the volume measurement is determined using in general an average liquid density for each range of pressure tanks. In practice, this density depends on the evolution of the internal pressures and the temperature. As a result, the higher the pressure in the tank, the greater the volume measurement error.
  • the total mass measurement depends less on the density, hence a better measurement. However, this measure is not sufficient to determine the mass capacity of a tank at the time of delivery (it would be necessary to know the density of the fluid as for the calculation of the volume).
  • a known technique for measuring the level of fluid in a cryogenic tank is to determine a liquid level using the principle of the pressure difference on the height of the tank.
  • volume measurement (% V) is simply obtained from the measurement of a pressure differential (DP) as follows:
  • the density of the liquid k 0 is in fact not constant. It varies over time depending on the internal pressure and the temperature of the liquid when it is not saturated, or when filling with a subcooled liquid. Therefore, the volume result given by the known technique is not always relevant (0-10% error approximately). This phenomenon is especially felt for high pressure (HP) tanks, that is to say greater than 10 bar and especially greater than 15 bar.
  • Liquid level estimates are usually based on regular measurements of the pressure differential (eg every hour), pressure, and reservoir geometry (diameter, height, and maximum level).
  • the estimation errors are therefore mainly due to: a poor conversion of a pressure differential measured in the corresponding liquid level, to the difference between the measured pressure differential and the real pressure differential (effects of the measuring conduits connecting tank sensors), calibration errors (density of the liquid).
  • the measurement of the pressure differential DP is the image of a mass, one thus obtains measurements independent of the evolution of the density. Nevertheless, two masses must be distinguished: the total mass before filling (gas + liquid) and the mass capacity (mass of liquid after filling). The first can be easily obtained by measuring the differential pressure DP.
  • the second requires knowing the characteristics and quantity of the product remaining in the tank in order to estimate the quantity of liquid actually deliverable by a truck.
  • a filling level at 100% that corresponding to the maximum mass that can be placed in the reservoir. This mass corresponds of course to a complete filling of the tank with the sub-cooled liquid. With such a calibration of 100%, we will be sure to never exceed 100%. However, we will never know what is the value corresponding to a full filling, this will depend on the conditions at the time of filling. Therefore we join the constraints of using the value of the volume.
  • the single level of liquid or the mass of liquid is not sufficient to meet both issues: trigger threshold of a delivery and quantity deliverable.
  • An object of the present invention is to overcome all or part of the disadvantages of the prior art noted above.
  • the method according to the invention is essentially characterized in that it calculates: by a thermal model at each time step the average temperatures of the liquid and the gas in the tank from the measured pressure differential and at least one of the pressures of said differential,
  • dP - pl.g.dhl from the calculated density of liquid (with dP the variation of pressure of the liquid, pi the density of the liquid, g the terrestrial acceleration, and dhl the variation in the height of the liquid).
  • embodiments of the invention may include one or more of the following features:
  • the method calculates the density of the gas at each instant from a part of the average temperature of the gas obtained and estimated at the previous instant and, secondly, from the pressure differential and from at least one of the measured pressures of said differential.
  • the thermal model calculates, at each time step, the average temperatures of the liquid and the gas in the tank from a part of the measured pressure differential and at least one of the pressures of said differential and from the liquid and gas temperatures known from the previous moment,
  • the model uses as initial value for the temperature of the liquid and the gas, the initial temperature values obtained during a complete filling of the tank, the level of the liquid at this known instant of complete filling being the known level of the weir of the tank,
  • the model uses the following calculation approximation: the gas after complete filling is at equilibrium liquid-vapor at the pressure of the tank,
  • the model uses the following calculation approximation: the liquid and the gas are constantly isothermal, each in its respective volume, but at respective temperatures which may be different, the model calculates at each time step the mean temperatures of the liquid and the gas in the tank from the mass and energy balances modeled and applied separately to the liquid and the gas contained in the reservoir,
  • the mass and energy balances modeled and applied separately to the liquid and the gas contained in the reservoir are made at a given instant on the basis of the density values and the volume of liquid and gas calculated from an estimate of the the instant and the model iterates the calculation of the average liquid and gas temperatures in the reservoir until the average liquid and a predetermined value towards the temperatures estimated at the preceding instant and in that after convergence the process resumes the steps of calculating the density and liquid level temperatures for the following instant,
  • - mass and energy budgets modeled and applied to gas use the enthalpy equation according to which the variation of the enthalpy of the gas corresponds to the thermal and mass exchanges applied to the gas, that is to say in taking into account at least one of the following exchanges: the exchange of heat between the liquid and the gas, the heat exchange between the outside of the tank and the gas, the heat exchange provided by a possible vaporization heater located usually below the tank, the vaporization of liquid in the tank,
  • the mass and energy balances modeled and applied to the liquid use the enthalpy equation according to which the variation of the enthalpy of the liquid corresponds to the heat and mass exchanges applied to the gases, that is to say in taking into account at least one of the following exchanges: the exchange of heat between the outside and the liquid, the heat exchange between the gas and the liquid, the heat exchange provided by a possible vaporization heater located in general below the tank, the vaporization of liquid in the tank, the consumption of liquid withdrawn by a user of the tank,
  • the method comprises measuring the temperature in the vicinity of the tank, to calculate the heat exchange between the outside and the tank,
  • the tank being of the jacketed type with inter-wall volume under vacuum, characterized in that the additional pressure difference value is calculated by adding or subtracting the gas and liquid levels in the measuring pipes by taking into account the calculated liquid level and neglecting the pressure influence of the pipe portions located in the inter-wall space,
  • the additional pressure difference value is calculated according to a formula of the type:
  • DPpipe pg_pipe.g. (2len_w + total_length)
  • pg_pipe the density of the gas at room temperature (outside the tank)
  • g the Earth's acceleration in m / s 2
  • len_w the thickness of the tank wall and totaljength the total height of the inner tank and in that, when the upper pipe is located inside the tank in the inter-wall space, the additional pressure difference value (DPpipe) is calculated according to a formula of the type:
  • DPpipe DPside_gas + DPsideJiq;
  • DPside_gas being the pressure difference in the portion of the pipe connected to the upper part of the tank and facing the gas in the tank
  • DPsideJiq being the pressure difference in the part of the upper pipe and facing the liquid in the tank
  • the method for dynamically determining the filling level of a cryogenic tank is intended to contain a diphasic liquid-gas mixture, according to any one of the preceding claims, comprising a step of measuring a pressure differential between levels situated respectively at the low and high ends of the tank, characterized in that it comprises a calculation of a volume and / or a mass of liquid in the tank from the pressure differential measured, the known or estimated geometry of the reservoir and at least one density of liquid in the tank, the method further comprising the steps of calculating the following quantities for a moment (t + ⁇ t):
  • a first step of calculating a density of the liquid in the reservoir from the pressure measurements at the level of the low and high ends of the reservoir respectively, and the value estimated at the preceding instant of a temperature Tl of the liquid in the tank, a second stage of calculating the level of liquid in the tank by applying the law in hydrostatic fluid type: dP - pl.g.dhl from the density of the liquid calculated in the previous step a third step of calculating the pressure level at the interface between the liquid phase and the gas phase in the reservoir from the calculated liquid level in the reservoir,
  • a seventh step of calculating the temperature Tg of the gas from the energy balance of the fifth and sixth steps and an eighth step of comparing the temperature T (t + ⁇ t) calculated for the following moment at the seventh step; with and the temperature T (t) estimated for the preceding instant, and when the difference between the temperature T (t + ⁇ t) calculated for the following instant at the seventh stage and the temperature T (t) estimated for the moment preceding is greater than a determined threshold: a step of returning to the second step and of iteration until convergence, when the difference between the temperature T (t + ⁇ t) calculated for the following instant in the seventh step and the temperature T (t) estimated for the previous instant is less than a threshold (convergence): reiteration of the above steps for the moment (t + 2 ⁇ t) with the pressure values measured for this moment.
  • FIG. 1 represents a schematic view illustrating a first example of a cryogenic reservoir for implementing the invention (conduits external to the walls of the tank),
  • FIG. 2 represents a schematic view illustrating a second example of a cryogenic reservoir for implementing the invention (conduits internal to the walls of the tank),
  • FIG. 3 is a simplified and partial representation of the steps implemented by the method according to the invention.
  • FIG. 4 is a simplified and partial representation of the initialization steps implemented by the method according to the invention.
  • the estimation method that will be described below may be implemented by a computer of a tank control system (local or remote).
  • This method includes a pressure measurement (differential DP), an estimate and may include a remote transmission.
  • the pressures are measured via conduits 11, 12 which may be in the inter-wall space of the reservoir (FIG. 2) or outside 11 (FIG. 1).
  • the tank 1 comprises a pressurizing device such as a vaporizer heater 3 able to take liquid to vaporize and reinject it into the tank.
  • This heater 3 conventionally regulates the pressure within the tank 1.
  • tank which stores the fluid
  • the liquid brought by truck when filling may also be considered in the steady state (temperature range of 0 10 K around the balance, for example from 77.2 to 87.9 ° K for nitrogen) .
  • the pressure of the liquid in the truck is chosen, according to the pressure of the tank, between 1 and 2 bar.
  • the liquid is introduced into the tank by pumping.
  • gas and liquid-specific temperatures are considered in the tank but without these temperatures being a function of the location in the tank. That is, in the following the gas temperatures Tg and liquid T1 are average temperatures.
  • the pressure at the interface between the liquid phase and the gaseous phase taken as the average of the pressures at the bottom and at the top (PB and PH), is considered.
  • the energy balance equations are applied separately to the gas and liquid phases of the reservoir.
  • the calculated liquid height hl1 (in m) is calculated according to the formula (equation 1):
  • VlX Ti R 2 M - -
  • the calculated liquid level hl1 is corrected taking into account an additional pressure difference value DPpipe created by the gas present in the measuring pipes 11, 12, both when the pipes 11 are located in the tank ( Figure 2) or out of the tank ( Figure 1). That is, the pressure sensors 4 are deported and "read" fluid-influenced pressures in the conduits 11, 12 connecting them to the upper and lower parts of the reservoir.
  • DPwall is the pressure differential between the two ends of the vertical duct crossing the inter-wall (top or bottom).
  • DPtotJength being the pressure difference due to the gas pressure in the pipe portion 11 connecting the highest point to the measuring member 4 (sensor) remote.
  • DPamb is the pressure difference due to the gas pressure in the pipe portion 11 connecting the lowest point to the remote sensor member 4 (sensor).
  • pg pipe being the density of the gas in the conduit thick_w being the thickness of the wall of the tank
  • T (x) being the temperature at point x
  • P the pressure.
  • the pressure differential DPwall between the two ends of the vertical duct passing through the inter-wall (top or bottom) can be considered substantially identical to the upper and lower parts (only the fact of gas in the duct).
  • the duct extends near the outer casing to "capture" the calories outside the tank and completely vaporize the fluid in the conduit 12 measurement. Between the upper and lower ends of this portion, the pressure is substantially the same (differential of 0.5 bar maximum).
  • pg_pipe is the density of gas at asbestos temperature (the outer conduits 11, 12 are preferably not isolated).
  • len_w being the thickness of the walls of the tank and totaljength being the total height of the tank forming the storage volume.
  • the upper duct 11 Since the upper duct 11 is close to the outer casing, it can not be considered that it contains only gas at ambient temperature. Two temperatures are to be considered:
  • Tside-gas the temperature of the gas in the upper part of the duct 11 (adjacent to the gas phase),
  • Tsidejiq the temperature of the gas in the lower part of the duct 11 (adjacent to the liquid phase of the tank).
  • dpipe being the distance (spacing) between the upper duct 11 and the wall of the inner tank.
  • ⁇ PH PHr + DPside _ gas + DPside _ Uq - DPamb
  • totaljength being the total height of the inner tank, hl and hg the actual heights of liquid and gas in that shell.
  • DPipe being the additional pressure difference value created by the gas present in the measuring pipes 11, 12.
  • equations 11 and 8 indicate that the measured pressure differential DP underestimates the real pressure differential DPreal. And since the liquid level in the internal reservoir and the real pressure differential DPreal are proportional, at low liquid levels this underestimation can be significant.
  • fold being the density of the liquid calibrated by default according to the procedure mentioned above.
  • pg being the density of the gas.
  • the members A, B and C of equation 15 correspond respectively to the effects of the calibration, the effect of the measuring ducts 11, 12 and the effect of the gas in the tank.
  • Equation 15 shows that the error in the liquid level due to the effect of the gas in the tank partially offsets the error due to the effect of the measuring ducts.
  • the members A, B and C of equation 16 here also correspond respectively to the effects of the calibration, the effect of the measuring ducts and the effect of the gas in the tank.
  • thermodynamic model makes it possible to calculate these values.
  • the mass and energy variation is the sum of the incoming masses minus the sum of the outgoing masses.
  • the tank can be divided into two working volumes: one for the gas phase and one for the liquid phase.
  • m vap being the mass flow (in m / s) of incoming gas generated by the
  • Mc m Hg (pec T) + ⁇ vap Hg (Teq) - ⁇ t g ai i e rhG (Tg) + Q - W o Jt prc ⁇ prc vap ⁇ ⁇ - * * "is""ga
  • Hg being the enthalpy of the gas.
  • Hg (T prc ) is the enthalpy of the gas coming from the vaporizer
  • T 8 ⁇ ⁇ prc the temperature of the gas from the vaporizer T prc is substantially equal to the temperature (average) of the gas Tg in the tank (is verified if the return gas duct vaporized the inner tank).
  • the flow m prc of incoming energy due to the vaporizer is relatively low.
  • the gas passing through the vaporizer receives from the outside an energy estimated by H 1 PrC (Hg (TaItIb) -Hl (Tl)) and redistributes this energy by heat exchange against the reservoir towards the liquid and gaseous phases.
  • the heat Q supplied to the gas phase can be evaluated as (equation 25):
  • NI m cons mass flow leaving consumption by the user of the tank.
  • equation 28 the mechanical work is neglected because the liquid is considered incompressible.
  • the power related to the vaporization is not taken into account because this equation is used only in the case where there is precisely no vaporization (liquid temperature below equilibrium temperature 7) ⁇ T eq ).
  • the total thermal power entering the liquid phase is composed of three terms I, II, III (equation 30):
  • hp the exchange coefficient between the ambient and the gas in the tank
  • hlg the coefficient of heat exchange between the gas and the liquid in the tank.
  • the first term (I) represents the natural convection between the outside ambient and the liquid located inside the tank.
  • thermodynamic behavior of the tank between two fillings can be established on the basis of the following conditions:
  • PH PH + pg_pipe.g (H max + of _ len - M)
  • PB PB B - pg pipe.gMl
  • Pr eal DP + pg_ pipe. g (H max-H of _ len)
  • the density of liquid in the duct 9S-P ⁇ e is considered at room temperature in the first case ( Figure 1) and at the equilibrium temperature in the second case ( Figure 2).
  • hp.Stot representing the power exchanged by degree of temperature difference can be calculated according to the mass flow rate of vaporized liquid and the latent heat of liquid oxygen at room temperature (298 ° K) according to the following equation 33:
  • Hlatent being the latent heat of vaporization of oxygen O 2 under pressure
  • fttioss ⁇ 2 being the mass flow rate of oxygen loss characteristic of the reservoir
  • Tamb being the ambient temperature relative to the
  • mlCpl hpStot (Tamb - Tl) - ⁇ ⁇ Tl ⁇ Teq dt
  • equation 38 The solution of equation 37 is given by equation 38 below:
  • the first effect is a condensation of the gas due to the cold liquid injected at the top.
  • the second effect is a gas evacuation through the safety valve.
  • the first effect is privileged while the second effect is privileged when filling from below.
  • Model with condensation the gas evacuated by the valve (model with condensation) is neglected, or it is estimated that all the mass of lost gas is evacuated by the valve (non-condensing model).
  • Hl _ affermi _ after + Hg _ after.mg _ after ml _ bef.Hl _ bef + ml _ deliv.Hl _ deliv + mg _ befilg _ bef
  • the working hypotheses are as follows: the filling is at 100%, that is to say that the liquid level is at most Hmax, the volumes of gas and liquid are known, the temperature of the gas is at equilibrium after filling With these two assumptions it is possible to determine the temperature after filling and the mass delivered during filling. Model without condensation:
  • Hlmix being the enthalpy of the liquid mixture after filling.
  • This first model underestimates the mass actually delivered during filling because it does not consider gas losses (venting through the valve).
  • the second model non-condensing will overestimate the mass actually delivered because the reduction in mass during filling is only due to venting, so it will replace this mass with liquid from the delivery truck.
  • the numerical estimation can be conducted mathematically, we use: the letter u for the measured data used as input
  • the liquid mass ml and the mass of mg gas in the tank are calculated at each time step on the basis of the measured input data.
  • the calculation can comprise three phases:
  • step 101 After initialization (101) the temperature and mass values are calculated by iteration for each time step since the last filling (step 102). If a filling is detected (step 103; "O", step 104), the method calculates the delivered liquid mass mdeliv and the liquid temperatures T1, gas Tg, masses of liquid ml and gas mg after filling according to the model above. As long as a filling is not detected (N, steps 103 and 105), the process continues to calculate by iteration for each time step the values of liquid temperature T1, gas temperature Tg, mass of liquid ml and mass of gas mg.
  • the first assumption is based on the fact that liquid is evacuated during an overflow (the operator must make sure during a filling).
  • the second hypothesis is based on the fact that during a filling a large part of the gas condenses (top filling in particular), the remaining gas can be considered at equilibrium.
  • the real pressure differential value DPreal makes it possible to calculate the density of the liquid at this instant pi according to the formula of step 201.
  • hg the height of gas
  • pg the density of the gas
  • hl the height of liquid
  • g being the terrestrial attraction.
  • a new temperature TM for the liquid can be calculated (step 202, for example applying an equation representing the liquid temperature T1 as a polynomial function of the pi density).
  • the difference between the new calculated temperature value TM and the previous TIo is calculated. If this difference Diff is greater than a threshold S, the procedure starts again at step 200 using the new value of the calculated temperature TM. Otherwise the new temperature TM is adopted as the initial temperature.

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Filling Or Discharging Of Gas Storage Vessels (AREA)

Abstract

Procédé de détermination en temps réel du niveau de remplissage d'un réservoir (1 ) cryogénique destiné à contenir un mélange diphasique liquide-gaz dans lequel est calculé à chaque pas de temps (t, t+ ?t...) pour le liquide et éventuellement pour le gaz au moins l'un parmi : le niveau, le volume ou la masse contenue dans le réservoir (1 ), le procédé mesurant à chaque pas de temps le différentiel de pression (DP=PB-PH) (en Pa) entre les parties haute et basse du réservoir et au moins une des pressions (PH, Pl) dudit différentiel, caractérisé en ce que le procédé : calcule par un modèle thermique à chaque pas de temps (t, t+ ?t...) les températures moyennes du liquide (Tl) et du gaz (Tg) dans le réservoir (1 ) à partir du différentiel de pression (PB-PH) mesuré et d'au moins une des pressions (PH, Pl) dudit différentiel, calcule l'évolution dans le temps (t, t+ ?t...) au moins la masse volumique du liquide (pl) (en kg/m3) à partir de la température moyenne du liquide (Tl) et des pressions (PH) (PB) dans le réservoir, calcule le niveau de liquide (hl) (en m) dans le réservoir (1 ) en appliquant la loi générale de l'hydrostatique au liquide du type : dP=- pl.g.dhl à partir de la masse volumique de liquide (pl) calculée (avec dP la variation de pression du liquide, pl la masse volumique du liquide, g l'accélération terrestre, et dhl la variation de hauteur du liquide).

Description

Procédé de détermination en temps réel du niveau de remplissage d'un réservoir cryogénique
La présente invention concerne un procédé de détermination en temps réel du niveau de remplissage d'un réservoir cryogénique. L'invention concerne également un procédé de détermination en temps réel de la quantité de fluide livrable à chaque instant.
Ce procédé nécessite de connaître les paramètres (géométrie) du réservoir. Ces paramètres peuvent être connus du système ou estimés selon une méthode d'estimation (rayon, hauteur ...) indépendante.
L'invention concerne plus particulièrement un procédé de détermination en temps réel du niveau de remplissage d'un réservoir cryogénique destiné à contenir un mélange diphasique liquide-gaz dans lequel est calculé à chaque pas de temps pour le liquide et éventuellement pour le gaz au moins l'un parmi : le niveau, le volume ou la masse contenue dans le réservoir, le procédé mesurant à chaque pas de temps le différentiel de pression entre les parties haute et basse du réservoir et au moins une des pressions dudit différentiel.
L'invention concerne ainsi 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ée dans un réservoir externe (ou enveloppe externe) avec un vide entre ces réservoirs. 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.
En général, les livraisons de liquide cryogénique par camion sont basées sur deux informations : l'instant de livraison et la quantité livrable. L'instant de livraison est basé sur le franchissement par le liquide d'un seuil fixe (seuil de livraison généralement de 30% de la capacité du réservoir). Ce seuil intègre l'incertitude sur la mesure pour éviter tout assèchement de l'utilisateur du réservoir. En fiabilisant la mesure de niveau de liquide dans le réservoir il est possible de baisser ce seuil de déclenchement et ainsi de diminuer la fréquence de livraison et donc les coûts associés. La connaissance de la quantité de fluide livrable dans chaque réservoir permet de déterminer les tournées des camions de livraison. Plus la quantité prévue d'être livrée sera proche de celle réellement livrée, plus le plan de la tournée sera respecté permettant ainsi d'exploiter pleinement les bénéfices potentiels engendrés par les nombreuses améliorations effectuées sur la partie logistique.
Cette connaissance du niveau de remplissage courant des réservoirs est obtenue par des mesures et calculs d'une grandeur physique représentative par exemple le volume (généralement utilisé France) ou la masse (généralement utilisé en Allemagne).
Les techniques connues d'estimation de ces grandeurs (volume ou masse) sont généralement satisfaisantes mais ne permettent pas d'avoir une précision ou une répétitivité suffisante, notamment pour les réservoirs de stockage à haute pression, ce qui réduit fortement l'efficacité de la chaîne logistique correspondante.
Quelque soit le type de mesure exploitée, il faut connaître l'état du liquide dans le réservoir pour améliorer la précision de la mesure. En effet, la mesure de volume est déterminée en utilisant en général une masse volumique de liquide moyenne pour chaque gamme de pression de réservoirs. En pratique, cette masse volumique dépend de l'évolution des pressions internes et de la température. De ce fait, plus la pression est élevée dans le réservoir, plus l'erreur de mesure de volume est importante.
La mesure de masse totale dépend moins de la masse volumique d'où une meilleure mesure. En revanche cette mesure n'est pas suffisante pour déterminer la capacité en masse d'un réservoir au moment de la livraison (il faudrait pour cela connaître la masse volumique du fluide comme pour le calcul du volume).
Une technique connue de mesure du niveau de fluide dans un réservoir cryogénique consiste à déterminer un niveau de liquide en utilisant le principe de la différence de pression sur la hauteur du réservoir.
Selon cette technique la mesure de volume (%V) est simplement obtenue à partir de la mesure d'un différentiel de pression (DP) comme suit :
Vol%= K.DeltaP Le facteur «K» dépend physiquement de la masse volumique k0 du fluide et de la géométrie du réservoir k (K= k / ko). En pratique, ce facteur est déterminé et utilisé dans la configuration de façon à avoir 100% du volume occupé par le liquide lors du débordement par le trop plein (au moment de la calibration du réservoir).
La masse volumique du liquide k0 n'est en fait pas constante. Elle varie dans le temps en fonction de la pression interne et de la température du liquide lorsqu'il n'est pas saturé, ou lors d'un remplissage avec un liquide sous-refroidi. C'est pourquoi, le résultat en volume donné par la technique connue n'est pas toujours pertinent (0-10% d'erreur environ). Ce phénomène est surtout ressenti pour les réservoirs à Haute Pression (HP), c'est-à-dire supérieur à 10 bar environ et notamment supérieur à 15 bar.
De plus, deux phénomènes supplémentaires conduisent à des dispersions de mesures importantes en haute pression. En effet, dans ce cas, la masse volumique du liquide se rapproche fortement de celle du gaz, or les appareils de mesure et d'estimation connus ne font pas de distinction entre le liquide et le gaz. De plus, le volume de gaz devient supérieur à celui du liquide lors du déclenchement de la livraison. En revanche, en basse pression, la masse volumique du gaz est négligeable, c'est pourquoi les techniques connues donnent de meilleurs résultats aux pressions relativement basses (quelques bar par exemple).
Les estimations de niveaux de liquide sont basées en général sur des mesures régulières du différentiel de pression (par exemple toutes les heures), d'une pression et de la géométrie du réservoir (diamètre, hauteur et niveau maximum).
Les erreurs d'estimation sont donc essentiellement dues : à une mauvaise conversion d'un différentiel de pression mesuré en niveau de liquide correspondant, à l'écart entre le différentiel de pression mesuré et le différentiel de pression réel (effets des conduits de mesure reliant les capteurs au réservoir), les erreurs de calibration (masse volumique du liquide). Dans le cas où c'est la masse du liquide qui doit être mesurée (estimée), la mesure du différentiel de pression DP est l'image d'une masse, on obtient ainsi des mesures indépendantes de l'évolution de la masse volumique. Néanmoins, il faut distinguer deux masses : la masse totale avant remplissage (gaz+liquide) et la capacité massique (masse de liquide après remplissage). La première peut être facilement obtenue par la mesure du différentiel de pression DP.
La seconde nécessite en revanche de connaître les caractéristiques et quantité du produit restant dans le réservoir afin d'estimer la quantité de liquide réellement livrable par un camion.
Comme la masse volumique varie, même en travaillant en masse, il est impossible de relier un niveau physique à un niveau géométrique.
En revanche, on peut définir par exemple un niveau de remplissage à 100% celui correspondant à la masse maximum que l'on puisse mettre dans le réservoir. Cette masse correspond bien sûr à un remplissage complet du réservoir avec du liquide sous-refroidi. Avec une telle calibration du 100%, on sera sûr de ne jamais dépasser le 100%. Cependant, on ne saura jamais quelle est la valeur correspondante à un remplissage complet, ceci dépendra des conditions au moment du remplissage. De ce fait on rejoint les contraintes de l'utilisation de la valeur du volume.
Ainsi, que l'on travaille en masse ou en volume, le seul niveau de liquide ou la masse de liquide n'est pas suffisante pour répondre à la fois aux deux problématiques : seuil de déclenchement d'une livraison et quantité livrable.
Il existe d'autres principes de mesures pour obtenir un niveau de liquide listés ci-après avec leurs inconvénients :
- Mesure par un flotteur : sa mise en œuvre est cependant délicate et la précision est faible avec de gros risques d'hystérésis (frottement, blocage par le givre),
- Mesure via un plongeur immergé : ce système est simple mais il y possibilité d'hystérésis et la mesure dépend de la masse volumique,
- Mesure par la pression hydrostatique différentielle : il est cependant nécessaire de connaître la masse volumique,
- Mesure par la technique du bullage : cependant la canne de mesure peut se boucher, et un débit de gaz excessif conduit à une erreur. De plus, la mesure dépend de la hauteur liquide au dessus de l'orifice inférieur de la canne et de la masse volumique et est de ce fait peu précise.
- Mesure via une jauge résistive : il y a cependant dissipation importante de chaleur et dérive des valeurs de résistance dans le temps,
- Mesure par jauge supraconductrice : cette technique est coûteuse et nécessite un temps de réponse important pour atteindre l'état d'équilibre,
- Mesure par capteur conductimétrique : celui-ci est cependant sensible à la conductivité du liquide,
- Mesure par sonde capacitive : cette technique est cependant sensible au dépôt dans le réservoir, à la constante diélectrique, nécessite un étalonnage lourd et contraignant sur le plan de l'instrumentation,
- Mesure par tube de Taconis : il y a cependant évaporation importante de liquide, et ceci ne fonctionne que pour de l'hélium liquide (de plus il n'y a pas de signal transmissible)
- Mesure par ultrasons : ici le résultat dépend de la surface du liquide (agitée, ou mousse), de la pression, de la température, de la nature du gaz, de poussières, de condensats...
- Mesure par la technique du peson : il y a dans ce cas nécessité de tenir compte des réactions des tuyauteries, de plus cette technique est très coûteuse,
- Mesure par radar : cette méthode est robuste et précise, mais coûteuse,
- Mesure par rayonnement Gamma : cette mesure est cependant sensible à la masse volumique du produit traversé et est contraignante au niveau de la sécurité (source radioactive).
En résumé, soit la méthode est précise mais très coûteuse (et difficile à mettre en œuvre), soit elle est simple, peu coûteuse mais trop sensible à l'environnement du réservoir pour être exploitable de façon suffisamment fiable.
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 calcule : - par un modèle thermique à chaque pas de temps les températures moyennes du liquide et du gaz dans le réservoir à partir du différentiel de pression mesuré et d'au moins une des pressions dudit différentiel,
- l'évolution dans le temps au moins la masse volumique du liquide à partir la température moyenne du liquide et des pressions dans le réservoir,
- le niveau de liquide dans le réservoir en appliquant la loi générale de l'hydrostatique au liquide du type : dP=- pl.g.dhl à partir de la masse volumique de liquide calculée (avec dP la variation de pression du liquide, pi la masse volumique du liquide, g l'accélération terrestre, et dhl la variation de hauteur du liquide).
Sauf indication contraire, les grandeurs physiques sont exprimées en unité SI (m pour les distances, m3 pour les volumes, kg pour les masses, kg/m3 pour les masses volumiques, Pa pour les pressions ou différentiels de pression etc.). Par ailleurs, des modes de réalisation de l'invention peuvent comporter l'une ou plusieurs des caractéristiques suivantes :
- le procédé calcule la masse volumique du gaz à chaque instant à partir d'une part de la température moyenne du gaz obtenue et estimée à l'instant précédent et, d'autre part, du différentiel de pression et d'au moins une des pressions mesurée dudit différentiel.,
- le modèle thermique calcule à chaque pas de temps les températures moyennes du liquide et du gaz dans le réservoir à partir d'une part du différentiel de pression mesuré et d'au moins une des pressions dudit différentiel et d'autre part à partir des températures du liquide et du gaz connues lors de l'instant précédent,
- le modèle utilise comme valeur de départ pour la température du liquide et du gaz, les valeurs de températures initiales obtenues lors d'un remplissage complet du réservoir, le niveau du liquide à cet instant connu de remplissage complet étant le niveau connu du déversoir du réservoir,
- le modèle utilise l'approximation de calcul suivante : le gaz après remplissage complet est à l'équilibre liquide-vapeur à la pression du réservoir,
- le modèle utilise l'approximation de calcul suivante : le liquide et le gaz sont constamment isothermes chacun dans son volume respectif mais à des températures respectives qui peuvent être différentes, - le modèle calcule à chaque pas de temps les températures moyennes du liquide et du gaz dans le réservoir à partir des bilans de masse et d'énergie modélisés et appliqués séparément au liquide et au gaz contenus dans le réservoir,
- les bilans de masse et d'énergie modélisés et appliqués séparément au liquide et au gaz contenus dans le réservoir sont réalisés à un instant sur la base des valeurs de masse volumique et volume de liquide et de gaz calculées à partir d'une estimation des températures de liquide et de gaz à l'instant précédent et en ce que le modèle itère le calcul des températures moyennes du liquide et du gaz dans le réservoir jusqu'à ce que les températures moyennes du liquide et du gaz calculées à un instant convergent d'une valeur prédéterminée vers les températures estimée à l'instant précédent et en ce que après convergence le procédé recommence les étapes de calcul des températures de masse volumique et de niveau de liquide pour l'instant suivant,
- les bilans de masse et d'énergie modélisés et appliqués au gaz utilisent l'équation d'enthalpie selon laquelle la variation de l'enthalpie du gaz correspond aux échanges thermiques et de masse appliqués au gaz, c'est-à-dire en prenant en compte au moins l'un des échanges suivants : l'échange de chaleur entre le liquide et le gaz, l'échange thermique entre l'extérieur du réservoir et le gaz, l'échange thermique fourni par un éventuel réchauffeur de vaporisation situé en général en dessous du réservoir, la vaporisation de liquide dans le réservoir,
- les bilans de masse et d'énergie modélisés et appliqués au liquide utilisent l'équation d'enthalpie selon laquelle la variation de l'enthalpie du liquide correspond aux échanges thermiques et de masse appliqués aux gaz, c'est-à-dire en prenant en compte au moins l'un des échange suivants : l'échange de chaleur entre l'extérieur et le liquide, l'échange thermique entre le gaz et le liquide, l'échange thermique fourni par un éventuel réchauffeur de vaporisation situé en général en dessous du réservoir, la vaporisation de liquide dans le réservoir , la consommation de liquide soutirée par un utilisateur du réservoir,
- le procédé comporte une mesure de température au voisinage extérieur du réservoir, pour calculer l'échange de chaleur entre l'extérieur et le réservoir,
- la mesure du différentiel de pression entre les parties haute et basse du réservoir 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, caractérisé en ce que le procédé corrige le niveau de liquide (hl) calculé et/ou le différentiel de pression mesuré en tenant compte d'une valeur de différence de pression additionnelle créée par le gaz présent dans les tuyaux de mesure, le différentiel de pression 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 » entre les parties hautes et basses du réservoir selon la formule : DP=PB-PH = DPreal- DPpipe,
- le réservoir étant du type à double enveloppe avec volume inter-paroi sous vide d'air, caractérisé en ce que la valeur de différence de pression additionnelle est calculée en ajoutant ou retranchant les niveaux de gaz et de liquide dans les tuyaux de mesure en tenant compte du niveau de liquide calculé et en négligeant l'influence en pression des portions de tuyaux situées dans l'espace inter-paroi,
- lorsque le tuyau supérieur de mesure de pression relié à la partie supérieure du réservoir est situé à l'extérieur du réservoir, la valeur de différence de pression additionnelle est calculée selon une formule du type :
DPpipe = pg_pipe.g.(2len_w+total_length) dans laquelle pg_pipe est la masse volumique du gaz à température ambiante (à l'extérieur au réservoir), g est l'accélération terrestre en m/s2, len_w est l'épaisseur de la paroi du réservoir et totaljength la hauteur totale du réservoir interne et en ce que, lorsque le tuyau supérieur est situé à l'intérieur du réservoir dans l'espace inter-parois, la valeur de différence de pression additionnelle (DPpipe) est calculée selon une formule du type : DPpipe =DPside_gas + DPsideJiq ;
DPside_gas étant 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, DPsideJiq étant la différence de pression dans la partie du tuyau supérieur et se trouvant en face du liquide dans le réservoir ; DPside_gas étant obtenu selon la formule DPside_gas =pSide_gas g hg dans laquelle pSide_gas est la masse volumique du gaz dans le tuyau calculée à une température représentative du tuyau, g étant l'accélération terrestre et hg la hauteur de gaz dans le réservoir ; DPsideJiq étant obtenu selon la formule DPsideJiq =pSide_iiq g hl dans laquelle Psidejiq est la masse volumique du gaz dans le tuyau supérieur se trouvant en face du liquide dans le réservoir calculée à une température représentative du tuyau, g étant l'accélération terrestre et hl la hauteur de liquide dans le réservoir,
- le procédé de détermination en dynamique du niveau de remplissage d'un réservoir cryogénique est destiné à contenir un mélange diphasique liquide-gaz, selon l'une quelconque des revendications précédente comprenant une étape de mesure d'un différentiel de pression entre des niveaux situés respectivement aux extrémités basse et haute du réservoir, caractérisé en ce qu'il comporte un calcul d'un volume et/ou d'une masse de liquide dans le réservoir à partir du différentiel de pression mesuré, de la géométrie connue ou estimée du réservoir et d'au moins une masse volumique de liquide dans le réservoir, le procédé comprenant en outre les étapes de calcul des grandeurs suivantes pour un instant (t+Δt) :
- une première étape de calcul d'une masse volumique du liquide dans le réservoir à partir des mesures de pression au niveau respectivement des extrémités basse et haute du réservoir, et de la valeur estimée à l'instant précédent d'une température Tl du liquide dans le réservoir, une seconde étape de calcul du niveau de liquide dans le réservoir en appliquant la loi en hydrostatique sur le liquide du type : dP=- pl.g.dhl à partir de la masse volumique du liquide calculée à l'étape précédente, une troisième étape de calcul du niveau de pression à l'interface entre la phase liquide et la phase gaz dans le réservoir à partir du niveau de liquide calculé dans le réservoir,
- une quatrième étape de calcul du volume de liquide et de gaz dans le réservoir à partir du niveau de liquide calculé à la seconde étape,
- une cinquième étape de calcul du bilan d'énergie appliqué au gaz du réservoir à partir de l'équation de variation de l'enthalpie du gaz et en prenant en compte au moins l'un des échanges thermiques suivants : l'échange de chaleur entre le liquide et le gaz, l'échange thermique entre l'extérieur du réservoir et le gaz, l'échange thermique fourni par un éventuel réchauffeur de vaporisation situé dans le réservoir, la vaporisation de liquide dans le réservoir, - une sixième étape de calcul du bilan d'énergie appliquée au liquide du réservoir à partir de l'équation de variation de l'enthalpie du liquide et en prenant en compte au moins l'un des échange suivants : l'échange de chaleur entre l'extérieur et le liquide, l'échange thermique entre le gaz et le liquide, l'échange thermique fourni par un éventuel réchauffeur de vaporisation situé dans le réservoir, la vaporisation de liquide dans le réservoir,
- une septième étape de calcul de la température Tg du gaz à partir du bilan d'énergie des cinquième et sixième étapes, et une huitième étape de comparaison de la température T(t+Δt) calculée pour l'instant suivant à la septième étape avec et la température T(t) estimée pour l'instant précédent, et lorsque la différence entre la température T(t+Δt) calculée pour l'instant suivant à la septième étape et la température T(t) estimée pour l'instant précédent est supérieure à un seuil déterminé : une étape de retour à la seconde étape et d'itération jusqu' à convergence, lorsque la différence entre de la température T(t+Δt) calculée pour l'instant suivant à la septième étape et la température T(t) estimée pour l'instant précédent est inférieure à un seuil (convergence) : réitération des étapes ci-dessus pour l'instant (t+2Δt) avec les valeurs de pression mesurées pour cet instant.
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 œuvre 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 œuvre de l'invention (conduits internes aux parois du réservoir),
- la figure 3 représente de façon simplifiée et partielle des étapes mises en œuvre par le procédé selon l'invention,
- la figure 4 représente de façon simplifiée et partielle des étapes d'initialisation mises en œuvre par le procédé selon l'invention. La méthode d'estimation qui va être décrite ci-après peut être mise en œuvre 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 (différentiel DP), une estimation et peut comprendre une transmission à distance. 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 comprend 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 lors des remplissages peut également être considéré à l'état d'équilibre (gamme de température de 100K autour de l'équilibre, par exemple 77,2 à 87,9°K pour de l'azote). La pression du liquide dans le camion est choisie, selon la pression du réservoir, entre 1 et 2 bar. Le liquide est introduit dans le réservoir par pompage.
Lors d'un transfert de liquide depuis un camion dans le réservoir, du liquide relativement plus froid est ajouté au liquide relativement plus chaud dans le réservoir. Lorsque du liquide sort du trop plein, le réservoir est rempli. Après remplissage le réchauffeur de vaporisation 3 modifie la pression si nécessaire. Ainsi, pendant le remplissage une partie du liquide est consommée pour augmenter et stabiliser la pression (de l'ordre de 4kg pour un réservoir de 10000 litres).
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 (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 change et 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 soit 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.
Selon une particularité avantageuse, on considère la pression à l'interface entre la phase liquide et la phase gazeuse, prise comme la moyenne des pressions en bas et en haut (PB et PH).
De préférence on applique séparément les équations de bilan d'énergie aux phases gazeuse et liquide du réservoir.
La méthode ne mesure pas directement le niveau de liquide dans le réservoir mais l'estime sur la base d'une mesure de différentiel DP=PB-PH de pression entre les parties hautes PH et basse PB du réservoir 1 et d'une pression (PB ou PH).
Le niveau de liquide estimé est basé sur le différentiel de pression mesuré DP= PB-PH entre les extrémités basses et hautes du réservoir.
Selon la méthode actuelle, la hauteur de liquide calculée hl1 (en m) est calculée selon la formule (équation 1 ) :
DP{enPa) hll = p/1. g
Avec pli une valeur de masse volumique de liquide de calibration constante en kg/m3 (mais modifiable par un opérateur) ; g étant l'accélération de la pesanteur terrestre en m/s2.
Du fait que le réservoir n'est pas un cylindre géométriquement parfait (extrémités elliptiques, cf. figures 1 et 2), la volume VU 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 hl1 est dans la zone F elliptique
alors VlX =Ti R2 M - -
sinon - hlif
R étant le rayon du réservoir (au niveau de sa portion cylindrique). Le volume maximum Vlmax de liquide dans le réservoir 1 est calculé en remplaçant la valeur de hauteur de liquide hl1 dans l'équation 2 par la hauteur maximale de liquide Hmax (équation 3):
,ι F Vl max = π R \ Hmax
Le pourcentage de liquide total %VI dans le réservoir est calculé en utilisant l'équation suivante (équation 4) :
%n = îoo FZl
Vlmax
La masse de liquide contenue dans le réservoir mil est déduite en utilisant la masse volumique du liquide pli (équationδ) :
»ι/l = p/1 VJl
En utilisant les équations 1 , 2 et 5 ci-dessus la masse de liquide mil peut être exprimée en fonction du différentiel de pression mesuré DP=PB-PH (en Pa).
Suivant que le niveau de liquide hl1 est dans la zone F elliptique (cf. cas (a) ci-dessous) ou en dessous (cf. cas (b) ci-dessous), l'équation devient (équation 6):
DP p li F cas (a) ml 1 = π R 2 g cas (b)
DP mi l = — p /lπ F R 2 - π p /lF - * g
De préférence selon une caractéristique avantageuse possible de l'invention, le niveau de liquide hl1 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) ou 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 DP= PB-PH 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 : PB-PH = DPreal-DPpipe
Cas où les conduits sont à l'extérieur de la paroi du réservoir (figure 1 ) : Dans ce cas, la relation entre les pressions « réelles » PHr, PBr, le différentiel de pression réel DPreal=PBr-PHr et le différentiel mesuré DP=PB-PH est (équation 7) : (cf. figure 5)
PH = PHr - DPwall + DPtot _ length PB = PBr + DPwall - DPamb
DPtot _ length = pg _ pipe . g (len_w + total _leng th - hl 1 ) avec i DPamb = pg _ pipe .g ( hl 1 + len_w) thιck_w
DPwall = g j pg(T(x),P) dx
DPwall étant le différentiel de pression entre les deux extrémités du conduit vertical traversant l'inter-paroi (en haut ou en bas).
DPtotJength étant 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é. pg pipe étant la masse volumique du gaz dans la conduit thick_w étant l'épaisseur de la paroi du réservoir, T(x) étant la température au point x et P la pression. Dans cette équation 7 le différentiel de pression DPwall 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).
Dans l'équation 7, pg_pipe est la masse volumique de gaz à température amiante (les conduits extérieures 11 , 12 sont de préférence non isolés).
Le différentiel de pression réel DPreal=PBr-PHr peut être déduit du différentiel de pression mesuré DP=PB-PH selon la formule :
DP = PB - PH
D Pr eal = DP - 2 DPwall + pg _ pipe. g (2len _ w + total _ length )
len_w étant l'épaisseur des parois du réservoir et totaljength étant la hauteur totale du réservoir formant le volume de stockage.
Cas où les conduits sont dans l'inter-paroi (figure 2) :
Le conduit supérieur 11 étant proche de l'enveloppe extérieure, il n'est pas possible de considérer qu'il contient uniquement du gaz à température ambiante. Deux températures sont à considérer :
Tside-gas = la température du gaz dans la partie supérieure du conduit 11 (adjacente à la phase gazeuse),
Tsidejiq = la température du gaz dans la partie inférieure du conduit 11 (adjacente à a phase liquide du réservoir).
Pour estimer ces températures, on fait l'hypothèse d'un profil de température linéaire entre l'ambiant et le liquide cryogénique du réservoir. Du fait que la température du liquide cryogénique Tl est plus froide que la température du gaz Tg, le profile de température dans le conduit près du liquide est différent du profil de température adjacent au gaz. On obtient l'approximation suivante (équation 9) :
I TsiMae gas = I Tg
I TsiMae l /i'q = a Jpi-peIl Tamb ~ Tl I + , 1 T1/ len w
dpipe étant la distance (espacement) entre le conduit supérieur 11 et la paroi du réservoir intérieur.
De plus, même lorsque le conduit 11 est situé à l'extérieur, il y a une chute de pression due au gaz à température ambiante DPamb.
Il y a donc trois termes correctifs pour la pression mesurée en partie supérieure, chacun de ces trois termes correctifs (I, II, III) correspondant à une masse volumique de gaz différente (équation 10, cf. figure 2) :
\PH = PHr + DPside _ gas + DPside _ Uq - DPamb
[PB = PBr + DPwall - DPamb avec
DPside _ gas = pside _ gas. g hg ; pside _ gas = pg(Tside _ gas, P)
DPside _ Hq = pside _ Hq. g hl ; pside _ Hq = pg(Tside _ Hq, P) hg = total length - hl
totaljength étant la hauteur totale du réservoir intérieur, hl et hg les hauteurs réelles de liquide et de gaz dans cette enveloppe.
Ainsi, la relation entre le différentiel de pression réel DPreal=PBr-PHr et le différentiel de pression mesuré DP=PB-PH est (équation 11 ):
D Pr eal = DP + DPside gas + DPside Uq - DPwall
On peut écrire la formule globale (équation 12) DP = DVreal - DPpipe= (PB - PH)- DPpipe avec: DPpipe=
[- 2 DPwall + pg _ pipe, g il len_w + total _ lengtfi) ( 1 er cas, figure 1 ) [DPside_ gas + DPside_ Hq - DPwall (2nd cas, figure 2))
DPipe étant la valeur de différence de pression additionnelle créée par le gaz présent dans les tuyaux de mesure 11 , 12.
Comme l'épaisseur de l'isolant du réservoir interne (len_w) est de quelques centimètres, on peut négliger le différentiel de pression DPwall entre les deux extrémités du conduit vertical traversant l'inter-paroi. On remarque ainsi que les équations 11 et 8 indiquent que le différentiel de pression mesuré DP sous-estime le différentiel de pression réel DPreal. Et comme le niveau de liquide dans le réservoir interne et le différentiel de pression réel DPreal sont proportionnels, à des niveaux de liquide faibles cette sous-estimation peut être significative.
Cette sous-estimation est également plus importante pour des réservoirs à haute pression et pour le cas de conduits de mesure 11 situés dans l'espace interparois.
Selon l'invention on prend en considération la phase gazeuse dans le réservoir. Le différentiel de pression réel DPreal est en fait (équation 13) :
DVr eal = [pg. kg. g + p/ hl. g]
Avec hg et hl, les hauteurs de gaz et de liquide dans le réservoir.
En substituant cette formue dans la formule 12 on obtient (équation 14) :
DP = g.pl.hl + g.pg.hg -DPpipe.
(les deux derniers membres de cet équation correspondant respectivement à l'effet du gaz dans le réservoir et à l'effet des conduits de mesure).
En utilisant les équations 1 et 14 on peut exprimer la hauteur du liquide hl1 calculée (mesurée selon la technique antérieure) en fonction de la hauteur de liquide réelle corrigée hl (équation15): hn = PL M _ DPpφe + pg p/1 p/.g p/1
pli étant masse volumique du liquide calibrée par défaut selon la procédure citée précédemment. pg étant masse volumique du gaz .
Les membres A, B et C de l'équation 15 correspondant respectivement aux effets de la calibration, à l'effet des conduits de mesure 11 , 12 et à l'effet du gaz dans le réservoir.
Ainsi, le fait d'opter pour une masse volumique fixe de calibration comme dans l'art antérieur génère une erreur proportionnelle à la différence entre d'une part la masse volumique réelle du liquide et d'autre part la masse volumique de liquide constante choisie lors de la calibration.
L'équation 15 montre que l'erreur sur le niveau de liquide due à l'effet du gaz dans le réservoir compense partiellement l'erreur due à l'effet des conduits de mesure.
En exprimant l'erreur relative sur la hauteur de liquide la relation s'écrit (équation 16) :
hll - hl _ pi - pli DPpipe pg.hg hl ~ pli DP pll hl
A B
Les membres A, B et C de l'équation 16 correspondant ici aussi respectivement aux effets de la calibration, à l'effet des conduits de mesure et à l'effet du gaz dans le réservoir.
De l'équation 4 on peut déterminer le volume de liquide avant remplissage VIa (équation 17):
Par conséquent la masse livrée md est égale à (équation 18) : md = (Vlmax- Vla) pll
Cette estimation de niveau ne considère pas que les températures du gaz et du liquide sont toujours à l'équilibre (ceci serait faux en particulier juste après un remplissage). Pour suivre l'évolution de la température au cours du temps, deux modèles peuvent être développés (un modèle dit « détaillé » et un modèle dit « à flux constant »).
Modèle détaillé :
Pour décrire l'évolution du niveau de liquide et de la masse il est nécessaire de connaître les masses volumiques du liquide et du gaz, qui elle-même, dépendent de la pression et de la température. Le modèle thermodynamique décrit ci-après permet de calculer ces valeurs.
En considérant les températures de gaz et de liquide comme homogènes dans leurs volumes respectifs et une isolation constante, on réalise les bilans de masse et d'énergie modélisés et appliqués séparément au liquide et au gaz contenus dans le réservoir.
Pour un volume ouvert, la variation de masse et énergies est la somme des masses entrantes moins la somme des masses sortantes.
Le réservoir peut être divisé en deux volumes de travail : un pour la phase gazeuse et un pour la phase liquide.
En appliquant le bilan de masse à la phase gazeuse (mg = masse gazeuse) on déduit l'équation suivante (équation 19):
dm
, — = " mιvap 4 ~- " mι prc — ' m"g air
Avec mvap étant le débit massique (en m/s) de gaz entrant généré par la
vaporisation de liquide, "*prc étant le flux de gaz entrant provenant du
vaporiseur; et m g_air étant le débit massique de gaz sortant par la valve de sécurité (mise à l'air). Concernant l'énergie de la phase gazeuse, en négligeant les énergies cinétiques et potentielles on obtient le bilan d'énergie suivant (équation 20) :
d(mgeg)
= m mcHg(Tpec ) + ήιvapHg(Teq) - ήtg aiierHg(Tg) + Q - W Jt prc o \ prc vap β ι -* *"i s "" g a
II IU
eg étant l'énergie du gaz, Hg étant l'enthalpie du gaz .
Dans lequel Hg(T prc) est l'enthalpie du gaz venant du vaporiseur, juste
avant de passer en phase gazeuse. Hg(Teq) l'enthalpie du gaz à l'équilibre (Teq
= température à l'équilibre) et Hg(Tg) est l'enthalpie du gaz dans le réservoir (à la température Tg du gaz).
D'après la dernière équation, la variation d'enthalpie est due à :
• m prc**g( 1 prc ) = l'énergie entrant due au vaporiseur.
• ^1 VUpHg(TCq) = l'énergie entrant due à la vaporisation de liquide.
m g _air**g( * g ) = l'énergie sortant lorsque la valve de mise à l'air est passante.
• Q = la puissance thermique à reçue par la phase gazeuse, c'est la somme de la contribution de l'interface gaz/liquide et de la surface du réservoir interne en contact avec la phase gazeuse.
• W = la puissance mécanique à travers les surfaces latérales de la phase gazeuse.
Pg
En exprimant l'énergie du gaz eg = Hg j (Pg étant la pression du gaz
et pg étant sa masse volumique) la partie gauche de la formule de bilan d'énergie devient (équation 21 ) : eg dmg + mg^- - Vg^- + MgPg dÇ>g dt dt dt 9g' dt
II est possible de trouver une relation exprimant la variation de masse dog volumique du gaz — — selon la formule (équation 22) :
dmg. = Vg *pg_+ 9g *vg dt dt dt
(Vg= volume de gaz).
Si la compressibilité du gaz est prise en compte, la seule puissance mécanique W à travers la surface est produite par l'augmentation ou la diminution du niveau de liquide (équation 23) :
dVg
W = Pg dt
Finalement, en substituant les équations 20, 21 , 22 dans l'équation 23, le bilan d'enthalpie de la phase gazeuse s'écrit (équation 24) :
mg^ = fkmc(Hg(T ) - Hg(Tg)) + m (Hg(Teq) - Hg(Tg)) + Q + Vg dPg dt dt
Pour simplifier le modèle on fait les approximations suivantes : - Cas où les conduits sont intra-parois :
T8 ≈ τ prc la température du gaz provenant du vaporiseur Tprc est sensiblement égale à la température (moyenne) du gaz Tg dans le réservoir (se vérifie si le conduit de retour de gaz vaporisé épouse le réservoir intérieur). - Cas où les conduits sont externes aux parois :
Tamb ≈ T prc Ainsi, le flux m prc d'énergie entrant due au vaporiseur est relativement faible. Le gaz passant dans le vaporiseur reçoit de l'extérieur une énergie estimée par H1 PrC (Hg(TaItIb) - Hl(Tl)) et redistribue cette énergie par échange thermique contre le réservoir vers les phases liquide et gazeuse.
La chaleur Q fournie à la phase gazeuse peut être évaluée comme (équation 25) :
hg hpSg(Tamb-Tg)+hlëSlëÇri-Tg)+-—rz m (Hg(Tamd) -Hl (Tl)) total _ length
III
Q = pour les conduits intra- parois hp.Sg(Taml^Tg)+h\gS\g(Tl-Tg) + m (Hg(Tamfy-Hl (Tl))
1 U III pour les conduits externes aux parois avec hp:le coefficient d'échange entre l'ambiant et le gaz dans le réservoir hlg:le coefficient d'échange thermique entre le gaz et le liquide dans le réservoir SIg=TlR2 = surface
Les termes I, II, et III correspondent respectivement à :
I = la chaleur venant de l'ambiant extérieur vers le gaz dans le réservoir à travers les parois
II = l'échange de chaleur entre le liquide et le gaz à l'interface et à travers les parois du réservoir
111= la puissance libérée par le gaz passant à travers le vaporiseur vers l'intérieur du réservoir à la phase gazeuse (proportionnel à la hauteur de gaz présent dans le réservoir pour les conduits intra-parois)
Ainsi, en substituant l'équation 25 dans l'équation 24, le bilan d'enthalpie pour la phase gazeuse Hg s'écrit (équation26) : mg— dHrg == hhppSSggTTaammbb--TTgg)) ++ h ItlIggSSlIgg((TTll--TTgg)) ++ —^ hg; -m prc( Hg(T amb)- Hl(Tl)) dt total length dPg
+ mvvmap(Hg(Teq) - Hg(Tg))+ Vg- dt
Les mêmes relations que ci-dessus appliquées à la phase liquide donnent les équations 27 et 28 suivantes (« g » remplacé par « I »): dml j . ~ ~ mvap ~ mprc ~ mcons
^Ml = _ m Hl(Tl) - mconsHl(Tl) + Q dt ^l ^ n m
Dans l'équation 27 de bilan de masse ci-dessus, les membres I, II et III définissent :
I = mvap le débit massique sortant de vaporisation
II = mprc le débit massique sortant vaporisé dans le vaporiseur
NI = mcons \e débit massique sortant de consommation par l'utilisateur du réservoir.
En ce qui concerne le bilan d'énergie appliqué à la phase liquide, les termes I, II, IN de l'équation 28 ci-dessus ont les significations suivantes :
I = mprcHl(Tl) = puissance sortant lorsque le vaporiseur fonctionne
N= mconsHl(Tl) = puissance sortant lorsque du liquide est consommé par le client
NI = Q = puissance thermique interne
Dans l'équation 28 le travail mécanique est négligé car le liquide est considéré incompressible. La puissance liée à la vaporisation n'est pas prise en compte car cette équation est utilisée uniquement dans le cas où il n'y a justement pas de vaporisation (température du liquide inférieure à la température d'équilibre 7} < Teq ).
En revanche, lorsque la température du liquide est supérieure à la température d'équilibre ^i - Teq la puissance entrant dans le système ne va pas augmenter la température du liquide mais le vaporiser. Il y a alors la relation suivante (équation 29) :
mvap(Hg(Teq)- Hl(Tl)) = Q
Tout comme pour la phase liquide, la puissance thermique totale entrant dans la phase liquide est composée de trois termes I, II, III (équation 30) :
hp.Sg (Tamb - Tl) + h \g S Ig(Tl - Tg ) + -^ m prc (Hg (Tamb ) - Hl (Tl)) v v ' v v ' total length i π ^ - ° v , m Q = \ pour les conduits intra - parois hp.Sg (Tamb - Tl) + h \g S Ig(Tl - Tg )
T ÎΓ pour les conduits externes aux parois
avec hp : le coefficient d'échange entre l'ambiant et le gaz dans le réservoir hlg : le coefficient d'échange thermique entre le gaz et le liquide dans le réservoir.
S Ig = TlR2 = surface
Le premier terme (I) représente la convection naturelle entre l'ambiant extérieur et le liquide situé à l'intérieur du réservoir.
Le second (II) représente la convection naturelle entre le liquide et le gaz au niveau de l'interface et le troisième (III) représente l'autre partie d'énergie libérée par le gaz passant à travers le vaporiseur 3 et qui n'est pas pris en compte dans l'équation appliquée au gaz. En utilisant les équations 27 à 30 ci-dessus, le bilan d'enthalpie pour la phase liquide peut s'exprimer comme suit (équation 31 ):
a) si Tl < Teq IIgg((TTgg -- TTll)) ++ —— -^——mprc (Hg(Tamb) - Hl(Tl)) b) si Tl ≥ Teq Hl = Hl(Teq) hl hpSg(Tam - Tl) + hlg S Ig(Tg - Tl)+ — - m prc (Hl(Tamb) - Hl(Tl)) total _ length
™vap =
(Hg( Teq) - Hl(Tl))
Modèle à flux constant :
Un modèle simplifié du comportement thermodynamique du réservoir entre deux remplissages peut être établi sur la base des conditions suivantes :
1 ) le taux de consommation de liquide prélevé est constant
2) le gaz est toujours à l'équilibre liquide-gaz, aussi bien dans le réservoir que dans les conduits,
3) la chaleur spécifique du liquide est constante C pi ,
4) l'effet de vaporiseur est négligeable
5) la totalité de la chaleur transitant à travers les parois est absorbée par le liquide.
Dans ces conditions, la correction du différentiel de pression mesuré DP=PB- PH peut être exprimée dans les deux cas présentés ci-dessus sous la forme suivante (équation 32) :
PH = PH + pg _ pipe.g(H max+ of _ len - M) PB = PBB - pg pipe.gMl D Pr eal = DP + pg _ pipe. g (H max-H of _ len)
Dans ce modèle, la masse volumique de liquide dans le conduit 9S - PΨe est considérée à la température ambiante dans le premier cas (figure 1 ) et à la température d'équilibre dans le second cas (figure 2). hp.Stot qui représente la puissance échangée par degré de différence de température peut être calculé en fonction du débit massique de liquide vaporisé et la chaleur latente de l'oxygène liquide à température ambiante (298°K) selon l'équation 33 suivante :
•
AHlatent02atm . m I0SsO2 hpStot = — ≈ — 2-
Tamb - Teq atm
Hlatent étant la chaleur latente de vaporisation de l'oxygène O2 à pression
atmosphérique, fttiossθ2 étant le débit massique de perte d'oxygène caractéristique du réservoir, Tamb étant la température ambiante relative à la
mesure de tn iossθτ .
L'évolution de la température du liquide est donnée par l'équation 34 suivante :
mlCpl = hpStot(Tamb - Tl)—^→Tl < Teq dt
Tl = Teq—?!→Tl ≥ Teq
Avec CpI = la chaleur spécifique du liquide est constante.
De la première condition la variation de masse du liquide peut être déduite par la formule 35 suivante :
ml = minit - ήlconst
Avec ml = masse du liquide ; minit = masse initiale et mcons = débit massique consommé par le client.
En substituant cela dans les équations 33 et 34 et en considérant la température ambiante Tamb constante, l'équation différentielle suivante est obtenue (équation 36) : dTl _ d(Tamb - Tl) _ hpStot(Tamb - Tl) dt dt (minit - mconst)Cpl
Avec Tl = température du liquide
Ceci est une équation différentielle séparable dont la solution peut être obtenue sans difficulté par intégration entre la température T à l'instant initial Tinit et le température à l'instant T(t) (équation 37):
T( ft) d d((TTaammbb -- TTll)) _ f f hpStotJt
Ti Jnit ( (TTaUmmhb -~ ΗTl)) ~ J 0 (tHinit ~ "1COHs* )Q>1
La solution de l'équation 37 est donnée par l'équation 38 ci-dessous :
T rli(st*) \ = T ramb u - I [Tramb u si >T rw,l, < T rpeq
Tl = Teq^→Tl ≥ Teq
Estimation de la température après remplissage :
Le remplissage par du liquide plus froid entraîne une baisse de température dans le réservoir. Dans les deux modèles (modèle à flux constant et modèle détaillé) décrits ci-dessus il y a réalisation d'un bilan de l'enthalpie entre le fluide restant dans le réservoir et le liquide livré. Ceci permet d'estimer la température du liquide après remplissage et la masse de liquide livrable.
Pendant un remplissage il y a une diminution de la masse de gaz ce qui peut être lié à deux effets distincts.
Le premier effet est une condensation du gaz due au liquide froid injecté en partie haute. Le second effet est une évacuation de gaz par la valve de sécurité.
En remplissant par le haut le premier effet est privilégié tandis que le second effet est privilégié lors d'un remplissage par le bas.
Dans les deux modèles d'estimation considérés ci-après soit on néglige le gaz évacué par la valve (modèle avec condensation), soit on estime que toute la masse de gaz perdue est évacué par la valve (modèle sans condensation). Modèle avec condensation :
En appliquant le bilan d'enthalpie entre le fluide restant dans le camion et le liquide du camion l'équation résultante s'écrit (équation 39) :
Hl _ affermi _ after + Hg _ after.mg _ after = ml _ bef.Hl _ bef + ml _ deliv.Hl _ deliv + mg _ befilg _ bef
Hl_after étant l'enthalpie de liquide après livraison, ml_after étant la masse de liquide après livraison, Hl_bef étant l'enthalpie de liquide avant livraison, ml_bef étant la masse de liquide avant livraison, Hl deliv étant l'enthalpie du liquide livrée, ml_deliv étant la masse de liquide livrée, les mêmes désignations s'appliquent au gaz en remplaçant « I » par « g ».
Dans cette équation la masse livrée dans le réservoir WlI _ dellV est (équation 40) : ml _ deliv = ml _ after + mg _ after - ml _ bef - mg _ bef
= Vl _ max . p/ _ after + Vg _ after. pg _ after - ml _ bef - mg _ bef
Pour obtenir une équation fonction uniquement de la température du liquide Tl après remplissage, les hypothèses de travail sont les suivantes : le remplissage est à 100%, c'est-à-dire que le niveau de liquide est au maximum Hmax , les volumes de gaz et de liquide sont connus, la température du gaz est à l'équilibre après remplissage Grâce à ces deux hypothèses il est possible de déterminer la température après remplissage et la masse livrée lors du remplissage. Modèle sans condensation :
Dans ce modèle le gaz n'est pas considéré, les équations résultantes sont (équation 41 ) :
Hl _ mix.ml _ after = ml _ bef.Hl _ bef + ml _ deliv. Hl _ deliv avec ml _ deliv = ml _ after - ml _ bef = Vl _ max . p/ _ after - ml _ bef
Hlmix étant l'enthalpie du mélange liquide après remplissage. On conserve les mêmes hypothèses pour calculer la température après remplissage et la masse livrable lors du remplissage. Ce premier modèle sous-estime la masse réellement livrée lors du remplissage car il ne considère pas les pertes de gaz (mise à l'air par la valve). Le second modèle (sans condensation) va surestimer la masse réellement livrée car il que la diminution de masse pendant le remplissage est uniquement due à la mise à l'air, il va donc remplacer cette masse par du liquide provenant du camion de livraison.
L'estimation numérique peut être conduite mathématiquement, on utilise : la lettre u pour les données mesurées utilisées comme en entrée
(différentiel de pression mesuré DP et une pression mesurée, par exemple PB), la lettre x pour les variables (température du liquide et du gaz (Tl
Tg), la lettre p pour les paramètres constants
Selon cette méthode on calcule à chaque pas de temps la masse de liquide ml et la masse de gaz mg dans le réservoir sur la base des données d'entrée mesurées. Après la mesure des grandeurs physiques, le calcul peut comporter trois phases :
1 ) initialisation
2) simulation du comportement du réservoir entre deux remplissages,
3) réinitialisation lorsqu'un remplissage est détecté.
La figure 3 illustre schématiquement le processus : après mesure des grandeurs physiques (u= différentiel de pression étape 100), un processus d'initialisation est effectué (étape 101 , figure 3).
Après initialisation (101) on calcule par itération pour chaque pas de temps les valeurs de températures et de masse depuis le dernier remplissage (étape 102). Si un remplissage est détecté (étape 103; « O », étape 104), le procédé calcule la masse de liquide livrée mdeliv et les températures de liquide Tl, de gaz Tg, les masses de liquide ml et de gaz mg après remplissage selon le modèle ci- dessus. Tant qu'un remplissage n'est pas détecté (N, étapes 103 et 105), le procédé continue de calculer par itération pour chaque pas de temps les valeurs de températures de liquide Tl, de température de gaz Tg, masse de liquide ml et de masse de gaz mg. Pour initialiser les premières valeurs de température du liquide Tl et du gaz Tg, le procédé fait les hypothèses suivantes : le premier point de calcul correspond au moment d'un remplissage complet (100%) lorsque le niveau de liquide hl =Hmax est maximum et le niveau de gaz minimum le gaz est à saturation à l'état d'équilibre.
La première hypothèse est basée sur le fait que du liquide est évacué lors d'un trop plein (l'opérateur doit s'en assurer lors d'un remplissage). La seconde hypothèse est basée sur le fait que lors d'un remplissage une grande partie du gaz se condense (remplissage par le haut notamment), le gaz restant peut être considéré à l'équilibre.
L'itération du calcul jusqu'à convergence permet de conduire à une valeur initiale fiable pour la température du liquide TIO (figure 4).
Le différentiel de pression réel DPreal peut alors être calculé par l'équation 12 en fonction de cette température initiale de liquide TIO, en fonction de la température du gaz Tg, en fonction du différentiel de pression mesuré DP=u et en fonction des constantes p (étape 200 figure 4).
La valeur de différentiel de pression réel DPreal permet de calculer la masse volumique du liquide à cet instant pi selon la formule de l'étape 201. Avec hg la hauteur de gaz, pg la masse volumique du gaz, hl, la hauteur de liquide, g étant l'attraction terrestre.
En fonction de la pression et la masse volumique de liquide obtenues, on peut calculer une nouvelle température TM pour le liquide (étape 202, application par exemple d'une équation représentant la température de liquide Tl comme une fonction polynomiale de la masse volumique pi).
La différence entre la nouvelle valeur de température calculée TM et la précédente TIo est calculée. Si cette différence Diff est supérieure à un seuil S, la procédure recommence à l'étape 200 en utilisant la nouvelle valeur de la température calculée TM . Sinon la nouvelle température TM est adoptée comme température initiale.
Ensuite, après initialisation (température de liquide initiale adoptée), l'équation différentielle pour les températures de gaz et de liquide est résolue par itération (par exemple par la méthode d'Euler).

Claims

Selon ces équations, la procédure d'estimation peut se dérouler comme suit :Initialisation de l'étape d'intégration avec les températures de liquide et de gaz obtenues à l'itération précédente et mesure du différentiel de pression DP,Correction du différentiel selon l'équation 12Calcul des hauteurs, volume, masse et masse volumique du liquide et du gazEstimation des nouvelles valeurs de température de liquide et de gaz par discrétisation des équations 26 et 31Calcul des différences entre les nouvelles températures obtenues et les températures précédentes,Si les différences sont supérieures à un seuil de convergence, on recommence l'itération avec les dernières valeurs de température obtenues.Ce processus est utilisé pour estimer l'évolution temporelle des températures tant qu'un nouveau remplissage n'est pas détecté. Un remplissage est détecté lors d'une augmentation de pression différentielle mesurée supérieure à un seuil fixé.La détermination du niveau de remplissage du réservoir permet de déterminer ainsi la quantité de fluide livrable dans le réservoir à chaque instant. REVENDICATIONS
1. Procédé de détermination en temps réel du niveau de remplissage d'un réservoir (1 ) cryogénique destiné à contenir un mélange diphasique liquide-gaz dans lequel est calculé à chaque pas de temps (t, t+ Δt...) pour le liquide et éventuellement pour le gaz au moins l'un parmi : le niveau, le volume ou la masse contenue dans le réservoir (1 ), le procédé mesurant à chaque pas de temps le différentiel de pression (DP=PB-PH) (en Pa) entre les parties haute et basse du réservoir et au moins une des pressions (PH, Pl) dudit différentiel, caractérisé en ce que le procédé : calcule par un modèle thermique à chaque pas de temps (t, t+ Δt...) les températures moyennes du liquide (Tl) et du gaz (Tg) dans le réservoir (1 ) à partir du différentiel de pression (PB-PH) mesuré et d'au moins une des pressions (PH, Pl) dudit différentiel ,
- calcule l'évolution dans le temps (t, t+ Δt...) au moins la masse volumique du liquide (pi) (en kg/m3) à partir la température moyenne du liquide (Tl) et des pressions (PH) (PB) dans le réservoir,
- calcule le niveau de liquide (hl) (en m) dans le réservoir (1 ) en appliquant la loi générale de l'hydrostatique au liquide du type : dP=- pl.g.dhl à partir de la masse volumique de liquide (pi) calculée (avec dP la variation de pression du liquide, pi la masse volumique du liquide, g l'accélération terrestre, et dhl la variation de hauteur du liquide).
2. Procédé selon la revendication 1 , caractérisé que ce que le procédé calcule la masse volumique du liquide (pi) et à chaque instant (t+ Δt...) à partir d'une part de la température moyenne du liquide (Tl) obtenue et estimée à l'instant précédent (t) et, d'autre part, du différentiel de pression (PB-PH) et d'au moins une des pressions (PH, Pl) mesurée dudit différentiel (DP).
3. Procédé selon la revendication 1 ou 2, caractérisé en ce que le procédé calcule la masse volumique du gaz (pg) à chaque instant (t+ Δt...) à partir d'une part de la température moyenne du gaz (Tg) obtenue et estimée à l'instant précédent (t) et, d'autre part, du différentiel de pression (DP=PB- PH) et d'au moins une des pressions (PH, Pl) mesurée dudit différentiel.
4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que le modèle thermique calcule à chaque pas de temps (t+ Δt) les températures moyennes du liquide et du gaz (Tl, Tg) dans le réservoir (1 ) à partir d'une part du différentiel de pression (DP=PB-PH) mesuré et d'au moins une des pressions (PH, Pl) dudit différentiel et d'autre part à partir des températures du liquide et du gaz connues (Tl, Tg) lors de l'instant précédent (t).
5. Procédé selon la revendication 4, caractérisé en ce que le modèle utilise comme valeur de départ pour la température du liquide (Tl) et du gaz (Tg), les valeurs de températures initiales obtenues lors d'un remplissage complet du réservoir (1 ), le niveau du liquide (hl) à cet instant connu de remplissage complet étant le niveau connu du déversoir du réservoir (1 ).
6. Procédé selon la revendication 5, caractérisé en ce que le modèle utilise l'approximation de calcul suivante : le gaz après remplissage complet est à l'équilibre liquide-vapeur à la pression du réservoir.
7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce que le modèle utilise l'approximation de calcul suivante : le liquide et le gaz sont constamment isothermes chacun dans son volume respectif mais à des températures respectives qui peuvent être différentes
8. Procédé selon l'une quelconque des revendication 1 à 7, caractérisé en ce que le modèle calcule à chaque pas de temps les températures moyennes du liquide et du gaz (Tl, Tg) dans le réservoir (1 ) à partir des bilans de masse (m) et d'énergie modélisés et appliqués séparément au liquide et au gaz contenus dans le réservoir (1 ).
9. Procédé selon la revendication 8, caractérisé en ce que les bilans de masse et d'énergie modélisés et appliqués séparément au liquide et au gaz contenus dans le réservoir (1 ) sont réalisés à un instant (t+ Δt) sur la base des valeurs de masse volumique et volume de liquide et de gaz calculées à partir d'une estimation des températures de liquide et de gaz à l'instant précédent (t) et en ce que le modèle itère le calcul des températures moyennes du liquide et du gaz (Tl, Tg) dans le réservoir (1 ) jusqu'à ce que les températures moyennes du liquide et du gaz (Tl, Tg) calculées à un instant (t+ Δt) convergent d'une valeur prédéterminée vers les températures estimée à l'instant précédent (t) et en ce que après convergence le procédé recommence les étapes de calcul des températures (Tl, Tg) de masse volumique (pi) et de niveau de liquide (hl) pour l'instant suivant (t+ 2Δt).
10. Procédé selon la revendication 8 ou 9, caractérisé en ce que les bilans de masse et d'énergie modélisés et appliqués au gaz utilisent l'équation d'enthalpie selon laquelle la variation de l'enthalpie (Hg) du gaz correspond aux échanges thermiques et de masse appliqués au gaz, c'est-à- dire en prenant en compte au moins l'un des échanges suivants : l'échange de chaleur entre le liquide et le gaz, l'échange thermique entre l'extérieur du réservoir (1 ) et le gaz, l'échange thermique fourni par un éventuel réchauffeur (3) de vaporisation situé en général en dessous du réservoir, la vaporisation de liquide dans le réservoir (1 ).
11. Procédé selon l'une quelconque des revendications 8 à 10, caractérisé en ce que les bilans de masse et d'énergie modélisés et appliqués au liquide utilisent l'équation d'enthalpie selon laquelle la variation de l'enthalpie du liquide (Hl) correspond aux échanges thermiques et de masse appliqués aux gaz, c'est-à-dire en prenant en compte au moins l'un des échange suivants : l'échange de chaleur entre l'extérieur et le liquide, l'échange thermique entre le gaz et le liquide, l'échange thermique fourni par un éventuel réchauffeur de vaporisation situé en général en dessous du réservoir, la vaporisation de liquide dans le réservoir (1 ), la consommation de liquide soutirée par un utilisateur du réservoir (1 ).
12. Procédé selon la revendication 10 ou 11 , caractérisé en ce qu'il comporte une mesure de température au voisinage extérieur du réservoir (1), pour calculer l'échange de chaleur entre l'extérieur et le réservoir (1 ).
13. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la mesure du différentiel de pression (DP=PB-PH) entre les parties haute et basse du réservoir 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, caractérisé en ce que le procédé corrige le niveau de liquide (hl) calculé et/ou le différentiel de pression (PB-PH) mesuré 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, le différentiel de pression (DP=PB-PH) 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 » entre les parties hautes et basses du réservoir (DPreal) selon la formule : DP=PB- PH = DPreal-DPpipe.
14. Procédé selon la revendication 13, le réservoir étant du type à double enveloppe avec volume inter-paroi sous vide d'air, caractérisé en ce que la valeur de différence de pression additionnelle (DPpipe) est calculée en ajoutant ou retranchant les niveaux de gaz et de liquide dans les tuyaux de mesure en tenant compte du niveau de liquide (hl) calculé et en négligeant l'influence en pression des portions de tuyaux situées dans l'espace interparoi.
15. Procédé selon la revendication 14, caractérisé en ce que lorsque le tuyau supérieur de mesure de pression relié à la partie supérieure du réservoir est situé à l'extérieur du réservoir, la valeur de différence de pression additionnelle (DPpipe) est calculée selon une formule du type : DPpipe = pg_pipe.g.(2len_w+total_length) dans laquelle pg_pipe est la masse volumique du gaz à température ambiante (à l'extérieur au réservoir), g est l'accélération terrestre, len_w est l'épaisseur de la paroi du réservoir et totaljength la hauteur totale du réservoir interne et en ce que, lorsque le tuyau supérieur est situé à l'intérieur du réservoir dans l'espace inter-parois, la valeur de différence de pression additionnelle (DPpipe) est calculée selon une formule du type : DPpipe =DPside_gas + DPsideJiq ;
DPside_gas étant 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, DPsideJiq étant la différence de pression dans la partie du tuyau supérieur et se trouvant en face du liquide dans le réservoir ; DPside_gas étant obtenu selon la formule DPside_gas =pSide_gas g hg dans laquelle pSide_gas est la masse volumique du gaz dans le tuyau calculée à une température représentative du tuyau, g étant l'accélération terrestre et hg la hauteur de gaz dans le réservoir ; DPsideJiq étant obtenu selon la formule DPsideJiq =pSide_iiq g hl dans laquelle Psidejiq est la masse volumique du gaz dans le tuyau supérieur se trouvant en face du liquide dans le réservoir calculée à une température représentative du tuyau, g étant l'accélération terrestre et hl la hauteur de liquide dans le réservoir.
16. Procédé de détermination en dynamique du niveau de remplissage d'un réservoir (1) cryogénique destiné à contenir un mélange diphasique liquide-gaz, selon l'une quelconque des revendications précédente comprenant une étape de mesure d'un différentiel de pression (DP=PB-PH) entre des niveaux situés respectivement aux extrémités basse et haute du réservoir (1 ), caractérisé en ce qu'il comporte un calcul d'un volume et/ou d'une masse de liquide dans le réservoir (1) à partir du différentiel de pression (DP) mesuré, de la géométrie connue ou estimée du réservoir (1 ) et d'au moins une masse volumique de liquide dans le réservoir (1 ), le procédé comprenant en outre les étapes de calcul des grandeurs suivantes pour un instant (t+Δt) :
- une première étape de calcul d'une masse volumique (pi) du liquide dans le réservoir (1 ) à partir des mesures de pression (PB, PH) au niveau respectivement des extrémités basse et haute du réservoir (1 ), et de la valeur estimée à l'instant précédent (t) d'une température Tl du liquide dans le réservoir (1 ), une seconde étape de calcul du niveau de liquide (hl) dans le réservoir (1) en appliquant la loi en hydrostatique sur le liquide du type : dP=- pl.g.dhl à partir de la masse volumique (pi) du liquide calculée à l'étape précédente, une troisième étape de calcul du niveau de pression (Pl) à l'interface entre la phase liquide et la phase gaz dans le réservoir (1 ) à partir du niveau de liquide (hl) calculé dans le réservoir (1 ),
- une quatrième étape de calcul du volume de liquide (hl) et de gaz dans le réservoir (1 ) à partir du niveau de liquide calculé à la seconde étape, une cinquième étape de calcul du bilan d'énergie appliqué au gaz du réservoir (1 ) à partir de l'équation de variation de l'enthalpie du gaz et en prenant en compte au moins l'un des échanges thermiques suivants : l'échange de chaleur entre le liquide et le gaz, l'échange thermique entre l'extérieur du réservoir (1) et le gaz, l'échange thermique fourni par un éventuel réchauffeur de vaporisation situé dans le réservoir, la vaporisation de liquide dans le réservoir (1 ), une sixième étape de calcul du bilan d'énergie appliquée au liquide du réservoir (1 ) à partir de l'équation de variation de l'enthalpie du liquide et en prenant en compte au moins l'un des échange suivants : l'échange de chaleur entre l'extérieur et le liquide, l'échange thermique entre le gaz et le liquide, l'échange thermique fourni par un éventuel réchauffeur de vaporisation situé dans le réservoir, la vaporisation de liquide dans le réservoir (1 ), une septième étape de calcul de la température Tg du gaz à partir du bilan d'énergie des cinquième et sixième étapes, et une huitième étape de comparaison de la température T(t+Δt) calculée pour l'instant suivant à la septième étape avec et la température T(t) estimée pour l'instant précédent, et lorsque la différence entre la température T(t+Δt) calculée pour l'instant suivant à la septième étape et la température T(t) estimée pour l'instant précédent est supérieure à un seuil déterminé : une étape de retour à la seconde étape et d'itération jusqu' à convergence, lorsque la différence entre de la température T(t+Δt) calculée pour l'instant suivant à la septième étape et la température T(t) estimée pour l'instant précédent est inférieure à un seuil (convergence) : réitération des étapes ci-dessus pour l'instant (t+2Δt) avec les valeurs de pression mesurées pour cet instant.
EP08840893A 2007-10-26 2008-10-16 Procede de determination en temps reel du niveau de remplissage d'un reservoir cryogenique Withdrawn EP2205899A2 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0758609A FR2922992B1 (fr) 2007-10-26 2007-10-26 Procede de determination en temps reel du niveau de remplissage d'un reservoir cryogenique
PCT/FR2008/051868 WO2009053648A2 (fr) 2007-10-26 2008-10-16 Procede de determination en temps reel du niveau de remplissage d'un reservoir cryogenique

Publications (1)

Publication Number Publication Date
EP2205899A2 true EP2205899A2 (fr) 2010-07-14

Family

ID=39427642

Family Applications (1)

Application Number Title Priority Date Filing Date
EP08840893A Withdrawn EP2205899A2 (fr) 2007-10-26 2008-10-16 Procede de determination en temps reel du niveau de remplissage d'un reservoir cryogenique

Country Status (4)

Country Link
US (1) US8370088B2 (fr)
EP (1) EP2205899A2 (fr)
FR (1) FR2922992B1 (fr)
WO (1) WO2009053648A2 (fr)

Families Citing this family (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2922991B1 (fr) * 2007-10-26 2015-06-26 Air Liquide Procede d'estimation de parametres caracteristiques d'un reservoir cryogenique et notamment de parametres geometriques du reservoir
SE534120C2 (sv) * 2009-09-16 2011-05-03 Firefly Ab Mätutrustning för att fastställa ett förhållande mellan innesluten gasvolym och innesluten vätskevolym i en ackumulatortank
US9869429B2 (en) * 2010-08-25 2018-01-16 Chart Industries, Inc. Bulk cryogenic liquid pressurized dispensing system and method
JP2012145406A (ja) * 2011-01-11 2012-08-02 Toshiba Corp 水位計測系の非凝縮性ガス排出装置
US10119657B2 (en) 2012-06-19 2018-11-06 Shailendra Suman Propane tank continuous monitoring system
US9453611B2 (en) * 2013-06-19 2016-09-27 Shailendra K Suman Propane tank continuous monitoring system
WO2014087251A2 (fr) * 2012-06-20 2014-06-12 Proyectos Y Generadores Libelula, S.A. De C.V. Systèmes et procédés pour production répartie de gaz naturel liquéfié
GB201214969D0 (en) * 2012-08-22 2012-10-03 Airbus Uk Ltd Fuel quantity measurement
FR2998643B1 (fr) * 2012-11-23 2015-11-13 Air Liquide Procede de remplissage d'un reservoir de gaz liquefie
FR2998642B1 (fr) * 2012-11-23 2015-10-30 Air Liquide Procede et dispositif de remplissage d'un reservoir de gaz liquefie
US9147018B2 (en) * 2013-01-10 2015-09-29 General Electric Company Method and system for use in controlling a pressure vessel
DE102013203187A1 (de) * 2013-02-26 2014-09-11 Bayerische Motoren Werke Aktiengesellschaft Verfahren und Vorrichtung zur Ermittlung der Füllmasse eines kryogen gespeicherten Gases in einem Behälter
US9377341B1 (en) * 2013-04-12 2016-06-28 Joe David Watson Electronic liquid level sensing device and gauge for liquid-immersed power transformers, reactors and similar equipment
US20140305201A1 (en) * 2013-04-12 2014-10-16 Joe David Watson Electronic liquid level sensing device and gauge for liquid-immersed power transformers, reactors and similar equipment
WO2015031646A1 (fr) * 2013-08-28 2015-03-05 Worthington Cylinders Corporation Jauge de niveau de liquide pour un cylindre de fluide cryogénique
US11085805B2 (en) 2013-10-30 2021-08-10 S1 Technologies, Inc. System and method for identifying a fuel loss
US11100456B2 (en) 2013-10-30 2021-08-24 S1 Technologies, Inc. System and method for determining volume of fluid in a tank
US12094269B2 (en) 2013-10-30 2024-09-17 S1 Technologies, Inc. System and method for identifying a fuel loss event
US9528872B2 (en) * 2013-10-30 2016-12-27 S1 Technologies, Inc. Redundant data communication system for confirming a fuel event and method therefor
CN103674163B (zh) * 2013-12-18 2016-08-24 中国西电电气股份有限公司 一种油浸式变压器储油系统的油位监测方法及监测系统
DE102014211503A1 (de) * 2014-06-16 2015-12-17 Bayerische Motoren Werke Aktiengesellschaft Überwachungsvorrichtung für einen Drucktank sowie Drucktank
CN104316134B (zh) * 2014-11-06 2017-08-01 山东理工大学 油浸式高压电气设备油位监测方法
DE102015203702A1 (de) * 2015-03-02 2016-09-08 Bayerische Motoren Werke Aktiengesellschaft Verfahren zur Bestimmung der Isolationsgüte eines kryogenen Druckbehälters sowie Kraftfahrzeug mit einem kryogenen Druckbehälter
DE102015206823A1 (de) * 2015-04-15 2016-10-20 Bayerische Motoren Werke Aktiengesellschaft Verfahren zum Betrieb eines Kraftfahrzeuges mit einem kryogenen Druckbehälter
JP6607796B2 (ja) 2016-01-27 2019-11-20 三菱重工業株式会社 液体供給設備の液位検出装置、液体供給設備の液位検出方法および当該液位検出装置を備えた液体供給設備
FR3053432B1 (fr) * 2016-06-30 2019-05-10 Engie Procede et systeme pour calculer en temps reel la quantite d'energie transportee dans une cuve de gaz naturel liquefie pressurisee et non refrigeree.
CN106959143A (zh) * 2017-03-22 2017-07-18 中国科学院上海微系统与信息技术研究所 一种自适应高精度压力式水位监测报警装置及方法
DE102018108214A1 (de) * 2018-04-06 2019-10-10 Samson Ag Tankanordnung und Verfahren zur Füllstandsregelung
CN112584785B (zh) * 2018-07-20 2024-09-24 阿特瑞克尔公司 低温外科系统
DE102018121267A1 (de) * 2018-08-31 2020-03-05 Bayerische Motoren Werke Aktiengesellschaft Verfahren zum Betrieb eines Kraftfahrzeugs sowie Kraftfahrzeug
FR3088406B1 (fr) * 2018-11-12 2021-05-07 Naval Group Procédé de contrôle d'un réservoir cryogénique, réservoir cryogénique et bâtiment sous-marin correspondants
US11009897B2 (en) 2018-12-28 2021-05-18 Rosemount Inc. Remote seal system with improved temperature compensation
CN110081948B (zh) * 2019-05-23 2024-01-30 江西制氧机有限公司 一种真空绝热罐体用差压液位计液相气化取压结构
CN111288296A (zh) * 2020-02-13 2020-06-16 蓝箭航天技术有限公司 一种氮气置换设备的控制装置和方法
CN111692334B (zh) * 2020-06-22 2022-04-15 中国核动力研究设计院 压力容器液位测量系统及方法
CN111928920B (zh) * 2020-07-16 2023-09-22 中广核核电运营有限公司 蒸发器液位检测方法、装置、计算机设备和存储介质
CN115560902A (zh) * 2022-09-23 2023-01-03 四川港通医疗设备集团股份有限公司 一种低温液化气体的计量方法

Family Cites Families (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2765326A (en) 1954-08-17 1956-10-02 Minnesota Mining & Mfg Process of making fluorocarbon carbonyl fluorides
US2765205A (en) 1954-10-12 1956-10-02 Primo J Capella Label dispensing system
FR1481242A (fr) * 1966-04-07 1967-05-19 Gaz De France Procédé et dispositif de mesure de la hauteur d'un liquide contenu dans une cuve ou analogue
US3371534A (en) * 1966-05-10 1968-03-05 Foxboro Co Level sensing apparatus
FR2554230A1 (fr) 1983-10-26 1985-05-03 Air Liquide Procede et appareil pour determiner le poids ou la masse d'un gaz liquefie contenu dans un reservoir
US4602344A (en) * 1984-10-25 1986-07-22 Air Products And Chemicals, Inc. Method and system for measurement of liquid level in a tank
US4736329A (en) * 1984-10-25 1988-04-05 Air Products And Chemicals, Inc. Method and system for measurement of liquid level in a tank
US4788648A (en) * 1984-10-25 1988-11-29 Air Products And Chemicals, Inc. Method and system for measurement of a liquid level in a tank
US4782451A (en) * 1984-11-30 1988-11-01 Union Carbide Corporation Process for maintaining liquid supply
US5333498A (en) * 1992-06-19 1994-08-02 W. L. Walker Co., Inc. Apparatus and method for measuring physical characteristics of a liquid
FR2765326B1 (fr) * 1997-06-26 1999-08-06 Air Liquide Procede de generation d'au moins un signal d'alerte relatif au niveau dans un reservoir de stockage d'un produit et dispositif pour la mise en oeuvre du procede
FR2765205B1 (fr) * 1997-06-26 1999-08-06 Air Liquide Dispositif de surveillance du niveau de stockage d'un produit stocke dans un reservoir
FR2811752B1 (fr) 2000-07-13 2002-10-18 Profroid Procede et dispositif de mesure du volume d'un liquide
US6782339B2 (en) 2000-07-31 2004-08-24 Chart Industries, Inc. Differential pressure gauge for cryogenic fluids which selects a density value based on pressure measurement
US6542848B1 (en) * 2000-07-31 2003-04-01 Chart Inc. Differential pressure gauge for cryogenic fluids
JP2003254811A (ja) * 2002-03-05 2003-09-10 Texcell Kk 差圧式液位測定装置
FR2841963B1 (fr) 2002-07-05 2005-07-01 Air Liquide Procede de regulation en pression d'un reservoir de fluide cryogenique, et reservoir correspondant
US6938635B2 (en) * 2002-07-26 2005-09-06 Exxonmobil Research And Engineering Company Level switch with verification capability
US6889508B2 (en) * 2002-10-02 2005-05-10 The Boc Group, Inc. High pressure CO2 purification and supply system
ES2235646B1 (es) * 2003-12-22 2006-03-16 Ros Roca Indox Equipos E Ingenieria, S.L. Planta movil de regasificacion de gnl.
US6932028B1 (en) * 2004-10-06 2005-08-23 Vogt Power International Inc. Apparatus and method for determining a liquid level in a steam drum
WO2007111052A1 (fr) * 2006-03-28 2007-10-04 Mitsui Mining & Smelting Co., Ltd. Appareil d'identification de fluide et méthode d'identification de fluide
US7539593B2 (en) * 2007-04-27 2009-05-26 Invensys Systems, Inc. Self-validated measurement systems
FR2922991B1 (fr) * 2007-10-26 2015-06-26 Air Liquide Procede d'estimation de parametres caracteristiques d'un reservoir cryogenique et notamment de parametres geometriques du reservoir

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO2009053648A2 *

Also Published As

Publication number Publication date
WO2009053648A2 (fr) 2009-04-30
US8370088B2 (en) 2013-02-05
WO2009053648A3 (fr) 2009-08-27
FR2922992B1 (fr) 2010-04-30
FR2922992A1 (fr) 2009-05-01
US20100241371A1 (en) 2010-09-23

Similar Documents

Publication Publication Date Title
WO2009053648A2 (fr) Procede de determination en temps reel du niveau de remplissage d&#39;un reservoir cryogenique
EP2217890A1 (fr) Procede d&#39;estimation de parametres caracteristiques d&#39;un reservoir cryogenique et notamment de parametres geometriques du reservoir
CA2753588C (fr) Appareil et methode d&#39;estimation du volume et de la masse d&#39;un fluide polyphase conserve a des temperatures cryogeniques
JP4760353B2 (ja) 液体水素タンク残量検知システム
US6595048B1 (en) Accurate cryogenic liquid dispenser
US5616838A (en) Metering apparatus for cryogenic liquids
EP2824378B1 (fr) Procédé de remplissage d&#39;un réservoir de gaz
FR3061285B1 (fr) Systeme et dispositif de determination de la valeur d&#39;un parametre de cuve de stockage de liquide, utilisation dudit dispositif
CN110291325A (zh) 用于确定双壁真空绝热容器的绝热质量的装置和方法
JP6234482B2 (ja) 容器内に極低温で貯蔵されたガスの充填程度を検出するための方法及び装置
EP3390893A1 (fr) Procede et systeme pour calculer en temps reel la duree d&#39;autonomie d&#39;une cuve non refrigeree contenant du gnl
CN111219598B (zh) 一种真空绝热储罐真空度检测方法及装置
FR2472177A1 (fr) Mesureur de niveau pour liquide, notamment pour l&#39;eau de refroidissement d&#39;un reacteur
WO2009074757A1 (fr) Procédé de détermination de la masse de fluide dans un réservoir cryogénique ainsi que du débit massique de fluide consommé
US6571624B1 (en) Low gravity liquid level sensor rake
JP2006160287A (ja) タンクローリ及びその質量管理システム
FR2811752A1 (fr) Procede et dispositif de mesure du volume d&#39;un liquide
JP5712037B2 (ja) 内容量推定装置
JP5656416B2 (ja) 液化ガス供給システム
FR3038052A1 (fr) Procede et dispositif d&#39;estimation de la masse de gaz contenue dans un reservoir de gaz sous pression
EP4166915A1 (fr) Systèmes et procédés de mesure de quantité de carburant
GB2546271A (en) A cylinder for pressurised liquefied gas and a method of calculating the liquid volume
FR3138181A1 (fr) Dispositif et procédé de remplissage de réservoir de gaz sous pression
FR2935112A3 (fr) Systeme et procede d&#39;estimation de la quantite de carburant contenu dans un reservoir de vehicule automobile

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20100526

AK Designated contracting states

Kind code of ref document: A2

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MT NL NO PL PT RO SE SI SK TR

RIN1 Information on inventor provided before grant (corrected)

Inventor name: AMMOURI, FOUAD

Inventor name: DONZEL, ALAIN

Inventor name: MACRON, JONATHAN

Inventor name: BOUTEMY, FLORENCE

DAX Request for extension of the european patent (deleted)
17Q First examination report despatched

Effective date: 20130814

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20170503