WO1993007513A1

WO1993007513A1 - Methods and apparatus for discrete-frequency tube-wave logging of boreholes

Info

Publication number: WO1993007513A1
Application number: PCT/US1991/007494
Authority: WO
Inventors: Shu-Kong Chang; Chaur-Jian Hsu
Original assignee: Chang Shu Kong; Hsu Chaur Jian
Priority date: 1991-10-11
Filing date: 1991-10-11
Publication date: 1993-04-15

Abstract

The present invention offers methods and apparatus for low-frequency sonic borehole logging with tube waves using discrete frequency measurements. A source transducer (16) emits sonic energy at a location in a borehole of at least one discrete frequency which induces propagation of tube waves in the borehole (14); sonic energy of the tube waves is received at multiple receiver locations (20, 22, 24, 26, 28, 30) in the borehole, the multiple receiver locations (20, 22, 24, 26, 28, 30) being spaced apart from one another and from the source transducer location (16); and at each said receiver location (20, 22, 24, 26, 28, 30), the complex pressure response (amplitude and phase, or real and imaginary parts) of a receiver to the received sonic energy of the tube waves is detected relative to a phase reference. The detected complex pressure response may be processed to determine phase velocity of the tube waves in the borehole as a function of frequency, to determine attenuation of the tube waves in the borehole as a function of frequency and/or to determine attenuation and amplitude as functions of frequency. In addition, tube-wave reflection coefficients near fractures, bed boundaries, etc., may be calculated. The tube-3ave data are useful for permeability correlation, fracture evaluation, rock rigidity studies and other purposes.

Description

METHODS AND APPARATUS FOR DISCRETE-FREQUENCY TUBE-WAVE LOGGING OF BOREHOLES

BACKGROUND OF THE INVENTION

1. Field of the Invention

This invention relates to logging of low-frequency tube waves in boreholes, and more particularly to logging of tube waves propagated in boreholes by emission of sonic energy at discrete frequencies.

2. Description of Related Art

Emission of broad-band sonic energy in a fluid-filled tube, such as a mud-filled borehole in the earth, induces propagation of sonic waves in a number of modes. These modes include compressional. shear, and tube wave propagation.

Tube waves are a special type of surface waves. Surface waves were first postulated by Lord Rayleigh in 1885. Lord Rayleigh considered the case of a flat plane separating some solid elastic material from a vacuum, and found that the motion of the material near the plane executed a sort of elliptical motion as the surface wave passed. In 1924. R. Stoneley examined the related case in which a second solid is substituted for the vacuum, and discovered a new type of surface wave. In 1948. J.G. Scholte considered the special case when one of the materials was a fluid. The case considered by Scholte is the one of interest with respect to a fluid-filled borehole, although Stoneley's name is more commoπlv used in describine it. If Scholte's plane is bent to form a cylinder, the geometry of the borehole is encountered. At very small wavelengths relative to a given borehole diameter, the cylinder will appear to the Stoneley wave as a flat plane. But with longer wavelengths, the Stoneley wave decays very little across the borehole. For example, at an operating frequency of, for example, 10 KHz. the wavelength in the mud is about 6 inches (comparable with a typical oilfield borehole diameter) and the wavelength in the formation surrounding the borehole may be as long as 2 feet. At the very low-frequency limit, the field strength of the Stoneley wave is practically flat across the borehole, since the borehole diameter is very small compared to a wavelength. Under these conditions the Stoneley wave is known as a tube wave. For purposes of the present invention and the disclosure given below, the terms "tube wave" and "Stoneley wave" will be considered synonymous. Propagation of the tube wave up and down the borehole is believed to be affected by fluid communication between the borehole and the formation. Low-frequency tube waves have found a number of uses in borehole exploration, such as for determining formation permeability, evaluating formation fractures, and understanding rock rigidity and stress. Correlation between tube wave parameters (velocity and attenuation) and formation permeability has been reported in theory and from field data. See, for example, J.E. WHITE, UNDERGROUND SOUND, pp. 139-191, Elsevier. 1983: and D.M. WILLIAMS, J. ZEMANEK, F.A. ANGONA, CL. DENNIS, and R.L. CALDWELL, THE LONG SPACED ACOUSTIC LOGGING TOOL (Transactions of the SPWLA 25th Annual Logging Symposium. Paper T. 1984). The tube wave has also been used as a fracture indicator. See. for example. K. Hsu. A. BRIE, and R. PLUMB. A NEW METHOD FOR FRACTURE IDENTIFICATION USING SONIC ARRAY TOOLS (Paper presented at the 1985 Annual Technical Conference and Exhibition of the Society of Petroleum Engineers of AIME. September 22 - 25, 1985, publication no. SPE 14397). Use of tube wave information to derive the formation shear modulus is also known. See, for example. O.Y. LIU, STONELEY WAVE-DERIVED Δt

SHEAR LOG (Paper presented at the SPWLA Twenty-Fifth Annual Logging Symposium. June 10-13. 1984); and J.L. STEVENS and S.M. DAY, SHEAR VELOCITY LOGGING IN SLOW FORMATIONS USING THE STONELEY WAVE (Borehole Geophysics Abstract No. BHG7 of the Extended Abstracts of the Annual Meeting of the Society of Exploration Geophysicists. September. 1983). Traditionally, sonic logging measurements are made in the time domain: a broad-band sonic energy source excites propagation of sonic waves in the borehole, and waveforms detected at an array of receiver locations spaced from the source in the borehole are recorded as functions of time. See. for example, D.M. WILLIAMS. J. ZEMANEK, F.A. ANGONA, CL. DENNIS, and R.L. CALDWELL, THE LONG SPACED ACOUSTIC LOGGING TOOL (Transactions of the SPWLA 25th Annual Logging Symposium. Paper T. 1984). Tube wave velocity and reflection coefficients at the dominant tube wave frequency can be estimated by waveform stacking techniques, such as semblance or radon transforms. However, quantitative interpretation and applications of the tube wave data require information about the tube waves in the frequency domain. That is, information such as tube wave velocity dispersion, tube wave attenuation and tube wave reflection coefficients are needed as functions of frequency. To obtain the frequency-dependent information, waveforms recorded in the time domain are transformed into the frequency domain before processing.

Time domain recording has several disadvantages if the frequency domain results are desired. The waveforms must be sampled in small time steps over a long period, resulting in a large set of data to transmit, store, and process. The truncation in time can cause interference in the processing. Also, the signal-to- noise ratio can be quite small for the wide-band transient measurements. The noise problem can be significant for low-frequency data (less than 500 Hz) because of the noises generated by the sonde traveling in the borehole. These problems may be avoided in accordance with the present invention by taking measurements in frequency domain.

U.S. Patent No. 3,330.375. issued July 11, 1967 to J.E. White proposes a form of acoustic well logging in which the propagation velocities of compressional, shear, mud and casing waves are determined from the expression velocity = frequency × wavelength by employing transmitter and/or receiver tuning techniques to determine the wavelength for a known frequency. A variety of such techniques are disclosed, involving wavelength tuning by varying frequency and/or phase shift. In all the techniques, frequency and/or phase shift is adjusted until an amplitude peak is observed, the frequencies which produce such amplitude peaks are noted, and wave propagation velocity is calculated from the expression given above. In one such technique, the wavelength of the transmitted signal is tuned to a fixed wavelength of the sending and receiving arrays by varying its frequency. In another, the wavelength of the arrays is tuned to a fixed value of the transmitted wavelength by varying the phase shift between adjacent transducer elements. In another, the frequency of the transmitted signal or the phase shift between adjacent receiving transducers, or both, are varied.

The technique of U.S. Patent No. 3.330.375. although not directed to tube wave logging, would have the disadvantage if used for tube wave logging that insufficient information is gathered to permit derivation of both tube wave velocity and tube wave attenuation as functions of frequency. Furthermore, the disclosed method requires downhole tuning of the frequency and/or the wavelength while acquiring the data. The method is not directly applicable to obtain the frequency dependent characteristics of dispersive waves, such as the tube wave in a borehole.

In another approach. U.S. Patent No. 4.419.748. issued December 6. 1983 to R.W. Siegfried, II., proposes a continuous wave sonic logging method in which a continuous sine wave at a single frequency is emitted and received, and a spatial Fourier transform is performed over the receiver array. The resulting spatial frequency component are then used to indicate the velocities of various sonic paths. The logging method of U.S. Patent No.4,419,748 would have several disadvantages for tube wave logging. For example, the method requires a large number of receivers in order to facilitate the spatial Fourier transforms. Further, the disclosed method records the instantaneous values of the received signal: that is. it is recorded in time domain. Therefore, the measurement is subject to noise interferences. No improvement in the signal to noise ratio is realized by the proposed method. Furthermore, the measurement can only be done one frequency at a time (due to the time domain recording). This would require numerous logging runs for the dispersive waves, for which the wave characteristics are functions of frequency. Since the logging time is a costly factor in wire iine logging services, the method is not practical for logging dispersive waves.

It is an object of the present invention to provide methods and apparatus for borehole logging with tube waves in which the aforementioned disadvantages of time-domain recording are avoided. It is a further object of the present invention to provide methods and apparatus for borehole logging with tube waves in which the complex pressure response (amplitude and phase, or real and imaginary parts) of tube waves propagated in the borehole is detected for use in determining parameters such as tube wave phase velocity and tube wave attenuation as functions of frequency.

SUMMARY OF THE INVENTION

The present invention offers methods and apparatus for low-frequency sonic borehole logging with tube waves using discrete frequency measurements. In a preferred form of the invention a source transducer emits sonic energy at a location in a borehole of at least one discrete frequency which induces propagation of tube waves in the borehole; sonic energy of the tube waves is received at multiple receiver locations in the borehole, the multiple receiver locations being spaced apart from one another and from the source transducer location: and for each said receiver location, the complex pressure response of the received sonic energy of the tube waves is detected relative to a phase reference. The source transducer may be driven continuously at one or a plurality of discrete frequencies, or may be driven intermittently, or may be driven alternately from a plurality of sources which provide one or a plurality of discrete frequencies. In one preferred form, the source transducer is driven from a square-wave signal source to emit sonic energy at a plurality of discrete frequencies. The phase reference is preferably derived from the same signal source which drives the source transducer.

The detected complex pressure response (amplitude and phase, or real and imaginary parts) of the received sonic energy of the tube waves may be processed to determine phase velocity of the tube waves in the borehole as a function of frequency, to determine attenuation of the tube waves in the borehole as a function of frequency and/or to determine attenuation and amplitude as functions of frequency. In addition, tube wave reflection coefficients near fractures, bed boundaries, etc.. can be calculated. The tube wave data are useful for permeability correlation, fracture evaluation, rock rigidity studies and other purposes. A number of different techniques may be used to process the tube wave data. Among these are: three receiver analytical solutions, least square fittings, two receiver phase difference and amplitude ratio. Prony's method, wave impedance methods and wall impedance methods.

BRIEF DESCRIPTION OF THE DRAWING

Figure 1 shows in schematic block form a first preferred embodiment of a discrete-frequency tube-wave logging apparatus in accordance with the invention:

Figure 2 shows in schematic block form a modification of the embodiment of Figure 1 which permits multiple, discrete frequencies to be emitted sequentially in bursts for tube-wave logging in accordance with the invention:

Figure 3 shows in schematic block form a further modification of the embodiment of Figure 1 which permits multiple, discrete frequencies to be emitted simultaneously for tube-wave logging in accordance with the invention:

Figure 4 shows amplitude of a periodic signal as a function of time, the signal having a discrete-frequency spectrum as shown in Figure 5: Figure 5 shows amplitude of the discrete frequencies of the signal of Figure 4 as a function of frequency; and

Figure 6 shows in schematic block form a further modification of the embodiment of Figure 1 which permits multiple, discrete frequencies to be emitted simultaneously for tube-wave logging in accordance with the invention.

DETAILED DESCRIPTION OF THE INVENTION

Figure 1 shows in schematic block form a first preferred embodiment of a logging apparatus in accordance with the present invention, useful for recording tube-wave data at discrete frequencies. Referring to Figure 1, acoustic logging sonde 10 is shown suspended by cable 12 in borehole 14 in the earth.

Borehole 14 is filled with drilling mud or other fluids (not illustrated). Logging sonde 10 includes a sending, or source, transducer 16 and a plurality of receiving transducers 20 - 30 arrayed at locations spaced apart from one another and from sending transducer 16 along sonde 10. While Figure 1 shows six receiving transducers, it will be appreciated that any number may be used.

Narrow-band waveform generator 32 generates sinusoidal signals at one or more discrete frequencies, the signals being passed through power amplifier 34 to drive sending transducer 16. The signals may be continuous, or may be transmitted in bursts of sufficient duration (for example, on the order of 10 to 20 cycles per burst) to approximate a continuous wave for purposes of exciting propagation of tube waves in the borehole at the desired frequency or frequencies. When so driven, sending transducer 16 emits sonic waves in the borehole at the prescribed frequencies. The sonic waves propagate in the borehole and are received at receiving transducers 20 - 30. The signal output from each of receiving transducers 20 - 30 is supplied to a respective lock-in amplifier which detects the complex pressure response (amplitude and phase, or real and imaginary pans) of the received waves at the frequencies of generator 32. The lock-in amplifiers (shown collectively at 36) also receive phase reference signals from generator 32. The phase reference signals are preferably the signals generated by generator 32 for driving sending transducer 16. Data detected by lock-in amplifiers 32 is passed to conventional telemetry circuitry 36 for transmission to an up-hole system for recording and processing the data. Generator 32, power amplifier 34, lock-in amplifiers 36 and telemetry circuitry 38 are preferably built into sonde 10. Telemetry circuitry 38 preferably transmits data output from lock-in amplifiers 36 to the up-hole recording and processing system 40 via cable 12 or any other suitable means. Sending transducer 16 may be any conventional type of monopole pressure wave source comprising, for example, one or more piezoelectric ceramic transducers, magnetostrictive transducers, or electrodynamic vibrators. Sending transducer 16 may be driven in any manner suitable for exciting propagation of tube waves in the borehole at one or more discrete frequencies, only some of which will be described herein. Receiving transducers 20 - 30 are of any type suitable for detecting pressure waves, such as hydrophones.

Figure 2 shows in schematic block form one embodiment of waveform generator 32 of Figure 1, which permits multiple, discrete frequencies to be emitted in sequential bursts for tube-wave logging in accordance with the invention. A plurality of resonators, for example three resonators 42. 44 and 46, generate sinusoidal signals at respective frequencies suitable for exciting tube waves in the borehole, for example at frequencies of 1 KHz, 3 KHz and 5 KHz. Suitable switching means 48 alternately passes multi-cycle bursts of the respective frequencies to power amplifier 34 to drive transducer 16 and to lock- in amplifiers 36 for use as a phase reference. The multi-cycle bursts are of sufficient length (for example, 10 to 20 cycles) to effectively simulate a continuous signal for purposes of exciting tube wave propagation in the borehole at the selected frequencies.

Figure 3 shows in schematic block form another embodiment of waveform generator 32 of Figure 1 which permits multiple, discrete frequencies to be emitted simultaneously for tube-wave logging in accordance with the invention. A plurality of resonators, for example three resonators 50, 52 and 54, generate sinusoidal signals at respective frequencies. The signals are supplied to mixer 56 which produces multiple frequencies suitable for exciting tube waves in the borehole, for example at frequencies of, for example 1 KHz. 3 KHz and 5 KHz. The multiple-discrete-frequency output of mixer 56 is passed through low-pass filter 58 to remove higher-frequency harmonics, and then through power amplifier 34 to drive sending transducer 16. Mixer 56 also supplies a phase-reference signal for each frequency to a lock-in amplifier. In this case, a lock-in amplifier is provided for each receiving transducer for each frequency. With the embodiment of Figure 3, the complex pressure response (phase and amplitude, or real and imaginary parts) of the received tube waves may be detected simultaneously for a plurality of discrete frequencies.

In yet another embodiment, the waveform generator 32 of Figure 1 drives the source transducer or transducers at a desired set of frequencies by using a repetitive signal sequence. If the period of the sequence is τ . then a discrete frequency spectrum at f = n/τ is obtained, where n = 1. 2. 3, etc. The amplitude distribution for the frequencies of interest can be designed by the shape of the waveform of the sequence. In other words, the shape of the waveform can be designed such that a desired discrete frequency spectrum is obtained. This may be achieved by first specifying the amplitude distribution for the discrete frequencies, and then by summing the Fourier series of all the discrete frequency terms with the pre-assigned amplitude distribution.

One example of such a repetitive signal sequence is given if Figure 4. which shows the amplitude of a periodic signal as a function of time. This periodic signal has a discrete frequency spectrum at f_n = 0.5, 1, 2. 3, 4, 5 kilohertz. The Fourier amplitude for each frequency is inversely proportional to f_n.

Thus, the Fourier spectrum of the signal of Figure 4 is as shown in Figure 5. For purposes of the present invention, the Fourier spectrum of the signal preferably has much higher energy in its lower frequency components than in its higher frequency components. This is because the upper frequency boundary for the energy of tube waves propagating in a borehole is about 5 kHz for boreholes of a diameter typically encountered in oilfield well logging. (It will be recalled that this upper frequency boundary is dependent on borehole diameter and the formation surrounding the borehole.) Greater energy at the low-frequency end of the spectrum will thus provide greater signal-to-noise ratio of the tube waves. Selecting minimum and maximum frequencies of f_min and f_max , respectively, then the square wave signal may be chosen such that r = 1 / f_min , if every f_n is an integer multiple of the minimum frequency. f_min . and N≈ f_max τ A low-pass filter can be applied to the periodic signal to eliminate the harmonics above f_max . Other (e.g., square-wave, saw-tooth, specially-designed pulse shape, etc.) periodic signals may also be used for exciting discrete frequency signals.

Figure 6 shows still another embodiment of the waveform generator 32 of Figure 1 which permits multiple, discrete frequencies to be emitted

simultaneously for tube-wave logging in accordance with the invention.

Referring now to Figure 6, the output signal from square-wave generator 60 is supplied to low-pass filter 62 which passes frequencies f≤ N/τ . The filtered signal is then passed through power amplifier 64 to drive sending transducer 16. The emitted sonic energy from sending transducer 16 excites tube wave propagation in the borehole and the tube waves are received at each of the receiving transducers, only one of which is shown (at 20). The output signal from each of the receiving transducers is passed through a respective band pass filter (one of which is shown at 66) having a pass-band frequency characteristic f_min < f < f_max . The filtered signal is then supplied to the input of a phase-sensitive detector circuit (one of which is shown at 68) which detects the phase of each frequency component of the filtered signal relative to a respective phase reference signal, as well as the amplitude of each frequency component. The phase reference signals may be obtained from the circuitry which drives sending transducer 16. Phase-sensitive detector 68 may be of any suitable analog, digital or hybrid construction. The phase and amplitude data for each frequency is passed by suitable means (for example, telemetry to an up-hole site) from phase-sensitive detector 68 to processor 70 which makes use of the data for any of a number of purposes, some of which are described below.

Those of skill in the an will recognize that suitable modification of the detection circuitry will permit detection of the real and imaginary parts of the complex pressure response of the receiving transducer, rather than the amplitude and phase of the complex pressure response.

As noted above, the upper frequency boundary for the energy of tube waves propagating in a borehole is about 5 KHz for boreholes of a diameter typically encountered in oilfield well logging. This upper frequency boundary is dependent on borehole diameter and the formation surrounding the borehole. Thus, the discrete frequencies at which tube waves are to be excited in the borehole preferably includes a number of discrete frequencies up to the upper frequency boundary, for example, 1 kHz, 2 kHz, 3 kHz, 4 kHz and 5 kHz. Spacing of the receiving transducers is not critical, but is preferably such that the length of the receiving transducer array is not substantially less than the wavelength of the tube wave energy in order to accurately determine the tube wave phase velocity.

Measuring tube waves in frequency domain offers several advantages. For example, measurements in frequency domain give higher signal-to-noise ratio than do time domain measurements. This is because the signal-to-noise ratio is inversely proportional to the bandwidth of the signal, and the bandwidth of the signal can be made extremely narrow when working with discrete frequencies. Also, as the present invention records the complex pressure response

(amplitude and phase, or real and imaginary pans) of the detected tube wave energy for discrete frequencies, the quantity of data recorded is less than that of time-domain logging techniques where full waveforms are recorded. This has at least two advantages. First, as the quantity of data to be transferred uphole is substantially reduced, the logging speed can be increased. Second, it may be possible to store all data in the sonde until logging is complete, thereby also increasing logging speed and efficiency. Further, as there is no need to transform the waveforms from time domain to frequency domain, the processing is simplified. In addition, transducer calibration at specific frequencies is more reliable than is transducer calibration for the broad-band time domain sondes. These features are believed to make the frequency domain logging methods and apparatus of the present invention more accurate than traditional time-domain sonic logging for obtaining the frequency- dependent information of tube waves.

The data obtained with the methods and apparatus of the present invention are the amplitude and phase (or real and imaginary parts) of the tube wave signal at each receiver for each frequency. In the following, selected prior art processing techniques are described which may be used to extract information from the measured data. Other processing techniques may also be employed. 1. Two Receiver Technique

As a tube wave travels in the borehole in one direction, the pressure function can be expressed in the following form: P(z, f) = A(f) e^{-α z} e^{i 2z(f/v) z} (1) where P is the pressure; A is the amplitude function which is independent of transmitter-receiver spacings; is the transmitter-receiver spacing; f is the frequency; v is the tube wave phase velocity: and α is the attenuation constant.

If the complex pressure function is measured at two receiver locations, z₁ andz₂ the tube wave velocity and attenuation constant mav be calculated according to the following formulas: v( f)^-1 = [ phase(P(z₁, f)) - phase(P(z₂, f)) ] / 2 πf (z₁ - z₂) (2a) α = In [ abs(P(z₂, f)/P(z₁, f) ) ] / (z₁ - z₂) (2b)

The two receiver method has one constraint: the tube wave must travel in only one direction. If there are significant reflection signals from tool housing structure, bed boundaries, or fractures, then other processing procedures should be used which can include both forward and backward traveling waves.

2. Three Receiver Technique

Tube waves can be reflected, either by the borehole/formation inhomogeneity or by the sonde housing structure. With the tube wave traveling up and down the borehole, the pressure function can be expressed as:

P(z,t) = A e^ikz + B e^-ikz (3a) where k is the complex wave-number of the tube wave. k = 2π (f/v) + i α (3b)

If P_n is defined as the complex pressure response of the n-th receiver located at z_n , and if the receivers are equally spaced at interval of Δz . then Equation

(3a) can be rewritten as

P_n = A e^ikzn + B e^-ikzn (4)

Taking the data.from only three receivers, and defining A' = A e^ikz2 , B' = B e^-ikz2 and G = e^ikΔz, then the following equations result: P₁ = A'G^-1 + B'G (5a)

P₂ = A' + B ' (5b)

P₃ = A'G + B'G^-1 (5c)

Eliminating A' and B', Equations (5a) - (5c) can be reduced to a quadratic equation of G. The solution of the quadratic equation is

Once G is solved, k can be obtained from the expression k = (In G)/(iΔz). The tube wave phase velocity, v , and the tube wave attenuation constant, α , can then be calculated from Equation (3b). Solving for A' and B' from

Equations (5a), (5b) and (5c), the following expressions are obtained:

The tube wave reflection coefficient can be obtained from the ratio of B' and A' with proper exponential factors.

3. Array Techniques

Any of a number of known techniques may be used to process data from an anay of receivers, such as: least square method, Prony's method, maximum likelihood method, and others.

In the least square method, A and B and k are found such that

is minimized. Here, the total number of receivers in the array is N . The least square method is preferable if the tube wave is the only arrival at the receivers (so that Equation (4) is valid). If there are other arrivals, such as compressional. fluid and other waves, then other spectral estimation techniques (Prony's method, maximum likelihood method, etc.) are preferred.

In applying Prony's method, the receiver pressure at a fixed frequency is first assumed to be

A linear predictor is defined as h(z) so as to satisfy the relationship

P(z)ⓩ h(z) = 0 (10) where the notation, ⓩ , means convolution in z. If H(k) is the Fourier transform of h(z) in the wave-number, domain. Equations (9) and (10) can be satisfied only when H(k) is zero at k = k_m.

In Prony's method, the spatial variables P(z) and h(z) are first truncated to finite intervals of z. The linear equations of Equation (10) are then used to solve for h(z) (assuming h(0) = 1 ). The function h(z) is Fourier-transformed to H(k). Solving the k's for the equation H(k) = 0 , the estimates of k_m's are obtained. The tube wave phase velocity and attenuation are then obtained from the k_m's according to Equation (3b). The amplitude functions A _m's can also be estimated by least square fitting. (This is different from the least square method described above. Here, the calculated values of k_m are used in the fitting.) Further information on Prony's method may be found, for example, at S.L. MARPLE, JR., DIGITAL SPECTRAL ANALYSIS, Prentice Hall, 1987 (incorporated herein by this reference for all purposes), at pages 303-349.

In the maximum likelihood method, the function L(k) is first calculated as

where and

The superscript denotes conjugate transpose of the matrix or the vector. The values of k_m are obtained by searching for the locations of the local minimum of L(k) . Further information on the maximum likelihood method may be found in Capon. Greenfield and Kolker, Multidimensional Maximum Likelihood Processing of a Large Aperture Seismic Array, 55 PROCEEDINGS OF THE IEEE. 1967 (incorporated herein by this reference for all purposes), at page 192-211.

Tube wave measurements made in the frequency .domain in accordance with the present invention have several applications. For example, the velocity and attenuation values may be correlated to the fluid flow permeability of the rock formations: the reflection coefficient may be used to evaluate fracture parameters such as permeability, height, depth, etc.: and/or the wave impedance and wall impedance may be used to obtain fracture parameters and

permeability. In addition, tube wave data as functions of frequency can be related to formation rock rigidity, as well as borehole alteration (caused by stress release or other factors ): and accurate tube wave measurements are useful for improving VSP processing and interpretation. The preferred embodiments described above are not intended to be limiting, but are instead intended as merely illustrative of the present invention. Those of skill in the an will recognize that many modifications may be made in the disclosed embodiments without departing from the spirit and scope of the present invention as defined by the following claims.

Claims

1. A method of borehole logging, comprising the steps of:

(a) emitting sonic energy at a first location in a borehole of at least one discrete frequency which induces propagation of tube waves in said borehole:

(b) receiving sonic energy of said tube waves at multiple receiver locations in said borehole, said multiple receiver locations being spaced apart from one another and from said first location; and

(c) for each said receiver location, detecting a complex pressure response of the received sonic energy of said tube waves.

2. The method of claim 1, further comprising the step of repeating steps (a) - (c) for each discrete frequency.

3. The method of claim 1, further comprising the step of deriving a phase reference from said emitted sonic energy, and wherein said step of detecting a complex pressure response of the received sonic energy of said tube waves comprises the step of detecting said complex pressure response of the received sonic energy of said tube waves in relation to said phase reference.

4. The method of claim 1. further comprising the step of processing said detected complex pressure response to determine attenuation of said tube waves in said borehole.

5. The method of claim 2. further comprising the step of processing said detected complex pressure response to determine attenuation of said tube waves in said borehole as a function of frequency.

6. The method of claim 1, further comprising the step of calculating tube wave phase velocity v and tube wave attenuation constant α from said detected complex pressure response at two of said receiver locations and from the expressions : v(f)^-1 = [ phase(P(z₁, f)) - phase(P(z₂, f)) ] / 2πf (z₁ - z₂) α = In [ abs(P(z₂, f)/P(z₁, f) ) ] / (z₁ - z₂) where the tube wave pressure function is expressed as:

P(z, f) = A(f) e^{-α z} e^{i 2π(f/v) z} and where P is pressure: A is amplitude as a function of frequency and is independent of spacing between said emitting location and said receiver locations: z is the spacing between said emitting location and a receiver location; f is the discrete frequency at which the tube wave energy is propagated: and z₁ and z₂ are respective said receiver Iocations.

7. The method of claim 1. further comprising the step of calculating tube wave phase velocity v and tube wave attenuation constant α from said detected complex pressure response at three of said receiver locations equally spaced at intervals Δz and from the expression : k = 2ττ (flv) + i α = (In G)/(iΔz)

where

and where k is the complex wave-number of the received tube wave, f is the frequency of the received tube wave, and P_n is the complex pressure response of the receiver located at an n-th one of said receiving locations. 0. The method of claim 7, further comprising the step of calculating a tube wave reflection coefficient from the ratio of B' and A' in which

9. The method of claim 1, further comprising the step of determining values for a variable A and a variable B and a variable k such that the expression:

is minimized and where N is a plural number of said receiver locations. k is the complex wave-number of the received tube wave, z_n is an n-th one of said receiving locations, and P_n is the complex pressure response of the receiver located at z_n and is expressed as:

P_n - A e^ikzn + B e^-ikzn

10. The method of claim 9, further comprising the step of calculating the tube wave reflection coefficient from the ratio of said determined values for A and B .

11. The method of claim 9, further comprising the step of calculating the tube wave phase velocity v from the expression: k = 2π (f/v) + i α

12. The method of claim 9. further comprising the step of calculating the tube wave attenuation constant a from the expression: k = 2π (f/v) + i α

13. The method of claim 1, further comprising the steps of

(a) defining the detected complex pressure response P(z) at a fixed frequency to be

(b) defining a linear predictor h(z) so as to satisfy the relationship

P(z)ⓩ h(z) = 0 where the notation. ⓩ , means convolution in z:

(c) defining H(k) as the Fourier transform of h(z) in the wave- number domain, where said definitions of receiver pressure P(z) and said linear predictor h(z) can be satisfied only when H(k) is zero at k = k_m ;

(d) truncating spatial variables P(z) and h(z) to finite intervals of z;

(e) solving the linear equations of the relationship

P(z)ⓩ h(z) = 0 for h(z) , assuming h(0) = 1 ;

(f) transforming the function h(z) to H(k);

(g) calculating values of k for the equation H(k) = 0 to obtain estimated vaiues for k_m .

14. The method of claim 13. further comprising the step of calculating the tube wave phase velocity and tube wave attenuation from said values for k_m from the relationship k = 20 (f/v) + i α

15. The method of claim 13. further comprising the step of estimating said amplitude functions A_m by least square fitting.

16. The method of claim 1, further comprising the steps of (a) calculating the function L(k) as

where

and

and where P_n is the complex pressure response of the receiver located at z_n , N is the number of receivers in the array, and the superscript denotes conjugate transpose of the matrix or the vector: and

(b) determining values of k_m by searching for the locations of the local minimum of L(k) .

17. Apparatus for borehole logging, comprising:

(a) means for emitting sonic energy at a first location in a borehole of at least one discrete frequency which induces propagation of tube waves in the borehole;

(b) means for receiving sonic energy of said tube waves at multiple receiver locations in said borehole, said multiple receiver locations being spaced apart from one another and from said first location: and

(c) means for detecting, for each said receiver location, a complex pressure response of the received sonic energy of said tube waves.

18. The apparatus of claim 17, wherein said emitting means comprises a source transducer and means for driving said source transducer to continuously emit sonic energy at said at least one discrete frequency.

19. The apparatus of claim 17. wherein said emitting means comprises a source transducer and means for driving said source transducer to intermittently emit sonic energy at said at least one discrete frequency.

20. The apparatus of claim 17. wherein said emitting means comprises a source transducer and means for driving said source transducer to emit sonic energy at a plurality of discrete frequencies.

21. The apparatus of claim 17, wherein said emitting means comprises a source transducer and means for driving said source transducer to emit sonic energy at in sequential bursts at a variety of discrete frequencies, and wherein, for each said burst, said detecting means detects for each said receiver location amplitude of the received sonic energy of said tube waves and phase of the received sonic energy of said tube waves relative to said predetermined phase reference.

22. The apparatus of claim 17, further comprising means for supplying a phase reference to said detecting means, and wherein said detecting means comprises means for detecting a complex pressure response of the received sonic energy of said tube waves in relation to said phase reference. 23. The apparatus of claim 17. further comprising means for processing said detected complex pressure response to determine attenuation of said tube waves in said borehole.

24. The apparatus of claim 17, further comprising means for processing said detected complex pressure response to determine attenuation of said tube waves in said borehole as a function of frequency.