CA2457117C  Method of determining formation anisotropy in deviated wells using separation of induction mode  Google Patents
Method of determining formation anisotropy in deviated wells using separation of induction modeInfo
 Publication number
 CA2457117C CA2457117C CA 2457117 CA2457117A CA2457117C CA 2457117 C CA2457117 C CA 2457117C CA 2457117 CA2457117 CA 2457117 CA 2457117 A CA2457117 A CA 2457117A CA 2457117 C CA2457117 C CA 2457117C
 Authority
 CA
 Grant status
 Grant
 Patent type
 Prior art keywords
 measurements
 conductivity
 horizontal
 method
 formation
 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.)
 Expired  Fee Related
Links
Classifications

 G—PHYSICS
 G01—MEASURING; TESTING
 G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS
 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
 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 welllogging
 G01V3/26—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for welllogging 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 welllogging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
Abstract
A combination of principal component measurements (301) is used to determine the horizontal resistivity (305) of the earth formations (201). The determined horizontal resistivities are used in a model for inversion of other components of the data to obtain the vertical formations resistivities (309). When multifrequency measurements are available, frequency focusing (301) is used.
Description
METHOD OF DETERMINING FORMATION ANISOTROPY IN DEVIATED
WELLS USING SEPARATION OF INDUCTION MODE
Leonty A. Tabarovsky, Mikhail I. Epov, Michael B. Rabinovich BACKGROUND OF THE INVENTION
1. Field of the Invention [0001] The invention is related generally to the field of interpretation of measurements made by well logging instruments for the purpose of determining the properties of earth formations. More specifically, the invention is related to a method for determination of anisotropic formation resistivity in a deviated wellbore using multifrequency, multicomponent resistivity data.
2. Background of the Art [0002] Electromagnetic induction and wave propagation logging tools are commonly used for determination of electrical properties of formations surrounding a borehole.
These logging tools give measurements of apparent resistivity (or conductivity) of the formation that when properly interpreted are diagnostic of the petrophysical properties of the formation and the fluids therein.
[0003] The physical principles of electromagnetic induction resistivity well logging are described, for example, in, H. G. Doll, Introduction to Induction Logging and Application to Logging of Wells Drilled with Qil Based Mud, Journal of Petroleum Technology, vol. 1, p.148, Society of Petroleum Engineers, Richardson Tex.
(1949).
Many improvements and modifications to electromagnetic induction resistivity instruments have been devised since publication of the Doll reference, supra.
Examples of such modifications and improvements can be found, for example, in U.S.
Pat. No. 4,837,517; U.S. Pat. No. 5,157,605 issued to Chandler et al, and U.S.
Pat.
No. 5,452,761 issued to Beard et al.
[0004] A limitation to the electromagnetic induction resistivity well logging instruments known in the art is that they typically include transmitter coils and receiver coils wound so that the magnetic moments of these coils are substantially parallel only to the axis of the instrument. Eddy currents are induced in the earth formations from the magnetic field generated by the transmitter coil, and in the induction instruments known in the art these eddy currents tend to flow in ground loops which are substantially perpendicular to the axis of the instrument.
Voltages are then induced in the receiver coils related to the magnitude of the eddy currents.
Certain earth formations, however, consist of thin layers of electrically conductive materials interleaved with thin layers of substantially nonconductive material. The response of the typical electromagnetic induction resistivity well logging instrument will be largely dependent on the conductivity of the conductive layers when the layers are substantially parallel to the flow path of the eddy currents. The substantially nonconductive layers will contribute only a small amount to the overall response of the instrument and therefore their presence will typically be masked by the presence of the conductive layers. The nonconductive layers, however, are the ones that are typically hydrocarbon bearing and are of the most interest to the instrument user.
Interpreting a well log made using the electromagnetic induction resistivity well logging instruments known in the art therefore may overlook some earth formations that might be of commercial interest.
[0005] The effect of formation anisotropy on resistivity logging measurements has long been recognized. Kunz and Moran studied the anisotropic effect on the response of a conventional logging device in a borehole perpendicular to the bedding plane of t thick anisotropic bed. Moran and Gianzero extended this work to accommodate an arbitrary orientation of the borehole to the bedding planes.
[0006] Rosthal (U.S. Patent No. 5,329,448) discloses a method for determining the horizontal and vertical conductivities from a propagation or induction well logging device. The method assumes 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.
[00071 Hagiwara (U.S. Patent No. 5,656,930) shows that the log response of an inductiontype logging tool can be described by an equation of the form Voc L3(2eikz(1ikL)+ikl(e'k e''`L)) (1) where fl2 = cosy 9 +,V sin 2 9 (2) and k2 = w2,u(Sh + 10h /w) (3) where L is the spacing between the transmitter and receiver, k is the wavenumber of the electromagnetic wave, 1a is the magnetic permeability of the medium, 0 is the deviation of the borehole angle from the normal to the formation, A is the anisotropy factor for the formation, c) is the angular frequency of the electromagnetic wave, o, is the horizontal conductivity of the medium and eh is the horizontal dielectric constant of the medium.
[00081 Eq. (3) is actually a pair of equations, one corresponding to the real part and one corresponding to the imaginary part of the measured signal, and has two unknowns. By making two measurements of the measured signal, the parameters k and,6can be determined. The two needed measurements can be obtained from (1) R
and X signals from induction logs, (2) phase and attenuation measurements from induction tools, (3) phase or attenuation measurements from induction tools with two different spacings, or (4) resistivity measurements at two different frequencies. In the low frequency limit, 6 can be neglected in Eq. (3) and from known values of cv and p, the conductivity o' can be determined from k, assuming a value of p equal to the permittivity of free space.
[0009] Wu (U.S. Patent 6,092,024) recognized that the solution to eqs. (1)(3) may be nonunique and showed how this ambiguity in the solution may be resolved using a plurality of measurements obtained from multiple spacings and/or multiple frequencies.
[0010] One solution to the limitation of the induction instruments known in the art is to include a transverse transmitter coil and a transverse receiver coil on the induction instrument, whereby the magnetic moments of these transverse coils is substantially perpendicular to the axis of the instrument. Such as solution was suggested in Tabarovsky and Epov, "Geometric and Frequency Focusing in Exploration of Anisotropic Seams", Nauka, USSR Academy of Science, Siberian Division, Novosibirsk, pp. 67129 (1972). Tabarovsky and Epov suggest various arrangements of transverse transmitter coils and transverse receiver coils, and present simulations of the responses of these transverse coil systems configured as shown therein.
Tabarovsky and Epov also describe a method of substantially reducing the effect on the voltage induced in transverse receiver coils which would be caused by eddy currents flowing in the wellbore and invaded zone. The wellbore is typically filled with a conductive fluid known as drilling mud. Eddy currents that flow in the drilling mud can substantially affect the magnitude of voltages induced in the transverse receiver coils. The wellbore signal reduction method described by Tabarovsky and Epov can be described as "frequency focusing", whereby induction voltage measurements are made at more than one frequency, and the signals induced in the transverse receiver coils are combined in a manner so that the effects of eddy currents flowing within certain geometries, such as the wellbore and invasion zone, can be substantially eliminated from the final result. Tabarovsky and Epov, however, do not suggest any configuration of signal processing circuitry which could perform the frequency focusing method suggested in their paper.
100111 Strack et al. (U.S. Patent 6,147,496) describe a multicomponent logging tool comprising a pair of 3component transmitters and a pair of 3component receivers.
Using measurements made at two different frequencies, a combined signal is generated having a reduced dependency on the electrical conductivity in the wellbore region. U.S. Patent 5,781,436 to Forgang et al discloses a suitable hardware configuration for multifrequency, multicomponent induction logging.
[0012] United States Patent 5,999,883 issued to Gupta et al, (the"Gupta patent") discloses a method for determination of an initial estimate of the horizontal and vertical conductivity of anisotropic earth formations. Electromagnetic induction signals induced by induction transmitters oriented along three mutually orthogonal axes are measured at a single frequency. One of the mutually orthogonal axes is substantially parallel to a logging instrument axis. The electromagnetic induction signals are measured using first receivers each having a magnetic moment parallel to one of the orthogonal axes and using second receivers each having a magnetic moment perpendicular to a one of the orthogonal axes which is also perpendicular to the instrument axis. A relative angle of rotation of the perpendicular one of the orthogonal axes is calculated from the receiver signals measured perpendicular to the instrument axis. An intermediate measurement tensor is calculated by rotating magnitudes of the receiver signals through a negative of the angle of rotation. A relative angle of inclination of one of the orthogonal axes that is parallel to the axis of the instrument is calculated, from the rotated magnitudes, with respect to a direction of the vertical conductivity.
The rotated magnitudes are rotated through a negative of the angle of inclination.
Horizontal conductivity is calculated from the magnitudes of the receiver signals after the second step of rotation. An anisotropy parameter is calculated from the receiver signal magnitudes after the second step of rotation Vertical conductivity is calculated from the horizontal conductivity and the anisotropy parameter. One drawback in the teachings of Gupta et al is that the step of determination of the relative angle of inclination of the required measurements of three components of data with substantially identical transmitterreceiver spacings. Because of limitations on the physical size of the tools, this condition is difficult to obtain;
consequently the estimates of resistivities are susceptible to error. In addition, due to the highly nonlinear character of the response of multicomponent tools, such inversion methods are time consuming at a single frequency and even more so at multiple frequencies.
100131 Analysis of the prior art leads to the conclusion that known methods of determining anisotropic resistivities in real time require very low frequencies; as a consequence of the low frequencies, the signaltonoise ratio in prior art methods is quite low.
100141 Copending United States Patent Application Ser. No. 09/825,104, published as US 2002/0173913 on November 21, 2002 (referred to hereinafter as the `104 application), and filed on April 3, 2001 teaches a computationally fast method of determination of horizontal and vertical conductivities of subsurface formations using a combination of data acquired with a transverse induction logging tool such as the 3DEXSM tool and data acquired with a conventional high definition induction logging tool (HDIL). The data are acquired in a vertical borehole. The HDIL
data are used to determine horizontal resistivities that are used in an isotropic model to obtain expected values of the transverse components of the 3DEXSM .
Differences between the model output and the acquired 3DEXSM data are indicative of anisotropy and this difference is used to derive an anisotropy factor. The method described therein has difficulties in deviated boreholes as the HDIL measurements used to derive the isotropic model are responsive to both horizontal and vertical resistivity.
[0015] There is a need for a fast and robust method of determination of anisotropic resistivity. Such a method should preferably be able to use high frequency measurements that are known to have better signaltonoise ratio than low frequency methods. The present invention satisfies this need.
SUMMARY OF THE INVENTION
[0016] A method of logging subsurface formations using data acquired with a transverse induction logging tool, the formation having a horizontal conductivity and a vertical conductivity, by obtaining a plurality of frequencies measurements indicative of vertical and horizontal conductivities in a tool referenced coordinate system. The data are transformed to a subsurface formation coordinate system.
Multifrequency focusing is applied to the measurements at a plurality of frequencies.
Horizontal formation conductivities are determined from a subset of the focused conductivity measurements. Vertical formation conductivities are determined from the focused conductivity measurements associated with the subsurface formation and the horizontal conductivities.
10016a] In one embodiment, a method of logging a subsurface formation comprising a plurality of layers each having a horizontal conductivity and a vertical conductivity is described. An electromagnetic logging tool is conveyed into a borehole in the subsurface formation, the logging tool including a plurality of transmitters and a plurality of receivers. At least one of the transmitters and at least one the receivers is inclined to an axis of the tool. The borehole has an axis inclined at a nonzero angle to a normal to the layers. The electromagnetic logging tool is used for obtaining a plurality of measurements with a plurality of pairs of the transmitters and receivers. A first subset of the plurality of measurements is used for determining a horizontal conductivity associated with each of the layers. Determined horizontal conductivities and a second subset of the plurality of measurements are used for determining a vertical conductivity associated with each of the layers.
[0017] In a preferred embodiment of the invention, a transformation independent of the formation azimuth may be used to determine the conductivity of the transversely anisotropic formation.
BRIEF DESCRIPTION OF THE FIGURES
[0018] FIG.1(Prior art) shows an induction instrument disposed in a wellbore penetrating earth formations.
FIG. 2 shows the arrangement of transmitter and receiver coils in a preferred embodiment of the present invention marketed under the name 3DEX
FIG. 3 shows an earth model example used in the present invention.
FIG. 4 is a flow chart illustrating steps comprising the present invention.
7a DETAILED DESCRIPTION OF THE INVENTION
[0019] FIG. 1 (prior art) shows an induction well logging instrument 10 disposed in a wellbore 2 penetrating earth formations. The earth formations are shown generally at 6, 8, 12 and 14. The instrument 10 is typically lowered into the wellbore 2 at one end of an armored electrical cable 22, by means of a winch 28 or similar device known in the art. An induction well logging instrument which will generate appropriate types of signals for performing the process of the present invention is described, for example, in U.S. Pat. No. 5,452,761 issued to Beard et al. The prior art induction logging tool includes a transmitter coil and a plurality of receiver coils 18A18F. The coils in the prior art device are oriented with the coil axes parallel to the axis of the tool and to the wellbore.
[0020] Turning now to Fig. 2, the configuration of transmitter and receiver coils in a preferred embodiment of the 3DExplorer TM induction logging instrument of Baker Hughes is disclosed. Such a logging instrument is an example of a transverse induction logging tool. Three orthogonal transmitters 101,103 and 105 that are referred to as the T,,, Tz, and T. transmitters are shown (the z axis is the longitudinal axis of the tool). Corresponding to the transmitters 101, 103 and 105 are receivers 107,109 and 111, referred to as the R., R., and RY receivers, for measuring the corresponding components (H,.,, 14n., Hoof induced signals. In addition, crosscomponents are also measured. These are denoted by H,ryõ H,y etc.
[0021] Fig. 3 is a schematic illustration of the model used in the present invention.
The subsurface of the earth is characterized by a plurality of layers 201a, 201b, ..., 2011. T h e l a y e r s h a v e t h i c k n e s s e s d e n o t e d by hl, h2, ... h,. The horizontal and vertical resistivities in the layers are denoted by Rw, Rh2, ...Rh; and R,1, RO, ...Ry, respectively. Equivalently, the model may be defined in terms of conductivities (reciprocal of resistivity). The borehole is indicated by 202 and associated with each of the layers are invaded zones in the vicinity of the borehole wherein borehole fluid has invaded the formation and altered is properties so that the electrical properties are not the same as in the uninvaded portion of the formation. The invaded zones have lengths L,,01, L.x02, = = =L 01 extending away from the borehole. The resistivities in the invaded zones are altered to values R.z01, R.x02, . . R 01. In the embodiment of the invention discussed here, the invaded zones are assumed to be isotropic while an alternate embodiment of the invention includes invaded zones that are anisotropic, i.e., they have different horizontal and vertical resistivities. It should further be noted that the discussion of the invention herein may be made in terms of resistivities or conductivities (the reciprocal of resistivity). The z component of the 3DEXSM
tool is oriented along the borehole axis and makes an angle 0 (not shown) with the normal to the bedding plane. The x component of the tool makes an angle q5 with the "up"
direction.
[0022] In the `104 application to Tabarovsky, et al. multifrequency, multicomponent induction data are obtained using, for example, the 3DEXTM tool, and a multifrequency focusing is applied to these data. As disclosed in United States Patent 5,703,773 to Tabarovsky et al., the contents of which are fully incorporated herein by reference, the response at multiple frequencies may be approximated by a Taylor series expansion of the form:
6 (0)1) CO 1/2 CO 3/2 ... w n/2 s a 1 1 1 0 1/2 3/2 n/1 0a(0)2) 1 002 001 ... w2 s1/2 (4) 1/2 3/2 n/2 a (&m_1) 1 ~m1 com1 ..= o) n/2 "S
6 (n1)/2 1/2 3/2 n/2 6a C0m 1 CVm 00m ... 0) m Sn/2 [0023] In a preferred embodiment of the invention of the `104 application, the number m of frequencies co is ten. In eq.(4), n is the number of terms in the Taylor series expansion. This can be any number less than or equal to in. The coefficient s312 of the cv312 term (cv being the square of k, the wave number) is generated by the primary field and is relatively unaffected by any inhomogeneities in the medium surround the logging instrument, i.e., it is responsive primarily to the formation parameters and not to the borehole and invasion zone. In fact, the coefficient s312 of the &3/2 term is responsive to the formation parameters as though there were no borehole in the formation. Specifically, these are applied to the H and HI, components. Those versed in the art would recognize that in a vertical borehole, the H and Hn, would be the same, with both being indicative of the vertical conductivity of the formation. In one embodiment of the invention, the sum of the H and H , is used so as to improve the signal to noise ratio (SNR). This multifrequency focused measurement is equivalent to the zero frequency value. As would be known to those versed in the art, the zero frequency value may also be obtained by other methods, such as by focusing using focusing electrodes in a suitable device.
[0024] Along with the 3DEXT`, the method of the `104 application also uses data from a prior art High Definition Induction Logging (HDIL) tool having transmitter and receiver coils aligned along the axis of the tool. These data are inverted using a method such as that taught by Tabarovsky et al, or by U.S. Patent 5,884,227 to Rabinovich et al., the contents of which are fully incorporated herein by reference, to give an isotropic model of the subsurface formation. Instead of, or in addition to the inversion methods, a focusing method may also be used to derive the initial model.
Such focusing methods would be known to those versed in the art and are not discussed further here. As discussed above, the HDIL tool is responsive primarily to the horizontal conductivity of the earth formations when run in a borehole that is substantially orthogonal to the bedding planes. The inversion methods taught by Tabarovsky et al and by Rabinovich et al are computationally fast and may be implemented in real time. These inversions give an isotropic model of the horizontal conductivities (or resistivities) [0025] Using the isotropic model derived, a forward modeling is used in the `104 application to calculate a synthetic response of the 3DE)(m tool at a plurality of frequencies. A suitable forward modeling program for the purpose is disclosed in Tabarovsky and Epov "Alternating Electromagnetic Field in an Anisotropic Layered Medium" Geol. Geoph., No. 1, pp. 101109. (1977). Multifrequency focusing is then applied to these synthetic data . In a preferred embodiment of the invention of the '104 application, the method taught by Tabarovsky is used for the purpose.
[0026] In the absence of anisotropy, the output from model output should be identical to the multifrequency focused measurements. Denoting by o;,, the multifrequency focused transverse component synthetic data from and by nQas the multifrequency focused field data from, the anisotropy factor A is then calculated in the `104 application.
[0027] The H,, for an anisotropic medium is given by H TL 7 + 3+ 1 8 (5) v h where 91 Fo ' gh 2 a /~ = 6h D/J Oh 6v For a threecoil subarray, H 1 (,+T,)( ev,u ~h M (6) 4Tr 3 2 Upon introducing the apparent conductivity for H this gives 3/2 3(1 1 3/2 ~meas = + h or (v.3/2 meas  6iso ) = 6h 4 + 4~  1  Ch 4A
which gives the result _ 1 A 3/2 3/2 (7) 4 iso  2 n:eas 1 3 a. 3/2 i where o, is the conductivity obtained from the HDIL data, i.e., the horizontal conductivity. The vertical conductivity is obtained by dividing u, by the anisotropy factor from eq. (5). An important aspect of the `104 application is that in a vertical borehole, the measurements made by a HDIL tool depend only on the horizontal conductivities and not on the vertical resistivities. The method of the invention disclosed there is to obtain an isotropic model from the HDIL data, use the isotropic model to predict measurements made on other components and to use a difference between the predicted and actual measurements to obtain the vertical conductivity.
[00281 In a similar manner, the method of the present invention can be viewed as finding a combination of 3DEXSM measurements (called modes of the induction measurements) that are responsive only to the horizontal conductivity, deriving a model of horizontal conductivity from this combination of measurements, predicting values of other components of 3DEXSM measurements and using a difference between these predicted measurements and the actual measurements to determine a vertical conductivity. In a particular embodiment of the present invention, the desired , combination includes only the principal component measurements, i.e., upon H, H1õ and HLZ . The flow chart of the method of the present invention is shown in Fig.
[00291 The method of the present invention starts with an estimate of the dip and azimuth of the formation relative to the borehole axis 300. These angles are defined below. In addition, a sensor on the logging tool also provides another angle measurement called the "toolface angle" that is also used in the analysis of the data.
Multifrequency focused data are derived from multifrequency measurements' 301.
The data are transformed 303 as discussed below to give measurements that are indicative only of horizontal conductivity. These data are inverted 305 to give a model of the horizontal conductivity of the data. These estimates of horizontal conductivity are used in an isotropic model as estimates of the vertical conductivity 307. The measured data are then inverted using this initial estimate of vertical conductivities 309. A check is made of the goodness of fit (difference between a model output based on the inverted model and the actual measurements) 311. If the difference is below a predetermined threshold, then the model is accepted. If the difference is excessive, a iterative procedure is carried out with an updated angle estimate 313 until the result is acceptable. Any suitable iterative procedure may be used such as that based on a gradient method or a method of steepest descent.
Such iterative methods would be known to those versed in the art and are not discussed further.
[0030]At this point we develop the principal component structure for measuring formation anisotropy in bedding planes when the borehole is not normal (perpendicular) to the bedding plane. For simplifying the notation, and to avoid confusion, the x, y and z components in the tool coordinates are called hereafter the (1,2,3) coordinate system. The x, y and z components in the earth coordinate system will be referred to as such.
[0031] In the tool coordinate system, the matrix of magnetic components, HT, may be represented in the following form:
X11 h12 13 HT = "21 122 ]23 (8) h1 "32 "33 For layered formations, the matrix, HT, is symmetric. We measure three diagonal elements, hl,, h22, and h33. The nondiagonal elements are not needed in the present invention.
[00321 In the earth coordinate system, {x, y, z}, associated with the plane formation boundaries (zaxis is perpendicular to the boundaries and directed downwards )the magnetic matrix may be presented as follows:
hxx hxy hxz HM = hX, hn, hyZ (9) hxh hyZ hZZ
The formation resistivity is described as a tensor, p. In the coordinate system associated with a formation, the resistivity tensor has only diagonal elements in the absence of azimuthal anisotropy:
pt 0 0 p= 0 p, 0 (10) 0 0 p, p, =P.. =p,,, Pn = PZZ
[00331 The "tool coordinate " system (1, 2, 3) can be obtained from the "formation coordinate" system (x, y, z) as a result of two sequential rotations:
(1) Rotation about the axis "2" by the angle 0, such that the axis "3" in a new position (let us call it "3' ") becomes parallel to the axis z of the "tool"
system;
(2) Rotation about the axis "3' " by the angle (p, such that the new axis "1"
(let us call it "1' ") becomes parallel to the axis x of the tool system.
In the present invention, an iterative procedure as shown in Fig. 4 is used for 0 and 4.
0 is the relative inclination of the borehole axis to the normal to the bedding while is the azimuth. An initial estimate for 0 and cb may be determined from borehole surveys or from resistivity imaging devices and from knowledge of the toolface angle.
:
[0034] The first rotation is described using matrices 8 and 0' C0 0 S C0 0  Se 8= 0 1 0 eT = 0 1 0 (11)  SB 0 CB so 0 CO
Here, CO = cos 0, So = sin 9 The second rotation is described using matrices ip and gyp':
(p = S9 C9 0 , CpT =  SP C9 0 (12) Here, C9 = cos cp, SP = sin cp Matrices HM (the formation coordinate system) and HT (the tool coordinate system) are related as follows:
HT = RT Hn,R (13) RT =CpT~T R =60 (14) It is worth noting that the matrix HM contains zero elements:
hXY = hyx = 0 (15) This is true for multifrequency focused measurements as described below. It is also important to note that the following three components of the matrix HM depend only on the horizontal resistivity.
h. = f (ot) , 11yz = fYz(Pt) I hzz = fzz(Pt) (16) [0035] Two remaining elements depend on both horizontal and vertical resistivities.
hxx =f (P, pa ), hyy  .f,Y (Pt 9 Pn) (17) Taking into account Equations (11), (12), (14) and (15), we can rewrite Equation (13) as follows:
1111 412 113 C9 Sq, 6 "'CO 0 S9h 0 hxz C0 0 SO C9 S9, 0 h11 k2 k3 = S, C9, 0 0 1 0 0 ho hyz 0 1 0 S. C. 0 (18) Ai 1132 1133 0 0 1S0 0 CO i hyz hzz SB 0 CB 0 0 1 The following expanded calculations are performed in order to present Equation (18) in a form more convenient for analysis.
C0 0 S Cc,  Sc, 0 CBCc  C Sc Se Al = 0 1 0 Sc, Cc, 0 = Sc, Cc, 0  so 0 CB 0 0 1  S Cc SBSc C
hxx 0 hxZ C Cc  CBSc so A2 = 0 hyy hyZ Sc Cc, 0 =
hxZ hyZ hZZ  SBCc SBSc CB
CBCc,hxx SBCq,hxZ CBSc,hxx +S Sc,hxZ S hxx +C hxZ
Sc,hyy SBCc,hyZ Cc,hyy +SBSc,hyZ C hyZ
C CphxZ +S ,hyz SBCc,h2Z CBSc,h. +Cc,hyZ +SBSch. Seh. +CehZZ
C0 0 S0 C Cc,hxx S CChxZ CBSc,hxx +SBSqhxZ Soh= +CehxZ
A3 = 0 1 0 Sc,hyy SoCv hyZ Cghyy +SBSc,hyZ C hyZ
SB 0 CB C Cc,hxZ +SghyZ SBCc,h. CBSc,hxZ +Cc,hyz +SBSc,h= SehxZ +C9h2Z
The components of A 3 are given as 41 CBCc,hxx C SBCc,hxZ CBS Cc,h,. S Sc,hyz +SeCc,hZZ
1 =C02 Cchxx 2CBSBCchxz S SchyZ +S Cchzj] (*) a12) =COSch,, +CBSBSchxZ +C SBCc,h: SBCc,hyz S Sc,hzz [a12 =CBSc,hxx +2C SBSc,hxZ S CchyZ S Sch2Z] (*) a13) = C9 Sehxx + C hZ  SehxZ  C9S9hZZ
{ai3 = CeSeh + (CB  Se )h.xZ  CeSehZZ } (*) 3) =S,h,,3, S C,h'] (*) [a( [a(3 =C9h3,y +SOSVh 22 (*) [a(3)=C h ] (*) 23 B yz a(3i) =C0S C9h. S Cqh.Z +CBCQ,h. +C Sehy7 C SeCLhZZ
[a3i) =CBSBCQ,h,, +(C 8 02 )Cq,hxz +CeSehyZ C SeC~hZZ~ (*) 3) = CeSeSPh + SBSphxZ  C02S9hZ + COC,hYZ + CBSeS 9hZZ
a( [a(3) =COSOS~hxx (CB SB)S,hXZ +CBC,hyZ +CBSOSq,hZZ] (*) a33) = Seh, + COSeh,Z + COSeh. + Cehz.
[a(3) = SB hxx + 2C0Seh +C02 h . ~ (*) [00361 Taking into account all the above calculations, we can represent Equation (18) in the following form:
h11 h12 h13 Cq, 'sip 0 a 3 ll a 3 12 a 3 h21 h22 h23 =  S~ CP 0 a21 a22 a23 h31 h32 h33 0 0 1 a 3 31 a 3 32 a 3 The method of the present invention involves defming a linear combination of the measurements that are responsive substantially to the horizontal conductivity and not responsive to the vertical conductivity. In a preferred embodiment of the invention, the linear combination is that of measurements h,,, h22, and h33 (i.e., the principal components only), although in alternate embodiments of the invention, a linear combination of any of the measurements may be used. The example given below is that of the preferred embodiment.
[0037] Let us consider expressions for the measured principal components, h,,, h22, and h33:
h = a(3)C + a(3)S
ii u 9 21 h22 = a1z)S + a23)C (19) h =a(3) More detailed representation yields:
hll = C02C2h  2CBSBC2hx. So CPSPhYZ + S02C2hZZ + S2h3ry  S0CPSq,I2 [hll = C02C2hXX + S2h3ry  2CBSBC2hxZ  2SBC9SPhyZ + S02C2hZZ 1 (20) h22 = C0 Soh,  2C SBSPhxZ + SBCPSph3,Z +S0 Phzz +C2 h3ry + S CPSPh3,Z
[h22 = C02S2h,, +C2 hYY  2COSOS2h z + 2SBCPSPh3Z + SBS2hZZ
(21) [h33 =S h.x +2C Sehz. +C0hZZ1 (22) [0038] Expressions for each component, h,,, h22, and h33, contain two types of functions: some depending only on pl, and some others depending on both, p, and p,,.
Equations (13)(15) may be represented in the following form:
h11 =C0Cgh +SVhY5, +f11(Pt) h22 =CBSPh +Cq,hyy +f22(Pt) (23) h33 = SB hx, + f33 (Pt ) Here, fll (Pt) _ 2C SBC2h  2S CTSPhyZ + S2C2hZZ
f22(A)=2C S,, S2hx, +2S CPSq,hyZ +SeS2hZ (24) f33(Pt)=2C S0h. +C2hZZ
A linear combination of Equations (23) is defined in the form:
h = ahll + Nh22 + h33 (25) [0039] Detailed consideration of Equation (25) yields:
h=oC CChh+aSp2hyy+c1(Pt)+/X'BS~h +/ ~hyy+/~J22(Pt)+Seh +f331F t) h=(aCBC. +N"0 s9 +SB)hxx (aS2 P +/"'PY`yy +a l1(Pt)+ff22(Pt)+f33(Pt [0040] The method of the present invention involves defining the coefficients, a and 8, in such a way that the resulting linear combination, h, does not depend on the vertical resistivity. To achieve that, we need to null the following part of the expression for h:
hf = (aC02C2 +/3C02B S2 + S )ham + (aS2 + fC2 )h~ , = 0 (26) Imposing the following conditions satisfies equation (26):
1aCeC2 + f3C02S2 +SB =0 V 4p as OP + 8cv = 0 (27) [00411 Let us calculate the coefficients, a and. The second Equation in (27) yields:
P  C 2 a (28) Ip After substitution of Equation (28) in the first Equation of (27), we obtain:
aC C~  C 2 a CB S2 + S = 0 = aC C~   , + SB = aCe ~C 2 ' _ SB
( 2 2 z 2 aC `C: + S XC:  SP S 2 aC C22 = SB
p C+2 S2 a a (29) C2q C0 To obtain the coefficient, /3, let us substitute Equation (29) in Equation (28):
Jo _ 0 = 0 (30) C:C29C0 2 2 [0042] Finally, a =  S0 z C2V C0 (31) C 2q C 0 [0043] It is convenient to normalize coefficients, a and /3. Let us introduce a normalization factor, x:
K = 1+a2 +/32 (32) Equation (25) may be presented in the form:
hf =a'h11 +18h22 +77133 (33) Here, h f, = h f /K, a' = a/K, fi' _ /3/K, Y' = Y /K = (34) Calculations yield:
LC4 A4 S4 S4 C4 +S4 S4 C2~C4 + (C: + s~ )se K2 =1+ ~~ c4 + 0V 4 1 + ~2 C4 = C2 C4 2p 0 2p 0 2q' B 2p 0 c2,c + C4 +S4 se K = 2 (35) C2,C0 Consequently, Y,  C2r,C9 Cep 0 C: + S~ 0 Cp S0 C2,C0 C
a = C2,, Ce VCq,C +C:+S: B Cp 0+C:+S:
2 2 2 s2 s2 Sp s0 C2~pC0 _ Sp 0 C2 , Ce V Cq,C + C~ + Sp e C2C + C4 + S: e [0044] Finally, we obtain:
a ' _  _ 9 K ' 9 (36) K ' Y~ = C2g) 0 K ' Here, K'= C 2 C + C+ S S 0 (37) The coefficient, K, degenerates under the following conditions:
0=0, (p=% =:> K'=0 (38) [0045] Using the derivation given above, for an estimated value of 0 and the conductivities may be derived. A difference between the model output and the measured values may then be used in the iterative procedure described with respect to Fig. 4.
[0046] The derivation above has been done for a single frequency data.
Multifrequency Focused (MFF) data is a linear combination of single frequency measurements so that the derivation given above is equally applicable to MFF
data. It can be proven that the three principle 3DEX' measurements, MFF (multifrequency focusing) processed, may be expressed in the following form:
MFF (ham ) MFF (hl,) a, a2 a3 a4 ) MFF (h22) = bt b2 b3 b4 MFF (h Y, 1lIF'F (h ) (39) MFF (h33) C1 C2 C3 C4 zz MFF (hx_) The matrix coefficients of Equation 39 depend on 9,, 0, and three trajectory measurements: deviation, azimuth and rotation.
[0047] The components of the vector in the right hand side of Equation 39 represent all nonzero field components generated by three orthogonal induction transmitters in the coordinate system associated with the formation. Only two of them depend on vertical resistivity: h and h.,Y. This allows us to build a linear combination of measurements, h11, h22, and h33 in such a way that the resulting transformation 'depends only on h~Z and hX;, or, in other words, only on horizontal resistivity. Let T be the transformation with coefficients aa, 8 and y:
T= )+/M'l("22)+ 1 '33) (40) [0048] The coefficients a, Band ymust satisfy the following system of equations:
ala+bl/3+c1y=0 a2a + b2fl + c2y = 0 (41) a2 +,132 + y2 = 1 [0049] From the above discussion it follows that a transformation may be developed that is independent of the formation azimuth. The formation azimuthindependent transformation may be expressed as:
T,, =(h,l +h22)Sin 20h33(1+cos20) (42) where 0 is the dip of the formation and To is the linear transformation to separate modes. With this transformation and the above series of equations we may determine the conductivity of the transversely anisotropic formation.
[0050] It is to be noted, however, that when the earth formation is uniform (i.e., there are no formation boundaries within the region of investigation of the tool), it is not possible to satisfy eq. (40). It is necessary then to have a measurement of at least one crosscomponent.
[0051]The present invention has been discussed above with respect to measurements made by a transverse inductionlogging tool conveyed on a wireline. This is not intended to be a limitation and the method is equally applicable to measurements made using a comparable tool conveyed on a measurementwhiledrilling (MWD) assembly or on coiled tubing.
[0052] While the foregoing disclosure is directed to the preferred embodiments of the invention, various modifications will be apparent to those skilled in the art.
It is intended that all variations within the scope and spirit of the appended claims be embraced by the foregoing disclosure.
Claims (14)
(a) conveying an electromagnetic logging tool into a borehole in the subsurface formation, said logging tool including a plurality of transmitters and a plurality of receivers at least one of said transmitters and at least one of said receivers inclined to an axis of the tool, said borehole having an axis inclined at a nonzero angle to a normal to said layers;
(b) using said electromagnetic logging tool for obtaining a plurality of measurements with a plurality of pairs of said transmitters and receivers;
(c) using a first subset of said plurality of measurements for determining a horizontal conductivity associated with each of said layers; and (d) using determined horizontal conductivities and a second subset of said plurality of measurements for determining a vertical conductivity associated with each of said layers.
Priority Applications (5)
Application Number  Priority Date  Filing Date  Title 

US31265501 true  20010815  20010815  
US60/312,655  20010815  
US10072173 US6636045B2 (en)  20010403  20020207  Method of determining formation anisotropy in deviated wells using separation of induction mode 
US10/072,173  20020207  
PCT/US2002/025787 WO2003016950A1 (en)  20010815  20020814  Method of determining formation anisotropy in deviated wells using separation of induction mode 
Publications (2)
Publication Number  Publication Date 

CA2457117A1 true CA2457117A1 (en)  20030227 
CA2457117C true CA2457117C (en)  20101019 
Family
ID=26753071
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CA 2457117 Expired  Fee Related CA2457117C (en)  20010403  20020814  Method of determining formation anisotropy in deviated wells using separation of induction mode 
Country Status (5)
Country  Link 

US (1)  US6636045B2 (en) 
EP (1)  EP1425611A4 (en) 
CA (1)  CA2457117C (en) 
GB (1)  GB2400672B (en) 
WO (1)  WO2003016950A1 (en) 
Families Citing this family (40)
Publication number  Priority date  Publication date  Assignee  Title 

FR2836557B1 (en) *  20020228  20040528  Schlumberger Services Petrol  Method and geophysical prospecting device of a porous earth formation containing at least one fluid add electrolytes 
US6998844B2 (en) *  20020419  20060214  Schlumberger Technology Corporation  Propagation based electromagnetic measurement of anisotropy using transverse or tilted magnetic dipoles 
US7038457B2 (en) *  20020729  20060502  Schlumberger Technology Corporation  Constructing colocated antennas by winding a wire through an opening in the support 
US7286091B2 (en)  20030613  20071023  Schlumberger Technology Corporation  Colocated antennas 
GB2410336B (en) *  20020927  20060308  Baker Hughes Inc  A method for resistivity anisotropy determination in conductive borehole environments 
CA2500337A1 (en) *  20020927  20040408  Baker Hughes Incorporated  A method for resistivity anisotropy determination in near vertical wells 
US7436183B2 (en) *  20020930  20081014  Schlumberger Technology Corporation  Replaceable antennas for wellbore apparatus 
US7031839B2 (en) *  20021115  20060418  Baker Hughes Incorporated  Multifrequency focusing for MWD resistivity tools 
US7274991B2 (en) *  20040615  20070925  Baker Hughes Incorporated  Geosteering in anisotropic formations using multicomponent induction measurements 
US7269515B2 (en) *  20040615  20070911  Baker Hughes Incorporated  Geosteering in anisotropic formations using multicomponent induction measurements 
US6819111B2 (en) *  20021122  20041116  Baker Hughes Incorporated  Method of determining vertical and horizontal resistivity, and relative dip in anisotropic earth formations having an arbitrary electromagnetic antenna combination and orientation with additional rotation and position measurements 
US7382135B2 (en) *  20030522  20080603  Schlumberger Technology Corporation  Directional electromagnetic wave resistivity apparatus and method 
US7027922B2 (en) *  20030825  20060411  Baker Hughes Incorporated  Deep resistivity transient method for MWD applications using asymptotic filtering 
US7043370B2 (en) *  20030829  20060509  Baker Hughes Incorporated  Real time processing of multicomponent induction tool data in highly deviated and horizontal wells 
US7737697B2 (en)  20031203  20100615  Baker Hughes Incorporated  Method and apparatus for use of the real component of a magnetic field of multicomponent resistivity measurements 
US7336080B2 (en) *  20031203  20080226  Baker Hughes Incorporated  Method and apparatus for use of the real component of a magnetic field of multicomponent resistivity measurements 
US8112227B2 (en) *  20040615  20120207  Baker Hughes Incorporated  Processing of multicomponent induction measurements in a biaxially anisotropic formation 
US7392137B2 (en) *  20040615  20080624  Baker Hughes Incorporated  Determination of formation anistrophy, dip and azimuth 
CA2586680A1 (en) *  20041104  20060518  Baker Hughes Incorporated  Multiscale multidimensional well log data inversion and deep formation imaging method 
US7496451B2 (en) *  20060306  20090224  Baker Hughes Incorporated  Real time data quality control and determination of formation angles from multicomponent induction measurements using neural networks 
US7268555B1 (en)  20060406  20070911  Baker Hughes Incorporated  Correction of crosscomponent induction measurements for misalignment using comparison of the XY formation response 
US7333891B2 (en) *  20060406  20080219  Baker Hughes Incorporated  Correction of crosscomponent induction measurements for misalignment using comparison of the XY formation response 
WO2007117631A3 (en)  20060406  20080502  Baker Hughes Inc  Correction of crosscomponent induction measurements for misalignment using comparison of the xy formation response 
US7379818B2 (en) *  20060406  20080527  Baker Hughes Incorporated  Correction of crosscomponent induction measurements for misalignment using comparison of the XY formation response 
US8931335B2 (en)  20060407  20150113  Baker Hughes Incorporation  Processing of multicomponent induction data in the presence of borehole abnormalities 
US7629791B2 (en) *  20060801  20091208  Baker Hughes Incorporated  Method and apparatus for making multicomponent measurements in deviated wells 
US7778778B2 (en) *  20060801  20100817  Baker Hughes Incorporated  Correction of multicomponent measurements for tool eccentricity in deviated wells 
US20080224706A1 (en) *  20061113  20080918  Baker Hughes Incorporated  Use of Electrodes and MultiFrequency Focusing to Correct Eccentricity and Misalignment Effects on Transversal Induction Measurements 
US7554328B2 (en) *  20061113  20090630  Baker Hughes Incorporated  Method and apparatus for reducing borehole and eccentricity effects in multicomponent induction logging 
KR100837910B1 (en) *  20061205  20080613  현대자동차주식회사  A device maintaining height of an active headrest 
US7898260B2 (en) *  20070410  20110301  Baker Hughes Incorporated  Method and apparatus for detecting borehole effects due to eccentricity of induction instruments 
US8046170B2 (en) *  20080903  20111025  Baker Hughes Incorporated  Apparatus and method for estimating eccentricity effects in resistivity measurements 
US8614578B2 (en) *  20090618  20131224  Schlumberger Technology Corporation  Attenuation of electromagnetic signals passing through conductive material 
US20110051552A1 (en) *  20090825  20110303  Schlumberger Technology Corporation  Methods and apparatus to calculate a distance from a borehole to a boundary of an anisotropic subterranean rock layer 
US9229125B2 (en) *  20090917  20160105  Baker Hughes Incorporated  TDEM forward focusing system for downhole use 
US8466682B2 (en) *  20090929  20130618  Schlumberger Technology Corporation  Apparatus and method for downhole electromagnetic measurement while drilling 
WO2014098838A1 (en) *  20121219  20140626  Halliburton Energy Services, Inc.  Method and apparatus for optimizing deep resistivity measurements with multicomponent antennas 
WO2014143015A1 (en) *  20130315  20140918  Halliburton Energy Services Inc.  Dip correction for array induction tool data 
EP3063569A1 (en) *  20131030  20160907  Halliburton Energy Services, Inc.  Apparatus and method of processing multicomponent induction data 
CN104597510B (en) *  20150129  20170308  中国科学院电子学研究所  A method for obtaining a frequencydomain electromagnetic system tensor rotational invariants 
Family Cites Families (32)
Publication number  Priority date  Publication date  Assignee  Title 

US4531187A (en) *  19821021  19850723  Uhland Joseph C  Game monitoring apparatus 
DE3505793A1 (en) *  19840221  19850822  Komatsu Mfg Co Ltd  Object identification method 
JPS62126487A (en) *  19851127  19870608  Shinko Electric Co Ltd  Edge detecting circuit for binary image 
DE3542484C2 (en) *  19851130  19870924  Ant Nachrichtentechnik Gmbh, 7150 Backnang, De  
JP2556484B2 (en) *  19861020  19961120  キヤノン株式会社  Graphic contour extracting device 
US4725879A (en) *  19870327  19880216  Honeywell Inc.  Chroma responsive inspection apparatus selectively producing analog voltage levels based on the luminance, the phase of the chrominance subcarrier, or the amplitude of the chrominance subcarrier 
US5115477A (en) *  19880331  19920519  Honeywell Inc.  Image recognition edge detection method and system 
US4969202A (en) *  19880331  19901106  Honeywell Inc.  Image recognition edge detection method and system 
CA1338342C (en) *  19880608  19960521  Mitsuru Kurita  Image processing apparatus 
US4991223A (en) *  19880630  19910205  American Innovision, Inc.  Apparatus and method for recognizing image features using color elements 
US5216498A (en) *  19890322  19930601  Konica Corporation  Image processing apparatus capable of detecting marked region 
JPH0372778A (en) *  19890811  19910327  Fuji Xerox Co Ltd  Area identification system for picture processing unit 
JPH07106646B2 (en) *  19891225  19951115  富士ゼロックス株式会社  Image processing apparatus 
US5227871A (en) *  19901130  19930713  Canon Kabushiki Kaisha  Image processing apparatus capable of discriminating a predetermined image 
US5258837A (en) *  19910107  19931102  Zandar Research Limited  Multiple security video display 
JP2585874B2 (en) *  19910318  19970226  富士通株式会社  Color image of the edge detection method 
US5329448A (en)  19910807  19940712  Schlumberger Technology Corporation  Method and apparatus for determining horizontal conductivity and vertical conductivity of earth formations 
US5375177A (en) *  19910927  19941220  E. I. Du Pont De Nemours And Company  Method of identifying and characterizing a valid object by color 
JP3073599B2 (en) *  19920422  20000807  本田技研工業株式会社  Image edge detection device 
US5343538A (en) *  19921002  19940830  International Remote Imaging Systems, Inc.  Method and an apparatus for identifying an object using quantile partitions 
US5757191A (en) *  19941209  19980526  Halliburton Energy Services, Inc.  Virtual induction sonde for steering transmitted and received signals 
US5656930A (en)  19950206  19970812  Halliburton Company  Method for determining the anisotropic properties of a subterranean formation consisting of a thinly laminated sand/shale sequence using an induction type logging tool 
US5703773A (en)  19960208  19971230  Western Atlas International, Inc.  Realtime 2dimensional inversion process and its application to induction resistivity well logging 
US5666057A (en)  19960229  19970909  Western Atlas International, Inc.  Method of skin effect correction and data quality verification for a multifrequency induction well logging instrument 
WO1998000733A1 (en)  19960701  19980108  Shell Internationale Research Maatschappij B.V.  Electrical logging of a laminated earth formation 
US6044325A (en) *  19980317  20000328  Western Atlas International, Inc.  Conductivity anisotropy estimation method for inversion processing of measurements made by a transverse electromagnetic induction logging instrument 
US5781436A (en)  19960726  19980714  Western Atlas International, Inc.  Method and apparatus for transverse electromagnetic induction well logging 
US5884227A (en)  19970401  19990316  Western Atlas International, Inc.  Method for interpreting induction logs in high resistivity contrast earth formations 
US6092024A (en)  19970403  20000718  Baker Hughes Incorporated  Method and apparatus for determining resistivity and dielectric anisotropy parameters of earth formations by using multifrequency and/or multispacing measurements 
US6466872B1 (en) *  19991108  20021015  Baker Hughes Incorporated  Method for determination of apparent resistivities of anisotropic reservoirs 
US6502036B2 (en) *  20000929  20021231  Baker Hughes Incorporated  2D inversion of multicomponent induction logging data to resolve anisotropic resistivity structure 
US6584408B2 (en)  20010626  20030624  Schlumberger Technology Corporation  Subsurface formation parameters from triaxial measurements 
Also Published As
Publication number  Publication date  Type 

EP1425611A1 (en)  20040609  application 
EP1425611A4 (en)  20121121  application 
US6636045B2 (en)  20031021  grant 
WO2003016950A1 (en)  20030227  application 
US20020186013A1 (en)  20021212  application 
GB2400672B (en)  20050622  grant 
CA2457117A1 (en)  20030227  application 
GB0403497D0 (en)  20040324  grant 
GB2400672A (en)  20041020  application 
Similar Documents
Publication  Publication Date  Title 

US6556015B1 (en)  Method and system for determining formation anisotropic resistivity with reduced borehole effects from tilted or transverse magnetic dipoles  
US6534986B2 (en)  Permanently emplaced electromagnetic system and method for measuring formation resistivity adjacent to and between wells  
US5703773A (en)  Realtime 2dimensional inversion process and its application to induction resistivity well logging  
US5260661A (en)  Calibrating and compensating influence of casing thickness variations on measurements of low frequency A.C. magnetic fields within cased boreholes to determine properties of geological formations  
Wang et al.  3D electromagnetic anisotropy modeling using finite differences  
US7057392B2 (en)  Method and apparatus for directional resistivity measurement while drilling  
US4278941A (en)  High frequency induction log for determining resistivity and dielectric constant of the earth  
US7046009B2 (en)  Method for measuring transient electromagnetic components to perform deep geosteering while drilling  
US5184079A (en)  Method and apparatus for correcting data developed from a well tool disposed at a dip angle in a wellbore to eliminate the effects of the dip angle on the data  
US20030229450A1 (en)  Induction well logging apparatus and method  
Li et al.  New directional electromagnetic tool for proactive geosteering and accurate formation evaluation while drilling  
US6219619B1 (en)  Inhomogeneous backgroundbased software focusing method for arraytype induction logging tools  
US6556016B2 (en)  Induction method for determining dip angle in subterranean earth formations  
US7382135B2 (en)  Directional electromagnetic wave resistivity apparatus and method  
US6809521B2 (en)  Apparatus and method for wellbore resistivity measurements in oilbased muds using capacitive coupling  
US6734675B2 (en)  Apparatus accurately measuring properties of a formation  
US6393364B1 (en)  Determination of conductivity in anisotropic dipping formations from magnetic coupling measurements  
US20090248308A1 (en)  Simultaneous inversion of induction data for dielectric permittivity and electric conductivity  
US7038457B2 (en)  Constructing colocated antennas by winding a wire through an opening in the support  
US20050083063A1 (en)  Electromagnetic method for determining dip angles independent of mud type and borehole environment  
US4258322A (en)  Electromagnetic subsoil prospecting process using an asymptotic low frequency range  
US6304086B1 (en)  Method and apparatus for evaluating the resistivity of formations with high dip angles or highcontrast thin layers  
US4837517A (en)  Spatial frequency method and apparatus for investigating earth conductivity with high vertical resolution by induction techniques  
US6819111B2 (en)  Method of determining vertical and horizontal resistivity, and relative dip in anisotropic earth formations having an arbitrary electromagnetic antenna combination and orientation with additional rotation and position measurements  
US20060043973A1 (en)  Colocated antennas 
Legal Events
Date  Code  Title  Description 

EEER  Examination request  
MKLA  Lapsed 
Effective date: 20160815 