EP1155342A1 - Determination de la saturation en eau et d'une fraction de sable a partir d'un outil d'imagerie de la resistivite des puits de forage, d'une diagraphie par induction transversale et d'un modele tensoriel de saturation en eau - Google Patents

Determination de la saturation en eau et d'une fraction de sable a partir d'un outil d'imagerie de la resistivite des puits de forage, d'une diagraphie par induction transversale et d'un modele tensoriel de saturation en eau

Info

Publication number
EP1155342A1
EP1155342A1 EP99968973A EP99968973A EP1155342A1 EP 1155342 A1 EP1155342 A1 EP 1155342A1 EP 99968973 A EP99968973 A EP 99968973A EP 99968973 A EP99968973 A EP 99968973A EP 1155342 A1 EP1155342 A1 EP 1155342A1
Authority
EP
European Patent Office
Prior art keywords
resistivity
shale
formation
sand
tool
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
EP99968973A
Other languages
German (de)
English (en)
Other versions
EP1155342A4 (fr
Inventor
Richard A. Mollison
Juergen H. Schoen
Otto N. Fanini
Berthold F. Kriegshauser
Milomir Pavlovic
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.)
Baker Hughes Holdings LLC
Original Assignee
Baker Hughes Inc
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 Baker Hughes Inc filed Critical Baker Hughes Inc
Priority to EP08170797.8A priority Critical patent/EP2026105A3/fr
Publication of EP1155342A1 publication Critical patent/EP1155342A1/fr
Publication of EP1155342A4 publication Critical patent/EP1155342A4/fr
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/26Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
    • G01V3/28Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Definitions

  • the invention is related generally to the field of interpretation of measurements made by well logging instruments for the purpose of determining the fluid content of earth formations. More specifically, the invention is related to methods for calculating fractional volumes of various fluids disposed in the pore spaces of earth formations where these earth formations include laminations of shale with reservoir rock that may include dispersed shales.
  • a significant number of hydrocarbon reservoirs include deep water turbidite deposits that consist of thin bedded, laminated sands and shales.
  • a common method for evaluating the hydrocarbon content of reservoirs is the use of resistivity measurements.
  • typically one or more types of porosity-related measurement will be combined with measurements of the electrical resistivity (or its inverse, electrical conductivity) of the earth formations to infer the fluid content within the pore spaces of the earth formations.
  • the fractional volumes of connate water and hydrocarbons can be inferred from empirical relationships of formation resistivity Rt with respect to porosity and connate water resistivity such as, for example, the well known Archie relationship.
  • the Archie relationship fractional volume of water in the pore space is represented, as shown in the following expression, by Sw - known as "water saturation”: w R, R t ⁇ " 1 ( l )
  • a and m are empirically determined factors which relate the porosity (represented by ⁇ ) to the resistivity of the porous rock formation when it is completely water-saturated (Ro), R w represents the resistivity of the connate water disposed in the pore spaces of the formation, and m represents an empirically determined "cementation" exponent, n is the saturation exponent.
  • Relationships such as the Archie formula shown in equation (1) do not work very well when the particular earth formation being analyzed includes some amount of extremely fme-grained, clay mineral-based components known in the art as "shale". Shale typically occurs, among other ways, in earth formations as "dispersed” shale, where particles of clay minerals occupy some of the pore spaces in the hydrocarbon- bearing earth formations, or as laminations (layers) of clay mineral-based rock interleaved with layers of reservoir-type rock in a particular earth formation.
  • fractional volume of pore space which is capable of containing movable (producible) hydrocarbons.
  • the fractional volume of such formations which is occupied by dispersed shale can be estimated using such well logging devices as natural gamma ray radiation detectors. See for example, M. H. Waxman et al, "Electrical Conductivities in Oil Bearing Shaly Sands", SPE Journal, vol. 8, no. 2, Society of Petroleum Engineers, Richardson, TX (1968).
  • the layers sometimes are thick enough to be within the vertical resolution of, and therefore are determinable by, well logging instruments such as a natural gamma ray detector.
  • the shale layers are determined not to be reservoir rock formation and are generally ignored for purposes of determining hydrocarbon content of the particular earth formation.
  • a problem in laminated shale reservoirs is where the shale laminations are not thick enough to be fully determined using gamma ray detectors and are not thick enough to have their electrical resistivity accurately determined by electrical resistivity measuring devices known in the art.
  • Sands that have high hydrocarbon saturation are typically more resistive than shales.
  • conventional induction logging tools greatly underestimate the resistivity of the reservoir: the currents induced in the formation by the logging tool flow preferentially through the conductive shale layers leading to an overestimate of the conductivity of the formation.
  • R t represents the electrical resistivity (inverse of conductivity) in the reservoir rock layers of the formation and R ⁇ represents the resistivity in the shale layers.
  • the reservoir exhibits an anisotropy in the resistivity.
  • This anisotropy may be detected by using a logging tool that has, in addition to the usual transmitter coil and receiver coil aligned along with the axis of the borehole, a receiver or a transmitter coil aligned at an angle to the borehole axis.
  • a logging tool that has, in addition to the usual transmitter coil and receiver coil aligned along with the axis of the borehole, a receiver or a transmitter coil aligned at an angle to the borehole axis.
  • Such devices have been well described in the past for dip determination. See, for example, United States patent 3,510,757 to Huston and United States patent 5,115,198 to Gianzero,
  • United States patent 5,656,930 issued to Hagiwara discloses a method of determining the horizontal resistivity, the vertical resistivity, and the anisotropy coefficient of a subterranean formation by means of an induction type logging tool positioned in a deviated borehole within the subterranean formation.
  • the induction type logging tool is first calibrated to determine a proportionality constant.
  • a predetermined relationship between the proportionality constant, the phase shift derived resistivity, the attenuation derived resistivity, the horizontal resistivity, the vertical resistivity, and the anisotropy coefficient is then generated and stored in the memory of a programmed central processing unit.
  • the phase shift derived resistivity and attenuation derived resistivity are then received and processed by the programmed central processing unit in accordance with the predetermined relationship to generate the horizontal resistivity, the vertical resistivity, and the anisotropy coefficient.
  • These measured values of horizontal and vertical resistivities when combined with a predetermined relationship between the horizontal resistivity, the vertical resistivity, the net/gross ratio, and the ratio of the sand layer resistivity to the shale layer resistivity make it possible to obtain a net/gross ratio.
  • the sands may include dispersed shales. Interpretation of formation water saturation in such reservoirs can be in error if the combined effects of laminations, dispersed shales within the sand, and possible intrinsic anisotropy of the shales is not considered.
  • the present invention is method of accounting for the distribution of shale in a reservoir including laminated shaly sands using vertical and horizontal conductivities derived from multi-component induction data.
  • data may also be acquired using a borehole resistivity imaging tool.
  • the data from the borehole resistivity imaging tool give measurements of the dip angle of the reservoir, and the resistivity and thickness of the layers on a fine scale.
  • the measurements made by the borehole resistivity imaging tool are calibrated with the data from the induction logging tool that gives measurements having a lower resolution than the borehole resistivity imaging tool.
  • a tensor petrophysical model determines the laminar shale volume and laminar sand conductivity from vertical and horizontal conductivities derived from the log data.
  • the volume of dispersed shale and the total and effective porosities of the laminar sand fraction are determined using a Thomas-
  • Stieber-Juhasz approach Removal of laminar shale conductivity and porosity effects reduces the laminated shaly sand problem to a single dispersed shaly sand model to which the Waxman-Smits equation can be applied.
  • FIG. 1 shows a resistivity imaging tool suspended in a borehole
  • FIG. 2 (PRIOR ART) is a mechanical schematic view of the imaging tool of Fig. 1;
  • FIG. 2 A is a detail view of an electrode pad for the tool of Figs. 1, 2;
  • FIG. 3 (PRIOR ART) is a pictorial view of a composite imaging log obtained by merging the resistivity image data shown in acoustic image data;
  • FIG. 4 is a flow chart illustrating the principal steps of the process of the invention.
  • FIG. 5 gives the steps of one subprocess of an embodiment of the invention for determination of water saturation from measured values of vertical an horizontal resistivities.
  • FIG. 6 is a schematic illustration of the components of the tensor petrophysical model of the present invention.
  • Fig. 4 is a schematic flowchart of the major steps of the process used in the present invention.
  • a borehole resistivity imaging tool 102 such as is described in United States Patent 5,502,686 issued to Dory et al., and the contents of which are fully incorporated here by reference. It should be noted that the Dory patent is an example of a device that can be used for obtaining measurements borehole resistivity measurements: any other suitable device could also be used.
  • Fig. 1 shows a composite imaging tool 10 suspended in a borehole 12, that penetrates earth formations such as 13, from a suitable cable 14 that passes over a sheave 16 mounted on drilling rig 18.
  • the cable 14 includes a stress member and seven conductors for transmitting commands to the tool and for receiving data back from the tool as well as power for the tool.
  • the tool 10 is raised and lowered by draw works 20.
  • Electronic module 22, on the surface 23, transmits the required operating commands downhole and in return, receives digital data back which may be recorded on an archival storage medium of any desired type for concurrent or later processing.
  • a data processor 24, such as a suitable computer may be provided for performing data analysis in the field in real time or the recorded data may be sent to a processing center or both for post processing of the data.
  • Fig. 2 is a schematic external view of the unified borehole sidewall imager system. This may be used to provide the data that may be used in an optional embodiment of the invention.
  • the tool 10 comprising the imager system includes four important components: 1) resistivity arrays 26; 2) electronics modules 28 and 38; 3) a mud cell 30; and 4) a circumferential acoustic televiewer 32. All of the components are mounted on a mandrel 34 in a conventional well-known manner. The outer diameter of the assembly is about 5.4 inches and about five feet long.
  • An orientation module 36 including a magnetometer and an inertial guidance system is mounted above the imaging assemblies 26 and 32.
  • the upper portion 38 of the tool 10 contains a telemetry module for sampling, digitizing and transmission of the data samples from the various components uphole to surface electronics 22 in a conventional manner.
  • the acoustic data are digitized although in an alternate arrangement, the data may be retained in analog form for transmission to the surface where it is later digitized by surface electronics 22.
  • each array includes 32 electrodes or buttons identified as 39 that are mounted on a pad such as 40 in four rows of eight electrodes each. Because of design considerations, the respective rows preferably are staggered as shown, to improve the spatial resolution. For reasons of clarity, less than eight buttons are shown in Fig. 2A.
  • each pad can be no more than about 4.0 inches wide. The pads are secured to extendable arms such as 42. Hydraulic or spring- loaded caliper-arm actuators (not shown) of any well-known type extend the pads and their electrodes against the borehole sidewall for resistivity measurements.
  • the extendable caliper arms 42 provide the actual measurement of the borehole diameter as is well known in the art.
  • the voltage drop and current flow is measured between a common electrode on the tool and the respective electrodes on each array to furnish a measure of the resistivity of the sidewall (or its inverse, conductivity) as a function of azimuth.
  • the acoustic imager that forms the circumferential borehole imaging system 32 provides 360 ° sampling of the sidewall acoustic reflectivity data from which a continuous acoustic imaging log or sonogram can be constructed to provide a display of the imaged data.
  • the borehole resistivity imaging tool arrays necessarily allow sampling only across preselected angular segments of the borehole sidewall. From those data, a resistivity imaging log, consisting of data strips, one strip per array, separated by gaps, can be constructed and displayed. The angular width of each data-scan strip is equal to 2 sin " ' (S/(2R) ⁇ , where S is the array width and R is the borehole radius.
  • the common data from the two imagers are merged together in a data processing operation to provide a substantially seamless display as shown in Fig. 3.
  • the merging incorporates equalizing the dynamic range of the resistivity measurements with respect to the acoustic measurements. That balance is essential in order that the continuity of a displayed textural feature is not distorted when scanning across a resistivity segment of the display, between adjacent acoustic segments.
  • the display in Fig. 3 incorporates measurements from directional sensors to align the resistivity measurements with geographical coordinates (North, East, South, West), with the resistivity image being "unfolded” to provide a flat image of the cylindrical surface of the borehole.
  • geographical coordinates North, East, South, West
  • the resistivity image being "unfolded” to provide a flat image of the cylindrical surface of the borehole.
  • the display in Fig. 3 shows many such sinusoids, some corresponding to bedding planes and others corresponding to fractures.
  • the dip angle and the dip direction corresponding to the various sinusoids are determined in the present invention using known methods.
  • the sinusoids have essentially zero amplitude.
  • the resistivity measurements are averaged circumferentially and vertically within each identified layer to give an average resistivity measurement for each layer identified above.
  • the subsurface may be characterized by a number of plane layers, each of which has a constant resistivity. With the resolution of the button-electrode tool, these layers may range in thickness from a few millimeters to a few centimeters.
  • the averaging described above is limited to electrodes in the strike direction: these measurements would be more likely representative of the true formation resistivity at the depth of measurement.
  • the resistivity measurements obtained by the averaging process correspond to layers that are beyond the resolution of electromagnetic induction logging tools or propagation resistivity tools. Accordingly, the resistivity measurements obtained at this point are averaged to give resistivities on a scale that would be measurable by an induction logging tool. This is depicted by 104 in Fig. 5.
  • a finely laminated sequence of layers having different resistivities exhibits a transverse isotropy on a larger scale where the wavelength of the electromagnetic wave is much greater than the layer thickness.
  • This condition is easily satisfied even for propagation resistivity tools that, e.g., operate at a frequency of 2MHz (with a wavelength ⁇ ⁇ 6 meters); for induction logging tools that have frequencies of the order of 50kHz to 200kHz, the wavelengths are even longer.
  • the layered medium is characterized by a horizontal resistivity R h * and a vertical resistivity R v * given by:
  • W is a window used to average the resistivities
  • Ah is the depth sampling interval of the electrodes
  • R is the measured resistivity for a given depth.
  • the terms “horizontal” and “vertical” are to be understood in terms of reference to the bedding planes and the anisotropy axes of the subsurface formations, i.e., “horizontal” refers to parallel to the bedding plane, and “vertical” refers to vertical to the bedding plane.
  • the anisotropy axis is taken to be the normal to the bedding plane.
  • data from the orientation module 36 in Fig. 1 may be used to correct the resistivity measurements made by the resistivity imaging tool to give measurements parallel to and perpendicular to the bedding planes.
  • the resistivity measurements made by the electrode-pad system described above may be in error and, in particular, may need to have a scaling factor applied to the data.
  • this data When this data is acquired, it may be calibrated by relating the values given by equations (3) and (4) to data from an induction logging tool or a propagation resistivity tool.
  • an induction or wave propagation tool is used to make measurements of the vertical and horizontal resistivity of the earth formations 104.
  • United States Patent 5781436 to Forgang et al discloses a method an apparatus for making measurements of horizontal and vertical resistivities of a transversely isotropic formation.
  • the method disclosed by Forgang et al comprises selectively passing an alternating current through transmitter coils inserted into the wellbore.
  • Each of the transmitter coils has a magnetic moment direction different from the magnetic moment direction of the other ones of the transmitter coils.
  • the alternating current includes a first and a second frequency.
  • the amplitude at the first frequency has a predetermined relationship to the amplitude at the second frequency. The relationship corresponds to the first and the second frequencies.
  • the method includes selectively receiving voltages induced in a receiver coil having a sensitive direction substantially parallel to the axis of the corresponding transmitter coil through which the alternating current is passed.
  • a difference in magnitudes between a component of the received voltage at the first frequency and a component of the voltage at the second frequency is measured, and conductivity is calculated from the difference in magnitudes of the components of the received voltage at the two frequencies.
  • Rosthal U.S. Patent 5,329,448 discloses a method for determining the horizontal and vertical conductivities from a propagation logging device. The method assumes that ⁇ , the angle between the borehole axis and the normal to the bedding plane, is known. Conductivity estimates are obtained by two methods. The first method measures the attenuation of the amplitude of the received signal between two receivers and derives a first estimate of conductivity from this attenuation. The second method measures the phase difference between the received signals at two receivers and derives a second estimate of conductivity from this phase shift.
  • Two estimates are used to give the starting estimate of a conductivity model and based on this model, an attenuation and a phase shift for the two receivers are calculated.
  • An iterative scheme is then used to update the initial conductivity model until a good match is obtained between the model output and the actual measured attenuation and phase shift.
  • model is generated of the axial distribution of the horizontal and vertical conductivities, from induction signals acquired by the instrument using two-frequency alternating current.
  • the model is generated by calculating an initial estimate of the conductivity distribution and axially inverting the estimate with respect to the measurements made using the two-frequency alternating current. Shoulder correction is applied to measurements made by the instrument using single- frequency alternating current.
  • An estimate of the radial distribution of the conductivities is generated from the shoulder corrected induction signals acquired using the single- frequency alternating current.
  • a 2-dimensional model is made of the conductivity distribution from the model of axial distribution and from the estimate of radial distribution.
  • the initial model for the inversion is based at least in part on data acquired by the resistivity imaging tool 105.
  • the resistivity imaging tool may need to be normalized in some way to correct the resistivity measurements
  • the layer boundaries determined by the resistivity imaging tool serve as a good starting point for the layers used in the inversion of the transverse induction logging tool data.
  • the initial layers for the model may be determined from other high resolution logging tools, such as a Laterolog TM or a gamma ray logging too.
  • the two-frequency induction signals are corrected for near wellbore effects using two-frequency whole space responses calculated using the 2-dimensional model.
  • the corrected two- frequency signals are then axially inverted to generate a 2-dimensional model.
  • all the previous steps are repeated until differences between the corrected two-frequency induction signals from successive repetitions (iterations) of the steps fall below a predetermined threshold.
  • the two-dimensional model extant when process is halted becomes the final two-dimensional model.
  • the horizontal and vertical resistivities obtained therefrom are analyzed using a petrophysical model 108.
  • the values of vertical and horizontal resistivity thus obtained are related to the fluid content and fractional volume of pore spaces in subsurface layers by expressions such as the following derived from the Patchett- Herrick water saturation model for shaly sand formations:
  • R, (i - R sh gives the horizontal resistivity in the reservoir-rock (non-shale) layers of the formation.
  • F * sd h in equation (5) represents the formation resistivity factor for the horizontal resistivity
  • B- Q v is a factor related to the resistivity of "dispersed" shale (shale located within the pore spaces of the reservoir rock).
  • V sh represents the fractional volume within the earth formation of interest of the layers of shale ("laminated shale volume"). The other terms represent the same quantities as described in the Background section herein. See for example, J. G. Patchett et al, "Introduction Section III.
  • the porosity can be determined by any one of a number of well known measurements, such as acoustic travel time, neutron porosity, bulk density, or combinations of measurements such as these as is well known in the art.
  • the porosity measurements just described are meant only as examples of measurements used to determine the porosity and are not meant to limit the invention in any way.
  • equation (5) is written in a form relating to conductivity (inverse of resistivity) rather than in a form related to resistivity because as is known in the art, the signal measured by an induction logging instrument, where eddy currents are induced substantially along layer perpendicular to the wellbore, is related in magnitude to the sum of the conductivities of the individual layers.
  • this can be thought of as current passing through a set of resistors connected in parallel.
  • equation (6) is expressed in terms of resistivity, because where eddy currents are induced in a direction perpendicular to the layers, the effect of layering on the magnitude of the induction signal is similar to passing electrical current through a set of resistors connected in series.
  • ⁇ sh represents an "anisotropy factor" relating the vertical and horizontal conductivities (or resistivities) in the same formation.
  • Equation (9) can be readily solved for S w to provide a calculation of the water saturation (and its complement, the hydrocarbon saturation) in the non-shale layers which does not require explicit determination of the resistivity (or conductivity) of the shale layers in the reservoir earth formation of interest.
  • the C, terms in equation (9) represent total conductivity (electrical conductivity of both the shale and reservoir rock portions of the earth formation of interest).
  • Fig. 5 a flow chart of the subprocess for determination of water saturation according to one embodiment of the invention using the Patchett- Herrick model is depicted.
  • a logging tool such as disclosed in the '436 patent measurements are made within a borehole.
  • the logging tool makes measurements of induction signals along and perpendicular to the axis of the instrument as well as cross-component signals. As described above, these measurements are processed to give a "horizontal' 'and “vertical" resistivity 220 at each depth in the borehole.
  • a measurement of the connate water resistivity R w is obtained 240.
  • Another embodiment of the invention uses an orthogonal tensor model based on electrical anisotropy (R v / R h ) instead of one based on single scalar parallel conductivity models.
  • the tensor model is easily implemented for isotropic and anisotropic shales with isotropic sands.
  • the laminar shale volume must be determined from some external model such as Thomas-Stieber or image log data.
  • True laminar sand porosity must also be derived from the Thomas-Stieber (1975) model and is essential to true laminar reservoir characterization.
  • the tensor model uses the Waxman-Smits equation to evaluate the laminar sand component and utilizes the Hill, Shirley, and Klein (1979) equation to derive Q v from the dispersed clay bound water fraction.
  • the final result of this two- step resistivity tensor model is consistent with the scalar model originally proposed by Patchett and Herrick.
  • the laminar sand resistivity component is derived directly from the tensor model and is implicitly linked to the correct laminar shale volume.
  • Figure 6 illustrates a laminated formation 300 comprising interbedded sand 301 and shale layers 302.
  • the sand may further comprise a clean Archie type sand 303a having poor sorting (i.e., different grain sizes) , a fine laminated sand 303b that includes laminae of well sorted sand with differences in grain size between the individual laminae, and a dispersed shaly sand 303c having shale dispersed therein.
  • the tensor model used herein focuses primarily on the determination of laminar sand properties and hydrocarbon saturation calculation of the interbedded sand layers.
  • the model assumes macroscopic anisotropy because the dimension of laminae is less than the vertical resolution of the measurement. Intrinsic anisotropy within the individual sand laminae is also considered in a general case.
  • the tensor model is two dimensional and only 'vertical' anisotropy is considered. Laminae or beds are assumed to be horizontally or 'laterally' isotropic.
  • the tensor model is implemented in a two-step process.
  • the first step of analysis permits the decoupling of the determination of laminar shale volume, laminae resistivities, and laminar sand properties including porosity, dispersed clay volume, Qv, etc.
  • the determination of saturation is made in a second step in order to allow the application of different sets of assumptions or models to the data
  • the second step of the analysis calculates the laminar sand component hydrocarbon saturation.
  • Sand properties can be described as a clean, homogeneous Archie (1942) sand or as a dispersed shaly sand using the Waxman-Smits (1968) equation.
  • V sh _ Volume fraction of the laminated shale 1-
  • V sh _ ⁇ /G Volume fraction of the sand
  • the forward modeling is derived from normalized parameters by the defined dimensionless anisotropy parameters: ⁇ , ⁇ sd , ⁇ sh
  • the sand normalized conductivities are:
  • the first step in using the tensor petrophysical model is the determination of laminar shale volume and sand conductivity from the measured anisotropy data and total conductivity in the ho ⁇ zontal and vertical directions, Q A and v .
  • the sand and shale conductivities in both horizontal and vertical directions (C sdh , C sd v , C sh h , and C sh v ) are calculated from these data.
  • the volumet ⁇ c sand and shale parameters i.e., V sh and N/G, must be derived. Depending upon the available data and any additional information about the sand and shale properties, three cases are considered.
  • CASE B Shale is anisotropic and the sand is isotropic.
  • the sand conductivity C sd and the laminated shale content V sM can be calculated from the composite formation conductivities, C t h and v , given the two shale conductivities, C sh h , and C sh v :
  • the second step in the application of the tensor model is the analysis of the laminar sand reservoir component.
  • Water saturation of the laminar sand is a function of the 'true' laminar sand porosity, electrical properties, and laminar sand conductivity and can be calculated using various published relationships or other models such as effective medium or electrical efficiency.
  • the Waxman- Smits equation is applied to quantitatively correct for dispersed clay conductivity in the sand. This equation reduces to Archie's equation when dispersed clay is not present.

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Electromagnetism (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

L'invention concerne un procédé permettant de déterminer (22) la porosité totale d'une formation (13), un volume fractionnel du schiste dans ladite formation (13) ainsi que la résistivité dudit schiste, dans un réservoir stratifié comprenant des sables dans lesquels du schiste peut être dispersé. Un modèle pétrophysique tensoriel permet de déterminer le volume de schiste laminaire et la conductivité du sable laminaire à partir des conductivités verticale et horizontale déduites des données (10) multi-composant de la diagraphie par induction. Le volume du schiste dispersé et les porosités totale et effective de la fraction de sable laminaire sont déterminés à l'aide de la technique Thomas-Stieber-Juhasz. L'élimination de la conductivité du schiste laminaire et des effets de la porosité permet de ramener les problèmes posés par le sable schisteux laminaire à un seul modèle de sable schisteux dispersé auquel on peut appliquer l'équation de Waxman-Smits.
EP99968973A 1998-12-30 1999-12-29 Determination de la saturation en eau et d'une fraction de sable a partir d'un outil d'imagerie de la resistivite des puits de forage, d'une diagraphie par induction transversale et d'un modele tensoriel de saturation en eau Withdrawn EP1155342A4 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
EP08170797.8A EP2026105A3 (fr) 1998-12-30 1999-12-29 Détermination de la saturation en eau et d'une fraction de sable à partir d'un outil d'imagerie de la résistivité des puits de forage, d'une diagraphie par induction transversale et d'un modèle tensoriel de saturation en eau

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US22296798A 1998-12-30 1998-12-30
US222967 1998-12-30
US16094399P 1999-10-22 1999-10-22
US160943P 1999-10-22
PCT/US1999/031104 WO2000039612A1 (fr) 1998-12-30 1999-12-29 Determination de la saturation en eau et d'une fraction de sable a partir d'un outil d'imagerie de la resistivite des puits de forage, d'une diagraphie par induction transversale et d'un modele tensoriel de saturation en eau

Related Child Applications (1)

Application Number Title Priority Date Filing Date
EP08170797.8A Division EP2026105A3 (fr) 1998-12-30 1999-12-29 Détermination de la saturation en eau et d'une fraction de sable à partir d'un outil d'imagerie de la résistivité des puits de forage, d'une diagraphie par induction transversale et d'un modèle tensoriel de saturation en eau

Publications (2)

Publication Number Publication Date
EP1155342A1 true EP1155342A1 (fr) 2001-11-21
EP1155342A4 EP1155342A4 (fr) 2006-05-24

Family

ID=46149834

Family Applications (1)

Application Number Title Priority Date Filing Date
EP99968973A Withdrawn EP1155342A4 (fr) 1998-12-30 1999-12-29 Determination de la saturation en eau et d'une fraction de sable a partir d'un outil d'imagerie de la resistivite des puits de forage, d'une diagraphie par induction transversale et d'un modele tensoriel de saturation en eau

Country Status (3)

Country Link
EP (1) EP1155342A4 (fr)
BR (1) BR9916714B1 (fr)
WO (1) WO2000039612A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104020276A (zh) * 2014-06-24 2014-09-03 中国石油大学(北京) 横观各向同性页岩储层岩石力学参数的确定方法

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6686736B2 (en) 2000-08-30 2004-02-03 Baker Hughes Incorporated Combined characterization and inversion of reservoir parameters from nuclear, NMR and resistivity measurements
US7363164B2 (en) * 2004-12-20 2008-04-22 Schlumberger Technology Corporation Method of evaluating fluid saturation characteristics in a geological formation
CN101012748B (zh) * 2007-02-02 2011-08-10 中国石化集团胜利石油管理局测井公司 微电阻率扫描井壁成像测井仪
US7983845B2 (en) 2008-05-15 2011-07-19 Schlumberger Technology Corporation Method and system for analyzing a laminated sand/shale formation
US9037415B2 (en) 2011-04-14 2015-05-19 Baker Hughes Incorporated OBM resistivity image enhancement using principal component analysis with first moment estimation
AU2013315927B2 (en) * 2012-09-13 2017-09-21 Chevron U.S.A. Inc. System and method for performing simultaneous petrophysical analysis of composition and texture of rock formations
US11008835B2 (en) 2015-09-03 2021-05-18 Halliburton Energy Services, Inc. Horizontal reservoir description systems
US10301935B2 (en) * 2016-10-18 2019-05-28 Halliburton Energy Services, Inc. MCI logging for processing downhole measurements
CN107797154B (zh) * 2017-09-22 2019-04-12 中国石油天然气股份有限公司 一种电成像测井图像刻度方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5355088A (en) * 1991-04-16 1994-10-11 Schlumberger Technology Corporation Method and apparatus for determining parameters of a transition zone of a formation traversed by a wellbore and generating a more accurate output record medium
US5671136A (en) * 1995-12-11 1997-09-23 Willhoit, Jr.; Louis E. Process for seismic imaging measurement and evaluation of three-dimensional subterranean common-impedance objects

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
No further relevant documents disclosed *
See also references of WO0039612A1 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104020276A (zh) * 2014-06-24 2014-09-03 中国石油大学(北京) 横观各向同性页岩储层岩石力学参数的确定方法

Also Published As

Publication number Publication date
BR9916714B1 (pt) 2012-08-07
WO2000039612A1 (fr) 2000-07-06
BR9916714A (pt) 2001-12-04
EP1155342A4 (fr) 2006-05-24

Similar Documents

Publication Publication Date Title
US6493632B1 (en) Water saturation and sand fraction determination from borehole resistivity imaging tool, transverse induction logging and a tensorial water saturation model
US6711502B2 (en) Water saturation and sand fraction determination from borehole resistivity imaging tool, transverse induction logging and a tensorial water saturation model
US6686736B2 (en) Combined characterization and inversion of reservoir parameters from nuclear, NMR and resistivity measurements
US7359800B2 (en) Determination of fracture orientation and length using multi-component and multi-array induction data
US7619540B2 (en) Apparatus and methods for determining isotropic and anisotropic formation resistivity in the presence of invasion
US6502036B2 (en) 2-D inversion of multi-component induction logging data to resolve anisotropic resistivity structure
US6466872B1 (en) Method for determination of apparent resistivities of anisotropic reservoirs
US7778778B2 (en) Correction of multi-component measurements for tool eccentricity in deviated wells
US6400148B1 (en) Use of redundant data for log quality measurements
US7505851B2 (en) Use of multi-component measurements in delineating geology of deep-water sediments
US7629791B2 (en) Method and apparatus for making multi-component measurements in deviated wells
US20040133351A1 (en) Method for resistivity anisotropy determination in near vertical wells
WO2011130587A2 (fr) Procédé et appareil pour déterminer une dépression structurelle géologique en utilisant des mesures d'induction multiaxiale
WO2000039612A1 (fr) Determination de la saturation en eau et d'une fraction de sable a partir d'un outil d'imagerie de la resistivite des puits de forage, d'une diagraphie par induction transversale et d'un modele tensoriel de saturation en eau
RU2242029C2 (ru) Определение водонасыщенности и фракции пласта песка с использованием инструмента формирования изображения удельного сопротивления в буровой скважине, инструмента поперечного индукционного каротажа и тензорной модели водонасыщенности

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: 20010725

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE CH CY DE DK ES FI FR GB GR IE IT LI LU MC NL PT SE

AX Request for extension of the european patent

Free format text: AL;LT;LV;MK;RO;SI

RIN1 Information on inventor provided before grant (corrected)

Inventor name: PAVLOVIC, MILOMIR

Inventor name: KRIEGSHAUSER, BERTHOLD, F.

Inventor name: FANINI, OTTO, N.

Inventor name: SCHOEN, JUERGEN, H.

Inventor name: MOLLISON, RICHARD, A.

RBV Designated contracting states (corrected)

Designated state(s): GB IT

REG Reference to a national code

Ref country code: DE

Ref legal event code: 8566

A4 Supplementary search report drawn up and despatched

Effective date: 20060406

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: BAKER HUGHES INCORPORATED

17Q First examination report despatched

Effective date: 20080610

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: 20160701