US20010039477A1 - Inhomogeneous background based focusing method for multiarray induction measurements in a deviated well - Google Patents

Inhomogeneous background based focusing method for multiarray induction measurements in a deviated well Download PDF

Info

Publication number
US20010039477A1
US20010039477A1 US09798867 US79886701A US2001039477A1 US 20010039477 A1 US20010039477 A1 US 20010039477A1 US 09798867 US09798867 US 09798867 US 79886701 A US79886701 A US 79886701A US 2001039477 A1 US2001039477 A1 US 2001039477A1
Authority
US
Grant status
Application
Patent type
Prior art keywords
focusing
response
formation
background
induction
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.)
Granted
Application number
US09798867
Other versions
US6442488B2 (en )
Inventor
Jaiqi Xiao
Ingo Geldmacher
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 Inc
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

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS
    • 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

Abstract

In highly deviated wells or when the relative dip angle between formation layering and wellbore axis is large, array induction measurements exhibit erratic spikes, misleading curve separations, and inaccurate resistivity values, preventing log analysts from accurately evaluating invasion and formation resistivities. To address these problems, various correction methods and inversion techniques have been developed. The correction methods, however, only yield satisfactory results when the relative dip angle is low to moderate, and inversion techniques are typically very time consuming. Conventionally, dipping bed effects are considered in terms of a charge effect and a volumetric effect. As even more complex earth models are considered, we find that formation anisotropy also exaggerates the dipping effect, manifested by misleading curve separations in the array instrument readings. Our newly developed deviated-well software focusing (DSF) method simultaneously accounts for all these dipping effects. The DSF method is derived from the Born approximation. The induction response is separated into two portions: a background response and a perturbation response. An inhomogeneous, anisotropic background formation model is used to calculate the background response, and the perturbation response is interpreted through a software focusing technique. The combination of the two solutions is the final result. We show that our method significantly reduces spurious separation between shallow and deep reading curves and minimizes the confusion between apparent and real invasion.

Description

    CROSS REFERENCES TO RELATED APPLICATIONS
  • [0001]
    This application is a Continuation-in-Part of U.S. patent application Ser. No. 09/264,105 filed on Mar. 8, 1999, now U.S. Pat. No.______ the contents of which are fully incorporated herein by reference.
  • FIELD OF THE INVENTION
  • [0002]
    The present invention relates to electromagnetic induction well logging. More specifically, the present invention relates to an improved method of focusing the measurements obtained with array induction tools in a deviated borehole.
  • BACKGROUND OF THE INVENTION
  • [0003]
    The production of hydrocarbons from subsurface formations typically commences by forming a borehole through the earth to a subsurface reservoir thought to contain hydrocarbons. From the borehole, various physical, chemical, and mechanical properties are “logged” for the purpose of determining the nature and characteristics, including for example, the porosity, permeability, saturation, and depth, of the subsurface formations encountered. One such logging technique commonly used in the industry is referred to as induction logging. Induction logging measures the conductivity or its inverse, the resistivity, of a formation. Formation conductivity is one possible indicator of the presence or absence of a significant accumulation of hydrocarbons because, generally speaking, hydrocarbons are relatively poor conductors of electricity. Formation water, on the other hand, typically salty, is a relatively good conductor of electricity. Thus, induction logging tools can obtain information that, properly interpreted, indicates the presence or absence of hydrocarbons.
  • [0004]
    These induction (also known as electromagnetic induction) well logging instruments were first introduced by Doll, H. G., “Introduction to Induction Logging and Application to Logging of Wells Drilled with Oil Based Mud,”Journal of Petroleum Technology, June, 1949, pp. 148-62. Induction well logging instruments typically include a sonde having one or more transmitter coils and one or more receiver coils at axially spaced apart locations. Induction well logging instruments also typically include a source of alternating current (AC) which is conducted through the transmitter coils. The AC passing through the transmitter coils induces a magnetic field within the surrounding formation, causing the flow of eddy currents within the earth formations. In general, the magnitude of the eddy currents is proportional to the electrical conductivity (the inverse of the electrical resistivity) of the earth formations surrounding the instrument. The eddy currents, in turn, induce a magnetic field that is coupled to the receiver coil, thereby inducing in the receiver coil a voltage signal with magnitude and phase dependent upon the electrical characteristics of the adjacent formation.
  • [0005]
    Typically, the signal from the receiver coil is applied to one or more phase detection circuits, each of which generates a signal proportional to the magnitude of that component of the receiver coil signal having a particular, predetermined phase. Thus, one such phase detector circuit senses the magnitude of the component of the receiver coil signal that is in-phase with the transmitter current in the transmitter coil. This component signal is commonly referred to as the real or inphase (R) component. A second phase detection circuit commonly used in induction logging tools detects the component of the receiver coil signal that is 90° out of phase with the transmitter current.
  • [0006]
    This latter component signal is commonly referred to as the quadrature-phase (X) component signal. Because the output signal from the receiver coil is not itself an absolute measure of conductivity, but rather is merely proportional to the true formation conductivity, the output signal must be processed to obtain a log or plot of the true formation conductivity as a function of axial depth in the borehole. Most modem theoretical analysis of induction log processing is based upon the work of H. G. Doll which is summarized in his 1949 article. According to Doll's analysis, the in-phase component of the signal induced in the receiver coil is directly proportional to the conductivity of the surrounding formation, and the constant of proportionality, referred to by Doll as the “geometrical factor,” is a function of the geometry of the tool as it relates to the portion of the formation being measured.
  • [0007]
    Doll calculated what he termed the “unit geometrical factor,” which defines the relationship between the conductivity of a so-called “unit ground loop,” a horizontal loop of homogeneous formational material having a circular shape with its center on the axis of the borehole and having a very small, square cross section, and the elementary voltage signal contributed by the unit ground loop to the total in-phase signal induced in the receiver coil. By integrating the unit geometrical factor across all unit ground loops lying within a horizontal plane spaced at some axial distance z from the center of the coil system, Doll obtained the geometrical factor for a “unit bed.” A plot of this geometrical factor as a function of the axial distance from the center of the coil system gives what is commonly referred to as the “vertical geometrical factor” for the tool. It is an accurate plot of the sonde response function relating formation conductivity to output voltage measurements for the tool, assuming no attenuation or phase shift of the induced magnetic field as a consequence of the conductivity of the surrounding formation.
  • [0008]
    Induction logging technology has evolved significantly since its introduction by Doll. In recent years, induction devices consisting of several complex coil combinations have been replaced by tools with multiple arrays (see, for example, Beard, D. R., et at., “A New, Fully Digital, Full-spectrum Induction Device for Determining Accurate Resistivity with Enhanced Diagnostics and Data Integrity Verification,” SPWLA 37th Annual Logging Symposium, June, 1996, Paper B; Beard, D. R., et at., “Practical Applications of a New Multichannel and Fully Digital Spectrum Induction System,” SPE Annual Technical Conference and Exhibition, 1996, Paper No. 36504; and Barber, T. D., et at., “A Multiarray Induction Tool optimized for Efficient Wellsite Operation,” SPE 7oth Annual Technical Conference and Exhibition, 1995, Paper No. 30583). Each array consists of one transmitter and a pair of receiver coils. These new induction devices are commonly referred to as array-type induction tools.
  • [0009]
    A simple induction array (two-coil array and three-coil array) responds to all its surrounding media, including formations, the borehole, and invasion zones if there are any. This response will be degraded by severe borehole effect and will suffer from low vertical and radial resolution. In order to avoid the weaknesses of the simple induction arrays, array combinations are used to increase the response contribution from the medium of interest, such as uninvaded formation, and at the same time to reduce the response contribution from the medium of disinterest, such as the borehole. This process by which the output of an induction logging instrument is made to effectively zoom in on a specific space of its surrounding medium and mute the unwanted peripherals is referred to as focusing.
  • [0010]
    The older style tools attempt to focus the tool response using carefully selected coil arrangements. The focusing therefore is fixed by the tool design, i.e. these tools are “hardware-focused”. In array-type induction tools, the measurements from various arrays are combined through a software algorithm to achieve focusing of the signal response, i.e. these tools are “software-focused”. This processing produces a set of curves with predetermined depth of investigation, vertical resolution and other optimized 2-D features.
  • [0011]
    Using software-based focusing provides greater flexibility for handling various logging environments and for creating more reliable induction logs. However, the quality and accuracy of the final focused logs are dependant on the accuracy of the software focusing method.
  • [0012]
    Prior art focusing methods are based on a method that was proposed by Barber and Zhou (see Barber, T. D. and Rosthal, R. A., “Using a Multiarray Induction Tool to Achieve High-Resolution Logs with Minimum Environmental Effects,” SPE 66th Annual Technical Conference and Exhibition, 1991, Paper No. 22725 and Zhou, Q., Beard, D. and Tabrovsky L., “Numerical Focusing of Induction Logging Measurements,” 12th Workshop in Electromagnetic Induction in Earth,” August, 1994) and is, for reference purposes, here referred to as the “conventional focusing method”. The conventional focusing method can be expressed mathematically as σ TRF ( z ) = i = 1 m ary z = z mm z max W i ( z ) σ ai ( z - z ) ( 1 )
    Figure US20010039477A1-20011108-M00001
  • [0013]
    where σai is the measured log from the ith array; Wi is the focusing filter; mary is the total number of arrays; and Zmin and Zmax define the depth window surrounding the output point.
  • [0014]
    Theoretically, the software focusing method described by eq. (1) can be traced back to the Born approximation (a linear approximation of the measured response of a medium) and then the condition for eq. (1) is an assumption of an homogeneous background. This focusing method produces good quality focused logs when the formation conductivity varies with small to moderate contrasts between adjacent formation beds. However, when the formation conductivity varies with very large conductivity contrasts, i.e. if the formation is very “inhomogeneous”, the focused logs are not as good as would be expected.
  • [0015]
    The root cause of this shortcoming is the nonlinearity of the induction response with respect to the formation conductivity. The basic assumption for the focusing algorithm expressed through eq. (1) is that the array measurements behave linearly with conductivity. The error due to the violation of this linearity assumption is referred to as the nonlinearity effect. The nonlinearity effect is formation-dependent: the larger the inhomogeneity, the stronger the nonlinearity effect. A focusing method based on a formation response with a homogeneous background, propagates or even amplifies the nonlinearity effect.
  • [0016]
    In a given logging environment, the inhomogeneity of a formation is described by numerous factors. Formation layering contributes to the vertical inhomogeneity, which is conventionally described by the Rt/Rs contrast, where Rt is the formation resistivity and Rs is the shoulder resistivity. Radial inhomogeneity is expressed through the Rt/Rx0 and the Rx0/Rm contrast, where Rx0 is the resistivity of the invaded zone and Rm is the mud resistivity. Other inhomogeneities are introduced through borehole irregularity, tool eccentricity, borehole deviation, etc. Despite the multitude of factors, the vertical inhomogeneity often dominates, particularly when the Rt/ Rs contrast is large.
  • [0017]
    The induction logging response function varies with formation inhomogeneity due to the nonlinearity of the induction measurements. The focusing filters are designed based on Born geometric factors, which equal the response functions under a homogeneous background. With such designed filters and eq. (1), the nonlinearity effect is propagated or even amplified through the focusing process, especially when the formation is inhomogeneous with a large conductivity contrast.
  • [0018]
    U.S. patent application Ser. No. 09/264,105 (now U.S. Pat. No.______) having the same assignee as the present application and the contents of which are fully incorporated herein by reference, teaches use of an an inhomogeneous background formation model. Using this inhomogeneous background formation model, the formation response of the induction tool measurement can be split into two portions: the response due to an inhomogeneous background conductivity, a background response, and a certain “response residue”. The background response is the computer simulated measurements of the inhomogeneous background model. The response residue is the difference between raw measurements and the background responses.
  • [0019]
    The '105 application uses an initial formation model estimated from raw array measurements or processed logs and used as the background conductivity model. For the background response, the focusing result can be directly obtained using focusing target functions instead of applying the conventional focusing processing. Therefore, the focusing result of the background response is ideal and free of any nonlinearity effect. The conventional focusing procedure is applied to the response residue. The final focusing response is obtained by adding the two focusing results. Due to the relatively small amplitude of the response residue, the nonlinearity effect introduced to the focused result is very small. Hence, the nonlinearity effect on the final focusing result is largely reduced. Thus, by introducing an inhomogeneous formation background model into the focusing algorithm, an improved focusing method is achieved having a reduction in the propagation of nonlinearity effects.
  • [0020]
    The teachings of the '105 application are, however, limited to near vertical wells (wells in which the borehole axis makes a small angle with the normal to the bedding planes). In highly deviated wells when the angle between the normal to the bedding planes and wellbore axis is large, array induction measurements exhibit erratic spikes, misleading curve separations, and inaccurate resistivity values, preventing log analysts from accurately evaluating invasion and formation resistivities. There is a need for a method of determination of formation resistivities in deviated boreholes. Such a method should preferably also account for any anisotropy in the resistivity of the formations. The present invention satisfies this need.
  • SUMMARY OF THE INVENTION
  • [0021]
    The present invention is a method and apparatus of using an induction logging tool in a deviated borehole within earth formations for determining an anisotropic conductivity of the formation. An anisotropic background model is estimated from borehole corrected data. An asymptotic conductivity of the anisotropic background model is determined based on the horizontal and vertical conductivities of the model and a known inclination of the tool to the normal to the bedding planes. A modeled response is calculated for this anisotropic background model using a forward modeling program and skin effect corrections are applied to this modeled response. A response residual is calculated as a difference between the skin-effect corrected response and the borehole corrected data. This residual is focused and focused residual is combined with a focused output of the background model. Where appropriate, the intermediate outputs are converted to a true vertical depth to account for the inclination of the borehole axis.
  • BRIEF DESCRIPTION OF THE FIGURES
  • [0022]
    [0022]FIG. 1 shows an induction instrument disposed in a wellbore penetrating earth formations.
  • [0023]
    [0023]FIG. 2 shows a model used to illustrate the method of the present invention
  • [0024]
    [0024]FIG. 3 is a flow chart illustrating steps of the present invention.
  • [0025]
    [0025]FIGS. 4a and 4 b shows the effect of anisotropy and deviation on resistivity measurements.
  • [0026]
    [0026]FIGS. 5a and 5 b show a comparison of results obtained using prior art methods and the method of the present invention on model data.
  • [0027]
    [0027]FIGS. 6a and 6 b illustrate the relative insensitivity of the method of the present invention to an incorrect inclination angle used in the analysis on model data..
  • [0028]
    FIGS. 7 shows a comparison of results obtained using prior art methods and the method of the present invention on field measurements.
  • DETAILED DESCRIPTION OF THE PRESENT INVENTION
  • [0029]
    Referring now to FIG. 1, an electromagnetic induction well logging instrument 10 is shown disposed in a wellbore 2 drilled through earth formations. The earth formations are shown generally at 4. The instrument 10 can be lowered into and withdrawn from the wellbore 2 by means of an armored electrical cable 6 or similar conveyance known in the art. The instrument 10 can be assembled from three subsections: an auxiliary electronics unit 14 disposed at one end of the instrument 10; a coil mandrel unit 8 attached to the auxiliary electronics unit 14; and a receiver/signal processing/telemetry electronics unit 12 attached to the other end of the coil mandrel unit 8, this unit 12 typically being attached to the cable 6.
  • [0030]
    The coil mandrel unit 8 includes induction transmitter and receiver coils, as will be further explained, for inducing electromagnetic fields in the earth formations 4 and for receiving voltage signals induced by eddy currents flowing in the earth formations 4 as a result of the electromagnetic fields induced therein.
  • [0031]
    The auxiliary electronics unit 14 can include a signal generator and power amplifiers (not shown) to cause alternating currents of selected frequencies to flow through transmitter coils in the coil mandrel unit 8.
  • [0032]
    The receiver/signal processing/telemetry electronics unit 12 can include receiver circuits (not shown) for detecting voltages induced in receiver coils in the coil mandrel unit 8, and circuits for processing these received voltages (not shown) into signals representative of the conductivities of various layers, shown as 4A through 4F of the earth formations 4. As a matter of convenience the receiver/signal processing/telemetry electronics unit 12 can include signal telemetry to transmit the conductivity-related signals to the earth's surface along the cable 6 for further processing, or alternatively can store the conductivity related signals in an appropriate recording device (not shown) for processing after the instrument 10 is withdrawn from the wellbore 2.
  • [0033]
    [0033]FIG. 2 shows the model used in the present invention. The model comprises layers 111 a, 111 b, 111 c. . . 111 i, . . . each having thickness hi. The horizontal and vertical resistivities are denoted by Rhl, Rvl, . . . Rhl, Rvl, or equivalently, by horizontal and vertical conductivities σvl, σvl, . . . σhl, σvl, with the normal to the layer boundaries inclined at an angle θ to the axis 321 of the borehole that defines the z- axis of the tool. In the method of the present invention, the angle θ is assumed to be a known quantity. In practice, this may be obtained by making measurements that give the relative dip of the formations. Micro-resistivity logs, dipmeters or other measurements may be used for the purpose. It is commonly observed that the normal to the layer boundaries defines a symmetry axis for the anisotropy of the formations. The horizontal and vertical resistivities (or conductivities) are defined with respect to directions parallel to the layer boundaries and normal to the layer boundaries.
  • [0034]
    The focusing filter for the conventional focusing method given by eq. (1) is obtained by solving the following problem: g p ( ρ , z ) = i = 1 m ary z = z max z max W i ( z ) g B i ( ρ , z - z ) ( 2 )
    Figure US20010039477A1-20011108-M00002
  • [0035]
    with certain constraints. In eq. (2) gp is a 2-D response function defining ideal 2-D features of the focusing and is often referred to as the “focusing 2-D target function”; gi B is the Born 2-D geometric factor of the i-th array. The basic assumption for the focusing algorithm expressed through eq. (1) is the Born-Doll approximation, σ a i ( z ) = 0 ρ - z g B ι ( ρ , z - z ) σ ( ρ , z ) ( 3 )
    Figure US20010039477A1-20011108-M00003
  • [0036]
    where σl a is the measurement of the i-th array, and σ(ρ′,z′) is the 2-D formation conductivity distribution.
  • [0037]
    The inhomogeneous background-based focusing (IBF) method taught in the '105 application stems from the Born Approximation and is mathematically expressed as σ a ( z ) = 0 ρ - z g B ( ρ , z - z ) [ σ ( ρ , z ) - σ b ( ρ , z ) ] + σ ba ( z ) ( 4 )
    Figure US20010039477A1-20011108-M00004
  • [0038]
    where σba is the induction response of the background, σb is the background conductivity. The difference between the true formation conductivity an d the background conductivity, σ−σb, is referred to as the “conductivity perturbation.” By introducing a background, the Born approximation separates the induction response into two portions: the background contribution and the perturbation contribution. If the background conductivity is properly chosen, the Born approximation can better represent the induction measurements than the Doll approximation and the Born-Doll approximation.
  • [0039]
    The IBF takes the following form:
  • σp(z)=σbp(z)+Δσp(z)  (5)
  • [0040]
    where σbp is the focusing result of the background model, something that is easily calculated from the background model, and Δσp is the focusing result of the perturbation, expressed as: Δ σ p ( z ) = i = 1 m ary z = z min z max W i ( z ) [ σ a i ( z - z ) - σ ba i ( z - z ) ] . ( 6 )
    Figure US20010039477A1-20011108-M00005
  • [0041]
    Because σbp is directly determined, the focusing filter W only applies to the residue, [σl a(Z)−σl ba(Z)]. Thus, the error propagation of the software focusing is largely reduced and the final focusing results are improved
  • [0042]
    The physics of induction measurements in a deviated well can be mathematically expressed by the following integral equation E ( r ) = E 0 ( r ) + r g 0 ( r , r ) [ k 2 ( r ) - k b 2 ] E ( r ) + r g 0 ( r , r ) { E ( r ) · [ k 2 ( r ) ] k 2 ( r ) } ( 7 )
    Figure US20010039477A1-20011108-M00006
  • [0043]
    where E is the electric field; r is the location of the observation point; k is the wave number of the medium; gO is the Green's function. From this equation two basic dipping effects are derived: the charge effect and the volumetric effect. The charge effect corresponds to the third term of the integral equation. It is caused by the charge build-up while the induced eddy current flows across the formation interfaces. The volumetric effect relates to the second term. When the dip angle is relatively large (>30°), the software focusing of array induction measurements amplifies these dipping effects.
  • [0044]
    The software focusing of array induction measurements in deviated wells also faces a conceptual dilemma. The major functionality of the software focusing of array induction measurements is to enhance the response contribution of the target area, which is predefined by the required vertical resolution and the depth of investigation. As the depth of investigation increases, the target area moves deeper into the formation in the direction perpendicular to the tool axis, referred to as the “focusing direction.” When the borehole is perpendicular to the formation layering, i.e., the relative dip angle equals zero, no other boundaries cross the focusing direction except for the invasion frontiers. In this case any separation of the focused curves indicates the presence of invasion. While in deviated wells, formation interfaces cross the focusing direction. As a result, the separation of the focused curves no longer represents invasion only, making the interpretation of the focused curves impossible.
  • [0045]
    Besides the charge effect and the volumetric effect, the anisotropy is also an important dipping effect contributor. When the tool is located far away from formation boundaries or the formation is homogeneous, the induction measurement of a two-coil configuration in an anisotropic formation can be expressed as σ a 1 ω μ L 2 [ ( 2 - ik h L ) e ik h L - ik h Le ik h L β ] ( 8 )
    Figure US20010039477A1-20011108-M00007
  • [0046]
    where ω is the angular frequency, μ is the magnetic permeability, L is the spacing between the transmitter and the receiver, and kh is defined as k h = ( 1 + i ) ω μ σ h 2
    Figure US20010039477A1-20011108-M00008
  • [0047]
    where σh is the horizontal conductivity and β is defined as β = sin 2 θ + λ 2 cos 2 θ λ ( 8 a )
    Figure US20010039477A1-20011108-M00009
  • [0048]
    and λ is defined as the square root of the ratio of the horizontal conductivity to the vertical conductivity.
  • [0049]
    Eq. (8) leads to the observation that the effect of anisotropy in an induction measurement is a function of the frequency and the transmitter-receiver spacing. Because the skin-effect correction and the standard software focusing assume the formation is isotropic, the anisotropy effect propagates into the final focusing results.
  • [0050]
    In terms of the amplitude of induction measurements, Moran and Gianzero (1979) derived that a co-axial induction tool reads an asymptotic conductivity σe given by:
  • σe=βσh  (9)
  • [0051]
    The charge effect, the volumetric effect, and the anisotropy effect all relate to the anisotropic, inhomogeneous background of the formation. Hence, with the Born approximation these dipping effects can be separated from the borehole- and invasion-related perturbation response if the background conductivities σbh σbv are properly determined and the relative dip angle is given.
  • [0052]
    Turning now to FIG. 3, a schematic flow chart of the method of the present invention is shown. Starting with the borehole corrected data 201, a background model is formulated 203. As would be recognized by those versed in the art, the better the background model, the better the results will be; however, the choice of the background formation model is not as critical in the method of the present invention as is the choice of an initial model in an iterative nonlinear optimization method.
  • [0053]
    Using a known inclination angle θ 205 (this may be obtained by high resolution acoustic or resistivity imaging tools that can determine the orientation of bedding planes relative to the borehole), an asymptotic conductivity model is obtained using eqs. (8 a) and (9). With the asymptotic conductivity, the background model may be treated as a deviated isotropic model. This model is then converted to true vertical depth 209 to give a conductivity model as a function of the borehole depth (TVD). This is a straightforward stretching related to the angle θ. A response to the background model is generated 211 and corrected for skin effects 213 using known methods.
  • [0054]
    The borehole corrected data are compared with the skin-corrected model output and a response residue is determined 217 as a difference between the former and the latter. This residue is focused using eq. (6).
  • [0055]
    In parallel, the ideal focusing results for the model corrected for true vertical depth are determined 215. When the output of 215 is summed 221 with the output from 219, the end result is the desired focusing result 223.
  • EXAMPLES
  • [0056]
    Turning now to FIGS. 4a and 4 b, numerical results of forward modeling are shown. The model consists of two semi-infinite formations. The upper one is an isotropic formation of 1 Ohm-m resistivity and the lower one is an anisotropic formation with vertical resistivity of 10 Ohm-m and horizontal resistivity of 0.54 Ohm-m. The relative dip angle is 60 degrees. The resistivities of the lower half space are specifically chosen to have an asymptotic resistivity of 1 Ohm-m. FIG. 4a shows the HDIL skin-effect-corrected curves for seven different source-receiver distances. The ‘bump’ and the separation of the curves in the vicinity of the formation interface present a pronounced boundary-related dipping effect. Because the anisotropy effect differs with frequency and array spacing, a constant curve separation is seen in the anisotropic formation. These two anisotropy-related features propagate through the HDIL software focusing. The HDIL vertical resolution matching (VRM) curves shown in FIG. 4b clearly demonstrate this fact. FIGS. 4a and 4 b also show that a co-axial induction tool tends to read the asymptotic conductivity in an anisotropic formation, which is determined by eq. (9). Vertical Resolution Matching (VRM) is a process that enhances the curves to predetermined vertical resolutions of, for example 0.5 foot, 1 foot, 2 feet, and/or 4 feet. In FIGS. 4a, 4 b and also in the examples below, only the equivalent or asymptotic resistivity values are shown. If the anisotropy factor and the angle of inclination are known, then the vertical and horizontal resistivities may be recovered.
  • [0057]
    [0057]FIG. 5a shows the result of processing data from a second, layered model using the method described in the '105 application, i.e., assuming that the borehole is vertical and ignoring the effect of anisotropy. In the actual model, the borehole had a 70° inclination to the normal to the bedding plane. Again, the data for seven different source-receiver distances are displayed after VRM and clearly show that ignoring the effects of anisotropy and dip. In comparison, FIG. 5b shows the results of processing the same model output with the method of the present invention, i.e., including the effects of the deviated borehole and of anisotropy. The results in FIG. 5b are clearly superior to the prior art method.
  • [0058]
    The method of the present invention is relatively insensitive to errors in the assumed inclination angle. This is shown in FIGS. 6a and 6 b where the same model output for which results are shown in FIGS. 5a, 5 b were processed using the method of the present invention with an incorrect assumption of the inclination angle. FIG. 6a was obtained by assuming an angle of 650 whereas FIG. 6b was obtained assuming an angle of 75°, when, as noted above, the actual angle was 70°. The method of the present invention is seen to be robust for errors in the inclination angle of 5°.
  • [0059]
    [0059]FIG. 7 shows a field example. Caliper and gamma-ray logs are presented in the left track. The HDIL 1-ft VRM curves produced with the standard focusing method are shown in the middle track. The curve separation over the pay zones (e.g., at depths 630-650) and can be easily interpreted as resistive invasion. This interpretation, however, conflicts with local experiences. With a relative dip angle of 45 degrees, the method of the present invention produced a set of more reliable curves (in the right track). The value of the dip angle was later confirmed by image tools.
  • [0060]
    The present invention has been discussed above with respect to measurements made by a transverse induction logging tool conveyed on a wireline. This is not intended to be a limitation and the method is equally applicable to measurements made using a comparable tool conveyed on a measurement-while-drilling (MWD) assembly conveyed on a drill string or on coiled tubing.
  • [0061]
    While the foregoing disclosure is directed to the preferred embodiments of the invention, various modifications will be apparent to those skilled in the art. It is intended that all variations within the scope and spirit of the appended claims be embraced by the foregoing disclosure.

Claims (14)

  1. 1. A method of using an induction logging tool in a borehole within earth formations for determining conductivities thereof, the borehole axis having an inclination angle to a symmetry axis of said earth formations, the method comprising:
    (a) estimating an anisotropic background formation model from borehole corrected data, said anisotropic background model comprising horizontal and vertical resistivities;
    (b) determining asymptotic conductivities of the anisotropic background model relative to the borehole axis using an assumed value for said inclination angle;
    (c) obtaining a modeled response for said anisotropic background model with said asymptotic conductivities;
    (c) making skin-effect corrections to said modeled response giving a skin-effect corrected response;
    (d) calculating a response residual as a difference between the skin-effect corrected response and the borehole corrected data;
    (e) focusing said residual giving a focused residual;
    (f) calculating an idealized focused result for said background model; and
    (g) combining the results of said focusing of said residual with said ideal focused results.
  2. 2. The method of
    claim 1
    wherein obtaining the asymptotic conductivity of the anisotropic background model further comprises using the relation σe=βσh wherein σe is the asymptotic conductivity, σh is the horizontal conductivity,
    β = sin 2 θ + λ 2 cos 2 θ λ 2
    Figure US20010039477A1-20011108-M00010
    where θ is the assumed angle of inclination and λ is the square root of the ratio of the horizontal conductivity and vertical conductivity.
  3. 3. The method of
    claim 1
    wherein calculating said response residual further comprises converting the skin-effect corrected response to true vertical depth.
  4. 4. The method of
    claim 1
    wherein calculating the idealized focusing result for said background model further comprises converting the background model to true vertical depth.
  5. 5. The method of
    claim 1
    , wherein said response residual is smaller than said modeled response.
  6. 6. The method of
    claim 1
    wherein focusing said residual further comprises using linear focusing.
  7. 7. The method of
    claim 2
    wherein calculating said response residual further comprises converting the skin-effect corrected response to true vertical depth.
  8. 8. The method of
    claim 7
    wherein calculating the idealized focusing result for said background model further comprises converting the background model to true vertical depth.
  9. 9. The method of
    claim 1
    wherein the assumed angle of inclination differs from the angle of inclination by less than 5°.
  10. 10. A system for making geophysical measurements in a deviated borehole having an axis with an inclination angle to a symmetry axis of earth formations, the system comprising:
    (a) a logging tool to make geophysical measurements having at least one transmitter and at least one receiver array; and
    (b) a computing device programmed to
    (i) estimate an background formation model from said geophysical measurements,
    (ii) determine an asymptotic conductivity of the anisotropic background model using an assumed inclination angle,
    (iii) obtain a modeled response for said background model, giving a modeled response,
    (iii) make skin-effect corrections for said modeled response giving a skin-effect corrected response,
    (iv) calculate a response residual as a difference between the skin-effect corrected response and borehole corrected geophysical measurements,
    (v) focus said residual giving a focused residual,
    (vi) calculate an idealized focused result for said background model, and
    (vii) combine the results of said focusing of said residual with said ideal focused results.
  11. 11. The system of
    claim 10
    wherein the computing device is further programmed to obtain said asymptotic conductivity using the relation σe=βσh wherein σe is the asymptotic conductivity, σh is the horizontal conductivity,
    β = sin 2 θ + λ 2 cos 2 θ λ 2
    Figure US20010039477A1-20011108-M00011
    where θ is the assumed angle of inclination, and λ is the square root of the ratio of the horizontal conductivity and vertical conductivity.
  12. 12. The system of
    claim 10
    wherein the computing device is further programmed to convert the skin-effect corrected response to true vertical depth.
  13. 13. The system of
    claim 10
    wherein the computing device is further programmed to convert the background model to true vertical depth.
  14. 14. The system of
    claim 10
    wherein the assumed angle of inclination differs from the angle of inclination by less than 5°.
US09798867 1999-03-08 2001-03-01 Inhomogeneous background based focusing method for multiarray induction measurements in a deviated well Active US6442488B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US09264105 US6219619B1 (en) 1999-03-08 1999-03-08 Inhomogeneous background-based software focusing method for array-type induction logging tools
US09798867 US6442488B2 (en) 1999-03-08 2001-03-01 Inhomogeneous background based focusing method for multiarray induction measurements in a deviated well

Applications Claiming Priority (6)

Application Number Priority Date Filing Date Title
US09798867 US6442488B2 (en) 1999-03-08 2001-03-01 Inhomogeneous background based focusing method for multiarray induction measurements in a deviated well
CA 2439619 CA2439619C (en) 2001-03-01 2002-01-25 An inhomogeneous background based focusing method for multi-array induction measurements in a deviated well
EP20020703251 EP1373939B1 (en) 2001-03-01 2002-01-25 Focussing method for multi-array induction measurements in a deviated borehole
PCT/US2002/002346 WO2002071101A1 (en) 2001-03-01 2002-01-25 An inhomogeneous background based focusing method for multi-array induction measurements in a deviated well
GB0322001A GB2391321B (en) 2001-03-01 2002-01-25 An inhomogeneous background based focusing method for multi-array induction measurements in a deviated well
DE2002637276 DE60237276D1 (en) 2001-03-01 2002-01-25 Fokusierungsmethode for multi-array induction logs in an inclined hole

Publications (2)

Publication Number Publication Date
US20010039477A1 true true US20010039477A1 (en) 2001-11-08
US6442488B2 US6442488B2 (en) 2002-08-27

Family

ID=25174469

Family Applications (1)

Application Number Title Priority Date Filing Date
US09798867 Active US6442488B2 (en) 1999-03-08 2001-03-01 Inhomogeneous background based focusing method for multiarray induction measurements in a deviated well

Country Status (6)

Country Link
US (1) US6442488B2 (en)
EP (1) EP1373939B1 (en)
CA (1) CA2439619C (en)
DE (1) DE60237276D1 (en)
GB (1) GB2391321B (en)
WO (1) WO2002071101A1 (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030105591A1 (en) * 2001-12-03 2003-06-05 Teruhiko Hagiwara Method for determining anisotropic resistivity and dip angle in an earth formation
US6614716B2 (en) * 2000-12-19 2003-09-02 Schlumberger Technology Corporation Sonic well logging for characterizing earth formations
US6795774B2 (en) 2002-10-30 2004-09-21 Halliburton Energy Services, Inc. Method for asymptotic dipping correction
US20080109167A1 (en) * 2005-07-27 2008-05-08 Baker Hughes Incorporated Method of Generating a Deep Resistivity Image in LWD Measurements
US20100026305A1 (en) * 2008-07-31 2010-02-04 Baker Hughes Incorporated Method and Apparatus for Imaging Boreholes
WO2013009847A2 (en) * 2011-07-14 2013-01-17 Thrubit B.V. Smart drop-off tool and hang-off tool for a logging string
WO2014143015A1 (en) * 2013-03-15 2014-09-18 Halliburton Energy Services Inc. Dip correction for array induction tool data
US20150032376A1 (en) * 2012-02-16 2015-01-29 Halliburton Energy Services, Inc. Apparatus and methods of skin effect correction
WO2015099773A1 (en) * 2013-12-27 2015-07-02 Halliburton Energy Services, Inc. Dip correction in invaded zones

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6885947B2 (en) 2001-03-08 2005-04-26 Baker Hughes Incorporated Method for joint interpretation of multi-array induction and multi-component induction measurements with joint dip angle estimation
CA2500337A1 (en) * 2002-09-27 2004-04-08 Baker Hughes Incorporated A method for resistivity anisotropy determination in near vertical wells
US7382135B2 (en) * 2003-05-22 2008-06-03 Schlumberger Technology Corporation Directional electromagnetic wave resistivity apparatus and method
GB2426845B (en) * 2004-01-20 2008-08-06 Saudi Arabian Oil Co Real time earth model for collaborative geosteering
US7269514B2 (en) * 2004-05-11 2007-09-11 Halliburton Energy Services, Inc System and method for correcting induction logging device measurements by alternately estimating geometry and conductivity parameters
US8812237B2 (en) * 2009-02-05 2014-08-19 Schlumberger Technology Corporation Deep-reading electromagnetic data acquisition method

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5210691A (en) 1990-05-08 1993-05-11 Schlumberger Technology Corporation Method and apparatus for producing a more accurate resistivity log from data recorded by an induction sonde in a borehole
US5329448A (en) 1991-08-07 1994-07-12 Schlumberger Technology Corporation Method and apparatus for determining horizontal conductivity and vertical conductivity of earth formations
US5339037A (en) 1992-10-09 1994-08-16 Schlumberger Technology Corporation Apparatus and method for determining the resistivity of earth formations
US5841281A (en) 1996-03-22 1998-11-24 Western Atlas International, Inc. Method for determining electrical conductivity of earth formations from electromagnetic induction voltage measurements
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
US5905657A (en) 1996-12-19 1999-05-18 Schlumberger Technology Corporation Performing geoscience interpretation with simulated data
US6219619B1 (en) * 1999-03-08 2001-04-17 Baker Hughes Incorporated Inhomogeneous background-based software focusing method for array-type induction logging tools

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6614716B2 (en) * 2000-12-19 2003-09-02 Schlumberger Technology Corporation Sonic well logging for characterizing earth formations
US20030105591A1 (en) * 2001-12-03 2003-06-05 Teruhiko Hagiwara Method for determining anisotropic resistivity and dip angle in an earth formation
US6760666B2 (en) * 2001-12-03 2004-07-06 Shell Oil Company Method for determining anisotropic resistivity and dip angle in an earth formation
US6795774B2 (en) 2002-10-30 2004-09-21 Halliburton Energy Services, Inc. Method for asymptotic dipping correction
US20080109167A1 (en) * 2005-07-27 2008-05-08 Baker Hughes Incorporated Method of Generating a Deep Resistivity Image in LWD Measurements
US7894990B2 (en) * 2005-07-27 2011-02-22 Baker Hughes Incorporated Method of generating a deep resistivity image in LWD measurements
US20100026305A1 (en) * 2008-07-31 2010-02-04 Baker Hughes Incorporated Method and Apparatus for Imaging Boreholes
WO2013009847A2 (en) * 2011-07-14 2013-01-17 Thrubit B.V. Smart drop-off tool and hang-off tool for a logging string
WO2013009847A3 (en) * 2011-07-14 2013-04-04 Thrubit B.V. Smart drop-off tool and hang-off tool for a logging string
US8844618B2 (en) 2011-07-14 2014-09-30 Schlumberger Technology Corporation Smart drop-off tool and hang-off tool for a logging string
US20150032376A1 (en) * 2012-02-16 2015-01-29 Halliburton Energy Services, Inc. Apparatus and methods of skin effect correction
WO2014143015A1 (en) * 2013-03-15 2014-09-18 Halliburton Energy Services Inc. Dip correction for array induction tool data
WO2015099773A1 (en) * 2013-12-27 2015-07-02 Halliburton Energy Services, Inc. Dip correction in invaded zones
US9354353B2 (en) 2013-12-27 2016-05-31 Halliburton Energy Services, Inc. Dip correction in invaded zones

Also Published As

Publication number Publication date Type
EP1373939B1 (en) 2010-08-11 grant
GB2391321B (en) 2004-11-03 grant
DE60237276D1 (en) 2010-09-23 grant
EP1373939A1 (en) 2004-01-02 application
GB0322001D0 (en) 2003-10-22 grant
US6442488B2 (en) 2002-08-27 grant
CA2439619A1 (en) 2002-09-12 application
EP1373939A4 (en) 2009-07-01 application
CA2439619C (en) 2006-03-21 grant
WO2002071101A1 (en) 2002-09-12 application
GB2391321A (en) 2004-02-04 application

Similar Documents

Publication Publication Date Title
US7046009B2 (en) Method for measuring transient electromagnetic components to perform deep geosteering while drilling
US5038107A (en) Method and apparatus for making induction measurements through casing
US7031839B2 (en) Multi-frequency focusing for MWD resistivity tools
US6509738B1 (en) Electromagnetic induction well logging instrument having azimuthally sensitive response
US6100696A (en) Method and apparatus for directional measurement of subsurface electrical properties
US4278941A (en) High frequency induction log for determining resistivity and dielectric constant of the earth
US20050083063A1 (en) Electromagnetic method for determining dip angles independent of mud type and borehole environment
US5184079A (en) Method and apparatus for correcting data developed from a well tool disposed at a dip angle in a wellbore to eliminate the effects of the dip angle on the data
US6584408B2 (en) Subsurface formation parameters from tri-axial measurements
US6294917B1 (en) Electromagnetic induction method and apparatus for the measurement of the electrical resistivity of geologic formations surrounding boreholes cased with a conductive liner
US20050256644A1 (en) System and method for correcting induction logging device measurements by alternately estimating geometry and conductivity parameters
US20040059513A1 (en) Simultaneous resolution enhancement and dip correction of resistivity logs through nonlinear iterative deconvolution
US6541975B2 (en) Integrated borehole system for reservoir detection and monitoring
US6885947B2 (en) Method for joint interpretation of multi-array induction and multi-component induction measurements with joint dip angle estimation
US7382135B2 (en) Directional electromagnetic wave resistivity apparatus and method
US20050256645A1 (en) Determination of fracture orientation and length using multi-component and multi-array induction data
US6686736B2 (en) Combined characterization and inversion of reservoir parameters from nuclear, NMR and resistivity measurements
US5389881A (en) Well logging method and apparatus involving electromagnetic wave propagation providing variable depth of investigation by combining phase angle and amplitude attenuation
US6393363B1 (en) Method and apparatus for the measurement of the electrical resistivity of geologic formations employing modeling data
US6998844B2 (en) Propagation based electromagnetic measurement of anisotropy using transverse or tilted magnetic dipoles
US20030222651A1 (en) System and method for evaluation of thinly laminated earth formations
US7030617B2 (en) System, apparatus, and method for conducting electromagnetic induction surveys
US5854991A (en) Method for inversion processing of transverse electromagnetic induction well logging measurements
US20040019427A1 (en) Method for determining parameters of earth formations surrounding a well bore using neural network inversion
US6643589B2 (en) Simultaneous determination of formation angles and anisotropic resistivity using multi-component induction logging data

Legal Events

Date Code Title Description
AS Assignment

Owner name: BAKER HUGHES INCORPORATED, TEXAS

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:XIAO, JIAQI;GELDMACHER, INGO;REEL/FRAME:011880/0202;SIGNING DATES FROM 20010524 TO 20010530

FPAY Fee payment

Year of fee payment: 4

FPAY Fee payment

Year of fee payment: 8

FPAY Fee payment

Year of fee payment: 12