Method and system for limited frequency seismic imaging

• US 6,859,734 B2
• Filed: 10/06/2003
• Issued: 02/22/2005
• Est. Priority Date: 10/04/2002
• Status: Active Grant

• Current Assignee: Paradigm Geophysical Corporation (Emerson Electric Company)
• Original Assignee: Paradigm Geophysical Corporation (Emerson Electric Company)
• Inventors: Bednar, J. Bee
• Primary Examiner(s): MCELHENY JR, DONALD E
• Application Number:
  US10/679,923
• Publication Number:
  US 20040111217A1
• Time in Patent Office: 505 Days
• Field of Search: 702/14, 702/16, 702/17, 702/18
• US Class Current: 702/14
• CPC Class Codes: G01V 1/28 Processing seismic data, e.g. for interpretation or for event detection G01V1/48 takes precedence

Abstract

The method and system of the present inventions provides a new approach for selecting frequencies comprising a minimal or limited frequency data set for imaging leading to an order of magnitude decrease in computation resources. Seismic data are acquired and transformed a frequency domain. The data may be gathered in any manner. A plurality of minimal frequency subsets are formed comprising selected frequencies. The frequencies may be selected so that more frequencies are selected for shallower image positions. The selected frequencies are imaged to form frequency images. The selected frequencies may be individually weighted to assure a full and balanced bandwidth. The frequency images are composited to form an intermediate or final composited migration image.

Claims

1. A method for processing seismic data comprising:

a) acquiring seismic data;

b) transforming said seismic data to a frequency domain;

c) forming a plurality of minimal frequency subsets comprising selected frequencies;

d) imaging said selected frequencies to form frequency images; and

e) compositing said frequency images to form a composited migration image.

2. The method of claim 1 wherein forming said plurality of minimal frequency subsets further comprises selecting frequencies for the minimal frequency subsets according to f_n=α

(z)<sup>−

n</sup>f₀ for n=0 to n=N where $N=\frac{\log <br><br>\text{\hspace{1em}<br><br>}<br><br>f_0-\log <br><br>\text{\hspace{1em}<br><br>}<br><br>f_1}{\log <br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{\alpha <br><br>}<br><br>\left(z\right)}.$

3. The method of claim 1 further comprising reducing the number of frequencies as depth increases.

4. The method of claim 1 further comprising reducing the number of frequencies as depth increases according to $\alpha$

⁡

(z)=11+(khkz)2=11+tan2⁢

 

⁢

γ

⁡

(z).

5. The method of claim 1 further comprising weighting at least one of i) a plurality of said selected frequencies at a depth level and ii) a plurality of said imaged frequencies at a depth level.

6. A digital computer programmed to utilize seismic data traces obtained over a region of the earth'"'"'s subsurface to perform a process comprising:

(a) forming a gather of seismic data traces;

(b) transforming said seismic data to a frequency domain;

(c) forming a plurality of minimal frequency subsets comprising selected frequencies;

(d) imaging said selected frequencies to form frequency images; and

(e) compositing said frequency images to form a composited migration image.

7. The digital computer of claim 6 further programmed to perform a process for forming said plurality of minimal frequency subsets by selecting frequencies for the minimal frequency subsets according to f_n=α

(z)<sup>−

n</sup>f₀ for n=0 to n=N where $N=\frac{\log <br><br>\text{\hspace{1em}<br><br>}<br><br>f_0-\log <br><br>\text{\hspace{1em}<br><br>}<br><br>f_1}{\log <br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{\alpha <br><br>}<br><br>\left(z\right)}.$

8. The digital computer of claim 6 further programmed to perform a process for reducing the number of frequencies as depth increases.

9. The digital computer of claim 6 further programmed to perform a process for reducing the number of frequencies as depth increases according to $\alpha$

⁡

(z)=11+(khkz)2=11+tan2⁢

 

⁢

γ

⁡

(z).

10. The digital computer of claim 6 further programmed to perform a process for weighting at least one of i) a plurality of said selected frequencies at a depth level and ii) a plurality of said imaged frequencies at a depth level.

11. A method for processing seismic data comprising:

a) acquiring seismic data;

b) transforming said seismic data to a frequency domain;

c) forming a plurality of minimal frequency subsets comprising selected frequencies;

said selected frequencies in the range (f₀, f₁) for the minimal frequency subsets selected according to f_n=α

(z)<sup>−

n</sup>f₀ for n=0 ti n=N where $N=\frac{\log <br><br>\text{\hspace{1em}<br><br>}<br><br>f_0-\log <br><br>\text{\hspace{1em}<br><br>}<br><br>f_1}{\log <br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{\alpha <br><br>}<br><br>\left(z\right)}.$

d) imaging said selected frequencies to form frequency images; and

e) compositing said frequency images to form a composited migration image.

12. The system of claim 11 further comprising reducing the number of frequencies as depth increases.

13. The system of claim 11 further comprising reducing the number of frequencies as depth increases according to $\alpha$

⁡

(z)=11+(khkz)2=11+tan2⁢

 

⁢

γ

⁡

(z).

14. A digital computer programmed to utilize seismic data traces obtained over a region of the earth'"'"'s subsurface to perform a process comprising:

a) forming a gather of seismic data traces;

b) transforming said seismic data to a frequency domain;

c) forming a plurality of minimal frequency subsets comprising selected frequencies, said selected frequencies in the range (f₀, f₁) for the minimal frequency subsets selected according to f_n=α

(z)<sup>−

n</sup>f₀ for n=0 to n=N where $N=\frac{\log <br><br>\text{\hspace{1em}<br><br>}<br><br>f_0-\log <br><br>\text{\hspace{1em}<br><br>}<br><br>f_1}{\log <br><br>\text{\hspace{1em}<br><br>}<br><br>\mathrm{\alpha <br><br>}<br><br>\left(z\right)};<br><br>$

d) imaging said selected frequencies to form frequency images; and

e) compositing said frequency images to form a composited migration image.

15. The digital computer of claim 14 further programmed to perform a process for reducing die number of frequencies as depth increases.

16. The digital computer of claim 14 further programmed to perform a process for reducing the number of frequencies as depth increases according to $\alpha$

⁡

(z)=11+(khkz)2=11+tan2⁢

 

⁢

γ

⁡

(z).