Indirect shearwave determination

• US 4,575,830 A
• Filed: 10/15/1982
• Issued: 03/11/1986
• Est. Priority Date: 10/15/1982
• Status: Expired due to Term

• Current Assignee: Schlumberger Technology Corporation (SLB)
• Original Assignee: Schlumberger Technology Corporation (SLB)
• Inventors: Murray, Josephine L., Ingram, John D.
• Primary Examiner(s): NOT, DEFINED
• Assistant Examiner(s): Kaiser, K. R.
• Application Number:
  US06/434,658
• Time in Patent Office: 1,243 Days
• Field of Search: 367/25, 367/28, 367/31, 367/32, 367/75, 364/422
• US Class Current: 367/26
• CPC Class Codes: G01H 5/00 Measuring propagation velocity of ultrasonic, sonic or infrasonic waves , e.g. of pressure waves
• CPC Class Codes: G01V 1/50 Analysing data

Abstract

A method and apparatus for estimating shear wave velocities where such velocities are difficult or impossible to measure because the shear energy transmission is attenuated or its velocity is equal to or less than direct mud velocity or because the shear wave may not propagate at all under certain conditions. The method and apparatus relate to an indirect technique for determining shear modulus by measuring the phase velocity of the Stoneley wave excitation and using it to estimate the shear modulus and the shear velocity. Using at least two receivers, a window is placed over the Stoneley arrival of the wave detected at each receiver and the Fourier transform of each wave is taken. The difference in phase of the Stoneley wave arrivals is determined. The relative shear modulus involving the Lame constants is estimated by a procedure which includes matching two points on Stoneley wave dispersion curves. From the relative shear modulus, the ratio of formation density to the density of mud, and the compressional velocity of the borehole mud is estimated the shear velocity.

Claims

1. A method of estimating parameters of shear waves included in acoustic waves arriving sequentially at a plurality of points spaced vertically along the length of a borehole comprising the steps of:

(a) generating acoustic energy in the borehole,

(b) receiving the acoustic energy at each of the points after refraction, reflection and direct transmission through and along a formation adjacent the borehole,

(c) establishing a window of known length with a known moveout between adjacent receiving points,

(d) positioning said window along the composite wave relative to the estimated arrival of a selected guided acoustic wave,

(e) multiplying the energy received at each point by the window,

(f) generating the Fourier transform of that portion of received guided wave energy multiplied by the window, to produce a plurality of complex signals in the frequency domain,

(g) determining the cross-spectral magnitude and phase for pairs of signals from adjacent receiving points,

(h) computing the phase velocity α

from the phase at each frequency,

(i) scanning the cross-spectral magnitudes for a peak in a selected frequency range and identifying the frequency at which the peak occurs,

(j) selecting from the phase velocity α

versus frequency a value of guided wave phase velocity α

_s at the identified frequency,

(k) estimating the ratio of Lame constants μ

₂ and λ

₁ in accordance with the expression ##EQU5## where A and γ

are values found by matching two points on known dispersion curves for said guided wave where α

, represents mud velocity,

(l) obtaining an estimate of shear wave velocity β

₂ in accordance with the expression ##EQU6## where ρ

₂ is formation density and ρ

₁ is mud density, and (m) plotting the estimate of β

₂ as a function of depth.

2. The method of claim 1 wherein the selected guided wave is a Stoneley wave.

3. The method of claim 1 wherein the window is a Blackman window.

4. The method of claim 3 wherein the Blackman window has a duration of 1200 microseconds.

5. The method of claim 1 including the step of time shifting the windowed energy prior to Fourier transformation to compensate for large phase differences.

6. The method of claim 2 including the step of time shifting the windowed energy prior to Fourier transformation to compensate for large phase differences.

7. The method of claim 3 including the step of time shifting the windowed energy prior to Fourier transformation to compensate for large phase differences.

8. The method of claim 5 wherein the amount of time shift introduced to the windowed energy is removed following Fourier transformation.

9. The method of claim 6 wherein the amount of time shift introduced to the windowed energy is removed following Fourier transformation.

10. The method of claim 7 wherein the amount of time shift introduced to the windowed energy is removed following Fourier transformation.

11. The method of claim 1 wherein the mud velocity α

, is assumed in a range of 5555 ft/sec to 5000 ft sec.

12. The method of claim 1 wherein the mud velocity α

, is approximated from a measured value of mud density.

13. The method of claim 1 wherein the ratio ρ

₂ /ρ

₁ is assumed a constant.

14. The method of claim 13 wherein the assumed value of the ratio is 2.5.

15. A system for estimating parameters of shear waves included in acoustic waves arriving sequentially at a plurality of points spaced vertically along the length of a borehole comprising the steps of:

(a) means for generating acoustic energy in the borehole,

(b) means for receiving the acoustic energy at each of the points after refraction, reflection and direct transmission through and along a formation adjacent the borehole,

(c) means for establishing a window of known length with a known moveout between adjacent points and for positioning said window along the composite wave relative to the estimated arrival of Stoneley waves,

(d) means for multiplying the energy received at each point by the window,

(e) means for generating the Fourier transform of that portion of received Stoneley wave energy multiplied by the window, to produce a plurality of complex signals in the frequency domain,

(f) means to determine the cross-spectral magnitude and phase for pairs of signals from adjacent points,

(g) means to compute the phase velocity α

from the phase at each frequency,

(h) means to scan the cross-spectral magnitudes for a peak in a selected frequency range and identify the frequency at which the peak occurs,

(i) means to select from the phase velocity α

versus frequency a value of Stoneley phase velocity α

_s at the identified frequency,

(j) means for estimating the ratio of Lame constant μ

₂ and λ

₁ in accordance with the expression ##EQU7## where α

, represents mud velocity, and A and γ

are value found by matching two points on known Stoneley wave dispersion curves,

(k) means for obtaining an estimate of shear wave velocity β

₂ in accordance with the expression ##EQU8## where ρ

₂ is formation density and ρ

₁ is mud density, and (l) means for plotting the estimate of ρ

₂ as a function of depth.