Wavefield reconstruction

• US 10,539,695 B2
• Filed: 04/28/2015
• Issued: 01/21/2020
• Est. Priority Date: 04/28/2014
• Status: Active Grant

• Current Assignee: WesternGeco LLC (Schlumberger NV)
• Original Assignee: WesternGeco LLC (Schlumberger NV)
• Inventors: Ozbek, Ali, Vassallo, Massimiliano, Halliday, David Fraser, Hopperstad, Jon-Fredrik, Laws, Robert Montgomery
• Primary Examiner(s): Barbee, Manuel L
• Assistant Examiner(s): Nimox, Raymond L
• Application Number:
  US15/307,742
• Publication Number:
  US 20170115415A1
• Time in Patent Office: 1,729 Days
• CPC Class Codes: G01V 1/282 Application of seismic models, synthetic seismograms
• CPC Class Codes: G01V 1/325 Transforming one representation into another
• CPC Class Codes: G01V 2210/42 Waveform, i.e. using raw or pre-filtered trace data

Abstract

A technique for reconstructing a seismic wavefield includes receiving data over one or more channels of a plurality of channels from a plurality of stations. The data is recorded by a plurality of seismic receivers and represent measurements of properties of the seismic wavefield. Each station includes a region in space including one or more seismic receivers. Each channel either measures a property of the seismic wavefield or a property of the seismic wavefield after the seismic wavefield has undergone a known transformation. At least one channel is derived as a function of one or more other channels. The technique includes using a processor based machine to process the data to model the seismic wavefield as a sum of basis functions; apply to the basis functions at least one forward transformation that describes the measurements received over the channel(s); and determine optimum basis functions based at least in part on the measurements.

Claims

1. A method for reconstructing a seismic wavefield, comprising:

receiving data over one or more channels of a plurality of channels from a plurality of stations, wherein;

the data are recorded by a plurality of seismic receivers and represent measurements of properties of the seismic wavefield;

each station of the plurality of stations comprises a region in space including one or more seismic receivers of the plurality of seismic receivers;

each channel of the plurality of channels either measures a property of the seismic wavefield or a property of the seismic wavefield after the seismic wavefield has undergone a known transformation; and

at least one channel of the plurality of channels is derived as a function of one or more other channels of the plurality of channels; and

using a processor based machine to process the data to;

model the seismic wavefield as a sum of basis functions;

apply to the basis functions forward transformations describing the measurements received over the one or more channels of the plurality of channels; and

determine optimum basis functions based at least in part on the measurements.

2. The method of claim 1 wherein at least one channel of the plurality of channels is derived as a function of at least two other channels of the plurality of channels.

3. The method of claim 1, wherein the task of processing data in the processor-based machine to determine the optimum basis functions comprises:

iteratively matching the transformed basis functions to the measurements received over the one or more channels; and

using the best matched transformed basis functions to reconstruct the seismic wavefield at one or more locations.

4. The method of claim 1, wherein the task of processing data in processor-based machine to determine the optimum basis functions comprises:

solving a sparsity-promoting optimization problem; and

using the transformed optimum basis functions to reconstruct the seismic wavefield at one or more locations.

5. The method of claim 1, wherein an average station interval associated with the stations in at least one spatial dimension is greater than half of the shortest wavelength of interest in the wavefield.

6. The method in claim 1, wherein the function comprises a finite-difference operation.

7. The method of claim 6, wherein the function comprises a finite-difference operation having an order of at least two.

8. The method of claim 7, wherein the finite-difference operation comprises a multi-dimensional finite-difference operation.

9. The method of claim 1, wherein the function comprises a finite difference operation, a given station has an associated average receiver interval within the given station over which finite-differences are computed according to the finite different operation is shorter than half of the shortest wavelength of interest in the wavefield.

10. A method for reconstructing a seismic wavefield, comprising:

receiving data from one or more stations in response to one or more seismic source arrays operating in one or more seismic source array configurations, wherein;

the data are recorded by one or more seismic receivers and represent measurements of properties of a seismic wavefield excited by the one or more seismic source arrays operating in the one or more seismic source array configurations;

each configuration of the plurality of seismic source array configurations is associated with a source location and is associated with a configuration of one or more seismic source elements; and

each seismic source array configuration of the plurality of seismic source array configurations emits energy having an associated radiation pattern; and

using a processor based machine to process the data to;

model the seismic wavefield as a sum of basis functions;

apply to the basis functions forward transformations describing the emitted radiation pattern associated with at least one source array configuration of the plurality of seismic source array configurations; and

determine optimum basis functions based at least in part on the measurements.

11. The method of claim 10, wherein the reconstructed seismic wavefield comprises the wavefield that would be received in response to one or more seismic source configurations with given radiation patterns from one or more locations.

12. The method of claim 11, wherein at least one of the given radiation patterns comprises a monopolar radiation pattern.

13. The method of claim 10, wherein at least one seismic source array configuration of the plurality of seismic source array configurations is derived as a function of one or more other seismic source array configurations.

14. The method of claim 10, wherein the task of processing data in processor-based machine to determine the optimum basis functions comprises:

iteratively matching the transformed basis functions to the measurements represented by the data in response to the one or more seismic source array configurations; and

using the matching to reconstruct the seismic wavefield emitted from the one or more source locations.

15. The method of claim 10, wherein the task of processing data in processor-based machine to determine the optimum basis functions comprises:

solving a sparsity-promoting optimization problem; and

using the transformed optimum basis functions to reconstruct the seismic wavefield emitted from the one or more source locations.

16. The method of claim 10, wherein an average source location interval associated with the source location in at least one spatial dimension is greater than half of the shortest wavelength of interest in the wavefield.

17. The method in claim 10, wherein the source array configuration comprises a source array comprising a subset of elements activated substantially in anti-phase.

18. The method in claim 10, wherein the source array configuration comprises a source array comprising elements activated in phase.

19. The method in claim 10, wherein the source array configuration at one or more locations comprises a source array comprising a subset of elements activated in anti-phase;

and another source array configuration at one or more locations comprises a source array comprising elements activated in phase.

20. The method of claim 10, wherein the source element comprises an airgun or a marine vibrator.

21. A method comprising:

receiving data over one or more channels of a plurality of channels from a plurality of stations, wherein;

the data are recorded by a plurality of receivers and represent measurements of properties of a signal to be reconstructed;

each station of the plurality of stations comprises a region in space including one or more seismic receivers of the plurality of receivers;

each channel of the plurality of channels either measures a property of a signal to be reconstructed or a property of the signal after the signal has undergone a known transformation; and

at least one channel of the plurality of channels is derived as a function of one or more other channels of the plurality of channels; and

using a processor based machine to process the data to;

model the signal as a sum of basis functions;

apply to the basis functions at least one forward transformation describing the measurements received over the one or more channels of the plurality of channels; and

determine optimum basis functions based at least in part on the measurements.

22. A system comprising:

an interface to receive data over one or more channels of a plurality of channels from a plurality of stations, wherein;

the data are recorded by a plurality of seismic receivers and represent measurements of properties of the seismic wavefield;

each station of the plurality of stations comprises a region in space including one or more seismic receivers of the plurality of seismic receivers;

each channel of the plurality of channels either measures a property of the seismic wavefield or a property of the seismic wavefield after the seismic wavefield has undergone a known transformation; and

at least one channel of the plurality of channels is derived as a function of one or more other channels of the plurality of channels; and

the processor processes the data to;

model the seismic wavefield as a sum of basis functions;

apply to the basis functions forward transformations describing the measurements received over the one or more channels of the plurality of channels; and

determine optimum basis functions based at least in part on the measurements.

23. A system comprising:

an interface to receive data from one or more stations in response to one or more seismic source arrays operating in one or more seismic source array configurations, wherein;

the data are recorded by one or more seismic receivers and represent measurements of properties of a seismic wavefield excited by the one or more seismic source arrays operating in the one or more seismic source array configurations;

each configuration of the plurality of seismic source array configurations is associated with a source location and is associated with a configuration of one or more seismic source elements; and

each seismic source array configuration of the plurality of seismic source array configurations emits energy having an associated radiation pattern; and

a processor to;

model a seismic wavefield to be reconstructed as a sum of basis functions;

apply to the basis functions at least one forward transformation describing the emitted radiation pattern associated with at least one source array configuration of the plurality of seismic source array configurations; and

determine optimum basis functions based at least in part on the measurements.