METHOD AND APPARATUS FOR DETERMINING SATURATION AND PERMEABILITY OF A ROCK SAMPLE USING MERCURY CAPILLARY PRESSURE CURVES

Abstract: A method and apparatus using mercury capillary pressure curves for determining various characteristics of formations. The mercury capillary pressure curve is amplified and small changes are used to compute residual saturation, permeability, Archie's exponent m and the resistivity factor F of the formation.

Description: The present invention relates to a method and apparatus for determining the characteristics of reservoir rocks and particularly, the characteristics of core samples that are removed from reservoir rocks. One method that has been used for determining the characteristics of reservoir rocks is the measurement of capillary pressure. In particular, the mercury pressure response of the reservoir rock has been measured from which estimates have been made of the permeability of the rock as well as the pore volume. One mercury capillary pressure method and apparatus is described in US patent specification 2,604,779. In this method the rock sample or core is placed in a closed vessel the vessel is evacuated and then partially filled with mercury to cover the core. Mercury is then forced into the vessel at a constant rate while plotting the volume of mercury versus the pressure in the vessel. This provides a plot of the mercury capillary pressure versus volume of the core from which various core characteristics may be determined. For example, total pore volume and an indication of the permeability of the core can be determined. In an article appearing in the Journal of Geophysical Research, Volume 71, No. 12, June 15, 1966, at pp. 2911, entitled "Flow of Fluids Through Porous Mediums", there is described experiments with air/water interfaces moving through simple synthesized permeable specimens. The authors observed that the interface did not move continuously through the synthesized sample but rather in a non-continuous manner from pore to pore. The abrupt movements of the surface resulted in pressure drops in the pressure measured in the sample holder. While the author observed the non-continuous manner in which the interface moved through the synthetic permeable specimen, they did not utilize actual rock samples. Furthermore, they relied upon air/water interfaces. For the results to be useful the interface must be conventional water/oil interface that occurs in reservoir rocks or mercury/air interface. Further, the pores that occur in reservoir rocks are smaller than those present in the authors' synthetic sample by a factor of at least 10 to 100. The equipment and procedures used by the authors would be incapable of measuring the pressure and volume changes occurring in an actual rock sample. It is therefore an object of the invention to provide a method and apparatus suitable for measuring characteristics of rock samples in an easy manner. The invention therefore provides a method for determining the characteristics of a rock sample using mercury capillary pressure data, characterized by the steps of: amplifying the mercury saturation versus mercury capillary pressure data; measuring the volume of each subison S v to obtain the total volume T s of the subison pore systems in the amplified data; measuring the total pore volume P v of the sample; and using the measurements of the total volume of the subison pore systems and total pore volume to calculate the sample characteristics. The invention also provides an apparatus for obtaining a saturation versus pressure curve from a rock sample, characterized by a sample holder, said sample holder being formed from two cup-shaped members, said two cup-shaped members being secured together to form a hollow sample holder, the mating surfaces of said cup-shaped members having ground surfaces to form a leak-tight seal; a piston type pump; drive means, said drive means being coupled to said piston pump; volume measuring means, said volume measuring means being coupled to said pump to digitally measure the volume displaced by said pump; a pressure transducer, said pressure transducer being coupled to said sample cell to measure the pressure in said sample cell and producing an electrical signal related thereto; amplifying means, said amplifying means being coupled to said pressure transducer to amplify said electrical signal; an analog-to-digital converter, said converter being coupled to a pressure transducer to convert said pressure signal to a digital quantity; and recording means, said volume measuring means and said digital converter being coupled to said recorder. The present invention is based on the discovery that by using extremely slow injection rates and making accurate measurements of both the volume and pressure of the displacing fluid that is injected, one can detect the pressure at which the displacing fluid breaks through to each subison pore system in the rock sample. The term "subison" will be explained below. The invention utilizes these measurements of the volumes of subison pore systems to yield the saturation of the non-wetting phase in the rock sample and the residual saturation from a given initial saturation. The data also is utilized to predict the absolute permeability of the sample as well as the value of Archie's exponent m and the formation resistivity factor F used in Archie's equations that have, since 1942, been used to calculate oil and/or water saturations in clean sands. Archie's original empirical relations were modified by Waxman-Smits equations described in "Electrical Conductivity in Oil-Bearing Shaly Sands", Waxman-Smits, 1968 SPE Journal, June, pp. 107-122. Both the Archie and the Waxman-Smits equations require determination of the formation resistivity factor F and Archie's exponent m. The apparatus of the present invention is an improvement over that previously used in that it permits very slow controlled injection rates in combination with means for accurately measuring both the injected volumes and pressures. Further, the sample holder has been modified to provide a leak-tight vessel without the use of gaskets or other sealing means. The elimination of gaskets removes a large source of error in previous equipment used for making capillary pressure measurements. The present invention will now be described by way of example in more detail with reference to the accompanying drawings, in which: Fig. 1 is a typical capillary pressure plot with a portion of the curve enlarged; Fig. 2 is the portion of the curve that is enlarged in Fig. 1 shown to a larger scale; Fig. 3 is a block diagram of an apparatus suitable for performing the present invention. Referring now to Figs. 1 and 2, there is shown a mercury capillary pressure curve and a portion of the curve enlarged. In Fig. 1 the mercury saturation (horizontal axis) increases to the left while the mercury capillary pressure (vertical axis) increases in a "Y" direction. This is a typical curve from a formation rock, for example, of Berea sandstone and shows for low saturations, the pressure increases very slowly but with increasing mercury saturation the pressure increases rapidly. Also shown in Fig. 1 is an enlargement of a small portion of the curve showing the pressure drops 10 and 11. The pressure drops are the result of the mercury passing through the throat of a pore and flowing into the pore space. As the mercury flows into the pore space the pressure decreases rapidly until the pore starts to fill at which time, the pressure then increases at a steady rate to approximately its initial value before breaking through the throat of the pore. The enlarged portion of Fig. 1 is shown to an even larger scale in Fig. 2 where the various portions of the curve are labeled. In particular, the rising portion A of the capillary pressure curve to previously unattained levels is referred to a rison. On the other hand, the rising portion B of the capillary pressure curve to previously attained levels is referred to as a subison, while the decreasing pressure portion indicated by C of the curve is referred to as a rheon. The volume between the point at which the pressure initially decreases and the point at which it again achieves its original value is the volume of a subison pore system indicated by D. A subison pore system consists of only a rheons and subisons. If the horizontal distance of each pore system were summed, one would obtain the volume of the pore system of the sample or more particularly the trapping pore volume. The trapping pore volume is the volume of the pore system that can hold hydrocarbons. From the above description it can be seen that the trapped saturation of the non-wetting phase S nwr fluid after a waterflood can be predicted from the following equation: where V sT = total subison pore system volume and V p = total pore volume of sample. The above relationship holds true for any strongly water-wet system which would be present after a waterflood type of recovery. In this type of system, the water present in the pores will be removed by the vacuum and thus the subison pore system volume that is measured by the method of the present invention will be result of the removal of the residual oil. One could also predict the residual saturation S r if one knows the initial saturation S I by summing the volume of all the given subison pore systems, that would be saturated at an initial saturation of S I . where S v = the volume of a single subison pore system and S I = initial saturation. Similarly the permeability K of the formation can be derived from the following expression: where  = porosity of the sample N r = number of measured risons V ri = volume of a given rison V rT = total volume of risons P r = entry pressure of rison The model equations for the permeability are derived by recognizing that the distribution of risons governs the flow of fluid, i.e. permeability. Each rison can be considered as a capillary tube of a diameter equal to the entry diameter (which can be calculated from the entry pressure) and with a volume equal to the rison volume. Thus, the permeability can be determined from the distribution of risons represented as capillary tubes in parallel over a wide range of different diameters. Likewise, the Archie's exponent m and the formation resistivity factor F can be calculated from the following equations: These equations for m and F are calculated, as with the permeability, from the risons represented as capillary tubes in parallel. The total conductivity is then the sum of the conductivities from each tube. The angular dependence disappears because each rison is assumed to have an isotropic contribution. Thus, the sum of the conductivities amounts by cancellation to just the total rison volume, T r . The average pore volume can be calculated from the following expression: where N s = number of subison pore systems. All of the above data can be obtained from accurate measurements of the volumes and the pressures during a single capillary pressure measurement using the apparatus of this invention. For example, if the core is first evacuated and then mercury forced through the core until it is completely saturated to obtain a complete mercury capillary pressure curve, both the total pore volume as well as the volume of each subison will be obtained. The mercury may be injected into the core at any suitable rate, for example less than 10⁻⁵ cubic centimeters per second. The pressure is measured with an accuracy of better than 1 part in 10⁴, using a 16-bit analog-to-digital converter, while the volume is measured with an accuracy of 1 part in 10⁶ and the measured values converted to a digital quantity, the required calculations can be performed in a small computer or specialized processor. Since the measurements are in digital form, it will be a simple matter for the processor to detect a decrease in pressure indicating a rheon as shown in Fig. 2. The processor can sum the volume until the pressure again reaches the value it had at the time that it entered the pore space to determine the total volume of all the subisons. The remaining calculations can then be readily performed. Referring to Fig. 3, there is shown an instrument suitable for carrying out the above measurements with the accuracy required. In particular, there is shown a sample holder 10 which is formed of two cupshaped members 11 and 12. The mating surfaces 13 and 14 of the cup-shaped surfaces are ground and lapped so that they will form a fluid-tight seal without the addition of any sealing material or gaskets. The presence of gaskets or sealing material will produce faulty results since the voids in gaskets, such as O-rings, or sealing materials can contribute to the volume of the pore spaces that are being measured. In order to effect the seal between the two mating surfaces, the two units should be bolted together with the tension on the bolts being accurately measured as by means of torque wrenches or similar devices. The sample holder is provided with a line 20 and a valve 21 which is coupled to a source suitable vacuum (arrow V). This permits the sample holder and the core material to be completely evacuated before the measurements are started. Prior to evacuating the chamber and the core material the chamber should be filled to a level 16 with a fluid which is to be used to obtain the capillary pressure curve. A pressure intensifier 30 is used for applying pressure to the fluid used in the capillary pressure measurements. The pressure intensifier is, for example, a positive displacement piston type, pump which is driven by actuating means 31. In an advantageous embodiment, the actuating means takes the rotary motion from the speed reducer 32 and converts it to linear motion for driving the piston of the pump unit. The speed reducer in turn is driven by a synchronous motor 33 whose rotation is determined by a shaft/encoder 34. The speed reducer should have a speed reduction of approximately 1000 to 1 when using the conventional 3600 rpm motor and the actuating means 31 should convert this rotary motion into a linear movement of the piston that provides a fluid displacement of approximately 1 × 10⁻⁵cc/sec or less. The pressure in the sample holder is detected by a pressure transducer 40 having an accuracy of better than 1 part in 10⁴ with the signal being amplified by an amplifier 41 at least 1000 times. The amplified signal is supplied to a 16-bit analog-to-digital converter 42 which converts the analog signal to a related digital signal. The digital signal is supplied to both a data acquisition and processing unit 43 as well as to a digital-to-analog converter 44 whose signal is recorded on a chart recorder 45. In particular, the digital-to-analog converter 44 uses only the first 10 bits of the data word from the converter 42 for display on the chart recorder. In particular, the 10-bit signal is displayed full scale on the chart recorder and upon the occurrence of the 11th bit in the data word from the converter 42, the converter 44 resets the chart recorder to 0 so that it can again be used full scale for displaying the first 4 bits of the analog-to-digital converter 42. This, in effect, provides an accuracy of 1 in 65,000 for the data recorder on the recorder 45. The data acquisition processing unit 43 also receives data from the shaft encoder 34 and can make the calculations described above. In particular, the shaft encoder should provide an output signal related to the position of the motor shaft which when correlated with the volume supplied by the pressure intensifier 30 provides an overall accuracy of at least 1 part in 10⁶. This system is operated in much the same manner as the system described in the patent referred to above. In particular, the sample holder is opened and the core material or other sample placed in the sample holder and the sample holder then assembled. The sample holder is partially filled with a displacing fluid, for example mercury, and then evacuated through line 20. After evacuation, the pressure on the mercury is increased by the intensifier 30 while the measurements are made. Once the measurements are obtained, the various characteristics of the core material described above can be readily computed. Various modifications of the invention will become apparent to those skilled in the art from the foregoing description and accompanying drawings. Such modifications are intended to fall within the scope of the appended claims.

Claims: 1. A method for determining the characteristics of a rock sample using mercury capillary pressure data, characterized by the steps of: amplifying the mercury saturation versus mercury capillary pressure data; measuring the volume of each subison Sv to obtain the total volume Ts of the subison pore system in the amplified data; measuring the total pore volume Pv of the sample; and using the measurements of the total volume of the subison pore systems and total pore volume to calculate the sample characteristics. 2. The method as claimed in claim 1 characterized in that the residual saturation of the non-wetting phase Snwr is determined from the expression: wherein VsT is the total subison pore system volume, and Vp is the total pore volume of sample. 3. The method as claimed in claim 1 characterized in that the residual saturation Sr is determined from a given initial saturation SI, from the expression: wherein Sv and Vp are as above defined. 4. The method as claimed in claim 1 characterized in that the permeability k is determined from the expression: wherein  is the porosity of the sample, Vri is the volume of a given rison, Nr the number of measured risons, Pr is the entry pressure of the rison, and VrT the total volume of risons. 5. The method as claimed in claim 1 characterized in that Archie's lithogical exponent m is determined from the expression: wherein , VrT and Vp are as above defined. 6. The method as claimed in claim 1 characterized in that the resistivity factor F is determined from the expression: wherein Vp and VrT are as above defined. 7. The method as claimed in claim 1 characterized in that the mercury is injected into the core at a rate of less than 10⁻⁵ cubic centimeters per second. 8. The method as claimed in claim 7 characterized in that the volume of injected fluid is measured with a precision of at least 1 part in 10⁶ cubic centimeters and the pressure is measured with a precision of at least 1 part in 10⁴. 9. The method as claimed in claim 8 characterized in that the measured pressure is converted to a digital quantity using a 16-bit analog-to-digital converter. 10. The method as claimed in claim 9 characterized in that the pressure is recorded versus time using a chart recorder, the pressure being recorded using only some of the first occurring bits for full scale and resetting the recorder to zero whenever a bit in the next position after said first occurring bit changes value. 11. The method as claimed in claim 10 characterized in that only the first 10 bits are recorded full scale and the recorder is reset to zero when a change occurs in the 11th bit. 12. An apparatus for obtaining a saturation versus pressure curve from a rock sample, characterized by a sample holder, said sample holder being formed from two cup-shaped members, said two cup-shaped members being secured together to form a hollow sample holder, the mating surfaces of said cup-shaped members having ground surfaces to form a leak-tight seal; a piston type pump; drive means, said drive means being coupled to said piston pump; volume measuring means, said volume measuring means being coupled to said pump to digitally measure the volume displaced by said pump; a pressure transducer, said pressure transducer being coupled to said sample cell to measure the pressure in said sample cell and producing an electrical signal related thereto; amplifying means, said amplifying means being coupled to said pressure transducer to amplify said electrical signal; an analog-to-digital converter, said converter being coupled to a pressure transducer to convert said pressure signal to a digital quantity; and recording means, said volume measuring means and said digital converter being coupled to said recorder. 13. The apparatus as claimed in claim 11 characterized in that the drive means includes a motor and lead screw for moving said piston pump, said motor being coupled to a shaft encoder for producing a digital signal related to the position of the piston pump. 14. The apparatus as claimed in claim 13 characterized in that said analog-to-digital converter is a sixteen bit converter. 15. The apparatus as claimed in claim 14 characterized in that said converter is coupled to a chart recorder to record said pressure with respect to volume. 16. The apparatus as claimed in claim 15 characterized in that said chart recorder records the first 10 bits of information in said converter and resets to zero on the occurrence of a change in the 11th bit of said converter.

Inventors: STAPLETON, ROBERT GENE, SWANSON, BENJAMIN FRITHIOF, YUAN, HERBERT HI-HWA
Applicants: SHELL INTERNATIONALE RESEARCH MAATSCHAPPIJ B.V.
Application year: 1986
Main classification: G01N15
IPC codes: G01N15