Abstract: Methods and apparatuses for directional designature in shot domain are provided. Azimuth and take-off angles are calculated for each record in the seismic data. Directional designature is then applied to the seismic data using a source signature dependent on the azimuth and take-off angles.

Abstract: Computing device, computer instructions and method for improving accuracy of an image of a surveyed subsurface. The method includes a step of receiving seismic data (D) recorded with seismic sensors; a step of deghosting the seismic data (D) using a sparseT-pinversion algorithm together with a low-rank optimization scheme and a high-cut filter to calculate ghost-free wave-fields (U0(t,x,y)) in the time-space domain; and a step of processing the ghost-free wave-fields (U0(t,x,y)) to generate the image of the surveyed subsurface.

Abstract: Computing device, computer instructions and method for jointly deghosting first and second 3-dimensional (3D) seismic vintages of a same subsurface. The method includes receiving the first 3D vintage; receiving the second 3D vintage, wherein the second 3D vintage is taken later in time than the first 3D vintage, over the same subsurface; jointly deghosting the first and second 3D vintages based on a common ghost-free model U0, a first vintage ghost-free model Ub, and a second vintage ghost-free model Um; and generating an image of the subsurface indicative of changes between the first and second 3D vintages based on the common ghost-free model U0.

Abstract: An input seismic data set is obtained, and a first inversion is conducted on a first frequency filtered input seismic data set across all slowness values in the plurality of slowness values using an inversion matrix with dipole modulation containing a plurality of Tau-p operators and a plurality of dipole terms to transform the first frequency filtered input seismic data set from a time domain to a Tau-p domain having a plurality of first inversion Tau-p coefficients. Each Tau-p coefficient associated with one of the plurality of slowness values. The first inversion Tau-p coefficients are used to identify a subset of slowness values from the plurality of slowness values, and iterative sparse Tau-p inversion using the identified subset of slowness values is performed. The input seismic data set is modified to generate a noise attenuated seismic data set based on the iterative sparse Tau-p inversion.

Abstract: Methods and systems for optimized receiver-based ghost filter generation are described. The optimized ghost filter self-determines its parameters based on an iterative calculation of recorded data transformed from a time-space domain to a Tau-P domain. An initial ghost filter prediction is made based on generating mirror data from the recorded data and using a least squares technique during a premigration stage.