CROSS-REFERENCE TO RELATED APPLICATIONSThis patent application claims the benefit of U.S. Provisional Patent Application Ser. No. 61/824,854, filed May 17, 2013, entitled “Methods for deriving relative permeability from resistivity measurements;” and U.S. Provisional Patent Application Ser. No. 61/817,160, filed Apr. 29, 2013, entitled “Correlation between resistivity and petrophysical parameters of a carbonate rock during two-phase flow displacements;” both of which are incorporated herein by reference in their entireties.
FIELDThe subject disclosure generally relates to production of oil and gas from subterranean reservoirs. More particularly, the subject disclosure relates to methods for deriving relative permeability from borehole resistivity measurements.
BACKGROUNDIn formation evaluation and reservoir engineering, resistivity index, relative permeability, and capillary pressure are often valuable parameters for estimating oil reserves and planning for production. They can be determined in the laboratory using conventional and Special Core Analysis, or SCAL techniques.
Several theoretical models have been proposed to infer relative permeability from capillary pressure. A few studies have been initiated to correlate relative permeability and/or capillary pressure with resistivity. For example, Pirson et al., found an empirical relationship between relative permeability and resistivity index (Pirson et al., “Prediction of relative permeability characteristics of intergranular reservoir rocks from electrical measurements,” J. Petrol. Technol., (1964), 561-570). Li et al., (Li, K. and Horne, R. N. “Experimental Verification of Methods to Calculate Relative Permeability Using Capillary Pressure Data,” SPE 76757, Proceedings of the 2002 SPE Western Region Meeting/AAPG Pacific Section Joint Meeting, Anchorage, Ak., May 20-22, 2002; Li, K., “A Semi-analytical Method to Calculate Relative Permeability from Resistivity Well Logs,” SPE 95575, SPE Annual Technical Conference and Exhibition, Dallas, USA, 9-12 Oct. 2005; and Li, K., “A New Method for Calculating Two-Phase Relative Permeability from Resistivity Data in Porous Media,” Transport in Porous Media (2007), DOI 10.1007/s11242-007-9178-4) has developed a semi-analytical model to infer relative permeability from resistivity, and confirmed it using experimental data. All of the above references are incorporated herein by reference in their entirety.
SUMMARYThis summary is provided to introduce a selection of concepts that are further described below in the detailed description. This summary is not intended to identify key or essential features of the claimed subject matter, nor is it intended to be used as an aid in limiting the scope of the claimed subject matter.
According to some embodiments a method is described for determining relative permeability for a subterranean rock formation. The method includes: receiving downhole resistivity data representative of resistivity measurements made in a wellbore penetrating the subterranean rock formation during a water flooding event of the subterranean rock formation; and determining one or more relative permeability values for the subterranean rock formation based at least in part on the downhole resistivity data. According to some embodiments, the permeability values include relative permeability for wetting (e.g., water) and non-wetting (e.g., oil) phases during an imbibition mode corresponding to the water flooding event.
According to some embodiments, the resistivity measurements are made using a logging-while-drilling tool. In such cases the aqueous wetting phase can be a water-based drilling mud and the water flooding event can be the introduction of the water-based drilling mud into the rock formation during a drilling process.
According to some embodiments, the water flooding event is caused by water being injected from an injection well. The resistivity measurements can be made from an injection well, an observer well and/or a production well. According to some other embodiments, the water flooding event is caused by a brine invasion from a second rock formation.
According to some embodiments, the relative permeability values are further based on porosity values of the rock formation. The porosity can be derived from, for example borehole measurements such as neutron density measurements, NMR measurements, dielectric measurements, and/or acoustic measurements.
According to some embodiments, the relative permeability values are further based on one or more derived Archie's law parameters, such as saturation exponent and/or cementation factor. The Archie's law parameters can be based on a laboratory core analysis procedure, and/or on borehole measurements.
According to some embodiments, a system is described for determining relative permeability for a subterranean rock formation. The system includes: a downhole resistivity measurement tool configured to be deployed in a borehole penetrating the rock formation and to take resistivity measurements during a water flooding event of the rock formation; and a processing system configured to determine one or more relative permeability values for the subterranean rock formation, based at least in part on downhole measurements made during the water flooding event made by the resistivity tool.
As used herein the term “determining” is to be broadly construed and includes like terms such as deriving, calculating, modeling, obtaining, acquiring, estimating and extracting.
As used herein the term “water” includes aqueous fluids such as brine and water-based drilling mud.
Further features and advantages of the subject disclosure will become more readily apparent from the following detailed description when taken in conjunction with the accompanying drawings.
BRIEF DESCRIPTION OF THE DRAWINGSThe subject disclosure is further described in the detailed description which follows, in reference to the noted plurality of drawings by way of non-limiting examples of the subject disclosure, in which like reference numerals represent similar parts throughout the several views of the drawings, and wherein:
FIG. 1 is a flow chart illustrating aspects of determining relative permeability from drainage-phase laboratory measurements, according to some embodiments;
FIGS. 2-1,2-2 and2-3 depict the experimental results in primary drainage, according to some embodiments;
FIG. 3 is a flow chart illustrating aspects of a method for determining relative permeability for imbibition following a drainage-phase from laboratory measurements, according to some embodiments;
FIGS. 4-1,4-2 and4-3 depict the experimental results in imbibition, according to some embodiments;
FIG. 5 is a flow chart illustrating aspects of a method for determining relative permeability on a log scale based on downhole resistivity measurements, according to some embodiments;
FIG. 6 is a diagram illustrating aspects of systems and methods for determining relative permeability from borehole resistivity measurements, according to some embodiments; and
FIG. 7 illustrates a wellsite system in which the relative permeability can be determined from resistivity and other measurements taken while drilling with a water-based mud, according to some embodiments.
DETAILED DESCRIPTIONThe particulars shown herein are by way of example and for purposes of illustrative discussion of the examples of the subject disclosure only, and are presented in the cause of providing what is believed to be the most useful and readily understood description of the principles and conceptual aspects of the subject disclosure. In this regard, no attempt is made to show structural details in more detail than is necessary, the description taken with the drawings making apparent to those skilled in the art how the several forms of the subject disclosure may be embodied in practice. Furthermore, like reference numbers and designations in the various drawings indicate like elements.
According to some embodiments, methods are disclosed for predicting relative permeability from resistivity measured in the lab using a steady-state method and from direct measurements in the well (i.e., downhole resistivity measurements). The steady-state flooding method measures resistivity and steady-state relative permeability at the equilibrium state.
According to some other embodiments, methods are disclosed to modify and improve Li's model in primary drainage. According to some embodiments, the pore size distribution index λ in Li's model is substituted by the saturation exponent n in the equation of the relative permeability of the wetting phase. This substitution can be justified by the fact that, when brine is the wetting phase, as is the case in water-wet conditions, the mobility of the wetting brine phase is highly affected by the capillary forces, whereas the mobility of the non-wetting oil phase is more governed by the viscous forces. The saturation exponent n accounts for this in the mobility of the wetting phase. According to some embodiments, for the non-wetting phase which will be less affected by the capillary effects (no capillary suction due to the water wettability of the rock), the factor λ is retained in the model.
According to some embodiments, methods are disclosed for determining Kr during an imbibition cycle. Experiments performed have shown good agreement between Kr determined from resistivity and Kr measured in the lab using the steady-state method during an imbibition cycle. The imbibition methods described are believed to more accurately describe what happens in the wellbore during brine invasion (water-based mud while drilling or brine while injecting). According to some embodiments, the above modeling methods for drainage and imbibition are used to determine relative permeability in a reservoir rock using borehole resistivity measurements made in the reservoir rock.
FIG. 1 is a flow chart illustrating aspects of determining relative permeability from drainage-phase laboratory measurements, according to some embodiments. More particularly, the described method is for determining relative permeability from resistivity measured on water-wet core plugs during a drainage steady-state experiment. Inblock110, conventional core analysis is used to determine porosity and absolute permeability. Inblock112, steady-state experiments are conducted with resistivity measurements carried out at equilibrium. Inblock114, a resistivity index (RI) curve is plotted to determine the saturation exponent (n).
Inblock116, the equivalence betweenEquations 1 and 2 below are used to determine λ.
Inblock118, Krw (relative permeability to the wetting phase or to water) and Kro (relative permeability to oil) are calculated using Equations 4 and 5 below.
Krw=(Sw)(2+n)/n (Equation 4)
Krnw=(1−Sw*)2[1−(Sw*)(2+λ)/λ] (Equation 5)
Inblock120, the Kr curves are un-normalized. The curve intervals will be Swi<Sw<1 for the X-axis, 0<Krw<1 and 0<Kro<Kro(Swi) for the Y-axis. At the end of the drainage steady-state experiment, Ko(Swi) is known.
The model using Equations 4 and 5 was validated in the lab on a water-wet rock. In general, for primary drainage, rocks are assumed to be water-wet because they are cleaned before starting the SCAL experiments at Sw=1 (water saturation).
FIGS. 2-1,2-2 and2-3 depict the experimental results in primary drainage, according to some embodiments.FIG. 2-1 depicts the steady-state Kr (relative permeability) withcurve210 showing Krw (water) andcurve212 showing Kro (oil).FIG. 2-2 depicts the RI (resistivity index)curve220 as plotted according to block114 inFIG. 1.FIG. 2-3 depicts the comparison between experimental steady state (SS) Kr (curves230 and232), Kr from Li's model (curves240 and242), and Kr from the methods in primary drainage shown inFIG. 1 (curve250), according to some embodiments. As can be seen by comparingcurves230,240 and250, a significant improvement over Li's model is obtained using the model described, according to some embodiments.
FIG. 3 is a flow chart illustrating aspects of a method for determining relative permeability for imbibition following a drainage-phase from laboratory measurements, according to some embodiments. More particularly, the method shown inFIG. 3 is for the imbibition steady-state cycle, following a primary drainage, to determine relative permeability from resistivity. Inblock310, in continuity with the primary drainage, a steady-state experiment is run with resistivity measurements taken at equilibrium. Inblock312, the RI (resistivity) curve is plotted to determine n (saturation exponent).
Inblock314, the equivalence betweenEquation 6 and 7 below is used to determine λ.
Inblock316, Krw (relative permeability to the wetting phase or to water) and Kro (relative permeability to oil) are calculated using the equations below, and Kr (relative permeability) is plotted on a normalized curve.
Krw=(Sw*)(2+λ)/λ (Equation 10)
Krnw=(1−Sw*)2[1−(Sw*)(2+λ)/λ] (Equation 11)
The Kr curves are un-normalized. Curve intervals will be Swi<Sw<Sor (Swi is irreducible water saturation, Sw is the water saturation and Sor is oil saturation resistivity) for the X-axis, 0<Kr0<1 (Ko(Swi) is used to normalize imbibition Kr) and 0<Krw<Krw(Sor) (Krw is the relative permeability to the wetting phase or to water) for the Y-axis. At the end of the imbibition steady-state experiment, Kw(Sor) is known.
FIGS. 4-1,4-2 and4-3 depict the experimental results in imbibition, according to some embodiments.FIG. 4-1 depicts the steady state Krw (curve410) and Kro (curve412) for the imbibition.FIG. 4-2 depicts theRI curve420.FIG. 4-3 depicts the comparison between experimental steady-state (SS) Krw (curve430) and Kro (curve432), and Kr from the model for imbibition mode described inFIG. 3, supra. As can be seen, an acceptable match is obtained for both wetting and non-wetting phase relative permeability values.
FIG. 5 is a flow chart illustrating aspects of a method for determining relative permeability on a log scale based on downhole resistivity measurements, according to some embodiments. For oil reservoirs, the imbibition model can be used since the reservoir is at Swi (irreducible water saturation). Inblock510, porosity is determined from neutron-density or other logging techniques. Examples include NMR or dielectric logging tools. Inblock512, Rt (resistivity) is determined from an electrical log during a brine invasion or a water injection process. Inblock514, m and n (or other suitable Archie's law parameters) are determined from core analysis or from other logging techniques. Inblock516, Sw (water saturation) is determined from the second Archie's law. Inblock518, λ is determined using the equivalence between Equation 11 above and Equation 12 below:
RIminis the resistivity index at the end of the water flooding (end of imbibition or Sor) when the resistivity reaches its minimal and constant value.
Inblock520, Krw and Kro are calculated using Equation 16 and 17 below, and Kr is plotted on a normalized curve.
Krw=(Sw*)(2+λ)/λ (Equation 16)
Krnw=(1−Sw*)2[1−(Sw*)(2+λ)/λ] (Equation 17)
Note that λ is not substituted by n in Krw since it is an imbibition cycle, and not a drainage cycle.
Inblock522, the Kr curve is un-normalized by determining the Kr end-points using resistivity and formation-tester measurements (See Zeybek, M., Ramakrishnan, T. S., Al-Otaibi, S. S., Salamy, S. P., Kuchuk, F. J., “Estimating Multiphase-Flow Properties from Dual-Packer Formation-Tester Interval Tests and Openhole Array Resistivity Measurements,” SPE Reservoir Eval. Eng., 40-46, February 2004, which is incorporated by reference herein).
FIG. 6 is a diagram illustrating aspects of systems and methods for determining relative permeability from borehole resistivity measurements, according to some embodiments. An injection well610 is used to inject water into a formation600, which is for example a hydrocarbon bearing rock formation. On the surface of injection well610,wellsite612 includes pumping and monitoring equipment for both injecting water and other fluid which can be stored intank614. Also located atwellsite612, according to some embodiments, is surfacedata monitoring unit616 that is in communication with a permanently or semi-permanently installedresistivity measuring unit624. The water injected is via well610 at a packer-isolatedinjection zone618.
The fluid produced from reservoir600 is collected by one or more producer wells, for example well producer well630. According to some embodiments, the resistivity in producer well630 is monitored by either a permanent (or semi-permanent)resistivity measuring unit628 or, according to some embodiments, adownhole resistivity tool636 is deployed viawireline634 andwireline truck632. According to some other embodiments,wireline tool636 is used to retrieve stored data from permanentresistivity measurement unit628. According to yet other embodiments, resistivity measurements are obtained using a permanent orsemi-permanent measurement unit626 that is positioned in anobserver well620. The data fromunit626 is transmitted tosurface station622. Using an observer well, according to some embodiments, is useful in some applications since it is less likely to be effected by “end effects” and so can provide an improved representation of the resistivity changes during the water injection event.
According to some embodiments, the resistivity data fromunits624,626,628 and/ortool636 are transmitted to adata processing unit650. The processing unit includes astorage system642, communications and input/output modules640, auser display646 and auser input system648. According to some embodiments, theprocessing unit650 may be located in thelogging truck632, or at another wellsite location, such as atwellsite612 or withinsurface station622.Data processing unit650 carries out the calculations that facilitate the determinations of relative permeability, such as described with respect to some or all ofFIGS. 1,3 and5 described supra.
According to some embodiments, the resistivity measurements fromunits624,626,628 and/ortool636 are taken during a water (e.g., brine) invasion event rather than a water injection procedure.
FIG. 7 illustrates a wellsite system in which the relative permeability can be determined from resistivity and other measurements taken while drilling with a water-based mud, according to some embodiments. The wellsite can be onshore or offshore. In this system, aborehole711 is formed in subsurface formations by rotary drilling in a manner that is well known. Embodiments of the invention can also use directional drilling, as will be described hereinafter.
Adrill string712 is suspended within theborehole711 and has abottom hole assembly700 that includes adrill bit705 at its lower end. The surface system includes platform andderrick assembly710 positioned over theborehole711, theassembly710 including a rotary table716,kelly717,hook718 androtary swivel719. Thedrill string712 is rotated by the rotary table716, energized by means not shown, which engages thekelly717 at the upper end of the drill string. Thedrill string712 is suspended from ahook718, attached to a traveling block (also not shown), through thekelly717 and arotary swivel719, which permits rotation of the drill string relative to the hook. As is well known, a top drive system could also be used.
In the example of this embodiment, the surface system further includes drilling fluid ormud726, stored in apit727 formed at the well site. Apump729 delivers thedrilling fluid726 to the interior of thedrill string712 via a port in theswivel719, causing the drilling fluid to flow downwardly through thedrill string712, as indicated by thedirectional arrow708. The drilling fluid exits thedrill string712 via ports in thedrill bit705, and then circulates upwardly through the annulus region between the outside of the drill string and the wall of the borehole, as indicated by thedirectional arrows709. In this well-known manner, the drilling fluid lubricates thedrill bit705 and carries formation cuttings up to the surface as it is returned to thepit727 for recirculation.
Thebottom hole assembly700 of the illustrated embodiment contains a logging-while-drilling (LWD)module720, a measuring-while-drilling (MWD)module730, a roto-steerable system and motor, anddrill bit705.
TheLWD module720 is housed in a special type of drill collar, as is known in the art, and can contain one or a plurality of known types of logging tools. It will also be understood that more than one LWD and/or MWD module can be employed, e.g., as represented at720A. (References, throughout, to a module at the position of720 can alternatively mean a module at the position of720A as well.) The LWD module includes capabilities for measuring, processing, and storing information, as well as for communicating with the surface equipment. In the present embodiment, the LWD module includes a resistivity measuring device as well as a number of other devices, such as a neutron-density measuring device.
TheMWD module730 is also housed in a special type of drill collar, as is known in the art, and can contain one or more devices for measuring characteristics of the drill string and drill bit. The MWD tool further includes an apparatus (not shown) for generating electrical power to the downhole system. This may typically include a mud turbine generator powered by the flow of the drilling fluid, it being understood that other power and/or battery systems may be employed. In the present embodiment, the MWD module includes one or more of the following types of measuring devices: a weight-on-bit measuring device, a torque measuring device, a vibration measuring device, a shock measuring device, a stick slip measuring device, a direction measuring device, and an inclination measuring device.
According to some embodiments, theresistivity data760 and other data such as neutron density data, is transmitted todata processing unit650, which can be located a the wellsite or in some remote location. According to some embodiments, measurements made during the water invasion in the form of water-based mud from the drilling process is used to determine relative permeability as described herein (e.g., with respect toFIG. 5, supra).
In view of the above description it will be appreciated that features of the subject disclosure may be implemented in computer programs stored on a computer readable medium and run by processors, application specific integrated circuits and other hardware. Moreover, the computer programs and hardware may be distributed across devices including but not limited to tooling which is inserted into the borehole and equipment which is located at the surface, whether onsite or elsewhere.
Although only a few examples have been described in detail above, those skilled in the art will readily appreciate that many modifications are possible in the examples without materially departing from this subject disclosure. Accordingly, all such modifications are intended to be included within the scope of this disclosure as defined in the following claims. In the claims, means-plus-function clauses are intended to cover the structures described herein as performing the recited function and not only structural equivalents, but also equivalent structures. Thus, although a nail and a screw may not be structural equivalents in that a nail employs a cylindrical surface to secure wooden parts together, whereas a screw employs a helical surface, in the environment of fastening wooden parts, a nail and a screw may be equivalent structures. It is the express intention of the applicant not to invoke 35 U.S.C. §112, paragraph 6 for any limitations of any of the claims herein, except for those in which the claim expressly uses the words ‘means for’ together with an associated function.