Patent US5057835A successfully refreshed from Google Patents.
← Back to Patent Explorer
US4575830A Expired / Lapsed US ⏱️ 0.0 years remaining

Indirect shearwave determination

Assignee: Schlumberger Technology Corp
•
Filing: 1982-10-15
•
Granted: 1986-03-11
•
Expires: 2002-10-15
📄 Original PDF 🌐 Google Patents 🔄 Sync Now

Abstract

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.

Inventors

  • 👤 John D. Ingram
  • 👤 Josephine L. Murray

Bibliographic Details

Application Number —
Jurisdiction US
Data Source google_patents
Last Synchronized 2026-10-24 21:46
CPC Classifications
G; G01; G01H; G01H5/00; G01V; G01V1/00; G01V1/40; G01V1/44; G01V1/48; G01V1/50

Claims Analysis

Total Claims: 17
Claim 1
Independent Claim 289 words
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 μ 2 and λ 1 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 β 2 in accordance with the expression ##EQU6## where ρ 2 is formation density and ρ 1 is mud density, and (m) plotting the estimate of β 2 as a function of depth.
Claim 2
Independent Claim 289 words
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 μ 2 and λ 1 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 β 2 in accordance with the expression ##EQU6## where ρ 2 is formation density and ρ 1 is mud density, and (m) plotting the estimate of β 2 as a function of depth.
Claim 3
Dependent Claim 14 words
The method of claim 1 wherein the selected guided wave is a Stoneley wave.
Claim 4
Dependent Claim 12 words
The method of claim 1 wherein the window is a Blackman window.
Claim 5
Dependent Claim 15 words
The method of claim 3 wherein the Blackman window has a duration of 1200 microseconds.
Claim 6
Dependent Claim 24 words
The method of claim 1 including the step of time shifting the windowed energy prior to Fourier transformation to compensate for large phase differences.
Claim 7
Dependent Claim 24 words
The method of claim 2 including the step of time shifting the windowed energy prior to Fourier transformation to compensate for large phase differences.
Claim 8
Dependent Claim 24 words
The method of claim 3 including the step of time shifting the windowed energy prior to Fourier transformation to compensate for large phase differences.
Claim 9
Dependent Claim 21 words
The method of claim 5 wherein the amount of time shift introduced to the windowed energy is removed following Fourier transformation.
Claim 10
Dependent Claim 21 words
The method of claim 6 wherein the amount of time shift introduced to the windowed energy is removed following Fourier transformation.
Claim 11
Dependent Claim 21 words
The method of claim 7 wherein the amount of time shift introduced to the windowed energy is removed following Fourier transformation.
Claim 12
Dependent Claim 22 words
The method of claim 1 wherein the mud velocity α, is assumed in a range of 5555 ft/sec to 5000 ft sec.
Claim 13
Dependent Claim 19 words
The method of claim 1 wherein the mud velocity α, is approximated from a measured value of mud density.
Claim 14
Dependent Claim 16 words
The method of claim 1 wherein the ratio ρ 2 /ρ 1 is assumed a constant.
Claim 15
Dependent Claim 14 words
The method of claim 13 wherein the assumed value of the ratio is 2.5.
Claim 16
Independent Claim 306 words
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 μ 2 and λ 1 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 β 2 in accordance with the expression ##EQU8## where ρ 2 is formation density and ρ 1 is mud density, and (l) means for plotting the estimate of ρ 2 as a function of depth.
Claim 17
Independent Claim 306 words
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 μ 2 and λ 1 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 β 2 in accordance with the expression ##EQU8## where ρ 2 is formation density and ρ 1 is mud density, and (l) means for plotting the estimate of ρ 2 as a function of depth.

References Cited (Backward Citations)

US4210965A Cited
US4210966A Cited

USPTO Annuity & Maintenance Windows

Annuity Stage Window Opens Payment Deadline Surcharge Period Until Estimated Fee Status
3.5 Years (Year 4) 1989-03-10 1989-09-09 1990-03-11 $2000.00 Window Closed / Lapsed
7.5 Years (Year 8) 1993-03-10 1993-09-09 1994-03-11 $3760.00 Window Closed / Lapsed
11.5 Years (Year 12) 1997-03-10 1997-09-09 1998-03-11 $7700.00 Window Closed / Lapsed

Assigned Categories & Tags

No categories assigned to this patent yet.

Assign Category Tag