Method and apparatus for true relative amplitude correction of seismic data for normal moveout stretch effects

• US 7,230,879 B2
• Filed: 02/12/2005
• Issued: 06/12/2007
• Est. Priority Date: 02/12/2005
• Status: Expired due to Fees

• Current Assignee: Chevron USA Incorporated (Chevron Corp.)
• Original Assignee: Chevron USA Incorporated (Chevron Corp.)
• Inventors: Herkenhoff, E. Frederic, Alford, Richard B., Martin, Harry L.
• Primary Examiner(s): Lobo; Ian J.
• Application Number:
  US11/056,640
• Publication Number:
  US 20060193205A1
• Time in Patent Office: 850 Days
• Field of Search: 367/47, 367/52, 702/14, 702/17
• US Class Current: 367/47
• CPC Class Codes: G01V 1/362 Effecting static or dynamic corrections; Stacking
• CPC Class Codes: G01V 2210/52 Move-out correction
• CPC Class Codes: G01V 2210/522 Dip move-out [DMO]

Abstract

The present invention provides a method and apparatus for arriving at true relative amplitude destretched seismic traces from stretched seismic traces. The method compensates for offset varying reflection interference effects due to normal moveout. Stretch factors β and also input spectra are determined for NMOR stretched seismic traces. Estimates are then made of stretched wavelet spectra from the input spectra. A destretched wavelet spectra is then obtained. Shaping correction factors are determined by taking the ratio of the destretched wavelet spectra to the stretched wavelet spectra and are applied to the input spectra of the stretched traces to arrive at a destretched trace spectra. True relative amplitude scaling factors are computed by taking the ratio of a true relative amplitude property of the destretched wavelet spectra to a corresponding true relative amplitude property of the stretched wavelet spectra. Finally, the true relative amplitude scaling factors are applied to the destretched trace spectra to arrive at true relative amplitude destretched seismic traces.

Claims

1. A method for arriving at true relative amplitude destretched seismic traces from stretched seismic traces, the method comprising the steps of:

a. acquiring stretched seismic traces;

b. determining stretch factors β

for the stretched seismic traces;

c. determining input spectra for the stretched seismic traces;

d. determining estimates of stretched wavelet spectra from the input spectra;

e. obtaining a destretched wavelet spectra;

f. determining shaping correction factors by taking the ratio of the destretched wavelet spectra to the stretched wavelet spectra;

g. applying the shaping correction factors to the input spectra of the stretched traces to arrive at a destretched trace spectra;

h. computing true relative amplitude scaling factors by taking the ratio of a true relative amplitude property of the destretched wavelet spectra to a corresponding true relative amplitude property of the stretched wavelet spectra; and

i. applying the true relative amplitude scaling factors to the destretched trace spectra to arrive at true relative amplitude destretched seismic trace whereby the true relative amplitude property of the stretched wavelet spectra is substantially preserved.

2. The method of claim 1 wherein:

the destretched wavelet spectra is obtained by mapping the stretched wavelet spectra utilizing the similarity theorem with the stretch factor β

into a destretched wavelet spectra utilizing the following mathematical formulas;

W(f)=W_n(f/β

)/|β

|;

where;

W(f)=destretched wavelet spectra;

W_n(f)=the stretched wavelet spectra; and

β

=the wavelet stretch factor.

3. The method of claim 2 wherein the stretched wavelet spectra are modified to correct for a non-white reflectivity utilizing the following mathematical formula:

W_s(f)=W_n(f)R(f)/R(fβ

);

where;

W_s(f)=the stretched wavelet spectra corrected for a non-white reflectivity spectrum;

W_n(f)=the stretched wavelet spectra;

R(f)=a user specified reflectivity spectrum in normal incidence time; and

β

=the wavelet stretch factor.

4. The method of claim 1 wherein:

the destretched wavelet spectra is obtained from a user specified target wavelet spectra;

W(f)=W_d(f);

where;

W(f)=destretched wavelet spectra; and

W_d(f)=a user specified target wavelet spectra.

5. The method of claim 4 wherein the target wavelet spectra are modified to correct for a non-white reflectivity utilizing the following mathematical formula:

W(f)=W_d(f)R(f);

where;

W(f)=destretched wavelet spectra;

W_d(f)=a user specified target wavelet spectra; and

R(f)=a user specified reflectivity spectrum in vertical time.

6. The method of claim 1 wherein:

the stretch factors β

are calculated deterministically as functions of offset, time, rms velocity, and interval velocity.

7. The method of claim 1 wherein:

the stretch factors β

are calculated deterministically as functions of offset, time, rms velocity, interval velocity, and at least one of overburden anisotropy and geologic dip.

8. The method of claim 7 wherein:

the stretched seismic traces are CDP gathers and β

is calculated in accordance with the following mathematical expression;

β

=(1−

sin²θ

_r)<sup>−

1/2</sup>=1/cos θ

_r;

where;

sin²θ

_r=(xv_i cos θ

_d)²/(v_a⁴(t₀²+(x cos θ

_d)²/v_a²));

v_a=v_rms(1+2ε

sin²α

))^1/2;

sin²α

=x²/(x²+v_rmst₀)²); and

where;

θ

_r is the reflection angle at the interface;

x is shot to group offset;

v_i is the subsurface interval velocity at the t₀;

v_rms is the rms velocity of overburden;

t₀ is the zero offset travel time;

ε

is the overburden anisotropy parameter; and

θ

_d is the geologic dip at the interface.

9. The method of claim 7 wherein:

the stretched seismic traces are dip moveout (DMO) gathers and β

is calculated in accordance with the following mathematical expression;

β

=(1−

sin²θ

_r)<sup>−

1/2</sup>;

where;

sin²θ

_r=(x cos θ

_d)²/((t_dmo²v_a⁴/v_i²)+(cos²θ

_dv_a²/v_i²));

v_a=vrms(1+2ε

sin²α

(2−

sin²α

))^1/2;

sin²α

=x²/(x²+v_rmst₀)²); and

where;

θ

_r is the reflection angle at the interface;

x is shot to group offset,

v_i is the subsurface interval velocity at the t₀;

v_rms is the rms velocity of overburden;

t_dmo is the zero offset travel time of DMO gathers;

ε

is the overburden anisotropy parameter; and

θ

_d is the geologic dip at the interface.

10. The method of claim 7 wherein:

the stretched seismic traces are CRP gathers and β

is calculated in accordance with the following mathematical expression;

β

=(1−

sin²θ

_r)<sup>−

1/2</sup>;

where;

sin²θ

_r=(x cos²θ

_d)²/((t_mig²v_a⁴/v_i²)+(cos⁴θ

_dv_a²/v_i²))

v_a=vrms(1+2ε

sin²α

(2−

sin²α

))^1/2;

sin²α

=x²/(x²+v_rmst₀)²); and

where;

θ

is the reflection angle at the interface;

x is shot to group offset,

v_i is the subsurface interval velocity at the t₀;

v_rms is the rms velocity of overburden;

t_mig is the zero offset time of migrated gathers;

ε

is the overburden anisotropy parameter; and

θ

_d is the geologic dip at the interface.

11. The method of claim 1 wherein:

the true relative amplitude property, that is preserved in the true relative amplitude destretched traces, is selected from one of the following;

zero time value of the wavelet, the area under the amplitude spectrum of the wavelet, and time average mean absolute value of the stretched traces.

12. A program storage device readable by a machine, tangibly embodying a program of instructions executable by the machine to perform method steps for arriving at true relative amplitude destretched seismic traces from stretched seismic traces, the method comprising the steps of:

a acquiring stretched seismic traces;

b determining stretch factors β

for the stretched seismic traces;

c. determining input spectra for the stretched seismic traces;

d. determining estimates of stretched wavelet spectra from the input spectra;

e. obtaining a destretched wavelet spectra;

f. determining shaping correction factors by taking the ratio of the destretched wavelet spectra to the stretched wavelet spectra;

g. applying the shaping correction factors to the input spectra of the stretched traces to arrive at a destretched trace spectra;

h. computing true relative amplitude scaling factors by taking the ratio of a true relative amplitude property of the destretched wavelet spectra to a corresponding true relative amplitude property of the stretched wavelet spectra; and

i. applying the true relative amplitude scaling factors to the destretched trace spectra to arrive at true relative amplitude destretched seismic trace whereby the true relative amplitude property of the stretched wavelet spectra is substantially preserved.