Suppressing noises in seismic data

1. A computer-implemented method for suppressing noises in seismic data, comprising:

receiving, at a data processing apparatus, a set of seismic data associated with a subsurface region;

flattening, by the data processing apparatus, the received set of seismic data according to an identified seismic event, wherein the identified seismic event is identified according to an interpretation of the subsurface region, and the identified seismic event comprises a fault, and wherein the flattening comprises:

selecting an anchor trace; and

squeezing or stretching other traces according to the anchor trace;

dividing, by the data processing apparatus, the flattened set of seismic data into a plurality of spatial windows according to spatial distances of the seismic data in a horizontal direction, wherein each of the plurality of spatial windows corresponds to a different horizontal distance index;

randomizing, by the data processing apparatus, the windowed set of seismic data according to a random sequential order, wherein the randomizing comprises:

assigning a new horizontal distance index to each of the plurality of spatial windows according to the random sequential order; and

reordering positions of the plurality of spatial windows according to the new horizontal distance indices;

filtering, by the data processing apparatus, the randomized seismic data, wherein a median filtering window size used in the filtering is smaller than a randomization window size used in the randomizing;

reorganizing, by the data processing apparatus, the filtered seismic data according to a pre-randomization order; and

generating, by the data processing apparatus, an image of the subsurface region based on the reorganized data.

2. The method of claim 1 , wherein the seismic data is post-stack data.

3. The method of claim 2 , wherein the seismic data includes at least one of pre-stack time migration results or pre-stack depth migration results.

4. The method of claim 1 , further comprising: subtracting, by the data processing apparatus, the set of reorganized seismic data from the seismic data.

5. The method of claim 1 , further comprising: filtering, by the data processing apparatus, the reorganized seismic data.

6. A non-transitory computer-readable medium storing instructions which, when executed, cause a computing device to perform operations comprising:

receiving, at a data processing apparatus, a set of seismic data associated with a subsurface region;

flattening, by the data processing apparatus, the received set of seismic data according to an identified seismic event, wherein the identified seismic event is identified according to an interpretation of the subsurface region, and the identified seismic event comprises a fault, and wherein the flattening comprises:

selecting an anchor trace; and

squeezing or stretching other traces according to the anchor trace;

dividing, by the data processing apparatus, the flattened set of seismic data into a plurality of spatial windows according to spatial distances of the seismic data in a horizontal direction, wherein each of the plurality of spatial windows corresponds to a different horizontal distance index;

randomizing, by the data processing apparatus, the windowed set of seismic data according to a random sequential order, wherein the randomizing comprises:

assigning a new horizontal distance index to each of the plurality of spatial windows according to the random sequential order; and

reordering positions of the plurality of spatial windows according to the new horizontal distance indices;

filtering, by the data processing apparatus, the randomized seismic data, wherein a median filtering window size used in the filtering is smaller than a randomization window size used in the randomizing;

reorganizing, by the data processing apparatus, the filtered seismic data according to a pre-randomization order; and

generating, by the data processing apparatus, an image of the subsurface region based on the reorganized data.

7. The non-transitory computer-readable medium of claim 6 , wherein the seismic data is post-stack data.

8. The non-transitory computer-readable medium of claim 7 , wherein the seismic data includes at least one of pre-stack time migration results or pre-stack depth migration results.

9. The non-transitory computer-readable medium of claim 6 , the operations further comprising: subtracting, by the data processing apparatus, the set of reorganized seismic data from the seismic data.

10. The non-transitory computer-readable medium of claim 6 , the operations further comprising: filtering, by the data processing apparatus, the reorganized seismic data.

11. A device, comprising:

at least one hardware processor; and

a non-transitory computer-readable storage medium coupled to the at least one hardware processor and storing programming instructions for execution by the at least one hardware processor, wherein the programming instructions, when executed, cause the at least one hardware processor to perform operations comprising:

receiving, at the least one hardware processor, a set of seismic data associated with a subsurface region;

flattening, by the least one hardware processor, the received set of seismic data according to an identified seismic event, wherein the identified seismic event is identified according to an interpretation of the subsurface region, and the identified seismic event comprises a fault, and wherein the flattening comprises:

selecting an anchor trace; and

squeezing or stretching other traces according to the anchor trace;

dividing, by the least one hardware processor, the flattened set of seismic data into a plurality of spatial windows according to spatial distances of the seismic data in a horizontal direction, wherein each of the plurality of spatial windows corresponds to a different horizontal distance index;

randomizing, by the least one hardware processor, the windowed set of seismic data according to a random sequential order, wherein the randomizing comprises:

assigning a new horizontal distance index to each of the plurality of spatial windows according to the random sequential order; and

reordering positions of the plurality of spatial windows according to the new horizontal distance indices;

filtering, by the least one hardware processor, the randomized seismic data, wherein a median filtering window size used in the filtering is smaller than a randomization window size used in the randomizing;

reorganizing, by the least one hardware processor, the filtered seismic data according to a pre-randomization order; and

generating, by the least one hardware processor, an image of the subsurface region based on the reorganized data.

12. The device of claim 11 , wherein the seismic data is post-stack data.

13. The device of claim 12 , wherein the seismic data includes at least one of pre-stack time migration results or pre-stack depth migration results.

14. The device of claim 11 , the operations further comprising: subtracting the set of reorganized seismic data from the seismic data.