A method for generating a geophysical image of a subsurface region includes defining a computational sub-volume for the geophysical image including a predetermined number of seismic traces of a plurality of seismic traces and a predetermined number of samples per each one of the plurality of seismic traces, generating a data matrix corresponding to a first sub-volume of the subsurface region based on the defined computational sub-volume, the data matrix comprising the predetermined number of samples for the predetermined number of traces of a portion of a seismic dataset corresponding to the first sub-volume. The method also includes estimating a coherence between the predetermined number of traces of the data matrix by performing a sum of a variance of the predetermined number of samples of the data matrix, and assigning the estimated coherence to a location in the geophysical image.

Seismic acquisition having high geophone densities is compressed based on Functional Quantization (FQ) for an infinite dimensional space. Using FQ, the entire sample path of the seismic waveform in a target function space is quantized. An efficient solution for the construction of a functional quantizer is given. It is based on Monte-Carlo simulation to circumvent the limitations of high dimensionality and avoids explicit construction of Voronoi regions to tessellate the function space of interest. The FQ architecture is then augmented with three different Vector Quantization (VQ) techniques which yield hybridized FQ strategies of 1) FQ-Classified VQ, 2) FQ-Residual/Multistage VQ and 3) FQ-Recursive VQ. Joint quantizers are obtained by replacing regular VQ codebooks in these hybrid quantizers by their FQ equivalents. Simulation results show that the FQ combined with any one of the different VQ techniques yields improved rate-distortion compared to either FQ or VQ techniques alone.

A method of compressing data from seismic waves using Gabor frames utilizes a plurality of geophones positioned within a region of interest. Each of the plurality of geophones is communicably coupled with at least one remote server. Thus, a plurality of reflected-seismic signals received through the plurality of geophones can be transmitted to the at least one remote server for analyzing and calculations. The plurality of reflected-seismic signals is converted into a set of Gabor frames, wherein the Gabor frames is generated via a plurality of prolate spheroidal wave functions (PSWF). A Gabor frame-generating calculation module utilizes the plurality of PSWF to generate the set of Gabor frames. A dual frame for each of the set of Gabor frames is derived and used for quantization purposes. Preferably, a tree structured vector quantization process is followed.

A system for seismic surveying, method for performing seismic and a non-transitory computer readable medium having instructions stored therein that, when executed by one or more processors, cause the one or more processors to perform a method for performing seismic surveying including emitting seismic waves into a substrate, receiving seismic waves reflected from discontinuities within the substrate, converting the seismic waves into seismic traces, and representing the seismic traces by superimposed multiple tone sinusoidal waves using a parameter estimation. An optimized residual of the modelling is compressed using entropy coding or quantization coding techniques, and the optimized residual and the parameter sets are transmitted to a remote processing station for reconstruction and analysis of the discontinuities.

The present method predicts and separates dispersive surface waves from seismic data using dispersion estimation and is completely data-driven and computer automated and no human intervention is needed. The method is capable of predicting and suppressing surface waves from recorded seismic data without damaging the reflections. Nonlinear signal comparison (NLSC) is used to obtain a high resolution and accurate dispersion. Based on the dispersion, surface waves are predicted from the field recorded seismic data. The predicted surface waves are then subtracted from the original data.

A method of compressing data from seismic waves using Gabor frames utilizes a plurality of geophones positioned within a region of interest. Each of the plurality of geophones is communicably coupled with at least one remote server. Thus, a plurality of reflected-seismic signals received through the plurality of geophones can be transmitted to the at least one remote server for analyzing and calculations. The plurality of reflected-seismic signals is converted into a set of Gabor frames, wherein the Gabor frames is generated via a plurality of prolate spheroidal wave functions (PSWF). A Gabor frame-generating calculation module utilizes the plurality of PSWF to generate the set of Gabor frames. A dual frame for each of the set of Gabor frames is derived and used for quantization purposes. Preferably, a tree structured vector quantization process is followed.

A system for seismic surveying, method for performing seismic and a non-transitory computer readable medium having instructions stored therein that, when executed by one or more processors, cause the one or more processors to perform a method for performing seismic surveying including emitting seismic waves into a substrate, receiving seismic waves reflected from discontinuities within the substrate, converting the seismic waves into seismic traces, and representing the seismic traces by superimposed multiple tone sinusoidal waves using a parameter estimation. An optimized residual of the modelling is compressed using entropy coding or quantization coding techniques, and the optimized residual and the parameter sets are transmitted to a remote processing station for reconstruction and analysis of the discontinuities.

Seismic acquisition having high geophone densities is compressed based on Functional Quantization (FQ) for an infinite dimensional space. Using FQ, the entire sample path of the seismic waveform in a target function space is quantized. An efficient solution for the construction of a functional quantizer is given. It is based on Monte-Carlo simulation to circumvent the limitations of high dimensionality and avoids explicit construction of Voronoi regions to tessellate the function space of interest. The FQ architecture is then augmented with three different Vector Quantization (VQ) techniques which yield hybridized FQ strategies of 1) FQ-Classified VQ, 2) FQ-Residual/Multistage VQ and 3) FQ-Recursive VQ. Joint quantizers are obtained by replacing regular VQ codebooks in these hybrid quantizers by their FQ equivalents. Simulation results show that the FQ combined with any one of the different VQ techniques yields improved rate-distortion compared to either FQ or VQ techniques alone.

A method is described for seismic imaging that may include receiving digital seismic data; processing the digital seismic data to create a digital seismic image in a seismic domain; flattening the digital seismic image to generate a digital flattened image; identifying artifacts in the digital flattened image; transforming the artifacts back into the seismic domain; and reprocessing the digital seismic data based on the artifacts in the seismic domain to generate a digital image with reduced artifacts. The method may be executed by a computer system.

A method for generating a geophysical image of a subsurface region includes defining a computational sub-volume for the geophysical image including a predetermined number of seismic traces of a plurality of seismic traces and a predetermined number of samples per each one of the plurality of seismic traces, generating a data matrix corresponding to a first sub-volume of the subsurface region based on the defined computational sub-volume, the data matrix comprising the predetermined number of samples for the predetermined number of traces of a portion of a seismic dataset corresponding to the first sub-volume. The method also includes estimating a coherence between the predetermined number of traces of the data matrix by performing a sum of a variance of the predetermined number of samples of the data matrix, and assigning the estimated coherence to a location in the geophysical image.