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 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.

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.

Techniques are disclosed relating to geophysical surveying. In various embodiments, a computer system may access seismic data for a geological formation, where the seismic data is recorded, using one or more sensors, during a seismic survey in which a first vibratory source was driven using a first digital code for at least a first time interval. The first digital code, in some embodiments, may include a first plurality of subsections corresponding to portions of the first time interval. In some embodiments, the computer system may image a first location of the geological formation using a correlation of only a first sub-section of the first plurality of sub-sections with the seismic data. Further, in some embodiments, the computer system may image a second location of the geological formation using a correlation of two or more of the first plurality of sub-sections with the seismic data.

Generating spectrally enhanced seismic data expresses seismic data as a convolution of reflectivity and a seismic source wavelet. This seismic source wavelet varies over a sampling interval and defining a total amount of energy over the sampling interval. An enhanced seismic source wavelet that is a single-valued energy spike that yields the total amount of energy over the sampling interval is generated. In addition, the reflectivity is modified to preserve amplitude variation with angle. The reflectivity is convoluted with the enhanced seismic source wavelet and residual energy is added to the convolution to generate the spectrally enhanced seismic data.

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.

Methods and apparatus for interferometric seismic imaging and creation of a high-resolution three-dimensional seismic volume in proximity to a wellbore are described. In contrast to current methods that deliver positions of microseismic events using a small fraction of the microseismic wavefield, the present invention provides for the analysis of a full microseismic wavefield. In a preferred embodiment, the method includes creating a planar image slice and/or corridor between the location of one or more microseismic events and one or several sensor arrays to produce a multitude of azimuthally orientated planar image slices and/or corridors. The method further includes adding the planar image slice and/or corridor contributions to create a single three-dimensional volume for analysis and rendering.

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 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.

Generating spectrally enhanced seismic data expresses seismic data as a convolution of reflectivity and a seismic source wavelet. This seismic source wavelet varies over a sampling interval and defining a total amount of energy over the sampling interval. An enhanced seismic source wavelet that is a single-valued energy spike that yields the total amount of energy over the sampling interval is generated. In addition, the reflectivity is modified to preserve amplitude variation with angle. The reflectivity is convoluted with the enhanced seismic source wavelet and residual energy is added to the convolution to generate the spectrally enhanced seismic data.