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.

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.

The present disclosure includes a method including determining a spatial region for analysis and selecting a segment of time for analysis, analyzing and correcting a plurality of traces from a plurality of receivers using an iterative non-linear inversion algorithm, wherein each iteration of the non-linear algorithm corrects the plurality of traces using at least one set of parameters defining a microseismic event, determining whether a final stack value of the plurality of traces corrected based on the at least one set of parameters of a final iteration of the iterative non-linear inversion algorithm exceeds a predetermined threshold and upon a determination that the final stack value exceeds the predetermined threshold, detecting a microseismic event defined by the at least one set of parameters of final iteration. The present disclosure also includes associated systems and computer-readable media.

Computing device, computer instructions and method for processing energy at a free-surface reflection relating to an air-water interface. The method includes receiving input seismic data recorded with seismic sensors; receiving wave-height data that describes an actual shape of a top surface of a body of water; processing up-going energy at a receiver and down-going energy following a reflection at the sea-surface, using the input seismic data and a linear operator modified to take into account the wave-height data; and generating an image of the subsurface based on the up-going energy or the down-going energy or a combination of the input seismic data and one of the up-going or down-going energy.

Computing device, computer instructions and method for processing energy at a free-surface reflection relating to an air-water interface. The method includes receiving input seismic data recorded with seismic sensors; receiving wave-height data that describes an actual shape of a top surface of a body of water; processing up-going energy at a receiver and down-going energy following a reflection at the sea-surface, using the input seismic data and a linear operator modified to take into account the wave-height data; and generating an image of the subsurface based on the up-going energy or the down-going energy or a combination of the input seismic data and one of the up-going or down-going energy.

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.

Synthetic survey data is generated using a two-way or one-way wave propagator based on a current model of a target structure. The current model is modified to reduce a difference between the synthetic survey data and observed survey data, while maintaining unchanged a velocity component of the current model, where the modifying of the current model produces a modified model. The modified model is used to reduce an adverse effect of multiples in the target structure, or to promote a favorable effect of multiples in the target structure.

Synthetic survey data is generated using a two-way or one-way wave propagator based on a current model of a target structure. The current model is modified to reduce a difference between the synthetic survey data and observed survey data, while maintaining unchanged a velocity component of the current model, where the modifying of the current model produces a modified model. The modified model is used to reduce an adverse effect of multiples in the target structure, or to promote a favorable effect of multiples in the target structure.

A seismic modeling system is provided which may include a seismic model data storage device, and a processor cooperating with the seismic model data storage device to identify an obstruction within a seismic model spatial domain data set having a lateral boundary, and define a suspect area including spatial domain data within the lateral boundary and directly below the obstruction. The processor may further inpaint the suspect area in the seismic model spatial domain data set based upon an exemplar inpainting algorithm.

Embodiments of dynamic gain adjustments in seismic surveys are described. One method of acquiring a seismic survey includes determining an arrival time at a seismic receiver of a downgoing seismic wavefield associated with a seismic source based at least in part on an estimated position of the seismic source, an estimated position of the seismic receiver, or combinations thereof. The method also includes adjusting a gain of the seismic receiver based at least in part on the determined arrival time of the downgoing seismic wavefield in order to, for example, help prevent overdriving or clipping of the seismic receiver when the downgoing seismic wavefield arrives at or passes by the seismic receiver.