AVA compliant pre-stack frequency spectrum enhancement of seismic data

Patent metadata

Field	Value
Publication number	US-10191168-B2
Application number	US-201715402830-A
Country	US
Kind code	B2
Filing date	Jan 10, 2017
Priority date	Jan 12, 2016
Publication date	Jan 29, 2019
Grant date	Jan 29, 2019

Abstract

Generating spectrally enhanced seismic data expresses seismic data as a convolution of reflectivity and a seismic source wavelet. This seismic source wavelet varies over a sampling interval and defining a total amount of energy over the sampling interval. An enhanced seismic source wavelet that is a single-valued energy spike that yields the total amount of energy over the sampling interval is generated. In addition, the reflectivity is modified to preserve amplitude variation with angle. The reflectivity is convoluted with the enhanced seismic source wavelet and residual energy is added to the convolution to generate the spectrally enhanced seismic data.

First claim

What is claimed is:

1. A method for generating spectrally enhanced seismic data, the method comprising:
obtaining seismic data for a subsurface;
modelling the seismic data as linear function of a reflectivity and a seismic source wavelet, the seismic source wavelet varying over a sampling interval;
determining an enhanced seismic source wavelet comprising a single-valued spike that preserves wavelet energy; and
convoluting the reflectivity with the enhanced seismic source wavelet to generate spectrally enhanced seismic data,
wherein the modelling the seismic data is constrained by an amplitude variation with reflection angle.

2. The method of claim 1 , wherein the method further comprises compensating the seismic data for subsurface absorption using a sparse inversion.

3. The method of claim 1 , wherein the method further comprises solving the expressed function of reflectivity and the seismic source wavelet using non-stationary estimates of time-dependent and angle-dependent wavelets and a sparse inversion algorithm to generate an estimated reflectivity.

4. The method of claim 1 , wherein the method further comprises including linearized amplitude variation with angle coefficients in the reflectivity to preserve amplitude variation with angle in the spectrally enhanced seismic data.

5. The method of claim 1 , wherein the seismic source wavelet comprises a time-dependent convolutional wavelet matrix comprising a plurality of time-dependent and angle-dependent wavelets.

6. The method of claim 1 , wherein the determining of the enhanced seismic source wavelet comprises:
calculating an area under the seismic source wavelet over the sampling interval; and
identifying the enhanced seismic source wavelet to enclose the area over the sampling interval.

7. The method of claim 1 , wherein convoluting the reflectivity with the enhanced seismic source wavelet to generate spectrally enhanced seismic data further comprises adding residual energy representing noise in the seismic data to the convolution of the reflectivity and the enhanced seismic wavelet to generate spectrally enhanced seismic data.

8. A non-transitory computer-readable medium containing computer-executable code that when read by a computer causes the computer to perform a method for generating spectrally enhanced seismic data, the method comprising:
obtaining seismic data for a subsurface;
modelling the seismic data as linear function of a reflectivity and a seismic source wavelet, the seismic source wavelet varying over a sampling interval;
determining an enhanced seismic source wavelet comprising a single-valued spike that preserves wavelet energy; and
convoluting the reflectivity with the enhanced seismic source wavelet to generate spectrally enhanced seismic data,
wherein the modelling the seismic data is constrained by an amplitude variation with reflection angle.

9. The non-transitory computer-readable medium of claim 8 , wherein the method further comprises solving the expressed function of reflectivity and the seismic source wavelet using non-stationary estimates of time-dependent and angle-dependent wavelets and a sparse inversion algorithm to generate an estimated reflectivity.

10. The non-transitory computer-readable medium of claim 8 , wherein the method further comprises including linearized amplitude variation with angle coefficients in the reflectivity to preserve amplitude variation with angle in the spectrally enhanced seismic data.

11. The non-transitory computer-readable medium of claim 8 , wherein the method further comprises modifying reflectivity in accordance with an applied structural consistency in the subsurface to reduce noise in the spectrally enhanced seismic data.

12. The non-transitory computer-readable medium of claim 8 , wherein the method further comprises compensating the seismic data for subsurface absorption using a sparse inversion.

13. The non-transitory computer-readable medium of claim 8 , wherein the determining of the enhanced seismic source wavelet comprises:
calculating an area under the seismic source wavelet over the sampling interval; and
identifying the enhanced seismic source wavelet to enclose the area over the sampling interval.

14. The non-transitory computer-readable medium of claim 8 wherein convoluting the reflectivity with the enhanced seismic source wavelet to generate spectrally enhanced seismic data further comprises adding residual energy representing noise in the seismic data to the convolution of the reflectivity and the enhanced seismic wavelet to generate spectrally enhanced seismic data.

15. A computing system for generating spectrally enhanced seismic data, the computing system comprising:
a storage device comprising recorded seismic data for a subsurface; and
a processor in communication with the storage device and configured to:
model the seismic data as linear function of a reflectivity and a seismic source wavelet, the seismic source wavelet varying over a sampling interval;
determine an enhanced seismic source wavelet comprising a single-valued spike that preserves wavelet energy; and
convolute the reflectivity with the enhanced seismic source wavelet to generate spectrally enhanced seismic data,

wherein modelling the seismic data is constrained by an amplitude variation with reflection angle.

16. The computing system of claim 15 , wherein the processor is further configured to include linearized amplitude variation with angle coefficients in the reflectivity to preserve amplitude variation with angle in the spectrally enhanced seismic data.

17. The computing system of claim 15 , wherein the processor is further configured to modify reflectivity in accordance with an applied structural consistency in the subsurface to reduce noise in the spectrally enhanced seismic data.

18. The computing system of claim 15 , wherein the seismic source wavelet comprises a time-dependent convolutional wavelet matrix comprising a plurality of time-dependent and angle-dependent wavelets.

19. The computing system of claim 15 , wherein the processor in determining the enhanced seismic source wavelet:
calculates an area under the seismic source wavelet over the sampling interval; and
identifies the enhanced seismic source wavelet to enclose the area over the sampling interval.

20. The computing system of claim 15 , wherein the processor after convoluting the reflectivity with the enhanced seismic source wavelet to generate spectrally enhanced seismic data adds residual energy representing noise in the seismic data to the convolution of the reflectivity and the enhanced seismic wavelet to generate spectrally enhanced seismic data.

Assignees

• Cgg Services Sas

Inventors

• Peng Can
• Jafargandomi Arash
• Hoeber Henning
• Cooper James
• Poole Gordon

Classifications

• G01V2210/632

  Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI] · CPC title

• G01V1/362 Primary

  Effecting static or dynamic corrections; Stacking · CPC title

• G01V1/375 Primary

  Correlating received seismic signals with the emitted source signal · CPC title

• G01V1/307

  for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity · CPC title

• G01V1/282

  Application of seismic models, synthetic seismograms · CPC title