Source: https://patents.google.com/patent/US8085050B2/en
Timestamp: 2019-04-18 13:23:30+00:00

Document:
2007-05-21 Assigned to HALLIBURTON ENERGY SERVICES, INC reassignment HALLIBURTON ENERGY SERVICES, INC ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: HU, GUOYU (DAVID), BITTAR, MICHAEL S.
Methods and systems for determining the horizontal resistivity, vertical resistivity, and relative dip angle of anisotropic earth formations. Some of the disclosed methods and systems measure sinusoidal variation of azimuthally sensitive resistivity logging tool measurements, determine parameters representative of the sinusoidal variation, and perform inversion based on the sinusoidal parameters. When cast in this manner, the inversion process may yield more accurate and consistent resistivity and dip angle estimates. The sinusoidal parameters preferably take the form of average and peak-to-peak measurements, but may also take other forms. Moreover, use of such sinusoidal parameters enables a condensed representation of the resistivity logging tool measurements, enabling significantly more efficient communication and storage of these measurements. The condensed representations continue to enable directional boundary detection and geosteering.
In the field of petroleum well drilling and logging, resistivity logging tools are frequently used to provide an indication of the electrical resistivity of rock formations surrounding an earth borehole. Such information regarding resistivity is useful in ascertaining the presence or absence of hydrocarbons. A typical resistivity logging tool includes a transmitter antenna and a pair of receiver antennas located at different distances from the transmitter antenna along the axis of the tool. The transmitter antenna is used to create electromagnetic fields in the surrounding formation. In turn, the electromagnetic fields in the formation induce an electrical voltage in each receiver antenna. Due to geometric spreading and absorption by the surrounding earth formation, the induced voltages in the two receiving antennas have different phases and amplitudes. Experiments have shown that the phase difference (Φ) and amplitude ratio (attenuation, A) of the induced voltages in the receiver antennas are indicative of the resistivity of the formation. The formation region (as defined by a radial distance from the tool axis) to which such a resistivity measurement pertains is a function of the frequency of the transmitter and the distance from the transmitter to the mid-point between the two receivers. Thus, one may achieve multiple radial depths of investigation of resistivity either by providing multiple transmitters at different distances from the receiver pair or by operating a single transmitter at multiple frequencies.
If a formation is electrically isotropic, the resistivities measured at the various depths of investigation by such a resistivity logging tool will be the same. However, if the resistivities corresponding to the various depths of investigation are different, such differences indicate that the formation being measured is electrically anisotropic. In electrically anisotropic formations, the anisotropy is generally attributable to extremely fine layering during the sedimentary build-up of the formation. Hence, in a formation coordinate system oriented such that the x-y plane is parallel to the formation layers and the z axis is perpendicular to the formation layers, resistivities Rx and Ry in directions x and y, respectively, are the same, but resistivity Rz in the z direction is different from Rx and Ry. Thus, the resistivity in a direction parallel to the plane of the formation (i.e., the x-y plane) is known as the horizontal resistivity, Rh, and the resistivity in the direction perpendicular to the plane of the formation (i.e., the z direction) is known as the vertical resistivity, Rv. The index of anisotropy, η, is defined as η=[Rv/Rh]1/2.
The relative dip angle, θ, is the angle between the tool axis and the normal to the plane of the formation. If the axis of a resistivity logging tool is perpendicular to the plane of an anisotropic formation (i.e., θ=0°), both the phase shift and amplitude attenuation measurements reflect only the horizontal resistivity. However, if the axis of the tool is inclined with respect to the normal of the formation plane (i.e., for non-zero relative dip angle), the rock anisotropy affects the resistivity derived from phase shift measurements (“phase shift resistivity” or RΦ) differently than it affects the resistivity derived from amplitude attenuation measurements (“amplitude attenuation resistivity” or RA). For small relative dip angles (e.g., θ less than about 45°), the difference between phase shift and amplitude attenuation resistivities is relatively small. However, this difference becomes significant for relative dip angles greater than about 50°, and the difference is large for horizontal boreholes (i.e., θ=90°).
Thus anisotropy and dip each have significant effects on resistivity logging tool measurements. As a result, resistivity logging systems should account for formation anisotropy and relative dip if accurate resistivity logs are to be obtained. In U.S. Pat. No. 6,163,155, Michael Bittar discloses one such approach that employs resistivity logging tools having tilted transmitter and/or receiver antennas. Improvements to this disclosed approach have been developed and are presented herein.
FIG. 20 is an illustrative graph of the difference between the tool responses of FIGS. 18 and 19.
Disclosed herein are various methods and systems for determining the horizontal resistivity, vertical resistivity, and relative dip angle of anisotropic earth formations. Some of the disclosed methods and systems measure sinusoidal variation of azimuthally sensitive resistivity logging tool measurements, determine parameters representative of the sinusoidal variation, and perform inversion based on the sinusoidal parameters. When cast in this manner, the inversion process may yield more accurate and consistent resistivity and dip angle estimates. The sinusoidal parameters preferably take the form of average and peak-to-peak measurements, but may also take other forms. Moreover, use of such sinusoidal parameters enables a condensed representation of the resistivity logging tool measurements, enabling significantly more efficient communication and storage of these measurements. The condensed representations continue to enable directional boundary detection and geosteering.
The disclosed tool configurations and operations are best understood in the context of the larger systems in which they operate. Accordingly, an illustrative logging while drilling (LWD) environment is shown in FIG. 1. A drilling platform 102 supports a derrick 104 having a traveling block 106 for raising and lowering a drill string 108. A kelly 110 supports the drill string 8 as it is lowered through a rotary table 12. A drill bit 14 is driven by a downhole motor and/or rotation of the drill string 8. As bit 14 rotates, it creates a borehole 16 that passes through various formations 18. A pump 20 circulates drilling fluid through a feed pipe 22 to kelly 10, downhole through the interior of drill string 8, through orifices in drill bit 14, back to the surface via the annulus around drill string 8, and into a retention pit 24. The drilling fluid transports cuttings from the borehole into the pit 24 and aids in maintaining the borehole integrity.
The disclosed tool configurations and operations are best understood in the context of the larger systems in which they operate. Accordingly, an illustrative logging while drilling (LWD) environment is shown in FIG. 1. A drilling platform 102 is equipped with a derrick 104 that supports a hoist 106 for raising and lowering a drill string 108. The hoist 106 suspends a top drive 110 that is used to rotate the drill string 108 and to lower the drill string through the well head 112. Connected to the lower end of the drill string 108 is a drill bit 114. As bit 114 rotates, it creates a borehole 120 that passes through various formations 118. A pump 116 circulates drilling fluid through a supply pipe 118 to top drive 110, downhole through the interior of drill string 108, through orifices in drill bit 114, back to the surface via the annulus around drill string 108, and into a retention pit 124. The drilling fluid transports cuttings from the borehole into the pit 124 and aids in maintaining the integrity of the borehole 120.
In wells employing acoustic telemetry for LWD, downhole sensors (including resistivity logging tool 126) are coupled to an acoustic telemetry transmitter 128 that transmits telemetry signals in the form of acoustic vibrations in the tubing wall of drill string 108. An acoustic telemetry receiver array 130 may be coupled to tubing below the top drive 110 to receive transmitted telemetry signals. One or more repeater modules 132 may be optionally provided along the drill string to receive and retransmit the telemetry signals.
An electromagnetic resistivity logging tool 126 is integrated into the bottom-hole assembly near the bit 114. As the bit extends the borehole through the formations, logging tool 126 collects measurements relating to various formation properties as well as the tool orientation and position and various other drilling conditions. The logging tool 126 may take the form of a drill collar, i.e., a thick-walled tubular that provides weight and rigidity to aid the drilling process. A telemetry sub 128 may be included to transfer tool measurements to a surface receiver 130 and to receive commands from the surface receiver. In some alternative embodiments, the telemetry sub 128 collects and stores tool measurements for later retrieval when the tool is brought back to the surface.
The orientation measurements may be performed using an azimuthal orientation indicator, which may include magnetometers, inclinometers, and/or accelerometers, though other sensor types such as gyroscopes may be used. Preferably, the tool includes a 3-axis fluxgate magnetometer and a 3-axis accelerometer. As is known in the art, the combination of those two sensor systems enables the measurement of the toolface angle, borehole inclination angle, and borehole azimuth angle. In some embodiments, the toolface and hole inclination angles are calculated from the accelerometer sensor output. The magnetometer sensor outputs are used to calculate the hole azimuth. With the toolface, the hole inclination, and the hole azimuth information, various resistivity logging tools disclosed herein can be used to steer the bit to the desirable bed. Specifically, an azimuthal response difference or the response ratio can be used effectively to enter a desired payzone or to stay within the payzone of interest.
At various times during the drilling process, the drill string 108 is removed from the borehole as shown in FIG. 2. Once the drill string has been removed, logging operations can be conducted using a wireline logging tool 134, i.e., a sensing instrument sonde suspended by a cable 142 having conductors for transporting power to the tool and telemetry from the tool to the surface. A resistivity imaging portion of the logging tool 134 may have centralizing arms 136 that center the tool within the borehole as the tool is pulled uphole. A logging facility 144 collects measurements from the logging tool 134, and includes computing facilities for processing and storing the measurements gathered by the logging tool.
FIGS. 1 and 2 show formations 118 that are not perpendicular to the borehole, a situation that may occur naturally or may occur due to directional drilling operations. When measuring formation resistivity and orientation, it is convenient to use the formation coordinate system shown in FIG. 3. FIG. 3 shows a portion of a sedimentary formation bed, with the z-axis oriented perpendicular to the plane of the formation in the direction of the sedimentary accretion. As previously mentioned, the formation resistivity when measured along this axis is often different than formation resistivity measured in the x-y plane. In a dipping bed, the x-axis is chosen to be oriented in the direction of deepest descent, i.e., “downhill”.
The relationship between the formation coordinate system and the borehole coordinate system is shown in FIG. 4. The z-axis of the borehole coordinate system is aligned with the borehole's long axis, and the x-axis of the borehole coordinate system is directed to the north side (or alternatively, the high side) of the hole. As shown in FIG. 4, the two coordinate systems are related by two rotations. Beginning with the formation coordinate system (x,y,z), a first rotation of angle β is made about the z axis. The resulting coordinate system is denoted (x′,y′,z′). Angle β is the relative strike angle, which indicates the direction of the formation dip relative to the borehole's coordinate system. A second rotation of angle θ is then made about the y′ axis. This aligns the formation coordinate system with the borehole coordinate system. Angle θ is the relative dip angle, which is the slope angle of the beds relative to the long axis of the borehole.
The vertical resistivity is generally found to be the resistivity as measured perpendicular to the plane of the formation, and the horizontal resistivity is the resistivity as measured within the plane of the formation. Determination of each of these parameters (dip angle, strike angle, vertical resistivity, and horizontal resistivity) is desirable. Accordingly, disclosed herein is an improved downhole method and apparatus for simultaneously determining the horizontal resistivity, vertical resistivity, and relative dip angle for anisotropic earth formations. Some system embodiments employ an electromagnetic logging tool having an antenna configuration in which a transmitter antenna and a receiver antenna are oriented in non-parallel planes such that the vertical resistivity and the relative dip angle are decoupled. Preferably, either the transmitter or the receiver is mounted in a conventional orientation in a first plane that is normal to the tool axis, and the other antenna is mounted in a second plane that is not parallel to the first plane. The disclosed embodiments are suitable for LWD applications, and are also applicable to wireline and possibly other applications.
Also disclosed herein are an improved downhole method and apparatus for steering a downhole tool during directional drilling operations in order to maintain the borehole within a desired geological formation by providing an advance indication of the resistivity of a given bed before entry into that bed. In some embodiments, this steering capability is achieved by providing transmitter and receiver antennas that are mounted in non-parallel planes and computing the azimuthal difference or the ratio of the phase-based or amplitude-based responses of the receiver antennas. With such an antenna arrangement, the azimuthal ratio or difference of the responses indicates whether the resistivity of an approaching bed is higher or lower than the resistivity of the present bed. With such information, the driller may steer the drilling apparatus in order to maintain the borehole in a desired geological bed.
FIG. 5 shows a coordinate system for specifying the orientation of a loop antenna. The tool coordinate system (x″, y″, z″) has a z-axis oriented parallel to the long axis of the logging tool. The x-axis is directed to the “tool face scribe line” that is used to specify the tool orientation. In the tool coordinate system, the orientation of the loop antenna is represented by a normal vector N that is perpendicular to the plane of the loop antenna windings. The tilt angle of the antenna is the angle ξ between the z-axis and the normal vector N. The azimuthal orientation α of the loop antenna is the angle α between the x-axis and the projection of the normal vector N in the x-y plane.
FIG. 6 shows an illustrative logging tool 10 having N coaxial transmitter antennas T1, T2, T3, . . . , TN spaced along the length of the logging tool. Illustrative logging tool 10 also has two tilted receiver antennas R1 and R2 that are axially spaced apart from the transmitter antennas and from each other. Preferably, the transmitter and receiver antennas are mounted in recesses in tool 10 and are protected by a non-conductive material or a material with non-conducting apertures. In some tool embodiments designed to operate at more than one frequency, for example, f1=2 MHz and f2=1 MHz, each receiver antenna includes a pair of coils, with one coil tuned to f1 and one coil tuned to f2. In such embodiments, the coils in each pair may be located side by side around the periphery of tool 10 or may be concentrically stacked. The transmitter and receiver antennas may be fabricated in accordance with the teachings of U.S. Pat. No. 4,940,943. It should be appreciated that the body of tool 10 is preferably made of steel in order to prevent tool 10 from becoming a weak link in the drill string 14.
It should be appreciated that logging tool 10 also has the requisite electronic circuitry for processing the signals received by the receivers antennas as disclosed further below, thereby converting the received signals into a log or another indication of formation resistivity. It should also be appreciated that the processed signals can be recorded within the electronics section of tool 10 or may be coneyed to the surface by a telemetry system for concurrent processing and readout at the surface.
Still referring to the illustrative logging tool of FIG. 6, amplifiers A1, A2, A3 . . . AN are coupled to the transmitter coils T1, T2, T3 . . . TN, respectively. Each of the amplifiers is in turn driven by an oscillators F1, F2, F3 . . . FN. The operating frequencies of the oscillators are preferably between about 0.5 MHz up to about 4 MHz. Because of power attenuation at greater depths of investigation, such as is the case with the longer spaced transmitters, the frequencies preferably conform to the relationship F1≧F2≧F3≧ . . . FN. The oscillators F1, F2, F3 . . . FN are controlled by a transmitter enable circuit 30, which interfaces with a microprocessor 32, which in turn interfaces with a communication interface circuit 34 and an analog-to-digital (A/D) converter 36. Communication interface circuit 34 enables communication between the microprocessor 32 and a tool controller, integrated memory, a telemetry unit, and/or an external port for communicating with operators or computers once the tool has been removed from the borehole.
The receiver antennas R1 and R2 are respectively connected to amplifiers 40 and 42, which are connected, respectively, to mixer circuits 44 and 46. Oscillators F1, F2, F3 . . . FN are coupled to an oscillator select circuit 48, the output of which is connected to the inputs of mixer circuits 44 and 46. Oscillator select circuit 48 interfaces to microprocessor 32 to determine which oscillater is to be coupled to the mixer circuits.
The respective outputs of mixer circuits 44 and 46 drive low pass filters 50 and 52, respectively, the outputs of which drive amplitude measurement circuits 54 and 56, respectively. The outputs of amplitude measurement circuits 54 and 56 are connected to a multiplexer circuit 60. The outputs of low pass filter circuits 50 and 52 are also connected to the inputs of a relative phase measurement circuit 62, the output of which is fed into multiplexer 60. Under control of the microprocessor 32, multiplexer 60 forwards a selected one of its inputs to A/D converter 36 for sampling and conversion into a digital value that the microprocessor can store and process. The microprocessor 32 further acquires tool orientation measurements from orientation sensors 70, and determines a tool orientation to be associated with each resistivity measurement derived from the sampled receiver signals.
In the operation of the device and circuitry illustrated in FIG. 6, it should be appreciated that it is desirable to process intermediate frequency signals rather than the nominal signals received by receiver pair R1 and R2. Accordingly, the oscillators F1′, F2′, F3′ . . . FN′ are designed to provide frequencies very near the corresponding frequencies F1, F2, F3 . . . FN. For example, FN can be set at a frequency of 1.998 MHz and thus provide an intermediate frequency coming out of mixer circuits 44 and 46 of 0.002 MHz (2 KHz). In a similar fashion, F2′ and F3′ can be set at 1.998 MHz and 0.998 MHz, respectively. Thus, the only signals that pass to low pass filters 50 and 52 will be the intermediate frequencies which are obtained by mixing the frequencies of F1, F2, F3 . . . FN with the frequencies F1′, F2′, F3′ . . . FN′, respectively. It should be appreciated that amplitude measurement circuit 54 provides a measure of the amplitude of the signal received by receiver antenna R1, whereas amplitude measurement circuit 56 measures the amplitude of the incoming signals received by receiver antenna R2. Similarly, relative phase measurement circuit 62 provides an indication of the phase difference between the signals received at receiver antenna R1 and the signals received at receiver antenna R2. The amplitude measurements (A) and the phase measurements (Φ) (alternatively, the amplitude ratio and phase difference measurements) are each indicative of formation resistivity.
It should also be appreciated that the frequencies F1, F2, F3 . . . FN could all be the same frequency, though some embodiments employ different frequencies to account for increased power loss in the formation for larger transmitter-receiver antenna spacings. Preferably, the individual transmitter antennas fire in sequence, though in some alternative embodiments, simultaneous operation of the transmitter antennas is possible. However, those skilled in the art will recognize that simultaneous transmission of all of the transmitter signals will usually require additional filters and processing circuitry to enable the instrument to properly discriminate between the different frequencies.
In addition to the antenna configuration of FIG. 6, FIGS. 7A-7N illustrate a variety of contemplated transmitter/receiver antenna configurations, each having azimuthal sensitivity due to a tilted transmitter antenna, a tilted receiver antenna, or both. FIG. 7A shows tool having a coaxial transmitter antenna T1 and a tilted receiver antenna R1. FIG. 7B shows a tool having a tilted transmitter antenna T1 and a coaxial receiver antenna R1. However, the antenna configuration does not necessarily include a coaxial antenna. FIG. 7C shows a tool having a transmitter antenna T1 tilted at an angle ξT1 and a receiver antenna R1 tilted at an angle ξR1.
While the antenna configurations of FIGS. 7A-7C provide useful resistivity measurements, such measurements generally suffer from asymmetric responses at bed boundaries, often making the resistivity logs unnecessarily difficult to interpret. To address this issue, symmetric antenna configurations may be employed. FIG. 7D shows a tool having a tilted receiver antenna R1 positioned midway between two coaxial transmitter antennas T1, T2. FIG. 7E shows a tool having a coaxial receiver antenna R1 positioned midway between two tilted transmitter antennas T1, T2. The titlted transmitter antennas of FIG. 7E are parallel, but this is optional, as shown by the transmitter antennas in FIG. 7F, which are tilted at opposite azimuthal angles. As before, the antenna configuration does not necessarily include a coaxial antenna. FIG. 7G shows a tool having a first transmitter antenna T1 tilted at an angle ξT1, a second transmitter antenna T2 tilted at an angle ξT2, and a receiver antenna R1 tilted at an angle ξR1.
The antenna configurations of FIGS. 7A-7G are suitable for making absolute (non-differential) attenuation and phase shift measurements. While such measurements are useful, they may not have sufficient spatial resolution. To address this issue, differential antenna configurations may be employed. FIG. 7H shows a tool having a coaxial transmitter antenna T1 and a pair of tilted receiver antennas R1 and R2. FIG. 7I shows a tool having a tilted transmitter antenna T1 and a pair of coaxial receiver antennas R1 and R2. This antenna configuration does not necessarily include a coaxial antenna, and the pair of receiver antennas are not necessarily parallel. FIG. 7J shows a tool having a transmitter antenna T1 tilted at an angle ξT1, a first receiver antenna R1 tilted at an angle ξR1, and a second receiver antenna R2 tilted at an angle ξR2.
Though not suitable for providing improved spatial resolution, FIG. 7K shows an antenna configuration having a pair of co-located receiver antennas R1 and R2 that are tilted in different azimuthal directions. This antenna configuration may be particularly suitable for deriving a directional steering signal.
By combining symmetric and differential antenna configuration designs, it is possible to create antenna designs that provide compensated measurements, i.e., measurements that are protected against errors attributable to temperature-induced drift in the electronic circuit components. FIG. 7L shows a tool having a pair of tilted receiver antennas R1 and R2 having a midpoint centered between two coaxial transmitter antennas T1 and T2. FIG. 7M shows a tool having a pair of coaxial receiver antennas R1 and R2 having a midpoint centered between two tilted transmitter antennas T1 and T2. Note that the tilted antennas may be parallel, but are not necessarily so. Moreover, in each of the foregoing antenna configurations, the transmitter and receiver roles may be exchanged in accordance with the principle of reciprocity. FIG. 7N, for example, shows a pair of coaxial transmitter antennas T1 and T2 having a midpoint centered between two tilted receiver antennas R1 and R2. Furthermore, in each of the foregoing antenna configurations, additional transmitter and receiver antennas may be present to enable measurements at additional transmitter-receiver antenna spacings.
FIG. 8 illustrates the amplitude and phase measurements that may be made by a compensated resistivity logging tool 802. In response to a signal from transmitter antenna 808, receiver antennas 810 and 812 measure signals having an attenuation A1 and A2, respectively, and having phase shifts φ1 and φ2, respectively. From these measurements a first differential attenuation (logA2-logA1) and differential phase (φ2-φ1) can be determined. Similarly, in response to a signal from transmitter antenna 814, receiver antennas 810 and 812 measure signals having an attenuation A4 and A3, respectively, and having phase shifts φ4 and φ3, respectively. From these measurements a second differential attenuation (logA4-logA3) and second differential phase (φ4-φ3) can be determined. The first and second differential attenuation measurements can then be averaged together (as can the first and second differential phase measurements) to obtain a compensated measurement, i.e., a measurement in which fixed biases in the electronics are canceled.
FIG. 9 shows an illustrative relationship between attenuation and isotropic formation resistivity for a tilted antenna resistivity logging tool. FIG. 10 shows an illustrative relationship between phase shift and isotropic formation resistivity for a tilted antenna resistivity logging tool. Such relationships are known in the art, and will be employed in the description below.
To enable simultaneous measurement of resistivity, anisotropy, and dip, the disclosed logging tools and methods employ one or more tilted antennas to obtain azimuthally-sensitive resistivity measurements. The circumference of the tool (or borehole) is divided into azimuthal bins as shown in FIG. 11, and as the tool rotates and the tool face scribe line passes through each azimuthal bin, a set of resistivity measurements is made within each bin at a given position in the borehole. In various contemplated embodiments, the number of bins ranges as low as 6 or 8 bins, and as high as 16 or 32 bins. Depending on the relative rotational and axial traversal rates, multiple measurement sets may be obtained within each bin and combined (e.g., by averaging). The set of resistivity measurements may be absolute phase and/or attenuation measurements, or differential phase and/or attenuation measurements. However, in at least some embodiments, the resistivity measurements are compensated attenuation resistivity and compensated phase resistivity measurements (e.g., determined using relationships like those of FIGS. 9-10 and ignoring for the moment any dip or anisotropy effects), while in other embodiments the resistivity measurements are compensated attenuation resistivity measurements (e.g., using a relationship like that of FIG. 9). The set of resistivity measurements may include measurements at different frequencies and/or different transmitter-receiver spacings.
FIG. 12 illustrates the azimuthal dependence of phase resistivity measurements at different frequencies in a dipping anisotropic formation. FIG. 13 illustrates the azimuthal dependence of phase resistivity measurements at different transmitter-receiver spacings. At each frequency and spacing, the resistivity measurements are approximately sinusoidal. The sinusoidal-like nature of these measurements is preferably exploited to condense the azimuthally-dependent measurements into a representative set of characteristic parameters. In some preferred embodiments, the characteristic parameters are (1) a mean value, (2) a max-to-min difference value, and (3) an identifier for the location of the peak of the sinusoid (e.g., an azimuthal bin number). Other sets of representative characteristic parameters may also be employed. For example, the mean value could be replaced with a maximum value, a minimum value, or a root-mean-square value. The max-to-min difference value could be replaced with a max-to-mean difference value or with a second extremum value that, when paired with the first parameter, indicates the amplitude of the sinusoidal-like oscillation. The location of the peak could be replaced with the location of the valley, or dropped entirely if the strike angle of the formation is already known. In some embodiments it may be sufficient to determine a subset of a representative set of sinusoidal parameters to characterize a selected aspect of the azimuthal variation.
The condensed representation of resistivity measurements at a given depth in the borehole may be useful for storage and telemetry communications. More importantly, however, is the use of these characteristic parameters in the inversion process for determining vertical and horizontal resistivities and formation dip. When an inversion process employing only the maximum value or mean value is compared to an inversion process that employs a mean value and a max-to-min difference value, a significantly improved accuracy is observed and the determined formation characteristics are more resistant to measurement noise.
θ=the relative dip angle (the angle between the tool axis and the normal to the formation).
Equation  shows that the induced voltage, V, depends on kh and β. In turn, kh depends on σh; and β depends on σh, σv, and θ. These relationships indicate that σv and θ are dependent, and this dependency prevents convergence of a simultaneous solution for σh, σv, and θ. To break this dependency and enable a solution for σh, σv, and θ, it is desirable to have the transmitter antenna and receiver antenna tilted at different angles. Although the above formulation is for an coaxial transmitter with a tilted receiver, the theory of reciprocity provides that the same result also applies to a tilted transmitter with an untilted receiver. Indeed, both the transmitter and the receiver may be tilted, provided that the respective angles of tilt are not the same, i.e., ξT≠ξR.
The sinusoidal-like variation in FIGS. 12-13 is evident. FIGS. 14-16 illustrate the dip dependence of the sinusoid peak values for a typical earth formation having a horizontal resistivity of 1 ohm-m and a vertical resistivity of 4 ohm-m. FIG. 14 depicts amplitude attenuation resistivity and phase shift resistivity as a function of relative dip angle using a single transmitter-receiver pair at a single frequency. FIG. 15 depicts phase shift resistivity as a function of relative dip angle using three transmitter-receiver pair spacings at a single frequency. FIG. 16 depicts phase shift resistivity as a function of relative dip angle using a single transmitter-receiver pair at three different frequencies.
FIG. 17 shows an illustrative resistivity logging process to determine horizontal resistivity, vertical resistivity, and relative dip angle as a function of position along the borehole. After the resistivity logging tool has been placed in a borehole and drilling commences, with accompanying rotation of the logging tool, the logging process begins in block 1702 with the selection of a transmitter antenna and a transmit frequency. As the process continues, the selection operation iterates through each transmit antenna and each transmit frequency that is suitable for that transmit antenna.
In block 1704, the logging tool drives the selected transmit antenna at the selected frequency, and measures the receiver responses. The logging tool further measures the tool position and orientation to be associated with the measured receiver responses. In various embodiments, the receiver responses may be absolute phase shift, absolute attenuation, differential phase shift, and/or differential attenuation. In block 1706, the position and orientation information is used to associate the measured receiver responses with an azimuthal bin. If multiple measurements are obtained for a given bin, the measurements may be combined, e.g., by averaging. In block 1708, a test is made to determine if there are more measurements to be made at this position in the borehole (e.g., a test to see if the measured tool position is still within a predetermined range). If so, blocks 1702-1708 are repeated.
In optional block 1710, the logging tool determines compensated phase and/or attenuation measurements as previously described with respect to FIG. 8. In block 1712, the logging tool determines a set of sinusoidal parameters to represent the measurements obtained at the current position in the borehole. Preferably, the set of characteristic parameters is indicative of the mean and the max-to-min difference, and may also indicate the position of the peak or valley of the sinusoid. In some embodiments, the parameters are determined from a parameterized sine function that fits the measured data with a minimum square error. The sinusoidal parameters may be communicated to the surface for the subsequent operations to be performed by surface computing facilities.
In block 1714, the characteristic parameters are inverted to determine one or more of the formation strike angle, the formation dip angle, the formation resistivity (horizontal or vertical), and the formation anisotropy. The determined formation measurements are stored in the form of a log, which may be optionally displayed and updated in block 1716. In block 1718, a test is made to determine whether the logging should continue (e.g., whether the drilling process is ongoing), and if so, blocks 1702-1718 are repeated.
The inversion process can take the form of a look-up table-based search with interpolation, but is preferably performed with a closed-form forward model of the tool responses that are expected in response to estimates of the formation dip, resistivity, and anisotropy. Beginning with random, arbitrary, or in some cases, predetermined, estimates of the formation dip, resistivity, and anisotropy, the inversion process repeatedly updates these estimates until the sinusoidal parameter values predicted by the forward model match the measured sinusoidal parameter values to within some predetermined threshold. The update routine may, for example, employ the Levenberg-Marquardt method discussed by Tianfei Zhu and Larry D. Brown, “Two-dimensional Velocity Inversion and Synthetic Seismogram Computation,” Geophysics, vol. 52, no. 1, January 1987, p. 37-50.
In some embodiments, the measured characteristic parameter values are the mean and the max-to-min difference values of the differential phase shift measured as a function of azimuth at each of three or more transmit frequencies. In other embodiments, the characteristic parameter values are the mean and max-to-min difference values of the differential attenuation measured as a function of azimuth at three or more transmitter-to-receiver spacings. In yet other embodiments, both phase and attenuation information are measured and used as the basis of the inversion. Other characteristic parameter values are possible and may be used.
It is noted here that the tool response to a boundary between isotropic formations exhibits a sinusoidal response much like those shown in FIGS. 12 and 13. The maximum value of the response occurs in the azimuthal direction towards nearby boundaries with regions of higher conductivity and in the azimuthal direction away from nearby boundaries with regions of lower conductivity. Conversely, the minimum value of the response occurs in the azimuthal direction away from boundaries with regions of higher conductivity and in the azimuthal direction toward nearby boundaries with regions of lower conductivity. In both cases, the direction of the steering signal is consistent as the tool crosses a boundary.
In embodiments that condense the azimuthal resistivity measurement information into sinusoidal parameters such as a mean, a max-to-min difference, and peak orientation indicator, the max-to-min difference and the peak orientation indicator can be used as the steering signal. A max-to-mean difference and peak orientation indicator could be employed to similar effect. In some alternative embodiments, a ratio of max-to-min value is employed. The drilling operator may combine the steering signal information with tool position and orientation measurements and knowledge derived from test holes or seismic surveys to formulate directional steering decisions. The steering signals described herein are expected to be particularly effective at enabling a drilling operator to effectively detect and enter a payzone and at enabling the drilling operator to maintain a borehole course that maximizes the borehole interval in the payzone.
Although the foregoing specific details describe a preferred embodiment of this invention, persons reasonably skilled in the art of petroleum well drilling and logging will recognize that various changes may be made in the details of the method and apparatus of this invention without departing from the spirit and scope of the invention as defined in the appended claims. Therefore, it should be understood that this invention is not to be limited to the specific details shown and described herein.
a processor coupled to the at least one receiver antenna and the rotational angle sensor to collect at least one formation measurement for at least three azimuthal angles at each of the first and second frequencies, and to determine based on said measurements at least one sinusoidal parameter characterizing the azimuthal variation of said at least one formation measurement at the first frequency and at least one additional sinusoidal parameter characterizing the azimuthal variation of said at least one formation measurement at the second frequency.
2. The logging tool of claim 1, wherein the parameter is indicative of a max-to-min difference of the formation measurement.
3. The logging tool of claim 1, wherein the at least one formation measurement comprises at least one of a complex voltage, an amplitude ratio, and a phase difference derived from the received electromagnetic wave.
4. The logging tool of claim 3, wherein the formation measurement is determined with reference to the electromagnetic wave from the transmitter antenna.
5. The logging tool of claim 3, wherein the formation measurement is determined with reference to the electromagnetic wave received by a second receiver antenna.
6. The logging tool of claim 1, wherein the at least one formation measurement is a resistivity component of the formation.
7. The logging tool of claim 6, wherein the at least one formation measurement is vertical resistivity.
8. The logging tool of claim 1, wherein the processor calculates vertical and horizontal resistivity of the formation.
9. The logging tool of claim 8, wherein the processor further calculates a relative dip of the formation.
10. The logging tool of claim 1, wherein the transmitter antenna is tilted at a first angle relative to the tool axis, the receiver antenna is tilted at a second angle relative to the tool axis, and a difference between the first and second angles is less than 90°.
11. The logging tool of claim 1, further comprising a second receiver antenna oriented parallel to the first receiver antenna to receive an electromagnetic wave from the formation.
12. The logging tool of claim 1, wherein the second transmitter antenna is oriented parallel to the first transmitter antenna.
storing at least one sinusoidal parameter to characterize the azimuthal variation of the at least one formation measurement at the first frequency and at least one additional sinusoidal parameter to characterize the azimuthal variation of the at least one formation measurement at the second frequency.
14. The method of claim 13, wherein the formation measurement is at least one of phase shift, attenuation, and complex voltage derived from the received electromagnetic wave.
displaying a resistivity log based at least in part on at least one of the sinusoidal parameters.
16. The method of claim 15, wherein the resistivity log represents at least one of a horizontal resistivity of the formation, a vertical resistivity of the formation, and an anisotropy coefficient of the formation.
17. The method of claim 16, further comprising determining a relative dip angle of the formation based at least in part on at least one of the sinusoidal parameters.
a processor coupled to the receiver antenna and the rotational angle indicator to determine a formation measurement at each frequency for at least three rotational angles, and to determine for each frequency at least one parameter of an azimuthally-dependent sinusoid that best fits the formation measurement at that frequency.
19. The tool of claim 18, wherein the rotational angle indicator comprises at least one of a fluxgate, a magnetometer, an accelerometer, and a gyroscope.
20. The tool of claim 18, wherein the at least one parameter comprises an indicator of an azimuthal orientation of a peak of the periodical signal.
MX343007B (en) * 2012-06-29 2016-10-21 Halliburton Energy Services Inc Full tensor micro-impedance imaging.
"Australian First Examiner's Report", dated Mar. 5, 2010, Appl No. 2007349251, "Robust Inversion Systems and Methods for Azimuthally Sensitive Resistivity Logging Tools" filed Oct. 3, 2008.
"EP Non-Final Office Action", dated Dec. 30, 2009, Appl No. 01 962 294.03, Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Geosteering Within a Desired Payzone, filed Jul. 1, 2001, 3 pgs.
"European Office Action", dated Jul. 17, 2006, Application No. 01 096 294.3 , Electromagnetic wave resistivity tool having a tilted antenna for determining the horizontal and vertical resistivities and relative dip angle in anisotropic earth formations, Jul. 10, 2001.
"European Office Action", dated Jul. 31, 2007, Application No. 00908351.0, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Relative Dip Angle in Anisotropic Earth Formations" filed Jan. 24, 2000.
"European Office Action", dated Sep. 23, 3008, Application No. 01096 294.3, Electromagnetic wave resistivity tool having a tilted antenna for determining the horizontal and vertical resistivities and relative dip angle in anisotropic earth formations Jul. 10, 2001.
"European Office Action", dated Sep. 3, 2010, Appl No. 09166567.9, Electromagnetic Wave Resistivity Tool Having a Tilted Antenna, filed Jul. 28, 2009, 4 pgs.
"European Search Report", dated Dec. 14, 2009, Appl No. 09166566.1, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Relative Dip Angle in Anisotropic Earth Formations", filed Jul. 28, 2009, 9 pgs.
"European Search Report", dated Nov. 17, 2009, Appl No. 09166404.5, "Electromagnetic Wave Resistivity tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Relative Dip Angle in Anisotropic Earth Formations", filed Jul. 24, 2009, 4 pgs.
"European Search Report", dated Nov. 20, 2009, Appl No. 09166405.2, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Related Dip Angle in Anisotropic Earth Formations", filed Jul. 24, 2009, 6 pgs.
"European Search Report", dated Oct. 16, 2009, Appl No. 09166567.9, Electromagnetic Wave Resistivity Tool Having a Tilted Antenna, filed Jul. 28, 2009, 5 pgs.
"European Search Report", dated Sep. 25, 2009, Appl No. 09166566.1, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Relative Dip Angle in Anisotropic Earth Formations", filed Jul. 28, 2009, 7 pgs.
"European Supplemental Search Report", dated Jun. 12, 2003 Application No. 00908351.0, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Relative Dip Angle in Anisotropic Earth Formations", filed Jan. 24, 2000.
"Indonisia First Office Action", dated Sep. 2009, Appl No. WOO200803740, Method and Apparatus for Building a Tilted Antenna, filed Nov. 19, 2008, 2 pgs.
"Malaysian Office Action", dated Nov. 30, 2010, Appl No. 032599, "Robust Inversion Systems and Methods for Azimuthally Sensitive Resistivity Logging Tools" filed Oct. 15, 2008, 2 pgs.
"Malaysian Search Report", dated Nov. 30, 2010, Appl No. 032599, "Robust Inversion Systems and Methods for Azimuthally Sensitive Resistivity Logging Tools" filed Oct. 15, 2008, 1 pg.
"No Office Action", dated Apr. 3, 2009, Application No. 2001 3707, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Relative Dip Angle in Anisotropic Earth Formations", filed Jan 24, 2000, 5 pgs.
"PCT International Preliminary Report on Patentability", dated Apr. 30, 2008, Appl No. PCT/US06/62149, "Antenna Coupling Component Measurement Tool Having a Rotating Antenna Configuration", filed Dec. 15, 2006, 6 pgs.
"PCT International Preliminary Report on Patentability", dated Oct. 1, 2009, Appl No. PCT/US07/064221, "Robust Inversion Systems and Methods for Azimuthally Sensitive Resistivity Logging Tools", filed Mar. 16, 2007, 2 pgs.
"PCT International Search Report and Written Opinion", dated Jun. 27, 2008, Appl No. PCT/US08/51447, EM-Guided Drilling Relative to an Existing Borehole, 8 pgs.
"PCT International Search Report and Written Opinion", dated May 15, 2000, Application No. PCT/US00/01693, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Relative Dip Angle in Anisotropic Earth Formations", 01/24/200.
"PCT Written Opinion of the International Searching Authority", dated Feb. 5, 2008, Appl No. PCT/US07/64221, "Robust Inversion systems and Methods for Azimuthally Sensitive Resistivity Logging Tools", filed Mar. 16, 2007, 4 pgs.
"Preliminary Report on Patentability", dated Jul. 29, 2010, Appl No. PCT/US08/51447, "EM-Guided Drilling Relative to an Existing Borehole", filed Jan. 18, 2008, 7 pgs.
"Russian Office Action", dated Jul. 9, 2009, Appl No. 2009104466, Modular Geosteering Tool Assembly, filed Feb. 10, 2009, 8 pgs.
"Singapore Examination Report", dated Oct. 28, 2010, Application No. 200808573-0, "Method and Apparatus for Building a Tilted Antenna", filed Nov. 19, 2008, 6 pgs.
"Singapore First Office Action", dated Apr. 11, 2009, Appl. No. 200808573-0, "Method and Apparatus for Building a Tilted Antenna", Filing Date Nov. 19, 2008, 8 pgs.
"US Final Office Action", dated Oct. 15, 2010, U.S. Appl. No. 12/467,427, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Relative Dip Angle in Anisotropic Earth Formations", May 18, 2009, 24 pgs.
"US Final Office Action", dated Oct. 22, 2010, U.S. Appl. No. 12/467,434, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Geosteering Within a Desired Payzone", filed May 18, 2009, 19 pgs.
"US Non-Final Office Action", dated Apr. 26, 2000, Application No. 09/23832, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Relative Dip Angle in Anisotropic Earth Formations", filed Jan. 28, 1999.
"US Non-Final Office Action", dated Dec. 23, 2010, U.S. Appl. No. 12/306,267, "Modular Geosteering Tool Assembly", filed Dec. 23, 2008, 18 pgs.
"US Non-Final Office Action", dated Jul. 28, 2003, Application No. 10/255,048, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Relative Dip Angle in Anisotropic Earth Formations", filed Sep. 25, 2002.
"US Non-Final Office Action", dated Jun. 3, 2010, U.S. Appl. No. 12/467,434, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Geosteering Within a Desired Payzone", filed May 18, 2009, 19 pgs.
"US Non-Final Office Action", dated Mar. 15, 2010, U.S. Appl. No. 12/467,427, "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Determining the Horizontal and Vertical Resistivities and Relative Dip Angle in Anisotropic Earth Formations", filed May 18, 2009, 25 pgs.
Bittar, Michael S., "Electromagnetic Wave Resistivity Tool Having a Tilted Antenna for Geosteering Within a Desired Payzone", U.S. Appl. No. 11/457,709, filed Jul. 14, 2006.
Bittar, Michale S., PCT Application US2007/075455, filed Aug. 8, 2007, Resistivity Logging with Reduced Dip Artifacts.
Gianzero, Stanley et al., "Method and Apparatus Having Antennas Configured to Measure Electrical Anisotropy", U.S. Appl. No. 10/957,415, filed Oct. 1, 2004.
Li, Qiming et al., "New Directional Electromagnetic Tool for Proactive Geosteering and Accurate Formation Evaluation While Drilling", SPWLA 46th Annual Logging Symposium, New Orleans, LA, USA, (Jun. 26, 2005), pp. 1-16.
UK Examination Report, dated Aug. 19, 2011, Appl No. 0914189.6, "EM-Guided Drilling Relative to an Existing Borehole", filed Aug. 13, 2009, 5 pgs.
US Non-Final Office Action, dated Jul. 18, 2011, U.S. Appl. No. 12/306,954, "Method and Apparatus for Building a Tilted Antenna", filed Dec. 30, 2008, 21 pgs.

References: Application No. 01
 Application No. 00908351
 Application No. 01096
 Application No. 00908351
 Application No. 2001
 Application No. 200808573
 Application No. 09
 Application No. 10