Processes and systems described herein are directed to imaging a subterranean formation from seismic data recorded in a marine survey with moving marine vibrators. The marine vibrators generate random sweeps with random sweep signatures. Processes and systems generate an up-going pressure wavefield from measured pressure and vertical velocity wavefield data recorded in the marine survey and obtain a downgoing vertical acceleration wavefield that records source wavefields, directivity, source ghosts, and random signatures of the random sweeps. The downgoing vertical acceleration wavefield data is deconvolved from the up-going pressure wavefield to obtain a subsurface reflectivity wavefield that is used to generate an image of the subterranean formation with reduced contamination from source wavefields, directivity, source ghosts, and random signatures of the random sweeps.

Disclosed are methods of marine 3D seismic data acquisition that do not require compensation for winds and currents.

Determination of an impulse response at a subsurface image level can include extrapolation of an up-going pressure wavefield to a subsurface image level, extrapolation of a down-going velocity wavefield to the subsurface image level, and determination of the impulse response at the subsurface image level from a hypothetical seismic source by spectral division of the extrapolated up-going pressure wavefield by the extrapolated down-going velocity wavefield.

The present disclosure describes methods and systems for estimating dispersion spectra for full waveform sonic (FWS) logging. One computer-implemented method includes receiving FWS data, performing frequency-spatial (FX) transform on the FWS data, using a nonstationary predictive error filtering (PEF) inversion on the transformed FWS data to estimate local matrix L and matrix P, calculating an inverse covariance matrix based on the estimated local matrix L and matrix P, and obtaining a nonstationary maximum likelihood method (MLM) spectra based on the inverse covariance matrix.

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.

Processes and systems described herein are directed to imaging a subterranean formation from seismic data recorded in a marine survey with moving marine vibrators. The marine vibrators generate random sweeps with random sweep signatures. Processes and systems generate an up-going pressure wavefield from measured pressure and vertical velocity wavefield data recorded in the marine survey and obtain a downgoing vertical acceleration wavefield that records source wavefields, directivity, source ghosts, and random signatures of the random sweeps. The downgoing vertical acceleration wavefield data is deconvolved from the up-going pressure wavefield to obtain a subsurface reflectivity wavefield that is used to generate an image of the subterranean formation with reduced contamination from source wavefields. directivity, source ghosts, and random signatures of the random sweeps.

The present disclosure describes methods and systems for estimating dispersion spectra for full waveform sonic (FWS) logging. One computer-implemented method includes receiving FWS data, performing frequency-spatial (FX) transform on the FWS data, using a nonstationary predictive error filtering (PEF) inversion on the transformed FWS data to estimate local matrix L and matrix P, calculating an inverse covariance matrix based on the estimated local matrix L and matrix P, and obtaining a nonstationary maximum likelihood method (MLM) spectra based on the inverse covariance matrix.

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.

Disclosed are methods of marine 3D seismic data acquisition that do not require compensation for winds and currents.

Methods and systems of generating seismic images from primaries and multiples are described. Methods separate pressure data into up-going pressure data and down-going pressure data from pressure data and vertical velocity data. Irregularly spaced receiver coordinates of the down-going and up-going pressure data are regularized to grid points of a migration grid and interpolation is used to fill in down-going and up-going pressure data at grid points of the migration grid. A seismic image is calculated at grid points of the migration grid based on the interpolated and regularized down-going pressure data and the interpolated and regularized up-going pressure data. The seismic images are high-resolution, have lower signal-to-noise ratio than seismic images generated by other methods, and have reduced acquisition artifacts and crosstalk effects.