Received shot gathers are sorted to a common receiver gather. A target three-dimensional (3D) amplitude spectrum of seismic wavefield direct arrivals is computed from synthetic data. A 3D amplitude spectrum of seismic wavefield direct arrivals in field data is computed for each receiver. A matched filter is calculated from the 3D amplitude spectrum of field data to target response and applied to downgoing seismic wavefields separated from the common receiver gather to generate filtered downgoing seismic wavefields. Time-dependent smoothing of the filtered downgoing seismic wavefields is performed to generate smoothed downgoing seismic wavefields. A cross-correlation is calculated between upgoing seismic wavefields separated from the common receiver gather and the smoothed downgoing seismic wavefields.

Processes and systems described herein are directed to performing marine surveys with a moving vibrational source that emits a continuous source wavefield into a body of water above a subterranean formation. The continuous source wavefield is formed from multiple sweeps in which each sweep is emitted from the moving vibrational source into the body of water with a randomized phase and/or with a randomized sweep duration. Reflections from the subterranean formation are continuously recorded in seismic data as the moving vibrational source travels above the subterranean formation. Processes and systems include iteratively deconvolving the source wavefield from the continuously recorded seismic data to obtain an earth response in the common receiver domain with little to no harmful effects from spatial aliasing and residual crosstalk noise. The earth response may be processed to generate an image of the subterranean formation.

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.

The present disclosure describes methods and systems, including computer-implemented methods, computer program products, and computer systems, for generating a reflectivity model for a subsurface area. One method includes: receiving a set of seismic data associated with the subsurface area; generating analytic source wavefields; generating analytic residual wavefields based on the set of seismic data and an initial reflectivity model; decomposing the analytic source wavefields and the analytic residual wavefields to obtain down-going and up-going components of the analytic source wavefields and the analytic residual wavefields; calculating a gradient vector using the down-going components of the analytic source wavefields and the up-going components of the analytic residual wavefields; calculating a source illumination factor using the down-going components of the analytic source wavefields; calculating a preconditioned gradient vector, based on the gradient vector and the source illumination factor; and generating an updated reflectivity model based on the preconditioned gradient vector.

A survey method includes generating a first amplitude modulation signal by amplitude-modulating a carrier wave repeating the same pattern at a predetermined cycle in each of a plurality of vibrators with a modulation signal whose cycle is 1/m times the predetermined period and is different for each of the vibrators, transmitting the seismic wave based on the first amplitude modulation signal, generating a second amplitude modulation signal in one or more receivers, the second amplitude modulation signal being identical to the first amplitude modulation signal generated by any one of the seismic vibrators, generating a reception signal in each of the one or more receivers by receiving a synthetic seismic wave in which the seismic waves generated by the seismic vibrators are synthesized, calculating a correlation value between the reception signal and the second amplitude modulation signal, and analyzing characteristics of the medium on the basis of the correlation value.

New methods for calculating acquisition illumination are computationally less expensive in comparison with conventional methods. In one such new method, source wavefield propagations are calculated and assigned to corresponding zero-offset receivers. Further, the number of non-zero-offset receivers within the coverage of the shot at the source location is decimated. Such a method is most advantages in reverse time migration, in which all source wavefield propagations are already calculated. The receiver-side illumination for each shot can be obtained by summing up all the source-side illumination with the source located within receiver coverage. All the source-side illumination and receiver-side illumination can be summed up to get the acquisition illumination for the survey. The acquisition illumination can be used to value the acquisition system and to compensate the migration images.

A method for echo detection may comprise recording one or more reflected waveforms, segmenting the one or more reflected waveforms based at least in part on a firing pulse length, applying a shaped filter to each segment of the one or more reflected waveforms, decoupling the one or more reflected waveforms into a time-frequency energy map, extracting a firing frequency band time domain plot from the decoupled time-frequency map, identifying a maximum amplitude in the extracted firing frequency band of the one or more reflected waveforms as an excitation, and identifying a second maximum amplitude in the extracted firing frequency band of the one or more reflected waveforms as an echo. A system for echo detection may comprise a digital signal processor, a transmitter, a transducer, a receiver, an analog to digital converter configured to digitize the measurement, and an information handling system.

A distributed acoustic sensor is positioned within a wellbore of a geologic formation. Seismic waves are detected using the distributed acoustic sensor. A raw seismic profile is generated based on the detected seismic waves. Resonant noise is identified and reduced in seismic data associated with the raw seismic profile.

Methods of geophysical modeling and inversion are disclosed. A sparse domain is defined for a geophysical model, over which a sparse model result is computed. A full model result is then resolved by interpolation over the sparse domain. The full model result may be used as the forward modeling result in a geophysical inversion process. Reconstruction error, or model error, or both may be used to adjust the sparse domain, the model, or the geophysical basis of the model.

A method for applying separated wave-field imaging onshore (1) by artificially creating up-going and down-going fields and (2) by using these fields in a migration algorithm. If there are any surface multiples in the data, the resulting image created using the migration algorithm will be distorted by the unknown free-surface reflection coefficient. In fact, the surface multiples may be generated with a complex series of reflection coefficients. The distortions found in the resulting image created using the migration algorithm are then removed.