Well logging apparatus and method with borehole compensation

Abstract: A well logging apparatus and technique for measuring formation resistivity is disclosed which provides advantages of borehole compensation without the requirement of having, for each transmitter to receiver pair spacing, an "equal and opposite" transmitter with the same spacing on the other side of the receiver pair. In one disclosed embodiment there are provided two transmitting antennas (T 1 ,T 2 ) at equal and opposite spacing from the midpoint of the receiving antenna pair (R 1 ,R 2 ), and a third transmitting antenna (T 3 ) at a different spacing from the midpoint of the receiving antenna pair (R 1 ,R 2 ), without the need for providing a fourth transmitting antenna to obtain the benefits of borehole compensation. This reduces the complexity and cost of adding additional depths of investigation to a logging tool, while retaining the advantages of borehole compensation. In another disclosed embodiment at least three transmitting antennas are provided with at least one of them (T 2 ) on one side of the receiving antenna pair (R 1 ,R 2 ) and at least two of them (T 1 ,T 3 ) on the other side of the receiving antenna pair (R 1 ,R 2 ). In this embodiment, none of the spacings between respective transmitting antennas and the midpoint of the receiving antenna pair are the same.

Description: BACKGROUND OF THE INVENTION FIELD OF THE INVENTION This invention relates to the field of well logging and, more particularly, to well logging techniques and apparatus for determining formation properties, such as resistivity, at several different radial depths of investigation and with compensation for factors such as borehole rugosity and drift of electronic components. The invention has general application in the well logging art, but is particularly useful in measuring while drilling. DESCRIPTION OF THE RELATED ART A commonly used technique for evaluating formations surrounding an earth borehole is resistivity logging. Porous formations having high resistivity generally indicate the presence of hydrocarbons, while porous formations with low resistivity are generally water saturated. However, the region immediately surrounding the borehole can be invaded by borehole fluid or mud filtrate and have a different resistivity than the virgin formation. If a resistivity logging device has only one radial depth of investigation, there is limited ability to measure resistivity of all zones of interest, and there may be difficulty in determining if the measured resistivity represents the invaded zone, the virgin zone, or some combination of the two zones. However, if the resistivity logging device has multiple radial depths of investigation, there is greater flexibility. In addition to the advantage of having, for example, a shallow measurement and a deep measurement individually, the combination of the two can provide additional information such as the extent of invasion. It is also possible to combine two or more measurements, for example a shallow measurement and a deeper measurement, to compute a better estimate of the true formation resistivity. Wireline resistivity logging tools have long been provided with two or more radial depths of investigation. Conventional wireline resistivity logging tools typically achieve two depths of investigation by using a short and a long vertical array of electrodes or coils. In general, a long vertical array provides a greater radial depth of investigation than does a short vertical array. More recently, as will be discussed momentarily, measuring while drilling tools have been provided with multiple radial depths of investigation. A type of well logging which is of interest herein is so-called electromagnetic propagation logging, which can be used to measure the resistivity of the formation surrounding a borehole. For example, U.S. Patent No. 3,551,797 describes a technique wherein electromagnetic energy is transmitted into the formation, and energy which returns to the borehole is measured at a receiver pair to determine the attenuation and/or the phase shift of the electromagnetic energy propagating in the formation. More than one vertical spacing between a transmitter and different receiver pairs may be used to obtain different radial depths of investigation. For example, a receiver pair relatively close to the transmitter can be used to obtain attenuation and/or phase shift information from which the properties of the invaded zone are determined, and measurements of the attenuation and/or phase shift from a receiver pair relatively far from the transmitter can be used to obtain the properties of the deeper uninvaded formations. Either attenuation or phase shift can be used to determine a bulk value of the electromagnetic skin depth for the formation, with the bulk conductivity then being determinable from the electromagnetic skin depth. Various other techniques also exist in the art for utilizing multiple transmitters and/or receivers to investigate resistivity at different depths of investigation. In U.S. Patent No. 4,899,112 there is disclosed a logging apparatus for determining resistivity of formations at two different radial depths of investigation using signals received at only a single receiver pair. The resistivity of formations at a relatively shallow depth of investigation around the receiver pair is determined as a function of the phase shift measured at the receiver pair, and the resistivity of formations at a relatively deep depth of investigation around the receiver pair is determined as a function of the attenuation measured at the receiver pair. The apparatus is particularly advantageous for measurement while-drilling, where it is desirable to obtain resistivity at multiple depths of investigation while minimizing the length and complexity of the logging device. u.s. Patent No. 4,899,112 also discloses a so-called borehole compensated embodiment, where the receiver pair is located between and equally spaced from, a pair of transmitting antennas that can be alternately energized. The signals received at the receiver pair can be averaged to obtain borehole compensated signals; i.e. signals from which (1) drift of electronic components, and (2) borehole rugosity, have been reduced or removed by cancellation. One or more additional transmitters can be added to the logging device just described, at different spacing(s) , to attain further depths of investigation. However, in order to provide borehole compensation for each additional spacing, it would be necessary to provide two transmitters (one on each side of the receivers) for each such additional spacing. U.S. Patent No. 5,187,661 describes a logging apparatus equipped with multiple transmitter an receiver antennas for determining various formation properties. A sequence of electromagnetic pulses is transmitted at multiple frequencies to make the subsurface measurements with this tool. The disclosed apparatus is configured with its transmitter antennas equally spaced above and below a designated mid-point. Mack S. G. et al., MWD tool accurately measures four resistivities , OIL AND GAS JOURNAL, vol. 90, no. 21, 25 May 1992, pages 42-46 describes a measurement while drilling sensor equipped with multiple transmitter and receiver antennas for determining logs at several depths of investigation. The disclosed apparatus is configured with a pair of receiver antennas and a series of transmitter antennas disposed at one end of the receiver pair. The transmitter antennas are disposed along the tool at varied spacings from the receiver pair. It is among the objects of the present invention to provide a logging device and technique that determines formation properties at several depths of investigation, with the advantage of borehole compensation and also with minimal use of antennas. SUMMARY OF THE INVENTION The apparatus and technique of the present invention provides advantages of borehole compensation without the requirement of having, for each transmitter to receiver pair spacing, an "equal and opposite" transmitter with the same spacing on the other side of the receiver pair. In one form of the invention, there is provided an apparatus for determining a property of earth formations surrounding a borehole. The apparatus comprises an elongated logging device movable through the borehole; a transmitter array of first, third, and second transmitting antennas longitudinally spaced on the device; a pair of longitudinally spaced receiving antennas on the device, the pair of receiving antennas being located between the second transmitting antenna and one of the first and third transmitting antennas; the first, third, and second transmitting antennas being spaced from the midpoint between the receiver pair by first, third, and second distances, respectively, and the first and second distances being substantially the same and being different from the third distance; means for energising the first, third, and second transmitting antennas with electromagnetic energy; means for receiving, at the pair of receiving antennas, electromagnetic energy from the transmitting antennas; means for determining, from electromagnetic energy received at the pair of receiving antennas, first, third, and second phase shift and/or attenuation measurements of electromagnetic energy transmitted from the first, third, and second transmitting antennas, respectively; and means for determining a property of said earth formations from the first, third, and second phase shift and/or attenuation measurements. In another form of the invention, an apparatus for determining a property of earth formations surrounding a borehole is provided. The apparatus comprising an elongated logging device movable through the borehole; the device including a sequence of spaced apart transmitting antennas with odd numbered designations 1, 3,. .. (2m-1) and at least one spaced apart transmitting antenna with an even numbered designation 2 ... (2m), where m is a positive integer; the device also including a pair of longitudinally spaced receiving antennas, said pair of receiving antennas being located between said sequence of transmitting antennas with odd numbered designations and said at least one antenna with an even numbered designation; the spacing of at least one of said transmitting antennas from the midpoint between said pair of receiving antennas being different than the spacing of every other transmitting antenna from the midpoint between said pair of receiving antennas; means for energising each of said transmitting antennas with electromagnetic energy; means for receiving, at said pair of receiving antennas, electromagnetic energy from said transmitting antennas; means for determining, from electromagnetic energy received at said pair of receiving antennas, phase shift and/or attenuation measurements of electromagnetic energy transmitted from each of the respective transmitting antennas; and means for determining a property of said earth formations from said phase shift and/or attenuation measurements. Another form of the invention provides a method for determining a property of earth formations surrounding a borehole. The method comprises providing an elongated logging device movable through the borehole, said device having thereon a sequence of spaced apart transmitting antennas with odd numbered designations 1, 3 ... (2m-1), and at least one spaced apart transmitting antenna with an even numbered designation 2 ... (2m), where m is a positive integer, and a pair of longitudinally spaced receiving antennas, said pair of receiving antennas being located between said sequence of transmitting antennas with odd numbered designations and said at least one transmitting antenna with an even numbered designation, the spacing of at least one of said transmitting antennas from the midpoint between said pair of receiving antennas being different than the spacing of every other transmitting antenna from the midpoint between said pair of receiving antennas; energising each of said transmitting antennas with electromagnetic energy; receiving, at said pair of receiving antennas, electromagnetic energy from said transmitting antennas; determining, from electromagnetic energy received at said pair of receiving antennas, phase shift and/or attenuation measurements of electromagnetic energy transmitted from each of the respective transmitting antennas; and determining a property of said earth formations from said phase shift and/or attenuation measurements. Further features and advantages of the invention will become more readily apparent from the following detailed description when taken in conjunction with the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS Figure 1 is a schematic diagram, partially in block form, of a system in accordance with an embodiment of the invention. Figure 2 is a schematic diagram, partially in block form, of an embodiment of the electronics of the Fig. 1 embodiment. Figure 3 illustrates a variation of the Figure 2 embodiment. Figure 4, which includes Figures 4A and 4B placed one below another, is a flow diagram of a routine for controlling a processor in practicing the embodiment of Figure 2 or 3. Figure 5 illustrates a logging device with five transmitting antennas in accordance with a further embodiment of the invention. Figure 6 is a diagram which illustrates the placement of antennas in accordance with a further embodiment of the invention. Figure 7 illustrates another logging device with three transmitters in accordance with a further embodiment of the invention. Figure 8, which includes Figures 8A and 8B placed one belowanother, is a flow diagram of a routine for controlling a processor in accordance with a further embodiment of the invention. DESCRIPTION OF THE PREFERRED EMBODIMENTS Referring to Fig. 1, there is illustrated an embodiment of the invention in the form of a logging-while-drilling apparatus and method. A platform and derrick 10 are positioned over a borehole 11 that is formed in the earth by rotary drilling. A drill string 12 is suspended within the borehole and includes a drill bit 15 at its lower end. The drill string 12, and the drill bit 15 attached thereto, is rotated by a rotating table 16 (energized by means not shown) which engages a kelly 17 at the upper end of the drill string. The drill string is suspended from a hook 18 attached to a travelling block (not shown). The kelly is connected to the hook through a rotary swivel 19 which permits rotation of the drill string relative to the hook. Drilling fluid or mud 26 is contained in a pit 27 in the earth. A pump 29 pumps the drilling fluid into the drill string 12 via a port in the swivel 19 to flow downward through the center of drill string 12. The drilling fluid exits the drill string via ports in the drill bit 15 and then circulates upward in the region between the outside of the drill string and the periphery of the borehole. As is well known, the drilling fluid thereby carries formation cuttings to the surface of the earth, and the drilling fluid is returned to the pit 27 for recirculation. The small arrows in the Figure illustrate the typical direction of flow of the drilling fluid. Mounted within the drill string 12, preferably near the drill bit 15, is a downhole sensing, processing, storing and transmitting subsystem 100. Subsystem 100 includes a measuring apparatus 200 which includes at least three transmitting antennas T 1 , T 2 , and T 3 , and receiving antennas R 1 and R 2 , and operates in the manner to be described below. [The antennas can be of the type described in the referenced U.S. Patent 4,899,112; that is coils wound on mounting material around a special section of metal drill collar which comprises part of the subsystem 100.] A communications transmitting portion of the downhole subsystem includes an acoustic transmitter 56, which generates an acoustic signal in the drilling fluid that is representative of the measured downhole conditions. One suitable type of acoustic transmitter, which is known in the art, employs a device known as a "mud siren" which includes a slotted stator and a slotted rotor that rotates and repeatedly interrupts the flow of drilling fluid to establish a desired acoustic wave signal in the drilling fluid. The generated acoustic mud wave travels upward in the fluid through the center of the drill string at the speed of sound in the fluid. The acoustic wave is received at the surface of the earth by transducers represented by reference numeral 31. The transducers, which are, for example, piezoelectric transducers, convert the received acoustic signals to electronic signals. The output of the transducers 31 is coupled to the uphole receiver subsystem 90 which is operative to demodulate the transmitted signals, which are then coupled to processor 85 and recorder 45. Transmitter 56 is controlled by transmitter control and driving electronics 57 which includes analog-to-digital (A/D) circuitry that converts the signals representative of downhole conditions into digital form. The control and driving electronics 57 may also include a suitable modulator, such as a phase shift keying (PSK) modulator, which conventionally produces driving signals for application to the transmitter 56. These driving signals can be used to apply appropriate modulation to the mud siren of transmitter 56. It will be understood that alternative techniques can be employed for communicating logging information to the surface of the earth. The downhole subsystem 100 further includes acquisition and processor electronics 58. These electronics include a microprocessor (with associated memory, clock circuitry, and interface circuitry), and processing circuitry. The acquisition and processor electronics 58 are coupled to the measuring apparatus 200 and obtain measurement information therefrom. The acquisition and processor electronics is capable of storing data from the measuring apparatus, processing the data and storing the results, and coupling any desired portion of the information it contains to the transmitter control and driving electronics 57 for transmission to the surface by transmitter 56. A battery 53 may provide downhole power. As known in the art, a downhole generator (not shown) such as a so-called « mud turbine » powered by the drilling fluid, can also be utilized to provide power during drilling. Fig. 2 is a block diagram of electronics in the measuring apparatus 200 and the acquisition and processor electronics 58. An oscillator 211 produces an electrical signal of frequency f [e.g. at 2 MHz], which is amplified by amplifier 212 and applied, via electronically controlled switch 275, to one of the transmitting antennas T 1 , T 2 , or T 3 , through respective impedance matching circuits 221, 222 and 223. The switch 275, under control of the processor 270, selects the transmitter T 1 , T 2 , or T 3 , that is to be energized. The signal from receiver R 1 is coupled, via impedance matching circuit 230 and preamplifier 231, to one input of an electronically controlled switch 250. Similarly, the signal from receiver R 2 is coupled via impedance matching circuit 240 and preamplifier 241 to the other input of the switch 250. The switch 250, under control of processor 270, selects the (R 1 ) or (R 2 ) receiver output. The selected signal is amplified by amplified 251 and subsequently translated to a lower frequency Δf using a known heterodyne technique. A local oscillator 213 is in a phase locked loop (represented at 215) with the master oscillator 211. The local oscillator 213 has a frequency f - Δf, where Δf is typically a few kilohertz or less. A signal from the local oscillator is mixed with the received signal by a mixer 255, and the mixer output is passed through a low pass filter 257 which blocks signals at f and f + Δf and passes the signal at Δf. The signal at frequency Δf contains the phase and amplitude information of the original signal at frequency f. The recognized advantage of this frequency translation is that it is easier to measure the phase and amplitude at kilohertz or lower frequencies than at megahertz frequencies. The signal at frequency Δf is measured with a phase meter 261 and with an amplitude meter 262, and the results are input to the processor 270. The phase meter 261 may utilize a reference signal from the phase locked loop 215. The phase and amplitude meters may also include sample-and-hold circuits, to implement comparison of the signal from the respective receivers. Alternatively, the processor may compute the relative phase and amplitude from the measurements it receives. Reference can be made to U.S. Patents 4,185,238 and 4,899,112. As noted above, the processor 270 has associated memory clocking, and interface circuitry (not shown), as is conventional. The processor implements storage of the measurements of phase and amplitude, processing of these measurements in the manner to be described, storage of the results of the processing, and/or coupling of the measurements and/or processing results to the control and driving electronics for transmission to the earth's surface. In the present embodiment, a downhole clock is utilized to keep track of time, which can subsequently be correlated with depth level by virtue of keeping a record of the drill string progress. The clock, which can typically be part of the system in conjunction with the processor 270, can be synchronized with the system before a trip downhole. Also, if desired, communication with the downhole subsystem, using mud pulse technology or other suitable communication means, can be utilized to convey timing synchronization and/or depth level information, as well as to communicate data in general. It will be understood that the invention can be utilized in conjunction with any suitable technique for keeping track of depth level. With the exception of transmitter T 3 and its associated circuitry (and the processing of the signals received therefrom), the Fig. 2 circuit is like that of the borehole compensated embodiment of U.S. Patent 4,899,112. U.S. Patent 4,899,112 describes in detail how resistivity at a relatively deep depth of investigation is obtained from attenuation measured at the receiver pair, and how resistivity at a relatively shallow depth of investigation is obtained from phase shift measured at the receiver pair. As an example, a table look-up can be used for this purpose, and reference can be made to said U.S. Patent 4,899,112 for further detail. Consider, first, the operation of borehole compensation in the referenced U.S. Patent 4,899,112. When T 1 is activated (and T 2 off), a downward propagating wave has a phase P 1D and an amplitude A 1D at the location of receiver R 1 and a phase P 2D and an amplitude A 2D at the location of receiver R 2 . The true phase shift for the downward propagating wave is PS=P 2D -P 1D , and the true attenuation is AT=20.log 10 (A 2D /A 1D ), where P 1D , P 2D , A 1D , and A 2D are the phases and amplitudes as measured by an ideal circuit with no phase shift, gain, or offset errors. The actual or measured voltage derived from R 1 and input to switch 250 has a net phase X 1D =P 1D +P 1 , where P 1 is the additional phase increment due to receiver R 1 , impedance matching circuit 230, and preamplifier 231. The actual or measured voltage derived from R 1 and input to switch 250 has an amplitude B 1D =A 1 *A 1D , where A 1 is the gain due to receiver R 1 , impedance matching circuit 230, and preamplifier 231. Similarly, the voltage derived from R 2 and input to switch 250 has a net phase X 2D =P 2D +P 2 and an amplitude B 2D =A 2 *A 2D , where P 2 and A 2 are related to receiver R 2 , impedance matching circuit 240, and preamplifier 241. The rest of the electronics chain (switch 250, amplifier 251, mixer 255, low pass filter 257, phase meter 261, amplitude meter 262, and processor 270) are common to signals from R 1 and R 2 , and do not introduce phase or amplitude differences. The phase shift of the downward propagating wave, as measured by the electronics, can be expressed as PS DM =X 2D -X 1D = (P 2D -P 1D ) + (P 2 -P 1 ) =PS + (P 2 -P 1 ) PS+PS ERR and the attenuation of the downward propagating wave as measured by the electronics can be expressed as: AT DM =20 . log 10 (B 2D /B 1D )=20 . log 10 (A 2D /A 1D )+20 . log 10 (A 2 /A 1 ) =AT+20 . log 10 (A 2 /A 1 )=AT+AT ERR The differences between the two receivers and between the two receiver electronics introduce a phase shift error (PS ERR =P 2 -P 1 ) and an attenuation error (AT ERR =20.log 10 (A 2 /A 1 )), which reduce the accuracy of the logging device measurements, particularly for relatively high formation resistivities where variations of PS and AT with respect to resistivity are small. These errors can change with temperature, pressure, borehole mud resistivity, and time. The inaccuracies which occur in uncompensated measurements are eliminated in the referenced U.S. Patent 4,899,112 with borehole compensation, which is obtained using the second transmitter T 2 . This can be seen by considering the case with T 2 activated and T 1 off. The upward propagating wave has phase P 1U and amplitude A 1U at the location of receiver R 1 and a phase P 2U and an amplitude A 2U at the location of receiver R 2 . The true phase shift for the upward propagating wave is PS=P 1U -P 2U , and the true attenuation is AT=20.log 10 (A 1U /A 2U ). These quantities will generally be the same for upward and downward propagating waves. The actual or measured voltage derived from R 1 and input to switch 250 has a net phase X 1U =P 1U +P 1 and an amplitude B 1U =A 1 *A 1U . The voltage derived from R 2 and input to switch 250 has a net phase X 2U =P 2U +P 2 and an amplitude B 2U =A 2 *A 2U . The phase shift of the upward propagating wave as measured by the electronics can be expressed as PS UM = (X 1U -X 2U ) = (P 1U -P 2U ) + (P 1 -P 2 ) =PS-PS ERR , and the attenuation of the upward propagating wave as measured by the electronics can be expressed as AT UM =20.log 10 (B 1U /B 2U ) =20 . log 10 (A 1U /A 2U ) +20 . log 10 (A 1 /A 2 )=AT-AT ERR The errors introduced by the differences in the receivers and receiver electronics have opposite signs for the upward and downward propagating waves and cancel in the borehole compensated quantities PS BHC and AT BHC : PS BHC =(PS UM +PS DM )/2=PS AT BHC =(AT UM +AT DM )/2=AT Although the above describes the elimination of electronic error, part of the error is due to "path" difference caused by borehole rugosity. Borehole compensation is also known to reduce the effects of borehole rugosity since rugosity causes a difference in the path lengths through the formation to the respective receivers. Without borehole compensation, this difference would cause transient phase shift and attenuation errors. In accordance with an embodiment of the invention, the transmitter T 3 is utilized to obtain two further depths of investigation, and the measurements thereof are borehole compensated, without the need for providing an "equal and opposite" transmitter; i.e., a corresponding transmitter at the same spacing from the receivers as T 3 , but on the other side of the receivers. Error signals are derived from the measurement obtained using T 1 and T 2 , and are used to correct (or compensate) measurements made with transmitter T 3 . The error terms for phase shift and attenuation that were obtained when T 1 and T 2 were active were PS ERR = (PS UM -PS DM ) /2 AT ERR = (AT UM -AT DM ) /2 With transmitter T 3 activated, and transmitters T 1 and T 2 off, the downward propagating wave has a phase P' 1D and an amplitude A' 1D at the location of receiver R 1 and a phase P' 2D and an amplitude A' 2D at the location of receiver R 2 . The true phase shift for the downward propagating wave is PS'=P' 2D -P' 1D , and the true attenuation is AT'=20 . log 10 (A' 2D /A' 1D ). Because transmitter T 3 is closer to the receivers than T 1 and T 2 , PS' and AT' will generally differ from PS and AT. However, the phase shift and attenuation errors introduced by differences between the receivers and the receiver electronics will be unchanged. The voltage derived from R 1 and input to switch 250 has a net phase X' 1D =P' 1D +P 1 and an amplitude B' 1D =A 1 *A' 1D . Similarly, the voltage derived from R 2 and input to switch 250 has a net phase X' 2D =P' 2D +P 2 and an amplitude B' 2D =A 2 *A' 2D . The downward phase shift as measured by the electronics is: PS' DM =X' 2D -X' 1D = (P' 2D -P' 1D ) + (P 2 -P 1 ) =PS'+PS ERR and the downward attenuation as measured by the electronics is: AT' DM =20 . log 10 (B' 2D /B' 1D )=20 . log(A' 2D /A' 1D )+20 . log 10 (A 2 /A 1 ) =AT'+AT ERR These quantities can be corrected by subtracting PS ERR and AT ERR as measured by T 1 and T 2 . That is, the corrected phase shift for T 3 will be PS' C =PS' DM -PS ERR =PS' and the corrected attenuation will be AT' C =AT' DM -AT ERR =AT' Thus, when these corrected measurements are generated, it is not necessary to add a fourth transmitter to obtain the benefits of borehole compensation, including correction for electronics errors and correction for borehole rugosity. This technique reduces the complexity and cost of adding additional depths of investigation to the logging tool, while retaining the advantages of borehole compensation. Significant further advantage is obtained for the situation when T 3 is placed above T 1 , as shown in Figure 3. [Processing will be the same as just described.] In this case, the tool length must only be increased by approximately the distance between T 1 and T 3 . In a conventionally borehole compensated long array, the tool length would have to be increased by twice such distance. It will be understood that further transmitters can be added, consistent with the principles set forth. Referring to Figure 4, there is shown a flow diagram of a routine for controlling the processor 270 in accordance with an embodiment of the invention for obtaining borehole compensated resistivity from attenuation and borehole compensated resistivity from phase shift for the spacing of transmitters (T 1 , T 2 ), and also for the spacing of transmitter T 3 . The block 410 represents initializing of a transmitter index to the first transmitter (T 1 ). The transmitter is activated, and the relative phase shift and attenuation are measured and stored, as represented by the block 412. Inquiry is then made (decision block 415) as to whether the last transmitter index has been reached. If not, the transmitter index is incremented (block 417), block 412 is re-entered, and the loop 419 is continued until measurements have been made and stored for each of the three transmitters (T 1 , T 2 , and T 3 , for the present embodiment). The block 425 is then entered, and, for the T 1 , T 2 spacing, borehole compensated phase shift [PS BHC ] and attenuation [AT BHC ] are computed using equation (5) and equation (6). These measured values can then be used to determine and store (block 427) the resistivity from attenuation [R AT ] and the resistivity from phase shift [R  ] at the spacing of T 1 , T 2 . This determination can be made, for example, in accordance with the teachings of U.S. Patent No. 4,899,112, by using the table look-up technique described in detail therein. [As noted in U.S. Patent No. 4,899,112, the look-up table can take into account dielectric constant, which can be treated as a function of conductivity.] The block 430 is then entered, and the phase shift error [PS ERR ] and the attenuation error [AT ERR ] are determined from equation (7) and equation (8), respectively. Then, in accordance with a feature of the invention, and as represented by the block 435, the corrected (that is, pseudo borehole compensated) phase shift [PS C '] and attenuation [AT C '], for the T 3 spacing are computed from equation (11) and equation (12). The same type of table look-up technique can then be used (block 440) to determine and store the resistivity from attenuation [R' AT ] and the resistivity from phase shift [R'  ] at the T 3 spacing. Thus, using the configuration of Figure 2 (or Figure 3), resistivity can be determined, with borehole compensation, at four effective depths of investigation [namely, R AT for the T 1 , T 2 spacing, R  for the T 1 , T 2 spacing, R' AT for the T 3 spacing, and R'  for the T 3 spacing]. The block 410 can then be re-entered for processing at further depth levels. Referring to Figure 5, there is shown a further embodiment of the invention, which may be employed, for example, as the logging device 200 of the measuring while drilling equipment illustrated in Figure 1. The embodiment of Figure 5 has a receiver pair comprising receiving antennas R 1 , R 2 , and five transmitting antennas, T 1 , T 2 , T 3 , T 4 , and T 5 . The transmitting antennas T 1 , T 3 , and T 5 are above the receiver antennas, and the transmitting antennas T 2 and T 4 are below the receiving antennas. The spacings of the transmitting antennas, with respect to the midpoint of the receiving antennas, (0), are all different. As one example, the spacing between receiving antennas R 1 and R 2 may be 6 inches, and the spacing from the receiver midpoint to each transmitter T 1 through T 5 may be 10 inches, 16 inches, 22 inches, 28 inches, and 34 inches, respectively. Each transmitter can be individually energized (for example, on a time multiplexed or a frequency multiplexed basis), and a relative phase shift measurement and an attenuation measurement can be determined from the signals received at the receiving antennas, for electromagnetic energy transmitted from each of the transmitting antennas T 1 through T 5 . The electronics for the Figure 5 embodiment can be similar to that of the Figure 3 embodiment. In this case,electronic switch 775 receives the signal from amplifier 212, and couples the signal to a selected one of the transmitting antennas T 1 , T 2 , T 3 , T 4 or T 5 , via impedance matching circuits 701, 702, 703, 704 or 705, under control of a signal from processor 270. The processor can energize the transmitting antennas in any desired order, for example, sequentially. The routine of Figure 8 (to be described below) can be used to program the processor 270 to control the energizing of the transmitting antennas and the storing and processing of relative phase shift and attenuation measurements from the receiving antennas. Before further describing the processing of the received signals of the Figure 5 embodiment, consider the general case of Figure 6, which illustrates a Z-axis (which can be a tool axis), with Z=0 being at the midpoint of receiving antenna locations represented by dots at Z R1 and - Z R2 , where Z R1 =Z R2 . Transmitting antennas are located along the Z axis, with the i-th transmitter, Z i , at: Z i =(-1) i-1 . [Z min +(i-1). Z inc ], i=1,2,3,....,n where Z min is the spacing of the closest transmitter, Z inc is the spacing increment (preferably at least one inch), and n is an integer that is at least 3. Accordingly, as seen in Figure 6, the transmitters will be located at Z 1 =Z min , Z 2 =-(Z min +Z inc ), Z 3 = (Z min +2Z inc ), Z 4 =-(Z min +3Z inc ), Z 5 =(Z min +4Z inc ), etc. A three transmitter device with transmitters at Z 1 =Z min , Z 2 =-(Z min +Z inc ) and Z 3 =(Z min +2Z inc ) is illustrated in Figure 7. In an illustrative example, the receivers are spaced 6 inches apart (Z R1 =-Z R2 = 3 inches), Z min is 10 inches, and Z inc is 6 inches. This means that the transmitter positions with respect to the receiver midpoint, for this example, are at Z 1 =10 inches, Z 2 =-16 inches, and Z 3 =22 inches. The electronics coupled with the transmitting and receiving antennas can be of the same type shown in Figure 2, except that the determinations of borehole compensated relative phase shift and attenuation are determined by processor 270 from the measured values in a manner to be described. For example, six effective depths of investigation (from three physical TR spacings), each with borehole compensation, can be obtained using the configuration of Figure 7 and applying appropriate weightings of the relative phase shift and/or attenuation measurements associated with each transmitter. If the relative phase shift or the attenuation measured at the receiving antennas when the transmitter T i is active is designated X(i), the following equivalent spacings, which illustrate the different types of borehole compensation, can be obtained: 1) Equivalent TR spacing of 10 inches: X (e.g. 10 inch) = .75 . X(1) + .5 . X(2) - .25 . X(3) The weighted mixing of measurements (relative phase shift or attenuation) for this type of case is referred to herein as "mix-pseudo", and has some similarity to the embodiment of Figure 2 or 3 (where T 1 and T 2 are symmetrically spaced, whereas here they are not), in that an imbalance function type of correction is used in obtaining borehole compensation. This type of mixing produces an asymmetric response to a thin bed, in varying degrees, depending on the spacing between transmitters and receivers. The response is more asymmetric for shorter spacings and less asymmetric for longer spacings. 2) Equivalent TR spacing of 13 inches: X (eg. 13 inch) = .5 . X(1) + .5 . X(2) The weighted (even weightings in this instance) mixing of measurements for this type of case is referred to herein as "mix-asymmetric". The depth of investigation substantially corresponds to the average of the two components spacings TR (10 inches and 16 inches, for this example), and borehole compensation is again achieved. These is no physical sensor at the equivalent TR spacing (13 inches, in this example), and the result is an interpolated one. This type of mixing also produces an asymmetric response to a thin bed, although the degree of asymmetry is less than that of the previously described mix-pseudo type. 3) Equivalent TR spacing of 16 inches: X (eg. 16 inch) = .25 . X(1) + .5 . X(2) + .25 . X(3) The weighted mixing of measurements for this type of case is referred to herein as "mix-symmetric". This type of mixing produces a symmetric response, with borehole compensation, at the equivalent TR spacing that has a physical sensor (at 16 inches, in this example). Table 1 shows, for the general case of up to n transmitters, the equivalent TR spacing, the mixing type, and the preferred weighting coefficients. A continuum of effective TR spacings may be obtained by adjusting the weighting factors of each measurement subject to the constraints that the sum of the weighting factors of the positive spacings should equal the sum of the weighting factors of the negative spacings, and the sum of all the factors should equal unity. However, the number of independent spacings is only equal to the number of transmitters at different spacings. Equivalent TR spacing type weighting coefficients Zmin mix-pseudo 0.75 . X(1)+0.5 . X(2)-0.25 . X(3) Zmin+0.5 . Zinc mix-asym. 0.5 . X(1)+0.5 . X(2) Zmin+Zinc mix-sym. 0.25 . X(1)+0.5 . X(2)+0.25 . X(3) Zmin+m.Zinc mix-sym. 0.25 . X(m)+0.5 . X(m+1)+0.25 . X(M+2) Zn-Zinc mix-sym. 0.25 . X(n-2)+0.5 . X(n-1)+0.25 . X(n) Zn-0.5 . Zinc mix-asym. 0.5 . X(n-1)+0.5 . X(n) Zn mix-pseudo 0.75 . X(n)+0.5 . X(n-1)-0.25 . X(n-2) Other combinations and weighting coefficients can be used, including, for example, a combination of measurements obtained from using four or more of the transmitters. For a particular combination, the measurements from transmitters not involved in the computation have weightings of zero. Accordingly, the general formula for a particular spacing will be X = w 1 . X(1) + w 2 . X(2) + w 3 . X(3) + ... w n . X(n) where n is the number of transmitters, and w 1 , w 2 , w 3 ... w n are weighting coefficients that can have positive, negative, or zero values, preferably subject to the above-indicated constraints. Reference can then again be made to the embodiment of Figure 5, with five transmitters at the indicated spacings: 10 inches, 16 inches, 22 inches, 28 inches and 34 inches [which, as above, results from the transmitter placement equation (13), with Z min = 10 inches and Z inc = 6 inches]. Table 2 shows equivalent TR spacings that can be obtained with the Figure 5 embodiment, the mixing types, the preferred weighting coefficients, and whether or not a physical sensor is present at each equivalent TR spacing. In this case only a maximum of five spacings are independent. The others are "dependent" spacings which can be linearly derived from combinations of measurements from the independent spacings. It will be understood that, if desired, weightings can be varied while operating downhole, either automatically or under operator control when communications capability is present. Equivalent TR spacing type weighting coefficients physical sensor 10" mix-pseudo 0.75 . X(1)+0.5 . X(2)-0.25 . X(3) Yes 13" mix-asym. 0.5 . X(1)+0.5 . X(2) No 16" mix-sym. 0.25 . X(1)+0.5 . X(2)+0.25 . X(3) Yes 22" mix-sym. 0.25 . X(2)+0.5 . X(3)+0.25 . X(4) Yes 28" mix-sym. 0.25 . X(3)+0.5 . X(4)+0.25 . X(5) Yes 31" mix-asym. 0.5 . X(4)+0.5 . X(5) No 34" mix-pseudo 0.75 . X(5)+0.5 . X(4)-0.25 . X(3) Yes It is seen that in the embodiment of Figure 5, fourteen different depths of investigation (one attenuation-based and one phase shift-based, for each of seven equivalent TR spacings) can be obtained, with borehole compensation, using the indicated preferred weightings. Referring to Figure 8, there is shown a flow diagram of a routine for controlling a processor (e.g. 270 of Figure 2) of the Figure 5 embodiment, with five transmitters, in accordance with an embodiment of the invention. [The same type of routine can be used for any number of transmitters.] The blocks 410, 412, 415, and 417, and the loop 419 operate as described in Figure 4, to sequentially cause energizing of each transmitter (T 1 through T 5 , in this case), and to measure and store the relative phase shift and attenuation of the received signals, associated with each transmitting antenna. The block 851 is then entered, this block representing the initializing of an equivalent spacing index to the first equivalent spacing of the Table 2; namely, the equivalent spacing of 10 inches. The mixed phase shift value is then computed (block 855) using the Table 2 weighting values. Also, the mixed attenuation value is computed, again, using the Table 2 weighting values (block 857). [The blocks 855A and 857A represent loading of modified combinations and/or weighting values, for situations where it is desirable to change and/or add combinations and/or weighting values.] The block 860 is then entered, this block representing, again, use of the table look-up to determine and store the resistivity from attenuation [R AT ] and the resistivity from phase shift [R  ] for the present spacing. A determination is then made (decision block 865) as to whether the last equivalent spacing index has been reached. [In this case, that would be the equivalent spacing of 34 inches.] If not, the equivalent spacing index is incremented (block 870), and the block 855 is re-entered. The loop 872 is then continued until the attenuation-based resistivity [R AT ] and the phase shift-based resistivity [R  ] have been computed for each equivalent spacing. The block 410 can then be re-entered for processing of subsequent depth levels. In all, for the example of this embodiment, borehole compensated resistivity can be obtained at fourteen depths of investigation; namely, two depths of investigation [R AT and R  ] for each of seven equivalent spacings. The invention has been described with regard to particular preferred embodiments, but variations within the spirit and scope of the invention will occur to those skilled in the art. For example, it will be understood that properties of earth formations, such as vertical and horizontal resistivity and dielectric constant, can be determined. Also, further receiving antennas can be used. Finally, while the invention is well suited for measurement-while-drilling application, it can also be applied in wireline logging when the drill string is removed from the borehole.

Claims: Apparatus for determining a property of earth formations surrounding a borehole (11), comprising: an elongated logging device (200) movable through the borehole; a first sequence of at least two spaced apart transmitting antennas (T 1 , T 2 ) on said device (200), and a second sequence of at least one transmitting antenna (T 3 ), on said device (200); a pair of longitudinally spaced receiving antennas (R 1 , R 2 ) on said device, said pair of receiving antennas (R 1 , R 2 ) being located between one of the transmitting antennas (T 2 ) of said first sequence and one of the transmitting antennas (T 3 ) of said second sequence ; means for energising (211) each of said transmitting antennas (T 1 , T 2 , T 3 ) with electromagnetic energy; means for receiving (58), at said pair of receiving antennas (R 1 , R 2 ), electromagnetic energy from said transmitting antennas (T 1 , T 2 , T 3 ); means for determining (58, 85), from electromagnetic energy received at said pair of receiving antennas (R 1 , R 2 ), phase shift and/or attenuation measurements of electromagnetic energy transmitted from each of the respective transmitting antennas (T 1 , T 2 , T 3 ); and means for determining (58, 85) a property of said earth formations from said phase shift and/or attenuation measurements, characterised in that the spacing of the transmitting antennas from said second sequence (T 3 ) from the midpoint (0) between said pair of receiving antennas (R 1 , R 2 ) is different from the spacing of every other transmitting antenna from said first sequence (T 1 , T 2 ) from said midpoint (0) between said pair of receiving antennas (R 1 , R 2 ). Apparatus as defined by claim 1, wherein the spacings of the transmitting antennas (T 1 , T 2 , T 3 ) from the midpoint (0) between the pair of receiving antennas (R 1 , R 2 ) are all different from each other. Apparatus as defined by claim 1 or 2, wherein said spacings of said transmitting antennas (T 1 , T 2 , T 3 ) each differ by at least 1 inch. Apparatus as defined by any of claims 1 to 3, wherein said property is resistivity. Apparatus as defined by claim 4, wherein said means for determining resistivity includes means (58, 85, 200, 270) for combining at least one weighted measurement of phase shift and/or attenuation of electromagnetic energy from at least one of said first transmitting antenna sequence (T 1 , T 2 ), and at least one weighted measurement of phase shift and/or attenuation of electromagnetic energy from at least one of said second transmitting antenna sequence (T 3 ). Apparatus as defined by claim 5, wherein said at least one weighted measurement of phase shift and/or attenuation of electromagnetic energy from at least one of said first transmitting antenna sequence (T 1 , T 3 , T 5 ) comprises weighted measurements of phase shift and/or attenuation of electromagnetic energy from two of said second transmitting antenna sequence (T 2 , T 4 ). Apparatus as defined by any of claims 4 to 6, wherein said means for determining resistivity (58, 85) comprises means for determining borehole compensated resistivity at several depths of investigation. Apparatus as defined by any of claims 1 to 7, wherein said logging device is mounted in a drill string (12) for logging while drilling operation. Apparatus as defined by any of claims 1 to 8, wherein said transmitting antennas (T 1 , T 2 , T 3 , T 4 , T 5 ) are energised in sequence. Apparatus as defined by any of claims 1 to 9, wherein said transmitting (T 1 , T 2 , T 3 , T 4 , T 5 ) and receiving antennas (R 1 , R 2 ) are coils wound on a generally cylindrical drill collar. A method for determining a property of earth formations surrounding a borehole (11), comprising the steps of: providing an elongated logging device (200) movable through the borehole, said device having thereon a first sequence of at least two spaced apart transmitting antennas (T 1 , T 2 ), and a second sequence of at least one spaced apart transmitting antenna (T 3 ), and a pair of longitudinally spaced receiving antennas (R 1 , R 2 ), said pair of receiving antennas being located between one transmitting antenna (T 2 ) of said first sequence (T 1 , T 2 ) and one transmitting antenna (T 3 ) of said second sequence (T 3 ), the spacing of the transmitting antennas (T 3 ) of the second sequence from the midpoint (0) between said pair of receiving antennas being different to the spacing of every other transmitting antenna (T 1 , T 2 ) from the midpoint between said pair of receiving antennas; energising each of said transmitting antennas with electromagnetic energy; receiving, at said pair of receiving antennas, electromagnetic energy from said transmitting antennas; determining, from electromagnetic energy received at said pair of receiving antennas, phase shift and/or attenuation measurements of electromagnetic energy transmitted from each of the respective transmitting antennas (T 1 , T 2 , T 3 ); and determining a property of said earth formations from said phase shift and/or attenuation measurements. The method as defined by claim 11, wherein the spacings of the transmitting antennas (T 1 , T 2 , T 3 ) from the midpoint (0) between the pair of receiving antennas (R 1 , R 2 ) are all different from each other. The method as defined by claim 11 or 12, wherein said property is resistivity, and wherein said step of determining resistivity includes combining at least one weighted measurement of phase shift and/or attenuation of electromagnetic energy from at least one of said first transmitting antenna sequence (T 1 , T 2 ), and at least one weighted measurement of phase shift and/or attenuation of electromagnetic energy from at least one of said second transmitting antenna sequence (T 3 ). The method as defined by claim 13, wherein said property is resistivity, and wherein said at least one weighted measurement of phase shift and/or attenuation of electromagnetic energy from at least one of said first transmitting antenna sequence (T 1 , T 3 , T 5 ) comprises weighted measurements of phase shift and/or attenuation of electromagnetic energy from two of said second transmitting antenna sequence (T 2 , T 4 ). The method as defined by claim 13 or 14, wherein said step of determining resistivity comprises determining borehole (11) compensated resistivity at several depths of investigation. The method as defined by any of claims 11 to 15, wherein said method is performed as a logging while drilling operation. The method as defined by claim 16, wherein said property is resistivity, and wherein said step of determining resistivity includes combining at least one weighted measurement of phase shift and/or attenuation of electromagnetic energy from at least one of said first transmitting antenna sequence (T 1 , T 2 ), and at least one weighted measurement of phase shift and/or attenuation of electromagnetic energy from at least one of said second transmitting antenna sequence (T 3 ), the method further comprising the step of modifying the weightings of said measurements as part of said logging while drilling operation.

Inventors: CLARK BRIAN, GRIJALVA VICTOR E, WU PETER T, CLARK, BRIAN, GRIJALVA, VICTOR E., WU, PETER T.
Applicants: ANADRILL INTERNAT SA, SCHLUMBERGER SERVICES PETROL, ANADRILL INTERNATIONAL SA, SERVICES PETROLIERS SCHLUMBERGER
Application year: 1995
Main classification: G01V3
IPC codes: G01V3