Coherent noise attenuation using statistical methods

• US 10,444,393 B2
• Filed: 03/27/2015
• Issued: 10/15/2019
• Est. Priority Date: 03/27/2015
• Status: Active Grant

• Current Assignee: Halliburton Energy Services Incorporated (Halliburton Company)
• Original Assignee: Halliburton Energy Services Incorporated (Halliburton Company)
• Inventors: Cheng, Dongjie, Laza, Laura Torres, Zhao, Xiaomin
• Primary Examiner(s): Ngon, Ricky
• Application Number:
  US15/301,629
• Publication Number:
  US 20170031047A1
• Time in Patent Office: 1,663 Days
• Field of Search: 702 6
• CPC Class Codes: G01V 1/325 Transforming one representation into another
• CPC Class Codes: G01V 1/364 Seismic filtering G01V1/37 takes precedence
• CPC Class Codes: G01V 1/48 Processing data
• CPC Class Codes: G01V 2210/121 Active source
• CPC Class Codes: G01V 2210/1295 Land surface
• CPC Class Codes: G01V 2210/1425 Land surface
• CPC Class Codes: G01V 2210/1429 Subsurface, e.g. in borehole or below weathering layer or mud line
• CPC Class Codes: G01V 2210/3246 Coherent noise, e.g. spatially coherent or predictable
• CPC Class Codes: G01V 2210/47 Slowness, e.g. tau-pi

Abstract

A system for attenuating coherent noise from seismic data comprises one or more sensors configured to sense waves generated by a seismic source and a coherent noise attenuation module communicably coupled to the one or more sensors and comprising a processor and memory. The coherent noise attenuation module is operable to receive a plurality of traces of seismic data from the one or more sensors and apply a first transformation to the plurality of traces, identify one or more outlier waveforms in the transformed traces, attenuate the identified outlier waveforms, and apply a second transformation to the plurality of traces that is the inverse of the first transformation.

Claims

1. A system for attenuating coherent noise from seismic data, comprising:

one or more sensors configured to sense waves generated by a seismic source; and

a coherent noise attenuation module comprising a processor and memory, the coherent noise attenuation module communicably coupled to the one or more sensors and operable to;

receive a plurality of traces of seismic data from the one or more sensors;

apply a first transformation to the plurality of traces of seismic data;

identify one or more outlier waveforms in the transformed traces;

attenuate the identified outlier waveforms; and

apply a second transformation to the plurality of transformed traces,

the second transformation being an inverse of the first transformation.

2. The system of claim 1, wherein

the coherent noise attenuation module operable to apply the first transformation is further operable to apply at least one of a time shift to each particular trace of the plurality of traces that is different from time shifts applied to each of the other traces of the plurality of traces and an inverse tau-p transformation to the plurality of traces of seismic data;

and

the coherent noise attenuation module operable to apply the second transformation is further operable to apply at least one of a time shift to each particular trace of the plurality of transformed traces that is the opposite of the time shift applied to the particular trace in the first transformation and an inverse tau-p transform to the plurality of transformed traces.

3. The system of claim 1, wherein the coherent noise attenuation module is further operable to isolate a portion of the plurality of traces of seismic data based on at least one of a time window and a moveout of a coherent noise event identified in the plurality of traces of seismic data.

4. The system of claim 1, wherein the coherent noise attenuation module operable to identify the one or more outlier waveforms is further operable to identify the one or more outlier waveforms using a threshold applied to at least a portion of the plurality of transformed traces.

5. The system of claim 4, wherein the threshold is based on at least one of mean and standard deviation of amplitudes of waveforms in the at least a portion of the plurality of transformed traces and percentiles of amplitudes of waveforms in the at least a portion of the plurality of transformed traces.

6. The system of claim 1, wherein the coherent noise attenuation module is further operable to determine the first transformation using an iterative estimation process.

7. The system of claim 1, wherein the coherent noise attenuation module operable to attenuate the identified outlier waveforms is further operable to at least one of multiply the waveforms by a coefficient that is less than one and replace each of the identified outlier waveforms with an interpolated waveform based on adjacent waveforms to the outlier waveform.

8. A method for attenuating coherent noise from seismic data, comprising:

receiving a plurality of traces of seismic data from one or more sensors configured to sense waves generated by a seismic source;

applying a first transformation to the plurality of traces of seismic data;

identifying one or more outlier waveforms in the transformed traces;

attenuating the identified outlier waveforms; and

applying a second transformation to the plurality of transformed traces, the second transformation being an inverse of the first transformation.

9. The method of claim 8, wherein:

applying the first transformation includes at least one of applying a time shift to each particular trace of the plurality of traces that is different from time shifts applied to each of the other traces of the plurality of traces and a tau-p transformation to the plurality of traces of seismic data; and

applying the second transformation includes applying at least one of a time shift to each particular trace of the plurality of transformed traces that is the opposite of the time shift applied to the particular trace in the first transformation and an inverse tau-p transform to the plurality of transformed traces.

10. The method of claim 8, further comprising isolating a portion of the plurality of traces of seismic data based on at least one of a time window and a moveout of a coherent noise event identified in the plurality of traces of seismic data.

11. The method of claim 8, wherein identifying the one or more outlier waveforms further comprises identifying the one or more outlier waveforms using a threshold applied to at least a portion of the plurality of transformed traces.

12. The method of claim 11, wherein the threshold is based on at least one of mean and standard deviation of amplitudes of waveforms in the at least a portion of the plurality of transformed traces and percentiles of amplitudes of waveforms in the at least a portion of the plurality of transformed traces.

13. The method of claim 8, further comprising determining the first transformation using an iterative estimation process.

14. The method of claim 8, wherein attenuating the identified outlier waveforms further comprises at least one of multiplying the waveforms by a coefficient that is less than one and replacing each of the identified outlier waveforms with an interpolated waveform based on adjacent waveforms to the outlier waveform.

15. A non-transitory computer-readable medium comprising storing one or more executable instructions that, when executed, by a processor, cause the processor to:

receive a plurality of traces of seismic data from one or more sensors configured to sense waves generated by a seismic source;

apply a first transformation to the plurality of traces of seismic data;

identify one or more outlier waveforms in the transformed traces;

attenuate the identified outlier waveforms; and

apply a second transformation to the plurality of transformed traces, the second transformation being an inverse of the first transformation.

16. The non-transitory computer-readable medium of claim 15, wherein:

applying the first transformation comprises applying at least one of a time shift to each particular trace of the plurality of traces that is different from time shifts applied to each of the other traces of the plurality of traces and a tau-p transformation to the plurality of traces of seismic data;

applying the second transformation comprises applying at least one of a time shift to each particular trace of the plurality of transformed traces that is the opposite of the time shift applied to the particular trace in the first transformation an inverse tau-p transform to the plurality of transformed traces; and

identifying the one or more outlier waveforms comprises identifying the one or more outlier waveforms using a threshold applied to at least a portion of the plurality of transformed traces.

17. The non-transitory computer-readable medium of claim 15, wherein the one or more instructions, when executed by the processor, further cause the processor to isolate a portion of the plurality of traces of seismic data based on at least one of a time window and a moveout of a coherent noise event identified in the plurality of traces of seismic data.

18. The non-transitory computer-readable medium of claim 15, wherein the threshold is based on at least one of mean and standard deviation of amplitudes of waveforms in the at least a portion of the plurality of transformed traces and percentiles of amplitudes of waveforms in the at least a portion of the plurality of transformed traces.

19. The non-transitory computer-readable medium of claim 15, wherein the one or more instructions, when executed by the processor, further cause the processor to determine the first transformation using an iterative estimation process.

20. The non-transitory computer-readable medium of claim 15, wherein the one or more instructions, when executed by the processor, further cause the processor to attenuate the identified outlier waveforms are further operable to cause the processor to at least one of multiply the waveforms by a coefficient that is less than one and replace each of the identified outlier waveforms with an interpolated waveform based on adjacent waveforms to the outlier waveform.