Enhancing Interface Waves Filtering by Hamiltonian Monte-Carlo Sampling
AssigneeSaudi Arabian Oil Company
InventorMamadou Sanou Diallo
Life of this patent
Estimating the term…
Abstract
Filtering interface waves from seismic data using Hamiltonian Monte Carlo (HMC) sampling to produce filtered seismic data. A time frequency transform up to a maximum frequency is performed on single component or multicomponent seismic data. An attribute for distinguishing interface waves is transformed to produce a normalized attribute. HMC sampling is applied to sample the time-frequency coefficients associated most likely associated with the interface waves and identify the interface waves component. The interface wave component may be subtracted from the raw seismic trace using an adaptive subtraction algorithm.
Claims21
A computer-implemented method for producing filtered seismic data from seismic data having a seismic trace generated from a seismic receiver station configured to sense seismic signals originating from a seismic source station, the method comprising:
- obtaining seismic data from the seismic receiver station;
- performing a time-frequency transform of the seismic data up to a maximum frequency (f_max) to produce transformed seismic data, wherein the maximum frequency defines a maximum of a frequency range of interface wave arrivals in the seismic data;
- performing a transformation of an attribute of the transformed seismic data to produce a normalized attribute;
- performing, on the normalized attribute, a Hamiltonian Monte Carlo sampling of a parameter associated with interface waves to determine an interface wave component of the transformed seismic data;
- and
- adaptively subtracting the interface wave component from the seismic trace to produce the filtered seismic data.
A non-transitory computer-readable storage medium having executable code stored thereon for producing filtered seismic data from seismic data having a seismic trace generated from a seismic receiver station configured to sense seismic signals originating from a seismic source station, the executable code comprising a set of instructions that causes a processor to perform operations comprising:
- obtaining seismic data from the seismic receiver station;
- performing a time-frequency transform of the seismic data up to a maximum frequency (f_max) to produce transformed seismic data, wherein the maximum frequency defines a maximum of a frequency range of interface wave arrivals in the seismic data;
- performing a transformation of an attribute of the transformed seismic data to produce a normalized attribute;
- performing, on the normalized attribute, a Hamiltonian Monte Carlo sampling of a parameter associated with interface waves to determine an interface wave component of the transformed seismic data;
- and
- adaptively subtracting the interface wave component from the seismic trace to produce the filtered seismic data.
A system, comprising:
- a seismic source station;
- a seismic receiver station configured to sense seismic signals originating from a seismic source station;
- a seismic data processor;
- a non-transitory computer-readable storage memory accessible by the seismic data processor and having executable code stored thereon for producing filtered seismic data from seismic data having a seismic trace generated from the seismic receiver station, the executable code comprising a set of instructions that causes the seismic data processor to perform operations comprising:
- performing a time-frequency transform of the seismic data up to a maximum frequency (f_max) to produce transformed seismic data, wherein the maximum frequency defines a maximum of a frequency range of interface wave arrivals in the seismic data;
- performing a transformation of an attribute of the transformed seismic data to produce a normalized attribute;
- performing, on the normalized attribute, a Hamiltonian Monte Carlo sampling of a parameter associated with interface waves to determine an interface wave component of the transformed seismic data;
- and
- adaptively subtracting the interface wave component from the seismic trace to produce the filtered seismic data.
Dependencies are inferred from the claim text.