Systems, methods, and apparatuses for generating a subsurface image using diffraction energy information are disclosed. The systems, methods, and apparatuses may include converting a shot gather into one or more plane-wave gather using a Radon transform. The plane-wave gathers may be extrapolated into source-side wavefields and receiver-side wavefields and further generate a pseudo dip-angle gather. The diffraction energy information may be extracted from the pseudo dip-angle gather, and an image containing subsurface features may be generated from the extracted diffraction energy information. The receiver-side wavefields may be decomposed using a recursive Radon transform.

Raw 3D seismic streamer wavefield data is received as a receiver-ghosted shot gather. The received receiver-ghosted shot gather shot gather is processed into a normalized form as normalized data. The normalized data is partitioned into a plurality of user-defined sub-gathers and processed to generate a complete receiver-deghosted shot gather. Output of the complete receiver-deghosted shot gather is initiated.

Computing device, computer instructions and method for processing input seismic data d. The method includes receiving the input seismic data d recorded in a data domain, solving a linear inversion problem constrained by input seismic data d to obtain a model domain and its energy, wherein the linear inversion problem is dependent on sparseness weights that are simultaneously a function of both time and frequency, reverse transforming the model domain energy to the data domain, and generating an image of a surveyed subsurface based on the reverse transformed model domain energy.

An example method for displaying sonic logging data associated with a formation surrounding a borehole can include acquiring sonic data at a plurality of depths using an acoustic array located in the borehole and transforming the acquired sonic data from a time-space domain to a frequency-wave number domain at a limited number of discrete frequencies. The method can also include estimating slowness values at the limited number of discrete frequencies from the transformed sonic data, interpolating the estimated slowness values to obtain a projection of one or more slowness-frequency dispersions of the acquired sonic data and displaying the projection of the slowness-frequency dispersions. The projection of the slowness-frequency dispersions can include a plurality of color bands corresponding to each of the limited number of discrete frequencies.

A method for deblending seismic signals includes entering as input to a computer recorded signals comprising seismic energy from a plurality of actuations of one or more seismic energy sources. A model of deblended seismic data and a blending matrix are initialized. A blending matrix inversion is performed using the initialized model. The inversion includes using a scaled objective function. The inversion is constrained by a thresholding operator. The thresholding operator is arranged to recover coefficients of the model of the deblended seismic data that are substantially nonzero, against a Gaussian white noise background. The thresholded model is projected into data space. Performing the blending matrix inversion is repeated if a data residual exceeds a selected threshold and the inversion is terminated if the data residual is below the selected threshold. At least one of storing and displaying an output of the blending matrix inversion is performed when the blending matrix inversion is terminated.

A method of generating geophysical data using at least one source. The method may include the steps of generating a geophysical wavefield with a varying signature using at least one source, wherein the signature is varied in a periodic pattern.

A zero-offset wavefield synthesis workflow to calculate a synthesized zero-offset wavefield output without the commitment to an rms velocity field output to circumvent velocity uncertainty. Said zero-offset wavefield synthesis workflow comprises calculating a migration cube output. Rendering a demigration cube output from said migration cube output with a demigration cube calculation. Rendering said synthesized zero-offset wavefield output from said demigration cube output with a zero-offset wavefield synthesis procedure.

A method of generating target-oriented acquisition-imprint-free prestack angle gathers using common focus point (CFP) operators includes receiving a plurality of seismic traces associated with a target point in a reservoir. A first angle domain common image gather (ADCIG) is generated based on the received plurality of seismic traces. A plurality of synthetic traces associated with the target point is generated. A second ADCIG is generated based on the synthetic traces. An enhanced ADCIG is generated using the first ADCIG and the second ADCIG.

Post-stack seismic data is received. Transformed seismic data is created from the received post-stack seismic data, including performing a super-resolution radon transform on the post-stack seismic data. Signal and noise regions are separated using the transformed seismic data, including using a defined muting function to remove unwanted noise. An inverse radon transform is performed using the separated signal and noise regions, outputting only signals.

Fast and a scalable algorithms and methods adaptable to available resources for computing (1) the DPRT on multicore CPUs by distributing the computation of the DPRT primary directions among the different cores, and (2) the DPRT on GPUs using parallel, distributed, and synchronized ray computations among the GPU cores with ray referring to one of the sums required for computing the DPRT or its inverse along a prime direction.