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.

A method may include providing a sensor in a first wellbore segment, providing a sensor in a second wellbore segment, observing upgoing acoustic waves or downgoing acoustic waves with the sensors, and separating the upgoing acoustic waves and/or the downgoing acoustic waves from a total wavefield. The first wellbore segment and the second wellbore segment may be separated by a distance. At least one of the wellbore segments may be non-vertical and/or the first wellbore segment may not be parallel to the second wellbore segment. The first wellbore segment may be part of a first set of wellbores and the second wellbore segment may be part of a second set of wellbores. The separated upgoing and downgoing acoustic waves may be used to generate deghosted data.

Processes and systems for deblending blended seismic data with attenuated source signatures and source ghost are described. Processes and systems compute blended upgoing pressure wavefield based on blended pressure wavefield and blended vertical velocity wavefield recorded in a marine survey of a subterranean formation. Downgoing vertical velocity wavefield is computed based on near-field pressure wavefields generated by source elements of sources activated in the marine survey. Deblended wavefield is computed based on the blended upgoing pressure wavefield and the downgoing vertical velocity source wavefield. The deblended wavefield may be used to generate an image of the subterranean formation with the source signatures and source ghosts contained in the blended pressure wavefield and blended vertical velocity wavefield.

Processes and systems for deblending blended seismic data with attenuated source signatures and free-surface effects are described. The blended seismic data may have been recorded in a marine survey in which multiple sources are activated in the body of water above a subterranean formation. Receivers record overlapping pressure and vertical velocity wavefield responses from the subterranean formation as corresponding blended pressure wavefield and blended vertical velocity wavefield. Processes and systems compute an upgoing pressure wavefield and a downgoing vertical velocity wavefield based on the blended pressure wavefield and blended vertical velocity wavefield. Deblended primary pressure wavefields are computed based on the upgoing pressure and downgoing vertical velocity. The deblended primary pressure wavefields may be used to generate images of the subterranean formation that are substantially free of source signatures and free-surface effects.

A method for geological formation analysis may include collecting time-lapsed well-based measurement data from a first borehole in a geological formation over a measurement time period, and collecting time-lapsed electromagnetic (EM) cross-well measurement data via a plurality of spaced-apart second boreholes in the geological formation over the measurement time period. The method may further include determining simulated changes to a hydrocarbon resource in the geological formation over the measurement time period based upon a geological model using a processor, and using the processor to determine if the simulated changes are within an error threshold of the time-lapsed well-based measurement data and the time-lapsed cross-well EM measurement data. If the simulated changes are not within the error threshold, then the geological model may be updated.

Processes and systems for imaging a subterranean formation using continuously recorded seismic data obtained during a marine seismic geophysical survey of the subterranean formation are described herein. The processes and systems compute upgoing pressure data at stationary-receiver locations, and low-frequency noise attenuation processes and systems are applied to the upgoing pressure wavefield data to obtain low-frequency noise attenuated upgoing pressure wavefield data. An image of the subterranean formation, or data indicative thereof, may be generated using the low-frequency noise attenuated upgoing pressure wavefield data at stationary-receiver locations.

Embodiments included herein are directed towards a marine seismic streamer. The seismic streamer may include an outer skin formed in longitudinally extending tubular shape, an inner surface of the outer skin defining an internal volume containing a gel substance. The seismic streamer may also include a plurality of micro-electro-mechanical (MEMS) sensors and plurality of hydrophones associated with the outer skin, wherein in the plurality of MEMS sensors are spaced non-uniformly in the seismic streamer along an axial direction of the streamer, such that not more than 100 MEMS sensors are located in the seismic streamer over a continuous 100 meter axial length of seismic streamer. The seismic streamer may further include an electronics system extending axially through an inside portion of the outer skin and a strength member core extending axially through an inside portion of the outer skin.

Prediction and subtraction of multiple diffractions may include transforming previously acquired seismic data from a time-space domain to a transformed domain using a dictionary of compressive basis functions and separating, within the transformed previously acquired seismic data, a first portion and a second portion of the transformed previously acquired seismic data. Prediction and subtraction of multiple diffractions may also include predicting a plurality of multiple diffractions based on the separated first and second portions and adaptively subtracting the predicted multiple diffractions from the previously acquired seismic data. Prediction and subtraction of multiple diffractions may also include inverse transforming a particular seismic data set from the transformed domain to the time-space domain.

Computing device, computer instructions and method for correcting an image, of a surveyed surface, due to a free-surface reflection. The method includes calculating a free surface reflection operator for a seismic source displaced in water based on a position of the source, and an air-water interface datum; receiving recorded seismic data d recorded with seismic sensors (r), wherein the recorded seismic data is associated with a pressure and/or a particle motion produced by a seismic wave in earth; correcting the recorded seismic data d based on the free surface reflection operator to obtain transformed seismic data; and generating an image of the surveyed subsurface, based on the transformed seismic data, wherein the image is indicative of various layers of the earth. The free surface reflection operator varies while a source signal is being emitted by the source.

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 sparse.sub.T-Pinversion algorithm together with a low-rank optimization scheme and a high-cut filter to calculate ghost-free wave-fields (U.sub.0(t,x,y)) in the time-space domain; and a step of processing the ghost-free wave-fields (U.sub.0(t,x,y)) to generate the image of the surveyed subsurface.