US7043370B2 - Real time processing of multicomponent induction tool data in highly deviated and horizontal wells - Google Patents

Real time processing of multicomponent induction tool data in highly deviated and horizontal wells Download PDF

Info

Publication number
US7043370B2
US7043370B2 US10/652,762 US65276203A US7043370B2 US 7043370 B2 US7043370 B2 US 7043370B2 US 65276203 A US65276203 A US 65276203A US 7043370 B2 US7043370 B2 US 7043370B2
Authority
US
United States
Prior art keywords
horizontal
resistivity
measurements
obtaining
processor
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 - Lifetime, expires
Application number
US10/652,762
Other versions
US20050049792A1 (en
Inventor
Liming Yu
Berthold Kriegshäuser
Otto Fanini
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Baker Hughes Holdings LLC
Original Assignee
Baker Hughes Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Baker Hughes Inc filed Critical Baker Hughes Inc
Priority to US10/652,762 priority Critical patent/US7043370B2/en
Assigned to BAKER HUGHES INCORPORATED reassignment BAKER HUGHES INCORPORATED ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: YU, LIMING, FANINI, OTTO, KRIEGSHAUSER, BERTHOLD
Priority to PCT/US2004/027751 priority patent/WO2005024467A1/en
Publication of US20050049792A1 publication Critical patent/US20050049792A1/en
Application granted granted Critical
Publication of US7043370B2 publication Critical patent/US7043370B2/en
Adjusted expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

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

Definitions

  • the invention is related generally to the field of electrical resistivity well logging methods. More specifically, the invention is related to methods for determining resistivity values in a multi-component induction measurement in highly-deviated wells.
  • Electromagnetic induction and wave propagation logging tools are commonly used for the 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, reasonably determine the petrophysical properties of the formation and the fluids therein.
  • a typical electrical induction instrument is described in U.S. Pat. No. 5,452,761 (Beard.)
  • the induction logging instrument described in Beard includes a number of receiver coils spaced at various axial distances from a transmitter coil. Alternating current is passed through the transmitter coil, which induces alternating electromagnetic fields in the earth formations. Voltages, , are induced in the receiver coils as a result of electromagnetic induction phenomena related to the alternating electromagnetic fields. A continuous record of the voltages form curves. These curves are also referred to as induction logs.
  • Induction instruments that are composed of multiple sets of receiver coils are referred to as multi-array induction instruments. Every set of receiver coils together with the transmitter is considered as a subarray.
  • a multi-array induction tool comprises numerous subarrays and acquires measurements with all the subarrays.
  • Voltages induced in axially more distal receiver coils are the result of electromagnetic induction phenomena occurring in a larger volume surrounding the instrument, and voltages induced in axially proximal receiver coils are the result of induction phenomena occurring more proximal to the instrument. Therefore, different receiver coils see a formation layer boundary with different shoulder-bed contributions, or shoulder-bed effects. The longer-spaced receiver coils see the formation layer boundary at further distance from the boundary than the shorter-spaced receiver coils do. As a result, the logs of longer-spaced receiver coils have longer shoulder-bed effects than the logs of shorter-spaced receiver coils. The logs of all the receiver coils form a certain pattern.
  • 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 non-conductive 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 non-conductive 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 non-conductive layers, layers which are typically hydrocarbon-bearing, are of the most interest to the user. Some earth formations which might be of commercial interest therefore may be overlooked by interpreting a well log made using conventional electromagnetic induction resistivity well logging instruments.
  • the thin layers may correspond to geologically distinct intervals, yet if the laminations are beyond the resolving capability of the logging instrument, it would appear as a single layer (for geophysical purposes).
  • the term “layer” could include both geological layers and geophysical layers.
  • 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 which is parallel to the axis of the instrument is calculated, from the rotated magnitudes, with respect to the 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 box of rotation.
  • An anisotropy parameter is calculated from the receiver signal magnitudes after the second box of rotation.
  • Vertical conductivity is calculated from the horizontal conductivity and the anisotropy parameter.
  • the method of U.S. Pat. No. 5,867,806 to Strickland, et al. corrects for shoulder bed effect in LWD resistivity logs through inversion.
  • the method selects one or more control depths at one or more locations of each of a plurality of detected beds in the formation.
  • the method estimates the resistivity of each bed only at the selected control depths to produce an estimated resistivity of the beds.
  • the method then computes a simulated log value at each control depth using a current estimate of the resistivity of the beds.
  • the computer-simulated log is then computed to the actual log data at each control depth, and the resistivity of each bed is adjusted using the difference between the actual and simulated values at the control depths.
  • the method iteratively repeats a plurality of times until the simulated log substantially matches the actual log at the control depths.
  • the relative formation dip angle is vital for proper and accurate interpretation of data acquired by the new multi-component induction instrument.
  • This newly developed induction instrument comprises three mutually orthogonal transmitter-receiver arrays. These configurations allow us to determine both horizontal and vertical resistivities for an anisotropic formation in vertical, deviated, and horizontal boreholes.
  • a description of the tool can be found in U.S. Pat. No. 6,147,496, to Strack, et al.
  • the transmitters induce currents in all three spatial directions and the receivers measure the corresponding magnetic fields (H xx , H yy , and H zz ,).
  • the first index indicates the direction of the transmitter and the second index denotes the receiver direction.
  • H zz is the magnetic field induced by a z-direction transmitter coil and measured by a z-directed receiver.
  • the z-direction is parallel to the borehole.
  • the instrument measures all other cross-components of the magnetic fields, i.e., H xy , H xz , H yx , H yz , H zx , and H zy .
  • the signals acquired by the main receiver coils (H xx , H yy , and H zz ) are used to determine both the horizontal and vertical resistivity of the formation. This is done by inverse processing techniques of the data. These inverse processing techniques automatically adjust formation parameters in order to optimize in a least-square sense the data match of the synthetic tool responses with measured data. Required inputs in this process are accurate information of the relative formation dip and relative formation azimuth. This information can be derived using, in addition to the main signals (H xx , H yy , and H zz ), the data from the cross-components. In highly deviated wells (70° or higher), skin-depth-corrected single frequency H xx +H yy measurements depend largely on the horizontal resistivity of the formation. The sum of dual frequency H xx +H yy measurements also mainly depends on the horizontal resistivity of the formation.
  • Prior methods are useful in determining resistivity values in wells where the angle of deviation from vertical is substantially less than 70 degrees. Typically, errors in conductivity values appear in highly-deviated wells. There is a need for a method that obtains a resistivity value in highly-deviated and horizontal wells.
  • the present invention fulfills that need.
  • the present invention is a method of logging a subsurface formation comprising a plurality of layers each having a horizontal conductivity and a vertical conductivity.
  • the method is primarily applicable to near horizontal boreholes (i.e., where the borehole axis is within 30 degrees of bedding.
  • Principal component resistivity measurements are obtained (or obtained by coordinate transformation of other suitable measurements).
  • the zz-measurements are used for deriving an initial model of horizontal and vertical resisitivities. An anisotropy factor may be assumed for the purposes.
  • an updated value of the horizontal resistivity is obtained.
  • an updated value of the vertical resistivity is obtained. The process may be iterated until convergence is reached.
  • Skin depth correction may be applied to the measurements. Multiple spacing measurements may be used, as could multifrequency measurements. These improve the convergence.
  • the iterative process includes only one variable at a time, greatly speeding up the computation.
  • FIG. 1 (Prior Art) shows a measurement-while-drilling tool suitable for use with the present invention.
  • FIG. 2 (Prior Art) depicts a multi-layer model of an anisotropic resistivity formation.
  • FIG. 3 (Prior Art) shows a horizontal well.
  • FIG. 4 shows a two-coil configuration of a transmission coil and a receiver coil.
  • FIG. 5 shows a three-coil configuration of a transmission coil and two receiver coils.
  • FIG. 6 shows a multi-component configuration for well exploration.
  • FIG. 7 shows a flowchart outlining a method of the invention.
  • the method of the present invention can be used in, among others, wireline and measurement-while-drilling (MWD) devices.
  • a typical wireline electromagnetic device is shown generally, for example, in U.S. Pat. No. 6,348,792 B1, issued to Beard, et al.
  • An MWD device suitable for use with the present invention is described generally below. Items of particular importance to wireline and MWD devices will be pointed out where appropriate.
  • FIG. 1 shows a schematic diagram of a drilling system 10 with a drillstring 20 carrying a drilling assembly 90 (also referred to as the bottom hole assembly, or “BHA”) conveyed in a “wellbore” or “borehole” 26 for drilling the wellbore.
  • the drilling system 10 includes a conventional derrick 11 erected on a floor 12 which supports a rotary table 14 that is rotated by a prime mover such as an electric motor (not shown) at a desired rotational speed.
  • the drillstring 20 includes a tubing such as a drill pipe 22 or a coiled-tubing extending downward from the surface into the borehole 26 . The drillstring 20 is pushed into the wellbore 26 when a drill pipe 22 is used as the tubing.
  • a tubing injector such as an injector (not shown), however, is used to move the tubing from a source thereof, such as a reel (not shown), to the wellbore 26 .
  • the drill bit 50 attached to the end of the drillstring breaks up the geological formations when it is rotated to drill the borehole 26 .
  • the drillstring 20 is coupled to a drawworks 30 via a Kelly joint 21 , swivel 28 , and line 29 through a pulley 23 .
  • the drawworks 30 is operated to control the weight on bit, which is an important parameter that affects the rate of penetration.
  • the operation of the drawworks is well known in the art and is thus not described in detail herein.
  • a suitable drilling fluid 31 from a mud pit (source) 32 is circulated under pressure through a channel in the drillstring 20 by a mud pump 34 .
  • the drilling fluid passes from the mud pump 34 into the drillstring 20 via a desurger (not shown), fluid line 28 and Kelly joint 21 .
  • the drilling fluid 31 is discharged at the borehole bottom 51 through an opening in the drill bit 50 .
  • the drilling fluid 31 circulates uphole through the annular space 27 between the drillstring 20 and the borehole 26 and returns to the mud pit 32 via a return line 35 .
  • the drilling fluid acts to lubricate the drill bit 50 and to carry borehole cutting or chips away from the drill bit 50 .
  • a sensor S 1 preferably placed in the line 38 provides information about the fluid flow rate.
  • a surface torque sensor S 2 and a sensor S 3 associated with the drillstring 20 respectively provide information about the torque and rotational speed of the drillstring.
  • a sensor (not shown) associated with line 29 is used to provide the hook load of the drillstring 20
  • the drill bit 50 is rotated by only rotating the drill pipe 22 .
  • a downhole motor 55 (mud motor) is disposed in the drilling assembly 90 to rotate the drill bit 50 and the drill pipe 22 is rotated usually to supplement the rotational power, if required, and to effect changes in the drilling direction.
  • the mud motor 55 is coupled to the drill bit 50 via a drive shaft (not shown) disposed in a bearing assembly 57 .
  • the mud motor rotates the drill bit 50 when the drilling fluid 31 passes through the mud motor 55 under pressure.
  • the bearing assembly 57 supports the radial and axial forces of the drill bit.
  • a stabilizer 58 coupled to the bearing assembly 57 acts as a centralizer for the lowermost portion of the mud motor assembly.
  • a drilling sensor module 59 is placed near the drill bit 50 .
  • the drilling sensor module contains sensors, circuitry, and processing software and algorithms relating to the dynamic drilling parameters. Such parameters preferably include bit bounce, stick-slip of the drilling assembly, backward rotation, torque, shocks, borehole and annulus pressure, acceleration measurements and other measurements of the drill bit condition.
  • a suitable telemetry or communication sub 72 using, for example, two-way telemetry, is also provided as illustrated in the drilling assembly 90 .
  • the drilling sensor module processes the sensor information and transmits it to the surface control unit 40 via the telemetry system 72 .
  • the communication sub 72 , a power unit 78 and an MWD tool 79 are all connected in tandem with the drillstring 20 . Flex subs, for example, are used in connecting the MWD tool 79 in the drilling assembly 90 . Such subs and tools form the bottom hole drilling assembly 90 between the drillstring 20 and the drill bit 50 .
  • the drilling assembly 90 makes various measurements including the pulsed nuclear magnetic resonance measurements while the borehole 26 is being drilled.
  • the communication sub 72 obtains the signals and measurements and transfers the signals, using two-way telemetry, for example, to be processed on the surface. Alternatively, the signals can be processed using a downhole processor in the drilling assembly 90 .
  • the surface control unit or processor 40 also receives signals from other downhole sensors and devices and signals from sensors S 1 –S 3 and other sensors used in the system 10 and processes such signals according to programmed instructions provided to the surface control unit 40 .
  • the surface control unit 40 displays desired drilling parameters and other information on a display/monitor 42 utilized by an operator to control the drilling operations.
  • the surface control unit 40 preferably includes a computer or a microprocessor-based processing system, memory for storing programs or models and data, a recorder for recording data, and other peripherals.
  • the control unit 40 is preferably adapted to activate alarms 44 when certain unsafe or undesirable operating conditions occur.
  • FIG. 2 shows an earth model usable with the present invention comprising multi-layers, labeled 201 a , 201 b , . . . , 201 n ⁇ 1 , 201 n .
  • Each layer is transversely isotropic, having horizontal and vertical resistivities.
  • the vertical resistivity is the resistivity for a current flowing along a line perpendicular to the bed layers.
  • the horizontal resistivity is the resistivity measurement for a current flowing in the plane of the bedding.
  • the vertical resistivity R v is always greater than or equal to the horizontal resistivity R h .
  • FIG. 3 a horizontal well 301 oriented substantially in the plane parallel with the horizontal resistivity R h is shown.
  • the angle of deviation ⁇ of the horizontal well is 90°.
  • FIG. 4 shows a two-coil configuration in a horizontal position of a measurement device.
  • the two-coil configuration comprises a transmitter 401 and a receiver 403 having a set separation r between them.
  • the transmitter and receiver are oriented substantially along the longitudinal axis of the device, thereby enabling measurement of the longitudinal component of the magnetic response.
  • the two-coil configuration could be directed perpendicular to the longitudinal axis, enabling measurement of x- and y-components of the magnetic response.
  • the z-axis is parallel to the tool axis and the borehole while the x- and y-axes are orthogonal to each other and to the z-axis.
  • FIG. 5 shows a typical three-coil configuration.
  • the three-coil configuration comprises an array having a transmitter T z 501 and two receivers R z1 503 and R z2 504 .
  • the array is configured such that the tool response in free space is zero.
  • Distance r 1 separates T z and R z1
  • distance r 2 separates T z and R z2 .
  • the transmitter and the two receivers are directed along the longitudinal axis, thereby enabling measurement of a longitudinal component (z-component) of the magnetic response.
  • This array can be directed perpendicular to the longitudinal axis, thereby enabling measurement of x- and y-components of the magnetic response.
  • FIG. 6 shows a typical multi-component induction device, comprising three mutually orthogonal transmitter-receiver coil arrays.
  • Transmitter T z 630 and receivers R z1 632 and R z2 633 are oriented along the longitudinal axis.
  • This coil array can be referred to as the ZZ array.
  • the ZZ array measures H zz components.
  • Similar arrays XX and YY are located axially along the drill tool.
  • Array XX comprises transmitter T x 610 and receivers R x1 612 and R x2 613 .
  • the coils of the XX array are directed substantially perpendicular to the axis of the drill tool.
  • Array YY comprises transmitter T y 620 and receivers R y1 622 and R y2 623 .
  • the coils of the YY array are directed substantially perpendicular to the axis of the drill tool, and substantially perpendicular to the direction of the XX array.
  • the Y-axis is defined to lie in a vertical plane through the tool axis. Measurements can be made for the conventional magnetic field component parallel to the borehole axis, H zz , and the two components perpendicular to the wellbore, H xx and H yy . It should be noted that the particular orientation of the coils (i.e., along the x-, y- and z-axes is a matter of convenience. Other orientations could also be used and the data recorded therein rotated to the components discussed below. Alternatively, the processing methodology could be used for other coil orientations with a more complicated analysis.
  • a multi-component tool in the manner of a two-coil configuration at a single frequency.
  • the low frequency approximations of the magnetic field components can be given by:
  • is the frequency of the applied field
  • r is the distance separating transmitter and receiver
  • R h is the horizontal conductivity of the formation layer
  • R v is the vertical conductivity of the formation layer.
  • M x , M y , and M z are the magnitudes of the magnetic dipole moments produced in the x-, y-, and z-directions, respectively, and IH xx , IH yy , and IH zz represents the imaginary component of the H xx , H yy , and H zz response.
  • Eqs. (1) and (2) demonstrate the characteristics found in whole space solutions. This shows that the H xx and H yy response depends only on R h .
  • Eqs. (3) and (4) measure an apparent conductivity corrected for skin depth.
  • IH xx , IH yy , and IH zz represents the imaginary component of the H xx , H yy , and H zz response.
  • the skin depth corrected apparent conductivity of the H xx and H yy responses depend on R h for a given spacing of r 1 and r 2 .
  • the skin depth corrected apparent conductivity of the zz array is related to the square root of the product of horizontal and vertical resistivities,
  • the responses of xx and yy arrays in a non-horizontal well are functions of borehole deviation ⁇ , tool rotation, and both horizontal and vertical resistivities.
  • the sum H xx +H yy is similarly dependent on borehole deviation and horizontal and vertical resistivities, but is independent of the tool rotation angle.
  • the response of the zz array is a function of the borehole deviation ⁇ and both horizontal and vertical resistivities, but is also independent of the angle of tool rotation.
  • equations (5) and (6) lose accuracy when the relative deviation ⁇ decreases from 90° (well deviates from horizontal).
  • the sum of H xx and H yy in Eq.(5) mainly depends on the horizontal resistivity R h at high deviation.
  • Transmitter and receivers spacings are given by r′ 1 and r′ 2 .
  • the relative error of the horizontal resistivity computed from equation (5) compared with the results using equation (3) is 17%.
  • the near zone effect on the H zz response is small.
  • the single frequency responses H xx and H yy are affected by the near zone environment, e.g. borehole mud and tool eccentricity.
  • the effect of a mud cake and an invaded zone can be significant.
  • the effect of tool eccentricity can be so large so as to make it impossible to use the single frequency xx and yy array data effectively.
  • a dual frequency response leads to the simplest method of processing multi-frequency data. This method is shown in equation (7):
  • H d H ⁇ ( f l ) - f l f h ⁇ H ⁇ ( f h ) ( 7 ) where f h is a high frequency of operation, f l is a low frequency of operation.
  • H xx and H yy are not simple functions of the horizontal resistivity. Rather, they are functions of tool rotation and both horizontal and vertical resistivities. However, the summation of xx and yy array dual frequency responses is independent of tool rotation, and the low frequency approximation depends mainly on the horizontal resistivity. For a three-coil configuration in horizontal wells, this summation of xx and yy arrays leads to the dual-frequency response:
  • FIG. 7 shows a flowchart of a preferred mode of the present invention. Analysis is performed on a level-by-level basis assuming a whole space model with transverse isotropy.
  • Box 701 comprises obtaining multi-component induction data. Box 701 can be performed using, for instance, a multi-component tool as shown in FIG. 6 , or a 3DEXTM induction tool of Baker Hughes Incorporated.
  • an initial value of R h and R v is chosen using the skin depth corrected apparent conductivity obtained from single frequency H zz data. As H zz is dependent on a product of the horizontal and vertical resistivities, an assumption is made of an isotropic medium for this calculation, thereby enabling an initial value of R h and R v .
  • an anisotropy factor may be assumed, giving an initial values for R h and R v . It is to be noted that the method of the present invention is quite robust, so that the method described above for determination of the initial model is not essential: almost any initial model can be used without detracting from the method of the invention.
  • a first inversion for R h is obtained preferably using a sum of dual frequency H xx and H yy values. Again, using the sum of H xx and H yy values is not essential to the invention.
  • the method of the present invention will also work using H zz and either one of H xx and H yy if the toolface orientation is known. It should also be pointed out that in a horizontal borehole with a sufficiently low frequency of operation of the tool, it is not even necessary to know the toolface orientation. In MWD tools, the tool will be rotating and an orientation sensor is necessary for defining the times at which the x- and y-axes are in a horizontal and vertical plane respectively.
  • An orientation sensor may also be used in wireline devices, though in deviated boreholes, the tool will be in a generally fixed orientation due to gravity.
  • R v is held fixed.
  • the horizontal resistivity can be inverted using single frequency values of H xx , H yy , or the dual frequency value of H xx +H yy . Because the sum of dual frequency H xx and H yy has much less near zone effects that the single frequency H xx and H yy measurements, it is preferred in the procedure; whereas, single frequency measurements will require more iterations to produce accurate results.
  • a second inversion is performed in Box 704 to obtain a value of R v .
  • the second inversion of Box 704 utilizes the single frequency values of H zz with the value of R h held fixed.
  • the orientation sensors may be selected from any of the commonly used sensors such as (i) a magnetometer, (ii) a gyroscope, (iii) an accelerometer, and, (iv) an inclinometer. It should also be pointed out that in a deviated borehole, a wireline logging tool will commonly orient itself in a preferred direction just from gravity alone, and hence an orientation sensor may not be necessary provided the xx- and yy-directions are known relative to the preferred direction for gravity orientation.
  • the inversions can be performed using a full solution synthetic model to obtain a full solution of the whole space model. Because only one unknown parameter is used in each inversion, an explicit formula can be written. Each subsequent inversion uses an initial R h determined using an iterative process, e.g:
  • R h m + 1 R h m + ⁇ ⁇ ⁇ R h m ( 9 )
  • R h m is the value of the horizontal resistivity at the m th iteration
  • ⁇ R h m is the horizontal resistivity correction at the m th iteration.
  • ⁇ R h m continually updates horizontal resistivity. In most cases, two iterations of the inversion process provide an acceptable degree of accuracy.
  • the advantages of the method of FIG. 7 are found in its extreme speed and the insensitivity of the accuracy to the angle of relative deviation.
  • the method of the present invention typically can process 1000 feet of data level by level on a Sun Microsystems ULTRA 5 computer in about 20 seconds.
  • a typical logging speed is 60 ft/minute, so it can be seen that the method of the present invention is substantially faster than real-time.
  • the whole-space inversion can be implemented as real-time processing software for 3DEXTM data in highly deviated and horizontal wells. This processing is preferably done using a downhole proessor, though any combination of surface and downhole processing may be used.
  • This processing is applicable for geosteering as the drill bit path is close to the delineating plane of the target-drilling horizon.
  • This processing can be combined with different frequency focusing techniques like dual frequency or multi-frequency processing.
  • An example of such processing can be found, for example in U.S. Pat. No. 6,147,496, issued to Strack ,et al which teaches a dual frequency method and apparatus.
  • the invention also can use frequency-focusing terms like the k 3 term of a wave Taylor series expansion, which has been shown to have little near zone effects.
  • An example of such processing can be found, for example in U.S. Pat. No. 6,574,562, issued to Tabarovsky, et al. having the same assignee and the contents of which are incorporated herein by reference.
  • frequency focusing as used in this invention is intended to include both the dual frequency method described by Strack et al, as well as the more general method taught by Tabarovsky et al. It is to be further noted that the method described in Tabarovsky is not applicable to the case where measurements are made with an electromagnetic logging tool that has a finite and non-zero conductivity.
  • U.S. patent application Ser. No. 10/295,969 of Tabarovsky having the same assignee as the present application and the contents of which are fully incorporated herein by reference, teaches a frequency focusing method and apparatus that is applicable to the case where the logging tool has a casing with a finite non-zero conductivity.
  • the induction sonde can also have multiple spacings for depth of investigation focusing and the coils can have an arbitrary orientation to capture the essential components of the electromagnetic signal in the formation.
  • a method of rotating to appropriate coordinate systems is discussed, for example, in U.S. patent application Ser. No. 10/302,239, filed on Nov. 22, 2002.
  • this data can be transformed into equivalent orthogonal coil measurements to be processed. Care should be taken to distinguish the shoulder bed effect. Error due to this effect typically occurs in the resistive layers that are embedded within the conductive layers above and below.
  • the present invention is applicable only when the inclination angle is large: typically, the borehole axis is within 30° of the bedding plane.
  • the method of the present invention is computationally fast, and can be used with only the principal component measurements (i.e., the xx, yy, and zz components). The analysis described above made use of these principal components. As noted above, any set of components could be used as long as the three principal components can be recovered using standard coordinate rotation methods.
  • Synthetic data was generated using both a whole-space model and a layered model.
  • Field data was obtained using a multicomponent induction tool in a customer well.
  • Table 1 shows the results of the applied method using the whole space model at a given deviation angle (80°).
  • the horizontal resistivity is 1 ⁇ -m and the vertical resistivity is varied over a range (1 ⁇ -m to 100 ⁇ -m).
  • the anisotropy ratio which is R v /R h , therefore varies from 1 to 100.
  • Table 2 shows results with the same set of horizontal and vertical resistivities taken at a different relative deviation angle of 70°.
  • Table 3 shows results at a deviation angle of 65°.
  • the relative errors of inverted horizontal and vertical resistivities are less than or approximately 1%.
  • the relative errors of inverted horizontal and vertical resistivities are less than 2%.
  • the relative errors increase with an increase in the anisotropy ratio. In most real situations, the anisotropy ratio will be less than 10. Similar results can be obtained through an increase in the horizontal resistivity.
  • the anisotropy ratio is a determining factor.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Electromagnetism (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

A method is discussed for obtaining resistivity values in a horizontal or highly-deviated wellbore penetrating an anisotropic formation layer. Multi-component data is obtained from the formation layer. An initial value of horizontal and vertical resistivity can be determined by using a whole space model and assuming an isotropic layer. Revised values of the vertical and horizontal resistivities are obtained through a first and second inversion process. A first inversion is performed on the horizontal resistivity, holding a value of vertical resistivity constant. A second inversion is performed on the vertical resistivity, holding a value of horizontal resistivity constant. These inversions are iterated until a desired convergence is achieved.

Description

BACKGROUND OF THE INVENTION
1. Field of the Invention
The invention is related generally to the field of electrical resistivity well logging methods. More specifically, the invention is related to methods for determining resistivity values in a multi-component induction measurement in highly-deviated wells.
2. Description of the Related Art
Electromagnetic induction and wave propagation logging tools are commonly used for the 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, reasonably determine the petrophysical properties of the formation and the fluids therein.
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 Oil-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,. Examples of such modifications and improvements can be found, for example, in U.S. Pat. No. 4,837,517; to Barber, 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.
A typical electrical induction instrument is described in U.S. Pat. No. 5,452,761 (Beard.) The induction logging instrument described in Beard includes a number of receiver coils spaced at various axial distances from a transmitter coil. Alternating current is passed through the transmitter coil, which induces alternating electromagnetic fields in the earth formations. Voltages, , are induced in the receiver coils as a result of electromagnetic induction phenomena related to the alternating electromagnetic fields. A continuous record of the voltages form curves. These curves are also referred to as induction logs. Induction instruments that are composed of multiple sets of receiver coils are referred to as multi-array induction instruments. Every set of receiver coils together with the transmitter is considered as a subarray. Hence, a multi-array induction tool comprises numerous subarrays and acquires measurements with all the subarrays.
Voltages induced in axially more distal receiver coils are the result of electromagnetic induction phenomena occurring in a larger volume surrounding the instrument, and voltages induced in axially proximal receiver coils are the result of induction phenomena occurring more proximal to the instrument. Therefore, different receiver coils see a formation layer boundary with different shoulder-bed contributions, or shoulder-bed effects. The longer-spaced receiver coils see the formation layer boundary at further distance from the boundary than the shorter-spaced receiver coils do. As a result, the logs of longer-spaced receiver coils have longer shoulder-bed effects than the logs of shorter-spaced receiver coils. The logs of all the receiver coils form a certain pattern.
Variations of measurements with azimuthal angle can be corrected. If the layers are not perpendicular to the axis of the instrument, the conductivity of the media surrounding the instrument can vary with azimuthal angle, causing any inferences about the conductivity from the measurements of the induction voltage to be in error. A method for correcting this error is described in U.S. Pat. No. 5,774,360 issued to Xiao and Zhou. The method requires the relative dip angle as a priori information. The relative dip angle is the angle between the borehole axis and the normal of the bedding plane. Because the formation layers can also be inclined, the relative dip angle is normally unknown even though the wellbore deviation is known. U.S. Pat. No. 6,049,209 issued to Xiao and Geldmacher teaches another method that has also been developed to interpret induction logs in environments of relative inclination and anisotropy. The method requires the relative dip angle and the anisotropy coefficient as a priori information. The anisotropy coefficient can be defined as the ratio between the resistivity perpendicular to bedding and the resistivity parallel to bedding.
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 non-conductive 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 non-conductive 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 non-conductive layers, layers which are typically hydrocarbon-bearing, are of the most interest to the user. Some earth formations which might be of commercial interest therefore may be overlooked by interpreting a well log made using conventional electromagnetic induction resistivity well logging instruments. The thin layers may correspond to geologically distinct intervals, yet if the laminations are beyond the resolving capability of the logging instrument, it would appear as a single layer (for geophysical purposes). Hence for the purposes of this invention, the term “layer” could include both geological layers and geophysical layers.
U.S. Pat. No. 5,999,883 issued to Gupta, et al., (the “Gupta patent”), the contents of which are fully incorporated herein by reference, discloses a method for determining the horizontal and vertical conductivity of anisotropic earth formations. Measurements are made of electromagnetic signals induced by induction transmitters oriented along three mutually orthogonal axes. 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 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 which is parallel to the axis of the instrument is calculated, from the rotated magnitudes, with respect to the 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 box of rotation. An anisotropy parameter is calculated from the receiver signal magnitudes after the second box of rotation. Vertical conductivity is calculated from the horizontal conductivity and the anisotropy parameter.
Shoulder bed corrections related to the effect of formations above and below the depth being evaluated also have to be applied to the data. Methods for making these corrections to data acquired with conventional logging tools are well known in the art. For example, U.S. Pat. No. 5,446,654 to Chemali teaches the conversion of a resistivity log as a function of well depth into a rectangularized curve so that the interfaces of the adjacent strata are located, and a suitable number of iterations, a correction factor is applied. The corrected rectangular log is obtained with a correction coefficient computed at each depth. For each computation, the impact of all the strata within a specified depth window is considered, while strata beyond that window are simplified by representing the strata beyond the window with single equivalent bed values to reduce the number of computations required. This then provides a resistivity log which is substantially free of shoulder bed effect.
The method of U.S. Pat. No. 5,867,806 to Strickland, et al. corrects for shoulder bed effect in LWD resistivity logs through inversion. The method selects one or more control depths at one or more locations of each of a plurality of detected beds in the formation. The method then estimates the resistivity of each bed only at the selected control depths to produce an estimated resistivity of the beds. The method then computes a simulated log value at each control depth using a current estimate of the resistivity of the beds. The computer-simulated log is then computed to the actual log data at each control depth, and the resistivity of each bed is adjusted using the difference between the actual and simulated values at the control depths. The method iteratively repeats a plurality of times until the simulated log substantially matches the actual log at the control depths.
Electrically anisotropic reservoirs are encountered frequently in hydrocarbon exploration. For accurate saturation estimation and optimum hydrocarbon recovery from these reservoirs, it is essential to detect and properly describe their electrical properties. For example, in laminated sand-shale sequences or sands with different grain size distributions, the vertical resistivity (perpendicular to the bedding) is more indicative of the hydrocarbon content than the horizontal resistivity (parallel to the bedding). However, the response measured by conventional induction tools with their transmitter-receiver coil moments oriented normal to bedding is dominated by the horizontal resistivity. Therefore, a petrophysical evaluation based on these data can either overlook hydrocarbons present in laminated sands or underestimate their productivity.
The relative formation dip angle is vital for proper and accurate interpretation of data acquired by the new multi-component induction instrument. This newly developed induction instrument comprises three mutually orthogonal transmitter-receiver arrays. These configurations allow us to determine both horizontal and vertical resistivities for an anisotropic formation in vertical, deviated, and horizontal boreholes. A description of the tool can be found in U.S. Pat. No. 6,147,496, to Strack, et al. The transmitters induce currents in all three spatial directions and the receivers measure the corresponding magnetic fields (Hxx, Hyy, and Hzz,). In this nomenclature of the field responses, the first index indicates the direction of the transmitter and the second index denotes the receiver direction. As an example, Hzz is the magnetic field induced by a z-direction transmitter coil and measured by a z-directed receiver. The z-direction is parallel to the borehole. In addition, the instrument measures all other cross-components of the magnetic fields, i.e., Hxy, Hxz, Hyx, Hyz, Hzx, and Hzy.
The signals acquired by the main receiver coils (Hxx, Hyy, and Hzz) are used to determine both the horizontal and vertical resistivity of the formation. This is done by inverse processing techniques of the data. These inverse processing techniques automatically adjust formation parameters in order to optimize in a least-square sense the data match of the synthetic tool responses with measured data. Required inputs in this process are accurate information of the relative formation dip and relative formation azimuth. This information can be derived using, in addition to the main signals (Hxx, Hyy, and Hzz), the data from the cross-components. In highly deviated wells (70° or higher), skin-depth-corrected single frequency Hxx+Hyy measurements depend largely on the horizontal resistivity of the formation. The sum of dual frequency Hxx+Hyy measurements also mainly depends on the horizontal resistivity of the formation.
Prior methods are useful in determining resistivity values in wells where the angle of deviation from vertical is substantially less than 70 degrees. Typically, errors in conductivity values appear in highly-deviated wells. There is a need for a method that obtains a resistivity value in highly-deviated and horizontal wells. The present invention fulfills that need.
SUMMARY OF THE INVENTION
The present invention is a method of logging a subsurface formation comprising a plurality of layers each having a horizontal conductivity and a vertical conductivity. The method is primarily applicable to near horizontal boreholes (i.e., where the borehole axis is within 30 degrees of bedding. Principal component resistivity measurements are obtained (or obtained by coordinate transformation of other suitable measurements). The zz-measurements are used for deriving an initial model of horizontal and vertical resisitivities. An anisotropy factor may be assumed for the purposes. Using the initial model and other principal components (xx- and yy-measurements), an updated value of the horizontal resistivity is obtained. Using this updated model, and the zz-measurements, an updated value of the vertical resistivity is obtained. The process may be iterated until convergence is reached.
Skin depth correction may be applied to the measurements. Multiple spacing measurements may be used, as could multifrequency measurements. These improve the convergence. In a preferred embodiment of the invention, the iterative process includes only one variable at a time, greatly speeding up the computation.
BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 (Prior Art) shows a measurement-while-drilling tool suitable for use with the present invention.
FIG. 2 (Prior Art) depicts a multi-layer model of an anisotropic resistivity formation.
FIG. 3 (Prior Art) shows a horizontal well.
FIG. 4 shows a two-coil configuration of a transmission coil and a receiver coil.
FIG. 5 shows a three-coil configuration of a transmission coil and two receiver coils.
FIG. 6 shows a multi-component configuration for well exploration.
FIG. 7 shows a flowchart outlining a method of the invention.
DESCRIPTION OF PREFERRED EMBODIMENT
The method of the present invention can be used in, among others, wireline and measurement-while-drilling (MWD) devices. A typical wireline electromagnetic device is shown generally, for example, in U.S. Pat. No. 6,348,792 B1, issued to Beard, et al. An MWD device suitable for use with the present invention is described generally below. Items of particular importance to wireline and MWD devices will be pointed out where appropriate.
FIG. 1 shows a schematic diagram of a drilling system 10 with a drillstring 20 carrying a drilling assembly 90 (also referred to as the bottom hole assembly, or “BHA”) conveyed in a “wellbore” or “borehole” 26 for drilling the wellbore. The drilling system 10 includes a conventional derrick 11 erected on a floor 12 which supports a rotary table 14 that is rotated by a prime mover such as an electric motor (not shown) at a desired rotational speed. The drillstring 20 includes a tubing such as a drill pipe 22 or a coiled-tubing extending downward from the surface into the borehole 26. The drillstring 20 is pushed into the wellbore 26 when a drill pipe 22 is used as the tubing. For coiled-tubing applications, a tubing injector, such as an injector (not shown), however, is used to move the tubing from a source thereof, such as a reel (not shown), to the wellbore 26. The drill bit 50 attached to the end of the drillstring breaks up the geological formations when it is rotated to drill the borehole 26. If a drill pipe 22 is used, the drillstring 20 is coupled to a drawworks 30 via a Kelly joint 21, swivel 28, and line 29 through a pulley 23. During drilling operations, the drawworks 30 is operated to control the weight on bit, which is an important parameter that affects the rate of penetration. The operation of the drawworks is well known in the art and is thus not described in detail herein.
During drilling operations, a suitable drilling fluid 31 from a mud pit (source) 32 is circulated under pressure through a channel in the drillstring 20 by a mud pump 34. The drilling fluid passes from the mud pump 34 into the drillstring 20 via a desurger (not shown), fluid line 28 and Kelly joint 21. The drilling fluid 31 is discharged at the borehole bottom 51 through an opening in the drill bit 50. The drilling fluid 31 circulates uphole through the annular space 27 between the drillstring 20 and the borehole 26 and returns to the mud pit 32 via a return line 35. The drilling fluid acts to lubricate the drill bit 50 and to carry borehole cutting or chips away from the drill bit 50. A sensor S1 preferably placed in the line 38 provides information about the fluid flow rate. A surface torque sensor S2 and a sensor S3 associated with the drillstring 20 respectively provide information about the torque and rotational speed of the drillstring. Additionally, a sensor (not shown) associated with line 29 is used to provide the hook load of the drillstring 20.
In one embodiment of the invention, the drill bit 50 is rotated by only rotating the drill pipe 22. In another embodiment of the invention, a downhole motor 55 (mud motor) is disposed in the drilling assembly 90 to rotate the drill bit 50 and the drill pipe 22 is rotated usually to supplement the rotational power, if required, and to effect changes in the drilling direction.
In the preferred embodiment of FIG. 1, the mud motor 55 is coupled to the drill bit 50 via a drive shaft (not shown) disposed in a bearing assembly 57. The mud motor rotates the drill bit 50 when the drilling fluid 31 passes through the mud motor 55 under pressure. The bearing assembly 57 supports the radial and axial forces of the drill bit. A stabilizer 58 coupled to the bearing assembly 57 acts as a centralizer for the lowermost portion of the mud motor assembly.
In one embodiment of the invention, a drilling sensor module 59 is placed near the drill bit 50. The drilling sensor module contains sensors, circuitry, and processing software and algorithms relating to the dynamic drilling parameters. Such parameters preferably include bit bounce, stick-slip of the drilling assembly, backward rotation, torque, shocks, borehole and annulus pressure, acceleration measurements and other measurements of the drill bit condition. A suitable telemetry or communication sub 72 using, for example, two-way telemetry, is also provided as illustrated in the drilling assembly 90. The drilling sensor module processes the sensor information and transmits it to the surface control unit 40 via the telemetry system 72.
The communication sub 72, a power unit 78 and an MWD tool 79 are all connected in tandem with the drillstring 20. Flex subs, for example, are used in connecting the MWD tool 79 in the drilling assembly 90. Such subs and tools form the bottom hole drilling assembly 90 between the drillstring 20 and the drill bit 50. The drilling assembly 90 makes various measurements including the pulsed nuclear magnetic resonance measurements while the borehole 26 is being drilled. The communication sub 72 obtains the signals and measurements and transfers the signals, using two-way telemetry, for example, to be processed on the surface. Alternatively, the signals can be processed using a downhole processor in the drilling assembly 90.
The surface control unit or processor 40 also receives signals from other downhole sensors and devices and signals from sensors S1–S3 and other sensors used in the system 10 and processes such signals according to programmed instructions provided to the surface control unit 40. The surface control unit 40 displays desired drilling parameters and other information on a display/monitor 42 utilized by an operator to control the drilling operations. The surface control unit 40 preferably includes a computer or a microprocessor-based processing system, memory for storing programs or models and data, a recorder for recording data, and other peripherals. The control unit 40 is preferably adapted to activate alarms 44 when certain unsafe or undesirable operating conditions occur.
FIG. 2 shows an earth model usable with the present invention comprising multi-layers, labeled 201 a, 201 b, . . . , 201 n1, 201 n. Each layer is transversely isotropic, having horizontal and vertical resistivities. The vertical resistivity is the resistivity for a current flowing along a line perpendicular to the bed layers. The horizontal resistivity is the resistivity measurement for a current flowing in the plane of the bedding. The vertical resistivity Rv is always greater than or equal to the horizontal resistivity Rh.
Whole-space analytical solutions can be found in various prior art (see, for example, Zhdanov, et al., 2001. “Foundations of tensor induction well-logging” Petrophysics, Vol. 42, No. 6, 588–610.). The low frequency response approximation can be treated as a skin-depth-corrected response. As shown below, in horizontal wells, the normalized imaginary components of magnetic field Hxx, Hyy and Hzz data are independent of tool rotation The Hxx and Hyy values depend only on the horizontal resistivity Rh, and the Hzz value depends on the square root of the product of horizontal and vertical resistivities √{square root over (RhRv)}.
Turning now to FIG. 3, a horizontal well 301 oriented substantially in the plane parallel with the horizontal resistivity Rh is shown. The angle of deviation θ of the horizontal well is 90°. Much of the discussion of the present invention is based on this geometry.
FIG. 4 shows a two-coil configuration in a horizontal position of a measurement device. The two-coil configuration comprises a transmitter 401 and a receiver 403 having a set separation r between them. As illustrated in the figure, the transmitter and receiver are oriented substantially along the longitudinal axis of the device, thereby enabling measurement of the longitudinal component of the magnetic response. Also, the two-coil configuration could be directed perpendicular to the longitudinal axis, enabling measurement of x- and y-components of the magnetic response. The z-axis is parallel to the tool axis and the borehole while the x- and y-axes are orthogonal to each other and to the z-axis.
FIG. 5 shows a typical three-coil configuration. The three-coil configuration comprises an array having a transmitter T z 501 and two receivers R z1 503 and R z2 504. The array is configured such that the tool response in free space is zero. Distance r1 separates Tz and Rz1, and distance r2 separates Tz and Rz2. In the figure, the transmitter and the two receivers are directed along the longitudinal axis, thereby enabling measurement of a longitudinal component (z-component) of the magnetic response. This array can be directed perpendicular to the longitudinal axis, thereby enabling measurement of x- and y-components of the magnetic response.
FIG. 6 shows a typical multi-component induction device, comprising three mutually orthogonal transmitter-receiver coil arrays. Transmitter Tz 630 and receivers Rz1 632 and R z2 633 are oriented along the longitudinal axis. This coil array can be referred to as the ZZ array. The ZZ array measures Hzz components. Similar arrays XX and YY are located axially along the drill tool. Array XX comprises transmitter T x 610 and receivers Rx1 612 and Rx2 613. The coils of the XX array are directed substantially perpendicular to the axis of the drill tool. Array YY comprises transmitter T y 620 and receivers Ry1 622 and Ry2 623. The coils of the YY array are directed substantially perpendicular to the axis of the drill tool, and substantially perpendicular to the direction of the XX array. The Y-axis is defined to lie in a vertical plane through the tool axis. Measurements can be made for the conventional magnetic field component parallel to the borehole axis, Hzz, and the two components perpendicular to the wellbore, Hxx and Hyy. It should be noted that the particular orientation of the coils (i.e., along the x-, y- and z-axes is a matter of convenience. Other orientations could also be used and the data recorded therein rotated to the components discussed below. Alternatively, the processing methodology could be used for other coil orientations with a more complicated analysis.
In a horizontal well, a multi-component tool can be used in the manner of a two-coil configuration at a single frequency. The low frequency approximations of the magnetic field components can be given by:
?? H xx M x = ?? H yy M y ω μ 0 8 π rR h ( 1 ) and ?? H zz M z = ω μ 0 4 π r R h × R v ( 2 )
where ω is the frequency of the applied field, r is the distance separating transmitter and receiver, Rh is the horizontal conductivity of the formation layer and Rv is the vertical conductivity of the formation layer. Mx, My, and Mz are the magnitudes of the magnetic dipole moments produced in the x-, y-, and z-directions, respectively, and ℑHxx, ℑHyy, and ℑHzz represents the imaginary component of the Hxx, Hyy, and Hzz response. Eqs. (1) and (2) demonstrate the characteristics found in whole space solutions. This shows that the Hxx and Hyy response depends only on Rh. These measurements can be taken using the multi-component induction tool of FIG. 6.
In a horizontal well, if a three-coil configuration is now employed at a single frequency, the low frequency approximations of the magnetic field responses obtained using this configuration are given by:
?? H xx M x = ?? H yy M y ω μ 0 8 π r 1 R h [ 1 - ( r 2 r 1 ) 2 ] , ( 3 ) and ?? H zz M z = ω μ 0 4 π r 1 R h × R v [ 1 - ( r 2 r 1 ) 2 ] . ( 4 )
where r1 is the distance separating Tz and Rz1, and r2 is the distance separating Tz and Rz2. Eqs. (3) and (4) measure an apparent conductivity corrected for skin depth. ℑHxx, ℑHyy, and ℑHzz represents the imaginary component of the Hxx, Hyy, and Hzz response. From Eq. (3), the skin depth corrected apparent conductivity of the Hxx and Hyy responses depend on Rh for a given spacing of r1 and r2. From Eq. (4), the skin depth corrected apparent conductivity of the zz array is related to the square root of the product of horizontal and vertical resistivities,
When the well is not horizontal, these relations become more complex. The responses of xx and yy arrays in a non-horizontal well are functions of borehole deviation θ, tool rotation, and both horizontal and vertical resistivities. The sum Hxx+Hyy is similarly dependent on borehole deviation and horizontal and vertical resistivities, but is independent of the tool rotation angle. The response of the zz array is a function of the borehole deviation θ and both horizontal and vertical resistivities, but is also independent of the angle of tool rotation. Using the three-coil configuration, a summation of the Hxx and Hyy responses can be given as
?? H xx M x + ?? H yy M y = ω μ 0 4 π r 1 R h [ 1 - ( r 2 r 1 ) 2 ] × [ 1 - cos 2 θ ( λ 2 - 1 ) λ sin 2 θ + λ 2 cos 2 θ ] ( 5 ) and ?? H zz M z = ω μ 0 4 π r 1 R h × R v [ 1 - ( r 2 r 1 ) 2 ] × sin 2 θ + λ 2 cos 2 θ , ( 6 ) where λ 2 = R v R h .
As the contribution of anisotropy is determined by the cos θ function, equations (5) and (6) lose accuracy when the relative deviation θ decreases from 90° (well deviates from horizontal). The sum of Hxx and Hyy in Eq.(5) mainly depends on the horizontal resistivity Rh at high deviation. Transmitter and receivers spacings are given by r′1 and r′2. As an example, using a the relative deviation of 70 degrees and a horizontal and vertical resistivities of 1 Ω-m and 5 Ω-m respectively, the relative error of the horizontal resistivity computed from equation (5) compared with the results using equation (3) is 17%.
Typically, the near zone effect on the Hzz response is small. In contrast, the single frequency responses Hxx and Hyy are affected by the near zone environment, e.g. borehole mud and tool eccentricity. In addition, in wireline implementations, the effect of a mud cake and an invaded zone can be significant. In highly deviated wells having conductive mud therein, the effect of tool eccentricity can be so large so as to make it impossible to use the single frequency xx and yy array data effectively. On the other hand, studies show that use of dual frequency and multifrequency data can significantly reduce the near zone effect. A dual frequency response leads to the simplest method of processing multi-frequency data. This method is shown in equation (7):
H d = H ( f l ) - f l f h H ( f h ) ( 7 )
where fh is a high frequency of operation, fl is a low frequency of operation.
The dual frequency responses of Hxx and Hyy are not simple functions of the horizontal resistivity. Rather, they are functions of tool rotation and both horizontal and vertical resistivities. However, the summation of xx and yy array dual frequency responses is independent of tool rotation, and the low frequency approximation depends mainly on the horizontal resistivity. For a three-coil configuration in horizontal wells, this summation of xx and yy arrays leads to the dual-frequency response:
?? H xx d M x + ?? H yy d M y = ( 2 π f l μ 0 ) 3 / 2 3 π R h 3 / 2 × [ 1 - ( r 2 r 1 ) 3 ] × ( 1 - f h f l ) × [ 1 - 3 8 ( 1 - R h R v ) ] . ( 8 )
where fh is a high frequency of operation, and fl is a low frequency of operation. As shown in Eq. (8), the maximum contribution from anisotropy is 37.5%.
FIG. 7 shows a flowchart of a preferred mode of the present invention. Analysis is performed on a level-by-level basis assuming a whole space model with transverse isotropy. Box 701 comprises obtaining multi-component induction data. Box 701 can be performed using, for instance, a multi-component tool as shown in FIG. 6, or a 3DEX™ induction tool of Baker Hughes Incorporated. In Box 702, an initial value of Rh and Rv is chosen using the skin depth corrected apparent conductivity obtained from single frequency Hzz data. As Hzz is dependent on a product of the horizontal and vertical resistivities, an assumption is made of an isotropic medium for this calculation, thereby enabling an initial value of Rh and Rv. Alternatively, an anisotropy factor may be assumed, giving an initial values for Rh and Rv. It is to be noted that the method of the present invention is quite robust, so that the method described above for determination of the initial model is not essential: almost any initial model can be used without detracting from the method of the invention.
In Box 703, a first inversion for Rh is obtained preferably using a sum of dual frequency Hxx and Hyy values. Again, using the sum of Hxx and Hyy values is not essential to the invention. The method of the present invention will also work using Hzz and either one of Hxx and Hyy if the toolface orientation is known. It should also be pointed out that in a horizontal borehole with a sufficiently low frequency of operation of the tool, it is not even necessary to know the toolface orientation. In MWD tools, the tool will be rotating and an orientation sensor is necessary for defining the times at which the x- and y-axes are in a horizontal and vertical plane respectively. An orientation sensor may also be used in wireline devices, though in deviated boreholes, the tool will be in a generally fixed orientation due to gravity. During this first inversion, Rv is held fixed. The horizontal resistivity can be inverted using single frequency values of Hxx, Hyy, or the dual frequency value of Hxx+Hyy. Because the sum of dual frequency Hxx and Hyy has much less near zone effects that the single frequency Hxx and Hyy measurements, it is preferred in the procedure; whereas, single frequency measurements will require more iterations to produce accurate results. After the first inversion of Box 703 has been performed, a second inversion is performed in Box 704 to obtain a value of Rv. The second inversion of Box 704 utilizes the single frequency values of Hzz with the value of Rh held fixed.
The orientation sensors may be selected from any of the commonly used sensors such as (i) a magnetometer, (ii) a gyroscope, (iii) an accelerometer, and, (iv) an inclinometer. It should also be pointed out that in a deviated borehole, a wireline logging tool will commonly orient itself in a preferred direction just from gravity alone, and hence an orientation sensor may not be necessary provided the xx- and yy-directions are known relative to the preferred direction for gravity orientation.
The inversions can be performed using a full solution synthetic model to obtain a full solution of the whole space model. Because only one unknown parameter is used in each inversion, an explicit formula can be written. Each subsequent inversion uses an initial Rh determined using an iterative process, e.g:
R h m + 1 = R h m + Δ R h m ( 9 )
where Rh m is the value of the horizontal resistivity at the mth iteration and ΔRh m is the horizontal resistivity correction at the mth iteration. ΔRh m is calculated using:
ΔR h m=(J T J+α)−1(J T Δd)  (10)
where J is the Jacobian vector comprising derivatives which relate the change in the data to small changes in the horizontal resistivity, α is a damping factor and Δd is the difference between the measured data and the synthetic data obtained from whole space forward modeling. Using Eq. (9), ΔRh m continually updates horizontal resistivity. In most cases, two iterations of the inversion process provide an acceptable degree of accuracy.
In Box 705, a check is made to determine the convergence of the values of Rh and Rv obtained in Box 703 and 704 respectively. If convergence has been sufficiently determined, the method stops (Box 706). Otherwise, the method returns to Box 703 for iteration of the inversion processes of Box 703 and 704. Additional iteration enables improved accuracy of the results. Subsequent initial Rh values are determined using the procedure of Eq. (9).
The advantages of the method of FIG. 7 are found in its extreme speed and the insensitivity of the accuracy to the angle of relative deviation. The method of the present invention typically can process 1000 feet of data level by level on a Sun Microsystems ULTRA 5 computer in about 20 seconds. A typical logging speed is 60 ft/minute, so it can be seen that the method of the present invention is substantially faster than real-time. The whole-space inversion can be implemented as real-time processing software for 3DEX™ data in highly deviated and horizontal wells. This processing is preferably done using a downhole proessor, though any combination of surface and downhole processing may be used.
This processing is applicable for geosteering as the drill bit path is close to the delineating plane of the target-drilling horizon. This processing can be combined with different frequency focusing techniques like dual frequency or multi-frequency processing. An example of such processing can be found, for example in U.S. Pat. No. 6,147,496, issued to Strack ,et al which teaches a dual frequency method and apparatus. The invention also can use frequency-focusing terms like the k3 term of a wave Taylor series expansion, which has been shown to have little near zone effects. An example of such processing can be found, for example in U.S. Pat. No. 6,574,562, issued to Tabarovsky, et al. having the same assignee and the contents of which are incorporated herein by reference. The term “frequency focusing” as used in this invention is intended to include both the dual frequency method described by Strack et al, as well as the more general method taught by Tabarovsky et al. It is to be further noted that the method described in Tabarovsky is not applicable to the case where measurements are made with an electromagnetic logging tool that has a finite and non-zero conductivity. U.S. patent application Ser. No. 10/295,969 of Tabarovsky, having the same assignee as the present application and the contents of which are fully incorporated herein by reference, teaches a frequency focusing method and apparatus that is applicable to the case where the logging tool has a casing with a finite non-zero conductivity. This too is intended to be within the scope of the term “frequency focusing.” The induction sonde can also have multiple spacings for depth of investigation focusing and the coils can have an arbitrary orientation to capture the essential components of the electromagnetic signal in the formation. A method of rotating to appropriate coordinate systems is discussed, for example, in U.S. patent application Ser. No. 10/302,239, filed on Nov. 22, 2002. Optionally, this data can be transformed into equivalent orthogonal coil measurements to be processed. Care should be taken to distinguish the shoulder bed effect. Error due to this effect typically occurs in the resistive layers that are embedded within the conductive layers above and below.
U.S. patent application Ser. No. 10/072,173 (now U.S. Pat. No. 6,636,045) of Tabarovsky et al having the same assignee as the present invention teaches a method in which a a subset of the measurements is used to determine the horizontal resistivity of the earth formations. The determined horizontal resistivities are used in a model for inversion of another subset of the data to obtain the vertical formations resistivities. The method of the Tabarovsky '173 application may be used for any borehole inclination and requires a full-fledged inversion of several components of data, including the cross components wherein the transmitter and receiver have different orientations. In contrast, the present invention is applicable only when the inclination angle is large: typically, the borehole axis is within 30° of the bedding plane. As noted above, the method of the present invention is computationally fast, and can be used with only the principal component measurements (i.e., the xx, yy, and zz components). The analysis described above made use of these principal components. As noted above, any set of components could be used as long as the three principal components can be recovered using standard coordinate rotation methods.
The preferred algorithm of the present invention has been tested on synthetic data and on field data. Synthetic data was generated using both a whole-space model and a layered model. Field data was obtained using a multicomponent induction tool in a customer well.
Table 1 shows the results of the applied method using the whole space model at a given deviation angle (80°). The horizontal resistivity is 1 Ω-m and the vertical resistivity is varied over a range (1 Ω-m to 100 Ω-m). The anisotropy ratio, which is Rv/Rh, therefore varies from 1 to 100. Table 2 shows results with the same set of horizontal and vertical resistivities taken at a different relative deviation angle of 70°. Finally, Table 3 shows results at a deviation angle of 65°.
TABLE 1
80° deviation
Relative Relative
True Rh Inverted Rh error % True Rv Inverted Rv error %
1.00 0.9994 0.0552 1.00 1.0122 −1.2213
1.00 0.9999 0.0078 4.00 4.0393 −0.9822
1.00 0.9971 0.2941 9.00 9.0818 −0.9087
1.00 0.9940 0.6010 16.00 16.1234 −0.7713
1.00 0.9914 0.8641 25.00 25.1666 −0.6666
1.00 0.9894 1.0624 36.00 36.2600 −0.7223
1.00 0.9877 1.2329 49.00 49.3255 −0.6643
1.00 0.9862 1.3839 64.00 64.4313 −0.6739
1.00 0.9848 1.5231 81.00 81.4893 −0.6041
1.00 0.9835 1.6496 100.00 100.5426 −0.5426
TABLE 2
70° deviation
Relative Relative
True Rh Inverted Rh error % True Rv Inverted Rv error %
1.00 0.9979 0.2149 1.00 1.0149 −1.4936
1.00 0.9987 0.1304 4.00 4.0361 −0.9016
1.00 0.9993 0.0724 9.00 9.0801 −0.8896
1.00 0.9998 0.0176 16.00 16.1265 −0.7908
1.00 1.0005 −0.0478 25.00 25.2483 −0.9931
1.00 1.0012 −0.1175 36.00 36.3221 −0.8948
1.00 1.0019 −0.1930 49.00 49.3018 −0.6159
1.00 1.0027 −0.2707 64.00 64.5275 −0.8243
1.00 1.0035 −0.3515 81.00 81.4058 −0.5010
1.00 1.0044 −0.4361 100.00 100.6267 −0.6267
TABLE 3
65° deviation
Relative Relative
True Rh Inverted Rh error % True Rv Inverted Rv error %
1.00 0.9964 0.3605 1.00 1.0174 −1.7446
1.00 1.0104 −1.0360 4.00 3.9736 0.6588
1.00 1.0243 −2.4329 9.00 8.8686 1.4606
1.00 1.0382 −3.8234 16.00 15.5536 2.7898
1.00 1.0523 −5.2268 25.00 23.8972 4.4112
1.00 1.0664 −6.6401 36.00 33.6467 6.5370
1.00 1.0807 −8.0667 49.00 44.7809 8.6104
1.00 1.0951 −9.5077 64.00 56.9753 10.9761
1.00 1.1096 −10.9643 81.00 69.9713 13.6157
1.00 1.1243 −12.4271 100.00 83.5903 16.4097
For the relative deviation of 80 degrees and above (Table 1), the relative errors of inverted horizontal and vertical resistivities are less than or approximately 1%. For the relative deviation of 70 degrees, (Table 2) the relative errors of inverted horizontal and vertical resistivities are less than 2%. For a relative deviation of 65 degrees (Table 3), the relative errors increase with an increase in the anisotropy ratio. In most real situations, the anisotropy ratio will be less than 10. Similar results can be obtained through an increase in the horizontal resistivity. The anisotropy ratio is a determining factor.
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 such variations within the scope and spirit of the appended claims be embraced by the foregoing disclosure.

Claims (40)

1. A method of logging a subsurface formation comprising a plurality of layers, each having a horizontal resistivity and a vertical resistivity, the method comprising;
(a) conveying an electromagnetic logging tool into a deviated borehole in the subsurface formation;
(b) obtaining principal component measurements from measurements made by said logging tool;
(c) obtaining an initial model of said horizontal and vertical resistivities;
(d) using said model and only a subset of the principal components measurements for determining an updated value of said horizontal resistivity;
(e) obtaining from only another subset of the principal component measurements and said updated value of said horizontal resistivity an updated value of said vertical resistivity.
2. The method of claim 1 wherein obtaining said principal component measurements further comprises:
(i) making measurements with said electromagnetic tool at a plurality of transmitter-receiver orientations, and
(ii) performing a coordinate transformation thereon.
3. The method of claim 1 wherein obtaining said initial model further comprises applying a skin effect correction.
4. The method of claim 1 wherein said obtaining said initial model further comprises assuming a value for a ratio of said horizontal and vertical conductivity.
5. The method of claim 1 wherein obtaining said principal component measurements further comprises making measurements with said electromagnetic tool at a plurality of transmitter-receiver spacings.
6. The method of claim 1 wherein obtaining said principal component measurements further comprises making measurements with said electromagnetic tool at a plurality of frequencies.
7. The method of claim 6 wherein obtaining said principal component measurements further comprises applying a frequency focusing.
8. The method of claim 6 wherein the number of said frequencies is two.
9. The method of claim 1 further comprising conveying said electromagnetic logging tool into the borehole on one of (i) a wireline, and (ii) a drilling tubular.
10. The method of claim 1 further comprising repeating (a)–(e) at a plurality of depths in the borehole.
11. The method of claim 7 wherein applying said frequency focusing further comprises using a Taylor series expansion.
12. The method of claim 1 wherein said borehole axis is deviated by about 65° or more to a normal to a layer boundary in said formation.
13. The method of claim 12 wherein said borehole is substantially horizontal.
14. The method of claim 1 wherein the subset comprises at least one of (i) an xx-measurement, and (ii) a yy-measurement, and the another subset comprises a zz-measurement.
15. The method of claim 1 further comprising using an orientation sensor for determining a toolface orientation of said logging tool.
16. The method of claim 14 wherein obtaining said initial model further comprises using said xx-measurement.
17. The method of claim 14 wherein obtaining said at least one additional measurement further comprises obtaining both xx- and y-measurements.
18. The method of claim 17 wherein determining said horizontal resistivity further comprises using a sum of said xx- and yy-measurements.
19. The method of claim 15 wherein said orientation sensor is selected from the group consisting of (i) a magnetometer, (ii) a gyroscope, (iii) an accelerometer, and (iv) an inclinometer.
20. The method of claim 1 further comprising using a processor for performing at least one of operations selected from (b),(c),(d) and (e).
21. The method of claim 19 wherein said processor is located at a downhole location.
22. The method of claim 19 wherein said processor performs said at least one operation substantially in real time.
23. A system for logging a subsurface formation comprising a plurality of layers each having a horizontal resistivity and a vertical resistivity, the system comprising:
(a) a logging tool conveyed in a deviated borehole in the earth formation, the logging tool having a plurality of transmitters and receivers, at least one of said plurality of transmitters and receivers having an axis inclined to a longitudinal axis of said logging tool, said transmitters conveying an electromagnetic signal into said subsurface formation, said receivers receiving resulting signals from an interaction of said transmitted signals with said subsurface formation;
(b) a processor which determines principal component measurements from said received signals;
(c) a processor which obtains an initial model of said horizontal and vertical resistivities;
(d) a processor which uses said model and only a subset of the principal component measurements and determines an updated value of said horizontal resistivity;
(e) a processor which determines only from another subset of the principal component measurements and said updated value of said horizontal resistivity an updated value of said vertical resistivity.
24. The system of claim 23 wherein the same processor is used for (b)–(e).
25. The system of claim 24 wherein said same processor is at a downhole location.
26. The system of claim 23 wherein said processor in (d) performs a coordinate transformation thereon.
27. The system of claim 23 wherein obtaining said initial model further comprises applying a skin effect correction.
28. The system of claim 23 wherein said received signals are associated with a plurality of transmitter-receiver spacings.
29. The system of claim 23 wherein said transmitters operate at a plurality of frequencies.
30. The system of claim 24 wherein said transmitters operate at a plurality of frequencies.
31. The system of claim 30 wherein said same processor further applies a frequency focusing.
32. The system of claim 30 wherein said plurality of frequencies is two.
33. The system of claim 23 further comprising a conveyance device selected from (i) a wireline, and (ii) a drilling tubular for conveying said logging tool into a borehole in said subsurface formation.
34. The system of claim 31 wherein applying said frequency focusing further comprises using a Taylor series expansion.
35. The system of claim 23 further comprising an orientation sensor which determines a toolface orientation of said logging tool.
36. The apparatus of claim 23 wherein the subset comprises at least one of (i) an xx-measurement, and (ii) a yy-measurement, and the another subset comprises a zz-measurement.
37. The system of claim 35 wherein said orientation sensor is selected from the group consisting of (i) a magnetometer, (ii) a gyroscope, (iii) an accelerometer, and, (iv) an inclinometer.
38. The system of claim 36 wherein determining said updated value of said horizontal resistivity further comprises using a sum of said xx- and yy-measurements.
39. The system of claim 24 wherein said same processor performs (b)–(e) substantially in real time.
40. A system for determining in real time a horizontal conductivity and a vertical conductivity of an subsurface formation having a plurality of layers, the system comprising:
(a) a logging tool having a plurality of transmitters and receivers, at least one of said plurality of transmitters and receivers having an axis inclined to a longitudinal axis of said logging tool, said transmitters conveying an electromagnetic signal into said subsurface formation, said receivers receiving signals resulting from an interaction of said transmitted signals with said subsurface formation;
(b) a processor which determines from said received signals zz-measurements and at least one additional measurement selected from an xx-measurement and a yy-measurement (principal components);
(c) a processor which determines in real time said horizontal and vertical resistivity in real time only from said principal component measurements.
US10/652,762 2003-08-29 2003-08-29 Real time processing of multicomponent induction tool data in highly deviated and horizontal wells Expired - Lifetime US7043370B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US10/652,762 US7043370B2 (en) 2003-08-29 2003-08-29 Real time processing of multicomponent induction tool data in highly deviated and horizontal wells
PCT/US2004/027751 WO2005024467A1 (en) 2003-08-29 2004-08-26 Real time processing of multicomponent induction tool data in highly deviated and horizontal wells

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US10/652,762 US7043370B2 (en) 2003-08-29 2003-08-29 Real time processing of multicomponent induction tool data in highly deviated and horizontal wells

Publications (2)

Publication Number Publication Date
US20050049792A1 US20050049792A1 (en) 2005-03-03
US7043370B2 true US7043370B2 (en) 2006-05-09

Family

ID=34217731

Family Applications (1)

Application Number Title Priority Date Filing Date
US10/652,762 Expired - Lifetime US7043370B2 (en) 2003-08-29 2003-08-29 Real time processing of multicomponent induction tool data in highly deviated and horizontal wells

Country Status (2)

Country Link
US (1) US7043370B2 (en)
WO (1) WO2005024467A1 (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050256645A1 (en) * 2004-05-11 2005-11-17 Baker Hughes Incorporated Determination of fracture orientation and length using multi-component and multi-array induction data
US20060173624A1 (en) * 2005-01-31 2006-08-03 Baker Hughes Incorporated Method for real-time well-site interpretation of array resistivity log data in vertical and deviated wells
US20070234787A1 (en) * 2006-04-07 2007-10-11 Baker Hughes Incorporated Processing of multi-component induction data in the presence of borehole abnormalities
US9927551B2 (en) 2013-12-27 2018-03-27 Halliburton Energy Services, Inc. Multifrequency processing to determine formation properties
US10317565B2 (en) 2016-03-10 2019-06-11 Halliburton Energy Services, Inc. Induction logging borehole correction for water-based mud
US11175430B1 (en) 2020-05-19 2021-11-16 Halliburton Energy Services, Inc. Processing nuclear magnetic resonance signals in a downhole environment
US11487042B2 (en) * 2008-04-17 2022-11-01 Richard H. Hardman Methods for producing a log of material properties

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7150316B2 (en) * 2004-02-04 2006-12-19 Baker Hughes Incorporated Method of eliminating conductive drill parasitic influence on the measurements of transient electromagnetic components in MWD tools
EP1718999A1 (en) * 2004-02-23 2006-11-08 Oxford Geoservice Limited Method and apparatus for gradient electromagnetic induction well logging
EP1577683B1 (en) * 2004-03-16 2008-12-17 Services Petroliers Schlumberger Characterizing properties of a geological formation by coupled acoustic and electromagnetic measurements
US7317991B2 (en) * 2005-01-18 2008-01-08 Baker Hughes Incorporated Multicomponent induction measurements in cross-bedded and weak anisotropy approximation
MX2007015622A (en) * 2005-06-09 2008-02-25 Exxonmobil Upstream Res Co Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys.
US20120209528A1 (en) * 2011-02-10 2012-08-16 Baker Hughes Incorporated Inversion-Based Method to Correct for the Pipe Residual Signal in Transient MWD Measurements
NO346993B1 (en) * 2013-04-02 2023-03-27 Schlumberger Technology Bv Extended 1D inversion of electromagnetic measurements for subterranean formation evaluation
US9611731B2 (en) 2013-10-04 2017-04-04 Halliburton Energy Services, Inc. Determination of formation dip/azimuth with multicomponent induction data
WO2015053741A1 (en) * 2013-10-07 2015-04-16 Halliburton Energy Services, Inc. Multi-component induction logging methods and systems having a trend-based data quality indicator
EP3069173A1 (en) * 2014-01-13 2016-09-21 Halliburton Energy Services, Inc. Dip correction using estimated formation layer resistivities
US10401529B2 (en) 2014-10-17 2019-09-03 Halliburton Energy Services, Inc. Fast-changing dip formation resistivity estimation
US10935690B2 (en) * 2015-11-12 2021-03-02 Schlumberger Technology Corporation Methods for inverting electromagnetic logging measurements
US20190086575A1 (en) * 2016-07-01 2019-03-21 Halliburton Energy Services, Inc. Real-Time Self-Consistency Quality Indicators for Multi-Component Induction Tools
CN111127243B (en) * 2018-10-31 2023-12-01 中国石油化工集团有限公司 Method for correcting resistivity transfer curve of electromagnetic wave while drilling

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4360777A (en) * 1979-12-31 1982-11-23 Schlumberger Technology Corporation Induction dipmeter apparatus and method
US4837517A (en) 1987-07-16 1989-06-06 Schlumberger Technology Corporation Spatial frequency method and apparatus for investigating earth conductivity with high vertical resolution by induction techniques
US5157605A (en) 1987-04-27 1992-10-20 Schlumberger Technology Corporation Induction logging method and apparatus including means for combining on-phase and quadrature components of signals received at varying frequencies and including use of multiple receiver means associated with a single transmitter
US5446654A (en) 1992-07-16 1995-08-29 Halliburton Company Shoulder effect logging method
US5452761A (en) 1994-10-31 1995-09-26 Western Atlas International, Inc. Synchronized digital stacking method and application to induction logging tools
US5774360A (en) 1996-03-26 1998-06-30 Western Atlas International, Inc. Method of correcting formation resistivity well logs for the effects of formation layer inclination with respect to the wellbore
US5842149A (en) * 1996-10-22 1998-11-24 Baker Hughes Incorporated Closed loop drilling system
US5867806A (en) 1996-03-13 1999-02-02 Halliburton Energy Services, Inc. System and method for performing inversion on LWD resistivity logs with enhanced resolution
US5999883A (en) 1996-07-26 1999-12-07 Western Atlas International, Inc. Conductivity anisotropy estimation method for inversion processing of measurements made by a transverse electromagnetic induction logging instrument
US6049209A (en) 1998-05-14 2000-04-11 Western Atlas International, Inc. Axial resolution enhancement method for resistivity well logs which accounts for lateral variations in the resistivity of earth formations penetrated by a wellbore
US6147496A (en) 1996-07-01 2000-11-14 Shell Oil Company Determining electrical conductivity of a laminated earth formation using induction logging
US6206108B1 (en) * 1995-01-12 2001-03-27 Baker Hughes Incorporated Drilling system with integrated bottom hole assembly
US6233524B1 (en) * 1995-10-23 2001-05-15 Baker Hughes Incorporated Closed loop drilling system
US6393364B1 (en) * 2000-05-30 2002-05-21 Halliburton Energy Services, Inc. Determination of conductivity in anisotropic dipping formations from magnetic coupling measurements
US6584408B2 (en) * 2001-06-26 2003-06-24 Schlumberger Technology Corporation Subsurface formation parameters from tri-axial measurements
US6643589B2 (en) * 2001-03-08 2003-11-04 Baker Hughes Incorporated Simultaneous determination of formation angles and anisotropic resistivity using multi-component induction logging data
US6727696B2 (en) * 1998-03-06 2004-04-27 Baker Hughes Incorporated Downhole NMR processing
US6794875B2 (en) * 2002-05-20 2004-09-21 Halliburton Energy Services, Inc. Induction well logging apparatus and method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6636045B2 (en) * 2001-04-03 2003-10-21 Baker Hughes Incorporated Method of determining formation anisotropy in deviated wells using separation of induction mode

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4360777A (en) * 1979-12-31 1982-11-23 Schlumberger Technology Corporation Induction dipmeter apparatus and method
US5157605A (en) 1987-04-27 1992-10-20 Schlumberger Technology Corporation Induction logging method and apparatus including means for combining on-phase and quadrature components of signals received at varying frequencies and including use of multiple receiver means associated with a single transmitter
US4837517A (en) 1987-07-16 1989-06-06 Schlumberger Technology Corporation Spatial frequency method and apparatus for investigating earth conductivity with high vertical resolution by induction techniques
US5446654A (en) 1992-07-16 1995-08-29 Halliburton Company Shoulder effect logging method
US5452761A (en) 1994-10-31 1995-09-26 Western Atlas International, Inc. Synchronized digital stacking method and application to induction logging tools
US6206108B1 (en) * 1995-01-12 2001-03-27 Baker Hughes Incorporated Drilling system with integrated bottom hole assembly
US6233524B1 (en) * 1995-10-23 2001-05-15 Baker Hughes Incorporated Closed loop drilling system
US5867806A (en) 1996-03-13 1999-02-02 Halliburton Energy Services, Inc. System and method for performing inversion on LWD resistivity logs with enhanced resolution
US5774360A (en) 1996-03-26 1998-06-30 Western Atlas International, Inc. Method of correcting formation resistivity well logs for the effects of formation layer inclination with respect to the wellbore
US6147496A (en) 1996-07-01 2000-11-14 Shell Oil Company Determining electrical conductivity of a laminated earth formation using induction logging
US5999883A (en) 1996-07-26 1999-12-07 Western Atlas International, Inc. Conductivity anisotropy estimation method for inversion processing of measurements made by a transverse electromagnetic induction logging instrument
US5842149A (en) * 1996-10-22 1998-11-24 Baker Hughes Incorporated Closed loop drilling system
US6727696B2 (en) * 1998-03-06 2004-04-27 Baker Hughes Incorporated Downhole NMR processing
US6049209A (en) 1998-05-14 2000-04-11 Western Atlas International, Inc. Axial resolution enhancement method for resistivity well logs which accounts for lateral variations in the resistivity of earth formations penetrated by a wellbore
US6393364B1 (en) * 2000-05-30 2002-05-21 Halliburton Energy Services, Inc. Determination of conductivity in anisotropic dipping formations from magnetic coupling measurements
US6643589B2 (en) * 2001-03-08 2003-11-04 Baker Hughes Incorporated Simultaneous determination of formation angles and anisotropic resistivity using multi-component induction logging data
US6584408B2 (en) * 2001-06-26 2003-06-24 Schlumberger Technology Corporation Subsurface formation parameters from tri-axial measurements
US6794875B2 (en) * 2002-05-20 2004-09-21 Halliburton Energy Services, Inc. Induction well logging apparatus and method

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050256645A1 (en) * 2004-05-11 2005-11-17 Baker Hughes Incorporated Determination of fracture orientation and length using multi-component and multi-array induction data
US7359800B2 (en) * 2004-05-11 2008-04-15 Baker Hughes Incorporated Determination of fracture orientation and length using multi-component and multi-array induction data
US20060173624A1 (en) * 2005-01-31 2006-08-03 Baker Hughes Incorporated Method for real-time well-site interpretation of array resistivity log data in vertical and deviated wells
US7313479B2 (en) * 2005-01-31 2007-12-25 Baker Hughes Incorporated Method for real-time well-site interpretation of array resistivity log data in vertical and deviated wells
US20070234787A1 (en) * 2006-04-07 2007-10-11 Baker Hughes Incorporated Processing of multi-component induction data in the presence of borehole abnormalities
US8931335B2 (en) * 2006-04-07 2015-01-13 Baker Hughes Incorporation Processing of multi-component induction data in the presence of borehole abnormalities
US11487042B2 (en) * 2008-04-17 2022-11-01 Richard H. Hardman Methods for producing a log of material properties
US9927551B2 (en) 2013-12-27 2018-03-27 Halliburton Energy Services, Inc. Multifrequency processing to determine formation properties
US10317565B2 (en) 2016-03-10 2019-06-11 Halliburton Energy Services, Inc. Induction logging borehole correction for water-based mud
US11175430B1 (en) 2020-05-19 2021-11-16 Halliburton Energy Services, Inc. Processing nuclear magnetic resonance signals in a downhole environment

Also Published As

Publication number Publication date
US20050049792A1 (en) 2005-03-03
WO2005024467A1 (en) 2005-03-17

Similar Documents

Publication Publication Date Title
US7043370B2 (en) Real time processing of multicomponent induction tool data in highly deviated and horizontal wells
US7612566B2 (en) Method and apparatus for the use of multicomponent induction tool for geosteering and formation resistivity data interpretation in horizontal wells
US7421345B2 (en) Geosteering in earth formations using multicomponent induction measurements
US7375530B2 (en) Method for signal enhancement in azimuthal propagation resistivity while drilling
US7269515B2 (en) Geosteering in anisotropic formations using multicomponent induction measurements
US6502036B2 (en) 2-D inversion of multi-component induction logging data to resolve anisotropic resistivity structure
US8049507B2 (en) Transient EM for geosteering and LWD/wireline formation evaluation
US10768336B2 (en) Formation logging using multicomponent signal-based measurement of anisotropic permittivity and resistivity
US6925384B2 (en) Method for resistivity anisotropy determination in conductive borehole environments
US7471088B2 (en) Elimination of the anisotropy effect in LWD azimuthal resistivity tool data
US8305081B2 (en) Cancellation of vibration noise in deep transient resistivity measurements while drilling
US20070216416A1 (en) Electromagnetic and Magnetostatic Shield To Perform Measurements Ahead of the Drill Bit
US8060310B2 (en) Geosteering in earth formations using multicomponent induction measurements
US20060197533A1 (en) Method for measuring transient electromagnetic components to perform deep geosteering while drilling
US20070024286A1 (en) Compensation for tool disposition in LWD resistivity measurements
WO2012129058A2 (en) Correction of deep azimuthal resistivity measurements for bending
EP1483602B1 (en) Use of a multicomponent induction tool for geosteering and formation resistivity data interpretation in horizontal wells
US20060192560A1 (en) Well placement by use of differences in electrical anisotropy of different layers
EP1875275A1 (en) Geosteering in anisotropic formations using multicomponent induction measurements

Legal Events

Date Code Title Description
AS Assignment

Owner name: BAKER HUGHES INCORPORATED, TEXAS

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:YU, LIMING;KRIEGSHAUSER, BERTHOLD;FANINI, OTTO;REEL/FRAME:014772/0184;SIGNING DATES FROM 20030902 TO 20030908

FEPP Fee payment procedure

Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

STCF Information on status: patent grant

Free format text: PATENTED CASE

FPAY Fee payment

Year of fee payment: 4

FPAY Fee payment

Year of fee payment: 8

MAFP Maintenance fee payment

Free format text: PAYMENT OF MAINTENANCE FEE, 12TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1553)

Year of fee payment: 12