US5349527A
Expired / Lapsed
US
⏱️ 0.0 years remaining
Method of seismic time migration using a massively parallel computer
Assignee: Schlumberger Technology Corp
•
Filing: 1991-12-20
•
Granted: 1994-09-20
•
Expires: 2011-12-20
Abstract
Abstract A method for 3 dimensional, one pass migration of post stack seismic data is implemented on a massively parallel computer. The quality of the migration and speed of execution makes possible iterative 3D post stack depth migrations. The depth migration method computes in the frequency domain and downward continues one or more frequency planes through all depth levels of interest. For a single frequency, at each depth level, the method extrapolates a shifted position for each x,y position by applying a filter for 2D convolution. A processing component is assigned to each subgrid of x,y spatial positions and the processing components operate concurrently in parallel to determine filter coefficients independently for each x,y spatial position and to extrapolate the x,y shifted positions. The filter coefficients are derived independently at each x,y position by storing a limited table of filter coefficients in local memory of each processing component.
Inventors
- 👤 Andrew Pieprzak
- 👤 Peter T. Highnam
Bibliographic Details
Application Number
—
Jurisdiction
US
Data Source
google_patents
Last Synchronized
2026-11-07 15:19
CPC Classifications
G; G01; G01V; G01V1/00; G01V1/28; G01V1/30; G01V1/301
Claims Analysis
Total Claims: 22
Claim 1
Independent Claim
196 words
A method for time migration of space-time domain seismic data to build a reflector time map using a parallel computer having a number of processing components, comprising the steps of: (1) transforming the seismic data from the space-time domain x,y,t to a space frequency domain x,y,f data cube using a Fourier transform process; (2) storing the space frequency data cube in computer memory; (3) grouping one or more frequency planes from said space frequency data cube into frequency chunks; (4) downwardly continuing the space frequency data cube for each frequency chunk through a number of time steps to obtain a reflector time map, where for each time step the method includes the substeps of (a) modeling a velocity profile, (b) loading a table of operator coefficients into memory, (c) for a single frequency plane within the frequency chunk, assigning processing components to respective subgrids of one or more x,y spatial positions in the x,y plane, (d) extrapolating all of the x,y spatial positions by concurrently operating several processing components in parallel, for each x,y position the extrapolation using the velocity profile and some of the operator coefficients from the table, (e) building the reflector time map.
Claim 2
Independent Claim
196 words
A method for time migration of space-time domain seismic data to build a reflector time map using a parallel computer having a number of processing components, comprising the steps of: (1) transforming the seismic data from the space-time domain x,y,t to a space frequency domain x,y,f data cube using a Fourier transform process; (2) storing the space frequency data cube in computer memory; (3) grouping one or more frequency planes from said space frequency data cube into frequency chunks; (4) downwardly continuing the space frequency data cube for each frequency chunk through a number of time steps to obtain a reflector time map, where for each time step the method includes the substeps of (a) modeling a velocity profile, (b) loading a table of operator coefficients into memory, (c) for a single frequency plane within the frequency chunk, assigning processing components to respective subgrids of one or more x,y spatial positions in the x,y plane, (d) extrapolating all of the x,y spatial positions by concurrently operating several processing components in parallel, for each x,y position the extrapolation using the velocity profile and some of the operator coefficients from the table, (e) building the reflector time map.
Claim 3
Dependent Claim
15 words
The method of claim 1, wherein the migration is performed using post stack seismic data.
Claim 4
Dependent Claim
25 words
The method of claim 1, wherein the downward continuation step includes the substep of controlling the processing components using a single instruction multiple data method.
Claim 5
Dependent Claim
29 words
The method of claim 1, the modeling a velocity profile substep (4)(a) comprising reading the velocity plane for the depth of interest into local memory of each processing component.
Claim 6
Dependent Claim
14 words
The method of claim 1, wherein the extrapolating substep includes applying a 2D filter.
Claim 7
Dependent Claim
28 words
The method of claim 5, wherein the 2D filter having a number of filter coefficients is expressed as a sum of fixed size convolution filters using Chebyshev recursion.
Claim 8
Dependent Claim
37 words
The method of claim 5, wherein the 2D filter is in part described by a number of filter coefficients retrieved from a table of filter coefficients, the table being stored in local memory of each processing component.
Claim 9
Dependent Claim
27 words
The method of claim 7, wherein an entry in a filter coefficient table is retrieved by indexing the table with a velocity value and a frequency value.
Claim 10
Dependent Claim
26 words
The method of claim 8, wherein each processing component includes a plurality of processors, and each processor can access the stored table with different index values.
Claim 11
Dependent Claim
17 words
The method of claim 5, wherein the filter coefficients are determined independently for each x,y spatial position.
Claim 12
Dependent Claim
18 words
The method of claim 5, wherein the 2D filter is in part described by a symmetrical operator G.
Claim 13
Dependent Claim
27 words
The method of claim 11, wherein the operator G is a 5×5 mask with some of the mask entries variant as a function of frequency and velocity.
Claim 14
Dependent Claim
20 words
The method of claim 11, wherein the operator G is a 5×5 mask with constant mask entries for all frequencies.
Claim 15
Independent Claim
160 words
A method for migrating seismic data using a computer system having a number of processing components operable in parallel, where the seismic data cube is in the space frequency domain x,y,f and is stored in computer memory, characterized by: grouping one or more frequency planes from the seismic data cube into frequency chunks; and downwardly continuing for each frequency chunk through a number of time steps to obtain a seismic reflector time map where for each step the method includes the substeps of (a) for a single frequency plane within the frequency chunk, assigning processing components to respective subgrids of one or more x,y spatial positions, (b) determining a velocity profile for respective subgrids, and (c) extrapolating data at all x,y spatial positions in the subgrids by concurrently operating several processing components in parallel, using said velocity profile and frequency associated with said frequency plane to determine the filter coefficients of a 2D filter applied to said single frequency plane.
Claim 16
Independent Claim
160 words
A method for migrating seismic data using a computer system having a number of processing components operable in parallel, where the seismic data cube is in the space frequency domain x,y,f and is stored in computer memory, characterized by: grouping one or more frequency planes from the seismic data cube into frequency chunks; and downwardly continuing for each frequency chunk through a number of time steps to obtain a seismic reflector time map where for each step the method includes the substeps of (a) for a single frequency plane within the frequency chunk, assigning processing components to respective subgrids of one or more x,y spatial positions, (b) determining a velocity profile for respective subgrids, and (c) extrapolating data at all x,y spatial positions in the subgrids by concurrently operating several processing components in parallel, using said velocity profile and frequency associated with said frequency plane to determine the filter coefficients of a 2D filter applied to said single frequency plane.
Claim 17
Dependent Claim
37 words
The method of claim 14, wherein the 2D filter is in part described by a number of filter coefficients retrieved from a table of filter coefficients, the table being stored in local memory of each processing component.
Claim 18
Dependent Claim
17 words
The method of claim 14, wherein the filter coefficients are determined independently for each x,y spatial position.
Claim 19
Dependent Claim
29 words
The method of claim 14, wherein the downward continuation step includes the substep of updating said seismic reflector time map with said extrapolated data at all x,y spatial positions.
Claim 20
Dependent Claim
24 words
The method of claim 14, wherein the downward continuation step is performed concurrently for a plurality of single frequency planes within said frequency chunk.
Claim 21
Independent Claim
112 words
A method of migrating a space-time seismic data cube to a space-time reflector map comprising the steps of: transforming the space-time data cube into a space-frequency data cube comprising a plurality of frequency planes; extrapolating some of the frequency planes of the space-frequency data cube through the time levels of interest by assigning for a frequency plane processing components of a computer system to subgrids of one or more x,y spatial positions in the frequency plane, applying a frequency and velocity dependent 2D filter by concurrently operating the processing components in parallel to obtain a wavefield extrapolation; and creating a space-time reflector map using the wavefield extrapolations for said some frequency planes.
Claim 22
Independent Claim
112 words
A method of migrating a space-time seismic data cube to a space-time reflector map comprising the steps of: transforming the space-time data cube into a space-frequency data cube comprising a plurality of frequency planes; extrapolating some of the frequency planes of the space-frequency data cube through the time levels of interest by assigning for a frequency plane processing components of a computer system to subgrids of one or more x,y spatial positions in the frequency plane, applying a frequency and velocity dependent 2D filter by concurrently operating the processing components in parallel to obtain a wavefield extrapolation; and creating a space-time reflector map using the wavefield extrapolations for said some frequency planes.
References Cited (Backward Citations)
US4598400A
Cited
US4773038A
Cited
WO1988000711A2
Cited
USH482H
Cited
US4827403A
Cited
US5148406A
Cited
GB2236393A
Cited
US5150331A
Cited
US5128899A
Cited
USPTO Annuity & Maintenance Windows
| Annuity Stage | Window Opens | Payment Deadline | Surcharge Period Until | Estimated Fee | Status |
|---|---|---|---|---|---|
| 3.5 Years (Year 4) | 1997-09-19 | 1998-03-21 | 1998-09-20 | $2000.00 | Window Closed / Lapsed |
| 7.5 Years (Year 8) | 2001-09-19 | 2002-03-21 | 2002-09-20 | $3760.00 | Window Closed / Lapsed |
| 11.5 Years (Year 12) | 2005-09-19 | 2006-03-21 | 2006-09-20 | $7700.00 | Window Closed / Lapsed |