Seismic data are obtained by recording simultaneously in seismic streamer, acquired by activating approximately simultaneously two or more seismic sources towed at two positions in the vicinity of seismic streamers. A residual is updated iteratively for an inversion solution for the activations of the two or more seismic sources. The iterative updating of the residuals utilizes a sequence of overlapping temporal windows containing reflection events and utilizes normal moveout corrections based on largest reflection events in each temporal window. A final updated residual is added to a final updated model result.

The disclosure provides systems and methods to enhance pre-stack data for seismic data analysis by: sorting the reflection seismic data acquired from cross-spread gathers into sets of data sections; performing data enhancement on the sets of data sections to generate enhanced traces by: (i) applying forward normal-moveout (NMO) corrections such that arrival times of primary reflection events become more flat, (ii) estimating beamforming parameters including a nonlinear traveltime surface and a summation aperture, (iii) generating enhanced traces that combine contributions from original traces in the sets of data sections, and (iv) applying inverse NMO corrections to the enhanced traces such that temporal rearrangements due to the forward NMO corrections are undone.

This specification describes workflows for, but is not limited to, performing full waveform inversion (FWI) to build high resolution velocity models to improve the accuracy of seismic imaging of a subterranean formation. This specification describes processes to automatically edit and enhance S/N quality of seismic data (such as land seismic data) to prepare the datasets for FWI. The methods for automatic corrections and pre-processing include: automatic iterative surface-consistent residual statics calculation, automatic rejection of anomalous traces (such as dead traces), and the automatic correction of surface-consistent amplitude anomalies (such as by scalar or deconvolution approaches). The operations include automatic “muting” of noise before first arrivals.

A sonic logging method is provided that transmits acoustic signals using a high order acoustic source and processes waveform data to identify a set of arrival events and time picks by automatic and/or manual methods. Ray tracing inversion is carried out for each arrival event over a number of possible raypath types that include at least one polarized shear raypath type to determine two-dimensional reflector positions and predicted inclination angles of the arrival event for the possible raypath types. One or more three-dimensional slowness-time coherence representations are generated for the arrival event and raypath type(s) and evaluated to determine azimuth, orientation and raypath type of a corresponding reflector. The method outputs a three-dimensional position and orientation for at least one reflector. The information derived from the method can be conveyed in various displays and plots and structured formats for reservoir understanding and also output for use in reservoir analysis and other applications.

Methods and systems including computer programs encoded on a computer storage medium, for utilizing equivalent linear velocity for first arrival picking of seismic refraction. In one aspect, a method includes receiving data for the shot gather record, generating a diving wave equation curve for a particular parameter pair of multiple parameter pairs, and integrating the shot gather record data corresponding to the diving wave equation curve over a selected range of offsets of the shot gather to generate an equivalent linear velocity value for the particular parameter pair and the shot gather record data, selecting, from the equivalent linear velocity values for the plurality of parameter pairs, a greatest equivalent linear velocity value of the equivalent linear velocity values, the greatest equivalent linear velocity value corresponding to a first-arrival parameter pair, and determining, using the first-arrival parameter pair, a set of first-arrival onsets for the selected sub-range of offsets.

A method for performing a formation-related operation based on corrected vertical seismic profile (VSP) data of an earth formation includes performing a VSP survey and applying a normal moveout (NMO) correction equation to the survey data that is a function of source offset to wellhead. The method also includes solving the NMO correction equation using a simulated annealing algorithm having an object function that is a coherence coefficient of semblance analysis of an NMO corrected reflection event within a time window to provide NMO corrected data. The method further includes performing the formation-related operation at at least one of a location, a depth and a depth interval based on the VSP NMO corrected data.

A method may include obtaining seismic data for a geological region of interest. The method may further include determining various velocity semblance values for the geological region of interest using a time window and the seismic data. The method may further include determining, automatically by a computer processor, one or more stacking velocities for the geological region of interest using a traced path based on the velocity semblance values and a path tracing algorithm. The path tracing algorithm may recursively determine an accumulated amplitude array based on the velocity semblance values. The path tracing algorithm may further determine the traced path among the velocity semblance values and the accumulated amplitude array, the traced path corresponding to the one or more stacking velocities. The method may further include generating a velocity model of the geological region of interest using the one or more stacking velocities.

The present invention provides a seismic signal processing method, device and system. The method comprises: obtaining an offset of a reflected seismic signal at a sampling point and the corresponding reflected wave arrival time; constructing a non-hyperbolic dynamic correction formula based on Pade approximation according to the offset of the reflected seismic signal at the sampling point and the corresponding reflected wave arrival time; extracting a vertical propagation velocity and anisotropy parameters of the reflected seismic signal according to the non-hyperbolic dynamic correction formula constructed based on Pade approximation.

Geophysical prospecting may be achieved using borehole seismic data and processing velocity seismic profiles using downward continuation to simulate the seismic source being at the depth of the borehole receivers. Such methods may involve collecting seismic data for a subterranean formation with at least one borehole receiver; grouping the seismic data into a one common receiver gather corresponding to each borehole receiver; performing a downward continuation on at least one of the common receiver gathers to produce corresponding downward continued common receiver gathers; performing a normal moveout analysis on at least one of the downward continued common receiver gathers to produce corresponding semblance velocity spectra; and analyzing at least one of the semblance velocity spectra for a zone of interest in the subterranean formation.

A method includes receiving, via a processor, input data based upon received seismic data, migrating, via the processor, the input data via a pre-stack depth migration technique to generate migrated input data, encoding, via the processor, the input data via an encoding function as a migration attribute to generate encoded input data having a migration function that is non-monotonic versus an attribute related to the input data, migrating, via the processor, the encoded input data via the pre-stack depth migration technique to generate migrated encoded input data, and generating an estimated common image gather based upon the migrated input data and the migrated encoded input data. The method also includes generating a seismic image utilizing the estimated common image gather, wherein the seismic image represents hydrocarbons in a subsurface region of the Earth or subsurface drilling hazards.