Abstract:
Method allowing to perform at least an evaluation of the inlet capillary pressure value of a porous medium. From a displacement experiment on a sample of the medium, a curve of the volume of saturating fluid expelled from the sample as a function of time is drawn. The differential pressure between the inlet face and at least one point located at a distance Li from the inlet face that is greater than the distance between the inlet face and the interface between the two fluids within the sample is then continuously measured as a function of time. At least one motive pressure gradient of the first fluid is thereafter calculated by means of distance Li and of the curve. Finally, at least one value of the inlet capillary pressure is determined by calculating the difference between the differential pressure and the value of the motive pressure gradient of the expelled fluid. Application: notably production of oil reservoirs for example.

Description:
FIELD OF THE INVENTION  
       [0001]     The present invention relates to a method allowing to determine at least one value for the inlet capillary pressure of a porous medium.  
         [0002]     The method can notably be applied under operating conditions representative of geologic formations, as regards the nature of the fluids as well as the thermodynamic pressure and temperature conditions.  
         [0003]     The potential applications of this method relate in particular to the characterization of low-permeability porous media such as the cap rocks of an underground reservoir, within the scope of the evaluation of geologic formations as storage sites for fluids such as hydrocarbons, CO 2  or other fluids.  
       BACKGROUND OF THE INVENTION  
       [0004]     The following documents, mentioned in the course of the description hereafter, illustrate the state of the art:  
         [0005]     Chiquet O., Broseta D. and Thibeau S.,            Capillary alteration of shaly caprocks by carbon dioxide           , SPE 94183, 14 th  Europec Conference, Madrid, Spain, 13-16 Jun. 2005,  
         [0006]     Hildenbrand, A., Schlömer S., and Krooss B.,            Gas breakthrough experiments on Fine - grained sedimentary rocks , Geofluids Vol. 2, 3-23, 2002,  
         [0007]     Monicard R.,            Caractéristiques des roches réservoir—Analyse des carottes , Paris, France, Éditions Technip, 1981,  
         [0008]     Zweigel P. et al.,            Towards a methodology for top seal efficacy assessment for underground CO 2  storage , 7 th  International Conference on Greenhouse Gas Control Technologies, Vancouver, 5-9 Sep. 2004.  
         [0009]     A geologic formation capable of keeping fluids (hydrocarbons for oil reservoirs, CO 2  or other gases for storage sites, . . . ) consists of a reservoir rock allowing to collect fluids coming from a source (mother rock or injection) and of an impervious cap rock located above (at the top) of the reservoir and allowing to prevent migration of the fluids from the reservoir to the surface. This geologic formation is then referred to as geologic trap.  
         [0010]     The capacity of a geologic formation to store fluids such as hydrocarbons or CO 2  for example mainly depends on the morphology of the geologic trap and on the petrophysical properties of the rocks that make up the cap layer.  
         [0011]     In general terms, the morphology of a geologic trap is evaluated by means of a geologic characterization based on geophysical data (seismic data for example), and also on data from the wells drilled in the zone (logs, cores and drill cuttings analysis, . . . ). If the geologic formation selected is not located in too hilly a zone, this type of study generally allows to precisely determine the shape and the extent of the geologic trap.  
         [0012]     Petrophysical characterization of the cap rocks requires specific laboratory experiments that can be very long considering the low permeability of such media (typically below 10 −4  mD). It can be reminded that the permeability of a porous medium corresponds to its capacity to allow a fluid (liquid or gas) to flow under the effect of a pressure gradient. Among all the petrophysical properties, it is by far the inlet capillary pressure that plays the most important part in the capacity of cap rocks to maintain the fluids in the reservoir since it controls the allowable maximum storage overpressure at the top of the reservoir. In the literature, the inlet capillary pressure is also referred to as threshold capillary pressure or breakthrough pressure.  
         [0013]     What is referred to as inlet capillary pressure is the minimum pressure difference to be imposed between a non-wetting phase and a wetting phase for the non-wetting phase to be able to start saturating the porous medium considered. The inlet capillary pressure is denoted by P c   E .  
         [0014]     The significance of the value of the inlet capillary pressure is illustrated on the reservoir scale within the scope of CO 2  storage: we consider the case of an aquifer wherein CO 2  is injected so as to reduce the emissions discharged to the atmosphere. During storage, the injected CO 2  whose density is under usual thermodynamic conditions lower than that of the water in place progressively forms a pocket located in the upper part of the reservoir. At the lower boundary of the pocket (water/CO 2  interface), the pressures in each one of the water and CO 2  phases are equal since the capillary pressure curve of the reservoir rocks does not generally exhibit a high inlet capillary pressure because these rocks have the property of readily accommodating fluids. Each one of the two phases having a different density, the pressure gradient is also different, which leads to the existence of a pressure difference between the two phases as shown in  FIG. 1 . This pressure difference is directly related to a height h measured above the water/CO 2  interface.  FIG. 1  illustrates the pressure difference for two heights h 1  and h 2  of the water/CO 2  interface. The abscissas represent pressures P, P CO2  for the CO 2  and P W  for the pressure in the water phase. The ordinates represent height h above the interface. In the first case, the interface is at a depth point C 11 . At depth point C x , height h above the interface is h 1 . The difference between the pressures is P C (h 1 ). In the second case, the interface is at a depth point C 12 . At depth point C x , height h above the interface is h 2 . The difference between the pressures is P C (h 2 ). Within the scope of a CO 2  storage operation, height h above the interface changes from h 2  to h 1 .  FIG. 1  thus shows an increase in the pressure difference imposed on the top with time, within the scope of a CO 2  storage operation.  
         [0015]     In general terms, we write: P CO2 (h)−P W (h)=(ρ W −ρ CO2 )gh=P c (h) with: 
        h: the height above the water/CO 2  interface     P CO2  (h): the pressure in the CO 2  phase for a height h     P W (h): the pressure in the water phase for a height h     ρ CO2 : density of the CO 2       ρP W : density of the water     g: gravity     P c (h): the capillary pressure corresponding to a height h.        
 
         [0023]     This pressure difference directly corresponds to the notion of capillary pressure. This capillary pressure P c (h) furthermore controls the value of the saturation for a given height h. P c (h) increases as a function of h as shown in  FIG. 1 : the greater h, the higher the CO 2  pressure. The maximum capillary pressure in the reservoir is thus reached at the reservoir top (h=H). This value is denoted by P c   T =P c (H). Since capillary continuity is provided at the reservoir/cap rocks interface, the value of the capillary pressure at the top P c   T  is therefore also imposed on the cap rock. Two cases can then arise:  
         [0000]     the capillary pressure at the reservoir top is lower than the inlet capillary pressure (P c   T &lt;P c   E ): the CO 2  remains confined;  
         [0000]     the capillary pressure at the reservoir top is higher than the inlet capillary pressure (P c   T &gt;P c   E ): the CO 2  starts circulating in the cap rock and the water saturation in the cap layer will tend towards the value corresponding to P c   T .  
         [0024]     In practice, within the scope of fluid injections in an underground reservoir, it is advisable to take a given margin in relation to the value of the capillary pressure at the reservoir top calculated from H because the injection itself can generate dynamic overpressures that can locally lead to capillary pressures at the top that are higher than the calculated capillary pressure at the top (P c   T ).  
         [0025]     The previous reminder shows how important it is to properly evaluate the value of the inlet capillary pressure of a porous medium, for example within the scope of the storage of fluids such as hydrocarbons, CO 2  or other fluids in geologic traps.  
         [0026]     There are various methods for evaluating the inlet capillary pressure P c   E  of a porous medium for the storage conditions (thermodynamic conditions and nature of the fluids).  
         [0027]     There is for example a known technique based on the mercury porosimetry method. This method consists in converting a capillary pressure curve obtained by mercury porosimetry for the reservoir conditions by means of the following conversion formula (Monicard, 1981):  
           P   c   E     ⁡     (   s   )       =         P   c   E     ⁡     (   m   )       ⁢         σ   s     ⁢   cos   ⁢           ⁢     θ   s           σ   m     ⁢   cos   ⁢           ⁢     θ   m               
 
 with: 
        σ m : mercury/air interfacial tension=480 mN/m     θ m : mercury/air contact angle=140°    σ s : interfacial tension for the fluids considered in the reservoir (typically CO 2 /brine within the scope of CO 2  storage)     θ s : contact angle for the fluids considered in the reservoir (typically CO 2 /brine within the scope of CO 2  storage)     P c   E (s): value of the inlet capillary pressure under storage conditions     P c   E (m): value of the inlet capillary pressure for the mercury porosimetry measurement conditions under ambient conditions.        
 
         [0034]     Although this method allows very fast estimation of a value for the inlet capillary pressure under storage conditions, the representativity thereof can be affected because of the uncertainty as regards the wettability phenomena (contact angle θ storage ). Since the value of the contact angle is generally not known, it is selected equal to zero, which corresponds to a perfect wettability of the fluid in place. Recent experimental measurements have shown that this hypothesis could turn out to be erroneous in particular in the case of geologic storage of CO 2  (Chiquet et al., 2005). This approach can lead to significant errors in the calculation of the inlet capillary pressure.  
         [0035]     There is also another approach, referred to as “conventional” approach, whose principle is based on the very definition of the inlet capillary pressure. This method is for example described by Monicard (1981).  
         [0036]     During this approach, the sample to be studied is first saturated and placed in a containment cell, which allows to work under imposed pressure and temperature conditions. The inlet end piece of the cell is then swept so as to bring the non-wetting fluid, such as CO 2  for example, for which the inlet capillary pressure is sought, just in contact with the face of the sample. A device allowing to measure small liquid productions (either by weighing the liquids produced or by direct measurement from a fine capillary tube) is then set at the level of the outlet end piece.  
         [0037]     The experiment then consists in increasing the pressure of the non-wetting fluid at the inlet face in successive increasing stages while monitoring the production level of the fluid saturating the sample at the outlet. The value of the inlet capillary pressure of the rock in relation to the two fluids used then corresponds to the imposed pressure for which production start of the fluid in place has been observed.  
         [0038]     Although the principle is simple, implementation of this type of experiment is however delicate within the scope of cap rock evaluation, for the following reasons:  
         [0039]     Length: the number of stages before the desired threshold is reached can be large since, in many cases, no realistic approximations are available before the experiment is started. On the other hand, the required waiting time for each pressure plateau is generally rather long to allow effective detection of the production at the outlet;  
         [0040]     Accuracy: in the vicinity of the inlet capillary pressure, the non-wetting fluid invasion kinetics is particularly slow because its mobility threshold is then reached, which makes its outlet flow rate extremely low and therefore difficult to detect, all the more so since the rock studied is of low permeability.  
         [0041]     The conventional approach thus leads to very long experiment times and rather to an overestimation of the inlet capillary pressure because of a wrong detection of the mobility threshold.  
         [0042]     Another technique is the residual capillary pressure approach. This approach was recently proposed (Hildenbrand, 2002) to provide a faster experimental alternative in relation to the conventional approach. The porous medium test cell is prepared in the same way as for the conventional method, but it is placed between two cells C 1  and C 2  arranged on either side and containing the non-wetting fluid. A valve initially separates cell C 1  from the sample which is however in contact with cell C 2 . A pressure P 1  is initially imposed in C 1  and a lower pressure P 2  is imposed in C 2  by seeing to it that the difference between the two pressures is greater than the estimated value of the capillary pressure sought.  
         [0043]     The valve is then opened while recording the evolution of P 1  and P 2 . A progressive decrease, then a stabilization is observed for P 1  in the course of time, which corresponds to a circulation of the non-wetting fluid towards C 2  through the sample. Similarly, pressure P 2  increases, then stabilizes. A residual differential pressure is experimentally observed between the two cells, which is interpreted as the inlet capillary pressure of the rock in relation to the fluids studied. Within the scope of this approach, it is also possible to keep P 1  constant and to monitor the evolution of P 2  only in the course of time.  
         [0044]     This approach aroused a real fad as it was published because of the rapidity thereof and of the ease of interpretation. However, recent work has shown that it is risky to directly interpret the differential pressure measured at the end of the experiment directly in terms of inlet capillary pressure (Zweigel et al., 2005). In fact, during the experiment, the upstream part of the sample undergoes an initial drainage stage during circulation of the wetting fluid, then an imbibition stage as the differential pressure progressively decreases. The measured residual pressure therefore corresponds to a pressure at the end of an imbibition stage and not to a pressure at the start of a drainage stage like the inlet capillary pressure. Now, many experimental findings show that these two pressures are generally not equal, the differential pressure at the end of the imbibition stage being systematically lower than the inlet capillary pressure. Although interpretation of the results of this method a priori affords many advantages, it leads, as it is currently considered, to a systematic underestimation of the inlet capillary pressure value.  
         [0045]     In relation to the methods currently used and described above, the method according to the invention allows to obtain a result rapidly, as regards acquisition of the necessary experimental data as well as their interpretation in terms of inlet capillary pressure.  
       SUMMARY OF THE INVENTION  
       [0046]     The invention relates to a method allowing to perform at least an evaluation of the inlet capillary pressure value of a porous medium, from a displacement experiment wherein a first fluid saturating a sample of said medium is caused to flow by injection of a second fluid at the level of a face of the sample referred to as “inlet” face.  
         [0047]     The method comprises the following stages:  
         [0000]     applying to the sample a constant pressure allowing said second fluid to flow into said sample,  
         [0000]     drawing a curve of the volume of the first fluid expelled from said sample as a function of time,  
         [0048]     continuously measuring as a function of time a local differential pressure DP t   i  between said inlet face and at least one point located at a distance L i  from the inlet face, this distance being greater than the distance between said inlet face and an interface between the two fluids within the sample,  
         [0000]     calculating at least one differential pressure of the first fluid DP w   i  by means of distance L i  and of said curve,  
         [0000]     determining at least one inlet capillary pressure value from local differential pressure DP t    i  and the differential pressure value DP w   i  of the first fluid.  
         [0049]     The inlet capillary pressure value can be defined as the difference between local differential pressure DP w   i  and the differential pressure value DP w   i  of the first fluid.  
         [0050]     According to the method, the differential pressure value DP w   i  of the first fluid can be calculated after a slope change of said curve and by determining this new slope.  
         [0051]     According to an embodiment, distance L i  is equal to the length of the sample and only the total differential pressure DP t  between said inlet face and an opposite face of the sample is then measured.  
         [0052]     According to another embodiment, distance L i  is smaller than the length of the sample, and the position of the interfaces between the two fluids can be determined by means of local saturation measurements along the sample. Several inlet capillary pressure values can thus be determined and an uncertainty on the value of this inlet capillary pressure can be deduced therefrom.  
         [0053]     According to the invention, it may be wise to leave a volume of first fluid upstream from the sample prior to starting injection of the second fluid. Furthermore, the total imposed differential pressure, the temperature of the fluids and the nature of the fluids can allow to restore conditions representative of oil reservoirs or of fluid storage sites.  
         [0054]     Within the scope of gas storage for example, and notably within the scope of CO 2  storage, the first fluid can be water and the second fluid a gas.  
         [0055]     Finally, the method can be applied to a porous medium of reservoir cap type so as to evaluate the storage capacity of a storage reservoir intended for a gaseous fluid, a hydrocarbon or CO 2  for example. 
     
    
     BRIEF DESCRIPTION OF THE FIGURES  
       [0056]     Other features and advantages of the method according to the invention will be clear from reading the description hereafter of embodiments given by way of non limitative example, with reference to the accompanying figures wherein:  
         [0057]      FIG. 1  illustrates the increase with time of the capillary pressure imposed at the top within the scope of a CO 2  storage operation,  
         [0058]      FIG. 2  shows the experimental device used to evaluate the dynamic inlet capillary pressures according to the invention,  
         [0059]      FIG. 3  shows a production curve obtained within the scope of the method proposed and underlines the flow rate change observed when the non-wetting phase starts saturating the sample (decrease in the slope of the production curve, therefore in the flow rate). 
     
    
     DETAILED DESCRIPTION  
       [0060]     The methods used in the profession and described above are essentially based on a “static” determination of the inlet capillary pressure value (semi-static for the residual capillary pressure method since we go through a transient flow phase), i.e. a fluid is injected, then stabilization of the outlet flow rate is monitored. From a more “dynamic” point of view, the inlet capillary pressure can also be considered to be a differential pressure between the two phases that does not contribute to the flow. Consider a sample initially saturated with a wetting fluid that is caused to flow by injecting a non-wetting fluid with a total differential pressure on the sample that is constant and equal to DP t . The pressure profile can then be split up into several parts according to the nature of the phases present in the various parts of the sample:
 
 DP   t   =P   g   am   =P   w   av   =P   g   am   −P   g   fr   +P   g   fr   −P   w   fr   +P   w   fr   −P   w   av 
 
 with: 
        DP t : the total differential pressure imposed on the sample     P g   am : the gas pressure upstream from the sample (inlet face)     P w   av : the water pressure downstream from the sample (outlet face)     P g   fr : the gas pressure at the front (gas/water interface)     P w   fr : the water pressure at the front (gas/water interface)        
 
         [0066]     In the case considered above, the total differential pressure imposed on the sample can thus be split up into three terms:
 
 DP   t   =DP   g   +P   c   fr   +DP   w 
 
 with: 
        DP g : the differential pressure (pressure drop) in the zone invaded by the gas     DP w : the differential pressure (pressure drop) in the zone that is not invaded by the gas     P c   fr : the capillary pressure drop at the front which corresponds to the inlet capillary pressure: P c   fr =P c   E .        
 
         [0070]     We can therefore write:
 
 DP   t   =DP   g   +P   c   E   +DP   w 
 
         [0071]     Consider the very beginning of the invasion of the medium by the gas (non-wetting fluid). We can then assume that:  
         [0000]     the differential pressure in the gas zone (DP g ) is negligible considering the limited extent of this zone (and also considering the low viscosity of the fluid injected, as it is the case for CO 2  even under storage conditions),  
         [0072]     the differential pressure in the water (DP w ) corresponds to the effective pressure difference that leads to the production of liquid in the buret at the outlet. This differential pressure DP w  can therefore be directly calculated from Darcy&#39;s law for the production flow rate measured at the outlet (Q w ) and a permeability value (K) measured otherwise by means of techniques known to specialists:  
         Q   w     =             K   ·   S       μ   w       ·       DP   w     L       ⇒     DP   w       =           μ   w     ·   L       K   ·   S       ⁢     Q   w             
 
 with: 
        S: the section of the sample (known)     L: the length of the sample (known)     μ w : the dynamic viscosity of the water (known).        
 
         [0076]     The method according to the invention thus allows to evaluate the inlet capillary pressure P c   E  by means of the expression as follows:
 
 P   c   E   =DP   t   −DP   w 
 
         [0077]     It is thus possible to determine the value of the inlet capillary pressure P c   E  by injecting directly the non-wetting fluid (gas) and by measuring the effective flow rate of the wetting fluid at the outlet (Q w ), which allows to calculate the differential pressure in this phase DP w . The inlet capillary pressure (P c   E ) is then simply obtained by calculating the difference between the total differential pressure imposed on the sample (DP t ) and the value of the differential pressure in the wetting phase (DP w ).  
         [0078]     Thus, the method according to the invention can be broken down into five major stages.  
         [0079]     1—A non-wetting fluid such as gas is directly injected into a sample of section S and of length L, saturated with a wetting fluid such as water, and constrained to a total differential pressure on the sample DP t .  
         [0080]     2—The total differential pressure imposed on the sample (DP t ) is measured and a curve representing the volume of the wetting fluid expelled at the sample outlet as a function of time is drawn. This curve is a line whose slope corresponds to the reference flow rate (Q ref ).  
         [0081]     3—The effective flow rate of the wetting fluid at the sample outlet (Q w ), which corresponds to the new slope of the curve, is calculated from the previous curve and after a slope change thereof.  
         [0082]     4—The differential pressure in the wetting phase DP w  is calculated, for example from Darcy&#39;s law and by means of the dimensions of the sample (L and S), the dynamic viscosity μ w , the permeability of the sample (K) and the effective flow rate of the wetting fluid at the sample outlet (Q w ).  
         [0083]     5—The value of the inlet capillary pressure P c   E  is determined by calculating the difference between the total differential pressure imposed on the sample (DP t ) and the value of the differential pressure in the wetting phase (DP w ): P c   E =DP t −DP w .  
         [0084]     Experimental implementation  
         [0085]     A sample of the porous medium whose inlet capillary pressure is to be evaluated is first saturated with a wetting fluid, then placed in a test cell known to the man skilled in the art and described for example in patent FR-2,708,742 (U.S. Pat. No. 5,679,885). What is understood to be a porous medium is any medium having pores through which a fluid can flow. The porosity can therefore have any value.  
         [0086]     This device allows multi-flow rate displacement experiments to be carried out on a sample (E) from a porous medium. This device is diagrammatically shown in  FIG. 2 . It comprises an elongate containment cell  1  that can be cylindrical, and which contains the sample to be tested E between two end pieces  7   a  and  7   b.  This cell is placed inside a thermostat-controlled enclosure  2  (not shown) so as to subject the sample to be tested to a predetermined temperature T. At the outlet of one of the faces of cell  1 , an outlet end piece  3  allows the expelled fluid to be sent to a system  4  allowing to determine the volume of wetting fluid expelled. At the opposite face, the device comprises an inlet end piece  5  allowing injection of a non-wetting fluid into sample (E). At the level of containment cell  1 , the sample to be tested E is placed within a deformable sheath  6  and the assembly is arranged within enclosure  2 . The annular space around sheath  6  is communicated with a source of fluid under pressure (not shown) so as to subject the sample to be tested to a predetermined pressure P. The device also includes at least one differential pressure detector CDP t  between the inlet and the outlet of the sample, for measuring the differential pressure produced on the sample itself DP t . The device can comprise several such detectors, referred to as CDP 1 , CDP 2 , . . . and arranged at a respective distance from the inlet face L 1 , L 2 , . . . .  
         [0087]     The experimental conditions are representative of the storage conditions insofar as the thermodynamic conditions, the nature of the fluids used and the state of mechanical stress (confining pressure in the test cell) are respected.  
         [0088]     Without local saturation and pressure measurement  
         [0089]     According to a first embodiment, the total pressure is measured without carrying out local pressure and saturation measurements.  
         [0090]     Unlike the “conventional” method for determining the inlet capillary pressure, a significant volume of wetting fluid, typically some cubic centimeters, can be kept in inlet end piece  5 . The non-wetting fluid is then injected at an imposed pressure so that the total differential pressure DP t  is higher than an “a priori” estimated inlet capillary pressure value. The wetting fluid saturating the sample is then caused to flow with two distinct periods as shown in  FIG. 3 :  
         [0091]     Period  1 : as long as the non-wetting fluid is located in the inlet end piece, the wetting fluid circulates with a motive gradient directly related to the imposed total differential pressure DP t . A linear production curve whose slope corresponds to the flow rate calculated with Darcy&#39;s law for the DP t  and which thus represents the reference flow rate (Q ref ) for the imposed total differential pressure during the test is thus obtained. The flow is a single-phase flow throughout the sample;  
         [0092]     Period P 2 : as soon as the non-wetting fluid reaches the inlet face, the differential pressure in the wetting fluid phase decreases by the value of the inlet capillary pressure and a decrease in the production slope at the outlet is immediately observed. This slope then corresponds to the effective flow rate of the wetting fluid at the sample outlet (Q w ) allowing to calculate the differential pressure in the wetting phase DP w .  
         [0093]      FIG. 3  shows a curve of the volume of fluid produced (V) as a function of time in hour (T) obtained within the scope of the method proposed. Two distinct slopes can be clearly seen, which correspond first to a strictly single-phase flow in the sample (P 1 ), then to the arrival (A) of the non-wetting fluid at the inlet face, producing a reduction in the slope (P 2 ) and therefore in the production flow rate (Q w ).  
         [0094]     The presence of the volume of wetting fluid initially located upstream from the sample is not a prerequisite within the scope of the interpretation but it contributes to the quality control of the test carried out since it allows to establish the reference flow rate (Q ref ) just before the non-wetting fluid starts flowing through the sample.  
         [0095]     It can furthermore be noted that if the value of DP t  was selected in a too conservative way (typically below the inlet capillary pressure value), the non-wetting fluid will stop at the inlet face without compromising the success of the test. In this case, the value of DP t  just has to be significantly increased to carry the test through to completion.  
         [0096]     With local saturation and pressure measurement  
         [0097]     According to another embodiment, a measurement of the total differential pressure and local pressure and saturation measurements are carried out. If the test cell is provided with an instrumentation allowing local pressure measurements to be performed along the sample and with a device allowing the saturation profile to be measured, it is possible to carry on with the degree of interpretation. This device is also shown in  FIG. 2 .  
         [0098]     In fact, as long as the gas front located by means of the saturation measurement, by X-rays for example (RX), is upstream from the local pressure measuring points, each one of these points can be used to apply the previous approach “locally”. The following formulas are then used, i representing the number of the local pressure detector used and L i  the distance from the inlet face and the position of the detector considered:  
         Q   w     =         K   ·   S       μ   w       ⁢       DP   w   i       L   i             
         P   c   E     =       DP   t   i     -     DP   w   i           
 
         [0099]     This approach thus allows, with a single experimentation, to obtain several values for the inlet capillary pressure of a sample, which improves determination of this parameter while giving an error range for the result obtained. This is particularly useful within the scope of larger-scale studies intended to evaluate the associated uncertainties and risks as regards a storage site.  
         [0100]     The table hereafter allows to compare the three methods mentioned for three different rocks (R1, R2 and R3):  
                                                                                                         Porosity   P c   E     P c   E     P c   E             Type   K (mD)   (%)   conventional   residual   dynamic                                    R1   Chalk   1.7   40   0.9   0.2   0.8       R2   Carbonate   0.016   14.5   6   3.1   6.2       R3   Sandstone   0.0014   13   10   7.7   9.6                  
 
         [0101]     It can be seen that the method according to the invention provides results that are comparable to those of the conventional method. These results are however more accurate and they were obtained more rapidly than those of the conventional method.