A system to detect and control noise in seismic surveys is provided. The system receives, responsive to a seismic wave generated by a source, seismic data detected by a sensor component of a seismic data acquisition unit. The system generates, for windows of the seismic data, Hough tensors for seismic data transforms in multiple dimensions. The system detects, based on a comparison of an eigenvector and eigenvalue of a canonical matrix of the Hough tensors with a historical eigenvector and eigenvalue of a historical canonical matrix of historical Hough tensors of historical seismic data, a first presence of noise in the seismic data. The first presence of noise can correspond to a noisy spectra pattern in a seismic data transform of the seismic data. The system provides, responsive to detection of the first presence of noise in the seismic data, a notification to adjust a characteristic of the seismic survey.

Techniques are disclosed relating to geophysical surveying. In various embodiments, a marine survey vessel may tow a plurality of streamers that each include a plurality of seismic sensors. Further, the survey vessel may tow a plurality of vibratory sources. In various embodiments, a first sweep may be performed, using one or more of the plurality of vibratory sources, for a first time interval. Further, in various embodiments, disclosed techniques may include recording, during the first time interval using the plurality of seismic sensors, seismic data on a tangible, computer-readable medium, thereby creating a geophysical data product.

Processes and systems described herein are directed to performing marine surveys with marine vibrators that emit orthogonal coded pseudo-random sweeps. In one aspect, coded pseudo-random signals are generated based on coded pseudo-random sequences. The coded pseudo-random sequences are used to activate the marine vibrators in a body of water above a subterranean formation. The activated marine vibrators generate orthogonal coded pseudo-random sweeps. A wavefield emitted from the subterranean formation in response to the orthogonal coded pseudo-random sweeps is detected at receivers located in a body of water. Seismic signals generated by the receivers may be cross-correlated with a signature of one of the orthogonal coded pseudo-random sweeps to obtain seismic data with incoherent residual noise.

A method of seismic acquisition using a dispersed-source array (DSA) comprising two or more sources. The method comprises determining, for each of the two or more sources of the DSA, an individual spectrally-banded waveform. For each of the two or more sources, a primary waveform is formed by repeating the individual spectrally-banded waveform. For each of the two or more sources, a secondary waveform is formed based on the primary waveform. The secondary waveform is spectrally shifted relative to the primary waveform such that secondary waveforms of any two of the two or more sources are spectrally non-overlapping. The blending operator based on the secondary waveform of each of the two or more sources is provided to the DSA. The method also includes performing deblended-data reconstruction of acquired seismic data using one or more properties of the blending operators of the two or more sources.

A computer implemented method for identifying reservoir facies in a subsurface region includes obtaining a set of seismic data points of both petrophysical and geophysical parameters relating to the subsurface region, identifying one or more correlated clusters of petrophysical parameters, generating, from the one or more correlated clusters of petrophysical parameters, one or more corresponding multi-dimensional clusters of seismic data points, storing, in a facies database, a multi-dimensional cluster center point for at least one multi-dimensional clusters, and recursively splitting the multi-dimensional clusters into distinct sub-clusters of seismic data points corresponding to facies types.

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.

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.

Techniques for determining a drill bit location includes identifying a plurality of acoustic energy signals received at a plurality of sets of acoustic receivers from a passive acoustic energy source that is part of a wellbore drilling system; processing the plurality of acoustic energy signals; determining a location of a drill bit of the wellbore drilling system based on the processed plurality of acoustic signals; and updating a geo-steering path of the drill bit based on the determined location of the drill bit.