CA2401898A1 - Method of interpreting induction logs in horizontal wells - Google Patents
Method of interpreting induction logs in horizontal wells Download PDFInfo
- Publication number
- CA2401898A1 CA2401898A1 CA002401898A CA2401898A CA2401898A1 CA 2401898 A1 CA2401898 A1 CA 2401898A1 CA 002401898 A CA002401898 A CA 002401898A CA 2401898 A CA2401898 A CA 2401898A CA 2401898 A1 CA2401898 A1 CA 2401898A1
- Authority
- CA
- Canada
- Prior art keywords
- signals
- properties
- subset
- zone
- borehole
- 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.)
- Abandoned
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/26—Electric 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/28—Electric 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Electromagnetism (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
An induction logging tool (100) is used in an inclined borehole (10, 110, 410) for determining properties of subsurface formations away from the borehole.
Measurements are made at a plurality of transmitter-receiver (T-R) (106-108F) distances. After correction (202) of the data for skin effects and optionally correcting (204) for eddy currents within the borehole, the shallow measurements (those from short T-R spacing or from high frequency data) are inverted to give a model (206) of the near borehole (invaded zone resistivity and diameter) and the resistivity of the formation outside the invaded zone.
Using this model, a prediction (208) is made of the data measured by the mid-level and deep sensors (long T-R spacings).
Measurements are made at a plurality of transmitter-receiver (T-R) (106-108F) distances. After correction (202) of the data for skin effects and optionally correcting (204) for eddy currents within the borehole, the shallow measurements (those from short T-R spacing or from high frequency data) are inverted to give a model (206) of the near borehole (invaded zone resistivity and diameter) and the resistivity of the formation outside the invaded zone.
Using this model, a prediction (208) is made of the data measured by the mid-level and deep sensors (long T-R spacings).
Description
FIELD OF~THE INVENTION
The invention is related to the field of induction logging of subsurface earth formations. Specifically, it relates to the use of induction logging for obtaining structural information of the subsurface away from a borehole.
BACKGROUND OF THE INVENTION
to Induction resistivity well logging devices are l~nown in the art for measuring the electrical resistivity of earth formations penetrated by wellbores. Electrical resistivity measurements are used for, among other things, to infer the fluid content of pore spaces of the earth formations.
15 Induction electrical resistivity instruments l~nown in the art, such as one described in Hunka, J. "A New Resistivity Measurement System for Deep Formation Imaging and High Resolution Formation Evaluation", paper no. 20559, Society of Petroleum Engineers, Richardson, Tex., 1990, typically measure voltages induced in receiver coils positioned at spaced apart locations along the instrument. The voltages are induced by magnetic fields 20 generated by eddy currents flowing in the earth formations around the instrument, in response to alternating currents of various frequencies being passed through a transmitter coil. The voltages induced in a particular one of the receiver coils are dependent on the electrical conductivity (which is inversely related to resistivity) of the media, including the earth formations, surrounding the instrument, and are dependent on the spacing of the particular receiver coil with respect to the transmitter coil. Generally, the greater the spacing between the transmitter coil and the particular receiver coil, the greater is the radial depth from the wellbore into the earth formation from which the measurement at that receiver coil corresponds. The vertical resolution of the measurements made by the closer-spaced receiver coils, however, is proportionately finer than is the vertical resolution of the measurements made by the to more distant receiver coils.
When the wellbore is drilled, fluids from the drilling process may be forced into the pore spaces of some of the earth formations, changing their fluid content and therefore their resistivity. The process of fluid being forced into the pore spaces is 15 generally referred to as "invasion". Wellbore can be drilled with fluids having differing electrical resistivities, and the wellbore itself can have widely varying diameters over its length. The effects of variable wellbore diameter and differing wellbore fluid resistivity can indeterminately affect the total voltage being induced in each one of the receiver coils, because eddy currents also can flow within the 2o conductive wellbore fluid and within that portion of the earth formations surrounding the wellbore subject to the invasion (the so-called "invaded zone") United States Patent 5,703,773 ( referred hereafter as the " '773 patent" ) having the same assignee as the present application and the contents of which are fully incorporated herein by reference, discloses a method of inversion processing signals from an induction well logging instrument. The instrument includes a transmitter and a plurality of receivers at axially spaced apart locations. The method includes shin effect correcting the responses of the receivers by extrapolating the receiver responses to zero frequency. A model is generated of the media surrounding said instnunent.
Conductivities of elements in the model are adjusted so that a measure of misfit l0 between the shin-effect corrected receiver responses and simulated receiver responses based on the model is minimized. The geometry of the model is adjusted so that the measure of misfit between the shin-effect corrected receiver responses and the simulated receiver responses based on the model is ftirther minimized.
15 The method of the '773 patent and other prior art processing methods for well logs are developed for vertical or near vertical wells. If the well deviation exceeds 30°
- 45° , traditional methods may fail. To interpret the data, 3-D
modeling and inversion must be carried out. However, 3-D modeling and inversion is rarely practical due to the computational time requirements and the unavailability of powerful computers.
20 As would be known to those versed in the art, some information about the subsurface is obtained using measurement-while-drilling ("MWD") methods in horizontal or near horizontal boreholes. In an MWD environment, there is relatively little invasion of the formation by the drilling mud and the relatively short transmitter to receiver spacings may be able to provide sufficient information about the formation resistivity to be useful in applications such as geosteering.
It is desirable to have an invention that is able to give information further away from the borehole in a near horizontal well. Fig. 1 is an end-on view of a borehole within a reservoir rock 20. Surrounding the well is an invaded zone 12. Also shown in Fig. 1 are interfaces 22 and 26 marking a resistivity contrast between the to zone 20 and zones 24 and 28. Typically, the interfaces 22 and 26 could be bed boundaries marking the separation between a reservoir formation 20 and non-reservoir formations 24 and 28. However, the boundary 26 could also be a fluid contact within a formation marling the separation between a water wet region 28 and a hydrocarbon wet region 20 wlule formation 22 is the caprock. Another possible configuration is that 20 is a water wet reservoir rock while 24 is the hydrocarbon saturated region with 28 being the seal beneath the reservoir.
As part of the reservoir evaluation process, it would be desirable to map the boundaries such as 22 and 26 accurately. Such information is useful in the estimation of reserves and in development of completion procedures for developing the reservoir.
The present invention satisfies this need for situations such as in Fig. 1 and also in other situations where information relating to formation properties away from the borehole is needed.
SUMMARY OF THE INVENTION
The present invention is a method for making induction log measurements in an inclined borehole and determining properties of subsurface formations formation away from the borehole. Measurements are made using an induction logging instr~unent having a plurality of transmitter-receiver (T-R) distances. After correction of the data for slcin effects and optionally correcting for eddy currents within the borehole, the shallow measurements (those from short T-R spacing or from high to frequency data) are inverted to give a model of the near borehole (invaded zone resistivity and diameter) and the resistivity of the formation outside the invaded zone.
Using this model, a prediction is made of the data measured by the mid-level and deep sensors (long T-R spacings). A discrepancy between these predicted values and the actual measurements made by the midlevel and deep sensors is indicative of 15 additional layer boundaries in the proximity of the wellbore. The midlevel and deep sensor data are corrected for invasion effects and inverted using the model derived from the shallow data to give an estimate of the distance to these additional layer boundaries and of the resistivity beyond these additional boundaries. In this manner, resistivity measurements in a near horizontal borehole may be used to map subsurface 20 layers away from the borehole. In one embodiment of the invention, the quadrature component of the signal induced in the receiver coils is used. This quadrature component is independent of T-R spacing and has a larger depth of penetration than the commonly used in-phase component of the receiver signal.
BRIEF DESCRIPTION OF THE DRAWINGS
Fig. l is a transverse cross section of a logging tool within a near horizontal borehole showing subsurface formations in proximity to the borehole.
Fig. 2 is a schematic illustration of a logging tool within a near horizontal borehole to showing the transmitter and receiver configuration.
Fig. 3 is a flow chart of the major steps of the method of the present invention for obtaining information about subsurface formations in the vicinity of a borehole.
15 Fig. 4 is an illustration of the geometry of an inclined borehole within the subsurface with uncomformable beds.
DESCRIPTION OF THE PREFERRED EMBODIMENT
2o The present invention is a method for determination of boundaries and material properties of regions having different electrical resisitivities away from a near horizontal borehole. The invention is best understood by reference to Figs. 1-4 and the related discussion hereunder.
FIG. 2 shows an induction well logging instrument 100 disposed in a generally horizontal wellbore 110 penetrating earth formations. The instrument 100 is typically lowered into the wellbore 100 at one end of an armored electrical cable 102 by suitable apparatus at the surface of the earth (not shown). An instrument which will generate appropriate induction signals for performing the process of the present invention is to described, for example, in U.S. Pat. No. 5,452,761 issued to Beard et al.
and having the same assignee as the present application. The instrmnent described in the Beard et al '761 patent is not an exclusive representation of well logging instruments which can generate signals for performing the process of the present invention and therefore is not meant as a limitation on the present invention.
The instrument 100 can include a telemetry/signal processing unit 104 (SPU).
The SPU 104 can include a source of alternating current (not shown separately) to be conducted through a transmitter coil 106 disposed on the instrument 100 Receiver coils 108A-108F can be disposed at axially spaced apart locations along the instrument 100 The SPU 104 include receiver circuits (not shown separately) connected to the receiver coils 108A-108F for detecting voltages induced in each of the receiver coils 108A-108F The SPU 104 can also impart signals to the cable corresponding to the magnitude of the voltages induced in each of the receiver coils 108A-108F. It is to be understood that the number of receiver coils and the arrangement of the coils shown in the instrument in FIG.1 is not meant to be a limitation on the present invention. It is to be further understood that the receiver coils shown in FIG.1 can either be single coils or a type of coil known in the art as "mutually balanced" wherein a second coil(not shown) is series connected in inverse polarity to the primary coil (not shown) so as to reduce the effect of direct iilduction from the transmitter coil 106.
As is understood by those slcilled in the art, the alternating current passing through the transmitter coil 106 induces eddy currents in the earth formations 120, 124 and 128. The eddy currents correspond in magnitude to the electrical conductivity of the earth formations 120,124 and 128 and to the relative position of the particular earth formation with respect to the transmitter coil 106. The eddy currents in turn induce voltages in the receiver coils l OSA-108F, the magnitude of which depends on both the eddy current magnitude and the relative position of the earth formation with respect to the individual receiver coil 108A-108F
A more closely spaced receiver coil such as IOSA would have mare of its voltage induced by eddy currents flowing from entirely within a closer formation such as 120 than would be the case for a longer spaced receiver coil such as 108F.
Conversely, the eddy currents which induce the voltages in receiver coil 108A
would more likely correspond to the conductivity within a zone such as shown at 112, which is affected by fluid "invasion" into its pore spaces from the liquid phase of a fluid used to drill the wellbore (commonly known as "drilling mud", the liquid phase known as "mud filtrate"). .
The drilling mud itself can be electrically conductive. Eddy currents can flow 1 o in such conductive drilling mud, and therefore the voltages induced in each of the receiver coils 108A-108F can also partially depend on the conductivity of the mud and the diameter of the wellbore. As is understood by those spilled in the art, the wellbore diameter is subject to variation as a result of "caving" or "washout".
15 The signals corresponding to the voltages in each receiver coil 108A-108F
can be transmitted along the cable 102 to surface electronics (not shown). The surface electronics can include processors (not shown) for analyzing the signals from the instrument 100 according to the present invention.
2o As is understood by those skilled in the art, the correspondence between the magnitude of the voltages induced in each receiver coil 108A-108F and the conductivity of the media surrounding the instrument 100 is affected by a phenomenon referred to as the "skin effect". In the present embodiment of the invention, the voltage signals induced in each receiver coil 108A-108F can be adjusted to compensate for the skin effect to more precisely correspond to the conductivity of the media surrounding the instrument 100. The present invention uses the method described in the '773 patent to compensate for the skin effect and produce a calculated response for the zero frequency limit of the apparent conductivity at each of the receiver coils 108A-108F.
l0 After the skin effect correction is performed, it is desirable to further adjust the signals from each one of the receiver coils 108A-108F for the effect of the eddy current signal which originates in the wellbore 110. Methods for correcting the receiver coils signals are known in the art, one of which is described, for example in G. P. Grove and G. N. Minerbo, "An Adaptive Borehole Correction Scheme for Array 15 Induction Tools", SPWLA 32nd Annual Logging Symposium, pp 1 25, 1991). It is to be understood that the step of correcting the receiver coil signals for the effects of the conductivity of the wellbore 110 is optional in performing the process of the present invention. Correction for wellbore 110 conductivity can improve the accuracy of the final result. Under some conditions, however, particularly when the drilling mud 20 comprises a substantially non conductive liquid phase, the diameter of the wellbore D is relatively small (typically in a range of 5 to 8 inches where the diameter of the instrument 100 is about l0 3 5/8 inches), and the instrument 100 is substantially centralized in the wellbore I10, the overall effect of the wellbore 110 on the receiver coil signals can be negligible, and the step of correction can be omitted without seriously affecting the accuracy of the result of the inversion process.
The '773 invention describes a method in which data from an instrument such as the one described above are used to determine a conductivity model for the to formation surrounding a vertical borehole. As would be known to those versed in the art, in a vertical borehole that intersects a number of substantially horizontal layers, the conductivity model of the earth formations is a lot more complicated than in the present situation. For a vertical borehole, each layer has a different conductivity associated with it and the response of each receiver coil depends upon the relative 15 position of the receiver coil to the individual layer boundaries.
The use of the method of the '773 invention here is best understood with reference to Fig. 3. The process starts at 200 with the data from all the receivers 108A -108F. The shin effect correction is applied 202 and optionally, a correction for 2o fluid in the borehole is applied 204. With the conf guration of the tool in a near horizontal borehole as shown in Figs. 1 and 2, those versed in the art would recognize that the near receivers, e.g., 108A -108B are responsive almost entirely to the invaded zone 12, 112, while the intermediate receivers 108D -108D are also affected by the material inunediately outside the invaded zone 20,120, referred to hereafter as the uninvaded zone. Finally, the far sensors 108E -108F are also affected by data from the shoulder zones described below. Using the method of the '773 patent, these data from the near and intermediate receivers are inverted to give parameters of the invaded zone and the uninvaded zone. These include the conductivity of the invaded zone and the minvaded zone and the thicl~ness of the invaded zone. Those versed in the art would recognize that what is important is the depth of penetration into the to formation of the induced electromagnetic signal. In the embodiment of the '773 patent, the near receivers are "shallow" sensors, the intermediate receivers are "midlevel" sensors and the far receivers are the "deep" sensors. Other embodiments exist in which shallow data is obtained by high frequency signals while the deep data is obtained by low frequency electromagnetic induction signals. Such variations are intended to be within the scope of the present application. Accordingly, the near receivers such as 108A -108B will hereafter be referred to as the shallow sensors, 108C -108D as sensors and the far receivers such as 108E -108Fwi11 be referred to as the deep sensors.
2o Using the conductivity derived for the near zone, a predicted values for corrected measurements in the midlevel and deep sensors axe obtained 208 assuming that there are no boundaries such as 122 and 126, i.e., that the formation is radially substantially homogenous with respect to the borehole axis. These predicted values at the midlevel and deep sensors are compared to the actual measurements 210. At 212, a comparison between the predicted and actual corrected measurements is made.
If this difference is less than a set threshold value, then the process stops 2I4: the midlevel and deep data are consistent with the model derived from the shallow data and there are no boundaries such as 122 and 126 present. If the difference between the predicted and actual corrected measurements exceeds the threshold 212 then the midlevel and deep data are not consistent with the shallow data. At this point, the Io data are corrected for invasion 215 and the inversion method of the '773 patent is used to derive an enhanced model 216 including therein a distance to a boundary such as 122,126 and a conductivity on the far side 124,128 (i.e., away from the borehole) of the boundary.
IS Those versed in the art would recognize that with azimuthally symmetrical tools such as used in the '773 patent, the inversion would not provide any information regarding the direction of the additional boundary. The additional boundary in the enhanced model is called a shoulder boundary and the layer on the side of the shoulder boundary away from the borehole is called a shoulder layer.
Using the enhanced model, corrected data for the midlevel and deep sensors is predicted and compared with the actual corrected midlevel and deep sensor data 218.
If this difference 220 is less than a set threshold value, then the process stops 214: the deep data are consistent with the enhanced model having a single shoulder layer. If the difference between the predicted and actual corrected data exceeds the threshold then the deep data are not consistent with the shallow data and a model having a single shoulder layer. At this point, the method of the '773 patent is used to derive an refined model 224 including therein a distance to a second shoulder and a conductivity for the second shoulder layer.
l0 In a preferred embodiment of the invention, the two shoulder layers are on opposite sides of the borehole and are parallel to each other, i.e., the inversion and modeling is one dimensional. As noted above, the enhanced model and the refined model by themselves camiot define the direction in which the one or two shoulders lie. However, based upon the relative values of the resistivity of the uninvaded zone 15 and the one or two shoulder layers, an inference can be drawn as to the nature and position of the shoulder layers.
As an example only, and without limitation, if the uninvaded zone has a low conductivity, consistent with a hydrocarbon rich reservoir rock and the first shoulder 20 layer has a very high conductivity, it suggests that the first shoulder layer is a water layer. This water layer would be positioned below the borehole. The second shoulder layer would then be a caprock above the borehole. Conversely, if the uninvaded zone has a high conductivity, indicative of a water saturated reservoir, then a shoulder zone with low conductivity could be indicative of a hydrocarbon/water boundary positioned above the borehole. Those versed in the art would recognize various combinations of the resistivities of the invaded zone, the first and second shoulder zones, as being diagnostic of various combinations of water saturated rocl, gas saturated rock oil saturated roclc and impermeable caproclc.
Fig. 4 illustrates a situation in which a near horizontal wellbore 410 intersects l0 a plurality of steeply dipping formations 412, 414, 416, 418. W version of induction log data acquired in such a wellbore males it possible to determine the subsurface structure in terms of layer tluclnesses. Data from the shallow sensors are inverted in a 2-D inversion to give properties of the formation in the irninediate vicinity of the borehole and, in particular, layer tluclnesses l1 , l2 . . . of the layers as traversed by 15 the borehole. These Iayer resistivities are then used as an initial estimate for a 1-D
inversion of the data from the deep sensors to give the shoulder boundaries 420, 422 and the shoulder beds 424, 426. The methodology methodology discussed above may be used to obtain a resistivity model of the subsurface.
20 Another aspect of the invention uses multifrequency and/or quadrature data as input to the inversion process. United States Patent 5,884,227 having the same inventors and assignee as the present application and the contents of which are fully incorporated herein by reference, shows that when data are acquired at m angular frequencies tvl , cv2 ,w3 ,. . . cvm , the in-phase components of voltages in the receiver coils 108A -108F in Fig. 2 generally correspond to the imaginary component of the magnetic field. This imaginary component of the magnetic field may be represented by a Taylor series of the form k=m ,5 H ) - ~ 0 k12~ k/2 1 k=2 A general representation of the apparent conductivity from the imaginary component to of the magnetic field is given by a Taylor series expansion of the form k=w ~a - ~ skl2 ~ k/2 ; SZ J-1 = O; j = h~~ ... . , ~
k=0 where s ~,2 are coefficients depending on the conductivity distribution and the tool's geometric configuration. They do not depend upon the frequency.
Considering the series for each of the measured frequencies, this gives CO 1 ~ 1/2 ~ 3/2 , . , s a( ~ n/2 1) 1 1 1 1/2 3/2 n/2 0'a(C02) 1 CV2 CUZ . . . (U2 S1/2 1/2 3/2 n/2 6a(~m-1) 1 ~m-1 ~m-1 ~ ~ ~ ~m-1 1/2 3/2 n/2 6a 1 CDm (.Um . . CUm ~n12 (lD~n ) (3) where m represents the number of component frequencies, and h can be selected from a Taylor series expansion for an apparent conductivity integral.
The ultimate objective of the expansion series in Eq. (3) is to solve for the coefficient at the zero frequency limit- sp. Solving a system according to Eq. (3) will calculate a receiver coil response approximately equal to that described in H. G. Doll, "Introduction to Induction Logging and Application to Logging Wells Drilled with Oil Base Mud", AIME
Transactions Vol. 1, No. 6, pp 148-162 (1949), for the zero frequency limit of the apparent conductivity at each one of the receiver coils (108A 108F in FIG.1).
The leading terms of the expansion given by Eq. (1) are of the form h(.fl) _ .flhl + f13/2h2 + flz y~+... (4) where f1 is the frequency. The component fl hl is substantially dependent on the conductivity in the wellbore region, i.e., this term depends mainly on the conductivity of the wellbore fluid. PCT Patent Publication WO 98!0073 shows that by making the magnetic field measurement at a second frequency f~ , a second field measurement h2 is obtained h(.fz ) _ .fzhi + .fi ~Zhz + .fi hs+... I (5) A linear combination of the measurements at the two frequencies h(.fi ~ .fz ) = h(.fi ) - (.fl ~ .fz ~h(,fz ) ( makes it possible to eliminate terms that are linear and frequency. As would be known to those versed in the art, the term that is linear in frequency is primarily responsive to the borehole. Having done this, it is possible to invert the multiple frequency data to obtain the bed boundaries for the configurations of Figs. l, 4 or 5.
Specifically, for the configuration of Fig. l, a two-dimensional (2-D) inversion of the shallow data is carried out to give the formation resistivity for the region 20 of Fig. 1 as at 206 in Fig. 3. The rest of the analysis proceeds as discussed above with reference to Fig. 3.
For the subsurface configuration of Fig. 4, single-frequency shallow data may be used to derive the borehole and formation resistivities as discussed above with respect to Fig. 1. Alternatively, multiple frequency data may be used for a one-dimensional inversion of the deep data to obtain the positions of the boundaries 420, 422 and the resistivities of the layers 424, 426. In a preferred embodiment of the invention, the multiple frequency data are acquired at two frequencies.
One embodiment of the present invention includes measurements of the quadrature component of the induced signals in the receiver coils (the quadrature to component is 90° out of phase with the phase of the transmitter signal). The quadrature component of eq. (2) gives the following equation:
O- CO 1 (~ 3/2 ~ 2 . . , a( 1) ~ n/2 -1 3/2 2 n/2 6a (~2 1 ~2 ~2 ~ ~ ~ ~2 'S3/2 ) ,S2 3/2 2 n/2 6a(~m-i), 1 ant-1 ~m-1 ' ' ' ~m-1 3/2 2 n/2 6a (CUm 1 CUm CUZ . . [pm snl2 ) In eq. (7), the 6a 's are the apparent conductivities measured at the corresponding frequencies. Solving the system of eq. (7) with respect to s3/2 gives the zero-frequency limit for quadrature measurements. It is known in literature that this component is independent of the tool spacing and it also has a large depth of investigation. These two properties are important for induction logging in near horizontal boreholes and inversion of the data makes it possible to identify bed boundaries a distance from the borehole.
Those versed in the art would recognize that if the number of frequencies m at which quadrature measurements are made is greater than h, the number of terms in the series expansion given by eq. (2), then eq. (~ is an overdetennined system of equations. Inversion of the values of s3/~ mares it possible to obtain the subsurface to model corresponding to any of Figs. l, 4 and 5. In a preferred embodiment of the invention, the in-phase components of the shallow data are used to obtain the resistivity of the formation near the borehole and the quadrature component data at a low frequency limit are used to obtain the distances and resistivities of the shoulder beds.
Those slcilled in the art will devise other embodiments of this invention which do not depart from the spirit of the invention as disclosed herein.
Accordingly, the invention should be limited in scope only by the attached claims.
The invention is related to the field of induction logging of subsurface earth formations. Specifically, it relates to the use of induction logging for obtaining structural information of the subsurface away from a borehole.
BACKGROUND OF THE INVENTION
to Induction resistivity well logging devices are l~nown in the art for measuring the electrical resistivity of earth formations penetrated by wellbores. Electrical resistivity measurements are used for, among other things, to infer the fluid content of pore spaces of the earth formations.
15 Induction electrical resistivity instruments l~nown in the art, such as one described in Hunka, J. "A New Resistivity Measurement System for Deep Formation Imaging and High Resolution Formation Evaluation", paper no. 20559, Society of Petroleum Engineers, Richardson, Tex., 1990, typically measure voltages induced in receiver coils positioned at spaced apart locations along the instrument. The voltages are induced by magnetic fields 20 generated by eddy currents flowing in the earth formations around the instrument, in response to alternating currents of various frequencies being passed through a transmitter coil. The voltages induced in a particular one of the receiver coils are dependent on the electrical conductivity (which is inversely related to resistivity) of the media, including the earth formations, surrounding the instrument, and are dependent on the spacing of the particular receiver coil with respect to the transmitter coil. Generally, the greater the spacing between the transmitter coil and the particular receiver coil, the greater is the radial depth from the wellbore into the earth formation from which the measurement at that receiver coil corresponds. The vertical resolution of the measurements made by the closer-spaced receiver coils, however, is proportionately finer than is the vertical resolution of the measurements made by the to more distant receiver coils.
When the wellbore is drilled, fluids from the drilling process may be forced into the pore spaces of some of the earth formations, changing their fluid content and therefore their resistivity. The process of fluid being forced into the pore spaces is 15 generally referred to as "invasion". Wellbore can be drilled with fluids having differing electrical resistivities, and the wellbore itself can have widely varying diameters over its length. The effects of variable wellbore diameter and differing wellbore fluid resistivity can indeterminately affect the total voltage being induced in each one of the receiver coils, because eddy currents also can flow within the 2o conductive wellbore fluid and within that portion of the earth formations surrounding the wellbore subject to the invasion (the so-called "invaded zone") United States Patent 5,703,773 ( referred hereafter as the " '773 patent" ) having the same assignee as the present application and the contents of which are fully incorporated herein by reference, discloses a method of inversion processing signals from an induction well logging instrument. The instrument includes a transmitter and a plurality of receivers at axially spaced apart locations. The method includes shin effect correcting the responses of the receivers by extrapolating the receiver responses to zero frequency. A model is generated of the media surrounding said instnunent.
Conductivities of elements in the model are adjusted so that a measure of misfit l0 between the shin-effect corrected receiver responses and simulated receiver responses based on the model is minimized. The geometry of the model is adjusted so that the measure of misfit between the shin-effect corrected receiver responses and the simulated receiver responses based on the model is ftirther minimized.
15 The method of the '773 patent and other prior art processing methods for well logs are developed for vertical or near vertical wells. If the well deviation exceeds 30°
- 45° , traditional methods may fail. To interpret the data, 3-D
modeling and inversion must be carried out. However, 3-D modeling and inversion is rarely practical due to the computational time requirements and the unavailability of powerful computers.
20 As would be known to those versed in the art, some information about the subsurface is obtained using measurement-while-drilling ("MWD") methods in horizontal or near horizontal boreholes. In an MWD environment, there is relatively little invasion of the formation by the drilling mud and the relatively short transmitter to receiver spacings may be able to provide sufficient information about the formation resistivity to be useful in applications such as geosteering.
It is desirable to have an invention that is able to give information further away from the borehole in a near horizontal well. Fig. 1 is an end-on view of a borehole within a reservoir rock 20. Surrounding the well is an invaded zone 12. Also shown in Fig. 1 are interfaces 22 and 26 marking a resistivity contrast between the to zone 20 and zones 24 and 28. Typically, the interfaces 22 and 26 could be bed boundaries marking the separation between a reservoir formation 20 and non-reservoir formations 24 and 28. However, the boundary 26 could also be a fluid contact within a formation marling the separation between a water wet region 28 and a hydrocarbon wet region 20 wlule formation 22 is the caprock. Another possible configuration is that 20 is a water wet reservoir rock while 24 is the hydrocarbon saturated region with 28 being the seal beneath the reservoir.
As part of the reservoir evaluation process, it would be desirable to map the boundaries such as 22 and 26 accurately. Such information is useful in the estimation of reserves and in development of completion procedures for developing the reservoir.
The present invention satisfies this need for situations such as in Fig. 1 and also in other situations where information relating to formation properties away from the borehole is needed.
SUMMARY OF THE INVENTION
The present invention is a method for making induction log measurements in an inclined borehole and determining properties of subsurface formations formation away from the borehole. Measurements are made using an induction logging instr~unent having a plurality of transmitter-receiver (T-R) distances. After correction of the data for slcin effects and optionally correcting for eddy currents within the borehole, the shallow measurements (those from short T-R spacing or from high to frequency data) are inverted to give a model of the near borehole (invaded zone resistivity and diameter) and the resistivity of the formation outside the invaded zone.
Using this model, a prediction is made of the data measured by the mid-level and deep sensors (long T-R spacings). A discrepancy between these predicted values and the actual measurements made by the midlevel and deep sensors is indicative of 15 additional layer boundaries in the proximity of the wellbore. The midlevel and deep sensor data are corrected for invasion effects and inverted using the model derived from the shallow data to give an estimate of the distance to these additional layer boundaries and of the resistivity beyond these additional boundaries. In this manner, resistivity measurements in a near horizontal borehole may be used to map subsurface 20 layers away from the borehole. In one embodiment of the invention, the quadrature component of the signal induced in the receiver coils is used. This quadrature component is independent of T-R spacing and has a larger depth of penetration than the commonly used in-phase component of the receiver signal.
BRIEF DESCRIPTION OF THE DRAWINGS
Fig. l is a transverse cross section of a logging tool within a near horizontal borehole showing subsurface formations in proximity to the borehole.
Fig. 2 is a schematic illustration of a logging tool within a near horizontal borehole to showing the transmitter and receiver configuration.
Fig. 3 is a flow chart of the major steps of the method of the present invention for obtaining information about subsurface formations in the vicinity of a borehole.
15 Fig. 4 is an illustration of the geometry of an inclined borehole within the subsurface with uncomformable beds.
DESCRIPTION OF THE PREFERRED EMBODIMENT
2o The present invention is a method for determination of boundaries and material properties of regions having different electrical resisitivities away from a near horizontal borehole. The invention is best understood by reference to Figs. 1-4 and the related discussion hereunder.
FIG. 2 shows an induction well logging instrument 100 disposed in a generally horizontal wellbore 110 penetrating earth formations. The instrument 100 is typically lowered into the wellbore 100 at one end of an armored electrical cable 102 by suitable apparatus at the surface of the earth (not shown). An instrument which will generate appropriate induction signals for performing the process of the present invention is to described, for example, in U.S. Pat. No. 5,452,761 issued to Beard et al.
and having the same assignee as the present application. The instrmnent described in the Beard et al '761 patent is not an exclusive representation of well logging instruments which can generate signals for performing the process of the present invention and therefore is not meant as a limitation on the present invention.
The instrument 100 can include a telemetry/signal processing unit 104 (SPU).
The SPU 104 can include a source of alternating current (not shown separately) to be conducted through a transmitter coil 106 disposed on the instrument 100 Receiver coils 108A-108F can be disposed at axially spaced apart locations along the instrument 100 The SPU 104 include receiver circuits (not shown separately) connected to the receiver coils 108A-108F for detecting voltages induced in each of the receiver coils 108A-108F The SPU 104 can also impart signals to the cable corresponding to the magnitude of the voltages induced in each of the receiver coils 108A-108F. It is to be understood that the number of receiver coils and the arrangement of the coils shown in the instrument in FIG.1 is not meant to be a limitation on the present invention. It is to be further understood that the receiver coils shown in FIG.1 can either be single coils or a type of coil known in the art as "mutually balanced" wherein a second coil(not shown) is series connected in inverse polarity to the primary coil (not shown) so as to reduce the effect of direct iilduction from the transmitter coil 106.
As is understood by those slcilled in the art, the alternating current passing through the transmitter coil 106 induces eddy currents in the earth formations 120, 124 and 128. The eddy currents correspond in magnitude to the electrical conductivity of the earth formations 120,124 and 128 and to the relative position of the particular earth formation with respect to the transmitter coil 106. The eddy currents in turn induce voltages in the receiver coils l OSA-108F, the magnitude of which depends on both the eddy current magnitude and the relative position of the earth formation with respect to the individual receiver coil 108A-108F
A more closely spaced receiver coil such as IOSA would have mare of its voltage induced by eddy currents flowing from entirely within a closer formation such as 120 than would be the case for a longer spaced receiver coil such as 108F.
Conversely, the eddy currents which induce the voltages in receiver coil 108A
would more likely correspond to the conductivity within a zone such as shown at 112, which is affected by fluid "invasion" into its pore spaces from the liquid phase of a fluid used to drill the wellbore (commonly known as "drilling mud", the liquid phase known as "mud filtrate"). .
The drilling mud itself can be electrically conductive. Eddy currents can flow 1 o in such conductive drilling mud, and therefore the voltages induced in each of the receiver coils 108A-108F can also partially depend on the conductivity of the mud and the diameter of the wellbore. As is understood by those spilled in the art, the wellbore diameter is subject to variation as a result of "caving" or "washout".
15 The signals corresponding to the voltages in each receiver coil 108A-108F
can be transmitted along the cable 102 to surface electronics (not shown). The surface electronics can include processors (not shown) for analyzing the signals from the instrument 100 according to the present invention.
2o As is understood by those skilled in the art, the correspondence between the magnitude of the voltages induced in each receiver coil 108A-108F and the conductivity of the media surrounding the instrument 100 is affected by a phenomenon referred to as the "skin effect". In the present embodiment of the invention, the voltage signals induced in each receiver coil 108A-108F can be adjusted to compensate for the skin effect to more precisely correspond to the conductivity of the media surrounding the instrument 100. The present invention uses the method described in the '773 patent to compensate for the skin effect and produce a calculated response for the zero frequency limit of the apparent conductivity at each of the receiver coils 108A-108F.
l0 After the skin effect correction is performed, it is desirable to further adjust the signals from each one of the receiver coils 108A-108F for the effect of the eddy current signal which originates in the wellbore 110. Methods for correcting the receiver coils signals are known in the art, one of which is described, for example in G. P. Grove and G. N. Minerbo, "An Adaptive Borehole Correction Scheme for Array 15 Induction Tools", SPWLA 32nd Annual Logging Symposium, pp 1 25, 1991). It is to be understood that the step of correcting the receiver coil signals for the effects of the conductivity of the wellbore 110 is optional in performing the process of the present invention. Correction for wellbore 110 conductivity can improve the accuracy of the final result. Under some conditions, however, particularly when the drilling mud 20 comprises a substantially non conductive liquid phase, the diameter of the wellbore D is relatively small (typically in a range of 5 to 8 inches where the diameter of the instrument 100 is about l0 3 5/8 inches), and the instrument 100 is substantially centralized in the wellbore I10, the overall effect of the wellbore 110 on the receiver coil signals can be negligible, and the step of correction can be omitted without seriously affecting the accuracy of the result of the inversion process.
The '773 invention describes a method in which data from an instrument such as the one described above are used to determine a conductivity model for the to formation surrounding a vertical borehole. As would be known to those versed in the art, in a vertical borehole that intersects a number of substantially horizontal layers, the conductivity model of the earth formations is a lot more complicated than in the present situation. For a vertical borehole, each layer has a different conductivity associated with it and the response of each receiver coil depends upon the relative 15 position of the receiver coil to the individual layer boundaries.
The use of the method of the '773 invention here is best understood with reference to Fig. 3. The process starts at 200 with the data from all the receivers 108A -108F. The shin effect correction is applied 202 and optionally, a correction for 2o fluid in the borehole is applied 204. With the conf guration of the tool in a near horizontal borehole as shown in Figs. 1 and 2, those versed in the art would recognize that the near receivers, e.g., 108A -108B are responsive almost entirely to the invaded zone 12, 112, while the intermediate receivers 108D -108D are also affected by the material inunediately outside the invaded zone 20,120, referred to hereafter as the uninvaded zone. Finally, the far sensors 108E -108F are also affected by data from the shoulder zones described below. Using the method of the '773 patent, these data from the near and intermediate receivers are inverted to give parameters of the invaded zone and the uninvaded zone. These include the conductivity of the invaded zone and the minvaded zone and the thicl~ness of the invaded zone. Those versed in the art would recognize that what is important is the depth of penetration into the to formation of the induced electromagnetic signal. In the embodiment of the '773 patent, the near receivers are "shallow" sensors, the intermediate receivers are "midlevel" sensors and the far receivers are the "deep" sensors. Other embodiments exist in which shallow data is obtained by high frequency signals while the deep data is obtained by low frequency electromagnetic induction signals. Such variations are intended to be within the scope of the present application. Accordingly, the near receivers such as 108A -108B will hereafter be referred to as the shallow sensors, 108C -108D as sensors and the far receivers such as 108E -108Fwi11 be referred to as the deep sensors.
2o Using the conductivity derived for the near zone, a predicted values for corrected measurements in the midlevel and deep sensors axe obtained 208 assuming that there are no boundaries such as 122 and 126, i.e., that the formation is radially substantially homogenous with respect to the borehole axis. These predicted values at the midlevel and deep sensors are compared to the actual measurements 210. At 212, a comparison between the predicted and actual corrected measurements is made.
If this difference is less than a set threshold value, then the process stops 2I4: the midlevel and deep data are consistent with the model derived from the shallow data and there are no boundaries such as 122 and 126 present. If the difference between the predicted and actual corrected measurements exceeds the threshold 212 then the midlevel and deep data are not consistent with the shallow data. At this point, the Io data are corrected for invasion 215 and the inversion method of the '773 patent is used to derive an enhanced model 216 including therein a distance to a boundary such as 122,126 and a conductivity on the far side 124,128 (i.e., away from the borehole) of the boundary.
IS Those versed in the art would recognize that with azimuthally symmetrical tools such as used in the '773 patent, the inversion would not provide any information regarding the direction of the additional boundary. The additional boundary in the enhanced model is called a shoulder boundary and the layer on the side of the shoulder boundary away from the borehole is called a shoulder layer.
Using the enhanced model, corrected data for the midlevel and deep sensors is predicted and compared with the actual corrected midlevel and deep sensor data 218.
If this difference 220 is less than a set threshold value, then the process stops 214: the deep data are consistent with the enhanced model having a single shoulder layer. If the difference between the predicted and actual corrected data exceeds the threshold then the deep data are not consistent with the shallow data and a model having a single shoulder layer. At this point, the method of the '773 patent is used to derive an refined model 224 including therein a distance to a second shoulder and a conductivity for the second shoulder layer.
l0 In a preferred embodiment of the invention, the two shoulder layers are on opposite sides of the borehole and are parallel to each other, i.e., the inversion and modeling is one dimensional. As noted above, the enhanced model and the refined model by themselves camiot define the direction in which the one or two shoulders lie. However, based upon the relative values of the resistivity of the uninvaded zone 15 and the one or two shoulder layers, an inference can be drawn as to the nature and position of the shoulder layers.
As an example only, and without limitation, if the uninvaded zone has a low conductivity, consistent with a hydrocarbon rich reservoir rock and the first shoulder 20 layer has a very high conductivity, it suggests that the first shoulder layer is a water layer. This water layer would be positioned below the borehole. The second shoulder layer would then be a caprock above the borehole. Conversely, if the uninvaded zone has a high conductivity, indicative of a water saturated reservoir, then a shoulder zone with low conductivity could be indicative of a hydrocarbon/water boundary positioned above the borehole. Those versed in the art would recognize various combinations of the resistivities of the invaded zone, the first and second shoulder zones, as being diagnostic of various combinations of water saturated rocl, gas saturated rock oil saturated roclc and impermeable caproclc.
Fig. 4 illustrates a situation in which a near horizontal wellbore 410 intersects l0 a plurality of steeply dipping formations 412, 414, 416, 418. W version of induction log data acquired in such a wellbore males it possible to determine the subsurface structure in terms of layer tluclnesses. Data from the shallow sensors are inverted in a 2-D inversion to give properties of the formation in the irninediate vicinity of the borehole and, in particular, layer tluclnesses l1 , l2 . . . of the layers as traversed by 15 the borehole. These Iayer resistivities are then used as an initial estimate for a 1-D
inversion of the data from the deep sensors to give the shoulder boundaries 420, 422 and the shoulder beds 424, 426. The methodology methodology discussed above may be used to obtain a resistivity model of the subsurface.
20 Another aspect of the invention uses multifrequency and/or quadrature data as input to the inversion process. United States Patent 5,884,227 having the same inventors and assignee as the present application and the contents of which are fully incorporated herein by reference, shows that when data are acquired at m angular frequencies tvl , cv2 ,w3 ,. . . cvm , the in-phase components of voltages in the receiver coils 108A -108F in Fig. 2 generally correspond to the imaginary component of the magnetic field. This imaginary component of the magnetic field may be represented by a Taylor series of the form k=m ,5 H ) - ~ 0 k12~ k/2 1 k=2 A general representation of the apparent conductivity from the imaginary component to of the magnetic field is given by a Taylor series expansion of the form k=w ~a - ~ skl2 ~ k/2 ; SZ J-1 = O; j = h~~ ... . , ~
k=0 where s ~,2 are coefficients depending on the conductivity distribution and the tool's geometric configuration. They do not depend upon the frequency.
Considering the series for each of the measured frequencies, this gives CO 1 ~ 1/2 ~ 3/2 , . , s a( ~ n/2 1) 1 1 1 1/2 3/2 n/2 0'a(C02) 1 CV2 CUZ . . . (U2 S1/2 1/2 3/2 n/2 6a(~m-1) 1 ~m-1 ~m-1 ~ ~ ~ ~m-1 1/2 3/2 n/2 6a 1 CDm (.Um . . CUm ~n12 (lD~n ) (3) where m represents the number of component frequencies, and h can be selected from a Taylor series expansion for an apparent conductivity integral.
The ultimate objective of the expansion series in Eq. (3) is to solve for the coefficient at the zero frequency limit- sp. Solving a system according to Eq. (3) will calculate a receiver coil response approximately equal to that described in H. G. Doll, "Introduction to Induction Logging and Application to Logging Wells Drilled with Oil Base Mud", AIME
Transactions Vol. 1, No. 6, pp 148-162 (1949), for the zero frequency limit of the apparent conductivity at each one of the receiver coils (108A 108F in FIG.1).
The leading terms of the expansion given by Eq. (1) are of the form h(.fl) _ .flhl + f13/2h2 + flz y~+... (4) where f1 is the frequency. The component fl hl is substantially dependent on the conductivity in the wellbore region, i.e., this term depends mainly on the conductivity of the wellbore fluid. PCT Patent Publication WO 98!0073 shows that by making the magnetic field measurement at a second frequency f~ , a second field measurement h2 is obtained h(.fz ) _ .fzhi + .fi ~Zhz + .fi hs+... I (5) A linear combination of the measurements at the two frequencies h(.fi ~ .fz ) = h(.fi ) - (.fl ~ .fz ~h(,fz ) ( makes it possible to eliminate terms that are linear and frequency. As would be known to those versed in the art, the term that is linear in frequency is primarily responsive to the borehole. Having done this, it is possible to invert the multiple frequency data to obtain the bed boundaries for the configurations of Figs. l, 4 or 5.
Specifically, for the configuration of Fig. l, a two-dimensional (2-D) inversion of the shallow data is carried out to give the formation resistivity for the region 20 of Fig. 1 as at 206 in Fig. 3. The rest of the analysis proceeds as discussed above with reference to Fig. 3.
For the subsurface configuration of Fig. 4, single-frequency shallow data may be used to derive the borehole and formation resistivities as discussed above with respect to Fig. 1. Alternatively, multiple frequency data may be used for a one-dimensional inversion of the deep data to obtain the positions of the boundaries 420, 422 and the resistivities of the layers 424, 426. In a preferred embodiment of the invention, the multiple frequency data are acquired at two frequencies.
One embodiment of the present invention includes measurements of the quadrature component of the induced signals in the receiver coils (the quadrature to component is 90° out of phase with the phase of the transmitter signal). The quadrature component of eq. (2) gives the following equation:
O- CO 1 (~ 3/2 ~ 2 . . , a( 1) ~ n/2 -1 3/2 2 n/2 6a (~2 1 ~2 ~2 ~ ~ ~ ~2 'S3/2 ) ,S2 3/2 2 n/2 6a(~m-i), 1 ant-1 ~m-1 ' ' ' ~m-1 3/2 2 n/2 6a (CUm 1 CUm CUZ . . [pm snl2 ) In eq. (7), the 6a 's are the apparent conductivities measured at the corresponding frequencies. Solving the system of eq. (7) with respect to s3/2 gives the zero-frequency limit for quadrature measurements. It is known in literature that this component is independent of the tool spacing and it also has a large depth of investigation. These two properties are important for induction logging in near horizontal boreholes and inversion of the data makes it possible to identify bed boundaries a distance from the borehole.
Those versed in the art would recognize that if the number of frequencies m at which quadrature measurements are made is greater than h, the number of terms in the series expansion given by eq. (2), then eq. (~ is an overdetennined system of equations. Inversion of the values of s3/~ mares it possible to obtain the subsurface to model corresponding to any of Figs. l, 4 and 5. In a preferred embodiment of the invention, the in-phase components of the shallow data are used to obtain the resistivity of the formation near the borehole and the quadrature component data at a low frequency limit are used to obtain the distances and resistivities of the shoulder beds.
Those slcilled in the art will devise other embodiments of this invention which do not depart from the spirit of the invention as disclosed herein.
Accordingly, the invention should be limited in scope only by the attached claims.
Claims (33)
1. A method of logging earth formations penetrated by a borehole using a well logging instrument conveyed therein, the method comprising:
(a) activating a transmitter on the instrument to transmit an electromagnetic (EM) wave thereby inducing currents in the earth formations, (b) measuring signals induced in a plurality of receivers on the instrument in response to the currents induced in said earth formations;
(c) processing with a processor a first subset of the measured signals for determining properties of the formation in an invaded zone proximate to the borehole and in an uninvaded zone immediately outside said invaded zone; and (d) using said determined properties and at least one additional subset of the measured signals to determine a distance to at least one boundary outside the invaded zone and properties of the earth formations on a side of the at least one boundary away from the borehole.
(a) activating a transmitter on the instrument to transmit an electromagnetic (EM) wave thereby inducing currents in the earth formations, (b) measuring signals induced in a plurality of receivers on the instrument in response to the currents induced in said earth formations;
(c) processing with a processor a first subset of the measured signals for determining properties of the formation in an invaded zone proximate to the borehole and in an uninvaded zone immediately outside said invaded zone; and (d) using said determined properties and at least one additional subset of the measured signals to determine a distance to at least one boundary outside the invaded zone and properties of the earth formations on a side of the at least one boundary away from the borehole.
2. The method of claim 1 wherein said properties of the invaded zone comprise a resistivity and radius of the invaded zone.
3. The method of claim 2 wherein said properties of said uninvaded zone comprise a resistivity thereof.
4. The method of claim 1 wherein said measured signals comprise shallow, midlevel and deep signals and said first subset comprises shallow and midlevel signals.
5. The method of claim 1 wherein determining properties of the invaded zone and of the uninvaded zone immediately outside the invaded zone further comprises performing an inversion of said first subset of signals.
6. The method of claim 1 wherein the at least one additional subset comprises deep signals.
7. The method of claim 4 wherein said shallow signals comprise signals corresponding to shoat transmitter-receiver spacings, said deep signals comprise corresponding to long transmitter-receiver spacings, and said midlevel signals comprise signals corresponding to intermediate transmitter-receiver spacings.
8. The method of claim 4 wherein said shallow signals comprise signals having a high frequency, said deep signals comprise signals having a low frequency, and said midlevel signals comprise signals having an intermediate frequency.
9. The method of claim 1 wherein the at least one boundary comprises two boundaries.
10. The method of claim 9 wherein the two boundaries are on opposite sides of the borehole.
11. The method of claim 10 further comprising using values of resistivities of the uninvaded zone and of the earth formations on sides of the boundaries away from the borehole for estimating petrophysical and fluid properties of said portions of earth formations.
12. The method of claim 1 wherein determining properties of the formation in the invaded zone further comprises determining and applying a skin effect correction
13. The method of claim 1, wherein the first subset of the measured signals have a phase substantially the same as a phase of said transmitted EM wave.
14. The method of claim 1 wherein the second subset of the measured signals have a phase that is one of (i) substantially the same as the phase of said transmitted EM wave, and, (ii) substantially orthogonal to the phase of said transmitted EM wave.
15. The method of claim 1 further comprising transmitting said electromagnetic wave at a plurality of frequencies.
16. The method of claim 5 further comprising transmitting said electromagnetic wave at a plurality of frequencies and wherein performing said inversion further comprises using a Taylor series expansion.
17. A method of determining a distance to at least one bed boundary in earth formations penetrated by a near horizontal borehole using an induction logging instrument conveyed therein, the method comprising:
(a) using a plurality of transmitter-receiver combinations on the logging instrument having different transmitter-receiver distances to obtain signals indicative of properties of the earth formations;
(b) determining from a first subset of said signals properties of the formation in an invaded zone proximate to the borehole and an uninvaded zone immediately outside said invaded zone, said at least one bed boundary being outside the invaded zone; and (c) using said determined properties and at least one additional subset of said signals for determining a distance to the at least one boundary.
(a) using a plurality of transmitter-receiver combinations on the logging instrument having different transmitter-receiver distances to obtain signals indicative of properties of the earth formations;
(b) determining from a first subset of said signals properties of the formation in an invaded zone proximate to the borehole and an uninvaded zone immediately outside said invaded zone, said at least one bed boundary being outside the invaded zone; and (c) using said determined properties and at least one additional subset of said signals for determining a distance to the at least one boundary.
18. The method of claim 17 wherein said properties of the invaded zone comprise a resistivity and radius of the invaded zone.
19. The method of claim 17 wherein said properties of said uninvaded zone comprise a resistivity thereof.
20. The method of claims 17 wherein said first subset of said signals have smaller transmitter-receiver distances than said second subset of said signals.
21. The method of claim 17 wherein determining properties of the invaded zone and of the uninvaded zone immediately outside the invaded zone further comprises performing an inversion of said first subset of signals.
22. The method of claim 17 wherein the at least one boundary comprises two boundaries disposed on opposite sides of the borehole.
23. The method of claim 17 wherein determining properties of the formation in the invaded zone further comprises determining and applying a skin effect correction
24. The method of claim 1 7 wherein the second subset of the measured signals have a phase that is one of (i) substantially the same as the phase of an electromagnetic wave transmitted by a transmitter, and, (ii) substantially orthogonal to a phase of an electromagnetic wave transmitted by a transmitter.
25. A method of determining a distance to at least one bed boundary in earth formations penetrated by a near horizontal borehole using an induction logging instrument conveyed therein, the method comprising:
(a) using at least ore transmitter and at least one receiver on the logging instrument at a plurality of frequencies to obtain signals indicative of properties of the earth formations;
(b) determining from a first subset of said signals properties of the formation in an invaded zone proximate to the borehole and an uninvaded zone immediately outside said invaded zone, said at least one bed boundary being outside the invaded zone; and (c) using said determined properties and at least one additional subset of said signals for determining a distance to the at least one boundary.
(a) using at least ore transmitter and at least one receiver on the logging instrument at a plurality of frequencies to obtain signals indicative of properties of the earth formations;
(b) determining from a first subset of said signals properties of the formation in an invaded zone proximate to the borehole and an uninvaded zone immediately outside said invaded zone, said at least one bed boundary being outside the invaded zone; and (c) using said determined properties and at least one additional subset of said signals for determining a distance to the at least one boundary.
26. The method of claim 25 wherein said properties of the invaded zone comprise a resistivity and radius of the invaded zone.
27. The method of claim 25 wherein said properties of said uninvaded zone comprise a resistivity thereof.
28. The method of claim 25 wherein said first subset of said signals have smaller frequencies than said second subset of said signals.
29. The method of claim 25 wherein determining properties of the invaded zone and of the uninvaded zone immediately outside the invaded zone further comprises, performing an inversion of said first subset of signals.
30. The method of claims 25 wherein the at least one boundary comprises two boundaries disposed on opposite sides of the borehole.
31. The method of claim 25 wherein determining properties of the formation in the invaded zone further comprises determining and applying a skin effect correction
32. The method of claim 25 wherein the second subset of the measured signals have a phase that is one of (i) substantially the same as the phase of an electromagnetic wave transmitted by a transmitter, and, (ii) substantially orthogonal to a phase of an electromagnetic wave transmitted by a transmitter.
33. The method of claim 25 further comprising representing said signals in a Taylor serves expansion.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US09/519,029 US6308136B1 (en) | 2000-03-03 | 2000-03-03 | Method of interpreting induction logs in horizontal wells |
| US09/519,029 | 2000-03-03 | ||
| PCT/US2001/006877 WO2001067135A1 (en) | 2000-03-03 | 2001-03-05 | Method of interpreting induction logs in horizontal wells |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CA2401898A1 true CA2401898A1 (en) | 2001-09-13 |
Family
ID=24066472
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002401898A Abandoned CA2401898A1 (en) | 2000-03-03 | 2001-03-05 | Method of interpreting induction logs in horizontal wells |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US6308136B1 (en) |
| EP (1) | EP1269222A4 (en) |
| CA (1) | CA2401898A1 (en) |
| NO (1) | NO20024207D0 (en) |
| WO (1) | WO2001067135A1 (en) |
Families Citing this family (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7375530B2 (en) * | 2002-03-04 | 2008-05-20 | Baker Hughes Incorporated | Method for signal enhancement in azimuthal propagation resistivity while drilling |
| US7463035B2 (en) * | 2002-03-04 | 2008-12-09 | Baker Hughes Incorporated | Method and apparatus for the use of multicomponent induction tool for geosteering and formation resistivity data interpretation in horizontal wells |
| US6832159B2 (en) | 2002-07-11 | 2004-12-14 | Schlumberger Technology Corporation | Intelligent diagnosis of environmental influence on well logs with model-based inversion |
| US7414391B2 (en) * | 2002-07-30 | 2008-08-19 | Schlumberger Technology Corporation | Electromagnetic logging tool calibration system |
| US6810331B2 (en) * | 2002-09-25 | 2004-10-26 | Halliburton Energy Services, Inc. | Fixed-depth of investigation log for multi-spacing multi-frequency LWD resistivity tools |
| GB2410336B (en) * | 2002-09-27 | 2006-03-08 | Baker Hughes Inc | A method for resistivity anisotropy determination in conductive borehole environments |
| BR0314800A (en) * | 2002-09-27 | 2005-08-02 | Baker Hughes Inc | Process for determination of resistivity anisotropy in nearby vertical wells |
| US7063174B2 (en) * | 2002-11-12 | 2006-06-20 | Baker Hughes Incorporated | Method for reservoir navigation using formation pressure testing measurement while drilling |
| US7031839B2 (en) * | 2002-11-15 | 2006-04-18 | Baker Hughes Incorporated | Multi-frequency focusing for MWD resistivity tools |
| US6906521B2 (en) * | 2002-11-15 | 2005-06-14 | Baker Hughes Incorporated | Multi-frequency focusing for MWD resistivity tools |
| US7093672B2 (en) * | 2003-02-11 | 2006-08-22 | Schlumberger Technology Corporation | Systems for deep resistivity while drilling for proactive geosteering |
| US6944546B2 (en) * | 2003-10-01 | 2005-09-13 | Halliburton Energy Services, Inc. | Method and apparatus for inversion processing of well logging data in a selected pattern space |
| US8060310B2 (en) * | 2004-06-15 | 2011-11-15 | Baker Hughes Incorporated | Geosteering in earth formations using multicomponent induction measurements |
| US7269515B2 (en) * | 2004-06-15 | 2007-09-11 | Baker Hughes Incorporated | Geosteering in anisotropic formations using multicomponent induction measurements |
| US7274991B2 (en) * | 2004-06-15 | 2007-09-25 | Baker Hughes Incorporated | Geosteering in anisotropic formations using multicomponent induction measurements |
| US7471088B2 (en) * | 2004-12-13 | 2008-12-30 | Baker Hughes Incorporated | Elimination of the anisotropy effect in LWD azimuthal resistivity tool data |
| WO2006091487A1 (en) * | 2005-02-21 | 2006-08-31 | Baker Hughes Incorporated | Well placement by use of differences in electrical anisotropy of different layers |
| US7254486B2 (en) * | 2005-04-12 | 2007-08-07 | Baker Hughes Incorporated | Method and apparatus for shale bed detection in deviated and horizontal wellbores |
| US7483793B2 (en) * | 2005-07-27 | 2009-01-27 | Baker Hughes Incorporated | Method of generating a deep resistivity image in LWD measurements |
| WO2007015992A1 (en) * | 2005-07-27 | 2007-02-08 | Baker Hughes Incorporated | Compensation for tool disposition in lwd resistivity measurements |
| KR100837910B1 (en) * | 2006-12-05 | 2008-06-13 | 현대자동차주식회사 | A device maintaining height of an active headrest |
| US8756015B2 (en) * | 2008-08-20 | 2014-06-17 | Baker Hughes Incorporated | Processing of azimuthal resistivity data in a resistivity gradient |
| US9869792B2 (en) * | 2011-10-31 | 2018-01-16 | Schlumberger Technology Corporation | Inversion-based workflow for processing nuclear density images in high-angle and horizontal wells |
| AU2013394401B2 (en) | 2013-07-18 | 2017-02-02 | Halliburton Energy Services, Inc. | Detecting boundary locations of multiple subsurface layers |
| CN104895561A (en) * | 2015-04-28 | 2015-09-09 | 中国海洋石油总公司 | Electromagnetic propagation resistivity logging joint inversion method based on parallel computing technology |
| WO2018052449A1 (en) * | 2016-09-19 | 2018-03-22 | Halliburton Energy Services, Inc. | Method of detecting substance saturation in a formation |
| CN110552689A (en) * | 2018-05-15 | 2019-12-10 | 中国石油化工股份有限公司 | method for determining distance from while-drilling instrument to stratum boundary |
| CN112677833B (en) * | 2021-01-18 | 2022-06-03 | 中车青岛四方机车车辆股份有限公司 | Rail vehicle shielding grounding system, wiring method thereof and rail vehicle |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB1588865A (en) * | 1978-05-13 | 1981-04-29 | Texaco Development Corp | Electromagnetic wave propagation well logging utilizing multiple phase shift measurement |
| US4360778A (en) * | 1978-10-30 | 1982-11-23 | Shell Oil Company | High frequency induction log for locating formation interfaces |
| US4472684A (en) * | 1980-07-24 | 1984-09-18 | Schlumberger Technology Corporation | Deep investigation induction logging with mirror image coil arrays |
| US5703773A (en) | 1996-02-08 | 1997-12-30 | Western Atlas International, Inc. | Real-time 2-dimensional inversion process and its application to induction resistivity well logging |
| US5666057A (en) | 1996-02-29 | 1997-09-09 | Western Atlas International, Inc. | Method of skin effect correction and data quality verification for a multi-frequency induction well logging instrument |
| AU707977B2 (en) | 1996-07-01 | 1999-07-22 | Shell Internationale Research Maatschappij B.V. | Electrical logging of a laminated earth formation |
| US5781436A (en) * | 1996-07-26 | 1998-07-14 | Western Atlas International, Inc. | Method and apparatus for transverse electromagnetic induction well logging |
| US5884227A (en) | 1997-04-01 | 1999-03-16 | Western Atlas International, Inc. | Method for interpreting induction logs in high resistivity contrast earth formations |
-
2000
- 2000-03-03 US US09/519,029 patent/US6308136B1/en not_active Expired - Lifetime
-
2001
- 2001-03-05 EP EP01914671A patent/EP1269222A4/en not_active Withdrawn
- 2001-03-05 WO PCT/US2001/006877 patent/WO2001067135A1/en active Application Filing
- 2001-03-05 CA CA002401898A patent/CA2401898A1/en not_active Abandoned
-
2002
- 2002-09-04 NO NO20024207A patent/NO20024207D0/en not_active Application Discontinuation
Also Published As
| Publication number | Publication date |
|---|---|
| EP1269222A4 (en) | 2008-10-01 |
| US6308136B1 (en) | 2001-10-23 |
| EP1269222A1 (en) | 2003-01-02 |
| NO20024207D0 (en) | 2002-09-04 |
| WO2001067135A1 (en) | 2001-09-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6308136B1 (en) | Method of interpreting induction logs in horizontal wells | |
| US7049821B2 (en) | Determination of borehole geometry inside cased wells with crosswell electromagnetics | |
| US6574562B2 (en) | Determination of formation anisotropy using multi-frequency processing of induction measurements with transverse induction coils | |
| US6703838B2 (en) | Method and apparatus for measuring characteristics of geological formations | |
| US6636045B2 (en) | Method of determining formation anisotropy in deviated wells using separation of induction mode | |
| GB2390432A (en) | Formation anisotropy determination while drilling, with tilted magnetic dipole antennas | |
| WO2002027357A1 (en) | Method for determination of apparent resistivities of anisotropic reservoirs | |
| US8931335B2 (en) | Processing of multi-component induction data in the presence of borehole abnormalities | |
| ZA200503381B (en) | A method for resistivity anisotropy determination in near vertical wells | |
| CA2439619C (en) | An inhomogeneous background based focusing method for multi-array induction measurements in a deviated well | |
| CA2500340A1 (en) | A method for resistivity anisotropy determination in conductive borehole environments | |
| NO324050B1 (en) | Method for determining a formation drop angle using a virtual controlled induction probe | |
| US8278930B2 (en) | Deep MWD resistivity measurements using EM shielding | |
| WO2008070239A2 (en) | Increasing the resolution of electromagnetic tools for resistivity evaluations in near borehole zones | |
| US20050083061A1 (en) | Methods and systems for estimating formation resistivity that are less sensitive to skin effects, shoulder-bed effects and formation dips | |
| WO2009052495A1 (en) | Methods for interpreting multi-component induction logs using the x-signal measurements | |
| US20050122116A1 (en) | Method and apparatus for use of the real component of a magnetic field of multicomponent resistivity measurements | |
| EP0872744B1 (en) | Determining conductivity of subterranean formations | |
| WO2008124774A1 (en) | Detection of borehole effects due to eccentricity on induction instruments with tilted transducers | |
| WO2024030779A1 (en) | Direct determination of formation apparent conductivity from em propagation measurements |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| FZDE | Discontinued |